<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
On 11/09/2011 04:50 PM, Clifford J Mugnier wrote:
<blockquote
cite="mid:B597ADFC46B67741B18EA6AA13613FF7026A6D00@email002.lsu.edu"
type="cite">
<div dir="ltr" id="idOWAReplyText38541">
<div dir="ltr"><font color="#000000" face="Times New Roman"
size="2">Vincenzo,</font></div>
</div>
</blockquote>
Dear Clifford,<br>
<blockquote
cite="mid:B597ADFC46B67741B18EA6AA13613FF7026A6D00@email002.lsu.edu"
type="cite">
<div dir="ltr" id="idOWAReplyText38541">
<div dir="ltr"> </div>
<div dir="ltr"><font size="2">Since you did not mention where
you are, I had to look up the coordinates to give you a
clue. Since your coordinates are near Rome, the local
coordinate system is a Gauss-Boaga Transverse Mercator
referenced to the Monte Mario Datum of 1940 on the
International ellipsoid. For further details, see my column
on the <strong>Italian Republic</strong> in the August 2005
issue of <em>Photogrammetric Engineering & Remote
Sensing</em>. The pdf of that column may be downloaded
from the society's website at <a moz-do-not-send="true"
href="http://www.ASPRS.org">www.ASPRS.org</a>. </font></div>
<div dir="ltr"> </div>
</div>
<div dir="ltr" id="idSignature34563">
<div><font color="#000000" face="Times New Roman" size="2"><span
style="FONT-SIZE: 10pt">
<div style="MARGIN: 0in 0in 0pt" class="MsoNormal"><span
style="FONT-SIZE: 10pt; mso-ansi-language: DE"
lang="DE">Clifford J. Mugnier, C.P., C.M.S.</span></div>
</span></font></div>
</div>
</blockquote>
based on your advice i surfed the web and found:<br>
<br>
<meta http-equiv="content-type" content="text/html;
charset=ISO-8859-1">
<a
href="http://www.rigacci.org/wiki/doku.php/tecnica/gps_cartografia_gis/gauss_boaga_wgs84">http://www.rigacci.org/wiki/doku.php/tecnica/gps_cartografia_gis/gauss_boaga_wgs84</a><br>
<br>
So i tried to use pyproj to do the computation:<br>
<br>
In [50]: p = Proj(proj="tmerc", lat_0="0", lon_0="9", k="0.9996",
x_0="1500000", y_0="0", ellps="intl", units="m",
towgs84="-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68")<br>
<br>
In [51]: x,y = p(45.463681, 9.1881714)<br>
<br>
In [52]: print x,y<br>
5790656.63279 1259107.51994<br>
<br>
But it's still wrong, i should obtain something close to X=2159820
Y=1457635<br>
<br>
Any other suggestion?<br>
<br>
Thanks for all.<br>
<pre class="moz-signature" cols="72">--
Vincenzo Ampolo
<a class="moz-txt-link-freetext" href="http://vincenzo-ampolo.net">http://vincenzo-ampolo.net</a>
<a class="moz-txt-link-freetext" href="http://goshawknest.wordpress.com">http://goshawknest.wordpress.com</a></pre>
</body>
</html>