Search NASA⌕ Search

SEARCH · Search NASA

Results for “maximum likelihood expectation maximization (MLEM)”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

6 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Performance of neutron spectrum unfolding using deuterated liquid scintillator

The performance of the maximum likelihood expectation maximization method for unfolding neutron energy spectra using deuterated liquid scintillator is evaluated for future utilization with rare isotope beams. High-resolution neutron energy spectra as well as the detector response matrix were measured at the Edwards Accelerator Laboratory at Ohio University. Maximum-likelihood expectation maximization (MLEM) unfolded neutron spectra are compared with spectra from neutron time-of-flight. In conclusion, the effects of the MLEM stopping criteria and spectrum statistics are also investigated.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Point Source Likelihood Framework (psl) v1.0

PSL is a software framework for localizing radiation point sources in a 3D environment using a freely moving radiation detector system. The framework assumes a single point source of radiation is in the environment and uses known detector positions and orientations along with their radiation measurements to find the most likely location of the point source. This approach to the problem enables analyses that improve upon the capability of a standard Maximum Likelihood Expectation Maximization (MLEM) approach, such as setting a source detection threshold, finding a likelihood interval for the position of the source, finding a likelihood interval for the activity (intensity) of the source, determining whether more than one point source may be present, and generating a map of Minimum Detectable Activity (MDA) for a given source type. Applications include radiological and nuclear security, consequence management, decontamination and decommissioning of nuclear facilities, and radiological emergency response.

Joshi, Tenzing↗
DOE OSTI2022

Maximum Likelihood Spectrum Decomposition for Isotope Identification and Quantification

A spectral decomposition method has been implemented to identify and quantify isotopic source terms in high-resolution gamma-ray spectroscopy in static geometry and shielding scenarios. Monte Carlo simulations were used to build the response matrix of a shielded high-purity germanium detector monitoring an effluent stream with a Marinelli configuration. The decomposition technique was applied to a series of calibration spectra taken with the detector using a multi-nuclide standard. These results are compared with decay-corrected values from the calibration certificate. For most nuclei in the standard ( 241 Am, 109 Cd, 137 Cs, and 60 Co), the deviations from the certificate values were generally no more than 6% with a few outliers as high as 10%. Furthermore, for 57 Co, the radionuclide with the lowest activity, the deviations from the standard reached as high as 25%, driven by the meager statistics in the calibration spectra. In addition, a complete treatment of error propagation for the technique is presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Monte Carlo Simulation and Reconstruction: Assessment of Myocardial Perfusion Imaging of Tracer Dynamics With Cardiac Motion Due to Deformation and Respiration Using Gamma Camera With Continuous Acquisition

Purpose: Myocardial perfusion imaging (MPI) with single photon emission computed tomography (SPECT) is routinely used for stress testing in nuclear medicine. Recently, our group extended its potential going from 3D visual qualitative image analysis to 4D spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration and the estimated myocardial blood flow (MBF) and coronary flow reserve (CFR). However, the quality of reconstructed image is compromised due to cardiac deformation and respiration. The work presented here develops an algorithm that reconstructs the dynamic sequence of separate respiratory and cardiac phases and evaluates the algorithm with data simulated with a Monte Carlo simulation for the continuous image acquisition and processing with a slowly rotating SPECT camera. Methods: A clinically realistic Monte Carlo (MC) simulation is developed using the 4D Extended Cardiac Torso (XCAT) digital phantom with respiratory and cardiac motion to model continuous data acquisition of dynamic cardiac SPECT with slowly rotating gamma cameras by incorporating deformation and displacement of the myocardium due to cardiac and respiratory motion. We extended our previously developed 4D maximum-likelihood expectation-maximization (MLEM) reconstruction algorithm for a data set binned from a continuous list mode (LM) simulation with cardiac and respiratory information. Our spatiotemporal image reconstruction uses splines to explicitly model the temporal change of the tracer for each cardiac and respiratory gate that delineates the myocardial spatial position as the tracer washes in and out. Unlike in a fully list-mode data acquisition and reconstruction the accumulated photons are binned over a specific but very short time interval corresponding to each cardiac and respiratory gate. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it continuously deforms. These results are then compared with the conventional 4D spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. Mean Stabilized Activity (MSA), signal to noise ratio (SNR) and Bias for the myocardium activities for three different target-to-background ratios (TBRs) are evaluated. Dynamic quantitative indices such as wash-in (K1) and wash-out (k2) rates at each gate were also estimated. Results: The MSA and SNR are higher with higher TBRs while biases were improved with higher TBRs to less than 10%. The correlation between exhalation-inhalation sequence with the ground truth during respiratory cycle was excellent. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Estimated values of K1 and k2 were also consistent with the ground truth. Conclusion: The continuous image acquisition for dynamic scan using conventional two-head gamma cameras can provide valuable information for MPI. Our study demonstrated the viability of using a continuous image acquisition method on a widely used clinical two-head SPECT system. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Precise implementation of reconstruction algorithms, better segmentation techniques by generating images of different tissue types and background activity would improve the feasibility of the method in real clinical environment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

3D Source Reconstruction Using Coded Aperture Gamma-Ray Imaging

Recent measurements with coded-aperture imagers demonstrate material mass determination in a holdup setting to an accuracy within a few percent. This capability is of particular interest to the Surplus Plutonium Disposition (SPD) project, which aims to dilute and dispose of surplus plutonium oxide. Gamma-ray imagers can be used to determine holdup without interrupting normal operations. In this work, we examine techniques for 3D source localization and mass determination using gamma-ray imagers. Coded-aperture imagers provide excellent source localization within the 2D image plane; however, multiple imagers operating in tandem are necessary to identify source location in 3D space. A Maximum Likelihood Expectation-Maximization (MLEM) method for fitting detector mappings is a powerful tool for accomplishing this task. MLEM allows 3D source localization to be simultaneously constrained using multiple gamma-ray imagers by constructing the basis for the MLEM fit using detector mappings from different detector locations stitched together. Each of these basis points represents a singular response from a source in 3D space and is generated using Monte Carlo simulations of sources placed individually at different locations throughout the imager’s field of view. Additionally, implementing knowledge of the physical equipment in the simulations of the glovebox used for the SPD project incorporates attenuation effects that are needed to calculate material holdup.

Laminack, Alex↗
DOE OSTI2021

Characterization of a collimated neutron imager for low-rate fast neutron imaging

Spatial localization of special nuclear materials (SNM) via their neutron signatures amidst background requires knowledge of the background neutron environment or a means of separating a source from background based on low amounts of information. This requirement has created the need for characterizing the spatial distribution of the cosmogenic neutron background. Neutron scatter cameras have been developed and optimized for rapid detection of high activity sources, but have low imaging efficiency, making it difficult to use them to characterize low rate diffuse sources, such as the neutron background. The Low Intensity Neutron Imaging System (LINIS) is a collimated neutron imager that has been designed and optimized for imaging diffuse cosmogenic neutron background in the energy range of 0.5–15 MeV. LINIS operates using 16 liquid scintillation detectors shielded by ultra-high molecular weight polyethylene cylindrical collimators in a staggered orientation and rotates to 7 discrete positions, giving it roughly 2π sensitivity. Finally, LINIS has been characterized using (α, n) and fission neutron sources using two imaging techniques for neutron source localization, simple backprojection and Maximum Likelihood Expectation Maximization.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗