Showing papers in "EJNMMI Physics in 2017"

Cerenkov luminescence imaging: physics principles and potential applications in biomedical sciences

Abstract: Cerenkov luminescence imaging (CLI) is a novel imaging modality to study charged particles with optical methods by detecting the Cerenkov luminescence produced in tissue. This paper first describes the physical processes that govern the production and transport in tissue of Cerenkov luminescence. The detectors used for CLI and their most relevant specifications to optimize the acquisition of the Cerenkov signal are then presented, and CLI is compared with the other optical imaging modalities sharing the same data acquisition and processing methods. Finally, the scientific work related to CLI and the applications for which CLI has been proposed are reviewed. The paper ends with some considerations about further perspectives for this novel imaging modality.

Variations in the practice of molecular radiotherapy and implementation of dosimetry: results from a European survey

Abstract: Currently, the implementation of dosimetry in molecular radiotherapy (MRT) is not well investigated, and in view of the Council Directive (2013/59/Euratom), there is a need to understand the current availability of dosimetry-based MRT in clinical practice and research studies. The aim of this study was to assess the current practice of MRT and dosimetry across European countries. An electronic questionnaire was distributed to European countries. This addressed 18 explicitly considered therapies, and for each therapy, a similar set of questions were included. Questions covered the number of patients and treatments during 2015, involvement of medical specialties and medical physicists, implementation of absorbed dose planning, post-therapy imaging and dosimetry, and the basis of therapy prescription. Responses were obtained from 26 countries and 208 hospitals, administering in total 42,853 treatments. The most common therapies were 131I-NaI for benign thyroid diseases and thyroid ablation of adults. The involvement of a medical physicist (mean over all 18 therapies) was reported to be either minority or never by 32% of the responders. The percentage of responders that reported that dosimetry was included on an always/majority basis differed between the therapies and showed a median value of 36%. The highest percentages were obtained for 177Lu-PSMA therapy (100%), 90Y microspheres of glass (84%) and resin (82%), 131I-mIBG for neuroblastoma (59%), and 131I-NaI for benign thyroid diseases (54%). The majority of therapies were prescribed based on fixed-activity protocols. The highest number of absorbed-dose based prescriptions were reported for 90Y microsphere treatments in the liver (64% and 96% of responses for resin and glass, respectively), 131I-NaI treatment of benign thyroid diseases (38% of responses), and for 131I-mIBG treatment of neuroblastoma (18% of responses). There is a wide variation in MRT practice across Europe and for different therapies, including the extent of medical-physicist involvement and the implementation of dosimetry-guided treatments.

Point-spread function reconstructed PET images of sub-centimeter lesions are not quantitative.

Abstract: Background PET image reconstruction methods include modeling of resolution degrading phenomena, often referred to as point-spread function (PSF) reconstruction. The aim of this study was to develop a clinically relevant phantom and characterize the reproducibility and accuracy of high-resolution PSF reconstructed images of small lesions, which is a prerequisite for using PET in the prediction and evaluation of responses to treatment. Sets of small homogeneous 18F-spheres (range 3–12 mm diameter, relevant for small lesions and lymph nodes) were suspended and covered by a 11C-silicone, which provided a scattering medium and a varying sphere-to-background ratio. Repeated measurements were made on PET/CT scanners from two vendors using a wide range of reconstruction parameters. Recovery coefficients (RCs) were measured for clinically used volume-of-interest definitions.

Dosimetry-based treatment planning for molecular radiotherapy: a summary of the 2017 report from the Internal Dosimetry Task Force

Abstract: The European directive on basic safety standards (Council directive 2013/59 Euratom) mandates dosimetry-based treatment planning for radiopharmaceutical therapies. The directive comes into operation February 2018, and the aim of a report produced by the Internal Dosimetry Task Force of the European Association of Nuclear Medicine is to address this aspect of the directive. A summary of the report is presented. A brief review of five of the most common therapy procedures is included in the current text, focused on the potential to perform patient-specific dosimetry. In the full report, 11 different therapeutic procedures are included, allowing additional considerations of effectiveness, references to specific literature on quantitative imaging and dosimetry, and existing evidence for absorbed dose-effect correlations for each treatment. Individualized treatment planning with tracer diagnostics and verification of the absorbed doses delivered following therapy is found to be scientifically feasible for almost all procedures investigated, using quantitative imaging and/or external monitoring. Translation of this directive into clinical practice will have significant implications for resource requirements. Molecular radiotherapy is undergoing a significant expansion, and the groundwork for dosimetry-based treatment planning is already in place. The mandated individualization is likely to improve the effectiveness of the treatments, although must be adequately resourced.

Accuracy of (177)Lu activity quantification in SPECT imaging: a phantom study.

Abstract: The 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. Quantification error for scans in air was better for TEW scatter correction ( 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. Following 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.

Saul Hertz, MD, and the birth of radionuclide therapy.

Abstract: The year, 2016, marked the 75th anniversary of Dr. Saul Hertz first using radioiodine to treat a patient with thyroid disease. In November of 1936, a luncheon was held of the faculty of Harvard Medical School where Karl Compton, PhD, president of the Massachusetts Institute of Technology was invited to give a presentation entitled "What Physics Can Do for Biology and Medicine." Saul Hertz who attended the luncheon spontaneously asked the very pertinent question that perhaps changed the course of treatment of thyroid disease, "Could iodine be made radioactive artificially?" We review the events leading up to the asking of this question, the preclinical investigations by Dr. Hertz and his colleague Arthur Roberts prior to the treatment of the first patient and what occurred in the years following this landmark event. This commentary seeks to set the record straight to the sequence of events leading to the first radioiodine therapy, so that those involved can be recognized with due credit.

[(123)I]FP-CIT ENC-DAT normal database: the impact of the reconstruction and quantification methods

Abstract: [123I]FP-CIT is a well-established radiotracer for the diagnosis of dopaminergic degenerative disorders. The European Normal Control Database of DaTSCAN (ENC-DAT) of healthy controls has provided age and gender-specific reference values for the [123I]FP-CIT specific binding ratio (SBR) under optimised protocols for image acquisition and processing. Simpler reconstruction methods, however, are in use in many hospitals, often without implementation of attenuation and scatter corrections. This study investigates the impact on the reference values of simpler approaches using two quantifications methods, BRASS and Southampton, and explores the performance of the striatal phantom calibration in their harmonisation. BRASS and Southampton databases comprising 123 ENC-DAT subjects, from gamma cameras with parallel collimators, were reconstructed using filtered back projection (FBP) and iterative reconstruction OSEM without corrections (IRNC) and compared against the recommended OSEM with corrections for attenuation and scatter and septal penetration (ACSC), before and after applying phantom calibration. Differences between databases were quantified using the percentage difference of their SBR in the dopamine transporter-rich striatum, with their significance determined by the paired t test with Bonferroni correction. Attenuation and scatter losses, measured from the percentage difference between IRNC and ACSC databases, were of the order of 47% for both BRASS and Southampton quantifications. Phantom corrections were able to recover most of these losses, but the SBRs remained significantly lower than the “true” values (p < 0.001). Calibration provided, in fact, “first order” camera-dependent corrections, but could not include “second order” subject-dependent effects, such as septal penetration from extra-cranial activity. As for the ACSC databases, phantom calibration was instrumental in compensating for partial volume losses in BRASS (~67%, p < 0.001), while for the Southampton method, inherently free from them, it brought no significant changes and solely corrected for residual inter-camera variability (−0.2%, p = 0.44). The ENC-DAT reference values are significantly dependent on the reconstruction and quantification methods and phantom calibration, while reducing the major part of their differences, is unable to fully harmonize them. Clinical use of any normal database, therefore, requires consistency with the processing methodology. Caution must be exercised when comparing data from different centres, recognising that the SBR may represent an “index” rather than a “true” value.

Minimum lesion detectability as a measure of PET system performance

Abstract: A phantom in combination with an imaging protocol was developed to measure the limit of small lesion detection on different PET systems. Seven small spheres with inner diameters ranging from 3.95 up to 15.43 mm were imaged in a Jaszczak ECT Phantom, in air, in a cold background, and with sphere to background contrast ratios of 15:1 down to 1.88:1. The imaging times varied from 1 to 16 min. The imaging protocol was performed on the Gemini TF and Vereos by Philips, the mCT and HRRT by Siemens, and the Discovery 710 by General Electric. For each scanning condition, the images were reconstructed with image voxel sizes of 1 to 4 mm cubic voxels. The reconstruction method used for each system was the one recommended by the manufacture to achieve best small image lesion detection results. A human observer study was performed to determine the smallest observable sphere for each scanning condition. All systems were able to image the smallest sphere of 3.95 mm inner diameter at the 15 to 1 signal to background ratio when imaged for 16 min. For a typical whole body per bed position scan time of 2 to 4 min, the smallest imaged sphere varied between 4.95 and 6.23 mm at the 15:1 contrast ratio and 12.43 and 15.43 mm at a contrast ratio of 1.88:1. In general, all systems were consistent with the Rose criteria when determining lesion detectability. Besides demonstrating that the current state of the art clinical PET/CT systems have the same lesion detection ability, the study demonstrates how sensitive scan time can be to detecting small lesions which have a relatively small contrast uptake in the range of just 2:1. This should help guide imaging protocols to use longer scan times over regions of the subject in which small lesions are suspect.

Comparison of machine learning and semi-quantification algorithms for (I123)FP-CIT classification: the beginning of the end for semi-quantification?

Abstract: Semi-quantification methods are well established in the clinic for assisted reporting of (I123) Ioflupane images. Arguably, these are limited diagnostic tools. Recent research has demonstrated the potential for improved classification performance offered by machine learning algorithms. A direct comparison between methods is required to establish whether a move towards widespread clinical adoption of machine learning algorithms is justified. This study compared three machine learning algorithms with that of a range of semi-quantification methods, using the Parkinson’s Progression Markers Initiative (PPMI) research database and a locally derived clinical database for validation. Machine learning algorithms were based on support vector machine classifiers with three different sets of features: Semi-quantification methods were based on striatal binding ratios (SBRs) from both putamina, with and without consideration of the caudates. Normal limits for the SBRs were defined through four different methods: Each machine learning and semi-quantification technique was evaluated with stratified, nested 10-fold cross-validation, repeated 10 times. The mean accuracy of the semi-quantitative methods for classification of local data into Parkinsonian and non-Parkinsonian groups varied from 0.78 to 0.87, contrasting with 0.89 to 0.95 for classifying PPMI data into healthy controls and Parkinson’s disease groups. The machine learning algorithms gave mean accuracies between 0.88 to 0.92 and 0.95 to 0.97 for local and PPMI data respectively. Classification performance was lower for the local database than the research database for both semi-quantitative and machine learning algorithms. However, for both databases, the machine learning methods generated equal or higher mean accuracies (with lower variance) than any of the semi-quantification approaches. The gain in performance from using machine learning algorithms as compared to semi-quantification was relatively small and may be insufficient, when considered in isolation, to offer significant advantages in the clinical context.

Model selection criteria for dynamic brain PET studies

Abstract: Several criteria exist to identify the optimal model for quantification of tracer kinetics. The purpose of this study was to evaluate the correspondence in kinetic model preference identification for brain PET studies among five model selection criteria: Akaike Information Criterion (AIC), AIC unbiased (AICC), model selection criterion (MSC), Schwartz Criterion (SC), and F-test. Six tracers were evaluated: [11C]FMZ, [11C]GMOM, [11C]PK11195, [11C]Raclopride, [18F]FDG, and [11C]PHT, including data from five subjects per tracer. Time activity curves (TACs) were analysed using six plasma input models: reversible single-tissue model (1T2k), irreversible two-tissue model (2T3k), and reversible two-tissue model (2T4k), all with and without blood volume fraction parameter (V B). For each tracer and criterion, the percentage of TACs preferring a certain model was calculated. For all radiotracers, strong agreement was seen across the model selection criteria. The F-test was considered as the reference, as it is a frequently used hypothesis test. The F-test confirmed the AIC preferred model in 87% of all cases. The strongest (but minimal) disagreement across regional TACs was found when comparing AIC with AICC. Despite these regional discrepancies, same preferred kinetic model was obtained using all criteria, with an exception of one FMZ subject. In conclusion, all five model selection criteria resulted in similar conclusions with only minor differences that did not affect overall model selection.

MLAA-based attenuation correction of flexible hardware components in hybrid PET/MR imaging

Abstract: Accurate PET quantification demands attenuation correction (AC) for both patient and hardware attenuation of the 511 keV annihilation photons. In hybrid PET/MR imaging, AC for stationary hardware components such as patient table and MR head coil is straightforward, employing CT-derived attenuation templates. AC for flexible hardware components such as MR-safe headphones and MR radiofrequency (RF) surface coils is more challenging. Registration-based approaches, aligning CT-based attenuation templates with the current patient position, have been proposed but are not used in clinical routine. Ignoring headphone or RF coil attenuation has been shown to result in regional activity underestimation values of up to 18%. We propose to employ the maximum-likelihood reconstruction of attenuation and activity (MLAA) algorithm to estimate the attenuation of flexible hardware components. Starting with an initial attenuation map not including flexible hardware components, the attenuation update of MLAA is applied outside the body outline only, allowing to estimate hardware attenuation without modifying the patient attenuation map. Appropriate prior expectations on the attenuation coefficients are incorporated into MLAA. The proposed method is investigated for non-TOF PET phantom and 18F-FDG patient data acquired with a clinical PET/MR device, using headphones or RF surface coils as flexible hardware components. Although MLAA cannot recover the exact physical shape of the hardware attenuation maps, the overall attenuation of the hardware components is accurately estimated. Therefore, the proposed algorithm significantly improves PET quantification. Using the phantom data, local activity underestimation when neglecting hardware attenuation was reduced from up to 25% to less than 3% under- or overestimation as compared to reference scans without hardware present or to CT-derived AC. For the patient data, we found an average activity underestimation of 7.9% evaluated in the full brain and of 6.1% for the abdominal region comparing the uncorrected case with MLAA. MLAA is able to provide accurate estimations of the attenuation of flexible hardware components and can therefore be used to significantly improve PET quantification. The proposed approach can be readily incorporated into clinical workflow.

Calculation of left ventricular volumes and ejection fraction from dynamic cardiac-gated 15O-water PET/CT: 5D-PET.

Abstract: Quantitative measurement of myocardial blood flow (MBF) is of increasing interest in the clinical assessment of patients with suspected coronary artery disease (CAD). 15O-water positron emission tomography (PET) is considered the gold standard for non-invasive MBF measurements. However, calculation of left ventricular (LV) volumes and ejection fraction (EF) is not possible from standard 15O-water uptake images. The purpose of the present work was to investigate the possibility of calculating LV volumes and LVEF from cardiac-gated parametric blood volume (V B) 15O-water images and from first pass (FP) images. Sixteen patients with mitral or aortic regurgitation underwent an eight-gate dynamic cardiac-gated 15O-water PET/CT scan and cardiac MRI. V B and FP images were generated for each gate. Calculations of end-systolic volume (ESV), end-diastolic volume (EDV), stroke volume (SV) and LVEF were performed with automatic segmentation of V B and FP images, using commercially available software. LV volumes and LVEF were calculated with surface-, count-, and volume-based methods, and the results were compared with gold standard MRI. Using V B images, high correlations between PET and MRI ESV (r = 0.89, p < 0.001), EDV (r = 0.85, p < 0.001), SV (r = 0.74, p = 0.006) and LVEF (r = 0.72, p = 0.008) were found for the volume-based method. Correlations for FP images were slightly, but not significantly, lower than those for V B images when compared to MRI. Surface- and count-based methods showed no significant difference compared with the volume-based correlations with MRI. The volume-based method showed the best agreement with MRI with no significant difference on average for EDV and LVEF but with an overestimation of values for ESV (14%, p = 0.005) and SV (18%, p = 0.004) when using V B images. Using FP images, none of the parameters showed a significant difference from MRI. Inter-operator repeatability was excellent for all parameters (ICC > 0.86, p < 0.001). Calculation of LV volumes and LVEF from dynamic 15O-water PET is feasible and shows good correlation with MRI. However, the analysis method is laborious, and future work is needed for more automation to make the method more easily applicable in a clinical setting.

3D absorbed dose distribution estimated by Monte Carlo simulation in radionuclide therapy with a monoclonal antibody targeting synovial sarcoma

Abstract: Radiolabeled OTSA101, a monoclonal antibody targeting synovial sarcoma (SS) developed by OncoTherapy Science, was used to treat relapsing SS metastases following a theranostic procedure: in case of significant 111In-OTSA101 tumor uptake and favorable biodistribution, patient was randomly treated with 370/1110 MBq 90Y-OTSA101. Monte Carlo-based 3D dosimetry integrating time-activity curves in VOI was performed on 111In-OTSA101 repeated SPECT/CT. Estimated absorbed doses (AD) in normal tissues were compared to biological side effects and to the admitted maximal tolerated absorbed dose (MTD) in normal organs. Results in the tumors were also compared to disease evolution. Biodistribution and tracer quantification were analyzed on repeated SPECT/CT acquisitions performed after injection of 111In-OTSA101 in 19/20 included patients. SPECT images were warped to a common coordinates system with deformable registration. Volumes of interest (VOI) for various lesions and normal tissues were drawn on the first CT acquisition and reported to all the SPECT images. Tracer quantification and residence time of 111In-OTSA101 in VOI were used to evaluate the estimated absorbed doses per MBq of 90Y-OTSA101 by means of Monte Carlo simulations (GATE). A visual scale analysis was applied to assess tumor uptake (grades 0 to 4) and results were compared to the automated quantification. Results were then compared to biological side effects reported in the selected patients treated with 90Y-OTSA101 but also to disease response to treatment. After screening, 8/20 patients were treated with 370 or 1110 MBq 90Y-OTSA101. All demonstrated medullary toxicity, only one presented with transient grade 3 liver toxicity due to disease progression, and two patients presented with transient grade 1 renal toxicity. Median absorbed doses were the highest in the liver (median, 0.64 cGy/MBq; [0.27 −1.07]) being far lower than the 20 Gy liver MTD, and the lowest in bone marrow (median, 0.09 cGy/MBq; [0.02 −0.18]) being closer to the 2 Gy bone marrow MTD. Most of the patients demonstrated progressive disease on RECIST criteria during patient follow-up. 111In-OTSA101 tumors tracer uptake visually appeared highly heterogeneous in inter- and intra-patient analyses, independently of tumor sizes, with variable kinetics. The majority of visual grades corresponded to the automated computed ones. Estimated absorbed doses in the 95 supra-centimetric selected lesions ranged from 0.01 to 0.71 cGy per injected MBq (median, 0.22 cGy/MBq). The maximal tumor AD obtained was 11.5 Gy. 3D dosimetry results can explain the observed toxicity and tumors response. Despite an intense visual 111In-OTSA101 liver uptake, liver toxicity was not the dose limiting factor conversely to bone marrow toxicity. Even though tumors 111In-OTSA101 avidity was visually obvious for treated patients, the low estimated tumors AD obtained by 3D dosimetry explain the lack of tumor response.

Accuracy and precision assessment for activity quantification in individualized dosimetry of 177 Lu-DOTATATE therapy

Abstract: In order to obtain a reliable 177Lu-DOTATATE therapy dosimetry, it is crucial to acquire accurate and precise activity measurements with the radionuclide calibrator, the SPECT/CT camera, and the NaI(Tl) well counter. The aim of this study was to determine, in a clinical context, the accuracy and the precision of their activity quantification over a range of activities and time. Ninety-three 177Lu sources from the manufacturer were measured in the radionuclide calibrator over 2.5 years to evaluate its calibration accuracy and precision compared to the manufacturer’s value. A NEMA 2012/IEC 2008 phantom was filled with a 177Lu activity concentration sphere-to-background ratio of five. It was acquired with the SPECT/CT camera to determine the reconstruction parameters offering the best compromise between partial volume effect and signal-to-noise ratio. The calibration factor was computed accordingly. The calibration quality was monitored over 2.5 years with 33 phantom acquisitions with activities ranging from 7040 to 0.6 MBq. Home-made sources were used to calibrate the well counter. Its reliability was evaluated with activities ranging from 150 to 0.2 kBq measured 34 times over 2.5 years. For the radionuclide calibrator, median [interquartile range] for the error on activity measurement was −0.99 [1.31] %. The optimal SPECT reconstruction parameters were obtained with 16 iterations, 16 subsets and a 12-mm Gaussian post-filter. The calibration factor was 9.87 cps/MBq with an error of −1.05 [2.12] %. The well counter was calibrated with 31.5 cps/kBq, and the error was evaluated to −12.89 [16.55] %. The accuracy and the precision of activity quantification using dedicated quality control were found to be sufficient for use in dosimetry implemented in clinical routine. The proposed methodology could be implemented in other centres to obtain reproducible 177Lu-based treatment dosimetry.

Guidelines for quality control of PET/CT scans in a multicenter clinical study

Abstract: To date, there is no published detailed checklist with parameters referencing the DICOM tag information with respect to the quality control (QC) of PET/CT scans. The aims of these guidelines are to provide the know-how for effectively controlling the quality of PET/CT scans in multicenter studies, to standardize the QC, to give sponsors and regulatory agencies a basis for justification of the data quality when using standardized uptake values as an imaging biomarker, to document the compliance with the imaging guidelines, to verify the per protocol population versus intent to treat population, and to safeguard the validity of multicenter study conclusions employing standardized uptake value (SUV) as an imaging biomarker which is paramount to the scientific community. Following the proposed guidelines will ensure standardized prospective imaging QC of scans applicable to most studies where SUVs are used as an imaging biomarker. The multitude of factors affecting SUV measurements when not controlled inflicts noise on the data. Decisions on patient management with substantial noise would be devastating to patients, ultimately undermine treatment outcome, and invalidate the utility of SUV as an imaging biomarker usefulness. Strict control of the data quality used for the validation of SUV as an imaging biomarker would ensure trust and reliability of the data.

Quantitative myocardial blood flow imaging with integrated time-of-flight PET-MR.

Abstract: The use of integrated PET-MR offers new opportunities for comprehensive assessment of cardiac morphology and function. However, little is known on the quantitative accuracy of cardiac PET imaging with integrated time-of-flight PET-MR. The aim of the present work was to validate the GE Signa PET-MR scanner for quantitative cardiac PET perfusion imaging. Eleven patients (nine male; mean age 59 years; range 46–74 years) with known or suspected coronary artery disease underwent 15O-water PET scans at rest and during adenosine-induced hyperaemia on a GE Discovery ST PET-CT and a GE Signa PET-MR scanner. PET-MR images were reconstructed using settings recommended by the manufacturer, including time-of-flight (TOF). Data were analysed semi-automatically using Cardiac VUer software, resulting in both parametric myocardial blood flow (MBF) images and segment-based MBF values. Correlation and agreement between PET-CT-based and PET-MR-based MBF values for all three coronary artery territories were assessed using regression analysis and intra-class correlation coefficients (ICC). In addition to the cardiac PET-MR reconstruction protocol as recommended by the manufacturer, comparisons were made using a PET-CT resolution-matched reconstruction protocol both without and with TOF to assess the effect of time-of-flight and reconstruction parameters on quantitative MBF values. Stress MBF data from one patient was excluded due to movement during the PET-CT scanning. Mean MBF values at rest and stress were (0.92 ± 0.12) and (2.74 ± 1.37) mL/g/min for PET-CT and (0.90 ± 0.23) and (2.65 ± 1.15) mL/g/min for PET-MR (p = 0.33 and p = 0.74). ICC between PET-CT-based and PET-MR-based regional MBF was 0.98. Image quality was improved with PET-MR as compared to PET-CT. ICC between PET-MR-based regional MBF with and without TOF and using different filter and reconstruction settings was 1.00. PET-MR-based MBF values correlated well with PET-CT-based MBF values and the parametric PET-MR images were excellent. TOF and reconstruction settings had little impact on MBF values.

Digital PET compliance to EARL accreditation specifications

Abstract: Our aim was to evaluate if a recently introduced TOF PET system with digital photon counting technology (Philips Healthcare), potentially providing an improved image quality over analogue systems, can fulfil EANM research Ltd (EARL) accreditation specifications for tumour imaging with FDG-PET/CT. We have performed a phantom study on a digital TOF PET system using a NEMA NU2-2001 image quality phantom with six fillable spheres. Phantom preparation and PET/CT acquisition were performed according to the European Association of Nuclear Medicine (EANM) guidelines. We made list-mode ordered-subsets expectation maximization (OSEM) TOF PET reconstructions, with default settings, three voxel sizes (4 × 4 × 4 mm3, 2 × 2 × 2 mm3 and 1 × 1 × 1 mm3) and with/without point spread function (PSF) modelling. On each PET dataset, mean and maximum activity concentration recovery coefficients (RCmean and RCmax) were calculated for all phantom spheres and compared to EARL accreditation specifications. The RCs of the 4 × 4 × 4 mm3 voxel dataset without PSF modelling proved closest to EARL specifications. Next, we added a Gaussian post-smoothing filter with varying kernel widths of 1–7 mm. EARL specifications were fulfilled when using kernel widths of 2 to 4 mm. TOF PET using digital photon counting technology fulfils EARL accreditation specifications for FDG-PET/CT tumour imaging when using an OSEM reconstruction with 4 × 4 × 4 mm3 voxels, no PSF modelling and including a Gaussian post-smoothing filter of 2 to 4 mm.

EORTC PET response criteria are more influenced by reconstruction inconsistencies than PERCIST but both benefit from the EARL harmonization program.

Abstract: This study evaluates the consistency of PET evaluation response criteria in solid tumours (PERCIST) and European Organisation for Research and Treatment of Cancer (EORTC) classification across different reconstruction algorithms and whether aligning standardized uptake values (SUVs) to the European Association of Nuclear Medicine acquisition (EANM)/EARL standards provides more consistent response classification Baseline (PET1) and response assessment (PET2) scans in 61 patients with non-small cell lung cancer were acquired in protocols compliant with the EANM guidelines and were reconstructed with point-spread function (PSF) or PSF + time-of-flight (TOF) reconstruction for optimal tumour detection and with a standardized ordered subset expectation maximization (OSEM) reconstruction known to fulfil EANM harmonizing standards Patients were recruited in three centres Following reconstruction, EQPET, a proprietary software solution was applied to the PSF ± TOF data (PSF ± TOFEQ) to harmonize SUVs to the EANM standards The impact of differing reconstructions on PERCIST and EORTC classification was evaluated using standardized uptake values corrected for lean body mass (SUL) Using OSEMPET1/OSEMPET2 (standard scenario), responders displayed a reduction of −575% ± 234 and −639% ± 224 for SULmax and SULpeak, respectively, while progressing tumours had an increase of +634% ± 265 and +607% ± 196 for SULmax and SULpeak respectively The use of PSF ± TOF reconstruction impacted the classification of tumour response For example, taking the OSEMPET1/PSF ± TOFPET2 scenario reduced the apparent reduction in SUL in responding tumours (−397% ± 313 and −555% ± 263 for SULmax and SULpeak, respectively) but increased the apparent increase in SUL in progressing tumours (+1300% ± 507 and +911% ± 396 for SULmax and SULpeak, respectively) Consequently, variation in reconstruction methodology (PSF ± TOFPET1/OSEMPET2 or OSEM PET1/PSF ± TOFPET2) led, respectively, to 11/61 (180%) and 10/61 (164%) PERCIST classification discordances and to 17/61 (289%) and 19/61 (311%) EORTC classification discordances An agreement was better for these scenarios with application of the propriety filter, with kappa values of 100 and 095 compared to 075 and 077 for PERCIST and kappa values of 093 and 095 compared to 061 and 055 for EORTC, respectively PERCIST classification is less sensitive to reconstruction algorithm-dependent variability than EORTC classification but harmonizing SULs within the EARL program is equally effective with either

Depth of interaction determination in monolithic scintillator with double side SiPM readout

Abstract: Monolithic scintillators read out by arrays of photodetectors represent a promising solution to obtain high spatial resolution and the depth of interaction (DOI) of the annihilation photon. We have recently investigated a detector geometry composed of a monolithic scintillator readout on two sides by silicon photomultiplier (SiPM) arrays, and we have proposed two parameters for the DOI determination: the difference in the number of triggered SiPMs on the two sides of the detector and the difference in the maximum collected signal on a single SiPM on each side. This work is focused on the DOI calibration and on the determination of the capability of our detector. For the DOI calibration, we studied a method which can be implemented also in detectors mounted in a full PET scanner. We used a PET detector module composed of a monolithic 20 × 20 × 10 mm3 LYSO scintillator crystal coupled on two opposite faces to two arrays of SiPMs. On each side, the scintillator was coupled to 6 × 6 SiPMs. In this paper, the two parameters previously proposed for the DOI determination were calibrated with two different methods. The first used a lateral scan of the detector with a collimated 511 keV pencil beam at steps of 0.5 mm to study the detector DOI capability, while the second used the background radiation of the 176Lu in the scintillator. The DOI determination capability was tested on different regions of the detector using each parameter and the combination of the two. With both parameters for the DOI determination, in the lateral scan, the bias between the mean reconstructed DOI and the real beam position was lower than 0.3 mm, and the DOI distribution had a standard deviation of about 1.5 mm. When using the calibration with the radioactivity of the LYSO, the mean bias increased of about 0.2 mm but with no degradation of the standard deviation of the DOI distribution. The two parameters allow to achieve a DOI resolution comparable with the state of the art, giving a continuous information about the three-dimensional interaction position of the scintillation. These results were obtained by using simple estimators and a detector scalable to a whole PET system. The DOI calibration obtained using lutetium natural radioactivity gives results comparable to the other standard method but appears more readily applicable to detectors mounted in a full PET scanner.

The effect of high count rates on cardiac perfusion quantification in a simultaneous PET-MR system using a cardiac perfusion phantom

Abstract: PET-MRI is under investigation as a new strategy for quantitative myocardial perfusion imaging. Consideration is required as to the maximum scanner count rate in order to limit dead-time losses resulting from administered activity in the scanner field of view during the first pass of the radiotracer. We performed a decaying-source experiment to investigate the high count-rate performance of a PET-MR system (Siemens mMR) over the expected range of activities during a clinical study. We also performed imaging of a cardiac perfusion phantom, which provides an experimental simulation of clinical transit of a simultaneous radiotracer (phantom injected activities range 252 to 997 MBq) and gadolinium-based contrast agent (GBCA). Time-activity and time-intensity curves of the aorta and myocardium compartments from PET and MR images were determined, and quantification of perfusion was then performed using a standard cardiac kinetic model. The decaying-source experiment showed a maximum noise equivalent count rate (NECRmax) of 286 kcps at a singles rate of 47.1 Mcps. NECR was maintained within 5% (NECR95%) of the NECRmax with a singles rate of 34.1 Mcps, corresponding to 310 MBq in the phantom. Count-rate performance was degraded above the singles rate of 64.9 Mcps due to the number of detection events impacting the quantitative accuracy of reconstructed images. A 10% bias in image activity concentration was observed between singles rates of 78.2 and 82.9 Mcps. Perfusion phantom experiments showed that image-based activity concentration and quantified values of perfusion were affected by count losses when the total singles rate was greater than 64.9 Mcps. This occurred during the peak arterial input function (AIF) phase of imaging for injected activities to the phantom of 600 MBq and greater. Care should be taken to avoid high count-rate losses in simultaneous PET-MRI studies. Based on our results in phantoms, bias in reconstructed images should be avoided by adhering to a singles rate lower than 64.9 Mcps on the mMR system. Quantification of perfusion values using singles rates higher than 64.9 Mcps on this system may be compromised and should be avoided.

Ejection fraction in myocardial perfusion imaging assessed with a dynamic phantom: comparison between IQ-SPECT and LEHR

Abstract: Developments in single photon emission tomography instrumentation and reconstruction methods present a potential for decreasing acquisition times. One of such recent options for myocardial perfusion imaging (MPI) is IQ-SPECT. This study was motivated by the inconsistency in the reported ejection fraction (EF) and left ventricular (LV) volume results between IQ-SPECT and more conventional low-energy high-resolution (LEHR) collimation protocols. IQ-SPECT and LEHR quantitative results were compared while the equivalent number of iterations (EI) was varied. The end-diastolic (EDV) and end-systolic volumes (ESV) and the derived EF values were investigated. A dynamic heart phantom was used to produce repeatable ESVs, EDVs and EFs. Phantom performance was verified by comparing the set EF values to those measured from a gated multi-slice X-ray computed tomography (CT) scan (EFTrue). The phantom with an EF setting of 45, 55, 65 and 70% was imaged with both IQ-SPECT and LEHR protocols. The data were reconstructed with different EI, and two commonly used clinical myocardium delineation software were used to evaluate the LV volumes. The CT verification showed that the phantom EF settings were repeatable and accurate with the EFTrue being within 1% point from the manufacture’s nominal value. Depending on EI both MPI protocols can be made to produce correct EF estimates, but IQ-SPECT protocol produced on average 41 and 42% smaller EDV and ESV when compared to the phantom’s volumes, while LEHR protocol underestimated volumes by 24 and 21%, respectively. The volume results were largely similar between the delineation methods used. The reconstruction parameters can greatly affect the volume estimates obtained from perfusion studies. IQ-SPECT produces systematically smaller LV volumes than the conventional LEHR MPI protocol. The volume estimates are also software dependent.

Design and implementation of a prototype head and neck phantom for the performance evaluation of gamma imaging systems

Abstract: A prototype anthropomorphic head and neck phantom has been designed to simulate the adult head and neck anatomy including some internal organs and tissues of interest, such as thyroid gland and sentinel lymph nodes (SLNs). The design of the head and neck phantom includes an inner jig holding the simulated SLNs and thyroid gland. The thyroid gland structure was manufactured using three-dimensional (3D) printing taking into consideration the morphology and shape of a healthy adult thyroid gland. The head and neck phantom was employed to simulate a situation where there are four SLNs distributed at two different vertical levels and at two depths within the neck. Contrast to noise ratio (CNR) calculations were performed for the detected SLNs at an 80 mm distance between both pinhole collimators (0.5 and 1.0 mm diameters) and the surface of the head and neck phantom with a 100 s acquisition time. The recorded CNR values for the simulated SLNs are higher when the hybrid gamma camera (HGC) was fitted with the 1.0 mm diameter pinhole collimator. For instance, the recorded CNR values for the superficially simulated SLN (15 mm depth) containing 0.1 MBq of 99mTc using 0.5 and 1.0 mm diameter pinhole collimators are 6.48 and 16.42, respectively (~87% difference). Gamma and hybrid optical images were acquired using the HGC for the simulated thyroid gland. The count profiles through the middle of the simulated thyroid gland images provided by both pinhole collimators were obtained. The HGC could clearly differentiate the individual peaks of both thyroid lobes in the gamma image produced by the 0.5-mm pinhole collimator. In contrast, the recorded count profile for the acquired image using the 1.0-mm-diameter pinhole collimator showed broader peaks for both lobes, reflecting the degradation of the spatial resolution with increasing the diameter of the pinhole collimator. This anthropomorphic head and neck phantom provides a valuable tool for assessing the imaging ability of gamma cameras used for imaging the head and neck region. The standardisation of test phantoms for SFOV gamma systems will provide an opportunity to collect data across various medical centres. The phantom described is cost effective, reproducible, flexible and anatomically representative.

Monte Carlo-based quantitative pinhole SPECT reconstruction using a ray-tracing back-projector

Abstract: Monte Carlo simulations provide accurate models of nuclear medicine imaging systems as they can properly account for the full physics of photon transport. The accuracy of the model included in the maximum-likelihood–expectation-maximization (ML-EM) reconstruction limits the overall accuracy of the reconstruction results. In this paper, we present a Monte Carlo-based ML-EM reconstruction method for pinhole single-photon emission computed tomography (SPECT) that has been incorporated into the SIMIND Monte Carlo program. The Monte Carlo-based model, which accounts for all of the physical and geometrical characteristics of the camera system, is used in the forward-projection step of the reconstruction, while a simpler model based on ray-tracing is used for back-projection. The aim of this work was to investigate the quantitative accuracy of this combination of forward- and back-projectors in the clinical pinhole camera GE Discovery NM 530c. The total activity was estimated in 99mTc-filled spheres with volumes between 0.5 and 16 mL. The total sphere activity was generally overestimated but remained within 10% of the reference activity defined by the phantom preparation. The recovered activity converged towards the reference activity as the number of iterations increased. Furthermore, the recovery of the activity concentrations within the physical boundaries of the spheres increased with increasing sphere volume. Additionally, the Monte Carlo-based reconstruction enabled recovery of the true activity concentration in the myocardium of a cardiac phantom mounted in a torso phantom regardless of whether the torso was empty or water-filled. A qualitative comparison to data reconstructed using the clinical reconstruction algorithm showed that the two methods performed similarly, although the images reconstructed using the clinical software were more uniform due to the incorporation of noise regularization and post-filtration in that reconstruction technique. We developed a Monte Carlo-based reconstruction method for pinhole SPECT and evaluated it using phantom measurements. The combination of a Monte Carlo-based forward-projector and a simplified analytical ray-tracing back-projector produced quantitative images of acceptable image quality. No explicit calibration is necessary in this method since the forward-projector model maintains a relationship between the number of counts and activity.

Cardiac PET/CT With Rb-82: Optimization of Image Acquisition and Reconstruction Parameters

Abstract: Our aim was to characterize the influence of time-of-flight (TOF) and point spread function (PSF) recovery corrections, as well as ordered subset expectation maximization (OSEM) reconstruction parameters, in 82Rb PET/CT quantification of myocardial blood flow (MBF) and myocardial flow reserve (MFR). Rest and stress list-mode dynamic 82Rb PET acquisition data from 10 patients without myocardial flow defects and 10 patients with myocardial blood flow defects were reconstructed retrospectively. OSEM reconstructions were performed with Gaussian filters of 4, 6, and 8 mm, different iterations, and subset numbers (2 × 24; 2 × 16; 3 × 16; 4 × 16). Rest and stress global, regional, and segmental MBF and MFR were computed from time activity curves with FlowQuant© software. Left ventricular segmentation using the 17-segment American Heart Association model was obtained. Whole left ventricle (LV) MBF at rest and stress were 0.97 ± 0.30 and 2.30 ± 1.00 mL/min/g, respectively, and MFR was 2.40 ± 1.13. Concordance was excellent and all reconstruction parameters had no significant impact on MBF, except for the exclusion of TOF which led to significantly decreased concordance in rest and stress MBF in patients with or without perfusion defects on a coronary artery basis and in MFR in patients with perfusion defects. Changes in reconstruction parameters in perfusion 82Rb PET/CT studies influence quantitative MBF analysis. The inclusion of TOF information in the tomographic reconstructions had significant impact in MBF quantification.

Comparison of PET/CT and whole-mount histopathology sections of the human prostate: a new strategy for voxel-wise evaluation.

Abstract: Implementation of PET/CT in diagnosis of primary prostate cancer (PCa) requires a profound knowledge about the tracer, preferably from a quantitative evaluation. Direct visual comparison of PET/CT slices to whole prostate sections is hampered by considerable uncertainties from imperfect coregistration and fundamentally different image modalities. In the current study, we present a novel method for advanced voxel-wise comparison of histopathology from excised prostates to pre-surgical PET. Resected prostates from eight patients who underwent PSMA-PET/CT were scanned (ex vivo CT) and thoroughly pathologically prepared. In vivo and ex vivo CT including histopathology were coregistered with three different methods (manual, semi−/automatic). Spatial overlap after CT-based registration was evaluated with dice similarity (DSC). Furthermore, we constructed 3D cancer distribution models from histopathologic information in various slices. Subsequent smoothing reflected the intrinsically limited spatial resolution of PSMA-PET. The resulting histoPET models were used for quantitative analysis of spatial histopathology-PET pattern agreement focusing on p values and coefficients of determination (R 2). We examined additional rigid mutual information (MI) coregistration directly based on PSMA-PET and histoPET. Mean DSC for the three different methods (ManReg, ScalFactReg, and DefReg) were 0.79 ± 0.06, 0.82 ± 0.04, and 0.90 ± 0.02, respectively, while quantification of PET-histopathology pattern agreement after CT-based registration revealed R 2 45.7, 43.2, and 41.3% on average with p < 10−5. Subsequent PET-based MI coregistration yielded R 2 61.3, 55.9, and 55.6%, respectively, while implying anatomically plausible transformations. Creating 3D histoPET models based on thorough histopathological preparation allowed sophisticated quantitative analyses showing highly significant correlations between histopathology and (PSMA-)PET. We recommend manual CT-based coregistration followed by a PET-based MI algorithm to overcome limitations of purely CT-based coregistrations for meaningful voxel-wise comparisons between PET and histopathology.

Predictive SIRT dosimetry based on a territorial model

Abstract: In the planning of selective internal radiation therapy (SIRT) for liver cancer treatment, one major aspect is to determine the prescribed activity and to estimate the resulting absorbed dose inside normal liver and tumor tissue. An optimized partition model for SIRT dosimetry based on arterial liver territories is proposed. This model is dedicated to characterize the variability of dose within the whole liver. For an arbitrary partition, the generalized absorbed dose is derived from the classical partition model. This enables to consider normal liver partitions for each arterial perfusion supply area and one partition for each tumor for activity and dose calculation. The proposed method excludes a margin of 11 mm emitting range around tumor volumes from normal liver to investigate the impact on activity calculation. Activity and dose calculation was performed for five patients using the body-surface-area (BSA) method, the classical and territorial partition model. The territorial model reaches smaller normal liver doses and significant higher tumor doses compared to the classical partition model. The exclusion of a small region around tumors has a significant impact on mean liver dose. Determined tumor activities for the proposed method are higher in all patients when limited by normal liver dose. Activity calculation based on BSA achieves in all cases the lowest amount. The territorial model provides a more local and patient-individual dose distribution in normal liver taking into account arterial supply areas. This proposed arterial liver territory-based partition model may be used for SPECT-independent activity calculation and dose prediction under the condition of an artery-based simulation for particle distribution.

An experimental phantom study of the effect of gadolinium-based MR contrast agents on PET attenuation coefficients and PET quantification in PET-MR imaging: application to cardiac studies

Abstract: Simultaneous cardiac perfusion studies are an increasing trend in PET-MR imaging. During dynamic PET imaging, the introduction of gadolinium-based MR contrast agents (GBCA) at high concentrations during a dual injection of GBCA and PET radiotracer may cause increased attenuation effects of the PET signal, and thus errors in quantification of PET images. We thus aimed to calculate the change in linear attenuation coefficient (LAC) of a mixture of PET radiotracer and increasing concentrations of GBCA in solution and furthermore, to investigate if this change in LAC produced a measurable effect on the image-based PET activity concentration when attenuation corrected by three different AC strategies. We performed simultaneous PET-MR imaging of a phantom in a static scenario using a fixed activity of 40 MBq [18 F]-NaF, water, and an increasing GBCA concentration from 0 to 66 mM (based on an assumed maximum possible concentration of GBCA in the left ventricle in a clinical study). This simulated a range of clinical concentrations of GBCA. We investigated two methods to calculate the LAC of the solution mixture at 511 keV: (1) a mathematical mixture rule and (2) CT imaging of each concentration step and subsequent conversion to LAC at 511 keV. This comparison showed that the ranges of LAC produced by both methods are equivalent with an increase in LAC of the mixed solution of approximately 2% over the range of 0–66 mM. We then employed three different attenuation correction methods to the PET data: (1) each PET scan at a specific millimolar concentration of GBCA corrected by its corresponding CT scan, (2) each PET scan corrected by a CT scan with no GBCA present (i.e., at 0 mM GBCA), and (3) a manually generated attenuation map, whereby all CT voxels in the phantom at 0 mM were replaced by LAC = 0.1 cm−1. All attenuation correction methods (1–3) were accurate to the true measured activity concentration within 5%, and there were no trends in image-based activity concentrations upon increasing the GBCA concentration of the solution. The presence of high GBCA concentration (representing a worst-case scenario in dynamic cardiac studies) in solution with PET radiotracer produces a minimal effect on attenuation-corrected PET quantification.

Comparison of novel multi-level Otsu (MO-PET) and conventional PET segmentation methods for measuring FDG metabolic tumor volume in patients with soft tissue sarcoma

Abstract: We have previously developed a novel and highly consistent PET segmentation algorithm using a multi-level Otsu method (MO-PET). The aim of this study was to evaluate the reliability of MO-PET compared to conventional PET segmentation methods for measuring 18F-FDG (FDG) PET metabolic tumor volume (MTV) in patients with soft tissue sarcoma (STS). Clinical and imaging data were obtained from the Cancer Imaging Archive. Forty-eight STS patients with FDG PET/CT and MR prior to therapy were analyzed. MTV of the tumor using MO-PET was compared to other conventional methods (absolute SUV threshold values of 2.0, 2.5, or 3.0 and percentage of tumor SUVmax values of 30, 40, 50, or 60%) and gradient-based method (PET Edge™). The reference volume was defined as an MR-based gross tumor volume (GTV). Spearman, intra-class correlation, and Bland-Altman analysis were performed to evaluate the correlation and agreement of MTV to GTV. MTVs obtained using each conventional SUV parameter, PET Edge™, and MO-PET were highly correlated with the GTV in Spearman and intra-class correlation analysis (p < 0.05). MO-PET and PET Edge™ showed high intra-class correlation coefficient of MTV to GTV (0.93 and 0.84, respectively). The Bland-Altman bias results showed the highest agreement for MTV using MO-PET with GTV (26.0 ± 489.6 cm3) compared to other methods (SUV 2.0 with − 69.3 ± 765.8, 30% SUVmax with − 255.0 ± 876.6, and PET Edge™ with − 26.46 ± 668.82 cm3). PET MTV segmented with MO-PET showed higher correlation and agreement with GTV in comparison to conventional percentage SUVmax and absolute SUV threshold-based PET segmentation methods. MO-PET is comparable to PET Edge™. MO-PET is a reliable and consistent method for measuring tumor MTV.

A CZT-based blood counter for quantitative molecular imaging

Abstract: Robust quantitative analysis in positron emission tomography (PET) and in single-photon emission computed tomography (SPECT) typically requires the time-activity curve as an input function for the pharmacokinetic modeling of tracer uptake. For this purpose, a new automated tool for the determination of blood activity as a function of time is presented. The device, compact enough to be used on the patient bed, relies on a peristaltic pump for continuous blood withdrawal at user-defined rates. Gamma detection is based on a 20 × 20 × 15 mm3 cadmium zinc telluride (CZT) detector, read by custom-made electronics and a field-programmable gate array-based signal processing unit. A graphical user interface (GUI) allows users to select parameters and easily perform acquisitions. This paper presents the overall design of the device as well as the results related to the detector performance in terms of stability, sensitivity and energy resolution. Results from a patient study are also reported. The device achieved a sensitivity of 7.1 cps/(kBq/mL) and a minimum detectable activity of 2.5 kBq/ml for 18F. The gamma counter also demonstrated an excellent stability with a deviation in count rates inferior to 0.05% over 6 h. An energy resolution of 8% was achieved at 662 keV. The patient study was conclusive and demonstrated that the compact gamma blood counter developed has the sensitivity and the stability required to conduct quantitative molecular imaging studies in PET and SPECT.

Evaluation of two intraoperative gamma detectors for assessment of 177 Lu activity concentration in vivo.

Abstract: Patients with somatostatin receptor-expressing neuroendocrine tumours can be treated with intravenously administered 177Lu-octreotate. Few patients are cured with the present protocol due to the current dose limitation of normal organs at risk, such as the kidneys. By locally administering 177Lu-octreotate to the liver for the purpose of treating liver metastases, a substantially reduced absorbed dose to organs at risk could be achieved. The development of such a technique requires the capability of measuring the 177Lu activity concentration in tissues in vivo. The aim of this study was to evaluate different performance parameters of two commercially available intraoperative gamma detectors in order to investigate whether intraoperative gamma detector measurements could be used to determine 177Lu activity concentration in vivo. Measurements were made using different sources containing 177Lu. Response linearity, sensitivity, spatial resolution and its depth dependence, organ thickness dependence of the measured count rate and tumour detectability were assessed for two intraoperative gamma detectors. The two detectors (a scintillation and a semiconductor detector) showed differences in technical performance. For example, the sensitivity was higher for the scintillation detector, while the spatial resolution was better for the semiconductor detector. Regarding organ thickness dependence and tumour detectability, similar results were obtained for both detectors, and even relatively small simulated tumours of low tumour-to-background activity concentration ratios could be detected. Acceptable results were obtained for both detectors, although the semiconductor detector proved more advantageous for our purpose. The measurements demonstrated factors that must be corrected for, such as organ thickness or dead-time effects. Altogether, intraoperative gamma detector measurements could be used to determine 177Lu activity concentration in vivo.