Scatter Events Research Articles

SiPM signal readout for inter-crystal scatter event identification in PET detectors

In positron emission tomography (PET) with pixelated detectors, a significant number of annihilation photons interact with scintillation crystals through single or multiple Compton scattering events. When these partial energy depositions occur across multiple crystal elements, we call them inter-crystal scatter (ICS) events. ICS events lead to incorrect localization of the annihilation photons, thereby degrading the PET image contrast, spatial resolution, and lesion detectability. The accurate identification of ICS events is the first essential step to improve the quality of PET images by rejecting ICS events or recovering ICS events without losing PET sensitivity. In this study, we propose a novel silicon photomultiplier (SiPM) readout method to identify ICS events in one-to-one coupled PET detectors with a reduced number of data acquisition channels. For concept verification, we assembled a PET detector that consists of a 16-channel SiPM array and 4 4 lutetium oxyorthosilicate (LSO) array with a 3.2 mm crystal pitch. The proposed SiPM readout scheme serializes the 16 SiPM anode signals into four pulse train outputs encoded with four increasing time-delays in steps of 250 ns intervals. A Sum signal of the 16 SiPM anodes provides the timing information for time-of-flight measurement and a trigger signal for coincidence detection. A time-over-threshold (TOT) method was applied for obtaining the energy information followed by a subsequent TOT-to-energy calibration. We successfully identified the ICS events and determined their interacted positions and deposited energies by analyzing the digital pulses from the four pulse train output channels. The occurrence rate of ICS events was 10.85% for the 4 × 4 PET detector module with 3.2 mm-pitch LSO crystals. The PET detector yielded an energy resolution of 10.9 0.6% and coincidence timing resolution of 285 12 ps FWHM. We expect that the proposed method can be a useful solution for alleviating the readout burden of SiPM-based PET scanners with ICS event identification capability.

Identifying multiply scattered wavepaths in strongly scattering and dispersive media.

The ability to extract information from scattered waves is usually limited to singly scattered energy even if multiple scattering might occur in the medium. As a result, the information in arrival times of higher-order scattered events is underexplored. This information is extracted using fingerprinting theory. This theory has never previously been applied successfully to real measurements, particularly when the medium is dispersive. The theory is used to estimate the arrival times and scattering paths of multiply scattered waves in a thin sheet using an automated scheme in a dispersive medium by applying an additional dispersion compensation method. Estimated times and paths are compared with predictions based on a sequence of straight ray paths for each scattering event given the known scatterer locations. Additionally, numerical modelling is performed to verify the interpretations of the compensated data. Since the source also acts as a scatterer in these experiments, initially, the predictions and the numerical results did not conform to the experimental observations. By reformulating the theory and the processing scheme and adding a source scatterer in the modelling, it is shown that predictions of all observed scattering events are possible with both prediction methods, verifying that the methods are both effective and practically achievable.

Fisher information analysis of list-mode SPECT emission data for joint estimation of activity and attenuation distribution

The potential to perform attenuation and scatter compensation (ASC) in single-photon emission computed tomography (SPECT) without a separate transmission scan is highly significant. In this context, attenuation in SPECT is primarily due to Compton scattering, where the probability of Compton scatter is proportional to the attenuation coefficient of the tissue and the energy of the scattered photon and the scattering angle are related. Based on this premise, we investigated whether the SPECT emission data, including the scatter-window data, acquired in list-mode (LM) format and including the energy information, can be used to estimate the attenuation map. For this purpose, we propose a Fisher-information-based method that yields the Cramér–Rao bound (CRB) for the task of jointly estimating the activity and attenuation distribution using only the SPECT emission data. In the process, a path-based formalism to process the LM SPECT emission data, including the scattered-photon data, is proposed. The Fisher information method was implemented on NVIDIA graphics processing units (GPUs) for acceleration. The method was applied to quantify the information content of SPECT LM emission data, which contains up to first-order scattered events, in a simulated SPECT system with parameters modeling a clinical system using realistic computational studies with 2D digital synthetic and anthropomorphic phantoms. Experiments with anthropomorphic phantoms simulated myocardial perfusion and dopamine transporter (DaT)-Scan SPECT studies. The method was also applied to LM data containing up to second-order scatter for a synthetic phantom. The results show that the CRB obtained for the attenuation and activity coefficients was typically much lower than the true value of these coefficients. An increase in the number of detected photons yielded lower CRB for both the attenuation and activity coefficients. Further, we observed that systems with better energy resolution yielded a lower CRB for the attenuation coefficient. Overall, the results provide evidence that LM SPECT emission data, including the scatter-window data, contains information to jointly estimate the activity and attenuation coefficients.

Iterative Algorithms for Joint Scatter and Attenuation Estimation From Broken Ray Transform Data

The single-scatter approximation is fundamental in many tomographic imaging problems including x-ray scatter imaging and optical scatter imaging for certain media. In all cases, noisy measurements are affected by both local scatter events and nonlocal attenuation. Prior works focus on reconstructing one of two images: scatter density or total attenuation. However, both images are media specific and useful for object identification. Nonlocal effects of the attenuation image on the data are summarized by the broken ray transform (BRT). While analytic inversion formulas exist, poor conditioning of the inverse problem is only exacerbated by noisy measurements and sampling errors. This has motivated interest in the related star transforms incorporating BRT measurements from multiple source-detector pairs. However, all analytic methods operate on the log of the data. For media comprising regions with no scatter a new approach is required. We are the first to present a joint estimation algorithm based on Poisson data models for a single-scatter measurement geometry. Monotonic reduction of the log-likelihood function is guaranteed for our iterative algorithm while alternating image updates. We also present a fast algorithm for computing the discrete BRT forward operator. Our generalized approach can incorporate both transmission and scatter measurements from multiple source-detector pairs. Transmission measurements resolve low-frequency ambiguity in the joint image estimation problem, while multiple scatter measurements resolve the attenuation image. The benefits of joint estimation, over single-image estimation, vary with problem scaling. Our results quantify these benefits and should inform design of future acquisition systems.

Computational model for detector timing effects in Compton-camera based prompt-gamma imaging for proton radiotherapy

This paper describes a realistic simulation of a Compton-camera (CC) based prompt-gamma (PG) imaging system for proton range verification for a range of clinical dose rates, and its comparison to PG measured data with a pre-clinical CC. We used a Monte Carlo plus Detector Effects (MCDE) model to simulate the production of prompt gamma-rays (PG) and their energy depositions in the CC. With Monte Carlo, we simulated PG emission resulting from irradiation of a high density polyethylene phantom with a 150 MeV proton pencil beam at dose rates of 5.0 × 108, 2.6 × 109, and 4.6 × 109 p+ s−1. Realistic detector timing effects (e.g. delayed triggering time, event-coincidence, dead time, etc,) were added in post-processing to allow for flexible count rate variations. We acquired PG emission measurements with our pre-clinical CC during irradiation with a clinical 150 MeV proton pencil beam at the same dose rates. For simulations and measurements, three primary changes could be seen in the PG emission data as the dose rate increased: (1) reduction in the total number of detected events due to increased dead-time percentage; (2) increase in false-coincidence events (i.e. multiple PGs interacting, rather than a single PG scatter); and (3) loss of distinct PG emission peaks in the energy spectrum. We used the MCDE model to estimate the quality of our measured PG data, primarily with regards to true and false double-scatters and triple-scatters recorded by the CC. The simulation results showed that of the recorded double-scatter PG interactions 22%, 57%, and 70% were false double-scatters and for triple-scatter interactions 3%, 21%, and 35% were false events at 5.0 × 108, 2.6 × 109, and 4.6 × 109 p+ s−1, respectively. These false scatter events represent noise in the data, and the high percentage of these events in the data represents a major limitation in our ability to produce usable PG images with our prototype CC.

Limitation of Rayleigh sky model for bioinspired polarized skylight navigation in three-dimensional attitude determination

Insects such as desert ants and drosophilae can sense polarized skylight for navigation. Inspired by insects, many researchers have begun to study how to use skylight polarization patterns for attitude determination. The Rayleigh sky model has become the most widely used skylight polarization model for bioinspired polarized skylight navigation due to its simplicity and practicality. However, this is an ideal model considering only single Rayleigh scatter events, and the limitation of this model in bio-inspired attitude determination has not been paid much attention and lacks strict inference proof. To address this problem, the rotational and plane symmetry of the Rayleigh sky model are analyzed in detail, and it is theoretically proved that this model contains only single solar vector information, which contains only two independent scalar pieces of attitude information, so it is impossible to determine three Euler angles simultaneously in real-time. To further verify this conclusion, based on a designed hypothetical polarization camera, we discuss what conditions different three-dimensional attitudes must satisfy so that the polarization images taken at different 3D attitudes are the same; this indicates that multiple solutions will appear when only using the Rayleigh sky model to determine 3D attitude. In conclusion, due to its single solar vector information and the existence of multiple solutions, it is fully proved that 3D attitude cannot be determined in real time based only upon the Rayleigh sky model. Code is available at: https://github.com/HuajuLiang/HypotheticalPolarizationCamera

Exploiting stilbene’s scintillation anisotropy for neutron source localization

A technique for neutron source localization that exploits the scintillation anisotropy of stilbene is presented. The light output anisotropy of stilbene and constraints on the direction of recorded scatter events imposed by a threshold were used to estimate the direction of a 252Cf neutron source relative to the crystal axes of a stilbene crystal. The neutron source location was determined via triangulation using source directions estimated from multiple detectors. Two measurements that illustrate the efficacy of the technique are presented. The first measurement was designed with a favorable geometry for triangulation, which resulted in low uncertainties in the estimate of the neutron source location. The second measurement was designed with both favorable and unfavorable geometries for triangulation, which exhibited low and high uncertainties in the estimates of neutron source location, respectively. Pair-wise combinations of three detectors were used to estimate the source location for both measurements. The neutron source was localized with errors of 5.0, 8.9, and 3.9 cm for the first measurement with source detector-distance of 89, 90, and 88 cm for the three detectors. The neutron source was localized with errors of 2.8, 14.2, 6.2 cm for the second measurement with source detector-distances of 233, 123, and 130 cm.

Performance evaluation of HiPET, a high sensitivity and high resolution preclinical PET tomograph

HiPET is a recently developed prototype preclinical PET scanner dedicated to high sensitivity and high resolution molecular imaging. The HiPET system employs a phoswich depth of interaction (DOI) detector design, which also allows identification of the large majority of the cross layer crystal scatter (CLCS) events. This work evaluates its performance characteristics following the National Electrical Manufacturers Association (NEMA) NU4-2008 protocol.The HiPET consists of twenty flat panel type detectors arranged in two rings. The inner diameter is 160 mm and the axial field of view (FOV) is 104 mm. Each detector is comprised of two layers of phoswich scintillator crystal arrays, a tapered, pixelated glass lightguide and a multi anode photomultiplier tube (MAPMT). The front (gamma ray entrance) layer is a 48 × 48 pixelated cerium doped lutetium yttrium orthosilicate (LYSO) scintillator array with individual crystals measuring 1.01 × 1.01 × 6.1 mm. The back (towards the PMT) layer is a 32 × 32 pixelated bismuth germanate (BGO) scintillator array with individual crystals measuring 1.55 × 1.55 × 8.9 mm.For energy windows of 250–650 keV and 350–650 keV, the peak absolute sensitivity at the center of the FOV was 13.5% and 10.4% including CLCS events, and 11.8% and 8.9% excluding CLCS events, respectively. The average detector energy resolution derived by averaging the individual crystal spectra was 11.7% ± 1.4% for LYSO and 17.0% ± 1.4% for BGO. The 3D ordered-subsets expectation maximization (OSEM) reconstructed image of a point source in air, ranged from 0.73 mm to 1.19 mm, with an average value of 0.93 ± 0.09 mm at all measured locations. The peak noise equivalent count rate (NECR) and scatter fraction were 179 kcps at 12.4 MBq and 6.9% for the mouse-sized phantom, and 63 kcps at 11.3 MBq and 18.3% for the rat-sized phantom. For the NEMA image quality phantom, the uniformity was 5.8%, and the spillover ratios measured in the water- and air-filled cold region chambers were 0.047 and 0.044, respectively. The recovery coefficients (RC) ranged from 0.31 to 0.92.These results and in vivo evaluation demonstrate that the HiPET can achieve high quality molecular imaging for biomedical applications.

Sputter deposition of NiW films from a rotatable target

A series of NiW thin films was magnetron sputter deposited from a rotatable Ni80W20 target at Ar pressures ranging from 0.25 to 1.00 Pa on a matrix of (1 0 0) Si substrates. All films were synthesized in static deposition mode. The chosen setup allowed to study the influence of lateral substrate position and Ar pressure on film microstructure and elemental composition. All films deposited opposite the target were dense and Ni-rich, whereas films grown at extended positions from the target center had an open voided morphology and were W-rich. The differences in elemental composition with respect to the lateral position on the substrate carrier are attributed to the respective emission profiles of Ni and W and to scatter events of sputtered particles with Ar gas. To assess the influence of the target surface on final film properties, a planar Ni81W19 target was sputter eroded and investigated. The observed roughening of the target surface stems from the interplay between the high and low sputter yield of Ni and W, influencing the emission angles of both elements. Additional SRIM simulations of a thin Ni top-layer atop a W under-layer provided evidence for the existence of sputter yield amplification active within the studied system.

A low energy rare event search with the MAJORANA DEMONSTRATOR

The MAJORANA DEMONSTRATOR is sensitive to rare events near its energy threshold, including bosonic dark matter, solar axions, and lightly ionizing particles. In this analysis, a novel training set of low energy small-angle Compton scatter events is used to determine the efficiency of pulse shape analysis cuts, and we present updated bosonic dark matter and solar axion results from an 11.17 kg-y dataset using a 5 keV analysis threshold.

Optimization of Scatter Correction Method in Samarium-153 Single-photon Emission Computed Tomography using Triple-Energy Window: A Monte Carlo Simulation Study.

Purpose:In single-photon emission computed tomography imaging, the presence of scatter degrades image quality. The goal of this study is to optimize the main- and sub-energy windows for triple-energy window (TEW) method using Monte Carlo SImulating Medical Imaging Nuclear Detectors (SIMIND) code for samarium-153 (Sm-153) imaging.Materials and Methods:The comparison is based on the Monte Carlo simulation data with the results estimated using TEW method. Siemens Symbia gamma-camera equipped with low-energy high-resolution collimator was simulated for Sm-153 point source located in seven positions in water cylindrical phantom. Three different main-energy window widths (10%, 15%, and 20%) and three different sub-energy window widths (2, 4, and 6 keV) were evaluated. We compared the true scatter fraction determined by SIMIND and scatter fraction estimated using TEW scatter correction method at each position. In order to evaluate the image quality, we used the full width at half maximum (FWHM) computed on the PSF and image contrast using Jaszczak phantom.Results:The scatter fraction using TEW method is similar to the true scatter fraction for 20% of the main-energy window and 6 keV sub-energy windows. For these windows, the results show that the resolution and contrast were improved.Conclusion:TEW method could be a useful scatter correction method to remove the scatter event in the image for Sm-153 imaging.

Интеграция украинских мигрантов на польском рынке труда

Ukrainian citizens constitute an absolute majority among all categories of foreigners working in Poland. Their employment is largely characterized by over-representation in jobs with poor working and hyper-intensive conditions; by concentration in certain sectors of the economy; by a mismatch between work performed and available qualifications; and by weak legal protection. Limited access to stable employment, work in the «shadow economy», failure of employers to meet their financial obligations, and legal insecurity discriminate labor immigrants. The integration policy consists mainly of scattered events and does not form a coherent, long-term national strategy. While NGOs in several Polish voivodships and individual employers may tackle the issue of migrant integration, only a small proportion of legally employed Ukrainians are covered by integration measures. Short-term migrants and those employed in the “grey zone” are virtually “invisible” to both Polish society and the authorities.

Primary, scatter, and penetration characterizations of parallel-hole and pinhole collimators for I-123 SPECT

Multi-pinhole (MPH) collimators are known to provide better trade-off between sensitivity and resolution for preclinical, as well as for smaller regions in clinical SPECT imaging compared to conventional collimators. In addition to this geometric advantage, MPH plates typically offer better stopping power for penetration than the conventional collimators, which is especially relevant for I-123 imaging. The I-123 emits a series of high-energy (>300 keV, ~2.5% abundance) gamma photons in addition to the primary emission (159 keV, 83% abundance). Despite their low abundance, high-energy photons penetrate through a low-energy parallel-hole (LEHR) collimator much more readily than the 159 keV photons, resulting in large downscatter in the photopeak window. In this work, we investigate the primary, scatter, and penetration characteristics of a single pinhole collimator that is commonly used for I-123 thyroid imaging and our two MPH collimators designed for I-123 DaTscan imaging for Parkinson’s Disease, in comparison to three different parallel-hole collimators through a series of experiments and Monte Carlo simulations. The simulations of a point source and a digital human phantom with DaTscan activity distribution showed that our MPH collimators provide superior count performance in terms of high primary counts, low penetration, and low scatter counts compared to the parallel-hole and single pinhole collimators. For example, total scatter, multiple scatter, and collimator penetration events for the LEHR were 2.5, 7.6 and 14 times more than that of MPH within the 15% photopeak window. The total scatter fraction for LEHR was 56% where the largest contribution came from the high-energy scatter from the back compartments (31%). For the same energy window, the total scatter for MPH was 21% with only 1% scatter from the back compartments. We therefore anticipate that using MPH collimators, higher quality reconstructions can be obtained in a substantially shorter acquisition time for I-123 DaTscan and thyroid imaging.

Testing the Radiation Pattern of Meteor Radio Afterglow

AbstractRadio emission from meteors or meteor radio afterglows (MRAs) was first detected using the all‐sky imaging capabilities of the first station of the Long Wavelength Array (LWA1). In this work, we use the recently commissioned LWA Sevilleta (LWA‐SV) station along with the LWA1 to carry out coordinated observations. The combined all‐sky observations with LWA1 and LWA‐SV have coobserved 32 MRAs and 21 transmitter reflections from meteors (meteor scatter events), which are believed to be specular reflections from overdense trails. The flux density of the events observed by each station was measured from the all‐sky images. Triangulating the angular direction of events from each station gave the physical location and the distance of the event to each station. The luminosity of the events in each station was calculated using the flux distance relation for an isotropic source. The luminosity distribution for MRAs and meteor scatter events observed by each station shows a clear distinction between these two types of events as the ratio of luminosities is closer to unity for MRAs than the meteor scatter events. Furthermore, we find that MRAs follow an isotropic radiation pattern. This suggests, either a complete incoherent emission mechanism or an incoherent addition of coherently emitting small regions within the meteor trail.

Effect of CZT system characteristics on Compton scatter event recovery.

Improving 511 keV photon detection sensitivity is a common goal for positron emission tomography system designers. One attractive approach to increase sensitivity is recovering events that are normally rejected. The kinematics of Compton scattering can be used to recover the line of response through direction difference angle (DDA). The uncertainty of DDA is determined by the energy and spatial resolution of a system. In this work, we evaluated the performance of small animal CZT-based positron emission tomography systems with energy resolution of 1%, 4%, and 6% and different spatial resolution based on prior work for guiding new design efforts. Designs with energy resolution limited by counting statistics and by electronic noise were considered. The influence of modifying the conventional energy window and uncertainty of DDA was investigated. For a system with 4% energy resolution and limited by electronic noise, the figure of merit of noise equivalent count increases by 65% as the lower energy bound increases from 471 keV to 493 keV. If the system-wide energy resolution becomes worse than 4% of the full width half maximum at 511 keV, going to a pixel size finer than 1 mm has very limited effect in reducing total angular uncertainty. For a system with 1% energy resolution, as the spatial resolution improves from 1 mm to 0.5 mm, the contrast-to-noise ratio increases by 9%.

Time efficient scatter correction for time-of-flight PET: the immediate scatter approximation

Utilization of time-of-flight (TOF) information allows us to improve image quality and convergence rate in iterative PET image reconstruction. In order to obtain quantitatively correct images accurate scatter correction (SC) is required that accounts for the non-uniform distribution of scatter events over the TOF bins. However, existing simplified TOF-SC algorithms frequently exhibit limited accuracy while the currently accepted reference method—the TOF extension of the single scatter simulation approach (TOF-SSS)—is computationally demanding and can substantially slow down the reconstruction. In this paper we propose and evaluate a new accelerated TOF-SC algorithm in order to improve this situation. The key idea of the algorithm is the use of an immediate scatter approximation (ISA) for scatter time distribution calculation which speeds up estimation of the required TOF scatter by a factor of up to five in comparison to TOF-SSS. The proposed approach was evaluated in dedicated phantom measurements providing challenging high activity contrast conditions as well as in representative clinical patient data sets. Our results show that ISA is a viable alternative to TOF-SSS. The reconstructed images are in excellent quantitative agreement with those obtained with TOF-SSS while overall reconstruction time can be reduced by a factor of two in whole-body studies, even when using a listmode reconstruction not optimized for speed.

Correlation between the velocity scattering angle and product relative translational energy for SN2 reactions. Comparison of experiments and direct dynamics simulations

In previous research, direct dynamics simulations have been used to provide atomistic information for the Cl− + CH3I, F− + CH3I, and OH− + CH3I SN2 and proton transfer (PT) reactions. An important component of these simulations is comparison with molecular beam, ion imaging experiments. From the simulations, comparisons may be made with the product translational energy distributions determined from the experiments and the simulations give quite good agreement with these distributions. Though the experiments provide the in-plane angular distribution for a reaction’s product translational energy distribution, the number of direct dynamics trajectories for a simulation is not sufficient to compare with this angular distribution. In the work presented here, the average percentage product translational energy partitioning for forward scattered events with scattering angle θ = 0–90° and for backward scattered events with θ = 90–180° are compared for the experiments and simulations. Overall good agreement is found, with a maximum difference as high as 5–10%. Additional atomistic details, regarding the dynamics of these reactions, are provided by scatter plots of their product relative translational energy versus the scattering angle θ.

Compton PET: a simulation study for a PET module with novel geometry and machine learning for position decoding

This paper describes a simulation study of a positron emission tomography (PET) detector module that can reconstruct the kinematics of Compton scattering within the scintillator. We used a layer structure, with which we could recover the positions and energies for the multiple interactions of a gamma ray in the different layers. Using the Compton scattering formalism, the sequence of interactions can be estimated. The true first interaction position extracted in the Compton scattering will help minimize the degradation of the reconstructed image resolution caused by intercrystal scatter events. Because of the layer structure, this module also has readily available user-defined resolution for the depth of interaction. The semi-monolithic crystals enable high light collection efficiency and an energy resolution of ∼10% has been achieved in the simulation. We used machine learning to decode the gamma ray interaction locations, achieving an average spatial resolution of 0.40 mm. Our proposed detector design provides a pathway to increase the sensitivity of a system without affecting other key performance features.

An improved method for liquid density measurement using NaI(Tl) detector and low-strength source

The current work focuses on using low-strength radiation source and NaI(Tl) detector in liquid density measurement. To do this, the detector collimator is enlarged, which has resulted in the increased multiple scattered events. The scattering spectrum analysis method was applied successfully in removing the multiple scattered component. The obtained single scattering component corresponding to different liquid densities was used to build the linear standard curve which shows the linear dependence between of them with the R2 values of 0.968 and 0.982 for the collimator diameters of 3.0 and 9.5 cm, respectively. Besides, the Monte Carlo method was used in the current work which the obtained results also support for experimental ones.

Building blocks of a multi-layer PET with time sequence photon interaction discrimination and double Compton camera

Current PET detectors have a very low sensitivity, of the order of a few percent. One of the reasons is the fact that Compton interactions are rejected. If an event involves multiple Compton scattering and the total deposited energy lays within the photoelectric peak, then an energy-weighted centroid is the given output for the coordinates of the reconstructed interaction point. This introduces distortion in the final reconstructed image. The aim of our work is to prove that Compton events are a very rich source of additional information as one can improve the resolution of the detector and implicitly the final reconstructed image. This could be a real breakthrough for PET detector technology as one should be able to obtain better results with less patient radiation. Using a PET as a double Compton camera, by means of Compton cone matching i.e., Compton cones coming from the same event should be compatible, is applied to discard randoms, patient scattered events and also, to perform a correct matching among events with multiple coincidences. In order to fully benefit experimentally from Compton events using monolithic scintillators a multi-layer configuration is needed and a good time-of-flight resolution.