Hello gdal,
I'm trying to convert GeoTIFF to TMS tiles using gdal2tiles.py script.
listgeo:
Geotiff_Information:
Version: 1
Key_Revision: 1.0
Tagged_Information:
ModelTiepointTag (2,3):
0 0 0
600449 4934034 0
ModelPixelScaleTag (1,3):
0.5 0.5 0
End_Of_Tags.
Keyed_Information:
GTModelTypeGeoKey (Short,1): ModelTypeProjected
GTRasterTypeGeoKey (Short,1): RasterPixelIsArea
GTCitationGeoKey (Ascii,257): "IMAGINE GeoTIFF Support\nCopyright
1991 - 2006 by Leica Geosystems Geospatial Imaging, LLC. All Rights
Reserved\n@(#)$RCSfile: egtf.c $ IMAGINE 9.1 $Revision: 22.0 $ $Date:
2006/05/24 17:16:00 EST $\nProjection Name = UTM\nUnits =
meters\nGeoTIFF Units = meters"
ProjectedCSTypeGeoKey (Short,1): PCS_WGS84_UTM_zone_36N
PCSCitationGeoKey (Ascii,244): "IMAGINE GeoTIFF
Support\nCopyright 1991 - 2006 by Leica Geosystems Geospatial Imaging,
LLC. All Rights Reserved\n@(#)$RCSfile: egtf.c $ IMAGINE 9.1 $Revision:
22.0 $ $Date: 2006/05/24 17:16:00 EST $\nUTM Zone 36N\nEllipsoid = WGS
84\nDatum = WGS 84"
ProjLinearUnitsGeoKey (Short,1): Linear_Meter
End_Of_Keys.
End_Of_Geotiff.
PCS = 32636 (WGS 84 / UTM zone 36N)
Projection = 16036 (UTM zone 36N)
Projection Method: CT_TransverseMercator
ProjNatOriginLatGeoKey: 0.000000 ( 0d 0' 0.00"N)
ProjNatOriginLongGeoKey: 33.000000 ( 33d 0' 0.00"E)
ProjScaleAtNatOriginGeoKey: 0.999600
ProjFalseEastingGeoKey: 500000.000000 m
ProjFalseNorthingGeoKey: 0.000000 m
GCS: 4326/WGS 84
Datum: 6326/World Geodetic System 1984
Ellipsoid: 7030/WGS 84 (6378137.00,6356752.31)
Prime Meridian: 8901/Greenwich (0.000000/ 0d 0' 0.00"E)
Projection Linear Units: 9001/metre (1.000000m)
Corner Coordinates:
Upper Left ( 600449.000, 4934034.000) ( 34d15'52.85"E, 44d33' 9.50"N)
Lower Left ( 600449.000, 4931945.500) ( 34d15'51.39"E, 44d32' 1.83"N)
Upper Right ( 602409.500, 4934034.000) ( 34d17'21.69"E, 44d33' 8.51"N)
Lower Right ( 602409.500, 4931945.500) ( 34d17'20.20"E, 44d32' 0.84"N)
Center ( 601429.250, 4932989.750) ( 34d16'36.53"E, 44d32'35.17"N)
I use bash script to convert:
#!/bin/sh
input=1.tif
gdal_translate -of VRT -a_srs EPSG:900913 -gcp 0 0 600449 4934034 -gcp
3921 4177 602409.5 4931945.5 -gcp 3921 0 602409.5 4934034 $input 2.vrt
rm 3.vrt
rm -rf sat
gdalwarp -of VRT -s_srs EPSG:900913 -t_srs EPSG:4326 2.vrt 3.vrt
./gdal2tiles.py -v -z 4-18 --s_srs EPSG:4326 3.vrt sat
The problem is tiles are generated but it looks like they lie far to the
west from place they should be.
I guess ProjNatOriginLongGeoKey tells about this shift but I can't find
out what it means and how use it?
Thanks in advance,
Eduard
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
