WO2011121762A1

WO2011121762A1 - 粒子線照射装置及び粒子線治療装置

Info

Publication number: WO2011121762A1
Application number: PCT/JP2010/055863
Authority: WO
Inventors: 高明岩田
Original assignee: 三菱電機株式会社
Priority date: 2010-03-31
Filing date: 2010-03-31
Publication date: 2011-10-06
Also published as: JPWO2011121762A1; TWI535470B; EP2554217A4; US8957393B2; US8829466B2; US20140350324A1; TW201406424A; US20110240875A1; TWI425962B; US8592778B2; EP2554217B1; TW201417854A; CN102686276B; JP4560142B1; CN102686276A; TWI526231B; TW201200195A; US20130345490A1; EP2554217A1

Abstract

　走査電磁石のヒステリシスの影響を排除し、ラスタースキャニングやハイブリッドスキャニングにおいて高精度なビーム照射を実現する粒子線照射装置を得ることを目的とする。　走査電磁石（３）の励磁電流を出力する走査電源（４）と、走査電源（４）を制御する照射制御装置（５）とを備え、照射制御装置（５）は、走査電源（４）に出力された励磁電流の指令値（Ｉ_ｋ）による一連の照射動作であるランスルーの結果を評価し、評価の結果が所定の条件を満たさない場合に励磁電流の指令値（Ｉ_ｋ）を更新してランスルーが実行され、評価の結果が所定の条件を満たした励磁電流の指令値（Ｉ_ｋ）を走査電源（４）に出力する走査電磁石指令値学習生成器（３７）を有した。

Description

粒子線照射装置及び粒子線治療装置

　この発明は、医療用や研究用に用いられる粒子線治療装置に関し、ラスタースキャニングのような走査型の粒子線照射装置及び粒子線治療装置に関する。

　一般に粒子線治療装置は、荷電粒子ビームを発生するビーム発生装置と、ビーム発生装置につながれ、発生した荷電粒子ビームを加速する加速器と、加速器で設定されたエネルギーまで加速された後に出射される荷電粒子ビームを輸送するビーム輸送系と、ビーム輸送系の下流に設置され、荷電粒子ビームを照射対象に照射するための粒子線照射装置とを備える。粒子線照射装置には大きく、荷電粒子ビームを散乱体で散乱拡大し、拡大した荷電粒子ビームを照射対象の形状にあわせて照射野を形成するブロード照射方式と、照射対象の形状に合わせるように、細いペンシル状のビームを走査して照射野形成するスキャニング照射方式（スポットスキャニング、ラスタースキャニング等）とがある。

　スポットスキャニングは、簡単に言えば、点描画のように粒子線ビームを小さなスポット状に照射し、照射野を形成する方法である。すなわち、ビーム供給（点を打つ）、ビーム停止、移動、を繰り返す。照射する線量をスポット位置ごとに変えられる自由度の高い照射方法であり、近年非常に注目をあびている。

　ラスタースキャニングは、簡単に言えば、一筆書きのように粒子線ビームを照射し続け、照射野を形成する方法である。すなわち、ビームを照射し続けながら、目標線量が一定な領域を等速で移動させる方法である。ビームの供給／停止を頻繁に繰り返さなくてよいので、治療時間が短くてすむメリットがある。

　スポットスキャニングとラスタースキャニングの中間の照射方法も提案されている。ラスタースキャニングのようにビームは照射し続け、スポットスキャニングのようにビーム照射位置はスポット位置間を次々と移動していく。スポットスキャニングとラスタースキャニングの両方のメリットを取り入れようとしたものである。本明細書において、この中間の照射方法をハイブリッドスキャニングと呼ぶことにする。

　ブロード照射方式は、コリメータやボーラスを用いて患部形状に合う照射野を形成する。患部形状に合う照射野を形成し、正常組織への不要な照射を防いでおり、最も汎用的に用いられている、優れた照射方式である。しかし、患者ごとにボーラスを製作したり、患部に合わせてコリメータを変形させたりする必要がある。

　一方、スキャニング照射方式は、コリメータやボーラスが不要といった自由度の高い照射方式である。しかし、患部以外の正常組織への照射を防ぐこれら部品を用いないため、ブロード照射方式以上に高いビーム照射位置精度が要求される。

　粒子線治療装置において、照射位置や照射線量の精度を上げる様々な発明がなされている。特許文献１には、正確に患部を照射することができる粒子線治療装置を提供することを目的とし、以下の発明が開示されている。特許文献１の発明は、走査装置による荷電粒子ビームの走査量とその際にビーム位置検出器により検出する荷電粒子ビームのビーム位置とを記憶装置に記憶し、この記憶された走査量及びビーム位置を用い、制御装置により治療計画情報に基づくビーム位置に応じて走査装置の走査量を設定する。実際に照射して得られた走査量とビーム位置との関係が記憶装置に記憶されているため、正確に患部を照射することが期待できるものである。

　特許文献２には、高い安全性を確保し、高精度で荷電粒子ビームを照射できる粒子治療装置を提供することを目的とし、以下の発明が開示されている。特許文献２の発明は、荷電粒子ビーム発生装置から出射された荷電粒子ビームを、ビーム進行方向と垂直な照射面上に走査する走査電磁石に供給し、この走査電磁石を通過した荷電粒子ビームの照射面上における位置及び線量に基づいて、荷電粒子ビーム発生装置からの荷電粒子ビームの出射量を制御する。具体的には、照射面上で分割して形成される複数の領域のうち、目標線量に達した領域への荷電粒子ビームの供給を停止し、目標線量に達していない他の領域に荷電粒子ビームを供給する。このように、各領域における照射線量を目標線量と比較し、荷電粒子ビームの出射量をＯＮ／ＯＦＦ制御（供給／停止）することによって、高い安全性を期待するものである。

　特許文献３には、走査電磁石の電流と磁場との間にあるヒステリシス特性がビーム照射位置の精度を低下させる課題に対して、以下の発明が開示されている。特許文献３の発明は、照射計画に基づくビーム照射位置に対応して、ヒステリシスの影響を考慮しない走査電磁石の電流値を演算する第１の演算手段と、第１の演算手段で演算した走査電磁石の電流値を、ヒステリシスの影響を考慮して補正演算する第２の演算手段とを有し、照射制御装置は、第２の演算手段の演算結果に基づいて走査電磁石の電流を制御する。このように第２の演算手段にヒステリシスの影響を排除するように補正演算を実施することで、すなわち、第２の演算手段にヒステリシス特性を表す数学モデルを有することにより、演算によりビーム照射位置の精度の向上を期待するものである。

特開２００５－２９６１６２号公報特開２００８－２７２１３９号公報特開２００７－１３２９０２号公報

　特許文献１に開示された発明においては、実際に照射をして得られた荷電粒子ビームの走査量とビーム位置との実データに基づいて変換テーブルを作成し、この変換テーブルを用いて走査電磁石の設定電流値を演算している。

　しかしながら、実際には特許文献３に示されているように、走査電磁石の電流と磁場との間にはヒステリシス特性が存在し、電流値が増加しているときと、電流値が減少しているときとでは、異なった磁場となる。すなわち、ある瞬間における走査電磁石の電流値が分かったとしても、その情報からだけでは、磁場の正確な値を特定できない。したがって、特許文献１に開示されている発明では、電磁石のヒステリシスの影響により正確に患部を照射することができない問題点があった。

　特許文献２に開示された発明においては、定義した各領域における照射線量が目標線量となるよう、荷電粒子ビームの出射量をＯＮ／ＯＦＦ制御（供給／停止）している。

　しかしながら、特許文献２に開示された発明に記載された照射面上で分割して形成される複数の領域は、対応する走査電磁石の励磁電流の範囲によって定義される励磁電流空間内の領域（励磁領域）であり、実際の照射空間内の領域（照射領域）とは一致しない。走査電磁石のヒステリシスを考慮しなければ、この励磁領域と照射領域が正確に１対１に対応することはないからである。したがって、このように励磁領域単位で照射線量を管理して安全性を高めようとした装置や方法においても、走査電磁石のヒステリシスの影響を排除しなければ、小領域で線量を管理する効果を発揮できない。すなわち、走査電磁石のヒステリシスによりビーム照射位置の精度が悪くなるという問題点があった。

　特許文献３に開示された発明においては、演算手段の内部にヒステリシスの数学的モデルを作り、演算によって走査電磁石の電流値を補正している。

　しかしながら、ヒステリシスを考慮していても、特許文献３の発明のような考慮の仕方では、まだいくつか問題点がある。第１の問題点は、ヒステリシス特性を演算的な手法を用いて高精度に補正することは、実際にはかなり難しいということである。例えば、電流と磁場とのヒステリシス特性を表す曲線は、入力（電流）の振幅だけではなく、入力（電流）を変化させる速度や、変化させるパターンによって様々な様態となる。この複雑なヒステリシス現象を演算的手法、すなわち数学的モデルで表すことは、多くの分野で長年いろいろと研究工夫がなされてきているが、やはり現実的には難しい。また、第２の問題点は、ビーム照射位置の検出方法にある。従来の多くは、この特許文献３に開示された発明のように、１台もしくは複数のビーム位置モニタのみによってビーム照射位置を検出しようとしている。ビーム位置モニタは、荷電粒子ビームが照射されて、はじめてビーム照射位置が分かるものである。したがって、ビームがターゲットからはずれて正常組織等を照射してしまったとき、単にビームを停止することしかできず、本来照射すべきだった正しい照射位置へビーム照射位置を制御するといったことができない問題点があった。

　この発明は上記のような課題を解決するためになされたものであり、走査電磁石のヒステリシスの影響を排除し、ラスタースキャニングやハイブリッドスキャニングにおいて高精度なビーム照射を実現する粒子線照射装置を得ることを目的とする。

　走査電磁石の励磁電流を出力する走査電源と、走査電源を制御する照射制御装置とを備えた。照射制御装置は、走査電源に出力された励磁電流の指令値による一連の照射動作であるランスルーの結果を評価し、評価の結果が所定の条件を満たさない場合に励磁電流の指令値を更新してランスルーが実行され、評価の結果が所定の条件を満たした励磁電流の指令値を走査電源に出力する走査電磁石指令値学習生成器を有した。

　この発明に係る粒子線照射装置は、走査電磁石指令値学習生成器がランスルーの結果に基づいて、走査電源への励磁電流の指令値を好適なものに学習するので、走査電磁石のヒステリシスの影響を排除し、ラスタースキャニングやハイブリッドスキャニングにおいて高精度なビーム照射を実現することができる。

この発明の実施の形態１における粒子線治療装置の概略構成図である。図１の照射制御装置の構成図である。図１の他の照射制御装置の構成図である。磁場空間で定義した複数領域を示す図ある。学習照射の際の点数表の例を示す図である。この発明の実施の形態２における照射制御装置の構成図である。この発明の実施の形態２における他の照射制御装置の構成図である。この発明の実施の形態３における指令電流を生成する数学的モデルの例である。この発明の実施の形態３における照射制御装置の構成図である。

実施の形態１．
　図１はこの発明の実施の形態１における粒子線治療装置の概略構成図である。粒子線治療装置は、ビーム発生装置５１と、加速器５２と、ビーム輸送装置５３と、粒子線照射装置５４と、治療計画装置５５と、データサーバ５６とを備える。なお、治療計画装置５５を粒子線治療装置の構成部品とは考えず、独立して準備することも可能である。ビーム発生装置５１は、イオン源で発生させた荷電粒子を加速して荷電粒子ビームを発生させる。加速器５２は、ビーム発生装置５１に接続され、発生した荷電粒子ビームを加速する。ビーム輸送装置５３は、加速器５２で設定されたエネルギーまで加速された後に出射される荷電粒子ビームを輸送する。粒子線照射装置５４は、ビーム輸送系５３の下流に設置され、荷電粒子ビームを照射対象１５に照射する。

　治療計画装置５５は、照射対象である患部１５の３次元データに基づいて、照射条件などいくつかの治療の計画を立案し、粒子線治療をシミュレートすることができる。粒子線治療を行う医者が最終的に選択した治療計画は、ここで粒子線治療装置を駆動するためのコードに変換される。駆動するためのコードとは、例えばスポットスキャニングの場合であれば、照射層（レイヤー）Ｚ_ｉ（下添え字ｉはレイヤー番号）ごとのスポットの座標Ｘ_ｊ，Ｙ_ｊ（下添え字ｊはスポット番号）及びそのスポットに照射する目標線量Ｄ_ｊであり、また、ラスタースキャニング等の場合であれば、照射層（レイヤー）Ｚ_ｉごとの照射軌道を、サンプリング周期ごとの照射位置Ｘ_ｋ，Ｙ_ｋ（下添え字ｋはシーケンス番号）で表した時系列データ等である。ここで、スポットに照射する目標線量、としたのは、ブラッグピークの特性から、下層の照射の影響を上層が受けてしまうからであり、その箇所に必要な目標線量ではない。データサーバ５６は、治療計画装置５５で患者毎に生成した治療計画データ及び駆動するためのコードを記憶する。

　粒子線照射装置５４は、ビーム輸送装置５３から入射された入射荷電粒子ビーム１ａを輸送するビーム輸送ダクト２と、入射荷電粒子ビーム１ａに垂直な方向であるＸ方向及びＹ方向に入射荷電粒子ビーム１ａを走査する走査電磁石３ａ、３ｂと、走査電磁石３ａ、３ｂが発生させる磁場を検出する磁場センサ２０ａ、２０ｂと、磁場データ変換器２１と、ビーム位置モニタ７と、位置データ変換器８と、線量モニタ１１と、線量データ変換器１２と、照射制御装置５と、走査電源４とを備える。磁場センサ２０ａ、２０ｂは例えばピックアップコイルを有する磁場センサである。なお、図１に示したように入射荷電粒子ビーム１ａの進行方向はＺ方向である。

　走査電磁石３ａは入射荷電粒子ビーム１ａをＸ方向に走査するＸ方向走査電磁石であり、走査電磁石３ｂは入射荷電粒子ビーム１ａをＹ方向に走査するＹ方向走査電磁石である。磁場センサ２０ａはＸ方向の磁場を検出するＸ方向磁場センサであり、磁場センサ２０ｂはＹ方向の磁場を検出するＹ方向磁場センサである。磁場データ変換器２１は磁場センサ２０ａ、２０ｂで検出した磁場を表すセンサの電気信号を、デジタルデータの測定磁場Ｂｓに変換する。ビーム位置モニタ７は走査電磁石３ａ、３ｂで偏向された出射荷電粒子ビーム１ｂの通過位置を検出する。位置データ変換器８はビーム位置モニタ７で検出した通過位置を表すセンサの電気信号から、照射層（レイヤー）での照射位置を計算し、デジタルデータの測定位置座標Ｐｓを生成する。線量モニタ１１は出射荷電粒子ビーム１ｂの線量を検出する。線量データ変換器１２は、線量モニタ１１で検出した線量を表すセンサの電気信号を、デジタルデータの測定線量Ｄｓに変換する。

　照射制御装置５は、走査電源４への励磁電流の指令値である指令電流Ｉｘ_ｋ，Ｉｙ_ｋを出力し、各照射層（レイヤー）Ｚ_ｉにおける照射位置を制御する。走査電源４は照射制御装置５から出力された指令電流Ｉｘ_ｋ，Ｉｙ_ｋに基づいて、走査電磁石３ａ、３ｂへ実際に駆動する励磁電流を出力する。

　図２は照射制御装置５の構成図である。照射制御装置５は、逆写像生成器３０と、逆写像演算器２２と、指令値出力器２５と、指令評価器３３と、指令更新器３４と、走査電磁石指令値系列生成器３５と、ビーム供給指令出力器２６とを有する。走査電磁石指令値系列生成器３５は指令値記憶装置３６を有する。指令値出力器２５と、指令評価器３３と、指令更新器３４と、走査電磁石指令値系列生成器３５は、走査電磁石指令値学習生成器３７を構成する。走査電磁石指令値学習生成器３７は、走査電源４に出力された指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）による一連の照射動作であるランスルーの結果を評価し、評価の結果が所定の条件を満たさない場合に指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）を更新してランスルーが実行され、評価の結果が所定の条件を満たした指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）を走査電源４に出力する。

　照射制御装置５の動作について説明する。粒子線治療装置の照射には、大きく分けるとキャリブレーション時の試し照射と、治療時の本照射とがある。一般的にキャリブレーション時の試し照射は、いわゆる校正のための照射であり、患者のいない状態で、校正が必要なときのみ行う。Ｘ方向走査電磁石３ａへの制御入力（電流Ｉｘ）とＹ方向走査電磁石３ｂへの制御入力（電流Ｉｙ）をいろいろな値にふって試し照射を行い、そのときのビーム照射位置を測定する。実施の形態１におけるキャリブレーション時の試し照射は、従来と同様に行うが、スポットスキャニング的に荷電粒子ビームの供給及び停止を行い、試し照射時において、ビームの測定位置座標Ｐｓ（ｘｓ，ｙｓ）だけではなく、磁場センサ２０ａ，２０ｂにより測定磁場Ｂｓ（Ｂｘｓ，Ｂｙｓ）をも測定する。試し照射をスポットスキャニング的に行うことで、ビームの測定位置座標Ｐｓ（ｘｓ，ｙｓ）を明確に測定することができる。このときの走査電磁石３の測定磁場Ｂｓ（Ｂｘｓ，Ｂｙｓ）とビームの測定位置座標Ｐｓ（ｘｓ，ｙｓ）との関係は、逆写像生成器３０によって生成された逆写像演算器２２の数式モデルとして実現する。

　あらかじめ、試し照射用に指令電流Ｉ_ｌ＝（Ｉｘ_ｌ，Ｉｙ_ｌ）（下添え字ｌは試し照射のスポット番号を表す）を用意する（ステップＳ００１）。試し照射用の指令電流Ｉ_ｌ＝（Ｉｘ_ｌ，Ｉｙ_ｌ）は、照射されるスポットが、粒子線照射装置５４の想定した照射範囲内で照射されるように用意する。指令値出力器２５は、指令電流Ｉ_ｌ＝（Ｉｘ_ｌ，Ｉｙ_ｌ）を走査電源４に出力する（ステップＳ００２）。走査電源４は、指令電流Ｉ_ｌ＝（Ｉｘ_ｌ，Ｉｙ_ｌ）に従って走査電磁石３を制御する（ステップＳ００３）。

　ビーム供給指令出力器２６は、指令値出力器２５と同期して、走査電磁石３ａ、３ｂの磁場が整定するに十分な整定時間を待って、ビーム発生装置５１にビームの発生を指示するビーム供給指令Ｓ_{ｓｔａｒｔ}を出力する。ビーム発生装置５１は荷電粒子ビームの照射を開始する。試し照射に必要な照射時間Ｔ_ｏｎ経過後、ビーム発生装置５１にビームの停止を指示するビーム停止指令Ｓ_ｓｔｏｐを出力し、ビーム発生装置５１は荷電粒子ビームの照射を停止する。

　指令電流Ｉ_ｌ＝（Ｉｘ_ｌ，Ｉｙ_ｌ）により制御された走査電磁石３ａ、３ｂの磁場を、磁場センサ２０ａ、２０ｂにより測定する。試し照射のスポット毎に測定された測定磁場Ｂ_ｌ＝（Ｂｘ_ｌ，Ｂｙ_ｌ）は、磁場データ変換器２１を介して逆写像生成器３０に入力される。

　ビーム位置モニタ７により、スポット毎に走査電磁石３で走査された出射荷電粒子ビーム１ｂの照射位置座標Ｐ_ｌ＝（Ｘ_ｌ，Ｙ_ｌ）を算出する。照射位置座標Ｐ_ｌ＝（Ｘ_ｌ，Ｙ_ｌ）は、位置データ変換器８を介して逆写像生成期３０に入力される。

逆写像生成器３０は、全スポットの測定磁場Ｂ_ｌ＝（Ｂｘ_ｌ，Ｂｙ_ｌ）及び照射位置座標Ｐ_ｌ＝（Ｘ_ｌ，Ｙ_ｌ）を、内蔵された記憶装置であるメモリに記憶する（ステップＳ００４）。

　逆写像生成器３０は、記憶した測定磁場Ｂ_ｌ＝（Ｂｘ_ｌ，Ｂｙ_ｌ）及び照射位置座標Ｐ_ｌ＝（Ｘ_ｌ，Ｙ_ｌ）に基づいて数式モデルを作成し、作成した数式モデルを逆写像演算器２２に格納する（ステップＳ００５）。

　逆写像演算器２２の数式モデルは、好適な一例として、多項式を用いて実現する。従来の変換テーブルとは異なり、逆写像演算器２２としたことについて説明をする。走査電磁石３の仕様、走査電源４の仕様、及び照射ビームの仕様（照射エネルギー、入射ビーム位置など）が一定であるという仮定のもとでは、走査電磁石３の磁場Ｂ（Ｂｘ，Ｂｙ）が決まれば、ビームの照射位置座標Ｐ（ｘ，ｙ）が一意に決まるので、磁場Ｂとビームの照射位置座標Ｐとの関係に関する物理現象は２入力２出力の正写像と捉えることができる。

　しかし、治療における本照射の際は、ビームの目標照射位置座標Ｐ_ｏｂｊ＝（Ｐｘ_ｏｂｊ，Ｐｙ_ｏｂｊ）が先に与えられて、そのビームの目標照射位置座標Ｐ_ｏｂｊ＝（Ｐｘ_ｏｂｊ，Ｐｙ_ｏｂｊ）を実現するように走査電磁石３の磁場Ｂ（Ｂｘ，Ｂｙ）を制御しなければならない。つまり、治療における本照射ではビームの目標照射位置座標Ｐ_ｏｂｊ＝（Ｐｘ_ｏｂｊ，Ｐｙ_ｏｂｊ）から、その目標照射位置座標Ｐ_ｏｂｊ＝（Ｐｘ_ｏｂｊ，Ｐｙ_ｏｂｊ）を実現するように走査電磁石３の磁場Ｂ（Ｂｘ，Ｂｙ）の推定値Ｂ_ｅｓｔ＝（▲Ｂｘ▼，▲Ｂｙ▼）を算出しなければならない（▲Ｂｘ▼及び▲Ｂｙ▼の説明は、数式（１）及び数式（２）参照）。このように、本照射では位置から磁場への、逆向きの写像が必要であることに注意する。したがって、磁場Ｂ（Ｂｘ，Ｂｙ）の推定値Ｂ_ｅｓｔを得るために逆写像演算器２２が必要となる。

　逆写像演算器２２の数式モデルを多項式にて実現する方法の概略を説明する。ここで多項式とは、数学で一般的に定義されている多項式（polynomial）であり、「定数及び不定元の和と積のみからなる式」などと定義されている。具体的には、例えば以下の数式に示すようなものである。

ここで、ｍ_００、ｍ_０１、ｍ_０２、ｍ_１０、ｍ_１１、ｍ_２０、ｎ_００、ｎ_０１、ｎ_０２、ｎ_１０、ｎ_１１、ｎ_２０は未知パラメータ定数である。また、Ｐｘ_ｏｂｊ、Ｐｙ_ｏｂｊは多項式の不定元に相当する。また、数式（１）の左辺（Ｂの上に小⌒が付いたＢｘ、これを▲Ｂｘ▼で示す。）はＢｘの推定値であり、数式（２）の左辺（Ｂの上に小⌒が付いたＢｙ、これを▲Ｂｙ▼で示す。）はＢｙの推定値であることを示す。磁場Ｂ（Ｂｘ，Ｂｙ）の推定値Ｂ_ｅｓｔは（▲Ｂｘ▼，▲Ｂｙ▼）である。

　多項式の未知パラメータ定数は、試し照射の測定磁場Ｂ_ｌ＝（Ｂｘ_ｌ，Ｂｙ_ｌ）及び照射位置座標Ｐ_ｌ＝（Ｘ_ｌ，Ｙ_ｌ）に基づいて、最小二乗法などにより求める。

　目標照射位置座標Ｐ_ｏｂｊ＝（Ｐｘ_ｏｂｊ，Ｐｙ_ｏｂｊ）を実現する磁場Ｂ＝（Ｂｘ，Ｂｙ）の推定値Ｂ_ｅｓｔは、求めた未知パラメータ定数を代入した数式（１）、数式（２）により求める。

　従来技術においては、キャリブレーションの走査電磁石３への制御入力（指令電流Ｉ_ｌ＝（Ｉｘ_ｌ，Ｉｙ_ｌ））とビームの照射位置座標Ｐ_ｌ＝（Ｘ_ｌ，Ｙ_ｌ）との関係を、変換テーブルとして作成し、走査電磁石指令値生成器６にこの変換テーブルを記憶しておく方式であった。

　すなわち、Ｘ方向走査電磁石３ａへの制御入力（指令電流Ｉｘ_ｌ）はビームの目標照射位置Ｐ_ｏｂｊのｘ座標（Ｐｘ_ｏｂｊ）から、Ｙ方向走査電磁石３ｂへの制御入力（指令電流Ｉｙ_ｌ）はビームの目標照射位置Ｐ_ｏｂｊのｙ座標（Ｐｙ_ｏｂｊ）から、それぞれ独立に求めていた。

　しかし、実際にはＸ方向走査電磁石３ａへの制御入力（指令電流Ｉｘ_ｌ）はビームの照射位置Ｐのｘ座標にもｙ座標にも、また、Ｙ方向走査電磁石３ｂへの制御入力（指令電流Ｉｙ_ｌ）もビーム照射位置Ｐのｘ座標にもｙ座標にも影響を与える、すなわち干渉項があるので、独立に求める変換テーブルによる手法では、照射位置精度は劣化する。

　実施の形態１の粒子線照射装置５４は、目標照射位置座標Ｐ_ｏｂｊ＝（Ｐｘ_ｏｂｊ，Ｐｙ_ｏｂｊ）を実現する磁場Ｂ＝（Ｂｘ，Ｂｙ）の推定値Ｂ_ｅｓｔを求めるにあたり、前記干渉項を考慮した数式モデルを逆写像演算器２２に実現したので、従来とは異なり、出射荷電粒子ビーム１ｂの照射位置精度を向上させることができる。

　次に実施の形態１の粒子線治療装置における本照射について説明する。本照射は、ビーム照射位置及び線量の制御を最適化する学習照射と、患者の照射対象１５にビームを照射する治療照射に分かれる。学習照射は以下の手順により行う。

　或る照射対象１５に対して、治療計画装置５５で作成された治療計画のうち医者が最終的に選択したものが、粒子線治療装置を駆動するためのコードに変換され、照射制御装置５に送られる（ステップＳ１０１）。ここでは、学習照射及び治療照射はラスタースキャニングを想定し、駆動するためのコードが照射層（レイヤー）Ｚ_ｉ（下添え字ｉはレイヤー番号）ごとの照射軌道をサンプリング周期ごとの照射位置Ｐ_ｋ（Ｘ_ｋ，Ｙ_ｋ）（下添え字ｋはシーケンス番号）で表した時系列データであるとして説明する。

　学習照射のための指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）（下添え字ｋはシーケンス番号）は後述する方法により作成する（ステップＳ１０２）。学習照射は、患者のいない状態で、学習照射のための指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）に従って、ランスルーを行う。

　学習照射のための指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）は、その候補としては、あるいは学習の初期値としては、極論すれば何でもよい。ここでは、従来技術による方法で求めたものを初期値とする。学習照射のための指令電流は、^ａＩ_ｋ＝（^ａＩｘ_ｋ，^ａＩｙ_ｋ）（ただし、ａは学習の回数であり、初期値の場合ａ＝０とする。）と表し、学習により指令電流が更新されることを明記する。

　また、走査電磁石指令値系列生成器３５は、指令値記憶装置３６に指令電流^ａＩ_ｋ＝（^ａＩｘ_ｋ，^ａＩｙ_ｋ）を記憶させる。指令電流^ａＩ_ｋ＝（^ａＩｘ_ｋ，^ａＩｙ_ｋ）は、下添え字ｋがシーケンス番号を表す、サンプリング周期ごとの制御入力である。患者のいない状態で初期値の指令電流^０Ｉ_ｋ＝（^０Ｉｘ_ｋ，^０Ｉｙ_ｋ）のランスルーを行う（ステップＳ１０３）。

　粒子線治療装置の操作者の、学習照射開始指示により、信号生成器２９が生成した本照射開始信号Ｓｔがビーム供給指令出力器２６及び走査電磁石指令値系列生成器３５に送られる。走査電磁石指令値系列生成器３５は、学習０回目の指令電流^０Ｉ_ｋ＝（^０Ｉｘ_ｋ，^０Ｉｙ_ｋ）を、サンプリング周期毎に、シーケンス番号順に出力する。

　ビーム供給指令出力器２６は、本照射開始信号Ｓｔを受けて、ビーム発生装置５１にビームを発生するビーム供給指令Ｓ_{ｓｔａｒｔ}を出力する。ビーム発生装置５１は荷電粒子ビームの照射を開始する。

　ビーム位置モニタ７により出射荷電粒子ビーム１ｂの通過位置は検出され、位置データ変換器８により計算された測定位置座標Ｐｓである^ａＰ_ｋは指令評価器３３に入力される。指令評価器３３は目標照射位置である照射位置Ｐ_ｋ（Ｘ_ｋ，Ｙ_ｋ）と測定位置座標Ｐｓである^ａＰ_ｋ（^ａＸ_ｋ，^ａＹ_ｋ）を比較し、初期値の指令電流^０Ｉ_ｋ＝（^０Ｉｘ_ｋ，^０Ｉｙ_ｋ）のランスルーの採点を行う（ステップＳ１０４）。ランスルーの採点方法については後述する。

　信号生成器２９は、本照射が終了する時間に本照射終了信号Ｓｅを送信する。本照射が終了する時間とは、本照射が開始された時間からサンプリング周期×ｋ（シーケンス総数）だけ経過した時間である。ビーム供給指令出力器２６は、本照射終了信号Ｓｅを受けて、ビーム発生装置５１にビームを停止するビーム停止指令Ｓ_ｓｔｏｐを出力する。ビーム発生装置５１はビーム停止指令Ｓ_ｓｔｏｐを受けて、荷電粒子ビーム１ａを停止する（ステップＳ１０５）。

　指令更新器３４は、指令評価器３３から出力された評価結果である採点結果に基づいて、学習ａ＝０回目の指令電流^ａＩ_ｋ＝（^ａＩｘ_ｋ，^ａＩｙ_ｋ）の一部のシーケンスを変更の対象に決定し、当該シーケンスを変更する（ステップＳ１０６）。例えば学習ａ＝０回目の指令電流^０Ｉ_ｋ＝（^０Ｉｘ_ｋ，^０Ｉｙ_ｋ）の第３シーケンス^０Ｉ_３＝（^０Ｉｘ_３，^０Ｉｙ_３）を、^０Ｉｘ_３から^０Ｉｘ_３＋ΔＩのように少し変更する。

　次に、患者のいない状態で、一部のシーケンスを少し変更した指令電流を用いて、再度、ランスルーを行う（ステップＳ１０７）。すなわち、ステップＳ１０３乃至ステップＳ１０５を第２のランスルーとして実行する。

　一部のシーケンスを少し変更したことによって、採点結果がＪ点からＪ＋ΔＪ点だけ変わったとする。すると、指令電流の第３のシーケンス^０Ｉｘ_３は、ΔＪ／ΔＩの情報を用いて更新すればよいことがわかる。一般的な学習機能と同様に、ΔＩが正の場合に採点結果が悪ければ、ΔＩを負にする等の更新を行えばよい。この作業は、採点結果が影響を受けるすべてのシーケンスについて実行するとよい。この更新で、学習が１回（ａをインクリメント（ａ＋１）する）とカウントする。

　指令更新器３４は、更新した学習ａ＝１回目の指令電流^ａＩ_ｋ＝（^ａＩｘ_ｋ，^ａＩｙ_ｋ）を作成する（ステップＳ１０８）。第２回のランスルーが前回のランスルーよりも採点結果がよいことを確認し、学習を続ける。学習は、あらかじめ設定した条件（合格点等）を満たすまで繰り返す。最終的に学習が終了した指令電流^ａＩ_ｋ＝（^ａＩｘ_ｋ，^ａＩｙ_ｋ）は、指令値記憶装置３６に記憶する。また、学習の結果、直前のランスルーよりも採点結果が悪くなった場合は、更新の速度（１回で更新する量）を少なくする、あるいは上記で説明したように指令電流を更新する等の工夫をする。

　治療照射は以下の手順により行う。粒子線治療装置の操作者の治療照射開始指示により、本照射開始信号Ｓｔがビーム供給指令出力器２６及び走査電磁石指令値系列生成器３５に送られる。査電磁石指令値系列生成器３５は、学習が終了した指令電流^ａＩ_ｋ＝（^ａＩｘ_ｋ，^ａＩｙ_ｋ）を、サンプリング周期毎に、シーケンス番号順に出力する。

　ビーム供給指令出力器２６は、本照射開始信号Ｓｔを受けて、ビーム発生装置５１にビームを発生するビーム供給指令Ｓ_{ｓｔａｒｔ}を出力する。ビーム発生装置５１は荷電粒子ビームの照射を開始する（ステップＳ１０９）。

　信号生成器２９は、本照射が終了する時間に本照射終了信号Ｓｅを送信する。本照射が終了する時間とは、本照射が開始された時間からサンプリング周期×ｋ（シーケンス総数）だけ経過した時間である。ビーム供給指令出力器２６は、本照射終了信号Ｓｅを受けて、ビーム発生装置５１にビームを停止するビーム停止指令Ｓ_ｓｔｏｐを出力する。ビーム発生装置５１はビーム停止指令Ｓ_ｓｔｏｐを受けて、荷電粒子ビーム１ａを停止する（ステップＳ１１０）。

　次にランスルーの採点方法について説明する。最も直接的な採点方法（第１の採点方法）は、ラスタースキャニング用の駆動するためのコード、すなわち、照射層（レイヤー）Ｚ_ｉごとの照射軌道を、サンプリング周期ごとの照射位置Ｐ_ｋ（Ｘ_ｋ，Ｙ_ｋ）（下添え字ｋはシーケンス番号）で表した時系列データと、実際にランスルーをしたときのサンプリング周期ごとの照射位置^ａＰ_ｋ（^ａＸ_ｋ，^ａＹ_ｋ）（ａは学習の回数）を比較して、以下の評価関数を考えることである。この評価関数の値が、所定の値になった場合（所定の条件を満たした場合）に学習を終了する。

　また、ランスルーの採点は、照射線量に着目して、以下の方法（第２の採点方法）で行うこともできる。ランスルーの採点は、図４に示すように磁場空間内で定義した複数の小領域毎に、目標線量Ｄｉと測定線量Ｄｓとを比較し、図５に示す点数表Ｔにしたがって採点を行う。ランスルーの評価は、例えば、磁場空間内で定義した複数の小領域毎の点数を加えたものを評価関数として定義し、この評価関数の点数により評価を行う。評価関数の点数が高くなったランスルーは点数が低いものより優れていると判定する。図４は磁場空間（Ｂｘ，Ｂｙ）で定義した磁場小領域Ｓ_ｉ，ｊを示す図であり、図５は学習照射の際の点数表Ｔの例を示す図である。なお、小領域に対応した目標線量Ｄｉは、治療計画装置により計算され与えられる。測定線量Ｄｓは、ビーム位置モニタ７の測定結果と、荷電粒子ビーム１ｂが前記小領域を通過した時間等から求める。

　なお、ランスルーの第２の採点方法においては、図３に示すように、指令評価器３３に入力される信号が異なる。これに伴いステップＳ１０４も異なる。図３はランスルーの第２の採点方法を採用した照射制御装置の構成図である。照射制御装置５ｂの走査電磁石指令値学習生成器３７ｂの指令評価器３３ｂには目標線量Ｄｉ及び測定線量Ｄｓ（^ａＤ_ｋ）が入力される。ステップＳ１０４において、線量モニタ１１により線量が検出され、線量データ変換器１２により変換された測定線量Ｄｓは指令評価器３３ｂに入力される。指令評価器３３ｂは目標線量Ｄｉと測定線量Ｄｓを比較し、初期値の指令電流^０Ｉ_ｋ＝（^０Ｉｘ_ｋ，^０Ｉｙ_ｋ）のランスルーの採点を行う。

　図４において、表の左列の（Ｂ_０，Ｂ_１）は、磁場ＢのＸ成分ＢｘがＢ_０≦Ｂｘ＜Ｂ_１の関係を満たすことを簡略して示したものであり、同様に（Ｂ_ｍ－１，Ｂ_ｍ）は、ＢｘがＢ_ｍ－１≦Ｂｘ＜Ｂ_ｍの関係を満たすことを簡略して示したものである。表の上段の（Ｂ_０，Ｂ_１）は、磁場ＢのＹ成分ＢｙがＢ_０≦Ｂｙ＜Ｂ_１の関係を満たすことを簡略して示したものであり、同様に（Ｂ_ｍ－１，Ｂ_ｍ）は、ＢｙがＢ_ｍ－１≦Ｂｙ＜Ｂ_ｍの関係を満たすことを簡略して示したものである。領域Ｓ_０，０はＢ_０≦Ｂｘ＜Ｂ_１及びＢ_０≦Ｂｙ＜Ｂ_１の関係を満たす領域であり、領域Ｓ_{ｍ－１，ｍ－１}はＢ_ｍ－１≦Ｂｘ＜Ｂ_ｍ及びＢ_ｍ－１≦Ｂｙ＜Ｂ_ｍの関係を満たす領域である。

　ランスルーの採点は、粒子線照射装置５４の照射可能範囲に対応する磁場空間の全領域で行う。これにより患者の照射対象１５は目標線量Ｄｉで管理され、照射対象１５ではない正常組織相当箇所は線量が零で管理されるので、照射対象１５及び非照射対象の荷電粒子ビームの線量管理を精度良く行うことができる。

　走査電磁石３の磁場が決まれば、荷電粒子ビームの照射位置は唯一に決まる。言い換えれば、走査電磁石３の磁場と、荷電粒子ビームの照射位置は１対１の関係にある。したがって、従来技術の指令値電流空間で定義された領域のような走査電磁石のヒステリシスの影響を受けることがない。また、キャリブレーションのときに荷電粒子ビームを照射して得られた磁場センサ２０による磁場とビーム位置モニタ７によるビーム位置との関係は、キャリブレーションのときと同様に荷電粒子を走査する場合の本照射の際の磁場とビーム位置との関係とは、極めてよく一致する。実照射空間は、荷電粒子ビームの出射位置、ビーム位置モニタ７における通過位置、粒子線照射装置５４と照射対象１５との位置関係により求めることができるので、実照射空間の領域は磁場空間で定義された領域と写像関係を有し、本照射の際にもこの写像関係はほとんど変わらない。したがって、ランスルーの採点を磁場空間内で定義した複数の磁場小領域Ｓ_ｉ，ｊ毎に行い、この磁場小領域Ｓ_ｉ，ｊ毎に荷電粒子ビームの線量管理を行うので、照射対象１５における実照射空間での線量管理を精度よく行うことができる。

　図５に示した学習照射の際の点数表Ｔは減点方式の例である。線量誤差ｄｅは領域測定線量Ｄｓｓから目標線量Ｄｉを引いた差である。領域測定線量Ｄｓｓは磁場空間（Ｂｘ，Ｂｙ）で定義した磁場小領域Ｓ_ｉ，ｊの実照射線量であり、磁場センサ２０で測定した測定磁場Ｂｓと、線量モニタ１１で測定した測定線量Ｄｓにより作成する。Δｄは線量誤差の幅であり、許容範囲である所定の値に設定される。測定線量Ｄｓが目標線量Ｄｉを超過した場合における点数の変化率の絶対値は、測定線量Ｄｓが目標線量Ｄｉに満たない場合における点数の変化率の絶対値よりも大きくしている。これにより測定線量Ｄｓが目標線量Ｄｉを超過した場合の修正を迅速に的確に行うことができる。

　実施の形態１の第１の採点方法を適用した粒子線照射装置５４は、実際の照射位置を評価関数としているため、走査電磁石３への指令電流^ａＩ_ｋ＝（^ａＩｘ_ｋ，^ａＩｙ_ｋ）を好適なものに学習することができる。また、実施の形態１の第２の採点方法を適用した粒子線照射装置５４は、磁場空間（Ｂｘ，Ｂｙ）で定義した磁場小領域Ｓ_ｉ，ｊ毎に目標線量Ｄｉと領域測定線量Ｄｓｓから点数をつけ、この点数により評価関数を定義したことにより、走査電磁石３への指令電流^ａＩ_ｋ＝（^ａＩｘ_ｋ，^ａＩｙ_ｋ）を好適なものに学習することができる。したがって、高精度で安全性の高い粒子線照射装置を提供することができる。

　粒子線照射装置５４は、第１の採点方法で、走査電磁石３の電流と磁場との間に生じるヒステリシスの影響が反映された実際の照射位置を評価関数として走査電磁石３への指令電流^ａＩ_ｋ＝（^ａＩｘ_ｋ，^ａＩｙ_ｋ）を好適なものに学習しているため、走査電磁石３の電流と磁場との間に生じるヒステリシスの影響を排除することができる。また、第２の採点方法で、粒子線照射装置５４は、磁場空間（Ｂｘ，Ｂｙ）で定義した磁場小領域Ｓ_ｉ，ｊ毎に荷電粒子ビームの線量管理を行うので、走査電磁石３の電流と磁場との間に生じるヒステリシスの影響を排除することができる。したがって走査電磁石のヒステリシスの影響を排除し、高精度なビーム照射を実現することができる。

　磁場センサ２０は、ホール素子を有する磁場センサでも構わない。ホール素子を用いることで走査電磁石３により発生された磁場の絶対値を測定でき、ピックアップコイルで計測された電圧を積分等の演算を行う必要がなくなる。したがって磁場データ変換器２１を簡略化、小型化することができる。

　磁場センサ２０として最も望ましいのは、ビックアップコイルとホール素子の両方を備えるものである。ホール素子の、磁場の絶対値が測定できるという利点と、ピックアップコイルの、磁場の変化分をヒステリシスがなく測定できるという利点との、両方を取り入れることができるからである。

　従来の粒子線照射装置は、１台もしくは複数のビーム位置モニタのみによってビーム照射位置を検出して、測定位置座標により荷電粒子ビームのフィードバック制御をしていた。位置モニタなどの荷電粒子ビームを遮るものを多く配置することは、ビームが散乱拡大することにつながり、所望のビームスポット径が得られないといった問題点があった。

　実施の形態１の粒子線照射装置５４は、治療照射の際には走査電磁石指令値系列生成器３５は学習が終了した指令電流^ａＩ_ｋ＝（^ａＩｘ_ｋ，^ａＩｙ_ｋ）を、サンプリング周期毎に、シーケンス番号順に出力して荷電粒子ビームの照射位置及び照射線量の制御を行ったので、本照射の際にビーム位置モニタ７を図示しない移動装置により移動させ、出射荷電粒子ビーム１ｂがビーム位置モニタ７を通過しないようにしてもよい。このようにすることでビーム位置モニタ７によって出射荷電粒子ビーム１ｂが散乱し拡大されることを防止できる。これにより、ビームスポット径を小さくできる。したがって、小さなビーム径で照射する方がよい場合には、適したスポット径で治療を行うことができる。

　また、指令値出力器２５が、走査電磁石３の磁場Ｂ（Ｂｘ，Ｂｙ）の推定値Ｂ_ｅｓｔに基づいて指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）を生成し、これを初期値として実行したランスルーにより学習することもできる。

　なお、照射制御装置５に逆写像生成器３０、逆写像演算器２２を有する例で説明したが、逆写像生成器３０、逆写像演算器２２がない場合であっても、走査電磁石のヒステリシスの影響を排除し、ラスタースキャニングやハイブリッドスキャニングにおいて高精度なビーム照射を実現することができることは言うまでもない。

　以上のように実施の形態１の粒子線照射装置５４によれば、走査電磁石３の励磁電流を出力する走査電源４と、走査電源４を制御する照射制御装置５とを備え、照射制御装置５は、走査電源４に出力された励磁電流の指令値Ｉ_ｋによる一連の照射動作であるランスルーの結果を評価し、評価の結果が所定の条件を満たさない場合に励磁電流の指令値Ｉ_ｋを更新してランスルーが実行され、評価の結果が所定の条件を満たした励磁電流の指令値Ｉ_ｋを走査電源４に出力する走査電磁石指令値学習生成器３７を有したので、走査電磁石指令値学習生成器３７がランスルーの結果に基づいて、走査電源４への励磁電流の指令値Ｉ_ｋを好適なものに学習でき、走査電磁石のヒステリシスの影響を排除し、ラスタースキャニングやハイブリッドスキャニングにおいて高精度なビーム照射を実現することができる。

　実施の形態１の粒子線治療装置によれば、荷電粒子ビームを発生させるビーム発生装置５１と、ビーム発生装置５１で発生された荷電粒子ビームを加速する加速器５２と、加速器５２により加速された荷電粒子ビームを輸送するビーム輸送装置５３と、ビーム輸送装置５３で輸送された荷電粒子ビームを走査電磁石３で走査して照射対象１５に照射する粒子線照射装置５４とを備え、粒子線照射装置５４は、走査電磁石３の励磁電流を出力する走査電源４と、走査電源４を制御する照射制御装置５とを備え、照射制御装置５は、走査電源４に出力された励磁電流の指令値Ｉ_ｋによる一連の照射動作であるランスルーの結果を評価し、評価の結果が所定の条件を満たさない場合に励磁電流の指令値Ｉ_ｋを更新してランスルーが実行され、評価の結果が所定の条件を満たした励磁電流の指令値Ｉ_ｋを走査電源４に出力する走査電磁石指令値学習生成器３７を有したので、走査電磁石指令値学習生成器３７がランスルーの結果に基づいて、走査電源４への励磁電流の指令値Ｉ_ｋを好適なものに学習でき、走査電磁石のヒステリシスの影響が排除され、ラスタースキャニングやハイブリッドスキャニングにおいて高精度なビーム照射を使用して高精度な粒子線治療を実現することができる。

実施の形態２．
　実施の形態１では、本照射の学習照射において荷電粒子ビームの照射を行ったが、荷電粒子ビームの照射を行わずに学習し、指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）を最適化することができる。以下に説明する。実施の形態２の粒子線治療装置では、学習照射のランスルーの採点として２種類のランスルーの採点方法（第３の採点方法、第４の採点方法）を用いることができる。第３の採点方法は、ラスタースキャニング用の駆動するためのコード、すなわち、照射層（レイヤー）Ｚ_ｉごとの照射軌道を、サンプリング周期ごとの照射位置Ｐ_ｋ（Ｘ_ｋ，Ｙ_ｋ）（下添え字ｋはシーケンス番号）に対応した磁場の推定値Ｂ_ｅｓｔであるＢ_ｋ＝（Ｂｘ_ｋ，Ｂｙ_ｋ）で表した時系列データと、実際にランスルーをしたときのサンプリング周期ごとの測定磁場^ａＢ_ｋ（^ａＢｘ_ｋ，^ａＢｙ_ｋ）（ａは学習の回数）を比較して、以下の評価関数を考えることである。なお、磁場の推定値Ｂ_ｋは逆写像演算器２２で演算したものである。

　第４の採点方法は、図５に示す点数表Ｔにしたがって採点を行う方法である。しかし、学習照射の際に荷電粒子ビームの照射が行われないので、線量誤差ｄｅは、実施の形態１とは異なり、領域線量計算値Ｄｓｃから目標線量Ｄｉを引いた差とする。領域線量計算値Ｄｓｃは磁場センサ２０で測定した測定磁場^ａＢ_ｋ（^ａＢｘ_ｋ，^ａＢｙ_ｋ）から、磁場空間（Ｂｘ，Ｂｙ）で定義した磁場小領域Ｓ_ｉ，ｊ毎に出射荷電粒子ビーム１ｂが滞在する時間分で線量を積分して求める。

　図６はこの発明の実施の形態２における照射制御装置の構成図であり、学習照射のランスルーの採点として第３の採点方法を採用したものである。実施の形態１の照射制御装置とは指令評価器３３に入力される信号が異なる。照射制御装置５ｃの走査電磁石指令値学習生成器３７ｃの指令評価器３３ｃには磁場の推定値Ｂ_ｋ＝（Ｂｘ_ｋ，Ｂｙ_ｋ）及び測定磁場^ａＢ_ｋ（^ａＢｘ_ｋ，^ａＢｙ_ｋ）が入力される。実施の形態２における照射制御装置５ｃの動作について説明する。キャリブレーション時の試し照射は実施の形態１のステップＳ００１乃至ステップＳ００５と同様である。学習照射は、基本的に実施の形態１のステップＳ１０１乃至ステップＳ１０８と同様であるが、荷電粒子ビームの照射を行わないので、ステップＳ１０４が異なる。実施の形態１のステップＳ１０４の代わりにステップＳ２０１を実行する。

　指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）により制御された走査電磁石３ａ、３ｂが発生する磁場を、磁場センサ２０ａ、２０ｂにより測定する。磁場センサ２０により測定され、磁場データ変換器２１で変換された測定磁場^ａＢ_ｋ（^ａＢｘ_ｋ，^ａＢｙ_ｋ）は指令評価器３３ｃに入力される。指令評価器３３ｃは、シーケンス番号毎に逆写像演算器２２で演算された磁場の推定値Ｂ_ｋ＝（Ｂｘ_ｋ，Ｂｙ_ｋ）と測定磁場^ａＢ_ｋ（^ａＢｘ_ｋ，^ａＢｙ_ｋ）を比較し、初期値の指令電流^０Ｉ_ｋ＝（^０Ｉｘ_ｋ，^０Ｉｙ_ｋ）のランスルーの採点を行う（ステップＳ２０１）。

　次に、学習照射のランスルーの採点として第４の採点方法を採用した場合の照射制御装置及び動作を説明する。図７はこの発明の実施の形態２における他の照射制御装置の構成図であり、学習照射のランスルーの採点として第４の採点方法を採用したものである。実施の形態１の照射制御装置とは指令評価器３３に入力される信号が異なる。照射制御装置５ｄの走査電磁石指令値学習生成器３７ｄの指令評価器３３ｄには目標線量Ｄｉ及び領域線量計算値Ｄｓｃが入力される。実施の形態２における照射制御装置５ｄの動作について説明する。キャリブレーション時の試し照射は実施の形態１のステップＳ００１乃至ステップＳ００５と同様である。学習照射は、基本的に実施の形態１のステップＳ１０１乃至ステップＳ１０８と同様であるが、荷電粒子ビームの照射を行わないので、ステップＳ１０４が異なる。実施の形態１のステップＳ１０４の代わりにステップＳ２０２を実行する。

　指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）により制御される荷電粒子ビーム１が磁場小領域Ｓ_ｉ，ｊ毎に滞在する時間分で積分され、領域線量計算値Ｄｓｃが計算され、指令評価器３３ｄに入力される。指令評価器３３ｃは、目標線量Ｄｉと領域線量計算値Ｄｓｃを比較し、初期値の指令電流^０Ｉ_ｋ＝（^０Ｉｘ_ｋ，^０Ｉｙ_ｋ）のランスルーの採点を行う（ステップＳ２０２）。

　実施の形態２の第３の採点方法を適用した粒子線照射装置５４は、測定磁場を評価関数としているため、荷電粒子ビームを照射しなくても走査電磁石３への指令電流^ａＩ_ｋ＝（^ａＩｘ_ｋ，^ａＩｙ_ｋ）を好適なものに学習することができる。また、実施の形態２の第４の採点方法を適用した粒子線照射装置５４は、磁場空間（Ｂｘ，Ｂｙ）で定義した磁場小領域Ｓ_ｉ，ｊ毎に計算された領域線量計算値Ｄｓｃと目標線量Ｄｉから点数をつけ、この点数により評価関数を定義したことにより、荷電粒子ビームを照射しなくても走査電磁石３への指令電流^ａＩ_ｋ＝（^ａＩｘ_ｋ，^ａＩｙ_ｋ）を好適なものに学習することができる。したがって、高精度で安全性の高い粒子線照射装置を提供することができる。

　粒子線照射装置５４は、第３の採点方法で、測定磁場を評価関数として走査電磁石３への指令電流^ａＩ_ｋ＝（^ａＩｘ_ｋ，^ａＩｙ_ｋ）を好適なものに学習しているため、走査電磁石３の電流と磁場との間に生じるヒステリシスの影響を排除することができる。また、第４の採点方法で、粒子線照射装置５４は、磁場空間（Ｂｘ，Ｂｙ）で定義した磁場小領域Ｓ_ｉ，ｊ毎に荷電粒子ビームの線量管理を行うので、走査電磁石３の電流と磁場との間に生じるヒステリシスの影響を排除することができる。したがって走査電磁石のヒステリシスの影響を排除し、高精度なビーム照射を実現することができる。

　また、学習照射のランスルーは荷電粒子ビームを照射しないで行うことができるので、無用なエネルギー消費を抑えることができる。

実施の形態３．
　学習機能とは、「ひとつの課題に対して、より理想的な解に近づく機能」というふうに解釈することができる。実施の形態１及び２では、治療計画装置５５で患者毎に生成したひとつの治療計画データに、いかに実際の照射を近づけるかという課題に対して、より理想的な照射を実現する指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）を生成する機能について説明した。

　制御工学におけるよりアドバンストな学習機能（或いはよりインテリジェントな学習機能）とは、「過去の経験の蓄積により、未来の未知な課題に対して、より理想的な解を導けるための機能」と解釈される。そこで、実施の形態３においては、実施の形態１や２で説明した学習機能をさらに応用し、よりアドバンストな学習機能を備えた粒子線照射装置及び粒子線治療装置について説明する。

　実施の形態１において、学習照射は、ランスルーの採点を行い、採点結果に基づいて、指令電流Ｉ_ｋ自身をより最適なものへと更新していった。この方法は、ひとつの治療計画データに対して、より理想的な照射を実現する指令電流Ｉ_ｋは生成できるが、ここでの経験は、他の治療計画データへは反映することができない。

　実施の形態３においては、学習照射は実施の形態１若しくは２と同様にランスルーの採点を行う。しかし、ランスルーの採点結果に基づいて、指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）自身を更新するのではなく、指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）を生成する数学的モデルを備えるようにし、その数学的モデルのパラメータをより最適なものへと更新するようにする。以下は、その詳細を説明したものである。

　図８は、実施の形態３における、前記指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）を生成する数学的モデルの例を表した図である。図８に基づいて、指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）を生成する仕組みを説明する。

　単層の隠れ層をもつフィードフォワード型ＮＮ（ニューラルネットワーク、以下「ＮＮ」。）６０は、前記指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）を生成する数学的モデルの一例である。入力層６１は、ＮＮ６０の入力部であり、本発明の実施の形態３では目標であるサンプリング周期ごとの照射位置Ｐ_ｋ（Ｘ_ｋ，Ｙ_ｋ）が該当する。入力層６１ａにＸ_ｋが入力され、入力層６１ｂにＹ_ｋが入力される。出力層６３は、ＮＮ６０の出力部であり、本発明の実施の形態３では指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）が該当する。出力層６３ａにＩｘ_ｋが出力され、出力層６３ｂにＩｙ_ｋが出力される。隠れ層６２は、ＮＮ６０の基底関数（活性化関数）であり、複数の隠れ層６２ａ～６２ｎにより、入力層６１ａ、６１ｂからの入力信号に対して重みを与え、出力層６３ａ、６３ｂへと出力する。

　粒子線治療装置においては、走査電磁石３により荷電粒子ビーム１が偏向され、その照射位置が決定される。すなわち、走査電磁石３への一連の指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）が決まれば、ビームの照射位置座標Ｐ＝（Ｘ，Ｙ）が一意に決まるので、一連の指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）からビームの照射位置座標Ｐへの物理現象は、ヒステリシス特性というダイナミクスを有する２入力２出力のシステムと捉えることができる。ＮＮ６０で実現しようとしているのは、まさにこのヒステリシス特性のダイナミクスを有する２入力２出力の逆システムを、数学的に実現しようとするものである。

　実施の形態１及び２では、学習照射における更新は、指令電流Ｉ_ｋ＝（Ｉｘ_ｋ，Ｉｙ_ｋ）自身であった。それに対し、実施の形態３では、ＮＮ６０の隠れ層６２における重みを更新するようにした。図９にＮＮ６０を有する走査電磁石指令値系列生成器３５ｂを備えた照射制御装置５（５ｅ）を示す。照射制御装置５ｅは、実施の形態１の照射制御装置５ａから、逆写像生成器３０、逆写像演算器２２、指令値出力器２５、指令更新器３４が削除され、走査電磁石指令値学習生成器３７ｅにパラメータ更新器３８を有し、走査電磁石指令値系列生成器３５にＮＮ６０が搭載された例である。照射制御装置５ｅは、学習照射のランスルーの採点方法を第１の採点方法とした場合である。パラメータ更新器３８は指令評価器３３の評価結果が所定の条件になるようにＮＮ６０のパラメータである隠れ層６２における重みを更新する。ＮＮ６０を有する走査電磁石指令値系列生成器３５ｂを備えた照射制御装置５は、学習照射のランスルーの採点方法を第２乃至４の採点方法とすることもできる。この場合は、図９の指令評価器３３ａに代えて、それぞれ指令評価器３３ｂ乃至３３ｄにすればよい。なお、第３の採点方法の場合、照射制御装置５は逆写像生成器３０、逆写像演算器２２を備え、磁場の推定値Ｂ_ｋを逆写像演算器２２で演算したものとすればよい。

　以上のように構成することにより、本発明の実施の形態３による粒子線治療装置においては、過去に行った治療計画データに対する学習照射のランスルー（学習ランスルー）の経験が蓄積される。多くの学習ランスルーの経験が蓄積された場合は、あらたな治療計画データに対し、毎回学習ランスルーをする必要がなく、効率的に、走査電磁石のヒステリシスの影響を受け難い、高精度な照射を行う粒子線治療装置を得られる。

　なお、実施の形態１乃至３で説明した学習アルゴリズムは、ほんの一例であり、最急降下法、遺伝的アルゴリズムなど、他の技術分野で実施されている別のアルゴリズムを応用してもよい。

　また、ランスルーの評価における評価関数は、ここで説明したものに限らず、他の評価関数を用いても構わない。数式（３）や数式（４）のシーケンス番号ｋ毎に重み付けを付けてもよい。また磁場空間内で定義した複数の小領域毎に点数を付ける場合に、小領域毎の点数Ｓｃ_ｉ，ｊに重み付けｗ_ｉ，ｊを行ったものを加えた評価関数でもよい。これらの場合は重要な位置や領域を精度よく最適化を行うことができる。

　この発明に係る粒子線照射装置および粒子線治療装置は、医療用や研究用に用いられる粒子線治療装置に好適に適用できる。

１　　　荷電粒子ビーム　　　　　　　　　１ａ　　入射荷電粒子ビーム
１ｂ　　出射荷電粒子ビーム　　　　　　　３　　　走査電磁石
３ａ　　Ｘ方向走査電磁石　　　　　　　　３ｂ　　Ｙ方向走査電磁石
４　　　走査電源　　　　　　　　　　　　７　　　ビーム位置モニタ
１１　　線量モニタ　　　　　　　　　　　１５　　照射対象
２０　　磁場センサ
２０ａ　Ｘ方向電磁石用磁場センサ
２０ｂ　Ｙ方向電磁石用磁場センサ　　　　２２　　逆写像演算器
３３　　指令評価器　　　　　　　　　　　３４　　指令更新器
３５　　走査電磁石指令値系列生成器
３５ａ　走査電磁石指令値系列生成器
３５ｂ　走査電磁石指令値系列生成器
３７　　走査電磁石指令値学習生成器
３７ａ　走査電磁石指令値学習生成器
３７ｂ　走査電磁石指令値学習生成器
３７ｃ　走査電磁石指令値学習生成器
３７ｄ　走査電磁石指令値学習生成器
３７ｅ　走査電磁石指令値学習生成器　　　３８　　パラメータ更新器
５１　　ビーム発生装置　　　　　　　　　５２　　加速器
５３　　ビーム輸送装置　　　　　　　　　５４　　粒子線照射装置
６０　　ＮＮ（ニューラルネットワーク）　Ｐｓ　　測定位置座標
^ａＰ_ｋ　　照射位置　　　　　　　　　　　　Ｐ_ｋ　　照射位置
Ｄｓ　　測定線量　　　　　　　　　　　　^ａＤ_ｋ　　測定線量
Ｄｉ　　目標線量　　　　　　　　　　　　Ｂｓ　　測定磁場
^ａＢ_ｋ　　測定磁場　　　　　　　　　　　　Ｂ_ｅｓｔ　磁場の推定値
Ｂ_ｋ　　磁場の推定値　　　　　　　　　　　Ｄｓｓ　領域測定線量
Ｄｓｃ　領域線量計算値　　　　　　　　　ｄｅ　　線量誤差
Ｔ　　　点数表　　　　　　　　　　　　　Ｓ_ｉ，ｊ　磁場小領域
Ｉ_ｋ　　指令電流

Claims

　加速器により加速された荷電粒子ビームを走査し、ヒステリシスを有する走査電磁石と、
前記走査電磁石を駆動する励磁電流を出力する走査電源と、
前記走査電源を制御する照射制御装置とを備えた粒子線照射装置であって、
前記照射制御装置は、
前記走査電源に出力された前記励磁電流の指令値による一連の照射動作であるランスルーの結果を評価し、前記評価の結果が所定の条件を満たさない場合に前記励磁電流の指令値を更新して前記ランスルーが実行され、前記評価の結果が所定の条件を満たした前記励磁電流の指令値を前記走査電源に出力する走査電磁石指令値学習生成器を有したことを特徴とする粒子線照射装置。
　前記走査電磁石指令値学習生成器は、
前記励磁電流の指令値による前記ランスルーの結果を評価する指令評価器と、前記指令評価器による評価結果が前記所定の条件を満たさない場合に前記励磁電流の指令値を更新する指令更新器と、前記指令更新器で更新された前記励磁電流の指令値により前記ランスルーが実行され、前記評価の結果が前記所定の条件を満たした前記励磁電流の指令値を前記走査電源に出力する走査電磁石指令値系列生成器を有したことを特徴とする請求項１記載の粒子線照射装置。
　前記走査電磁石指令値学習生成器は、
前記励磁電流の指令値を生成しパラメータを有する数学モデルを備え、
前記励磁電流の指令値によるランスルーの結果を評価する指令評価器と、前記指令評価器による評価結果が前記所定の条件を満たさない場合に前記数学モデルのパラメータを更新するパラメータ更新器と、前記パラメータが更新された前記数学モデルが生成した前記励磁電流の指令値により前記ランスルーが実行され、前記評価の結果が前記所定の条件を満たした前記励磁電流の指令値を前記走査電源に出力する走査電磁石指令値系列生成器を有したことを特徴とする請求項１記載の粒子線照射装置。
　前記荷電粒子ビームは照射対象に向けて照射され、前記荷電粒子ビームの通過位置を検出するビーム位置モニタを備え、
前記指令評価器は、前記通過位置から計算された前記照射対象での照射位置の座標である測定位置座標と、前記照射対象における前記荷電粒子ビームの目標照射位置座標とを比較して評価を行うことを特徴とする請求項２または３に記載の粒子線照射装置。
　前記荷電粒子ビームの線量を測定する線量モニタを備え、
前記指令評価器は、前記線量モニタにより測定された測定線量と前記荷電粒子ビームが照射される照射対象に対する目標線量とを比較して評価を行うことを特徴とする請求項２または３に記載の粒子線照射装置。
　前記励磁電流によって駆動された前記走査電磁石の磁場を測定する磁場センサを備え、前記指令評価器は、前記磁場センサにより測定された測定磁場と、前記荷電粒子ビームが照射される照射対象における前記荷電粒子ビームの目標照射位置座標に対応する前記走査電磁石の磁場の推定値とを比較して評価を行うことを特徴とする請求項２または３に記載の粒子線照射装置。
　前記荷電粒子ビームを照射する照射対象は、前記励磁電流によって駆動される前記走査電磁石が発生するＸ方向及びＹ方向の磁場で定義された磁場小領域に分けられ、
前記指令評価器は、前記磁場小領域毎に前記測定線量と前記目標線量とを比較して評価を行うことを特徴とする請求項５記載の粒子線照射装置。
　前記励磁電流によって駆動された前記走査電磁石の磁場を測定する磁場センサを備え、前記荷電粒子ビームを照射する照射対象は、前記励磁電流によって駆動される前記走査電磁石が発生するＸ方向及びＹ方向の磁場で定義された磁場小領域に分けられ、
前記指令評価器は、前記荷電粒子ビームが前記磁場小領域に滞在する滞在時間に基づいて計算された領域線量計算値と、前記荷電粒子ビームが照射される照射対象に対する目標線量とを前記磁場小領域毎に比較して評価を行うことを特徴とする請求項２または３に記載の粒子線照射装置。
　前記指令評価器は、前記磁場小領域における前記測定線量である領域測定線量と前記目標線量との差である線量誤差を点数化する点数表を有し、前記点数表により評価を行うことを特徴とする請求項７記載の粒子線照射装置。
　前記点数表において、前記測定線量が前記目標線量を超過した場合における点数の変化率の絶対値は、前記測定線量が前記目標線量に満たない場合における点数の変化率の絶対値よりも大きいことを特徴とする請求項９記載の粒子線照射装置。
　前記指令評価器は、前記領域線量計算値と前記目標線量との差である線量誤差を点数化する点数表を有し、前記点数表により評価を行うことを特徴とする請求項８記載の粒子線照射装置。
　前記点数表において、前記領域線量計算値が前記目標線量を超過した場合における点数の変化率の絶対値は、前記領域線量計算値が前記目標線量に満たない場合における点数の変化率の絶対値よりも大きいことを特徴とする請求項１１記載の粒子線照射装置。
　前記ランスルーの評価は、前記荷電粒子ビームが照射される照射対象を含む照射可能範囲にて行うことを特徴とする請求項１乃至１２のいずれか１項に記載の粒子線照射装置。
　前記荷電粒子ビームが照射される照射対象における前記荷電粒子ビームの目標照射位置座標から目標磁場を演算する逆写像演算器を有し、
前記ランスルーの初回において前記走査電源に出力された前記励磁電流の指令値は、前記目標磁場に基づいて生成されたことを特徴とする請求項１または２に記載の粒子線照射装置。
　前記逆写像演算器は、前記目標照射位置座標の多項式からなることを特徴とする請求項１４記載の粒子線照射装置。
　荷電粒子ビームを発生させるビーム発生装置と、前記ビーム発生装置で発生された前記荷電粒子ビームを加速する加速器と、前記加速器により加速された荷電粒子ビームを輸送するビーム輸送装置と、前記ビーム輸送装置で輸送された荷電粒子ビームを走査電磁石で走査して照射対象に照射する粒子線照射装置とを備え、
前記粒子線照射装置は、請求項１乃至１５記載の粒子線照射装置であることを特徴とする粒子線治療装置。