[postgis-devel] ST_3DIntersects segfaults

[Tom van Tilburg]

The query below segfaults with me on
"POSTGIS="2.5.0rc1dev r16693" [EXTENSION] PGSQL="100"
GEOS="3.7.0beta2-CAPI-1.11.0 e7d54545" SFCGAL="1.3.5" PROJ="Rel. 5.1.0,
June 1st, 2018" GDAL="GDAL 2.4.0dev-3ea4d4ed14-dirty, released 2018/07/09"
LIBXML="2.9.4" LIBJSON="0.12.1" LIBPROTOBUF="1.2.1" RAST (...)"

I don't have another postgis version at hand to test, can others confirm
that?
The polyhedral surface seems valid, removing decimal precision makes the
query run without segfaulting.

WITH data as (
SELECT
ST_GeometryFromText(
'POLYGON Z ((
122395.299 489126.697 8.61546664325712,
122389.298 489128.73 8.55588025324629,
122391.489 489135.198 8.5526708028059,
122397.49 489133.165 8.61225719281685,
122395.299 489126.697 8.61546664325712))'
) as geoma,
ST_GeometryFromText(
'POLYHEDRALSURFACE Z (((
122390.998245685 489133.068537491 0,
122391.003145022 489133.066423547 0,
122391.003145022 489133.066423547 10,
122390.998245685 489133.068537491 10,
122390.998245685 489133.068537491 0
)),((
122391.003145022 489133.066423547 0,
122383.269575402 489114.842869866 0,
122383.269575402 489114.842869866 10,
122391.003145022 489133.066423547 10,
122391.003145022 489133.066423547 0
)))'
) as geomb
)
SELECT ST_3DIntersection(geoma, geomb)
FROM data

Best,
 Tom
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/postgis-devel/attachments/20180823/4292b6f2/attachment.html>