/***************************************************************************** * * This file is part of Mapnik (c++ mapping toolkit) * * Copyright (C) 2021 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 "process_geojson_file_x3.hpp" #if defined(MAPNIK_MEMORY_MAPPED_FILE) #include MAPNIK_DISABLE_WARNING_PUSH #include #include #include #include MAPNIK_DISABLE_WARNING_POP #include #else #include #endif #include #include #include #include #include #include namespace { constexpr mapnik::json::well_known_names feature_properties[] = {mapnik::json::well_known_names::type, mapnik::json::well_known_names::geometry, mapnik::json::well_known_names::properties}; // sorted constexpr mapnik::json::well_known_names geometry_properties[] = { mapnik::json::well_known_names::type, mapnik::json::well_known_names::coordinates}; // sorted constexpr mapnik::json::well_known_names geometry_collection_properties[] = { mapnik::json::well_known_names::type, mapnik::json::well_known_names::geometries}; // sorted template std::string join(Keys const& keys) { std::string result; bool first = true; for (auto const& key : keys) { if (!first) result += ","; result += "\"" + std::string(mapnik::json::wkn_to_string(key)) + "\""; first = false; } return result; } template bool has_keys(Iterator first1, Iterator last1, Keys const& keys) { auto first2 = std::begin(keys); auto last2 = std::end(keys); for (; first2 != last2; ++first1) { if (first1 == last1 || *first2 < std::get<0>(*first1)) return false; if (!(std::get<0>(*first1) < *first2)) ++first2; } return true; } template bool validate_geojson_feature(mapnik::json::geojson_value& value, Keys const& keys, bool verbose) { if (!value.is()) { if (verbose) std::clog << "Expecting an GeoJSON object" << std::endl; return false; } mapnik::json::geojson_object& feature = mapnik::util::get(value); std::sort(feature.begin(), feature.end(), [](auto const& e0, auto const& e1) { return std::get<0>(e0) < std::get<0>(e1); }); if (!has_keys(feature.begin(), feature.end(), feature_properties)) { if (verbose) std::clog << "Expecting one of " << join(feature_properties) << std::endl; return false; } for (auto& elem : feature) { auto const key = std::get<0>(elem); if (key == mapnik::json::well_known_names::geometry) { auto& geom_value = std::get<1>(elem); if (!geom_value.is()) { if (verbose) std::clog << "\"geometry\": xxx <-- expecting an JSON object here" << std::endl; return false; } auto& geometry = mapnik::util::get(geom_value); // sort by property name std::sort(geometry.begin(), geometry.end(), [](auto const& e0, auto const& e1) { return std::get<0>(e0) < std::get<0>(e1); }); if (!has_keys(geometry.begin(), geometry.end(), geometry_properties) && !has_keys(geometry.begin(), geometry.end(), geometry_collection_properties)) { if (verbose) std::clog << "\"geometry\": xxx <-- expecting one of " << join(geometry_properties) << " or " << join(geometry_collection_properties) << std::endl; return false; } mapnik::geometry::geometry_types geom_type; mapnik::json::positions const* coordinates = nullptr; for (auto& elem2 : geometry) { auto const key2 = std::get<0>(elem2); if (key2 == mapnik::json::well_known_names::type) { auto const& geom_type_value = std::get<1>(elem2); if (!geom_type_value.is()) { if (verbose) std::clog << "\"type\": xxx <-- expecting an GeoJSON geometry type here" << std::endl; return false; } geom_type = mapnik::util::get(geom_type_value); if (geom_type == mapnik::geometry::geometry_types::GeometryCollection) { if (verbose) std::clog << "GeometryCollections are not allowed" << std::endl; ; return false; } } else if (key2 == mapnik::json::well_known_names::coordinates) { auto const& coordinates_value = std::get<1>(elem2); if (!coordinates_value.is()) { if (verbose) std::clog << "\"coordinates\": xxx <-- expecting an GeoJSON positions here" << std::endl; return false; } coordinates = &mapnik::util::get(coordinates_value); } } if (geom_type == mapnik::geometry::geometry_types::Point) { // expecting single position if (!coordinates->is()) { if (verbose) std::clog << "Expecting single position in Point" << std::endl; return false; } } else if (geom_type == mapnik::geometry::geometry_types::LineString) { // expecting if (!coordinates->is()) { if (verbose) std::clog << "Expecting sequence of positions (ring) in LineString" << std::endl; return false; } else { auto const& ring = mapnik::util::get(*coordinates); if (ring.size() < 2) { if (verbose) std::clog << "Expecting at least two coordinates in LineString" << std::endl; return false; } } } else if (geom_type == mapnik::geometry::geometry_types::Polygon) { // expecting if (!coordinates->is()) { if (verbose) std::clog << "Expecting an array of rings in Polygon" << std::endl; return false; } else { auto const& rings = mapnik::util::get(*coordinates); if (rings.size() < 1) { if (verbose) std::clog << "Expecting at least one ring in Polygon" << std::endl; return false; } for (auto const& ring : rings) { if (ring.size() < 4) { if (verbose) std::clog << "Expecting at least four coordinates in Polygon ring" << std::endl; return false; } } } } } } return true; }; using box_type = mapnik::box2d; using boxes_type = std::vector>>; using base_iterator_type = char const*; auto const& geojson_value = mapnik::json::grammar::geojson_value; } // namespace namespace mapnik { namespace detail { template std::pair process_geojson_file_x3(T& boxes, std::string const& filename, bool validate_features, bool verbose) { using box_type = typename T::value_type::first_type; box_type extent; #if defined(MAPNIK_MEMORY_MAPPED_FILE) mapnik::mapped_region_ptr mapped_region; const auto memory = mapnik::mapped_memory_cache::instance().find(filename, true); if (!memory.has_value()) { std::clog << "Error : cannot memory map " << filename << std::endl; return std::make_pair(false, extent); } else { mapped_region = *memory; } base_iterator_type start = reinterpret_cast(mapped_region->get_address()); base_iterator_type end = start + mapped_region->get_size(); #else mapnik::util::file file(filename); if (!file) { std::clog << "Error : cannot open " << filename << std::endl; return std::make_pair(false, extent); } std::string file_buffer; file_buffer.resize(file.size()); auto count = std::fread(&file_buffer[0], file.size(), 1, file.get()); base_iterator_type start = file_buffer.c_str(); base_iterator_type end = (count == 1) ? start + file_buffer.length() : start; #endif base_iterator_type itr = start; // make a copy to preserve `start` iterator state try { mapnik::json::extract_bounding_boxes(itr, end, boxes); } catch (boost::spirit::x3::expectation_failure const& ex) { std::clog << ex.what() << std::endl; std::clog << "Expected: " << ex.which(); std::clog << " Got: \"" << std::string(ex.where(), ex.where() + 200) << '"' << std::endl; return std::make_pair(false, extent); } catch (std::exception const& ex) { std::clog << "mapnik-index (GeoJSON) : could not extract bounding boxes from : '" << filename << "'" << std::endl; return std::make_pair(false, extent); } using namespace boost::spirit; using space_type = mapnik::json::grammar::space_type; auto keys = mapnik::json::get_keys(); #if BOOST_VERSION >= 106700 auto feature_grammar = x3::with(keys)[geojson_value]; #else auto feature_grammar = x3::with(std::ref(keys))[geojson_value]; #endif for (auto const& item : boxes) { if (item.first.valid()) { if (!extent.valid()) extent = item.first; else extent.expand_to_include(item.first); if (validate_features) { base_iterator_type feat_itr = start + item.second.first; base_iterator_type feat_end = feat_itr + item.second.second; mapnik::json::geojson_value feature_value; try { bool result = x3::phrase_parse(feat_itr, feat_end, feature_grammar, space_type(), feature_value); if (!result || feat_itr != feat_end) { if (verbose) std::clog << "Failed to parse: offset=" << item.second.first << " size=" << item.second.second << std::endl; return std::make_pair(false, extent); } } catch (x3::expectation_failure const& ex) { if (verbose) std::clog << ex.what() << std::endl; return std::make_pair(false, extent); } catch (...) { if (verbose) std::clog << "Failed to parse: offset=" << item.second.first << " size=" << item.second.second << std::endl; return std::make_pair(false, extent); } if (!validate_geojson_feature(feature_value, keys, verbose)) { if (verbose) std::clog << "Failed to validate: [" << std::string(start + item.second.first, feat_end) << "]" << std::endl; return std::make_pair(false, extent); } } } else if (validate_features) { if (verbose) std::clog << "Invalid bbox encountered " << item.first << std::endl; return std::make_pair(false, extent); } } return std::make_pair(true, extent); } template std::pair process_geojson_file_x3(boxes_type&, std::string const&, bool, bool); } // namespace detail } // namespace mapnik