/***************************************************************************** * * This file is part of Mapnik (c++ mapping toolkit) * * Copyright (C) 2016 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 * *****************************************************************************/ // mapnik #include "geojson_index_featureset.hpp" #include #include #include #include #include #include #include #include // stl #include #include #include geojson_index_featureset::geojson_index_featureset(std::string const& filename, mapnik::filter_in_box const& filter) : #if defined(MAPNIK_MEMORY_MAPPED_FILE) // #elif defined _WINDOWS file_(_wfopen(mapnik::utf8_to_utf16(filename).c_str(), L"rb"), std::fclose), #else file_(std::fopen(filename.c_str(),"rb"), std::fclose), #endif ctx_(std::make_shared()) { #if defined (MAPNIK_MEMORY_MAPPED_FILE) boost::optional memory = mapnik::mapped_memory_cache::instance().find(filename, true); if (memory) { mapped_region_ = *memory; } else { throw std::runtime_error("could not create file mapping for " + filename); } #else if (!file_) throw std::runtime_error("Can't open " + filename); #endif std::string indexname = filename + ".index"; std::ifstream index(indexname.c_str(), std::ios::binary); if (!index) throw mapnik::datasource_exception("GeoJSON Plugin: can't open index file " + indexname); mapnik::util::spatial_index::query(filter, index, positions_); std::sort(positions_.begin(), positions_.end(), [](value_type const& lhs, value_type const& rhs) { return lhs.first < rhs.first;}); itr_ = positions_.begin(); } geojson_index_featureset::~geojson_index_featureset() {} mapnik::feature_ptr geojson_index_featureset::next() { while( itr_ != positions_.end()) { auto pos = *itr_++; #if defined(MAPNIK_MEMORY_MAPPED_FILE) char const* start = (char const*)mapped_region_->get_address() + pos.first; char const* end = start + pos.second; #else std::fseek(file_.get(), pos.first, SEEK_SET); std::vector record; record.resize(pos.second); std::fread(record.data(), pos.second, 1, file_.get()); auto const* start = record.data(); auto const* end = start + record.size(); #endif static const mapnik::transcoder tr("utf8"); mapnik::feature_ptr feature(mapnik::feature_factory::create(ctx_, feature_id_++)); using mapnik::json::grammar::iterator_type; mapnik::json::parse_feature(start, end, *feature, tr); // throw on failure // skip empty geometries if (mapnik::geometry::is_empty(feature->get_geometry())) continue; return feature; } return mapnik::feature_ptr(); }