(57)【特許請求の範囲】
【請求項１】
医用画像撮影装置により撮影された画像であって被検体の注目部位を示す医用画像と、前記注目部位の疾患をシミュレーションするためのモデルであって前記注目部位を示す部位モデルとを位置合わせをして合成し、前記医用画像及び前記部位モデルの合成画像を生成して表示する手段を備え、
前記部位モデルは、モデルデータベースに記憶されており、機械バネ定数、静止膜電位及びイオン濃度の情報を含んでいることを特徴とする医用画像処理装置。
【請求項２】
前記生成して表示する手段は、
前記医用画像の複数の特異点及び前記部位モデルの複数の特異点を抽出する手段と、
抽出した前記医用画像の複数の特異点及び前記部位モデルの複数の特異点に基づいて前記医用画像の注目部位と前記部位モデルの注目部位との位置合わせを行う手段と、
を具備していることを特徴とする請求項１記載の医用画像処理装置。
【請求項３】
前記位置合わせを行う手段は、
抽出した前記医用画像の複数の特異点及び前記部位モデルの複数の特異点に基づいて、前記医用画像の注目部位と前記部位モデルの注目部位との大きさを一致させる手段と、
抽出した前記医用画像の複数の特異点及び前記部位モデルの複数の特異点に基づいて、前記医用画像の注目部位と前記部位モデルの注目部位との形状を一致させる手段と、
を具備していることを特徴とする請求項２記載の医用画像処理装置。
【請求項４】
前記大きさを一致させる手段は、抽出した前記医用画像の複数の特異点から各々の特異点間の離間距離及び抽出した前記部位モデルの複数の特異点から各々の特異点間の離間距離を求め、求めた前記医用画像の各々の前記離間距離と求めた前記部位モデルの各々の前記離間距離とがそれぞれ同じになるように前記医用画像の注目部位と前記部位モデルの注目部位との大きさを一致させることを特徴とする請求項３記載の医用画像処理装置。
【請求項５】
前記形状を一致させる手段は、前記医用画像の注目部位の輪郭及び前記部位モデルの注目部位の輪郭を抽出し、抽出した前記医用画像の複数の特異点及び前記部位モデルの複数の特異点に基づいて、抽出した前記医用画像の注目部位の輪郭と前記部位モデルの注目部位の輪郭とが同じになるように前記医用画像の注目部位と前記部位モデルの注目部位との形状を一致させることを特徴とする請求項３又は４記載の医用画像処理装置。
【請求項６】
生成した前記合成画像に前記注目部位の異常部分を示す手段を備えることを特徴とする請求項１乃至５のいずれか一に記載の医用画像処理装置。
【請求項７】
前記生成して表示する手段は、生成した前記合成画像を表示する表示部を備えることを特徴とする請求項１乃至６のいずれか一に記載の医用画像処理装置。
【請求項８】
医用画像を撮影する医用画像撮影装置と、
請求項１乃至７のいずれか一に記載の医用画像処理装置と、
を備えることを特徴とする医用画像診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、医用画像を処理する医用画像処理装置及びその医用画像処理装置を備える医用画像診断装置に関する。
【背景技術】
【０００２】
医用画像診断装置は、医用画像を撮影する医用画像撮影装置やその医用画像を処理する医用画像処理装置等を備えている。医用画像撮影装置としては、例えば、Ｘ線断層撮影装置（Ｘ線ＣＴ装置）、核磁気共鳴装置（ＭＲＩ装置）及び超音波診断装置等が挙げられる。Ｘ線断層撮影装置は、被検体にＸ線を照射し、その被検体を通過したＸ線を検出し、被検体の内部（注目部位）をＣＴ画像として表示する装置である。超音波診断装置は、被検体に超音波を送波し、その反射波（エコー波）を受波し、被検体の内部（注目部位）を超音波画像として表示する装置である。
【０００３】
前述の医用画像を処理する画像処理技術としては、様々な画像処理技術が提案されている。例えば、Ｘ線ＣＴ装置及び超音波診断装置により得られた同一関心部位（注目部位）についてのＣＴ画像と超音波画像を合成して表示する複合型医用画像診断装置や、磁気共鳴画像の特定的な特徴部に対して超音波画像の対応する特徴部を位置合わせし、超音波画像を重ね合わせた磁気共鳴画像を生成する画像診断支援システムが提案されている（例えば、特許文献１及び特許文献２参照）。
【０００４】
一方、被検体の注目部位の疾患をシミュレーションする技術も提案されている。例えば、心臓の形状をモデル化し、そのモデルに心筋の機械的パラメータを組み込むことによって、心臓の動き（壁運動）や心腔内の血液の動態をシミュレーションする技術や、同様に、心臓の形状をモデル化し、そのモデルに心筋の電気パラメータを組み込むことによって、興奮伝播過程や電位変化をシミュレーションする技術が提案されている。このようなシミュレーションにより得られた部位モデル（例えば心臓モデル）からは、注目部位のある種の機能情報が得られる。また、医用画像撮影装置により得られた医用画像（例えば心臓画像）からは、注目部位の形態情報や別の機能情報が得られる。
【特許文献１】特開平１０−１２７６２３号公報
【特許文献２】特開２００３−１５３８７７号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、医用画像撮影装置により得られた医用画像と、シミュレーションにより得られた部位モデルとは別画像となるため、注目部位の形態情報とその機能情報との関係認識が難しく、医師等の診断者が被検体の注目部位の状況を的確に把握することは困難である。
【０００６】
本発明は、上記に鑑みてなされたものであり、被検体の注目部位の形態情報とその機能情報との関係認識を容易にし、診断者が注目部位の状況を的確に把握することができる医用画像処理装置及び医用画像診断装置を提供することである。
【課題を解決するための手段】
【０００７】
本発明の実施の形態に係る第１の特徴は、医用画像処理装置において、医用画像撮影装置により撮影された画像であって被検体の注目部位を示す医用画像と、注目部位の疾患をシミュレーションするためのモデルであって注目部位を示す部位モデルとを位置合わせをして合成し、医用画像及び部位モデルの合成画像を生成して表示する手段を備え、部位モデルは、モデルデータベースに記憶されており、機械バネ定数、静止膜電位及びイオン濃度の情報を含んでいることである。
【０００８】
本発明の実施の形態に係る第２の特徴は、医用画像診断装置において、医用画像を撮影する医用画像撮影装置と、前述の第１の特徴に係る医用画像処理装置とを備えることである。
【発明の効果】
【０００９】
本発明によれば、被検体の注目部位の形態情報とその機能情報との関係認識を容易にし、診断者が注目部位の状況を的確に把握することができる医用画像処理装置及び医用画像診断装置を提供することができる。
【発明を実施するための最良の形態】
【００１０】
本発明の実施の一形態について図面を参照して説明する。
【００１１】
図１に示すように、本発明の実施の形態に係る医用画像診断装置１は、被検体の注目部位（例えば、心臓）を示す医用画像を撮影する医用画像撮影装置２と、その撮影した医用画像を保管する医用画像データベース３と、注目部位の疾患をシミュレーションするための部位モデル（例えば、心臓モデル）を保管するモデルデータベース４と、医用画像及び部位モデルの合成画像を生成する合成処理を行う医用画像処理装置５とを備えている。これらの各部はＬＡＮ（Local Area Network）等のネットワーク６により接続されている。
【００１２】
医用画像撮影装置２は、被検体の注目部位（関心部位）を示す医用画像を撮影する撮影装置である。この医用画像撮影装置２としては、例えば、Ｘ線断層撮影装置（Ｘ線ＣＴ装置）や超音波診断装置等を用いる。Ｘ線断層撮影装置を用いた場合には、医用画像はＣＴ画像であり、超音波診断装置を用いた場合には、医用画像は超音波画像（超音波エコー画像）である。
【００１３】
医用画像データベース３は、医用画像撮影装置２により得られた医用画像（例えば、被検体の各部位に関するＣＴ画像であるボリュームデータや超音波画像等）を記憶する記憶装置である。ここで、被検体の各部位としては、例えば、心臓、肺及び胃等が挙げられる。なお、ＣＴ画像（ボリュームデータ）や超音波画像は、医用画像撮影装置２により得られ、ネットワーク６を介して医用画像データベース３に格納される。
【００１４】
モデルデータベース４は、注目部位の疾患をシミュレーションするためのモデルであって、その注目部位を示すシミュレーションモデルを部位モデルとして記憶する記憶装置である。注目部位が心臓である場合には、心臓の形状がモデル化され、心臓モデルが作成される。その心臓モデルの情報がモデルデータベース４に記憶されている。この心臓モデルの情報としては、例えば、形態、機械バネ定数、静止膜電位及びイオン濃度等の設定値情報が含まれる。
【００１５】
ここで、心臓の形状モデルとしては、セルモデル（数万個のセルを心臓の形状に配置したモデル）や四面体分割モデル（色々な形状の四面体の組み合わせにより心臓を表現したモデル）等が用いられる。例えば、Ｘ線断層撮影装置により撮影された心臓の解剖画像から心臓領域を抽出し、自動的にセル分割又は四面体分割を行うことにより、個人の心臓形状が表現され、個人の心臓形状を有する心臓モデルが用いられる。また、個人の心臓形状ではなく、標準の心臓形状を有する心臓モデルが用いられる場合もある。
【００１６】
医用画像処理装置５は、図２に示すように、各部を集中的に制御するＣＰＵ（Central Processing Unit）等の制御部１１と、ＲＯＭ（Read Only Memory）やＲＡＭ（Random Access Memory）等のメモリ１２と、医用画像や部位モデル等の各種画像を表示する表示部１３と、操作者からの入力操作を受け付ける操作部１４と、各種プログラムや各種データ等を記憶する記憶部１５と、外部装置との通信を行う通信部１６と、医用画像や部位モデル等の各種画像を処理する画像処理部１７とを備えている。これらの各部は、バスライン１８により電気的に接続されている。
【００１７】
制御部１１は、記憶部１５に記憶された各種プログラムや各種データ等に基づいて、各部を制御する。特に、制御部１１は、各種のプログラムやデータに基づいて、データの計算又は加工等を行う一連のデータ処理、医用画像及び部位モデルの合成画像を生成する合成処理及び画像を表示する表示処理等を実行する。
【００１８】
メモリ１２は、制御部１１が実行する起動プログラム等を記憶するメモリであって、制御部１１のワークエリアとしても機能するメモリである。なお、起動プログラムは、医用画像処理装置５の起動時に制御部１１により読み出されて実行される。
【００１９】
表示部１３は、二次元画像や三次元画像等の各種画像をカラー表示する表示装置である。この表示部１３としては、例えば、液晶ディスプレイやＣＲＴ（Cathode Ray Tube）ディスプレイ等を用いる。
【００２０】
操作部１４は、操作者により入力操作される入力部であって、画像表示の開始や画像の切り替え、設定の変更等の各種の入力操作を受け付ける入力部である。この操作部１４としては、例えば、マウスやキーボード等の入力デバイスを用いる。
【００２１】
記憶部１５は、各種のプログラムやデータ等を記憶する記憶装置であって、特に、医用画像と部位モデルとの合成（位置合わせ）に関する合成用データＤ１を記憶する記憶装置である。この記憶部１５としては、例えば、磁気ディスク装置や半導体ディスク装置（フラッシュメモリ）等を用いる。なお、合成用データＤ１は、医用画像と部位モデルとの合成に必要な情報、例えば伸縮量や変形量等の情報である。この合成用データＤ１は、画像処理部１７により得られ、バスライン１８を介して記憶部１５に格納される。
【００２２】
通信部１６は、ＬＡＮやインターネット等のネットワーク６を介して外部装置との通信を行う装置である。この通信部１６としては、ＬＡＮカードやモデム等を用いる。外部装置としては、医用画像撮影装置２、医用画像データベース３及びモデルデータベース４等が挙げられる。
【００２３】
画像処理部１７は、図３に示すように、時相確認手段１７ａと、医用画像情報取得手段１７ｂと、モデル情報取得手段１７ｃと、特異点抽出手段１７ｄ、１７ｅと、内壁及び外壁の輪郭抽出手段１７ｆ、１７ｇと、画像及びモデル間の特異点フィッティング手段１７ｈと、画像及びモデル間の輪郭フィッティング手段１７ｉと、画像及びモデルの位置合わせ手段１７ｊとを具備している。なお、画像処理部１７は、ソフトウェア又はハードウェア（回路）あるいはそれらの両方により構成されている。ここから、注目部位が心臓である場合について説明する。すなわち、医用画像が心臓画像（心臓形態画像）Ｇ１となり、部位モデルが心臓モデルＭ１となる（図４参照）。
【００２４】
時相確認手段１７ａは、心臓画像Ｇ１と心臓モデルＭ１との心電図における時相を確認し、それらの時相が同じ時相となるように時相を調整する手段である。同じ時相に調整することにより、心臓画像Ｇ１と心臓モデルＭ１の形状を精度良く合わせこむことが可能になる。時相としては、例えば、拡張期の特定の時相、特に拡張末期を用いる。心臓の形状モデルを拡張末期の時相で構築した場合、拡張末期では、心臓画像Ｇ１と心臓モデルＭ１とが同じ状態（略静止状態）になるため、この拡張末期の心臓画像Ｇ１及び心臓モデルＭ１を用いた場合、それらの形状をより精度良く合わせこむことが可能である。
【００２５】
医用画像情報取得手段１７ｂは、特定の時相（例えば、拡張末期）における心臓画像Ｇ１を医用画像データベース３から取得する手段である。同様に、モデル情報取得手段１７ｃは、特定の時相（例えば、拡張末期）における心臓モデルＭ１をモデルデータベース４から取得する手段である。なお、心臓モデルＭ１の情報としては、例えば、形態、機械バネ定数、静止膜電位及びイオン濃度等の設定値情報が取得される。
【００２６】
この医用画像情報取得手段１７ｂにより、例えば、図４及び図５に示すような心臓画像Ｇ１が医用画像データベース３から取得される。ここで、図５では、図４の心臓画像Ｇ１を２次元的に示している。また、モデル情報取得手段１７ｃにより、例えば、図４及び図６に示すような心臓モデルＭ１がモデルデータベース４から取得される。ここで、図６では、図４の心臓モデルＭ１を２次元的に示している。
【００２７】
特異点抽出手段１７ｄは、心臓画像Ｇ１から特異点を選択して抽出する手段である。同様に、特異点抽出手段１７ｅは、心臓モデルＭ１から特異点を選択して抽出する手段である。この特異点としては、例えば、心尖部、弁の位置、大動脈、大静脈、冠動脈、心室中隔及び心房中隔等が挙げられる。なお、心臓画像Ｇ１上においては、特異点は自動抽出され、心臓モデルＭ１上においては、特異点は予め設定されおり、前述の設定値情報から得られる。
【００２８】
この特異点抽出手段１７ｄにより、図５に示すように、心臓画像Ｇ１では、例えば、複数の特異点Ａ１〜Ａ４が抽出される。また、この特異点抽出手段１７ｅにより、図６に示すように、心臓モデルＭ１では、例えば、複数の特異点Ｂ１〜Ｂ４が得られる。ここで、各特異点Ａ１〜Ａ４、Ｂ１〜Ｂ４は、三次元の座標（ｘ、ｙ、ｚ）で示される。なお、特異点Ａ２及び特異点Ｂ２がそれぞれ原点とされている。
【００２９】
内壁及び外壁の輪郭抽出手段１７ｆは、心臓画像Ｇ１から心臓の内壁及び外壁の輪郭を抽出する手段である。同様に、内壁及び外壁の輪郭抽出手段１７ｇは、心臓モデルＭ１から心臓の内壁及び外壁の輪郭を抽出する手段である。なお、心臓画像Ｇ１上においては、それらの輪郭は自動抽出され、心臓モデルＭ１上においては、それらの輪郭は予め設定されおり、前述の設定値情報から得られる。
【００３０】
画像及びモデル間の特異点フィッティング手段１７ｈは、心臓画像Ｇ１及び心臓モデルＭ１間において、各特異点Ａ１〜Ａ４、Ｂ１〜Ｂ４をフィッティングする手段である。このとき、心臓画像Ｇ１の各特異点Ａ１〜Ａ４が動かされても良いし、心臓モデルＭ１の各特異点Ｂ１〜Ｂ４が動かされても良い。これは、操作部１４に対する操作者の入力操作に応じて予め設定されている。このフィッティング手段により、特異点Ａ１〜Ａ４、Ｂ１〜Ｂ４の移動情報（例えば、回転量や伸縮量等を示す変換係数）が求められる。この特異点Ａ１〜Ａ４、Ｂ１〜Ｂ４の移動情報は、合成用データＤ１として記憶部１５に格納される。
【００３１】
この画像及びモデル間の特異点フィッティング手段１７ｈにより、図５及び図６に示すように、例えば、心臓画像Ｇ１の各特異点Ａ１〜Ａ４間の離間距離（特異点間距離）ｍ１〜ｍ５と、心臓モデルＭ１の各特異点Ｂ１〜Ｂ４間の離間距離（特異点間距離）ｎ１〜ｎ５とが一致するように（ｍ１＝ｎ１、ｍ２＝ｎ２、ｍ３＝ｎ３、ｍ４＝ｎ４、ｍ５＝ｎ５）、心臓画像Ｇ１の各特異点Ａ１〜Ａ４又は心臓モデルＭ１の各特異点Ｂ１〜Ｂ４が動かされる。このときの各伸縮量が各セルに線形に分配され（線形分配）、すなわち各セルの伸縮量が線形に調整され、各セルの移動量（特異点の移動情報）が求められる。これらのセルの移動量は、合成用データＤ１として記憶部１５に格納される。
【００３２】
画像及びモデル間の輪郭フィッティング手段１７ｉは、心臓画像Ｇ１及び心臓モデルＭ１間において、内壁や外壁等の輪郭をフィッティングする手段である。このとき、心臓画像Ｇ１の輪郭が動かされても良いし、心臓モデルＭ１の輪郭が動かされても良い。これは、操作部１４に対する操作者の入力操作に応じて予め設定されている。このフィッティング手段により、輪郭の移動情報（例えば、回転量や伸縮量等を示す変換係数）が求められる。この輪郭の移動情報は、合成用データＤ１として記憶部１５に格納される。
【００３３】
この画像及びモデル間の輪郭フィッティング手段１７ｉにより、図５及び図６に示すように、例えば、各特異点間距離ｍ１〜ｍ５、ｎ１〜ｎ５を固定しながら（各特異点Ａ１〜Ａ４、Ｂ１〜Ｂ４を固定しながら）、心臓画像Ｇ１の内壁における複数の中点（図５中の×印）から１つの中点ａ１が選択され、選択した中点ａ１に相関のある中点ｂ１が心臓モデルＭ１から選択される。次いで、特異点Ａ１−中点ａ１−特異点Ａ４の形状と特異点Ｂ１−中点ｂ１−特異点Ｂ４の形状とが一致するように心臓画像Ｇ１の中点ａ１又は心臓モデルＭ１の中点ｂ１が動かされる。その動かされた中点ａ１又は中点ｂ１に相当するセルも一緒に動かされる。このときの各伸縮量（セルが動いた分だけ）が各セルに線形に分配され（線形分配）、すなわち各セルの伸縮量が線形に調整され、各セルの移動量（輪郭の移動情報）が求められる。これらのセルの移動量は、合成用データＤ１として記憶部１５に格納される。このような処理が全ての中点について行われ、また、心臓画像Ｇ１の外壁でも同様に行われる。
【００３４】
画像及びモデルの位置合わせ手段１７ｊは、記憶部１５に格納された合成用データＤ１に基づいて、心臓画像Ｇ１と心臓モデルＭ１との位置合わせを行ってそれらを合成し、心臓画像Ｇ１及び心臓モデルＭ１の合成画像Ｇ２を生成する手段である。これにより、それらの合成画像Ｇ２が得られる。この合成画像Ｇ２は表示部１３に表示され、加えて、必要に応じて記憶部１５に保存される。
【００３５】
この画像及びモデルの位置合わせ手段１７ｊにより、合成用データＤ１、すなわち各セルの移動量を用いて、図４に示すように、心臓画像Ｇ１及び心臓モデルＭ１の位置合わせ（大きさ及び形状等）が合わせられ、心臓画像Ｇ１及び心臓モデルＭ１が合成される。これにより、心臓画像Ｇ１及び心臓モデルＭ１の合成画像Ｇ２が得られ、その合成画像Ｇ２は表示部１３に表示される。このとき、合成画像Ｇ２上には、壁運動異常部分Ｒ１や壁厚異常部分Ｒ２もそれぞれ色を変えて表示される。
【００３６】
画像及びモデル間の輪郭フィッティング手段１７ｉと、画像及びモデルの位置合わせ手段１７ｊとは、一般的に次のように構成することができる。フィッティング前において、心臓画像Ｇ１の各特異点Ａ１〜Ａ４は心臓画像上の座標系で表されており、心臓モデルＭ１の各特異点Ｂ１〜Ｂ４は心臓モデルの座標系で表されている。まず、これら心臓画像Ｇ１の特異点Ａ１〜Ａ４と心臓モデルＭ１の各特異点Ｂ１〜Ｂ４のそれぞれの座標系での座標値を用いて、両座標系の間の関係を求める。一般に、それらの関係はａ＝ｆ（ｂ）の様に表される。ａは心臓画像の座標値であり、ｂは心臓モデルの座標値である。ｆは両者の座標の写像を表す関数であり、代表的な例としては線形関数による線形写像が挙げられる。他の例としては多項式関数やｓｐｌｉｎｅ関数等が挙げられる。これらの後者の場合、変形を含むフィッティングを行うことができる。いずれの場合も関数ｆの具体的な形状は有限個のパラメータで表される。フィッティング処理では、心臓画像Ｇ１の各特異点Ａ１〜Ａ４と心臓モデルＭ１の各特異点Ｂ１〜Ｂ４との各々の座標値をａ＝ｆ（ｂ）に代入して得られる連立方程式を解くことにより、関数ｆのパラメータを決定する。得られたパラメータセットは心臓画像Ｇ１の座標系と心臓モデルＭ１の座標系との点の対応関係を具体的に表すものであり、これが決定された段階で、全ての点の対応関係が得られたことになる。心臓画像Ｇ１にて抽出された輪郭点ａ１、ａ２、・・・をａ＝ｆ（ｂ）に代入し、これらの心臓モデルＭ１内での対応する位置のｂ１、ｂ２、・・・を求めることもできる。その逆に、心臓モデルＭ１の輪郭点ｂ１、ｂ２、・・・をａ＝ｆ（ｂ）に代入し、これらの心臓モデルＭ１内での対応する位置のａ１、ａ２、・・・を求めることもできる。
【００３７】
このような画像処理部１７では、シミュレーションとリアルタイム心エコーとの心電同期技術や呼吸動除去技術、加えて、心疾患シミュレーションとＣＴ画像とをレジストレーションする技術、心疾患シミュレーションモデルとリアルタイム超音波画像の取得断面及び領域との位置合わせ技術、そのモデルから対応断面及び対応領域を切り出すＭＰＲ（Multi Planer Reformation）やＶＲ（Volume Reformation）技術（バーチャルソノグラフィー技術）等が用いられる。
【００３８】
次に、このような医用画像診断装置１、特に、医用画像処理装置５の合成表示処理について図７及び図８を参照して説明する。医用画像処理装置５の制御部１１が、画像処理部１７により合成処理を行い、合成処理後の合成画像Ｇ２を表示する表示処理を実行する。ここでは、一例として、心臓画像Ｇ１を心臓モデルＭ１に合わす合成処理を行う。
【００３９】
図７に示すように、制御部１１は、画像処理部１７により、同じ位相時間の心臓画像Ｇ１及び心臓モデルＭ１を選択し（ステップＳ１）、選択した心臓画像Ｇ１の各特異点Ａ１〜Ａ４（図５参照）を抽出し（ステップＳ２）、選択した心臓モデルＭ１の各特異点Ｂ１〜Ｂ４（図６参照）を抽出する（ステップＳ３）。
【００４０】
次いで、制御部１１は、画像処理部１７により、心臓画像Ｇ１の各特異点Ａ１〜Ａ４間の離間距離（特異点間距離）ｍ１〜ｍ５を算出し（ステップＳ４）、心臓モデルＭ１の各特異点Ｂ１〜Ｂ４間の離間距離（特異点間距離）ｎ１〜ｎ５を算出する（ステップＳ５）。
【００４１】
次に、制御部１１は、画像処理部１７により、心臓画像Ｇ１及び心臓モデルＭ１の各特異点間距離ｍ１〜ｍ５、ｎ１〜ｎ５が同じになるように、心臓画像Ｇ１において原点となる特異点Ａ２以外の特異点Ａ１、Ａ３、Ａ４の１つを移動させ（ステップＳ６）、このときの伸縮方向と伸縮量に合わせて、その伸縮量を各セルに線形に分配し、すなわち各セルの伸縮量を線形に調整する（ステップＳ７）。その後、制御部１１は、各セルの移動量及び変形量（特異点Ａ１〜Ａ４、Ｂ１〜Ｂ４の移動情報）を合成用データＤ１として記憶部１５に保存する（ステップＳ８）。
【００４２】
さらに、制御部１１は、特異点Ａ１〜Ａ４、Ｂ１〜Ｂ４の位置が完全に一致しているか否かを判断し（ステップＳ９）、それらの位置が完全に一致するまでステップＳ６〜Ｓ８を繰り返す（ステップＳ９のＮＯ）。それらの位置が、設定された精度（例えば、１ｍｍ）の範囲内で完全に一致していると判断した場合には（ステップＳ９のＹＥＳ）、図８に示すように、心臓画像Ｇ１の内壁を選択し（ステップＳ１０）、選択した内壁の輪郭を抽出、すなわち内壁をトレースし（ステップＳ１１）、トレースした線分中に複数の中点（図５及び図６参照）を設定する（ステップＳ１２）。このときの中点の設定は、等間隔でもよく、偏っていてもよい。なお、前述の精度の設定は、操作部１４に対する操作者の入力操作に応じて予め設定されている。
【００４３】
次に、制御部１１は、画像処理部１７により、複数の中点から１つの中点ａ１を選択し（ステップＳ１３）、各特異点Ａ１〜Ａ４を固定させた状態で、その中点ａ１を移動させる（ステップＳ１４）。このとき、例えば、制御部１１は、心臓画像Ｇ１での異なる２つの特異点と中点とからなる形状（特異点Ａ１−中点ａ１−特異点Ａ４の形状）と、それらの特異点及び中点に対応する心臓モデルＭ１での異なる２つの特異点と中点とからなる形状（特異点Ｂ１−中点ｂ１−特異点Ｂ４の形状）とが一致するように中点ａ１を動かす。このとき、中点ａ１に相当するセルもしくは要素（モデルを構成する要素）も一緒に動かされる。
【００４４】
次いで、制御部１１は、このときの伸縮量を各セルに線形に分配し、すなわち各セルの伸縮量を線形に調整し、各セルの移動量（内壁の輪郭の移動情報）を合成用データＤ１として記憶部１５に保存する（ステップＳ１５）。各特異点Ａ１〜Ａ４と中点との距離及び方向が必要な精度内に達したか否かを判断し（ステップＳ１６）、全ての中点においてその距離及び方向が必要な精度内に達するまで、ステップＳ１３〜Ｓ１５を繰り返す（ステップＳ１６のＮＯ）。
【００４５】
その距離及び方向が必要な精度内に達したと判断した場合には（ステップＳ１６のＹＥＳ）、心臓画像Ｇ１の外壁を選択し（ステップＳ１７）、選択した外壁の輪郭を抽出、すなわち外壁をトレースし（ステップＳ１８）、トレースした線分中に複数の中点を設定する（ステップＳ１９）。このときの中点の設定は、等間隔でもよく、偏っていてもよい。
【００４６】
次に、制御部１１は、画像処理部１７により、複数の中点から１つの中点を選択し（ステップＳ２０）、特異点Ａ１〜Ａ４を固定させた状態で、その中点を移動させる（ステップＳ２１）。このときも、前述同様、例えば、制御部１１は、心臓画像Ｇ１での異なる２つの特異点と中点とからなる形状と、それらの特異点及び中点に対応する心臓モデルＭ１での異なる２つの特異点と中点とからなる形状とが一致するように中点を動かす。
【００４７】
次いで、制御部１１は、このときの伸縮量を各セルに線形に分配し、すなわち各セルの伸縮量を線形に調整し、各セルの移動量（外壁の輪郭の移動情報）を合成用データＤ１として記憶部１５に保存する（ステップＳ２２）。各特異点と中点との距離及び方向が必要な精度内に達したか否かを判断し（ステップＳ２３）、全ての中点においてその距離及び方向が必要な精度内に達するまで、ステップＳ２０〜Ｓ２２を繰り返す（ステップＳ２３のＮＯ）。
【００４８】
その距離及び方向が必要な精度内に達したと判断した場合には（ステップＳ２３のＹＥＳ）、記憶部１５に保存されている合成用データＤ１、すなわち各セルの移動量を用いて、心臓画像Ｇ１を変形し、心臓モデルＭ１に一致させる（ステップＳ２４）。これにより、心臓画像Ｇ１と心臓モデルＭ１との合成画像Ｇ２が得られる（図４参照）。制御部１１は、その合成画像Ｇ２を表示部１３に表示し、さらに、必要に応じて、その合成画像Ｇ２を記憶部１５に保存する。このとき、制御部１１は、合成画像Ｇ２上に壁運動異常部分Ｒ１や壁厚異常部分Ｒ２をそれぞれの色を変えて表示する。
【００４９】
このような合成表示処理により、合成画像Ｇ２が表示部１３に表示され、医師等の診断者により視認される。このとき、心臓画像Ｇ１と心臓モデルＭ１とが１つの合成画像Ｇ２として表示され、被検体の心臓の形態情報とその機能情報との関係認識が容易になるので、診断者は被検体の心臓の状況（状態）を的確に把握することができる。さらに、心臓モデルＭ１上でシミュレーションするための疾患部位の情報が心臓画像Ｇ１上に重畳して表示されるので、注目部位である心臓の疾患に対する新たな質的診断情報を診断者に与えることができる。特に、診断者の経験や熟練度等に依存することなく、患者の心臓の状況（形態及び機能等）を的確に把握することが可能になるので、正確な診断を行うことができる。
【００５０】
さらに、合成画像Ｇ２に壁運動異常部分Ｒ１や壁厚異常部分Ｒ２等の異常部分を表示することによって、その異常部分が合成画像Ｇ２に表示され、心臓の状況に加え、その異常部分を即座にかつ正確に把握することができる。すなわち、壁運動異常がない状態のＣＴ画像でＡ氏（個人）の心臓モデルＭ１を作成し、その心臓モデルＭ１を用いた心疾患シミュレーションによって、壁運動異常等の異常を合成画像Ｇ２上に再現することが可能になるので、心臓の疾患に対する新たな質的診断情報が医師等に与えられ、医師等の診断者は、被検体の心臓の形態情報とその機能情報との関係認識を容易に行うことができる。さらに、ＣＴ画像や超音波画像等の心臓画像Ｇ１と、心臓モデルＭ１の心疾患シミュレーションとの様々な比較検討及び比較診断を容易に行うことが可能になるので、患者の心臓の状況（形態及び機能等）を的確に把握することができ、その結果、正確な診断を行うことができる。
【００５１】
以上説明したように、本発明の実施の形態によれば、医用画像撮影装置２により撮影された画像であって被検体の注目部位を示す医用画像（例えば心臓画像Ｇ１）と、注目部位の疾患をシミュレーションするためのモデルであって注目部位を示す部位モデル（例えば心臓モデルＭ１）とを合成し、医用画像及び部位モデルの合成画像Ｇ２を生成することによって、医用画像と部位モデルとが１つの合成画像Ｇ２として表示され、被検体の注目部位の形態情報とその機能情報との関係認識が容易になるので、医師等の診断者は被検体の心臓の状況（状態）を的確に把握することができる。また、部位モデル上でシミュレーションするための疾患部位の情報を医用画像（診断画像）上に重畳し、合成画像Ｇ２として表示することが可能になるので、注目部位の疾患に対する新たな質的診断情報を与えることができる。その結果、医師等の診断者は被検体の注目部位の状況を的確に把握することができる。
【００５２】
特に、注目部位が心臓である場合には、ＣＴ画像や超音波画像等の心臓画像Ｇ１と心疾患シミュレーションの心臓モデルＭ１との複合的な診断情報（画像情報）を比較表示、すなわち重畳表示することが可能になるので、心疾患の総合診断環境を提供することができる。すなわち、虚血性心疾患や狭心症等に関する質的診断情報、例えば、心筋細胞のイオン電流の活性度や心筋細胞の機械的な伸縮度（動き）等を提示することができる。これにより、診断者は、心筋の梗塞状態や虚血状態を推測するための新たな質的診断情報を得ることができる。加えて、超音波画像だけではカバーすることができない広い視野領域で心機能を観察することができる。
【００５３】
また、医用画像（例えば心臓画像Ｇ１）と部位モデル（例えば心臓モデルＭ１）との合成を行う場合には、医用画像の複数の特異点Ａ１〜Ａ４及び部位モデルの複数の特異点Ｂ１〜Ｂ４を抽出し、抽出した医用画像の複数の特異点Ａ１〜Ａ４及び部位モデルの複数の特異点Ｂ１〜Ｂ４に基づいて医用画像の注目部位と部位モデルの注目部位との位置合わせを行うことによって、精度良く医用画像及び部位モデルの合成画像Ｇ２を生成することができる。
【００５４】
さらに、位置合わせを行う場合には、抽出した医用画像（例えば心臓画像Ｇ１）の複数の特異点Ａ１〜Ａ４及び部位モデル（例えば心臓モデルＭ１）の複数の特異点Ｂ１〜Ｂ４に基づいて、医用画像の注目部位と部位モデルの注目部位との大きさを一致させ、抽出した医用画像の複数の特異点Ａ１〜Ａ４及び部位モデルの複数の特異点Ｂ１〜Ｂ４に基づいて、医用画像の注目部位と部位モデルの注目部位との形状を一致させることによって、より精度良く医用画像及び部位モデルの合成画像Ｇ２を生成することができる。
【００５５】
加えて、大きさを一致させる場合には、抽出した医用画像（例えば心臓画像Ｇ１）の複数の特異点Ａ１〜Ａ４から各々の特異点間の離間距離ｍ１〜ｍ５及び抽出した部位モデル（例えば心臓モデルＭ１）の複数の特異点Ｂ１〜Ｂ４から各々の特異点間の離間距離ｎ１〜ｎ５を求め、求めた医用画像の各々の離間距離ｍ１〜ｍ５と求めた部位モデルの各々の離間距離ｎ１〜ｎ５とがそれぞれ同じになるように医用画像の注目部位と部位モデルの注目部位との大きさを一致させることから、容易な処理により医用画像及び部位モデルの合成画像Ｇ２を生成することが可能になり、さらに、医用画像及び部位モデルの大きさを精度良く一致させることが可能になるので、処理時間が長くなることを防止しつつ精度が高い合成画像Ｇ２を得ることができる。
【００５６】
また、形状を一致させる場合には、医用画像の注目部位（例えば心臓画像Ｇ１）の輪郭及び部位モデル（例えば心臓モデルＭ１）の注目部位の輪郭を抽出し、抽出した医用画像の複数の特異点及び抽出した部位モデルの複数の特異点に基づいて、抽出した医用画像の注目部位の輪郭と抽出した部位モデルの注目部位の輪郭とが同じになるように医用画像の注目部位と部位モデルの注目部位との形状を一致させることから、容易な処理により医用画像及び部位モデルの合成画像Ｇ２を生成することが可能になり、さらに、医用画像及び部位モデルの形状を精度良く一致させることが可能になるので、処理時間が長くなることを防止しつつ精度が高い合成画像Ｇ２を得ることができる。
【００５７】
さらに、生成した合成画像Ｇ２に注目部位の異常部分（例えば、壁運動異常部分Ｒ１や壁厚異常部分Ｒ２）を示すことによって、その異常部分が合成画像Ｇ２に表示されるので、注目部位の疾患に対する新たな質的診断情報を与えることが可能になるので、診断者は心臓の状況に加え、その異常部分を即座にかつ正確に把握することができ、さらに、心臓の状況（形態及び機能等）に対する異常部分の位置関係等も容易に把握することができる。
【００５８】
（他の実施の形態）
なお、本発明は、前述の実施の形態に限るものではなく、その要旨を逸脱しない範囲において種々変更可能である。
【００５９】
例えば、前述の実施の形態においては、医用画像情報（心臓画像Ｇ１）やモデル情報（心臓モデルＭ１）を取得する場合、時相として、拡張期の特定の時相（拡張末期）を用いているが、これに限るものではなく、他の時相を用いてもよく、さらに、時相毎に医用画像情報やモデル情報を取得して合成画像Ｇ２を生成するようにしてもよい。この場合には、生成した複数の合成画像Ｇ２を時系列に並べて表示部１３に連続表示させることによって、心臓の状況の時系列的な変化を直感的に認識することが可能になるので、心臓の状況を容易に把握することができる。
【００６０】
また、前述の実施の形態においては、Ｘ線断層撮影装置（Ｘ線ＣＴ装置）や超音波診断装置等を用いて医用画像撮影装置２を構成しているが、これに限るものではなく、例えば、他の撮影装置を用いるようにしてもよく、さらに、Ｘ線断層撮影装置や超音波診断装置等の医用画像撮影装置２に医用画像処理装置５を組み込み、医用画像診断装置１を構成するようにしてもよい。
【図面の簡単な説明】
【００６１】
【図１】本発明の実施の一形態に係る医用画像診断装置の概略構成を示すブロック図である。
【図２】図１に示す医用画像診断装置が備える医用画像処理装置の概略構成を示すブロック図である。
【図３】図２に示す医用画像処理装置が備える画像処理部の概略構成を示すブロック図である。
【図４】図３に示す画像処理部が行う画像合成を説明するための説明図である。
【図５】図３に示す画像処理部が行う画像合成を説明するために心臓画像を示す模式図である。
【図６】図３に示す画像処理部が行う画像合成を説明するために心臓モデルを示す模式図である。
【図７】図１に示す医用画像診断装置、特に、図２に示す医用画像処理装置が行う合成処理の流れを示すフローチャートである。
【図８】図１に示す医用画像診断装置、特に、図２に示す医用画像処理装置が行う合成処理の流れを示すフローチャートである。
【符号の説明】
【００６２】
１          医用画像診断装置
２          医用画像撮影装置
５          医用画像処理装置
Ａ１〜Ａ４  特異点
Ｂ１〜Ｂ４  特異点
Ｇ１        医用画像（心臓画像）
Ｇ２        合成画像
Ｍ１        部位モデル（心臓モデル）
ｍ１〜ｍ５  離間距離（特異点離間距離）
ｎ１〜ｎ５  離間距離（特異点離間距離）
【図１】
0005173303.tif 000002
【図２】
0005173303.tif 000003
【図３】
0005173303.tif 000004
【図４】
0005173303.tif 000005
【図５】
0005173303.tif 000006
【図６】
0005173303.tif 000007
【図７】
0005173303.tif 000008
【図８】
0005173303.tif 000009