[postgis-users] line_interpolate_point does not return a point that intersects the line

2010-05-03 Thread Francis Markham
Hi all, I think I've got a bit of a problem with snapping points to lines. I had thought that using line_interpolate_point(..., line_locate_point(...)) would do the trick, but this seems not to be the case: gis=# SELECT ST_AsText(v.the_geom), ST_SRID(v.the_geom), ST_AsText(r.the_geom),

Re: [postgis-users] line_interpolate_point does not return a point that intersects the line

2010-05-03 Thread strk
On Mon, May 03, 2010 at 04:41:09PM +1000, Francis Markham wrote: Hi all, I think I've got a bit of a problem with snapping points to lines. I had thought that using line_interpolate_point(..., line_locate_point(...)) would do the trick, but this seems not to be the case: Francis, thing is

Re: [postgis-users] line_interpolate_point does not return a point that intersects the line

2010-05-03 Thread Francis Markham
Mmm, so it is a floating point issue - the distance is 9.57711783686133e-12 meters. Is there anyway to specify a tolerance for ST_Intersects and similar predicates? -Francis On 3 May 2010 17:53, strk s...@keybit.net wrote: On Mon, May 03, 2010 at 04:41:09PM +1000, Francis Markham wrote: Hi

Re: [postgis-users] line_interpolate_point does not return a point that intersects the line

2010-05-03 Thread strk
On Mon, May 03, 2010 at 06:15:14PM +1000, Francis Markham wrote: Mmm, so it is a floating point issue - the distance is 9.57711783686133e-12 meters. Is there anyway to specify a tolerance for ST_Intersects and similar predicates? You may use DWithin() for a tolerance-based intersection

Re: [postgis-users] line_interpolate_point does not return a point that intersects the line

2010-05-03 Thread Francis Markham
D'oh. Totally forgot about DWithin! Problem, solved. Thanks very much! -Francis On 3 May 2010 18:24, strk s...@keybit.net wrote: On Mon, May 03, 2010 at 06:15:14PM +1000, Francis Markham wrote: Mmm, so it is a floating point issue - the distance is 9.57711783686133e-12 meters. Is