我试图将3个不同的空间数据框架叠加在一起 - 到目前为止 - 只有前两个有效 - pts bndry
每当我向混合添加polygon()时，我都会收到以下错误：
> plot(bndry, xlim = c(590900, 595900), ylim = c(4785786, 4793200))
> points(pts, col=colcode, pch=16)
> polygon(polys)
Error in as.double(y) :
cannot coerce type 'S4' to vector of type 'double'
确保pts和polys投射在完全相同的CRS中：
> polys@proj4string
CRS arguments:
+init=epsg:32617 +proj=utm +zone=17 +datum=WGS84 +units=m
+no_defs +ellps=WGS84 +towgs84=0,0,0
> pts@proj4string
CRS arguments:
+init=epsg:32617 +proj=utm +zone=17 +datum=WGS84 +units=m
+no_defs +ellps=WGS84 +towgs84=0,0,0
这个错误是否有原因以及已知的解决方法？