0
我想有R讀取一個簡單的景觀分類(貼here)的shape文件和我得到一個錯誤Shapefile環不是有限的,並且不被R讀取 - 如何修復?
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)座標。
現在我知道了解決問題,我更接近解決方案。但是,我不知道如何去做有限的座標。提前致謝。
謝謝。