R2-14-第九次地图

lucyylx 2018-01-02 22:56:21 阅读: 1362

library(maptools)

library(rgdal)

library(ggplot2)

library(plyr)

china_map=rgdal::readOGR(“/Users/LUCY/Desktop/科研狗R语言第二期/2-9/china-province-border-data/bou2_4p.shp”)

china_map1<-fortify(china_map)

x<-china_map@data

xs<-data.frame(x,id=seq(0:924)-1)

china_map_data <- join(china_map1, xs, type = "full")

hunan<-subset(china_map_data,NAME=="湖南省")

map1<gplot(hunan,aes(x=long,y=lat,group=group,fill=NAME))+geom_polygon(fill="beige")+geom_path(colour="grey40")+geom_point(x=112.58,y=28.11,fill="black")+annotate("text",x=112.7,y=28,label="长沙市")+annotate("text",x=111,y=29,label="湖南省",size=8,colour="red")+ggtitle("R2-14-湖南省-长沙市")

mytheme<- theme(panel.grid = element_blank(), panel.background = element_blank(),axis.text = element_blank(),axis.ticks = element_blank(),axis.title = element_blank())

map1+mytheme

library(mapproj)

mypop<-read.csv("/Users/LUCY/Desktop/科研狗R语言第二期/2-9/map.csv")

china_data<-join(china_map_data,pop,type="full")

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

map2+mytheme