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:
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.
quelle
Antworten:
Mit den folgenden GDAL-Befehlen können Sie die Daten an die richtige Position verschieben:
(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:
quelle
(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.
quelle
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.
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.
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
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.
quelle