/usr/include/osmium/geometry/ogr_multipolygon.hpp is in libosmium-dev 0.0~20140910-9a069af-1+b1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 | #ifndef OSMIUM_GEOMETRY_OGR_MULTIPOLYGON_HPP
#define OSMIUM_GEOMETRY_OGR_MULTIPOLYGON_HPP
/*
Copyright 2012 Jochen Topf <jochen@topf.org> and others (see README).
This file is part of Osmium (https://github.com/joto/osmium).
Osmium is free software: you can redistribute it and/or modify it under the
terms of the GNU Lesser General Public License or (at your option) the GNU
General Public License as published by the Free Software Foundation, either
version 3 of the Licenses, or (at your option) any later version.
Osmium 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 and the GNU
General Public License for more details.
You should have received a copy of the Licenses along with Osmium. If not, see
<http://www.gnu.org/licenses/>.
*/
#define OSMIUM_COMPILE_WITH_CFLAGS_OGR `gdal-config --cflags`
#define OSMIUM_LINK_WITH_LIBS_OGR `gdal-config --libs`
#include <ogr_geometry.h>
#include <geos/geom/LineString.h>
#include <geos/geom/CoordinateSequence.h>
#include <geos/geom/Polygon.h>
#include <geos/geom/GeometryCollection.h>
#include <osmium/geometry/multipolygon.hpp>
namespace Osmium {
namespace Geometry {
namespace {
inline void add_ring(OGRPolygon* ogrpolygon, const geos::geom::LineString* geosring) {
OGRLinearRing* ogrring = new OGRLinearRing;
const geos::geom::CoordinateSequence* cs = geosring->getCoordinatesRO();
ogrring->setNumPoints(cs->getSize());
for (size_t i = 0; i < cs->getSize(); ++i) {
ogrring->setPoint(i, cs->getX(i), cs->getY(i));
}
ogrpolygon->addRingDirectly(ogrring);
}
inline OGRPolygon* make_polygon(const geos::geom::Polygon* geospolygon) {
OGRPolygon* ogrpolygon = new OGRPolygon;
add_ring(ogrpolygon, geospolygon->getExteriorRing());
for (size_t i=0; i < geospolygon->getNumInteriorRing(); ++i) {
add_ring(ogrpolygon, geospolygon->getInteriorRingN(i));
}
return ogrpolygon;
}
}
/**
* Create OGR geometry of a MultiPolygon.
*
* Caller takes ownership.
*/
inline OGRMultiPolygon* create_ogr_geometry(const Osmium::Geometry::MultiPolygon& multipolygon) {
OGRMultiPolygon* ogrmp = new OGRMultiPolygon;
const geos::geom::GeometryCollection* geosgeom = dynamic_cast<const geos::geom::GeometryCollection*>(multipolygon.borrow_geos_geometry());
for (geos::geom::GeometryCollection::const_iterator it = geosgeom->begin(); it != geosgeom->end(); ++it) {
OGRPolygon* ogrpolygon = make_polygon(dynamic_cast<const geos::geom::Polygon*>(*it));
OGRErr result = ogrmp->addGeometryDirectly(ogrpolygon);
if (result != OGRERR_NONE) {
// delete ogrpolygon; XXX are we supposed to delete this ourselves?
delete ogrmp;
throw Osmium::Geometry::IllegalGeometry();
}
}
return ogrmp;
}
} // namespace Geometry
} // namespace Osmium
#endif // OSMIUM_GEOMETRY_OGR_MULTIPOLYGON_HPP
|