With the development of the compressive sensing theory, the image reconstruction from the projections viewed in limited angles is one of the hot problems in the research of computed tomography technology. This paper d...With the development of the compressive sensing theory, the image reconstruction from the projections viewed in limited angles is one of the hot problems in the research of computed tomography technology. This paper develops an iterative algorithm for image reconstruction, which can fit the most cases. This method gives an image reconstruction flow with the difference image vector, which is based on the concept that the difference image vector between the reconstructed and the reference image is sparse enough. Then the l1-norm minimization method is used to reconstruct the difference vector to recover the image for flat subjects in limited angles. The algorithm has been tested with a thin planar phantom and a real object in limited-view projection data. Moreover, all the studies showed the satisfactory results in accuracy at a rather high reconstruction speed.展开更多
Objective To evaluate the feasibility of using a low concentration of contrast medium （Visipaque 270 mgl/mL）, low tube voltage, and an advanced image reconstruction algorithm in head and neck computed tomography ang...Objective To evaluate the feasibility of using a low concentration of contrast medium （Visipaque 270 mgl/mL）, low tube voltage, and an advanced image reconstruction algorithm in head and neck computed tomography angiography （CTA）. Methods Forty patients （22 men and 18 women; average age 48.7 ± 14.25 years; average body mass index 23.9 ± 3.7 kg/m^2） undergoing CTA for suspected vascular diseases were randomly assigned into two groups. Group A （n = 20） was administered 370 mgl/mL contrast medium, and group B （n = 20） was administered 270 mgl/mL contrast medium. Both groups were administered at a rate of 4.8 mL/s and an injection volume of 0.8 mL/kg. Images of group A were obtained with 120 kVp and filtered back projection （FBP） reconstruction, whereas images of group B were obtained with 80 kVp and 80% adaptive iterative statistical reconstruction algorithm （ASiR）. The CT values and standard deviations of intracranial arteries and image noise on the corona radiata were measured to calculate the contrast-to-noise ratio （CNR） and signal-to-noise ratio （SNR）. The beam-hardening artifacts （BHAs） around the skull base were calculated. Two readers evaluated the image quality with volume rendered images using scores from 1 to 5. The values between the two groups were statistically compared. Results The mean CT value of the intracranial arteries in group B was significantly higher than that in group A （P 〈 0.001）. The CNR and SNR values in group B were also statistically higher than those in group A （P 〈 0.001）. Image noise and BHAs were not significantly different between the two groups. The image quality score of VR images of in group B was significantly higher than that in group A （P = 0.001）. However, the quality scores of axial enhancement images in group B became significantly smaller than those in group A （P〈 0.001）. The CT dose index volume and dose-length product were decreased by 63.8% and 64%, respectively, in group B （P 〈 0.001 for both）. Conclusion Visipaque combined with 80 kVp and 80% ASiR provided similar image quality in intracranial CTA with 64% radiation dose reduction compared with the use of lopamidol, 120 kVp, and FBP reconstruc-tion.展开更多
The projection matrix model is used to describe the physical relationship between reconstructed object and projection.Such a model has a strong influence on projection and backprojection,two vital operations in iterat...The projection matrix model is used to describe the physical relationship between reconstructed object and projection.Such a model has a strong influence on projection and backprojection,two vital operations in iterative computed tomographic reconstruction.The distance-driven model（DDM） is a state-of-the-art technology that simulates forward and back projections.This model has a low computational complexity and a relatively high spatial resolution;however,it includes only a few methods in a parallel operation with a matched model scheme.This study introduces a fast and parallelizable algorithm to improve the traditional DDM for computing the parallel projection and backprojection operations.Our proposed model has been implemented on a GPU（graphic processing unit） platform and has achieved satisfactory computational efficiency with no approximation.The runtime for the projection and backprojection operations with our model is approximately 4.5 s and 10.5 s per loop,respectively,with an image size of 256×256×256 and 360 projections with a size of 512×512.We compare several general algorithms that have been proposed for maximizing GPU efficiency by using the unmatched projection/backprojection models in a parallel computation.The imaging resolution is not sacrificed and remains accurate during computed tomographic reconstruction.展开更多
In application of tomography imaging, limited-angle problem is a quite practical and important issue.In this paper, an iterative reprojection-reconstruction(IRR) algorithm using a modified Papoulis-Gerchberg(PG)iterat...In application of tomography imaging, limited-angle problem is a quite practical and important issue.In this paper, an iterative reprojection-reconstruction(IRR) algorithm using a modified Papoulis-Gerchberg(PG)iterative scheme is developed for reconstruction from limited-angle projections which contain noise. The proposed algorithm has two iterative update processes, one is the extrapolation of unknown data, and the other is the modification of the known noisy observation data. And the algorithm introduces scaling factors to control the two processes, respectively. The convergence of the algorithm is guaranteed, and the method of choosing the scaling factors is given with energy constraints. The simulation result demonstrates our conclusions and indicates that the algorithm proposed in this paper can obviously improve the reconstruction quality.展开更多
Although emission spectral tomography （EST） combines emission spectral measurement with optical computed tomography （OCT）, it is difficult to gain transient emission data from a large number of views, therefore, h...Although emission spectral tomography （EST） combines emission spectral measurement with optical computed tomography （OCT）, it is difficult to gain transient emission data from a large number of views, therefore, high precision OCT algorithms with few views ought to be studied for EST application. To improve the reconstruction precision in the case of few views, a new computed tomography reconstruction algorithm based on multipurpose optimal criterion and simulated annealing theory （multi-criterion simulated annealing reconstruction technique, MCSART） is proposed. This algorithm can suffice criterion of least squares, criterion of most uniformity, and criterion of most smoothness synchronously. We can get global optimal solution by MCSART algorithm with simulated annealing theory. The simulating experiment result shows that this algorithm is superior to the traditional algorithms under various noises.展开更多
The paper presents the implementation of a parallel version of FDK （Felkamp, David e Kress） algorithm using graphics processing units. Discussion was briefly some elements the computed tomographic scan and FDK algor...The paper presents the implementation of a parallel version of FDK （Felkamp, David e Kress） algorithm using graphics processing units. Discussion was briefly some elements the computed tomographic scan and FDK algorithm; and some ideas about GPUs （Graphics Processing Units） and its use in general purpose computing were presented. The paper shows a computational implementation of FDK algorithm and the process of parallelization of this implementation. Compare the parallel version of the algorithm with the sequential version, used speedup as a performance metric. To evaluate the performance of parallel version, two GPUs, GeForce 9400GT （16 cores） a low capacity GPU and Quadro 2000 （192 cores） a medium capacity GPU was reached speedup of 3.37.展开更多
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.展开更多
The order of the projection in the algebraic reconstruction technique(ART)method has great influence on the rate of the convergence.Although many scholars have studied the order of the projection,few theoretical proof...The order of the projection in the algebraic reconstruction technique(ART)method has great influence on the rate of the convergence.Although many scholars have studied the order of the projection,few theoretical proofs are given.Thomas Strohmer and Roman Vershynin introduced a randomized version of the Kaczmarz method for consistent,and over-determined linear systems and proved whose rate does not depend on the number of equations in the systems in 2009.In this paper,we apply this method to computed tomography(CT)image reconstruction and compared images generated by the sequential Kaczmarz method and the randomized Kaczmarz method.Experiments demonstrates the feasibility of the randomized Kaczmarz algorithm in CT image reconstruction and its exponential curve convergence.展开更多
目的比较不同水平的高级建模迭代重建算法(ADMIRE)对颞骨CT图像质量的影响。方法收集2023年1月至2月在空军军医大学唐都医院行颞骨CT扫描的患者24例,其中男性13例,女性11例;年龄21~61岁,平均年龄43.3岁。采用德国西门子SOMATOM FORCE C...目的比较不同水平的高级建模迭代重建算法(ADMIRE)对颞骨CT图像质量的影响。方法收集2023年1月至2月在空军军医大学唐都医院行颞骨CT扫描的患者24例,其中男性13例,女性11例;年龄21~61岁,平均年龄43.3岁。采用德国西门子SOMATOM FORCE CT扫描仪对24例患者行颞骨CT扫描,分别采用传统滤波反投影(FBP)和AD MIRE 1、3、5对扫描后的图像进行重建,比较4组图像客观评价(包括CT值、噪声、信噪比)、主观评价的差异。结果4组图像的客观评价、主观评价总分差异均有统计学意义(P<0.05);相对于FBP,ADMIRE重建后,图像的噪声减低,信噪比升高,ADMIRE 5噪声最小,信噪比最高;ADMIRE 1与FBP组图像主观评价差异无统计学意义,ADMIRE 3图像质量最好,5分占比最高(75%)。结论采用高级建模迭代重建算法后,图像噪声减低,信噪比提高,但并非迭代比例越高,图像质量越好。展开更多
Metal objects in X-ray computed tomography can cause severe artifacts.The state-of-the-art metal artifact reduction methods are in the sinogram inpainting category and are iterative methods.This paper proposes a proje...Metal objects in X-ray computed tomography can cause severe artifacts.The state-of-the-art metal artifact reduction methods are in the sinogram inpainting category and are iterative methods.This paper proposes a projectiondomain algorithm to reduce the metal artifacts.In this algorithm,the unknowns are the metal-affected projections,while the objective function is set up in the image domain.The data fidelity term is not utilized in the objective function.The objective function of the proposed algorithm consists of two terms:the total variation of the metalremoved image and the energy of the negative-valued pixels in the image.After the metal-affected projections are modified,the final image is reconstructed via the filtered backprojection algorithm.The feasibility of the proposed algorithm has been verified by real experimental data.展开更多
基金Project supported by the National Basic Research Program of China(Grant No.2006CB7057005)the National High Technology Research and Development Program of China(Grant No.2009AA012200)the National Natural Science Foundation of China (Grant No.60672104)
文摘With the development of the compressive sensing theory, the image reconstruction from the projections viewed in limited angles is one of the hot problems in the research of computed tomography technology. This paper develops an iterative algorithm for image reconstruction, which can fit the most cases. This method gives an image reconstruction flow with the difference image vector, which is based on the concept that the difference image vector between the reconstructed and the reference image is sparse enough. Then the l1-norm minimization method is used to reconstruct the difference vector to recover the image for flat subjects in limited angles. The algorithm has been tested with a thin planar phantom and a real object in limited-view projection data. Moreover, all the studies showed the satisfactory results in accuracy at a rather high reconstruction speed.
文摘Objective To evaluate the feasibility of using a low concentration of contrast medium （Visipaque 270 mgl/mL）, low tube voltage, and an advanced image reconstruction algorithm in head and neck computed tomography angiography （CTA）. Methods Forty patients （22 men and 18 women; average age 48.7 ± 14.25 years; average body mass index 23.9 ± 3.7 kg/m^2） undergoing CTA for suspected vascular diseases were randomly assigned into two groups. Group A （n = 20） was administered 370 mgl/mL contrast medium, and group B （n = 20） was administered 270 mgl/mL contrast medium. Both groups were administered at a rate of 4.8 mL/s and an injection volume of 0.8 mL/kg. Images of group A were obtained with 120 kVp and filtered back projection （FBP） reconstruction, whereas images of group B were obtained with 80 kVp and 80% adaptive iterative statistical reconstruction algorithm （ASiR）. The CT values and standard deviations of intracranial arteries and image noise on the corona radiata were measured to calculate the contrast-to-noise ratio （CNR） and signal-to-noise ratio （SNR）. The beam-hardening artifacts （BHAs） around the skull base were calculated. Two readers evaluated the image quality with volume rendered images using scores from 1 to 5. The values between the two groups were statistically compared. Results The mean CT value of the intracranial arteries in group B was significantly higher than that in group A （P 〈 0.001）. The CNR and SNR values in group B were also statistically higher than those in group A （P 〈 0.001）. Image noise and BHAs were not significantly different between the two groups. The image quality score of VR images of in group B was significantly higher than that in group A （P = 0.001）. However, the quality scores of axial enhancement images in group B became significantly smaller than those in group A （P〈 0.001）. The CT dose index volume and dose-length product were decreased by 63.8% and 64%, respectively, in group B （P 〈 0.001 for both）. Conclusion Visipaque combined with 80 kVp and 80% ASiR provided similar image quality in intracranial CTA with 64% radiation dose reduction compared with the use of lopamidol, 120 kVp, and FBP reconstruc-tion.
基金supported by the National High Technology Research and Development Program of China(Grant No.2012AA011603)the National Natural Science Foundation of China(Grant No.61372172)
文摘The projection matrix model is used to describe the physical relationship between reconstructed object and projection.Such a model has a strong influence on projection and backprojection,two vital operations in iterative computed tomographic reconstruction.The distance-driven model（DDM） is a state-of-the-art technology that simulates forward and back projections.This model has a low computational complexity and a relatively high spatial resolution;however,it includes only a few methods in a parallel operation with a matched model scheme.This study introduces a fast and parallelizable algorithm to improve the traditional DDM for computing the parallel projection and backprojection operations.Our proposed model has been implemented on a GPU（graphic processing unit） platform and has achieved satisfactory computational efficiency with no approximation.The runtime for the projection and backprojection operations with our model is approximately 4.5 s and 10.5 s per loop,respectively,with an image size of 256×256×256 and 360 projections with a size of 512×512.We compare several general algorithms that have been proposed for maximizing GPU efficiency by using the unmatched projection/backprojection models in a parallel computation.The imaging resolution is not sacrificed and remains accurate during computed tomographic reconstruction.
文摘In application of tomography imaging, limited-angle problem is a quite practical and important issue.In this paper, an iterative reprojection-reconstruction(IRR) algorithm using a modified Papoulis-Gerchberg(PG)iterative scheme is developed for reconstruction from limited-angle projections which contain noise. The proposed algorithm has two iterative update processes, one is the extrapolation of unknown data, and the other is the modification of the known noisy observation data. And the algorithm introduces scaling factors to control the two processes, respectively. The convergence of the algorithm is guaranteed, and the method of choosing the scaling factors is given with energy constraints. The simulation result demonstrates our conclusions and indicates that the algorithm proposed in this paper can obviously improve the reconstruction quality.
基金This work was supported by the Chinese Natural Science Foundation of China(No.60577016)the Foundation(No. 0512034)of Jiangxi Natural Science+1 种基金the Science and Technology Program(No. 2006-164)of Jiangxi Provincial Department of Educationthe Program(No.2005-314)of Key Laboratory of Nondestructive Testing Technology,Ministry of Education.
文摘Although emission spectral tomography （EST） combines emission spectral measurement with optical computed tomography （OCT）, it is difficult to gain transient emission data from a large number of views, therefore, high precision OCT algorithms with few views ought to be studied for EST application. To improve the reconstruction precision in the case of few views, a new computed tomography reconstruction algorithm based on multipurpose optimal criterion and simulated annealing theory （multi-criterion simulated annealing reconstruction technique, MCSART） is proposed. This algorithm can suffice criterion of least squares, criterion of most uniformity, and criterion of most smoothness synchronously. We can get global optimal solution by MCSART algorithm with simulated annealing theory. The simulating experiment result shows that this algorithm is superior to the traditional algorithms under various noises.
文摘The paper presents the implementation of a parallel version of FDK （Felkamp, David e Kress） algorithm using graphics processing units. Discussion was briefly some elements the computed tomographic scan and FDK algorithm; and some ideas about GPUs （Graphics Processing Units） and its use in general purpose computing were presented. The paper shows a computational implementation of FDK algorithm and the process of parallelization of this implementation. Compare the parallel version of the algorithm with the sequential version, used speedup as a performance metric. To evaluate the performance of parallel version, two GPUs, GeForce 9400GT （16 cores） a low capacity GPU and Quadro 2000 （192 cores） a medium capacity GPU was reached speedup of 3.37.
文摘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.
基金National Natural Science Foundation of China(No.61171179,No.61171178)Natural Science Foundation of Shanxi Province(No.2010011002-1,No.2010011002-2and No.2012021011-2)
文摘The order of the projection in the algebraic reconstruction technique(ART)method has great influence on the rate of the convergence.Although many scholars have studied the order of the projection,few theoretical proofs are given.Thomas Strohmer and Roman Vershynin introduced a randomized version of the Kaczmarz method for consistent,and over-determined linear systems and proved whose rate does not depend on the number of equations in the systems in 2009.In this paper,we apply this method to computed tomography(CT)image reconstruction and compared images generated by the sequential Kaczmarz method and the randomized Kaczmarz method.Experiments demonstrates the feasibility of the randomized Kaczmarz algorithm in CT image reconstruction and its exponential curve convergence.
文摘目的比较不同水平的高级建模迭代重建算法(ADMIRE)对颞骨CT图像质量的影响。方法收集2023年1月至2月在空军军医大学唐都医院行颞骨CT扫描的患者24例,其中男性13例,女性11例;年龄21~61岁,平均年龄43.3岁。采用德国西门子SOMATOM FORCE CT扫描仪对24例患者行颞骨CT扫描,分别采用传统滤波反投影(FBP)和AD MIRE 1、3、5对扫描后的图像进行重建,比较4组图像客观评价(包括CT值、噪声、信噪比)、主观评价的差异。结果4组图像的客观评价、主观评价总分差异均有统计学意义(P<0.05);相对于FBP,ADMIRE重建后,图像的噪声减低,信噪比升高,ADMIRE 5噪声最小,信噪比最高;ADMIRE 1与FBP组图像主观评价差异无统计学意义,ADMIRE 3图像质量最好,5分占比最高(75%)。结论采用高级建模迭代重建算法后,图像噪声减低,信噪比提高,但并非迭代比例越高,图像质量越好。
基金This research is partially supported by NIH,No.R15EB024283.
文摘Metal objects in X-ray computed tomography can cause severe artifacts.The state-of-the-art metal artifact reduction methods are in the sinogram inpainting category and are iterative methods.This paper proposes a projectiondomain algorithm to reduce the metal artifacts.In this algorithm,the unknowns are the metal-affected projections,while the objective function is set up in the image domain.The data fidelity term is not utilized in the objective function.The objective function of the proposed algorithm consists of two terms:the total variation of the metalremoved image and the energy of the negative-valued pixels in the image.After the metal-affected projections are modified,the final image is reconstructed via the filtered backprojection algorithm.The feasibility of the proposed algorithm has been verified by real experimental data.
