这一期呢，小编给回应读者的需求，绘制在地图上添加饼图的可视化作品，这类图表在GIS或生态方面涉及的比较多哈。本期内容主要如下：
-
R-scatterpie包主要函数介绍
-
R-scatterpie包实例演示
R-scatterpie包主要函数介绍
要想在地图上绘制饼图，使用常规的ggplot2较为繁琐，得益于丰富第三方包，我们可以使用scatterpie包轻松实现该需求。scatterpie包主要有：
-
geom_scatterpie()
-
geom_scatterpie_legend()
前者主要用于绘制饼图绘制，后者则对应绘制饼图大小的图例，具体可参考如下链接：https://cran.r-project.org/web/packages/scatterpie/vignettes/scatterpie.html。
接下来，小编将结合具体实际例子用于帮助小伙伴们理解scatterpie绘制函数，快速绘制带饼图的空间地图可视化作品。
R-scatterpie包实例演示
这一部分我们分别通过虚构数据绘制美国和中国的饼图地图作品，所使用的数据也会在文末给出免费获取方式。
-
以美国地图为例 数据处理代码如下(较为简单且都有注释)：
#导入包
library(readxl)
library(albersusa)
library(biscale)
library(sf)
library(tidyverse)
library(hrbrthemes)
library(ggtext)
library(scatterpie)
library("ggsci")#读入虚构数据
data_test <- read_excel("usa_lon_lat.xlsx")
#转换投影
data_test_pro <- data_test %>% sf::st_as_sf(coords = c("lon", "lat"),crs = 4326)%>% st_transform(crs = 2163)
#将结果经纬度信息单独列出并准换成DF类型
data_test_pro <- data_test_pro %>% dplyr::mutate(lon = sf::st_coordinates(.)[,1],lat = sf::st_coordinates(.)[,2])%>% select(values01,values02,values03,total,lon,lat)%>%data.frame(stringsAsFactors = F) #可视化绘制
US_map <- ggplot()+geom_sf(data = us_sf,size=.125,aes(geometry = geometry)) +scatterpie::geom_scatterpie(data=data_select,aes(x = lon,y = lat,r=total*40),cols=c("values01","values02","values03"),alpha=.9) +geom_scatterpie_legend(data_select$total*40, x=2068885.4, y=-2242789,n = 3,labeller=function(x) x/40) +scale_fill_jco(name="Type") +labs(x="",y="",title = "Example of <span style='color:#D20F26'>scatterpie::geom_scatterpie function</span>",subtitle = "processed charts with <span style='color:#1A73E8'>usa map()</span>",caption = "Visualization by <span style='color:#DD6449'>DataCharm</span>") +theme_ipsum(base_family = "Roboto Condensed") +theme(plot.title = element_markdown(hjust = 0.5,vjust = .5,color = "black",size = 20, margin = margin(t = 1, b = 12)),plot.subtitle = element_markdown(hjust = 0,vjust = .5,size=15),plot.caption = element_markdown(face = 'bold',size = 12))
注意：geom_scatterpie_legend()中x,y的位置为转换后的数据哈。可视化结果如下：
Scatterpie Map Example of USA
-
以中国地图为例 同样的，我们还是虚构数据构建关于中国的类似地图的绘制，详细绘制代码如下：
china_data <- read_excel("China_lon_lat.xlsx")
china_data_pro <- china_data %>% sf::st_as_sf(coords = c("lon", "lat"),crs = 4326)%>% st_transform(crs = 2343)%>% dplyr::mutate(lon = sf::st_coordinates(.)[,1],lat = sf::st_coordinates(.)[,2])%>% select(data01,data02,data03,total,lon,lat) %>% data.frame(stringsAsFactors = F)
#读取中国地图文件信息
china_shp <- "中国省级地图GS（2019）1719号.geojson"
nine <- "九段线GS（2019）1719号.geojson"
china <- sf::read_sf(china_shp)
nine_line <- sf::read_sf(nine)
# 可视化绘制
china_map <- ggplot()+geom_sf(data = china,fill="gray90",size=.125,color="black") + geom_sf(data = nine_line,size=.125) + coord_sf(crs = 2343) +scatterpie::geom_scatterpie(data=china_data_select,aes(x = lon,y = lat,r=total*500),cols=c("data01","data02","data03"),alpha=.9) +geom_scatterpie_legend(china_data_select$total*500,x = 25000,y=1500000,n = 3,labeller=function(x) x/500) +scale_fill_nejm(name="Type") +labs(x="",y="",title = "Example of <span style='color:#D20F26'>scatterpie::geom_scatterpie function</span>",subtitle = "processed charts with <span style='color:#1A73E8'>China map</span>",caption = "Visualization by <span style='color:#DD6449'>DataCharm</span>") +theme_ipsum(base_family = "Roboto Condensed") +theme(plot.title = element_markdown(hjust = 0.5,vjust = .5,color = "black",size = 20, margin = margin(t = 1, b = 12)),plot.subtitle = element_markdown(hjust = 0,vjust = 0,size=15),plot.caption = element_markdown(face = 'bold',size = 12))
可视化结果如下(小编使用了ggsci包的不同颜色主题哦)：
Scatterpie Map Example of China
以上，就是本期推文的全部内容，较为简单且常用，希望小伙伴们能够喜欢，别忘了获取练习数据进行复现哦~
总结
本期算是解决一个读者的提问了哈，希望大家有所收获