R語言怎么實現(xiàn)地圖上的迷你條形圖

本篇內(nèi)容介紹了“R語言怎么實現(xiàn)地圖上的迷你條形圖”的有關(guān)知識，在實際案例的操作過程中，不少人都會遇到這樣的困境，接下來就讓小編帶領(lǐng)大家學習一下如何處理這些情況吧！希望大家仔細閱讀，能夠?qū)W有所成！

加載包：

library(maptools)

library(ggplot2)

library(plyr)

library(ggmap)

導(dǎo)入地圖素材及省會城市經(jīng)緯度數(shù)據(jù)

china_shp <-readShapePoly("c:/rstudy/bou2_4p.shp")    

china_map <- fortify(china_shp)     

province_city <- read.csv("c:/rstudy/chinaprovincecity.csv") #城市經(jīng)緯度數(shù)據(jù) 

#新建三個指標，2015年和2016年的某經(jīng)濟度量指標（虛構(gòu)）N15、N16和環(huán)比增長率Ratio。

province_city1<-mutate(province_city,N15=runif(34,min=500,max=1000),N16=runif(34,600,1100),Ratio=round((N16-N15)/N15,3))

province_data<-province_city1[sample(nrow(province_city1),10),]

合并后的數(shù)據(jù)結(jié)構(gòu)如下，隨即選取了10個城市的14、15年度某項經(jīng)濟指標（虛構(gòu)）以及同比增長率。

ggplot的現(xiàn)有圖層圖形中是沒有直接根據(jù)點坐標生成條形圖、柱形圖的，所以這里我們只能曲線救國，使用線條圖和誤差線來進行模擬。

首先來畫底圖：

ggplot()+

geom_polygon(aes(x=long, y=lat, group=group), data=china_map, fill="white", colour="grey60")

接下來使用geom_linerange函數(shù)（也就是條線圖函數(shù)）進行各個坐標點的模擬柱形繪制。

線畫出14年的指標值：

ggplot()+

geom_polygon(aes(x=long, y=lat,group=group),data=china_map, fill="white", colour="grey60")+

geom_linerange(data=province_data,aes(x=jd-0.5,ymin=wd,ymax=wd+N15/max(N15,N16)*5),size=5,color="orange",alpha=0.8)

再疊加一次16年的指標值

ggplot()+

geom_polygon(aes(x=long, y=lat,group=group),data=china_map, fill="white", colour="grey60")+

geom_linerange(data=province_data,aes(x=jd-0.5,ymin=wd,ymax=wd+0.7*N15/max(N15,N16)*5),size=3,color="#5B88A0",alpha=0.8)+

geom_linerange(data=province_data,aes(x=jd+0.5,ymin=wd,ymax=wd+0.7*N16/max(N15,N16)*5),size=3,color="#FB882C",alpha=0.8)

接下來我們來對已經(jīng)塑造好的雙柱條形圖進行美化操作，

windowsFonts(myFont = windowsFont("微軟雅黑")) 

ggplot()+

geom_polygon(aes(x=long, y=lat,group=group),data=china_map, fill="white", colour="grey60")+

geom_linerange(data=province_data,aes(x=jd-0.5,ymin=wd,ymax=wd+0.7*N15/max(N15,N16)*5),size=3,color="#5B88A0",alpha=0.8)+

geom_linerange(data=province_data,aes(x=jd+0.5,ymin=wd,ymax=wd+0.7*N16/max(N15,N16)*5),size=3,color="#FB882C",alpha=0.8)+

geom_text(aes(x=jd,y=wd-0.6,label=paste0(province_data$province,ifelse(Ratio>0,"▲","▼"),Ratio*100,"%")), data=province_data,family="myFont",size=2.5)+

theme_nothing()+

annotate("text", x=105, y=52, label="● 2015", color= "#5B88A0", size=8)+ 

annotate("text", x=105, y=49, label="● 2016", color= "#FB882C", size=8)

其實最初在網(wǎng)上看到該圖（大概是劉萬祥老師的博客，在此感謝老師的啟發(fā)），他的做法是使用geom_errbar圖層函數(shù)，也就是誤差線圖層函數(shù)，其實理念是一致的。這里給出原方法。

ggplot()+

geom_polygon(aes(x=long, y=lat,group=group),data=china_map, fill="white", colour="grey60")+

geom_errorbar(data=province_data,aes(x=jd-0.5,ymin=wd,ymax=wd+0.7*N15/max(N15,N16)*5),size=3,color="#5B88A0",width=0,alpha=0.8)+

geom_errorbar(data=province_data,aes(x=jd+0.5,ymin=wd,ymax=wd+0.7*N16/max(N15,N16)*5),size=3,color="#FB882C",width=0,alpha=0.8)+

geom_text(aes(x=jd,y=wd-0.6,label=paste0(province_data$province,ifelse(Ratio>0,"▲","▼"),Ratio*100,"%")), data=province_data,family="myFont",size=2.5)+

theme_nothing()+

annotate("text", x=105, y=52, label="● 2014", color= "#5B88A0", size=8)+ 

annotate("text", x=105, y=49, label="● 2015", color= "#FB882C", size=8)

其實如果能換個思路，使用geom_errorh函數(shù)，想必一定了以做出橫向的條形圖。

ggplot()+

geom_polygon(aes(x=long, y=lat,group=group),data=china_map, fill="white", colour="grey60")+

geom_errorbarh(data=province_data,aes(x=jd,y=wd,xmin=jd-3,xmax=jd+3*N15/max(N15,N16)),size=3,color="#5B88A0",height=0,alpha=0.8)+

geom_errorbarh(data=province_data,aes(x=jd,y=wd-0.8,xmin=jd-3,xmax=jd+3*N16/max(N15,N16)),size=3,color="#FB882C",height=0,alpha=0.8)+

geom_text(aes(x=jd+0.2,y=wd+1,label=paste0(province_data$province,ifelse(Ratio>0,"▲","▼"),Ratio*100,"%")), data=province_data,family="myFont",size=2.5)+

theme_nothing()+

annotate("text", x=105, y=52, label="● 2015", color= "#5B88A0", size=7)+ 

annotate("text", x=105, y=50, label="● 2016", color= "#FB882C", size=7)

“R語言怎么實現(xiàn)地圖上的迷你條形圖”的內(nèi)容就介紹到這里了，感謝大家的閱讀。如果想了解更多行業(yè)相關(guān)的知識可以關(guān)注億速云網(wǎng)站，小編將為大家輸出更多高質(zhì)量的實用文章！