2016-04-07 18 views
2

Ich versuche, eine netCDF-Datei als Raster in R zu lesen. Die netCDF-Datei beschreibt die mittlere Jahrestemperatur im Ozean als Funktion von Länge, Breite und Tiefe. Ich interessiere mich für den Oberflächenozean (d. H. Die erste Ebene in der netCDF-Datei) und t_an ist der Name der Variablen in der netCDF-Datei. Deshalb habe ich den folgenden Code verwenden:R: Lesen einer netCDF-Datei als Raster

MyRast <- raster("Temperature.nc", level = 1, varname = "t_an") 

Das gibt mir die folgende Warnung:

Warning message: 
In .getCRSfromGridMap4(atts) : cannot process these parts of the CRS: epsg_code=EPSG:4326 

Wie Sie sehen können, hat die netCDF-Datei CRS EPSG 4326 (oder 84 WGS), noch das Raster das erstellt wird, hat die folgende CRS:

+proj=longlat +lon_0=0 +a=6378137 +rf=298.257232666016 

Irgendwelche Ideen, wie ich in der Datei netCDF mit dem richtigen CRS lesen?

+0

können Sie einen Link zu dem netCDF geben? –

+0

Von diesem Link heruntergeladen https://www.nodc.noaa.gov/cgi-bin/OC5/woa13/woa13.pl Wählen Sie netCDF als Format und 1/4º als verfügbares Raster und die Datei heißt t00_04v2.nc unter die jährliche Überschrift. – Lyngbakr

+0

bist du tot mit netCDF eingestellt? – MikeJewski

Antwort

3

Alles, was Sie tun müssen, ist die Projektion einstellen, nachdem Sie in den Daten lesen:

r <- raster("Temperature.nc", varname = "t_an") 
proj4string(r)=CRS("+init=EPSG:4326") 
+1

oder 'crs (r) <-" + init = EPSG: 4326 "' – RobertH