2014-04-12 16 views
5

Làm cách nào để xác định được khu vực của shapefile?Xác định khu vực của Shapefiles

Khi tôi sử dụng summary(data), nó cung cấp cho khu vực, nhưng sau đó tôi không biết khu vực của mỗi đa giác shapefile là gì.

Ngoài ra, khu vực có thể được chuyển đổi thành mét không?

+0

nhìn trong 'rgeos' thư viện cho hàm 'gArea()' –

+0

Phụ thuộc vào hệ tọa độ. Cho chúng tôi biết 'tóm tắt (dữ liệu) 'thực sự mang lại và có thể chúng tôi có thể giúp nhiều hơn. – Spacedman

+0

Tôi đang làm việc với dữ liệu DIVA-GIS. Và không thể tìm ra các đơn vị của khu vực. Nếu thông tin nằm trong vĩ độ dài nếu tôi sử dụng hàm 'gArea()', làm cách nào để biết các đơn vị của vùng này? – GabrielMontenegro

Trả lời

8
# load the rgeos library 
library(rgeos) 

# make a polygon (borrowed from ref manual for package) 
sample_poly <- readWKT("POLYGON((0 1,0.95 0.31,0.59 -0.81,-0.59 -0.81,-0.95 0.31,0 1))") 

# and calculate the area 
gArea(sample_poly) 
[1] 2.3803 
+0

Từ một shapefile của các quốc gia khác nhau ví dụ từ DIVA-GIS, làm thế nào bạn có thể tìm ra các đơn vị của khu vực đó? – GabrielMontenegro

+0

@ JavierM88 Tôi sẽ kiểm tra tài liệu liên quan đến tập dữ liệu của bạn. Không biết vị trí mà các đơn vị được nhúng vào chính bản thân shapefile. –

2
[email protected][[1]]@Polygons[[1]]@area 
[email protected][[2]]@Polygons[[1]]@area 

sẽ cung cấp cho bạn các lĩnh vực của đa giác đầu tiên và thứ hai, bạn có thể figue một vòng lặp for để cung cấp cho bạn tất cả hoặc chỉ cho chúng ta những cái tên với:

names(data) 
+1

'sapply' có thể thuận tiện hơn' for': 'sapply (dữ liệu @ đa giác, hàm (x) x @ Đa giác [[1]] @ area)' – jbaums

+1

Phương pháp này nguy hiểm trong trường hợp bạn có đa giác có lỗ, như nó sẽ không xem xét lỗ và do đó ước lượng quá mức diện tích. – Matifou

Các vấn đề liên quan