[PostGIS] #5786: Side-location conflict: new edge starts in face 42 and ends in face 10

PostGIS trac at osgeo.org
Thu Oct 3 07:39:29 PDT 2024


#5786: Side-location conflict: new edge starts in face 42 and ends in face 10
--------------------------------+---------------------------
  Reporter:  Lars Aksel Opsahl  |      Owner:  strk
      Type:  defect             |     Status:  new
  Priority:  medium             |  Milestone:  PostGIS 3.5.1
 Component:  topology           |    Version:  master
Resolution:                     |   Keywords:  robustness
--------------------------------+---------------------------
Comment (by strk):

 Edge 2 is found to be ABOVE edge 1 before the snapping:
 {{{
 =# select edge_id, st_azimuth(st_pointn(geom,1),st_pointn(geom,2)) from
 t5786.edge where edge_id = 1;
  edge_id |    st_azimuth
 ---------+-------------------
        1 | 4.688117807826012

 =# select edge_id, st_azimuth(st_pointn(geom,3),st_pointn(geom,2)) from
 t5786.edge where edge_id = 2;
  edge_id |    st_azimuth
 ---------+-------------------
        2 | 4.688117807826215
 }}}

 Once both edges are snapped to the new node, edge 2 is found to be BELOW
 edge 1:
 {{{
 =# update t5786.edge set geom = ST_Snap(geom, 'POINT(11.812029186127067
 59.7793864213727)', 1e-6);
 UPDATE 2
 =# select edge_id, st_azimuth(st_pointn(geom,2),st_pointn(geom,3)) from
 t5786.edge where edge_id = 1;
  edge_id |    st_azimuth
 ---------+-------------------
        1 | 4.688117807809833
 (1 row)

 strk at strk=# select edge_id,
 st_azimuth(st_pointn(geom,3),st_pointn(geom,2)) from t5786.edge where
 edge_id = 2;
  edge_id |    st_azimuth
 ---------+-------------------
        2 | 4.688117807781473
 (1 row)
 }}}
-- 
Ticket URL: <https://trac.osgeo.org/postgis/ticket/5786#comment:16>
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.


More information about the postgis-tickets mailing list