Opened 14 years ago

Closed 5 years ago

#3038 closed defect (wontfix)

incorrect result from gdalwarp - NAD83 Albers to Michigan GeoRef

Reported by: ctoney Owned by: warmerdam
Priority: normal Milestone: closed_because_of_github_migration
Component: default Version: 1.6.1
Severity: normal Keywords:
Cc: Kyle Shannon


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",




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",




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)",



SPHEROID["GRS 1980",6378137,298.257222101,


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






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

Attachments (1) (12.2 KB ) - added by ctoney 14 years ago.

Download all attachments as: .zip

Change History (4)

by ctoney, 14 years ago

Attachment: added

comment:1 by Kyle Shannon, 14 years ago

Cc: Kyle Shannon added

comment:2 by Even Rouault, 14 years ago

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

comment:3 by Even Rouault, 5 years ago

Milestone: closed_because_of_github_migration
Resolution: wontfix
Status: newclosed

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 if it is not already reported there.

Note: See TracTickets for help on using tickets.