Noise aliasing and the 3D NEQ of flat-panel cone-beam CT: effect of 2D/3D apertures and sampling

Technical note: Does scan resolution or downsampling impact the analysis of trabecular bone architecture?

The "gold standard" for the assessment of trabecular bone structure is high-resolution micro-CT. In this technical note, we test the influence of initial scan resolution and post hoc downsampling on the quantitative and qualitative analysis of trabecular bone in a Gorilla tibia. We analyzed trabecular morphology in the right distal tibia of one Gorilla gorilla individual to investigate the impact of variation in voxel size on measured trabecular variables. For each version of the micro-CT volume, trabecular bone was segmented using the medical image analysis method. Holistic morphometric analysis was then used to analyze bone volume (BV/TV), anisotropy (DA), trabecular thickness (Tb.Th), spacing (Tb.Sp), and number (Tb.N). Increasing voxel size during initial scanning was found to have a strong impact on DA and Tb.Th measures, while BV/TV, Tb.Sp, and Tb.N were found to be less sensitive to variations in initial scan resolution. All tested parameters were not substantially influenced by downsampling up to 90 μm resolution. Color maps of BV/TV and DA also retained their distribution up to 90 μm. This study is the first to examine the effect of variation in micro-CT voxel size on the analysis of trabecular bone structure using whole epiphysis approaches. Our results indicate that microstructural variables may be measured for most trabecular parameters up to a voxel size of 90 μm for both scan and downsampled resolutions. Moreover, if only BV/TV, Tb.Sp or Tb.N is measured, even larger voxel sizes might be used without substantially affecting the results.

Task-based image quality assessment of an intraoperative CBCT for spine surgery compared with conventional CT

PURPOSE

To analyze the image quality of a novel, state-of-the art platform for CBCT image-guided spine surgery, focusing particularly on the dose-effectiveness compared with conventional CT (the gold standard for postoperative assessment).

METHODS

The ClarifEye platform (Philips Healthcare) with integrated augmented-reality surgical navigation, has been compared with a GE Revolution CT (GE Healthcare). The 3D spatial resolution (TTF) and noise (NPS) were evaluated considering relevant feature contrasts (200-900 HU) and background noise for differently sized patients (200-300 mm water-equivalent diameter). These measures were used to determine the noise equivalent quanta (NEQ) and observer model detectability.

RESULTS

The CBCT system exhibited a linear response with 50% TTF at 5.7 cycles/cm (10% TTF at 9.2 cycles/cm), and the axial noise power peaking at about 3.6 cycles/cm (average frequency of 4.1 cycles/cm). The noise magnitude and texture differed markedly compared to iteratively reconstructed CT images (GE ASiR-V). The CBCT system had 26% lower detectability for a high-frequency task (related to edge detection) compared with CT images reconstructed using the Bone kernel combined with ASiR-V 50%. Likewise, it had 18% lower detectability for low- and mid-frequency tasks compared with CT images reconstructed using the Standard kernel. This difference translates to 50%-80% higher CBCT imaging doses required to match the CT image quality.

CONCLUSIONS

The ClarifEye platform demonstrates intraoperative CBCT-imaging capabilities that under certain circumstances are comparable with conventional CT. However, due to limited dose-effectiveness, a trade-off between timeliness and radiation exposure must be considered if end-of-procedure CBCT is to replace postoperative CT.

Noise power spectrum (NPS) in computed tomography: Enabling local NPS measurement without stationarity and ergodicity assumptions

BACKGROUND

Conventional methods for estimating the noise power spectrum (NPS) often necessitate multiple computed tomography (CT) data acquisitions and are required to satisfy stringent stationarity and ergodicity conditions, which prove challenging in CT imaging systems.

PURPOSE

The aim was to revisit the conventional NPS estimation method, leading to a new framework that estimates local NPS without relying on stationarity or ergodicity, thus facilitating experimental NPS estimations.

METHODS

The scientific foundation of the conventional CT NPS measurement method, based on the Wiener-Khintchine theorem, was reexamined, emphasizing the critical conditions of stationarity and ergodicity. This work proposes an alternative framework, characterized by its independence from stationarity and ergodicity, and its ability to facilitate local NPS estimations. A spatial average of local NPS over a Region of Interest (ROI) yields the conventional NPS for that ROI. The connections and differences between the proposed alternative method and the conventional method are discussed. Experimental studies were conducted to validate the new method.

RESULTS

(1) The NPS estimated using the conventional method was demonstrated to correspond to the spatial average of pointwise NPS from the proposed NPS estimation framework. (2) The NPS estimated over an ROI with the conventional method was shown to be the sum of the NPS estimated from the proposed method and a contribution from measurement uncertainty. (3) Local NPS estimations from the proposed method in this work elucidate the impact of surrounding image content on local NPS variations.

CONCLUSION

The NPS estimation method proposed in this work allows for the estimation of local NPS without relying on stationarity and ergodicity conditions, offering local NPS estimations with significantly improved precision.

Experimental measurement of local noise power spectrum (NPS) in photon counting detector-CT (PCD-CT) using a single data acquisition

BACKGROUND

Accurate noise power spectra (NPS) measurement in clinical X-ray CT exams is challenging due to the need for repeated scans, which expose patients to high radiation risks. A reliable method for single CT acquisition NPS estimation is thus highly desirable.

PURPOSE

To develop a method for estimating local NPS from a single photon counting detector-CT (PCD-CT) acquisition.

METHODS

A novel nearly statistical bias-free estimator was constructed from the raw counts data of PCD-CT scan to estimate the variance of sinogram projection data. An analytical algorithm is employed to reconstruct point-wise covariancecov ( x i , x j ) $\text{cov}({\bf x}_i,{\bf x}_j)$ between any two image pixel/voxel locationsx i ${\bf x}_i$ andx j ${\bf x_j}$ . A Fourier transform is applied to obtain the desired point-wise NPS for any chosen locationx i ${\bf x}_i$ . The method was validated using experimental data acquired from a benchtop PCD-CT system with various physical phantoms, and the results were compared with the conventional local NPS measurement method using repeated scans and statistical ensemble averaging.

RESULTS

The experimental results demonstrate that (1) the proposed method can achieve pointwise/local NPS measurement for a region of interest (ROI) located at any chosen position, accurately characterizing the NPS with spatial structures resulting from image content heterogeneity; (2) the local NPS measured using the proposed method show a higher precision in the measured NPS compared to the conventional measurement method; (3) spatial averaging of the local NPS yields the conventional NPS for a given local ROI.

CONCLUSION

A new method was developed to enable local NPS from a single PCD-CT acquisition.

Technical note: Characterization, validation, and spectral optimization of a dedicated breast CT system for contrast-enhanced imaging

BACKGROUND

The development of a new imaging modality, such as 4D dynamic contrast-enhanced dedicated breast CT (4D DCE-bCT), requires optimization of the acquisition technique, particularly within the 2D contrast-enhanced imaging modality. Given the extensive parameter space, cascade-systems analysis is commonly used for such optimization.

PURPOSE

To implement and validate a parallel-cascaded model for bCT, focusing on optimizing and characterizing system performance in the projection domain to enhance the quality of input data for image reconstruction.

METHODS

A parallel-cascaded system model of a state-of-the-art bCT system was developed and model predictions of the presampled modulation transfer function (MTF) and the normalized noise power spectrum (NNPS) were compared with empirical data collected in the projection domain. Validation was performed using the default settings of 49 kV with 1.5 mm aluminum filter and at 65 kV and 0.257 mm copper filter. A 10 mm aluminum plate was added to replicate the breast attenuation. Air kerma at the isocenter was measured at different tube current levels. Discrepancies between the measured projection domain metrics and model-predicted values were quantified using percentage error and coefficient of variation (CoV) for MTF and NNPS, respectively. The optimal filtration was for a 5 mm iodine disk detection task at 49, 55, 60, and 65 kV. The detectability index was calculated for the default aluminum filtration and for copper thicknesses ranging from 0.05 to 0.4 mm.

RESULTS

At 49 kV, MTF errors were +5.1% and -5.1% at 1 and 2 cycles/mm, respectively; NNPS CoV was 5.3% (min = 3.7%; max = 8.5%). At 65 kV, MTF errors were -0.8% and -3.2%; NNPS CoV was 13.1% (min = 11.4%; max = 16.9%). Air kerma output was linear, with 11.67 µGy/mA (R2 = 0.993) and 19.14 µGy/mA (R2 = 0.996) at 49 and 65 kV, respectively. For iodine detection, a 0.25 mm-thick copper filter at 65 kV was found optimal, outperforming the default technique by 90%.

CONCLUSION

The model accurately predicts bCT system performance, specifically in the projection domain, under varied imaging conditions, potentially contributing to the enhancement of 2D contrast-enhanced imaging in 4D DCE-bCT.

Image quality models for 2D and 3D x-ray imaging systems: A perspective vignette

Image quality models based on cascaded systems analysis and task-based imaging performance were an important aspect of the emergence of 2D and 3D digital x-ray systems over the last 25 years. This perspective vignette offers cursory review of such developments and personal insights that may not be obvious within previously published scientific literature. The vignette traces such models to the mid-1990s, when flat-panel x-ray detectors were emerging as a new base technology for digital radiography and benefited from the rigorous, objective characterization of imaging performance gained from such models. The connection of models for spatial resolution and noise to spatial-frequency-dependent descriptors of imaging task provided a useful framework for system optimization that helped to accelerate the development of new technologies to first clinical use. Extension of the models to new technologies and applications is also described, including dual-energy imaging, photon-counting detectors, phase contrast imaging, tomosynthesis, cone-beam CT, 3D image reconstruction, and image registration.

3D vertebrae labeling in spine CT: an accurate, memory-efficient (Ortho2D) framework

Purpose.Accurate localization and labeling of vertebrae in computed tomography (CT) is an important step toward more quantitative, automated diagnostic analysis and surgical planning. In this paper, we present a framework (called Ortho2D) for vertebral labeling in CT in a manner that is accurate and memory-efficient.Methods. Ortho2D uses two independent faster R-convolutional neural network networks to detect and classify vertebrae in orthogonal (sagittal and coronal) CT slices. The 2D detections are clustered in 3D to localize vertebrae centroids in the volumetric CT and classify the region (cervical, thoracic, lumbar, or sacral) and vertebral level. A post-process sorting method incorporates the confidence in network output to refine classifications and reduce outliers. Ortho2D was evaluated on a publicly available dataset containing 302 normal and pathological spine CT images with and without surgical instrumentation. Labeling accuracy and memory requirements were assessed in comparison to other recently reported methods. The memory efficiency of Ortho2D permitted extension to high-resolution CT to investigate the potential for further boosts to labeling performance.Results. Ortho2D achieved overall vertebrae detection accuracy of 97.1%, region identification accuracy of 94.3%, and individual vertebral level identification accuracy of 91.0%. The framework achieved 95.8% and 83.6% level identification accuracy in images without and with surgical instrumentation, respectively. Ortho2D met or exceeded the performance of previously reported 2D and 3D labeling methods and reduced memory consumption by a factor of ∼50 (at 1 mm voxel size) compared to a 3D U-Net, allowing extension to higher resolution datasets than normally afforded. The accuracy of level identification increased from 80.1% (for standard/low resolution CT) to 95.1% (for high-resolution CT).Conclusions. The Ortho2D method achieved vertebrae labeling performance that is comparable to other recently reported methods with significant reduction in memory consumption, permitting further performance boosts via application to high-resolution CT.

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.

Effects of kV, filtration, dose, and object size on soft tissue and iodine contrast in dedicated breast CT

PURPOSE

Clinical use of dedicated breast computed tomography (bCT) requires relatively short scan times necessitating systems with high frame rates. This in turn impacts the x-ray tube operating range. We characterize the effects of tube voltage, beam filtration, dose, and object size on contrast and noise properties related to soft tissue and iodine contrast agents as a way to optimize imaging protocols for soft tissue and iodine contrast at high frame rates.

METHODS

This study design uses the signal-difference-to-noise ratio (SDNR), noise-equivalent quanta (NEQ), and detectability (d´) as measures of imaging performance for a prototype breast CT scanner that utilizes a pulsed x-ray tube (with a 4 ms pulse width) at 43.5 fps acquisition rate. We assess a range of kV, filtration, breast phantom size, and mean glandular dose (MGD). Performance measures are estimated from images of adipose-equivalent breast phantoms machined to have a representative size and shape of small, medium, and large breasts. Water (glandular tissue equivalent) and iodine contrast (5 mg/ml) were used to fill two cylindrical wells in the phantoms.

RESULTS

Air kerma levels required for obtaining an MGD of 6 mGy ranged from 7.1 to 9.1 mGy and are reported across all kV, filtration, and breast phantom sizes. However, at 50 kV, the thick filters (0.3 mm of Cu or Gd) exceeded the maximum available mA of the x-ray generator, and hence, these conditions were excluded from subsequent analysis. There was a strong positive association between measurements of SDNR and d' (R²  > 0.97) within the range of parameters investigated in this work. A significant decrease in soft tissue SDNR was observed for increasing phantom size and increasing kV with a maximum SDNR at 50 kV with 0.2 mm Cu or 0.2 mm Gd filtration. For iodine contrast SDNR, a significant decrease was observed with increasing phantom size, but a decrease in SDNR for increasing kV was only observed for 70 kV (50 and 60 kV were not significantly different). Thicker Gd filtration (0.3 mm Gd) resulted in a significant increase in iodine SDNR and decrease in soft tissue SDNR but requires significantly more tube current to deliver the same MGD.

CONCLUSIONS

The choice of 60 kV with 0.2 mm Gd filtration provides a good trade-off for maximizing both soft tissue and iodine contrast. This scanning technique takes advantage of the ~50 keV Gd k-edge to produce contrast and can be achieved within operating range of the x-ray generator used in this work. Imaging at 60 kV allows for a greater range in dose delivered to the large breast sizes when uniform image quality is desired across all breast sizes. While imaging performance metrics (i.e., detectability index and SDNR) were shown to be strongly correlated, the methodologies presented in this work for the estimation of NEQ (and subsequently d') provides a meaningful description of the spatial resolution and noise characteristics of this prototype bCT system across a range of beam quality, dose, and object sizes.

Task-driven optimization of the non-spectral mode of photon counting CT for intracranial hemorrhage assessment

Non-contrast CT (NCCT) is widely employed as the first-line imaging test to evaluate intracranial hemorrhage (ICH). Advances in mutidetector CT (MDCT) technology have greatly improved the image quality of NCCT for the detection of established, relatively large, and acute ICHs. Meanwhile, the reliability of MDCT in detecting microbleeds and chronic hemorrhage, and in predicting hemorrhagic transformation needs to be further improved. The purpose of this work was to investigate the potential use of non-spectral photon counting CT (PCCT) to address these challenges in ICH imaging. Towards this goal, the NCCT protocol of an experimental PCCT system that simulates the geometry of a general-purpose MDCT was optimized. The optimization was driven by three imaging tasks: detection of a 4.0 mm intraparenchymal hemorrhage, detection of a 1.5 mm subarachnoid hemorrhage, and discrimination of a sulcus in the insular cortex from the parenchymal background. These imaging tasks were custom-built into an anthropomorphic head phantom. Under the guidance of the frequency-dependent noise equivalent quanta and the ideal observer model detectability index [Formula: see text], the optimal PCD detection mode, energy threshold, and reconstruction kernel were found to be the anti-charge sharing mode, 15 keV, and an apodized ramp kernel, respectively. Compared with a clinical MDCT operated with an ICH protocol and at a matched dose level, the PCCT system provided at least 20% improvements in [Formula: see text] for all three ICH imaging tasks. These results demonstrated the potential benefits of non-spectral PCCT in ICH assessment.

Dose Efficiency of Quarter-Millimeter Photon-Counting Computed Tomography: First-in-Human Results

PURPOSE

The aim of this study was to assess the clinical feasibility, image quality, and radiation dose implications of 0.25-mm imaging mode in a cohort of humans, achieved by dividing the photon-counting detector (PCD) size in half compared with standard-resolution photon-counting computed tomography (CT) (0.5 mm).

METHODS

In this technical feasibility study, a whole-body prototype PCD-CT scanner was studied in the 0.25 mm detector mode (measured at isocenter). A high-resolution PCD-CT protocol was first tested in phantom and canine studies in terms of image noise and spatial resolution. Then, 8 human subjects (mean age, 58 ± 8 years; 2 men) underwent axial PCD 0.25-mm scans of the brain, the thorax, and at the level of the upper left kidney. Filtered backprojection reconstruction was performed with a sharp kernel (B70) for standard-resolution and high-resolution data at 0.5-mm isotropic image voxel. High-resolution data, in addition, were reconstructed with an ultrasharp kernel (U70) at 0.25-mm isotropic voxels.

RESULTS

Image reconstructions from the PCD 0.25-mm detector system led to an improvement in resolution from 9 to 18 line pairs/cm in a line pair phantom. Modulation transfer function improved from 9.5 to 15.8 line pairs/cm at 10% modulation transfer function. When fully exploiting this improvement, image noise increased by 75% compared with dose-matched 0.5-mm slice PCD standard-resolution acquisition. However, when comparing with standard-resolution data at same in-plane resolution and slice thickness, the PCD 0.25-mm detector mode showed 19% less image noise in phantom, animal, and human scans.

CONCLUSION

High-resolution photon-counting CT in humans showed improved image quality in terms of spatial resolution and image noise compared with standard-resolution photon-counting.

Quarter-millimeter spectral coronary stent imaging with photon-counting CT: Initial experience

PURPOSE

To evaluate the performance and clinical feasibility of 0.25 mm resolution mode of a dual-energy photon-counting detector (PCD) computed tomography (CT) system for coronary stent imaging and to compare the results to state-of-the-art dual-energy energy-integrating detector (EID) CT.

MATERIALS AND METHODS

Coronary stents with different diameters (2.0-4.0 mm) were examined inside a coronary artery phantom consisting of plastic tubes filled with iodine-based and gadolinium-based contrast material diluted to approximate clinical concentrations (n = 18). EID images were acquired using 2nd and 3rd generation dual-source CT systems (SOMATOM Flash and SOMATOM Force, Siemens Healthcare) at 0.60 mm (defined as standard-resolution (SR)) isotropic voxel size. Radiation-dose matched PCD images were acquired using a human prototype PCD system (Siemens Healthcare) at 0.50 mm (SR) and 0.25 mm (HR) imaging modes. Images were reconstructed using optimized convolution kernels.

RESULTS

Dual-energy HR PCD images significantly better stent lumen visualization (median: 69.5%, IQR: 61.2-78.9%) over dual-energy EID, and standard-resolution PCD images (median: 53.2-57.4%, all P < 0.01). HR PCD acquisitions reconstructed at SR image voxel size showed 25.3% lower image noise compared to SR PCD acquisitions (P < 0.001). High-resolution iodine and gadolinium maps, as well as virtual monoenergetic images, were calculated from the PCD data and enabled estimation of contrast agent concentration in the lumen without interference from the coronary stent.

CONCLUSION

HR spectral PCD imaging significantly improves coronary stent lumen visibility over dual-energy EID. When the PCD-HR data was reconstructed into standard voxel sizes (0.5 mm isotropic) the image noise decreased by 25% compared to SR acquisition of PCD. Both dual-energy systems were consistent in estimating contrast agent concentrations.

Impact of anti-charge sharing on the zero-frequency detective quantum efficiency of CdTe-based photon counting detector system: cascaded systems analysis and experimental validation

Inter-pixel communication and anti-charge sharing (ACS) technologies have been introduced to photon counting detector (PCD) systems to address the undesirable charge sharing problem. In addition to improving the energy resolution of PCD, ACS may also influence other aspects of PCD performance such as detector multiplicity (i.e. the number of pixels triggered by each interacted photon) and detective quantum efficiency (DQE). In this work, a theoretical model was developed to address how ACS impacts the multiplicity and zero-frequency DQE [DQE(0)] of PCD systems. The work focused on cadmium telluride (CdTe)-based PCD that often involves the generation and transport of K-fluorescence photons. Under the parallel cascaded systems analysis framework, the theory takes both photoelectric and scattering effects into account, and it also considers both the reabsorption and escape of photons. In a new theoretical treatment of ACS, it was considered as a modified version of the conventional single pixel (i.e. non-ACS) mode, but with reduced charge spreading distance and K-fluorescence travel distance. The proposed theoretical model does not require prior knowledge of the detailed ACS implementation method for each specific PCD, and its parameters can be experimentally determined using a radioisotope without invoking any Monte-Carlo simulation. After determining the model parameters, independent validation experiments were performed using a diagnostic x-ray tube and four different polychromatic beams (from 50 to 120 kVp). Both the theoretical and experimental results demonstrate that ACS increased the first and second moments of multiplicity for a majority of the x-ray energy and threshold levels tested, except when the threshold level was much lower than the x-ray energy level. However, ACS always improved DQE(0) at all energy and threshold levels tested.

Development of a chest digital tomosynthesis R/F system and implementation of low-dose GPU-accelerated compressed sensing (CS) image reconstruction

PURPOSE

This work describes the hardware and software developments of a prototype chest digital tomosynthesis (CDT) R/F system. The purpose of this study was to validate the developed system for its possible clinical application on low-dose chest tomosynthesis imaging.

METHODS

The prototype CDT R/F system was operated by carefully controlling the electromechanical subsystems through a synchronized interface. Once a command signal was delivered by the user, a tomosynthesis sweep started to acquire 81 projection views (PVs) in a limited angular range of ±20°. Among the full projection dataset of 81 images, several sets of 21 (quarter view) and 41 (half view) images with equally spaced angle steps were selected to represent a sparse view condition. GPU-accelerated and total-variation (TV) regularization strategy-based compressed sensing (CS) image reconstruction was implemented. The imaged objects were a flat-field using a copper filter to measure the noise power spectrum (NPS), a Catphan^® CTP682 quality assurance (QA) phantom to measure a task-based modulation transfer function (MTF_Task ) of three different cylinders' edge, and an anthropomorphic chest phantom with inserted lung nodules. The authors also verified the accelerated computing power over CPU programming by checking the elapsed time required for the CS method. The resultant absorbed and effective doses that were delivered to the chest phantom from two-view digital radiographic projections, helical computed tomography (CT), and the prototype CDT system were compared.

RESULTS

The prototype CDT system was successfully operated, showing little geometric error with fast rise and fall times of R/F x-ray pulse less than 2 and 10 ms, respectively. The in-plane NPS presented essential symmetric patterns as predicted by the central slice theorem. The NPS images from 21 PVs were provided quite different pattern against 41 and 81 PVs due to aliased noise. The voxel variance values which summed all NPS intensities were inversely proportional to the number of PVs, and the CS method gave much lower voxel variance by the factors of 3.97-6.43 and 2.28-3.36 compared to filtered backprojection (FBP) and 20 iterations of simultaneous algebraic reconstruction technique (SART). The spatial frequencies of the f₅₀ at which the MTF_Task reduced to 50% were 1.50, 1.55, and 1.67 cycles/mm for FBP, SART, and CS methods, respectively, in the case of Bone 20% cylinder using 41 views. A variety of ranges of TV reconstruction parameters were implemented during the CS method and we could observe that the NPS and MTF_Task preserved best when the regularization and TV smoothing parameters α and τ were in a range of 0.001-0.1. For the chest phantom data, the signal difference to noise ratios (SDNRs) were higher in the proposed CS scheme images than in the FBP and SART, showing the enhanced rate of 1.05-1.43 for half view imaging. The total averaged reconstruction time during 20 iterations of the CS scheme was 124.68 s, which could match-up a clinically feasible time (<3 min). This computing time represented an enhanced speed 386 times greater than CPU programming. The total amounts of estimated effective doses were 0.12, 0.53 (half view), and 2.56 mSv for two-view radiographs, the prototype CDT system, and helical CT, respectively, showing 4.49 times higher than conventional radiography and 4.83 times lower than a CT exam, respectively.

CONCLUSIONS

The current work describes the development and performance assessment of both hardware and software for tomosynthesis applications. The authors observed reasonable outcomes by showing a potential for low-dose application in CDT imaging using GPU acceleration.

Modeling and evaluation of a high-resolution CMOS detector for cone-beam CT of the extremities

PURPOSE

Quantitative assessment of trabecular bone microarchitecture in extremity cone-beam CT (CBCT) would benefit from the high spatial resolution, low electronic noise, and fast scan time provided by complementary metal-oxide semiconductor (CMOS) x-ray detectors. We investigate the performance of CMOS sensors in extremity CBCT, in particular with respect to potential advantages of thin (<0.7 mm) scintillators offering higher spatial resolution.

METHODS

A cascaded systems model of a CMOS x-ray detector incorporating the effects of CsI:Tl scintillator thickness was developed. Simulation studies were performed using nominal extremity CBCT acquisition protocols (90 kVp, 0.126 mAs/projection). A range of scintillator thickness (0.35-0.75 mm), pixel size (0.05-0.4 mm), focal spot size (0.05-0.7 mm), magnification (1.1-2.1), and dose (15-40 mGy) was considered. The detectability index was evaluated for both CMOS and a-Si:H flat-panel detector (FPD) configurations for a range of imaging tasks emphasizing spatial frequencies associated with feature size aobj. Experimental validation was performed on a CBCT test bench in the geometry of a compact orthopedic CBCT system (SAD = 43.1 cm, SDD = 56.0 cm, matching that of the Carestream OnSight 3D system). The test-bench studies involved a 0.3 mm focal spot x-ray source and two CMOS detectors (Dalsa Xineos-3030HR, 0.099 mm pixel pitch) - one with the standard CsI:Tl thickness of 0.7 mm (C700) and one with a custom 0.4 mm thick scintillator (C400). Measurements of modulation transfer function (MTF), detective quantum efficiency (DQE), and CBCT scans of a cadaveric knee (15 mGy) were obtained for each detector.

RESULTS

Optimal detectability for high-frequency tasks (feature size of ~0.06 mm, consistent with the size of trabeculae) was ~4× for the C700 CMOS detector compared to the a-Si:H FPD at nominal system geometry of extremity CBCT. This is due to ~5× lower electronic noise of a CMOS sensor, which enables input quantum-limited imaging at smaller pixel size. Optimal pixel size for high-frequency tasks was <0.1 mm for a CMOS, compared to ~0.14 mm for an a-Si:H FPD. For this fine pixel pitch, detectability of fine features could be improved by using a thinner scintillator to reduce light spread blur. A 22% increase in detectability of 0.06 mm features was found for the C400 configuration compared to C700. An improvement in the frequency at 50% modulation (f50 ) of MTF was measured, increasing from 1.8 lp/mm for C700 to 2.5 lp/mm for C400. The C400 configuration also achieved equivalent or better DQE as C700 for frequencies above ~2 mm-1 . Images of cadaver specimens confirmed improved visualization of trabeculae with the C400 sensor.

CONCLUSIONS

The small pixel size of CMOS detectors yields improved performance in high-resolution extremity CBCT compared to a-Si:H FPDs, particularly when coupled with a custom 0.4 mm thick scintillator. The results indicate that adoption of a CMOS detector in extremity CBCT can benefit applications in quantitative imaging of trabecular microstructure in humans.

Impact of bowtie filter and object position on the two-dimensional noise power spectrum of a clinical MDCT system

PURPOSE

Noise characteristics of clinical multidetector CT (MDCT) systems can be quantified by the noise power spectrum (NPS). Although the NPS of CT has been extensively studied in the past few decades, the joint impact of the bowtie filter and object position on the NPS has not been systematically investigated. This work studies the interplay of these two factors on the two dimensional (2D) local NPS of a clinical CT system that uses the filtered backprojection algorithm for image reconstruction.

METHODS

A generalized NPS model was developed to account for the impact of the bowtie filter and image object location in the scan field-of-view (SFOV). For a given bowtie filter, image object, and its location in the SFOV, the shape and rotational symmetries of the 2D local NPS were directly computed from the NPS model without going through the image reconstruction process. The obtained NPS was then compared with the measured NPSs from the reconstructed noise-only CT images in both numerical phantom simulation studies and experimental phantom studies using a clinical MDCT scanner. The shape and the associated symmetry of the 2D NPS were classified by borrowing the well-known atomic spectral symbols s, p, and d, which correspond to circular, dumbbell, and cloverleaf symmetries, respectively, of the wave function of electrons in an atom. Finally, simulated bar patterns were embedded into experimentally acquired noise backgrounds to demonstrate the impact of different NPS symmetries on the visual perception of the object.

RESULTS

(1) For a central region in a centered cylindrical object, an s-wave symmetry was always present in the NPS, no matter whether the bowtie filter was present or not. In contrast, for a peripheral region in a centered object, the symmetry of its NPS was highly dependent on the bowtie filter, and both p-wave symmetry and d-wave symmetry were observed in the NPS. (2) For a centered region-ofinterest (ROI) in an off-centered object, the symmetry of its NPS was found to be different from that of a peripheral ROI in the centered object, even when the physical positions of the two ROIs relative to the isocenter were the same. (3) The potential clinical impact of the highly anisotropic NPS, caused by the interplay of the bowtie filter and position of the image object, was highlighted in images of specific bar patterns oriented at different angles. The visual perception of the bar patterns was found to be strongly dependent on their orientation.

CONCLUSIONS

The NPS of CT depends strongly on the bowtie filter and object position. Even if the location of the ROI with respect to the isocenter is fixed, there can be different symmetries in the NPS, which depend on the object position and the size of the bowtie filter. For an isolated off-centered object, the NPS of its CT images cannot be represented by the NPS measured from a centered object.

[Evaluation for Basic Image Qualities Dependence on the Position in XYZ Directions and Acquisition Parameters of the Cone Beam CT for Angiography System with Flat Panel Detector]

The purpose of this study was to investigate the effect of the position in XYZ directions and acquisition parameters on the basic image qualities of for cone beam computed tomography (CBCT) in an angiography system with flat panel detector. The resolution property (modulation transfer function: MTF) and the noise property (Wiener spectrum: WS) of CBCT images in X-Y plane were measured with different acquisition parameters (scan matrix number and projection number) and the effect of the position in XYZ directions. The MTFs with 1024×1024 matrix were higher than those of 512×512 matrix and decreased in the peripheral areas due to the reduction of projection number. The highest and the lowest MTFs were measured at the X-ray tube side and on the detector side of the position in X-Y plane, respectively. The WS-doubled projection number showed about 50% lesser noise level. There were differences in the Wiener spectra (WS) at the position in XYZ directions. We conclude that the resolution and the noise property of CBCT image in X-Y plane showed dependences on the position in XYZ directions and acquisition parameters of the CBCT.

Technical Note: spektr 3.0-A computational tool for x-ray spectrum modeling and analysis

PURPOSE

A computational toolkit (spektr 3.0) has been developed to calculate x-ray spectra based on the tungsten anode spectral model using interpolating cubic splines (TASMICS) algorithm, updating previous work based on the tungsten anode spectral model using interpolating polynomials (TASMIP) spectral model. The toolkit includes a matlab (The Mathworks, Natick, MA) function library and improved user interface (UI) along with an optimization algorithm to match calculated beam quality with measurements.

METHODS

The spektr code generates x-ray spectra (photons/mm(2)/mAs at 100 cm from the source) using TASMICS as default (with TASMIP as an option) in 1 keV energy bins over beam energies 20-150 kV, extensible to 640 kV using the TASMICS spectra. An optimization tool was implemented to compute the added filtration (Al and W) that provides a best match between calculated and measured x-ray tube output (mGy/mAs or mR/mAs) for individual x-ray tubes that may differ from that assumed in TASMICS or TASMIP and to account for factors such as anode angle.

RESULTS

The median percent difference in photon counts for a TASMICS and TASMIP spectrum was 4.15% for tube potentials in the range 30-140 kV with the largest percentage difference arising in the low and high energy bins due to measurement errors in the empirically based TASMIP model and inaccurate polynomial fitting. The optimization tool reported a close agreement between measured and calculated spectra with a Pearson coefficient of 0.98.

CONCLUSIONS

The computational toolkit, spektr, has been updated to version 3.0, validated against measurements and existing models, and made available as open source code. Video tutorials for the spektr function library, UI, and optimization tool are available.

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.

Investigation on slice direction dependent detectability of volumetric cone beam CT images

We investigate the detection performance of transverse and longitudinal planes for various signal sizes (i.e., 1 mm to 8 mm diameter spheres) in cone beam computed tomography (CBCT) images. CBCT images are generated by computer simulation and images are reconstructed using an FDK algorithm. For each slice direction and signal size, a human observer study is conducted with a signal-known-exactly/background-known-exactly (SKE/BKE) binary detection task. The detection performance of human observers is compared with that of a channelized Hotelling observer (CHO). The detection performance of an ideal linear observer is also calculated using a CHO with Laguerre-Gauss (LG) channels. The detectability of high contrast small signals (i.e., up to 4-mm-diameter spheres) is higher in the longitudinal plane than the transverse plane. It is also shown that CHO performance correlates well with human observer performance in both transverse and longitudinal plane images.

Characteristic performance evaluation of a photon counting Si strip detector for low dose spectral breast CT imaging

PURPOSE

The possible clinical applications which can be performed using a newly developed detector depend on the detector's characteristic performance in a number of metrics including the dynamic range, resolution, uniformity, and stability. The authors have evaluated a prototype energy resolved fast photon counting x-ray detector based on a silicon (Si) strip sensor used in an edge-on geometry with an application specific integrated circuit to record the number of x-rays and their energies at high flux and fast frame rates. The investigated detector was integrated with a dedicated breast spectral computed tomography (CT) system to make use of the detector's high spatial and energy resolution and low noise performance under conditions suitable for clinical breast imaging. The aim of this article is to investigate the intrinsic characteristics of the detector, in terms of maximum output count rate, spatial and energy resolution, and noise performance of the imaging system.

METHODS

The maximum output count rate was obtained with a 50 W x-ray tube with a maximum continuous output of 50 kVp at 1.0 mA. A109Cd source, with a characteristic x-ray peak at 22 keV from Ag, was used to measure the energy resolution of the detector. The axial plane modulation transfer function (MTF) was measured using a 67 μm diameter tungsten wire. The two-dimensional (2D) noise power spectrum (NPS) was measured using flat field images and noise equivalent quanta (NEQ) were calculated using the MTF and NPS results. The image quality parameters were studied as a function of various radiation doses and reconstruction filters. The one-dimensional (1D) NPS was used to investigate the effect of electronic noise elimination by varying the minimum energy threshold.

RESULTS

A maximum output count rate of 100 million counts per second per square millimeter (cps/mm2) has been obtained (1 million cps per 100×100 μm pixel). The electrical noise floor was less than 4 keV. The energy resolution measured with the 22 keV photons from a 109Cd source was less than 9%. A reduction of image noise was shown in all the spatial frequencies in 1D NPS as a result of the elimination of the electronic noise. The spatial resolution was measured just above 5 line pairs per mm (lp/mm) where 10% of MTF corresponded to 5.4 mm(-1). The 2D NPS and NEQ shows a low noise floor and a linear dependence on dose. The reconstruction filter choice affected both of the MTF and NPS results, but had a weak effect on the NEQ.

CONCLUSIONS

The prototype energy resolved photon counting Si strip detector can offer superior imaging performance for dedicated breast CT as compared to a conventional energy-integrating detector due to its high output count rate, high spatial and energy resolution, and low noise characteristics, which are essential characteristics for spectral breast CT imaging.

Task-based measures of image quality and their relation to radiation dose and patient risk

The theory of task-based assessment of image quality is reviewed in the context of imaging with ionizing radiation, and objective figures of merit (FOMs) for image quality are summarized. The variation of the FOMs with the task, the observer and especially with the mean number of photons recorded in the image is discussed. Then various standard methods for specifying radiation dose are reviewed and related to the mean number of photons in the image and hence to image quality. Current knowledge of the relation between local radiation dose and the risk of various adverse effects is summarized, and some graphical depictions of the tradeoffs between image quality and risk are introduced. Then various dose-reduction strategies are discussed in terms of their effect on task-based measures of image quality.

Statistical model based iterative reconstruction (MBIR) in clinical CT systems: experimental assessment of noise performance

PURPOSE

To reduce radiation dose in CT imaging, the statistical model based iterative reconstruction (MBIR) method has been introduced for clinical use. Based on the principle of MBIR and its nonlinear nature, the noise performance of MBIR is expected to be different from that of the well-understood filtered backprojection (FBP) reconstruction method. The purpose of this work is to experimentally assess the unique noise characteristics of MBIR using a state-of-the-art clinical CT system.

METHODS

Three physical phantoms, including a water cylinder and two pediatric head phantoms, were scanned in axial scanning mode using a 64-slice CT scanner (Discovery CT750 HD, GE Healthcare, Waukesha, WI) at seven different mAs levels (5, 12.5, 25, 50, 100, 200, 300). At each mAs level, each phantom was repeatedly scanned 50 times to generate an image ensemble for noise analysis. Both the FBP method with a standard kernel and the MBIR method (Veo(®), GE Healthcare, Waukesha, WI) were used for CT image reconstruction. Three-dimensional (3D) noise power spectrum (NPS), two-dimensional (2D) NPS, and zero-dimensional NPS (noise variance) were assessed both globally and locally. Noise magnitude, noise spatial correlation, noise spatial uniformity and their dose dependence were examined for the two reconstruction methods.

RESULTS

(1) At each dose level and at each frequency, the magnitude of the NPS of MBIR was smaller than that of FBP. (2) While the shape of the NPS of FBP was dose-independent, the shape of the NPS of MBIR was strongly dose-dependent; lower dose lead to a "redder" NPS with a lower mean frequency value. (3) The noise standard deviation (σ) of MBIR and dose were found to be related through a power law of σ ∝ (dose)(-β) with the component β ≈ 0.25, which violated the classical σ ∝ (dose)(-0.5) power law in FBP. (4) With MBIR, noise reduction was most prominent for thin image slices. (5) MBIR lead to better noise spatial uniformity when compared with FBP. (6) A composite image generated from two MBIR images acquired at two different dose levels (D1 and D2) demonstrated lower noise than that of an image acquired at a dose level of D1+D2.

CONCLUSIONS

The noise characteristics of the MBIR method are significantly different from those of the FBP method. The well known tradeoff relationship between CT image noise and radiation dose has been modified by MBIR to establish a more gradual dependence of noise on dose. Additionally, some other CT noise properties that had been well understood based on the linear system theory have also been altered by MBIR. Clinical CT scan protocols that had been optimized based on the classical CT noise properties need to be carefully re-evaluated for systems equipped with MBIR in order to maximize the method's potential clinical benefits in dose reduction and/or in CT image quality improvement.

Comparative power law analysis of structured breast phantom and patient images in digital mammography and breast tomosynthesis

PURPOSE

This work characterizes three candidate mammography phantoms with structured background in terms of power law analysis in the low frequency region of the power spectrum for 2D (planar) mammography and digital breast tomosynthesis (DBT).

METHODS

The study was performed using three phantoms (spheres in water, Voxmam, and BR3D CIRS phantoms) on two DBT systems from two different vendors (Siemens Inspiration and Hologic Selenia Dimensions). Power spectra (PS) were calculated for planar projection, DBT projection, and reconstructed images and curve fitted in the low frequency region from 0.2 to 0.7 mm(-1) with a power law function characterized by an exponent β and magnitude κ. The influence of acquisition dose and tube voltage on the power law parameters was first explored. Then power law parameters were calculated from images acquired with the same anode∕filter combination and tube voltage for the three test objects, and compared with each other. Finally, PS curves for automatic exposure controlled acquisitions (anode∕filter combination and tube voltages selected by the systems based on the breast equivalent thickness of the test objects) were compared against PS analysis performed on patient data (for Siemens 80 and for Hologic 48 mammograms and DBT series). Dosimetric aspects of the three test objects were also examined.

RESULTS

The power law exponent (β) was found to be independent of acquisition dose for planar mammography but varied more for DBT projections of the sphere-phantom. Systematic increase of tube voltage did not affect β but decreased κ, both in planar and DBT projection phantom images. Power spectra of the BR3D phantom were closer to those of the patients than these of the Voxmam phantom; the Voxmam phantom gave high values of κ compared to the other phantoms and the patient series. The magnitude of the PS curves of the BR3D phantom was within the patient range but β was lower than the average patient value. Finally, PS magnitude for the sphere-phantom coincided with the patient curves for Siemens but was lower for the Hologic system. Close agreement of doses for all three phantoms with patient doses was found.

CONCLUSIONS

Power law parameters of the phantoms were close to those of the patients but no single phantom matched in terms of both magnitude (κ) and texture (β) for the x-ray systems in this work. PS analysis of structured phantoms is feasible and this methodology can be used to suggest improvements in phantom design.

Fundamental relationship between the noise properties of grating-based differential phase contrast CT and absorption CT: theoretical framework using a cascaded system model and experimental validation

PURPOSE

Using a grating interferometer, a conventional x-ray cone beam computed tomography (CT) data acquisition system can be used to simultaneously generate both conventional absorption CT (ACT) and differential phase contrast CT (DPC-CT) images from a single data acquisition. Since the two CT images were extracted from the same set of x-ray projections, it is expected that intrinsic relationships exist between the noise properties of the two contrast mechanisms. The purpose of this paper is to investigate these relationships.

METHODS

First, a theoretical framework was developed using a cascaded system model analysis to investigate the relationship between the noise power spectra (NPS) of DPC-CT and ACT. Based on the derived analytical expressions of the NPS, the relationship between the spatial-frequency-dependent noise equivalent quanta (NEQ) of DPC-CT and ACT was derived. From these fundamental relationships, the NPS and NEQ of the DPC-CT system can be derived from the corresponding ACT system or vice versa. To validate these theoretical relationships, a benchtop cone beam DPC-CT/ACT system was used to experimentally measure the modulation transfer function (MTF) and NPS of both DPC-CT and ACT. The measured three-dimensional (3D) MTF and NPS were then combined to generate the corresponding 3D NEQ.

RESULTS

Two fundamental relationships have been theoretically derived and experimentally validated for the NPS and NEQ of DPC-CT and ACT: (1) the 3D NPS of DPC-CT is quantitatively related to the corresponding 3D NPS of ACT by an inplane-only spatial-frequency-dependent factor 1∕f (2), the ratio of window functions applied to DPC-CT and ACT, and a numerical factor C(g) determined by the geometry and efficiency of the grating interferometer. Note that the frequency-dependent factor is independent of the frequency component f(z) perpendicular to the axial plane. (2) The 3D NEQ of DPC-CT is related to the corresponding 3D NEQ of ACT by an f (2) scaling factor and numerical factors that depend on both the attenuation and refraction properties of the image object, as well as C(g) and the MTF of the grating interferometer.

CONCLUSIONS

The performance of a DPC-CT system is intrinsically related to the corresponding ACT system. As long as the NPS and NEQ of an ACT system is known, the corresponding NPS and NEQ of the DPC-CT system can be readily estimated using additional characteristics of the grating interferometer.

Scaling-law for the energy dependence of anatomic power spectrum in dedicated breast CT

PURPOSE

To determine the x-ray photon energy dependence of the anatomic power spectrum of the breast when imaged with dedicated breast computed tomography (CT).

METHODS

A theoretical framework for scaling the empirically determined anatomic power spectrum at one x-ray photon energy to that at any given x-ray photon energy when imaged with dedicated breast CT was developed. Theory predicted that when the anatomic power spectrum is fitted with a power curve of the form k f(-β), where k and β are fit coefficients and f is spatial frequency, the exponent β would be independent of x-ray photon energy (E), and the amplitude k scales with the square of the difference in energy-dependent linear attenuation coefficients of fibroglandular and adipose tissues. Twenty mastectomy specimens based numerical phantoms that were previously imaged with a benchtop flat-panel cone-beam CT system were converted to 3D distribution of glandular weight fraction (f(g)) and were used to verify the theoretical findings. The 3D power spectrum was computed in terms of f(g) and after converting to linear attenuation coefficients at monoenergetic x-ray photon energies of 20-80 keV in 5 keV intervals. The 1D power spectra along the axes were extracted and fitted with a power curve of the form k f(-β). The energy dependence of k and β were analyzed.

RESULTS

For the 20 mastectomy specimen based numerical phantoms used in the study, the exponent β was found to be in the range of 2.34-2.42, depending on the axis of measurement. Numerical simulations agreed with the theoretical predictions that for a power-law anatomic spectrum of the form k f(-β), β was independent of E and k(E) = k(1)[μ(g)(E) - μ(a)(E)](2), where k(1) is a constant, and μ(g)(E) and μ(a)(E) represent the energy-dependent linear attenuation coefficients of fibroglandular and adipose tissues, respectively.

CONCLUSIONS

Numerical simulations confirmed the theoretical predictions that in dedicated breast CT, the spatial frequency dependence of the anatomic power spectrum will be independent of x-ray photon energy, and the amplitude of the anatomic power spectrum scales by the square of difference in linear attenuation coefficients of fibroglandular and adipose tissues.

To bin or not to bin? The effect of CT system limiting resolution on noise and detectability

We examine the noise advantages of having a computed tomography (CT) detector whose spatial resolution is significantly better (e.g. a factor of 2) than needed for a desired resolution in the reconstructed images. The effective resolution of detectors in x-ray CT is sometimes degraded by binning cells because the small cell size and fine sampling are not needed to achieve the desired resolution (e.g. with flat panel detectors). We studied the effect of the binning process on the noise in the reconstructed images and found that while the images in the absence of noise can be made identical for the native and the binned system, for the same system MTF in the presence of noise, the binned system always results in noisier reconstructed images. The effect of the increased noise in the reconstructed images on lesion detection is scale (frequency content) dependent with a larger difference between the high resolution and binned systems for imaging fine structure (small objects). We show simulated images reconstructed with both systems for representative objects and quantify the impact of the noise on the detection of the lesions based on mathematical observers. Through both subjective assessment of the reconstructed images and the quantification using mathematical observers, we show that for a CT system where the photon noise is dominant, higher resolution in the detectors leads to better noise performance in the reconstructed images at any resolution.

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.

A dedicated cone-beam CT system for musculoskeletal extremities imaging: design, optimization, and initial performance characterization

PURPOSE

This paper reports on the design and initial imaging performance of a dedicated cone-beam CT (CBCT) system for musculoskeletal (MSK) extremities. The system complements conventional CT and MR and offers a variety of potential clinical and logistical advantages that are likely to be of benefit to diagnosis, treatment planning, and assessment of therapy response in MSK radiology, orthopaedic surgery, and rheumatology.

METHODS

The scanner design incorporated a host of clinical requirements (e.g., ability to scan the weight-bearing knee in a natural stance) and was guided by theoretical and experimental analysis of image quality and dose. Such criteria identified the following basic scanner components and system configuration: a flat-panel detector (FPD, Varian 3030+, 0.194 mm pixels); and a low-power, fixed anode x-ray source with 0.5 mm focal spot (SourceRay XRS-125-7K-P, 0.875 kW) mounted on a retractable C-arm allowing for two scanning orientations with the capability for side entry, viz. a standing configuration for imaging of weight-bearing lower extremities and a sitting configuration for imaging of tensioned upper extremity and unloaded lower extremity. Theoretical modeling employed cascaded systems analysis of modulation transfer function (MTF) and detective quantum efficiency (DQE) computed as a function of system geometry, kVp and filtration, dose, source power, etc. Physical experimentation utilized an imaging bench simulating the scanner geometry for verification of theoretical results and investigation of other factors, such as antiscatter grid selection and 3D image quality in phantom and cadaver, including qualitative comparison to conventional CT.

RESULTS

Theoretical modeling and benchtop experimentation confirmed the basic suitability of the FPD and x-ray source mentioned above. Clinical requirements combined with analysis of MTF and DQE yielded the following system geometry: a -55 cm source-to-detector distance; 1.3 magnification; a 20 cm diameter bore (20 x 20 x 20 cm3 field of view); total acquisition arc of -240 degrees. The system MTF declines to 50% at -1.3 mm(-1) and to 10% at -2.7 mm(-1), consistent with sub-millimeter spatial resolution. Analysis of DQE suggested a nominal technique of 90 kVp (+0.3 mm Cu added filtration) to provide high imaging performance from -500 projections at less than -0.5 kW power, implying -6.4 mGy (0.064 mSv) for low-dose protocols and -15 mGy (0.15 mSv) for high-quality protocols. The experimental studies show improved image uniformity and contrast-to-noise ratio (without increase in dose) through incorporation of a custom 10:1 GR antiscatter grid. Cadaver images demonstrate exquisite bone detail, visualization of articular morphology, and soft-tissue visibility comparable to diagnostic CT (10-20 HU contrast resolution).

CONCLUSIONS

The results indicate that the proposed system will deliver volumetric images of the extremities with soft-tissue contrast resolution comparable to diagnostic CT and improved spatial resolution at potentially reduced dose. Cascaded systems analysis provided a useful basis for system design and optimization without costly repeated experimentation. A combined process of design specification, image quality analysis, clinical feedback, and revision yielded a prototype that is now awaiting clinical pilot studies. Potential advantages of the proposed system include reduced space and cost, imaging of load-bearing extremities, and combined volumetric imaging with real-time fluoroscopy and digital radiography.

Task-based modeling and optimization of a cone-beam CT scanner for musculoskeletal imaging

PURPOSE

This work applies a cascaded systems model for cone-beam CT imaging performance to the design and optimization of a system for musculoskeletal extremity imaging. The model provides a quantitative guide to the selection of system geometry, source and detector components, acquisition techniques, and reconstruction parameters.

METHODS

The model is based on cascaded systems analysis of the 3D noise-power spectrum (NPS) and noise-equivalent quanta (NEQ) combined with factors of system geometry (magnification, focal spot size, and scatter-to-primary ratio) and anatomical background clutter. The model was extended to task-based analysis of detectability index (d') for tasks ranging in contrast and frequency content, and d' was computed as a function of system magnification, detector pixel size, focal spot size, kVp, dose, electronic noise, voxel size, and reconstruction filter to examine trade-offs and optima among such factors in multivariate analysis. The model was tested quantitatively versus the measured NPS and qualitatively in cadaver images as a function of kVp, dose, pixel size, and reconstruction filter under conditions corresponding to the proposed scanner.

RESULTS

The analysis quantified trade-offs among factors of spatial resolution, noise, and dose. System magnification (M) was a critical design parameter with strong effect on spatial resolution, dose, and x-ray scatter, and a fairly robust optimum was identified at M ∼ 1.3 for the imaging tasks considered. The results suggested kVp selection in the range of ∼65-90 kVp, the lower end (65 kVp) maximizing subject contrast and the upper end maximizing NEQ (90 kVp). The analysis quantified fairly intuitive results-e.g., ∼0.1-0.2 mm pixel size (and a sharp reconstruction filter) optimal for high-frequency tasks (bone detail) compared to ∼0.4 mm pixel size (and a smooth reconstruction filter) for low-frequency (soft-tissue) tasks. This result suggests a specific protocol for 1 × 1 (full-resolution) projection data acquisition followed by full-resolution reconstruction with a sharp filter for high-frequency tasks along with 2 × 2 binning reconstruction with a smooth filter for low-frequency tasks. The analysis guided selection of specific source and detector components implemented on the proposed scanner. The analysis also quantified the potential benefits and points of diminishing return in focal spot size, reduced electronic noise, finer detector pixels, and low-dose limits of detectability. Theoretical results agreed quantitatively with the measured NPS and qualitatively with evaluation of cadaver images by a musculoskeletal radiologist.

CONCLUSIONS

A fairly comprehensive model for 3D imaging performance in cone-beam CT combines factors of quantum noise, system geometry, anatomical background, and imaging task. The analysis provided a valuable, quantitative guide to design, optimization, and technique selection for a musculoskeletal extremities imaging system under development.

Cone-Beam CT with a Flat-Panel Detector: From Image Science to Image-Guided Surgery

The development of large-area flat-panel x-ray detectors (FPDs) has spurred investigation in a spectrum of advanced medical imaging applications, including tomosynthesis and cone-beam CT (CBCT). Recent research has extended image quality metrics and theoretical models to such applications, providing a quantitative foundation for the assessment of imaging performance as well as a general framework for the design, optimization, and translation of such technologies to new applications. For example, cascaded systems models of Fourier domain metrics, such as noise-equivalent quanta (NEQ), have been extended to these modalities to describe the propagation of signal and noise through the image acquisition and reconstruction chain and to quantify the factors that govern spatial resolution, image noise, and detectability. Moreover, such models have demonstrated basic agreement with human observer performance for a broad range of imaging conditions and imaging tasks. These developments in image science have formed a foundation for the knowledgeable development and translation of CBCT to new applications in image-guided interventions - for example, CBCT implemented on a mobile surgical C-arm for intraoperative 3D imaging. The ability to acquire high-quality 3D images on demand during surgical intervention overcomes conventional limitations of surgical guidance in the context of preoperative images alone. A prototype mobile C-arm developed in academic-industry partnership demonstrates CBCT with low radiation dose, sub-mm spatial resolution, and soft-tissue visibility potentially approaching that of diagnostic CT. Integration of the 3D imaging system with real-time tracking, deformable registration, endoscopic video, and 3D visualization offers a promising addition to the surgical arsenal in interventions ranging from head-and-neck / skull base surgery to spine, orthopaedic, thoracic, and abdominal surgeries. Cadaver studies show the potential for significant boosts in surgical performance under CBCT guidance, and early clinical trials demonstrate feasibility, workflow, and image quality within the surgical theatre.

Local and global 3D noise power spectrum in cone-beam CT system with FDK reconstruction

PURPOSE

The authors examine the nonstationary noise behavior of a cone-beam CT system with FDK reconstruction.

METHODS

To investigate the nonstationary noise behavior, an analytical expression for the NPS of local volumes and an entire volume was derived and quantitatively compared to the NPS estimated from experimental air and water images.

RESULTS

The NPS of local volumes at different locations along the z-axis showed radial symmetry in the f(x)-f(y) plane and different missing cone regions in the f(z) direction depending on the tilt angle of rays through the local volumes. For local volumes away from the z-axis, the NPS of air and water images showed sharp transitions in the f(x)-f(y) and f(y)-f(z) planes and lack of radial symmetry in the f(x)-f(y) plane. These effects are mainly caused by varying magnification and different noise levels from view to view. In the NPS of the entire volume, the f(x)-f(y) plane showed radial symmetry because the nonstationary noise behaviors of local volumes were averaged out. The nonstationary sharp transitions were manifested as a high-frequency roll-off.

CONCLUSIONS

The results from noise power analysis for local volumes and an entire volume demonstrate the spatially varying noise behavior in the reconstructed cone-beam CT images.

The noise power spectrum in CT with direct fan beam reconstruction

The noise power spectrum (NPS) is a useful metric for understanding the noise content in images. To examine some unique properties of the NPS of fan beam CT, the authors derived an analytical expression for the NPS of fan beam CT and validated it with computer simulations. The nonstationary noise behavior of fan beam CT was examined by analyzing local regions and the entire field-of-view (FOV). This was performed for cases with uniform as well as nonuniform noise across the detector cells and across views. The simulated NPS from the entire FOV and local regions showed good agreement with the analytically derived NPS. The analysis shows that whereas the NPS of a large FOV in parallel beam CT (using a ramp filter) is proportional to frequency, the NPS with direct fan beam FBP reconstruction shows a high frequency roll off. Even in small regions, the fan beam NPS can show a sharp transition (discontinuity) at high frequencies. These effects are due to the variable magnification and therefore are more pronounced as the fan angle increases. For cases with nonuniform noise, the NPS can show the directional dependence and additional effects.

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.