[gdal-dev] Re: [GDAL_TRANSLATE] Gray scale

Solimyr giannicristian at msn.com
Thu Jun 30 10:44:37 EDT 2011


I already tried this step:

$ gdal_translate -of ENVI -ot Byte -scale  File.tif Prova.raw
Input file size is 1828, 1646
0...10...20...30...40...50...60...70...80...90...100 - done.

$ gdalinfo -stats prova.raw
Driver: ENVI/ENVI .hdr Labelled
Files: prova.raw
       prova.raw.aux.xml
       prova.hdr
Size is 1828, 1646
Coordinate System is:
PROJCS["UTM Zone 47, Southern Hemisphere",
    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"]],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",99],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",10000000],
    UNIT["Meter",1]]
Origin = (650301.875000000000000,9895384.375000000000000)
Pixel Size = (1.250000000000000,-1.250000000000000)
Metadata:
  Band_1=Band 1
  TIFFTAG_XRESOLUTION=1
  TIFFTAG_YRESOLUTION=1
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  650301.875, 9895384.375) (100d21' 2.82"E,  0d56'46.41"S)
Lower Left  (  650301.875, 9893326.875) (100d21' 2.84"E,  0d57'53.40"S)
Upper Right (  652586.875, 9895384.375) (100d22'16.73"E,  0d56'46.38"S)
Lower Right (  652586.875, 9893326.875) (100d22'16.76"E,  0d57'53.37"S)
Center      (  651444.375, 9894355.625) (100d21'39.79"E,  0d57'19.89"S)
Band 1 Block=1828x1 Type=Byte, ColorInterp=Undefined
  Description = Band 1
  Minimum=0.000, Maximum=255.000, Mean=2.364, StdDev=3.013
  Metadata:
    STATISTICS_MINIMUM=0
    STATISTICS_MAXIMUM=255
    STATISTICS_MEAN=2.3638978918458
    STATISTICS_STDDEV=3.0127939567847

And the results from ENVI is:

$ gdalinfo -stats fromENVI.raw
Driver: ENVI/ENVI .hdr Labelled
Files: fromENVI.raw
       fromENVI.raw.hdr
Size is 1828, 1646
Coordinate System is:
PROJCS["UTM Zone 47, Southern Hemisphere",
    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"]],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",99],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",10000000],
    UNIT["Meter",1]]
Origin = (650301.875000000000000,9895384.375000000000000)
Pixel Size = (1.250000000000000,-1.250000000000000)
Metadata:
  Band_1=Gray Scale (Band 1:File.tif)
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  650301.875, 9895384.375) (100d21' 2.82"E,  0d56'46.41"S)
Lower Left  (  650301.875, 9893326.875) (100d21' 2.84"E,  0d57'53.40"S)
Upper Right (  652586.875, 9895384.375) (100d22'16.73"E,  0d56'46.38"S)
Lower Right (  652586.875, 9893326.875) (100d22'16.76"E,  0d57'53.37"S)
Center      (  651444.375, 9894355.625) (100d21'39.79"E,  0d57'19.89"S)
Band 1 Block=1828x1 Type=Byte, ColorInterp=Undefined
  Description = Gray Scale (Band 1:File.tif)
  Minimum=0.000, Maximum=255.000, Mean=53.375, StdDev=48.454
  Metadata:
    STATISTICS_MINIMUM=0
    STATISTICS_MAXIMUM=255
    STATISTICS_MEAN=53.374992023631
    STATISTICS_STDDEV=48.454042778451

Looking at the histogram seems that ENVI do a conversion using a kind of
equalize scale. Do you know where I can find the method that ENVI use? I
have searched a lot without solution and seems that I'm the only one with
this problem! Very strange!

--
View this message in context: http://osgeo-org.1803224.n2.nabble.com/GDAL-TRANSLATE-Gray-scale-tp6533237p6533558.html
Sent from the GDAL - Dev mailing list archive at Nabble.com.


More information about the gdal-dev mailing list