Multi-Stage Adaptive Spline Autofocus (MASA) with a Learned Metric for Deformable Motion Compensation in Interventional Cone-Beam CT

Purpose

Cone-beam CT (CBCT) is widespread in abdominal interventional imaging, but its long acquisition time makes it susceptible to patient motion. Image-based autofocus has shown success in CBCT deformable motion compensation, via deep autofocus metrics and multi-region optimization, but it is challenged by the large parameter dimensionality required to capture intricate motion trajectories. This work leverages the differentiable nature of deep autofocus metrics to build a novel optimization strategy, Multi-Stage Adaptive Spine Autofocus (MASA), for compensation of complex deformable motion in abdominal CBCT.

Methods

MASA poses the autofocus problem as a multi-stage adaptive sampling strategy of the motion trajectory, sampled with Hermite spline basis with variable amplitude and knot temporal positioning. The adaptive method permits simultaneous optimization of the sampling phase, local temporal sampling density, and time-dependent amplitude of the motion trajectory. The optimization is performed in a multi-stage schedule with increasing number of knots that progressively accommodates complex trajectories in late stages, preconditioned by coarser components from early stages, and with minimal increase in dimensionality. MASA was evaluated in controlled simulation experiments with two types of motion trajectories: i) combinations of slow drifts with sudden jerk (sigmoid) motion; and ii) combinations of periodic motion sources of varying frequency into multi-frequency trajectories. Further validation was obtained in clinical data from liver CBCT featuring motion of contrast-enhanced vessels, and soft-tissue structures.

Results

The adaptive sampling strategy provided successful motion compensation in sigmoid trajectories, compared to fixed sampling strategies (mean SSIM increase of 0.026 compared to 0.011). Inspection of the estimated motion showed the capability of MASA to automatically allocate larger sampling density to parts of the scan timeline featuring sudden motion, effectively accommodating complex motion without increasing the problem dimension. Experiments on multi-frequency trajectories with 3-stage MASA (5, 10, and 15 knots) yielded a twofold SSIM increase compared to single-stage autofocus with 15 knots (0.076 vs 0.040, respectively). Application of MASA to clinical datasets resulted in simultaneous improvement on the delineation of both contrast-enhanced vessels and soft-tissue structures in the liver.

Conclusion

A new autofocus framework, MASA, was developed including a novel multi-stage technique for adaptive temporal sampling of the motion trajectory in combination with fully differentiable deep autofocus metrics. This novel adaptive sampling approach is a crucial step for application of deformable motion compensation to complex temporal motion trajectories.

A nonconvex model-based combined geometric calibration scheme for micro cone-beam CT with irregular trajectories

BACKGROUND

Many dedicated cone-beam CT (CBCT) systems have irregular scanning trajectories. Compared with the standard CBCT calibration, accurate calibration for CBCT systems with irregular trajectories is a more complex task, since the geometric parameters for each scanning view are variable. Most of the existing calibration methods assume that the intrinsic geometric relationship of the fiducials in the phantom is precisely known, and rarely delve deeper into the issue of whether the phantom accuracy is adapted to the calibration model.

PURPOSE

A high-precision phantom and a highly robust calibration model are interdependent and mutually supportive, and they are both important for calibration accuracy, especially for the high-resolution CBCT. Therefore, we propose a calibration scheme that considers both accurate phantom measurement and robust geometric calibration.

METHODS

Our proposed scheme consists of two parts: (1) introducing a measurement model to acquire the accurate intrinsic geometric relationship of the fiducials in the phantom; (2) developing a highly noise-robust nonconvex model-based calibration method. The measurement model in the first part is achieved by extending our previous high-precision geometric calibration model suitable for CBCT with circular trajectories. In the second part, a novel iterative method with optimization constraints based on a back-projection model is developed to solve the geometric parameters of each view.

RESULTS

The simulations and real experiments show that the measurement errors of the fiducial ball bearings (BBs) are within the subpixel level. With the help of the geometric relationship of the BBs obtained by our measurement method, the classic calibration method can achieve good calibration based on far fewer BBs. All metrics obtained in simulated experiments as well as in real experiments on Micro CT systems with resolutions of 9 and 4.5 μm show that the proposed calibration method has higher calibration accuracy than the competing classic method. It is particularly worth noting that although our measurement model proves to be very accurate, the classic calibration method based on this measurement model can only achieve good calibration results when the resolution of the measurement system is close to that of the system to be calibrated, but our calibration scheme enables high-accuracy calibration even when the resolution of the system to be calibrated is twice that of the measurement system.

CONCLUSIONS

The proposed combined geometrical calibration scheme does not rely on a phantom with an intricate pattern of fiducials, so it is applicable in Micro CT with high resolution. The two parts of the scheme, the "measurement model" and the "calibration model," prove to be of high accuracy. The combination of these two models can effectively improve the calibration accuracy, especially in some extreme cases.

3D-printed large-area focused grid for scatter reduction in cone-beam CT

BACKGROUND

Cone-beam computed tomography (CBCT) systems acquire volumetric data more efficiently than fan-beam or multislice CT, particularly when the anatomy of interest resides within the axial field-of-view of the detector and data can be acquired in one rotation. For such systems, scattered radiation remains a source of image quality degradation leading to increased noise, image artifacts, and CT number inaccuracies.

PURPOSE

Recent advances in metal additive manufacturing allow the production of highly focused antiscatter grids (2D-ASGs) that can be used to reduce scatter intensity, while preserving primary radiation transmission. We present the first implementation of a large-area, 2D-ASG for flat-panel CBCT, including grid-line artifact removal and related improvements in image quality.

METHODS

A 245 × 194 × 10 mm 2D-ASG was manufactured from chrome-cobalt alloy using laser powder-bed fusion (LPBF) (AM-400; Renishaw plc, New Mills Wotton-under-Edge, UK). The 2D-ASG had a square profile with a pitch of 9.09 lines/cm and 10:1 grid-ratio. The nominal 0.1 mm grid septa were focused to a 732 mm x-ray source to optimize primary x-ray transmission and reduce grid-line shadowing at the detector. Powder-bed fusion ensured the structural stability of the ASG with no need for additional interseptal support. The 2D-ASG was coupled to a 0.139-mm element pitch flat-panel detector (DRX 3543, Carestream Health) and proper alignment was confirmed by consistent grid-line shadow thickness across the whole detector array. A 154-mm diameter CBCT image-quality-assurance phantom was imaged using a rotary stage and a ceiling-mounted, x-ray unit (Proteus XR/a, GE Medical Systems, 80kVp, 0.5mAs). Grid-line artifacts were removed using a combination of exposure-dependent gain correction and spatial-frequency, Fourier filtering. Projections were reconstructed using a Parker-weighted, FDK algorithm and voxels were spatially averaged to 357 × 357 × 595 µm to improve the signal-to-noise characteristics of the CBCT reconstruction. Finally, in order to compare image quality with and without scatter, the phantom was scanned again under the same CBCT conditions but with no 2D-ASG. No additional antiscatter (i.e., air-gap, bowtie filtration) strategies were used to evaluate the effects in image quality caused by the 2D-ASG alone.

RESULTS

The large-area, 2D-ASG prototype was successfully designed and manufactured using LPBF. CBCT image-quality improvements using the 2D-ASG included: an overall 14.5% CNR increase across the volume; up to 48.8% CNR increase for low-contrast inserts inside the contrast plate of the QA phantom; and a 65% reduction of cupping artifact in axial profiles of water-filled cross sections of the phantom. Advanced image processing strategies to remove grid line artifacts did not affect the spatial resolution or geometric accuracy of the system.

CONCLUSIONS

LPBF can be used to manufacture highly efficient, 2D-focused ASGs that can be easily coupled to clinical, flat-panel detectors. The implementation of ASGs in CBCT leads to reduced scatter-related artifacts, improved CT number accuracy, and enhanced CNR with no increased equivalent dose to the patient. Further improvements to image quality might be achieved with a combination of scatter-correction algorithms and iterative-reconstruction strategies. Finally, clinical applications where other scatter removal strategies are unfeasible might now achieve superior soft-tissue visualization and quantitative capabilities.

CBCT-based Prostate IGRT With and Without Implanted Markers: Assessment of Geometric Corrections and Time for Completion

BACKGROUND/AIM

Cone-beam computed tomography (CBCT) is the most commonly used system in modern radiotherapy of prostate cancer for daily positioning verification. The use of intraprostatic radiopaque fiducials (FMs) may be added to CBCT. We wanted to investigate the possible advantage of using FMs in daily CBCT repositioning.

MATERIALS AND METHODS

We selected three CBCTs for each treatment course for 13 patients (seven with and six without use of FMs) treated at our centre. Seven experienced Radiation Oncologists retrospectively reviewed the CBCTs, recording couch movements for correct patient positioning, and time spent to do it. Analysis of variance and t-test were carried out for comparison of different groups and for differences in mean values of the movements recorded (with p<0.05 as significance level).

RESULTS

No statistically significant difference was found between operators in the analysis of images with FMs nor of images without them. A difference was only found in the mean corrections in couch rotation and pitch angle, which were higher in the FM group, and in the mean time for image analysis, which was shorter in this group. Using the van Herk formula, we found a possible reduction of clinical target volume and planning target volume margins for the FM group.

CONCLUSION

According to our study, the use of intraprostatic FMs in daily CBCT seems useful for better detection of and correction for non-negligible rotational errors. Furthermore, FMs reduced the time to treatment start, which is very important in reducing the risk of intrafraction organ motion. These results need to be confirmed by further studies.

Machine Learning Analysis of the Anatomical Parameters of the Upper Airway Morphology: A Retrospective Study from Cone-Beam CT Examinations in a French Population

The objective of this study is to assess, using cone-beam CT (CBCT) examinations, the correlation between hard and soft anatomical parameters and their impact on the characteristics of the upper airway using symbolic regression as a machine learning strategy. Methods: On each CBCT, the upper airway was segmented, and 24 anatomical landmarks were positioned to obtain six angles and 19 distances. Some anatomical landmarks were related to soft tissues and others were related to hard tissues. To explore which variables were the most influential to explain the morphology of the upper airway, principal component and symbolic regression analyses were conducted. Results: In total, 60 CBCT were analyzed from subjects with a mean age of 39.5 ± 13.5 years. The intra-observer reproducibility for each variable was between good and excellent. The horizontal soft palate measure mostly contributed to the reduction of the airway volume and minimal section area with a variable importance of around 50%. The tongue and the position of the hyoid bone were also linked to the upper airway morphology. For hard anatomical structures, the anteroposterior position of the mandible and the maxilla had some influence. Conclusions: Although the volume of the airway is not accessible on all CBCT scans performed by dental practitioners, this study demonstrates that a small number of anatomical elements may be markers of the reduction of the upper airway with, potentially, an increased risk of obstructive sleep apnea. This could help the dentist refer the patient to a suitable physician.

3D-2D image registration in the presence of soft-tissue deformation in image-guided transbronchial interventions

Purpose. Target localization in pulmonary interventions (e.g. transbronchial biopsy of a lung nodule) is challenged by deformable motion and may benefit from fluoroscopic overlay of the target to provide accurate guidance. We present and evaluate a 3D-2D image registration method for fluoroscopic overlay in the presence of tissue deformation using a multi-resolution/multi-scale (MRMS) framework with an objective function that drives registration primarily by soft-tissue image gradients.Methods. The MRMS method registers 3D cone-beam CT to 2D fluoroscopy without gating of respiratory phase by coarse-to-fine resampling and global-to-local rescaling about target regions-of-interest. A variation of the gradient orientation (GO) similarity metric (denotedGO') was developed to downweight bone gradients and drive registration via soft-tissue gradients. Performance was evaluated in terms of projection distance error at isocenter (PDE_iso). Phantom studies determined nominal algorithm parameters and capture range. Preclinical studies used a freshly deceased, ventilated porcine specimen to evaluate performance in the presence of real tissue deformation and a broad range of 3D-2D image mismatch.Results. Nominal algorithm parameters were identified that provided robust performance over a broad range of motion (0-20 mm), including an adaptive parameter selection technique to accommodate unknown mismatch in respiratory phase. TheGO'metric yielded median PDE_iso= 1.2 mm, compared to 6.2 mm for conventionalGO.Preclinical studies with real lung deformation demonstrated median PDE_iso= 1.3 mm with MRMS +GO'registration, compared to 2.2 mm with a conventional transform. Runtime was 26 s and can be reduced to 2.5 s given a prior registration within ∼5 mm as initialization.Conclusions. MRMS registration via soft-tissue gradients achieved accurate fluoroscopic overlay in the presence of deformable lung motion. By driving registration via soft-tissue image gradients, the method avoided false local minima presented by bones and was robust to a wide range of motion magnitude.

Image quality models for 2D and 3D x-ray imaging systems: A perspective vignette

Image quality models based on cascaded systems analysis and task-based imaging performance were an important aspect of the emergence of 2D and 3D digital x-ray systems over the last 25 years. This perspective vignette offers cursory review of such developments and personal insights that may not be obvious within previously published scientific literature. The vignette traces such models to the mid-1990s, when flat-panel x-ray detectors were emerging as a new base technology for digital radiography and benefited from the rigorous, objective characterization of imaging performance gained from such models. The connection of models for spatial resolution and noise to spatial-frequency-dependent descriptors of imaging task provided a useful framework for system optimization that helped to accelerate the development of new technologies to first clinical use. Extension of the models to new technologies and applications is also described, including dual-energy imaging, photon-counting detectors, phase contrast imaging, tomosynthesis, cone-beam CT, 3D image reconstruction, and image registration.

[Automatic Delineation of Clinical Target Volume and Organ at Risk by Deep Learning for Prostate Cancer Adaptive Radiotherapy]

Adaptive radiotherapy can modify the treatment plan online based on the clinical target volume (CTV) and organ at risk (OAR) contours on the cone-beam CT (CBCT) before treatment, improving the accuracy of radiotherapy. However, manual delineation of CTV and OAR on CBCT is time-consuming. In this study, a deep neural network-based method based on U-Net was purposed. CBCT images and corresponding mask were used for model training and validation, showing superior performance in terms of the segmentation accuracy. The proposed method could be used in the clinic to support rapid CTV and OAR contouring for prostate adaptive radiotherapy.

Effect of subject motion and gantry rotation speed on image quality and dose delivery in CT-guided radiotherapy

PURPOSE

To investigate the effects of subject motion and gantry rotation speed on computed tomography (CT) image quality over a range of image acquisition speeds for fan-beam (FB) and cone-beam (CB) CT scanners, and quantify the geometric and dosimetric errors introduced by FB and CB sampling in the context of adaptive radiotherapy.

METHODS

Images of motion phantoms were acquired using four CT scanners with gantry rotation speeds of 0.5 s/rotation (denoted FB-0.5), 1.9 s/rotation (FB-1.9), 16.6 s/rotation (CB-16.6), and 60.0 s/rotation (CB-60.0). A phantom presenting various tissue densities undergoing motion with 4-s period and ranging in amplitude from ±0.5 to ±10.0 mm was used to characterize motion artifacts (streaks), motion blur (edge-spread function, ESF), and geometric inaccuracy (excursion of insert centroids and distortion of known shape). An anthropomorphic abdomen phantom undergoing ±2.5-mm motion with 4-s period was used to simulate an adaptive radiotherapy workflow, and relative geometric and dosimetric errors were compared between scanners.

RESULTS

At ±2.5-mm motion, phantom measurements demonstrated mean ± SD ESF widths of 0.6 ± 0.0, 1.3 ± 0.4, 2.0 ± 1.1, and 2.9 ± 2.0 mm and geometric inaccuracy (excursion) of 2.7 ± 0.4, 4.1 ± 1.2, 2.6 ± 0.7, and 2.0 ± 0.5 mm for the FB-0.5, FB-1.9, CB-16.6, and CB-60.0 scanners, respectively. The results demonstrated nonmonotonic trends with scanner speed for FB and CB geometries. Geometric and dosimetric errors in adaptive radiotherapy plans were largest for the slowest (CB-60.0) scanner and similar for the three faster systems (CB-16.6, FB-1.9, and FB-0.5).

CONCLUSIONS

Clinically standard CB-60.0 demonstrates strong image quality degradation in the presence of subject motion, which is mitigated through faster CBCT or FBCT. Although motion blur is minimized for FB-0.5 and FB-1.9, such systems suffer from increased geometric distortion compared to CB-16.6. Each system reflects tradeoffs in image artifacts and geometric inaccuracies that affect treatment delivery/dosimetric error and should be considered in the design of next-generation CT-guided radiotherapy systems.

Evaluation of intrafraction couch shifts for proton treatment delivery in head-and-neck cancer patients: Toward optimal imaging frequency

PURPOSE

Treatment planning for head-and-neck (H&N) cancer, in particular oropharynx, nasopharynx, and paranasal sinus cases, at our center requires noncoplanar proton beams due to the complexity of the anatomy and target location. Targeting accuracy for all beams is carefully evaluated by using image guidance before delivering proton beam therapy (PBT). In this study, we analyzed couch shifts to evaluate whether imaging is required before delivering each field with different couch angles.

METHODS

After the Institutional Review Board approval, a retrospective analysis was performed on data from 28 H&N patients treated with PBT. Each plan was made with two-to-three noncoplanar and two-to-three coplanar fields. Cone-beam computed tomography and orthogonal kilovoltage (kV) images were acquired for setup and before delivering each field, respectively. The Cartesian (longitudinal, vertical, and lateral) and angular (pitch and roll) shifts for each field were recorded from the treatment summary on the first two fractions and every subsequent fifth fraction. A net magnitude of the three-dimensional (3D) shift in Cartesian coordinates was calculated, and a 3D vector was created from the 6 degrees of freedom coordinates for transforming couch shifts in the system coordinate to the beam's-eye view.

RESULTS

A total of 3219 Cartesian and 2146 angular shift values were recorded for 28 patients. Of the Cartesian shifts, 2069 were zero (64.3%), and 1150 (35.7%) were nonzero (range, -7 to 11 mm). Of the angular shifts, 1034 (48.2%) were zero, and 1112 (51.8%) were nonzero (range, -3.0° to 3.2°). For 17 patients, the couch shifts increased toward the end of the treatment course. We also found that patients with higher body mass index (BMI) presented increased net couch shifts (p < 0.001). With BMI < 27, all overall net shift averages were <2 mm, and overall maximum net shifts were <6 mm.

CONCLUSIONS

These results confirm the need for orthogonal kV imaging before delivering each field of H&N PBT at our center, where a couch rotation is involved.

Dosimetric and feasibility evaluation of a CBCT-based daily adaptive radiotherapy protocol for locally advanced cervical cancer

PURPOSE

Evaluate a cone-beam computed tomography (CBCT)-based daily adaptive platform in cervical cancer for multiple endpoints: (1) physics contouring accuracy of daily CTVs, (2) CTV coverage with adapted plans and reduced PTV margins versus non-adapted plans with standard-of-care (SOC) margins, (3) dosimetric improvements to CTV and organs-at-risk (OARs), and (4) on-couch time.

METHODS AND MATERIALS

Using a Varian Ethos™ emulator and KV-CBCT scans, we simulated the doses 15 retrospective cervical cancer patients would have received with/without online adaptation for five fractions. We compared contours and doses from SOC plans (5-15 mm CTV-to-PTV margins) to adapted plans (3 mm margins). Auto-segmented CTVs and OARs were reviewed and edited by trained physicists. Physics-edited targets were evaluated by an oncologist. Time spent reviewing and editing auto-segmented structures was recorded. Metrics from the CTV (D99%), bowel (V45Gy, V40Gy), bladder (D50%), and rectum (D50%) were compared.

RESULTS

The physician approved the physics-edited CTVs for 55/75 fractions; 16/75 required reductions, and 4/75 required CTV expansions. CTVs were encapsulated by unadapted, SOC PTVs for 56/75 (72%) fractions-representative of current clinical practice. CTVs were completely covered by adapted 3 mm PTVs for 71/75 (94.6%) fractions. CTV D99% values for adapted plans were comparable to non-adapted SOC plans (average difference of -0.9%), while all OAR metrics improved with adaptation. Specifically, bowel V45Gy and V40Gy decreased on average by 87.6 and 109.4 cc, while bladder and rectum D50% decreased by 37.7% and 35.8%, respectively. The time required for contouring and calculating an adaptive plan for 65/75 fractions was less than 20 min (range: 1-29 min).

CONCLUSIONS

Improved dose metrics with daily adaption could translate to reduced toxicity while maintaining tumor control. Training physicists to perform contouring edits could minimize the time physicians are required at adaptive sessions improving clinical efficiency. All emulated adaptive sessions were completed within 30 min however extra time will be required for patient setup, image acquisition, and treatment delivery.

Cone-beam computed tomography guided nusinersen administrations in adult spinal muscular atrophy patients with challenging access: a single- center experience

BACKGROUND

The challenging anatomic predispositions in adult patients with spinal muscular atrophy (SMA) preclude the conventional lumbar punctures. Consequently, an introduction of alternative method for intrathecal delivery of nusinersen is required. Cone-beam CT (CBCT) allows volumetric display of the area of interest, pre-procedural planning and real time needle guidance which results in accurate anatomic navigation. The aim of the study was to evaluate technical success, safety, and feasibility of CBCT lumbar intrathecal delivery of nusinersen in the adult SMA patients with challenging anatomical access.

PATIENTS AND METHODS

Thirty-eight adult SMA patients were treated in our institution. Patients with challenging access were selected by multidisciplinary board for image guided administration of nusinersen either due to implantation of the posterior fusion instrumentation, severe scoliosis defined as Cobb's angle > 40º or body mass index over 35. Technical success, radiation exposure and occurrence of adverse events were assessed.

RESULTS

Twenty patients were selected, and 108 CBCT-guided procedures were performed. Each patient underwent at least 4 administrations. Transforaminal approach was performed in 82% of patients. The technical success was 100%, with primary success of 93.5%. The median radiation effective dose of the administrations was 5 mSv, the mean value equalled 10 mSv. Only mild adverse events were reported in the study.

CONCLUSIONS

CBCT-guided lumbar intrathecal administrations of nusinersen in an adult SMA population with challenging access was feasible and safe image guided method.

Cone Beam CT to Guide Transorbital Treatment of a Cavernous Sinus Dural Arteriovenous Fistula in a Patient with Middle Meningeal Artery Origin of the Ophthalmic Artery

Transorbital puncture to embolize cavernous sinus (CS) dural arteriovenous fistulas (DAVF) is a useful strategy when conventional transvenous routes are inaccessible. We report a case of bilateral CS DAVF associated with bilateral middle meningeal artery (MMA) origin of ophthalmic arteries (OA) who had recently undergone transvenous coil embolization. She presented with persistent symptoms of conjunctival congestion and proptosis in the left eye. Angiogram revealed residual left CS DAVF with dilated SOV. Inferior petrosal sinus or facial vein access was not possible. Transorbital access of the SOV was planned. Cone-beam CT (CBCT) angiography was used to delineate the relationship between the variant OA and SOV and also to plan a safe trajectory. Using fluoroscopy guidance, the SOV was punctured and embolization was done using Onyx-18. CBCT is a valuable tool in planning and executing transorbital treatment of CS DAVF, especially in the setting of variant OA.

Reference-free learning-based similarity metric for motion compensation in cone-beam CT

Purpose. Patient motion artifacts present a prevalent challenge to image quality in interventional cone-beam CT (CBCT). We propose a novel reference-free similarity metric (DL-VIF) that leverages the capability of deep convolutional neural networks (CNN) to learn features associated with motion artifacts within realistic anatomical features. DL-VIF aims to address shortcomings of conventional metrics of motion-induced image quality degradation that favor characteristics associated with motion-free images, such as sharpness or piecewise constancy, but lack any awareness of the underlying anatomy, potentially promoting images depicting unrealistic image content. DL-VIF was integrated in an autofocus motion compensation framework to test its performance for motion estimation in interventional CBCT.Methods. DL-VIF is a reference-free surrogate for the previously reported visual image fidelity (VIF) metric, computed against a motion-free reference, generated using a CNN trained using simulated motion-corrupted and motion-free CBCT data. Relatively shallow (2-ResBlock) and deep (3-Resblock) CNN architectures were trained and tested to assess sensitivity to motion artifacts and generalizability to unseen anatomy and motion patterns. DL-VIF was integrated into an autofocus framework for rigid motion compensation in head/brain CBCT and assessed in simulation and cadaver studies in comparison to a conventional gradient entropy metric.Results. The 2-ResBlock architecture better reflected motion severity and extrapolated to unseen data, whereas 3-ResBlock was found more susceptible to overfitting, limiting its generalizability to unseen scenarios. DL-VIF outperformed gradient entropy in simulation studies yielding average multi-resolution structural similarity index (SSIM) improvement over uncompensated image of 0.068 and 0.034, respectively, referenced to motion-free images. DL-VIF was also more robust in motion compensation, evidenced by reduced variance in SSIM for various motion patterns (σ_DL-VIF = 0.008 versusσ_{gradient entropy} = 0.019). Similarly, in cadaver studies, DL-VIF demonstrated superior motion compensation compared to gradient entropy (an average SSIM improvement of 0.043 (5%) versus little improvement and even degradation in SSIM, respectively) and visually improved image quality even in severely motion-corrupted images.Conclusion: The studies demonstrated the feasibility of building reference-free similarity metrics for quantification of motion-induced image quality degradation and distortion of anatomical structures in CBCT. DL-VIF provides a reliable surrogate for motion severity, penalizes unrealistic distortions, and presents a valuable new objective function for autofocus motion compensation in CBCT.

Accuracy of Mini-Implant Placement Using a Computer-Aided Designed Surgical Guide, with Information of Intraoral Scan and the Use of a Cone-Beam CT

Background: The purpose of the study was to investigate the accuracy of mini-implant placement with the use of a computer designed surgical guide derived by intraoral scanning alongside Cone-Beam Computed Tomography (CBCT) or the use of a 2D radiograph. Methods: Thirty-five mini-implants (Aarhus^® System: n = 20, Ø: 1.5 mm and AbsoAnchor^®: n = 15, Ø: 1.3–1.4 mm) were placed in the maxilla and mandible of 15 orthodontic patients for anchorage purposes in cases where a CBCT was needed due to impacted teeth or for safety reasons. All were placed with the help of a computer designed surgical guide. One implant became loose and was excluded from the study. For 18 mini-implants (study group), CBCT and intraoral scanning were used for guide design, while for 16 (control group) only intraoral scanning and panoramic imaging information were used. Mini-implant position was recorded by angular and linear measurements on digital models created by combining Digital Imaging and Communications in Medicine (DICOM) and stereolithography (.stl) files. Accuracy in positioning was determined by comparing corresponding measurements for virtual and real positioned mini-implants on digital models before and after operation. The results were statistically analyzed with t-tests and the Mann-Whitney test. Results: No significant statistical differences were found for pre- and post-operational angular measurements in the study group, while significant statistical differences occurred on the same measurements for the control group (coronal angle 13.6° pre-op and 22.7° post-op, p-value = 0.002, axial angle 13.4° pre-op and 15.9° post-op, p-value = 0.034). Linear measurements pre- and post-operational for either group presented no significant statistical differences. Conclusions: A 3D designed and manufactured surgical guide with information concerning CBCT and intraoral scanning ensures accuracy on mini-implant placement while design of the guide without the use of a CBCT is less accurate, especially on inclination of the implant.

Dual-Energy Cone-Beam CT with Three-Material Decomposition for Bone Marrow Edema Imaging

We investigate the feasibility of bone marrow edema (BME) detection using a kV-switching Dual-Energy (DE) Cone-Beam CT (CBCT) protocol. This task is challenging due to unmatched x-ray paths in the low-energy (LE) and high-energy (HE) spectral channels, CBCT non-idealities such as x-ray scatter, and narrow spectral separation between fat (bone marrow) and water (BME). We propose a comprehensive DE decomposition framework consisting of projection interpolation onto matching LE and HE view angles, fast Monte Carlo scatter correction with low number of tracked photons and Gaussian denoising, and two-stage three-material decompositions involving two-material (fat-Aluminium) Projection-Domain Decomposition (PDD) followed by image-domain three-material (fat-water-bone) base-change. Performance in BME detection was evaluated in simulations and experiments emulating a kV-switching CBCT wrist imaging protocol on a robotic x-ray system with 60 kV LE beam, 120 kV HE beam, and 0.5° angular shift between the LE and HE views. Cubic B-spline interpolation was found to be adequate to resample HE and LE projections of a wrist onto common view angles required by PDD. The DE decomposition maintained acceptable BME detection specificity (<0.2 mL erroneously detected BME volume compared to 0.85 mL true BME volume) over +/-10% range of scatter magnitude errors, as long as the scatter shape was estimated without major distortions. Physical test bench experiments demonstrated successful discrimination of ~20% change in fat concentrations in trabecular bone-mimicking solutions of varying water and fat content.

Diagnostic Performance of Ultrasonography for Evaluation of Osteoarthritis of Ankle Joint: Comparison With Radiography, Cone-Beam CT, and Symptoms

OBJECTIVES

To determine the diagnostic performance of ultrasonography (US) for evaluation of the ankle joint osteoarthritic (OA) changes. Cone-beam computed tomography (CT) was used as the gold standard and US performance was compared with conventional radiography (CR). As a secondary aim, associations between the imaging findings and ankle symptoms were assessed.

METHODS

US was performed to 51 patients with ankle OA. Every patient had prior ankle CR and underwent cone-beam CT during the same day as US examination. On US, effusion/synovitis, osteophytes, talar cartilage damage, and tenosynovitis were evaluated. Comparison to respective imaging findings on CR and cone-beam CT was then performed. Single radiologist blinded to other modalities assessed all the imaging studies. Symptoms questionnaire, the Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC), was available for 48 patients.

RESULTS

US detected effusion/synovitis of the talocrural joint with 45% sensitivity and 90% specificity. For the detection of anterior talocrural osteophytes, US sensitivity was 78% and specificity 79%. For the medial talocrural osteophytes, they were 39 and 83%, and for the lateral talocrural osteophytes 54 and 100%, respectively. Considering cartilage damage of the talus, US yielded a low sensitivity of 18% and high specificity of 97%. Overall, the performance of US was only moderate and comparable to CR. The imaging findings showed only weak associations with ankle symptoms.

CONCLUSIONS

The ability of US to detect ankle OA is only moderate. Interestingly, performance of CR also remained moderate. The associations between imaging findings and WOMAC score seem to be weak in ankle OA.

Dosimetric evaluation of cone-beam CT-based synthetic CTs in pediatric patients undergoing intensity-modulated proton therapy

PURPOSE

To evaluate dosimetric changes detected using synthetic computed tomography (sCT) derived from online cone-beam CTs (CBCT) in pediatric patients treated using intensity-modulated proton therapy (IMPT).

METHODS

Ten pediatric patients undergoing IMPT and aligned daily using proton gantry-mounted CBCT were identified for retrospective analysis with treated anatomical sites fully encompassed in the CBCT field of view. Dates were identified when the patient received both a CBCT and a quality assurance CT (qCT) for routine dosimetric evaluation. sCTs were generated based on a deformable registration between the initial plan CT (pCT) and CBCT. The clinical IMPT plans were re-computed on the same day qCT and sCT, and dosimetric changes due to tissue change or response from the initial plan were computed using each image. Linear regression analysis was performed to determine the correlation between dosimetric changes detected using the qCT and the sCT. Gamma analysis was also used to compare the dose distributions computed on the qCT and sCT.

RESULTS

The correlation coefficients (p-values) between qCTs and sCTs for changes detected in target coverage, overall maximum dose, and organ at risk dose were 0.97 (< .001), 0.84 (.002) and 0.91 (< .001), respectively. Mean ± SD gamma pass rates of the sCT-based dose compared to the qCT-based dose at 3%/3 mm, 3%/2 mm, and 2%/2 mm criteria were 96.5%±4.5%, 93.2%±6.3%, and 91.3%±7.8%, respectively. Pass rates tended to be lower for targets near lung.

CONCLUSION

While insufficient for re-planning, sCTs provide approximate dosimetry without administering additional imaging dose in pediatric patients undergoing IMPT. Dosimetric changes detected using sCTs are correlated with changes detected using clinically-standard qCTs; however, residual differences in dosimetry remain a limitation. Further improvements in sCT image quality may both improve online dosimetric evaluation and reduce imaging dose for pediatric patients by reducing the need for routine qCTs.

Cone-beam breast CT using an offset detector: effect of detector offset and image reconstruction algorithm

Objective.A dedicated cone-beam breast computed tomography (BCT) using a high-resolution, low-noise detector operating in offset-detector geometry has been developed. This study investigates the effects of varying detector offsets and image reconstruction algorithms to determine the appropriate combination of detector offset and reconstruction algorithm.Approach.Projection datasets (300 projections in 360°) of 30 breasts containing calcified lesions that were acquired using a prototype cone-beam BCT system comprising a 40 × 30 cm flat-panel detector with 1024 × 768 detector pixels were reconstructed using Feldkamp-Davis-Kress (FDK) algorithm and served as the reference. The projection datasets were retrospectively truncated to emulate cone-beam datasets with sinograms of 768×768 and 640×768 detector pixels, corresponding to 5 cm and 7.5 cm lateral offsets, respectively. These datasets were reconstructed using the FDK algorithm with appropriate weights and an ASD-POCS-based Fast, total variation-Regularized, Iterative, Statistical reconstruction Technique (FRIST), resulting in a total of 4 offset-detector reconstructions (2 detector offsets × 2 reconstruction methods). Signal difference-to-noise ratio (SDNR), variance, and full-width at half-maximum (FWHM) of calcifications in two orthogonal directions were determined from all reconstructions. All quantitative measurements were performed on images in units of linear attenuation coefficient (1/cm).Results.The FWHM of calcifications did not differ (P > 0.262) among reconstruction algorithms and detector formats, implying comparable spatial resolution. For a chosen detector offset, the FRIST algorithm outperformed FDK in terms of variance and SDNR (P < 0.0001). For a given reconstruction method, the 5 cm offset provided better results.Significance.This study indicates the feasibility of using the compressed sensing-based, FRIST algorithm to reconstruct sinograms from offset-detectors. Among the reconstruction methods and detector offsets studied, FRIST reconstructions corresponding to a 30 cm × 30 cm with 5 cm lateral offset, achieved the best performance. A clinical prototype using such an offset geometry has been developed and installed for clinical trials.

Technical assessment of 2D and 3D imaging performance of an IGZO-based flat-panel X-ray detector

BACKGROUND

Indirect detection flat-panel detectors (FPDs) consisting of hydrogenated amorphous silicon (a-Si:H) thin-film transistors (TFTs) are a prevalent technology for digital x-ray imaging. However, their performance is challenged in applications requiring low exposure levels, high spatial resolution, and high frame rate. Emerging FPD designs using metal oxide TFTs may offer potential performance improvements compared to FPDs based on a-Si:H TFTs.

PURPOSE

This work investigates the imaging performance of a new indium gallium zinc oxide (IGZO) TFT-based detector in 2D fluoroscopy and 3D cone-beam CT (CBCT).

METHODS

The new FPD consists of a sensor array combining IGZO TFTs with a-Si:H photodiodes and a 0.7-mm thick CsI:Tl scintillator. The FPD was implemented on an x-ray imaging bench with system geometry emulating intraoperative CBCT. A conventional FPD with a-Si:H TFTs and a 0.6-mm thick CsI:Tl scintillator was similarly implemented as a basis of comparison. 2D imaging performance was characterized in terms of electronic noise, sensitivity, linearity, lag, spatial resolution (modulation transfer function, MTF), image noise (noise-power spectrum, NPS), and detective quantum efficiency (DQE) with entrance air kerma (EAK) ranging from 0.3 to 1.2 μGy. 3D imaging performance was evaluated in terms of the 3D MTF and noise-equivalent quanta (NEQ), soft-tissue contrast-to-noise ratio (CNR), and image quality evident in anthropomorphic phantoms for a range of anatomical sites and dose, with weighted air kerma, K w ${K_w}$ , ranging from 0.8 to 4.9 mGy.

RESULTS

The 2D imaging performance of the IGZO-based FPD exhibited up to ∼1.7× lower electronic noise than the a-Si:H FPD at matched pixel pitch. Furthermore, the IGZO FPD exhibited ∼27% increase in mid-frequency DQE (1 mm^-1 ) at matched pixel size and dose (EAK ≈ 1.0 μGy) and ∼11% increase after adjusting for differences in scintillator thickness. 2D spatial resolution was limited by the scintillator for each FPD. The IGZO-based FPD demonstrated improved 3D NEQ at all spatial frequencies in both head (≥25% increase for all dose levels) and body (≥10% increase for K w ${K_w}$ ≤2 mGy) imaging scenarios. These characteristics translated to improved low-contrast visualization in anthropomorphic phantoms, demonstrating ≥10% improvement in CNR and extension of the low-dose range for which the detector is input-quantum limited.

CONCLUSION

The IGZO-based FPD demonstrated improvements in electronic noise, image lag, and NEQ that translated to measurable improvements in 2D and 3D imaging performance compared to a conventional FPD based on a-Si:H TFTs. The improvements are most beneficial for 2D or 3D imaging scenarios involving low-dose and/or high-frame rate.

Head and neck synthetic CT generated from ultra-low-dose cone-beam CT following Image Gently Protocol using deep neural network

PURPOSE

Image guidance is used to improve accuracy of radiation therapy delivery but results in increased dose to patients. This is of particular concern in children who need be treated per Pediatric Image Gently Protocols due to long term risks from radiation exposure. The purpose of this study is to design a deep neural network (DNN) architecture and loss function for improving soft-tissue contrast and preserving small anatomical features in ultra-low-dose cone-beam CTs (CBCT) of head and neck cancer (HNC) imaging.

METHODS

A 2-D compound U-Net architecture (modified U-Net++) with different depths was proposed to enhance the network capability of capturing small-volume structures. A mask weighted loss function (Mask-Loss) was applied to enhance soft-tissue contrast. Fifty-five paired CBCT and CT images of HNC patients were retrospectively collected for network training and testing. The output enhanced CBCT images from the present study were evaluated with quantitative metrics including mean absolute error (MAE), signal-to-noise ratio (SNR), and structural similarity (SSIM), and compared with those from the previously proposed network architectures (U-Net and wide U-Net) using MAE loss functions. A visual assessment of ten selected structures in the enhanced CBCT images of each patient was performed to evaluate image quality improvement, blindly scored by an experienced radiation oncologist specialized in HN cancer.

RESULTS

All the enhanced CBCT images showed reduced artifactual distortion and image noise. U-Net++ outperformed the U-Net and wide U-Net in terms of MAE, contrast near structure boundaries, and small structures. The proposed Mask-Loss improved image contrast and accuracy of the soft-tissue regions. The enhanced CBCT images predicted by U-Net++ and Mask-Loss demonstrated improvement compared to the U-Net in terms of average MAE (52.41 vs. 42.85 HU), SNR (14.14 vs. 15.07 dB), and SSIM (0.84 vs. 0.87), respectively (p < 0.01, in all paired t-tests). The visual assessment showed that the proposed U-Net++ and Mask-Loss significantly improved original CBCTs (p < 0.01), compared to the U-Net and MAE loss.

CONCLUSIONS

The proposed network architecture and loss function effectively improved image quality in soft-tissue contrast, organ boundary, and small structures preservation for ultra-low-dose CBCT following Image Gently Protocol. This method has potential to provide sufficient anatomical representation on the enhanced CBCT images for accurate treatment delivery and potentially fast online-adaptive re-planning for HN cancer patients. This article is protected by copyright. All rights reserved.

Improved Diagnostic Accuracy of Ameloblastoma and Odontogenic Keratocyst on Cone-Beam CT by Artificial Intelligence

OBJECTIVE

The purpose of this study was to utilize a convolutional neural network (CNN) to make preoperative differential diagnoses between ameloblastoma (AME) and odontogenic keratocyst (OKC) on cone-beam CT (CBCT).

METHODS

The CBCT images of 178 AMEs and 172 OKCs were retrospectively retrieved from the Hospital of Stomatology, Wuhan University. The datasets were randomly split into a training dataset of 272 cases and a testing dataset of 78 cases. Slices comprising lesions were retained and then cropped to suitable patches for training. The Inception v3 deep learning algorithm was utilized, and its diagnostic performance was compared with that of oral and maxillofacial surgeons.

RESULTS

The sensitivity, specificity, accuracy, and F1 score were 87.2%, 82.1%, 84.6%, and 85.0%, respectively. Furthermore, the average scores of the same indexes for 7 senior oral and maxillofacial surgeons were 60.0%, 71.4%, 65.7%, and 63.6%, respectively, and those of 30 junior oral and maxillofacial surgeons were 63.9%, 53.2%, 58.5%, and 60.7%, respectively.

CONCLUSION

The deep learning model was able to differentiate these two lesions with better diagnostic accuracy than clinical surgeons. The results indicate that the CNN may provide assistance for clinical diagnosis, especially for inexperienced surgeons.

Feasibility of Dual-Energy Cone-Beam CT of Bone Marrow Edema Using Dual-Layer Flat Panel Detectors

Purpose

We investigated the feasibility of detection and quantification of bone marrow edema (BME) using dual-energy (DE) Cone-Beam CT (CBCT) with a dual-layer flat panel detector (FPD) and three-material decomposition.

Methods

A realistic CBCT system simulator was applied to study the impact of detector quantization, scatter, and spectral calibration errors on the accuracy of fat-water-bone decompositions of dual-layer projections. The CBCT system featured 975 mm source-axis distance, 1,362 mm source-detector distance and a 430 × 430 mm2 dual-layer FPD (top layer: 0.20 mm CsI:Tl, bottom layer: 0.55 mm CsI:Tl; a 1 mm Cu filter between the layers to improve spectral separation). Tube settings were 120 kV (+2 mm Al, +0.2 mm Cu) and 10 mAs per exposure. The digital phantom consisted of a 160 mm water cylinder with inserts containing mixtures of water (volume fraction ranging 0.18 to 0.46) - fat (0.5 to 0.7) - Ca (0.04 to 0.12); decreasing fractions of fat indicated increasing degrees of BME. A two-stage three-material DE decomposition was applied to DE CBCT projections: first, projection-domain decomposition (PDD) into fat-aluminum basis, followed by CBCT reconstruction of intermediate base images, followed by image-domain change of basis into fat, water and bone. Sensitivity to scatter was evaluated by i) adjusting source collimation (12 to 400 mm width) and ii) subtracting various fractions of the true scatter from the projections at 400 mm collimation. The impact of spectral calibration was studied by shifting the effective beam energy (± 2 keV) when creating the PDD lookup table. We further simulated a realistic BME imaging framework, where the scatter was estimated using a fast Monte Carlo (MC) simulation from a preliminary decomposition of the object; the object was a realistic wrist phantom with an 0.85 mL BME stimulus in the radius.

Results

The decomposition is sensitive to scatter: approx. <20 mm collimation width or <10% error of scatter correction in a full field-of-view setting is needed to resolve BME. A mismatch in PDD decomposition calibration of ± 1 keV results in ~25% error in fat fraction estimates. In the wrist phantom study with MC scatter corrections, we were able to achieve ~0.79 mL true positive and ~0.06 mL false positive BME detection (compared to 0.85 mL true BME volume).

Conclusions

Detection of BME using DE CBCT with dual-layer FPD is feasible, but requires scatter mitigation, accurate scatter estimation, and robust spectral calibration.

Targeted Deformable Motion Compensation for Vascular Interventional Cone-Beam CT Imaging

Purpose

Cone-beam CT has become commonplace for 3D guidance in interventional radiology (IR), especially for vascular procedures in which identification of small vascular structures is crucial. However, its long image acquisition time poses a limit to image quality due to soft-tissue deformable motion that hampers visibility of small vessels. Autofocus motion compensation has shown promising potential for soft-tissue deformable motion compensation, but it lacks specific target to the imaging task. This work presents an approach for deformable motion compensation targeted at imaging of vascular structures.

Methods

The proposed method consists on a two-stage framework for: i) identification of contrast-enhanced blood vessels in 2D projection data and delineation of an approximate region covering the vascular target in the volume space, and, ii) a novel autofocus approach including a metric designed to promote the presence of vascular structures acting solely in the region of interest. The vesselness of the image is quantified via evaluation of the properties of the 3D image Hessian, yielding a vesselness filter that gives larger values to voxels candidate to be part of a tubular structure. A cost metric is designed to promote large vesselness values and spatial sparsity, as expected in regions of fine vascularity. A targeted autofocus method was designed by combining the presented metric with a conventional autofocus term acting outside of the region of interest. The resulting method was evaluated on simulated data including synthetic vascularity merged with real anatomical features obtained from MDCT data. Further evaluation was obtained in two clinical datasets obtained during TACE procedures with a robotic C-arm (Artis Zeego, Siemens Healthineers).

Results

The targeted vascular autofocus effectively restored the shape and contrast of the contrast-enhanced vascularity in the simulation cases, resulting in improved visibility and reduced artifacts. Segmentations performed with a single threshold value on the target vascular regions yielded a net increase of up to 42% in DICE coefficient computed against the static reference. Motion compensation in clinical datasets resulted in improved visibility of vascular structures, observed in maximum intensity projections of the contrast-enhanced liver vessel tree.

Conclusion

Targeted motion compensation for vascular imaging showed promising performance for increased identification of small vascular structures in presence of motion. The development of autofocus metrics and methods tailored to vascular imaging opens the way for reliable compensation of deformable motion while preserving the integrity of anatomical structures in the image.

Comparison of an in-house hybrid DIR method to NiftyReg on CBCT and CT images for head and neck cancer

An in-house hybrid deformable image registration (DIR) method, which combines free-form deformation (FFD) and the viscous fluid registration method, is proposed. Its results on the planning computed tomography (CT) and the day 1 treatment cone-beam CT (CBCT) image from 68 head and neck cancer patients are compared with the results of NiftyReg, which uses B-spline FFD alone. Several similarity metrics, the target registration error (TRE) of annotated points, as well as the Dice similarity coefficient (DSC) and Hausdorff distance (HD) of the propagated organs at risk are employed to analyze their registration accuracy. According to quantitative analysis on mutual information, normalized cross-correlation, and the absolute pixel value differences, the results of the proposed DIR are more similar to the CBCT images than the NiftyReg results. Smaller TRE of the annotated points is observed in the proposed method, and the overall mean TRE for the proposed method and NiftyReg was 2.34 and 2.98 mm, respectively (p < 0.001). The mean DSC in the larynx, spinal cord, oral cavity, mandible, and parotid given by the proposed method ranged from 0.78 to 0.91, significantly higher than the NiftyReg results (ranging from 0.77 to 0.90), and the HD was significantly lower compared to NiftyReg. Furthermore, the proposed method did not suffer from unrealistic deformations as the NiftyReg did in the visual evaluation. Meanwhile, the execution time of the proposed method was much higher than NiftyReg (96.98 ± 11.88 s vs. 4.60 ± 0.49 s). In conclusion, the in-house hybrid method gave better accuracy and more stable performance than NiftyReg.