[英]geom_polygon in ggplot2 - Facets for different fillings?
有沒有辦法使用構面組合這三個地圖?
emi <- readShapePoly('prov2011_g.shp')
names(emi) <- c('REG', 'ID', 'NOME', 'SHAPE', 'AREA', 'MIG1', 'MIG2', 'MIG3')
emi_geom <- poly_coords(emi)
map <- qplot(PolyCoordsY, PolyCoordsX, data = emi_geom, group = Poly_Name, geom = 'polygon', fill = MIG1)
map1 <- qplot(PolyCoordsY, PolyCoordsX, data = emi_geom, group = Poly_Name, geom = 'polygon', fill = MIG2)
map2 <- qplot(PolyCoordsY, PolyCoordsX, data = emi_geom, group = Poly_Name, geom = 'polygon', fill = MIG3)
我遵循了這個很棒的教程( http://www.r-bloggers.com/maps-with-ggplot2/ )。 poly_coords函數從shapefile中提取坐標。 我曾嘗試使用melt(reshape2),但我不知道如何將其應用於這種情況-使用其他geom_ *會容易得多。 希望你能幫助我。
發布答案,這樣將來的SO搜索者都不會認為qplot()
甚至是個好主意。 無需接受或+1。 請停止從博客文章中剪切和粘貼代碼,甚至不嘗試使用它。
library(rgdal)
library(maptools)
library(ggplot2)
library(viridis)
prov <- readOGR("prov2011_g.shp", "prov2011_g", stringsAsFactors=FALSE, verbose=FALSE)
longlat <- "+init=epsg:4121 +proj=longlat +ellps=GRS80 +datum=GGRS87 +no_defs +towgs84=-199.87,74.79,246.62"
prov <- SpatialPolygonsDataFrame(spTransform(prov, CRS(longlat)), prov@data)
prov_map <- fortify(prov, region="NOME_PRO")
fac <- data.frame(area=rep(prov@data$NOME_PRO, 3),
measure=rep(c("A", "B", "C"), each=nrow(prov@data)),
val=sample(100, nrow(prov@data)*3, replace=TRUE))
gg <- ggplot()
gg <- gg + geom_map(data=prov_map, map=prov_map,
aes(x=long, y=lat, map_id=id),
color="#2b2b2b", size=0.1, fill=NA)
gg <- gg + geom_map(data=fac, map=prov_map,
aes(fill=val, group=measure, map_id=area))
gg <- gg + scale_fill_viridis()
gg <- gg + facet_wrap(~measure)
gg <- gg + coord_map()
gg <- gg + ggthemes::theme_map()
gg
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.