THE NEW DECADE: DIGITS AND DIVERSITY

The President of the Photogrammetric Society, Mr. I Newton, chaired a panel discussion concerned with those developments in photogrammetry which may take place during the 1990s.     

INSTRUMENTAL REQUIREMENTS FOR MODERN MAP REVISION

An examination of the need for map revision and the methods currently used allows the identification of problems specific to this operation. It is suggested that many instruments, old and new, suitable for new map production do not solve the problems associated with map revision. The main characteristics of instruments designed for map revision in the future are outlined, with the suggestion that the time has now come for instruments to be designed specifically for the map revision task, rather than practitioners having to devise revision procedures to suit existing equipment.     

A technique for modeling the Earth’s gravity field on the basis of satellite accelerations

A technique is proposed for Earths gravity field modeling on the basis of satellite accelerations that are derived from precise orbit data. The functional model rests on Newtons second law. The computational procedure is based on the pre-conditioned conjugate-gradient (PCCG) method. The data are treated as weighted average accelerations rather than as point-wise ones. As a result, a simple three-point numerical differentiation scheme can be used to derive them. Noise in the orbit-derived accelerations is strongly dependent on frequency. Therefore, the key element of the proposed technique is frequency-dependent data weighting. Fast convergence of the PCCG procedure is ensured by a block-diagonal pre-conditioner (approximation of the normal matrix), which is derived under the so-called Colombo assumptions. Both uninterrupted data sets and data with gaps can be handled. The developed technique is compared with other approaches: (1) the energy balance approach (based on the energy conservation law) and (2) the traditional approach (based on the integration of variational equations). Theoretical considerations, supported by a numerical study, show that the proposed technique is more accurate than the energy balance approach and leads to approximately the same results as the traditional one. The former finding is explained by the fact that the energy balance approach is only sensitive to the along-track force component. Information about the cross-track and the radial component of the gravitational potential gradient is lost because the corresponding force components do no work and do not contribute to the energy balance. Furthermore, it is shown that the proposed technique is much (possibly, orders of magnitude) faster than the traditional one because it does not require the computation of the normal matrix. Hints are given on how the proposed technique can be adapted to the explicit assembling of the normal matrix if the latter is needed for the computation of the model covariance matrix.Acknowledgments. Professor R. Klees is thanked for support of the project and for numerous fruitful discussions. The authors are also thankful to Dr. J. Kusche for useful remarks and to Dr. E. Schrama, his solid background in satellite geodesy proved to be very helpful. A large number of valuable comments were made by Dr. S.-C. Han, Dr. P. Schwintzer, and an anonymous reviewer; their contribution is greatly acknowledged. The satellite orbits used in the numerical study were kindly provided by Dr. P. Visser (Aerospace Department, Delft University of Technology). Access to the SGI Origin 3800 computer was provided by Stichting Nationale Computerfaciliteiten (NCF), grant SG-027.     

On the reliability of the VCV Matrix: A case study based on GAMIT and Bernese GPS Software

Commonly, the variance-covariance (VCV) matrix derived from GPS processing software underestimates the magnitude of the error, mainly due to the fact that physical correlations are normally neglected. The GAMIT and Bernese software packages serve the scientific community as important tools for GPS measurement processing and analyzing, especially in precise applications. Therefore, the reliability of the VCV matrices derived by the GAMIT and Bernese packages is of great importance. Formal accuracies derived from both software need to be scaled by applying a scaling factor (SF) that multiplies the software-derived formal errors. However, to the best of our knowledge, no standard approach approved by the GPS community exists. In this report, an analysis is carried out in order to test the reliability and the validity of the VCV matrices in both software, and to provide SFs needed to calculate the realistic accuracies reflecting the actual error levels. The method applied in this study allows deriving SFs for formal accuracies obtained from GAMIT and Bernese. The results attained from the time series of eight days for eight baselines (lengths of 20–415 km) indicate that the overall SF for GAMIT is more than 10 times smaller than for Bernese (1.9 and 23.0, respectively). Although no distance-dependent SF was detected in either case, the session-duration dependence was detected for the Bernese software, while no clear session-duration dependence was observed for the GAMIT. Furthermore, no receiver/antenna dependence could be deduced from the results of this analysis.     

Wet tropospheric effects on precise relative GPS height determination

Summary Considerable interest has been generated recently in the use of the Global Positioning System (GPS) for precise height determination. A major error source in these measurements is the propagation delay due to atmospheric water vapour. In order to achieve the high precisions required for such applications as absolute sea-level monitoring improvement of wet delay modelling is necessary. Results from a GPS campaign show a significant correlation (0.91) between the variability of the wet delay measured using a water vapour radiometer (WVR) at the Onsala site and the absolute value of the residual error in the height determination of a 134 km baseline from Onsala to Jönköping. This correlation indicates that the atmosphericvariability as inferred from the WVR data includes information on the quality of the GPS height estimate. During periods of high atmospheric activity, e.g., during the passage of a weather front, the use of a six-parameter gradient model reduces the spread for the vertical coordinate from 40 mm to 20 mm (with standard deviations of 17 mm and 9 mm respectively) over the 134 km baseline (less than 1 × 10^–7) using 8 hour data spans on 11 different days over a six month period.     

Spectrum estimation using maximum entropy and multiresolution considerations

Summary Given a sample autocovariance sequence of finite length for some observed random process, the spectrum estimation problem involves the extension of this sequence for the required Fourier transformation. The maximum entropy approach which is based on the optimal use of information contents, leads to a dual sequence of reflection coefficients with reciprocal spectrum of the process. The estimation of the maximum entropy spectrum implies results identical to those using autoregressive modeling in one dimension under appropriate white noise assumptions. In cases of a non-white noise component, the approach is generalized to an autoregressive-moving-average model. Recent developments in multiresolution analysis with spectral domain decompositions also offer possibilities of subband spectrum estimation for specific applications. Using a simulated data sequence with two close frequencies, the estimated spectrum from a two-level decomposition with autoregressive modeling shows better resolution than with conventional processing. Geodetic and geophysical applications are briefly indicated.     

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     

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     

The hierarchical tessellation model and its use in spatial analysis

The hierarchid tessellation model belongs to a class of spatial data models based on the recursive decomposition of space. The quadtree is one such tessellation and is characterized by square cells and a 1:4 decomposition ratio. To relax these constraints in the tessellation, a generalized hierarchical tessellation data model, called Adaptive Recursive Tessellations (ART), has been proposed. ART increases flexibility in the tessellation by the use of rectangular cells and variable decomposition ratios. In ART, users can specify cell sizes which are intuitively meaningful to their applications, or which can reflect the scales of data. ART is implemented in a data structure called Adaptive Recursive Run-Encoding (ARRE), which is a variant of two-dimensional run-encoding whose running path can vary with the different tessellation structures incorporated in an ART model. Given the recognition of the benefits of implementing statistical spatial analysis in GIS, the use of hierarchical tessellation models such as ART in spatial analysis is discussed. Three examples are introduced to show how ART can: (1) be applied to solve the quadrat size problem in quadrat analysis of point patterns; (2) act as the data model in the variable resolution block kriging technique for geostatistical data to reduce variation in kriging error; and (3) facilitate the evaluation of spatial autocorrelation for area data at multiple map resolutions via the construction of a connectivity matrix for calculating spatial autocorrelation indices based on ARRE.