A new method for receiver function inversion by wavelet transformation is presented in this paper. Receiver func-tion is expanded to different scales with different resolution by wavelet transformation. After an initi...A new method for receiver function inversion by wavelet transformation is presented in this paper. Receiver func-tion is expanded to different scales with different resolution by wavelet transformation. After an initial model be-ing taken, a generalized least-squares inversion procedure is gradually carried out for receiver function from low to high scale, with the inversion result for low order receiver function as the initial model for high order. A neighborhood containing the global minimum is firstly searched from low scale receiver function, and will gradu-ally focus at the global minimum by introducing high scale information of receiver function. With the gradual ad-dition of high wave-number to smooth background velocity structure, wavelet transformation can keep the inver-sion result converge to the global minimum, reduce to certain extent the dependence of inversion result on the initial model, overcome the nonuniqueness of generalized least-squares inversion, and obtain reliable crustal and upper mantle velocity with high resolution.展开更多
In this paper, we present a new image compression method based on the direct and inverse F1-transform, a generalization of the concept of fuzzy transform. Under weak compression rates, this method improves the quality...In this paper, we present a new image compression method based on the direct and inverse F1-transform, a generalization of the concept of fuzzy transform. Under weak compression rates, this method improves the quality of the images with respect to the classical method based on the fuzzy transform.展开更多
A method for reconstructing crustal velocity structure using the optimization of stacking receiver function amplitude in the depth domain,named common conversion amplitude(CCA)inversion,is presented.The conversion amp...A method for reconstructing crustal velocity structure using the optimization of stacking receiver function amplitude in the depth domain,named common conversion amplitude(CCA)inversion,is presented.The conversion amplitude in the depth domain,which represents the impedance change in the medium,is obtained by assigning the receiver function amplitude to the corresponding conversion position where the P-to-S conversion occurred.Utilizing the conversion amplitude variation with depth as an optimization objective,imposing reliable prior constraints on the structural model frame and velocity range,and adopting a stepwise search inversion technique,this method efficiently weakens the tendency of easily falling into the local extremum in conventional receiver function inversion.Synthetic tests show that the CCA inversion can reconstruct complex crustal velocity structures well and is especially suitable for revealing crustal evolution by estimating diverse velocity distributions.Its performance in reconstructing crustal structure is superior to that of the conventional receiver function imaging method.展开更多
The local reconstruction from truncated projection data is one area of interest in image reconstruction for com- puted tomography (CT), which creates the possibility for dose reduction. In this paper, a filtered-bac...The local reconstruction from truncated projection data is one area of interest in image reconstruction for com- puted tomography (CT), which creates the possibility for dose reduction. In this paper, a filtered-backprojection (FBP) algorithm based on the Radon inversion transform is presented to deal with the three-dimensional (3D) local recon- struction in the circular geometry. The algorithm achieves the data filtering in two steps. The first step is the derivative of projections, which acts locally on the data and can thus be carried out accurately even in the presence of data trun- cation. The second step is the nonlocal Hilbert filtering. The numerical simulations and the real data reconstructions have been conducted to validate the new reconstruction algorithm. Compared with the approximate truncation resistant algorithm for computed tomography (ATRACT), not only it has a comparable ability to restrain truncation artifacts, but also its reconstruction efficiency is improved. It is about twice as fast as that of the ATRACT. Therefore, this work provides a simple and efficient approach for the approximate reconstruction from truncated projections in the circular cone-beam CT.展开更多
Based on anisotropic total variation regularization(ATVR), a nonnegativity and support constraints recursive inverse filtering(NAS-RIF) blind restoration method is proposed to enhance the quality of optical coherence ...Based on anisotropic total variation regularization(ATVR), a nonnegativity and support constraints recursive inverse filtering(NAS-RIF) blind restoration method is proposed to enhance the quality of optical coherence tomography(OCT) image. First, ATVR is introduced into the cost function of NAS-RIF to improve the noise robustness and retain the details in the image.Since the split Bregman iterative is used to optimize the ATVR based cost function, the ATVR based NAS-RIF blind restoration method is then constructed. Furthermore, combined with the geometric nonlinear diffusion filter and the Poisson-distribution-based minimum error thresholding, the ATVR based NAS-RIF blind restoration method is used to realize the blind OCT image restoration. The experimental results demonstrate that the ATVR based NAS-RIF blind restoration method can successfully retain the details in the OCT images. In addition, the signal-to-noise ratio of the blind restored OCT images can be improved, along with the noise robustness.展开更多
Radio astronomy radio telescope plays the role of a linear operator, affecting the function that describes the object of research, formation of image of a monitored object. This paper presents methods for reconstructi...Radio astronomy radio telescope plays the role of a linear operator, affecting the function that describes the object of research, formation of image of a monitored object. This paper presents methods for reconstruction and correction of solar radio images using the algorithm of rejections, the updated Weiner-filter, and the method CLEAN designed by Hegbomom (Pseudonym, 2009) for point sources. It is the process of numerical convolution in signal handling, an algorithm for separating weak-contrast formations on the solar which represents most points of the actual limb by using the ellipse equation. Consequently, the filling algorithm is applied by moving from the center to the ellipse points and filling each point by solar image data. Finally, a linear limb-darkening expression is used to remove the limb darkening. Different examples of the intermediate and final results are presented in addition to the developed algorithm.展开更多
To improve the global search ability and imaging quality of electrical resistivity imaging(ERI) inversion, a two-stage learning ICPSO algorithm of radial basis function neural network(RBFNN) based on information crite...To improve the global search ability and imaging quality of electrical resistivity imaging(ERI) inversion, a two-stage learning ICPSO algorithm of radial basis function neural network(RBFNN) based on information criterion(IC) and particle swarm optimization(PSO) is presented. In the proposed method, IC is applied to obtain the hidden layer structure by calculating the optimal IC value automatically and PSO algorithm is used to optimize the centers and widths of the radial basis functions in the hidden layer. Meanwhile, impacts of different information criteria to the inversion results are compared, and an implementation of the proposed ICPSO algorithm is given. The optimized neural network has one hidden layer with 261 nodes selected by AKAIKE's information criterion(AIC) and it is trained on 32 data sets and tested on another 8 synthetic data sets. Two complex synthetic examples are used to verify the feasibility and effectiveness of the proposed method with two learning stages. The results show that the proposed method has better performance and higher imaging quality than three-layer and four-layer back propagation neural networks(BPNNs) and traditional least square(LS) inversion.展开更多
We survey fundamental concepts for inverse programming and then present the Universal Resolving Algorithm, an algorithm for inverse computation in a first order, functional programming language. We discuss the key co...We survey fundamental concepts for inverse programming and then present the Universal Resolving Algorithm, an algorithm for inverse computation in a first order, functional programming language. We discuss the key concepts of the algorithm, including a three step approach based on the notion of a perfect process tree, and demonstrate our implementation with several examples of inverse computation.展开更多
Reflection imaging results generally reveal large-scale continuous geological information,and it is difficult to identify small-scale geological bodies such as breakpoints,pinch points,small fault blocks,caves,and fra...Reflection imaging results generally reveal large-scale continuous geological information,and it is difficult to identify small-scale geological bodies such as breakpoints,pinch points,small fault blocks,caves,and fractures,etc.Diffraction imaging is an important method to identify small-scale geological bodies and it has higher resolution than reflection imaging.In the common-offset domain,reflections are mostly expressed as smooth linear events,whereas diffractions are characterized by hyperbolic events.This paper proposes a diffraction extraction method based on double sparse transforms.The linear events can be sparsely expressed by the high-resolution linear Radon transform,and the curved events can be sparsely expressed by the Curvelet transform.A sparse inversion model is built and the alternating direction method is used to solve the inversion model.Simulation data and field data experimental results proved that the diffractions extraction method based on double sparse transforms can effectively improve the imaging quality of faults and other small-scale geological bodies.展开更多
In this paper,we first construct generalized q^2-cosine,q^2-sine and q^2-exponential functions.We then use q^2-exponential function in order to define and investigate a q^2-Fourier transform.We establish q-analogues o...In this paper,we first construct generalized q^2-cosine,q^2-sine and q^2-exponential functions.We then use q^2-exponential function in order to define and investigate a q^2-Fourier transform.We establish q-analogues of inversion and Plancherel theorems.展开更多
the technique of image processing and analysis of gravity and magnetic data is one of themost effective ways to extract geological information from gravity and msanetic data. The presentpaper investigates, from an ang...the technique of image processing and analysis of gravity and magnetic data is one of themost effective ways to extract geological information from gravity and msanetic data. The presentpaper investigates, from an angle of generalized joint inversion, thc methods and procedures ofcomprehensive processing of multi-source geological image , and a specific example in Huai Nan coalfield is given here as well.展开更多
The adaptive reconstruction for the lost information of the rectangular image area is very important for the robust transmission and restoration of the image. In this paper, a new reconstruction method based on the Di...The adaptive reconstruction for the lost information of the rectangular image area is very important for the robust transmission and restoration of the image. In this paper, a new reconstruction method based on the Discrete Cosine Transform (DCT) domain has been put forward. According to the low pass character of the human visual system and the energy distribution of the DCT coefficients on the rectangular boundary, the DCT coefficients of the rectangular image area are adaptively selected and recovered. After the Inverse Discrete Cosine Transform (IDCT), the lost information of the rectangular image area can be reconstructed. The experiments have demonstrated that the subjective and objective qualities of the reconstructed images are enhanced greatly than before.展开更多
基金National Nature Science Foundation of China (49974021).
文摘A new method for receiver function inversion by wavelet transformation is presented in this paper. Receiver func-tion is expanded to different scales with different resolution by wavelet transformation. After an initial model be-ing taken, a generalized least-squares inversion procedure is gradually carried out for receiver function from low to high scale, with the inversion result for low order receiver function as the initial model for high order. A neighborhood containing the global minimum is firstly searched from low scale receiver function, and will gradu-ally focus at the global minimum by introducing high scale information of receiver function. With the gradual ad-dition of high wave-number to smooth background velocity structure, wavelet transformation can keep the inver-sion result converge to the global minimum, reduce to certain extent the dependence of inversion result on the initial model, overcome the nonuniqueness of generalized least-squares inversion, and obtain reliable crustal and upper mantle velocity with high resolution.
文摘In this paper, we present a new image compression method based on the direct and inverse F1-transform, a generalization of the concept of fuzzy transform. Under weak compression rates, this method improves the quality of the images with respect to the classical method based on the fuzzy transform.
基金financially supported by the National Natural Science Foundation of China(Grant 91755214).
文摘A method for reconstructing crustal velocity structure using the optimization of stacking receiver function amplitude in the depth domain,named common conversion amplitude(CCA)inversion,is presented.The conversion amplitude in the depth domain,which represents the impedance change in the medium,is obtained by assigning the receiver function amplitude to the corresponding conversion position where the P-to-S conversion occurred.Utilizing the conversion amplitude variation with depth as an optimization objective,imposing reliable prior constraints on the structural model frame and velocity range,and adopting a stepwise search inversion technique,this method efficiently weakens the tendency of easily falling into the local extremum in conventional receiver function inversion.Synthetic tests show that the CCA inversion can reconstruct complex crustal velocity structures well and is especially suitable for revealing crustal evolution by estimating diverse velocity distributions.Its performance in reconstructing crustal structure is superior to that of the conventional receiver function imaging method.
基金Project supported by the National High Technology Research and Development Program of China (Grant No. 2012AA011603)
文摘The local reconstruction from truncated projection data is one area of interest in image reconstruction for com- puted tomography (CT), which creates the possibility for dose reduction. In this paper, a filtered-backprojection (FBP) algorithm based on the Radon inversion transform is presented to deal with the three-dimensional (3D) local recon- struction in the circular geometry. The algorithm achieves the data filtering in two steps. The first step is the derivative of projections, which acts locally on the data and can thus be carried out accurately even in the presence of data trun- cation. The second step is the nonlocal Hilbert filtering. The numerical simulations and the real data reconstructions have been conducted to validate the new reconstruction algorithm. Compared with the approximate truncation resistant algorithm for computed tomography (ATRACT), not only it has a comparable ability to restrain truncation artifacts, but also its reconstruction efficiency is improved. It is about twice as fast as that of the ATRACT. Therefore, this work provides a simple and efficient approach for the approximate reconstruction from truncated projections in the circular cone-beam CT.
基金Supported by National Key Research and Development Program of China(2016YFF0201005)。
文摘Based on anisotropic total variation regularization(ATVR), a nonnegativity and support constraints recursive inverse filtering(NAS-RIF) blind restoration method is proposed to enhance the quality of optical coherence tomography(OCT) image. First, ATVR is introduced into the cost function of NAS-RIF to improve the noise robustness and retain the details in the image.Since the split Bregman iterative is used to optimize the ATVR based cost function, the ATVR based NAS-RIF blind restoration method is then constructed. Furthermore, combined with the geometric nonlinear diffusion filter and the Poisson-distribution-based minimum error thresholding, the ATVR based NAS-RIF blind restoration method is used to realize the blind OCT image restoration. The experimental results demonstrate that the ATVR based NAS-RIF blind restoration method can successfully retain the details in the OCT images. In addition, the signal-to-noise ratio of the blind restored OCT images can be improved, along with the noise robustness.
文摘Radio astronomy radio telescope plays the role of a linear operator, affecting the function that describes the object of research, formation of image of a monitored object. This paper presents methods for reconstruction and correction of solar radio images using the algorithm of rejections, the updated Weiner-filter, and the method CLEAN designed by Hegbomom (Pseudonym, 2009) for point sources. It is the process of numerical convolution in signal handling, an algorithm for separating weak-contrast formations on the solar which represents most points of the actual limb by using the ellipse equation. Consequently, the filling algorithm is applied by moving from the center to the ellipse points and filling each point by solar image data. Finally, a linear limb-darkening expression is used to remove the limb darkening. Different examples of the intermediate and final results are presented in addition to the developed algorithm.
基金Project(41374118)supported by the National Natural Science Foundation,ChinaProject(20120162110015)supported by Research Fund for the Doctoral Program of Higher Education,China+3 种基金Project(2015M580700)supported by the China Postdoctoral Science Foundation,ChinaProject(2016JJ3086)supported by the Hunan Provincial Natural Science Foundation,ChinaProject(2015JC3067)supported by the Hunan Provincial Science and Technology Program,ChinaProject(15B138)supported by the Scientific Research Fund of Hunan Provincial Education Department,China
文摘To improve the global search ability and imaging quality of electrical resistivity imaging(ERI) inversion, a two-stage learning ICPSO algorithm of radial basis function neural network(RBFNN) based on information criterion(IC) and particle swarm optimization(PSO) is presented. In the proposed method, IC is applied to obtain the hidden layer structure by calculating the optimal IC value automatically and PSO algorithm is used to optimize the centers and widths of the radial basis functions in the hidden layer. Meanwhile, impacts of different information criteria to the inversion results are compared, and an implementation of the proposed ICPSO algorithm is given. The optimized neural network has one hidden layer with 261 nodes selected by AKAIKE's information criterion(AIC) and it is trained on 32 data sets and tested on another 8 synthetic data sets. Two complex synthetic examples are used to verify the feasibility and effectiveness of the proposed method with two learning stages. The results show that the proposed method has better performance and higher imaging quality than three-layer and four-layer back propagation neural networks(BPNNs) and traditional least square(LS) inversion.
文摘We survey fundamental concepts for inverse programming and then present the Universal Resolving Algorithm, an algorithm for inverse computation in a first order, functional programming language. We discuss the key concepts of the algorithm, including a three step approach based on the notion of a perfect process tree, and demonstrate our implementation with several examples of inverse computation.
基金supported by National Natural Science Foundation of China(41974166)Natural Science Foundation of Hebei Province(D2019403082,D2021403010)+1 种基金Hebei Province“three-threethree talent project”(A202005009)Funding for the Science and Technology Innovation Team Project of Hebei GEO University(KJCXTD202106)
文摘Reflection imaging results generally reveal large-scale continuous geological information,and it is difficult to identify small-scale geological bodies such as breakpoints,pinch points,small fault blocks,caves,and fractures,etc.Diffraction imaging is an important method to identify small-scale geological bodies and it has higher resolution than reflection imaging.In the common-offset domain,reflections are mostly expressed as smooth linear events,whereas diffractions are characterized by hyperbolic events.This paper proposes a diffraction extraction method based on double sparse transforms.The linear events can be sparsely expressed by the high-resolution linear Radon transform,and the curved events can be sparsely expressed by the Curvelet transform.A sparse inversion model is built and the alternating direction method is used to solve the inversion model.Simulation data and field data experimental results proved that the diffractions extraction method based on double sparse transforms can effectively improve the imaging quality of faults and other small-scale geological bodies.
文摘In this paper,we first construct generalized q^2-cosine,q^2-sine and q^2-exponential functions.We then use q^2-exponential function in order to define and investigate a q^2-Fourier transform.We establish q-analogues of inversion and Plancherel theorems.
文摘the technique of image processing and analysis of gravity and magnetic data is one of themost effective ways to extract geological information from gravity and msanetic data. The presentpaper investigates, from an angle of generalized joint inversion, thc methods and procedures ofcomprehensive processing of multi-source geological image , and a specific example in Huai Nan coalfield is given here as well.
文摘The adaptive reconstruction for the lost information of the rectangular image area is very important for the robust transmission and restoration of the image. In this paper, a new reconstruction method based on the Discrete Cosine Transform (DCT) domain has been put forward. According to the low pass character of the human visual system and the energy distribution of the DCT coefficients on the rectangular boundary, the DCT coefficients of the rectangular image area are adaptively selected and recovered. After the Inverse Discrete Cosine Transform (IDCT), the lost information of the rectangular image area can be reconstructed. The experiments have demonstrated that the subjective and objective qualities of the reconstructed images are enhanced greatly than before.