<div dir="ltr">Hi,<br><br>I solved, it was because epsg:3003 was missing in epsg file and proj env variable were not set. I have add<br><pre>+proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=1500000 +y_0=0 +ellps=intl +units=m +no_defs</pre>
in epsg file<br>now I have same result as you.<br><br>Thanks a lot for the helping me!<br><br>cheers <br>Antonio<br><br><br><div class="gmail_quote">2008/9/15 Antonio Parrotta <span dir="ltr"><<a href="mailto:antonioparrotta@gmail.com">antonioparrotta@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div dir="ltr">I'm trying proj4 too, but after downloaded binary and <br><br>set PROJ_LIB=<nad_intallation_dir><br>
set set PATH=%PATH%;<bin_intallation_path><div class="Ih2E3d"><br><br>>cs2cs -v +proj=latlong +datum=WGS84 +to +init=epsg:3003 +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68<br>
<br></div>I obtain :<br>Using to definition: init=epsg:3003 towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68<br>Rel. 4.4.6, 3 March 2003<br><cs2cs>:<br>projection initialization failure<br>cause: no options found in 'init' file<br>
program abnormally terminated<br><br>maybe I'm missing some file, is it?<br><br>thanks<br>Antonio<br><br><div class="gmail_quote">2008/9/15 Antonio Parrotta <span dir="ltr"><<a href="mailto:antonioparrotta@gmail.com" target="_blank">antonioparrotta@gmail.com</a>></span><div>
<div></div><div class="Wj3C7c"><br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div dir="ltr">I downloaded the source from <a href="http://download.osgeo.org/gdal/gdal152.zip" target="_blank">http://download.osgeo.org/gdal/gdal152.zip</a>, built on windows XP with VC2008 with and PROJ4 dll from <a href="http://download.osgeo.org/proj/proj446_win32_bin.zip" target="_blank">http://download.osgeo.org/proj/proj446_win32_bin.zip</a><br>
<br>any clue<br><br>Antonio<br><br><div class="gmail_quote">2008/9/15 Markus Neteler <span dir="ltr"><<a href="mailto:neteler@osgeo.org" target="_blank">neteler@osgeo.org</a>></span><div><div></div><div>
<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div><div></div><div>On Mon, Sep 15, 2008 at 11:02 AM, Antonio Parrotta<br>
<<a href="mailto:antonioparrotta@gmail.com" target="_blank">antonioparrotta@gmail.com</a>> wrote:<br>
> Hi All,<br>
><br>
> I'm facing in strange result using gdaltransform. What I'm doing is to try<br>
> to trasform coordinates from WGS84 to Gauss-Boaga:<br>
><br>
> gdaltransform.exe -s_srs "WGS84" -t_srs "+proj=tmerc +ellps=intl +lat_0=0<br>
> +lon_0=9 +k=0.9996 +x_0=1500000 +y_0=0 +units=m<br>
> +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68"<br>
> entering lon and lat<br>
>>7,69746 45,106977<enter><br>
> it gives<br>
> 1342393.88434212 4984915.97865081 -47.0342488819733<br>
> doing the same transformation with qgis (same parameters), I have<br>
> 1397559,446 4995677,02<br>
> that is quite near to the target expected also checking out in map.<br>
><br>
> Do we know the reason of this delta? Do I have to set any other correction<br>
> parameter?<br>
<br>
<br>
</div></div>Using cs2cs from PROJ4 Rel. 4.6.1, 21 July 2008, I get:<br>
<br>
cs2cs -v +proj=latlong +datum=WGS84 +to +init=epsg:3003<br>
<div>+towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68<br>
</div># ---- From Coordinate System ----<br>
#Lat/long (Geodetic alias)<br>
#<br>
# +proj=latlong +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0<br>
# ---- To Coordinate System ----<br>
#Transverse Mercator<br>
# Cyl, Sph&Ell<br>
# +init=epsg:3003 +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68<br>
# +no_defs +proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=1500000 +y_0=0<br>
# +ellps=intl +units=m<br>
#--- following specified but NOT used<br>
# +no_defs<br>
7.69746 45.106977<br>
1397556.77 4995679.28 -46.70<br>
<br>
Using gdaltransform on Linux (almost current SVN), I get:<br>
gdaltransform -s_srs "WGS84" -t_srs "+proj=tmerc +ellps=intl +lat_0=0 \<br>
<div> +lon_0=9 +k=0.9996 +x_0=1500000 +y_0=0 +units=m \<br>
+towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68"<br>
</div>7.69746 45.106977<br>
1397556.76994996 4995679.27916302 -46.6971342787147<br>
<br>
which is the same as the cs2cs output.<br>
<br>
How old is your gdaltransform.exe and where does it originate from?<br>
<font color="#888888"><br>
Markus<br>
</font></blockquote></div></div></div><br></div>
</blockquote></div></div></div><br></div>
</blockquote></div><br></div>