地理绘图网与R程序.docVIP

国家基础地理信息系统：/。（提供国家SHP文件数据下载服务） /091019/5920591.html提供mapinfo数据下载，经过mapinfo软件可得到各省份SHP数据，可研究其它各省地理绘图。 /地信网论坛，提供软件学习和省份SHP数据。 绘图代码： library(maptools) x=readShapePoly(bou2_4p.shp)###读入地理数据 plot(x) 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); col=c(othercol,provcol)[colIndex+1]; return(col); } qgx-c(6859.6,6791.004,6852.123036,7078.243096,7191.494986,7299.367411,7627.838944,8054.997925,7982.502944) qgx1-c(6859.6,7780.606061,8481.444775,9130.581569,10008.73712,11050.97451,12397.35209,13438.8583,14758.68283 ) qgx-round(qgx,digits=2) qgx1-round(qgx1,digits=2) shuju-read.table(2009全国可支配收入数据.csv,sep=,,header = TRUE) as.character(na.omit(unique(x$NAME)));按各省顺序给经济数据赋值 as.character(na.omit(unique(x1$NAME)));按各省顺序给经济数据赋值 provname=c(北京市,天津市,河北省,山西省,内蒙古自治区, 辽宁省,吉林省,黑龙江省,上海市,江苏省, 浙江省,安徽省,福建省,江西省,山东省, 河南省,湖北省,湖南省,广东省, 广西壮族自治区,海南省,重庆市,四川省,贵州省, 云南省,西藏自治区,陕西省,甘肃省,青海省, 宁夏回族自治区,新疆维吾尔自治区,台湾省, 香港特别行政区); #################### ##按各年自然划分为基础## #################### ####2001年 pop-shuju[,2] provcol=rep(white,33); for(i in 1:33){if(pop[i]0)provcol[i]=red} for(i in 1:33){if(pop[i]quantile(shuju[,2], probs=seq(0,1,0.2))[[2]])provcol[i]=yellow} for(i in 1:33){if(pop[i]quantile(shuju[,2], probs=seq(0,1,0.2))[[3]])provcol[i]=green} for(i in 1:33){if(pop[i]quantile(shuju[,2], probs=seq(0,1,0.2))[[4]])provcol[i]=blue} for(i in 1:33){if(pop[i]quantile(shuju[,2], probs=seq(0,1,0.2))[[5]])provcol[i]=purple} provcol brks-round(quantile(shuju[,2], probs=seq(0,1,0.2)), digits=2) colours-c(red,yellow,green,blue,purple) plot(x,col=getColor(x,provname,provcol,white),xlab=,ylab=); legend(x=92,y=59,legend=leglabs(brks),fill=colours, bty=n)####图例绘图 title(main = list(2001年各省份城镇可支配收入, cex=1.5, col=red, font=3)) ##图标##用R画中国地图并标注城市位置 par(mar=rep(0,4)) dat = read.csv(textConnection(城市,jd,wd 北 京,118.4666667,39.9 上 海,118.4833333,31 天 津,116.5833333,38.65