ST_Difference drops z-dimension on empty result

[lkajan]

When calling ST_Difference with two 3D geometries the result may be a 2D geometry in case when geometries are equal and result is an EMPTY geometry.

Expected behavior:
As the documentation states: "This function supports 3d and will not drop the z-index."

To reproduce:

SELECT ST_AsText(ST_Difference('LINESTRINGZ(0 0 0,1 0 0)'::geometry, 'LINESTRINGZ(0 0 0,1 0 0)'::geometry));
    st_astext
------------------
 LINESTRING EMPTY
(1 row)

The result is a LINESTRING EMPTY and not a LINESTRINGZ EMPTY as one would expect.

PostGIS version:
POSTGIS="3.4.0 0874ea3" [EXTENSION] PGSQL="160" 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" LIBXML="2.9.10" LIBJSON="0.15" LIBPROTOBUF="1.3.3" WAGYU="0.5.0 (Internal)" TOPOLOGY

[pramsey]

Hum, this has been fixed at some point, either in PostGIS or GEOS (I'm guessing the latter).

# SELECT ST_AsText(ST_Difference('LINESTRINGZ(0 0 0,1 0 0)'::geometry, 'LINESTRINGZ(0 0 0,1 0 0)'::geometry));
     st_astext      
--------------------
 LINESTRING Z EMPTY

POSTGIS="3.5.0dev"
GEOS="3.13.0dev"

[robe]

Works for me too on GEOS 3.12.2dev LINESTRING Z EMPTY

Going to close this out as GEOS.

Ikajan your GEOS is a very old micro. The latest release GEOS for 3.9 is 3.9.5 ​https://libgeos.org/usage/download/. If you can upgrade that that would be the best.