mapnik/plugins/input/shape/shape_index_featureset.cpp

185 lines
6.4 KiB
C++
Raw Normal View History

2006-03-31 12:32:02 +02:00
/*****************************************************************************
*
2006-03-31 12:32:02 +02:00
* This file is part of Mapnik (c++ mapping toolkit)
2005-06-14 17:06:59 +02:00
*
* Copyright (C) 2011 Artem Pavlenko
2005-06-14 17:06:59 +02:00
*
2006-03-31 12:32:02 +02:00
* 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,
2005-06-14 17:06:59 +02:00
* but WITHOUT ANY WARRANTY; without even the implied warranty of
2006-03-31 12:32:02 +02:00
* 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
2005-06-14 17:06:59 +02:00
*
2006-03-31 12:32:02 +02:00
*****************************************************************************/
2005-06-14 17:06:59 +02:00
// stl
#include <fstream>
2005-06-14 17:06:59 +02:00
// mapnik
2012-04-08 02:20:56 +02:00
#include <mapnik/debug.hpp>
#include <mapnik/feature_factory.hpp>
// boost
2014-10-22 02:44:47 +02:00
#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wunused-parameter"
#pragma GCC diagnostic ignored "-Wunused-local-typedef"
#include <boost/algorithm/string.hpp>
#ifdef SHAPE_MEMORY_MAPPED_FILE
#include <boost/interprocess/streams/bufferstream.hpp>
#endif
2014-10-22 02:44:47 +02:00
#pragma GCC diagnostic pop
#include "shape_index_featureset.hpp"
#include "shape_utils.hpp"
2013-01-14 10:49:35 +01:00
#include "shp_index.hpp"
using mapnik::feature_factory;
using mapnik::geometry_type;
2005-06-14 17:06:59 +02:00
template <typename filterT>
shape_index_featureset<filterT>::shape_index_featureset(filterT const& filter,
std::unique_ptr<shape_io> && shape_ptr,
std::set<std::string> const& attribute_names,
std::string const& encoding,
std::string const& shape_name,
int row_limit)
: filter_(filter),
ctx_(std::make_shared<mapnik::context_type>()),
shape_ptr_(std::move(shape_ptr)),
tr_(new mapnik::transcoder(encoding)),
row_limit_(row_limit),
count_(0),
feature_bbox_()
2005-06-14 17:06:59 +02:00
{
shape_ptr_->shp().skip(100);
setup_attributes(ctx_, attribute_names, shape_name, *shape_ptr_,attr_ids_);
2012-04-17 07:22:35 +02:00
auto index = shape_ptr_->index();
2010-09-06 21:20:59 +02:00
if (index)
2005-06-14 17:06:59 +02:00
{
#ifdef SHAPE_MEMORY_MAPPED_FILE
//shp_index<filterT,stream<mapped_file_source> >::query(filter, index->file(), offsets_);
shp_index<filterT,boost::interprocess::ibufferstream>::query(filter, index->file(), offsets_);
#else
shp_index<filterT,std::ifstream>::query(filter, index->file(), offsets_);
#endif
2005-06-14 17:06:59 +02:00
}
std::sort(offsets_.begin(), offsets_.end());
2012-02-02 02:37:35 +01:00
MAPNIK_LOG_DEBUG(shape) << "shape_index_featureset: Query size=" << offsets_.size();
2012-02-02 02:37:35 +01:00
itr_ = offsets_.begin();
2005-06-14 17:06:59 +02:00
}
template <typename filterT>
feature_ptr shape_index_featureset<filterT>::next()
{
if (row_limit_ && count_ >= row_limit_)
{
return feature_ptr();
}
while ( itr_ != offsets_.end())
2005-06-14 17:06:59 +02:00
{
shape_ptr_->move_to(*itr_++);
shape_file::record_type record(shape_ptr_->reclength_ * 2);
shape_ptr_->shp().read_record(record);
int type = record.read_ndr_integer();
feature_ptr feature(feature_factory::create(ctx_,shape_ptr_->id_));
switch (type)
2010-07-25 21:55:40 +02:00
{
case shape_io::shape_point:
case shape_io::shape_pointm:
case shape_io::shape_pointz:
{
double x = record.read_double();
double y = record.read_double();
std::unique_ptr<geometry_type> point(new geometry_type(mapnik::geometry_type::types::Point));
point->move_to(x, y);
feature->paths().push_back(point.release());
break;
2005-06-14 17:06:59 +02:00
}
case shape_io::shape_multipoint:
case shape_io::shape_multipointm:
case shape_io::shape_multipointz:
{
shape_io::read_bbox(record, feature_bbox_);
if (!filter_.pass(feature_bbox_)) continue;
int num_points = record.read_ndr_integer();
for (int i = 0; i < num_points; ++i)
{
double x = record.read_double();
double y = record.read_double();
std::unique_ptr<geometry_type> point(new geometry_type(mapnik::geometry_type::types::Point));
point->move_to(x, y);
feature->paths().push_back(point.release());
}
break;
}
case shape_io::shape_polyline:
case shape_io::shape_polylinem:
case shape_io::shape_polylinez:
2005-06-14 17:06:59 +02:00
{
shape_io::read_bbox(record, feature_bbox_);
if (!filter_.pass(feature_bbox_)) continue;
shape_io::read_polyline(record,feature->paths());
break;
2005-06-14 17:06:59 +02:00
}
case shape_io::shape_polygon:
case shape_io::shape_polygonm:
case shape_io::shape_polygonz:
{
shape_io::read_bbox(record, feature_bbox_);
if (!filter_.pass(feature_bbox_)) continue;
shape_io::read_polygon(record,feature->paths());
break;
}
default :
MAPNIK_LOG_DEBUG(shape) << "shape_index_featureset: Unsupported type" << type;
return feature_ptr();
}
2013-07-23 22:00:16 +02:00
// FIXME: https://github.com/mapnik/mapnik/issues/1020
feature->set_id(shape_ptr_->id_);
if (attr_ids_.size())
{
shape_ptr_->dbf().move_to(shape_ptr_->id_);
std::vector<int>::const_iterator itr = attr_ids_.begin();
std::vector<int>::const_iterator end = attr_ids_.end();
try
{
for (; itr!=end; ++itr)
{
shape_ptr_->dbf().add_attribute(*itr, *tr_, *feature);
}
}
catch (...)
{
MAPNIK_LOG_ERROR(shape) << "Shape Plugin: error processing attributes";
}
}
++count_;
return feature;
}
MAPNIK_LOG_DEBUG(shape) << "shape_index_featureset: " << count_ << " features";
return feature_ptr();
2005-06-14 17:06:59 +02:00
}
template <typename filterT>
shape_index_featureset<filterT>::~shape_index_featureset() {}
template class shape_index_featureset<mapnik::filter_in_box>;
template class shape_index_featureset<mapnik::filter_at_point>;