Bắt các khu vực đa giác bằng cách sử dụng geopandas?

Đưa ra geopandas GeoDataFramemột loạt các đa giác, tôi muốn có được diện tích tính bằng km vuông của mỗi tính năng trong danh sách của tôi.

Đây là một vấn đề khá phổ biến và giải pháp được đề xuất thông thường trong quá khứ là sử dụng shapelyvà pyprojtrực tiếp (ví dụ ở đây và đây ).

Có cách nào để làm điều này trong sạch geopandas?

Nếu các crs của GeoDataFrame được biết đến (EPSG: 4326 đơn vị = độ, ở đây), bạn không cần Shapely, cũng không phải pyproj trong tập lệnh của mình vì GeoPandas sử dụng chúng).

import geopandas as gpd
test = gpd.read_file("test_wgs84.shp")
print test.crs
test.head(2)

Bây giờ hãy sao chép GeoDataFrame của bạn và thay đổi phép chiếu sang hệ thống Cartesian (EPSG: 3857, unit = m như trong câu trả lời của ResMar)

tost = test.copy()
tost= tost.to_crs({'init': 'epsg:3857'})
print tost.crs
tost.head(2)

Bây giờ diện tích bằng km vuông

tost["area"] = tost['geometry'].area/ 10**6
tost.head(2)

Nhưng các bề mặt trong phép chiếu Mercator là không chính xác, do đó, với các phép chiếu khác tính bằng mét.

tost= tost.to_crs({'init': 'epsg:32633'})
tost["area"] = tost['geometry'].area/ 10**6
tost.head(2)

Tôi tin là có. Sau đây nên làm việc:

gdf['geometry'].to_crs({'init': 'epsg:3395'})\
               .map(lambda p: p.area / 10**6)

Điều này chuyển đổi hình học thành hình chiếu có diện tích bằng nhau, tìm nạp shapelykhu vực (được trả về trong m ^ 2) và ánh xạ tới km ^ 2 (bước cuối cùng này là tùy chọn).