以下に示すように、ポイントの強度の一種のヒートマップを生成するために、density.ppp を使用して GPS ポイントのデータセットを分析しました。
ただし、以下のように、画像をシェープファイルの境界に限定したいと思います。
最初の画像は次のように呼ばれます
x <- readShapePoly("dk.shp")
xlim<-c(min(912),max(920))
ylim<-c(min(8023),max(8030))
a<-ppp(cases@coords[,1], cases@coords[,2], xlim, ylim, unitname=c("km"))
plot(density.ppp(a, 0.1), col=COLORS)
plot(x, add=T, border="white")
ここで、cases@coords は各対象ポイントの GPS 座標であり、x は地理単位のアウトラインを提供するシェープファイルです。
2 番目のイメージは、次のコードを使用して呼び出されます。
plot(x, axes=T, col=COLORS, border="White")
これがどのように行われるか知っている人はいますか?おそらく plot() では不可能であり、別のパッケージが必要になります。
余談ですが、私が計画している次のステップは、この画像を GoogleEarth からインポートした地図の上に重ねることです。それを行う方法もまだわかりませんが、解決したら答えを投稿します
どうもありがとう