/***************************************************************************** * * This file is part of Mapnik (c++ mapping toolkit) * * Copyright (C) 2012 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 * *****************************************************************************/ #include "topojson_datasource.hpp" #include "topojson_featureset.hpp" #include #include // boost #include #include // mapnik #include #include #include #include #include #include #include using mapnik::datasource; using mapnik::parameters; DATASOURCE_PLUGIN(topojson_datasource) struct attr_value_converter : public mapnik::util::static_visitor { mapnik::eAttributeType operator() (mapnik::value_integer /*val*/) const { return mapnik::Integer; } mapnik::eAttributeType operator() (double /*val*/) const { return mapnik::Double; } mapnik::eAttributeType operator() (float /*val*/) const { return mapnik::Double; } mapnik::eAttributeType operator() (bool /*val*/) const { return mapnik::Boolean; } mapnik::eAttributeType operator() (std::string const& /*val*/) const { return mapnik::String; } mapnik::eAttributeType operator() (mapnik::value_unicode_string const& /*val*/) const { return mapnik::String; } mapnik::eAttributeType operator() (mapnik::value_null const& /*val*/) const { return mapnik::String; } }; struct geometry_type_visitor : public mapnik::util::static_visitor { int operator() (mapnik::topojson::point const&) const { return static_cast(mapnik::datasource::Point); } int operator() (mapnik::topojson::multi_point const&) const { return static_cast(mapnik::datasource::Point); } int operator() (mapnik::topojson::linestring const&) const { return static_cast(mapnik::datasource::LineString); } int operator() (mapnik::topojson::multi_linestring const&) const { return static_cast(mapnik::datasource::LineString); } int operator() (mapnik::topojson::polygon const&) const { return static_cast(mapnik::datasource::Polygon); } int operator() (mapnik::topojson::multi_polygon const&) const { return static_cast(mapnik::datasource::Polygon); } int operator() (mapnik::topojson::invalid const&) const { return -1; } }; struct collect_attributes_visitor : public mapnik::util::static_visitor { mapnik::layer_descriptor & desc_; collect_attributes_visitor(mapnik::layer_descriptor & desc): desc_(desc) {} void operator() (mapnik::topojson::invalid const& g) {} template void operator() (GeomType const& g) { if (g.props) { for (auto const& p : *g.props) { desc_.add_descriptor(mapnik::attribute_descriptor(std::get<0>(p), mapnik::util::apply_visitor(attr_value_converter(), std::get<1>(p)))); } } } }; topojson_datasource::topojson_datasource(parameters const& params) : datasource(params), type_(datasource::Vector), desc_(topojson_datasource::name(), *params.get("encoding","utf-8")), filename_(), inline_string_(), extent_(), tr_(new mapnik::transcoder(*params.get("encoding","utf-8"))), #if BOOST_VERSION >= 105600 tree_() #else tree_(16,1) #endif { boost::optional inline_string = params.get("inline"); if (inline_string) { inline_string_ = *inline_string; } else { boost::optional file = params.get("file"); if (!file) throw mapnik::datasource_exception("TopoJSON Plugin: missing parameter"); boost::optional base = params.get("base"); if (base) filename_ = *base + "/" + *file; else filename_ = *file; } if (!inline_string_.empty()) { std::istringstream in(inline_string_); parse_topojson(in); } else { #if defined (_WINDOWS) std::ifstream in(mapnik::utf8_to_utf16(filename_),std::ios_base::in | std::ios_base::binary); #else std::ifstream in(filename_.c_str(),std::ios_base::in | std::ios_base::binary); #endif if (!in.is_open()) { throw mapnik::datasource_exception("TopoJSON Plugin: could not open: '" + filename_ + "'"); } parse_topojson(in); in.close(); } } template void topojson_datasource::parse_topojson(T & stream) { using base_iterator_type = std::istreambuf_iterator; boost::spirit::multi_pass begin = boost::spirit::make_default_multi_pass(base_iterator_type(stream)); boost::spirit::multi_pass end = boost::spirit::make_default_multi_pass(base_iterator_type()); static const mapnik::topojson::topojson_grammar > g; boost::spirit::standard_wide::space_type space; bool result = boost::spirit::qi::phrase_parse(begin, end, g, space, topo_); if (!result) { throw mapnik::datasource_exception("topojson_datasource: Failed parse TopoJSON file '" + filename_ + "'"); } std::size_t count = 0; for (auto const& geom : topo_.geometries) { mapnik::box2d box = mapnik::util::apply_visitor(mapnik::topojson::bounding_box_visitor(topo_), geom); if (box.valid()) { if (count == 0) { extent_ = box; collect_attributes_visitor assessor(desc_); mapnik::util::apply_visitor( std::ref(assessor), geom); } else { extent_.expand_to_include(box); } #if BOOST_VERSION >= 105600 tree_.insert(std::make_pair(box_type(point_type(box.minx(),box.miny()),point_type(box.maxx(),box.maxy())),count++)); #else tree_.insert(box_type(point_type(box.minx(),box.miny()),point_type(box.maxx(),box.maxy())),count++); #endif } } } topojson_datasource::~topojson_datasource() { } const char * topojson_datasource::name() { return "topojson"; } boost::optional topojson_datasource::get_geometry_type() const { boost::optional result; int multi_type = 0; std::size_t num_features = topo_.geometries.size(); for (std::size_t i = 0; i < num_features && i < 5; ++i) { mapnik::topojson::geometry const& geom = topo_.geometries[i]; int type = mapnik::util::apply_visitor(geometry_type_visitor(),geom); if (type > 0) { if (multi_type > 0 && multi_type != type) { result.reset(mapnik::datasource::Collection); return result; } else { result.reset(static_cast(type)); } multi_type = type; } } return result; } mapnik::datasource::datasource_t topojson_datasource::type() const { return type_; } mapnik::box2d topojson_datasource::envelope() const { return extent_; } mapnik::layer_descriptor topojson_datasource::get_descriptor() const { return desc_; } mapnik::featureset_ptr topojson_datasource::features(mapnik::query const& q) const { // if the query box intersects our world extent then query for features mapnik::box2d const& b = q.get_bbox(); if (extent_.intersects(b)) { box_type box(point_type(b.minx(),b.miny()),point_type(b.maxx(),b.maxy())); #if BOOST_VERSION >= 105600 topojson_featureset::array_type index_array; tree_.query(boost::geometry::index::intersects(box),std::back_inserter(index_array)); return std::make_shared(topo_, *tr_, std::move(index_array)); #else index_array_ = tree_.find(box); return std::make_shared(topo_, *tr_, tree_.find(box)); #endif } // otherwise return an empty featureset pointer return mapnik::featureset_ptr(); } mapnik::featureset_ptr topojson_datasource::features_at_point(mapnik::coord2d const& pt, double tol) const { mapnik::box2d query_bbox(pt, pt); query_bbox.pad(tol); mapnik::query q(query_bbox); for (auto const& attr_info : desc_.get_descriptors()) { q.add_property_name(attr_info.get_name()); } return features(q); }