Xóa một số đa giác nhất định khỏi SpatialPolygonDataFrame

Tôi đang sử dụng thư viện R maptools để phân tích shapefile thành một danh sách các đa giác. Hàm readShapeSpatial

shp <- readShapeSpatial("<path to my shapefile>")

cung cấp cho tôi một đối tượng SpatailPolygonDataFrame . Trong ví dụ của tôi, SpatialPolygonsDataFrame của tôi có các cột sau:

> names(shp)
[1] "AREA"       "PERIMETER"  "COMAREA_"   "COMAREA_ID" "AREA_NUMBE"
[6] "COMMUNITY"  "AREA_NUM_1" "SHAPE_AREA" "SHAPE_LEN

Tôi biết rằng tôi có thể loại bỏ một số đa giác bằng cách của họ row.id, ví dụ

shp.dropI  <- shp[-i, ]
shp.subset <- shp[i %in% c(1,2,3),]

Bây giờ tôi muốn thả (các) đa giác nhất định với một điều kiện, nói AREA > 10. Làm thế nào để tôi thực hiện điều này một cách thanh lịch? Phương pháp duy nhất tôi có bây giờ là lặp qua tất cả các hàng và tìm tương ứng row.id.

Đầu tiên, tôi rất khuyên bạn nên sử dụng readOGR, từ thư viện rgdal, để đọc shapefile của bạn. Nó sẽ giữ lại thông tin chiếu (proj4 chuỗi) và lưu lại nhiều vấn đề đau đầu, khi khớp chuỗi, sử dụng các chức năng khác.

Hai cách nhanh chóng để thực hiện những gì bạn đang có sau khi sử dụng một chỉ mục hoặc sử dụng tập hợp con. Điều này sẽ giữ lại các đa giác có diện tích <10 (giảm những người> 10).

shp.sub <- shp[shp$AREA < 10,] 
shp.sub <- subset(shp, AREA < 10)