<html><head></head><body><div>Yes, I can confirm crashing </div><div><br></div><div>postgis_full_version</div><div>POSTGIS="2.4.4 r16526" GEOS="3.7.0dev-CAPI-1.11.0 084542e5" SFCGAL="1.3.0" PROJ="Rel. 4.9.3, 15 August 2016" GDAL="GDAL 2.1.2, released 2016/10/24" LIBXML="2.9.4" LIBJSON="0.12.99" LIBPROTOBUF="1.2.1" (core procs from "2.4.0dev r15351" need upgrade) RASTER (raster procs from "2.4.0dev r15351" need upgrade) (sfcgal procs from "2.4.0dev r15351" need upgrade)</div><div><br></div><div>That is handled in sfcgal I guess.</div><div><br></div><div>/Nicklas</div><div><br></div><div><br></div><div>On Thu, 2018-08-23 at 14:04 +0200, Tom van Tilburg wrote:</div><blockquote type="cite"><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>
<pre>_______________________________________________
postgis-devel mailing list
<a href="mailto:postgis-devel@lists.osgeo.org">postgis-devel@lists.osgeo.org</a>
<a href="https://lists.osgeo.org/mailman/listinfo/postgis-devel">https://lists.osgeo.org/mailman/listinfo/postgis-devel</a></pre></blockquote></body></html>