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

[Vincenzo Ampolo]

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:

http://www.rigacci.org/wiki/doku.php/tecnica/gps_cartografia_gis/gauss_boaga_wgs84

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", 
towgs84="-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68")

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
http://vincenzo-ampolo.net
http://goshawknest.wordpress.com

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20111109/8cdb521c/attachment.html>