<div dir="ltr"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Other ideas about performance improvements are welcome.  For instance, I have had some ideas about only noding geometry in the region of actual overlap in the case of an intersection computation.  That probably has implications for robustness, however.  Needs some research...</div></div></blockquote><div><br><br>I tried to implement this speed up in PostGIS via box clipping.</div><div>Current GEOS recangle clipping is unfortunately not robust, but helps quite a lot when not crashing.<br><br>I tried redoing clipping by doing just O(N) clamping to the box. Result retains many required properties (ie, raycasting will give correct inside/outside result for it, area is okay) but has holes touching edges and edges touching themselves. GEOSIntersection can't handle that. <br><br>Is there a simple way to recover the boundary for such geometry to feed further to GEOSIntersection? Buffer(,0) does not help in many cases, ST_MakeValid too.<br><br>Can we make GEOSIntersection handle such invalid input, or is it completely incompatible with monotone chain conception?</div></div></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Darafei Praliaskouski<br>Support me: <a href="http://patreon.com/komzpa">http://patreon.com/komzpa</a></div></div>