Efficient spatial-spectral computation of local planar gravimetric terrain corrections from high-resolution digital elevation models
Computation of gravimetric terrain corrections (TCs) is a numerical challenge, especially when using very high-resolution (say, ∼30 m or less) digital elevation models (DEMs). TC computations can use spatial or/and spectral techniques: Spatial domain methods are more exact but can be very time-consu...
| Main Authors: | , , , |
|---|---|
| Format: | Journal Article |
| Language: | English |
| Published: |
OXFORD UNIV PRESS
2020
|
| Subjects: | |
| Online Access: | http://hdl.handle.net/20.500.11937/82547 |
| _version_ | 1848764517819875328 |
|---|---|
| author | Goyal, Ropesh Featherstone, Will Tsoulis, D. Dikshit, O. |
| author_facet | Goyal, Ropesh Featherstone, Will Tsoulis, D. Dikshit, O. |
| author_sort | Goyal, Ropesh |
| building | Curtin Institutional Repository |
| collection | Online Access |
| description | Computation of gravimetric terrain corrections (TCs) is a numerical challenge, especially when using very high-resolution (say, ∼30 m or less) digital elevation models (DEMs). TC computations can use spatial or/and spectral techniques: Spatial domain methods are more exact but can be very time-consuming; the discrete/fast Fourier transform (D/FFT) implementation of a binomial expansion is efficient, but fails to achieve a convergent solution for terrain slopes >45°. We show that this condition must be satisfied for each and every computation-roving point pair in the whole integration domain, not just at or near the computation points. A combination of spatial and spectral methods has been advocated by some through dividing the integration domain into inner and outer zones, where the TC is computed from the superposition of analytical mass-prism integration and the D/FFT. However, there remain two unresolved issues with this combined approach: (1) deciding upon a radius that best separates the inner and outer zones and (2) analytical mass-prism integration in the inner zone remains time-consuming, particularly for high-resolution DEMs. This paper provides a solution by proposing: (1) three methods to define the radius separating the inner and outer zones and (2) a numerical solution for near-zone TC computations based on the trapezoidal and Simpson's rules that is sufficiently accurate w.r.t. the exact analytical solution, but which can reduce the computation time by almost 50 per cent. |
| first_indexed | 2025-11-14T11:20:37Z |
| format | Journal Article |
| id | curtin-20.500.11937-82547 |
| institution | Curtin University Malaysia |
| institution_category | Local University |
| language | English |
| last_indexed | 2025-11-14T11:20:37Z |
| publishDate | 2020 |
| publisher | OXFORD UNIV PRESS |
| recordtype | eprints |
| repository_type | Digital Repository |
| spelling | curtin-20.500.11937-825472021-03-05T04:15:05Z Efficient spatial-spectral computation of local planar gravimetric terrain corrections from high-resolution digital elevation models Goyal, Ropesh Featherstone, Will Tsoulis, D. Dikshit, O. Science & Technology Physical Sciences Geochemistry & Geophysics Gravity anomalies and Earth structure Fourier analysis Numerical approximations and analysis Numerical solutions GRAVITY GRADIOMETRY DENSITY PRISM FIELD FFT Computation of gravimetric terrain corrections (TCs) is a numerical challenge, especially when using very high-resolution (say, ∼30 m or less) digital elevation models (DEMs). TC computations can use spatial or/and spectral techniques: Spatial domain methods are more exact but can be very time-consuming; the discrete/fast Fourier transform (D/FFT) implementation of a binomial expansion is efficient, but fails to achieve a convergent solution for terrain slopes >45°. We show that this condition must be satisfied for each and every computation-roving point pair in the whole integration domain, not just at or near the computation points. A combination of spatial and spectral methods has been advocated by some through dividing the integration domain into inner and outer zones, where the TC is computed from the superposition of analytical mass-prism integration and the D/FFT. However, there remain two unresolved issues with this combined approach: (1) deciding upon a radius that best separates the inner and outer zones and (2) analytical mass-prism integration in the inner zone remains time-consuming, particularly for high-resolution DEMs. This paper provides a solution by proposing: (1) three methods to define the radius separating the inner and outer zones and (2) a numerical solution for near-zone TC computations based on the trapezoidal and Simpson's rules that is sufficiently accurate w.r.t. the exact analytical solution, but which can reduce the computation time by almost 50 per cent. 2020 Journal Article http://hdl.handle.net/20.500.11937/82547 10.1093/gji/ggaa107 English OXFORD UNIV PRESS fulltext |
| spellingShingle | Science & Technology Physical Sciences Geochemistry & Geophysics Gravity anomalies and Earth structure Fourier analysis Numerical approximations and analysis Numerical solutions GRAVITY GRADIOMETRY DENSITY PRISM FIELD FFT Goyal, Ropesh Featherstone, Will Tsoulis, D. Dikshit, O. Efficient spatial-spectral computation of local planar gravimetric terrain corrections from high-resolution digital elevation models |
| title | Efficient spatial-spectral computation of local planar gravimetric terrain corrections from high-resolution digital elevation models |
| title_full | Efficient spatial-spectral computation of local planar gravimetric terrain corrections from high-resolution digital elevation models |
| title_fullStr | Efficient spatial-spectral computation of local planar gravimetric terrain corrections from high-resolution digital elevation models |
| title_full_unstemmed | Efficient spatial-spectral computation of local planar gravimetric terrain corrections from high-resolution digital elevation models |
| title_short | Efficient spatial-spectral computation of local planar gravimetric terrain corrections from high-resolution digital elevation models |
| title_sort | efficient spatial-spectral computation of local planar gravimetric terrain corrections from high-resolution digital elevation models |
| topic | Science & Technology Physical Sciences Geochemistry & Geophysics Gravity anomalies and Earth structure Fourier analysis Numerical approximations and analysis Numerical solutions GRAVITY GRADIOMETRY DENSITY PRISM FIELD FFT |
| url | http://hdl.handle.net/20.500.11937/82547 |