使用重采样和。在R中汇总以生成高级栅格

时间:2019-07-14 10:57:08

标签: r aggregate raster

我使用双线性方法尝试了helm templateaggregate函数。功能结果在50 km分辨率下相等,但在100 km分辨率下不同。专家可以告诉我这是什么错误吗? 谢谢

resample

汇总

setwd()
trmm25 = raster("ppt_200204.tif")
class       : RasterLayer 
dimensions  : 400, 1440, 576000  (nrow, ncol, ncell)
resolution  : 0.25, 0.25  (x, y)
extent      : -180, 180, -50, 50  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : E:/PhD research paper and Data/Data/TRMM_3B43 v7 DATA_nc/TRMM .tif monthly/ppt_200204.tif 
names       : ppt_200204 
values      : 0, 1.564059  (min, max)

重新采样

trmm50km = aggregate(trmm25, fact = 2)
class       : RasterLayer 
dimensions  : 200, 720, 144000  (nrow, ncol, ncell)
resolution  : 0.5, 0.5  (x, y)
extent      : -180, 180, -50, 50  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : in memory
names       : ppt_200204 
values      : 0, 1.482134  (min, max)
trmm100km = aggregate(trmm25, fact = 4)
class       : RasterLayer 
dimensions  : 100, 360, 36000  (nrow, ncol, ncell)
resolution  : 1, 1  (x, y)
extent      : -180, 180, -50, 50  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : in memory
names       : ppt_200204 
values      : 0, 1.369626  (min, max)

0 个答案:

没有答案