所以我试图在地球上绘制一堆坐标并跟踪每个国家有多少个坐标。我已经很好地绘制了地图和坐标,但是当我尝试使用交集来计算每个国家(多边形)内有多少个坐标时,会导致错误。我尝试使用 st_make_valid 函数来修复地球形状文件,但它弄乱了几何形状。我是 R 的新手,因此我们将不胜感激。
我使用以下代码来绘制地球形状文件和顶部的坐标:
library(tidyverse)
library(sf)
library(rmapshaper)
library(rnaturalearth)
library(rnaturalearthdata)
library(sp)
library(raster)
###############
# Load Data
###############
# Read in data from .csv file
MeteoriteData <- read.csv("C:/Users/ChaseDickson_/Desktop/College/AERO 689/Semester Project/Meteorite Landings.csv")
# Convert these points to an SF object, specifying the X and Y
# column names, and supplying the CRS as 4326 (which is WGS84)
MeteoriteData.sf <- st_as_sf(MeteoriteData, coords=c('long', 'lat'), crs=4326)
world <- (ne_countries(scale = "medium", returnclass = "sf"))
MeteoriteMap <- ggplot(data = world) +
geom_sf() +
geom_sf(data = …Run Code Online (Sandbox Code Playgroud)