ggplot2在一个颜色中映射县边界,在同一个地图上映射另一个颜色和州边界

时间:2013-05-19 04:50:49

标签: r ggplot2

我正在创建一个带有灰色边框的等值区县地图,我还希望将州边界包含在黑色中。有谁知道如何在现有的县地图中添加第二层状态映射?

以下是data set以及我最终使用的代码:

#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")

谢谢!

1 个答案:

答案 0 :(得分:12)

只需在代码中添加geom_path ...

我使用红色突出显示边界,但您可以轻松地将其设置为黑色。

    ggplot( wip_map, aes( x = long , y = lat , group=group ) ) +
        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( ) +
        geom_path( data = state_map , colour = "red")

enter image description here