In seismic data processing, blind deconvolution is a key technology. Introduced in this paper is a flow of one kind of blind deconvolution. The optimal precondition conjugate gradients (PCG) in Kyrlov subspace is als...In seismic data processing, blind deconvolution is a key technology. Introduced in this paper is a flow of one kind of blind deconvolution. The optimal precondition conjugate gradients (PCG) in Kyrlov subspace is also used to improve the stability of the algorithm. The computation amount is greatly decreased.展开更多
Successful blind image deconvolution algorithms require the exact estimation of the Point Spread Function size, PSF. In the absence of any priori information about the imagery system and the true image, this estimatio...Successful blind image deconvolution algorithms require the exact estimation of the Point Spread Function size, PSF. In the absence of any priori information about the imagery system and the true image, this estimation is normally done by trial and error experimentation, until an acceptable restored image quality is obtained. This paper, presents an exact estimation of the PSF size, which yields the optimum restored image quality for both noisy and noiseless images. It is based on evaluating the detail energy of the wave packet decomposition of the blurred image. The minimum detail energies occur at the optimum PSF size. Having accurately estimated the PSF, the paper also proposes a fast double updating algorithm for improving the quality of the restored image. This is achieved by the least squares minimization of a system of linear equations that minimizes some error functions derived from the blurred image. Moreover, a technique is also proposed to improve the sharpness of the deconvolved images, by constrained maximization of some of the detail wavelet packet energies. Simulation results of several examples have verified that the proposed technique manages to yield a sharper image with higher PSNR than classical approaches.展开更多
Former frequency-domain blind devolution algorithms need to consider a large number of frequency bins and recover the sources in different orders and with different amplitudes in each frequency bin,so they suffer from...Former frequency-domain blind devolution algorithms need to consider a large number of frequency bins and recover the sources in different orders and with different amplitudes in each frequency bin,so they suffer from permutation and amplitude indeterminacy troubles. Based on sliding discrete Fourier transform,the presented deconvolution algorithm can directly recover time-domain sources from frequency-domain convolutive model using single frequency bin. It only needs to execute blind sepa-ration of instantaneous mixture once there are no permutation and amplitude indeterminacy troubles. Compared with former algorithms,the algorithm greatly reduces the computation cost as only one frequency bin is considered. Its good and robust per-formance is demonstrated by simulations when the signal-to-noise-ratio is high.展开更多
As an ill-posed problem, multiframe blind super resolution imaging recovers a high resolution image from a group of low resolution images with some degradations when the information of blur kernel is limited. Note tha...As an ill-posed problem, multiframe blind super resolution imaging recovers a high resolution image from a group of low resolution images with some degradations when the information of blur kernel is limited. Note that the quality of the recovered image is influenced more by the accuracy of blur estimation than an advanced regularization. We study the traditional model of the multiframe super resolution and modify it for blind deblurring. Based on the analysis, we proposed two algorithms. The first one is based on the total variation blind deconvolution algorithm and formulated as a functional for optimization with the regularization of blur. Based on the alternating minimization and the gradient descent algorithm, the high resolution image and the unknown blur kernel are estimated iteratively. By using the median shift and add operator, the second algorithm is more robust to the outlier influence. The MSAA initialization simplifies the interpolation process to reconstruct the blurred high resolution image for blind deblurring and improves the accuracy of blind super resolution imaging. The experimental results demonstrate the superiority and accuracy of our novel algorithms.展开更多
An algorithm is presented for image prior combinations based blind deconvolution and applied to astronomical images.Using a hierarchical Bayesian framework, the unknown original image and all required algorithmic para...An algorithm is presented for image prior combinations based blind deconvolution and applied to astronomical images.Using a hierarchical Bayesian framework, the unknown original image and all required algorithmic parameters are estimated simultaneously. Through utilization of variational Bayesian analysis,approximations of the posterior distributions on each unknown are obtained by minimizing the Kullback-Leibler(KL) distance, thus providing uncertainties of the estimates during the restoration process. Experimental results on both synthetic images and real astronomical images demonstrate that the proposed approaches compare favorably to other state-of-the-art reconstruction methods.展开更多
In this paper, we propose a novel seismic blind deconvolution approach based on the Spearman’s rho in the case of band-limited seismic data with a low dominant frequency and short data records. The Spearman’s rho is...In this paper, we propose a novel seismic blind deconvolution approach based on the Spearman’s rho in the case of band-limited seismic data with a low dominant frequency and short data records. The Spearman’s rho is a measure of the dependence between two continuous random variables without the influence of the marginal distributions, by which a new criterion for blind deconvolution is constructed. The optimization program for new criterion of blind deconvolution is performed by applying Neidell’s wavelet model to the inverse filter. The noise-free and noisy synthetic data, onshore seismic trace in the Ordos Basin, and offshore stacked section in the Bohai Bay Basin examples show good results of the method.展开更多
In order to alleviate the shortcomings of most blind deconvolution algorithms,this paper proposes an improved fast algorithm for blind deconvolution based on decorrelation technique and broadband block matrix.Althougt...In order to alleviate the shortcomings of most blind deconvolution algorithms,this paper proposes an improved fast algorithm for blind deconvolution based on decorrelation technique and broadband block matrix.Althougth the original algorithm can overcome the shortcomings of current blind deconvolution algorithms,it has a constraint that the number of the source signals must be less than that of the channels.The improved algorithm deletes this constraint by using decorrelation technique.Besides,the improved algorithm raises the separation speed in terms of improving the computing methods of the output signal matrix.Simulation results demonstrate the validation and fast separation of the improved algorithm.展开更多
An iterative separation approach, i.e. source signals are extracted and removed one by one, is proposed for multichannel blind deconvolution of colored signals. Each source signal is extracted in two stages: a filtere...An iterative separation approach, i.e. source signals are extracted and removed one by one, is proposed for multichannel blind deconvolution of colored signals. Each source signal is extracted in two stages: a filtered version of the source signal is first obtained by solving the generalized eigenvalue problem, which is then followed by a single channel blind deconvolution based on ensemble learning. Simulation demonstrates the capability of the approach to perform efficient mutichannel blind deconvolution.展开更多
Recently we have developed an eigenvector method (EVM) which can achieve the blind deconvolution (BD) for MIMO systems. One of attractive features of the proposed algorithm is that the BD can be achieved by calculatin...Recently we have developed an eigenvector method (EVM) which can achieve the blind deconvolution (BD) for MIMO systems. One of attractive features of the proposed algorithm is that the BD can be achieved by calculating the eigenvectors of a matrix relevant to it. However, the performance accuracy of the EVM depends highly on computational results of the eigenvectors. In this paper, by modifying the EVM, we propose an algorithm which can achieve the BD without calculating the eigenvectors. Then the pseudo-inverse which is needed to carry out the BD is calculated by our proposed matrix pseudo-inversion lemma. Moreover, using a combination of the conventional EVM and the modified EVM, we will show its performances comparing with each EVM. Simulation results will be presented for showing the effectiveness of the proposed methods.展开更多
Vehicle reidentification is an elegant solution for gathering several pieces of valuable traffic information, e.g., space mean speed, travel time, vehicle tracking, and origin/destination data. Recently, a number of v...Vehicle reidentification is an elegant solution for gathering several pieces of valuable traffic information, e.g., space mean speed, travel time, vehicle tracking, and origin/destination data. Recently, a number of vehiclereidentification algorithms utilizing inductive loop signals have been proposed to take advantage of the widespread availability of loop detectors. These algorithms, however, all directly utilize the raw inductance signals for pattern matching and feature extraction without deconvolution. The raw loop signals are essentially a convolved output between the true vehicle inductance signature and the loop system function, and thus a deconvolution is needed in order to expose the detailed features of individual vehicles. The purpose of this paper is to present a recent investigation on restoration of true inductance signatures by applying a blind deconvolution process. The main advantage of blind deconvolution over the conventional deconvolution is that the computation does not require modeling of a precise loop-detector system function. Experimental results show that the proposed blind deconvolution reveals much more detailed features of inductance signals and, as a result, increases the vehicle reidentification accuracy.展开更多
In this paper, the inverse problem of reconstructing reflectivity function of a medium is examined within a blind deconvolution framework. The ultrasound pulse is estimated using higher-order statistics, and Wiener fi...In this paper, the inverse problem of reconstructing reflectivity function of a medium is examined within a blind deconvolution framework. The ultrasound pulse is estimated using higher-order statistics, and Wiener filter is used to obtain the ultrasonic reflectivity function through wavelet-based models. A new approach to the parameter estimation of the inverse filtering step is proposed in the nondestructive evaluation field, which is based on the theory of Fourier-Wavelet regularized deconvolution (ForWaRD). This new approach can be viewed as a solution to the open problem of adaptation of the ForWaRD framework to perform the convolution kernel estimation and deconvolution interdependently. The results indicate stable solutions of the esti- mated pulse and an improvement in the radio-frequency (RF) signal taking into account its signal-to-noise ratio (SNR) and axial resolution. Simulations and experiments showed that the proposed approach can provide robust and optimal estimates of the reflectivity function.展开更多
In the communication field, during transmission, a source signal undergoes a convolutive distortion between its symbols and the channel impulse response. This distortion is referred to as Intersymbol Interference (ISI...In the communication field, during transmission, a source signal undergoes a convolutive distortion between its symbols and the channel impulse response. This distortion is referred to as Intersymbol Interference (ISI) and can be reduced significantly by applying a blind adaptive deconvolution process (blind adaptive equalizer) on the distorted received symbols. But, since the entire blind deconvolution process is carried out with no training symbols and the channel’s coefficients are obviously unknown to the receiver, no actual indication can be given (via the mean square error (MSE) or ISI expression) during the deconvolution process whether the blind adaptive equalizer succeeded to remove the heavy ISI from the transmitted symbols or not. Up to now, the output of a convolution and deconvolution process was mainly investigated from the ISI point of view. In this paper, the output of a convolution and deconvolution process is inspected from the leading digit point of view. Simulation results indicate that for the 4PAM (Pulse Amplitude Modulation) and 16QAM (Quadrature Amplitude Modulation) input case, the number “1” is the leading digit at the output of a convolution and deconvolution process respectively as long as heavy ISI exists. However, this leading digit does not follow exactly Benford’s Law but follows approximately the leading digit (digit 1) of a Gaussian process for independent identically distributed input symbols and a channel with many coefficients.展开更多
基金With the support of the key project of Knowledge Innovation, CAS(KZCX1-y01, KZCX-SW-18), Fund of the China National Natural Sciences and the Daqing Oilfield with Grant No. 49894190
文摘In seismic data processing, blind deconvolution is a key technology. Introduced in this paper is a flow of one kind of blind deconvolution. The optimal precondition conjugate gradients (PCG) in Kyrlov subspace is also used to improve the stability of the algorithm. The computation amount is greatly decreased.
文摘Successful blind image deconvolution algorithms require the exact estimation of the Point Spread Function size, PSF. In the absence of any priori information about the imagery system and the true image, this estimation is normally done by trial and error experimentation, until an acceptable restored image quality is obtained. This paper, presents an exact estimation of the PSF size, which yields the optimum restored image quality for both noisy and noiseless images. It is based on evaluating the detail energy of the wave packet decomposition of the blurred image. The minimum detail energies occur at the optimum PSF size. Having accurately estimated the PSF, the paper also proposes a fast double updating algorithm for improving the quality of the restored image. This is achieved by the least squares minimization of a system of linear equations that minimizes some error functions derived from the blurred image. Moreover, a technique is also proposed to improve the sharpness of the deconvolved images, by constrained maximization of some of the detail wavelet packet energies. Simulation results of several examples have verified that the proposed technique manages to yield a sharper image with higher PSNR than classical approaches.
基金Project (No. 2005EB040486) supported by the National Torch Program of China
文摘Former frequency-domain blind devolution algorithms need to consider a large number of frequency bins and recover the sources in different orders and with different amplitudes in each frequency bin,so they suffer from permutation and amplitude indeterminacy troubles. Based on sliding discrete Fourier transform,the presented deconvolution algorithm can directly recover time-domain sources from frequency-domain convolutive model using single frequency bin. It only needs to execute blind sepa-ration of instantaneous mixture once there are no permutation and amplitude indeterminacy troubles. Compared with former algorithms,the algorithm greatly reduces the computation cost as only one frequency bin is considered. Its good and robust per-formance is demonstrated by simulations when the signal-to-noise-ratio is high.
基金Supported by the National Natural Science Foundation of China(No.61340034)the Research Program of Application Foundation and Advanced Technology of Tianjin(No.13JCYBJC15600)
文摘As an ill-posed problem, multiframe blind super resolution imaging recovers a high resolution image from a group of low resolution images with some degradations when the information of blur kernel is limited. Note that the quality of the recovered image is influenced more by the accuracy of blur estimation than an advanced regularization. We study the traditional model of the multiframe super resolution and modify it for blind deblurring. Based on the analysis, we proposed two algorithms. The first one is based on the total variation blind deconvolution algorithm and formulated as a functional for optimization with the regularization of blur. Based on the alternating minimization and the gradient descent algorithm, the high resolution image and the unknown blur kernel are estimated iteratively. By using the median shift and add operator, the second algorithm is more robust to the outlier influence. The MSAA initialization simplifies the interpolation process to reconstruct the blurred high resolution image for blind deblurring and improves the accuracy of blind super resolution imaging. The experimental results demonstrate the superiority and accuracy of our novel algorithms.
文摘An algorithm is presented for image prior combinations based blind deconvolution and applied to astronomical images.Using a hierarchical Bayesian framework, the unknown original image and all required algorithmic parameters are estimated simultaneously. Through utilization of variational Bayesian analysis,approximations of the posterior distributions on each unknown are obtained by minimizing the Kullback-Leibler(KL) distance, thus providing uncertainties of the estimates during the restoration process. Experimental results on both synthetic images and real astronomical images demonstrate that the proposed approaches compare favorably to other state-of-the-art reconstruction methods.
文摘In this paper, we propose a novel seismic blind deconvolution approach based on the Spearman’s rho in the case of band-limited seismic data with a low dominant frequency and short data records. The Spearman’s rho is a measure of the dependence between two continuous random variables without the influence of the marginal distributions, by which a new criterion for blind deconvolution is constructed. The optimization program for new criterion of blind deconvolution is performed by applying Neidell’s wavelet model to the inverse filter. The noise-free and noisy synthetic data, onshore seismic trace in the Ordos Basin, and offshore stacked section in the Bohai Bay Basin examples show good results of the method.
基金Natural Science Fund of Anhui Province of China (050420101)
文摘In order to alleviate the shortcomings of most blind deconvolution algorithms,this paper proposes an improved fast algorithm for blind deconvolution based on decorrelation technique and broadband block matrix.Althougth the original algorithm can overcome the shortcomings of current blind deconvolution algorithms,it has a constraint that the number of the source signals must be less than that of the channels.The improved algorithm deletes this constraint by using decorrelation technique.Besides,the improved algorithm raises the separation speed in terms of improving the computing methods of the output signal matrix.Simulation results demonstrate the validation and fast separation of the improved algorithm.
基金Supported by the National Natural Science Foundation of China(No.60072048)the Doctoral Program Fund(No.20010561007)
文摘An iterative separation approach, i.e. source signals are extracted and removed one by one, is proposed for multichannel blind deconvolution of colored signals. Each source signal is extracted in two stages: a filtered version of the source signal is first obtained by solving the generalized eigenvalue problem, which is then followed by a single channel blind deconvolution based on ensemble learning. Simulation demonstrates the capability of the approach to perform efficient mutichannel blind deconvolution.
文摘Recently we have developed an eigenvector method (EVM) which can achieve the blind deconvolution (BD) for MIMO systems. One of attractive features of the proposed algorithm is that the BD can be achieved by calculating the eigenvectors of a matrix relevant to it. However, the performance accuracy of the EVM depends highly on computational results of the eigenvectors. In this paper, by modifying the EVM, we propose an algorithm which can achieve the BD without calculating the eigenvectors. Then the pseudo-inverse which is needed to carry out the BD is calculated by our proposed matrix pseudo-inversion lemma. Moreover, using a combination of the conventional EVM and the modified EVM, we will show its performances comparing with each EVM. Simulation results will be presented for showing the effectiveness of the proposed methods.
文摘Vehicle reidentification is an elegant solution for gathering several pieces of valuable traffic information, e.g., space mean speed, travel time, vehicle tracking, and origin/destination data. Recently, a number of vehiclereidentification algorithms utilizing inductive loop signals have been proposed to take advantage of the widespread availability of loop detectors. These algorithms, however, all directly utilize the raw inductance signals for pattern matching and feature extraction without deconvolution. The raw loop signals are essentially a convolved output between the true vehicle inductance signature and the loop system function, and thus a deconvolution is needed in order to expose the detailed features of individual vehicles. The purpose of this paper is to present a recent investigation on restoration of true inductance signatures by applying a blind deconvolution process. The main advantage of blind deconvolution over the conventional deconvolution is that the computation does not require modeling of a precise loop-detector system function. Experimental results show that the proposed blind deconvolution reveals much more detailed features of inductance signals and, as a result, increases the vehicle reidentification accuracy.
基金Project (No. PRC 03-41/2003) supported by the Ministry of Con-struction of Cuba
文摘In this paper, the inverse problem of reconstructing reflectivity function of a medium is examined within a blind deconvolution framework. The ultrasound pulse is estimated using higher-order statistics, and Wiener filter is used to obtain the ultrasonic reflectivity function through wavelet-based models. A new approach to the parameter estimation of the inverse filtering step is proposed in the nondestructive evaluation field, which is based on the theory of Fourier-Wavelet regularized deconvolution (ForWaRD). This new approach can be viewed as a solution to the open problem of adaptation of the ForWaRD framework to perform the convolution kernel estimation and deconvolution interdependently. The results indicate stable solutions of the esti- mated pulse and an improvement in the radio-frequency (RF) signal taking into account its signal-to-noise ratio (SNR) and axial resolution. Simulations and experiments showed that the proposed approach can provide robust and optimal estimates of the reflectivity function.
文摘In the communication field, during transmission, a source signal undergoes a convolutive distortion between its symbols and the channel impulse response. This distortion is referred to as Intersymbol Interference (ISI) and can be reduced significantly by applying a blind adaptive deconvolution process (blind adaptive equalizer) on the distorted received symbols. But, since the entire blind deconvolution process is carried out with no training symbols and the channel’s coefficients are obviously unknown to the receiver, no actual indication can be given (via the mean square error (MSE) or ISI expression) during the deconvolution process whether the blind adaptive equalizer succeeded to remove the heavy ISI from the transmitted symbols or not. Up to now, the output of a convolution and deconvolution process was mainly investigated from the ISI point of view. In this paper, the output of a convolution and deconvolution process is inspected from the leading digit point of view. Simulation results indicate that for the 4PAM (Pulse Amplitude Modulation) and 16QAM (Quadrature Amplitude Modulation) input case, the number “1” is the leading digit at the output of a convolution and deconvolution process respectively as long as heavy ISI exists. However, this leading digit does not follow exactly Benford’s Law but follows approximately the leading digit (digit 1) of a Gaussian process for independent identically distributed input symbols and a channel with many coefficients.