  -Purpose: 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. ralvarez@aprendtech.com Networks trained with
high noise data have both large bias and large variance. Conclusion: With a
properly chosen level of added training data noise, a neural network estimator
for photon counting data with pileup can have variance close to the CRLB with
negligible bias.
