The piecewise-linear dynamic attenuator reduces the impact of count rate loss with photon-counting detectors

Photon count rates estimated from 1980 clinical CT scans

BACKGROUND

All photon counting detectors have a characteristic count rate over which their performance degrades. Degradation in the clinical setting takes the form of increased noise, reduced material quantification accuracy, and image artifacts. Count rate is a function of patient attenuation, beam filtration, scanner geometry, and X-ray technique.

PURPOSE

To guide protocol and technology development in the photon counting space, knowledge of clinical count rates spanning the complete range of clinical indications and patient sizes is needed. In this paper, we use clinical data to characterize the range of computed tomography (CT) count rates.

METHODS

We retrospectively gathered 1980 patient exams spanning the entire body (head/neck/chest/abdomen/extremity) and sampled 36 951 axial image slices. We assigned the tissue labels air/lung/fat/soft tissue/bone to each voxel for each slice using CT number thresholds. We then modeled four different bowtie filters, 70/80/100/120/140 kV spectra, and a range of mA values. We forward-projected each slice to obtain detector-incident count rates, using the geometry of a GE Revolution Apex scanner. Our analysis divided the detector into thirds: the central one-third, one-third of the detector split into two equal regions adjacent to the central third, and the final one-third divided equally between the outer detector edges. We report the 99th percentile of counts to mimic the upper limits of count rates making passing through a patient as a function of patient water equivalent diameter. We also report the percentage of patient scans, by body region, over different count rate thresholds for all combinations of bowtie and beam energy.

RESULTS

For routine exam types, we recorded count rates of approximately 3.5 × 10⁸  counts/mm² /s in the torso, extremities, and brain. For neck scans, we observed count rates near 6 × 10⁸  counts/mm² /s. Our simulations of 1000 mA, appropriately mimicking the mA needs for fast pediatric, fast thoracic, and cardiac scanning, resulted in count rates of over 10 × 10⁸  counts/mm² /s for the torso, extremities, and brain. At 1000 mA, for the neck region, we observed count rates close to 2 × 10⁹  counts/mm² /s. Importantly, we saw only a small change in maximum count rate needs over patient size, which we attribute to patient mis-positioning with respect to the bowtie filters. As expected, combinations of kV and bowtie filter with higher beam energies and wider/less attenuating bowtie fluence profiles lead to higher count rates relative to lower energies. The 99th-50th percentile count rate changed the most for the torso region, with a maximum variation of 3.9 × 10⁸ to 1.2 × 10⁷  counts/mm² /s. The head/neck/extremity regions had less than a 50% change in count rate from the 99th to 50th percentiles.

CONCLUSIONS

Our results are the first to use a large patient cohort spanning all body regions to characterize count rates in CT. Our results should be useful in helping researchers understand count rates as a function of body region and mA for various combinations of bowtie filter designs and beam energies. Our results indicate clinical rates >1 × 10⁹  counts/mm² /s, but they do not predict the image quality impact of using a detector with lower characteristic count rates.

Reconstruction of three-dimensional tomographic patient models for radiation dose modulation in CT from two scout views using deep learning

BACKGROUND

A tomographic patient model is essential for radiation dose modulation in x-ray computed tomography (CT). Currently, two-view scout images (also known as topograms) are used to estimate patient models with relatively uniform attenuation coefficients. These patient models do not account for the detailed anatomical variations of human subjects, and thus, may limit the accuracy of intraview or organ-specific dose modulations in emerging CT technologies.

PURPOSE

The purpose of this work was to show that 3D tomographic patient models can be generated from two-view scout images using deep learning strategies, and the reconstructed 3D patient models indeed enable accurate prescriptions of fluence-field modulated or organ-specific dose delivery in the subsequent CT scans.

METHODS

CT images and the corresponding two-view scout images were retrospectively collected from 4214 individual CT exams. The collected data were curated for the training of a deep neural network architecture termed ScoutCT-NET to generate 3D tomographic attenuation models from two-view scout images. The trained network was validated using a cohort of 55 136 images from 212 individual patients. To evaluate the accuracy of the reconstructed 3D patient models, radiation delivery plans were generated using ScoutCT-NET 3D patient models and compared with plans prescribed based on true CT images (gold standard) for both fluence-field-modulated CT and organ-specific CT. Radiation dose distributions were estimated using Monte Carlo simulations and were quantitatively evaluated using the Gamma analysis method. Modulated dose profiles were compared against state-of-the-art tube current modulation schemes. Impacts of ScoutCT-NET patient model-based dose modulation schemes on universal-purpose CT acquisitions and organ-specific acquisitions were also compared in terms of overall image appearance, noise magnitude, and noise uniformity.

RESULTS

The results demonstrate that (1) The end-to-end trained ScoutCT-NET can be used to generate 3D patient attenuation models and demonstrate empirical generalizability. (2) The 3D patient models can be used to accurately estimate the spatial distribution of radiation dose delivered by standard helical CTs prior to the actual CT acquisition; compared to the gold-standard dose distribution, 95.0% of the voxels in the ScoutCT-NET based dose maps have acceptable gamma values for 5 mm distance-to-agreement and 10% dose difference. (3) The 3D patient models also enabled accurate prescription of fluence-field modulated CT to generate a more uniform noise distribution across the patient body compared to tube current-modulated CT. (4) ScoutCT-NET 3D patient models enabled accurate prescription of organ-specific CT to boost image quality for a given body region-of-interest under a given radiation dose constraint.

CONCLUSION

3D tomographic attenuation models generated by ScoutCT-NET from two-view scout images can be used to prescribe fluence-field-modulated or organ-specific CT scans with high accuracy for the overall objective of radiation dose reduction or image quality improvement for a given imaging task.

On the potential of ROI imaging in x-ray CT - A comparison of novel dynamic beam attenuators with current technology

PURPOSE

In this work, we explore the potential of region-of-interest (ROI) imaging in x-ray computed tomography (CT). Using two dynamic beam attenuator (DBA) concepts for fluence field modulation (FFM) previously developed, we investigate and evaluate the potential dose savings in comparison with current FFM technology.

METHODS

ROI imaging is a special application of FFM where the bulk of x-ray radiation is propagated toward a certain anatomical target (ROI), specified by the imaging task, while the surrounding tissue is spared from radiation. We introduce a criterion suitable to quantitatively describe the balance between image quality inside an ROI and total radiation dose with respect to a given ROI imaging task. It accounts for the mean image variance at the ROI and the effective patient dose calculated from Monte Carlo simulations. The criterion is further used to compile task-specific DBA trajectories determining the primary x-ray fluence, and eventually used for comparing different FFM techniques, namely the sheet-based dynamic beam attenuator (sbDBA), the z-aligned sbDBA (z-sbDBA), and an adjustable static operation mode of the z-sbDBA. Furthermore, two static bowtie filters and the influence of tube current modulation (TCM) are included in the comparison.

RESULTS

Our findings demonstrate by simulations that the presented trajectory optimization method determines reasonable DBA trajectories. The influence of TCM is strongly depending on the imaging task. The narrow bowtie filter allows for dose reductions of about 10% compared to the regular bowtie filter in the considered ROI imaging tasks. The DBAs are shown to realize substantially larger dose reductions. In our cardiac imaging scenario, the DBAs can reduce the effective dose by about 30% (z-sbDBA) or 60% (sbDBA). We can further verify that the noise characteristics are not adversely affected by the DBAs.

CONCLUSION

Our research demonstrates that ROI imaging using the presented DBA concepts is a promising technique toward a more patient- and task-specific CT imaging requiring lower radiation dose. Both the sbDBA and the z-sbDBA are potential technical solutions for realizing ROI imaging in x-ray CT.

Anomalous edge response of cadmium telluride-based photon counting detectors jointly caused by high-flux radiation and inter-pixel communication

This work reports an edge enhancing effect experimentally observed in cadmium telluride (CdTe)-based photon counting detector (PCD) systems operated under the charge summing (CS) mode and irradiated by high-flux x-rays. Experimental measurements of the edge spread functions (ESFs) of a PCD system (100μm pixel size, 88 ns deadtime) were performed at different input flux levels from 4.5 × 105count per second (cps) mm-2to 1.5 × 109cps mm-2for the single pixel mode (SP) and the CS mode. A theoretical model that incorporates the impacts of inter-pixel communications and the arbitration process involved in the CS mode was developed to help explain the physical origin of the observed edge enhancing effect. Compared with the monotonically increasing ESF of the SP mode, the ESF of the CS mode measured at high-flux levels shows a peak at an intermediate location (50μm from the edge). The peak became more pronounced with increasing flux levels. The theoretically calculated ESFs agreed well with experimental results with relative errors less than 5% at all flux levels and tested. These results indicate that the anomalous edge enhancing effect is jointly caused by the pileup effect and the CS circuit that introduces negative correlations between adjacent pixels. When the input flux is high enough to deliver photons to multiple adjacent pixels within the same deadtime period, the CS mode may treat the coincident x-rays as shared charges and thus introduce count losses in addition to the well-known pileup count loss. When a high contrast object partially blocks certain pixels from x-rays, the adjacent unblocked pixels have an increased probability of registering counts as a result of the negative correlation. This leads to a peak on the ESF at a pixel-to-edge distance half of the pixel pitch.

Generation and analysis of clinically relevant breast imaging x-ray spectra

PURPOSE

The purpose of this work was to develop and make available x-ray spectra for some of the most widely used digital mammography (DM), breast tomosynthesis (BT), and breast CT (bCT) systems in North America.

METHODS

The Monte Carlo code MCNP6 was used to simulate minimally filtered (only beryllium) x-ray spectra at 8 tube potentials from 20 to 49 kV for DM/BT, and 9 tube potentials from 35 to 70 kV for bCT. Vendor-specific anode compositions, effective anode angles, focal spot sizes, source-to-detector distances, and beryllium filtration were simulated. For each 0.5 keV energy bin in all simulated spectra, the fluence was interpolated using cubic splines across the range of simulated tube potentials to produce spectra in 1 kV increments from 20 to 49 kV for DM/BT and from 35 to 70 kV for bCT. The HVL of simulated spectra with conventional filtration (at 35 kV for DM/BT and 49 kV for bCT) was used to assess spectral differences resulting from variations in: (a) focal spot size (0.1 and 0.3 mm IEC), (b) solid angle at the detector (i.e., small and large FOV size), and (c) geometrical specifications for vendors that employ the same anode composition.

RESULTS

Averaged across all DM/BT vendors, variations in focal spot and FOV size resulted in HVL differences of 2.2% and 0.9%, respectively. Comparing anode compositions separately, the HVL differences for Mo (GE, Siemens) and W (Hologic, Philips, and Siemens) spectra were 0.3% and 0.6%, respectively. Both the commercial Koning and prototype "Doheny" (UC Davis) bCT systems utilize W anodes with a 0.3 mm focal spot. Averaged across both bCT systems, variations in FOV size resulted in a 2.2% difference in HVL. In addition, the Koning spectrum was slightly harder than Doheny with a 4.2% difference in HVL. Therefore to reduce redundancy, a generic DM/BT system and a generic bCT system were used to generate the new spectra reported herein. The spectral models for application to DM/BT were dubbed the Molybdenum, Rhodium, and Tungsten Anode Spectral Models using Interpolating Cubic Splines (MASMICS_M-T , RASMICS_M-T , and TASMICS_M-T ; subscript "M-T" indicating mammography and tomosynthesis). When compared against reference models (MASMIP_M , RASMIP_M , and TASMIP_M ; subscript "M" indicating mammography), the new spectral models were in close agreement with mean differences of 1.3%, -1.3%, and -3.3%, respectively, across tube potential comparisons of 20, 30, and 40 kV with conventional filtration. TASMICS_bCT -generated bCT spectra were also in close agreement with the reference TASMIP model with a mean difference of -0.8%, across tube potential comparisons of 35, 49, and 70 kV with 1.5 mm Al filtration.

CONCLUSIONS

The Mo, Rh, and W anode spectra for application in DM and BT (MASMICS_M-T , RASMICS_M-T , and TASMICS_M-T ) and the W anode spectra for bCT (TASMICS_bCT ) as described in this study should be useful for individuals interested in modeling the performance of modern breast x-ray imaging systems including dual-energy mammography which extends to 49 kV. These new spectra are tabulated in spreadsheet form and are made available to any interested party.

An Investigation of Low-Dose 3D Scout Scans for Computed Tomography

PURPOSE

Commonly 2D scouts or topograms are used prior to CT scan acquisition. However, low-dose 3D scouts could potentially provide additional information for more effective patient positioning and selection of acquisition protocols. We propose using model-based iterative reconstruction to reconstruct low exposure tomographic data to maintain image quality in both low-dose 3D scouts and reprojected topograms based on those 3D scouts.

METHODS

We performed tomographic acquisitions on a CBCT test-bench using a range of exposure settings from 16.6 to 231.9 total mAs. Both an anthropomorphic phantom and a 32 cm CTDI phantom were scanned. The penalized-likelihood reconstructions were made using Matlab and CUDA libraries and reconstruction parameters were tuned to determine the best regularization strength and delta parameter. RMS error between reconstructions and the highest exposure reconstruction were computed, and CTDI_W values were reported for each exposure setting. RMS error for reprojected topograms were also computed.

RESULTS

We find that we are able to produce low-dose (0.417 mGy) 3D scouts that show high-contrast and large anatomical features while maintaining the ability to produce traditional topograms.

CONCLUSIONS

We demonstrated that iterative reconstruction can mitigate noise in very low exposure CT acquisitions to enable 3D CT scout. Such additional 3D information may lead to improved protocols for patient positioning and acquisition refinements as well as a number of advanced dose reduction strategies that require localization of anatomical features and quantities that are not provided by simple 2D topograms.

Technical Note: spektr 3.0-A computational tool for x-ray spectrum modeling and analysis

PURPOSE

A computational toolkit (spektr 3.0) has been developed to calculate x-ray spectra based on the tungsten anode spectral model using interpolating cubic splines (TASMICS) algorithm, updating previous work based on the tungsten anode spectral model using interpolating polynomials (TASMIP) spectral model. The toolkit includes a matlab (The Mathworks, Natick, MA) function library and improved user interface (UI) along with an optimization algorithm to match calculated beam quality with measurements.

METHODS

The spektr code generates x-ray spectra (photons/mm(2)/mAs at 100 cm from the source) using TASMICS as default (with TASMIP as an option) in 1 keV energy bins over beam energies 20-150 kV, extensible to 640 kV using the TASMICS spectra. An optimization tool was implemented to compute the added filtration (Al and W) that provides a best match between calculated and measured x-ray tube output (mGy/mAs or mR/mAs) for individual x-ray tubes that may differ from that assumed in TASMICS or TASMIP and to account for factors such as anode angle.

RESULTS

The median percent difference in photon counts for a TASMICS and TASMIP spectrum was 4.15% for tube potentials in the range 30-140 kV with the largest percentage difference arising in the low and high energy bins due to measurement errors in the empirically based TASMIP model and inaccurate polynomial fitting. The optimization tool reported a close agreement between measured and calculated spectra with a Pearson coefficient of 0.98.

CONCLUSIONS

The computational toolkit, spektr, has been updated to version 3.0, validated against measurements and existing models, and made available as open source code. Video tutorials for the spektr function library, UI, and optimization tool are available.

A prototype piecewise-linear dynamic attenuator

The piecewise-linear dynamic attenuator has been proposed as a mechanism in CT scanning for personalizing the x-ray illumination on a patient- and application-specific basis. Previous simulations have shown benefits in image quality, scatter, and dose objectives. We report on the first prototype implementation. This prototype is reduced in scale and speed and is integrated into a tabletop CT system with a smaller field of view (25 cm) and longer scan time (42 s) compared to a clinical system. Stainless steel wedges were machined and affixed to linear actuators, which were in turn held secure by a frame built using rapid prototyping technologies. The actuators were computer-controlled, with characteristic noise of about 100 microns. Simulations suggest that in a clinical setting, the impact of actuator noise could lead to artifacts of only 1 HU. Ring artifacts were minimized by careful design of the wedges. A water beam hardening correction was applied and the scan was collimated to reduce scatter. We scanned a 16 cm water cylinder phantom as well as an anthropomorphic pediatric phantom. The artifacts present in reconstructed images are comparable to artifacts normally seen with this tabletop system. Compared to a flat-field reference scan, increased detectability at reduced dose is shown and streaking is reduced. Artifacts are modest in our images and further refinement is possible. Issues of mechanical speed and stability in the challenging clinical CT environment will be addressed in a future design.

A dynamic attenuator improves spectral imaging with energy-discriminating, photon counting detectors

Energy-discriminating, photon counting (EDPC) detectors have high potential in spectral imaging applications but exhibit degraded performance when the incident count rate approaches or exceeds the characteristic count rate of the detector. In order to reduce the requirements on the detector, we explore the strategy of modulating the X-ray flux field using a recently proposed dynamic, piecewise-linear attenuator. A previous paper studied this modulation for photon counting detectors but did not explore the impact on spectral applications. In this work, we modeled detection with a bipolar triangular pulse shape (Taguchi et al., 2011) and estimated the Cramer-Rao lower bound (CRLB) of the variance of material selective and equivalent monoenergetic images, assuming deterministic errors at high flux could be corrected. We compared different materials for the dynamic attenuator and found that rare earth elements, such as erbium, outperformed previously proposed materials such as iron in spectral imaging. The redistribution of flux reduces the variance or dose, consistent with previous studies on benefits with conventional detectors. Numerical simulations based on DICOM datasets were used to assess the impact of the dynamic attenuator for detectors with several different characteristic count rates. The dynamic attenuator reduced the peak incident count rate by a factor of 4 in the thorax and 44 in the pelvis, and a 10 Mcps/mm (2) EDPC detector with dynamic attenuator provided generally superior image quality to a 100 Mcps/mm (2) detector with reference bowtie filter for the same dose. The improvement is more pronounced in the material images.