JP3689610B2

JP3689610B2 - Ridge filter design method

Info

Publication number: JP3689610B2
Application number: JP2000006562A
Authority: JP
Inventors: たか文中川
Original assignee: Mitsubishi Electric Corp
Current assignee: Mitsubishi Electric Corp
Priority date: 2000-01-14
Filing date: 2000-01-14
Publication date: 2005-08-31
Anticipated expiration: 2020-01-14
Also published as: JP2001194499A

Description

【０００１】
【発明の属する技術分野】
本発明は、透過する粒子線の線量を変化させてターゲット内に所定の線量分布を作るリッジフィルタの設計方法に関する。
【０００２】
【従来の技術】
癌を治療する方法の１つに、粒子線を癌の腫瘍に照射する治療方法がある。この治療方法は、癌腫瘍へ粒子線を照射することにより癌細胞を死滅させるものである。図１３は、粒子線１０を照射される腫瘍１３０の例を示す。粒子線１０は以下のように照射される。腫瘍１３０は粒子線１０に垂直な平面上にある程度の広がりを持つので、粒子線１０は、コリメータ１３２および補正用媒質１３４によりその広がりに応じて整形される。腫瘍１３０以外の正常な組織へ粒子線１０を照射しないようにするためである。腫瘍１３０は深さ方向にも厚みＤを持つので、その厚みＤの部分に線量が所望の値で一様な粒子線１０を照射する必要がある。
【０００３】
深さ方向の所定の範囲で線量が一様な粒子線１０は、複数の異なるエネルギーを持つ粒子線から得られる。例えば、図２の（ｂ）は、水ターゲットに入射した粒子線の線量分布（ブラッグ曲線）を示すところ、線量が一様な領域ＳＯＢＰ（Spread-Out Bragg Peak）は、ブラッグ曲線２２、２４’’のそれぞれを重ね合わせることにより得られる。ブラッグ曲線２２、２４’’のそれぞれは、図２の（ａ）におけるブラッグ曲線２２、２４に重みをかけたものである。
【０００４】
このようなＳＯＢＰを得るための複数の異なるエネルギーを持つ粒子線は、リッジフィルタを用いて生成される。リッジフィルタは、入射したあるエネルギーを持つ粒子線を、それぞれ重みをかけた複数の単色エネルギーの粒子線に変換するフィルタである。すなわちリッジフィルタは、ブラッグ曲線２２で表されるエネルギーの粒子線を受け、ブラッグ曲線２４’’によって表される複数の単色エネルギーの粒子線に変換する。これらの粒子線が重ね合わされることにより、リッジフィルタから出力された粒子線は、照射対象（ターゲット）に対してブラッグ曲線２６で表されるＳＯＢＰを持つ粒子線を形成する。図３は、リッジフィルタの例を示す。リッジフィルタの入射面３０からの階段の高さｈに応じて、粒子線は異なるエネルギーを持つ複数の粒子線に変換される。すなわち、階段の高さｈが小さいほどリッジフィルタ内で減衰する量が小さいので、深い位置まで到達できるエネルギーを持つ粒子線が透過される。よって、階段の高さｈは図２のグラフにおける横軸方向のブラッグ曲線の位置を規定する。さらに階段の幅ｗを調整して通過する粒子数を制御することにより、線量レベルを調整することができる。よって、階段の幅ｗは、図２のグラフにおける縦軸方向のブラッグ曲線の重みに相当する。換言すれば、この階段の幅ｗは粒子束に相当する。
【０００５】
粒子線の進行距離に関して所望の距離で、かつ所望の線量値で一様なＳＯＢＰを得るために、適切な階段の高さｈおよび幅ｗを持つリッジフィルタを設計する必要がある。従来のリッジフィルタの設計手順は、例えば「MONTE CARLO SIMULATION OF A PROTONTHERAPY SYSTEM FOR THE CALCULATION OF THE DOSE DISTRIBUTION IN A PATIENT」TERA, 95, No.4, pp25 (1995)に記載されている。図１４は、従来のリッジフィルタの設計フローを示す。（１）まず目標とするＳＯＢＰ幅を決定する。（２）異なるエネルギーの複数のブラッグ曲線（図２のブラッグ曲線２４）を、モンテカルロ計算により求める（ステップＳ１４２）。（３）目標とするＳＯＢＰ上の評価点と、その位置での目標値に対して、全てのブラッグ曲線の重みｘを、式Ａｘ＝ｂから決定する（ステップＳ１４４）。ここで、行列Ａには各ブラッグ曲線２４上のサンプル値を格納している。その１行目は最大エネルギーの値、第２行は２番目のエネルギーの値等であり、ｂは目標値である。（４）求めた重みｘに基づいてリッジフィルタの階段の幅ｗを決定する（ステップＳ１４６）。このとき、階段の高さの差（階段の各ステップの高さ）は一定である。
【０００６】
【発明が解決しようとする課題】
上の（２）で述べたように、従来は多数回モンテカルロ計算を行い、エネルギーの異なる複数のブラッグ曲線２４（図２（ａ））を得なければならなかった。一般にモンテカルロ計算には時間がかかるため、これでは設計を迅速に行うことができない。
【０００７】
さらに、従来の方法では目標とするＳＯＢＰ上の評価点の位置によっては、平坦度のよいＳＯＢＰが得られないという問題がある。例えば、評価点は飛び飛びの位置を採用するため評価点間のＳＯＢＰ値の挙動が考慮されず、したがって評価点以外の点で目標値との差が大きくなる場合がある。特に最大エネルギーのブラッグ曲線のピークが評価点として採用されない場合には重みの影響が大きく、ＳＯＢＰグラフの右端に突起が出ることがある（後述する図５の横軸２８５ｍｍ付近の値）。
【０００８】
一旦作成したリッジフィルタを使って理論上の値を計算したとき、理論値と実際に必要な値とが一致しない場合は適当な補償方法がなく、再度設計を行わなければならなかった。または実測値に基づいた複数のブラッグ曲線を得て、それに基づいてリッジフィルタを設計する必要があった。
【０００９】
本発明の目的は、均一なＳＯＢＰを実現するリッジフィルタを迅速に設計する方法を提供することである。
【００１０】
【課題を解決するための手段】
本発明のリッジフィルタの設計方法は、粒子線の進行距離と線量との関係を表すブラッグ曲線を用いて前記粒子線が透過する透過部の高さと幅とを調整し、透過後の前記線量が前記進行距離の所定の範囲において所望の値で一様なＳＯＢＰを実現するリッジフィルタの設計方法であって、モンテカルロ計算により、前記粒子線の第１のブラッグ曲線を取得するステップと、前記第１のブラッグ曲線をｎ回（ｎ：自然数）平行移動することにより、ｎ本の第２のブラッグ曲線を取得するステップと、前記ＳＯＢＰの前記所望の値と、前記第１のブラッグ曲線と、前記第２のブラッグ曲線とに基づいて前記第１のブラッグ曲線と前記第２のブラッグ曲線の各々とにかける重みを計算するステップであって、前記重みをかけた前記第１のブラッグ曲線と前記第２のブラッグ曲線の各々とを重ね合わせると、前記所望の値の前記ＳＯＢＰが得られるよう重みを計算するステップと、前記重みに基づいて前記透過部の幅を決定するステップと、からなる、リッジフィルタ設計方法であり、これにより上記目的が達成される。
【００１１】
前記リッジフィルタ設計方法は、前記重みに基づいてモンテカルロ計算を行い、前記ＳＯＢＰの理論値を計算するステップと、前記理論値と、前記所望の値との差を計算するステップと、前記差と、前記第１のブラッグ曲線と、前記第２のブラッグ曲線とに基づいて前記差を補正するための前記第１のブラッグ曲線と前記第２のブラッグ曲線の各々とにかける補正重みを計算するステップであって、前記補正重みをかけた前記第１のブラッグ曲線と前記第２のブラッグ曲線の各々とを、さらに前記重みをかけた前記第１のブラッグ曲線と前記第２のブラッグ曲線の各々とに重ね合わせると、前記差に相当する前記線量が最も小さくなるよう前記補正重みを計算するステップと、をさらに含み、前記決定するステップは、さらに前記補正重みに基づいて前記透過部の幅を決定してもよい。
【００１２】
前記リッジフィルタ設計方法は、前記重みに基づいてモンテカルロ計算を行い、前記ＳＯＢＰの理論値を計算するステップと、前記所望の値と、前記理論値との差を計算するステップと、前記所望の値に前記差を加えた値と、前記第１のブラッグ曲線と、前記第２のブラッグ曲線とに基づいて、前記重みを修正した修正重みを計算するステップであって、前記修正重みをかけた前記第１のブラッグ曲線と前記第２のブラッグ曲線の各々とを重ね合わせると、前記所望の値の前記ＳＯＢＰが得られるよう重みを計算するステップと、をさらに含み、前記決定するステップは、前記修正重みに基づいて前記透過部の幅を修正してもよい。
【００１３】
前記リッジフィルタ設計方法は、前記重みと前記補正重みとに基づいて得られるＳＯＢＰの値と、前記ＳＯＢＰの前記所望の値との誤差を計算するステップと、前記誤差が所定の許容範囲を超えている位置を選択するステップと、をさらに含み、選択された前記位置の誤差を小さくするように、前記重みを計算するステップと前記補正重みを計算するステップとを再度実行してもよい。
【００１４】
前記リッジフィルタは、前記粒子線の入射方向に平行な高さと前記粒子線の入射方向に垂直な幅とを有する前記透過部を備えたリッジフィルタであってもよい。
【００１５】
前記リッジフィルタは、前記入射方向と平行な軸を中心として回転する、少なくとも１つの前記透過部を備えた回転形リッジフィルタであり、前記透過部の幅は、前記透過部の回転方向の幅を変化させることによって調整されてもよい。
【００１６】
前記回転形リッジフィルタは複数の前記透過部を有し、前記複数の透過部の各々は、それぞれ異なる高さを有していてもよい。
【００１７】
【発明の実施の形態】
（実施の形態１）
以下では、本発明により設計されるリッジフィルタを利用する装置の説明をした後、そのリッジフィルタを設計するための手順を説明する。本明細書において粒子線とは、例えば陽子線である。また、粒子線を照射する対象は深さ方向に厚みＤを持つ腫瘍１３０（図１３）とする。なお「深さ方向」とは粒子線が照射され、進行する方向をいう。
【００１８】
図１は、リッジフィルタを用いたワブラー法による粒子線照射装置１の模式図である（「Instrumentation for treatment of cancer using proton and light-iron beams」Rev. Sci. Instrum. Vol. 64 No.8 pp.2055-2091 (1993)）。粒子線照射装置１は、腫瘍１３０（図１３）の形状に応じて粒子線１０を照射する領域（深さ方向も含む）を調整し、必要な線量の粒子線１０を腫瘍１３０（図１３）に照射する装置である。粒子線照射装置１は、１対の水平方向揺動用電磁石１１と、１対の垂直方向揺動用電磁石１２と、粒子のエネルギーを変えるためのレンジシフタ１３と、リッジフィルタ１４と、照射ボリューム１６とを含む。粒子線源（図示せず）から放射された粒子線１０は、１対の水平方向揺動用電磁石１１と、１対の垂直方向揺動用電磁石１２とを用いて揺動され、照射面積が広げられる。揺動された粒子線１０はレンジシフタ１３を通り、粒子線１０のエネルギーが適当な量だけ小さくされる。その後、粒子線１０はリッジフィルタ１４に入射する。リッジフィルタ１４は、透過した媒質の長さに応じて粒子線のエネルギーが減衰することを利用して、入射した粒子線をそれぞれ重みをかけた複数の単色エネルギーの粒子線に変換するフィルタである。リッジフィルタ１４から出力された粒子線は、照射対象（ターゲット）に対してＳＯＢＰ（Spread-Out Bragg Peak）を持つ粒子線を形成する。ＳＯＢＰとは粒子線進行距離の所定の範囲で線量が一様な領域であり、例えば、ブラッグ曲線２６（図２（ａ））で表される。図３は、リッジフィルタの例を示す。リッジフィルタ１４は、入射面からの高さｈが異なる、複数の段差のある透過部を持つものが多く利用されている。リッジフィルタ１４の入射面３０からの階段の高さｈに応じて、粒子線は、異なるエネルギーを持つ複数の粒子線に変換される。すなわち、階段の高さｈが小さいほどリッジフィルタでの減衰量が小さいので、深い位置まで到達可能なエネルギーを持つ粒子線が透過される。さらに階段の幅ｗを調整して通過する粒子数を制御することにより、線量レベルを調整することができる。したがってこの階段の幅ｗは、各粒子線に対してかける重みであるといえる。重みをかけることにより粒子線（図２の（ｂ）のブラッグ曲線２４’’）を得ることができる。なお、リッジフィルタ１４の透過部において、各エネルギーの粒子線が出射する部分の形状は任意であるが、本実施の形態では長方形である。リッジフィルタ１４の階段は、最終的に得られたブラッグ曲線２４’’（図２）の間隔（例えば、ピークの間隔）が狭くなれば、滑らかになる。間隔が狭いことは、すなわちエネルギーレベルが近接するブラッグ曲線２４’’（図２）を利用することを意味し、よって隣り合う階段の高さの差ｄが小さくなるためである。リッジフィルタ１４から出力された粒子線１０が、揺動され広がった粒子束１５として示されている。なお、以下の説明では、揺動され広がった粒子束１５についても粒子線１０として言及する。
【００１９】
粒子線１０は、続いて照射ボリューム１６によって適切な範囲に絞られ、腫瘍１３０（図１３）に照射される。照射ボリューム１６により調整されるのは、図１に示す座標系ではＸＹ平面内の照射領域である。ｚ軸方向、すなわち腫瘍１３０（図１３）の深さ方向に照射すべき量は、リッジフィルタ１４によって既に調整されている。すなわち腫瘍１３０（図１３）の厚み方向に一様かつ所望の線量のＳＯＢＰが形成されるよう、粒子線の照射量はリッジフィルタ１４によって既に調整されている。
【００２０】
先に述べたように、リッジフィルタ１４（図１）の入射面からの透過部の高さｈと階段の幅ｗとにより透過する粒子線の線量が決定されるので、それらを調整することにより、ＳＯＢＰを形成することができる。以下、より具体的に説明する。まず、図２の（ａ）は、水ターゲットに入射した粒子線１０（図１）の水の深さに対する線量分布を示す。粒子線の進行距離と線量との関係（線量分布）を表す曲線は、一般に「ブラッグ曲線」と呼ばれる。複数のブラッグ曲線に適当な重みをかけて重ね合わせると、所定の範囲（図２の場合は約１００ｍｍ〜３００ｍｍ）で線量が一様な領域、すなわちＳＯＢＰを持つブラッグ曲線２６を得ることができる。ＳＯＢＰが腫瘍１３０（図１３）の厚みＤに一致するよう、ブラッグ曲線２６により表される粒子線を照射することで、厚みＤの部分の線量を所望の値で一様にできる。図３を参照してこのようなＳＯＢＰを形成するためのリッジフィルタ１４をさらに説明する。リッジフィルタ１４は、図の奥行き方向に同じ断面を有している。粒子線のエネルギーは透過した媒質の長さに応じて減衰するので、リッジフィルタ１４の入射面３０からの透過部の高さｈによって、異なるエネルギーの粒子線１０を複数同時に得ることができる。線量レベルの差は、透過部の各階段の幅ｗに差異があるために生じる。換言すれば、リッジフィルタ１４は、透過部の各階段の幅ｗを変化させることにより粒子の透過面積を変化させ、透過する粒子数を制御する。これによりブラッグ曲線の線量レベル（図２の縦軸方向のレベル）を調整できる。例えば、各階段の幅ｗを２倍にして透過部の面積を２倍にすると、透過する粒子数も２倍になり、ブラッグ曲線は図２の縦軸の正方向に平行移動する。したがってリッジフィルタ１４の階段の幅ｗは、ブラッグ曲線２２、２４を平行移動させる因子である。この因子は、以下、ブラッグ曲線にかける「重み」としても言及される。なお各階段の高さの差ｄの大きさは特に限定されないが、より小さい方が深さ方向（図２の横軸方向）に密接した複数のブラッグ曲線を得ることができ、より容易にかつ均一なＳＯＢＰを得ることができる。
【００２１】
続いて、リッジフィルタの入射面３０（図３）からの階段の高さｈと階段の幅ｗを決定するための手順を説明する。本発明では、ＳＯＢＰを持つブラッグ曲線２６（図２）（以下、「ＳＯＢＰブラッグ曲線２６」という）を得るために行っていた全てのブラッグ曲線についてのモンテカルロ計算を、最大のエネルギーを持つ粒子線について１回だけ行う。すなわち本発明では、モンテカルロ計算により求めるのはブラッグ曲線２２（図２の（ａ））のみである。続いて、ブラッグ曲線２２（図２）を平行移動することによりブラッグ曲線２４’が取得される。ブラッグ曲線２４’が、ＳＯＢＰブラッグ曲線２６（図２）を得るために必要な別のエネルギーのブラッグ曲線となる。平行移動は、図２の横軸方向に行う。このように、すべてのブラッグ曲線の取得に時間のかかるモンテカルロ計算を採用しないことで、リッジフィルタの設計に要する計算時間を短縮することができる。
【００２２】
以下、その手順をより詳しく説明する。図４は、本実施の形態によるリッジフィルタの設計フローを示す。まず初めに、最大エネルギーのブラッグ曲線２２（図２）をモンテカルロ計算により求める（ステップＳ４０）。モンテカルロ計算は周知の計算手法であるのでその説明は省略する。次に、得られたブラッグ曲線２２（図２）をｎ回（ｎ：自然数）平行移動して、ｎ本のブラッグ曲線２４（図２）を得る（ステップＳ４１）。このステップではブラッグ曲線２２（図２）の平行移動をするのみであるから、ｎ本のブラッグ曲線を得るためにモンテカルロ計算を行うよりは短時間ですむ。
【００２３】
次に、各ブラッグ曲線２２、２４’（図２）へかける重みを計算する（ステップＳ４２）。ここで「重み」は、各ブラッグ曲線２２、２４’に（図２）その重みをかけて重ね合わせると、ＳＯＢＰの必要な線量レベルと所望の位置での所望の長さが得られるよう求められる。重みは以下のようにして求める。まず、各ブラッグ曲線２２、２４’（図２）を複数個の点（例えば、各ブラッグ曲線についてｍ個ずつ）でサンプリングする。次に所望のＳＯＢＰを有するＳＯＢＰブラッグ曲線２６（図２）も複数個（ｍ個）の点でサンプリングする。この点は以下、「評価点」として言及される。評価点の位置は任意であるが、各ブラッグ曲線２２、２４’（図２）およびＳＯＢＰブラッグ曲線２６（図２）のｋ番目の評価点は、いずれも横軸方向の同じ位置ｐ_kにおける各ブラッグ曲線２２、２４’（図２）およびＳＯＢＰブラッグ曲線２６（図２）上の点が採用される。ＳＯＢＰブラッグ曲線２６（図２）のＳＯＢＰの位置以外の位置では、線量レベルを低くする方が好ましい。その位置は腫瘍１３０（図１３）ではない正常な組織がある位置であり、そのような組織に照射される粒子線レベルはより低い方がよいからである。なお、所望のＳＯＢＰを得ることだけを考える場合には、ＳＯＢＰブラッグ曲線２６（図２）の所望のＳＯＢＰ位置ですべてのサンプリングを行えばよい。続いて得られたサンプル値を利用して、ブラッグ曲線２２、２４’（図２）のサンプル値を成分に持つｍ×（ｎ＋１）行列Ａと、重みを表す（ｎ＋１）×１重みベクトルｘと、ＳＯＢＰブラッグ曲線２６（図２）上のｍ×１サンプル値列ベクトルｂとを作る。ブラッグ曲線２２および２４’は計（ｎ＋１）本であるので、成分も（ｎ＋１）行、または列だけ必要となる。これらの行列とベクトルの間には、
Ａｘ＝ｂ・・・（１）
の関係が成り立つ。重みベクトルｘの各成分の値は、例えば、逐次代入法または最適化手法を利用して求められる。なお、これらの手法は周知であるので詳細な説明は省略する。ここで、行列Ａの第１列はブラッグ曲線２２（図２）、第２列はブラッグ曲線２２、２４’（図２）中で２番目のエネルギーを持つブラッグ曲線、・・・（ｎ＋１）列目は（ｎ＋１）番目のエネルギーを持つブラッグ曲線を表す。この結果、重みベクトルｘの第ｋ行（ｋ：１〜ｎ＋１の自然数）はブラッグ曲線２２、２４’（図２）のうちｋ番目のエネルギーを持つブラッグ曲線の重みとなる。以上のようにして重みが計算される。
【００２４】
求めた重みとブラッグ曲線２２、２４’（図２）とから、実際にＳＯＢＰを持つブラッグ曲線を計算する。図５は、本発明の手法により得た重みから計算したＳＯＢＰブラッグ曲線のグラフを示す。図２では、横軸は粒子線の進行距離（ｍｍ）、縦軸はＳＯＢＰが得られているか否かの観点から、線量の最大レベルで正規化している。図５によれば、ＳＯＢＰが約１４０ｍｍ〜２９０ｍｍの長さ得られていることがわかる。粒子線のエネルギーは２３０ＭｅＶ、ブラッグ曲線のサンプリング間隔は２．５ｍｍである。なお「粒子線の進行距離」とは、先に「水の深さ」として言及した、粒子線照射装置１（図１）から照射された粒子線が進んだ距離である。
【００２５】
重みが求められると、リッジフィルタの階段の幅ｗを決定することができる（ステップＳ４３）。例えば、あるエネルギーのブラッグ曲線に対する重みが０．２である場合を考える。平行移動して得られたブラッグ曲線が、リッジフィルタの幅ｗ₁のときに実現されるとすると、重みをかけた後のブラッグ曲線は、リッジフィルタの幅０．２×ｗ₁のときに実現される。これは、幅を０．２倍にして面積を０．２倍にすることにより、リッジフィルタを通過する粒子線１０（図３）の粒子数を０．２倍したことを意味する。このようにしてリッジフィルタの階段の幅ｗを決定できる。リッジフィルタの入射面３０（図３）からの階段の高さｈは、ブラッグ曲線２２、２４（図２）の位置に基づいて決定すればよいが、例えば、図２の横軸方向への平行移動の移動量により決定してもよい。この場合は、最大エネルギーのブラッグ曲線２２（図２）を与える入射面３０（図３）からの階段の高さｈ₁が得られれば、他のブラッグ曲線２４（図２）を与える入射面３０（図３）からの階段の高さｈ_kも容易に得ることができる。なお、幅ｗ₁と階段の高さｈ₁は、粒子数および減衰量を考慮して周知の手法により容易に得ることができるので、その説明は省略する。
【００２６】
このようにして得られた階段の幅ｗおよび入射面３０（図３）からの高さｈを用いることによって、リッジフィルタ１４（図３）を設計することができる。この結果、全てのブラッグ曲線についてモンテカルロ計算をしなくとも、ＳＯＢＰを得るためのリッジフィルタを得ることができる。
【００２７】
続いて、以上のようにして求めたＳＯＢＰ計算結果の評価を行い、より一層精度の高いＳＯＢＰを求めるための方法を説明する。本方法では、求めた重みに基づいて設計されたリッジフィルタを、ある粒子線が透過するとしたときのＳＯＢＰの理論値ｂ’を計算する（ステップＳ４４の「はい」に分岐）。すなわち、得られた幅ｗおよび高さｈをもつリッジフィルタにある粒子線を通過させたときの、ＳＯＢＰの理論上の値ｂ’をモンテカルロ計算により計算する（ステップＳ４５）。このモンテカルロ計算は、得られたリッジフィルタと透過させる粒子線のエネルギーとに基づいて行われる計算である。図６は、ＳＯＢＰブラッグ曲線の理論値のグラフを示す。縦軸および横軸等の説明は図５と同じである。この理論値のグラフは、ブラッグ曲線の平行移動により求めたグラフ（図５）と比較して、一様なレベルのＳＯＢＰが得られていないことがわかる。その理由は、エネルギーの低い粒子線はエネルギーの高い粒子線よりも比較的散乱が大きいことが知られているが、ブラッグ曲線の平行移動によって取得された低エネルギー粒子線のブラッグ曲線では粒子線の散乱が考慮されていないからである。これは、実際のブラッグ曲線は、平行移動して得られたブラッグ曲線よりも線量が小さいことを意味する。したがって、より一層精度よく一様なＳＯＢＰを得るためには、エネルギーが低い粒子線の散乱をも考慮して重みを計算し、リッジフィルタを設計しなければならない。よって、上で得られた重みを用いてモンテカルロ計算を行い、ＳＯＢＰの理論値ｂ’と所望のＳＯＢＰ値ｂとの差を補償する必要がある。
【００２８】
再び図４を参照して、得られたＳＯＢＰの理論値ｂ’と、ＳＯＢＰの所望の値ｂとの差を計算し（ステップＳ４６）、その差を最も小さくする（補償する）ための補正重みを計算する（ステップＳ４７）。すなわち「補正重み」は、その補正重みをかけた各ブラッグ曲線２２、２４’（図２）を、前に得られた、重みをかけた各ブラッグ曲線２２、２４’（図２）にさらに重ね合わせて、その差に相当する線量を最も小さくするよう求められる。具体的には、補正重みは以下のようにして求める。前に述べた行列Ａと、補正重みを表す（ｎ＋１）×１補正重みベクトルΔｘと、ＳＯＢＰの理論値ｂ’とＳＯＢＰの所望の値ｂとの差を表す、ｍ×１差分サンプル値列ベクトルΔｂとを用ると、
ＡΔｘ＝Δｂ・・・（２）
が成り立ち、式（１）と同様の周知の計算手法により補正重みベクトルΔｘを求められる。評価点の位置に関しては先の説明と同様に、ｋ番目の評価点は、いずれも同じ位置ｐ_kにおける各ブラッグ曲線およびＳＯＢＰブラッグ曲線上の点が採用される。
【００２９】
式（２）から得られた補正重みΔｘをブラッグ曲線２２、２４’（図２）に加重し、モンテカルロ計算により既に得ている、重みｘをかけたブラッグ曲線にさらに重ね合わせると、その差を補正したＳＯＢＰブラッグ曲線を得ることができる。図７は、先に求めた重みと補正重みとに基づいてモンテカルロ計算で求めたＳＯＢＰブラッグ曲線を示す。補正重みを加える前のＳＯＢＰブラッグ曲線（図６）と比較して、より平坦かつ十分均一なＳＯＢＰブラッグ曲線を得ることができる。したがって、この重みおよび補正重みとをかけたブラッグ曲線２２、２４’それぞれの形状は、ブラッグ曲線２４’’ （図２の（ｂ））の形状にほぼ一致する。このようにして各ブラッグ曲線に対する重みと補正重みとが求められたので、補正重みをさらに考慮してリッジフィルタの階段の幅ｗを決定することができる（図４のステップＳ４３）。
【００３０】
以上、実施の形態１を説明した。本実施の形態によれば、迅速に均一なＳＯＢＰを与えるリッジフィルタを設計することができる。また本実施の形態によれば、理論値と実測値が違う場合であっても、補正重みを用いて実測値と理論値を十分等しくできるので、精度のよいリッジフィルタを設計できる。この場合でも、モンテカルロ計算を２回行うだけで十分な精度のＳＯＢＰを得ることができる。
【００３１】
（実施の形態２）
本実施の形態では、所望のＳＯＢＰレベルｂを得るための他の方法を説明する。実施の形態１では、ＳＯＢＰの理論値ｂ’と所望のＳＯＢＰ値ｂとの差をΔｂとして、Ａ△ｘ＝△ｂから補正重みを決定した。本実施の形態では、所望のＳＯＢＰ値ｂからのＳＯＢＰの理論値ｂ’の差を所望のＳＯＢＰ値ｂに加え、その値を得られるような重みを求める。これは、所望のＳＯＢＰ値ｂを求めようとして生じた差を、あらかじめ所望のＳＯＢＰ値ｂに加えておいて重みを求めることを意味する。このようにして求めた重みによれば、当初から誤差のない所望のＳＯＢＰ値を実現できる。
【００３２】
図８は、本実施の形態によるリッジフィルタの設計フローを示す。ステップＳ８０〜Ｓ８６のそれぞれは、図４のステップＳ４０〜Ｓ４６のそれぞれに対応するので、その説明は省略する。なおステップＳ８６では、所望のＳＯＢＰ値ｂからのＳＯＢＰの理論値ｂ’の差Δｂが計算されている。ステップＳ８７では、この差Δｂを本来の所望のＳＯＢＰ値ｂに加算し、その値を用いて、既に求めた重みｘを修正した修正重みｘ’を求める。すなわち、修正重みｘ’は
Ａｘ’＝ｂ＋Δｂ・・・（３）
を満たすベクトルとして与えられる。これは、所望のＳＯＢＰ値ｂを求めようとして重みｘを求めても、モンテカルロ計算を行って得られた理論値ｂ’は所望のＳＯＢＰ値ｂからΔｂだけ差が生じていたことに基づく。すなわち、所望のＳＯＢＰ値ｂにあらかじめΔｂを加えた値を実現するような修正重みｘ’を求めると、その修正重みｘ’に基づいて計算した理論値は、ｂ＋ΔｂよりもΔｂだけ小さい、もとの所望のＳＯＢＰ値ｂになる。なお、ｘ’を求める手法は、実施の形態１で簡単に述べたような逐次代入法または最適化手法を用いればよい。このようにして得られた修正重みｘ’を新たな重みとすると、この重みは、すでに求めた重みｘに代わるものである。よって、この重みｘ’を各ブラッグ曲線２２、２４’（図２の（ａ））に加重すると、各ブラッグ曲線２２、２４’’（図２の（ｂ））を得ることができ、それらを重ね合わせることによりＳＯＢＰブラッグ曲線（図２（ｂ））を得ることができる。
【００３３】
このようにして各ブラッグ曲線に対する新たな重みが求められたので、その重みに基づいて、既に求めていたリッジフィルタの階段の幅ｗを修正することができる（図８のステップＳ８３）。その手順は実施の形態１で既に説明したので、説明は省略する。
【００３４】
以上、実施の形態２を説明した。本実施の形態によれば、迅速に均一なＳＯＢＰを与えるリッジフィルタを設計することができる。
【００３５】
（実施の形態３）
本実施の形態は、より均一な所望のＳＯＢＰを得るための方法を説明する。すなわち、所望のＳＯＢＰを得る過程で、ＳＯＢＰの右端にピークが生じる場合がある（図５の２８５ｍｍ付近におけるＳＯＢＰブラッグ曲線を参照）。これは、サンプリングした評価位置が連続的でなく離散的であることに起因する。図２の（ｂ）を参照してより具体的に説明する。ブラッグ曲線２２は、評価点として通常採用されるピーク以外の位置、例えば、ピークから少しずれた水の深さ３００ｍｍの位置でもある線量値を持つ。この位置では、他のブラッグ曲線２４’’も何らかの線量値を持つ。このような位置では各ブラッグ曲線にかけられた重みの影響が大きく出ているので、重ね合わされたブラッグ曲線のその位置の線量値の和が所望のＳＯＢＰのレベルよりも大きくなる場合がある。この最大値が突出した形でＳＯＢＰの右端に出現する（シフトする）。これがＳＯＢＰの右端に生じるピークである。これはすなわち複数のブラッグカーブを重ね合わせた結果、ピーク位置が評価点の位置と一致していないことを意味する。計算の都合上、評価点が離散的になることは避けられないが、このような最大値を与える位置が評価点として採用されない場合には、ＳＯＢＰの均一度が悪化する。
【００３６】
そこで、このようなピークが生じたと判断する場合には、新しく生じた突起の位置をＳＯＢＰブラッグ曲線２６上の評価点の１つとして採用し、重みの計算を再度行う。ＳＯＢＰの突起部分だけでなく、別の評価点間の誤差を計算してその誤差が大きい点をさらに探してもよい。この誤差を定量的に判断する方法として、本実施の形態では、ある一定の許容誤差範囲を設け、その範囲を超える場合にはその位置の誤差が大きいと判断する。許容誤差範囲は、例えば所望のＳＯＢＰレベルｂの±０．２％などのように設定できる。したがって一般的には、得られたＳＯＢＰと、所望のＳＯＢＰレベルｂとの誤差を計算し、その誤差が許容誤差範囲を超えている位置を選択して新たな評価位置とし、その位置で再び重みを計算すればよい。それにより、選択された新たな評価位置における誤差を小さくできる。図９は、再度計算した重みに基づいて得られたＳＯＢＰブラッグ曲線を示す。突起が生じていた位置（図５の２８５ｍｍ付近）にはピークはなく、より均一なＳＯＢＰブラッグ曲線が得られていることがわかる。これは、設計したリッジフィルタが高精度であることを意味する。
【００３７】
以上、本発明によるリッジフィルタの設計手順を説明した。続いて、実施の形態１〜３において設計されるリッジフィルタの具体的な例を説明する。以下では、長方形でない、扇形の回転型リッジフィルタを用いても、上述の実施の形態がそのまま適用できることを説明する。
【００３８】
図１０は、回転型リッジフィルタ１００を示す。回転型リッジフィルタ１００は、粒子線の入射方向と平行な軸を中心に回転して、粒子線が回転する透過部１０２を透過する際にエネルギーを変化させ、複数の異なるエネルギーを持つ粒子線に変換するフィルタである。回転型リッジフィルタ１００は、複数の扇形の透過部１０２がねじ等の固定具１０４で固定されている。複数の扇形の透過部１０２は積み重ねられ、階段状の透過部を形成している。この構成は、これまでの説明のリッジフィルタ透過部と同様である。回転型リッジフィルタ１００では、透過部１０２の交換により弧の幅ｗ（階段の幅ｗ）を変えることができる。上で述べたように階段の幅ｗを変化させると粒子の透過面積が変化するので、階段の幅ｗを変化させることは、所定のエネルギーを持つブラッグ曲線にかける重みを自由に変化させられることをあらわす。すなわち、回転型リッジフィルタも実施の形態１〜３で求めた重みを反映して設計できる。例えば、あるエネルギーのブラッグ曲線に対する重みが０．２である場合を考える。透過する粒子数を０．２倍にするためには、扇形の弧の幅ｗを０．２倍にすればよい。
【００３９】
なお、階段の幅ｗは、透過部１０２の交換のみならず回転方向にずらしてもよい。例えば回転方向に透過部１０２の幅を狭くするようずらすと、回転方向と逆側の透過部１０２は幅が広くなる。したがってこの場合には、粒子線は、リッジフィルタ１００の回転方向に幅を狭く調整した領域にのみ粒子線が透過するよう照射する必要がある。
【００４０】
続いて、回転型リッジフィルタ１００（図１０）とは別の回転型リッジフィルタについて説明する。上で説明したように、この回転型リッジフィルタも上述の実施の形態はそのまま適用できる。
【００４１】
図１１の（ａ）は、回転型リッジフィルタ１１０を示す。回転型リッジフィルタ１１０は、回転型リッジフィルタ１００（図１０）とは異なるのは、４つの透過部１１２、１１４、１１６、１１８から構成されている点である。各透過部内の階段の高さは一定であるが、４種類の階段の高さは異なる。プロペラが回転することによって、４枚の透過部１１２、１１４、１１６、１１８を持つリッジフィルタ、すなわち４種類の高さを持つリッジフィルタを等価的に実現できる。リッジフィルタは切削で製作するので、各透過部の製作段階では階段を切削する工程を減らすことができ、階段の高さが高いとより容易に製作できる。図１２は、４つの透過部をもつ回転型リッジフィルタ１２０のより具体的な外観を示す。直線の矢印が指す方向が粒子線が入射する方向である。回転型リッジフィルタ１２０は、粒子線の入射方向に垂直な平面内で回転する。
【００４２】
図１１の（ｂ）は、透過部１１２および１１４の断面図を示す。粒子線の入射面からの高さｈは、透過部１１２と透過部１１４とでは異なることがわかる。より一般的には、回転型リッジフィルタがそれぞれ高さが異なるｎ種類の透過部を備えていれば、等価的にｎ段の階段を持つリッジフィルタとして扱うことができるので、上述した効果を得ることができる。なお各階段の高さの差ｄはどの透過部で同じであってもよいし、必要に応じて異ならせてもよい。
【００４３】
【発明の効果】
本発明によれば、迅速にＳＯＢＰを与えるリッジフィルタを設計することができる。さらに本発明によれば、理論値と実測値が違う場合であっても、補正重みを用いて、またはあらかじめその差を考慮して実測値を理論値に近づけることができるので、精度のよいリッジフィルタを設計できる。
【００４４】
さらに本発明によれば、均一なＳＯＢＰを得ることができるので、精度のよいリッジフィルタを設計することができる。
【図面の簡単な説明】
【図１】リッジフィルタを用いたワブラー法による粒子線照射装置１の模式図である。
【図２】水ターゲットに入射した粒子線１０（図１）の水の深さに対する線量分布を示す図である。（ａ）はモンテカルロ計算によるブラッグ曲線２２および２４と、本発明による平行移動したブラッグ曲線２４’と、ＳＯＢＰブラッグ曲線２６とを示すグラフである。（ｂ）は重みをかけたブラッグ曲線２２および２４’’と、ＳＯＢＰブラッグ曲線２６とを示すグラフである。
【図３】リッジフィルタの例を示す図である。
【図４】実施の形態１によるリッジフィルタの設計フローを示す。
【図５】本発明の手法により得た重みから計算したＳＯＢＰブラッグ曲線のグラフである。
【図６】ＳＯＢＰブラッグ曲線の理論値のグラフである。
【図７】重みと補正重みとに基づいてモンテカルロ計算したＳＯＢＰブラッグ曲線を示す図である。
【図８】実施の形態２によるリッジフィルタの設計フローを示す図である。
【図９】再度計算した重みに基づいて得られたＳＯＢＰブラッグ曲線を示す図である。
【図１０】回転型リッジフィルタを示す図である。
【図１１】（ａ）は、回転型リッジフィルタ１１０を示す図である。（ｂ）は、透過部１１２および１１４の断面図を示す図である。
【図１２】４つの透過部をもつ回転型リッジフィルタのより具体的な外観を示す。
【図１３】粒子線を照射される腫瘍の例を示す図である。
【図１４】従来のリッジフィルタの設計フローを示す図である。
【符号の説明】
１粒子線照射装置、１０粒子線、１１１対の水平方向揺動用電磁石、１２１対の垂直方向揺動用電磁石、１３レンジシフタ、１４リッジフィルタ、１６照射ボリューム[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a method for designing a ridge filter that creates a predetermined dose distribution in a target by changing the dose of a penetrating particle beam.
[0002]
[Prior art]
One of the methods for treating cancer is a treatment method in which a particle beam is irradiated to a cancer tumor. In this treatment method, cancer cells are killed by irradiating a cancer tumor with a particle beam. FIG. 13 shows an example of a tumor 130 irradiated with the particle beam 10. The particle beam 10 is irradiated as follows. Since the tumor 130 has a certain extent on a plane perpendicular to the particle beam 10, the particle beam 10 is shaped according to the extent by the collimator 132 and the correction medium 134. This is to prevent the normal beam other than the tumor 130 from being irradiated with the particle beam 10. Since the tumor 130 has a thickness D also in the depth direction, it is necessary to irradiate the portion of the thickness D with the uniform particle beam 10 with a desired dose.
[0003]
The particle beam 10 having a uniform dose in a predetermined range in the depth direction is obtained from a plurality of particle beams having different energies. For example, FIG. 2B shows the dose distribution (Bragg curve) of the particle beam incident on the water target. The region SOBP (Spread-Out Bragg Peak) with a uniform dose is represented by the Bragg curves 22, 24 ′. Obtained by superimposing each of '. Each of the Bragg curves 22, 24 ″ is obtained by applying a weight to the Bragg curves 22, 24 in FIG.
[0004]
A plurality of particle beams having different energies for obtaining such SOBP are generated using a ridge filter. The ridge filter is a filter that converts an incident particle beam having a certain energy into a plurality of monochromatic energy particle beams each weighted. That is, the ridge filter receives the energy particle beam represented by the Bragg curve 22 and converts it into a plurality of monochromatic energy particle beams represented by the Bragg curve 24 ″. By superimposing these particle beams, the particle beam output from the ridge filter forms a particle beam having SOBP represented by the Bragg curve 26 with respect to the irradiation target (target). FIG. 3 shows an example of a ridge filter. Depending on the height h of the staircase from the incident surface 30 of the ridge filter, the particle beam is converted into a plurality of particle beams having different energies. That is, since the amount of attenuation in the ridge filter is smaller as the step height h is smaller, a particle beam having energy that can reach a deep position is transmitted. Therefore, the height h of the stairs defines the position of the Bragg curve in the horizontal axis direction in the graph of FIG. Furthermore, the dose level can be adjusted by adjusting the width w of the stairs to control the number of particles passing through. Therefore, the width w of the stairs corresponds to the weight of the Bragg curve in the vertical axis direction in the graph of FIG. In other words, the width w of this step corresponds to the particle bundle.
[0005]
In order to obtain a uniform SOBP at a desired distance with respect to the traveling distance of the particle beam and at a desired dose value, it is necessary to design a ridge filter having an appropriate step height h and width w. A conventional ridge filter design procedure is described, for example, in “MONTE CARLO SIMULATION OF A PROTON THERAPY SYSTEM FOR THE CALCULATION OF THE DOSE DISTRIBUTION IN A PATIENT” TERA, 95, No. 4, pp25 (1995). FIG. 14 shows a design flow of a conventional ridge filter. (1) First, the target SOBP width is determined. (2) A plurality of Bragg curves (Bragg curves 24 in FIG. 2) having different energies are obtained by Monte Carlo calculation (step S142). (3) The weight x of all the Bragg curves is determined from the expression Ax = b for the target evaluation point on the SOBP and the target value at that position (step S144). Here, the matrix A stores sample values on each Bragg curve 24. The first row is the maximum energy value, the second row is the second energy value, etc., and b is the target value. (4) The step width w of the ridge filter is determined based on the obtained weight x (step S146). At this time, the difference in the height of the stairs (the height of each step of the stairs) is constant.
[0006]
[Problems to be solved by the invention]
As described in (2) above, conventionally, a plurality of Monte Carlo calculations had to be performed to obtain a plurality of Bragg curves 24 (FIG. 2A) having different energies. In general, Monte Carlo calculation takes time, and this makes it impossible to design quickly.
[0007]
Furthermore, the conventional method has a problem that an SOBP with good flatness cannot be obtained depending on the target position of the evaluation point on the SOBP. For example, since the evaluation point adopts a jump position, the behavior of the SOBP value between the evaluation points is not taken into consideration, and therefore the difference from the target value may be large at points other than the evaluation point. In particular, when the peak of the Bragg curve with the maximum energy is not adopted as the evaluation point, the influence of the weight is large, and a protrusion may appear at the right end of the SOBP graph (a value around 285 mm in the horizontal axis in FIG. 5 described later).
[0008]
When the theoretical value was calculated using the ridge filter once created, if the theoretical value did not match the actually required value, there was no appropriate compensation method, and the design had to be performed again. Alternatively, it is necessary to obtain a plurality of Bragg curves based on actually measured values and design a ridge filter based on the obtained Bragg curves.
[0009]
An object of the present invention is to provide a method for quickly designing a ridge filter that achieves uniform SOBP.
[0010]
[Means for Solving the Problems]
The ridge filter design method of the present invention adjusts the height and width of the transmission part through which the particle beam passes using a Bragg curve representing the relationship between the traveling distance of the particle beam and the dose, and the dose after transmission is A method of designing a ridge filter that realizes uniform SOBP with a desired value within a predetermined range of the travel distance, the first Bragg curve of the particle beam being obtained by Monte Carlo calculation, and the first N times (n: natural number) of the Bragg curve are translated to obtain n second Bragg curves, the desired value of the SOBP, the first Bragg curve, and the first Calculating a weight to be applied to each of the first Bragg curve and the second Bragg curve based on two Bragg curves, the weighted first Bragg curve And calculating each of the second Bragg curves to calculate a weight so as to obtain the SOBP having the desired value, and determining a width of the transmission portion based on the weight. Thus, the above object is achieved.
[0011]
The ridge filter design method performs a Monte Carlo calculation based on the weight, calculates a theoretical value of the SOBP, calculates a difference between the theoretical value and the desired value, and the difference; Calculating a correction weight to be applied to each of the first Bragg curve and the second Bragg curve for correcting the difference based on the first Bragg curve and the second Bragg curve; Each of the first Bragg curve and the second Bragg curve to which the correction weight is applied, and each of the first Bragg curve and the second Bragg curve to which the weight is further applied. Calculating the correction weight so that the dose corresponding to the difference becomes the smallest when superimposed, and the determining step further adds the correction weight to the correction weight. It may determine the width of the transmissive portion Zui.
[0012]
The ridge filter design method includes performing a Monte Carlo calculation based on the weight, calculating a theoretical value of the SOBP, calculating a difference between the desired value and the theoretical value, and the desired value. Calculating a correction weight obtained by correcting the weight based on the value obtained by adding the difference to the first Bragg curve and the second Bragg curve, and applying the correction weight Calculating a weight so as to obtain the SOBP of the desired value upon superimposing the first Bragg curve and each of the second Bragg curves, wherein the determining step comprises the correction You may correct the width | variety of the said permeation | transmission part based on a weight.
[0013]
The ridge filter design method includes a step of calculating an error between the SOBP value obtained based on the weight and the correction weight and the desired value of the SOBP, and the error exceeds a predetermined allowable range. And a step of calculating the weight and the step of calculating the correction weight so as to reduce an error in the selected position.
[0014]
The ridge filter may be a ridge filter including the transmission portion having a height parallel to the incident direction of the particle beam and a width perpendicular to the incident direction of the particle beam.
[0015]
The ridge filter is a rotary ridge filter having at least one transmission part that rotates about an axis parallel to the incident direction, and the transmission part has a width in the rotation direction of the transmission part. It may be adjusted by changing.
[0016]
The rotary ridge filter may include a plurality of transmission parts, and each of the plurality of transmission parts may have a different height.
[0017]
DETAILED DESCRIPTION OF THE INVENTION
(Embodiment 1)
In the following, after describing an apparatus using a ridge filter designed according to the present invention, a procedure for designing the ridge filter will be described. In this specification, the particle beam is, for example, a proton beam. Moreover, the object irradiated with the particle beam is a tumor 130 (FIG. 13) having a thickness D in the depth direction. The “depth direction” means a direction in which the particle beam is irradiated and travels.
[0018]
FIG. 1 is a schematic diagram of a particle beam irradiation apparatus 1 by a wobbler method using a ridge filter (“Instrumentation for treatment of cancer using proton and light-iron beams” Rev. Sci. Instrum. Vol. 64 No. 8 pp. .2055-2091 (1993)). The particle beam irradiation apparatus 1 adjusts the region (including the depth direction) to which the particle beam 10 is irradiated according to the shape of the tumor 130 (FIG. 13), and the particle beam 10 having a necessary dose is transferred to the tumor 130 (FIG. 13). It is an apparatus that irradiates. The particle beam irradiation apparatus 1 includes a pair of horizontally oscillating electromagnets 11, a pair of vertically oscillating electromagnets 12, a range shifter 13 for changing the energy of particles, a ridge filter 14, and an irradiation volume 16. Including. The particle beam 10 radiated from a particle beam source (not shown) is swung by using a pair of horizontal swinging electromagnets 11 and a pair of vertical swinging electromagnets 12 to expand the irradiation area. . The oscillated particle beam 10 passes through the range shifter 13 and the energy of the particle beam 10 is reduced by an appropriate amount. Thereafter, the particle beam 10 enters the ridge filter 14. The ridge filter 14 is a filter that converts the incident particle beam into a plurality of monochromatic energy particle beams each weighted by utilizing the attenuation of the particle beam energy according to the length of the transmitted medium. . The particle beam output from the ridge filter 14 forms a particle beam having SOBP (Spread-Out Bragg Peak) with respect to the irradiation target (target). The SOBP is a region where the dose is uniform within a predetermined range of the particle beam traveling distance, and is represented by a Bragg curve 26 (FIG. 2A), for example. FIG. 3 shows an example of a ridge filter. As the ridge filter 14, a filter having a plurality of steps with different heights h from the incident surface is often used. Depending on the height h of the staircase from the incident surface 30 of the ridge filter 14, the particle beam is converted into a plurality of particle beams having different energies. That is, the smaller the staircase height h, the smaller the amount of attenuation in the ridge filter, so that a particle beam having energy that can reach a deep position is transmitted. Furthermore, the dose level can be adjusted by adjusting the width w of the stairs to control the number of particles passing through. Therefore, it can be said that the width w of the stairs is a weight applied to each particle beam. By applying a weight, a particle beam (Bragg curve 24 ″ in FIG. 2B) can be obtained. In the transmission part of the ridge filter 14, the shape of the part where the particle beam of each energy is emitted is arbitrary, but in the present embodiment, it is a rectangle. The staircase of the ridge filter 14 becomes smooth if the interval (for example, peak interval) of the Bragg curve 24 ″ (FIG. 2) finally obtained becomes narrow. The narrow distance means that the Bragg curve 24 ″ (FIG. 2) having close energy levels is used, and thus the difference d between the heights of adjacent stairs becomes small. The particle beam 10 output from the ridge filter 14 is shown as a particle bundle 15 that is swung and spread. In the following description, the particle bundle 15 that is swung and spread is also referred to as the particle beam 10.
[0019]
The particle beam 10 is subsequently focused to an appropriate range by the irradiation volume 16 and irradiated to the tumor 130 (FIG. 13). What is adjusted by the irradiation volume 16 is an irradiation area in the XY plane in the coordinate system shown in FIG. The amount to be irradiated in the z-axis direction, that is, the depth direction of the tumor 130 (FIG. 13) has already been adjusted by the ridge filter 14. That is, the dose of the particle beam has already been adjusted by the ridge filter 14 so that a uniform and desired dose of SOBP is formed in the thickness direction of the tumor 130 (FIG. 13).
[0020]
As described above, since the dose of the particle beam to be transmitted is determined by the height h of the transmission part from the incident surface of the ridge filter 14 (FIG. 1) and the width w of the staircase, by adjusting them, , SOBP can be formed. More specific description will be given below. First, (a) of FIG. 2 shows the dose distribution with respect to the depth of water of the particle beam 10 (FIG. 1) incident on the water target. A curve representing the relationship (dose distribution) between the traveling distance of a particle beam and a dose (dose distribution) is generally called a “Bragg curve”. By superimposing a plurality of Bragg curves with appropriate weights, a Bragg curve 26 having a uniform dose in a predetermined range (about 100 mm to 300 mm in the case of FIG. 2), that is, a SOBP can be obtained. By irradiating the particle beam represented by the Bragg curve 26 so that the SOBP matches the thickness D of the tumor 130 (FIG. 13), the dose of the thickness D portion can be made uniform at a desired value. The ridge filter 14 for forming such SOBP will be further described with reference to FIG. The ridge filter 14 has the same cross section in the depth direction of the drawing. Since the energy of the particle beam is attenuated according to the length of the transmitted medium, a plurality of particle beams 10 having different energies can be obtained simultaneously depending on the height h of the transmission part from the incident surface 30 of the ridge filter 14. The difference in dose level occurs because there is a difference in the width w of each step of the transmission part. In other words, the ridge filter 14 changes the transmission area of particles by changing the width w of each step of the transmission part, and controls the number of transmitted particles. Thereby, the dose level (level in the vertical axis direction of FIG. 2) of the Bragg curve can be adjusted. For example, if the width w of each staircase is doubled and the area of the transmission part is doubled, the number of transmitted particles is also doubled, and the Bragg curve is translated in the positive direction of the vertical axis in FIG. Therefore, the step width w of the ridge filter 14 is a factor for translating the Bragg curves 22 and 24. This factor is also referred to below as the “weight” applied to the Bragg curve. Although the magnitude of the height difference d between the steps is not particularly limited, a smaller one can obtain a plurality of Bragg curves closely in the depth direction (horizontal axis direction in FIG. 2). Uniform SOBP can be obtained.
[0021]
Next, a procedure for determining the height h of the staircase from the incident surface 30 (FIG. 3) and the width w of the staircase will be described. In the present invention, the Monte Carlo calculation for all the Bragg curves performed to obtain the Bragg curve 26 having SOBP (FIG. 2) (hereinafter referred to as “SOBP Bragg curve 26”) is performed for the particle beam having the maximum energy. Perform only once. That is, in the present invention, only the Bragg curve 22 ((a) of FIG. 2) is obtained by Monte Carlo calculation. Subsequently, a Bragg curve 24 ′ is obtained by translating the Bragg curve 22 (FIG. 2). The Bragg curve 24 'becomes another Bragg curve of energy necessary to obtain the SOBP Bragg curve 26 (FIG. 2). The parallel movement is performed in the horizontal axis direction of FIG. In this way, the calculation time required for designing the ridge filter can be shortened by not adopting the Monte Carlo calculation which takes time to acquire all the Bragg curves.
[0022]
Hereinafter, the procedure will be described in more detail. FIG. 4 shows a design flow of the ridge filter according to the present embodiment. First, the Bragg curve 22 (FIG. 2) with the maximum energy is obtained by Monte Carlo calculation (step S40). Since the Monte Carlo calculation is a well-known calculation method, its description is omitted. Next, the obtained Bragg curve 22 (FIG. 2) is translated n times (n: natural number) to obtain n Bragg curves 24 (FIG. 2) (step S41). In this step, since only the Bragg curve 22 (FIG. 2) is translated, it takes less time to perform Monte Carlo calculation in order to obtain n Bragg curves.
[0023]
Next, the weight applied to each Bragg curve 22, 24 '(FIG. 2) is calculated (step S42). Here, the “weight” is obtained so that the required dose level of SOBP and the desired length at the desired position can be obtained by superimposing each weight on the Bragg curves 22 and 24 ′ (FIG. 2). . The weight is obtained as follows. First, each Bragg curve 22, 24 ′ (FIG. 2) is sampled at a plurality of points (eg, m for each Bragg curve). Next, the SOBP Bragg curve 26 (FIG. 2) having a desired SOBP is also sampled at a plurality (m) of points. This point is hereinafter referred to as “evaluation point”. The position of the evaluation point is arbitrary, but the k-th evaluation point of each Bragg curve 22, 24 ′ (FIG. 2) and the SOBP Bragg curve 26 (FIG. 2) is the same position p in the horizontal axis direction. _k The points on each Bragg curve 22, 24 ′ (FIG. 2) and SOBP Bragg curve 26 (FIG. 2) are taken. At a position other than the SOBP position on the SOBP Bragg curve 26 (FIG. 2), it is preferable to lower the dose level. This is because the position is a position where there is a normal tissue that is not the tumor 130 (FIG. 13), and it is better that the particle beam level irradiated to such tissue is lower. When only considering obtaining a desired SOBP, all sampling may be performed at a desired SOBP position on the SOBP Bragg curve 26 (FIG. 2). Subsequently, by using the obtained sample values, an m × (n + 1) matrix A having the sample values of the Bragg curves 22, 24 ′ (FIG. 2) as components, and an (n + 1) × 1 weight vector x representing weights, , And an m × 1 sample value sequence vector b on the SOBP Bragg curve 26 (FIG. 2). Since there are a total of (n + 1) Bragg curves 22 and 24 ', only (n + 1) rows or columns are required for the components. Between these matrices and vectors,
Ax = b (1)
The relationship holds. The value of each component of the weight vector x is obtained using, for example, a sequential substitution method or an optimization method. Since these methods are well known, detailed description thereof is omitted. Here, the first column of the matrix A is the Bragg curve 22 (FIG. 2), the second column is the Bragg curve having the second energy in the Bragg curves 22 and 24 ′ (FIG. 2), (n + 1) columns. The eye represents a Bragg curve having the (n + 1) th energy. As a result, the k-th row (k: natural number from 1 to n + 1) of the weight vector x is the weight of the Bragg curve having the k-th energy among the Bragg curves 22, 24 ′ (FIG. 2). The weight is calculated as described above.
[0024]
A Bragg curve having an actual SOBP is calculated from the obtained weight and the Bragg curves 22, 24 '(FIG. 2). FIG. 5 shows a graph of the SOBP Bragg curve calculated from the weights obtained by the method of the present invention. In FIG. 2, the horizontal axis is normalized by the maximum dose level from the viewpoint of whether or not the traveling distance (mm) of the particle beam is obtained, and the vertical axis is whether SOBP is obtained. According to FIG. 5, it can be seen that the SOBP has a length of about 140 mm to 290 mm. The particle beam energy is 230 MeV, and the Bragg curve sampling interval is 2.5 mm. The “traveling distance of the particle beam” is a distance traveled by the particle beam irradiated from the particle beam irradiation apparatus 1 (FIG. 1), which has been referred to as “depth of water”.
[0025]
When the weight is obtained, the step width w of the ridge filter can be determined (step S43). For example, consider a case where the weight for a Bragg curve of a certain energy is 0.2. The Bragg curve obtained by translation is the width w of the ridge filter. ₁ Is realized, the Bragg curve after weighting has a width of 0.2 × w of the ridge filter. ₁ Realized when. This means that the number of particles of the particle beam 10 (FIG. 3) passing through the ridge filter is increased by 0.2 by increasing the width by 0.2 and the area by 0.2. In this way, the step width w of the ridge filter can be determined. The height h of the staircase from the incident surface 30 (FIG. 3) of the ridge filter may be determined based on the positions of the Bragg curves 22 and 24 (FIG. 2). For example, it is parallel to the horizontal axis direction of FIG. You may determine by the movement amount of a movement. In this case, the height h of the staircase from the incident surface 30 (FIG. 3) that gives the Bragg curve 22 (FIG. 2) of maximum energy. ₁ Is obtained, the height h of the staircase from the incident surface 30 (FIG. 3) giving another Bragg curve 24 (FIG. 2). _k Can also be easily obtained. Width w ₁ And the height of the stairs h ₁ Can be easily obtained by a known method in consideration of the number of particles and the amount of attenuation, and the description thereof will be omitted.
[0026]
The ridge filter 14 (FIG. 3) can be designed by using the step width w and the height h from the incident surface 30 (FIG. 3) thus obtained. As a result, a ridge filter for obtaining SOBP can be obtained without performing Monte Carlo calculation for all the Bragg curves.
[0027]
Subsequently, a method for evaluating the SOBP calculation result obtained as described above and obtaining SOBP with higher accuracy will be described. In this method, the SOBP theoretical value b ′ when a particle beam passes through the ridge filter designed based on the obtained weight is calculated (branch to “Yes” in step S44). That is, the theoretical value b ′ of SOBP when passing through the particle beam in the obtained ridge filter having width w and height h is calculated by Monte Carlo calculation (step S45). This Monte Carlo calculation is performed based on the obtained ridge filter and the energy of the particle beam to be transmitted. FIG. 6 shows a graph of the theoretical values of the SOBP Bragg curve. Description of the vertical and horizontal axes is the same as in FIG. It can be seen that the graph of this theoretical value does not provide a uniform level of SOBP, compared to the graph obtained by the parallel movement of the Bragg curve (FIG. 5). The reason is that a low energy particle beam is known to scatter more than a high energy particle beam, but the low energy particle beam obtained by the translation of the Bragg curve is This is because scattering is not considered. This means that the actual Bragg curve has a smaller dose than the Bragg curve obtained by translation. Therefore, in order to obtain a uniform SOBP with higher accuracy, the ridge filter must be designed by calculating weights in consideration of scattering of particle beams having low energy. Therefore, it is necessary to perform a Monte Carlo calculation using the weights obtained above to compensate for the difference between the SOBP theoretical value b ′ and the desired SOBP value b.
[0028]
Referring to FIG. 4 again, the difference between the obtained SOBP theoretical value b ′ and the desired SOBP value b is calculated (step S46), and the correction weight for minimizing (compensating) the difference is calculated. Is calculated (step S47). That is, the “correction weight” is obtained by superimposing the Bragg curves 22, 24 ′ (FIG. 2) to which the correction weight is applied on the weighted Bragg curves 22, 24 ′ (FIG. 2) obtained previously. In addition, it is required to minimize the dose corresponding to the difference. Specifically, the correction weight is obtained as follows. M × 1 difference sample value sequence vector representing the difference between the matrix A described above, the (n + 1) × 1 correction weight vector Δx representing the correction weight, and the SOBP theoretical value b ′ and the desired SOBP value b Using Δb,
AΔx = Δb (2)
Thus, the correction weight vector Δx can be obtained by a well-known calculation method similar to the equation (1). As for the position of the evaluation point, the k-th evaluation point is the same position p as in the previous description. _k Each Bragg curve and a point on the SOBP Bragg curve are adopted.
[0029]
When the correction weight Δx obtained from Equation (2) is weighted to the Bragg curves 22, 24 ′ (FIG. 2) and further superimposed on the Bragg curve multiplied by the weight x already obtained by Monte Carlo calculation, the difference is calculated. A corrected SOBP Bragg curve can be obtained. FIG. 7 shows an SOBP Bragg curve obtained by Monte Carlo calculation based on the previously obtained weight and the correction weight. Compared to the SOBP Bragg curve before applying the correction weight (FIG. 6), a flatter and sufficiently uniform SOBP Bragg curve can be obtained. Therefore, the shape of each of the Bragg curves 22 and 24 ′ multiplied by the weight and the correction weight substantially matches the shape of the Bragg curve 24 ″ (FIG. 2B). Since the weight and the correction weight for each Bragg curve are obtained in this way, the step width w of the ridge filter can be determined by further considering the correction weight (step S43 in FIG. 4).
[0030]
The first embodiment has been described above. According to the present embodiment, it is possible to design a ridge filter that provides uniform SOBP quickly. Further, according to the present embodiment, even if the theoretical value and the actually measured value are different, the actually measured value and the theoretical value can be made sufficiently equal using the correction weight, so that a highly accurate ridge filter can be designed. Even in this case, a sufficiently accurate SOBP can be obtained by performing the Monte Carlo calculation only twice.
[0031]
(Embodiment 2)
In the present embodiment, another method for obtaining a desired SOBP level b will be described. In the first embodiment, the difference between the SOBP theoretical value b ′ and the desired SOBP value b is Δb, and the correction weight is determined from AΔx = Δb. In the present embodiment, a difference between the SOBP theoretical value b ′ from the desired SOBP value b is added to the desired SOBP value b, and a weight capable of obtaining the value is obtained. This means that the difference generated when the desired SOBP value b is obtained is added to the desired SOBP value b in advance to obtain the weight. According to the weights thus obtained, a desired SOBP value with no error can be realized from the beginning.
[0032]
FIG. 8 shows a design flow of the ridge filter according to the present embodiment. Each of steps S80 to S86 corresponds to each of steps S40 to S46 in FIG. In step S86, a difference Δb between the SOBP theoretical value b ′ and the desired SOBP value b is calculated. In step S87, the difference Δb is added to the original desired SOBP value b, and a corrected weight x ′ obtained by correcting the already determined weight x is obtained using the value. That is, the correction weight x ′ is
Ax ′ = b + Δb (3)
Given as a vector that satisfies. This is based on the fact that the theoretical value b ′ obtained by performing the Monte Carlo calculation is different from the desired SOBP value b by Δb even if the weight x is obtained in order to obtain the desired SOBP value b. That is, when a correction weight x ′ that realizes a value obtained by adding Δb in advance to a desired SOBP value b is obtained, a theoretical value calculated based on the correction weight x ′ is smaller than b + Δb by Δb. Is the desired SOBP value b. As a method for obtaining x ′, a sequential substitution method or an optimization method as briefly described in the first embodiment may be used. Assuming that the modified weight x ′ obtained in this way is a new weight, this weight replaces the already obtained weight x. Therefore, when this weight x ′ is weighted to each Bragg curve 22, 24 ′ (FIG. 2A), each Bragg curve 22, 24 ″ (FIG. 2B) can be obtained. By superimposing, an SOBP Bragg curve (FIG. 2B) can be obtained.
[0033]
Thus, since the new weight with respect to each Bragg curve was calculated | required, the width | variety w of the staircase of the ridge filter already calculated | required can be corrected based on the weight (step S83 of FIG. 8). Since the procedure has already been described in Embodiment 1, the description thereof is omitted.
[0034]
The second embodiment has been described above. According to the present embodiment, it is possible to design a ridge filter that provides uniform SOBP quickly.
[0035]
(Embodiment 3)
This embodiment describes a method for obtaining a more uniform desired SOBP. That is, in the process of obtaining a desired SOBP, a peak may occur at the right end of the SOBP (see the SOBP Bragg curve near 285 mm in FIG. 5). This is because the sampled evaluation positions are not continuous but discrete. A more specific description will be given with reference to FIG. The Bragg curve 22 has a dose value that is also a position other than the peak that is normally adopted as an evaluation point, for example, a position at a depth of 300 mm of water slightly deviated from the peak. In this position, the other Bragg curve 24 '' also has some dose value. Since the influence of the weight applied to each Bragg curve is significant at such a position, the sum of the dose values at that position of the superimposed Bragg curve may be greater than the desired SOBP level. This maximum value appears (shifts) at the right end of the SOBP in a protruding form. This is the peak that occurs at the right end of SOBP. This means that as a result of superimposing a plurality of Bragg curves, the peak position does not coincide with the position of the evaluation point. For convenience of calculation, it is inevitable that the evaluation points become discrete. However, if such a position that gives the maximum value is not adopted as the evaluation point, the uniformity of SOBP deteriorates.
[0036]
Therefore, when it is determined that such a peak has occurred, the position of the newly generated protrusion is adopted as one of the evaluation points on the SOBP Bragg curve 26, and the weight is calculated again. In addition to the SOBP protrusion, an error between other evaluation points may be calculated to further search for a point having a large error. As a method for quantitatively determining this error, in the present embodiment, a certain allowable error range is provided, and when the error exceeds this range, it is determined that the error in the position is large. The allowable error range can be set, for example, as ± 0.2% of the desired SOBP level b. Therefore, in general, an error between the obtained SOBP and the desired SOBP level b is calculated, a position where the error exceeds the allowable error range is selected as a new evaluation position, and the weight is again set at that position. Should be calculated. Thereby, the error in the selected new evaluation position can be reduced. FIG. 9 shows the SOBP Bragg curve obtained based on the recalculated weight. It can be seen that there is no peak at the position where the protrusion occurred (near 285 mm in FIG. 5), and a more uniform SOBP Bragg curve is obtained. This means that the designed ridge filter is highly accurate.
[0037]
The design procedure of the ridge filter according to the present invention has been described above. Subsequently, a specific example of the ridge filter designed in the first to third embodiments will be described. Hereinafter, it will be described that the above-described embodiment can be applied as it is even when a fan-shaped rotary ridge filter that is not rectangular is used.
[0038]
FIG. 10 shows the rotary ridge filter 100. The rotary ridge filter 100 rotates about an axis parallel to the incident direction of the particle beam, changes the energy when the particle beam passes through the rotating transmission part 102, and converts the particle beam into a plurality of particle beams having different energies. This is the filter to convert. In the rotary ridge filter 100, a plurality of fan-shaped transmission portions 102 are fixed by a fixture 104 such as a screw. The plurality of fan-shaped transmission parts 102 are stacked to form a step-shaped transmission part. This configuration is the same as the ridge filter transmission portion described so far. In the rotary ridge filter 100, the arc width w (step width w) can be changed by exchanging the transmission part 102. As described above, changing the step width w changes the particle transmission area, so changing the step width w can freely change the weight applied to the Bragg curve with a given energy. Is expressed. That is, the rotary ridge filter can also be designed reflecting the weights obtained in the first to third embodiments. For example, consider a case where the weight for a Bragg curve of a certain energy is 0.2. In order to increase the number of transmitted particles by 0.2, the width w of the fan-shaped arc may be increased by 0.2.
[0039]
Note that the width w of the stairs may be shifted in the rotation direction as well as the replacement of the transmission part 102. For example, when the width of the transmission part 102 is shifted in the rotation direction so as to be narrowed, the transmission part 102 on the opposite side to the rotation direction becomes wider. Therefore, in this case, it is necessary to irradiate the particle beam so that the particle beam is transmitted only to the region whose width is adjusted to be narrow in the rotation direction of the ridge filter 100.
[0040]
Next, a rotary ridge filter different from the rotary ridge filter 100 (FIG. 10) will be described. As described above, the above-described embodiment can be applied to this rotary ridge filter as it is.
[0041]
FIG. 11A shows the rotary ridge filter 110. The rotary ridge filter 110 is different from the rotary ridge filter 100 (FIG. 10) in that it includes four transmission portions 112, 114, 116, and 118. The height of the stairs in each transmission part is constant, but the heights of the four types of stairs are different. By rotating the propeller, a ridge filter having four transmission portions 112, 114, 116, and 118, that is, a ridge filter having four kinds of heights can be equivalently realized. Since the ridge filter is manufactured by cutting, the steps for cutting the stairs can be reduced at the manufacturing stage of each transmission part, and the ridge filter can be manufactured more easily when the height of the stairs is high. FIG. 12 shows a more specific appearance of the rotary ridge filter 120 having four transmission parts. The direction indicated by the straight arrow is the direction in which the particle beam is incident. The rotary ridge filter 120 rotates in a plane perpendicular to the incident direction of the particle beam.
[0042]
FIG. 11B shows a cross-sectional view of the transmission parts 112 and 114. It can be seen that the height h of the particle beam from the incident surface is different between the transmission part 112 and the transmission part 114. More generally, if the rotary ridge filter has n types of transmission portions having different heights, it can be treated as a ridge filter having n steps in equivalent, so that the above-described effect is obtained. be able to. Note that the difference d in the height of each staircase may be the same in any transmission part, or may be different as necessary.
[0043]
【The invention's effect】
According to the present invention, it is possible to design a ridge filter that provides SOBP quickly. Furthermore, according to the present invention, even when the theoretical value and the actual measurement value are different, the actual measurement value can be brought close to the theoretical value by using the correction weight or by taking the difference into consideration in advance. You can design filters.
[0044]
Furthermore, according to the present invention, since a uniform SOBP can be obtained, a highly accurate ridge filter can be designed.
[Brief description of the drawings]
FIG. 1 is a schematic diagram of a particle beam irradiation apparatus 1 by a wobbler method using a ridge filter.
FIG. 2 is a diagram showing a dose distribution with respect to the depth of water of a particle beam 10 (FIG. 1) incident on a water target. (A) is a graph showing Bragg curves 22 and 24 by Monte Carlo calculation, a Bragg curve 24 'translated by the present invention, and an SOBP Bragg curve 26. (B) is a graph showing weighted Bragg curves 22 and 24 '' and an SOBP Bragg curve 26. FIG.
FIG. 3 is a diagram illustrating an example of a ridge filter.
FIG. 4 shows a design flow of a ridge filter according to the first embodiment.
FIG. 5 is a graph of an SOBP Bragg curve calculated from weights obtained by the method of the present invention.
FIG. 6 is a graph of theoretical values of SOBP Bragg curves.
FIG. 7 is a diagram showing an SOBP Bragg curve calculated by Monte Carlo based on weights and correction weights.
FIG. 8 is a diagram illustrating a design flow of a ridge filter according to the second embodiment.
FIG. 9 is a diagram showing an SOBP Bragg curve obtained based on weights calculated again.
FIG. 10 is a diagram showing a rotary ridge filter.
11A is a view showing a rotary ridge filter 110. FIG. FIG. 6B is a cross-sectional view of the transmission parts 112 and 114.
FIG. 12 shows a more specific appearance of a rotary ridge filter having four transmission parts.
FIG. 13 is a diagram showing an example of a tumor irradiated with a particle beam.
FIG. 14 is a diagram showing a design flow of a conventional ridge filter.
[Explanation of symbols]
DESCRIPTION OF SYMBOLS 1 Particle beam irradiation apparatus, 10 Particle beam, 11 1 pair of electromagnet for horizontal direction swing, 12 1 pair of electromagnet for vertical direction swing, 13 Range shifter, 14 Ridge filter, 16 Irradiation volume

Claims

The Bragg curve representing the relationship between the travel distance of the particle beam and the dose is used to adjust the height and width of the transmission part through which the particle beam passes, and the dose after transmission is desired within the predetermined range of the travel distance A method for designing a ridge filter that realizes uniform SOBP by value,
Obtaining a first Bragg curve of the particle beam by Monte Carlo calculation;
Obtaining n second Bragg curves by translating the first Bragg curve n times (n: natural number);
Calculating a weight to be applied to each of the first Bragg curve and the second Bragg curve based on the desired value of the SOBP, the first Bragg curve, and the second Bragg curve; And calculating the weight so that the SOBP of the desired value is obtained by superimposing the weighted first Bragg curve and each of the second Bragg curves;
Determining the width of the transmission part based on the weight;
A ridge filter design method comprising:

Performing a Monte Carlo calculation based on the weights to calculate a theoretical value of the SOBP;
Calculating a difference between the theoretical value and the desired value;
A correction weight to be applied to each of the first Bragg curve and the second Bragg curve for correcting the difference based on the difference, the first Bragg curve, and the second Bragg curve. And calculating each of the first Bragg curve and the second Bragg curve to which the correction weight is applied, and further applying the weight to the first Bragg curve and the second Bragg curve. Further calculating the correction weight so that the dose corresponding to the difference is minimized when superimposed on each of
The ridge filter design method according to claim 1, wherein the determining step further determines a width of the transmission part based on the correction weight.

Performing a Monte Carlo calculation based on the weights to calculate a theoretical value of the SOBP;
Calculating a difference between the desired value and the theoretical value;
Calculating a corrected weight obtained by correcting the weight based on a value obtained by adding the difference to the desired value, the first Bragg curve, and the second Bragg curve, the corrected weight Calculating a weight so that the SOBP of the desired value is obtained by superimposing each of the first Bragg curve and each of the second Bragg curves multiplied by
The ridge filter design method according to claim 1, wherein the determining step corrects a width of the transmission part based on the correction weight.

Calculating an error between the SOBP value obtained based on the weight and the correction weight and the desired value of the SOBP;
Selecting a position where the error exceeds a predetermined tolerance, and
The ridge filter design method according to claim 2, wherein the step of calculating the weight and the step of calculating the correction weight are performed again so as to reduce an error of the selected position.

5. The ridge filter according to claim 1, wherein the ridge filter includes the transmission unit having a height parallel to an incident direction of the particle beam and a width perpendicular to the incident direction of the particle beam. Ridge filter design method.

The ridge filter is a rotary ridge filter having at least one transmission part that rotates about an axis parallel to the incident direction, and the width of the transmission part is the width of the transmission part in the rotation direction. The ridge filter design method according to claim 5, wherein the ridge filter design method is adjusted by changing.

The ridge filter design method according to claim 6, wherein the rotary ridge filter has a plurality of transmission parts, and each of the plurality of transmission parts has a different height.