Georeferencing Mrsid raster file with a world file

Mon Nov 12 04:04:07 EST 2007

I've renamed world file and the extent is stil the same. So I assume there
is a native georeferencing data as you call. But the thing is the dataset
already has sdw file so there shouldn't be any conflict with gdal. I've
googled a little and find out that sdw world files are recognized by gdal
and wld files are recognized by mapserver so i tried both and failed to
change the extent. 

MrSID tools change the sdw file not the embedded data I'm afraid.

I cannot change the format because of the performance issues. The file that
big might cause some problems. 

I'm not a GIS expert but these layers work fine with GIS software which
means their projection and registries match. But GIS software take sdw file
into account. Mapserver should do that too. This is either a bug on
mapserver - gdal side or I don't have a correctly configured gdal -
mapserver which I doubt because I can see the raster image when I set my
mapfile to its extent. 

Layer configuration in mapfile might also be wrong since there are 3 files
and 2 different methods to open them ( I believe i tried all the
combinations though)

  _____  

From: Rahkonen Jukka [mailto:Jukka.Rahkonen at mmmtike.fi] 
Sent: Monday, November 12, 2007 9:56 AM
To: Volkan Ulukut; MAPSERVER-USERS at LISTS.UMN.EDU
Subject: Re: [UMN_MAPSERVER-USERS] Georeferencing Mrsid raster file with a
world file

Hi,

I do not know anything about MrSID, but I believe you have some conflict
between MrSID native georeference data and world file and MapServer (or
DGAL) gives priority for the native data.  To check this you could rename
the world file to .BAK or whatever and see if the situation remains the
same. If this is the case three alternatives comes into my mind:

- ask if there is some hidden switch to tell GDAL to give priority for the
world file by sending a question to GDAL-dev mailing list

- try to find some MrSID tool that could be used for changing the embedded
metadata

- try to convert your image to some other format and drop the georeferencing
information totally; then you can use your world file with the new image.

I know that it is possible to tell GDAL to write clean tiff files without
Geotiff tags so tiff could be used with the last approach, but unfortunately
your image is far too big to be uncompressed to tiff.  And probably you
would prefer keeping the image compressen anyway. Are you sure that the
world file extensions and your default projection are correct for the MrSID
image, and not the embedded data?

-Jukka Rahkonen-

  _____  

Lähettäjä: UMN MapServer Users List [mailto:MAPSERVER-USERS at LISTS.UMN.EDU]
Puolesta Volkan Ulukut
Lähetetty: 12. marraskuuta 2007 7:57
Vastaanottaja: MAPSERVER-USERS at LISTS.UMN.EDU
Aihe: [UMN_MAPSERVER-USERS] Georeferencing Mrsid raster file with a world
file

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/e76a4ba3/attachment.html