Perfusion measurements by micro-CT using prior image constrained compressed sensing (PICCS): initial phantom results

Iterative reconstruction for low-dose cerebral perfusion computed tomography using prior image induced diffusion tensor

Cerebral perfusion computed tomography (CPCT) can depict the functional status of cerebral circulation at the tissue level; hence, it has been increasingly used to diagnose patients with cerebrovascular disease. However, there is a significant concern that CPCT scanning protocol could expose patients to excessive radiation doses. Although reducing the x-ray tube current when acquiring CPCT projection data is an effective method for reducing radiation dose, this technique usually results in degraded image quality. To enhance the image quality of low-dose CPCT, we present a prior image induced diffusion tensor (PIDT) for statistical iterative reconstruction, based on the penalized weighted least-squares (PWLS) criterion, which we referred to as PWLS-PIDT, for simplicity. Specifically, PIDT utilizes the geometric features of pre-contrast scanned high-quality CT image as a structure prior for PWLS reconstruction; therefore, the low-dose CPCT images are enhanced while preserving important features in the target image. An effective alternating minimization algorithm is developed to solve the associated objective function in the PWLS-PIDT reconstruction. We conduct qualitative and quantitative studies to evaluate the PWLS-PIDT reconstruction with a digital brain perfusion phantom and patient data. With this method, the noise in the reconstructed CPCT images is more substantially reduced than that of other competing methods, without sacrificing structural details significantly. Furthermore, the CPCT sequential images reconstructed via the PWLS-PIDT method can derive more accurate hemodynamic parameter maps than those of other competing methods.

Iterative quality enhancement via residual-artifact learning networks for low-dose CT

Radiation exposure and the associated risk of cancer for patients in computed tomography (CT) scans have been major clinical concerns. The radiation exposure can be reduced effectively via lowering the x-ray tube current (mA). However, this strategy may lead to excessive noise and streak artifacts in the conventional filtered back-projection reconstructed images. To address this issue, some deep convolutional neural network (ConvNet) based approaches have been developed for low-dose CT imaging inspired by the recent development of machine learning. Nevertheless, some of the image textures reconstructed by the ConvNet could be corrupted by the severe streaks, especially in ultra-low-dose cases, which could be close to prostheses and hamper diagnosis. Therefore, in this work, we propose an iterative residual-artifact learning ConvNet (IRLNet) approach to improve the reconstruction performance over the ConvNet based approaches. Specifically, the proposed IRLNet estimates the high-frequency details within the noise and then removes them iteratively; after eliminating severe streaks in the low-dose CT images, the residual low-frequency details can be processed through the conventional network. Moreover, the proposed IRLNet scheme can be extended for robust handling of quantitative dual energy CT/cerebral perfusion CT imaging, and statistical iterative reconstruction. Real patient data are used to evaluate the proposed IRLNet, and the experimental results demonstrate that the proposed IRLNet approach outperforms the previous ConvNet based approaches in reducing the image noise and streak artifacts efficiently at the same time as preserving edge details well, suggesting that the proposed IRLNet approach can be used to improve the CT image quality, especially in ultra-low-dose cases.

Low-Dose CT Perfusion of the Liver using Reconstruction of Difference

Liver CT perfusion (CTP) is used in the detection, staging, and treatment response analysis of hepatic diseases. Unfortunately, CTP radiation exposures is significant, limiting more widespread use. Traditional CTP data processing reconstructs individual temporal samples, ignoring a large amount of shared anatomical information between temporal samples, suggesting opportunities for improved data processing. We adopt a prior-image-based reconstruction approach called Reconstruction of Difference (RoD) to enable low-exposure CTP acquisition. RoD differs from many algorithms by directly estimating the attenuation changes between the current patient state and a prior CT volume. We propose to use a high-fidelity unenhanced baseline CT image to integrate prior anatomical knowledge into subsequent data reconstructions. Using simulation studies based on a 4D digital anthropomorphic phantom with realistic time-attenuation curves, we compare RoD with conventional filtered-backprojection, penalized-likelihood estimation, and prior image penalized-likelihood estimation. We evaluate each method in comparisons of reconstructions at individual time points, accuracy of estimated time-attenuation curves, and in an analysis of common perfusion metric maps including hepatic arterial perfusion, hepatic portal perfusion, perfusion index, and time-to-peak. Results suggest that RoD enables significant exposure reductions, outperforming standard and more sophisticated model-based reconstruction, making RoD a potentially important tool to enable low-dose liver CTP.

Radiation Dose Reduction in CT Myocardial Perfusion Imaging Using SMART-RECON

In this paper, a newly developed statistical model-based image reconstruction [referred to as Simultaneous Multiple Artifacts Reduction in Tomographic RECONstruction (SMART-RECON)] is applied to low dose computer tomography (CT) myocardial perfusion imaging (CT-MPI). This method uses the nuclear norm of the spatial-temporal image matrix of the CT-MPI images as a regularizer, rather than a conventional spatial regularizer that incorporates image smoothness, edge preservation, or spatial sparsity into the reconstruction. In addition to providing the needed noise reduction for low-dose CT-MPI, SMART-RECON provides images with spatial resolution and noise power spectrum (NPS) properties, which are independent of contrast and dose levels. Both numerical simulations and in vivo animal studies were performed to validate the proposed method. In these studies, it was found that: 1) quantitative accuracy of perfusion maps in CT-MPI was well maintained for radiation dose level as low as 10 mAs per image frame, compared with the reference standard of 200 mAs for conventional filtered backprojection; 2) flow-occluded myocardium in the porcine heart was well delineated by SMART-RECON at 10 mAs per frame when compared with model-based image reconstruction using spatial total variation (TV) as the regularizer (referred to as TV-SIR) or spatial-temporal TV (ST-TV-SIR); the CT-MPI results were confirmed with positron-emission tomography imaging; 3) image sharpness in SMART-RECON images was nearly independent of image contrast level and radiation dose level, in stark contrast to TV-SIR and ST-TV-SIR, which displayed a strong dependence on both image contrast and radiation dose levels; and 4) the structure of the dose-normalized NPS for the SMART-RECON method did not depend on dose, while the TV-SIR and ST-TV-SIR NPS structure was dose-dependent.

Local attenuation curve optimization framework for high quality perfusion maps in low-dose cerebral perfusion CT

PURPOSE

Cerebral perfusion x-ray computed tomography (PCT) is a powerful tool for noninvasive imaging of hemodynamic information throughout the brain. Conventional PCT requires the brain to be imaged multiple times during the perfusion process, and hence radiation dose is a major concern. The authors propose a PCT reconstruction algorithm that allows for lowering the dose while maintaining a high quality of the perfusion maps. It relies on an accurate estimation of the arterial input function (AIF), which in turn depends on the quality of the attenuation curves in the arterial region.

METHODS

The authors propose the local attenuation curve optimization (LACO) framework. It accurately models the attenuation curves inside the vessel and arterial regions and optimizes its shape directly based on the acquired x-ray projection data.

RESULTS

The LACO algorithm is extensively validated with simulation and real clinical experiments. Quantitative and qualitative results show that our proposed approach accurately estimates the vessel and arterial attenuation curves from only few x-ray projections. In contrast to conventional approaches, where the AIF is estimated based on the reconstructed images, our method computes an optimal AIF directly based on the projection data, resulting in far more accurate perfusion maps.

CONCLUSIONS

The LACO algorithm allows estimating high quality perfusion maps in low dose scanning protocols.

Iterative reconstruction for CT perfusion with a prior-image induced hybrid nonlocal means regularization: Phantom studies

PURPOSE

In computed tomography perfusion (CTP) imaging, an initial phase CT acquired with a high-dose protocol can be used to improve the image quality of later phase CT acquired with a low-dose protocol. For dynamic regions, signals in the later low-dose CT may not be completely recovered if the initial CT heavily regularizes the iterative reconstruction process. The authors propose a hybrid nonlocal means (hNLM) regularization model for iterative reconstruction of low-dose CTP to overcome the limitation of the conventional prior-image induced penalty.

METHODS

The hybrid penalty was constructed by combining the NLM of the initial phase high-dose CT in the stationary region and later phase low-dose CT in the dynamic region. The stationary and dynamic regions were determined by the similarity between the initial high-dose scan and later low-dose scan. The similarity was defined as a Gaussian kernel-based distance between the patch-window of the same pixel in the two scans, and its measurement was then used to weigh the influence of the initial high-dose CT. For regions with high similarity (e.g., stationary region), initial high-dose CT played a dominant role for regularizing the solution. For regions with low similarity (e.g., dynamic region), the regularization relied on a low-dose scan itself. This new hNLM penalty was incorporated into the penalized weighted least-squares (PWLS) for CTP reconstruction. Digital and physical phantom studies were performed to evaluate the PWLS-hNLM algorithm.

RESULTS

Both phantom studies showed that the PWLS-hNLM algorithm is superior to the conventional prior-image induced penalty term without considering the signal changes within the dynamic region. In the dynamic region of the Catphan phantom, the reconstruction error measured by root mean square error was reduced by 42.9% in PWLS-hNLM reconstructed image.

CONCLUSIONS

The PWLS-hNLM algorithm can effectively use the initial high-dose CT to reconstruct low-dose CTP in the stationary region while reducing its influence in the dynamic region.

Low-dose cerebral perfusion computed tomography image restoration via low-rank and total variation regularizations

Cerebral perfusion x-ray computed tomography (PCT) is an important functional imaging modality for evaluating cerebrovascular diseases and has been widely used in clinics over the past decades. However, due to the protocol of PCT imaging with repeated dynamic sequential scans, the associative radiation dose unavoidably increases as compared with that used in conventional CT examinations. Minimizing the radiation exposure in PCT examination is a major task in the CT field. In this paper, considering the rich similarity redundancy information among enhanced sequential PCT images, we propose a low-dose PCT image restoration model by incorporating the low-rank and sparse matrix characteristic of sequential PCT images. Specifically, the sequential PCT images were first stacked into a matrix (i.e., low-rank matrix), and then a non-convex spectral norm/regularization and a spatio-temporal total variation norm/regularization were then built on the low-rank matrix to describe the low rank and sparsity of the sequential PCT images, respectively. Subsequently, an improved split Bregman method was adopted to minimize the associative objective function with a reasonable convergence rate. Both qualitative and quantitative studies were conducted using a digital phantom and clinical cerebral PCT datasets to evaluate the present method. Experimental results show that the presented method can achieve images with several noticeable advantages over the existing methods in terms of noise reduction and universal quality index. More importantly, the present method can produce more accurate kinetic enhanced details and diagnostic hemodynamic parameter maps.

Exploitation of temporal redundancy in compressed sensing reconstruction of fMRI studies with a prior-based algorithm (PICCS)

PURPOSE

Compressed sensing is a technique used to accelerate magnetic resonance imaging (MRI) acquisition without compromising image quality. While it has proven particularly useful in dynamic imaging procedures such as cardiac cine, very few authors have applied it to functional magnetic resonance imaging (fMRI). The purpose of the present study was to check whether the prior image constrained compressed sensing (PICCS) algorithm, which is based on an available prior image, can improve the statistical maps in fMRI better than other strategies that also exploit temporal redundancy.

METHODS

PICCS was compared to spatiotemporal total variation (TTV) and k-t FASTER, since they have already demonstrated high performance and robustness in other MRI applications, such as cardiac cine MRI and resting state fMRI, respectively. The prior image for PICCS was the average of all undersampled data. Both PICCS and TTV were solved using the split Bregman formulation. K-t FASTER algorithm relies on matrix completion to reconstruct the undersampled k-spaces. The three algorithms were evaluated using two datasets with high and low signal-to-noise ratio (SNR)-BOLD contrast-acquired in a 7 T preclinical MRI scanner and retrospectively undersampled at various rates (i.e., acceleration factors). The authors evaluated their performance in terms of the sensitivity/specificity of BOLD detection through receiver operating characteristic curves and by visual inspection of the statistical maps.

RESULTS

With high SNR studies, PICCS performed similarly to the state-of-the-art algorithms TTV and k-t FASTER and provided consistent BOLD signal at the ROI. In scenarios with low SNR and high acceleration factors, PICCS still provided consistent maps and higher sensitivity/specificity than TTV, whereas k-t FASTER failed to provide significant maps.

CONCLUSIONS

The authors performed a comparison between three reconstructions (PICCS, TTV, and k-t FASTER) that exploit temporal redundancy in fMRI. The prior-based algorithm, PICCS, preserved BOLD activation and sensitivity/specificity better than TTV and k-t FASTER in noisy scenarios. The PICCS algorithm can potentially reach an acceleration factor of ×8 and still provide BOLD contrast in the ROI with an area under the curve over 0.99.

A Novel Prior- and Motion-Based Compressed Sensing Method for Small-Animal Respiratory Gated CT

Low-dose protocols for respiratory gating in cardiothoracic small-animal imaging lead to streak artifacts in the images reconstructed with a Feldkamp-Davis-Kress (FDK) method. We propose a novel prior- and motion-based reconstruction (PRIMOR) method, which improves prior-based reconstruction (PBR) by adding a penalty function that includes a model of motion. The prior image is generated as the average of all the respiratory gates, reconstructed with FDK. Motion between respiratory gates is estimated using a nonrigid registration method based on hierarchical B-splines. We compare PRIMOR with an equivalent PBR method without motion estimation using as reference the reconstruction of high dose data. From these data acquired with a micro-CT scanner, different scenarios were simulated by changing photon flux and number of projections. Methods were evaluated in terms of contrast-to-noise-ratio (CNR), mean square error (MSE), streak artefact indicator (SAI), solution error norm (SEN), and correction of respiratory motion. Also, to evaluate the effect of each method on lung studies quantification, we have computed the Jaccard similarity index of the mask obtained from segmenting each image as compared to those obtained from the high dose reconstruction. Both iterative methods greatly improved FDK reconstruction in all cases. PBR was prone to streak artifacts and presented blurring effects in bone and lung tissues when using both a low number of projections and low dose. Adopting PBR as a reference, PRIMOR increased CNR up to 33% and decreased MSE, SAI and SEN up to 20%, 4% and 13%, respectively. PRIMOR also presented better compensation for respiratory motion and higher Jaccard similarity index. In conclusion, the new method proposed for low-dose respiratory gating in small-animal scanners shows an improvement in image quality and allows a reduction of dose or a reduction of the number of projections between two and three times with respect to previous PBR approaches.

Iterative image reconstruction for limited-angle inverse helical cone-beam computed tomography

Helical trajectory satisfying the condition of exact reconstruction, has been widely utilized in the commercial computed tomography (CT). While limited by the scanning environment in some practical applications, the conventional helical cone-beam CT imaging is hard to complete, thus, developing an imaging system suited for long-object may be valuable. Three-dimensional C-arm CT is an innovative imaging technique which has been greatly concerned. Since there is a high degree of freedom of C-arm, more flexible image acquisition trajectories for 3D imaging can be achieved. In this work, a fast iterative reconstruction algorithm based on total variation minimization is developed for a trajectory of limited-angle inverse helical cone-beam CT, which can be applied to detect long-object without slip-ring technology. The experimental results show that the developed algorithm can yield reconstructed images of low noise level and high image quality.

Investigation of different sparsity transforms for the PICCS algorithm in small-animal respiratory gated CT

Respiratory gating helps to overcome the problem of breathing motion in cardiothoracic small-animal imaging by acquiring multiple images for each projection angle and then assigning projections to different phases. When this approach is used with a dose similar to that of a static acquisition, a low number of noisy projections are available for the reconstruction of each respiratory phase, thus leading to streak artifacts in the reconstructed images. This problem can be alleviated using a prior image constrained compressed sensing (PICCS) algorithm, which enables accurate reconstruction of highly undersampled data when a prior image is available. We compared variants of the PICCS algorithm with different transforms in the prior penalty function: gradient, unitary, and wavelet transform. In all cases the problem was solved using the Split Bregman approach, which is efficient for convex constrained optimization. The algorithms were evaluated using simulations generated from data previously acquired on a micro-CT scanner following a high-dose protocol (four times the dose of a standard static protocol). The resulting data were used to simulate scenarios with different dose levels and numbers of projections. All compressed sensing methods performed very similarly in terms of noise, spatiotemporal resolution, and streak reduction, and filtered back-projection was greatly improved. Nevertheless, the wavelet domain was found to be less prone to patchy cartoon-like artifacts than the commonly used gradient domain.

Accelerated barrier optimization compressed sensing (ABOCS) for CT reconstruction with improved convergence

Recently, we proposed a new algorithm of accelerated barrier optimization compressed sensing (ABOCS) for iterative CT reconstruction. The previous implementation of ABOCS uses gradient projection (GP) with a Barzilai-Borwein (BB) step-size selection scheme (GP-BB) to search for the optimal solution. The algorithm does not converge stably due to its non-monotonic behavior. In this paper, we further improve the convergence of ABOCS using the unknown-parameter Nesterov (UPN) method and investigate the ABOCS reconstruction performance on clinical patient data. Comparison studies are carried out on reconstructions of computer simulation, a physical phantom and a head-and-neck patient. In all of these studies, the ABOCS results using UPN show more stable and faster convergence than those of the GP-BB method and a state-of-the-art Bregman-type method. As shown in the simulation study of the Shepp-Logan phantom, UPN achieves the same image quality as those of GP-BB and the Bregman-type methods, but reduces the iteration numbers by up to 50% and 90%, respectively. In the Catphan©600 phantom study, a high-quality image with relative reconstruction error (RRE) less than 3% compared to the full-view result is obtained using UPN with 17% projections (60 views). In the conventional filtered-backprojection reconstruction, the corresponding RRE is more than 15% on the same projection data. The superior performance of ABOCS with the UPN implementation is further demonstrated on the head-and-neck patient. Using 25% projections (91 views), the proposed method reduces the RRE from 21% as in the filtered backprojection (FBP) results to 7.3%. In conclusion, we propose UPN for ABOCS implementation. As compared to GP-BB and the Bregman-type methods, the new method significantly improves the convergence with higher stability and fewer iterations.

Iterative image reconstruction for sparse-view CT using normal-dose image induced total variation prior

X-ray computed tomography (CT) iterative image reconstruction from sparse-view projection data has been an important research topic for radiation reduction in clinic. In this paper, to relieve the requirement of misalignment reduction operation of the prior image constrained compressed sensing (PICCS) approach introduced by Chen et al, we present an iterative image reconstruction approach for sparse-view CT using a normal-dose image induced total variation (ndiTV) prior. The associative objective function of the present approach is constructed under the penalized weighed least-square (PWLS) criteria, which contains two terms, i.e., the weighted least-square (WLS) fidelity and the ndiTV prior, and is referred to as "PWLS-ndiTV". Specifically, the WLS fidelity term is built based on an accurate relationship between the variance and mean of projection data in the presence of electronic background noise. The ndiTV prior term is designed to reduce the influence of the misalignment between the desired- and prior- image by using a normal-dose image induced non-local means (ndiNLM) filter. Subsequently, a modified steepest descent algorithm is adopted to minimize the associative objective function. Experimental results on two different digital phantoms and an anthropomorphic torso phantom show that the present PWLS-ndiTV approach for sparse-view CT image reconstruction can achieve noticeable gains over the existing similar approaches in terms of noise reduction, resolution-noise tradeoff, and low-contrast object detection.

A comparative analysis of water-soluble and blood-pool contrast agents for in vivo vascular imaging with micro-CT

RATIONALE AND OBJECTIVES

In recent years, micro-computed tomography (micro-CT) has emerged as a high-resolution modality for vascular exploration in vivo. Several x-ray contrast agents for in vivo imaging are on the market and are based on different formulations. The objective of this study was to compare contrast-related and pharmacokinetic properties of a water-soluble compound containing iomeprol (Iomeron 400) and blood-pool agents (eXIA160XL, AuroVist 15 nm, and ExiTron nano 12000) for the identification of suitable in vivo vascular imaging applications.

MATERIALS AND METHODS

Forty-four healthy C57BL/6J mice were used in this study. Iomeprol was administered with a continuous infusion protocol; the other agents as a bolus. Anatomical micro-CT was applied at the head, neck, and lower hind limb before (baseline) and immediately after contrast injection, and used to quantify contrast-related properties of the agents. Dynamic micro-CT was applied at the same regions to characterize the agents pharmacokinetics.

RESULTS

All contrast media revealed safe, except for eXIA160XL, which caused death in four of eight tested animals and was therefore excluded early from the study. AuroVist 15 nm provided the highest attenuation (2.33/mm) as compared to iomeprol (1.97/mm) and ExiTron nano 12000 (1.58/mm) and a maximum temporal variation of contrast of 20% after 30 minutes, but the appearance of a dark skin staining did not allow multiple injections of the agent. Iomeprol passively diffused across capillary membranes, and after 30 minutes doubled the tissue contrast with respect to its initial levels. ExiTron nano 12000 revealed temporal variations of contrast below 10% and significantly reduced clearance rates after the third consecutive injection.

CONCLUSION

AuroVist 15 nm is best suited for anatomical investigation of the vascular network, while the high extravasation levels of iomeprol can be exploited for perfusion analysis. ExiTron nano 12000 is indicated for use in longitudinal monitoring with repeated injections.

Characterization of statistical prior image constrained compressed sensing (PICCS): II. Application to dose reduction

PURPOSE

The ionizing radiation imparted to patients during computed tomography exams is raising concerns. This paper studies the performance of a scheme called dose reduction using prior image constrained compressed sensing (DR-PICCS). The purpose of this study is to characterize the effects of a statistical model of x-ray detection in the DR-PICCS framework and its impact on spatial resolution.

METHODS

Both numerical simulations with known ground truth and in vivo animal dataset were used in this study. In numerical simulations, a phantom was simulated with Poisson noise and with varying levels of eccentricity. Both the conventional filtered backprojection (FBP) and the PICCS algorithms were used to reconstruct images. In PICCS reconstructions, the prior image was generated using two different denoising methods: a simple Gaussian blur and a more advanced diffusion filter. Due to the lack of shift-invariance in nonlinear image reconstruction such as the one studied in this paper, the concept of local spatial resolution was used to study the sharpness of a reconstructed image. Specifically, a directional metric of image sharpness, the so-called pseudopoint spread function (pseudo-PSF), was employed to investigate local spatial resolution.

RESULTS

In the numerical studies, the pseudo-PSF was reduced from twice the voxel width in the prior image down to less than 1.1 times the voxel width in DR-PICCS reconstructions when the statistical model was not included. At the same noise level, when statistical weighting was used, the pseudo-PSF width in DR-PICCS reconstructed images varied between 1.5 and 0.75 times the voxel width depending on the direction along which it was measured. However, this anisotropy was largely eliminated when the prior image was generated using diffusion filtering; the pseudo-PSF width was reduced to below one voxel width in that case. In the in vivo study, a fourfold improvement in CNR was achieved while qualitatively maintaining sharpness; images also had a qualitatively more uniform noise spatial distribution when including a statistical model.

CONCLUSIONS

DR-PICCS enables to reconstruct CT images with lower noise than FBP and the loss of spatial resolution can be mitigated to a large extent. The introduction of statistical modeling in DR-PICCS may improve some noise characteristics, but it also leads to anisotropic spatial resolution properties. A denoising method, such as the directional diffusion filtering, has been demonstrated to reduce anisotropy in spatial resolution effectively when it was combined with DR-PICCS with statistical modeling.

C-arm CT measurement of cerebral blood volume and cerebral blood flow using a novel high-speed acquisition and a single intravenous contrast injection

BACKGROUND AND PURPOSE

Assessment of perfusion parameters is important in the selection of patients who are most likely to benefit from revascularization after an acute ischemic stroke. The aim of this study was to evaluate the feasibility of measuring cerebral perfusion parameters with the use of a novel high-speed C-arm CT acquisition in conjunction with a single intravenous injection of contrast.

MATERIALS AND METHODS

Seven canines had experimentally induced focal ischemic regions confirmed by CT perfusion imaging. Four hours after ischemic injury creation, each subject underwent cerebral perfusion measurements with the use of standard perfusion CT, immediately followed by the use of C-arm CT. Cerebral blood flow and cerebral blood volume maps measured by C-arm CT were quantitatively and qualitatively compared with those measured by perfusion CT for 6 of the 7 canine subjects.

RESULTS

Results from independent observer evaluations of perfusion CT and C-arm perfusion maps show strong agreement between observers for identification of ischemic lesion location. Significant percentage agreement between observers for lesion detection and identification of perfusion mismatch between CBV and CBF maps indicate that the maps for both perfusion CT and C-arm are easy to interpret. Quantitative region of interest-based evaluation showed a strong correlation between the perfusion CT and C-arm CBV and CBF maps (R(2) = 0.68 and 0.85). C-arm measurements for both CBV and CBF were consistently overestimated when compared with perfusion CT.

CONCLUSIONS

Qualitative and quantitative measurements of CBF and CBV with the use of a C-arm CT acquisition and a single intravenous injection of contrast agent are feasible. Future improvements in flat detector technology and software algorithms probably will enable more accurate quantitative perfusion measurements with the use of C-arm CT.

Characterization of statistical prior image constrained compressed sensing. I. Applications to time-resolved contrast-enhanced CT

PURPOSE

Prior image constrained compressed sensing (PICCS) is an image reconstruction framework that takes advantage of a prior image to improve the image quality of CT reconstructions. An interesting question that remains to be investigated is whether or not the introduction of a statistical model of the photon detection in the PICCS reconstruction framework can improve the performance of the algorithm when dealing with high noise projection datasets. The goal of the research presented in this paper is to characterize the noise properties of images reconstructed using PICCS with and without statistical modeling. This paper investigates these properties in the clinical context of time-resolved contrast-enhanced CT.

METHODS

Both numerical phantom studies and an Institutional Review Board approved human subject study were used in this research. The conventional filtered backprojection (FBP), and PICCS with and without the statistical model were applied to each dataset. The prior image used in PICCS was generated by averaging over FBP reconstructions from different time frames of the time-resolved CT exam, thus reducing the noise level. Numerical studies were used to evaluate if the noise characteristics are altered for varying levels of noise, as well as for different object shapes. The dataset acquired in vivo was used to verify that the conclusions reached from numerical studies translate adequately to a clinical case. The results were analyzed using a variety of qualitative and quantitative metrics such as the universal image quality index, spatial maps of the noise standard deviations, the noise uniformity, the noise power spectrum, and the model-observer detectability.

RESULTS

The noise characteristics of PICCS were shown to depend on the noise level contained in the data, the level of eccentricity of the object, and whether or not the statistical model was applied. Most differences in the characteristics were observed in the regime of low incident x-ray fluence. No substantial difference was observed between PICCS with and without statistics in the high fluence domain. Objects with a semi-major axis ratio below 0.85 were more accurately reconstructed with lower noise using the statistical implementation. Above that range, for mostly circular objects, the PICCS implementation without the statistical model yielded more accurate images and a lower noise level. At all levels of eccentricity, the noise spatial distribution was more uniform and the model-observer detectability was greater for PICCS with the statistical model. The human subject study was consistent with the results obtained using numerical simulations.

CONCLUSIONS

For mildly eccentric objects in the low noise regime, PICCS without the noise model yielded equal or better noise level and image quality than the statistical formulation. However, in a vast majority of cases, images reconstructed using statistical PICCS have a noise power spectrum that facilitated the detection of model lesions. The inclusion of a statistical model in the PICCS framework does not always result in improved noise characteristics.

Noise spatial nonuniformity and the impact of statistical image reconstruction in CT myocardial perfusion imaging

PURPOSE

To achieve high temporal resolution in CT myocardial perfusion imaging (MPI), images are often reconstructed using filtered backprojection (FBP) algorithms from data acquired within a short-scan angular range. However, the variation in the central angle from one time frame to the next in gated short scans has been shown to create detrimental partial scan artifacts when performing quantitative MPI measurements. This study has two main purposes. (1) To demonstrate the existence of a distinct detrimental effect in short-scan FBP, i.e., the introduction of a nonuniform spatial image noise distribution; this nonuniformity can lead to unexpectedly high image noise and streaking artifacts, which may affect CT MPI quantification. (2) To demonstrate that statistical image reconstruction (SIR) algorithms can be a potential solution to address the nonuniform spatial noise distribution problem and can also lead to radiation dose reduction in the context of CT MPI.

METHODS

Projection datasets from a numerically simulated perfusion phantom and an in vivo animal myocardial perfusion CT scan were used in this study. In the numerical phantom, multiple realizations of Poisson noise were added to projection data at each time frame to investigate the spatial distribution of noise. Images from all datasets were reconstructed using both FBP and SIR reconstruction algorithms. To quantify the spatial distribution of noise, the mean and standard deviation were measured in several regions of interest (ROIs) and analyzed across time frames. In the in vivo study, two low-dose scans at tube currents of 25 and 50 mA were reconstructed using FBP and SIR. Quantitative perfusion metrics, namely, the normalized upslope (NUS), myocardial blood volume (MBV), and first moment transit time (FMT), were measured for two ROIs and compared to reference values obtained from a high-dose scan performed at 500 mA.

RESULTS

Images reconstructed using FBP showed a highly nonuniform spatial distribution of noise. This spatial nonuniformity led to large fluctuations in the temporal direction. In the numerical phantom study, the level of noise was shown to vary by as much as 87% within a given image, and as much as 110% between different time frames for a ROI far from isocenter. The spatially nonuniform noise pattern was shown to correlate with the source trajectory and the object structure. In contrast, images reconstructed using SIR showed a highly uniform spatial distribution of noise, leading to smaller unexpected noise fluctuations in the temporal direction when a short scan angular range was used. In the numerical phantom study, the noise varied by less than 37% within a given image, and by less than 20% between different time frames. Also, the noise standard deviation in SIR images was on average half of that of FBP images. In the in vivo studies, the deviation observed between quantitative perfusion metrics measured from low-dose scans and high-dose scans was mitigated when SIR was used instead of FBP to reconstruct images.

CONCLUSIONS

(1) Images reconstructed using FBP suffered from nonuniform spatial noise levels. This nonuniformity is another manifestation of the detrimental effects caused by short-scan reconstruction in CT MPI. (2) Images reconstructed using SIR had a much lower and more uniform noise level and thus can be used as a potential solution to address the FBP nonuniformity. (3) Given the improvement in the accuracy of the perfusion metrics when using SIR, it may be desirable to use a statistical reconstruction framework to perform low-dose dynamic CT MPI.

Improved compressed sensing-based cone-beam CT reconstruction using adaptive prior image constraints

Volumetric cone-beam CT (CBCT) images are acquired repeatedly during a course of radiation therapy and a natural question to ask is whether CBCT images obtained earlier in the process can be utilized as prior knowledge to reduce patient imaging dose in subsequent scans. The purpose of this work is to develop an adaptive prior image constrained compressed sensing (APICCS) method to solve this problem. Reconstructed images using full projections are taken on the first day of radiation therapy treatment and are used as prior images. The subsequent scans are acquired using a protocol of sparse projections. In the proposed APICCS algorithm, the prior images are utilized as an initial guess and are incorporated into the objective function in the compressed sensing (CS)-based iterative reconstruction process. Furthermore, the prior information is employed to detect any possible mismatched regions between the prior and current images for improved reconstruction. For this purpose, the prior images and the reconstructed images are classified into three anatomical regions: air, soft tissue and bone. Mismatched regions are identified by local differences of the corresponding groups in the two classified sets of images. A distance transformation is then introduced to convert the information into an adaptive voxel-dependent relaxation map. In constructing the relaxation map, the matched regions (unchanged anatomy) between the prior and current images are assigned with smaller weight values, which are translated into less influence on the CS iterative reconstruction process. On the other hand, the mismatched regions (changed anatomy) are associated with larger values and the regions are updated more by the new projection data, thus avoiding any possible adverse effects of prior images. The APICCS approach was systematically assessed by using patient data acquired under standard and low-dose protocols for qualitative and quantitative comparisons. The APICCS method provides an effective way for us to enhance the image quality at the matched regions between the prior and current images compared to the existing PICCS algorithm. Compared to the current CBCT imaging protocols, the APICCS algorithm allows an imaging dose reduction of 10-40 times due to the greatly reduced number of projections and lower x-ray tube current level coming from the low-dose protocol.

Dynamic cone beam CT angiography of carotid and cerebral arteries using canine model

PURPOSE

This research is designed to develop and evaluate a flat-panel detector-based dynamic cone beam CT system for dynamic angiography imaging, which is able to provide both dynamic functional information and dynamic anatomic information from one multirevolution cone beam CT scan.

METHODS

A dynamic cone beam CT scan acquired projections over four revolutions within a time window of 40 s after contrast agent injection through a femoral vein to cover the entire wash-in and wash-out phases. A dynamic cone beam CT reconstruction algorithm was utilized and a novel recovery method was developed to correct the time-enhancement curve of contrast flow. From the same data set, both projection-based subtraction and reconstruction-based subtraction approaches were utilized and compared to remove the background tissues and visualize the 3D vascular structure to provide the dynamic anatomic information.

RESULTS

Through computer simulations, the new recovery algorithm for dynamic time-enhancement curves was optimized and showed excellent accuracy to recover the actual contrast flow. Canine model experiments also indicated that the recovered time-enhancement curves from dynamic cone beam CT imaging agreed well with that of an IV-digital subtraction angiography (DSA) study. The dynamic vascular structures reconstructed using both projection-based subtraction and reconstruction-based subtraction were almost identical as the differences between them were comparable to the background noise level. At the enhancement peak, all the major carotid and cerebral arteries and the Circle of Willis could be clearly observed.

CONCLUSIONS

The proposed dynamic cone beam CT approach can accurately recover the actual contrast flow, and dynamic anatomic imaging can be obtained with high isotropic 3D resolution. This approach is promising for diagnosis and treatment planning of vascular diseases and strokes.

Prior image constrained compressed sensing: implementation and performance evaluation

PURPOSE

Prior image constrained compressed sensing (PICCS) is an image reconstruction framework which incorporates an often available prior image into the compressed sensing objective function. The images are reconstructed using an optimization procedure. In this paper, several alternative unconstrained minimization methods are used to implement PICCS. The purpose is to study and compare the performance of each implementation, as well as to evaluate the performance of the PICCS objective function with respect to image quality.

METHODS

Six different minimization methods are investigated with respect to convergence speed and reconstruction accuracy. These minimization methods include the steepest descent (SD) method and the conjugate gradient (CG) method. These algorithms require a line search to be performed. Thus, for each minimization algorithm, two line searching algorithms are evaluated: a backtracking (BT) line search and a fast Newton-Raphson (NR) line search. The relative root mean square error is used to evaluate the reconstruction accuracy. The algorithm that offers the best convergence speed is used to study the performance of PICCS with respect to the prior image parameter α and the data consistency parameter λ. PICCS is studied in terms of reconstruction accuracy, low-contrast spatial resolution, and noise characteristics. A numerical phantom was simulated and an animal model was scanned using a multirow detector computed tomography (CT) scanner to yield the projection datasets used in this study.

RESULTS

For λ within a broad range, the CG method with Fletcher-Reeves formula and NR line search offers the fastest convergence for an equal level of reconstruction accuracy. Using this minimization method, the reconstruction accuracy of PICCS was studied with respect to variations in α and λ. When the number of view angles is varied between 107, 80, 64, 40, 20, and 16, the relative root mean square error reaches a minimum value for α ≈ 0.5. For values of α near the optimal value, the spatial resolution of the reconstructed image remains relatively constant and the noise texture is very similar to that of the prior image, which was reconstructed using the filtered backprojection (FBP) algorithm.

CONCLUSIONS

Regarding the performance of the minimization methods, the nonlinear CG method with NR line search yields the best convergence speed. Regarding the performance of the PICCS image reconstruction, three main conclusions can be reached. (1) The performance of PICCS is optimal when the weighting parameter of the prior image parameter is selected to be near α = 0.5. (2) The spatial resolution measured for static objects in images reconstructed using PICCS from undersampled datasets is not degraded with respect to the fully-sampled reconstruction for α near its optimal value. (3) The noise texture of PICCS reconstructions is similar to that of the prior image, which was reconstructed using the conventional FBP method.

Interventional 4-D C-arm CT perfusion imaging using interleaved scanning and partial reconstruction interpolation

Tissue perfusion measurement during catheter-guided stroke treatment in the interventional suite is currently not possible. In this work, we present a novel approach that uses a C-arm angiography system capable of computed tomography (CT)-like imaging (C-arm CT) for this purpose. With C-arm CT one reconstructed volume can be obtained every 4-6 s which makes it challenging to measure the flow of an injected contrast bolus. We have developed an interleaved scanning (IS) protocol that uses several scan sequences to increase temporal sampling. Using a dedicated 4-D reconstruction approach based on partial reconstruction interpolation (PRI) we can optimally process our data. We evaluated our combined approach (IS-PRI) with simulations and a study in five healthy pigs. In our simulations, the cerebral blood flow values (unit: ml/100 g/min) were 60 (healthy tissue) and 20 (pathological tissue). For one scan sequence the values were estimated with standard deviations of 14.3 and 2.9, respectively. For two interleaved sequences the standard deviations decreased to 3.6 and 1.5, respectively. We used perfusion CT to validate the in vivo results. With two interleaved sequences we achieved promising correlations ranging from r=0.63 to r=0.94. The results suggest that C-arm CT tissue perfusion imaging is feasible with two interleaved scan sequences.

Dynamic contrast-enhanced micro-CT on mice with mammary carcinoma for the assessment of antiangiogenic therapy response

OBJECTIVE

To evaluate the potential of in vivo dynamic contrast-enhanced micro-computed tomography (DCE micro-CT) for the assessment of antiangiogenic drug therapy response of mice with mammary carcinoma.

METHODS

20 female mice with implanted MCF7 tumours were split into control group and therapy group treated with a known effective antiangiogenic drug. All mice underwent DCE micro-CT for the 3D analysis of functional parameters (relative blood volume [rBV], vascular permeability [K], area under the time-enhancement curve [AUC]) and morphology. All parameters were determined for total, peripheral and central tumour volumes of interest (VOIs). Immunohistochemistry was performed to characterise tumour vascularisation. 3D dose distributions were determined.

RESULTS

The mean AUCs were significantly lower in therapy with P values of 0.012, 0.007 and 0.023 for total, peripheral and central tumour VOIs. K and rBV showed significant differences for the peripheral (P(per)(K) = 0.032, P(per) (rBV) = 0.029), but not for the total and central tumour VOIs (P(total)(K) = 0.108, P(central)(K) = 0.246, P(total) (rBV) = 0.093, P(central) (rBV) = 0.136). Mean tumour volume was significantly smaller in therapy (P (in vivo) = 0.001, P (ex vivo) = 0.005). Histology revealed greater vascularisation in the controls and central tumour necrosis. Doses ranged from 150 to 300 mGy.

CONCLUSIONS

This study indicates the great potential of DCE micro-CT for early in vivo assessment of antiangiogenic drug therapy response.

KEY POINTS

Dynamic contrast enhanced micro-CT (computed tomography) is a new experimental laboratory technique. DCE micro-CT allows early in vivo assessment of antiangiogenic drug therapy response. Pharmaceutical drugs can be tested before translation to clinical practice. Both morphological and functional parameters can be obtained using DCE micro-CT. Antiangiogenic effects can be visualised with DCE micro-CT.

Nonconvex prior image constrained compressed sensing (NCPICCS): theory and simulations on perfusion CT

PURPOSE

To present and evaluate a new image reconstruction method for dynamic CT based on a nonconvex prior image constrained compressed sensing (NCPICCS) algorithm. The authors systematically compared the undersampling potential, functional information recovery, and solution convergence speed of four compressed sensing (CS) based image reconstruction methods using perfusion CT data: Standard l1-based CS, nonconvex CS (NCCS), and l1-based and nonconvex CS, including an additional constraint based on a prior image (PICCS and NCPICCS, respectively).

METHODS

The Shepp-Logan phantom was modified such that its uppermost ellipses changed attenuation through time, simulating both an arterial input function (AIF) and a homogeneous tissue perfusion region. Data were simulated with and without Poisson noise added to the projection data and subsequently reconstructed with all four CS-based methods at four levels of undersampling: 20, 12, 6, and 4 projections. Root mean squared (RMS) error of reconstructed images and recovered time attenuation curves (TACs) were assessed as well as convergence speed. The performance of both PICCS and NCPICCS methods were also evaluated using a kidney perfusion animal experiment data set.

RESULTS

All four CS-based methods were able to reconstruct the phantoms with 20 projections, with similar results on the RMS error of the recovered TACs. NCCS allowed accurate reconstructions with as few as 12 projections, PICCS with as few as six projections, and NCPICCS with as few as four projections. These results were consistent for noise-free and noisy data. NCPICCS required the fewest iterations to converge across all simulation conditions, followed by PICCS, NCCS, and then CS. On animal data, at the lowest level of undersampling tested (16 projections), the image quality of NCPICCS was better than PICCS with fewer streaking artifacts, while the TAC accuracy on the selected region of interest was comparable.

CONCLUSIONS

The authors have presented a novel method for image reconstruction using highly undersampled dynamic CT data. The NCPICCS method takes advantage of the information provided by a prior image, as in PICCS, but employs a more general nonconvex sparsity measure [such as the l(p)-norm (0 < p < or = 1)] rather than the conventional convex l1-norm. Despite the lack of guarantees of a globally optimal solution, the proposed nonconvex extension of PICCS consistently allowed for image reconstruction from fewer samples than the analogous l1-based PICCS method. Both nonconvex sparsity measures as well as prior image information (when available) significantly reduced the number of iterations required for convergence, potentially providing computational advantages for practical implementation of CS-based image reconstruction techniques.

4D micro-CT for cardiac and perfusion applications with view under sampling

Micro-CT is commonly used in preclinical studies to provide anatomical information. There is growing interest in obtaining functional measurements from 4D micro-CT. We report here strategies for 4D micro-CT with a focus on two applications: (i) cardiac imaging based on retrospective gating and (ii) pulmonary perfusion using multiple contrast injections/rotations paradigm. A dual source micro-CT system is used for image acquisition with a sampling rate of 20 projections per second. The cardiac micro-CT protocol involves the use of a liposomal blood pool contrast agent. Fast scanning of free breathing mice is achieved using retrospective gating. The ECG and respiratory signals are used to sort projections into ten cardiac phases. The pulmonary perfusion protocol uses a conventional contrast agent (Isovue 370) delivered by a micro-injector in four injections separated by 2 min intervals to allow for clearance. Each injection is synchronized with the rotation of the animal, and each of the four rotations is started with an angular offset of 22.5 from the starting angle of the previous rotation. Both cardiac and perfusion protocols result in an irregular angular distribution of projections that causes significant streaking artifacts in reconstructions when using traditional filtered backprojection (FBP) algorithms. The reconstruction involves the use of the point spread function of the micro-CT system for each time point, and the analysis of the distribution of the reconstructed data in the Fourier domain. This enables us to correct for angular inconsistencies via deconvolution and identify regions where data is missing. The missing regions are filled with data from a high quality but temporally averaged prior image reconstructed with all available projections. Simulations indicate that deconvolution successfully removes the streaking artifacts while preserving temporal information. 4D cardiac micro-CT in a mouse was performed with adequate image quality at isotropic voxel size of 88 µm and 10 ms temporal resolution. 4D pulmonary perfusion images were obtained in a mouse at 176 µm and 687 ms temporal resolution. Compared with FBP reconstruction, the streak reduction ratio is 70% and the contrast to noise ratio is 2.5 times greater in the deconvolved images. The radiation dose associated with the proposed methods is in the range of a typical micro-CT dose (0.17 Gy for the cardiac study and 0.21 Gy for the perfusion study). The low dose 4D micro-CT imaging presented here can be applied in high-throughput longitudinal studies in a wide range of applications, including drug safety and cardiopulmonary phenotyping.

64-slice CT perfusion imaging of pancreatic adenocarcinoma and mass-forming chronic pancreatitis

RATIONALE AND OBJECTIVES

To investigate 64 computed tomography (CT) perfusion imaging features of patients with pancreatic cancer and mass-forming chronic pancreatitis.

MATERIALS AND METHODS

Between January 2003 and April 2010, 234 patients with pancreatic mass underwent 64-CT perfusion imaging. Among them, the histopathological results of 64 patients were proven to be pancreatic adenocarcinoma and 15 patients were proven to be mass-forming chronic pancreatitis. Additionally, CT perfusion imaging was performed in 33 healthy volunteers served as controls. The slice data were processed using CT perfusion software. Perfusion parameters including time density curve, blood flow, blood volume, permeability, peak enhancement, and time to peak were recorded.

RESULTS

Blood flow was 77% lower in patients with pancreatic adenocarcinoma than in controls, 48% lower in patients with mass-forming chronic pancreatitis than in controls, and 56% lower in patients with pancreatic adenocarcinoma than with mass-forming chronic pancreatitis (P < .016). Blood volume was 65% lower in pancreatic adenocarcinoma than in controls, 27% lower in mass-forming chronic pancreatitis than in controls, and 53% lower in cancer than mass-forming chronic pancreatitis (P < .016). Permeability was 559% higher in pancreatic adenocarcinoma than in controls, 821% higher in mass-forming chronic pancreatitis than in controls, and 28% lower in cancer than mass-forming chronic pancreatitis (P < .016). Peak enhancement was 27% lower and time to peak 23% longer in pancreatic adenocarcinoma than mass-forming chronic pancreatitis (P < .016). Time-density curve showed the peak of mass-forming chronic pancreatitis is earlier and higher than that of pancreatic adenocarcinoma, and the peak of mass-forming chronic pancreatitis is later and lower than that of controls.

CONCLUSION

CT perfusion imaging can provide additional quantitative hemodynamic information of pancreatic adenocarcinoma and mass-forming chronic pancreatitis.