netcdf file
Verfasst: Di Mai 15, 2018 1:05 pm
Hallo,
ich arbeite gerade für meine Masterarbeit mit einem netcdf file mit täglichen Niederschlagsdaten. Ich habe noch nie mit diesem Format gearbeitet und bin gerade am verzweifeln... In der Datei sind 3 Dimensionen (Longitude, Latitude, Time) vorhanden, aber ich habe Schwierigkeiten die Dimension Time zu extrahieren...
Wenn ich die Datei in ArcMap importiere zeigt es mir aber auch ein falsches Datum an, deswegen habe ich nun zusätzliche Bedenken, dass etwas mit den Daten nicht stimmt. (Habe die Datei nochmals heruntergeladen, um Fehler beim Download ausschließen zu können)
Wenn ich wie untenstehend vorgehe, dann bekomme ich bei time einen Zahlencode und nicht das wirkliche Datum.
Wäre toll, wenn mir jemand helfen könnte
Vielen Dank schon mal!
cpath <- "/E-OBS/scPDSI.cru.3.24.bams.2016.GLOBAL.1901.2015_with.times.nc/"
ncname <- "scPDSI.cru.3.24.bams.2016.GLOBAL.1901.2015_with.times"
ncfname <- paste(ncpath, ncname, ".nc", sep="")
dname <- "rr" # note: rainfall
# open a NetCDF file
ncin <- nc_open(ncfname)
print(ncin)
attributes(ncin$var)
dat <- ncvar_get(ncin, attributes(ncin$var)$names[1])
dat
# get longitude and latitude
lon <- ncvar_get(ncin,"longitude")
nlon <- dim(lon)
head(lon)
lat <- ncvar_get(ncin,"latitude")
nlat <- dim(lat)
head(lat)
print(c(nlon,nlat))
# get time
time <- ncvar_get(ncin,"time")
time
tunits <- ncatt_get(ncin,"time","units")
nt <- dim(time)
nt
ich arbeite gerade für meine Masterarbeit mit einem netcdf file mit täglichen Niederschlagsdaten. Ich habe noch nie mit diesem Format gearbeitet und bin gerade am verzweifeln... In der Datei sind 3 Dimensionen (Longitude, Latitude, Time) vorhanden, aber ich habe Schwierigkeiten die Dimension Time zu extrahieren...
Wenn ich die Datei in ArcMap importiere zeigt es mir aber auch ein falsches Datum an, deswegen habe ich nun zusätzliche Bedenken, dass etwas mit den Daten nicht stimmt. (Habe die Datei nochmals heruntergeladen, um Fehler beim Download ausschließen zu können)
Wenn ich wie untenstehend vorgehe, dann bekomme ich bei time einen Zahlencode und nicht das wirkliche Datum.
Wäre toll, wenn mir jemand helfen könnte
Vielen Dank schon mal!
cpath <- "/E-OBS/scPDSI.cru.3.24.bams.2016.GLOBAL.1901.2015_with.times.nc/"
ncname <- "scPDSI.cru.3.24.bams.2016.GLOBAL.1901.2015_with.times"
ncfname <- paste(ncpath, ncname, ".nc", sep="")
dname <- "rr" # note: rainfall
# open a NetCDF file
ncin <- nc_open(ncfname)
print(ncin)
attributes(ncin$var)
dat <- ncvar_get(ncin, attributes(ncin$var)$names[1])
dat
# get longitude and latitude
lon <- ncvar_get(ncin,"longitude")
nlon <- dim(lon)
head(lon)
lat <- ncvar_get(ncin,"latitude")
nlat <- dim(lat)
head(lat)
print(c(nlon,nlat))
# get time
time <- ncvar_get(ncin,"time")
time
tunits <- ncatt_get(ncin,"time","units")
nt <- dim(time)
nt