Test procedures for verification of an electron pencil beam algorithm implemented for treatment planning

The calculation of an electron dose distribution in a patient is a difficult problem because of the presence of tissue and surface inhomogeneities. Verification of the dose planning system is therefore essential. In this investigation, a novel method is used to evaluate a commercially available system (Helax-TMS), at electron energies between 10 and 50 MeV, both for a conventional treatment unit and an MLC-collimated scanned beam unit with a helium-filled treatment head. First, the experiments were designed to verify the local beam database and some fundamental characteristics of the electron beam calculations. Secondly, a number of generalised situations that would be encountered in the clinical treatment planning were evaluated oblique incidence, field shaping with multi-leaf collimator, bolus edges, and air cavities. Dose distributions in two generalised anatomical phantoms simulating a neck and a nose were also analysed. The results have, when so possible, been presented as the dose ratio within the 'flattened area' for dose profiles and down to the 'treatment depth' (80% dose level) for depth doses. In the penumbra region and in the dose fall-off region, the comparison has been represented by the distance deviation between calculated and measured dose profiles or depth doses. A new tool, 'volume integration', was used to evaluate the deviations from a more clinical point of view. Most results were within +/- 2% in dose for volumes larger than a sphere with a diameter of 15 mm, or +/- 2 mm in position. Dose deviations were generally found for oblique incidences and below heterogeneities such as small air cavities and bolus edges in limited volumes.     

The simple excision

A method for the fully computerized determination and optimization of positions of target points and collimator sizes in convergent beam irradiation is presented. In conventional interactive trial and error methods, which are very time consuming, the treatment parameters are chosen according to the operator's experience and improved successively. This time is reduced significantly by the use of a computerized procedure. After the definition of target volume and organs at risk in the CT or MR scans, an initial configuration is created automatically. In the next step the target point positions and collimator diameters are optimized by the program. The aim of the optimization is to find a configuration for which a prescribed dose at the target surface is approximated as close as possible. At the same time dose peaks inside the target volume are minimized and organs at risk and tissue surrounding the target are spared. To enhance the speed of the optimization a fast method for approximate dose calculation in convergent beam irradiation is used. A possible application of the method for calculating the leaf positions when irradiating with a micromultileaf collimator is briefly discussed. The success of the procedure has been demonstrated for several clinical cases with up to six target points.     

铁粉生产计划的佩特里网优化方法

讨论了有序加工类生产过程优化计划的自动生成,并选用模拟退火算法解决该大规模组合优化问题。为减少算法中目标函数的计算量,不仅整个生产过程被分解成若干个子过程,过程模型也被简化。为了补偿模型简化导致的精度下降,引入佩特里网验证上述算法所得到的优化计划,实例证明该方法是有效的     

Fast iterative algorithms for three-dimensional inverse treatment planning

Three types of iterative algorithms, algebraic inverse treatment planning (AITP), simultaneous iterative inverse treatment planning (SIITP), and iterative least-square inverse treatment planning (ILSITP), differentiated according to their updating sequences, were generalized to three dimension with true beam geometry and dose model. A rapid ray-tracing approach was developed to optimize the primary beam components. Instead of recalculating the dose matrix at each iteration, the dose distribution was generated by scaling up or down the dose matrix elements of the previous iteration. This significantly increased the calculation speed. The iterative algorithms started with an initial intensity profile for each beam, specified by a two-dimensional pixel beam map of M elements. The calculation volume was divided into N voxels, and the calculation was done by repeatedly comparing the calculated and desired doses and adjusting the values of the beam map elements to minimize an objective function. In AITP, the iteration is performed voxel by voxel. For each voxel, the dose discrepancy was evaluated and the contributing pencil beams were updated. In ILSITP and SIITP, the iteration proceeded pencil beam by pencil beam instead of voxel by voxel. In all cases, the iteration procedure was repeated until the best possible dose distribution was achieved. The algorithms were applied to two examples and the results showed that the iterative techniques were able to produce superior isodose distributions.     

Monte Carlo and analytical calculation of proton pencil beams for computerized treatment plan optimization

Proton pencil beams in water, in a format suitable for treatment planning algorithms and covering the radiotherapy energy range (50-250 MeV), have been calculated using a modified version of the Monte Carlo code PTRAN. A simple analytical model has also been developed for calculating proton broad-beam dose distributions which is in excellent agreement with the Monte Carlo calculations. Radial dose distributions are also calculated analytically and narrow proton pencil-beam dose distributions derived. The physical approximations in the Monte Carlo code and in the analytical model together with their limitations are discussed. Examples showing the use of the calculated set of proton pencil beams as input to an existing photon treatment planning algorithm based on biological optimization are given for fully 3D scanned proton pencil beams; these include intensity modulated beams with range shift and scanning in the transversal plane.     

Description and evaluation of a new 3-D computerized treatment planning dose compensator system

Evaluation has been performed of compensators generated by means of a computerized three-dimensional treatment planning system that can utilize either digitized slice profiles or CT scans. Two methods of calculating compensator thickness are used: the modified Batho power law (dSAR) method for digitized profiles and the equivalent TAR (eqTAR) method for CT scans. This system not only compensates for patient surface contours but also compensates for internal inhomogeneities. In addition, any required wedging will be incorporated in the compensator generation. This system has been tested for a number of extreme cases with inhomogeneities and sloping contours. Good agreement was obtained between the measured and computer calculated dose profiles especially along the central axis of the beam. A "Profile Uniformity Index" was defined to quantify the goodness of dose compensation in three dimensions. Compensation using this system can achieve good dose uniformity within the target volume in all clinical cases and is definitely an improvement over systems based solely on tissue deficit.     

Optimal Design with Probabilistic Objective and Constraints

Significant challenges are associated with solving optimal structural design problems involving the failure probability in the objective and constraint functions. In this paper, we develop gradient-based optimization algorithms for estimating the solution of three classes of such problems in the case of continuous design variables. Our approach is based on a sequence of approximating design problems, which is constructed and then solved by a semiinfinite optimization algorithm. The construction consists of two steps: First, the failure probability terms in the objective function are replaced by auxiliary variables resulting in a simplified objective function. The auxiliary variables are determined automatically by the optimization algorithm. Second, the failure probability constraints are replaced by a parametrized first-order approximation. The parameter values are determined in an adaptive manner based on separate estimations of the failure probability. Any computational reliability method, including first-order reliability and second-order reliability methods and Monte Carlo simulation, can be used for this purpose. After repeatedly solving the approximating problem, an approximate solution of the original design problem is found, which satisfies the failure probability constraints at a precision level corresponding to the selected reliability method. The approach is illustrated by a series of examples involving optimal design and maintenance planning of a reinforced concrete bridge girder.     

A new genetic algorithm technique in optimization of permanent 125I prostate implants

Real time optimized treatment planning at the time of the implant is desirable for ultrasound-guided transperineal 125I permanent prostate implants. Currently available optimization algorithms are too slow to be used in the operating room. The goal of this work is to develop a robust optimization algorithm, which is suitable for such application. Three different genetic algorithms (sGA, sureGA and securGA) were developed and compared in terms of the number of function evaluations and the corresponding fitness. The optimized dose distribution was achieved by searching the best seed distribution through the minimization of a cost function. The cost function included constraints on the periphery dose of the planned target volume, the dose uniformity within the target volume, and the dose to the critical structure. Adjustment between the peripheral dose, the dose uniformity and critical structure dose can be achieved by varying the weighting factors in the cost function. All plans were evaluated in terms of the dose nonuniformity ratio, the conformation number and the dose volume histograms. Among these three GA algorithms, the securGA provided the best performance. Within 2500 function evaluations, the near optimum results were obtained. For a large target volume (5 cm x 4 cm x 4.5 cm) including urethra with 20 needles, the computer time needed for the optimization was less than 5 min on a HP735 workstation. The results showed that once the best set of parameters was found, they were applicable for all sizes of prostate volume. For a fixed needle geometry, the optimized plan showed much better dose distribution than that of nonoptimized plan. If the critical structure was considered in the optimization, the dose to the critical structure could be minimized. In the cases of irregular and skewed needle geometry, the optimized treatment plans were almost as good as ideal needle geometry. It is concluded that this new genetic algorithm (securGA) allows for an efficient and rapid optimization of dose distribution, which is suitable for real time treatment planning optimization for ultrasound-guided prostate implant.     

步进炉炉底机械运动速度的合理设定与电液控制

轧管车间步进式再加热炉炉底机械运动速度的突变不仅直接影响到步进炉炉底机械设备和炉底水冷系统的使用寿命,还会对炉内加热的荒管造成损伤,影响成品质量。所以,合理设计炉底机械液压系统、正确设定并计算炉底机械运动速度曲线、优化各种液压元件的配置,对减少炉底机械运动速度突变是十分重要的。     

An on-line computer graphics terminal for radiotherapy treatment planning

A programmed graphics terminal has been connected on-line to a large time-shared computer for calculating dose distributions in radiotherapy treatment planning and provides a viable alternative to dedicated systems and batch working. The terminal equipment is based on a mini-computer and includes a function-key devise for outline input, a large-screen refresh oscilloscope for viewing results and an X-Y plotter for hard-copy. Radiotherapy dose computation programs in standard Fortran are stored and run on a large remote computer with graphical interaction at the terminal. External beam programs can calculate dose distributions for most commonly used treatment situations and can compute in off-axis planes. Data input is fully interactive and easy to understnad. Dose distributions are displayed as isodose contours. Advantages of the system include accuracy, speed, ease of use and maintenance, and transferability of the programs between different host computers.     

Intensity modulation with electrons: calculations, measurements and clinical applications

Intensity modulation of electron beams is one step towards truly conformal therapy. This can be realized with the MM50 racetrack microtron that utilizes a scanning beam technique. By adjusting the scan pattern it is possible to obtain arbitrary fluence distributions. Since the monitor chambers in the treatment head are segmented in both x- and y-directions it is possible to verify the fluence distribution to the patient at any time during the treatment. Intensity modulated electron beams have been measured with film and a plane parallel chamber and compared with calculations. The calculations were based on a pencil beam method. An intensity distribution at the multileaf collimator (MLC) level was calculated by superposition of measured pencil beams over scan patterns. By convolving this distribution with a Gaussian pencil beam, which has propagated from the MLC to the isocentre, a fluence distribution at isocentre level was obtained. The agreement between calculations and measurements was within 2% in dose or 1 mm in distance in the penumbra zones. A standard set of intensity modulated electron beams has been developed. These beams have been implemented in a treatment planning system and are used for manual optimization. A clinical example (prostate) of such an application is presented and compared with a standard irradiation technique.     

A three-dimensional algorithm for optimizing beam weights and wedge filters

An essential step towards optimizing and automating radiation therapy treatment planning is to develop an effective algorithm to find the optimal beam weights and wedge filters for a given set of beam directions and modalities. This problem is solved by introducing a variable transformation based on the universal and omni wedge principles. Instead of directly optimizing an objective function with respect to wedge angles and orientations, each field is first decomposed into a superposition of an open field and two orthogonal wedged fields. This transforms the problem of finding J beam weights, wedge angles, and orientations to that of optimizing a system with 3J beam weights (J open beams and 2J nominal wedged beams), where J is the total number of incident beam directions. An iterative algorithm based on a method originally developed for image reconstruction is used to find the 3J beam weights. The technique is applied to a few clinical cases. Treatment plans are improved compared to those obtained through the conventional manual trial and error planning process. In addition, planning time and effort are greatly reduced.     

Dose accuracy check of the 3D electron beam algorithm in a treatment planning system

The accuracy of the recently implemented three-dimensional electron beam dose calculating algorithm in CADPLAN version 2.62 manufactured by Varian Dosetek was investigated. The algorithm uses a generalized Gaussian pencil beam model and the dose distributions are calculated as the sum of three weighted Gaussians. To use the calculating program in an optimum way, one needs to know the dose calculation accuracy of the algorithm as well as its limitations. This investigation includes comparisons of measured relative dose distributions with calculated dose distributions and also comparisons of measured and calculated monitor units. The geometries tested were quadratic fields, irregularly shaped fields, oblique fields, irregularly shaped phantom surfaces and internal heterogeneities and were most often irradiated with 8 and 20 MeV electrons. The results indicate that the algorithm is well suited for clinical three-dimensional dose planning. Some deviations occurred but they were most often within the limits of international criteria of acceptability.     

混沌人工鱼群的鲁棒保性能控制权值矩阵优化方法

针对鲁棒保性能控制中的权值矩阵依赖经验选取,无法最大限度的减小系统保守性的问题,提出了一种基于混沌人工鱼群算法的鲁棒保性能控制权值矩阵优化方法.该方法中,将保性能控制鲁棒界作为优化的目标函数来寻找最优权值矩阵是整个算法实现的关键.该种改进的人工鱼群优化算法融合了混沌搜索与自适应步长和视野的人工鱼群优化算法,有效的解决了基本人工鱼群算法的后期收敛速度慢、易陷入局部最优等缺点.通过测试函数对比验证了该种改进人工鱼群优化算法的优越性,并通过应用实例验证了该权值矩阵优化方法的有效性.      

Paleolithic and Neolithic lineages in the European mitochondrial gene pool

PURPOSE: This article presents a new optimization method for stereotactic radiosurgery treatment planning for gamma unit treatment system. METHODS AND MATERIALS: The gamma unit has been utilized in stereotactic radiosurgery for about 30 years, but the usual procedure for a physician-physicist team to design a treatment plan is a trial-and-error approach. Isodose curves are viewed on two-dimensional computed tomography (CT) or magnetic resonance (MR) image planes, which is not only time consuming but also seldom achieves the optimal treatment plan, especially when the isocenter weights are regarded. We developed a treatment-planning system on a computer workstation in which Powell's optimization method is realized. The optimization process starts with the initial parameters (the number of isocenters as well as corresponding 3D isocenters' coordinates, collimator sizes, and weight factors) roughly determined by the physician-physicist team. The objective function can be changed to consider protection of sensitive tissues. RESULTS: We use the plan parameters given by a well-trained physician-physicist team, or ones that the author give roughly as the initial parameters for the optimization procedure. Dosimetric results of optimization show a better high dose-volume conformation to the target volume compared to the doctor's plan. CONCLUSION: This method converges quickly and is not sensitive to the initial parameters. It achieves an excellent conformation of the estimated isodose curves with the contours of the target volume. If the initial parameters are varied, there will be a little difference in parameters' configuration, but the dosimetric results proved almost to be the same.     

Adaptive Parametric Identification Scheme for a Class of Nondeteriorating and Deteriorating Nonlinear Hysteretic Behavior

The adaptive parametric identification of deteriorating and nondeteriorating nonlinear hysteretic phenomena is considered using a generalization of Masing model based on the observed memory behavior of distributed element models. The model permits a parametric identification to be performed using nonlinear optimization techniques for arbitrary response time histories. A changing objective function, defined as the normalized force estimation error over a shifting window of recent data, is employed so that classic nonlinear optimization techniques can be used for the adaptive identification problem. A variation of the steepest descent method is used with significant modifications. To achieve the best performance for any given problem, a set of a priori numeric tests are suggested to design the identification scheme. The design identification scheme exhibits a very good performance in identifying the correct values of the parameters and is rather robust in dealing with noise. The proposed approach has applications to adaptive identification of much wider types of nonlinear rate-dependent hysteretic behavior. Also, the set of guidelines proposed by the authors is a contribution toward having more effective autonomous identification schemes, using minimal information about the model and input.     

A patient-specific Monte Carlo dose-calculation method for photon beams

A patient-specific, CT-based, Monte Carlo dose-calculation method for photon beams has been developed to correctly account for inhomogeneity in the patient. The method employs the EGS4 system to sample the interaction of radiation in the medium. CT images are used to describe the patient geometry and to determine the density and atomic number in each voxel. The user code (MCPAT) provides the data describing the incident beams, and performs geometry checking and energy scoring in patient CT images. Several variance reduction techniques have been implemented to improve the computation efficiency. The method was verified with measured data and other calculations, both in homogeneous and inhomogeneous media. The method was also applied to a lung treatment, where significant differences in dose distributions, especially in the low-density region, were observed when compared with the results using an equivalent pathlength method. Comparison of the DVHs showed that the Monte Carlo calculated plan predicted an underdose of nearly 20% to the target, while the maximum doses to the cord and the heart were increased by 25% and 33%, respectively. These results suggested that the Monte Carlo method may have an impact on treatment designs, and also that it can be used as a benchmark to assess the accuracy of other dose calculation algorithms. The computation time for the lung case employing five 15-MV wedged beams, with an approximate field size of 13 X 13 cm and the dose grid size of 0.375 cm, was less than 14 h on a 175-MHz computer with a standard deviation of 1.5% in the high-dose region.     

Simple look-up table of optimized dwell time intervals using a high-dose-rate remote afterloader for endovascular irradiation

PURPOSE: This article's objective is to develop a simple methodology deliver a uniform radiation dose to the wall of a narrow peripheral artery for preventing restenosis using a high-dose-rate (HDR) 192Ir remote afterloader. METHODS AND MATERIALS: Based upon published two-dimensional data such as anisotropy factors of an HDR 192Ir source calculated from the Monte-Carlo method, arterial wall doses at a close range from an HDR source may be easily calculated using the special formula suggested in Task Group Report No. 43 published by the American Association of Physicists in Medicine. An optimization procedure was used to calculate the optimized dwell times for delivering a uniform dose along arterial walls for various arterial diameters and lengths of lesions. RESULTS: Based on lengths of the stenosis and diameters of arteries or angioplasty balloons, a set of simple look-up tables for optimal dwell time intervals of endovascular radiation treatment have been developed for the MicroSelectron HDR remote afterloader. CONCLUSION: Doses for endovascular irradiation have been accurately calculated with anisotropy factors. For delivering uniform doses along the arterial wall, a set of look-up tables listed for optimal dwell times is available for the HDR remote afterloader.     

A dual source photon beam model used in convolution/superposition dose calculations for clinical megavoltage x-ray beams

A realistic photon beam model based on Monte Carlo simulation of clinical linear accelerators was implemented in a convolution/superposition dose calculation algorithm. A primary and an extra-focal sources were used in this beam model to represent the direct photons from the target and the scattered photons from other head structures, respectively. The effect of the finite size of the extra-focal source was modeled by a convolution of the source fluence distribution with the collimator aperture function. Relative photon output in air (Sc) and in phantom (Scp) were computed using the convolution method with this new photon beam model. Our results showed that in a 10 MV photon beam, the Sc, Sp (phantom scatter factor), and Scp factors increased by 11%, 10%, and 22%, respectively, as the field size changed from 3 x 3 cm2 to 40 x 40 cm2. The variation of the Sc factor was contributed mostly by an increase of the extra-focal radiation with field size. The radiation backscattered into the monitor chamber inside the accelerator head affected the Sc by about 2% in the same field range. The output factors in elongated fields, asymmetric fields, and blocked fields were also investigated in this study. Our results showed that if the effect of the backscattered radiation was taken into account, output factors in these treatment fields can be predicted accurately by our convolution algorithm using the dual source photon beam model.