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

library(maptools)
library(rgeos)
library(sp)

world_map <- maptools::getMap(resolution = "high")

world_map_sp <- maptools::map2SpatialPolygons(world_map, IDs = world_map$names)

bbox <- bbox(world_map_sp)

# 假设你有一个数据框df，其中包含经度(lon)和纬度(lat)信息
df <- data.frame(lon = c(120, 130, 140), lat = c(30, 40, 50))

# 创建一个空的逻辑向量，用于存储数据点是否在边界框内的信息
in_bbox <- logical(nrow(df))

# 判断数据点是否在边界框内，并将结果存储在in_bbox中
for (i in 1:nrow(df)) {
  point <- SpatialPoints(matrix(c(df$lon[i], df$lat[i]), ncol = 2))
  in_bbox[i] <- gIntersects(point, bbox, byid = TRUE)
}

# 保留在边界框内的数据点
df_filtered <- df[in_bbox, ]