• 한국방사선학회논문지
• /
• 제17권5호
• /
• pp.633-640
• /
• 2023

최근 근접 치료에서 방사선 차폐막을 사용하여 선량 분포를 변조하여 선량을 전달하는 정적 및 동적 변조 근접 치료 방법이 개발됨에 따라 새로운 방향성 빔 세기 변조 근접 치료에 적합한 역방향 치료 계획 및 치료 계획 최적화 알고리즘에서 선량 계산에 필요한 파라미터 및 데이터의 양이 증가하고 있다. 세기 변조 근접 치료는 방사선의 정확한 선량 전달이 가능하지만, 파라미터와 데이터의 양이 증가하기 때문에 선량 계산에 필요한 경과 시간이 증가한다. 본 연구에서는 선량 계산 경과 시간의 증가를 줄이기 위해 그래픽 카드 기반의 CUDA 가속 선량 계산 알고리즘을 구축하였다. 계산 과정의 가속화 방법은 관심 체적의 시스템 행렬 계산 및 선량 계산의 병렬화를 이용하여 진행하였다. 개발된 알고리즘은 모두 인텔(3.7GHz, 6코어) CPU와 단일 NVIDIA GTX 1080ti 그래픽 카드가 장착된 동일한 컴퓨팅 환경에서 수행하였으며, 선량 계산 시간은 디스크에서 데이터를 불러오고 전처리를 위한 작업 등의 추가 적으로 필요한 시간은 제외하고 선량 계산 시간만 측정하여 평가하였다. 그 결과 가속화된 알고리즘은 CPU로만 계산할 때보다 선량 계산 시간이 약 30배 단축된 것으로 나타났다. 가속화된 선량 계산 알고리즘은 적응방사선치료와 같이 매일 변화되는 어플리케이터의 움직임을 고려하여 새로운 치료 계획을 수립해야 하는 경우나 동적 변조 근접 치료와 같이 선량 계산에 변화되는 파라미터를 고려해야 하는 경우 치료 계획 수립 속도를 높일 수 있을 것으로 판단된다.

• 대한방사선기술학회지:방사선기술과학
• /
• 제43권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.

• 한국의학물리학회지:의학물리
• /
• 제31권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.

• Radiation Oncology Journal
• /
• 제2권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.

• 대한방사선기술학회지:방사선기술과학
• /
• 제43권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.

• 대한방사선치료학회지
• /
• 제26권2호
• /
• pp.171-176
• /
• 2014

목 적 : Lung SABR plan 에서 AAA의 calculation grid를 변화시켜 선량변화를 분석하고 그에 따른 영향을 연구하여 적절한 적용 방안에 대해 고찰한다. 대상 및 방법 : 모든 plan에 이용된 4D CT image는 Brilliance Big Bore CT(Philips, Netherlands)에서 촬영되었으며 10 건의 Lung SABR plan($Eclipse^{TM}$ ver 10.0.42, Varian, the USA)에서 anisotropic analytic algorithm (AAA, ver. 10, Varian Medical Systems, Palo Alto, CA, USA)을 이용하여 각각 1.0, 3.0, 5.0 mm의 calculation grid로 계산하였다. 결 과 : 10 건의 Lung SABR plan에서 1.0 mm calculation grid를 사용한 경우 $V_{98}$이 각각 처방선량의 약 $99.5{\pm}1.5%$ 였으며 Dmin이 각각 처방선량의 약 $92.5{\pm}1.5%$ 였고 Homogeneity Index(HI)는 약 $1.0489{\pm}0.0025$로 나타났다. 3.0 mm calculation grid를 사용한 경우 $V_{98}$이 각각 처방선량의 약 $90{\pm}4.5%$였으며, Dmin이 각각 처방선량의 약 $87.5{\pm}3%$ 였고 HI가 약 $1.07{\pm}1$로 나타났다. 5.0 mm calculation grid를 사용한 경우 $V_{98}$이 각각 처방선량의 약 $63{\pm}15%$ 였으며, Dmin이 각각 처방선량의 약 $83{\pm}4%$ 였고 HI가 약 $1.13{\pm}0.2$로 나타났다. 결 론 : 1.0 mm calculation grid의 계산 시간이 3.0 mm, 5.0 mm 보다 오래 걸렸지만 grid의 간격이 좁을수록 상대적으로 작은 PTV를 갖는 plan의 정확성을 향상시키는 것으로 나타났다. 또한 Lung과 같이 비교적 넓게 퍼져 있으며 밀도가 낮은 장기의 작은 PTV를 치료해야 하는 경우에는 1.0 mm의 calculation grid를 사용하는 것이 좋을 것으로 사료된다.

• 대한방사선기술학회지:방사선기술과학
• /
• 제39권3호
• /
• pp.305-312
• /
• 2016

본 연구에서는 진단용 X선 검사에서 환자에게 피폭되는 두부 및 사지를 다양한 선량 계산법을 통해 실측 선량과 비교 실험하였다. 또한 촬영 장비의 형태, 장비 설정조건, X선의 용량, X선관과 환자와의 거리, X선 후방산란차이 등을 고려한 새로운 계산 방법을 제시하여 피폭선량을 산출하였다. 그 결과 피부입사선량이 기존의 선량 계산법보다 실측과의 오차가 줄어들었으며, 환자가 피폭되는 선량을 쉽게 계산할 수 있었고 의료선량 평가가 이루어지게 되어 방사선 관련 종사자들의 의료 선량 관리가 더욱 수월해지는 계기가 될 것으로 사료된다.

• Nuclear Engineering and Technology
• /
• 제36권1호
• /
• pp.1-11
• /
• 2004

Two approximate methods for a cosmic radiation shielding calculation in low earth orbits were developed and assessed. Those are a sectoring method and a chord-length distribution method. In order to simulate a change in cosmic radiation environments along the satellite mission trajectory, IGRF model and AP(E)-8 model were used. When the approximate methods were applied, the geometrical model of satellite structure was approximated as one-dimensional slabs, and a pre-calculated dose-depth conversion function was introduced to simplify the dose calculation process. Verification was performed with mission data of KITSAT-1 and the calculated results were also compared with detailed 3-dimensional calculation results using Monte Carlo calculation. Dose results from the approximate methods were conservatively higher than Monte Carlo results, but were lower than experimental data in total dose rate. Differences between calculation and experimental data seem to come from the AP-8 model, for which it is reported that fluxes of proton are underestimated. We confirmed that the developed approximate method can be applied to commercial satellite shielding calculations. It is also found that commercial products of semi-conductors can be damaged due to total ionizing dose under LEO radiation environment. An intensive shielding analysis should be taken into account when commercial devices are used.

• Nuclear Engineering and Technology
• /
• 제54권5호
• /
• pp.1769-1780
• /
• 2022

A new method of dose calculation algorithm, called GPU-accelerated Monte Carlo and collapsed cone Convolution (GMCC) was developed to improve the calculation speed of BNCT treatment planning system. The GPU-accelerated Monte Carlo routine in GMCC is used to simulate the neutron transport over whole energy range and the Collapsed Cone Convolution method is to calculate the gamma dose. Other dose components due to alpha particles and protons, are calculated using the calculated neutron flux and reaction data. The mathematical principle and the algorithm architecture are introduced. The accuracy and performance of the GMCC were verified by comparing with the FLUKA results. A water phantom and a head CT voxel model were simulated. The neutron flux and the absorbed dose obtained by the GMCC were consistent well with the FLUKA results. In the case of head CT voxel model, the mean absolute percentage error for the neutron flux and the absorbed dose were 3.98% and 3.91%, respectively. The calculation speed of the absorbed dose by the GMCC was 56 times faster than the FLUKA code. It was verified that the GMCC could be a good candidate tool instead of the Monte Carlo method in the BNCT dose calculations.

• 한국의학물리학회:학술대회논문집
• /
• 한국의학물리학회 2003년도 제27회 추계학술대회
• /
• pp.41-41
• /
• 2003

Purpose: Even if the wedge filter is widely used for the radiation therapy to modify the photon beam intensity, the wedged photon beam dose calculation is not so easy. Radiation therapy planning systems (RTPS) have been used the empirical or semi-analytical methods such as attenuation method using wedge filter parameters or wedge filter factor obtained from measurement. However, these methods can cause serious error in penumbra region as well as in edge region. In this study, we propose the dose calculation algorithm for wedged field to minimize the error especially in the outer beam region. Materials and Method: Modified intensity by wedge filter was calculated using tissue-maximum ratio (TMR) and scatter-maximum ratio (SMR) of wedged field. Profiles of wedged and non-wedged direction was also used. The result of new dose calculation was compared with measurement and the result from attenuation method. Results: Proposed algorithm showed the good agreement with measurement in the high dose-gradient region as well as in the inner beam region. The error was decreased comparing to attenuation method. Conclusion: Although necessary beam data for the RTPS commissioning was increased, new algorithm would guarantee the improved dose calculation accuracy for wedged field. In future, this algorithm could be adopted in RTPS.