使用ggplot在R中绘制组合形状文件时出现问题，未给出定义的错误

使用的图书馆：

library(sp)
library(sf) 
library(ggplot2) 
library(ggmap)

创建了一个名为“coordinate.data”的数据框，其中经度和纬度作为列名，气象站位置作为行名

longitude <- c(-73.964482,-73.953678,-73.893522,-73.815856,-74.148499)
latitude <- c(40.767544,40.631762,40.872481,40.734335,40.604014)
coordinate.data <- data.frame(longitude,latitude)
rownames(coordinate.data) <- c("MANH","BKLN","BRON","QUEE","STAT")

---

经度我不太确定你从那些形状文件中删除了什么。我没有删除任何内容。我也没有结合任何东西。在我的输出中有四个独立的投影层：底图、NY形状文件、NJ形状文件，以及下面创建和描述的名为sites的形状文件

此地图是使用3个现有形状文件和在以下步骤中创建的场地形状文件创建的。打印并显示从每个步骤创建的变量以及每个步骤的说明

首先，创建空间几何体变量：

MULTIPOINT ((-73.96448 40.76754), (-73.95368 40.63176), (-73.89352
 40.87248), (-73.81586 40.73434), (-74.1485 40.60401))

然后，创建几何体列（将多点变量与crs组合）

然后创建一个由一列或一列站点名称组成的常规df

然后创建简单要素对象-（将几何体列与df组合）。输出显示：

A sf: 5 × 2
site    pts.sfc
<fct>   <MULTIPOINT [m]>
MANH    MULTIPOINT ((-73.96448 40.7...
BKLN    MULTIPOINT ((-73.96448 40.7...
BRON    MULTIPOINT ((-73.96448 40.7...
QUEE    MULTIPOINT ((-73.96448 40.7...
STAT    MULTIPOINT ((-73.96448 40.7...

输出显示了纽约州的5个站点以及新泽西州和纽约州。NY形状文件的包容性似乎较低，仅包含一些有限的NY数据

要使整个映射正确，需要按顺序执行四个步骤。几何体、几何体列、sf对象、写入形状文件夹

除了一些小问题，您的代码在我的机器上运行良好。
mesonet.map <-ggplot() +
  ggtitle("NY Mesonet Site Locations") +
  xlab("Longitude") +
  ylab("Latitude") +
  geom_point(data=coordinate.data,aes(x=longitude,y=latitude))+
  geom_text(aes(x=longitude,y=latitude,label=rownames(coordinate.data)),size=3.25,nudge_y=0.02)+
  geom_sf(data=combined.data,fill='darkgreen') +

mesonet.map + theme(
  panel.background=element_rect(fill="lightblue",color="lightblue",size=0.5,linetype="solid"),
  panel.grid.major=element_line(size=0.5,linetype='solid',color="white"), 
  panel.grid.minor=element_line(size=0.25,linetype='solid',color="white")
)

MULTIPOINT ((-73.96448 40.76754), (-73.95368 40.63176), (-73.89352
 40.87248), (-73.81586 40.73434), (-74.1485 40.60401))

MULTIPOINT ((-73.96448 40.76754), (-73.95368 40...

Geometry set for 1 feature 
geometry type:  MULTIPOINT
dimension:      XY
bbox:           xmin: -74.1485 ymin: 40.60401 xmax: -73.81586 ymax: 40.87248
projected CRS:  WGS 84 / Pseudo-Mercator

A sf: 5 × 2
site    pts.sfc
<fct>   <MULTIPOINT [m]>
MANH    MULTIPOINT ((-73.96448 40.7...
BKLN    MULTIPOINT ((-73.96448 40.7...
BRON    MULTIPOINT ((-73.96448 40.7...
QUEE    MULTIPOINT ((-73.96448 40.7...
STAT    MULTIPOINT ((-73.96448 40.7...

Writing layer `Weather_sites' to data source `C:/pathR' using driver `ESRI Shapefile'
Writing 5 features with 1 fields and geometry type Multi Point.