A small problem about soil particle regularization and contacts but essential to geotechnical engineering was studied.The soils sourced from Guangzhou and Xiamen were sieved into five different particle scale ranges(d...A small problem about soil particle regularization and contacts but essential to geotechnical engineering was studied.The soils sourced from Guangzhou and Xiamen were sieved into five different particle scale ranges(d<0.075 mm,0.075 mm≤d<0.1 mm,0.1 mm≤d<0.2 mm,0.2 mm≤d<0.5 mm and 0.5 mm≤d<1.0 mm)to study the structures and particle contacts of granite residual soil.The X-ray micro computed tomography method was used to reconstruct the microstructure of granite residual soil.The particle was identified and regularized using principal component analysis(PCA).The particle contacts and geometrical characteristics in 3D space were analyzed and summarized using statistical analyses.The results demonstrate that the main types of contact among the particles are face-face,face-angle,face-edge,edge-edge,edge-angle and angle-angle contacts for particle sizes less than 0.2 mm.When the particle sizes are greater than 0.2 mm,the contacts are effectively summarized as face-face,face-angle,face-edge,edge-edge,edge-angle,angle-angle,sphere-sphere,sphere-face,sphere-edge and sphere-angle contacts.The differences in porosity among the original sample,reconstructed sample and regularized sample are closely related to the water-swelling and water-disintegrable characteristics of granite residual soil.展开更多
The geometric and spatial characteristics of pore structures determine the permeability and water retention of soils, which have important effects on soil functional diversity and ecological restoration. Until recentl...The geometric and spatial characteristics of pore structures determine the permeability and water retention of soils, which have important effects on soil functional diversity and ecological restoration. Until recently, there have not been tools and methods to visually and quantitatively describe the characteristics of soil pores. To solve this problem, this research reconstructs the geometry and spatial distribution of soil pores by the marching cubes method, texture mapping method and the ray casting method widely used in literature. The objectives were to explore an optimal method for three-dimensional visualization of soil pore structure by comparing the robustness of the three methods on soil CT images with single pore structure and porosity ranging from low (2–5%) to high (12–18%), and to evaluate the reconstruction performance of the three methods with different geometric features. The results demonstrate that there are aliases (jagged edges) and deficiency at the boundaries of the model reconstructed by the marching cubes method and pore volumes are smaller than the ground truth, whereas the results of the texture mapping method lack the details of pore structures. For all the soil images, the ray casting method is preferable since it better preserves the pore characteristics of the ground truth. Furthermore, the ray casting method produced the best soil pore model with higher rendering speed and lower memory consumption. Therefore, the ray casting method provides a more advanced method for visualization of pore structures and provides an optional technique for the study of the transport of moisture and the exchange of air in soil.展开更多
The relationship between noise variance and spatial resolution in grating-based x-ray phase computed tomography（PCT） imaging is investigated with reverse projection extraction method, and the noise variances of the ...The relationship between noise variance and spatial resolution in grating-based x-ray phase computed tomography（PCT） imaging is investigated with reverse projection extraction method, and the noise variances of the reconstructed absorption coefficient and refractive index decrement are compared. For the differential phase contrast method, the noise variance in the differential projection images follows the same inverse-square law with spatial resolution as in conventional absorption-based x-ray imaging projections. However, both theoretical analysis and simulations demonstrate that in PCT the noise variance of the reconstructed refractive index decrement scales with spatial resolution follows an inverse linear relationship at fixed slice thickness, while the noise variance of the reconstructed absorption coefficient conforms with the inverse cubic law. The results indicate that, for the same noise variance level, PCT imaging may enable higher spatial resolution than conventional absorption computed tomography（ACT）, while ACT benefits more from degraded spatial resolution. This could be a useful guidance in imaging the inner structure of the sample in higher spatial resolution.展开更多
In conventional computed tomography （CT） reconstruction based on fixed voltage, the projective data often ap- pear overexposed or underexposed, as a result, the reconstructive results are poor. To solve this problem...In conventional computed tomography （CT） reconstruction based on fixed voltage, the projective data often ap- pear overexposed or underexposed, as a result, the reconstructive results are poor. To solve this problem, variable voltage CT reconstruction has been proposed. The effective projective sequences of a structural component are obtained through the variable voltage. The total variation is adjusted and minimized to optimize the reconstructive results on the basis of iterative image using algebraic reconstruction technique （ART）. In the process of reconstruction, the reconstructive image of low voltage is used as an initial value of the effective proiective reconstruction of the adjacent high voltage, and so on until to the highest voltage according to the gray weighted algorithm. Thereby the complete structural information is reconstructed. Simulation results show that the proposed algorithm can completely reflect the information of a complicated structural com- ponent, and the pixel values are more stable than those of the conventional.展开更多
X-ray grating interferometer has attracted widely attention in the past years due to its capability in achieving x-ray phase contrast imaging with low brilliance source. However, the widely used phase stepping informa...X-ray grating interferometer has attracted widely attention in the past years due to its capability in achieving x-ray phase contrast imaging with low brilliance source. However, the widely used phase stepping information extraction method reduces system stability and prolongs data acquisition time by several times compared with conventional x-ray absorption- based imaging. The mechanical stepping can be avoided by using a staggered grating, but at the cost of low vertical spatial resolution. In this paper, employing a modified staggered grating and the angular signal radiography, we proposed a single-shot grating-based x-ray differential phase contrast imaging with decent vertical spatial resolution. The theoretical framework was deduced and proved by numerical experiments. Absorption, phase, and scattering computed tomography can be performed without phase stepping. Therefore, we believe this fast and highly stable imaging method with decent resolution would be widely applied in x-ray grating-based phase contrast imaging.展开更多
We present a new iterative reconstruction algorithm to improve the algebraic reconstruction technique （ART） for the Single-Photon Emission Computed Tomography. Our method is a generalization of the Kaczmarz iterativ...We present a new iterative reconstruction algorithm to improve the algebraic reconstruction technique （ART） for the Single-Photon Emission Computed Tomography. Our method is a generalization of the Kaczmarz iterative algorithm for solving linear systems of equations and introduces exact and implicit attenuation correction derived from the attenuated Radon transform operator at each step of the algorithm. The performances of the presented algorithm have been tested upon various numerical experiments in presence of both strongly non-uniform attenuation and incomplete measurements data. We also tested the ability of our algorithm to handle moderate noisy data. Simulation studies demonstrate that the proposed method has a significant improvement in the quality of reconstructed images over ART. Moreover, convergence speed was improved and stability was established, facing noisy data, once we incorporate filtration procedure in our algorithm.展开更多
基金Projects(41572277,41877229) supported by the National Natural Science Foundation of ChinaProject(2015A030313118) supported by the Natural Science Foundation of Guangdong Province,ChinaProject(201607010023) supported by the Science and Technology Program of Guangzhou,China
文摘A small problem about soil particle regularization and contacts but essential to geotechnical engineering was studied.The soils sourced from Guangzhou and Xiamen were sieved into five different particle scale ranges(d<0.075 mm,0.075 mm≤d<0.1 mm,0.1 mm≤d<0.2 mm,0.2 mm≤d<0.5 mm and 0.5 mm≤d<1.0 mm)to study the structures and particle contacts of granite residual soil.The X-ray micro computed tomography method was used to reconstruct the microstructure of granite residual soil.The particle was identified and regularized using principal component analysis(PCA).The particle contacts and geometrical characteristics in 3D space were analyzed and summarized using statistical analyses.The results demonstrate that the main types of contact among the particles are face-face,face-angle,face-edge,edge-edge,edge-angle and angle-angle contacts for particle sizes less than 0.2 mm.When the particle sizes are greater than 0.2 mm,the contacts are effectively summarized as face-face,face-angle,face-edge,edge-edge,edge-angle,angle-angle,sphere-sphere,sphere-face,sphere-edge and sphere-angle contacts.The differences in porosity among the original sample,reconstructed sample and regularized sample are closely related to the water-swelling and water-disintegrable characteristics of granite residual soil.
基金supported by the National Natural Science Foundation Project(41501283)Beijing Science and Technology Plan Project(Z161100000916012)+2 种基金the National Key Research and Development Program(2017YFD0600901)Special Fund for Beijing Common Construction Projectthe Fundamental Research Funds for the Central Universities(2015ZCQ-GX-04)
文摘The geometric and spatial characteristics of pore structures determine the permeability and water retention of soils, which have important effects on soil functional diversity and ecological restoration. Until recently, there have not been tools and methods to visually and quantitatively describe the characteristics of soil pores. To solve this problem, this research reconstructs the geometry and spatial distribution of soil pores by the marching cubes method, texture mapping method and the ray casting method widely used in literature. The objectives were to explore an optimal method for three-dimensional visualization of soil pore structure by comparing the robustness of the three methods on soil CT images with single pore structure and porosity ranging from low (2–5%) to high (12–18%), and to evaluate the reconstruction performance of the three methods with different geometric features. The results demonstrate that there are aliases (jagged edges) and deficiency at the boundaries of the model reconstructed by the marching cubes method and pore volumes are smaller than the ground truth, whereas the results of the texture mapping method lack the details of pore structures. For all the soil images, the ray casting method is preferable since it better preserves the pore characteristics of the ground truth. Furthermore, the ray casting method produced the best soil pore model with higher rendering speed and lower memory consumption. Therefore, the ray casting method provides a more advanced method for visualization of pore structures and provides an optional technique for the study of the transport of moisture and the exchange of air in soil.
基金Project supported by the National Basic Research Program of China(Grant No.2012CB825800)the Science Fund for Creative Research Groups,the Knowledge Innovation Program of the Chinese Academy of Sciences(Grant Nos.KJCX2-YW-N42 and Y4545320Y2)the National Natural Science Foundation of China(Grant Nos.11475170,11205157,11305173,11205189,11375225,11321503,11179004,and U1332109)
文摘The relationship between noise variance and spatial resolution in grating-based x-ray phase computed tomography（PCT） imaging is investigated with reverse projection extraction method, and the noise variances of the reconstructed absorption coefficient and refractive index decrement are compared. For the differential phase contrast method, the noise variance in the differential projection images follows the same inverse-square law with spatial resolution as in conventional absorption-based x-ray imaging projections. However, both theoretical analysis and simulations demonstrate that in PCT the noise variance of the reconstructed refractive index decrement scales with spatial resolution follows an inverse linear relationship at fixed slice thickness, while the noise variance of the reconstructed absorption coefficient conforms with the inverse cubic law. The results indicate that, for the same noise variance level, PCT imaging may enable higher spatial resolution than conventional absorption computed tomography（ACT）, while ACT benefits more from degraded spatial resolution. This could be a useful guidance in imaging the inner structure of the sample in higher spatial resolution.
文摘In conventional computed tomography （CT） reconstruction based on fixed voltage, the projective data often ap- pear overexposed or underexposed, as a result, the reconstructive results are poor. To solve this problem, variable voltage CT reconstruction has been proposed. The effective projective sequences of a structural component are obtained through the variable voltage. The total variation is adjusted and minimized to optimize the reconstructive results on the basis of iterative image using algebraic reconstruction technique （ART）. In the process of reconstruction, the reconstructive image of low voltage is used as an initial value of the effective proiective reconstruction of the adjacent high voltage, and so on until to the highest voltage according to the gray weighted algorithm. Thereby the complete structural information is reconstructed. Simulation results show that the proposed algorithm can completely reflect the information of a complicated structural com- ponent, and the pixel values are more stable than those of the conventional.
基金supported by the National Natural Science Foundation of China(Grant Nos.11275204,11475175,and 11405175)the China Postdoctoral Science Foundation(Grant No.2017M612097)the Fundamental Research Funds for the Central Universities(Grant No.WK2310000065)
文摘X-ray grating interferometer has attracted widely attention in the past years due to its capability in achieving x-ray phase contrast imaging with low brilliance source. However, the widely used phase stepping information extraction method reduces system stability and prolongs data acquisition time by several times compared with conventional x-ray absorption- based imaging. The mechanical stepping can be avoided by using a staggered grating, but at the cost of low vertical spatial resolution. In this paper, employing a modified staggered grating and the angular signal radiography, we proposed a single-shot grating-based x-ray differential phase contrast imaging with decent vertical spatial resolution. The theoretical framework was deduced and proved by numerical experiments. Absorption, phase, and scattering computed tomography can be performed without phase stepping. Therefore, we believe this fast and highly stable imaging method with decent resolution would be widely applied in x-ray grating-based phase contrast imaging.
文摘We present a new iterative reconstruction algorithm to improve the algebraic reconstruction technique （ART） for the Single-Photon Emission Computed Tomography. Our method is a generalization of the Kaczmarz iterative algorithm for solving linear systems of equations and introduces exact and implicit attenuation correction derived from the attenuated Radon transform operator at each step of the algorithm. The performances of the presented algorithm have been tested upon various numerical experiments in presence of both strongly non-uniform attenuation and incomplete measurements data. We also tested the ability of our algorithm to handle moderate noisy data. Simulation studies demonstrate that the proposed method has a significant improvement in the quality of reconstructed images over ART. Moreover, convergence speed was improved and stability was established, facing noisy data, once we incorporate filtration procedure in our algorithm.
