[gdal-dev] gdal2tiles.py and ProjNatOriginLongGeoKey

[Eduard Bondarenko]

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