A local shift-variant Fourier model and experimental validation of circular cone-beam computed tomography artifacts

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.

AAPM Task Group Report 238: 3D C-arms with volumetric imaging capability

This report reviews the image acquisition and reconstruction characteristics of C-arm Cone Beam Computed Tomography (C-arm CBCT) systems and provides guidance on quality control of C-arm systems with this volumetric imaging capability. The concepts of 3D image reconstruction, geometric calibration, image quality, and dosimetry covered in this report are also pertinent to CBCT for Image-Guided Radiation Therapy (IGRT). However, IGRT systems introduce a number of additional considerations, such as geometric alignment of the imaging at treatment isocenter, which are beyond the scope of the charge to the task group and the report. Section 1 provides an introduction to C-arm CBCT systems and reviews a variety of clinical applications. Section 2 briefly presents nomenclature specific or unique to these systems. A short review of C-arm fluoroscopy quality control (QC) in relation to 3D C-arm imaging is given in Section 3. Section 4 discusses system calibration, including geometric calibration and uniformity calibration. A review of the unique approaches and challenges to 3D reconstruction of data sets acquired by C-arm CBCT systems is give in Section 5. Sections 6 and 7 go in greater depth to address the performance assessment of C-arm CBCT units. First, Section 6 describes testing approaches and phantoms that may be used to evaluate image quality (spatial resolution and image noise and artifacts) and identifies several factors that affect image quality. Section 7 describes both free-in-air and in-phantom approaches to evaluating radiation dose indices. The methodologies described for assessing image quality and radiation dose may be used for annual constancy assessment and comparisons among different systems to help medical physicists determine when a system is not operating as expected. Baseline measurements taken either at installation or after a full preventative maintenance service call can also provide valuable data to help determine whether the performance of the system is acceptable. Collecting image quality and radiation dose data on existing phantoms used for CT image quality and radiation dose assessment, or on newly developed phantoms, will inform the development of performance criteria and standards. Phantom images are also useful for identifying and evaluating artifacts. In particular, comparing baseline data with those from current phantom images can reveal the need for system calibration before image artifacts are detected in clinical practice. Examples of artifacts are provided in Sections 4, 5, and 6.

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.

Flat-panel conebeam CT in the clinic: history and current state

Research into conebeam CT concepts began as soon as the first clinical single-slice CT scanner was conceived. Early implementations of conebeam CT in the 1980s focused on high-contrast applications where concurrent high resolution ( < 200    μ m ), for visualization of small contrast-filled vessels, bones, or teeth, was an imaging requirement that could not be met by the contemporaneous CT scanners. However, the use of nonlinear imagers, e.g., x-ray image intensifiers, limited the clinical utility of the earliest diagnostic conebeam CT systems. The development of consumer-electronics large-area displays provided a technical foundation that was leveraged in the 1990s to first produce large-area digital x-ray detectors for use in radiography and then compact flat panels suitable for high-resolution and high-frame-rate conebeam CT. In this review, we show the concurrent evolution of digital flat panel (DFP) technology and clinical conebeam CT. We give a brief summary of conebeam CT reconstruction, followed by a brief review of the correction approaches for DFP-specific artifacts. The historical development and current status of flat-panel conebeam CT in four clinical areas-breast, fixed C-arm, image-guided radiation therapy, and extremity/head-is presented. Advances in DFP technology over the past two decades have led to improved visualization of high-contrast, high-resolution clinical tasks, and image quality now approaches the soft-tissue contrast resolution that is the standard in clinical CT. Future technical developments in DFPs will enable an even broader range of clinical applications; research in the arena of flat-panel CT shows no signs of slowing down.

DeepRegularizer: Rapid Resolution Enhancement of Tomographic Imaging Using Deep Learning

Optical diffraction tomography measures the three-dimensional refractive index map of a specimen and visualizes biochemical phenomena at the nanoscale in a non-destructive manner. One major drawback of optical diffraction tomography is poor axial resolution due to limited access to the three-dimensional optical transfer function. This missing cone problem has been addressed through regularization algorithms that use a priori information, such as non-negativity and sample smoothness. However, the iterative nature of these algorithms and their parameter dependency make real-time visualization impossible. In this article, we propose and experimentally demonstrate a deep neural network, which we term DeepRegularizer, that rapidly improves the resolution of a three-dimensional refractive index map. Trained with pairs of datasets (a raw refractive index tomogram and a resolution-enhanced refractive index tomogram via the iterative total variation algorithm), the three-dimensional U-net-based convolutional neural network learns a transformation between the two tomogram domains. The feasibility and generalizability of our network are demonstrated using bacterial cells and a human leukaemic cell line, and by validating the model across different samples. DeepRegularizer offers more than an order of magnitude faster regularization performance compared to the conventional iterative method. We envision that the proposed data-driven approach can bypass the high time complexity of various image reconstructions in other imaging modalities.

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.

Dynamical machine learning volumetric reconstruction of objects' interiors from limited angular views

Limited-angle tomography of an interior volume is a challenging, highly ill-posed problem with practical implications in medical and biological imaging, manufacturing, automation, and environmental and food security. Regularizing priors are necessary to reduce artifacts by improving the condition of such problems. Recently, it was shown that one effective way to learn the priors for strongly scattering yet highly structured 3D objects, e.g. layered and Manhattan, is by a static neural network [Goy et al. Proc. Natl. Acad. Sci. 116, 19848-19856 (2019)]. Here, we present a radically different approach where the collection of raw images from multiple angles is viewed analogously to a dynamical system driven by the object-dependent forward scattering operator. The sequence index in the angle of illumination plays the role of discrete time in the dynamical system analogy. Thus, the imaging problem turns into a problem of nonlinear system identification, which also suggests dynamical learning as a better fit to regularize the reconstructions. We devised a Recurrent Neural Network (RNN) architecture with a novel Separable-Convolution Gated Recurrent Unit (SC-GRU) as the fundamental building block. Through a comprehensive comparison of several quantitative metrics, we show that the dynamic method is suitable for a generic interior-volumetric reconstruction under a limited-angle scheme. We show that this approach accurately reconstructs volume interiors under two conditions: weak scattering, when the Radon transform approximation is applicable and the forward operator well defined; and strong scattering, which is nonlinear with respect to the 3D refractive index distribution and includes uncertainty in the forward operator.

Dataset of patient-derived digital breast phantoms for in silico studies in breast computed tomography, digital breast tomosynthesis, and digital mammography

PURPOSE

To present a dataset of computational digital breast phantoms derived from high-resolution three-dimensional (3D) clinical breast images for the use in virtual clinical trials in two-dimensional (2D) and 3D x-ray breast imaging.

ACQUISITION AND VALIDATION METHODS

Uncompressed computational breast phantoms for investigations in dedicated breast CT (BCT) were derived from 150 clinical 3D breast images acquired via a BCT scanner at UC Davis (California, USA). Each image voxel was classified in one out of the four main materials presented in the field of view: fibroglandular tissue, adipose tissue, skin tissue, and air. For the image classification, a semi-automatic software was developed. The semi-automatic classification was compared via manual glandular classification performed by two researchers. A total of 60 compressed computational phantoms for virtual clinical trials in digital mammography (DM) and digital breast tomosynthesis (DBT) were obtained from the corresponding uncompressed phantoms via a software algorithm simulating the compression and the elastic deformation of the breast, using the tissue's elastic coefficient. This process was evaluated in terms of glandular fraction modification introduced by the compression procedure. The generated cohort of 150 uncompressed computational breast phantoms presented a mean value of the glandular fraction by mass of 12.3%; the average diameter of the breast evaluated at the center of mass was 105 mm. Despite the slight differences between the two manual segmentations, the resulting glandular tissue segmentation did not consistently differ from that obtained via the semi-automatic classification. The difference between the glandular fraction by mass before and after the compression was 2.1% on average. The 60 compressed phantoms presented an average glandular fraction by mass of 12.1% and an average compressed thickness of 61 mm.

DATA FORMAT AND ACCESS

The generated digital breast phantoms are stored in DICOM files. Image voxels can present one out of four values representing the different classified materials: 0 for the air, 1 for the adipose tissue, 2 for the glandular tissue, and 3 for the skin tissue. The generated computational phantoms datasets were stored in the Zenodo public repository for research purposes (http://doi.org/10.5281/zenodo.4529852, http://doi.org/10.5281/zenodo.4515360).

POTENTIAL APPLICATIONS

The dataset developed within the INFN AGATA project will be used for developing a platform for virtual clinical trials in x-ray breast imaging and dosimetry. In addition, they will represent a valid support for introducing new breast models for dose estimates in 2D and 3D x-ray breast imaging and as models for manufacturing anthropomorphic physical phantoms.

Design of a Computed Tomography Automation Architecture

This paper presents a literature review on techniques related to the computed tomography procedure that incorporate automation elements in their research investigations or industrial applications. Computed tomography (CT) is a non-destructive testing (NDT) technique in that the imaging and inspection are performed without damaging the sample, allowing for additional or repeated analysis if necessary. The reviewed literature is organized based on the steps associated with a general NDT task in order to define an end-to-end computed tomography automation architecture. The process steps include activities prior to image collection, during the scan, and after the data are collected. It further reviews efforts related to repeating this process based on a previous scan result. By analyzing the multiple existing but disparate efforts found in the literature, we present a framework for fully automating NDT procedures and discuss the remaining technical gaps in the developed framework.

Accelerated 3D image reconstruction with a morphological pyramid and noise-power convergence criterion

Model-based iterative reconstruction (MBIR) for cone-beam CT (CBCT) offers better noise-resolution tradeoff and image quality than analytical methods for acquisition protocols with low x-ray dose or limited data, but with increased computational burden that poses a drawback to routine application in clinical scenarios. This work develops a comprehensive framework for acceleration of MBIR in the form of penalized weighted least squares optimized with ordered subsets separable quadratic surrogates. The optimization was scheduled on a set of stages forming a morphological pyramid varying in voxel size. Transition between stages was controlled with a convergence criterion based on the deviation between the mid-band noise power spectrum (NPS) measured on a homogeneous region of the evolving reconstruction and that expected for the converged image, computed with an analytical model that used projection data and the reconstruction parameters. A stochastic backprojector was developed by introducing a random perturbation to the sampling position of each voxel for each ray in the reconstruction within a voxel-based backprojector, breaking the deterministic pattern of sampling artifacts when combined with an unmatched Siddon forward projector. This fast, forward and backprojector pair were included into a multi-resolution reconstruction strategy to provide support for objects partially outside of the field of view. Acceleration from ordered subsets was combined with momentum accumulation stabilized with an adaptive technique that automatically resets the accumulated momentum when it diverges noticeably from the current iteration update. The framework was evaluated with CBCT data of a realistic abdomen phantom acquired on an imaging x-ray bench and with clinical CBCT data from an angiography robotic C-arm (Artis Zeego, Siemens Healthineers, Forchheim, Germany) acquired during a liver embolization procedure. Image fidelity was assessed with the structural similarity index (SSIM) computed with a converged reconstruction. The accelerated framework provided accurate reconstructions in 60 s (SSIM = 0.97) and as little as 27 s (SSIM = 0.94) for soft-tissue evaluation. The use of simple forward and backprojectors resulted in 9.3× acceleration. Accumulation of momentum provided extra ∼3.5× acceleration with stable convergence for 6-30 subsets. The NPS-driven morphological pyramid resulted in initial faster convergence, achieving similar SSIM with 1.5× lower runtime than the single-stage optimization. Acceleration of MBIR to provide reconstruction time compatible with clinical applications is feasible via architectures that integrate algorithmic acceleration with approaches to provide stable convergence, and optimization schedules that maximize convergence speed.

Differentiated Backprojection Domain Deep Learning for Conebeam Artifact Removal

Conebeam CT using a circular trajectory is quite often used for various applications due to its relative simple geometry. For conebeam geometry, Feldkamp, Davis and Kress algorithm is regarded as the standard reconstruction method, but this algorithm suffers from so-called conebeam artifacts as the cone angle increases. Various model-based iterative reconstruction methods have been developed to reduce the cone-beam artifacts, but these algorithms usually require multiple applications of computational expensive forward and backprojections. In this paper, we develop a novel deep learning approach for accurate conebeam artifact removal. In particular, our deep network, designed on the differentiated backprojection domain, performs a data-driven inversion of an ill-posed deconvolution problem associated with the Hilbert transform. The reconstruction results along the coronal and sagittal directions are then combined using a spectral blending technique to minimize the spectral leakage. Experimental results under various conditions confirmed that our method generalizes well and outperforms the existing iterative methods despite significantly reduced runtime complexity.

Dynamic focal plane estimation for dental panoramic radiography

PURPOSE

The digital panoramic radiography is widely used in dental clinics and provides the anatomical information of the intraoral structure along the predefined arc-shaped path. Since the intraoral structure varies depending on the patient, however, it is nearly impossible to design a common and static focal path or plane fitted to the dentition of all patients. In response, we introduce an imaging algorithm for digital panoramic radiography that can provide a focused panoramic radiographic image for all patients, by automatically estimating the best focal plane for each patient.

METHODS

The aim of this study is to improve the image quality of dental panoramic radiography based on a three-dimensional (3D) dynamic focal plane. The plane is newly introduced to represent the arbitrary 3D intraoral structure of each patient. The proposed algorithm consists of three steps: preprocessing, focal plane estimation, and image reconstruction. We first perform preprocessing to improve the accuracy of focal plane estimation. The 3D dynamic focal plane is then estimated by adjusting the position of the image plane so that object boundaries in the neighboring projection data are aligned or focused on the plane. Finally, a panoramic radiographic image is reconstructed using the estimated dynamic focal plane.

RESULTS

The proposed algorithm is evaluated using a numerical phantom dataset and four clinical human datasets. In order to examine the image quality improvement owing to the proposed algorithm, we generate panoramic radiographic images based on a conventional static focal plane and estimated 3D dynamic focal planes, respectively. Experimental results show that the image quality is dramatically improved for all datasets using the 3D dynamic focal planes that are estimated from the proposed algorithm.

CONCLUSIONS

We propose an imaging algorithm for digital panoramic radiography that provides improved image quality by estimating dynamic focal planes fitted to each individual patient's intraoral structure.

Non-Destructive Testing of Ceramic Knee Implants Using Micro-Computed Tomography

BACKGROUND

Because of the accumulated numbers and the increasing rate of knee replacement surgeries larger numbers of revision cases are likely. Although the success rate of knee arthroplasties is high, complications like the loosening of the implant necessitate subsequent treatments. Therefore, new concepts such as metal-free ceramic implants are necessary, for example, using zirconium dioxide (ZrO₂). Several studies showed that the strength of ceramic ZrO₂ implants is equivalent to cobalt-chromium components.

METHODS

Non-destructive testing remains challenging due to the high density (6 g/cm³) of ZrO₂. In this feasibility study, we investigated 8 tibial and 8 femoral implants respectively using an industrial X-ray micro-computed tomography (XCT) system at a voxel size of 100 μm. We established a non-destructive testing protocol for ceramic knee implants optimizing scanning parameters and sample orientation using CT simulations. Finally, we used an iterative artifact reduction procedure for beam hardening correction.

RESULTS

The results show that corrected image data enable the non-destructive inspection of high-density components. In this sample, none of the investigated components show any internal defects like pores or cracks. In general, XCT is a major imaging method that is able to provide a 3-dimensional representation of higher dense objects that allows the inspection of metal or ceramic knee implants. Even though we established an optimized scanning routine for tibial and femoral ceramic components, it is not possible to completely eliminate scanning artifacts of XCT.

CONCLUSION

Altogether, after visual inspection, none of the beam-hardening corrected XCT data sets for femoral and tibial implants showed any defects, that is, no inclusions, cracks, or pores were detected. XCT test is therefore an essential addition to the fatigue testing since it is the only non-destroying testing method.

Quantification of Tomographic Incompleteness in Cone-Beam Reconstruction

For situations of cone-beam scanning where the measurements are incomplete, we propose a method to quantify the severity of the missing information at each voxel. This incompleteness metric is geometric; it uses only the relative locations of all cone-beam vertices with respect to the voxel in question, and does not apply global information such as the object extent or the pattern of incompleteness of other voxels. The values are non-negative, with zero indicating "least incompleteness," i.e. minimal danger of incompleteness artifacts. The incompleteness value can be related to the severity of the potential reconstruction artifact at the voxel location, independent of reconstruction algorithm. We performed a computer simulation of x-ray sources along a circular trajectory, and used small multi-disk test-objects to examine the local effects of data incompleteness. The observed behavior of the reconstructed test-objects quantitatively matched the precalculated incompleteness values. A second simulation of a hypothetical SPECT breast imaging system used only 12 pinholes. Reconstructions were performed using analytic and iterative methods, and five reconstructed test-objects matched the behavior predicted by the incompleteness model. The model is based on known sufficiency conditions for data incompleteness, and provides strong predictive guidance for what can go wrong with incomplete cone-beam data.

A method to measure slice sensitivity profiles of CT images under low-contrast and high-noise conditions

Noise reduction features of iterative reconstruction (IR) methods in computed tomography might accompany the sacrifice of the longitudinal resolution, or slice sensitivity profile (SSP), at low contrast-to-noise ratio (CNR) conditions. To assess the benefit of IR methods correctly, the difference of SSP between IR methods and filtered-backprojection (FBP) must be taken into account. Therefore, SSP measurement under low-CNR conditions is necessary. Although edge methods are predominantly used, their performance under low-CNR conditions appears to be not fully established. We developed a method that is compatible with extremely low-CNR conditions. Thin plastic disk-shaped sheets embedded in acrylic resin were used as low-contrast test objects. The lowest peak contrast used was approximately 17 [HU]. We assessed the performance of our method by using FBP images. We identified a source of measurement instability aside from noise: the measured thin-slice SSP is dependent on the orbital phase of helical scan, presumably because of cone-beam artifacts. This impediment to high accuracy is manageable using phase-controlled scans. We confirmed that table position repeatability is much better than the value of the specifications, and therefore the ensemble-averaged images of multiple scans can be used for SSP measurement. Accurate measurement of SSP under extremely low-CNR conditions is possible, even when the test object is visually indiscernible from the noisy background. Low-contrast SSP behavior is elucidated for IR methods (AIDR-3D, FIRST, and AiSR-V) by using this measurement method.

A data-efficient method for local noise power spectrum (NPS) estimation in FDK-reconstructed 3D cone-beam CT

PURPOSE

For computed tomography (CT) systems in which noise is nonstationary, a local noise power spectrum (NPS) is often needed to characterize its noise property. We have previously developed a data-efficient radial NPS method to estimate the two-dimensional (2D) local NPS for filtered back projection (FBP)-reconstructed fan-beam CT utilizing the polar separability of CT NPS. In this work, we extend this method to estimate three-dimensional (3D) local NPS for feldkamp-davis-kress (FDK)-reconstructed cone-beam CT (CBCT) volumes.

METHODS

Starting from the 2D polar separability, we analyze the CBCT geometry and FDK image reconstruction process to derive the 3D expression of the polar separability for CBCT local NPS. With the polar separability, the 3D local NPS of CBCT can be decomposed into a 2D radial NPS shape function and a one-dimensional (1D) angular amplitude function with certain geometrical transforms. The 2D radial NPS shape function is a global function characterizing the noise correlation structure, while the 1D angular amplitude function is a local function reflecting the varying local noise amplitudes. The 3D radial local NPS method is constructed from the polar separability. We evaluate the accuracy of the 3D radial local NPS method using simulated and real CBCT data by comparing the radial local NPS estimates to a reference local NPS in terms of normalized mean squared error (NMSE) and a task-based performance metric (lesion detectability).

RESULTS

In both simulated and physical CBCT examples, a very small NMSE (<5%) was achieved by the radial local NPS method from as few as two scans, while for the traditional local NPS method, about 20 scans were needed to reach this accuracy. The results also showed that the detectability-based system performances computed using the local NPS estimated with the NPS method developed in this work from two scans closely reflected the actual system performance.

CONCLUSIONS

The polar separability greatly reduces the data dimensionality of the 3D CBCT local NPS. The radial local NPS method developed based on this property is shown to be capable of estimating the 3D local NPS from only two CBCT scans with acceptable accuracy. The minimum data requirement indicates the potential utility of local NPS in CBCT applications even for clinical situations.

Tomographic phase microscopy: principles and applications in bioimaging [Invited]

Tomographic phase microscopy (TPM) is an emerging optical microscopic technique for bioimaging. TPM uses digital holographic measurements of complex scattered fields to reconstruct three-dimensional refractive index (RI) maps of cells with diffraction-limited resolution by solving inverse scattering problems. In this paper, we review the developments of TPM from the fundamental physics to its applications in bioimaging. We first provide a comprehensive description of the tomographic reconstruction physical models used in TPM. The RI map reconstruction algorithms and various regularization methods are discussed. Selected TPM applications for cellular imaging, particularly in hematology, are reviewed. Finally, we examine the limitations of current TPM systems, propose future solutions, and envision promising directions in biomedical research.

Task-Driven Orbit Design and Implementation on a Robotic C-Arm System for Cone-Beam CT

PURPOSE

This work applies task-driven optimization to the design of non-circular orbits that maximize imaging performance for a particular imaging task. First implementation of task-driven imaging on a clinical robotic C-arm system is demonstrated, and a framework for orbit calculation is described and evaluated.

METHODS

We implemented a task-driven imaging framework to optimize orbit parameters that maximize detectability index d'. This framework utilizes a specified Fourier domain task function and an analytical model for system spatial resolution and noise. Two experiments were conducted to test the framework. First, a simple task was considered consisting of frequencies lying entirely on the f_z-axis (e.g., discrimination of structures oriented parallel to the central axial plane), and a "circle + arc" orbit was incorporated into the framework as a means to improve sampling of these frequencies, and thereby increase task-based detectability. The orbit was implemented on a robotic C-arm (Artis Zeego, Siemens Healthcare). A second task considered visualization of a cochlear implant simulated within a head phantom, with spatial frequency response emphasizing high-frequency content in the (f_y , f_z ) plane of the cochlea. An optimal orbit was computed using the task-driven framework, and the resulting image was compared to that for a circular orbit.

RESULTS

For the f_z -axis task, the circle + arc orbit was shown to increase d' by a factor of 1.20, with an improvement of 0.71 mm in a 3D edge-spread measurement for edges located far from the central plane and a decrease in streak artifacts compared to a circular orbit. For the cochlear implant task, the resulting orbit favored complementary views of high tilt angles in a 360° orbit, and d' was increased by a factor of 1.83.

CONCLUSIONS

This work shows that a prospective definition of imaging task can be used to optimize source-detector orbit and improve imaging performance. The method was implemented for execution of non-circular, task-driven orbits on a clinical robotic C-arm system. The framework is sufficiently general to include both acquisition parameters (e.g., orbit, kV, and mA selection) and reconstruction parameters (e.g., a spatially varying regularizer).

Fourier Analysis of Noise Characteristics in Cone-Beam Microtomography Laboratory Scanners

GOAL

We investigate the signal and noise performance of an x-ray microtomography system that incorporates a complementary metal-oxide-semiconductor flat-panel detector as a projection image receptor.

METHODS

Signal and noise performance is analyzed in the Fourier domain using modulation-transfer function (MTF), noise-power spectrum (NPS), and noise-equivalent number of quanta (NEQ) with respect to magnification and different convolution kernels for image reconstruction.

RESULTS

Higher magnification provides lower NPS, and thus, higher NEQ performance in the transaxial planes from microtomography. A window function capable of smoothing the ramp filter edge to below one-half of the Nyquist limit results in better performance in terms of NPS and NEQ. The characteristics of convolution kernels do not affect signal and noise performance in longitudinal planes; hence, MTF performance mainly dominates the NEQ performance. The signal and noise performances investigated in this study are demonstrated with images obtained from the contrast phantom and postmortem mouse.

CONCLUSION

The results of our study could be helpful in developing x-ray microtomography systems based on flat-panel detectors.

Accelerated statistical reconstruction for C-arm cone-beam CT using Nesterov's method

PURPOSE

To accelerate model-based iterative reconstruction (IR) methods for C-arm cone-beam CT (CBCT), thereby combining the benefits of improved image quality and/or reduced radiation dose with reconstruction times on the order of minutes rather than hours.

METHODS

The ordered-subsets, separable quadratic surrogates (OS-SQS) algorithm for solving the penalized-likelihood (PL) objective was modified to include Nesterov's method, which utilizes "momentum" from image updates of previous iterations to better inform the current iteration and provide significantly faster convergence. Reconstruction performance of an anthropomorphic head phantom was assessed on a benchtop CBCT system, followed by CBCT on a mobile C-arm, which provided typical levels of incomplete data, including lateral truncation. Additionally, a cadaveric torso that presented realistic soft-tissue and bony anatomy was imaged on the C-arm, and different projectors were assessed for reconstruction speed.

RESULTS

Nesterov's method provided equivalent image quality to OS-SQS while reducing the reconstruction time by an order of magnitude (10.0 ×) by reducing the number of iterations required for convergence. The faster projectors were shown to produce similar levels of convergence as more accurate projectors and reduced the reconstruction time by another 5.3 ×. Despite the slower convergence of IR with truncated C-arm CBCT, comparison of PL reconstruction methods implemented on graphics processing units showed that reconstruction time was reduced from 106 min for the conventional OS-SQS method to as little as 2.0 min with Nesterov's method for a volumetric reconstruction of the head. In body imaging, reconstruction of the larger cadaveric torso was reduced from 159 min down to 3.3 min with Nesterov's method.

CONCLUSIONS

The acceleration achieved through Nesterov's method combined with ordered subsets reduced IR times down to a few minutes. This improved compatibility with clinical workflow better enables broader adoption of IR in CBCT-guided procedures, with corresponding benefits in overcoming conventional limits of image quality at lower dose.

Extended ellipse-line-ellipse trajectory for long-object cone-beam imaging with a mounted C-arm system

Recent reports show that three-dimensional cone-beam (CB) imaging with a floor-mounted (or ceiling-mounted) C-arm system has become a valuable tool in interventional radiology. Currently, a circular short scan is used for data acquisition, which inevitably yields CB artifacts and a short coverage in the direction of the patient table. To overcome these two limitations, a more sophisticated data acquisition geometry is needed. This geometry should be complete in terms of Tuy's condition and should allow continuous scanning, while being compatible with the mechanical constraints of mounted C-arm systems. Additionally, the geometry should allow accurate image reconstruction from truncated data. One way to ensure such a feature is to adopt a trajectory that provides full R-line coverage within the field-of-view (FOV). An R-line is any segment of line that connects two points on a source trajectory, and the R-line coverage is the set of points that belong to an R-line. In this work, we propose a novel geometry called the extended ellipse-line-ellipse (ELE) for long-object imaging with a mounted C-arm system. This trajectory is built from modules consisting of two elliptical arcs connected by a line. We demonstrate that the extended ELE can be configured in many ways so that full R-line coverage is guaranteed. Both tight and relaxed parametric settings are presented. All results are supported by extensive mathematical proofs provided in appendices. Our findings make the extended ELE trajectory attractive for axially-extended FOV imaging in interventional radiology.

Comparative study of iterative reconstruction algorithms for missing cone problems in optical diffraction tomography

In optical tomography, there exist certain spatial frequency components that cannot be measured due to the limited projection angles imposed by the numerical aperture of objective lenses. This limitation, often called as the missing cone problem, causes the under-estimation of refractive index (RI) values in tomograms and results in severe elongations of RI distributions along the optical axis. To address this missing cone problem, several iterative reconstruction algorithms have been introduced exploiting prior knowledge such as positivity in RI differences or edges of samples. In this paper, various existing iterative reconstruction algorithms are systematically compared for mitigating the missing cone problem in optical diffraction tomography. In particular, three representative regularization schemes, edge preserving, total variation regularization, and the Gerchberg-Papoulis algorithm, were numerically and experimentally evaluated using spherical beads as well as real biological samples; human red blood cells and hepatocyte cells. Our work will provide important guidelines for choosing the appropriate regularization in ODT.

Task-based detectability in CT image reconstruction by filtered backprojection and penalized likelihood estimation

PURPOSE

Nonstationarity is an important aspect of imaging performance in CT and cone-beam CT (CBCT), especially for systems employing iterative reconstruction. This work presents a theoretical framework for both filtered-backprojection (FBP) and penalized-likelihood (PL) reconstruction that includes explicit descriptions of nonstationary noise, spatial resolution, and task-based detectability index. Potential utility of the model was demonstrated in the optimal selection of regularization parameters in PL reconstruction.

METHODS

Analytical models for local modulation transfer function (MTF) and noise-power spectrum (NPS) were investigated for both FBP and PL reconstruction, including explicit dependence on the object and spatial location. For FBP, a cascaded systems analysis framework was adapted to account for nonstationarity by separately calculating fluence and system gains for each ray passing through any given voxel. For PL, the point-spread function and covariance were derived using the implicit function theorem and first-order Taylor expansion according to Fessler ["Mean and variance of implicitly defined biased estimators (such as penalized maximum likelihood): Applications to tomography," IEEE Trans. Image Process. 5(3), 493-506 (1996)]. Detectability index was calculated for a variety of simple tasks. The model for PL was used in selecting the regularization strength parameter to optimize task-based performance, with both a constant and a spatially varying regularization map.

RESULTS

Theoretical models of FBP and PL were validated in 2D simulated fan-beam data and found to yield accurate predictions of local MTF and NPS as a function of the object and the spatial location. The NPS for both FBP and PL exhibit similar anisotropic nature depending on the pathlength (and therefore, the object and spatial location within the object) traversed by each ray, with the PL NPS experiencing greater smoothing along directions with higher noise. The MTF of FBP is isotropic and independent of location to a first order approximation, whereas the MTF of PL is anisotropic in a manner complementary to the NPS. Task-based detectability demonstrates dependence on the task, object, spatial location, and smoothing parameters. A spatially varying regularization "map" designed from locally optimal regularization can improve overall detectability beyond that achievable with the commonly used constant regularization parameter.

CONCLUSIONS

Analytical models for task-based FBP and PL reconstruction are predictive of nonstationary noise and resolution characteristics, providing a valuable framework for understanding and optimizing system performance in CT and CBCT.

Soft-tissue imaging with C-arm cone-beam CT using statistical reconstruction

The potential for statistical image reconstruction methods such as penalized-likelihood (PL) to improve C-arm cone-beam CT (CBCT) soft-tissue visualization for intraoperative imaging over conventional filtered backprojection (FBP) is assessed in this work by making a fair comparison in relation to soft-tissue performance. A prototype mobile C-arm was used to scan anthropomorphic head and abdomen phantoms as well as a cadaveric torso at doses substantially lower than typical values in diagnostic CT, and the effects of dose reduction via tube current reduction and sparse sampling were also compared. Matched spatial resolution between PL and FBP was determined by the edge spread function of low-contrast (∼ 40-80 HU) spheres in the phantoms, which were representative of soft-tissue imaging tasks. PL using the non-quadratic Huber penalty was found to substantially reduce noise relative to FBP, especially at lower spatial resolution where PL provides a contrast-to-noise ratio increase up to 1.4-2.2 × over FBP at 50% dose reduction across all objects. Comparison of sampling strategies indicates that soft-tissue imaging benefits from fully sampled acquisitions at dose above ∼ 1.7 mGy and benefits from 50% sparsity at dose below ∼ 1.0 mGy. Therefore, an appropriate sampling strategy along with the improved low-contrast visualization offered by statistical reconstruction demonstrates the potential for extending intraoperative C-arm CBCT to applications in soft-tissue interventions in neurosurgery as well as thoracic and abdominal surgeries by overcoming conventional tradeoffs in noise, spatial resolution, and dose.

Cascaded systems analysis of noise and detectability in dual-energy cone-beam CT

PURPOSE

Dual-energy computed tomography and dual-energy cone-beam computed tomography (DE-CBCT) are promising modalities for applications ranging from vascular to breast, renal, hepatic, and musculoskeletal imaging. Accordingly, the optimization of imaging techniques for such applications would benefit significantly from a general theoretical description of image quality that properly incorporates factors of acquisition, reconstruction, and tissue decomposition in DE tomography. This work reports a cascaded systems analysis model that includes the Poisson statistics of x rays (quantum noise), detector model (flat-panel detectors), anatomical background, image reconstruction (filtered backprojection), DE decomposition (weighted subtraction), and simple observer models to yield a task-based framework for DE technique optimization.

METHODS

The theoretical framework extends previous modeling of DE projection radiography and CBCT. Signal and noise transfer characteristics are propagated through physical and mathematical stages of image formation and reconstruction. Dual-energy decomposition was modeled according to weighted subtraction of low- and high-energy images to yield the 3D DE noise-power spectrum (NPS) and noise-equivalent quanta (NEQ), which, in combination with observer models and the imaging task, yields the dual-energy detectability index (d(')). Model calculations were validated with NPS and NEQ measurements from an experimental imaging bench simulating the geometry of a dedicated musculoskeletal extremities scanner. Imaging techniques, including kVp pair and dose allocation, were optimized using d(') as an objective function for three example imaging tasks: (1) kidney stone discrimination; (2) iodine vs bone in a uniform, soft-tissue background; and (3) soft tissue tumor detection on power-law anatomical background.

RESULTS

Theoretical calculations of DE NPS and NEQ demonstrated good agreement with experimental measurements over a broad range of imaging conditions. Optimization results suggest a lower fraction of total dose imparted by the low-energy acquisition, a finding consistent with previous literature. The selection of optimal kVp pair reveals the combined effect of both quantum noise and contrast in the kidney stone discrimination and soft-tissue tumor detection tasks, whereas the K-edge effect of iodine was the dominant factor in determining kVp pairs in the iodine vs bone task. The soft-tissue tumor task illustrated the benefit of dual-energy imaging in eliminating anatomical background noise and improving detectability beyond that achievable by single-energy scans.

CONCLUSIONS

This work established a task-based theoretical framework that is predictive of DE image quality. The model can be utilized in optimizing a broad range of parameters in image acquisition, reconstruction, and decomposition, providing a useful tool for maximizing DE-CBCT image quality and reducing dose.

Beyond noise power in 3D computed tomography: the local NPS and off-diagonal elements of the Fourier domain covariance matrix

PURPOSE

To investigate the correlation and stationarity of noise in volumetric computed tomography (CT) using the local discrete noise-power spectrum (NPS) and off-diagonal elements of the covariance matrix of the discrete Fourier transform of noise-only images (denoted Σ(DFT)). Experimental conditions were varied to affect noise correlation and stationarity, the effects were quantified in terms of the NPS and Σ(DFT), and practical considerations in CT performance characterization were identified.

METHODS

Cone-beam CT (CBCT) images were acquired using a benchtop system comprising an x-ray tube and flat-panel detector for a range of acquisition techniques (e.g., dose and x-ray scatter) and three phantom configurations hypothesized to impart distinct effects on the NPS and Σ(DFT): (A) air, (B) a 20-cm-diameter water cylinder with a bowtie filter, and (C) the cylinder without a bowtie filter. The NPS and off-diagonal elements of the Σ(DFT) were analyzed as a function of position within the reconstructions.

RESULTS

The local NPS varied systematically throughout the axial plane in a manner consistent with changes in fluence transmitted to the detector and view sampling effects. Variability in fluence was manifest in the NPS magnitude-e.g., a factor of ~2 variation in NPS magnitude within the axial plane for case C (cylinder without bowtie), compared to nearly constant NPS magnitude for case B (bowtie filter matched to the cylinder). View sampling effects were most prominent in case A (air) where the variance increased at greater distance from the center of reconstruction and in case C (cylinder) where the NPS exhibited correlations in the radial direction. The effects of detector lag were observed as azimuthal correlation. The cylinder (without bowtie) had the strongest nonstationarity because of the larger variability in fluence transmitted to the detector. The diagonal elements of the Σ(DFT) were equivalent to the NPS estimated from the periodogram, and the average off-diagonal elements of the Σ(DFT) exhibited amplitude of ~1% of the NPS for the experimental conditions investigated. Furthermore, the off-diagonal elements demonstrated fairly long tails of nearly constant amplitude, with magnitude somewhat reduced for experimental conditions associated with greater stationarity (viz., lower Σ(DFT) tails for cases A and B in comparison to case C).

CONCLUSIONS

Volumetric CT exhibits nonstationarity in the NPS as hypothesized in relation to fluence uniformity and view sampling. Measurement of the NPS should seek to minimize such changes in noise correlations and include careful reporting of experimental conditions (e.g., phantom design and use of a bowtie filter) and spatial dependence (e.g., analysis at fixed radius within a phantom). Off-diagonal elements of the Σ(DFT) similarly depend on experimental conditions and can be readily computed from the same data as the NPS. This work begins to check assumptions in NPS analysis examine the extent to which NPS is an appropriate descriptor of noise correlations, and investigate the magnitude of off-diagonal elements of the Σ(DFT). While the magnitude of such off-diagonal elements appears to be low, their cumulative effect on space-variant detectability remains to be investigated-e.g., using task-specific figures of merit.

Frameless angiogram-based stereotactic radiosurgery for treatment of arteriovenous malformations

PURPOSE

Stereotactic radiosurgery (SRS) is an effective alternative to microsurgical resection or embolization for definitive treatment of arteriovenous malformations (AVMs). Digital subtraction angiography (DSA) is the gold standard for pretreatment diagnosis and characterization of vascular anatomy, but requires rigid frame (skull) immobilization when used in combination with SRS. With the advent of advanced proton and image-guided photon delivery systems, SRS treatment is increasingly migrating to frameless platforms, which are incompatible with frame-based DSA. Without DSA as the primary image, target definition may be less than optimal, in some cases precluding the ability to treat with a frameless system. This article reports a novel solution.

METHODS AND MATERIALS

Fiducial markers are implanted into the patient's skull before angiography. Angiography is performed according to the standard clinical protocol, but, in contrast to the previous practice, without the rigid frame. Separate images of a specially designed localizer box are subsequently obtained. A target volume projected on DSA can be transferred to the localizer system in three dimensions, and in turn be transferred to multiple CT slices using the implanted fiducials. Combined with other imaging modalities, this "virtual frame" approach yields a highly precise treatment plan that can be delivered by frameless SRS technologies.

RESULTS

Phantom measurements for point and volume targets have been performed. The overall uncertainty of placing a point target to CT is 0.4 mm. For volume targets, deviation of the transformed contour from the target CT image is within 0.6 mm. The algorithm and software are robust. The method has been applied clinically, with reliable results.

CONCLUSIONS

A novel and reproducible method for frameless SRS of AVMs has been developed that enables the use of DSA without the requirement for rigid immobilization. Multiple pairs of DSA can be used for better conformality. Further improvement, including using nonimplanted fiducials, is potentially feasible.

A new method to combine 3D reconstruction volumes for multiple parallel circular cone beam orbits

PURPOSE

This article presents a new reconstruction method for 3D imaging using a multiple 360 degrees circular orbit cone beam CT system, specifically a way to combine 3D volumes reconstructed with each orbit. The main goal is to improve the noise performance in the combined image while avoiding cone beam artifacts.

METHODS

The cone beam projection data of each orbit are reconstructed using the FDK algorithm. When at least a portion of the total volume can be reconstructed by more than one source, the proposed combination method combines these overlap regions using weighted averaging in frequency space. The local exactness and the noise performance of the combination method were tested with computer simulations of a Defrise phantom, a FORBILD head phantom, and uniform noise in the raw data.

RESULTS

A noiseless simulation showed that the local exactness of the reconstructed volume from the source with the smallest tilt angle was preserved in the combined image. A noise simulation demonstrated that the combination method improved the noise performance compared to a single orbit reconstruction.

CONCLUSIONS

In CT systems which have overlap volumes that can be reconstructed with data from more than one orbit and in which the spatial frequency content of each reconstruction can be calculated, the proposed method offers improved noise performance while keeping the local exactness of data from the source with the smallest tilt angle.

Quantitative imaging in breast tomosynthesis and CT: comparison of detection and estimation task performance

PURPOSE

This work investigates a framework for modeling volumetric breast imaging to compare detection and estimation task performance and optimize quantitative breast imaging.

METHODS

Volumetric reconstructions of a breast phantom, which incorporated electronic, quantum, and anatomical noise with embedded spherical lesions, were simulated over a range of acquisition angles varying from 4 degrees to 204 degrees with a constant total acquisition dose of 1.5 mGy. A maximum likelihood estimator was derived in terms of the noise power spectrum, which yielded figures of merit for quantitative imaging performance in terms of accuracy and precision. These metrics were computed for estimation of lesion area, volume, and location. Estimation task performance was optimized as a function of acquisition angle and compared to the performance of a more conventional lesion detection task.

RESULTS

Results revealed tradeoffs between electronic, quantum, and anatomical noise. The detection of a 4 mm sphere was optimal at an acquisition angle of 84 degrees, where reconstructed images using a smaller acquisition angle exhibited increased anatomical noise and reconstructed images using a larger acquisition angle exhibited increased quantum and electronic noise. For all estimation tasks, accuracy was found to be fairly constant as a function acquisition angle indicating adequate system calibration, whereas a more significant dependence on acquisition angle was observed for precision performance. Precision for the 2D area estimation task was optimal at approximately 104 degrees, while precision of the 3D volume estimation task was optimal at larger angles (approximately 124 degrees). Precision for the localization task showed orientation dependence where localization was significantly inferior in the depth direction. Overall, precision for localization was optimal at larger angles (i.e., > 125 degrees) compared to the size estimation tasks. Results suggested that for quantitative imaging tasks, the acquisition angle should be larger than currently used in conventional breast tomosynthesis for lesion detection.

CONCLUSIONS

Analysis of quantitative imaging performance using Fourier-based metrics highlights the difference between estimation and detection task in volumetric breast imaging and provides a meaningful framework for optimizing the performance of breast imaging systems for quantitative imaging applications.

Anatomical background and generalized detectability in tomosynthesis and cone-beam CT

PURPOSE

Anatomical background presents a major impediment to detectability in 2D radiography as well as 3D tomosynthesis and cone-beam CT (CBCT). This article incorporates theoretical and experimental analysis of anatomical background "noise" in cascaded systems analysis of 2D and 3D imaging performance to yield "generalized" metrics of noise-equivalent quanta (NEQ) and detectability index as a function of the orbital extent of the (circular arc) source-detector orbit.

METHODS

A physical phantom was designed based on principles of fractal self-similarity to exhibit power-law spectral density (kappa/Fbeta) comparable to various anatomical sites (e.g., breast and lung). Background power spectra [S(B)(F)] were computed as a function of source-detector orbital extent, including tomosynthesis (approximately 10 degrees -180 degrees) and CBCT (180 degrees + fan to 360 degrees) under two acquisition schemes: (1) Constant angular separation between projections (variable dose) and (2) constant total number of projections (constant dose). The resulting S(B) was incorporated in the generalized NEQ, and detectability index was computed from 3D cascaded systems analysis for a variety of imaging tasks.

RESULTS

The phantom yielded power-law spectra within the expected spatial frequency range, quantifying the dependence of clutter magnitude (kappa) and correlation (beta) with increasing tomosynthesis angle. Incorporation of S(B) in the 3D NEQ provided a useful framework for analyzing the tradeoffs among anatomical, quantum, and electronic noise with dose and orbital extent. Distinct implications are posed for breast and chest tomosynthesis imaging system design-applications varying significantly in kappa and beta, and imaging task and, therefore, in optimal selection of orbital extent, number of projections, and dose. For example, low-frequency tasks (e.g., soft-tissue masses or nodules) tend to benefit from larger orbital extent and more fully 3D tomographic imaging, whereas high-frequency tasks (e.g., microcalcifications) require careful, application-specific selection of orbital extent and number of projections to minimize negative effects of quantum and electronic noise.

CONCLUSIONS

The complex tradeoffs among anatomical background, quantum noise, and electronic noise in projection imaging, tomosynthesis, and CBCT can be described by generalized cascaded systems analysis, providing a useful framework for system design and optimization.

Implementation and characterization of a 320-slice volumetric CT scanner for simulation in radiation oncology

PURPOSE

Effective target definition and broad employment of treatment response assessment with dynamic contrast-enhanced CT in radiation oncology requires increased speed and coverage for use within a single bolus injection. To this end, a novel volumetric CT scanner (Aquilion One, Toshiba, Tochigi Pref., Japan) has been installed at the Princess Margaret Hospital for implementation into routine CT simulation. This technology offers great advantages for anatomical and functional imaging in both scan speed and coverage. The aim of this work is to investigate the system's imaging performance and quality as well as CT quantification accuracy which is important for radiotherapy dose calculations.

METHODS

The 320-slice CT scanner uses a 160 mm wide-area (2D) solid-state detector design which provides the possibility to acquire a volumetric axial length of 160 mm without moving the CT couch. This is referred to as "volume" and can be scanned with a rotation speed of 0.35-3 s. The scanner can also be used as a 64-slice CT scanner and perform conventional (axial) and helical acquisitions with collimation ranges of 1-32 and 16-32 mm, respectively. Commissioning was performed according to AAPM Reports TG 66 and 39 for both helical and volumetric imaging. Defrise and other cone-beam image analysis tests were performed.

RESULTS

Overall, the imaging spatial resolution and geometric efficiency (GE) were found to be very good (>10 lp/mm, <1 mm spatial integrity and GE160 mm=85%) and within the AAPM guidelines as well as IEC recommendations. Although there is evidence of some cone-beam artifacts when scanning the Defrise phantom, image quality was found to be good and sufficient for treatment planning (soft tissue noise <10 HU). Measurements of CT number stability and contrast-to-noise values across the volume indicate clinically acceptable scan accuracy even at the field edge.

CONCLUSIONS

Initial experience with this exciting new technology confirms its accuracy for routine CT simulation within radiation oncology and allows for future investigations into specialized dynamic volumetric imaging applications.