In this study, we present a practical technique of transforming cross-hole EM data into the inter-well resistivity distribution. The a priori information constraint is incorporated into an iterative regularized invers...In this study, we present a practical technique of transforming cross-hole EM data into the inter-well resistivity distribution. The a priori information constraint is incorporated into an iterative regularized inversion procedure and a variable roughness is added into the inversion process. Finite element approximation based on a two and a half-dimensional （2.5D） model has been developed for the forward problem and the ＂pseudo-forward＂ problem needed for constructing the sensitivity matrix and synthetic data set. The regularized least-squares inversion scheme, constrained with the a priori information obtained from well logs, was adopted to reconstruct the inter-well resistivity profile from two synthetic electromagnetic data sets and field data acquired in the Gudao Oil Field, East China. The partial derivatives of the sensitivity matrix were computed by the adjoint equation based on the reciprocity principle. Inversion results of the synthetic and field data examples suggest that our method is robust and stable in the presence of random noise in the field data and can be used for cross-hole EM field data interpretation.展开更多
To solve the problem that metal artifacts severely damage the clarity of the organization structure in computed tomography（CT） images, a sinogram fusion-based metal artifact correction method is proposed. First, the...To solve the problem that metal artifacts severely damage the clarity of the organization structure in computed tomography（CT） images, a sinogram fusion-based metal artifact correction method is proposed. First, the metal image is segmented from the original CT image by the pre-set threshold. The original CT image and metal image are forward projected into the original projection sinogram and metal projection sinogram, respectively. The interpolation-based correction method and mean filter are used to correct the original CT image and preserve the edge of the corrected CT image, respectively. The filtered CT image is forward projected into the filtered image sinogram. According to the position of the metal sinogram in the original sinogram and filtered image sinogram, the corresponding sinograms PM^D （ in the original sinogram） and PM^C （ in the filtered image sinogram）can be acquired from the original sinogram and filtered image sinogram, respectively. Then, PM^D and PM^C are fused into the fused metal sinogram PM^F according to a certain proportion.The final sinogram can be acquired by fusing PM^F , PM^D and the original sinogram P^O. Finally, the final sinogram is reconstructed into the corrected CT image and metal information is compensated into the corrected CT image.Experiments on clinical images demonstrate that the proposed method can effectively reduce metal artifacts. A comparison with classical metal artifacts correction methods shows that the proposed metal artifacts correction method performs better in metal artifacts suppression and tissue feature preservation.展开更多
Linear scan computed tomography （LCT） is of great benefit to online industrial scanning and security inspection due to its characteristics of straight-line source trajectory and high scanning speed. However, in prac...Linear scan computed tomography （LCT） is of great benefit to online industrial scanning and security inspection due to its characteristics of straight-line source trajectory and high scanning speed. However, in practical applications of LCT, there are challenges to image reconstruction due to limited-angle and insufficient data. In this paper, a new reconstruction algorithm based on total-variation （TV） minimization is developed to reconstruct images from limited-angle and insufficient data in LCT. The main idea of our approach is to reformulate a TV problem as a linear equality constrained problem where the objective function is separable, and then minimize its augmented Lagrangian function by using alternating direction method （ADM） to solve subproblems. The proposed method is robust and efficient in the task of reconstruction by showing the convergence of ADM. The numerical simulations and real data reconstructions show that the proposed reconstruction method brings reasonable performance and outperforms some previous ones when applied to an LCT imaging problem.展开更多
In recent years,karst construction projects in the built-up area of Wuhan(capital of Hubei Province,China)are increasing,and the karst geological disasters have aroused social concerns.The actual engineering projects ...In recent years,karst construction projects in the built-up area of Wuhan(capital of Hubei Province,China)are increasing,and the karst geological disasters have aroused social concerns.The actual engineering projects usually use shallow geophysical exploration methods to explore karst.This paper uses Spatial Auto-Correlation Method(SPAC)and electromagnetic Computerized Tomography(CT)to detect karst in urban built-up areas.Depending on the different physical properties of rock and soil,the SPAC method can better reveal the interface between soil and rock strata and the interface between soil layers.The electromagnetic CT method can identify strata according to the apparent absorption coefficient,which can better reveal the interface between soil and rock,the interface between the more intact and weathered rock.The SPAC method is mainly qualitative to measure the low-speed area,namely,the wrong geological body i.e.,karst cave,but also can detect the fracture zone or filling mode of karst cave,and at the same time,cannot use exploration holes or logging observation.The electromagnetic CT method can accurately detect the location and scale of the karst caves and has a higher accuracy detecting karst bands.In addition,exploration holes or well logging observations are also expected to be conducted,and their detection effect is greatly affected by lithology.展开更多
Current output dose measurement in CT is based on (CTDI). The conventional methodology of CT dosimetric performance characterization is not appropriate to modern CT scanners with helical scanning modes, dose modulatio...Current output dose measurement in CT is based on (CTDI). The conventional methodology of CT dosimetric performance characterization is not appropriate to modern CT scanners with helical scanning modes, dose modulation, array detectors and multiple slice planes or cone-beam irradiation geometries. AAPM TG 111 report recognizes the shortfall of the CTDI methods and recommends a new technique that more accurately characterizes the dose profile from modern CT scanners, which utilizes a short conventional ion chamber rather than a pencil chamber. We developed and characterize a in-house phantom design using a three separate anatomical regions of clinical scan sequences (Head, chest and abdomen), and determined the equilibrium dose in our dose equilibrium phantom, measured if the attenuation of the beam is the equal to that of CTDI Perspex phantom and compare CTDI dose estimations using a standard pencil chamber to the dose equilibrium phantom measurements. This methodology allows measurements of the accumulated dose for any clinical scan length and allowing measurement of the equilibrium dose. Using the new methodology, we determined that the CTDI approach can underestimate the dose by 25% to 35% and all of our dose values from the water phantom and farmer chamber were independently verified with TLD measurements.展开更多
Various kinds of heritages, such as architectures, statues, grave posts and towers, are made by stone, and they are facing the crisis of weathering. Therefore, it is necessary to give appropriate treatments to keep th...Various kinds of heritages, such as architectures, statues, grave posts and towers, are made by stone, and they are facing the crisis of weathering. Therefore, it is necessary to give appropriate treatments to keep them in good conditions. Kyushu Research Institute for Cultural Properties Inc. and Kumamoto University introduced a new method, Aquo-Siloxane Method, in order to protect the stone heritages. In this study, preservation effect by Aquo-Siloxane method towards water permeation and material diffusion was verified. Here one-dimensional permeation and diffusion tests were conducted, and the intrinsic permeability and diffusion coefficient of rock samples with and without Aquo-Siloxane treatments were evaluated. As rock samples, 3 types of sandstones and concrete were applied. It was found that the permeability decreased to less than 1/10 to 1/100 of without Aquo-Siloxane treatment, and that the effect gradually developed during more than one year. One-dimensional diffusion tests were also conducted, and the diffusion process in rock samples are visualized by X-ray CT scanner system. It was confirmed internal structures of samples are clearly visualized, and that the diffusion process was also visualized as X-ray CT images. In order to extract the necessary information due to diffusion, image subtraction method was applied to image data. Then, by comparing obtained CT image data and numerical solutions, diffusion coefficients of rock samples with and without Aquo-Siloxane treatments were evaluated. As a result, diffusion coefficients also became smaller by applying Aquo-Siloxane treatments. It is revealed that material movement due to diffusion was also efficiently suppressed by applying Aquo-Siloxane method. Finally, Aquo-Siloxane method was applied to a stone heritage. It was found that no more chippings and cracks were observed and that the heritage has been kept in good condition for at least four years.展开更多
Purpose: To present an application of the anisotropic diffusion (AD) method to improve the accuracy of the functional images of perfusion parameters such as cerebral blood flow (CBF), cerebral blood volume (CBV) and m...Purpose: To present an application of the anisotropic diffusion (AD) method to improve the accuracy of the functional images of perfusion parameters such as cerebral blood flow (CBF), cerebral blood volume (CBV) and mean transit time (MTT) generated from cerebral CT perfusion studies using multi-detector row CT (MDCT). Materials and Methods: Continuous scans (1 sec/rotation ×60 sec) consisting of four 5-mm-thick contiguous slices were acquired after an intravenous injection of iodinated contrast material in 6 patients with cerebrovascular disease using an MDCT scanner with a tube voltage of 80 kVp and a tube current of 200 mA. New image data were generated by thinning out the above original images at an interval of 2 sec or 3 sec. The thinned-out images were then interpolated by linear interpolation to generate the same number of images as originally acquired. The CBF, CBV and MTT images were generated using deconvolution analysis based on singular value decomposition. Results: When using the AD method, the correlation coefficient between the MTT values obtained from the original and thinned-out images was significantly improved. Furthermore, the coefficients of variation of the CBF, CBV and MTT values in the white matter significantly decreased as compared to not using the AD method. Conclusion: Our results suggest that the AD method is useful for improving the accuracy of the functional images of perfusion parameters and for reducing radiation exposure in cerebral CT perfusion studies using MDCT.展开更多
Computed tomography(CT) blurring caused by point spread function leads to errors in quantification and visualization. In this paper, multichannel blind CT image restoration is proposed to overcome the effect of point ...Computed tomography(CT) blurring caused by point spread function leads to errors in quantification and visualization. In this paper, multichannel blind CT image restoration is proposed to overcome the effect of point spread function. The main advantage from multichannel blind CT image restoration is to exploit the diversity and redundancy of information in different acquisitions. The proposed approach is based on a variable splitting to obtain an equivalent constrained optimization formulation, which is addressed with the alternating direction method of multipliers and simply implemented in the Fourier domain. Numerical experiments illustrate that our method obtains a higher average gain value of at least 1.21 d B in terms of Q metric than the other methods, and it requires only 7 iterations of alternating minimization to obtain a fast convergence.展开更多
The construction of modern livable cities faces challenges in karst areas,including ground collapse and engineering problems.Wuhan,with a population of 13.74×10^(6) and approximately 1161 km^(2)of soluble rocks i...The construction of modern livable cities faces challenges in karst areas,including ground collapse and engineering problems.Wuhan,with a population of 13.74×10^(6) and approximately 1161 km^(2)of soluble rocks in the urban area of 8569.15 km^(2),predominantly consists of concealed karst areas where occasional ground collapse events occur,posing significant threats to underground engineering projects.To address these challenges,a comprehensive geological survey was conducted in Wuhan,focusing on major karstrelated issues.Geophysical methods offer advantages over drilling in detecting concealed karst areas due to their efficiency,non-destructiveness,and flexibility.This paper reviewed the karst geological characteristics in Wuhan and the geophysical exploration methods for karst,selected eight effective geophysical methods for field experimentation,evaluated their suitability,and proposed method combinations for different karst scenarios.The results show that different geophysical methods have varying applicability for karst detection in Wuhan,and combining multiple methods enhances detection effectiveness.The specific recommendations for method combinations provided in this study serve as a valuable reference for karst detection in Wuhan.展开更多
文摘In this study, we present a practical technique of transforming cross-hole EM data into the inter-well resistivity distribution. The a priori information constraint is incorporated into an iterative regularized inversion procedure and a variable roughness is added into the inversion process. Finite element approximation based on a two and a half-dimensional （2.5D） model has been developed for the forward problem and the ＂pseudo-forward＂ problem needed for constructing the sensitivity matrix and synthetic data set. The regularized least-squares inversion scheme, constrained with the a priori information obtained from well logs, was adopted to reconstruct the inter-well resistivity profile from two synthetic electromagnetic data sets and field data acquired in the Gudao Oil Field, East China. The partial derivatives of the sensitivity matrix were computed by the adjoint equation based on the reciprocity principle. Inversion results of the synthetic and field data examples suggest that our method is robust and stable in the presence of random noise in the field data and can be used for cross-hole EM field data interpretation.
基金Open Research Fund of the Key Laboratory of Computer Netw ork and Information Integration of Ministry of Education of Southeast University(No.K93-9-2014-10C)the Scientific Research Foundation of Education Department of Anhui Province(No.KJ2014A186,SK2015A433)the National Basic Research Program of China(973 Program)(No.2010CB732503)
文摘To solve the problem that metal artifacts severely damage the clarity of the organization structure in computed tomography（CT） images, a sinogram fusion-based metal artifact correction method is proposed. First, the metal image is segmented from the original CT image by the pre-set threshold. The original CT image and metal image are forward projected into the original projection sinogram and metal projection sinogram, respectively. The interpolation-based correction method and mean filter are used to correct the original CT image and preserve the edge of the corrected CT image, respectively. The filtered CT image is forward projected into the filtered image sinogram. According to the position of the metal sinogram in the original sinogram and filtered image sinogram, the corresponding sinograms PM^D （ in the original sinogram） and PM^C （ in the filtered image sinogram）can be acquired from the original sinogram and filtered image sinogram, respectively. Then, PM^D and PM^C are fused into the fused metal sinogram PM^F according to a certain proportion.The final sinogram can be acquired by fusing PM^F , PM^D and the original sinogram P^O. Finally, the final sinogram is reconstructed into the corrected CT image and metal information is compensated into the corrected CT image.Experiments on clinical images demonstrate that the proposed method can effectively reduce metal artifacts. A comparison with classical metal artifacts correction methods shows that the proposed metal artifacts correction method performs better in metal artifacts suppression and tissue feature preservation.
基金the National High Technology Research and Development Program of China(Grant No.2012AA011603)
文摘Linear scan computed tomography （LCT） is of great benefit to online industrial scanning and security inspection due to its characteristics of straight-line source trajectory and high scanning speed. However, in practical applications of LCT, there are challenges to image reconstruction due to limited-angle and insufficient data. In this paper, a new reconstruction algorithm based on total-variation （TV） minimization is developed to reconstruct images from limited-angle and insufficient data in LCT. The main idea of our approach is to reformulate a TV problem as a linear equality constrained problem where the objective function is separable, and then minimize its augmented Lagrangian function by using alternating direction method （ADM） to solve subproblems. The proposed method is robust and efficient in the task of reconstruction by showing the convergence of ADM. The numerical simulations and real data reconstructions show that the proposed reconstruction method brings reasonable performance and outperforms some previous ones when applied to an LCT imaging problem.
文摘In recent years,karst construction projects in the built-up area of Wuhan(capital of Hubei Province,China)are increasing,and the karst geological disasters have aroused social concerns.The actual engineering projects usually use shallow geophysical exploration methods to explore karst.This paper uses Spatial Auto-Correlation Method(SPAC)and electromagnetic Computerized Tomography(CT)to detect karst in urban built-up areas.Depending on the different physical properties of rock and soil,the SPAC method can better reveal the interface between soil and rock strata and the interface between soil layers.The electromagnetic CT method can identify strata according to the apparent absorption coefficient,which can better reveal the interface between soil and rock,the interface between the more intact and weathered rock.The SPAC method is mainly qualitative to measure the low-speed area,namely,the wrong geological body i.e.,karst cave,but also can detect the fracture zone or filling mode of karst cave,and at the same time,cannot use exploration holes or logging observation.The electromagnetic CT method can accurately detect the location and scale of the karst caves and has a higher accuracy detecting karst bands.In addition,exploration holes or well logging observations are also expected to be conducted,and their detection effect is greatly affected by lithology.
文摘Current output dose measurement in CT is based on (CTDI). The conventional methodology of CT dosimetric performance characterization is not appropriate to modern CT scanners with helical scanning modes, dose modulation, array detectors and multiple slice planes or cone-beam irradiation geometries. AAPM TG 111 report recognizes the shortfall of the CTDI methods and recommends a new technique that more accurately characterizes the dose profile from modern CT scanners, which utilizes a short conventional ion chamber rather than a pencil chamber. We developed and characterize a in-house phantom design using a three separate anatomical regions of clinical scan sequences (Head, chest and abdomen), and determined the equilibrium dose in our dose equilibrium phantom, measured if the attenuation of the beam is the equal to that of CTDI Perspex phantom and compare CTDI dose estimations using a standard pencil chamber to the dose equilibrium phantom measurements. This methodology allows measurements of the accumulated dose for any clinical scan length and allowing measurement of the equilibrium dose. Using the new methodology, we determined that the CTDI approach can underestimate the dose by 25% to 35% and all of our dose values from the water phantom and farmer chamber were independently verified with TLD measurements.
文摘Various kinds of heritages, such as architectures, statues, grave posts and towers, are made by stone, and they are facing the crisis of weathering. Therefore, it is necessary to give appropriate treatments to keep them in good conditions. Kyushu Research Institute for Cultural Properties Inc. and Kumamoto University introduced a new method, Aquo-Siloxane Method, in order to protect the stone heritages. In this study, preservation effect by Aquo-Siloxane method towards water permeation and material diffusion was verified. Here one-dimensional permeation and diffusion tests were conducted, and the intrinsic permeability and diffusion coefficient of rock samples with and without Aquo-Siloxane treatments were evaluated. As rock samples, 3 types of sandstones and concrete were applied. It was found that the permeability decreased to less than 1/10 to 1/100 of without Aquo-Siloxane treatment, and that the effect gradually developed during more than one year. One-dimensional diffusion tests were also conducted, and the diffusion process in rock samples are visualized by X-ray CT scanner system. It was confirmed internal structures of samples are clearly visualized, and that the diffusion process was also visualized as X-ray CT images. In order to extract the necessary information due to diffusion, image subtraction method was applied to image data. Then, by comparing obtained CT image data and numerical solutions, diffusion coefficients of rock samples with and without Aquo-Siloxane treatments were evaluated. As a result, diffusion coefficients also became smaller by applying Aquo-Siloxane treatments. It is revealed that material movement due to diffusion was also efficiently suppressed by applying Aquo-Siloxane method. Finally, Aquo-Siloxane method was applied to a stone heritage. It was found that no more chippings and cracks were observed and that the heritage has been kept in good condition for at least four years.
文摘Purpose: To present an application of the anisotropic diffusion (AD) method to improve the accuracy of the functional images of perfusion parameters such as cerebral blood flow (CBF), cerebral blood volume (CBV) and mean transit time (MTT) generated from cerebral CT perfusion studies using multi-detector row CT (MDCT). Materials and Methods: Continuous scans (1 sec/rotation ×60 sec) consisting of four 5-mm-thick contiguous slices were acquired after an intravenous injection of iodinated contrast material in 6 patients with cerebrovascular disease using an MDCT scanner with a tube voltage of 80 kVp and a tube current of 200 mA. New image data were generated by thinning out the above original images at an interval of 2 sec or 3 sec. The thinned-out images were then interpolated by linear interpolation to generate the same number of images as originally acquired. The CBF, CBV and MTT images were generated using deconvolution analysis based on singular value decomposition. Results: When using the AD method, the correlation coefficient between the MTT values obtained from the original and thinned-out images was significantly improved. Furthermore, the coefficients of variation of the CBF, CBV and MTT values in the white matter significantly decreased as compared to not using the AD method. Conclusion: Our results suggest that the AD method is useful for improving the accuracy of the functional images of perfusion parameters and for reducing radiation exposure in cerebral CT perfusion studies using MDCT.
基金Supported by the National Natural Science Foundaton of China(No.61340034)China Postdoctoral Science Foundation(No.2013M530873)the Research Program of Application Foundation and Advanced Technology of Tianjin(No.13JCYBJC15600)
文摘Computed tomography(CT) blurring caused by point spread function leads to errors in quantification and visualization. In this paper, multichannel blind CT image restoration is proposed to overcome the effect of point spread function. The main advantage from multichannel blind CT image restoration is to exploit the diversity and redundancy of information in different acquisitions. The proposed approach is based on a variable splitting to obtain an equivalent constrained optimization formulation, which is addressed with the alternating direction method of multipliers and simply implemented in the Fourier domain. Numerical experiments illustrate that our method obtains a higher average gain value of at least 1.21 d B in terms of Q metric than the other methods, and it requires only 7 iterations of alternating minimization to obtain a fast convergence.
基金jointly supported by the project of Chinese National Natural Science Foundation(42107485)National Key R&D Program(2020YFC1512400,2018YFC800804)China Geological Survey(DD20190282,DD20221734,and DD20230323)。
文摘The construction of modern livable cities faces challenges in karst areas,including ground collapse and engineering problems.Wuhan,with a population of 13.74×10^(6) and approximately 1161 km^(2)of soluble rocks in the urban area of 8569.15 km^(2),predominantly consists of concealed karst areas where occasional ground collapse events occur,posing significant threats to underground engineering projects.To address these challenges,a comprehensive geological survey was conducted in Wuhan,focusing on major karstrelated issues.Geophysical methods offer advantages over drilling in detecting concealed karst areas due to their efficiency,non-destructiveness,and flexibility.This paper reviewed the karst geological characteristics in Wuhan and the geophysical exploration methods for karst,selected eight effective geophysical methods for field experimentation,evaluated their suitability,and proposed method combinations for different karst scenarios.The results show that different geophysical methods have varying applicability for karst detection in Wuhan,and combining multiple methods enhances detection effectiveness.The specific recommendations for method combinations provided in this study serve as a valuable reference for karst detection in Wuhan.
