Geophysical data sparse reconstruction based on L0-norm minimization

Missing data are a problem in geophysical surveys, and interpolation and reconstruction of missing data is part of the data processing and interpretation. Based on the sparseness of the geophysical data or the transform domain, we can improve the accuracy and stability of the reconstruction by transforming it to a sparse optimization problem. In this paper, we propose a mathematical model for the sparse reconstruction of data based on the LO-norm minimization. Furthermore, we discuss two types of the approximation algorithm for the LO- norm minimization according to the size and characteristics of the geophysical data： namely, the iteratively reweighted least-squares algorithm and the fast iterative hard thresholding algorithm. Theoretical and numerical analysis showed that applying the iteratively reweighted least-squares algorithm to the reconstruction of potential field data exploits its fast convergence rate, short calculation time, and high precision, whereas the fast iterative hard thresholding algorithm is more suitable for processing seismic data, moreover, its computational efficiency is better than that of the traditional iterative hard thresholding algorithm.

Efficient Concurrent L1-Minimization Solvers on GPUs

Given that the concurrent L1-minimization(L1-min)problem is often required in some real applications,we investigate how to solve it in parallel on GPUs in this paper.First,we propose a novel self-adaptive warp implementation of the matrix-vector multiplication(Ax)and a novel self-adaptive thread implementation of the matrix-vector multiplication(ATx),respectively,on the GPU.The vector-operation and inner-product decision trees are adopted to choose the optimal vector-operation and inner-product kernels for vectors of any size.Second,based on the above proposed kernels,the iterative shrinkage-thresholding algorithm is utilized to present two concurrent L1-min solvers from the perspective of the streams and the thread blocks on a GPU,and optimize their performance by using the new features of GPU such as the shuffle instruction and the read-only data cache.Finally,we design a concurrent L1-min solver on multiple GPUs.The experimental results have validated the high effectiveness and good performance of our proposed methods.

Synthetic aperture radar imaging based on attributed scatter model using sparse recovery techniques

The sparse recovery algorithms formulate synthetic aperture radar （SAR） imaging problem in terms of sparse representation （SR） of a small number of strong scatters＇ positions among a much large number of potential scatters＇ positions, and provide an effective approach to improve the SAR image resolution. Based on the attributed scatter center model, several experiments were performed with different practical considerations to evaluate the performance of five representative SR techniques, namely, sparse Bayesian learning （SBL）, fast Bayesian matching pursuit （FBMP）, smoothed 10 norm method （SL0）, sparse reconstruction by separable approximation （SpaRSA）, fast iterative shrinkage-thresholding algorithm （FISTA）, and the parameter settings in five SR algorithms were discussed. In different situations, the performances of these algorithms were also discussed. Through the comparison of MSE and failure rate in each algorithm simulation, FBMP and SpaRSA are found suitable for dealing with problems in the SAR imaging based on attributed scattering center model. Although the SBL is time-consuming, it always get better performance when related to failure rate and high SNR.

A Fast High Order Iterative Solver for the Electromagnetic Scattering by Open Cavities Filled with the Inhomogeneous Media

The scattering of the open cavity filled with the inhomogeneous media is studied.The problem is discretized with a fourth order finite difference scheme and the immersed interfacemethod,resulting in a linear system of equations with the high order accurate solutions in the whole computational domain.To solve the system of equations,we design an efficient iterative solver,which is based on the fast Fourier transformation,and provides an ideal preconditioner for Krylov subspace method.Numerical experiments demonstrate the capability of the proposed fast high order iterative solver.

Surface registration algorithm for rapid detection of surface thermal deformation of paraboloid antennas

In order to obtain and master the surface thermal deformation of paraboloid antennas,a fast iterative closest point( FICP) algorithm based on design coordinate guidance is proposed,which can satisfy the demands of rapid detection for surface thermal deformation. Firstly,the basic principle of the ICP algorithm for registration of a free surface is given,and the shortcomings of the ICP algorithm in the registration of surface are analysed,such as its complex computation,long calculation time,low efficiency,and relatively strict initial registration position. Then an improved FICP algorithm based on design coordinate guidance is proposed. Finally,the FICP algorithm is applied to the fast registration test for the surface thermal deformation of a paraboloid antenna. Results indicate that the approach offers better performance with regard to fast surface registration and the algorithm is more simple,efficient,and easily realized in practical engineering application.