R：在经纬度图上覆盖州地图

我已经在经度和纬度的网格上绘制了一个zipcodes图。然而，我想在这幅图上加上宾夕法尼亚州的轮廓。我不知道该怎么做。以下是我目前掌握的代码：

library(tidyverse)
library(zipcode)
data(zipcode)


toplot <- merge(zipcode, nhmzipuse, sort = FALSE, by.x = "zip", by.y = "zip")

ggplot(toplot, aes(longitude, latitude)) + geom_point(aes(color=count))

库（tidyverse）
库（zipcode）
数据（zipcode）
这是另一个很好的起点。答案中的链接提供了一个非常好的开始指南