r - 在 R 中打开栅格!和一些统计操作

Question

我想打开栅格文件(ASCII 或 TIFF 格式)，聚合它们的单元格，并在此操作后计算这个新栅格文件中的值与另一个栅格文件中的值之间的相关性。不幸的是，我不知道我的命令有什么问题 - 我收到一条错误消息:

x <- GDAL.open('~/Pulpit/dods/karol/TVDI 113_121/TVDI_kamp_evi_TRANSF.asc') 

CPL ERROR 4: `~/Pulpit/dods/karol/TVDI 113_121/TVDI_kamp_evi_TRANSF.asc' does not exist in the file system, and is not recognised as a supported dataset name.

Error in .local(.Object, ...) : `~/Pulpit/dods/karol/TVDI 113_121/TVDI_kamp_evi_TRANSF.asc' does not exist in the file system, and is not recognised as a supported dataset name.

Answer 1

如果您在获取文件名时遇到问题，您可以这样做:

my_asc_files = dir("../somepath", pattern="*.asc", recursive=T, full.names=T)
files_I_want = my_asc_files[c(1,12,32,33)]

然后你可以像这样加载你的文件

library(raster)
my_rasters = lapply(files_I_want, raster)

那么你可以这样做:

pairs(my_rasters) 

还有这个:

for(i in 1:length(my_rasters)) 
  for(j in i:length(my_rasters))   
    if(i != j) {
      df = na.omit(data.frame(values(my_rasters[[i]]), values(my_rasters[[j]])))
      cor(df[,1], df[,2])
    }

虽然如果栅格太大以至于您不能同时在内存中保存两个，您将会遇到问题。如果没有更好的问题，就很难给你更好的建议。