<div dir="ltr">Hey Nicklas,<div>thanks for your answer.<br><br>I'm simply trying to solve a problem in current ST_Intersects behavior regarding a point and a line. IN most of the case, a point on a line(to double limit precision) is not on the line for st_intersects !</div>
<div><br></div><div>it works as intended in very very few cases (ie, all the points that intersects a line is a very sparse set and is not a line at all !).</div><div><br></div><div>Here is a simple illustration : <a href="http://hpics.li/c3485b5">http://hpics.li/c3485b5</a> represents a line with max zoom : it is all the point that follow the line equation and are rounded to the(16 or 17 digits)</div>
<div><a href="http://hpics.li/f30aeb7">http://hpics.li/f30aeb7</a> now the big points represents the point that are on the line for ST_Intersect, in a favorable case (bad case could mean no points at all).<br></div><div><br>
</div><div>This is without speaking of custom precision model, ie working with the postgis 'infinite' precision (double precision, 15 to 17 digits)</div><div><br></div><div>I tried several workaround (like several random walk along the line to find a point that st_intersects like, </div>
<div>and also reverse engineer to find which points are seen on the line by st_intersects (related to prime factorization of slope it appears) ). </div><div><br></div><div>The issue is there is a flaw in ST_Intersects and in bad cases there is no workaround.</div>
<div><br></div><div>Could you elaborate about 'native' functions please?</div><div><br></div><div>Cheers,</div><div><br>Rémi-C</div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/11/6 Nicklas Avén <span dir="ltr"><<a href="mailto:nicklas.aven@jordogskog.no" target="_blank">nicklas.aven@jordogskog.no</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div>I think postgis has native functions for st_intersects when at least one of the geometries is a point. And there is two paths in the code. One that prepares the "not point" geometry if it is to be used multiple times and one that doesn't.</div>
<div><br></div><div>I do not get the idea do far. Is it a tolerance you are implementing? </div><div><br></div><div><br></div><div><br></div><div>Best regards</div><div>Nicklad</div><div><br></div><div><div style="font-size:75%;color:#575757">
Skickat från min Samsung Mobil</div></div><br><br><br>-------- Originalmeddelande --------<br>Från: Rémi Cura <<a href="mailto:remi.cura@gmail.com" target="_blank">remi.cura@gmail.com</a>> <br>Datum: 06-11-2013  8:59  (GMT+01:00) <br>
Till: PostGIS Development Discussion <<a href="mailto:postgis-devel@lists.osgeo.org" target="_blank">postgis-devel@lists.osgeo.org</a>> <br>Rubrik: [postgis-devel] SNapPointToLine : GEOS : where is the code for st_intersects between line and point? <br>
<div class="im"> <br><br><div dir="ltr">Hey,<div>I've been working on a ST_SnapPointToLine for the past week,</div><div>by reverse engineering when points are considered to be on line by st_intersect.</div><div><br></div>
<div>The (classic) problem is purely a numerical issue, and I suspect a problem in design of st_intersects.</div>
<div><br></div><div><br></div><div>So far my solution doesn't work for all line, can someone please point me in the right direction in GEOS so I can see if something can be improved regarding precision.</div><div><br>

</div><div>I know where the API is, but I'm looking for the exact part that decide if a point and a line intersects.</div><div><br></div><div>Cheers,</div><div><br></div><div>Rémi-C</div></div>
</div></div><br>_______________________________________________<br>
postgis-devel mailing list<br>
<a href="mailto:postgis-devel@lists.osgeo.org">postgis-devel@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-devel" target="_blank">http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-devel</a><br></blockquote></div><br></div>