Physics in medicine and biology最新文献

•Objective:In this study, we explore positron emission tomography(PET)/magnetic resonance imaging (MRI) joint reconstruction within a deeplearning (DL) framework, introducing a novel synergistic method. •Approach:We propose a new approach based on a variational autoencoder (VAE)constraint combined with the alternating direction method of multipliers (ADMM)optimization technique. We compare several VAE architectures, including jointVAE, mixture of experts (MoE) and product of experts (PoE), to determine theoptimal latent representation for the two modalities. We trained then evaluatedthe architectures on a brain PET/MRI dataset. •Main results:We showed that our approach takes advantage of each modalitysharing information to each other, which results in improved peak signal-to-noiseratio (PSNR) and structural similarity (SSIM) as compared with traditionalreconstruction methods, particularly for short acquisition times. We find that theone particular architecture, MMJSD, is the most effective for our methodology. •Significance:The proposed method outperforms classical approaches especiallyin noisy and undersampled conditions by making use of the two modalities together to compensate for the missing information.

Objective. Bragg peak measurements play a key role in the beam quality assurance in proton therapy. Used as base data for the treatment planning softwares, the accuracy of the data is crucial when defining the range of the protons in the patient.Approach. In this paper a protocol to reconstruct a Pristine Bragg Peak exploring the direct correlation between the particle flux and the dose deposited by particles is presented. Proton flux measurements at the HollandPTC and FLUKA Monte Carlo simulations are used for this purpose. This new protocol is applicable to plastic scintillator detectors developed for Quality Assurance applications. In order to obtain the Bragg curve using a plastic fiber detector, a PMMA phantom with a decoupled and moveable stepper was designed. The step phantom allows to change the depth of material in front of the fiber detector during irradiations. The Pristine Bragg Peak reconstruction protocol uses the measured flux of particles at each position and multiplies it by the average dose obtained from the Monte Carlo simulation at each position.Main results. The results show that with this protocol it is possible to reconstruct the Bragg Peak with an accuracy of about 470µm, which is in accordance with the tolerances set by the AAPM.Significance. It has the advantage to be able to overcome the quenching problem of scintillators in the high ionization density region of the Bragg peak.

Objective. To compare in reproducible and equalized conditions the performance of two independent proton range verification systems based on prompt gamma-ray detectors from two different proton therapy centers.Approach. An anthropomorphic head phantom with calibrated stopping power, serving as ground truth, was irradiated with comparable treatment plans, spot positions and energies in both facilities. Clinical beam current, tumor contour and dose were used. The absolute range measurement was compared to the expected value according to the ground truth. The statistical precision was assessed by repeating each measurement ten times. Sensitivity to relative range shifts was evaluated by introducing 2 mm and 5 mm plastic slabs on half of the field.Main results. The resulting absolute range accuracy was within 2.4 mm in all cases. Relative range shifts were detected with deviations lower than 14%.Significance. The performance of both systems was deemed worthy of clinical application for the detection of range deviations. This study represents the first comparison of independent prompt gamma-ray-based proton range verification systems under equalized conditions with realistic treatment fields and beam currents.

Objective. Image-guided and adaptive proton therapy rely on daily CBCT or CT imaging, which increases radiation dose and radiation-induced cancer risk. Online adaptation however also reduces setup uncertainty, and the additional risk might be compensated by reducing the setup robustness margin. This work developed a framework to investigate how much this robustness margin should be reduced to offset the secondary cancer risk from additional imaging dose and applied it to proton therapy for head-and-neck cancer.Approach. For five patients with head-and-neck cancer, voxel-wise CT and CBCT imaging doses were estimated with Monte Carlo radiation transport simulations, calibrated with air and PMMA phantom measurements. The total dose of several image-guided and adaptive treatments protocols was calculated by summing the planning CT dose, daily and weekly CBCT or CT dose, and therapy dose, robustly optimized with setup margins between 0 and 4 mm. These were compared to a reference protocol with 4 mm setup margin without daily imaging. All plans further used 3% range robustness. Organ-wise excess absolute risk (EAR) of cancer was calculated with three models to determine at which setup margin the total EAR of image-guided and adaptive treatment protocols was equal to the total EAR of the reference.Results. The difference between the simulated and measured CT and CBCT doses was within 10%. Using the Monte Carlo models, we found that a 1 mm setup margin reduction was sufficient for most patients, treatment protocols, and cancer risk models to compensate the additional risk imposed by daily and weekly imaging. For some protocols, even a smaller reduction sufficed, depending on the imaging frequency and type. The risk reduction by reducing the margin was mainly due to reducing the risk for carcinomas in the brain and, for some patients, the oral cavity.Significance. Our framework allows to compare an increased imaging dose with the reduced treatment dose from margin reductions in terms of radiation-induced cancer risk. It is extendable to different treatment sites, modalities, and imaging protocols, in clinic-specific or even patient-specific assessments.

Objective.Deep learning (DL) is becoming increasingly important in generating attenuation maps for accurate attenuation correction (AC) in cardiac perfusion SPECT imaging. Typically, DL models take inputs from initial reconstructed SPECT images, which are performed on the photopeak window and often also on scatter windows. While prior studies have demonstrated improvements in DL performance when scatter window images are incorporated into the DL input, the comprehensive analysis of the impact of employing different scatter windows remains unassessed. Additionally, existing research mainly focuses on applying DL to SPECT scans obtained at clinical standard count levels. This study aimed to assess utilities of DL from two aspects: (1) investigating the impact when different scatter windows were used as input to DL, and (2) evaluating the performance of DL when applied on SPECT scans acquired at a reduced count level.Approach.We utilized 1517 subjects, with 386 subjects for testing and the remaining 1131 for training and validation.Main results.The results showed that as scatter window width increased from 4% to 30%, a slight improvement was observed in DL estimated attenuation maps. The application of DL models to quarter-count (¼-count) SPECT scans, compared to full-count scans, showed a slight reduction in performance. Nonetheless, discrepancies across different scatter window configurations and between count levels were minimal, with all normalized mean square error (NMSE) values remaining within 2.1% when comparing the different DL attenuation maps to the reference CT maps. For attenuation corrected SPECT slices using DL estimated maps, NMSE values were within 0.5% when compared to CT correction.Significance.This study, leveraging an extensive clinical dataset, showed that the performance of DL seemed to be consistent across the use of varied scatter window settings. Moreover, our investigation into reduced count studies indicated that DL could provide accurate AC even at a ¼-count level.

Objective: Linear energy transfer (LET) verification was conducted using a silicon-on-insulator (SOI) microdosimeter during the commissioning of LET-optimized carbon-ion radiotherapy. This advanced treatment technique is expected to improve local control rates, especially in hypoxic tumors. Approach: An SOI microdosimeter with a cylindrical sensitive volume of 30 μm diameter and 5 μm thickness was used. Simple cubic plans and patient plans using the carbon-ion beams were created by treatment planning system, and the calculated LETd values were compared with the measured LETd values obtained by the SOI microdosimeter. Main results: Reasonable agreement between the measured and calculated LETd was seen in the plateau region of depth LETd profile, whereas the measured LETd were below the calculated LETd in the peak region, specifically where LETd exceeds 75 keV/μm. The discrepancy in the peak region may arise from the uncertainties in the calibration process of the SOI microdosimeter. Excluding the peak region, the average ratio and standard deviation between measured and calculated LETd values were 0.996 and 7%, respectively. Significance: This verification results in the initiation of clinical trials for LET-optimized carbon-ion radiotherapy at QST Hospital, National Institutes for Quantum Science and Technology.

Objective: Spot size stability is crucial in pencil beam scanning (PBS) proton therapy, and variations in spot size can disrupt dose distributions. Recently, a novel proton beam delivery method known as DynamicARC PBS scanning has been introduced. The current study investigates the dosimetric impact of spot size errors in DynamicARC proton therapy for head and neck (HNC), prostate, and lung cancers. Approach: Robustly optimized DynamicARC proton therapy plans were created for HNC (n=4), prostate (n=4), and lung (n=4) cancer patients. Spot size errors of ±10%, ±15%, and ±20% were introduced, and their effects on target coverage (D95% and D99%), homogeneity index (HI), and organ-at-risk (OAR) doses were analyzed across different cancer sites. Main Results: HNC and lung cancer plans showed greater vulnerability to spot size errors, with reductions in target coverage of up to 4.8% under -20% spot size errors. Dose homogeneity was also more affected in these cases, with HI degrading by 0.12 in lung cancer. Prostate cancer demonstrated greater resilience to spot size variations, even under errors of ±20%. For spot size errors ±10%, the oral cavity, parotid glands, and constrictor muscles experienced Dmean deviations within ±1.2%, while deviations were limited to ±0.5% for D10% of the bladder and rectum and ±0.3% for V20Gy(RBE) of the lungs. The robustness analysis indicated that lung cancer plans were most susceptible to robustness reductions caused by spot size errors, while HNC plans demonstrated moderate sensitivity. Conversely, prostate cancer plans demonstrated high robustness, experiencing only minimal reductions in target coverage. Significance: While the ±10% spot size tolerance is appropriate in the majority of the cases, lung cancer plans may require more stringent criteria. As DynamicARC becomes clinically available, measuring spot size errors in practice will be essential to validate these findings and refine tolerance thresholds for clinical use.

Cone beam CT is integral to most modern radiotherapy treatments. The application of daily and repeat CBCT imaging can lead to high imaging doses over a large volume of tissue that extends beyond the treatment site. Hence, it is important to ensure exposures are optimised to keep doses as low as reasonably achievable, whilst ensuring images are suitable for the clinical task. This IPEM topical report presents the results of the first UK survey of dose indices in radiotherapy CBCT. Dose measurements, as defined by the cone beam dose index (CBDI_w), were collected along with protocol information for seven treatment sites. Where a range of optimised protocols were available in a centre, a sample of patient data demonstrating the variation in protocol use were requested. Protocol CBDI_wvalues were determined from the average dosimetry data for each type of linear accelerator, and median CBDI_wand scan length were calculated for each treatment site at each centre. Median CBDI_wvalues were compared and summary statistics derived that enable the setting of national dose reference levels (DRLs). A total of 63 UK radiotherapy centres contributed data. The proposed CBDI_wDRLs are; prostate 20.6 mGy, gynaecological 20.8 mGy, breast 5.0 mGy, 3D-lung 6.0 mGy, 4D-lung 11.8 mGy, brain 3.5 mGy and head/neck 4.2 mGy. However, large differences between models of imaging system were noted. Where centres had pro-active optimisation strategies in place, such as sized based protocols with selection criteria, dose reductions on the 'average' patient were possible compared with vendor defaults. Optimisation of scan length was noted in some clinical sites, with Elekta users tending to fit different collimators for prostate imaging (relatively short) compared with gynaecological treatments (longest). This contrasts with most Varian users who apply the default scan length in most cases.

Objective.The highly conformal carbon-ion radiotherapy is associated with an increased sensitivity of the dose distributions to internal changes in the patient during the treatment course. Hence, monitoring methodologies capable of detecting such changes are of vital importance. We established experimental setup conditions to address the sensitivity of a monitoring approach based on secondary-fragment tracking for detecting clinically motivated air cavity dimensions in a homogeneous head-sized PMMA phantom in 40 mm depth.Approach.The air cavities were positioned within the entrance channel of a treatment field of 50 mm diameter at three lateral positions. The measured secondary-fragment emission profiles were compared to a reference measurement without cavities. The experiments were conducted at the Heidelberg Ion-Beam Therapy Center in Germany at typical doses and dose rates.Main results.Significances above a detectability threshold of 2σfor the larger cavities (20 mm diameter and 4 mm thickness, and 20 mm diameter and 2 mm thickness) across the entire treatment field. The smallest cavity of 10 mm diameter and 2 mm thickness, which is on the lower limit of clinical interest, could not be detected at any position. We also demonstrated that it is feasible to reconstruct the lateral position of the cavity on average within 2.8 mm, once the cavity is detected. This is sufficient for the clinicians to estimate medical effects of such a cavity and to decide about the need for a control imaging CT.Significance.This investigation defines well-controlled reference conditions for the evaluation of the performance of any kind of treatment monitoring method and its capability to detect internal changes within head-sized objects. Four air cavities with volumes between 0.31 cm³and 1.26 cm³were narrowed down around the detectability threshold of this secondary-fragment-based monitoring method.

Objective.This study describes the development, validation, and integration of a detector response model that accounts for the combined effects of x-ray crosstalk, charge sharing, and pulse pileup in photon-counting detectors.Approach.The x-ray photon transport was simulated using Geant4, followed by analytical charge sharing simulation in MATLAB. The analytical simulation models charge clouds with Gaussian-distributed charge densities, which are projected on a 3×3 pixel neighborhood of interaction location to compute detected counts. For pulse pileup, a prior analytical method for redistribution of energy-binned counts was implemented for delta pulses. The x-ray photon transport and charge sharing components were validated using experimental data acquired on the CdTe-based Pixirad-1/Pixie-III detector using monoenergetic beams at 26, 33, 37, and 50 keV. The pulse pileup implementation was verified with a comparable Monte Carlo simulation. The model output without pulse pileup was used to generate spatio-energetic response matrices for efficient simulation of scanner-specific photon-counting CT (PCCT) images with DukeSim, with pulse pileup modeled as a post-processing step on simulated projections. For analysis, images for the Gammex multi-energy phantom and the XCAT chest phantom were simulated at 120 kV, both with and without pulse pileup for a range of doses (27-1344 mAs). The XCAT images were evaluated qualitatively at 120 mAs, while images for the Gammex phantom were evaluated quantitatively for all doses using measurements of attenuation coefficients and Calcium concentrations.Main results.Reasonable agreement was observed between simulated and experimental spectra with Mean Absolute Percentage Error Values (MAPE) between 10%and 31%across all incident energies and detector modes. The increased pulse pileup from increased dose affected attenuation coefficients and calcium concentrations, with an effect on calcium quantification as high as MAPE of 28%.Significance.The presented approach demonstrates the viability of the model for enabling VITs to assess and optimize the clinical performance of PCCT.