2

我进行了地理分析

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

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

4

0 回答 0