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

[Chris Toney]

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?
 
Thanks.
 
Chris Toney
 
 
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]],
    PROJECTION["Albers_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]],
    PROJECTION["Hotine_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"]],
    PROJECTION["Hotine_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