[mapserver-users] Using Spot DIMAP Format in Mapserver

[Eric Francois]

Hello,

I'm trying to use Dimap format from Spot Image( metadata.dim) in my mapfile
but my raster is always bad located.
METADATA.DIM opens directly in Qgis is well placed and well oriented. The
correspondant IMAGERY.TIF is well placed but not well oriented.
I'm not able to create the worldfile (pb to catch orientation info) for my
tif file so I want to use directly dimap format but I don't understand why
metadata.dim is not well interpreted by Mapserver.

My Mapfile contains:
TYPE RASTER
DATA "country/METADATA.DIM"
PROJECTION
	'proj=longlat'
	'ellps=WGS84'
	'datum=WGS84'
	'no_defs'
END

Gdalinfo on IMAGERY.TIF:
Driver: GTiff/GeoTIFF
Files: IMAGERY.TIF
Size is 3000, 3000
Coordinate System is:
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.2572235630016,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4326"]]
Origin = (11.795597991857200,-4.699265933053320)
Pixel Size = (0.000070782000000,-0.000070782000000)
Metadata:
  AREA_OR_POINT=Point
  TIFFTAG_IMAGEDESCRIPTION=B3 B2 B1 B4
  TIFFTAG_DATETIME=20090428 15:04:18
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  11.7955980,  -4.6992659)
Lower Left  (  11.7955980,  -4.9116119)
Upper Right (  12.0079440,  -4.6992659)
Lower Right (  12.0079440,  -4.9116119)
Center      (  11.9017710,  -4.8054389)
Band 1 Block=3000x500 Type=Byte, ColorInterp=Red Band 2 Block=3000x500
Type=Byte, ColorInterp=Green Band 3 Block=3000x500 Type=Byte,
ColorInterp=Blue Band 4 Block=3000x500 Type=Byte, ColorInterp=Undefined

Gdalinfo on the METADATA.DIM of IMAGERY.TIF

Driver: DIMAP/SPOT DIMAP
Files: METADATA.DIM
Size is 3000, 3000
Coordinate System is `'
GCP Projection = GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS
84",6378137,298.
257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenw
ich"
,0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EP
SG",
"9122"]],AUTHORITY["EPSG","4326"]]
GCP[  0]: Id=1, Info=
          (0.5,0.5) -> (11.795633,-4.699301,0) GCP[  1]: Id=2, Info=
          (2999.5,0.5) -> (12.010877,-4.734618,0) GCP[  2]: Id=3, Info=
          (0.5,2999.5) -> (11.745949,-4.913964,0) GCP[  3]: Id=4, Info=
          (2999.5,2999.5) -> (11.961267,-4.949299,0)
Metadata:
  JOB_ID=R2_CAT_090428123250353_1
  PRODUCT_INFO=FORMOSAT-2 SCENE level 1A
  PRODUCT_TYPE=FORMOSAT-2 SCENE
  DATASET_PRODUCER_NAME=NSPO / Spot Image
  DATASET_PRODUCTION_DATE=2009-04-28 15:04:20.000000
  GRID_REFERENCE=0067-0481
  SHIFT_VALUE=0
  IMAGING_DATE=2006-02-28
  IMAGING_TIME=08:57:15.197455
  MISSION=FORMOSAT
  MISSION_INDEX=2
  INSTRUMENT=RSI
  INSTRUMENT_INDEX=1
  IMAGING_MODE=MS
  SCENE_PROCESSING_LEVEL=0
  VIEWING_ANGLE_ALONG_TRACK=0.273794
  VIEWING_ANGLE_ACROSS_TRACK=-3.151896
  SATELLITE_INCIDENCE_ANGLE=3.611694
  SATELLITE_AZIMUTH_ANGLE=98.138828
  SUN_AZIMUTH=98.126829
  SUN_ELEVATION=52.988591
  REVOLUTION_NUMBER=09085
  PROCESSING_LEVEL=1A
  GEOMETRIC_PROCESSING=SYSTEM
  RADIOMETRIC_PROCESSING=SYSTEM
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0, 3000.0)
Upper Right ( 3000.0,    0.0)
Lower Right ( 3000.0, 3000.0)
Center      ( 1500.0, 1500.0)
Band 1 Block=3000x500 Type=Byte, ColorInterp=Red Band 2 Block=3000x500
Type=Byte, ColorInterp=Green Band 3 Block=3000x500 Type=Byte,
ColorInterp=Blue Band 4 Block=3000x500 Type=Byte, ColorInterp=Undefined


Thanks for all 

ERIC