GDAL Wish: Projection information may print UTM zone

[Markus Neteler]

Frank,

may I ask you to add the UTM *zone* information to the printout
of the projection information?
E.g. 

gdalinfo /cdrom/L71173073_07320011008_HTM.FST
Driver: FAST/EOSAT FAST Format
Size is 3707, 3457
Coordinate System is:
PROJCS["UTM",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            TOWGS84[0,0,0,0,0,0,0],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9108"]],
        AXIS["Lat","NORTH"],
        AXIS["Long","EAST"],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",27],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",10000000]]
Origin = (234345.000000,8025915.000000)
Pixel Size = (60.000000,-60.000000)
Metadata:
  BIAS1=       0.000000000000000

  GAIN1=       0.066823534667492

  BIAS2=       3.199999809265137

  GAIN2=       0.037058822810650

Corner Coordinates:
Upper Left  (  234345.000, 8025915.000) ( 24d29'36.59"E, 17d50'18.15"S)
Lower Left  (  234345.000, 7818495.000) ( 24d27'56.30"E, 19d42'40.30"S)
Upper Right (  456765.000, 8025915.000) ( 26d35'30.95"E, 17d51'14.56"S)
Lower Right (  456765.000, 7818495.000) ( 26d35'14.60"E, 19d43'43.08"S)
Center      (  345555.000, 7922205.000) ( 25d32'4.48"E, 18d47'9.92"S)
Band 1 Block=3707x1 Type=Byte, ColorInterp=Undefined
Band 2 Block=3707x1 Type=Byte, ColorInterp=Undefined

From looking into a UTM map I found that central meridian 27E
refers to UTM 35. Perhaps you could print the EPSG code for
the projection found and/or something like

PROJCS["UTM zone 35S",

Thanks in advance,
 merry christmas

   Markus

[warmerdam]

Andrey,

Could you try removing the SetProjCS() call in fastdataset.cpp?  If
you do this, I think the SetUTM() call will set a more detailed name
for the projcs including the UTM zone. 

[neteler@…]

Hi,

the same *may* apply also to these formats?

find . -type f -name "*.cpp" -exec grep -l SetProjCS {} \;
./frmts/hdf4/hdf4imagedataset.cpp
./frmts/hfa/hfadataset.cpp

Best regards

 Markus

[dron]

Hi, Markus,

I did it for FAST and HDF formats, now UTM zone reported correctly. I can't just
remove SetProjCS() call from Erdas Imagine driver (hfadataset.cpp), it will
broke the driver logic.

            Happy New Year!
            Andrey

[neteler@…]

Hi Andrey,

... works perfectly (have tested FAST only).

Thanks and happy new year,

 Markus

[dron]

Ok, considered as fixed.