A MATLAB toolbox for attitude determination with GPS multi-antenna systems

In this paper a MATLAB toolbox for determining the attitude of a rigid platform by means of multiple non-dedicated antennas using global positioning system is presented. The programs embedded in this toolbox cover the RINEX data analysis, single point positioning, differential positioning, coordinate conversion, attitude determination, and other auxiliary functions. After forming the baselines through double-differenced (carrier phase smoothed) code observables, the attitude parameters are obtained by applying the direct attitude computation and the least squares attitude estimation. The theoretical background is summarized, and some hints regarding the software implementation are given in the paper. Moreover, improvements yielding an expanded functionality are proposed.

Network-based geometry-free three carrier ambiguity resolution and phase bias calibration

Continuously operating reference stations (CORS) are increasingly used to deliver real-time and near-real-time precise positioning services on a regional basis. A CORS network-based data processing system uses either or both of the two types of measurements: (1) ambiguity-resolved double-differenced (DD) phase measurements, and (2) phase bias calibrated zero-differenced (ZD) phase measurements. This paper describes generalized, network-based geometry-free models for three carrier ambiguity resolution (TCAR) and phase bias estimation with DD and ZD code and phase measurements. First, the geometry-free TCAR models are constructed with two Extra-Widelane (EWL)/Widelane (WL) virtual observables to allow for rapid ambiguity resolution (AR) for DD phase measurements without distance constraints. With an ambiguity-resolved WL phase measurement and the ionospheric estimate derived from the two EWL observables, an additional geometry-free equation is formed for the third virtual observable linearly independent of the previous two. AR with the third geometry-free model requires a longer period of observations for averaging than the first two, but is also distance-independent. A more general formulation of the geometry-free model for a baseline or network is also introduced, where all the DD ambiguities can be more rigorously resolved using the LAMBDA method. Second, the geometry-free models for calibration of three carrier phase biases of ZD phase measurements are similarly defined for selected virtual observables. A network adjustment procedure is then used to improve the ZD phase biases with known DD integer constraints. Numerical results from experiments with 24-h dual-frequency GPS data from three US CORS stations baseline lengths of 21, 56 and 74 km confirm the theoretical predictions concerning AR reliability of the network-based geometry-free algorithms.

Critical network infrastructure analysis: interdiction and system flow

Effective management of critical network infrastructure requires the assessment of potential interdiction scenarios. Optimization approaches have been essential for identifying and evaluating such scenarios in networked systems. Although a primary function of any network is the distribution of flow between origins and destinations, the complexity and difficulty of mathematically abstracting interdiction impacts on connectivity or flow has been a challenge for researchers. This paper presents an optimization approach for identifying interdiction bounds with respect to connectivity and/or flow associated with a system of origins and destinations. Application results for telecommunications flow are presented, illustrating the capabilities of this approach.

---

An ANN–RTS smoother scheme for accurate INS/GPS integrated attitude determination

Digital mobile mapping, the method that integrates digital imaging with direct geo-referencing, has developed rapidly over the past 15 years. The Kalman filter (KF) is considered an optimal estimation tool for real-time INS/GPS integrated kinematic positioning and orientation determination. However, the accuracy requirements of general mobile mapping applications cannot be easily achieved even when using the KF scheme. Therefore, this study proposes an intelligent scheme combining ANN and RTS backward smoother to overcome the limitations of KF and to enhance the overall accuracy of attitude determination for tactical grade and MEMS INS/GPS integrated systems.

Impact of systematic errors on precise long-baseline kinematic GPS positioning

Many kinematic GPS applications rely on high accuracy, which usually requires the ambiguities to be fixed. Normally, a reference station in the rover’s vicinity is needed for successful ambiguity resolution. Alternatively, a network surrounding the rover and allowing one to derive area correction parameters is needed. Unfortunately, both approaches are not feasible in certain situations. This paper is a contribution to precise kinematic positioning over long baselines. Atmospheric refraction becomes critical in the error budget, but progress has been made to use numerical weather models to derive tropospheric corrections, for instance. The spatial correlation of both ionospheric and tropospheric propagation delays is investigated in this paper and special attention is paid on the systematic error behavior of tropospheric refraction. The principles developed are applied to an extended reliability test of the ambiguities. Finally, it is demonstrated in positioning experiments that kinematic positioning retrieval with fixed ambiguities is actually possible for baselines between 150 and 300 km with an accuracy of approximately 2 cm in post-mission processing.

---

Ultra-tight GPS/INS/PL integration: a system concept and performance analysis

The architecture of the ultra-tight GPS/INS/PL integration is the key to its successful performance; the main feature of this architecture is the Doppler feedback to the GPS receiver tracking loops. This Doppler derived from INS, when integrated with the carrier tracking loops, removes the Doppler due to vehicle dynamics from the GPS/PL signal thereby achieving a significant reduction in the carrier tracking loop bandwidth. The bandwidth reduction provides several advantages such as: improvement in anti-jamming performance, and increase in post correlated signal strength which in turn increases the dynamic range and accuracy of measurements. Therefore, any degradation in the derived Doppler estimates will directly affect the tracking loop bandwidth and hence its performance. The quadrature signals from the receiver correlator, I (in-phase) and Q (quadrature), form the measurements, whereas the inertial sensor errors, position, velocity and attitude errors form the states of the complementary Kalman filter. To specify a reliable measurement model of the filter for this type of integrated system, a good understanding of GPS/PL signal characteristics is essential. It is shown in this paper that phase and frequency errors are the variables that relate the measurements and the states in the Kalman filter. The main focus of this paper is to establish the fundamental mathematical relationships that form the measurement model, and to show explicitly how the system error states are related to the GPS/PL signals. The derived mathematical relationships encapsulated in a Kalman filter, are tested by simulation and shown to be valid.

---

GPS/陀螺组合测姿中整周模糊度的快速解算

采用方向余弦矩阵描述姿态,建立GPS/陀螺组合姿态确定系统模型,由矩阵Kalman滤波方法解算整周模糊度的浮点解,然后再利用MCLambda方法得到整周模糊度固定解。仿真实验结果表明,附加方向余弦矩阵约束的Kalman滤波方法可以有效地提高整周模糊度浮点解的精度,使得整周模糊度的固定成功率和效率均得到提高,尤其是在GPS观测条件较差的情况下。     

Ionospheric delay corrections for single-frequency GPS receivers over Europe using tomographic mapping

The majority of navigation satellites receivers operate on a single frequency and experience a positioning error due to the ionospheric delay. This can be compensated for using a variety of approaches that are compared in this paper. The study focuses on the last solar maximum. A 4D tomographic imaging technique is used to map the ionospheric electron density over the European region during 2002 and 2003. The electron density maps are then used to calculate the excess propagation delay on the L1 frequency experienced by GPS receivers at selected locations across Europe. The excess delay is applied to correct the pseudo-range single frequency observations at each location and the improvements to the resulting positioning are calculated. The real-time tomographic technique is shown to give navigation solutions that are better than empirical modelling methods and approach the accuracy of the full dual-frequency solution. The improvements in positioning accuracy vary from day to day depending on ionospheric conditions but can be up to 25 m during mid-day during these solar maximum conditions at European mid-latitudes.

Rectilinear minimax hub location problems

The minimax hub location problem sites a facility to minimize the maximum weighted interaction cost between pairs of fixed nodes. In this paper, distances are represented by a rectilinear norm and may be suited to factory layout or street network problems. The problem is already well known (in 2-D) as the round trip location problem and is extended to 3-D in this paper. One rationale for the solution method is based on an extension of the geometric arguments used to solve the minimax single facility location problem. Suppose a budget is provided for interactions, and that each interaction must be accomplished for no more than this cost. The algorithm uses a bi-section search for the feasible budget until it finds the expenditure needed to provide for these flows. The extension in the present paper is that the nodes are permitted to be on different layers (levels). This 3-D version of the problem appears to be a new variant of the hub model. The models and solution techniques developed in the paper are illustrated using a small 55 node problem. Because of a relatively efficient implementation of the bi-section search, the algorithm in 2-D and 3-D is also applied successfully to a 550 node problem.

Pitch and roll attitude estimation of a small-scaled helicopter using single antenna GPS with gyroscopes

A method is presented for estimating the roll and pitch attitude of a small-scaled unmanned helicopter based on the velocity measurements of the global positioning system (GPS). The small-scaled helicopter is a radio controlled (RC) model which is readily available and affordable for academic laboratories as a research platform. Only one single antenna GPS receiver is equipped on the RC helicopter to acquire the velocity measurements needed for the attitude estimation. The velocity information is recorded by the onboard computer for post-processing. An attitude and heading reference system (AHRS) is used to provide the reference attitudes. The required angular rates and heading for this study are also given by the gyroscopes and compass of the AHRS for the sake of system’s simplification. The Kalman filter is applied to estimate the helicopter’s accelerations by using the GPS velocity measurements. The estimated accelerations form the fundamental elements of synthesizing the pseudo-roll and the pseudo-pitch. With some legitimate simplifications and assumptions, the relation between the helicopter’s attitudes and the accelerations estimated from the GPS velocity measurements can be developed. Furthermore, to enhance the accuracy of the pseudo-attitudes, the angular rates acquired from the gyroscopes are incorporated into the estimation algorithm of pseudo-attitudes by using a complementary filter.

Incorporating spatial variation in housing attribute prices: a comparison of geographically weighted regression and the spatial expansion method

Hedonic house price models typically impose a constant price structure on housing characteristics throughout an entire market area. However, there is increasing evidence that the marginal prices of many important attributes vary over space, especially within large markets. In this paper, we compare two approaches to examine spatial heterogeneity in housing attribute prices within the Tucson, Arizona housing market: the spatial expansion method and geographically weighted regression (GWR). Our results provide strong evidence that the marginal price of key housing characteristics varies over space. GWR outperforms the spatial expansion method in terms of explanatory power and predictive accuracy.

Correction of apparent position shifts caused by GNSS antenna changes

Antenna changes at GNSS reference stations frequently produce discontinuities in the coordinate time series. These apparent position shifts are mainly caused by changes of carrier-phase multipath effects and different errors in the antenna phase center corrections. A monitoring method was developed and successfully tested, which requires additional GNSS observations from a local, temporary reference station. Changes of carrier-phase measurement errors due to the antenna change are determined and stored in L1 and L2 phase maps. These phase maps provide corrections to be applied either to the observation data obtained before the antenna change or to the observation data obtained after the antenna change. The observation corrections are able to remove coordinate discontinuities independent of the selected coordinate estimation algorithm.

A weighted difference barrier method in landscape genetics

Identifying barriers of species and characterize their effects on spatial distribution provide essential information to research in landscape genetics. We propose a weighted difference barrier (WDB) method as an alternative to maximum difference barriers (MDB), and to initiate and integrate more spatial modeling and methods into the problem solving process. Overall, WDB provides quick and straightforward improvements to the drawbacks of MDB. WDB integrates more sample location relationships into the barrier construction and reveals potential barriers that would otherwise go undetected. WDB incorporates both within group and between group genetic information, and delineates the barriers as a more complex pattern.

Simplified equivalent multiple baseline solutions with elevation-dependent weights

Since the assumption of all stations tracking the same satellites with identical weights was previously employed by Shen and Xu (GPS Solut 12:99–108, 2008) to derive the simplified GNSS single- and double-differenced equivalent equations, this supplementary paper expands these simplified equations in the case of each station tracking different satellites with elevation-dependent weights. Numerical experiments are performed to demonstrate the computational efficiency of the simplified equivalent algorithm relative to the traditional method in various scenarios of multi-baseline solutions with tracking different satellites. The fast computational speed of the simplified equivalent algorithm will potentially benefit the local, regional and even global GNSS multi-baseline solutions as well as the combined GNSS application.

Specifying a joint space- and time-lag using a bivariate Poisson distribution

Separate space- or time-lags have been considered regularly in data analyses; as space–time models are more recently being studied extensively in data analytic fashion, joint estimation of both lags has to be considered explicitly. This paper addresses this issue, taking into special consideration parametric parsimony together with specification richness; use of the bivariate Poisson frequency distribution is advocated and applied to an empirical case. The relation of this approach to random effects specifications is investigated. Data for Belgian regional products constitute the empirical case study.

Hierarchical spatial interaction among the Italian regions: a nonlinear relative dynamics approach

Regional interactions and spillover effects should be considered as important factors in growth analysis of regional economies. Using modified versions of the Dendrinos–Sonis model, this paper analyses the spatial hierarchical system of Italy. The interaction among Italian regions is considered at three different levels of spatial aggregation, the NUTS-1, NUTS-2 and NUTS-3 levels. Compared to more popular spatial econometric approaches, the Dendrinos–Sonis model and its extensions provide greater flexibility in the way interaction between regions is handled but the results strongly depend on the choice of a reference region.

Latent lifestyle preferences and household location decisions

Lifestyle, indicating preferences towards a particular way of living, is a key driver of the decision of where to live. We employ latent class choice models to represent this behavior, where the latent classes are the lifestyles and the choice model is the choice of residential location. Thus, we simultaneously estimate lifestyle groups and how lifestyle impacts location decisions. Empirical results indicate three latent lifestyle segments: suburban dwellers, urban dwellers, and transit-riders. The suggested lifestyle segments have intriguing policy implications. Lifecycle characteristics are used to predict lifestyle preferences, although there remain significant aspects that cannot be explained by observable variables.

The ecological fallacy in a time series context: evidence from Spanish regional unemployment rates

The ecological fallacy (EF) is a common problem regional scientists have to deal with when using aggregated data in their analyses. Although there is a wide number of studies considering different aspects of this problem, little attention has been paid to the potential negative effects of the EF in a time series context. Using Spanish regional unemployment data, this paper shows that EF effects are not only observed at the cross-section level, but also in a time series framework. The empirical evidence obtained shows that analytical regional configurations are the least susceptible to time effects relative to both normative and random regional configurations, while normative configurations are an improvement over random ones.

A sampling approach to estimate the log determinant used in spatial likelihood problems

Likelihood-based methods for modeling multivariate Gaussian spatial data have desirable statistical characteristics, but the practicality of these methods for massive georeferenced data sets is often questioned. A sampling algorithm is proposed that exploits a relationship involving log-pivots arising from matrix decompositions used to compute the log determinant term that appears in the model likelihood. We demonstrate that the method can be used to successfully estimate log-determinants for large numbers of observations. Specifically, we produce an log-determinant estimate for a 3,954,400 by 3,954,400 matrix in less than two minutes on a desktop computer. The proposed method involves computations that are independent, making it amenable to out-of-core computation as well as to coarse-grained parallel or distributed processing. The proposed technique yields an estimated log-determinant and associated confidence interval.

Show me the code: spatial analysis and open source

This paper considers the intersection of academic spatial analysis with the open source revolution. Its basic premise is that the potential for cross-fertilization between the two is rich, yet some misperceptions about these two communities pose challenges to realizing these opportunities. The paper provides a primer on the open source movement for academicians with an eye towards correcting these misperceptions. It identifies a number of ways in which increased adoption of open source practices in spatial analysis can enhance the development of the next generation of tools and the wider practice of scientific research and education.