[GRASS-user] geocentric coordinates in GRASS

[Martin Landa]

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 *