1、
library(maps)
library(mapdata)
par(mai = c(0,0,0,0),mgp = c(3,1,0))
map("china",col = "black",xlim=c(70,140),ylim = c(15, 60),
lwd=2.5,panel.first = grid())
axis(1,lwd=3,font = 2,cex.axis=1.5)
axis(2,lwd=3,font = 2,cex.axis=1.5)
axis(3,lwd=3,font = 2,cex.axis=1.5)
axis(4,lwd=3,font = 2,cex.axis=1.5)
box(which = "plot",lwd=3)
library(REmap)
options(remap.ak = "输入个人秘钥")
city_vec = c("Beijing","Xian","Haerbin","Shanghai","Guangzhou")
pos <- get_geo_position (city_vec)
pos
points(pos$lon,pos$lat, pch= c(15,16,17,18,25),bg = "purple",cex=2,
col=c("red","cyan","blue","green","purple"))
legend("bottomright",c("aa","bb","cc","dd","ee"),inset = 0.04,
pch = c(15,16,17,18,25),cex=1.8,
col=c("red","cyan","blue","green","purple"),
pt.bg ="purple",box.lwd =3)
原文:https://www.cnblogs.com/liujiaxin2018/p/14323053.html