Wo wird dieser Rasterversatz in die QGIS-Verarbeitungs-Toolbox aufgenommen?

8

Ich verwende die (fantastische!) Verarbeitungs-Toolbox in QGIS 2.0, um einige NoData-Werte in einem Raster von -9999 auf 0 neu zu klassifizieren, aber irgendwann wird ein Versatz von etwa 0,2 Grad in die Ausgabeebene eingeführt.

Ich verwende das Tool "Reclassify (Simple)" (ich glaube, es ist Teil von SAGA) mit einigen ziemlich grundlegenden Parametern:

Parameter neu klassifizieren

. . . und wo sich das Eingabe-Raster im Vergleich zu diesem Hex-Raster korrekt befand:

Eingabe-Raster

. . . Das Ausgabe-Raster wird von zwei Hex-Zellen oder etwa 0,2 Grad nach Südosten gestoßen:

Geben Sie hier die Bildbeschreibung ein

Der Eingabeursprung ist -180, 83.623. . . Der Ausgabeursprung ist jedoch -179.05,83.05. Normalerweise würde ich annehmen, dass dies das Ergebnis einer fehlerhaften WGS84-Konvertierung ist, aber die Eingabe, Ausgabe und das Projekt sind alle EPSG: 4326.

Woher kommt dieser Ursprungsversatz und lohnt es sich, ein SAGA-Ticket einzureichen, oder ist es nur etwas, das ich übersehen habe?

Aktualisieren Das Prozessprotokoll zeigt an, dass eine Transformation für den Eingabeursprung durchgeführt wird:

Transformation x' = -179.999989 + x * 0.100000 + y * 0.000000
Transformation y' = 83.623125 + x * 0.000000 + y * -0.100000

und ein paar Zeilen später werden die neuen Parameter angezeigt, die jetzt falsch sind:

Grid system: 0.1; 3600x 1735y; -179.949989x -89.826875y
Bill Morris
quelle
1
Sie sollten das native SAGA-Tool über die GUI oder die CLI testen und prüfen, ob Sie das gleiche Ergebnis erhalten.
Giovanni Manghi
1
Zu viele Sackgassen aufgrund dieser Installationsdokumente auf meinem Computer - Ich bin nicht daran interessiert, meine aktuelle GDAL / PROJ-Konfiguration zu löschen, nur um SAGA nativ zu testen.
Bill Morris
1
Wenn Sie SAGA in QGIS verwenden, ist SAGA bereits auf Ihrem Computer installiert.
Giovanni Manghi
Nach mehr als 2 Stunden Durchstöbern der Dokumente habe ich nicht die geringste Ahnung, wie man die GUI in OSX startet. Die Dokumentation erinnert in all den schlechten Dingen an R. Keine Variation von "saga" (-gis, -gui usw.) erzeugt irgendetwas in der Befehlszeile, und das Kompilieren aus dem Quellcode liefert völlig nicht hilfreiche Fehler. Dem barmherzigen Himmel (also den Kernentwicklern) sei Dank, dass dieses Chaos jetzt für Menschen in einer QGIS-Seitenleiste zugänglich ist. Ich hoffe aufrichtig, dass SAGA nicht den Ursprung meiner Rasterausgabe ändert, da ich ihn niemals debuggen kann.
Bill Morris
Können Sie die Saga saga_guistattdessen mit einem Befehl starten ?
SaultDon

Antworten:

1

Ihr Versatz beträgt 1/2 Pixel (0,1 Grad).

Längengrad:

-179,999989 - -179,949989 = 0,05

Der Breitengrad ist weniger offensichtlich, da ein Raster von links oben und das andere von links unten referenziert wird:

83,623125 - (1735 * 0,1) - -89,826875 = 0,05

Ich vermute, dies hängt mit dem Problem RasterPixelIsPoint oder RasterPixelIsArea zusammen:

http://www.remotesensing.org/geotiff/faq.html#PixelIsPoint

Verwenden Sie gdal 1.7 oder 1.8?

Marc Pfister
quelle
Systemweit habe ich GDAL 1.10.1; Aus irgendeinem Grund kann ich die QGIS-Python-Konsole nicht dazu bringen, die dort verwendete Version
Bill Morris,
es sieht weit mehr als 1/2 Pixel auf dem Bild aus.
Radouxju
Es ist, aber die Pixelgröße ist auch aus. Hier ist die gdalinfo für die Eingabe- und Ausgabe-Tiffs
Bill Morris
Ah, wenn man die gdalinfo sieht, sieht es so aus, als ob die Verschiebung .95 und -0.5731265 ist. Die obige Verschiebung um 1/2 Pixel ist wahrscheinlich eine korrekte Umwandlung von Bereich zu Punkt. Ihr ursprüngliches Raster hat einige merkwürdige Grenzen, vielleicht vermasseln sie die Konvertierung.
Marc Pfister
Guter Punkt. Seltsam, dass es bis zu einem gewissen Grad am Rande der Welt hängt. Das Raster ist das Ergebnis einer Konvertierung vom Vektor. Lassen Sie mich sehen, ob ich diese Konvertierung mit den richtigen Grenzen wiederholen kann.
Bill Morris