Evaluation of the geometric stability and the accuracy potential of digital cameras — Comparing mechanical stabilisation versus parameterisation

Recent tests on the geometric stability of several digital cameras that were not designed for photogrammetric applications have shown that the accomplished accuracies in object space are either limited or that the accuracy potential is not exploited to the fullest extent. A total of 72 calibrations were calculated with four different software products for eleven digital camera models with different hardware setups, some with mechanical fixation of one or more parts. The calibration procedure was chosen in accord to a German guideline for evaluation of optical 3D measuring systems [VDI/VDE, VDI/VDE 2634 Part 1, 2002. Optical 3D Measuring Systems–Imaging Systems with Point-by-point Probing. Beuth Verlag, Berlin]. All images were taken with ringflashes which was considered a standard method for close-range photogrammetry. In cases where the flash was mounted to the lens, the force exerted on the lens tube and the camera mount greatly reduced the accomplished accuracy. Mounting the ringflash to the camera instead resulted in a large improvement of accuracy in object space. For standard calibration best accuracies in object space were accomplished with a Canon EOS 5D and a 35 mm Canon lens where the focusing tube was fixed with epoxy (47  μm maximum absolute length measurement error in object space). The fixation of the Canon lens was fairly easy and inexpensive resulting in a sevenfold increase in accuracy compared with the same lens type without modification. A similar accuracy was accomplished with a Nikon D3 when mounting the ringflash to the camera instead of the lens (52  μm maximum absolute length measurement error in object space). Parameterisation of geometric instabilities by introduction of an image variant interior orientation in the calibration process improved results for most cameras. In this case, a modified Alpa 12 WA yielded the best results (29  μm maximum absolute length measurement error in object space). Extending the parameter model with FiBun software to model not only an image variant interior orientation, but also deformations in the sensor domain of the cameras, showed significant improvements only for a small group of cameras. The Nikon D3 camera yielded the best overall accuracy (25  μm maximum absolute length measurement error in object space) with this calibration procedure indicating at the same time the presence of image invariant error in the sensor domain. Overall, calibration results showed that digital cameras can be applied for an accurate photogrammetric survey and that only a little effort was sufficient to greatly improve the accuracy potential of digital cameras.     

Nonlinear Adjustment of GPS Observations of Type Pseudo-Ranges

The nonlinear adjustment of GPS observations of type pseudo-ranges is performed in two steps. In step one a combinatorial minimal subset of observations is constructed which is rigorously converted into station coordinates by means of Groebner basis algorithm or the multipolynomial resultant algorithm. The combinatorial solution points in a polyhedron are reduced to their barycentric in step two by means of their weighted mean. Such a weighted mean of the polyhedron points in ℝ³ is generated via the Error Propagation law/variance-covariance propagation. The Fast Nonlinear Adjustment Algorithm (FNon Ad Al) has been already proposed by Gauss whose work was published posthumously and Jacobi (1841). The algorithm, here referred to as the Gauss-Jacobi Combinatorial algorithm, solves the over-determined GPS pseudo-ranging problem without reverting to iterative or linearization procedure except for the second moment (Variance-Covariance propagation). The results compared well with the solutions obtained using the linearized least squares approach giving legitimacy to the Gauss-Jacobi combinatorial procedure. ? 2002 Wiley Periodicals, Inc.     

Use of the L2C signal for inversions of GPS radio occultation data in the neutral atmosphere

Results from processing FORMOSAT-3/COSMIC radio occultations (RO) with the new GPS L2C signal acquired both in phase locked loop (PLL) and open loop (OL) modes are presented. Analysis of L2P, L2C, and L1CA signals acquired in PLL mode shows that in the presence of strong ionospheric scintillation not only L2P tracking, but also L1CA tracking often fails, while L2C tracking is most stable. The use of L2C improves current RO processing in the neutral atmosphere mainly by increasing the number of processed occultations (due to significant reduction in the number of L2 tracking failures) and marginally by a reduction in noise in statistics. The latter is due to the combination of reduced L2C noise (compared to L2P) and increased L1CA noise in those occultations where L2P would have failed. This result suggests application of OL tracking for L1CA and L2C signals throughout an entire occultation to optimally acquire RO data. Two methods of concurrent processing of L1CA and L2C RO signals are considered. Based on testing of individual occultations, these methods allow: (1) reduction in uncertainty of bending angles retrieved by wave optics in the lower troposphere and (2) reduction in small-scale residual errors of the ionospheric correction in the stratosphere.     

Motion Analysis in SAR Images of Unfocused Objects Using Time–Frequency Methods

Synthetic aperture radar (SAR) image formation processing assumes that the scene is stationary, and to focus an object, one coherently sums a large number of independent returns. Any target motion introduces phases that distort and/or translate the target's image. Target motion produces a smear primarily in the azimuth direction of the SAR image. Time-frequency (TF) modeling is used to analyze and correct the residual phase distortions. An interactive focusing algorithm based on TF modeling demonstrates how to correct the phase and to rapidly focus the mover. This is demonstrated on two watercraft observed in a SAR image. Then, two time-frequency representations (TFRs) are applied to estimate the motion parameters of the movers or refocus them or both. The first is the short-time Fourier transform, from which a velocity profile is constructed based on the length of the smear. The second TFR is the time-frequency distribution series, which is a robust derivative of the Wigner-Ville distribution that works well in this SAR environment. The smear is a modulated chirp, from which a velocity profile is plotted and the phase corrections are integrated to focus the movers. The relationship between these two methods is discussed. Both methods show good agreement on the example.     

TAZ-level variation in work trip mode choice between 1990 and 2000 and the presence of rail transit

Rail transit continues to be a popular alternative for cities as a tool for alleviating automobile congestion and for redeveloping areas into transit and pedestrian-friendly environments. Ideally, rail transit will draw trips away from cars, but the quantitative work that tests this notion has often been either case studies of neighborhoods, in which conclusions are tough to generalize, or citywide comparisons where important spatial variation is lost in aggregation. This study seeks to narrow this gap in the research by using multivariate analysis of covariance to isolate the effect of covariates and cities on changes in work trip mode choice at the traffic analysis zone (TAZ) level for nine cities between 1990 and 2000. The results suggest differences by city in the change at the TAZ level of the proportion of people driving alone and taking transit. Increases in transit usage were associated with cities that built rail transit, while increases in automobile commuting and decreases in transit usage were associated with cities that did not.     

Evaluation of the first GOCE static gravity field models using terrestrial gravity,vertical deflections and EGM2008 quasigeoid heights

Recently, four global geopotential models (GGMs) were computed and released based on the first 2 months of data collected by the Gravity field and steady-state Ocean Circulation Explorer (GOCE) dedicated satellite gravity field mission. Given that GOCE is a technologically complex mission and different processing strategies were applied to real space-collected GOCE data for the first time, evaluation of the new models is an important aspect. As a first assessment strategy, we use terrestrial gravity data over Switzerland and Australia and astrogeodetic vertical deflections over Europe and Australia as ground-truth data sets for GOCE model evaluation. We apply a spectral enhancement method (SEM) to the truncated GOCE GGMs to make their spectral content more comparable with the terrestrial data. The SEM utilises the high-degree bands of EGM2008 and residual terrain model data as a data source to widely bridge the spectral gap between the satellite and terrestrial data. Analysis of root mean square (RMS) errors is carried out as a function of (i) the GOCE GGM expansion degree and (ii) the four different GOCE GGMs. The RMS curves are also compared against those from EGM2008 and GRACE-based GGMs. As a second assessment strategy, we compare global grids of GOCE GGM and EGM2008 quasigeoid heights. In connection with EGM2008 error estimates, this allows location of regions where GOCE is likely to deliver improved knowledge on the Earth’s gravity field. Our ground truth data sets, together with the EGM2008 quasigeoid comparisons, signal clear improvements in the spectral band ~160–165 to ~180–185 in terms of spherical harmonic degrees for the GOCE-based GGMs, fairly independently of the individual GOCE model used. The results from both assessments together provide strong evidence that the first 2 months of GOCE observations improve the knowledge of the Earth’s static gravity field at spatial scales between ~125 and ~110 km, particularly over parts of Asia, Africa, South America and Antarctica, in comparison with the pre-GOCE-era.     

SIMD correlator library for GNSS software receivers

Software receivers have had a discernable impact on the GNSS research community. Often such receivers are implemented in a compiled programming language, such as C or C++. A software receiver must emulate the digital signal processing (DSP) algorithms executed on dedicated hardware in a traditional receiver. The DSP algorithms, most notably correlation, have a high computational cost; this burden precludes many software receivers from running in real time. However, the computational cost can be lessened by utilizing single instruction multiple data (SIMD) operations found on modern ×86 processors. The following demonstrates how C/C++ compatible code can be written to directly utilize the SIMD instructions. First, an analysis is carried out to demonstrate why real time operation is not possible when using traditional C/C++ code is carried out. Secondly a tutorial outlines how to write and insert ×86 assembly, with SIMD operations, into C/C++ code. Performance gains achieved via SIMD operations are then demonstrated, and pseudo code outlines how SIMD operations can be employed to perform correlation. Finally, a C/C++ compatible SIMD enabled arithmetic library is added to the GPS Toolbox for use in software receivers.     

Illumination effects on the differenced Normalized Burn Ratio's optimality for assessing fire severity

The influence of illumination effects on the optimality of the dNBR (differenced Normalized Burn Ratio) was evaluated for the case of the 2007 Peloponnese (Greece) wildfires using a pre/post-fire Landsat TM (Thematic Mapper) image couple. Well-illuminated pixels (south and south-east facing slopes) exhibited more optimal displacements in the bi-spectral feature space than more shaded pixels (north and north-west exposed slopes). Moreover, pixels experiencing a small image-to-image difference in illumination obtained a higher optimality than pixels with a relatively large difference in illumination. To correct for illumination effects, the c-correction method and a modified c-correction technique were applied. The resulting median dNBR optimality of uncorrected, c-corrected and modified c-correction data was respectively 0.58, 0.60 and 0.71 (differences significant for p < 0.001). The original c-correction method improved the optimality of badly illuminated pixels while deteriorating the optimality of well-illuminated pixels. In contrast, the modified c-correction technique improved the optimality of all the pixels while retaining the prime characteristic of topographic correction techniques, i.e. detrending the illumination–reflectance relationship. For a minority of the data, for shaded pixels and/or pixels with a high image-to-image difference in illumination, the original c-correction outperformed the modified c-correction technique. In this study conducted in rugged terrain and with a bi-temporal image acquisition scheme that deviated up to two months from the ideal anniversary date scheme the modified c-correction technique resulted in a more reliable change detection.