## 1 Introduction

Geodetic deformation analysis is about change of form and size of the earth’s surface or of objects on, below or above it, and also of the relative position and orientation of the objects. The objects to be analysed are represented by points that constitute a three-dimensional geodetic network. It is nowadays common practice to use total stations, GPS receivers and other devices for the analysis. If processing the measurements results in three-dimensional *x*, *y*, *z* coordinates, these can be presented in two-dimensional graphs, showing the displacements in time or in space or both. It is, however, generally difficult to come to statistically sound conclusions by analysing the graphs. Computational methods to test statistical hypotheses are desirable. For two epochs methods are available to perform an adjustment of the coordinates, taking into account their covariance matrix, and to perform hypothesis testing [4, p. 494 ff.]. The analysis is generally not invariant for the used S-basis [10].

In this paper an adjustment model is proposed that analyses a time series of 3D coordinates, taking account of the covariance matrices and analysing the deformations of all points and all epochs simultaneously, by computing statistics of deformation patterns and testing them. The model can be applied to any 3D geodetic network, observed quasi-continuously (i. e. with permanently installed sensors measuring frequently). Examples are the monitoring of the movement of a subset of points through all epochs, or the periodic oscillation of a subset of points.

In the next section the problem is defined. After describing existing approaches in section 3, section 4 describes the solution set-up. The adjustment model is treated in section 5. The adjustment itself, the deformation testing and the S-basis invariance are handled in section 6. Section 7 gives an experimental validation of the model.

## 2 Problem Definition

The problem addressed in this paper is the adjustment and testing for deformations of a time series of three-dimensional coordinates of a geodetic network, with a covariance matrix of the coordinates that is full and generally singular, because each epoch of the time series is adjusted as a free network.

The described problem will be handled by constructing a least squares adjustment model. As a practical application to show the usability of the model, the continuous monitoring by a total station of points, situated on built structures that are prone to deformations, is analysed.

## 3 Existing Solutions

Heunecke et al. give an overview of existing approaches for deformation analysis of two time epochs of deformation measurements [4, p. 521]. The general approach is to compute displacement vectors between coordinates of two epochs and their covariance matrix. Different approaches exists to analyse the displacement vectors, e. g. by using 95 %-confidence ellipsoids after a least squares adjustment [2, 5] or a L1-norm adjustment [1, 3], or using constraints on common points and analysing the quadratic form of the weighted estimated least squares residuals [4, p. 500 ff.].

Heunecke et al. give also methods to analyse time series [4, p. 548]. They do not take advantage of the covariance matrices of the coordinates and do not perform the analysis for all three dimensions (*x*, *y* and *z*) simultaneously. As a consequence the choice of datum definition and the solution’s invariance for it, are not addressed.

A comprehensive 3D multi-epoch model is treated by [2, p. 164 ff.]. It takes care of singular covariance matrices and incorporates deterministic deformation models. It assumes all epochs to be defined relative to the same S-basis, which has to be defined by points, measured in all epochs. Testing is treated for the sequential adjustment case. Quality description of the tests is not treated.

## 4 Solution Set-up

### 4.1 Form and Size, Position and Orientation

The subject of geodetic deformation analysis is the change in time of the form and size of objects, and also of the relative position and orientation of the objects. Form, size, relative position and orientation can be recorded by Euclidian *x*, *y*, *z* coordinates. It is assumed that the coordinates are normally distributed with a probability density function, which is fully described by a known covariance matrix, except for the first moments. This matrix may be singular, e. g. because it stems from a free network adjustment. If there are reference points, i. e. points that are considered not to be influenced by the deformation to be analysed, they are part of the geodetic network, and are analysed simultaneously with the object points.

The Euclidian coordinates describe the position and orientation of the network relative to the coordinate origin and axes as well. These, however, are *not* subject of the analysis. Their uncertaintity, as it is reflected in the covariance matrix, has to be eliminated from the analysis. This is realised in the adjustment model by a congruence or similarity transformation of the coordinates of each epoch to the coordinate system of the reference epoch. It is shown that after these transformations, testing of deformation hypotheses can be done independently from the S-bases chosen for the individual epochs. The first epoch is chosen in this paper as reference epoch, but any other epoch as reference epoch would give the same analysis results.

The choice between a congruence and a similarity transformation depends on the question, whether the scale (unit of length) is considered stable between epochs and essential for the analysis.

The set-up of the adjustment model, with transformations between the epochs incorporated into it, not only removes the influence of origin, axes and scale of the reference system on the analysis. It also makes it possible to test, without additional S-transformations, for deformations of all kinds of subsets of points, independent of their being reference or object points, or being part of the S-basis or not. It is possible to include in one hypothesis that is to be tested, both reference and object points, and both points within and outside the S-basis.

### 4.2 Nonstochastic Observations

The adjustment model is built as a model of observation equations with the coordinates as observations, arranged according to the epochs. The parameters are the expectations of the coordinates of all epochs. Each point has for each epoch different coordinates in the parameter vector. Also the transformation parameters of each epoch relative to the previous one appear in the parameter vector.

Constraints are stated concerning the coordinates of all epochs. In the case of stability analysis the constraints state that the expectations of coordinates of the same points in different epochs are equal. These constraints are added to the observation vector as nonstochastic observations, following the approach of [9]. If coordinates are assumed to be subject to some kind of deformation, for example a linear movement of one or more points, or a deformation pattern with a periodic character for a subset of points, the constraints add extra unknown parameters, for example the linear rate of movement, or the coefficients of the periodic pattern, to the parameter vector.

The advantage of using nonstochastic observations is that testing of deformation hypotheses is done in the same way as testing of one- or multi-dimensional hypotheses on biases in the other observations. Least squares estimates of the deformations are determined using standard formulas. Also minimal detectable biases can be computed with standard formulas, giving information on the deformation sizes that can be detected with the tests.

### 4.3 Full, Singular Covariance Matrices

Observations, for example direction and distance observations of total stations, and their stochastic model are used for a deformation analysis, which is performed in two phases. In the first phase the direction and distance observations are adjusted for each epoch separately. The results are coordinates and their covariance matrices for all epochs. The second step is the subject of this paper: the deformation analysis of the coordinates of many epochs. The covariance matrices of all epochs have to be used [7, p. 154]. These matrices are generally full matrices (no or few zeros) and singular, because each epoch is adjusted as a free network, not connected to control points. The adjustment model of section 5 can handle full, singular covariance matrices.

### 4.4 Solution Characteristics

An overview of the solution characteristics can now be given. The most relevant terms are listed below and the solution procedure is illustrated by a Nassi-Schneidermann diagram (figure 1).

**A geodetic network per epoch** is a set of points on, above, or under the earth’s surface, in this paper assumed to be represented by 3D Euclidian coordinates.

**Form and size** of a geodetic network (and their changes in time) are of interest, not the position and orientation relative to the reference system. Transformations are therefore included in the adjustment model.

**Stability** assumes the expectations of coordinates to be equal through all epochs, except for the above mentioned transformations.

**A deformation pattern** is the relation between the geodetic networks per epoch, formulated by giving the expectations of coordinates through the epochs using mathematical functions, described by nonstochastic observations, which depend on unknown deformation parameters, for example a linear movement rate, or the coefficients of a series expansion of a periodic function.

**Inside or outside** the adjustment model we put the description of the deformation. If it is *inside*, nonstochastic observations describe the deformation pattern, and extra deformation parameters are included in the parameter vector. If it is *outside*, the adjustment model assumes *stability* and no extra deformation parameters are in the parameter vector. The hypothesis of stability is tested against alternative hypotheses, describing deformation patterns, by appropiate test statistics, which make use of the nonstochastic observations to determine matrix **Z**_{∇} of equation (16).

**Singular, full covariance matrices** result from the free network adjustments of each epoch, and are used in the adjustment model. If only coordinates are available for each epoch, a substitute matrix, for example a unit matrix, is used as covariance matrix, which yields sub optimal adjustment and testing results.

## 5 Adjustment Model

### 5.1 Observations and Parameters

The adjustment model is built taking as:

- –observations:
- i.cartesian 3D point coordinates of a geodetic network and their covariance matrix, available for at least two epochs. For the first epoch they are assembled in vector
**a**_{1}(an underlined variable indicates a stochastic variable) with the covariance matrix**D**{**a**_{1}}. For the second and later epochs they are assembled in vectors**b**_{i}, with*i*the epoch number, which runs from 2 to*p*with*p*the number of epochs. Each**b**_{i}has a covariance matrix**D**{**b**_{i}}; - ii.nonstochastic observations
**z**_{f}, describing constraints on the transformation parameters; their covariance matrix is the zero matrix; - iii.nonstochastic observations
**z**_{d}, describing the deformation pattern; their covariance matrix is the zero matrix.

- i.cartesian 3D point coordinates of a geodetic network and their covariance matrix, available for at least two epochs. For the first epoch they are assembled in vector
- –unknown parameters:
- i.expectations of cartesian 3D network point coordinates, for each epoch assembled in vector
**c**_{i}of epoch*i*. Vector**c**takes all epochs together:$c={\left({c}_{1}\mathrm{,}\cdots \mathrm{,}{c}_{p}\right)}^{T}.$ - ii.vector of transformation parameters
**f**, subdivided in subvectors**f**_{i, i–1}for the transformation in each epoch interval between epoch*i*and*i*– 1, with*i*= 2, . . .,*p*. - iii.additonal parameters ∇ to describe the trend function of the deformation, see section 4.2.

- i.expectations of cartesian 3D network point coordinates, for each epoch assembled in vector

### 5.2 Nonlinear Adjustment Model

In the adjustment model the expectations of all point coordinates are expressed in the reference system of the first epoch, and are *parameters* in vector **c**. The *observed* coordinates in the first epoch are taken together in vector **a**_{1}. We have:

with *E*{.} the expectation operator, and **P**_{1} the matrix that selects the points of the first epoch from **c**. **P**_{1} has only ones and zeros. The observed coordinates **b**_{i} in a following epoch *i* (*i* = 2, . . ., *p*) are assumed to be in a separate reference system, indicated by a superindex (*i*):

These coordinates are transformed with a vector function *ϕ*_{i, i–1} to the reference system of epoch (*i *– 1):

then with *ϕ*_{i–1, i–2} to the reference system of epoch (*i *– 2), and so on, and we get the transformed coordinates **a**_{i}:

and:

**P**_{i} selects **c**_{i} from **c**. It follows that:

Hopping from epoch *i* through all intermediate epochs to the first one, is chosen, and not a direct transformation to the first epoch, because it is assumed that in general more common points are available for successive epochs.

Following the approach of [8], the transformation *ϕ*_{i,j} between epoch *i* and *j*, is of a general form, for example an affine transformation, which is changed to another type of transformation, for example a similarity transformation, by the use of constraints. These constraints are formulated as nonstochastic observations:

for which zeros are assumed as observed values.

The deformation pattern is described by a vector of nonstochastic observations **z**_{d} and a deformation function ζ_{d}, which gives a relation between the elements of **c** and the elements of a vector of deformation parameters ∇:

For **z**_{d} we assume zeros as the observed values.

From equations (2), (7), (8) and (9) follows the following system for *p* epochs:

with

A point may be present in an epoch, but missing in one or more other epochs. This is handled by matrix **P**_{i}. The S-basis definition of an epoch is arbitrary (see section 6.3) and may be realised by only a few points, by many points, or by all. There can even be no S-basis, i. e. the covariance matrix is regular, and the S-basis can be considered to lie outside the geodetic network. The fact that a point is missing, be it in the first or in any other epoch, does therefore not pose any problem for the deformation analysis with model (10).

### 5.3 Transformations

### 5.3.1 Affine Transformation

As general form of transformation *ϕ*_{i, j} the affine transformation is taken, written as:

The column vectors **u**, **v**, **w** contain resp. the *x*˗, *y*˗, *z*˗ coordinates of **b**_{i} before transformation. The vectors **x**, **y**, **z** contain the coordinates after transformation. **R** describes the rotation, shear and scale change of the affine transformation, **t** the translation.

### 5.3.2 Congruence Transformation

Equation (11) describes a congruence (or rigid body) transformation, if the nine coefficients of matrix **R** meet the following six constraints:

In the following sections a linearised adjustment model is derived. The linearised constraints are:

where **0** is the (1 × 3) zero vector and **a**_{i}. ∆ indicates the difference of the quantity concerned and its approximate value.

### 5.3.3 Similarity Transformation

For the similarity transformation the affine transformation is constrained with five constraints. Three constraints state that the three rows of **R** are perpendicular to each other. The two remaining constraints state that the lengths of the first and second row, and those of the second and third row are equal. The linearised constraints are:

### 5.3.4 Approximate Transformation

Before the adjustment, **b**_{i} is approximately transformed to **a**_{1}, using equation (5). Likewise **D**{**b**_{i}} is transformed to *a*_{ij} in the constraints of the congruence or similarity transformation we can take

In the following sections, if **b**_{i} is written,

### 5.4 Linearised Adjustment Model

A linearised adjustment model is built for the deformation analysis. Linearisation of all equations of system (10) is done with implicit differentiation relative to the observed vectors **a**_{1}, **b**_{i} (*i* = 2, . . ., *p*), **z**_{d} and **z**_{f}, and the unknown parameter vectors **f**, **c** and ∇.

The first and last two equations of system (10) are linearised as:

We define for later use:

The partial derivatives of the vectors ζ_{f} and ζ_{d} with respect to the vectors **f**, **c** and ∇ are matrices. The parentheses with zero (.)_{0} indicate that approximate values of the parameters have to be used to get the values in the matrices.

For the equations with **b**_{i} (*i* = 2, . . ., *p*) in system (10) the linearised equations are:

with the matrices **B**_{i} defined as follows:

and with ( *j* = 2, . . ., *i *– 1):

**F**_{i} is defined for *i* = 2, . . . *p* as follows:

with (*p*-*i*) matrices **0** of zeros, which have the same number of rows as **F**_{2,1 }, and the partioning of **F**_{i} in columns in accordance with the partitioning of ∆**f**:

For **F**_{i, i–1} (*i* = 2, . . ., *p*) we have:

Matrix **B**_{i, i–1} for an affine transformation is given by [8] as follows:

with *a*_{ij} and **I** the (*n* × *n*) unit matrix and *n* the amount of points in **b**_{i}.

As explained in section 5.3.4, we can take **B**_{i, i–1}, from which follows, see equations (18) and (22):

Matrix **F**_{i, i–1} for an affine transformation is given by [8] as follows:

where *β*_{i}, *ε*_{1}, *ε*_{2}, *ε*_{3} and **0** are all (*n* × 3) matrices, as follows:

*β*_{i} = (**u**_{0}, **v**_{0}, **w**_{0}); **u**_{0}, **v**_{0}, **w**_{0} are approximate values of **u**, **v**, **w**, (the *x*, *y*, z coordinates in **b**_{i}), which can be transformed to make the barycentre the origin.

*ε*_{2} and *ε*_{3} are analogous matrices as *ε*_{1} with ones in the second, resp. third column, **0** is the (*n* × 3) zero matrix.

We define **F**_{1} as the null matrix **0** and put it together with the **F**_{i}, *i* = 2 . . . *p* of equation (20) into matrix **F**. Analogously we take all **P**_{i} together in a matrix **P**:

We define vector ∆**b** as:

We can now formulate the linearised equivalent of system (10):

The covariance matrix of the observation vector on the left-hand side consists of the covariance matrices of **a**_{1} and **b**_{i}, *i *= 2, . . .,*p,* as described in section 5.1, approximately transformed as described in section 5.3.4, and zero matrices for the remainder if no correlation between the epochs is assumed (which is, however, not necessary to solve the model).

The model takes each epoch as a separate geodetic network: each point has a different point number for each epoch, for example point *A* is called *A*_{1} in epoch 1, *A*_{2} in epoch 2, etc. The hypothesis that no deformation has occured is formulated by stating that

The separate geodetic networks are linked together in this way. Equation (29) gives the nonstochastic observation equations (the zeros constitute together vector **z**_{d} and have a standard deviation of zero). The number of rows of matrix **Z**_{d} is three times the number of points. In each row there are zeros and one 1 and one −1 for respectively the coordinate of epoch 2 and epoch 1 (which are separate unknowns in the parameter vector). There are no parameters ∇ and no matrix **Z**_{∇}.

Let us now assume that a deformation is present for point *A*. Let it be a linear movement for which we write:

The *a*_{x}*, a*_{y}*, a*_{z} are unknown parameters, which enter the parameter vector ∇, and for which a least squares estimate is determined in the adjustment. *t*_{12} is the time interval between epoch 1 and 2. The matrix **Z**_{∇} is in this case a matrix with three columns and three elements *t*_{12} on the rows of the three nonstochastic observations mentioned, and with zeros on all other positions.

We can also leave *a*_{x}*, a*_{y}*, a*_{z} out of the adjustment. Then the last column of the coefficient matrix of equation (28) disappears. The null hypothesis states now that there is no deformation. We test for a linear movement by using **Z**_{∇} in the test statistic of equation (34).

Generally the transformation between epoch *i* and *i – *1 is a similarity or congruence, not an affine transformation. Matrix **F**_{i} is constructed according to equation (20) from matrices **F**_{i, i –1} as given in equation (22) for the affine transformation. Matrix **Z**_{f} is the matrix that describes the constraints for a congruence or similarity transformation. The coefficient matrix of equation (13) or (14) is used to construct matrix **Z**_{f}.

## 6 Adjustment and Testing

### 6.1 Adjustment

System (28) is a linear system of observation equations and can be solved by least squares. If sufficient points are available in all epochs to determine the transformation parameters, the coefficient matrix is of full rank.

Because of the nonstochastic observations, and because of possible singularities of the covariance matrices of **a**_{1} and the vectors **b**_{i}, the covariance matrix of the observation vector of system (28) is singular. To get a least squares solution of the system, at least five methods are available that make it possible to test nonstochastic observations in the same way as stochastic observations [9]:

- i.The adjustment model is split into two parts for the stochastic and the nonstochastic observations respectively, and a sequential adjustment is applied.
- ii.A switch is made from the model of observation equations to the model of condition equations.
- iii.The covariance matrix is regularised.
- iv.The standard deviations in the covariance matrix that are zero, are replaced by values that are very small.
- v.The observations are orthogonalised and the nonstochastic observations eliminated. A follow-up adjustment determines the test quantities.

Because the system is linearised, iteration is needed to find the least squares solution. To start the iteration good approximate values for all observations and all parameters are needed, which have to satisfy the non-linear equations (10) and the non-linear constraints of section 5.3.2 or 5.3.3. As described by [8], in each iteration step the approximate transformation parameters are updated, using singular value decomposition. Also, in each iteration step, all **b**_{i} and their covariance matrices **D**{**b**_{i}} are transformed with the new approximate transformation parameters to new coordinates **a**_{1} and in the reference system of **a**_{1}.

In each iteration step the approximate values of all observations, and of all parameters have to comply again with the non-linear equations (10) and the non-linear constraints of section 5.3.2 or 5.3.3.

### 6.2 Deformation Testing

If one of the five methods, mentioned in the previous section, is used, standard methods for testing can be applied with the formulas given by [9]. Also the nonstochastic observations can be tested with the same formulas, which means that a method of testing deformation patterns is provided.

If it is not sure whether there is any deformation, or what type of deformation happens, a null hypothesis *H*_{0} is formulated, where no deformation is assumed (∇ is missing in system (28)), and an alternative hypothesis *H*_{a}:

where __∆ y__,

**A**and ∆

**x**are respectively the observation vector, the coefficient matrix and the parameter vector of system (28). In

**A**the last column of the coefficient matrix is missing and in ∆

**x**the parameters ∆∇.

The alternative hypothesis is tested against the null hypothesis, without the need to perform a complete adjustment of (32), by using test statistic *T*_{q} [6, p. 77]:

*q* is the number of columns in ^{2} is the variance factor of unit weight, and __ ê__ the usual least squares residuals and

**Q**

_{y}the cofactor matrix of the observations

__:__

**y**The probability density function of *T*_{q} is a *χ*^{2}-distribution with an expected value of *q*. The test is to choose a significance level *α*, to compute the critical value and to test, whether the computed value of *T*_{q} exceeds the critical value. If this happens, the null hypothesis is rejected [6, p. 78].

### 6.3 S-basis Invariance

In [9] and [10] it is shown that the test statistic of equation (34) is invariant for a change of S-basis of the parameter vector **x**. It is evident from the fact that **x** has been eliminated, and therefore a change of S-basis of **x** doesn’t influence

The test statistic of equation (34) is also invariant for changes of S-bases of the observed coordinate vectors **a**_{1} and **b**_{i}, i = 2, . . ., *p,* if deformation patterns are tested. To see this, model (28) is simplified and reduced. To do this, we assume all observed vectors **a**_{1} and **b**_{i}, and also the parameter vectors **c**_{i} to contain coordinates of the same points in the same order, from which follows:

We also assume stability of all points, and therefore:

and we reduce ∆**c** to a vector **c****I**_{p} defined with unit matrices **I** as:

we get:

This means that the nonstochastic observations **z**_{d} disappear. Furthermore we assume that the nonstochastic observations **z**_{f} are elimated. This can be done by noting that the equation:

means that ∆**f** lies in the nullspace of **Z**_{f}. If **N** is a base matrix that spans this nullspace, we have:

with **f****f** contains 12 parameters of an affine transformation and there are 5 nonstochastic observations to constrain the transformation into a similarity transformation, **f****F****= FN,** it follows that:

So if we use **F****f****F** ∆**f** in model (28), we can omit the nonstochastic observations **z**_{f}.

Because of the stability assumption, no parameters ∇ exist and no matrix **Z**_{∇}.

With (36), (39) and (42), model (28) is written as:

To eliminate **c****H** as:

and the vector **d**, containing the difference vectors of all epoch intervals, as:

Premultiplying equation (43) with **H*** ^{T}*, we get:

This model has the same redundancy as model (43) and yields the same least squares solution.

Let the vectors **a**_{1} and **b**_{i}, *i = *2*, . . .,p* all have been S-transformed to other S-bases. It means that we have new vectors __ b__′:

where *ψ* is the vector of the differential transformations of the coordinate vectors of all epochs. These relate, however, to the same degrees of freedom as the transformations in **f****S** = **F**

A proof for two epochs that test statistic (34) is invariant for changes of S-bases of **a**_{1} and **b**_{2}, by proving that **G**, which is chosen to fulfil:

with *R*(**G**) the complementary space of *R*(**H**^{T}**F****a**_{1} and **b**_{i}, *i* = 2, . . ., *p.*

The conclusion is that if hypotheses concerning deformation patterns are formulated in terms of the original model (28), and they can be reformulated in terms of model (46), which is generally possible, test statistic (34) is invariant for changes to other S-bases of the coordinate vectors **a**_{1} and **b**_{i}, *i* = 2, . . ., *p*.

A deformation hypothesis may concern a point that is part of the S-basis definition and whose coordinates are fixed with a zero standard deviation. No S-transformation is needed to test such a point for deformation. It is demonstrated by the example of [9].

## 7 Experimental Validation

The proposed model can be applied to the 3D monitoring by GPS and total stations of deformations of buildings, harbour quays, bridges, tunnels, land slides, etc. The model gives the possibility to compute statistics and to test hypotheses that describe complex deformation patterns, like the abnormal movement of a subset of points through many epochs, or the periodic oscillation of a subset of points, for example caused by changes of temperature.

To validate experimentally the model, the monitoring of some buildings is taken. To be able to judge effectively the performance of the model, observations have been generated with known standard deviations, to which artificially deformations have been added. Fifteen points have been measured with direction and distance observations from a total station during five epochs. The fifteen points are positioned on three buildings (figure 2), which are monitored because of construction works. The instrument point is not fixed (not monumented). The observations are adjusted using the software package MOVE3 (http://www.move3.com), resulting in *x*, *y*, *z* coordinates and their covariance matrix. The network is not attached to a control network.

A Matlab programme has been written to do the computations. The observations have been generated with the following standard deviations:

- –directions: 0.3 mgon;
- –distances: 1 mm;
- –zenith angles: 0.3 mgon.

The precision with which a point is defined (idealisation precision) is supposed to be 0.5 mm, indicating the precision by which a removable prism can be put on a point.

First no deformation is put in the observations. The adjustment model to test stability of all points is created by adding for each epoch interval, for each point and for each coordinate direction a nonstochastic observation, i. e. 4 × 15 × 3 = 180 nonstochastic observations:

From these equations follows matrix **Z**_{d} in equation (28).

The epochs are joined together with similarity transformations, which are realised for each epoch interval by five constraints on the affine transformation parameters.

The model contains 425 observations (5 × 15 × 3 = 225 coordinates, 4 × 5 = 20 transformation constraints, 180 point constraints), and 273 parameters (225 coordinates, 4 × 12 = 48 transformation parameters), which yields a redundancy of 152. Adjusting the model leads to an overall model test of 0.60. With a critical value of 1.004, based on the use of the B-method of testing with a significance level of a one-dimensional test of 0.1 % and of 47 % for the overall model test, the null hypothesis is accepted.

Then a movement of 1 mm in each epoch interval, in the direction of each coordinate axis of one point (point 101), is added and observations are generated. The same adjustment model as before is used and leads to rejection of the null hypothesis. The movement is then modelled with twelve nonstochastic observations. Assume vector **c** that contains the *x*, *y,* and *z* coordinates of point 101 in epoch *i*. Let vector **a** contain the movements in *x*, *y* and *z* direction between epoch *i* and *j*, and ∆**a** the difference of **a** with its approximate value, necessary for the linearised model. The following 12 nonstochastic observations describe the deformation.

From these equations the matrix

Test of linear point movement over 5 epochs; values are per epoch interval.

Stat. | Est. def. (mm) | |||
---|---|---|---|---|

Pnt. | x | y | z | |

101 | 70.9 | 0.9 | 1.1 | 1.0 |

103 | 16.9 | –0.4 | –0.5 | –0.3 |

102 | 7.0 | –0.1 | –0.1 | –0.5 |

104 | 4.7 | –0.2 | –0.3 | –0.1 |

Minimal detectable deformations (m. d. d.); values are per epoch interval.

M. d. d. (mm) | |||
---|---|---|---|

Pnt. | axis 1 | axis 2 | axis 3 |

101 | 1.55 | 0.80 | 0.76 |

103 | 1.49 | 0.73 | 0.73 |

102 | 1.52 | 0.76 | 0.74 |

104 | 1.48 | 0.72 | 0.72 |

with *λ*_{0} the non-centrality parameter of the *χ*^{2}-distribution and ∇_{0} describing the minimal detectable deformations. They give the deformations that can be detected with the three-dimensional point test of five epochs with a power of 80 %.

Finally five points (101, . . .,105) are given a movement of 1 mm in both the *x* and *y* direction and −0.7 mm in the *z* direction in each epoch interval. It is modelled by 60 nonstochastic observations. Let vector **c** that contains the *x*, *y*, *z* coordinates of the five points in epoch *i*. Let **k** = (1, 1, 1, 1, 1)^{T}, **I**_{3} the (3 × 3)-unit matrix, and **E** = **I**_{3} ⊗ **k**, with ⊗ denoting the kronecker product. From the following nonstochastic observations the matrix **Z**_{∇} is deduced:

with ∆**a** as defined before.

The null hypothesis is rejected again. The test of the hypothesis that the five points have shifted gives a very large test statistic (74.2 with a critical value of 12.6, if the significance level is 0.6 %), indicating that it is a very good hypothesis. The estimated deformation and the minimal detectable deformations are given in table 3. The length of the deformation vector is 1.6 mm, which is exactly the length of the vector that has been put intentionally into the coordinates.

Linear movement of points 101–105 over five epochs; values are per epoch interval.

Est. def. (mm) | M. d. d. (mm) | |||||
---|---|---|---|---|---|---|

Pnt. | x | y | z | axis 1 | axis 2 | axis 3 |

101–105 | 1.0 | 0.7 | –1.1 | 0.93 | 0.71 | 0.60 |

Point movements that are nonlinear in time are modelled by nonstochastic observations that are nonlinear functions of the deformation parameters. To be used in the model, the functions have to be linearised.

If the deformation pattern to be expected is not known, a search has to be performed for the best alternative hypothesis. A strategy is described by [9] for two epochs. Extending it to more than two epochs, one could for example systematically test for a constant linear movement through all epochs of each point individually, of combinations of two points close together, of combinations of three points close together, etc. Because it is not needed to solve a complete adjustment model, only to compute test statistic (34), its degrees of freedom *q* determines the computational burden of testing many hypotheses.

## 8 Conclusions

A model has been built for the adjustment of a time series of 3D coordinates in a geodetic point field. The covariance matrices of the coordinates of all epochs of the time series are used and they may be full and singular. Deformation patterns, or their absence, are modelled as nonstochastic observations. To make the testing of the model invariant for S-transformations, transformations between all epochs are built into the model. The transformations can be similarity or congruence transformations, and are modelled as affine transformations, subject to constraints.

The constraints are implemented as nonstochastic observations. The model is first built as a nonlinear one, and then linearised. The approximate parameter values and their updates in the iteration steps (needed because of the linearisation) have to comply with all nonstochastic observations. For the rotation parameters this is accomplished with singular value decomposition.

In many cases it is a sound deformation analysis procedure to formulate a null hypothesis that assumes no deformation. The nonstochastic observation equations state that the coordinate differences between the epochs are expected to be zero after the transformations. Alternative hypotheses are formulated that describe movements of one or many points over one or many epoch intervals. Standard hypothesis testing is used to test the alternative hypothesis against the null hypothesis. The quality of the tests is described by the sizes of the minimal detectable deformations.

The point movements are formulated as nonstochastic observation equations, which give the matrices to be used in the testing equations.

The model and its adjustment and testing have been verified experimentally with a geodetic network, where 15 points are measured by a total station during five epochs. The results show that 3D deformation analysis of time series of coordinates is possible with the model proposed.

The support of research programme Maps4Society is gratefully acknowledged.

## References

- [1]↑
Caspary, W. and Borutta, H., Robust estimation in deformation models, Surv Rev 29 (1987), 29–45.

- [2]↑
Caspary, W. F., Concepts of Network and Deformation Analysis, Monograph 11, School of Surveying, The University of New South Wales, Australia, Report, 2000.

- [3]↑
Chen, Y.Q., Analysis of Deformation Surveys – A generalized method, University of New Brunswick, Department of Geodesy and Geomatics Engineering, Canada, Report, 1983.

- [4]↑
Heunecke, O., Kuhlmann, H., Welsch, W. M., Eichhorn, A. and Neuner, H., Handbuch Ingenieurgeodäsie: Auswertung geodätischer Überwachungsmessungen, second ed, Wichmann, H, 2013.

- [5]↑
Kami´nski, W. and Nowel, K., Local variance factors in deformation analysis of non-homogenous monitoring networks, Surv Rev 45 (2013), 44–50.

- [8]↑
Velsink, H., Extendable linearised adjustment model for deformation analysis, Survey Review, 47 (2015), 397–410.

- [10]↑
Velsink, H., On the deformation analysis of point fields, Journal of Geodesy 89 (2015), 1071–1087.