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
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
