[Gdal-dev] Geotiff metadata conversion?

[Don Tidrow]

I have a set of geotiff images that I need to be able to reproject into 
lat/long space in order to slice them up into tiles of specific dimensions 
(in arcseconds).  The output of gdalinfo for the files looks like this:

Driver: GTiff/GeoTIFF
Size is 7876, 12796
Coordinate System is `'
GCP Projection = GEOGCS["unnamed",DATUM["WGS_1984",SPHEROID["WGS 
84",6378137,298.2572235629972,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]]
GCP[  0]: Id=1, Info=
           (0,0) -> (-76.6208,37.1811,0)
GCP[  1]: Id=2, Info=
           (7876,0) -> (-76.5321,37.1823,0)
GCP[  2]: Id=3, Info=
           (0,12796) -> (-76.6184,37.0658,0)
Metadata:
   TIFFTAG_SOFTWARE=Sensor Systems RemoteView
   TIFFTAG_DATETIME=2004:00:14 11:44:23
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,12796.0)
Upper Right ( 7876.0,    0.0)
Lower Right ( 7876.0,12796.0)
Center      ( 3938.0, 6398.0)
Band 1 Block=1024x1024 Type=Byte, ColorInterp=Gray


I have other imagery, downloaded from the Landsat repository, that works 
fine in my image-processing program.  The output of gdalinfo for these 
files looks like this:

Driver: GTiff/GeoTIFF
Size is 7882, 7517
Coordinate System is:
PROJCS["WGS 84 / UTM zone 38N",
     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"]],
     PROJECTION["Transverse_Mercator"],
     PARAMETER["latitude_of_origin",0],
     PARAMETER["central_meridian",45],
     PARAMETER["scale_factor",0.9996],
     PARAMETER["false_easting",500000],
     PARAMETER["false_northing",0],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]],
     AUTHORITY["EPSG","32638"]]
Origin = (390022.500000,3779841.000000)
Pixel Size = (28.500000,-28.500000)
Metadata:
   TIFFTAG_XRESOLUTION=72
   TIFFTAG_YRESOLUTION=72
   TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
Corner Coordinates:
Upper Left  (  390022.500, 3779841.000) ( 43d48'25.11"E, 34d 9'13.32"N)
Lower Left  (  390022.500, 3565606.500) ( 43d49'58.39"E, 32d13'18.10"N)
Upper Right (  614659.500, 3779841.000) ( 46d14'37.70"E, 34d 9'11.50"N)
Lower Right (  614659.500, 3565606.500) ( 46d13'0.45"E, 32d13'16.41"N)
Center      (  502341.000, 3672723.750) ( 45d 1'30.41"E, 33d11'36.11"N)
Band 1 Block=7882x1 Type=Byte, ColorInterp=Gray
   Overviews: 986x940


Is there a gdal utility/function(s) that will convert from the first format 
into the second, or is there simply not enough metadata in the first 
geotiff to convert it?

Don

-- 

+----------------------------------------+
|  Don Tidrow                            |
|  Visualization Software Engineer       |
|  Object Sciences Corporation           |
|  ph: 703-253-1094                      |
+----------------------------------------+