如何在光栅上绘制分组点和不同形状？

如何解决如何在光栅上绘制分组点和不同形状？

我想为我的点绘制不同的形状，并将它们相对于 year 这样的整数分组，并用图例显示这一点。

我可以用点绘制地图并操纵 pch，尽管我不知道如何将这些与年份相关联，这样在使用 {{1} 时，图例中会显示带有点的年份} 在 points 上。

这是带点的地图：

plot

获取地图的代码如下，对于范围和 CRS 投影，可以从上面的栅格层复制（因为我有一个栅格来执行此操作）：

class      : RasterLayer 
dimensions : 1637,3203,5243311  (nrow,ncol,ncell)
resolution : 2280,2060  (x,y)
extent     : -2821420,4481420,-1449450,1922770  (xmin,xmax,ymin,ymax)
crs        : +proj=aea +lat_0=37.5 +lon_0=-96 +lat_1=29.5 +lat_2=45.5 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs 
source     : memory
names      : layer 
values     : 1,1  (min,max)

我尝试过的：

us.map <- getData('GADM',country='USA',level=1)
us.map.1 <- st_as_sf(us.map)

这是在地图上模仿 plot(us.map.1) points(xy[,1:2],add=TRUE,pch=c(1,2)[xy$year]) 的可重现代码：

black points