Transformation of Lambert Conic Conformal Coordinates from a Global Datum to a Local Datum

This paper treats the problem of how to transform from global datum, for example, from the International Terrestrial Reference System (ITRS), to a local datum, for example, regional or national, for the practical case of the Lambert projection of the sphere or the ellipsoid-of-revolution to the cone. We design the two projection constants n(ϕ₁, ϕ₂) and m(ϕ₁) for the Universal Lambert Conic projection of the ellipsoid-of-revolution. The task to transform from a global datum with respect to the ellipsoid-of-revolution E_A,B² to local datum with respect to the alternative ellipsoid-of-revolution E_a,b², without local ellipsoidal height, is solved by an extended numerical example. Ideas in this paper could be of interest to those working with maps and coordinates transformation from global geodetic datum to local geodetic datum and vice versa, under the Universal Lambert Conic projection, and applicable to precise positioning and navigation, boundary demarcation and determination in the marine environment.     

Methodology and use of tensor invariants for satellite gravity gradiometry

Although its use is widespread in several other scientific disciplines, the theory of tensor invariants is only marginally adopted in gravity field modeling. We aim to close this gap by developing and applying the invariants approach for geopotential recovery. Gravitational tensor invariants are deduced from products of second-order derivatives of the gravitational potential. The benefit of the method presented arises from its independence of the gradiometer instrument’s orientation in space. Thus, we refrain from the classical methods for satellite gravity gradiometry analysis, i.e., in terms of individual gravity gradients, in favor of the alternative invariants approach. The invariants approach requires a tailored processing strategy. Firstly, the non-linear functionals with regard to the potential series expansion in spherical harmonics necessitates the linearization and iterative solution of the resulting least-squares problem. From the computational point of view, efficient linearization by means of perturbation theory has been adopted. It only requires the computation of reference gravity gradients. Secondly, the deduced pseudo-observations are composed of all the gravitational tensor elements, all of which require a comparable level of accuracy. Additionally, implementation of the invariants method for large data sets is a challenging task. We show the fundamentals of tensor invariants theory adapted to satellite gradiometry. With regard to the GOCE (Gravity field and steady-state Ocean Circulation Explorer) satellite gradiometry mission, we demonstrate that the iterative parameter estimation process converges within only two iterations. Additionally, for the GOCE configuration, we show the invariants approach to be insensitive to the synthesis of unobserved gravity gradients.     

The gravitational field of topographic-isostatic masses and the hypothesis of mass condensation II-the topographic-isostatic geoid

In Part I we focussed on a convergent representation of the gravitational potential generated bytopographic masses on top of the equipotential surface atMean Sea Level, thegeoid, and by those masses which compensate topography. Topographic masses have also been condensated, namely represented by a single layer. Part II extends the computation of the gravitational field of topographic-isostatic masses by a detailed analysis of itsforce field in terms ofvector-spherical harmonic functions. In addition, the discontinuous mass-condensated topographic gravitational force vector (head force) is given. Once we identify theMoho discontinuity asone interface of isostatically compensated topographical masses, we have computed the topographic potential and the gravitational potential which is generated by isostatically compensated masses atMean Sea Level, the geoid, and illustrated by various figures of geoidal undulations. In comparison to a data oriented global geoid computation ofJ. Engels (1991) the conclusion can be made that the assumption of aconstant crustal mass density, the basic condition for isostatic modeling, does not apply. Insteaddensity variations in the crust, e.g. betweenoceanic and continental crust densities, have to be introduced in order to match the global real geoid and its topographic-isostatic model. The performed analysis documents that thestandard isostatic models based upon aconstant crustal density areunreal.     

Probability distribution of eigenspectra and eigendirections of a twodimensional,symmetric rank two random tensor

Let there be given a twodimensional symmetric rank two tensor of random type (examples:strain, stress) which is either directly observed or indirectly estimated from observations by an adjustment procedure. Under the assumption of normalityof tensor components we compute the joint probability density functionas well as the marginal probability density functionsof its eigenspectra (eigenvalues) and eigendirections (orientation parameters). Due to the nonlinearity of the relation between eigenspectra-eigendirections and the random tensor components, via the inverse nonlinear error propagationbiases and aliases of their first and centralized second moments (mean value, variance-covariance) are expressed in terms of Jacobianand Hessianmatrices. The joint probability density function and the first and second moments thus form the fundamental of hypothesis testing and qualify control of eigenspectra (eigenvalues, principal components) and eigendirections (orientation parameters, eigenvectors, principial direction) of a twodimensional, symmetric rank two random tensor.     

Closed-form solution of P4P or the three-dimensional resection problem in terms of Möbius barycentric coordinates

 The perspective 4 point (P4P) problem - also called the three-dimensional resection problem - is solved by means of a new algorithm: At first the unknown Cartesian coordinates of the perspective center are computed by means of M?bius barycentric coordinates. Secondly these coordinates are represented in terms of observables, namely space angles in the five-dimensional simplex generated by the unknown point and the four known points. Substitution of M?bius barycentric coordinates leads to the unknown Cartesian coordinates (2.8)–(2.10) of Box 2.2. The unknown distances within the five-dimensional simplex are determined by solving the Grunert equations, namely by forward reduction to one algebraic equation (3.8) of order four and backward linear substitution. Tables 1.–4. contain a numerical example. Finally we give a reference to the solution of the 3 point (P3P) problem, the two-dimensional resection problem, namely to the Ansermet barycentric coordinates initiated by C.F. Gau? (1842), A. Schreiber (1908) and A.␣Ansermet (1910). Received: 05 March 1996; Accepted: 15 October 1996     

Closed form solution to the twin P4P or the combined three dimensional resection-intersection problem in terms of Möbius barycentric coordinates

The twin perspective 4 point (twin P4P) problem – also called the combined three dimensional resection-intersection problem – is the problem of finding the position of a scene object from 4 correspondence points and a scene stereopair. While the perspective centers of the left and right scene image are positioned by means of a double three dimensional resection, the position of the scene object imaged on the left and right photograph is determined by a three dimensional intersection based upon given resected perspective centers. Here we present a new algorithm solving the twin P4P problem by means of M?bius barycentric coordinates. In the first algorithmic step we determine the distances between the perspective centers and the unknown intersected point by solving a linear system of equations. Typically, area elements of the left and right image build up the linear equation system. The second algorithmic step allows for the computation of the M?bius barycentric coordinates of the unknown intersected point which are thirdly converted into three dimensional object space coordinates {X,Y,Z} of the intersected point. Typically, this three-step algorithm based upon M?bius barycentric coordinates takes advantage of the primary double resection problem from which only distances from four correspondence points to the left and right perspective centre are needed. No orientation parameters and no coordinates of the left and right perspective center have to be made available. Received 1 May 1996; Accepted 13 September 1996     

Construction of Green's function to the external Dirichlet boundary-value problem for the Laplace equation on an ellipsoid of revolution

Green's function to the external Dirichlet boundary-value problem for the Laplace equation with data distributed on an ellipsoid of revolution has been constructed in a closed form. The ellipsoidal Poisson kernel describing the effect of the ellipticity of the boundary on the solution of the investigated boundary-value problem has been expressed as a finite sum of elementary functions which describe analytically the behaviour of the ellipsoidal Poisson kernel at the singular point ψ = 0. We have shown that the degree of singularity of the ellipsoidal Poisson kernel in the vicinity of its singular point is of the same degree as that of the original spherical Poisson kernel. Received: 4 June 1996 / Accepted: 7 April 1997     

The spacetime gravitational field of a deformable body

The high-resolution analysis of orbit perturbations of terrestrial artificial satellites has documented that the eigengravitation of a massive body like the Earth changes in time, namely with periodic and aperiodic constituents. For the space-time variation of the gravitational field the action of internal and external volume as well as surface forces on a deformable massive body are responsible. Free of any assumption on the symmetry of the constitution of the deformable body we review the incremental spatial (“Eulerian”) and material (“Lagrangean”) gravitational field equations, in particular the source terms (two constituents: the divergence of the displacement field as well as the projection of the displacement field onto the gradient of the reference mass density function) and the `jump conditions' at the boundary surface of the body as well as at internal interfaces both in linear approximation. A spherical harmonic expansion in terms of multipoles of the incremental Eulerian gravitational potential is presented. Three types of spherical multipoles are identified, namely the dilatation multipoles, the transport displacement multipoles and those multipoles which are generated by mass condensation onto the boundary reference surface or internal interfaces. The degree-one term has been identified as non-zero, thus as a “dipole moment” being responsible for the varying position of the deformable body's mass centre. Finally, for those deformable bodies which enjoy a spherically symmetric constitution, emphasis is on the functional relation between Green functions, namely between Fourier-/ Laplace-transformed volume versus surface Love-Shida functions (h(r),l(r) versus h ^′(r),l ^′(r)) and Love functions k(r) versus k ^′(r). The functional relation is numerically tested for an active tidal force/potential and an active loading force/potential, proving an excellent agreement with experimental results. Received: December 1995 / Accepted: 1 February 1997