def find_tiles_inside(domain, corners):
    'Determine which tiles are inside `domain`\n\n    The function uses `corners` the list of corners for each tile\n    '
    p = Polygon(domain)
    tileslist = []
    for (tile, c) in corners.items():
        q = Polygon(c)
        if (p.overlaps(q) or p.contains(q)):
            tileslist += [tile]
    return tileslist
Determine which tiles are inside `domain`

The function uses `corners` the list of corners for each tile
def find_tiles_inside(domain, corners):
    'Determine which tiles are inside `domain`\n\n    The function uses `corners` the list of corners for each tile\n    '
    p = Polygon(domain)
    tileslist = []
    for (tile, c) in corners.items():
        q = Polygon(c)
        if (p.overlaps(q) or p.contains(q)):
            tileslist += [tile]
    return tileslist