[postgis-tickets] [PostGIS] #3072: ST_Intersects sfcgal implementation gives error with polyhedral surface

[PostGIS]

#3072: ST_Intersects sfcgal implementation gives error with polyhedral surface
--------------------+-------------------------------------------------------
 Reporter:  robe    |       Owner:  colivier     
     Type:  defect  |      Status:  new          
 Priority:  medium  |   Milestone:  PostGIS 2.2.0
Component:  sfcgal  |     Version:  2.1.x        
 Keywords:          |  
--------------------+-------------------------------------------------------
 I can't compare to the built in ST_Intersects since that one throws a
 polyhedralsurface not supported error.

 I have this:


 {{{
 set postgis.backend = sfcgal;
 WITH pa AS ( SELECT
 ST_Extrude(ST_Force3D(ST_GeomFromText('MULTIPOLYGON(((6011980.6923649
 2117710.26819432,6011983.40466805 2117690.92097183,
 6011988.97899699 2117691.72493859,6012014.67481019
 2117512.89736198,6011838.71290818 2117490.73975961,6011813.97754787
 2117666.48492323,
 6011961.83622488 2117687.8102276,6011958.87617183
 2117707.03845887,6011980.6923649 2117710.26819432)))',2227)),0,0,860) As
 geom ),
 drone As ( SELECT
     ST_Translate(
         ST_Force3D(
             ST_Transform(ST_GeomFromText('POINT(-122.4028
 37.7952)',4326),2227) ),
         0, 0, z) As geom, z
  FROM generate_series(500,1000, 500) As z )
 SELECT z
 FROM drone INNER JOIN pa ON ST_Intersects(drone.geom, pa.geom);

 }}}


 I get this error:


 {{{
 NOTICE:  During intersects(A,B) :
 NOTICE:    with A: POINT(6455232916477051/1073741824
 2273757146426363/1073741824 500/1)
 NOTICE:     and B:
 GEOMETRYCOLLECTION(POLYHEDRALSURFACE(((6455315114472671/1073741824
 4547748172148997/2147483648 0/1,3227659013393001/536870912
 4547706624305065/2147483648 0/1,3227662006088063/536870912
 568463543851317/268435456 0/1,1613837900711365/268435456 227


 }}}

 This is using:


 {{{
 POSTGIS="2.2.0dev r13298" GEOS="3.5.0dev-CAPI-1.9.0 r4048" SFCGAL="1.0.5"
 PROJ="Rel. 4.8.0, 6 March 2012" GDAL="GDAL 1.11.1, released 2014/09/24"
 LIBXML="2.7.8" LIBJSON="0.12" RASTER
 }}}


 I was expecting this to give me the footprint intersection.

 On another note -- I thought if I used ST_Extrude in 2.2, that the
 ST_3DInteresects would give me a solid intersection (bit would be set to
 solid), it seems not to.  I'll bring this up on dev list.

-- 
Ticket URL: <http://trac.osgeo.org/postgis/ticket/3072>
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.