'china'数据库的经纬度数据把全国分成了712个regions,麻烦的就是怎样把这712个polygon对应到31个省市上去(实际上是30个),而且重庆的问题仍然不好解决……除非有它的地理数据
x=map('china',col='lightgray')<br />
dat=cbind(x$x,x$y)<br />
y=c(1,which(is.na(x$x)),length(x$x))<br />
for(i in 2:length(y)){<br />
lines(dat[y[i-1]:y[i],])<br />
Sys.sleep(1) #or locator(1)<br />
}
要是有足够的时间和耐心,其实也可以手工把各省的数据“挑”出来:上面的程序每作一条线,都手工记下来这条线是哪(两)个省的边界,最后把712组数据分好类,得到每个省的经纬度数据。