Tôi đang cố đọc tệp raster ở định dạng .DEM trên windows bằng gói 'raster' trong R.
Tôi gặp vấn đề với các giá trị NA, khi tải dữ liệu vào R trong Windows 7, nhưng tôi không gặp vấn đề gì với máy Mac với OSX Lion. Trên cửa sổ, các giá trị NA dường như không được đọc chính xác. Câu hỏi là tại sao điều này xảy ra?
Tệp raster được sử dụng đã được tải xuống từ USGS với mã R sau:
download.file('http://edcftp.cr.usgs.gov/pub/data/gtopo30/global/e020n90.tar.gz', 'e020n90.tar.gz')
untar('e020n90.tar.gz')
Sau đó, tôi đọc raster vào R bằng cách sử dụng gói 'raster'. Trong OSX Lion và R64 phiên bản 2.13.1, các giá trị NA được nhận dạng:
> 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
Nhưng trên Windows 7 (64Bit, cùng phiên bản R), nó chuyển đổi các giá trị ô nên là NA thành số:
> 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
Tại sao không có giá trị NA trong raster khi tôi đọc nó trên Windows? Làm thế nào tôi có thể làm việc xung quanh nó? Tôi đoán là nó có liên quan đến cách lưu trữ số, rất nhiều giá trị NA được chuyển đổi thành 55540.
Thông tin từ Windows (sau khi tải raster):
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
Thông tin từ OSX (sau khi tải raster):
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()
trong bài viết của bạn?