查看原文
其他

R绘图模板——世界地图的绘制及拓展!!!

王志山 科研后花园 2023-09-08

】 

点击上方

“科研后花园”

关注我们

主要代码如下:

1、设置工作环境:
rm(list = ls())setwd("D:/桌面/世界地图绘制")
#加载R包library(ggplot2)2、加载绘图世界地图的数据(这里按理说应该加载maps包,但是好像加载ggplot2包即可导入数据,没搞明白):
df_map <- map_data("world")3、绘制基础世界地图:
##绘制基础世界地图p <- ggplot(df_map,aes(long, lat, group = group)) + geom_polygon()+ borders("world",regions = ".",color = "grey50",fill="white",size=0.5)+ #主题设置 theme_void()p

4、给不同国家和海水部分添加不同颜色:
p+geom_polygon(aes(fill=region),show.legend = F)+ theme(plot.background = element_rect(fill="#76daff"))

5、连续变量——这里直接使用数据中的order列数据:
p+geom_polygon(aes(fill=order))

6、转变为极地地图:
p+coord_map("ortho")

7、添加采样点信息:
#构造采样点的经纬度df_sample1 <- data.frame( lon=runif(12, min = 0, max = 120), lat=runif(12, min = 25, max = 45), group=rep(c('1980','2000','2005','2010'),times=c(3,3,3,3)))df_sample2 <- data.frame( lon=runif(12, min = -125, max = -75), lat=runif(12, min = 30, max = 65), group=rep(c('1980','2000','2005','2010'),times=c(1,5,4,2)))df_sample3 <- data.frame( lon=runif(12, min = -5, max = 20), lat=runif(12, min = 10, max = 25), group=rep(c('1980','2000','2005','2010'),times=c(2,5,3,2)))#合并数据df_sample <- rbind(df_sample1,df_sample2,df_sample3)#顺便统计出各组的数量data <- data.frame(group=c('1980','2000','2005','2010'), value=c(6,13,10,7))##绘制子图p2 <- ggplot(data,aes(group,value,fill=group))+ geom_col()+ scale_fill_manual(values = c("#00c700", "#da0000","#0054da","#6a3d00"))+ theme_classic()+ theme(legend.position = "none", axis.text.x = element_blank(), axis.text.y = element_text(color = "black",size=12), axis.line.x = element_blank(), axis.ticks.x = element_blank(), panel.background = element_blank(), plot.background = element_blank())+ labs(x='', y='')+  scale_y_continuous(expand = c(0,0))
#添加采样点信息p+geom_point(df_sample, mapping=aes(lon, lat, fill = group), size=3,shape=21,color="black", show.legend = T)+ scale_fill_manual(name=NULL,values = c("#00c700", "#da0000","#0054da","#6a3d00"))+ annotation_custom(grob=ggplotGrob(p2), ymin = -75, ymax = -30, xmin=-190, xmax=-80)+ theme(legend.position = c(0.1,0.5))

PS: 以上内容是小编个人学习代码笔记分享,仅供参考学习,欢迎大家一起交流学习。


完整代码(附带注释)及原始数据获取:

夸克网盘(手机端下载夸克APP进行搜索提取):

链接:https://pan.quark.cn/s/053edd38709c

提取码:Rznq

更多推荐

基于MicrobiotaProcess包进行PCoA分析!

值得推荐的豆瓣高分书单,百万读者共同的选择!

跟着Nature学绘图——柱状图+散点+配对连线+显著性!!!

R绘图模板——中国地图+散点+柱状图!!!

R绘图模板——热图和柱状堆积图组合图!!!

扩增子测序数据分析还不会?小编整理的全套R语言代码助您轻松解决问题!

R绘图模板源代码获取!

不会用代码做RDA分析?没关系,Canoco软件帮你实现!

R绘图模板——热图+网络图展示mantel test相关性!!!

R绘图模板——箱线图+点线图+显著性+分组!!!

跟着Nature学绘图——绘制一张不一样的点线图!

R绘图模板——散点图+拟合曲线+边际组合图形!!!

跟着Nature学绘图——热图+显著性+间隔+注释+柱状图!

R绘图模板——散点+箱线图+小提琴图+辅助线+显著性!

跟着Nature学绘图——双向柱状图!

您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存