Reduction in x-ray scatter and radiation dose for volume-of-interest (VOI) cone-beam breast CT--a phantom study

Novel Technology Allowing Cone Beam Computed Tomography in 6 Seconds: A Patient Study of Comparative Image Quality

PURPOSE

The goal of this study was to evaluate the image quality provided by a novel cone beam computed tomography (CBCT) platform (HyperSight, Varian Medical Systems), a platform with enhanced reconstruction algorithms as well as rapid acquisition times. Image quality was compared with both status quo CBCT for image guidance, and to fan beam CT (FBCT) acquired on a CT simulator (CTsim).

METHODS AND MATERIALS

In a clinical study, 30 individuals were recruited for whom either deep inspiration (DIBH) or deep exhalation breath hold (DEBH) was used during imaging and radiation treatment of tumors involving liver, lung, breast, abdomen, chest wall, and pancreatic sites. All subjects were imaged during breath hold with CBCT on a standard image guidance platform (TrueBeam 2.7, Varian Medical Systems) and FBCT CT (CTsim, GE Optima). HyperSight imaging with both breath hold (HSBH) and free breathing (HSFB) was performed in a single session. The 4 image sets thus acquired were registered and compared using metrics quantifying artifact index, image nonuniformity, contrast, contrast-to-noise ratio, and difference of Hounsfield unit (HU) from CTsim.

RESULTS

HSBH provided less severe artifacts compared with both HSFB and TrueBeam. The severity of artifacts in HSBH images was similar to that in CTsim images, with statistically similar artifact index values. CTsim provided the best image uniformity; however, HSBH provided improved uniformity compared with both HSFB and TrueBeam. CTsim demonstrated elevated contrast compared with HyperSight imaging, but both HSBH and HSFB imaging showed superior contrast-to-noise ratio characteristics compared with TrueBeam. The median HU difference of HSBH from CTsim was within 1 HU for muscle/fat tissue, 12 HU for bone, and 14 HU for lung.

CONCLUSIONS

The HyperSight system provides 6-second CBCT acquisition with image artifacts that are significantly reduced compared with TrueBeam and comparable to those in CTsim FBCT imaging. HyperSight breath hold imaging was of higher quality compared with free breathing imaging on the same system. The median HU value in HyperSight breath hold imaging is within 15 HU of that in CTsim imaging for muscle, fat, bone, and lung tissue types, indicating the utility of image data for direct dose calculation in adaptive workflows.

Comparison and evaluation of different deep learning models of synthetic CT generation from CBCT for nasopharynx cancer adaptive proton therapy

BACKGROUND

Cone-beam computed tomography (CBCT) scanning is used for patient setup in image-guided radiotherapy. However, its inaccurate CT numbers limit its applicability in dose calculation and treatment planning.

PURPOSE

This study compares four deep learning methods for generating synthetic CT (sCT) to determine which method is more appropriate and offers potential for further clinical exploration in adaptive proton therapy for nasopharynx cancer.

METHODS

CBCTs and deformed planning CT (dCT) from 75 patients (60/5/10 for training, validation and testing) were used to compare cycle-consistent Generative Adversarial Network (cycleGAN), Unet, Unet+cycleGAN and conditionalGenerative Adversarial Network (cGAN) for sCT generation. The sCT images generated by each method were evaluated against dCT images using mean absolute error (MAE), structural similarity (SSIM), peak signal-to-noise ratio (PSNR), spatial non-uniformity (SNU) and radial averaging in the frequency domain. In addition, dosimetric accuracy was assessed through gamma analysis, differences in water equivalent thickness (WET), and dose-volume histogram metrics.

RESULTS

The cGAN model has demonstrated optimal performance in the four models across various indicators. In terms of image quality under global condition, the average MAE has been reduced to 16.39HU, SSIM has increased to 95.24%, and PSNR has increased to 28.98. Regarding dosimetric accuracy, the gamma passing rate (2%/2 mm) has reached 99.02%, and the WET difference is only 1.28 mm. The D95 value of CTVs coverage and Dmax value of spinal cord, brainstem show no significant differences between dCT and sCT generated by cGAN model.

CONCLUSIONS

The cGAN model has been shown to be a more suitable approach for generating sCT using CBCT, considering its characteristics and concepts. The resulting sCT has the potential for application in adaptive proton therapy.

Streaking artifact reduction for CBCT-based synthetic CT generation in adaptive radiotherapy

BACKGROUND

Cone-beam computed tomography (CBCT) is widely used for daily image guidance in radiation therapy, enhancing the reproducibility of patient setup. However, its application in adaptive radiotherapy (ART) is limited by many imaging artifacts and inaccurate Hounsfield units (HUs). The correction of CBCT image is necessary and of great value for CBCT-based ART.

PURPOSE

To explore the synthetic CT (sCT) generation from CBCT images of thorax and abdomen patients, which usually surfer from serious artifacts duo to organ state changes. In this study, a streaking artifact reduction network (SARN) is proposed to reduce artifacts and combine with cycleGAN to generate high-quality sCT images from CBCT and achieve an accurate dose calculation.

METHODS

The proposed SARN was trained in a self-supervised manner. Artifact-CT images were generated from planning CT by random deformation and projection replacement, and SARN was trained based on paired artifact-CT and CT images. The planning CT and CBCT images of 260 patients with cancer, including 120 thoracic and 140 abdominal CT scans, were used to train and evaluate neural networks. The CBCT images of another 12 patients in late treatment fractions, which contained large anatomy changes, were also tested by trained models. The trained models include commonly used U-Net, cycleGAN, attention-gated cycleGAN (cycAT), and cascade models combined SARN with cycleGAN or cycAT. The generated sCT images were compared in terms of image quality and dose calculation accuracy.

RESULTS

The sCT images generated by SARN combined with cycleGAN and cycAT showed the best image quality, removed the most artifacts, and retained the normal anatomical structure. The SARN+cycleGAN performed best in streaking artifacts removal with the maximum percent integrity uniformity (PIU_m ) of 91.0% and minimum standard deviation (SD) of 35.4 HU for delineated artifact regions among all models. The mean absolute error (MAE) of CBCT images in the thorax and abdomen were 71.6 and 55.2 HU, respectively, using planning CT images after deformable registration as ground truth. Compared with CBCT, the thoracic and abdominal sCT images generated by each model had significantly improved image quality with smaller MAE (p < 0.05). The SARN+cycAT obtained the minimum MAEs of 42.5 HU in the thorax while SARN+cycleGAN got the minimum MAEs of 32.0 HU in the abdomen. The sCT generated by U-Net had a remarkably lower anatomical structure accuracy compared with the other models. The thoracic and abdominal sCT images generated by SARN+cycleGAN showed optimal dose calculation accuracy with gamma passing rates (2 mm/2%) of 98.2% and 96.9%, respectively.

CONCLUSIONS

The proposed SARN can reduce serious streaking artifacts in CBCT images. The SARN combined with cycleGAN can generate high-quality sCT images with fewer artifacts, high-accuracy HU values, and accurate anatomical structures, thus providing reliable dose calculation in ART.

Region-of-interest intra-arc MV imaging to facilitate sub-mm positional accuracy with minimal imaging dose during treatment deliveries of small cranial lesions

Purpose

To automate the generation of region‐of‐interest (ROI) apertures for use with megavoltage imaging for online positional corrections during cranial stereotactic radiosurgery.

Materials and methods

Digitally reconstructed radiographs (DRRs) were created for a 3D‐printed skull phantom at 5 degree gantry angle increments for a three‐arc beam arrangement. At each angle, 3000 random rectangular apertures were generated, and 100 shifts on a grid were applied to the anatomy within the frame. For all shifts, the mutual information (MI) between the shifted and unshifted DRR was calculated to derive an average MI gradient. The top 10% of apertures that minimized registration errors were overlaid and discretely thresholded to generate imaging plans. Imaging was acquired with the skull while implementing simulated patient motion on a linac. Control point‐specific couch motions were derived to align the skull to its planned positioning.

Results

Apertures with a range of repositioning errors less than 0.1 mm possessed a 42% larger average MI gradient when compared with apertures with a range greater than 1 mm. Dose calculations with Monte Carlo exhibited an 84% reduction in the dose received by 50% of the skull with the 50% thresholded plan when compared to a constant 22 × 22 cm² imaging plan. For all different imaging plans (with and without motion), the calculated median 3D‐errors with respect to the tracking of a metal‐BB fiducial positioned at isocenter in the skull were sub‐mm except for the 80% thresholded plan.

Conclusions

Sub‐mm positional errors are achievable with couch motions derived from control point–specific ROI imaging. Smaller apertures that conform to an anatomical ROI can be utilized to minimize the imaging dose incurred at the expense of larger errors.

A novel hardware duo of beam modulation and shielding to reduce scatter acquisition and dose in cone-beam breast CT

PURPOSE

In cone-beam breast CT, scattered photons form a large portion of the acquired signal, adversely impacting image quality throughout the frequency response of the imaging system. Prior simulation studies provided proof of concept for utilization of a hardware solution to prevent scatter acquisition. Here, we report the design, implementation, and characterization of an auxiliary apparatus of fluence modulation and scatter shielding that does indeed lead to projections with a reduced level of scatter.

METHODS

An apparatus was designed for permanent installation within an existing cone-beam CT system. The apparatus is composed of two primary assemblies: a "Fluence Modulator" (FM) and a "Scatter Shield" (SS). The design of the assemblies enables them to operate in synchrony during image acquisition, converting the sourced x-rays into a moving narrow beam. During a projection, this narrow beam sweeps the entire fan angle coverage of the imaging system. As the two assemblies are contingent on one another, their joint implementation is described in the singular as apparatus FM-SS. The FM and the SS assemblies are each comprised a metal housing, a sensory system, and a robotic system. A controller unit handles their relative movements. A series of comparative studies were conducted to evaluate the performance of a cone-beam CT system in two "modes" of operation: with and without FM-SS installed, and to compare the results of physical implementation with those previously simulated. The dynamic range requirements of the utilized detector in the cone-beam CT imaging system were first characterized, independent of the mode of operation. We then characterized and compared the spatial resolution of the imaging system with, and without, FM-SS. A physical breast phantom, representative of an average size breast, was developed and imaged. Actual differences in signal level obtained with, versus without, FM-SS were then compared to the expected level gains based on previously reported simulations. Following these initial assessments, the scatter acquisition in each projection in both modes of operation was investigated. Finally, as an initial study of the impact of FM-SS on radiation dose in an average size breast, a series of Monte Carlo simulations were coupled with physical measurements of air kerma, with and without FM-SS.

RESULTS

With implementation of FM-SS, the detector's required dynamic range was reduced by a factor of 5.5. Substantial reduction in the acquisition of the scattered rays, by a factor of 5.1 was achieved. With the implementation of FM-SS, deposited dose was reduced by 27% in the studied breast.

CONCLUSIONS

The disclosed implementation of FM-SS, within a cone-beam breast CT system, results in reduction of scatter-components in acquired projections, reduction of dose deposit to the breast, and relaxation of requirements for the detector's dynamic range. Controlling or correcting for patient motion occurring during image acquisition remains an open problem to be solved prior to practical clinical usage of FM-SS cone-beam breast CT.

Intraoperative cone-beam and slot-beam CT: 3D image quality and dose with a slot collimator on the O-arm imaging system

PURPOSE

To characterize the 3D imaging performance and radiation dose for a prototype slot-beam configuration on an intraoperative O-arm™ Surgical Imaging System (Medtronic Inc., Littleton, MA) and identify potential improvements in soft-tissue image quality for surgical interventions.

METHODS

A slot collimator was integrated with the O-arm^™ system for slot-beam axial CT. The collimator can be automatically actuated to provide 1.2° slot-beam longitudinal collimation. Cone-beam and slot-beam configurations were investigated with and without an antiscatter grid (12:1 grid ratio, 60 lines/cm). Dose, scatter, image noise, and soft-tissue contrast resolution were evaluated in quantitative phantoms for head and body configurations over a range of exposure levels (beam energy and mAs), with reconstruction performed via filtered-backprojection. Qualitative imaging performance across various anatomical sites and imaging tasks was assessed with anthropomorphic head, abdomen, and pelvis phantoms.

RESULTS

The dose for a slot-beam scan varied from 0.02-0.06 mGy/mAs for head protocols to 0.01-0.03 mGy/mAs for body protocols, yielding dose reduction by ∼1/5 to 1/3 compared to cone-beam, owing to beam collimation and reduced x-ray scatter. The slot-beam provided an ∼6-7× reduction in scatter-to-primary ratio (SPR) compared to the cone-beam, yielding SPR ∼20-80% for head and body without the grid and ∼7-30% with the grid. Compared to cone-beam scans at equivalent dose, slot-beam images exhibited an ∼2.5× increase in soft-tissue contrast-to-noise ratio (CNR) for both grid and gridless configurations. For slot-beam scans, a further ∼10-30% improvement in CNR was achieved when the grid was removed. Slot-beam imaging could benefit certain interventional scenarios in which improved visualization of soft tissues is required within a fairly narrow longitudinal region of interest ( ± 7 mm in z )--for example, checking the completeness of tumor resection, preservation of adjacent anatomy, or detection of complications (e.g., hemorrhage). While preserving existing capabilities for fluoroscopy and cone-beam CT, slot-beam scanning could enhance the utility of intraoperative imaging and provide a useful mode for safety and validation checks in image-guided surgery.

CONCLUSIONS

The 3D imaging performance and dose of a prototype slot-beam CT configuration on the O-arm^™ system was investigated. Substantial improvements in soft-tissue image quality and reduction in radiation dose are evident with the slot-beam configuration due to reduced x-ray scatter.

Scatter Correction for Spectral CT Using a Primary Modulator Mask

The problem of scattered radiation correction in computed tomography (CT) is well known because scatter induces a bias, a loss of contrast and artifacts. Numerous strategies have been proposed in conventional CT (using energy-integrating detectors) but the problem is still open in the field of spectral CT, a new imaging technique based on energy-selective photon counting detectors. The aim of the present study is to introduce a scatter correction method adapted to multi-energy imaging and based on the use of a primary modulator mask. The main contributions are a correction matrix, which compensates for the effect of the mask, a scatter model based on B-splines and a cost function based on the mask structures and robust to the object structures. The performances of the method have been evaluated on both simulated and experimental data. The mean relative error was reduced from 20% in the lower energy-bins without correction to 4% with the proposed technique, which is close to the error caused by statistical noise.

CBCT-based synthetic CT generation using deep-attention cycleGAN for pancreatic adaptive radiotherapy

PURPOSE

Current clinical application of cone-beam CT (CBCT) is limited to patient setup. Imaging artifacts and Hounsfield unit (HU) inaccuracy make the process of CBCT-based adaptive planning presently impractical. In this study, we developed a deep-learning-based approach to improve CBCT image quality and HU accuracy for potential extended clinical use in CBCT-guided pancreatic adaptive radiotherapy.

METHODS

Thirty patients previously treated with pancreas SBRT were included. The CBCT acquired prior to the first fraction of treatment was registered to the planning CT for training and generation of synthetic CT (sCT). A self-attention cycle generative adversarial network (cycleGAN) was used to generate CBCT-based sCT. For the cohort of 30 patients, the CT-based contours and treatment plans were transferred to the first fraction CBCTs and sCTs for dosimetric comparison.

RESULTS

At the site of abdomen, mean absolute error (MAE) between CT and sCT was 56.89 ± 13.84 HU, comparing to 81.06 ± 15.86 HU between CT and the raw CBCT. No significant differences (P > 0.05) were observed in the PTV and OAR dose-volume-histogram (DVH) metrics between the CT- and sCT-based plans, while significant differences (P < 0.05) were found between the CT- and the CBCT-based plans.

CONCLUSIONS

The image similarity and dosimetric agreement between the CT and sCT-based plans validated the dose calculation accuracy carried by sCT. The CBCT-based sCT approach can potentially increase treatment precision and thus minimize gastrointestinal toxicity.

Shading artifact correction in breast CT using an interleaved deep learning segmentation and maximum-likelihood polynomial fitting approach

PURPOSE

The purpose of this work was twofold: (a) To provide a robust and accurate method for image segmentation of dedicated breast CT (bCT) volume data sets, and (b) to improve Hounsfield unit (HU) accuracy in bCT by means of a postprocessing method that uses the segmented images to correct for the low-frequency shading artifacts in reconstructed images.

METHODS

A sequential and iterative application of image segmentation and low-order polynomial fitting to bCT volume data sets was used in the interleaved correction (IC) method. Image segmentation was performed through a deep convolutional neural network (CNN) with a modified U-Net architecture. A total of 45 621 coronal bCT images from 111 patient volume data sets were segmented (using a previously published segmentation algorithm) and used for neural network training, validation, and testing. All patient data sets were selected from scans performed on four different prototype breast CT systems. The adipose voxels for each patient volume data set, segmented using the proposed CNN, were then fit to a three-dimensional low-order polynomial. The polynomial fit was subsequently used to correct for the shading artifacts introduced by scatter and beam hardening in a method termed "flat fielding." An interleaved utilization of image segmentation and flat fielding was repeated until a convergence criterion was satisfied. Mathematical and physical phantom studies were conducted to evaluate the dependence of the proposed algorithm on breast size and the distribution of fibroglandular tissue. In addition, a subset of patient scans (not used in the CNN training, testing or validation) were used to investigate the accuracy of the IC method across different scanner designs and beam qualities.

RESULTS

The IC method resulted in an accurate classification of different tissue types with an average Dice similarity coefficient > 95%, precision > 97%, recall > 95%, and F1-score > 96% across all tissue types. The flat fielding correction of bCT images resulted in a significant reduction in either cupping or capping artifacts in both mathematical and physical phantom studies as measured by the integral nonuniformity metric with an average reduction of 71% for cupping and 30% for capping across different phantom sizes, and the Uniformity Index with an average reduction of 53% for cupping and 34% for capping.

CONCLUSION

The validation studies demonstrated that the IC method improves Hounsfield Units (HU) accuracy and effectively corrects for shading artifacts caused by scatter contamination and beam hardening. The postprocessing approach described herein is relevant to the broad scope of bCT devices and does not require any modification in hardware or existing scan protocols. The trained CNN parameters and network architecture are available for interested users.

Characterization of X-ray scattering for various phantoms and clinical breast geometries using breast CT on a dedicated hybrid system

OBJECTIVE

The purpose of this study was to utilize a dedicated breast CT system using a 2D beam stop array to physically evaluate the scatter to primary ratios (SPRs) of different geometric phantoms and prospectively acquired clinical patient data.

METHODS

Including clinically unrealizable compositions of 100% glandular and 100% fat, projection images were acquired using three geometrically different phantoms filled with fluids simulating breast tissue. The beam stop array method was used for measuring scatter in projection space, and creating the scatter corrected primary images. 2D SPRs were calculated. Additionally, a new figure of merit, the 3D normalized scatter contribution (NSC) volumes were calculated.

RESULTS

The 2D SPR values (0.52-1.10) were primarily dependent on phantom geometry; a secondary dependence was due to their uniform density; 2D SPRs were low frequency and smoothly varying in the uniformly filled phantoms. SPRs of clinical patient data followed similar trends as phantoms, but with noticeable deviations and high frequency components due to the heterogeneous distribution of glandular tissue. The maximum measured patient 2D SPRs were all <0.6, even for the largest diameter breast. These results demonstrate modest scatter components with changing object geometries and densities; the 3D NSC volumes with higher frequency components help visualize scatter distribution throughout the reconstructed image volumes. Furthermore, the SPRs in the heterogeneous clinical breast cases were underestimated by the equivalent density, uniformly filled phantoms.

CONCLUSIONS

These results provide guidance on the use of uniformly distributed density and differently shaped phantoms when considering simulations. They also clearly demonstrate that results from patients can vary considerably from 2D SPRs of uniformly simulated phantoms.

Interior tomographic imaging of mouse heart in a carbon nanotube micro-CT

BACKGROUND

The relatively high radiation dose from micro-CT is a cause for concern in preclinical research involving animal subjects. Interior region-of-interest (ROI) imaging was proposed for dose reduction, but has not been experimentally applied in micro-CT.

OBJECTIVE

Our aim is to implement interior ROI imaging in a carbon nanotube (CNT) x-ray source based micro-CT, and present the ROI image quality and radiation dose reduction for interior cardiac micro-CT imaging of a mouse heart in situ.

METHODS

An aperture collimator was mounted at the source-side to induce a small-sized cone beam (10 mm width) at the isocenter. Interior in situ micro-CT scans were conducted on a mouse carcass and several micro-CT phantoms. A GPU-accelerated hybrid iterative reconstruction algorithm was employed for volumetric image reconstruction. Radiation dose was measured for the same system operated at the interior and global micro-CT modes.

RESULTS

Visual inspection demonstrated comparable image quality between two scan modes. Quantitative evaluation demonstrated high structural similarity index (up to 0.9614) with improved contrast-noise-ratio (CNR) on interior micro-CT mode. Interior micro-CT mode yielded significant reduction (up to 83.9%) for dose length product (DLP).

CONCLUSIONS

This work demonstrates the applicability of using CNT x-ray source based interior micro-CT for preclinical imaging with significantly reduced radiation dose.

Radiation doses in volume-of-interest breast computed tomography--A Monte Carlo simulation study

PURPOSE

Cone beam breast computed tomography (breast CT) with true three-dimensional, nearly isotropic spatial resolution has been developed and investigated over the past decade to overcome the problem of lesions overlapping with breast anatomical structures on two-dimensional mammographic images. However, the ability of breast CT to detect small objects, such as tissue structure edges and small calcifications, is limited. To resolve this problem, the authors proposed and developed a volume-of-interest (VOI) breast CT technique to image a small VOI using a higher radiation dose to improve that region's visibility. In this study, the authors performed Monte Carlo simulations to estimate average breast dose and average glandular dose (AGD) for the VOI breast CT technique.

METHODS

Electron-Gamma-Shower system code-based Monte Carlo codes were used to simulate breast CT. The Monte Carlo codes estimated were validated using physical measurements of air kerma ratios and point doses in phantoms with an ion chamber and optically stimulated luminescence dosimeters. The validated full cone x-ray source was then collimated to simulate half cone beam x-rays to image digital pendant-geometry, hemi-ellipsoidal, homogeneous breast phantoms and to estimate breast doses with full field scans. 13-cm in diameter, 10-cm long hemi-ellipsoidal homogeneous phantoms were used to simulate median breasts. Breast compositions of 25% and 50% volumetric glandular fractions (VGFs) were used to investigate the influence on breast dose. The simulated half cone beam x-rays were then collimated to a narrow x-ray beam with an area of 2.5 × 2.5 cm(2) field of view at the isocenter plane and to perform VOI field scans. The Monte Carlo results for the full field scans and the VOI field scans were then used to estimate the AGD for the VOI breast CT technique.

RESULTS

The ratios of air kerma ratios and dose measurement results from the Monte Carlo simulation to those from the physical measurements were 0.97 ± 0.03 and 1.10 ± 0.13, respectively, indicating that the accuracy of the Monte Carlo simulation was adequate. The normalized AGD with VOI field scans was substantially reduced by a factor of about 2 over the VOI region and by a factor of 18 over the entire breast for both 25% and 50% VGF simulated breasts compared with the normalized AGD with full field scans. The normalized AGD for the VOI breast CT technique can be kept the same as or lower than that for a full field scan with the exposure level for the VOI field scan increased by a factor of as much as 12.

CONCLUSIONS

The authors' Monte Carlo estimates of normalized AGDs for the VOI breast CT technique show that this technique can be used to markedly increase the dose to the breast and thus the visibility of the VOI region without increasing the dose to the breast. The results of this investigation should be helpful for those interested in using VOI breast CT technique to image small calcifications with dose concern.

Volume-of-interest imaging of the inner ear in a human temporal bone specimen using a robot- driven C-arm flat panel detector CT system

VOI imaging can provide higher image quality at a reduced dose for a subregion. In this study with a robot-driven C-arm FDCT system, the goals were proof of feasibility for inner ear imaging, higher flexibility during data acquisition, and easier processing during reconstruction. First a low-dose OV scan was acquired allowing an orientation and enabling the selection of the VOI. The C-arm was then moved by the robotic system without a need for patient movement and the VOI was scanned with adapted parameters. Uncompromised artifact-free image quality was achieved by the 2-scan approach and the dose was reduced by 80%-90% in comparison with conventional MSCT and FPCT scans.

High resolution dual detector volume-of-interest cone beam breast CT--Demonstration with a bench top system

PURPOSE

In this study, we used a small field high resolution detector in conjunction with a full field flat panel detector to implement and investigate the dual detector volume-of-interest (VOI) cone beam breast computed tomography (CBCT) technique on a bench-top system. The potential of using this technique to image small calcifications without increasing the overall dose to the breast was demonstrated. Significant reduction of scatter components in the high resolution projection image data of the VOI was also shown.

METHODS

With the regular flat panel based CBCT technique, exposures were made at 80 kVp to generate an air kerma of 6 mGys at the isocenter. With the dual detector VOI CBCT technique, a high resolution small field CMOS detector was used to scan a cylindrical VOI (2.5 cm in diameter and height, 4.5 cm off-center) with collimated x-rays at four times of regular exposure level. A flat panel detector was used for full field scan with low x-ray exposures at half of the regular exposure level. The low exposure full field image data were used to fill in the truncated space in the VOI scan data and generate a complete projection image set. The Feldkamp-Davis-Kress (FDK) filtered backprojection algorithm was used to reconstruct high resolution images for the VOI. Two scanning techniques, one breast centered and the other VOI centered, were implemented and investigated. Paraffin cylinders with embedded thin aluminum (Al) wires were imaged and used in conjunction with optically stimulated luminescence (OSL) dose measurements to demonstrate the ability of this technique to image small calcifications without increasing the mean glandular dose (MGD).

RESULTS

Using exposures that produce an air kerma of 6 mGys at the isocenter, the regular CBCT technique was able to resolve the cross-sections of Al wires as thin as 254 μm in diameter in the phantom. For the specific VOI studied, by increasing the exposure level by a factor of 4 for the VOI scan and reducing the exposure level by a factor of 2 for the full filed scan, the dual-detector CBCT technique was able to resolve the cross-sections of Al wires as thin as 152 μm in diameter. The CNR evaluated for the entire Al wire cross-section was found to be improved from 5.5 in regular CBCT to 14.4 and 16.8 with the breast centered and VOI centered scanning techniques, respectively. Even inside VOI center, the VOI scan resulted in significant dose saving with the dose reduced by a factor of 1.6 at the VOI center. Dose saving outside the VOI was substantial with the dose reduced by a factor of 7.3 and 7.8 at the breast center for the breast centered and VOI centered scans, respectively, when compared to full field scan at the same exposure level. The differences between the two dual detector techniques in terms of dose saving and scatter reduction were small with VOI scan at 4× exposure level and full field scan at 0.5 × exposure level. The MGDs were only 94% of that from the regular CBCT scan.

CONCLUSIONS

For the specific VOI studied, the dual detector VOI CBCT technique has the potential to provide high quality images inside the VOI with MGD similar to or even lower than that of full field breast CBCT. It was also found that our results were compromised by the use of inadequate detectors for the VOI scan. An appropriately selected detector would better optimize the image quality improvement that can be achieved with the VOI CBCT technique.

Region-of-interest material decomposition from truncated energy-resolved CT

PURPOSE

Energy-resolved CT using photon-counting detectors has the potential to provide improved material decomposition compared to dual-kVp approaches. However, available photon-counting detectors are susceptible to pulse-pileup artifacts, especially at the periphery of the field of view (FOV) where the object attenuation is low compared to the center of the FOV. Pulse pileup may be avoided by imaging a region-of-interest (ROI) where the dynamic range is expected to be limited. This work investigated performing material decomposition and reconstructing ROI basis images from truncated energy-resolved data.

METHODS

A method is proposed to reconstruct images of basis functions primarily contained within the ROI, such as targeted or localized K-edge contrast agents. Material decomposition is performed independently for each ray in the sinogram, followed by filtered backprojection from the truncated data encompassing the ROI. A second method is proposed that uses a prior conventional energy-integrating image to estimate energy-resolved data outside the ROI. The measured and estimated energy-resolved data are decomposed into basis projections and merged into basis sinograms of the full FOV. Basis images of the ROI are then reconstructed through filtered backprojection. This method is most easily applied to objects that do not contain K-edge contrast agents outside the ROI. Simulations of a voxelized thorax phantom with iodine in the blood pool and a detector with five energy bins were performed. Full FOV, truncated, and truncated data merged with data estimated from the prior energy-integrating image were decomposed into Compton, photoelectric, and iodine basis functions. An empirical weighting factor was determined to blend the merged sinogram at the boundary of the truncated data. The effects of noise and misalignment in the prior image were also quantified. Basis images of the central 15 cm × 15 cm ROI containing the heart were reconstructed via filtered backprojection. Basis image accuracy was quantified relative to gold-standard basis images reconstructed from full FOV energy-resolved data.

RESULTS

The error in the iodine basis image reconstructed from truncated energy-resolved data without prior information was less than 1% for the central 7 cm of the 7.5-cm-radius ROI and 3% at the edge of the ROI. When the truncated and estimated basis sinograms were blended, the error was below 1% throughout the ROI for photoelectric basis images and ranged from 1% at the center of the ROI to 4% at the edge for the Compton basis image.

CONCLUSIONS

The density of localized K-edge contrast agents can be estimated to within 1% error using filtered back projection without prior information. For noncontrast and localized-contrast scans, ROI images of general basis functions can be reconstructed to within a few percent error using a prior energy-integrating image. The ability to perform material decomposition for a limited ROI may facilitate energy-resolved CT with available photon-counting detectors.

A general framework and review of scatter correction methods in x-ray cone-beam computerized tomography. Part 1: Scatter compensation approaches

Since scattered radiation in cone-beam volume CT implies severe degradation of CT images by quantification errors, artifacts, and noise increase, scatter suppression is one of the main issues related to image quality in CBCT imaging. The aim of this review is to structurize the variety of scatter suppression methods, to analyze the common structure, and to develop a general framework for scatter correction procedures. In general, scatter suppression combines hardware techniques of scatter rejection and software methods of scatter correction. The authors emphasize that scatter correction procedures consist of the main components scatter estimation (by measurement or mathematical modeling) and scatter compensation (deterministic or statistical methods). The framework comprises most scatter correction approaches and its validity also goes beyond transmission CT. Before the advent of cone-beam CT, a lot of papers on scatter correction approaches in x-ray radiography, mammography, emission tomography, and in Megavolt CT had been published. The opportunity to avail from research in those other fields of medical imaging has not yet been sufficiently exploited. Therefore additional references are included when ever it seems pertinent. Scatter estimation and scatter compensation are typically intertwined in iterative procedures. It makes sense to recognize iterative approaches in the light of the concept of self-consistency. The importance of incorporating scatter compensation approaches into a statistical framework for noise minimization has to be underscored. Signal and noise propagation analysis is presented. A main result is the preservation of differential-signal-to-noise-ratio (dSNR) in CT projection data by ideal scatter correction. The objective of scatter compensation methods is the restoration of quantitative accuracy and a balance between low-contrast restoration and noise reduction. In a synopsis section, the different deterministic and statistical methods are discussed with respect to their properties and applications. The current paper is focused on scatter compensation algorithms. The multitude of scatter estimation models will be dealt with in a separate paper.

Reconstruction of a cone-beam CT image via forward iterative projection matching

PURPOSE

To demonstrate the feasibility of reconstructing a cone-beam CT (CBCT) image by deformably altering a prior fan-beam CT (FBCT) image such that it matches the anatomy portrayed in the CBCT projection data set.

METHODS

A prior FBCT image of the patient is assumed to be available as a source image. A CBCT projection data set is obtained and used as a target image set. A parametrized deformation model is applied to the source FBCT image, digitally reconstructed radiographs (DRRs) that emulate the CBCT projection image geometry are calculated and compared to the target CBCT projection data, and the deformation model parameters are adjusted iteratively until the DRRs optimally match the CBCT projection data set. The resulting deformed FBCT image is hypothesized to be an accurate representation of the patient's anatomy imaged by the CBCT system. The process is demonstrated via numerical simulation. A known deformation is applied to a prior FBCT image and used to create a synthetic set of CBCT target projections. The iterative projection matching process is then applied to reconstruct the deformation represented in the synthetic target projections; the reconstructed deformation is then compared to the known deformation. The sensitivity of the process to the number of projections and the DRR/CBCT projection mismatch is explored by systematically adding noise to and perturbing the contrast of the target projections relative to the iterated source DRRs and by reducing the number of projections.

RESULTS

When there is no noise or contrast mismatch in the CBCT projection images, a set of 64 projections allows the known deformed CT image to be reconstructed to within a nRMS error of 1% and the known deformation to within a nRMS error of 7%. A CT image nRMS error of less than 4% is maintained at noise levels up to 3% of the mean projection intensity, at which the deformation error is 13%. At 1% noise level, the number of projections can be reduced to 8 while maintaining CT image and deformation errors of less than 4% and 13%, respectively. The method is sensitive to contrast mismatch between the simulated projections and the target projections when the soft-tissue contrast in the projections is low.

CONCLUSIONS

By using prior knowledge available in a FBCT image, the authors show that a CBCT image can be iteratively reconstructed from a comparatively small number of projection images, thus saving acquisition time and reducing imaging dose. This will enable more frequent daily imaging during radiation therapy. Because the process preserves the CT numbers of the FBCT image, the resulting 3D image intensities will be more accurate than a CBCT image reconstructed via conventional backprojection methods. Reconstruction errors are insensitive to noise at levels beyond what would typically be found in CBCT projection data, but are sensitive to contrast mismatch errors between the CBCT projection data and the DRRs.

Volume-of-interest (VOI) imaging in C-arm flat-detector CT for high image quality at reduced dose

PURPOSE

A novel method for flat-detector computed tomography was developed to enable volume-of-interest (VOI) imaging at high resolution, low noise, and reduced dose. For this, a full low-dose overview (OV) scan and a local high-dose scan of a VOI are combined.

METHODS

The first scan yields an overview of the whole object and enables the selection of an arbitrary VOI. The second scan of that VOI assures high image quality within the interesting volume. The combination of the two consecutive scans is based on a forward projection of the reconstructed OV volume that was registered to the VOI. The artificial projection data of the OV scan are combined with the measured VOI data in the raw data domain. Different projection values are matched by an appropriate transformation and weighting. The reconstruction is performed with a standard Feldkamp-type algorithm. In simulations, the combination of OV scan and VOI scan was investigated on a mathematically described phantom. In measurements, spatial resolution and noise were evaluated with image quality phantoms. Modulation transfer functions and noise values were calculated. Measurements of an anthropomorphic head phantom were used to validate the proposed method for realistic applications, e.g., imaging stents. In Monte Carlo simulations, 3D dose distributions were calculated and dose values were assessed quantitatively.

RESULTS

By the proposed combination method, an image is generated which covers the whole object and provides the VOI at high image quality. In the OV image, a resolution of 0.7 lp/mm (line pairs per millimeter) and noise of 63.5 HU were determined. Inside the VOI, resolution was increased to 2.4 lp/mm and noise was decreased to 18.7 HU. For the performed measurements, the cumulative dose was significantly reduced in comparison to conventional scans by up to 93%. The dose of a high-quality scan, for example, was reduced from 97 to less than 7 mGy, while keeping image quality constant within the VOI.

CONCLUSIONS

The proposed VOI application with two scans is an effective way to ensure high image quality within the VOI while simultaneously reducing the cumulative patient dose.