[Mapserver-users] Raster: rotation params in WLD file - Resampling

[Frank Warmerdam]

Marian,

OK, I tracked through and it seems that since MapServer 3.6 I broket the
support for resampling raster files with non-zero rotation/shear coefficients
in the world file.  I have corrected the MapServer 3.7 code in CVS and added
a test to the auto-test suite (msautotest/gdal/wld_rotate.map).

I now get correct operation with the following mapfile and world file. Thank
you for preparing the concise example for me to work with.  The key issue a
user needs to keep in mind is to specify a projection in the MAP and LAYER (even
if they are just a dummy one like here) or the resampling logic can't kick in.

try.map:
NAME EX1
EXTENT 0 -558 880 322
SIZE 645 645

PROJECTION
    "proj=utm"
    "zone=12"
    "datum=WGS84"
END

IMAGECOLOR 0 200 200

LAYER
   NAME "map1"
   DATA "map1.tif"
   TYPE RASTER
   STATUS DEFAULT
   DEBUG ON
   PROJECTION
     "proj=utm"
     "zone=12"
     "datum=WGS84"
   END
END # end of layer object
END # end of map file

map1.wld:
0.8660254
0.5
0.5
-0.8660254
0
0

Best regards,

-- 
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up   | Frank Warmerdam, warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Programmer for Rent