共查询到20条相似文献,搜索用时 78 毫秒
1.
While there is a large literature on the form of epidemic waves in the time domain, models of their structure and shape in the spatial domain remain poorly developed. This paper concentrates on the changing spatial distribution of an epidemic wave over time and presents a simple method for identifying the leading and trailing edges of the spatial advance and retreat of such waves. Analysis of edge characteristics is used to (a) disaggregate waves into ‘swash’ and ‘backwash’ stages, (b) measure the phase transitions of areas from susceptible, S, through infective, I, to recovered, R, status (S → I → R) as dimensionless integrals and (c) estimate a spatial version of the basic reproduction number, R
0. The methods used are illustrated by application to measles waves in Iceland over a 60-year period from 1915 to 1974. Extensions of the methods for use with more complex waves are possible through modifying the threshold values used to define the start and end points of an event.
相似文献
Peter HaggettEmail: |
2.
This paper presents a continuous version of the model of distribution dynamics to analyse the transition dynamics and implied
long-run behaviour of the EU-27 NUTS-2 regions over the period 1995–2003. It departs from previous research in two respects:
first, by introducing kernel estimation and three-dimensional stacked conditional density plots as well as highest density regions
plots for the visualisation of the transition function, based on Hyndman et al. (J Comput Graph Stat 5(4):315–336, 1996), and second, by combining Getis’ spatial filtering view with kernel estimation to explicitly account for the spatial dimension of the
growth process. The results of the analysis indicate a very slow catching-up of the poorest regions with the richer ones,
a process of shifting away of a small group of very rich regions, and highlight the importance of geography in understanding
regional income distribution dynamics.
相似文献
Manfred M. FischerEmail: |
3.
Computer algebra solution of the GPS N-points problem 总被引:1,自引:0,他引:1
A computer algebra solution is applied here to develop and evaluate algorithms for solving the basic GPS navigation problem: finding a point position using four or more pseudoranges at one
epoch (the GPS N-points problem). Using Mathematica 5.2 software, the GPS N-points problem is solved numerically, symbolically, semi-symbolically, and with Gauss–Jacobi, on a work station. For the case of N > 4, two minimization approaches based on residuals and distance norms are evaluated for the direct numerical solution and their computational duration is compared. For N = 4, it is demonstrated that the symbolic computation is twice as fast as the iterative direct numerical method. For N = 6, the direct numerical solution is twice as fast as the semi-symbolic, with the residual minimization requiring less computation
time compared to the minimization of the distance norm. Gauss–Jacobi requires eight times more computation time than the direct
numerical solution. It does, however, have the advantage of diagnosing poor satellite geometry and outliers. Besides offering
a complete evaluation of these algorithms, we have developed Mathematica 5.2 code (a notebook file) for these algorithms (i.e., Sturmfel’s resultant, Dixon’s resultants, Groebner basis, reduced
Groebner basis and Gauss–Jacobi). These are accessible to any geodesist, geophysicist, or geoinformation scientist via the
GPS Toolbox () website or the Wolfram Information Center ().
相似文献
Erik W. GrafarendEmail: |
4.
L. E. Sjöberg 《Journal of Geodesy》1999,73(3):118-124
With access to dual-frequency pseudorange and phase Global Positioning System (GPS) data, the wide-lane ambiguity can easily
be fixed. Advantage is taken of this information in the linear combination of the above four observables for base ambiguity
estimation (i.e. of N
1 and N
2). Starting points for our analysis are the Best Linear Unbiased Estimators BLUE1 and BLUE2. BLUE1 is the best one (with minimum mean square error, MSE) if the ionosphere effect is negligible. If this is not the case, BLUE2 has the smallest variance, but not necessarily the least mean square error. Hence, both estimators may suffer from a non-optimal
treatment of the ionosphere bias. BLUE1 ignores possible ionosphere bias, while BLUE2 compensates for this bias in a less favourable way by eliminating it at the price of increased noise. As an alternative,
linear estimators are derived, which make a compromise between the ionosphere bias and the random observation errors. This
leads to the derivation of the Best Linear Estimator (BLE) and the Restricted Best Linear Estimator (RBLE) with minimum MSE.
The former is generally not very useful, while the RBLE is recommended for practical use. It is shown that the MSE of the
RBLE is limited by the variances of BLUE1 and BLUE2, i.e.
However, as is always the case with a BLE, it cannot be used strictly: some parameter (in this case the ionosphere bias) must
be approximately known.
Received: 25 June 1997 / Accepted: 16 November 1998 相似文献
5.
Alan T. Murray Timothy C. Matisziw Tony H. Grubesic 《Journal of Geographical Systems》2007,9(2):103-117
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.
相似文献
Alan T. Murray (Corresponding author)Email: |
Timothy C. MatisziwEmail: |
Tony H. GrubesicEmail: |
6.
Lars Sjöberg 《Journal of Geodesy》1979,53(3):227-230
The method of Bjerhammar is studied in the continuous case for a sphere. By varying the kernel function, different types of
unknowns (u*) are obtained at the internal sphere (the Bjerhammar sphere). It is shown that a necessary condition for the existence of
u* is that the degree variances (σ
n
2
) of the observations are of an order less than n−2. According to Kaula’s rule this condition is not satisfied for the earth’s gravity anomaly field (σ
n
2
=n−1) but well for the geopotential (σ
n
2
=n−3). 相似文献
7.
Engineering projects that require deformation monitoring frequently utilize geodetic sensors to measure displacements of target
points located in the deformation zone. In situations where control stations and targets are separated by a kilometer or more,
GPS can offer higher precision position updates at more frequent intervals than can normally be achieved using total station
technology. For large-scale deformation projects requiring the highest precision, it is therefore advisable to use a combination
of the two sensors. In response to the need for high precision, continuous GPS position updates in harsh deformation monitoring
environments, a software has been developed that employs triple-differenced carrier-phase measurements in a delayed-state
Kalman filter. Two data sets were analyzed to test the capabilities of the software. In the first test, a GPS antenna was
displaced using a translation stage to mimic slow deformation. In the second test, data collected at a large open pit mine
were processed. It was shown that the delayed-state Kalman filter developed could detect millimeter-level displacements of
a GPS antenna. The actual precision attained depends upon the amount of process noise infused at each epoch to accommodate
the antenna displacements. Higher process noise values result in quicker detection times, but at the same time increase the
noise in the solutions. A slow, 25 mm displacement was detected within 30 min of the full displacement with sigma values in
E, N and U of ±10 mm or better. The same displacement could also be detected in less than 5 h with sigma values in E, N and U of ±5 mm or better. The software works best for detecting long period deformations (e.g., 20 mm per day or less) for which
sigma values of 1–2 mm are attained in all three solution components. It was also shown that the triple-differenced carrier-phase
observation can be used to significantly reduce the effects of residual tropospheric delay that would normally plague double-differenced
observations in harsh GPS environments.
相似文献
Don KimEmail: |
8.
Gianfranco Piras Kieran P. Donaghy Giuseppe Arbia 《Journal of Geographical Systems》2007,9(4):311-344
This paper presents an innovative approach to the study of regional economic dynamics within a nonlinear continuous-time econometric
framework—a generalized specification of the Lotka–Volterra system of equations. This specification, which accounts for interdependent
behavior of three industrial sectors and spillover effects of activities in neighboring regions, is employed in an analysis
of five Italian regions between 1980 and 2003. For these regions, we report estimation results, characterize the varying systems
dynamics, analyze the models’ local and global stability properties, and determine via sensitivity analyses which structural
features appear to exert the greatest influence on these properties.
相似文献
Kieran P. DonaghyEmail: |
9.
Hierarchical spatial interaction among the Italian regions: a nonlinear relative dynamics approach 总被引:2,自引:2,他引:0
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.
相似文献
Geoffrey J. D. HewingsEmail: |
10.
The resolution of a nonlinear parametric adjustment model is addressed through an isomorphic geometrical setup with tensor
structure and notation, represented by a u-dimensional “model surface” embedded in a flat n-dimensional “observational space”.
Then observations correspond to the observational-space coordinates of the pointQ, theu initial parameters correspond to the model-surface coordinates of the “initial” pointP, and theu adjusted parameters correspond to the model-surface coordinates of the “least-squares” point
. The least-squares criterion results in a minimum-distance property implying that the vector
Q must be orthogonal to the model surface. The geometrical setup leads to the solution of modified normal equations, characterized
by a positive-definite matrix. The latter contains second-order and, optionally, thirdorder partial derivatives of the observables
with respect to the parameters. This approach significantly shortens the convergence process as compared to the standard (linearized)
method. 相似文献
11.
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.
相似文献
Daniel A. GriffithEmail: |
12.
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.
相似文献
Guochang XuEmail: |
13.
A data archive of GPS navigation messages 总被引:2,自引:1,他引:1
Since 18 June 2007 navigation data messages transmitted by the GPS constellation are recorded by five receivers within GeoForschungsZentrum’s
global groundstation network. We describe the recording, processing, validation, analysis and archiving of the navigation
data. During the 197 days between 18 June 2007 and 31 December 2007 a total of 125,723,666 subframes were collected. By taking
into consideration that the same data set frequently is observed by two or more receivers concurrently, 65,153,955 unique
subframes could be extracted from the observations. With an estimated 88,099,200 subframes transmitted by the constellation
during this time period a data yield of about 74% was achieved. Simulation studies suggest that with two additional GPS receivers,
which are scheduled for addition to the network in 2008, about 95% of the transmitted subframes will be retrieved. The message
data archive is open to the scientific community for non-commercial purposes and may be accessed through GFZ’s Information
System and Data Center ().
相似文献
G. BeyerleEmail: |
14.
E.W. Grafarend 《Journal of Geodesy》2005,78(10):594-615
Harmonic maps are generated as a certain class of optimal map projections. For instance, if the distortion energy over a meridian strip of the International Reference Ellipsoid is minimized, we are led to the Laplace–Beltrami vector-valued partial differential equation. Harmonic functions x(L,B), y(L,B) given as functions of ellipsoidal surface parameters of Gauss ellipsoidal longitude L and Gauss ellipsoidal latitude B, as well as x(,q), y(,q) given as functions of relative isometric longitude =L–L0 and relative isometric latitude q=Q–Q0 gauged to a vector-valued boundary condition of special symmetry are constructed. The easting and northing {x(b,),y(b,)} of the new harmonic map is then given. Distortion energy analysis of the new harmonic map is presented, as well as case studies for (1) B[–40°,+40°], L[–31°,+49°], B0= ±30°, L0=9° and (2) B[46°,56°], L{[4.5°, 7.5°]; [7.5°, 10.5°]; [10.5°,13.5°]; [13.5°,16.5°]}, B0= 51°, L0 {6°,9°,12°,15°}. 相似文献
15.
Burkhard Schaffrin 《Journal of Geodesy》1987,61(3):276-280
The Bayesian estimates
b of the standard deviation σ in a linear model—as needed for the evaluation of reliability—is well known to be proportional
to the square root of the Bayesian estimate (s
2)
b
of the variance component σ2 by a proportionality factor
involving the ratio of Gamma functions. However, in analogy to the case of the respective unbiased estimates, the troublesome
exact computation ofa
b may be avoided by a simple approximation which turns out to be good enough for most applications even if the degree of freedom
ν is rather small.
Paper presented to the Int. Conf. on “Practical Bayesian Statistics”, Cambridge (U.K.), 8.–11. July 1986. 相似文献
16.
A constrained LAMBDA method for GPS attitude determination 总被引:4,自引:0,他引:4
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.
相似文献
Bo WangEmail: |
17.
Urho A. Rauhala 《Journal of Geodesy》1979,53(4):317-342
Array algebra forms the general base of fast transforms and multilinear algebra making rigorous solutions of a large number
(millions) of parameters computationally feasible. Loop inverses are operators solving the problem of general matrix inverses.
Their derivation starts from the inconsistent linear equations
by a parameter exchangeX→L
0, where
X is a set of unknown observables,A
0 forming a basis of the so called “problem space”. The resulting full rank design matrix of parameters L0 and its ℓ-inverse reveal properties speeding the computational least squares solution
expressed in observed values
. The loop inverses are found by the back substitution expressing ∧X in terms ofL through
. Ifp=rank (A) ≤n, this chain operator creates the pseudoinverseA
+. The idea of loop inverses and array algebra started in the late60's from the further specialized case,p=n=rank (A), where the loop inverse A
0
−1
(AA
0
−1
)ℓ reduces into the ℓ-inverse Aℓ=(ATA)−1AT. The physical interpretation of the design matrixA A
0
−1
as an interpolator, associated with the parametersL
0, and the consideration of its multidimensional version has resulted in extended rules of matrix and tensor calculus and mathematical
statistics called array algebra. 相似文献
18.
Modeling spatial autocorrelation in spatial interaction data: empirical evidence from 2002 Germany journey-to-work flows 总被引:3,自引:3,他引:0
Daniel A. Griffith 《Journal of Geographical Systems》2009,11(2):117-140
Since before the inception of work by Okabe, the intermingling of spatial autocorrelation (i.e., local distance and configuration)
and distance decay (i.e., global distance) effects has been suspected in spatial interaction data. This convolution was first
treated conceptually because technology and methodology did not exist at the time to easily or fully address spatial autocorrelation
effects within spatial interaction model specifications. Today, however, sufficient computer power coupled with eigenfunction-based
spatial filtering offers a means for accommodating spatial autocorrelation effects within a spatial interaction model for
modest-sized problems. In keeping with Okabe’s more recent efforts to dissemination spatial analysis tools, this paper summarizes
how to implement the methodology utilized to analyze a particular empirical flows dataset.
相似文献
Daniel A. GriffithEmail: |
19.
The regularized solution of the external sphericalStokes boundary value problem as being used for computations of geoid undulations and deflections of the vertical is based upon theGreen functions S
1(0, 0, , ) ofBox 0.1 (R = R
0) andV
1(0, 0, , ) ofBox 0.2 (R = R
0) which depend on theevaluation point {0, 0} S
R0
2
and thesampling point {, } S
R0
2
ofgravity anomalies
(, ) with respect to a normal gravitational field of typegm/R (free air anomaly). If the evaluation point is taken as the meta-north pole of theStokes reference sphere S
R0
2
, theStokes function, and theVening-Meinesz function, respectively, takes the formS() ofBox 0.1, andV
2() ofBox 0.2, respectively, as soon as we introduce {meta-longitude (azimuth), meta-colatitude (spherical distance)}, namely {A, } ofBox 0.5. In order to deriveStokes functions andVening-Meinesz functions as well as their integrals, theStokes andVening-Meinesz functionals, in aconvolutive form we map the sampling point {, } onto the tangent plane T0S
R0
2
at {0, 0} by means ofoblique map projections of type(i) equidistant (Riemann polar/normal coordinates),(ii) conformal and(iii) equiareal.Box 2.1.–2.4. andBox 3.1.– 3.4. are collections of the rigorously transformedconvolutive Stokes functions andStokes integrals andconvolutive Vening-Meinesz functions andVening-Meinesz integrals. The graphs of the correspondingStokes functions S
2(),S
3(r),,S
6(r) as well as the correspondingStokes-Helmert functions H
2(),H
3(r),,H
6(r) are given byFigure 4.1–4.5. In contrast, the graphs ofFigure 4.6–4.10 illustrate the correspondingVening-Meinesz functions V
2(),V
3(r),,V
6(r) as well as the correspondingVening-Meinesz-Helmert functions Q
2(),Q
3(r),,Q
6(r). The difference between theStokes functions / Vening-Meinesz functions andtheir first term (only used in the Flat Fourier Transforms of type FAST and FASZ), namelyS
2() – (sin /2)–1,S
3(r) – (sinr/2R
0)–1,,S
6(r) – 2R
0/r andV
2() + (cos /2)/2(sin2 /2),V
3(r) + (cosr/2R
0)/2(sin2
r/2R
0),,
illustrate the systematic errors in theflat Stokes function 2/ or flatVening-Meinesz function –2/2. The newly derivedStokes functions S
3(r),,S
6(r) ofBox 2.1–2.3, ofStokes integrals ofBox 2.4, as well asVening-Meinesz functionsV
3(r),,V
6(r) ofBox 3.1–3.3, ofVening-Meinesz integrals ofBox 3.4 — all of convolutive type — pave the way for the rigorousFast Fourier Transform and the rigorousWavelet Transform of theStokes integral / theVening-Meinesz integral of type equidistant, conformal and equiareal. 相似文献
20.
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.
相似文献
Yun-Wen Huang (Corresponding author)Email: |