• Progress in Medical Physics
• /
• v.31 no.3
• /
• pp.54-62
• /
• 2020

Dose calculation algorithms play an important role in radiation therapy and are even the basis for optimizing treatment plans, an important feature in the development of complex treatment technologies such as intensity-modulated radiation therapy. We reviewed the past and current status of dose calculation algorithms used in the treatment planning system for radiation therapy. The radiation-calculating dose calculation algorithm can be broadly classified into three main groups based on the mechanisms used: (1) factor-based, (2) model-based, and (3) principle-based. Factor-based algorithms are a type of empirical dose calculation that interpolates or extrapolates the dose in some basic measurements. Model-based algorithms, represented by the pencil beam convolution, analytical anisotropic, and collapse cone convolution algorithms, use a simplified physical process by using a convolution equation that convolutes the primary photon energy fluence with a kernel. Model-based algorithms allowing side scattering when beams are transmitted to the heterogeneous media provide more precise dose calculation results than correction-based algorithms. Principle-based algorithms, represented by Monte Carlo dose calculations, simulate all real physical processes involving beam particles during transportation; therefore, dose calculations are accurate but time consuming. For approximately 70 years, through the development of dose calculation algorithms and computing technology, the accuracy of dose calculation seems close to our clinical needs. Next-generation dose calculation algorithms are expected to include biologically equivalent doses or biologically effective doses, and doctors expect to be able to use them to improve the quality of treatment in the near future.

• Journal of the Korean Society of Radiology
• /
• v.17 no.5
• /
• pp.633-640
• /
• 2023

With the recent development of static and dynamic modulated brachytherapy methods in brachytherapy, which use radiation shielding to modulate the dose distribution to deliver the dose, the amount of parameters and data required for dose calculation in inverse treatment planning and treatment plan optimization algorithms suitable for new directional beam intensity modulated brachytherapy is increasing. Although intensity-modulated brachytherapy enables accurate dose delivery of radiation, the increased amount of parameters and data increases the elapsed time required for dose calculation. In this study, a GPU-based CUDA-accelerated dose calculation algorithm was constructed to reduce the increase in dose calculation elapsed time. The acceleration of the calculation process was achieved by parallelizing the calculation of the system matrix of the volume of interest and the dose calculation. The developed algorithms were all performed in the same computing environment with an Intel (3.7 GHz, 6-core) CPU and a single NVIDIA GTX 1080ti graphics card, and the dose calculation time was evaluated by measuring only the dose calculation time, excluding the additional time required for loading data from disk and preprocessing operations. The results showed that the accelerated algorithm reduced the dose calculation time by about 30 times compared to the CPU-only calculation. The accelerated dose calculation algorithm can be expected to speed up treatment planning when new treatment plans need to be created to account for daily variations in applicator movement, such as in adaptive radiotherapy, or when dose calculation needs to account for changing parameters, such as in dynamically modulated brachytherapy.

• Journal of radiological science and technology
• /
• v.43 no.6
• /
• pp.431-441
• /
• 2020

The purpose of this study was to construct a model of MVCT(Megavoltage Computed Tomography) dose calculation by using Dosimetry Check™, a program that radiation treatment dose verification, and establish a protocol that can be accumulated to the radiation treatment dose distribution. We acquired sinogram of MVCT after air scan in Fine, Normal, Coarse mode. Dosimetry Check™(DC) program can analyze only DICOM(Digital Imaging Communications in Medicine) format, however acquired sinogram is dat format. Thus, we made MVCT RC-DICOM format by using acquired sinogram. In addition, we made MVCT RP-DICOM by using principle of generating MLC(Multi-leaf Collimator) control points at half location of pitch in treatment RP-DICOM. The MVCT imaging dose in fine mode was measured by using ionization chamber, and normalized to the MVCT dose calculation model, the MVCT imaging dose of Normal, Coarse mode was calculated by using DC program. As a results, 2.08 cGy was measured by using ionization chamber in Fine mode and normalized based on the measured dose in DC program. After normalization, the result of MVCT dose calculation in Normal, Coarse mode, each mode was calculated 0.957, 0.621 cGy. Finally, the dose resulting from the process for acquisition of MVCT can be accumulated to the treatment dose distribution for dose evaluation. It is believed that this could be contribute clinically to a more realistic dose evaluation. From now on, it is considered that it will be able to provide more accurate and realistic dose information in radiation therapy planning evaluation by using Tomotherapy.

• The Journal of Korean Society for Radiation Therapy
• /
• v.24 no.2
• /
• pp.137-147
• /
• 2012

Purpose: Analyze the Effectiveness of Radiation Treatment Planning by dose calculation and optimization algorithm, apply consideration of actual treatment planning, and then suggest the best way to treatment planning protocol. Materials and Methods: The treatment planning system use Eclipse 10.0. (Varian, USA). PBC (Pencil Beam Convolution) and AAA (Anisotropic Analytical Algorithm) Apply to Dose calculation, DVO (Dose Volume Optimizer 10.0.28) used for optimized algorithm of Intensity Modulated Radiation Therapy (IMRT), PRO II (Progressive Resolution Optimizer V 8.9.17) and PRO III (Progressive Resolution Optimizer V 10.0.28) used for optimized algorithm of VAMT. A phantom for experiment virtually created at treatment planning system, $30{\times}30{\times}30$ cm sized, homogeneous density (HU: 0) and heterogeneous density that inserted air assumed material (HU: -1,000). Apply to clinical treatment planning on the basis of general treatment planning feature analyzed with Phantom planning. Results: In homogeneous density phantom, PBC and AAA show 65.2% PDD (6 MV, 10 cm) both, In heterogeneous density phantom, also show similar PDD value before meet with low density material, but they show different dose curve in air territory, PDD 10 cm showed 75%, 73% each after penetrate phantom. 3D treatment plan in same MU, AAA treatment planning shows low dose at Lung included area. 2D POP treatment plan with 15 MV of cervical vertebral region include trachea and lung area, Conformity Index (ICRU 62) is 0.95 in PBC calculation and 0.93 in AAA. DVO DVH and Dose calculation DVH are showed equal value in IMRT treatment plan. But AAA calculation shows lack of dose compared with DVO result which is satisfactory condition. Optimizing VMAT treatment plans using PRO II obtained results were satisfactory, but lower density area showed lack of dose in dose calculations. PRO III, but optimizing the dose calculation results were similar with optimized the same conditions once more. Conclusion: In this study, do not judge the rightness of the dose calculation algorithm. However, analyzing the characteristics of the dose distribution represented by each algorithm, especially, a method for the optimal treatment plan can be presented when make a treatment plan. by considering optimized algorithm factors of the IMRT or VMAT that needs to optimization make a treatment plan.

• Radiation Oncology Journal
• /
• v.2 no.1
• /
• pp.129-137
• /
• 1984

In brachytherapy, it is important to determine the positions of the radiation sources which are inserted into a patient and to estimate the dose resulting from the treatment. Calculation of the dose distribution throughout an implant is so laborious that it is rarely done by manual methods except for model cases. It is possible to calculate isodose distributions and tumor doses for individual patients by the use of a microcomputer. In this program, the dose rate and dose distributions are calculated by numerical integration of point source and the localization of radiation sources are obtained from two radiographs at right angles taken by a simulator developed for the treatment planning. By using microcomputer for brachytherapy, we obtained the result as following 1. Dose calculation and irradiation time for tumor could be calculated under one or five seconds after input data. 2. It was same value under$\pm2\%$ error between dose calculation by computer program and measurement dose. 3. It took about five minutes to reconstruct completely dose distribution for intracavitary irradiation. 4. Calculating by computer made remarkly reduction of dose errors compared with Quimby's calculation in interstitial radiation implantation. 5. It could calculate the biological isoffect dose for high and low dose rate activities.

• Journal of radiological science and technology
• /
• v.43 no.4
• /
• pp.265-272
• /
• 2020

The purpose of this study was to investigate the dose-volume indices and radiobiological indices according to the change in dose calculation grid size during the planning of nasopharyngeal cancer VMAT treatment. After performing the VMAT treatment plan using the 3.0 mm dose calculation grid size, dose calculation from 1.0 mm to 5.0 mm was performed repeatedly to obtain a dose volume histogram. The dose volume index and radiobiological index were evaluated using the obtained dose volume histogram. The smaller the dose calculation grid size, the smaller the mean dose for CTV and the larger the mean dose for PTV. For OAR of spinal cord, brain stem, lens and parotid gland, the mean dose did not show a significant difference according to the change in dose calculation grid size. The smaller the grid size, the higher the conformity of the dose distribution as the CI of the PTV increases. The CI and HI showed the best results at 3.0 mm. The smaller the dose calculation grid size, the higher the TCP of the PTV. The smaller the dose calculation grid size, the lower the NTCP of lens and parotid. As a result, when performing the nasopharynx cancer VMAT plan, it was found that the dose calculation grid size should be determined in consideration of dose volume index, radiobiological index, and dose calculation time. According to the results of various experiments, it was determined that it is desirable to apply a grid size of 2.0 - 3.0 mm.

• Progress in Medical Physics
• /
• v.10 no.2
• /
• pp.103-114
• /
• 1999

In this study, as a preliminary study for developing a full 3D photon dose calculation algorithm, We developed 2.5D photon dose calculation algorithm by extending 2D calculation algorithm to allow non-coplanar configurations of photon beams. For this purpose, we defined the 3d patient coordinate system and the 3d beam coordinate system, which are appropriate to 3d treatment planning and dose calculation. and then, calculate a transformation matrix between them. For dose calculation, we extended 2d "Clarkson-Cunningham" model to 3d one, which can calculate wedge fields as well as regular and irregular fields on arbitrary plane. The simple Batho's power-law method was implemented as an inhomogeneity correction. We evaluated the accuracy of our dose model following procedures of AAPM TG#23; radiation treatment planning dosimetry verifications for 4MV of Varian Clinac-4. As results, PDDs (percent depth dose) of cubic fields, the accuracy of calculation are within 1% except buildup region, and $\pm$3% for irregular fields and wedge fields. And for 45$^{\circ}$ oblique incident beam, the deviations between measurements and calculations are within $\pm$4%. In the case of inhomogeneity correction, the calculation underestimate 7% at the lung/water boundary and overestimate 3% at the bone/water boundary. At the conclusions, we found out our model can predict dose with 5% accuracy at the general condition. we expect our model can be used as a tool for educational and research purpose.. purpose..

• Proceedings of the Korean Society of Medical Physics Conference
• /
• 2002.09a
• /
• pp.291-293
• /
• 2002

Accurate dose calculation in radiation treatment planning is most important for successful treatment. Since human body is composed of various materials and not an ideal shape, it is not easy to calculate the accurate effective dose in the patients. Many methods have been proposed to solve the inhomogeneity and surface contour problems. Monte Carlo simulations are regarded as the most accurate method, but it is not appropriate for routine planning because it takes so much time. Pencil beam kernel based convolution/superposition methods were also proposed to correct those effects. Nowadays, many commercial treatment planning systems, including Pinnacle and Helax-TMS, have adopted this algorithm as a dose calculation engine. The purpose of this study is to verify the accuracy of the dose calculated from pencil beam kernel based treatment planning system Helax-TMS comparing to Monte Carlo simulations and measurements especially in inhomogeneous region. Home-made inhomogeneous phantom, Helax-TMS ver. 6.0 and Monte Carlo code BEAMnrc and DOSXYZnrc were used in this study. Dose calculation results from TPS and Monte Carlo simulation were verified by measurements. In homogeneous media, the accuracy was acceptable but in inhomogeneous media, the errors were more significant.

• Progress in Medical Physics
• /
• v.5 no.1
• /
• pp.67-74
• /
• 1994

Dose distribution of point source represents an inverse square law as the distance, Difference of measurement value and calculation value according to moving distance of radiation source show very large error in dose calculation of Brachytherapy. Therefore, in RALS of high dose rate, dose calculation have an important effect in treatment of uterine cervix cancer and recurrent rate. In this paper, authors measured moving distance of radiation source carrying out RALS. And we measured Rectum dose compared with calculationdose.

• Proceedings of the Korean Society of Medical Physics Conference
• /
• 2002.09a
• /
• pp.68-73
• /
• 2002

In standard teletherapy, a treatment plan is generated with the aid of a treatment planning system, but it is common to perform an independent monitor unit verification calculation (MUVC). In exact analogy, we propose and demonstrate that a simple and accurate MUVC in Intensity Modulated Radiotherapy (IMRT) is possible. We introduce a concept of Modified Clarkson Integration (MCI). In MCI, we exploit the rotational symmetry of scattering to simplify the dose calculation. For dose calculation along a central axis (CAX), we first replace the incident IMRT fluence by an azimuthally averaged fluence. Second, the Clarkson Integration is carried over annular sectors instead of over pie sectors. We wrote a computer code, implementing the MCI technique, in order to perform a MUVC for IMRT purposes. We applied the code to IMRT plans generated by CORVUS. The input to the code consists of CORVUS plan data (e.g., DMLC files, jaw settings, MU for each IMRT field, depth to isocenter for each IMRT field), and the output is dose contribution by individual IMRT field to the isocenter. The code uses measured beam data for Sc, Sp, TPR, (D/Mu)$\_$ref/ and includes effects from MLC transmission, and radiation field offset. On a 266 MHZ desktop computer, the code takes less than 15 sec to calculate a dose. The doses calculated with MCI algorithm agreed within +/- 3% with the doses calculated by CORVUS, which uses a 1cm x 1cm pencil beam in dose calculation. In the present version of MCI, skin contour variations and inhomogeneities were neglected.