Opened 15 years ago
Closed 15 years ago
#3037 closed defect (duplicate)
incorrect result from gdalwarp - NAD83 Albers to Michigan GeoRef
Reported by: | ctoney | Owned by: | warmerdam |
---|---|---|---|
Priority: | normal | Milestone: | |
Component: | default | Version: | 1.6.1 |
Severity: | normal | Keywords: | |
Cc: |
Description
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