Near optimal neural network estimator for spectral x-ray photon counting data with pileup

Abstract

A neural network estimator to process x-ray spectral measurements from photon counting detectors with pileup. The estimator is used with an expansion of the attenuation coefficient as a linear combination of functions of energy multiplied by coefficients that depend on the material composition at points within the object [R.E. Alvarez and A. Macovski, Phys. Med. Biol., 1976, 733-744]. The estimator computes the line integrals of the coefficients from measurements with different spectra. Neural network estimators are trained with measurements of a calibration phantom with the clinical x-ray system. One estimator uses low noise training data and another network is trained with data computed by adding random noise to the low noise data. The performance of the estimators is compared to each other and to the Cramer-Rao lower bound (CRLB). Methods: The estimator performance is measured using a Monte Carlo simulation with an idealized model of a photon counting detector that includes only pileup and quantum noise. Transmitted x-ray spectra are computed for a calibration phantom. The transmitted spectra are used to compute random data for photon counting detectors with pileup. Detectors with small and large dead times are considered. Neural network training data with extremely low noise are computed by averaging the random detected data with pileup for a large numbers of exposures of the phantom. Each exposure is equivalent to a projection image or one projection of a computed tomography scan. Training data with high noise are computed by using data from one exposure. Finally, training data are computed by adding random data to the low noise data. The added random data are multivariate normal with zero mean and covariance equal to the sample covariance of data for an object with properly chosen attenuation. To test the estimators, random data are computed for different thicknesses of three test objects with different compositions. These are used as inputs to the neural network estimators. The mean squared errors (MSE), variance and square of the bias of the neural networks' outputs with the random object data are each compared to the CRLB. Results: The MSE for a network trained with low noise data and added noise is close to the CRLB for both the low and high pileup cases. Networks trained with very low noise data have low bias but large variance for both pileup cases.