Mendapatkan area poligon menggunakan geopanda?

Diberikan geopandas GeoDataFrameserangkaian poligon, saya ingin mendapatkan luas dalam km persegi dari setiap fitur dalam daftar saya.

Ini adalah masalah yang cukup umum, dan solusi yang disarankan di masa lalu adalah menggunakan shapelydan pyprojsecara langsung (misalnya di sini dan di sini ).

Apakah ada cara untuk melakukan ini dengan murni geopandas?

Jika crs GeoDataFrame diketahui (EPSG: 4326 unit = degree, di sini), Anda tidak perlu Shapely, atau pyproj dalam skrip Anda karena GeoPandas menggunakannya).

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

Sekarang salin GeoDataFrame Anda dan ubah proyeksi ke sistem Cartesian (EPSG: 3857, unit = m seperti pada jawaban ResMar)

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

Sekarang area dalam kilometer persegi

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

Tetapi permukaan dalam proyeksi Mercator tidak benar, demikian juga dengan proyeksi lainnya dalam meter.

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

Saya percaya ya. Berikut ini harus bekerja:

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

Ini mengubah geometri menjadi proyeksi dengan luas yang sama, mengambil shapelyarea tersebut (dikembalikan dalam m ^ 2), dan memetakannya menjadi km ^ 2 (langkah terakhir ini opsional).