R 将ShapeFile中的地理分析结果打印到Google地图分幅上
我进行了地理分析R 将ShapeFile中的地理分析结果打印到Google地图分幅上,r,google-maps,R,Google Maps,我进行了地理分析 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) 有人知道如何覆盖这些吗?据我所
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)
有人知道如何覆盖这些吗?据我所知,产生颜色的数值存储在dens.cas$v中,并存储为128*128矩阵
非常感谢任何能帮忙的人
shp <- importShapefile('map.shp', projection="LL");
bb <- qbbox(lat = shp[,"Y"], lon = shp[,"X"]);
MyMap <- GetMap.bbox(bb$lonR, bb$latR, destfile = "map.jpg");
PlotPolysOnStaticMap(MyMap, shp, lwd=.5, col = rgb(1,1,1,0.2), add = F)