Dual-energy material decomposition for cone-beam computed tomography in image-guided radiotherapy

Technical note: TIGRE-DE for the creation of virtual monoenergetic images from dual-energy cone-beam CT

BACKGROUND

Dual-energy (DE)-CBCT represents a promising imaging modality that can produce virtual monoenergetic (VM) CBCT images. VM images, which provide enhanced contrast and reduced imaging artifacts, can be used to assist in soft-tissue visualization during image-guided radiotherapy.

PURPOSE

This work reports the development of TIGRE-DE, a module in the open-source TIGRE toolkit for the performance of DE-CBCT and the production of VM CBCT images. This module is created to make DE-CBCT tools accessible in a wider range of clinical and research settings.

METHODS

We developed an add-on (TIGRE-DE) to the TIGRE toolkit that performs DE material decomposition. To verify its performance, sequential CBCT scans at 80 and 140 kV of a Catphan 604 phantom were decomposed into equivalent thicknesses of aluminum (Al) and polymethyl-methylacrylate (PMMA) basis materials. These basis material projections were used to synthesize VM projections for a range of x-ray energies, which were then reconstructed using the Feldkamp-Davis-Kress (FDK) algorithm. Image quality was assessed by computing Hounsfield units (HU) and contrast-to-noise ratios (CNR) for the material inserts of the phantom and comparing with the constituent 80 and 140 kV images.

RESULTS

All VM images generated using TIGRE-DE showed good general agreement with the theoretical HU values of the material inserts of the phantom. Apart from the highest-density inserts imaged at the extremes of the energy range, the measured HU values agree with theoretical HUs within the clinical tolerance of ±50 HU. CNR measurements for the various inserts showed that, of the energies selected, 60 keV provided the highest CNR values. Moreover, 60 keV VM images showed average CNR enhancements of 63% and 66% compared to the 80 and 140 kV full-fan protocols.

CONCLUSIONS

TIGRE-DE successfully implements DE-CBCT material decomposition and VM image creation in an accessible, open-source platform.

Feasibility of dual-energy CBCT material decomposition in the human torso with 2D anti-scatter grids and grid-based scatter sampling

BACKGROUND

Dual-energy (DE) imaging techniques in cone-beam computed tomography (CBCT) have potential clinical applications, including material quantification and improved tissue visualization. However, the performance of DE CBCT is limited by the effects of scattered radiation, which restricts its use to small object imaging.

PURPOSE

This study investigates the feasibility of DE CBCT material decomposition by reducing scatter with a 2D anti-scatter grid and a measurement-based scatter correction method. Specifically, the investigation focuses on iodine quantification accuracy and virtual monoenergetic (VME) imaging in phantoms that mimic head, thorax, abdomen, and pelvis anatomies.

METHODS

A 2D anti-scatter grid prototype was utilized with a residual scatter correction method in a linac-mounted CBCT system to investigate the effects of robust scatter suppression in DE CBCT. Scans were acquired at 90 and 140 kVp using phantoms that mimic head, thorax, and abdomen/pelvis anatomies. Iodine vials with varying concentrations were placed in each phantom, and CBCT images were decomposed into iodine and water basis material images. The effect of a 2D anti-scatter grid with and without residual scatter correction on iodine concentration quantification and contrast visualization in VME images was evaluated. To benchmark iodine concentration quantification accuracy, a similar set of experiments and DE processing were also performed with a conventional multidetector CT scanner.

RESULTS

In CBCT images, a 2D grid with or without scatter correction can differentiate iodine and water after DE processing in human torso-sized phantom images. However, iodine quantification errors were up to 10 mg/mL in pelvis phantoms when only the 2D grid was used. Adding scatter correction to 2D-grid CBCT reduced iodine quantification errors below 1.5 mg/mL in pelvis phantoms, comparable to iodine quantification errors in multidetector CT. While a noticeable contrast-to-noise ratio improvement was not observed in VME CBCT images, contrast visualization was substantially better in 40 keV VME images in visual comparisons with 90 and 140 kVp CBCT images across all phantom sizes investigated.

CONCLUSIONS

This study indicates that accurate DE decomposition is potentially feasible in DE CBCT of the human torso if robust scatter suppression is achieved with 2D anti-scatter grids and residual scatter correction. This approach can potentially enable better contrast visualization and tissue and contrast agent quantification in various CBCT applications.

Split-filter dual energy computed tomography radiotherapy: From calibration to image guidance

Background and purpose

Dual-energy computed tomography (DECT) is an emerging technology in radiotherapy (RT). Here, we investigate split-filter DECT throughout the RT treatment chain as compared to single-energy CT (SECT).

Materials and methods

DECT scans were acquired with a tin-gold split-filter at 140 kV resulting in a low- and high-energy CT reconstruction (recon). Ten cancer patients (four head-and-neck (HN)​, three rectum​, two anal/pelvis and one abdomen) were DECT scanned without and with iodine administered. A cylindrical and an anthropomorphic HN phantom were scanned with DECT and 120 kV SECT. The DECT images generated were: 120 kV SECT-equivalent (CTmix), virtual monoenergetic images (VMIs), iodine map, virtual non-contrast (VNC), effective atomic number (Zeff), and relative electron density (ρe,w). The clinical utility of these recons was investigated for calibration, delineation, dose calculation and image-guided RT (IGRT).

Results

A calibration curve for 75 keV VMI had a root-mean-square-error (RMSE) of 34 HU in closest agreement with the RSME of SECT calibration. This correlated with a phantom-based dosimetric agreement to SECT of γ1%1mm > 98%. A 40 keV VMI recon was most promising to improve tumor delineation accuracy with an average evaluation score of 1.6 corresponding to "partial improvement". The dosimetric impact of iodine was in general < 2%. For this setup, VNC vs. non-contrast CTmix based dose calculations are considered equivalent. SECT- and DECT-based IGRT was in agreement within the setup uncertainty.

Conclusions

DECT-based RT could be a feasible alternative to SECT providing additional recons to support the different steps of the RT workflow.

A generalizable new figure of merit for dose optimization in dual energy cone beam CT scanning protocols

Objective. This study proposes and evaluates a new figure of merit (FOMn) for dose optimization of Dual-energy cone-beam CT (DE-CBCT) scanning protocols based on size-dependent modeling of radiation dose and multi-scale image quality.Approach. FOMn was defined using Z-score normalization and was proportional to the dose efficiency providing better multi-scale image quality, including comprehensive contrast-to-noise ratio (CCNR) and electron density (CED) for CatPhan604 inserts of various materials. Acrylic annuluses were combined with CatPhan604 to create four phantom sizes (diameters of the long axis are 200 mm, 270 mm, 350 mm, and 380 mm, respectively). DE-CBCT was decomposed using image-domain iterative methods based on Varian kV-CBCT images acquired using 25 protocols (100 kVp and 140 kVp combined with 5 tube currents).Main results. The accuracy of CED was approximately 1% for all protocols, but degraded monotonically with the increased phantom sizes. Combinations of lower voltage + higher current and higher voltage + lower current were optimal protocols balancing CCNR and dose. The most dose-efficient protocols for CED and CCNR were inconsistent, underlining the necessity of including multi-scale image quality in the evaluation and optimization of DE-CBCT. Pediatric and adult anthropomorphic phantom tests confirmed dose-efficiency of FOMn-recommended protocols.Significance. FOMn is a comprehensive metric that collectively evaluates radiation dose and multi-scale image quality for DE-CBCT. The models and data can also serve as lookup tables, suggesting personalized dose-efficient protocols for specific clinical imaging purposes.

TIME-Net: Transformer-Integrated Multi-Encoder Network for limited-angle artifact removal in dual-energy CBCT

Dual-energy cone-beam computed tomography (DE-CBCT) is a promising imaging technique with foreseeable clinical applications. DE-CBCT images acquired with two different spectra can provide material-specific information. Meanwhile, the anatomical consistency and energy-domain correlation result in significant information redundancy, which could be exploited to improve image quality. In this context, this paper develops the Transformer-Integrated Multi-Encoder Network (TIME-Net) for DE-CBCT to remove the limited-angle artifacts. TIME-Net comprises three encoders (image encoder, prior encoder, and transformer encoder), two decoders (low- and high-energy decoders), and one feature fusion module. Three encoders extract various features for image restoration. The feature fusion module compresses these features into more compact shared features and feeds them to the decoders. Two decoders perform differential learning for DE-CBCT images. By design, TIME-Net could obtain high-quality DE-CBCT images using two complementary quarter-scans, holding great potential to reduce radiation dose and shorten the acquisition time. Qualitative and quantitative analyses based on simulated data and real rat data have demonstrated the promising performance of TIME-Net in artifact removal, subtle structure restoration, and reconstruction accuracy preservation. Two clinical applications, virtual non-contrast (VNC) imaging and iodine quantification, have proved the potential utility of the DE-CBCT images provided by TIME-Net.

Planning CT-guided robust and fast cone-beam CT scatter correction using a local filtration technique

PURPOSE

Cone-beam CT (CBCT) has been widely utilized in image-guided radiotherapy. Planning CT (pCT)-aided CBCT scatter correction could further enhance image quality and extend CBCT application to dose calculation and adaptive planning. Nevertheless, existing pCT-based approaches demand accurate registration between pCT and CBCT, leading to limited imaging performance and increased computational cost when large anatomical discrepancies exist. In this work, we proposed a robust and fast CBCT scatter correction method using local filtration technique and rigid registration between pCT and CBCT (LF-RR).

METHODS

First of all, the pCT was rigidly registered with CBCT, then forward projection was performed on registered pCT to create scatter-free projections. The raw scatter signals were obtained via subtracting the scatter-free projections from the measured CBCT projections. Based on frequency and intensity threshold criteria, reliable scatter signals were selected from the raw scatter signals, and further filtered for global scatter estimation via local filtration technique. Finally, corrected CBCT was reconstructed with the projections generated by subtracting the scatter estimation from the raw CBCT projections using FDK algorithm. The LF-RR method was evaluated via comparison with another pCT-based scatter correction method based on Median and Gaussian filters (MG method).

RESULTS

Proposed method was first validated on an anthropomorphic pelvis phantom, and showed satisfied performance on scatter removal even when anatomical mismatches were intentionally created on pCT. The quantitative analysis was further performed on four clinical CBCT images. Compared with the uncorrected CBCT, CBCT corrected by MG with rigid registration (MG-RR), MG with deformable registration (MG-DR), and LF-RR reduced the CT number error from 79 ± 35 to 25 ± 18 , 17 ± 13 and 7 ± 3 HU for adipose and from 115 ± 61 to 36 ± 22 , 30 ± 24 , 7 ± 3 HU for muscle, respectively. After correction, the spatial non-uniformity (SNU) of CBCT corrected with MG-RR, MG-DR and LF-RR was 51 ± 13 , 60 ± 21 , and 21 ± 9 HU for adipose, and 50 ± 22 , 57 ± 41 , and 25 ± 6 HU for muscle, respectively. Meanwhile, the contrast-to-noise ratio (CNR) between muscle and adipose was increased by a factor of 2.70, 2.89 and 2.56, respectively. Using the LF-RR method, the scatter correction of 656 projections can be finished within 10 s and the corrected volumetric images (200 slices) can be obtained within 2 min.

CONCLUSION

We developed a fast and robust pCT-based CBCT scatter correction method which exploits the local-filtration technique to promote the accuracy of scatter estimation and is resistant to pCT-to-CBCT registration uncertainties. Both phantom and patient studies showed the superiority of the proposed correction in imaging accuracy and computational efficiency, indicating promisingfuture clinical application.

Performance characterization of a prototype dual-layer cone-beam computed tomography system

PURPOSE

Conventional cone-beam computed tomography CT (CBCT) provides limited discrimination between low-contrast tissues. Furthermore, it is limited to full-spectrum energy integration. A dual-energy CBCT system could be used to separate photon energy spectra with the potential to increase the visibility of clinically relevant features and acquire additional information relevant in a multitude of clinical imaging applications. In this work, the performance of a novel dual-layer dual-energy CBCT (DL-DE-CBCT) C-arm system is characterized for the first time.

METHODS

A prototype dual-layer detector was fitted into a commercial interventional C-arm CBCT system to enable DL-DE-CBCT acquisitions. DL-DE reconstructions were derived from material-decomposed Compton scatter and photoelectric base functions. The modulation transfer function (MTF) of the prototype DL-DE-CBCT was compared to that of a commercial CBCT. Noise and uniformity characteristics were evaluated using a cylindrical water phantom. Effective atomic numbers and electron densities were estimated in clinically relevant tissue substitutes. Iodine quantification was performed (for 0.5-15 mg/ml concentrations) and virtual noncontrast (VNC) images were evaluated. Finally, contrast-to-noise ratios (CNR) and CT number accuracies were estimated.

RESULTS

The prototype and commercial CBCT showed similar spatial resolution, with a mean 10% MTF of 5.98 cycles/cm and 6.28 cycles/cm, respectively, using a commercial standard reconstruction. The lowest noise was seen in the 80 keV virtual monoenergetic images (VMI) (7.40 HU) and the most uniform images were seen at VMI 60 keV (4.74 HU) or VMI 80 keV (1.98 HU), depending on the uniformity measure used. For all the tissue substitutes measured, the mean accuracy in effective atomic number was 98.2% (SD 1.2%) and the mean accuracy in electron density was 100.3% (SD 0.9%). Iodine quantification images showed a mean difference of -0.1 (SD 0.5) mg/ml compared to the true iodine concentration for all blood and iodine-containing objects. For VNC images, all blood substitutes containing iodine averaged a CT number of 43.2 HU, whereas a blood-only substitute measured 44.8 HU. All water-containing iodine substitutes measured a mean CT number of 2.6 in the VNC images. A noise-suppressed dataset showed a CNR peak at VMI 40 keV and low at VMI 120 keV. In the same dataset without noise suppression applied, a peak in CNR was obtained at VMI 70 keV and a low at VMI 120 keV. The estimated CT numbers of various clinically relevant objects were generally very close to the calculated CT number.

CONCLUSIONS

The performance of a prototype dual-layer dual-energy C-arm CBCT system was characterized. Spatial resolution and noise were comparable with a commercially available C-arm CBCT system, while offering dual-energy capability. Iodine quantifications, effective atomic numbers, and electron densities were in good agreement with expected values, indicating that the system can be used to reliably evaluate the material composition of clinically relevant tissues. The VNC and monoenergetic images indicate a consistent ability to separate clinically relevant tissues. The results presented indicate that the system could find utility in diagnostic, interventional, and radiotherapy planning settings.