Robust Low-Dose CT Perfusion Deconvolution via Tensor Total-Variation Regularization

• Published in: IEEE transactions on medical imaging 2015-07, Vol.34 (7), p.1533-1548
• Main Authors: Ruogu Fang, Shaoting Zhang, Tsuhan Chen, Sanelli, Pina C.
• Format: Article
• Language: English
• Subjects: Algorithms; Brain; Computation; Computed tomography; Computed tomography perfusion; Deconvolution; Diseases; Estimation; Image reconstruction; low-dose; Mathematical analysis; Mathematical models; Mortality; Noise; Phantoms; Radiation dosage; radiation dose safety; regularization; Tensile stress; tensor total variation; Tensors
• ISSN: 0278-0062; 1558-254X
• DOI: 10.1109/TMI.2015.2405015

Acute brain diseases such as acute strokes and transit ischemic attacks are the leading causes of mortality and morbidity worldwide, responsible for 9% of total death every year. "Time is brain" is a widely accepted concept in acute cerebrovascular disease treatment. Efficient and accurate computational framework for hemodynamic parameters estimation can save critical time for thrombolytic therapy. Meanwhile the high level of accumulated radiation dosage due to continuous image acquisition in CT perfusion (CTP) raised concerns on patient safety and public health. However, low-radiation leads to increased noise and artifacts which require more sophisticated and time-consuming algorithms for robust estimation. In this paper, we focus on developing a robust and efficient framework to accurately estimate the perfusion parameters at low radiation dosage. Specifically, we present a tensor total-variation (TTV) technique which fuses the spatial correlation of the vascular structure and the temporal continuation of the blood signal flow. An efficient algorithm is proposed to find the solution with fast convergence and reduced computational complexity. Extensive evaluations are carried out in terms of sensitivity to noise levels, estimation accuracy, contrast preservation, and performed on digital perfusion phantom estimation, as well as in vivo clinical subjects. Our framework reduces the necessary radiation dose to only 8% of the original level and outperforms the state-of-art algorithms with peak signal-to-noise ratio improved by 32%. It reduces the oscillation in the residue functions, corrects over-estimation of cerebral blood flow (CBF) and under-estimation of mean transit time (MTT), and maintains the distinction between the deficit and normal regions.