def METHOD_NAME(geometries: Iterator[MultiPolygon]) -> Optional[MultiPolygon]:
    total_area = GEOSGeometry("POINT EMPTY", srid=4326)
    for geom in geometries:
        total_area = total_area.union(geom)
    if isinstance(total_area, Point):
        return None
    return to_multipolygon(total_area)