열기 래스터에서 R! 일부 통계 작업

Question

래스터 파일 (ASCII 또는 TIFF 형식)을 열고 셀을 집계하고이 작업 이후에이 새 래스터 파일의 값과 다른 작업 사이의 상관 관계를 계산합니다. 불행하게도 내 명령에 잘못된 무엇인지 - 나는 오류 메시지가 얻을 :

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

CPL 오류 4 :`~/설교는 dods은/캐롤/TVDI 113_121이/TVDI_kamp_evi_TRANSF.asc '는 존재하지 않습니다/파일 시스템에서 지원되며 지원되는 데이터 세트 이름으로 인식되지 않습니다.

'~/Pulpit/dods/karol/TVDI 113_121/TVDI_kamp_evi_TRANSF.asc'파일 시스템에 존재하지 않으며 지원되는 데이터 세트 이름으로 인식되지 않습니다. .

Answer 1

, 당신은이 작업을 수행 할 수 있습니다

pairs(my_rasters) 

:

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) 

처럼 파일을로드 할 수 있습니다

및이 :

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]) 
    } 

래스터가 너무 커서 메모리에 두 개를 동시에 저장할 수 없으면 문제가 발생합니다. 더 좋은 질문이 없으면 더 나은 조언을하는 것이 어려울 것입니다.

Answer 2

는 래스터 (오픈) 읽으려면 한 가지 방법은 사용하는 것입니다 readGDAL :

library(rgdal) 
r <- readGDAL("~/myhome/thisdir/IhaveaFile.asc") 

이건 내 개인적인 취향이고, 유일한 이유는 그렇지 않으면 내 컴퓨터가 '아무튼 경우 GDAL.open 또는 raster가 사용하는 문제의 데이터 세트를 처리하기 위해 RAM (+ abit)이 없습니다. 당신은 파일 이름을 얻는 데 문제가있는 경우