[postgis-users] St_makevalid issue

[Imre Samu]

> Doesn't anybody have an idea? Or is at least curious?
>"POSTGIS="2.5.1 r17027" [EXTENSION] PGSQL="110" GEOS="3.7.0-CAPI-1.11.0
3.7.1"

imho: try to update your environment with the latest patches/fixes  and
please re-check.
Postgis:  from 2.5.1  -> to 2.5.4 (
https://postgis.net/2020/02/28/postgis-2.5.4/
https://git.osgeo.org/gitea/postgis/postgis/raw/tag/2.5.4/NEWS )
Geos:  from 3.7.0    ->  to the latest (3.8.1 )  (
https://github.com/libgeos/geos/blob/master/NEWS )

Imre

Paul van der Linden <paul.doskabouter at gmail.com> ezt írta (időpont: 2020.
márc. 26., Cs, 15:54):

> Doesn't anybody have an idea? Or is at least curious?
>
> Especially the select st_isvalid(geo) where not st_valid returning a
> "true" row?
>
> >I have a query that gives inconsistent results.
> >It almost seems like st_makevalid doesn't always produce the same results
> >given the same input, and sometimes even returns invalid geometries...
> >
> >Results are
> >f;1159320407;"Bangladesh"
> >f;1159321073;"Mozambique"
> >f;1159316737;"Nord"
> >f;1159314827;"Uusimaa"
> >
> >second run:
> >f;1159321073;"Mozambique"
> >t;1159315455;"Syddanmark"
> >t;1159314827;"Uusimaa"
> >t;1159317795;"Gerona"
> >
> >third run:
> >f;1159321073;"Mozambique"
> >t;1159315455;"Syddanmark"
> >
> >note the 't' in the second and third run...
> >
> >Query:
> >
> >SELECT ST_isvalid(simplifiedgeo2),* FROM (
> >SELECT
> >    foo2.ne_id,
> >    foo2.name,
> >    CASE
> >        WHEN ST_isvalid(foo2.simplifiedgeo) THEN foo2.simplifiedgeo
> >        ELSE ST_makevalid(foo2.simplifiedgeo)
> >    END AS simplifiedgeo2
> >   FROM (
> >     SELECT
> >       ST_simplifypreservetopology(foo.geom,
> >         LEAST(
> >           ST_xmax(foo.shiftedgeo::box3d) -
> ST_xmin(foo.shiftedgeo::box3d),
> >ST_ymax(foo.shiftedgeo::box3d) - ST_ymin(foo.shiftedgeo::box3d))
> >         ) AS simplifiedgeo,
> >       foo.ne_id,
> >       foo.name,
> >       foo.geom
> >     FROM (
> >        SELECT
> >          natearth.ne_id,
> >          natearth.geom,
> >          natearth.name,
> >          CASE
> >              WHEN (
> >                ST_xmax(ST_shiftlongitude(natearth.geom)::box3d) -
> >ST_xmin(ST_shiftlongitude(natearth.geom)::box3d))
> >                <
> >                (ST_xmax(natearth.geom::box3d) -
> >ST_xmin(natearth.geom::box3d)
> >              ) THEN ST_shiftlongitude(natearth.geom)
> >              ELSE natearth.geom
> >          END AS shiftedgeo
> >          FROM (
> >             SELECT
> >                ne_10m_admin_0_countries_lakes.ne_id,
> >                ne_10m_admin_0_countries_lakes.geom,
> >                ne_10m_admin_0_countries_lakes.name
> >              FROM natural_earth.ne_10m_admin_0_countries_lakes
> >          UNION ALL
> >             SELECT
> >                ne_10m_admin_1_states_provinces_lakes.ne_id,
> >                ne_10m_admin_1_states_provinces_lakes.geom,
> >                ne_10m_admin_1_states_provinces_lakes.name
> >             FROM natural_earth.ne_10m_admin_1_states_provinces_lakes
> >          ) natearth
> >     ) foo
> >   ) foo2
> >) as fff
> >WHERE NOT ST_isvalid(simplifiedgeo2)
> >
> >Any advice on how to further debug this?
> >
> >"POSTGIS="2.5.1 r17027" [EXTENSION] PGSQL="110" GEOS="3.7.0-CAPI-1.11.0
> >3.7.1" PROJ="Rel. 4.9.3, 15 August 2016" GDAL="GDAL 2.2.4, released
> >2018/03/19 GDAL_DATA not found" LIBXML="2.7.8" LIBJSON="0.12"
> >LIBPROTOBUF="1.2.1" (core procs from "2.5.0 r16836" nee (...)"
> >"PostgreSQL 11.2, compiled by Visual C++ build 1914, 64-bit"
> >
> >Running on windows 64-bit.
> _______________________________________________
> postgis-users mailing list
> postgis-users at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/postgis-users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/postgis-users/attachments/20200326/17f39f55/attachment.html>