カーネル密度を使用してRの道路密度を計算しますか?[閉まっている]
閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか?この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 閉じた2年前。 道路の大きな(〜70MB)シェープファイルがあり、これを各セルの道路密度を持つラスターに変換したい。理想的には、必要に応じてGDALコマンドラインツールとともにRでこれを実行したいと思います。 私の最初のアプローチは、このスレッドごとに各セルの線分セグメントの長さを直接計算することでした。これにより目的の結果が得られますが、シェープファイルが私のものよりはるかに小さい場合でも非常に遅くなります。正しいセル値が明らかな非常に単純化された例は次のとおりです。 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 …