# Optimization of brain PET imaging for a multicentre trial: the French CATI experience

Marie-Odile Habert, Sullivan Marie, Hugo Bertin, Moana Reynal, Jean-Baptiste Martini, Mamadou Diallo, Aurélie Kas, Régine Trébossen

► To cite this version:

Marie-Odile Habert, Sullivan Marie, Hugo Bertin, Moana Reynal, Jean-Baptiste Martini, et al.. Optimization of brain PET imaging for a multicentre trial: the French CATI experience. EJNMMI Physics, 2016, 3 (1), 10.1186/s40658-016-0141-8. hal-01300073

HAL Id: hal-01300073

https://hal.sorbonne-universite.fr/hal-01300073

Submitted on 8 Apr 2016

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

# Optimization of brain PET imaging for a multicentre trial: the French CATI experience

Marie-Odile Habert<fnref n="1" />,<fnref n="2" />,³,⁵*, Sullivan Marie<fnref n="1" />,<fnref n="2" />, Hugo Bertin<fnref n="1" />,<fnref n="2" />, Moana Reynal<fnref n="1" />,<fnref n="2" />, Jean-Baptiste Martini<fnref n="1" />,<fnref n="2" />, Mamadou Diallo<fnref n="1" />,<fnref n="2" />, Aurélie Kas<fnref n="1" />,<fnref n="2" />,³, the CATI Nuclear Medicine network and Régine Trébossen⁴

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

CrossMark

**Abstract**

**Background:** CATI is a French initiative launched in 2010 to handle the neuroimaging of a large cohort of subjects recruited for an Alzheimer's research program called MEMENTO. This paper presents our test protocol and results obtained for the 22 PET centres (overall 13 different scanners) involved in the MEMENTO cohort. We determined acquisition parameters using phantom experiments prior to patient studies, with the aim of optimizing PET quantitative values to the highest possible per site, while reducing, if possible, variability across centres.

**Methods:** Jaszczak's and 3D-Hoffman's phantom measurements were used to assess image spatial resolution (ISR), recovery coefficients (RC) in hot and cold spheres, and signal-to-noise ratio (SNR). For each centre, the optimal reconstruction parameters were chosen as those maximizing ISR and RC without a noticeable decrease in SNR. Point-spread-function (PSF) modelling reconstructions were discarded. The three figures of merit extracted from the images reconstructed with optimized parameters and routine schemes were compared, as were volumes of interest ratios extracted from Hoffman acquisitions. The net effect of the 3D-OSEM reconstruction parameter optimization was investigated on a subset of 18 scanners without PSF modelling reconstruction.

**Results:** Compared to the routine parameters of the 22 PET centres, average RC in the two smallest hot and cold spheres and average ISR remained stable or were improved with the optimized reconstruction, at the expense of slight SNR degradation, while the dispersion of values was reduced.

For the subset of scanners without PSF modelling, the mean RC of the smallest hot sphere obtained with the optimized reconstruction was significantly higher than with routine reconstruction. The putamen and caudate-to-white matter ratios measured on 3D-Hoffman acquisitions of all centres were also significantly improved by the optimization, while the variance was reduced.

**Conclusions:** This study provides guidelines for optimizing quantitative results for multicentric PET neuroimaging trials.

**Keywords:** Neurology, PET/CT, Quality assurance, Standardization, Multicentre trials

<footnote>* Correspondence: habert@upmc.fr</footnote>

<footnote n="1">¹Sorbonne Universités, UPMC Univ Paris 06, CNRS, INSERM, Laboratoire d'Imagerie Biomédicale, F-75013 Paris, France</footnote>

<footnote n="2">²Centre pour l'Acquisition et le Traitement des Images (www.cati-neuroimaging.com), Paris, France</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
## Background

CATI is a French platform funded in 2010 with the aim of supporting multicentre clinical trials involving neuroimaging (cati-neuroimaging.com). The main project devoted to CATI was to handle the multimodality imaging aspect of a large cohort of subjects included in a research program on Alzheimer's disease (AD). This program, called MEMENTO, is a longitudinal study aiming at better understanding the natural history of Alzheimer's disease. It has enrolled 2300 subjects with either light cognitive deficits or cognitive complaints that will be followed up for at least 5 years (http://clinicaltrials.gov/show/NCT01926249).

The imaging aspect of the MEMENTO study includes the acquisition of PET cerebral glucose metabolism and amyloid distribution, with MRI studies on anatomy, diffusion, and resting state. In this paper, we focus on the optimization of PET quantitative values to the highest possible per site, while reducing, if possible, the variability across centres for the MEMENTO cohort. Phantom data were acquired during a technical visit to the PET centres, prior to patients' participation in the trial. The same phantom data were also acquired during follow-up visits after 18 months. The present study only addresses the data acquired during the site technical set-up.

Several large cohort studies on AD have been launched since the early 2000s. Among them, ADNI-1 was a forward-thinking, imaging-based multicentre clinical trial that involved 50 centres in North America, with the aim of identifying biomarkers of AD. For the PET imaging carried out in this study, the ADNI PET core determined reconstruction parameters for each scanner model and 3D-Hoffman phantoms were acquired with a standard protocol. However, discrepancies between the image characteristics from different PET centres remained high and were accounted for by degrading the spatial resolution to the lowest value among centres [1].

Other multicentre studies specifically focusing on AD have been conducted [2, 3]. In contrast to the ADNI trial, no phantom acquisition was performed prior to the study, and no post-filtering was applied to the data. Differences in spatially normalized FDG-PET scans obtained with scanners of different resolutions were minimized by the following measures: (i) restricting the analysis to voxels with intensity 80 % greater than the whole-brain mean, and (ii) excluding voxels from the uppermost 10 slices (i.e., from the top 22.5 mm of the brain) and from the lowermost 5 slices, where significant inter-scanner effects due to different fields of view have been reported [3].

The most advanced domain for PET image acquisition harmonization across multiple centres is oncological imaging. Guidelines for PET centres involved in clinical trials have been written by a group of experts, under the umbrella of the European society of Nuclear Medicine [4, 5], or the American College of Radiology Imaging Network [6]. These guidelines are based on quality checks with phantoms reproducing the imaging conditions encountered in the abdomen and the thorax in glucose metabolism studies. In addition to these tomography quality checks, to ensure the highest reproducibility of the measures extracted from the image for patient follow-up, disease evolution, and recovery after treatments, these groups have formulated recommendations for data analysis. Overall, these recommendations enable the cross-comparison of measurements from different centres. According to our knowledge, the latest guidelines for FDG-PET brain imaging were published in 2009 [7] and do not meet our objectives.
This paper presents our test protocol and the results obtained for the 22 PET centres involved in the MEMENTO cohort. French PET centres were equipped with systems set up between 2003 and 2012. No brain-dedicated high-resolution system or 2D tomography systems were included. Therefore, we chose to optimize acquisition and reconstruction parameters using phantom experiments prior to patient studies, with the aim of optimizing PET quantitative values to the highest possible per site, while reducing, if possible, the variability across centres. We also checked the impact of this harmonization on signal-to-noise ratios, as well as the exclusion of reconstruction with point-spread-function (PSF) modelling.

# Methods

## Phantom studies

Two phantom studies were acquired for the qualification process of the 22 centres, the Jaszczak phantom and the Hoffman 3D brain phantom. Standardized uptake value (SUV) measurements were also checked to assess the cross-calibration between the PET scanner and the dose calibration system.

The Deluxe Jaszczak phantom (model ECT/DLX/P) is a Plexiglas cylinder partially composed of cold rods (diameters 4.8 to 12.7 mm) and partially of six inserted hollow spheres. The main cylinder of 6060 mL was filled with a 5 kBq mL⁻¹ FDG solution. Four spheres (internal diameters 7.86, 12.43, 15.43, and 24.82 mm) were filled with a 15 kBq mL⁻¹ FDG solution. The remaining two spheres (internal diameters 9.89 and 31.27 mm) were filled with cold water.

The Hoffman 3D brain phantom (model BR/3D/P) was filled with a 37-to-55 kBq mL⁻¹ FDG solution, taking great care to avoid the presence of bubbles.

A dynamic acquisition of 3 × 5 min was performed for each phantom. CT acquisition parameters were set such that the effective dose would be low (≤0.3 mSv), in accordance with the recommendations from the European Association of Nuclear Medicine [7].

The Hoffman phantom was positioned on a head holder when available, with the Deluxe Jaszczak phantom (20–30 kBq mL⁻¹ FDG solution) below it for the simulation of a diffusion medium, as in subjects. Different reconstruction parameters were tested with varying reconstruction algorithms, number of iterations and subsets, reconstruction diameters, matrix sizes, or filters, always including those routinely used in brain studies. For the first 5 pilot centres, up to 10 sets of reconstructions were tested. Attenuation (derived from a CT scan) and scatter were corrected using standard software supplied by the scanner manufacturers.

Acquisition from a uniform cylinder filled with 5 kBq mL⁻¹ FDG solution was also performed to verify cross-calibration.

The same team of two technologists certified for radiation protection prepared all the phantoms. The CATI PET project manager supervised the whole procedure and the setting of reconstruction parameters.

## Analysis of the Deluxe Jaszczak phantom studies

### Jaszczak phantom studies

The different image reconstruction schemes were compared based on recovery coefficients (RC) computed for each sphere and spatial resolution estimated from the cold
rods. RC were calculated from measurements in volumes of interest (VOIs) defined on each sphere and on the background. In-house software was developed for this purpose. It includes the following steps (all displayed for quality check):

1) For all sets of reconstructed PET volumes, the computation of the spheres' centre of gravity, based on an automatic segmentation of the spheres.

2) The computation of the spheres' centre of gravity based on the automatic segmentation of the spheres on CT scan volume, using the Hough algorithm [8], optimized with the incorporation of the spheres' specifications.

3) VOIs were defined using both the CT spheres' centre of gravity (determined in step 2) and the known sphere diameters, corrected if necessary for a mismatch between PET and CT spheres, then eroded to exclude the sphere wall. The background activity concentration was measured in VOIs obtained after the rotation of the VOIs drawn on the spheres.

4) The measurement of the mean activity in each sphere and background (BG), and calculation of RC (RC<sub>H</sub> and RC<sub>C</sub>, respectively, for hot and cold spheres):

$$RC_H = \frac{(M_S - M_{BG}) / (M_S + M_{BG})}{(A_S - A_{BG}) / (A_S + A_{BG})}$$

$$RC_C = \frac{(M_{BG} - M_S) / (M_S + M_{BG})}{(A_{BG} - A_S) / (A_S + A_{BG})}$$

where M<sub>S</sub> is the activity measured in the sphere, M<sub>BG</sub> is the activity measured in the background, A<sub>S</sub> is the actual activity injected in the sphere, and A<sub>BG</sub> is the actual activity injected in the background.

5) Signal-to-noise ratio was calculated as follows:

$$SNR = \frac{M_{S_{12.43}} - M_{BG}}{SD_{BG}}$$

where M<sub>S<sub>12.43</sub></sub> is the mean activity measured in the hot sphere of 12.43 mm, M<sub>BG</sub> is the mean activity measured in the background, and SD<sub>BG</sub> is the standard deviation in the background.

6) Display of the cold rod images and spheres for visual assessment of noise and spatial resolution, and of tables and graphics for all obtained RC according to the different reconstruction parameters (Fig. 1).

We also computed the image spatial resolution (ISR) using a method developed by Prieto et al. [9] to estimate the full width at half maximum (FWHM) of each PET scanner from the three largest hot spheres. The process can be summarized in three main steps:

1) Calculation of a ratio R<sub>mes</sub> as follows:
<figure><img src="image_3.png" /></figure>

$$R_{\text{mes}} = \frac{\text{Max}_S}{7A_S}$$

where Max_S is the maximum activity measured in the sphere and A_S is the actual activity injected in the hot spheres. Each R_mes was normalized by the corresponding value of the largest sphere in order to limit the bias introduced by attenuation and scatter corrections.

2) Computation of a mask of the three spheres at the same position in the scanner field of view as during the acquisition; this mask was then convolved with a three-dimensional isotropic Gaussian function with different FWHM values ranging from 4 to 10 mm and a step of 0.1 mm. A theoretical ratio R_the was then calculated for each of the three spheres as follows:

$$R_{\text{the}} = \frac{\text{Max}_S}{A_S}$$

where Max_S is the maximum value measured in the sphere after convolution and A_S is the actual value.
3) The optimal ISR value was selected as the FWHM value that best fitted the experimental values, according to minimization of the normalized sum of squared differences between observed $R_{mes}$ and simulated $R_{the}$ over the three spheres.

## Analysis of the Hoffman 3D brain phantom studies

All reconstructed volumes were co-registered to the highest quality CT scan of the phantom (acquired at centre no. 16) using SPM8. The grey matter ribbon was obtained by thresholding CT image intensity at 60 % of the maximum value. The CT scan was also manually segmented in volumes of interest (Fig. 2) in order to compute right-to-left (R/L) and anterior-to-posterior (L/P) cortical ratios, as well as caudate-to-white matter (C/ WM), putamen-to-white matter (P/WM), and grey-to-white matter (GM/WM) ratios from PET images. The upper and lower slices were discarded for the analyses.

## Image reconstruction

First, the three frames were summed. We then systematically compared phantom measurements obtained from images reconstructed with the algorithm selected for diagnostic purposes by each centre with the images reconstructed using the parameters optimized by CATI. Reconstruction algorithms incorporating the modelling of the spatial resolution of the tomographs were used for diagnostic purposes in four (4/22) centres. These algorithms were discarded by CATI to avoid possible additional centre

<figure><img src="image_4.png" /><figcaption>**Fig. 2** Volumes of interest used for analysis of the Hoffman 3D brain phantom acquisitions. **a** White matter. **b** Caudate (**blue**) and putamen (**red**). **c** Right and left cortex. **d** Anterior and posterior cortex</figcaption></figure>
effects, and 3D-OSEM statistical image reconstruction algorithms, or FORE 2D-OSEM when the 3D-OSEM was not available, were selected.

The optimization was based on three figures of merit extracted from measurements on a Jaszczak phantom, the RC of spheres of various diameters, the signal-to-background-noise ratio (SNR) and the ISR. The image reconstruction parameters were chosen at each site as a compromise between maximized RC values, ISR, and SNR. Volumes of interest ratios extracted from Hoffman acquisitions were also compared to confirm the choice of reconstruction parameters for each site. Our approach can be described as follows:

1) Reconstruction matrix was set to obtain pixel spacing inferior to 3 mm

2) The number of total iterations multiplied by the number of subsets was chosen such as the product iterations × subsets was superior to 50, and optimized with post-reconstruction smoothing

3) RC and ISR obtained with different combinations of reconstruction parameters were compared first: the parameters giving highest values with acceptable SNR were chosen.

4) Hoffmann images were then quantitatively and visually checked for the best compromise between spatial resolution and noise.

5) The optimized quantitative values for both phantoms obtained at each centre were finally compared to the routine values.

The optimized parameters chosen according to the model of scanner are presented in Table 1.

## Statistical analyses

The Wilcoxon signed rank test was used to compare routine and optimized values obtained from phantoms' studies. Variances were compared with a Pitman test [10]. A significance threshold of 0.05 was adopted for all statistical analyses. Statistical analyses were performed for all centres and also for a subset of 18 centres where PSF modelling reconstruction was not available in order to assess the impact of PSF modelling on the final optimization.

## Results

All tomographs were 3D PET/CT, installed between 2003 and 2012, and consisted of 9 GE, 3 Philips, and 10 Siemens systems, with a total of 13 different models of scanners (Table 1).

In four centres (nos. 8, 13, 19, and 22), an error greater than 10 % was found for cross-calibration measurements and was later corrected by the centre's physicist.

Four other scanners (nos. 10, 11, 14, and 20) presented a mismatch of more than 3 mm between CT and PET images, which required fixing by servicing the scanners.

### Jaszczak phantom results

For all centres, the optimized RC of the two smallest hot spheres ranged from 0.17 to 0.52 (mean ± SD 0.33 ± 0.09) and from 0.56 to 0.86 (mean ± SD 0.71 ± 0.08), respectively
<table><caption>**Table 1** Optimized parameters chosen according to the model of the scanner</caption><thead><tr><th>GE models</th><th>Number</th><th>Installation year</th><th>Slice thickness (mm)</th><th>Recons. method</th><th>Number of iterations (iterations x subsets)</th><th>Post-reconstruction filter type and FWHM if relevant (mm)</th><th>Matrix</th><th>Pixel spacing (mm x mm)</th><th>Others</th></tr></thead><tbody><tr><td>Discovery 690</td><td>3</td><td>2009/2010/2011</td><td>3.27</td><td>VPHD</td><td>10 x 36</td><td>Gaussian 3.5</td><td>256 x 256</td><td>1.17 x 1.17</td><td>Z-filter: standard</td></tr><tr><td>Discovery RX</td><td>2</td><td>2009</td><td>3.27</td><td>FORE + OSEM 2D</td><td>7 x 35</td><td>Gaussian 3.5</td><td>128 x 128</td><td>2.34 x 2.34</td><td>Loop filter: 2.34 mm<br>Z-filter: standard</td></tr><tr><td>Discovery DST - E</td><td>1</td><td>2004</td><td>3.27</td><td>FORE + OSEM 2D</td><td>5 x 35</td><td>Gaussian 2.57</td><td>128 x 128</td><td>2.34 x 2.34</td><td>Loop filter: 2.34 mm<br>Z-filter: standard</td></tr><tr><td>Discovery ST 4</td><td>1</td><td>2004</td><td>3.27</td><td>FORE + OSEM 2D</td><td>7 x 35</td><td>Gaussian 3</td><td>256 x 256</td><td>1.17 x 1.17</td><td>Loop filter: 2 mm<br>Z-filter: standard</td></tr><tr><td>Discovery ST</td><td>2</td><td>2003/2004</td><td>3.27</td><td>FORE + OSEM 2D</td><td>7 x 32</td><td>Gaussian 2 or 2.57</td><td>128 x 128</td><td>2.34 x 2.34</td><td>Loop filter: 2.34 mm</td></tr><tr><th>Philips models</th><th>Number</th><th>Installation year</th><th>Slice thickness (mm)</th><th>Recons. method</th><th>Number of iterations (iterations x subsets)</th><th>Smoothing (mm)</th><th>Matrix</th><th>Pixel spacing (mm x mm)</th><th>Others</th></tr><tr><td>Gemini TF</td><td>2</td><td>2008/2010</td><td>2</td><td>LOR - RAMLA</td><td>10 x 33</td><td>Smooth B</td><td>128 x 128</td><td>2 x 2</td><td>-</td></tr><tr><td>Gemini GXL</td><td>1</td><td>2006</td><td>2</td><td>LOR - RAMLA</td><td>10 x NA</td><td>Smooth</td><td>128 x 128</td><td>2 x 2</td><td>-</td></tr><tr><th>Siemens models</th><th>Number</th><th>Installation year</th><th>Slice thickness (mm)</th><th>Recons. method</th><th>Number of iterations (iterations x subsets)</th><th>Post-reconstruction filter type and FWHM if relevant (mm)</th><th>Matrix</th><th>Pixel spacing (mm x mm)</th><th>Others</th></tr><tr><td>Biograph mCT 40/64</td><td>3</td><td>2009/2012</td><td>2.027</td><td>OSEM 3D</td><td>12 x 24</td><td>Gaussian 3 or 4</td><td>256 x 256</td><td>1.59 x 1.59</td><td>-</td></tr><tr><td>Biograph Hirez TruePoint</td><td>1</td><td>2008</td><td>3</td><td>OSEM 3D</td><td>8 x 21</td><td>Gaussian 4</td><td>336 x 336</td><td>1.02 x 1.02</td><td>-</td></tr><tr><td>Biograph 6 VB 20B True V</td><td>1</td><td>2006</td><td>2</td><td>FORE + OSEM 2D</td><td>8 x 21</td><td>Gaussian 4</td><td>256 x 256</td><td>1.34 x 1.34</td><td>-</td></tr><tr><td>Biograph 6</td><td>2</td><td>2004/2005</td><td>2</td><td>FORE + OSEM 2D</td><td>8 x 24</td><td>Gaussian 2 or 3</td><td>256 x 256</td><td>1.33 x 1.33</td><td>-</td></tr><tr><td>Biograph 16 VB40B</td><td>1</td><td>2004</td><td>2</td><td>FORE + OSEM 2D</td><td>8 x 24</td><td>Gaussian 3</td><td>336 x 336</td><td>1.02 x 1.02</td><td>-</td></tr><tr><td>Biograph LSO DUO</td><td>2</td><td>2003/2004</td><td>3.375</td><td>FORE + OSEM 2D</td><td>8 x 16</td><td>Gaussian 2</td><td>256 x 256</td><td>1.33 x 1.33</td><td>-</td></tr></tbody></table>

NA non applicable, OSEM Ordered Subset Expectation Maximization, VPHD VUE Point High Definition, LOR-RAMLA Line-of-Response Row-Action Maximum Likelihood Algorithm, FORE Fourier Rebinning
(Fig. 3). Mean RCs were not significantly different for routine and optimized reconstructions, but the variance was significantly reduced (*p* = 0.008 and *p* = 0.002, respectively). The optimized RC of the two cold spheres ranged from 0.22 to 0.52 (mean ± SD 0.41 ± 0.08) for the smallest and from 0.60 to 0.84 for the largest (mean ± SD 0.77 ± 0.05). Mean RC significantly increased with optimized reconstructions (*p* = 0.003 and *p* = 0.004, respectively) (Fig. 4). Variance was also significantly reduced for the largest sphere (*p* = 0.002). All RC values obtained with routine and optimized parameters are presented in Table 2, with the exception of four missing data from one centre, because of an operating error.

For the 18 centres where PSF modelling reconstruction was not available, the optimized RC of the two smallest hot spheres ranged from 0.17 to 0.52 (mean ± SD 0.31 ± 0.09) and from 0.56 to 0.86 (mean ± SD 0.70 ± 0.08), respectively (Fig. 3). Mean RC of the smallest hot sphere was significantly higher for optimized reconstructions *p* = 0.04), while the dispersion of values remained equivalent for both spheres. The optimized RC of the two cold spheres ranged from 0.22 to 0.52 (mean ± SD 0.40 ± 0.08) for the smallest and from 0.60 to 0.84 for the largest (mean ± SD 0.76 ± 0.05) (Fig. 4). Mean RC significantly increased with optimized reconstructions (*p* = 0.00002 and *p* = 0.001, respectively), but no significant difference was found for variances.

No significant difference could be detected (*p* = 0.61) between the values of the ratio “measured background/actual background” for routine reconstruction parameters (0.97 ± 0.13) and optimized reconstruction parameters (0.99 ± 0.08). In all centres, the SNR decreased between routine and optimized parameters, from 14.2 ± 6.16 to 9.74 ±

<figure><img src="image_5.png" /><figcaption>**Fig. 3** Jaszczak phantom results. Recovery coefficients (RC) obtained with routine and optimized acquisition and reconstruction parameters for the two smallest hot spheres, for all centres (upper row) and for the subset of 18 centres without PSF modelling (lower row). *P* values represent the significant test results either for comparison of means (Wilcoxon test) or for comparison of standard deviations (Pitman test)</figcaption></figure>
<figure><img src="image_6.png" /></figure>

2.41 (mean ± SD). In the subset of 18 centres, it also decreased from 12.67 ± 4.52 to 9.70 ± 2.5 (mean ± SD).

## Image spatial resolution measurements

For all centres, the ISR ranged from 5.9 to 7.7 mm (mean ± SD 6.76 ± 0.6) with optimized parameters and from 5.5 to 7.7 mm (mean ± SD 6.98 ± 0.6) with routine parameters and did not differ significantly.

## 3D Hoffman phantom results

For all centres, a significant increase of the mean C/WM ($p = 0.005$) and P/WM ($p = 0.04$) ratios was observed between routine and optimized reconstructions, but not for the other ratios. The variance was significantly reduced for the GM/WM, C/WM, and P/WM ratios ($p = 0.001$; $p = 0.005$; $p = 0.006$, respectively) (Fig. 5). In the subset of 18 centres, a significant increase of both mean C/WM ($p = 0.0002$) and P/WM ratios ($p = 0.006$) was observed. Variances were not significantly different. The ratio values for both reconstructions and all centres are presented in Table 3.

## Discussion

Multicentre clinical trials involving imaging require procedures for image acquisition and reconstruction to be optimized to account for inter-subject profile variability and to ensure the robustness of the analysis after pooling the data and for patient follow-
<table><caption>**Table 2** Recovery coefficient values for both reconstructions and all centres obtained from Jaszczak acquisitions</caption><thead><tr><th rowspan="2">Centre number</th><th colspan="2">7.86 mm Hot sphere</th><th colspan="2">12.43 mm Hot sphere</th><th colspan="2">14.43 mm Hot sphere</th><th colspan="2">24.82 mm Hot sphere</th><th colspan="2">9.89 mm Cold sphere</th><th colspan="2">24.82 mm Cold sphere</th></tr><tr><th>RP</th><th>OP</th><th>RP</th><th>OP</th><th>RP</th><th>OP</th><th>RP</th><th>OP</th><th>RP</th><th>OP</th><th>RP</th><th>OP</th></tr></thead><tbody><tr><td>1</td><td>0.17</td><td>0.25</td><td>0.52</td><td>0.57</td><td>0.56</td><td>0.59</td><td>0.74</td><td>0.74</td><td>0.20</td><td>0.34</td><td>0.57</td><td>0.81</td></tr><tr><td>2</td><td>0.25</td><td>0.29</td><td>0.64</td><td>0.68</td><td>0.76</td><td>0.77</td><td>0.95</td><td>0.95</td><td>0.24</td><td>0.32</td><td>0.64</td><td>0.79</td></tr><tr><td>3</td><td>0.24</td><td>0.26</td><td>0.59</td><td>0.66</td><td>0.70</td><td>0.75</td><td>0.91</td><td>0.92</td><td>0.40</td><td>0.50</td><td>0.74</td><td>0.74</td></tr><tr><td>4</td><td>0.23</td><td>0.25</td><td>0.55</td><td>0.56</td><td>0.74</td><td>0.73</td><td>0.93</td><td>0.91</td><td>0.25</td><td>0.29</td><td>0.76</td><td>0.75</td></tr><tr><td>5</td><td>0.29</td><td>0.25</td><td>0.76</td><td>0.70</td><td>0.78</td><td>0.74</td><td>0.94</td><td>0.92</td><td>0.33</td><td>0.37</td><td>0.67</td><td>0.76</td></tr><tr><td>6</td><td>0.49</td><td>0.52</td><td>0.85</td><td>0.86</td><td>0.90</td><td>0.91</td><td>1.01</td><td>1.01</td><td>0.33</td><td>0.46</td><td>0.64</td><td>0.77</td></tr><tr><td>7</td><td>0.37</td><td>0.46</td><td>0.79</td><td>0.78</td><td>0.89</td><td>0.82</td><td>1.04</td><td>0.97</td><td>0.33</td><td>0.52</td><td>0.70</td><td>0.78</td></tr><tr><td>8</td><td>0.20</td><td>0.24</td><td>0.61</td><td>0.62</td><td>0.70</td><td>0.70</td><td>0.92</td><td>0.91</td><td>0.22</td><td>0.22</td><td>0.64</td><td>0.60</td></tr><tr><td>9</td><td>0.30</td><td>0.41</td><td>NA</td><td>NA</td><td>0.87</td><td>0.82</td><td>1.06</td><td>0.99</td><td>NA</td><td>NA</td><td>0.73</td><td>0.80</td></tr><tr><td>10</td><td>0.15</td><td>0.24</td><td>0.71</td><td>0.80</td><td>0.80</td><td>0.85</td><td>0.99</td><td>0.99</td><td>0.29</td><td>0.42</td><td>0.64</td><td>0.78</td></tr><tr><td>11</td><td>0.15</td><td>0.38</td><td>0.63</td><td>0.75</td><td>0.68</td><td>0.78</td><td>0.92</td><td>0.95</td><td>0.26</td><td>0.44</td><td>0.68</td><td>0.84</td></tr><tr><td>12</td><td>0.17</td><td>0.17</td><td>0.71</td><td>0.71</td><td>0.73</td><td>0.73</td><td>0.91</td><td>0.91</td><td>0.41</td><td>0.41</td><td>0.74</td><td>0.74</td></tr><tr><td>13</td><td>0.67</td><td>0.41</td><td>1.02</td><td>0.77</td><td>0.97</td><td>0.79</td><td>1.03</td><td>0.94</td><td>0.57</td><td>0.41</td><td>0.89</td><td>0.78</td></tr><tr><td>14</td><td>0.44</td><td>0.30</td><td>0.78</td><td>0.64</td><td>0.72</td><td>0.70</td><td>0.93</td><td>0.94</td><td>0.40</td><td>0.52</td><td>0.78</td><td>0.80</td></tr><tr><td>15</td><td>0.26</td><td>0.24</td><td>0.70</td><td>0.69</td><td>0.75</td><td>0.73</td><td>0.93</td><td>0.92</td><td>0.35</td><td>0.50</td><td>0.72</td><td>0.79</td></tr><tr><td>16</td><td>0.65</td><td>0.45</td><td>0.97</td><td>0.79</td><td>0.95</td><td>0.84</td><td>0.98</td><td>0.98</td><td>0.60</td><td>0.48</td><td>0.92</td><td>0.79</td></tr><tr><td>17</td><td>0.43</td><td>0.48</td><td>0.74</td><td>0.73</td><td>0.78</td><td>0.77</td><td>0.95</td><td>0.94</td><td>0.35</td><td>0.43</td><td>0.66</td><td>0.72</td></tr><tr><td>18</td><td>0.43</td><td>0.41</td><td>0.76</td><td>0.74</td><td>0.78</td><td>0.76</td><td>0.87</td><td>0.86</td><td>0.39</td><td>0.43</td><td>0.74</td><td>0.77</td></tr><tr><td>19</td><td>0.25</td><td>0.37</td><td>0.66</td><td>0.72</td><td>0.68</td><td>0.73</td><td>0.92</td><td>0.93</td><td>0.21</td><td>0.37</td><td>0.54</td><td>0.78</td></tr><tr><td>20</td><td>0.27</td><td>0.30</td><td>0.58</td><td>0.60</td><td>0.75</td><td>0.76</td><td>0.97</td><td>0.97</td><td>0.38</td><td>0.39</td><td>0.79</td><td>0.79</td></tr><tr><td>21</td><td>0.34</td><td>0.30</td><td>0.82</td><td>0.77</td><td>0.84</td><td>0.80</td><td>1.00</td><td>0.99</td><td>0.31</td><td>0.33</td><td>0.70</td><td>0.76</td></tr><tr><td>22</td><td>0.22</td><td>0.33</td><td>0.63</td><td>0.72</td><td>0.76</td><td>0.82</td><td>0.96</td><td>0.97</td><td>0.27</td><td>0.39</td><td>0.59</td><td>0.72</td></tr><tr><td>Mean</td><td>0.32</td><td>0.33**</td><td>0.72</td><td>0.71</td><td>0.78</td><td>0.77</td><td>0.95</td><td>0.94</td><td>0.34</td><td>0.41*,**</td><td>0.70</td><td>0.77*,**</td></tr><tr><td>S.D.</td><td>0.14</td><td>0.09***</td><td>0.13</td><td>0.08***</td><td>0.09</td><td>0.06***</td><td>0.07</td><td>0.06</td><td>0.10</td><td>0.08</td><td>0.09</td><td>0.05***</td></tr></tbody></table>

RP routine parameters, OP optimized parameters

up. This phantom study proposed and validated a strategy for optimizing PET scanners for multicentre cerebral studies while reducing, but not minimizing, variability across centres. It involved 22 PET centres with 13 different PET/CT scanner models, most of which were of recent generation. We ensured the reproducibility of phantom measurements by sending the same technologists to the PET centres for scanner set-up.

We checked the cross-calibration between the tomograph and the dose calibration device and found a difference of more than 10 % for four centres, which was immediately corrected. The alignment between the CT and PET scanners was also systematically checked and a correction was performed if necessary. Such a correction was needed in four additional centres.

## Choice of reconstruction parameters

Harmonization across scanners and centres for multicentre cerebral imaging trials was one of the achievements of a previous study by the ADNI [11]. For that study, which

<footnote>*Wilcoxon signed rank test, *p* < 0.05 for 22 centres; **Wilcoxon signed rank test, *p* < 0.05 for 18 centres; ***Pitman test, *p* < 0.05</footnote>
<figure><img src="image_7.png" /><figcaption>**Fig. 5** 3D Hoffman phantom results. Ratio values obtained with routine and optimized acquisition and reconstruction parameters in all centres. GM grey matter, WM white matter. *P* values represent the significant test results either for comparison of means (Wilcoxon test) or for comparison of standard deviations (Pitman test)</figcaption></figure>

included 50 centres and 17 different PET scanners, the PET centres were asked to acquire two 3D-Hoffman studies with recommended parameters. The ADNI quality-check team then checked the phantom images. For the analysis of the pooled images, a post-reconstruction smoothing filter, determined from phantom measurements, was applied to the images. This filter aimed at homogenizing the spatial resolution of the images across centres, and its application translated to a degradation of the resolution to the lowest one encountered [1].

In the present study, we chose to optimize the reconstruction parameters (with a product iterations × subsets superior to 50) and the post-reconstruction filter so that the recovery coefficients in the small cold and hot spheres would reach an optimized mean value and present limited dispersion around this optimal value. To this end, we reconstructed the images using a conventional 3D algorithm with a description of the statistics of the recorded data only, although PSF modelling reconstructions were available on the scanners that were of the more recent generations.

As expected, the reconstructions with PSF modelling provided recovery coefficients closer to 1 in the two smallest hot and cold spheres than the reconstructions without resolution modelling. However, Gibbs artefacts [12] were detected on the images at the edges of spherical objects. We therefore discarded them to avoid these artefacts and also to limit the discrepancies in RC in images reconstructed with and without resolution modelling across centres. On the other hand, we allowed algorithms with time-of-flight (TOF) modelling in two centres. It is known that state-of-the-art TOF values
<table><caption>**Table 3** Ratio values for both reconstructions and all centres obtained from Hoffman-3D acquisitions</caption><thead><tr><th rowspan="2">Centre number</th><th colspan="2">Right/left ratio</th><th colspan="2">Anterior/posterior ratio</th><th colspan="2">Putamen/WM ratio</th><th colspan="2">Caudate/WM ratio</th><th colspan="2">GM/WM ratio</th></tr><tr><th>RP</th><th>OP</th><th>RP</th><th>OP</th><th>RP</th><th>OP</th><th>RP</th><th>OP</th><th>RP</th><th>OP</th></tr></thead><tbody><tr><td>1</td><td>1.01</td><td>1.00</td><td>1.00</td><td>1.01</td><td>2.39</td><td>2.49</td><td>1.68</td><td>1.95</td><td>2.92</td><td>2.98</td></tr><tr><td>2</td><td>1.00</td><td>1.00</td><td>1.01</td><td>1.01</td><td>2.50</td><td>2.53</td><td>1.84</td><td>2.00</td><td>2.98</td><td>3.02</td></tr><tr><td>3</td><td>1.01</td><td>1.01</td><td>1.00</td><td>1.00</td><td>2.37</td><td>2.42</td><td>1.73</td><td>1.85</td><td>2.84</td><td>2.92</td></tr><tr><td>4</td><td>1.01</td><td>1.01</td><td>1.01</td><td>1.01</td><td>2.57</td><td>2.57</td><td>1.92</td><td>1.94</td><td>3.09</td><td>3.09</td></tr><tr><td>5</td><td>0.99</td><td>0.99</td><td>1.02</td><td>1.02</td><td>2.50</td><td>2.45</td><td>1.89</td><td>1.93</td><td>3.01</td><td>2.91</td></tr><tr><td>6</td><td>1.02</td><td>1.02</td><td>1.03</td><td>1.03</td><td>2.48</td><td>2.52</td><td>1.91</td><td>2.06</td><td>2.96</td><td>2.98</td></tr><tr><td>7</td><td>1.00</td><td>1.00</td><td>1.05</td><td>1.01</td><td>2.37</td><td>2.56</td><td>1.71</td><td>2.03</td><td>2.85</td><td>2.98</td></tr><tr><td>8</td><td>1.02</td><td>1.02</td><td>1.04</td><td>1.05</td><td>2.24</td><td>2.25</td><td>1.57</td><td>1.61</td><td>2.73</td><td>2.70</td></tr><tr><td>9</td><td>1.01</td><td>1.01</td><td>1.01</td><td>1.01</td><td>2.66</td><td>2.77</td><td>1.78</td><td>2.20</td><td>3.40</td><td>3.32</td></tr><tr><td>10</td><td>1.03</td><td>1.03</td><td>1.00</td><td>1.00</td><td>2.65</td><td>2.72</td><td>1.86</td><td>2.08</td><td>3.33</td><td>3.37</td></tr><tr><td>11</td><td>1.03</td><td>1.03</td><td>0.98</td><td>0.97</td><td>2.33</td><td>2.64</td><td>1.57</td><td>2.11</td><td>2.83</td><td>3.18</td></tr><tr><td>12</td><td>1.01</td><td>1.01</td><td>1.03</td><td>1.03</td><td>2.32</td><td>2.32</td><td>1.74</td><td>1.74</td><td>2.69</td><td>2.69</td></tr><tr><td>13</td><td>1.00</td><td>1.00</td><td>1.05</td><td>1.04</td><td>2.96</td><td>2.62</td><td>2.48</td><td>2.05</td><td>3.75</td><td>3.09</td></tr><tr><td>14</td><td>1.01</td><td>1.00</td><td>1.05</td><td>0.99</td><td>2.15</td><td>2.37</td><td>1.75</td><td>1.81</td><td>2.34</td><td>2.82</td></tr><tr><td>15</td><td>1.00</td><td>1.00</td><td>1.01</td><td>1.01</td><td>2.42</td><td>2.51</td><td>1.83</td><td>2.05</td><td>2.79</td><td>2.84</td></tr><tr><td>16</td><td>0.99</td><td>1.00</td><td>1.02</td><td>1.03</td><td>3.06</td><td>2.78</td><td>2.60</td><td>2.19</td><td>3.77</td><td>3.28</td></tr><tr><td>17</td><td>1.00</td><td>1.00</td><td>1.01</td><td>1.02</td><td>2.57</td><td>2.62</td><td>1.89</td><td>2.08</td><td>3.06</td><td>3.06</td></tr><tr><td>18</td><td>1.00</td><td>1.00</td><td>1.05</td><td>1.05</td><td>2.52</td><td>2.50</td><td>2.00</td><td>2.03</td><td>2.99</td><td>2.94</td></tr><tr><td>19</td><td>0.99</td><td>0.99</td><td>1.01</td><td>1.02</td><td>2.46</td><td>2.64</td><td>1.68</td><td>2.02</td><td>3.00</td><td>3.18</td></tr><tr><td>20</td><td>1.00</td><td>1.00</td><td>1.01</td><td>1.01</td><td>2.48</td><td>2.50</td><td>1.89</td><td>1.93</td><td>2.93</td><td>2.95</td></tr><tr><td>21</td><td>1.00</td><td>1.00</td><td>1.02</td><td>1.02</td><td>2.58</td><td>2.51</td><td>1.95</td><td>1.93</td><td>3.06</td><td>2.93</td></tr><tr><td>22</td><td>1.03</td><td>1.02</td><td>1.00</td><td>1.01</td><td>2.57</td><td>2.67</td><td>1.82</td><td>2.04</td><td>3.23</td><td>3.28</td></tr><tr><td>Mean</td><td>1.01</td><td>1.01</td><td>1.02</td><td>1.02</td><td>2.51</td><td>2.54<sup>***</sup></td><td>1.87</td><td>1.98<sup>***</sup></td><td>3.03</td><td>3.02</td></tr><tr><td>S.D.</td><td>0.01</td><td>0.01</td><td>0.02</td><td>0.02</td><td>0.20</td><td>0.13<sup>***</sup></td><td>0.24</td><td>0.14<sup>***</sup></td><td>0.32</td><td>0.18<sup>***</sup></td></tr></tbody></table>

RP routine parameters, OP optimized parameters, WM white matter, GM grey matter

*Wilcoxon signed rank test, *p* < 0.05 for all centres; **Wilcoxon signed rank test, *p* < 0.05 for 18 centres; ***Pitman test, *p* < 0.05

have an important effect on signal-to-noise ratio of whole body imaging, but that effect is negligible on smaller objects such as brain structures [13].

Conversely, in images where both spatial resolution and RC were too low, we chose to use more iterations of the algorithm in order to enhance the spatial resolution of the images and to apply a Gaussian (FWHM between 2 and 4 mm) post-reconstruction smoothing filter to the images. The pixel spacing was between 1 and 3 mm in all optimized images.

## Improving contrast recovery and dispersion of RC values

With optimized parameters, the RC significantly improved for the cold spheres, but not for the hot spheres, of close diameter. That difference between cold and hot spheres is partly related to the presence of the sphere walls, which are intrinsically cold. These walls affect the quantification to a greater extent in hot spheres than in cold spheres. Such a cold wall is specific to the phantom. One should also note that the optimized RC was higher in the hot spheres than in the cold spheres of similar diameter. The
quantification in cold objects is complex, depending not only on spatial resolution but also on scatter correction and spatial sampling [14, 15], and of the non-negativity constraint of the statistical reconstruction algorithm MLEM without a specific description [16].

We also significantly reduced the variability of RC in four of the six spheres of the Jaszczak phantom. As shown in Figs. 3 and 4, this reduction of variability was mostly due to the suppression of outliers by discarding PSF modelling reconstruction algorithms.

It should be emphasized that the RC measurements were obtained with in-house software, which was developed with great care to ensure reproducibility and precision. In particular, the VOIs were not drawn directly on the PET images but were automatically segmented on CT images and optimized with knowledge from the phantom specifications.

## Improving quantification and variability in small brain structures and image uniformity

We checked the quantification in the striatum of the Hoffman brain phantom and observed an increase in activity concentration measured in the putamen and caudate, when compared to non-optimized images. We also significantly reduced after the optimization the dispersion of values for putamen, caudate, and grey-matter-to-white-matter ratios.

Finally, the good uniformity of the images across the field of view was confirmed by the ratio of right-to-left activity concentration and anterior-to-posterior activity concentration measured on the Hoffman brain images. Additionally, as expected, the choice of the reconstruction mode parameters had no impact on the ratio of activity concentration in the right-to-left and frontal-to-occipital regions of interest. Those ratios should not be affected by the spatial resolution. However, they are likely affected by low frequency variation resulting from scatter, uniformity, and attenuation corrections. Attempts to take into account these effects before pooling data sets from multiple centres were suggested by Joshi et al. [1], but the results were not considered as convincing by the authors. In the present work, the compensation for scatter, random coincidences, and attenuation was performed using the manufacturer's latest techniques, with slight differences in the implementation for the different scanner models. The Hoffman phantom also enabled us to verify that the quantification in typical small structures of the brain, such as the caudate, was improved after the reconstruction optimization.

## Image spatial resolution and signal-to-noise ratio

On average, image spatial resolution was slightly improved, despite a marked decrease for the centres equipped with reconstruction algorithms including TOF and PSF modelling. This improvement is likely due to the increase in the number of iterations in the other centres. As expected, this improvement is counterbalanced by a decrease of the signal-to-noise ratio. We checked with Hoffman phantom acquisitions that this reduction did not impact the visual aspect of the images.

## Variability across centres

The variability across centres was reduced but not minimized by this work. For that, it would have been necessary to apply post-filtering tailored per site according to ISR. We
deliberately chose not to do this in order to preserve the spatial resolution. This latter
goal was achieved for the 18 centres without PSF modelling reconstruction.

In order to account for the residual variability across centres, we are currently design-
ing a statistical model that takes into account the centre effect (as assessed by the
phantom studies). Several other parameters have been shown to influence the variabil-
ity of the measurements across centres, including injected dose per kilogram, delay be-
tween the injection and the scan acquisition, patient positioning, weight, and
glycaemia, and these have been extensively discussed by Boellard [17, 18]. We will add
to the statistical model these additional acquisition and patient-related parameters (as
assessed by an ancillary monocentric study involving 300 subjects), and we will test it
on the MEMENTO cohort.

If the statistical model happens not to be sufficiently effective, our approach could be
updated in a second phase according to the strategy proposed by Lasnon et al. [19].
The latest advances in reconstruction schemes could be used to reconstruct images to-
gether with an adequate post-reconstruction filtering, offering the best compromise in
image quantification for multicentre studies.

Conclusions

This work was undertaken by CATI's team in charge of harmonization, quality check, and analysis of multicentre brain PET or SPECT acquisitions. The proposed procedure for PET imaging optimization enabled the production of images with a more homogenous spatial resolution across the centres. The spatial resolution was also preserved thanks to the use of statistical image reconstructions with iterations x subsets of at least 50, a Gaussian smoothing post-reconstruction filter with a FWHM between 2 and 4 mm, and a pixel spacing between 1 and 3 mm. The uniformity across the fields of view of the scanners was good.

We provided recommendations to each centre for minimizing the influence of other
factors such as injected dose and interval post-injection. In addition, centres were
qualified after the analysis of the images of a first test patient. CATI has now received
more than 3000 PET images.

This step of optimization of the image characteristics acquired at different centres
will thus allow us to account for the residual variability, which could be handled in the
final statistical analysis using adequate modelling.

Ethical approval

This article does not contain any studies with human participants or with animals per-
formed by any of the authors.

Competing interests

MO Habert received honoraria from GE healthcare as a consultant and from PIRAMAL and ROCHE as a speaker. S
Marie has no potential conflict of interest to declare. H Bertin has no potential conflict of interest to declare. JB Martini
has no potential conflict of interest to declare. M Reynal has no potential conflict of interest to declare. M Diallo has
no potential conflict of interest to declare. A Kas received honoraria from PIRAMAL and ROCHE as a speaker. R Trebossen
has no potential conflict of interest to declare.

Authors' contributions

MOH participated in the study design and coordination, and wrote the draft. SM supervised fantom acquisitions and
participated in the analyses. HB participated in the data phantom analyses and carried out figures and tables of the
manuscript. MR and JBM developed in-house software for the study. MD carried out the computation of the spatial
image resolution. AK participated in the study design. The CATI Nuclear medicine network contributed to data
acquisition. RT conceived of the study and participated in the study design. All authors read and approved the final manuscript.

### Acknowledgements

We sincerely thank the staff from all departments of Nuclear Medicine who participated in the MEMENTO cohort study for their availability and their patience in spite of our numerous requests.

We also are grateful to Alain Giron and Marion Houot for their help with statistical analyses, and to the staff from Esprimed (www.esprimed.net), who were in charge of phantom acquisitions, for their efficiency and professionalism. Finally, we greatly appreciate the support of Philips, GE Healthcare, and Siemens engineers in determining and setting up the optimized reconstruction parameters.

### Funding

The French Foundation Plan Alzheimer supported this work.

### The CATI Nuclear Medicine network includes as follows:

Pascal Bailly (CHU Sud, Amiens), Sébastien Balduyck (Hôpital Purpan, Toulouse), Francis Bouchet (CHU Angers), Hatem Boulahdour (Hôpital Jean Minjoz, Besançon), David Bourhis (CHRU Morvan, Brest), Thomas Carlier (Hôpital Nord Laennec, Nantes), Jacques Darcourt (Centre Antoine Lacassagne, Nice), Marielle Decousus (Hôpital Nord, Saint-Etienne), Delphine De Verbier Lonjon (Hôpital Gui de Chauliac, Montpellier), Inna Dygai-Cochet (Centre George François Leclerc, Dijon), Élise Enderlin (Hôpital de Hautepierre, Strasbourg), Karim Farid (Centre Cardiologique du Nord, Saint-Denis), Bardia Farman (Hôpital La Timone, Marseille), Philippe Fernandez (Hôpital Xavier Arnozan, Pessac), Pacôme Fosse (CHU Angers), Marjolaine Fourcade (Hôpital Gui de Chauliac, Montpellier), Christelle Gallais (CHU Poitiers - TEP GIE), Céline Gallazzini-Crepin (Hôpital Albert Michallon, Grenoble), Éric Guedj (Hôpital La Timone, Marseille), Sébastien Hapdey (Centre Henri Becquerel, Rouen), Anne Hitzel (Hôpital Purpan, Toulouse), Claude Hossein-Foucher (Hôpital Huriez, Lille), Laetitia Imbert (CHU Nancy-Brabois), Anthony Kelly (Centre Jean Perrin, Clermont-Ferrand), Pierre Malick Koulibaly (Centre Antoine Lacassagne, Nice), Fabien Maurel (CHU Nice), Marc-Etienne Meyer (CHU Sud, Amiens), Corinne Millardet (Centre Jean Perrin, Clermont-Ferrand), Izzy Jacques Namer (Hôpital de Hautepierre, Strasbourg), Amandine Pallardy (Hôpital Nord Laennec, Nantes), Catherine Pérault (Hôpital Bretonneau, Tours), Rémy Perdrisot (CHU Poitiers - TEP GIE), Yolande Petegnief (Hôpital Jean Minjoz, Besançon), Gabriel Pop (Hôpital Avicenne, Bobigny), Nathalie Prevot-Bitot (Hôpital Nord, Saint-Etienne), Mathieu Queneau (Centre Cardiologique du Nord, Saint-Denis), Solenne Querellou (CHRU Morvan, Brest), Micheline Razzouk-Cadet (CHU Nice), Ghislaine Reboulet (Hôpital Albert Michallon, Grenoble), Maria-Joao Santiago-Ribeiro (Hôpital Bretonneau, Tours), Pierre Vera (Centre Henri Becquerel, Rouen), Antoine Verger (CHU Nancy-Brabois), Maximilien Vermandel (Hôpital Huriez, Lille), Jean-Marc Vrigneaud (Centre George François Leclerc, Dijon), Marie-Joséphine Waryn (Hôpital Avicenne, Bobigny), and Michel Zanca (Hôpital Gui de Chauliac, Montpellier).

### Author details

¹Sorbonne Universités, UPMC Univ Paris 06, CNRS, INSERM, Laboratoire d'Imagerie Biomédicale, F-75013 Paris, France.

²Centre pour l'Acquisition et le Traitement des Images (www.cati-neuroimaging.com), Paris, France. ³AP-HP, Hôpital Pitié-Salpétrière, Département de Médecine Nucléaire, F-75013 Paris, France. ⁴CEA-Service Hospitalier Frédéric Joliot, F-91401 Orsay, France. ⁵Faculté de Médecine - Laboratoire d'Imagerie Biomédicale, 91 Boulevard de l'Hôpital, F-75013 Paris, France.

Received: 11 November 2015 Accepted: 9 February 2016

Published online: 05 April 2016

### References

1. Joshi A, Koeppé RA, Fessler JA. Reducing between scanner differences in multi-center PET studies. NeurolImage. 2009;46(1):154-9.

2. Herholz K, Salmon E, Perani D, Baron C, Holthoff V, Frölich L, Schönknecht P, Ito K, Mielke R, and Kalbe E. Discrimination between Alzheimer dementia and controls by automated analysis of multicenter FDG PET. NeurolImage. 2002;17(1):302-16.

3. Mosconi L, Tsui WH, Herholz K, Pupi A, Drzezga A, Lucignani G, Reiman EM, Holthoff V, Kalbe E, Sorbi S, Diehl-Schmid J, Perneczky R, Clerici F, Caselli R, Beuthien-Baumann B, Kurz A, Minoshima S, and de Leon MJ. Multicenter standardized 18 F-FDG PET diagnosis of mild cognitive impairment, Alzheimer's disease, and other dementias. J Nucl Med. 2008;49(3):390-8.

4. Geworski L, Knoop BO, de Wit M, Ivancević V, Bares R, Munz DL. Multicenter comparison of calibration and cross calibration of PET scanners. J Nucl Med. 2002;43(5):635-9.

5. Westerterp M, Pruim J, Oyen W, Hoekstra O, Paans A, Visser E, van Lanschot J, Sloof G, and Boellaard R. Quantification of FDG PET studies using standardised uptake values in multi-centre trials: effects of image reconstruction, resolution and ROI definition parameters. Eur J Nucl Med Mol Imaging. 2007;34(3):392-404.

6. Scheuermann JS, Saffer JR, Karp JS, Levering AM, Siegel BA. Qualification of PET scanners for use in multicenter cancer clinical trials: the American College of Radiology Imaging Network experience. J Nucl Med. 2009;50(7): 1187-93.

7. Varrone A, Asenbaum S, Vander Borght T, Booi J, Nobilli F, Någren K, Darcourt J, Kapucu OL, Tatsch K, Bartenstein P, Van Laere K, and European Association of Nuclear Medicine Neuroimaging Committee. EANM procedure guidelines for PET brain imaging using [18 F]FDG, version 2. Eur J Nucl Med Mol Imaging. 2009;36(12):2103-10.

8. Ballard DH. Generalizing the Hough transform to detect arbitrary shapes. Pattern Recogn. 1981;13(2):111-22.

9. Prieto E, Martí-Climent JM, Arbizu J, Garrastachu P, Domínguez I, Quincoces G, García-Velloso MJ, Lecumberri P, Gómez-Fernández M, Richter JA. Evaluation of spatial resolution of a PET scanner through the simulation and experimental measurement of the recovery coefficient. Comput Biol Med. 2010;40(1):75-80.
10. Pitman EJG. A note on normal correlation. Biometrika. 2013;31:9-12.

11. Weiner MW, Veitch DP, Aisen PS, Beckett LA, Cairns NJ, Green RC, Harvey D, Jack CR, Jagust W, Liu E, Morris JC, Petersen RC, Saykin AJ, Schmidt ME, Shaw L, Shen L, Siuciak JA, Soares H, Toga AW, Trojanowski JQ, and Alzheimer's Disease Neuroimaging Initiative. The Alzheimer's Disease Neuroimaging Initiative: a review of papers published since its inception. Alzheimers Dement. 2013;9(5):e111-94.

12. Zeng GL. Gibbs artifact reduction by nonnegativity constraint. J Nucl Med Technol. 2011;39(3):213-9.

13. Conti M. Focus on time-of-flight PET: the benefits of improved time resolution. Eur J Nucl Med Mol Imaging. 2011; 38:1147-57.

14. Stute S, Comtat C. Practical considerations for image-based PSF and blobs reconstruction in PET. Phys Med Biol. 2013;58(11):3849-70.

15. Bendriem B, Dewey SL, Schlyer DJ, Wolf AP, Volkow ND. Quantitation of the human basal ganglia with positron emission tomography: a phantom study of the effect of contrast and axial positioning. IEEE Trans Med Imaging. 1991;10(2):216-22.

16. Van Slambrouck K, Stute S, Comtat C, Sibomana M, van Velden F, Boellaard R, and Nuyts J. Bias reduction for low-statistics PET: maximum likelihood reconstruction with a modified Poisson distribution. IEEE Trans Med Imaging. 2015;34(1):126-36.

17. Boellaard R, O'Doherty MJ, Weber WA, Mottaghy FM, Lonsdale MN, Stroobants SG, Oyen WJ, Kotzerke J, Hoekstra OS, Pruim J, Marsden PK, Tatsch K, Hoekstra CJ, Visser EP, Arends B, Verzijlbergen FJ, Zijlstra JM, Comans EF, Lammertsma AA, Paans AM, Willemse AT, Beyer T, Bockisch A, Schaefer-Prokop C, Delbeke D, Baum RP, Chiti A, and Krause BJ. FDG PET and PET/CT: EANM procedure guidelines for tumour PET imaging: version 1.0. Eur J Nucl Med Mol Imaging. 2010; 37(1):181-200.

18. Boellaard R. Need for standardization of 18 F-FDG PET/CT for treatment response assessments. J Nucl Med. 2011; 52 Suppl 2:93S-100S.

19. Lasnon C, Desmonts C, Quak E, Gervais R, Do P, Dubos-Arvis C, and Aide N. Harmonizing SUVs in multicentre trials when using different generation PET systems: prospective validation in non-small cell lung cancer patients. Eur J Nucl Med Mol Imaging. 2013;40(7):985-96.

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

► Convenient online submission

► Rigorous peer review

► Immediate publication on acceptance

► Open access: articles freely available online

► High visibility within the field

► Retaining the copyright to your article