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