2006-03-31 10:32:02 +00:00
|
|
|
/*****************************************************************************
|
2011-11-14 03:37:50 +00:00
|
|
|
*
|
2006-03-31 10:32:02 +00:00
|
|
|
* This file is part of Mapnik (c++ mapping toolkit)
|
2005-06-14 15:06:59 +00:00
|
|
|
*
|
2011-10-22 13:27:28 +00:00
|
|
|
* Copyright (C) 2011 Artem Pavlenko
|
2005-06-14 15:06:59 +00:00
|
|
|
*
|
2006-03-31 10:32:02 +00: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 15:06:59 +00:00
|
|
|
* but WITHOUT ANY WARRANTY; without even the implied warranty of
|
2006-03-31 10:32:02 +00: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 15:06:59 +00:00
|
|
|
*
|
2006-03-31 10:32:02 +00:00
|
|
|
*****************************************************************************/
|
2005-06-14 15:06:59 +00:00
|
|
|
|
2011-10-22 13:27:28 +00:00
|
|
|
// stl
|
|
|
|
#include <iostream>
|
|
|
|
|
2011-05-16 23:41:34 +00:00
|
|
|
// mapnik
|
2012-04-08 00:20:56 +00:00
|
|
|
#include <mapnik/debug.hpp>
|
2011-05-16 23:41:34 +00:00
|
|
|
#include <mapnik/feature_factory.hpp>
|
2013-01-04 17:23:06 +00:00
|
|
|
#include <mapnik/unicode.hpp>
|
|
|
|
|
|
|
|
// boost
|
|
|
|
#include <boost/make_shared.hpp>
|
2011-05-16 23:41:34 +00:00
|
|
|
|
2006-10-04 11:22:18 +00:00
|
|
|
#include "shape_featureset.hpp"
|
2012-01-12 12:18:50 +00:00
|
|
|
#include "shape_utils.hpp"
|
2005-06-14 15:06:59 +00:00
|
|
|
|
2011-05-16 23:41:34 +00:00
|
|
|
using mapnik::geometry_type;
|
|
|
|
using mapnik::feature_factory;
|
2012-01-12 10:04:08 +00:00
|
|
|
using mapnik::context_ptr;
|
2012-02-02 01:37:35 +00:00
|
|
|
|
2005-06-14 15:06:59 +00:00
|
|
|
template <typename filterT>
|
2012-01-12 12:18:50 +00:00
|
|
|
shape_featureset<filterT>::shape_featureset(filterT const& filter,
|
|
|
|
std::string const& shape_name,
|
|
|
|
std::set<std::string> const& attribute_names,
|
2007-02-14 21:55:04 +00:00
|
|
|
std::string const& encoding,
|
2011-08-26 07:02:36 +00:00
|
|
|
long file_length,
|
|
|
|
int row_limit)
|
2011-11-14 03:37:50 +00:00
|
|
|
: filter_(filter),
|
|
|
|
shape_(shape_name, false),
|
|
|
|
query_ext_(),
|
2012-12-11 09:37:53 +00:00
|
|
|
feature_bbox_(),
|
2011-11-14 03:37:50 +00:00
|
|
|
tr_(new transcoder(encoding)),
|
|
|
|
file_length_(file_length),
|
2012-01-12 12:28:53 +00:00
|
|
|
row_limit_(row_limit),
|
2012-02-02 16:58:20 +00:00
|
|
|
count_(0),
|
2013-09-20 13:00:11 +00:00
|
|
|
ctx_(std::make_shared<mapnik::context_type>())
|
2005-06-14 15:06:59 +00:00
|
|
|
{
|
2010-01-20 15:25:58 +00:00
|
|
|
shape_.shp().skip(100);
|
2012-01-12 12:18:50 +00:00
|
|
|
setup_attributes(ctx_, attribute_names, shape_name, shape_,attr_ids_);
|
2005-06-14 15:06:59 +00:00
|
|
|
}
|
|
|
|
|
|
|
|
template <typename filterT>
|
2005-09-08 13:09:04 +00:00
|
|
|
feature_ptr shape_featureset<filterT>::next()
|
2005-06-14 15:06:59 +00:00
|
|
|
{
|
2012-11-28 22:15:53 +00:00
|
|
|
if (row_limit_ && count_ >= row_limit_)
|
2011-10-22 13:27:28 +00:00
|
|
|
{
|
2011-08-26 07:02:36 +00:00
|
|
|
return feature_ptr();
|
2011-10-22 13:27:28 +00:00
|
|
|
}
|
|
|
|
|
2011-08-26 07:02:36 +00:00
|
|
|
|
2012-12-11 09:37:53 +00:00
|
|
|
while (shape_.shp().pos() < std::streampos(file_length_ * 2))
|
2011-08-12 16:43:28 +00:00
|
|
|
{
|
2012-12-11 09:37:53 +00:00
|
|
|
shape_.move_to(shape_.shp().pos());
|
|
|
|
shape_file::record_type record(shape_.reclength_ * 2);
|
|
|
|
shape_.shp().read_record(record);
|
|
|
|
int type = record.read_ndr_integer();
|
2011-11-14 03:37:50 +00:00
|
|
|
|
2012-12-11 09:37:53 +00:00
|
|
|
// skip null shapes
|
|
|
|
if (type == shape_io::shape_null) continue;
|
2011-11-14 03:37:50 +00:00
|
|
|
|
2012-12-11 09:37:53 +00:00
|
|
|
feature_ptr feature(feature_factory::create(ctx_, shape_.id_));
|
|
|
|
switch (type)
|
2010-01-29 01:54:15 +00:00
|
|
|
{
|
2012-12-11 09:37:53 +00:00
|
|
|
case shape_io::shape_point:
|
|
|
|
case shape_io::shape_pointm:
|
|
|
|
case shape_io::shape_pointz:
|
2010-01-29 01:54:15 +00:00
|
|
|
{
|
2012-12-11 09:37:53 +00:00
|
|
|
double x = record.read_double();
|
|
|
|
double y = record.read_double();
|
|
|
|
if (!filter_.pass(mapnik::box2d<double>(x,y,x,y)))
|
|
|
|
continue;
|
2013-09-03 11:15:31 +00:00
|
|
|
std::unique_ptr<geometry_type> point(new geometry_type(mapnik::geometry_type::types::Point));
|
2011-10-22 13:27:28 +00:00
|
|
|
point->move_to(x, y);
|
2013-04-19 08:50:50 +00:00
|
|
|
feature->paths().push_back(point.release());
|
2012-12-11 09:37:53 +00:00
|
|
|
break;
|
2010-01-29 01:54:15 +00:00
|
|
|
}
|
2012-12-11 09:37:53 +00:00
|
|
|
case shape_io::shape_multipoint:
|
|
|
|
case shape_io::shape_multipointm:
|
|
|
|
case shape_io::shape_multipointz:
|
2010-01-29 01:54:15 +00:00
|
|
|
{
|
2012-12-11 09:37:53 +00:00
|
|
|
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)
|
2010-01-29 01:59:50 +00:00
|
|
|
{
|
2012-12-11 09:37:53 +00:00
|
|
|
double x = record.read_double();
|
|
|
|
double y = record.read_double();
|
2013-09-03 11:15:31 +00:00
|
|
|
std::unique_ptr<geometry_type> point(new geometry_type(mapnik::geometry_type::types::Point));
|
2012-12-11 09:37:53 +00:00
|
|
|
point->move_to(x, y);
|
2013-04-19 08:50:50 +00:00
|
|
|
feature->paths().push_back(point.release());
|
2010-01-29 01:59:50 +00:00
|
|
|
}
|
2012-12-11 09:37:53 +00:00
|
|
|
break;
|
2010-01-29 01:54:15 +00:00
|
|
|
}
|
2011-10-22 13:27:28 +00:00
|
|
|
|
2012-12-11 09:37:53 +00:00
|
|
|
case shape_io::shape_polyline:
|
|
|
|
case shape_io::shape_polylinem:
|
|
|
|
case shape_io::shape_polylinez:
|
|
|
|
{
|
|
|
|
shape_io::read_bbox(record, feature_bbox_);
|
|
|
|
if (!filter_.pass(feature_bbox_)) continue;
|
|
|
|
shape_io::read_polyline(record, feature->paths());
|
|
|
|
break;
|
2010-01-29 01:54:15 +00:00
|
|
|
}
|
2012-12-11 09:37:53 +00: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_featureset: Unsupported type" << type;
|
|
|
|
return feature_ptr();
|
|
|
|
}
|
|
|
|
|
2012-01-17 18:33:35 +00:00
|
|
|
// FIXME: https://github.com/mapnik/mapnik/issues/1020
|
2011-05-20 01:41:28 +00:00
|
|
|
feature->set_id(shape_.id_);
|
2010-01-29 01:54:15 +00:00
|
|
|
if (attr_ids_.size())
|
|
|
|
{
|
|
|
|
shape_.dbf().move_to(shape_.id_);
|
2011-10-22 13:27:28 +00:00
|
|
|
std::vector<int>::const_iterator itr = attr_ids_.begin();
|
|
|
|
std::vector<int>::const_iterator end = attr_ids_.end();
|
2011-11-14 03:37:50 +00:00
|
|
|
try
|
2010-01-29 01:54:15 +00:00
|
|
|
{
|
2011-10-22 13:27:28 +00:00
|
|
|
for (; itr != end; ++itr)
|
2011-11-14 03:37:50 +00:00
|
|
|
{
|
2011-10-22 13:27:28 +00:00
|
|
|
shape_.dbf().add_attribute(*itr, *tr_, *feature); //TODO optimize!!!
|
2010-07-25 19:55:47 +00:00
|
|
|
}
|
2011-08-11 11:04:19 +00:00
|
|
|
}
|
|
|
|
catch (...)
|
|
|
|
{
|
2012-04-09 01:00:51 +00:00
|
|
|
MAPNIK_LOG_ERROR(shape) << "Shape Plugin: error processing attributes";
|
2010-01-29 01:54:15 +00:00
|
|
|
}
|
|
|
|
}
|
2012-12-11 09:37:53 +00:00
|
|
|
++count_;
|
2010-01-29 01:54:15 +00:00
|
|
|
return feature;
|
2010-01-20 15:25:58 +00:00
|
|
|
}
|
2012-04-09 10:05:49 +00:00
|
|
|
|
2012-12-11 09:37:53 +00:00
|
|
|
MAPNIK_LOG_DEBUG(shape) << "shape_featureset: Total shapes read=" << count_;
|
|
|
|
return feature_ptr();
|
2005-06-14 15:06:59 +00:00
|
|
|
}
|
|
|
|
|
|
|
|
template <typename filterT>
|
|
|
|
shape_featureset<filterT>::~shape_featureset() {}
|
|
|
|
|
2007-03-16 10:11:37 +00:00
|
|
|
template class shape_featureset<mapnik::filter_in_box>;
|
|
|
|
template class shape_featureset<mapnik::filter_at_point>;
|