将shapefile内的地理分析结果绘制到Google地图图块上

Question

我已经进行了地理分析

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)

产生以下图像：

我想在一些谷歌地图图块上叠加这个，我已经能够告诉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");  

但是，我似乎无法找到任何解决方案来将分析结果放在地图图块上。我知道可以使用此解决方案覆盖shapefile：

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矩阵。

非常感谢任何可以提供帮助的人