[Gdal-dev] gdalwarp: bad output header?

[Michael Tan]

Hello Frank-

The files are from OpenEV_FW_154.zip

Here's the info you requested:
> proj
Rel. 4.4.5, 9 January 2002
usage: bin\proj [ -beEfiIlormsStTvVwW [args] ] [ +opts[=arg] ] [ files ]

> gdalinfo p174r38_5t870814_nn2.tif
Driver: GTiff/GeoTIFF
Size is 7857, 7777
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.2572235630016,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    AUTHORITY["EPSG","32636"],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",33],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (621670.500000,3625029.000000)
Pixel Size = (28.500000,-28.500000)
Corner Coordinates:
Upper Left  (  621670.500, 3625029.000) ( 34d17'55.89"E, 32d45'23.07"N)
Lower Left  (  621670.500, 3403384.500) ( 34d16'16.43"E, 30d45'25.92"N)
Upper Right (  845595.000, 3625029.000) ( 36d41'11.06"E, 32d42'32.04"N)
Lower Right (  845595.000, 3403384.500) ( 36d36'29.44"E, 30d42'47.65"N)
Center      (  733632.750, 3514206.750) ( 35d27'58.26"E, 31d44'22.16"N)
Band 1 Block=7857x1 Type=Byte, ColorInterp=Gray

-----Original Message-----
From: gdal-dev-admin at remotesensing.org
[mailto:gdal-dev-admin at remotesensing.org]On Behalf Of Frank Warmerdam
Sent: Friday, April 11, 2003 6:04 PM
To: gdal-dev at remotesensing.org
Subject: Re: [Gdal-dev] gdalwarp: bad output header?


Michael Tan wrote:
> Hi-
>
> I am trying to reproject a UTM GeoTIFF to an LatLon GeoTIFF. gdalwarp
> seems to do it, but does not seem to produce a correct header (which I
> need to determine the lat-lon coords for the corners).

Michael,

I am not sure what is going wrong.  This _can_ work since I just tried it
without a problem:


warmerda at gdal2200[18]% gdalwarp -t_srs "+proj=latlong +datum=WGS84"
openev/utm.tif out.tif
Creating output file is that 558P x 466L.
:0...10...20...30...40...50...60...70...80...90...100 - done.
warmerda at gdal2200[19]% gdalinfo out.tif
Driver: GTiff/GeoTIFF
Size is 558, 466
Coordinate System is:
GEOGCS["WGS 84",
     DATUM["WGS_1984",
         SPHEROID["WGS 84",6378137,298.2572235629972,
             AUTHORITY["EPSG","7030"]],
         AUTHORITY["EPSG","6326"]],
     PRIMEM["Greenwich",0],
     UNIT["degree",0.0174532925199433],
     AUTHORITY["EPSG","4326"]]
Origin = (-117.642054,33.903651)
Pixel Size = (0.000597,-0.000597)
Corner Coordinates:
Upper Left  (-117.6420540,  33.9036511) (117d38'31.39"W, 33d54'13.14"N)
Lower Left  (-117.6420540,  33.6252533) (117d38'31.39"W, 33d37'30.91"N)
Upper Right (-117.3086936,  33.9036511) (117d18'31.30"W, 33d54'13.14"N)
Lower Right (-117.3086936,  33.6252533) (117d18'31.30"W, 33d37'30.91"N)
Center      (-117.4753738,  33.7644522) (117d28'31.35"W, 33d45'52.03"N)
Band 1 Block=558x14 Type=Byte, ColorInterp=Gray
warmerda at gdal2200[20]%

What version of GDAL are you using?  What did the coordinate system info
on
the source file look like?  The gdalinfo report on the source should
be sufficient.

Best regards,

--
---------------------------------------+---------------------------------
-----
I set the clouds in motion - turn up   | Frank Warmerdam,
warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Programmer for Rent


_______________________________________________
Gdal-dev mailing list
Gdal-dev at remotesensing.org
http://remotesensing.org/mailman/listinfo/gdal-dev