本节来介绍如何使用sf包来绘制中国地图,本文由阿猹同学完成,我们提供数据可视化的有偿服务,有需要的小伙伴欢迎咨询
加载R包
rm(list = ls())
pacman::p_load(tidyverse, sf, raster, ggspatial, stars)
设置投影
crs_84 <- st_crs("EPSG:4326") ## WGS 84 大地坐标
crs_al <- st_crs("+proj=aea +lat_1=25 +lat_2=47 +lon_0=105") ## Albers Equal Area Conic投影
地图获取
http://datav.aliyun.com/portal/school/atlas/area_selector#&lat=31.769817845138945&lng=104.29901249999999&zoom=4
china_all <-
sf::st_read("https://geo.datav.aliyun.com/areas_v3/bound/100000_full.json") %>%
st_transform(crs_al)
hainan <-
sf::st_read("https://geo.datav.aliyun.com/areas_v3/bound/460000_full.json") %>%
st_transform(crs_al)
截取地图
版权声明:本文为Mrrunsen原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接和本声明。