Threshold segmentation for PET target volume delineation in radiation treatment planning: the role of target-to-background ratio and target size.

A multivariable approach was adopted to study the dependence of the percentage threshold [TH(%)] used to define the boundaries of 18F-FDG positive tissue on emission scan duration (ESD) and activity at the start of acquisition (Aacq) for different target sizes and target-to-background (T/B) ratios. An anthropomorphic model, at least for counting rate characteristics, was used to study this dependence in conditions resembling the ones that can be encountered in the clinical studies. An annular ring of water bags of 3 cm thickness was fitted over an International Electro-technical Commission (IEC) phantom in order to obtain counting rates similar to those found in average patients. The scatter fraction of the modified IEC phantom was similar to the mean scatter fraction measured on patients, with a similar scanner. A supplemental set of microhollow spheres was positioned inside the phantom. The NEMA NU 2-2001 scatter phantom was positioned at the end of the IEC phantom to approximate the clinical situation of having activity that extends beyond the scanner field of view. The phantoms were filled with a solution of water and 18F (12 kBq/mL) and the spheres with various T/B ratios of 22.5, 10.3, and 3.6. Sequential imaging was performed to acquire PET images with varying background activity concentrations of about 12, 9, 6.4, 5.3, and 3.1 kBq/mL. The ESD was set to 60, 120, 180, and 240 s/bed. Data were fitted using two distinct multiple linear regression models for sphere ID < or = 10 mm and sphere ID > 10 mm. The fittings of both models were good with an R2 of 0.86 in both cases. Neither ESD nor Aacq resulted as significant predictors of the TH(%). For sphere ID < or =10 mm the target size was the most significant predictor of the TH(%), followed by the T/B ratio, while for sphere ID > 10 mm the explanatory power of the target size and T/B ratio were reversed, the T/B ratio being now the most important predictor of the TH(%). Both the target size and T/B ratio play a major role in explaining the variance of the TH(%), throughout the whole range of target sizes and T/B ratios examined. Thus, algorithms aimed at automatic threshold segmentation should incorporate both variables with a relative weight which critically depends on target size.

The use of molecular sieves to simulate hot lesions in (18)F-fluorodeoxyglucose--positron emission tomography imaging.

We investigated the use of a kind of zeolite, the Bowie chabazite, to produce radioactive sources of different shapes, dimensions and activity concentrations that can be used for lesion simulation in positron emission tomography (PET) imaging. The (18)F-fluorodeoxyglucose ((18)F-FDG) uptake of a group of 12 zeolites was studied as a function of their weight (120-1,520 mg) and of the activity concentration of the (18)F-FDG solution (1-37 MBq ml(-1)), using a multiple linear regression model. The reproducibility, homogeneity and stability over time of the (18)F-FDG uptake were assessed. The fit of the regression model is good (r(2) = 0.83). This relation allows the production of zeolites of a desired (18)F-FDG activity using knowledge of the concentration of the soaking solution and the weight of the zeolite. The reproducibility of the (18)F-FDG uptake after heating the zeolites is elevated (CV% = 3.68). The almost complete regeneration of the zeolites allows us to reuse them in successive experiments. The stability of the (18)F-FDG uptake on zeolites is far from ideal. When placed in a saline solution the 'activated' zeolites release the (18)F-FDG with an effective half-time of 53 min. The sealing of the zeolites in plastic film bags has been demonstrated to be effective in preventing any release of (18)F-FDG. These features, together with their variable dimensions and shapes, make them ideal (18)F-FDG sources with a fixed target-to-background ratio that can be placed anywhere in a phantom to study lesion detectability in PET imaging.

Impact of target-to-background ratio, target size, emission scan duration, and activity on physical figures of merit for a 3D LSO-based whole body PET/CT scanner.

The aim of our work is to describe the way in which physical figures of merit such as contrast-to-noise ratio (CNR) behave when varying acquisition parameters such as emission scan duration (ESD) or activity at the start of acquisition (A(acq)) that in clinical practice can be selected by the user, or object properties such as target dimensions or target-to-background (T/B) ratio, which depend uniquely on the intrinsic characteristics of the object being imaged. Figures of merit, used to characterize image quality and quantitative accuracy for a 3D-LSO based PET/CT scanner, were studied as a function of ESD and A(acq) for different target sizes and T/B ratios using a multivariate approach in a wide range of conditions approaching the ones that can be encountered in clinical practice. An annular ring of water bags of 3 cm thickness was fitted over an IEC phantom in order to obtain counting rates similar to those found in average patients. The average scatter fraction (SF) of the modified IEC phantom was similar to the mean SF measured on patients with a similar scanner. A supplemental set of micro-hollow spheres was positioned inside the phantom. The NEMA NU 2-2001 scatter phantom was positioned at the end of the IEC phantom to approximate the clinical situation of having activity that extends beyond the scanner. The phantoms were filled with a solution of water and 18F (12 kBq/mL) and the spheres with various T/B ratios of 22.5, 10.3, and 3.6. Sequential imaging was performed to acquire PET images with varying background activity concentrations of about 12, 9, 6.4, 5.3, and 3.1 kBq/mL, positioned on the linear portion of the phantom's NECR curve, well below peak NECR of 61.2 kcps that is reached at 31.8 kBq/mL. The ESD was set to 1, 2, 3, and 4 min/bed. With T/B ratios of 3.6, 10.3, and 22.5, the 13.0, 8.1, and 6.5 mm spheres were detectable for the whole ranges of background activity concentration and ESD, respectively. The ESD resulted as the most significant predictor of CNR variance, followed by T/B ratio and the cross sectional area of the given sphere. Only last comes A(acq) with a weight more than halved with respect to ESD. Thus, raising ESD seems to be much more effective than raising A(acq) in order to obtain higher CNR, which is the physical figure of merit closely related with target detectability, at least in the simple task of the signal known exactly background known exactly model.

Conversion factors of effective and equivalent organ doses with the air kerma area product in patients undergoing coronary angiography and percutaneous coronary interventions.

To derive effective dose (E), organ dose (HT) and conversion factors with the air kerma area product (KAP) in coronary angiography (CA) and percutaneous coronary intervention (PCI) by the radial route, using the ICRP 103 tissue weighting factors. The study included 34 patients referred for CA and 31 for PCI. E and HT were derived from in-the-field KAP measurements using Montecarlo methods. Median KAP of 23.2 and 56.8Gycm2 and E of 6.9 and 20.0mSv were found for CA and PCI, respectively. Mean KAP and E were significantly higher in males than in females (52.4±40.0 vs 32.3±16.6Gycm2; p=0.02) and (16.8±13.6 vs 10.7±5.8mSv; p=0.04). KAP (r=0.39; p=0.001) and E (r=0.34; p=0.005) showed a significant correlation with the patient's weight. Conversion factors between KAP and E (E/KAP) were 0.30±0.04mSvGy-1cm-2 for CA and 0.33±0.05mSvGy-1cm-2 for PCI. No significant differences in the E/KAP between males and females were found (0.31±0.05 vs 0.33±0.05; p=0.08). Again, no significant correlation was found between E/KAP and patient's weight (r=0.23; p=0.07). The correlation between E and KAP was excellent for CA (r=0.99) and PCI (r=0.96). The correlation between HT and KAP ranged from r=0.87 to r=1 and from r=0.71 to r=0.98 for CA and PCI, respectively. A single factor, the total KAP, could be used for a specific acquisition protocol to reliably estimate E and HT without the need of a patient's specific analysis. Conversion factors might be installation, X-ray beam quality or protocol dependent.

Performance comparison of two resolution modeling PET reconstruction algorithms in terms of physical figures of merit used in quantitative imaging.

PURPOSE: Resolution modeling (RM) of PET systems has been introduced in iterative reconstruction algorithms for oncologic PET. The RM recovers the loss of resolution and reduces the associated partial volume effect. While these methods improved the observer performance, particularly in the detection of small and faint lesions, their impact on quantification accuracy still requires thorough investigation. The aim of this study was to characterize the performances of the RM algorithms under controlled conditions simulating a typical (18)F-FDG oncologic study, using an anthropomorphic phantom and selected physical figures of merit, used for image quantification. METHODS: Measurements were performed on Biograph HiREZ (B_HiREZ) and Discovery 710 (D_710) PET/CT scanners and reconstructions were performed using the standard iterative reconstructions and the RM algorithms associated to each scanner: TrueX and SharpIR, respectively. RESULTS: RM determined a significant improvement in contrast recovery for small targets (≤17 mm diameter) only for the D_710 scanner. The maximum standardized uptake value (SUVmax) increased when RM was applied using both scanners. The SUVmax of small targets was on average lower with the B_HiREZ than with the D_710. Sharp IR improved the accuracy of SUVmax determination, whilst TrueX showed an overestimation of SUVmax for sphere dimensions greater than 22 mm. The goodness of fit of adaptive threshold algorithms worsened significantly when RM algorithms were employed for both scanners. CONCLUSIONS: Differences in general quantitative performance were observed for the PET scanners analyzed. Segmentation of PET images using adaptive threshold algorithms should not be undertaken in conjunction with RM reconstructions.

An Adaptive Thresholding Method for BTV Estimation Incorporating PET Reconstruction Parameters: A Multicenter Study of the Robustness and the Reliability.

OBJECTIVE: The aim of this work was to assess robustness and reliability of an adaptive thresholding algorithm for the biological target volume estimation incorporating reconstruction parameters. METHOD: In a multicenter study, a phantom with spheres of different diameters (6.5-57.4 mm) was filled with (18)F-FDG at different target-to-background ratios (TBR: 2.5-70) and scanned for different acquisition periods (2-5 min). Image reconstruction algorithms were used varying number of iterations and postreconstruction transaxial smoothing. Optimal thresholds (TS) for volume estimation were determined as percentage of the maximum intensity in the cross section area of the spheres. Multiple regression techniques were used to identify relevant predictors of TS. RESULTS: The goodness of the model fit was high (R(2): 0.74-0.92). TBR was the most significant predictor of TS. For all scanners, except the Gemini scanners, FWHM was an independent predictor of TS. Significant differences were observed between scanners of different models, but not between different scanners of the same model. The shrinkage on cross validation was small and indicative of excellent reliability of model estimation. CONCLUSIONS: Incorporation of postreconstruction filtering FWHM in an adaptive thresholding algorithm for the BTV estimation allows obtaining a robust and reliable method to be applied to a variety of different scanners, without scanner-specific individual calibration.

Dosimetry study in patients with autonomous thyroid nodule who are candidates for radioiodine therapy.

UNLABELLED: A dosimetry study was performed on 26 patients with an autonomous thyroid nodule and suppressed serum thyroid-stimulating hormone, to determine the dose to extranodular tissue when the nodule receives 300 Gy for 131I therapy. METHODS: Parameters of radioiodine turnover to be used in the dosimetry formula were separately obtained for the nodule and the contralateral lobe, as a measurable example of the extranodular tissue, using 55 MBq 123I and a computer-assisted gamma camera. The biologic half-life of 123I was then converted into the effective half-life of 131I, and the volumes of the nodule and the lobe were obtained by scintigraphy or sonography. RESULTS: The mean dose to the contralateral lobe from uptake and irradiation by the nodule was calculated to be 32 Gy, and that to the ipsilateral lobe was estimated to be 34 Gy. CONCLUSION: During radioiodine therapy for autonomous thyroid nodules, the extranodular tissue receives a higher dose than is generally assumed, which explains the relatively high rate of post-treatment hypothyroidism reported in the literature.

Changes in the energy response of a dedicated gamma camera after exposure to a high-flux irradiation.

This work reports the effects of the gain variation of the photomultiplier tubes (PMTs) observed on a cardiac dedicated gamma camera after accidental high-flux irradiation. One detector of this dual-headed 90 degrees-fixed gamma camera was accidentally left uncollimated during a quality assurance procedure on the other detector with a 57Co flood source (259 MBq) and received a non-uniform high flux of 1.9-0.6 Mcps over 25000 mm2 areas for about 30 min. To evaluate the severity and the duration of the perturbation effect on the energy response of the detector, the photopeak position was monitored for about 1 month with a 99mTc point source. The 140 keV photopeak shifted to 158 keV soon after irradiation, reached the correct position after 9 days and moved to a stable value of 132 keV after 15 days. Afterwards, a new energy calibration reset the photopeak position at 140 keV and the correct energy response of the gamma camera. This experience suggests that particular care should be taken to avoid exposures to high radiation fluxes that induce persistent gain shifts on the PMTs of this system.

Minimum acceptable sensitivity of intraoperative gamma probes used for sentinel lymph node detection in melanoma patients.

The aim of this study was to determine the suspension level for the sensitivity of an intraoperative scintillation gamma probe in the detection of the sentinel lymph node (SLN) in melanoma patients. Thirty-eight consecutive patients with melanoma were enrolled in the study during a 12-month period and underwent lymphatic scintigraphy after the peritumoral intradermal administration of about 14 MBq of (99m)Tc-nanocolloids. The SLNs were successfully removed during the surgical intervention about 4 h later. To identify and localize the SLN, a scintillation NaI(Tl) collimated probe was used. Predictably, the probe sensitivity decreased as the photopeak energy window was progressively narrowed, from 6.9 ± 0.7 counts per second (cps)/kBq (designated as the 'optimum,' or 'OPT,' sensitivity) to 2.5 ± 0.3 cps/kBq (LOW sensitivity) and to 1.4 ± 0.2 cps/kBq (VLOW sensitivity). Maximum lymph node count rates (cps) were determined for the foregoing energy windows prior to skin incision (PREOPT, PRELOW, PREVLOW, respectively) and in vivo after incision (INVOPT, INVLOW, INVVLOW). Forty-three SLNs were removed with a mean source-to-detector distance of 46 ± 24 mm (min 12 mm, max 92 mm). Four SLNs could not have been detected using PRELOW. This figure would have decreased to 34, with nine undetectable lymph nodes, with PREVLOW. One SLN could not have been identified using INVLOW and four could not have be identified using INVVLOW. In the clinical scenario of SLN detection in melanoma patients, a system sensitivity of 2.5 cps/kBq represents a suspension level, that is, a level under which the equipment must be suspended from clinical use and the poor performance must be investigated.

Permeability of gloves used in nuclear medicine departments to [(99m)Tc]-pertechnetate and [(18)F]-fluorodeoxyglucose: radiation protection considerations.

In order to evaluate the safety of the individual protection devices, the permeability of four different types of disposable gloves, commonly used in hospitals, was tested in relation to [(99m)Tc]-pertechnetate and to [(18)F]-fluorodeoxyglucose ([(18)F]-FDG). From these radiopharmaceutical solutions, a drop was deposited on the external surface of the glove which was opened and stretched with the external surface placed upward. The smear test technique permitted to evaluate the activity onto the inner surface of the glove at different times. The smear tests were measured in a well sodium iodide detector calibrated in efficiency for (99m)Tc and (18)F. The permeability was tested on ten samples of each type of gloves and was expressed as the ratio of the activity onto the inner surface at each time interval to the activity deposited on the external surface of the glove. For each type of gloves and for each sampling time, mean value, standard deviation and percentage coefficient of variation of permeability were evaluated. One type of gloves showed a low resistance to permeation of both radiopharmaceuticals, while another one only to pertechnetate. The other gloves were good performers. The results of this study suggest to test permeability for gloves used for handling radiopharmaceuticals, before their adoption in the clinical routine. This practice will provide a more careful service of radiation protection for nuclear medicine department staff.

Five-year experience of quality control for a 3D LSO-based whole-body PET scanner: results and considerations.

PET scanners require routine monitoring and quality control (QC) to ensure proper scanner performance. QC helps to ensure that PET equipment performs as specified by the manufacturer and that there have not been significant changes in the system response since acceptance. In this work we describe the maintenance history and we report on the results obtained from the PET system QC testing program over 5 years at two centers, both utilizing a Siemens Biograph 16 HiRez PET/CT system. QC testing programs were based on international standards and included the manufacturer's daily QC, monthly uniformity and sensitivity, quarterly cross-calibration and annual resolution and image quality. For the Winnipeg and Novara sites, two and one PET detector blocks have been replaced, respectively. Neither system has had other significant PET system related hardware replacements. The manufacturer's suggested daily QC was sensitive to detecting problems in the function of PET detector elements. The same test was not sensitive for detecting long term drifts in the systems: the Novara system observed a significant deterioration over five years of testing in the sensitivity which exhibited a decrease of 16% as compared to its initial value measured at system installation. The measure of the energy spectrum, showed that the 511 keV photopeak had shifted to a position of 468 keV. This shift was corrected by having service personnel perform a complete system calibration and detector block setup. We recommend including tests of system energy response and of sensitivity as part of a QC program since they can provide useful information on the actual performance of the scanner. A modification of the daily QC test by the manufacturer is suggested to monitor the long term stability of the system. Image quality and spatial resolution tests have proven to be of limited value for monitoring the system over time.

Reproducibility and variability of global and regional dyssynchrony parameters derived from phase analysis of gated myocardial perfusion SPECT.

AIM: The aim of the study was to assess the reproducibility (R) and variability of quantitative parameters derived from Cedars-Sinai QGS software for phase analysis of G-SPECT in 39 coronary artery disease (CAD) patients with or without severe heart-failure (HF). METHODS: Phase standard deviation (s) and entropy (e) are global parameters that quantify dyssynchrony. Regional measures of dyssynchrony were also computed including differences between the lateral wall versus the septal wall (DmW) and differences between the averages of midventricular lateral versus midventricular septal segments (DmS). RESULTS: Global parameters e and s and regional parameters DmW and DmS exhibit excellent values of R=0.92, 0.99, 0.99 and 0.96, respectively. In regional parameters DmW and DmS there is a significant variability in individual scoring assessed by a standard error of measurement of =9.17 and 21.7, respectively. The box plots of e in patients with or without HF do not show any significant superimposition, while the box plots of s and DmW show a partial overlap mainly due to the significant variability of s and DmW within patients with HF. Conversely the box plots of DmS in patients with or without HF show a significant overlap due to the great variability of DmS within patients with HF. CONCLUSION: Regional parameters derived from phase analysis of G-SPECT studies are not useful in the individual assessment of dyssynchrony in CAD patients either due to large variability in individual scoring and to a large heterogeneity in HF patients. Global parameters (e and s) exhibit both an excellent reproducibility. Nonetheless, e seems to perform better than s in individual assessment of dyssynchrony due to a better separation between HF and non HF patients. Using the QGS software approach no manual intervention is necessary to ensure a good reproducibility of global parameters.

Characterization of the [(153)Sm]Sm-EDTMP pharmacokinetics and estimation of radiation absorbed dose on an individual basis.

[(153)Sm]Sm-EDTMP is a radiopharmaceutical used in palliation cares of bone metastases. The purpose of this study is to provide an explicit description of [(153)Sm]Sm-EDTMP pharmacokinetics, adopting a simple three-compartmental model with the analytical expressions calculating the rate constants and determining biodistribution parameters, like radiopharmaceutical uptake and clearance. This biokinetic model allowed us to calculate on an individual basis the dose to bone surface and to red bone marrow and to assess the degree of variability in dosimetric parameters using a fixed administered activity based only on patient weight. In this study twenty patients were enrolled and were treated with [(153)Sm]Sm-EDTMP, administering a fixed activity per kilogram (37 MBq/kg); blood and urine samples were collected during 24 h post treatment. The median value of the administered activity was 2.7 GBq. Blood clearance confirmed that an aliquot of [(153)Sm]Sm-EDTMP rapidly localizes and is retained in bone, while the remainder is rapidly cleared from the blood pool by the urinary system. Our data show a bi-exponential clearance from blood: the rapid component has a half life median value of 6 min (range: 2-24 min), while the slow one has a half life median value of 1.4 h (range: 0.6-5.8 h). Median value of the urinary excretion is 40 (range: 3-75) % of the administered activity. Our model shows the behaviour of a tracer which is distributed in the extracellular space of the body, localized in the skeleton and excreted via glomerular filtration. Half life median values of [(153)Sm]Sm-EDTMP transferring between compartments, T(1/2) (blood→ECF), T(1/2) (ECF→blood) are 7.4 (range: 1.9-37) and 48 (range: 8-408) min, respectively. Median values of half lives of [(153)Sm]Sm-EDTMP clearance through the urine and of uptake into bone are 1.0 (range: 0.1-6.0) and 1.6 (range: 0.6-9.0) h, respectively. Median value of red marrow absorbed dose is 2.1 (range: 0.7-3.5) Gy and 0.8 (range: 0.3-2.1) Gy/GBq, while median value of bone surface absorbed dose is 11.5 Gy (range: 5.0-18.4) and 4.4 (range: 2.3-14.3) Gy/GBq. It is remarkable that there is a really great biological variability within patients, especially considering the excreted activity. The cumulated activity in bone and red marrow doses were significantly higher in prostate cancer, where metastatic bone lesions are osteoblastic, than in breast cancer where metastatic bone lesions are osteolytic or mixed (lytic/blastic). The relevant biological variability in biodistribution and metabolism of [(153)Sm]Sm-EDTMP suggests that the fixed administered activity based on patient weight is not sufficient to optimize the treatment and a better optimization would be reached by using a predictive dosimetry tailored to individual patient characteristics.

The effect of activity outside the field of view on image quality for a 3D LSO-based whole body PET/CT scanner.

The purpose of this study was to quantify the influence of outside field of view (FOV) activity concentration (A(c)(,out)) on the noise equivalent count rate (NECR), scatter fraction (SF) and image quality of a 3D LSO whole-body PET/CT scanner. The contrast-to-noise ratio (CNR) was the figure of merit used to characterize the image quality of PET scans. A modified International Electrotechnical Commission (IEC) phantom was used to obtain SF and counting rates similar to those found in average patients. A scatter phantom was positioned at the end of the modified IEC phantom to simulate an activity that extends beyond the scanner. The modified IEC phantom was filled with (18)F (11 kBq mL(-1)) and the spherical targets, with internal diameter (ID) ranging from 10 to 37 mm, had a target-to-background ratio of 10. PET images were acquired with background activity concentrations into the FOV (A(c)(,bkg)) about 11, 9.2, 6.6, 5.2 and 3.5 kBq mL(-1). The emission scan duration (ESD) was set to 1, 2, 3 and 4 min. The tube inside the scatter phantom was filled with activities to provide A(c)(,out) in the whole scatter phantom of zero, half, unity, twofold and fourfold the one of the modified IEC phantom. Plots of CNR versus the various parameters are provided. Multiple linear regression was employed to study the effects of A(c)(,out) on CNR, adjusted for the presence of variables (sphere ID, A(c)(,bkg) and ESD) related to CNR. The presence of outside FOV activity at the same concentration as the one inside the FOV reduces peak NECR of 30%. The increase in SF is marginal (1.2%). CNR diminishes significantly with increasing outside FOV activity, in the range explored. ESD and A(c)(,out) have a similar weight in accounting for CNR variance. Thus, an experimental law that adjusts the scan duration to the outside FOV activity can be devised. Recovery of CNR loss due to an elevated A(c)(,out) activity seems feasible by modulating the ESD in individual bed positions according to A(c)(,out).