根据道路shapefile计算道路密度栅格

Question

我希望将美国东半部的道路(包括每条道路的一列长度)的shapefile转换为1x1公里道路密度的栅格，使用R。

我在Arcmap中找不到直接的方法(Line density与单元格中心的半径一起工作，而不仅仅是单元格)。

Answer 1

这是一个从栅格单元创建多边形的解决方案(改编自我的答案here)。您可能需要对数据集的子集执行此操作，然后进行合并。

示例数据

library(terra)
v <- vect(system.file("ex/lux.shp", package="terra"))
roads <- as.lines(v)
rs <- rast(v)

解决方案

values(rs) <- 1:ncell(rs)
names(rs) <- "rast"    
rsp <- as.polygons(rs)

rp <- intersect(roads, rsp)

rp$length <- perim(rp) / 1000 #km
x <- tapply(rp$length, rp$rast, sum)

r <- rast(rs)
r[as.integer(names(x))] <- as.vector(x)

plot(r)
lines(roads)