> I don't understand all of your caveats here, but I assume
> PointInPolygon(box2d(<polygon>), the_geom doesn't cut it.

If you are determining whether 100,000 points are within a polygon, you can 
save yourself a ton of work by "preparing" a suitable data structure from your 
polygon for use in the point-in-poly test. If you build the right structure, 
each PiP test takes almost no time at all.

Chris
 





