Next: , Previous: SQL schema SQLite and gama-local, Up: Top

## 4 Network adjustment with gama-local

Adjustment of local geodetic network is a classical case of adjustment of indirect observations. After estimation of approximate values of unknown parameters (coordinates of points) and linearization of functions describing relations between observations and parameters we solve linear system of equations

     (1)      Ax = b + v,


where A is coefficient matrix, b is vector of absolute terms (right hand side) and v is vector of residuals. This system is (generally) overdetermined and we seek the solution x satisfying the basic criterion of Least Squares

     (2)      v'Pv = min,


where P is weight matrix. This criterion unambiguously defines the shape of adjusted network.

In the case of free network the system (1) is singular (matrix A has linearly dependent columns) and we have to define second regularization criterion

     (3)      \sum x_i^2 = min,      for all selected i


stating that at the same time we demand that the sum of squares corrections of selected parameters is minimal (corrections of unknown parameters with indexes from the set of all selected unknowns. Geometrically this criterion is equivalent to adjustment of the network according to (2) with simultaneous transformation to the selected set of fiducial points. This transformation does not change the shape of adjusted network.

Often it is advantageous to work with a homogenized system, ie. with the system of project equations in which coefficient of each row and absolute term are multiplied by square root of the weight of corresponding observation.

     (4)      ~A x = ~b,


where ~A = P^{1/2} A, ~b = P^{1/2} A. Symbol P^{1/2} denotes diagonal matrix of square roots of observation weights (or Cholesky decomposition of covariance matrix in the case of correlated observations). To criterion (2) corresponds in the case of homogenized system criterion

     (5)      ~v'v = min.


Normal equations are clearly equivalent for both systems.