r - 从栅格计算缓冲区内的线密度

Question

有没有办法计算空间线周围缓冲区内的道路密度(km/km²)？道路由栅格中的像素(1 像素 = 625 平方米)表示。于是我开始使用函数rasterToContour将道路像素转换为折线。 (package raster) .然后，我正在考虑计算缓冲区(以公里为单位)和缓冲区(以平方公里为单位)内的线的总长度。

r <- raster(rast_path)
x <- rasterToContour(r)

这是一个可重现的示例:

## To create raster:
library(raster)
library(rgeos)
r <- raster(ncols=90, nrows=50)
values(r) <- sample(1:10, ncell(r), replace=TRUE)

## Road raster
r[r[] < 10] <- 0
r[r[] >= 10] <- 1
plot(r)

## To create spatial lines 
line1 <- rbind(c(-125,0), c(0,60))
line2 <- rbind(c(0,60), c(40,5))
line3 <- rbind(c(40,5), c(15,-45))
line1_sp <- spLines(line1)
line2_sp <- spLines(line2)
line3_sp <- spLines(line3)

## To create buffer around lines
line2_buff <- gBuffer(line2_sp, width=20)
plot(line2_sp,add=T)
plot(line2_buff,add=T)

Answer 1

您正在寻找道路长度(公里)除以缓冲区(平方公里)，对吗？要计算道路长度，您可以使用您的方法 rasterToContour() ，尽管这在您提供的示例中无法重现。

但是要计算以平方公里为单位的缓冲区面积，您可以执行以下操作:n <- length(extract(r, line2_buff))获取 n缓冲区中的像素数，每个像素为 625 m^2。您需要的转换系数是 1 km^2 = 1,000,000 m^2。总而言之，以 km^2 为单位的缓冲区面积由下式给出:

length(extract(r, line2_buff)) * 625 / 1000000