Ich versuche, eine Rasterdatei im .DEM-Format unter Windows mit dem 'Raster'-Paket in R zu lesen.
Ich bekomme Probleme mit NA-Werten, wenn ich die Daten in R in Windows 7 lade, aber ich habe das Problem auf einem Mac mit OSX Lion nicht. Unter Windows scheinen die NA-Werte nicht richtig gelesen zu werden. Die Frage ist, warum das passiert?
Die verwendete Rasterdatei wurde von USGS mit dem folgenden R-Code heruntergeladen:
download.file('http://edcftp.cr.usgs.gov/pub/data/gtopo30/global/e020n90.tar.gz', 'e020n90.tar.gz')
untar('e020n90.tar.gz')
Dann habe ich das Raster mit dem 'Raster'-Paket in R eingelesen. In OSX Lion und R64 Version 2.13.1 werden die NA-Werte erkannt:
> onMac <- raster('E020N90.DEM')
> onMac
class : RasterLayer
dimensions : 6000, 4800, 28800000 (nrow, ncol, ncell)
resolution : 0.008333333, 0.008333333 (x, y)
extent : 20, 60, 40, 90 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs
values : /Users/Tam/Desktop/E020N90.DEM
min value : -9999
max value : 5483
> summary(values(onMac))
Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
-137 85 148 213 213 5483 13046160
Unter Windows 7 (64 Bit, dieselbe R-Version) werden die Zellenwerte, die NAs sein sollten, in Zahlen konvertiert:
> onWindows <- raster('E020N90.DEM')
> onWindows
class : RasterLayer
dimensions : 6000, 4800, 28800000 (nrow, ncol, ncell)
resolution : 0.008333333, 0.008333333 (x, y)
extent : 20, 60, 40, 90 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0
values : E:/WorldDegreeDays/gsoddata/gtopo/E020N90.DEM
min value : -9999
max value : 5483
> summary(values(onWindows))
Min. 1st Qu. Median Mean 3rd Qu. Max.
1 150 946 27190 55540 65540
Warum enthält das Raster keine NA-Werte, wenn ich es unter Windows lese? Wie könnte ich das umgehen? Ich vermute, es hat mit der Art und Weise zu tun, wie Zahlen gespeichert werden. Viele der NA-Werte werden in 55540 konvertiert.
Info von Windows (nach dem Laden des Rasters):
SessionInfo()
R version 2.13.1 (2011-07-08)
Platform: x86_64-pc-mingw32/x64 (64-bit)
locale:
[1] LC_COLLATE=English_United States.1252
[2] LC_CTYPE=English_United States.1252
[3] LC_MONETARY=English_United States.1252
[4] LC_NUMERIC=C
[5] LC_TIME=English_United States.1252
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] rgdal_0.7-1 raster_1.9-12 sp_0.9-88
loaded via a namespace (and not attached):
[1] grid_2.13.1 lattice_0.19-30
Infos von OSX (nach dem Laden des Rasters):
R version 2.13.1 (2011-07-08)
Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit)
locale:
[1] en_US.UTF-8/en_US.UTF-8/C/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods
[7] base
other attached packages:
[1] rgdal_0.6-33 raster_1.9-12 sp_0.9-88
loaded via a namespace (and not attached):
[1] grid_2.13.1 lattice_0.19-33
sessionInfo()
in Ihren Beitrag aufnehmen?