Differential D2007 Diff 4944 tools/osm-simplify/ShpCoastlineProcessor.cpp

Changeset View

Standalone View

tools/osm-simplify/ShpCoastlineProcessor.cpp

This file was added.

1		//
2		// This file is part of the Marble Virtual Globe.
3		//
4		// This program is free software licensed under the GNU LGPL. You can
5		// find a copy of this license in LICENSE.txt in the top directory of
6		// the source code.
7		//
8		// Copyright 2016 David Kolozsvari <freedawson@gmail.com>
9		//
10
11		#include "ShpCoastlineProcessor.h"
12
13		#include "BaseClipper.h"
14
15		#include "GeoDataPlacemark.h"
16		#include "OsmPlacemarkData.h"
17
18		#include <QDebug>
19
20		ShpCoastlineProcessor::ShpCoastlineProcessor(GeoDataDocument* document) :
21		PlacemarkFilter(document)
22		{
23
24		}
25
26		void ShpCoastlineProcessor::process()
27		{
28		OsmPlacemarkData marbleLand;
29		marbleLand.addTag("marble_land","landmass");
30
31		foreach (GeoDataObject* object, m_objects) {
32		GeoDataPlacemark* placemark = static_cast<GeoDataPlacemark*>(object);
33		if(placemark->geometry()->nodeType() == GeoDataTypes::GeoDataPolygonType) {
34		placemark->setOsmData(marbleLand);
35		}
36		}
37		}
38
39		GeoDataDocument *ShpCoastlineProcessor::cutToTiles(unsigned int zoomLevel, unsigned int tileX, unsigned int tileY)
40		{
41		unsigned int N = pow(2, zoomLevel);
42
43		GeoDataDocument* tile = new GeoDataDocument();
44		QString tileName = QString("%1/%2/%3").arg(zoomLevel).arg(tileX).arg(tileY);
45		tile->setName(tileName);
46
47		GeoDataLatLonBox tileBoundary;
48		qreal north = BaseClipper::tileY2lat(tileY, N);
49		qreal south = BaseClipper::tileY2lat(tileY+1, N);
50		qreal west = BaseClipper::tileX2lon(tileX, N);
51		qreal east = BaseClipper::tileX2lon(tileX+1, N);
52
53		// qDebug() << tileName;
54		// qDebug() << "west: " << westRAD2DEG << "\t east: " << eastRAD2DEG;
55		// qDebug() << "north: " << northRAD2DEG << "\t south: " << southRAD2DEG;
56		// qDebug() << "\n";
57
58		tileBoundary.setBoundaries(north, south, east, west);
59
60		foreach (GeoDataObject* object, m_objects) {
61		GeoDataPlacemark* placemark = static_cast<GeoDataPlacemark*>(object);
62
63		if(placemark->geometry()->nodeType() == GeoDataTypes::GeoDataPolygonType) {
64		GeoDataPolygon* marblePolygon = static_cast<GeoDataPolygon*>(placemark->geometry());
65
66		if(tileBoundary.intersects(marblePolygon->latLonAltBox())) {
67		BaseClipper clipper;
68		clipper.initClipRect(tileBoundary);
69
70		QVector<QPolygonF> clippedPolygons;
71
72		QPolygonF outerBoundary = *BaseClipper::linearRing2Qpolygon(marblePolygon->outerBoundary());
73		qDebug() << "Size before:" << outerBoundary.size();
74
75		qDebug() << "Clipping...";
76		clipper.clipPolyObject(outerBoundary, clippedPolygons, true);
77		qDebug() << "done.";
78
79		qDebug() << "Number of polygons after clipping: " << clippedPolygons.size();
80		qDebug() << "Size(s) after:";
81		foreach(const QPolygonF& polygon, clippedPolygons) {
82		qDebug() << polygon.size();
83
84		// Dirty fix to close the polygons.
85		GeoDataLinearRing outerBoundary = *BaseClipper::qPolygon2linearRing(polygon);
	nienhueserUnsubmitted Not Done Manual polygon closing should not be needed anymore with http://commits.kde.org/marble/39fb2c303bc182d1d513c453c6233e5d1a341635 nienhueser: Manual polygon closing should not be needed anymore with http://commits.kde.
86		if(outerBoundary.first() != outerBoundary.last()) {
87		outerBoundary.append(outerBoundary.first());
88		}
89
90		GeoDataPolygon* newMarblePolygon = new GeoDataPolygon();
91		newMarblePolygon->setOuterBoundary(outerBoundary);
92
93		GeoDataPlacemark* newPlacemark = new GeoDataPlacemark();
94		newPlacemark->setGeometry(newMarblePolygon);
95		newPlacemark->setVisualCategory(GeoDataFeature::Landmass);
96		newPlacemark->setOsmData(placemark->osmData());
97
98		tile->append(newPlacemark);
99		}
100		qDebug() << "";
101		}
102		}
103		}
104
105
106		return tile;
107		}