2015-01-17 2 views
1

저는 R 사용자이고 18 Geotiff 파일을 R로 읽고 관심 영역에 따라 tiff를 마스크했습니다. 이제 18 개의 마스크 된 지오 티프 파일을 쓸 수 없습니다.R에 여러 GeoTiff 파일 쓰기

require(raster) 
raster_data <- list.files(path=getwd()) 
s <- stack(raster_data) 
spf<-readShapePoly("basin.shp") 
rc<-crop(s, extent(spf)) 
rm<-mask(rc, spf) 
rf <- writeRaster(rm, filename=outputFile, overwrite=TRUE) 

나는 18 Geotiff 파일 내 Shape 파일을 형성 마스크 원하지만 출력은 하나의 TIF 파일이며, 그것은 열리지 않습니다. 인터넷에서 검색 옵션을 사용하지만 내 쿼리에 대한 적절한 답변을 찾을 수 없습니다.

고맙습니다.

+0

같은

뭔가 당신은 당신이 당신의 접근 방식에있어 오류를 설명 할 수? –

+0

'rf <- writeRaster (u, filename = outputFile, overwrite = TRUE)' –

+0

또는'writeRaster (rm, filename = outputFile) '앞에'u <- unstack (rm) , overwrite = TRUE, bylayer = T)' –

답변

1

당신은 확실히 귀하의 질문에 대한 답변을 찾을 수 있습니다. herehere 전에 문의했으며 대답은 here뿐 아니라 here이고 bylayer 옵션은 writeRaster입니다.

bylayer = T을 목록에 스택 해제 할 필요가없는 경우 시도하십시오.

library(raster) 
r1 <- raster(ncol=10, nrow=10) 
r1[] <- 1:100 
s <- stack(r1, r1) 
writeRaster(s, '~:/r.tif', bylayer = T) 
+0

이전에 요청되었다는 것을 알고 있다면 게시물을 중복으로 플래그해야합니다. – jbaums

+0

오류 메시지는 .local (x, filename, ...) 오류입니다. 파일명의 수가 1보다 크지 만 레이어 수와 같지 않습니다. – user3356873

+0

@ user3356873 - 'writeRaster (s, names (들), bylayer = T, 형식 = 'GTiff')'. – jbaums