[Proj] load arcSDE raw data from database and convert to lat/long

Vincenzo Ampolo vincenzo.ampolo at gmail.com
Wed Nov 9 11:15:52 EST 2011

On 11/09/2011 04:50 PM, Clifford J Mugnier wrote:
> Vincenzo,
Dear Clifford,
> 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 *Italian 
> Republic* in the August 2005 issue of /Photogrammetric Engineering & 
> Remote Sensing/.  The pdf of that column may be downloaded from the 
> society's website at www.ASPRS.org <http://www.ASPRS.org>.
> Clifford J. Mugnier, C.P., C.M.S.
based on your advice i surfed the web and found:


So i tried to use pyproj to do the computation:

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", 

In [51]: x,y = p(45.463681, 9.1881714)

In [52]: print x,y
5790656.63279 1259107.51994

But it's still wrong, i should obtain something close to X=2159820 Y=1457635

Any other suggestion?

Thanks for all.

Vincenzo Ampolo

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.maptools.org/pipermail/proj/attachments/20111109/8cdb521c/attachment.htm 

More information about the Proj mailing list