0

我正在尝试让 R 读取一个简单景观分类的 shapefile(在此处发布)并且出现错误

shape1 <-readShapePoly("JustForest_BahiaSul_1990_Test.shp", force_ring=TRUE, delete_null_obj=TRUE)
Error in Polygon(coords = crds) : ring not closed

当我尝试使用 rgdal 作为替代方案时,同样的错误结果:

shape1 <-readOGR("JustForest_BahiaSul_1990_Test.shp", layer="JustForest_BahiaSul_1990_Test")OGR data source with driver: ESRI Shapefile
Source: "JustForest_BahiaSul_1990_Test.shp", layer: "JustForest_BahiaSul_1990_Test"
with 19122 features and 6 fields Feature type: wkbPolygon with 2 dimensionsError in stopifnot(is.list(srl)) : ring not closed

一位非常有帮助的 R 同事指出,R 没有正确读取 shapefile 的原因是 shapefile 在实体 15755 中具有无限 (NaN) 坐标。

现在我知道了问题所在,我离解决方案更近了一点。但是,我不知道如何使坐标有限。提前致谢。

谢谢你。

4

1 回答 1

-1

通常使用参数“force_ring=TRUE”,错误不再出现。仅尝试使用该参数,它对我有用。

问候!

于 2013-07-31T00:27:32.533 回答