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

[rouault]

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