使用核密度计算R中的道路密度？[关闭]
已关闭。这个问题需要更加集中。它当前不接受答案。 想改善这个问题吗？更新问题，使其仅通过编辑此帖子来关注一个问题。 2年前关闭。 我有一个很大的道路图形文件（约70MB），并希望将其转换为每个单元格中具有道路密度的栅格。理想情况下，如有必要，我想在R中与GDAL命令行工具一起使用。 我最初的方法是按照此线程直接计算每个单元格中线段的长度。这样可以产生所需的结果，但是即使对于比我的小得多的shapefile也非常慢。这是一个非常简单的示例，其正确的单元格值显而易见： require(sp) require(raster) require(rgeos) require(RColorBrewer) # Create some sample lines l1 <- Lines(Line(cbind(c(0,1),c(.25,0.25))), ID="a") l2 <- Lines(Line(cbind(c(0.25,0.25),c(0,1))), ID="b") sl <- SpatialLines(list(l1,l2)) # Function to calculate lengths of lines in given raster cell lengthInCell <- function(i, r, l) { r[i] <- 1 rpoly <- rasterToPolygons(r, na.rm=T) lc …