Physics in medicine and biology最新文献

Objective.This paper introduces a novel PET imaging methodology called 3-dimensional positron imaging (3Dπ), which integrates total-body coverage, time-of-flight (TOF) technology, ultra-low dose imaging capabilities, and ultra-fast readout electronics inspired by emerging technology from the DarkSide collaboration.Approach.The study evaluates the performance of 3Dπusing Monte Carlo simulations based on NEMA NU 2-2018 protocols. The methodology employs a homogenous, monolithic scintillator composed of liquid argon (LAr) doped with xenon (Xe) with silicon photomultipliers (SiPMs) operating at cryogenic temperatures.Main results.Substantial improvements in system performance are observed, with the 3Dπsystem achieving a noise equivalent count rate of 3.2 Mcps at 17.3 kBq ml^-1, continuing to increase up to 4.3 Mcps at 40 kBq ml^-1. Spatial resolution measurements show an average FWHM of 2.7 mm across both axial positions. The system exhibits superior sensitivity, with values reaching 373 kcps MBq^-1with a line source at the center of the field of view. Additionally, 3Dπachieves a TOF resolution of 151 ps at 5.3 kBq ml^-1, highlighting its potential to produce high-quality images with reduced noise levels.Significance.The study underscores the potential of 3Dπin improving PET imaging performance, offering the potential for shorter scan times and reduced radiation exposure for patients. The Xe-doped LAr offers advantages such as fast scintillation, enhanced light yield, and cost-effectiveness. Future research will focus on optimizing system geometry and further refining reconstruction algorithms to exploit the strengths of 3Dπfor clinical applications.

Objective: In particle therapy (PT), several methods are being investigated to help reduce range margins and identify deviations from the original treatment plan, such as prompt-gamma (PG) imaging with Compton cameras (CC). To reconstruct the images, the Origin Ensemble (OE) algorithm is commonly used. In the context of PT, artifacts and strong noise often affect CC images. To improve the ability of OE to identify range shifts, and also to enhance image quality, we propose to regularize OE using beam a-priori knowledge (beam prior). Approach: We implemented the beam prior to OE using the class of Gibbs' distribution functions. For evaluation, Monte-Carlo simulations of centered and off-center beams with therapeutic energies impinging on a PMMA target were conducted in GATE. To introduce range shifts, air layers were introduced into the target. In addition, the effect of a bone layer, closer to a realistic scenario, was investigated. OE with the beam prior (BP-OE) and conventional OE (reference) were compared using the spill-over-ratio (SOR) as well as shifts in the distal falloff in projections using cubic splines with Chebyshev nodes. Main results: BP-OE improved the shift estimates by up to 11% compared to conventional OE for centered and up to 250% with off-centered beams. BP-OE decreased the image noise level, improving the SOR significantly by up to 96%. Significance: BP-OE applied to CC data can improve shift estimations compared to conventional OE. The developed Gibbs-based regularization framework also allows further prior functions to be included into OE, for instance, smoothing or edge-preserving priors. BP-OE could be extended to PET-based range verification or multiple-beam scenarios.

Objective: Utilizing prompt gammas in preclinical pinhole-collimated PET avoids image degradation due to positron range blurring and photon down scatter, enables multi-isotope PET and can improve counting statistics for low-abundance positron emitters. This was earlier reported for 124I, 89Zr and simultaneous 124I -18F PET using the VECTor scanner (MILabs, The Netherlands), demonstrating sub-mm resolution despite long positron ranges. The aim of the present study is to investigate if such sub-mm PET imaging is also feasible for a large variety of other isotopes including those with extremely high energy prompt gammas (>1 MeV) or with complex emission spectra of prompt gammas.

Approach: We use Monte Carlo simulations to assess achievable image resolutions and uniformity across a broad range of spectrum types and emitted prompt gamma energies (603 keV - 2.2 MeV), using 52Mn, 94Tc, 89Zr, 44Sc, 86Y, 72As, 124I, 38K, and 66Ga.

Main results: Our results indicate that sub-millimeter resolution imaging may be feasible for almost all isotopes investigated, with the currently used cluster pinhole collimators. At prompt gamma energies of 603 keV of 124I, an image resolution of ~0.65 mm was achieved, while for emissions at 703, 744, 834, and 909 keV of 94Tc, 52Mn, 72As, and 89Zr, respectively, ~0.7 mm resolution was obtained. Finally, at ultra-high energies of 1.2 (44Sc) and 1.4 MeV (52Mn) resolutions of ~0.75 mm and ~0.8 mm could still be achieved although ring artifacts were observed at the highest energies (1.4 MeV). For 38K (2.2 MeV), an image resolution of 1.2 mm was achieved utilizing its 2.2 MeV prompt emission.

Significance: This work shows that current cluster pinhole collimators are suitable for sub-mm resolution prompt PET up till at least 1.4 MeV. This may open up new avenues to developing new tracer applications and therapies utilizing these PET isotopes.

Objective: Clinical implementation of in-beam PET monitoring in proton therapy requires the integration of an online fast and reliable dose calculation engine. This manuscript reports on the achievement of real-time reconstruction of 3D dose and activity maps with proton range verification from experimental in-beam PET measurements. Approach: Several cylindrical homogeneous PMMA phantoms were irradiated with a monoenergetic 70-MeV proton beam in a clinical facility. Additionally, PMMA range-shifting foils of varying thicknesses were placed at the proximal surface of the phantom to investigate range shift prediction capabilities. PET activity was measured using a state-of-the-art in-house developed six-module PET scanner equipped with online PET reconstruction capabilities. For real-time dose estimation, we integrated this system with an in-beam dose estimation (IDE) algorithm, which combines a GPU-based 3D reconstruction algorithm with a dictionary-based software, capable of estimating deposited doses from the 3D PET activity images. The range shift prediction performance has been quantitatively studied in terms of the minimum dose to be delivered and the maximum acquisition time. Main results: With this framework, 3D dose maps were accurately reconstructed and displayed with a delay as short as one second. For a dose fraction of 8.4 Gy at the Bragg peak maximum, range shifts as small as 1 mm could be detected. The quantitative analysis shows that accumulating 20 seconds of statistics from the start of the irradiation, doses down to 1 Gy could be estimated online with total uncertainties smaller than 2 mm. Significance. The hardware and software combination employed in this work can deliver dose maps and accurately predict range shifts after short acquisition times and small doses, suggesting that real-time monitoring and dose reconstruction during proton therapy are within reach. Future work will focus on testing the methodology in more complex clinical scenarios and on upgrading the PET prototype for increased sensitivity. .

Objective: Radiotherapy planning requires significant expertise to balance tumor control and organ-at-risk (OAR) sparing. Automated planning can improve both efficiency and quality. This study introduces GPT-Plan, a novel multi-agent system powered by the GPT-4 family of large language models (LLMs), for automating the iterative radiotherapy plan optimization.Approach: GPT-Plan uses LLM-driven agents, mimicking the collaborative clinical workflow of a dosimetrist and physicist, to iteratively generate and evaluate text-based radiotherapy plans based on predefined criteria. Supporting tools assist the agents by leveraging historical plans, mitigating LLM hallucinations, and balancing exploration and exploitation. Performance was evaluated on 12 lung (IMRT) and 5 cervical (VMAT) cancer cases, benchmarked against the ECHO auto-planning method and manual plans. The impact of historical plan retrieval on efficiency was also assessed.Results: For IMRT lung cancer cases, GPT-Plan generated high-quality plans, demonstrating superior target coverage and homogeneity compared to ECHO while maintaining comparable or better OAR sparing. For VMAT cervical cancer cases, plan quality was comparable to a senior physicist and consistently superior to a junior physicist, particularly for OAR sparing. Retrieving historical plans significantly reduced the number of required optimization iterations for lung cases (p < 0.01) and yielded iteration counts comparable to those of the senior physicist for cervical cases (p=0.313). Occasional LLM hallucinations have been mitigated by self-reflection mechanisms. One limitation was the inaccuracy of vision-based LLMs in interpreting dose images.Significance: This pioneering study demonstrates the feasibility of automating radiotherapy planning using LLM-powered agents for complex treatment decision-making tasks. While challenges remain in addressing LLM limitations, ongoing advancements hold potential for further refining and expanding GPT-Plan's capabilities. .

Internationally, reference dosimetry for clinical proton beams largely follows the guidelines published by the International Atomic Energy Agency (IAEA TRS-398 Rev. 1 (2024). This approach yields a relative standard uncertainty of 1.7% (k= 1) on the absorbed dose to water determined under reference conditions. The new IPEM code of practice presented here, enables the relative standard uncertainty on the absorbed dose to water measured under reference conditions to be reduced to 1.0% (k= 1). This improvement is based on the absorbed dose to water calibration service for proton beams provided by the National Physical Laboratory (NPL), the UK's primary standards laboratory. This significantly reduced uncertainty is achieved through the use of a primary standard level graphite calorimeter to derive absorbed dose to water directly in the clinical department's beam. This eliminates the need for beam quality correction factors (kQ,Q0) as required by the IAEA TRS-398 approach. The portable primary standard level graphite calorimeter, developed over a number of years at the NPL, is sufficiently robust to be useable in the proton beams of clinical facilities both in the UK and overseas. The new code of practice involves performing reference dosimetry measurements directly traceable to the primary standard level graphite calorimeter in a clinical proton beam. Calibration of an ionisation chamber is performed in the centre of a standard test volume (STV) of dose, defined here to be a 10 × 10 × 10 cm volume in water, centred at a depth of 15 cm. Further STVs at reduced and increased depths are also utilised. The designated ionisation chambers are Roos-type plane-parallel chambers. This article provides all the necessary background material, formalism, and specifications of reference conditions required to implement reference dosimetry according to this new code of practice. The Annexes provide a detailed review of ion recombination and how this should be assessed (Annex A1) and detailed work instructions for creating and delivering the STVs (Annex A2).

Objective.Fast computation of daily reoptimization is key for an efficient online adaptive proton therapy workflow. Various approaches aim to expedite this process, often compromising daily dose. This study compares Massachusetts General Hospital's (MGH's) online dose reoptimization approach, Paul Scherrer Institute's (PSI's) online replanning workflow and a full reoptimization adaptive workflow for head and neck cancer (H&N) patients.Approach.Ten H&N patients (PSI:5, MGH:5) with daily cone beam computed tomographys (CBCTs) were included. Synthetic CTs were created by deforming the planning CT to each CBCT. Targets and organs at risk (OARs) were deformed on daily images. Three adaptive approaches were investigated: (i) an online dose reoptimization approach modifying the fluence of a subset of beamlets, (ii) full reoptimization adaptive workflow modifying the fluence of all beamlets, and (iii) a full online replanning approach, allowing the optimizer to modify both fluence and position of all beamlets. Two non-adapted (NA) scenarios were simulated by recalculating the original plan on the daily image using: Monte Carlo for NA_MGHand raycasting algorithm for NA_PSI.Main results.All adaptive scenarios from both institutions achieved the prescribed daily target dose, with further improvements from online replanning. For all patients, low-dose CTV D_98%shows mean daily deviations of -2.2%, -1.1%, and 0.4% for workflows (i), (ii), and (iii), respectively. For the online adaptive scenarios, plan optimization averages 2.2 min for (iii) and 2.4 for (i) while the full dose reoptimization requires 72 min. The OA_MGH20%dose reoptimization approach produced results comparable to online replanning for most patients and fractions. However, for one patient, differences up to 11% in low-dose CTV D_98%occurred.Significance.Despite significant anatomical changes, all three adaptive approaches ensure target coverage without compromising OAR sparing. Our data suggests 20% dose reoptimization suffices, for most cases, yielding comparable results to online replanning with a marginal time increase due to Monte Carlo. For optimal daily adaptation, a rapid online replanning is preferable.

Objective. X-ray computed tomography employing low-dose x-ray source is actively researched to reduce radiation exposure. However, the reduced photon count in low-dose x-ray sources leads to severe noise artifacts in analytic reconstruction methods like filtered backprojection. Recently, deep learning (DL)-based approaches employing uni-domain networks, either in the image-domain or projection-domain, have demonstrated remarkable effectiveness in reducing image noise and Poisson noise caused by low-dose x-ray source. Furthermore, dual-domain networks that integrate image-domain and projection-domain networks are being developed to surpass the performance of uni-domain networks. Despite this advancement, dual-domain networks require twice the computational resources of uni-domain networks, even though their underlying network architectures are not substantially different.Approach. The U-Net architecture, a type of Hourglass network, comprises encoder and decoder modules. The encoder extracts meaningful representations from the input data, while the decoder uses these representations to reconstruct the target data. In dual-domain networks, however, encoders and decoders are redundantly utilized due to the sequential use of two networks, leading to increased computational demands. To address this issue, this study proposes a cross-domain DL approach that leverages analytical domain transfer functions. These functions enable the transfer of features extracted by an encoder trained in input domain to target domain, thereby reducing redundant computations. The target data is then reconstructed using a decoder trained in the corresponding domain, optimizing resource efficiency without compromising performance.Main results. The proposed cross-domain network, comprising a projection-domain encoder and an image-domain decoder, demonstrated effective performance by leveraging the domain transfer function, achieving comparable results with only half the trainable parameters of dual-domain networks. Moreover, the proposed method outperformed conventional iterative reconstruction techniques and existing DL approaches in reconstruction quality.Significance. The proposed network leverages the transfer function to bypass redundant encoder and decoder modules, enabling direct connections between different domains. This approach not only surpasses the performance of dual-domain networks but also significantly reduces the number of required parameters. By facilitating the transfer of primal representations across domains, the method achieves synergistic effects, delivering high quality reconstruction images with reduced radiation doses.

Objective.Computational uncertainty and variability of power absorption and temperature rise in humans for radiofrequency (RF) exposure is a critical factor in ensuring human protection. This aspect has been emphasized as a priority. However, accurately modeling head tissue composition and assigning tissue dielectric and thermal properties remains a challenging task. This study investigated the impact of segmentation-based versus segmentation-free models for assessing localized RF exposure.Approach.Two computational head models were compared: one employing traditional tissue segmentation and the other leveraging deep learning to estimate tissue dielectric and thermal properties directly from magnetic resonance images. The finite-difference time-domain method and the bioheat transfer equation was solved to assess temperature rise for local exposure. Inter-subject variability and dosimetric uncertainties were analyzed across multiple frequencies.Main results.The comparison between the two methods for head modeling demonstrated strong consistency, with differences in peak temperature rise of 7.6 ± 6.4%. The segmentation-free model showed reduced inter-subject variability, particularly at higher frequencies where superficial heating dominates. The maximum relative standard deviation in the inter-subject variability of heating factor was 15.0% at 3 GHz and decreased with increasing frequencies.Significance.This study highlights the advantages of segmentation-free deep-learning models for RF dosimetry, particularly in reducing inter-subject variability and improving computational efficiency. While the differences between the two models are relatively small compared to overall dosimetric uncertainty, segmentation-free models offer a promising approach for refining individual-specific exposure assessments. These findings contribute to improving the accuracy and consistency of human protection guidelines against RF electromagnetic field exposure.

Objective: Cone beam computed tomography (CBCT) has become an essential tool in head and neck cancer (HNC) radiotherapy (RT) treatment delivery. Automatic segmentation of the organs at risk (OARs) on CBCT can trigger and accelerate treatment replanning but is still a challenge due to the poor soft tissue contrast, artifacts, and limited field-of-view of these images, alongside the lack of large, annotated datasets to train deep learning models. This study aims to develop a comprehensive framework to segment 25 HN OARs on CBCT to facilitate treatment replanning. Approach. The proposed framework was developed in three steps: (i) refining an in-house framework to segment 25 OARs on computed tomography (CT); (ii) training a deep learning model to segment the same OARs on synthetic CT (sCT) images derived from CBCT using contours propagated from CT as ground truth, integrating high-contrast information from CT and texture features of sCT; and (iii) validating the clinical relevance of sCT segmentations through a dosimetric analysis on an external cohort. Main results. Most OARs achieved a Dice Score Coefficient over 70%, with mean Average Surface Distances of 1.30 mm for CT and 1.27 mm for sCT. The dosimetric analysis demonstrated a strong agreement in the mean dose and D2 (%) values, with most OARs showing non-significant differences between automatic CT and sCT segmentations. Significance. These results support the feasibility and clinical relevance of using deep learning models for OAR segmentation on both CT and CBCT for HNC RT. .