将栅格空间分辨率从 5 角分重新采样到 1km2

Question

我有以下栅格图层：

> raster
class      : RasterLayer 
dimensions : 2160, 4320, 9331200  (nrow, ncol, ncell)
resolution : 0.08333333, 0.08333333  (x, y)
extent     : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs 

如何将分辨率从 5 角分重新采样到 1km2？

我的理解是 1km2 空间分辨率相当于 30 弧秒，因此我会使用以下解决方案： 分解（栅格，因子 = 10）。

Answer 1

library(raster)
#> Warning: package 'raster' was built under R version 4.0.5
#> Loading required package: sp
#> Warning: package 'sp' was built under R version 4.0.5
r <- raster(ncol=10, nrow=10)
values(r) <- runif(ncell(r))
par(mfrow=c(1,2))
plot(r)
rd <- disaggregate(r, fact=c(10,10))
plot(rd

它们看起来一样，但分辨率被修改了：

ncol(r)
#> [1] 10
ncol(rd)
#> [1] 100
Created on 2021-09-23 by the reprex package (v2.0.1)