距离栅格子集

Question

我计算距离栅格以计算从栅格像元到红色边框的距离。结果显示在下图中：

尽管如此，我现在想计算从所有单元格到边框内部（边框最右边的部分）的距离。我试图“栅格化”边界，但很难只 select 边界的右侧部分。

有什么想法吗？

提前致谢，

为了可重复性（所有可用数据）：

library(maptools) #To get the polygon data 
data(wrld_simpl)
colven <- c("Colombia", "Venezuela")
colven_map <- wrld_simpl[wrld_simpl$NAME %in% colven, ]

#Create a raster object
library(raster)
raster <- raster(colven_map, nrow=100, ncol=100)
raster[] <- 1:length(raster)
raster_colven <- mask(raster, colven_map)

#Calculate distance raster
col <- wrld_simpl[wrld_simpl$NAME == "Colombia", ]
ven <- wrld_simpl[wrld_simpl$ NAME == "Venezuela", ]
ven_l <- as(ven, "SpatialLines")
ven_p <- as(ven_l, "SpatialPoints")
distance_raster_colven <- distanceFromPoints(raster, ven_p)
distance_raster_colven <- mask(distance_raster_colven, colven_map)

Answer 1

您可以手动select您想要的片段：

plot(ven_l)
x <- crop(ven_l, drawPoly())
# draw a polygon on the map

现在

y <- rasterize(x, raster)
d <- distance(y)

距离栅格子集

Distance raster subset

r

distance

raster