Premium
An iterative method for tomographic x‐ray perfusion estimation in a decomposition model‐based approach
Author(s) -
Neukirchen Christoph,
Giordano Marco,
Wiesner Steffen
Publication year - 2010
Publication title -
medical physics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.473
H-Index - 180
eISSN - 2473-4209
pISSN - 0094-2405
DOI - 10.1118/1.3495818
Subject(s) - iterative reconstruction , imaging phantom , attenuation , computer science , projection (relational algebra) , tomographic reconstruction , tomography , correction for attenuation , radon transform , algorithm , regularization (linguistics) , noise (video) , computer vision , artificial intelligence , optics , physics , image (mathematics)
Purpose: X‐ray based tomographic blood perfusion imaging requires recovery of contrast time‐attenuation‐curves from dynamic projection data. When using slowly rotating imaging systems, this task is challenging due to nonsimultaneous projection acquisition. A dynamic reconstruction method is proposed that aims at compensating the lack of simultaneously acquired information by incorporating prior knowledge about the expected temporal contrast dynamics. Methods: A decomposition model using temporal basis functions to approximate time‐attenuation‐curves is integrated into an iterative tomographic reconstruction method. The computationally efficient implementation of the proposed approach makes use of standard forward‐projections and backprojections, as well as scalar products in image space. The critical issue of projection noise propagation is tackled by the application of regularization which is realized by the early stopping of iteration cycles and by the proper selection of smooth temporal basis functions. The performance of the proposed dynamic reconstruction approach is evaluated in a simulation study concerning various aspects: Noise propagation and regularization, specification of the temporal model, and type of acquisition mode. Results: The evaluation based on dynamic phantom data indicates that tomographic recovery of contrast time‐attenuation‐curves in tissue can be achieved with an average range of accuracy of ∼ 2 % (with respect to dynamic peak attenuation) under ideal noise‐free conditions. The relative estimation error for arterial time‐attenuation‐curves is in the range of 8%, which is due to faster contrast dynamics in the artery. In general, performance depends on the level of acquired information contained in the projection data, which is mainly influenced by the type of rotational acquisition mode; restrictions in angular range and speed can lead to limited accuracy. The analysis of propagated projection noise in a statistical bias‐variance framework reveals relative noise levels in estimated time‐attenuation‐curves of 3%–4% in tissue regions and below 1% in vessels when using optimized settings for regularization. Here, the effect of noise suppression depends on the interrelation between the number of iteration cycles and the constraints imposed by the temporal decomposition model. Conclusions: For usage with slowly rotating imaging systems, the presented model‐based iterative dynamic reconstruction method is capable of recovering contrast time‐attenuation‐curves related to tissue perfusion. The proposed regularization framework is an effective means to limit the impact of projection noise, which is a factor dominating estimation accuracy in tissue regions.