mapnik/plugins/input/geos/geos_featureset.cpp

/*****************************************************************************
 * 
 * This file is part of Mapnik (c++ mapping toolkit)
 *
 * Copyright (C) 2010 Artem Pavlenko
 *
 * This library is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2.1 of the License, or (at your option) any later version.
 *
 * This library is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
 * License along with this library; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA
 *
 *****************************************************************************/
//$Id$

#include <iostream>
#include <fstream>

// mapnik
#include <mapnik/global.hpp>
#include <mapnik/datasource.hpp>
#include <mapnik/box2d.hpp>
#include <mapnik/geometry.hpp>
#include <mapnik/feature.hpp>
#include <mapnik/feature_layer_desc.hpp>
#include <mapnik/wkb.hpp>
#include <mapnik/unicode.hpp>
#include <mapnik/feature_factory.hpp>

// ogr
#include "geos_featureset.hpp"

using std::clog;
using std::endl;

using mapnik::query;
using mapnik::box2d;
using mapnik::coord2d;
using mapnik::CoordTransform;
using mapnik::Feature;
using mapnik::feature_ptr;
using mapnik::geometry_utils;
using mapnik::transcoder;
using mapnik::feature_factory;


geos_featureset::geos_featureset(GEOSGeometry* geometry,
                                 GEOSGeometry* extent,
                                 int identifier,
                                 const std::string& field,
                                 const std::string& field_name,
                                 const std::string& encoding,
                                 bool multiple_geometries)
   : geometry_(geometry),
     tr_(new transcoder(encoding)),
     extent_(extent),
     identifier_(identifier),
     field_(field),
     field_name_(field_name),
     multiple_geometries_(multiple_geometries),
     already_rendered_(false)
{
}

geos_featureset::~geos_featureset() 
{
}

feature_ptr geos_featureset::next()
{
    if (! already_rendered_)
    {
        already_rendered_ = true;
        
        if (GEOSisValid(geometry_) && ! GEOSisEmpty(geometry_))
        {
            bool render_geometry = true;
            
            if (*extent_ != NULL && GEOSisValid(*extent_) && !GEOSisEmpty(*extent_))
            {
                const int type = GEOSGeomTypeId(*extent_);
                render_geometry = false;

                switch ( type )
                {
                case GEOS_POINT:
                    if (GEOSIntersects(*extent_, geometry_))
                    {
                        render_geometry = true;
                    }
                    break;
                case GEOS_POLYGON:
                    if (GEOSContains(*extent_, geometry_)
                        || GEOSWithin(geometry_, *extent_)
                        || GEOSEquals(geometry_, *extent_))
                    {
                        render_geometry = true;
                    }
                    break;
               default:
#ifdef MAPNIK_DEBUG
                    clog << "GEOS Plugin: unknown extent geometry_type=" << type << endl;
#endif
                    break;
               }
            }

            if (render_geometry)
            {
                geos_wkb_ptr wkb(geometry_);
                if (wkb.is_valid())
                {
                    feature_ptr feature(feature_factory::create(identifier_));

                    geometry_utils::from_wkb(feature->paths(),
                                             wkb.data(),
                                             wkb.size(),
                                             multiple_geometries_);

                    if (field_ != "")
                    {
                        boost::put(*feature, field_name_, tr_->transcode(field_.c_str()));
                    }
                    
                    return feature;
                }
            }
        }
    }

    return feature_ptr();
}