[PostGIS] #6071: ST_intersects fails for srid 3035

[PostGIS]

#6071: ST_intersects fails for srid 3035
-------------------------------+---------------------------
 Reporter:  Lars Aksel Opsahl  |      Owner:  pramsey
     Type:  defect             |     Status:  new
 Priority:  medium             |  Milestone:  PostGIS 3.6.3
Component:  postgis            |    Version:  3.6.x
 Keywords:                     |
-------------------------------+---------------------------
 We have two objects that intersect in SRID 4258. The intersection area is
 not very small, but approximately 21.96131633594632.

 {{{
 SELECT ST_area(ST_intersection(
 ST_transform(ST_SetSRID('0103000020DB0B000001000000050000001DBF0EB1DF4F514104447475EF7D5341112F8BDD854E51412DDFB3CF638D53411BE3F8F2525651411CF8786E038E5341E0AB8694B95751419CB58506907E53411DBF0EB1DF4F514104447475EF7D5341'::geometry,
 3035),4258),
 ST_transform(ST_SetSRID('0103000020DB0B0000010000001E000000086879984555514154EEF8B95D7E5341F3BB2D3E465551419A5D6AB75C7E5341CDFEF858465551418E08FC6B5B7E5341D464B6F045555141F1BD3C2F5A7E5341D1D9AE48455551415361814C597E5341F0E7FA7944555141EED77BA6587E5341F0E7FA7944555141EED77BA6587E5341C1A0F1AF44555141223FC4CC577E5341C1A0F1AF44555141223FC4CC577E5341FDEF3EEE475551414CA0B94F587E5341056A0B4B4955514188FF7256587E534137DCC1DC4955514105356B91567E534137DCC1DC4955514105356B91567E5341312E151A4A555141FA90CC84557E5341312E151A4A555141FA90CC84557E5341E22B94504D5551411D733750567E53419C118C3452555141C434AF60577E534166D164EF55555141B09040ED577E534171EE750F59555141CC8126DD587E53413A4DD0615F55514194AA387C5A7E53413A4DD0615F55514194AA387C5A7E5341E2DAF2265F5551418CFD050A5B7E5341F8E4D65C5E55514173DEA6955C7E53413A9657E05D555141C20DD6E45D7E5341DAA090085C5551418E8FB4515E7E534108A0C74B5A5551412A4B86415E7E53419ED99C7A57555141369F0DCE5D7E5341D3E490795355514155680BEB5D7E5341F2D1AC1E4F555141A69091AA5D7E5341086879984555514154EEF8B95D7E5341'::geometry,
 3035),4258)
 ),true) ;
       st_area
 -------------------
  21.96131633594632
 (1 row)

 }}}

 Then try the same with SRID 3035 and zero area

 {{{
 SELECT ST_area(ST_intersection(
 ST_SetSRID('0103000020DB0B000001000000050000001DBF0EB1DF4F514104447475EF7D5341112F8BDD854E51412DDFB3CF638D53411BE3F8F2525651411CF8786E038E5341E0AB8694B95751419CB58506907E53411DBF0EB1DF4F514104447475EF7D5341'::geometry,
 3035),
 ST_SetSRID('0103000020DB0B0000010000001E000000086879984555514154EEF8B95D7E5341F3BB2D3E465551419A5D6AB75C7E5341CDFEF858465551418E08FC6B5B7E5341D464B6F045555141F1BD3C2F5A7E5341D1D9AE48455551415361814C597E5341F0E7FA7944555141EED77BA6587E5341F0E7FA7944555141EED77BA6587E5341C1A0F1AF44555141223FC4CC577E5341C1A0F1AF44555141223FC4CC577E5341FDEF3EEE475551414CA0B94F587E5341056A0B4B4955514188FF7256587E534137DCC1DC4955514105356B91567E534137DCC1DC4955514105356B91567E5341312E151A4A555141FA90CC84557E5341312E151A4A555141FA90CC84557E5341E22B94504D5551411D733750567E53419C118C3452555141C434AF60577E534166D164EF55555141B09040ED577E534171EE750F59555141CC8126DD587E53413A4DD0615F55514194AA387C5A7E53413A4DD0615F55514194AA387C5A7E5341E2DAF2265F5551418CFD050A5B7E5341F8E4D65C5E55514173DEA6955C7E53413A9657E05D555141C20DD6E45D7E5341DAA090085C5551418E8FB4515E7E534108A0C74B5A5551412A4B86415E7E53419ED99C7A57555141369F0DCE5D7E5341D3E490795355514155680BEB5D7E5341F2D1AC1E4F555141A69091AA5D7E5341086879984555514154EEF8B95D7E5341'::geometry,
 3035)
 )) ;
  st_area
 ---------
        0
 (1 row)

 }}}

 A even simple intersects fails here.


 {{{
 SELECT ST_intersects(
 ST_SetSRID('0103000020DB0B000001000000050000001DBF0EB1DF4F514104447475EF7D5341112F8BDD854E51412DDFB3CF638D53411BE3F8F2525651411CF8786E038E5341E0AB8694B95751419CB58506907E53411DBF0EB1DF4F514104447475EF7D5341'::geometry,
 3035),
 ST_SetSRID('0103000020DB0B0000010000001E000000086879984555514154EEF8B95D7E5341F3BB2D3E465551419A5D6AB75C7E5341CDFEF858465551418E08FC6B5B7E5341D464B6F045555141F1BD3C2F5A7E5341D1D9AE48455551415361814C597E5341F0E7FA7944555141EED77BA6587E5341F0E7FA7944555141EED77BA6587E5341C1A0F1AF44555141223FC4CC577E5341C1A0F1AF44555141223FC4CC577E5341FDEF3EEE475551414CA0B94F587E5341056A0B4B4955514188FF7256587E534137DCC1DC4955514105356B91567E534137DCC1DC4955514105356B91567E5341312E151A4A555141FA90CC84557E5341312E151A4A555141FA90CC84557E5341E22B94504D5551411D733750567E53419C118C3452555141C434AF60577E534166D164EF55555141B09040ED577E534171EE750F59555141CC8126DD587E53413A4DD0615F55514194AA387C5A7E53413A4DD0615F55514194AA387C5A7E5341E2DAF2265F5551418CFD050A5B7E5341F8E4D65C5E55514173DEA6955C7E53413A9657E05D555141C20DD6E45D7E5341DAA090085C5551418E8FB4515E7E534108A0C74B5A5551412A4B86415E7E53419ED99C7A57555141369F0DCE5D7E5341D3E490795355514155680BEB5D7E5341F2D1AC1E4F555141A69091AA5D7E5341086879984555514154EEF8B95D7E5341'::geometry,
 3035)
 ) ;
  st_intersects
 ---------------
  f
 (1 row)

 }}}

 This tested on


 {{{
 POSTGIS="3.4.2 c19ce56" [EXTENSION] PGSQL="160" GEOS="3.12.1-CAPI-1.18.1"
 (compiled against GEOS 3.10.2) SFCGAL="SFCGAL 1.4.1, CGAL 5.3.1, BOOST
 1.74.0" PROJ="8.2.1 NETWORK_ENABLED=OFF URL_ENDPOINT=https://cdn.proj.org
 USER_WRITABLE_DIRECTORY=/tmp/proj DATABASE_PATH=/usr/share/proj/proj.db"
 GDAL="GDAL 3.4.1, released 2021/12/27" LIBXML="2.9.13" LIBJSON="0.15"
 LIBPROTOBUF="1.3.3" WAGYU="0.5.0 (Internal)" TOPOLOGY RASTER

 POSTGIS="3.7.0dev 3.6.0rc2-305-g3b3488ddf" [EXTENSION] PGSQL="160"
 GEOS="3.14.0-CAPI-1.20.4" PROJ="9.7.0 NETWORK_ENABLED=OFF
 URL_ENDPOINT=https://cdn.proj.org
 USER_WRITABLE_DIRECTORY=/Users/lop/Library/Application Support/proj
 DATABASE_PATH=/opt/homebrew/Cellar/proj/9.7.0/share/proj/proj.db"
 (compiled against PROJ 9.7.0) LIBXML="2.13.0" LIBJSON="0.18" TOPOLOGY

 POSTGIS="3.7.0dev 3.6.0rc2-435-g2c9277d3d" [EXTENSION] PGSQL="160"
 GEOS="3.12.1-CAPI-1.18.1" PROJ="8.2.1 NETWORK_ENABLED=OFF
 URL_ENDPOINT=https://cdn.proj.org USER_WRITABLE_DIRECTORY=/tmp/proj
 DATABASE_PATH=/usr/share/proj/proj.db" (compiled against PROJ 8.2.1)
 LIBXML="2.9.13" LIBJSON="0.15" LIBPROTOBUF="1.3.3" WAGYU="0.5.0
 (Internal)" TOPOLOGY

 }}}
-- 
Ticket URL: <https://trac.osgeo.org/postgis/ticket/6071>
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.