SPECT image segmentation for estimation of tumour volume and activity concentration in 177Lu-DOTATATE radionuclide therapy

Background: Dosimetry in radionuclide therapy has the potential to allow for a treatment tailored to the individual patient. One therapeutic radiopharmaceutical where patient-specific dosimetry is feasible is 177Lu-DOTATATE, used for the treatment of neuroendocrine tumours. The emission of gamma photons by 177Lu allows for imaging with SPECT (single photon emission computed tomography). One important step for dosimetry using this imaging technique is the SPECT image segmentation, which needs to be robust and accurate for the estimated quantities to be reliable. This work investigates different methods for automatic tumour delineation in 177Lu-DOTATATE SPECT images. Three segmentation methods are considered: a fixed 42% threshold (FT), the Otsu method (OM) and a method based on Fourier surfaces (FS). Effects of including resolution compensation in the iterative SPECT image reconstruction are also studied. Evaluation is performed based on Monte Carlo-simulated SPECT images from 24 h and 336 h post injection (p.i.), for determination of the volume, activity concentration and dice similarity coefficient. In addition, patient data are used to investigate the correspondence of tumour volumes when delineated in SPECT or morphological CT or MR images. Patient data are also used to examine the sensitivity to the operator-dependent initialization. Results: For simulated images from 24 h p.i. reconstructed without resolution compensation, a volume and activity-concentration root-mean-square error below 15% is typically obtained for tumours above approximately 10 cm3 when using OM or FS, while FT performs considerably worse. When including resolution compensation, the tumour volume becomes underestimated and the activity concentration overestimated. The FS method appears to be robust to noise, as seen for the 336 h images. The differences between the tumour volumes estimated from the SPECT images and the volumes estimated from morphological images are generally larger than the discrepancies seen for the simulated data sets. Conclusions: Segmentation results are encouraging for future dosimetry of tumours with volumes above approximately 10 cm3. Using resolution compensation in the reconstruction may have a negative effect on volume estimation

Biologically effective dose in fractionated molecular radiotherapy-application to treatment of neuroblastoma with (131)I-mIBG.

In this work, the biologically effective dose (BED) is investigated for fractionated molecular radiotherapy (MRT). A formula for the Lea-Catcheside G-factor is derived which takes the possibility of combinations of sub-lethal damage due to radiation from different administrations of activity into account. In contrast to the previous formula, the new G-factor has an explicit dependence on the time interval between administrations. The BED of tumour and liver is analysed in MRT of neuroblastoma with (131)I-mIBG, following a common two-administration protocol with a mass-based activity prescription. A BED analysis is also made for modified schedules, when due to local regulations there is a maximum permitted activity for each administration. Modifications include both the simplistic approach of delivering this maximum permitted activity in each of the two administrations, and also the introduction of additional administrations while maintaining the protocol-prescribed total activity. For the cases studied with additional (i.e. more than two) administrations, BED of tumour and liver decreases at most 12% and 29%, respectively. The decrease in BED of the tumour is however modest compared to the two-administration schedule using the maximum permitted activity, where the decrease compared to the original schedule is 47%

Testis dosimetry in individual patients by combining a small-scale dosimetry model and pharmacokinetic modeling-application of (111)In-Ibritumomab Tiuxetan (Zevalin(®)).

A heterogeneous distribution of radionuclides emitting low-energy electrons in the testicles may result in a significant difference between an absorbed dose to the radiosensitive spermatogonia and the mean absorbed dose to the whole testis. This study focused on absorbed dose distribution in patients at a finer scale than normally available in clinical dosimetry, which was accomplished by combining a small-scale dosimetry model with patient pharmacokinetic data. The activity in the testes was measured and blood sampling was performed for patients that underwent pre-therapy imaging with (111)In-Zevalin(®). Using compartment modeling, testicular activity was separated into two components: vascular and extravascular. The uncertainty of absorbed dose due to geometry variations between testicles was explored by an assumed activity micro-distribution and by varying the radius of the interstitial tubule. Results showed that the absorbed dose to germ cells might be strongly dependent on the location of the radioactive source, and may exceed the absorbed dose to the whole testis by as much as a factor of two. Small-scale dosimetry combined with compartmental analysis of clinical data proved useful for gauging tissue dosimetry and interpreting how intrinsic geometric variation influences the absorbed dose

EANM dosimetry committee recommendations for dosimetry of 177Lu-labelled somatostatin-receptor- and PSMA-targeting ligands

The purpose of the EANM Dosimetry Committee is to provide recommendations and guidance to scientists and clinicians on patient-specific dosimetry. Radiopharmaceuticals labelled with lutetium-177 (177Lu) are increasingly used for therapeutic applications, in particular for the treatment of metastatic neuroendocrine tumours using ligands for somatostatin receptors and prostate adenocarcinoma with small-molecule PSMA-targeting ligands. This paper provides an overview of reported dosimetry data for these therapies and summarises current knowledge about radiation-induced side effects on normal tissues and dose-effect relationships for tumours. Dosimetry methods and data are summarised for kidneys, bone marrow, salivary glands, lacrimal glands, pituitary glands, tumours, and the skin in case of radiopharmaceutical extravasation. Where applicable, taking into account the present status of the field and recent evidence in the literature, guidance is provided. The purpose of these recommendations is to encourage the practice of patient-specific dosimetry in therapy with 177Lu-labelled compounds. The proposed methods should be within the scope of centres offering therapy with 177Lu-labelled ligands for somatostatin receptors or small-molecule PSMA.</p

Inter-comparison of quantitative imaging of lutetium-177 (177Lu) in European hospitals

Background: This inter-comparison exercise was performed to demonstrate the variability of quantitative SPECT/CT imaging for lutetium-177 (177Lu) in current clinical practice. Our aim was to assess the feasibility of using international inter-comparison exercises as a means to ensure consistency between clinical sites whilst enabling the sites to use their own choice of quantitative imaging protocols, specific to their systems. Dual-compartment concentric spherical sources of accurately known activity concentrations were prepared and sent to seven European clinical sites. The site staff were not aware of the true volumes or activity within the sources—they performed SPECT/CT imaging of the source, positioned within a water-filled phantom, using their own choice of parameters and reported their estimate of the activities within the source. Results: The volumes reported by the participants for the inner section of the source were all within 29% of the true value and within 60% of the true value for the outer section. The activities reported by the participants for the inner section of the source were all within 20% of the true value, whilst those reported for the outer section were up to 83% different to the true value. Conclusions: A variety of calibration and segmentation methods were used by the participants for this exercise which demonstrated the variability of quantitative imaging across clinical sites. This paper presents a method to assess consistency between sites using different calibration and segmentation methods

EANM practical guidance on uncertainty analysis for molecular radiotherapy absorbed dose calculations

A framework is proposed for modelling the uncertainty in the measurement processes constituting the dosimetry chain that are involved in internal absorbed dose calculations. The starting point is the basic model for absorbed dose in a site of interest as the product of the cumulated activity and a dose factor. In turn, the cumulated activity is given by the area under a time–activity curve derived from a time sequence of activity values. Each activity value is obtained in terms of a count rate, a calibration factor and a recovery coefficient (a correction for partial volume effects). The method to determine the recovery coefficient and the dose factor, both of which are dependent on the size of the volume of interest (VOI), are described. Consideration is given to propagating estimates of the quantities concerned and their associated uncertainties through the dosimetry chain to obtain an estimate of mean absorbed dose in the VOI and its associated uncertainty. This approach is demonstrated in a clinical example

Whole-remnant and maximum-voxel SPECT/CT dosimetry in ¹³¹I-NaI treatments of differentiated thyroid cancer.

Purpose To investigate the possible differences between SPECT/CT based whole-remnant and maximum-voxel dosimetry in patients receiving radio-iodine ablation treatment of differentiated thyroid cancer (DTC).Methods Eighteen DTC patients were administered 1.11 GBq of 131I-NaI after near-total thyroidectomy and rhTSH stimulation. Two patients had two remnants, so in total dosimetry was performed for 20 sites. Three SPECT/CT scans were performed for each patient at 1, 2, and 3-7 days after administration. The activity, the remnant mass, and the maximum-voxel activity were determined from these images and from a recovery-coefficient curve derived from experimental phantom measurements. The cumulated activity was estimated using trapezoidal-exponential integration. Finally, the absorbed dose was calculated using S-values for unit-density spheres in whole-remnant dosimetry and S-values for voxels in maximum-voxel dosimetry.Results The mean absorbed dose obtained from whole-remnant dosimetry was 40 Gy (range 2-176 Gy) and from maximum-voxel dosimetry 34 Gy (range 2-145 Gy). For any given patient, the activity concentrations for each of the three time-points were approximately the same for the two methods. The effective half-lives varied (R = 0.865), mainly due to discrepancies in estimation of the longer effective half-lives. On average, absorbed doses obtained from whole-remnant dosimetry were 1.2 ± 0.2 (1 SD) higher than for maximum-voxel dosimetry, mainly due to differences in the S-values. The method-related differences were however small in comparison to the wide range of absorbed doses obtained in patients.Conclusions Simple and consistent procedures for SPECT/CT based whole-volume and maximum-voxel dosimetry have been described, both based on experimentally determined recovery coefficients. Generally the results from the two approaches are consistent, although there is a small, systematic difference in the absorbed dose due to differences in the S-values, and some variability due to differences in the estimated effective half-lives, especially when the effective half-life is long. Irrespective of the method used, the patient absorbed doses obtained span over two orders of magnitude

Image Processing for Quantitative Scintillation-Camera Imaging. Application to Radionuclide Therapy.

Individual-based determinations of the absorbed dose in radionuclide therapy largely rely on absolute measurement of the activity distribution and its redistribution over time. Scintillation-camera imaging is the most commonly employed measuring technique, applied in planar or SPECT mode, sometimes in combination with structural images from CT. In this thesis, methods for processing and analysis of scintillation-camera images have been studied, with the intent to compensate for the physical limitations involved in the imaging process and thereby improve the quantification accuracy. In particular, methods for image registration and segmentation have been developed. Evaluation of the earlier SPECT based quantification program revealed limitations for the measurement of small volumes, their activity and contrast, in situations where the object contrast deviated from 100%. A 3D deformable segmentation method was developed and applied for SPECT images. In preliminary evaluations, this method gave an adequate boundary, also in situations with low contrast and high levels of noise. For registration of SPECT and CT images, a method was developed where the possibility to use images acquired in the Compton scatter region was explored. The registration accuracy was obtained to about 10 mm, where the largest deviations occurred for slices with heterogeneous activity distributions. To achieve automation, this method was further developed by applying mutual information as similarity measure and spatial transformations in 3D. This registration method was applied for activity quantification and absorbed dose calculations in 3D, using the combined information from CT and SPECT. For planar image activity quantification, a 2D method for registration of whole-body emission and transmission images was developed, including a tailored spatial transformation. The accuracy was obtained to below pixel level (< 3.6 mm) for simulated images, to 9 mm in comparison to point markers for patients. A quantification method based on conjugate scintillation-camera images was also developed, in which a registered CT image was used for attenuation correction, background compensation and segmentation of organs. Evaluation using Monte Carlo simulated images showed an accuracy of within 10% for organ activity quantification. Throughout this thesis, the major application has been images acquired in association with radionuclide therapy using an 131I labeled monoclonal antibody