Optical-resolution photoacoustic microscopy is a novel imaging technique that combines the advantages of optical and ultrasound imaging,enabling high-resolution visualization of biological tissues at the micrometer sc...Optical-resolution photoacoustic microscopy is a novel imaging technique that combines the advantages of optical and ultrasound imaging,enabling high-resolution visualization of biological tissues at the micrometer scale.However,the divergence of the excited Gaussian beam limits the depth-of-field of the system to less than 100μm,which hinders accurate three-dimensional imaging of living tissues and restrictsits applicability in biological research.Therefore,there is an urgent need for an effective method to enhance the depth-of-field without altering the hardware configuration.This paper presents a photoacoustic microscopy depth-of-field extension method and system based on three-dimensional continuity and sparsity deconvolution.This method utilizes a depth-varying point spread function and incorporates continuity and sparsity con-straints into the deconvolution process to mitigate the effect of background noise,enhancing the stability and accuracy of the depth-of-field extension.Experimental results using tungsten wire phantoms suggest that the depth-of-field of system can be extended to 650 pm,which is 7.2 times greater than conventional system,while improving the resolution of the defocused region by an average factor of 3.5.Furthermore,experiments on zebrafish and nude mouse ears with irregular topologies demonstrate that the proposed method successfully overcomes image blurring and the loss of structural information due to limited depth-of-field.All the results suggest that the system with higher lateral resolution and enhanced depth-of-field has significant potential for a wide range of practical biomedical applications.展开更多
Traditional deconvolution methods based on single-channel inversion do not consider the spatial structural relation between channels,and hence,they yield high-resolution results with the existing transverse inconsiste...Traditional deconvolution methods based on single-channel inversion do not consider the spatial structural relation between channels,and hence,they yield high-resolution results with the existing transverse inconsistency or discontinuity.Therefore,in this study,the local dip angle was used to obtain the structural information and construct the spatial structurally constraint operator.This operator is then introduced into multichannel deconvolution as a regularization operator to improve the resolution and maintain the transverse continuity of seismic data.Model tests and actual seismic data processing have demonstrated the effectiveness and practicability of this method.展开更多
The conventional nonstationary convolutional model assumes that the seismic signal is recorded at normal incidence. Raw shot gathers are far from this assumption because of the effects of offsets. Because of such prob...The conventional nonstationary convolutional model assumes that the seismic signal is recorded at normal incidence. Raw shot gathers are far from this assumption because of the effects of offsets. Because of such problems, we propose a novel prestack nonstationary deconvolution approach. We introduce the radial trace （RT） transform to the nonstationary deconvolution, we estimate the nonstationary deconvolution factor with hyperbolic smoothing based on variable-step sampling （VSS） in the RT domain, and we obtain the high-resolution prestack nonstationary deconvolution data. The RT transform maps the shot record from the offset and traveltime coordinates to those of apparent velocity and traveltime. The ray paths of the traces in the RT better satisfy the assumptions of the convolutional model. The proposed method combines the advantages of stationary deconvolution and inverse Q filtering, without prior information for Q. The nonstationary deconvolution in the RT domain is more suitable than that in the space-time （XT） domain for prestack data because it is the generalized extension of normal incidence. Tests with synthetic and real data demonstrate that the proposed method is more effective in compensating for large-offset and deep data.展开更多
The Robinson convolution model is mainly restricted by three inappropriate assumptions, i.e., statistically white reflectivity, minimum-phase wavelet, and stationarity. Modern reflectivity inversion methods（e.g., spa...The Robinson convolution model is mainly restricted by three inappropriate assumptions, i.e., statistically white reflectivity, minimum-phase wavelet, and stationarity. Modern reflectivity inversion methods（e.g., sparsity-constrained deconvolution） generally attempt to suppress the problems associated with the first two assumptions but often ignore that seismic traces are nonstationary signals, which undermines the basic assumption of unchanging wavelet in reflectivity inversion. Through tests on reflectivity series, we confirm the effects of nonstationarity on reflectivity estimation and the loss of significant information, especially in deep layers. To overcome the problems caused by nonstationarity, we propose a nonstationary convolutional model, and then use the attenuation curve in log spectra to detect and correct the influences of nonstationarity. We use Gabor deconvolution to handle nonstationarity and sparsity-constrained deconvolution to separating reflectivity and wavelet. The combination of the two deconvolution methods effectively handles nonstationarity and greatly reduces the problems associated with the unreasonable assumptions regarding reflectivity and wavelet. Using marine seismic data, we show that correcting nonstationarity helps recover subtle reflectivity information and enhances the characterization of details with respect to the geological record.展开更多
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.展开更多
Sparsity constrained deconvolution can improve the resolution of band-limited seismic data compared to conventional deconvolution. However, such deconvolution methods result in nonunique solutions and suppress weak re...Sparsity constrained deconvolution can improve the resolution of band-limited seismic data compared to conventional deconvolution. However, such deconvolution methods result in nonunique solutions and suppress weak reflections. The Cauchy function, modified Cauchy function, and Huber function are commonly used constraint criteria in sparse deconvolution. We used numerical experiments to analyze the ability of sparsity constrained deconvolution to restore reflectivity sequences and protect weak reflections under different constraint criteria. The experimental results demonstrate that the performance of sparsity constrained deconvolution depends on the agreement between the constraint criteria and the probability distribution of the reflectivity sequences; furthermore, the modified Cauchy- constrained criterion protects the weak reflections better than the other criteria. Based on the model experiments, the probability distribution of the reflectivity sequences of carbonate and clastic formations is statistically analyzed by using well-logging data and then the modified Cauchy-constrained deconvolution is applied to real seismic data much improving the resolution.展开更多
The predictive deconvolution algorithm （PD）, which is based on second-order statistics, assumes that the primaries and the multiples are implicitly orthogonal. However, the seismic data usually do not satisfy this a...The predictive deconvolution algorithm （PD）, which is based on second-order statistics, assumes that the primaries and the multiples are implicitly orthogonal. However, the seismic data usually do not satisfy this assumption in practice. Since the seismic data （primaries and multiples） have a non-Gaussian distribution, in this paper we present an improved predictive deconvolution algorithm （IPD） by maximizing the non-Gaussianity of the recovered primaries. Applications of the IPD method on synthetic and real seismic datasets show that the proposed method obtains promising results.展开更多
A deconvolution algorithm is proposed to account for the distortions of impulse shape introduced by propagation process. By finding the best correlation of the received waveform with the multiple templates, the number...A deconvolution algorithm is proposed to account for the distortions of impulse shape introduced by propagation process. By finding the best correlation of the received waveform with the multiple templates, the number of multipath components is reduced as the result of eliminating the ＂phantom paths＂, and the captured energy increases. Moreover, it needs only a single reference measurement in real measurement environment （do not need the anechoic chamber）, which by far simplifies the templates acquiring procedure.展开更多
In wave-equation migration and demigration,the cross-correlation imaging/forwarding step implicitly injects an additional copy of the source wavelet,so that the amplitude spectrum of the wavelet is applied redundantly...In wave-equation migration and demigration,the cross-correlation imaging/forwarding step implicitly injects an additional copy of the source wavelet,so that the amplitude spectrum of the wavelet is applied redundantly(effectively imposing a wavelet-spectrum weighting,often akin to an amplitude-squared bias).This redundancy degrades structural fidelity and amplitude balance yet is frequently overlooked.We(i)formalize the mechanism by which cross-correlation duplicates the source-wavelet amplitude effect in both migration and demigration,and(ii)introduce a source-equalized operator that removes the redundancy by deconvolving(or dividing by)the wavelet amplitude spectrum in the imaging condition and its demigration counterpart,while leaving phase/kinematics intact.Using a band-limited Ricker wavelet on a two-layer model and on Marmousi,we show that,if unmanaged,the redundant wavelet spectrum broadens main lobes,introduces ringing,and suppresses vertical resolution in migrated images,and inflates spectrum mismatches between demigrated and observed data even when peak times agree.With our correction,images recover observed-data-consistent bandwidth and sharpened interfaces,and demigrated data also exhibit improved spectrum conformity and reduced amplitude misfit.The results clarify when source amplitudes matter,why cross-correlation makes them redundantly matter,and how a lightweight spectral correction restores physically meaningful amplitude behavior in wave-equation migration/demigration.展开更多
The Paper Presents a deconvolution method based on the double passthrough the human eye By using an optical-digital-computer. S system, the methodsolves the point-spread function and the optical transfer function of t...The Paper Presents a deconvolution method based on the double passthrough the human eye By using an optical-digital-computer. S system, the methodsolves the point-spread function and the optical transfer function of the optical systemof eye from the reflective image of retina, reconstructs the retinal image.As the double pass feature says, in the Positive direction and negative direction ofthe light propagation, the amplitude spread function is the same along the same media ofeye. If He-Ne laser beam is used as a short-time point impulse source, the incident lightcan be referred to as coherent light. Because the surface of retina is not fine smooth,usually the reflection is in company with the scattering, the reflective light beam is considered as incoherent light. After many times of flashings, samplings and averagings,the reflective ratio of retina will approximate to a constant and be neglected. Therefore,the conclusion is: the light intensity distribution outside of eye is the auto-convolution ofthe light intensity distribution of the renal image, i. e. the auto-convolution of thePOint-spread function. By the convolution theorem in the spatial frequency domain,two-dimensional Fourier transform of the point-spread runction, i. e. the optical transfer runction (OTF), can be calculated. Sequentially, by inverse Fourier transform, thePOint-spead function (PSF) is obtained. Thus the de-auto-convolution is accomplished.If a real picture replaces the point light source, the picture image on the fovea near thelight axis and then is reflected incoherently, the intensity distribution of the receivedimage outside of eye is equal to the crossconvolution between the intensity distribution ofretina and the point-spread function. By doing deconvolution once more. i. e. once decross-convolution, the retinal image of the picture can be reconstructed. The results ofthe experiments show that the hybrid system is advanced at the objectivity, auto-calibration and dynamic recording.展开更多
基金supported by the National Key R&D Program of China[Grant No.2022YFC2402400]the National Natural Science Foundation of China[Grant No.62275062]+2 种基金Project of Shandong Innovation and Startup Community of High-end Medical Apparatus and Instruments[Grant Nos.2023-SGTTXM-002 and 2024-SGTTXM-005]the Shandong Province Technology Innovation Guidance Plan(Central Leading Local Science and Technology Development Fund)[Grant No.YDZX2023115]the Taishan Scholar Special Funding Project of Shandong Province,and the Shandong Laboratory of Advanced Biomaterials and Medical Devices in Weihai[Grant No.ZL202402].
文摘Optical-resolution photoacoustic microscopy is a novel imaging technique that combines the advantages of optical and ultrasound imaging,enabling high-resolution visualization of biological tissues at the micrometer scale.However,the divergence of the excited Gaussian beam limits the depth-of-field of the system to less than 100μm,which hinders accurate three-dimensional imaging of living tissues and restrictsits applicability in biological research.Therefore,there is an urgent need for an effective method to enhance the depth-of-field without altering the hardware configuration.This paper presents a photoacoustic microscopy depth-of-field extension method and system based on three-dimensional continuity and sparsity deconvolution.This method utilizes a depth-varying point spread function and incorporates continuity and sparsity con-straints into the deconvolution process to mitigate the effect of background noise,enhancing the stability and accuracy of the depth-of-field extension.Experimental results using tungsten wire phantoms suggest that the depth-of-field of system can be extended to 650 pm,which is 7.2 times greater than conventional system,while improving the resolution of the defocused region by an average factor of 3.5.Furthermore,experiments on zebrafish and nude mouse ears with irregular topologies demonstrate that the proposed method successfully overcomes image blurring and the loss of structural information due to limited depth-of-field.All the results suggest that the system with higher lateral resolution and enhanced depth-of-field has significant potential for a wide range of practical biomedical applications.
基金supported by the basic and forward-looking project(No.2023YQX302)。
文摘Traditional deconvolution methods based on single-channel inversion do not consider the spatial structural relation between channels,and hence,they yield high-resolution results with the existing transverse inconsistency or discontinuity.Therefore,in this study,the local dip angle was used to obtain the structural information and construct the spatial structurally constraint operator.This operator is then introduced into multichannel deconvolution as a regularization operator to improve the resolution and maintain the transverse continuity of seismic data.Model tests and actual seismic data processing have demonstrated the effectiveness and practicability of this method.
基金financially supported by the National Science and Technology Major Project of China(No.2011ZX05023-005-005)the National Natural Science Foundation of China(No.41274137)
文摘The conventional nonstationary convolutional model assumes that the seismic signal is recorded at normal incidence. Raw shot gathers are far from this assumption because of the effects of offsets. Because of such problems, we propose a novel prestack nonstationary deconvolution approach. We introduce the radial trace （RT） transform to the nonstationary deconvolution, we estimate the nonstationary deconvolution factor with hyperbolic smoothing based on variable-step sampling （VSS） in the RT domain, and we obtain the high-resolution prestack nonstationary deconvolution data. The RT transform maps the shot record from the offset and traveltime coordinates to those of apparent velocity and traveltime. The ray paths of the traces in the RT better satisfy the assumptions of the convolutional model. The proposed method combines the advantages of stationary deconvolution and inverse Q filtering, without prior information for Q. The nonstationary deconvolution in the RT domain is more suitable than that in the space-time （XT） domain for prestack data because it is the generalized extension of normal incidence. Tests with synthetic and real data demonstrate that the proposed method is more effective in compensating for large-offset and deep data.
基金funded by the National Basic Research Program of China(973 Program)(Grant No.2011CB201100)Major Program of the National Natural Science Foundation of China(Grant No.2011ZX05004003)
文摘The Robinson convolution model is mainly restricted by three inappropriate assumptions, i.e., statistically white reflectivity, minimum-phase wavelet, and stationarity. Modern reflectivity inversion methods（e.g., sparsity-constrained deconvolution） generally attempt to suppress the problems associated with the first two assumptions but often ignore that seismic traces are nonstationary signals, which undermines the basic assumption of unchanging wavelet in reflectivity inversion. Through tests on reflectivity series, we confirm the effects of nonstationarity on reflectivity estimation and the loss of significant information, especially in deep layers. To overcome the problems caused by nonstationarity, we propose a nonstationary convolutional model, and then use the attenuation curve in log spectra to detect and correct the influences of nonstationarity. We use Gabor deconvolution to handle nonstationarity and sparsity-constrained deconvolution to separating reflectivity and wavelet. The combination of the two deconvolution methods effectively handles nonstationarity and greatly reduces the problems associated with the unreasonable assumptions regarding reflectivity and wavelet. Using marine seismic data, we show that correcting nonstationarity helps recover subtle reflectivity information and enhances the characterization of details with respect to the geological record.
基金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.
基金supported by the Major Basic Research Development Program of China (973 Program)(No.2013CB228606)the National Science foundation of China (No.41174117)+1 种基金the National Major Science-Technology Project (No.2011ZX05031-001)Innovation Fund of PetroChina (No.2010D-5006-0301)
文摘Sparsity constrained deconvolution can improve the resolution of band-limited seismic data compared to conventional deconvolution. However, such deconvolution methods result in nonunique solutions and suppress weak reflections. The Cauchy function, modified Cauchy function, and Huber function are commonly used constraint criteria in sparse deconvolution. We used numerical experiments to analyze the ability of sparsity constrained deconvolution to restore reflectivity sequences and protect weak reflections under different constraint criteria. The experimental results demonstrate that the performance of sparsity constrained deconvolution depends on the agreement between the constraint criteria and the probability distribution of the reflectivity sequences; furthermore, the modified Cauchy- constrained criterion protects the weak reflections better than the other criteria. Based on the model experiments, the probability distribution of the reflectivity sequences of carbonate and clastic formations is statistically analyzed by using well-logging data and then the modified Cauchy-constrained deconvolution is applied to real seismic data much improving the resolution.
基金National 863 Foundation of China(No.2006AA09A102-10)National Natural Science Foundation of China(No.40874056)NCET Fund
文摘The predictive deconvolution algorithm （PD）, which is based on second-order statistics, assumes that the primaries and the multiples are implicitly orthogonal. However, the seismic data usually do not satisfy this assumption in practice. Since the seismic data （primaries and multiples） have a non-Gaussian distribution, in this paper we present an improved predictive deconvolution algorithm （IPD） by maximizing the non-Gaussianity of the recovered primaries. Applications of the IPD method on synthetic and real seismic datasets show that the proposed method obtains promising results.
基金the Key Program of the National Natural Science Foundation of China (60432040)the China Postdoctors Science Foundation (20060390792).
文摘A deconvolution algorithm is proposed to account for the distortions of impulse shape introduced by propagation process. By finding the best correlation of the received waveform with the multiple templates, the number of multipath components is reduced as the result of eliminating the ＂phantom paths＂, and the captured energy increases. Moreover, it needs only a single reference measurement in real measurement environment （do not need the anechoic chamber）, which by far simplifies the templates acquiring procedure.
基金supported by the National Natural Science Foundation of China(42430303)Strategy Priority Research Program(Category B)of the Chinese Academy of Sciences(XDB0710000)+2 种基金National Natural Science Foundation of China(42288201)the National Key R&D Program of China(2023YFF0803203)the IGGCAS start-up funding(Grant No.E251510101).
文摘In wave-equation migration and demigration,the cross-correlation imaging/forwarding step implicitly injects an additional copy of the source wavelet,so that the amplitude spectrum of the wavelet is applied redundantly(effectively imposing a wavelet-spectrum weighting,often akin to an amplitude-squared bias).This redundancy degrades structural fidelity and amplitude balance yet is frequently overlooked.We(i)formalize the mechanism by which cross-correlation duplicates the source-wavelet amplitude effect in both migration and demigration,and(ii)introduce a source-equalized operator that removes the redundancy by deconvolving(or dividing by)the wavelet amplitude spectrum in the imaging condition and its demigration counterpart,while leaving phase/kinematics intact.Using a band-limited Ricker wavelet on a two-layer model and on Marmousi,we show that,if unmanaged,the redundant wavelet spectrum broadens main lobes,introduces ringing,and suppresses vertical resolution in migrated images,and inflates spectrum mismatches between demigrated and observed data even when peak times agree.With our correction,images recover observed-data-consistent bandwidth and sharpened interfaces,and demigrated data also exhibit improved spectrum conformity and reduced amplitude misfit.The results clarify when source amplitudes matter,why cross-correlation makes them redundantly matter,and how a lightweight spectral correction restores physically meaningful amplitude behavior in wave-equation migration/demigration.
文摘The Paper Presents a deconvolution method based on the double passthrough the human eye By using an optical-digital-computer. S system, the methodsolves the point-spread function and the optical transfer function of the optical systemof eye from the reflective image of retina, reconstructs the retinal image.As the double pass feature says, in the Positive direction and negative direction ofthe light propagation, the amplitude spread function is the same along the same media ofeye. If He-Ne laser beam is used as a short-time point impulse source, the incident lightcan be referred to as coherent light. Because the surface of retina is not fine smooth,usually the reflection is in company with the scattering, the reflective light beam is considered as incoherent light. After many times of flashings, samplings and averagings,the reflective ratio of retina will approximate to a constant and be neglected. Therefore,the conclusion is: the light intensity distribution outside of eye is the auto-convolution ofthe light intensity distribution of the renal image, i. e. the auto-convolution of thePOint-spread function. By the convolution theorem in the spatial frequency domain,two-dimensional Fourier transform of the point-spread runction, i. e. the optical transfer runction (OTF), can be calculated. Sequentially, by inverse Fourier transform, thePOint-spead function (PSF) is obtained. Thus the de-auto-convolution is accomplished.If a real picture replaces the point light source, the picture image on the fovea near thelight axis and then is reflected incoherently, the intensity distribution of the receivedimage outside of eye is equal to the crossconvolution between the intensity distribution ofretina and the point-spread function. By doing deconvolution once more. i. e. once decross-convolution, the retinal image of the picture can be reconstructed. The results ofthe experiments show that the hybrid system is advanced at the objectivity, auto-calibration and dynamic recording.
