时间戳

首页 > 解决方案 > 如何在没有精确光栅蒙版的情况下光栅化形状？

r - 如何在没有精确光栅蒙版的情况下光栅化形状？

问题描述

我有一个我想要栅格化的多边形，我有一个带有投影和像素大小等的栅格，我想用于多边形，但栅格小于多边形。

我需要一个具有多边形范围的新栅格，还需要像素大小和像素顺序以及像较小栅格一样的位置。

标签： r

解决方案

示例数据

library(raster)
p <- readRDS(system.file('external/lux.rds', package='raster'))
r <- raster(p[4,])
res(r) <- 0.05

假设您要栅格化所有 p。r具有所需的分辨率，但范围太小。

要获得更大的栅格，您可以执行

 rr1 <- setExtent(r, extent(p)+res(r), keepres=TRUE, snap=TRUE)
 #or
 rr2 <- extend(r, extent(p)+res(r))

其次是

 x <- rasterize(p, rr1)

你也可以对此做一些变化

 r2 <- raster(xmn=5.7, xmx=6.6, ymn=49.4, ymx=50.2)
 res(r2) <- 0.03

 y <- rasterize(p, r2)

推荐阅读