问题描述
我想为我的点绘制不同的形状，并将它们相对于 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解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）