Cone-beam Computed Tomography System Research Articles

Objective image quality assurance in cone-beam CT: Test methods, analysis, and workflow in longitudinal studies.

Standards for image quality evaluation in multi-detector CT (MDCT) and cone-beam CT (CBCT) are evolving to keep pace with technological advances. A clear need is emerging for methods that facilitate rigorous quality assurance (QA) with up-to-date metrology and streamlined workflow suitable to a range of MDCT and CBCT systems. To evaluate the feasibility and workflow associated with image quality (IQ) assessment in longitudinal studies for MDCT and CBCT with a single test phantom and semiautomated analysis of objective, quantitative IQ metrology. A test phantom (CorgiTM Phantom, The Phantom Lab, Greenwich, New York, USA) was used in monthly IQ testing over the course of 1 year for three MDCT scanners (one of which presented helical and volumetric scan modes) and four CBCT scanners. Semiautomated software analyzed image uniformity, linearity, contrast, noise, contrast-to-noise ratio (CNR), 3D noise-power spectrum (NPS), modulation transfer function (MTF) in axial and oblique directions, and cone-beam artifact magnitude. The workflow was evaluated using methods adapted from systems/industrial engineering, including value stream process modeling (VSPM), standard work layout (SWL), and standard work control charts (SWCT) to quantify and optimize test methodology in routine practice. The completeness and consistency of DICOM data from each system was also evaluated. Quantitative IQ metrology provided valuable insight in longitudinal quality assurance (QA), with metrics such as NPS and MTF providing insight on root cause for various forms of system failure-for example, detector calibration and geometric calibration. Monthly constancy testing showed variations in IQ test metrics owing to system performance as well as phantom setup and provided initial estimates of upper and lower control limits appropriate to QA action levels. Rigorous evaluation of QA workflow identified methods to reduce total cycle time to ∼10min for each system-viz., use of a single phantom configuration appropriate to all scanners and Head or Body scan protocols. Numerous gaps in the completeness and consistency of DICOM data were observed for CBCT systems. An IQ phantom and test methodology was found to be suitable to QA of MDCT and CBCT systems with streamlined workflow appropriate to busy clinical settings.

Super Resolution Dual-Energy Cone-Beam CT Imaging With Dual-Layer Flat-Panel Detector.

In flat-panel detector (FPD) based cone-beam computed tomography (CBCT) imaging, the native receptor array is usually binned into a smaller matrix size. By doing so, the signal readout speed could be increased by 4-9 times at the expense of a spatial resolution loss of 50%-67%. Clearly, such manipulation poses a key bottleneck in generating high spatial and high temporal resolution CBCT images at the same time. In addition, the conventional FPD is also difficult in generating dual-energy CBCT images. In this paper, we propose an innovative super resolution dual-energy CBCT imaging method, named as suRi, based on dual-layer FPD (DL-FPD) to overcome these aforementioned difficulties at once. With suRi, specifically, a 1D or 2D sub-pixel (half pixel in this study) shifted binning is applied instead of the conventionally aligned binning to double the spatial sampling rate during the dual-energy data acquisition. As a result, the suRi approach provides a new strategy to enable high spatial resolution CBCT imaging while at high readout speed. Moreover, a penalized likelihood material decomposition algorithm is developed to directly reconstruct the high resolution bases from these dual-energy CBCT projections containing sub-pixel shifts. Numerical and physical experiments are performed to validate this newly developed suRi method with phantoms and biological specimen. Results demonstrate that suRi can significantly improve the spatial resolution of the CBCT image. We believe this developed suRi method would greatly enhance the imaging performance of the DL-FPD based dual-energy CBCT systems in future.

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

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. 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 trajectoryscans. 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 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 with 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 by , 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). When using the saddle trajectory, the region without CB artifacts was increased from 43 to 190mm in the SI direction compared to the circular trajectory. This region coincided with low values for . When was larger than 0.02, we found there was a linear relationship between the FWHM and . For the saddle, IR tracking allowed the increase of MTF10 from 0.37 to 0.98 lp/mm. 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-circulartrajectories.

High-density objects in exomass affect the volume of high-density objects inside the field of view.

To evaluate the effect of the presence and the number of high-density objects in the exomass on the volume of a high-density object in cone-beam CT (CBCT). Cylinders of cobalt-chromium (Co-Cr), titanium (Ti), and zirconium (Zi) were inserted into a polymethylmethacrylate phantom in five different combinations of number and position: 1-no cylinder; 2-one cylinder in a posterior region; 3-one cylinder in an anterior region; 4-two cylinders in posterior regions; and 5-three cylinders in anterior and posterior regions. The phantom underwent CBCT scanning using OP300 and X800 systems, with the afore mentioned cylinders of the same composition placed in the exomass and an additional high-density cylinder placed in the centre of the field of view (FOV), corresponding to the left-anterior region. The tomographic volume of the cylinder inside the FOV was measured using semi-automatic segmentation. The volumetric alteration (VA) between the segmented and physical volumes, in percentage, was compared among the experimental groups using repeated measures ANOVA and Tukey post-hoc (α = 5%). The factors material, combination, and their interaction affected the volume or both CBCT systems. In OP300, more cylinders in the exomass reduced the VA, mainly for Co-Cr. In X800, more cylinders in the exomass tended to increase the VA inside the FOV, except for Zi. In general, the presence of high-density objects in the exomass influences the VA of the object inside the FOV, although this oscillates according to object composition, number and position in the exomass, and CBCT system.

Joint denoising and interpolating network for low-dose cone-beam CT reconstruction under hybrid dose-reduction strategy

Cone-beam computed tomography (CBCT) is generally reconstructed with hundreds of two-dimensional X-Ray projections through the FDK algorithm, and its excessive ionizing radiation of X-Ray may impair patients’ health. Two common dose-reduction strategies are to either lower the intensity of X-Ray, i.e., low-intensity CBCT, or reduce the number of projections, i.e., sparse-view CBCT. Existing efforts improve the low-dose CBCT images only under a single dose-reduction strategy. In this paper, we argue that applying the two strategies simultaneously can reduce dose in a gentle manner and avoid the extreme degradation of the projection data in a single dose-reduction strategy, especially under ultra-low-dose situations. Therefore, we develop a Joint Denoising and Interpolating Network (JDINet) in projection domain to improve the CBCT quality with the hybrid low-intensity and sparse-view projections. Specifically, JDINet mainly includes two important components, i.e., denoising module and interpolating module, to respectively suppress the noise caused by the low-intensity strategy and interpolate the missing projections caused by the sparse-view strategy. Because FDK actually utilizes the projection information after ramp-filtering, we develop a filtered structural similarity constraint to help JDINet focus on the reconstruction-required information. Afterward, we employ a Postprocessing Network (PostNet) in the reconstruction domain to refine the CBCT images that are reconstructed with denoised and interpolated projections. In general, a complete CBCT reconstruction framework is built with JDINet, FDK, and PostNet. Experiments demonstrate that our framework decreases RMSE by approximately 8 %, 15 %, and 17 %, respectively, on the 1/8, 1/16, and 1/32 dose data, compared to the latest methods. In conclusion, our learning-based framework can be deeply imbedded into the CBCT systems to promote the development of CBCT. Source code is available at https://github.com/LianyingChao/FusionLowDoseCBCT.

Full 3-D modulation transfer function estimation of tomosynthesis system using modified Richardson-Lucy deconvolution.

Tomosynthesis systems are three-dimensional (3-D) medical imaging devices that operate over limited acquisition angles using low radiation dosages. To measure the spatial resolution performance of a tomosynthesis system, the modulation transfer function (MTF) is widely used as a quantitative evaluationmetric. We previously introduced a method to estimate the full 3-D MTF of a cone-beam computed tomography system using two-dimensional (2-D) Richardson-Lucy (RL) deconvolution with Tikhonov-Miller regularization. However, this method can not be applied directly to estimate the 3-D MTF of a tomosynthesis system, since the unique artifacts (i.e., shadow artifacts, spreading tails, directional blurring, and high-level noise) of the system produce several errors that lower the estimation performance. Varying positions of the negative pixels due to shadow artifacts and spreading tails cause inconsistent deconvolution performances at each of the directional projections, and the severe noise in the reconstructed images cause noise amplification during estimation. This work proposes several modifications to the previous method to resolve the inconsistent performance and noise amplification errors to increase the full 3-D MTF estimationaccuracy. Three modifications were introduced to the 2-D RL deconvolution to prevent estimation errors and improve MTF estimation performance: non-negativity relaxation function, cost function to terminate the iterative process of RL deconvolution, and regularization strength for noise control. To validate the effectiveness of the proposed modifications, we reconstructed sphere phantoms from simulation and experimental tomosynthesis studies in the iso-center and offset-center positions as well as estimated the full 3-D MTFs using the previous and proposed methods. We compared the 3-D render images, central plane images, and center profiles of the estimated 3-D MTFs and calculated the full widths at half and tenth maximum for quantitativeevaluation. The previous method cannot estimate the full 3-D MTF of a tomosynthesis system; its inaccurate negative pixel relaxation produces circular-shaped errors, and the mean squared error based simple cost function for termination causes inconsistent estimation at each directional projection to diminish the clear edges of the low-frequency drop and missing sample regions. Noise amplification from lack of noise regularization is also seen in the previous method results. Compared to the previous method, the proposed method shows superior estimation performance at reducing errors in both the simulation and experimental studies regardless of object position. The proposed method preserves the low-frequency drop, missing sample regions from the limited acquisition angles, and missing cone region from the offset-center position; the estimated MTFs also show FWHM and FWTM values close to those of the ideal MTFs than with the previousmethod. This work presents a method to estimate the full 3-D MTF of a tomosynthesis system. The proposed modifications prevent circular-shaped errors and noise amplification due to the geometry for limited acquisition angles and high noise levels. Compared to our previous method, the proposed scheme show better performance for estimating the 3-D MTF of the tomosynthesissystem.

A simulation study to evaluate the effect of 2D antiscatter grid primary transmission on flat panel detector based CBCT image quality

Abstact Purpose. Two-dimensional antiscatter grids’ (2D-ASGs) septal shadows and their impact on primary transmission play a critical role in cone-beam computed tomography (CBCT) image noise and artifact characteristics. Therefore, a numerical simulation platform was developed to evaluate the effect of 2D-ASG’s primary transmission on image quality, as a function of grid geometry and CBCT system properties. Methods. To study the effect of 2D-ASG’s septal shadows on primary transmission and CBCT image quality, two new methods were introduced; one to simulate projection signal gradients in septal shadows, and the other to simulate septal shadow variations due to gantry flex. Signal gradients in septal shadows were simulated by generating a system point spread function that was directly extracted from projection images of 2D-ASG prototypes in experiments. Variations in septal shadows due to gantry flex were simulated by generating oversampled shadow profiles extracted from experiments. Subsequently, the effect of 2D-ASG’s septal shadows on primary transmission and image quality was evaluated. Results. For an apparent septal thickness of 0.15 mm, the primary transmission of 2D-ASG varied between 72%–90% for grid pitches 1–3 mm. In low-contrast phantoms, the effect of 2D-ASG’s radiopaque footprint on information loss was subtle. At high spatial frequencies, information loss manifested itself as undersampling artifacts, however, its impact on image quality is subtle when compared to quantum noise. Effects of additive electronic noise and gantry flex induced ring artifacts on image quality varied as a function of grid pitch and septal thickness. Such artifacts were substantially less in lower resolution images. Conclusion. The proposed simulation platform allowed successful evaluation of CBCT image quality variations as a function of 2D-ASG primary transmission properties and CBCT system characteristics. This platform can be potentially used for optimizing 2D-ASG design properties based on the imaging task and properties of the CBCT system.

Proton Online Adaptation Using Novel Cone-Beam Computed Tomography System

Cone beam computed tomography (CBCT) has been used in clinic frequently to provide quick, online, and handy three-dimensional images. However, due to its significant artifacts and inaccurate Hounsfield values (HU), CBCT based online proton adaptation is infeasible. Recently, a novel in-room imaging solution was developed to allow larger images, better contrast and faster CBCT imaging acquisition. Our work demonstrated for the first time the feasibility of using this novel CBCT images for direct intensity modulated proton planning. Three patients and three CIRS phantoms were scanned using novel imaging technique (CBCTp) for this work. CT curves were acquired by scanning a CIRS electron density phantom. Stopping power ratio was computed using the stoichiometric method. Proton plans were made on treatment simulation CT and further evaluated on CBCTp to compare the differences. The table below shows three patients and three phantoms plan comparisons. Multiple sites and dose levels were studied. The planning target volume (PTV) coverage (D95%) and mean dose difference between simulation CT and CBCTp are 0.8% and 0.3%, correspondingly. Dosimetrically, phantom and patient plans are almost identical between two imaging techniques. Lung patient plan show the largest variation due to patient tumor change and the quality of breathing. The novel and high quality of CBCT is feasible for direct proton planning with accurate CT calibration. It provides a reliable alternative solution for proton online adaptive therapy.

Implementation and evaluation of the AAPM TG 111 for radiation dosimetry in a state-of-the-art C-arm cone beam CT system

The aim of this study was to assess the dosimetric characteristics of a state-of-the-art C-arm cone beam computed tomography (CBCT) system using the methodology proposed by the American Association of Physicists in Medicine (AAPM) Task Group (TG) 111. The dose measurement methodology described in AAPM TG 111 for wide cone beam acquisitions without table translation was employed to estimate equilibrium beam length () and equilibrium dose () in various interventional task-specific protocols with different tube arc projection geometries. Dose profiles were derived from point dose measurements in the centre and peripheral locations of the ICRU/AAPM and standard polymethyl methacrylate (PMMA) body phantom. Dose measurements were performed in phantom and free-in-air using a solid-state point detector. Monte Carlo (MC) based simulation dosimetry was performed to quantify the inhomogeneous dose patterns imparted in the phantoms. Estimated and on the ICRU/AAPM phantom was up to 49.4 cm and 6.17 mGy/100 mAs, respectively. Corresponding values determined on the PMMA phantom were 139 cm and 8.8 mGy/100 mAs, respectively. Free-in-air dose measurement ranged from 1.43 mGy/100 mAs to 5.93 mGy/100 mAs. Per cent difference in and between MC simulation and solid-state point detector measurement methods in the ICRU/AAPM phantom were within 16% and 18%, respectively. Manufacturers can use the presented methodology to characterize the dosimetric properties of C-arm CBCT systems. Clinical medical physicists may follow this methodology to verify corresponding data provided by the manufacturer and check for C-arm CBCT system performance dosimetric consistency.

Comparison of intraoperative CT- and cone beam CT-based spinal navigation for the treatment of atlantoaxial instability

BACKGROUND CONTEXTDue to the complexity of neurovascular structures in the atlantoaxial region, spinal navigation for posterior C1–C2 instrumentation is nowadays a helpful tool to increase accuracy of surgery and safety of patients. Many available intraoperative navigation devices have proven their reliability in this part of the spine. Two main imaging techniques are used: intraoperative CT (iCT) and cone beam computed tomography (CBCT). PURPOSEComparison of iCT- and CBCT-based technologies for navigated posterior instrumentation in C1-C2 instability. STUDY DESIGNRetrospective study. PATIENT SAMPLEA total of 81 consecutive patients from July 2014 to April 2020. OUTCOME MEASURESScrew accuracy and operating time. METHODSPatients with C1–C2 instability received posterior instrumentation using C2 pedicle screws, C1 lateral mass or pedicle screws. All screws were inserted using intraoperative imaging either using iCT or CBCT systems and spinal navigation with autoregistration technology. Following navigated screw insertion, a second intraoperative scan was performed to assess the accuracy of screw placement. Accuracy was defined as the percentage of correctly placed screws or with minor cortical breach (<2 mm) as graded by an independent observer compared to misplaced screws. RESULTSA total of 81 patients with C1–C2 instability were retrospectively analyzed. Of these, 34 patients were operated with the use of iCT and 47 with CBCT. No significant demographic difference was found between groups. In the iCT group, 97.7% of the C1–C2 screws were correctly inserted; 2.3% showed a minor cortical breach (<2 mm); no misplacement (>2 mm). In the CBCT group, 98.9% of screws were correctly inserted; no minor pedicle breach; 1.1% showed misplacement >2 mm. Accuracy of screw placement demonstrated no significant difference between groups. Both technologies allowed sufficient identification of screw misplacement intraoperatively leading to two screw revisions in the iCT and three in the CBCT group. Median time of surgery was significantly shorter using CBCT technology (166.5 minutes [iCT] vs 122 minutes [CBCT]; p<.01). CONCLUSIONSSpinal navigation using either iCT- or CBCT-based systems with autoregistration allows safe and reliable screw placement and intraoperative assessment of screw positioning. Using the herein presented procedural protocols, CBCT systems allow shorter operating time.

A method for sensitivity analysis of deviation factor for geometric correction of cone-beam CT system

To propose a sensitivity test method for geometric correction position deviation of cone-beam CT systems. We proposed the definition of center deviation and its derivation. We analyzed the influence of the variation of the three-dimensional spatial center of the steel ball point, the projection center and the size of the steel ball point on the deviation of geometric parameters and the reconstructed image results by calculating the geometric correction parameters based on the Noo analytical method using the FDK reconstruction algorithm for image reconstruction. The radius of the steel ball point was within 3 mm. The deviation of the center of the calibration parameter was within the order of magnitude and negligible. A 10% Gaussian perturbation of a single pixel in the 3D spatial coordinates of the steel ball point produced a deviation of about 3 pixel sizes, while the same Gaussian perturbation of the 2D projection coordinates of the steel ball point produced a deviation of about 2 pixel sizes. The geometric correction is more sensitive to the deviation generated by the three-dimensional spatial coordinates of the steel ball point with limited sensitivity to the deviation generated by the two-dimensional projection coordinates of the steel ball point. The deviation sensitivity of a small diameter steel ball point can be ignored.

Evaluation of 2- and 3-dimensional anatomic parameters of C-shaped root canals with cone beam computed tomography, microcomputed tomography, and nanocomputed tomography

The objective of this study was to evaluate 2-dimensional (2D) and 3D morphometric parameters of C-shaped root canals on cone beam computed tomography (CBCT) and microcomputed tomography (microCT) images using nanocomputed tomography (nanoCT) as the reference standard. Sixty mandibular molars with C-shaped canals were individually scanned using nanoCT and microCT. Cone beam computed tomography acquisitions were then performed with 4 CBCT systems, using high and standard resolutions. The 2D parameters of roundness and major and minor diameters were obtained in the cross sections of the root canals at 1, 2, and 3 mm from the root apex. The 3D parameters of surface area, volume, and structure model index were measured considering the entire extension of the root canals. Absolute error (AE) in measurement was calculated against the nanoCT values. Data were statistically analyzed with the Shapiro-Wilk test and analysis of variance (α=0.05). No significant differences in AE were discovered for the 2D parameters among microCT and the CBCT scans. The AE values for the 3D parameters of volume and surface area were significantly smaller in microCT compared to all CBCT units (P < .05). Significantly lower AE values for surface area were observed in high resolution compared to standard resolution for all CBCT units (P < .05). Structure model index did not differ significantly among microCT and all CBCT protocols. Cone beam computed tomography images showed accuracy for evaluating 2D parameters and over- and underestimation for 3D parameters.

CBCT Evaluation of Lingual Foramen and its Anatomic Variations in Northeast Indian Population.

The identification of nearby neurovascular structures is a vital part of dental implant treatment planning. Commonly two implants are placed in the mandibular inter foraminal region to support the overdenture. Placing dental implants in this region is considered relatively safe. The purpose of this study was to employ cone beam computed tomography (CBCT) on a representative sample of Northeast Indians to find out the location and distribution of lingual foramen in the anterior mandible. A total of 104 CBCT scans were performed using the Cranex 3D CBCT system. Three Oral and Maxillofacial Radiologists agreed on a visibility grade for assessing the lingual foramen using a 3-point scale. The number and presence of foramina were noted. On cross-sectional pictures, the distance in millimeters between the crest of the alveolar ridge and the superior cortical margin of the foramen was determined. Furthermore, the analysis of the distribution of lingual foramina in terms of frequency and percentages was performed. Scan results showed that the lingual foramen was clearly visible in 65.4% of cases, somewhat visible in 12.5% of scans, and not visible at all in 22.1% of cases. Half of the people in the study had two foramina, whereas 27.9 percent had one and 22.1 percent had none. The mean distance from the crest of the alveolar ridge was 10.84 mm in males whereas 9.91 mm in females. The Northeast Indian population exhibits an unusually high prevalence of lingual foramen in the anterior portion of the mandible. This research will aid in the precise pre-operative identification of the position and number of lingual foramina prior to surgical treatments involving the midline of the mandible.

Weight-bearing gantry-free cone-beam CT of the lumbar spine: Image quality analysis and dose efficiency

PurposeThe effect of static forces under load limits the prognostic value of lumbar spine CT in a horizontal position. Using a gantry-free scanner architecture, this study was designed to assess the feasibility of weight-bearing cone-beam CT (CBCT) of the lumbar spine and to establish the most dose-effective combination of scan parameters. MethodsEight formalin-fixated cadaveric specimens were examined with a gantry-free CBCT system in upright position with the aid of a dedicated positioning backstop. Cadavers were scanned with eight combinations of tube voltage (102 or 117 kV), detector entrance dose level (high or low), and frame rates (16 or 30 fps). Five radiologists independently analyzed datasets for overall image quality and posterior wall assessability. Additionally, image noise and signal-to-noise ratio (SNR) were compared based on region-of-interest (ROI) measurements in the gluteal muscles. ResultsRadiation dose ranged from 6.8 ± 1.6 (117 kV, dose level low, 16 fps) to 24.3 ± 6.3 mGy (102 kV, dose level high, 30 fps). Both image quality and posterior wall assessability were favored with 30 over 16 fps (all p ≤ 0.008). In contrast, both tube voltage (all p > 0.999) and dose level (all p > 0.096) did not significantly impact reader assessment. Image noise decreased considerably with higher frame rates (all p ≤ 0.040), while SNR ranged from 0.56 ± 0.03 to 1.11 ± 0.30 without a significant difference between scan protocols (all p ≥ 0.060). ConclusionsEmploying an optimized scan protocol, weight-bearing gantry-free CBCT of the lumbar spine allows for diagnostic imaging at reasonable radiation dose.

Nonlinear ill-posed problem in low-dose dental cone-beam computed tomography

Abstract This paper describes the mathematical structure of the ill-posed nonlinear inverse problem of low-dose dental cone-beam computed tomography (CBCT) and explains the advantages of a deep learning-based approach to the reconstruction of computed tomography images over conventional regularization methods. This paper explains the underlying reasons why dental CBCT is more ill-posed than standard computed tomography. Despite this severe ill-posedness, the demand for dental CBCT systems is rapidly growing because of their cost competitiveness and low radiation dose. We then describe the limitations of existing methods in the accurate restoration of the morphological structures of teeth using dental CBCT data severely damaged by metal implants. We further discuss the usefulness of panoramic images generated from CBCT data for accurate tooth segmentation. We also discuss the possibility of utilizing radiation-free intra-oral scan data as prior information in CBCT image reconstruction to compensate for the damage to data caused by metal implants.

Image quality and radiation dose in spinal surgery: a comparison of three imaging systems in phantom.

Using optimal settings for x-ray scans is crucial for obtaining three-dimensional images of high quality while keeping the patient dose low. Our work compares dose and image quality (IQ) of three intraoperative imaging systems [O-arm cone-beam computed tomography (CBCT), ClarifEye C-arm CBCT, and Airo computed tomography] used for spinal surgery. Patients of 70, 90, and 110kg were simulated with an anthropomorphic phantom by adding tissue-equivalent material. Titanium inserts were placed in the phantom spine for reproducing metal artifacts in the images. Organ dose was measured with thermo-luminescent dosimeters for effective dose () calculation. Subjective IQ was assessed by ranking the images acquired with the manufacturer-defined imaging protocols. Objective IQ was assessed with a customized Catphan phantom. The ClarifEye protocols resulted in the lowest ranging from 1.4 to 5.1mSv according to phantom size and protocol. The highest was measured for the high-definition protocol of O-arm ( 2.2 to 9mSv) providing the best subjective IQ for imaging of the spine without titanium inserts. For the images with metal, the best IQ was obtained with ClarifEye. Airo ( 5.5 to 8.4mSv) was ranked with the lowest IQ for images without metal while the rank improved for images with metal. Airo images had better uniformity, noise, and contrast sensitivity compared with CBCTs but worse high-contrast resolution. The values of these parameters were comparable between the CBCT systems. Both CBCT systems provided better IQ compared with Airo for navigation of lumbar spinal surgery for the original phantom. Metal artifacts particularly affect O-arm images decreasing the subjective IQ. The high spatial resolution of CBCT systems resulted in a relevant parameter for the visibility of anatomical features important for spine navigation. Low dose protocols were enough to obtain a clinically acceptable contrast-to-noise ratio in the bones.

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

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 ] and related to an empirical FOM (denoted ) for measurement of cone-beam artifact magnitude in a test phantom. A previously proposed analytical FOM [, 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 -axis) at various locations throughout the field of view, quantifying cone-beam artifact magnitude in terms of (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 360deg circular orbit, (b)tilted and untilted semi-circular (196deg) orbits, (c)multi-source (three x-ray sources distributed along the axis) semi-circular orbits, and (d)a non-circular (sine-on-sphere, SoS) orbit. The incompleteness of sampling [] and magnitude of cone-beam artifacts () were evaluated for each system and orbit. The results show visually and quantitatively the effect of system geometry and scan orbit on cone-beam sampling effects, demonstrating the relationship between analytical and empirical . 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 metric were sensitive to variations in CBCT system geometry and scan orbit and provided a surrogate measure of underlying sampling completeness. 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 Novel Reconstruction of the Sparse-View CBCT Algorithm for Correcting Artifacts and Reducing Noise

X-ray tomography is often affected by noise and artifacts during the reconstruction process, such as detector offset, calibration errors, metal artifacts, etc. Conventional algorithms, including FDK and SART, are unable to satisfy the sampling theorem requirements for 3D reconstruction under sparse-view constraints, exacerbating the impact of noise and artifacts. This paper proposes a novel 3D reconstruction algorithm tailored to sparse-view cone-beam computed tomography (CBCT). Drawing upon compressed sensing theory, we incorporate the weighted Schatten p-norm minimization (WSNM) algorithm for 2D image denoising and the adaptive steepest descent projection onto convex sets (ASD-POCS) algorithm, which employs a total variation (TV) regularization term. These inclusions serve to reduce noise and ameliorate artifacts. Our proposed algorithm extends the WSNM approach into three-dimensional space and integrates the ASD-POCS algorithm, enabling 3D reconstruction with digital brain phantoms, clinical medical data, and real projections from our portable CBCT system. The performance of our algorithm surpasses traditional methods when evaluated using root mean square error (RMSE), peak signal-to-noise ratio (PSNR), and structural similarity index measure (SSIM) metrics. Furthermore, our approach demonstrates marked enhancements in artifact reduction and noise suppression.

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

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. 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. 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. 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.

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.