[postgis-devel] [PostGIS] #1618: [raster] ST_Transform to SRID 3978 produce wrong results

[PostGIS]

#1618: [raster] ST_Transform to SRID 3978 produce wrong results
---------------------+------------------------------------------------------
 Reporter:  pracine  |       Owner:  pracine      
     Type:  defect   |      Status:  new          
 Priority:  blocker  |   Milestone:  PostGIS 2.0.0
Component:  raster   |     Version:  trunk        
 Keywords:           |  
---------------------+------------------------------------------------------
 I have a raster in srid=4269:

 {{{
 -- Raster unprojected
 SELECT ST_AsBinary(ST_AddBand(ST_MakeEmptyRaster(10, 10, -168.0, 85.0,
 0.083, -0.083, 0, 0, 4269), 1, '8BUI', 1, 0)::geometry);
 }}}

 I project the geometric extent of this raster to srid=3978:

 {{{
 -- Raster extent projected
 SELECT ST_AsBinary(ST_Transform(ST_AddBand(ST_MakeEmptyRaster(10, 10,
 -168.0, 85.0, 0.083, -0.083, 0, 0, 4269), 1, '8BUI', 1, 0)::geometry,
 3978));
 }}}

 When I project the raster (not the vectorization of its extent) to the
 same srid=3978, it is wrongly projected:

 {{{
 --Raster projected
 SELECT ST_AsBinary(ST_Transform(ST_AddBand(ST_MakeEmptyRaster(10, 10,
 -168.0, 85.0, 0.083, -0.083, 0, 0, 4269), 1, '8BUI', 1, 0),
 3978)::geometry);
 }}}

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