2016-09-06 5 views
0

나는 세부있는 전 세계의 래스터 파일이 있습니다하나의 래스터 형식을 다른 형식으로 변환하려면 어떻게해야합니까?

class  : RasterLayer 
dimensions : 18000, 43200, 777600000 (nrow, ncol, ncell) 
resolution : 0.008333333, 0.008333333 (x, y) 
extent  : -180, 180, -60, 90 (xmin, xmax, ymin, ymax) 
coord. ref. : +proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs 
data source : E:\bc26pr50\bc26pr501.tif 
names  : bc26pr501 
values  : 0, 1100 (min, max) 

을 그리고 미국만을위한 또 다른 래스터 파일 모방을 부분 집합하려면 :

class  : RasterLayer 
dimensions : 3500, 7000, 24500000 (nrow, ncol, ncell) 
resolution : 0.01, 0.01 (x, y) 
extent  : -130, -60, 20, 55 (xmin, xmax, ymin, ymax) 
coord. ref. : NA 
data source : E:\dem.tif 
names  : dem 
values  : -85.25208, 4385.539 (min, max) 

가 어떻게 변환 할 수 있습니다 해상도뿐만 아니라 데이터 집합의 하위 집합? RGDAL을 사용해야합니까? 모든 포인터가 도움이 될 것입니다.

답변

0

첫째, 래스터 패키지에서 작물() 함수를 사용

r1.crop <- crop(r1,r2) # probably not obligatory, but will speed up the next action

그런 다음 당신은 당신의 목표에 매우 가까이 있어야 리샘플링 기능

r1.resamp <- resample(r1.crop, r2) # can be slow

를 사용합니다. 보다 효율적인 방법이 있지만 먼저이 방법을 시도한 다음 최적화하십시오.

1

기능을 패키지 gdalUtils에 사용하는 것이 좋습니다. "te"및 "tr"매개 변수를 설정해야하며 '탭'으로 비트를 재생해야하지만 원하는 것을 쉽게 얻을 수 있습니다.

HTH,

로렌조