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

[Noel Zinn (cc)]

Vincenzo,

Here’s a radical suggestion for your consideration.  Approach the source of these coordinates and get projection, grid, ellipsoid and datum shift parameters.  Forensic sleuthing is great entertainment, but without definitive documentation it’s just an informed guess.

Noel

Noel Zinn, Principal, Hydrometronics LLC
+1-832-539-1472 (office), +1-281-221-0051 (cell)
noel.zinn at hydrometronics.com (email)
http://www.hydrometronics.com (website)


From: Vincenzo Ampolo 
Sent: Wednesday, November 09, 2011 10:15 AM
To: PROJ.4 and general Projections Discussions 
Subject: Re: [Proj] load arcSDE raw data from database and convert tolat/long

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. 

  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

--------------------------------------------------------------------------------
_______________________________________________
Proj mailing list
Proj at lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20111109/916c2c44/attachment.html>