Hi all,<br><br>needed to create linestrings between a point relation and a line relation. <br>used makeline, line_interpolate_point and line_locate_point with good results in the following query:<br>_______________________________________________________________________________________
<br>create table tmp.jumper6 with oids as <br>select distinct on (<a href="http://b.id">b.id</a>) <a href="http://b.id">b.id</a>,distance(a.geo,b.geo),<br>    makeline(line_interpolate_point(a.geo, line_locate_point(a.geo
, b.geo)),b.geo)<br>from wkh a<br>inner join xfmr_mapping.xfmr_deadend_ b on expand(a.geo,15) && b.geo AND distance(a.geo,b.geo) < 20<br>where b.category = 'TRANSFORMER'<br>order by <a href="http://b.id">
b.id</a>, distance<br>________________________________________________________________________________________<br><br>the new table has 387 rows. <br><br>When the following query is run:<br> ________________________________________________________________________________________
<br>select distance(a.makeline,b.geo),intersects(a.makeline,b.geo),a.*,b.gid from tmp.jumper6 a<br>inner join wkh b on a.makeline && expand(b.geo,10) AND intersects(a.makeline,b.geo) <br>order by id<br>________________________________________________________________________________________
<br><br>94 rows are returned, if I change "intersects(a.makeline,b.geo)" to "distance(a.makeline,b.geo) = 0", 240 rows are returned.<br><br>So i guess the question is this: the distance between two geometries can be 0 and they do not intersect? Is this a matter of precision? or is this a bug?
<br><br>Regards,<br><br>Rhys<br>