Statistical reconstruction for cosmic ray muon tomography.
Highly penetrating cosmic ray muons constantly shower the earth at a rate of about 1 muon per cm2 per minute. We have developed a technique which exploits the multiple Coulomb scattering of these particles to perform nondestructive inspection without the use of artificial radiation. In prior work [1]-[3], we have described heuristic methods for processing muon data to create reconstructed images. In this paper, we present a maximum likelihood/expectation maximization tomographic reconstruction algorithm designed for the technique. This algorithm borrows much from techniques used in medical imaging, particularly emission tomography, but the statistics of muon scattering dictates differences. We describe the statistical model for multiple scattering, derive the reconstruction algorithm, and present simulated examples. We also propose methods to improve the robustness of the algorithm to experimental errors and events departing from the statistical model.
Schultz
Larry J
LJ
Los Alamos National Laboratory, Los Alamos, NM 87545, USA. schultz@lanl.gov
Blanpied
Gary S
GS
Borozdin
Konstantin N
KN
Fraser
Andrew M
AM
Hengartner
Nicolas W
NW
Klimenko
Alexei V
AV
Morris
Christopher L
CL
Orum
Chris
C
Sossong
Michael J
MJ
D000465:Algorithms; D003198:Computer Simulation; D003359:Cosmic Radiation; D003627:Data Interpretation, Statistical; D007089:Image Enhancement / Q000379:methods; D007090:Image Interpretation, Computer-Assisted / Q000379:methods*; D015233:Models, Statistical; D015203:Reproducibility of Results; D012680:Sensitivity and Specificity; D014054:Tomography / Q000379:methods*