评论 - 1129使用R绘制中国地图
由于绘制地图的包不是R的基础包，所以要下载安装：
install.packages("maps");
install.packages("mapdata")
install.packages("maptools")
library（maps）
library(mapdata)
map("china")
即可载入中国地图；
此处载入的中国地图版本较老，四川省还包括重庆。
******************
library(maptools);
x=read.shape('bou2_4p.shp');#下文中会继续用到x这个变量，
#如果你用的是其它的名称，
#请在下文的程序中也进行相应的改动。
新版本的maptools包不再提供read.shape()函数，请用readShapePoly()代替。
*******************************************************************
plot(x,fg=gray(924:0/924));
新版本的maptools包的绘图参数也有所改变，请将fg换成col。
++++++++++++++++++++++++++++++++++++
如何获取某一个特定地区的ID，进而设置我们想要的颜色？事实上，
在变量x中，就已经存储了我们想要的信息。在R中输入“x[[2]]”或
“x$att.data”，会得到一个925行7列的数据框，这其实是
bou2_4p.dbf这个文件中存储的信息，之前的read.shape()函数虽然读
取的是bou2_4p.shp文件，但在默认情况下会把dbf文件的信息也放到变
量之中。对于这个数据框，其行名就是每一个区域的ID编号，第一列和
第二列分别是面积和周长，最后一列是该区域所属的行政区名，其它的
列应该也是一些编号性质的变量。于是，通过查找相应的行政区对应的
行名，就可以对fg参数进行赋值了。
**************************************************************
getColor=function(mapdata,provname,provcol,othercol)
f=function(x,y) ifelse(x %in% y,which(y==x),0);
colIndex=sapply(mapdata$att.data$NAME,f,provname);
fg=c(othercol,provcol)[colIndex+1];
return(fg);
地图数据的组织形式有所变化，上面函数中的mapdata$att.data$NAME需要替换为。
**********************************************************
provname=c("北京市","天津市","上海市","重庆市");
provcol=c("red","green","yellow","purple");
plot(x,fg=getColor(x,provname,provcol,"white"));
++++++++++++++++++++++++++++++
其中mapdata是存放地图数据的变量，在上面的例子中就是x，
provname是需要改变颜色的地区的名称，provcol是对应于
provname的代表颜色的向量（名称和数字均可），othercol是其
它地区的颜色。
*************************************************************
provname=c("北京市","天津市","上海市","重庆市");
provcol=c("red","green","yellow","purple");
plot(x,fg=getColor(x,provname,provcol,"white"));
as.character(na.omit(unique(x$att.data$NAME)));
+++++++++++++++++++++++++++++++++++++++++++++++++
由此生成的向量有33个元素，少了澳门特别行政区，
这是这个数据中的一块瑕疵。在x$att.data的第899行有一个NA，
不知道它代表的是否就是澳门。
*****************************************************
provname=c("北京市","天津市","河北省","山西省","内蒙古自治区",
&&"辽宁省","吉林省","黑龙江省","上海市","江苏省",
&&"浙江省","安徽省","福建省","江西省","山东省",
&&"河南省","湖北省","湖南省","广东省",
&&"广西壮族自治区","海南省","重庆市","四川省","贵州省",
&&"云南省","西藏自治区","陕西省","甘肃省","青海省",
&&"宁夏回族自治区","新疆维吾尔自治区","台湾省",
&&"香港特别行政区");
pop=c(43,98,58,7625,
&&81,60,49,
&&6,14,284,,
&&552,610,3);
provcol=rgb(red=1-pop/max(pop)/2,green=1-pop/max(pop)/2,blue=0);
plot(x,fg=getColor(x,provname,provcol,"white"),xlab="",ylab="");
++++++++++++++++++++++++++++++++++++++++++++++++++++
在绘制地图的过程中，还有一个比较有用的参数是recs，它是一个
由多边形ID组成的向量，表示在地图中只画出这些ID所代表的区域。
利用这个参数，就可以画出某一部分的地图，例如下面的例子是我国
中部六省的地图：
++++++++++++++++++++++++++++++++++++++
getID=function(mapdata,provname)
index=mapdata$att.data$NAME %in%
ids=rownames(mapdata$att.data[index,]);
return(as.numeric(ids));
midchina=c("河南省","山西省","湖北省","安徽省","湖南省","江西省");
plot(x,recs=getID(x,midchina),fg="green",ol="white",xlab="",
&&ylab="");
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
新版本的maptools包的绘图函数已经取消了recs这个参数，现在要实现这个功能，
可以在颜色上把不需要的省份变成白色，其中填充色用col参数，
边界颜色用border参数。例如上面的例子可以用下面的函数来实现：
###################################################
plot(x, col = getColor(x, midchina, rep("green", 6),
& & "white"), border = "white",
xlab = "", ylab = "")
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
最后要说的是，在画出的图上仍然可以用points()函数和
text()函数加上点和文字，而maptools包中还提供了一个
pointLabel()函数，用来解决文本标签的重叠问题。
这部分内容请参阅博文：用R画中国地图并标注城市位置，
以及避免文本标签重叠：maptools中的pointLabel()。
______________________________________________
par(mar=rep(0,4))
dat = read.csv(textConnection("城市,jd,wd
津,117..15
哈尔滨,126..75
阳,123.4,41.
呼和浩特,111.8,40.
石家庄,114..
& & 济 南,117,36.
& & 郑 州,113.7,34.8
安,108.9,34.
州,103..05
宁,101.75,36.
& & 乌鲁木齐,87.6,43.8
& & 合 肥,117.3,31.85
州,120.15,30.
& & 长 沙,113,28.
汉,114.35,30.
都,104..65
阳,106.7,26.
州,119.3,26.
北,121..05
州,113.25,23.
& & 昆 明,102.
& & 澳门,113.5,22.2"))
library(maps)
library(mapdata)
map("china", col = "darkgray", ylim = c(18, 54), panel.first =
points(dat$jd, dat$wd, pch = 19, col = rgb(0, 0, 0, 0.5))
text(dat$jd, dat$wd, dat[, 1], cex = 0.9, col = rgb(0,
& & 0, 0, 0.7), pos = c(2, 4, 4,
4, 3, 4, 2, 3, 4, 2, 4, 2, 2,
& & 4, 3, 2, 1, 3, 1, 1, 2, 3, 2,
2, 1, 2, 4, 3, 1, 2, 2, 4, 4, 2))
axis(1, lwd = 0); axis(2, lwd = 0); axis(3, lwd = 0); axis(4, lwd =
***********************************************
已投稿到：
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。公司需要做一个软件，就是允许用户在一张地图上使用鼠标圈划将地图划分成若干区域，请问怎么做
[问题点数：20分，结帖人CSDN]
公司需要做一个软件，就是允许用户在一张地图上使用鼠标圈划将地图划分成若干区域，请问怎么做
[问题点数：20分，结帖人CSDN]
不显示删除回复
显示所有回复
显示星级回复
显示得分回复
只显示楼主
相关帖子推荐：
本帖子已过去太久远了，不再提供回复功能。浅析城市影像地图与普通线划地图在设计制作上的差异_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
评价文档：
&&¥2.00
&&¥0.50
&&¥2.00
&&¥2.00
&&¥1.00
喜欢此文档的还喜欢
浅析城市影像地图与普通线划地图在设计制作上的差异
城​市​影​像​地​图​ ​ ​普​通​线​划​地​图
阅读已结束，如果下载本文需要使用
想免费下载本文？
把文档贴到Blog、BBS或个人站等：
普通尺寸(450*500pix)
较大尺寸(630*500pix)
你可能喜欢}