[postgis-devel] [PostGIS] #1402: [raster] ST_Resample producing incorrect transformation

[PostGIS]

#1402: [raster] ST_Resample producing incorrect transformation
----------------------------+-----------------------------------------------
 Reporter:  dustymugs       |       Owner:  dustymugs    
     Type:  defect          |      Status:  new          
 Priority:  blocker         |   Milestone:  PostGIS 2.0.0
Component:  postgis raster  |     Version:  trunk        
 Keywords:                  |  
----------------------------+-----------------------------------------------

Comment(by dustymugs):

 Fascinating.  Depending on the format passed to GDAL's gdaltransform, the
 resulting answer varies

 {{{
 gdaltransform -s_srs 'PROJCS["Amersfoort / RD
 New",GEOGCS["Amersfoort",DATUM["Amersfoort",SPHEROID["Bessel
 1841",6377397.155,299.1528128,AUTHORITY["EPSG","7004"]],AUTHORITY["EPSG","6289"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4289"]],UNIT["metre",1,AUTHORITY["EPSG","9001"]],PROJECTION["Oblique_Stereographic"],PARAMETER["latitude_of_origin",52.15616055555555],PARAMETER["central_meridian",5.38763888888889],PARAMETER["scale_factor",0.9999079],PARAMETER["false_easting",155000],PARAMETER["false_northing",463000],AUTHORITY["EPSG","28992"],AXIS["X",EAST],AXIS["Y",NORTH]]'
 -t_srs 'PROJCS["Popular Visualisation CRS / Mercator
 (deprecated)",GEOGCS["Popular Visualisation
 CRS",DATUM["Popular_Visualisation_Datum",SPHEROID["Popular Visualisation
 Sphere",6378137,0,AUTHORITY["EPSG","7059"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6055"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4055"]],UNIT["metre",1,AUTHORITY["EPSG","9001"]],PROJECTION["Mercator_1SP"],PARAMETER["central_meridian",0],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],AUTHORITY["EPSG","3785"],AXIS["X",EAST],AXIS["Y",NORTH]]'
 195000 450000
 664595.85121482 6773064.24012838 -13206.6219004256
 }}}

 {{{
 gdaltransform -s_srs "+proj=sterea +lat_0=52.15616055555555
 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel
 +towgs84=565.237,50.0087,465.658,-0.406857,0.350733,-1.87035,4.0812
 +units=m +no_defs" -t_srs "+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0
 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs"
 195000 450000
 664595.859810406 6806797.32452805 43.5179952728836
 }}}

 {{{
 gdaltransform -s_srs EPSG:28992 -t_srs 'PROJCS["Popular Visualisation CRS
 / Mercator (deprecated)",GEOGCS["Popular Visualisation
 CRS",DATUM["Popular_Visualisation_Datum",SPHEROID["Popular Visualisation
 Sphere",6378137,0,AUTHORITY["EPSG","7059"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6055"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4055"]],UNIT["metre",1,AUTHORITY["EPSG","9001"]],PROJECTION["Mercator_1SP"],PARAMETER["central_meridian",0],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],AUTHORITY["EPSG","3785"],AXIS["X",EAST],AXIS["Y",NORTH]]'
 195000 450000
 664595.85126645 6773064.23984736 -13206.6222764212
 }}}

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