Home > Products > GeoLab > Articles > Using Observation Groups in GeoLab
Welcome Products Services Contact Us Links
 
Products
GeoLab
Explorer
Closure
Smoother
Surface Area Calculator
Color Picker
File Format Tools
 
Services
Adjustments
Training
Software
Geoids

 

Using Observation Groups in GeoLab  

This article describes the various options and formatting requirements for groups of observations in GeoLab. The different types of observation groups covered are:

Important: You should read the Getting Started with GeoLab article before reading this one.

General Layout of Observation Groups
All observation groups in GeoLab should begin with a GRP record, which effectively names the group of measurements. The group name has no effect on the network adjustment, but it simplifies the revisions to your IOB file due to flagged residuals in the output listing. The observations named by a GRP record are all of those measurements following that GRP record, up to the next GRP record encountered.

The general structure of observation groups is as follows (please note that we will use square brackets, [], to enclose a description of a group of one or more records):

 GRP  Observation Group 001
[Observations in Group 001]
 GRP  Observation Group 002
[Observations in Group 002]
etc...

The general layout of records for the 2DC, 2DD, 3DC, and 3DD observation groups is as follows:

[Header record (2DC, 2DD, 3DC, or 3DD)]
[Set of measurement value records]
[Covariance matrix header record]
[Matrix elements records]

The "Covariance matrix header record" and the "Matrix elements records" are formatted the same for all types of header records (2DC, 2DD, 3DC, and 3DD), so we will cover their formatting here.

Covariance Matrix Formatting
The specific formatting of records for the "covariance matrix" of 2DC, 2DD, 3DC, and 3DD groups of observations depends on the form of covariance information you have for the observation group. The most common type of covariance matrix data provided for coordinate and coordinate difference measurements are the actual elements of a covariance matrix (in which case you must use a COV matrix header record). Other forms handled by GeoLab are weight matrix elements (WGT header record), correlation matrix elements (CORR header record), or you can generate a covariance matrix (e.g. for simulations) using the GENC header record (which requires no ELEM records).

The required size of the matrix (which determines the number of ELEM records in the "Matrix elements records" section) is determined by the number N of records in the "Set of measurement value records" as follows:

Group type Obs. record type(s) Matrix size
2DC PL, NE, XY 2 * N
2DD PL, NE, XY 2 * (N - 1)
3DC PLH, PLO, NEH, NEO, XYZ 3 * N
3DD PLH, PLO, NEH, NEO, XYZ 3 * (N - 1)
3DD DXYZ 3 * N

For example, if we have the following three (N = 3) sets of coordinate difference observations:

Station 1 to station 2: d1 = X2 - X1; d2 = Y2 - Y1; d3 = Z2 - Z1;
Station 2 to station 3: d4 = X3 - X2; d5 = Y3 - Y2; d6 = Z3 - Z2;
Station 3 to station 1: d7 = X3 - X1; d8 = Y3 - Y1; d9 = Z3 - Z1;

for a total of nine observations (3 sets of dX, dY, and dZ), we would input them in three DXYZ records as follows:

 3DD
 DXYZ12 d1 d2 d3
 DXYZ23 d4 d5 d6
 DXYZ31 d7 d8 d9

If we denote a covariance matrix element as sij (where s11 is the variance of d1, and s12 is the covariance between d1 and d2, etc) the upper-triangular portion of the covariance matrix for these observations is as follows (note that the order of the covariance elements depends only on the order in which the observations are given as DXYZ records):

s11 s12 s13 s14 s15 s16 s17 s18 s19
  s22 s23 s24 s25 s26 s27 s28 s29
    s33 s34 s35 s36 s37 s38 s39
      s44 s45 s46 s47 s48 s49
        s55 s56 s57 s58 s59
          s66 s67 s68 s69
            s77 s78 s79
              s88 s89
                s99

The COV and ELEM records for this observation would therefore be:

 COV  CT UPPR
 ELEM s11 s12 s13
 ELEM s14 s15 s16
 ELEM s17 s18 s19
 ELEM s22 s23 s24
 ELEM s25 s26 s27
 ELEM s28 s29
 ELEM s33 s34 s35
 ELEM s36 s37 s38
 ELEM s39
 ELEM s44 s45 s46
 ELEM s47 s48 s49
 ELEM s55 s56 s57
 ELEM s58 s59
 ELEM s66 s67 s68
 ELEM s69
 ELEM s77 s78 s79
 ELEM s88 s89
 ELEM s99

As you can see, the general approach is to start the entry of each row of the upper-triangular portion of the matrix with a new ELEM record, and you use as many ELEM records as required for the elements of each row (an ELEM record can hold up to three elements). For the UPPR matrix form, each row is entered starting with the diagonal element. For the DIAG matrix form, you simply enter the diagonal elements (three at a time) in ELEM records from the upper-left element to the lower-right element.

The structure of the ELEM records is similar when you use the WGT matrix header (except, of course, weight matrix elements are entered instead of covariance matrix elements).

The structure of the ELEM records is slightly different when you use the CORR matrix header record. In this case the diagonal elements must always be 1.0, and the off-diagonal elements are the correlation coefficients. Immediately following these elements, you must also provide ELEM records that specify the standard deviations of the coordinate or coordinate difference measurements. If the matrix size is N, then N standard deviations must be entered in as many ELEM records as required (an ELEM record can only contain up to three elements).

 
2DC: 2D Coordinate Observation Groups
The format for a 2D coordinate observation group is as follows (the use of the optional sigma-ID is described below):
 
 2DC  [optional sigma ID]
[Set of 2D coordinate records]
[Covariance matrix header record]
[Matrix elements records]

The "Set of 2D coordinate records" must contain only PL, NE, or XY coordinate records. It may also contain HI, HT, 4PAR, and 7PAR records.
 

2DD: 2D Coordinate Difference Observation Groups
The format for a 2D coordinate difference observation group is as follows:
 
 2DD  [optional sigma ID]
[Set of 2D coordinate records]
[Covariance matrix header record]
[Matrix elements records]

The "Set of 2D coordinate records" must contain only PL, NE, or XY coordinate records (it may also contain HI, HT, 4PAR, and 7PAR records). GeoLab will automatically subtract the first of these records from the others to arrive at the actual coordinate difference measurement values. Note that the covariance matrix elements entered must correspond to these computed coordinate differences.
 

3DC: 3D Coordinate Observation Groups
The format for a 3D coordinate observation group is as follows:
 
 3DC  [optional sigma ID]
[Set of 3D coordinate records]
[Covariance matrix header record]
[Matrix elements records]

The "Set of 3D coordinate records" must contain only PLH, PLO, NEH, NEO, or XYZ coordinate records (it may also contain HI, HT, 4PAR, and 7PAR records).
 

3DD: 3D Coordinate Difference Observation Groups
The format for a 3D coordinate observation group is as follows:
 
 3DD  [optional sigma ID]
[Set of 3D coordinate records or DXYZ records]
[Covariance matrix header record]
[Matrix elements records]

The "Set of 3D coordinate records or DXYZ records" must contain only PLH, PLO, NEH, NEO, XYZ coordinate records, or DXYZ coordinate difference records (it may also contain HI, HT, 4PAR, and 7PAR records). If DXYZ records are not used, GeoLab will automatically subtract the first of these records from the others to arrive at the actual coordinate difference measurement values (in which case the covariance matrix elements entered must correspond to these computed coordinate differences).
 

Using the Covariance Matrix Header Record and Sigma Records
Each of the coordinate and coordinate-difference observation group header records (2DC, 2DD, 3DC, and 3DD) takes an optional "sigma ID". As you know, a "sigma-ID" refers to a previously given SIGM record. The elements of the sigma-record are used in these groups of measurements as follows:
 
SIGM Element Description
Standard deviation This square of this value is added to the "Addition constant (diagonal)" value specified in the matrix header record.
PPM This value is added to the "PPM (diagonal)" value specified in the matrix header record.
Centering errors For 2DC and 3DC groups, the square of the at-centering value is added to the "Factor (diagonal)" value specified in the matrix header record.  For 2DD and 3DD groups, the sum of the squares of the from-centering and to-centering values is added to the "Factor (diagonal)" value specified in the matrix header record.
Auxiliary parameter Adds the specified auxiliary parameter to the observation group.

Note that if centering errors are used in the sigma-record, the covariance matrix coordinate system must be LG (local geodetic).
 
The covariance matrix header record can be one of the following record types:

Record Type Description
COV Specifies that covariance matrix elements follow in a set of ELEM records.
CORR Specifies that correlation matrix elements and standard deviation values follow in a set of ELEM records.
GENC Specifies that the covariance matrix will be built from the fields of the GENC record (see below) and that no matrix elements are specified (no ELEM records follow).
WGT Specifies that weight matrix elements follow in a set of ELEM records.

The "matrix revision" fields of the matrix header record (COV, CORR, GENC, and WGT) are as follows:

Columns Description
015-024 Addition constant (entire matrix)
026-035 Factor (entire matrix)
037-046 Addition constant (diagonal)
048-057 Factor (diagonal)
059-068 PPM (diagonal)
070-079 Addition constant (height/Z)
081-090 Factor (height/Z)

For the GENC matrix header, all elements of the covariance matrix are initially set to zero before the following revisions are made. For all other matrix header types, the elements of the covariance matrix are initially given (either directly or indirectly) by the ELEM records.
 
The "Addition constant (entire matrix)" and the "Addition constant (diagonal)" are first added directly to the covariance matrix. The "PPM (diagonal)" value is used to calculate a standard deviation (using either a geocentric radius for 2DC and 3DC groups, or the vector lengths for 2DD and 3DD groups) and that value is squared and added to the diagonal of the covariance matrix.

Next, for 3DC and 3DD groups, the "Addition constant (height/Z)" and "Factor (height/Z)" values are used to revise the covariance matrix. If the matrix coordinate system is CT then the Z or Z-difference diagonal elements are revised. If the matrix coordinate system is LG then the height or height-difference diagonal elements are revised. The revision is performed as follows (E denotes the current matrix element, A the  "Addition constant (entire matrix)" value, and F the "Factor (height/Z)" value:

   E = (E + A) * F.

Finally, the "Factor (entire matrix)" and the "Factor (diagonal)" are applied to the covariance matrix.
 

 
DSET: Direction Observation Groups
A (horizontal) direction set must be in the following format:
 
 DSET [optional sigma ID]
[Set of DIR observation records]

Each DIR record in a DSET must have the same from-station specified.

If a sigma-record is used for a DSET, the standard deviation, the from-centering, and to-centering values are used to modify or define the standard deviations of the direction measurements in the DSET group.
 

Please let us know if you would like to see this article expanded in any way. Thanks.
 

 
 
©2006 BitWise Ideas Inc. All Rights Reserved.
All names are (registered) trademarks of their respective companies.