Iterative reconstruction of detector response of an Anger gamma camera

Experimental validation of the ANTS2 code for modelling optical photon transport in monolithic LYSO crystals

In this work the scintillation energy spectra originating from the background radioactivity from polished monolithic lutetium yttrium oxyorthosilicate coupled to position-sensitive silicon photomultipliers (SiPM) was studied using the open source Monte Carlo simulation package ANTS2. Two crystal sizes, fully and partially covering the photosensor area, three surface crystal wrappings (black, specular or diffuse) and the full signal formation process in the photosensor were considered. The simulation results were validated with experimental data acquired under the same geometric and detector operating conditions. In all cases ANTS2 simulated spectra have very good agreement with experimental results, reproducing the expected shape, with correct onset and end at 88 and 1190 keV, respectively, as well as sharp edges at the reference energies of 88, 88 + 202, 88 + 307 and 88 + 202 + 307 keV. The normalized root-mean square error between simulated and measured spectra varied between 4.3% and 10.4%.

Characterization of highly multiplexed monolithic PET / gamma camera detector modules

PET detectors use signal multiplexing to reduce the total number of electronics channels needed to cover a given area. Using measured thin-beam calibration data, we tested a principal component based multiplexing scheme for scintillation detectors. The highly-multiplexed detector signal is no longer amenable to standard calibration methodologies. In this study we report results of a prototype multiplexing circuit, and present a new method for calibrating the detector module with multiplexed data. A [Formula: see text] mm3 LYSO scintillation crystal was affixed to a position-sensitive photomultiplier tube with [Formula: see text] position-outputs and one channel that is the sum of the other 64. The 65-channel signal was multiplexed in a resistive circuit, with 65:5 or 65:7 multiplexing. A 0.9 mm beam of 511 keV photons was scanned across the face of the crystal in a 1.52 mm grid pattern in order to characterize the detector response. New methods are developed to reject scattered events and perform depth-estimation to characterize the detector response of the calibration data. Photon interaction position estimation of the testing data was performed using a Gaussian Maximum Likelihood estimator and the resolution and scatter-rejection capabilities of the detector were analyzed. We found that using a 7-channel multiplexing scheme (65:7 compression ratio) with 1.67 mm depth bins had the best performance with a beam-contour of 1.2 mm FWHM (from the 0.9 mm beam) near the center of the crystal and 1.9 mm FWHM near the edge of the crystal. The positioned events followed the expected Beer-Lambert depth distribution. The proposed calibration and positioning method exhibited a scattered photon rejection rate that was a 55% improvement over the summed signal energy-windowing method.

Iterative reconstruction of SiPM light response functions in a square-shaped compact gamma camera

Compact gamma cameras with a square-shaped monolithic scintillator crystal and an array of silicon photomultipliers (SiPMs) are actively being developed for applications in areas such as small animal imaging, cancer diagnostics and radiotracer guided surgery. Statistical methods of position reconstruction, which are potentially superior to the traditional centroid method, require accurate knowledge of the spatial response of each photomultiplier. Using both Monte Carlo simulations and experimental data obtained with a camera prototype, we show that the spatial response of all photomultipliers (light response functions) can be parameterized with axially symmetric functions obtained iteratively from flood field irradiation data. The study was performed with a camera prototype equipped with a 30  ×  30  ×  2 mm³ LYSO crystal and an 8  ×  8 array of SiPMs for 140 keV gamma rays. The simulations demonstrate that the images, reconstructed with the maximum likelihood method using the response obtained with the iterative approach, exhibit only minor distortions: the average difference between the reconstructed and the true positions in X and Y directions does not exceed 0.2 mm in the central area of 22  ×  22 mm² and 0.4 mm at the periphery of the camera. A similar level of image distortions is shown experimentally with the camera prototype.