R语言绘制地图

library(geoviz)
library(tidyverse)
library(sf)
library(terra)
library(rasterVis)
library(ggspatial)
library(rgdal)
library(rnaturalearth)
library(rnaturalearthdata)
library(raster)

download.file("https://geo.datav.aliyun.com/areas_v3/bound/530000_full.json",destfile="yunnan.json")
yunnan_shp=sf::read_sf("yunnan.json") ## 
ggplot(yunnan_shp) +
  geom_sf(aes(fill=name)) + 
  theme_bw() -> p

export::graph2ppt(p, file = "云南省地图.pptx")
# ggsave("yunnan.pdf",height = 9,width = 9)

sf::read_sf("./data/sanqimetagenome/data/530000_云南省_full.json") %>% 
  dplyr::mutate(name = "") %>% 
  ggplot() +
  geom_sf(aes(geometry = geometry),
          show.legend = FALSE) +
  geom_sf_text(aes(label = name), size = 2) +
  geom_point(data = df.sample.num, aes(x = longitude, y = latitude, size = n), color = "#1380A1") +
  labs(x = "Longitude", y = "Latitude") +
  pac4xiang::mytheme() +
  theme(legend.position = c(0.1,0.15),
        legend.title = element_blank())-> p.map

ggsave(p.map,
       file = "./data/sanqimetagenome/results/Figure 1/samples.map.png",
       width = 6, height = 6, dpi = 500)