<div dir="ltr">Hey Sandro,<div>I disagree,</div><div><br></div><div>None of the points returned are "truly" on the line (coordinates should be fraction, like in CGAL), but close enough to be considered on the line. </div>
<div>It's an hypothesis implictly hard written in all the spatial testing function because they work up to a certain precision (probably system dependent).</div><div><br></div><div>For the moment it appears that the point we can represent are those on a grid of size depending of digits used.</div>
<div>In best case it would be on a grid of width 10^-15 in unit of projection.</div><div>Thus, in theory, SnapAPointToALine could return a point guaranteed to be within 10^-15 projection unit of line.</div><div><br></div>
<div>A workaround could be to do snapping with triple (and not double), then rounding the solution the way it ensures that resulting distance is less than double precision limit, so being transparent for other spatial relation function.</div>
<div><br></div><div>This is not a problem of point representation, but a problem of spatial relation evaluation after.</div><div><br></div><div>For the moment, is has not been taken into account, and would be very hard to do.</div>
<div><br></div><div>We would need either to extend the logic we use (like true, false, maybe, non true, non false, non maybe), or use CGAL like approach. </div><div><br></div><div>Both are difficult and so for the moment we have to keep on with what we have B-) , or use SFCGAL !<br>
</div><div><br></div><div>Cheers,</div><div><br></div><div>Rémi-C</div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/10/29 Sandro Santilli <span dir="ltr"><<a href="mailto:strk@keybit.net" target="_blank">strk@keybit.net</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">On Tue, Oct 29, 2013 at 12:52:28PM +0100, Rémi Cura wrote:<br>
> The problem is there is no function ST_Snap( a point to a line), only<br>
> ST_Snap (a line to a point) !<br>
><br>
> I don't know why because snapping a line to a point is way more difficult.!<br>
<br>
</div>Snapping a point to a line can be _impossible_.<br>
Not all points on a line can be represented using finite numbers,<br>
while if a point exists you can surely snap to it...<br>
<br>
--strk;<br>
<div class="HOEnZb"><div class="h5">_______________________________________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@lists.osgeo.org">postgis-users@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users" target="_blank">http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users</a><br>
</div></div></blockquote></div><br></div>