[postgis-tickets] [PostGIS] #2586: [raster] ST_Clip with geometry smaller than pixel size

[PostGIS]

#2586: [raster] ST_Clip with geometry smaller than pixel size
-----------------------+----------------------------------------------------
 Reporter:  dustymugs  |       Owner:  dustymugs    
     Type:  defect     |      Status:  new          
 Priority:  medium     |   Milestone:  PostGIS 2.1.2
Component:  raster     |     Version:  2.1.x        
 Keywords:             |  
-----------------------+----------------------------------------------------

Comment(by dustymugs):

 This occurs with a geometry that is sufficiently smaller than the pixel

 {{{
 ST_Multi(ST_Buffer('SRID=910001;POINT(624772 1534404)', 10))
 }}}

 The following geometry works.

 {{{
 ST_Multi(ST_Buffer('SRID=910001;POINT(624772 1534404)', 85))
 }}}

 {{{
 WITH foo AS (
         SELECT ST_Multi(ST_Buffer('SRID=910001;POINT(624772 1534404)',
 80)) AS geom
 )
 SELECT
         ST_DumpValues(ST_Clip(r.rast, foo.geom))
 FROM test_raster r
 JOIN foo
         ON ST_Intersects(r.rast, foo.geom)
 }}}

 The test raster was loaded with the following. The SRID is something user-
 specified for the projection of the raster.

 {{{
 raster2pgsql -s 910001 -t 1073x1 -R -C /tmp/ds.maxt.bin test_raster | psql
 }}}

-- 
Ticket URL: <http://trac.osgeo.org/postgis/ticket/2586#comment:1>
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.