[Gdal-dev] Confusion with gdalwarp

[R. Piepul]

I seem to have a misunderstanding about the use of gdalwarp and am
hoping that someone can tell me what is wrong.  I am trying to prepare
a GeoTiff image for use in mapserver.

I'm using gdalwarp to project a GeoTiff raster image.  The coordinates
of the input file are in lat long.  After processing the file with
gdalwarp, I used gdalinfo to examine the output file.  It appears that
the coordinate system information is no longer there. Also, the corner
coordinates in the output file are in pixels.

I was expecting gdalinfo to display information about the projected
coordinate system.  Also, I expected the corner coordinates to be
displayed in projected units (meters).  Is this the normal behavior
for gdalwarp?  I'm concerned that I may have a configuration or build
issue.

Thanks for your help,
Rob



$ gdalinfo map.tif
ERROR 1: Undefined symbol "pj_get_def"
Driver: GTiff/GeoTIFF
Size is 5183, 3887
Coordinate System is:
GEOGCS["NAD83",
    DATUM["North_American_Datum_1983",
        SPHEROID["GRS 1980",6378137,298.2572221010002,
            AUTHORITY["EPSG","7019"]],
        AUTHORITY["EPSG","6269"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4269"]]
Origin = (-96.932100,36.761300)
Pixel Size = (0.000278,-0.000278)
Metadata:
  TIFFTAG_SOFTWARE=IMAGINE TIFF Support
Copyright 1991 - 1999 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: etif.c $ $Revision: 1.9.1.2 $ $Date: 2001/12/05 00:33:12Z $
Corner Coordinates:
Upper Left  ( -96.9321000,  36.7613000) ( 96d55'55.56"W, 36d45'40.68"N)
Lower Left  ( -96.9321000,  35.6813695) ( 96d55'55.56"W, 35d40'52.93"N)
Upper Right ( -95.4921000,  36.7613000) ( 95d29'31.56"W, 36d45'40.68"N)
Lower Right ( -95.4921000,  35.6813695) ( 95d29'31.56"W, 35d40'52.93"N)
Center      ( -96.2121000,  36.2213347) ( 96d12'43.56"W, 36d13'16.81"N)
Band 1 Block=5183x1 Type=Byte, ColorInterp=Palette
  Color Table (RGB with 256 entries)
    0: 0,0,0,255
    . . .
    [snipped]


$ gdalwarp -t_srs 'proj=utm zone=14' map.tif temp1.tif
:0...10...20...30...40...50...60...70...80...90...100 - done.


$ gdalinfo temp1.tif
Driver: GTiff/GeoTIFF
Size is 4879, 4536
Coordinate System is `'
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0, 4536.0)
Upper Right ( 4879.0,    0.0)
Lower Right ( 4879.0, 4536.0)
Center      ( 2439.5, 2268.0)
Band 1 Block=4879x1 Type=Byte, ColorInterp=Gray
$


$ gdal-config --version 
1.1.8
$ gdal-config --ogr-enabled
yes
$ gdal-config --formats
gxf gtiff hfa aigrid aaigrid ceos ceos2 iso8211 bsb xpm sdts raw dted
mem jdem envisat elas fit vrt usgsdem l1b nitf bmp hdf4 gif jpeg png
grass