The introduction of photon-counting detectors in x-ray computed tomography raises the question of how reconstruction algorithms should be adapted to photon-counting measurement data. The transition from energy-integrating to photon-counting detectors introduces new effects into the data model, such as pure Poisson statistics and increased cross talk between detector pixels, (e.g. due to charge sharing), but it is still not known in detail how these effects can be treated accurately by the reconstruction algorithm. In this work, we propose a new reconstruction method based on penalized-likelihood reconstruction that incorporates these effects. By starting from a simple, easily-solved reconstruction problem and adding correction terms for the additional physical effects, we obtain a series expansion for the solution to the image reconstruction problem. This approach serves the twofold purpose of (1) yielding a new, potentially faster method of incorporating complex detector models in the reconstruction process and (2) providing insight into the impact of the non-ideal physical effects on the reconstructed image. We investigate the potential for reconstructing images from simulated photon-counting energy-resolving CT data with the new algorithm by including correction terms representing pure Poisson statistics and interpixel cross talk; and we investigate the impact of these physical effects on the reconstructed images. Results indicate that using two correction terms gives good agreement with the converged solution, suggesting that the new method is feasible in practice. This new approach to image reconstruction can help in developing improved reconstruction algorithms for photon-counting CT.
Purpose In computed tomography (CT), miscalibrated or imperfect detector elements produce stripe artifacts in the sinogram. The stripe artifacts in Radon space are responsible for concentric ring artifacts in the reconstructed images. In this work, a novel optimization model is proposed to remove the ring artifacts in an iterative reconstruction procedure. Method In the proposed optimization model, a novel ring total variation (RTV) regularization is developed to penalize the ring artifacts in the image domain. Moreover, to correct the sinogram, a new correcting vector is proposed to compensate for malfunctioning of detectors in the projection domain. The optimization problem is solved by using the alternating minimization scheme (AMS). In each iteration, the fidelity term along with the RTV regularization is solved using the alternating direction method of multipliers (ADMM) to find the image, and then the correcting coefficient vector is updated for certain detectors according to the obtained image. Because the sinogram and the image are simultaneously updated, the proposed method basically performs in both image and sinogram domains. Results The proposed method is evaluated using both simulated and physical phantom datasets containing different ring artifact patterns. In the simulated datasets, the Shepp-Logan phantom, a real chest scan image and a noisy low-contrast phantom are considered for the performance evaluation of our method. We compare the quantitative root mean square error (RMSE) and structural similarity (SSIM) results of our algorithm with wavelet-Fourier sinogram filtering method by Munch et al., the ring artifact reduction method by Brun et al., and the TV-based ring correction method by Paleo and Mirone. Our proposed method is also evaluated using a physical phantom dataset where strong ring artifacts are manifest due to the miscalibration of a large number of detectors. Our proposed method outperforms the competing methods in terms of both qualitative and quantitative evaluation results. Conclusion The experimental results in both simulated and physical phantom datasets show that the proposed method achieves the state-of-the-art ring artifact reduction performance in terms of RMSE, SSIM, and subjective visual quality.
Purpose To compare the CT imaging performance of a carboxybetaine zwitterionic-coated tantalum oxide (TaCZ) nanoparticle CT contrast agent with that of a conventional iodinated contrast agent in a swine model meant to simulate overweight and obese patients. Materials and Methods Four swine were evaluated inside three different-sized adipose-equivalent encasements emulating abdominal girths of 102, 119, and 137 cm. Imaging was performed with a 64-detector row CT scanner at six scan delays after intravenous injection of 240 mg element (Ta or I) per kilogram of body weight of TaCZ or iopromide. For each time point, contrast enhancement of the aorta and liver were measured by using regions of interest. Two readers independently recorded the clarity of vasculature using a five-point Likert scale. Findings were compared by using paired t tests and Wilcoxon signed-rank tests. Results Mean peak enhancement was higher for TaCZ than for iopromide in the aorta (270 HU [σ = 24.5] vs 199 HU [σ = 10.2], P < .001) and liver (61.3 HU [σ = 11.7] vs 45.2 HU [σ = 8], P < .001). Vascular clarity was higher for TaCZ than for iopromide in 63% (132 of 208), 82% (170 of 208), and 86% (178 of 208) of the individual vessels at the 102-, 119-, and 137-cm girths, respectively (P < .01). Arterial clarity scores were higher for TaCZ than for iopromide in 62% (208 of 336) of vessels. Venous clarity scores were higher for TaCZ than for iopromide in 89% (128 of 144) of the veins in the venous phase and in 100% (144 of 144) of veins in the delayed phase (P < .01). No vessel showed higher clarity score with iopromide than with TaCZ. Conclusion An experimental tantalum nanoparticle-based contrast agent showed greater contrast enhancement compared with iopromide in swine models meant to simulate overweight and obese patients. © RSNA, 2018.
Fractional Flow Reserve (FFR), the ratio of arterial pressure distal to a coronary lesion to the proximal pressure, is indicative of its hemodynamic significance. This quantity can be determined from invasive measurements made with a catheter, or by using computational methods incorporating models of the the coronary vasculature. One of the inputs needed by a model-based approach for estimating FFR from Computed Tomography Angiography (CTA) images (denoted FFR-CT) is the geometry of the coronary arteries, which requires segmentation of the coronary lumen. Several algorithms have been proposed for coronary lumen segmentation, including the recent application of machine learning techniques. For evaluating these algorithms or for training machine learning algorithms, manual segmentation of the lumen has been considered as ground truth. However, since there is inter-subject variability in manual segmentation, it would be useful to first assess the extent to which this variability affects the predicted FFR values. In the current study, we evaluated the impact of inter-subject variability in manual segmentation on computed FFR, using datasets with three different manual segmentations provided as part of the Rotterdam Coronary Artery Evaluation Framework. FFR was computed using a coronary blood flow model. Our results indicate that variability in manual segmentations on FFR estimates depend on the FFR value. For FFR ≥ 0.97, variability in manual segmentations does not impact FFR estimates, while, for lower FFR values, the variability in manual segmentations leads to significant variability in FFR. The results of this study indicate that researchers should exercise caution when treating manual segmentations as ground truth for estimating FFR from CTA images.
OBJECTIVES:The aim of this study was to assess whether the low- to high-kVp computed tomography (CT) number ratio at dual-energy CT is affected by changes in patient diameter.METHODS:Seven contrast-producing elements were housed sequentially within an abdomen phantom. Fat rings enlarged the phantom diameter from 26 to 44 cm. The phantom was scanned using single-energy CT at tube potentials of 80 and 140 kVp and rapid-kVp-switching dual-energy CT.RESULTS:CT numbers decreased proportionally (∼20% CT number reduction for smallest to largest phantom diameters) for low- and high-energy acquisitions but resulted in consistent dual-energy ratios for each contrast element. For 17 of 21 material pair combinations, the dual-energy ratio ranges of the two elements did not overlap, implying that discrimination should remain possible for these material pairs at all patient sizes.CONCLUSIONS:The dual-energy ratio for different contrast materials is largely unaffected by changes in phantom diameter. This should allow for robust separation of most contrast material combinations irrespective of patient size.
PurposeTo quantitatively evaluate and compare six proposed system architectures for cardiac CT scanning.MethodsStarting from the clinical requirements for cardiac CT, we defined six dedicated cardiac CT architectures. We selected these architectures based on a previous screening study and defined them in sufficient detail to comprehensively analyze their cost and performance. We developed rigorous comparative evaluation methods for the most important aspects of performance and cost, and we applied these evaluation criteria to the defined cardiac CT architectures.ResultsWe found that CT system architectures based on the third‐generation geometry provide nearly linear performance improvement versus the increased cost of additional beam lines (i.e., source–detector pairs), although similar performance improvement could be achieved with advanced motion‐correction algorithms. The third‐generation architectures outperform even the most promising of the proposed architectures that deviate substantially from the traditional CT system architectures.ConclusionThis work confirms the validity of the current trend in commercial CT scanner design. However, we anticipate that over time, CT hardware and software technologies will evolve, the relative importance of the performance criteria will change, the relative costs of components will vary, some of the remaining challenges will be addressed, and perhaps new candidate architectures will be identified; therefore, the conclusion of a comparative analysis like this may change. The evaluation methods that we used can provide a framework for other researchers to analyze their own proposed CT architectures.
Purpose To develop a simple method for producing liquid‐tissue‐surrogate ( LTS ) materials that accurately represent human soft tissues in terms of density and X‐ray attenuation coefficient. Methods and materials We evaluated hypothetical mixtures of water, glycerol, butanol, methanol, sodium chloride, and potassium nitrate; these mixtures were intended to emulate human adipose, blood, brain, kidney, liver, muscle, pancreas, and skin. We compared the hypothetical densities, effective atomic numbers (Z eff ), and calculated discrete‐energy CT attenuation [Hounsfield Units ( HU )] of the proposed materials with those of human tissue elemental composition as specified in International Commission on Radiation Units ( ICRU ) Report 46. We then physically produced the proposed LTS materials for adipose, liver, and pancreas tissue, and we measured the polyenergetic CT attenuation (also expressed as HU ) of these materials within a 32 cm phantom using a 64‐slice clinical CT scanner at 80 kV p, 100 kV p, 120 kV p, and 140 kV p. Results The predicted densities, Z eff , and calculated discrete‐energy CT attenuation of our proposed formulations generally agreed with those of ICRU within < 1% or < 10 HU . For example, the densities of our hypothetical materials agreed precisely with ICRU 's reported values and were 0.95 g/ mL for adipose tissue, 1.04 g/ mL for pancreatic tissue, and 1.06 g/ mL for liver tissue; the discrete‐energy CT attenuation at 60 keV of our hypothetical materials (and ICRU ‐specified compositions) were −107 HU (−113 HU ) for adipose #3, −89 HU (−90 HU ) for adipose #2, 56 HU (55 HU ) for liver tissue, and 31 HU (31 HU ) for pancreatic tissue. The densities of our physically produced materials (compared to ICRU ‐specified compositions) were 0.947 g/ mL (0.0%) for adipose #2, 1.061 g/ mL (+2.0%) for pancreatic tissue, and 1.074 g/ mL (+1.3%) for liver tissue. The empirical polyenergetic CT attenuation measurements of our LTS materials (and the discrete‐energy HU of the ICRU compositions at the mean energy of each spectrum) at 80 kV p were −104 HU (−113 HU ) for adipose #3, −87 HU (−90 HU ) for adipose #2, 59 HU (55 HU ) for liver tissue, and 33 HU (31 HU ) for pancreatic tissue; at 120 kV p, these were −83 HU (−83 HU ) for adipose #3, −68 HU (−63 HU ) for adipose #2, 55 HU (52 HU ) for liver tissue, and 35 HU (33 HU ) for pancreatic tissue. Conclusion Our method for formulating tissue surrogates allowed straightforward production of solutions with CT attenuation that closely matched the target tissues' expected CT attenuation values and trends with kV p. The LTS s' inexpensive and widely available constituent chemicals, combined with their liquid state, should enable rapid production and versatile use among different phantom and experiment types. Further study is warranted, such as the inclusion of contrast agents. These liquid tissue surrogates may potentially accelerate development and testing of advanced CT imaging techniques and technologies.
Coronary Artery Disease (CAD) is the leading cause of death globally [1]. Modern cardiac computed tomography angiography (CCTA) is highly effective at identifying and assessing coronary blockages associated with CAD. The diagnostic value of this anatomical information can be substantially increased in combination with a non-invasive, low dose, correlative, quantitative measure of blood supply to the myocardium. While CT perfusion has shown promise of providing such indications of ischemia, artifacts due to motion, beam hardening, and other factors confound clinical findings and can limit quantitative accuracy. In this paper, we investigate the impact of applying a novel motion correction algorithm to correct for motion in the myocardium. This motion compensation algorithm (originally designed to correct for the motion of the coronary arteries in order to improve CCTA images) has been shown to provide substantial improvements in both overall image quality and diagnostic accuracy of CCTA. We have adapted this technique for application beyond the coronary arteries and present an assessment of its impact on image quality and quantitative accuracy within the context of dual-energy CT perfusion imaging. We conclude that motion correction is a promising technique that can help foster the routine clinical use of dual-energy CT perfusion. When combined, the anatomical information of CCTA and the hemodynamic information from dual-energy CT perfusion should facilitate better clinical decisions about which patients would benefit from treatments such as stent placement, drug therapy, or surgery and help other patients avoid the risks and costs associated with unnecessary, invasive, diagnostic coronary angiography procedures.
BACKGROUND: We are interested in exploring dedicated, high-performance cardiac CT systems optimized to provide the best tradeoff between system cost, image quality, and radiation dose. OBJECTIVE: We sought to identify and evaluate a broad range of CT architectures that could provide an optimal, dedicated cardiac CT solution. METHODS: We identified and evaluated thirty candidate architectures using consistent design choices. We defined specific evaluation metrics related to cost and performance. We then scored the candidates versus the defined metrics. Lastly, we applied a weighting system to combine scores for all metrics into a single overall score for each architecture. CT experts with backgrounds in cardiovascular radiology, x-ray physics, CT hardware and CT algorithms performed the scoring and weighting. RESULTS: We found nearly a twofold difference between the most and the least promising candidate architectures. Architectures employed by contemporary commercial diagnostic CT systems were among the highest-scoring candidates. We identified six architectures that show sufficient promise to merit further in-depth analysis and comparison. CONCLUSION: Our results suggest that contemporary diagnostic CT system architectures outperform most other candidates that we evaluated, but the results for a few alternatives were relatively close. We selected six representative high-scoring candidates for more detailed design and further comparative evaluation.
PURPOSE This paper summarizes the development of a high-power distributed x-ray source, or "multisource," designed for inverse-geometry computed tomography (CT) applications [see B. De Man et al., "Multisource inverse-geometry CT. Part I. System concept and development," Med. Phys. 43, 4607-4616 (2016)]. The paper presents the evolution of the source architecture, component design (anode, emitter, beam optics, control electronics, high voltage insulator), and experimental validation. METHODS Dispenser cathode emitters were chosen as electron sources. A modular design was adopted, with eight electron emitters (two rows of four emitters) per module, wherein tungsten targets were brazed onto copper anode blocks-one anode block per module. A specialized ceramic connector provided high voltage standoff capability and cooling oil flow to the anode. A matrix topology and low-noise electronic controls provided switching of the emitters. RESULTS Four modules (32 x-ray sources in two rows of 16) have been successfully integrated into a single vacuum vessel and operated on an inverse-geometry computed tomography system. Dispenser cathodes provided high beam current (>1000 mA) in pulse mode, and the electrostatic lenses focused the current beam to a small optical focal spot size (0.5 × 1.4 mm). Controlled emitter grid voltage allowed the beam current to be varied for each source, providing the ability to modulate beam current across the fan of the x-ray beam, denoted as a virtual bowtie filter. The custom designed controls achieved x-ray source switching in <1 μs. The cathode-grounded source was operated successfully up to 120 kV. CONCLUSIONS A high-power, distributed x-ray source for inverse-geometry CT applications was successfully designed, fabricated, and operated. Future embodiments may increase the number of spots and utilize fast read out detectors to increase the x-ray flux magnitude further, while still staying within the stationary target inherent thermal limitations.
The introduction of spectral CT imaging in the form of fast clinical dual-energy CT enabled contrast material to be differentiated from other radiodense materials, improved lesion detection in contrast-enhanced scans, and changed the way that existing iodine and barium contrast materials are used in clinical practice. More profoundly, spectral CT can differentiate between individual contrast materials that have different reporter elements such that high-resolution CT imaging of multiple contrast agents can be obtained in a single pass of the CT scanner. These spectral CT capabilities would be even more impactful with the development of contrast materials designed to complement the existing clinical iodine- and barium-based agents. New biocompatible high-atomic number contrast materials with different biodistribution and X-ray attenuation properties than existing agents will expand the diagnostic power of spectral CT imaging without penalties in radiation dose or scan time.
A direct-conversion photon-counting detector has great potential in CT imaging and has been intensively studied in recently years. However, the optimization of the detector configuration based on clinically-relevant metrics is still lacking. In this study, we developed an optimization scheme for CZT-based photon-counting detector design. The spectral response functions characterizing different detector designs were generated using a Monte Carlo-based approach and evaluated with a dual-energy analysis tool. The results demonstrate the degradation in dose efficiency caused by charge sharing when the detector pixel size is below ∼500 μm. The results also show that thinner CZT sensor material has degraded performance due to the lower absorption efficiency. Based on this analysis, ∼500 μm pixel size and 2.5 mm thickness seem to be a good balance between dose efficiency and spatial resolution. The results also show that electronic noise has limited impact on dose efficiency and is not a major risk to the detector development.
Modern diagnostic X-ray computed tomography (CT) usually provides a limited number of standard X-ray tube operating voltages (kVps) that are well calibrated and optimized. In certain scenarios, however, a spectrum using a non-standard kVp may have potential of achieving improved image quality and dose efficiency. For a CT scan using a non-standard kVp, estimating an effective spectrum for calibration is a problem to be solved. An efficient but not exact method of spectrum estimation in practice is to decompose the non-standard spectrum into a linear combination of spectra generated using standard kVps, e.g. spectrum decomposition. This approach works well when the spectrum generated at a non-standard kVp is reasonably close to that generated by a standard kVp. If it is too far away, the residual error (systematic bias) could be significant enough to cause artifact/inaccuracy in CT images. To improve spectrum decomposition, we present a bias-corrected spectrum decomposition approach. A lookup table (LUT) is established in advance to store the spectrum estimation errors using a set of known spectra generated using known kVps from simulation that do not have to be well calibrated. A reasonable and key assumption here is that the spectral estimation errors in the LUT generation would be close to that of actual spectrum estimation errors. Simulations of a dual-energy CT scan at non-standard kVps demonstrate the effectiveness of this approach, leading to more quantitative CT images.
PURPOSE:Traditionally, 2D radiographic preparatory scan images (scout scans) are used to plan diagnostic CT scans. However, a 3D CT volume with a full 3D organ segmentation map could provide superior information for customized scan planning and other purposes. A practical challenge is to design the volumetric scout acquisition and processing steps to provide good image quality (at least good enough to enable 3D organ segmentation) while delivering a radiation dose similar to that of the conventional 2D scout.METHODS:The authors explored various acquisition methods, scan parameters, postprocessing methods, and reconstruction methods through simulation and cadaver data studies to achieve an ultralow dose 3D scout while simultaneously reducing the noise and maintaining the edge strength around the target organ.RESULTS:In a simulation study, the 3D scout with the proposed acquisition, preprocessing, and reconstruction strategy provided a similar level of organ segmentation capability as a traditional 240 mAs diagnostic scan, based on noise and normalized edge strength metrics. At the same time, the proposed approach delivers only 1.25% of the dose of a traditional scan. In a cadaver study, the authors' pictorial-structures based organ localization algorithm successfully located the major abdominal-thoracic organs from the ultralow dose 3D scout obtained with the proposed strategy.CONCLUSIONS:The authors demonstrated that images with a similar degree of segmentation capability (interpretability) as conventional dose CT scans can be achieved with an ultralow dose 3D scout acquisition and suitable postprocessing. Furthermore, the authors applied these techniques to real cadaver CT scans with a CTDI dose level of less than 0.1 mGy and successfully generated a 3D organ localization map.
Metal artifacts have been a problem associated with computed tomography (CT) since its introduction. Recent techniques to mitigate this problem have included utilization of high-energy (keV) virtual monochromatic spectral (VMS) images, produced via dual-energy CT (DECT). A problem with these high-keV images is that contrast enhancement provided by all commercially available contrast media is severely reduced. Contrast agents based on higher atomic number elements can maintain contrast at the higher energy levels where artifacts are reduced. This study evaluated three such candidate elements: bismuth, tantalum, and tungsten, as well as two conventional contrast elements: iodine and barium. A water-based phantom with vials containing these five elements in solution, as well as different artifact-producing metal structures, was scanned with a DECT scanner capable of rapid operating voltage switching. In the VMS datasets, substantial reductions in the contrast were observed for iodine and barium, which suffered from contrast reductions of 97% and 91%, respectively, at 140 versus 40 keV. In comparison under the same conditions, the candidate agents demonstrated contrast enhancement reductions of only 20%, 29%, and 32% for tungsten, tantalum, and bismuth, respectively. At 140 versus 40 keV, metal artifact severity was reduced by 57% to 85% depending on the phantom configuration. (C) 2015 Society of Photo-Optical Instrumentation Engineers (SPIE)
PURPOSETo quantify the computed tomographic (CT) image contrast produced by potentially useful contrast material elements in clinically relevant imaging conditions.MATERIALS AND METHODSEqual mass concentrations (grams of active element per milliliter of solution) of seven radiodense elements, including iodine, barium, gadolinium, tantalum, ytterbium, gold, and bismuth, were formulated as compounds in aqueous solutions. The compounds were chosen such that the active element dominated the x-ray attenuation of the solution. The solutions were imaged within a modified 32-cm CT dose index phantom at 80, 100, 120, and 140 kVp at CT. To simulate larger body sizes, 0.2-, 0.5-, and 1.0-mm-thick copper filters were applied. CT image contrast was measured and corrected for measured concentrations and presence of chlorine in some compounds.RESULTSEach element tested provided higher image contrast than iodine at some tube potential levels. Over the range of tube potentials that are clinically practical for average-sized and larger adults-that is, 100 kVp and higher-barium, gadolinium, ytterbium, and tantalum provided consistently increased image contrast compared with iodine, respectively demonstrating 39%, 56%, 34%, and 24% increases at 100 kVp; 39%, 66%, 53%, and 46% increases at 120 kVp; and 40%, 72%, 65%, and 60% increases at 140 kVp, with no added x-ray filter.CONCLUSIONThe consistently high image contrast produced with 100-140 kVp by tantalum compared with bismuth and iodine at equal mass concentration suggests that tantalum could potentially be favorable for use as a clinical CT contrast agent.
3D volumetric CT images hold the potential to become a rich source of information for 3D organ segmentation and far exceed that made available through 2D radiograph images. Acquiring and generating 3D volumetric images for scan preparation purposes, i.e. 3D scout, while delivering radiation dose equivalent to conventional 2D radiograph is challenging. We explore various acquisition parameters and post-processing methods to reduce dose of a 3D scout while reducing the noise and maintaining the edge strength around the target organ. We demonstrated that similar edge strength and noise to the conventional dose CT scan can be achieved with 3D scout acquisition and post-processing while being dose neutral to a 2D scout acquisition.
Distributed X-ray sources open the way to innovative system concepts in X-ray and computed tomography. They offer promising opportunities in terms of system performance, but they pose unique challenges in terms of source and system technologies. Several academic and industrial teams have proposed a variety of concepts and developed some unique prototypes. We present a broad review of multisource systems. We also discuss X-ray source components and challenges. We close with our perspective on the future prospects of multisource imaging.