地图上的 R spplot 图答案

【问题标题】：R spplot diagrams on the map地图上的 R spplot 图
【发布时间】：2015-03-18 19:09:02
【问题描述】：

大家好，任何人都可以帮助我升级我的代码，并可能在我的地图中插入额外的数据。这是为我绘制迁移强度地图的代码，我正在尝试添加每个区域的 ehtnic 信息（许多小饼图）。

绘制地图

con <- url("http://biogeo.ucdavis.edu/data/gadm2/R/UKR_adm1.RData")
print(load(con))
close(con)
name<-gadm$VARNAME_1
value<-c(4,2,5,2,1,2,4,2,2,4,1,1,1,4,3,3,1,1,3,1,2,4,5,3,4,2,1)
gadm$VARNAME_1<-as.factor(value)
col<- colorRampPalette(c('cadetblue4','cadetblue1','mediumseagreen','tan2','tomato3'))(260) 
spplot(gadm, "VARNAME_1", main="Ukraine", scales = list(draw = TRUE), col.regions=col)

sp.label <- function(x, label) {
        list("sp.text", coordinates(x), label)
}

NAME.sp.label <- function(x) {
        sp.label(x, x$NAME_1)
}

draw.sp.label <- function(x) {
        do.call("list", NAME.sp.label(x))
}

spplot(gadm, 'VARNAME_1', sp.layout = draw.sp.label(gadm), col.regions=col,
       colorkey = list(labels = list( labels = c("Very low","Low", "Average",
                                                 "High","Very high"),        
                                      width = 1, cex = 1)))

这是 df 的一部分，我正在尝试将其作为饼图或条形图添加到该地图中，每个纬度 (lat) 和经度 (long) 都可以定位 mu 条形图或饼图。

df<-data.frame(region=c('Kiev oblast', 'Donezk oblast'), 
               rus=c(45,35), ukr=c(65,76), mold=c(11,44),long=c(50.43,48),
               lat=c(30.52, 37.82))

我找到了一个example 和another，但是...无法弄清楚如何在 ma 的情况下使用它。希望得到您的帮助，谢谢。

只有我现在发现的那个解决方案，但它没有升级我的代码（（（

mapPies( df,nameX="lat", nameY="long", nameZs=c('rus','ukr','mold'),
xlim=c(30,33), ylim=c(44,53), symbolSize = 2)

【问题讨论】：

标签： r plot charts geospatial

【解决方案1】：

也许这会有所帮助：

pieSP该函数根据几个属性提供SpatialPolygonsDataFrame，可用于plotGoogleMaps或spplot。

library(plotGoogleMaps)

data(meuse)
coordinates(meuse)<-~x+y
proj4string(meuse) <- CRS('+init=epsg:28992')

pies <- pieSP(meuse,zcol=c('zinc','lead','copper'), max.radius=120)
pies$pie  <- rep(c('zinc','lead','copper'),155)
pies$pie2 <- rep(1:3,155)

spplot(pies, 'pie2')

【讨论】：