Längsschnittprobleme mit NetCDF-Daten

8

Ich suche Rat bei der Arbeit mit Weltdatensätzen, bei denen das GCS einen positiven Bereich aufweist (0 bis 360 Grad). Die Daten, mit denen ich arbeite, sind ozeanografische NetCDF- Daten und haben wie erwähnt positive Koordinatenwerte. Wenn es in einem normalen GCS WGS84 in ArcGIS 10 von ESRI angezeigt wird, wird es von anderen Daten im Bereich von -180 bis 180 versetzt. Wenn ich es neu projiziere, treten Lücken um den Nullmeridian auf (hauptsächlich aufgrund seiner Herkunft, manchmal bis zu diesem Punkt) aus als 25 Grad W). Mein Gedanke ist es, ein benutzerdefiniertes Projektions- / Koordinatensystem zu erstellen, das die positiven Werte aufweist, und zu hoffen, dass andere Weltdatensätze im laufenden Betrieb problemlos neu projiziert werden. Hat jemand Probleme damit oder andere Lösungen? Hoffe das macht Sinn. Ich freue mich über jede Hilfe im Voraus.

** Update ** Dies ist ein Screenshot eines Länderdatensatzes mit einem Raster oben, um das normale CRS von GCS zu zeigen. Wenn Sie die NetCDF-Daten so wie sie sind importieren, werden sie eingegeben, aber die Werte beginnen erst 20,5 Grad östlich des Nullmeridians.

Normales GCS http://grafa.co/rnd/img/GCS_Normal.png

Wenn ich mich dafür entscheide, das CRS des NetCDF zu verwenden (das eigentlich dasselbe ist), werden die Weltdaten im laufenden Betrieb neu projiziert, wenn sie definiert sind. Beachten Sie, dass das Raster undefiniert ist und nicht neu projiziert wird. Art von neu projiziertem GCS http://grafa.co/rnd/img/GCS_Other.png

Wenn ich dann versuche, alles auf ein normales GCS mit negativen Werten zu projizieren, ist es so, als würde man eine flache Karte der Daten auf der ganzen Welt einwickeln und sie verschwindet, wenn sie auf den Nullmeridian trifft. Art von neu projiziertem GCS http://grafa.co/rnd/img/reprojected.png

Jetzt weiß ich, dass es keine Werte von 0 bis 20,5 gibt, wie in den Metadaten angegeben. Aber warum können die Werte aus dem Rest der Daten nicht in der Lücke angezeigt werden? Ich habe sogar erfolglos versucht, die Raster-Tools zu verschieben.

geografa
quelle
geografa - Wenn Sie eine benutzerdefinierte Projektion definieren, können Sie das Dokument nicht auf ein anderes crs (einen der anderen Weltdatensätze) festlegen und müssen diese nicht "on-the-fly" ausführen. Aber das netcdf würde?
Brad Nesom
Bitte posten Sie irgendwo eine Kopie des NetCDF. Es gibt vielleicht eine Möglichkeit, das zu tun, was Sie wollen, aber es wäre einfacher, es Ihnen anhand eines realen Beispiels zu beschreiben. Im Wesentlichen könnte es möglich sein, etwas wie im Abschnitt "Alles zusammenbinden" dieser Präsentation zu tun: download.osgeo.org/gdal/presentations/… . ODER es ist möglicherweise möglich, einfach eine neue Weltdatei mit Standardkoordinatenbereichen zu erstellen und die in der NetCDF integrierte Georeferenzierung zu überladen. In jedem Fall würde eine Datenprobe helfen.
Lagerratrobe
Ich habe festgestellt, dass ich das Raster neu projizieren kann. Ich musste ihm nur den richtigen Versatz geben. Ich werde in den nächsten ein oder zwei Tagen eine prägnantere Antwort veröffentlichen. Daten finden Sie unter data-portal.ecmwf.int/data/d/edit/personal/temporary/netcdf29 .
Geografa
Ich wäre daran interessiert, hier den richtigen Versatz einzustellen. Eine Sache zu beachten ist, dass, glaube ich, 0 und 360 gleich sind.
Haarige

Antworten:

1

OK, ich denke, das Hauptproblem hier ist, dass ESRI die Mitte des Gitters verwendet, um die Position zu berechnen, und wenn ich die Daten geschrieben habe, verwende ich den unteren linken Wert der Zelle, um die Daten zu berechnen. Ich habe dies gelöst, indem ich meine netCDF-Datei mit einem X-Raster-Ursprung in der Mitte der ersten Zelle erstellt und dann jede Zelle als Ausgangspunkt berechnet habe. Schreiben Sie mir im Chat eine Nachricht, wenn Sie dies weiter diskutieren möchten, da es jetzt funktioniert.

Behaart
quelle
0

Dies ist jetzt ein bekannter Fehler, da wir das gleiche Problem festgestellt haben. Mit ein bisschen zu viel und zu wenig, was dazu führte, dass ich ihnen Details zum Panopol schickte, um die netCDF-Daten ohne die fehlende Zeile anzuzeigen.

Interessanterweise ist dies nicht nur auf das netCDF beschränkt. Wenn Sie Punkte mithilfe der Gitter in der netCDF erstellen und dann die Punkte aggregieren, um Polygone zu erstellen, wird kein Polygon über dem Nullmeridian erstellt, sondern es handelt sich um separate Punkte.

Der Fehler ist NIM070286: NetCDF Raster Layer zeigt keine Daten an, die sich entlang des Nullmeridians befinden

Ich habe einen Hotfix von ESRI Inc. angefordert. Ich werde dies als Lieblingsfrage beibehalten und alles aktualisieren, sobald ich mehr davon höre. Ich bin bald im Urlaub, also hoffentlich wird dies in den nächsten 10 Tagen geschehen.

Behaart
quelle