A platform for research: civil engineering, architecture and urbanism
A platform for research: civil engineering, architecture and urbanism
On the inversion of geodetic integrals defined over the sphere using 1-D FFT
Abstract An iterative method is presented which performs inversion of integrals defined over the sphere. The method is based on one-dimensional fast Fourier transform (1-D FFT) inversion and is implemented with the projected Landweber technique, which is used to solve constrained least-squares problems reducing the associated 1-D cyclic-convolution error. The results obtained are as precise as the direct matrix inversion approach, but with better computational efficiency. A case study uses the inversion of Hotine’s integral to obtain gravity disturbances from geoid undulations. Numerical convergence is also analyzed and comparisons with respect to the direct matrix inversion method using conjugate gradient (CG) iteration are presented. Like the CG method, the number of iterations needed to get the optimum (i.e., small) error decreases as the measurement noise increases. Nevertheless, for discrete data given over a whole parallel band, the method can be applied directly without implementing the projected Landweber method, since no cyclic convolution error exists.
On the inversion of geodetic integrals defined over the sphere using 1-D FFT
Abstract An iterative method is presented which performs inversion of integrals defined over the sphere. The method is based on one-dimensional fast Fourier transform (1-D FFT) inversion and is implemented with the projected Landweber technique, which is used to solve constrained least-squares problems reducing the associated 1-D cyclic-convolution error. The results obtained are as precise as the direct matrix inversion approach, but with better computational efficiency. A case study uses the inversion of Hotine’s integral to obtain gravity disturbances from geoid undulations. Numerical convergence is also analyzed and comparisons with respect to the direct matrix inversion method using conjugate gradient (CG) iteration are presented. Like the CG method, the number of iterations needed to get the optimum (i.e., small) error decreases as the measurement noise increases. Nevertheless, for discrete data given over a whole parallel band, the method can be applied directly without implementing the projected Landweber method, since no cyclic convolution error exists.
On the inversion of geodetic integrals defined over the sphere using 1-D FFT
García, R. V. (author) / Alejo, C. A. (author)
Journal of Geodesy ; 79
2005
Article (Journal)
Electronic Resource
English
On the inversion of geodetic integrals defined over the sphere using 1-D FFT
| Online Contents | 2005
Thin-plate spline quadrature of geodetic integrals
| Online Contents | 1994
Fast space-domain evaluation of geodetic surface integrals
| Online Contents | 1997
On the accurate numerical evaluation of geodetic convolution integrals
| Online Contents | 2011
On the accurate numerical evaluation of geodetic convolution integrals
| Online Contents | 2011