Welches Koordinatensystem ist dieses netCDF-Raster?

10

Ich habe eine netCDF-Rasterdatei erhalten, konnte jedoch keine Metadaten abrufen, um den Namen des Koordinatensystems abzurufen, auf dem das Raster aufgebaut wurde. In das Raster selbst ist kein Koordinatensystem eingebettet. Ich dachte, es wäre nur ein WGS84, und auf den ersten Blick sah es so aus, aber bei weiteren Untersuchungen in ArcMap sah ich, dass es sich um ein eher ungewöhnliches System handelt. So wird es angezeigt: Geben Sie hier die Bildbeschreibung ein

Das orangefarbene Raster ist ein normales Raster in WGS84, das ich hier zu Vergleichszwecken eingefügt habe. Das lila ist das Raster mit dem unbekannten Koordinatensystem. Haben Sie eine Ahnung, was das sein könnte?

Einige Updates: Hier ist das netCDF-Raster: https://www.dropbox.com/s/nottbl9yt6dwss6/sic_average_nclimate.nc?dl=0 Ich konnte auch einige Metadaten vom Bildanbieter abrufen:

netcdf sic_average_nclimate {
dimensions:
    nlon = 361 ;
    nlat = 90 ;
    nseas = 4 ;
variables:
    float SIC_Change(nlat, nlon) ;
        SIC_Change:Title = "Gridded Multi-Model Ensemble Mean Annual Mean Change in Ice Concentration 21C-20C" ;
    float SIC_Season_Change(nseas, nlat, nlon) ;
        SIC_Season_Change:Title = "Gridded Multi-Model Ensemble Mean Seasonal Mean Change in Ice Concentration 21C-20C" ;
    float SIC_Change_STD(nlat, nlon) ;
        SIC_Change_STD:Title = "Gridded Multi-Model Standard Deviation of the Annual Mean Change in Ice Concentration 21C-20C" ;
    float SIC_Season_Change_STD(nseas, nlat, nlon) ;
        SIC_Season_Change_STD:Title = "Gridded Multi-Model Standard Deviation of the Seasonal Mean Change in Ice Concentration 21C-20C" ;
    float LAT(nlat) ;
        LAT:Title = "Latitude" ;
    float LON(nlon) ;
        LON:Title = "Longitude" ;

// global attributes:
        :Title = "Ice Concentration metrics for Model subset as in Figure 1 of NCLIMATE paper" ;

Sie zeigen die Grenze lat-lang, aber anscheinend keine Informationen bezüglich des Koordinatensystems.

multigoodverse
quelle
1
Die Antarktis befindet sich normalerweise in polarer stereografischer Projektion, sodass alles vom Südpol nach Norden verläuft. nsidc.org/data/polar_stereo/ps_grids.html
Mapperz
1
Seltsam, sieht aus wie 0-360, aber warum zeichnet es auf der Nordhalbkugel?
Mkennedy
1
Wenn die gemeldeten Extents im NetCDF-Header wie lat / lon aussehen, aber lat: 0 bis ?, Lon: 0 bis 360. Können Sie sie bis lat: -90 bis (? - 90) und lon: -180 bis + hacken 180?
Mkennedy
@Mapperz & mkennedy, ich habe versucht, ihm eine polare stereografische Projektion zuzuweisen, aber es hat nicht funktioniert. Ich kann es wahrscheinlich vertikal nach unten spiegeln, aber die rechte Hälfte wird immer noch über 180E Längen liegen. Ich habe die netCDF mit dem Editor geöffnet, konnte aber keine Informationen zu lat, long finden. Ich füge das Raster der Frage bei, wenn Sie neugierig sind, es auf Ihrer Seite zu öffnen.
Multigoodverse
Dieser Link funktioniert nicht - Dropbox sagt "Der Eigentümer hat keinen autorisierten Zugriff".
BradHards

Antworten:

6

Mit den folgenden GDAL-Befehlen können Sie die Daten an die richtige Position verschieben:

gdal_translate -a_srs EPSG:4326 -a_ullr 0 0 360 -90 NETCDF:"sic_average_nclimate.nc":SIC_Change change360.tif
gdalwarp -t_srs WGS84 change360.tif change180.tif  -wo SOURCE_EXTRA=1000 --config CENTER_LONG 0

(mit ein wenig Hilfe von Frank Warmerdam: Wie man ein Raster von 0 360 auf -180 180 projiziert, indem man 180 Meridian schneidet )

und das Ergebnis sieht in einer stereografischen Projektion so aus:

Geben Sie hier die Bildbeschreibung ein

AndreJ
quelle
1

(Dies ist keine Antwort, sondern ein Prozessvorschlag, der zu lang ist, um einen Kommentar abzugeben.)

Sie können es in QGIS mit dem GeoReferencer-Tool georeferenzieren, da Sie anscheinend in der Lage sind, gemeinsame Punkte auf beiden Rastern zu identifizieren. Dies würde ein GDAL-Skript mit einer Reihe von GCP-Punkten und Transformationsparametern generieren. Sie können dieses Skript mit einem Beispiel dieser Punkte zu Ihrer Frage hinzufügen, da es möglicherweise den Unterschied zwischen den beiden Projektionen erklärt.

Im schlimmsten Fall hätten Sie eine grobe Schätzung der Werte basierend auf dem neu positionierten Raster, die Ihnen bei Ihrem Ziel helfen könnten.

djq
quelle
+ 1 und ja, Georeferenzierung ist möglicherweise die letzte Kugel, aber der Prozess würde einige Fehler ergeben, da es schwierig ist, Kontrollpunkte zuzuordnen. Außerdem müsste das netCDF-Raster aufgeteilt werden, da es nicht denselben Nullmeridian hat. Zum Beispiel ist der Schwanz der Antarktis rechts, während das orangefarbene Raster ihn links zeigt.
Multigoodverse
0

Die untere linke Ecke liegt bei (0,0) Koordinaten, aber für mich ist (0,0) die Koordinate der oberen rechten Ecke in Ihrem Fall, und ich kenne kein geografisches Koordinatensystem mit einem Ursprung bei -90.

Während Sie in ArcGIS arbeiten, können Sie die Georeferenzierungssymbolleiste verwenden , um Ihre Daten an den richtigen Ort zu verschieben. Bewegen Sie die obere rechte Seite nach (0,0) und die untere linke nach (-360, -90). Dies sollte eine genaue Transformation ergeben, aber Sie hätten immer noch ein Problem mit der Datenlinie.

Geben Sie hier die Bildbeschreibung ein

Sie können also diesen Link ausprobieren , der Ihnen beim Erstellen eines benutzerdefinierten geografischen Koordinatensystems helfen soll. Um es zusammenzufassen, erstellen Sie in Ihrem Referenzsystem ein benutzerdefiniertes WGS 84 mit dem Nullmeridian bei -180 und anschließend eine benutzerdefinierte Transformation mit Längengradrotation.

Geben Sie hier die Bildbeschreibung ein

Georeferenzierung in diesem neuen System würde darin bestehen, den Ursprung auf (-180, -90) zu setzen. Ein Kontrollpunkt reicht aus, da Sie nur eine Übersetzung benötigen. Wählen Sie das untere linke Pixel sorgfältig aus, klicken Sie mit der rechten Maustaste und fügen Sie die genauen XY-Werte hinzu. Dann sollte es so aussehen

Geben Sie hier die Bildbeschreibung ein

Als letzte Bemerkung gibt es keine Gewissheit, dass das Koordinatensystem auf WGS 84 basiert. Es handelt sich jedoch mit Sicherheit um ein geografisches Koordinatensystem in Grad, und die Auflösung Ihres Rasters ist wahrscheinlich gröber als der potenzielle Fehler aufgrund einer falschen Bezugsauswahl.

Radouxju
quelle