首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
Wavelet evaluation of the Stokes and Vening Meinesz integrals   总被引：1，自引：0，他引：1  
The wavelet transform is a powerful tool in evaluating some singular geodetic integrals. Due to its localization properties in both the time (space) and frequency (scale) domains, and because the kernels of some geodetic integrals have singular points and decay smoothly and quickly away from the singularities, many wavelet transform coefficients of the kernels become zeros or negligible, and only a small number of wavelet transform coefficients are significant. It is thus possible to significantly compress the kernels of these integrals on a wavelet basis by neglecting the zero coefficients and the small coefficients below a certain threshold. Therefore, wavelets provide a convenient way of efficiently evaluating these integrals in terms of fast computation and savings of computer memory. A modified algorithm for the wavelet evaluation of Stokes' integral is presented. The same modified algorithm is applied to the evaluation of the Vening Meinesz integral, whose kernel has a stronger singularity than does Stokes' kernel. Numerical examples illustrate the efficiency and accuracy of the wavelet methods.Acknowledgments.The author express their sincere thanks to Dr. Salamonwicz for providing his PhD thesis. E-mail correspondence between the authors and Dr. Barthelmes and Dr. Benciolini contributed to the work. R. Benciolini and the other two anonymous reviewers are thanked for their constructive comments. Support for this research was provided by research grants to Dr. Sideris from the Natural Sciences and Engineering Reserch Council of Canada (NSERC) and the Geomatics for Informed Decisions (GEOIDE) Network of Centres of Excellence. The MATLAB Wavelet Toolbox package was used as the platform to develop the software in this project.  相似文献   

2.
The previous work of Xu on discrete nonlinear filtering is extended to continuous systems. The new results are summarized as follows: (1) a second-order unbiased prediction of the true state governed by a vector stochastic differential equation is worked out; (2) a set of coupled differential equations for a new truncated second-order nonlinear filter and its variance–covariance matrix are derived from the frequentist point of view. The new filter is proved to be unbiased to the second-order approximation; and, most importantly, (3) comparison of the new filtering and accuracy results with the literature on nonlinear filtering has indicated that more than 40 years of nonlinear filtering of continuous systems may have foundational problems.Acknowledgments.This work is supported by a Grant-in-Aid for Scientific Research (C13640422). The author thanks Prof J.A.R. Blais, Prof A. Dermanis and Prof B. Schaffrin for their constructive comments.  相似文献   

3.
An algorithm for the determination of the spherical harmonic coefficients of the terrestrial gravitational field representation from the analysis of a kinematic orbit solution of a low earth orbiting GPS-tracked satellite is presented and examined. A gain in accuracy is expected since the kinematic orbit of a LEO satellite can nowadays be determined with very high precision, in the range of a few centimeters. In particular, advantage is taken of Newton's Law of Motion, which balances the acceleration vector with respect to an inertial frame of reference (IRF) and the gradient of the gravitational potential. By means of triple differences, and in particular higher-order differences (seven-point scheme, nine-point scheme), based upon Newton's interpolation formula, the local acceleration vector is estimated from relative GPS position time series. The gradient of the gravitational potential is conventionally given in a body-fixed frame of reference (BRF) where it is nearly time independent or stationary. Accordingly, the gradient of the gravitational potential has to be transformed from spherical BRF to Cartesian IRF. Such a transformation is possible by differentiating the gravitational potential, given as a spherical harmonics series expansion, with respect to Cartesian coordinates by means of the chain rule, and expressing zero- and first-order Ferrer's associated Legendre functions in terms of Cartesian coordinates. Subsequently, the BRF Cartesian coordinates are transformed into IRF Cartesian coordinates by means of the polar motion matrix, the precession–nutation matrices and the Greenwich sidereal time angle (GAST). In such a way a spherical harmonic representation of the terrestrial gravitational field intensity with respect to an IRF is achieved. Numerical tests of a resulting Gauss–Markov model document not only the quality and the high resolution of such a space gravity spectroscopy, but also the problems resulting from noise amplification in the acceleration determination process.  相似文献   

4.
The three-dimensional (3-D) resection problem is usually solved by first obtaining the distances connecting the unknown point P{X,Y,Z} to the known points Pi{Xi,Yi,Zi}i=1,2,3 through the solution of the three nonlinear Grunert equations and then using the obtained distances to determine the position {X,Y,Z} and the 3-D orientation parameters {,, }. Starting from the work of the German J. A. Grunert (1841), the Grunert equations have been solved in several substitutional steps and the desire as evidenced by several publications has been to reduce these number of steps. Similarly, the 3-D ranging step for position determination which follows the distance determination step involves the solution of three nonlinear ranging (`Bogenschnitt') equations solved in several substitution steps. It is illustrated how the algebraic technique of Groebner basis solves explicitly the nonlinear Grunert distance equations and the nonlinear 3-D ranging (`Bogenschnitt') equations in a single step once the equations have been converted into algebraic (polynomial) form. In particular, the algebraic tool of the Groebner basis provides symbolic solutions to the problem of 3-D resection. The various forward and backward substitution steps inherent in the classical closed-form solutions of the problem are avoided. Similar to the Gauss elimination technique in linear systems of equations, the Groebner basis eliminates several variables in a multivariate system of nonlinear equations in such a manner that the end product normally consists of a univariate polynomial whose roots can be determined by existing programs e.g. by using the roots command in Matlab.Acknowledgments.The first author wishes to acknowledge the support of JSPS (Japan Society of Promotion of Science) for the financial support that enabled the completion of the write-up of the paper at Kyoto University, Japan. The author is further grateful for the warm welcome and the good working atmosphere provided by his hosts Professors S. Takemoto and Y. Fukuda of the Department of Geophysics, Graduate School of Science, Kyoto University, Japan.  相似文献   

5.
  相似文献   

6.
Mapping and monitoring impervious surface dynamic change in a complex urban-rural frontier with medium or coarse spatial resolution images is a challenge due to the mixed pixel problem and the spectral confusion between impervious surfaces and other non-vegetation land covers. This research selected Lucas do Rio Verde County in Mato Grosso State, Brazil as a case study to improve impervious surface estimation performance by the integrated use of Landsat and QuickBird images and to monitor impervious surface change by analyzing the normalized multitemporal Landsat-derived fractional impervious surfaces. This research demonstrates the importance of two-step calibrations. The first step is to calibrate the Landsat-derived fraction impervious surface values through the established regression model based on the QuickBird-derived impervious surface image in 2008. The second step is to conduct the normalization between the calibrated 2008 impervious surface image with other dates of impervious surface images. This research indicates that the per-pixel based method overestimates the impervious surface area in the urban-rural frontier by 50%-60%. In order to accurately estimate impervious surface area, it is necessary to map the fractional impervious surface image and further calibrate the estimates with high spatial resolution images. Also normalization of the multitemporal fractional impervious surface images is needed to reduce the impacts from different environmental conditions, in order to effectively detect the impervious surface dynamic change in a complex urban-rural frontier. The procedure developed in this paper for mapping and monitoring impervious surface area is especially valuable in urban-rural frontiers where multitemporal Landsat images are difficult to be used for accurately extracting impervious surface features based on traditional per-pixel based classification methods as they cannot effectively handle the mixed pixel problem.  相似文献   

7.
The role of compatibility conditions for altimetry–gravimetry boundary value problems (AGBVPs) is quite complicated by way of: guaranteeing the existence of the solution; smoothing the data (boundary conditions) at the boundary; and providing a higher level of regularization and smoothness to the solution. Pseudo-differential operators (PDOs) can be applied not only for the reformulation of AGBVPs in a simple way, but also for imposing compatibility conditions at the coastline. It is shown here that spherical PDOs can be combined with the theory of spherical harmonics and spherical wavelets to obtain a numerical solution for AGBVPs with compatibility conditions along the coastline. The role of compatibility conditions, is summarized and emphasized and compatibility conditions are given in an explicit form based on the combined representation of functionals of the disturbing potential in terms of spherical harmonics, spherical PDOs and wavelets. Numerical aspects of the application of the compatibility conditions to a case-study area in Canada are discussed.  相似文献   

8.
 采用Landsat TM数据分析了阿克苏河—塔里木河断面水质污染状况,通过波段的DN值和常规监测数据建立能反映水质状况的 污染物监测模型。结果发现,将2000年常规监测数据代入模型后,与遥感数据的结果基本吻合| 重建阿克苏河—塔里木河的连续水体 污染变化曲线,得出污染物浓度随着远离上游而增加。  相似文献   

9.
How to handle topography in practical geoid determination: three examples   总被引：3，自引：1，他引：2  
 Three different methods of handling topography in geoid determination were investigated. The first two methods employ the residual terrain model (RTM) remove–restore technique, yielding the quasigeoid, whereas the third method uses the classical Helmert condensation method, yielding the geoid. All three methods were used with the geopotential model Earth Gravity Model (1996) (EGM96) as a reference, and the results were compared to precise global positioning system (GPS) levelling networks in Scandinavia. An investigation of the Helmert method, focusing on the different types of indirect effects and their effects on the geoid, was also carried out. The three different methods used produce almost identical results at the 5-cm level, when compared to the GPS levelling networks. However, small systematic differences existed. Received: 18 March 1999 / Accepted: 21 March 2000  相似文献   

10.
介绍了Pleiades卫星影像数据及基于Erdas 2013的Pleiades卫星影像应用的可行性和区域网平差技术,为该卫星影像应用于相关测绘工程提供了技术依据。  相似文献   

11.
 Several pre-analysis measures which help to expose the behavior of L 1 -norm minimization solutions are described. The pre-analysis measures are primarily based on familiar elements of the linear programming solution to L 1-norm minimization, such as slack variables and the reduced-cost vector. By examining certain elements of the linear programming solution in a probabilistic light, it is possible to derive the cumulative distribution function (CDF) associated with univariate L 1-norm residuals. Unlike traditional least squares (LS) residual CDFs, it is found that L 1-norm residual CDFs fail to follow the normal distribution in general, and instead are characterized by both discrete and continuous (i.e. piecewise) segments. It is also found that an L 1 equivalent to LS redundancy numbers exists and that these L 1 equivalents are a byproduct of the univariate L 1 univariate residual CDF. Probing deeper into the linear programming solution, it is found that certain combinations of observations which are capable of tolerating large-magnitude gross errors can be predicted by comprehensively tabulating the signs of slack variables associated with the L 1 residuals. The developed techniques are illustrated on a two-dimensional trilateration network. Received: 6 July 2001 / Accepted: 21 February 2002  相似文献   

12.
 Two iterative vector methods for computing geodetic coordinates (φ, h) from rectangular coordinates (x, y, z) are presented. The methods are conceptually simple, work without modification at any latitude and are easy to program. Geodetic latitude and height can be calculated to acceptable precision in one iteration over the height range from −106 to +109 m. Received: 13 December 2000 / Accepted: 13 July 2001  相似文献   

13.
关于我国基础测绘发展战略问题研究   总被引：5，自引：0，他引：5  
对我国基础测绘发展现状和趋势进行了分析,并就我国基础测绘发展的战略目标和实现这些目标需要开展的主要任务进行了深入探讨。  相似文献   

14.
 The use of GPS for height control in an area with existing levelling data requires the determination of a local geoid and the bias between the local levelling datum and the one implicitly defined when computing the local geoid. If only scarse gravity data are available, the heights of new data may be collected rapidly by determining the ellipsoidal height by GPS and not using orthometric heights. Hence the geoid determination has to be based on gravity disturbances contingently combined with gravity anomalies. Furthermore, existing GPS/levelling data may also be used in the geoid determination if a suitable general gravity field modelling method (such as least-squares collocation, LSC) is applied. A comparison has been made in the Aswan Dam area between geoids determined using fast Fourier transform (FFT) with gravity disturbances exclusively and LSC using only the gravity disturbances and the disturbances combined with GPS/levelling data. The EGM96 spherical harmonic model was in all cases used in a remove–restore mode. A total of 198 gravity disturbances spaced approximately 3 km apart were used, as well as 35 GPS/levelling points in the vicinity and on the Aswan Dam. No data on the Nasser Lake were available. This gave difficulties when using FFT, which requires the use of gridded data. When using exclusively the gravity disturbances, the agreement between the GPS/levelling data were 0.71 ± 0.17 m for FFT and 0.63 ± 0.15 for LSC. When combining gravity disturbances and GPS/levelling, the LSC error estimate was ±0.10 m. In the latter case two bias parameters had to be introduced to account for a possible levelling datum difference between the levelling on the dam and that on the adjacent roads. Received: 14 August 2000 / Accepted: 28 February 2001  相似文献   

15.
 Considering a GPS satellite and two terrestrial stations, two types of equations are derived relating the heights of the two stations to the measured data (frequency ratio or clock rate differences) and the coordinates and velocity components of all three participating objects. The potential possibilities of using such relations for the determination of heights (in terms of geopotential numbers or orthometric heights) are discussed. Received: 6 December 2000 / Accepted: 9 July 2001  相似文献   

16.
Deconvolution with wavelets and vaguelettes   总被引：2，自引：0，他引：2  
The use of wavelets for the solution of convolution equations is studied as a possible alternative to the well-established Fast Fourier Transform (FFT) technique. Two possible solution strategies are investigated: (1) The use of wavelets for the representation of both the given data and the unknown solution. This leads to an algorithm with good de-noising and data-compression properties. In terms of computational efficiency this algorithm is inferior to FFT. (2) The use of wavelets for the representation of the unknown solution and of so-called vaguelettes for the representations of the given data. This leads to an algorithm which is even faster than FFT. Received: 14 October 1998 / Accepted: 30 November 1999  相似文献   

17.
New solutions for the geodetic coordinate transformation   总被引：5，自引：2，他引：5  
 The Cartesian-to-geodetic-coordinate transformation is approached from a new perspective. Existence and uniqueness of geodetic representation are presented, along with a clear geometric picture of the problem and the role of the ellipse evolute. A new solution is found with a Newton-method iteration in the reduced latitude; this solution is proved to work for all points in space. Care is given to error propagation when calculating the geodetic latitude and height. Received: 9 August 2001 / Accepted: 27 March 2002 Acknowledgments. The author would like to thank the Clifford W.␣Tompson scholarship fund, Dr. Brian DeFacio, the University of Missouri College of Arts &Sciences, and the United States Air Force. He also thanks a reviewer for suggesting and providing a prototype MATLAB code. A MATLAB program for the iterative sequence is presented at the end of the paper (Appendix A).  相似文献   

18.
 A new method called Trop_NetAdjust is described to predict in real time the residual tropospheric delays on the GPS carrier phase observables using the redundant measurements from a network of GPS reference stations. This method can not only enhance the effectiveness and reliability of real-time kinematic users within the network, but also provide a valid approach to tropospheric parameter variation forecasting. Trop_NetAdjust is theoretically based upon LS prediction criteria and enables the prediction of residual tropospheric delays remaining after a standard model has been applied to the raw GPS measurements. Two cases are analyzed, namely a first case when the delay is required for an existing satellite at a new point within the network and a second case when the delay is required for a new satellite. Field tests were conducted using data collected in a network of 11 reference stations covering a 400×600 km region in southern Norway. The results were analyzed in the measurement domain (ionospheric-free double-difference residuals) and showed improvements of 20 to 65% RMS errors using Trop_NetAdjust. The estimates of the Trop_NetAdjust prediction accuracy were also obtained using the covariance analysis method. The agreement was consistently better than 30% when compared with data from a real network. Received: 28 February 2000 / Accepted: 9 January 2001  相似文献   

19.
  The Western Alps are among the best studied collisional belts with both detailed structural mapping and also crustal geophysical investigations such as the ECORS and EGT seismic profile. By contrast, the present-day kinematics of the belt is still largely unknown due to small relative motions and the insufficient accuracy of the triangulation data. As a consequence, several tectonic problems still remain to be solved, such as the amount of N–S convergence in the Occidental Alps, the repartition of the deformation between the Alpine tectonic units, and the relation between deformation and rotation across the Alpine arc. In order to address these problems, the GPS ALPES group, made up of French, Swiss and Italian research organizations, has achieved the first large-scale GPS surveys of the Western Alps. More than 60 sites were surveyed in 1993 and 1998 with a minimum observation of 3 days at each site. GPS data processing has been done by three independent teams using different software. The different solutions have horizontal repeatabilities (N–E) of 4–7 mm in 1993 and 2–3 mm in 1998 and compare at the 3–5-mm level in position and 2-mm/yr level in velocity. A comparison of 1993 and 1998 coordinates shows that residual velocities of the GPS marks are generally smaller than 2 mm/yr, precluding a detailed tectonic interpretation of the differential motions. However, these data seem to suggest that the N–S compression of the Western Alps is quite mild (less than 2 mm/yr) compared to the global convergence between the African and Eurasian plate (6 mm/yr). This implies that the shortening must be accomodated elsewhere by the deformation of the Maghrebids and/or by rotations of Mediterranean microplates. Also, E–W velocity components analysis supports the idea that E–W extension exists, as already suggested by recent structural and seismotectonic data interpretation. Received: 27 November 2000 / Accepted: 17 September 2001  相似文献   

20.
GPS-assisted GLONASS orbit determination   总被引：1，自引：0，他引：1  
 Using 1 week of data from a network of GPS/GLONASS dual-tracking receivers, 15-cm accurate GLONASS orbit determination is demonstrated with an approach that combines GPS and GLONASS data. GPS data are used to define the reference frame, synchronize receiver clocks and determine troposphere delay for the GLONASS tracking network. GLONASS tracking data are then processed separately, with the GPS-defined parameters held fixed, to determine the GLONASS orbit. The quality of the GLONASS orbit determination is currently limited by the size and distribution of the tracking network, and by the unavailability of a sufficiently refined solar pressure model. Temporal variations in the differential clock bias of the dual-tracking receivers are found to have secondary impact on the orbit determination accuracy. Received: 5 January 2000 / Accepted: 15 February 2001  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号