Next Article in Journal
An Acceptance Test for Assistive Robots
Next Article in Special Issue
Bearing Damage Detection of a Reinforced Concrete Plate Based on Sensitivity Analysis and Chaotic Moth-Flame-Invasive Weed Optimization
Previous Article in Journal
Low Cost and Compact FMCW 24 GHz Radar Applications for Snowpack and Ice Thickness Measurements
Previous Article in Special Issue
Weighted Structured Sparse Reconstruction-Based Lamb Wave Imaging Exploiting Multipath Edge Reflections in an Isotropic Plate
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation and Prediction of Vertical Deformations of Random Surfaces, Applying the Total Least Squares Collocation Method

by
Zbigniew Wiśniewski
1 and
Waldemar Kamiński
2,*
1
Faculty of Geoengineering, University of Warmia and Mazury in Olsztyn, 10-957 Olsztyn, Poland
2
Faculty of Civil and Environmental Engineering, Gdańsk University of Technology, 80-233 Gdańsk, Poland
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(14), 3913; https://doi.org/10.3390/s20143913
Submission received: 20 June 2020 / Revised: 10 July 2020 / Accepted: 11 July 2020 / Published: 14 July 2020

Abstract

:
This paper proposes a method for determining the vertical deformations treated as random fields. It is assumed that the monitored surfaces are subject not only to deterministic deformations, but also to random fluctuations. Furthermore, the existence of random noise coming from surface’s vibrations is also assumed. Such noise disturbs the deformation’s functional models. Surface monitoring with the use of the geodetic levelling network of a free control network class is carried out. Assuming that, in some cases, the control networks are insufficient in surface’s deformation analysis, additional and non–measurable reference points have been provided. The prediction of these points’ displacements and estimation of the free control network points’ displacement are carried out using the collocation method applying the total least squares adjustment. The proposed theoretical solutions were verified by the simulation methods and on the example of a real control network.

1. Introduction

Engineering structures require monitoring and estimation of their technical conditions. Various systems and measurement techniques are used for this purpose, for example, those constructed with the use of specialized sensors and detectors [1,2]. Geodetic methods used for accurate determination of positions of monitored structure selected points have a significant role in the stability monitoring of buildings, bridges, and dams, among others. Geodetic monitoring is conducted using a variety of techniques and measurement methods. The most common are global navigation satellite systems [3,4,5], photogrammetry, and remote sensing [6], as well as laser scanning [7]. However, the classic geodetic control networks [8,9] still play a basic role in the engineering structures’ deformation analysis with the use of the geodetic methods. Structural deformation is then determined on the basis of changes of points’ positioning (displacements) within this network over time. Such changes are determined on the basis of periodic measurements (measurement epochs).
Geodetic control networks are comprised of controlled points and reference points. Reference points require a stability analysis carried out applying appropriate procedures and statistical tests [10,11]. On the basis of the networks measurements conducted in two or more epochs, the displacements of controlled points and potential reference points are then determined. The displacements of these points can be obtained by determining the differences between the estimators of the networks coordinates in different measurement epochs. The displacements can be also treated as parameters in models of the observations differences [12,13]. In this case, in the literature of the subject, the displacement of individual points of the network is sometimes treated as the shift of the model’s parameters [14,15]. The resignation from the reference points leads to solutions that apply the principles of free adjustment [16,17,18,19]. Such networks belong to the free control networks (FCNs) class. The principles of free adjustment [16,17,19,20] are applied to determine these networks points’ displacements. These principles derive from the free adjustment theory presented in numerous papers [21,22,23,24,25].
It is generally assumed that the deformation of a structure in time is deterministic. There are, however, structures for which this assumption is inadequate, particularly if they are situated in a dynamically changing environment or are a subject to dynamic loads [26]. For example, when analysing the bridges’ deformations, not only deterministic deformations can be considered, but also the deformations caused by vehicles traffic and other factors as well as random strains and vibrations of the examined object [27,28]. This also applies to such structures as river dams, tunnels, port quays, and so on [29,30]. In many cases, the object of interest is the deformation of selected surfaces of engineering structures. The results of the object’s random deformations are the fluctuations of its surface, whereas the result of the vibration is the random noise. In this approach, the monitored surfaces can be treated as random fields represented by appropriate random functions.
In practice, there are situations when the subject of interest is focused on the positioning changes not only of the controlled points (CPs), but also of the additional points that are not directly related to the geodetic network. These points, hereinafter called the extended controlled points (ECPs), are necessary for examining the engineering object condition. However, it is not possible to determine their displacements based on the direct geodetic observations (e.g., covering, lack of access, damage). It is assumed that ECPs are situated in close proximity to actual controlled points. Assuming that the monitored surface is a random field, then the use of statistical interpolation and prediction methods could be proposed to determine ECPs’ displacement. It is important to remember that there are limited possibilities of using ECPs determined in this way. However, in certain cases, the solutions proposed below can be the way to approach this issue (or to complete other methods) in difficult and unconventional situations.
The least-squares collocation (LSC) is the most popular method of interpolating and predicting the values of random functions measured in discrete points. The idea of this method was proposed and developed in the following studies inter alia [31,32,33,34,35]. The LSC was applied in solving various geodetic problems, including, inter alia, determination of gravimetric measurements results [36,37,38] and coordinates transformation [39,40,41,42]. In the paper [43], the use of collocations in the mobile LiDAR (Light Detection and Ranging) system was proposed. Tscherning [44,45] applied the LSC method in analysing the observations gathered during the GOCE (Gravity field and steady-state Ocean Circulation Explorer) mission. The principles of this collocation were also the basis of the analysis of deformations caused by earthquakes [46]. Yang [47] presented a way for the robustness of the collocation method against the observations’ gross errors.
In this paper, it is assumed that FCN is a levelling network. All of its points serve as CPs and are located on a surface that is characterised by deterministic and random deformations. The general LSC principles are applied in determining the vertical displacements of CPs as well as in displacements of non-measurable ECPs. In the paper, it will be presented that surface noise leads to random disturbances of the special matrices occurring in functional observation models. For this reason, the conventional solutions used in LSC will be replaced by solutions proposed in the total least squares (TLS) theory. The theory and basic applications of TLS were discussed, inter alia, in [48,49,50]. The collocation that uses the TLS for solving the optimisation problem is called TLSC (total least-squares collocation). The application of TLSC will allow not only the determination of the deterministic and random displacements of CPs and ECPs, but also the estimation of the value of random disturbance in these points.

2. Basic Assumptions and Models

2.1. Surface Fluctuation

Assume that the surface π is the subject to the time-variable t T vertical displacements. These deformations are expressed by changes in the heights H ( t ) of its points (in the adopted reference system). Assume now that vertical deformations are subject to time- and position-dependent random fluctuations. When fluctuations are described by the random function ξ t ( ω ) , where ω Ω is the position parameter (e.g., ω = ( x , y ) ), then the heights of the points on the surface can be expressed as
( t , ω ) = H ( t ) + ξ t ( ω )
For function ξ t ( ω ) , it is assumed that E { ξ t ( ω ) } = 0 , hence E { ( t , ω ) } = H ( t ) , where E ( ) is an expected value. In the simplest case, the fluctuations can be treated as signals s t ( ω ) with assumed covariance function. In this case, ξ t ( ω ) = s t ( ω ) . This approach is applied in solving physical geodesy problems, for instance, in processes of filtering and predicting the gravimetric anomalies [32,33,38].
When surface π is an element of the engineering object, then not only its deterministic deformations and random fluctuations, but also additional noise caused by structural vibrations [28,29], should be projected. This noise should be marked with e t and called the primary surface noise. Assume that these are random values with an expected value of E ( e t ) = 0 and variance of σ e 2 (for each t T ). Furthermore, similarly to [51], it is assumed that surface noise is Gaussian.
At each point P f π with coordinates ω f , the primary surface noise can be modified (amplified or damped) by signals s t ( ω f ) occurring at these points. Treating standardised signals s t * ( ω f ) = s t ( ω f ) / σ s as noise-modifying functions, secondary noise ϑ t at point P f is presented as
ϑ t ( ω f ) = e t , f s t * ( ω f )
( σ s —signal’s standard deviation identical for each t T ). Equation (2) results in the fact that secondary noise does not exist in two cases: (1) no object vibrations, hence no primary noise; and (2) there are vibrations generating noise, but their secondary version decays in the absence of signals.
Assume now that the secondary noise generated from other points on the surface reaches each specific point of surface π . At point P i π , the total (resultant) noise presented in the following way can be expected:
ϑ ¯ i , t = f = 1 ϑ i , t ( ω f ) = f = 1 e i , t , f s t * ( ω f )
where ϑ i , t ( ω f ) = e i , t , f s t * ( ω f ) is the secondary noise generated at point P f reaching point P i . The resultant of secondary noise ϑ ¯ i , t and the signal s t ( ω i ) form the total fluctuations of the P i point’s height on this surface, that is,
ξ t ( ω i ) = s t ( ω i ) + ϑ ¯ i , t
In practice, the surface vertical deformations are determined on the basis of changes in heights H ( t ) of CPs on the control levelling network (in this paper, it is assumed that this is the FCN). Points P k of this network create a set P = { P 1 , , P r } , whereas their positions are determined by coordinates ω k = ( x k , y k ) included in this set Ω = { ω 1 , , ω r } . When surface points π are limited to the set P , then in accordance with Equation (3), the resultant of secondary noise reaching point P k P (CP) can be presented as
ϑ ¯ k , t = f = 1 r e k , t , f s t * ( ω f ) = e k , t T s t *
where e k , t = [ e k , t , 1 , , e k , t , r ] T , s t * = [ s 1 , t * , , s r , t * ] T , s f , t * = s t * ( ω f ) = s t ( ω f ) / σ s , f = 1 , , r .
In deformation analysis, an interesting issue along with displacement and signals at CPs is also the prediction of displacements and signals at other non-measured points on the surface π . These points (relative to FCN) are referred to as ECPs. These points and their coordinates create the sets P = { P 1 , , P z } and Ω = { ω 1 , , ω z } , respectively. Moreover, assume that s t * ( ω f ) = s t ( ω f ) / σ s is a standardised signal at point P f P . Taking into account the ECPs, it leads to the resultant of secondary noise reaching the measured point P k P in the following, extended form
ϑ ¯ k , t = f = 1 r e k , t , f s t * ( ω f ) + g = 1 z e k , t , g s t * ( ω f ) = e k , t T s t * + ( e k , t ) T ( s t ) *
where e k , t = [ e k , t , 1 , , e k , t , z ] T , ( s t ) * = [ ( s 1 , t ) * , , ( s z , t ) * ] T , ( s g , t ) * = s t * ( ω g ) = s t ( ω g ) / σ s , g = 1 , , z On the basis of Equations (4) and (6), total fluctuation of the heights of point P k P can be presented as
ξ t ( ω k ) = s t ( ω k ) + ϑ ¯ k , t = s t ( ω k ) + e k , t T s t * + ( e k , t ) T ( s t ) *
For the signals s t ( ω f ) , ω f Ω Ω , it is assumed that they form a stationary and isotropic process. In this case, the identical following covariance function for each t T can be assigned to these magnitudes [52]
C ( ω i , ω j ) = C ( ω i , ω i + Δ ω ) = C ( Δ ω ) = C ( d i , j ) = σ s 2 ρ ( d i , j )
where ω i , ω j Ω Ω , Δ ω = ω j ω i , and d i , j = Δ ω 1 / 2 are distance between any random pair of points belonging to the set P P ( ρ ( d i , j ) —correlation function). The following values can be obtained using Function (8): covariances cov ( s k , s l ) = C ( d k , l ) of signals s k and s l at points P k , P l P (CPs), covariances cov ( s q , s p ) = C ( d q , p ) of signals s q and s p at points P q , P p P (ECPs), and covariances cov ( s q , s k ) = C ( d q , k ) among signals at points P q P and signals at points P k P .

2.2. Heights and Displacements Models

Considering the total surface fluctuations, the height of point P k P lying on this surface is presented in the following form (based on Equations (1) and (7))
( t , ω k ) = H k ( t ) + ξ t ( ω k ) = H k ( t ) + s t ( ω k ) + e k , t T s t * + ( e k , t ) T ( s t ) * = H k , t + s k , t + e k , t T s t * + ( e k , t ) T ( s t ) *
where H k , t = H k ( t ) and s k , t = s t ( ω k ) . After introducing a surface noise vector e k , t = [ e k , t T , e k , t T ] T , joint for CPs and ECPs, and their joint vector of standardised signals s t * = [ ( s t * ) T , ( ( s t ) * ) T ] T , Equation (9) can also be written as
( t , ω k ) = H k , t + s k , t + e k , t T s t *
In case of signals’ absence at all CPs and ECPs, model (10) is reduced to its classic form ( t , ω k ) = H k , t . This means that, in this case, even if there are surface vibrations and there is random noise coming from these vibrations, they do not contribute to the height model for point P k . Consequently, similarly to random measurement errors, surface noise will disturb the observation (these magnitudes cannot be separated without additional assumptions). Whereas when f : e k , t , f = 0 ,   e k , t , f = 0 , then surface fluctuation is composed only of signals, this means that ξ t ( ω k ) = s t ( ω k ) takes place at point P k P and
( t , ω k ) = H k , t + s k , t
In classic levelling networks forming the FCN, heights H k , t are not directly measured. They are treated as unknown parameters in functional observations models (the observations of heights differences between pairs of points). In this case, separating the signals from non-random parameters is a problem. Although it is possible to conduct the estimation of random parameters ( t , ω k ) = H k , t + s k , t = H k , t r a n d (e.g., applying the Bayes method, such as in [53,54]), the determination of the height estimator H k , t and the estimation of the signal s k , t related to this height remain an unsolved problem. In the proposed model (10), there is not only a signal at point P k , but there are also signals at other CPs (and in ECPs) that, together with surface noise, form the resultant of secondary noise. Therefore, the direct connection of deterministic parameters and random signals is violated here. This creates an opportunity for separating them in the estimation process. Model (10) can be also applied in estimating the surfaces noise contained in vectors e k , t .
Assume now that
C s = [ C s C s s C s s C s ]
is a covariance matrix of the joint signal vector s t = [ s t T , ( s t ) T ] T determined on the basis of the covariance function under Equation (8), where ( C s ) k , l = cov ( s k , s l ) , ( C s ) q , p = cov ( s q , s p ) , ( C s s ) q , k = cov ( s q , s k ) , and C s , s = C s , s T ( ( C ) i , j i, jth element of the matrix C). If the covariance matrix C s is positively defined, then there is a lower triangular matrix R for which C s = R R T (Cholesky decomposition, for example, [54]). The vector created on its basis
s t * = R 1 s t
has a covariance matrix C s * = R 1 C s ( R 1 ) T = I r + z identical for each t T ( I r + z —unit matrix of r + z ). After Equation (13) is introduced into Equation (10), the result is
( t , ω k ) = H k , t + s k , t + e k , t T R 1 s t
The levelling geodetic networks CPs’ displacements are calculated on the basis of the observations of the heights’ differences between these points, conducted during a few measurement epochs. Assume that the height difference Δ ( k , l ) = ( t , ω l ) ( t , ω k ) between points P k , P l P is the object of measurement in epochs t α , t β T . On the basis of Equation (14), the heights of these points in epoch t α can be written as
( t α , ω k ) = H k , α + s k , α + e k , α T R 1 s α ,   ( t α , ω l ) = H l , α + s l , α + e l , α T R 1 s α
In this case, the heights’ differences observation model Δ ( k , l ) is formed as
h i , α = h ( k , l ) , α = Δ ( k , l ) α + v i , α = ( t α , ω l ) ( t α , ω k ) + v i , α = H l , α H k , α + s l , α s k , α + e l , α T R 1 s α e k , α T R 1 s α + v i , α = H l , α H k , α + s l , α s k , α + ε i , α T R 1 s α + v i , α
where ε i , α = e l , α e k , α , i = ( k , l ) . The random error of observation h i , α is marked with v i , α (with expected value E ( v i , α ) = 0 and variance σ v 2 ). The elements of vector ε i , α are differential noise with identical variances σ E 2 = 2 σ e 2 . The magnitudes of H k , α and H l , α , k , l = 1 , , r are the parameters of the functional model (16), determined by estimation. For total FCN, these parameters form vector X α = [ H 1 , α , , H r , α ] T . After using this vector, the model of ith observation can be presented as
h i , α = a i X α + a i s α + ε i , α T R 1 s α + v i , α
where a i = [ 0 , , 0 , 1 k , 0 , , 0 , 1 l , 0 , , 0 ] . When surface vibrations are absent, then ε i , α = 0 , hence
h i , α = a i X α + a i s α + v i , α = a i ( X α + s α ) + v i , α = a i X α r a n d + v i , α
where X α r a n d = X α + s α . In this specific case, the random signals are combined with parameters, which consequently leads to a random parameter X α r a n d with a covariance matrix of C X α r a n d = C s .
Until this point, the signals in ECPs were treated only as the components of secondary noise reaching individual points within the FCN. Assume that these magnitudes, that is, signals at points P 1 , , P z P , will also be determined by estimation process. Then, according to the principles of collocation [32,33,34], model (17) can be supplemented to the following form
h i , α = a i X α + a i s α + ε i , α T R 1 s α + 0 s α ( ω 1 ) + + 0 s α ( ω z ) + v i , α = a i X α + a i s α + ε i , α T R 1 s α + 0 z T s α + v i , α
where 0 z = [ 0 1 , , 0 z ] T is the zero vector (the zero component 0 z T s α added to model (17) does not change the primary sense of this model). After further transformation of Equation (19), the formula is as follows
h i , α = a i X α + a i s α + 0 z T s α + ε i , α T R 1 s α + v i , α = a i X α + [ a i , 0 z T ] [ s α s α ] + ε i , α T R 1 s α + v i , α = a i X α + b i s α + ε i , α T R 1 s α + v i , α = a i X α + ( b i + ε i , α T R 1 ) s α + v i , α
where b i = [ a i , 0 z T ] . For total FCN, that is, including all observations h i , α , i = ( k , l ) = 1 , , n , based on Equation (20), the following model is obtained
y α = A X α + ( B + E α R 1 ) s α + v α = A X α + ( B + E α ) s α + v α
where y α = [ h 1 , α , , h n , α ] T is the observations vector and v α = [ v 1 , α , , v n , α ] T is the vector of random observation errors with a covariance matrix of C v = σ v 2 I n . Matrix A = [ a 1 T , , a n T ] T is a classic coefficient matrix in functional models of levelling networks. What is notable, in a levelling FCN, is that the A matrix is the matrix of vertically incomplete rank, that is, with rank ( A ) = u = r d = r 1 ( d = 1 —deficient rank). Furthermore, B = [ A , 0 n , z ] ( 0 n , z —zero matrix of n × z dimensions) and E α = [ ε 1 , α , , ε n , α ] T = [ E α , E α ] . The E α matrix is a matrix of random noise as a result of surface vibrations, whereas E α = [ ε 1 , α , , ε n , α ] T applies to noise at CPs and E α = [ ε 1 , α , , ε n , α ] T is the noise matrix at ECPs. All elements of these matrices have common variance σ E 2 = 2 σ e 2 . Matrix E α = E α R 1 = [ E α , E α ] , where E α = E α R 1 , E α = E α R 1 , is a matrix that disturbs the corresponding blocks of B matrix ( E α disturbs the A matrix, whereas E α disturbs the zero matrix 0 n , z ).
Similarly to above, the determination of the model of observation vector y β in epoch t β , that is, based on Equation (21), gives
y β = A X β + ( B + E β ) s β + v β
(index α is substituted here with index β indicating the current measurement epoch).
Assume that the network points’ displacements are determined on the basis of difference δ y = y β y α of observation vectors in epochs t α and t β . Using Equations (21) and (22), the following model of this difference is obtained
δ y = A u + ( B + E β ) s β ( B + E α ) s α + υ
where u = X β X α is the CPs’ vertical displacement in these epochs. The magnitude u will be also treated as the deterministic parameter shift (e.g., [14]). Vector υ = v β v α is the vector of the differences of the observations errors with a covariance matrix of C υ = 2 σ v 2 I n .
In model (23), it is assumed that disturbances of the B matrix differ in epochs t α and t β . In practice, however, it is possible to adopt the assumption that object vibrations in both of these epochs are similar. In this case, it can be expected that primary noise e t will have similar values coming from these vibrations. It should be noted that, even with this simplification, fluctuations of heights of the points in epochs t α and t β can differ from one another. Signals s t ( ω f ) modifying primary noise (see Equation (2)) remain dependent on time. The assumptions made are not general in nature (i.e., they do not correspond to each practical situation), but they can significantly simplify model (23). Without this simplification, it may be difficult to estimate surface noise. Applying equation E = E α = E β , and thus adopting identical disturbs of B matrix in both measurement epochs, that is, E = E α = E β , the model of the vector of the observations’ differences of heights’ differences in epochs t α and t β can be expressed as
δ y = A u + ( B + E ) ( s β s α ) + υ = A u + ( B + E ) η + υ
Vector η = s β s α = [ η T , ( η ) T ] T , where η = s β s α , η = s β s α is a vector of differential signals with zero expected values and covariance matrix of C η = 2 C s . Its elements can be also treated as signals for points‘ displacements in epochs t α and t β . Vector η is thus the vector of CPs’ random displacement, which can be used to determine the total displacements u η = u + η of these points, in conjunction with deterministic displacements u . Meanwhile, the η vector is the vector of predictions of random ECPs’ displacements. The deterministic displacements of these points are not parameters in model (24), hence they are not determined by estimation process. However, if it is assumed that u E C P is the vector of interpolated ECPs’ displacements (calculated on the basis of the displacements of CPs located in their proximity), then their total displacements can be expressed as u η , E C P = u E C P + η .
For the purpose of further optimisation process, the E matrix can be conveniently substituted with a vector created from its subsequent columns, that is, with vector
v E = vec ( E ) = [ vec ( E ) T , vec ( E ) T ] T
with a covariance matrix of
C v E = σ E 2 I n ( r + z ) = σ E 2 ( I n + z I n )
( Kronecker product). Applying the following, general dependence (e.g., [49,51,55])
vec ( A B C ) = ( C T A ) vec ( B )
the expression can be converted to E η = E R 1 η as
E η = E R 1 η = I n E ( R 1 η ) = ( η T R 1 I n ) v E
Introducing the above transformation into model (24), the result is as follows
δ y = A u + ( B + E ) η + υ = A u + B η + ( η T R 1 I n ) v E + υ

3. Optimisation Problem and Its Solution

The unknown elements of model (29) can be determined according to the principles of the TLS method (e.g., [49]). The following optimisation criterion is adopted
φ ( υ , v E , η , u , λ ) = υ T C υ 1 υ + v E T C v E 1 v E + η T C η 1 η 2 λ T ( W ( υ , v E , η , u ) ) = min
where
W ( υ , v E , η , u ) = A u + ( B + E ) η + υ δ y = A u + B η + ( η T R 1 I n ) v E + υ δ y = 0
is the equation that links the determined magnitudes ( λ —vector of Lagrange multipliers).
The optimisation problem φ ( υ , v E , η , u , λ ) = min is solved by such magnitudes as υ ^ , v ^ E , η ^ , u ^ , for which the following Euler–Lagrange conditions are fulfilled (e.g., [48,49,50])
φ υ ^ | υ ^ , v ^ E , η ^ , u ^ , λ ^ = 2 C υ 1 υ ^ 2 λ ^ = 0
φ v ^ E | υ ^ , v ^ E , η ^ , u ^ , λ ^ = 2 C v E 1 v ^ E 2 ( R 1 η ^ I n ) λ ^ = 0
φ η | υ ^ , v ^ E , η ^ , u ^ , λ ^ = 2 C η 1 η ^ 2 ( B + E ^ ) T λ ^ = 0
φ u | υ ^ , v ^ E , η ^ , u ^ , λ ^ = 2 A T λ ^ = 0
φ λ | υ ^ , v ^ E , η ^ , u ^ , λ ^ = A u ^ + ( B + E ^ ) η ^ + υ ^ δ y = 0
Furthermore, the necessary conditions used to positively define second derivatives must be fulfilled here
2 φ υ ^ υ ^ T = 2 C υ 1 , 2 φ v E v E T = 2 C v E 1 , 2 φ η η T = 2 C η 1
The following residual vector is a result of Equation (32) and concerns the vector of observation differences in epochs t α and t β
υ ^ = C υ λ ^
Meanwhile, applying Equation (33) and taking into consideration the covariance matrix (26), the following equation will be obtained
v ^ E = C v E ( R 1 η ^ I n ) λ ^ = σ E 2 ( I n + z I n ) ( R 1 η ^ I n ) λ ^ = σ E 2 ( R 1 η ^ I n ) λ ^
Because v ^ E = vec ( E ^ ) and λ ^ = vec ( λ ^ ) , and based on Equation (27), the following matrix is determined
E ^ = σ E 2 λ ^ η ^ T ( R 1 ) T = [ E ^ , E ^ ]
which constitutes the assessment of the noise matrix in CPs ( E ^ ) and in ECPs ( E ^ ). Therefore, the matrix that disturbs corresponding matrix blocks B = [ A , 0 n , z ] is
E ^ = E ^ R 1 = σ E 2 λ ^ η ^ T ( R 1 ) T R 1 = σ E 2 λ ^ η ^ T C s 1 = [ E ^ , E ^ ]
Using Equation (35), the following vector can be determined
η ^ = C η ( B + E ^ ) T λ ^ = [ η ^ T , ( η ^ ) T ] T
Its corresponding blocks are the estimations of CPs’ random displacements (vector η ^ ) and prediction of ECPs’ random displacements (vector η ^ ). Inserting Equations (38) and (42) into Equation (36), the result is
( B + E ^ ) C η ( B + E ^ ) T λ ^ + C υ λ ^ + A u ^ δ y = D λ ^ + A u ^ δ y = 0
where
D = ( B + E ^ ) C η ( B + E ^ ) T + C υ
If rank ( D ) = n , then the solution of Equation (43) is the following vector of Lagrange multipliers
λ ^ = D 1 ( A u ^ δ y )
After inserting Equation (45) into Equation (35), the following equation is obtained
A T D 1 A u ^ A T D 1 δ y = 0
which can be used to determine the estimator u ^ of the vector of deterministic displacements u .
Because rank ( A T D 1 A ) = rank ( A ) = r 1 , then any of the g-inverse of matrix A T D 1 A , should be used to solve Equation (46). In free networks adjustment, the following g-inverse of minimum norm, is particularly significant (e.g., [21,22,23,24,25,56])
( A T D 1 A ) + = A T D 1 A ( A T D 1 A A T D 1 A )
The estimator u ^ determined with the use of this inverse fulfils the optimisation criterion u 2 = u T u = min (other criteria are also used in similar problems, see, for example, [17]). Because the rank of A matrix is known, the method based on the presentation of this matrix in form of block A = [ A a , A b ] can be used to calculate ( A T D 1 A ) + . The A a matrix is a matrix with r 1 columns, for which rank ( A a ) = rank ( A ) , whereas A b is a single-column matrix (e.g., [14,25,57]). Then,
( A T D 1 A ) + = [ A a T D 1 A a ( G G T ) 1 0 r 1 , 1 A b T D 1 A a ( G G T ) 1 0 ]
where G = [ A a T D 1 A a A a T D 1 A b ] . After applying Equation (48), the solution of Equation (46) is the following estimator of CPs’ deterministic displacements
u ^ = ( A T D 1 A ) + A T D 1 δ y = G T ( G G T ) 1 A a T D 1 δ y
The expressions obtained above indicate that the solution to the problem φ ( υ , v E , η , u , λ ) = min is non-linear. It is required to use an iterative procedure with a properly selected starting step and a rational stopping criterion. Owing to Equations (40) and (42), the starting matrix of random noise E plays a primary role in initiating iteration. Assume that E ( 0 ) = 0 is this matrix. Then, E ( 0 ) = E ( 0 ) R 1 = 0 and
η ( 0 ) = C η ( B + E ( 0 ) ) T λ = C η B T λ = C η [ A , 0 n , z ] T λ
Because A T λ = 0 (see Equation (35)) and 0 n , z T λ = 0 , then it is obtained that η ( 0 ) = 0 for matrix E ( 0 ) = 0 . This means that, for this initial noise matrix, the iteration process will not start. A different, non-zero matrix E ( 0 ) should be thus adopted for the initiation of this process. On the other hand, there can be cases where there is no object’s noise, in that case, E = 0 . For this reason, the authors propose that the elements of the starting matrix E ( 0 ) should be small enough not to disturb the structure of observation (in the case of the absence of noise), and large enough to ensure that the iteration process is able to start (if noise exists). The authors propose that, according to the adopted nature of the noise matrix E , matrix E ( 0 ) should be simulated with the use of a generator of Gaussian’s random matrices with a standard deviation common to all of its elements σ E ( 0 ) = c σ E . The c coefficient should be selected in such way that E ( 0 ) fulfils the aforementioned requirements (e.g., c = 0.01 ).
For the adopted starting noise matrix and starting disturbance matrix E ( 0 ) = E ( 0 ) R 1 , the following iterative algorithm ( i = 0 , 1 , 2 , , m ) is proposed:
D ( i + 1 ) = ( B + E ( i ) ) C η ( B + E ( i ) ) T + C υ
G ( i + 1 ) = [ A a T ( D ( i + 1 ) ) 1 A a A a T ( D ( i + 1 ) ) 1 A b ]
u ( i + 1 ) = ( G ( i + 1 ) ) T ( G ( i + 1 ) G ( i + 1 ) T ) 1 A a T ( D ( i + 1 ) ) 1 δ y
λ ( i + 1 ) = ( D ( i + 1 ) ) 1 ( A u ( i + 1 ) δ y )
υ ( i + 1 ) = C υ λ ( i + 1 )
η ( i + 1 ) = C η ( B + E ( i ) ) T λ ( i + 1 )
E ( i + 1 ) = σ E 2 λ ( i + 1 ) ( η ( i + 1 ) ) T
E i + 1 = E ( i + 1 ) R 1
Iteration ends for such m that u ( m ) u ( m 1 ) < ξ u , λ ( m ) λ ( m 1 ) < ξ λ , η ( m ) η ( m 1 ) < ξ η and E ( m ) E ( m 1 ) < ξ E (for given ξ u , ξ λ , ξ η , ξ E ).

4. Numerical Tests

4.1. Simulated Levelling Network

Basic numerical tests were carried out on the example of simulated levelling network (Figure 1) constituting the FCN. This network is composed of a set of P = { P 1 , , P r } , r = 25 , CPs and a set of P = { P 1 , , P z } , z = 5 , ECPs. The positions of these points are defined by their coordinates, ω k = ( x k , y k ) Ω and ω l = ( x l , y l ) Ω (for k = 1 , , 25 and l = 26 , , 30 ), respectively. It was assumed that the network is measured in two epochs, t α and t β . In each of these epochs, the observations of heights differences h i , i = 1 , , 40 between CPs were simulated by adding the observations errors v i of normal distribution N ( 0 , σ v 2 ) to their theoretical values, where σ v is the assumed standard deviation of these errors. These errors were generated using a random number generator σ v r a n d n ( n , 1 ) , n = 40 included in the MatLab package. Observation vectors y α = [ h 1 , α , , h 40 , α ] T and y β = [ h 1 , β , , h 40 , β ] T were obtained this way and then used to calculate the vector of observations differences δ y = y β y α in epochs t α and t β .
The simulation of differential signals in CPs and in ECPs, that is, the random displacements of these points, required the application of the following procedure:
  • Generation of mutually independent elements of vectors, s t * = [ ( s t * ) T , ( ( s t ) * ) T ] T ~ N ( 0 , I 30 ) , t = t α , t β using generator r a n d n ( r + z , 1 ) .
  • Creation of a covariance matrix C s of signals vector s t = [ s t T , ( s t ) T ] T for the adopted covariance function C ( d i , j ) = σ s 2 ρ ( d i , j ) , i , j = 1 , , 30 .
  • Calculation of signals vectors s t = R s t * ~ N ( 0 , C s ) , t = t α , t β , based on the R matrix of C ( d i , j ) = σ s 2 ρ ( d i , j ) = σ s 2 exp ( k d i , j 2 ) , C s = R R T distribution and simulated vectors s t * .
  • Calculation of simulated random displacements vector of η = ( s β s α ~ N ( 0 , 2 C s ) .
Matrix of random noise resulting from surface vibrations, that is, the E = [ E , E ] matrix was simulated using generator σ E r a n d n ( n , r + z ) . This matrix is the basis for the calculation of the E = E R 1 matrix that disturbs the B matrix in model (24).
The Gaussian’s covariance function of signals was adopted in the tests [41,42,54,58] in the following form
C ( d i , j ) = σ s 2 ρ ( d i , j ) = σ s 2 exp ( k d i , j 2 )
The k coefficient in the correlation function ρ ( d i , j ) = exp ( k d i , j 2 ) was determined in such way that, for the most distant points of network d max = 230   m , the correlation coefficient between the signals at these points has the assumed minimum value ρ min . The value of this coefficient must be selected in such way that the covariance matrix C s is positively defined (too high a value of ρ min leads to the negative definitiveness). For instance, k = ln ( ρ min ) / d max 2 = 8.7 10 5 is obtained for ρ min = 0.01 . There are also functions that belong to other classes of theoretical covariance functions and can be applied in practice. The chosen function should, however, correspond to the unique engineering object’s specification as well as its monitoring programme and technical conditions.
The influence of not recognising random displacements and surface noise on estimators of the u parameter in functional observation models was verified. For this reason, an analysis applying the Crude version of the Monte Carlo method was carried out. For N = 5000 of independent simulations of vectors v α , s i m , v β , s i m , s α , s i m , and s β , s i m and matrices E s i m , the theoretical vector of differential observations δ ¯ y = A u ¯ was loaded with a random vector ξ s i m i = ( B + E i ) η s i m i + υ s i m i , where E i = E s i m i R 1 , η s i m i = s β , s i m i s α , s i m i , and υ s i m i = v β , s i m i v α , s i m i (for each i = 1 , , N ). The u ¯ vector is the adopted theoretical displacement of CPs. LS-estimators u ^ i of u parameter were determined on the basis of the model δ y = A u + υ , for the simulated observations δ y , s i m i = A u ¯ + ξ s i m i . It should be noted that the model adopted in this paper does not take into account the signals and noise, although observations tend to be burdened by them. For comparison, the u ^ p i LS-estimators determined on the basis of observations burdened only with the vector of differential observation errors υ s i m i (without signals and noise) were also calculated. The estimators u ^ i and u ^ p i , i = 1 , , N , were the basis for calculating the magnitudes of
J u ^ M C = i = 1 N ( u ¯ u ^ i 2 / r ) 1 / 2 / N , J u ^ p M C = i = 1 N ( u ¯ u ^ p i 2 / r ) 1 / 2 / N
which constitute the Monte Carlo assessment of the following comparative parameters
J u ^ = ( u ¯ u ^ i 2 / r ) 1 / 2 , J u ^ p = ( u ¯ u ^ p i 2 / r ) 1 / 2
Calculations were carried out for three versions of the vector u ¯
u ¯ = [ u ¯ 1 ,   u ¯ 2 ,   u ¯ 3 ,   ,   u ¯ 25 ] T = { [ 0 ,   0 ,   0 ,     ,   0 ] T mm [ 5 ,   0 ,   0 ,     ,   0 ] T mm [ 10 ,   0 ,   0 ,     ,   0 ] T mm
Estimators u ^ i = [ u ^ 1 i ,   ,   u ^ 25 i ] T and u ^ p i = [ u ^ p , 1 i ,   ,   u ^ p , 25 i ] T were determined for different values of signal’s standard deviations (the same for each measurement epoch): σ s = 0   mm (no signals), σ s = 0.5   mm , and σ s = 1.0   mm , and for different values of surface noise’s standard deviations: σ E = 0   mm (no noise), σ E = 0.1   mm , σ E = 0.2   mm , and σ E = 0.3   mm . In each case, an identical standard deviation of observations errors σ v = 0.3   mm was adopted. The values of J u ^ M C and J u ^ p M C obtained using the Monte Carlo (MC) method are presented in Table 1 (all values are expressed in mm). Figure 2 presents the courses of estimators u ^ k i and u ^ p , k i obtained in the sequence i = 1 , , N of simulations (for variants u ¯ 1 = 10   mm , σ v = 0.3   mm , σ s = 1.0   mm , and σ E = 0.1   mm ).
The obtained results indicate that the differential signals (random displacement) omitted in the functional model can significantly affect the displacements values of the control network points that are determined with the use of this model. According to the theoretical assumptions, this influence grows together with the increase of noise values produced by surface vibrations. For instance, for noise with standard deviations σ E = 0.1   mm and σ E = 0.3   mm , the difference between comparative parameters J u ^ M C and J u ^ p M C will be 0.60   mm and 1.18   mm (for the variant, where u ¯ 1 = 10   mm and σ s = 1.0   mm ), respectively. What is noteworthy is that the influence of signals is present also when there is no noise. In the discussed variant, J u ^ M C J u ^ p M C = 0.49 was obtained for σ E = 0   mm . However, in this and similar examples, there is an issue with separating random signals from deterministic displacements (which was discussed in the theoretical part of the paper).
Both differential signals and surface noise are random values that are difficult to foresee (similarly to observations errors). In practice, it is thus impossible to eliminate their influence on the vector of differential observations δ y , and thus it is impossible to actually determine the u ^ p estimator. Extending the observation model to the form of Equation (24) and applying the proposed TLSC method lead to different solutions. In this case, the calculation results can go beyond TLSC estimators of the u deterministic displacements, but also include estimations of the η random displacements (with their prediction for ECPs), and estimations of the E noise matrices. The example of this solution is presented below.
The example assumes the same simulated FCN as above. The simulation was carried out for the following vector of CPs’ theoretical displacements of u ¯ = [ u ¯ 1 ,   u ¯ 2 ,   u ¯ 3 ,   ,   u ¯ 25 ] T = [ 5 ,   0 ,   0 ,     ,   0 ] T mm . The following values σ v = 0.3   mm , σ s = 1.0   mm , and σ E = 0.1   mm were adopted for each of the measurement epochs t α and t β . Similarly to the previous analysis, the correlation coefficient between the signals in the most distant points is ρ min = 0.01 . The vector of observations differences δ y , s i m = A u ¯ + ξ s i m , where ξ s i m = ( B + E s i m ) η s i m + υ s i m , was determined in the same way. The vector was the basis for iterative calculations performed with the use of algorithm defined by Equations (51)–(58). The starting matrix of random noise E ( 0 ) was generated by adopting σ E ( 0 ) = 0.01 ,   σ E = 0.001   mm . The objects of these calculations are as follows: the u ^ estimator of CPs’ deterministic displacements, the η ^ estimation of these points’ random displacements, the η ^ prediction of ECPs random displacements, the total CPs displacement u ^ η = u ^ + η ^ , prediction u ^ η = u ^ + η ^ of ECPs’ total displacements, and estimation E ^ of the noise matrix. By comparison, the classic least squares estimator u ^ L S of the u parameter was also determined. Model δ y = A u + υ was the basis for determining estimator u ^ L S and was used for the same observations’ differences vector δ y , s i m as above. Simulated and determined values are presented in Table 2. The table also presents the values of the following comparative parameters
J u ^ = ( u ¯ u ^ 2 / r ) 1 / 2 , J u ^ L S = ( u ¯ u ^ L S 2 / r ) 1 / 2 , J η ^ = ( η η 2 / r ) 1 / 2 , J u ^ η = ( ( u ¯ η u ^ η ) 2 / r ) 1 / 2
where u ¯ η = u ¯ + η s i m is the simulated total displacement.
The random displacements’ prediction η ^ allows for prediction of ECPs’ total displacements. However, an interpolation of deterministic displacements u of these points must be carried out first. Assume that point P i P lies between points P k P l P , and for these points, the displacements estimators u ^ k and u ^ l have been determined. Then, in the simplest case, the linear interpolation expressed by the following formula can be used to determine the displacement u ^ i
u ^ i = u ^ k + u ^ l u ^ k d k , l d k , i
where d k , l is the distance between points P k and P l , whereas d k , i is the distance between point P k and point P i . Interpolated displacements of points P 1 , P 5 (the u ^ vector) and total displacements of these points (vector u ^ η = u ^ + η ^ ) determined using Equation (53) are also presented in Table 2. Final information about random, deterministic, and total displacements of all points (CPs and ECPs) within a simulated free control network is presented in Figure 3.
While analysing the obtained results, attention should be paid to the values of comparative parameters J u ^ L S = 0.92   mm and J u ^ = 0.79   mm . These values mean that the deterministic displacement estimators designated using this method are generally closer to the theoretical values of these displacements, compared with corresponding LS estimators. This particularly applies to displacement u 1 , for which the absolute values of differences between its relatively high theoretical value u ¯ 1 = 5   mm and the designated estimators are | u ¯ 1 u ^ 1 , L S | = 0.57   mm and | u ¯ 1 u ^ 1 | = 0.05   mm , respectively. The theoretical displacements of the remaining points are equal to zero. Assume that such displacements are the u 0 parameter that is common to these points. Then, LS and TLSC estimators of parameter u 0 can be treated as the realisations u ^ 0 , k , L S and u ^ 0 , k , k = 1 , , 25 of two random variables, u ^ 0 , L S and u ^ 0 , respectively. In this approach, the comparisons regarding the quality of identification of zero displacement can be made on the basis of the determined estimators’ histograms (Figure 4). These histograms indicate that, compared with LS estimators, TLSC estimators are more focused around the theoretical value u ¯ 0 = 0   mm .
Taking into account the issue of determining the network points’ displacements, it is interesting that total displacement estimators agree with their simulated values. This is displayed by parameter J u ^ η = 0.55   mm , which differs insignificantly from parameter J η ^ = 0.62   mm that defines the quality of random displacements estimation. These estimations are generally similar to the simulated values (Figure 5a), whereas this similarity is also applicable for predictions at ECPs.
Owing to the random character of theoretical signals, a better match can hardly be expected here, particularly with reference to specific, individual signals. Attention should also be paid to the random noise estimations at CPs and ECPs, presented in Figure 5b. These general estimations have lower values compared with simulated noise. This is primarily caused by the adopted, relatively small standard deviation of these values. In this situation, although the adopted model should prevent them, some noise becomes part of observations errors. The consequence of this autonomous integration is that the assessments of the E matrix elements are underestimated. From the theoretical point of view, this could be prevented by adopting a higher value for σ E , for example, comparable with the value of the observation’s standard deviation. The conducted numerical tests have, however, indicated that the iteration under Equations (51)–(58) is not convergent in many cases. Substituting the proposed, convenient iteration procedure with a different procedure in relation to such cases is an important issue, which is, however, beyond the scope of this paper. In the presented example, the iteration process for end criteria ξ u = ξ η = ξ E = 0.001   mm and ξ λ = 0.001   mm ended after m = 10 steps. The values of norms u ( i ) u ( i 1 ) , λ ( i ) λ ( i 1 ) , η ( i ) η ( i 1 ) , and E ( i ) E ( i 1 ) obtained in individual steps are presented in Figure 6.

4.2. Real Free Control Network

The simulations presented in the previous section have been supplemented with the results of real FCN calculations. The network was established to monitor the loading quays (Figure 7a). It is comprised of CPs P 1 , , P 14 , among which heights differences were measured (Figure 7b). In difficult to measure but critical points of the quay, points P 1 , , P 4 were created, to be treated as ECPs.
The calculations were performed for two measurement epochs: the year t α = 1998 and the year t β = 2008 . The specific nature of operations carried out on port quays justifies the application of the proposed method. The ongoing operation of gantries and other loading equipment, the stacking and lifting containers, and the docking of vessels can cause vibrations and elastic displacements. In this case, not only the object’s deterministic deformations, but also random fluctuations consisting of signals and surface noise can be expected.
Observations of heights differences h ( k , l ) among points P k ,   P l , k = 1 , , 14 , l = 2 , , 14 , 1 and differences δ y i = h ( k , l ) β h ( k , l ) α , i = ( k , l ) among these heights differences, obtained in both measurement epochs, are presented in Table 3. In each measurement epoch, observations were obtained assuming standard deviation σ v = 0.2   mm . In practice, the selection of adequate signals values and surface noise’s standard deviation is an important issue. In this example, calculations were performed in two versions: for signal’s standard deviation, clearly exceeding observations’ standard deviation ( σ s = 3 σ ν = 0.6   mm ), and when these standard deviations are equal to each other ( σ s = σ ν = 0.2   mm ). In both versions, identical standard deviations were adopted for surface noise σ E =   0.1   mm . The subject of the calculations was the deterministic, random, and total CPs’ displacements (vectors u ^ , η ^ , and u ^ η = u ^ + η ^ , respectively). For comparative reasons, a classic estimator of displacement u ^ L S for these points was also determined. The calculations also produced the prediction η ^ of random displacements regarding ECPs (i.e., points P 1 , , P 4 ). The η ^ vector and interpolated displacements u ^ of ECPs were then used to determine their total displacements u ^ η = u ^ + η ^ . The results are shown in Table 4. A graphic presentation of determined random displacements and surface noise is presented in Figure 8a,b. The iteration process leading to the presented results ends after eight steps. The course of the most important norms u ( i ) u ( i 1 ) and η ( i ) η ( i 1 ) essential to this process is presented in Figure 9.
The TLSC estimators of deterministic displacements for both versions of signals’ standard deviations vary from each other. This is shown by the ( u ^ 0.6 u ^ 0.2 2 / r ) 1 / 2 = 0.55   mm parameter, where u ^ 0.2 and u ^ 0.6 are heights’ difference estimators obtained for signals’ standard deviations, σ s = 0.2   mm and σ s = 0.6   mm , respectively. The estimations of random displacements, for which ( η ^ 0.6 η ^ 0.2 2 / r ) 1 / 2 = 0.58   mm is obtained, vary in a similar way. With reference to real FCNs, it is particularly important to compare TLSC estimators of deterministic displacements with their classic LS estimators. Calculating relevant comparative indicators, the results are as follows: ( u ^ 0.2 u ^ L S 2 / r ) 1 / 2 = 0.28   mm and ( u ^ 0.6 u ^ L S 2 / r ) 1 / 2 = 0.83   mm . For the variant where σ s = 0.6   mm , the difference between TLSC and LS estimators of deterministic displacements can have a practical significance. However, paying attention to the total displacements of the analysed network points, for these displacements, where σ s = 0.6   mm ( u ^ η , 0.6 u ^ L S 2 / r ) 1 / 2 = 0.28   mm is obtained. This parameter shows that the total displacement estimations obtained using the TLSC method hardly differ from the LS estimators of deterministic displacements (for σ s = 0.2   mm , this parameter is 0.09   mm ). This means that the displacements that were treated as deterministic in the LS method have been split into deterministic and random in the TLSC method.
Additionally, in this practical example, the prediction of ECPs’ random displacements (for the version in which σ s = 0.6   mm ) was analysed. For instance, point P 1 lies between points P 4 and P 5 at similar distances, d P 1 P 4 = 23   m and d P 1 P 5 = 26   m . The random displacements recorded for these points are η ^ P 4 = 1.99   mm and η ^ P 5 = 0.24   mm . At point P 1 , the random displacement prediction determined using the TLSC method is η ^ P 1 = 0.83   mm , which is the approximate average of the displacements recorded in proximate points (this average is 0.875   mm ). At this moment, it should be noted that the closer CP has a bigger (but still unclear) influence on the random displacement’s prediction at ECP. A different ECP, that is, point P 4 , lies at distances from points P 10 and P 11 of d P 4 P 10 = 35   m and d P 4 P 11 = 69   m . For these points, the assessments of random vertical displacements are η ^ P 10 = 1.03   mm and η ^ P 11 = 0.02   mm . The determined prediction of displacement of point P 4 is equal to η ^ P 4 = 0.81   mm . In this case, the impact of closer CP on the displacement prediction at ECP is definitely larger.

5. Conclusions

In this paper, the free control network was limited to the class of levelling networks. Such networks play the important role in vertical surface deformations analyses. Splitting the points’ displacements within an FCN into deterministic and random parts allows for expansion of the analysis, primarily in terms of the character of monitored surface deformations. Identification of this aspect can be of high practical importance. The different object’s reactions to deterministic deformations as well as to random surface fluctuations can be expected. There are different ways of interpreting these deformations and methods of preventing their effects.
The total displacement estimations obtained using the TLSC method hardly differ from the LS estimators of deterministic displacements. When these methods are applied, the general conclusions pertaining to network deformations can be similar. Therefore, if the character of deformations is not the object of interest, TLSC may be treated as an alternative to the LS method. The proposed method, however, opens the path to more detailed interpretation of displacements, and to these magnitudes’ prediction, for additional points of the monitored surface that are not connected to the control network.
The solution of the optimisation problem in TLSC is iterative. The iterative process is generally convergent and ended after 8–10 steps in the tests conducted as part of the study. The problem with the convergence of the iterative process occurs only when the surface noise’s standard deviations are too large (equal to or exceeding the values of observations’ standard deviations). Selecting the surface noise’s standard deviations, adequate to the specific, monitored object, as well as signals’ standard deviations, is a separate problem requiring additional analyses. This also applies to the selection of adequate covariance functions and their parameters. In this study, the role of empirical examples and analyses was limited to demonstrating the results that could be expected after applying the proposed method. For this reason, the most popular Gaussian’s covariance function using the parameter coming from the adopted, minimum value of correlation coefficient was applied. In specific instances, this function can be replaced by a different, more justified covariance function.

Author Contributions

Conceptualization, Z.W. and W.K.; Methodology, Z.W. and W.K.; Software, Z.W.; Validation, W.K.; Formal analysis, Z.W.; Investigation, Z.W.; Resources, W.K.; Writing—original draft preparation, Z.W.; Writing—review and editing, W.K.; Visualization, Z.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Institute of Geodesy, University of Warmia and Mazury in Olsztyn, statutory research grant no. 28.610.002–300 and from funds of the Faculty of Civil and Environmental Engineering, Gdansk University of Technology allocated for scientific research.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Meo, M.; Zumpano, G. On the optimal sensor placement techniques for a bridge structure. Eng. Struct. 2005, 27, 1488–1497. [Google Scholar] [CrossRef]
  2. Glišić, B.; Inaudi, D. Fibre Optic Methods for Structural Health Monitoring; John Wiley & Sons: Chichester, UK, 2007; p. 276. [Google Scholar]
  3. Cai, J.; Wang, J.; Wu, J.; Hu, C.; Grafarend, E.; Chen, J. Horizontal deformation rate analysis based on multiepoch GPS measurements in Shanghai. J. Surv. Eng. 2008, 134, 132–137. [Google Scholar] [CrossRef]
  4. Li, L.; Kuhlmann, H. Deformation detection in the GPS real–time series by the multiple Kalman filters model. J. Surv. Eng. 2010, 136, 157–164. [Google Scholar] [CrossRef]
  5. Yi, T.H.; Li, H.N.; Gu, M. Experimental assessment of high–rate GPS receivers for deformation monitoring of bridge. Measurement 2013, 46, 420–432. [Google Scholar] [CrossRef]
  6. Xiao, Z.; Liang, J.; Yu, D.; Asundi, A. Large field–of–view deformation measurement for transmission tower based on close–range photogrammetry. Measurement 2011, 44, 1705–1712. [Google Scholar] [CrossRef]
  7. Lindenbergh, R.; Pietrzyk, P. Change detection and deformation analysis using static and mobile laser scanning. Appl. Geomat. 2015, 7, 65–74. [Google Scholar] [CrossRef]
  8. Pingue, F.; Petrazzuoli, S.M.; Obrizzo, F.; Tammaro, U.; De Martino, P.; Zuccaro, G. Monitoring system of buildings with high vulnerability in presence of slow ground deformations (The Campi Flegrei, Italy, case). Measurement 2011, 44, 1628–1644. [Google Scholar] [CrossRef]
  9. Baselga, S.; García–Asenjo, L.; Garrigues, P. Deformation monitoring and the maximum number of stable points method. Measurement 2015, 70, 27–35. [Google Scholar] [CrossRef]
  10. Chen, Y.Q.; Chrzanowski, A.; Secord, J.M. A strategy for the analysis of the stability of reference points in deformation surveys. CISM JACSGG 1990, 44, 141–149. [Google Scholar] [CrossRef]
  11. Amiri–Simkooei, A.R.; Alaei–Tabatabaei, S.M.; Zangeneh–Nejad, F.; Voosoghi, B. Stability Analysis of Deformation–Monitoring Network Points Using Simultaneous Observation Adjustment of Two Epochs. J. Surv. Eng. 2017, 143, 1–12. [Google Scholar] [CrossRef] [Green Version]
  12. Koch, K.R. Parameter Estimation and Hypothesis Testing in Linear Models, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 1999; p. 334. [Google Scholar] [CrossRef]
  13. Durdag, U.M.; Hekimoglu, S.; Erdogan, B. Reliability of Models in Kinematic Deformation Analysis. J. Surv. Eng. 2018, 144, 1–15. [Google Scholar] [CrossRef]
  14. Duchnowski, R.; Wiśniewski, Z. Estimation of the shift between parameters of functional models of geodetic observations by applying Msplit estimation. J. Surv. Eng. 2012, 138, 1–8. [Google Scholar] [CrossRef]
  15. Wiśniewski, Z.; Zienkiewicz, M.H. Shift–M*split estimation in deformation analyses. J. Surv. Eng. 2016, 142, 1–13. [Google Scholar] [CrossRef]
  16. Chen, Y.Q. Analysis of Deformation Surveys—A Generalized Method; Technical Report No. 94; University of New Brunswick: Fredericton, NB, USA, 1983; pp. 54–72. [Google Scholar]
  17. Caspary, W.F. Concepts of Network and Deformation Analysis; The University of New South Wales: Utrecht, The Netherlands, 2000; p. 183. Available online: https://trove.nla.gov.au/version/39909771 (accessed on 27 February 2020).
  18. Even-Tzur, G. Deformation analysis by means of extended free network adjustment constraints. J. Surv. Eng. 2011, 137, 47–52. [Google Scholar] [CrossRef]
  19. Nowel, K.; Kamiński, W. Robust estimation from observation differences for free control networks. J. Geod. 2014, 88, 749–764. [Google Scholar] [CrossRef] [Green Version]
  20. Denli, H.H. Stable point research on deformation networks. Surv. Rev. 2008, 40, 74–82. [Google Scholar] [CrossRef]
  21. Mittermayer, E. Zur Ausgleichung freier Netze. Zeitschrift für Vermessungswesen (ZFV) 1972, 11, 481–489. [Google Scholar]
  22. Mittermayer, E. A generalization of the least–squares method for the adjustment of free networks. Bull. Géodésique 1972, 104, 139–157. [Google Scholar] [CrossRef]
  23. Wolf, H. Helmers Lösung zum Problem der freien Netze mit singulärer Normal–gleichungsmatrix. Zeitschrift für Vermessungswesen (ZFV) 1972, 5, 189–192. [Google Scholar]
  24. Perelmuter, A. Adjustment of free networks. Bull. Géodésique 1979, 53, 291–296. [Google Scholar] [CrossRef]
  25. Papo, H.B.; Perelmuter, A. Datum definition by free net adjustment. Bull. Géodésique 1981, 55, 218–226. [Google Scholar] [CrossRef]
  26. Caetano, E.; Silva, S.; Bateira, J. A vision system for vibration monitoring of civil engineering structures. Exp. Tech. 2011, 35, 74–82. [Google Scholar] [CrossRef]
  27. Flesch, R.; Kahmen, H. Building structures as kinematic systems—Dynamic monitoring and system analysis. Int. Assoc. Geod. Symp. 1999, 121, 313–318. [Google Scholar] [CrossRef]
  28. Reynolds, P. Dynamic testing and monitoring of civil engineering structures. Exp. Tech. 2008, 32. [Google Scholar] [CrossRef]
  29. Rogozhin, E.A.; Kapustian, N.K.; Antonovskaya, G.N. A New System for Monitoring Hydraulic Engineering Structures. Seism. Instrum. 2018, 54, 531–535. [Google Scholar] [CrossRef]
  30. Amezquita–Sanchez, J.P.; Adeli, H. Signal Processing Techniques for Vibration–Based Health Monitoring of Smart Structures. Arch. Comput. Methods Eng. 2016, 23, 1–15. [Google Scholar] [CrossRef]
  31. Krarup, T. A Contribution to the Mathematical Foundation of Physical Geodesy; Publish no. 44; Danish Geodetic Institute: Copenhagen, Denmark, 1969; p. 80. [Google Scholar]
  32. Moritz, H. A generalized least–squares model. Studia Geophysica Geodaetica 1970, 14, 353–362. [Google Scholar] [CrossRef]
  33. Moritz, H. Advanced least–squares methods. Department of Geodetic Science, Report; No. 175; Ohio State University: Columbus, OH, USA, 1972. [Google Scholar]
  34. Moritz, H. Least–Squares Collocation; Deutsche Geodatische Kommission Reihe a Heft: München, Germany, 1973; p. 75. [Google Scholar]
  35. Koch, K.R. Least squares adjustment and collocation. Bull. Géodésique 1977, 51, 127–135. [Google Scholar] [CrossRef]
  36. Tscherning, C.C. Collocation and least squares methods as a tool for handling gravity field dependent data obtained through space research techniques. Bull. Géodésique 1978, 52, 199–212. [Google Scholar] [CrossRef]
  37. Kotsakis, C. The multiresolution character of collocation. J. Geod. 2000, 74, 275–290. [Google Scholar] [CrossRef]
  38. Kotsakis, C. Least–squares collocation with covariance–matching constraints. J. Geod. 2007, 81, 661–677. [Google Scholar] [CrossRef]
  39. Deakin, R.E.; Collier, P.A.; Leahy, F.J. Transformation of coordinates using least squares collocation. Aust. Surv. 1994, 396–418. [Google Scholar] [CrossRef]
  40. Kutoglu, H.; Ayan, T.; Mekik, C. Integrating GPS with national networks by collocation method. Appl. Math. Comput. 2006, 177, 508–514. [Google Scholar] [CrossRef]
  41. You, R.J.; Hwang, H.W. Coordinate transformation between two geodetic datums of Taiwan by least–squares collocation. J. Surv. Eng. 2006, 132, 64–70. [Google Scholar] [CrossRef]
  42. Yang, Y.; Zeng, A.; Zhang, J. Adaptive collocation with application in height system transformation. J. Geod. 2009, 83, 403–410. [Google Scholar] [CrossRef]
  43. Mao, O.; Zhang, L.; Hu, Q.; Yu, J.; Feng, S.; Washington, O.; Gong, H. A least squares collocation method for accuracy improvement of mobile LiDAR systems. Remote Sens. 2015, 7, 7402–7424. [Google Scholar] [CrossRef] [Green Version]
  44. Tscherning, C.C. Testing frame transformation, gridding and filtering of GOCE gradiometer data by least–squares collocation using simulated data. IAG Proc. 2004, 128, 277–282. [Google Scholar] [CrossRef]
  45. Tscherning, C.C. The use of least–squares collocation for the processing of GOCE data. Oesterreichische Z. Fuer Vermess. Geoinf. (Austrian Mag. Surv. Geoinf.) 2010, 98, 21–26. [Google Scholar]
  46. Gómez, D.; Piñón, D.; Smalley, R.; Bevis, M.; Cimbaro, S.; Lenzano, L.; Barón, J. Reference frame access under the effects of great earthquakes: A least squares collocation approach for non–secular post–seismic evolution. J. Geod. 2016, 90, 263–273. [Google Scholar] [CrossRef]
  47. Yang, Y. Robustifying collocation. Manuscr. Geod. 1992, 17, 21–28. [Google Scholar]
  48. Schaffrin, B.; Wieser, A. On weighted total least–squares adjustment for linear regression. J. Geod. 2008, 82, 415–421. [Google Scholar] [CrossRef]
  49. Shen, Y.; Li, B.; Chen, Y. An iterative solution of weighted total least–squares adjustment. J. Geod. 2011, 85, 229–238. [Google Scholar] [CrossRef]
  50. Mahboub, V. On weighted total least–squares for geodetic transformations. J. Geod. 2012, 86, 359–367. [Google Scholar] [CrossRef]
  51. Ditlevsen, O. Traffic loads on large bridges modeled as white–noise fields. ASCE J. Eng. Mech. 1994, 120, 681–694. [Google Scholar] [CrossRef]
  52. Lifshits, M.A. Gaussian Random Functions; Spinger Science+Business Media: Dordrecht, The Netherlands, 1995; p. 337. [Google Scholar] [CrossRef]
  53. Koch, K.R. Bayesian Inference with Geodetic Applications; Springer: Berlin/Heidelberg, Germany, 1990; p. 199. [Google Scholar] [CrossRef]
  54. Osborne, M. Bayesian Gaussian Processes for Sequential Prediction, Optimization and Quadrature; New College, University of Oxford: London, UK, 2010. [Google Scholar]
  55. Magnus, J.R. Linear Structures; Charles Griffin and Company Ltd.: London, UK; Oxford University Press: London, UK, 1988. [Google Scholar]
  56. Bock, Y. Estimating crustal deformations from a combination of baseline measurements and geophysical models. J. Geod. 1983, 57, 294–311. [Google Scholar] [CrossRef]
  57. Rao, C.R. Linear Statistical Inference and Its Applications; John Wiley & Sons Inc.: New York, NY, USA, 1973. [Google Scholar] [CrossRef]
  58. Moritz, H. Statistical foundations of collocation. Boll. Geod. Sci. Affini 1980, 2, 131–150. [Google Scholar]
Figure 1. Simulated free control network. CP, controlled point; ECP, extended CP.
Figure 1. Simulated free control network. CP, controlled point; ECP, extended CP.
Sensors 20 03913 g001
Figure 2. Estimators u ^ k i and u ^ p , k i , k = 1 , , 25 of deterministic displacements of CPs, obtained in the sequence i = 1 , , 5000 of simulations (for u ¯ 1 = 10   mm , σ v = 0.3   mm , σ s = 1.0   mm , and σ E = 0.1   mm ).
Figure 2. Estimators u ^ k i and u ^ p , k i , k = 1 , , 25 of deterministic displacements of CPs, obtained in the sequence i = 1 , , 5000 of simulations (for u ¯ 1 = 10   mm , σ v = 0.3   mm , σ s = 1.0   mm , and σ E = 0.1   mm ).
Sensors 20 03913 g002
Figure 3. Random, deterministic, and total displacements of points within the simulated free control network (in mm).
Figure 3. Random, deterministic, and total displacements of points within the simulated free control network (in mm).
Sensors 20 03913 g003
Figure 4. The histograms of least squares (LS) and total least-squares collocation (TLSC) estimators regarding zero CPs displacements.
Figure 4. The histograms of least squares (LS) and total least-squares collocation (TLSC) estimators regarding zero CPs displacements.
Sensors 20 03913 g004
Figure 5. Graphic illustration of estimation results: (a) simulated and determined random displacements of CPs and ECPs, (b) simulated and determined surface noise in these points.
Figure 5. Graphic illustration of estimation results: (a) simulated and determined random displacements of CPs and ECPs, (b) simulated and determined surface noise in these points.
Sensors 20 03913 g005
Figure 6. The course of controlled values of norms in individual iteration steps.
Figure 6. The course of controlled values of norms in individual iteration steps.
Sensors 20 03913 g006
Figure 7. Monitored loading quay: (a) object’s photo (source: Google Maps), (b) free control network (FCN) diagram with ECPs marked.
Figure 7. Monitored loading quay: (a) object’s photo (source: Google Maps), (b) free control network (FCN) diagram with ECPs marked.
Sensors 20 03913 g007
Figure 8. Graphic illustration of (a) random displacement and (b) surface noise.
Figure 8. Graphic illustration of (a) random displacement and (b) surface noise.
Sensors 20 03913 g008
Figure 9. The courses of norms values in the iterative process: (a) random displacement (b) surface noise.
Figure 9. The courses of norms values in the iterative process: (a) random displacement (b) surface noise.
Sensors 20 03913 g009
Table 1. Comparative parameters J u ^ M C and J u ^ p M C obtained with the Monte Carlo (MC) method.
Table 1. Comparative parameters J u ^ M C and J u ^ p M C obtained with the Monte Carlo (MC) method.
σ s = 0
u ¯ 1 = 0 u ¯ 1 = 5 u ¯ 1 = 10
σ E σ E σ E
0.00.10.20.30.00.10.20.30.00.10.20.3
J u ^ M C 0.240.240.240.240.290.290.290.290.400.400.400.40
J u ^ p M C 0.240.240.240.240.290.290.290.290.400.400.400.40
σ s = 0.5
u ¯ 1 = 0 u ¯ 1 = 5 u ¯ 1 = 10
σ E σ E σ E
0.00.10.20.30.00.10.20.30.00.10.20.3
J u ^ M C 0.460.630.991.410.490.661.001.400.570.711.031.44
J u ^ p M C 0.240.240.240.240.290.290.290.290.400.400.400.40
σ s = 1.0
u ¯ 1 = 0 u ¯ 1 = 5 u ¯ 1 = 10
σ E σ E σ E
0.00.10.20.30.00.10.20.30.00.10.20.3
J u ^ M C 0.810.931.211.550.840.941.201.560.891.001.241.58
J u ^ p M C 0.240.240.240.240.290.290.290.290.400.400.400.40
Table 2. The results of estimation and prediction (in mm) for the simulated free control network (FCN).
Table 2. The results of estimation and prediction (in mm) for the simulated free control network (FCN).
CP u ¯ u ^ L S u ^ η s i m η ^ u ¯ η = u ¯ + η s i m u ^ η = u ^ + η ^
P 1 −5−4.43−5.050.870.41−4.13−4.64
P 2 00.860.050.570.600.570.64
P 3 01.090.270.090.610.090.87
P 4 01.080.510.110.360.110.86
P 5 02.502.420.78−0.130.782.29
P 6 01.441.840.45−0.610.451.23
P 7 00.290.63−0.07−0.56−0.070.08
P 8 01.401.500.11−0.310.111.19
P 9 01.721.360.710.140.711.51
P 10 00.880.140.810.520.810.66
P 11 0−0.73−1.77−0.050.82−0.05−0.95
P 12 0−0.11−0.21−0.26−0.12−0.26−0.33
P 13 00.331.10−0.78−0.99−0.780.11
P 14 0−0.400.57−0.61−1.19−0.61−0.62
P 15 00.691.290.14−0.820.140.48
P 16 00.440.67−0.11−0.45−0.110.22
P 17 0−1.48−0.64−1.15−1.05−1.15−1.69
P 18 0−1.47−0.65−1.58−1.04−1.58−1.69
P 19 0−1.13−1.07−1.04−0.28−1.04−1.35
P 20 0−1.55−0.47−0.730.70−0.73−1.77
P 21 00.270.18−0.86−0.14−0.860.05
P 22 0−0.50−0.02−0.98−0.69−0.98−0.72
P 23 0−1.88–1.19–1.38–0.91–1.38–2.09
P 24 0–0.85– 0.55–0.87–0.52–0.87–1.07
P 25 01.531.090.220.230.221.32
J u ^ L S = 0.92 J u ^ = 0.79 J η ^ = 0.62 J u ^ η = 0.55
ECP u ^ η s i m η ^ u ^ η = u ^ + η ^
P 1 0.160.190.610.77
P 2 0.390.100.510.90
P 3 1.490.410.211.70
P 4 −0.86−1.16−0.75−1.61
P 5 0.29−0.19−0.080.21
Table 3. Observation of heights’ differences within a real FCN (in mm).
Table 3. Observation of heights’ differences within a real FCN (in mm).
h ( 1 , 2 ) h ( 2 , 3 ) h ( 3 , 4 ) h ( 4 , 5 ) h ( 5 , 6 ) h ( 6 , 7 ) h ( 7 , 8 ) h ( 8 , 9 ) h ( 9 , 10 ) h ( 10 , 11 ) h ( 11 , 12 ) h ( 12 , 13 ) h ( 13 , 14 ) h ( 14 , 1 )
t α −6.911.6−13.2−25.3−3.6−14.920.6−3.48.5−15.78.20.112.1−29.7
t β −3.58.9−6.623.4−1.2−13.923.2−2.710.8−17.57.9−0.414.8−38.2
δ y i 3.4−2.76.6−1.92.41.02.60.72.3−1.8−0.3−0.52.7−8.5
Table 4. The results of estimation and prediction for the real FCN (in mm).
Table 4. The results of estimation and prediction for the real FCN (in mm).
CP u ^ L S u ^ η ^ u ^ η = u ^ + η ^
σ s = 0.2 σ s = 0.6 σ s = 0.2 σ s = 0.6 σ s = 0.2 σ s = 0.6
P 1 −5.97−6.10−6.350.030.10−6.07−6.24
P 2 −2.98−2.91−2.76−0.16−0.49−3.07−3.25
P 3 −6.12−5.60−4.58−0.61−1.82−6.21−6.39
P 4 0.050.621.76−0.66−1.99−0.04−0.23
P 5 −2.27−2.44−2.780.080.24−2.36−2.54
P 6 −0.29−0.72−1.570.331.00−0.39−0.57
P 7 0.270.210.09−0.03−0.090.18−0.00
P 8 2.452.332.090.030.082.362.17
P 9 2.722.702.67−0.08−0.232.622.44
P 10 4.604.163.300.341.034.504.32
P 11 2.352.252.060.010.022.262.07
P 12 1.591.721.96−0.21−0.641.511.32
P 13 0.660.801.09−0.24−0.710.560.38
P 14 2.952.973.02−0.12−0.352.852.67
ECP u ^ η ^ u ^ η = u ^ + η ^
P 1 −0.81−0.37−0.28−0.83−1.09−1.20
P 2 −0.13−0.510.070.21−0.06−0.30
P 3 2.482.320.010.042.492.36
P 4 3.522.870.270.813.793.68

Share and Cite

MDPI and ACS Style

Wiśniewski, Z.; Kamiński, W. Estimation and Prediction of Vertical Deformations of Random Surfaces, Applying the Total Least Squares Collocation Method. Sensors 2020, 20, 3913. https://doi.org/10.3390/s20143913

AMA Style

Wiśniewski Z, Kamiński W. Estimation and Prediction of Vertical Deformations of Random Surfaces, Applying the Total Least Squares Collocation Method. Sensors. 2020; 20(14):3913. https://doi.org/10.3390/s20143913

Chicago/Turabian Style

Wiśniewski, Zbigniew, and Waldemar Kamiński. 2020. "Estimation and Prediction of Vertical Deformations of Random Surfaces, Applying the Total Least Squares Collocation Method" Sensors 20, no. 14: 3913. https://doi.org/10.3390/s20143913

APA Style

Wiśniewski, Z., & Kamiński, W. (2020). Estimation and Prediction of Vertical Deformations of Random Surfaces, Applying the Total Least Squares Collocation Method. Sensors, 20(14), 3913. https://doi.org/10.3390/s20143913

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop