1
|
Gustafsson J, Ljungberg M, Alm Carlsson G, Larsson E, Warfvinge CF, Asp P, Sjögreen Gleisner K. Averaging of absorbed doses: How matter matters. Med Phys 2023; 50:6600-6613. [PMID: 37272586 DOI: 10.1002/mp.16528] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2022] [Revised: 04/05/2023] [Accepted: 05/10/2023] [Indexed: 06/06/2023] Open
Abstract
BACKGROUND Dosimetry in radionuclide therapy often requires the calculation of average absorbed doses within and between spatial regions, for example, for voxel-based dosimetry methods, for paired organs, or across multiple tumors. Formation of such averages can be made in different ways, starting from different definitions. PURPOSE The aim of this study is to formally specify different averaging strategies for absorbed doses, and to compare their results when applied to absorbed dose distributions that are non-uniform within and between regions. METHODS For averaging within regions, two definitions of the average absorbed dose are considered: the simple average over the region (the region average) and the average when weighting by the mass density (density-weighted region average). The latter is shown to follow from the definition of mean absorbed dose according to the ICRU, and to be consistent with the MIRD formalism. For averaging between different spatial regions, three definitions follow: the volume-weighted, the mass-weighted, and the unweighted average. With respect to characterizing non-uniformity, the different average definitions lead to the use of dose-volume histograms (DVHs) (region average), dose-mass histograms (DMHs) (density-weighted region average), and unweighted histograms (unweighted average). Average absorbed doses are calculated for three worked examples, starting from the different definitions. The first, schematic, example concerns the calculation of the average absorbed dose between two regions with different volumes or mass densities. The second, stylized, example concerns voxel-based dosimetry, for which the average absorbed-dose rate within a region is calculated. The geometries studied include three 177 Lu-filled voxelized spheres, where the sphere masses are held constant while the material compositions, densities, and volumes are varied. For comparison, the mean absorbed-dose rates obtained using unit-density sphere S-values are also included. The third example concerns SPECT/CT-based tumor dosimetry for five patients undergoing therapy with 177 Lu-PSMA and six patients undergoing therapy with 177 Lu-DOTA-TATE, for which the average absorbed-dose rates across multiple tumors are calculated. For the second and third examples, analyses also include representations by histograms. RESULTS Example 1 shows that the average absorbed doses, calculated using different definitions, can differ considerably if the masses and absorbed doses for two regions are markedly different. From example 2 it is seen that the density-weighted region average is stable under different activity and density distributions and is also in line with results using S-values. In contrast, the region average varies as function of the activity distribution. In example 3, the absorbed dose rates for individual tumors differ by (1.1 ± 4.3)% and (-0.1 ± 0.4)% with maximum deviations of +34.4% and -1.4% for 177 Lu-PSMA and 177 Lu-DOTA-TATE, respectively, when calculated as region averages or density-weighted region averages, with largest deviations obtained when the density is non-uniform. The average absorbed doses calculated across all tumors are similar when comparing mass-weighted and volume-weighted averages but these differ substantially from unweighted averages. CONCLUSION Different strategies for averaging of absorbed doses within and between regions can lead to substantially different absorbed-dose estimates. At reporting of radionuclide therapy dosimetry, it is important to specify the averaging strategy applied.
Collapse
Affiliation(s)
| | | | - Gudrun Alm Carlsson
- Department of Radiation Physics, Faculty of Health Sciences, Linköping University, Linköping, Sweden
| | - Erik Larsson
- Radiation Physics, Skåne University Hospital, Lund, Sweden
| | - Carl Fredrik Warfvinge
- Division of Oncology and Pathology, Department of Clinical Sciences Lund, Lund University, Lund, Sweden
| | - Pernilla Asp
- Division of Oncology and Pathology, Department of Clinical Sciences Lund, Lund University, Lund, Sweden
| | | |
Collapse
|
2
|
Lyu Y, Chen G, Lu Z, Chen Y, Mok GSP. The effects of mismatch between SPECT and CT images on quantitative activity estimation - A simulation study. Z Med Phys 2023; 33:54-69. [PMID: 35644776 PMCID: PMC10082378 DOI: 10.1016/j.zemedi.2022.03.004] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2021] [Revised: 03/19/2022] [Accepted: 03/25/2022] [Indexed: 01/11/2023]
Abstract
BACKGROUND Quantitative activity estimation is essential in nuclear medicine imaging. Mismatch between SPECT and CT images at the same imaging time point due to patient movement degrades accuracy in both diagnostic studies and target radionuclide therapy dosimetry. This work aims to study the mismatch effects between CT and SPECT data on attenuation correction (AC), volume-of-interest (VOI) delineation, and registration for activity estimation. METHODS Nine 4D XCAT phantoms were generated at 1, 24, and 144 h post In-111 Zevalin injection, varying in activity distributions, body sizes, and organ sizes. Realistic noisy SPECT projections were generated by an analytical projector and reconstructed with a quantitative OS-EM method. CT images were shifted, corresponding to SPECT images at each imaging time point, from -5 to 5 voxels and also according to a clinical reference. The effect of mismatched AC maps was evaluated using mismatched CT images for AC in SPECT reconstruction while VOIs were mapped out from matched CTs. The effect of mismatched VOI drawings was evaluated using mismatched CTs to map out target organs while using matched CTs for AC. The effect of mismatched CT images for registration was evaluated by registering sequential mismatched CTs to align corresponding SPECT images, with no AC and VOI mismatch. Bi-exponential curve fitting was performed to obtain time-integrated activity (TIA). Organ activity errors (%OAE) and TIA errors (%TIAE) were calculated. RESULTS According to the clinical reference, %OAE was larger for organs near ribs for AC effect. For VOI effect, %OAE was larger for small and low uptake organs. For registration effect, %TIAE were larger when mismatch existed in more numbers of SPECT/CT images, while no substantial difference was observed when using mismatched CT at different imaging time points as registration reference. %TIAE was highest for VOI, followed by registration and AC, e.g., 20.62%±8.61%, 9.33%±4.66% and 1.13%±0.90% respectively for kidneys. CONCLUSIONS The mismatch between CT and SPECT images poses a significant impact on the accuracy of quantitative activity estimation, attributed particularly from VOI delineation errors. It is recommended to perform registration between emission and transmission images at the same time point to ensure diagnostic and dosimetric accuracy.
Collapse
Affiliation(s)
- Yingqing Lyu
- Biomedical Imaging Laboratory (BIG), Department of Electrical and Computer Engineering, Faculty of Science and Technology, University of Macau, Taipa, Macau SAR, China
| | - Gefei Chen
- Biomedical Imaging Laboratory (BIG), Department of Electrical and Computer Engineering, Faculty of Science and Technology, University of Macau, Taipa, Macau SAR, China
| | - Zhonglin Lu
- Biomedical Imaging Laboratory (BIG), Department of Electrical and Computer Engineering, Faculty of Science and Technology, University of Macau, Taipa, Macau SAR, China
| | - Yue Chen
- Department of Nuclear Medicine, The Affiliated Hospital of Southwest Medical University, Nuclear Medicine and Molecular Imaging Key Laboratory of Sichuan Province, No. 25, Taiping St., Luzhou, Sichuan, China.
| | - Greta S P Mok
- Biomedical Imaging Laboratory (BIG), Department of Electrical and Computer Engineering, Faculty of Science and Technology, University of Macau, Taipa, Macau SAR, China; Center for Cognitive and Brain Sciences, Institute of Collaborative Innovation, University of Macau, Taipa, Macau SAR, China; Ministry of Education Frontiers Science Center for Precision Oncology, Faculty of Health Science, University of Macau, Taipa, Macau SAR, China.
| |
Collapse
|
3
|
EANM dosimetry committee recommendations for dosimetry of 177Lu-labelled somatostatin-receptor- and PSMA-targeting ligands. Eur J Nucl Med Mol Imaging 2022; 49:1778-1809. [PMID: 35284969 PMCID: PMC9015994 DOI: 10.1007/s00259-022-05727-7] [Citation(s) in RCA: 79] [Impact Index Per Article: 39.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2021] [Accepted: 02/13/2022] [Indexed: 12/25/2022]
Abstract
The purpose of the EANM Dosimetry Committee is to provide recommendations and guidance to scientists and clinicians on patient-specific dosimetry. Radiopharmaceuticals labelled with lutetium-177 (177Lu) are increasingly used for therapeutic applications, in particular for the treatment of metastatic neuroendocrine tumours using ligands for somatostatin receptors and prostate adenocarcinoma with small-molecule PSMA-targeting ligands. This paper provides an overview of reported dosimetry data for these therapies and summarises current knowledge about radiation-induced side effects on normal tissues and dose-effect relationships for tumours. Dosimetry methods and data are summarised for kidneys, bone marrow, salivary glands, lacrimal glands, pituitary glands, tumours, and the skin in case of radiopharmaceutical extravasation. Where applicable, taking into account the present status of the field and recent evidence in the literature, guidance is provided. The purpose of these recommendations is to encourage the practice of patient-specific dosimetry in therapy with 177Lu-labelled compounds. The proposed methods should be within the scope of centres offering therapy with 177Lu-labelled ligands for somatostatin receptors or small-molecule PSMA.
Collapse
|
4
|
Pettinato C, Richetta E, Cremonesi M. Dosimetry with single photon emission tomography (SPECT). Nucl Med Mol Imaging 2022. [DOI: 10.1016/b978-0-12-822960-6.00173-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022] Open
|
5
|
SCreg: a registration-based platform to compare unicondylar knee arthroplasty SPECT/CT scans. BMC Musculoskelet Disord 2020; 21:162. [PMID: 32164663 PMCID: PMC7066757 DOI: 10.1186/s12891-020-3185-9] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/05/2019] [Accepted: 03/03/2020] [Indexed: 11/11/2022] Open
Abstract
Background A combination of conventional computed tomography and single photon emitted computed tomography (SPECT/CT) provides simultaneous data on the intensity and location of osteoblastic activity. Currently, since SPECT/CT scans are not spatially aligned, scans following knee arthroplasty are compared by extracting average and maximal values of osteoblastic activity intensity from large subregions of the structure of interest, which leads to a loss of resolution, and hence, information. Therefore, this paper describes the SPECT/CT registration platform (SCreg) based on the principle of image registration to spatially align SPECT/CT scans following unicondylar knee arthroplasty (UKA) and allow full resolution intra-subject and inter-subject comparisons. Methods SPECT-CT scans of 20 patients were acquired before and 1 year after UKA. Firstly, scans were pre-processed to account for differences in voxel sizes and divided in volumes of interest. This was followed by optimization of registration parameters according to their volumetric agreement, and alignment using a combination of rigid, affine and non-rigid registration. Finally, radiotracer uptakes were normalized, and differences between pre-operative and post-operative activity were computed for each voxel. Wilcoxon signed rank sum test was performed to compare Dice similarity coefficients pre- and post-registration. Results Qualitative and quantitative validation of the platform assessing the correct alignment of SPECT/CT scans resulted in Dice similarity coefficient values over 80% and distances between predefined anatomical landmarks below the fixed threshold of (2;2;0) voxels. Locations of increased and decreased osteoblastic activity obtained during comparisons of osteoblastic activity before and after UKA were mainly consistent with literature. Conclusions Thus, a full resolution comparison performed on the platform could assist surgeons and engineers in optimizing surgical parameters in view of bone remodeling, thereby improving UKA survivorship.
Collapse
|
6
|
Gustafsson J, Brolin G, Ljungberg M. Monte Carlo-based SPECT reconstruction within the SIMIND framework. Phys Med Biol 2018; 63:245012. [PMID: 30523946 DOI: 10.1088/1361-6560/aaf0f1] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
This paper presents the development and validation of a Monte Carlo-based singe photon emission computed tomography reconstruction program for parallel-hole collimation contained within the SIMIND Monte Carlo framework. The Monte Carlo code is used as an accurate forward-projector and is combined with a simplified back-projector to perform iterative tomographic reconstruction using the Maximum Likelihood Expectation Maximization and Ordered Subsets Expectation Maximization algorithms, together forming a program called SIMREC. The Monte Carlo simulation transforms the estimated source distribution directly from activity to counts in its projections. Hence, the reconstructed image is expressed in activity without reference to an external calibration. The program is tested using phantom measurements of spheres filled with 99mTc, 177Lu and 131I placed in air and centrally and peripherally in a water-filled elliptical phantom. The feasibility of applying the reconstruction to patients is also demonstrated for a range of radiopharmaceuticals. The deviation in total activity in the spheres ranged between -4.1% and 6.2% compared with the activity determined when preparing the phantom. The SIMREC program was found to be accurate with respect to activity estimation and to reconstruct visually acceptable images within a few hours when applied to patient examples.
Collapse
Affiliation(s)
- Johan Gustafsson
- Department of Medical Radiation Physics, Lund University, Lund, Sweden. Author to whom any correspondence should be addressed
| | | | | |
Collapse
|
7
|
Ljungberg M, Sjogreen Gleisner K. 3-D Image-Based Dosimetry in Radionuclide Therapy. IEEE TRANSACTIONS ON RADIATION AND PLASMA MEDICAL SCIENCES 2018. [DOI: 10.1109/trpms.2018.2860563] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
8
|
Li T, Wu NY, Song N, Mok GSP. Evaluation of sequential SPECT and CT for targeted radionuclide therapy dosimetry. Ann Nucl Med 2017; 32:34-43. [PMID: 29143283 DOI: 10.1007/s12149-017-1218-8] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2017] [Accepted: 11/08/2017] [Indexed: 12/22/2022]
Abstract
PURPOSE In targeted radionuclide therapy (TRT), a prior knowledge of the absorbed dose biodistribution is essential for pre-therapy treatment planning. Previously, we showed that non-rigid organ-by-organ registration in sequential quantitative SPECT images improved dose estimation. This study aims to investigate if sequential CT can further improve TRT dosimetric accuracy. METHODS We simulated SPECT/CT acquisitions at 1, 12, 24, 72 and 144 h In-111 Zevalin post-injection using an analytical MEGP projector, modeling attenuation, scatter and collimator-detector response. We later recruited a patient injected with 222 MBq In-111 DTPAOC imaged at 3 SPECT/CT sessions for clinical evaluations. Four registration schemes were evaluated: whole-body-based registration performed on sequential (1) SPECT (WB-SPECT) or (2) CT (WB-CT) images; organ-based registration applied on organs individually segmented from sequential (3) SPECT (O-SPECT) or (4) CT (O-CT) images. Voxel-by-voxel integration was performed followed by Y-90 voxel-S-kernel convolution. Organ-absorbed doses, iso-dose curves, dose-volume histograms (DVHs) were generated for targeted organs for analysis. RESULTS In simulation study, organ-absorbed dose errors were (- 8.66 ± 2.83)%, (- 2.51 ± 3.69)%, (- 9.23 ± 3.28)%, (- 7.17 ± 2.53)% for liver, (- 14.81 ± 4.91)%, (- 3.60 ± 4.37)%, (- 18.13 ± 4.44)%, (- 11.34 ± 4.22)% for spleen, for O-SPECT, O-CT, WB-SPECT and WB-CT registrations, respectively. For all organs, O-CT showed superior results. Results of iso-dose contour, DVHs were in accordance with the organ-absorbed doses. In clinical studies, the results were also consistent which showed O-CT method deviated the most from the result with no registration. CONCLUSIONS We conclude that if both sequential SPECT/CT scans are available, CT organ-based registration method can more effectively improve the 3D dose estimation. Sequential low-dose CT scans might be considered to be included in the standard TRT protocol.
Collapse
Affiliation(s)
- Tiantian Li
- Biomedical Imaging Laboratory, Department of Electrical and Computer Engineering, Faculty of Science and Technology, University of Macau, Macau SAR, China
| | - Nien-Yun Wu
- Department of Biomedical Imaging and Radiological Sciences, National Yang Ming University, Taipei, Taiwan, Republic of China.,Department of Nuclear Medicine, Taipei Veterans General Hospital, Taipei, Taiwan, Republic of China
| | - Na Song
- Department of Nuclear Medicine, Albert Einstein College of Medicine, Yeshiva University, Bronx, New York, 10461, USA
| | - Greta S P Mok
- Biomedical Imaging Laboratory, Department of Electrical and Computer Engineering, Faculty of Science and Technology, University of Macau, Macau SAR, China. .,Faculty of Health Sciences, University of Macau, Macau SAR, China.
| |
Collapse
|
9
|
Sarrut D, Halty A, Badel JN, Ferrer L, Bardiès M. Voxel-based multimodel fitting method for modeling time activity curves in SPECT images. Med Phys 2017; 44:6280-6288. [PMID: 28921613 DOI: 10.1002/mp.12586] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2017] [Revised: 08/22/2017] [Accepted: 08/22/2017] [Indexed: 11/10/2022] Open
Abstract
PURPOSE Estimating the biodistribution and the pharmacokinetics from time-sequence SPECT images on a per-voxel basis is useful for studying activity nonuniformity or computing absorbed dose distributions by convolution of voxel kernels or Monte-Carlo radiation transport. Current approaches are either region-based, thus assuming uniform activity within the region, or voxel-based but using the same fitting model for all voxels. METHODS We propose a voxel-based multimodel fitting method (VoMM) that estimates a fitting function for each voxel by automatically selecting the most appropriate model among a predetermined set with Akaike criteria. This approach can be used to compute the time integrated activity (TIA) for all voxels in the image. To control fitting optimization that may fail due to excessive image noise, an approximated version based on trapezoid integration, named restricted method, is also studied. From this comparison, the number of failed fittings within images was estimated and analyzed. Numerical experiments were used to quantify uncertainties and feasibility was demonstrated with real patient data. RESULTS Regarding numerical experiments, root mean square errors of TIA obtained with VoMM were similar to those obtained with bi-exponential fitting functions, and were lower (< 5% vs. > 10%) than with single model approaches that consider the same fitting function for all voxels. Failure rates were lower with VoMM and restricted approaches than with single-model methods. On real clinical data, VoMM was able to fit 90% of the voxels and led to less failed fits than single-model approaches. On regions of interest (ROI) analysis, the difference between ROI-based and voxel-based TIA estimations was low, less than 4%. However, the computation of the mean residence time exhibited larger differences, up to 25%. CONCLUSIONS The proposed voxel-based multimodel fitting method, VoMM, is feasible on patient data. VoMM leads organ-based TIA estimations similar to conventional ROI-based method. However, for pharmacokinetics analysis, studies of spatial heterogeneity or voxel-based absorbed dose assessment, VoMM could be used preferentially as it prevents model overfitting.
Collapse
Affiliation(s)
- David Sarrut
- Université de Lyon, CREATIS, CNRS UMR5220, Inserm U1044, INSA-Lyon, Université Lyon 1, Centre Léon Bérard, Lyon, France
| | - Adrien Halty
- Université de Lyon, CREATIS, CNRS UMR5220, Inserm U1044, INSA-Lyon, Université Lyon 1, Centre Léon Bérard, Lyon, France
| | - Jean-Noel Badel
- Université de Lyon, CREATIS, CNRS UMR5220, Inserm U1044, INSA-Lyon, Université Lyon 1, Centre Léon Bérard, Lyon, France
| | - Ludovic Ferrer
- Medical Physics Department, ICO René Gauducheau, St Herblain, France.,CRCINA UMR 892, Nantes, France
| | - Manuel Bardiès
- Inserm UMR1037, CRCT, Université Toulouse III-Paul Sabatier, F-31000, Toulouse, France
| |
Collapse
|
10
|
Li T, Ao ECI, Lambert B, Brans B, Vandenberghe S, Mok GSP. Quantitative Imaging for Targeted Radionuclide Therapy Dosimetry - Technical Review. Theranostics 2017; 7:4551-4565. [PMID: 29158844 PMCID: PMC5695148 DOI: 10.7150/thno.19782] [Citation(s) in RCA: 50] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2017] [Accepted: 07/25/2017] [Indexed: 01/06/2023] Open
Abstract
Targeted radionuclide therapy (TRT) is a promising technique for cancer therapy. However, in order to deliver the required dose to the tumor, minimize potential toxicity in normal organs, as well as monitor therapeutic effects, it is important to assess the individualized internal dosimetry based on patient-specific data. Advanced imaging techniques, especially radionuclide imaging, can be used to determine the spatial distribution of administered tracers for calculating the organ-absorbed dose. While planar scintigraphy is still the mainstream imaging method, SPECT, PET and bremsstrahlung imaging have promising properties to improve accuracy in quantification. This article reviews the basic principles of TRT and discusses the latest development in radionuclide imaging techniques for different theranostic agents, with emphasis on their potential to improve personalized TRT dosimetry.
Collapse
Affiliation(s)
- Tiantian Li
- Biomedical Imaging Laboratory, Department of Electrical and Computer Engineering, Faculty of Science and Technology, University of Macau, Macau SAR, China
| | - Edwin C. I. Ao
- Biomedical Imaging Laboratory, Department of Electrical and Computer Engineering, Faculty of Science and Technology, University of Macau, Macau SAR, China
| | - Bieke Lambert
- Dept of Radiology and Nuclear medicine, Ghent University, De Pintelaan 185 9000 Gent, Belgium
- AZ Maria Middelares, Buiten-Ring-Sint-Denijs 30, 9000 Gent, Belgium
| | - Boudewijn Brans
- Dept of Nuclear Medicine, UZ Ghent-Ghent University, St-Pietersnieuwstraat 41, 9000 Gent, Belgium
| | - Stefaan Vandenberghe
- MEDISIP-ELIS-IBITECH-IMEC, Ghent University, St-Pietersnieuwstraat 41, 9000 Gent, Belgium
| | - Greta S. P. Mok
- Biomedical Imaging Laboratory, Department of Electrical and Computer Engineering, Faculty of Science and Technology, University of Macau, Macau SAR, China
- Faculty of Health Sciences, University of Macau, Macau SAR, China
| |
Collapse
|
11
|
Gustafsson J, Sundlöv A, Sjögreen Gleisner K. SPECT image segmentation for estimation of tumour volume and activity concentration in 177Lu-DOTATATE radionuclide therapy. EJNMMI Res 2017; 7:18. [PMID: 28233160 PMCID: PMC5323339 DOI: 10.1186/s13550-017-0262-7] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2016] [Accepted: 01/31/2017] [Indexed: 11/23/2022] Open
Abstract
Background Dosimetry in radionuclide therapy has the potential to allow for a treatment tailored to the individual patient. One therapeutic radiopharmaceutical where patient-specific dosimetry is feasible is 177Lu-DOTATATE, used for the treatment of neuroendocrine tumours. The emission of gamma photons by 177Lu allows for imaging with SPECT (single photon emission computed tomography). One important step for dosimetry using this imaging technique is the SPECT image segmentation, which needs to be robust and accurate for the estimated quantities to be reliable. This work investigates different methods for automatic tumour delineation in 177Lu-DOTATATE SPECT images. Three segmentation methods are considered: a fixed 42% threshold (FT), the Otsu method (OM) and a method based on Fourier surfaces (FS). Effects of including resolution compensation in the iterative SPECT image reconstruction are also studied. Evaluation is performed based on Monte Carlo-simulated SPECT images from 24 h and 336 h post injection (p.i.), for determination of the volume, activity concentration and dice similarity coefficient. In addition, patient data are used to investigate the correspondence of tumour volumes when delineated in SPECT or morphological CT or MR images. Patient data are also used to examine the sensitivity to the operator-dependent initialization. Results For simulated images from 24 h p.i. reconstructed without resolution compensation, a volume and activity-concentration root-mean-square error below 15% is typically obtained for tumours above approximately 10 cm3 when using OM or FS, while FT performs considerably worse. When including resolution compensation, the tumour volume becomes underestimated and the activity concentration overestimated. The FS method appears to be robust to noise, as seen for the 336 h images. The differences between the tumour volumes estimated from the SPECT images and the volumes estimated from morphological images are generally larger than the discrepancies seen for the simulated data sets. Conclusions Segmentation results are encouraging for future dosimetry of tumours with volumes above approximately 10 cm3. Using resolution compensation in the reconstruction may have a negative effect on volume estimation.
Collapse
Affiliation(s)
- Johan Gustafsson
- Department of Medical Radiation Physics, Clinical Sciences Lund, Lund University, Lund, Sweden.
| | - Anna Sundlöv
- Department of Oncology and Pathology, Clinical Sciences Lund, Lund University, Lund, Sweden.,Department of Oncology, Skåne University Hospital, Lund, Sweden
| | | |
Collapse
|
12
|
Sarrut D, Badel JN, Halty A, Garin G, Perol D, Cassier P, Blay JY, Kryza D, Giraudet AL. 3D absorbed dose distribution estimated by Monte Carlo simulation in radionuclide therapy with a monoclonal antibody targeting synovial sarcoma. EJNMMI Phys 2017; 4:6. [PMID: 28101733 PMCID: PMC5241566 DOI: 10.1186/s40658-016-0172-1] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2016] [Accepted: 12/28/2016] [Indexed: 11/10/2022] Open
Abstract
Backround Radiolabeled OTSA101, a monoclonal antibody targeting synovial sarcoma (SS) developed by OncoTherapy Science, was used to treat relapsing SS metastases following a theranostic procedure: in case of significant 111In-OTSA101 tumor uptake and favorable biodistribution, patient was randomly treated with 370/1110 MBq 90Y-OTSA101. Monte Carlo-based 3D dosimetry integrating time-activity curves in VOI was performed on 111In-OTSA101 repeated SPECT/CT. Estimated absorbed doses (AD) in normal tissues were compared to biological side effects and to the admitted maximal tolerated absorbed dose (MTD) in normal organs. Results in the tumors were also compared to disease evolution. Results Biodistribution and tracer quantification were analyzed on repeated SPECT/CT acquisitions performed after injection of 111In-OTSA101 in 19/20 included patients. SPECT images were warped to a common coordinates system with deformable registration. Volumes of interest (VOI) for various lesions and normal tissues were drawn on the first CT acquisition and reported to all the SPECT images. Tracer quantification and residence time of 111In-OTSA101 in VOI were used to evaluate the estimated absorbed doses per MBq of 90Y-OTSA101 by means of Monte Carlo simulations (GATE). A visual scale analysis was applied to assess tumor uptake (grades 0 to 4) and results were compared to the automated quantification. Results were then compared to biological side effects reported in the selected patients treated with 90Y-OTSA101 but also to disease response to treatment. After screening, 8/20 patients were treated with 370 or 1110 MBq 90Y-OTSA101. All demonstrated medullary toxicity, only one presented with transient grade 3 liver toxicity due to disease progression, and two patients presented with transient grade 1 renal toxicity. Median absorbed doses were the highest in the liver (median, 0.64 cGy/MBq; [0.27 −1.07]) being far lower than the 20 Gy liver MTD, and the lowest in bone marrow (median, 0.09 cGy/MBq; [0.02 −0.18]) being closer to the 2 Gy bone marrow MTD. Most of the patients demonstrated progressive disease on RECIST criteria during patient follow-up. 111In-OTSA101 tumors tracer uptake visually appeared highly heterogeneous in inter- and intra-patient analyses, independently of tumor sizes, with variable kinetics. The majority of visual grades corresponded to the automated computed ones. Estimated absorbed doses in the 95 supra-centimetric selected lesions ranged from 0.01 to 0.71 cGy per injected MBq (median, 0.22 cGy/MBq). The maximal tumor AD obtained was 11.5 Gy. Conclusions 3D dosimetry results can explain the observed toxicity and tumors response. Despite an intense visual 111In-OTSA101 liver uptake, liver toxicity was not the dose limiting factor conversely to bone marrow toxicity. Even though tumors 111In-OTSA101 avidity was visually obvious for treated patients, the low estimated tumors AD obtained by 3D dosimetry explain the lack of tumor response.
Collapse
Affiliation(s)
- David Sarrut
- Univ Lyon, INSA-Lyon, Université Lyon 1, CNRS, Inserm, CREATIS UMR 5220, U1206, Lyon, F-69008, France. .,Univ Lyon, Centre Léon Bérard, Lyon, 69008, France.
| | | | - Adrien Halty
- Univ Lyon, INSA-Lyon, Université Lyon 1, CNRS, Inserm, CREATIS UMR 5220, U1206, Lyon, F-69008, France.,Univ Lyon, Centre Léon Bérard, Lyon, 69008, France
| | | | - David Perol
- Univ Lyon, Centre Léon Bérard, Lyon, 69008, France
| | | | | | - David Kryza
- Univ Lyon, Université Lyon 1, CNRS, LAGEP UMR 5007, Lyon, F-69008, France.,Hospices Civils de Lyon, Imthernat plateform, Lyon, F-69008, France
| | | |
Collapse
|
13
|
Personalized Dosimetry for Radionuclide Therapy Using Molecular Imaging Tools. Biomedicines 2016; 4:biomedicines4040025. [PMID: 28536392 PMCID: PMC5344265 DOI: 10.3390/biomedicines4040025] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2016] [Revised: 11/07/2016] [Accepted: 11/08/2016] [Indexed: 11/16/2022] Open
Abstract
For treatment of systemic malignancies, when external radiation therapy is not applicable, radionuclide therapy can be an alternative. In this form of therapy, radionuclides are administered to the patient, often in a form where the radionuclide is labelled to a molecule that plays the active part in the localization of the tumor. Since the aim is to impart lethal damage to tumor cells while maintaining possible side-effects to normal tissues at tolerable levels, a proper and accurate personalized dosimetry should be a pre-requisite. In radionuclide therapy, there is a need to measure the distribution of the radiopharmaceutical in vivo, as well as its re-distribution over time, in order estimate the total energy released in radioactive decays and subsequent charged-particle interactions, governing the absorbed dose to different organs and tumors. Measurements are usually performed by molecular imaging, more specifically planar and SPECT (Single-Photon Emission Computed Tomography) imaging, combined with CT. This review describes the different parts in the dosimetry chain of radionuclide therapy. Emphasis is given to molecular imaging tools and the requirements for determining absorbed doses from quantitative planar and SPECT images. As example solutions to the different problems that need to be addressed in such a dosimetric chain, we describe our tool, Lundadose, which is a set of methods that we have developed for personalized dosimetry.
Collapse
|
14
|
Ljungberg M, Celler A, Konijnenberg MW, Eckerman KF, Dewaraja YK, Sjögreen-Gleisner K. MIRD Pamphlet No. 26: Joint EANM/MIRD Guidelines for Quantitative 177Lu SPECT Applied for Dosimetry of Radiopharmaceutical Therapy. J Nucl Med 2015; 57:151-62. [DOI: 10.2967/jnumed.115.159012] [Citation(s) in RCA: 176] [Impact Index Per Article: 19.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2015] [Accepted: 10/01/2015] [Indexed: 11/16/2022] Open
|
15
|
Gustafsson J, Brolin G, Cox M, Ljungberg M, Johansson L, Gleisner KS. Uncertainty propagation for SPECT/CT-based renal dosimetry in177Lu peptide receptor radionuclide therapy. Phys Med Biol 2015; 60:8329-46. [DOI: 10.1088/0031-9155/60/21/8329] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
|
16
|
Ljungberg M, Gleisner KS. Hybrid Imaging for Patient-Specific Dosimetry in Radionuclide Therapy. Diagnostics (Basel) 2015; 5:296-317. [PMID: 26854156 PMCID: PMC4665601 DOI: 10.3390/diagnostics5030296] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2015] [Revised: 06/29/2015] [Accepted: 07/01/2015] [Indexed: 11/16/2022] Open
Abstract
Radionuclide therapy aims to treat malignant diseases by systemic administration of radiopharmaceuticals, often using carrier molecules such as peptides and antibodies. The radionuclides used emit electrons or alpha particles as a consequence of radioactive decay, thus leading to local energy deposition. Administration to individual patients can be tailored with regards to the risk of toxicity in normal organs by using absorbed dose planning. The scintillation camera, employed in planar imaging or single-photon emission computed tomography (SPECT), generates images of the spatially and temporally varying activity distribution. Recent commercially available combined SPECT and computed tomography (CT) systems have dramatically increased the possibility of performing accurate dose planning by using the CT information in several steps of the dose-planning calculation chain. This paper discusses the dosimetry chain used for individual absorbed-dose planning and highlights the areas where hybrid imaging makes significant contributions.
Collapse
Affiliation(s)
- Michael Ljungberg
- Department of Medical Radiation Physics, Lund University, 221 85 Lund, Sweden.
| | | |
Collapse
|
17
|
Kost SD, Dewaraja YK, Abramson RG, Stabin MG. VIDA: a voxel-based dosimetry method for targeted radionuclide therapy using Geant4. Cancer Biother Radiopharm 2015; 30:16-26. [PMID: 25594357 PMCID: PMC4322792 DOI: 10.1089/cbr.2014.1713] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
We have developed the Voxel-Based Internal Dosimetry Application (VIDA) to provide patient-specific dosimetry in targeted radionuclide therapy performing Monte Carlo simulations of radiation transport with the Geant4 toolkit. The code generates voxel-level dose rate maps using anatomical and physiological data taken from individual patients. Voxel level dose rate curves are then fit and integrated to yield a spatial map of radiation absorbed dose. In this article, we present validation studies using established dosimetry results, including self-dose factors (DFs) from the OLINDA/EXM program for uniform activity in unit density spheres and organ self- and cross-organ DFs in the Radiation Dose Assessment Resource (RADAR) reference adult phantom. The comparison with reference data demonstrated agreement within 5% for self-DFs to spheres and reference phantom source organs for four common radionuclides used in targeted therapy ((131)I, (90)Y, (111)In, (177)Lu). Agreement within 9% was achieved for cross-organ DFs. We also present dose estimates to normal tissues and tumors from studies of two non-Hodgkin Lymphoma patients treated by (131)I radioimmunotherapy, with comparison to results generated independently with another dosimetry code. A relative difference of 12% or less was found between methods for mean absorbed tumor doses accounting for tumor regression.
Collapse
Affiliation(s)
- Susan D. Kost
- Department of Physics and Astronomy, Vanderbilt University, Nashville, Tennessee
| | - Yuni K. Dewaraja
- Department of Radiology, University of Michigan, Ann Arbor, Michigan
| | - Richard G. Abramson
- Department of Radiology/Radiological Science, Vanderbilt University, Nashville, Tennessee
| | - Michael G. Stabin
- Department of Radiology/Radiological Science, Vanderbilt University, Nashville, Tennessee
| |
Collapse
|
18
|
Ao ECI, Wu NY, Wang SJ, Song N, Mok GSP. Improved dosimetry for targeted radionuclide therapy using nonrigid registration on sequential SPECT images. Med Phys 2015; 42:1060-70. [DOI: 10.1118/1.4906242] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
19
|
Grimes J, Uribe C, Celler A. JADA: a graphical user interface for comprehensive internal dose assessment in nuclear medicine. Med Phys 2014; 40:072501. [PMID: 23822450 DOI: 10.1118/1.4810963] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE The main objective of this work was to design a comprehensive dosimetry package that would keep all aspects of internal dose calculation within the framework of a single software environment and that would be applicable for a variety of dose calculation approaches. METHODS Our MATLAB-based graphical user interface (GUI) can be used for processing data obtained using pure planar, pure SPECT, or hybrid planar/SPECT imaging. Time-activity data for source regions are obtained using a set of tools that allow the user to reconstruct SPECT images, load images, coregister a series of planar images, and to perform two-dimensional and three-dimensional image segmentation. Curve fits are applied to the acquired time-activity data to construct time-activity curves, which are then integrated to obtain time-integrated activity coefficients. Subsequently, dose estimates are made using one of three methods. RESULTS The organ level dose calculation subGUI calculates mean organ doses that are equivalent to dose assessment performed by OLINDA/EXM. Voxelized dose calculation options, which include the voxel S value approach and Monte Carlo simulation using the EGSnrc user code DOSXYZnrc, are available within the process 3D image data subGUI. CONCLUSIONS The developed internal dosimetry software package provides an assortment of tools for every step in the dose calculation process, eliminating the need for manual data transfer between programs. This saves times and minimizes user errors, while offering a versatility that can be used to efficiently perform patient-specific internal dose calculations in a variety of clinical situations.
Collapse
Affiliation(s)
- Joshua Grimes
- Department of Physics and Astronomy, University of British Columbia, Vancouver V5Z 1M9, Canada.
| | | | | |
Collapse
|
20
|
Gleisner KS, Ljungberg M. Patient-specific whole-body attenuation correction maps from a CT system for conjugate-view-based activity quantification: method development and evaluation. Cancer Biother Radiopharm 2013; 27:652-64. [PMID: 23210927 DOI: 10.1089/cbr.2011.1082] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
For activity quantification based on planar scintillation camera measurements, photon attenuation is an important factor that needs to be corrected for in a patient- and organ-specific manner. One possibility for obtaining attenuation correction maps is to use X-ray CT scout images. Since the intensity of scout images is in relative numbers, their image values need to be multiplied by a factor to become quantitative and thus useful for attenuation correction. The calibration factor can for our current imaging system be obtained from a scanner system file, but is generally not available. For this purpose, a method based on the patient weight has been developed. Results based on 79 patient scout images show that the calibration factor thus determined correlates well with values that, in this case, are independently specified by the system. The accuracy of attenuation correction factors (ACFs) derived from the scout-based attenuation correction maps is evaluated by comparison to ACFs derived from three-dimensional CT studies. For photon energies of 208, 245, and 364 keV, scout-based ACFs are on average 1.2% and 0.5% from the CT-derived values, using the system-based and the weight-based values of the scout-image calibration factor, respectively. The imprecision is somewhat higher for the weight-based method, due to variability in the delineation of the patient contour used as a part of this method. In conclusion, X-ray scouts are found useful for attenuation correction with a satisfactory accuracy obtained, both using the new, weight-based method, and using the previous, system-based method, for determining the required calibration factor.
Collapse
|
21
|
Dewaraja YK, Frey EC, Sgouros G, Brill AB, Roberson P, Zanzonico PB, Ljungberg M. MIRD pamphlet No. 23: quantitative SPECT for patient-specific 3-dimensional dosimetry in internal radionuclide therapy. J Nucl Med 2012; 53:1310-25. [PMID: 22743252 DOI: 10.2967/jnumed.111.100123] [Citation(s) in RCA: 241] [Impact Index Per Article: 20.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
In internal radionuclide therapy, a growing interest in voxel-level estimates of tissue-absorbed dose has been driven by the desire to report radiobiologic quantities that account for the biologic consequences of both spatial and temporal nonuniformities in these dose estimates. This report presents an overview of 3-dimensional SPECT methods and requirements for internal dosimetry at both regional and voxel levels. Combined SPECT/CT image-based methods are emphasized, because the CT-derived anatomic information allows one to address multiple technical factors that affect SPECT quantification while facilitating the patient-specific voxel-level dosimetry calculation itself. SPECT imaging and reconstruction techniques for quantification in radionuclide therapy are not necessarily the same as those designed to optimize diagnostic imaging quality. The current overview is intended as an introduction to an upcoming series of MIRD pamphlets with detailed radionuclide-specific recommendations intended to provide best-practice SPECT quantification-based guidance for radionuclide dosimetry.
Collapse
Affiliation(s)
- Yuni K Dewaraja
- Department of Radiology, University of Michigan Medical School, Ann Arbor, MI 48109, USA
| | | | | | | | | | | | | |
Collapse
|
22
|
Minarik D, Sjögreen-Gleisner K, Linden O, Wingårdh K, Tennvall J, Strand SE, Ljungberg M. 90Y Bremsstrahlung Imaging for Absorbed-Dose Assessment in High-Dose Radioimmunotherapy. J Nucl Med 2010; 51:1974-8. [DOI: 10.2967/jnumed.110.079897] [Citation(s) in RCA: 49] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
|