Radiation doses and accuracy of navigated pedicle screw placement in cervical and thoracic spine surgery: a comparison of sliding gantry CT and mobile cone-beam CT in a homogeneous cohort

OBJECTIVE

Multiple solutions for navigation-guided pedicle screw placement are currently available. Intraoperative imaging techniques are invaluable for spinal surgery, but often there is little attention paid to patient radiation exposure. This study aimed to compare the applied radiation doses of sliding gantry CT (SGCT)- and mobile cone-beam CT (CBCT)-based pedicle screw placement for spinal instrumentation.

METHODS

The authors retrospectively analyzed 183 and 54 patients who underwent SGCT- or standard CBCT-based pedicle screw placement, respectively, for spinal instrumentation at their department between June 2019 and January 2020. SGCT uses an automated radiation dose adjustment.

RESULTS

Baseline characteristics, including the number of screws per patient and the number of instrumented levels, did not significantly differ between the two groups. Although the accuracy of screw placement according to Gertzbein-Robbins classification did not differ between the two groups, more screws had to be revised intraoperatively in the CBCT group (SGCT 2.7% vs CBCT 6.0%, p = 0.0036). Mean (± SD) radiation doses for the first (SGCT 484.0 ± 201.1 vs CBCT 687.4 ± 188.5 mGy*cm, p < 0.0001), second (SGCT 515.8 ± 216.3 vs CBCT 658.3 ± 220.1 mGy*cm, p < 0.0001), third (SGCT 531.3 ± 237.5 vs CBCT 641.6 ± 177.3 mGy*cm, p = 0.0140), and total (SGCT 1216.9 ± 699.3 vs CBCT 2000.3 ± 921.0 mGy*cm, p < 0.0001) scans were significantly lower for SGCT. This was also true for radiation doses per scanned level (SGCT 461.9 ± 429.3 vs CBCT 1004.1 ± 905.1 mGy*cm, p < 0.0001) and radiation doses per screw (SGCT 172.6 ± 110.1 vs CBCT 349.6 ± 273.4 mGy*cm, p < 0.0001).

CONCLUSIONS

The applied radiation doses were significantly lower using SGCT for navigated pedicle screw placement in spinal instrumentation. A modern CT scanner on a sliding gantry leads to lower radiation doses, especially through automated 3D radiation dose adjustment.

Three-dimensional carpal alignment: computer-aided CT analysis of carpal axes and normal ranges

Assessment of carpal alignment traditionally uses carpal bone axes measured on lateral radiographs. In this study, three-dimensional axes were defined for carpal bones using segmentation and numerical modelling of CT data of 121 neutrally positioned, asymptomatic wrists. The geometric axis was used for radius, scaphoid and capitate, whereas the axis based on a line perpendicular to the articular surface was used for the other carpal bones. Normal values of radiocarpal angles in the radial coordinate and the reliability of the computer-aided analysis are reported. The mean sagittal radiocarpal angles (positive in palmar direction) were as follows: scaphoid 58° (SD 10°), lunate 0° (SD 11°), triquetrum 12° (SD 8°), trapezium 17° (SD 8°), trapezoid -10° (SD 7°), capitate -17° (SD 9°) and hamate 2° (SD 7°). The mean coronal radiocarpal angles (positive in ulnar direction) were -42° (SD 9°), -20° (SD 4°), -49° (SD 4°), -32° (SD 6°), -16° (SD 5°), 2° (SD 7°) and 8° (SD 6°), respectively. The intra-observer reliability of the measurements was excellent (mean intraclass correlations coefficient 0.98). This study provides guidelines on how to measure and quantify carpal alignment three-dimensionally, and a database for the normal values. Together, these may be useful when analysing various wrist pathologies and kinematics of the wrist.

CBCT-Based Dose Monitoring and Adaptive Planning Triggers in Head and Neck PBS Proton Therapy

PURPOSE

To investigate the feasibility of using cone-beam computed tomography (CBCT)-derived synthetic CTs to monitor the daily dose and trigger a plan review for adaptive proton therapy (APT) in head and neck cancer (HNC) patients.

METHODS

For 84 HNC patients treated with proton pencil-beam scanning (PBS), same-day CBCT and verification CT (vfCT) pairs were retrospectively collected. The ground truth CT (gtCT) was created by deforming the vfCT to the same-day CBCT, and it was then used as a dosimetric baseline and for establishing plan review trigger recommendations. Two different synthetic CT algorithms were tested; the corrected CBCT (corrCBCT) was created using an iterative image correction method and the virtual CT (virtCT) was created by deforming the planning CT to the CBCT, followed by a low-density masking process. Clinical treatment plans were recalculated on the image sets for evaluation.

RESULTS

Plan review trigger criteria for adaptive therapy were established after closely reviewing the cohort data. Compared to the vfCT, the corrCBCT and virtCT reliably produced dosimetric data more similar to the gtCT. The average discrepancy in D99 for high-risk clinical target volumes (CTV) was 1.1%, 0.7%, and 0.4% and for standard-risk CTVs was 1.8%, 0.5%, and 0.5% for the vfCT, corrCBCT, and virtCT, respectively.

CONCLUSION

Streamlined APT has been achieved with the proposed plan review criteria and CBCT-based synthetic CT workflow.

Cone-Beam CT With Enhanced Needle Guidance and Augmented Fluoroscopy Overlay: Applications in Interventional Radiology

Digital flat-panel detector cone-beam CT (CBCT), introduced in the early 2000s, was historically used in interventional radiology primarily for liver-directed therapies. However, contemporary advanced imaging applications, including enhanced needle guidance and augmented fluoroscopy overlay, have evolved substantially over the prior decade and now work synergistically with CBCT guidance to overcome limitations encountered with other imaging modalities. CBCT with advanced imaging applications has become increasingly used to facilitate a broad range of minimally invasive procedures, particularly relating to pain and musculoskeletal interventions. Potential advantages of CBCT with advanced imaging applications include greater accuracy for complex needle paths, improved targeting in the presence of metal artifact, enhanced visualization during injection of contrast medium or cement, increased ease when space in the gantry is limited, and reduced radiation doses versus conventional CT guidance. Nonetheless, CBCT guidance remains underutilized, partly relating to lack of familiarity with the technique. This article describes the practical implementation of CBCT with enhanced needle guidance and augmented fluoroscopy overlay and depicts the technique's application for an array of interventional radiology procedures, including epidural steroid injections, celiac plexus block and neurolysis, pudendal block, spine ablation, percutaneous osseous ablation fixation and osteoplasty, biliary recanalization, and transcaval type II endoleak repair.

Evaluating the effect of setup uncertainty reduction and adaptation to geometric changes on normal tissue complication probability using online adaptive head and neck intensity modulated proton therapy

Objective. To evaluate the impact of setup uncertainty reduction (SUR) and adaptation to geometrical changes (AGC) on normal tissue complication probability (NTCP) when using online adaptive head and neck intensity modulated proton therapy (IMPT).Approach.A cohort of ten retrospective head and neck cancer patients with daily scatter corrected cone-beam CT (CBCT) was studied. For each patient, two IMPT treatment plans were created: one with a 3 mm setup uncertainty robustness setting and one with no explicit setup robustness. Both plans were recalculated on the daily CBCT considering three scenarios: the robust plan without adaptation, the non-robust plan without adaptation and the non-robust plan with daily online adaptation. Online-adaptation was simulated using an in-house developed workflow based on GPU-accelerated Monte Carlo dose calculation and partial spot-intensity re-optimization. Dose distributions associated with each scenario were accumulated on the planning CT, where NTCP models for six toxicities were applied. NTCP values from each scenario were intercompared to quantify the reduction in toxicity risk induced by SUR alone, AGC alone and SUR and AGC combined. Finally, a decision tree was implemented to assess the clinical significance of the toxicity reduction associated with each mechanism.Main results. For most patients, clinically meaningful NTCP reductions were only achieved when SUR and AGC were performed together. In these conditions, total reductions in NTCP of up to 30.48 pp were obtained, with noticeable NTCP reductions for aspiration, dysphagia and xerostomia (mean reductions of 8.25, 5.42 and 5.12 pp respectively). While SUR had a generally larger impact than AGC on NTCP reductions, SUR alone did not induce clinically meaningful toxicity reductions in any patient, compared to only one for AGC alone.SignificanceOnline adaptive head and neck proton therapy can only yield clinically significant reductions in the risk of long-term side effects when combining the benefits of SUR and AGC.

Image-guided radiotherapy (IGRT) in Lombardy, Italy: a survey by the Lombardy section of the Italian Association of Radiotherapy and Clinical Oncology (AIRO-Lombardy)

BACKGROUND

Image-guided radiation therapy (IGRT) has changed clinical practice. We proposed a survey to radiotherapy centers in Lombardy to picture the current clinical practice of its use.

RESEARCH DESIGN AND METHODS

The survey consisted of 32 multiple-choice questions, divided into five topics: type of hospital, patients treated in 2019, number of LINACs; presence of protocols and staff involved in IGRT; IGRT in stereotaxis; IGRT in non-stereotactic treatments; availability of medical and technical staff.

RESULTS

Twenty-seven directors answered (77%). Most centers (74%) have produced protocols to ensure uniformity in the IGRT process. The most widely used IGRT modality (92%) is cone-beam CT. Daily IGRT control is favored for prostate (100%), head and neck (87%), and lung (78%) neoplasms. The resident doctors can always perform supervised IGRT matching in only six centers. Radiation therapists perform IGRT controls only for some sites in 12 cases (44%) and always in 9 cases (33%). Radiation oncologists are present in real time, in most cases.

CONCLUSIONS

Today, IGRT can be considered standard practice but at the price of more time-consuming procedures. A balance between a fully physician-controlled process and an increased role for specifically trained RTTs is actively being sought.

Combining physics-based models with deep learning image synthesis and uncertainty in intraoperative cone-beam CT of the brain

BACKGROUND

Image-guided neurosurgery requires high localization and registration accuracy to enable effective treatment and avoid complications. However, accurate neuronavigation based on preoperative magnetic resonance (MR) or computed tomography (CT) images is challenged by brain deformation occurring during the surgical intervention.

PURPOSE

To facilitate intraoperative visualization of brain tissues and deformable registration with preoperative images, a 3D deep learning (DL) reconstruction framework (termed DL-Recon) was proposed for improved intraoperative cone-beam CT (CBCT) image quality.

METHODS

The DL-Recon framework combines physics-based models with deep learning CT synthesis and leverages uncertainty information to promote robustness to unseen features. A 3D generative adversarial network (GAN) with a conditional loss function modulated by aleatoric uncertainty was developed for CBCT-to-CT synthesis. Epistemic uncertainty of the synthesis model was estimated via Monte Carlo (MC) dropout. Using spatially varying weights derived from epistemic uncertainty, the DL-Recon image combines the synthetic CT with an artifact-corrected filtered back-projection (FBP) reconstruction. In regions of high epistemic uncertainty, DL-Recon includes greater contribution from the FBP image. Twenty paired real CT and simulated CBCT images of the head were used for network training and validation, and experiments evaluated the performance of DL-Recon on CBCT images containing simulated and real brain lesions not present in the training data. Performance among learning- and physics-based methods was quantified in terms of structural similarity (SSIM) of the resulting image to diagnostic CT and Dice similarity metric (DSC) in lesion segmentation compared to ground truth. A pilot study was conducted involving seven subjects with CBCT images acquired during neurosurgery to assess the feasibility of DL-Recon in clinical data.

RESULTS

CBCT images reconstructed via FBP with physics-based corrections exhibited the usual challenges to soft-tissue contrast resolution due to image non-uniformity, noise, and residual artifacts. GAN synthesis improved image uniformity and soft-tissue visibility but was subject to error in the shape and contrast of simulated lesions that were unseen in training. Incorporation of aleatoric uncertainty in synthesis loss improved estimation of epistemic uncertainty, with variable brain structures and unseen lesions exhibiting higher epistemic uncertainty. The DL-Recon approach mitigated synthesis errors while maintaining improvement in image quality, yielding 15%-22% increase in SSIM (image appearance compared to diagnostic CT) and up to 25% increase in DSC in lesion segmentation compared to FBP. Clear gains in visual image quality were also observed in real brain lesions and in clinical CBCT images.

CONCLUSIONS

DL-Recon leveraged uncertainty estimation to combine the strengths of DL and physics-based reconstruction and demonstrated substantial improvements in the accuracy and quality of intraoperative CBCT. The improved soft-tissue contrast resolution could facilitate visualization of brain structures and support deformable registration with preoperative images, further extending the utility of intraoperative CBCT in image-guided neurosurgery.

Orthodontic Treatment of Palatally Impacted Maxillary Canines with the Use of a Digitally Designed and 3D-Printed Metal Device

The purpose of this article is to present a computer designed and 3D-printed metal device, which was used for the surgical exposure and orthodontic treatment of maxillary palatally impacted canines. In two cases which presented a palatally impacted canine, a Cone-Beam Computed Tomography (CBCT) was acquired and an intraoral scanning was performed, to determine the exact location of the canine. Based on a digital model, a device leaning on the teeth and mucosa was designed to serve as a guiding tool for the oral surgeon to expose the crown of the canine and help the orthodontist to provide proper traction. The device was then 3D-printed in biocompatible dental alloy and placed in the patients' mouth. After the surgical exposure of the canine's crown in both cases, a gold chain apparatus was bonded on and it was mounted on the metal projection of the device through an elastic chain. Within 3 months of traction, the crown of the canines appeared in the patients' palate to the exact location that was predicted and guided. A 3D-designed and manufactured metal device, with information acquired by CBCT and intraoral scanning, can be used for the exposure and traction of palatally impacted canines.

TOPAS-imaging: extensions to the TOPAS simulation toolkit for medical imaging systems

Objective. The TOol for PArticle Simulation (TOPAS) is a Geant4-based Monte Carlo software application that has been used for both research and clinical studies in medical physics. So far, most users of TOPAS have focused on radiotherapy-related studies, such as modeling radiation therapy delivery systems or patient dose calculation. Here, we present the first set of TOPAS extensions to make it easier for TOPAS users to model medical imaging systems.Approach. We used the extension system of TOPAS to implement pre-built, user-configurable geometry components such as detectors (e.g. flat-panel and multi-planar detectors) for various imaging modalities and pre-built, user-configurable scorers for medical imaging systems (e.g. digitizer chain).Main results. We developed a flexible set of extensions that can be adapted to solve research questions for a variety of imaging modalities. We then utilized these extensions to model specific examples of cone-beam CT (CBCT), positron emission tomography (PET), and prompt gamma (PG) systems. The first of these new geometry components, the FlatImager, was used to model example CBCT and PG systems. Detected signals were accumulated in each detector pixel to obtain the intensity of x-rays penetrating objects or prompt gammas from proton-nuclear interaction. The second of these new geometry components, the RingImager, was used to model an example PET system. Positron-electron annihilation signals were recorded in crystals of the RingImager and coincidences were detected. The simulated data were processed using corresponding post-processing algorithms for each modality and obtained results in good agreement with the expected true signals or experimental measurement.Significance. The newly developed extension is a first step to making it easier for TOPAS users to build and simulate medical imaging systems. Together with existing TOPAS tools, this extension can help integrate medical imaging systems with radiotherapy simulations for image-guided radiotherapy.

A Cone Beam CT Study on the Correlation between Crestal Bone Loss and Periapical Disease

The aim of this study was to determine whether the degree of bone loss around teeth can be linked to the loss of vitality of adjacent teeth and periapical disease, which necessitates root canal treatments. Three hundred and twenty-one full maxilla cone-beam computed tomography scans were examined. The parameters investigated included the degree of crestal bone loss in relation to the cementoenamel junction, the presence/absence of apical periodontitis, and the presence/absence of root canal treatments. Out of the 2001 teeth examined, 696 (34.8%) showed evidence of crestal bone loss. The degree of crestal bone loss was classified as mild, moderate, or severe. A significant association (p < 0.001) was found between the presence of crestal bone loss around a tooth and root canal treatment of that tooth. It was found that it is more likely for teeth with crestal bone loss to be root canal treated compared to teeth with existing root canal treatment and healthy crestal bone levels. Furthermore, teeth with buccal or lingual crestal bone loss were significantly associated with a higher rate of periapical disease than teeth without crestal bone loss (p < 0.001). CBCT identified the severity of bone loss on all surfaces of the teeth, and the most common presentation was bone loss to the mid-root level. Teeth with crestal bone loss were significantly more likely to be associated with a higher rate of periapical disease. Teeth with crestal bone loss were more likely to be root treated than teeth with healthy crestal bone levels.

Advanced Imaging for Robotic Bronchoscopy: A Review

Recent advances in navigational platforms have led bronchoscopists to make major strides in diagnostic interventions for pulmonary parenchymal lesions. Over the last decade, multiple platforms including electromagnetic navigation and robotic bronchoscopy have allowed bronchoscopists to safely navigate farther into the lung parenchyma with increased stability and accuracy. Limitations persist, even with these newer technologies, in achieving a similar or higher diagnostic yield when compared to the transthoracic computed tomography (CT) guided needle approach. One of the major limitations to this effect is due to CT-to-body divergence. Real-time feedback that better defines the tool-lesion relationship is vital and can be obtained with additional imaging using radial endobronchial ultrasound, C-arm based tomosynthesis, cone-beam CT (fixed or mobile), and O-arm CT. Herein, we describe the role of this adjunct imaging with robotic bronchoscopy for diagnostic purposes, describe potential strategies to counteract the CT-to-body divergence phenomenon, and address the potential role of advanced imaging for lung tumor ablation.

A correlated sampling-based Monte Carlo simulation for fast CBCT iterative scatter correction

BACKGROUND

In recent years, cone-beam computed tomography (CBCT) has played an important role in medical imaging. However, the applications of CBCT are limited due to the severe scatter contamination. Conventional Monte Carlo (MC) simulation can provide accurate scatter estimation for scatter correction, but the expensive computational cost has always been the bottleneck of MC method in clinical application.

PURPOSE

In this work, an MC simulation method combined with a variance reduction technique called correlated sampling is proposed for fast iterative scatter correction.

METHODS

Correlated sampling exploits correlation between similar simulation systems to reduce the variance of interest quantities. Specifically, conventional MC simulation is first performed on the scatter-contaminated CBCT to generate the initial scatter signal. In the subsequent correction iterations, scatter estimation is then updated by applying correlated MC sampling to the latest corrected CBCT images by reusing the random number sequences of the task-related photons in conventional MC. Afterward, the corrected projections obtained by subtracting the scatter estimation from raw projections are utilized for FDK reconstruction. These steps are repeated until an adequate scatter correction is obtained. The performance of the proposed framework is evaluated by the accuracy of the scatter estimation, the quality of corrected CBCT images and efficiency.

RESULTS

Overall, the difference in mean absolute percentage error between scatter estimation with and without correlated sampling is 0.25% for full-fan case and 0.34% for half-fan case, respectively. In simulation studies, scatter artifacts are substantially eliminated, where the mean absolute error value is reduced from 15 to 2 HU in full-fan case and from 53 to 13 HU in half-fan case. Scatter-to-primary ratio is reduced to 0.02 for full-fan and 0.04 for half-fan, respectively. In phantom study, the contrast-to-noise ratio (CNR) is increased by a factor of 1.63, and the contrast is increased by a factor of 1.77. As for clinical studies, the CNR is improved by 11% and 14% for half-fan and full-fan, respectively. The contrast after correction is increased by 19% for half-fan and 44% for full-fan. Furthermore, root mean square error is also effectively reduced, especially from 78 to 4 HU for full-fan. Experimental results demonstrate that the figure of merit is improved between 23 and 43 folds when using correlated sampling. The proposed method takes less than 25 s for the whole iterative scatter correction process.

CONCLUSIONS

The proposed correlated sampling-based MC simulation method can achieve fast and accurate scatter correction for CBCT, making it suitable for real-time clinical use.

Dynamic cone-beam CT reconstruction using spatial and temporal implicit neural representation learning (STINR)

Objective. Dynamic cone-beam CT (CBCT) imaging is highly desired in image-guided radiation therapy to provide volumetric images with high spatial and temporal resolutions to enable applications including tumor motion tracking/prediction and intra-delivery dose calculation/accumulation. However, dynamic CBCT reconstruction is a substantially challenging spatiotemporal inverse problem, due to the extremely limited projection sample available for each CBCT reconstruction (one projection for one CBCT volume).Approach. We developed a simultaneous spatial and temporal implicit neural representation (STINR) method for dynamic CBCT reconstruction. STINR mapped the unknown image and the evolution of its motion into spatial and temporal multi-layer perceptrons (MLPs), and iteratively optimized the neuron weightings of the MLPs via acquired projections to represent the dynamic CBCT series. In addition to the MLPs, we also introduced prior knowledge, in the form of principal component analysis (PCA)-based patient-specific motion models, to reduce the complexity of the temporal mapping to address the ill-conditioned dynamic CBCT reconstruction problem. We used the extended-cardiac-torso (XCAT) phantom and a patient 4D-CBCT dataset to simulate different lung motion scenarios to evaluate STINR. The scenarios contain motion variations including motion baseline shifts, motion amplitude/frequency variations, and motion non-periodicity. The XCAT scenarios also contain inter-scan anatomical variations including tumor shrinkage and tumor position change.Main results. STINR shows consistently higher image reconstruction and motion tracking accuracy than a traditional PCA-based method and a polynomial-fitting-based neural representation method. STINR tracks the lung target to an average center-of-mass error of 1-2 mm, with corresponding relative errors of reconstructed dynamic CBCTs around 10%.Significance. STINR offers a general framework allowing accurate dynamic CBCT reconstruction for image-guided radiotherapy. It is a one-shot learning method that does not rely on pre-training and is not susceptible to generalizability issues. It also allows natural super-resolution. It can be readily applied to other imaging modalities as well.

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.

Quantitative Dual-Energy Imaging of Bone Marrow Edema Using Multisource Cone-Beam CT with Model-Based Decomposition

Purpose

We investigated the feasibility of dual-energy (DE) detection of bone marrow edema (BME) using a dedicated extremity cone-beam CT (CBCT) with a unique three-source x-ray unit. The sources can be operated at different energies to enable single-scan DE acquisitions. However, they are arranged parallel to the axis of rotation, resulting in incomplete sampling and precluding the application of DE projection-domain decompositions (PDD) for beam-hardening reduction. Therefore, we propose a novel combination of a model-based "one-step" DE two-material decomposition followed by a constrained image-domain change-of-basis to obtain virtual non-calcium (VNCa) images for BME detection.

Methods

DE projections were obtained using an "alternating-kV" protocol by operating the peripheral two sources of the CBCT system at low-energy (60 kV, 0.105 mAs/frame) and the central source at high-energy (100 kV, 0.028 mAs/frame), for a total of 600 frames over 216° of gantry rotation. Projections were processed with detector lag, glare and fast Monte Carlo (MC)-based iterative scatter corrections. Model-based material decomposition (MBMD) was then implemented to obtain aluminum (Al) and polyethylene (PE) volume fraction images with minimal beam-hardening. Statistical ray weights in MBMD were modified to account for regions with highly oblique sampling by the peripheral sources. To generate the VNCa maps, image-domain decomposition (IDD) constrained by the volume conservation principle (VCP) was performed to convert the Al and PE MBMD images into volume fractions of water, fat and cortical bone. Accuracy of BME detection was evaluated using physical phantom data acquired on the multi-source extremity CBCT scanner.

Results

The proposed framework estimated the volume of BME with ~10% error. The MC-based scatter corrections and the modified MBMD ray weights were essential to achieve such performance - the error without MC scatter corrections was >30%, whereas the uniformity of estimated VNCa images was 3x improved using the modified weights compared to the conventional weights.

Conclusions

The proposed DE decomposition framework was able to overcome challenges of high scatter and incomplete sampling to achieve BME detection on a CBCT system with axially-distributed x-ray sources.

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.