C-arm orbits for metal artifact avoidance (MAA) in cone-beam CT

Effect of the cone-beam CT acquisition trajectory on image quality in spine surgery: experimental cadaver study

BACKGROUND

Intraoperative 3D imaging with cone-beam CT (CBCT) improves assessment of implant position and reduces complications in spine surgery. It is also used for image-guided surgical techniques, resulting in improved quality of care. However, in some cases, metal artifacts can reduce image quality and make it difficult to assess pedicle screw position and reduction.

PURPOSE

The objective of this study was to investigate whether a change in CBCT acquisition trajectory in relation to pedicle screw position during dorsal instrumentation can reduce metal artifacts and consequently improve image quality and clinical assessability.

STUDY DESIGN

Experimental cadaver study.

METHODS

A human cadaver was instrumented with pedicle screws in the thoracic and lumbar spine region (Th11 to L5). Then, the acquisition trajectory of the CBCT (Cios Spin, Siemens, Germany) to the pedicle screws was systematically changed in 5° steps in angulation (-30° to +30°) and swivel (-25° to +25°). Subsequently, radiological evaluation was performed by 3 blinded, qualified raters on image quality using 9 questions (including anatomical structures, implant position, appearance of artifacts) with a score (1-5 points). For statistical evaluation, the image quality of the different acquisition trajectories was compared to the standard acquisition trajectory and checked for significant differences.

RESULTS

The angulated acquisition trajectory significantly increased the score for subjective image quality (p<.001) as well as the clinical assessability of pedicle screw position (p<.001) with particularly strong effects on subjective image quality in the vertebral pedicle region (d=1.61). Swivel of the acquisition trajectory significantly improved all queried domains of subjective image quality (p<.001) as well as clinical assessability of pedicle screw position (p<.001).

CONCLUSIONS

In this cadaver study, the angulation as well as the swivel of the acquisition trajectory led to a significantly improved image quality in intraoperative 3D imaging (CBCT) with a constant isocenter. The data show that maximizing the angulation/swivel angle towards 30°/25° provides the best tested subjective image quality and enhances clinical assessability. Therefore, a correct adjustment of the acquisition trajectory can help to make intraoperative revision decisions more reliably.

CLINICAL SIGNIFICANCE

The knowledge of enhanced image quality by changing the acquisition trajectory in intraoperative 3D imaging can be used for the assessment of critical screw positions in spine surgery. The implementation of this knowledge requires only a minor change of the current intraoperative imaging workflow without additional technical equipment and could further reduce the need for revision surgery.

An interactive task-based method for the avoidance of metal artifacts in CBCT

PURPOSE

Intraoperative cone-beam CT imaging enables 3D validation of implant positioning and fracture reduction for orthopedic and trauma surgeries. However, the emergence of metal artifacts, especially in the vicinity of metallic objects, severely degrades the clinical value of the imaging modality. In previous works, metal artifact avoidance (MAA) methods have been shown to reduce metal artifacts by adapting the scanning trajectory. Yet, these methods fail to translate to clinical practice due to remaining methodological constraints and missing workflow integration.

METHODS

In this work, we propose a method to compute the spatial distribution and calibrated strengths of expected artifacts for a given tilted circular trajectory. By visualizing this as an overlay changing with the C-Arm's tilt, we enable the clinician to interactively choose an optimal trajectory while factoring in the procedural context and clinical task. We then evaluate this method in a realistic human cadaver study and compare the achieved image quality to acquisitions optimized using global metrics.

RESULTS

We assess the effectiveness of the compared methods by evaluation of image quality gradings of depicted pedicle screws. We find that both global metrics as well as the proposed visualization of artifact distribution enable a drastic improvement compared to standard non-tilted scans. Furthermore, the novel interactive visualization yields a significant improvement in subjective image quality compared to the state-of-the-art global metrics. Additionally we show that by formulating an imaging task, the proposed method allows to selectively optimize image quality and avoid artifacts in the region of interest.

CONCLUSION

We propose a method to spatially resolve predicted artifacts and provide a calibrated measure for artifact strength grading. This interactive MAA method proved practical and effective in reducing metal artifacts in the conducted cadaver study. We believe this study serves as a crucial step toward clinical application of an MAA system to improve image quality and enhance the clinical validation of implant placement.

Cone-beam CT with a noncircular (sine-on-sphere) orbit: imaging performance of a clinical system for image-guided interventions

Purpose

We aim to compare the imaging performance of a cone-beam CT (CBCT) imaging system with noncircular scan protocols (sine-on-sphere) to a conventional circular orbit.

Approach

A biplane C-arm system (ARTIS Icono; Siemens Healthineers) capable of circular and noncircular CBCT acquisition was used, with the latter orbit (sine-on-sphere, "Sine Spin") executing a sinusoidal motion with ± 10    deg tilt amplitude over the half-scan orbit. A test phantom was used for the characterization of image uniformity, noise, noise-power spectrum (NPS), spatial resolution [modulation transfer function (MTF) in axial and oblique directions], and cone-beam artifacts. Findings were interpreted using an anthropomorphic head phantom with respect to pertinent tasks in skull base neurosurgery.

Results

The noncircular scan protocol exhibited several advantages associated with improved 3D sampling-evident in the NPS as filling of the null cone about thef z spatial frequency axis and reduction of cone-beam artifacts. The region of support at the longitudinal extrema was reduced from 16 to ∼ 12    cm at a radial distance of 6.5 cm. Circular and noncircular orbits exhibited nearly identical image uniformity and quantum noise, demonstrating cupping of - 16.7 % and overall noise of ∼ 27    HU . Although both the radially averaged axial MTF (f x , y ) and 45 deg oblique MTF (f x , y , z ) were ∼ 20 % lower for the noncircular orbit compared with the circular orbit at the default full reconstruction field of view (FOV), there was no difference in spatial resolution for the medium reconstruction FOV (smaller voxel size). Differences in the perceptual image quality for the anthropomorphic phantom reinforced the objective, quantitative findings, including reduced beam-hardening and cone-beam artifacts about structures of interest in the skull base.

Conclusions

Image quality differences between circular and noncircular CBCT orbits were quantitatively evaluated on a clinical system in the context of neurosurgery. The primary performance advantage for the noncircular orbit was the improved sampling and elimination of cone-beam artifacts.

Purse-string suture with nylon cords and metal clips for the treatment of duodenal fistulae under the endoscope: a case report

Purse-string suture with nylon cords and metal clips under the endoscope is a novel therapeutic technique which is minimally invasive and it is particularly indicated for the closure and repair of gastrointestinal fistula or perforations such as duodenal fistulae. Duodenal fistulae are often caused by medical manipulation, disease progression or trauma. Once this occurs, it leads to a series of pathophysiologic changes and a variety of complications. In most cases, these complications will exacerbate the damage to the organism, and the complications are difficult to treat and can lead to infections, nutrient loss, multi-organ dysfunction and many other adverse effects. In this case report, the use of endoscopic nylon cords combined with purse-string suture and metal clips in the treatment of duodenal fistula is presented and discussed. The patient was treated with endoscopic purse-string suture and the duodenal fistula was significantly improved. The results indicate that endoscopic purse-string suture is an effective strategy for the treatment of duodenal fistulae.

Deformable registration of preoperative MR and intraoperative long-length tomosynthesis images for guidance of spine surgery via image synthesis

PURPOSE

Improved integration and use of preoperative imaging during surgery hold significant potential for enhancing treatment planning and instrument guidance through surgical navigation. Despite its prevalent use in diagnostic settings, MR imaging is rarely used for navigation in spine surgery. This study aims to leverage MR imaging for intraoperative visualization of spine anatomy, particularly in cases where CT imaging is unavailable or when minimizing radiation exposure is essential, such as in pediatric surgery.

METHODS

This work presents a method for deformable 3D-2D registration of preoperative MR images with a novel intraoperative long-length tomosynthesis imaging modality (viz., Long-Film [LF]). A conditional generative adversarial network is used to translate MR images to an intermediate bone image suitable for registration, followed by a model-based 3D-2D registration algorithm to deformably map the synthesized images to LF images. The algorithm's performance was evaluated on cadaveric specimens with implanted markers and controlled deformation, and in clinical images of patients undergoing spine surgery as part of a large-scale clinical study on LF imaging.

RESULTS

The proposed method yielded a median 2D projection distance error of 2.0 mm (interquartile range [IQR]: 1.1-3.3 mm) and a 3D target registration error of 1.5 mm (IQR: 0.8-2.1 mm) in cadaver studies. Notably, the multi-scale approach exhibited significantly higher accuracy compared to rigid solutions and effectively managed the challenges posed by piecewise rigid spine deformation. The robustness and consistency of the method were evaluated on clinical images, yielding no outliers on vertebrae without surgical instrumentation and 3% outliers on vertebrae with instrumentation.

CONCLUSIONS

This work constitutes the first reported approach for deformable MR to LF registration based on deep image synthesis. The proposed framework provides access to the preoperative annotations and planning information during surgery and enables surgical navigation within the context of MR images and/or dual-plane LF images.

Marker-based C-arm self-calibration with unknown calibration pattern

BACKGROUND

Accurate tomographic reconstructions require the knowledge of the actual acquisition geometry. Many mobile C-arm CT scanners have poorly reproducible acquisition geometries and thus need acquisition-specific calibration procedures. Most of geometric self-calibration methods based on projection data either need prior information or are limited to the estimation of a low number of geometric calibration parameters. Other self-calibration methods generally use a calibration pattern with known geometry and are hardly implementable in practice for clinical applications.

PURPOSE

We present a three-step marker based self-calibration method which does not require the prior knowledge of the calibration pattern and thus enables the use of calibration patterns with arbitrary markers positions.

METHODS

The first step of the method aims at detecting the set of markers of the calibration pattern in each projection of the CT scan and is performed using the YOLO (You Only Look Once) Convolutional Neural Network. The projected marker trajectories are then estimated by a sequential projection-wise marker association scheme based on the Linear Assignment Problem which uses Kalman filters to predict the markers 2D positions in the projections. The acquisition geometry is finally estimated from the marker trajectories using the Bundle-adjustment algorithm.

RESULTS

The calibration method has been tested on realistic simulated images of the ICRP (International Commission on Radiological Protection) phantom, using calibration patterns with 10 and 20 markers. The backprojection error was used to evaluate the self-calibration method and exhibited sub-millimeter errors. Real images of two human knees with 10 and 30 markers calibration patterns were then used to perform a qualitative evaluation of the method, which showed a remarkable artifacts reduction and bone structures visibility improvement.

CONCLUSIONS

The proposed calibration method gave promising results that pave the way to patient-specific geometric self-calibrations in clinics.

Simulation-driven training of vision transformers enables metal artifact reduction of highly truncated CBCT scans

BACKGROUND

Due to the high attenuation of metals, severe artifacts occur in cone beam computed tomography (CBCT). The metal segmentation in CBCT projections usually serves as a prerequisite for metal artifact reduction (MAR) algorithms.

PURPOSE

The occurrence of truncation caused by the limited detector size leads to the incomplete acquisition of metal masks from the threshold-based method in CBCT volume. Therefore, segmenting metal directly in CBCT projections is pursued in this work.

METHODS

Since the generation of high quality clinical training data is a constant challenge, this study proposes to generate simulated digital radiographs (data I) based on real CT data combined with self-designed computer aided design (CAD) implants. In addition to the simulated projections generated from 3D volumes, 2D x-ray images combined with projections of implants serve as the complementary data set (data II) to improve the network performance. In this work, SwinConvUNet consisting of shift window (Swin) vision transformers (ViTs) with patch merging as encoder is proposed for metal segmentation.

RESULTS

The model's performance is evaluated on accurately labeled test datasets obtained from cadaver scans as well as the unlabeled clinical projections. When trained on the data I only, the convolutional neural network (CNN) encoder-based networks UNet and TransUNet achieve only limited performance on the cadaver test data, with an average dice score of 0.821 and 0.850. After using both data II and data I during training, the average dice scores for the two models increase to 0.906 and 0.919, respectively. By replacing the CNN encoder with Swin transformer, the proposed SwinConvUNet reaches an average dice score of 0.933 for cadaver projections when only trained on the data I. Furthermore, SwinConvUNet has the largest average dice score of 0.953 for cadaver projections when trained on the combined data set.

CONCLUSIONS

Our experiments quantitatively demonstrate the effectiveness of the combination of the projections simulated under two pathways for network training. Besides, the proposed SwinConvUNet trained on the simulated projections performs state-of-the-art, robust metal segmentation as demonstrated on experiments on cadaver and clinical data sets. With the accurate segmentations from the proposed model, MAR can be conducted even for highly truncated CBCT scans.

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.

Reduction of cone-beam CT artifacts in a robotic CBCT device using saddle trajectories with integrated infrared tracking

BACKGROUND

Cone beam computed tomography (CBCT) is widely used in many medical fields. However, conventional CBCT circular scans suffer from cone beam (CB) artifacts that limit the quality and reliability of the reconstructed images due to incomplete data.

PURPOSE

Saddle trajectories in theory might be able to improve the CBCT image quality by providing a larger region with complete data. Therefore, we investigated the feasibility and performance of saddle trajectory CBCT scans and compared them to circular trajectory scans.

METHODS

We performed circular and saddle trajectory scans using a novel robotic CBCT scanner (Mobile ImagingRing (IRm); medPhoton, Salzburg, Austria). For the saddle trajectory, the gantry executed yaw motion up to± 10 ∘ $\pm 10^{\circ }$ using motorized wheels driving on the floor. An infrared (IR) tracking device with reflective markers was used for online geometric calibration correction (mainly floor unevenness). All images were reconstructed using penalized least-squares minimization with the conjugate gradient algorithm from RTK with0.5 × 0.5 × 0.5 mm 3 $0.5 \times 0.5\times 0.5 \text{ mm}^3$ voxel size. A disk phantom and an Alderson phantom were scanned to assess the image quality. Results were correlated with the local incompleteness value represented bytan ( ψ ) $\tan (\psi)$ , which was calculated at each voxel as a function of the source trajectory and the voxel's 3D coordinates. We assessed the magnitude of CB artifacts using the full width half maximum (FWHM) of each disk profile in the axial center of the reconstructed images. Spatial resolution was also quantified by the modulation transfer function at 10% (MTF10).

RESULTS

When using the saddle trajectory, the region without CB artifacts was increased from 43 to 190 mm in the SI direction compared to the circular trajectory. This region coincided with low values fortan ( ψ ) $\tan (\psi)$ . Whentan ( ψ ) $\tan (\psi)$ was larger than 0.02, we found there was a linear relationship between the FWHM andtan ( ψ ) $\tan (\psi)$ . For the saddle, IR tracking allowed the increase of MTF10 from 0.37 to 0.98 lp/mm.

CONCLUSIONS

We achieved saddle trajectory CBCT scans with a novel CBCT system combined with IR tracking. The results show that the saddle trajectory provides a larger region with reliable reconstruction compared to the circular trajectory. The proposed method can be used to evaluate other non-circular trajectories.

Effect of changing the acquisition trajectory of the 3D C-arm (CBCT) on image quality in spine surgery: experimental study using an artificial bone model

BACKGROUND

Intraoperative 3D imaging using cone-beam CT (CBCT) provides improved assessment of implant position and reduction in spine surgery, is used for navigated surgical techniques, and therefore leads to improved quality of care. However, in some cases the image quality is not sufficient to correctly assess pedicle screw position and reduction, especially due to metal artifacts. The aim of this study was to investigate whether changing the acquisition trajectory of the CBCT in relation to the pedicle screw position during dorsal instrumentation of the spine can reduce metal artifacts and consequently improve image quality as well as clinical assessability on the artificial bone model.

METHODS

An artificial bone model was instrumented with pedicle screws in the thoracic and lumbar spine region (Th10 to L5). Then, the acquisition trajectory of the CBCT (Cios Spin, Siemens, Germany) to the pedicle screws was systematically changed in 5° steps in angulation (- 30° to + 30°) and swivel (- 30° to + 30°). Subsequently, radiological evaluation was performed by three blinded, qualified raters on image quality using 9 questions (including anatomical structures, implant position, appearance of artifacts) with a score (1-5 points). For statistical evaluation, the image quality of the different acquisition trajectories was compared to the standard acquisition trajectory and checked for significant differences.

RESULTS

The angulated acquisition trajectory increased the score for subjective image quality (p < 0.001) as well as the clinical assessability of pedicle screw position (p < 0.001) highly significant with particularly strong effects on subjective image quality in the vertebral pedicle region (d = 1.06). Swivel of the acquisition trajectory significantly improved all queried domains of subjective image quality (p < 0.001) as well as clinical assessability of pedicle screw position (p < 0.001). The data show that maximizing the angulation or swivel angle toward 30° provides the best tested subjective image quality. Angulation and swivel of the acquisition trajectory result in a clinically relevant improvement in image quality in intraoperative 3D imaging (CBCT) during dorsal instrumentation of the spine.

Cone-beam CT sampling incompleteness: analytical and empirical studies of emerging systems and source-detector orbits

Purpose

Motivated by emerging cone-beam computed tomography (CBCT) systems and scan orbits, we aim to quantitatively assess the completeness of data for 3D image reconstruction-in turn, related to "cone-beam artifacts." Fundamental principles of cone-beam sampling incompleteness are considered with respect to an analytical figure-of-merit [FOM, denoted tan ( ψ min ) ] and related to an empirical FOM (denoted z mod ) for measurement of cone-beam artifact magnitude in a test phantom.

Approach

A previously proposed analytical FOM [tan ( ψ min ) , defined as the minimum angle between a point in the 3D image reconstruction and the x-ray source over the scan orbit] was analyzed for a variety of CBCT geometries. A physical test phantom was configured with parallel disk pairs (perpendicular to the z -axis) at various locations throughout the field of view, quantifying cone-beam artifact magnitude in terms of z mod (the relative signal modulation between the disks). Two CBCT systems were considered: an interventional C-arm (Cios Spin 3D; Siemens Healthineers, Forcheim Germany) and a musculoskeletal extremity scanner; Onsight3D, Carestream Health, Rochester, United States)]. Simulations and physical experiments were conducted for various source-detector orbits: (a) a conventional 360 deg circular orbit, (b) tilted and untilted semi-circular (196 deg) orbits, (c) multi-source (three x-ray sources distributed along the z axis) semi-circular orbits, and (d) a non-circular (sine-on-sphere, SoS) orbit. The incompleteness of sampling [tan ( ψ min ) ] and magnitude of cone-beam artifacts (z mod ) were evaluated for each system and orbit.

Results

The results show visually and quantitatively the effect of system geometry and scan orbit on cone-beam sampling effects, demonstrating the relationship between analytical tan ( ψ min ) and empirical z mod . Advanced source-detector orbits (e.g., three-source and SoS orbits) exhibited superior sampling completeness as quantified by both the analytical and the empirical FOMs. The test phantom and z mod metric were sensitive to variations in CBCT system geometry and scan orbit and provided a surrogate measure of underlying sampling completeness.

Conclusion

For a given system geometry and source-detector orbit, cone-beam sampling completeness can be quantified analytically (in terms arising from Tuy's condition) and/or empirically (using a test phantom for quantification of cone-beam artifacts). Such analysis provides theoretical and practical insight on sampling effects and the completeness of data for emerging CBCT systems and scan trajectories.

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

BACKGROUND

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Computed Tomography: State-of-the-Art Advancements in Musculoskeletal Imaging

ABSTRACT

Although musculoskeletal magnetic resonance imaging (MRI) plays a dominant role in characterizing abnormalities, novel computed tomography (CT) techniques have found an emerging niche in several scenarios such as trauma, gout, and the characterization of pathologic biomechanical states during motion and weight-bearing. Recent developments and advancements in the field of musculoskeletal CT include 4-dimensional, cone-beam (CB), and dual-energy (DE) CT. Four-dimensional CT has the potential to quantify biomechanical derangements of peripheral joints in different joint positions to diagnose and characterize patellofemoral instability, scapholunate ligamentous injuries, and syndesmotic injuries. Cone-beam CT provides an opportunity to image peripheral joints during weight-bearing, augmenting the diagnosis and characterization of disease processes. Emerging CBCT technologies improved spatial resolution for osseous microstructures in the quantitative analysis of osteoarthritis-related subchondral bone changes, trauma, and fracture healing. Dual-energy CT-based material decomposition visualizes and quantifies monosodium urate crystals in gout, bone marrow edema in traumatic and nontraumatic fractures, and neoplastic disease. Recently, DE techniques have been applied to CBCT, contributing to increased image quality in contrast-enhanced arthrography, bone densitometry, and bone marrow imaging. This review describes 4-dimensional CT, CBCT, and DECT advances, current logistical limitations, and prospects for each technique.

Source-detector trajectory optimization in cone-beam computed tomography: a comprehensive review on today’s state-of-the-art

Cone-beam computed tomography (CBCT) imaging is becoming increasingly important for a wide range of applications such as image-guided surgery, image-guided radiation therapy as well as diagnostic imaging such as breast and orthopaedic imaging. The potential benefits of non-circular source-detector trajectories was recognized in early work to improve the completeness of CBCT sampling and extend the field of view (FOV). Another important feature of interventional imaging is that prior knowledge of patient anatomy such as a preoperative CBCT or prior CT is commonly available. This provides the opportunity to integrate such prior information into the image acquisition process by customized CBCT source-detector trajectories. Such customized trajectories can be designed in order to optimize task-specific imaging performance, providing intervention or patient-specific imaging settings. The recently developed robotic CBCT C-arms as well as novel multi-source CBCT imaging systems with additional degrees of freedom provide the possibility to largely expand the scanning geometries beyond the conventional circular source-detector trajectory. This recent development has inspired the research community to innovate enhanced image quality by modifying image geometry, as opposed to hardware or algorithms. The recently proposed techniques in this field facilitate image quality improvement, FOV extension, radiation dose reduction, metal artifact reduction as well as 3D imaging under kinematic constraints. Because of the great practical value and the increasing importance of CBCT imaging in image-guided therapy for clinical and preclinical applications as well as in industry, this paper focuses on the review and discussion of the available literature in the CBCT trajectory optimization field. To the best of our knowledge, this paper is the first study that provides an exhaustive literature review regarding customized CBCT algorithms and tries to update the community with the clarification of in-depth information on the current progress and future trends.

Accuracy Assessment of Percutaneous Pedicle Screw Placement Using Cone Beam Computed Tomography with Metal Artifact Reduction

Metal artifact reduction (MAR) algorithms are used with cone beam computed tomography (CBCT) during augmented reality surgical navigation for minimally invasive pedicle screw instrumentation. The aim of this study was to assess intra- and inter-observer reliability of pedicle screw placement and to compare the perception of baseline image quality (NoMAR) with optimized image quality (MAR). CBCT images of 24 patients operated on for degenerative spondylolisthesis using minimally invasive lumbar fusion were analyzed retrospectively. Images were treated using NoMAR and MAR by an engineer, thus creating 48 randomized files, which were then independently analyzed by 3 spine surgeons and 3 radiologists. The Gertzbein and Robins classification was used for screw accuracy rating, and an image quality scale rated the clarity of pedicle screw and bony landmark depiction. Intra-class correlation coefficients (ICC) were calculated. NoMAR and MAR led to similarly good intra-observer (ICC > 0.6) and excellent inter-observer (ICC > 0.8) assessment reliability of pedicle screw placement accuracy. The image quality scale showed more variability in individual image perception between spine surgeons and radiologists (ICC range 0.51−0.91). This study indicates that intraoperative screw positioning can be reliably assessed on CBCT for augmented reality surgical navigation when using optimized image quality. Subjective image quality was rated slightly superior for MAR compared to NoMAR.

Universal orbit design for metal artifact elimination

Objective.Metal artifacts are a persistent problem in CT and cone-beam CT. In this work, we propose to reduce or even eliminate metal artifacts by providing better sampling of data using non-circular orbits.Approach.We treat any measurements intersecting metal as missing data, and aim to design a universal orbit that can generally accommodate arbitrary metal shapes and locations. We adapted a local sampling completeness metric based on Tuy's condition to quantify the extent of sampling in the presence of metal. A maxi-min objective over all possible metal locations was used for orbit design. A simple class of sinusoidal orbits was evaluated as a function of frequencies, maximum tilt angles, and orbital extents. Experimental implementation of these orbits were performed on an imaging bench and evaluated on two phantoms, one containing metal balls and the other containing a pedicle screw assembly for spine fixation. Metal artifact reduction (MAR) performance was compared amongst three approaches: non-circular orbits only, algorithmic correction only, and a combined approach.Main results.Theoretical evaluations of the objective favor sinusoidal orbits with large tilt angles and large orbital extents. Furthermore, orbits that leverage redundant azimuthal angles to sample non-redundant data have better performance, e.g. even or non-integer frequency sinusoids for a 360° acquisition. Experimental data support the trends observed in theoretical evaluations. Reconstructions using even or non-integer frequency orbits present less streaking artifacts and background details with finer resolution, even when multiple metal objects are present and even in the absence of MAR algorithms. The combined approach of non-circular orbits and MAR algorithm yields the best performance. The observed trend in image quality is supported by quantitative measures of sampling and severity of streaking artifact.Significance.This work demonstrates that sinusoidal orbits are generally robust against metal artifacts and can provide an avenue for improved image quality in interventional imaging.

Theory, method, and test tools for determination of 3D MTF characteristics in cone-beam CT

PURPOSE

The modulation transfer function (MTF) is widely used as an objective metric of spatial resolution of medical imaging systems. Despite advances in capability for three-dimensional (3D) isotropic spatial resolution in computed tomography (CT) and cone-beam CT (CBCT), MTF evaluation for such systems is typically reported only in the axial plane, and practical methodology for assessment of fully 3D spatial resolution characteristics is lacking. This work reviews fundamental theoretical relationships of two-dimensional (2D) and 3D spread functions and reports practical methods and test tools for analysis of 3D MTF in CBCT.

METHODS

Fundamental aspects of 2D and 3D MTF measurement are reviewed within a common notational framework, and three MTF test tools with analysis code are reported and made available online (https://istar.jhu.edu/downloads/): (a) a multi-wire tool for measurement of the axial plane MTF [denoted as M T F ( f r ; φ = 0 ∘ ) , where φ is the measurement angle out of the axial plane] as a function of position in the axial plane; (b) a wedge tool for measurement of the MTF in any direction in the 3D Fourier domain [e.g.,  φ  = 45°, denoted as M T F ( f r ; φ = 45 ∘ ) ]; and (c) a sphere tool for measurement of the MTF in any or all directions in the 3D Fourier domain. Experiments were performed on a mobile C-arm with CBCT capability, showing that M T F ( f r ; φ = 45 ∘ ) yields an informative one-dimensional (1D) representation of the overall 3D spatial resolution characteristics, capturing important characteristics of the 3D MTF that might be missed in conventional analysis. The effects of anisotropic filters and detector readout mode were investigated, and the extent to which a system can be said to provide "isotropic" resolution was evaluated by quantitative comparison of MTF at various φ .

RESULTS

All three test tools provided consistent measurement of M T F ( f r ; φ = 0 ∘ ) , and the wedge and sphere tools demonstrated how measurement of the MTF in directions outside the axial plane ( φ > 0 ∘ ) can reveal spatial resolution characteristics to which conventional axial MTF measurement is blind. The wedge tool was shown to reduce statistical measurement error compared to the sphere tool due to improved sampling, and the sphere tool was shown to provide a basis for measurement of the MTF in any or all directions (outside the null cone) from a single scan. The C-arm system exhibited non-isotropic spatial resolution with conventional non-isotropic 1D apodization filters (i.e., frequency cutoff filters) - which is common in CBCT - and implementation of isotropic 2D apodization yielded quantifiably isotropic MTF. Asymmetric pixel binning modes were similarly shown to impart non-isotropic effects on the 3D MTF, and the overall 3D MTF characteristics were evident in each case with a single, 1D measurement of the 1D M T F ( f r ; φ = 45 ∘ ).

CONCLUSION

Three test tools and corresponding MTF analysis methods were presented within a consistent framework for analysis of 3D spatial resolution characteristics in a manner amenable to routine, practical measurements. Experiments on a CBCT C-arm validated many intuitive aspects of 3D spatial resolution and quantified the extent to which a CBCT system may be considered to have isotropic resolution. Measurement of M T F ( f r ; φ = 45 ∘ ) provided a practical 1D measure of the underlying 3D MTF characteristics and is extensible to other CT or CBCT systems offering high, isotropic spatial resolution.

Sinogram + image domain neural network approach for metal artifact reduction in low-dose cone-beam computed tomography

Purpose: Cone-beam computed tomography (CBCT) is commonly used in the operating room to evaluate the placement of surgical implants in relation to critical anatomical structures. A particularly problematic setting, however, is the imaging of metallic implants, where strong artifacts can obscure visualization of both the implant and surrounding anatomy. Such artifacts are compounded when combined with low-dose imaging techniques such as sparse-view acquisition. Approach: This work presents a dual convolutional neural network approach, one operating in the sinogram domain and one in the reconstructed image domain, that is specifically designed for the physics and setting of intraoperative CBCT to address the sources of beam hardening and sparse view sampling that contribute to metal artifacts. The networks were trained with images from cadaver scans with simulated metal hardware. Results: The trained networks were tested on images of cadavers with surgically implanted metal hardware, and performance was compared with a method operating in the image domain alone. While both methods removed most image artifacts, superior performance was observed for the dual-convolutional neural network (CNN) approach in which beam-hardening and view sampling effects were addressed in both the sinogram and image domain. Conclusion: The work demonstrates an innovative approach for eliminating metal and sparsity artifacts in CBCT using a dual-CNN framework which does not require a metal segmentation.

Long-length tomosynthesis and 3D-2D registration for intraoperative assessment of spine instrumentation

PURPOSE

A system for long-length intraoperative imaging is reported based on longitudinal motion of an O-arm gantry featuring a multi-slot collimator. We assess the utility of long-length tomosynthesis and the geometric accuracy of 3D image registration for surgical guidance and evaluation of long spinal constructs.

METHODS

A multi-slot collimator with tilted apertures was integrated into an O-arm system for long-length imaging. The multi-slot projective geometry leads to slight view disparity in both long-length projection images (referred to as 'line scans') and tomosynthesis 'slot reconstructions' produced using a weighted-backprojection method. The radiation dose for long-length imaging was measured, and the utility of long-length, intraoperative tomosynthesis was evaluated in phantom and cadaver studies. Leveraging the depth resolution provided by parallax views, an algorithm for 3D-2D registration of the patient and surgical devices was adapted for registration with line scans and slot reconstructions. Registration performance using single-plane or dual-plane long-length images was evaluated and compared to registration accuracy achieved using standard dual-plane radiographs.

RESULTS

Longitudinal coverage of ∼50-64 cm was achieved with a single long-length slot scan, providing a field-of-view (FOV) up to (40 × 64) cm², depending on patient positioning. The dose-area product (reference point air kerma × x-ray field area) for a slot scan ranged from ∼702-1757 mGy·cm², equivalent to ∼2.5 s of fluoroscopy and comparable to other long-length imaging systems. Long-length scanning produced high-resolution tomosynthesis reconstructions, covering ∼12-16 vertebral levels. 3D image registration using dual-plane slot reconstructions achieved median target registration error (TRE) of 1.2 mm and 0.6° in cadaver studies, outperforming registration to dual-plane line scans (TRE = 2.8 mm and 2.2°) and radiographs (TRE = 2.5 mm and 1.1°). 3D registration using single-plane slot reconstructions leveraged the ∼7-14° angular separation between slots to achieve median TRE ∼2 mm and <2° from a single scan.

CONCLUSION

The multi-slot configuration provided intraoperative visualization of long spine segments, facilitating target localization, assessment of global spinal alignment, and evaluation of long surgical constructs. 3D-2D registration to long-length tomosynthesis reconstructions yielded a promising means of guidance and verification with accuracy exceeding that of 3D-2D registration to conventional radiographs.