<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">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. <div><br></div><div>I downloaded some aerial imagery from here:<br><a href="http://www.michigan.gov/dnr/0,1607,7-153-10371_14546-30211--,00.html">http://www.michigan.gov/dnr/0,1607,7-153-10371_14546-30211--,00.html</a><br><br>I want to import it into a UTM zone 17 project. So I did the following:<br>gdalwarp -t_srs '+proj=utm +zone=17 +datum=WGS84' input.sid output_utm.sid<br><br>gdapwarp responds with an error, but then continues with the conversion:<br>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 '.<br>lat_0 = 0 or 90 or alpha = 90<br>0...10...20...30...40...50...60...70...80...90...100 - done.<br><div><br></div><div>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. </div><div><br></div><div>Any suggestions for what I might be doing wrong? </div><div><br></div><div><br></div><div>If I do gdalinfo on the original file here is what I get:</div><div><br></div><div>Driver: MrSID/Multi-resolution Seamless Image Database (MrSID)<br>Files: oakgrove_sw.sid<br> oakgrove_sw.sdw<br> oakgrove_sw.sid.aux.xml<br>Size is 6415, 8061<br>Coordinate System is:<br>PROJCS["IMAGINE GeoTIFF Support<br>Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved<br>@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $<br>Projection = Oblique Mercator (Hotine)|IMAGINE",<br> GEOGCS["NAD83",<br> DATUM["North_American_Datum_1983",<br> SPHEROID["GRS 1980",6378137,298.2572221010002,<br> AUTHORITY["EPSG","7019"]],<br> AUTHORITY["EPSG","6269"]],<br> PRIMEM["Greenwich",0],<br> UNIT["degree",0.0174532925199433],<br> AUTHORITY["EPSG","4269"]],<br> PROJECTION["Hotine_Oblique_Mercator"],<br> PARAMETER["latitude_of_center",0],<br> PARAMETER["longitude_of_center",0],<br> PARAMETER["azimuth",0],<br> PARAMETER["rectified_grid_angle",90],<br> PARAMETER["scale_factor",1],<br> PARAMETER["false_easting",0],<br> PARAMETER["false_northing",0],<br> UNIT["metre",1,<br> AUTHORITY["EPSG","9001"]]]<br>Origin = (663144.500000000000000,240035.500000000000000)<br>Pixel Size = (1.000000000000000,-1.000000000000000)<br>Metadata:<br> IMAGE__INPUT_NAME=K:\blk18\oakgrove_sw.tif<br> IMAGE__INPUT_FILE_SIZE=155199578.000000<br> GEOTIFF_NUM__1024__GTModelTypeGeoKey=1<br> GEOTIFF_CHAR__GTModelTypeGeoKey=ModelTypeProjected<br> GEOTIFF_NUM__1025__GTRasterTypeGeoKey=1<br> GEOTIFF_CHAR__GTRasterTypeGeoKey=RasterPixelIsArea<br> GEOTIFF_NUM__1026__GTCitationGeoKey=IMAGINE GeoTIFF Support<br>Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved<br>@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $<br>Projection = Oblique Mercator (Hotine)|IMAGINE GeoTIFF Support<br>Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved<br>@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $<br>Projection Name = Oblique Mercator (Hotine)<br>Units = meters<br>GeoTIFF Units = meters|<br> GEOTIFF_NUM__2048__GeographicTypeGeoKey=4269<br> GEOTIFF_CHAR__GeographicTypeGeoKey=GCS_NAD83<br> GEOTIFF_NUM__2060__GeogAzimuthUnitsGeoKey=9102<br> GEOTIFF_CHAR__GeogAzimuthUnitsGeoKey=Angular_Degree<br> GEOTIFF_NUM__3072__ProjectedCSTypeGeoKey=32767<br> GEOTIFF_CHAR__ProjectedCSTypeGeoKey=User-Defined<br> GEOTIFF_NUM__3073__PCSCitationGeoKey=IMAGINE GeoTIFF Support<br>Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved<br>@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $<br>Projection = Oblique Mercator (Hotine)|IMAGINE GeoTIFF Support<br>Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved<br>@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $<br>Projection Name = Oblique Mercator (Hotine)<br>Units = meters<br>GeoTIFF Units = meters|<br> GEOTIFF_NUM__3074__ProjectionGeoKey=32767<br> GEOTIFF_CHAR__ProjectionGeoKey=User-Defined<br> GEOTIFF_NUM__3075__ProjCoordTransGeoKey=3<br> GEOTIFF_CHAR__ProjCoordTransGeoKey=CT_ObliqueMercator<br> GEOTIFF_NUM__3076__ProjLinearUnitsGeoKey=9001<br> GEOTIFF_CHAR__ProjLinearUnitsGeoKey=Linear_Meter<br> GEOTIFF_NUM__3088__ProjCenterLongGeoKey=0.999600,45.309167,2546731.496000,-4354009.816000,337.255560,-86.000000<br> GEOTIFF_NUM__3089__ProjCenterLatGeoKey=0.999600,45.309167,2546731.496000,-4354009.816000,337.255560,-86.000000<br> GEOTIFF_NUM__3090__ProjCenterEastingGeoKey=0.999600,45.309167,2546731.496000,-4354009.816000,337.255560,-86.000000<br> GEOTIFF_NUM__3091__ProjCenterNorthingGeoKey=0.999600,45.309167,2546731.496000,-4354009.816000,337.255560,-86.000000<br> GEOTIFF_NUM__3093__ProjScaleAtCenterGeoKey=0.999600,45.309167,2546731.496000,-4354009.816000,337.255560,-86.000000<br> GEOTIFF_NUM__3094__ProjAzimuthAngleGeoKey=0.999600,45.309167,2546731.496000,-4354009.816000,337.255560,-86.000000<br> IMAGE__Z_RESOLUTION=0.000000<br> GEO__ModelTypeGeoKey=1<br> GEO__ProjectedCSTypeGeoKey=32767<br> GEO__PCSCitationGeoKey=IMAGINE GeoTIFF Support<br>Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved<br>@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $<br>Projection = Oblique Mercator (Hotine)|IMAGINE GeoTIFF Support<br>Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved<br>@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $<br>Projection Name = Oblique Mercator (Hotine)<br>Units = meters<br>GeoTIFF Units = meters|<br> IMAGE__Z_ORIGIN=0.000000<br> IMAGE__COMPRESSION_VERSION=1,6,1<br> IMAGE__TARGET_COMPRESSION_RATIO=20.000000<br> IMAGE__COMPRESSION_NLEV=4<br> IMAGE__COMPRESSION_WEIGHT=4.000000<br> IMAGE__COMPRESSION_GAMMA=2.000000<br> IMAGE__COMPRESSION_BLOCK_SIZE=512<br> IMAGE__CREATION_DATE=Wed Mar 10 17:11:37 2004<br><br> VERSION=MG2<br> IMAGE__BITS_PER_SAMPLE=8<br> IMAGE__WIDTH=6415<br> IMAGE__HEIGHT=8061<br> IMAGE__COLOR_SCHEME=0<br> IMAGE__DATA_TYPE=0<br> IMAGE__XY_ORIGIN=663145.000000,240035.000000<br> IMAGE__X_RESOLUTION=1.000000<br> IMAGE__Y_RESOLUTION=1.000000<br>Image Structure Metadata:<br> INTERLEAVE=PIXEL<br>Corner Coordinates:<br>Upper Left ( 663144.500, 240035.500) <br>Lower Left ( 663144.500, 231974.500) <br>Upper Right ( 669559.500, 240035.500) <br>Lower Right ( 669559.500, 231974.500) <br>Center ( 666352.000, 236005.000) <br>Band 1 Block=1024x128 Type=Byte, ColorInterp=Red<br> Minimum=0.000, Maximum=255.000, Mean=97.677, StdDev=49.163<br> Overviews: 3208x4031, 1604x2016, 802x1008, 401x504<br>Band 2 Block=1024x128 Type=Byte, ColorInterp=Green<br> Minimum=0.000, Maximum=255.000, Mean=97.555, StdDev=49.371<br> Overviews: 3208x4031, 1604x2016, 802x1008, 401x504<br>Band 3 Block=1024x128 Type=Byte, ColorInterp=Blue<br> Minimum=0.000, Maximum=255.000, Mean=97.402, StdDev=49.024<br> Overviews: 3208x4031, 1604x2016, 802x1008, 401x504<br></div><div><div><br class="webkit-block-placeholder"></div><div>
<div>Thanks much!</div><div><br></div><div>--Adam</div><div><br class="khtml-block-placeholder"></div><div><br class="khtml-block-placeholder"></div>
</div>
<br></div></div></body></html>