incorrect result from gdalwarp - NAD83 Albers to Michigan GeoRef

[ctoney]

Results from gdalwarp seem incorrect when reprojecting from USGS North American Albers (NAD83) to Michigan Oblique Mercator (NAD83). The output raster is not in the right place.

Command line:

gdalwarp -t_srs EPSG:3078 -tr 30 30 -srcnodata -9999 -dstnodata -9999 -of HFA evt_1_Alb83.img evt_1_MIGeoRef_gdal.img

Test data sets are attached:
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-

gdalinfo -norat -noct evt_1_alb83.img Driver: HFA/Erdas Imagine Images (.img) Files: evt_1_alb83.img Size is 30, 24 Coordinate System is: PROJCS["NAD_1983_Albers",

GEOGCS["GCS_North_American_1983",

DATUM["North_American_Datum_1983",

SPHEROID["GRS_1980",6378137.0,298.257222101]],

PRIMEM["Greenwich",0.0], UNIT["Degree",0.0174532925199433]],

PROJECTIONAlbers_Conic_Equal_Area, PARAMETER["False_Easting",0.0], PARAMETER["False_Northing",0.0], PARAMETER["longitude_of_center",-96.0], PARAMETER["Standard_Parallel_1",29.5], PARAMETER["Standard_Parallel_2",45.5], PARAMETER["latitude_of_center",23.0], UNIT["Meter",1.0]]

Origin = (633135.000000000000000,2629425.000000000000000)

Pixel Size = (30.000000000000000,-30.000000000000000)

Corner Coordinates:

Upper Left ( 633135.000, 2629425.000) ( 87d46'39.79"W, 46d24'25.53"N)

Lower Left ( 633135.000, 2628705.000) ( 87d46'42.70"W, 46d24'2.24"N)

Upper Right ( 634035.000, 2629425.000) ( 87d45'57.92"W, 46d24'23.00"N)

Lower Right ( 634035.000, 2628705.000) ( 87d46'0.83"W, 46d23'59.71"N)

Center ( 633585.000, 2629065.000) ( 87d46'20.31"W, 46d24'12.62"N)

Band 1 Block=30x24 Type=Int16, ColorInterp=Palette

Description = Band_1 Min=21.000 Max=2481.000 Minimum=21.000, Maximum=2481.000, Mean=2266.315, StdDev=395.461 NoData Value=-9999

gdalinfo -norat -noct evt_1_migeoref_gdal.img Driver: HFA/Erdas Imagine Images (.img) Files: evt_1_migeoref_gdal.img Size is 32, 27 Coordinate System is: PROJCS["NAD83_Michigan_Oblique_Mercator",

GEOGCS["GCS_North_American_1983",

DATUM["North_American_Datum_1983",

SPHEROID["GRS_1980",6378137,298.257222101]],

PRIMEM["Greenwich",0], UNIT["Degree",0.017453292519943295]],

PROJECTIONHotine_Oblique_Mercator, PARAMETER["latitude_of_center",45.30916666666666], PARAMETER["longitude_of_center",-86], PARAMETER["azimuth",337.25556], PARAMETER["scale_factor",0.9996], PARAMETER["false_easting",2546731.496], PARAMETER["false_northing",-4354009.816], UNIT["Meter",1]]

Origin = (2410019.541034979300000,-4230498.913706443300000)

Pixel Size = (30.000000000000000,-30.000000000000000)

Corner Coordinates:

Upper Left ( 2410019.541,-4230498.914) ( 87d46'43.45"W, 46d24'25.47"N)

Lower Left ( 2410019.541,-4231308.914) ( 87d46'42.60"W, 46d23'59.24"N)

Upper Right ( 2410979.541,-4230498.914) ( 87d45'58.50"W, 46d24'26.17"N)

Lower Right ( 2410979.541,-4231308.914) ( 87d45'57.66"W, 46d23'59.93"N)

Center ( 2410499.541,-4230903.914) ( 87d46'20.55"W, 46d24'12.70"N)

Band 1 Block=64x64 Type=Int16, ColorInterp=Palette

Description = Layer_1 Min=21.000 Max=2481.000 Minimum=21.000, Maximum=2481.000, Mean=2266.046, StdDev=395.452 NoData Value=-9999

gdalinfo -norat -noct evt_1_migeoref_imagine.img Driver: HFA/Erdas Imagine Images (.img) Files: evt_1_migeoref_imagine.img Size is 32, 27 Coordinate System is: PROJCS["Oblique Mercator (Hotine)",

GEOGCS["NAD83",

DATUM["North_American_Datum_1983",

SPHEROID["GRS 1980",6378137,298.257222101,

AUTHORITY["EPSG","7019"]],

TOWGS84[0,0,0,0,0,0,0], AUTHORITY["EPSG","6269"]],

PRIMEM["Greenwich",0,

AUTHORITY["EPSG","8901"]],

UNIT["degree",0.0174532925199433,

AUTHORITY["EPSG","9108"]],

AUTHORITY["EPSG","4269"]],

PROJECTIONHotine_Oblique_Mercator, PARAMETER["latitude_of_center",45.30916666666666], PARAMETER["longitude_of_center",-86], PARAMETER["azimuth",337.25556], PARAMETER["rectified_grid_angle",0], PARAMETER["scale_factor",0.9996], PARAMETER["false_easting",2546731.496], PARAMETER["false_northing",-4354009.816], UNIT["meters",1]]

Origin = (363137.000000000000000,652110.000000000000000)

Pixel Size = (30.000000000000000,-30.000000000000000)

Corner Coordinates:

Upper Left ( 363137.000, 652110.000) (178d20'9.95"E, 72d16'45.90"N)

Lower Left ( 363137.000, 651300.000) (178d21'35.77"E, 72d16'45.30"N)

Upper Right ( 364097.000, 652110.000) (178d20'12.26"E, 72d17'16.88"N)

Lower Right ( 364097.000, 651300.000) (178d21'38.13"E, 72d17'16.28"N)

Center ( 363617.000, 651705.000) (178d20'54.03"E, 72d17'1.09"N)

Band 1 Block=32x27 Type=Int16, ColorInterp=Palette

Description = Band_1 Min=21.000 Max=2481.000 Minimum=21.000, Maximum=2481.000, Mean=2265.965, StdDev=395.718 NoData Value=-9999

[Even Rouault]

Might be another instance of a more general problem when dealing with Hotine Oblique Mercator, as in #3311 or #3038

[Even Rouault]

This ticket has been automatically closed because Trac is no longer used for GDAL bug tracking, since the project has migrated to GitHub. If you believe this ticket is still valid, you may file it to ​https://github.com/OSGeo/gdal/issues if it is not already reported there.