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.展开更多
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.展开更多
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.展开更多
Objective To compare the impact of different reconstruction algorithms on the image quality of 60 kVp head and neck CT angiography(CTA)using subjective and objective metrics,with a focus on vessel edge sharpness.Metho...Objective To compare the impact of different reconstruction algorithms on the image quality of 60 kVp head and neck CT angiography(CTA)using subjective and objective metrics,with a focus on vessel edge sharpness.Methods This prospective study enrolled 45 patients who underwent ultra-low-voltage(60 kVp)head and neck CTA.Image datasets were reconstructed with filtered back-projection(FBP),ClearView(CV)and ClearInfinity(CI)algorithms at low(30%),medium(50%),and high(70%)strengths.Image quality was assessed subjectively and objectively via the Kruskal‒Wallis test for multiple comparisons.Objective parameters,including edge rise slope(ERS)and edge rise distance(ERD),were analyzed via the Friedman test of multiple comparisons statistics.Results Subjective assessments favored the CI50 reconstruction algorithm,demonstrating superior or satisfactory results compared to the other algorithms,with significantly better vessel delineation,edge definition and diagnostic confidence(all P<0.05).Objective analysis revealed that the CV50 and CV70 algorithms significantly reduced ERS and/or elevated ERD(both P<0.05).However,the CI50 algorithm maintained comparable vessel edge sharpness(P>0.05)across all evaluated head and neck vascular segments when compared with the FBP algorithm.Conclusions The CI50 reconstruction algorithm optimizes image quality in 60 kVp head and neck CTA.It provides vessel edge sharpness comparable to FBP while offering superior vessel delineation,edge definition,and diagnostic confidence compared to FBP and CV algorithm.These findings suggest that CI50 has the potential to improve diagnostic accuracy in low-dose vascular imaging.展开更多
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.展开更多
In this paper, an approximate analytical algorithm in the form of direct Fourier reconstruction is obtained for the recon- struction of data functions arisen from ^-scheme short-scan sin- gle-photon emission computed ...In this paper, an approximate analytical algorithm in the form of direct Fourier reconstruction is obtained for the recon- struction of data functions arisen from ^-scheme short-scan sin- gle-photon emission computed tomography(SPECT) with uniform attenuation, and the modified central slice theorem is developed. Numerical simulations are conducted to demonstrate the effec- tiveness of the developed method.展开更多
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.展开更多
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.展开更多
文摘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.
基金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.
基金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 Grant from the National Key Research and Development Program of China(No.2024YFC2419300)the National Natural Science Foundation of China(No.82471967)+1 种基金the Hubei Provincial Key Research and Development Program(No.2024BCB008)the Hubei Provincial Natural Science Foundation of China(No.2025AFB733).
文摘Objective To compare the impact of different reconstruction algorithms on the image quality of 60 kVp head and neck CT angiography(CTA)using subjective and objective metrics,with a focus on vessel edge sharpness.Methods This prospective study enrolled 45 patients who underwent ultra-low-voltage(60 kVp)head and neck CTA.Image datasets were reconstructed with filtered back-projection(FBP),ClearView(CV)and ClearInfinity(CI)algorithms at low(30%),medium(50%),and high(70%)strengths.Image quality was assessed subjectively and objectively via the Kruskal‒Wallis test for multiple comparisons.Objective parameters,including edge rise slope(ERS)and edge rise distance(ERD),were analyzed via the Friedman test of multiple comparisons statistics.Results Subjective assessments favored the CI50 reconstruction algorithm,demonstrating superior or satisfactory results compared to the other algorithms,with significantly better vessel delineation,edge definition and diagnostic confidence(all P<0.05).Objective analysis revealed that the CV50 and CV70 algorithms significantly reduced ERS and/or elevated ERD(both P<0.05).However,the CI50 algorithm maintained comparable vessel edge sharpness(P>0.05)across all evaluated head and neck vascular segments when compared with the FBP algorithm.Conclusions The CI50 reconstruction algorithm optimizes image quality in 60 kVp head and neck CTA.It provides vessel edge sharpness comparable to FBP while offering superior vessel delineation,edge definition,and diagnostic confidence compared to FBP and CV algorithm.These findings suggest that CI50 has the potential to improve diagnostic accuracy in low-dose vascular imaging.
文摘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.
基金Supported by the National Natural Science Foundation of China(61271398)the Natural Science Foundation of Ningbo(2012A610031)
文摘In this paper, an approximate analytical algorithm in the form of direct Fourier reconstruction is obtained for the recon- struction of data functions arisen from ^-scheme short-scan sin- gle-photon emission computed tomography(SPECT) with uniform attenuation, and the modified central slice theorem is developed. Numerical simulations are conducted to demonstrate the effec- tiveness of the developed method.
基金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.
基金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.