Dataset Information

Low-Dose Dynamic Cerebral Perfusion Computed Tomography Reconstruction via Kronecker-Basis-Representation Tensor Sparsity Regularization.

ABSTRACT: Dynamic cerebral perfusion computed tomography (DCPCT) has the ability to evaluate the hemodynamic information throughout the brain. However, due to multiple 3-D image volume acquisitions protocol, DCPCT scanning imposes high radiation dose on the patients with growing concerns. To address this issue, in this paper, based on the robust principal component analysis (RPCA, or equivalently the low-rank and sparsity decomposition) model and the DCPCT imaging procedure, we propose a new DCPCT image reconstruction algorithm to improve low-dose DCPCT and perfusion maps quality via using a powerful measure, called Kronecker-basis-representation tensor sparsity regularization, for measuring low-rankness extent of a tensor. For simplicity, the first proposed model is termed tensor-based RPCA (T-RPCA). Specifically, the T-RPCA model views the DCPCT sequential images as a mixture of low-rank, sparse, and noise components to describe the maximum temporal coherence of spatial structure among phases in a tensor framework intrinsically. Moreover, the low-rank component corresponds to the "background" part with spatial-temporal correlations, e.g., static anatomical contribution, which is stationary over time about structure, and the sparse component represents the time-varying component with spatial-temporal continuity, e.g., dynamic perfusion enhanced information, which is approximately sparse over time. Furthermore, an improved nonlocal patch-based T-RPCA (NL-T-RPCA) model which describes the 3-D block groups of the "background" in a tensor is also proposed. The NL-T-RPCA model utilizes the intrinsic characteristics underlying the DCPCT images, i.e., nonlocal self-similarity and global correlation. Two efficient algorithms using alternating direction method of multipliers are developed to solve the proposed T-RPCA and NL-T-RPCA models, respectively. Extensive experiments with a digital brain perfusion phantom, preclinical monkey data, and clinical patient data clearly demonstrate that the two proposed models can achieve more gains than the existing popular algorithms in terms of both quantitative and visual quality evaluations from low-dose acquisitions, especially as low as 20 mAs.

SUBMITTER: Zeng D

PROVIDER: S-EPMC5711606 | biostudies-literature | 2017 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Low-Dose Dynamic Cerebral Perfusion Computed Tomography Reconstruction via Kronecker-Basis-Representation Tensor Sparsity Regularization.

Zeng Dong D Xie Qi Q Cao Wenfei W Lin Jiahui J Zhang Hao H Zhang Shanli S Huang Jing J Bian Zhaoying Z Meng Deyu D Xu Zongben Z Liang Zhengrong Z Chen Wufan W Ma Jianhua J

IEEE transactions on medical imaging 20170904 12

Dynamic cerebral perfusion computed tomography (DCPCT) has the ability to evaluate the hemodynamic information throughout the brain. However, due to multiple 3-D image volume acquisitions protocol, DCPCT scanning imposes high radiation dose on the patients with growing concerns. To address this issue, in this paper, based on the robust principal component analysis (RPCA, or equivalently the low-rank and sparsity decomposition) model and the DCPCT imaging procedure, we propose a new DCPCT image r ...[more]

PMID: 28880164

Similar Datasets

Project description:BackgroundMulti-energy computed tomography (CT) provides multiple channel-wise reconstructed images, and they can be used for material identification and k-edge imaging. Nonetheless, the projection datasets are frequently corrupted by various noises (e.g., electronic, Poisson) in the acquisition process, resulting in lower signal-noise-ratio (SNR) measurements. Multi-energy CT images have local sparsity, nonlocal self-similarity in spatial dimension, and correlation in spectral dimension.MethodsIn this paper, we propose an image-spectral decomposition extended-learning assisted by sparsity (IDEAS) method to fully exploit these intrinsic priors for multi-energy CT image reconstruction. Particularly, a nonlocal low-rank Tucker decomposition (TD) is employed to utilize the correlation and nonlocal self-similarity priors. Moreover, considering the advantages of multi-task tensor dictionary learning (TDL) in sparse representation, an adaptive spatial dictionary and an adaptive spectral dictionary are trained during the iterative reconstruction process. Furthermore, a weighted total variation (TV) regularization term is employed to encourage local sparsity.ResultsNumerical simulation, physical phantom, and preclinical mouse experiments are performed to validate the proposed IDEAS algorithm. Specifically, in the simulation experiments, the proposed IDEAS reconstructed high-quality images that are very close to the references. For example, the root mean square error (RMSE) of IDEAS image in energy bin 1 is as low as 0.0672, while the RMSE of other methods are higher than 0.0843. Besides, the structural similarity (SSIM) of IDEAS reconstructed image in energy bin 1 is greater than 0.98. For material decomposition, the RMSE of IDEAS bone component is as low as 0.0152, and other methods are higher than 0.0199. In addition, the computational cost of IDEAS is as low as 98.8 s for one iteration, and the competing tensor decomposition method is higher than 327 s.ConclusionsTo further improve the quality of the reconstructed multi-energy CT images, multiple prior regularizations are introduced to the multi-energy CT reconstructed model, leading to an IDEAS method. Both qualitative and quantitative evaluation of our results confirm the outstanding performance of the proposed algorithm compared to the state-of-the-arts.

Project description:The recent advent of tensor tomography techniques has enabled tomographic investigations of the 3D nanostructure organization of biological and material science samples. These techniques extended the concept of conventional X-ray tomography by reconstructing not only a scalar value such as the attenuation coefficient per voxel, but also a set of parameters that capture the local anisotropy of nanostructures within every voxel of the sample. Tensor tomography data sets are intrinsically large as each pixel of a conventional X-ray projection is substituted by a scattering pattern, and projections have to be recorded at different sample angular orientations with several tilts of the rotation axis with respect to the X-ray propagation direction. Currently available reconstruction approaches for such large data sets are computationally expensive. Here, a novel, fast reconstruction algorithm, named iterative reconstruction tensor tomography (IRTT), is presented to simplify and accelerate tensor tomography reconstructions. IRTT is based on a second-rank tensor model to describe the anisotropy of the nanostructure in every voxel and on an iterative error backpropagation reconstruction algorithm to achieve high convergence speed. The feasibility and accuracy of IRTT are demonstrated by reconstructing the nanostructure anisotropy of three samples: a carbon fiber knot, a human bone trabecula specimen and a fixed mouse brain. Results and reconstruction speed were compared with those obtained by the small-angle scattering tensor tomography (SASTT) reconstruction method introduced by Liebi et al. [Nature (2015), 527, 349-352]. The principal orientation of the nanostructure within each voxel revealed a high level of agreement between the two methods. Yet, for identical data sets and computer hardware used, IRTT was shown to be more than an order of magnitude faster. IRTT was found to yield robust results, it does not require prior knowledge of the sample for initializing parameters, and can be used in cases where simple anisotropy metrics are sufficient, i.e. the tensor approximation adequately captures the level of anisotropy and the dominant orientation within a voxel. In addition, by greatly accelerating the reconstruction, IRTT is particularly suitable for handling large tomographic data sets of samples with internal structure or as a real-time analysis tool during the experiment for online feedback during data acquisition. Alternatively, the IRTT results might be used as an initial guess for models capturing a higher complexity of structural anisotropy such as spherical harmonics based SASTT in Liebi et al. (2015), improving both overall convergence speed and robustness of the reconstruction.

Project description:PurposeComputed tomography (CT)-derived ventilation methods compute respiratory induced volume changes as a surrogate for pulmonary ventilation. Currently, there are no known methods to derive perfusion information from noncontrast CT. We introduce a novel CT-Perfusion (CT-P) method for computing the magnitude mass changes apparent on dynamic noncontrast CT as a surrogate for pulmonary perfusion.MethodsCT-Perfusion is based on a mass conservation model which describes the unknown mass change as a linear combination of spatially corresponding inhale and exhale HU estimated voxel densities. CT-P requires a deformable image registration (DIR) between the inhale/exhale lung CT pair, a preprocessing lung volume segmentation, and an estimate for the Jacobian of the DIR transformation. Given this information, the CT-P image, which provides the magnitude mass change for each voxel within the lung volume, is formulated as the solution to a constrained linear least squares problem defined by a series of subregional mean magnitude mass change measurements. Similar to previous robust CT-ventilation methods, the amount of uncertainty in a subregional sample mean measurement is related to measurement resolution and can be characterized with respect to a tolerance parameter τ . Spatial Spearman correlation between single photon emission CT perfusion (SPECT-P) and the proposed CT-P method was assessed in two patient cohorts via a parameter sweep of τ . The first cohort was comprised of 15 patients diagnosed with pulmonary embolism (PE) who had SPECT-P and 4DCT imaging acquired within 24 h of PE diagnosis. The second cohort was comprised of 15 nonsmall cell lung cancer patients who had SPECT-P and 4DCT images acquired prior to radiotherapy. For each test case, CT-P images were computed for 30 different uncertainty parameter values, uniformly sampled from the range [0.01, 0.125], and the Spearman correlation between the SPECT-P and the resulting CT-P images were computed.ResultsThe median correlations between CT-P and SPECT-P taken over all 30 test cases ranged between 0.49 and 0.57 across the parameter sweep. For the optimal tolerance τ = 0.0385, the CT-P and SPECT-P correlations across all 30 test cases ranged between 0.02 and 0.82. A one-sample sign test was applied separately to the PE and lung cancer cohorts. A low Spearmen correlation of 15% was set as the null median value and two-sided alternative was tested. The PE patients showed a median correlation of 0.57 (IQR = 0.305). One-sample sign test was statistically significant with 96.5 % confidence interval: 0.20-0.63, P < 0.00001. Lung cancer patients had a median correlation of 0.57(IQR = 0.230). Again, a one-sample sign test for median was statistically significant with 96.5 percent confidence interval: 0.45-0.71, P < 0.00001.ConclusionCT-Perfusion is the first mechanistic model designed to quantify magnitude blood mass changes on noncontrast dynamic CT as a surrogate for pulmonary perfusion. While the reported correlations with SPECT-P are promising, further investigation is required to determine the optimal CT acquisition protocol and numerical method implementation for CT-P imaging.

Project description:See Saver (doi:10.1093/awx020) for a scientific commentary on this article.Stroke shortens an individual's disability-free life. We aimed to assess the relative prognostic influence of pre- and post-treatment perfusion computed tomography imaging variables (e.g. ischaemic core and penumbral volumes) compared to standard clinical predictors (such as onset-to-treatment time) on long-term stroke disability in patients undergoing thrombolysis. We used data from a prospectively collected international, multicentre, observational registry of acute ischaemic stroke patients who had perfusion computed tomography and computed tomography angiography before treatment with intravenous alteplase. Baseline perfusion computed tomography and follow-up magnetic resonance imaging were analysed to derive the baseline penumbra volume, baseline ischaemic core volume, and penumbra salvaged from infarction. The primary outcome measure was the effect of imaging and clinical variables on Disability-Adjusted Life Year. Clinical variables were age, sex, National Institutes of Health Stroke Scale score, and onset-to-treatment time. Age, sex, country, and 3-month modified Rankin Scale were extracted from the registry to calculate disability-adjusted life-year due to stroke, such that 1 year of disability-adjusted life-year equates to 1 year of healthy life lost due to stroke. There were 772 patients receiving alteplase therapy. The number of disability-adjusted life-year days lost per 1 ml of baseline ischaemic core volume was 17.5 (95% confidence interval, 13.2-21.9 days, P < 0.001). For every millilitre of penumbra salvaged, 7.2 days of disability-adjusted life-year days were saved (β = -7.2, 95% confidence interval, -10.4 to -4.1 days, P < 0.001). Each minute of earlier onset-to-treatment time resulted in a saving of 4.4 disability-free days after stroke (1.3-7.5 days, P = 0.006). However, after adjustment for imaging variables, onset-to-treatment time was not significantly associated with savings in disability-adjusted life-year days. Pretreatment perfusion computed tomography can (independently of clinical variables) predict significant gains, or loss, of disability-free life in patients undergoing reperfusion therapy for stroke. The effect of earlier treatment on disability-free life appears explained by salvage of penumbra, particularly when the ischaemic core is not too large.

Dataset Information

Low-Dose Dynamic Cerebral Perfusion Computed Tomography Reconstruction via Kronecker-Basis-Representation Tensor Sparsity Regularization.

Publications

Low-Dose Dynamic Cerebral Perfusion Computed Tomography Reconstruction via Kronecker-Basis-Representation Tensor Sparsity Regularization.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets