ฉันพยายามอ่านไฟล์แรสเตอร์ในรูปแบบ. DEM บน windows โดยใช้แพ็คเกจ 'raster' ใน R
ฉันพบปัญหากับค่า NA เมื่อโหลดข้อมูลลงใน R ใน Windows 7 แต่ฉันไม่มีปัญหาใน Mac ที่มี OSX Lion บน windows ค่า NA ดูเหมือนจะอ่านไม่ถูกต้อง คำถามคือทำไมสิ่งนี้เกิดขึ้น?
ไฟล์แรสเตอร์ที่ใช้ถูกดาวน์โหลดจาก USGS ด้วยรหัส R ต่อไปนี้:
download.file('http://edcftp.cr.usgs.gov/pub/data/gtopo30/global/e020n90.tar.gz', 'e020n90.tar.gz')
untar('e020n90.tar.gz')
จากนั้นฉันก็อ่านแรสเตอร์เป็น R โดยใช้แพ็คเกจ 'raster' ใน OSX Lion และ R64 เวอร์ชัน 2.13.1 ค่า NA จะได้รับการยอมรับ:
> 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
แต่ใน Windows 7 (64 บิตรุ่น R เดียวกัน) จะแปลงค่าเซลล์ที่ควรเป็น NA เป็นตัวเลข:
> 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
ทำไมไม่มีค่า NA ในแรสเตอร์เมื่อฉันอ่านบน Windows ฉันจะหลีกเลี่ยงมันได้อย่างไร ฉันเดาว่ามันเกี่ยวข้องกับวิธีจัดเก็บตัวเลขค่า NA จำนวนมากถูกแปลงเป็น 55540
ข้อมูลจาก Windows (หลังจากโหลดแรสเตอร์):
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
ข้อมูลจาก OSX (หลังจากโหลดแรสเตอร์):
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()
ในโพสต์ของคุณ?