Assignment 1: Histogram, Descriptive Statistics, and Scatter Plots
No abstract available
Newman, Jessica A.↗
SEARCH · Search NASA
Results for “histogram”
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.
High-Performance Zonal Histogramming on Large-Scale Geospatial Rasters Using GPUs and GPU-Accelerated Clusters
Explore the source record for details and available documents.
Zhang, Jianting↗
Gamma ray spectroscopy employing divalent europium-doped alkaline earth halides and digital readout for accurate histogramming
A scintillator radiation detector system according to one embodiment includes a scintillator; and a processing device for processing pulse traces corresponding to light pulses from the scintillator, wherein pulse digitization is used to improve energy resolution of the system. A scintillator radiation detector system according to another embodiment includes a processing device for fitting digitized scintillation waveforms to an algorithm based on identifying rise and decay times and performing a direct integration of fit parameters. A method according to yet another embodiment includes processing pulse traces corresponding to light pulses from a scintillator, wherein pulse digitization is used to improve energy resolution of the system. A method in a further embodiment includes fitting digitized scintillation waveforms to an algorithm based on identifying rise and decay times; and performing a direct integration of fit parameters. Additional systems and methods are also presented.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Gamma ray spectroscopy employing divalent europium-doped alkaline earth halides and digital readout for accurate histogramming
According to one embodiment, a scintillator radiation detector system includes a scintillator, and a processing device for processing pulse traces corresponding to light pulses from the scintillator, where the processing device is configured to: process each pulse trace over at least two temporal windows and to use pulse digitization to improve energy resolution of the system. According to another embodiment, a scintillator radiation detector system includes a processing device configured to: fit digitized scintillation waveforms to an algorithm, perform a direct integration of fit parameters, process multiple integration windows for each digitized scintillation waveform to determine a correction factor, and apply the correction factor to each digitized scintillation waveform.
36 MATERIALS SCIENCE↗
A Formula for the Maximum Number of Bins for Approximately Smooth Histograms with Application to the Pearson Chi2 Test and Related Topics.
Abstract not provided.
Dalbey, Keith↗
Quantifying cell migration distance as a contributing factor to the development of rectal toxicity after prostate radiotherapy
Purpose: Spatial information is usually neglected in mathematical models of radiation-induced toxicity. In the presence of inhomogeneous dose distributions produced by intensity modulated radiation therapy (IMRT) and volumetric arc therapy, this may be a limitation. We present a model to quantify the spatial characteristics of the dose distribution on the rectum through the quantification of the distribution of distances between dose points on the surface of the rectum in three-dimensions. The method allows us to evaluate the hypothesis that distances between lower and higher dose regions on the rectum influence radiation damage repair due to the migration of normal cells into damaged areas, and consequently, the development of radiation-induced toxicity in patients treated with radiation for prostate cancer. Methods: We present a method to compute distances between dose points on the surface of the rectum in three dimensions (3D) and to generate distance maps representing the distances between specific dose regions on the rectum. We introduce the concept of the distance dose surface histogram (DDSH), which is computed from the distance maps. The DDSH is a 2D histogram of rectum area on a grid defined by pairwise combinations of dose and distance. Each bin in the DDSH quantifies the area of the rectum exposed to a given dose and at a given distance from other another dose region on the rectum. By summing across the columns and rows of the DDSH, we can generate the dose surface histogram (DSH) and distance surface histogram (DiSH) for a particular dose region. The DiSH is a marginal histogram showing the distribution of distances for the dose points in a specific dose region from another region. We computed the DDSH, DiSH, and DSH for 33 patients treated with IMRT for prostate cancer, nine of whom developed late Grade 2 or higher late rectal toxicity. Results: We show how even though the total area of the rectum exposed to a given dose may be the same for different patients, the distribution of distances for points receiving that dose can be different depending on the shape and contiguity of the region(s) containing those dose points. We also show that area of the rectum in the region receiving more than 75 Gy and at a distance of 16 to 22 mm from the 50 Gy isodose line was significantly correlated to the development of toxicity (p = 0.004, two sided t-test). This suggests that examining the distance distribution of points in specific dose regions could provide valuable additional information in predicting the risk of a patient developing radiation-induced rectal toxicity. Conclusions: We present a new method to quantify the spatial distribution of points in a given region relative to other regions on the rectum. The method provides a means to evaluate the hypothesis that distances between lower and higher dose regions on the rectum influence radiation damage repair due to the migration of normal cells into damaged areas, and may be a contributing factor to the development of radiation-induced toxicity in patients treated with radiation for prostate cancer.
62 RADIOLOGY AND NUCLEAR MEDICINE↗
A Bayesian Modeling Approach for Estimation of a Shape-Free Groundwater Age Distribution using Multiple Tracers
The mixing of groundwaters with different ages in aquifers, groundwater age is more appropriately represented by a distribution rather than a scalar number. To infer a groundwater age distribution from environmental tracers, a mathematical form is often assumed for the shape of the distribution and the parameters of the mathematical distribution are estimated using deterministic or stochastic inverse methods. We found that the prescription of the mathematical form limits the exploration of the age distribution to the shapes that can be described by the selected distribution. In this paper, the use of freeform histograms as groundwater age distributions is evaluated. A Bayesian Markov Chain Monte Carlo approach is used to estimate the fraction of groundwater in each histogram bin. This method was able to capture the shape of a hypothetical gamma distribution from the concentrations of four age tracers. The number of bins that can be considered in this approach is limited based on the number of tracers available. The histogram method was also tested on tracer data sets from Holten (The Netherlands; 3 H, 3 He, 85 Kr, 39 Ar) and the La Selva Biological Station (Costa-Rica; SF 6 , CFCs, 3 H, 4 He and 14 C), and compared to a number of mathematical forms. According to standard Bayesian measures of model goodness, the best mathematical distribution performs better than the histogram distributions in terms of the ability to capture the observed tracer data relative to their complexity. Among the histogram distributions, the four bin histogram performs better in most of the cases. The Monte Carlo simulations showed strong correlations in the posterior estimates of bin contributions, indicating that these bins cannot be well constrained using the available age tracers. The fact that mathematical forms overall perform better than the freeform histogram does not undermine the benefit of the freeform approach, especially for the cases where a larger amount of observed data is available and when the real groundwater distribution is more complex than can be represented by simple mathematical forms.
54 ENVIRONMENTAL SCIENCES↗
Boundary detection evaluation
Illustrative embodiments are directed to a method and apparatus for evaluating boundary detection in an image. A processed image is received, wherein a detected boundary of an image of an object is identified in the processed image. A Radon transform is applied to the processed image for a plurality of angles to form a processed image histogram corresponding to the detected boundary for each of the plurality of angles. The processed image histogram for each of the plurality of angles and a corresponding ground truth histogram for each of the plurality of angles is normalized to provide a normalized processed image histogram and a normalized ground truth histogram for each of the plurality of angles, wherein the ground truth histogram corresponds to a ground truth boundary of the object for a corresponding angle. An indication of the edges of the normalized processed image histogram for each of the plurality of angles is plotted to form a boundary detection evaluation visualization.
Wantuch, Andrew C.↗
Boundary detection evaluation
Illustrative embodiments are directed to a method and apparatus for evaluating boundary detection in an image. A processed image is received, wherein a detected boundary of an image of an object is identified in the processed image. A Radon transform is applied to the processed image for a plurality of angles to form a processed image histogram corresponding to the detected boundary for each of the plurality of angles. The processed image histogram for each of the plurality of angles and a corresponding ground truth histogram for each of the plurality of angles is normalized to provide a normalized processed image histogram and a normalized ground truth histogram for each of the plurality of angles, wherein the ground truth histogram corresponds to a ground truth boundary of the object for a corresponding angle. An indication of the edges of the normalized processed image histogram for each of the plurality of angles is plotted to form a boundary detection evaluation visualization.
Wantuch, Andrew C.↗
GPU Acceleration of Kernel Density Estimators in Monte Carlo Neutron Transport Simulations
Kernel Density Estimators (KDEs) have recently been explored for use in Monte Carlo radiation transport simulations as an alternative to histogram tallies for capturing spatially-resolved quantities such as neutron flux and reaction rates. Histogram tallies suffer from large uncertainties when detailed spatial resolution of the quantity of interest is required; increasing the resolution of the underlying mesh increases the uncertainty in each histogram bin. KDEs obtain estimates of underlying densities on user-defined points, with the uncertainty at those points being independent of the spatial resolution desired since collisions and particle-track segments are allowed to contribute to the score at multiple tally points. Thus, KDEs show potential for obtaining estimates of a smoother spatially-resolved with reduced variance when compared to a histogram. However, this ability that enables KDEs to obtain lower uncertainties per particle history is more computationally expensive than using a traditional histogram bins on a structured mesh. While KDEs are capable of producing lower tally uncertainties compared to histogram binning for a given number of histories, their increased computational cost may make them less efficient. Even so, it is possible to reduce this additional cost by exporting the KDE tally routines onto Graphics Processing Units (GPUs). GPUs are commonly found on high-performance computing machines as a means of increasing performance without adding additional compute nodes. While it is difficult for existing Monte Carlo radiation transport codes to take advantage of GPUs without rewriting large portions of the code, it is possible to leverage GPUs through heterogeneous computing. This work uses a modified version of OpenMC and exports the KDE tally routines onto a GPU to obtain problem-dependent speedups ranging between 1.6 and 5.0 for the problem studied in this paper. While this methodology is applied to KDEs, the same ideas can be applied to other computationally expensive vectorizable portions of Monte Carlo algorithms. (authors)
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Texture analysis of small renal cell carcinomas at MDCT for predicting relevant histologic and protein biomarkers
Purpose: To assess CT texture features of small renal cell carcinomas (≤ 4cm) for association with key pathologic features including protein biomarkers. Methods: Quantitative CT texture analysis (CTTA) of small renal cancers (≤ 4cm) was performed on non-contrast and portal venous phase abdominal MDCT scans with an ROI drawn at the largest cross-sectional diameter of the tumor using commercially available software. Texture parameters including mean pixel attenuation, the standard deviation (SD) of the pixel distribution histogram, entropy, the mean of positive pixels, the skewness (i.e., asymmetry) of the pixel histogram, kurtosis (i.e., peakness) of the pixel histogram, and the percentage of positive pixels were correlated with pathologic data from surgical resection, including histology and nuclear grade, as well as microarray analysis in a subset (n = 40) including Ki67 index, CRP, and neovascularization (CD105/CD31). Results: Portal venous phase images were available in 249 patients (105 women, 144 men; mean age, 56.7 years) with tumors ≤ 4cm (mean, median, range, ± SD; 2.66, 2.60, 0.3–4.0 ± 0.85 cm). CT texture features of standard deviation, mean of the positive pixels, and entropy of the pixel histogram were significantly associated with histologic cell type (clear vs. non-clear; p < 0.001). Entropy and mean of the positive pixels also showed an association with nuclear grade, although not statistically significant. In the microarray analysis subset, kurtosis of the pixel histogram was associated with CD105/CD31 (p = 0.05). SD also showed some association with CD 105 positivity (p = 0.02) and CAIX expression (p = 0.01). Non-contrast CT images were available in 174 patients (72 women, 102 men; mean age, 57.5 years). Although the association with histology was not as strong as on the portal venous phase, in the subset of patients with microarray data, SD was found to correlate with CRP (p = 0.08), kurtosis with CRP (p = 0.004), CD105/CD31 (p = 0.002), and with Ki 67 index (p < 0.001). Conclusion: CT texture features were significantly associated with important histopathologic features in small renal cancers. These non-invasive measures can be performed retrospectively and may provide useful information when determining follow-up and treatment of small renal cancers.
62 RADIOLOGY AND NUCLEAR MEDICINE↗
Cosmic-ray exposure history of ordinary chondrites
The exposure age histograms for H-, L-, and LL-chondrites are discussed. None of these histograms is consistent with a continuous delivery of asteroidal material to the earth, as the observed T(e) histograms clearly disagree with expected exponential distributions for a variety of orbital lifetimes. It is concluded that T(e) histograms are dominated by stochastic events and that the continuous supply of asteroidal material can account only for a minor background of the T(e) histograms. An attempt is made to identify major collisional levels among the major classes of ordinary chondrites in order to estimate the frequency of stochastic events. Orbital maturity in the inner solar system is documented by the p.m./total fall ratio among observed meteorite falls. All chondrite classes exhibit a uniform ratio of 2/3 except type H5 chondrites, which reveal a 0.5 or lower p.m./total fall ratio. This shift in the time of fall statistics suggests a strongly evolved orbit for the H5 parent at the time of collision about 7 Ma ago.
Marti, K.↗
Insights into Long-Term Global Precipitation from IMERG and GPCP
Multiple satellite-based datasets provide estimates of the long-term record of global precipitation. Each necessarily includes both the real atmospheric behavior and a collection of artifacts driven by the input data sources and design choices in the retrievals and dataset construction. For this presentation, the Integrated Multi-satellitE Retrievals for GPM (IMERG) products from the Global Precipitation Measurement (GPM) mission’s U.S. Science Team are designed as a modern High-Resolution Precipitation Product (HRPP), while the Global Precipitation Climatology Project’s (GPCP) products are designed as a modern Climate Data Record (CDR). Although they share some common input data sources, the computational process for each is rather different. Comparing the precipitation estimates from each algorithm gives us insights into likely artifacts and the natural variations that are common to both. The goal of an HRPP is to give the best estimate of precipitation at each time step in the dataset, generally meaning that “all possible available” data are used. In contrast, the CDR is intended to provide a precipitation record that has relatively homogeneous statistics, necessary for climate analysis. Each analysis strives for both goals, of course, but this means that the HRPP uses the disparate collection of satellites whose statistics don’t quite match the ideal record, while the CDR is computed from fewer satellites that provide a relatively homogeneous set of precipitation statistics. This presentation employs the current versions of the IMERG and GPCP products, V07B and V3.2, respectively. A variety of analyses have demonstrated that IMERG has more artifacts, as expected, which tend to result from changes in the source of intercalibration for the various input satellites, namely the Combined Radar-Radiometer Algorithm using TMI and PR during the TRMM era and using GMI and DPR-Ku during the GPM era. As well, a shift occurs when the altitudes of the TRMM and GPM satellites change as a result of orbit boosts. For the most part the mean precipitation shows good continuity across these boundaries, but there are noticeable changes in the respective histograms. The latter has taken on more importance in recent years due to increased scrutiny on extremes which implicitly focus on the upper end of the precipitation histograms. The GPCP product is more homogeneous, with only one major change in calibrator, transitioning from the SSMI series of satellites to the SSMIS series of satellites in 2009. The GPCP analysis does a good job of minimizing artifacts in the means at this boundary, while the histograms show a smaller, but still noticeable shift in the histograms, a result that is similar to data boundaries in IMERG. The GPCP Daily, which is a month-by-month rescaling of the daily accumulated IMERG Final product, shows more consistency than the equivalent daily IMERG, but inherits the histogram shifts from IMERG. Given this overlay of artifacts, there is enough consistency between the products to illustrate some important long-term variations, including interannual variations, not all of which are easily attributable to ENSO events, and trends that are large regionally, but comparatively small when averaged across the globe.
precipitation↗
Validation of the TOPAS Monte Carlo toolkit for HDR brachytherapy simulations
The goal of this work is to validate the user-friendly Geant4-based Monte Carlo toolkit TOol for PArticle Simulation (TOPAS) for brachytherapy applications. Brachytherapy simulations performed with TOPAS were systematically compared with published TG-186 reference data. The photon emission energy spectrum, the air-kerma strength, and the dose-rate constant of the model-based dose calculation algorithm (MBDCA)-WG generic Ir-192 source were extracted. For dose calculations, a track-length estimator was implemented. Additionally, the four Joint AAPM/ESTRO/ABG MBDCA-WG test cases were evaluated through histograms of the local and global dose difference volumes. A prostate, a palliative lung, and a breast case were simulated. For each case, the dose ratio map, the histogram of the global dose difference volume, and cumulative dose-volume histograms were calculated. The air-kerma strength was (9.772 ± 0.001) × 10 -8 U Bq -1 (within 0.3% of the reference value). The dose-rate constant was 1.1107 ± 0.0005 cGy h -1 U -1 (within 0.01% of the reference value). For all cases, at least 96.9% of voxels had a local dose difference within [-1%, 1%] and at least 99.9% of voxels had a global dose difference within [-0.1%, 0.1%]. The implemented track-length estimator scorer was more efficient than the default analog dose scorer by a factor of 237. For all clinical cases, at least 97.5% of voxels had a global dose difference within [-1%, 1%]. Dose-volume histograms were consistent with the reference data. TOPAS was validated for high-dose-rate brachytherapy simulations following the TG-186 recommended approach for MBDCAs. Built on top of Geant4, TOPAS provides broad access to a state-of-the-art Monte Carlo code for brachytherapy simulations.
59 BASIC BIOLOGICAL SCIENCES↗
Modeling the dosimetry of organ-at-risk in head and neck IMRT planning: An intertechnique and interinstitutional study
Purpose: To build a statistical model to quantitatively correlate the anatomic features of structures and the corresponding dose-volume histogram (DVH) of head and neck (HN) Tomotherapy (Tomo) plans. To study if the model built upon one intensity modulated radiation therapy (IMRT) technique (such as conventional Linac) can be used to predict anticipated organs-at-risk (OAR) DVH of patients treated with a different IMRT technique (such as Tomo). To study if the model built upon the clinical experience of one institution can be used to aid IMRT planning for another institution. Methods: Forty-four Tomotherapy intensity modulate radiotherapy plans of HN cases (Tomo-IMRT) from Institution A were included in the study. A different patient group of 53 HN fixed gantry IMRT (FG-IMRT) plans was selected from Institution B. The analyzed OARs included the parotid, larynx, spinal cord, brainstem, and submandibular gland. Two major groups of anatomical features were considered: the volumetric information and the spatial information. The volume information includes the volume of target, OAR, and overlapped volume between target and OAR. The spatial information of OARs relative to PTVs was represented by the distance-to-target histogram (DTH). Important anatomical and dosimetric features were extracted from DTH and DVH by principal component analysis. Two regression models, one for Tomotherapy plan and one for IMRT plan, were built independently. The accuracy of intratreatment-modality model prediction was validated by a leave one out cross-validation method. The intertechnique and interinstitution validations were performed by using the FG-IMRT model to predict the OAR dosimetry of Tomo-IMRT plans. The dosimetry of OARs, under the same and different institutional preferences, was analyzed to examine the correlation between the model prediction and planning protocol. Results: Significant patient anatomical factors contributing to OAR dose sparing in HN Tomotherapy plans have been analyzed and identified. For all the OARs, the discrepancies of dose indices between the model predicted values and the actual plan values were within 2.1%. Similar results were obtained from the modeling of FG-IMRT plans. The parotid gland was spared in a comparable fashion during the treatment planning of two institutions. The model based on FG-IMRT plans was found to predict the median dose of the parotid of Tomotherapy plans quite well, with a mean error of 2.6%. Predictions from the FG-IMRT model suggested the median dose of the larynx, median dose of the brainstem and D2 of the brainstem could be reduced by 10.5%, 12.8%, and 20.4%, respectively, in the Tomo-IMRT plans. This was found to be correlated to the institutional differences in OAR constraint settings. Re-planning of six Tomotherapy patients confirmed the potential of optimization improvement predicted by the FG-IMRT model was correct. Conclusions: The authors established a mathematical model to correlate the anatomical features and dosimetric indexes of OARs of HN patients in Tomotherapy plans. The model can be used for the setup of patient-specific OAR dose sparing goals and quality control of planning results. The institutional clinical experience was incorporated into the model which allows the model from one institution to generate a reference plan for another institution, or another IMRT technique.
62 RADIOLOGY AND NUCLEAR MEDICINE↗
TH-E-BRF-09: Gaussian Mixture Model Analysis of Radiation-Induced Parotid-Gland Injury: An Ultrasound Study of Acute and Late Xerostomia in Head-And-Neck Radiotherapy
Purpose: Xerostomia (dry mouth), secondary to parotid-gland injury, is a distressing side-effect in head-and-neck radiotherapy (RT). This study's purpose is to develop a novel ultrasound technique to quantitatively evaluate post-RT parotid-gland injury. Methods: Recent ultrasound studies have shown that healthy parotid glands exhibit homogeneous echotexture, whereas post-RT parotid glands are often heterogeneous, with multiple hypoechoic (inflammation) or hyperechoic (fibrosis) regions. We propose to use a Gaussian mixture model to analyze the ultrasonic echo-histogram of the parotid glands. An IRB-approved clinical study was conducted: (1) control-group: 13 healthy-volunteers, served as the control; (2) acutetoxicity group − 20 patients (mean age: 62.5 ± 8.9 years, follow-up: 2.0±0.8 months); and (3) late-toxicity group − 18 patients (mean age: 60.7 ± 7.3 years, follow-up: 20.1±10.4 months). All patients experienced RTOG grade 1 or 2 salivary-gland toxicity. Each participant underwent an ultrasound scan (10 MHz) of the bilateral parotid glands. An echo-intensity histogram was derived for each parotid and a Gaussian mixture model was used to fit the histogram using expectation maximization (EM) algorithm. The quality of the fitting was evaluated with the R-squared value. Results: (1) Controlgroup: all parotid glands fitted well with one Gaussian component, with a mean intensity of 79.8±4.9 (R-squared>0.96). (2) Acute-toxicity group: 37 of the 40 post-RT parotid glands fitted well with two Gaussian components, with a mean intensity of 42.9±7.4, 73.3±12.2 (R-squared>0.95). (3) Latetoxicity group: 32 of the 36 post-RT parotid fitted well with 3 Gaussian components, with mean intensities of 49.7±7.6, 77.2±8.7, and 118.6±11.8 (R-squared>0.98). Conclusion: RT-associated parotid-gland injury is common in head-and-neck RT, but challenging to assess. This work has demonstrated that the Gaussian mixture model of the echo-histogram could quantify acute and late toxicity of the parotid glands. This study provides meaningful preliminary data from future observational and interventional clinical research.
60 APPLIED LIFE SCIENCES↗
Postmortem validation of breast density using dual-energy mammography
Purpose: Mammographic density has been shown to be an indicator of breast cancer risk and also reduces the sensitivity of screening mammography. Currently, there is no accepted standard for measuring breast density. Dual energy mammography has been proposed as a technique for accurate measurement of breast density. The purpose of this study is to validate its accuracy in postmortem breasts and compare it with other existing techniques. Methods: Forty postmortem breasts were imaged using a dual energy mammography system. Glandular and adipose equivalent phantoms of uniform thickness were used to calibrate a dual energy basis decomposition algorithm. Dual energy decomposition was applied after scatter correction to calculate breast density. Breast density was also estimated using radiologist reader assessment, standard histogram thresholding and a fuzzy C-mean algorithm. Chemical analysis was used as the reference standard to assess the accuracy of different techniques to measure breast composition. Results: Breast density measurements using radiologist reader assessment, standard histogram thresholding, fuzzy C-mean algorithm, and dual energy were in good agreement with the measured fibroglandular volume fraction using chemical analysis. The standard error estimates using radiologist reader assessment, standard histogram thresholding, fuzzy C-mean, and dual energy were 9.9%, 8.6%, 7.2%, and 4.7%, respectively. Conclusions: The results indicate that dual energy mammography can be used to accurately measure breast density. The variability in breast density estimation using dual energy mammography was lower than reader assessment rankings, standard histogram thresholding, and fuzzy C-mean algorithm. Improved quantification of breast density is expected to further enhance its utility as a risk factor for breast cancer.
60 APPLIED LIFE SCIENCES↗
Evaluation of Matrix9 silicon photomultiplier array for small-animal PET
Purpose: The MatrixSL-9-30035-OEM (Matrix9) from SensL is a large-area silicon photomultiplier (SiPM) photodetector module consisting of a 3 × 3 array of 4 × 4 element SiPM arrays (total of 144 SiPM pixels) and incorporates SensL’s front-end electronics board and coincidence board. Each SiPM pixel measures 3.16 × 3.16 mm{sup 2} and the total size of the detector head is 47.8 × 46.3 mm{sup 2}. Using 8 × 8 polished LSO/LYSO arrays (pitch 1.5 mm) the performance of this detector system (SiPM array and readout electronics) was evaluated with a view for its eventual use in small-animal positron emission tomography (PET). Methods: Measurements of noise, signal, signal-to-noise ratio, energy resolution, flood histogram quality, timing resolution, and array trigger error were obtained at different bias voltages (28.0–32.5 V in 0.5 V intervals) and at different temperatures (5 °C–25 °C in 5 °C degree steps) to find the optimal operating conditions. Results: The best measured signal-to-noise ratio and flood histogram quality for 511 keV gamma photons were obtained at a bias voltage of 30.0 V and a temperature of 5 °C. The energy resolution and timing resolution under these conditions were 14.2% ± 0.1% and 4.2 ± 0.1 ns, respectively. The flood histograms show that all the crystals in the 1.5 mm pitch LSO array can be clearly identified and that smaller crystal pitches can also be resolved. Flood histogram quality was also calculated using different center of gravity based positioning algorithms. Improved and more robust results were achieved using the local 9 pixels for positioning along with an energy offset calibration. To evaluate the front-end detector readout, and multiplexing efficiency, an array trigger error metric is introduced and measured at different lower energy thresholds. Using a lower energy threshold greater than 150 keV effectively eliminates any mispositioning between SiPM arrays. Conclusions: In summary, the Matrix9 detector system can resolve high-resolution scintillator arrays common in small-animal PET with adequate energy resolution and timing resolution over a large detector area. The modular design of the Matrix9 detector allows it to be used as a building block for simple, low channel-count, yet high performance, small animal PET or PET/MRI systems.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗