2013-05-07 40 views
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)座標。

現在我知道了解決問題,我更接近解決方案。但是,我不知道如何去做有限的座標。提前致謝。

謝謝。

回答

-1

通常使用參數'force_ring = TRUE',錯誤不會再出現。只嘗試使用該參數,它爲我工作。

問候!

相關問題