[PostGIS] #6009: Error of ST_Polygon(raster) function in PostGIS 3.4.3 and 3.5.4

[PostGIS]

#6009: Error of ST_Polygon(raster) function in PostGIS 3.4.3 and 3.5.4
-------------------------+-------------------------------------------------
  Reporter:              |      Owner:  robe
  GISuser5432            |
      Type:  defect      |     Status:  new
  Priority:  high        |  Milestone:  PostGIS 3.5.5
 Component:  raster      |    Version:  3.5.x
Resolution:              |   Keywords:  PostGIS_Raster; Raster;
                         |  ST_Polygon(Raster)
-------------------------+-------------------------------------------------
Description changed by GISuser5432:

Old description:

> I tried to run ST_Polygon(Raster) on a raster table and for 3 rows, it
> produced an error.
>
> This is the simple query:
>
> {{{#!div style="font-size: 80%"
> Code highlighting:
>   {{{#!sql
>   SELECT ST_Polygon(raster)
>   FROM error_sounding_2m
>   }}}
> }}}
>

>
> This is the error message:
>
> {{{
> ERROR:  rt_raster_surface: Could not union the pixel polygons using
> GEOSUnaryUnion()
>
> SQL state: XX000
> }}}
>

> This is the version of the PostGIS:
>
> {{{
> POSTGIS="3.4.3 e365945" [EXTENSION] PGSQL="170" GEOS="3.9.0-CAPI-1.16.2"
> PROJ="7.2.1 NETWORK_ENABLED=OFF URL_ENDPOINT=https://cdn.proj.org
> USER_WRITABLE_DIRECTORY=/var/lib/postgresql/.local/share/proj
> DATABASE_PATH=/usr/share/proj/proj.db" GDAL="GDAL 3.2.2, released
> 2021/03/05" LIBXML="2.9.10" LIBJSON="0.15" LIBPROTOBUF="1.3.3"
> WAGYU="0.5.0 (Internal)" RASTER
>
> }}}
>

>
> I also checked it for 3.5 and I had the same error. Only in 3.6 the error
> was fixed!
>

> The DB dump is attached for 3 raw of raster in a DB called enav.

New description:

 I tried to run ST_Polygon(Raster) on a raster table and for 3 rows, it
 produced an error.

 This is the simple query:

 {{{#!div style="font-size: 80%"
 Code highlighting:
   {{{#!sql
   SELECT ST_Polygon(raster)
   FROM error_sounding_2m
   }}}
 }}}



 This is the error message:

 {{{
 ERROR:  rt_raster_surface: Could not union the pixel polygons using
 GEOSUnaryUnion()

 SQL state: XX000
 }}}


 This is the version of the PostGIS:

 {{{
 POSTGIS="3.4.3 e365945" [EXTENSION] PGSQL="170" GEOS="3.9.0-CAPI-1.16.2"
 PROJ="7.2.1 NETWORK_ENABLED=OFF URL_ENDPOINT=https://cdn.proj.org
 USER_WRITABLE_DIRECTORY=/var/lib/postgresql/.local/share/proj
 DATABASE_PATH=/usr/share/proj/proj.db" GDAL="GDAL 3.2.2, released
 2021/03/05" LIBXML="2.9.10" LIBJSON="0.15" LIBPROTOBUF="1.3.3"
 WAGYU="0.5.0 (Internal)" RASTER

 }}}



 I also checked it for 3.5 and I had the same error. Only in 3.6 the error
 was fixed!

 I must also say that all tests were done in Docker containers of PostGIS.

 The DB dump is attached for 3 raw of raster in a DB called enav.

--
-- 
Ticket URL: <https://trac.osgeo.org/postgis/ticket/6009#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.