Development and tuning of models for accurate simulation of CT spatial resolution using CatSim.

Objective. We sought to systematically evaluate CatSim's ability to accurately simulate the spatial resolution produced by a typical 64-detector-row clinical CT scanner in the projection and image domains, over the range of clinically used x-ray techniques.Approach.Using a 64-detector-row clinical scanner, we scanned two phantoms designed to evaluate spatial resolution in the projection and image domains. These empirical scans were performed over the standard clinically used range of x-ray techniques (kV, and mA). We extracted projection data from the scanner, and we reconstructed images. For the CatSim simulations, we developed digital phantoms to represent the phantoms used in the empirical scans. We developed a new, realistic model for the x-ray source focal spot, and we empirically tuned a published model for the x-ray detector temporal response. We applied these phantoms and models to simulate scans equivalent to the empirical scans, and we reconstructed the simulated projections using the same methods used for the empirical scans. For the empirical and simulated scans, we qualitatively and quantitatively compared the projection-domain and image-domain point-spread functions (PSFs) as well as the image-domain modulation transfer functions. We reported four quantitative metrics and the percent error between the empirical and simulated results.Main Results.Qualitatively, the PSFs matched well in both the projection and image domains. Quantitatively, all four metrics generally agreed well, with most of the average errors substantially less than 5% for all x-ray techniques. Although the errors tended to increase with decreasing kV, we found that the CatSim simulations agreed with the empirical scans within limits required for the anticipated applications of CatSim.Significance.The new focal spot model and the new detector temporal response model are significant contributions to CatSim because they enabled achieving the desired level of agreement between empirical and simulated results. With these new models and this validation, CatSim users can be confident that the spatial resolution represented by simulations faithfully represents results that would be obtained by a real scanner, within reasonable, known limits. Furthermore, users of CatSim can vary parameters including but not limited to system geometry, focal spot size/shape and detector parameters, beyond the values available in physical scanners, and be confident in the results. Therefore, CatSim can be used to explore new hardware designs as well as new scanning and reconstruction methods, thus enabling acceleration of improved CT scan capabilities.

Semiempirical, parameterized spectrum estimation for x-ray computed tomography.

PURPOSE: To develop a tool to produce accurate, well-validated x-ray spectra for standalone use or for use in an open-access x-ray/CT simulation tool. Spectrum models will be developed for tube voltages in the range of 80 kVp through 140 kVp and for anode takeoff angles in the range of 5° to 9°. METHODS: Spectra were initialized based on physics models, then refined using empirical measurements, as follows. A new spectrum-parameterization method was developed, including 13 spline knots to represent the bremsstrahlung component and 4 values to represent characteristic lines. Initial spectra at 80, 100, 120, and 140 kVp and at takeoff angles from 5° to 9° were produced using physics-based spectrum estimation tools XSPECT and SpekPy. Empirical experiments were systematically designed with careful selection of attenuator materials and thicknesses, and by reducing measurement contamination from scatter to <1%. Measurements were made on a 64-row CT scanner using the scanner's detector and using multiple layers of polymethylmethacrylate (PMMA), aluminum, titanium, tin, and neodymium. Measurements were made at 80, 100, 120, and 140 kVp and covering the entire 64-row detector (takeoff angles from 5° to 9°); a total of 6,144 unique measurements were made. After accounting for the detector's energy response, parameterized representations of the initial spectra were refined for best agreement with measurements using two proposed optimization schemes: based on modulation and based on gradient descent. X-ray transmission errors were computed for measurements vs calculations using the nonoptimized and optimized spectra. Half-value, tenth-value, and hundredth-value layers for PMMA, Al, and Ti were calculated. RESULTS: Spectra before and after parameterization were in excellent agreement (e.g., R2 values of 0.995 and 0.997). Empirical measurements produced smoothly varying curves with x-ray transmission covering a range of up to 3.5 orders of magnitude. Spectra from the two optimization schemes, compared with the unoptimized physic-based spectra, each improved agreement with measurements by twofold through tenfold, for both postlog transmission data and for fractional value layers. CONCLUSION: The resulting well-validated spectra are appropriate for use in the open-access x-ray/CT simulator under development, the x-ray-based Cancer Imaging Toolkit (XCIST), or for standalone use. These spectra can be readily interpolated to produce spectra at arbitrary kVps over the range of 80 to 140 kVp and arbitrary takeoff angles over the range of 5° to 9°. Furthermore, interpolated spectra over these ranges can be obtained by applying the standalone Matlab function available at https://github.com/xcist/documentation/blob/master/XCISTspectrum.m.