<div dir="ltr"><div>The query below segfaults with me on <br></div><div><span style="font-family:monospace,monospace">"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 (...)"</span><br></div><div></div><div></div><div><br></div><div>I don't have another postgis version at hand to test, can others confirm that?</div><div></div><div>The polyhedral surface seems valid, removing decimal precision makes the query run without segfaulting.<br></div><div><br></div><div><span style="font-family:monospace,monospace">WITH data as (<br>SELECT<br>ST_GeometryFromText(<br>'POLYGON Z ((<br>122395.299 489126.697 8.61546664325712,<br>122389.298 489128.73 8.55588025324629,<br>122391.489 489135.198 8.5526708028059,<br>122397.49 489133.165 8.61225719281685,<br>122395.299 489126.697 8.61546664325712))'<br>) as geoma,<br>ST_GeometryFromText(<br>'POLYHEDRALSURFACE Z (((<br>122390.998245685 489133.068537491 0,<br>122391.003145022 489133.066423547 0,<br>122391.003145022 489133.066423547 10,<br>122390.998245685 489133.068537491 10,<br>122390.998245685 489133.068537491 0<br>)),((<br>122391.003145022 489133.066423547 0,<br>122383.269575402 489114.842869866 0,<br>122383.269575402 489114.842869866 10,<br>122391.003145022 489133.066423547 10,<br>122391.003145022 489133.066423547 0<br>)))'<br>) as geomb<br>)<br>SELECT ST_3DIntersection(geoma, geomb)<br>FROM data</span><br></div><div><br></div><div>Best,</div><div> Tom<br></div></div>