Rstudio绘制中国地图,代码运行没有错误,但不显示颜色是怎么回事?
x=readShapePoly("C:/Users/AppData/Local/R/win-library/4.2/maptools/shapes/bou2_4p.shp")
library(maptools)
library(sp)
plot(x,col=gray(924:0/924))
getColor=function(mapdata,provname,provcol,othercol){f=function(x,y) return(ifelse (x%in%y,which(y==x),0));colIndex=sapply(mapdata@data$NAME,f,provname);col=c(othercol,provcol)[colIndex+1];return(col);}
provname=c("北京市","天津市","上海市","重庆市")
provcol=c("red","green","yellow","purple")
plot(x,col=getColor(x,provname,provcol,"white"))