# Clinical implementation of PLANET® Dose for dosimetric assessment after [177Lu]Lu-DOTA-TATE: comparison with Dosimetry Toolkit® and OLINDA/EXM® V1.0

Lore Santoro, L. Pitalot, D. Trauchessec, E. Mora-Ramirez, P. Kotzki, M.

Bardiès, E. Deshayes

► To cite this version:

Lore Santoro, L. Pitalot, D. Trauchessec, E. Mora-Ramirez, P. Kotzki, et al. Clinical imple-
mentation of PLANET® Dose for dosimetric assessment after [177Lu]Lu-DOTA-TATE: compari-
son with Dosimetry Toolkit® and OLINDA/EXM® V1.0. EJNMMI Research, 2021, 11 (1), pp.1.
10.1186/s13550-020-00737-8. hal-03619967

HAL Id: hal-03619967

https://hal.umontpellier.fr/hal-03619967

Submitted on 25 Mar 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ORIGINAL RESEARCH

Open Access

<figure><img src="image_2.png" /></figure>

Check for
updates

# Clinical implementation of PLANET® Dose for dosimetric assessment after [¹⁷⁷Lu] Lu-DOTA-TATE: comparison with Dosimetry Toolkit® and OLINDA/EXM® V1.0

Lore Santoro<fnref n="1" />*, L. Pitalot<fnref n="1" />, D. Trauchessec<fnref n="1" />, E. Mora-Ramirez²,³,⁴, P. O. Kotzki<fnref n="1" />,⁵, M. Bardiès²,³ and E. Deshayes<fnref n="1" />,⁵

**Abstract**

**Background:** The aim of this study was to compare a commercial dosimetry workstation (PLANET® Dose) and the dosimetry approach (GE Dosimetry Toolkit® and OLINDA/EXM® V1.0) currently used in our department for quantification of the absorbed dose (AD) to organs at risk after peptide receptor radionuclide therapy with [¹⁷⁷Lu]Lu-DOTA-TATE.

**Methods:** An evaluation on phantom was performed to determine the SPECT calibration factor variations over time and to compare the Time Integrated Activity Coefficients (TIACs) obtained with the two approaches. Then, dosimetry was carried out with the two tools in 21 patients with neuroendocrine tumours after the first and second injection of 7.2 ± 0.2 GBq of [¹⁷⁷Lu]Lu-DOTA-TATE (40 dosimetry analyses with each software). SPECT/CT images were acquired at 4 h, 24 h, 72 h and 192 h post-injection and were reconstructed using the Xeleris software (General Electric). The liver, spleen and kidneys masses and TIACs were determined using Dosimetry Toolkit® (DTK) and PLANET® Dose. The ADs were calculated using OLINDA/EXM® V1.0 and the Local Deposition Method (LDM) or Dose voxel-Kernel convolution (DK) on PLANET® Dose.

**Results:** With the phantom, the 3D calibration factors showed a slight variation (0.8% and 3.3%) over time, and TIACs of 225.19 h and 217.52 h were obtained with DTK and PLANET® Dose, respectively. In patients, the root mean square deviation value was 8.9% for the organ masses, 8.1% for the TIACs, and 9.1% and 7.8% for the ADs calculated with LDM and DK, respectively. The Lin's concordance correlation coefficient was 0.99 and the Bland-Altman plot analysis estimated that the AD value difference between methods ranged from – 0.75 to 0.49 Gy, from – 0.20 to 0.64 Gy, and from – 0.43 to 1.03 Gy for 95% of the 40 liver, kidneys and spleen dosimetry analyses. The dosimetry method had a minor influence on AD differences compared with the image registration and organ segmentation steps.

**Conclusions:** The ADs to organs at risk obtained with the new workstation PLANET® Dose are concordant with those calculated with the currently used software and in agreement with the literature. These results validate the use of PLANET® Dose in clinical routine for patient dosimetry after targeted radiotherapy with [¹⁷⁷Lu]Lu-DOTA-TATE.

**Keywords:** Dosimetry workstation, Peptide receptor radionuclide therapy, [¹⁷⁷Lu]Lu-DOTA-TATE, 3D calibration factor, MIRD, Voxel-based dosimetry

**Background**

Dosimetry applications are expanding and the number of nuclear medicine departments performing patient dosimetry is growing, especially for patients with neuroendocrine tumours receiving peptide receptor radionuclide therapy [1]. This is possible thanks to the

<figure><img src="image_3.png" /></figure>

<footnote>*Correspondence: lore.santoro@icm.unicancer.fr</footnote>

<footnote n="1">¹ Nuclear Medicine Department, Montpellier Cancer Institute (ICM), Univ. Montpellier, 208 Avenue des Apothicaires, 34298 Montpellier Cedex 5, France</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
multidisciplinary collaboration between nuclear medicine physicians, medical physicists, and nuclear medicine technologists [2].

In the clinic, many medical teams have developed their own methodology using the tools available in their department and according to their own organizational possibilities [3–5]. This has led to the local implementation of in-house developed dosimetry software and programs [6–10]. However, the legislation on Medical Devices restricts their use to clinical trials. For routine clinical use, a software package must have received the CE mark. Therefore, many medical teams are now acquiring commercial packages, because obtaining the CE label is usually beyond the missions/capabilities of academic structures and due to the increasing availability of commercial software tools. In a previous article we evaluated some commercial packages already on the market [11]. An updated table (Table 1) is presented below, but the field is rapidly and constantly evolving. Only CE-marked software tools are presented (therefore, STRATOS from the Philips research station Imalytics is not included), and some features that are still under development may not have been approved yet.

In our department, dosimetry is integrated in the clinical routine for patients with neuroendocrine tumours treated by peptide receptor radionuclide therapy with [¹⁷⁷Lu]Lu-DOTA-TATE. Up to now, dosimetry analyses of organs at risk (OAR) have been performed using the combination of Dosimetry Toolkit® (GE Healthcare,

Milwaukee, USA) and OLINDA/EXM® V1.0 to calculate the Time Integrated Activity Coefficient (TIAC) and organ-level absorbed dose (AD), respectively [12]. Recently, our department acquired PLANET® Dose (DOSIsoft SA, Cachan, France), a new CE-marked commercial dosimetry workstation. The initial motivation for changing was to use a vendor-neutral solution to support multicentric trials and to put in place a central dosimetry system. However, before its clinical routine implementation, we wanted to compare this new dosimetric package with our internal reference. Our validation plan involved comparing the results obtained on phantom (calibration factor, TIAC) and in patients (mean AD to OARs, TIACs and organ volumes), using similar parameters in terms of segmentation, registration, and time activity curve fitting. This allowed assessing the consistency between the dosimetric results obtained using Dosimetry Toolkit® and OLINDA/EXM® V1.0 (our reference) and with PLANET® Dose.

## Material and methods

### Dosimetry software packages

The characteristics of the two dosimetry software packages used in this study are presented in Table 2.

Dosimetry Toolkit® is an application of the Xeleris® software (GE Healthcare, Milwaukee, USA) [13, 15]. GE Healthcare recommends to use a procedure based on planar acquisition to determine the calibration factor (CF; in counts.s⁻¹.MBq⁻¹). For clinical dosimetry, different

<table><caption>**Table 1** List of commercial packages with CE marking</caption><thead><tr><td>Name</td><td>Manufacturer</td><td>Image format</td><td>VOI/voxel dosimetry</td><td>Absorbed dose calculation</td><td>References</td></tr></thead><tbody><tr><td>Dosimetry Toolkit®</td><td>GE Healthcare, Waukesha, WI, USA</td><td>SPECT/CT Hybrid Planar</td><td>VOI (volume of interest)</td><td>OLINDA/EXM® V1/V2</td><td>https://www.gehealthcare.com/products/molecular-imaging/nuclear-medicine/xeleris-4-dr<br>Stabin et al. [14]<br>Kupitz et al. [15]</td></tr><tr><td>PLANET® Dose</td><td>DOSIsoft, Cachan, France</td><td>SPECT/CT Hybrid</td><td>Voxel-based absorbed dose rates calculation with integration over the VOI</td><td>Local energy deposition<br>Convolution<br>+/- density correction</td><td>https://www.dosisoft.com/products/planet-dose/Huizing et al. [16]</td></tr><tr><td>MIM SurePlan™ MRT</td><td>MIM Software Inc., Cleveland, OH, USA</td><td>SPECT/CT Hybrid</td><td>Voxel-based</td><td>Convolution</td><td>https://www.mimsoftwar.com/mim_sureplan_mrt<br>Maughan et al. [17]</td></tr><tr><td>Organ Dosimetry™<br>Voxel Dosimetry™</td><td>Hermes Medical Solutions, Stockholm, Sweden</td><td>SPECT/CT Hybrid Planar</td><td>VOI/Voxel-based</td><td>OLINDA/EXM® V2<br>Monte Carlo method</td><td>https://www.hermesmedical.com/dosimetry/Hippeläinen et al. [18]</td></tr><tr><td>QDOSE®</td><td>ABX-CRO, Dresden, Germany</td><td>SPECT/CT Hybrid Planar</td><td>VOI</td><td>IDAC-Dose 2.1<br>Convolution</td><td>https://www.quantitati vedose.com<br>Barna et al. [19]</td></tr></tbody></table>

GE, general electric; VOI, volume of interest; SPECT/CT, single photon emission computed tomography/computed tomography; IDAC, internal dose assessment by computer
<table><caption>Table 2 General characteristics for each step of the dosimetry workflow in the currently used dosimetry “Dosimetry Toolkit + OLINDA” software (the reference in this study), and the new PLANET® Dose package</caption><thead><tr><td></td><td>Dosimetry Toolkit® + OLINDA/EXM® V1.0</td><td>PLANET® Dose</td></tr></thead><tbody><tr><td>Calibration procedure</td><td>GE recommendations: planar acquisitions<br>Calibration factor: counts.s<sup>-1</sup>.MBq<sup>-1</sup></td><td>No manufacturer's recommendations<br>Calibration factor: Bq.Count<sup>-1</sup></td></tr><tr><td>Clinical imaging expected</td><td>Whole body scans, SPECT/CT or hybrid<br>5 time points maximum<br>Only from GE devices</td><td>SPECT/CT or hybrid<br>Unlimited time points<br>Import in DICOM format from all devices</td></tr><tr><td>Reconstruction</td><td>« Preparation for Dosimetry Toolkit» application</td><td>Non applicable</td></tr><tr><td>Registration</td><td>« Preparation for Dosimetry Toolkit» application<br>Automatic rigid registration between CT images (Full FOV)</td><td>Rigid or elastic registration<br>Using CT data</td></tr><tr><td>Segmentation and propagation</td><td>Manual, semi-automatic or automatic segmentation using the first images (NM or CT)<br>Rigid propagation (constant volumes)<br>Segmentation adjusted by translation or rotation on others CT images</td><td>Manual, semi-automatic or automatic segmentation using the first images (NM or CT)<br>Rigid propagation (constant volumes)<br>No adjustment on others CT images</td></tr><tr><td>TIA fitting</td><td>Mono-exponential exclusively</td><td>Mono-exponential, Bi-exponential, Xexp, Trapezoidal, Tri-exponential, ...</td></tr><tr><td>Absorbed dose calculation</td><td>Not available on Dosimetry Toolkit®<br>TIAC exported to OLINDA/EXM® V1.0</td><td>Local Deposition Method (LDM)<br>Dose Kernel (DK) approach with / without density correction</td></tr><tr><td>Import/export</td><td>Not applicable</td><td>DICOM format (DICOM-RT-Structure and RT-Dose)</td></tr></tbody></table>

GE, general electric; FOV, field of view; NM, nuclear medicine; CT, computed tomography; TIA, Time integrated activity; TIAC, time integrated activity coefficient; DICOM-RT, digital imaging and communications in medicine-radiotherapy

scenarios are available: whole-body, hybrid, or multi-SPECT/CT image acquisition. It includes two steps. The first, “Preparation for Dosimetry Toolkit”, is used for the reconstruction of SPECT/CT raw data and registration of the CT or planar whole body images. The second, “Dosimetry Toolkit”, is used to segment the different organs, create the time activity curves fitted by a mono-exponential function, and calculate the TIAC for each of them. These TIACs are then uploaded in OLINDA/EXM® V1.0 [14] to calculate the organ mass-adjusted ADs. At the time of the study, only OLINDA V1 was available. The future version of Dosimetry Toolkit should propose OLINDA V2, but any model-based dosimetry software that includes S values (IDAC-Dose for example) could be used [20].

PLANET® Dose is a treatment planning system from DOSIsoft. The calibration procedure is left to the user's discretion and a CF (in Bq.count<sup>-1</sup>) is required. This dosimetry package does not reconstruct SPECT/CT data, but accepts reconstructed data in DICOM format from all devices. It provides multi-time point registration (using rigid and elastic algorithms), organ segmentation (manual and automatic), and TIAC calculation with a wide choice of interpolation methods of the time-activity curve (linear, trapezoidal, mono-, X-, bi-, tri-exponential...). The mean AD can be calculated by assuming the local energy deposition (local deposition model—LDM) or by convolution of AD voxel kernels (dose kernels—DK), with or without media density correction [21–23].

For this study, the segmentation, registration and TIAC steps were carried out on PLANET® Dose using parameters similar to those defined in Dosimetry Toolkit® to allow the comparison.

## Dosimetry imaging protocol

All imaging acquisitions were performed with a SPECT/CT Discovery NM/CT 670 (General Electric [GE] Healthcare), including a BrightSpeed 16 CT scanner and a 3/8-inch NaI(Tl) crystal, according to the acquisition protocol previously described [12]. Briefly, nuclear medicine images were acquired using a medium-energy general purpose parallel-hole collimator. A 20% energy window centred on the 208 keV photopeak and a 10% scatter correction window centred on 177 keV were applied. NM acquisitions were performed using a body contour option, rotation of 180° per detector, total of 60 projections and 45 s each. For attenuation correction, CT images were acquired (120 kV, automatic mA regulation with a max at 200 mA, noise index at 6.43, slice thickness of 5 mm, rotation time of 0.8 s, pitch 1.375, 512 × 512 pixels matrix), with standard reconstruction.

The application “Preparation for Dosimetry Toolkit” was used for SPECT/CT image reconstruction for both dosimetry approaches. The Ordered Subset Expectation Maximization iterative reconstruction algorithm was used with 6 iterations and 10 subsets, attenuation, scatter, recovery resolution corrections and a Gaussian post-filter of 0.11 cm [12].
Phantom study

Calibration factor and time integrated activity coefficient

A NEMA IEC body phantom (Body Phantom NU2-2001/2007) that contains two bottles of 250 mL filled with 200 mL of 82.2±4.1 (i.e. maximum activity measurement error of 5%) MBq [¹⁷⁷Lu]Lu-DOTA-TATE was chosen with the aim of mimicking the size of kidneys. The background was filled with non-radioactive water. SPECT/CT images were acquired at different time points to evaluate CF variations over time. The CF was estimated using one of the two bottles filled with [¹⁷⁷Lu]Lu-DOTA-TATE (Fig. 1).

Dosimetry toolkit A CT rigid registration based on
the full phantom was performed with “Preparation for
Dosimetry Toolkit”. Then, using the “Dosimetry Toolkit”

application, an isocontour representing a volume of
200 mL was automatically segmented on the first nuclear
medicine image and was replicated for the images at 24 h,
72 h, 120 h and 216 h. For each time-point, the segmented
volume was kept constant. To determine the CF in counts.
s⁻¹.MBq⁻¹ at each time point, the number of events in the
volume was divided by the acquisition time provided by
the DICOM data (2700 s) and by the activity. For radio-
activity decay correction, a physical half-life of 6.647 days
[24] was applied and the activity at each acquisition time
point was corrected for the phantom preparation time.
The CF was calculated as the mean of the CF values
obtained at the different time points. To obtain the time-
activity curve fitted by a mono-exponential function and
the TIAC (h), information about the radionuclide and the
CF were entered in the appropriate interface.

<figure><img src="image_4.png" /></figure>
**PLANET® Dose** For the body phantom, CT image rigid registration based on the volume that tightly enclosed the bottle was performed with PLANET® Dose. Similarly, an isocontour that represented a volume of 200 mL was automatically segmented on the first functional image and strictly transferred to the SPECT/CT images acquired at the other time points using a rigid contour propagation algorithm. The segmented volume was maintained over time. To determine the CF in Bq.count⁻¹ at each time point, the number of events in the segmented volume was divided by the activity in Bq, by taking into account the radioactivity decay. The CF was calculated as the mean of the CF values obtained at the different time points.

A mono-exponential fitting function, similar to the Dosimetry Toolkit® approach, was used to calculate the TIA Ã, and the TIAC τ, as [25]:

$$ \tilde{A} = \left[ \frac{A_0 \times T_e}{\ln 2} \right] \quad \text{and} \quad \tau = \left[ \frac{1}{\ln 2} \right] \times T_e = 230.15 \text{ h} $$

where A₀ is the initial activity, and T_e is the experimental half-life assessed using a mono-exponential fit.

As in that situation, only physical decay is observed (T=6.647 days) [24], it is possible to derive the reference TIAC τ ref as:

$$ \tau_{\text{ref}} = \left[ \frac{1}{\ln 2} \right] \times T = 230.15 \text{ h} $$

The reference TIAC τ ref was compared to the experimental τ values obtained with the two dosimetry platforms.

## Clinical study

### Patients and treatment

Twenty-one patients (5 women and 16 men; median age 68 years, range 41–82 years) with a neuroendocrine tumour and treated with [¹⁷⁷Lu]Lu-DOTA-TATE, Lutathera® (Advanced Accelerator Applications, Saint Genis Pouilly, France) were evaluated (Additional file 1: Table S1). The treatment consisted in 7.2±0.2 GBq activity (four infusions in total) injected every 8 weeks. Amino acids (lysine + arginine) were administered concomitantly to ensure renal protection by reducing tubular reabsorption of the radiolabelled peptides. All patients were hospitalized in specialized radioprotection rooms for 24 h after injection. Patients were then discharged and had to come back for the other post-infusion imaging sessions. Dosimetry to liver, kidneys and spleen (i.e. the OARs) was performed after the first and second infusion of [¹⁷⁷Lu]Lu-DOTA-TATE. Dosimetry data after the first injection were not evaluable in one patient, and another patient died before the second infusion. In total,

40 dosimetry analyses were performed with each dosimetry package.

### Dosimetry workflow

The dosimetry workflow for the two packages is presented in Table 3. SPECT/CT images were acquired at 4 h, 24 h, 72 h and 192 h after infusion. For some patients, due to health problems, technical issues or calendar reasons, SPECT/CT images were acquired at only three time points after injection. As dosimetry for the first two injections of [¹⁷⁷Lu]Lu-DOTA-TATE is performed routinely in our department, currently with Dosimetry Toolkit + OLINDA, we carried out a retrospective additional analysis of the already available dosimetric data with PLANET® Dose.

**Reference dosimetry method** For both infusions, after the last SPECT/CT image acquisition at 192 h, all SPECT/CT data were loaded on the “Preparation for Dosimetry Toolkit” application. Imaging data were reconstructed and an automatic rigid registration of the CT scans was performed. The results were loaded on the “Dosimetry Toolkit” application. Liver, kidneys and spleen were manually segmented using the CT images collected at 4 h post-injection, and then the segmented contour was rigidly propagated to the 24 h, 72 h, and 192 h images. For each time point, the segmented volume was maintained, but sometimes it was adjusted by translation or rotation. The partial volume effect correction was considered negligible. The administered activity, date and time of administration, radionuclide and CF values (in counts.s⁻¹.MBq⁻¹) were entered. To obtain the TIAC, the time-activity curves were fitted using a mono-exponential function, the only fitting model available in the “Dosimetry Toolkit” application. Then, the TIAC values were exported to OLINDA/EXM® V1.0 to calculate the ADs to liver, kidneys and spleen (i.e. the OARs). The organ masses included in this software were determined from the organ volume defined on the CT images using “Dosimetry Toolkit” and the biological tissue density (1.06 g.cm⁻³ for liver and spleen; 1.05 g.cm⁻³ for kidneys).

**PLANET® Dose** The transversal slices reconstructed using the “Preparation for Dosimetry Toolkit” application and the corresponding CT images were uploaded on PLANET® Dose. Automatic and rigid registration was performed iteratively for each organ, based on a volume tightly enclosing the organ of interest. Then, the organs of interest on the first CT image were manually segmented, and the volumes were transferred from the reference CT image to the SPECT/CT images acquired at the other time points by rigid propagation. The volume of each OAR remained constant at all time points. As mentioned
<table><caption>**Table 3** Dosimetry workflow in patients using the reference dosimetry approach Dosimetry Toolkit + OLINDA and the new workstation PLANET® Dose</caption><thead><tr><th></th><th>Dosimetry Toolkit® + OLINDA/EXM® V1.0</th><th>PLANET® Dose</th></tr></thead><tbody><tr><td>Calibration procedure</td><td>SPECT/CT acquisitions<br>Calibration factor (counts.s<sup>-1</sup>.MBq<sup>-1</sup>)</td><td>SPECT/CT acquisitions<br>Calibration factor (Bq.Count<sup>-1</sup>)</td></tr><tr><td>Clinical imaging</td><td>4 SPECT/CT at 4 h, 24 h, 72 h and 192 h after injection—60 projections of 45 s</td><td></td></tr><tr><td>Reconstruction</td><td>« Preparation for Dosimetry Toolkit» 6 iterations 10 subsets-AC, SC, RR, Gaussian post filter 0.11 cm</td><td></td></tr><tr><td>Registration</td><td>« Preparation for Dosimetry Toolkit» application<br>Automatic rigid registration using CT images (Full FOV)</td><td>Organ-based rigid registration using CT images</td></tr><tr><td>Segmentation and propagation</td><td>Manual segmentation using the first CT image<br>Rigid propagation (constant volumes)<br>Segmentation adjusted by translation or rotation for other images</td><td>Manual segmentation using the first CT image<br>Rigid propagation (constant volumes)<br>No adjustment for other images</td></tr><tr><td>TIA fitting</td><td>Mono-exponential</td><td>Mono-exponential</td></tr><tr><td>Absorbed dose calculation</td><td>OLINDA/EXM® V1.0<br>Patient-adapted organ masses</td><td>Local Deposition Method (LDM)<br>Dose Kernel (DK) approach with/without density correction</td></tr></tbody></table>

AC, attenuation correction; SC, scatter correction; RR, recovery resolution; FOV, field of view; CT, computed tomography; TIA, time integrated activity

above, the partial volume effect correction was considered negligible for the OARs under study.

The administered activity, date and time of administration, radionuclide used, and CF value (in Bq.counts⁻¹) were entered. The time-activity curve was fitted using a similar approach as the one used in Dosimetry Toolkit® (i.e. mono-exponential function) to provide the TIAC (h) and the TIA (MBq.s). The mean ADs were calculated using the LDM and DK methods, with correction of density.

### Mean absorbed dose calculation

The mean AD was calculated using the first two approaches presented in Fig. 2:

* The “Dosimetry Toolkit” approach used TIAC provided by Dosimetry Toolkit®, entered in OLINDA/EXM® V1.0, to derive the mean ADs to liver, kidneys and spleen (Fig. 2a).

* The PLANET® Dose approach used reconstructed images (Preparation for Dosimetry Toolkit) and full processing (registration, segmentation, TIAC and AD calculation using the LDM and DK methods, with density correction) on PLANET® Dose (Fig. 2b).

A third approach was also used to evaluate independently the previous two AD calculation methods. This approach was similar to the PLANET® Dose method, but after the TIAC step, the ADs were computed with OLINDA/EXM® V1.0 (Fig. 2c) using the same phantom model than in the first approach (adult male or female).

### Statistical analyses

Besides the relative difference (in %), the root mean-square deviation (RMSD) of organ masses, TIACs and ADs per cycle obtained with PLANET® Dose and Dosimetry Toolkit + OLINDA (taken as reference) was calculated as follows:

$$RMSD = \sqrt{\frac{\sum_i \left[ \frac{(X_{\text{PLANET dose i}} - X_{\text{DTK OLinda i}})}{X_{\text{DTK OLinda i}}} \right]}{\text{number of dosimetry analysis}}}$$

where $X_{(i)}$ was the organ masses, TIACs or ADs obtained for the dosimetry analysis $i$.

The Lin's concordance correlation coefficient [26] was used to evaluate the agreement between PLANET® Dose LDM and the reference method (Dosimetry Toolkit + OLINDA). This analysis was performed using the mean values for all patients and organs after the two infusions. Moreover, the absolute differences between the mean AD values obtained with the two approaches were assessed for each organ using the Bland-Altman plot analysis [27]. The 95% limits of agreement, from -1.96 to +1.96 SD, were calculated for each organ.

The paired Student's t-test was used to compare the liver, kidneys and spleen ADs calculated with Dosimetry Toolkit + OLINDA and PLANET® Dose LDM ($n=40$). This analysis was performed using the mean values for all patients and organs after the two infusions.

To evaluate independently the AD calculation methods, the mean AD obtained with PLANET® Dose and with “PLANET® Dose + OLINDA/EXM® V1.0” (taken as
<figure><img src="image_5.png" /></figure>

reference) were compared. A bland–Altman analysis and a correlation evaluation were performed.

# Results

## Phantom-based study

SPECT/CT CF values of 5.60±0.04 counts.s⁻¹.MBq⁻¹ and 5.53±0.19 counts.s⁻¹.MBq⁻¹ were obtained with Dosimetry Toolkit® and PLANET® Dose, respectively. These values did not vary significantly over time (0.8% and 3.3% of variation, respectively) (Fig. 3). For PLANET® Dose, the CF value was converted to 67±2.2 Bq.counts⁻¹ to be uploaded on the software.

These CF values were then used in the clinical study to calculate the time activity curves and TIACs for liver, kidneys, and spleen.

TIACs of 217.52 h and 225.19 h were obtained with Dosimetry Toolkit and PLANET® Dose, respectively. These values showed a deviation of –5.5% and –2.2% relative to the reference residence time of 230.15 h.

# Clinical results

## Package comparison

The mean organ masses, TIACs and ADs for liver, kidneys and spleen (i.e. OARs) obtained using each dosimetry package are summarized in Table 4a. The relative differences of the organ mass values and TIACs between Dosimetry Toolkit + OLINDA/EXM® V1.0 (the reference) and PLANET® Dose are presented in

<figure><img src="image_6.png" /></figure>

Table 4b. These results highlighted RMSD values lower than 10% for organ mass values and TIACs, but for the spleen TIAC (RMSD = 10.4%).

The comparison of the OAR AD values obtained with Dosimetry Toolkit + OLINDA and with PLANET® Dose, LDM and DK, are presented in Table 5a–c and Fig. 4. Again, the mean difference and RMSD values were lower than 10% for all organs, except for spleen (RMSD = 10.9%). For kidneys and
spleen, the AD values obtained with PLANET® Dose were slightly, but significantly higher (p < 0.05).

The mean difference between the AD values calculated with LDM and DK (PLANET® Dose) was 2.2%. The values obtained with the LDM method were always higher than those obtained with the DK method (Fig. 4d). Moreover, whatever the software used, liver presented the highest AD variability, with the highest values reaching 16 Gy (Fig. 4a).

### Concordance between packages

For the 40 dosimetry evaluations, the estimated Lin's concordance correlation coefficient was 0.99 (95% CI 0.99; 0.99; R²=0.9736) (Fig. 5a). This result suggests an excellent concordance between our current dosimetry method and PLANET® Dose. According to the Bland-Altman plot method (Fig. 5b-d), the “bias” value (i.e. the mean of the AD differences between PLANET® Dose and Dosimetry Toolkit+OLINDA) was -0.13 Gy for liver, 0.22 Gy for kidneys, and 0.30 Gy for spleen. Moreover, this approach estimated that the difference of AD values between methods ranged from -0.75 to 0.49 Gy, from -0.20 to 0.64 Gy, and from -0.43 to 1.03 Gy for 95% of the 40 liver, kidneys and spleen dosimetry analyses, respectively. For liver, the Bland-Altman analysis showed that the mean absolute dose differences progressively increased with the dose (R²=0.5742).

### Calculation method assessment

The comparison of the mean AD obtained with PLANET® Dose LDM and DK and using the TIACs from PLANET® Dose uploaded on OLINDA/EXM® V1.0 (i.e. the reference in this comparison) (Additional file 2: Fig. S1 and Table 6) showed that the mean relative difference and RMSD values were lower than 5% for all organs, but for spleen (maximum RMSD value of 6.5%).

The estimated Lin's concordance correlation coefficient was 1 (R²=0.9966) (Additional file 3: Fig. S2a). The “bias” value of the Bland-Altman plot analysis (i.e. the average of the differences between PLANET® Dose LDM and PLANET® Dose+OLINDA) was -0.16 Gy for liver, -0.06 Gy for kidneys, and -0.04 Gy for spleen. The difference of AD values between methods ranged from -0.57 to 0.24 Gy, from -0.18 to 0.06 Gy, and from -0.34 to 0.26 Gy for 95% of the 40 liver, kidneys and spleen dosimetry analyses, respectively (Additional file 3: Fig. S2b-d). For liver, the negative trend (Bland-Altman analysis) showed an R² value of 0.8358.

### Discussion

Due to the need of implementing central dosimetric data processing for multicentre trials, our department acquired a vendor-neutral solution, PLANET® Dose, for dosimetry assessments. However, to take advantage of the experience and data already acquired with Dosimetry Toolkit® and OLINDA/EXM® V1.0, we needed to

<table><caption>**Table 4 (a) Mean and standard deviation (n=40 dosimetry analyses) of OAR organ masses, TIACs and ADs calculated using Dosimetry Toolkit+OLINDA/EXM® V1.0 and PLANET® Dose with LDM and DK and density correction. (b) Relative differences (%) of the organ mass values and TIACs between PLANET® Dose and Dosimetry Toolkit+OLINDA, DTK, Dosimetry Toolkit**</caption><caption>(a)</caption><thead><tr><th rowspan="2">ORGANS</th><th colspan="2">Mass (g)</th><th colspan="2">TIAC (h)</th><th colspan="3">Absorbed dose (Gy)</th></tr><tr><th>DTK + OLINDA/EXM® V1.0</th><th>PLANET® Dose</th><th>DTK + OLINDA/EXM® V1.0</th><th>PLANET® Dose</th><th>DTK + OLINDA/EXM® V1.0</th><th>PLANET® Dose</th><th>LDM</th><th>DK</th></tr></thead><tbody><tr><td>Liver</td><td>2141.6 ± 1213.3</td><td>2191.7 ± 1205</td><td>14.9 ± 24.1</td><td>15.1 ± 24.2</td><td>3.40 ± 3.9</td><td>3.27 ± 3.7</td><td>3.21 ± 3.6</td><td></td></tr><tr><td>Kidneys</td><td>478.3 ± 111.4</td><td>461.5 ± 108.3</td><td>2.2 ± 0.7</td><td>2.4 ± 0.8</td><td>3.01 ± 0.9</td><td>3.23 ± 0.9</td><td>3.16 ± 0.9</td><td></td></tr><tr><td>Spleen</td><td>290.4 ± 181.3</td><td>281.7 ± 178.5</td><td>2.0 ± 1.5</td><td>2.0 ± 1.4</td><td>4.15 ± 1</td><td>4.45 ± 1.2</td><td>4.36 ± 1.2</td><td></td></tr></tbody></table><table><caption>(b)</caption><thead><tr><th colspan="10">Relative difference (%)/DTK+OLINDA</th></tr><tr><th rowspan="2"></th><th colspan="3">Mass</th><th colspan="3">TIAC</th><th rowspan="2">Mean</th><th rowspan="2">Min</th><th rowspan="2">Max</th><th rowspan="2">RMSD</th></tr><tr><th>Liver</th><th>Kidneys</th><th>Spleen</th><th>Liver</th><th>Kidneys</th><th>Spleen</th></tr></thead><tbody><tr><td>Mean</td><td>2.9</td><td>-3.4</td><td>-2.4</td><td>5.7</td><td>6.4</td><td>4.9</td><td></td><td></td><td></td><td></td></tr><tr><td>Min</td><td>-4.0</td><td>-18.8</td><td>-23.4</td><td>-6.9</td><td>-9.1</td><td>-18.9</td><td></td><td></td><td></td><td></td></tr><tr><td>Max</td><td>19.7</td><td>16.4</td><td>17.3</td><td>17.3</td><td>25.7</td><td>23.5</td><td></td><td></td><td></td><td></td></tr><tr><td>RMSD</td><td>5.2</td><td>7.2</td><td>8.9</td><td>8.1</td><td>9.4</td><td>10.4</td><td></td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>Table 5 Comparison of the AD values to each OAR: liver (a), kidneys (b) and spleen (c) calculated with Dosimetry Toolkit + OLINDA and PLANET® Dose with density correction. DTK, Dosimetry Toolkit</caption><thead><tr><th rowspan="3">Dosimetry analyses</th><th colspan="3">Absorbed dose to liver (Gy)</th><th colspan="2">Relative difference (%)/DTK + OLINDA</th></tr><tr><th rowspan="2">DTK + OLINDA/EXM® V1.0</th><th rowspan="2">PLANET® Dose LDM</th><th rowspan="2">PLANET® Dose DK</th><th colspan="2">PLANET® Dose LDM</th></tr><tr><th>PLANET® Dose DK</th></tr></thead><tbody><tr><td>1</td><td>1.04</td><td>1.04</td><td>1.02</td><td>-0.1</td><td>-1.9</td></tr><tr><td>2</td><td>1.59</td><td>1.64</td><td>1.61</td><td>3.2</td><td>1.2</td></tr><tr><td>3</td><td>0.92</td><td>0.86</td><td>0.85</td><td>-6.7</td><td>-8.3</td></tr><tr><td>4</td><td>16.44</td><td>15.89</td><td>15.56</td><td>-3.4</td><td>-5.4</td></tr><tr><td>5</td><td>14.19</td><td>13.31</td><td>13.06</td><td>-6.2</td><td>-8.0</td></tr><tr><td>6</td><td>3.27</td><td>3.30</td><td>3.24</td><td>0.9</td><td>-1.0</td></tr><tr><td>7</td><td>4.59</td><td>4.40</td><td>4.32</td><td>-4.1</td><td>-5.9</td></tr><tr><td>8</td><td>0.57</td><td>0.61</td><td>0.60</td><td>7.0</td><td>4.4</td></tr><tr><td>9</td><td>2.05</td><td>2.04</td><td>2.00</td><td>-0.3</td><td>-2.2</td></tr><tr><td>10</td><td>1.26</td><td>1.09</td><td>1.07</td><td>-13.2</td><td>-14.7</td></tr><tr><td>11</td><td>1.12</td><td>1.11</td><td>1.10</td><td>-0.9</td><td>-2.2</td></tr><tr><td>12</td><td>2.40</td><td>2.28</td><td>2.23</td><td>-5.0</td><td>-7.0</td></tr><tr><td>13</td><td>4.02</td><td>3.83</td><td>3.76</td><td>-4.7</td><td>-6.4</td></tr><tr><td>14</td><td>0.91</td><td>0.93</td><td>0.92</td><td>2.1</td><td>0.7</td></tr><tr><td>15</td><td>0.79</td><td>0.80</td><td>0.78</td><td>1.3</td><td>-1.0</td></tr><tr><td>16</td><td>9.58</td><td>9.00</td><td>8.83</td><td>-6.0</td><td>-7.8</td></tr><tr><td>17</td><td>9.32</td><td>8.88</td><td>8.72</td><td>-4.7</td><td>-6.4</td></tr><tr><td>18</td><td>1.47</td><td>1.49</td><td>1.44</td><td>1.5</td><td>-1.9</td></tr><tr><td>19</td><td>1.92</td><td>1.89</td><td>1.86</td><td>-1.8</td><td>-3.3</td></tr><tr><td>20</td><td>0.94</td><td>0.86</td><td>0.84</td><td>-8.7</td><td>-10.7</td></tr><tr><td>21</td><td>1.13</td><td>1.18</td><td>1.16</td><td>4.1</td><td>2.4</td></tr><tr><td>22</td><td>1.65</td><td>1.68</td><td>1.65</td><td>1.8</td><td>0.1</td></tr><tr><td>23</td><td>0.88</td><td>0.93</td><td>0.91</td><td>5.7</td><td>3.5</td></tr><tr><td>24</td><td>6.95</td><td>6.89</td><td>6.75</td><td>-0.9</td><td>-2.8</td></tr><tr><td>25</td><td>2.09</td><td>2.10</td><td>2.06</td><td>0.7</td><td>-1.1</td></tr><tr><td>26</td><td>2.34</td><td>2.29</td><td>2.25</td><td>-2.1</td><td>-3.8</td></tr><tr><td>27</td><td>7.57</td><td>6.82</td><td>6.69</td><td>-9.9</td><td>-11.7</td></tr><tr><td>28</td><td>0.62</td><td>0.65</td><td>0.63</td><td>5.2</td><td>2.7</td></tr><tr><td>29</td><td>1.61</td><td>1.64</td><td>1.61</td><td>1.9</td><td>0.1</td></tr><tr><td>30</td><td>1.14</td><td>1.14</td><td>1.12</td><td>-0.4</td><td>-2.3</td></tr><tr><td>31</td><td>1.39</td><td>1.31</td><td>1.29</td><td>-5.6</td><td>-7.0</td></tr><tr><td>32</td><td>2.24</td><td>2.19</td><td>2.15</td><td>-2.1</td><td>-4.1</td></tr><tr><td>33</td><td>4.13</td><td>4.72</td><td>4.25</td><td>14.3</td><td>2.9</td></tr><tr><td>34</td><td>0.78</td><td>0.84</td><td>0.83</td><td>7.7</td><td>6.5</td></tr><tr><td>35</td><td>0.81</td><td>0.81</td><td>0.80</td><td>0.0</td><td>-1.6</td></tr><tr><td>36</td><td>10.43</td><td>9.19</td><td>9.02</td><td>-11.9</td><td>-13.5</td></tr><tr><td>37</td><td>8.15</td><td>7.51</td><td>7.38</td><td>-7.9</td><td>-9.5</td></tr><tr><td>38</td><td>1.42</td><td>1.36</td><td>1.34</td><td>-4.1</td><td>-5.7</td></tr><tr><td>39</td><td>1.57</td><td>1.49</td><td>1.46</td><td>-5.1</td><td>-6.8</td></tr><tr><td>40</td><td>0.92</td><td>0.94</td><td>0.93</td><td>2.0</td><td>0.5</td></tr><tr><td>Mean</td><td>3.40</td><td>3.27</td><td>3.20</td><td>-1.4</td><td>-3.5</td></tr><tr><td>SD</td><td>3.89</td><td>3.65</td><td>3.58</td><td>5.5</td><td>4.9</td></tr><tr><td>Min</td><td>0.57</td><td>0.61</td><td>0.60</td><td>-13.2</td><td>-14.7</td></tr><tr><td>Max</td><td>16.44</td><td>15.89</td><td>15.56</td><td>14.3</td><td>6.5</td></tr><tr><td>RMSD</td><td></td><td></td><td></td><td>5.6</td><td>6.0</td></tr></tbody></table>
<table><caption>Table 5 (continued)</caption><caption>(b)</caption><thead><tr><th rowspan="2">Dosimetry analyses</th><th colspan="3">Absorbed dose to kidneys</th><th colspan="2">Relative difference (%)/DTK + OLINDA</th></tr><tr><th>DTK + OLINDA/EXM® V1.0</th><th>PLANET® Dose LDM</th><th>PLANET® Dose DK</th><th>PLANET® Dose LDM</th><th>PLANET® Dose DK</th></tr></thead><tbody><tr><td>1</td><td>2.78</td><td>3.14</td><td>3.07</td><td>12.8</td><td>10.3</td></tr><tr><td>2</td><td>2.79</td><td>3.10</td><td>3.03</td><td>11.1</td><td>8.5</td></tr><tr><td>3</td><td>3.26</td><td>3.41</td><td>3.33</td><td>4.6</td><td>2.1</td></tr><tr><td>4</td><td>3.68</td><td>3.45</td><td>3.38</td><td>-6.3</td><td>-8.2</td></tr><tr><td>5</td><td>2.05</td><td>2.03</td><td>1.91</td><td>-0.7</td><td>-6.9</td></tr><tr><td>6</td><td>3.06</td><td>3.66</td><td>3.57</td><td>19.6</td><td>16.8</td></tr><tr><td>7</td><td>3.56</td><td>3.84</td><td>3.75</td><td>7.7</td><td>5.1</td></tr><tr><td>8</td><td>3.48</td><td>3.83</td><td>3.70</td><td>10.0</td><td>6.2</td></tr><tr><td>9</td><td>2.53</td><td>2.40</td><td>2.35</td><td>-5.2</td><td>-7.2</td></tr><tr><td>10</td><td>3.07</td><td>3.29</td><td>3.22</td><td>7.2</td><td>5.0</td></tr><tr><td>11</td><td>2.10</td><td>2.40</td><td>2.34</td><td>14.3</td><td>11.5</td></tr><tr><td>12</td><td>2.51</td><td>2.61</td><td>2.56</td><td>3.9</td><td>1.8</td></tr><tr><td>13</td><td>3.24</td><td>3.35</td><td>3.28</td><td>3.4</td><td>1.1</td></tr><tr><td>14</td><td>2.44</td><td>2.72</td><td>2.66</td><td>11.5</td><td>8.8</td></tr><tr><td>15</td><td>2.20</td><td>2.50</td><td>2.45</td><td>13.6</td><td>11.1</td></tr><tr><td>16</td><td>1.98</td><td>2.10</td><td>2.06</td><td>5.8</td><td>3.8</td></tr><tr><td>17</td><td>2.96</td><td>3.49</td><td>3.41</td><td>18.0</td><td>15.4</td></tr><tr><td>18</td><td>2.90</td><td>3.13</td><td>3.05</td><td>8.1</td><td>5.4</td></tr><tr><td>19</td><td>4.69</td><td>5.37</td><td>5.24</td><td>14.5</td><td>11.7</td></tr><tr><td>20</td><td>2.48</td><td>2.46</td><td>2.40</td><td>-0.8</td><td>-3.2</td></tr><tr><td>21</td><td>3.31</td><td>3.73</td><td>3.65</td><td>12.7</td><td>10.2</td></tr><tr><td>22</td><td>2.85</td><td>3.15</td><td>3.12</td><td>10.7</td><td>9.5</td></tr><tr><td>23</td><td>3.93</td><td>4.00</td><td>3.91</td><td>1.8</td><td>-0.5</td></tr><tr><td>24</td><td>2.60</td><td>2.60</td><td>2.54</td><td>0.0</td><td>-2.3</td></tr><tr><td>25</td><td>3.54</td><td>3.91</td><td>3.82</td><td>10.5</td><td>8.0</td></tr><tr><td>26</td><td>3.31</td><td>3.50</td><td>3.42</td><td>5.7</td><td>3.3</td></tr><tr><td>27</td><td>7.36</td><td>7.09</td><td>6.93</td><td>-3.7</td><td>-5.8</td></tr><tr><td>28</td><td>3.35</td><td>3.86</td><td>3.76</td><td>15.2</td><td>12.3</td></tr><tr><td>29</td><td>2.51</td><td>2.55</td><td>2.49</td><td>1.6</td><td>-0.8</td></tr><tr><td>30</td><td>3.32</td><td>3.34</td><td>3.27</td><td>0.6</td><td>-1.6</td></tr><tr><td>31</td><td>2.62</td><td>2.82</td><td>2.75</td><td>7.5</td><td>5.0</td></tr><tr><td>32</td><td>2.31</td><td>2.55</td><td>2.49</td><td>10.5</td><td>8.0</td></tr><tr><td>33</td><td>3.37</td><td>3.92</td><td>3.83</td><td>16.3</td><td>13.5</td></tr><tr><td>34</td><td>2.53</td><td>2.68</td><td>2.62</td><td>5.9</td><td>3.5</td></tr><tr><td>35</td><td>2.30</td><td>2.71</td><td>2.64</td><td>17.8</td><td>14.9</td></tr><tr><td>36</td><td>2.15</td><td>2.22</td><td>2.17</td><td>3.4</td><td>1.2</td></tr><tr><td>37</td><td>2.08</td><td>2.10</td><td>2.06</td><td>1.2</td><td>-0.8</td></tr><tr><td>38</td><td>2.87</td><td>3.11</td><td>3.04</td><td>8.5</td><td>6.0</td></tr><tr><td>39</td><td>4.57</td><td>5.02</td><td>4.90</td><td>9.8</td><td>7.2</td></tr><tr><td>40</td><td>2.52</td><td>2.80</td><td>2.73</td><td>11.1</td><td>4.6</td></tr><tr><td>Mean</td><td>3.03</td><td>3.25</td><td>3.17</td><td>7.5</td><td>4.9</td></tr><tr><td>SD</td><td>0.95</td><td>0.97</td><td>0.95</td><td>6.5</td><td>6.4</td></tr><tr><td>Min</td><td>1.98</td><td>2.03</td><td>1.91</td><td>-6.3</td><td>-8.2</td></tr><tr><td>Max</td><td>7.36</td><td>7.09</td><td>6.93</td><td>19.6</td><td>16.8</td></tr><tr><td>RMSD</td><td></td><td></td><td></td><td>9.9</td><td>8.0</td></tr></tbody></table>
<table><caption>Table 5 (continued)</caption><caption>(c)</caption><thead><tr><th rowspan="2">Dosimetry analyses</th><th colspan="3">Absorbed dose to spleen</th><th colspan="2">Relative difference (%)/DTK + OLINDA</th></tr><tr><th>DTK + OLINDA/EXM® V1.0</th><th>PLANET® Dose LDM</th><th>PLANET® Dose DK</th><th>PLANET® Dose LDM</th><th>PLANET® Dose DK</th></tr></thead><tbody><tr><td>1</td><td>3.86</td><td>4.12</td><td>4.03</td><td>6.8</td><td>4.5</td></tr><tr><td>2</td><td>4.67</td><td>5.41</td><td>5.28</td><td>15.8</td><td>13.1</td></tr><tr><td>3</td><td>3.83</td><td>4.40</td><td>4.31</td><td>15.0</td><td>12.5</td></tr><tr><td>4</td><td>4.01</td><td>3.68</td><td>3.59</td><td>-8.2</td><td>-10.4</td></tr><tr><td>5</td><td>4.35</td><td>3.93</td><td>3.83</td><td>-9.7</td><td>-12.0</td></tr><tr><td>6</td><td>2.98</td><td>3.30</td><td>3.23</td><td>10.9</td><td>8.3</td></tr><tr><td>7</td><td>5.03</td><td>5.37</td><td>5.28</td><td>6.7</td><td>5.0</td></tr><tr><td>8</td><td>2.07</td><td>2.36</td><td>2.31</td><td>14.3</td><td>11.7</td></tr><tr><td>9</td><td>3.67</td><td>3.40</td><td>3.32</td><td>-7.2</td><td>-9.4</td></tr><tr><td>10</td><td>4.87</td><td>4.59</td><td>4.37</td><td>-5.7</td><td>-10.2</td></tr><tr><td>11</td><td>3.69</td><td>3.83</td><td>3.74</td><td>3.8</td><td>1.4</td></tr><tr><td>12</td><td>4.27</td><td>4.35</td><td>4.25</td><td>1.8</td><td>-0.5</td></tr><tr><td>13</td><td>4.17</td><td>4.47</td><td>4.37</td><td>7.2</td><td>4.7</td></tr><tr><td>14</td><td>3.22</td><td>3.30</td><td>3.22</td><td>2.5</td><td>0.0</td></tr><tr><td>15</td><td>2.91</td><td>3.33</td><td>3.26</td><td>14.6</td><td>12.0</td></tr><tr><td>16</td><td>2.94</td><td>3.11</td><td>3.05</td><td>5.9</td><td>3.8</td></tr><tr><td>17</td><td>2.57</td><td>2.60</td><td>2.54</td><td>1.1</td><td>-1.2</td></tr><tr><td>18</td><td>4.75</td><td>5.19</td><td>5.07</td><td>9.3</td><td>6.8</td></tr><tr><td>19</td><td>5.50</td><td>6.17</td><td>6.02</td><td>12.2</td><td>9.5</td></tr><tr><td>20</td><td>3.77</td><td>3.94</td><td>3.86</td><td>4.6</td><td>2.4</td></tr><tr><td>21</td><td>4.83</td><td>5.64</td><td>5.51</td><td>16.8</td><td>14.2</td></tr><tr><td>22</td><td>3.88</td><td>4.38</td><td>4.35</td><td>12.9</td><td>12.1</td></tr><tr><td>23</td><td>4.91</td><td>5.55</td><td>5.42</td><td>13.0</td><td>10.4</td></tr><tr><td>24</td><td>4.05</td><td>3.81</td><td>3.71</td><td>-5.9</td><td>-8.3</td></tr><tr><td>25</td><td>6.38</td><td>7.30</td><td>7.13</td><td>14.4</td><td>11.8</td></tr><tr><td>26</td><td>4.18</td><td>4.53</td><td>4.42</td><td>8.3</td><td>5.7</td></tr><tr><td>27</td><td>5.57</td><td>5.12</td><td>5.02</td><td>-8.1</td><td>-9.9</td></tr><tr><td>28</td><td>2.50</td><td>2.72</td><td>2.66</td><td>8.8</td><td>6.4</td></tr><tr><td>29</td><td>3.43</td><td>4.20</td><td>4.10</td><td>22.4</td><td>19.3</td></tr><tr><td>30</td><td>4.75</td><td>4.82</td><td>4.70</td><td>1.5</td><td>-1.0</td></tr><tr><td>31</td><td>5.52</td><td>6.51</td><td>5.67</td><td>18.0</td><td>2.8</td></tr><tr><td>32</td><td>4.09</td><td>4.46</td><td>4.36</td><td>9.0</td><td>6.6</td></tr><tr><td>33</td><td>5.02</td><td>5.75</td><td>5.63</td><td>14.5</td><td>12.0</td></tr><tr><td>34</td><td>2.71</td><td>2.80</td><td>2.73</td><td>3.3</td><td>0.6</td></tr><tr><td>35</td><td>3.44</td><td>4.09</td><td>3.98</td><td>18.8</td><td>15.7</td></tr><tr><td>36</td><td>3.81</td><td>4.41</td><td>4.31</td><td>15.8</td><td>13.3</td></tr><tr><td>37</td><td>3.18</td><td>3.25</td><td>3.18</td><td>2.2</td><td>-0.1</td></tr><tr><td>38</td><td>5.86</td><td>6.46</td><td>6.32</td><td>10.2</td><td>7.8</td></tr><tr><td>39</td><td>6.03</td><td>6.57</td><td>6.40</td><td>9.0</td><td>6.1</td></tr><tr><td>40</td><td>4.84</td><td>4.85</td><td>4.75</td><td>0.2</td><td>-1.9</td></tr><tr><td>Mean</td><td>4.15</td><td>4.45</td><td>4.33</td><td>7.2</td><td>4.4</td></tr><tr><td>SD</td><td>1.03</td><td>1.19</td><td>1.14</td><td>8.3</td><td>8.0</td></tr><tr><td>Min</td><td>2.07</td><td>2.36</td><td>2.31</td><td>-9.7</td><td>-12.0</td></tr><tr><td>Max</td><td>6.38</td><td>7.30</td><td>7.13</td><td>22.4</td><td>19.3</td></tr><tr><td>RMSD</td><td></td><td></td><td></td><td>10.9</td><td>9.1</td></tr></tbody></table>
<figure><img src="image_7.png" /><figcaption>**Fig. 4** Box-and-whisker plots showing the ADs to liver (a), kidneys (b) and spleen (c) calculated using Dosimetry Toolkit + OLINDA and the PLANET® Dose packages, as well as the relative differences of AD values between PLANET® Dose and Dosimetry Toolkit + OLINDA (d). DTK, Dosimetry Toolkit</figcaption></figure>

compare their results. Overall, the two systems gave similar dosimetry results in the phantom study (CF and TIAC) and also when using imaging data from patients who received two injections of [¹⁷⁷Lu]Lu-DOTA-TATE (mean organ masses, TIACs and ADs for liver, kidneys and spleen).

In a recent study, Huizing et al. [16] compared dosimetry results obtained with Hybrid Viewer Dosimetry and PLANET® Dose in ten patients using hybrid imaging data obtained from 0.5 h to 72 h post-injection of [¹⁷⁷Lu]Lu-DOTA-TATE. Although our study also concerned PLANET® Dose, our reference method was Dosimetry Toolkit®. Furthermore, our study was based on complete 3D imaging data from 4 to 192 h post-injection in 21 patients undergoing [¹⁷⁷Lu]Lu-DOTA-TATE treatment, thus representing a total of 40 dosimetry analysis.

In our study, the reconstruction step was performed with the GE application “Preparation for Dosimetry Toolkit” because PLANET® Dose does not include this functionality. The latter accepts reconstructed data

supplied by others workstations, unlike the “Dosimetry Toolkit” application.

In a clinical dosimetry study, the preliminary step of calibration is crucial to obtain accurate activity quantification [28, 29]. A detailed calibration protocol is the cornerstone to achieve accurate and reliable image quantification in multicentric studies [30–33]. According to the GE recommendations, CF determined by planar acquisition of a 15 cm diameter Petri dish partially filled with a solution of ¹⁷⁷Lu is sufficient for Dosimetry Toolkit®. However as our first aim was to monitor the AD to OARs, we used a large phantom (i.e. a bottle with a volume relatively close to that of kidney). The same methodology based on SPECT/CT imaging was followed for both calibration and clinical imaging. As described by Gustafsson et al. [34], SPECT image segmentation is essential to determine the activity concentration. For our phantom-based study, we selected a fixed volume threshold method based on an automatically drawn isocontour around the bottle with a volume of 200 mL. This volume was segmented on the SPECT images acquired at the
<figure><img src="image_8.png" /></figure>

first time point and rigidly copied to the others (i.e. rigid propagation). This step is directly affected by the accuracy of the registration of all SPECT/CT images using a rigid algorithm [11, 35, 36]. The obtained CF was similar (within 10%) to the one determined by Peters et al. [31] using a comparable gamma camera model.

Moreover, CFs are expressed in counts.s⁻¹. MBq⁻¹ by Dosimetry Toolkit® and in Bq.count⁻¹ by PLANET® Dose. This implies that the CF must be modified with the acquisition duration for PLANET® Dose. This also highlights a risk of errors because of the lack of a standardized method for introducing CFs in dosimetry software tools. The crucial recommendation at this step is that the conditions used in clinical studies in terms of acquisition and reconstruction parameters must be similar to those used for calibration. We observed a negligible variation of the CF values over time, from T=0–216 h. This means that the same CF can be used for each time point. This observation is particularly interesting for Dosimetry Toolkit® in which a single CF must be

entered, whereas a different CF can be used at each time point with PLANET® Dose.

In the patient study, the liver masses obtained with the two packages were more similar than those obtained for smaller organs, such as kidneys and spleen, and the variability increased with the OAR decreasing size (Table 4b). Actually, as the DICOM-RT-Structure import could not be used in Dosimetry Toolkit®, organ delineation was done manually for each dosimetry package, leading to operator-induced variability. Thus, a small difference between contours could generate a more important relative deviation in smaller than larger organs. Moreover, the results presented in Table 5a–c seem to indicate a volume/mass effect. The RMSD of the AD values increased with the OAR decreasing size. For spleen, this could be explained by our package comparison methodology that used similar parameters in terms of registration and segmentation with constant volumes over time. Indeed, when using Dosimetry Toolkit®, the volume delineated on the first image was maintained, but adjusted by translation or rotation at each successive
<table><caption>Table 6 Mean, SD, range, RMSD value of the AD values to liver, kidneys and spleen calculated with PLANET® Dose + OLINDA and PLANET® Dose with density correction</caption><thead><tr><th colspan="10">Absorbed dose (Gy)</th></tr><tr><th rowspan="2">Liver</th><th colspan="4">Kidneys</th><th colspan="4">Spleen</th></tr><tr><th>PLANET+OLINDA/ PLANET®Dose EXM®V1.0</th><th>LDM</th><th>PLANET®Dose DK</th><th>PLANET+OLINDA</th><th>PLANET+OLINDA/ PLANET®Dose EXM®V1.0</th><th>LDM</th><th>PLANET®Dose DK</th><th>PLANET+OLINDA/ PLANET®Dose EXM®V1.0</th><th>LDM</th><th>PLANET®Dose DK</th></tr></thead><tbody><tr><td>Mean</td><td>3.43</td><td>3.27</td><td>3.20</td><td>-4.5</td><td>-6.2</td><td>3.31</td><td>3.25</td><td>3.17</td><td>-1.9</td><td>-4.2</td><td>4.47</td><td>4.45</td><td>4.33</td><td>-1.1</td><td>-3.4</td></tr><tr><td>SD</td><td>3.84</td><td>3.65</td><td>3.58</td><td>3.6</td><td>3.6</td><td>3.58</td><td>3.58</td><td>3.58</td><td>3.6</td><td>3.6</td><td>3.58</td><td>3.58</td><td>3.58</td><td>3.6</td><td>3.6</td></tr><tr><td>Min</td><td>0.62</td><td>0.61</td><td>0.60</td><td>-6.6</td><td>-8.5</td><td>1.97</td><td>1.94</td><td>1.91</td><td>-5.8</td><td>-7.7</td><td>2.40</td><td>2.36</td><td>2.31</td><td>-5.5</td><td>-7.6</td></tr><tr><td>Max</td><td>16.40</td><td>15.89</td><td>15.56</td><td>0.9</td><td>-1.0</td><td>7.36</td><td>7.09</td><td>6.93</td><td>2.9</td><td>-0.6</td><td>7.58</td><td>7.30</td><td>7.13</td><td>6.5</td><td>4.2</td></tr><tr><td>RMSD</td><td></td><td></td><td></td><td>4.9</td><td>6.5</td><td></td><td></td><td></td><td>2.5</td><td>4.4</td><td>4.6</td><td></td><td></td><td>3.3</td><td>4.6</td></tr></tbody></table>
time point, when necessary. However, this step was not available in PLANET® Dose. Therefore, to use similar approaches with both packages, in PLANET® Dose we chose to perform an organ-based registration followed by rigid propagation (i.e. the contour registered in the first CT image was exactly propagated to the images acquired at the other time points). Thus, for some patients, the initial spleen contour did not fully match the organ contour at the other time points, and included tissues with different density. This implied an important deviation when calculating the AD with density correction. For instance, in the dosimetry analysis n°31, part of the spleen contour was moved to the left lung at other time points, and the AD was overestimated because of the density correction. For kidneys, AD variations between packages could be explained by the fact that in PLANET® Dose, left and right kidney are considered separately, while in OLINDA/EXM, the two kidneys are considered as a single organ. This may affect pharmacokinetic assessments, especially when rigid registration is considered. Therefore, we tested an additional approach that uses the segmentations and TIACs provided by PLANET® Dose and the dosimetric results (ADs) obtained with PLANET® Dose LDM or DK and with OLINDA/EXM® V1.0 (Table 6). The differences between approaches were reduced when the comparison considered only the AD calculation. This is in agreement with the idea that registration and segmentation are major steps, and probably induce more variability than the AD calculation step. Grassi et al. [36] showed that ADs to organs are significantly affected by the registration algorithm used. Indeed, due to respiratory motion during SPECT/CT and organ deformation, registration errors can happen. Although this may not influence the results of our comparison, it is clear that additional studies on registration and time activity curve fitting are needed, by taking advantage of the additional possibilities available in PLANET® Dose. Moreover, for accurate registration, position reproducibility during the four SPECT/CT acquisitions is crucial. Thus, patient set-up and immobilization devices are strongly recommended.

The Bland–Altman plots showed that the biases were quite low for all organs. The limits of agreement were rather tight with a maximum value of approximately 1 Gy for spleen. Moreover, the results for liver highlighted a trend for slightly higher AD values obtained with Dosimetry Toolkit + OLINDA than PLANET® Dose when the AD increased. This trend was more pronounced for patients with liver metastases, i.e. high activity gradients within the liver. One hypothesis is that cross-absorbed doses (from photons) may be different for heterogeneous vs. homogeneous activity distributions. Another point to consider is the density (homogeneous for OLINDA

vs. voxel-based and coming from the CT images for PLANET® Dose). This may influence the AD calculation, but it is not clear to which extent. This certainly deserves to be thoroughly investigated, probably using Monte Carlo modelling to take into account also the possible impact of local density corrections [37].

The concordance evaluation (Lin's coefficient value of 0.99) highlighted an excellent agreement between methods. Moreover, the dosimetry results obtained using PLANET® Dose (AD to liver, kidneys and spleen of 0.45±0.50 Gy/GBq, 0.45±0.13 Gy/GBq and 0.62±0.17 Gy/GBq respectively) are in agreement with those of the literature [38].

As proposed by Gear et al. [39] in a practical guidance paper, the uncertainties at each step of the dosimetry analysis should be determined to express the accuracy of the dosimetry results. Currently, PLANET® Dose allows evaluating the relative proportion of interpolation (between time points) and extrapolation (after the last time point). Regardless of the goodness of fit, this is a good indication of the relevance of time sampling, but it is not sufficient to fully characterize the uncertainty associated with the dosimetric workflow. Such study requires important developments and will be implemented in the future.

From a qualitative point of view, PLANET® Dose is a user-friendly commercial solution that proposes a wide range of tools for segmentation, and several analytic fit functions. The time necessary for a dosimetry analysis is significantly reduced. Therefore, considering the good agreement with our reference dosimetry method, the concordance of the dosimetry results with the literature, the added value of this software (easy contouring, wide choice of time activity curve fitting models, time saving), and the fact that the observed differences were explainable and clinically acceptable, we think that the PLANET® Dose software can replace our current dosimetry package without any correction for dosimetry analysis. We can now start to investigate the different methodological possibilities offered by PLANET® Dose, such as elastic registration and propagation, time integration activity with multiple exponentials, AD rate computation at the voxel level, DICOM-RT import and export of structures. Its potential can now be fully explored, particularly for the determination of the tumour AD and for investigating the AD-response correlations.

## Conclusion

In this work we compared the dosimetric results obtained with the software currently used in our department and with a new dosimetry software package available on the market. The ADs to OARs (liver, kidney and
spleen) obtained with PLANET® Dose were concordant with those calculated with GE Dosimetry Toolkit® and OLINDA/EXM® V1.0, and in agreement with the literature. These results allow us to use PLANET® Dose in clinical routine for patient dosimetry after targeted radiotherapy with [177Lu]Lu-DOTA-TATE.

**Supplementary information**

The online version contains supplementary material available at https://doi.org/10.1186/s13550-020-00737-8.

Additional file 1. Table S1: Patients' characteristics.

Additional file 2. Fig. S1: Box-and-whisker plots showing the ADs to liver (a), kidneys (b) and spleen (c) calculated using PLANET® Dose and PLANET® Dose+OLINDA, as well as the relative AD differences between PLANET® Dose and PLANET® Dose+OLINDA (d).

Additional file 3. Fig. S2: Dispersion around the 45° line of the AD pairs obtained with PLANET+OLINDA and PLANET® Dose LDM with density correction for all organs combined (a). Bland-Altman plots of the AD to liver (b), kidneys (c) and spleen (d) calculated with PLANET+OLINDA and PLANET® Dose LDM with density correction.

**Abbreviations**

AD: Absorbed dose; CF: Calibration factor; DICOM-RT: Digital imaging and communications in medicine-radiotherapy; DK: Dose kernel; DTK: Dosimetry toolkit®; FOV: Field of view; GE: General Electric; IDAC: Internal dose assessment by computer; LDM: Local deposition method; MIRD: Medical internal radiation dose; NETS: Neuroendocrine tumours; NM: Nuclear medicine; OAR: Organ at risk; OSEM: Ordered subset expectation maximization; PRRT: Peptide receptor radionuclide therapy; RMSD: Root mean square deviation; SD: Standard deviation; SPECT/CT: Single photon emission computed tomography/ computed tomography; TACS: Time activity curves; TIA: Time integrated activity; TIAC: Time integrated activity coefficient; VOI: Volume of interest.

**Authors' contributions**

LS, LP, DT, EMR carried out the phantom measurements, image reconstructions and data analysis. MB and ED participated to the design of the study. POK and ED were responsible of the patient acquisitions and treatments. LS was responsible of the dosimetry software implementation and performed patient dosimetry evaluations with LP. LS wrote the manuscript. All authors read and approved the final manuscript.

**Funding**

No funding was received.

**Availability of data and materials**

Please contact the author for data requests.

**Ethics approval and consent to participate**

All procedures performed in this study involving human participants were in accordance with the ethical standards of the COMERE research committee ("commission interne de méthodologie en recherche clinique", ICMART201801).

**Consent for publication**

Consent has been obtained from participants to publish this work.

**Competing interests**

ED was invited in 2018 to attend an international meeting by AAA who sponsored registration fees and travel expenses. Otherwise, there are no potential conflicts of interest to disclose for any of the authors.

**Author details**

Received: 19 June 2020 Accepted: 23 November 2020

Published online: 04 January 2021

**References**

1. Huizing DMV, de Wit-van der Veen BJ, Verheij M, Stokkel MPM. Dosimetry methods and clinical applications in peptide receptor radionuclide therapy for neuroendocrine tumours: a literature review. EJNMMI Res (Internet). 2018. https://doi.org/10.1186/s13550-018-0443-z.

2. Flux GD, Sjogreen Gleisner K, Chiesa C, Lassmann M, Chouin N, Gear J, et al. From fixed activities to personalized treatments in radionuclide therapy: lost in translation? Eur J Nucl Med Mol Imaging. 2018;45(1):152-4.

3. Sundlöv A, Gustafsson J, Brolin G, Mortensen N, Hermann R, Bernhardt P, et al. Feasibility of simplifying renal dosimetry in 177Lu peptide receptor radionuclide therapy. EJNMMI Phys (Internet). 2018. https://doi.org/10.1186/s40658-018-0210-2.

4. Marin G, Vanderlinden B, Karfis I, Guiot T, Wimana Z, Reynaert N, et al. A dosimetry procedure for organs-at-risk in 177Lu peptide receptor radionuclide therapy of patients with neuroendocrine tumours. Phys Med. 2018;56:41-9.

5. Del Prete M, Arsenault F, Saighi N, Zhao W, Buteau F-A, Celler A, et al. Accuracy and reproducibility of simplified QSPECT dosimetry for personalized 177Lu-octreotate PRRT. EJNMMI Phys. 2018;5(1):25.

6. Johnson TK, McClure D, McCourt S. MABDOSE II: validation of a general purpose dose estimation code. Med Phys. 1999;26(7):1396-403.

7. Gardin I, Bouchet LG, Assié K, Caron J, Lisbona A, Ferrer L, et al. Voxeldose: a computer program for 3-D dose calculation in therapeutic nuclear medicine. Cancer Biother Radiopharm. 2003;18(1):109-15.

8. Grimes J, Uribe C, Celler A. JADA: a graphical user interface for comprehensive internal dose assessment in nuclear medicine: a GUI for comprehensive internal dose assessment. Med Phys. 2013;40(7):072501.

9. Kletting P, Schimmel S, Hänscheid H, Luster M, Fernández M, Nosske D, et al. The NUKDOS software for treatment planning in molecular radiotherapy. Z Für Med Phys. 2015;25(3):264-74.

10. Li T, Zhu L, Lu Z, Song N, Lin K-H, Mok GSP. BIGDOSE: software for 3D personalized targeted radionuclide therapy dosimetry. Quant Imaging Med Surg. 2020;10(1):160-70.

11. Mora-Ramirez E, Santoro L, Cassol E, Ocampo-Ramos JC, Clayton N, Kayal G, et al. Comparison of commercial dosimetric software platforms in patients treated with 177Lu-DOTATATE for peptide receptor radionuclide therapy. Med Phys. 2020;1:mp.14375.

12. Santoro L, Mora-Ramirez E, Trauchessec D, Chouaf S, Eustache P, Pouget J-P, et al. Implementation of patient dosimetry in the clinical practice after targeted radiotherapy using [177Lu-[DOTAO, Tyr3]-octreotate. EJNMMI Res. 2018;8(1):103.

13. GE Healthcare. Organ dose estimates for radio-isotope therapy treatment planning purposes. Dosimetry toolkit package. White paper. 2011.

14. Stabin MG, Sparks RB, Crowe E. OLINDA/EXM: the second-generation personal computer software for internal dose assessment in nuclear medicine. J Nucl Med. 2005;46(6):1023-7.

15. Kupitz D, Wetz C, Wissel H, Wedel F, Apostolova I, Wallbaum T, et al. Software-assisted dosimetry in peptide receptor radionuclide therapy with 177Lutetium-DOTATATE for various imaging scenarios. Mosley RL, editor. PLoS ONE. 2017;12(11):e0187570.

16. Huizing DMV, Peters SMB, Versleijen MWJ, Martens E, Verheij M, Sinaasappeł M, et al. A head-to-head comparison between two commercial software packages for hybrid dosimetry after peptide receptor radionuclide therapy. EJNMMI Phys. 2020;7(1):36.

17. Maughan NM, Garcia-Ramirez J, Arpidone M, Swallen A, Laforest R, Goddu SM, et al. Validation of post-treatment PET-based dosimetry software for hepatic radioembolization of Yttrium-90 microspheres. Med Phys. 2019;46(5):2394-402.

18. Hippeläinen ET, Tenhunen MJ, Mäenpää HO, Heikkonen JJ, Sohlberg AO. Dosimetry software Hermes Internal Radiation Dosimetry: from

<footnote n="1">¹ Nuclear Medicine Department, Montpellier Cancer Institute (ICM), Univ. Montpellier, 208 Avenue des Apothicaires, 34298 Montpellier Cedex 5, France.</footnote>

<footnote n="2">² Centre de Recherche en Cancérologie de Toulouse, Toulouse, France.</footnote>

<footnote n="3">³ INSERM, UMR 1037, Toulouse III Paul Sabatier University, Toulouse, France.</footnote>

<footnote n="4">⁴ University of Costa Rica, Physics School, CICANUM, San José, Costa Rica.</footnote>

<footnote n="5">⁵ Montpellier Cancer Research Institute, UMR 1194, Univ. Montpellier, Montpellier, France.</footnote>
quantitative image reconstruction to voxel-level absorbed dose distribu-
tion. Nucl Med Commun. 2017;38(5):357–65.

19. Barna S, Haug A, Hartenbach M, Rasul S, Grubmüller B, Kramer G, et al. Dose calculations and dose-effect relationships in 177Lu-PSMA I&T radionuclide therapy for metastatic castration-resistant prostate cancer. Clin Nucl Med. 2020;45(9):661–7.

20. Andersson M, Johansson L, Eckerman K, Mattsson S. IDAC-Dose 2.1, an internal dosimetry program for diagnostic nuclear medicine based on the ICRP adult reference voxel phantoms. EJNMMI Res. 2017;7(1):88.

21. Pasciak AS, Bourgeois AC, Bradley YC. A comparison of techniques for 90Y PET/CT image-based dosimetry following radioembolization with resin microspheres. Front Oncol (Internet). 2014. https://doi.org/10.3389/fonc.2014.00121/abstract.

22. Dieudonne A, Hobbs RF, Lebtahi R, Maurel F, Baechler S, Wahl RL, et al. Study of the impact of tissue density heterogeneities on 3-dimensional abdominal dosimetry: comparison between dose kernel convolution and direct Monte Carlo methods. J Nucl Med. 2013;54(2):236–43.

23. Dieudonne A, Hobbs RF, Bolch WE, Sgouros G, Gardin I. Fine-resolution voxel S values for constructing absorbed dose distributions at variable voxel size. J Nucl Med. 2010;51(10):1600–7.

24. Kondev FG. Lu-177\_comments on evaluation of decay data for beta-decay; 2002.

25. Simpkin DJ. Radiation Interactions and Internal Dosimetry in Nuclear Medicine; 1999;19(1):13.

26. Lin LI-K. A Concordance Correlation Coefficient to Evaluate Reproducibility. Biometrics. 1989;45(1):255.

27. Altman DG, Bland JM. Measurement in medicine: the analysis of method comparison studies. Statistician. 1983;32(3):307.

28. Zhao W, Esquinas PL, Hou X, Uribe CF, Gonzalez M, Beauregard J-M, et al. Determination of gamma camera calibration factors for quantitation of therapeutic radioisotopes. EJNMMI Phys (Internet). 2018. https://doi.org/10.1186/s40658-018-0208-9.

29. Uribe CF, Esquinas PL, Tanguay J, Gonzalez M, Gaudin E, Beauregard J-M, et al. Accuracy of 177Lu activity quantification in SPECT imaging: a phantom study. EJNMMI Phys (Internet). 2017. https://doi.org/10.1186/s40658-016-0170-3.

30. Zimmerman BE, Grošev D, Buvat I, Coca Pérez MA, Frey EC, Green A, et al. Multi-centre evaluation of accuracy and reproducibility of planar and

SPECT image quantification: an IAEA phantom study. Z Für Med Phys. 2017;27(2):98–112.

31. Peters SMB, Meyer Viol SL, van der Werf NR, de Jong N, van Velden FHP, Meeuwis A, et al. Variability in lutetium-177 SPECT quantification between different state-of-the-art SPECT/CT systems. EJNMMI Phys. 2020;7(1):9.

32. Lassmann M, Eberlein U. The relevance of dosimetry in precision medi-
cine. J Nucl Med. 2018;2018:1.

33. Tran-Gia J, Lassmann M. Optimizing image quantification for 177 Lu SPECT/CT based on a 3D printed 2-compartment kidney phantom. J Nucl Med. 2018;59(4):616–24.

34. Gustafsson J, Sundlöv A, Sjögreen GK. SPECT image segmentation for estimation of tumour volume and activity concentration in 177Lu-DOTA-TATE radionuclide therapy. EJNMMI Res. 2017;7(1):18.

35. Grassi E, Fioroni F, Ferri V, Mezzenga E, Sarti MA, Paulus T, et al. Quantitative comparison between the commercial software STRATOS® by Philips and a homemade software for voxel-dosimetry in radiopeptide therapy. Phys Med. 2015;31(1):72–9.

36. Grassi E, Fioroni F, Berenato S, Patterson N, Ferri V, Braglia L, et al. Effect of image registration on 3D absorbed dose calculations in 177Lu-DOTATOC peptide receptor radionuclide therapy. Phys Med. 2018;45:177–85.

37. Hippeläinen E, Tenhunen M, Sohlberg A. Fast voxel-level dosimetry for 177 Lu labelled peptide treatments. Phys Med Biol. 2015;60(17):6685–700.

38. Sandstrom M, Garske-Roman U, Granberg D, Johansson S, Widstrom C, Eriksson B, et al. Individualized dosimetry of kidney and bone marrow in patients undergoing 177Lu-DOTA-octreotate treatment. J Nucl Med. 2013;54(1):33–41.

39. Gear JI, Cox MG, Gustafsson J, Gleisner KS, Murray I, Glatting G, et al. EANM practical guidance on uncertainty analysis for molecular radiotherapy absorbed dose calculations. Eur J Nucl Med Mol Imaging (Internet). 2018. https://doi.org/10.1007/s00259-018-4136-7.

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in pub-
lished maps and institutional affiliations.

Submit your manuscript to a SpringerOpen®
journal and benefit from:

► Convenient online submission

► Rigorous peer review

► Open access: articles freely available online

► High visibility within the field

► Retaining the copyright to your article

Submit your next manuscript at ► springeropen.com