[postgis-tickets] [PostGIS] #4898: ST_SetZ seems to hang indefinitely if a vertex is not in raster

[PostGIS]

#4898: ST_SetZ seems to hang indefinitely if a vertex is not in raster
---------------------+---------------------------
 Reporter:  robe     |      Owner:  pramsey
     Type:  defect   |     Status:  new
 Priority:  blocker  |  Milestone:  PostGIS 3.2.0
Component:  postgis  |    Version:  master
 Keywords:           |
---------------------+---------------------------
 Note this example is very close to what is in the manual but with an added
 vertex not in raster


 {{{
 WITH test_raster AS (
 SELECT
 ST_SetValues(
   ST_AddBand(
     ST_MakeEmptyRaster(width => 2, height => 2,
       upperleftx => 0, upperlefty => 2,
       scalex => 1.0, scaley => -1.0,
       skewx => 0, skewy => 0, srid => 4326),
     index => 1, pixeltype => '16BSI',
     initialvalue => 0,
     nodataval => -999),
   1,1,1,
   newvalueset =>ARRAY[ARRAY[10.0::float8, 50.0::float8],
 ARRAY[40.0::float8, 20.0::float8]]) AS rast
 )
 SELECT
 ST_AsText(
   ST_SetZ(
     rast,
     band => 1,
     geom => 'SRID=4326;LINESTRING(1.0 1.9, 1.0 0.2, 10 100)'::geometry,
     resample => 'bilinear'
 ))
 FROM test_raster
 }}}



 After loosing my patience waiting forf an answer, I canceled, and saw tons
 of these notices:


 {{{
 NOTICE:  Attempting to get pixel value with out of range raster
 coordinates: (10, -98)
 NOTICE:  Attempting to get pixel value with out of range raster
 coordinates: (10, -98)
 NOTICE:  Attempting to get pixel value with out of range raster
 coordinates: (10, -98)
 NOTICE:  Attempting to get pixel value with out of range raster
 coordinates: (10, -98)
 }}}

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