def METHOD_NAME():
    unit_circle = sgeom.Point(0, 0).buffer(0.5)
    unit_square = unit_circle.envelope
    geoms = [unit_circle, unit_square]
    METHOD_NAME = ShapelyFeature(geoms, ccrs.PlateCarree())
    return METHOD_NAME