匹配两个光栅的分辨率

Question

我正在处理两个光栅，每个光栅的分辨率都不同。我想知道是否有更有效的方法将较粗的光栅分辨率与较精细的光栅分辨率相匹配。现在我正在使用遮罩功能来节省一些时间，剪辑到正确的范围并更改分辨率：

library(raster)
#the raster template with the desired resolution        
r <- raster(extent(-180, 180, -64, 84), res=0.04166667) 
# set some pixels to values, others to NA
r <- setValues(r, sample(c(1:3, NA), ncell(r), replace=TRUE))

#load the raster 
lc_r1 <- raster(r)
res(lc_r1) <- 0.5
values(lc_r1) <- 1:ncell(lc_r1)
lc_r1
##class       : RasterLayer 
##dimensions  : 296, 720, 213120  (nrow, ncol, ncell)
##resolution  : 0.5, 0.5  (x, y)
##extent      : -180, 180, -64, 84  (xmin, xmax, ymin, ymax)
##coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
##data source : in memory
##names       : layer 
##values      : 1, 213120  (min, max)

#create the new finer resolution raster.
lc_r2 <- mask (lc_r1, r2)
Error in compareRaster(x, mask) : different number or columns

我也在尝试 raster 中的 disaggregate 函数，但我得到了这个奇怪的错误！

lc_r2 <- disaggregate (lc_r1, nrows=3600 )
Error: !is.null(fact) is not TRUE

这似乎暂时有效，但不确定是否正确：

lc_r2 <- disaggregate (lc_r1, fact=c(12,12 ), method='bilinear')

Answer 1

为什么这个 Error: !is.null(fact) is not TRUE 会很奇怪？如果您查看 ?disaggregate，您会发现没有参数 nrows，但是有一个您没有提供的必需参数 fact。

你可以做到

lc_r2a <- disaggregate (lc_r1, fact=12)

或

lc_r2b <- disaggregate(lc_r1, fact=12, method='bilinear')

相当于

lc_r2c <- resample(lc_r1, r)

您为什么不确定这是正确的？

但是，鉴于您想要屏蔽 lc_r1，合乎逻辑的方法是反其道而行之，更改屏蔽的分辨率，r，

ra <- aggregate(r, fact=12, na.rm=TRUE) 
lcm <- mask(lc_r1, ra)

匹配两个光栅的分辨率

Matching the resolution of two rasters

r

resolution

raster