<div>Hi </div>
<div>&nbsp;</div>
<div>I am doing this ....</div>
<div>&nbsp;</div>
<div>I have lat/lon of Top Left Corner of the Map and&nbsp; Bottom Right Corner of The Map</div>
<div>iam converting&nbsp;them to projected X,Y using this projection defintion </div>
<div>&nbsp;</div>
<div>World Mercator &lt;54004&gt; +proj=merc +latts=0 +lon0=0 +k=1.000000 +x0=0<br>+y0=0 +ellps=WGS84 +datum=WGS84 +units=m no_defs &lt;&gt;</div>
<div>&nbsp;</div>
<div>( I have found on some site that Google Maps use this projection !!)<br>&nbsp;</div>
<div>&nbsp;</div>
<div>For example: </div>
<div>&nbsp;</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;tl_lat = 42.8483;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tl_lon = -74.2935;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; br_lat = 42.3813;&nbsp;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; br_lon = -73.6509;<br>&nbsp;</div>
<div>&nbsp;</div>
<div>using pj_transform of proj4.. iam converting them to X,Y.</div>
<div>&nbsp;</div>
<div>
<p>top left xy </p>
<p>-8270314.5893 , 5259843.0071</p>
<p>btm right xy </p>
<p>-8198780.6845 , 5189458.4882<br></p>
<p>I calculate the Centre cX,cY using the top left X,Y and Btm right X,Y...</p>
<p>The Size of the Image I have generated from the Map is (1031,801)...</p>
<p>width = absolute value of (btm right X - top left X)</p>
<p>height&nbsp;= absolute value of (btm right&nbsp;Y - top left Y)</p>
<p>Now I am trying to&nbsp;get the pixel x,y from&nbsp;the projected X,Y&nbsp;using this formula</p>
<p>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;px = 515.5 +&nbsp; (1031/width)*(X-cX);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; py = 400.5&nbsp; -&nbsp; (801/height)*(Y-cY);<br></p>
<p>This gives correct pixels if i give top left or bottom right or the centre X,Y...</p>
<p>But for some point inside the Image ... </p>
<p>the calcualted pixel X&nbsp;is some 2/3 pixels away from actual pixel X&nbsp;</p>
<p>the calcualted pixel&nbsp;Y&nbsp;is some 8/10 pixels away from actual pixel Y</p>
<p>where am I going wrong ... is the projection iam using correct for such calculations ..???</p>
<p>&nbsp;</p>
<p>Thanks in advance.</p>
<p>sandeep</p></div><br><br>
<div><span class="gmail_quote">On 10/6/06, <b class="gmail_sendername">Gerald I. Evenden</b> &lt;<a href="mailto:gerald.evenden@verizon.mail">gerald.evenden@verizon.mail</a>&gt; wrote:</span>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">On Friday 06 October 2006 4:38 am, Hugo Tyson wrote:<br>&gt; <a href="mailto:proj-bounces@lists.maptools.org">
proj-bounces@lists.maptools.org</a> &lt;&gt; wrote:<br>&gt; &gt; I want to transform ... this lat/lon 39.739/-104.985...<br>&gt; &gt; to world mercator project ..with the defintition<br>&gt; &gt;<br>&gt; &gt; World Mercator &lt;54004&gt; +proj=merc +latts=0 +lon0=0
<br>&gt; &gt; +k=1.000000 +x0=0 +y0=0 +ellps=WGS84 +datum=WGS84 +units=m no_defs &lt;&gt;<br>&gt; &gt;<br>&gt; &gt; but the method gives an &quot;error&quot;....<br>&gt; &gt;<br>&gt; &gt; this method works fine if i give lat/lon say 
42.4182,-74.2643<br>&gt;<br>&gt; I've not heard of &quot;world mercator&quot;, but looking at the parameters<br>&gt; it's projecting into a plane tangent to the sphere at lat 0, lon 0.<br>&gt; So it only works for longitudes from -90 to +90 (though these map
<br>&gt; to infinities in the projected plane).<br>&gt;<br>&gt; Your example of -104 longitude is around the back of the earth,<br>&gt; on the &quot;far side&quot; from the projection plane.&nbsp;&nbsp; Therefore a<br>&gt; line from the origin through the point never intersects the
<br>&gt; projection plane, so no result can be calculated.&nbsp;&nbsp;Hence your<br>&gt; error.&nbsp;&nbsp;Lon -74 is in range, so it's OK.<br><br>Excuse me for jumping into the middle of this discussion, but Mercator and<br>normal proj functioning span the range of +-180 degrees longitude.&nbsp;&nbsp;The plane
<br>of the cylindrical projection projection is &quot;unwrapped&quot; from a cut on the<br>projection cylinder opposite the longitude of the central meridian (CM) or<br>180 degrees from the CM.&nbsp;&nbsp;There is no such thing as the '&quot;far side&quot;' of the
<br>projection.<br><br>The following is an example:<br><br>gie@charon:~$ lproj +proj=merc +ellps=WGS84<br>-104.985 39.739<br>-11686876.74&nbsp;&nbsp;&nbsp;&nbsp;4800765.45<br>-74.2643 42.4182<br>-8267064.06&nbsp;&nbsp;&nbsp;&nbsp; 5195000.64<br>gie@charon:~$<br>
<br>All the other parameters in the original problem are superfluous as they are<br>default values of proj/lproj.&nbsp;&nbsp;The datum option is not part of the proj/lproj<br>system.<br>_______________________________________________
<br>Proj mailing list<br><a href="mailto:Proj@lists.maptools.org">Proj@lists.maptools.org</a><br><a href="http://lists.maptools.org/mailman/listinfo/proj">http://lists.maptools.org/mailman/listinfo/proj</a><br></blockquote>
</div><br><br clear="all"><br>-- <br>Sandeep Kumar Jakkaraju <br>WeBlog:<br><a href="http://jakkarajus.blogspot.com">http://jakkarajus.blogspot.com</a>