如何在 R 传单中将 SP 对象绘制为 SF

Question

我刚开始使用 GIS，所以请不要吝啬。

我正在尝试绘制不同的 shp 文件，但是这个特定的文件没有绘制到我使用传单包构建的交互式地图上。

文件本身可以很好地加载到 R 中，但是当我将它插入我的地图时，它看起来超级缩小并且看起来像 this

所以我运行了这段代码“属性（几何）” 这出现了-

$generic
[1] "geometry"
attr(,"package")
[1] "sp"

$package
[1] "sp"

$group
list()

$valueClass
character(0)

$signature
[1] "obj"

$default
`\001NULL\001`

$skeleton
(function (obj) 
stop("invalid call in method dispatch to 'geometry' (no default method)", 
    domain = NA))(obj)

$class
[1] "standardGeneric"
attr(,"package")
[1] "methods"

在我看来，我正在处理“sp”而不是“sf”.. 需要对此做些什么吗？

另外这部分“stop(”invalid call in method dispatch to 'geometry' (no default method)”看起来也很奇怪。

你能帮我解决这个问题吗？

我正在尝试使用传单中的“addPolygons”函数进行绘图。

非常感谢！

Answer 1

您似乎使用rgdal::readOGR() 来读取shapefile；它产生一个Spatial 对象（在 {sp} 领域，而不是 {sf} 领域）。

我建议您考虑使用sf::st_read() 将 shapefile 导入 R；那么它应该表现得更好。