Hi,
I'm having problems warping an image to WGS 1984 Web Mercator. The lat/long
center point moved south by 10 arc-minutes. Does anyone have some suggestion
on how to fix this?
I'm pretty sure the USA_Contiguous_Equidistant_Conic transform worked fine
because I looked at the lag/long in Google Earth and saw the same buildings.
I've also done the transform in two steps Source SRS to WGS84 and WGS84 to WGS
1984 Web Mercator but get the same answer.
Bellow is the command line I used to warp the image and info dumps of the
before and after images.
$ gdalwarp -t_srs ESRI::wm.prj Athens_RGB_CEC1FT_200303.tif
Athens_RGB_CEC1FT_200303-wm.tif
Creating output file that is 7455P x 6882L.
Processing input file Athens_RGB_CEC1FT_200303.tif.
0...10...20...30...40...50...60...70...80...90...100 - done.
Files: Athens_RGB_CEC1FT_200303.tif
Size is 6683, 6024
Coordinate System is:
PROJCS["USA_Contiguous_Equidistant_Conic",
GEOGCS["NAD83",
DATUM["North_American_Datum_1983",
SPHEROID["GRS 1980",6378137,298.2572221010002,
AUTHORITY["EPSG","7019"]],
AUTHORITY["EPSG","6269"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4269"]],
PROJECTION["Equidistant_Conic"],
PARAMETER["standard_parallel_1",33],
PARAMETER["standard_parallel_2",45],
PARAMETER["latitude_of_center",39],
PARAMETER["longitude_of_center",-96],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (1157874.742488616100000,-478125.149393160190000)
Pixel Size = (0.169893858067214,-0.169893858067214)
Corner Coordinates:
Upper Left ( 1157874.742, -478125.149) ( 83d24'42.40"W, 33d58'12.44"N)
Lower Left ( 1157874.742, -479148.590) ( 83d24'47.89"W, 33d57'39.54"N)
Upper Right ( 1159010.143, -478125.149) ( 83d23'58.53"W, 33d58'7.37"N)
Lower Right ( 1159010.143, -479148.590) ( 83d24'4.02"W, 33d57'34.47"N)
Center ( 1158442.443, -478636.870) ( 83d24'23.21"W, 33d57'53.46"N)
Files: Athens_RGB_CEC1FT_200303-wm.tif
Size is 7455, 6882
Coordinate System is:
PROJCS["WGS_1984_Web_Mercator",
GEOGCS["GCS_WGS_1984_Major_Auxiliary_Sphere",
DATUM["WGS_1984_Major_Auxiliary_Sphere",
SPHEROID["WGS_1984_Major_Auxiliary_Sphere",6378137,0]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Mercator_1SP"],
PARAMETER["central_meridian",0],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (-9285526.535261498800000,4000908.274406374400000)
Pixel Size = (0.204760916615493,-0.204760916615493)
Corner Coordinates:
Upper Left (-9285526.535, 4000908.274) ( 83d24'47.89"W, 33d47'31.24"N)
Lower Left (-9285526.535, 3999499.110) ( 83d24'47.89"W, 33d46'53.36"N)
Upper Right (-9284000.043, 4000908.274) ( 83d23'58.53"W, 33d47'31.24"N)
Lower Right (-9284000.043, 3999499.110) ( 83d23'58.53"W, 33d46'53.36"N)
Center (-9284763.289, 4000203.692) ( 83d24'23.21"W, 33d47'12.30"N)
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
