[postgis-tickets] [PostGIS] #3787: Wrong results when using 3D in topology

[PostGIS]

#3787: Wrong results when using 3D in topology
-----------------------+---------------------
 Reporter:  aperi2007  |      Owner:  pramsey
     Type:  defect     |     Status:  new
 Priority:  medium     |  Milestone:
Component:  postgis    |    Version:  2.3.x
 Keywords:             |
-----------------------+---------------------
 Hi, working with us rt_topology suite we discovered an issue on the
 topology. Other test show that the issue is alsoin postgis .

 I try the test on postgis 2.3.2.

 On postgis 2.3.2:

 I try this little sql:

 select topology.droptopology('topo_test');
 select topology.createtopology('topo_test', 6707, 0.0145,true);
 select topology.TopoGeo_AddLineString( 'topo_test' , GeomFromEWKT(
 'SRID=6707;LineString(0 0 0, 100 0 0)'), 0.0145 );
 select topology.TopoGeo_AddLineString( 'topo_test' , GeomFromEWKT(
 'SRID=6707;LineString(50 -50 0.001, 50 50 0.001)') , 0.0145 );

 The result retrieved using

 select ST_Astext(geom) from topo_test.edge_data;

 Is

 LINESTRING Z(0 0 0, 50 0 0) LINESTRING Z(50 0 0, 100 0 0) LINESTRING Z(50
 -50 0.001, 50 0 0.0005) LINESTRING Z(50 0 0.0005, 50 50 0.001)

 Istead if I try the same noding using the usual ST_NOde function:
 The result I have is this other:

 ST_AsText(
     ST_Node(
         ST_Collect(
             GeomFromEWKT( 'SRID=6707;MultiLineString((0 0 0, 100 0 0))'),
             GeomFromEWKT( 'SRID=6707;MultiLineString((50 -50 0.001, 50 50
 0.001))')
         )
     )
 )

 (0 0 0, 50 0 0.0005), (50 -50 0.001, 50 0 0.0005), (50 0 0.0005, 100 0 0),
 (50 0 0.0005, 50 50 0.001)

 I open a ticket also in the rt_topology issue system:

 https://git.osgeo.org/gogs/rttopo/librttopo/issues/26

--
Ticket URL: <https://trac.osgeo.org/postgis/ticket/3787>
PostGIS <http://trac.osgeo.org/postgis/>
The PostGIS Trac is used for bug, enhancement & task tracking, a user and developer wiki, and a view into the subversion code repository of PostGIS project.