在R中加载ESRI ArcGrid导出文件(e00)时出现问题
我正在尝试将ESRI ArcGrid导出文件加载到R中。该文件位于-。该文件的文档说明,“上述变量的网格化数据在压缩的ESRI ArcGRID导出文件中可用”,其他内容不多。链接到文档。我的主要问题是这是一个在R中加载ESRI ArcGrid导出文件(e00)时出现问题,r,geospatial,raster,R,Geospatial,Raster,我正在尝试将ESRI ArcGrid导出文件加载到R中。该文件位于-。该文件的文档说明,“上述变量的网格化数据在压缩的ESRI ArcGRID导出文件中可用”,其他内容不多。链接到文档。我的主要问题是这是一个e00文件,一种旧的Arc输出格式 正常加载光栅最直接的方法(如果这不是e00文件)是: require(raster) require(rgdal) test <- raster('/path/to/n_tw-2013.e00') 当我运行这个命令时,它会打印NULL,并创
e00
文件,一种旧的Arc输出格式
正常加载光栅最直接的方法(如果这不是e00
文件)是:
require(raster)
require(rgdal)
test <- raster('/path/to/n_tw-2013.e00')
当我运行这个命令时,它会打印NULL
,并创建两个目录,test1
,和info
,但是,这两个目录都是空的。我不确定这里出了什么问题。对于如何将其作为光栅加载到R中,以便我能够将数据提取到特定的lat/long,我将不胜感激
粘贴在此处的sessionInfo
的输出:
> sessionInfo()
R version 3.2.2 (2015-08-14)
Platform: x86_64-redhat-linux-gnu (64-bit)
Running under: CentOS release 6.7 (Final)
locale:
[1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
[3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
[5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
[7] LC_PAPER=en_US.UTF-8 LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods
[7] base
other attached packages:
[1] rgdal_1.2-5 raster_2.5-8 sp_1.2-4
loaded via a namespace (and not attached):
[1] tools_3.2.2 Rcpp_0.12.3 grid_3.2.2 lattice_0.20-33
使用
库(rgdal)
在您的示例中,应该可以:
n_tw_2013 <- readGDAL("path/to/the/file/n_tw-2013.e00")
n_tw_2013我使用rgdal 1.2-3
和光栅2.5-8
读取文件没有问题:
require(raster)
require(rgdal)
test <- raster('n_tw-2013.e00')
test
class : RasterLayer
dimensions : 775, 1440, 1116000 (nrow, ncol, ncell)
resolution : 4134.354, 4134.354 (x, y)
extent : -2950369, 3003101, 115686.8, 3319811 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=aea +lat_1=0 +lat_2=29.5 +lat_0=45.5 +lon_0=0 +x_0=0 +y_0=-96 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
data source : /Users/Martin/Desktop/n_tw-2013.e00
names : n_tw.2013
values : 0.8068619, 70.83445 (min, max)
有关读取中的.e00
文件的信息,请参阅软件包r@ImranAli见更新的问题。由于某些原因,RArcInfo
包仍然无法解决我的问题。签出我可以使用raster()
读取它。可能会更新软件包。奇怪的是,我刚刚更新了光栅
和rgdal
,但这对我来说仍然不起作用…我添加了sessionInfo。也许重新启动R,这样你就有了一个干净的环境。谢谢。您提示我在本地机器上测试此功能,在本地机器上只需调用光栅
。更新我的软件包或重新启动R都不走运,但这肯定会进一步推动我。我会尽快给你奖金。我仍然对它为什么不起作用感兴趣。您在不喜欢此代码段的计算机上使用的操作系统是什么?sessionInfo已粘贴到线程中。
require(raster)
require(rgdal)
test <- raster('n_tw-2013.e00')
test
class : RasterLayer
dimensions : 775, 1440, 1116000 (nrow, ncol, ncell)
resolution : 4134.354, 4134.354 (x, y)
extent : -2950369, 3003101, 115686.8, 3319811 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=aea +lat_1=0 +lat_2=29.5 +lat_0=45.5 +lon_0=0 +x_0=0 +y_0=-96 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
data source : /Users/Martin/Desktop/n_tw-2013.e00
names : n_tw.2013
values : 0.8068619, 70.83445 (min, max)
R version 3.3.1 (2016-06-21)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: OS X 10.12.2 (Sierra)
locale:
[1] de_DE.UTF-8/de_DE.UTF-8/de_DE.UTF-8/C/de_DE.UTF-8/de_DE.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] raster_2.5-8 sp_1.2-3
loaded via a namespace (and not attached):
[1] Rcpp_0.12.7 plotly_4.5.2 magrittr_1.5 munsell_0.4.3 colorspace_1.2-7
[6] viridisLite_0.1.3 lattice_0.20-34 R6_2.2.0 httr_1.2.1 plyr_1.8.4
[11] dplyr_0.5.0 tools_3.3.1 rgdal_1.2-3 grid_3.3.1 gtable_0.2.0
[16] DBI_0.5-1 htmltools_0.3.5 lazyeval_0.2.0 assertthat_0.1 digest_0.6.10
[21] tibble_1.2 purrr_0.2.2 ggplot2_2.2.0 tidyr_0.6.0 base64enc_0.1-3
[26] htmlwidgets_0.8 scales_0.4.1 jsonlite_1.1