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

Posted 2023-02-25

【中文标题】如何在没有精确光栅蒙版的情况下光栅化形状？

【英文标题】：How to rasterize shape without exact raster mask?

【发布时间】：2019-08-24 06:05:31

【问题描述】：

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

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

【问题讨论】：

您能否让您的问题可重现并展示您迄今为止所做的尝试（参见here）？另外，你知道有一个GIS SE site吗？我可能会做的是create a larger raster 基于你的小，然后使用多边形掩盖它。

【参考方案1】：

示例数据

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)