(57)【特許請求の範囲】
【請求項１】
患者の患部に放射線を照射して治療を行う放射線治療時に患者を載せる治療台を位置決めする患者位置決めシステムにおいて、
治療計画時の患者のデータである治療計画時患者データを、異なる種類の測定装置による複数の上記治療計画時患者データとして取得する治療計画時データ取得装置と、
上記複数の治療計画時患者データのうち少なくとも一つの治療計画時患者データを用いて治療計画を行う治療計画装置と、
治療時の患者のデータである治療時患者データを取得する治療時データ取得装置と、
上記複数の治療計画時患者データからそれぞれ特徴量を取得して、上記複数の治療計画時患者データをそれぞれ治療計画時患者特徴量データとし、統一の空間を定義して、上記それぞれの治療計画時患者特徴量データを上記統一の空間に統合化して上記統一の空間で表現した一つの統合化治療計画時患者データとして生成するとともに、上記治療時患者データから特徴量を取得して治療時患者特徴量データとし、この治療時患者特徴量データを上記統一の空間で表現した治療時患者データとして生成し、上記統合化治療計画時患者データと上記統一の空間で表現した治療時患者データとにより、治療計画時の患者の形状と治療時の患者の形状を比較して解析する形状解析装置と、
この形状解析装置により治療計画時の患者の形状と治療時の患者の形状を比較して解析した結果を基に上記治療台の移動量を演算して上記治療台を移動させる治療台位置決め処理装置と
を備えたことを特徴とする患者位置決めシステム。
【請求項２】
上記治療時データ取得装置は、上記治療時患者データを、異なる種類の測定装置による複数の治療時患者データとして取得し、上記形状解析装置は、上記複数の治療時患者データからそれぞれ特徴量を取得して上記複数の治療時患者データをそれぞれ治療時患者特徴量データとし、このそれぞれの治療時患者特徴量データを上記統一の空間に統合化して上記統一の空間で表現した一つの統合化治療時患者データとして生成して、上記統合化治療計画時患者データと上記統合化治療時患者データとにより、治療計画時の患者の形状と治療時の患者の形状を比較して解析することを特徴とする請求項１に記載の患者位置決めシステム。
【請求項３】
上記形状解析装置は、統一の空間において患者の物理モデルを生成し、この物理モデルを用いて患者の形状を解析することを特徴とする請求項１または２に記載の患者位置決めシステム。
【請求項４】
上記物理モデルは、ばねモデルを含むことを特徴とする請求項３に記載の患者位置決めシステム。
【請求項５】
上記統一の空間は電子密度空間であることを特徴とする請求項１から４のいずれか１項に記載の患者位置決めシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、放射線や粒子線などを患者の病巣部に照射して治療を行う放射線治療における、好適な患者位置決めシステムに関する。
【背景技術】
【０００２】
患者位置決めシステムでは、始めに、断層像撮影装置（例えば、Ｘ線ＣＴ(Computed Tomography)装置）を用いて患者の病巣部を撮影した治療計画用３次元ＣＴデータを取得し、このＣＴデータの診断結果に基づいて治療計画を立てる。このとき３次元ＣＴデータに基づいて腫瘍患部の位置や形状を特定し、放射線を照射する方向や照射線量などを決める。次に、決定した治療計画に基づいて放射線治療を行うことになる。しかし、ＣＴ撮影時から放射線治療までの間にかなりの時間が経過していると、治療時における治療台の患者の位置や体位が治療計画作成時の患者の位置や体位と異なっていることが多い。そのため放射線治療を行う前に、現在の患者位置と治療計画時の患者位置のずれを補正する必要がある。
【０００３】
このずれの補正量を算出するために必要な基準画像を、治療計画時の３次元治療計画用データから再構成して生成し、基準ＤＲＲ（Digitally Reconstructed Radiograph：デジタル再構成ラジオグラフ）画像とする。一方、現在の患者位置は、Ｘ線ＴＶ画像撮影装置を用いて取得する。そして、取得したＸ線ＴＶ画像と再構成して得られた基準ＤＲＲ画像とを比較し、画像処理を施すことによって補正量を算出する。治療ビームが患部の適切な位置を照射するように、算出した補正量に基づいて治療台の３次元位置と姿勢を調整する。以上の処理を行う装置が患者位置決めシステムである。こうした患者位置決めシステムでは、患者位置決めの精度向上や速度向上が要望されている。
【０００４】
特に、近年、たとえば粒子線など、体内で線量を集中させることが可能な放射線治療法がある。この治療法では、粒子線ビームのエネルギーを調節して腫瘍の深さ方向の位置に合わせることによって、高い線量部分を腫瘍患部に一致させることが可能である。つまり、腫瘍だけに高い線量を照射しつつ、周囲の正常組織に対する影響を低減できる。この性質を活かすためには、腫瘍患部のみに粒子線を照射するための、高精度な患者の位置決め技術が重要になってくる。
【０００５】
このような患者位置決めシステムにおいては、精度向上や速度向上のために、様々な補助的な装置や方式が用いられている。例えば、特許文献１では、治療計画時に患者の特定の人体箇所の基準座標を決定し、治療時に、まず患者の人体箇所の相対位置決め（例えば姿勢）が、基準座標の相対位置決めと合致するように行われ、その後患者単体を基準座標の絶対位置と合致するよう移動させて位置決めする、という技術が開示されている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特表２００２−５２８１６８号公報
【非特許文献】
【０００７】
【非特許文献１】矢野、外３名、「ＣＴ値から求めた電子密度線量計算への利用について」、日本放射線技術学会誌、昭和５６年１１月、第３７巻第６号、ｐｐ７２９〜７３７
【発明の概要】
【発明が解決しようとする課題】
【０００８】
特許文献１では、患者の回転方向を含む姿勢が、治療時に治療計画時と同じになるように位置決めすれば常に患者内部、すなわち患部も治療計画時と同じ位置になる、という考えで位置決めされる。しかしながら、このような位置決め装置において、治療計画時の腫瘍患部の位置や形状が、時間や患者体位のわずかなずれにより、治療時において異なる可能性がある。その結果、治療計画で決められた放射線を照射する方向や照射線量を精確に実現することが出来ず、治療時に周囲の正常組織に影響を与えてしまう可能性がある。
【０００９】
この発明は、上記のような問題点を解決するためになされたものであり、時間や体位の変化により体内の腫瘍患部が変形しても、周囲正常組織への影響を低減し、精度良く位置合わせできる患者位置決めシステムを得ることを目的としている。
【課題を解決するための手段】
【００１０】
この発明に係る患者位置決めシステムは、治療計画時の患者のデータである治療計画時患者データを、異なる種類の測定装置による複数の治療計画時患者データとして取得する治療計画時データ取得装置と、複数の治療計画時患者データのうち少なくとも一つの治療計画時患者データを用いて治療計画を行う治療計画装置と、治療時の患者のデータである治療時患者データを取得する治療時データ取得装置と、複数の治療計画時患者データからそれぞれ特徴量を取得して、複数の治療計画時患者データをそれぞれ治療計画時患者特徴量データとし、統一の空間を定義して、それぞれの治療計画時患者特徴量データを統一の空間に統合化して統一の空間で表現した一つの統合化治療計画時患者データとして生成するとともに、治療時患者データから特徴量を取得して治療時患者特徴量データとし、この治療時患者特徴量データを統一の空間で表現した治療時患者データとして生成し、統合化治療計画時患者データと統一の空間で表現した治療時患者データとにより、治療計画時の患者の形状と治療時の患者の形状を比較して解析する形状解析装置と、この形状解析装置により治療計画時の患者の形状と治療時の患者の形状を比較して解析した結果を基に上記治療台の移動量を演算して上記治療台を移動させる治療台位置決め処理装置とを備えたものである。
【発明の効果】
【００１１】
この発明によれば、患者の内部組織の形状変化を推定することにより、治療計画時と治療時の腫瘍患部の違いが推定でき、効率的に位置合わせができるとともに、治療における正常組織への影響を低減できる。
【図面の簡単な説明】
【００１２】
【図１】本発明の実施の形態１による患者位置決めシステムの概略構成を示すブロック図である。
【図２】本発明の実施の形態１による患者位置決めシステムの動作の概略を示すチャート図である。
【図３】本発明の実施の形態１による患者位置決めシステムの形状解析装置の動作を説明するチャート図である。
【図４】本発明の実施の形態1による患者位置決めシステムの統一の空間を示すイメージ図である。
【図５】本発明の実施の形態1による患者位置決めシステムのばねモデルの基本図である。
【図６】本発明の実施の形態1による患者位置決めシステムのばねモデルの詳細図である。
【図７】本発明の実施の形態1による患者位置決めシステムのばねモデルを３次元に展開する説明図である。
【図８】本発明の実施の形態1による患者位置決めシステムのラベル化を説明する図である。
【発明を実施するための形態】
【００１３】
実施の形態１．
図１は、本発明の実施の形態１に係る患者位置決めシステムの概略構成を示すブロック図である。患者位置決めシステムは、治療計画時データ取得装置１、治療計画装置２、治療時データ取得装置３、形状解析装置４、患者位置決め処理装置５、治療台６などで構成
される。
【００１４】
治療計画時データ取得装置１は、例えば、Ｘ線ＣＴ（Computed Tomography）装置、ス
テレオカメラで得られる体表面の形状モデルデータ生成装置、体表マーカの位置データを取得する装置、体内に埋め込む事により位置決めに用いられる体内マーカの位置データ取得装置、体内の透視画像であるＸ線ＴＶ画像生成装置、体内の組織状態を取得するための超音波装置、ＭＲＩ（Magnetic Resonance Imaging）、超音波を改良して得られる体内の弾性情報などを測定する装置、ＰＥＴ（ポジトロン放出型断層撮影：射線経路などが画像として表示できるシステム）など、患者の体表や内部組織の状態などが取得できる装置から選択され、患部の部位などに応じて一つ以上の測定装置により構成される。なお、治療計画時データ取得装置１により治療計画時データ取得工程が形成される。
【００１５】
治療計画装置２は、医師などが治療計画時データ取得装置１から得られたデータを用いて、放射線を照射する患部、照射角、照射線量などを計画するための装置である。治療計画装置２により治療計画工程が形成される。
【００１６】
治療時データ取得装置３は、治療計画時データ取得装置１と同様、治療時に患者の状態を取得する装置であり、必ずしも治療計画時データ取得装置１で用いられる装置と同じ構成である必要はない。例えば、治療計画時のデータ取得は治療計画室で行い、治療時のデータ取得は粒子線などの放射線治療を行うための治療室で行う、というように、治療計画時と治療時でデータ取得環境が異なる場合が多く、それぞれで備えられた装置が異なる場合が多い。治療時データ取得装置３により治療時データ取得工程が形成される。
【００１７】
形状解析装置４は、治療計画時データ取得装置１により得られたデータからデータに応じて体表面、病巣部、骨などの位置や形状を特定しモデル化する。得られるデータの点群や、エッジ処理を行う事により得られるデータをモデル化する。また、治療時データ取得装置３により得られるデータも同様にモデル化する。そして、このモデル化されたデータ群（情報）をもとに、患部の病巣部などの形状変化を解析する。形状解析装置４により形状解析工程が形成される。
【００１８】
患者位置決め処理装置５は、形状解析装置４で得られる解析結果を用いて治療台６の移動量を演算するものである。患者位置決め処理装置５により患者位置決め処理工程が形成される。
治療台６は、放射線治療の際に、放射線や粒子線など治療用のビームが患部の適切な位置を照射されるように、３次元位置および姿勢が調整可能な機構を備える。
【００１９】
なお、治療計画装置２、形状解析装置４、患者位置決め処理装置５は、計算機によってそれぞれの機能が実行される。これらの装置は、それぞれが別々の計算機によって実行されても良いし、一つの計算機が複数の装置の機能を実行するように構成されても良い。
【００２０】
次に、以上の患者位置決めシステムの動作の概要を図２に従って説明する。まず、治療計画室で治療計画時データ取得装置１により患者のデータを取得する（ＳＴ１０１）。この取得したデータを治療計画時患者データと呼ぶことにする。次に医師などが、治療計画時データ取得装置１、特にＸ線ＣＴにより取得されたデータを用いて、放射線を照射する患部、照射角、照射線量などを決定する治療計画を実施する（ＳＴ１０２）。以上で治療計画時の工程が終了する。
【００２１】
次に、放射線を照射して治療するに際して、患者を治療台６に載せて治療計画に従って放射線を照射する位置に患者を位置決めする工程に入る。まず患者を治療台６に載せた状態で治療時データ取得装置３により患者データを取得する（ＳＴ１０３）。この取得した
データを治療時患者データと呼ぶことにする。形状解析装置４において、治療計画時患者データと治療時患者データを、それぞれ一つの統一の空間に統合化処理する（ＳＴ１０４）。なお、治療計画時患者データは、治療計画時データ取得装置１により取得した後治療時までの間に統合化処理しておいても良いのは言うまでもない。
【００２２】
次に、形状解析装置４により統合化された治療計画時患者データと治療時患者データをモデル化して比較解析する（ＳＴ１０５）。モデル化および比較解析の内容については、後で詳細に説明する。この比較解析した結果が患者位置決め処理装置５に出力される（ＳＴ１０６）。患者位置決め処理装置５は、受け取った解析結果を基に治療台６の移動量を演算する（ＳＴ１０７）。この演算結果の移動量分治療台６を移動する（ＳＴ１０８）。治療台６の移動後、医師などが、患部が治療計画通りに放射線が照射される位置にあるかどうか、Ｘ線ＴＶ画面などにより確認し、患部の位置ずれが許容範囲内かどうかを判断する（ＳＴ１０９）。位置ずれが許容範囲内であれば放射線治療を実施する。（ＳＴ１１０）。もし位置ずれが許容範囲を超えていれば、治療台６を微調整して（ＳＴ１１１）、再度Ｘ線ＴＶ画面などにより位置ずれを確認して許容範囲内かどうかを判断する（ＳＴ１０９）。
【００２３】
以下で詳細を説明するように、本発明によれば、治療計画時患者データ取得時と治療時患者データ取得時において、患者の姿勢が異なる、あるいは外部力が異なる、などにより、患部近傍の形状が変化してもその変化の原因を推定して、患部の形状が治療計画時と同じになるように治療台６の姿勢などの位置決めを行うため、上記ＳＴ１１０の判断工程において患部の位置ずれが許容範囲を超えることはほとんどない。したがって医師などによるＸ線ＴＶ画面での確認回数は１回で済むため、患者のＸ線被ばくが抑制できる。
【００２４】
以下、形状解析装置４における比較解析の処理、すなわち図２におけるＳＴ１０４〜ＳＴ１０６の処理について説明する。形状解析装置４の内部処理は概略図３のように、データ統合手段７（図２のＳＴ１０４に相当）、データ解析手段８（図２のＳＴ１０５に相当）、データ出力手段９（図２のＳＴ１０６に相当）で構成されている。
【００２５】
まず、データ統合手段７において、異なる測定値のデータを統一した空間により定義する。この統一の空間のイメージ図を図４に示す。例えば、Ｘ線ＣＴ１０により測定された３次元データであるＣＴ値を電子密度で再定義する。ＣＴ値から電子密度を求める方法は例えば非特許文献１に記載された方法などを用いる。この電子密度で再定義した３次元空間を電子密度空間１３と呼ぶ。この電子密度空間１３によれば、体内の脂肪や骨などの密度を３次元空間上の位置ごとに定義できる。ここでは、この電子密度空間１３を統一の空間として用いる。
【００２６】
また、ＭＲＩ１２やＸ線ＴＶ画像１１など、Ｘ線ＣＴとは異なる種類の測定値（モダリティ）の２次元データに対しても、フィルタ処理をすることにより２次元データから３次元情報の候補となる特徴量を取得し、その特徴量を電子密度空間１３として定義した統一の空間の３次元位置に対応した情報として保存する。超音波１５のデータに対しては、プローブに３次元センサを搭載しておくことにより、プローブの位置と取得データの位置情報などが取得でき、電子密度空間１３の位置との一対一対応がとれる。さらにステレオカメラで得られる体表面１４のデータでは、データの３次元位置が分かっているために、電子密度空間１３に対して体表面のデータを容易に対応付けできる。また、上記以外の測定装置であっても患者の体表面や内部組織の２次元データや３次元データが得られる測定装置であれば使用することができる。
【００２７】
このように、統一の空間、ここでは電子密度空間１３、を用いて異なる測定装置の特徴量データを３次元位置に対応付けできる。このようにして、例えば超音波では超音波で取
得しやすい体内の内部情報、ＭＲＩではＭＲＩで取得しやすい乳腺の位置情報など、異なる測定装置によって得られる、それぞれの測定装置での最適な情報を一つの空間で定義することが可能となる。さらに、Ｘ線ＣＴでは乳腺を特定することが難しいが、ＭＲＩによるデータを統一の空間で定義することにより、乳腺の位置がＸ線ＣＴのデータから得られた電子密度空間１３上でも予測が可能となる。
【００２８】
上記の説明では、種々のデータ取得装置を用いたが、実際には、治療部位や準備できるデータ取得装置の種類に応じて、Ｘ線ＣＴとＭＲＩ、Ｘ線ＣＴと超音波とステレオカメラ、といった様に、データ取得装置の種類を決定し、それぞれの組み合わせに応じて適切な処理を行う。
【００２９】
次に、データ解析手段８においては、この統一の空間すなわち電子密度空間１３でデータ解析を行う。例えばこの電子密度空間１３に対してフィルタ処理を行う。また異なるモダリティにおけるデータに対しても、ボリュームデータ（３次元）の特徴解析、２次元データの特徴解析、２次元３次元特徴解析を、それぞれのデータに応じて実施する。例えば３次元のＣＴ、ＭＲＩボリュームデータであれば、表面のエッジデータ、骨データの抽出、内部臓器の形状抽出などをおこなう。Ｘ線ＴＶ画像のような２次元画像データであれば、画像処理による一般的なエッジ処理等による特徴点抽出処理を行う。２次元３次元特徴解析では、２次元データＤＲＲから３次元データの推定を特徴点による解析と解剖学的情報（各臓器や骨などの一般的な位置関係の情報）から算出する。この結果、測定装置で得られるＣＴやＭＲＩデータを電子密度空間１３上に変換できる。この変換により、３次元ボリュームデータ内の表面のエッジデータ、骨データ、内部臓器の形状などを電子密度空間上に対応付けでき、電子密度空間上で表現（定義）できる。例えば、ＣＴとＭＲＩを測定する領域を一致させておくことにより、一方では取得が困難な情報を補間する事が可能となる。そして、他のモダリティである２次元データのＸ線ＴＶ画像から得られる特徴と、２次元３次元情報が得られるＤＲＲでの特徴を、ＣＴやＭＲＩ等で得られる特徴と一致するようにデータ取得することにより、電子密度空間上に様々な測定装置の情報を集約できる。このマルチモダリティ（Ｘ線ＣＴ、ＭＲＩ・・・など複数のデータから得られる特徴を含んだ情報）の空間を潜在空間と定義する。
【００３０】
この得られた特徴点を用いてモデル化を実施する。例えば図５に示すようなばねモデルを定義する。すなわち、体表面の特徴点１６（座標(x1,y1,z1)）および骨の特徴点１７（座標(x2,y2,z2)）と腫瘍の特徴点１９（座標(x3,y3,z3)）の間をそれぞればね１８と仮定し、そのばね係数をそれぞれＫ１３、Ｋ２３とする。例えば、体表面の特徴点１６と腫瘍の特徴点１９の間のばねによる力ＦＫ１,３は、ばね係数Ｋ１３に変位を乗じたものとな
る。また、モデル化の方法は、ここで示す力学モデルの一種であるばねモデルではなく、他の力学モデルや力学モデル以外の物理モデル、例えば血液の流れなどを用いた血管などの流体モデル、電気信号等による電磁気モデルのようなものでも良く、物理モデルであればなんでもよい。ばねモデルを用いる場合では、例えば数（１）式のように外部からの力をＦall（重力、形状変化など全ての因子）としたときに、表面の力Ｆ表面と内部の力Ｆ
内部の和として定義ができる。ここでの内部の力とは呼吸、脈、心臓などの力を示す。一方表面の力は、外部から体表面が押される力などを示している。従って、数（１）式の定常状態（変化がない状態）では左辺が０であると定義ができる。
【００３１】
この数（１）式をばねモデル定義式数（２）式で表す。数（２）式の左辺、右辺の行列の上の矢印は、左辺、右辺がベクトル量であることを示している。左辺の(x,y,z)は体表
面の座標と内部の座標両方を示し、その座標上の力をＦi,jと定義する。ｉ＝ｊ、すなわ
ちＦi,iはｉ番目の特徴点での力を表し、ｉ≠ｊの場合はｉ番目の特徴点とｊ番目の特徴
点の間に働く力を表す。この左辺の力ベクトルは、外部の力が与えられたら、内部の特徴点間ごとの力ベクトルのトータルが外部の力ベクトルのトータルと一致するように、体内
のばねモデルとして反映させる事を示している。数（２）式の右辺では、ばねモデルの定義により、外部の力が与えられる特徴点や領域から、その周辺の特徴点に力が加わり（派生し）、周辺の特徴点の力に反映されるように定義する。すなわち、数（２）式の右辺はモデル化したばねによる力をＦＫi,jとして表しており、例えば図５の腫瘍（ｉ＝３）と
体表面（ｉ＝１）の間に働くばねの力をＦＫ1,3として表している。右辺の対角要素ＦＫi,iは自己相関でありばねの力が無いため通常は０となる。あるいは、次のように対角要素だけで表してもよい。すなわち、対角要素におけるばねの力としての要素であるＦＫ3,3の座標(x3,y3,z3）が、図６の腫瘍の位置であるとすると、ＦＫ3,3は近傍のばねの集合として数（３）式のように定義してもよい。数（３）式における引数は移動量（変位）を表している。このように、力を対角要素に集合することで、式を簡単化できるとともに、非剛体の変形する部位に対して、ある領域に外力が加わったときに周辺に与える影響をモデル化することが可能となる。
【００３２】
このように定義することにより、一旦ばね係数Ｋが求まれば、外部の力やずれ量（例えば体の向きが変わる事により重力の向きが変わり、体表面の形状が変るという外部の力などを重力モデルより算出する）を算出することにより内部の形状が推定できる。
【００３３】
【数１】
0005495886.tif 000002
【数２】
0005495886.tif 000003
【数３】
0005495886.tif 000004
【００３４】
ここでばね係数Ｋの求め方の一例を説明する。まず、治療計画で式（１）の右辺＝０という状態の式を算出する。この＝０というのはトータルＦall＝０を意味し、なにも力が
加わっていない状態である定常状態を表す。つまり治療計画で得られる患者の参照データであり、ベッドに横たわっている状態を表す。この状態で、特徴などを算出し、未知のパラメータであるばね係数Ｋを計算する。計算には未知のパラメータ数の方程式が必要となる。たとえば、定常状態から少し体を動かして（重心を動かして）ある特徴位置(x,y,z)
が(x’,y’,z’)へと移動した事が分かれば、方程式を得ることができる。未知のパラメ
ータ数の方程式を得ることができると、内部の未知のパラメータＫを算出することができ
る。この方法によれば、未知のパラメータ数の方程式を得るため、何回もデータを取る必要があり、ＣＴやＭＲＩの場合は、実用化が難しい。しかし、ＣＴを一回取得し、他のセンサ、例えば、超音波や体表面のデータを時間軸で取得することにより、体位の変化、呼吸の振幅などを取得することができ、複数の特徴点の移動を算出できる。このような方法で、例えば体表面だけのデータを複数取得することにより、ばねパラメータＫの未知数分の方程式が得られ、ばねパラメータＫの値が算出できる。
【００３５】
また、勾配法やScale Conjugate Gradient(パラメータを含んだ勾配法の最適化手法)、ＩＫ（Inverse Kinematics）などを用いることにより、重心などを動かさずに、位置移動のデータを取得することも可能である。このような方法で、治療計画データから数（１）式＝０を満たす方程式を生成して、ばね係数などの値を算出することもできる。
【００３６】
以上では、方程式を得ることでばね係数を求めたが、ばね係数は、ユーザが定義したり、解剖学的な情報をも用いて定義することも可能である。あるいは、ニューラルネットワークなどのモデルと、他のモダリティで得る情報からＫを最適化してもよい。または、Ｘ線ＣＴより変換されて得られる電子密度空間１３では密度に応じてばね係数を定義することもできる。例えば、電子密度が高い領域ではばね定数を大きくし、電子密度が小さい領域では、ばね定数を小さくする。このことにより、ばねモデルの係数が、電子密度に応じた最適なパラメータで評価できるとともに、異なったモダリティにおける情報を統一して評価が可能となる。そして図６の骨や体表面などは特徴解析で得られた特徴点を用いる。またここでのばね方程式だけでなく、重力、体位の変化などの物理モデルで表現してもよい。この処理により、患者の体位の違いにより変形する内部組織をモデル化できる。この結果容易に内部の形状を推定できる。
【００３７】
次に、治療時にデータを取得する。このデータから体の重心の位置などが取得できるとともに、特徴が取得できる。これを先ほど算出した数（１）式に代入することにより、取得できなかった内部の状態の特徴の位置を求めることができる。たとえば、図６で示すように、ｎ点の特徴点が求まる（図６では特徴点は７点を図示しているが、図示しない特徴点があるとして７点目をｎで表記している。）とする。そのｎ点の座標を(x1,y1,z1)・・・(xn,yn,zn)と定義し、その間にあるばね係数をそれぞれＫi,jとすると（i≦n,ｊ≦n）
となる。そこで数（２）式の左辺では、実際は対角の要素のみ値が入り対角以外の要素の値は０となる。しかし、治療計画で得られた特徴の位置に対して、現在の値の特徴の位置が移動したとすると、この移動量∇分のばねの力をＦＫi,jにそれぞれ代入して数（２）
式の右辺の力の要素を求める。
【００３８】
ここで問題になるのは、現在の特徴点と治療計画の特徴点が実際に特徴と一致しているのかである。そのため、容易に一致できる体表面の特徴点などを予め例えばＩＣＰアルゴリズムなどを用いて算出する必要がある。その結果、大まかな、特徴点の移動が算出できるとともに詳細な内部臓器の移動量などを、数（２）式などから、最適化を行う事により算出できる。また計算を簡易化するために、数（２）式の右辺の対角上の力を数（３）式のように簡易的に近似してもよい。数（２）式の対角要素は０である可能性があり、最適化するのに難しい問題があるが、以上のように処理すればこの問題が解決できる。
【００３９】
以上の処理を治療計画時データ取得装置１と治療時データ取得装置３で得られるデータ両方におこなう。その結果得られる図７のような３次元空間を用いて腫瘍や、体の変化を推定する。たとえば潜在空間の体表面のデータ１４を用いてＩＣＰ(Iterative Closet Points)マッチングを行う。その結果、電子密度空間１３ではその対応点が算出でき、ばねモデルが定義できる。このように他の潜在空間（ここでの潜在空間とは電子密度を生成するＸ線ＴＶ画像、Ｘ線ＣＴ、超音波などの測定装置における特徴空間）で変化を評価する事により様々な内部臓器を評価できる。これにより、体内の内部の特徴点の間をばねモデルで表現する事により、治療時に体表面や大まかな情報しかえられない状況においても内部の臓器の状態を容易に推定できるようになる。そしてデータ出力手段９で臓器の変化結果について算出する。
【００４０】
この内部の情報をもとに、患者位置決め処理装置５にて移動量を算出する。例えば、臓器の変形がなければ、治療計画時の状態と治療時の状態に変化が少ない。そのため特許文献１のような従来の位置合わせ方法を用いて治療計画通りに治療が可能である。しかし、照射ターゲットによっては、図７のように治療計画時から治療時にかけて形状が変化する図７の格子ばねモデルでは、表面しか図示していないが、体表面だけでなく内部もばねモデルとして表現する。そこで、例えば解析結果が図７の治療計画時格子ばねモデル２２と治療時格子ばねモデル２３のように得られたら、治療時格子ばねモデル２３から治療計画時格子ばねモデル２２のモデルになるように、外部から力を与える。あるいは、外部から力を与えるのではなく、治療台６を移動させることにより姿勢を変化（体の重心移動）させて、治療時格子ばねモデル２３が治療計画時格子ばねモデル２２のモデルになるようにする方法もある。また、これら外部から力を与えるのと姿勢を変化させるのとを組み合わせることにより治療時格子ばねモデル２３を治療計画時格子ばねモデルに一致させる方法もある。
【００４１】
このように、内部臓器をモデル化することにより、外部から与える力や、治療台６の姿勢などによりどのように内部の形状が変化するかを推定でき、この推定結果を基に、患者位置決め処理装置５にて、治療時格子ばねモデルが治療計画時格子ばねモデルになるような治療台６の移動量を演算することができる。この結果、患部が治療計画によって決めた放射線を照射するための所定の位置になるよう位置決めをすることができる。
【００４２】
乳がんなど形状変化が大きくなる可能性がある対象の場合、変形する対象に対してカップなどを用いて固定する方法を用いてもよい。さらにカップに吸引する機能などがあれば、ＩＫで得られた制御値を、その吸引の角度や吸引力等に反映させ、治療計画時の形状になるようにする。この内部の状態を確認するためには、ＤＲＲ(Digital Reconstruct Radiograph)をＸ線ＣＴから生成して確認を行う。ただし、乳がんなどの形状をＸ線ＣＴから生成したＤＲＲで表示する事は難しい。そこで、他のモダリティ、例えば、乳腺などの情報が得られるＭＲＩ１２の測定データから、統一の空間である電子密度空間を用いて位置情報を取得する。取得する方法として、図８で示すように、例えば、ＭＲＩで得られる対象部位の３次元位置情報（ラベル２４）を保存し、保存した位置が電子密度空間１３に写像した時にその写像空間上の点（ラベル２５）がＭＲＩ上での座標の特徴点を保持する。この処理により電子密度空間やＭＲＩでないモダリティ上で、例えばＤＲＲなどで、乳腺などの位置情報が取得できる。この位置情報を電子密度空間で得るための処理をラベル化と定義する。そのラベル化された特徴点を通過する領域が強調されるようにＤＲＲを生成して形状変化を確認する。このような処理を行うことにより、内部の形状が一致するまで、最適計算をおこなう。以上の結果から、移動量（ターゲットの形状を治療計画に一致できる移動量）を治療台６の姿勢と位置の移動量に反映させる。
【００４３】
以上の処理により、治療計画時と治療時との腫瘍患部近傍の形状を比較解析し、その変化から姿勢などの変化が推定でき、効率的に治療計画で決定した放射線照射ができるように患部の位置合わせができるとともに、治療における正常組織への影響を低減できる。
【００４４】
実施の形態２．
実施の形態１においては、治療時患者データも、治療計画時患者データと同じ統一の空間で処理を行った。必ずしもこのような処理を行う必要はなく、例えば治療計画時患者データに対して図４の統一の空間を生成する。一方、治療時データは統一の空間で処理を行うのではなく、取得した潜在空間のみで治療計画時データと比較してずれ量を算出する。
この算出方法は、テンプレートマッチングや、ＩＣＰなど一般的に用いられている手法で行う。そして、このずれ量は特徴点であるノードに対応して算出されるため、統一の空間である電子密度空間におけるばねモデルを用いて形状の変化を推定でき、治療台の移動量を演算できる。この方法によれば、計算時間が短縮できるとともに、治療時データが少ないモダリティであっても、当該モダリティからずれ量を算出することで形状の変化を推定することができ、治療台の移動量を演算することができる。
【００４５】
実施の形態３．
実施の形態１では、統一の空間を、治療計画時データ取得装置１にて得られるデータを用いて生成した。実施の形態３では、直接治療計画時取得データから統一の空間を生成するのではなく、治療計画装置から得られるデータを用いて統一の空間を生成する。すなわち、治療計画装置においては、治療計画時患者データに、例えばがん治療領域などを特定してその癌治療領域の情報を加えて治療計画データとする。このため、治療計画装置から得られるデータを用いて統一の空間を生成すれば、治療計画の情報を付加した状態で形状の変化が推定できる。これにより、患部の特定が容易になり、患部付近の変形をより精度良く推定することができるため、実施の形態１よりもさらに精度良く治療台の移動量を演算できる。
【００４６】
実施の形態４．
実施の形態１ではＸ線ＣＴから得られる電子密度の空間を統一の空間としたが、統一の空間を、電子密度空間ではなく、治療の部位に合わせて別の空間で定義してもよい。治療の部位によっては、Ｘ線ＣＴのデータから生成された電子密度空間よりも、他のモダリティのデータから生成した方が適当な場合がある。たとえば、ＭＲＩによる測定データによれば乳腺などの情報が良好に得られるため、乳がんの治療においてＭＲＩ１２のモダリティのデータを統一の空間とすれば患部の特定が容易になる。このように、電子密度空間に限らず、それぞれの目的に合わせて統一の空間を定義することができる。この結果、腫瘍が特定し易いモダリティを基準とする事ができ、最適な内部推定を可能とすることができる。
【００４７】
実施の形態５．
実施の形態１ではＸ線ＣＴの測定データを加工して例えば電子密度のような特定の特徴量により統一の空間を定義したが、測定データを加工して統一の空間を定義するのではなく、Image Fusion（画像融合）を用いてもよい。Image Fusionとは、異なる性質の画像を融合して一つの画像とするものである。例えば、Ｘ線ＣＴにより得られるデータを画像データとして、ＭＲＩにより得られるデータも画像データとして、これら二つの画像データを画像データのまま融合して一つの画像データとすることにより統一の空間を定義することができる。この融合した画像から特徴点を設定してばねモデルとしてモデル化を行う。この結果、複数のモダリティで得られる異なる特徴量のデータを統合する事ができ、データ量の削減が可能である。
【００４８】
実施の形態６．
得られる全ての特徴点に対してSelf Organization Map（ＳＯＭ、自己組織化マップ）
等によって自動的にグループ化を行ったり、骨や各臓器の一般的な相対位置関係を表す情報である解剖学的情報を用いて最適化されたモデルを用いてもよい。この結果計算量を削減することが可能となり、実用的な計算速度での形状の変化の推定を可能とする。
【００４９】
実施の形態７．
実施の形態１におけるばねモデル（数（１）式、数（２）式）において、単純なばねモデルではなく、各位置におけるばね係数以外の情報を含めた情報を用いてモデル化してもよい。たとえば、重力や摩擦の影響を考慮したモデル化や、呼吸、心臓、脈による影響を
含めてモデル化するために、時間による変化を含めたモデル化を行ってもよい。時間による変化を含めたモデル化とは、４次元ＣＴなど時間をも一つのパラメータとして動画のように取得したデータによりモデル化を行うことである。さらに解剖学的情報を数式化してもよい。
【００５０】
例えば、呼吸による影響を含める場合、数（１）式においてばね係数Ｋが位置によって決まるだけではなく、呼吸サイクルの時間ｔを含み、時間ｔの変化により特徴点(x,y,z)
の位置が変り、変位が生じることによりＦＫという力の値が得られる。この結果、患者内部の詳細なモデル化が時間軸を含んだ形で推定が可能となる。
【００５１】
実施の形態８．
実施の形態１においては、形状が変化するターゲットに対してＤＲＲ画像を用いて変形を確認したが、ＤＲＲ画像を用いずに他のモダリティにより変形を確認する処理を行うことができる。例えばＭＲＩ１２や超音波１５などにより、確認処理を行ってもよい。ＭＲＩや超音波によって得られるデータも統一の空間に対応して位置付けすることができるため、これらＸ線画像以外のモダリティによっても変形を確認することができる。Ｘ線画像により確認を行うのは、確認のためだけで患者にＸ線被ばくを与えてしまうが、ＭＲＩや超音波を用いれば、確認処理のために患者にＸ線被ばくを与えることがないというメリットがある。また、確認処理のために持ち込める計測装置に応じて位置合わせが可能となる。
【００５２】
実施の形態９．
実施の形態１では、治療計画時データ取得装置１により治療計画時患者データを取得して、この治療計画時患者データを用いて治療計画装置２により治療計画を実施した。この治療計画を、形状解析装置４から得られるばねモデル２３等のデータを生成した後、この生成したデータを用いて実施してもよい。この結果、例えば患者のわずかな向きのずれによるターゲット患部の変形を推定する事ができ、予めこのずれを反映して治療計画が立てられる。
【００５３】
実施の形態１０．
実施の形態１においては、治療時格子ばねモデル２３を治療計画時格子ばねモデル２２に一致するように力を加えたり、治療台を移動したりした。しかし本実施の形態１０ではこの変形の処理を、治療計画時格子ばねモデル２２から治療時格子ばねモデル２３に一致するようにして行う。すなわち、実施の形態１と逆の処理を行う。この処理により、例えばＤＲＲを生成するために必要なＸ線ＣＴ１０を、治療時に取得する必要がなくなり、Ｘ線ＴＶ画像１１などで行えるようになる。これにより、取得データの処理や装置の数を削減する事が可能である。ただし本実施の形態１０を実施するためには、実施の形態９の処理をおこない、患者の向きのずれ量が治療計画に考慮されているようにしなければならない。
【００５４】
実施の形態１１．
実施の形態１において、治療計画時格子ばねモデル２２と治療時格子ばねモデル２３に用いられる特徴点の数は計算コストに依存する。そこで、それぞれの、ばねモデルに用いる特徴点は一般的な画像圧縮方法である画像ピラミッドを生成して特徴量を圧縮して得られる特徴点をばねモデルのノードに用いてもよい。または、得られる特徴点を全て算出し、その特徴点なかで、特徴点らしい特徴点を一般的に用いられる固有値計算により算出し固有値の高い特徴点を用いて評価を行ってもよい。
【００５５】
以上各実施の形態で説明した処理は、ＧＰＵ（Graphic Processing Unit）などを用い
て並列計算処理を行うことができ、高速に処理を行うことができる。
【符号の説明】
【００５６】
１：治療計画時データ取得装置１        ２：治療計画装置２
３：治療時データ取得装置３            ４：形状解析装置４
５：患者位置決め処理装置５            ６：治療台
１３：統一の空間
【図１】
0005495886.tif 000005
【図２】
0005495886.tif 000006
【図３】
0005495886.tif 000007
【図４】
0005495886.tif 000008
【図５】
0005495886.tif 000009
【図６】
0005495886.tif 000010
【図７】
0005495886.tif 000011
【図８】
0005495886.tif 000012