Приведенные выше ответы кажутся правильными, за исключением того, что в последнее время параметры lat1 и lat2 в коде pyproj были переименованы с подчеркиванием: lat_1 и lat_2 (см. Https://stackoverflow.com/a/55259718/1538758 ). У меня недостаточно представителей, чтобы комментировать, поэтому я делаю новый ответ (извините, не извините)
import pyproj
import shapely
import shapely.ops as ops
from shapely.geometry.polygon import Polygon
from functools import partial
geom = Polygon([(0, 0), (0, 10), (10, 10), (10, 0), (0, 0)])
geom_area = ops.transform(
partial(
pyproj.transform,
pyproj.Proj(init='EPSG:4326'),
pyproj.Proj(
proj='aea',
lat_1=geom.bounds[1],
lat_2=geom.bounds[3])),
geom)
# Print the area in m^2
print geom_area.area