Ist es möglich, die Goode Homolosin-Projektion in QGIS zu verwenden?

11

Ist es möglich, die Goode Homolosin-Projektion in QGIS zu verwenden? Sowohl QGIS als auch GRASS scheinen Probleme mit den Unterbrechungen zu haben und die höheren Breiten zusammenzudrücken. Ich habe versucht, + proj = goode + ellps = sphere + datum = wgs84 als benutzerdefiniertes CRS.

qgis coordinate-system grass homolosine Mattwigway
quelle

Antworten:

0

Es gab schon einmal eine Frage zu unterbrochenen Projektionen: Diskontinuierliche oder unterbrochene Kartenprojektionen?

Die Schlussfolgerung war, dass proj4 (das für die Projektionsbehandlung in QGIS und GRASS verwendet wird) unterbrochene Projektionen nicht unterstützt.

Unterdunkel
quelle
Ich hätte ein Ticket auf dem Proj.4-Tracker für die Unterstützung von Goode Homolosin erwähnen sollen (das erstellt wurde, nachdem diese Frage gestellt wurde): trac.osgeo.org/proj/ticket/106 . Sieht so aus, als hätte ich kein so neues Projekt wie ich dachte; Ich werde versuchen, mein Projekt zu aktualisieren und sehen, wie es geht. Entschuldigung für die doppelte Frage.
Mattwigway
1
Alle Projektionen werden unterbrochen ... dh die Gegendomäne ist ein endlicher Raum, nicht unendlich;)
Luís de Sousa
6

Mit QGIS Lisboa gibt es jetzt eine grundlegende Unterstützung für die Goode Homolosin-Projektion.

Ich habe IGH-Daten von http://edc2.usgs.gov/glcc/fao/data_specifications.php heruntergeladen

und richten Sie diesen Benutzer CRS ein:

+proj=goode +lon_0=-100 +x_0=-11119487.43 +y_0=0 +a=6370997 +b=6370997 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs

für Nordamerika.

lon_0 und x_0 stammen aus Anhang 03 dieses Dokuments: http://proceedings.esri.com/library/userconf/proc98/proceed/TO850/PAP844/P844.HTM

Wenn ich die natürlichen Grenzen der Erdwelt überlagere, bekomme ich dieses Bild:

Geben Sie hier die Bildbeschreibung ein

Passend in Nordamerika (ca. 20 km verlegt) und in anderen Regionen verzerrt.

BEARBEITEN

Ich habe gerade bemerkt, dass proj inzwischen unterbrochenes gutes Homolosin kennt:

+proj=igh +ellps=sphere +towgs84=0,0,0 +lon_0=100w +x_0=-11119487.43

gibt mir dieses Bild der Weltgrenzen der natürlichen Erde und ein 10-Grad-Vektorgitter:

Geben Sie hier die Bildbeschreibung ein

Ich bin mir nicht sicher, welche Parameter außer +proj=ighnoch notwendig sind.

AndreJ
quelle
Mit der obigen Definition erhalte ich eine Pixelverschiebung. Für mich hat die Definition gut +proj=igh +towgs84=0,0,0funktioniert
Thomas
3

QGIS 2.18 unterstützt das Homolosin für Raster nahezu vollständig. Sie müssen sich nur des Problems Nr. 959 bewusst sein : Führen Sie Raster-Warping direkt mit GDAL unter Verwendung der CHECK_WITH_INVERT_PROJParameter durch, und Sie sind gut.

In Bezug auf Vektoren sind die Dinge immer noch komplex. Einige der Fehler, die AndreJ in seiner Antwort zeigt, überwiegen. Die Ursache ist einfach: QGis ist sich nicht bewusst, dass es sich um ein kartografisches Programm handelt, und betrachtet die Co-Domäne einer Projektion als einen endlichen Raum.

Dies hat eine einfache Lösung: Schneiden Sie jede in QGis importierte Vektorschicht mit der Homolosin-Co-Domäne. Vor einiger Zeit stellte ich einen diskreten Co-Domain-Vektor zur Verfügung , der in Intervallen von einem Längen- und Breitengrad definiert wurde.

DOI

Geben Sie hier die Bildbeschreibung ein

Luís de Sousa
quelle