<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv=Content-Type content="text/html; charset=us-ascii">
<meta name=Generator content="Microsoft Word 12 (filtered medium)">
<style>
<!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;}
@page Section1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.Section1
{page:Section1;}
-->
</style>
<!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang=EN-US link=blue vlink=purple>
<div class=Section1>
<p class=MsoNormal>Hi,<o:p></o:p></p>
<p class=MsoNormal> 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?<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal> 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.<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>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.<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>Bellow is the command line I used to warp the image and info
dumps of the before and after images.<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>$ gdalwarp -t_srs ESRI::wm.prj Athens_RGB_CEC1FT_200303.tif
Athens_RGB_CEC1FT_200303-wm.tif<o:p></o:p></p>
<p class=MsoNormal>Creating output file that is 7455P x 6882L.<o:p></o:p></p>
<p class=MsoNormal>Processing input file Athens_RGB_CEC1FT_200303.tif.<o:p></o:p></p>
<p class=MsoNormal>0...10...20...30...40...50...60...70...80...90...100 - done.<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>Files: Athens_RGB_CEC1FT_200303.tif<o:p></o:p></p>
<p class=MsoNormal>Size is 6683, 6024<o:p></o:p></p>
<p class=MsoNormal>Coordinate System is:<o:p></o:p></p>
<p class=MsoNormal>PROJCS["USA_Contiguous_Equidistant_Conic",<o:p></o:p></p>
<p class=MsoNormal> GEOGCS["NAD83",<o:p></o:p></p>
<p class=MsoNormal>
DATUM["North_American_Datum_1983",<o:p></o:p></p>
<p class=MsoNormal>
SPHEROID["GRS 1980",6378137,298.2572221010002,<o:p></o:p></p>
<p class=MsoNormal>
AUTHORITY["EPSG","7019"]],<o:p></o:p></p>
<p class=MsoNormal>
AUTHORITY["EPSG","6269"]],<o:p></o:p></p>
<p class=MsoNormal>
PRIMEM["Greenwich",0],<o:p></o:p></p>
<p class=MsoNormal>
UNIT["degree",0.0174532925199433],<o:p></o:p></p>
<p class=MsoNormal>
AUTHORITY["EPSG","4269"]],<o:p></o:p></p>
<p class=MsoNormal> PROJECTION["Equidistant_Conic"],<o:p></o:p></p>
<p class=MsoNormal>
PARAMETER["standard_parallel_1",33],<o:p></o:p></p>
<p class=MsoNormal>
PARAMETER["standard_parallel_2",45],<o:p></o:p></p>
<p class=MsoNormal>
PARAMETER["latitude_of_center",39],<o:p></o:p></p>
<p class=MsoNormal>
PARAMETER["longitude_of_center",-96],<o:p></o:p></p>
<p class=MsoNormal> PARAMETER["false_easting",0],<o:p></o:p></p>
<p class=MsoNormal> PARAMETER["false_northing",0],<o:p></o:p></p>
<p class=MsoNormal> UNIT["metre",1,<o:p></o:p></p>
<p class=MsoNormal>
AUTHORITY["EPSG","9001"]]]<o:p></o:p></p>
<p class=MsoNormal>Origin = (1157874.742488616100000,-478125.149393160190000)<o:p></o:p></p>
<p class=MsoNormal>Pixel Size = (0.169893858067214,-0.169893858067214)<o:p></o:p></p>
<p class=MsoNormal>Corner Coordinates:<o:p></o:p></p>
<p class=MsoNormal>Upper Left ( 1157874.742, -478125.149) (
83d24'42.40"W, 33d58'12.44"N)<o:p></o:p></p>
<p class=MsoNormal>Lower Left ( 1157874.742, -479148.590) (
83d24'47.89"W, 33d57'39.54"N)<o:p></o:p></p>
<p class=MsoNormal>Upper Right ( 1159010.143, -478125.149) (
83d23'58.53"W, 33d58'7.37"N)<o:p></o:p></p>
<p class=MsoNormal>Lower Right ( 1159010.143, -479148.590) ( 83d24'4.02"W,
33d57'34.47"N)<o:p></o:p></p>
<p class=MsoNormal>Center ( 1158442.443,
-478636.870) ( 83d24'23.21"W, 33d57'53.46"N)<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>Files: Athens_RGB_CEC1FT_200303-wm.tif<o:p></o:p></p>
<p class=MsoNormal>Size is 7455, 6882<o:p></o:p></p>
<p class=MsoNormal>Coordinate System is:<o:p></o:p></p>
<p class=MsoNormal>PROJCS["WGS_1984_Web_Mercator",<o:p></o:p></p>
<p class=MsoNormal> GEOGCS["GCS_WGS_1984_Major_Auxiliary_Sphere",<o:p></o:p></p>
<p class=MsoNormal>
DATUM["WGS_1984_Major_Auxiliary_Sphere",<o:p></o:p></p>
<p class=MsoNormal>
SPHEROID["WGS_1984_Major_Auxiliary_Sphere",6378137,0]],<o:p></o:p></p>
<p class=MsoNormal>
PRIMEM["Greenwich",0],<o:p></o:p></p>
<p class=MsoNormal>
UNIT["degree",0.0174532925199433]],<o:p></o:p></p>
<p class=MsoNormal> PROJECTION["Mercator_1SP"],<o:p></o:p></p>
<p class=MsoNormal> PARAMETER["central_meridian",0],<o:p></o:p></p>
<p class=MsoNormal> PARAMETER["scale_factor",1],<o:p></o:p></p>
<p class=MsoNormal> PARAMETER["false_easting",0],<o:p></o:p></p>
<p class=MsoNormal> PARAMETER["false_northing",0],<o:p></o:p></p>
<p class=MsoNormal> UNIT["metre",1,<o:p></o:p></p>
<p class=MsoNormal>
AUTHORITY["EPSG","9001"]]]<o:p></o:p></p>
<p class=MsoNormal>Origin = (-9285526.535261498800000,4000908.274406374400000)<o:p></o:p></p>
<p class=MsoNormal>Pixel Size = (0.204760916615493,-0.204760916615493)<o:p></o:p></p>
<p class=MsoNormal>Corner Coordinates:<o:p></o:p></p>
<p class=MsoNormal>Upper Left (-9285526.535, 4000908.274) (
83d24'47.89"W, 33d47'31.24"N)<o:p></o:p></p>
<p class=MsoNormal>Lower Left (-9285526.535, 3999499.110) ( 83d24'47.89"W,
33d46'53.36"N)<o:p></o:p></p>
<p class=MsoNormal>Upper Right (-9284000.043, 4000908.274) (
83d23'58.53"W, 33d47'31.24"N)<o:p></o:p></p>
<p class=MsoNormal>Lower Right (-9284000.043, 3999499.110) (
83d23'58.53"W, 33d46'53.36"N)<o:p></o:p></p>
<p class=MsoNormal>Center (-9284763.289,
4000203.692) ( 83d24'23.21"W, 33d47'12.30"N)<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
</div>
</body>
</html>