服务热线
178 0020 3020
R2-06 地图作业
part1
library(maptools)
library(ggplot2)
library(plyr)
library(rgdal)
china_map<-readShapePoly("bou2_4p.shp")
jingweidu<-china_map@data
xingzhengquyu<-data.frame(jignweidu,id=seq(0:924)-1)
china_map1<-fortify(china_map)
china_mapdata<-join(china_map1, xingzhengquyu, type = "full")
hebei<-subset(china_mapdata,NAME=="河北省")
ggplot(hebei, aes(x = long, y = lat, group = group,fill=NAME))+
geom_polygon(fill="yellow" )+
geom_path(colour = "grey40")+
ggtitle("R2-06-dean-SJZ")+
geom_text(mapping = aes(x=117,y=35,label=NAME),color="red",cex=5)+
geom_point(aes(x=114.48,y=38.03),colour="blue",size=1.8)+
annotate("text",x=115.48,y=38.03,label="石家庄")+
coord_map("polyconic")
part2
pop<-read.table(file.choose(),header = T,sep = ",")
china_pop<-join(china_mapdata, pop, type = "full")
ggplot(china_pop, aes(x = long, y = lat, group = group, fill=pop))+
geom_polygon(colour="green")+
scale_fill_gradient(low="white",high="red") +
coord_map("polyconic") +
ggtitle("R2-06-dean")
附件