Georeferencing Mrsid raster file with a world file

Volkan Ulukut arthan at DSL.TTNET.NET.TR
Sun Nov 11 21:56:54 PST 2007


Greetings list,

I've encountered a problem while trying to publish a sid file. The sid file
does not show in the extent I'm trying to show instead it has a fixed extent
and shows there. 

I got 3 files with the same prefix

kayseri_lokal.sid

kayseri_lokal.sdw

kayseri_lokal.aux

 

I've added this to my mapfile

 

 

LAYER

    NAME "kay_uydu7"

    STATUS ON

    TYPE RASTER

    MAXSCALE 24999

    DATA "uydu/kayseri_lokal.sid"

    GROUP "kay_uydu7"

   METADATA

     "queryable" "true"

    "imageformat" "alpha"

   END

    TRANSPARENCY alpha

    CLASS

    END

END

 

The extent of my mapfile is:

EXTENT 440773 4282202 467369 4295904

 

 

This is the contents of my sdw file which should put the raster on the right
extent:

 

             1.000000000000000

             0.000000000000000

             0.000000000000000

            -1.000000000000000

        440773.000000000000000

        467369.000000000000000

 

 

I've also tried making this a wld file in case gdal does not recognize sdw
file. 

When I type gdalinfo on sid file I get:

 

Driver: MrSID/Multi-resolution Seamless Image Database (MrSID)

Size is 55564, 66276

Coordinate System is `'

Origin = (75393.500000000000000,137653.500000000000000)

Pixel Size = (1.000000000000000,-1.000000000000000)

Metadata:

  IMAGE__INPUT_NAME=j:/kayseri_ik_lokal_v5/mrsid_lokal.img

  IMAGE__INPUT_FILE_SIZE=11047678992.000000

  GEOTIFF_NUM__1024__GTModelTypeGeoKey=0

  GEOTIFF_CHAR__GTModelTypeGeoKey=Undefined

  GEOTIFF_NUM__1025__GTRasterTypeGeoKey=1

  GEOTIFF_CHAR__GTRasterTypeGeoKey=RasterPixelIsArea

  IMAGE__DYNAMIC_RANGE_LEVEL=127.500000

  GEOTIFF_NUM__3076__ProjLinearUnitsGeoKey=9014

  GEOTIFF_CHAR__ProjLinearUnitsGeoKey=Linear_Fathom

  IMAGE__DYNAMIC_RANGE_WINDOW=255.000000

  GEO__ModelTypeGeoKey=0

  IMG__HORIZONTAL_UNITS=other

  IMAGE__COMPRESSION_VERSION=1,6,-3

  IMAGE__TARGET_COMPRESSION_RATIO=20.000000

  IMAGE__COMPRESSION_NLEV=8

  IMAGE__COMPRESSION_WEIGHT=4.000000

  IMAGE__COMPRESSION_GAMMA=2.000000

  IMAGE__COMPRESSION_BLOCK_SIZE=512

  IMAGE__CREATION_DATE=Fri Aug 03 13:21:10 2007

 

  IMAGE__WIDTH=55564

  IMAGE__HEIGHT=66276

  IMAGE__COLOR_SCHEME=0

  IMAGE__DATA_TYPE=0

  IMAGE__BITS_PER_SAMPLE=8

  IMAGE__XY_ORIGIN=75394.000000,137653.000000

  IMAGE__X_RESOLUTION=1.000000

  IMAGE__Y_RESOLUTION=1.000000

Corner Coordinates:

Upper Left  (   75393.500,  137653.500)

Lower Left  (   75393.500,   71377.500)

Upper Right (  130957.500,  137653.500)

Lower Right (  130957.500,   71377.500)

Center      (  103175.500,  104515.500)

Band 1 Block=1024x128 Type=Byte, ColorInterp=Red

  Overviews: 27782x33138, 13891x16569, 6946x8285, 3473x4143, 1737x2072,
869x1036, 435x518, 218x259

Band 2 Block=1024x128 Type=Byte, ColorInterp=Green

  Overviews: 27782x33138, 13891x16569, 6946x8285, 3473x4143, 1737x2072,
869x1036, 435x518, 218x259

Band 3 Block=1024x128 Type=Byte, ColorInterp=Blue

  Overviews: 27782x33138, 13891x16569, 6946x8285, 3473x4143, 1737x2072,
869x1036, 435x518, 218x259

 

 

When I change the mapfile extent to:

EXTENT 75393 71377 137653 130957

I can see the raster layer normally.But I need to move it to my initial
mapfile extent (since all the other data is there)

Thanks in advance

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20071112/207f253a/attachment.htm>


More information about the MapServer-users mailing list