我想让R读取简单景观分类的shapefile(发布here)并收到错误
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)坐标。
既然我知道这个问题,我就更接近解决方案了。但是,我不知道如何使坐标有限。提前谢谢。
谢谢。
答案 0 :(得分:-1)
通常使用参数'force_ring = TRUE',错误不再出现。只尝试使用该参数,它对我有用。
问候!