怎么画中国各省地图？
先拿出一张纸。
再画中国轮廓，比如说这张。
请点击输入图片描述
再画出省份，如果还要画省会的话，建议画大些。
请点击输入图片描述
一张完好的中国地图就出来了，如下图。
请点击输入图片描述
END
甘肃不要画于过宽，会影响宁夏、陕西、山西、河北、北京、天津。
不要画于过低，影响海南面积。
注意事项
简介
关于地理地图绘制已经由师弟师妹系统出了几期。具体推文可见：
- 注意：这里的 dem_data 可以是连续型或离散型变量（各市 GDP 排名），以下例子使用 runif() 随机生成连续型数据。library(leaflet) #加载包
- library(leafletCN)
- region = regionNames("浙江")
- dem_data = runif(length(region)) # 这里使用自己的数据
- dat = data.frame(region,dem_data)
- map = leafletGeo("浙江", dat)
本推文主要介绍，如何使用 leaflet 包绘制省/市级地图并以浙江省、温州市作为案例。读者可以使用函数保存该图片，并将其使用到论文中。
但是本推文的省市级绘图还存在几个问题，小编还没有能力通过代码解决：1. 如何凸显各市/县的名字到填充颜色上方？2. 有些市、县发生变化（从县变为区），但 leaflet 包没有及时更新，应该如何处理？
以上是小编在实际科研中存在的问题，我的“笨”办法是： 。或者使用 ppt，手动修改。如果读者们有什么好办法，可以留言分享下，十分感谢。
如果觉得推文有用，求关注，收藏，分享，点赞。各位的支持是我创作的最大动力。
网上有很多关于中国地图绘制教程，但是关于省市级地图绘制非常少，本推文就是来解决这个问题。下面推文主要以浙江省、温州市为例，使用 leaflet 包绘制省/市级地图。
绘制省级地图
首先，使用regionNames()导入浙江省各市的名字。dem_data可以是读者想要填充在地图上的数据（例如：各市的 GDP，空气质量指数等数据）。使用leafletGeo()创建一个 sp 对象的数据框。
构建好数据集之后，进行绘图。载入高德地图 amap，设置各市边界及颜色并加入图例，得到以下 html 形式图形。
#涂色环节
pal <- colorNumeric(
palette = "Blues",
domain = map$value)
#载入高德地图amap
leaflet(map) %>% amap() %>%
#加入框边界及颜色
addPolygons(stroke = TRUE,
smoothFactor = 2,
fillOpacity = 0.7,
weight = 1,
color = ~pal(value),
popup = ~htmltools::htmlEscape(popup)
) %>%
#加入右下角边框
addLegend("bottomright", pal = pal, values = ~value,
title = "legendTitle",
labFormat = leaflet::labelFormat(prefix = ""),
opacity = 1)
绘制市级地图
绘制市级地图与绘制省级地图类似，只需在regionNames()中进行变化即可，其他几乎相同。下面绘制的温州市地图做了一些小小的拓展。1. 使用真实案例数据；2. 填充颜色变化。数据处理方面主要是用 dplyr 包的full_join()，具体细节就不做过多介绍了。
#======================== 温州市地图绘制 =======================
library(leafletCN)
library(dplyr)
library(leaflet)
#demomap("温州")
#主函数——geojsonmap：地图标色函数
data = read.csv("各区县经营效率.csv",header=T)
data$mean = apply(data[,2:13],1,mean)
colnames(data) = c('region',as.character(2009:2020),"mean")
dim(data)
region = regionNames("温州市")
dat = data.frame(region,runif(length(region)))
data1 = full_join(dat,data[,c(1,8)])
data2 = data1[,-2]
map = leafletGeo("温州市", data2)
#涂色环节heat.colors(6,rev = T)
# c("#FFFF00FF","#FF0000FF","#00FF00FF","#00FFFFFF","#FF00FFFF","#0000FFFF")
pal <- colorNumeric(
palette = c("purple","blue","lightblue","green","yellow","orange",'red'),
domain = map$value)
#载入高德地图amap
leaflet(map) %>% amap() %>%
#加入框边界及颜色
addPolygons(stroke = TRUE,
smoothFactor = 1,
fillOpacity = 1,
weight = 1,
color = ~pal(value),
popup = ~htmltools::htmlEscape(popup)
) %>%
#加入右下角边框
addLegend("bottomright", pal = pal, values = ~value,
title = "效率值",
labFormat = leaflet::labelFormat(prefix = ""),
opacity = 2)
