(57)【特許請求の範囲】
【請求項１】
治療計画時に取得した３次元ＣＴ画像データと位置決め時に取得する２次元Ｘ線透視像のレジストレーションを行う際に、前記３次元ＣＴ画像データと前記２次元Ｘ線透視像との位置合わせを行うために、コンピュータ上にＸ線撮影のジオメトリを再現し、ＣＴ画像データに対して仮想的に透視投影を行うことで、２次元ＤＲＲを得る手段と、
該ＤＲＲと前記Ｘ線透視像の一致度が最大となるように、ＣＴ画像データの回転・平行移動に関する透視投影パラメータを、最短経路を探索して最適化し、患者の現在位置と治療計画時位置とのずれ量を算出する手段とを備えた放射線治療における患者自動位置決め装置において、
前記ＤＲＲと前記Ｘ線透視像とのずれ量を、相互情報量を評価関数として用いる相互情報量アルゴリズムで計算する手段と、
前記ＤＲＲと前記Ｘ線透視像とのずれ量を、微分値でエッジを強調させるＧｒａｄｉｅｎｔ  Ｄｉｆｆｅｒｅｎｃｅを評価関数として用いるエッジ検出アルゴリズムで計算する手段と、
前記相互情報量アルゴリズムにより計算されたずれ量と、前記エッジ検出アルゴリズムにより計算されたずれ量とが異なる場合、それぞれを重み付けて統合する手段と、
を備えたことを特徴とする放射線治療における患者自動位置決め装置。
【請求項２】
前記相互情報量が規格化相互情報量であることを特徴とする請求項１に記載の放射線治療における患者自動位置決め装置。
【請求項３】
治療計画時に取得した３次元ＣＴ画像データと位置決め時に取得する２次元Ｘ線透視像のレジストレーションを位置決め計算機で行う際に、前記３次元ＣＴ画像データと前記２次元Ｘ線透視像との位置合わせを行うために、コンピュータ上にＸ線撮影のジオメトリを再現させ、ＣＴ画像データに対して仮想的に透視投影させることで、コンピュータで２次元ＤＲＲを得て、
該ＤＲＲと前記Ｘ線透視像の一致度が最大となるように、ＣＴ画像データの回転・平行移動に関する透視投影パラメータを、コンピュータが最短経路を探索して最適化し、患者の現在位置と治療計画時位置とのずれ量をコンピュータが算出するようにした放射線治療における患者自動位置決め方法において、
コンピュータが、前記ＤＲＲと前記Ｘ線透視像とのずれ量を、相互情報量を評価関数として用いる相互情報量アルゴリズムで計算すると共に、
コンピュータが、前記ＤＲＲと前記Ｘ線透視像とのずれ量を、微分値でエッジを強調させるＧｒａｄｉｅｎｔ  Ｄｉｆｆｅｒｅｎｃｅを評価関数として用いるエッジ検出アルゴリズムで計算し、
コンピュータが、前記相互情報量アルゴリズムにより計算されたずれ量と、前記エッジ検出アルゴリズムにより計算されたずれ量とが異なる場合、それぞれを重み付けて統合することを特徴とする放射線治療における患者自動位置決め方法。
【請求項４】
治療計画時に取得した３次元ＣＴ画像データと位置決め時に取得する２次元Ｘ線透視像のレジストレーションを行う際に、前記３次元ＣＴ画像データと前記２次元Ｘ線透視像との位置合わせを行うために、コンピュータ上にＸ線撮影のジオメトリを再現し、ＣＴ画像データに対して仮想的に透視投影を行うことで、２次元ＤＲＲを得て、
該ＤＲＲと前記Ｘ線透視像の一致度が最大となるように、ＣＴ画像データの回転・平行移動に関する透視投影パラメータを、最短経路を探索して最適化し、患者の現在位置と治療計画時位置とのずれ量を算出するようにした放射線治療における患者自動位置決め方法をコンピュータに行なわせるための患者自動位置決め用プログラムであって、
前記ＤＲＲと前記Ｘ線透視像とのずれ量を、相互情報量を評価関数として用いる相互情報量アルゴリズムでコンピュータに計算させるステップと、
前記ＤＲＲと前記Ｘ線透視像とのずれ量を、微分値でエッジを強調させるＧｒａｄｉｅｎｔ  Ｄｉｆｆｅｒｅｎｃｅを評価関数として用いるエッジ検出アルゴリズムでコンピュータに計算させるステップと、
前記相互情報量アルゴリズムにより計算されたずれ量と、前記エッジ検出アルゴリズムにより計算されたずれ量とが異なる場合、それぞれを重み付けてコンピュータに統合させるステップと、
を備えたことを特徴とする放射線治療における患者自動位置決め用プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、放射線治療における患者自動位置決め装置及び方法並びに患者自動位置決め用プログラムに係り、特に、複数の評価関数を用いることで、骨が多い領域や臓器が多い領域など部位毎の特性に合った位置決めを行い、患者の現在位置と治療計画時位置とのずれ量誤差を軽減することが可能な放射線治療における患者自動位置決め装置及び方法並びに患者自動位置決め用プログラムに関する。
【背景技術】
【０００２】
放射線治療においては、通常、治療計画と照射が別に行われる。従って、治療計画時と照射時の患者体位の再現性を保つよう、治療の際に患者の位置決めを正確に行うことで治療精度が向上する。それ故、患者位置決めは放射線治療の過程では欠かせなくなっている。
【０００３】
従来、レーザを使用して患者位置決めが行われていたが、最近では高精度化のためにＸ線画像を使用して患者位置決めを行う技術が開発されている（特許文献１、２）。
【０００４】
放射線治療は、通常、分割して行われ、現在は、図１に示す如く、治療の度に技師によって手動による正確な位置決めが行われている。これを図２に示す如く自動化するため、治療計画時に取得した３次元ＣＴ画像データと位置決め時に取得する２次元Ｘ線透視像のレジストレーションが行われている（非特許文献１）。
【０００５】
具体的には、３次元ＣＴ画像データと、２次元Ｘ線透視像の位置合わせを行うために、コンピュータ上にＸ線撮影のジオメトリを再現し、ＣＴ画像データに対して仮想的に透視投影を行うことで、２次元ＤＲＲを得る。ＤＲＲとＸ線透視像の一致度が最大となるように、ＣＴ画像データの回転・平行移動に関する透視投影パラメータを最適化し、患者の現在位置と治療計画時位置とのずれ量を算出する。
【０００６】
更に、２方向Ｘ線画像と治療計画時ＣＴ画像データから、患者ずれ量を自動計算する方法も提案されている。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００７−２８２８７７号公報
【特許文献２】特開２０１０−５７８１０号公報
【非特許文献】
【０００８】
【非特許文献１】G. P. Penney，J. Weese，et al.：A comparison of similarity measures for use in 2D-3D medical image registration，IEEE Trans. Inform Technol.Biomed．17，4，586−595，1998.
【発明の概要】
【発明が解決しようとする課題】
【０００９】
非特許文献１に記載の方法では、パラメータの最適化のために、反復してＤＲＲを作成する必要があるが、ＤＲＲ作成は多大な計算コストを要するため、短時間で位置決めを行うことが困難であるという問題点を有していた。即ち、位置決めの高精度化と、それにかかる時間はトレードオフとなり、実際の位置決め時間に１５〜３０分かかっており、その間、同じ体位を維持する患者にとっては苦痛である。
【００１０】
又、自動位置決めの従来技術は、１つの位置決めアルゴリズムにより位置ずれ量を算出していた。体の部位により、また患者による個人差により、最適な位置決めアルゴリズムは異なっている。これをユーザが判別することは難しい。従って、位置決め部位によっては、高い精度で計算結果を提示できなかった。
【００１１】
又、従来は、自動位置決め計算により算出された結果を、治療台移動へ反映するときに考慮すべき種々の不確定要素まで考慮されていない。
【００１２】
これにより、手動位置決めではユーザの指示誤差や、目視判断による精度低下が生じていた。
【００１３】
又、従来の患者位置決め画像は静止画像を入力としていたため、特に呼吸性移動が影響する体幹部では、参照画像と位置決め画像の呼吸位相が異なり、本質的に患者形状が異なっていたことも問題であった。
【００１４】
本発明は、前記従来の問題点を解決するべくなされたもので、複数の評価関数を用いることで、骨が多い領域や臓器が多い領域など部位毎の特性に合った位置決めを行い、患者の現在位置と治療計画時位置とのずれ量誤差を軽減することを課題とする。
【課題を解決するための手段】
【００１５】
現状の患者位置決めは、Ｘ線透視像を直交する２方向で取得するため、ＤＲＲも同様に２方向で作成し、各方向の評価値を足し合わせることで全体の評価値としている。しかし、投影軸に沿った方向に対する平行移動量パラメータを変化させた場合は、投影面においてＤＲＲが拡大・縮小という変化となり、画像の変化が小さい。これは評価値においても同様な変化が生じると考えられ、投影軸に対する平行移動量パラメータの最適化の際には、これに直交する１方向でのみＤＲＲと実透視像の一致度を評価値の算出を行うことで、十分最適化可能ではないかと予想される。また、位置決め対象物体の、投影軸まわりの回転移動量パラメータを変化させた場合は、投影面においてもＤＲＲが回転移動するため、画像に生じる変化が大きく、１方向でのみＤＲＲと実透視像の一致度を評価値の算出を行うことで、十分最適化可能ではないかと予想される。
【００１６】
又、ＣＴ画像とＸ線画像という異なるモダリティー間の画像レジストレーションは、相互情報量を用いたレジストレーションが用いられるが、最適化計算時にlocal minimaに落ちる場合がある。また、高吸収物である骨形状に注目し、比較的計算量の少ないエッジ検出アルゴリズムを用いてレジストレーションを行う場合もある。しかし、異なるモダリティーにより取得された画像間のレジストレーションにおいては、誤差量が多くなるという欠点がある。そこで、これら２つのアルゴリズムを用いて計算することで、計算時間と精度のバランスが取れた結果を出すことができる。このように、複数の評価関数で位置ずれ量を算出した場合、これらは独立に計算される。そのため、それぞれの評価関数による位置ずれ量が異なる値を示す場合もある。それぞれの評価関数ごとの結果に、重み付け処理を行い最終的な、位置ずれ量補正値を算出することができる。
【００１７】
又、最近では、４次元ＣＴや透視撮影により時系列のデータを取得できるようになった。そこで、本発明法を時間軸に展開することで、これら時系列のデータを入力させることで、時間位相もあわせることが可能となる。
【００１８】
即ち、４次元ＣＴ撮影による時系列の参照ＣＴ画像、並びに透視装置による時系列の撮影画像を扱う場合が想定される。参照ＣＴ画像の位相を指定するか、撮影画像の位相を指定するかの２つが考えられるが、どちらの方法でも、静止画像の場合と同じプロセスに、位相パラメータを加えて、順次位相を変えながら、評価関数による計算を行えば良い。
【００１９】
本発明は、このような知見に基づいてなされたもので、治療計画時に取得した３次元ＣＴ画像データと位置決め時に取得する２次元Ｘ線透視像のレジストレーションを行う際に、前記３次元ＣＴ画像データと前記２次元Ｘ線透視像との位置合わせを行うために、コンピュータ上にＸ線撮影のジオメトリを再現し、ＣＴ画像データに対して仮想的に透視投影を行うことで、２次元ＤＲＲを得る手段と、該ＤＲＲと前記Ｘ線透視像の一致度が最大となるように、ＣＴ画像データの回転・平行移動に関する透視投影パラメータを、最短経路を探索して最適化し、患者の現在位置と治療計画時位置とのずれ量を算出する手段とを備えた放射線治療における患者自動位置決め装置において、前記ＤＲＲと前記Ｘ線透視像とのずれ量を、相互情報量を評価関数として用いる相互情報量アルゴリズムで計算する手段と、前記ＤＲＲと前記Ｘ線透視像とのずれ量を、微分値でエッジを強調させるＧｒａｄｉｅｎｔ  Ｄｉｆｆｅｒｅｎｃｅを評価関数として用いるエッジ検出アルゴリズムで計算する手段と、前記相互情報量アルゴリズムにより計算されたずれ量と、前記エッジ検出アルゴリズムにより計算されたずれ量とが異なる場合、それぞれを重み付けて統合する手段と、を備えることにより、前記課題を解決したものである。
【００２１】
ここで、前記相互情報量を規格化相互情報量とすることができる。
【００２２】
本発明は、又、治療計画時に取得した３次元ＣＴ画像データと位置決め時に取得する２次元Ｘ線透視像のレジストレーションを位置決め計算機で行う際に、前記３次元ＣＴ画像データと前記２次元Ｘ線透視像との位置合わせを行うために、コンピュータ上にＸ線撮影のジオメトリを再現させ、ＣＴ画像データに対して仮想的に透視投影させることで、コンピュータで２次元ＤＲＲを得て、該ＤＲＲと前記Ｘ線透視像の一致度が最大となるように、ＣＴ画像データの回転・平行移動に関する透視投影パラメータを、コンピュータが最短経路を探索して最適化し、患者の現在位置と治療計画時位置とのずれ量をコンピュータが算出するようにした放射線治療における患者自動位置決め方法において、コンピュータが、前記ＤＲＲと前記Ｘ線透視像とのずれ量を、相互情報量を評価関数として用いる相互情報量アルゴリズムで計算すると共に、コンピュータが、前記ＤＲＲと前記Ｘ線透視像とのずれ量を、微分値でエッジを強調させるＧｒａｄｉｅｎｔ  Ｄｉｆｆｅｒｅｎｃｅを評価関数として用いるエッジ検出アルゴリズムで計算し、コンピュータが、前記相互情報量アルゴリズムにより計算されたずれ量と、前記エッジ検出アルゴリズムにより計算されたずれ量とが異なる場合、それぞれを重み付けて統合することを特徴とする放射線治療における患者自動位置決め方法を提供するものである。
【００２３】
本発明は、又、治療計画時に取得した３次元ＣＴ画像データと位置決め時に取得する２次元Ｘ線透視像のレジストレーションを行う際に、前記３次元ＣＴ画像データと前記２次元Ｘ線透視像との位置合わせを行うために、コンピュータ上にＸ線撮影のジオメトリを再現し、ＣＴ画像データに対して仮想的に透視投影を行うことで、２次元ＤＲＲを得て、該ＤＲＲと前記Ｘ線透視像の一致度が最大となるように、ＣＴ画像データの回転・平行移動に関する透視投影パラメータを、最短経路を探索して最適化し、患者の現在位置と治療計画時位置とのずれ量を算出するようにした放射線治療における患者自動位置決め方法をコンピュータに行なわせるための患者自動位置決め用プログラムであって、前記ＤＲＲと前記Ｘ線透視像とのずれ量を、相互情報量を評価関数として用いる相互情報量アルゴリズムでコンピュータに計算させるステップと、前記ＤＲＲと前記Ｘ線透視像とのずれ量を、微分値でエッジを強調させるＧｒａｄｉｅｎｔ  Ｄｉｆｆｅｒｅｎｃｅを評価関数として用いるエッジ検出アルゴリズムでコンピュータに計算させるステップと、前記相互情報量アルゴリズムにより計算されたずれ量と、前記エッジ検出アルゴリズムにより計算されたずれ量とが異なる場合、それぞれを重み付けてコンピュータに統合させるステップと、を備えたことを特徴とする放射線治療における患者自動位置決め用プログラムを提供するものである。
【発明の効果】
【００２４】
本発明によれば、ＤＲＲとＸ線透視像との一致度を相互情報量アルゴリズムを用いたレジストレーションとエッジ検出アルゴリズムを用いたレジストレーションとで評価することによって、骨が多い領域や臓器が多い領域など部位毎の特性に合った位置決めができる。従って、患者の現在位置と治療計画時位置のずれ量誤差を軽減するだけでなく、位置決め時間を短縮することにより、患者の苦痛を軽減させ、患者数のスループットも向上させることができる。
【００２５】
更に、前記位置決めに際して、独立に計算した、前記相互情報量アルゴリズムを用いて計算したずれ量と、前記エッジ検出アルゴリズムを用いて計算したずれ量とが異なる場合、それぞれを重み付けて統合し、最終的な位置ずれ量を算出するようにしたので、従来の自動位置決め法よりも高精度化し、患者の患部に依存しない位置決めを行うことができる。
【図面の簡単な説明】
【００２６】
【図１】従来の放射線治療の一例の手順を示すフローチャート
【図２】従来の放射線治療の他の例の手順を示すフローチャート
【図３】本発明が適用される放射線治療システムの実施形態の全体構成を示すブロック図
【図４】同じく患者周辺を示す断面図
【図５】同じく照射ポート周辺を示す側面図
【図６】同じく正面図
【図７】変形例のレンジシフタ周辺を示す側面図
【図８】前記実施形態で用いられているスキップ機能の概略を示す図
【図９】同じくコプラナ照射で１８０°治療台位置が異なる場合のスキップ機能を示す平面図
【図１０】本発明の第１実施形態における自動位置決めの処理手順を示すフローチャート
【図１１】同じく参照画像を取得している状態を示す正面図
【図１２】同じくＤＲＲ計算方法を示す斜視図
【図１３】同じく計算領域を示す図
【図１４】本発明の第２実施形態における自動位置決めの処理手順を示すフローチャート
【図１５】本発明の実施形態においてリハーサルを実施しない場合の参照画像の作成手順を示すフローチャート
【図１６】同じくリハーサルを実施する場合の参照画像の作成手順を示すフローチャート
【図１７】本発明の実施形態においてダウンサンプリング法の概略を示す図
【図１８】同じく自動位置決め計算時の座標計算方法を示す図
【図１９】同じく自動位置決め時のリハーサルずれ量補正方法を示す図
【図２０】実施例における自動位置決め前の、参照ＤＲＲ画像とＸ線画像の重畳表示を示す図
【図２１】同じく正規化相互情報量を評価関数とした自動位置決め後、参照ＤＲＲ画像とＸ線画像の（Ａ）垂直方向、（Ｂ）水平方向を示す重畳表示図
【図２２】同じくGradient Differenceを評価関数とした自動位置決め後、参照ＤＲＲ画像とＸ線画像の（Ａ）垂直方向、（Ｂ）水平方向を示す重畳表示図
【図２３】同じく規格化相互情報量とGradient Differenceを評価関数とした自動位置決め後、参照ＤＲＲ画像とＸ線画像の（Ａ）垂直方向、（Ｂ）水平方向を示す重畳表示図
【発明を実施するための形態】
【００２７】
以下図面を参照して、本発明の実施形態を詳細に説明する。
【００２８】
本発明が適用される放射線治療システムの実施形態は、図３に示す如く、患者１０が乗せられる患者ベッド１２の、例えば側方及び上方にそれぞれ配設された、水平及び垂直治療ビームをそれぞれ照射するための水平及び垂直照射ポート２０、２２と、Ｘ線画像を取得するための、例えば患者ベッド１２の側方及び下方にそれぞれ配置される水平方向Ｘ線管３０及び垂直方向Ｘ線管３２と、該Ｘ線管３０、３２の対向位置にそれぞれ配設される２次元検出器３４、３６と、予めＣＴ画像を撮影するためのＣＴ撮像装置５０と、該ＣＴ撮像装置５０の撮像結果に基づいて治療計画を行うための治療計画装置５２と、患者データベース（ＤＢ）５４と、患者ファイルサーバ５６と、前記２次元検出器３４、３６の出力に基づいて位置決め計算を行う位置決め計算機６０と、位置決め結果に基づいて機器を制御するための機器制御用計算機６２とを備えている。
【００２９】
ここで、前記垂直方向Ｘ線管３２は、図４に示す如く、治療室４０の床４２に埋めることで、床上に突起物がなくなり、又、可働部分がなくなるようにすることができる。一方、前記水平方向Ｘ線管３０は図の上下方向に移動可能とし、必要なときにのみ患者１０の視野に入ってくるようにして、患者への不安を低減することができる。
【００３０】
前記２次元検出器３４、３６は、例えばフラットパネル検出器（ＦＰＤ）やイメージインテンシファイヤ（Ｉ．Ｉ．）とし、いずれも照射ポート２０、２２のカバー内に設置して、患者との距離を十分保つとともに、接触事故を防ぐようにすることができる。前記２次元検出器３４、３６は、図５（２次元検出器３４の例）に詳細に示す如く、照射ポート内に設置されたレンジシフタ２４の患者側に配設され、患者位置決め時は図６（ａ）に示す如く、照射口２０Ａの面を覆い、一方、治療時は図６（ｂ）に示す如く、照射口２０Ａを開放して治療ビームの照射を妨げないようにされている。なお、この例では２次元検出器３４がレンジシフタ２４と独立して設けられているが、図７に示す変形例のように、レンジシフタ２４の一番患者側に２次元検出器３４を配設し、レンジシフタ２４の駆動機構２６を利用して出し入れするように構成することで、構成を簡略化することもできる。なお、２次元画像を取得できる物であれば、ＦＰＤやＩ．Ｉ．以外の２次元検出器を用いることも可能である。
【００３１】
以下、本実施形態における処理を詳細に説明する。
【００３２】
本実施形態においては、セットアップ番号が同じ場合など、すでに取得した位置決め結果データ（画像、ＵＩＤ、治療台位置、サムネイル画像）を、現在行っている位置決め結果へ複製する位置決めスキップを行う。
【００３３】
ノンコプラナ照射時の例を図８に示す。図８（ａ）では、Ｂｅａｍ１の位置決め画像と照射確認画像の取得が終了した状況を示している。次に（図８（ｂ））、Ｂｅａｍ２の位置決め画像取得時に、Ｂｅａｍ１の位置決めスキップを実施すると、位置決め画像データ（バイナリファイル、情報ファイル）、ＤＩＣＯＭ（規格）データ、そして治療台位置情報（治療台位置１）が、Ｂｅａｍ２の位置決め結果として作業フォルダにコピーされる。ＴＭＳ（治療計画装置）通知コマンドは、作業フォルダに置いてあるＤＩＣＯＭファイルのＵＩＤ（識別子）を取得送信するためＤＩＣＯＭデータもコピーしている。次に、照射確認画像取得時は、実際に撮影する必要がある。
【００３４】
右側照射／左側照射のコプラナ照射では、それぞれの照射方向で位置決めを行っているが、２回目の照射方向の位置決めスキップへの対応も可能である。たとえば、Ｂｅａｍ１における「位置決め画像」と「照射確認画像」取得が終了し（図９（ａ））、Ｂｅａｍ２を開始する状況とする。Ｂｅａｍ２は、コプラナ照射であるがＢｅａｍ１の位置決め時天板座標値から１８０°回転したベッド天板（以下、単に天板とも称する）位置である。１８０°天板を回転させてから「位置決め画像」取得、「照射確認画像」取得を行う場合は問題ない（図９（ｂ）（ｃ））。しかし、「位置決め画像」取得に対してスキップ機能を実行した場合（図９（ｄ））、Ｂｅａｍ１の位置決め結果情報がＢｅａｍ２の位置決め結果として登録される。次に、天板を１８０°回転させて（図９（ｅ））、照射確認画像を取得する。
【００３５】
次に、自動位置決め機能について説明する。
【００３６】
２Ｄ−３Ｄレジストレーションは、位置決め参照画像と撮影画像のずれ量算出を自動計算で行うモードである。位置決め参照画像は、治療計画ＣＴ画像となる。
【００３７】
自動レジストレーションを選択すると、ＣＴ画像データを読み込み、２方向参照画像を表示する。ＣＴ画像データは、１つのボリュームデータとなる。
【００３８】
スクリーン構成は、例えば、左上：垂直参照画像、右上：垂直撮影画像、左下：水平参照画像、右下：水平撮影画像の構成とすることができる。
【００３９】
２方向ＦＰＤ画像を撮影するが、すでに撮影されていれば、その画像が表示される（モード切替時に自動読込）。
【００４０】
レジストレーションプリセット（詳細後述）による簡便な設定値指定が可能である。
【００４１】
以下、図１０を参照して、自動位置決め手順を説明する。
【００４２】
１．参照ＣＴ画像取得。これは、あらかじめ撮影しておいた基準となるＣＴ画像データである。
【００４３】
２．入室した時の患者を２方向からＸ線撮影を行い、Ｘ線２次元画像（以下、「撮影画像」と呼ぶ。それぞれの撮影画像を「撮影画像Ａ」、「撮影画像Ｂ」と呼ぶ。）を取得する（図１１）。この撮影画像は、イメージインテンシファイヤ(I.I.)やフラットパネル検出器(ＦＰＤ)により取得される。
【００４４】
３．参照ＣＴ画像を用いて、撮影画像１と撮影画像２を取得したＸ線照射角度と同じ方向から、２次元画像(Ｄigital Ｒeconstructed Ｒadiography、以下「ＤＲＲ」と呼ぶ。
図１２）を計算し、それぞれ「参照画像Ａ」、「参照画像Ｂ」と呼ぶ。ここでは、Ｘ線管焦点から２次元検出器(Image detector)へ、Ｘ線パスに沿ってＣＴボクセル値を積算していく。
【００４５】
４．自動位置決め計算に使用する画像領域を、参照画像Ａ、Ｂ上に指定する（図１３）。ここでは、腹部領域の画像を使用し、自動位置決め計算を実施した。
【００４６】
５．Ａ、Ｂそれぞれの参照画像と撮影画像に対して、評価関数の計算結果を算出する。
【００４７】
６．ＣＴ画像を画像空間上で移動や回転を行い、参照画像を作成し、評価関数の計算結果を算出する。ＤＲＲを作成することで、擬似的に患者を移動回転することと、同じ状況を作っている。
【００４８】
７．参照画像を計算する位置角度を変えつつ、評価関数の計算を行う。この計算結果が最小になるまで上記６を繰り返す。参照画像の計算方向は、それぞれの位置パラメータ(移動３軸、回転３軸)を総当りで計算すると、莫大な繰り返し数と時間がかかるため、評価関数計算結果が最小値になるような最短経路を探索しつつ（最適化処理と言われている）、位置パラメータを変えていく。評価関数計算結果が最小となるときの位置パラメータが、撮影画像と参照画像の位置ずれ量の計算値となる。
【００４９】
８．この位置ずれ量補正値分、治療天板を移動させると、参照画像を撮影した位置と同じ位置に患者をセットすることができる。
【００５０】
Ａ．患者位置決め精度向上１
身体部位毎に、骨が多い領域や臓器が多い領域があるため、１つの評価関数ですべての部位の自動位置決め計算を行うと、部位により計算誤差が大きくなる場合がある。そこで、このため、評価関数のパラメータを調整するなどの方法もあるが、それだけでは、高精度なずれ量を計算することが難しい。そこで、複数の評価関数を使用することで、ずれ量誤差を軽減する。複数の評価関数で位置ずれ量を算出した場合、これらは独立に計算される。そのため、それぞれの評価関数による位置ずれ量が異なる値を示す場合もある。それぞれの評価関数ごとの結果に、重み付け処理を行い最終的な、位置ずれ量補正値を算出する（ずれ量補正計算）。
【００５１】
Ｂ．患者位置決め精度向上２
本発明の「自動位置決め法」は、参照ＣＴ画像からＤＲＲ画像を作成し、Ｘ線画像とのずれ量を算出する方法である。一般に、患者位置決めに使用する参照画像は、一般には治療計画で使用したＣＴデータを元に作成した2次元画像(Digital Reconstructed Radiography: ＤＲＲ)が一般的である。この場合、参照ＤＲＲ画像は、元を正せばＣＴ画像であるため、この状況で自動位置決め計算結果のずれ量分、治療台を移動することで、参照ＤＲＲとＸ線画像の位置関係がほぼ一致する。
【００５２】
しかし、粒子線施設では、照射前（照射直前ではなく、多くは一日以上前に実施）に実施するリハーサルというプロセスにて、ＣＴ画像から作成したＤＲＲ画像を参照画像とし、Ｘ線画像を取得して位置決めを実施する。最終的に、リハーサルにて、ＤＲＲ参照画像とそこで取得したＸ線画像を用いて位置決めが実施される。最終的に、この二つの画像の位置関係が一致するＸ線画像が作成され、これが照射時参照画像となる。照射時に自動位置決めを実行する場合の計算上の参照画像は、ＣＴ画像となる。この場合、自動位置決めは、参照ＣＴ画像と撮影Ｘ線画像を用いて、これらのずれ量を算出する。このずれ量分治療台を移動しても、参照Ｘ線画像と撮影Ｘ線画像は一致しない場合がある。これは、リハーサル時に自動位置決めを実施したときに、自動位置決め結果がユーザが納得いく結果でなかったか、またはユーザの判断により、手動微調整（これをΔＲと呼ぶ）を実施したことが原因である。よって、照射時自動位置決めを実施する際には、このΔＲを補正した自動位置決め結果を治療台移動量に反映することで、参照Ｘ線画像と位置関係を一致させることができる。
【００５３】
Ｃ．計算時間短縮化
自動位置決め計算を行う際、数百から数千回のＤＲＲ画像作成を行う。また、このＤＲＲ画像計算は、計算量が多いため、画像1枚作成のための計算時間もかかる。また、複数の評価関数を実行させる場合、評価関数の計算量も増加する。そこで、高速計算させるために、以下の方法を用いた。
・最適化法（勾配法、凸最適化など）の最適化アルゴリズム
・Graphics Processing Unit(GPU)の使用によるマルチスレッド計算
【００５４】
ＤＲＲ画像作成が、位置決め計算の大幅を占めている。このＤＲＲ画像作成の高速化は上記GPUの利用で改善されるが、他のアプローチとして、画像作成数の減少の工夫があげられる。本発明では、画像面内の回転、並びに移動に関しては、ＤＲＲ画像作成時間よりもかなり少ない計算時間であることから、1枚のＤＲＲ画像自体を回転、移動させることで、ＤＲＲ画像作成枚数を減少させている。しかし、画像面内の移動に関しては、大幅に移動を行う場合には、Ｘ線焦点と被写体の位置関係がＤＲＲ画像に影響するため、ＤＲＲ画像を再作成する必要が出てくる。計算最終段階では、ＤＲＲ画像作成は、微小距離、微小角度回転であることから、画像面内の移動、回転に対してＤＲＲ画像再作成を行う必要がなくなり、更に計算時間を短縮化することが可能となる。
【００５５】
Ｄ．時間軸への拡張
４次元ＣＴ撮影による時系列の参照ＣＴ画像、並びに透視装置による時系列の撮影画像を扱う場合が想定される。
【００５６】
参照ＣＴ画像の位相を指定するか、撮影画像の位相を指定するかの2つが考えられるが、どちらの方法でも図１０に示した、静止画像の場合と同じプロセスに、図１４に示す第２実施形態のように位相パラメータを加えて、順次位相を変えながら、評価関数による計算を行う。
【００５７】
計算した結果は次の通りである。
【００５８】
図１３に示した腹部領域の患者データを使用した。図１３中、四角部分が、ずれ量の計算領域である。
【００５９】
１．相互情報量アルゴリズムでずれ量計算
相互情報量の評価関数は、規格化相互情報量を用いた。
【００６０】
２つの画像ＡとＢの離散確率変数をH(A)とH(B)と表すと、
【数１】
0006281849.tif 000002
pTAとpTBは、画像ＡとＢにおける周辺確率を示す。ΩTA,Bは、画像AとBの重複領域を示す。
【００６１】
このときの相互情報量は、次式で表せられる。
【数２】
0006281849.tif 000003
ｐTAB(a,b)は、同時確率密度関数である。
【００６２】
離散確率変数A,Bの同時分布関数H(A,B)で規格化した下式を規格化相互情報量という。
【００６３】
【数３】
0006281849.tif 000004
移動量(Ｘ軸：0.5mm, Ｙ軸：-3.0mm, Ｚ軸：5mm)、回転量(Ｘ軸：5.0度, Ｙ軸：-5.5度, Ｚ軸：11度)の誤差を含む。
【００６４】
２．エッジ検出アルゴリズムでずれ量計算
1/(1+x)の微分値でエッジを強調させる。同様な物としてsobelフィルタ、robertsフィルタなどがある。
【００６５】
移動量(X軸：0.2mm, Y軸：0.0mm, Z軸：0.0mm)、回転量(X軸：0.0度, Y軸：-0.1度, Z軸：0.1度)の誤差を含む。
【００６６】
３．規格化相互情報量アルゴリズムとエッジ検出アルゴリズムを併用し、ずれ量計算
移動量(X軸：0.0mm, Y軸：0.0mm, Z軸：0.0mm)、回転量(X軸：0.0度, Y軸：0.0度, Z軸：0.0度)の誤差を含む。
【００６７】
評価関数部分を含まずＤＲＲ計算のみで、CPUとGPUを用いた場合の計算時間は、ＣＰＵで計算：１１６秒、ＧＰＵで計算：２．８秒となった。ＣＴデータは、５１２×５１２×２５６マトリックスサイズとした。
【００６８】
リハーサルを実施しない場合の参照画像の作成手順を図１５に、リハーサルを実施する場合の参照画像の作成手順を図１６に示す。
【００６９】
次に、ダウンサンプリング法について説明する。
【００７０】
画像レジストレーションの解を高速に求めるため、画像データを低解像度から高解像度へと順次上げていき、計算領域を限定していく。解像度を下げると情報量が減少するため、得られる精度も低下する。ダウンサンプリングを行わない場合、高い精度は得られるが、位置決定までに時間がかかることが想定される。精度を保持したまま、高速に位置決めを行うことができるアルゴリズムが求められる。最適化手法にステージ計算を導入することで、高精度かつ、高速に位置決め可能な手法を実現できる。最適化手法であるＰｏｗｅｌｌ法では、イタレーションを繰り返して徐々に最適解に近づけていく。その過程で、序盤のイタレーションでは大まかな位置合わせが行われ、終盤のイタレーションでは詳細な位置合わせが行われる。そこで、Ｐｏｗｅｌｌ法にステージ計算を導入し、図１７に示す如く、序盤（ステージ１）では高いダウンサンプリングレベルでイタレーションを行い、ステージが進むごとにダウンサンプリングレベルを下げていくように改良を行った。本アプリケーションでは、ＤＲＲ作成後にダウンサンプリングするのではなく、ＤＲＲ作成時にダウンサンプリングして計算することで、計算コスト削減を行っている。
【００７１】
次に、座標計算方法について説明する。
【００７２】
手動位置決め時は、撮影画像を参照画像に一致するように、撮影画像を動かす計算を行う。一方、自動位置決め時は、撮影画像に合うように参照画像（ＣＴ画像データ）を動かすため、手動位置決め時と反対方向の画像合わせとなる。
【００７３】
通常、部屋座標系で計算を行うとき、回転（部屋アイソセンタ＝ターゲットアイソセンタ）→並進の順序である。しかし、自動位置決め時、参照先となる撮影画像は、患者の横たわる位置により、部屋アイソセンタにターゲットアイソセンタが位置しているとは限らない。それ故、回転を行ったときの回転中心が、部屋アイソセンタ＝ターゲットアイソセンタの位置以外で行われることとなり、その結果並進値に誤差が出る。
【００７４】
そこで、自動位置決めでは、図１８に示す如く、並進→回転（並進先位置）とした。また、手動位置決め時と参照画像−撮影画像の参照元の関係が、反対であることから、回転順序は、θ→ψ→Φの順序の反対Φ→ψ→θの順序とする。
【００７５】
このようにすることで、自動位置決め結果を反映するときの回転中心位置は、ターゲットアイソセンタ＝部屋アイソセンタとなる。
【００７６】
次に、参照ＣＴ画像データとのずれ量補正について説明する。
【００７７】
手動位置決めでは、撮影画像を参照画像に合わせる作業を行う。この参照画像は多くはリハーサル時に取得した画像である。この画像は、ＤＲＲ画像（治療計画ＣＴ画像データから作成）とＸ線画像の合わせ込みを行ったものである。ユーザは、ＤＲＲ画像に限りなく近いＸ線画像を取得するが、ＣＴ撮影時とリハーサル時に患者の形状が全く同じではない（Interfractional change）。そして手動位置決めを実施している以上、人の指示誤差が影響することは避けられない（この誤差量をΔＲと呼ぶ。）。一方、ＣＴ画像データを参照画像として使用する自動位置決めでは、ＣＴ画像をＸ線画像に合わせ込む計算を行っている。自動計算の結果を反映して、患者位置を変えＸ線画像を取得しても、このままでは参照画像と撮影画像は完全に一致しない。
【００７８】
そこで、図１９に示す如く、照射時の自動位置決め計算結果に、参照画像とＣＴ画像のずれ量（ΔＲ）を加えることで、照射時自動位置決め結果を治療台に正確に返すことができる。
【実施例】
【００７９】
骨盤領域の患者データを使用して、本件の実施例を示す。
【００８０】
自動位置決め実施前の、参照ＤＲＲ画像とＸ線撮影画像の位置関係を図２０に示す。
【００８１】
正規化相互情報量を評価関数として自動計算した結果を図２１に示す。正面像側面像ともに、位置が合っていないことがわかる。
【００８２】
Gradient Differenceを評価関数として自動位置決め計算を実施した結果を図２２に記す。Gradient Differenceは、異なる画像モダリティー間のずれ量算出には不向きであるため、X線画像の輝度により、計算精度が影響する。
【００８３】
規格化相互情報量とGradient Differenceを評価関数として自動位置決め計算を行った結果を図２３に記す。参照画像とX線撮影画像ともに、位置関係が良好であることがわかる。
【００８４】
本実施形態においては、２次元検出器を照射ポート内に設置したので、構成が簡略であるとともに、患者への心理的負担が少ない。なお、照射ポート以外の場所に２次元検出器を設けることも可能である。Ｘ線画像の撮影方向も垂直方向と水平方向に限定されない。角度も直交２方向に限定されず、任意で良い。
【００８５】
又、前記実施形態においては、照射ポートが垂直方向と水平方向に固定されていたが、本発明の適用対象はこれに限定されず、回転ガントリ内を照射ポートが移動するようにされた放射線治療システムにも同様に適用できる。Dual energy Ｘ線による骨強調画像や組織強調画像にも適応可能である。放射線も、陽子線等の粒子線を含む。
【符号の説明】
【００８６】
１０…患者
１２…患者ベッド
２０、２２…照射ポート
３０、３２…Ｘ線管
３４、３６…２次元検出器
５０…ＣＴ撮像装置
５２…治療計画装置
５４…患者データベース
５６…患者ファイルサーバ
６０…位置決め計算機
６２…機器制御用計算機
【図１】
0006281849.tif 000005
【図２】
0006281849.tif 000006
【図３】
0006281849.tif 000007
【図４】
0006281849.tif 000008
【図５】
0006281849.tif 000009
【図６】
0006281849.tif 000010
【図７】
0006281849.tif 000011
【図８】
0006281849.tif 000012
【図９】
0006281849.tif 000013
【図１０】
0006281849.tif 000014
【図１１】
0006281849.tif 000015
【図１２】
0006281849.tif 000016
【図１３】
0006281849.tif 000017
【図１４】
0006281849.tif 000018
【図１５】
0006281849.tif 000019
【図１６】
0006281849.tif 000020
【図１７】
0006281849.tif 000021
【図１８】
0006281849.tif 000022
【図１９】
0006281849.tif 000023
【図２０】
0006281849.tif 000024
【図２１】
0006281849.tif 000025
【図２２】
0006281849.tif 000026
【図２３】
0006281849.tif 000027