figure1.jpg

基于shp格式的底图，R包sf进行矢量空间数据的转换，R包tmap进行可视化

加载工具包

library(sf)
library(tmap)
library(ggplot2)
library(tmaptools)
sf_use_s2(FALSE)
library(scales)

设置颜色变量

# forest, grass, shrub
arid = c('#8E0152', '#C51B7D','#DE77AE','#F1B6DA')
tropical = c('#40004B','#762A83','#9970AB','#C2A5CF')
temperate = c("#003C30",'#01665E','#35978F','#80CDC1')
continental = c("#543005",'#8C510A','#BF812D','#DFC27D')
polar = '#FFFF33'
show_col(c(arid, tropical, temperate, continental, polar))
show_col(grep('^grey', colors(), value = T))

设置R语言中的地理/投影坐标系统

proj = "+proj=robin +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"

导入样点数据

包括坐标信息和分组信息

sites <- read.table("./sites.csv", sep = ",", header = T)
names(sites)[4] = 'Biomes'
sites.sf <- st_as_sf(sites, coords = c("lon", "lat"), crs=4326)#25831
sites.sf.robin <- st_transform(sites.sf, proj)

导入底图

# land = read_sf('10m_physical/ne_10m_land.shp') %>%
#   # subset(CONTINENT != 'Antarctica') %>%
#   st_geometry() %>%
#   st_transform(proj)

country = read_sf('10m_cultural/ne_10m_admin_0_countries.shp')%>%
  # subset(name != 'Antarctica') %>%
  st_geometry() %>%
  st_transform(proj)

绘图

fig = tm_shape(country)+
  tm_graticules(col = 'grey80',labels.show = F) +
  tm_borders(col = 'white', lwd = .01) +
  tm_fill(col = 'grey',border.col = 'black') +
  tm_layout(bg.color="white",
            outer.margins = c(0.001,0.001,0.001,0.1),
            inner.margins=c(.02, .02, .04, .02), #bottom, 
            earth.boundary = T,
            earth.boundary.color = 'grey80',
            space.color="grey90",
            legend.show = T,
            outer.bg.color = 'white')