我使用双线性方法尝试了helm template
和aggregate
函数。功能结果在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)