Hi,
I have problem with converting geocentric WGS84 coordinates in GRASS.
Source location:
-PROJ_INFO-------------------------------------------------
name : EPSG Projection 4979 - WGS 84
proj : geocent
ellps : wgs84
-PROJ_UNITS------------------------------------------------
unit : metre
units : metres
meters : 1
Target location:
-PROJ_INFO-------------------------------------------------
name : Krovak
proj : krovak
datum : hermannskogel
ellps : bessel
lat_0 : 49.5
lon_0 : 24.83333333333333
alpha : 30.28813972222222
k : 0.9999
x_0 : 0
y_0 : 0
no_defs : defined
towgs84 : 570.8,85.7,462.8,4.998,1.587,5.261,3.56
-PROJ_UNITS------------------------------------------------
unit : metre
units : metres
meters : 1
Test point:
v.out.ascii fixed_points where="cat = 9291" fs=space
3979316.43738 1050312.25981 4857066.89549 9291
after reprojecting:
v.out.ascii fixed_points where="cat = 9291" fs=space
-1420272.16059322 -7048695.05085213 4857066.89549 9291
While manual reprojecting works:
echo "3979316.43738 1050312.25981 4857066.89549" | cs2cs +proj=geocent
+a=6378137 +rf=298.257223563 +no_defs +to +proj=krovak +lat_0=49.5
+lon_0=24.83333333333333 +k=0.9999 +x_0=0 +y_0=0 +no_defs
+a=6377397.155 +rf=299.1528128
+towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56
-719603.84 -1065734.53 592.60
Thanks for any pointer in advance, Martin
--
Martin Landa <landa.martin gmail.com> * http://gama.fsv.cvut.cz/~landa *