Tôi đang tạo một bản đồ hạt choropleth với đường viền màu xám và tôi cũng muốn bao gồm các ranh giới trạng thái màu đen. Có ai biết làm thế nào tôi đi về việc thêm một lớp thứ hai của bản đồ nhà nước đến một bản đồ quận hiện có?ggplot2 lập bản đồ các ranh giới hạt trong một màu và các ranh giới trạng thái khác trên cùng một bản đồ
Đây là data set và mã tôi đã kết thúc bằng:
#load libraries
library(ggplot2)
library(ggmap)
library(maps)
library(plyr)
#get wif file
wip <- read.csv("wip.csv")
#get map data for US counties and states
county_map <- map_data("county")
state_map <- map_data("state")
#merge wip and county_map
wip_map <- merge(county_map, wip, by.x=c("region", "subregion"),
by.y=c("region","subregion"), all.x=TRUE)
#resort merged data
wip_map <- arrange(wip_map, group, order)
#relpace NA with 0's
wip_map[is.na(wip_map)] <- 0
#generate a disctrete color pallette
pal <- c("#F7FCF5","#74C476","#41AB5D","#238B45","#006D2C","#00441B")
theme_clean <- function(base_size = 12) {
require(grid)
theme_grey(base_size) %+replace%
theme(
axis.title = element_blank(),
axis.text = element_blank(),
panel.background = element_blank(),
panel.grid = element_blank(),
axis.ticks.length = unit(0,"cm"),
axis.ticks.margin = unit(0,"cm"),
panel.margin = unit(0,"lines"),
plot.margin = unit(c(0,0,0,0),"lines"),
complete = TRUE
)
}
final_map <- ggplot(wip_map, aes(x=long, y=lat, group=group, fill=factor(CATEGORY))) +
geom_polygon(colour="grey", aes(fill=factor(CATEGORY))) +
scale_fill_manual(values=pal) +
expand_limits(x = wip_map$long, y = wip_map$lat) +
coord_map("polyconic") +
labs(fill="Number Per\nCounty") +
theme_clean()
final_map + geom_path(data = state_map , colour = "red")
Cảm ơn!
Các tập tin mà bạn liên kết cho bộ dữ liệu của bạn trong thùng rác của Google Drive của bạn. Nó sẽ biến mất sớm thôi. Bạn có thể cập nhật nội dung cho một liên kết vĩnh viễn hơn được lưu trữ trên một thứ gì đó như pastebin không? – csalvato