Comparative study on gamma-ray detectors for in-situ ocean radiation monitoring system.

Large-sized crystals and state-of-the-art photosensors are desirable to cope with low environmental radioactivity (e.g., 1-2 Bq∙m-3137Cs in surface seawater) for homeland security purposes. We compared the performances of two different gamma-ray detector assemblies, GAGG crystal + silicon photomultiplier (SiPM) and NaI(Tl) crystal + photomultiplier tube, for our mobile in-situ ocean radiation monitoring system. We performed energy calibration, followed by water tank experiments with varying the depth of a137Cs point source. Experimental energy spectra were compared with MCNP-simulated spectra with identical setup and the consistency was validated. We finally assessed the detection efficiency and minimum detectable activity (MDA) of the detectors. Both GAGG and NaI detectors exhibited favorable energy resolutions (7.98 ± 0.13% and 7.01 ± 0.58% at 662 keV, respectively) and MDAs (33.1 ± 0.0645 and 13.5 ± 0.0327 Bq∙m-3 for 24-h 137Cs measurement, respectively). Matching the geometry of the GAGG crystal with that of the NaI crystal, the GAGG detector outperformed the NaI detector. The results demonstrated that the GAGG detector is potentially advantageous over the NaI detector in detection efficiency and compactness.

Underwater Optical-Sonar Image Fusion Systems.

Unmanned underwater operations using remotely operated vehicles or unmanned surface vehicles are increasing in recent times, and this guarantees human safety and work efficiency. Optical cameras and multi-beam sonars are generally used as imaging sensors in underwater environments. However, the obtained underwater images are difficult to understand intuitively, owing to noise and distortion. In this study, we developed an optical and sonar image fusion system that integrates the color and distance information from two different images. The enhanced optical and sonar images were fused using calibrated transformation matrices, and the underwater image quality measure (UIQM) and underwater color image quality evaluation (UCIQE) were used as metrics to evaluate the performance of the proposed system. Compared with the original underwater image, image fusion increased the mean UIQM and UCIQE by 94% and 27%, respectively. The contrast-to-noise ratio was increased six times after applying the median filter and gamma correction. The fused image in sonar image coordinates showed qualitatively good spatial agreement and the average IoU was 75% between the optical and sonar pixels in the fused images. The optical-sonar fusion system will help to visualize and understand well underwater situations with color and distance information for unmanned works.

Temporal Prediction of Paralytic Shellfish Toxins in the Mussel Mytilus galloprovincialis Using a LSTM Neural Network Model from Environmental Data.

Paralytic shellfish toxins (PSTs) are produced mainly by Alexandrium catenella (formerly A. tamarense). Since 2000, the National Institute of Fisheries Science (NIFS) has been providing information on PST outbreaks in Korean coastal waters at one- or two-week intervals. However, a daily forecast is essential for immediate responses to PST outbreaks. This study aimed to predict the outbreak timing of PSTs in the mussel Mytilus galloprovincialis in Jinhae Bay and along the Geoje coast in the southern coast of the Korea Peninsula. We used a long-short-term memory (LSTM) neural network model for temporal prediction of PST outbreaks from environmental data, such as water temperature (WT), tidal height, and salinity, measured at the Geojedo, Gadeokdo, and Masan tidal stations from 2006 to 2020. We found that PST outbreaks is gradually accelerated during the three years from 2018 to 2020. Because the in-situ environmental measurements had many missing data throughout the time span, we applied LSTM for gap-filling of the environmental measurements. We trained and tested the LSTM models with different combinations of environmental factors and the ground truth timing data of PST outbreaks for 5479 days as input and output. The LSTM model trained from only WT had the highest accuracy (0.9) and lowest false-alarm rate. The LSTM-based temporal prediction model may be useful as a monitoring system of PSP outbreaks in the coastal waters of southern Korean.

High Spatial-Resolution Red Tide Detection in the Southern Coast of Korea Using U-Net from PlanetScope Imagery.

Red tides caused by Margalefidinium polykrikoides occur continuously along the southern coast of Korea, where there are many aquaculture cages, and therefore, prompt monitoring of bloom water is required to prevent considerable damage. Satellite-based ocean-color sensors are widely used for detecting red tide blooms, but their low spatial resolution restricts coastal observations. Contrarily, terrestrial sensors with a high spatial resolution are good candidate sensors, despite the lack of spectral resolution and bands for red tide detection. In this study, we developed a U-Net deep learning model for detecting M. polykrikoides blooms along the southern coast of Korea from PlanetScope imagery with a high spatial resolution of 3 m. The U-Net model was trained with four different datasets that were constructed with randomly or non-randomly chosen patches consisting of different ratios of red tide and non-red tide pixels. The qualitative and quantitative assessments of the conventional red tide index (RTI) and four U-Net models suggest that the U-Net model, which was trained with a dataset of non-randomly chosen patches including non-red tide patches, outperformed RTI in terms of sensitivity, precision, and F-measure level, accounting for an increase of 19.84%, 44.84%, and 28.52%, respectively. The M. polykrikoides map derived from U-Net provides the most reasonable red tide patterns in all water areas. Combining high spatial resolution images and deep learning approaches represents a good solution for the monitoring of red tides over coastal regions.

FPGA-Based Interface of Digital DAQ System for Double-Scattering Compton Camera.

PURPOSE: The double-scattering Compton camera (DSCC) is a radiation imaging system that can provide both unknown source energy spectra and 3D spatial source distributions. The energies and detection locations measured in coincidence with three CdZnTe (CZT) detectors contribute to reconstructing emission energies and a spatial image based on conical surface integrals. In this study, we developed a digital data acquisition (DAQ) board to support our research into coincidence detection in the DSCC. METHODS: The main components of the digital DAQ board were 12 ADCs and one field programmable gate array (FPGA). The ADCs digitized the analog 96-channel CZT signals at a sampling rate of 50 MHz and transferred the serialized ADC samples and the bit and frame clocks to the FPGA. In order to correctly capture the ADC sample bits in the FPGA, we conducted individual sync calibrations for all the ADC channels to align the bit and frame clocks to the right positions of the ADC sample bits. The FPGA logic design was composed of IDELAY and IDDR components, six shift registers, and bit slip buffer resources. RESULTS: Using a Deskew test pattern, the delay value of the IDELAY component was determined to align the bit clock to the center of each sample bit. We determined the bit slip in the 12-bit ADC sample using an MSB test pattern by checking where the MSB value of one is located in the captured parallel data. CONCLUSIONS: After sync calibration, we tested the interface between the ADCs and the FPGA with a synthetic analog Gaussian signal. The 96 ADC channels yielded a mean R 2 goodness-of-fit value of 0.95 between the Gaussian curve and the captured 12-bit parallel data.

Analytic simulator and image generator of multiple-scattering Compton camera for prompt gamma ray imaging.

For prompt gamma ray imaging for biomedical applications and environmental radiation monitoring, we propose herein a multiple-scattering Compton camera (MSCC). MSCC consists of three or more semiconductor layers with good energy resolution, and has potential for simultaneous detection and differentiation of multiple radio-isotopes based on the measured energies, as well as three-dimensional (3D) imaging of the radio-isotope distribution. In this study, we developed an analytic simulator and a 3D image generator for a MSCC, including the physical models of the radiation source emission and detection processes that can be utilized for geometry and performance prediction prior to the construction of a real system. The analytic simulator for a MSCC records coincidence detections of successive interactions in multiple detector layers. In the successive interaction processes, the emission direction of the incident gamma ray, the scattering angle, and the changed traveling path after the Compton scattering interaction in each detector, were determined by a conical surface uniform random number generator (RNG), and by a Klein-Nishina RNG. The 3D image generator has two functions: the recovery of the initial source energy spectrum and the 3D spatial distribution of the source. We evaluated the analytic simulator and image generator with two different energetic point radiation sources (Cs-137 and Co-60) and with an MSCC comprising three detector layers. The recovered initial energies of the incident radiations were well differentiated from the generated MSCC events. Correspondingly, we could obtain a multi-tracer image that combined the two differentiated images. The developed analytic simulator in this study emulated the randomness of the detection process of a multiple-scattering Compton camera, including the inherent degradation factors of the detectors, such as the limited spatial and energy resolutions. The Doppler-broadening effect owing to the momentum distribution of electrons in Compton scattering was not considered in the detection process because most interested isotopes for biomedical and environmental applications have high energies that are less sensitive to Doppler broadening. The analytic simulator and image generator for MSCC can be utilized to determine the optimal geometrical parameters, such as the distances between detectors and detector size, thus affecting the imaging performance of the Compton camera prior to the development of a real system.

Direct Reconstruction of CT-based Attenuation Correction Images for PET with Cluster-Based Penalties.

Extremely low-dose CT acquisitions used for PET attenuation correction have high levels of noise and potential bias artifacts due to photon starvation. This work explores the use of a priori knowledge for iterative image reconstruction of the CT-based attenuation map. We investigate a maximum a posteriori framework with cluster-based multinomial penalty for direct iterative coordinate decent (dICD) reconstruction of the PET attenuation map. The objective function for direct iterative attenuation map reconstruction used a Poisson log-likelihood data fit term and evaluated two image penalty terms of spatial and mixture distributions. The spatial regularization is based on a quadratic penalty. For the mixture penalty, we assumed that the attenuation map may consist of four material clusters: air+background, lung, soft tissue, and bone. Using simulated noisy sinogram data, dICD reconstruction was performed with different strengths of the spatial and mixture penalties. The combined spatial and mixture penalties reduced the RMSE by roughly 2 times compared to a weighted least square and filtered backprojection reconstruction of CT images. The combined spatial and mixture penalties resulted in only slightly lower RMSE compared to a spatial quadratic penalty alone. For direct PET attenuation map reconstruction from ultra-low dose CT acquisitions, the combination of spatial and mixture penalties offers regularization of both variance and bias and is a potential method to reconstruct attenuation maps with negligible patient dose. The presented results, using a best-case histogram suggest that the mixture penalty does not offer a substantive benefit over conventional quadratic regularization and diminishes enthusiasm for exploring future application of the mixture penalty.

Mixed Confidence Estimation for Iterative CT Reconstruction.

Dynamic (4D) CT imaging is used in a variety of applications, but the two major drawbacks of the technique are its increased radiation dose and longer reconstruction time. Here we present a statistical analysis of our previously proposed Mixed Confidence Estimation (MCE) method that addresses both these issues. This method, where framed iterative reconstruction is only performed on the dynamic regions of each frame while static regions are fixed across frames to a composite image, was proposed to reduce computation time. In this work, we generalize the previous method to describe any application where a portion of the image is known with higher confidence (static, composite, lower-frequency content, etc.) and a portion of the image is known with lower confidence (dynamic, targeted, etc). We show that by splitting the image space into higher and lower confidence components, MCE can lower the estimator variance in both regions compared to conventional reconstruction. We present a theoretical argument for this reduction in estimator variance and verify this argument with proof-of-principle simulations. We also propose a fast approximation of the variance of images reconstructed with MCE and confirm that this approximation is accurate compared to analytic calculations of and multi-realization image variance. This MCE method requires less computation time and provides reduced image variance for imaging scenarios where portions of the image are known with more certainty than others allowing for potentially reduced radiation dose and/or improved dynamic imaging.

Mixed Confidence Estimation for Iterative CT Reconstruction.

We present a statistical analysis of our previously proposed Constrain-Static Target-Kinetic algorithm for 4D CT reconstruction. This method, where framed iterative reconstruction is only performed on the dynamic regions of each frame, while static regions are fixed across frames to a composite image, was proposed to reduce computation time. In this work, we generalize the previous method to describe any application where a portion of the image is known with higher confidence (static, composite, lower-frequency content, etc.) and a portion of the image is known with lower confidence (dynamic, targeted, etc). We show that by splitting the image space into higher and lower confidence components, CSTK can lower the estimator variance in both regions compared to conventional reconstruction. We present a theoretical argument for this reduction in estimator variance and verify this argument with proof-of-principle simulations. This method allows for reduced computation time and improved image quality for imaging scenarios where portions of the image are known with more certainty than others.

Resolution recovery reconstruction for a Compton camera.

The spatial resolution from Compton cameras suffers from measurement uncertainties in interaction positions and energies. The degree of degradation in spatial resolution is shift-variant (SV) over the field-of-view (FOV) because the imaging principle is based on the conical surface integration. In our study, the shift-variant point spread function (SV-PSF) is derived from point source measurements at various positions in the FOV and is incorporated into the system matrix of a fully three-dimensional, accelerated reconstruction, i.e. the listmode ordered subset expectation maximization (LMOSEM) algorithm, for resolution recovery. Simulation data from point sources were used to estimate SV and asymmetric parameters for Gaussian, Cauchy, and general parametric PSFs. Although little difference in the fitness accuracy between Gaussian and general parametric PSFs was observed, the general parametric model showed greater flexibility over the FOV in shaping the curve between that for Gaussian and Cauchy functions. The estimated asymmetric SV-PSFs were incorporated into the LMOSEM for resolution recovery. For simulation data from a single point source at the origin, all LMOSEM-SV-PSFs improved the spatial resolution by 2.6 times over the standard LMOSEM. For two point-source simulations, reconstructions also gave a two-fold improvement in spatial resolution and resulted in a greater recovered activity ratio at different positions in the FOV.

Direct Reconstruction of CT-based Attenuation Correction Images for PET with Cluster-Based Penalties.

Extremely low-dose CT acquisitions for the purpose of PET attenuation correction will have a high level of noise and biasing artifacts due to factors such as photon starvation. This work explores a priori knowledge appropriate for CT iterative image reconstruction for PET attenuation correction. We investigate the maximum a posteriori (MAP) framework with cluster-based, multinomial priors for the direct reconstruction of the PET attenuation map. The objective function for direct iterative attenuation map reconstruction was modeled as a Poisson log-likelihood with prior terms consisting of quadratic (Q) and mixture (M) distributions. The attenuation map is assumed to have values in 4 clusters: air+background, lung, soft tissue, and bone. Under this assumption, the MP was a mixture probability density function consisting of one exponential and three Gaussian distributions. The relative proportion of each cluster was jointly estimated during each voxel update of direct iterative coordinate decent (dICD) method. Noise-free data were generated from NCAT phantom and Poisson noise was added. Reconstruction with FBP (ramp filter) was performed on the noise-free (ground truth) and noisy data. For the noisy data, dICD reconstruction was performed with the combination of different prior strength parameters (ß and γ) of Q- and M-penalties. The combined quadratic and mixture penalties reduces the RMSE by 18.7% compared to post-smoothed iterative reconstruction and only 0.7% compared to quadratic alone. For direct PET attenuation map reconstruction from ultra-low dose CT acquisitions, the combination of quadratic and mixture priors offers regularization of both variance and bias and is a potential method to derive attenuation maps with negligible patient dose. However, the small improvement in quantitative accuracy relative to the substantial increase in algorithm complexity does not currently justify the use of mixture-based PET attenuation priors for reconstruction of CT images for PET attenuation correction.

Gap compensation during PET image reconstruction by constrained, total variation minimization.

PURPOSE: Positron emission tomography (PET) is a noninvasive molecular imaging tool with various clinical and preclinical applications. The polygonal structure of small-diameter PET scanners that are designed for specific purposes can lead to gaps between the detector modules and result in loss of PET data during measurement. In the current study, the authors applied the compressed sensing (CS)-based total variation (TV) minimization method to PET image reconstructions to reduce the artifacts caused by gaps in small-diameter PET systems. METHODS: The first step in each iteration estimates whether an image is consistent with the measured PET data using the existing common reconstruction algorithms (ART, OSEM, and RAMLA). The second step recovers sparsity in the gradient domain of the image by minimizing the TV of an estimated image. The authors evaluated the gap-compensable reconstruction algorithms with uniform disk and Shepp-Logan phantoms by simulating sinograms which contained Poisson random noise and a data loss due to detector gaps. In addition, these methods were applied to real high resolution research tomography (HRRT)-like sinograms of human brain and uniform phantom. A comparison with other methods for gap compensation prior to or during image reconstruction was also made. Quantitative evaluations were performed by computing the uniformity, root mean squared error, and difference between the reconstructed images of nongapped and gapped sinograms. RESULTS: The simulation results showed that the gap-compensable methods incorporating TV minimization could control gap artifacts, as well as Poisson random noise. In particular, OSEM-TV and RAMLA-TV showed distinct potential via the properties of convergence and robustness to different noise levels and gap angle. CONCLUSIONS: A TV minimization strategy incorporated into commonly used PET reconstruction algorithms was useful for reducing the occurrence of artifacts due to gaps between detector modules in small-diameter PET scanners.

Fully three-dimensional OSEM-based image reconstruction for Compton imaging using optimized ordering schemes.

Although the ordered subset expectation maximization (OSEM) algorithm does not converge to a true maximum likelihood solution, it is known to provide a good solution if the projections that constitute each subset are reasonably balanced. The Compton scattered data can be allocated to subsets using scattering angles (SA) or detected positions (DP) or a combination of the two (AP (angles and positions)). To construct balanced subsets, the data were first arranged using three ordering schemes: the random ordering scheme (ROS), the multilevel ordering scheme (MLS) and the weighted-distance ordering scheme (WDS). The arranged data were then split into J subsets. To compare the three ordering schemes, we calculated the coefficients of variation (CVs) of angular and positional differences between the arranged data and the percentage errors between mathematical phantoms and reconstructed images. All ordering schemes showed an order-of-magnitude acceleration over the standard EM, and their computation times were similar. The SA-based MLS and the DP-based WDS led to the best-balanced subsets (they provided the largest angular and positional differences for SA- and DP-based arrangements, respectively). The WDS exhibited minimum CVs for both the SA- and DP-based arrangements (the deviation in mean angular and positional differences between the ordered subsets was smallest). The combination of AP and WDS yielded the best results with the lowest percentage errors by providing larger and more uniform angular and positional differences for the SA and DP arrangements, and thus, is probably optimal Compton camera reconstruction using OSEM.