**用R编辑地图**

编辑地图**"/>

**用R编辑地图**

用R编辑地图
1、地图准备
下载最新的ArcGIS矢量地图数据，这是画地图的基础数据。
下载地址：
2、安装maptools包
install packages(“maptools”)
3、编程
#显示4个直辖市
#设置工作空间
setwd(“D:/R/地图文件”)
#加载安装包
library(maptools);
library(plyr)
library(sp)
#读取数据
x=readShapePoly(‘bou2_4p.shp’)
plot(x);
#绘制数据图形
plot(x,col=gray(924:0/924))
#定义功能块
getColor = function(mapdata, provname, provcol, othercol){
f = function(x, y) ifelse(x %in% y, which(y == x), 0);
colIndex = sapply(mapdata@data$NAME, f, provname);
fg = c(othercol, provcol)[colIndex + 1];
return(fg);
}
#命名
provname=c(“北京市”,“天津市”,“上海市”,“重庆市”)
#定义颜色
provcol=c(“red”,“green”,“yellow”,“purple”)
#绘制图形
plot(x,col=getColor(x,provname,provcol,“white”))