一篇舊文，解決一個困擾已經的小技術問題，權當是學習ggplot2以來的整理回顧與查漏補缺。windows
今天這一篇是昨天推送的基礎上進行了進一步的深化，主要講如何在離散顏色填充的地圖上進行氣泡圖圖層疊加。
微信
爲了使得案例先後一致，仍然使用昨天的數據集。ide
加載包：學習
library("ggplot2")ui
library("plyr")spa
library("maptools").net
library("sp")orm
library("ggthemes")對象
導入中國省界地圖：blog
china_map<-readShapePoly("D:/R/rstudy/CHN_adm/bou2_4p.shp")
data1<- china_map@data
data2<- data.frame(id=row.names(data1),data1)
數據格式轉化及業務數據合併：
china_map1 <- fortify(china_map)
china_map_data <- join(china_map1,data2, type = "full")
mydata <- read.csv("D:/R/rstudy/Province/geshengzhibiao.csv")
china_data <- join(china_map_data, mydata, type="full")
各省省會城市經緯度數據：
province_city <- read.csv("D:/R/rstudy/Province/chinaprovincecity.csv")
###根據本身的數據量級和具體業務須要設置分割點
mydata<-mydata[,-c(5,6)]
mydata$zhibiao<-rnorm(33,100,50)
mydata$zhibiao<-abs(mydata$zhibiao)
mydata$zhibiao2<-round(mydata$zhibiao,0)
mydata$fau <- cut(mydata$zhibiao, breaks = c(0,50,100,150,200,250))
###將轉換的分段因子變量從新命名爲咱們須要的分段閥值：
mydata$fam<-factor(mydata$fau,levels=c('(0,50]','(50,100]','(100,150]','(150,200]','(200,250]'),labels=c('0~50','50~100','100~150','150~200','200~250'),order=TRUE)
#將業務數據與地理信息數據合併：
china_data <- join(china_map_data, mydata, type="full")
windowsFonts(myFont = windowsFont("微軟雅黑"))
離散標度填充：
ggplot(china_data, aes(x = long, y = lat, group = group,fill =fam)) +
geom_polygon(colour="white")+
scale_fill_brewer(palette="Blues") + ###Blues&Greens
coord_map("polyconic") +
ggtitle("某公司2015~2016年度營業情況分佈圖")+
guides(fill=guide_legend(reverse=TRUE,title=NULL))+
theme_map() %+replace% theme(legend.position = c(0.08,0.4),legend.text.align=1)
在離散顏色標度的基礎上添加各省份散點圖：
ggplot() +
geom_polygon(data=china_data, aes(x=long,y=lat,group=group,fill=fam),colour="white")+
geom_point(data=province_city,aes(x=jd,y=wd),colour="red")+
scale_fill_brewer(palette="Blues") + ###Blues&Greens
coord_map("polyconic") +
ggtitle("某公司2015~2016年度營業情況分佈圖")+ #寫入標題
guides(fill=guide_legend(reverse=TRUE,title=NULL))+
theme_map() %+replace% theme(legend.position = c(0.08,0.4),legend.text.align=1)
更該散點圖爲氣泡圖：
province_city$PerforamA<-round(rnorm(34,100,30),0)
province_city$PerforamB<-round(rnorm(34,100,30),0)
ggplot() +
geom_polygon(data=china_data, aes(x=long,y=lat,group=group,fill=fam),colour="white")+
geom_point(data=province_city,aes(x=jd,y=wd,size=PerforamA),shape=21,fill="#8E0F2E",colour="black",alpha=0.6)+
scale_fill_brewer(palette="Blues") + ###Blues&Greens
scale_size_area(max_size=6)+
coord_map("polyconic") +
ggtitle("某公司2015~2016年度營業情況分佈圖")+ #寫入標題
guides(fill=guide_legend(reverse=TRUE,title=NULL),size=guide_legend(reverse=TRUE,title=NULL))+
theme_map() %+replace% theme(legend.position = c(0.08,0.4),legend.text.align=1)
原本打算再繼續在氣泡圖的基礎上進行顏色漸變填充呢，能夠試了一下，這樣的話前面的底圖離散顏色標度填充的的時候已經使用過了一個fill屬性設置選項，而要對氣泡圖進行顏色漸變填充就要再使用一次fill屬性，但是目前爲止我還不知道如何在多圖層中出現多個fill屬性的時候如何識別並分別進行標度設置，當我運行以下代碼時，軟件沒法識別兩個顏色標度設置分別對應的指標，於是圖表沒法跑出來：
ggplot() +
geom_polygon(data=china_data, aes(x=long,y=lat,group=group,fill=fam),colour="white")+
scale_fill_brewer(palette="Blues") + ###Blues&Greens geom_point(data=province_city,aes(x=jd,y=wd,size=PerforamA,fill=PerforamB),shape=21,fill="#8E0F2E",colour="black",alpha=0.6)+
scale_size_area(max_size=6)+
scale_fill_gradient2(low="DarkCyan", mid="Azure", high="Sienna", midpoint=median(province_city$PerforamB))+
coord_map("polyconic") +
ggtitle("某公司2015~2016年度營業情況分佈圖")+ #寫入標題
guides(fill=guide_legend(reverse=TRUE,title=NULL),size=guide_legend(reverse=TRUE,title=NULL))+
theme_map() %+replace% theme(legend.position = c(0.08,0.4),legend.text.align=1)
Scale for 'fill' is already present. Adding another scale for 'fill', which will replace the existing scale.
Error: Discrete value supplied to continuous scale
實在是太遺憾了，不知道哪位大神知道如何處理多圖層相同的fill屬性，能夠告知在下，這裏先行謝過了。
那麼這個問題真的那麼難解決嗎，假若放在半年之前，確實如此，那個時候爲此糾結了好久，一直耿耿於懷，但是現在的我，水平已經今非昔比了哈哈哈（容我傲嬌一下~）
其實問題是這樣的，該圖表對象的須要三層顏色標度映射。
最底層的離散填充標度；
氣泡圖的大小標度：
氣泡圖的填充標度；
困擾個人問題是，底層的多邊形填充使用了一次fill屬性，而氣泡的顏色填充又使用了一次fill屬性，因此兩個顏色標度——
撞車了……
軟件沒法識別兩個標度參數：
scale_fill_brewer(palette="Blues")
scale_fill_gradient2(low="DarkCyan", mid="Azure", high="Sienna")
分別是做用於那個圖層的：
geom_polygon()
geom_point()
因此本身就崩潰罷工了~
當時想的方法是，是否能夠給兩個fill標度進行組合或者命名，最後也沒搞定該如何進行區分。
前不久跟個人一個朋友聊天，聊起這個問題，也沒有直接的解決方法， 能夠聊天過程當中忽然發現了一條線索，R語言環境中的形狀一共有25種，其中1~20種僅有colour屬性而沒有fill屬性，21~25種既有colour屬性又有fill屬性，而我平時做圖使用氣泡圖的時候，用慣了21號形狀（主要是擁有兩個顏色調整參數，比較利於後期屬性調整）。
因此針對這個案例，使用21號形狀的fill屬性進行氣泡內顏色標度映射，必然與多邊形的fill映射衝撞，那麼如今問題就好辦了，放棄使用21號形狀，選擇一個只有colour屬性映射的形狀編號進行映射。
其中1六、19號都符合要求，我隨表挑選一個16號吧……
ggplot() +
geom_polygon(data=china_data, aes(x=long,y=lat,group=group,fill=fam),colour="white")+
scale_fill_brewer(palette="Greens") + ###Blues&Greens
geom_point(data=province_city,aes(x=jd,y=wd,size=PerforamA,colour=PerforamB),shape=16)+
scale_size_area(max_size=6)+
scale_colour_gradient(low="white",high="#D73434")+
coord_map("polyconic") +
ggtitle("某公司2015~2016年度營業情況分佈圖")+ #寫入標題
guides(fill=guide_legend(reverse=TRUE,title=NULL),size=guide_legend(reverse=TRUE,title=NULL))+
theme_map() %+replace% theme(legend.position = c(0.08,0.4),legend.text.align=1)
OK，完美的解決了標度重疊問題，如今該地圖已經用了三個可用的顏色標度了！
歡迎關注魔方學院QQ羣
本文分享自微信公衆號 - 數據小魔方（datamofang）。
若有侵權，請聯繫 support@oschina.cn 刪除。
本文參與「OSC源創計劃」，歡迎正在閱讀的你也加入，一塊兒分享。