用R画中国地图...

# install.packages("mpas")# install.packages('mapdata')# install.packages('ggmap')# install.packages('plyr')# install.packages('maptools')# library(ggmap)library(maps)library(mapdata)library(sp)library(ggplot2)library(plyr)library(maptools);library(mapproj)setwd('d:/software/workplace')cm=readShapePoly('bou2_4p.shp');dat = read.csv(text = "provname,jd,wd               北京,116.4666667,39.9               上海,121.4833333,31.23333333               天津,117.1833333,39.15               重庆,106.5333333,29.53333333               哈尔滨 ,126.6833333,45.75               长春,125.3166667,43.86666667               沈阳,123.4,41.83333333               呼和浩特,111.8,40.81666667               石家庄 ,114.4666667,38.03333333               太原,112.5666667,37.86666667               济南,117,36.63333333               郑州,113.7,34.8               西安,108.9,34.26666667               兰州,103.8166667,36.05               银川,106.2666667,38.33333333               西宁,101.75,36.63333333               乌鲁木齐,87.6,43.8               合肥,117.3,31.85               南京,118.8333333,32.03333333               杭州,120.15,30.23333333               长沙,113,28.18333333               南昌,115.8666667,28.68333333               武汉,114.35,30.61666667               成都,104.0833333,30.65               贵阳,106.7,26.58333333               福州,119.3,26.08333333               台北,121.5166667,25.05               广州,113.25,23.13333333               海口,110.3333333,20.03333333               南宁,108.3333333,22.8               昆明,102.6833333,25               拉萨,91.16666667,29.66666667               香港,114.1666667,22.3               澳门,113.5,22.2", sep=",")dat$provname = gsub(' ', '', dat$provname)NAME <- c("北京市","广东省","福建省","广西壮族自治区","河北省","湖南省","浙江省","山东省","江苏省","河南省","江西省","辽宁省","香港特别行政区","四川省","湖北省","上海市","山西省","海南省","吉林省","黑龙江省","陕西省","云南省","贵州省","安徽省","内蒙古自治区","台湾省","重庆市","天津市","宁夏回族自治区","甘肃省","新疆维吾尔自治区","西藏自治区","青海省")provname<-c("北京","广东","福建","广西","河北","湖南",  "浙江","山东","江苏","河南","江西","辽宁","香港","四川","湖北","上海","山西","海南","吉林","黑龙江","陕西","云南","贵州","安徽","内蒙古","台湾","重庆","天津","宁夏","甘肃","新疆","西藏","青海")pop=c(1633,1115,6943,3393,2405,4298,2730,3824,1858,7625,      5060,6118,3581,4368,9367,9360,5699,6355,9449,      4768,845,2816,8127,3762,4514,284,3748,2617,      552,610,2095,2296,693);df <- data.frame(NAME, provname, pop)xs <- data.frame(cm@data ,id=seq(0:924)-1)          #含岛屿共925个形状china_map1 <- fortify(cm)           #转化为数据框china_map_data <- join(china_map1, xs, type = "full") china_map <- join(china_map_data, df, type="full")head(china_map)midpos <- function(x) mean(range(x,na.rm=TRUE)) #取形状内的平均坐标centres <- ddply(china_map,.(NAME),colwise(midpos,.(long,lat)))p <- ggplot(china_map,aes(long,lat))+  geom_polygon(aes(group=group,fill=pop),colour="grey60")+  scale_fill_gradient(low="white",high="steelblue") +  coord_map("polyconic") +  # geom_text(aes(x = long,y = lat,label = NAME), data =centres,position = "identity")+  geom_point(data = dat, aes(x=jd, y=wd), colour = "red") +  geom_text(data=dat,aes( x=jd, y=wd), label = dat$provname, size =3) +  theme(    panel.grid = element_blank(),    panel.background = element_blank(),    axis.text = element_blank(),    axis.ticks = element_blank(),    axis.title = element_blank()  )ggsave(p, filename = 'd:/chinamap2.png', width = 16, height = 12)merge(df, dat )