Наведені вище відповіді здаються правильними, ПРИКЛАДНО, що в якийсь момент нещодавно параметри 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