Nest center should be inside the polygon even if it is concave about nestwatcher HOT 5 CLOSED

#pip install python-polylabel
#pip install shapely
#pip install geojson

from geojson import (
    Feature,
    FeatureCollection,
    dumps
)

from polylabel import polylabel
from shapely.ops import polylabel as polylabel2
from shapely.geometry import shape
import json

NEST_NAMES = ['Foobar Park', 'Foobar2 Park', 'Foobar3 Park', 'Foobar4']
NEST_IDS = [605XXXX44]

with open("nests.json") as json_file:
    json_data = json.load(json_file)

for item in json_data['features']:
    if item['properties']['name'] in NEST_NAMES:
        print(item['properties']['name'])
        print(polylabel(item['geometry']['coordinates']))
        print(polylabel2(shape(item['geometry']), tolerance=10))
    if item['id'] in NEST_IDS:
        print(polylabel(item['geometry']['coordinates']))

from geojson import (
    Feature,
    FeatureCollection,
    dumps
)

from polylabel import polylabel
from shapely.ops import polylabel as polylabel2
from shapely.geometry import shape
import json

with open("nests.json") as json_file:
    json_data = json.load(json_file)

for item in json_data['features']:
    print(item['properties']['name'])
    print(polylabel(item['geometry']['coordinates'], 1e-6))
    point = polylabel2(shape(item['geometry']), tolerance=1e-6)
    print(point)
    assert shape(item['geometry']).contains(point)
    print()