<div dir="ltr">I downloaded the source from <a href="http://download.osgeo.org/gdal/gdal152.zip">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">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">neteler@osgeo.org</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><div></div><div class="Wj3C7c">On Mon, Sep 15, 2008 at 11:02 AM, Antonio Parrotta<br>
<<a href="mailto:antonioparrotta@gmail.com">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 class="Ih2E3d">+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 class="Ih2E3d"> +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><br></div>