docs: def get_poly_area_geo(poly):
      """
      Calculates the area in meters squared of the individual polygon
      """
      minx, miny, maxx, maxy = poly.bounds
      # reproject polygon to get area
      reprojected_for_area = Proj("+proj=aea +lat_1={0} +lat_1={1} "
                                  "+lat_0={2} +lon_0={3}"
                                  .format(miny,
                                          maxy,
                                          (miny + maxy) / 2.0,
                                          (minx + maxx) / 2.0))
      geographic_proj = Proj(init='epsg:4326')
      project_func = partial(transform,
                             geographic_proj,
                             reprojected_for_area)
      reprojected_poly = shapely_transform(project_func, poly)
      return reprojected_poly.area