Mehr

Tilecache positioniert meine gdal_translated Karte 20km südlich von wo ich sie haben möchte

Tilecache positioniert meine gdal_translated Karte 20km südlich von wo ich sie haben möchte


Ich habe über gdal_translate eine Geotiff-Datei erstellt:

gdal_translate -a_ullr 547812.2384341886 6834734.336439242 547967.5011478929 6834622.35081697 -a_srs EPSG:900913 pdz3.jpg">

Ich schlage vor zu verwendenEPSG:3857als CRS für das Raster während des Kachelprozesses.

EPSG:900913 ist in aktuellen GDAL und QGIS nicht mehr enthalten, daher wird die spezielle sphärische Mercator-Behandlung nicht angewendet. Die proj-Definition für EPSG:3857 ist

+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1,0 +units=m [email protected] +wktext +no_defs

was ein bisschen widersprüchlich ist zu

SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]]

Dies verursacht typischerweise den Versatz von 20 bis 30 km nach Norden.


20 km falsch in Nord-Süd-Richtung werden sehr wahrscheinlich durch die Vermischung von sphärischen und ellipsoiden Darstellungen der Erde verursacht.

Google Maps und Open Street Maps (OSM) verwenden einen optimierten Mercator (der anscheinend "Web Mercator" heißt), der schneller zu berechnen ist. Die Kachel-Software hinter OSM verwendet wahrscheinlich auch das kugelförmige schnellere Modell, während Ihre Eingabedaten wahrscheinlich eine ellipsoide Erddarstellung verwenden.

Das Problem kann möglicherweise behoben werden, indem ein Parameter semi_minor = 6378137 hinzugefügt wird und alle Daten auf ein einziges gemeinsames Koordinatensystem neu projiziert werden.

Für eine lange Diskussion siehe (etwas schwer zu verdauen) http://alastaira.wordpress.com/2011/01/23/the-google-maps-bing-maps-spherical-mercator-projection/