Renal perfusion and hemodynamics: accurate in vivo determination at CT with a 10-fold decrease in radiation dose and HYPR noise reduction

Reliable Assessment of Swine Renal Fibrosis Using Quantitative Magnetization Transfer Imaging

OBJECTIVES

Quantitative magnetization transfer (qMT) is useful for measurement of murine renal fibrosis at high and ultrahigh field strengths. However, its utility at clinical field strengths and in human-like kidneys remains unknown. We tested the hypothesis that qMT would successfully detect fibrosis in swine kidneys with unilateral renal artery stenosis (RAS) at 3.0 T.

METHODS

The qMT protocol is composed of MT scans with variable flip angles and offset frequencies, and of B0, B1, and T1 mapping. Pigs were scanned 10 weeks after RAS or control. A 2-pool model was used to fit the bound pool fraction f of the renal cortex (CO) and outer medulla (OM). Then qMT-derived f in 5 normal and 10 RAS pigs was compared with histological fibrosis determined using Masson's trichrome staining and to renal perfusion assessed with computed tomography.

RESULTS

The qMT 2-pool model provided accurate fittings of data collected on swine kidneys. Stenotic kidneys showed significantly elevated f in both the CO (9.8% ± 2.7% vs 6.4% ± 0.9%, P = 0.002) and OM (7.6% ± 2.2% vs 4.7% ± 1.1%, P = 0.002), as compared with normal kidneys. Histology-measured renal fibrosis and qMT-derived f correlated directly in both the cortex (Pearson correlation coefficient r = 0.93, P < 0.001) and OM (r = 0.84, P = 0.002), and inversely with stenotic kidney perfusion (r = 0.85, P = 0.002).

CONCLUSIONS

This study demonstrates the feasibility of qMT for measuring fibrosis in human-like swine kidneys, and the association between tissue macromolecule content and renal perfusion. Therefore, qMT may be useful as a tool for noninvasive assessment of renal fibrosis in subjects with RAS at clinical field strengths.

Cubic-Spline Interpolation for Sparse-View CT Image Reconstruction With Filtered Backprojection in Dynamic Myocardial Perfusion Imaging

We investigated a projection interpolation method for reconstructing dynamic contrast-enhanced (DCE) heart images from undersampled x-ray projections with filtered backprojecton (FBP). This method may facilitate the application of sparse-view dynamic acquisition for ultralow-dose quantitative computed tomography (CT) myocardial perfusion (MP) imaging. We conducted CT perfusion studies on 5 pigs with a standard full-view acquisition protocol (984 projections). We reconstructed DCE heart images with FBP from all and a quarter of the measured projections evenly distributed over 360°. We interpolated the sparse-view (quarter) projections to a full-view setting using a cubic-spline interpolation method before applying FBP to reconstruct the DCE heart images (synthesized full-view). To generate MP maps, we used 3 sets of DCE heart images, and compared mean MP values and biases among the 3 protocols. Compared with synthesized full-view DCE images, sparse-view DCE images were more affected by streak artifacts arising from projection undersampling. Relative to the full-view protocol, mean bias in MP measurement associated with the sparse-view protocol was 10.0 mL/min/100 g (95%CI: −8.9 to 28.9), which was >3 times higher than that associated with the synthesized full-view protocol (3.3 mL/min/100 g, 95% CI: −6.7 to 13.2). The cubic-spline-view interpolation method improved MP measurement from DCE heart images reconstructed from only a quarter of the full projection set. This method can be used with the industry-standard FBP algorithm to reconstruct DCE images of the heart, and it can reduce the radiation dose of a whole-heart quantitative CT MP study to <2 mSv (at 8-cm coverage).

Renal Adiposity Does not Preclude Quantitative Assessment of Renal Function Using Dual-Energy Multidetector CT in Mildly Obese Human Subjects

RATIONALE AND OBJECTIVES

Multidetector computed tomography (MDCT) is useful for measuring in the research setting single-kidney perfusion and function using iodinated contrast time-attenuation curves. Obesity promotes deposition of intrarenal fat, which might decrease tissue attenuation and thereby interfere with quantification of renal function using MDCT. The purpose of this study was to test the hypothesis that background subtraction adequately accounts for intrarenal fat deposition in mildly obese human subjects during renal contrast enhanced dynamic CT.

MATERIALS AND METHODS

We prospectively recruited seventeen human subjects stratified as lean or mildly obese based on body mass index below or over 30 kg/m², respectively. Renal perfusion was quantified from CT-derived indicator-dilution curves after background subtraction. Dual-energy MDCT images were postprocessed to generate iodine and virtual-noncontrast datasets, and the ratios between kidney/aorta CT numbers and iodine values calculated as surrogates of renal function.

RESULTS

Subcutaneous adipose tissue was increased in obese subjects. Virtual-noncontrast maps revealed in obese patients a decrease in basal cortical and medullary attenuation. Overall, basal attenuation inversely correlated with body mass index, in line with renal fat deposition. Contrarily, the kidney/aorta CT attenuation (after background subtraction) and kidney/aorta iodine ratios were similar between lean and obese subjects and correlated directly. These observations show that following background subtraction, the CT number reliably reflects basal tissue attenuation.

CONCLUSION

Therefore, our findings support our hypothesis that background subtraction enables reliable assessment of kidney function in mildly obese subjects using MDCT, despite decreased basal attenuation due to renal adiposity.

A robust noise reduction technique for time resolved CT

PURPOSE

To develop a noise reduction method for time resolved CT data, especially those with significant patient motion.

METHODS

PArtial TEmporal Nonlocal (PATEN) means is a technique that uses the redundant information in time-resolved CT data to achieve noise reduction. In this method, partial temporal profiles are used to determine the similarity (or weight) between pixels, and the similarity search makes use of both spatial and temporal information, providing robustness to patient motion. The performance of the PATEN filter was qualitatively and quantitatively evaluated with nine cardiac CT patient data sets and five CT brain perfusion patient data sets. In cardiac CT, PATEN was applied to reduce noise primarily in the reduced-dose phases created with electrocardiographic (ECG) pulsing. CT number accuracy and noise reduction were evaluated in both full-dose phases and reduced-dose phases between filtered backprojection images and PATEN filtered images. In CT brain perfusion, simulated quarter dose data were obtained by adding noise to the raw data of a routine dose scan. PATEN was applied to the simulated low-dose images. Image noise, time-intensity profile accuracy, and perfusion parameter maps were compared among low-dose, low-dose+PATEN filter, and full-dose images. The noise reduction performance of PATEN was compared to a previously proposed noise reduction method, time-intensity profile similarity (TIPS) bilateral filtering.

RESULTS

In 4D cardiac CT, after PATEN filtering, the image noise in the reduced-dose phases was greatly reduced, making anatomical structures easier to identify. The mean decreases in noise values between the original and PATEN images were 11.0% and 53.8% for the full and reduced-dose phases of the cardiac cycle, respectively. TIPS could not achieve effective noise reduction. In CT brain perfusion, PATEN achieved a 55.8%-66.3% decrease in image noise in the low-dose images. The contrast to noise ratio in the axial images was increased and was comparable to the full-dose images. Differentiation of anatomical structure in the PATEN images and corresponding quantitative perfusion parameter maps were preferred by two neuroradiologists compared to the simulated low-dose and TIPS results. The mean perfusion parameters calculated from the PATEN images agreed with those determined from full-dose data to within 12% and 20% for normal and diseased regions.

CONCLUSIONS

In ECG-gated cardiac CT, where the dose had already been reduced by a factor of 5 in the reduced-dose phases, PATEN achieved a 53.8% noise reduction, which decreased the noise level in the reduced-dose phases close to that of the full-dose phases. In CT brain perfusion, a fourfold dose reduction was demonstrated to be achievable by PATEN filtering, which improved quantitative perfusion analysis. PATEN can be used to effectively reduce image noise to improve image quality, even when significant motion occurred between temporal samples.

Improved accuracy of quantitative parameter estimates in dynamic contrast-enhanced CT study with low temporal resolution

PURPOSE

A previously proposed method to reduce radiation dose to patient in dynamic contrast-enhanced (DCE) CT is enhanced by principal component analysis (PCA) filtering which improves the signal-to-noise ratio (SNR) of time-concentration curves in the DCE-CT study. The efficacy of the combined method to maintain the accuracy of kinetic parameter estimates at low temporal resolution is investigated with pixel-by-pixel kinetic analysis of DCE-CT data.

METHODS

The method is based on DCE-CT scanning performed with low temporal resolution to reduce the radiation dose to the patient. The arterial input function (AIF) with high temporal resolution can be generated with a coarsely sampled AIF through a previously published method of AIF estimation. To increase the SNR of time-concentration curves (tissue curves), first, a region-of-interest is segmented into squares composed of 3 × 3 pixels in size. Subsequently, the PCA filtering combined with a fraction of residual information criterion is applied to all the segmented squares for further improvement of their SNRs. The proposed method was applied to each DCE-CT data set of a cohort of 14 patients at varying levels of down-sampling. The kinetic analyses using the modified Tofts' model and singular value decomposition method, then, were carried out for each of the down-sampling schemes between the intervals from 2 to 15 s. The results were compared with analyses done with the measured data in high temporal resolution (i.e., original scanning frequency) as the reference.

RESULTS

The patients' AIFs were estimated to high accuracy based on the 11 orthonormal bases of arterial impulse responses established in the previous paper. In addition, noise in the images was effectively reduced by using five principal components of the tissue curves for filtering. Kinetic analyses using the proposed method showed superior results compared to those with down-sampling alone; they were able to maintain the accuracy in the quantitative histogram parameters of volume transfer constant [standard deviation (SD), 98th percentile, and range], rate constant (SD), blood volume fraction (mean, SD, 98th percentile, and range), and blood flow (mean, SD, median, 98th percentile, and range) for sampling intervals between 10 and 15 s.

CONCLUSIONS

The proposed method of PCA filtering combined with the AIF estimation technique allows low frequency scanning for DCE-CT study to reduce patient radiation dose. The results indicate that the method is useful in pixel-by-pixel kinetic analysis of DCE-CT data for patients with cervical cancer.

Assessment of renal function in patients with unilateral ureteral obstruction using whole-organ perfusion imaging with 320-detector row computed tomography

Background

Obstructed nephropathy is a common complication of several disease processes. Accurate evaluation of the functional status of the obstructed kidney is important to achieve a good outcome. The purpose of this study was to investigate renal cortical and medullary perfusion changes associated with unilateral ureteral obstruction (UUO) using whole-organ perfusion imaging with 320-detector row computed tomography (CT).

Methodology/Principle Findings

Sixty-four patients with UUO underwent whole-organ CT perfusion imaging. Patients were divided into 3 groups, mild, moderate, and severe, based on hydronephrosis severity. Twenty sex- and age-matched patients without renal disease, who referred to abdominal CT, were chosen as control subjects. Mean cortical and medullary perfusion parameters of obstructed and contralateral kidneys were compared, and mean perfusion ratios between obstructed and contralateral kidneys were calculated and compared. Mean cortical or medullary blood flow (BF) and blood volume (BV) of the obstructed kidneys in the moderate UUO and BF, BV, and clearance (CL) in the severe UUO were significantly lower than those of the contralateral kidneys (p < 0.05). The mean cortical or medullary BF of the obstructed kidney in the moderate UUO, and BF, BV, and CL in the severe UUO were significantly lower than those of the kidneys in control subjects (p < 0.05). Mean cortical or medullary BF of the non-obstructed kidneys in the severe UUO were statistically greater than that of normal kidneys in control subjects (p < 0.05). An inverse correlation was observed between cortical and medullary perfusion ratios and grades of hydronephosis (p < 0.01).

Conclusions/Significance

Perfusion measurements of the whole kidney can be obtained with 320-detector row CT, and estimated perfusion ratios have potential for quantitatively evaluating UUO renal injury grades.

Assessment of renal function; clearance, the renal microcirculation, renal blood flow, and metabolic balance

Historically, tools to assess renal function have been developed to investigate the physiology of the kidney in an experimental setting, and certain of these techniques have utility in evaluating renal function in the clinical setting. The following work will survey a spectrum of these tools, their applications and limitations in four general sections. The first is clearance, including evaluation of exogenous and endogenous markers for determining glomerular filtration rate, the adaptation of estimated glomerular filtration rate in the clinical arena, and additional clearance techniques to assess various other parameters of renal function. The second section deals with in vivo and in vitro approaches to the study of the renal microvasculature. This section surveys a number of experimental techniques including corticotomy, the hydronephrotic kidney, vascular casting, intravital charge coupled device videomicroscopy, multiphoton fluorescent microscopy, synchrotron-based angiography, laser speckle contrast imaging, isolated renal microvessels, and the perfused juxtamedullary nephron microvasculature. The third section addresses in vivo and in vitro approaches to the study of renal blood flow. These include ultrasonic flowmetry, laser-Doppler flowmetry, magnetic resonance imaging (MRI), phase contrast MRI, cine phase contrast MRI, dynamic contrast-enhanced MRI, blood oxygen level dependent MRI, arterial spin labeling MRI, x-ray computed tomography, and positron emission tomography. The final section addresses the methodologies of metabolic balance studies. These are described for humans, large experimental animals as well as for rodents. Overall, the various in vitro and in vivo topics and applications to evaluate renal function should provide a guide for the investigator or physician to understand and to implement the techniques in the laboratory or clinic setting.

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.

Reduction of image noise in low tube current dynamic CT myocardial perfusion imaging using HYPR processing: a time-attenuation curve analysis

PURPOSE

This study describes a HighlY constrained backPRojection (HYPR) image processing method for the reduction of image noise in low tube current time-resolved CT myocardial perfusion scans. The effect of this method on myocardial time-attenuation curve noise and fidelity is evaluated in an animal model, using varying levels of tube current.

METHODS

CT perfusion scans of four healthy pigs (42-59 kg) were acquired at 500, 250, 100, 50, 25, and 10 mA on a 64-slice scanner (4 cm axial coverage, 120 kV, 0.4 s∕rotation, 50 s scan duration). For each scan a sequence of ECG-gated images centered on 75% R-R was reconstructed using short-scan filtered back projection (FBP). HYPR processing was applied to the scans acquired at less than 500 mA using parameters designed to maintain the voxel noise level in the 500-mA FBP images. The processing method generates a series of composite images by averaging over a sliding time window and then multiplies the composite images by weighting images to restore temporal fidelity to the image sequence. HYPR voxel noise relative to FBP noise was measured in AHA myocardial segment numbers 1, 5, 6, and 7 at each mA. To quantify the agreement between HYPR and FBP time-attenuation curves (TACs), Bland-Altman analysis was performed on TACs measured in full myocardial segments. The relative degree of TAC fluctuation in smaller subvolumes was quantified by calculating the root mean square deviation of a TAC about the gamma variate curve fit to the TAC data.

RESULTS

HYPR image sequences were produced using 2, 7, and 20 beat composite windows for the 250, 100, and 50 mA scans, respectively. At 25 and 10 mA, all available beats were used in the composite (41-60; average 50). A 7-voxel-wide 3D cubic filter kernel was used to form weighting images. The average ratio of HYPR voxel noise to 500-mA FBP voxel noise was 1.06, 1.10, 0.97, 1.11, and 2.15 for HYPR scans at 250, 100, 50, 25, and 10 mA. The average limits-of-agreement between HYPR and FBP TAC values measured 0.02+∕-0.91, 0.04+∕-1.92, 0.19+∕-1.59, 1.13+∕-4.22, and 1.07+∕-6.37 HU (mean difference +∕-1.96 SD). The HYPR image subvolume that yielded a fixed level of TAC fluctuations was smaller, on average, than the FBP subvolume determined at the same mA.

CONCLUSIONS

HYPR processing is a feasible method for generating low noise myocardial perfusion data from a low-mA time-resolved CT myocardial perfusion scan. The method is applicable to current clinical scanners and uses conventional image reconstructions as input data.

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.

Diagnosis of hepatocellular carcinoma: newer radiological tools

With the recent dramatic advances in diagnostic modalities, the diagnosis of hepatocellular carcinoma (HCC) is primarily based on imaging. Ultrasound (US) plays a crucial role in HCC surveillance. Dynamic multiphasic multidetector-row CT (MDCT) and magnetic resonance imaging (MRI) are the standard diagnostic methods for the noninvasive diagnosis of HCC, which can be made based on hemodynamic features (arterial enhancement and delayed washout). The technical development of MDCT and MRI has made possible the fast scanning with better image quality and resolution, which enables an accurate CT hemodynamic evaluation of hepatocellular tumor, as well as the application of perfusion CT and MRI in clinical practice. Perfusion CT and MRI can measure perfusion parameters of tumor quantitatively and can be used for treatment response assessment to anti-vascular agents. Besides assessing the hemodynamic or perfusion features of HCC, new advances in MRI can provide a cellular information of HCC. Liver-specific hepatobiliary contrast agents, such as gadoxetic acid, give information regarding hepatocellular function or defect of the lesion, which improves lesion detection and characterization. Diffusion-weighted imaging (DWI) of the liver provides cellular information of HCC and also has broadened its role in lesion detection, lesion characterization, and treatment response assessment to chemotherapeutic agents. In this article, we provide an overview of the state-of-the art imaging techniques of the liver and their clinical role in management of HCC.

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.

A strategy to decrease partial scan reconstruction artifacts in myocardial perfusion CT: phantom and in vivo evaluation

PURPOSE

Partial scan reconstruction (PSR) artifacts are present in myocardial perfusion imaging using dynamic multidetector computed tomography (MDCT). PSR artifacts appear as temporal CT number variations due to inconsistencies in the angular data range used to reconstruct images and compromise the quantitative value of myocardial perfusion when using MDCT. The purpose of this work is to present and evaluate a technique termed targeted spatial frequency filtration (TSFF) to reduce CT number variations due to PSR when applied to myocardial perfusion imaging using MDCT.

METHODS

The TSFF algorithm requires acquiring enough X-ray projections to reconstruct both partial (π + fan angle α) and full (2π) scans. Then, using spatial linear filters, the TSFF-corrected image data are created by superimposing the low spatial frequency content of the full scan reconstruction (containing no PSR artifacts, but having low spatial resolution and worse temporal resolution) with the high spatial frequency content of the partial scan reconstruction (containing high spatial frequencies and better temporal resolution). The TSFF method was evaluated both in a static anthropomorphic thoracic phantom and using an in vivo porcine model and compared with a previously validated reference standard technique that avoids PSR artifacts by pacing the animal heart in synchrony with the gantry rotation. CT number variations were quantified by measuring the range and standard deviation of CT numbers in selected regions of interest (ROIs) over time. Myocardial perfusion parameters such as blood volume (BV), mean transit time (MTT), and blood flow (BF) were quantified and compared in the in vivo study.

RESULTS

Phantom experiments demonstrated that TSFF reduced PSR artifacts by as much as tenfold, depending on the location of the ROI. For the in vivo experiments, the TSFF-corrected data showed two- to threefold decrease in CT number variations. Also, after TSFF, the perfusion parameters had an average difference of 13.1% (range 4.5%-25.6%) relative to the reference method, in contrast to an average difference of 31.8% (range 0.3%-54.0%) between the non-TSFF processed data with the reference method.

CONCLUSIONS

TSFF demonstrated consistent reduction in CT number variations due to PSR using controlled phantom and in vivo experiments. TSFF-corrected data provided quantitative measures of perfusion (BV, MTT, and BF) with better agreement to a reference method compared to noncorrected data. Practical implementation of TSFF is expected to incur in an additional radiation exposure of 14%, when tube current is modulated to 20% of its maximum, to complete the needed full scan reconstruction.

Porcine ex vivo liver phantom for dynamic contrast-enhanced computed tomography: development and initial results

OBJECTIVES

: To demonstrate the feasibility of developing a fixed, dual-input, biologic liver phantom for dynamic contrast-enhanced computed tomography (CT) imaging and to report initial results of use of the phantom for quantitative CT perfusion imaging.

MATERIALS AND METHODS

: Porcine livers were obtained from completed surgical studies and perfused with saline and fixative. The phantom was placed in a body-shaped, CT-compatible acrylic container and connected to a perfusion circuit fitted with a contrast injection port. Flow-controlled contrast-enhanced imaging experiments were performed using 128-slice and 64-slice dual-source multidetector CT scanners. CT angiography protocols were used to obtain portal venous and hepatic arterial vascular enhancement, reproduced over a period of 4 to 6 months. CT perfusion protocols were used at different input flow rates to correlate input flow with calculated tissue perfusion, to test reproducibility, and to determine the feasibility of simultaneous dual-input liver perfusion. Histologic analysis of the liver phantom was also performed.

RESULTS

: CT angiogram 3-dimensional reconstructions demonstrated homogenous tertiary and quaternary branching of the portal venous system to the periphery of all lobes of the liver as well as enhancement of the hepatic arterial system to all lobes of the liver and gallbladder throughout the study period. For perfusion CT, the correlation between the calculated mean tissue perfusion in a volume of interest and input pump flow rate was excellent (R = 0.996) and color blood flow maps demonstrated variations in regional perfusion in a narrow range. Repeat perfusion CT experiments demonstrated reproducible time-attenuation curves, and dual-input perfusion CT experiments demonstrated that simultaneous dual input liver perfusion is feasible. Histologic analysis demonstrated that the hepatic microvasculature and architecture appeared intact and well preserved at the completion of 4 to 6 months of laboratory experiments and contrast-enhanced imaging.

CONCLUSIONS

: We have demonstrated successful development of a porcine liver phantom using a flow-controlled extracorporeal perfusion circuit. This phantom exhibited reproducible dynamic contrast-enhanced CT of the hepatic arterial and portal venous system over a 4- to 6-month period.

Noise reduction in spectral CT: reducing dose and breaking the trade-off between image noise and energy bin selection

PURPOSE

Our purpose was to reduce image noise in spectral CT by exploiting data redundancies in the energy domain to allow flexible selection of the number, width, and location of the energy bins.

METHODS

Using a variety of spectral CT imaging methods, conventional filtered backprojection (FBP) reconstructions were performed and resulting images were compared to those processed using a Local HighlY constrained backPRojection Reconstruction (HYPR-LR) algorithm. The mean and standard deviation of CT numbers were measured within regions of interest (ROIs), and results were compared between FBP and HYPR-LR. For these comparisons, the following spectral CT imaging methods were used:(i) numerical simulations based on a photon-counting, detector-based CT system, (ii) a photon-counting, detector-based micro CT system using rubidium and potassium chloride solutions, (iii) a commercial CT system equipped with integrating detectors utilizing tube potentials of 80, 100, 120, and 140 kV, and (iv) a clinical dual-energy CT examination. The effects of tube energy and energy bin width were evaluated appropriate to each CT system.

RESULTS

The mean CT number in each ROI was unchanged between FBP and HYPR-LR images for each of the spectral CT imaging scenarios, irrespective of bin width or tube potential. However, image noise, as represented by the standard deviation of CT numbers in each ROI, was reduced by 36%-76%. In all scenarios, image noise after HYPR-LR algorithm was similar to that of composite images, which used all available photons. No difference in spatial resolution was observed between HYPR-LR processing and FBP. Dual energy patient data processed using HYPR-LR demonstrated reduced noise in the individual, low- and high-energy images, as well as in the material-specific basis images.

CONCLUSIONS

Noise reduction can be accomplished for spectral CT by exploiting data redundancies in the energy domain. HYPR-LR is a robust method for reducing image noise in a variety of spectral CT imaging systems without losing spatial resolution or CT number accuracy. This method improves the flexibility to select energy bins in the manner that optimizes material identification and separation without paying the penalty of increased image noise or its corollary, increased patient dose.

Newer CT applications and their alternatives: what is appropriate in children?

Innovations in image acquisition and reconstruction technologies have greatly expanded the range of CT applications available in the routine clinical setting. CT images of sub-millimeter resolution can now be acquired of entire body regions in a few seconds or even sub-second time, allowing depiction of fine anatomical detail uncompromised by motion artifact. With sophisticated visualization software, image data can be processed into multiplanar, volume-rendered, cine and other formats to better display anatomical abnormalities and facilitate newer applications such as CT angiography, enterography, urography, tracheobronchography and cardiac CT. Newer applications including dual-energy material decomposition CT are furthering the transition of CT from a purely morphological to a combined anatomical, functional and metabolic imaging technique. These newer applications have largely been pioneered in adult populations, and heightened concern of the risk of carcinogenesis from ionizing radiation tempers dissemination of their use in children. Similar information can often be gleaned from alternative imaging modalities without ionizing radiation exposure, such as MRI and US, and what is most appropriate in children will depend on relative diagnostic efficacy, cost, availability and local expertise.

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.

Renal perfusion: noninvasive measurement with multidetector CT versus fluorescent microspheres in a pig model

PURPOSE

To validate the measurement of renal perfusion with multidetector computed tomography (CT) with a low-rate injection of contrast medium (ie, 3 mL/sec) through a catheter placed peripherally with gamma variate extended modeling in a pig model, compared with a reference method of fluorescent microspheres.

MATERIALS AND METHODS

This study was approved by the Institutional Animal Care and Use Committee. Renal perfusion was measured in 10 anesthetized pigs simultaneously with multidetector CT and with fluorescent microspheres, which are the reference standard for measuring regional renal perfusion. In each pig, measurements were obtained under three conditions. These were dopamine infusion, dopamine infusion with vascular expansion, and angiotensin II infusion. Aortic and cortical time-attenuation curves were modeled to measure renal perfusion with the gamma variate model. The renal perfusion measurements with the multidetector CT and that with microspheres were compared with least squares regression analysis and Bland-Altman plots.

RESULTS

Perfusion as measured with multidetector CT and that as measured with microspheres were strongly correlated (ρ = 0.93, P < .0001). Multidetector CT renal perfusion with dopamine infusion (3.13 mL/min/g ± 0.53) was not changed after volume expansion (3.37 mL/min/g ± 0.75, P = .35) but was significantly decreased after angiotensin II injection (2.01 mL/min/g ± 0.57, P = .0001).

CONCLUSION

Multidetector CT provides reliable measurements of single-kidney perfusion with peripheral low-rate contrast medium injection.

Noise reduction and image quality improvement of low dose and ultra low dose brain perfusion CT by HYPR-LR processing

PURPOSE

To evaluate image quality and signal characteristics of brain perfusion CT (BPCT) obtained by low-dose (LD) and ultra-low-dose (ULD) protocols with and without post-processing by highly constrained back-projection (HYPR)-local reconstruction (LR) technique.

METHODS AND MATERIALS

Simultaneous BPCTs were acquired in 8 patients on a dual-source-CT by applying LD (80 kV, 200 mAs, 14×1.2 mm) on tube A and ULD (80 kV, 30 mAs, 14×1.2 mm) on tube B. Image data from both tubes was reconstructed with identical parameters and post-processed using the HYPR-LR. Correlation coefficients between mean and maximum (MAX) attenuation values within corresponding ROIs, area under attenuation curve (AUC), and signal to noise ratio (SNR) of brain parenchyma were assessed. Subjective image quality was assessed on a 5-point scale by two blinded observers (1: excellent, 5: non-diagnostic).

RESULTS

Radiation dose of ULD was more than six times lower compared to LD. SNR was improved by HYPR: ULD vs. ULD+HYPR: 1.9±0.3 vs. 8.4±1.7, LD vs. LD+HYPR: 5.0±0.7 vs. 13.4±2.4 (both p<0.0001). There was a good correlation between the original datasets and the HYPR-LR post-processed datasets: r = 0.848 for ULD and ULD+HYPR and r = 0.933 for LD and LD+HYPR (p<0.0001 for both). The mean values of the HYPR-LR post-processed ULD dataset correlated better with the standard LD dataset (r = 0.672) than unprocessed ULD (r = 0.542), but both correlations were significant (p<0.0001). There was no significant difference in AUC or MAX. Image quality was rated excellent (1.3) in LD+HYPR and non-diagnostic (5.0) in ULD. LD and ULD+HYPR images had moderate image quality (3.3 and 2.7).

CONCLUSION

SNR and image quality of ULD-BPCT can be improved to a level similar to LD-BPCT when using HYPR-LR without distorting attenuation measurements. This can be used to substantially reduce radiation dose. Alternatively, LD images can be improved by HYPR-LR to higher diagnostic quality.

Radiation dose reduction in computed tomography: techniques and future perspective

Despite universal consensus that computed tomography (CT) overwhelmingly benefits patients when used for appropriate indications, concerns have been raised regarding the potential risk of cancer induction from CT due to the exponentially increased use of CT in medicine. Keeping radiation dose as low as reasonably achievable, consistent with the diagnostic task, remains the most important strategy for decreasing this potential risk. This article summarizes the general technical strategies that are commonly used for radiation dose management in CT. Dose-management strategies for pediatric CT, cardiac CT, dual-energy CT, CT perfusion and interventional CT are specifically discussed, and future perspectives on CT dose reduction are presented.