[gdal-dev] gdaltransform results differs from qgis

[Antonio Parrotta]

Hi,

I solved, it was because epsg:3003 was missing in epsg file and proj env
variable were not set. I have add

+proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=1500000 +y_0=0
+ellps=intl +units=m +no_defs

in epsg file
now I have same result as you.

Thanks a lot for the helping me!

cheers
Antonio


2008/9/15 Antonio Parrotta <antonioparrotta at gmail.com>

> I'm trying proj4 too, but after downloaded binary and
>
> set PROJ_LIB=<nad_intallation_dir>
> set set PATH=%PATH%;<bin_intallation_path>
>
> >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
>
> I obtain :
> Using to definition: init=epsg:3003
> towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68
> Rel. 4.4.6, 3 March 2003
> <cs2cs>:
> projection initialization failure
> cause: no options found in 'init' file
> program abnormally terminated
>
> maybe I'm missing some file,  is it?
>
> thanks
> Antonio
>
> 2008/9/15 Antonio Parrotta <antonioparrotta at gmail.com>
>
> I downloaded the source from http://download.osgeo.org/gdal/gdal152.zip,
>> built on windows XP with VC2008 with and PROJ4 dll from
>> http://download.osgeo.org/proj/proj446_win32_bin.zip
>>
>> any clue
>>
>> Antonio
>>
>> 2008/9/15 Markus Neteler <neteler at osgeo.org>
>>
>> On Mon, Sep 15, 2008 at 11:02 AM, Antonio Parrotta
>>> <antonioparrotta at gmail.com> wrote:
>>> > Hi All,
>>> >
>>> > I'm facing in  strange result using gdaltransform. What I'm doing is to
>>> try
>>> > to trasform coordinates from WGS84 to Gauss-Boaga:
>>> >
>>> > gdaltransform.exe -s_srs "WGS84" -t_srs "+proj=tmerc +ellps=intl
>>> +lat_0=0
>>> > +lon_0=9 +k=0.9996 +x_0=1500000 +y_0=0 +units=m
>>> > +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68"
>>> > entering lon and lat
>>> >>7,69746 45,106977<enter>
>>> > it gives
>>> > 1342393.88434212 4984915.97865081 -47.0342488819733
>>> > doing the same transformation with qgis (same parameters), I have
>>> > 1397559,446 4995677,02
>>> > that is quite near to the target expected also checking out in map.
>>> >
>>> > Do we know the reason of this delta? Do I have to set any other
>>> correction
>>> > parameter?
>>>
>>>
>>> Using cs2cs from PROJ4 Rel. 4.6.1, 21 July 2008, I get:
>>>
>>> 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
>>> # ---- From Coordinate System ----
>>> #Lat/long (Geodetic alias)
>>> #
>>> # +proj=latlong +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
>>> # ---- To Coordinate System ----
>>> #Transverse Mercator
>>> #       Cyl, Sph&Ell
>>> # +init=epsg:3003 +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68
>>> # +no_defs +proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=1500000 +y_0=0
>>> # +ellps=intl +units=m
>>> #--- following specified but NOT used
>>> # +no_defs
>>> 7.69746 45.106977
>>> 1397556.77      4995679.28 -46.70
>>>
>>> Using gdaltransform on Linux (almost current SVN), I get:
>>> gdaltransform -s_srs "WGS84" -t_srs "+proj=tmerc +ellps=intl +lat_0=0 \
>>>   +lon_0=9 +k=0.9996 +x_0=1500000 +y_0=0 +units=m \
>>>  +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68"
>>> 7.69746 45.106977
>>> 1397556.76994996 4995679.27916302 -46.6971342787147
>>>
>>> which is the same as the cs2cs output.
>>>
>>> How old is your gdaltransform.exe and where does it originate from?
>>>
>>> Markus
>>>
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20080915/57c41fe7/attachment.html