#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)
)
颜色参考
效果展示