放射線治療システム及び放射線治療プログラム
【課題】再構成処理において、散乱線が照射された領域についての条件を再構成に導入し、得られる散乱源分布の精度を向上することができる放射線治療システム等を提供すること。
【解決手段】治療放射線の線量を散乱線によりモニタする方式の放射線治療システムにおける再構成方法に関する。一般に治療用放射線は細い放射線錐を用いるが、従来の断層撮影法はこのような情報を利用しない。本発明では、反復再構成法を採用し、その反復中に測定した散乱線画像に基づいて、散乱線画像の値の小さい領域に対応する散乱源分布を減少させるステップを含める。また、散乱源分布が負の値の領域に対し散乱源分布を０に置き換えるステップを含める。特に散乱線画像が０に近い値を示す領域に対応する被検体の位置において、散乱源分布が０に近い値であるという条件を付加する。
【解決手段】治療放射線の線量を散乱線によりモニタする方式の放射線治療システムにおける再構成方法に関する。一般に治療用放射線は細い放射線錐を用いるが、従来の断層撮影法はこのような情報を利用しない。本発明では、反復再構成法を採用し、その反復中に測定した散乱線画像に基づいて、散乱線画像の値の小さい領域に対応する散乱源分布を減少させるステップを含める。また、散乱源分布が負の値の領域に対し散乱源分布を０に置き換えるステップを含める。特に散乱線画像が０に近い値を示す領域に対応する被検体の位置において、散乱源分布が０に近い値であるという条件を付加する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、放射線治療システムにおいて実際に照射した線量を計測し表示する機能を持つ放射線治療システムに関する。なかでも、散乱線画像から線量分布を３次元再構成する再構成処理に関する。
【背景技術】
【０００２】
Ｘ線外照射治療に代表される放射線治療では、治療前に、患者画像上で照射計画（病変部に対してどの方向から、どれだけの線量を照射するか）が立案され、これに基づいて患者への照射行われる。しかし、現在のところ実際に計画通りの位置、線量が患者に照射されているか否かを確認する手段がなく、病変部への過少照射や正常組織への過剰照射が起こっても気づかれないのが現状である。照射前にファントムとＸ線検出器を用いて、計画通りの照射が行えることが確認されることもあるが、簡便に持ち運びができ、自由に位置を調整できるファントムと異なり、患者を寝台上の、照射計画どおりの位置に置くことは困難であり、これらの照射前確認は、患者への計画通りの照射を完全に保証するものではない。
【０００３】
放射線照射装置で被検体の一部（患者の患部）を照射中に、放射線通過領域にて散乱するＸ線を散乱線検出器にて多方向で測定し、断層撮影により散乱源の分布を３次元再構成する方法が提案されている。
【０００４】
放射線治療では通常、幹部の形状に合わせた細い放射線ビームを用いる。散乱線は放射線が照射された領域にて主に発生するが、その領域は被検体の中のきわめて一部分である。一般のＣＴ再構成または断層撮影の再構成処理は解の存在領域などを仮定せずに再構成を行うが、放射線治療の照射領域のモニタに応用する場合、散乱源の場所に関して一部のみに存在するという条件を付加することが可能である。一般のＣＴ再構成または断層撮影の再構成処理を用いるより解の精度を向上できると考えられる。
【特許文献１】特開平５−１４６４２６号公報 この特許文献が開示する技術は、Ｘ線被写体の散乱Ｘ線を検出し、被写体の断層像を得るものである。ペンシル状ビームを走査することにより被写体の３次元散乱線像を再構成して得ることが特徴である。すなわち、本技術はペンシル状ビームのみを想定しており、Ｘ線治療で用いられる、有限の幅をもったビームが通過した領域の散乱像（治療ビームによる線量の空間分布）を得るものではない。また、エネルギーの高い治療ビーム(数MeV)の被写体内での散乱は前方散乱が優位となるため、入射Ｘ線方向に検出器を配置すると散乱線と透過線の区別が難しく、散乱線の検出に補正処理を必須としている。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
放射線治療では通常、幹部の形状に合わせた細い放射線ビームを用いる。散乱線は放射線が照射された領域にて主に発生するが、その領域は被検体の中のきわめて一部分である。一般のＣＴ再構成または断層撮影の再構成処理は解の存在領域などを仮定せずに再構成を行うが、放射線治療の照射領域を映像化する場合、散乱源の領域に関して一部のみに存在するという条件を付加することが可能である。従って、一般のＣＴ再構成または断層撮影の再構成処理を用いるよりも、放射線治療の照射領域の映像化固有の再構成手法を確立することができれば、画像診断の精度を向上できると考えられる。
【０００６】
本発明は、上記事情を鑑みてなされたもので、再構成処理において、散乱線が照射された領域についての条件を再構成に導入して、得られる散乱源分布の精度を向上することができる放射線治療システム及び放射線治療プログラムを提供することを目的としている。
【課題を解決するための手段】
【０００７】
本発明は、上記目的を達成するため、次のような手段を講じている。
【０００８】
請求項１に記載の発明は、被検体に対して治療用放射線ビームを照射する照射手段と、前記治療用放射線ビームに基づいて発生する前記被検体内からの散乱線を検出し散乱線データを発生する検出手段と、前記検出された各散乱線データ上における、前記被検体内の散乱源存在領域を決定する決定手段と、前記散乱源存在領域が決定された前記各散乱線データを入力とし所定の条件に従う再構成処理を繰り返し実行することで再構成画像を生成し、当該再構成画像を用いて、前記被検体内における散乱線発生密度の三次元的分布を示す散乱線ボリュームデータを生成する再構成手段と、前記散乱線ボリュームデータを、吸収された放射線量の三次元分布を示す吸収線量ボリュームデータに変換する変換手段と、前記吸収線量ボリュームデータを用いて、前記被検体内における吸収線量画像を生成する画像生成手段と、前記吸収線量画像を表示する表示手段と、を具備することを特徴とする放射線治療システムである。
【０００９】
請求項６に記載の発明は、コンピュータに、被検体に対して照射された治療用放射線ビームに基づいて発生する前記被検体内からの散乱線を検出させる検出機能と、前記検出された散乱線に基づいて複数の散乱線データを生成させるデータ生成機能と、前記検出された各散乱線データ上における、前記被検体内の散乱源存在領域を決定させ決定機能と、前記散乱源存在領域が決定された前記各散乱線データを入力とし所定の条件に従う再構成処理を繰り返し実行することで再構成画像を生成させ、当該再構成画像を用いて、前記被検体内における散乱線発生密度の三次元的分布を示す散乱線ボリュームデータを生成させる再構成機能と、前記散乱線ボリュームデータを、吸収された放射線量の三次元分布を示す吸収線量ボリュームデータに変換させる変換機能と、前記吸収線量ボリュームデータを用いて、前記被検体内における吸収線量画像を生成させる画像生成機能と、前記吸収線量画像を表示させる表示機能と、を実現させることを特徴とする放射線治療プログラムである。
【発明の効果】
【００１０】
以上本発明によれば、再構成処理において、散乱線が照射された領域についての条件を再構成に導入して、得られる散乱源分布の精度を向上することができる放射線治療システム及び放射線治療プログラムを実現することができる。
【発明を実施するための最良の形態】
【００１１】
（第１実施形態）
［原理と方法］
本実施形態に係る放射線治療システムは、被検体に対して照射した放射線に基づく当該被検体からの散乱線を計測し、これに基づいて被検体のどの部位に、どれだけの線量が照射されたかを客観的に示す情報を取得するものである。その原理と方法は、次の様である。
【００１２】
図１は、本放射線治療システムの治療用放射線に基づく被検体からの散乱線計測の原理、方法を説明するための図である。
【００１３】
外照射Ｘ線照射による治療効果は、主として患者体内で起こるＸ線の散乱によってもたらされる。すなわち、治療Ｘ線ビームが患者体内の電子によって散乱される際、エネルギーを受け取った電子は組織内を飛行したのち、停止する。このとき、電子は停止するまでに、組織内の分子をラジカル化し、細胞内のDNAに損傷を与える。そして、損傷を受け、修復することができなかった細胞は最終的に死に至る。これがＸ線照射による治療効果である。反跳電子が多く発生すればするほど組織を構成する細胞が死に至る確率が高くなるため、治療効果は、散乱反応が起こる回数に比例する。
【００１４】
上述から、組織内で起こった散乱の回数が分かれば、治療効果（=組織がどれだけ損傷を受けたか）を知ることができる。そして起こった散乱の回数は、散乱線の数を測定することで知ることができる。散乱されたＸ線の多くは、電子に進行方向を変えられた後、患者体外に出てくるため、患者体外に設置したＸ線検出器で測定することができる。
【００１５】
本実施形態の第１の実施例に係る放射線治療システムでは、治療Ｘ線ビームに対して特定の角度をなす位置にコリメータを備えた検出器を設置し、その方向に来た散乱線のみを選択的に検出する。コンプトン散乱で、どの角度に、どれだけＸ線が散乱されるかは理論的に分かるため、ある角度での散乱線を検出できれば、他の角度への散乱線の数も推定できる。さらに、患者体内の、散乱の起こった場所の分布を３次元的に得るために、照射中に検出器を回転させ、すべての方向から散乱線の測定を行う（例えば、図４参照）。その後、再構成処理を行い、被検体内部の散乱線の発生分布を３次元的に画像化する。
【００１６】
また、本実施形態の第２の実施例に係る放射線治療システムでは、治療Ｘ線ビームに対して所定の角度（散乱角）をなす位置にコリメータを備えた検出器を設置し、その方向に来た散乱線のみを選択的に検出し、この検出を照射部から照射される治療用Ｘ線ビームの軸と検出器の検出面とのなす角を維持しつつ治療用Ｘ線ビームと検出面とを移動させながら実行することで、被検体内の３次元領域をスキャンする。得られた所定の散乱角に関する３次元散乱線データを用いて、散乱線ボリュームデータを再構成すると共に、当該散乱線ボリュームデータを吸収された放射線量の３次元分布を示す吸収線量ボリュームデータに変換し、吸収線量画像を生成する。
【００１７】
［構成］
図２は、本実施形態に係る放射線治療システム１のブロック構成図を示している。同図に示すように、本放射線治療システム１は、放射線照射システム２、散乱線検出システム３、データ取得制御部４、データ処理システム５、表示部６、記憶部７、操作部８、ネットワークＩ／Ｆ９を具備している。放射線照射システム２及び散乱線検出システム３は架台（ガントリ）に設置され、架台を移動、回転させることで、被検体に対して任意の位置に配置することができる。また、データ取得制御部４、データ処理システム５、表示部６、記憶部７、操作部８、ネットワークＩ／Ｆ９は、例えば放射線治療システム１の本体（筐体）に設置される。
【００１８】
［放射線照射システム］
放射線照射システム２は、電力供給部２０１、照射部２０３、タイミング制御部２０５、ガントリ制御部２０７を有している。
【００１９】
電力供給部２０１は、データ取得制御部４からの制御に従って照射部２０３に電力を供給する。
【００２０】
照射部２０３は、例えば線形加速器（ライナック）等の機構を有する放射線照射装置である。当該照射部２０３では、加速管の一端に設けられた電子銃により、陰極から放射された熱電子は数１００keVになるまで加速される。次に、クライストロンで発生したマイクロ波は導波管を使って加速管まで導かれ、そこでこの熱電子は数MeVのエネルギーに達するまで加速される。この加速された熱電子は磁石によってその方向を変えられ、透過型ターゲットに衝突する。このとき制動放射により、数MeVのエネルギーのＸ線が発生する。照射部２０３は、コリメータによってこのＸ線を所定の形状（例えば、円錐形状、或いは薄い平面形状）に成形し、寝台上に配置された被検体の三次元領域に照射する。
【００２１】
タイミング制御部２０５は、データ取得制御部４からの制御に従って所定のタイミングで照射部２０３に電力が供給されるように、電力供給部２０１を制御する。
【００２２】
ガントリ制御部２０７は、例えば操作部８やデータ取得制御部４からの制御指示に従って、ガントリの移動位置・回転位置を制御する。
【００２３】
［散乱線検出システム］
散乱線検出システム３は、検出器３０１、コリメータ３０３、移動機構部３０５、位置検出部３０７を有している。
【００２４】
検出器３０１は、数1００keVのＸ線を検出できる半導体検出器や、イメージング・プレート等であり、被検体に対して照射した放射線に基づく当該被検体からの散乱線を検出する。この検出器の好ましいサイズ、照射ビーム軸に対する配置角度、画素数等については、後述する。
【００２５】
コリメータ３０３は、特定の方向に来た散乱線のみを選択的に検出するための絞り装置である。
【００２６】
移動機構部３０５は、照射部２０３の照射ビーム軸に対する検出器３０１の検出面の角度（すなわち、照射ビーム軸と検出器３０１の検出面の法線との角度）、放射線ビーム軸を中心とした検出器３０１の回転角、被検体と検出器３０１の検出面との距離等を制御するために、検出器３０１の位置や角度を移動させるための移動機構部である。
【００２７】
位置検出部３０７は、検出器３０１の位置を検出するためのエンコーダである。
【００２８】
［データ取得制御部］
データ取得制御部４は、放射線治療時における散乱線計測に関する総合的な制御を行う。例えば、データ取得制御部４は、放射線照射システム２のタイミング制御部２０５からの信号を得て、散乱線検出システム３に対して散乱線計測開始トリガーや検出データの伝送トリガーを送信する等、放射線照射、散乱線計測、データ処理、画像表示、ネットワーク通信等について、本放射線治療システム１を静的又は動的に制御する。また、データ取得制御部４は、必要に応じて、ネットワークを介して放射線治療計画装置から受け取った治療計画に基づいて、各照射の照射時間に合わせてスキャン時間を最適化する。
【００２９】
［データ処理システム］
データ処理システム５は、補正処理部５０１、散乱源存在領域決定部５０２、再構成処理部５０３、変換処理部５０５、データ処理部５０７を有している。
【００３０】
補正処理部５０１は、必要に応じてデータのキャリブレーション処理やノイズを除去するための補正処理等を行う。当該補正処理５０１が実行する補正処理の内容については、後で詳しく説明する。
【００３１】
散乱源存在領域決定部５０２は、散乱線検出器の各位置での測定結果を用いたバックプロジェクション処理により、被検体内に局在する散乱線源に対応する領域を決定する。この散乱源存在領域を決定する処理の内容については、後で詳しく説明する。
【００３２】
再構成処理部５０３は、散乱線検出システム３において検出された散乱線画像データと各散乱線画像データを検出した位置を示す位置情報とを用いて画像再構成処理を実行し、散乱イベント回数（散乱発生回数）の密度の三次元的分布を示す散乱線ボリュームデータを取得する。再構成法式としては、例えば、コリメータの方向がスキャン軸と直交していればＣＴの再構成手法を、一方直交していなければ、断層撮影の再構成手法を用いる。特に、本再構成処理部５０３は、照射領域以外に分布する散乱源を０と仮定することで散乱性存在領域のみを対象とし、プロジェクション画像が測定画像と一致するまで解を修正しながらバックプロジェクションとプロジェクションを反復する再構成処理を実行する。この反復再構成処理の内容については、後で詳しく説明する。
【００３３】
変換処理部５０５は、画像再構成処理によって得られた三次元画像データを、吸収された放射線量（吸収線量）の三次元分布を示す吸収線量ボリュームデータに変換する。
【００３４】
［表示部、記憶部、操作部、ネットワークＩ／Ｆ］
表示部６は、ＬＣＤ等のディスプレイで構成される。表示部６は、上記データ処理部５０７の各種モジュールにより出力されるデータをもとに、例えば、吸収線量画像を計画画像や照射直前、照射中に得た画像とフュージョンして表示を行う。
【００３５】
記憶部７は、照射する放射線ビームの軸を中心として検出器３０１を回転させながら散乱線データを取得（スキャン）するための所定のスキャンシーケンス７０１、補正処理、画像再構成処理、変換処理、表示処理等の実行、および治療計画を当該システムで表示、編集するための制御プログラム７０２や、当該放射線治療システム１によって取得された散乱線ボリュームデータ７０３、吸収線量ボリュームデータ７０４、Ｘ線コンピュータ断層撮影装置等の他のモダリティによって取得された形態画像データ７０５等を記憶する。当該記憶部７に記憶されているこれらのデータは、ネットワークＩ／Ｆ９０を経由して外部装置へ転送することも可能となっている。
【００３６】
操作部８は、オペレータからの各種指示、条件、関心領域（ＲＯＩ）の設定指示、種々の画質条件設定指示等を装置本体１１にとりこむための各種スイッチ、ボタン、トラックボール１３ｓ、マウス１３ｃ、キーボード１３ｄ等を有している。
【００３７】
ネットワークＩ／Ｆ９は、当該放射線治療システム１によって得られた吸収線量画像データ等をネットワーク経由で他の装置に転送し、また、例えば放射線治療計画装置において作成された治療計画等をネットワーク経由で取得する。
【００３８】
（散乱線ボリュームデータ等の生成方法）
（第１の実施例）
次に、第１の実施例に係る放射線治療システム１を用いた散乱線ボリュームデータ等の生成方法について説明する。本実施形態に係る放射線治療システムでは、治療Ｘ線ビームに対して特定の角度をなす位置にコリメータを備えた検出器を設置し、その方向に来た散乱線のみを選択的に検出する。さらに、患者体内の、散乱の起こった場所の分布を３次元的に得るために、照射中に検出器を回転させ、すべての方向から散乱線の測定を行う（例えば、図６参照）。その後、再構成処理を行い、被検体内部の散乱線の発生分布を３次元的に画像化する。
【００３９】
図３は、本実施形態に係る吸収線量画像データの生成処理を含む放射線治療時における処理の流れを示したフローチャートである。以下、各ステップの処理内容について説明する。
【００４０】
[被検体の配置等：ステップＳ１ａ]
まず、データ取得制御部４は、例えばネットワークを介して当該被検体に関する治療計画情報を取得し、表示部６に表示する。術者は、表示された治療計画に従って寝台上に被検体を配置すると共に、操作部８を介して、放射線照射時間の設定、散乱線計測を行う回転角度の設定、スキャンシーケンスの選択等を行う（ステップＳ１ａ）。なお、放射線照射時間の設定等については、取得した治療計画情報に基づいて、自動的に行うようにしてもよい。
【００４１】
[放射線照射／多方向における散乱線画像データの取得：ステップＳ２ａ]
図４は、本放射線治療システム１の散乱線の測定形態を示した図である。同図に示すように、放射線照射システム２は被検体に対して、三次元領域を照射するための治療用放射線を所定のタイミングで発生する。また、散乱線検出システム２は、当該照射放射線に基づいて被検体外に出てくる散乱線を照射される放射線ビームの軸を中心とした複数の回転角において検出する（ステップＳ２ａ）。例えば、ある１つの方向から３分間照射が行える場合、１方向につき１０秒ずつ、１８方向のデータを収集する。このとき、１８方向はビーム軸を中心として等角度間隔であることが好ましい。検出器３０３が各方向で検出した散乱線のカウント数及び位置検出部３０７で計測した散乱線検出時における検出器３０３の位置情報は、データ処理システム５に伝送される。
【００４２】
なお、本実施形態では、検出器３０１の配置角度を、散乱角θが１２０°≦θ≦１６５°の範囲のいずれか（例えば、１５５°）である後方散乱線を検出するように、検出器３０１の配置角度を設定するものとする。
【００４３】
また、上記の例において、例えば２Gyの照射が３方向から行われる場合、１方向あたりのカウント数は、1.24×1０5× 1/3となり、約4×1０4[counts/cm2]である。１方向あたり１８０秒で照射されるとして１０秒間測定すると、4×1０4×1０/18０=2×1０3[counts/cm2]となるが、Ｓ／Ｎ比に問題はない。
【００４４】
また、散乱線の検出は、少なくとも２つ以上の方向が必要であるが、現実にはできる限り多くの方向において検出することが好ましい。また、各検出位置は、照射ビームの軸を中心として等角度間隔に配置されていることが好ましい。
【００４５】
[前処理（補正処理等）：ステップＳ３ａ]
収集されたデータは、検出器設置角度方向に散乱されたＸ線のみカウントしている。しかし実際には、Ｘ線はあらゆる方向への散乱が起こっている。データ処理システム５の補正処理部５０１は、検出器のカウント値を補正し、所定の計算しきに従って、すべての方向への散乱数を取得する（ステップＳ３ａ）。
【００４６】
[画像再構成処理：ステップＳ４ａ]
次に、データ処理システム５の画像再構成処理部５０３は、多方向の投影データを用いて画像再構成処理を実行し、散乱線ボリュームデータを取得する（ステップＳ４ａ）。このとき、検出器３０１の回転軸とコリメータの方向が直交しており、１８０度（＋α）以上の角度範囲で画像を撮影する場合はＣＴの再構成方法を用いればよいが、その他の場合は断層撮影の再構成方法を用いる。断層撮影の手法として、例えば投影画像にフィルタ処理を適用した後バックプロジェクション処理を行うfiltered backprojection法を用いる。filterの構成方法としては古典的なShepp-Logan filterや、特願２００６−２８４３２５, 特願２００７−２６９４４７に開示されているフィルタを用いる。特に、特願２００６−２８４３２５, 特願２００７−２６９４４７に記載されている方法を用いれば、物理的意味が明確な散乱源分布画像を生成することができる。
【００４７】
検出器画像にフィルタ処理を施し、バックプロジェクションを行って得られる画像は、単位体積あたりの散乱線発生密度（単位体積あたりの散乱回数）である。上記の再構成処理の全ステップ（各種補正処理、フィルタ処理、バックプロジェクション処理）をとおして、治療用放射線が被検体を通過する近傍での散乱線発生密度の３次元分布（散乱線ボリュームデータ）を取得することができる。
【００４８】
なお、本実施形態に係る放射線治療システムでは、本ステップにおいて後述する条件付き反復再構成処理が実行されることになる。
【００４９】
[変換処理：ステップＳ５ａ]
次に、データ処理システム５の変換処理部５０７は、ボクセル（voxel）ごとに算出された単位体積あたりの散乱回数nを、吸収線量に換算することで、散乱線ボリュームデータを吸収された放射線量（吸収線量）の三次元分布を示す吸収線量ボリュームデータに変換する（ステップＳ５ａ）。
【００５０】
［吸収線量画像データの生成／画像データの表示：ステップＳ６ａ、Ｓ７ａ］
次に、画像処理部５０７は、吸収線量ボリュームデータ等を用いて、被検体の所定部位に関する吸収された放射線量（吸収線量）の分布を示す吸収線量画像データを生成し、例えばフュージョン表示するためにＣＴ画像と合成する（ステップＳ６ａ）。表示部６は、所定の形態にて吸収線量画像を表示する（ステップＳ７ａ）。
【００５１】
（第２の実施例）
次に、第２の実施例に係る放射線治療システム１を用いた散乱線ボリュームデータ等の生成方法について説明する。本実施形態に係る放射線治療システムでは、治療Ｘ線ビームに対して所定の角度（散乱角）をなす位置にコリメータを備えた検出器を設置し、その方向に来た散乱線のみを選択的に検出し、この検出を照射部から照射される治療用Ｘ線ビームの軸と検出器の検出面とのなす角を維持しつつ治療用Ｘ線ビームと検出面とを移動させながら実行することで、被検体内の３次元領域をスキャンする。得られた所定の散乱角に関する３次元散乱線データを用いて、散乱線ボリュームデータを再構成すると共に、当該散乱線ボリュームデータを吸収された放射線量の３次元分布を示す吸収線量ボリュームデータに変換し、吸収線量画像を生成する。
【００５２】
図５は、本実施形態に係る吸収線量画像データの生成処理を含む放射線治療時における処理の流れを示したフローチャートである。以下、各ステップの処理内容について説明する。
【００５３】
[被検体の配置等：ステップＳ１ｂ]
まず、第１の実施形態と同様に、被検体の配置等が実行される（ステップＳ１ｂ）。
【００５４】
[放射線照射（散乱線データの取得）：ステップＳ２]
図６は、本放射線治療システム１の散乱線の測定形態の一例を示した図である。同図に示すように、放射線照射システム２は、被検体に対して薄い平面状に整形されたＸ線ビームＢ２を所定のタイミングで照射し、放射線検出システム２は、当該照射放射線に基づいて被検体外に出てくる所定の散乱角の散乱線を検出する。また、データ取得制御部４は、照射部２０３から照射される治療用のＸ線ビームＢ２の軸と検出器３０１の視線方向とのなす角を維持しながらＸ線ビームＢ２による励起断面を移動させ、当該被検体内の３次元領域を走査（スキャン）するように、ガントリ制御部２０７或いは移動機構部３０５を制御する（ステップＳ２）。この治療用のＸ線ビームＢ２を用いた３次元領域のスキャンにより、Ｘ線ビームＢ２の平面に対応する複数の二次元散乱線データからなる３次元散乱線データが取得される。
【００５５】
なお、図６は、散乱線の測定形態の一例である。従って、本実施形態に係る散乱線の測定形態は、当該例に拘泥されない。例えば、図７に示すように、検出器３０１の検出面（及びコリメータ３０３の開口面）を、治療用放射線ビームの照射方向に対する検出面のなす角度を一定に保ちながら、治療用放射線ビームの軸の位置の移動に連動して移動させることによっても、複数の二次元散乱線データからなる３次元散乱線データを取得することができる。
【００５６】
[前処理（補正処理等）：ステップＳ３ｂ]
次に、データ処理システム５の補正処理部５０１は、減弱補正を含む前処理を実行し、投影データを取得する（ステップＳ３）。ここで、減弱補正とは、治療用放射線や散乱線が被検体内を伝播することに起因する信号減弱に関する補正処理である。
【００５７】
[画像再構成処理：ステップＳ４ｂ]
次に、データ処理システム５の画像再構成処理部５０３は、取得された投影データを用いて画像再構成処理を実行し、散乱線ボリュームデータを取得する（ステップＳ４）。なお、本実施形態に係る放射線治療システムでは、本ステップにおいて後述する条件付き反復再構成処理が実行されることになる。
【００５８】
[変換処理：ステップＳ５ｂ]
次に、データ処理システム５の変換処理部５０７は、第１の実施形態と同様に、散乱線ボリュームデータを吸収された放射線量（吸収線量）の３次元分布を示す吸収線量ボリュームデータに変換する（ステップＳ５）。
【００５９】
［吸収線量画像データの生成／画像データの表示：ステップＳ６ｂ、Ｓ７ｂ］
次に、画像処理部５０７は、吸収線量ボリュームデータ等を用いて、被検体の所定部位に関する吸収された放射線量（吸収線量）の分布を示す吸収線量画像データを生成し、例えばフュージョン表示するためにＣＴ画像と合成する（ステップＳ６ｂ）。表示部６は、所定の形態にて吸収線量画像を表示する（ステップＳ７ｂ）。
【００６０】
（条件付き反復再構成機能）
次に、本放射線治療システム１が有する条件付き反復再構成機能について説明する。一般のＸ線断層撮影やＸ線ＣＴと異なり、放射線治療システムのＸ線照射領域（散乱源分布）を画像化する場合は放射線を照射した領域が限局していることがわかっている。本条件付き反復再構成機能では、まず、図８（ａ）に示すように、散乱線検出器の各位置での測定結果を用いたバックプロジェクション処理により散乱線源の存在領域を決定する。次に、図８（ｂ）に示すように、再構成処理において、照射領域以外に分布する散乱源を０と仮定することで散乱性存在領域のみを対象とし、プロジェクション画像が測定画像と一致するまで解を修正しながらバックプロジェクションとプロジェクションを反復する再構成処理を実行し、被検体において局在する放射線の被照射領域における散乱源の空間分布（又は、吸収線量の空間分布）を好適に映像化するものである。なお、この条件付き反復再構成機能は、例えば図３のステップＳ４ａ、図４のステップＳ４ｂにおいて実行される。
【００６１】
図９は、本条件付き反復再構成機能に従う処理（条件付き反復再構成処理）の流れを示したフローチャートである。同図に示すように、本条件付き反復再構成処理は、大きく２つのステップからなる。以下、各ステップにおける処理の内容について説明する。
【００６２】
［散乱源存在可能領域の決定処理：ステップＳ４１］
図１０は、散乱源存在可能領域の決定処理の流れを示したフローチャートである。
【００６３】
まず、散乱性存在領域決定部５０２は、Np個のプロジェクションの各散乱線画像において、図１１（ａ）に示すように測定値が０に近い領域とそれ以外の領域に分割する（ステップＳ４１１）。領域分割には閾値法を用いればよいが、閾値を決定するためには例えば（ヒストグラムのピーク値＋定数)などを用いることができる。閾値を用いて2値化した後、ノイズの影響を除くため孤立した小さな前景領域を除去する処理を実施する。この時点で、プロジェクションごとに、２値化投影画像が得られる。2値化投影画像は散乱線画像が０に近い領域では値が０で、それ以外の領域では値が1となる。
【００６４】
次に、散乱性存在領域決定部５０２は、図１１（ｂ）に示すように2値化投影画像をバックプロジェクションする（ステップＳ４１２）。散乱性存在領域決定部５０２は、バックプロジェクションの結果得られた画像にて値がプロジェクション数と等しい領域を、散乱源存在可能領域と決定する（ステップＳ４１３）。なお、本実施形態の例では、図１１（ｃ）に示すように、値が4の領域が散乱源存在可能領域Ｓとして決定されることになる。
【００６５】
［反復再構成処理：ステップＳ４２］
図１２は、散乱源存在可能領域の決定処理の流れを示したフローチャートである。
【００６６】
まず、再構成処理部５０３は、暫定的な解画像を設定する（ステップＳ４２１）。本実施形態では、暫定解画像の初期値として、ステップＳ４１において得られた解を用いる。
【００６７】
次に、再構成処理部５０３は、暫定解画像の散乱源存在可能領域以外の領域を０に置き換えた後（ステップＳ４２２）、プロジェクション処理によりプロジェクション画像を生成し、測定値としての投影画像との差を求め（ステップＳ４２３）、この差を投影画像として用いて再度再構成する（ステップＳ４２４、ステップＳ４２５）。
【００６８】
次に、再構成処理部５０３は、差分の再構成画像を、1回目の暫定解画像（０に置き換えたもの）から引き算し（ステップＳ４２６）、当該再構成処理の反復回数が所定回数に到達したか否かを判定する。所定回数に到達していないと判定した場合には、再構成処理部５０３は、ステップＳ４２５で得られた最新の投影画像を暫定的な解画像として、ステップＳ４２２〜ステップＳ４２６間での処理を繰り返すことその暫定解画像を反復修正し最終結果に向けて収束させていく。一方、所定回数に到達したと判定した場合には、再構成処理部５０３は、反復再構成処理を終了する。
【００６９】
（効果）
以上述べた構成によれば、次の効果を実現することができる。
【００７０】
本放射線治療システムでは、散乱線検出器の各位置での測定結果を用いたバックプロジェクション処理により、被検体内の散乱線源存在領域を決定し、当該散乱源存在可能領域外の暫定解を０に置き換えて初回の再構成を実行する。これにより、誤差分の投影画像は、１回目の反復の再構成でその大部分が散乱源存在可能領域に再構成され、散乱源存在可能領域の暫定解画像の値は大きく減少する。また、次の再構成処理においては、前回の再構成処理で得られた暫定解画像を入力画像とし、その散乱源存在可能領域外の暫定解を０に置き換えて再び再構成処理を実行し、これを繰り返す。この様な反復再構成処理により、散乱源存在可能領域に対応する信号を忠実に映像化することができると共に、散乱源存在可能領域に対応しない信号を除外することができ、断層撮像特有のアーチファクトを低減させ、画質を向上させることができる。
【００７１】
なお、本条件付き反復再構成処理は、少なくとも２解の繰り返しでも再構成誤差を改善するために大きな効果がある。さらに反復回数を増やせばより効果があるが、最大でも1０回の反復すれば十分である。
【００７２】
（第２の実施形態）
第１の実施形態では、散乱源が存在すべきではない領域に対し散乱源の分布が０であるという条件を導入することで、解の精度の向上を図るものであった。これに対し、本実施形態では、この散乱源の分布が０という条件をより一般化した方法について説明する。散乱源が存在すべきではない領域に対し散乱源の分布が０という条件の導入は2つの効果がある。
【００７３】
・従来の再構成方法では、散乱源が本来０になるべき領域で０より大きい解が得られることがある。本来０であるはずの領域なので、この領域について散乱源＝０という条件を付加することで解の精度の向上を図ることができる。
【００７４】
・従来の再構成方法では、散乱源が小さな値を持つ領域で負の値の解が得られることがある。特に、散乱源が本来０になるべき領域では負の値の解が得られやすい。原理的に散乱源の値は負の値にはならないので、この領域について散乱源=０という条件を付加することで解の精度の向上図ることができる。
【００７５】
本実施形態では、上記２つの置き換えを下記のように一般化し、解の上限および下限についての条件に従う反復再構成処理を行う。
【００７６】
・解を投影した画像が、再構成処理において入力となる投影画像（再構成入力投影画像）を超えない。
【００７７】
・解は負の値にはならない（非負条件）。
【００７８】
上の条件を考慮した上で、「解を投影した画像が再構成入力投影画像にできるだけ近くなる」ように、解を選ぶという方針をとる。なお、第１の実施形態では、散乱源が本来０になるべき領域についてだけ条件を付加することができたが、本第２の実施形態では、上の２つの条件で解の上限と下限を定めているため、全領域について条件を付加することが可能になっている。
【００７９】
再構成処理部５０３において実行される具体的な処理の内容は、以下の様である。
【００８０】
［基本となる方法］
再構成空間の画像ｆのプロジェクション処理をＡｆと表記する。この処理の結果は投影空間のベクトルである。また、投影画像ｇのバックプロジェクション処理をＡＴｇと表記する。この処理の結果は、再構成空間のベクトルとなる。
【００８１】
まず、本手法のベースとなる従来方法について説明する。次の式（１）を満たすｆを求める場合を考える。
【数１】
【００８２】
Ｐｆは解ｆのペナルティーベクトルであり、望ましくないｆに対してｆＴＰｆが大きな値を持つ。常にｆＴＰｆ≧０である。上記式（１）より、次の式（２）が得られる。
【数２】
【００８３】
また、式（２）の右辺を導関数とする関数ｑ（ｆ）を次の式（３）によって定義する。
【数３】
【００８４】
すなわち、式（１）を解くには、式（３）を目的関数として最小化すればよい。式（３）式を最小化するには、共役勾配法などの最小化手法が利用できる。共役勾配法を含む多くの最小化手法では評価関数ｑの計算とその導関数の計算が必要であるが、それらの処理要素であるＡｆおよびＡＴｇは前述のとおりプロジェクション処理とバックプロジェクション処理なので、それらプロジェクション処理とバックプロジェクション処理を実装することが必要である。参考のため、ｑ（ｆ）の導関数を次の式（４）に示す。
【数４】
【００８５】
ある方向ｈに修正しようとしたとき（具体的には解ｆにλｈを加える）、方向ｈにどれだけの距離を進めば最小になるかを決定することが必要になる。これはｑ（ｆ＋λｈ）を最小にするλを決めることであり、次の式（５）を解くことで求めることができる。
【数５】
【００８６】
以下に、式（５）を解くための計算を示す。
【数６】
【００８７】
以上を具体的に実装した従来法の処理手順を図１３に示す。説明上の例として最急降下法の例を示してしるが、実際には、共役勾配法等、より優れた方法を用いるべきである。
【００８８】
［非負条件の導入］
本実施形態では、解は負の値にならないという条件を導入するため、図１３に示した処理手順を、図１４に示す処理手順に変更する（ステップ９ｂが新たに追加されたステップである）。
【００８９】
［解を投影した画像が、再構成入力投影画像を超えないという条件の導入］
式（３）に変えて、次の式（６）に示す目的関数を用いる。
【数７】
【００９０】
ここで、Ｐｏｓｉｔｉｖｅ（ｖ）は要素νｉの値が正の値のときｉ番目の要素の値がνｉ、０以下の値のときｉ番目の要素の値が０になる関数である。Ｎｅｇａｔｉｖｅ（ｖ）は逆に要素νｉの値が０以下のときｉ番目の要素の値がνｉ、正の値のときｉ番目の要素の値が０になる関数である。式（６）は、αｇ＝αｓ＝１のとき、式（３）と等価である。典型的にはαｇ＞αｓとなる値を用いる。
【００９１】
式（６）の導関数として次の式（７）を用いる。
【数８】
【００９２】
式（６）のＡｆ−ｇは暫定解の投影と入力投影画像の差である。αｇ＞αｓとなる値を用いたとき、入力投影画像を暫定解の投影が超えている場合、より大きな係数がかかるため、目的関数（６）はより大きな値になる。
【００９３】
図１５は、αｇ＝２、αｓ＝１の例である。最適化の過程では目的関数を小さくする解を選ぶので、結果として、暫定解の投影が入力投影画像を大きく超えない解が選ばれるようになる。
【００９４】
なお、処理手順は、図１６に示すようになる。この処理手順の特徴的な点は、ステップ6bでの導関数の算出において、暫定解ｈｉの投影と入力投影画像を比較し、Ａｆ−ｇの各要素の符号に応じ、正の値の方が大きい係数をＡｆ−ｇの各要素に乗じる処理が追加される点である。ステップ６ｂの変更に伴い、ステップ８で算出されるλが不正確なものになるので、ステップ８ｂを追加してλをより正確な値に修正することが望ましい。
【００９５】
（効果）
以上述べた本実施形態に係る条件付き反復再構成処理によれば、解が非負であるという下限に関する条件と、解を投影した画像が、再構成入力投影画像を超えないという上限に関する条件の導入が可能になり、解を制限することができる。例えば、特に散乱線画像が０に近い値を示す領域に対応する被検体の位置において、散乱源分布が０に近い値であるという条件が付加される。従って、第２の実施形態は、第１の実施形態と同様の効果を持っている。さらに、散乱線画像が０でないある一定値を示す領域においても、散乱源分布を投影した値がそれより小さい値になるという条件が付加される。このことは、第２の実施形態では第１の実施形態より強い制限事項を再構成に課している。従って、第２の実施形態によれば、第１の実施形態に比して、再構成の精度をさらに向上できるという効果がある。
【００９６】
（第３の実施形態）
第２の実施形態では、解全体について非負の条件を課すことができたが、このような方法では非負条件が有効に機能しないことがある。例えば、図１７（ａ）に示すように散乱源が一部だけ大きな値を持つような場合、図１７（ｂ）に示すようにその周囲では負の値を持つ誤った解が得られやすい。この場合、非負条件により解が改善することができる。しかし、図１７（ａ）に対して図１７（ｃ）に示すように一様な散乱源分布が加わる場合、図１７（ｂ）と同様の再構成方式では図１７（ｄ）に示すように解が負の値にならない。つまり、解に図１７（ｂ）と同じ間違いが含まれていても、一様成分が加わるために解が負の値にならず間違いを修正することができない。
【００９７】
そこで、本実施形態では、モデル化された複数の分布の合成で解（散乱源の分布）を表わす。すなわち、本来は一つである解を複数のモデル化された解に分解して考える。そのうえで、分解した解の各々について第２の実施形態の条件を課すことで、例えば図１７（ｃ）に示すような一様な散乱分布が存在する場合であっても、好適な画像再構成を実現するものである。
【００９８】
解の合成は、次の式（８）で表すことができる。
【数９】
【００９９】
上記式（８）では、解がＭ個の関数の合成で表されているが、そのうち、先頭のＲ個のみで近似した合成式が、次の式（９）で算出できるようにしておく。
【数１０】
【０１００】
分解された解分布ｆｊのうち小さなｊに対応する分布は広域的な分布とし、ｊが大きくなるほど局所分布となるようにｆｊを定める。とくに、ｆ０は一様分布である。ｆｊの各要素は全て０以上の値を持つように設定するのが良い。
【０１０１】
目的関数は、式（６）に式（８）を代入することで、次の式（１０）の様になる。
【数１１】
【０１０２】
また、導関数として次の式（１１）を用いる。
【数１２】
【０１０３】
ＦＲｍＲの値はＲ個の解分布ｆｊの合成なので、目的関数および導関数はＲの値により計算処理が異なる。Ｒを決めれば、第２の実施形態の方法でｍＲを求め、ＦＲｍＲを計算することで解分布を求めることができる。その際、第２の実施形態の記載中のＡは本実施形態のＡＦＲに置き換える必要がある。
【０１０４】
本実施形態では、Ｒ＝１から開始して、Ｒを徐々に大きくしながら、図１９に示すように反復的に再構成処理を実行する。投影画像ｇは最初再構成入力画像の値を設定しておくが、反復ごとに解ｆｋの投影Ａｆｋを引き算していく。反復中のステップｋにおいては、それより前の全ステップで差し引かれた投影画像からの再構成を行う。ｆｊは小さなｊに対する分布がより広い分布を持っており、反復が進むほどＲが増加していくので、ステップが進むほど細かい構造が再構成されるようになる。一例として最初は一様分布が推定され、その後、局所の分布が推定される。２ステップ目以降は、投影画像から一様分布に対する投影画像が除去されているので、最初に示した概念図で説明したように以降の各ステップで局所の分布を推定する際も非負条件が有効に働く。
【０１０５】
以上述べた本実施形態に係る再構成処理によれば、図１８（ｃ）の例（図１７（ｃ）の例と同じ）の場合、解を図１８（ｃ１）に示す一様成分と図１８（ｃ２）に示す局所成分に分解し、各々に非負の条件を課して再構成する。このとき、非負条件が加えられるので、解に含まれる誤差がそれぞれ図１８（ｅ１）、図１８（ｅ２）に示すように低減される。最終的に得られる解はそれぞれを合成であるが、もちろん、最終的な解に含まれる誤差も図１８（ｅ）に示すように低減されている。
【０１０６】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。具体的な変形例としては、例えば次のようなものがある。
【０１０７】
例えば、本実施形態に係る各機能は、当該処理を実行するプログラムをワークステーション等のコンピュータにインストールし、これらをメモリ上で展開することによっても実現することができる。このとき、コンピュータに当該手法を実行させることのできるプログラムは、磁気ディスク（フロッピー（登録商標）ディスク、ハードディスクなど）、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤなど）、半導体メモリなどの記録媒体に格納して頒布することも可能である。
【０１０８】
また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【産業上の利用可能性】
【０１０９】
以上本発明によれば、再構成処理において、散乱線が照射された領域についての条件を再構成に導入して、得られる散乱源分布の精度を向上することができる放射線治療システム及び放射線治療プログラムを実現することができる。
【図面の簡単な説明】
【０１１０】
【図１】図１は、本放射線治療用線量分布測定装置の治療用放射線に基づく被検体からの散乱線計測の原理、方法を説明するための図である。
【図２】図２は、本実施形態に係る放射線治療用線量分布測定装置１のブロック構成図を示している。
【図３】図３は、吸収線量画像データの生成処理を含む放射線治療時における処理の流れを示したフローチャートである。
【図４】図４は、本放射線治療システム１の散乱線の測定形態を示した図である。
【図５】図５は、本実施形態に係る吸収線量画像データの生成処理を含む放射線治療時における処理の流れを示したフローチャートである。
【図６】図６は、本放射線治療システム１の散乱線の測定形態の一例を示した図である。
【図７】図７は、本放射線治療システム１の散乱線の測定形態の他の例を示した図である。
【図８】図８（ａ）、（ｂ）は、条件付き反復再構成機能の概念を説明するための図である。
【図９】図９は、本条件付き反復再構成機能に従う処理の流れを示したフローチャートである。
【図１０】図１０は、散乱源存在可能領域の決定処理の流れを示したフローチャートである。
【図１１】図１１（ａ）、（ｂ）、（ｃ）は、散乱源存在可能領域の決定処理を説明するための図である。
【図１２】図１２は、散乱源存在可能領域の決定処理の流れを示したフローチャートである。
【図１３】図１３は、最急降下法の処理手順を示したフローチャートである。
【図１４】図１４は、非負条件が導入された最急降下法の処理手順を示したフローチャートである。
【図１５】図１５は、αｇ＝２、αｓ＝１の場合の入力投影画像及び暫定解の投影画像の例を示した図である。
【図１６】図１６は、解を投影した画像が、再構成入力投影画像を超えないという条件が導入された場合の最急降下法の処理手順を示したフローチャートである。
【図１７】図１７は、第３の実施形態に係る再構成処理の概念を説明するための図である。
【図１８】図１８は、第３の実施形態に係る再構成処理の概念を説明するための図である。
【図１９】図１９は、第３の実施形態に係る再構成処理の流れを示したフローチャートである。
【符号の説明】
【０１１１】
１…放射線治療用線量分布測定装置、２…放射線照射システム、３…散乱線検出システム、４…データ取得制御部、５…データ処理システム、６…表示部、７…記憶部、８…操作部、９…ネットワークＩ／Ｆ、２０１…電力供給部、２０３…照射部、２０５…タイミング制御部、２０７…ガントリ制御部、３０１…検出器、３０３…コリメータ、３０５…移動機構部、３０７…位置検出部、５０１…補正処理部、５０２…散乱源存在領域決定部、５０３…再構成処理部、５０５…変換処理部、５０７…画像処理部
【技術分野】
【０００１】
本発明は、放射線治療システムにおいて実際に照射した線量を計測し表示する機能を持つ放射線治療システムに関する。なかでも、散乱線画像から線量分布を３次元再構成する再構成処理に関する。
【背景技術】
【０００２】
Ｘ線外照射治療に代表される放射線治療では、治療前に、患者画像上で照射計画（病変部に対してどの方向から、どれだけの線量を照射するか）が立案され、これに基づいて患者への照射行われる。しかし、現在のところ実際に計画通りの位置、線量が患者に照射されているか否かを確認する手段がなく、病変部への過少照射や正常組織への過剰照射が起こっても気づかれないのが現状である。照射前にファントムとＸ線検出器を用いて、計画通りの照射が行えることが確認されることもあるが、簡便に持ち運びができ、自由に位置を調整できるファントムと異なり、患者を寝台上の、照射計画どおりの位置に置くことは困難であり、これらの照射前確認は、患者への計画通りの照射を完全に保証するものではない。
【０００３】
放射線照射装置で被検体の一部（患者の患部）を照射中に、放射線通過領域にて散乱するＸ線を散乱線検出器にて多方向で測定し、断層撮影により散乱源の分布を３次元再構成する方法が提案されている。
【０００４】
放射線治療では通常、幹部の形状に合わせた細い放射線ビームを用いる。散乱線は放射線が照射された領域にて主に発生するが、その領域は被検体の中のきわめて一部分である。一般のＣＴ再構成または断層撮影の再構成処理は解の存在領域などを仮定せずに再構成を行うが、放射線治療の照射領域のモニタに応用する場合、散乱源の場所に関して一部のみに存在するという条件を付加することが可能である。一般のＣＴ再構成または断層撮影の再構成処理を用いるより解の精度を向上できると考えられる。
【特許文献１】特開平５−１４６４２６号公報 この特許文献が開示する技術は、Ｘ線被写体の散乱Ｘ線を検出し、被写体の断層像を得るものである。ペンシル状ビームを走査することにより被写体の３次元散乱線像を再構成して得ることが特徴である。すなわち、本技術はペンシル状ビームのみを想定しており、Ｘ線治療で用いられる、有限の幅をもったビームが通過した領域の散乱像（治療ビームによる線量の空間分布）を得るものではない。また、エネルギーの高い治療ビーム(数MeV)の被写体内での散乱は前方散乱が優位となるため、入射Ｘ線方向に検出器を配置すると散乱線と透過線の区別が難しく、散乱線の検出に補正処理を必須としている。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
放射線治療では通常、幹部の形状に合わせた細い放射線ビームを用いる。散乱線は放射線が照射された領域にて主に発生するが、その領域は被検体の中のきわめて一部分である。一般のＣＴ再構成または断層撮影の再構成処理は解の存在領域などを仮定せずに再構成を行うが、放射線治療の照射領域を映像化する場合、散乱源の領域に関して一部のみに存在するという条件を付加することが可能である。従って、一般のＣＴ再構成または断層撮影の再構成処理を用いるよりも、放射線治療の照射領域の映像化固有の再構成手法を確立することができれば、画像診断の精度を向上できると考えられる。
【０００６】
本発明は、上記事情を鑑みてなされたもので、再構成処理において、散乱線が照射された領域についての条件を再構成に導入して、得られる散乱源分布の精度を向上することができる放射線治療システム及び放射線治療プログラムを提供することを目的としている。
【課題を解決するための手段】
【０００７】
本発明は、上記目的を達成するため、次のような手段を講じている。
【０００８】
請求項１に記載の発明は、被検体に対して治療用放射線ビームを照射する照射手段と、前記治療用放射線ビームに基づいて発生する前記被検体内からの散乱線を検出し散乱線データを発生する検出手段と、前記検出された各散乱線データ上における、前記被検体内の散乱源存在領域を決定する決定手段と、前記散乱源存在領域が決定された前記各散乱線データを入力とし所定の条件に従う再構成処理を繰り返し実行することで再構成画像を生成し、当該再構成画像を用いて、前記被検体内における散乱線発生密度の三次元的分布を示す散乱線ボリュームデータを生成する再構成手段と、前記散乱線ボリュームデータを、吸収された放射線量の三次元分布を示す吸収線量ボリュームデータに変換する変換手段と、前記吸収線量ボリュームデータを用いて、前記被検体内における吸収線量画像を生成する画像生成手段と、前記吸収線量画像を表示する表示手段と、を具備することを特徴とする放射線治療システムである。
【０００９】
請求項６に記載の発明は、コンピュータに、被検体に対して照射された治療用放射線ビームに基づいて発生する前記被検体内からの散乱線を検出させる検出機能と、前記検出された散乱線に基づいて複数の散乱線データを生成させるデータ生成機能と、前記検出された各散乱線データ上における、前記被検体内の散乱源存在領域を決定させ決定機能と、前記散乱源存在領域が決定された前記各散乱線データを入力とし所定の条件に従う再構成処理を繰り返し実行することで再構成画像を生成させ、当該再構成画像を用いて、前記被検体内における散乱線発生密度の三次元的分布を示す散乱線ボリュームデータを生成させる再構成機能と、前記散乱線ボリュームデータを、吸収された放射線量の三次元分布を示す吸収線量ボリュームデータに変換させる変換機能と、前記吸収線量ボリュームデータを用いて、前記被検体内における吸収線量画像を生成させる画像生成機能と、前記吸収線量画像を表示させる表示機能と、を実現させることを特徴とする放射線治療プログラムである。
【発明の効果】
【００１０】
以上本発明によれば、再構成処理において、散乱線が照射された領域についての条件を再構成に導入して、得られる散乱源分布の精度を向上することができる放射線治療システム及び放射線治療プログラムを実現することができる。
【発明を実施するための最良の形態】
【００１１】
（第１実施形態）
［原理と方法］
本実施形態に係る放射線治療システムは、被検体に対して照射した放射線に基づく当該被検体からの散乱線を計測し、これに基づいて被検体のどの部位に、どれだけの線量が照射されたかを客観的に示す情報を取得するものである。その原理と方法は、次の様である。
【００１２】
図１は、本放射線治療システムの治療用放射線に基づく被検体からの散乱線計測の原理、方法を説明するための図である。
【００１３】
外照射Ｘ線照射による治療効果は、主として患者体内で起こるＸ線の散乱によってもたらされる。すなわち、治療Ｘ線ビームが患者体内の電子によって散乱される際、エネルギーを受け取った電子は組織内を飛行したのち、停止する。このとき、電子は停止するまでに、組織内の分子をラジカル化し、細胞内のDNAに損傷を与える。そして、損傷を受け、修復することができなかった細胞は最終的に死に至る。これがＸ線照射による治療効果である。反跳電子が多く発生すればするほど組織を構成する細胞が死に至る確率が高くなるため、治療効果は、散乱反応が起こる回数に比例する。
【００１４】
上述から、組織内で起こった散乱の回数が分かれば、治療効果（=組織がどれだけ損傷を受けたか）を知ることができる。そして起こった散乱の回数は、散乱線の数を測定することで知ることができる。散乱されたＸ線の多くは、電子に進行方向を変えられた後、患者体外に出てくるため、患者体外に設置したＸ線検出器で測定することができる。
【００１５】
本実施形態の第１の実施例に係る放射線治療システムでは、治療Ｘ線ビームに対して特定の角度をなす位置にコリメータを備えた検出器を設置し、その方向に来た散乱線のみを選択的に検出する。コンプトン散乱で、どの角度に、どれだけＸ線が散乱されるかは理論的に分かるため、ある角度での散乱線を検出できれば、他の角度への散乱線の数も推定できる。さらに、患者体内の、散乱の起こった場所の分布を３次元的に得るために、照射中に検出器を回転させ、すべての方向から散乱線の測定を行う（例えば、図４参照）。その後、再構成処理を行い、被検体内部の散乱線の発生分布を３次元的に画像化する。
【００１６】
また、本実施形態の第２の実施例に係る放射線治療システムでは、治療Ｘ線ビームに対して所定の角度（散乱角）をなす位置にコリメータを備えた検出器を設置し、その方向に来た散乱線のみを選択的に検出し、この検出を照射部から照射される治療用Ｘ線ビームの軸と検出器の検出面とのなす角を維持しつつ治療用Ｘ線ビームと検出面とを移動させながら実行することで、被検体内の３次元領域をスキャンする。得られた所定の散乱角に関する３次元散乱線データを用いて、散乱線ボリュームデータを再構成すると共に、当該散乱線ボリュームデータを吸収された放射線量の３次元分布を示す吸収線量ボリュームデータに変換し、吸収線量画像を生成する。
【００１７】
［構成］
図２は、本実施形態に係る放射線治療システム１のブロック構成図を示している。同図に示すように、本放射線治療システム１は、放射線照射システム２、散乱線検出システム３、データ取得制御部４、データ処理システム５、表示部６、記憶部７、操作部８、ネットワークＩ／Ｆ９を具備している。放射線照射システム２及び散乱線検出システム３は架台（ガントリ）に設置され、架台を移動、回転させることで、被検体に対して任意の位置に配置することができる。また、データ取得制御部４、データ処理システム５、表示部６、記憶部７、操作部８、ネットワークＩ／Ｆ９は、例えば放射線治療システム１の本体（筐体）に設置される。
【００１８】
［放射線照射システム］
放射線照射システム２は、電力供給部２０１、照射部２０３、タイミング制御部２０５、ガントリ制御部２０７を有している。
【００１９】
電力供給部２０１は、データ取得制御部４からの制御に従って照射部２０３に電力を供給する。
【００２０】
照射部２０３は、例えば線形加速器（ライナック）等の機構を有する放射線照射装置である。当該照射部２０３では、加速管の一端に設けられた電子銃により、陰極から放射された熱電子は数１００keVになるまで加速される。次に、クライストロンで発生したマイクロ波は導波管を使って加速管まで導かれ、そこでこの熱電子は数MeVのエネルギーに達するまで加速される。この加速された熱電子は磁石によってその方向を変えられ、透過型ターゲットに衝突する。このとき制動放射により、数MeVのエネルギーのＸ線が発生する。照射部２０３は、コリメータによってこのＸ線を所定の形状（例えば、円錐形状、或いは薄い平面形状）に成形し、寝台上に配置された被検体の三次元領域に照射する。
【００２１】
タイミング制御部２０５は、データ取得制御部４からの制御に従って所定のタイミングで照射部２０３に電力が供給されるように、電力供給部２０１を制御する。
【００２２】
ガントリ制御部２０７は、例えば操作部８やデータ取得制御部４からの制御指示に従って、ガントリの移動位置・回転位置を制御する。
【００２３】
［散乱線検出システム］
散乱線検出システム３は、検出器３０１、コリメータ３０３、移動機構部３０５、位置検出部３０７を有している。
【００２４】
検出器３０１は、数1００keVのＸ線を検出できる半導体検出器や、イメージング・プレート等であり、被検体に対して照射した放射線に基づく当該被検体からの散乱線を検出する。この検出器の好ましいサイズ、照射ビーム軸に対する配置角度、画素数等については、後述する。
【００２５】
コリメータ３０３は、特定の方向に来た散乱線のみを選択的に検出するための絞り装置である。
【００２６】
移動機構部３０５は、照射部２０３の照射ビーム軸に対する検出器３０１の検出面の角度（すなわち、照射ビーム軸と検出器３０１の検出面の法線との角度）、放射線ビーム軸を中心とした検出器３０１の回転角、被検体と検出器３０１の検出面との距離等を制御するために、検出器３０１の位置や角度を移動させるための移動機構部である。
【００２７】
位置検出部３０７は、検出器３０１の位置を検出するためのエンコーダである。
【００２８】
［データ取得制御部］
データ取得制御部４は、放射線治療時における散乱線計測に関する総合的な制御を行う。例えば、データ取得制御部４は、放射線照射システム２のタイミング制御部２０５からの信号を得て、散乱線検出システム３に対して散乱線計測開始トリガーや検出データの伝送トリガーを送信する等、放射線照射、散乱線計測、データ処理、画像表示、ネットワーク通信等について、本放射線治療システム１を静的又は動的に制御する。また、データ取得制御部４は、必要に応じて、ネットワークを介して放射線治療計画装置から受け取った治療計画に基づいて、各照射の照射時間に合わせてスキャン時間を最適化する。
【００２９】
［データ処理システム］
データ処理システム５は、補正処理部５０１、散乱源存在領域決定部５０２、再構成処理部５０３、変換処理部５０５、データ処理部５０７を有している。
【００３０】
補正処理部５０１は、必要に応じてデータのキャリブレーション処理やノイズを除去するための補正処理等を行う。当該補正処理５０１が実行する補正処理の内容については、後で詳しく説明する。
【００３１】
散乱源存在領域決定部５０２は、散乱線検出器の各位置での測定結果を用いたバックプロジェクション処理により、被検体内に局在する散乱線源に対応する領域を決定する。この散乱源存在領域を決定する処理の内容については、後で詳しく説明する。
【００３２】
再構成処理部５０３は、散乱線検出システム３において検出された散乱線画像データと各散乱線画像データを検出した位置を示す位置情報とを用いて画像再構成処理を実行し、散乱イベント回数（散乱発生回数）の密度の三次元的分布を示す散乱線ボリュームデータを取得する。再構成法式としては、例えば、コリメータの方向がスキャン軸と直交していればＣＴの再構成手法を、一方直交していなければ、断層撮影の再構成手法を用いる。特に、本再構成処理部５０３は、照射領域以外に分布する散乱源を０と仮定することで散乱性存在領域のみを対象とし、プロジェクション画像が測定画像と一致するまで解を修正しながらバックプロジェクションとプロジェクションを反復する再構成処理を実行する。この反復再構成処理の内容については、後で詳しく説明する。
【００３３】
変換処理部５０５は、画像再構成処理によって得られた三次元画像データを、吸収された放射線量（吸収線量）の三次元分布を示す吸収線量ボリュームデータに変換する。
【００３４】
［表示部、記憶部、操作部、ネットワークＩ／Ｆ］
表示部６は、ＬＣＤ等のディスプレイで構成される。表示部６は、上記データ処理部５０７の各種モジュールにより出力されるデータをもとに、例えば、吸収線量画像を計画画像や照射直前、照射中に得た画像とフュージョンして表示を行う。
【００３５】
記憶部７は、照射する放射線ビームの軸を中心として検出器３０１を回転させながら散乱線データを取得（スキャン）するための所定のスキャンシーケンス７０１、補正処理、画像再構成処理、変換処理、表示処理等の実行、および治療計画を当該システムで表示、編集するための制御プログラム７０２や、当該放射線治療システム１によって取得された散乱線ボリュームデータ７０３、吸収線量ボリュームデータ７０４、Ｘ線コンピュータ断層撮影装置等の他のモダリティによって取得された形態画像データ７０５等を記憶する。当該記憶部７に記憶されているこれらのデータは、ネットワークＩ／Ｆ９０を経由して外部装置へ転送することも可能となっている。
【００３６】
操作部８は、オペレータからの各種指示、条件、関心領域（ＲＯＩ）の設定指示、種々の画質条件設定指示等を装置本体１１にとりこむための各種スイッチ、ボタン、トラックボール１３ｓ、マウス１３ｃ、キーボード１３ｄ等を有している。
【００３７】
ネットワークＩ／Ｆ９は、当該放射線治療システム１によって得られた吸収線量画像データ等をネットワーク経由で他の装置に転送し、また、例えば放射線治療計画装置において作成された治療計画等をネットワーク経由で取得する。
【００３８】
（散乱線ボリュームデータ等の生成方法）
（第１の実施例）
次に、第１の実施例に係る放射線治療システム１を用いた散乱線ボリュームデータ等の生成方法について説明する。本実施形態に係る放射線治療システムでは、治療Ｘ線ビームに対して特定の角度をなす位置にコリメータを備えた検出器を設置し、その方向に来た散乱線のみを選択的に検出する。さらに、患者体内の、散乱の起こった場所の分布を３次元的に得るために、照射中に検出器を回転させ、すべての方向から散乱線の測定を行う（例えば、図６参照）。その後、再構成処理を行い、被検体内部の散乱線の発生分布を３次元的に画像化する。
【００３９】
図３は、本実施形態に係る吸収線量画像データの生成処理を含む放射線治療時における処理の流れを示したフローチャートである。以下、各ステップの処理内容について説明する。
【００４０】
[被検体の配置等：ステップＳ１ａ]
まず、データ取得制御部４は、例えばネットワークを介して当該被検体に関する治療計画情報を取得し、表示部６に表示する。術者は、表示された治療計画に従って寝台上に被検体を配置すると共に、操作部８を介して、放射線照射時間の設定、散乱線計測を行う回転角度の設定、スキャンシーケンスの選択等を行う（ステップＳ１ａ）。なお、放射線照射時間の設定等については、取得した治療計画情報に基づいて、自動的に行うようにしてもよい。
【００４１】
[放射線照射／多方向における散乱線画像データの取得：ステップＳ２ａ]
図４は、本放射線治療システム１の散乱線の測定形態を示した図である。同図に示すように、放射線照射システム２は被検体に対して、三次元領域を照射するための治療用放射線を所定のタイミングで発生する。また、散乱線検出システム２は、当該照射放射線に基づいて被検体外に出てくる散乱線を照射される放射線ビームの軸を中心とした複数の回転角において検出する（ステップＳ２ａ）。例えば、ある１つの方向から３分間照射が行える場合、１方向につき１０秒ずつ、１８方向のデータを収集する。このとき、１８方向はビーム軸を中心として等角度間隔であることが好ましい。検出器３０３が各方向で検出した散乱線のカウント数及び位置検出部３０７で計測した散乱線検出時における検出器３０３の位置情報は、データ処理システム５に伝送される。
【００４２】
なお、本実施形態では、検出器３０１の配置角度を、散乱角θが１２０°≦θ≦１６５°の範囲のいずれか（例えば、１５５°）である後方散乱線を検出するように、検出器３０１の配置角度を設定するものとする。
【００４３】
また、上記の例において、例えば２Gyの照射が３方向から行われる場合、１方向あたりのカウント数は、1.24×1０5× 1/3となり、約4×1０4[counts/cm2]である。１方向あたり１８０秒で照射されるとして１０秒間測定すると、4×1０4×1０/18０=2×1０3[counts/cm2]となるが、Ｓ／Ｎ比に問題はない。
【００４４】
また、散乱線の検出は、少なくとも２つ以上の方向が必要であるが、現実にはできる限り多くの方向において検出することが好ましい。また、各検出位置は、照射ビームの軸を中心として等角度間隔に配置されていることが好ましい。
【００４５】
[前処理（補正処理等）：ステップＳ３ａ]
収集されたデータは、検出器設置角度方向に散乱されたＸ線のみカウントしている。しかし実際には、Ｘ線はあらゆる方向への散乱が起こっている。データ処理システム５の補正処理部５０１は、検出器のカウント値を補正し、所定の計算しきに従って、すべての方向への散乱数を取得する（ステップＳ３ａ）。
【００４６】
[画像再構成処理：ステップＳ４ａ]
次に、データ処理システム５の画像再構成処理部５０３は、多方向の投影データを用いて画像再構成処理を実行し、散乱線ボリュームデータを取得する（ステップＳ４ａ）。このとき、検出器３０１の回転軸とコリメータの方向が直交しており、１８０度（＋α）以上の角度範囲で画像を撮影する場合はＣＴの再構成方法を用いればよいが、その他の場合は断層撮影の再構成方法を用いる。断層撮影の手法として、例えば投影画像にフィルタ処理を適用した後バックプロジェクション処理を行うfiltered backprojection法を用いる。filterの構成方法としては古典的なShepp-Logan filterや、特願２００６−２８４３２５, 特願２００７−２６９４４７に開示されているフィルタを用いる。特に、特願２００６−２８４３２５, 特願２００７−２６９４４７に記載されている方法を用いれば、物理的意味が明確な散乱源分布画像を生成することができる。
【００４７】
検出器画像にフィルタ処理を施し、バックプロジェクションを行って得られる画像は、単位体積あたりの散乱線発生密度（単位体積あたりの散乱回数）である。上記の再構成処理の全ステップ（各種補正処理、フィルタ処理、バックプロジェクション処理）をとおして、治療用放射線が被検体を通過する近傍での散乱線発生密度の３次元分布（散乱線ボリュームデータ）を取得することができる。
【００４８】
なお、本実施形態に係る放射線治療システムでは、本ステップにおいて後述する条件付き反復再構成処理が実行されることになる。
【００４９】
[変換処理：ステップＳ５ａ]
次に、データ処理システム５の変換処理部５０７は、ボクセル（voxel）ごとに算出された単位体積あたりの散乱回数nを、吸収線量に換算することで、散乱線ボリュームデータを吸収された放射線量（吸収線量）の三次元分布を示す吸収線量ボリュームデータに変換する（ステップＳ５ａ）。
【００５０】
［吸収線量画像データの生成／画像データの表示：ステップＳ６ａ、Ｓ７ａ］
次に、画像処理部５０７は、吸収線量ボリュームデータ等を用いて、被検体の所定部位に関する吸収された放射線量（吸収線量）の分布を示す吸収線量画像データを生成し、例えばフュージョン表示するためにＣＴ画像と合成する（ステップＳ６ａ）。表示部６は、所定の形態にて吸収線量画像を表示する（ステップＳ７ａ）。
【００５１】
（第２の実施例）
次に、第２の実施例に係る放射線治療システム１を用いた散乱線ボリュームデータ等の生成方法について説明する。本実施形態に係る放射線治療システムでは、治療Ｘ線ビームに対して所定の角度（散乱角）をなす位置にコリメータを備えた検出器を設置し、その方向に来た散乱線のみを選択的に検出し、この検出を照射部から照射される治療用Ｘ線ビームの軸と検出器の検出面とのなす角を維持しつつ治療用Ｘ線ビームと検出面とを移動させながら実行することで、被検体内の３次元領域をスキャンする。得られた所定の散乱角に関する３次元散乱線データを用いて、散乱線ボリュームデータを再構成すると共に、当該散乱線ボリュームデータを吸収された放射線量の３次元分布を示す吸収線量ボリュームデータに変換し、吸収線量画像を生成する。
【００５２】
図５は、本実施形態に係る吸収線量画像データの生成処理を含む放射線治療時における処理の流れを示したフローチャートである。以下、各ステップの処理内容について説明する。
【００５３】
[被検体の配置等：ステップＳ１ｂ]
まず、第１の実施形態と同様に、被検体の配置等が実行される（ステップＳ１ｂ）。
【００５４】
[放射線照射（散乱線データの取得）：ステップＳ２]
図６は、本放射線治療システム１の散乱線の測定形態の一例を示した図である。同図に示すように、放射線照射システム２は、被検体に対して薄い平面状に整形されたＸ線ビームＢ２を所定のタイミングで照射し、放射線検出システム２は、当該照射放射線に基づいて被検体外に出てくる所定の散乱角の散乱線を検出する。また、データ取得制御部４は、照射部２０３から照射される治療用のＸ線ビームＢ２の軸と検出器３０１の視線方向とのなす角を維持しながらＸ線ビームＢ２による励起断面を移動させ、当該被検体内の３次元領域を走査（スキャン）するように、ガントリ制御部２０７或いは移動機構部３０５を制御する（ステップＳ２）。この治療用のＸ線ビームＢ２を用いた３次元領域のスキャンにより、Ｘ線ビームＢ２の平面に対応する複数の二次元散乱線データからなる３次元散乱線データが取得される。
【００５５】
なお、図６は、散乱線の測定形態の一例である。従って、本実施形態に係る散乱線の測定形態は、当該例に拘泥されない。例えば、図７に示すように、検出器３０１の検出面（及びコリメータ３０３の開口面）を、治療用放射線ビームの照射方向に対する検出面のなす角度を一定に保ちながら、治療用放射線ビームの軸の位置の移動に連動して移動させることによっても、複数の二次元散乱線データからなる３次元散乱線データを取得することができる。
【００５６】
[前処理（補正処理等）：ステップＳ３ｂ]
次に、データ処理システム５の補正処理部５０１は、減弱補正を含む前処理を実行し、投影データを取得する（ステップＳ３）。ここで、減弱補正とは、治療用放射線や散乱線が被検体内を伝播することに起因する信号減弱に関する補正処理である。
【００５７】
[画像再構成処理：ステップＳ４ｂ]
次に、データ処理システム５の画像再構成処理部５０３は、取得された投影データを用いて画像再構成処理を実行し、散乱線ボリュームデータを取得する（ステップＳ４）。なお、本実施形態に係る放射線治療システムでは、本ステップにおいて後述する条件付き反復再構成処理が実行されることになる。
【００５８】
[変換処理：ステップＳ５ｂ]
次に、データ処理システム５の変換処理部５０７は、第１の実施形態と同様に、散乱線ボリュームデータを吸収された放射線量（吸収線量）の３次元分布を示す吸収線量ボリュームデータに変換する（ステップＳ５）。
【００５９】
［吸収線量画像データの生成／画像データの表示：ステップＳ６ｂ、Ｓ７ｂ］
次に、画像処理部５０７は、吸収線量ボリュームデータ等を用いて、被検体の所定部位に関する吸収された放射線量（吸収線量）の分布を示す吸収線量画像データを生成し、例えばフュージョン表示するためにＣＴ画像と合成する（ステップＳ６ｂ）。表示部６は、所定の形態にて吸収線量画像を表示する（ステップＳ７ｂ）。
【００６０】
（条件付き反復再構成機能）
次に、本放射線治療システム１が有する条件付き反復再構成機能について説明する。一般のＸ線断層撮影やＸ線ＣＴと異なり、放射線治療システムのＸ線照射領域（散乱源分布）を画像化する場合は放射線を照射した領域が限局していることがわかっている。本条件付き反復再構成機能では、まず、図８（ａ）に示すように、散乱線検出器の各位置での測定結果を用いたバックプロジェクション処理により散乱線源の存在領域を決定する。次に、図８（ｂ）に示すように、再構成処理において、照射領域以外に分布する散乱源を０と仮定することで散乱性存在領域のみを対象とし、プロジェクション画像が測定画像と一致するまで解を修正しながらバックプロジェクションとプロジェクションを反復する再構成処理を実行し、被検体において局在する放射線の被照射領域における散乱源の空間分布（又は、吸収線量の空間分布）を好適に映像化するものである。なお、この条件付き反復再構成機能は、例えば図３のステップＳ４ａ、図４のステップＳ４ｂにおいて実行される。
【００６１】
図９は、本条件付き反復再構成機能に従う処理（条件付き反復再構成処理）の流れを示したフローチャートである。同図に示すように、本条件付き反復再構成処理は、大きく２つのステップからなる。以下、各ステップにおける処理の内容について説明する。
【００６２】
［散乱源存在可能領域の決定処理：ステップＳ４１］
図１０は、散乱源存在可能領域の決定処理の流れを示したフローチャートである。
【００６３】
まず、散乱性存在領域決定部５０２は、Np個のプロジェクションの各散乱線画像において、図１１（ａ）に示すように測定値が０に近い領域とそれ以外の領域に分割する（ステップＳ４１１）。領域分割には閾値法を用いればよいが、閾値を決定するためには例えば（ヒストグラムのピーク値＋定数)などを用いることができる。閾値を用いて2値化した後、ノイズの影響を除くため孤立した小さな前景領域を除去する処理を実施する。この時点で、プロジェクションごとに、２値化投影画像が得られる。2値化投影画像は散乱線画像が０に近い領域では値が０で、それ以外の領域では値が1となる。
【００６４】
次に、散乱性存在領域決定部５０２は、図１１（ｂ）に示すように2値化投影画像をバックプロジェクションする（ステップＳ４１２）。散乱性存在領域決定部５０２は、バックプロジェクションの結果得られた画像にて値がプロジェクション数と等しい領域を、散乱源存在可能領域と決定する（ステップＳ４１３）。なお、本実施形態の例では、図１１（ｃ）に示すように、値が4の領域が散乱源存在可能領域Ｓとして決定されることになる。
【００６５】
［反復再構成処理：ステップＳ４２］
図１２は、散乱源存在可能領域の決定処理の流れを示したフローチャートである。
【００６６】
まず、再構成処理部５０３は、暫定的な解画像を設定する（ステップＳ４２１）。本実施形態では、暫定解画像の初期値として、ステップＳ４１において得られた解を用いる。
【００６７】
次に、再構成処理部５０３は、暫定解画像の散乱源存在可能領域以外の領域を０に置き換えた後（ステップＳ４２２）、プロジェクション処理によりプロジェクション画像を生成し、測定値としての投影画像との差を求め（ステップＳ４２３）、この差を投影画像として用いて再度再構成する（ステップＳ４２４、ステップＳ４２５）。
【００６８】
次に、再構成処理部５０３は、差分の再構成画像を、1回目の暫定解画像（０に置き換えたもの）から引き算し（ステップＳ４２６）、当該再構成処理の反復回数が所定回数に到達したか否かを判定する。所定回数に到達していないと判定した場合には、再構成処理部５０３は、ステップＳ４２５で得られた最新の投影画像を暫定的な解画像として、ステップＳ４２２〜ステップＳ４２６間での処理を繰り返すことその暫定解画像を反復修正し最終結果に向けて収束させていく。一方、所定回数に到達したと判定した場合には、再構成処理部５０３は、反復再構成処理を終了する。
【００６９】
（効果）
以上述べた構成によれば、次の効果を実現することができる。
【００７０】
本放射線治療システムでは、散乱線検出器の各位置での測定結果を用いたバックプロジェクション処理により、被検体内の散乱線源存在領域を決定し、当該散乱源存在可能領域外の暫定解を０に置き換えて初回の再構成を実行する。これにより、誤差分の投影画像は、１回目の反復の再構成でその大部分が散乱源存在可能領域に再構成され、散乱源存在可能領域の暫定解画像の値は大きく減少する。また、次の再構成処理においては、前回の再構成処理で得られた暫定解画像を入力画像とし、その散乱源存在可能領域外の暫定解を０に置き換えて再び再構成処理を実行し、これを繰り返す。この様な反復再構成処理により、散乱源存在可能領域に対応する信号を忠実に映像化することができると共に、散乱源存在可能領域に対応しない信号を除外することができ、断層撮像特有のアーチファクトを低減させ、画質を向上させることができる。
【００７１】
なお、本条件付き反復再構成処理は、少なくとも２解の繰り返しでも再構成誤差を改善するために大きな効果がある。さらに反復回数を増やせばより効果があるが、最大でも1０回の反復すれば十分である。
【００７２】
（第２の実施形態）
第１の実施形態では、散乱源が存在すべきではない領域に対し散乱源の分布が０であるという条件を導入することで、解の精度の向上を図るものであった。これに対し、本実施形態では、この散乱源の分布が０という条件をより一般化した方法について説明する。散乱源が存在すべきではない領域に対し散乱源の分布が０という条件の導入は2つの効果がある。
【００７３】
・従来の再構成方法では、散乱源が本来０になるべき領域で０より大きい解が得られることがある。本来０であるはずの領域なので、この領域について散乱源＝０という条件を付加することで解の精度の向上を図ることができる。
【００７４】
・従来の再構成方法では、散乱源が小さな値を持つ領域で負の値の解が得られることがある。特に、散乱源が本来０になるべき領域では負の値の解が得られやすい。原理的に散乱源の値は負の値にはならないので、この領域について散乱源=０という条件を付加することで解の精度の向上図ることができる。
【００７５】
本実施形態では、上記２つの置き換えを下記のように一般化し、解の上限および下限についての条件に従う反復再構成処理を行う。
【００７６】
・解を投影した画像が、再構成処理において入力となる投影画像（再構成入力投影画像）を超えない。
【００７７】
・解は負の値にはならない（非負条件）。
【００７８】
上の条件を考慮した上で、「解を投影した画像が再構成入力投影画像にできるだけ近くなる」ように、解を選ぶという方針をとる。なお、第１の実施形態では、散乱源が本来０になるべき領域についてだけ条件を付加することができたが、本第２の実施形態では、上の２つの条件で解の上限と下限を定めているため、全領域について条件を付加することが可能になっている。
【００７９】
再構成処理部５０３において実行される具体的な処理の内容は、以下の様である。
【００８０】
［基本となる方法］
再構成空間の画像ｆのプロジェクション処理をＡｆと表記する。この処理の結果は投影空間のベクトルである。また、投影画像ｇのバックプロジェクション処理をＡＴｇと表記する。この処理の結果は、再構成空間のベクトルとなる。
【００８１】
まず、本手法のベースとなる従来方法について説明する。次の式（１）を満たすｆを求める場合を考える。
【数１】
【００８２】
Ｐｆは解ｆのペナルティーベクトルであり、望ましくないｆに対してｆＴＰｆが大きな値を持つ。常にｆＴＰｆ≧０である。上記式（１）より、次の式（２）が得られる。
【数２】
【００８３】
また、式（２）の右辺を導関数とする関数ｑ（ｆ）を次の式（３）によって定義する。
【数３】
【００８４】
すなわち、式（１）を解くには、式（３）を目的関数として最小化すればよい。式（３）式を最小化するには、共役勾配法などの最小化手法が利用できる。共役勾配法を含む多くの最小化手法では評価関数ｑの計算とその導関数の計算が必要であるが、それらの処理要素であるＡｆおよびＡＴｇは前述のとおりプロジェクション処理とバックプロジェクション処理なので、それらプロジェクション処理とバックプロジェクション処理を実装することが必要である。参考のため、ｑ（ｆ）の導関数を次の式（４）に示す。
【数４】
【００８５】
ある方向ｈに修正しようとしたとき（具体的には解ｆにλｈを加える）、方向ｈにどれだけの距離を進めば最小になるかを決定することが必要になる。これはｑ（ｆ＋λｈ）を最小にするλを決めることであり、次の式（５）を解くことで求めることができる。
【数５】
【００８６】
以下に、式（５）を解くための計算を示す。
【数６】
【００８７】
以上を具体的に実装した従来法の処理手順を図１３に示す。説明上の例として最急降下法の例を示してしるが、実際には、共役勾配法等、より優れた方法を用いるべきである。
【００８８】
［非負条件の導入］
本実施形態では、解は負の値にならないという条件を導入するため、図１３に示した処理手順を、図１４に示す処理手順に変更する（ステップ９ｂが新たに追加されたステップである）。
【００８９】
［解を投影した画像が、再構成入力投影画像を超えないという条件の導入］
式（３）に変えて、次の式（６）に示す目的関数を用いる。
【数７】
【００９０】
ここで、Ｐｏｓｉｔｉｖｅ（ｖ）は要素νｉの値が正の値のときｉ番目の要素の値がνｉ、０以下の値のときｉ番目の要素の値が０になる関数である。Ｎｅｇａｔｉｖｅ（ｖ）は逆に要素νｉの値が０以下のときｉ番目の要素の値がνｉ、正の値のときｉ番目の要素の値が０になる関数である。式（６）は、αｇ＝αｓ＝１のとき、式（３）と等価である。典型的にはαｇ＞αｓとなる値を用いる。
【００９１】
式（６）の導関数として次の式（７）を用いる。
【数８】
【００９２】
式（６）のＡｆ−ｇは暫定解の投影と入力投影画像の差である。αｇ＞αｓとなる値を用いたとき、入力投影画像を暫定解の投影が超えている場合、より大きな係数がかかるため、目的関数（６）はより大きな値になる。
【００９３】
図１５は、αｇ＝２、αｓ＝１の例である。最適化の過程では目的関数を小さくする解を選ぶので、結果として、暫定解の投影が入力投影画像を大きく超えない解が選ばれるようになる。
【００９４】
なお、処理手順は、図１６に示すようになる。この処理手順の特徴的な点は、ステップ6bでの導関数の算出において、暫定解ｈｉの投影と入力投影画像を比較し、Ａｆ−ｇの各要素の符号に応じ、正の値の方が大きい係数をＡｆ−ｇの各要素に乗じる処理が追加される点である。ステップ６ｂの変更に伴い、ステップ８で算出されるλが不正確なものになるので、ステップ８ｂを追加してλをより正確な値に修正することが望ましい。
【００９５】
（効果）
以上述べた本実施形態に係る条件付き反復再構成処理によれば、解が非負であるという下限に関する条件と、解を投影した画像が、再構成入力投影画像を超えないという上限に関する条件の導入が可能になり、解を制限することができる。例えば、特に散乱線画像が０に近い値を示す領域に対応する被検体の位置において、散乱源分布が０に近い値であるという条件が付加される。従って、第２の実施形態は、第１の実施形態と同様の効果を持っている。さらに、散乱線画像が０でないある一定値を示す領域においても、散乱源分布を投影した値がそれより小さい値になるという条件が付加される。このことは、第２の実施形態では第１の実施形態より強い制限事項を再構成に課している。従って、第２の実施形態によれば、第１の実施形態に比して、再構成の精度をさらに向上できるという効果がある。
【００９６】
（第３の実施形態）
第２の実施形態では、解全体について非負の条件を課すことができたが、このような方法では非負条件が有効に機能しないことがある。例えば、図１７（ａ）に示すように散乱源が一部だけ大きな値を持つような場合、図１７（ｂ）に示すようにその周囲では負の値を持つ誤った解が得られやすい。この場合、非負条件により解が改善することができる。しかし、図１７（ａ）に対して図１７（ｃ）に示すように一様な散乱源分布が加わる場合、図１７（ｂ）と同様の再構成方式では図１７（ｄ）に示すように解が負の値にならない。つまり、解に図１７（ｂ）と同じ間違いが含まれていても、一様成分が加わるために解が負の値にならず間違いを修正することができない。
【００９７】
そこで、本実施形態では、モデル化された複数の分布の合成で解（散乱源の分布）を表わす。すなわち、本来は一つである解を複数のモデル化された解に分解して考える。そのうえで、分解した解の各々について第２の実施形態の条件を課すことで、例えば図１７（ｃ）に示すような一様な散乱分布が存在する場合であっても、好適な画像再構成を実現するものである。
【００９８】
解の合成は、次の式（８）で表すことができる。
【数９】
【００９９】
上記式（８）では、解がＭ個の関数の合成で表されているが、そのうち、先頭のＲ個のみで近似した合成式が、次の式（９）で算出できるようにしておく。
【数１０】
【０１００】
分解された解分布ｆｊのうち小さなｊに対応する分布は広域的な分布とし、ｊが大きくなるほど局所分布となるようにｆｊを定める。とくに、ｆ０は一様分布である。ｆｊの各要素は全て０以上の値を持つように設定するのが良い。
【０１０１】
目的関数は、式（６）に式（８）を代入することで、次の式（１０）の様になる。
【数１１】
【０１０２】
また、導関数として次の式（１１）を用いる。
【数１２】
【０１０３】
ＦＲｍＲの値はＲ個の解分布ｆｊの合成なので、目的関数および導関数はＲの値により計算処理が異なる。Ｒを決めれば、第２の実施形態の方法でｍＲを求め、ＦＲｍＲを計算することで解分布を求めることができる。その際、第２の実施形態の記載中のＡは本実施形態のＡＦＲに置き換える必要がある。
【０１０４】
本実施形態では、Ｒ＝１から開始して、Ｒを徐々に大きくしながら、図１９に示すように反復的に再構成処理を実行する。投影画像ｇは最初再構成入力画像の値を設定しておくが、反復ごとに解ｆｋの投影Ａｆｋを引き算していく。反復中のステップｋにおいては、それより前の全ステップで差し引かれた投影画像からの再構成を行う。ｆｊは小さなｊに対する分布がより広い分布を持っており、反復が進むほどＲが増加していくので、ステップが進むほど細かい構造が再構成されるようになる。一例として最初は一様分布が推定され、その後、局所の分布が推定される。２ステップ目以降は、投影画像から一様分布に対する投影画像が除去されているので、最初に示した概念図で説明したように以降の各ステップで局所の分布を推定する際も非負条件が有効に働く。
【０１０５】
以上述べた本実施形態に係る再構成処理によれば、図１８（ｃ）の例（図１７（ｃ）の例と同じ）の場合、解を図１８（ｃ１）に示す一様成分と図１８（ｃ２）に示す局所成分に分解し、各々に非負の条件を課して再構成する。このとき、非負条件が加えられるので、解に含まれる誤差がそれぞれ図１８（ｅ１）、図１８（ｅ２）に示すように低減される。最終的に得られる解はそれぞれを合成であるが、もちろん、最終的な解に含まれる誤差も図１８（ｅ）に示すように低減されている。
【０１０６】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。具体的な変形例としては、例えば次のようなものがある。
【０１０７】
例えば、本実施形態に係る各機能は、当該処理を実行するプログラムをワークステーション等のコンピュータにインストールし、これらをメモリ上で展開することによっても実現することができる。このとき、コンピュータに当該手法を実行させることのできるプログラムは、磁気ディスク（フロッピー（登録商標）ディスク、ハードディスクなど）、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤなど）、半導体メモリなどの記録媒体に格納して頒布することも可能である。
【０１０８】
また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【産業上の利用可能性】
【０１０９】
以上本発明によれば、再構成処理において、散乱線が照射された領域についての条件を再構成に導入して、得られる散乱源分布の精度を向上することができる放射線治療システム及び放射線治療プログラムを実現することができる。
【図面の簡単な説明】
【０１１０】
【図１】図１は、本放射線治療用線量分布測定装置の治療用放射線に基づく被検体からの散乱線計測の原理、方法を説明するための図である。
【図２】図２は、本実施形態に係る放射線治療用線量分布測定装置１のブロック構成図を示している。
【図３】図３は、吸収線量画像データの生成処理を含む放射線治療時における処理の流れを示したフローチャートである。
【図４】図４は、本放射線治療システム１の散乱線の測定形態を示した図である。
【図５】図５は、本実施形態に係る吸収線量画像データの生成処理を含む放射線治療時における処理の流れを示したフローチャートである。
【図６】図６は、本放射線治療システム１の散乱線の測定形態の一例を示した図である。
【図７】図７は、本放射線治療システム１の散乱線の測定形態の他の例を示した図である。
【図８】図８（ａ）、（ｂ）は、条件付き反復再構成機能の概念を説明するための図である。
【図９】図９は、本条件付き反復再構成機能に従う処理の流れを示したフローチャートである。
【図１０】図１０は、散乱源存在可能領域の決定処理の流れを示したフローチャートである。
【図１１】図１１（ａ）、（ｂ）、（ｃ）は、散乱源存在可能領域の決定処理を説明するための図である。
【図１２】図１２は、散乱源存在可能領域の決定処理の流れを示したフローチャートである。
【図１３】図１３は、最急降下法の処理手順を示したフローチャートである。
【図１４】図１４は、非負条件が導入された最急降下法の処理手順を示したフローチャートである。
【図１５】図１５は、αｇ＝２、αｓ＝１の場合の入力投影画像及び暫定解の投影画像の例を示した図である。
【図１６】図１６は、解を投影した画像が、再構成入力投影画像を超えないという条件が導入された場合の最急降下法の処理手順を示したフローチャートである。
【図１７】図１７は、第３の実施形態に係る再構成処理の概念を説明するための図である。
【図１８】図１８は、第３の実施形態に係る再構成処理の概念を説明するための図である。
【図１９】図１９は、第３の実施形態に係る再構成処理の流れを示したフローチャートである。
【符号の説明】
【０１１１】
１…放射線治療用線量分布測定装置、２…放射線照射システム、３…散乱線検出システム、４…データ取得制御部、５…データ処理システム、６…表示部、７…記憶部、８…操作部、９…ネットワークＩ／Ｆ、２０１…電力供給部、２０３…照射部、２０５…タイミング制御部、２０７…ガントリ制御部、３０１…検出器、３０３…コリメータ、３０５…移動機構部、３０７…位置検出部、５０１…補正処理部、５０２…散乱源存在領域決定部、５０３…再構成処理部、５０５…変換処理部、５０７…画像処理部
【特許請求の範囲】
【請求項１】
被検体に対して治療用放射線ビームを照射する照射手段と、
前記治療用放射線ビームに基づいて発生する前記被検体内からの散乱線を検出し散乱線データを発生する検出手段と、
前記検出された各散乱線データ上における、前記被検体内の散乱源存在領域を決定する決定手段と、
前記散乱源存在領域が決定された前記各散乱線データを入力とし所定の条件に従う再構成処理を繰り返し実行することで再構成画像を生成し、当該再構成画像を用いて、前記被検体内における散乱線発生密度の三次元的分布を示す散乱線ボリュームデータを生成する再構成手段と、
前記散乱線ボリュームデータを、吸収された放射線量の三次元分布を示す吸収線量ボリュームデータに変換する変換手段と、
前記吸収線量ボリュームデータを用いて、前記被検体内における吸収線量画像を生成する画像生成手段と、
前記吸収線量画像を表示する表示手段と、
を具備することを特徴とする放射線治療システム。
【請求項２】
前記再構成手段は、前記各再構成処理において生成される前記再構成画像が、前記各再構成処理において入力とされる前記各散乱線データを超えないとする条件に従って、前記再構成処理を繰り返し実行することを特徴とする請求項１記載の放射線治療システム。
【請求項３】
前記再構成手段は、入力とされる前記各散乱線データについて、散乱線が負の値である領域についてはその値を０とする条件に従って、前記再構成処理を繰り返し実行することを特徴とする請求項１又は２記載の放射線治療システム。
【請求項４】
前記再構成手段は、入力とされる前記各散乱線データについて、前記散乱源存在領域に対応しない領域については散乱線を０とする条件に従って、前記再構成処理を繰り返し実行することを特徴とする請求項１又は２記載の放射線治療システム。
【請求項５】
前記再構成手段は、入力とされる前記各散乱線データを複数のモデル解に分解し、各モデル解について前記再構成処理を実行することを特徴とする請求項３又は４記載の放射線治療システム。
【請求項６】
コンピュータに、
被検体に対して照射された治療用放射線ビームに基づいて発生する前記被検体内からの散乱線を検出させる検出機能と、
前記検出された散乱線に基づいて複数の散乱線データを生成させるデータ生成機能と、
前記検出された各散乱線データ上における、前記被検体内の散乱源存在領域を決定させ決定機能と、
前記散乱源存在領域が決定された前記各散乱線データを入力とし所定の条件に従う再構成処理を繰り返し実行することで再構成画像を生成させ、当該再構成画像を用いて、前記被検体内における散乱線発生密度の三次元的分布を示す散乱線ボリュームデータを生成させる再構成機能と、
前記散乱線ボリュームデータを、吸収された放射線量の三次元分布を示す吸収線量ボリュームデータに変換させる変換機能と、
前記吸収線量ボリュームデータを用いて、前記被検体内における吸収線量画像を生成させる画像生成機能と、
前記吸収線量画像を表示させる表示機能と、
を実現させることを特徴とする放射線治療プログラム。
【請求項７】
前記再構成機能においては、前記各再構成処理において生成される前記再構成画像が、前記各再構成処理において入力とされる前記各散乱線データを超えないとする条件に従って、前記再構成処理を繰り返し実行させることを特徴とする請求項６記載の放射線治療プログラム。
【請求項８】
前記再構成機能においては、入力とされる前記各散乱線データについて、散乱線が負の値である領域についてはその値を０とする条件に従って、前記再構成処理を繰り返し実行させることを特徴とする請求項６又は７記載の放射線治療プログラム。
【請求項９】
前記再構成機能においては、入力とされる前記各散乱線データについて、前記散乱源存在領域に対応しない領域については散乱線を０とする条件に従って、前記再構成処理を繰り返し実行させることを特徴とする請求項６又は７記載の放射線治療プログラム。
【請求項１０】
前記再構成機能においては、入力とされる前記各散乱線データを複数のモデル解に分解させ、各モデル解について前記再構成処理を実行させることを特徴とする請求項８又は９記載の放射線治療プログラム。
【請求項１】
被検体に対して治療用放射線ビームを照射する照射手段と、
前記治療用放射線ビームに基づいて発生する前記被検体内からの散乱線を検出し散乱線データを発生する検出手段と、
前記検出された各散乱線データ上における、前記被検体内の散乱源存在領域を決定する決定手段と、
前記散乱源存在領域が決定された前記各散乱線データを入力とし所定の条件に従う再構成処理を繰り返し実行することで再構成画像を生成し、当該再構成画像を用いて、前記被検体内における散乱線発生密度の三次元的分布を示す散乱線ボリュームデータを生成する再構成手段と、
前記散乱線ボリュームデータを、吸収された放射線量の三次元分布を示す吸収線量ボリュームデータに変換する変換手段と、
前記吸収線量ボリュームデータを用いて、前記被検体内における吸収線量画像を生成する画像生成手段と、
前記吸収線量画像を表示する表示手段と、
を具備することを特徴とする放射線治療システム。
【請求項２】
前記再構成手段は、前記各再構成処理において生成される前記再構成画像が、前記各再構成処理において入力とされる前記各散乱線データを超えないとする条件に従って、前記再構成処理を繰り返し実行することを特徴とする請求項１記載の放射線治療システム。
【請求項３】
前記再構成手段は、入力とされる前記各散乱線データについて、散乱線が負の値である領域についてはその値を０とする条件に従って、前記再構成処理を繰り返し実行することを特徴とする請求項１又は２記載の放射線治療システム。
【請求項４】
前記再構成手段は、入力とされる前記各散乱線データについて、前記散乱源存在領域に対応しない領域については散乱線を０とする条件に従って、前記再構成処理を繰り返し実行することを特徴とする請求項１又は２記載の放射線治療システム。
【請求項５】
前記再構成手段は、入力とされる前記各散乱線データを複数のモデル解に分解し、各モデル解について前記再構成処理を実行することを特徴とする請求項３又は４記載の放射線治療システム。
【請求項６】
コンピュータに、
被検体に対して照射された治療用放射線ビームに基づいて発生する前記被検体内からの散乱線を検出させる検出機能と、
前記検出された散乱線に基づいて複数の散乱線データを生成させるデータ生成機能と、
前記検出された各散乱線データ上における、前記被検体内の散乱源存在領域を決定させ決定機能と、
前記散乱源存在領域が決定された前記各散乱線データを入力とし所定の条件に従う再構成処理を繰り返し実行することで再構成画像を生成させ、当該再構成画像を用いて、前記被検体内における散乱線発生密度の三次元的分布を示す散乱線ボリュームデータを生成させる再構成機能と、
前記散乱線ボリュームデータを、吸収された放射線量の三次元分布を示す吸収線量ボリュームデータに変換させる変換機能と、
前記吸収線量ボリュームデータを用いて、前記被検体内における吸収線量画像を生成させる画像生成機能と、
前記吸収線量画像を表示させる表示機能と、
を実現させることを特徴とする放射線治療プログラム。
【請求項７】
前記再構成機能においては、前記各再構成処理において生成される前記再構成画像が、前記各再構成処理において入力とされる前記各散乱線データを超えないとする条件に従って、前記再構成処理を繰り返し実行させることを特徴とする請求項６記載の放射線治療プログラム。
【請求項８】
前記再構成機能においては、入力とされる前記各散乱線データについて、散乱線が負の値である領域についてはその値を０とする条件に従って、前記再構成処理を繰り返し実行させることを特徴とする請求項６又は７記載の放射線治療プログラム。
【請求項９】
前記再構成機能においては、入力とされる前記各散乱線データについて、前記散乱源存在領域に対応しない領域については散乱線を０とする条件に従って、前記再構成処理を繰り返し実行させることを特徴とする請求項６又は７記載の放射線治療プログラム。
【請求項１０】
前記再構成機能においては、入力とされる前記各散乱線データを複数のモデル解に分解させ、各モデル解について前記再構成処理を実行させることを特徴とする請求項８又は９記載の放射線治療プログラム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５】
【図１６】
【図１７】
【図１８】
【図１９】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５】
【図１６】
【図１７】
【図１８】
【図１９】
【公開番号】特開２００９−１８９５５９（Ｐ２００９−１８９５５９Ａ）
【公開日】平成２１年８月２７日（２００９．８．２７）
【国際特許分類】
【出願番号】特願２００８−３３３６２（Ｐ２００８−３３３６２）
【出願日】平成２０年２月１４日（２００８．２．１４）
【出願人】（０００００３０７８）株式会社東芝 (54,554)
【出願人】（５９４１６４５４２）東芝メディカルシステムズ株式会社 (4,066)
【Ｆターム（参考）】
【公開日】平成２１年８月２７日（２００９．８．２７）
【国際特許分類】
【出願日】平成２０年２月１４日（２００８．２．１４）
【出願人】（０００００３０７８）株式会社東芝 (54,554)
【出願人】（５９４１６４５４２）東芝メディカルシステムズ株式会社 (4,066)
【Ｆターム（参考）】
[ Back to top ]