地理的分析を実行しました
library(PBSmapping); library(RgoogleMaps); library(maptools)
cases<-read.table("cases.csv")
cases.ppp<-ppp(cases@coords[,1], cases@coords[,2], unitname=c("km"))
dens.cas <- density.ppp(cases.ppp, bw)
plot(dens.cas)
次の画像が得られます。
これをいくつかの Google マップ タイルにオーバーレイしたいのですが、次のように R に適切なタイルをフェッチするように指示できました。
shp <- importShapefile('map.shp', projection="LL");
bb <- qbbox(lat = shp[,"Y"], lon = shp[,"X"]);
MyMap <- GetMap.bbox(bb$lonR, bb$latR, destfile = "map.jpg");
ただし、分析結果をマップ タイル上に配置するための解決策が見つからないようです。このソリューションを使用してシェープファイルをオーバーレイできることはわかっています。
https://stackoverflow.com/a/2105912/1540735
これは機能し、(以下のように) 地理位置情報に関しては完全に調整されていますが、分析の結果が失われます。
PlotPolysOnStaticMap(MyMap, shp, lwd=.5, col = rgb(1,1,1,0.2), add = F)
これらをオーバーレイする方法を知っている人はいますか？色を生成する数値は dens.cas$v に格納されており、私が知る限り、128*128 の行列として格納されています。
助けてくれる人に感謝します