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.

---

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.

A copula-based closed-form binary logit choice model for accommodating spatial correlation across observational units

This study focuses on accommodating spatial dependency in data indexed by geographic location. In particular, the emphasis is on accommodating spatial error correlation across observational units in binary discrete choice models. We propose a copula-based approach to spatial dependence modeling based on a spatial logit structure rather than a spatial probit structure. In this approach, the dependence between the logistic error terms of different observational units is directly accommodated using a multivariate logistic distribution based on the Farlie-Gumbel-Morgenstein (FGM) copula. The approach represents a simple and powerful technique that results in a closed-form analytic expression for the joint probability of choice across observational units, and is straightforward to apply using a standard and direct maximum likelihood inference procedure. There is no simulation machinery involved, leading to substantial computation gains relative to current methods to address spatial correlation. The approach is applied to teenagers’ physical activity participation levels, a subject of considerable interest in the public health, transportation, sociology, and adolescence development fields. The results indicate that failing to accommodate heteroscedasticity and spatial correlation can lead to inconsistent and inefficient parameter estimates, as well as incorrect conclusions regarding the elasticity effects of exogenous variables.

An exploratory analysis of hierarchical spatial interaction: the case of regional income shares in Indonesia

Regional interaction is generally understood as implying movement between regions at the same level of an hierarchy. This paper extends the notion to include an hierarchical system, thus facilitating the consideration of vertical interaction in the analysis of regional interaction. Obviously, vertical interaction is not altogether a new concept. One could find this concept in many analyses related to national-local or federal-state relationships. This paper treats hierarchy in a strict sense; spatial units are interacting one with another within, but not between, presumed super regions. A particular example drawn in this paper is the province-region relationships for the case of Indonesia. In this framework, provinces form a region, and regions form the nation. The Dendrinos–Sonis model is then used as the basis for measuring the hierarchical spatial interaction in Indonesia. The analysis will explore the degree to which complementarity and competitive interaction revealed at one level in the hierarchy persist at lower or higher levels.

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.

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.

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.

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.

Total electron content processing from GPS observations to facilitate ionospheric modeling

With the increasing global distribution of high rate dual-frequency global positioning system (GPS) receivers, the production of a real-time atmospheric constituent definition, total electron content (TEC), has become a beneficial contributor to the modeling applications used in the assessment of GPS position accuracy and the composition of the ionosphere, plasmasphere, and troposphere. Historically, TEC measurements have been obtained through post processing techniques to produce the quality of data necessary for modeling applications with rigorous error estimate requirements. These procedures necessitated the collection of large volumes of data to address the various abnormalities in the computation of TEC associated with the use of greater data quality controls and source selection while real-time modeling environments must rely on autonomous controls and filtration techniques to prevent the production of erroneous model results. In this paper we present methods for processing TEC in real time, which utilize several procedures including the application of an ionospheric model to automatically perform quality control on the TEC output and the computational techniques used to address receiver multipath, faulty receiver observations, cycle-slips, segmented processing, and receiver calibrations. The resulting TEC measurements are provided with rigorous error estimates validated using the vertical TEC from the Jason satellite mission.

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 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.

A constrained LAMBDA method for GPS attitude determination

An improved method to obtain fixed integer ambiguity in GPS attitude determination is presented. Known conditions are utilized as constraints to acquire attitude information when the float solution and its variance–covariance matrix are not accurate enough. The searching ellipsoidal region is first expanded to compensate for errors caused by the inaccurate float solution. Then the constraints are used to shrink the region to a proper size, which maintains the true integer ambiguity. Experimental results demonstrate that this scheme gives a fast search time and a higher success rate in determining the fixed integer ambiguity than the unconstrained method. The accuracy of attitude angles is also improved.

A spatiotemporal analysis of aggregate labour force behaviour by sex and age across the European Union

This study investigates the causes of variation in age-specific male and female labour force participation rates using annual data from 154 regions across ten European Union member states for the period 1983–1997. Regional participation rates appear to be strongly correlated in time, weakly correlated in space and to parallel their national counterparts. An econometric model is designed consistent with these empirical findings. To control for potential endogeneity of the explanatory variables, we use an instrumental variables estimation scheme based on a matrix exponential spatial specification of the error terms. Many empirical studies of aggregate labour force behaviour have ignored population distribution effects, relying instead on the representative-agent paradigm. In order for representative-agent models to accurately describe aggregate behaviour, all marginal reactions of individuals to changes in aggregate variables must be identical. It turns out that this condition cannot apply to individuals across different sex/age groups.

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.

---

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.

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.

Modeling network autocorrelation within migration flows by eigenvector spatial filtering

Although the assumption of independence among interaction flows frequently is engaged in spatial interaction modeling, in many circumstances it leads to misspecified models and incorrect inferences. An informed approach is to explicitly incorporate an assumed relationship structure among the interaction flows, and to explicitly model the network autocorrelation. This paper illustrates such an approach in the context of U.S. interstate migration flows. Behavioral assumptions, similar to those of the intervening opportunities or the competing destinations concepts, exemplify how to specify network flows that are related to particular origin–destination combinations. The stepwise incorporation of eigenvectors, which are extracted from a network link matrix, captures the network autocorrelation in a Poisson regression model specification context. Spatial autocorrelation in Poisson regression is measured by the test statistic of Jacqmin-Gadda et al. (Stat Med 16(11):1283–1297, 1997). Results show that estimated regression parameters in the spatial filtering interaction model become more intuitively interpretable.

Kalman-filter-based GPS clock estimation for near real-time positioning

In this article, an algorithm for clock offset estimation of the GPS satellites is presented. The algorithm is based on a Kalman-filter and processes undifferenced code and carrier-phase measurements of a global tracking network. The clock offset and drift of the satellite clocks are estimated along with tracking station clock offsets, tropospheric zenith path delay and carrier-phase ambiguities. The article provides a brief overview of already existing near-real-time and real-time clock products. The filter algorithm and data processing scheme is presented. Finally, the accuracy of the orbit and clock product is assessed with a precise orbit determination of the MetOp satellite and compared to results gained with other real-time products.