Nonlinear least-squares modeling of 3D interaction position in a monolithic scintillator block

This paper presents a study of possible models to describe the relation between the scintillation light point-of-origin and the measured photo detector pixel signals in monolithic scintillation crystals. From these models the X, Y and depth of interaction (DOI) coordinates can be estimated simultaneously by nonlinear least-square fitting. The method depends only on the information embedded in the signals of individual events, and therefore does not need any prior position training or calibration. Three possible distributions of the light sources were evaluated: an exact solid-angle-based distribution, an approximate solid-angle distribution and an extended approximate solid-angle-based distribution which includes internal reflection at side and bottom surfaces. The performance of the general model using these three distributions was studied using Monte Carlo simulated data of a 20 x 20 x 10 mm lutetium oxyorthosilicate (Lu?SiO? or LSO) block read out by 2 Hamamatsu S8550 avalanche photo diode arrays. The approximate solid-angle-based model had the best compromise between resolution and simplicity. This model was also evaluated using experimental data by positioning a narrow 1.2 mm full width at half maximum (FWHM) beam of 511 keV photons at known positions on the 20 x 20 x 10 mm LSO block. An average intrinsic resolution in the X-direction of 1.4 mm FWHM was obtained for positions covering the complete block. The intrinsic DOI resolution was estimated at 2.6 mm FWHM.     

A method for determining the information capacity of x-ray imaging scintillator detectors by means of luminescence and modulation transfer function measurements

A method to determine the information capacity of x-ray phosphor screens used in the detectors of medical imaging systems is described. Information capacity was determined via x-ray luminescence efficiency (XLE), modulation transfer function (MTF) and emission spectrum measurements. The method was applied to laboratory prepared screens from commonly employed phosphor materials. The screen coating weight varied from 50 mg cm⁻² to 140 mg cm⁻². Results indicated that information capacity decreased with screen coating thickness but also depended on intrinsic phosphor properties (density, effective atomic number, intrinsic conversion efficiency, light wavelength). The Gd₂O₂S:Tb phosphor, exhibiting high density and effective atomic number, was found to be superior to La₂O₂S:Tb and Y₂O₂S:Tb.     

Parametric positioning of a continuous crystal PET detector with depth of interaction decoding

Here we demonstrate a parametric positioning method on a continuous crystal detector. Three different models for the light distribution were tested. Diagnosis of the residues showed that the parametric model fits the experimental data better than Gaussian and Cauchy models in our particular experimental setup. Based on the correlation between the spread and the peak value of the light distribution model with the depth of interaction (DOI), we were able to estimate the three-dimensional position of a scintillation event. On our continuous miniature crystal element (cMiCE) detector module with 8 mm thick LYSO crystal, the intrinsic spatial resolution is 1.06 mm at the center and 1.27 mm at the corner using a maximum-likelihood estimation (MLE) method and the parametric model. The DOI resolution (full width at half maximum) is estimated to be approximately 3.24 mm. The positioning method using the parametric model outperformed the Gaussian and Cauchy models, in both MLE and weighted least-squares (WLS) fitting methods. The key feature of this technique is that it requires very little calibration of the detector, but still retains high resolution and high sensitivity.     

Assessment of new small-field detectors against standard-field detectors for practical stereotactic beam data acquisition.

Two new detectors (0.015 cm3 ion chamber from PTW, 0.6 mm diameter diode from Scanditronix AB) designed specifically for use in small stereotactic fields were compared against similar, more routine, detectors (0.125 cm3 ion chamber, parallel plate chamber, shielded and unshielded diodes and film). Percentage depth doses, tissue maximum ratios, off-axis ratios and relative output factors were compared for circular fields in the 40-12.5 mm diameter range, with a view to identifying the optimum detector for stereotactic beam data acquisition. Practical suggestions for beam data collection and analysis are made, with an emphasis on what is achievable practically in radiotherapy departments where the primary demand is to provide a routine service. No single detector was found to be ideal, and neither of the two new measurement devices had any significant advantages over more routine devices, in the situations measured. Although the new 0.015 cm3 ion chamber was an improvement on a 0.125 cm3 ion chamber in the measurement of profiles, it was still too large when compared with a diode. The new small diode had a low signal to noise ratio which made reliable data difficult to extract and its only advantage is possibly improved resolution in fields smaller than the range tested. The use of a larger unshielded diode is recommended for all measurements, with the additional cross-checking of data against at least one small ion chamber and film. A simple method of obtaining reliable output data from the detectors used is explained.     

A new method for preventing pulse pileup in scintillation detectors

A new method for preventing pulse pileup in scintillation detectors is proposed. In the new method (G-INT), the energy of an event is calculated from the 'gated integral' of the pulse signal and the period of integration. The period of integration is not fixed but is shortened by the arrival of the succeeding pulse so as to avoid post-pulse pileup. The effect of pre-pulse pileup is corrected by subtracting the remnant energy of the preceding pulses, which is calculated from the gated integral of the preceding pulse. To avoid error due to short pulse intervals, pre- and post-pulse deadtimes are imposed. The method is similar to Wong's method (W-SUM) that depicts the energy by the 'weighted sum' of the current signal and the integrated signal. The performance of G-INT has been studied by Monte Carlo simulation in comparison with W-SUM, the variable sampling-time technique and simple delay-line clipping. It is shown that G-INT provides the smallest degradation in pulse height resolution for a given count rate capability. The difference between G-INT and W-SUM is explained by the difference in the amount of statistical noise involved in the gated integral and in the weighted sum.     

A count-rate model for PET scanners using pixelated Anger-logic detectors with different scintillators

A high count-rate simulation (HCRSim) model has been developed so that all results are derived from fundamental physics principles. Originally developed to study the behaviour of continuous sodium iodide (NaI(Tl)) detectors, this model is now applied to PET scanners based on pixelated Anger-logic detectors using lanthanum bromide (LaBr(3)), gadolinium orthosilicate (GSO) and lutetium orthosilicate (LSO) scintillators. This simulation has been used to study the effect on scanner deadtime and pulse pileup at high activity levels due to the scintillator stopping power (mu), decay time (tau) and energy resolution. Simulations were performed for a uniform 20 cm diameter x 70 cm long cylinder (NEMA NU2-2001 standard) in a whole-body scanner with an 85 cm ring diameter and a 25 cm axial field-of-view. Our results for these whole-body scanners demonstrate the potential of a pixelated Anger-logic detector and the relationship of its performance with the scanner NEC rate. Faster signal decay and short coincidence timing window lead to a reduction in deadtime and randoms fraction in the LaBr(3) and LSO scanners compared to GSO. The excellent energy resolution of LaBr(3) leads to the lowest scatter fraction for all scanners and helps compensate for reduced sensitivity compared to the GSO and LSO scanners, leading to the highest NEC values at high activity concentrations. The LSO scanner has the highest sensitivity of all the scanner designs investigated here, therefore leading to the highest peak NEC value but at a lower activity concentration than that of LaBr(3).     

Image blurring due to light-sharing in PET block detectors

The spatial resolution in PET is poorer than that of CT or MRI. All modern PET scanners use block detectors, i.e., clusters of scintillation crystals coupled to four photomultiplier tubes. Some of the loss of spatial resolution in PET is attributed to the use of block detectors, because a photon that interacts with one crystal in the cluster may be incorrectly positioned, resulting in blurring of the reconstructed image. This is called the "block effect." The block effect was measured for detectors from the CTI HR+ scanner, and the GE Advance scanner; two popular clinical PET scanners. The effect of changing the depth of first interaction of a gamma ray in the scintillation crystals was also studied to determine if it may be a contributor to the block effect. The block effect was found to be 1.2 +/- 0.5 mm for the central crystals and negligible for the edge crystals in the CTI HR+ block. It was 0.9 +/- 0.3 mm in the central crystals of the GE Advance detector, and 0.7 +/- 0.3 mm in the edge crystals of the GE Advance detector. In the CTI HR+ detector, a depth dependence on the positioning of the event was observed, as was a dependence on the crystal location (edge versus center). In the GE Advance detector events that occurred at different interaction depths were positioned consistently. The percentage of events that may be positioned inaccurately was also calculated for both detectors. In the CTI HR+ detector as many as 16% of all events in the block detector may be positioned incorrectly. In the GE Advance detector as many as 13% of all events in the block detector may be positioned inaccurately. These results suggest that the depth of interaction of an annihilation photon may contribute to the block effect in detectors that use crystals cut from a single scintillation crystal (pseudodiscrete crystals), and is less dominant a factor for detectors that use discrete crystals with light sharing between the crystals. Investigating the effect of changing photon interaction depth in PET detectors can lead to better detector design, and an intuitive explanation of what sources of blurring may exist in the detector examined.     

一种简化的PET绝对定量分析方法

目的 探讨在PET绝对定量分析中,如何在使用无损输入函数的同时,减PET动态显像时间.方法 通过勾画PET图像兔左心室感兴趣区(ROI)得到左心室血池时间活度曲线(PTAC),验证用血池感兴趣区法与标准输入函数模型结合的方法,确定无损输入函数可以替代动脉连续采血.利用输入函数模型拟合PTAC,输出函数模型拟合组织活度时间曲线(TTAC),再利用微分方程数值解法求解葡萄糖代谢房室模型.对该算法在缩短PET动态采集时间方面的应用进行了计算机仿真和动物实验上的验证.结果 (1)无创输入函数验证部分,两者通过非线性最小二乘法得到的ki相关系数r=0.9824(P<0.01).(2)计算机模拟结果表明在中低噪声情况下,可以将PET动态采集时间缩短为25min.动物实验结果表明,当PET动态采集时间缩短为25 min时,ki_(60)=0.9127×ki_(25)-0.0013,相关系数r=0.9182(P<0.01).结论 本方法可以通过选取合适的血池来获得无损输入函数,避免连续动脉采血,同时减少PET动态采集的时间.     

Optimal whole-body PET scanner configurations for different volumes of LSO scintillator: a simulation study

The axial field of view (AFOV) of the current generation of clinical whole-body PET scanners range from 15-22?cm, which limits sensitivity and renders applications such as whole-body dynamic imaging or imaging of very low activities in whole-body cellular tracking studies, almost impossible. Generally, extending the AFOV significantly increases the sensitivity and count-rate performance. However, extending the AFOV while maintaining detector thickness has significant cost implications. In addition, random coincidences, detector dead time, and object attenuation may reduce scanner performance as the AFOV increases. In this paper, we use Monte Carlo simulations to find the optimal scanner geometry (i.e. AFOV, detector thickness and acceptance angle) based on count-rate performance for a range of scintillator volumes ranging from 10 to 93 l with detector thickness varying from 5 to 20?mm. We compare the results to the performance of a scanner based on the current Siemens Biograph mCT geometry and electronics. Our simulation models were developed based on individual components of the Siemens Biograph mCT and were validated against experimental data using the NEMA NU-2 2007 count-rate protocol. In the study, noise-equivalent count rate (NECR) was computed as a function of maximum ring difference (i.e. acceptance angle) and activity concentration using a 27?cm diameter, 200?cm uniformly filled cylindrical phantom for each scanner configuration. To reduce the effect of random coincidences, we implemented a variable coincidence time window based on the length of the lines of response, which increased NECR performance up to 10% compared to using a static coincidence time window for scanners with a large maximum ring difference values. For a given scintillator volume, the optimal configuration results in modest count-rate performance gains of up to 16% compared to the shortest AFOV scanner with the thickest detectors. However, the longest AFOV of approximately 2?m with 20?mm thick detectors resulted in performance gains of 25-31?times higher NECR relative to the current Siemens Biograph mCT scanner configuration.     

A motion-incorporated reconstruction method for gated PET studies

Cardiac and respiratory motion artefacts in PET imaging have been traditionally resolved by acquiring the data in gated mode. However, gated PET images are usually characterized by high noise content due to their low photon statistics. In this paper, we present a novel 4D model for the PET imaging system, which can incorporate motion information to generate a motion-free image with all acquired data. A computer simulation and a phantom study were conducted to test the performance of this approach. The computer simulation was based on a digital phantom that was continuously scaled during data acquisition. The phantom study, on the other hand, used two spheres in a tank of water, all of which were filled with (18)F water. One of the spheres was stationary while the other moved in a sinusoidal fashion to simulate tumour motion in the thorax. Data were acquired using both 4D CT and gated PET. Motion information was derived from the 4D CT images and then used in the 4D PET model. Both studies showed that this 4D PET model had a good motion-compensating capability. In the phantom study, this approach reduced quantification error of the radioactivity concentration by 95% when compared to a corresponding static acquisition, while signal-to-noise ratio was improved by 210% when compared to a corresponding gated image.     

Two-dimensional phased arrays of sources and detectors for depth discrimination in diffuse optical imaging

We present a multisource, multidetector phased-array approach to diffuse optical imaging that is based on postprocessing continuous-wave data. We previously showed that this approach enhances the spatial resolution of diffuse optical imaging. We now demonstrate the depth discrimination capabilities of this approach and its potential to perform tomographic sectioning of turbid media. The depth discrimination results from the dependence of the sensitivity function on the depth coordinate z. To demonstrate the potential of this approach, we perform an experimental study of a turbid medium containing cylindrical inhomogeneities that are placed 2.0, 3.0, and 4.0 cm from a seven-element, 2-D source array. A single detector element is placed at a distance of 6.0 cm from the source array, and the measurement is repeated after switching the positions of the detector and the source array to simulate the case where both sources and detectors consist of a 2-D array of elements. We find that the proposed phased-array method is able to separate cylinders at different depths, thus showing cross-sectioning capabilities.     

A method for small-animal PET/CT alignment calibration

Small-animal positron-emission tomography/computed tomography (PET/CT) scanners provide anatomical and molecular imaging, which enables the joint visualization and analysis of both types of data. A proper alignment calibration procedure is essential for small-animal imaging since resolution is much higher than that in human devices. This work presents an alignment phantom and two different calibration methods that provide a reliable and repeatable measurement of the spatial geometrical alignment between the PET and the CT subsystems of a hybrid scanner. The phantom can be built using laboratory materials, and it is meant to estimate the rigid spatial transformation that aligns both modalities. It consists of three glass capillaries filled with a positron-emitter solution and positioned in a non-coplanar triangular geometry inside the system field of view. The calibration methods proposed are both based on automatic line detection, but with different approaches to calculate the transformation of the lines between both modalities. Our results show an average accuracy of the alignment estimation of 0.39 mm over the whole field of view.     

Analog signal multiplexing for PSAPD-based PET detectors: simulation and experimental validation

A 1 mm(3) resolution clinical positron emission tomography (PET) system employing 4608 position-sensitive avalanche photodiodes (PSAPDs) is under development. This paper describes a detector multiplexing technique that simplifies the readout electronics and reduces the density of the circuit board design. The multiplexing scheme was validated using a simulation framework that models the PSAPDs and front-end multiplexing circuits to predict the signal-to-noise ratio and flood histogram performance. Two independent experimental setups measured the energy resolution, time resolution, crystal identification ability and count rate both with and without multiplexing. With multiplexing, there was no significant degradation in energy resolution, time resolution and count rate. There was a relative 6.9 ± 1.0% and 9.4 ± 1.0% degradation in the figure of merit that characterizes the crystal identification ability observed in the measured and simulated ceramic-mounted PSAPD module flood histograms, respectively.     

A convolution-subtraction scatter correction method for 3D PET

3D acquisition and reconstruction in positron emission tomography (PET) produce data with improved signal-to-noise ratios compared with conventional 2D slice-oriented methods. However, the sensitivity increase is accompanied by an increase in the number of scattered photons and random coincidences detected. This paper presents a scatter correction technique for 3D PET data where an estimate of the scattered photon distribution is subtracted from the data before reconstruction. The scatter distribution is estimated by iteratively convolving the photopeak projections with a mono-exponential kernel. The method accounts for the 3D acquisition geometry and nature of scatter by performing the scatter estimation on 2D projections. The assumptions of the method have been investigated by measuring the variation in the scatter fraction and the scatter function at different positions in a cylinder. Both parameters were found to vary by up to 50% from the centre to the edge of a large water-filled cylinder. Despite this, in a uniform cylinder containing water with different concentrations of radioactivity, scatter was reduced from 25% in a non-radioactive region to less than 5% using the convolution-subtraction method. In addition, the relative concentration of a cylinder containing an increased concentration, which was underestimated by almost 50% without scatter correction, was within 5% of the true concentration after correction.     

Multilayer Gafchromic film detectors for breast skin dose determination in vivo

Assessment of skin dose delivered to patients from radiotherapy x-ray beams should be performed both inside and outside the prescribed treatment fields. A multilayer Gafchromic film detector which has high sensitivity for detection of radiation can be used to measure skin dose in a two-dimensional map over the skin surface if required. This is an advantage over other detectors, which only provide point dose estimates. A study of 25 patients undergoing breast irradiation was performed to analyse the ability of the multilayer detector to analyse skin dose and to assess both in-field and out-of-field radiation doses delivered during tangent field breast irradiation. Results show that the main contributor to total skin dose within the treatment field was delivered by exit dose. However, outside the field, most dose was delivered by entry beams. Patients with smaller breast separations where found, in general, to receive a higher total skin dose from entry and exiting beams at the central axis. Results also showed that a significant skin dose was delivered outside the treatment field and the main cause of this dose was from electron contamination from entry beams. The multilayer Gafchromic film detector provided adequate skin dose assessment within one fraction of treatment for in vivo results.