R의 '래스터'패키지를 사용하여 Windows에서 .DEM 형식의 래스터 파일을 읽으려고합니다.
Windows 7에서 R로 데이터를로드 할 때 NA 값에 문제가 발생하지만 OSX Lion이 설치된 Mac에서는 문제가 없습니다. Windows에서 NA 값이 올바르게 읽히지 않는 것 같습니다. 문제는 왜 이런 일이 발생합니까?
사용 된 래스터 파일은 다음 R 코드와 함께 USGS에서 다운로드되었습니다.
download.file('http://edcftp.cr.usgs.gov/pub/data/gtopo30/global/e020n90.tar.gz', 'e020n90.tar.gz')
untar('e020n90.tar.gz')
그런 다음 '래스터'패키지를 사용하여 래스터를 R로 읽습니다. 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 (64Bit, 동일한 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
Windows에서 래스터를 읽을 때 NA 값이없는 이유는 무엇입니까? 어떻게 해결할 수 있습니까? 내 생각에 그것은 숫자가 저장되는 방식과 관련이 있으며 많은 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()
게시물에 게시물 을 포함시킬 수 있습니까 ?