[PostGIS] #5713: GEOS3.12 - "TopologyException: side location conflict" on valid input geometries
PostGIS
trac at osgeo.org
Thu Apr 18 15:59:58 PDT 2024
#5713: GEOS3.12 - "TopologyException: side location conflict" on valid input
geometries
-------------------------+---------------------------
Reporter: dannytoone | Owner: pramsey
Type: defect | Status: closed
Priority: medium | Milestone: PostGIS 3.4.3
Component: postgis | Version: 3.4.x
Resolution: wontfix | Keywords:
-------------------------+---------------------------
Changes (by pramsey):
* status: new => closed
* resolution: => wontfix
Old description:
> This is from a recent version of PostGIS running on AWS RDS:
> {{{
> POSTGIS="3.4.0 0874ea3" [EXTENSION] PGSQL="160" GEOS="3.12.0-CAPI-1.18.0"
> PROJ="8.0.1 NETWORK_ENABLED=OFF URL_ENDPOINT=https://cdn.proj.org
> USER_WRITABLE_DIRECTORY=/tmp/proj
> DATABASE_PATH=/rdsdbbin/postgres-16.1.R1/share/proj/proj.db" GDAL="GDAL
> 3.4.3, released 2022/04/22" LIBXML="2.9.1" LIBJSON="0.15"
> LIBPROTOBUF="1.3.2" WAGYU="0.5.0 (Internal)" RASTER
> }}}
>
> ST_Equals as well as ST_Touches fail on these input geometries. I will
> attach a csv output containing two geometry columns which fail.
New description:
This is from a recent version of PostGIS running on AWS RDS:
{{{
POSTGIS="3.4.0 0874ea3" [EXTENSION] PGSQL="160" GEOS="3.12.0-CAPI-1.18.0"
PROJ="8.0.1 NETWORK_ENABLED=OFF URL_ENDPOINT=https://cdn.proj.org
USER_WRITABLE_DIRECTORY=/tmp/proj
DATABASE_PATH=/rdsdbbin/postgres-16.1.R1/share/proj/proj.db" GDAL="GDAL
3.4.3, released 2022/04/22" LIBXML="2.9.1" LIBJSON="0.15"
LIBPROTOBUF="1.3.2" WAGYU="0.5.0 (Internal)" RASTER
}}}
ST_Equals as well as ST_Touches fail on these input geometries. I will
attach a csv output containing two geometry columns which fail.
--
Comment:
Confirmed. Since the Relate code is likely to be removed and replaced from
GEOS in the next cycle, this particular issue will likely not be fixed. A
work around it to apply a very small precision reduction, which should
shake out the very small nearly parallel lines in your original that are
causing the problem.
{{{
select st_equals(ST_ReducePrecision(mark_geom,0.00000001),
st_reduceprecision(fips_geom, 0.00000001)) from g5713
}}}
--
Ticket URL: <https://trac.osgeo.org/postgis/ticket/5713#comment:1>
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