R2-27 第九期作业-地图绘制

whk 2018-01-02 16:59:55 阅读: 949

setwd("D:/R")

library(maptools)

library(rgdal)

library(ggplot2)

library(plyr)

library(mapproj)

D <- theme_classic() +

theme(panel.grid.major = element_blank(),

panel.grid.minor = element_blank(),

axis.line = element_blank(),

axis.text = element_blank(),

axis.title = element_blank(),

legend.position = "none")

china_map<-readShapePoly("bou2_4p.shp")

data1<-china_map@data

data2<-data.frame(data1,id=seq(0:924)-1)

china_map1<-fortify(china_map)

china_mapdata<-join(china_map1, data2, type = "full")

P<-subset(china_mapdata,NAME=="北京市")

ggplot(P,aes(x=long,y=lat,group=group,fill=NAME))+geom_polygon(fill="yellow" )+geom_path(colour = "grey40")+ggtitle("R2-27-北京-房山区")+geom_text(mapping = aes(x=116,y=39.5,label=NAME),color="red",cex=6)+geom_point(aes(x=115.98,y=39.72),colour="red",size=4,pch=19)+annotate("text",x=116,y=39.8,color="red",label="房山区")+D

P<-theme( panel.background=element_blank(),panel.grid=element_blank(),axis.text=element_blank(),axis.ticks=element_blank(),axis.title=element_blank(),legend.position = c(0.9,0.4))

population<-read.csv("R2-9-1.csv")

mydata<-as.data.frame(population)

pop<-join(china_mapdata, mydata, type = "full")

ggplot(pop,aes(x=long,y=lat,group=group,fill=pop))+geom_polygon()+coord_map("polyconic")+geom_polygon(colour="green")+scale_fill_gradient(low="white",high="red") +ggtitle("R2-27")+P