[postgis-devel] [PostGIS] #925: [raster] ST_Transform

[PostGIS]

#925: [raster] ST_Transform
----------------------------+-----------------------------------------------
 Reporter:  dustymugs       |       Owner:  dustymugs    
     Type:  task            |      Status:  assigned     
 Priority:  medium          |   Milestone:  PostGIS 2.0.0
Component:  postgis raster  |     Version:  trunk        
 Keywords:  history         |  
----------------------------+-----------------------------------------------

Comment(by robe):

 I'm using a fairly recent GDAL trunk -- 1.9 SVN r22508 (so that I could
 get the changes that you and Jorge made to GDAL)


 Tried one of your new ones.  Sadly that doesn't seem to work either -- so
 it seems it's not transforming on mine but silently failing and returning
 the same thing.

 {{{
 SELECT
         ST_SRID(rast),
         ST_Width(rast),
         ST_Height(rast),
         ST_NumBands(rast),
         round(ST_ScaleX(rast)::numeric, 3),
         round(ST_ScaleY(rast)::numeric, 3),
         ST_Count(rast),
         round(ST_Mean(rast)::numeric, 3)
 FROM (SELECT ST_Transform(
         ST_SetValue(
                 ST_SetValue(
                         ST_SetValue(
                                 ST_AddBand(
                                         ST_MakeEmptyRaster(10, 10,
 -500000, 600000, 1000, 1000, 0, 0, 1003310)
                                         , 1, '64BF', 0, 0
                                 )
                                 , 1, 1, 1, -10
                         )
                         , 1, 5, 4, 0
                 )
                 , 1, 5, 5, 3.14159
         )
         , 1004269) AS rast
 ) foo;
 }}}

 Gives pretty much what I started with:

 {{{
 1004269 10      10      1       1000.000        -1000.000       2
 -3.429
 }}}

 I'll try compiling proj 4.7 and see if it behaves differently or it still
 needs those files for some reason.  Perhaps I'll set the proj data path on
 my box to see if it changes the behavior first.  I have been putting off
 compiling proj 4.7 just because it's cumbersome to compile under mingw
 windows -- always have to hack it (come to think of it only been succesful
 compiling it under my mingw-64 install though haven't tried for a while on
 my regular mingw.

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