将网格气候数据内插到R中的不同分辨率 [英] Interpolate gridded climate data to different resolutions in R

问题描述

我有0.5度分辨率的栅格数据。我想. 使用最近邻插值法将这些值内插到0.125度和1度

如何在R中实现这一点？以下是一些示例数据：

library(raster)
ras <- raster(res=0.5)
#get polygons for Iceland
iceland <- getData('GADM', country="Iceland", level=1)
cr <- crop(ras, iceland, snap="out")
values(cr) <- 1:ncell(cr)
fr <- mask(cr, iceland)

plot(fr)
lines(iceland)

推荐答案

数据示例

library(raster)
ras <- raster(res=0.5)
#get polygons for Iceland
iceland <- getData('GADM', country="Iceland", level=1)
cr <- crop(ras, iceland, snap="out")
values(cr) <- 1:ncell(cr)
fr <- mask(cr, iceland)

1度

a <- aggregate(fr, 2, mean, na.rm=TRUE)
 

0.125度

d <- disaggregate(fr, 4)

像这样分解数据可能是可行的，但可能不合适。这显然不是更改气候数据分辨率的好方法(但您问了如何做这件事，这就是本网站的目的)

如果不能用整数乘除得到新的分辨率，则有resample

target <- raster(fr)
res(target) <- 1/13
r <- resample(fr, target, "ngb")

#or alternatively
rb <- resample(fr, target, "bilinear")

这篇关于将网格气候数据内插到R中的不同分辨率的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持IT屋！