RESUMO
BACKGROUND: Stereotactic MR-guided Adaptive Radiation Therapy (SMART) dose painting for hypoxia has potential to improve treatment outcomes, but clinical implementation on low-field MR-Linac faces substantial challenges due to dramatically lower signal-to-noise ratio (SNR) characteristics. While quantitative MRI and T1 mapping of hypoxia biomarkers show promise, T1-to-noise ratio (T1NR) optimization at low fields is paramount, particularly for the clinical implementation of oxygen-enhanced (OE)-MRI. The 3D Magnetization Prepared (2) Rapid Gradient Echo (MP2RAGE) sequence stands out for its ability to acquire homogeneous T1-weighted contrast images with simultaneous T1 mapping. PURPOSE: To optimize MP2RAGE for low-field T1 mapping; conduct experimental validation in a ground-truth phantom; establish feasibility and reproducibility of low-field MP2RAGE acquisition and T1 mapping in healthy volunteers. METHODS: The MP2RAGE optimization was performed to maximize the contrast-to-noise ratio (CNR) of T1 values in white matter (WM) and gray matter (GM) brain tissues at 0.35T. Low-field MP2RAGE images were acquired on a 0.35T MR-Linac (ViewRay MRIdian) using a multi-channel head coil. Validation of T1 mapping was performed with a ground-truth Eurospin phantom, containing inserts of known T1 values (400-850 ms), with one and two average (1A and 2A) MP2RAGE scans across four acquisition sessions, resulting in eight T1 maps. Mean (± SD) T1 relative error, T1NR, and intersession coefficient of variation (CV) were determined. Whole-brain MP2RAGE scans were acquired in 5 healthy volunteers across two sessions (A and B) and T1 maps were generated. Mean (± SD) T1 values for WM and GM were determined. Whole-brain T1 histogram analysis was performed, and reproducibility was determined with the CV between sessions. Voxel-by-voxel T1 difference maps were generated to evaluate 3D spatial variation. RESULTS: Low-field MP2RAGE optimization resulted in parameters: MP2RAGETR of 3250 ms, inversion times (TI1/TI2) of 500/1200 ms, and flip angles (α1/α2) of 7/5°. Eurospin T1 maps exhibited a mean (± SD) relative error of 3.45% ± 1.30%, T1NR of 20.13 ± 5.31, and CV of 2.22% ± 0.67% across all inserts. Whole-brain MP2RAGE images showed high anatomical quality with clear tissue differentiation, resulting in mean (± SD) T1 values: 435.36 ± 10.01 ms for WM and 623.29 ± 14.64 ms for GM across subjects, showing excellent concordance with literature. Whole-brain T1 histograms showed high intrapatient and intersession reproducibility with characteristic intensity peaks consistent with voxel-level WM and GM T1 values. Reproducibility analysis revealed a CV of 0.46% ± 0.31% and 0.35% ± 0.18% for WM and GM, respectively. Voxel-by-voxel T1 difference maps show a normal 3D spatial distribution of noise in WM and GM. CONCLUSIONS: Low-field MP2RAGE proved effective in generating accurate, reliable, and reproducible T1 maps with high T1NR in phantom studies and in vivo feasibility established in healthy volunteers. While current work is focused on refining the MP2RAGE protocol to enable clinically efficient OE-MRI, this study establishes a foundation for TOLD T1 mapping for hypoxia biomarkers. This advancement holds the potential to facilitate a paradigm shift toward MR-guided biological adaptation and dose painting by leveraging 3D hypoxic spatial distributions and improving outcomes in conventionally challenging-to-treat cancers.
RESUMO
Objective.Proton therapy allows for highly conformal dose deposition, but is sensitive to range uncertainties. Several approaches currently under development measure composition-dependent secondary radiation to monitor the delivered proton rangein-vivo. To fully utilize these methods, an estimate of the elemental composition of the patient's tissue is often needed.Approach.A published dual-energy computed tomography (DECT)-based composition-extraction algorithm was validated against reference compositions obtained with two independent methods. For this purpose, a set of phantoms containing either fresh porcine tissue or tissue-mimicking samples with known, realistic compositions were imaged with a CT scanner at two different energies. Then, the prompt gamma-ray (PG) signal during proton irradiation was measured with a PG detector prototype. The PG workflow used pre-calculated Monte Carlo simulations to obtain an optimized estimate of the sample's carbon and oxygen contents. The compositions were also assessed with chemical combustion analysis (CCA), and the stopping-power ratio (SPR) was measured with a multi-layer ionization chamber. The DECT images were used to calculate SPR-, density- and elemental composition maps, and to assign voxel-wise compositions from a selection of human tissues. For a more comprehensive set of reference compositions, the original selection was extended by 135 additional tissues, corresponding to spongiosa, high-density bones and low-density tissues.Results.The root-mean-square error for the soft tissue carbon and oxygen content was 8.5 wt% and 9.5 wt% relative to the CCA result and 2.1 wt% and 10.3 wt% relative to the PG result. The phosphorous and calcium content were predicted within 0.4 wt% and 1.1 wt% of the CCA results, respectively. The largest discrepancies were encountered in samples whose composition deviated the most from tabulated compositions or that were more inhomogeneous.Significance.Overall, DECT-based composition estimations of relevant elements were in equal or better agreement with the ground truth than the established SECT-approach and could contribute toin-vivodose verification measurements.
Assuntos
Imagens de Fantasmas , Tomografia Computadorizada por Raios X , Animais , Suínos , Humanos , Método de Monte CarloRESUMO
Objective.We aim to: (1) quantify the benefits of lung sparing using non-adaptive magnetic resonance guided stereotactic body radiotherapy (MRgSBRT) with advanced motion management for peripheral lung cancers compared to conventional x-ray guided SBRT (ConvSBRT); (2) establish a practical decision-making guidance metric to assist a clinician in selecting the appropriate treatment modality.Approach.Eleven patients with peripheral lung cancer who underwent breath-hold, gated MRgSBRT on an MR-guided linear accelerator (MR linac) were studied. Four-dimensional computed tomography (4DCT)-based retrospective planning using an internal target volume (ITV) was performed to simulate ConvSBRT, which were evaluated against the original MRgSBRT plans. Metrics analyzed included planning target volume (PTV) coverage, various lung metrics and the generalized equivalent unform dose (gEUD). A dosimetric predictor for achievable lung metrics was derived to assist future patient triage across modalities.Main results.PTV coverage was high (median V100% > 98%) and comparable for both modalities. MRgSBRT had significantly lower lung doses as measured by V20 (median 3.2% versus 4.2%), mean lung dose (median 3.3 Gy versus 3.8 Gy) and gEUD. Breath-hold, gated MRgSBRT resulted in an average reduction of 47% in PTV volume and an average increase of 19% in lung volume. Strong correlation existed between lung metrics and the ratio of PTV to lung volumes (RPTV/Lungs) for both modalities, indicating that RPTV/Lungsmay serve as a good predictor for achievable lung metrics without the need for pre-planning. A threshold value of RPTV/Lungs< 0.035 is suggested to achieve V20 < 10% using ConvSBRT. MRgSBRT should otherwise be considered if the threshold cannot be met.Significance.The benefits of lung sparing using MRgSBRT were quantified for peripheral lung tumors; RPTV/Lungswas found to be an effective predictor for achievable lung metrics across modalities. RPTV/Lungscan assist a clinician in selecting the appropriate modality without the need for labor-intensive pre-planning, which has significant practical benefit for a busy clinic.
Assuntos
Tomografia Computadorizada Quadridimensional , Neoplasias Pulmonares , Pulmão , Imageamento por Ressonância Magnética , Radiocirurgia , Dosagem Radioterapêutica , Planejamento da Radioterapia Assistida por Computador , Humanos , Radiocirurgia/métodos , Neoplasias Pulmonares/radioterapia , Neoplasias Pulmonares/diagnóstico por imagem , Planejamento da Radioterapia Assistida por Computador/métodos , Imageamento por Ressonância Magnética/métodos , Pulmão/diagnóstico por imagem , Estudos Retrospectivos , Tomografia Computadorizada Quadridimensional/métodos , Masculino , Feminino , Radioterapia Guiada por Imagem/métodos , Suspensão da Respiração , Idoso , Pessoa de Meia-Idade , Tratamentos com Preservação do Órgão/métodos , Órgãos em RiscoRESUMO
AGuIX, a novel gadolinium-based nanoparticle, has been deployed in a pioneering double-blinded Phase II clinical trial aiming to assess its efficacy in enhancing radiotherapy for tumor treatment. This paper moves towards this goal by analyzing AGuIX uptake patterns in 23 patients. A phantom was designed to establish the relationship between AGuIX concentration and longitudinal ( T 1 ) relaxation. A 3T MRI and MP2RAGE sequence were used to generate patient T 1 maps. AGuIX uptake in tumors was determined based on longitudinal relaxivity. AGuIX (or placebo) was administered to 23 patients intravenously at 100 mg/kg 1-5 hours pre-imaging. Each of 129 brain metastases across 23 patients were captured in T 1 maps and examined for AGuIX uptake and distribution. Inferred AGuIX recipients had average tumor uptakes between 0.012 and 0.17 mg/ml, with a mean of 0.055 mg/ml. Suspected placebo recipients appeared to have no appreciable uptake. Tumors presented with varying spatial AGuIX uptake distributions, suspected to be related to differences in accumulation time and patient-specific bioaccumulation factors. This research demonstrates AGuIX's ability to accumulate in brain metastases, with quantifiable uptake via T 1 mapping. Future analyses will extend these methods to complete clinical trial data (~ 134 patients) to evaluate the potential relationship between nanoparticle uptake and possible tumor response following radiotherapy.Clinical Trial Registration Number: NCT04899908.Clinical Trial Registration Date: 25/05/2021.
Assuntos
Neoplasias Encefálicas , Gadolínio , Imageamento por Ressonância Magnética , Humanos , Neoplasias Encefálicas/secundário , Neoplasias Encefálicas/diagnóstico por imagem , Neoplasias Encefálicas/metabolismo , Neoplasias Encefálicas/tratamento farmacológico , Gadolínio/metabolismo , Gadolínio/administração & dosagem , Imageamento por Ressonância Magnética/métodos , Feminino , Pessoa de Meia-Idade , Masculino , Nanopartículas/química , Meios de Contraste/farmacocinética , Imagens de Fantasmas , Idoso , Adulto , Método Duplo-CegoRESUMO
Magnetic Resonance Imaging (MRI) is increasingly being used in treatment planning due to its superior soft tissue contrast, which is useful for tumor and soft tissue delineation compared to computed tomography (CT). However, MRI cannot directly provide mass density or relative stopping power (RSP) maps, which are required for calculating proton radiotherapy doses. Therefore, the integration of artificial intelligence (AI) into MRI-based treatment planning to estimate mass density and RSP directly from MRI has generated significant interest. A deep learning (DL) based framework was developed to establish a voxel-wise correlation between MR images and mass density as well as RSP. To facilitate the study, five tissue substitute phantoms were created, representing different tissues such as skin, muscle, adipose tissue, 45% hydroxyapatite (HA), and spongiosa bone. The composition of these phantoms was based on information from ICRP reports. Additionally, two animal tissue phantoms, simulating pig brain and liver, were prepared for DL training purposes. The phantom study involved the development of two DL models. The first model utilized clinical T1 and T2 MRI scans as input, while the second model incorporated zero echo time (ZTE) MRI scans. In the patient application study, two more DL models were trained: one using T1 and T2 MRI scans as input, and another model incorporating synthetic dual-energy computed tomography (sDECT) images to provide accurate bone tissue information. The DECT empirical model was used as a reference to evaluate the proposed models in both phantom and patient application studies. The DECT empirical model was selected as the reference for evaluating the proposed models in both phantom and patient application studies. In the phantom study, the DL model based on T1, and T2 MRI scans demonstrated higher accuracy in estimating mass density and RSP for skin, muscle, adipose tissue, brain, and liver. The mean absolute percentage errors (MAPE) were 0.42%, 0.14%, 0.19%, 0.78%, and 0.26% for mass density, and 0.30%, 0.11%, 0.16%, 0.61%, and 0.23% for RSP, respectively. The DL model incorporating ZTE MRI further improved the accuracy of mass density and RSP estimation for 45% HA and spongiosa bone, with MAPE values of 0.23% and 0.09% for mass density, and 0.19% and 0.07% for RSP, respectively. These results demonstrate the feasibility of using an MRI-only approach combined with DL methods for mass density and RSP estimation in proton therapy treatment planning. By employing this approach, it is possible to obtain the necessary information for proton radiotherapy directly from MRI scans, eliminating the need for additional imaging modalities.
Assuntos
Aprendizado Profundo , Imageamento por Ressonância Magnética , Imagens de Fantasmas , Terapia com Prótons , Imageamento por Ressonância Magnética/métodos , Terapia com Prótons/métodos , Humanos , Animais , Suínos , Planejamento da Radioterapia Assistida por Computador/métodos , Tomografia Computadorizada por Raios X/métodos , Dosagem RadioterapêuticaRESUMO
With the availability of MRI linacs, online adaptive intensity modulated radiotherapy (IMRT) has become a treatment option for liver cancer patients, often combined with hypofractionation. Intensity modulated proton therapy (IMPT) has the potential to reduce the dose to healthy tissue, but it is particularly sensitive to changes in the beam path and might therefore benefit from online adaptation. This study compares the normal tissue complication probabilities (NTCPs) for liver and duodenal toxicity for adaptive and non-adaptive IMRT and IMPT treatments of liver cancer patients. Adaptive and non-adaptive IMRT and IMPT plans were optimized to 50 Gy (RBE = 1.1 for IMPT) in five fractions for 10 liver cancer patients, using the original MRI linac images and physician-drawn structures. Three liver NTCP models were used to predict radiation-induced liver disease, an increase in albumin-bilirubin level, and a Child-Pugh score increase of more than 2. Additionally, three duodenal NTCP models were used to predict gastric bleeding, gastrointestinal (GI) toxicity with grades >3, and duodenal toxicity grades 2-4. NTCPs were calculated for adaptive and non-adaptive IMRT and IMPT treatments. In general, IMRT showed higher NTCP values than IMPT and the differences were often significant. However, the differences between adaptive and non-adaptive treatment schemes were not significant, indicating that the NTCP benefit of adaptive treatment regimens is expected to be smaller than the expected difference between IMRT and IMPT.
RESUMO
OBJECTIVE: Mapping CT number to material property dominates the proton range uncertainty. This work aims to develop a physics-constrained deep learning-based multimodal imaging (PDMI) framework to integrate physics, deep learning, MRI, and advanced dual-energy CT (DECT) to derive accurate patient mass density maps. METHODS: Seven tissue substitute MRI phantoms were used for validation including adipose, brain, muscle, liver, skin, spongiosa, 45% hydroxyapatite (HA) bone. MRI images were acquired using T1 weighted Dixon and T2 weighted short tau inversion recovery sequences. Training inputs are from MRI and twin-beam dual-energy images acquired at 120 kVp with gold/tin filters. The feasibility investigation included an empirical model and four residual networks (ResNet) derived from different training inputs and strategies by PDMI framework. PRN-MR-DE and RN-MR-DE denote ResNet (RN) trained with and without a physics constraint (P) using MRI (MR) and DECT (DE) images. PRN-DE stands for RN trained with a physics constraint using only DE images. A retrospective study using institutional patient data was also conducted to investigate the feasibility of the proposed framework. RESULTS: For the tissue surrogate study, PRN-MR-DE, PRN-DE, and RN-MR-DE result in mean mass density errors: -0.72%/2.62%/-3.58% for adipose; -0.03%/-0.61%/-0.18% for muscle; -0.58%/-1.36%/-4.86% for 45% HA bone. The retrospective patient study indicated that PRN-MR-DE predicted the densities of soft tissue and bone within expected intervals based on the literature survey, while PRN-DE generated large density deviations. CONCLUSION: The proposed PDMI framework can generate accurate mass density maps using MRI and DECT images. The supervised learning can further enhance model efficacy, making PRN-MR-DE outperform RN-MR-DE. The patient investigation also shows that the framework can potentially improve proton range uncertainty with accurate patient mass density maps. ADVANCES IN KNOWLEDGE: PDMI framework is proposed for the first time to inform deep learning models by physics insights and leverage the information from MRI to derive accurate mass density maps.
Assuntos
Aprendizado Profundo , Terapia com Prótons , Humanos , Processamento de Imagem Assistida por Computador/métodos , Estudos Retrospectivos , Prótons , Tomografia Computadorizada por Raios X/métodos , Imagem Multimodal/métodos , Imageamento por Ressonância Magnética/métodos , ObesidadeRESUMO
Objective. Range uncertainty in proton therapy is an important factor limiting clinical effectiveness. Magnetic resonance imaging (MRI) can measure voxel-wise molecular composition and, when combined with kilovoltage CT (kVCT), accurately determine mean ionization potential (Im), electron density, and stopping power ratio (SPR). We aimed to develop a novel MR-based multimodal method to accurately determine SPR and molecular compositions. This method was evaluated in tissue-mimicking andex vivoporcine phantoms, and in a brain radiotherapy patient.Approach. Four tissue-mimicking phantoms with known compositions, two porcine tissue phantoms, and a brain cancer patient were imaged with kVCT and MRI. Three imaging-based values were determined: SPRCM(CT-based Multimodal), SPRMM(MR-based Multimodal), and SPRstoich(stoichiometric calibration). MRI was used to determine two tissue-specific quantities of the Bethe Bloch equation (Im, electron density) to compute SPRCMand SPRMM. Imaging-based SPRs were compared to measurements for phantoms in a proton beam using a multilayer ionization chamber (SPRMLIC).Main results. Root mean square errors relative to SPRMLICwere 0.0104(0.86%), 0.0046(0.45%), and 0.0142(1.31%) for SPRCM, SPRMM, and SPRstoich, respectively. The largest errors were in bony phantoms, while soft tissue and porcine tissue phantoms had <1% errors across all SPR values. Relative to known physical molecular compositions, imaging-determined compositions differed by approximately ≤10%. In the brain case, the largest differences between SPRstoichand SPRMMwere in bone and high lipids/fat tissue. The magnitudes and trends of these differences matched phantom results.Significance. Our MR-based multimodal method determined molecular compositions and SPR in various tissue-mimicking phantoms with high accuracy, as confirmed with proton beam measurements. This method also revealed significant SPR differences compared to stoichiometric kVCT-only calculation in a clinical case, with the largest differences in bone. These findings support that including MRI in proton therapy treatment planning can improve the accuracy of calculated SPR values and reduce range uncertainties.
Assuntos
Neoplasias Encefálicas , Terapia com Prótons , Animais , Suínos , Prótons , Tomografia Computadorizada por Raios X/métodos , Imagens de Fantasmas , Imageamento por Ressonância Magnética , Calibragem , Planejamento da Radioterapia Assistida por Computador/métodosRESUMO
PURPOSE: The role of biliary stents in image-guided localization for pancreatic cancer has been inconclusive. To date, stent accuracy has been largely evaluated against implanted fiducials on cone beam computed tomography. We aim to use magnetic resonance (MR) soft tissue as a direct reference to examine the geometric and dosimetric impacts of stent-based localization on the newly available MR linear accelerator. METHODS: Thirty pancreatic cancer patients (132 fractions) treated on our MR linear accelerator were identified to have a biliary stent. In our standard adaptive workflow, patients were set up to the target using soft tissue for image registration and structures were re-contoured on daily MR images. The original plan was then projected on treatment anatomy and dose predicted, followed by plan re-optimization and treatment delivery. These online predicted plans were soft tissue-based and served as reference plans. Retrospective image registration to the stent was performed offline to simulate stent-based localization and the magnitude of shifts was taken as the geometric accuracy of stent localization. New predicted plans were generated based on stent-alignment for dosimetric comparison. RESULTS: Shifts were within 3 mm for 90% of the cases (mean = 1.5 mm); however, larger shifts up to 7.2 mm were observed. Average PTV coverage dropped by 1.1% with a maximum drop of 26.8%. The mean increase in V35Gy was 0.15, 0.05, 0.02, and 0.02 cc for duodenum, stomach, small bowel and large bowel, respectively. Stent alignment was significantly worse for all metrics except for small bowel (p = 0.07). CONCLUSIONS: Overall discrepancy between stent- and soft tissue-alignment was modest; however, large discrepancies were observed for select cases. While PTV coverage loss may be compensated for by using a larger margin, the increase in dose to gastrointestinal organs at risk may limit the role of biliary stents in image-guided localization.
Assuntos
Neoplasias Pancreáticas , Radiocirurgia , Radioterapia Guiada por Imagem , Humanos , Radiocirurgia/métodos , Estudos Retrospectivos , Neoplasias Pancreáticas/diagnóstico por imagem , Neoplasias Pancreáticas/radioterapia , Neoplasias Pancreáticas/cirurgia , Stents , Espectroscopia de Ressonância Magnética , Planejamento da Radioterapia Assistida por Computador/métodos , Dosagem Radioterapêutica , Radioterapia Guiada por Imagem/métodos , Neoplasias PancreáticasRESUMO
BACKGROUND: Paramagnetic species such as O2 and free radicals can enhance T1 and T2 relaxation times. If the change in relaxation time is sufficiently large, the contrast will be generated in magnetic resonance images. Since radiation is known to be capable of altering the concentration of O2 and free radicals during water radiolysis, it may be possible for radiation to induce MR signal change. PURPOSE: We present the first reported instance of x-ray-induced MR signal changes in water phantoms and investigate potential paramagnetic relaxation enhancement mechanisms associated with radiation chemistry changes in oxygen and free radical concentrations. METHODS: Images of water and 10 mM coumarin phantoms were acquired on a 0.35 T MR-linac before, during, and after a dose delivery of 80 Gy using an inversion-recovery dual-echo sequence with water nullified. Radiation chemistry simulations of these conditions were performed to calculate changes in oxygen and free radical concentrations. Published relaxivity values were then applied to calculate the resulting T1 change, and analytical MR signal equations were used to calculate the associated signal change. RESULTS: Compared to pre-irradiation reference images, water phantom images taken during and after irradiation showed little to no change, while coumarin phantom images showed a small signal loss in the irradiated region with a contrast-to-noise ratio (CNR) of 1.0-2.5. Radiation chemistry simulations found oxygen depletion of -11 µM in water and -31 µM in coumarin, resulting in a T1 lengthening of 24 ms and 68 ms respectively, and a simulated CNR of 1.0 and 2.8 respectively. This change was consistent with observations in both direction and magnitude. Steady-state superoxide, hydroxyl, hydroperoxyl, and hydrogen radical concentrations were found to contribute less than 1 ms of T1 change. CONCLUSION: Observed radiation-induced MR signal changes were dominated by an oxygen depletion mechanism. Free radicals were concluded to play a minor secondary role under steady-state conditions. Future applications may include in vivo FLASH treatment verification but would require an MR sequence with a better signal-to-noise ratio and higher temporal resolution than the one used in this study.
Assuntos
Imageamento por Ressonância Magnética , Oxigênio , Razão Sinal-Ruído , Radicais Livres , Imageamento por Ressonância Magnética/métodos , Imagens de FantasmasRESUMO
Current MRI-guided adaptive radiotherapy (MRgART) workflows require fraction-specific electron and/or mass density maps, which are created by deformable image registration (DIR) between the simulation CT images and daily MR images. Manual density overrides may also be needed where DIR-produced results are inaccurate. This approach slows the adaptive radiotherapy workflow and introduces additional dosimetric uncertainties, especially in the presence of the magnetic field. This study investigated a method based on a conditional generative adversarial network (cGAN) with a multi-planar method to generate synthetic CT images from low-field MR images to improve efficiency in MRgART workflows for prostate cancer. Fifty-seven male patients, who received MRI-guided radiation therapy to the pelvis using the ViewRay MRIdian Linac, were selected. Forty-five cases were randomly assigned to the training cohort with the remaining twelve cases assigned to the validation/testing cohort. All patient datasets had a semi-paired DIR-deformed CT-sim image and 0.35T MR image acquired using a true fast imaging with steady-state precession (TrueFISP) sequence. Synthetic CT images were compared with deformed CT images to evaluate image quality and dosimetric accuracy. To evaluate the dosimetric accuracy of this method, clinical plans were recalculated on synthetic CT images in the MRIdian treatment planning system. Dose volume histograms for planning target volumes (PTVs) and organs-at-risk (OARs) and dose distributions using gamma analyses were evaluated. The mean-absolute-errors (MAEs) in CT numbers were 30.1 ± 4.2 HU, 19.6 ± 2.3 HU and 158.5 ± 26.0 HU for the whole pelvis, soft tissue, and bone, respectively. The peak signal-to-noise ratio was 35.2 ± 1.7 and the structural index similarity measure was 0.9758 ± 0.0035. The dosimetric difference was on average less than 1% for all PTV and OAR metrics. Plans showed good agreement with gamma pass rates of 99% and 99.9% for 1%/1 mm and 2%/2 mm, respectively. Our study demonstrates the potential of using synthetic CT images created with a multi-planar cGAN method from 0.35T MRI TrueFISP images for the MRgART treatment of prostate radiotherapy. Future work will validate the method in a large cohort of patients and investigate the limitations of the method in the adaptive workflow.
RESUMO
Objective. Proton therapy of cancer improves dose conformality to the target and sparing of surrounding healthy tissues compared to conventional photon treatments. However, proton therapy's advantage could be even larger if proton range uncertainties were reduced. Sources of range uncertainties include computed tomography treatment planning images and variations in patient anatomy and setup. To reduce range uncertainties, we have developed a system for real-timein vivorange monitoring. The system is based on spectroscopy of prompt gamma-rays emitted through proton-nuclear interactions during irradiation. We validated the performance of our prompt gamma-ray spectroscopy detector prototype using tissue-mimicking and porcine samples.Approach. Measurements were performed in water, four tissue-mimicking samples (spongiosa, muscle, adipose tissue, and cortical bone), and two porcine samples (liver and brain). A dose of 0.9 Gy was delivered to a target at a depth of 12.5-17.5 cm. Multi-layer ionization chamber measurements were performed to determine stopping power ratios relative to water and ground truth proton ranges. Ground truth elemental compositions were determined using combustion analysis. Proton ranges and elemental compositions measured using prompt gamma-ray spectroscopy were compared to the ground truth.Main results. For all samples, the mean measured range over all pencil-beam spots differed from the ground truth by less than 1.2 mm. The mean standard deviation was 0.9 mm (range: 0.4-1.6 mm). The mean difference between ground truth and measured elemental compositions was 0.06gcm3(range: 0.00gcm3to 0.12gcm3).Significance. We verified the performance of our prompt gamma-ray spectroscopy detector prototype for proton range verification using tissue-mimicking and porcine samples. Measured proton ranges and elemental sample compositions were in good agreement with the ground truth. These measurements confirm the system's reliability for a variety of tissues and bridge the gap between previously-reported experiments and ongoingin vivopatient measurements.
Assuntos
Terapia com Prótons , Animais , Imagens de Fantasmas , Terapia com Prótons/métodos , Prótons , Planejamento da Radioterapia Assistida por Computador/métodos , Reprodutibilidade dos Testes , Análise Espectral , Suínos , Água/químicaRESUMO
Currently, adaptive strategies require time- and resource-intensive manual structure corrections. This study compares different strategies: optimization without manual structure correction, adaptation with physician-drawn structures, and no adaptation. Strategies were compared for 16 patients with pancreas, liver, and head and neck (HN) cancer with 1-5 repeated images during treatment: 'reference adaptation', with structures drawn by a physician; 'single-DIR adaptation', using a single set of deformably propagated structures; 'multi-DIR adaptation', using robust planning with multiple deformed structure sets; 'conservative adaptation', using the intersection and union of all deformed structures; 'probabilistic adaptation', using the probability of a voxel belonging to the structure in the optimization weight; and 'no adaptation'. Plans were evaluated using reference structures and compared using a scoring system. The reference adaptation with physician-drawn structures performed best, and no adaptation performed the worst. For pancreas and liver patients, adaptation with a single DIR improved the plan quality over no adaptation. For HN patients, integrating structure uncertainties brought an additional benefit. If resources for manual structure corrections would prevent online adaptation, manual correction could be replaced by a fast 'plausibility check', and plans could be adapted with correction-free adaptation strategies. Including structure uncertainties in the optimization has the potential to make online adaptation more automatable.
RESUMO
BACKGROUND: Megavoltage computed tomography (MVCT) has been implemented on many radiotherapy treatment machines for on-board anatomical visualization, localization, and adaptive dose calculation. Implementing an MR-only workflow by synthesizing MVCT from magnetic resonance imaging (MRI) would offer numerous advantages for treatment planning and online adaptation. PURPOSE: In this work, we sought to synthesize MVCT (sMVCT) datasets from MRI using deep learning to demonstrate the feasibility of MRI-MVCT only treatment planning. METHODS: MVCTs and T1-weighted MRIs for 120 patients treated for head-and-neck cancer were retrospectively acquired and co-registered. A deep neural network based on a fully-convolutional 3D U-Net architecture was implemented to map MRI intensity to MVCT HU. Input to the model were volumetric patches generated from paired MRI and MVCT datasets. The U-Net was initialized with random parameters and trained on a mean absolute error (MAE) objective function. Model accuracy was evaluated on 18 withheld test exams. sMVCTs were compared to respective MVCTs. Intensity-modulated volumetric radiotherapy (IMRT) plans were generated on MVCTs of four different disease sites and compared to plans calculated onto corresponding sMVCTs using the gamma metric and dose-volume-histograms (DVHs). RESULTS: MAE values between sMVCT and MVCT datasets were 93.3 ± 27.5, 78.2 ± 27.5, and 138.0 ± 43.4 HU for whole body, soft tissue, and bone volumes, respectively. Overall, there was good agreement between sMVCT and MVCT, with bone and air posing the greatest challenges. The retrospective dataset introduced additional deviations due to sinus filling or tumor growth/shrinkage between scans, differences in external contours due to variability in patient positioning, or when immobilization devices were absent from diagnostic MRIs. Dose distributions of IMRT plans evaluated for four test cases showed close agreement between sMVCT and MVCT images when evaluated using DVHs and gamma dose metrics, which averaged to 98.9 ± 1.0% and 96.8 ± 2.6% analyzed at 3%/3 mm and 2%/2 mm, respectively. CONCLUSIONS: MVCT datasets can be generated from T1-weighted MRI using a 3D deep convolutional neural network with dose calculation on a sample sMVCT in close agreement with the MVCT. These results demonstrate the feasibility of using MRI-derived sMVCT in an MR-only treatment planning workflow.
Assuntos
Planejamento da Radioterapia Assistida por Computador , Tomografia Computadorizada por Raios X , Humanos , Imageamento por Ressonância Magnética/métodos , Redes Neurais de Computação , Dosagem Radioterapêutica , Planejamento da Radioterapia Assistida por Computador/métodos , Estudos Retrospectivos , Tomografia Computadorizada por Raios X/métodosRESUMO
Proton therapy requires accurate dose calculation for treatment planning to ensure the conformal doses are precisely delivered to the targets. The conversion of CT numbers to material properties is a significant source of uncertainty for dose calculation. The aim of this study is to develop a physics-informed deep learning (PIDL) framework to derive accurate mass density and relative stopping power maps from dual-energy computed tomography (DECT) images. The PIDL framework allows deep learning (DL) models to be trained with a physics loss function, which includes a physics model to constrain DL models. Five DL models were implemented including a fully connected neural network (FCNN), dual-FCNN (DFCNN), and three variants of residual networks (ResNet): ResNet-v1 (RN-v1), ResNet-v2 (RN-v2), and dual-ResNet-v2 (DRN-v2). An artificial neural network (ANN) and the five DL models trained with and without physics loss were explored to evaluate the PIDL framework. Two empirical DECT models were implemented to compare with the PIDL method. DL training data were from CIRS electron density phantom 062M (Computerized Imaging Reference Systems, Inc., Norfolk, VA). The performance of DL models was tested by CIRS adult male, adult female, and 5-year-old child anthropomorphic phantoms. For density map inference, the physics-informed RN-v2 was 3.3%, 2.9% and 1.9% more accurate than ANN for the adult male, adult female, and child phantoms. The physics-informed DRN-v2 was 0.7%, 0.6%, and 0.8% more accurate than DRN-v2 without physics training for the three phantoms, respectfully. The results indicated that physics-informed training could reduce uncertainty when ANN/DL models without physics training were insufficient to capture data structures or derived significant errors. DL models could also achieve better image noise control compared to the empirical DECT parametric mapping methods. The proposed PIDL framework can potentially improve proton range uncertainty by offering accurate material properties conversion from DECT.
Assuntos
Aprendizado Profundo , Terapia com Prótons , Pré-Escolar , Feminino , Humanos , Masculino , Imagens de Fantasmas , Física , Tomografia Computadorizada por Raios X/métodosRESUMO
Objective. Kilovoltage computed tomography (kVCT) is the cornerstone of radiotherapy treatment planning for delineating tissues and towards dose calculation. For the former, kVCT provides excellent contrast and signal-to-noise ratio. For the latter, kVCT may have greater uncertainty in determining relative electron density (ρe) and proton stopping power ratio (SPR). Conversely, megavoltage CT (MVCT) may result in superior dose calculation accuracy. The purpose of this work was to convert kVCT HU to MVCT HU using deep learning to obtain higher accuracyρeand SPR.Approach. Tissue-mimicking phantoms were created to compare kVCT- and MVCT-determinedρeand SPR to physical measurements. Using 100 head-and-neck datasets, an unpaired deep learning model was trained to learn the relationship between kVCTs and MVCTs, creating synthetic MVCTs (sMVCTs). Similarity metrics were calculated between kVCTs, sMVCTs, and MVCTs in 20 test datasets. An anthropomorphic head phantom containing bone-mimicking material with known composition was scanned to provide an independent determination ofρeand SPR accuracy by sMVCT.Main results. In tissue-mimicking bone,ρeerrors were 2.20% versus 0.19% and SPR errors were 4.38% versus 0.22%, for kVCT versus MVCT, respectively. Compared to MVCT,in vivomean difference (MD) values were 11 and 327 HU for kVCT and 2 and 3 HU for sMVCT in soft tissue and bone, respectively.ρeMD decreased from 1.3% to 0.35% in soft tissue and 2.9% to 0.13% in bone, for kVCT and sMVCT, respectively. SPR MD decreased from 1.8% to 0.24% in soft tissue and 6.8% to 0.16% in bone, for kVCT and sMVCT, respectively. Relative to physical measurements,ρeand SPR error in anthropomorphic bone decreased from 7.50% and 7.48% for kVCT to <1% for both MVCT and sMVCT.Significance. Deep learning can be used to map kVCT to sMVCT, suggesting higher accuracyρeand SPR is achievable with sMVCT versus kVCT.
Assuntos
Terapia com Prótons , Prótons , Elétrons , Aprendizado de Máquina , Imagens de Fantasmas , Planejamento da Radioterapia Assistida por ComputadorRESUMO
Objective.In MRI-based radiation therapy planning, mitigating patient-specific distortion with standard high bandwidth scans can result in unnecessary sacrifices of signal to noise ratio. This study investigates a technique for distortion detection and mitigation on a patient specific basis.Approach.Fast B0 mapping was performed using a previously developed technique for high-resolution, large dynamic range field mapping without the need for phase unwrapping algorithms. A phantom study was performed to validate the method. Distortion mitigation was validated by reducing geometric distortion with increased acquisition bandwidth and confirmed by both the B0 mapping technique and manual measurements. Images and contours from 25 brain stereotactic radiosurgery patients and 95 targets were analyzed to estimate the range of geometric distortions expected in the brain and to estimate bandwidth required to keep all treatment targets within the ±0.5 mm iso-distortion contour.Main Results.The phantom study showed, at 3 T, the technique can measure distortions with a mean absolute error of 0.12 mm (0.18 ppm), and a maximum error of 0.37 mm (0.6 ppm). For image acquisition at 3 T and 1.0 mm resolution, mean absolute distortion under 0.5 mm in patients required bandwidths from 109 to 200 Hz px-1for patients with the least and most distortion, respectively. Maximum absolute distortion under 0.5 mm required bandwidths from 120 to 390 Hz px-1.Significance.The method for B0 mapping was shown to be valid and may be applied to assess distortion clinically. Future work will adapt the readout bandwidth to prospectively mitigate distortion with the goal to improve radiosurgery treatment outcomes by reducing healthy tissue exposure.
Assuntos
Radiocirurgia , Algoritmos , Encéfalo , Humanos , Imageamento por Ressonância Magnética/métodos , Imagens de Fantasmas , Radiocirurgia/métodosRESUMO
PURPOSE: To perform a propensity-score matched analysis comparing stereotactic body radiation therapy (SBRT) boost and high-dose-rate (HDR) boost for localized prostate cancer. METHODS AND MATERIALS: A single-institution retrospective chart review was conducted of men treated with pelvic external beam radiation therapy (EBRT) and SBRT boost (21 Gy and 19 Gy in 2 fractions) to the prostate for prostate cancer. A cohort treated at the same institution with HDR brachytherapy boost (19 Gy in 2 fractions) was compared. Propensity-score (PS) matching and multivariable Cox regression were used for analysis. Outcomes were biochemical recurrence freedom (BCRF) and metastasis freedom (MF). RESULTS: One hundred thirty-one men were treated with SBRT boost and 101 with HDR boost with median follow-up of 73.4 and 186.0 months, respectively. In addition, 68.8% of men had high-risk and 26.0% had unfavorable-intermediate disease, and 94.3% received androgen deprivation therapy. Five- and 10-year unadjusted BCRF was 88.8% and 85.3% for SBRT and 91.8% and 74.6% for HDR boost (log-rank P = .3), and 5- and 10-year unadjusted MF was 91.7% and 84.3% for SBRT and 95.8% and 82.0% for HDR (log-rank P = .8). After adjusting for covariates, there was no statistically significant difference in BCRF (hazard ratio [HR] 0.81; 95% confidence interval [CI], 0.37-1.79; P = .6) or MF (HR 1.07; 95% CI, 0.44-2.57; P = .9) between SBRT and HDR boost. Similarly, after PS matching, there was no statistically significant difference between SBRT and HDR (BCRF: HR 0.66, 0.27-1.62, P = .4; MF: HR 0.84, 0.31-2.26, P = .7). Grade 3+ genitourinary and gastrointestinal toxicity in the SBRT cohort were 4.6% and 1.5%, and 3.0% and 0.0% in the HDR cohorts (P = .4, Fisher exact test). CONCLUSIONS: SBRT boost plus pelvic EBRT for prostate cancer resulted in similar BCRF and MF to HDR boost in this single institution, PS matched retrospective analysis. Toxicity was modest. Prospective evaluation of SBRT boost for the treatment of unfavorable-intermediate and high-risk prostate cancer is warranted.
Assuntos
Braquiterapia/métodos , Neoplasias da Próstata/radioterapia , Radiocirurgia , Radioterapia de Intensidade Modulada , Idoso , Antagonistas de Androgênios/uso terapêutico , Anilidas/uso terapêutico , Braquiterapia/efeitos adversos , Estudos de Coortes , Terapia Combinada/métodos , Intervalos de Confiança , Fracionamento da Dose de Radiação , Humanos , Leuprolida/uso terapêutico , Masculino , Pessoa de Meia-Idade , Nitrilas/uso terapêutico , Pontuação de Propensão , Antígeno Prostático Específico/sangue , Neoplasias da Próstata/sangue , Neoplasias da Próstata/tratamento farmacológico , Neoplasias da Próstata/patologia , Radiocirurgia/efeitos adversos , Radioterapia de Intensidade Modulada/efeitos adversos , Análise de Regressão , Estudos Retrospectivos , Compostos de Tosil/uso terapêuticoRESUMO
PURPOSE: Proton therapy is becoming an increasingly popular cancer treatment modality due to the proton's physical advantage in that it deposits the majority of its energy at the distal end of its track where the tumor is located. The proton range in a material is determined from the stopping power ratio (SPR) of the material. However, SPR is typically estimated based on a computed tomography (CT) scan which can lead to range estimation errors due to the difference in x-ray and proton interactions in matter, which can preclude the ability to utilize protons to their full potential. Applications of magnetic resonance imaging (MRI) in radiotherapy have increased over the past decade and using MRI to calculate SPR directly could provide numerous advantages. The purpose of this study was to develop a practical implementation of a novel multimodal imaging method for estimating SPR and compare the results of this method to physical measurements in which values were computed directly using tissue substitute materials fabricated to mimic skin, muscle, adipose, and spongiosa bone. METHODS: For both the multimodal imaging method and physical measurements, SPR was calculated using the Bethe-Bloch equation from values of relative electron density and mean ionization potential determined for each tissue. Parameters used to estimate SPR using the multimodal imaging method were extracted from Dixon water-only and (1 H) proton density-weighted zero echo time MRI sequences and CT, with both kVCT and MVCT used separately to evaluate the performance of each. For comparison, SPR was also computed from kVCT using the stoichiometric method, the current clinical standard. RESULTS: Results showed that our multimodal imaging approach using MRI with either kVCT or MVCT was in close agreement to SPR calculated from physical measurements for the four tissue substitutes evaluated. Using MRI and MVCT, SPR values estimated using our method were within 1% of physical measurements and were more accurate than the stoichiometric method for the tissue types studied. CONCLUSIONS: We have demonstrated the methodology for improved estimation of SPR using the proposed multimodal imaging framework.