用R画中国地图
准备工作
- R版本:4.0.5
- RStudio
- 中国国界线文件,各省市的经纬度文件
链接:https://pan.baidu.com/s/1ANa4rMZxuXLgle-j-mruVw
提取码:1010 - 代码
install.packages("maptools")
install.packages("rgdal")
library(maptools)
library(rgdal)
str.crs <- " +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"
china <- readShapePoly("M:/temp/geo/china_border/bou2_4p.shp", proj4string = CRS(str.crs))
par(mar = c(2, 2, 0.1, 0.1))
plot(china,col=c("peachpuff","papayawhip","tan1","tan2","tomato"),
ylim = c(18, 54), panel.first = grid());
city=read.csv("M:/temp/geo/china_city/city.csv",header=T, sep=",",encoding = "UTF-8")
points(city$jd,city$wd,pch=19,col="black")
text(city$jd, city$wd, city[,1], cex = 0.6, col = rgb(0,0, 0, 0.7), pos = c(2, 4, 4, 4, 3, 4, 2, 3, 4, 2, 4, 2, 2,
4, 3, 2, 1, 3, 1, 1, 2, 3, 2, 2, 1, 2, 4, 3, 1, 2, 2, 4, 4, 2))
结果显示
![在这里插入图片描述](https://img-blog.csdnimg.cn/20210412161353639.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L1NhbmR5eno=,size_16,color_FFFFFF,t_70)