Manipulieren von azimutalen äquidistanten Projektionen in QGIS

8

Ich habe ein Shapefile der Weltländer, das als Azimuthal Equidistant mit Chicago, USA als Zentrum projiziert wird. Ich möchte die Projektion so ändern, dass Edinburgh, Großbritannien, das Zentrum ist. Man würde denken, dass dies so einfach sein sollte wie das Ändern von 2 Zahlen in den Koordinatensystemparametern ...

Das Problem, das ich zu haben scheine, ist, dass QGIS es schwierig, wenn nicht unmöglich macht, den vollständigen proj4-Code der Koordinatensystemparameter eines bestimmten Shapefiles anzuzeigen.

Folgendes habe ich bisher versucht:

  • Shapefile der als Azimuthal Equidistant projizierten Länder der Welt mit Chicago als Zentrum hinzugefügt
  • Geöffnete Eigenschaften, Registerkarte Metadaten, kopieren Sie den proj4-Parametercode
  • Ging zu Einstellungen, benutzerdefinierte Projektion, Einfügen in proj4-Code, änderte das lat / lon von dem von Chicago zu dem von Edinburgh
  • Gespeichertes Shapefile der Weltländer, spezifizierte meine neue benutzerdefinierte Projektion, fügte der Karte eine neue Datei hinzu
  • Meine neue Datei sieht aus wie eine Karte von Pangaea, auf der alle Kontinente zusammengewürfelt sind.

Ich denke, es muss andere Parameter geben als die, die im proj4-Code der Registerkarte Metadaten angezeigt werden. Das einzige, was ich im proj4-Code ändere, ist die Mitte der Karte, aber die Änderung, die ich sehe, ist viel mehr als das. Meine Beobachtung legt nahe, dass sich mehr als eine unabhängige Variable ändert.

Gerne sende ich mein Shapefile an alle, die damit spielen möchten.

Vielen Dank,

Daniel Wolf Umwelt- / Geodaten-Enthusiast

Dan Wolf
quelle

Antworten:

8

Verwenden einer einfachen Proj-Zeichenfolge:

+proj=aeqd +lat_0=51 +lon_0=7

Ich bekomme dieses Bild von den Weltgrenzen der natürlichen Erde und einem 10-Grad-Gitter, das sich von -150 ° bis 150 ° O erstrecken sollte:

aeqd auf einem WGS84-Ellipsoid

Es sieht so aus, als würde alles außerhalb von lon_0 +/- 90 ° abgeschnitten und auch Teile der südlichen Hemisphäre werden gebogen.


BEARBEITEN

Laut diesem Ticket: http://trac.osgeo.org/proj/ticket/211 proj hat ein Problem mit der Neuprojektion von aeqd-Daten auf einem Ellipsoid.

Eine einfache Lösung: Verwenden Sie eine Kugel (und nicht das Standardellipsoid), indem Sie einen Radius definieren:

+proj=aeqd  +R=6371000 +lat_0=51 +lon_0=7

Dies gibt mir das folgende Bild:

aeqd auf einer Kugel

Gitterlinien werden verdichtet, um schöne Kurven zu zeigen.

AndreJ
quelle