library(maptools)
axx<-readShapePoly("bou2_4p.shp", proj4string=CRS("+proj=longlat +ellps=clrk66"))
x=axx
getColor=function(mapdata,provname,provcol,othercol)
{
f=function(x,y) ifelse(x %in% y,which(y==x),0);
colIndex=sapply(mapdata$NAME,f,provname);
col=c(othercol,provcol)[colIndex+1];
return(col);
}
provname=c("湖北省","甘肃省","河北省","辽宁省","黑龙江省","内蒙古自治区");
provcol=c("green","green","green","green","green","green","green");
plot(x,col=getColor(x,provname,provcol,"white"))
怎么在各个省加上省的名字上去