Evaluation of detector readout gain mode and bowtie filters for cone-beam CT imaging of the head

Algorithm development for automatic laser alignment assessment on an ACR CT phantom and its evaluation on sixteen CT scanners

Purpose. The aim of this study is to develop software to automatically assess the laser alignment on the ACR CT phantom and evaluate its accuracy on sixteen CT scanners.Methods. Software for an automated method of laser alignment assessment on the ACR CT phantom was developed. Laser alignment assessment was based on the positions of the ball-bearing markers at the edge of the ACR CT phantom. The automatic assessment was performed using several steps, including segmentation to acquire the coordinates of the ball-bearing markers and determination of the distances between lines connecting them with lines through the center of the image. A comparison of the results from the automatic method with those from the manual method was performed. The manual measurements were carried out using MicroDicom Viewer. A Mann-Whitney U test was performed to determine the statistical difference between both methods. The evaluation was performed on images of the ACR CT phantom scanned with 16 CT scanners from 5 different CT manufacturers.Results. The results confirmed that our software successfully segments the ball-bearing markers and determines the laser alignment assessment on the ACR CT phantom. Evaluation of the algorithm with images from the 16 CT scanners revealed that the difference between the results from automatic and manual methods were about 0.2 mm with apvalue of around 0.7 (no statistical difference). Misalignment in they-axis was larger than the misalignment in the x-axisfor the majority of the scanners tested. It was found that the phantom tended to be placed 2 mm higher than the iso-center.Conclusions. Software to automatically assess CT laser alignment with the ACR CT phantom was successfully developed and evaluated. The automatic assessment was comparable to manual assessment. In addition, the automatic method was user independent and fast.

Scatter-to-primary ratio in dentomaxillofacial cone-beam CT: effect of field of view and beam energy

OBJECTIVE

The aim of this study was to evaluate the effect of field of view (FOV) and beam energy on the scatter-to-primary ratio (SPR) in dental cone-beam CT (CBCT).

METHODS

An anthropomorphic phantom representing an adult male (ATOM Max 711-HN, Norfolk, VA, USA) was scanned using the 3D Accuitomo 170 CBCT (J. Morita, Kyoto, Japan) using 11 FOVs. During each scan, half of the X-ray beam was blocked. Each scan was performed at three exposure settings with varying beam energy and equal radiation dose: 90 kV 5 mA, 77 kV 7.5 mA and 69 kV 10 mA. The SPR was estimated by measuring the grey values in the blocked and non-blocked regions of the RAW data. The effect of FOV on SPR was evaluated using Dunn's multiple comparison test, and the effect of the exposure settings was compared using a Wilcoxon signed rank test.

RESULTS

Larger FOVs showed increased scatter. FOVs with a shorter isocenter-detector distance showed a particularly high SPR. Most intercomparisons between FOVs were statistically significant. The largest difference was found between 17 × 12 cm and 6 × 6 cm (lower jaw), with the former showing a 4.9-fold higher SPR. The effect of beam energy was relatively small and varied between FOV sizes and positions.

CONCLUSION

While the choice of FOV size and position is determined by the diagnostic region of interest, the image quality deterioration for large FOVs due to scatter provides another incentive to limit the FOV size as much as possible.

3D Non-Rigid Alignment of Low-Dose Scans Allows to Correct for Saturation in Lower Extremity Cone-Beam CT

Detector saturation in cone-beam computed tomography occurs when an object of highly varying shape and material composition is imaged using an automatic exposure control (AEC) system. When imaging a subject's knees, high beam energy ensures the visibility of internal structures but leads to overexposure in less dense border regions. In this work, we propose to use an additional low-dose scan to correct the saturation artifacts of AEC scans. Overexposed pixels are identified in the projection images of the AEC scan using histogram-based thresholding. The saturation-free pixels from the AEC scan are combined with the skin border pixels of the low-dose scan prior to volumetric reconstruction. To compensate for patient motion between the two scans, a 3D non-rigid alignment of the projection images in a backward-forward-projection process based on fiducial marker positions is proposed. On numerical simulations, the projection combination improved the structural similarity index measure from 0.883 to 0.999. Further evaluations were performed on two in vivo subject knee acquisitions, one without and one with motion between the AEC and low-dose scans. Saturation-free reference images were acquired using a beam attenuator. The proposed method could qualitatively restore the information of peripheral tissue structures. Applying the 3D non-rigid alignment made it possible to use the projection images with inter-scan subject motion for projection image combination. The increase in radiation exposure due to the additional low-dose scan was found to be negligibly low. The presented methods allow simple but effective correction of saturation artifacts.

Cone-beam imaging with tilted rotation axis: Method and performance evaluation

PURPOSE

The recently introduced robotic x-ray systems provide the flexibility to acquire cone-beam computed tomography (CBCT) data using customized, application-specific source-detector trajectories. We exploit this capability to mitigate the effects of x-ray scatter and noise in CBCT imaging of weight-bearing foot and cervical spine (C-spine) using scan orbits with a tilted rotation axis.

METHODS

We used an advanced CBCT simulator implementing accurate models of x-ray scatter, primary attenuation, and noise to investigate the effects of the orbital tilt angle in upright foot and C-spine imaging. The system model was parameterized using a laboratory version of a three-dimensional (3D) robotic x-ray system (Multitom RAX, Siemens Healthineers). We considered a generalized tilted axis scan configuration, where the detector remained parallel to patient's long body axis during the acquisition, but the elevation of source and detector was changing. A modified Feldkamp-Davis-Kress (FDK) algorithm was developed for reconstruction in this configuration, which departs from the FDK assumption of a detector that is perpendicular to the scan plane. The simulated foot scans involved source-detector distance (SDD) of 1386 mm, orbital tilt angles ranging 10° to 40°, and 400 views at 1 mAs/view and 0.5° increment; the C-spine scans involved -25° to -45° tilt angles, SDD of 1090 mm, and 202 views at 1.3 mAs and 1° increment The imaging performance was assessed by projection-domain measurements of the scatter-to-primary ratio (SPR) and by reconstruction-domain measurements of contrast, noise and generalized contrast-to-noise ratio (gCNR, accounting for both image noise and background nonuniformity) of the metatarsals (foot imaging) and cervical vertebrae (spine imaging). The effects of scatter correction were also compared for horizontal and tilted scans using an ideal Monte Carlo (MC)-based scatter correction and a frame-by-frame mean scatter correction.

RESULTS

The proposed modified FDK, involving projection resampling, mitigated streak artifacts caused by the misalignment between the filtering direction and the detector rows. For foot imaging (no grids), an optimized 20° tilted orbit reduced the maximum SPR from ~1.5 in a horizontal scan to <0.5. The gCNR of the second metatarsal was enhanced twofold compared to a horizontal orbit. For the C-spine (with vertical grids), imaging with a tilted orbit avoided highly attenuating x-ray paths through the lower cervical vertebrae and shoulders. A -35° tilted orbit yielded improved image quality and visualization of the lower cervical spine: the SPR of lower cervical vertebrae was reduced from ~10 (horizontal orbit) to <6 (tilted orbit), and the gCNR for C5-C7 increased by a factor of 2. Furthermore, tilted orbits showed potential benefits over horizontal orbits by enabling scatter correction with a simple frame-by-frame mean correction without substantial increase in noise-induced artifacts after the correction.

CONCLUSIONS

Tilted scan trajectories, enabled by the emerging robotic x-ray system technology, were optimized for CBCT imaging of foot and cervical spine using an advanced simulation framework. The results demonstrated the potential advantages of tilted axis orbits in mitigation of scatter artifacts and improving contrast-to-noise ratio in CBCT reconstructions.

Cone-beam CT for imaging of the head/brain: Development and assessment of scanner prototype and reconstruction algorithms

PURPOSE

Our aim was to develop a high-quality, mobile cone-beam computed tomography (CBCT) scanner for point-of-care detection and monitoring of low-contrast, soft-tissue abnormalities in the head/brain, such as acute intracranial hemorrhage (ICH). This work presents an integrated framework of hardware and algorithmic advances for improving soft-tissue contrast resolution and evaluation of its technical performance with human subjects.

METHODS

Four configurations of a CBCT scanner prototype were designed and implemented to investigate key aspects of hardware (including system geometry, antiscatter grid, bowtie filter) and technique protocols. An integrated software pipeline (c.f., a serial cascade of algorithms) was developed for artifact correction (image lag, glare, beam hardening and x-ray scatter), motion compensation, and three-dimensional image (3D) reconstruction [penalized weighted least squares (PWLS), with a hardware-specific statistical noise model]. The PWLS method was extended in this work to accommodate multiple, independently moving regions with different resolution (to address both motion compensation and image truncation). Imaging performance was evaluated quantitatively and qualitatively with 41 human subjects in the neurosciences critical care unit (NCCU) at our institution.

RESULTS

The progression of four scanner configurations exhibited systematic improvement in the quality of raw data by variations in system geometry (source-detector distance), antiscatter grid, and bowtie filter. Quantitative assessment of CBCT images in 41 subjects demonstrated: ~70% reduction in image nonuniformity with artifact correction methods (lag, glare, beam hardening, and scatter); ~40% reduction in motion-induced streak artifacts via the multi-motion compensation method; and ~15% improvement in soft-tissue contrast-to-noise ratio (CNR) for PWLS compared to filtered backprojection (FBP) at matched resolution. Each of these components was important to improve contrast resolution for point-of-care cranial imaging.

CONCLUSIONS

This work presents the first application of a high-quality, point-of-care CBCT system for imaging of the head/ brain in a neurological critical care setting. Hardware configuration iterations and an integrated software pipeline for artifacts correction and PWLS reconstruction mitigated artifacts and noise to achieve image quality that could be valuable for point-of-care detection and monitoring of a variety of intracranial abnormalities, including ICH and hydrocephalus.

Dynamic fluence field modulation for miscentered patients in computed tomography

Traditional CT image acquisition uses bowtie filters to reduce dose, x-ray scatter, and detector dynamic range requirements. However, accurate patient centering within the bore of the CT scanner takes time and is often difficult to achieve precisely. Patient miscentering combined with a static bowtie filter can result in significant increases in dose, reconstruction noise, and CT number variations, and consequently raise overall exposure requirements. Approaches to estimate the patient position from scout scans and perform dynamic spatial beam filtration during acquisition are developed and applied in physical experiments on a CT test bench using different beam filtration strategies. While various dynamic beam modulation strategies have been developed, we focus on two approaches: (1) a simple approach using attenuation-based beam modulation using a translating bowtie filter and (2) dynamic beam modulation using multiple aperture devices (MADs)-an emerging beam filtration strategy based on binary filtration of the x-ray beam using variable width slits in a high-density beam blocker. Improved dose utilization and more consistent image performance with respect to an unmodulated baseline (static filter) are demonstrated for miscentered objects and dynamic beam filtration in physical experiments. For a homogeneous object miscentered by 4 cm, the dynamic filter reduced the maximum regional noise and dose penalties (compared with a centered object) from 173% to 16% and 42% to 14%, respectively, for a traditional bowtie, 29% to 8% and 24% to 15%, respectively, for a single MAD, and 275% to 11% and 56% to 18%, respectively, for a dual-MAD filter. The proposed methodology has the potential to relax patient centering requirements within the scanner, reduce setup time, and facilitate additional CT dose reduction.

Task-based statistical image reconstruction for high-quality cone-beam CT

Task-based analysis of medical imaging performance underlies many ongoing efforts in the development of new imaging systems. In statistical image reconstruction, regularization is often formulated in terms to encourage smoothness and/or sharpness (e.g. a linear, quadratic, or Huber penalty) but without explicit formulation of the task. We propose an alternative regularization approach in which a spatially varying penalty is determined that maximizes task-based imaging performance at every location in a 3D image. We apply the method to model-based image reconstruction (MBIR-viz., penalized weighted least-squares, PWLS) in cone-beam CT (CBCT) of the head, focusing on the task of detecting a small, low-contrast intracranial hemorrhage (ICH), and we test the performance of the algorithm in the context of a recently developed CBCT prototype for point-of-care imaging of brain injury. Theoretical predictions of local spatial resolution and noise are computed via an optimization by which regularization (specifically, the quadratic penalty strength) is allowed to vary throughout the image to maximize local task-based detectability index ([Formula: see text]). Simulation studies and test-bench experiments were performed using an anthropomorphic head phantom. Three PWLS implementations were tested: conventional (constant) penalty; a certainty-based penalty derived to enforce constant point-spread function, PSF; and the task-based penalty derived to maximize local detectability at each location. Conventional (constant) regularization exhibited a fairly strong degree of spatial variation in [Formula: see text], and the certainty-based method achieved uniform PSF, but each exhibited a reduction in detectability compared to the task-based method, which improved detectability up to ~15%. The improvement was strongest in areas of high attenuation (skull base), where the conventional and certainty-based methods tended to over-smooth the data. The task-driven reconstruction method presents a promising regularization method in MBIR by explicitly incorporating task-based imaging performance as the objective. The results demonstrate improved ICH conspicuity and support the development of high-quality CBCT systems.

Dual-rotation C-arm cone-beam computed tomography to increase low-contrast detection

PURPOSE

This paper investigates the capabilities of a dual-rotation C-arm cone-beam computed tomography (CBCT) framework to improve non-contrast-enhanced low-contrast detection for full volume or volume-of-interest (VOI) brain imaging.

METHOD

The idea is to associate two C-arm short-scan rotational acquisitions (spins): one over the full detector field of view (FOV) at low dose, and one collimated to deliver a higher dose to the central densest parts of the head. The angular sampling performed by each spin is allowed to vary in terms of number of views and angular positions. Collimated data is truncated and does not contain measurement of the incoming X-ray intensities in air (air calibration). When targeting full volume reconstruction, the method is intended to act as a virtual bow-tie. When targeting VOI imaging, the method is intended to provide the minimum full detector FOV data that sufficiently corrects for truncation artifacts. A single dedicated iterative algorithm is described that handles all proposed sampling configurations despite truncation and absence of air calibration.

RESULTS

Full volume reconstruction of dual-rotation simulations and phantom acquisitions are shown to have increased low-contrast detection for less dose, with respect to a single-rotation acquisition. High CNR values were obtained on 1% inserts of the Catphan® 515 module in 0.94 mm thick slices. Image quality for VOI imaging was preserved from truncation artifacts even with less than 10 non-truncated views, without using the sparsity a priori common to such context.

CONCLUSION

A flexible dual-rotation acquisition and reconstruction framework is proposed that has the potential to improve low-contrast detection in clinical C-arm brain soft-tissue imaging.

Technical assessment of a prototype cone-beam CT system for imaging of acute intracranial hemorrhage

PURPOSE

A cone-beam CT scanner has been developed for detection and monitoring of traumatic brain injury and acute intracranial hemorrhage (ICH) at the point of care. This work presents a technical assessment of imaging performance and dose for the scanner in phantom and cadaver studies as a prerequisite to clinical translation.

METHODS

The scanner incorporates a compact, rotating-anode x-ray source and a flat-panel detector (43 × 43 cm²) on a mobile U-arm gantry with source-axis distance = 550 mm and source-detector distance = 1000 mm. Central and peripheral doses were measured in 16 cm diameter CTDI phantoms using a 0.6 cm³ Farmer ionization chamber for various scan techniques and as a function of longitudinal position, including out of field. Spatial resolution, contrast, noise, and image uniformity were assessed in quantitative and anthropomorphic head phantoms. Two reconstruction protocols were evaluated, including filtered backprojection (FBP) for high-resolution bone imaging and penalized weighted least squares (PWLS) reconstruction for low-contrast soft tissue (ICH) visualization. A fresh cadaver was imaged with and without simulated ICH using the scanner as well as a diagnostic multidetector CT (MDCT) scanner using a standard head protocol. Images were interpreted by a fellowship-trained neuroradiologist for imaging tasks of ICH detection, gray-white-CSF differentiation, detection of midline shift, and fracture detection.

RESULTS

The nominal scan protocol involved 720 projections acquired over a 360° orbit at 100 kV and 216 mAs, giving a dose (weighted CTDI) of 22.8 mGy (∼1.2 mSv effective dose). Out-of-field dose decreased to <10% within 6 cm of the field edge (approximate to the thyroid position). Image uniformity demonstrated <1% variation between the edge of the field (near the cranium) and center of the image. The high-resolution FBP reconstruction protocol showed ∼0.9 mm point spread function (PSF) full-width at half-maximum (FWHM). The smooth PWLS reconstruction protocol yielded ∼1.2 mm PSF FWHM and contrast-to-noise ratio exceeding 5.7 in ∼50 HU spherical ICH, resulting in conspicuous depiction of ICH down to ∼2 mm (the smallest diameter investigated). Cadaver images demonstrated good differentiation of brain and CSF (sufficient, but inferior to MDCT, recognizing that the CBCT dose was one-third that of MDCT), excellent visualization of cranial sutures and fracture (potentially superior to MDCT), clear detection of midline shift, and conspicuous detection of ICH.

CONCLUSIONS

Technical assessment of the prototype demonstrates dose characteristics and imaging performance consistent with point-of-care detection and monitoring of head injury-most notably, conspicuous detection of ICH-and supports translation of the system to clinical studies.

Multi-resolution statistical image reconstruction for mitigation of truncation effects: application to cone-beam CT of the head

A prototype cone-beam CT (CBCT) head scanner featuring model-based iterative reconstruction (MBIR) has been recently developed and demonstrated the potential for reliable detection of acute intracranial hemorrhage (ICH), which is vital to diagnosis of traumatic brain injury and hemorrhagic stroke. However, data truncation (e.g. due to the head holder) can result in artifacts that reduce image uniformity and challenge ICH detection. We propose a multi-resolution MBIR method with an extended reconstruction field of view (RFOV) to mitigate truncation effects in CBCT of the head. The image volume includes a fine voxel size in the (inner) nontruncated region and a coarse voxel size in the (outer) truncated region. This multi-resolution scheme allows extension of the RFOV to mitigate truncation effects while introducing minimal increase in computational complexity. The multi-resolution method was incorporated in a penalized weighted least-squares (PWLS) reconstruction framework previously developed for CBCT of the head. Experiments involving an anthropomorphic head phantom with truncation due to a carbon-fiber holder were shown to result in severe artifacts in conventional single-resolution PWLS, whereas extending the RFOV within the multi-resolution framework strongly reduced truncation artifacts. For the same extended RFOV, the multi-resolution approach reduced computation time compared to the single-resolution approach (viz. time reduced by 40.7%, 83.0%, and over 95% for an image volume of 600³, 800³, 1000³ voxels). Algorithm parameters (e.g. regularization strength, the ratio of the fine and coarse voxel size, and RFOV size) were investigated to guide reliable parameter selection. The findings provide a promising method for truncation artifact reduction in CBCT and may be useful for other MBIR methods and applications for which truncation is a challenge.