绘制分省地图

#install.packages("maps")  
#install.packages("mapdata")  
#install.packages("maptools") 
library(maps)  
library(mapdata)  
library(maptools)  
#install.packages("mapproj")
library("mapproj")
library("maptools")
setwd("F:\\R\\地图")  
china_map<-readShapePoly('bou2_4p.shp')  

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

library(ggplot2)
china_map1 <- fortify(china_map)                

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

unique(china_map@data$NAME) 
mydata <- read.csv("F:\\R\\data\\BeiJing.csv")                  
china_data <- join(china_map_data, mydata, type="full")     

midpos <- function(x) mean(range(x,na.rm=TRUE)) 
centres <- ddply(china_data,.(NAME),colwise(midpos,.(long,lat)))

province_city <- read.csv("F:\\R\\地图\\chinaprovincecity.csv")  

ggplot(china_data,aes(long,lat))+
     geom_polygon(aes(group=group,fill=corr1),colour="grey40")+
     scale_fill_gradient(low="white",high="aquamarine4") +  #可自己调整颜色
     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(),
          legend.position=c(0.2,0.3)
          )

颜色参考
效果展示
这里写图片描述这里写图片描述

猜你喜欢

转载自blog.csdn.net/weixin_39777626/article/details/80213395