ترغب بنشر مسار تعليمي؟ اضغط هنا

The fixed-point iteration method for IMRT optimization with truncated dose deposition coefficient matrix

160   0   0.0 ( 0 )
 نشر من قبل Zhen Tian
 تاريخ النشر 2013
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

In the treatment plan optimization for intensity modulated radiation therapy (IMRT), dose-deposition coefficient (DDC) matrix is often pre-computed to parameterize the dose contribution to each voxel in the volume of interest from each beamlet of unit intensity. However, due to the limitation of computer memory and the requirement on computational efficiency, in practice matrix elements of small values are usually truncated, which inevitably compromises the quality of the resulting plan. A fixed-point iteration scheme has been applied in IMRT optimization to solve this problem, which has been reported to be effective and efficient based on the observations of the numerical experiments. In this paper, we aim to point out the mathematics behind this scheme and to answer the following three questions: 1) whether the fixed-point iteration algorithm converges or not? 2) when it converges, whether the fixed point solution is same as the original solution obtained with the complete DDC matrix? 3) if not the same, whether the fixed point solution is more accurate than the naive solution of the truncated problem obtained without the fixed-point iteration? To answer these questions, we first performed mathematical analysis and deductions using a simplified fluence map optimization (FMO) model. Then we conducted numerical experiments on a head-and-neck patient case using both the simplified and the original FMO model. Both our mathematical analysis and numerical experiments demonstrate that with proper DDC matrix truncation, the fixed-point iteration can converge. Even though the converged solution is not the one that we obtain with the complete DDC matrix, the fixed-point iteration scheme could significantly improve the plan accuracy compared with the solution to the truncated problem obtained without the fixed-point iteration.

قيم البحث

اقرأ أيضاً

Purpose: This study aims to optimize and characterize the response of a mPSD for in vivo dosimetry in HDR brachytherapy. Methods: An exhaustive analysis was carried out in order to obtain an optimized mPSD design that maximize the scintillation light collection produced by the interaction of ionizing photons. Several mPSD prototypes were built and tested in order to determine the appropriate order of scintillators relative to the photodetector, as well as their length as a function of the scintillation light emitted. Scintillators BCF-60, BCF-12 and BCF-10 constituted the mPSD sensitive volume.Each scintillator contribution to the total spectrum was determined by irradiations in the low energy range.For the best mPSD design, a numerical optimization was done in order to select the optical components that better match the light emission profile. The optimized dosimetric system was used for HDR brachytherapy dose determination. The system performance was quantified in term of signal to noise ratio and signal to background ratio. Results: It was determined that BCF-60 should be placed at the distal position, BCF-12 in the center and BCF-10 at proximal position with respect to the photodetector.This configuration allowed for optimized light transmission through the collecting fiber, avoiding inter-scintillator excitation and self-absorption effects.The optimized luminescence system allowed for signal deconvolution using a multispectral approach, extracting the dose to each element while taking into account Cerenkov stem effect.Differences between the mPSD measurements and TG-43 remain below 5%. In all measurement conditions, the system was able to properly differentiate the produced scintillation signal from the background one. Conclusions: A mPSD was constructed and optimized for HDR brachytherapy dosimetry, enabling real time dose determination, up to 6.5cm from the 192Ir source.
A new variant of the pencil-beam (PB) algorithm for dose distribution calculation for radiotherapy with protons and heavier ions, the grid-dose spreading (GDS) algorithm, is proposed. The GDS algorithm is intrinsically faster than conventional PB alg orithms due to approximations in convolution integral, where physical calculations are decoupled from simple grid-to-grid energy transfer. It was effortlessly implemented to a carbon-ion radiotherapy treatment planning system to enable realistic beam blurring in the field, which was absent with the broad-beam (BB) algorithm. For a typical prostate treatment, the slowing factor of the GDS algorithm relative to the BB algorithm was 1.4, which is a great improvement over the conventional PB algorithms with a typical slowing factor of several tens. The GDS algorithm is mathematically equivalent to the PB algorithm for horizontal and vertical coplanar beams commonly used in carbon-ion radiotherapy while dose deformation within the size of the pristine spread occurs for angled beams, which was within 3 mm for a single proton pencil beam of $30^circ$ incidence, and needs to be assessed against the clinical requirements and tolerances in practical situations.
The aim of this study was to investigate the impact of decay data provided by the newly developed stochastic atomic relaxation model BrIccEmis on dose point kernels (DPKs - radial dose distribution around a unit point source) and S-values (absorbed d ose per unit cumulated activity) of 14 Auger electron (AE) emitting radionuclides, namely 67Ga, 80mBr, 89Zr, 90Nb, 99mTc, 111In, 117mSn, 119Sb, 123I, 124I, 125I, 135La, 195mPt and 201Tl. Radiation spectra were based on the nuclear decay data from the medical internal radiation dose (MIRD) RADTABS program and the BrIccEmis code, assuming both an isolated-atom and condensed-phase approach. DPKs were simulated with the PENELOPE Monte Carlo (MC) code using event-by-event electron and photon transport. S-values for concentric spherical cells of various sizes were derived from these DPKS using appropriate geometric reduction factors. The number of Auger and Coster-Kronig (CK) electrons and x-ray photons released per nuclear decay (yield) from MIRD-RADTABS were consistently higher than those calculated using BrIccEmis. DPKs for the electron spectra from BrIccEmis were considerably different from MIRD-RADTABS in the first few hundred nanometres from a point source where most of the Auger electrons are stopped. S-values were, however, not significantly impacted as the differences in DPKS in the sub-micrometre dimension were quickly diminished in larger dimensions. Overestimation in the total AE energy output by MIRD-RADTABS leads to higher predicted energy deposition by AE emitting radionuclides, especially in the immediate vicinity of the decaying radionuclides. This should be taken into account when MIRD-RADTABS data are used to simulate biological damage at nanoscale dimensions.
For the unsorted database quantum search with the unknown fraction $lambda$ of target items, there are mainly two kinds of methods, i.e., fixed-point or trail-and-error. (i) In terms of the fixed-point method, Yoder et al. [Phys. Rev. Lett. 113, 2105 01 (2014)] claimed that the quadratic speedup over classical algorithms has been achieved. However, in this paper, we point out that this is not the case, because the query complexity of Yoders algorithm is actually in $O(1/sqrt{lambda_0})$ rather than $O(1/sqrt{lambda})$, where $lambda_0$ is a known lower bound of $lambda$. (ii) In terms of the trail-and-error method, currently the algorithm without randomness has to take more than 1 times queries or iterations than the algorithm with randomly selected parameters. For the above problems, we provide the first hybrid quantum search algorithm based on the fixed-point and trail-and-error methods, where the matched multiphase Grover operations are trialed multiple times and the number of iterations increases exponentially along with the number of trials. The upper bound of expected queries as well as the optimal parameters are derived. Compared with Yoders algorithm, the query complexity of our algorithm indeed achieves the optimal scaling in $lambda$ for quantum search, which reconfirms the practicality of the fixed-point method. In addition, our algorithm also does not contain randomness, and compared with the existing deterministic algorithm, the query complexity can be reduced by about 1/3. Our work provides an new idea for the research on fixed-point and trial-and-error quantum search.
294 - Lin Ma , Mingli Chen , Xuejun Gu 2021
Purpose: To develop a model to generate volumetric dose distribution from two isodose surfaces (iso-surfaces), and to interactively tune dose distribution by iso-surface dragging. Methods: We model volumetric dose distribution as analytical extension of two iso-surfaces with the extension variables as distances to iso-surfaces. We built a 3D lookup table (LUT) which are generated based on clinical dose distributions. Two LUT tables store the mean and standard deviation of voxel dose values of clinical doses and binned as distance to 100% iso-surface, reference iso-surface and reference dose level. The process of interactive tuning starts from a given base plan. A user drags iso-surface for a desired carving. Our method responds with tuned dose. The derivation of tuned dose follows two steps. Dose is extended from the two user-desired iso-surfaces (eg.100% and 50%) to the whole patient volume by table lookup, using distances to two iso-surfaces and reference dose level as keys. Then we fine tune the extended dose by a correction strategy utilizing the information of base plan. Results: We validated this method on coplanar VMAT doses of post-operative prostate plans. The LUT was populated by dose distributions of 27 clinical plans. We optimized two plans with different rectum sparing for an independent case to mimic the process of dose tuning. The plan with less rectum sparing is set as base plan. The 50% iso-surface of the more-sparing plan is defined as the desired iso-surface input. The dose output by our method (expansion and correction) agrees with the more-sparing plan obtained by optimization, in terms of gamma (97.2%), DVH and profiles. The overall dose reconstruction time is within two seconds. Conclusion: We developed a distance-to-isosurface based volumetric dose reconstruction method, and applied it to interactive tuning with iso-surface dragging.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا