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 * _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
