Abstract:When using traditional Euler deconvolution optimization strategies, it is difficult to distinguish between anomalies and their corresponding Euler tails (those solutions are often distributed outside the anomaly source, forming "tail"-shaped spurious solutions, i.e., misplaced Euler solutions, which must be removed or marked) with only the structural index. The nonparametric estimation method based on the normalized B-spline probability density (BSS) is used to separate the Euler solution clusters and mark different anomaly sources according to the similarity and density characteristics of the Euler solutions. For display purposes, the BSS needs to map the samples onto the estimation grid at the points where density will be estimated in order to obtain the probability density distribution. However, if the size of the samples or the estimation grid is too large, this process can lead to high levels of memory consumption and excessive computation times. To address this issue, a fast linear binning approximation algorithm is introduced in the BSS to speed up the computation process and save time. Subsequently, the sample data are quickly projected onto the estimation grid to facilitate the discrete convolution between the grid and the density function using a fast Fourier transform. A method involving multivariate B-spline probability density estimation based on the FFT (BSSFFT), in conjunction with fast linear binning appropriation, is proposed in this paper. The results of two random normal distributions show the correctness of the BSS and BSSFFT algorithms, which is verified via a comparison with the true probability density function (pdf) and Gaussian kernel smoothing estimation algorithms. Then, the Euler solutions of the two synthetic models are analyzed using the BSS and BSSFFT algorithms. The results are consistent with their theoretical values, which verify their correctness regarding Euler solutions. Finally, the BSSFFT is applied to Bishop 5X data, and the numerical results show that the comprehensive analysis of the 3D probability density distributions using the BSSFFT algorithm, derived from the Euler solution subset of x0,y0,z0, can effectively separate and locate adjacent anomaly sources, demonstrating strong adaptability.

Utilizing Bochners Theorem for Constrained Evaluation of Missing Fourier Data

Fast and Quantitative Reconstruction of 4D SPECT for Fan-beam Geometry Based on Novikov's Inverse Formula

Reducing the Sampling Burden of Fourier Sensing with a Non-rectangular Field-of-View

Uncertainty quantification for sparse Fourier recovery

With or Without Replacement? Improving Confidence in Fourier Imaging

Mathematical Equivalence of Zero-Padding Interpolation and Circular Sampling Theorem Interpolation with Implications for Direct Fourier Image Reconstruction.

Nonlinear Fourier analysis for discontinuous conductivities: computational results

Image Reconstruction From Sparse Data In Echo-Planar Imaging

FFT-Based Probability Density Imaging of Euler Solutions

Optimal density compensation factors for the reconstruction of the Fourier transform of bandlimited functions

Kaczmarz Iterative Projection and Nonuniform Sampling with Complexity Estimates

Frequency Sub-Sampling of Ultrasound Non-Destructive Measurements: Acquisition, Reconstruction and Performance

A Fourier representation of kernel Stein discrepancy with application to Goodness-of-Fit tests for measures on infinite dimensional Hilbert spaces

High-resolution imaging using partial separability of spatiotemporal signals with a novel data sampling scheme

Pseudo-inverse reconstruction of bandlimited signals from nonuniform generalized samples with orthogonal kernels

Accelerating Non-Cartesian MRI Reconstruction Convergence using k-space Preconditioning

Recovery of Point Clouds on Surfaces: Application to Image Reconstruction

Improved Reconstruction for Fourier-Sparse Signals

Reconstruction of large, irregularly sampled multidimensional images. A tensor-based approach

Signal reconstruction using determinantal sampling

A Riemann--Hilbert approach to computing the inverse spectral map for measures supported on disjoint intervals