Mathematical Problems in Engineering

Volume 2019, Article ID 5173956, 10 pages

https://doi.org/10.1155/2019/5173956

## Robust Prediction Algorithm Based on a General EIV Model for Multiframe Transformation

^{1}NASG Key Laboratory of Land Environment and Disaster Monitoring, China University of Mining and Technology, Xuzhou 221116, Jiangsu, China^{2}School of Environment Science and Spatial Informatics, China University of Mining & Technology, Xuzhou 221116, Jiangsu, China^{3}College of Civil Engineering, Anhui Jianzhu University, Hefei 230000, Anhui, China^{4}School of Geodesy and Geomatics, Anhui University of Science and Technology, Huainan 232001, Anhui, China

Correspondence should be addressed to Jingxiang Gao; nc.ude.tmuc@oagxj

Received 15 November 2018; Revised 6 January 2019; Accepted 30 January 2019; Published 11 February 2019

Academic Editor: Jixiang Yang

Copyright © 2019 Lihui Yao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

In modern geodesy, there are cases in which the target frame is unique and there is more than one source frame. Helmert transformations, which are extensively used to solve transformation parameters, can be separately solved between the target frame and one of the source frames. However, this is not globally optimal, even though each transformation is locally optimal on its own. Additionally, this also generates the problem of multiple solutions in the noncommon station of the target frame. Moreover, least squares solutions can cause estimation value distortion, with a gross error existing in observations. Thus, in this paper, Helmert transformations among three frames, that is, one target frame and two source frames, are studied as an example. A robust prediction algorithm based on the general errors-in-variables prediction algorithm and the robust estimation is derived in detail and is applied to achieve multiframe total transformation. Furthermore, simulation experiments were conducted and the results validated the superiority of the proposed total transformation method over classical separate approaches.

#### 1. Introduction

The precision and reliability of terrestrial reference frames have improved with the rapid development of the Global Navigation Satellite System (GNSS) [1]. At present, common terrestrial reference frames include WGS-84, PE-90, CGCS-2000, and ITRF2008. The three-dimensional (3D) space coordinates of ground control points in the source coordinate system are obtained using GNSS observations, whereas local coordinate systems are often used for geodetic surveying. Therefore, it is necessary to use coordinate transformation to unify the data. Similarity datum transformation is widely used in 3D coordinate transformation. Common point coordinates are typically used to solve the transformation parameters, and then the coordinates of noncommon points in the old coordinate system are converted to the new coordinate system according to the transformation parameter solution.

Seven parameters of 3D datum transformations can be used to relate two reference frames: three for translation, one for scale, and three for rotation. Using two sets of coordinates measured in two frames, these parameters can be solved using the Gauss–Markov (GM) model of least squares (LS) [2]. The assumption of the GM model is that there are no random errors in the coefficient matrix, but the observed coordinates of the reference station in the source coordinate system of the coefficient matrix have random errors; thus, the GM model cannot be used to describe the coordinate transformation problem accurately. For the case in which the coefficient matrix contains random errors, the error-in-variables (EIV) model can be used to describe the GM model, and the LS solution of the EIV model is called the total least squares (TLS).

In recent years, TLS has become a research topic of great interest in the field of geodesy and many algorithms have been developed and examined [3–6]. When the small rotation angle and scale are close to one, the datum transformation model can adopt the Bursa–Wolf model. Compared with classical LS, the accuracy of estimation results using TLS is greatly improved [7–10]. In the case of a large rotation angle and arbitrary scale, Felus et al. proposed an algorithm for a multiple EIV model [11]. Fang transformed the EIV model into a nonlinear Gauss–Markov model and proposed the 3D datum transformation algorithm under any rotation angle by taking the coordinates of common points in the source coordinate system as virtual observation values [12]. Chang and Li et al. proposed the analytical and closed-form solutions of 3D datum transformation, respectively [13, 14]. Because the rotation matrix of the datum transformation is a completely orthogonal matrix, there are strict orthogonal constraints. Fang proposed the universal constraint model of weighted TLS [5]. Lin et al. proposed the constraint solution of 3D datum transformation based on the nonlinear Gauss–Helmert model [15].

However, most TLS algorithms for datum transformation that focus on how to solve transformation parameters do not combine the two independent processes of transformation parameter estimation and noncommon station transformation; that is, noncommon station transformation is performed at the same time as the estimation of transformation parameters. In practice, there is a statistically significant correlation between the common and noncommon stations [16]; for example, the coordinates of common and noncommon stations come from the same control network adjustment results (e.g., GPS baseline network), but the existing model algorithm does not take into consideration the correlation between reference stations, so the estimation results cannot obtain the optimal solution in theory. To overcome this problem, Li et al. [17, 18] proposed the concept of seamless frame datum transformation, which combines the estimation parameters and coordinate transformation of noncommon stations. Compared with the traditional model, this transformation effectively improves the accuracy of the transformation coordinates. Wang et al. proposed the estimation model of 3D datum transformation under a large rotation angle and arbitrary scale conditions based on the generalized EIV model [19]. Lin et al. proposed the prediction model of 3D datum transformation of source coordinates in the form of geodetic coordinates based on the nonlinear Gauss–Helmert model [20].

Most of the current literature mainly concentrates on the datum transformation between two frames. When multiple (two or more) frames are involved, the transformation between the multiple source frames and a target frame, as mentioned earlier, clearly can be used to convert each source frame in turn with the target frame. The advantage is that we can draw lessons from the existing research achievements of Helmert transformation between the two frames directly. However, the existing methods cannot obtain the optimal solution. In recent years, a few researchers have studied the problem of simultaneous Helmert transformation with multiple frames. For example, Han et al. and Aktug proposed the closed solutions of this problem from the perspectives of geometry and algebra, respectively [21, 22], but the solutions all lack statistical optimality; in other words, the solutions are not LS.

The parameter estimation results of TLS are seriously affected when the observation of source frames contains gross errors. Mahboub et al. proposed a method of robust TLS based on selecting weight iteration for the EIV model [23]. Lu et al. and Tao et al. studied the robust TLS (RTLS) algorithm for 3D datum transformation and GPS height fitting, respectively, based on the Gauss–Helmert model [24, 25]. Jin et al. and Zhao et al. studied the total minimum L1 norm estimation method, taking into account the error of all variables [26, 27]. Wang et al. proposed the antidifference solution of weighted TLS (WTLS) by standardizing the residual to construct the equal-weight function and using the median method [6]. To solve the problem of processing outliers in universal 3D symmetrical Helmert transformation (probably with large rotation angles), Chang et al. derived the M-estimator by introducing the reweighted LS scheme [28]; in addition, the data snooping procedure was extended to the generalized EIV model by Wang et al. [29], consequently solving this problem by taking another approach.

Continuing from the analysis given above, we propose a robust predicted estimation for multiframe transformation, which constructs the weight factor function with standardized residuals and uses the median method to obtain the unit weight mean error. We need to consider the random errors of all observations and their correlation to finally achieve multiframe total transformation based on the generalized EIV model.

#### 2. Extrapolation Estimation Method Based on the Generalized EIV Model

The function model for the extrapolation estimation method based on the generalized EIV model [19, 20] is as follows:where the n × 1 vector denotes measurements and is the corresponding measurement error; the m × 1 vector denotes estimated values; the n × 1 vector and m × 1 vector are the measurements of the common and noncommon points, respectively, in the old frame; and are the corresponding measurement errors; is the t × 1 vector of estimated parameters; and is an abstract function.

The corresponding stochastic model can be described aswhere ,, and are the cofactor matrices of , , and , respectively; , , and are the intercofactor matrices of , , and , respectively.

Taylor series expansion is applied to the right-hand side of (1) and (2) at the approximate values , , and of , , and , respectively. Thus, (1) and (2) are now expressed aswhere denotes the correction vector of parameters, ,Then the standard LS form of (5) and (6) is described as follows: where , , , .

On the basis of the weighted LS error criterion, the LS solutions of the parameters and the estimation of the error vector can be calculated using (8):where .

#### 3. Robust Estimation Method Based on the Standardized Residuals and Median

The accurate estimated parameters are difficult to obtain using (8) if there are gross errors in the observation vector or coefficient matrix . For this, we propose a robust prediction algorithm based on the generalized EIV model. We use the robust estimation introduced by Wang et al. [6]; it is based on the standardized residuals and median robust model. The robust prediction model is more reasonable in theory than the existing TLS robust methods and has two advantages: the weight factor function constructed by the standardized residuals can take into account the robustness of both the observation space and structure space; the estimation of the unit weight mean error obtained by the median method has greater robustness in the iterative process. The estimation criterion of the robust prediction algorithm is described as follows:where where , , and are the equivalent cofactor matrices of observation vectors , , and , respectively; , , and are the equivalent intercofactor matrices of observation vectors , , and , respectively.

Then, the solutions of the parameters and the estimation of the error vector are as follows:

Because the correlation between the observations is taken into account in the stochastic model of the extrapolation estimation method, it is necessary to construct the equivalent cofactor matrices in the form of double factors as follows [6, 30]:where the subscripts and represent the serial numbers of the observed values; and denote the elements in cofactor matrix Q; and denote the elements in equivalent cofactor matrix ; and and are the reciprocals of the weight factors, which are called the cofactor factors [30]. Additionally, because the generalized EIV model is adopted in this paper, the prior information is directly the weight matrix/cofactor matrix of the observation vector, and there is no weight reduction problem of the standard EIV for the fixed elements of the coefficient matrix.

The corresponding cofactor factor function can be obtained by inverting the IGGIII weight factor function as follows.

It should be noted that the theoretical value of the corresponding cofactor should be infinite because the weight factor of the third line in (16) is zero. However, to facilitate calculation, a large number is used to satisfy the requirement in the numerical calculation, and the abnormal observation information has little effect on parameter estimation.

Because of the accuracy difference between variables, the cofactors are calculated with standardized residual to guarantee the antierror effect. The formula of the standardized residual is as follows:where is the th element of residual vector , is the cofactor of , and is the standard error of unit weight, which is expressed as follows:According to the law of cofactor propagation, we can easily obtain cofactor matrices of residual errors; see Wang et al. [6] for the detailed derivation process. Therefore, cofactor matrices can be estimated aswhere and .

It should be pointed out that if the number of redundancies is low, outliers would significantly affect the residuals of normal observations and the observation with larger residuals may not contain outliers. In other words, almost all methods that process outliers will break down when there are only a few redundant observations.

#### 4. Robust Prediction Algorithms for Multiframe Transformation

Without loss of generality, in this paper, we consider the transformation between two source frames and one target frame. The two source frames are called the first and second source frames. Assume the following four observation datasets: (i) the common stations of three frames *Π*_{1} = ; (ii) the common stations of the first and second source frames *Π*_{2} = ; (iii) the noncommon stations of the first source frame *Π*_{3} = ; and (iv) the noncommon stations of the second source frame *Π*_{4} = . The terms , , , and are used to represent the 3D coordinate vectors of the common stations in the above four observation datasets, respectively. Considering the first dataset as an example, , , and denote the 3D coordinate vectors of the target frame, first source frame, and second source frame, respectively. The observation or estimation equations with respect to the ith element in each dataset can be expressed aswhere denotes the error vector that corresponds to the coordinate observation vector; and are the 3 × 1 translation parameters between the first source frame or second source frame and the target frame, respectively; and are the scale parameters between the first source frame or second source frame and the target frame, respectively; and and are the 3 × 3 rotation matrices between the first source frame or second source frame and the target frame, respectively, and can be described as follows:where , , and are the rotation parameters between the first source frame and the target frame; , , and are the rotation parameters between the second source frame and the target frame.

The corresponding random model is as follows:

In a similar manner, (20) and (21) are expanded using the Taylor series expansion method. Then,wherewhere the format of and is the same as that of ; the format of and is the same as that of .

All observations can be linked to obtain where

Then the standard LS form of (29) is described as follows: where

Following (26) and (31), the robust LS solution of (31) can be obtained using the robust prediction estimation method in the previous section.

#### 5. Experimental Analysis

Assume that (i) the transformation parameters between the first source frame and the target frame are = , = 0.83, = 0.45, = 0.37, = 0.54; (ii) the transformation parameters between the second source frame and the target frame are = , = 0.79, = 0.23, = 0.48, = 0.33. The four dimensions of the corresponding observation datasets are = 10, = 6, = 5, and =4, respectively.

Considering the WGS-84 frame as the target frame, the geodetic coordinates of the target frame, including the latitude, longitude, and geodetic height, which are randomly generated in the form of a uniform distribution within the ranges , , and , are transformed into the corresponding 3D rectangular coordinates. The observation station coordinates of the first and second source frames are obtained from the true value of the given transformation parameters. The 3D distribution of the simulation data is shown in Figure 1.