710 likes | 720 Views
This symposium paper discusses the algebraic and kinematic approaches to the ITRF formulation problem, which involves combining time series of sub-network coordinates from various techniques into coordinates for the whole network.
E N D
The VII Hotine-Marussi Symposium Rome, July 6–10, 2009 Theoretical foundations of ITRF determination The algebraic and the kinematic approach Zuheir Altamimi1 & Athanasios Dermanis2 (1) IGN-LAREG - (2) Aristotle University of Thessaloniki
THE ITRF FORMULATION PROBLEM t given a time sequence of sub-network coordinates (one from each technique T = VLBI, SLR, GPS, DORIS) combine them into coordinates for the whole network obeying a time-evolution model Essentially: Determine the model parameters for each network point i
THE ITRF FORMULATION PROBLEM t given a time sequence of sub-network coordinates (one from each technique T = VLBI, SLR, GPS, DORIS) combine them into coordinates for the whole network obeying a time-evolution model Essentially: Determine the model parameters for each network point i
THE COORDINATE-FREE APPROACH (ALMOST) TO THE ITRF FORMULATION PROBLEM t Given a time sequence of sub-network shapes (one from each technique: VLBI, SLR, GPS, DORIS)
THE COORDINATE-FREE APPROACH (ALMOST) TO THE ITRF FORMULATION PROBLEM t Replace them with a smooth sequence of shapes
THE COORDINATE-FREE APPROACH (ALMOST) TO THE ITRF FORMULATION PROBLEM t Replace them with a smooth sequence of shapes
THE COORDINATE-FREE APPROACH (ALMOST) TO THE ITRF FORMULATION PROBLEM t Note that although shape variation is insignificant coordinates may vary significantly due to temporal instability in reference system maintenance
THE COORDINATE-FREE APPROACH (ALMOST) TO THE ITRF FORMULATION PROBLEM t To remove coordinate variation assign a different reference system at each epoch
THE COORDINATE-FREE APPROACH (ALMOST) TO THE ITRF FORMULATION PROBLEM t To remove coordinate variation assign a different reference system at each epoch
THE COORDINATE-FREE APPROACH (ALMOST) TO THE ITRF FORMULATION PROBLEM t To remove coordinate variation assign a different reference system at each epoch such that when networks are viewed in the “same” system
THE COORDINATE-FREE APPROACH (ALMOST) TO THE ITRF FORMULATION PROBLEM t To remove coordinate variation assign a different reference system at each epoch such that when networks are viewed in the “same” system coordinates vary in a smooth way
THE COORDINATE-FREE APPROACH (ALMOST) TO THE ITRF FORMULATION PROBLEM t To remove coordinate variation assign a different reference system at each epoch such that when networks are viewed in the “same” system coordinates vary in a smooth way in conformance with a coordinate time-variation model
THE COORDINATE-FREE APPROACH (ALMOST) TO THE ITRF FORMULATION PROBLEM t To remove coordinate variation assign a different reference system at each epoch such that when networks are viewed in the “same” system coordinates vary in a smooth way in conformance with a coordinate time-variation model currently:
data: STACKING FOR EACH PARTICULAR TECHNIQUE coordinate transformation parameters: t t coordinate variation model: model parameters: t
THE ITRF FORMULATION PROBLEM = SIMULTANEOUS STACKING FOR ALL TECHNIQUES VLBI SLR GPS DORIS ITRF t
THE ITRF FORMULATION PROBLEM IN AN OPERATIONALLY CONVENIENT COMPROMISE (1) Separate stackings one for each technique: Provides initial coordinates and velocities for the subnetwork of each technique Separation into 2 steps:
THE ITRF FORMULATION PROBLEM IN AN OPERATIONALLY CONVENIENT COMPROMISE (1) Separate stackings one for each technique: Provides initial coordinates and velocities for the subnetwork of each technique Separation into 2 steps: (2) Combination of initial coordinates and velocities: Provides initial coordinates and velocities for the whole ITRF network
THE MODEL FOR TIME EVOLUTION OF COORDINATES The (general) model: Point Pi coordinates: ai = point Pi parameters The current model:
WHAT THE MODEL DOES SMOOTHING Replaces observes time sequences of sub-network shapes with a single smooth time sequence for the whole ITRF network IMPOSES THE USE OF A REFERENCE SYSTEM so that network shapes are represented by coordinates INTERPOLATION Provides shapes expressed by coordinates for epochs other than observation ones
WHAT THE MODEL DOES NOT DO SMOOTHING Replaces observes time sequences of sub-network shapes with a single smooth time sequence for the whole ITRF network IMPOSES THE USE OF A REFERENCE SYSTEM so that network shapes are represented by coordinates It does not resolve the problem of the choice of the reference system INTERPOLATION Provides shapes expressed by coordinates for epochs other than observation ones MAIN ITRF FORMULATION PROBLEM Assign a reference system for each epoch
WHAT THE MODEL DOES NOT DO MAIN ITRF FORMULATION PROBLEM Assign a reference system for each epoch It does not resolve the problem of the choice of the reference system
WHAT THE MODEL DOES NOT DO MAIN ITRF FORMULATION PROBLEM Assign a reference system for each epoch It does not resolve the problem of the choice of the reference system PROBLEM SOLUTION:Introduce additional minimal constraintsin the Least-Square data analysis problem Minimal constraints: At any epoch t they determine the reference system without affecting the optimal network shape uniquely determined by the least-squares principle for the determination of ITRF parameters
WHAT THE MODEL DOES NOT DO MAIN ITRF FORMULATION PROBLEM Assign a reference system for each epoch It does not resolve the problem of the choice of the reference system PROBLEM SOLUTION:Introduce additional minimal constraintsin the Least-Square data analysis problem Minimal constraints: At any epoch t they determine the reference system without affecting the optimal network shape uniquely determined by the least-squares principle for the determination of ITRF parameters How to choose the minimal inner constraints? 2 approaches: (1) The algebraic approach (classical Meissl inner constraints) (2) The kinematic approach (new!)
THE ALGEBRAIC APPROACH THE KINEMATIC APPROACH Formulation of Least Squares problem with infinite solutions for different choices of reference system
THE ALGEBRAIC APPROACH THE KINEMATIC APPROACH Formulation of Least Squares problem with infinite solutions for different choices of reference system Choice of unique solution by inner constraints: or partial inner constraints:
THE ALGEBRAIC APPROACH THE KINEMATIC APPROACH Formulation of Least Squares problem Choice of reference system by minimization of apparent variation of coordinate for network points with infinite solutions for different choices of reference system Choice of unique solution by inner constraints: or partial inner constraints:
THE ALGEBRAIC APPROACH THE KINEMATIC APPROACH Formulation of Least Squares problem Choice of reference system by minimization of apparent variation of coordinate for network points with infinite solutions for different choices of reference system Measures of coordinate variation: Choice of unique solution by inner constraints: (1) Minimum relative kinetic energy = = vanishing relative angular momentum Discrete Tisserand Reference System (2) constant network barycenter or partial inner constraints: (3) constant mean quadratic scale
THE ALGEBRAIC APPROACH – INNER CONSTRAINTS Inner constraints determined from the linear variation of unknown parameters x when coordinate system changes with small transformation parameters p rotation angles translation vector scale parameter
THE ALGEBRAIC APPROACH – INNER CONSTRAINTS Inner constraints determined from the linear variation of unknown parameters x when coordinate system changes with small transformation parameters p rotation angles translation vector scale parameter Determine the parameter variation equations
THE ALGEBRAIC APPROACH – INNER CONSTRAINTS Inner constraints determined from the linear variation of unknown parameters x when coordinate system changes with small transformation parameters p rotation angles translation vector scale parameter Determine the parameter variation equations Then the (total) inner constraints are
MODEL PRESERVING APPROXIMATE TRANSFORMATIONS OF INITIAL COORDINATES AND VELOCITIES Transformation of coordinates in first order approximation
MODEL PRESERVING APPROXIMATE TRANSFORMATIONS OF INITIAL COORDINATES AND VELOCITIES Transformation of coordinates in first order approximation Model preserving transformations
MODEL PRESERVING APPROXIMATE TRANSFORMATIONS OF INITIAL COORDINATES AND VELOCITIES Transformation of coordinates in first order approximation Model preserving transformations Transformation of model parameters
THE ALGEBRAIC APPROACH – INNER CONSTRAINTS PER STATION For each station Pi determine the parameter variation equations The inner constraints per station are The (total) inner constraints are
MODEL PRESERVING APPROXIMATE TRANSFORMATIONS OF INITIAL COORDINATES AND VELOCITIES The use of model preserving transformations instead of arbitrary transformations leads to a sub-optimal solution: No matter what the optimality criterion, there exist an arbitrary transformation leading to a better solution which does not conform with the chosen model Strict optimality leads the solution OUTSIDE the adopted model !
MODEL PRESERVING APPROXIMATE TRANSFORMATIONS OF INITIAL COORDINATES AND VELOCITIES Transformation of model parameters in terms of corrections to approximate values Transformation of corrections to model parameters
MODEL PRESERVING APPROXIMATE TRANSFORMATIONS OF INITIAL COORDINATES AND VELOCITIES Transformation of model parameters in terms of corrections to approximate values Transformation of corrections to model parameters
MODEL PRESERVING APPROXIMATE TRANSFORMATIONS OF INITIAL COORDINATES AND VELOCITIES Transformation of model parameters in terms of corrections to approximate values Transformation of corrections to model parameters inner constraints sub-matrix
THE STACKING PROBLEM Original observation model GIVEN SOUGHT Observed coordinates in particular technique at epoch tk ITRF model coordinates at epoch tk NUISANCE Transformation parameters from ITRF system to technique-system at epoch tk
THE STACKING PROBLEM Original observation model GIVEN SOUGHT Observed coordinates in particular technique at epoch tk ITRF model coordinates at epoch tk NUISANCE Transformation parameters from ITRF system to technique-system at epoch tk In first order approximation In terms of corrections to approximate values
INNER CONSTRAINTS FOR THE STACKING PROBLEM Change of ITRF reference system
INNER CONSTRAINTS FOR THE STACKING PROBLEM (Total) inner constraints initial orientation initialtranslation initialscale orientation rate translation rate scale rate
INNER CONSTRAINTS FOR THE STACKING PROBLEM Partial inner constraints – Coordinates & velocities initialorientation initialtranslation initialscale orientation rate translation rate scale rate
INNER CONSTRAINTS FOR THE STACKING PROBLEM Partial inner constraints – Transformation parameters initialorientation initialtranslation initialscale orientation rate translation rate scale rate
THE COMBINATION PROBLEM Observation model Initial coordinates and velocities from each technique T Unknown ITRF initial coordinates and velocities Transformation parameters from ITRF system to technique (stacking) system GIVEN SOUGHT NUISANCE
INNER CONSTRAINTS FOR THE COMBINATION PROBLEM Change of ITRF reference system
INNER CONSTRAINTS FOR THE COMBINATION PROBLEM (Total) inner constraints initial orientation initialtranslation initialscale orientation rate translation rate scale rate
INNER CONSTRAINTS FOR THE COMBINATION PROBLEM Partial inner constraints – Coordinates & velocities initialorientation initialtranslation initialscale orientation rate translation rate scale rate
INNER CONSTRAINTS FOR THE COMBINATION PROBLEM Partial inner constraints – Coordinates & velocities initialorientation initialtranslation initialscale orientation rate translation rate scale rate Same as for the stacking problem !
INNER CONSTRAINTS FOR THE COMBINATION PROBLEM Partial inner constraints – Transformation parameters initialorientation initialtranslation initialscale orientation rate translation rate scale rate