docs: def contains(polygon, point):
          """
          Tests whether point lies within the polygon
          """
          in_hole = functools.reduce(
              lambda P, Q: P and Q,
              [interior.covers(point) for interior in polygon.interiors]
          ) if polygon.interiors else False
          return polygon.covers(point) and not in_hole