[postgis-devel] [PostGIS] #1536: [raster] There should be an option for certain variants of ST_Intersection(raster, raster)

[PostGIS]

#1536: [raster] There should be an option for certain variants of
ST_Intersection(raster, raster)
---------------------+------------------------------------------------------
 Reporter:  pracine  |       Owner:  pracine      
     Type:  defect   |      Status:  new          
 Priority:  low      |   Milestone:  PostGIS 2.0.0
Component:  raster   |     Version:  trunk        
 Keywords:           |  
---------------------+------------------------------------------------------

Comment(by pracine):

 I would expect, if I intersect a raster representing a circle surrounded
 with nodata values with another raster, that the area surrounding the
 circle would be filled with nodata values in BOTH bands...

 Otherwise we can not say that ST_Intersection is taking nodata value into
 account. It is just an extent intersection.

 Actually I guess this is what is happening based on the way you are
 calling ST_MapAlgebraExpr: since you are not providing alternate
 expressions, any time one pixel from any band is nodata both resulting
 bands gets nodata. This should be the expected behavior of ST_MapAlgebra.
 I check...

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