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:[EMAIL PROTECTED] Puolesta Volkan Ulukut Lähetetty: 12. marraskuuta 2007 7:57 Vastaanottaja: MAPSERVER-USERS@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
