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.

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.


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.


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

Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.