Noise Simulations For an Inverse-Geometry Volumetric CT System

This paper examines the noise performance of an inverse-geometry volumetric CT (IGCT) scanner through simulations. The IGCT system uses a large area scanned source and a smaller array of detectors to rapidly acquire volumetric data with negligible cone-beam artifacts. The first investigation compares the photon efficiency of the IGCT geometry to a 2D parallel ray system. The second investigation models the photon output of the IGCT source and calculates the expected noise. For the photon efficiency investigation. the same total number of photons was modeled in an IGCT acquisition and a comparable multi-slice 2D parallel ray acquisition. For both cases noise projections were simulated and the central axial slice reconstructed. In the second study. to investigate the noise in an IGCT system, the expected x-ray photon flux was modeled and projections simulated through ellipsoid phantoms. All simulations were compared to theoretical predictions. The results of the photon efficiency simulations verify that the IGCT geometry is as efficient in photon utilization as a 2D parallel ray geometry. For a 10 cm diameter 4 cm thick ellipsoid water phantom and for reasonable system parameters, the calculated standard deviation was approximately 15 HU at the center of the ellipsoid. For the same size phantom with maximum attenuation equivalent to 30 cm of water, the calculated noise was approximately 131 HU. The theoretical noise predictions for these objects were 15 HU and 112 HU respectively. These results predict acceptable noise levels for a system with a 0.16 second scan time and 12 lp/cm isotropic resolution

Three-Dimensional Reconstruction Algorithm for a Reverse-Geometry Volumetric CT System With a Large-Array Scanned Source

We have proposed a CT system design to rapidly produce volumetric images with negligible cone beam artifacts. The investigated system uses a large array scanned source with a smaller array of fast detectors. The x-ray source is electronically steered across a 2D target every few milliseconds as the system rotates. The proposed reconstruction algorithm for this system is a modified 3D filtered backprojection method. The data are rebinned into 2D parallel ray projections, most of which are tilted with respect to the axis of rotation. Each projection is filtered with a 2D kernel and backprojected onto the desired image matrix. To ensure adequate spatial resolution and low artifact level, we rebin the data onto an array that has sufficiently fine spatial and angular sampling. Due to finite sampling in the real system, some of the rebinned projections will be sparse, but we hypothesize that the large number of views will compensate for the data missing in a particular view. Preliminary results using simulated data with the expected discrete sampling of the source and detector arrays suggest that high resolution

Geometry Analysis of an Inverse-Geometry Volumetric CT System With Multiple Detector Arrays

An inverse-geometry volumetric CT (IGCT) system for imaging in a single fast rotation without cone-beam artifacts is being developed. It employs a large scanned source array and a smaller detector array. For a single-source/single-detector implementation, the FOV is limited to a fraction of the source size. Here we explore options to increase the FOV without increasing the source size by using multiple detectors spaced apart laterally to increase the range of radial distances sampled. We also look at multiple source array systems for faster scans. To properly reconstruct the FOV, Radon space must be sufficiently covered and sampled in a uniform manner. Optimal placement of the detectors relative to the source was determined analytically given system constraints (5cm detector width, 25cm source width, 45cm source-to-isocenter distance). For a 1x3 system (three detectors per source) detector spacing (DS) was 18deg and source-to-detector distances (SDD) were 113, 100 and 113cm to provide optimum Radon sampling and a FOV of 44cm. For multiple-source systems, maximum angular spacing between sources cannot exceed 125deg since detectors corresponding to one source cannot be occluded by a second source. Therefore, for 2x3 and 3x3 systems using the above DS and SDD, optimum spacing between sources is 115deg and 61deg respectively, requiring minimum scan rotations of 115deg and 107deg. Also, a 3x3 system can be much faster for full 360deg dataset scans than a 2x3 system (120deg vs. 245deg). We found that a significantly increased FOV can be achieved while maintaining uniform radial sampling as well as a substantial reduction in scan time using several different geometries. Further multi-parameter optimization is underway

Histogram Analysis of ADC in Brain Tumor Patients

At various stage of progression, most brain tumors are not homogenous. In this presentation, we retrospectively studied the distribution of ADC values inside tumor volume during the course of tumor treatment and progression for a selective group of patients who underwent an anti-VEGF trial. Complete MRI studies were obtained for this selected group of patients including pre- and multiple follow-up, post-treatment imaging studies. In each MRI imaging study, multiple scan series were obtained as a standard protocol which includes T1, T2, T1-post contrast, FLAIR and DTI derived images (ADC, FA etc.) for each visit. All scan series (T1, T2, FLAIR, post-contrast T1) were registered to the corresponding DTI scan at patient\u27s first visit. Conventionally, hyper-intensity regions on T1-post contrast images are believed to represent the core tumor region while regions highlighted by FLAIR may overestimate tumor size. Thus we annotated tumor regions on the T1-post contrast scans and ADC intensity values for pixels were extracted inside tumor regions as defined on T1-post scans. We fit a mixture Gaussian (MG) model for the extracted pixels using the Expectation-Maximization (EM) algorithm, which produced a set of parameters (mean, various and mixture coefficients) for the MG model. This procedure was performed for each visits resulting in a series of GM parameters. We studied the parameters fitted for ADC and see if they can be used as indicators for tumor progression. Additionally, we studied the ADC characteristics in the peri-tumoral region as identified by hyper-intensity on FLAIR scans. The results show that ADC histogram analysis of the tumor region supports the two compartment model that suggests the low ADC value subregion corresponding to densely packed cancer cell while the higher ADC value region corresponding to a mixture of viable and necrotic cells with superimposed edema. Careful studies of the composition and relative volume of the two compartments in tumor region may provide some insights in the early assessment of tumor response to therapy for recurrence brain cancer patients

Spectral optimization using fast kV switching and filtration for photon counting CT with realistic detector responses: a simulation study.

PurposePhoton counting CT (PCCT) provides spectral measurements for material decomposition. However, the image noise (at a fixed dose) depends on the source spectrum. Our study investigates the potential benefits from spectral optimization using fast kV switching and filtration to reduce noise in material decomposition.ApproachThe effect of the input spectra on noise performance in both two-basis material decomposition and three-basis material decomposition was compared using Cramer-Rao lower bound analysis in the projection domain and in a digital phantom study in the image domain. The fluences of different spectra were normalized using the CT dose index to maintain constant dose levels. Four detector response models based on Si or CdTe were included in the analysis.ResultsFor single kV scans, kV selection can be optimized based on the imaging task and object size. Furthermore, our results suggest that noise in material decomposition can be substantially reduced with fast kV switching. For two-material decomposition, fast kV switching reduces the standard deviation (SD) by ∼10% . For three-material decomposition, greater noise reduction in material images was found with fast kV switching (26.2% for calcium and 25.8% for iodine, in terms of SD), which suggests that challenging tasks benefit more from the richer spectral information provided by fast kV switching.ConclusionsThe performance of PCCT in material decomposition can be improved by optimizing source spectrum settings. Task-specific tube voltages can be selected for single kV scans. Also, our results demonstrate that utilizing fast kV switching can substantially reduce the noise in material decomposition for both two- and three-material decompositions, and a fixed Gd filter can further enhance such improvements for two-material decomposition

Study of Increased Radiation When an X-ray Tube is Placed in a Strong Magnetic Field

When a fixed anode x-ray tube is placed in a magnetic field (B) that is parallel to the anode-cathode axis, the x-ray exposure increases with increasing B. It was hypothesized that the increase was caused by backscattered electrons which were constrained by B and reaccelerated by the electric field onto the x-ray tube target. We performed computer simulations and physical experiments to study the behavior of the backscattered electrons in a magnetic field, and their effects on the radiation output, x-ray spectrum, and off-focal radiation. A Monte Carlo program (EGS4) was used to generate the combined energy and angular distribution of the backscattered electrons. The electron trajectories were traced and their landing locations back on the anode were calculated. Radiation emission from each point was modeled with published data (IPEM Report 78), and thus the exposure rate and x-ray spectrum with the contribution of backscattered electrons could be predicted. The point spread function for a pencil beam of electrons was generated and then convolved with the density map of primary electrons incident on the anode as simulated with a finite element program (Opera-3d, Vector Fields, UK). The total spatial distribution of x-ray emission could then be calculated. Simulations showed that for an x-ray tube working at 65 kV, about 54% of the electrons incident on the target were backscattered. In a magnetic field of 0.5 T, although the exposure would be increased by 33%, only a small fraction of the backscattered electrons landed within the focal spot area. The x-ray spectrum was slightly shifted to lower energies and the half value layer (HVL) was reduced by about 6%. Measurements of the exposure rate, half value layer and focal spot distribution were acquired as functions of B. Good agreement was observed between experimental data and simulation results. The wide spatial distribution of secondary x-ray emission can degrade the MTF of the x-ray system at low spatial frequencies for B {le} 0.5 T