मैं R में 'raster' पैकेज का उपयोग करके विंडोज़ पर .DEM प्रारूप में एक रैस्टर फ़ाइल पढ़ने की कोशिश कर रहा हूँ।
विंडोज 7 में आर में डेटा लोड करते समय मुझे एनए मानों के साथ समस्याएं मिलती हैं, लेकिन मुझे ओएसएक्स लायन के साथ मैक पर समस्या नहीं है। खिड़कियों पर, एनए मान सही ढंग से पढ़ा नहीं लगता है। सवाल यह है कि ऐसा क्यों होता है?
उपयोग की गई रेखापुंज फ़ाइल यूएसजीएस से निम्नलिखित आर कोड के साथ डाउनलोड की गई थी:
download.file('http://edcftp.cr.usgs.gov/pub/data/gtopo30/global/e020n90.tar.gz', 'e020n90.tar.gz')
untar('e020n90.tar.gz')
तब मैंने 'रेखापुंज' पैकेज का उपयोग करके आर को रेखापुंज में पढ़ा। ओएसएक्स लायन और आर 64 संस्करण 2.13.1 में, एनए मूल्यों को मान्यता दी गई है:
> 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
लेकिन विंडोज 7 (64 बिट, एक ही आर संस्करण) पर यह सेल मानों को परिवर्तित करता है जो 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 मान क्यों नहीं हैं? मैं इसके आसपास कैसे काम कर सकता था? मेरा अनुमान है कि यह संख्याओं को संग्रहीत करने के तरीके के साथ करना है, बहुत सारे एनए मान 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()
अपने पोस्ट में शामिल कर सकते हैं ?