User Tools

Site Tools


etrs89

Funkce v jazyku C pro výpočet transformace mezi WGS-84 a ETRS-89

void ITRF05tc_ETRS89tc(long double *x, long double *y, long double *z, long double *xE, long double *yE, long double *zE) 
{ 
 long double T1, T2, T3, R1, R2, R3;

 T1 = 5.6*0.01; 
 T2 = 4.8*0.01; 
 T3 = -3.7*0.01; 
 R1 = (0.054*0.001)/3600*M_PI/180; 
 R2 = (0.518*0.001)/3600*M_PI/180; 
 R3 = (-0.781*0.001)/3600*M_PI/180;

 *xE = *x + T1 + ((0 - R3 * *y + R2 * *z)*(tc - 1989)); 
 *yE = *y + T2 + ((R3 * *x + 0 - R1 * *z)*(tc - 1989)); 
 *zE = *z + T3 + ((-R2 * *x + R1 * *y + 0)*(tc - 1989));

 fprintf (v_ETRS, "%Lf\t", *xE); 
 fprintf (v_ETRS, "%Lf\t", *yE); 
 fprintf (v_ETRS, "%Lf\t", *zE);
etrs89.txt · Last modified: 2018/01/31 14:37 by Jan Forman