install.packages("maptools")
install.packages("plyr")
library(maptools)
library(ggplot2)
library(plyr)

mytheme<-theme( panel.background=element_blank(),panel.grid=element_blank(),axis.text=element_blank(),axis.ticks=element_blank(),axis.title=element_blank(),legend.position = "none")
china_map<-readShapePoly("bou2_4p.shp")
data<-china_map@data
d<-data.frame(data,id=seq(0:924)-1)
china_map1<-fortify(china_map)
china_mapdata<-join(china_map1, d, type = "full")
BJ<-subset(china_mapdata,NAME=="北京市")
p<-ggplot(BJ,aes(x=long,y=lat,group=group,fill=NAME))+geom_polygon(fill="yellow" )+geom_path(colour = "grey40")+ggtitle("R2-35-北京-海淀")+geom_text(mapping = aes(x=116.7,y=40.5,label=NAME),color="red",cex=6)+geom_point(aes(x=116.17,y=39.97),colour="red",size=4,pch=19)+annotate("text",x=116.15,y=40,color="red",label="海淀")+ mytheme
p

install.packages("mapproj")
library(mapproj)
mytheme<-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))scale_fill_gradient(low="white",high="red")+geom_path(colour = "yellow")+ggtitle("R2-35")+mytheme
population<-read.csv("R2-9-1.csv")
mydata<-as.data.frame(population)
mydata
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-35")+mytheme