• 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..

• 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.

• 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.

• Progress in Medical Physics
• /
• v.10 no.3
• /
• pp.133-140
• /
• 1999

In this paper, as a preliminary study for developing a full 3D electron dose calculation algorithm, We developed 2.5D electron dose calculation algorithm by extending 2D pencil-beam model to consider three dimensional geometry such as air-gap and obliquity appropriately. The dose calculation algorithm was implemented using the IDL5.2(Research Systems Inc., USA), For calculation of the Hogstrom's pencil-beam algorithm, the measured data of the central-axis depth-dose for 12 MeV(Siemens M6740) and the linear stopping power and the linear scattering power of water and air from ICRU report 35 was used. To evaluate the accuracy of the implemented program, we compared the calculated dose distribution with the film measurements in the three situations; the normal incident beam, the 45$^{\circ}$ oblique incident beam, and the beam incident on the pit-shaped phantom. As results, about 120 seconds had been required on the PC (Pentium III 450MHz) to calculate dose distribution of a single beam. It needs some optimizing methods to speed up the dose calculation. For the accuracy of dose calculation, in the case of the normal incident beam of the regular and irregular shaped field, at the rapid dose gradient region of penumbra, the errors were within $\pm$3 mm and the dose profiles were agreed within 5%. However, the discrepancy between the calculation and the measurement were about 10% for the oblique incident beam and the beam incident on the pit-shaped phantom. In conclusions, we expended 2D pencil-beam algorithm to take into account the three dimensional geometry of the patient. And also, as well as the dose calculation of irregular field, the irregular shaped body contour and the air-gap could be considered appropriately in the implemented program. In the near future, the more accurate algorithm will be implemented considering inhomogeneity correction using CT, and at that time, the program can be used as a tool for educational and research purpose. This study was supported by a grant (＃HMP-98-G-1-016) of the HAN(Highly Advanced National) Project, Ministry of Health & Welfare, R.O.K.

• Journal of Radiation Protection and Research
• /
• v.40 no.1
• /
• pp.55-64
• /
• 2015

Dose calculations which are a crucial requirement for radiotherapy treatment planning systems require accuracy and rapid calculations. The conventional radiotherapy treatment planning dose algorithms are rapid but lack precision. Monte Carlo methods are time consuming but the most accurate. The new combined system that Monte Carlo methods calculate part of interesting domain and the rest is calculated by neural can calculate the dose distribution rapidly and accurately. The preliminary study showed that neural networks can map functions which contain discontinuous points and inflection points which the dose distributions in inhomogeneous media also have. Performance results between scaled conjugated gradient algorithm and Levenberg-Marquardt algorithm which are used for training the neural network with a different number of neurons were compared. Finally, the dose distributions of homogeneous phantom calculated by a commercialized treatment planning system were used as training data of the neural network. In the case of homogeneous phantom;the mean squared error of percent depth dose was 0.00214. Further works are programmed to develop the neural network model for 3-dimensinal dose calculations in homogeneous phantoms and inhomogeneous phantoms.

• Proceedings of the Korean Society of Medical Physics Conference
• /
• 2003.09a
• /
• pp.33-33
• /
• 2003

목적 : 현재 국내에서 사용되고 있는 방사선치료계획 장비는 거의 모든 방사선종양학과에서 외국 회사의 제품을 사용하고 있는 실정이다. 본 연구는 국내 기술로 개발된 3 차원 조형 방사선치료계획 장비인 Core Plan의 선량계산 알고리즘의 정확성 및 방사선 치료계획 장비로서의 유용성에 관해 평가하고 본 장비의 특징에 대해 간략히 소개하고자 한다. 재료 및 방법 : 본 연구는 2002 년 11 월에 가톨릭대학교 성모병원 방사선종양학과와 서울씨앤제이의 연구계약에 의해 시행되었다. 본 장비에 대한 평가는 방사선분포 및 계산상의 정확성과 임상적용시의 유용성의 관점에서 시행되었다. 본 장비에서 이용된 광자선 선량계산은 Clarkson-Cunningham 모델이며 전자선 선량계산은 2.5D Hogstrom 알고리즘이다. 방사선분포 및 계산상의 정확성 평가를 위하여 방사선치료 장비는 본 병원이 보유한 Clinac 2100CD (Varian, USA)를 이용하였고 폴리스티렌 팬텀과 필름 및 이온 전리함을 이용하여 방사선분포 및 계산상의 정확성을 평가하였다. 방사선분포의 평가 방법은 필름을 이용한 방사선분포의 중심부단면 선량분포와 CorePlan 에서 재현된 방사선분포의 동일면 선량분포를 비교하였다. 임상적용은 2003 년 3 월부터 7 월까지 방사선치료를 받은 50 명의 환자를 대상으로 분석하였다. 본 시험에 적용된 환자는 본 병원이 보유한 3 차원 방사선치료계획 장비인 Prowess 3D (SSGI Inc., USA)를 이용하여 실제 치료된 환자들이며 이 환자를 대상으로 CorePlan에서 동일하게 재현하여 비교하였다. 결과 및 결론 : 방사선분포 및 계산상의 정확성 평가에서는 실제 측정된 결과와 CorePlan에서 재현된 결과가 모두 $\pm$3% 이내로 평가되었다. 50 명의 환자를 대상으로 시행한 임상시험 결과 Prowess 3D에서 나타낸 결과와 비교하여 두경부에서는 1.678%, 흉부 1.578%, 복부 1.271%의 선량값의 오차를 보였다. 본 연구를 진행하는 과정에서 많은 부분의 프로그램 수정이 있었으며 실제임상에 필요한 부분에 대한 추가 및 보완이 이루어졌다. 앞으로 진행될 과정은 실제임상의 사용에 필요한 부분의 계속적인 업그레이드와 전자선에 대한 방사선분포 및 계산정확성 평가, 임상적용에 있어서 Prowess3D 뿐만 아니라 다양한 방사선 치료계획 장비와의 비교를 할 예정이다.

• Progress in Medical Physics
• /
• v.23 no.2
• /
• pp.106-113
• /
• 2012

The purpose of this study is to evaluate the variation of radiation dose distribution for liver tumor located in liver dome and for the interest organs(normal liver, kidney, stomach) with the pencil beam convolution (PBC) algorithm versus anisotropic Analyticalal algorithm (AAA) of the Varian Eclipse treatment planning system, The target volumes from 20 liver cancer patients were used to create treatment plans. Treatment plans for 10 patients were performed in Stereotactic Body Radiation Therapy (SBRT) plan and others were performed in 3 Dimensional Conformal Radiation Therapy (3DCRT) plan. dose calculation was recalculated by AAA algorithm after dose calculation was performed by PBC algorithm for 20 patients. Plans were optimized to 100% of the PTV by the Prescription Isodose in Dose Calculation with the PBC algorithm. Plans were recalculated with the AAA, retaining identical beam arrangements, monitor units, field weighting and collimator condition. In this study, Total PTV was to be statistically significant (SRS: p=0.018, 3DCRT: p=0.006) between PBC and AAA algorithm. and in the case of PTV, ITV in liver dome, plans for 3DCRT were to be statistically significant respectively (p=0.013, p=0.024). normal liver and kidney were to be statistically significant (p=0.009, p=0.037). For the predictive index of dose variation, CVF ratio was to be statistically significant for PTV in the liver dome versus PTV (SRS r=0.684, 3DCRT r=0.732, p<0.01) and CVF ratio for Tumor size was to be statistically significant (SRS r=-0.193, p=0.017, 3DCRT r=0.237, p=0.023).

• The Journal of Korean Society for Radiation Therapy
• /
• v.24 no.1
• /
• pp.11-14
• /
• 2012

Purpose: Anisotropic Analytical Algorithm (AAA) provides more accurate dose calculation regarding impact on scatter and tissue inhomogeneity in comparison to Pencil Beam Convolution (PBC) algorithm. This study tries to analyze the difference of dose distribution according to PBC algorithm and dose calculation algorithm of AAA on breast cancer tangential plan. Materials and Methods: Computerized medical care plan using Eclipse treatment planning system (version 8.9, VARIAN, USA) has been established for the 10 breast cancer patients using 6 MV energy of Linac (CL-6EX, VARIAN, USA). After treatment plan of Conventional Radiation Therapy plan (Conventional plan) and Field-in-Field Intensity Modulated Radiation Therapy plan (FiF plan) using PBC algorithm has been established, MU has been fixed, implemented dose calculation after changing it to AAA, and compared and analyzed treatment plan using Dose Volume Histogram (DVH). Results: Firstly, as a result of evaluating PBC algorithm of Conventional plan and the difference according to AAA, the average difference of CI value on target volume has been highly estimated by 0.295 on PBC algorithm and as a result of evaluating dose of lung, $V_{47Gy}$ and $V_{45Gy}$ has been highly evaluated by 5.83% and 4.04% each, Mean dose, $V_{20Gy}$, $V_{5Gy}$, $V_{3Gy}$ has been highly evaluated 0.6%, 0.29%, 6.35%, 10.23% each on AAA. Secondly, in case of FiF plan, the average difference of CI value on target volume has been highly evaluated on PBC algorithm by 0.165, and dose on ipsilateral lung, $V_{47Gy}$, $V_{45Gy}$, Mean dose has been highly evaluated 6.17%, 3.80%, 0.15% each on PBC algorithm, $V_{20Gy}$, $V_{5Gy}$, $V_{3Gy}$ has been highly evaluated 0.14%, 4.07%, 4.35% each on AAA. Conclusion: When calculating with AAA on breast cancer tangential plan, compared to PBC algorithm, Conformity on target volume of Conventional plan, FiF plan has been less evaluated by 0.295, 0.165 each. For the reason that dose of high dose region of ipsilateral lung has been showed little amount, and dose of low dose region has been showed much amount, features according to dose calculation algorithm need to be considered when we evaluate dose for the lungs.

• Radiation Oncology Journal
• /
• v.20 no.3
• /
• pp.274-282
• /
• 2002

Purpose : Measurement of transmission dose is useful for in vivo dosimetry. In this study, the algorithm for estimating the transmission dose for open radiation fields was modified for application to partially blocked radiation fields. Materials and Methods : The beam data was measured with a flat solid phantom with various blocked fields. A new correction algorithm for partially blocked radiation field was developed from the measured data. This algorithm was tested in some settings simulating clinical treatment with an irregular field shape. Results : The correction algorithm for the beam block could accurately reflect the effect of the beam block, with an error within ${\pm}1.0\%$, with both square fields and irregularly shaped fields. Conclusion : This algorithm can accurately estimate the transmission dose in most radiation treatment settings, including irregularly shaped field.

• Progress in Medical Physics
• /
• v.23 no.1
• /
• pp.48-53
• /
• 2012

The pencil beam convolution (PBC) algorithms in radiation treatment planning system have been widely used to calculate the radiation dose. A new photon dose calculation algorithm, referred to as the anisotropic analytical algorithm (AAA), was released for use by the Varian medical system. The aim of this paper was to investigate the difference in dose calculation between the AAA and PBC algorithm using the intensity modulated radiation therapy (IMRT) plan for lung cancer cases that were inhomogeneous in the low density. We quantitatively analyzed the differences in dose using the eclipse planning system (Varian Medical System, Palo Alto, CA) and I'mRT matirxx (IBA, Schwarzenbruck, Germany) equipment to compare the gamma evaluation. 11 patients with lung cancer at various sites were used in this study. We also used the TLD-100 (LiF) to measure the differences in dose between the calculated dose and measured dose in the Alderson Rando phantom. The maximum, mean, minimum dose for the normal tissue did not change significantly. But the volume of the PTV covered by the 95% isodose curve was decreased by 6% in the lung due to the difference in the algorithms. The difference dose between the calculated dose by the PBC algorithms and AAA algorithms and the measured dose with TLD-100 (LiF) in the Alderson Rando phantom was -4.6% and -2.7% respectively. Based on the results of this study, the treatment plan calculated using the AAA algorithms is more accurate in lung sites with a low density when compared to the treatment plan calculated using the PBC algorithms.