code: def area_from_lon_lat_poly(geometry):
      
  
      import pyproj
      from shapely.ops import transform
      from functools import partial
  
  
      project = partial(
          pyproj.transform,
          pyproj.Proj(init='epsg:4326'), # Source: Lon-Lat
          pyproj.Proj(proj='aea')) # Target: Albers Equal Area Conical https://en.wikipedia.org/wiki/Albers_projection
  
      new_geometry = transform(project, geometry)
  
      #default area is in m^2
      return new_geometry.area/1e6
docs: Compute the area in km^2 of a shapely geometry, whose points are in longitude and latitude.
  
      Parameters
      ----------
      geometry: shapely geometry
          Points must be in longitude and latitude.
  
      Returns
      -------
      area:  float
          Area in km^2.