JP2015163309A - Charged particle beam irradiation device - Google Patents

Charged particle beam irradiation device Download PDF

Info

Publication number
JP2015163309A
JP2015163309A JP2015121259A JP2015121259A JP2015163309A JP 2015163309 A JP2015163309 A JP 2015163309A JP 2015121259 A JP2015121259 A JP 2015121259A JP 2015121259 A JP2015121259 A JP 2015121259A JP 2015163309 A JP2015163309 A JP 2015163309A
Authority
JP
Japan
Prior art keywords
dose distribution
charged particle
particle beam
function
irradiation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2015121259A
Other languages
Japanese (ja)
Other versions
JP6094906B2 (en
Inventor
健蔵 佐々井
Kenzo SASAI
健蔵 佐々井
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sumitomo Heavy Industries Ltd
Original Assignee
Sumitomo Heavy Industries Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sumitomo Heavy Industries Ltd filed Critical Sumitomo Heavy Industries Ltd
Priority to JP2015121259A priority Critical patent/JP6094906B2/en
Publication of JP2015163309A publication Critical patent/JP2015163309A/en
Application granted granted Critical
Publication of JP6094906B2 publication Critical patent/JP6094906B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radiation-Therapy Devices (AREA)

Abstract

PROBLEM TO BE SOLVED: To provide a charged particle beam irradiation device and a dose distribution calculation method of a charged particle beam that achieve accurate calculation of a dose distribution of the charged particle beam at an isocenter position.SOLUTION: A charged particle beam irradiation device 100 according to the invention irradiates an irradiated body at an isocenter position with a charged particle beam P, and comprises: an irradiation unit 1 that irradiates the irradiated body with the charged particle beam P; a dose distribution monitor 6 that is provided at the irradiation unit 1 and measures a dose distribution of the charged particle beam; and a dose distribution calculation unit 11 that calculates the dose distribution of the charged particle beam P at the isocenter position on the basis of a measured result by the dose distribution monitor 6. This charged particle irradiation device 100 makes it possible to calculate the dose distribution of the charged particle beam P at the isocenter position on the basis of the measured result of the dose distribution of the charged particle beam P passing through the irradiation unit 1.

Description

本発明は、荷電粒子線照射装置及び荷電粒子線の線量分布算出方法に関する。   The present invention relates to a charged particle beam irradiation apparatus and a charged particle beam dose distribution calculation method.

従来、患者の腫瘍に対して荷電粒子線の照射により治療を施す荷電粒子線照射装置が知られている。このような荷電粒子線照射装置では、アイソセンター位置に配置された腫瘍の形状に応じて荷電粒子線の線量分布を正確に設定する必要がある。   Conventionally, a charged particle beam irradiation apparatus for treating a patient's tumor by irradiation with a charged particle beam is known. In such a charged particle beam irradiation apparatus, it is necessary to accurately set the dose distribution of the charged particle beam according to the shape of the tumor arranged at the isocenter position.

例えば、特許文献１には、ＣＴ装置が得た患者の断面画像を用いて腫瘍の輪郭を決定し、腫瘍の形状から荷電粒子線の線量分布を計画する治療計画装置が開示されている。この治療計画装置では、腫瘍付近の断面画像から立体画像を構成し、計画した荷電粒子線の線量分布を立体画像へ投影することで線量分布計画の確認を行っている。   For example, Patent Literature 1 discloses a treatment planning apparatus that determines a contour of a tumor using a cross-sectional image of a patient obtained by a CT apparatus and plans a dose distribution of a charged particle beam from the shape of the tumor. In this treatment planning apparatus, a stereoscopic image is constructed from cross-sectional images near the tumor, and the dose distribution plan is confirmed by projecting the planned dose distribution of charged particle beams onto the stereoscopic image.

特開平５−３３７２０８号公報JP-A-5-337208

しかしながら、前述した治療計画装置においては、画像データ上で線量分布計画の確認を行っているが、実際に照射された荷電粒子線の線量分布を測定しておらず、確認に対する信頼性が高いとは言えない。荷電粒子線の径路上に線量分布モニタを設けることも考えられるが、この場合、測定できる線量分布は線量分布モニタ位置のものであり、アイソセンター位置における線量分布とは異なる。荷電粒子線の線量分布の確認は治療の信頼性を確保する上で重要であることから、アイソセンター位置における荷電粒子線の線量分布を精度良く確認するための技術が強く求められている。   However, in the treatment planning apparatus described above, the dose distribution plan is confirmed on the image data, but the dose distribution of the actually irradiated charged particle beam is not measured, and the confirmation is highly reliable. I can't say that. Although it is conceivable to provide a dose distribution monitor on the path of the charged particle beam, in this case, the measurable dose distribution is at the dose distribution monitor position and is different from the dose distribution at the isocenter position. Since confirmation of the dose distribution of the charged particle beam is important for ensuring the reliability of treatment, there is a strong demand for a technique for accurately confirming the dose distribution of the charged particle beam at the isocenter position.

そこで、本発明は、アイソセンター位置における荷電粒子線の線量分布を精度良く算出できる荷電粒子線照射装置及び荷電粒子線の線量分布算出方法を提供することを目的とする。   Therefore, an object of the present invention is to provide a charged particle beam irradiation apparatus and a charged particle beam dose distribution calculation method capable of accurately calculating a dose distribution of a charged particle beam at an isocenter position.

上記課題を解決するため、本発明は、アイソセンター位置に配置された被照射体に対して荷電粒子線を照射する荷電粒子線照射装置であって、被照射体に荷電粒子線を照射する照射部と、照射部に設けられ、荷電粒子線の線量分布を測定する線量分布測定手段と、線量分布測定手段の測定結果に基づいて、アイソセンター位置における荷電粒子線の線量分布を算出する線量分布算出手段と、を備えることを特徴とする。   In order to solve the above-described problems, the present invention is a charged particle beam irradiation apparatus that irradiates a charged particle beam to an irradiation object arranged at an isocenter position, and irradiates the irradiation object with the charged particle beam. A dose distribution measurement means for measuring a dose distribution of a charged particle beam provided in a radiation section, and a dose distribution for calculating a dose distribution of a charged particle beam at an isocenter position based on a measurement result of the dose distribution measurement means And a calculating means.

本発明に係る荷電粒子線照射装置によれば、照射部を通過する荷電粒子線の線量分布の測定結果に基づいて、アイソセンター位置における荷電粒子線の線量分布を算出することができる。従って、この荷電粒子線照射装置によれば、実際に照射することでアイソセンター位置における線量分布を治療前に確認することが可能になる。しかも、この荷電粒子線照射装置では、実際に照射された荷電粒子線の線量分布の測定結果を用いることで、実際の測定結果を用いない場合と比べて、アイソセンター位置における荷電粒子線の線量分布を精度良く算出することができる。   According to the charged particle beam irradiation apparatus according to the present invention, the dose distribution of the charged particle beam at the isocenter position can be calculated based on the measurement result of the dose distribution of the charged particle beam passing through the irradiation unit. Therefore, according to this charged particle beam irradiation apparatus, the dose distribution at the isocenter position can be confirmed before treatment by actually irradiating. Moreover, in this charged particle beam irradiation apparatus, the dose of the charged particle beam at the isocenter position is compared with the case where the actual measurement result is not used by using the measurement result of the dose distribution of the actually irradiated charged particle beam. The distribution can be calculated with high accuracy.

本発明に係る荷電粒子線照射装置において、照射部に設けられ、荷電粒子線を走査する電磁石を更に備え、線量分布算出手段は、照射部に入射した荷電粒子線の線量分布を荷電粒子線の照射方向と直交する方向の位置座標の関数とし、当該関数に対して電磁石の走査に対応する関数を畳み込み積分して得られる線量分布測定手段の位置における荷電粒子線の線量分布の関数と、線量分布測定手段の測定結果から得られる線量分布の中心位置座標と、に基づいて、アイソセンター位置における荷電粒子線の線量分布を算出してもよい。   The charged particle beam irradiation apparatus according to the present invention further includes an electromagnet that is provided in the irradiation unit and scans the charged particle beam, and the dose distribution calculating means calculates the dose distribution of the charged particle beam incident on the irradiation unit of the charged particle beam. A function of the position coordinate in the direction orthogonal to the irradiation direction, and a function of the dose distribution of the charged particle beam at the position of the dose distribution measuring means obtained by convolving and integrating the function corresponding to the scanning of the electromagnet with the function, and the dose The dose distribution of the charged particle beam at the isocenter position may be calculated based on the center position coordinates of the dose distribution obtained from the measurement result of the distribution measuring means.

この場合、線量分布測定手段の位置における線量分布を表す関数に対して実際の荷電粒子線の線量分布の測定結果を利用することで、アイソセンター位置における荷電粒子線の線量分布を算出することができる。   In this case, the dose distribution of the charged particle beam at the isocenter position can be calculated by using the measurement result of the actual charged particle beam dose distribution with respect to the function representing the dose distribution at the position of the dose distribution measuring means. it can.

本発明に係る荷電粒子線照射装置において、照射方向と直交する方向をＸ軸方向とした場合に、照射部に入射した荷電粒子線の線量分布のＸ軸位置座標の関数ｆ（ｘ）が下記の式（１）で表され、電磁石の走査に対応する関数ｇ（ｘ）が下記の式（２）で表され、線量分布測定手段の位置における荷電粒子線の線量分布の関数ＤＭ（ｘ）が下記の式（３）で表されても良い。

Figure 2015163309
Figure 2015163309
Figure 2015163309
 但し、式（１）におけるＡは規格化定数、σは線量分布のＸ軸方向の広がり、μは線量分布のＸ軸上の中心位置座標である。また、式（２）におけるＲはワブリング半径である。 In the charged particle beam irradiation apparatus according to the present invention, when the direction orthogonal to the irradiation direction is the X-axis direction, the function f (x) of the X-axis position coordinate of the dose distribution of the charged particle beam incident on the irradiation unit is The function g (x) corresponding to the scanning of the electromagnet is expressed by the following expression (2), and the function D M (x of the dose distribution of the charged particle beam at the position of the dose distribution measuring means ) May be represented by the following formula (3).
Figure 2015163309
Figure 2015163309
Figure 2015163309
 In Equation (1), A is a normalization constant, σ is the spread of the dose distribution in the X-axis direction, and μ is the center position coordinate on the X-axis of the dose distribution. Moreover, R in Formula (2) is a wobbling radius.

この場合、線量分布を表す関数ＤＭ（ｘ）に対して実際の荷電粒子線の線量分布の測定結果を利用することで、アイソセンター位置における荷電粒子線の線量分布を精度良く算出することができる。 In this case, it is possible to accurately calculate the dose distribution of the charged particle beam at the isocenter position by using the measurement result of the actual dose distribution of the charged particle beam with respect to the function D M (x) representing the dose distribution. it can.

本発明は、アイソセンター位置に配置された被照射体に対して照射される荷電粒子線の線量分布算出方法であって、被照射体に荷電粒子線を照射する照射部に設けられた線量分布モニタによって荷電粒子線の線量分布を測定する線量分布測定ステップと、線量分布測定ステップの測定結果に基づいて、アイソセンター位置における荷電粒子線の線量分布を算出する線量分布算出ステップと、を有することを特徴とする。   The present invention relates to a method for calculating a dose distribution of a charged particle beam irradiated to an irradiated object arranged at an isocenter position, the dose distribution provided in an irradiation unit for irradiating the irradiated object with the charged particle beam. A dose distribution measuring step for measuring the dose distribution of the charged particle beam by a monitor; and a dose distribution calculating step for calculating the dose distribution of the charged particle beam at the isocenter position based on the measurement result of the dose distribution measuring step. It is characterized by.

本発明に係る荷電粒子線の線量分布算出方法において、線量分布測定ステップでは、走査された荷電粒子線の線量分布を測定し、線量分布算出ステップでは、照射部に入射した荷電粒子線の線量分布を荷電粒子線の照射方向と直交する方向の位置座標の関数の関数とし、当該関数に対して電磁石の走査に対応する関数を畳み込み積分して得られる線量分布測定位置における荷電粒子線の線量分布の関数と、線量分布測定ステップの測定結果から得られる線量分布の中心位置座標と、に基づいて、アイソセンター位置における荷電粒子線の線量分布を算出しても良い。   In the charged particle beam dose distribution calculating method according to the present invention, in the dose distribution measuring step, the dose distribution of the scanned charged particle beam is measured, and in the dose distribution calculating step, the dose distribution of the charged particle beam incident on the irradiation unit. Is a function of a position coordinate function in a direction perpendicular to the irradiation direction of the charged particle beam, and the dose distribution of the charged particle beam at the dose distribution measurement position obtained by convolving and integrating the function corresponding to the scanning of the electromagnet with the function The dose distribution of the charged particle beam at the isocenter position may be calculated on the basis of the above function and the center position coordinates of the dose distribution obtained from the measurement result of the dose distribution measurement step.

この場合、線量分布測定位置における線量分布を表す関数に対して実際の荷電粒子線の線量分布の測定結果を利用することで、アイソセンター位置における荷電粒子線の線量分布を算出することができる。   In this case, the dose distribution of the charged particle beam at the isocenter position can be calculated by using the measurement result of the actual dose distribution of the charged particle beam with respect to the function representing the dose distribution at the dose distribution measurement position.

本発明に係る荷電粒子線の線量分布算出方法において、照射方向と直交する方向をＸ軸方向とした場合に、照射部に入射した荷電粒子線の線量分布のＸ軸位置座標の関数ｆ（ｘ）が下記の式（１）で表され、電磁石の走査に対応する関数ｇ（ｘ）が下記の式（２）で表され、線量分布測定位置における荷電粒子線の線量分布の関数ＤＭ（ｘ）が下記の式（３）で表されても良い。

Figure 2015163309
Figure 2015163309
Figure 2015163309
 但し、式（１）におけるＡは規格化定数、σは線量分布のＸ軸方向の広がり、μは線量分布のＸ軸上の中心位置座標である。また、式（２）におけるＲはワブリング半径である。 In the charged particle beam dose distribution calculating method according to the present invention, when the direction orthogonal to the irradiation direction is the X-axis direction, the function f (x) of the X-axis position coordinate of the dose distribution of the charged particle beam incident on the irradiation unit is used. ) Is expressed by the following formula (1), a function g (x) corresponding to scanning of the electromagnet is expressed by the following formula (2), and the function D M ( x) may be represented by the following formula (3).
Figure 2015163309
Figure 2015163309
Figure 2015163309
 In Equation (1), A is a normalization constant, σ is the spread of the dose distribution in the X-axis direction, and μ is the center position coordinate on the X-axis of the dose distribution. Moreover, R in Formula (2) is a wobbling radius.

この場合、線量分布を表す関数ＤＭ（ｘ）に対して実際の荷電粒子線の線量分布の測定結果を利用することで、アイソセンター位置における荷電粒子線の線量分布を精度良く算出することができる。 In this case, it is possible to accurately calculate the dose distribution of the charged particle beam at the isocenter position by using the measurement result of the actual dose distribution of the charged particle beam with respect to the function D M (x) representing the dose distribution. it can.

本発明によれば、アイソセンター位置における荷電粒子線の線量分布を精度良く算出できる。   According to the present invention, the dose distribution of the charged particle beam at the isocenter position can be calculated with high accuracy.

本発明に係る荷電粒子線照射装置の一実施形態を示す斜視図である。1 is a perspective view showing an embodiment of a charged particle beam irradiation apparatus according to the present invention. 図１に示すワブラー方式の照射部を示す斜視図である。It is a perspective view which shows the irradiation part of the wobbler system shown in FIG. 照射部における荷電粒子線の線量分布の変化を示す概略図である。It is the schematic which shows the change of the dose distribution of the charged particle beam in an irradiation part. 線量分布モニタにおける測定結果の一例を示すグラフである。It is a graph which shows an example of the measurement result in a dose distribution monitor. 本発明に係る荷電粒子線の線量分布算出方法を示すフローチャートである。It is a flowchart which shows the dose distribution calculation method of the charged particle beam which concerns on this invention.

以下、本発明の好適な実施形態について、図面を参照して詳細に説明する。   DESCRIPTION OF EMBODIMENTS Hereinafter, preferred embodiments of the present invention will be described in detail with reference to the drawings.

図１に示されるように、照射部１は、荷電粒子線照射装置１００において、治療台１０５を取り囲むように設けられた回転ガントリ１０３に取り付けられ、この回転ガントリ１０３によって治療台１０５の回りに回転可能とされている。そして、治療台１０５に寝かされた患者の腫瘍等の被照射体に対して、荷電粒子線Ｐ（図２参照）を照射する。荷電粒子線Ｐは、電荷をもった粒子を高速に加速したものであり、荷電粒子線Ｐとしては、例えば陽子線、重粒子（重イオン）線等が挙げられる。   As shown in FIG. 1, the irradiation unit 1 is attached to a rotating gantry 103 provided so as to surround the treatment table 105 in the charged particle beam irradiation apparatus 100, and is rotated around the treatment table 105 by the rotation gantry 103. It is possible. Then, a charged particle beam P (see FIG. 2) is irradiated to an irradiated object such as a tumor of a patient who is laid on the treatment table 105. The charged particle beam P is obtained by accelerating charged particles at a high speed. Examples of the charged particle beam P include proton beams and heavy particle (heavy ion) beams.

なお、図１には示されていないが、荷電粒子線照射装置１００は、イオン源で生成した荷電粒子線を加速して荷電粒子線Ｐを出射するサイクロトロンを治療台１０５及び回転ガントリ１０３から離れた位置に備えている。サイクロトロンから出射された荷電粒子線Ｐはビーム輸送系を介して照射部１に供給される。この荷電粒子線照射装置１００は、制御装置１０によって統括的に制御されている。   Although not shown in FIG. 1, the charged particle beam irradiation apparatus 100 separates the cyclotron that accelerates the charged particle beam generated by the ion source and emits the charged particle beam P away from the treatment table 105 and the rotating gantry 103. Prepared in a different position. The charged particle beam P emitted from the cyclotron is supplied to the irradiation unit 1 through the beam transport system. The charged particle beam irradiation apparatus 100 is comprehensively controlled by the control apparatus 10.

図２に示されるように、照射部１は、荷電粒子線Ｐの照射方向に順に配列された四極磁石２、ワブリング磁石３、散乱体４、リッジフィルタ５、線量分布モニタ（線量分布測定手段）６、マルチリーフコリメータ７、及びスノート８を備えている。以下、照射部１を通る荷電粒子線Ｐの照射方向をＺ軸方向、Ｚ軸に直交する二方向をＸ軸方向及びＹ軸方向として説明に用いる。Ｘ軸方向及びＹ軸方向は互いに直交する方向である。   As shown in FIG. 2, the irradiation unit 1 includes a quadrupole magnet 2, a wobbling magnet 3, a scatterer 4, a ridge filter 5 and a dose distribution monitor (dose distribution measuring means) arranged in order in the irradiation direction of the charged particle beam P. 6, a multi-leaf collimator 7 and a snout 8 are provided. Hereinafter, the irradiation direction of the charged particle beam P passing through the irradiation unit 1 will be used as an explanation in the Z axis direction, and two directions orthogonal to the Z axis will be used as an X axis direction and a Y axis direction. The X-axis direction and the Y-axis direction are directions orthogonal to each other.

照射部１の入口に最も近い四極磁石２は、ビーム輸送ラインを介して入力された荷電粒子線Ｐが発散するのを抑えて収束させる電磁石である。四極磁石２の下流には、ワブリング磁石３が配置されている。ワブリング磁石３は、正規分布状の荷電粒子線Ｐを円軌道で走査することで一様な照射野を形成するための電磁石である。   The quadrupole magnet 2 closest to the entrance of the irradiating unit 1 is an electromagnet that converges while suppressing the divergence of the charged particle beam P input through the beam transport line. A wobbling magnet 3 is disposed downstream of the quadrupole magnet 2. The wobbling magnet 3 is an electromagnet for forming a uniform irradiation field by scanning a normally distributed charged particle beam P in a circular orbit.

散乱体４は、ワブリング磁石３で走査された荷電粒子線Ｐの拡散を行う。このような散乱体４は、例えば鉛製の板から構成されている。散乱体４は、荷電粒子線Ｐを照射方向と直交する方向（ＸＹ平面内の方向）に広がりを持つ幅広のビームへと拡散させる。   The scatterer 4 diffuses the charged particle beam P scanned by the wobbling magnet 3. Such a scatterer 4 is composed of, for example, a lead plate. The scatterer 4 diffuses the charged particle beam P into a wide beam having a spread in a direction orthogonal to the irradiation direction (direction in the XY plane).

リッジフィルタ５は、散乱体４で散乱された荷電粒子線Ｐの線量分布の調整を行う。具体的には、リッジフィルタ５は、患者の体内の腫瘍の厚さ（照射方向における腫瘍の長さ）に対応するように、荷電粒子線Ｐに拡大ブラッグピーク（ＳＯＢＰ）を与える。   The ridge filter 5 adjusts the dose distribution of the charged particle beam P scattered by the scatterer 4. Specifically, the ridge filter 5 gives an enlarged Bragg peak (SOBP) to the charged particle beam P so as to correspond to the thickness of the tumor in the patient's body (the length of the tumor in the irradiation direction).

線量分布モニタ６は、通過する荷電粒子線Ｐの線量分布の測定を行う。線量分布モニタ６は、Ｘ軸方向又はＹ軸方向で延在して格子状に配列された複数の金属ワイヤを有しており、荷電粒子線Ｐが照射された金属ワイヤから発生する電子を検出する。これらの複数の金属線にはＸ軸及びＹ軸の位置座標に応じたチャネル番号が付されており、各チャネル番号の金属ワイヤから発生した電子を検出することで、荷電粒子線Ｐの線量分布の測定が行われる。   The dose distribution monitor 6 measures the dose distribution of the charged particle beam P that passes therethrough. The dose distribution monitor 6 has a plurality of metal wires extending in the X-axis direction or the Y-axis direction and arranged in a lattice pattern, and detects electrons generated from the metal wires irradiated with the charged particle beam P. To do. Channel numbers corresponding to the position coordinates of the X axis and the Y axis are assigned to the plurality of metal wires, and the dose distribution of the charged particle beam P is detected by detecting electrons generated from the metal wires of the respective channel numbers. Is measured.

マルチリーフコリメータ７は、照射方向と垂直な平面（ＸＹ平面）における荷電粒子線Ｐの形状の成形を行う。マルチリーフコリメータ７は、Ｘ軸方向で対向する二列の櫛歯を有しており、対向する櫛歯間の開口を調整することで、開口を通過する荷電粒子線Ｐの成形を行う。なお、マルチリーフコリメータ７に代えて、ブロックコリメータを用いても良い。ブロックコリメータとしては、患者の腫瘍形状に合わせて穴を形成された金属柱等が用いられる。   The multi-leaf collimator 7 forms the shape of the charged particle beam P on a plane (XY plane) perpendicular to the irradiation direction. The multi-leaf collimator 7 has two rows of comb teeth facing each other in the X-axis direction, and the charged particle beam P passing through the opening is formed by adjusting the opening between the facing comb teeth. Instead of the multi-leaf collimator 7, a block collimator may be used. As the block collimator, a metal column or the like in which a hole is formed according to the tumor shape of the patient is used.

スノート８は、照射方向（Ｚ軸方向）の奥行きについて荷電粒子線Ｐを腫瘍の形状に合わせて成形するボーラスを保持する部材である。ボーラスは、荷電粒子線Ｐの最大到達深さの部分の立体形状を、腫瘍の最大深さ部分の形状に合わせて成形する。   The snout 8 is a member that holds a bolus for shaping the charged particle beam P in accordance with the shape of the tumor with respect to the depth in the irradiation direction (Z-axis direction). The bolus is formed so that the three-dimensional shape of the portion where the charged particle beam P reaches the maximum depth is matched with the shape of the maximum depth portion of the tumor.

図１及び図２に示されるように、荷電粒子線照射装置１００は、その全体的な制御を行う制御装置１０を備えている。制御装置１０は、患者に対する荷電粒子線Ｐの照射制御やマルチリーフコリメータ７の位置制御を行う。また、制御装置１０は、アイソセンター位置における荷電粒子線Ｐの線量分布を算出するための線量分布算出部（線量分布算出手段）１１を有している。線量分布算出部１１は、線量分布モニタ６の測定結果に基づいて、アイソセンター位置における荷電粒子線Ｐの線量分布を算出する。   As shown in FIGS. 1 and 2, the charged particle beam irradiation apparatus 100 includes a control device 10 that performs overall control thereof. The control device 10 performs irradiation control of the charged particle beam P to the patient and position control of the multi-leaf collimator 7. Further, the control device 10 has a dose distribution calculation unit (dose distribution calculation means) 11 for calculating the dose distribution of the charged particle beam P at the isocenter position. The dose distribution calculation unit 11 calculates the dose distribution of the charged particle beam P at the isocenter position based on the measurement result of the dose distribution monitor 6.

次に、本実施形態に係る荷電粒子線照射装置１００における線量分布算出方法について説明する。   Next, a dose distribution calculation method in the charged particle beam irradiation apparatus 100 according to the present embodiment will be described.

図３に示されるように、本実施形態に係る線量分布算出方法は、線量分布モニタ６が荷電粒子線Ｐの線量分布を測定する線量分布測定ステップＳ１と、線量分布算出部１１がアイソセンター位置における荷電粒子線Ｐの線量分布を算出する線量分布算出ステップＳ２と、を有している。   As shown in FIG. 3, the dose distribution calculation method according to the present embodiment includes a dose distribution measurement step S1 in which the dose distribution monitor 6 measures the dose distribution of the charged particle beam P, and the dose distribution calculation unit 11 in the isocenter position. A dose distribution calculating step S2 for calculating the dose distribution of the charged particle beam P at.

ここで、図４は、ＸＺ平面における荷電粒子線Ｐの線量分布の変化を説明するための概略図である。図４において、四極磁石２、リッジフィルタ５、及びスノート８の描写を省略する。照射部１における中心軸（照射部１を構成する機器２〜８の略中心を通る軸）をＣとして示す。   Here, FIG. 4 is a schematic diagram for explaining a change in the dose distribution of the charged particle beam P in the XZ plane. In FIG. 4, the depiction of the quadrupole magnet 2, the ridge filter 5, and the snout 8 is omitted. A central axis in the irradiation unit 1 (an axis passing through the approximate center of the devices 2 to 8 constituting the irradiation unit 1) is indicated as C.

図４に示されるように、ビーム輸送ラインから入射した荷電粒子線Ｐは、四極磁石２により収束され、その線量分布はＷ１のような山なりの状態となる。この場合の線量分布Ｗ１は、位置座標ｘの関数ｆ（ｘ）として下記の式（１）で表すことができる。式（１）におけるＡは規格化定数である。また、μは実際の線量分布の中心位置座標であり、中心軸Ｃからの距離（Ｘ軸方向における距離）として表すことができる。σは線量分布の広がりである。この関数ｆ（ｘ）はガウシアン（ガウス関数）となる。

Figure 2015163309
As shown in FIG. 4, the charged particle beam P incident from the beam transport line is converged by the quadrupole magnet 2, and the dose distribution is in a mountainous state like W <b> 1. The dose distribution W1 in this case can be expressed by the following equation (1) as a function f (x) of the position coordinate x. A in Formula (1) is a normalization constant. Further, μ is the center position coordinate of the actual dose distribution, and can be expressed as a distance from the center axis C (a distance in the X-axis direction). σ is the spread of the dose distribution. This function f (x) is a Gaussian (Gaussian function).
Figure 2015163309

線量分布Ｗ１の荷電粒子線Ｐがワブリング磁石３によって走査されると、線量分布はＷ２の状態（二つの山からなる状態）となる。この場合のワブリング半径をＲとする。ワブリング半径Ｒは、入射した荷電粒子線Ｐの粒子エネルギーＥ、ワブリング磁石３の磁場強度Ｂ、ワブリング磁石３から線量分布Ｗ２の現在位置までの距離Ｌｗに依存するパラメータである。ここで、線形分布Ｗ１に対するワブリング磁石３の走査の影響を表す関数ｇ（ｘ）を検討すると、下記の式（２）として表すことができる。この場合、線量分布Ｗ２の関数ＤＭ（ｘ）は、下記の式（３）として表すことができる。すなわち、関数ＤＭ（ｘ）は、関数ｆ（ｘ）及び関数ｇ（ｘ）を畳み込み積分することで求めることができる。

Figure 2015163309
Figure 2015163309
 When the charged particle beam P of the dose distribution W1 is scanned by the wobbling magnet 3, the dose distribution is in the state of W2 (a state consisting of two peaks). Let R be the wobbling radius in this case. The wobbling radius R is a parameter that depends on the particle energy E of the incident charged particle beam P, the magnetic field strength B of the wobbling magnet 3, and the distance Lw from the wobbling magnet 3 to the current position of the dose distribution W2. Here, when the function g (x) representing the influence of the scanning of the wobbling magnet 3 on the linear distribution W1 is examined, it can be expressed as the following equation (2). In this case, the function D M (x) of the dose distribution W2 can be expressed as the following formula (3). That is, the function D M (x) can be obtained by convolving and integrating the function f (x) and the function g (x).
Figure 2015163309
Figure 2015163309

その後、荷電粒子線Ｐが散乱体４によってＸ軸方向及びＹ軸方向に拡散されると、線量分布はＷ３の状態となる。線量分布Ｗ３は、線量分布モニタ６の位置（線量分布測定位置）における荷電粒子線Ｐの線量分布である。この線量分布Ｗ３を表す関数は、上記式（２）において、関数ｆ（ｘ）に含まれるσを線量分布モニタ６の位置における広がりσＭに置き換えることで求められる。このσＭは、入射時における線形分布Ｗ１の広がりσ０、散乱体４の厚さｄ、散乱体４から線量分布モニタ６までの距離ＬＳＭに依存するパラメータである。この線量分布Ｗ３が線量分布測定ステップＳ１において線量分布モニタ６に測定される。 Thereafter, when the charged particle beam P is diffused by the scatterer 4 in the X-axis direction and the Y-axis direction, the dose distribution becomes W3. The dose distribution W3 is a dose distribution of the charged particle beam P at the position of the dose distribution monitor 6 (dose distribution measurement position). The function representing the dose distribution W3 is obtained by replacing σ included in the function f (x) with the spread σ M at the position of the dose distribution monitor 6 in the above equation (2). This σ M is a parameter that depends on the spread σ 0 of the linear distribution W 1 at the time of incidence, the thickness d of the scatterer 4, and the distance L SM from the scatterer 4 to the dose distribution monitor 6. This dose distribution W3 is measured by the dose distribution monitor 6 in the dose distribution measurement step S1.

図５は、線量分布モニタ６の測定結果の一例を示すグラフである。横軸は、線量分布モニタ６のＸ軸座標に対応するチャネル番号、縦軸は線量相当値を示している。線量相当値とは、荷電粒子線Ｐの線量に比例する値であり、所定の係数を乗じることで線量が得られる。   FIG. 5 is a graph showing an example of the measurement result of the dose distribution monitor 6. The horizontal axis represents the channel number corresponding to the X-axis coordinate of the dose distribution monitor 6, and the vertical axis represents the dose equivalent value. The dose equivalent value is a value proportional to the dose of the charged particle beam P, and the dose can be obtained by multiplying by a predetermined coefficient.

線量分布算出部１１は、図５に示す測定結果の線量分布データ（折れ線状のデータ）に対してフィッティングを行う。このようなフィッティングには、Levenberg−Marquardt法やNelder−Mead Simplex法等を利用することができる。フィッティング後の線量分布データの一例を図５に一点鎖線として示す。   The dose distribution calculation unit 11 performs fitting on the dose distribution data (string-line data) of the measurement results shown in FIG. For such fitting, the Levenberg-Marquardt method, the Nelder-Mead Simplex method, or the like can be used. An example of the dose distribution data after fitting is shown as a one-dot chain line in FIG.

ここで、線量モニタ６の位置における線量分布の広がりσＭと線量モニタ６の位置におけるワブリング半径ＲＭは、予め計算や実測で決定することができる。従って、上記の式（３）におけるフリーパラメータは、規格化定数Ａと線量モニタ６の位置における線量分布の中心位置座標μＭのみである。Ａは規格化定数なので、物理的な意味はない。このため、フィッティングによってμＭが求められる。線量分布算出部１１は、フィッティングによって求めたμＭを、アイソセンター位置における線量分布Ｗ４の中心位置座標μＩとして用いる。 Here, wobbling radius R M in the spread sigma M and the position of the dose monitor 6 of the dose distribution at the position of the dose monitor 6 can be determined in advance by calculation or actual measurement. Thus, the free parameters in the above equation (3) is only the center position coordinates mu M of dose distribution at the position of the normalization constant A and a dose monitor 6. Since A is a normalization constant, there is no physical meaning. For this reason, the μ M is obtained by fitting. Dose distribution calculation unit 11, a mu M as determined by fitting, used as the center position coordinates mu I dose distribution W4 at isocenter position.

アイソセンター位置における線量分布Ｗ４は、上記の式（３）にマルチリーフコリメータ７の効果を乗じることで求められる。マルチリーフコリメータ７の効果は、下記の式（４），（５）として表わすことができる。下記の式（４），（５）において、ａはマルチリーフコリメータ７の開口の幅ｑ（Ｘ軸方向における幅）等に依存するパラメータである。また、ｓは、入射した荷電粒子線Ｐの粒子エネルギーＥ、ワブリング磁石３の磁場強度Ｂ、マルチリーフコリメータ７からアイソセンター位置までの距離ＬＣＩ等に依存するパラメータである。

Figure 2015163309
The dose distribution W4 at the isocenter position is obtained by multiplying the above expression (3) by the effect of the multi-leaf collimator 7. The effect of the multi-leaf collimator 7 can be expressed as the following equations (4) and (5). In the following formulas (4) and (5), a is a parameter depending on the opening width q (width in the X-axis direction) of the multi-leaf collimator 7 and the like. S is a parameter that depends on the particle energy E of the incident charged particle beam P, the magnetic field strength B of the wobbling magnet 3, the distance LCI from the multi-leaf collimator 7 to the isocenter position, and the like.
Figure 2015163309

アイソセンター位置における荷電粒子線Ｐの線量分布Ｗ４を表す関数ＤＩ（ｘ）は、ｆ（ｘ）とｇ（ｘ）とを畳み込み積分をしたものにｈ（ｘ）を乗じたものとなる。従って、関数ＤＩ（ｘ）を下記の式（６）として表すことができる。式（６）において、アイソセンター位置における線量分布Ｗ４の中心位置座標μＩにはフィッティングで求めたμＭを用いる。また、アイソセンター位置における線量分布Ｗ４の広がりσＩ及びアイソセンター位置におけるワブリング半径ＲＩには計算や実測で決定した値を利用する。

Figure 2015163309
The function D I (x) representing the dose distribution W4 of the charged particle beam P at the isocenter position is obtained by convolving and integrating f (x) and g (x) with h (x). Therefore, the function D I (x) can be expressed as the following equation (6). In the formula (6), the center position coordinates mu I dose distribution W4 at isocenter position using mu M obtained in fitting. Further, values determined by calculation or measurement are used for the spread σ I of the dose distribution W4 at the isocenter position and the wobbling radius R I at the isocenter position.
Figure 2015163309

線量分布算出部１１は、線量分布算出ステップＳ２において上記の式（６）に基づいて、アイソセンター位置における荷電粒子線Ｐの線量分布Ｗ４の関数ＤＩ（ｘ）を算出する。線量分布算出部１１は、線量モニタ６の測定結果のフィッティングデータと上記の式（６）とに基づいて、アイソセンター位置における荷電粒子線Ｐの線量分布Ｗ４を算出することができる。なお、線量分布をＸ軸位置座標の関数とした場合の算出について説明したが、線量分布をＹ軸位置座標の関数とした場合も同様に算出することができる。 The dose distribution calculation unit 11 calculates the function D I (x) of the dose distribution W4 of the charged particle beam P at the isocenter position based on the above formula (6) in the dose distribution calculation step S2. The dose distribution calculation unit 11 can calculate the dose distribution W4 of the charged particle beam P at the isocenter position based on the fitting data of the measurement result of the dose monitor 6 and the above equation (6). Although the calculation when the dose distribution is a function of the X-axis position coordinates has been described, the calculation can be similarly performed when the dose distribution is a function of the Y-axis position coordinates.

以上説明した荷電粒子線照射装置１００及び荷電粒子線の線量分布算出方法によれば、照射部１を通過する荷電粒子線Ｐの線量分布の測定結果に基づいて、アイソセンター位置における荷電粒子線Ｐの線量分布を算出することができる。従って、この荷電粒子線照射装置１００及び線量分布算出方法によれば、実際に照射することでアイソセンター位置における線量分布を治療前に確認することが可能になる。しかも、この荷電粒子線照射装置１００及び線量分布算出方法では、実際に照射部１を通過する荷電粒子線Ｐの線量分布の測定結果を用いるので、線量分布の測定結果を用いない場合と比べて、アイソセンター位置における荷電粒子線Ｐの線量分布を精度良く算出することができる。   According to the charged particle beam irradiation apparatus 100 and the charged particle beam dose distribution calculation method described above, the charged particle beam P at the isocenter position is based on the measurement result of the dose distribution of the charged particle beam P passing through the irradiation unit 1. The dose distribution can be calculated. Therefore, according to the charged particle beam irradiation apparatus 100 and the dose distribution calculation method, the dose distribution at the isocenter position can be confirmed before treatment by actually irradiating. Moreover, in the charged particle beam irradiation apparatus 100 and the dose distribution calculation method, since the measurement result of the dose distribution of the charged particle beam P that actually passes through the irradiation unit 1 is used, the measurement result of the dose distribution is not used. The dose distribution of the charged particle beam P at the isocenter position can be calculated with high accuracy.

また、この荷電粒子線照射装置１００及び線量分布算出方法によれば、上述した式（１）〜（６）を用い、実際の荷電粒子線Ｐの線量分布の測定結果を踏まえて線量分布を算出することで、アイソセンター位置における荷電粒子線Ｐの線量分布を一層精度良く算出することができる。   In addition, according to the charged particle beam irradiation apparatus 100 and the dose distribution calculation method, the dose distribution is calculated based on the measurement result of the actual dose distribution of the charged particle beam P using the above-described equations (1) to (6). Thus, the dose distribution of the charged particle beam P at the isocenter position can be calculated with higher accuracy.

本発明は、上述した実施形態に限定されるものではない。例えば、アイソセンター位置における線量分布の中心位置座標μＩは、フィッティングにより求められたμＭと同じにする必要はない。照射部１の入口側のプロファイルモニタによって当該プロファイルモニタの位置における線量分布の中心位置座標μｐを測定した場合には、μｐを用いてアイソセンター位置における線量分布の中心位置座標μＩを算出することもできる。具体的には、μｐと、プロファイルモニタから線量分布モニタ６までの距離ＬＰＭと、線量分布モニタ６からアイソセンター位置までの距離ＬＭＩと、を用いて下記の式（７）から求めることができる。

Figure 2015163309
The present invention is not limited to the embodiment described above. For example, the center position coordinates mu I of dose distribution in isocenter position need not be the same as mu M obtained by the fitting. When measuring the center position coordinates μp dose distribution at the position of the profile monitored by the inlet side of the profile monitor the irradiation unit 1 is to calculate the center position coordinates mu I of dose distribution in the isocenter position with μp You can also. Specifically, it can be obtained from the following equation (7) using μp, the distance L PM from the profile monitor to the dose distribution monitor 6, and the distance L MI from the dose distribution monitor 6 to the isocenter position. it can.
Figure 2015163309

また、本発明は、ワブラー法による荷電粒子線の線量分布算出に限られない。例えば、ワブリング磁石を使わずに、複数の散乱体により荷電粒子線を拡散させる二重散乱体法においても本発明を有効に適用することができる。また、スキャニング法による荷電粒子線の線量分布算出にも適用することができる。更に、本発明は、上述したような関数を用いた算出方法に限られず、様々な方法によりアイソセンター位置の線量分布の算出を行うことができる。   The present invention is not limited to the calculation of the dose distribution of charged particle beams by the wobbler method. For example, the present invention can be effectively applied to a double scatterer method in which a charged particle beam is diffused by a plurality of scatterers without using a wobbling magnet. It can also be applied to the calculation of the dose distribution of charged particle beams by the scanning method. Furthermore, the present invention is not limited to the calculation method using the function as described above, and the dose distribution at the isocenter position can be calculated by various methods.

１…照射部 ２…四極磁石 ３…ワブリング磁石 ４…散乱体 ５…リッジフィルタ ６…線量分布モニタ ７…マルチリーフコリメータ ８…スノート １０…制御装置 １１…線量分布算出部（線量分布算出手段） １００…荷電粒子線照射装置 １０３…回転ガントリ １０５…治療台 Ｃ…中心軸 Ｐ…荷電粒子線 ｑ…開口幅 Ｒ…ワブリング半径 Ｗ１-Ｗ４…線形分布   DESCRIPTION OF SYMBOLS 1 ... Irradiation part 2 ... Quadrupole magnet 3 ... Wobbling magnet 4 ... Scattering body 5 ... Ridge filter 6 ... Dose distribution monitor 7 ... Multi-leaf collimator 8 ... Snout 10 ... Control apparatus 11 ... Dose distribution calculation part (dose distribution calculation means) 100 ... charged particle beam irradiation device 103 ... rotating gantry 105 ... treatment table C ... central axis P ... charged particle beam q ... aperture width R ... wobbling radius W1-W4 ... linear distribution

本発明は、荷電粒子線照射装置に関する。 The present invention relates to a charged particle beam irradiation equipment.

従来、患者の腫瘍に対して荷電粒子線の照射により治療を施す荷電粒子線照射装置が知られている。このような荷電粒子線照射装置では、アイソセンター位置に配置された腫瘍の形状に応じて荷電粒子線の線量分布を正確に設定する必要がある。   Conventionally, a charged particle beam irradiation apparatus for treating a patient's tumor by irradiation with a charged particle beam is known. In such a charged particle beam irradiation apparatus, it is necessary to accurately set the dose distribution of the charged particle beam according to the shape of the tumor arranged at the isocenter position.

例えば、特許文献１には、ＣＴ装置が得た患者の断面画像を用いて腫瘍の輪郭を決定し、腫瘍の形状から荷電粒子線の線量分布を計画する治療計画装置が開示されている。この治療計画装置では、腫瘍付近の断面画像から立体画像を構成し、計画した荷電粒子線の線量分布を立体画像へ投影することで線量分布計画の確認を行っている。   For example, Patent Literature 1 discloses a treatment planning apparatus that determines a contour of a tumor using a cross-sectional image of a patient obtained by a CT apparatus and plans a dose distribution of a charged particle beam from the shape of the tumor. In this treatment planning apparatus, a stereoscopic image is constructed from cross-sectional images near the tumor, and the dose distribution plan is confirmed by projecting the planned dose distribution of charged particle beams onto the stereoscopic image.

特開平５−３３７２０８号公報JP-A-5-337208

しかしながら、前述した治療計画装置においては、画像データ上で線量分布計画の確認を行っているが、実際に照射された荷電粒子線の線量分布を測定しておらず、確認に対する信頼性が高いとは言えない。荷電粒子線の径路上に線量分布モニタを設けることも考えられるが、この場合、測定できる線量分布は線量分布モニタ位置のものであり、アイソセンター位置における線量分布とは異なる。荷電粒子線の線量分布の確認は治療の信頼性を確保する上で重要であることから、アイソセンター位置における荷電粒子線の線量分布を精度良く確認するための技術が強く求められている。   However, in the treatment planning apparatus described above, the dose distribution plan is confirmed on the image data, but the dose distribution of the actually irradiated charged particle beam is not measured, and the confirmation is highly reliable. I can't say that. Although it is conceivable to provide a dose distribution monitor on the path of the charged particle beam, in this case, the measurable dose distribution is at the dose distribution monitor position and is different from the dose distribution at the isocenter position. Since confirmation of the dose distribution of the charged particle beam is important for ensuring the reliability of treatment, there is a strong demand for a technique for accurately confirming the dose distribution of the charged particle beam at the isocenter position.

そこで、本発明は、アイソセンター位置における荷電粒子線の線量分布を精度良く算出できる荷電粒子線照射装置を提供することを目的とする。 The present invention aims to provide a charged particle beam irradiation equipment which can accurately calculate the dose distribution of the charged particle beam at the isocenter position.

上記課題を解決するため、本発明は、患者に対して荷電粒子線を照射する荷電粒子線照射装置であって、患者が載置される治療台と、治療台のまわりを回転可能な回転ガントリと、回転ガントリに取り付けられ、荷電粒子線を照射する照射部と、荷電粒子線の線量分布を算出する線量分布算出部と、を備え、照射部は、通過する荷電粒子線の線量分布を測定する線量分布モニタと、線量分布モニタよりも荷電粒子線の照射方向における下流側に設けられ、荷電粒子線の照射方向と垂直な平面における荷電粒子線の形状の成形を行うコリメータと、を有し、線量分布算出部は、線量分布モニタの測定結果とコリメータの開口の形状のデータとに基づき、線量分布モニタの測定結果を補正してコリメータよりも荷電粒子線の照射方向における下流側における荷電粒子線の線量分布を算出することを特徴とする。In order to solve the above-described problems, the present invention provides a charged particle beam irradiation apparatus that irradiates a patient with a charged particle beam, a treatment table on which the patient is placed, and a rotating gantry that can rotate around the treatment table. And an irradiation unit that is attached to the rotating gantry and irradiates a charged particle beam, and a dose distribution calculation unit that calculates a dose distribution of the charged particle beam, and the irradiation unit measures a dose distribution of the charged particle beam that passes therethrough. And a collimator provided downstream of the dose distribution monitor in the irradiation direction of the charged particle beam and for shaping the shape of the charged particle beam in a plane perpendicular to the irradiation direction of the charged particle beam. The dose distribution calculation unit corrects the measurement result of the dose distribution monitor based on the measurement result of the dose distribution monitor and the shape data of the opening of the collimator, and is downstream of the collimator in the irradiation direction of the charged particle beam. And calculating a dose distribution of the charged particle beam in.

本発明に係る荷電粒子線照射装置によれば、照射部を通過する荷電粒子線の線量分布の測定結果に基づいて、アイソセンター位置における荷電粒子線の線量分布を算出することができる。従って、この荷電粒子線照射装置によれば、実際に照射することでアイソセンター位置における線量分布を治療前に確認することが可能になる。しかも、この荷電粒子線照射装置では、実際に照射された荷電粒子線の線量分布の測定結果を用いることで、実際の測定結果を用いない場合と比べて、アイソセンター位置における荷電粒子線の線量分布を精度良く算出することができる。   According to the charged particle beam irradiation apparatus according to the present invention, the dose distribution of the charged particle beam at the isocenter position can be calculated based on the measurement result of the dose distribution of the charged particle beam passing through the irradiation unit. Therefore, according to this charged particle beam irradiation apparatus, the dose distribution at the isocenter position can be confirmed before treatment by actually irradiating. Moreover, in this charged particle beam irradiation apparatus, the dose of the charged particle beam at the isocenter position is compared with the case where the actual measurement result is not used by using the measurement result of the dose distribution of the actually irradiated charged particle beam. The distribution can be calculated with high accuracy.

本発明に係る荷電粒子線照射装置において、照射部を通る荷電粒子線の照射方向をＺ軸方向、Ｚ軸に直行する方向をＸ方向とし、コリメータは、Ｘ軸方向で対向する二列の櫛歯を有し、対向する櫛歯の間の開口を調整することで開口を通過する荷電粒子線の成形を行うマルチリーフコリメータであり、線量分布算出部は、線量分布モニタの測定結果とマルチリーフコリメータの対向する櫛歯の間の開口のＸ軸方向における幅とに基づき、線量分布モニタの測定結果を補正してコリメータよりも荷電粒子線の照射方向における下流側における荷電粒子線の線量分布を算出してもよい。In the charged particle beam irradiation apparatus according to the present invention, the irradiation direction of the charged particle beam passing through the irradiation unit is the Z-axis direction, the direction orthogonal to the Z-axis is the X direction, and the collimator is two rows of combs facing each other in the X axis direction This is a multi-leaf collimator that has teeth and adjusts the opening between the opposing comb teeth to form a charged particle beam that passes through the opening. Based on the width in the X-axis direction of the opening between the comb teeth facing each other of the collimator, the measurement result of the dose distribution monitor is corrected to determine the dose distribution of the charged particle beam downstream in the irradiation direction of the charged particle beam from the collimator. It may be calculated.

本発明によれば、アイソセンター位置における荷電粒子線の線量分布を精度良く算出できる。   According to the present invention, the dose distribution of the charged particle beam at the isocenter position can be calculated with high accuracy.

本発明に係る荷電粒子線照射装置の一実施形態を示す斜視図である。1 is a perspective view showing an embodiment of a charged particle beam irradiation apparatus according to the present invention. 図１に示すワブラー方式の照射部を示す斜視図である。It is a perspective view which shows the irradiation part of the wobbler system shown in FIG. 照射部における荷電粒子線の線量分布の変化を示す概略図である。It is the schematic which shows the change of the dose distribution of the charged particle beam in an irradiation part. 線量分布モニタにおける測定結果の一例を示すグラフである。It is a graph which shows an example of the measurement result in a dose distribution monitor. 本発明に係る荷電粒子線の線量分布算出方法を示すフローチャートである。It is a flowchart which shows the dose distribution calculation method of the charged particle beam which concerns on this invention.

以下、本発明の好適な実施形態について、図面を参照して詳細に説明する。   DESCRIPTION OF EMBODIMENTS Hereinafter, preferred embodiments of the present invention will be described in detail with reference to the drawings.

図１に示されるように、照射部１は、荷電粒子線照射装置１００において、治療台１０５を取り囲むように設けられた回転ガントリ１０３に取り付けられ、この回転ガントリ１０３によって治療台１０５の回りに回転可能とされている。そして、治療台１０５に寝かされた患者の腫瘍等の被照射体に対して、荷電粒子線Ｐ（図２参照）を照射する。荷電粒子線Ｐは、電荷をもった粒子を高速に加速したものであり、荷電粒子線Ｐとしては、例えば陽子線、重粒子（重イオン）線等が挙げられる。   As shown in FIG. 1, the irradiation unit 1 is attached to a rotating gantry 103 provided so as to surround the treatment table 105 in the charged particle beam irradiation apparatus 100, and is rotated around the treatment table 105 by the rotation gantry 103. It is possible. Then, a charged particle beam P (see FIG. 2) is irradiated to an irradiated object such as a tumor of a patient who is laid on the treatment table 105. The charged particle beam P is obtained by accelerating charged particles at a high speed. Examples of the charged particle beam P include proton beams and heavy particle (heavy ion) beams.

なお、図１には示されていないが、荷電粒子線照射装置１００は、イオン源で生成した荷電粒子線を加速して荷電粒子線Ｐを出射するサイクロトロンを治療台１０５及び回転ガントリ１０３から離れた位置に備えている。サイクロトロンから出射された荷電粒子線Ｐはビーム輸送系を介して照射部１に供給される。この荷電粒子線照射装置１００は、制御装置１０によって統括的に制御されている。   Although not shown in FIG. 1, the charged particle beam irradiation apparatus 100 separates the cyclotron that accelerates the charged particle beam generated by the ion source and emits the charged particle beam P away from the treatment table 105 and the rotating gantry 103. Prepared in a different position. The charged particle beam P emitted from the cyclotron is supplied to the irradiation unit 1 through the beam transport system. The charged particle beam irradiation apparatus 100 is comprehensively controlled by the control apparatus 10.

図２に示されるように、照射部１は、荷電粒子線Ｐの照射方向に順に配列された四極磁石２、ワブリング磁石３、散乱体４、リッジフィルタ５、線量分布モニタ（線量分布測定手段）６、マルチリーフコリメータ７、及びスノート８を備えている。以下、照射部１を通る荷電粒子線Ｐの照射方向をＺ軸方向、Ｚ軸に直交する二方向をＸ軸方向及びＹ軸方向として説明に用いる。Ｘ軸方向及びＹ軸方向は互いに直交する方向である。   As shown in FIG. 2, the irradiation unit 1 includes a quadrupole magnet 2, a wobbling magnet 3, a scatterer 4, a ridge filter 5 and a dose distribution monitor (dose distribution measuring means) arranged in order in the irradiation direction of the charged particle beam P. 6, a multi-leaf collimator 7 and a snout 8 are provided. Hereinafter, the irradiation direction of the charged particle beam P passing through the irradiation unit 1 will be used as an explanation in the Z axis direction, and two directions orthogonal to the Z axis will be used as an X axis direction and a Y axis direction. The X-axis direction and the Y-axis direction are directions orthogonal to each other.

照射部１の入口に最も近い四極磁石２は、ビーム輸送ラインを介して入力された荷電粒子線Ｐが発散するのを抑えて収束させる電磁石である。四極磁石２の下流には、ワブリング磁石３が配置されている。ワブリング磁石３は、正規分布状の荷電粒子線Ｐを円軌道で走査することで一様な照射野を形成するための電磁石である。   The quadrupole magnet 2 closest to the entrance of the irradiating unit 1 is an electromagnet that converges while suppressing the divergence of the charged particle beam P input through the beam transport line. A wobbling magnet 3 is disposed downstream of the quadrupole magnet 2. The wobbling magnet 3 is an electromagnet for forming a uniform irradiation field by scanning a normally distributed charged particle beam P in a circular orbit.

散乱体４は、ワブリング磁石３で走査された荷電粒子線Ｐの拡散を行う。このような散乱体４は、例えば鉛製の板から構成されている。散乱体４は、荷電粒子線Ｐを照射方向と直交する方向（ＸＹ平面内の方向）に広がりを持つ幅広のビームへと拡散させる。   The scatterer 4 diffuses the charged particle beam P scanned by the wobbling magnet 3. Such a scatterer 4 is composed of, for example, a lead plate. The scatterer 4 diffuses the charged particle beam P into a wide beam having a spread in a direction orthogonal to the irradiation direction (direction in the XY plane).

リッジフィルタ５は、散乱体４で散乱された荷電粒子線Ｐの線量分布の調整を行う。具体的には、リッジフィルタ５は、患者の体内の腫瘍の厚さ（照射方向における腫瘍の長さ）に対応するように、荷電粒子線Ｐに拡大ブラッグピーク（ＳＯＢＰ）を与える。   The ridge filter 5 adjusts the dose distribution of the charged particle beam P scattered by the scatterer 4. Specifically, the ridge filter 5 gives an enlarged Bragg peak (SOBP) to the charged particle beam P so as to correspond to the thickness of the tumor in the patient's body (the length of the tumor in the irradiation direction).

線量分布モニタ６は、通過する荷電粒子線Ｐの線量分布の測定を行う。線量分布モニタ６は、Ｘ軸方向又はＹ軸方向で延在して格子状に配列された複数の金属ワイヤを有しており、荷電粒子線Ｐが照射された金属ワイヤから発生する電子を検出する。これらの複数の金属線にはＸ軸及びＹ軸の位置座標に応じたチャネル番号が付されており、各チャネル番号の金属ワイヤから発生した電子を検出することで、荷電粒子線Ｐの線量分布の測定が行われる。   The dose distribution monitor 6 measures the dose distribution of the charged particle beam P that passes therethrough. The dose distribution monitor 6 has a plurality of metal wires extending in the X-axis direction or the Y-axis direction and arranged in a lattice pattern, and detects electrons generated from the metal wires irradiated with the charged particle beam P. To do. Channel numbers corresponding to the position coordinates of the X axis and the Y axis are assigned to the plurality of metal wires, and the dose distribution of the charged particle beam P is detected by detecting electrons generated from the metal wires of the respective channel numbers. Is measured.

マルチリーフコリメータ７は、照射方向と垂直な平面（ＸＹ平面）における荷電粒子線Ｐの形状の成形を行う。マルチリーフコリメータ７は、Ｘ軸方向で対向する二列の櫛歯を有しており、対向する櫛歯間の開口を調整することで、開口を通過する荷電粒子線Ｐの成形を行う。なお、マルチリーフコリメータ７に代えて、ブロックコリメータを用いても良い。ブロックコリメータとしては、患者の腫瘍形状に合わせて穴を形成された金属柱等が用いられる。   The multi-leaf collimator 7 forms the shape of the charged particle beam P on a plane (XY plane) perpendicular to the irradiation direction. The multi-leaf collimator 7 has two rows of comb teeth facing each other in the X-axis direction, and the charged particle beam P passing through the opening is formed by adjusting the opening between the facing comb teeth. Instead of the multi-leaf collimator 7, a block collimator may be used. As the block collimator, a metal column or the like in which a hole is formed according to the tumor shape of the patient is used.

スノート８は、照射方向（Ｚ軸方向）の奥行きについて荷電粒子線Ｐを腫瘍の形状に合わせて成形するボーラスを保持する部材である。ボーラスは、荷電粒子線Ｐの最大到達深さの部分の立体形状を、腫瘍の最大深さ部分の形状に合わせて成形する。   The snout 8 is a member that holds a bolus for shaping the charged particle beam P in accordance with the shape of the tumor with respect to the depth in the irradiation direction (Z-axis direction). The bolus is formed so that the three-dimensional shape of the portion where the charged particle beam P reaches the maximum depth is matched with the shape of the maximum depth portion of the tumor.

図１及び図２に示されるように、荷電粒子線照射装置１００は、その全体的な制御を行う制御装置１０を備えている。制御装置１０は、患者に対する荷電粒子線Ｐの照射制御やマルチリーフコリメータ７の位置制御を行う。また、制御装置１０は、アイソセンター位置における荷電粒子線Ｐの線量分布を算出するための線量分布算出部（線量分布算出手段）１１を有している。線量分布算出部１１は、線量分布モニタ６の測定結果に基づいて、アイソセンター位置における荷電粒子線Ｐの線量分布を算出する。   As shown in FIGS. 1 and 2, the charged particle beam irradiation apparatus 100 includes a control device 10 that performs overall control thereof. The control device 10 performs irradiation control of the charged particle beam P to the patient and position control of the multi-leaf collimator 7. Further, the control device 10 has a dose distribution calculation unit (dose distribution calculation means) 11 for calculating the dose distribution of the charged particle beam P at the isocenter position. The dose distribution calculation unit 11 calculates the dose distribution of the charged particle beam P at the isocenter position based on the measurement result of the dose distribution monitor 6.

次に、本実施形態に係る荷電粒子線照射装置１００における線量分布算出方法について説明する。   Next, a dose distribution calculation method in the charged particle beam irradiation apparatus 100 according to the present embodiment will be described.

図３に示されるように、本実施形態に係る線量分布算出方法は、線量分布モニタ６が荷電粒子線Ｐの線量分布を測定する線量分布測定ステップＳ１と、線量分布算出部１１がアイソセンター位置における荷電粒子線Ｐの線量分布を算出する線量分布算出ステップＳ２と、を有している。   As shown in FIG. 3, the dose distribution calculation method according to the present embodiment includes a dose distribution measurement step S1 in which the dose distribution monitor 6 measures the dose distribution of the charged particle beam P, and the dose distribution calculation unit 11 in the isocenter position. A dose distribution calculating step S2 for calculating the dose distribution of the charged particle beam P at.

ここで、図４は、ＸＺ平面における荷電粒子線Ｐの線量分布の変化を説明するための概略図である。図４において、四極磁石２、リッジフィルタ５、及びスノート８の描写を省略する。照射部１における中心軸（照射部１を構成する機器２〜８の略中心を通る軸）をＣとして示す。   Here, FIG. 4 is a schematic diagram for explaining a change in the dose distribution of the charged particle beam P in the XZ plane. In FIG. 4, the depiction of the quadrupole magnet 2, the ridge filter 5, and the snout 8 is omitted. A central axis in the irradiation unit 1 (an axis passing through the approximate center of the devices 2 to 8 constituting the irradiation unit 1) is indicated as C.

図４に示されるように、ビーム輸送ラインから入射した荷電粒子線Ｐは、四極磁石２により収束され、その線量分布はＷ１のような山なりの状態となる。この場合の線量分布Ｗ１は、位置座標ｘの関数ｆ（ｘ）として下記の式（１）で表すことができる。式（１）におけるＡは規格化定数である。また、μは実際の線量分布の中心位置座標であり、中心軸Ｃからの距離（Ｘ軸方向における距離）として表すことができる。σは線量分布の広がりである。この関数ｆ（ｘ）はガウシアン（ガウス関数）となる。

Figure 2015163309
As shown in FIG. 4, the charged particle beam P incident from the beam transport line is converged by the quadrupole magnet 2, and the dose distribution is in a mountainous state like W <b> 1. The dose distribution W1 in this case can be expressed by the following equation (1) as a function f (x) of the position coordinate x. A in Formula (1) is a normalization constant. Further, μ is the center position coordinate of the actual dose distribution, and can be expressed as a distance from the center axis C (a distance in the X-axis direction). σ is the spread of the dose distribution. This function f (x) is a Gaussian (Gaussian function).
Figure 2015163309

線量分布Ｗ１の荷電粒子線Ｐがワブリング磁石３によって走査されると、線量分布はＷ２の状態（二つの山からなる状態）となる。この場合のワブリング半径をＲとする。ワブリング半径Ｒは、入射した荷電粒子線Ｐの粒子エネルギーＥ、ワブリング磁石３の磁場強度Ｂ、ワブリング磁石３から線量分布Ｗ２の現在位置までの距離Ｌｗに依存するパラメータである。ここで、線形分布Ｗ１に対するワブリング磁石３の走査の影響を表す関数ｇ（ｘ）を検討すると、下記の式（２）として表すことができる。この場合、線量分布Ｗ２の関数ＤＭ（ｘ）は、下記の式（３）として表すことができる。すなわち、関数ＤＭ（ｘ）は、関数ｆ（ｘ）及び関数ｇ（ｘ）を畳み込み積分することで求めることができる。

Figure 2015163309
Figure 2015163309
 When the charged particle beam P of the dose distribution W1 is scanned by the wobbling magnet 3, the dose distribution is in the state of W2 (a state consisting of two peaks). Let R be the wobbling radius in this case. The wobbling radius R is a parameter that depends on the particle energy E of the incident charged particle beam P, the magnetic field strength B of the wobbling magnet 3, and the distance Lw from the wobbling magnet 3 to the current position of the dose distribution W2. Here, when the function g (x) representing the influence of the scanning of the wobbling magnet 3 on the linear distribution W1 is examined, it can be expressed as the following equation (2). In this case, the function D M (x) of the dose distribution W2 can be expressed as the following formula (3). That is, the function D M (x) can be obtained by convolving and integrating the function f (x) and the function g (x).
Figure 2015163309
Figure 2015163309

その後、荷電粒子線Ｐが散乱体４によってＸ軸方向及びＹ軸方向に拡散されると、線量分布はＷ３の状態となる。線量分布Ｗ３は、線量分布モニタ６の位置（線量分布測定位置）における荷電粒子線Ｐの線量分布である。この線量分布Ｗ３を表す関数は、上記式（２）において、関数ｆ（ｘ）に含まれるσを線量分布モニタ６の位置における広がりσＭに置き換えることで求められる。このσＭは、入射時における線形分布Ｗ１の広がりσ０、散乱体４の厚さｄ、散乱体４から線量分布モニタ６までの距離ＬＳＭに依存するパラメータである。この線量分布Ｗ３が線量分布測定ステップＳ１において線量分布モニタ６に測定される。 Thereafter, when the charged particle beam P is diffused by the scatterer 4 in the X-axis direction and the Y-axis direction, the dose distribution becomes W3. The dose distribution W3 is a dose distribution of the charged particle beam P at the position of the dose distribution monitor 6 (dose distribution measurement position). The function representing the dose distribution W3 is obtained by replacing σ included in the function f (x) with the spread σ M at the position of the dose distribution monitor 6 in the above equation (2). This σ M is a parameter that depends on the spread σ 0 of the linear distribution W 1 at the time of incidence, the thickness d of the scatterer 4, and the distance L SM from the scatterer 4 to the dose distribution monitor 6. This dose distribution W3 is measured by the dose distribution monitor 6 in the dose distribution measurement step S1.

図５は、線量分布モニタ６の測定結果の一例を示すグラフである。横軸は、線量分布モニタ６のＸ軸座標に対応するチャネル番号、縦軸は線量相当値を示している。線量相当値とは、荷電粒子線Ｐの線量に比例する値であり、所定の係数を乗じることで線量が得られる。   FIG. 5 is a graph showing an example of the measurement result of the dose distribution monitor 6. The horizontal axis represents the channel number corresponding to the X-axis coordinate of the dose distribution monitor 6, and the vertical axis represents the dose equivalent value. The dose equivalent value is a value proportional to the dose of the charged particle beam P, and the dose can be obtained by multiplying by a predetermined coefficient.

線量分布算出部１１は、図５に示す測定結果の線量分布データ（折れ線状のデータ）に対してフィッティングを行う。このようなフィッティングには、Levenberg−Marquardt法やNelder−Mead Simplex法等を利用することができる。フィッティング後の線量分布データの一例を図５に一点鎖線として示す。   The dose distribution calculation unit 11 performs fitting on the dose distribution data (string-line data) of the measurement results shown in FIG. For such fitting, the Levenberg-Marquardt method, the Nelder-Mead Simplex method, or the like can be used. An example of the dose distribution data after fitting is shown as a one-dot chain line in FIG.

ここで、線量モニタ６の位置における線量分布の広がりσＭと線量モニタ６の位置におけるワブリング半径ＲＭは、予め計算や実測で決定することができる。従って、上記の式（３）におけるフリーパラメータは、規格化定数Ａと線量モニタ６の位置における線量分布の中心位置座標μＭのみである。Ａは規格化定数なので、物理的な意味はない。このため、フィッティングによってμＭが求められる。線量分布算出部１１は、フィッティングによって求めたμＭを、アイソセンター位置における線量分布Ｗ４の中心位置座標μＩとして用いる。 Here, wobbling radius R M in the spread sigma M and the position of the dose monitor 6 of the dose distribution at the position of the dose monitor 6 can be determined in advance by calculation or actual measurement. Thus, the free parameters in the above equation (3) is only the center position coordinates mu M of dose distribution at the position of the normalization constant A and a dose monitor 6. Since A is a normalization constant, there is no physical meaning. For this reason, the μ M is obtained by fitting. Dose distribution calculation unit 11, a mu M as determined by fitting, used as the center position coordinates mu I dose distribution W4 at isocenter position.

アイソセンター位置における線量分布Ｗ４は、上記の式（３）にマルチリーフコリメータ７の効果を乗じることで求められる。マルチリーフコリメータ７の効果は、下記の式（４），（５）として表わすことができる。下記の式（４），（５）において、ａはマルチリーフコリメータ７の開口の幅ｑ（Ｘ軸方向における幅）等に依存するパラメータである。また、ｓは、入射した荷電粒子線Ｐの粒子エネルギーＥ、ワブリング磁石３の磁場強度Ｂ、マルチリーフコリメータ７からアイソセンター位置までの距離ＬＣＩ等に依存するパラメータである。

Figure 2015163309
The dose distribution W4 at the isocenter position is obtained by multiplying the above expression (3) by the effect of the multi-leaf collimator 7. The effect of the multi-leaf collimator 7 can be expressed as the following equations (4) and (5). In the following formulas (4) and (5), a is a parameter depending on the opening width q (width in the X-axis direction) of the multi-leaf collimator 7 and the like. S is a parameter that depends on the particle energy E of the incident charged particle beam P, the magnetic field strength B of the wobbling magnet 3, the distance LCI from the multi-leaf collimator 7 to the isocenter position, and the like.
Figure 2015163309

アイソセンター位置における荷電粒子線Ｐの線量分布Ｗ４を表す関数ＤＩ（ｘ）は、ｆ（ｘ）とｇ（ｘ）とを畳み込み積分をしたものにｈ（ｘ）を乗じたものとなる。従って、関数ＤＩ（ｘ）を下記の式（６）として表すことができる。式（６）において、アイソセンター位置における線量分布Ｗ４の中心位置座標μＩにはフィッティングで求めたμＭを用いる。また、アイソセンター位置における線量分布Ｗ４の広がりσＩ及びアイソセンター位置におけるワブリング半径ＲＩには計算や実測で決定した値を利用する。

Figure 2015163309
The function D I (x) representing the dose distribution W4 of the charged particle beam P at the isocenter position is obtained by convolving and integrating f (x) and g (x) with h (x). Therefore, the function D I (x) can be expressed as the following equation (6). In the formula (6), the center position coordinates mu I dose distribution W4 at isocenter position using mu M obtained in fitting. Further, values determined by calculation or measurement are used for the spread σ I of the dose distribution W4 at the isocenter position and the wobbling radius R I at the isocenter position.
Figure 2015163309

線量分布算出部１１は、線量分布算出ステップＳ２において上記の式（６）に基づいて、アイソセンター位置における荷電粒子線Ｐの線量分布Ｗ４の関数ＤＩ（ｘ）を算出する。線量分布算出部１１は、線量モニタ６の測定結果のフィッティングデータと上記の式（６）とに基づいて、アイソセンター位置における荷電粒子線Ｐの線量分布Ｗ４を算出することができる。なお、線量分布をＸ軸位置座標の関数とした場合の算出について説明したが、線量分布をＹ軸位置座標の関数とした場合も同様に算出することができる。 The dose distribution calculation unit 11 calculates the function D I (x) of the dose distribution W4 of the charged particle beam P at the isocenter position based on the above formula (6) in the dose distribution calculation step S2. The dose distribution calculation unit 11 can calculate the dose distribution W4 of the charged particle beam P at the isocenter position based on the fitting data of the measurement result of the dose monitor 6 and the above equation (6). Although the calculation when the dose distribution is a function of the X-axis position coordinates has been described, the calculation can be similarly performed when the dose distribution is a function of the Y-axis position coordinates.

以上説明した荷電粒子線照射装置１００及び荷電粒子線の線量分布算出方法によれば、照射部１を通過する荷電粒子線Ｐの線量分布の測定結果に基づいて、アイソセンター位置における荷電粒子線Ｐの線量分布を算出することができる。従って、この荷電粒子線照射装置１００及び線量分布算出方法によれば、実際に照射することでアイソセンター位置における線量分布を治療前に確認することが可能になる。しかも、この荷電粒子線照射装置１００及び線量分布算出方法では、実際に照射部１を通過する荷電粒子線Ｐの線量分布の測定結果を用いるので、線量分布の測定結果を用いない場合と比べて、アイソセンター位置における荷電粒子線Ｐの線量分布を精度良く算出することができる。   According to the charged particle beam irradiation apparatus 100 and the charged particle beam dose distribution calculation method described above, the charged particle beam P at the isocenter position is based on the measurement result of the dose distribution of the charged particle beam P passing through the irradiation unit 1. The dose distribution can be calculated. Therefore, according to the charged particle beam irradiation apparatus 100 and the dose distribution calculation method, the dose distribution at the isocenter position can be confirmed before treatment by actually irradiating. Moreover, in the charged particle beam irradiation apparatus 100 and the dose distribution calculation method, since the measurement result of the dose distribution of the charged particle beam P that actually passes through the irradiation unit 1 is used, the measurement result of the dose distribution is not used. The dose distribution of the charged particle beam P at the isocenter position can be calculated with high accuracy.

また、この荷電粒子線照射装置１００及び線量分布算出方法によれば、上述した式（１）〜（６）を用い、実際の荷電粒子線Ｐの線量分布の測定結果を踏まえて線量分布を算出することで、アイソセンター位置における荷電粒子線Ｐの線量分布を一層精度良く算出することができる。   In addition, according to the charged particle beam irradiation apparatus 100 and the dose distribution calculation method, the dose distribution is calculated based on the measurement result of the actual dose distribution of the charged particle beam P using the above-described equations (1) to (6). Thus, the dose distribution of the charged particle beam P at the isocenter position can be calculated with higher accuracy.

本発明は、上述した実施形態に限定されるものではない。例えば、アイソセンター位置における線量分布の中心位置座標μＩは、フィッティングにより求められたμＭと同じにする必要はない。照射部１の入口側のプロファイルモニタによって当該プロファイルモニタの位置における線量分布の中心位置座標μｐを測定した場合には、μｐを用いてアイソセンター位置における線量分布の中心位置座標μＩを算出することもできる。具体的には、μｐと、プロファイルモニタから線量分布モニタ６までの距離ＬＰＭと、線量分布モニタ６からアイソセンター位置までの距離ＬＭＩと、を用いて下記の式（７）から求めることができる。

Figure 2015163309
The present invention is not limited to the embodiment described above. For example, the center position coordinates mu I of dose distribution in isocenter position need not be the same as mu M obtained by the fitting. When measuring the center position coordinates μp dose distribution at the position of the profile monitored by the inlet side of the profile monitor the irradiation unit 1 is to calculate the center position coordinates mu I of dose distribution in the isocenter position with μp You can also. Specifically, it can be obtained from the following equation (7) using μp, the distance L PM from the profile monitor to the dose distribution monitor 6, and the distance L MI from the dose distribution monitor 6 to the isocenter position. it can.
Figure 2015163309

また、本発明は、ワブラー法による荷電粒子線の線量分布算出に限られない。例えば、ワブリング磁石を使わずに、複数の散乱体により荷電粒子線を拡散させる二重散乱体法においても本発明を有効に適用することができる。また、スキャニング法による荷電粒子線の線量分布算出にも適用することができる。更に、本発明は、上述したような関数を用いた算出方法に限られず、様々な方法によりアイソセンター位置の線量分布の算出を行うことができる。   The present invention is not limited to the calculation of the dose distribution of charged particle beams by the wobbler method. For example, the present invention can be effectively applied to a double scatterer method in which a charged particle beam is diffused by a plurality of scatterers without using a wobbling magnet. It can also be applied to the calculation of the dose distribution of charged particle beams by the scanning method. Furthermore, the present invention is not limited to the calculation method using the function as described above, and the dose distribution at the isocenter position can be calculated by various methods.

１…照射部 ２…四極磁石 ３…ワブリング磁石 ４…散乱体 ５…リッジフィルタ ６…線量分布モニタ ７…マルチリーフコリメータ ８…スノート １０…制御装置 １１…線量分布算出部（線量分布算出手段） １００…荷電粒子線照射装置 １０３…回転ガントリ １０５…治療台 Ｃ…中心軸 Ｐ…荷電粒子線 ｑ…開口幅 Ｒ…ワブリング半径 Ｗ１-Ｗ４…線形分布
DESCRIPTION OF SYMBOLS 1 ... Irradiation part 2 ... Quadrupole magnet 3 ... Wobbling magnet 4 ... Scattering body 5 ... Ridge filter 6 ... Dose distribution monitor 7 ... Multi-leaf collimator 8 ... Snout 10 ... Control apparatus 11 ... Dose distribution calculation part (dose distribution calculation means) 100 ... charged particle beam irradiation device 103 ... rotating gantry 105 ... treatment table C ... central axis P ... charged particle beam q ... aperture width R ... wobbling radius W1-W4 ... linear distribution

Claims (6)

アイソセンター位置に配置された被照射体に対して荷電粒子線を照射する荷電粒子線照射装置であって、
前記被照射体に荷電粒子線を照射する照射部と、
前記照射部に設けられ、荷電粒子線の線量分布を測定する線量分布測定手段と、
前記線量分布測定手段の測定結果に基づいて、前記アイソセンター位置における荷電粒子線の線量分布を算出する線量分布算出手段と、
を備えることを特徴とする荷電粒子線照射装置。 A charged particle beam irradiation apparatus that irradiates an object to be irradiated arranged at an isocenter position with a charged particle beam,
An irradiation unit for irradiating the irradiated body with a charged particle beam;
A dose distribution measuring means for measuring a dose distribution of a charged particle beam provided in the irradiation unit;
A dose distribution calculating means for calculating a dose distribution of the charged particle beam at the isocenter position based on the measurement result of the dose distribution measuring means;
A charged particle beam irradiation apparatus comprising: 前記照射部に設けられ、荷電粒子線を走査する電磁石を更に備え、
前記線量分布算出手段は、前記照射部に入射した荷電粒子線の線量分布を荷電粒子線の照射方向と直交する方向の位置座標の関数とし、当該関数に対して前記電磁石の走査に対応する関数を畳み込み積分して得られる前記線量分布測定手段の位置における荷電粒子線の線量分布の関数と、前記線量分布測定手段の測定結果から得られる線量分布の中心位置座標と、に基づいて、前記アイソセンター位置における荷電粒子線の線量分布を算出することを特徴とする請求項１に記載の荷電粒子線照射装置。 An electromagnet that is provided in the irradiation unit and scans the charged particle beam,
The dose distribution calculating means uses the dose distribution of the charged particle beam incident on the irradiation unit as a function of position coordinates in a direction orthogonal to the irradiation direction of the charged particle beam, and a function corresponding to scanning of the electromagnet with respect to the function Based on the function of the dose distribution of the charged particle beam at the position of the dose distribution measuring means obtained by convolution integration and the center position coordinates of the dose distribution obtained from the measurement result of the dose distribution measuring means. The charged particle beam irradiation apparatus according to claim 1, wherein the dose distribution of the charged particle beam at the center position is calculated. 前記照射方向と直交する方向をＸ軸方向とした場合に、前記照射部に入射した荷電粒子線の線量分布のＸ軸位置座標の関数ｆ（ｘ）が下記の式（１）で表され、前記電磁石の走査に対応する関数ｇ（ｘ）が下記の式（２）で表され、前記線量分布測定手段の位置における荷電粒子線の線量分布の関数ＤＭ（ｘ）が下記の式（３）で表されることを特徴とする請求項２に記載の荷電粒子線照射装置。
Figure 2015163309
Figure 2015163309
Figure 2015163309
 但し、式（１）におけるＡは規格化定数、σは線量分布のＸ軸方向の広がり、μは線量分布のＸ軸上の中心位置座標である。また、式（２）におけるＲはワブリング半径である。 When the direction orthogonal to the irradiation direction is the X-axis direction, the function f (x) of the X-axis position coordinate of the dose distribution of the charged particle beam incident on the irradiation unit is expressed by the following equation (1): The function g (x) corresponding to the scanning of the electromagnet is expressed by the following formula (2), and the function D M (x) of the dose distribution of the charged particle beam at the position of the dose distribution measuring means is expressed by the following formula (3) The charged particle beam irradiation apparatus according to claim 2, wherein
Figure 2015163309
Figure 2015163309
Figure 2015163309
 In Equation (1), A is a normalization constant, σ is the spread of the dose distribution in the X-axis direction, and μ is the center position coordinate on the X-axis of the dose distribution. Moreover, R in Formula (2) is a wobbling radius. アイソセンター位置に配置された被照射体に対して照射される荷電粒子線の線量分布算出方法であって、
前記被照射体に荷電粒子線を照射する照射部に設けられた線量分布モニタによって荷電粒子線の線量分布を測定する線量分布測定ステップと、
前記線量分布測定ステップの測定結果に基づいて、前記アイソセンター位置における荷電粒子線の線量分布を算出する線量分布算出ステップと、
を有することを特徴とする荷電粒子線の線量分布算出方法。 A method for calculating a dose distribution of a charged particle beam irradiated to an irradiation object arranged at an isocenter position,
A dose distribution measuring step of measuring a dose distribution of the charged particle beam by a dose distribution monitor provided in an irradiation unit that irradiates the irradiated body with the charged particle beam; and
A dose distribution calculating step for calculating a dose distribution of the charged particle beam at the isocenter position based on the measurement result of the dose distribution measuring step;
A method for calculating a dose distribution of a charged particle beam, comprising: 前記線量分布測定ステップでは、走査された荷電粒子線の線量分布を測定し、
前記線量分布算出ステップでは、前記照射部に入射した荷電粒子線の線量分布を荷電粒子線の照射方向と直交する方向の位置座標の関数とし、当該関数に対して電磁石の走査に対応する関数を畳み込み積分して得られる線量分布測定位置における荷電粒子線の線量分布の関数と、前記線量分布測定ステップの測定結果から得られる線量分布の中心位置座標と、に基づいて、前記アイソセンター位置における荷電粒子線の線量分布を算出すること
を特徴とする請求項４に記載の荷電粒子線の線量分布算出方法。 In the dose distribution measurement step, the dose distribution of the scanned charged particle beam is measured,
In the dose distribution calculating step, the dose distribution of the charged particle beam incident on the irradiation unit is a function of position coordinates in a direction orthogonal to the irradiation direction of the charged particle beam, and a function corresponding to the scanning of the electromagnet with respect to the function Based on the function of the dose distribution of the charged particle beam at the dose distribution measurement position obtained by convolution integration and the center position coordinate of the dose distribution obtained from the measurement result of the dose distribution measurement step, the charge at the isocenter position The dose distribution calculation method of the charged particle beam according to claim 4, wherein the dose distribution of the particle beam is calculated. 前記照射方向と直交する方向をＸ軸方向とした場合に、前記照射部に入射した荷電粒子線の線量分布のＸ軸位置座標の関数ｆ（ｘ）が下記の式（１）で表され、前記電磁石の走査に対応する関数ｇ（ｘ）が下記の式（２）で表され、前記線量分布測定位置における荷電粒子線の線量分布の関数ＤＭ（ｘ）が下記の式（３）で表されることを特徴とする請求項５に記載の荷電粒子線の線量分布算出方法。
Figure 2015163309
Figure 2015163309
Figure 2015163309
 但し、式（１）におけるＡは規格化定数、σは線量分布のＸ軸方向の広がり、μは線量分布のＸ軸上の中心位置座標である。また、式（２）におけるＲはワブリング半径である。 When the direction orthogonal to the irradiation direction is the X-axis direction, the function f (x) of the X-axis position coordinate of the dose distribution of the charged particle beam incident on the irradiation unit is expressed by the following equation (1): The function g (x) corresponding to the scanning of the electromagnet is expressed by the following formula (2), and the function D M (x) of the dose distribution of the charged particle beam at the dose distribution measurement position is expressed by the following formula (3). The method of calculating a dose distribution of a charged particle beam according to claim 5, wherein:
Figure 2015163309
Figure 2015163309
Figure 2015163309
 In Equation (1), A is a normalization constant, σ is the spread of the dose distribution in the X-axis direction, and μ is the center position coordinate on the X-axis of the dose distribution. Moreover, R in Formula (2) is a wobbling radius.
JP2015121259A 2015-06-16 2015-06-16 Charged particle beam irradiation equipment Active JP6094906B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2015121259A JP6094906B2 (en) 2015-06-16 2015-06-16 Charged particle beam irradiation equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2015121259A JP6094906B2 (en) 2015-06-16 2015-06-16 Charged particle beam irradiation equipment

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP2011227078A Division JP6099183B2 (en) 2011-10-14 2011-10-14 Charged particle beam irradiation equipment

Publications (2)

Publication Number Publication Date
JP2015163309A true JP2015163309A (en) 2015-09-10
JP6094906B2 JP6094906B2 (en) 2017-03-15

Family

ID=54186508

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2015121259A Active JP6094906B2 (en) 2015-06-16 2015-06-16 Charged particle beam irradiation equipment

Country Status (1)

Country Link
JP (1) JP6094906B2 (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS55136069A (en) * 1979-04-11 1980-10-23 Nippon Electric Co Linear particle accelerator for medical treatment
JPH10314323A (en) * 1997-05-20 1998-12-02 Mitsubishi Electric Corp Irradiation method
JPH11169469A (en) * 1997-12-11 1999-06-29 Hitachi Ltd Irradiation of charged particle beam and device therefor
JPH11169471A (en) * 1997-09-25 1999-06-29 Siemens Medical Syst Inc Irradiation of radioactive ray and device therefor
JP2004041292A (en) * 2002-07-09 2004-02-12 Mitsubishi Electric Corp Dose calculation method
JP2009072287A (en) * 2007-09-19 2009-04-09 Hitachi-Ge Nuclear Energy Ltd Particle beam irradiation apparatus

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS55136069A (en) * 1979-04-11 1980-10-23 Nippon Electric Co Linear particle accelerator for medical treatment
JPH10314323A (en) * 1997-05-20 1998-12-02 Mitsubishi Electric Corp Irradiation method
JPH11169471A (en) * 1997-09-25 1999-06-29 Siemens Medical Syst Inc Irradiation of radioactive ray and device therefor
JPH11169469A (en) * 1997-12-11 1999-06-29 Hitachi Ltd Irradiation of charged particle beam and device therefor
JP2004041292A (en) * 2002-07-09 2004-02-12 Mitsubishi Electric Corp Dose calculation method
JP2009072287A (en) * 2007-09-19 2009-04-09 Hitachi-Ge Nuclear Energy Ltd Particle beam irradiation apparatus

Also Published As

Publication number Publication date
JP6094906B2 (en) 2017-03-15

Similar Documents

Publication Publication Date Title
US7920675B2 (en) Producing a radiation treatment plan
US8232536B2 (en) Particle beam irradiation system and method for controlling the particle beam irradiation system
EP2313155B1 (en) Device and method for particle therapy monitoring and verification
EP3006084B1 (en) Particle therapy device and method for setting dose calibration factor
EP2277591B1 (en) Corpuscular radiation treatment apparatus
US10376713B2 (en) Radiation therapy planning system, radiation therapy planning method, and radiation therapy system
US10661101B2 (en) Dose distribution calculation device, particle beam therapy system, and dose distribution calculation method
US20120313002A1 (en) Particle beam treatment device and irradiation dose setting method of the particle beam treatment device
DE19907771A1 (en) Method for checking the radiation control unit of an ion beam therapy system
DE19907774A1 (en) Method for verifying the calculated radiation dose of an ion beam therapy system
DE19907138A1 (en) Method for checking the beam generating means and the beam accelerating means of an ion beam therapy system
DE19907064A1 (en) Procedure for checking an emergency shutdown of an ion beam therapy system
DE19907065A1 (en) Method for checking an isocenter and a patient positioning device of an ion beam therapy system
Parodi et al. The influence of lateral beam profile modifications in scanned proton and carbon ion therapy: a Monte Carlo study
DE19907205A1 (en) Method for operating an ion beam therapy system while monitoring the beam position
EP2186542A1 (en) Method and device for determining a radiation isocenter in a particle therapy room
TWI600451B (en) Particle irradiation treatment device and treatment planning corrective method
JP2018513746A (en) Phantom and method for quality assurance of particle beam therapy system
JP5813207B2 (en) Particle beam scanning irradiation system
EP2442317A1 (en) Particle beam radiation device
JP6099183B2 (en) Charged particle beam irradiation equipment
JP6094906B2 (en) Charged particle beam irradiation equipment
JPH11169469A (en) Irradiation of charged particle beam and device therefor
JP6196912B2 (en) Treatment planning device and program for creating treatment plan information
JP7220403B2 (en) Particle beam therapy system, measurement particle beam CT image generation method, and CT image generation program

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150715

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20150716

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20160624

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20160628

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160826

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20170131

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20170202

R150 Certificate of patent or registration of utility model

Ref document number: 6094906

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150