Jawaban di atas tampaknya benar, KECUALI bahwa pada beberapa titik baru-baru ini, parameter lat1 dan lat2 dalam kode pyproj diganti nama dengan garis bawah: lat_1 dan lat_2 (lihat /programming//a/55259718/1538758 ). Saya tidak punya cukup perwakilan untuk berkomentar, jadi saya membuat jawaban baru (maaf bukan maaf)
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