geom_scatterpie绘制带饼图的地图异常问题求助

环境信息 #

• R版本：4.1.1
• scatterpie版本：0.2.1
• ggplot2版本：3.4.2

问题说明 #

尝试绘制以太平洋为中心、叠加饼图的世界地图，但生成的map_try_1和map_try_2图像显示异常，现附上运行脚本及输入文件内容，请求排查问题原因。

运行脚本 #

library(ggplot2)
library(scatterpie)

##############################################################################

lat_lon_province <- read.csv("lat_lon_province.txt", sep ="	")
pies <- read.csv("pies.txt", sep ="	")
colnames(pies)[colnames(pies) == "Province"] ="PROVINCE"

pies_ecoregion_lat_lon <- merge(pies, lat_lon_province, by = "PROVINCE")
pies_ecoregion_lat_lon[is.na(pies_ecoregion_lat_lon)] <- 0
colnms=c("Copepoda", "Isopoda", "Nematoda", "Cestoda", "Trematoda", "Monogenea", "Myxozoa")
pies_ecoregion_lat_lon$Total_parasites<-rowSums(pies_ecoregion_lat_lon[,colnms])

n <- nrow(pies_ecoregion_lat_lon)
pies_ecoregion_lat_lon$radius <- 6 * abs(rnorm(n))

lon_min <- 80
lon_max <- 270
lat_min <- -50
lat_max <- 50

# 使用map_data()
pacific_world <- map_data("world2")

# world2MapEnv或world2
# 这是世界数据库的替代版本，基于纬度范围[0, 360)，将太平洋置于地图中心。
mapplot_pacific <- ggplot(pacific_world) + 
  geom_map(data = pacific_world, map = pacific_world, 
           aes(x=long, y=lat, map_id=region), col = "white", fill = "gray50") + theme_bw()
mapplot_pacific + 
  coord_map(xlim=c(lon_min, lon_max), ylim=c(lat_min, lat_max)) +
  geom_scatterpie(aes(x=Longitude, y=Latitude, r=radius), 
                  data = pies_ecoregion_lat_lon, cols = colnames(pies_ecoregion_lat_lon[,c(3:9)]), 
                  color=NA, alpha=.8) +
  geom_scatterpie_legend(pies_ecoregion_lat_lon$radius, x=240, y=35)
ggsave("map_try_1.png", width = 10)

###############################################################################

pacific_map <- fortify(pacific_world, 
                       xlim = c(lon_min, lon_max), 
                       ylim = c(lat_min, lat_max), 
                       fill = T, plot = F)

pacific <- ggplot(pacific_map)

ggplot(pacific_map, aes(x = long, y = lat, group = group)) + geom_polygon() + 
  coord_map(xlim=c(lon_min, lon_max), ylim=c(lat_min, lat_max)) +
  geom_scatterpie(aes(x=Longitude, y=Latitude, group = PROVINCE, r=radius), 
                  data = pies_ecoregion_lat_lon, cols = colnames(pies_ecoregion_lat_lon[,c(3:9)]), 
                  color=NA, alpha=.8) +
  geom_scatterpie_legend(pies_ecoregion_lat_lon$radius, x=240, y=35) + 
  theme_bw()
ggsave("map_try_2.png", width = 10)

输入文件内容 #

pies.txt #

lat_lon_province.txt #

内容的提问来源于stack exchange，提问作者Christina P