[postgis-devel] [PostGIS] #1567: [raster] The SRID of the empty raster returned by ST_MapAlgebraExpr(rast, rast, 'intersection') should be 0

[PostGIS]

#1567: [raster] The SRID of the empty raster returned by ST_MapAlgebraExpr(rast,
rast, 'intersection') should be 0
----------------------+-----------------------------------------------------
  Reporter:  pracine  |       Owner:  pracine      
      Type:  defect   |      Status:  closed       
  Priority:  low      |   Milestone:  PostGIS 2.0.0
 Component:  raster   |     Version:  trunk        
Resolution:  invalid  |    Keywords:               
----------------------+-----------------------------------------------------
Changes (by pracine):

  * status:  new => closed
  * resolution:  => invalid


Comment:

 As a last test I would check the SRID of the empty geometry resulting from
 an intersection. If it's the same as the input geometry then no. If it's 0
 then yes.

 SELECT ST_SRID(ST_Intersection(ST_Buffer(ST_SetSRID(ST_MakePoint(1.5,
 -1.5), 4269), 0.2), ST_Buffer(ST_SetSRID(ST_MakePoint(4.5, -1.5), 4269),
 0.2)));

 Gives 4269... so it should be the same as the input raster...

 I set it to invalid.

 However:

 {{{
 SELECT ST_SRID(ST_MakeEmptyRaster(4, 4, 0, 0, 1, -1, 0, 0, -10))
 }}}

 should return 0 no? Since:

 {{{
 SELECT ST_SRID(ST_SetSRID(ST_MakePoint(4.5, -1.5), -4269));
 }}}

 Does...

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