> > The other place in the overlay algorithm where point-in-polygon is used is > during the polygon-extraction phase [1]. This is to assign the holes in > the overlay result to their parent shell. In this particular test case it > might be time-consuming, since all 4262 holes are tested against the > original 130,327 vertex parent ring, using a linear scan. This could be > improved by using an index. It will need a new class or a new entry point > to IndexedPointInAreaLocator however, since that class currently only > accepts Geometry inputs as the target. >
@Paul - do you think this is another hotspot? > > [1] https://github.com/libgeos/geos/blob/8bc1ea0f8019bcb83c60126e7883ae428ea7c677/src/operation/overlay/PolygonBuilder.cpp#L352
_______________________________________________ geos-devel mailing list geos-devel@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/geos-devel