Georeferencing Mrsid raster file with a world file

[Rahkonen Jukka]

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