%pylab --no-import-all inline from geo.generator import randomConvexPolygon, randomConcaveTiling from geo.drawer import plot, plotPoints from kirkpatrick import Locator # Create a random region initial = randomConvexPolygon(100, k=100) polygons = randomConcaveTiling(initial) plot(polygons) # Run Kirkpatrick's locator = Locator(polygons) # Visualize triangulated boundary plot(locator.boundary, style='g--') plot(locator.regions, style='g-') # Locate a point target = initial.smartInteriorPoint() region = locator.locate(target) plotPoints(target, style='bo') plot(locator.regions, style='g-') plot(region, style='ro-')