[Gdal-dev] GeoTIFF Georeferencing Problem

[Meek, Allan (Contractor)]

Hi,

I'm trying to work with a set of GeoTIFFs that are giving me
georeferencing headaches.

One example is available at:

http://www7333.nrlssc.navy.mil/browse/lvl4/modis/5.2/METOC-Halifax/7dayl
p/2007/aug/MODAM20072072007213.L4_LP_NOAA_HFX_sst.tiff

Gdalinfo produces the following output (I've removed some irrelevant
lines):

Driver: GTiff/GeoTIFF
Size is 1749, 869
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["unnamed",
        DATUM["unknown",
            SPHEROID["unnamed",6378137,298.257223560493]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Mercator_1SP"],
    PARAMETER["latitude_of_origin",39],
    PARAMETER["central_meridian",-59.5],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (1080355.008799013900000,4727540.859266426400000)
Pixel Size = (-1237.520055898069100,-1237.520055898069400)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_DOCUMENTNAME=MODAM20072072007213.L4_LP_NOAA_HFX_sst.tiff
  TIFFTAG_IMAGEDESCRIPTION=
createTime = "Wed Aug  1 20:07:38 2007"
createSoftware = "APS v3.5.3"
createPlatform = "i686-pc-linux-gnu"
createAgency = "Naval Research Laboratory, Stennis Space Center"
mapProjectionSystem = "NRL(USGS)"
mapProjection = "METOC-Halifax"
mapUpperLeft = [48.0003,-72]
mapUpperRight = [48.0003,-47.0286]
mapLowerLeft = [39,-72]
mapLowerRight = [39,-47.0286]
  TIFFTAG_SOFTWARE=imgConvertTIFF v3.5.3 (compiled Feb  5 2007, 12:26:26
for 
i686-pc-linux-gnu)
  TIFFTAG_DATETIME=2007:08:01 15:07:43
  TIFFTAG_XRESOLUTION=72
  TIFFTAG_YRESOLUTION=72
  TIFFTAG_RESOLUTIONUNIT=1 (unitless)
Corner Coordinates:
Upper Left  ( 1080355.009, 4727540.859) ( 49d47'42.02"W, 39d13'45.43"N)
Lower Left  ( 1080355.009, 3652135.931) ( 49d47'42.02"W, 31d19'11.08"N)
Upper Right (-1084067.569, 4727540.859) ( 69d14'18.04"W, 39d13'45.43"N)
Lower Right (-1084067.569, 3652135.931) ( 69d14'18.04"W, 31d19'11.08"N)
Center      (   -1856.280, 4189838.395) ( 59d31'0.03"W, 35d22'19.83"N)
Band 1 Block=1749x1 Type=Float32, ColorInterp=Gray

The reported corner coordinates are wrong -- you can see the actual
coordinates in mapUpperLeft, mapUpperRight, etc. Also, GEOGCS, DATUM,
and SPHEROID are blank.

When I open the file in ENVI, I get a window with the message:

GeoTIFF Datum: <blank>
  Unrecognized GeoTIFF datum name.
  Please select ENVI datum name from list.

Choosing WGS-84 loads the image properly.

Is there anything I can do to fix this -- to have GDAL report the
correct corner coordinates?

Thanks,
Allan