标签：sheng 技术 join axis new 加载 使用 for word
1、绘制中国地图
library("maptools")
china_map = readShapePoly("D:/Astatistics/R/画热力图数据材料/中国省级地图空间数据文件/bou2_4p.shp") # 读取地图空间数据
plot(china_map)#可以不画
2*、改变投影方式
library(ggplot2)
ggplot(china_map,aes(x=long,y=lat,group=group)) +
geom_polygon(fill="white",colour="grey") +
coord_map("polyconic")
3、加载和拼接行政信息
x <- china_map@data #读取行政信息
xs <- data.frame(x,id=seq(0:924)-1) #含岛屿共925个形状
library(ggplot2)
china_map1 <- fortify(china_map) #转化为数据框
library(plyr)
china_map_data <- join(china_map1, xs, type = "full") #合并两个数据框
4、制作业务工作表，保存为csv。需要NAME行（一定命名为NAME)，和绘图指标C1
其中，查看数据包里各省名称
unique(china_map@data$NAME)
5、读取业务指标，与地理数据合并
mydata <- read.csv("c:/rstudy/geshengzhibiao.csv") #读取指标数据，csv格式
china_data <- join(china_map_data, mydata, type="full") #合并两个数据框
提示：Joining by: NAME
6、绘制地图
7*、计算经纬度平均值添加标签
midpos <- function(x) mean(range(x,na.rm=TRUE)) #取形状内的平均坐标
centres <- ddply(china_data,.(province),colwise(midpos,.(long,lat)))
ggplot(china_data,aes(long,lat))+ #此处语法与前面不同，参考ggplot2一书P85
geom_polygon(aes(group=group,fill=zhibiao),colour="black")+
scale_fill_gradient(low="white",high="steelblue") +
coord_map("polyconic") +
geom_text(aes(label=province),data=centres) +
theme(
panel.grid = element_blank(),
panel.background = element_blank(),
axis.text = element_blank(),
axis.ticks = element_blank(),
axis.title = element_blank()
)
8*、读取省会城市坐标添加标签
province_city <- read.csv("D:/Astatistics/R/画热力图数据材料/省会坐标.csv") #读取省会城市坐标
ggplot(china_data,aes(long,lat))+
geom_polygon(aes(group=group,fill=zhibiao),colour="grey60")+
scale_fill_gradient(low="white",high="steelblue") +
coord_map("polyconic") +
geom_text(aes(x = jd,y = wd,label = province), data =province_city)+
theme(
panel.grid = element_blank(),
panel.background = element_blank(),
axis.text = element_blank(),
axis.ticks = element_blank(),
axis.title = element_blank()
)
标签：sheng 技术 join axis new 加载 使用 for word
原文地址：http://www.cnblogs.com/EXODUS1917/p/7581521.html
