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.     

Reviewers of the Journal of Geodesy for Volume 79

Acknowledgement for Referees

Reviewers of the Journal of Geodesy for Volume 79     

Effect of signal contamination in matched-filter detection of the signal on a cluttered background

To derive a matched filter for detecting a weak target signal in a hyperspectral image, an estimate of the band-to-band covariance of the target-free background scene is required. We investigate the effects of including some of the target signal in the background scene. Although the covariance is contaminated by the presence of a target signal (there is increased variance in the direction of the target signature), we find that the matched filter is not necessarily affected. In fact, if the variation in plume strength is strictly uncorrelated with the variation in background spectra, the matched filter and its signal-to-clutter ratio (SCR) performance will not be impaired. While there is little a priori reason to expect significant correlation between the plume and the background, there usually is some residual correlation, and this correlation leads to a suppressing effect that limits the SCR obtainable even for strong plumes. These effects are described and quantified analytically, and the crucial role of this correlation is illustrated with some numerical examples using simulated plumes superimposed on real hyperspectral imagery. In one example, we observe an order-of-magnitude loss in SCR for a matched filter based on the contaminated covariance.     

Plate kinematics of Nubia–Somalia using a combined DORIS and GPS solution

We have used up to 12 years of data to assess DORIS performance for geodynamics applications. We first examine the noise characteristics of the DORIS time-series of weekly station coordinates to derive realistic estimates of velocity uncertainties. We find that a combination of white and flicker noise best explains the DORIS time-series noise characteristics. Second, weekly solutions produced by the Institut Géographique National/Jet Propulsion Laboratory (IGN/JPL) DORIS Analysis Centre are combined to derive a global velocity field. This solution is combined with two independent GPS solutions, including 11 sites on Nubia and 5 on the Somalia plate. The combination indicates that DORIS horizontal velocities have an average accuracy of 3 mm/year, with best-determined sites having velocity accuracy better than 1 mm/year (one-sigma levels). Using our combined velocity field, we derive an updated plate kinematics model with a focus on the Nubia–Somalia area. Including DORIS data improves the precision of the angular velocity vector for Nubia by 15%. Our proposed model provides robust bounds on the maximum opening rates along the East African Rift (4.7–6.7 mm/year). It indicates opening rates 15 and 7% slower than values predicted by NUVEL-1A for the southern Atlantic Ocean and Indian Ocean, respectively. These differences are likely to arise from the fact that NUVEL-1A considered Africa as a single non-deforming plate, while here we use a more refined approach.     

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.     

Study of the optimum time span for distinguishing <Emphasis Type="Italic">Avena ludoviciana</Emphasis> from wheat crop based on their spectral reflectance characteristics

The study to establish the optimum time span for distinguishing Avena ludoviciana from wheat crop based on their spectral signatures was carried out at Student’s Research Farm, Department of Agronomy during 2006–07 and 2007–08. The experimental sites during both the seasons were sandy loam in texture, with normal soil reaction and electrical conductivity, low in organic carbon and available nitrogen and medium in available phosphorus and potassium. The experiment was laid out in randomized block design with four replications and consisting of twelve treatments comprising 0, 10, 15, 25, 50, 75, 100, 125, 150, 200, 250 plants m⁻² and a pure Avena ludoviciana plot (T_max). The results revealed that in all the treatments irrespective of wheat and weeds, the red reflectance (%) value decreased from 34 to 95 DAS (days after sowing) in 2006–07 and 45 DAS to 100 DAS during 2007–08, and thereafter a sharp increase was observed in all the treatments. This trend might be due to increased chlorophyll index after 34 DAS as red reflectance was reduced by chlorophyll absorption. Among all the treatments, T_max (Pure Avena ludoviciana plot) had the highest red reflectance and T₀ (Pure wheat plot) had a lowest value of red reflectance during both the years. The highest value of IR reflectance was obtained at 95 DAS (2006–07) and 70 DAS (2007–08) in all the treatments. IR reflectance of wheat crop ranged between 24.61 and 61.21 per cent during 2006–07 and 27.33 and 67.3 per cent during 2007–08. However, IR reflectance values declined after 95 DAS and 70 DAS up to harvesting during 2006–07 and 2007–08. This lower reflectance may have been due to the onset of senescence. The highest RR and NDVI values were recorded under pure wheat treatment and minimum under pure weed plots. This may be due to dark green colour and better vigor of the wheat as compared to Avena ludoviciana. It was observed that by using RR and NDVI, pure wheat can be distinguished from pure populations of Avena ludoviciana after 34 DAS and different levels of weed populations can be discriminated amongst themselves from 68 DAS up to 107 DAS during both the years of investigation.     

Combining EGM2008 and SRTM/DTM2006.0 residual terrain model data to improve quasigeoid computations in mountainous areas devoid of gravity data

A global geopotential model, like EGM2008, is not capable of representing the high-frequency components of Earth’s gravity field. This is known as the omission error. In mountainous terrain, omission errors in EGM2008, even when expanded to degree 2,190, may reach amplitudes of 10 cm and more for height anomalies. The present paper proposes the utilisation of high-resolution residual terrain model (RTM) data for computing estimates of the omission error in rugged terrain. RTM elevations may be constructed as the difference between the SRTM (Shuttle Radar Topography Mission) elevation model and the DTM2006.0 spherical harmonic topographic expansion. Numerical tests, carried out in the German Alps with a precise gravimetric quasigeoid model (GCG05) and GPS/levelling data as references, demonstrate that RTM-based omission error estimates improve EGM2008 height anomaly differences by 10 cm in many cases. The comparisons of EGM2008-only height anomalies and the GCG05 model showed 3.7 cm standard deviation after a bias-fit. Applying RTM omission error estimates to EGM2008 reduces the standard deviation to 1.9 cm which equates to a significant improvement rate of 47%. Using GPS/levelling data strongly corroborates these findings with an improvement rate of 49%. The proposed RTM approach may be of practical value to improve quasigeoid determination in mountainous areas without sufficient regional gravity data coverage, e.g., in parts of Asia, South America or Africa. As a further application, RTM omission error estimates will allow refined validation of global gravity field models like EGM2008 from GPS/levelling data.     

Kinematic precise point positioning at remote marine platforms

Precise kinematic differential positioning using the global positioning system (GPS) at a marine platform usually requires a relatively short distance (e.g. <500 km) to a land-based reference station. As an alternative, precise point positioning (PPP) is normally considered free from this limiting requirement. However, due to the prerequisite of network-based satellite products, PPP at a remote marine platform may still be affected by its distance to the reference network. Hence, this paper investigates this scenario by configuring rings of reference stations with different radii centered on a to-be-positioned marine platform. Particularly, we applied ambiguity resolution at a single station to PPP by estimating uncalibrated phase delays (UPDs). We used three rings of reference stations centered on a vessel, with radii of roughly 900, 2,000 and 3,600 km, to determine satellite clocks and UPDs independently. For comparison, we also performed differential positioning based on a single reference station with baseline lengths of about 400, 1,700 and 2,800 km. We demonstrate that, despite the increasing ring-network radius to a few 1,000 km, the overall change in accuracy of the satellite clocks that are used at the vessel is smaller than 0.02 ns, and the RMS values of differences between the three sets of narrow-lane UPD estimates are around 0.05 cycles only. Moreover, the kinematic positioning accuracy of PPP is affected by the increasing ring-network radius, but can still achieve several centimeters after ambiguity resolution when the vessel is over a few 1,000 km away from the ring network, showing better performance than that of differential positioning. Therefore, we propose that ambiguity-fixed PPP can be used at remote marine platforms that support precise oceanographic and geophysical applications in open oceans.     

Comparison of approaches of geographic partitioning for data anonymization

Journal of Geographical Systems - Given the large volumes of detailed data now being collected, there is a high demand for the release of this data for research purposes. In particular,...