(57)【要約】【課題】粒子線の照射において、外部に放出された放射線の計測画像から的確な線量分布を予測するための線量分布予測システム、線量分布予測方法及び線量分布予測プログラムを提供する。【解決手段】支援サーバ２０は、ＣＴ画像情報記憶部２３、照射条件情報記憶部２４を備える。そして、制御部２１は、粒子線の照射条件、照射対象の元素組成に基づいて、照射対象に対する粒子線の照射時の放射能分布と、外部に放出される放射線の計測画像とをシミュレーションにより算出し、照射条件、放射能分布、計測画像を用いて機械学習を行なうことにより、放射能分布学習モデルを生成する学習処理と、実際の照射条件において、照射対象に対する粒子線の照射時に外部に放出された放射線の実測画像を取得し、照射条件、実測画像を放射能分布学習モデルに適用し、放射能分布を予測し、予測した放射能分布により、線量分布を予測する予測処理とを実行する。【選択図】図１
2019030648.tif 000002
【特許請求の範囲】
【請求項１】
粒子線の照射条件、照射対象の元素組成に関する情報を記憶する記憶部に接続された制御部を備えた線量分布予測システムであって、
前記制御部が、
粒子線の照射条件、照射対象の元素組成に基づいて、前記照射対象に対する粒子線の照射時の放射能分布と、外部に放出される放射線の計測画像とをシミュレーションにより算出し、
前記照射条件、前記放射能分布、前記計測画像を用いて機械学習を行なうことにより、放射能分布学習モデルを生成する学習処理と、
実際の照射条件において、照射対象に対する粒子線の照射時に外部に放出された放射線の実測画像を取得し、
前記照射条件、前記実測画像を前記放射能分布学習モデルに適用し、放射能分布を予測し、
前記予測した放射能分布により、線量分布を予測する予測処理とを実行することを特徴とする線量分布予測システム。
【請求項２】
前記制御部が、
液管が配置された照射対象において、液流がある照射対象及び前記液流がない照射対象に粒子線を照射して取得した放射線の計測画像と、前記液管の配置情報とを用いて、機械学習を行なうことにより、流体拡散学習モデルを生成し、
前記流体拡散学習モデルを用いて、前記放射線の実測画像の高画質化処理を行なうことを特徴とする請求項１に記載の線量分布予測システム。
【請求項３】
前記制御部が、
放射線の計測画像と、前記放射線の計測画像にノイズを付加したノイズ加算計測画像とを用いて、機械学習を行なうことにより、ノイズ学習モデルを生成し、
前記ノイズ学習モデルを用いて、前記放射線の実測画像の高画質化処理を行なうことを特徴とする請求項１又は２に記載の線量分布予測システム。
【請求項４】
前記制御部が、
先行照射において予測した放射能分布に基づいて、後続の計測画像に影響する先行影響画像を生成し、
前記放射線の実測画像の計測点から前記先行影響画像の計測点を差し引くことを特徴とする請求項１〜３のいずれか一項に記載の線量分布予測システム。
【請求項５】
前記制御部が、
放射線の計測画像と、前記放射線の計測画像を圧縮した圧縮計測画像とを用いて、機械学習を行なうことにより、超解像学習モデルを生成し、
前記超解像学習モデルを用いて、前記放射線の実測画像の高画質化処理を行なうことを特徴とする請求項１〜４のいずれか一項に記載の線量分布予測システム。
【請求項６】
前記制御部が、
前記予測した放射能分布をペンシルビームの配置に変換し、前記ペンシルビームによる近似線量分布を予測することを特徴とする請求項１〜５のいずれか一項に記載の線量分布予測システム。
【請求項７】
前記制御部が、
照射条件により生成した本来のペンシルビームに対して、位置をずらしたペンシルビームによりバリエーションを生成し、
前記バリエーションにより、第１の放射能分布を予測し、
前記予測した第１の放射能分布に基づいて計測画像を予測し、
前記予測した計測画像を得られるように、最適化した第２の放射能分布を予測し、
前記照射条件、前記予測した計測画像、前記第１、第２の放射能分布を用いて機械学習を行なうことにより、前記放射能分布学習モデルを生成することを特徴とする請求項１〜６のいずれか一項に記載の線量分布予測システム。
【請求項８】
前記制御部が、
前記予測した放射能分布をペンシルビームの配置に変換し、前記ペンシルビームによる近似線量分布と、照射条件に基づいて生成した３次元線量分布と、照射条件とを用いて、線量分布推定モデルを生成し、
前記線量分布推定モデルを用いて、前記放射能分布から３次元線量分布を予測することを特徴とする請求項１〜７のいずれか一項に記載の線量分布予測システム。
【請求項９】
前記制御部が、
粒子線の照射後の第１時刻の第１実測画像と、前記第１時刻よりも後の第２時刻の第２実測画像とを用いて、画像予測モデルを生成し、
前記画像予測モデルを用いて、粒子線の照射後の第１時刻の実測画像から、第２時刻の第２画像を予測することを特徴とする請求項１〜８のいずれか一項に記載の線量分布予測システム。
【請求項１０】
粒子線の照射条件、照射対象の元素組成に関する情報を記憶する記憶部に接続された制御部を備えた線量分布予測システムを用いて、線量分布予測を行なう方法であって、
前記制御部が、
粒子線の照射条件、照射対象の元素組成に基づいて、前記照射対象に対する粒子線の照射時の放射能分布と、外部に放出される放射線の計測画像とをシミュレーションにより算出し、
前記照射条件、前記放射能分布、前記計測画像を用いて機械学習を行なうことにより、放射能分布学習モデルを生成する学習処理と、
実際の照射条件において、照射対象に対する粒子線の照射時に外部に放出された放射線の実測画像を取得し、
前記照射条件、前記実測画像を前記放射能分布学習モデルに適用し、放射能分布を予測し、
前記予測した放射能分布により、線量分布を予測する予測処理とを実行することを特徴とする線量分布予測方法。
【請求項１１】
粒子線の照射条件、照射対象の元素組成に関する情報を記憶する記憶部に接続された制御部を備えた線量分布予測システムを用いて、線量分布予測を行なうプログラムであって、
前記制御部を、
粒子線の照射条件、照射対象の元素組成に基づいて、前記照射対象に対する粒子線の照射時の放射能分布と、外部に放出される放射線の計測画像とをシミュレーションにより算出し、
前記照射条件、前記放射能分布、前記計測画像を用いて機械学習を行なうことにより、放射能分布学習モデルを生成する学習処理と、
実際の照射条件において、照射対象に対する粒子線の照射時に外部に放出された放射線の実測画像を取得し、
前記照射条件、前記実測画像を前記放射能分布学習モデルに適用し、放射能分布を予測し、
前記予測した放射能分布により、線量分布を予測する予測処理とを実行する手段として機能させることを特徴とする線量分布予測プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、照射した粒子線の線量分布を予測するための線量分布予測システム、線量分布予測方法及び線量分布予測プログラムに関する。
【背景技術】
【０００２】
Ｘ線放射線治療に用いるＸ線のエネルギーは、体の表面では大きく、体内の深いところでは小さくなってしまう。一方、粒子線の一つである陽子線は、入射陽子が体内で停止する寸前の場所で大きなエネルギーを損失し、その場所に「ブラッグピーク」と呼ばれる高線量領域を形成する。このため、正常な領域のダメージを減らし、体内の患部に強い放射線を集中的に照射することができる。このような粒子線を用いた治療においては、患者の体内での線量分布の把握が大切である。
【０００３】
陽子線を用いた治療において、実臨床で利用可能なポジトロン放出核種分布のシミュレーションも検討されている（例えば、特許文献１参照）。この文献に記載された技術では、特定の人体組織の構成元素を含む標準物質に陽子線を照射することにより生成されたポジトロン放出核種のアクティビティ分布の実測値を用いて人体組織の各構成元素のポジトロン放出核種のアクティビティ分布を計算する。そして、このアクティビティ分布の計算値を用い、人体組織の元素組成比に基づいて人体組織に陽子線を照射した場合のポジトロン放出核種のアクティビティ分布を計算する。
【０００４】
図１８を用いて、線量分布を予測する方法を説明する。
治療において、陽子線照射を行なう（ステップＳ００１）。この陽子線照射により、体内に線量分布が生じる。陽子線を物体に照射すると、陽子と物体の原子核との相互作用により、原子核破砕反応が起こり、陽電子放出核（放射化された原子核）が生成される（ステップＳ００２）。この陽電子放出核の分布を３次元アクティビティ分布と呼ぶ。
【０００５】
そして、陽電子放出核は、核種毎の半減期でβ＋崩壊する（ステップＳ００３）。この場合、新たな原子核や陽電子やニュートリノが放出される（ステップＳ００４）。そして、陽電子は電子と対消滅して、その際に電子と陽電子のエネルギーに等価なγ線（放射線）を放出する（ステップＳ００５）。
【０００６】
このγ線を計測することにより（ステップＳ００６）、２次元アクティビティ計測画像を取得する（ステップＳ００７）。そして、この取得した２次元アクティビティ計測画像を用いて、３次元アクティビティ分布、線量分布を予測することができる。
【０００７】
また、第１の情報に基づいて、第２の情報を予測するために、いわゆる「カップリング学習」技術も検討されている（例えば、特許文献２参照）。この文献に記載された技術では、画像処理システムは、学習対象を記録する画像記憶部と、学習結果を記録する学習結果記憶部と、入力部に接続された制御部とを備える。そして、制御部は、学習対象毎に、算出した第１の特徴ベクトル及び第２の特徴ベクトルの主成分係数ベクトルを、主成分係数ベクトルのばらつきを表わす指標で除算した行ベクトルを生成し、主成分分析を行なうことにより、結合主成分基底ベクトル、平均ベクトルを生成し、学習結果記憶部に記録する。そして、処理対象の第１の特徴ベクトルの主成分係数ベクトル、学習結果記憶部に記録された結合主成分基底ベクトルを用いて、第２の特徴ベクトルの主成分係数ベクトルを算出し、処理対象の第２の特徴ベクトルを算出する。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開２０１０−１７２４２７号公報
【特許文献２】特開２０１６−９５６５１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００９】
従来、陽子線照射された体内の３次元線量分布は、２次元のγ線測定結果（２次元アクティビティ計測画像）を用いて確認されていた。このγ線計測においては、照射領域を囲むようにして計測面を並べ、トモグラフィを利用し、体内の３次元アクティビティ分布を求める陽電子断層法（ＰＥＴ：Positron Emission Tomography）を用いる場合がある。しかしながら、計測した２次元アクティビティ計測画像から、的確な線量分布の予測は困難であった。例えば、陽子線照射装置内では、多くの計測面を設置することができないため、全方位に放出されるγ線について限られた計測面から３次元アクティビティ分布を推定する必要がある。
【課題を解決するための手段】
【００１０】
上記課題を解決する線量分布予測システムは、粒子線の照射条件、照射対象の元素組成に関する情報を記憶する記憶部に接続された制御部を備える。そして、制御部が、粒子線の照射条件、照射対象の元素組成に基づいて、前記照射対象に対する粒子線の照射時の放射能分布と、外部に放出される放射線の計測画像とをシミュレーションにより算出し、前記照射条件、前記放射能分布、前記計測画像を用いて機械学習を行なうことにより、放射能分布学習モデルを生成する学習処理と、実際の照射条件において、照射対象に対する粒子線の照射時に外部に放出された放射線の実測画像を取得し、前記照射条件、前記実測画像を前記放射能分布学習モデルに適用し、放射能分布を予測し、前記予測した放射能分布により、線量分布を予測する予測処理とを実行する。
【発明の効果】
【００１１】
本発明によれば、粒子線の照射において、外部に放出された放射線の計測画像から的確な線量分布を予測することができる。
【図面の簡単な説明】
【００１２】
【図１】第１の実施形態の線量分布予測システムの説明図。
【図２】第１の実施形態の処理手順の説明図であって、（ａ）は学習処理、（ｂ）は予測処理の説明図。
【図３】第１の実施形態の処理手順の説明図。
【図４】第１の実施形態の処理手順の説明図。
【図５】第１の実施形態の処理手順の説明図。
【図６】第１の実施形態の処理手順の説明図。
【図７】第１の実施形態の処理手順の説明図。
【図８】第１の実施形態の処理手順の説明図。
【図９】第１の実施形態の処理手順の説明図。
【図１０】第１の実施形態の処理手順の説明図。
【図１１】第１の実施形態の処理手順の説明図。
【図１２】第１の実施形態の処理手順の説明図。
【図１３】第１の実施形態の処理手順の説明図。
【図１４】他の実施形態の処理手順の説明図。
【図１５】第２の実施形態の処理手順の説明図。
【図１６】第２の実施形態の処理手順の説明図。
【図１７】第３の実施形態の処理手順の説明図。
【図１８】従来技術における陽子線照射時のγ線計測の説明図。
【発明を実施するための形態】
【００１３】
（第１の実施形態）
以下、図１〜図１３に従って、線量分布予測システム、線量分布予測方法及び線量分布予測プログラムを具体化した一実施形態を説明する。本実施形態では、粒子線としての陽子を患者の患部に照射した場合に計測したアクティビティ計測画像から線量分布を算出する場合を説明する。
【００１４】
アクティビティ計測画像は可能な限り高画質であることが望ましい。しかし、アクティビティ計測画像の画質には、以下の影響があると想定される。
（ａ）血流による陽電子放出核の拡散：陽電子放出核が、特有の半減期でβ＋崩壊するまでに、血流により拡散する場合がある。
【００１５】
（ｂ）ノイズ：原子核の崩壊はポアソン分布に従うため、計測には熱雑音とポアソンノイズが加算される。また、γ線は人体内における散乱により画像が不明瞭になる。
（ｃ）低い解像度：計測面の解像度が低いため、情報量が少ない。
（ｄ）先行照射の影響：複数の角度から順番に照射した場合、先行照射の影響が後続の計測に影響を与える場合がある。
そこで、本実施形態では、これらの要因を機械学習により解決する。
【００１６】
図１に示すように、本実施形態における線量分布予測方法においては、ネットワークを介して接続されたユーザ端末１０、画像取得部１５及び支援サーバ２０を用いる。
【００１７】
ユーザ端末１０は、ユーザが利用するコンピュータ端末である。このユーザ端末１０は、入力された情報を取得する入力部（キーボードやポインティングデバイス等）や、各種情報の出力を行なう出力部（ディスプレイ等）を備える。
画像取得部１５は、患者に対する陽子線照射において、γ線を２次元で計測したアクティビティ計測画像を取得する。
【００１８】
支援サーバ２０は、線量分布予測処理を行なうコンピュータシステム（線量分布予測システム）である。この支援サーバ２０は、制御部２１、学習結果記憶部２２、ＣＴ画像情報記憶部２３、照射条件情報記憶部２４を備えている。
【００１９】
制御部２１は、制御手段（ＣＰＵ、ＲＡＭ、ＲＯＭ等）を備え、後述する処理（シミュレーション段階、学習段階、予測段階、血管情報取得段階、残存影響処理段階等の各処理）を行なう。そのための線量分布予測プログラムを実行することにより、制御部２１は、シミュレータ２１０、学習部２１１、予測部２１２、血管情報取得部２１３、残存影響処理部２１４として機能する。
【００２０】
シミュレータ２１０は、機械学習に用いる学習データを生成する。
学習部２１１は、学習データを用いた機械学習により、予測を行なうための各種モデル（学習結果）を生成する。
【００２１】
予測部２１２は、学習部２１１によって生成された各種モデル（学習結果）を用いて、各種予測を行なう。ここでは、本来の高画質なアクティビティ計測画像、３次元アクティビティ分布、線量分布の予測を行なう。
【００２２】
本実施形態では、機械学習の手法として、特許文献２に開示された「カップリング学習」を用いる。
ここで、図２を用いて、この「カップリング学習」について説明する。
図２（ａ）に示すように、学習部２１１は、説明変数ベクトル５０１（予測時に入力する情報）と目的変数ベクトル５０２（予測する情報）を結合した行ベクトルを作成し、データ数分、縦に並べて行列５０３を作成する。そして、この行列を主成分分析することで、説明変数と目的変数を同時に表現する基底５０４を求める。
次に、学習部２１１は、前半部分（説明変数部分）について直交化して基底Ｄi,jとする。後半部分（目的変数部分）についても、直交化時と同じ線形結合係数で線形変換して基底Ｅi,kとする。前半部分を直交化しておくことで、予測時に線形結合係数を内積で求めることができる。
【００２３】
図２（ｂ）に示すように、予測部２１２は、入力された説明変数ベクトルＸjと、直交基底i,jとの内積を求めて線形結合係数ｂiを求める。そして、予測部２１２は、この線形結合係数ｂiと基底Ｅi,kの線形結合により、目的変数ベクトルＹkを再構成する。
【００２４】
図１に示す血管情報取得部２１３は、患者の血管（照射対象内の液管）の配置情報を取得する。
残存影響処理部２１４は、先行照射の影響を排除する処理を実行する。
【００２５】
学習結果記憶部２２には、カップリング学習によって生成された学習結果（主成分基底ベクトルと平均ベクトル）が記録される。この学習結果は、機械学習を行なった場合に記録される。本実施形態では、計測高速化モデル、血流拡散モデル、血管配置モデル、ノイズ加算モデル、超解像モデル、アクティビティ分布推定モデル(放射能分布学習モデル)、線量分布推定モデルが記録される。
【００２６】
計測高速化モデルは、第１時刻におけるアクティビティ計測画像に基づいて、第２時刻におけるアクティビティ計測画像を予測するためのモデルである。なお、以下では、第２時刻は第１時刻より遅い時刻という意味で用いている。
血流拡散モデルは、実測したアクティビティ計測画像において、陽電子放出核の血流拡散の影響を排除したアクティビティ計測画像を予測するための流体拡散学習モデルである。
【００２７】
血管配置モデルは、実測した３次元ＣＴ画像から３次元血管造影画像を予測するための学習モデルである。
ノイズ加算モデルは、ノイズが重畳した実測したアクティビティ計測画像から、ノイズを除去したアクティビティ計測画像を予測するためのノイズ学習モデルである。
【００２８】
超解像モデルは、解像度が低い実測したアクティビティ計測画像から、高解像度のアクティビティ計測画像を予測するための超解像学習モデルである。
アクティビティ分布推定モデルは、治療計画情報（ＣＴ画像、照射条件）、アクティビティ計測画像、３次元アクティビティ分布（最適化）から３次元アクティビティ分布（正解）を推定するための放射能分布学習モデルである。
【００２９】
線量分布推定モデルは、治療計画情報（ＣＴ画像、照射条件）、３次元線量分布（ペンシルビーム近似）から３次元線量分布（正解）を推定するための学習モデルである。
【００３０】
ＣＴ画像情報記憶部２３には、患者の治療計画についてのＣＴ画像を管理するためのＣＴ画像管理レコード２３０が記録される。このＣＴ画像管理レコード２３０は、各患者の治療前に３次元コンピュータ断層（ＣＴ）検査を行なった場合に記録される。ＣＴ画像管理レコード２３０には、患者コードに関連付けてＣＴ画像が記録される。
【００３１】
患者コードデータ領域には、各患者を特定するための識別子に関するデータが記録される。
ＣＴ画像データ領域には、この患者の３次元コンピュータ断層画像（３次元ＣＴ画像）データが記録される。この３次元ＣＴ画像により、組織組成を特定したり、血管の配置を予測したりすることができる。
【００３２】
照射条件情報記憶部２４には、患者の治療計画を管理するための照射条件管理レコード２４０が記録される。この照射条件管理レコード２４０は、各患者の治療時の陽子線の照射条件を決定した場合に記録される。照射条件管理レコード２４０には、患者コード、治療予定日に関連付けて照射条件が記録される。
【００３３】
患者コードデータ領域には、各患者を特定するための識別子に関するデータが記録される。
治療予定日データ領域には、この患者に対して、治療計画における陽子線照射による治療の予定日（年月日）に関するデータが記録される。
【００３４】
照射条件データ領域には、この患者に対して、治療予定日に照射する陽子線の照射条件に関するデータが記録される。照射条件としては、陽子線の照射位置、照射方向、照射エネルギー、ビーム照射法等がある。ここで、ビーム照射法には、「拡大ビーム照射法」と「スキャニング照射法」の２手法がある。
【００３５】
（ペンシルビーム近似処理）
次に、図３を用いて、ペンシルビーム近似処理を説明する。陽子線照射のシミュレーションにおいては、同一照射位置・同一照射方向・同一エネルギーの複数の陽子線を束ねて計算するペンシルビーム近似を用いる。
【００３６】
まず、支援サーバ２０の制御部２１は、ペンシルビーム生成処理を実行する（ステップＳ０１１）。具体的には、シミュレータ２１０は、ユーザ端末１０に照射条件入力画面を出力する。この照射条件入力画面においては、患者コード、治療予定日を入力する。この場合、シミュレータ２１０は、照射条件情報記憶部２４から、患者コード、治療予定日に関連付けられた照射条件（陽子線の位置、方向、エネルギー、ビーム照射法等）を取得する。そして、照射条件に応じたペンシルビームを生成する。なお、陽子線と人体組織の核反応は、元素毎に異なるため、ペンシルビームも元素（酸素、水素、炭素、カルシウム等）毎に用意する。
【００３７】
次に、支援サーバ２０の制御部２１は、組成密度分布の取得処理を実行する（ステップＳ０１２）。具体的には、シミュレータ２１０は、ＣＴ画像情報記憶部２３から、患者コードに関連付けられたＣＴ画像を取得する。そして、シミュレータ２１０は、ＣＴ画像から各元素の３次元組成密度分布に変換する。
【００３８】
そして、支援サーバ２０の制御部２１は、元素に応じたペンシルビーム毎に以下の処理を繰り返す。
ここでは、支援サーバ２０の制御部２１は、線量分布の算出処理を実行する（ステップＳ０１３）。具体的には、シミュレータ２１０は、処理対象のペンシルビームによる線量分布を算出し、メモリに記憶する。
【００３９】
次に、支援サーバ２０の制御部２１は、３次元アクティビティ分布の算出処理を実行する（ステップＳ０１４）。具体的には、シミュレータ２１０は、処理対象のペンシルビームによる３次元アクティビティ分布を算出し、メモリに記憶する。
すべてのペンシルビームについて処理を終了するまで繰り返す。
【００４０】
次に、支援サーバ２０の制御部２１は、すべての元素毎、ペンシルビーム毎の算出結果の合成処理を実行する（ステップＳ０１５）。具体的には、シミュレータ２１０は、メモリに記憶したすべての線量分布、３次元アクティビティ分布を、それぞれ積算する。
【００４１】
（高速化処理）
図４を用いて、高速化処理を説明する。アクティビティ計測では、陽電子放出核が放出した陽電子が電子と対消滅した際に生じるγ線を計測する。各陽電子放出核は、半減期に従ってβ＋崩壊するため、通常、その計測には２分程度の時間がかかる。この高速化処理では、計測時間を短縮した初期段階の画像（例えば、計測時間を１分にした画像）に基づいて、将来の画像を機械学習により予測する。
【００４２】
支援サーバ２０の制御部２１は、機械学習処理を実行する（ステップＳ０２１）。具体的には、制御部２１の学習部２１１は、時系列に計測されたアクティビティ計測画像の実測データから、学習用のアクティビティ計測画像（第１時刻）、アクティビティ計測画像（第２時刻）を取得する。このアクティビティ計測画像（第１時刻）は、実測データにおいて、観測開始から６０秒以内の計測のみを集計して作成する。また、アクティビティ計測画像（第２時刻）は、実測データにおいて、観測開始から１２０秒以内の計測のみを集計して作成する。
【００４３】
そして、学習部２１１は、教師データとして、アクティビティ計測画像（第１時刻）から生成した説明変数ベクトル、アクティビティ計測画像（第２時刻）から生成した目的変数ベクトルを用いて、機械学習を行なう。そして、学習部２１１は、算出した計測高速化モデルを生成し、学習結果記憶部２２に記録する。
【００４４】
アクティビティ計測時には、支援サーバ２０の制御部２１は、予測処理を実行する（ステップＳ０２２）。具体的には、制御部２１の予測部２１２は、計測したアクティビティ計測画像（第１時刻）を取得し、このアクティビティ計測画像（第１時刻）を説明変数ベクトルとして、学習結果記憶部２２に記録された計測高速化モデルに適用する。そして、予測部２１２は、目的変数ベクトルとして、アクティビティ計測画像（第２時刻）を算出する。
【００４５】
（血管配置の学習処理）
次に、図５を用いて、血管配置の学習処理を説明する。血流拡散の影響を予測するためには、患者における血管の配置情報が必要である。本実施形態では、機械学習を利用してＣＴ画像から血管配置の推測を行なう。
【００４６】
ここで、支援サーバ２０の制御部２１は、機械学習処理を実行する（ステップＳ０３１）。具体的には、制御部２１の学習部２１１は、複数の患者の血管情報（血管造影画像）と、３次元ＣＴ画像とを教師データとして取得する。そして、学習部２１１は、３次元ＣＴ画像から生成した説明変数ベクトル、血管造影画像から生成した目的変数ベクトルを用いて、機械学習を行なう。そして、学習部２１１は、算出した血管配置モデルを学習結果記憶部２２に記録する。
【００４７】
（血流拡散対応処理）
次に、図６、図７を用いて、血流拡散対応処理を説明する。ここでは、血流による陽電子放出核の拡散によるアクティビティ分布の復元を行なう。血管の配置情報を利用し、事前に血流によるアクティビティの拡散についてシミュレーションを行ない、拡散前後のアクティビティ計測画像を仮想的に作成する。作成した拡散前後のアクティビティ計測画像を教師データとして機械学習を用いた学習処理を実行し、拡散後のアクティビティ計測画像から、拡散前のアクティビティ計測画像を推定する予測処理を実行する。
【００４８】
（血流拡散対応の学習処理）
図６を用いて、血流拡散対応の学習処理を説明する。
まず、支援サーバ２０の制御部２１は、３次元拡散シミュレーション処理を実行する（ステップＳ０４１）。具体的には、制御部２１のシミュレータ２１０は、上述したペンシルビーム近似処理により、各照射条件に応じて、３次元アクティビティ分布を算出する。次に、この患者の血管情報（血管造影画像）に基づいて、血流拡散をシミュレーションする。そして、シミュレータ２１０は、血流拡散前、第１時刻（例えば、３０秒）後、第２時刻（例えば、４０秒）後の各アクティビティ計測画像を生成する。
【００４９】
次に、支援サーバ２０の制御部２１は、機械学習処理を実行する（ステップＳ０４２）。具体的には、制御部２１の学習部２１１は、第１時刻及び第２時刻のアクティビティ計測画像から生成した説明変数ベクトル、血流拡散前のアクティビティ計測画像から生成した目的変数ベクトルを用いて、機械学習を行なう。そして、学習部２１１は、算出した血流拡散モデルを学習結果記憶部２２に記録する。
【００５０】
（血流拡散対応の予測処理）
次に、図７を用いて、血流拡散対応の予測処理を説明する。ここでは、血流による拡散の復元を行なう。２つ以上の時間差情報を利用して、拡散前の画像予測を行なう。
【００５１】
まず、支援サーバ２０の制御部２１は、血管情報の取得処理を実行する（ステップＳ０５１）。具体的には、予め患者の血管情報（血管の配置情報）が準備されている場合には、制御部２１の血管情報取得部２１３は、ユーザ端末１０から血管情報を取得する。一方、ユーザ端末１０から患者の血管情報を取得できない場合には、血管情報取得部２１３は、ＣＴ画像情報記憶部２３から患者のＣＴ画像を取得し、上述した血管配置モデルにより、血管情報を取得する。
【００５２】
次に、支援サーバ２０の制御部２１は、予測処理を実行する（ステップＳ０５２）。具体的には、制御部２１の予測部２１２は、画像取得部１５から、計測した第１時刻（例えば、３０秒）のアクティビティ計測画像を取得する。そして、予測部２１２は、前述した高速化処理により、第１時刻（例えば、３０秒）のアクティビティ計測画像から第２時刻（例えば、４０秒）のアクティビティ計測画像を予測する。次に、予測部２１２は、第１時刻、第２時刻のアクティビティ計測画像、血管情報を説明変数ベクトルとして、学習結果記憶部２２の血流拡散モデルに適用する。そして、予測部２１２は、目的変数ベクトルとして、血流拡散前のアクティビティ計測画像を算出する。
【００５３】
（ノイズ対応処理）
次に、図８を用いて、ノイズ対応処理を説明する。ここでは、アクティビティ計測画像をシミュレーションにより作成し、人工的にノイズを加算し、ノイズ加算前後のアクティビティ計測画像を教師データとして機械学習を行なう。そして、ノイズが加算されたアクティビティ計測画像から、ノイズを除去したアクティビティ計測画像を推定する。
【００５４】
まず、支援サーバ２０の制御部２１は、ノイズシミュレーション処理を実行する（ステップＳ０６１）。具体的には、制御部２１のシミュレータ２１０は、シミュレーションによりアクティビティ計測画像を作成し、熱雑音やポアソンノイズ、散乱等、ノイズの統計的・物理的な性質に基づいて、アクティビティ計測画像にノイズを与える。これにより、アクティビティ計測画像（ノイズ加算前）及びアクティビティ計測画像（ノイズ加算後）を取得する。
【００５５】
次に、支援サーバ２０の制御部２１は、機械学習処理を実行する（ステップＳ０６２）。具体的には、制御部２１の学習部２１１は、説明変数ベクトルとしてアクティビティ計測画像（ノイズ加算後）、目的変数ベクトルとしてアクティビティ計測画像（ノイズ加算前）を用いて、機械学習を行なう。そして、学習部２１１は、算出したノイズ加算モデルを学習結果記憶部２２に記録する。
【００５６】
アクティビティ計測時には、支援サーバ２０の制御部２１は、予測処理を実行する（ステップＳ０６３）。具体的には、予測部２１２は、画像取得部１５から、実測したアクティビティ計測画像を取得する。このアクティビティ計測画像には、ノイズが含まれている。次に、予測部２１２は、アクティビティ計測画像（ノイズを含む）を説明変数ベクトルとして、学習結果記憶部２２のノイズ加算モデルに適用する。そして、予測部２１２は、目的変数ベクトルとして、アクティビティ計測画像（ノイズを含まない）を算出する。
【００５７】
（超解像処理）
次に、図９を用いて、超解像処理を説明する。ここでは、アクティビティ計測画像をシミュレーションにより作成し、縮小することにより、縮小前後のアクティビティ計測画像を教師データとして機械学習を行なう。そして、解像度の低いアクティビティ計測画像から、高精細なアクティビティ計測画像を推定する。
【００５８】
まず、支援サーバ２０の制御部２１は、画像の縮小処理を実行する（ステップＳ０７１）。具体的には、制御部２１のシミュレータ２１０は、シミュレーションによりアクティビティ計測画像を作成し、画像を縮小する。これにより、アクティビティ計測画像（縮小前）及びアクティビティ計測画像（縮小後）を取得する。
【００５９】
次に、支援サーバ２０の制御部２１は、機械学習処理を実行する（ステップＳ０７２）。具体的には、制御部２１の学習部２１１は、説明変数ベクトルとしてアクティビティ計測画像（縮小後）、目的変数ベクトルとしてアクティビティ計測画像（縮小前）を用いて、機械学習を行なう。そして、学習部２１１は、算出した超解像モデルを学習結果記憶部２２に記録する。
【００６０】
アクティビティ計測時には、支援サーバ２０の制御部２１は、予測処理を実行する（ステップＳ０７３）。具体的には、制御部２１の予測部２１２は、処理対象のアクティビティ計測画像を取得する。このアクティビティ計測画像は、解像度が低い画像になっている。次に、予測部２１２は、このアクティビティ計測画像（低解像度）を説明変数ベクトルとして、学習結果記憶部２２の超解像モデルに適用する。そして、予測部２１２は、目的変数ベクトルとして、アクティビティ計測画像（高解像度）を算出する。
【００６１】
（先行照射の影響排除処理）
図１０を用いて、先行照射の影響排除処理を説明する。陽子線を複数の角度から照射する場合、同時に複数の角度から照射はできないため、角度を変えて照射を繰り返す。この場合、アクティビティ計測画像には、先行照射における前の角度でのアクティビティ（影響）が残存する場合がある。このような場合には、計測画像から、前の角度での影響を差し引く必要がある。
【００６２】
まず、支援サーバ２０の制御部２１は、先行照射の３次元アクティビティ分布の推定処理を実行する（ステップＳ０８１）。具体的には、制御部２１の残存影響処理部２１４は、後述する３次元アクティビティ分布の推定処理により、先行照射の３次元アクティビティ分布を算出する。
【００６３】
次に、支援サーバ２０の制御部２１は、先行照射によるアクティビティ計測画像の予測処理を実行する（ステップＳ０８２）。具体的には、制御部２１のシミュレータ２１０は、算出した３次元アクティビティ分布に基づいて、先行影響画像であるアクティビティ計測画像（先行照射）を算出する。この場合、３次元アクティビティ分布における血流拡散の影響の追加、アクティビティ計測画像におけるノイズ加算、縮小化を行なう。
【００６４】
次に、支援サーバ２０の制御部２１は、差分算出処理を実行する（ステップＳ０８３）。具体的には、制御部２１のシミュレータ２１０は、処理対象のアクティビティ計測画像の計測点から、アクティビティ計測画像（先行照射）の計測点を差し引くことにより、先行照射の影響を排除したアクティビティ計測画像を算出する。
【００６５】
（学習データ生成処理）
次に、図１１を用いて、学習データ生成処理を説明する。ここでは、高画質化したアクティビティ計測画像と、治療計画情報（ＣＴ画像、照射条件）から、３次元アクティビティ分布を推定する。治療計画時にシミュレーションした３次元アクティビティ分布と、実際の３次元アクティビティ分布は、以下の変動要因で異なる場合がある。
【００６６】
（ａ）陽子線照射条件の変動（照射装置の位置ずれ等）
（ｂ）患者の位置決め条件のずれ変動
（ｃ）患者の体型変動
（ｄ）患者の体内の対象臓器の状態変化（患部の形状変化等）
上記原因で生じた３次元アクティビティ分布のずれは、アクティビティ計測画像に反映される。このため、アクティビティ計測画像から最適化計算により、３次元アクティビティ分布を推定する。そして、推定した３次元アクティビティ分布は、機械学習による本来の３次元アクティビティ分布（正解）を求める際の入力値（初期値）として利用する。
【００６７】
このため、治療計画のビーム配置に対して、上述した４種類の変動要因をランダムに一部のビームに与えることにより、進行方向前後に位置オフセット（位置ずれ）を与えた３次元アクティビティ分布（正解）と、アクティビティ計測画像のセットを大量にシミュレーションにより作成する。
【００６８】
ここでは、ペンシルビームの変動バリエーション毎に以下の処理を繰り返す。ここで、変動バリエーションとして、上述した４種類の変動要因（ａ）〜（ｄ）を想定する。
まず、支援サーバ２０の制御部２１は、３次元線量分布シミュレーション処理を実行する（ステップＳ０９１）。具体的には、シミュレータ２１０は、治療計画情報（ＣＴ画像、照射条件）において、変動バリエーションに基づいて、上述したペンシルビーム近似を行なう。そして、シミュレータ２１０は、変動バリエーションに応じた３次元線量分布（正解）を算出する。
【００６９】
次に、支援サーバ２０の制御部２１は、３次元アクティビティ分布シミュレーション処理を実行する（ステップＳ０９２）。具体的には、シミュレータ２１０は、治療計画情報（ＣＴ画像、照射条件）において、変動バリエーションに基づいて、上述したペンシルビーム近似を行なう。そして、変動バリエーションに応じた３次元アクティビティ分布（正解）を算出する。
【００７０】
次に、支援サーバ２０の制御部２１は、アクティビティ計測画像の作成処理を実行する（ステップＳ０９３）。具体的には、シミュレータ２１０は、算出した３次元アクティビティ分布を２次元計測面に投影することにより、アクティビティ計測画像を作成する。
【００７１】
次に、支援サーバ２０の制御部２１は、最適化処理を実行する（ステップＳ０９４）。具体的には、シミュレータ２１０は、ペンシルビーム配置を、アクティビティ計測画像との差異が最も少なくなるように再配置を行なう。ここでは、ビーム進行方向をＺ軸とすると、各ペンシルビームは進行方向（Ｚ軸）には独立に移動させることができるが、ＸＹ方向の移動は全体の平行移動のみを許容する。すなわち、各ペンシルビームの照射開始位置を初期値とし、最急降下法により、組成毎の各ペンシルビームの照射開始位置を、進行方向に前後させて、アクティビティ計測画像との差異が最も少ないビーム位置を求める。なお、ビーム照射開始位置を前後させる際には、改めて、その照射開始位置におけるペンシルビームの照射シミュレーションを行ない、３次元アクティビティ分布を求める。
以上の処理を、予め準備したペンシルビームの変動バリエーションについて繰り返す。
【００７２】
（３次元アクティビティ分布の推定処理）
次に、図１２を用いて、３次元アクティビティ分布の推定処理を説明する。ここでは、前述した学習データ生成処理において作成したデータを教師データとして用いる。第１の放射能分布としてのアクティビティ計測画像、第２の放射能分布としての３次元アクティビティ分布（最適化）、治療計画情報（照射条件、ＣＴ画像）に基づいて、機械学習を行なうことにより、高精度な３次元アクティビティ分布を求める。ここでは、治療計画情報とシミュレーションとにより作成したアクティビティ計測画像から、３次元アクティビティ分布（最適化）を作成する。そして、アクティビティ計測画像、３次元アクティビティ分布（最適化）、治療計画情報から、３次元アクティビティ分布（正解）を推定する予測モデルを機械学習により求める。
【００７３】
（３次元アクティビティ分布の学習処理）
まず、支援サーバ２０の制御部２１は、機械学習処理を実行する（ステップＳ１０１）。具体的には、制御部２１の学習部２１１は、治療計画情報として、ＣＴ画像情報記憶部２３から患者のＣＴ画像、照射条件情報記憶部２４から照射条件を取得する。次に、学習部２１１は、学習データ生成処理において作成した３次元アクティビティ分布（最適化）、アクティビティ計測画像、３次元アクティビティ分布（正解）を取得する。そして、学習部２１１は、治療計画情報（ＣＴ画像、照射条件）、アクティビティ計測画像、３次元アクティビティ分布（最適化）から生成した説明変数ベクトル、３次元アクティビティ分布（正解）から生成した目的変数ベクトルを用いて、機械学習を行なう。そして、学習部２１１は、算出したアクティビティ分布推定モデルを学習結果記憶部２２に記録する。
【００７４】
（３次元アクティビティ分布の予測処理）
次に、支援サーバ２０の制御部２１は、予測処理を実行する（ステップＳ１０２）。具体的には、制御部２１の予測部２１２は、治療計画情報として、ＣＴ画像情報記憶部２３から患者のＣＴ画像、照射条件情報記憶部２４から照射条件を取得する。次に、予測部２１２は、照射条件に応じたペンシルビーム近似に対して、アクティビティ計測画像（実測、高画質化後）を用いて、ステップＳ０９４と同様に最適化した３次元アクティビティ分布（最適化）を算出する。そして、予測部２１２は、アクティビティ計測画像（実測、高画質化後）、３次元アクティビティ分布（最適化）、治療計画情報（照射条件、ＣＴ画像）を説明変数ベクトルとして、学習結果記憶部２２に記録されたアクティビティ分布推定モデルに適用する。そして、予測部２１２は、目的変数ベクトルとして、３次元アクティビティ分布（予測結果）を算出する。
【００７５】
（３次元線量分布の推定処理）
次に、図１３を用いて、３次元線量分布の推定処理を説明する。ここでも、前述した学習データ生成処理において作成したデータを教師データとして用いる。ここでは、３次元アクティビティ分布から３次元線量分布を求める。３次元アクティビティ分布から３次元線量分布に直接変換することはできないため、ペンシルビーム近似を通して変換を行なう。
【００７６】
（３次元線量分布の学習処理）
まず、支援サーバ２０の制御部２１は、３次元アクティビティ分布の取得処理を実行する（ステップＳ１１１）。具体的には、制御部２１のシミュレータ２１０は、３次元アクティビティ分布の推定処理で算出した３次元アクティビティ分布を取得する。
【００７７】
次に、支援サーバ２０の制御部２１は、ペンシルビームへの変換処理を実行する（ステップＳ１１２）。具体的には、制御部２１のシミュレータ２１０は、取得した３次元アクティビティ分布を、各ペンシルビームの配置情報（照射開始位置）に変換する。ここでは、アクティビティ計測画像との差異が最も少なくなるように再配置した各ペンシルビームの配置情報（照射開始位置）を初期値とし、最急降下法により、組成毎に各ペンシルビームの照射開始位置を、進行方向に前後させて、算出した３次元アクティビティ分布との差異が最も少ないビーム位置を求める。なお、ビーム照射開始位置を前後させる際には、改めて、その照射開始位置におけるペンシルビームの照射シミュレーションを行ない、３次元アクティビティ分布を求める。
【００７８】
次に、支援サーバ２０の制御部２１は、ペンシルビーム近似の３次元線量分布の算出処理を実行する（ステップＳ１１３）。具体的には、制御部２１のシミュレータ２１０は、求めたビーム照射開始位置におけるペンシルビーム近似の照射シミュレーションを行ない、３次元線量分布を求める。
【００７９】
次に、支援サーバ２０の制御部２１は、機械学習処理を実行する（ステップＳ１１４）。具体的には、制御部２１の学習部２１１は、治療計画情報として、ＣＴ画像情報記憶部２３から患者のＣＴ画像、照射条件情報記憶部２４から照射条件を取得する。そして、学習部２１１は、ペンシルビーム近似の３次元線量分布から生成した説明変数ベクトル、学習データ生成処理により生成した３次元線量分布（正解）から生成した目的変数ベクトルを用いて、機械学習を行なう。そして、学習部２１１は、算出した線量分布推定モデルを学習結果記憶部２２に記録する。
【００８０】
（３次元線量分布の予測処理）
まず、支援サーバ２０の制御部２１は、３次元アクティビティ分布の取得処理を実行する（ステップＳ１２１）。具体的には、制御部２１のシミュレータ２１０は、画像取得部１５から取得したアクティビティ計測画像に基づいて、３次元アクティビティ分布の推定処理において算出した３次元アクティビティ分布（予測結果）を取得する。
【００８１】
次に、支援サーバ２０の制御部２１は、ステップＳ１１２と同様に、ペンシルビームへの変換処理を実行する（ステップＳ１２２）。
次に、支援サーバ２０の制御部２１は、ステップＳ１１３と同様に、ペンシルビーム近似の３次元線量分布の算出処理を実行する（ステップＳ１２３）。
【００８２】
次に、支援サーバ２０の制御部２１は、予測処理を実行する（ステップＳ１２４）。具体的には、制御部２１の予測部２１２は、ペンシルビーム近似の３次元線量分布、治療計画情報を説明変数ベクトルとして、学習結果記憶部２２に記録された線量分布推定モデルに適用する。そして、予測部２１２は、目的変数ベクトルとして、３次元線量分布（予測結果）を算出する。
【００８３】
以上、本実施形態によれば、以下に示す効果を得ることができる。
（１−１）本実施形態では、支援サーバ２０の制御部２１は、ペンシルビーム近似処理を実行する。これにより、計算量を削減し、計算の高速化を図ることができる。また、人体内の各元素の密度（組成密度分布）は一定ではないが、元素毎のペンシルビームシミュレーションを行ない、その結果を合計して、３次元線量分布や３次元アクティビティ分布を求めることができる。
【００８４】
（１−２）本実施形態では、支援サーバ２０の制御部２１は、高速化処理を実行する。これにより、短時間の計測により、将来のアクティビティ計測画像を予測し、計測の高速化を図ることができる。
【００８５】
（１−３）本実施形態では、支援サーバ２０の制御部２１は、血流対応処理を実行する。人体に陽子線を照射した場合、生成された陽電子放出核が血流によって、時間とともに拡散する可能性があるが、拡散前の線量分布を予測することができる。
【００８６】
（１−４）本実施形態では、支援サーバ２０の制御部２１は、ノイズ対応処理を実行する。これにより、γ線の計測時に重畳したノイズを考慮して、高画質なアクティビティ計測画像を予測することができる。
【００８７】
（１−５）本実施形態では、支援サーバ２０の制御部２１は、超解像処理を実行する。これにより、低解像度のアクティビティ計測画像に対して、高解像度のアクティビティ計測画像を予測することができる。
【００８８】
（１−６）本実施形態では、支援サーバ２０の制御部２１は、先行照射の影響排除処理を実行する。治療等において、連続して陽電子照射を行なう場合、先行照射が後続の計測に影響を与える可能性がある。このような場合にも、先行照射の影響を排除して、今回の陽電子照射における線量分布を的確に予測することができる。
【００８９】
（１−７）本実施形態では、支援サーバ２０の制御部２１は、学習データ生成処理を実行する。状況によって、照射するペンシルビームの配置がずれる可能性がある。このような状況においても、アクティビティ計測画像に応じたペンシルビームの配置を考慮し、３次元アクティビティ分布や線量分布を予測するための教師データを生成することができる。すなわち、全体の位置ずれ等のバリエーションは考慮する必要がなくなり、個別のペンシルビームの位置ずれのバリエーションに限定することが可能となるため、学習時間の削減と高精度化が期待できる。
【００９０】
（１−８）本実施形態では、支援サーバ２０の制御部２１は、３次元アクティビティ分布の推定処理を実行する。これにより、治療計画情報に基づいて、３次元アクティビティ分布を予測することができる。
（１−９）本実施形態では、支援サーバ２０の制御部２１は、３次元線量分布の推定処理を実行する。これにより、治療計画情報に基づいて、３次元線量分布を予測することができる。
【００９１】
また、上記実施形態は、以下の態様に変更してもよい。
・上記実施形態では、高速化〜３次元線量分布において機械学習を行ない、この学習結果に基づいて予測を行なう。すべてにおいて機械学習を用いる必要はない。
・上記実施形態では、粒子線として陽子線を用いる。ここで、粒子線は陽子線に限定されるものではなく、例えば、炭素線等を用いることも可能である。
・上記実施形態では、血管配置に基づく血流拡散について拡散前を予測する。拡散の要因は血流に限定されるものではなく、液体等の流体の所在情報（液管情報）を用いて、この流体によるアクティビティ分布の拡散に対応するようにしてもよい。
【００９２】
・上記実施形態では、血管配置学習処理において、機械学習を利用してＣＴ画像から推測を行なう。血管配置の予測方法は、機械学習に限定されるものではない。アトラス（解剖学）情報から推測するようにしてもよい。
【００９３】
・上記実施形態では、支援サーバ２０の制御部２１は、３次元アクティビティ分布の推定処理を実行する。ここでは、３次元アクティビティ分布を予測するが、予測対象は３次元アクティビティ分布そのものに限定されるものではない。例えば、推定した分布との差分「３次元アクティビティ分布（正解）−３次元アクティビティ分布（最適化）」を予測するようにしてもよい。
【００９４】
・上記実施形態では、支援サーバ２０の制御部２１は、３次元線量分布の推定処理を実行する。ここでは、３次元線量分布を予測するが、予測対象は３次元線量分布そのものに限定されるものではない。例えば、推定した分布との差分「３次元線量分布（正解）−３次元線量分布（ペンシルビーム近似）」を予測するようにしてもよい。
【００９５】
・上記実施形態では、ノイズ対応処理、超解像処理、先行照射の影響排除処理において、画質劣化要因について個別に高画質化するための機械学習を個別に行ない、それぞれの推定モデルを生成する。これらをまとめて、機械学習するようにしてもよい。ここでは、統合処理により機械学習を行なう。
【００９６】
図１４を用いて、統合処理を説明する。
まず、支援サーバ２０の制御部２１は、３次元アクティビティ分布シミュレーション処理を実行する（ステップＳ１３１）。具体的には、制御部２１のシミュレータ２１０は、血管情報（血管造影画像）を用いて、ペンシルビーム近似シミュレーションにより、アクティビティ計測画像（血流拡散前）、アクティビティ計測画像（血流拡散後）を生成する。
【００９７】
次に、支援サーバ２０の制御部２１は、低画質化シミュレーション処理を実行する（ステップＳ１３２）。具体的には、制御部２１のシミュレータ２１０は、ノイズシミュレーション、画像の縮小、先行照射の影響を加算して、ノイズ・縮小・先行照射の影響を受けた低画質化アクティビティ計測画像を生成する。
【００９８】
次に、支援サーバ２０の制御部２１は、機械学習処理を実行する（ステップＳ１３３）。具体的には、制御部２１の学習部２１１は、低画質化アクティビティ計測画像から生成した説明変数ベクトル、アクティビティ計測画像（血流拡散前）から生成した目的変数ベクトルを用いて、機械学習を行なう。そして、学習部２１１は、算出した統合推定モデルを学習結果記憶部２２に記録する。推定処理においては、予測部２１２は、アクティビティ計測画像を説明変数ベクトルとして、この統合推定モデルに適用する。そして、予測部２１２は、目的変数ベクトルとして、アクティビティ計測画像（血流拡散前）を算出する。これにより、低画質化の要因をまとめて学習し、効率的に高画質化を行なうことができる。
【００９９】
（第２の実施形態）
次に、図１５、１６に従って、線量分布予測システム、線量分布予測方法及び線量分布予測プログラムを具体化した第２の実施形態を説明する。なお、第２の実施形態は、第１の実施形態の機械学習処理、予測処理の手順を変更したのみの構成であるため、同様の部分についてはその詳細な説明を省略する。第１の実施形態では、粒子線の照射条件、照射対象の元素組成に基づいて、照射対象に対する粒子線の照射時の線量分布と、外部に放出される放射線の計測画像とをシミュレーションにより算出する。照射対象の患部配置、線量分布、外部に放出される放射線の計測画像（放射能分布に応じたアクティビティ計測画像）は相互に影響する。そこで、第２の実施形態では、アクティビティ計測画像から予測した患部配置を用いて、線量分布を算出する。
【０１００】
この場合には、支援サーバ２０に、カルテ情報記憶部を設ける。このカルテ情報記憶部には、患者コードに関連付けてカルテ情報が記録される。このカルテ情報には、患者の属性情報（年齢、性別等）、疾患情報（患部部位等）等の診療に関する情報が記録される。そして、支援サーバ２０の制御部２１は、治療支援処理を実行する。このために、制御部２１に、患部のシュリンク（患部の縮小）を含む患部状態を確認する治療支援部を設ける。
【０１０１】
図１５を用いて、治療支援処理を説明する。
まず、支援サーバ２０の制御部２１は、学習データの取得処理を実行する（ステップＳ１４１）。具体的には、制御部２１のシミュレータ２１０は、カルテ情報、陽子線照射条件、アクティビティ計測画像、線量分布、患部配置のデータセットを取得する。この場合、データセットのアクティビティ計測画像には、陽子線照射条件で陽子線の照射を行なった場合に実測したアクティビティ計測画像を用いる。なお、データセットのアクティビティ計測画像には、第１の実施形態と同様に、「解像度が低い実測したアクティビティ計測画像」から、「超解像モデルを用いて予測した高解像度のアクティビティ計測画像」を用いてもよい。また、データセットの線量分布には、第１の実施形態と同様に、陽子線照射条件に基づいて、ＣＴ画像から算出した組成密度分布に応じたペンシルビームを用いたシミュレーションにより算出した線量分布を用いる。データセットの患部配置には、ＣＴ画像において特定した患部領域を用いる。なお、学習データのデータセットについては、複数の患者のデータを用いることができる。
【０１０２】
更に、学習データのデータセット数を増やす場合には、シミュレーションを用いてもよい。
例えば、図１６に示すように、同一患者における異なる時期（ｔ１，ｔ２）のデータセット（６０１，６０２）が存在する場合には、これらを用いて新たなデータセット６０３を生成して、データセット数を増やしてもよい。
この場合、支援サーバ２０の制御部２１は、補間処理を実行する（ステップＳ１５１）。具体的には、データセット６０３のカルテ情報、陽子線照射条件については、制御部２１のシミュレータ２１０は、データセット（６０１，６０２）と同じデータを用いる。
また、患部配置については、制御部２１のシミュレータ２１０は、時期（ｔ１、ｔ２）のＣＴ画像の画像変換（例えば、モーフィング技術）を用いて、時期（ｔ１、ｔ２）の中間状態の画像（中間画像）を作成して用いる。ここで、時期ｔ２のＣＴ画像として、患部が完治している場合を想定して、ＣＴ画像において患部が消滅している画像を用いてもよい。
また、線量分布については、シミュレータ２１０は、陽子線照射条件に基づいて、中間画像から算出した組成密度分布に応じたペンシルビームを用いたシミュレーションにより算出した線量分布を用いる。
また、アクティビティ計測画像については、学習データ生成処理（図１１）を用いて、線量分布及び患部配置に基づいたシミュレーションにより作成する。なお、ここで、３次元アクティビティ分布の推定処理（図１２）を用いて生成した３次元アクティビティ分布（放射能分布）を用いてもよい。
上述した処理により、データセット６０３を作成し、学習データのデータセット数を増やすことができる。
【０１０３】
次に、図１５に示すように、支援サーバ２０の制御部２１は、機械学習処理を実行する（ステップＳ１４２）。具体的には、制御部２１の学習部２１１は、取得した学習データを用いて、患部配置を予測するための患部配置推定モデルを作成する。この場合には、学習データのデータセットの中で、アクティビティ分布（アクティビティ計測画像）、カルテ情報を入力層（説明変数ベクトル）として用い、患部配置（ＣＴ画像）を出力層（目的変数ベクトル）として機械学習を行なうことにより、患部配置推定モデルを生成し、学習結果記憶部２２に記録する。
【０１０４】
そして、支援サーバ２０の制御部２１は、予測処理を実行する（ステップＳ１４３）。具体的には、制御部２１の予測部２１２は、新たに予測用データを取得する。ここでは、予測用データとして、アクティビティ分布（アクティビティ計測画像）及びカルテ情報を取得する。アクティビティ計測画像には、患者の治療時に実測した画像を用いる。カルテ情報は、治療時の患者の患者コードを用いて、カルテ情報記憶部から取得した情報である。そして、予測部２１２は、アクティビティ計測画像、カルテ情報からなる予測用データを、患部配置推定モデルに適用して、患部配置を予測する。なお、学習処理において、アクティビティ計測画像の加工（例えば、３次元アクティビティ分布の推定処理（図１２））を行なった場合には、予測処理においても同様の加工を行なったアクティビティ分布を用いる。
【０１０５】
次に、支援サーバ２０の制御部２１は、線量分布の計算処理を実行する（ステップＳ１４４）。具体的には、制御部２１のシミュレータ２１０は、予測した患部配置を用いて、学習結果記憶部２２に記録された線量分布推定モデルにより線量分布を算出する。
【０１０６】
また、支援サーバ２０の制御部２１は、シュリンク確認処理を実行する（ステップＳ１４５）。具体的には、制御部２１の治療支援部は、患者コードを用いて、ＣＴ画像情報記憶部２３からＣＴ画像を取得する。そして、治療支援部は、ＣＴ画像の患部配置と、予測した患部配置とを比較し、比較結果をユーザ端末１０に出力する。この場合、異なる時期の患部配置と、現在の患部配置とにより、患部のシュリンク状態を確認することができる。
【０１０７】
以上、本実施形態によれば、上記（１−１）、（１−４）、（１−５）、（１−９）等に記載の効果と同様な効果に加えて、以下に示す効果を得ることができる。
（２−１）本実施形態では、支援サーバ２０の制御部２１は、学習データの取得処理を実行する（ステップＳ１４１）。この場合、支援サーバ２０の制御部２１は、補間処理を実行する（ステップＳ１５１）。これにより、学習データの不足を補うことができる。
（２−２）本実施形態では、支援サーバ２０の制御部２１は、機械学習処理（ステップＳ１４２）、予測処理（ステップＳ１４３）、線量分布の計算処理（ステップＳ１４４）を実行する。これにより、アクティビティ分布（放射能分布）、カルテ情報に基づいて、患部配置を予測し、予測した患部配置を用いて線量分布を算出することができる。
（２−３）本実施形態では、支援サーバ２０の制御部２１は、シュリンク確認処理を実行する（ステップＳ１４５）。これにより、患部の形状変化を把握し、治療の効果を確認することができる。また、患部の形状変化に応じて、粒子線の照射条件を変更することができる。
【０１０８】
（第３の実施形態）
次に、図１７に従って、線量分布予測システム、線量分布予測方法及び線量分布予測プログラムを具体化した第３の実施形態を説明する。なお、第３の実施形態は、第２の実施形態の機械学習処理、予測処理の手順を変更したのみの構成であるため、同様の部分についてはその詳細な説明を省略する。第２の実施形態では、アクティビティ画像から予測した患部領域を用いて、線量分布を算出する。第３の実施形態では、アクティビティ画像から患部領域や線量分布を予測する。この場合にも、上記第２実施形態と同様に、支援サーバ２０に、カルテ情報記憶部を設ける。
【０１０９】
図１７を用いて、予測処理を説明する。
支援サーバ２０の制御部２１は、ステップＳ１４１と同様に、学習データの取得処理を実行する（ステップＳ１６１）。
次に、支援サーバ２０の制御部２１は、機械学習処理を実行する（ステップＳ１６２）。具体的には、制御部２１の学習部２１１は、学習データを用いて、患部配置及び患部配置を予測するための推定モデルを作成する。この場合には、データセットの中で、アクティビティ分布（放射能分布に応じたアクティビティ計測画像）、カルテ情報、治療計画線量分布（第１の線量分布）を入力層（説明変数ベクトル）として用いる。この第１の線量分布は、治療計画の陽子線照射条件に基づいて算出される線量分布である。そして、患部配置（ＣＴ画像）、第２の線量分布を出力層（目的変数ベクトル）として用いる。この第２の線量分布は、患部配置を考慮して線量分布を想定した分布である。学習部２１１は、この入力層及び出力層を用いた機械学習を行なうことにより、患部配置及び線量分布の推定モデルを生成し、学習結果記憶部２２に記録する。
【０１１０】
次に、支援サーバ２０の制御部２１は、予測処理を実行する（ステップＳ１６３）。具体的には、制御部２１の予測部２１２は、新たに予測用データを取得する。ここでは、予測用データとして、アクティビティ計測画像、カルテ情報、治療計画線量分布（第１の線量分布）を取得する。アクティビティ計測画像は、患者の治療時に計測した画像である。カルテ情報は、治療時の患者の患者コードを用いて、カルテ情報記憶部から取得した情報である。第１の線量分布は、治療計画の陽子線照射条件に基づいて算出する。そして、予測部２１２は、アクティビティ計測画像、カルテ情報、第１の線量分布からなる予測用データを、患部配置及び線量分布の推定モデルに適用して、患部配置及び、この患部配置を考慮した第２の線量分布を予測する。
次に、支援サーバ２０の制御部２１は、ステップＳ１４５と同様に、シュリンク確認処理を実行する（ステップＳ１６４）。
【０１１１】
以上、本実施形態によれば、上記（２−１）、（２−３）に記載の効果と同様な効果に加えて、以下に示す効果を得ることができる。
（３−１）本実施形態では、アクティビティ分布（放射能分布）、カルテ情報、治療計画に基づいて、線量分布や患部配置を予測することができる。
【０１１２】
上記実施形態は、以下のように変更して実施することができる。上記実施形態及び以下の変更例は、技術的に矛盾しない範囲で互いに組み合わせて実施することができる。
・第２、第３の実施形態では、粒子線として陽子線を用いる。ここで、粒子線は陽子線に限定されるものではなく、例えば、炭素線等を用いることも可能である。
【０１１３】
次に、上記実施形態及び別例から把握できる技術的思想について、それらの効果とともに以下に追記する。
（ａ）放射能分布から患部配置を予測する推定モデルを記憶する記憶部に接続された制御部を備えた線量分布予測システムであって、
前記制御部が、
患部配置の患部画像、放射能分布を含む教師データを用いて、前記放射能分布から患部配置を予測する推定モデルを作成する学習処理と、
新たな放射能分布を取得した場合に、前記新たな放射能分布を、前記記憶部に記憶された前記推定モデルに適用して、患部配置を予測し、前記患部配置に応じた線量分布を算出する予測処理を実行することを特徴とする線量分布予測システム。
これにより、放射能分布から患部配置を予測し、この患部配置から線量分布を算出することができる。
【０１１４】
（ｂ）治療計画の線量分布及び放射能分布から、患部配置を考慮した線量分布を予測するための推定モデルを記憶する記憶部に接続された制御部を備えた線量分布予測システムであって、
前記制御部が、
患部画像、放射能分布、治療計画の線量分布を含む教師データを用いて、前記患部配置に応じた実際の線量分布を予測する推定モデルを作成し、記憶部に記録する学習処理と、
新たな放射能分布を取得した場合に、治療計画の線量分布及び前記新たな放射能分布を、前記記憶部に記憶された前記推定モデルに適用して、患部配置を考慮した線量分布を予測する予測処理を実行することを特徴とする線量分布予測システム。
これにより、放射能分布、治療計画の線量分布から、患部配置における線量分布を予測することができる。
【０１１５】
（ｃ）前記線量分布の予測時に患部配置を出力する処理を実行することを特徴とする（ｂ）に記載の線量分布予測システム
これにより、患部配置を予測することができる。
（ｄ）学習処理において、患者のカルテ情報を更に用いることを特徴とする（ａ）〜（ｃ）に記載の線量分布予測システム。
これにより、患者のカルテ情報を用いて、より的確な予測を行なうことができる。
（ｅ）治療前の患部配置に基づいて、治療後の患部配置を想定し、
前記推定モデルを用いて予測した患部配置と、前記想定した患部配置等が異なる場合には、注意喚起を出力することを特徴とする（ａ）〜（ｄ）に記載の線量分布予測システム。
これにより、患部配置における変化を確認することができる。
【符号の説明】
【０１１６】
１０…ユーザ端末、１５…画像取得部、２０…支援サーバ、２１…制御部、２１０…シミュレータ、２１１…学習部、２１２…予測部、２１３…血管情報取得部、２１４…残存影響処理部、２２…学習結果記憶部、２３…ＣＴ画像情報記憶部、２４…照射条件情報記憶部。
【図１】
2019030648.tif 000003
【図２】
2019030648.tif 000004
【図３】
2019030648.tif 000005
【図４】
2019030648.tif 000006
【図５】
2019030648.tif 000007
【図６】
2019030648.tif 000008
【図７】
2019030648.tif 000009
【図８】
2019030648.tif 000010
【図９】
2019030648.tif 000011
【図１０】
2019030648.tif 000012
【図１１】
2019030648.tif 000013
【図１２】
2019030648.tif 000014
【図１３】
2019030648.tif 000015
【図１４】
2019030648.tif 000016
【図１５】
2019030648.tif 000017
【図１６】
2019030648.tif 000018
【図１７】
2019030648.tif 000019
【図１８】
2019030648.tif 000020