Polygonflächen mit Geopandas erhalten?

Angesichts einer geopandas GeoDataFrameReihe von Polygonen möchte ich die Fläche der einzelnen Elemente in meiner Liste in km² anzeigen.

Dies ist ein ziemlich häufiges Problem, und die in der Vergangenheit übliche Lösung bestand darin, shapelyund pyprojdirekt zu verwenden (z. B. hier und hier ).

Gibt es eine Möglichkeit, dies in rein zu tun geopandas?

Wenn die crs des GeoDataFrame bekannt sind (EPSG: 4326 unit = degree, hier), benötigen Sie weder Shapely noch Pyproj in Ihrem Skript, da GeoPandas sie verwendet.

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

Kopieren Sie nun Ihren GeoDataFrame und ändern Sie die Projektion in ein kartesisches System (EPSG: 3857, Einheit = m wie in der Antwort von ResMar)

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

Jetzt ist die Fläche in Quadratkilometern

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

Die Flächen in der Mercator-Projektion sind jedoch nicht korrekt, also bei anderen Projektionen in Metern.

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

Ich glaube ja Folgendes sollte funktionieren:

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

Dadurch wird die Geometrie in eine flächengleiche Projektion konvertiert, die Fläche abgerufen shapely(zurückgegeben in m ^ 2) und auf km ^ 2 abgebildet (dieser letzte Schritt ist optional).