Pozyskiwanie obszarów wielokątów za pomocą geopandaz?

Biorąc pod uwagę geopandas GeoDataFrameserię wielokątów, chciałbym uzyskać powierzchnię w km2 każdej funkcji na mojej liście.

Jest to dość powszechny problem, a zwykle proponowane rozwiązanie w przeszłości był w użyciu shapelyi pyprojbezpośrednio (np tutaj i tutaj ).

Czy istnieje sposób, aby to zrobić w czystej postaci geopandas?

Jeśli znany jest crs GeoDataFrame (EPSG: 4326 jednostka = stopień, tutaj), nie potrzebujesz Shapely ani pyproj w swoim skrypcie, ponieważ używa go GeoPandas).

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

Teraz skopiuj GeoDataFrame i zmień rzut na system kartezjański (EPSG: 3857, jednostka = m jak w odpowiedzi ResMar)

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

Teraz obszar w kilometrach kwadratowych

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

Ale powierzchnie w rzucie Mercatora są niepoprawne, więc w przypadku innych rzutów w metrach.

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

Wierzę tak. Następujące powinny działać:

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

Konwertuje geometrię na rzut o jednakowym obszarze, pobiera shapelyobszar (zwrócony wm ^ 2) i odwzorowuje go na km ^ 2 (ten ostatni krok jest opcjonalny).