Stable reduction to the pole at the magnetic equator

Geophysics ◽  
2001 ◽  
Vol 66 (2) ◽  
pp. 571-578 ◽  
Author(s):  
Yaoguo Li ◽  
Douglas W. Oldenburg

The solution of reduction to the pole (RTP) of magnetic data in the wavenumber domain faces a long standing difficulty of instability when the observed data are acquired at low magnetic latitudes or at the equator. We develop a solution to this problem that allows stable reconstruction of the RTP field with a high fidelity even at the magnetic equator. The solution is obtained by inverting the Fourier transform of the observed magnetic data in the wavenumber domain with explicit regularization. The degree of regularization is chosen according to the estimated error level in the data. The Fourier transform of the RTP field is constructed as a model that is maximally smooth and, at the same time, has a power‐spectral decay common to all fields produced by the same source. The applied regularization alleviates the singularity associated with the wavenumber‐domain RTP operator, and the imposed power spectral decay ensures that the constructed RTP field has the correct spectral content. As a result, the algorithm can perform the reduction to the pole stably at any magnetic latitude, and the constructed RTP field yields a good representation of the true field at the pole even when the reduction is carried out at the equator.

GEOMATIKA ◽  
2018 ◽  
Vol 23 (2) ◽  
pp. 65 ◽  
Author(s):  
Mila Apriani ◽  
Admiral Musa Julius ◽  
Mahmud Yusuf ◽  
Damianus Tri Heryanto ◽  
Agus Marsono

<p align="center"><strong>ABSTRAK</strong></p><p> </p><p>Penelitian dengan analisis <em>power spectral</em> data anomali gayaberat telah banyak dilakukan untuk estimasi ketebalan sedimen. Dalam studi ini penulis melakukan analisis spektral data anomali gayaberat wilayah DKI Jakarta untuk mengetahui kedalaman sumber anomali yang bersesuaian dengan ketebalan sedimen. Data yang digunakan berupa data gayaberat dari BMKG tahun 2014 dengan 197 lokasi titik pengukuran yang tersebar di koordinat 6,08º-6,36º LU dan 106,68º-106,97º BT. Studi ini menggunakan metode <em>power spectral</em>  dengan mentransformasikan data dari domain jarak ke dalam domain bilangan gelombang memanfaatkan transformasi <em>Fourier</em>. Hasil penelitian dengan menggunakan metode transformasi <em>Fourier  </em>menunjukkan bahwa ketebalan sedimen di Jakarta dari arah selatan ke utara semakin besar, di sekitar Babakan ketebalan diperkirakan 92 meter, sekitar Tongkol, Jakarta Utara diperkirakan 331 meter.</p><p><strong> </strong></p><p><strong>Kata kunci</strong>: <em>power spectral</em>, anomali gayaberat, ketebalan sedimen</p><p align="center"><strong><em> </em></strong></p><p align="center"><strong><em>ABSTRACT</em></strong></p><p><em> </em></p><p><em>Studies of spectral analysis of gravity anomaly data have been carried out to estimate the thickness of sediment. In this study the author did spectral analysis of gravity anomaly data of DKI Jakarta area to know the depth of anomaly source which corresponds to the thickness of sediment. The data used in the form of gravity data from BMKG 2014 with 197 locations of measurement points spread in coordinates 6.08º - 6.36º N and 106.68º - 106.97º E. This study used the power spectral method by transforming the data from the distance domain into the wavenumber domain utilizing the Fourier transform. The result of the research using Fourier transform method shows that the thickness of sediment in Jakarta from south to north is getting bigger, in Babakan the thickness of the sediment is around 92 meter, in Tongkol, North Jakarta is around 331 meter.</em></p><p><strong><em> </em></strong></p><p><strong><em>Keywords</em></strong><em>: </em><em>power spectral, gravity anomaly, sediment thickness</em><em></em></p>


2021 ◽  
Author(s):  
David Howe

Statistical imputation is a field of study that attempts to fill missing data. It is commonly applied to population statistics whose data have no correlation with running time. For a time series, data is typically analyzed using the autocorrelation function (ACF), the Fourier transform to estimate power spectral densities (PSD), the Allan deviation (ADEV), trend extensions, and basically any analysis that depends on uniform time indexes. We explain the rationale for an imputation algorithm that fills gaps in a time series by applying a backward, inverted replica of adjacent live data. To illustrate, four intentional massive gaps that exceed 100% of the original time series are recovered. The L(f) PSD with imputation applied to the gaps is nearly indistinguishable from the original. Also, the confidence of ADEV with imputation falls within 90% of the original ADEV with mixtures of power-law noises. The algorithm in Python is included for those wishing to try it.


1981 ◽  
Vol 14 (6) ◽  
pp. 383-386 ◽  
Author(s):  
C. R. Worthington

A method of deconvolution using the modulus of the continuous Fourier transform of the unit cell is described. This method differs from previous deconvolution methods in membrane diffraction in that calculations are carried out in reciprocal space. The modulus profile is obtained from the continuous intensity transform which is itself the Fourier transform of the autocorrelation function. Sampling theorem methods are used to reconstruct the continuous Fourier transform of the unit cell. The various phase choices are examined and compared. In membrane diffraction, the autocorrelation function is derived in two distinct situations: when the interference function is broad as in the case of a few unit cells and when the membrane systems contain wide regions of constant electron density. It is concluded that, in the first situation, the derived autocorrelation function contains missing information and is incorrect. On the other hand, in the second situation, the derived autocorrelation function is a good representation of the true autocorrelation function.


Author(s):  
Jeremy J. Dawkins ◽  
David M. Bevly ◽  
Robert L. Jackson

This paper investigates the use of the Fourier transform and Wavelet transform as methods to supplement the more common root mean squared elevation and power spectral density methods of terrain characterization. Two dimensional terrain profiles were generated using the Weierstrass-Mandelbrot fractal equation. The Fourier and Wavelet transforms were used to decompose these terrains into a parameter set. A two degree of freedom quarter car model was used to evaluate the vehicle response before and after the terrain characterization. It was determined that the Fourier transform can be used to reduce the profile into the key frequency components. The Wavelet transform can effectively detect discontinuities of the profile and changes in the roughness of the profile. These two techniques can be added to current methods to yield a more robust terrain characterization.


Author(s):  
H. P. Hong ◽  
X. Z. Cui ◽  
D. Qiao

AbstractWe proposed a new iterative power and amplitude correction (IPAC) algorithm to simulate nonstationary and non-Gaussian processes. The proposed algorithm is rooted in the concept of defining the stochastic processes in the transform domain, which is elaborated and extend. The algorithm extends the iterative amplitude adjusted Fourier transform algorithm for generating surrogate and the spectral correction algorithm for simulating stationary non-Gaussian process. The IPAC algorithm can be used with different popular transforms, such as the Fourier transform, S-transform, and continuous wavelet transforms. The targets for the simulation are the marginal probability distribution function of the process and the power spectral density function of the process that is defined based on the variables in the transform domain for the adopted transform. The algorithm is versatile and efficient. Its application is illustrated using several numerical examples.


2021 ◽  
Vol 11 (6) ◽  
pp. 2582
Author(s):  
Lucas M. Martinho ◽  
Alan C. Kubrusly ◽  
Nicolás Pérez ◽  
Jean Pierre von der Weid

The focused signal obtained by the time-reversal or the cross-correlation techniques of ultrasonic guided waves in plates changes when the medium is subject to strain, which can be used to monitor the medium strain level. In this paper, the sensitivity to strain of cross-correlated signals is enhanced by a post-processing filtering procedure aiming to preserve only strain-sensitive spectrum components. Two different strategies were adopted, based on the phase of either the Fourier transform or the short-time Fourier transform. Both use prior knowledge of the system impulse response at some strain level. The technique was evaluated in an aluminum plate, effectively providing up to twice higher sensitivity to strain. The sensitivity increase depends on a phase threshold parameter used in the filtering process. Its performance was assessed based on the sensitivity gain, the loss of energy concentration capability, and the value of the foreknown strain. Signals synthesized with the time–frequency representation, through the short-time Fourier transform, provided a better tradeoff between sensitivity gain and loss of energy concentration.


2021 ◽  
Vol 262 ◽  
pp. 117928
Author(s):  
Shusaku Nakajima ◽  
Shuhei Horiuchi ◽  
Akifumi Ikehata ◽  
Yuichi Ogawa

Sign in / Sign up

Export Citation Format

Share Document