Bookmark this page: Add coord transformation with Proj4 to Yahoo MyWeb Add coord transformation with Proj4 to Google Bookmarks Add coord transformation with Proj4 to Windows Live Add coord transformation with Proj4 to Del.icio.us Digg coord transformation with Proj4! Add coord transformation with Proj4 to Netscape
  •  
  • Subject
  • Author
  • Date
If you were  Registered and logged in, you could reply and use other advanced thread options
Posted by Zoran on February 24, 2009, 3:28 pm


Hello!
Does anyone use proj4 with 7-parameter for coordinate transformation.
I try to convert wgs84 to Croatian HDKS but unsuccesfull.
I use this command line
cs2cs.exe -v -E +proj=latlong +ellps=GRS80
+towgs84=-528.6481,-119.5964,-467.1873,6.061424,2.749619,-9.525554,1.064918
+to +proj=tmerc +ellps=bessel +lat_0=0 +lon_0=15 +k=0.9999
+x_0=5500000 +y_0=0 inwgs > out.txt

Inwgs:
14 45 100

out.txt:
# ---- From Coordinate System ----
#Lat/long (Geodetic)
#        
# +proj=latlong +ellps=GRS80
#
+towgs84=-528.6481,-119.5964,-467.1873,6.061424,2.749619,-9.525554,1.064918
# ---- To Coordinate System ----
#Transverse Mercator
#        Cyl, Sph&Ell
# +proj=tmerc +ellps=bessel +lat_0=0 +lon_0=15 +k=0.9999 +x_0=5500000
+y_0=0
14 45 100        5421169.76        4984360.80 806.66


But this is not correct result.
I have correct result when i transform with alltrans using same 7
parameters
5421533.2467 4984445.3383 100.1539

Any idea?



Makni _BRISI_ iz mail adrese

Posted by miso@sushi.com on February 26, 2009, 11:34 pm


> Hello!
> Does anyone use proj4 with 7-parameter for coordinate transformation.
> I try to convert wgs84 to Croatian HDKS but unsuccesfull.
> I use this command line
> cs2cs.exe -v -E =A0+proj=3Dlatlong +ellps=3DGRS80
> +towgs84=3D-528.6481,-119.5964,-467.1873,6.061424,2.749619,-9.525554,1.06=
4918
> +to +proj=3Dtmerc +ellps=3Dbessel +lat_0=3D0 +lon_0=3D15 +k=3D0.9999
> +x_0=3D5500000 +y_0=3D0 inwgs > out.txt
> Inwgs:
> 14 45 100
> out.txt:
> # ---- From Coordinate System ----
> #Lat/long (Geodetic)
> # =A0 =A0 =A0
> # +proj=3Dlatlong +ellps=3DGRS80
> #
> +towgs84=3D-528.6481,-119.5964,-467.1873,6.061424,2.749619,-9.525554,1.06=
4918
> # ---- To Coordinate System ----
> #Transverse Mercator
> # =A0 =A0 =A0 Cyl, Sph&Ell
> # +proj=3Dtmerc +ellps=3Dbessel +lat_0=3D0 +lon_0=3D15 +k=3D0.9999 +x_0=
=3D5500000
> +y_0=3D0
> 14 45 100 =A0 =A0 =A0 5421169.76 =A0 =A0 =A04984360.80 806.66
> But this is not correct result.
> I have correct result when i transform with alltrans using same 7
> parameters
> 5421533.2467 4984445.3383 100.1539
> Any idea?
> Makni _BRISI_ iz mail adrese

Proj4 had a dedicate mailing list.

Posted by Zoran on March 2, 2009, 3:16 pm


Thank you...i looking for it.


On Thu, 26 Feb 2009 20:34:43 -0800 (PST), "miso@sushi.com"

>> Hello!
>> Does anyone use proj4 with 7-parameter for coordinate transformation.
>> I try to convert wgs84 to Croatian HDKS but unsuccesfull.
>> I use this command line
>> cs2cs.exe -v -E  +proj=latlong +ellps=GRS80
>> +towgs84=-528.6481,-119.5964,-467.1873,6.061424,2.749619,-9.525554,1.064918
>> +to +proj=tmerc +ellps=bessel +lat_0=0 +lon_0=15 +k=0.9999
>> +x_0=5500000 +y_0=0 inwgs > out.txt
>> Inwgs:
>> 14 45 100
>> out.txt:
>> # ---- From Coordinate System ----
>> #Lat/long (Geodetic)
>> #      
>> # +proj=latlong +ellps=GRS80
>> #
>> +towgs84=-528.6481,-119.5964,-467.1873,6.061424,2.749619,-9.525554,1.064918
>> # ---- To Coordinate System ----
>> #Transverse Mercator
>> #       Cyl, Sph&Ell
>> # +proj=tmerc +ellps=bessel +lat_0=0 +lon_0=15 +k=0.9999 +x_0=5500000
>> +y_0=0
>> 14 45 100       5421169.76      4984360.80 806.66
>> But this is not correct result.
>> I have correct result when i transform with alltrans using same 7
>> parameters
>> 5421533.2467 4984445.3383 100.1539
>> Any idea?
>> Makni _BRISI_ iz mail adrese
>Proj4 had a dedicate mailing list.



Makni _BRISI_ iz mail adrese

Posted by Nick on March 3, 2009, 4:04 am


Hi Zoran,

Not a solution to your problem I'm afraid, but there is an easier way
of testing Proj.4 parameters, using QGIS (Open Source). If you go
Settings -> Custom CRS... it brings up a utility where you can test
parameters by cutting and pasting them into a text box. It's much
easier to edit the parameters and test them using QGIS than to work
from a terminal.

Regards, N.