r - 如何将州和县线添加到 R 中的 geom_raster ggplot？

Question

我在加利福尼亚州南部有一个数据集(纬度、经度、值)，我可以使用 ggplot(见下文)绘制它。纬度范围从 31.5 到 35.5，经度范围从 -121 到 -113。 ggplot map 看起来不错，但我想将州和县边界线添加到我的 map 中。

我曾尝试添加 geom_polygon() 但我一直在获取整个美国西南部的边界，而不是仅在我的地块区域(加利福尼亚州南部)。

我想知道是否有人会对如何做到这一点有任何建议/想法。非常感谢。

“df 是南加州的 3D 数据框(纬度、经度、值)”
“纬度范围从 31 到 35”
“lon 范围从 -121 到 -113”
“值范围从 3000 到 5000”

library(maps)
library(fields)
library(ggplot2)    
state_map <- map_data("state")
bmap=map_data("state", xlim=c(-121,-113), ylim=c(31.5,35.5))

p1 = ggplot() + geom_raster(data = df, aes(x=Lon, y = Lat, fill=Value)) + 
     coord_fixed(ratio = 1) +
     geom_polygon(data=bmap,aes(x=long,y=lat,group=group), inherit.aes=F, 
     colour='black', fill=NA) +      
     scale_fill_gradientn(na.value="white",limits = c(3500,5000), 
                          colours=c("yellow","orange","green","blue")) +
     theme(panel.background = element_rect(fill = 'white',color="black"),
           legend.key = element_blank()) 
print(p1)

我正在获得整个西南地区的边界，但我只希望我的数据区域:

Answer 1

您可以使用 coord_cartesian 缩放 map :

library(maps)
library(fields)
library(ggplot2)    
state_map <- map_data("state")
bmap=map_data("state", xlim=c(-121,-113), ylim=c(31.5,35.5))

# Create the df dataset with random numbers for Values
n <- 80
df <- expand.grid(seq(31,35,length.out=n),
            seq(-121,-113,length.out=n))
df$Vaule <- 3000+2000*runif(nrow(df))
names(df) <- c("Lat","Lon","Value")

p1 = ggplot() + 
     coord_fixed(ratio = 1) +
     geom_raster(data = df, aes(x=Lon, y = Lat, fill=Value), alpha=0.2) +   
     geom_polygon(data=bmap,aes(x=long,y=lat,group=group), inherit.aes=F, 
     colour='black', fill=NA, lwd=1) +      
     scale_fill_gradientn(na.value="white",limits = c(3500,5000), 
                          colours=c("yellow","orange","green","blue")) +
     scale_y_continuous(expand = c(0,0)) + 
     scale_x_continuous(expand = c(0,0)) + 
     theme(panel.background = element_rect(fill = 'white',color="black"),
           legend.key = element_blank()) +
     coord_cartesian(xlim=c(-121, -113), ylim=c(31,35)) +
     theme(panel.border = element_rect(colour = "gray50", fill=NA, size=1))
print(p1)