# Normal Dual Isotope V/Q SPECT Model for Monte-Carlo Studies

David Bourhis, Laura Wagner, Marine Essayan, Philippe Robin, Romain Le Pennec, Pierre Yves Salaun, Pierre Yves Le Roux

► To cite this version:

David Bourhis, Laura Wagner, Marine Essayan, Philippe Robin, Romain Le Pennec, et al. Normal Dual Isotope V/Q SPECT Model for Monte-Carlo Studies. Frontiers in Medicine, 2020, 7, pp.461. 10.3389/fmed.2020.00461 . hal-05413497

HAL Id: hal-05413497

https://hal.univ-brest.fr/hal-05413497v1

Submitted on 12 Dec 2025

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
<figure><img src="image_2.png" /></figure>

Check for
updates

# Normal Dual Isotope V/Q SPECT Model for Monte-Carlo Studies

David Bourhis<fnref n="1" />,²*, Laura Wagner<fnref n="1" />, Marine Essayan<fnref n="1" />, Philippe Robin<fnref n="1" />,², Romain Le Pennec<fnref n="1" />, Pierre Yves Salaun<fnref n="1" />,² and Pierre Yves Le Roux<fnref n="1" />,²

**Background:** There is currently no reliable or validated tool to delineate and quantify functional lung volumes with ventilation/perfusion (V/Q) SPECT/CT. The main challenges encountered include the physiological non-uniformity of lung function, such as the anterior-to-posterior gradient on perfusion images, and the lack of ground truth to assess the accuracy of delineation algorithms. In that respect, Monte-Carlo simulations would be an interesting tool. Thus, the aim of this study was to develop a realistic model of dual-isotope lung V/Q SPECT-CT Monte-Carlo simulations, integrating the anterior to posterior gradient on perfusion.

**Methods:** Acquisitions and simulations parameters were set in accordance to nuclear medicine guidelines for V/Q lung SPECT-CT. Projections were acquired and simulated, then the reconstructions [with and without attenuation correction (AC)] were compared. A model was built from a patient's CT scan. To model the anterior to posterior gradient, the lungs were divided into sixteen coronal planes, where a rising radioactivity concentration was set. To assess the realism of simulations, they were compared to a normal co-registered normal cases database in terms of pixelwize Z-score map.

**Results:** For ventilation images, mean (SD) Zscores on Zscore maps were -0.2 (0.7) and -0.2 (0.7) for AC and noAC images, respectively. For perfusion images, mean (SD) Zscores were -0.2 (0.6) and -0.1 (0.6) for AC and noAC images, respectively.

**Conclusion:** We developed a model for dual isotopes lung V/Q SPECT-CT, integrating the anterior-to-posterior gradient on perfusion images. This model could be used to build a catalog of clinical scenarios, in order to test delineation methods of functional lung volumes.

**Keywords:** V/Q SPECT, lung function, PE, simulation, Z-score

## INTRODUCTION

Lung function evaluation mainly relies on pulmonary function tests, which provide information about global lung function. Lung ventilation/perfusion scintigraphy is an imaging modality that provides complementary information about the regional distribution of lung function. Ventilation is assessed using inert gases or radio-labeled aerosols while local pulmonary blood flow is assessed after administration of ⁹⁹mTc labeled albumin macro-aggregates (MAA) (1, 2). Quantification of lung functional volumes demonstrated a potential utility in various clinical scenarios. This includes the quantification of the pulmonary vascular obstruction index (PVOI) in patients with pulmonary

**OPEN ACCESS**

**Edited by:**

Thomas Carlier,
Centre Hospitalier Universitaire (CHU)
de Nantes, France

**Reviewed by:**

Dale L. Bailey,
The University of Sydney, Australia
Robin De Nijs,
Copenhagen University
Hospital, Denmark

**\*Correspondence:**

David Bourhis
david.bourhis@chu-brest.fr

**Specialty section:**

This article was submitted to
Nuclear Medicine,
a section of the journal
Frontiers in Medicine

**Received:** 03 June 2020

**Accepted:** 10 July 2020

**Published:** 18 August 2020

**Citation:**

Bourhis D, Wagner L, Essayan M, Robin P, Le Pennec R, Salaun PY and Le Roux PY (2020) Normal Dual Isotope V/Q SPECT Model for Monte-Carlo Studies. *Front. Med.* 7:461. doi: 10.3389/fmed.2020.00461

<footnote n="1">¹ Service de Médecine Nucléaire, Centre Hospitalier Régional Universitaire de Brest, Brest, France, ² EA3878 GETBO, Université de Bretagne Occidentale, Brest, France</footnote>
embolism (PE), which has been shown to be predictive of PE recurrence, either at diagnosis or after anticoagulant therapy (3–5). Another indication is radiation oncology (6). The last 20 years have seen the emergence of new methods of irradiation (IMRT, VMAT, stereostatic radiotherapy), increasing the possibility of individualizing radiotherapy treatment plans, in order to increase the intensity to the tumor and/or to preserve lung function and limit the risk of pneumonitis (7, 8).

Although lung scintigraphy is a very attractive tool for regional lung function assessment, the main factor limiting a wider use of the test in clinical practice is the lack of reliable and validated tools to delineate and quantify functional lung volumes. Functional lung delineation is not trivial because of the physiological non-uniformity of lung ventilation and blood flow, mainly because of anatomy and gravity (9), leading to an anterior-to-posterior intensity gradient on perfusion SPECT. This makes delineation methods based on pixel intensities inaccurate, as an intensity value could be normal in the ventral area or pathological in the dorsal area, even if a quantified reconstruction is used. Several original segmentation methods have been proposed (10–13), but studies are often limited by the lack of ground truth.

Monte-Carlo simulations would be an interesting tool for assessing and validating a delineation algorithm for lung functional volumes with V/Q SPECT/CT. First, the advantage of simulating realistic lung V/Q SPECT-CT would be to know the ground truth for the evaluation and comparison of delineation methods (e.g., the exact size of a perfusion defect for PE, the exact pulmonary vascular obstruction index). Second, it would allow to test a much greater number of clinical scenarios (e.g., PE of various size and location). Finally, simulations may integrate physiological parameters, such as gravity, which make conventional delineation methods inaccurate in some situations.

Monte-Carlo simulations have been widely proposed in SPECT for multiple objectives (14). However, only a few studies have been published in the area of lung V/Q SPECT (12, 13, 15–17), and modeled the lung segmentation (18–22), but none have modeled the anterior to posterior gradient on perfusion images. Furthermore, no study has modeled a dual energy 81mKr gas ventilation associated with 99mTc-MAA perfusion SPECT.

The aim of this study was to develop a realistic model of dual-isotopes lung V/Q SPECT-CT.

# MATERIALS AND METHODS

## Data Acquisition and Simulation Parameters

All acquisitions were performed on a Symbia T6 system (Siemens, Erlangen) equipped with a medium energy low penetration (MELP) collimator, as it is used for dual-isotopes lung V/Q SPECT-CT (1). SPECT were acquired with 4.7 mm pixels, 128 projections, 10 s/projections, step&shoot mode, non-circular orbit, in accordance with nuclear medicine clinical guidelines (1, 23). Energy windows for 99mTc scatter, 99mTc photopeak, 81mKr scatter, and 81mKr photopeak were [109.9, 129.5], [129.5, 150.5], [150.7, 177.6], [177.6, 206.4],

respectively. CT acquisitions parameters were 110kV, 16 mAs with automatic exposure, pitch1. Attenuation correction CT (ACCT) reconstruction parameters were 512² matrix, 0.98 mm pixel, 5 mm slice width and B08 filter. SPECT reconstructions were performed on Siemens MI-Apps software with FLASH3D (OSEM3D with collimator detector response modeling), 4 iterations, 8 subsets, 8.4 mm gaussian post-filtering, scatter correction [dual energy window method (DEW)], with and without attenuation correction.

The Monte-Carlo package used to run SPECT simulations was Simind (6.0) (24). The camera modeling parameters were set in order to correspond to a Symbia T6 system (Siemens, Erlangen) (25) equipped with a MELP collimator. All measurements were first acquired on a SymbiaT6 gamma camera, then simulated with the same parameters on Simind. ACCT reconstruction was used to define the simulation geometries. Simulation geometries were Zubal-like phantoms (26, 27), built from the ACCT. ACCT data was segmented according to hounsfield units using MiM software (7.0, Cleveland). Images bit depth was set to 8 bits, and a unique value was attributed to each segmented areas using ImageJ software (28). Those values were used in Simind to set the desired value of density and radioactivity concentration in the defined areas. SPECT were simulated with the same parameters than the acquisitions. Photons emitted from 99mTc decay were simulated with a 140 keV Energy and 88.5% abundance. As Krypton gas is continuously inspired and expired and has a very fast decay (half-life is 13 s), it was simulated as a stationary gas without significant decay, with homogeneous concentration, with a 190 keV energy and 100% abundance.

All SPECT reconstructions were performed on Siemens MI-Apps software, with FLASH3D, four iterations, eight subsets, 8.4 mm gaussian post-filtering, scatter correction (DEW), with and without attenuation correction.

In order to assess the simulated data, image quality tests were run on both simulated and acquired data, and compared in terms of Root Mean Square Deviation (RMSD) for spatial resolution and activity recovery (see Appendix 1).

## Lung V/Q SPECT Realistic Model
### Simulation Parameters

Using the ACCT, five representative tissues were delineated using MiM Software (MiM 7.0, Cleveland), including outside air, bones, fat, soft tissues, lungs, and bronchi. A code was assigned to each area. Digital phantom was sub-sampled in a 128² matrix to accelerate the simulation calculation, and the simulation grid was 128 × 128 × 108 matrix, corresponding to 3.92 × 3.92 × 3.59 mm voxels. With regards to ventilation, the simulated radioactivity was evenly set to 55 kBq·mL⁻¹ in the lungs and the airways. With regards to perfusion, in order to model anterior to posterior gradient of the distribution of radioactivity on perfusion SPECT, lungs were divided into sixteen coronal planes. For each coronal plane a relative to maximum radioactivity concentration value was assigned. The values were set to fit the mean anterior to posterior intensity gradient measured on a normal cases database (29), from 43 kBq·mL⁻¹ inside the first coronal plane to 65 kBq·mL⁻¹ in the last coronal plane. As the source map is different when simulating ventilation and perfusion, simulations were
not run simultaneously. Scatter data was stored at each energy window, and 81mKr scatter was added to 99mTc lower scatter and primary energy windows.

## Reconstructions Analysis and Z-Score Maps

In order to assess the realism and consistency of simulated SPECT, reconstructed simulations were compared to a 73 co-registered normal cases database (29). Z-score maps were generated as follows: First, SPECT simulated reconstructions were normalized according to the mean value. Z-score map was then calculated at each voxel coordinate as follows:

$$Z - score_{SPECT Reconstruction}(x, y, z) = \frac{[Pixelvalue(x, y, z) - MEANmapValue(x, y, z)]}{SDmapValue(x, y, z)};$$

Z score histograms were computed in a ROI built 0.5 cm inside the lungs CT boundaries to avoid partial volume effect artifacts. Simulations were evaluated as the mean Zscore ± SD, and as the pourcentage of voxels below or above 1 SD.

## Example of Application

As an example of the possible applications, a simulation with a segmental perfusion defect was performed. A segmental area was

delineated on ACCT and the relative radioactivity concentration was set to 50% in this zone.

## RESULTS

### Realistic Simulated Lung V/Q SPECT

Digital phantom, radioactivity source map and the corresponding AC and NoAC perfusion reconstructions are shown in Figure 1. An axial slice and the corresponding Zscore maps and 3D histogram for ventilation, perfusion, AC and NoAC reconstructions are shown in Figure 2. Zscores distribution for all reconstructions are summarized in Table 1. Mean Zscores ranged from -0.1 to -0.2. An example of PE simulation, its radioactive source map and the corresponding perfusion AC reconstruction is shown in Figure 3.

## DISCUSSION

In this study we created a realistic model for dual isotopes V/Q SPECT-CT Monte-Carlo simulations. The gamma camera was correctly modeled, as shown by spatial resolution and contrast recovery measurements. Reconstructed simulations of ventilation and perfusion scans were very close to a normal cases database including 73 normal co-registered V/Q SPECT/CT (29).

<figure><img src="image_3.png" /><figcaption>FIGURE 1 | Digital phantom, radioactivity source map definition, AC and NoAC reconstructions of the simulated projections.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | AC and NoAC reconstructions for ventilation and perfusion (A), corresponding Z-score map (B) and pixels Z-score histograms (C), and statistics measured 0.5 cm inside CT boundaries.</figcaption></figure>

<table><caption>**TABLE 1** | Zscores distribution for ventilation AC, NoAC, and perfusion AC and NoAC.</caption><thead><tr><th></th><th>Ventilation<br>AC</th><th>Ventilation<br>NoAC</th><th>Perfusion<br>AC</th><th>Perfusion<br>NoAC</th></tr></thead><tbody><tr><td>Mean Zscore (SD)</td><td>-0.2 (0.7)</td><td>-0.2 (0.7)</td><td>-0.2 (0.6)</td><td>-0.1 (0.6)</td></tr><tr><td>%pixels > 1</td><td>4.6</td><td>5.2</td><td>3.5</td><td>4.9</td></tr><tr><td>%pixels < -1</td><td>1.9</td><td>1.3</td><td>1.9</td><td>3.8</td></tr></tbody></table>

The Monte-Carlo code used in this study was simind. As compared with Gate (30) or MCNPX (31) packages, this program is known to be fast, thanks to variance reduction techniques and analytic calculation of the collimator response. Taheri et al. showed that the accuracy was comparable to other packages (32). The use of a user friendly graphical user interface, CHANGE, and the ability to write scripts to run automatically simulations with different parameters makes Simind an easy package to run multiple SPECT simulations.

The system performances comparison showed slight differences between simulations and acquisitions. The root mean square deviation was lower than 1 mm for spatial resolution and lower than 1% for contrast recovery. Thus, the system was correctly modeled and usable to run realistic simulations, in

accordance with published studies (33, 34). Energy resolution was not verified, as it is an input of the simulation. Sensitivity was not verified either, as images are scaled before adding Poisson noise. This may have been of concern if the simulations were run in a dosimetric purpose, but not for producing images. In order to model the anterior to posterior intensity gradient, we defined 16 horizontal planes in the patient orientation, so that each plane's width was close to 1 cm (3 pixels). This is far under the spatial resolution for SPECT with a 30 cm radius, and thus is enough to model precisely the gradient. The radioactivity values in each plane was defined based on a parametric mean perfusion AC map obtained from the normal cases database (29). Based on Z-score analysis, V/Q SPECT modeling was satisfactory, with a mean Z value close to zero, and at least 91.3% of the pixels ranging from -1 to 1 SD. We measured Z scores in ROI defined 0.5 cm inside CT boundaries corresponding to a 1-pixel width in order to avoid false measurements due to partial volume effects on the edges of the Zscore map.

A current challenge in PE management is the ability to estimate the Pulmonary Vascular Obstruction Index (PVOI), which has been shown to be predictive of PE recurrence (3-5). In order to measure precisely the PVOI with V/Q SPECT, several methods have been proposed (11-13, 35), but an accurate evaluation remains limited by the lack of ground truth. We illustrated the interest of a realistic model of lung V/Q
<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Application exemple: PE simulation. Radioactivity source map representing the ground Truth (A) and the corresponding perfusion AC reconstruction (B).</figcaption></figure>

SPECT-CT Monte-Carlo simulations in **Figure 3**. A knowledge of the real volume and geometry of perfusion defects, for various morphologies and anterior to posterior gradient intensities, should allow to develop and test new delineation methods, such as statistical map threshold.

Our study has some limitations. First, we used a static Zubal-like phantom instead of dynamic X-CAT phantom. Although respiratory motion is a challenge in PET imaging, it has less impact on SPECT images because of the lower spatial resolution of the technique. Indeed, based on the analysis of 73 normal co-registrated cases database (29), the impact of respiratory motion on the uptake variability in the basis areas was very low, especially as compared with the anterior to posterior gradient. Second, the DEW scatter correction method is not the more efficient, since both scatter (from ⁹⁹mTc itself) and downscatter (from ⁸¹mKr) are present in the main Tc-window around 140 keV. It has been shown that triple energy windows method (TEW) is better in the case of multiple energy isotopes (36), or here in multiple isotopes SPECT. However, given that the normal cases database was created with the DEW method, the reconstructions of the simulations were performed with the same method to avoid a bias in Z-score analysis. Third, we did not model the mild gradient described by some authors in the caudo-cranial direction (35). Similarly, this gradient was negligible in the 73

normal co-registrated cases database (29). Fourth, some more singularities have been described, such as the fissure sign or the segmental contour sign (37). The distribution of ⁹⁹mTc-MAA often does not extend to the periphery of a segment or lobe, often finishing 1 cm or less before the anatomical boundary. This is usually attributed to the lack of pulmonary artery supplied perfusion to the peripheral surface of the lung which is supplied by the bronchial circulation. This was not simulated in our model. Finally, Walrand et al. showed that it was possible to simulate regional heterogeneity of liver perfusion with ⁹⁰Y-microspheres, taking into account the diameter of arteries (38). Similarly, simulating the physics of MAA particles inside the pulmonary arterial tree may have improved the realism of simulations. However, these simulations were not developed for dosimetry purposes.

## CONCLUSION

We developed a realistic model for dual isotopes lung V/Q SPECT-CT, integrating the anterior posterior gradient on perfusion images. This model can be used to build a catalog of clinical scenarios, in order to test delineation methods of functional lung volumes. In the context of PE, this could help to develop new delineation methods for PVOI estimation.
# DATA AVAILABILITY STATEMENT

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

# AUTHOR CONTRIBUTIONS

DB, PL, and PS contributed to designing the study. DB, LW, ME, PR, RL, PL, and PS contributed to managing the imaging

procedures. DB, LW, ME, PL, and PS contributed to analyzing the data. All authors contributed to writing the manuscript, read, and approved the final manuscript.

# SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2020.00461/full#supplementary-material

# REFERENCES

1. Bajc M, Schümichen C, Grüning T, Lindqvist A, Le Roux PY, Alatri A, et al. EANM guideline for ventilation/perfusion single-photon emission computed tomography (SPECT) for diagnosis of pulmonary embolism and beyond. *Eur J Nucl Med Mol Imaging*. (2019) 46:2429–51. doi: 10.1007/s00259-019-04450-0

2. Le Roux PY, Robin P, Salaun PY. New developments and future challenges of nuclear medicine and molecular imaging for pulmonary embolism. *Thromb Res*. (2018) 163:236–41. doi: 10.1016/j.thromres.2017.06.031

3. Pesavento R, Filippi L, Palla A, Visonà A, Bova C, Marzolo M, et al. Impact of residual pulmonary obstruction on the long-term outcome of patients with pulmonary embolism. *Eur Respir J*. (2017) 49:1601980. doi: 10.1183/13993003.01980-2016

4. Planquette B, Ferré A, Peron J, Vial-Dupuy A, Pastre J, Mourin G, et al. Residual pulmonary vascular obstruction and recurrence after acute pulmonary embolism. A single center cohort study. *Thromb Res*. (2016) 148:70–5. doi: 10.1016/j.thromres.2016.10.030

5. Tromeur C, Sanchez O, Presles E, Pernod G, Bertoletti L, Jego P, et al. Risk factors for recurrent venous thromboembolism after unprovoked pulmonary embolism: the PADIS-PE randomised trial. *Eur Respir J*. (2018) 51:1701202. doi: 10.1183/13993003.01202-2017

6. De Jaeger K, Seppenwoolde Y, Boersma LJ, Muller SH, Baas P, Belderbos JS, et al. Pulmonary function following high-dose radiotherapy of non-small-cell lung cancer. *Int J Radiat Oncol Biol Phys*. (2003) 55:1331–40. doi: 10.1016/S0360-3016(02)04389-4

7. De Bari B, Filippi AR, Mazzola R, Bonomo P, Trovò M, Livi L, et al. Available evidence on re-irradiation with stereotactic ablative radiotherapy following high-dose previous thoracic radiotherapy for lung malignancies. *Cancer Treat Rev*. (2015) 41:511–8. doi: 10.1016/j.ctrv.2015.04.002

8. De Bari B, Deantonio L, Bourhis J, Prior JO, Ozsahin M. Should we include SPECT lung perfusion in radiotherapy treatment plans of thoracic targets? Evidences from the literature. *Crit Rev Oncol Hematol*. (2016) 102:111–7. doi: 10.1016/j.critrevonc.2016.04.009

9. Galvin I, Drummond G, Nirmalan M. Distribution of blood flow and ventilation in the lung: gravity is not the only factor. *Br J Anaesth*. (2007) 98:420–8. doi: 10.1093/bja/aem036

10. Grigorios-Aris C, Mariam AM, Kostas H, Anthony AH, Jonas J, Marika B, et al. Automatic segmentation of lungs in SPECT images using active shape model trained by meshes delineated in CT images. In: *2016 IEEE 38th Annual International Conference of the Engineering in Medicine and Biology Society (EMBC)*. Orlando, FL: IEEE (2016). doi: 10.1109/EMBC.2016.7590940

11. He Y, Chen P, Chen Y. Perfusion-ventilation lung SPECT image analysis system based on minimum cross-entropy threshold and watershed segmentation. In: *2008 ISECS International Colloquium on Computing, Communication, Control, and Management*. Guangzhou: IEEE (2008). doi: 10.1109/CCCM.2008.235

12. Wang A, Yan H. SPECT lung delineation via true 3D active contours. *IAENG Int J Comput Sci*. (2010) 37:3.

13. Wang A, Yan H. Delineating low-count defective-contour SPECT lung scans for PE diagnosis using adaptive dual exponential thresholding and active contours. *Int J Imaging Syst Technol*. (2010) 20:149–54. doi: 10.1002/ima.20222

14. Buvat I, Lazaro D. Monte Carlo simulations in emission tomography and GATE: an overview. *Nucl Instr Methods Phys Res A Acceler Spectrom Detect Assoc Equip*. (2006) 569:323–9. doi: 10.1016/j.nima.2006.08.039

15. Norberg P, Olsson A, Carlsson GA, Sandborg M, Gustafsson A. Optimisation of quantitative lung SPECT applied to mild COPD: a software phantom simulation study. *EJNMMI Res*. (2015) 5:16. doi: 10.1186/s13550-015-0086-2

16. Norberg P, Persson HL, Carlsson GA, Bake B, Kentson M, Sandborg M, et al., Quantitative lung SPECT applied on simulated early COPD humans with advanced COPD. *EJNMMI Res*. (2013) 3:28. doi: 10.1186/2191-219X-3-28

17. Holstensson M, Smedby Ö, Poludniowski G, Sanchez-Crespo A, Savitcheva I, Öberg M, et al. Comparison of acquisition protocols for ventilation/perfusion SPECT—a Monte Carlo study. *Phys Med Biol*. (2019) 64:235018. doi: 10.1088/1361-6560/ab36ee

18. Chicco P, Magnussen J, Mackey D, Murray I, Van der Wall H. Recognition of subsegmental scintigraphic defects in virtual lung scintigraphy. *Nucl Med Commun*. (1997) 18:728–33. doi: 10.1097/00006231-199708000-00006

19. Magnussen JS, Chicco P, Palmer AW, Bush V, Mackey DW, Storey G, et al. Single-photon emission tomography of a computerised model of pulmonary embolism. *Eur J Nucl Med*. (1999) 26:1430–8. doi: 10.1007/s002590050475

20. Magnussen JS, Chicco P, Palmer AW, Mackey DW. Optimization of the scintigraphic segmental anatomy of the lungs. *J Nucl Med*. (1997) 38:1987.

21. Magnussen JS, Chicco P, Palmer AW, Mackey DW, Magee M, Murray IPC, et al. Variability of perceived defect size in virtual lung scintigraphy. *J Nucl Med*. (1998) 39:361–5.

22. Magnussen JS, Chicco P, Palmer AW, Van der Wall H, Vu DH. Creation of a three-dimensional model of human segmental lung anatomy. *Am J Roentgenol*. (2000) 174:1333–6. doi: 10.2214/ajr.174.5.1741333

23. Parker JA, Coleman RE, Grady E, Royal HD, Siegel BA, Stabin MG, et al. SNM practice guideline for lung scintigraphy 4.0. *J Nucl Med Technol*. (2012) 40:57–65. doi: 10.2967/jnmt.111.101386

24. Ljungberg M, Strand S, King M. *The SIMIND Monte Carlo Program. Monte Carlo Calculation in Nuclear Medicine: Applications in Diagnostic Imaging*. Bristol: IOP Publishing (1998).

25. Toossi MB, Islamistan JP, Momennezhad M. Ljungberg M, Naseri S. SIMIND Monte Carlo simulation of a single photon emission CT. *J Med Phys Assoc Med Phys India*. (2010) 35:42. doi: 10.4103/0971-6203.55967

26. Ljungberg M. *The SIMIND Monte Carlo Program Manual*. Lund: Lund University (2017). Available online at: http://www2.msf.lu.se/simind/download/Simind_manual.pdf

27. Zubal IG, Harrell CR, Smith EO, Rattner Z, Gindi G, Hoffer PB. Computerized three-dimensional segmented human anatomy. *Med Phys*. (1994) 21:299–302. doi: 10.1118/1.597290

28. Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. *Nat Methods*. (2012) 9:671–5. doi: 10.1038/nmeth.2089

29. Bourhis D, Robin P, Essayan M, Abgral R, Querellou S, Tromeur C, et al. V/Q SPECT for the assessment of regional lung function: generation of normal mean and standard deviation 3-D maps. *Front Med*. (2020) 7:143. doi: 10.3389/fmed.2020.00143

30. Jan S, Santin G, Strul D, Staelens S, Assie K, Autret D, et al. GATE: a simulation toolkit for PET SPECT. *Phys Med Biol*. (2004) 49:4543. doi: 10.1088/0031-9155/49/19/007

31. Brown F. MCNP—A General Monte Carlo N-Particle Transport Code. Los Alamos, NM: Los Alamos National Laboratory (2003).

32. Taheri A, Heidary S, Shahrabi H. Monte Carlo simulation of a SPECT system: GATE. MCNPX or SIMIND? (a comparative study). *J Instrum*. (2017) 12:P12022. doi: 10.1088/1748-0221/12/12/P12022

33. Ejeh JE, van Staden JA, du Raan H. Validation of SIMIND Monte Carlo simulation software for modelling a siemens symbia T SPECT scintillation
camera. In: World Congress on Medical Physics and Biomedical Engineering 2018. Prague: Springer (2019). doi: 10.1007/978-981-10-9035-6_106

34. Ramonaheng K. SIMIND Monte Carlo validation for Lutetium-177 imaging. *J Nucl Med.* (2019) 60:1349.

35. Lau EM, Bailey DL, Bailey EA, Torzillo PJ, Roach PJ, Schembri GP, et al. Pulmonary hypertension leads to a loss of gravity dependent redistribution of regional lung perfusion: a SPECT/CT study. *Heart*. (2014) 100:47–53. doi: 10.1136/heartjnl-2013-304254

36. de Nijs R, Lagerburg V, Klausen TL, Holm S. Improving quantitative dosimetry in 177Lu-DOTATATE SPECT by energy window-based scatter corrections. *Nucl Med Commun.* (2014) 35:522. doi: 10.1097/MNM.0000000000000079

37. Freeman LM, Krynyckyi B, Zuckier LS. Enhanced lung scan diagnosis of pulmonary embolism with the use of ancillary scintigraphic findings and clinical correlation. *Semin Nucl Med.* (2001) 31:143–57. doi: 10.1053/snuc.2001.21273

38. Walrand S, Hesse M, Chiesa C, Lhommel R, Jamar F. The low hepatic toxicity per Gray of 90Y glass microspheres is linked to their transport in the arterial tree favoring a nonuniform trapping as observed in posttherapy PET imaging. *J Nucl Med.* (2014) 55:135–40. doi: 10.2967/jnumed.113.126839

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2020 Bourhis, Wagner, Essayan, Robin, Le Pennec, Salaun and Le Roux. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.