Multi-Purpose Particle and Heavy Ion Transport code System

Slides:



Advertisements
Similar presentations
P HI T S 陽子ビームで雪だるまを溶かそう Multi-Purpose Particle and Heavy Ion Transport code System title 年 3 月改訂.
Advertisements

EGS5 の導入 KEK 波戸芳仁 Last modified on 例題1 ベータ線を物質に打ち込 む ベータ線は物質で止まってしまうか?通 り抜けるか? 物質の内部でどのような反応が起こる か? ベータ線.
ユーザーコードの導入 2010 年 7 月 20 日 KEK 波戸. 例題1 ベータ線を物質に打ち込 む ベータ線 ベータ線は物質で止まってしまうか?通り抜けるか? 物質の内部でどのような反応が起こるか?
Image J を用いた粒子径分析 目次 1.ImageJ ダウンード・開始 2. 粒子径分析の基本操作.
P HI T S α 線, β 線, γ 線,中性子線を止めるに は? Multi-Purpose Particle and Heavy Ion Transport code System title 年 5 月改訂.
P HI T S PSFC4PHITS の使い方 Multi-Purpose Particle and Heavy Ion Transport code System Title 年 2 月改訂.
P HI T S スクリプト言語を用いた PHITS の連続 実行 Multi-Purpose Particle and Heavy Ion Transport code System title 年 2 月改訂.
計算効率を上げるためのvariance reduction(粒子のウェイト)の利用
7-1.WEKOコンテンツ 一括登録 マニュアル Version2.5
Cornell Spectrum Imagerの使い方
ボクセルファントムを用いた X線治療のシミュレーション
PHITS 講習会 基礎実習(III): 計算条件の設定
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
PHITS講習会 基礎実習(I): 体系及び線源の定義
PHITS講習会 基礎実習(II): Tally(タリー)の定義
このスペクトルをPC上で処理したい 1.
強度変調回転放射線治療(VMAT)の シミュレーション
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
第3回 ファイルとフォルダ 伊藤 高廣 計算機リテラシーM 第3回 ファイルとフォルダ 伊藤 高廣
DICOM医療画像を使用したPHITSシミュレーション
デジタルポートフォリオ作成支援ツール PictFolio 使用マニュアル
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
PHITSによるX線治療シミュレーション 基本操作復習編
ParaViewを用いたPHITS 計算結果の3次元表示
EGSに対応した粒子軌跡と 計算体系の3次元表示ソフト - CGVIEW -
XAS マニュアル 暫定版.
X線天文衛星用CCDカメラの 放射線バックグランドの評価
スクリプト言語を用いたPHITSの連続実行
PHITS 講習会 基礎実習(III): 計算条件の設定
誘導放射能の評価方法: PHITSとDCHAIN-SPの接続計算
α線,β線,γ線,中性子線を止めるには?
京都大学大学院医学研究科 画像応用治療学・放射線腫瘍学 石原 佳知
計算効率を上げるためのvariance reduction(粒子のウェイト)の利用 A
X線CTにおけるファントム中の エネルギー変化についての検討
Multi-Purpose Particle and Heavy Ion Transport code System
応用実習用資料 Environmental radioactivity
Multi-Purpose Particle and Heavy Ion Transport code System
東京海洋大産学官連携研究員/技術コンサルタント 高須 知二 Tomoji TAKASU
MeV internal meeting Oct. 2, 2015
計算効率を上げるためのvariance reduction(粒子のウェイト)の利用 A
PHITS 講習会 基礎実習(III): 計算条件の設定
応用実習用資料 Neutron target
実 習 4 2次元テーブルの利用.
IAEA phase space fileを用いた X線治療シミュレーション
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
治療用フィルムによる線量分布測定の 基礎的検討Ⅱ
Multi-Purpose Particle and Heavy Ion Transport code System
EGSに対応した粒子軌跡と 計算体系の3次元表示ソフト - CGVIEW -
フレアの非熱的成分とサイズ依存性    D1 政田洋平      速報@太陽雑誌会(10/24).
COBOLプログラムの概観 プログラミング資料.
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
μ+N→τ+N反応探索実験の ためのシミュレーション計算
ユーザーコードに記述する事項の概要 2010年7月21日 KEK 波戸.
α線,β線,γ線,中性子線を止めるには?
Multi-Purpose Particle and Heavy Ion Transport code System
平山 英夫、波戸 芳仁 KEK, 高エネルギー加速器研究機構
BLJ2013 BentleyArchitecture
Simulation study for drift region
計算効率を上げるためのvariance reduction(粒子のウェイト)の利用 B
Molecular Devices Japan
OPEN文、READ文、PERFORM文
目的とするユーザーコードを 作成するために
Presentation transcript:

Multi-Purpose Particle and Heavy Ion Transport code System PHITS Multi-Purpose Particle and Heavy Ion Transport code System X線治療シミュレーション実習 PHITS講習会 入門実習 2017年2月改訂 Title 1

実習内容 頭部ファントム(模型)をCT(スライス厚2mm)で撮影したDICOM画像を用いて、X線治療シミュレーションをPHITSで実行し、治療計画装置の結果と比較する 照射野8cm x 8 cm の4MV X線 Tumor 頭部ファントム Contents 2

DICOMイメージの読み込み dicom2phits.inpを以下のように変更してDICOM2PHITSを実行 DICOM2PHITS input "data/HumanVoxelTable.data" ! File for conversion of human voxel data "work_Kanazawa/DICOM/” ! DICOM file directory "work_Kanazawa/PHITSinput/" ! Filename for PHITS input 1 120 ! Minimum slice number, Maximum slice number 110 400 120 430 ! Clipping: Nxmin, Nxmax, Nymin, Nymax 4 4 2 ! Coarse graining: Nxc, Nyc, Nzc 1 ! Coordinates 0:Zero at center 1:DICOM system 1 ! PHITS parameter: 1:PhotonTherapy 2:ParticleTherapy [ Source ] s-type = 4 $ cylindrical source with energy spectrum … proj = photon $ Varian 6MV photon spectrum based on D. Sheikh-Bagheri et al. PMB 29 (2002) 391 e-type = 1 ne = 24 0.0000E+00 2.1400E-03 2.5000E-01 1.2600E-02 5.0000E-01 1.3100E-02 6MV X線スペクトル(初期設定) DICOM2PHITS 3

巨大なボクセルデータをあらかじめバイナリー化して読込時間短縮! 読込の高速化 目的 PHITSでは一度インプットファイルを全てバイナリー化してから再読込 巨大なボクセルデータをあらかじめバイナリー化して読込時間短縮! 手順 ① [Parameters]セクションのivoxelを有効にする($を消す) ivoxel = 2 # LatticeのFill部分をバイナリー化としてfile(18)に出力させるオプション file(18) = voxel.bin # 出力するバイナリーファイルのファイル名 ② PHITSを実行する → Binary file was successfully generated!!  ③ ivoxel = 1に変更し,Latticeを定義するinflコマンドをコメントアウト ivoxel = 1 # LatticeのFill部分をfile(18)から読み込むオプション 高速化! $ infl:{voxel.inp} 112行目 Binary generation 4

詳しくはマニュアルのEGS5用パラメータの項目を参照 マテリアル数が多い場合の注意点 マテリアル数が多い場合にEGSを使用すると、EGS用断面積データを用意するの に非常に時間がかかる(PHITS実行の最初に全てのマテリアルについて用意)。 複数回PHITSを実行する場合には、PHITS実行時にEGS用断面積データを毎回 準備するのではなく、一度作成した断面積データを使いまわすことができる。 実行方法 ① ipegs=-1 をパラメータセクションで設定し、PHITSを実行    ⇒ EGS用断面積データを生成し、輸送計算を実行せずに終了 ② ipegs=2 をパラメータセクションで設定し、PHITSを実行    ⇒ 既存のEGS用断面積データを使用して、輸送計算を実行 詳しくはマニュアルのEGS5用パラメータの項目を参照 Many materials with EGS 5

PHITS実行による体型の確認 icntl=8 deposit-xy.eps deposit-xz.eps icntl=0 Geometry check 6

線源の変更1 Adjustment of source 7 円柱スペクトル線源 直方体スペクトル線源 [ Transform ] $ Transform source directing beam along x-axis tr400 -0.29115 -10.88558 0.00000 … [ Source ] s-type = 4 $ Create radiation field centered on 0 $ Adjust center of radiation field in trcl=400 x0 = 0.0 y0 = 0.0 z0 = 0.0 z1 = 0.0 r0 = 5.0 dir = 1.0 trcl = 400 [ Transform ] $ Transform source directing beam along x-axis tr400 -0.29115 -10.88558 0.00000 … [ Source ] s-type = 5 $ Create radiation field centered on 0 $ Adjust center of radiation field in trcl=400 x0 = -0.049-4.0 x1 = -0.049+4.0 y0 = -1.266-4.0 y1 = -1.266+4.0 z0 = -0.3 z1 = -0.3 dir = 1.0 trcl = 400 円柱スペクトル線源 直方体スペクトル線源 Trcl=400でy軸向きにビームを回転 (-0.049, -1.266, -0.3)[cm] を中心に、xy平面で8x8cm2の照射野を設定 icntl=0でPHITSを実行し、deposit-xz.epsでビーム形状を確認 Adjustment of source 7

線源の変更2 Adjustment of source 8 4MV X線スペクトル に変更する 6MV X線スペクトル s-type = 4 … $ Varian 6MV photon spectrum based on D. Sheikh-Bagheri et al. PMB 29 (2002) 391 e-type = 1 ne = 24 0.0000E+00 2.1400E-03 2.5000E-01 1.2600E-02 5.0000E-01 1.3100E-02 7.5000E-01 1.1400E-02 1.0000E+00 9.7600E-03 [ Source ] s-type = 4 … $ Varian 6MV photon spectrum based on D. Sheikh-Bagheri et al. PMB 29 (2002) 391 e-type = 1 ne = 16 infl:{spectrum-4MV.inp} 4MV X線スペクトル に変更する 6MV X線スペクトル (初期設定) 初期設定のデータは削除 Adjustment of source 8

三次元線量分布出力と変換 PHITS2DICOM 9 ①PHITSで三次元線量分布を計算 ・ phits.inpの一つ目と二つ目のt-depositタリーにOFF、最後のt-depositタリーのOFFを削除 ・ icntl=0でPHITSを実行 DICOMフォルダにあるファイルをどれでも良いので1つ指定 ②入力ファイルを作成(phits2dicom.inp) PHITS2DICOM input PHITS2DICOM用のインプットであることを明示 data/sample.dcm DICOM RT-doseサンプルファイルを指定 work_Kanazawa/DICOM/2.16.840…………….dcm DICOM2PHITSの変換に使用したDICOMファイル work_Kanazawa/PHITSinput/deposit-3D.out PHITSで計算した三次元線量分布データを指定 work_Kanazawa/PHITSinput/phits.out PHITS実行で出力されたphits.outを指定 work_Kanazawa/PHITSinput/RTD.deposit-3D.dcm RT-dose形式で出力するファイル名を指定 0 線量の規格化を設定 0:しない(最大値が1で規格化) 1:する ③ PHITS2DICOM実行 (Windows) dicom2phits_win.batにphits2dicom.inpをドラッグ&ドロップ     (Mac) dicom2phits_mac.commnadをダブルクリック、現れる窓にphits2dicom.inpと入力 ④ DICOM viewer (例:dicompyler)で確認 ・ work_Tokushimaのフォルダを指定 ・ Rx値は1 cGyを選択 PHITS2DICOM 9

三次元線量分布 統計量が足りず、滑らかな分布が得られていない maxcasを10倍に増やして確認してみましょう 3D dose distribution 10

Improvement of statistics 統計精度 10倍程度ではまだ分布が滑らかにならない ⇒ 1000倍以上の計算量が必要 phits2dicom.inpを次の様に変更し、PHITS2DICOMを実行 PHITS2DICOM input PHITS2DICOM用のインプットであることを明示 data/sample.dcm DICOM RT-doseサンプルファイルを指定 work_Kanazawa/DICOM/2.16.840…………….dcm DICOM2PHITSの変換に使用したDICOMファイル work_Kanazawa/PHITSinput/LONGdeposit-3D.out PHITSで計算した三次元線量分布データを指定 work_Kanazawa/PHITSinput/LONGphits.out PHITS実行で出力されたphits.outを指定 work_Kanazawa/PHITSinput/RTD.deposit-3D.dcm RT-dose形式で出力するファイル名を指定 0 線量の規格化を設定 0:しない(最大値が1で規格化) 1:する 1000倍の統計量の計算結果 DICOM viewer で確認 (Rx値は100cGyを選択) 分布が滑らかに Improvement of statistics 11

治療計画装置の出力との比較 DICOM viewer でTP_Kanazawaフォルダを指定 規格化は? RT-imageと治療計画装置の RT-plan, RT-doseを内包 規格化は? Comparison with TP 12

治療計画装置の出力の規格化 DICOM Treeのタグ選択 プルダウンからRT Plan選択 規格化の位置 規格値 Comparison with TP 13

線量分布比較のための規格化 Comparison with TP 14 phits2dicom.inpを次の様に変更し、PHITS2DICOMを実行 PHITS2DICOM input PHITS2DICOM用のインプットであることを明示 data/sample.dcm DICOM RT-doseサンプルファイルを指定 work_Kanazawa/DICOM/2.16.840…………….dcm DICOM2PHITSの変換に使用したDICOMファイル work_Kanazawa/PHITSinput/LONGdeposit-3D.out PHITSで計算した三次元線量分布データを指定 work_Kanazawa/PHITSinput/LONGphits.out PHITS実行で出力されたphits.outを指定 work_Kanazawa/PHITSinput/RTD.deposit-3D.dcm RT-dose形式で出力するファイル名を指定 1 線量の規格化を設定 0:しない(最大値が1で規格化) 1:する -0.5 -12.6 -2.0 規格化する位置座標(x,y,z) [mm] 2.0 上記位置での線量値 治療計画装置のRT-planヘッダーから読み取った情報 DICOM viewer で比較 (Rx値は200cGyを選択)   治療計画装置の出力結果 (TP_Kanazawa) PHITSシミュレーション結果 (work_Kanazawa) Comparison with TP 14

線量分布の比較 Z=-3mm Z=-3mm 治療計画装置の出力 (TP_Kanazawa) PHITSシミュレーション (work_Kanazawa) Comparison with TP 15

(DICOM2PHITS→PHITS→PHITS2DICOM) まとめ DICOMイメージを使用したPHITSシミュレーション (DICOM2PHITS→PHITS→PHITS2DICOM) 線源を状況に合わせて変更 PHITS2DICOMを使用することで、PHITS出力の三次元線量分布をRT-dose形式に変換 十分な統計精度を得るためには入射粒子数を増やす必要あり 線量分布の比較のためには規格化が重要 Summary 16