Fully automatic online geometric calibration for non-circular cone-beam CT orbits using fiducials with unknown placement.

BACKGROUND: Cone-beam CT (CBCT) with non-circular scanning orbits can improve image quality for 3D intraoperative image guidance. However, geometric calibration of such scans can be challenging. Existing methods typically require a prior image, specialized phantoms, presumed repeatable orbits, or long computation time. PURPOSE: We propose a novel fully automatic online geometric calibration algorithm that does not require prior knowledge of fiducial configuration. The algorithm is fast, accurate, and can accommodate arbitrary scanning orbits and fiducial configurations. METHODS: The algorithm uses an automatic initialization process to eliminate human intervention in fiducial localization and an iterative refinement process to ensure robustness and accuracy. We provide a detailed explanation and implementation of the proposed algorithm. Physical experiments on a lab test bench and a clinical robotic C-arm scanner were conducted to evaluate spatial resolution performance and robustness under realistic constraints. RESULTS: Qualitative and quantitative results from the physical experiments demonstrate high accuracy, efficiency, and robustness of the proposed method. The spatial resolution performance matched that of our existing benchmark method, which used a 3D-2D registration-based geometric calibration algorithm. CONCLUSIONS: We have demonstrated an automatic online geometric calibration method that delivers high spatial resolution and robustness performance. This methodology enables arbitrary scan trajectories and should facilitate translation of such acquisition methods in a clinical setting.

Faster and lower dose imaging: evaluating adaptive, constant gantry velocity and angular separation in fast low-dose 4D cone beam CT imaging.

BACKGROUND: The adoption of four-dimensional cone beam computed tomography (4DCBCT) for image-guided lung cancer radiotherapy is increasing, especially for hypofractionated treatments. However, the drawbacks of 4DCBCT include long scan times (∼240 s), inconsistent image quality, higher imaging dose than necessary, and streaking artifacts. With the emergence of linear accelerators that can acquire 4DCBCT scans in a short period of time (9.2 s) there is a need to examine the impact that these very fast gantry rotations have on 4DCBCT image quality. PURPOSE: This study investigates the impact of gantry velocity and angular separation between x-ray projections on image quality and its implication for fast low-dose 4DCBCT with emerging systems, such as the Varian Halcyon that provide fast gantry rotation and imaging. Large and uneven angular separation between x-ray projections is known to reduce 4DCBCT image quality through increased streaking artifacts. However, it is not known when angular separation starts degrading image quality. The study assesses the impact of constant and adaptive gantry velocity and determines the level when angular gaps impair image quality using state-of-the-art reconstruction methods. METHODS: This study considers fast low-dose 4DCBCT acquisitions (60-80 s, 200-projection scans). To assess the impact of adaptive gantry rotations, the angular position of x-ray projections from adaptive 4DCBCT acquisitions from a 30-patient clinical trial were analyzed (referred to as patient angular gaps). To assess the impact of angular gaps, variable and static angular gaps (20°, 30°, 40°) were introduced into evenly separated 200 projections (ideal angular separation). To simulate fast gantry rotations, which are on emerging linacs, constant gantry velocity acquisitions (9.2 s, 60 s, 120 s, 240 s) were simulated by sampling x-ray projections at constant intervals using the patient breathing traces from the ADAPT clinical trial (ACTRN12618001440213). The 4D Extended Cardiac-Torso (XCAT) digital phantom was used to simulate projections to remove patient-specific image quality variables. Image reconstruction was performed using Feldkamp-Davis-Kress (FDK), McKinnon-Bates (MKB), and Motion-Compensated-MKB (MCMKB) algorithms. Image quality was assessed using Structural Similarity-Index-Measure (SSIM), Contrast-to-Noise-Ratio (CNR), Signal-to-Noise-Ratio (SNR), Tissue-Interface-Width-Diaphragm (TIW-D), and Tissue-Interface-Width-Tumor (TIW-T). RESULTS: Patient angular gaps and variable angular gap reconstructions produced similar results to ideal angular separation reconstructions, while static angular gap reconstructions produced lower image quality metrics. For MCMKB-reconstructions, average patient angular gaps produced SSIM-0.98, CNR-13.6, SNR-34.8, TIW-D-1.5 mm, and TIW-T-2.0 mm, static angular gap 40° produced SSIM-0.92, CNR-6.8, SNR-6.7, TIW-D-5.7 mm, and TIW-T-5.9 mm and ideal produced SSIM-1.00, CNR-13.6, SNR-34.8, TIW-D-1.5 mm, and TIW-T-2.0 mm. All constant gantry velocity reconstructions produced lower image quality metrics than ideal angular separation reconstructions regardless of the acquisition time. Motion compensated reconstruction (MCMKB) produced the highest contrast images with low streaking artifacts. CONCLUSION: Very fast 4DCBCT scans can be acquired provided that the entire scan range is adaptively sampled, and motion-compensated reconstruction is performed. Importantly, the angular separation between x-ray projections within each individual respiratory bin had minimal effect on the image quality of fast low-dose 4DCBCT imaging. The results will assist the development of future 4DCBCT acquisition protocols that can now be achieved in very short time frames with emerging linear accelerators.

Revealing pelvic structures in the presence of metal hip prothesis via non-circular CBCT orbits.

As the expansion of Cone Beam CT (CBCT) to new interventional procedures continues, the burdensome challenge of metal artifacts remains. Photon starvation and beam hardening from metallic implants and surgical tools in the field of view can result in the anatomy of interest being partially or fully obscured by imaging artifacts. Leveraging the flexibility of modern robotic CBCT imaging systems, implementing non-circular orbits designed for reducing metal artifacts by ensuring data-completeness during acquisition has become a reality. Here, we investigate using non-circular orbits to reduce metal artifacts arising from metallic hip prostheses when imaging pelvic anatomy. As a first proof-of-concept, we implement a sinusoidal and a double-circle-arc orbit on a CBCT test bench, imaging a physical pelvis phantom, with two metal hip prostheses, housing a 3D-printed iodine-filled radial line-pair target. A standard circular orbit implemented with the CBCT test bench acted as comparator. Imaging data collection and processing, geometric calibration and image reconstruction was completed using in-house developed software programs. Imaging with the standard circular orbit, image artifacts were observed in the pelvic bones and only 33 out of the possible 45 line-pairs of the radial line-pair target were partially resolvable in the reconstructed images. Comparatively, imaging with both the sinusoid and double-circle-arc orbits reduced artifacts in the surrounding anatomy and enabled all 45 line-pairs to be visibly resolved in the reconstructed images. These results indicate the potential of non-circular orbits to assist in revealing previously obstructed structures in the pelvic region in the presence of metal hip prosthesis.

X-ray source arrays for volumetric imaging during radiotherapy treatment.

This work presents a novel hardware configuration for radiotherapy systems to enable fast 3D X-ray imaging before and during treatment delivery. Standard external beam radiotherapy linear accelerators (linacs) have a single X-ray source and detector located at ± 90° from the treatment beam respectively. The entire system can be rotated around the patient acquiring multiple 2D X-ray images to create a 3D cone-beam Computed Tomography (CBCT) image before treatment delivery to ensure the tumour and surrounding organs align with the treatment plan. Scanning with a single source is slow relative to patient respiration or breath holds and cannot be performed during treatment delivery, limiting treatment delivery accuracy in the presence of patient motion and excluding some patients from concentrated treatment plans that would be otherwise expected to have improved outcomes. This simulation study investigated whether recent advances in carbon nanotube (CNT) field emission source arrays, high frame rate (60 Hz) flat panel detectors and compressed sensing reconstruction algorithms could circumvent imaging limitations of current linacs. We investigated a novel hardware configuration incorporating source arrays and high frame rate detectors into an otherwise standard linac. We investigated four potential pre-treatment scan protocols that could be achieved in a 17 s breath hold or 2-10 1 s breath holds. Finally, we demonstrated for the first time volumetric X-ray imaging during treatment delivery by using source arrays, high frame rate detectors and compressed sensing. Image quality was assessed quantitatively over the CBCT geometric field of view as well as across each axis through the tumour centroid. Our results demonstrate that source array imaging enables larger volumes to be imaged with acquisitions as short as 1 s albeit with reduced image quality arising from lower photon flux and shorter imaging arcs.

Technical note: Extended longitudinal and lateral 3D imaging with a continuous dual-isocenter CBCT scan.

BACKGROUND: The clinical benefits of intraoperative cone beam CT (CBCT) during orthopedic procedures include (1) improved accuracy for procedures involving the placement of hardware and (2) providing immediate surgical verification. PURPOSE: Orthopedic interventions often involve long and wide anatomical sites (e.g., lower extremities). Therefore, in order to ensure that the clinical benefits are available to all orthopedic procedures, we investigate the feasibility of a novel imaging trajectory to simultaneously expand the CBCT field-of-view longitudinally and laterally. METHODS: A continuous dual-isocenter imaging trajectory was implemented on a clinical robotic CBCT system using additional real-time control hardware. The trajectory consisted of 200° circular arcs separated by alternating lateral and longitudinal table translations. Due to hardware constraints, the direction of rotation (clockwise/anticlockwise) and lateral table translation (left/right) was reversed every 400°. X-ray projections were continuously acquired at 15 frames/s throughout all movements. A whole-body phantom was used to verify the trajectory. As comparator, a series of conventional large volume acquisitions were stitched together. Image quality was quantified using Root Mean Square Deviation (RMSD), Mean Absolute Percentage Deviation (MAPD), Structural Similarity Index Metric (SSIM) and Contrast-to-Noise Ratio (CNR). RESULTS: The imaging volume produced by the continuous dual-isocenter trajectory had dimensions of L = 95 cm × W = 45 cm × H = 45 cm. This enabled the hips to the feet of the whole-body phantom to be captured in approximately half the imaging dose and acquisition time of the 11 stitched conventional acquisitions required to match the longitudinal and lateral imaging dimensions. Compared to the stitched conventional images, the continuous dual-isocenter acquisition had RMSD of 4.84, MAPD of 6.58% and SSIM of 0.99. The CNR of the continuous dual-isocenter and stitched conventional acquisitions were 1.998 and 1.999, respectively. CONCLUSION: Extended longitudinal and lateral intraoperative volumetric imaging is feasible on clinical robotic CBCT systems.

Real-time radial reconstruction with domain transform manifold learning for MRI-guided radiotherapy.

BACKGROUND: MRI-guidance techniques that dynamically adapt radiation beams to follow tumor motion in real time will lead to more accurate cancer treatments and reduced collateral healthy tissue damage. The gold-standard for reconstruction of undersampled MR data is compressed sensing (CS) which is computationally slow and limits the rate that images can be available for real-time adaptation. PURPOSE: Once trained, neural networks can be used to accurately reconstruct raw MRI data with minimal latency. Here, we test the suitability of deep-learning-based image reconstruction for real-time tracking applications on MRI-Linacs. METHODS: We use automated transform by manifold approximation (AUTOMAP), a generalized framework that maps raw MR signal to the target image domain, to rapidly reconstruct images from undersampled radial k-space data. The AUTOMAP neural network was trained to reconstruct images from a golden-angle radial acquisition, a benchmark for motion-sensitive imaging, on lung cancer patient data and generic images from ImageNet. Model training was subsequently augmented with motion-encoded k-space data derived from videos in the YouTube-8M dataset to encourage motion robust reconstruction. RESULTS: AUTOMAP models fine-tuned on retrospectively acquired lung cancer patient data reconstructed radial k-space with equivalent accuracy to CS but with much shorter processing times. Validation of motion-trained models with a virtual dynamic lung tumor phantom showed that the generalized motion properties learned from YouTube lead to improved target tracking accuracy. CONCLUSION: AUTOMAP can achieve real-time, accurate reconstruction of radial data. These findings imply that neural-network-based reconstruction is potentially superior to alternative approaches for real-time image guidance applications.

CArdiac and REspiratory adaptive Computed Tomography (CARE-CT): a proof-of-concept digital phantom study.

Current respiratory 4DCT imaging for high-dose rate thoracic radiotherapy treatments are negatively affected by the complex interaction of cardiac and respiratory motion. We propose an imaging method to reduce artifacts caused by thoracic motion, CArdiac and REspiratory adaptive CT (CARE-CT), that monitors respiratory motion and ECG signals in real-time, triggering CT acquisition during combined cardiac and respiratory bins. Using a digital phantom, conventional 4DCT and CARE-CT acquisitions for nineteen patient-measured physiological traces were simulated. Ten respiratory bins were acquired for conventional 4DCT scans and ten respiratory bins during cardiac diastole were acquired for CARE-CT scans. Image artifacts were quantified for 10 common thoracic organs at risk (OAR) substructures using the differential normalized cross correlation between axial slices (ΔNCC), mean squared error (MSE) and sensitivity. For all images, on average, CARE-CT improved the ΔNCC for 18/19 and the MSE and sensitivity for all patient traces. The ΔNCC was reduced for all cardiac OARs (mean reduction 21%). The MSE was reduced for all OARs (mean reduction 36%). In the digital phantom study, the average scan time was increased from 1.8 ± 0.4 min to 7.5 ± 2.2 min with a reduction in average beam on time from 98 ± 28 s to 45 s using CARE-CT compared to conventional 4DCT. The proof-of-concept study indicates the potential for CARE-CT to image the thorax in real-time during the cardiac and respiratory cycle simultaneously, to reduce image artifacts for common thoracic OARs.

System requirements to improve adaptive 4-dimensional computed tomography (4D CT) imaging.

Four-Dimensional Computed Tomography (4D CT) is of increasing importance in stereotactic body radiotherapy (SBRT) treatments affected by respiratory motion. However, 4D CT images are commonly impacted by irregular breathing, causing image artifacts that can propagate through to treatment, negatively effecting local control. REspiratory Adaptive CT (REACT) is a real-time gating method demonstrated to reduce motion artifacts by avoiding imaging during irregular respiration. Despite artifact reduction seen throughin silicoand clinical phantom-based studies, REACT has not been able to remove all artifacts. Here, we explore several hardware and software latencies (gantry rotation time, couch shifts, acquisition delays and phase calculation method) inherently linked to REACT and 4D CT in general and investigate their contribution to artifacts beyond those caused by irregular breathing. Imaging was simulated using the digital extended cardiac-torso (XCAT) phantom for fifty patient-measured respiratory traces. Imaging protocols included conventional cine 4D CT and five REACT scans with systematically varied parameters to test the effect of different latencies on artifacts. Artifacts were quantified by comparing the image normalized cross correlation across couch transition points and determining the volume error compared to a static phantom ground truth both as a total error and individually across pixel rows in the main plane of motion. Artifacts were determined for each lung, the whole heart and lung tumour and were compared back to conventional 4D CT and REACT with standard clinical scanning parameters. The gantry rotation time and acquisition delay were found to have the largest impact on reducing image artifacts and should be the focus of future development. The phase calculation method was also found to influence motion artifacts and should potentially be assessed on a patient-to-patient basis. Finally, the correlation between an increase in artifacts and baseline drift suggests that longer scan times allowing drift to occur may impact image quality.

Fast CBCT Reconstruction using Convolutional Neural Networks for Arbitrary Robotic C-arm Orbits.

Cone-beam CT (CBCT) with non-circular acquisition orbits has the potential to improve image quality, increase the field-of view, and facilitate minimal interference within an interventional imaging setting. Because time is of the essence in interventional imaging scenarios, rapid reconstruction methods are advantageous. Model-Based Iterative Reconstruction (MBIR) techniques implicitly handle arbitrary geometries; however, the computational burden for these approaches is particularly high. The aim of this work is to extend a previously proposed framework for fast reconstruction of non-circular CBCT trajectories. The pipeline combines a deconvolution operation on the backprojected measurements using an approximate, shift-invariant system response prior to processing with a Convolutional Neural Network (CNN). We trained and evaluated the CNN for this approach using 1800 randomized arbitrary orbits. Noisy projection data were formed from 1000 procedurally generated tetrahedral phantoms as well as anthropomorphic data in the form of 800 CT and CBCT images from the Lung Image Database Consortium Image Collection (LIDC). Using this proposed reconstruction pipeline, computation time was reduced by 90% as compared to MBIR with only minor differences in performance. Quantitative comparisons of nRMSE, FSIM and SSIM are reported. Performance was consistent for projection data simulated with acquisition orbits the network has not previously been trained on. These results suggest the potential for fast processing of arbitrary CBCT trajectory data with reconstruction times that are clinically relevant and applicable - facilitating the application of non-circular orbits in CT image-guided interventions and intraoperative imaging.

Extended Intraoperative Longitudinal 3-Dimensional Cone Beam Computed Tomography Imaging With a Continuous Multi-Turn Reverse Helical Scan.

OBJECTIVES: Cone beam computed tomography (CBCT) imaging is becoming an indispensable intraoperative tool; however, the current field of view prevents visualization of long anatomical sites, limiting clinical utility. Here, we demonstrate the longitudinal extension of the intraoperative CBCT field of view using a multi-turn reverse helical scan and assess potential clinical utility in interventional procedures. MATERIALS AND METHODS: A fixed-room robotic CBCT imaging system, with additional real-time control, was used to implement a multi-turn reverse helical scan. The scan consists of C-arm rotation, through a series of clockwise and anticlockwise rotations, combined with simultaneous programmed table translation. The motion properties and geometric accuracy of the multi-turn reverse helical imaging trajectory were examined using a simple geometric phantom. To assess potential clinical utility, a pedicle screw posterior fixation procedure in the thoracic spine from T1 to T12 was performed on an ovine cadaver. The multi-turn reverse helical scan was used to provide postoperative assessment of the screw insertion via cortical breach grading and mean screw angle error measurements (axial and sagittal) from 2 observers. For all screw angle measurements, the intraclass correlation coefficient was calculated to determine observer reliability. RESULTS: The multi-turn reverse helical scans took 100 seconds to complete and increased the longitudinal coverage by 370% from 17 cm to 80 cm. Geometric accuracy was examined by comparing the measured to actual dimensions (0.2 ± 0.1 mm) and angles (0.2 ± 0.1 degrees) of a simple geometric phantom, indicating that the multi-turn reverse helical scan provided submillimeter and degree accuracy with no distortion. During the pedicle screw procedure in an ovine cadaver, the multi-turn reverse helical scan identified 4 cortical breaches, confirmed via the postoperative CT scan. Directly comparing the screw insertion angles (n = 22) measured in the postoperative multi-turn reverse helical and CT scans revealed an average difference of 3.3 ± 2.6 degrees in axial angle and 1.9 ± 1.5 degrees in the sagittal angle from 2 expert observers. The intraclass correlation coefficient was above 0.900 for all measurements (axial and sagittal) across all scan types (conventional CT, multi-turn reverse helical, and conventional CBCT), indicating excellent reliability between observers. CONCLUSIONS: Extended longitudinal field-of-view intraoperative 3-dimensional imaging with a multi-turn reverse helical scan is feasible on a clinical robotic CBCT imaging system, enabling long anatomical sites to be visualized in a single image, including in the presence of metal hardware.

Non-circular CBCT orbit design and realization on a clinical robotic C-arm for metal artifact reduction.

Metal artifacts have been a difficult challenge for cone-beam CT (CBCT), especially for intraoperative imaging. Metal surgical tools and implants are often present in the field of view and can attenuate X-rays so heavily that they essentially create a missing-data problem. Recently, an increasing number of intra-operative imaging systems such as robotic C-arms are capable of non-circular orbits for data acquisition. Such trajectories can potentially improve sampling and the degree of data completeness to solve the metal-induced missing-data problem, thereby reducing or eliminating the associated image artifacts. In this work, we extend our prior theoretical and experimental work and implement non-circular orbits for metal artifact reduction on a clinical robotic C-arm (Siemens Artis zeego). To maximize the potential for clinical translation, we restrict our implementation to standard built-in motion and data collection functions, also available on other zeego systems, and work within the physical constraints and limitations on positioning and motion. Customized software tools for data extraction, processing, calibration, and reconstruction are used. We demonstrate example non-circular orbits and the resulting image quality using a phantom containing pedicle screws for spine fixation. As compared with a standard circular CBCT orbit, these non-circular orbits exhibit significantly reduced metal artifacts. These results suggest a high potential for image quality improvements for intraoperative CBCT imaging when metal tools or implants are present in the field-of-view.

Generating patient-matched 3D-printed pedicle screw and laminectomy drill guides from Cone Beam CT images: Studies in ovine and porcine cadavers.

BACKGROUND: The emergence of robotic Cone Beam Computed Tomography (CBCT) imaging systems in trauma departments has enabled 3D anatomical assessment of musculoskeletal injuries, supplementing conventional 2D fluoroscopic imaging for examination, diagnosis, and treatment planning. To date, the primary focus has been on trauma sites in the extremities. PURPOSE: To determine if CBCT images can be used during the treatment planning process in spinal instrumentation and laminectomy procedures, allowing accurate 3D-printed pedicle screw and laminectomy drill guides to be generated for the cervical and thoracic spine. METHODS: The accuracy of drill guides generated from CBCT images was assessed using animal cadavers (ovine and porcine). Preoperative scans were acquired using a robotic CBCT C-arm system, the Siemens ARTIS pheno (Siemens Healthcare, GmbH, Germany). The CBCT images were imported into 3D-Slicer version 4.10.2 (www.slicer.org) where vertebral models and specific guides were developed and subsequently 3D-printed. In the ovine cadaver, 11 pedicle screw guides from the T1-T5 and T7-T12 vertebra and six laminectomy guides from the C2-C7 vertebra were planned and printed. In the porcine cadaver, nine pedicle screw guides from the C3-T4 vertebra were planned and printed. For the pedicle screw guides, accuracy was assessed by three observers according to pedicle breach via the Gertzbein-Robbins grading system as well as measured mean axial and sagittal screw error via postoperative CBCT and CT scans. For the laminectomies, the guides were designed to leave 1 mm of lamina. The average thickness of the lamina at the mid-point was used to assess the accuracy of the guides, measured via postoperative CBCT and CT scans from three observers. For all measurements, the intraclass correlation coefficient (ICC) was calculated to determine observer reliability. RESULTS: Compared with the planned screw angles for both the ovine and porcine procedures (n = 32), the mean axial and sagittal screw error measured on the postoperative CBCT scans from three observers were 3.9 ± 1.9° and 1.8 ± 0.8°, respectively. The ICC among the observes was 0.855 and 0.849 for the axial and sagittal measurements, respectively, indicating good reliability. In the ovine cadaver, directly comparing the measured axial and sagittal screw angle of the visible screws (n = 14) in the postoperative CBCT and conventional CT scans from three observers revealed an average difference 1.9 ± 1.0° in axial angle and 1.8 ± 1.0° in the sagittal angle. The average thickness of the lamina at the middle of each vertebra, as measured on-screen in the postoperative CBCT scans by three observes was 1.6 ± 0.2 mm. The ICC among observers was 0.693, indicating moderate reliability. No lamina breaches were observed in the postoperative images. CONCLUSION: Here, CBCT images have been used to generate accurate 3D-printed pedicle screw and laminectomy drill guides for use in the cervical and thoracic spine. The results demonstrate sufficient precision compared with those previously reported, generated from standard preoperative CT and MRI scans, potentially expanding the treatment planning capabilities of robotic CBCT imaging systems in trauma departments and operating rooms.

Reducing 4DCBCT imaging dose and time: exploring the limits of adaptive acquisition and motion compensated reconstruction.

This study investigates the dose and time limits of adaptive 4DCBCT acquisitions (adaptive-acquisition) compared with current conventional 4DCBCT acquisition (conventional-acquisition). We investigate adaptive-acquisitions as low as 60 projections (∼25 s scan, 6 projections per respiratory phase) in conjunction with emerging image reconstruction methods. 4DCBCT images from 20 patients recruited into the adaptive CT acquisition for personalized thoracic imaging clinical study (NCT04070586) were resampled to simulate faster and lower imaging dose acquisitions. All acquisitions were reconstructed using Feldkamp-Davis-Kress (FDK), McKinnon-Bates (MKB), motion compensated FDK (MCFDK), motion compensated MKB (MCMKB) and simultaneous motion estimation and image reconstruction (SMEIR) algorithms. All reconstructions were compared against conventional-acquisition 4DFDK-reconstruction using Structural SIMilarity Index (SSIM), signal-to-noise ratio (SNR), contrast-to-noise-ratio (CNR), tissue interface sharpness diaphragm (TIS-D), tissue interface sharpness tumor (TIS-T) and center of mass trajectory (COMT) for difference in diaphragm and tumor motion. All reconstruction methods using 110-projection adaptive-acquisition (11 projections per respiratory phase) had a SSIM of greater than 0.92 relative to conventional-acquisition 4DFDK-reconstruction. Relative to conventional-acquisition 4DFDK-reconstruction, 110-projection adaptive-acquisition MCFDK-reconstructions images had 60% higher SNR, 10% higher CNR, 30% higher TIS-T and 45% higher TIS-D on average. The 110-projection adaptive-acquisition SMEIR-reconstruction images had 123% higher SNR, 90% higher CNR, 96% higher TIS-T and 60% higher TIS-D on average. The difference in diaphragm and tumor motion compared to conventional-acquisition 4DFDK-reconstruction was within submillimeter accuracy for all acquisition reconstruction methods. Adaptive-acquisitions resulted in faster scans with lower imaging dose and equivalent or improved image quality compared to conventional-acquisition. Adaptive-acquisition with motion compensated-reconstruction enabled scans with as low as 110 projections to deliver acceptable image quality. This translates into 92% lower imaging dose and 80% less scan time than conventional-acquisition.

Practical workflow for arbitrary non-circular orbits for CT with clinical robotic C-arms.

Non-circular orbits in cone-beam CT (CBCT) imaging are increasingly being studied for potential benefits in field-of-view, dose reduction, improved image quality, minimal interference in guided procedures, metal artifact reduction, and more. While modern imaging systems such as robotic C-arms are enabling more freedom in potential orbit designs, practical implementation on such clinical systems remains challenging due to obstacles in critical stages of the workflow, including orbit realization, geometric calibration, and reconstruction. In this work, we build upon previous successes in clinical implementation and address key challenges in the geometric calibration stage with a novel calibration method. The resulting workflow eliminates the need for prior patient scans or dedicated calibration phantoms, and can be conducted in clinically relevant processing times.

Reducing 4DCBCT scan time and dose through motion compensated acquisition and reconstruction.

Conventional 4DCBCT captures 1320 projections across 4 min. Adaptive 4DCBCT has been developed to reduce imaging dose and scan time. This study investigated reconstruction algorithms that best complement adaptive 4DCBCT acquisition for reducing imaging dose and scan time whilst maintaining or improving image quality compared to conventional 4DCBCT acquisition using real patient data from the first 10 adaptive 4DCBCT patients. Adaptive 4DCBCT was implemented in the ADaptive CT Acquisition for Personalized Thoracic imaging clinical trial. Adaptive 4DCBCT modulates gantry rotation speed and kV acquisition rate in response to the patient's real-time respiratory signal, ensuring even angular spacing between projections at each respiratory phase. We examined the first 10 lung cancer radiotherapy patients that received adaptive 4DCBCT. Fast, 200-projection scans over 60-80 s, and slower, 600-projection scans over ∼240 s, were obtained after routine patient treatment and compared against conventional 4DCBCT acquisition. Adaptive 4DCBCT acquisitions were reconstructed using Feldkamp-Davis-Kress (FDK), McKinnon-Bates (MKB), Motion Compensated FDK (MCFDK) and Motion Compensated MKB (MCMKB) algorithms. Reconstructions were assessed via, Structural SIMilarity (SSIM), Signal-to-Noise-Ratio (SNR), Contrast-to-Noise-Ratio (CNR), Tissue Interface Sharpness of Diaphragm (TIS-D) and Tumor (TIS-T). The 200- and 600-projection adaptive 4DCBCT acquisition corresponded to 85% and 55% reduction in imaging dose, shorter and similar scan times of approximately 90 s and 236 s respectively, compared to conventional 4DCBCT acquisition. 200- and 600-projection adaptive 4DCBCT reconstructions achieved more than 0.900 SSIM relative to conventional 4DCBCT acquisition. Compared to conventional 4DCBCT acquisition, 200-projection adaptive 4DCBCT reconstructions achieved higher SNR, CNR, TIS-T, TIS-D with motion compensated algorithms, MCFDK (208%, 159%, 174%, 247%) and MCMKB (214%, 173%, 266%, 245%) respectively. The 200-projection adaptive 4DCBCT MCFDK- and MCMKB-reconstruction results show image quality improvements are possible even with 85% fewer projections acquired. We established acquisition-reconstruction protocols that provide substantial reductions in imaging time and dose whilst improving image quality.

Adaptive CaRdiac cOne BEAm computed Tomography (ACROBEAT): Developing the next generation of cardiac cone beam CT imaging.

PURPOSE: An important factor when considering the use of interventional cone beam computed tomography (CBCT) imaging during cardiac procedures is the trade-off between imaging dose and image quality. Accordingly, Adaptive CaRdiac cOne BEAm computed Tomography (ACROBEAT) presents an alternative acquisition method, adapting the gantry velocity and projection rate of CBCT imaging systems in accordance with a patient's electrocardiogram (ECG) signal in real-time. The aim of this study was to experimentally investigate that ACROBEAT acquisitions deliver improved image quality compared to conventional cardiac CBCT imaging protocols with fewer projections acquired. METHODS: The Siemens ARTIS pheno (Siemens Healthcare, GmbH, Germany), a robotic CBCT C-arm system, was used to compare ACROBEAT with a commercially available conventional cardiac imaging protocol that utilizes multisweep retrospective ECG-gated acquisition. For ACROBEAT, real-time control of the gantry position was enabled through the Siemens Test Automation Control system. ACROBEAT and conventional image acquisitions of the CIRS Dynamic Cardiac Phantom were performed, using five patient-measured ECG traces. The traces had average heart rates of 56, 64, 76, 86, and 100 bpm. The total number of acquired projections was compared between the ACROBEAT and conventional acquisition methods. The image quality was assessed via the contrast-to-noise ratio (CNR), structural similarity index (SSIM), and root-mean square error (RMSE). RESULTS: Compared to the conventional protocol, ACROBEAT reduced the total number of projections acquired by 90%. The visual image quality provided by the ACROBEAT acquisitions, across all traces, matched or improved compared to conventional acquisition and was independent of the patient's heart rate. Across all traces, ACROBEAT averaged 1.4 times increase in the CNR, a 23% increase in the SSIM and a 29% decrease in the RMSE compared to conventional and was independent of the patient's heart rate. CONCLUSION: Adaptive patient imaging is feasible on a clinical robotic CBCT system, delivering higher quality images while reducing the number of projections acquired by 90% compared to conventional cardiac imaging protocols.

Minimizing 4DCBCT imaging dose and scan time with Respiratory Motion Guided 4DCBCT: a pre-clinical investigation.

Current conventional 4D Cone Beam Computed Tomography (4DCBCT) imaging is hampered by inconsistent patient breathing that leads to long scan times, reduced image quality and high imaging dose. To address these limitations, Respiratory Motion Guided 4D cone beam computed tomography (RMG-4DCBCT) uses mathematical optimization to adapt the gantry rotation speed and projection acquisition rate in real-time in response to changes in the patient's breathing rate. Here, RMG-4DCBCT is implemented on an Elekta Synergy linear accelerator to determine the minimum achievable imaging dose. 8 patient-measured breathing traces were programmed into a 1D motion stage supporting a 3D-printed anthropomorphic thorax phantom. The respiratory phase and current gantry position were calculated in real-time with the RMG-4DCBCT software, which in turn modulated the gantry rotation speed and suppressed projection acquisition. Specifically, the effect of acquiring 20, 25, 30, 35 and 40 projections/respiratory phase bin RMG scans on scan time and image quality was assessed. Reconstructed image quality was assessed via the contrast-to-noise ratio (CNR) and the Edge Response Width (ERW) metrics. The performance of the system in terms of gantry control accuracy was also assessed via an analysis of the angular separation between adjacent projections. The median CNR increased linearly from 5.90 (20 projections/bin) to 8.39 (40 projections/bin). The ERW did not significantly change from 1.08 mm (20 projections/bin) to 1.07 mm (40 projections/bin), indicating the sharpness is not dependent on the total number of projections acquired. Scan times increased with increasing total projections and slower breathing rates. Across all 40 RMG-4DCBCT scans performed, the average difference in the acquired and desired angular separation between projections was 0.64°. RMG-4DCBCT provides the opportunity to enable fast low-dose 4DCBCT (∼70 s, 200 projections), without compromising on current clinical image quality.

Toward improved 3D carotid artery imaging with Adaptive CaRdiac cOne BEAm computed Tomography (ACROBEAT).

PURPOSE: Interventional treatments of aneurysms in the carotid artery are increasingly being supplemented with three-dimensional (3D) x-ray imaging. The 3D imaging provides additional information on device sizing and stent malapposition during the procedure. Standard 3D x-ray image acquisition is a one-size fits all model, exposing patients to additional radiation and results in images that may have cardiac-induced motion blur around the artery. Here, we investigate the potential of a novel dynamic imaging technique Adaptive CaRdiac cOne BEAm computed Tomography (ACROBEAT) to personalize image acquisition by adapting the gantry velocity and projection rate in real-time to changes in the patient's electrocardiogram (ECG) trace. METHODS: We compared the total number of projections acquired, estimated carotid artery widths and image quality between ACROBEAT and conventional (single rotation fixed gantry velocity and acquisition rate, no ECG-gating) scans in a simulation study and a proof-of-concept physical phantom experimental study. The simulation study dataset consisted of an XCAT digital software phantom programmed with five patient-measured ECG traces and artery motion curves. The ECG traces had average heart rates of 56, 64, 76, 86, and 100 bpm. To validate the concept experimentally, we designed and manufactured the physical phantom from an 8-mm diameter silicon rubber tubing cast into Phytagel. An artery motion curve and the ECG trace with an average heart rate of 56 bpm was passed through the phantom. To implement ACROBEAT on the Siemens ARTIS pheno angiography system for the proof-of-concept experimental study, the Siemens Test Automation Control System was used. The total number of projections acquired and estimated carotid artery widths were compared between the ACROBEAT and conventional scans. As the ground truth was available for the simulation studies, the image quality metrics of Root Mean Square Error (RMSE) and Structural Similarity Index (SSIM) were also utilized to assess image quality. RESULTS: In the simulation study, on average, ACROBEAT reduced the number of projections acquired by 63%, reduced carotid width estimation error by 65%, reduced RMSE by 11% and improved SSIM by 27% compared to conventional scans. In the proof-of-concept experimental study, ACROBEAT enabled a 60% reduction in the number of projections acquired and reduced carotid width estimation error by 69% compared to a conventional scan. CONCLUSION: A simulation and proof-of-concept experimental study was completed applying a novel dynamic imaging protocol, ACROBEAT, to imaging the carotid artery. The ACROBEAT results showed significantly improved image quality with fewer projections, offering potential applications to intracranial interventional procedures negatively affected by cardiac motion.

Dual cardiac and respiratory gated thoracic imaging via adaptive gantry velocity and projection rate modulation on a linear accelerator: A Proof-of-Concept Simulation Study.

PURPOSE: Cardiac motion is typically not accounted for during pretreatment imaging for central lung and mediastinal tumors. However, cardiac induced tumor motion averages 5.8 mm for esophageal tumors and 3-5 mm for some lung tumors, which can result in positioning errors. Our aim is to reduce both cardiac- and respiratory-induced motion artifacts in thoracic cone beam computed tomography (CBCT) images through gantry velocity and projection rate modulation on a standard linear accelerator (linac). METHODS: The acquisition of dual cardiac and respiratory gated CBCT thoracic images was simulated using the XCAT phantom with patient-measured respiratory and ECG traces. The gantry velocity and projection rate were modulated based on the cardiac and respiratory signals to maximize the angular consistency between adjacent projections in the gated cardiac-respiratory bin. The mechanical limitations of a gantry-mounted CBCT system were investigated. For our protocol, images were acquired during the 60%-80% window of cardiac phase and 20% displacement either side of peak exhale of the respiratory cycle. The comparator method was the respiratory-only gated CBCT acquisition with constant gantry speed and projection rate in routine use for respiratory correlated four-dimensional (4D) CBCT. All images were reconstructed using the Feldkamp-Davis-Kress (FDK) algorithm. The methods were compared in terms of image sharpness as measured using the edge response width (ERW) and contrast as measured using the contrast to noise ratio (CNR). The effects of the total number of projections acquired and magnitude of cardiac motion on scan time and image quality were also investigated. RESULTS: Median total scan times with our protocol ranged from 117 s (40 projections) through to 296 s (100 projections), compared with 240 s for the conventional protocol (1320 projections). The scan times were dictated by the number of projections acquired, heart rate, length of the respiratory cycle and mechanical constraints of the CBCT system. Our protocol was able to provide between 8% and 43% decrease in the median value of the ERW in the anterior/posterior (AP) direction across the 17 traces when there was 0.5 cm of cardiac motion and between 35% and 64% decrease when there was 1.0 cm of cardiac motion over conventional acquisition. In the superior-inferior (SI) direction, our protocol was able to provide between 22% and 26% decrease in the median value of the ERW across the 17 traces when there was 0.5 cm of cardiac motion and between 30% and 35% decrease when there was 1.0 cm of cardiac motion over conventional acquisition. The magnitude of the cardiac motion did not have an observable effect on the median value of the CNR. Across all 17 traces, our adaptive protocol produced noticeably more consistent, albeit lower CNR values compared with conventional acquisition. CONCLUSION: For the first time, the potential of adapting CBCT image acquisition to changes in the patient's cardiac and respiratory rates simultaneously for applications in radiotherapy was investigated. This work represents a step towards thoracic imaging that reduces the effects of both cardiac and respiratory motion on image quality.

Towards patient connected imaging with ACROBEAT: Adaptive CaRdiac cOne BEAm computed Tomography.

Robotic C-arm cone beam computed tomography (CBCT) systems are playing an increasingly pivotal role in interventional cardiac procedures and high precision radiotherapy treatments. One of the main challenges in any form of cardiac imaging is mitigating the intrinsic motion of the heart, which causes blurring and artefacts in the 3D reconstructed image. Most conventional 3D cardiac CBCT acquisition techniques attempt to combat heart motion through retrospective gating techniques, whereby acquired projections are sorted into the desired cardiac phase after the completion of the scan. However, this results in streaking artefacts and unnecessary radiation exposure to the patient. Here, we present our Adaptive CaRdiac cOne BEAm computed Tomography (ACROBEAT) acquisition protocol that uses the patient's electrocardiogram (ECG) signal to adaptively regulate the gantry velocity and projection time interval in real-time. It enables prospectively gated patient connected imaging in a single sweep of the gantry. The XCAT digital software phantom was used to complete a simulation study to compare ACROBEAT to a conventional multi-sweep retrospective ECG gated acquisition, under a variety of different acquisition conditions. The effect of location and length of the acquisition window and total number of projections acquired on image quality and total scan time were examined. Overall, ACROBEAT enables up to a 5 times average improvement in the contrast-to-noise ratio, a 40% reduction in edge response width and an 80% reduction in total projections acquired compared to conventional multi-sweep retrospective ECG gated acquisition.