在R中删除国家地图边界之外的数据

我知道这很简单，但无法让它发挥作用。我想删除下面地图上多余的数据点。我该怎么做？ 下面的代码给了我结果。

ggplot() +
geom_polygon(data = rwa2, aes(x = long, y = lat, group= group),
             colour = "black", size = 0.5, fill = "white") +
geom_tile(data = df, aes(x = Lon, y = Lat, z = z, fill = z), alpha = 0.8) +
ggtitle("State Data") +
xlab("Longitude") +
ylab("Latitude") +
scale_fill_distiller(type = "div", palette = "Spectral")+
theme_bw() +
theme(plot.title = element_text(size = 25, face = "bold"),
      legend.title = element_text(size = 15),
      axis.text = element_text(size = 15),
      axis.title.x = element_text(size = 20, vjust = -0.5),
      axis.title.y = element_text(size = 20, vjust = 0.2),
      legend.text = element_text(size = 10)) +
coord_map()

我想删除州边界之外的所有数据。 边界坐标是从 Rdata 文件中获取的readRDS，ID_2 代表州，ID_3 代表区，名称也是如此。请引导我到这里。

由于我们没有您的数据，因此很难处理您的案件。但是，我想给你留下一个方法。只要我从你的代码中看到，你有一个名为df。您想要创建一个临时的 SpatialPointsDataFrame。我们就这样称呼它吧spdf。您还有名为的多边形数据rwa2，这也是一个数据框。如果rwa2来自空间类对象（即 SpatialPolygonsDataFrame），您希望使用它来对保留在多边形内的数据点进行子集化。假设您有一个名为的空间对象sp.rwa2.

library(sp)
library(ggplot2)

步骤1：使用创建一个SpatialPointsDataFramedf确保您分配了 proj4stringsp.rwa2 to spdf。我在这段代码中只有一个示例 proj4string 。

spdf <- SpatialPointsDataFrame(coords = df[, c("Lon", "Lat")], data = df,
                               proj4string = CRS("+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"))

步骤 2：对 sf.rwa2 内的数据点进行子集化。

whatever <- spdf[!is.na(over(spdf, as(sp.rwa2, "SpatialPolygons"))), ]

步骤 3：将 spdf 转换为数据框

whatever <- as.data.frame(whatever)

然后，您将运行 ggplot 代码。rwa2是一个数据框。

ggplot() +
geom_polygon(data = rwa2, aes(x = long, y = lat, group = group),
             colour = "black", size = 0.5, fill = "white") +
geom_tile(data = whatever, aes(x = Lon, y = Lat, fill = z), alpha = 0.8) +
labs(title = "State Data", x = "Longitude", y = "Latitude") +
scale_fill_distiller(type = "div", palette = "Spectral") +
theme_bw() +
theme(plot.title = element_text(size = 25, face = "bold"),
      legend.title = element_text(size = 15),
      axis.text = element_text(size = 15),
      axis.title.x = element_text(size = 20, vjust = -0.5),
      axis.title.y = element_text(size = 20, vjust = 0.2),
      legend.text = element_text(size = 10)) +
coord_map()