Neuprojektion des Höhendatensatzes für "gekippte Erde"

Ich möchte Karten einer "geneigten Erde" erstellen, wie zum Beispiel das folgende Beispiel, in dem Australien am Nordpol liegt:

Geben Sie hier die Bildbeschreibung ein http://i1048.photobucket.com/albums/s379/laskaris_mods/tilted_earth.jpg

Mit der von mir verwendeten Software Geocart 3 kann ich Karten in Bildformaten wie .jpeg wie diesen neu projizieren. Ich muss jedoch die Höhendatensätze im netcdf-Format (die ETOPO .grd-Daten) auf diese Weise neu projizieren Ich habe genaue Höhendaten für meine "geneigte Erde".

Gibt es eine GIS-Software, die dies kann?

coordinate-system

— Laskaris
quelle

Auch auf Cartotalk gefragt: cartotalk.com/index.php?showtopic=9050, wo daan strebe berichtete, dass Geocart 3 TIFFs unterstützt (Sie müssen die TIFFs CRS in der Software identifizieren).

— Mkennedy

Aber ich muss die tatsächlichen Höhendaten des TIFF neu projizieren, und ich glaube nicht, dass Geocart das tut. Ich würde es auch vorziehen, wenn ich es mit dem netcdf-Format machen könnte.

— Laskaris

Da sich die Zellengrößen / -positionen geändert haben, müssen die Zellenwerte neu abgetastet / neu berechnet werden, aber das sollte passieren. Warum denkst du, ist es nicht? Das Neuprojektieren eines Punkts hat keinen Einfluss auf seinen Elevation / z / H-Wert.

— Mkennedy

Sie haben Recht, es funktioniert für TIFF mit Geocart. Ich habe mich dort geirrt. Vielen Dank! Wenn ich das gleiche mit netcdf machen könnte, wäre es sogar noch besser.

— Laskaris

Dies kann mit Python NetCDF4, Projection4 Library Pyproj und Numpy erfolgen. Vorausgesetzt, Sie kennen die netcdf-Datei CRS und die gekippte Erde CRS, wenn EPSG-Codes verfügbar sind, ist es super einfach. Die Schritte sind

Importieren Sie zuerst die netcdf-Datei mithilfe der netCDF4-Bibliothek in Python.
Fragen Sie die lat, long der netcdf-Datei ab und speichern Sie sie als Array.
Definieren Sie die aktuelle CRS der Netcdf-Datei und die erforderliche CRS für gekippte Erde. Wenn keine CRS vorhanden ist, kann die Projektion definiert werden. Ein Beispiel für die Konvertierung zwischen einer benutzerdefinierten lcc-Projektion in EPSG: 4326 (siehe unten). import pyproj tc={'proj':'lcc','width':'width_meters','height':'height_meters','lat_0':cen_lat,'lon_0':cen_lon,'lat_1':truelat1,'lat_2':truelat2} proj1=pyproj.Proj(tc) proj2_out='+init=EPSG:4326' proj2=pyproj.Proj(proj2_out) lat2,lon2=[],[] for k, l in zip(lon,lat): lat1,lon1=pyproj.transform(proj1,proj2,k,l) lat2.append(lon1) lon2.append(lat1)
Die Arrays lat2 und lon2 enthalten neu projizierte lat long-Werte. Verwenden Sie diese Werte, um eine Netcdf-Datei oder ein tiff mithilfe der Python-Ogr-Bibliothek neu zu erstellen.

Diese Antwort kann einen Starter für die Schritte 1 bis 2 geben. Oder lesen Sie diese erweiterte Einführung für die Projektionskonvertierung mit der Python-Bibliothek pyproj.

— nish
quelle