ggplot2优雅的绘制高端地形图

下载第一个JSON文件即可，复制链接可通过IDM等下载工具下载

shp <- sf::read_sf("shanxi.json") # 读入陕西省地图
shanxi<- raster::mask(dem,shp) %>% # 将地图与DEM数据结合
  crop(.,extent(shp))

df_shanxi <- as.data.frame(as(shanxi,"Raster"),xy=T) #格式转换

核心点3(获取河流信息)

通过ne_download下载世界河流信息，后续要将河流信息添加到ggplot2图中，因此通过 st_as_sf将其转换为sf

rivers <- ne_download(scale = 10, type = 'rivers_lake_centerlines',category = 'physical')

设置河流经纬度信息

此处的xmin，ymin代表河流的经纬度信息，应该根据后续绘制出来的地图判读，为了代码的整洁性我将其提前放置了

rivers_cropped <- st_crop(st_as_sf(rivers), xmin = 107, xmax =111,
                          ymin = 32, ymax =39)

完成以上3大核心步骤可以说我们已经就成功了，ggplot2绘制地形图关键点就在于前3点，此前参考气象水文科研猫的文档，作者只介绍了绘图上面三大核心点真是一笔带过，在此感谢作者的简单介绍

自定义颜色

colors <- c("#33A02C","#B2DF8A","#FDBF6F","#1F78B4","#999999",
             "#E31A1C","#E6E6E6","#A6CEE3")

ggplot2数据可视化

ggplot()+
  geom_sf(data=shp,aes(fill = NULL))+
  annotation_scale(location = "bl") + # 设置距离刻度尺
  annotation_north_arrow(location="tl",
                         style = north_arrow_nautical(
                           fill = c("grey40","white"),
                           line_col = "grey20"))+  # 添加指北针
  geom_tile(data=df_shanxi,aes(x=x,y=y,fill=layer),show.legend = F)+
  geom_sf(data = rivers_cropped,col='blue',size=1)+ # 添加河流信息
  scale_fill_gradientn(colours=colors,na.value="transparent")+ 
  labs(x=NULL,y=NULL)+
  geom_sf(data=shp,fill="NA",size=1,color="black")+ # 添加地图边界
  theme_bw()+
  theme(panel.grid.major=element_blank(),
        panel.grid.minor=element_blank(),
        panel.background = element_blank(),
        legend.title = element_blank())

喜欢的小伙伴欢迎关注我的公众号

R语言数据分析指南，持续分享数据可视化的经典案例及一些生信知识，希望对大家有所帮助