BackgroundThe aim of the study is to assess accuracy of activity quantification of 177Lu studies performed according to recommendations provided by the committee on Medical Internal Radiation Dose (MIRD) pamphlets 23 and 26. The performances of two scatter correction and three segmentation methods were compared. Additionally, the accuracy of tomographic and planar methods for determination of the camera normalization factor (CNF) was evaluated.Eight phantoms containing inserts of different sizes and shapes placed in air, water, and radioactive background were scanned using a Siemens SymbiaT SPECT/CT camera. Planar and tomographic scans with 177Lu sources were used to measure CNF. Images were reconstructed with our SPEQToR software using resolution recovery, attenuation, and two scatter correction methods (analytical photon distribution interpolated (APDI) and triple energy window (TEW)). Segmentation was performed using a fixed threshold method for both air and cold water scans. For hot water experiments three segmentation methods were compared as folows: a 40% fixed threshold, segmentation based on CT images, and our iterative adaptive dual thresholding (IADT). Quantification error, defined as the percent difference between experimental and true activities, was evaluated.ResultsQuantification error for scans in air was better for TEW scatter correction (<6%) than for APDI (<11%). This trend was reversed for scans in water (<10% for APDI and <14% for TEW). For hot water, the best results (<18% for small objects and <5% for objects >100 ml) were obtained when APDI and IADT were used for scatter correction and segmentation, respectively. Additionally, we showed that planar acquisitions with scatter correction and tomographic scans provide similar CNF values. This is an important finding because planar acquisitions are easier to perform than tomographic scans. TEW and APDI resulted in similar quantification errors with APDI showing a small advantage for objects placed in medium with non-uniform density.ConclusionsFollowing the MIRD recommendations for data acquisition and reconstruction resulted in accurate activity quantification (errors <5% for large objects). However, techniques for better organ/tumor segmentation must still be developed.Electronic supplementary materialThe online version of this article (doi:10.1186/s40658-016-0170-3) contains supplementary material, which is available to authorized users.

BackgroundCamera calibration, which translates reconstructed count map into absolute activity map, is a prerequisite procedure for quantitative SPECT imaging. Both planar and tomographic scans using different phantom geometries have been proposed for the determination of the camera calibration factor (CF). However, there is no consensus on which approach is the best. The aim of this study is to evaluate all these calibration methods, compare their performance, and propose a practical and accurate calibration method for SPECT quantitation of therapeutic radioisotopes. Twenty-one phantom experiments (Siemens Symbia SPECT/CT) and 12 Monte Carlo simulations (GATE v6.1) using three therapy isotopes (131I, 177Lu, and 188Re) have been performed. The following phantom geometries were used: (1) planar scans of point source in air (PS), (2) tomographic scans of insert(s) filled with activity placed in non-radioactive water (HS + CB), (3) tomographic scans of hot insert(s) in radioactive water (HS + WB), and (4) tomographic scans of cylinders uniformly filled with activity (HC). Tomographic data were reconstructed using OSEM with CT-based attenuation correction and triple energy window (TEW) scatter correction, and CF was determined using total counts in the reconstructed image, while for planar scans, the photopeak counts, corrected for scatter and background with TEW, were used. Additionally, for simulated data, CF obtained from primary photons only was analyzed.ResultsFor phantom experiments, CF obtained from PS and HS + WB agreed to within 6% (below 3% if experiments performed on the same day are considered). However, CF from HS + CB exceeded those from PS by 4–12%. Similar trend was found in simulation studies. Analysis of CFs from primary photons helped us to understand this discrepancy. It was due to underestimation of scatter by the TEW method, further enhanced by attenuation correction. This effect becomes less important when the source is distributed over the entire phantom volume (HS + WB and HC).ConclusionsCamera CF could be determined using planar scans of a point source, provided that the scatter and background contributions are removed, for example using the clinically available TEW method. This approach is simple and yet provides CF with sufficient accuracy (~ 5%) to be used in clinics for radiotracer quantification.

177Lu-DOTATATE therapy has been shown to produce encouraging results in treatment of neuroendocrine tumours (NETs). Unfortunately, since dosimetry for radionuclide therapy is considered to be challenging, typically similar amount of radiopharmaceutical is administered to every patient. There is growing evidence that the efficacy of this therapy can be significantly improved by employing personalized protocols, based on the organ-at-risk dosimetry. The aim of this study is to propose a practical and accurate dosimetry protocol based on the simplified acquisition schedules. Data from fifty-three therapy cycles in thirty-nine NET patients were analyzed. Three SPECT/CT scans, acquired at 4 h (D0), 23 h (D1) and 70 h (D3) after injection, were performed. The kidney volume was determined using CT and the activity was determined from quantitative SPECT using an iterative thresholding method. For each dataset, four methods were used to model the time-activity-curves (TAC): M1—two trapezoid segments (0 to D0 and D0 to D1), followed by monoexponential fit to D1  +  D3 data; M2—monoexponential fit to D0  +  D1  +  D3 data; M3 and M4—monoexponential fit to D0  +  D3 and D1  +  D3 data, respectively. Additionally, kidney doses obtained from single time point method using a monoexponential curve with the population mean effective half-life, normalized to activities at D0 or D1 or D3 points, were calculated. The accuracy of simplified dosimetry methods was assessed as the percentage difference relative to doses calculated from M1. The major contribution to the absorbed dose estimate comes from the area under the late time portion of the TAC (D1 to infinity). Therefore, information from the late scan (D3) is crucial for the determination of kidney absorbed doses. Single time point method using monoexponential TAC, with the population mean effective half-life normalized to the late data point (48–72 h for kidneys) produces  <10% deviation in the absorbed dose estimation, thus is recommended for clinical use.

BackgroundThe aim of this study was to investigate the deadtime (DT) effects that are present in 177Lu images acquired after radionuclide therapy injection, assess differences in DT based on the full spectrum and the photopeak-only measurements, and design a method to correct for the deadtime losses.A Siemens SymbiaT SPECT/CT camera with a medium energy collimator was used. A 295-mL bottle was placed off-center inside a large cylinder filled with water, and 177Lu activity was sequentially added up to a maximum of 9.12 GBq. The true count rates vs. observed count rates were plotted and fitted to the DT paralyzable model. This analysis was performed using counts recorded in the full spectrum and in other energy windows. The DT correction factors were calculated using the percentage difference between the true and the observed count rates.ResultsThe DT values of 5.99 ± 0.02 μs, 4.60 ± 0.052 μs, and 0.19 ± 0.18 μs were obtained for the primary photons (PP) recorded in the 113- and 208-keV photopeaks and for the full spectrum, respectively. For the investigated range of count rates, the DT correction factors of up to 23% were observed for PP corresponding to the 113-keV photopeak, while for the 208-keV photopeak values of up to 20% were obtained. These values were almost three times higher than the deadtime correction factors derived from the full spectrum.ConclusionsThe paralyzable model showed to be appropriate for the investigated range of counts, which were five to six times higher than those observed in the patient post-therapy imaging. Our results suggest that the deadtime corrections should be based on count losses in the scatter-corrected photopeak window and not on the deadtime determined from the full spectrum. Finally, a general procedure that can be followed to correct patient images for deadtime is presented.

PET images acquired after liver 90Y radioembolization therapies are typically very noisy, which significantly challenges both visualization and quantification of activity distributions. To improve their noise characteristics, regularized iterative reconstruction algorithms such as block sequential regularized expectation maximization (Q.Clear for GE Healthcare, USA) have been proposed. In this study, we aimed to investigate the effects which different reconstruction algorithms may have on patient images, with reconstruction parameters initially narrowed down using phantom studies. Moreover, we evaluated the impact of these reconstruction methods on voxel-based dose distribution in phantom and patient studies (lesions and healthy livers). The International Electrotechnical Commission (IEC)/NEMA phantom, containing six spheres, was filled with 90Y and imaged using a GE Discovery 690 PET/CT scanner with time-of-flight enabled. The images were reconstructed using Q.Clear (with β parameter ranging from 0 to 8000) and ordered subsets expectation maximization. The image quality and quantification accuracy were evaluated by computing the hot () and cold () contrast recovery coefficients, background variability (BV) and activity bias. Next, dose distributions and dose volume histograms were generated using MIM® software’s SurePlan LiverY90 toolbox. Subsequently, parameters optimized in these phantom studies were applied to five patient datasets. Dose parameters, such as max, mean, 70, and 100 , were estimated, and their variability for different reconstruction methods was investigated. Based on phantom studies, the parameter values optimized for image quality and quantification accuracy were 2500 and 300, respectively. When all investigated reconstructions were applied to patient studies, mean, D50, 70, and 100 showed coefficients of variation below 8%; whereas the variability of max was up to 30% for both phantom and patient images. Although = 300−1000 would provide accurate activity quantification for a region of interest, when considering activity/dose voxelized distribution, higher value (e.g. 4000–5000) would provide the greatest accuracy for dose distributions. In this 90Y radioembolization PET/CT study, the parameter in regularized iterative (Q.Clear) reconstruction was investigated for image quality, accurate quantification and dose distributions based on phantom experiments and then applied to patient studies. Our results indicate that more accurate dose distribution can be achieved from smoother PET images, reconstructed with larger values than those yielding the best activity quantifications but noisy images. Most importantly, these results suggest that quantitative measures, which are commonly used in clinics, such as SUVmax or SUVpeak( equivalent of max), should not be employed for 90Y PET images, since their values would hig...