[gdal-dev] gdapwarp problem

Adam Dershowitz, Ph.D., P.E. adershowitz at exponent.com
Fri Feb 11 13:20:41 EST 2011 

I am having a problem trying to use gdalwarp.  Perhaps I am just missing something, but I figured I would see if someone can help.  

I downloaded some aerial imagery from here:
http://www.michigan.gov/dnr/0,1607,7-153-10371_14546-30211--,00.html

I want to import it into a UTM zone 17 project.  So I did the following:
gdalwarp -t_srs '+proj=utm +zone=17 +datum=WGS84' input.sid output_utm.sid

gdapwarp responds with an error, but then continues with the conversion:
ERROR 6: Failed to initialize PROJ.4 with `+proj=omerc +lat_0=0 +lonc=0 +alpha=0 +k=1 +x_0=0 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs '.
lat_0 = 0 or 90 or alpha = 90
0...10...20...30...40...50...60...70...80...90...100 - done.

It does complete, but the area of interest is in the range of UTM (254500,4734000) and the output of gdalwarp is around 663144,240035.  

Any suggestions for what I might be doing wrong?  


If I do gdalinfo on the original file here is what I get:

Driver: MrSID/Multi-resolution Seamless Image Database (MrSID)
Files: oakgrove_sw.sid
      oakgrove_sw.sdw
      oakgrove_sw.sid.aux.xml
Size is 6415, 8061
Coordinate System is:
PROJCS["IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection = Oblique Mercator (Hotine)|IMAGINE",
   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"]],
   PROJECTION["Hotine_Oblique_Mercator"],
   PARAMETER["latitude_of_center",0],
   PARAMETER["longitude_of_center",0],
   PARAMETER["azimuth",0],
   PARAMETER["rectified_grid_angle",90],
   PARAMETER["scale_factor",1],
   PARAMETER["false_easting",0],
   PARAMETER["false_northing",0],
   UNIT["metre",1,
       AUTHORITY["EPSG","9001"]]]
Origin = (663144.500000000000000,240035.500000000000000)
Pixel Size = (1.000000000000000,-1.000000000000000)
Metadata:
 IMAGE__INPUT_NAME=K:\blk18\oakgrove_sw.tif
 IMAGE__INPUT_FILE_SIZE=155199578.000000
 GEOTIFF_NUM__1024__GTModelTypeGeoKey=1
 GEOTIFF_CHAR__GTModelTypeGeoKey=ModelTypeProjected
 GEOTIFF_NUM__1025__GTRasterTypeGeoKey=1
 GEOTIFF_CHAR__GTRasterTypeGeoKey=RasterPixelIsArea
 GEOTIFF_NUM__1026__GTCitationGeoKey=IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection = Oblique Mercator (Hotine)|IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection Name = Oblique Mercator (Hotine)
Units = meters
GeoTIFF Units = meters|
 GEOTIFF_NUM__2048__GeographicTypeGeoKey=4269
 GEOTIFF_CHAR__GeographicTypeGeoKey=GCS_NAD83
 GEOTIFF_NUM__2060__GeogAzimuthUnitsGeoKey=9102
 GEOTIFF_CHAR__GeogAzimuthUnitsGeoKey=Angular_Degree
 GEOTIFF_NUM__3072__ProjectedCSTypeGeoKey=32767
 GEOTIFF_CHAR__ProjectedCSTypeGeoKey=User-Defined
 GEOTIFF_NUM__3073__PCSCitationGeoKey=IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection = Oblique Mercator (Hotine)|IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection Name = Oblique Mercator (Hotine)
Units = meters
GeoTIFF Units = meters|
 GEOTIFF_NUM__3074__ProjectionGeoKey=32767
 GEOTIFF_CHAR__ProjectionGeoKey=User-Defined
 GEOTIFF_NUM__3075__ProjCoordTransGeoKey=3
 GEOTIFF_CHAR__ProjCoordTransGeoKey=CT_ObliqueMercator
 GEOTIFF_NUM__3076__ProjLinearUnitsGeoKey=9001
 GEOTIFF_CHAR__ProjLinearUnitsGeoKey=Linear_Meter
 GEOTIFF_NUM__3088__ProjCenterLongGeoKey=0.999600,45.309167,2546731.496000,-4354009.816000,337.255560,-86.000000
 GEOTIFF_NUM__3089__ProjCenterLatGeoKey=0.999600,45.309167,2546731.496000,-4354009.816000,337.255560,-86.000000
 GEOTIFF_NUM__3090__ProjCenterEastingGeoKey=0.999600,45.309167,2546731.496000,-4354009.816000,337.255560,-86.000000
 GEOTIFF_NUM__3091__ProjCenterNorthingGeoKey=0.999600,45.309167,2546731.496000,-4354009.816000,337.255560,-86.000000
 GEOTIFF_NUM__3093__ProjScaleAtCenterGeoKey=0.999600,45.309167,2546731.496000,-4354009.816000,337.255560,-86.000000
 GEOTIFF_NUM__3094__ProjAzimuthAngleGeoKey=0.999600,45.309167,2546731.496000,-4354009.816000,337.255560,-86.000000
 IMAGE__Z_RESOLUTION=0.000000
 GEO__ModelTypeGeoKey=1
 GEO__ProjectedCSTypeGeoKey=32767
 GEO__PCSCitationGeoKey=IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection = Oblique Mercator (Hotine)|IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection Name = Oblique Mercator (Hotine)
Units = meters
GeoTIFF Units = meters|
 IMAGE__Z_ORIGIN=0.000000
 IMAGE__COMPRESSION_VERSION=1,6,1
 IMAGE__TARGET_COMPRESSION_RATIO=20.000000
 IMAGE__COMPRESSION_NLEV=4
 IMAGE__COMPRESSION_WEIGHT=4.000000
 IMAGE__COMPRESSION_GAMMA=2.000000
 IMAGE__COMPRESSION_BLOCK_SIZE=512
 IMAGE__CREATION_DATE=Wed Mar 10 17:11:37 2004

 VERSION=MG2
 IMAGE__BITS_PER_SAMPLE=8
 IMAGE__WIDTH=6415
 IMAGE__HEIGHT=8061
 IMAGE__COLOR_SCHEME=0
 IMAGE__DATA_TYPE=0
 IMAGE__XY_ORIGIN=663145.000000,240035.000000
 IMAGE__X_RESOLUTION=1.000000
 IMAGE__Y_RESOLUTION=1.000000
Image Structure Metadata:
 INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  663144.500,  240035.500) 
Lower Left  (  663144.500,  231974.500) 
Upper Right (  669559.500,  240035.500) 
Lower Right (  669559.500,  231974.500) 
Center      (  666352.000,  236005.000) 
Band 1 Block=1024x128 Type=Byte, ColorInterp=Red
 Minimum=0.000, Maximum=255.000, Mean=97.677, StdDev=49.163
 Overviews: 3208x4031, 1604x2016, 802x1008, 401x504
Band 2 Block=1024x128 Type=Byte, ColorInterp=Green
 Minimum=0.000, Maximum=255.000, Mean=97.555, StdDev=49.371
 Overviews: 3208x4031, 1604x2016, 802x1008, 401x504
Band 3 Block=1024x128 Type=Byte, ColorInterp=Blue
 Minimum=0.000, Maximum=255.000, Mean=97.402, StdDev=49.024
 Overviews: 3208x4031, 1604x2016, 802x1008, 401x504

Thanks much!

--Adam



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20110211/6c7166e6/attachment.html

More information about the gdal-dev mailing list