[gdal-dev] incorrect result from gdalwarp to Michigan GeoRef

[Frank Warmerdam]

Chris Toney wrote:
> Hi,
> gdalwarp seems to give an incorrect result when reprojecting from USGS North American Albers (NAD83) to Michigan Oblique Mercator (NAD83). The output raster is not in the right place.
>  
> Tried with gdal 1.5.1, FWTools 2.2.8 and FWTools 2.3.0. Command line used:
>  
> gdalwarp -t_srs EPSG:3078 -tr 30 30 -srcnodata -9999 -dstnodata -9999 -of HFA evt_1_Alb83.img evt_1_MIGeoRef_gdal.img
>  
> The EPSG definition in pcs.csv is correct as far as I can tell.
>  
> Test data sets are available at:
> http://www.spatialsimulation.org/test/gdalwarp_MI_georef.zip
> evt_1_Alb83.img is the source raster
> evt1_MIGeoRef_gdal.img is the output from gdalwarp
> evt1_MIGeoRef_imagine.img is the output from Erdas Imagine based on specifying EPSG 3078. This raster is in the correct location.
>  
> gdalinfo for these rasters is copied below.
>  
> Am I missing anything with the gdalwarp usage?

Chris,

I have not reviewed this in detail, but I would suggest filing a ticket on
the issue.  It isn't immediately obvious to me what is going wrong.

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