(57)【特許請求の範囲】
【請求項１】
医用画像診断装置によって取得された所定の被検体に関する医用情報を用いて、当該被検体に関する複数のパラメータを抽出すると共に、予め取得された健常者及び非健常者に関する医用情報を用いて、前記健常者に関する複数の第１のパラメータと、前記非健常者における複数の疾患のそれぞれに関する前記複数の第２のパラメータとを抽出する抽出手段と、
前記複数の第１のパラメータを項目因子としてＭＴシステムを実行することで生成される前記健常者の状態空間における第１の物差しを計算し、前記非健常者の複数の第２のパラメータを項目因子としてＭＴシステムを実行することで生成される前記疾患症例毎の状態空間における複数の第２の物差しを計算し、前記被検体に関する複数のパラメータを項目因子としてＭＴシステムを実行することで、前記被検体と前記健常者及び前記複数の疾患の状態空間との間のそれぞれの距離を計算する計算手段と、
前記第１の前記物差しおよび前記疾患症例毎の前記複数の第２の物差しと前記被検体の前記複数の疾患症例毎の状態空間からの前記それぞれの距離とを用いて、前記被検体の健康状態が正常であるか異常であるいかを判定し、及び健康状態が異常である場合の疾患種別を判定する判定手段と、
前記判定の結果に基づいて、診断支援情報を生成する支援情報生成手段と、
前記診断支援情報を所定の形態で表示する表示手段と、
を具備することを特徴とする自動診断支援装置。
【請求項２】
前記医用画像診断装置は、超音波診断装置であることを特徴とする請求項１記載の自動診断支援装置。
【請求項３】
前記所定のパラメータは、Ｂモード画像に関するパターン認識データ、ドプラ波形による波形パターン認識データ、ドプラ自動計測による診断指標、心電図に基づく診断指標のうちの少なくとも一つを含むことを特徴とする請求項１又は２記載の自動診断支援装置。
【請求項４】
前記計算手段は、前記各疾患に関する物差し及び前記被検体の前記各疾患の非健常者の状態空間からの距離を正規化し、
前記判定手段は、
正規化された前記各疾患症例に関する物差し及び正規化された前記被検体の前記各疾患の非健常者の状態空間からの距離を用いて、疾患判別の相対的距離を計算し、
前記疾患判別の相対的距離に基づいて、前記被検体の疾患種別を判定すること、
を特徴とする請求項１乃至３のうちいずれか一項記載の自動診断支援装置。
【請求項５】
前記支援情報生成手段は、前記診断支援情報として、各疾患を変数するレーダーチャートを生成することを特徴とする請求項１乃至４のうちいずれか一項記載の自動診断支援装置。
【請求項６】
前記支援情報生成手段は、各疾患につき外周に行くほど当該疾患の可能性が高くなり、且つ正常と推定される範囲を参照表示する前記レーダーチャートを生成することを特徴とする請求項５記載の自動診断支援装置。
【請求項７】
前記支援情報生成手段は、前記診断支援情報として、前記各疾患の可能性を示すための定量値又はヒストグラムを生成することを特徴とする請求項１乃至６のうちいずれか一項記載の自動診断支援装置。
【請求項８】
前記医用画像診断装置による前記医用情報の取得の完了に連動して、当該装置の動作タイミングを制御する制御手段をさらに具備することを特徴とする請求項１乃至７のうちいずれか一項記載の自動診断支援装置。
【請求項９】
所定の被検体を超音波走査し、当該被検体に関する超音波データを取得する超音波データ取得手段と、
前記超音波データを含む前記被検体に関する医用情報を用いて、当該被検体に関する複数のパラメータを抽出すると共に、予め取得された健常者及び非健常者に関する医用情報を用いて、前記健常者に関する複数の第１のパラメータと、前記非健常者における複数の疾患のそれぞれに関する前記複数の第２のパラメータとを抽出する抽出手段と、
前記複数の第１のパラメータを項目因子としてＭＴシステムを実行することで生成される前記健常者の状態空間における第１の物差しを計算し、前記非健常者の複数の第２のパラメータを項目因子としてＭＴシステムを実行することで生成される前記疾患症例毎の状態空間における複数の第２の物差しを計算し、前記被検体に関する複数のパラメータを項目因子としてＭＴシステムを実行することで、前記被検体と前記健常者及び前記複数の疾患の状態空間との間のそれぞれの距離を計算する計算手段と、
前記第１の前記物差しおよび前記疾患症例毎の前記複数の第２の物差しと前記被検体の前記複数の疾患症例毎の状態空間からの前記それぞれの距離とを用いて、前記被検体の健康状態が正常であるか異常であるいかを判定し、及び健康状態が異常である場合の疾患種別を判定する判定手段と、
前記判定の結果に基づいて、診断支援情報を生成する支援情報生成手段と、
前記診断支援情報を所定の形態で表示する表示手段と、
を具備することを特徴とする超音波診断装置。
【請求項１０】
コンピュータに、
医用画像診断装置によって取得された所定の被検体に関する医用情報を用いて、当該被検体に関する複数のパラメータを抽出すると共に、予め取得された健常者及び非健常者に関する医用情報を用いて、前記健常者に関する複数の第１のパラメータと、前記非健常者における複数の疾患のそれぞれに関する前記複数の第２のパラメータとを抽出する抽出機能と、
前記複数の第１のパラメータを項目因子としてＭＴシステムを実行することで生成される前記健常者の状態空間における第１の物差しを計算し、前記非健常者の複数の第２のパラメータを項目因子としてＭＴシステムを実行することで生成される前記疾患症例毎の状態空間における複数の第２の物差しを計算し、前記被検体に関する複数のパラメータを項目因子としてＭＴシステムを実行することで、前記被検体と前記健常者及び前記複数の疾患の状態空間との間のそれぞれの距離を計算する計算機能と、
前記第１の前記物差しおよび前記疾患症例毎の前記複数の第２の物差しと前記被検体の前記複数の疾患症例毎の状態空間からの前記それぞれの距離とを用いて、前記被検体の健康状態が正常であるか異常であるいかを判定し、及び健康状態が異常である場合の疾患種別を判定する判定機能と、
前記判定の結果に基づいて、診断支援情報を生成する支援情報生成機能と、
前記診断支援情報を所定の形態で表示する表示機能と、
を実現させることを特徴とする自動診断支援プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像を用いたコンピュータ支援診断（computer-aided diagnosis： CAD）に関するもので、特に、異なるアルゴリズムに従う複数の読影結果を用いるコンピュータ支援診断に関する。
【背景技術】
【０００２】
コンピュータ支援診断とは、診断における最終判断は医師が行うことを前提としつつ、医師や技師による画像診断の負荷を低減させるため、画像パターン認識技術等を用いて、入力したデジタル画像から検出した病巣候補の位置や、コンピュータを用いて解析した定量的な結果を画像上にわかりやすく提示するものである。この技術は、例えば胸部X線写真や集団検診用CT画像による肺癌の検診、乳癌マンモグラフィの自動スクリーニングシステム、胃癌の検出システム、超音波画像診断における循環器ドプラ診断等に用いられる。
【０００３】
この様なコンピュータ支援診断を用いた自動診断支援装置では、近年、ＭＴ(Mahalanobis（マハラノビス）・Taguchi（タグチ）)システムと呼ばれる手法が用いられている。ＭＴシステムとは、品質工学において多変の項目（多変量）を総合して新しい尺度を構成する手法であり、その類型として、ＭＴ（Mahalanobis Taguchi）法、Ｔ（Taguchi）法、ＭＴＡ（Mahalanobis Taguchi Ajoint）法、ＴＳ（Taguchi schmitt）法等がある。このＭＴシステムが被検体の正常／異常の識別判断に応用され始めて約５年経ち、近年では、特殊な応用として肝疾患種類の判別にＳＮ比の要因効果パターンによる識別が試みられている。
【０００４】
なお、本願発明に関連する公知文献としては、次の様なものがある。
【先行技術文献】
【特許文献】
【０００５】
【非特許文献１】“心臓超音波検査法トレーニングプロトコールの評価  ―正確性と再現性について―  ” J Med Ultrasonics Vol.29 No.6 p537-544 (2002)
【非特許文献２】“我国における心エコー図のルーチン検査の現状（アンケート調査  ―第２報―  ” Jpn J Med Ultrasonics Vol.32 No.3 p329-337 (2005)
【非特許文献３】“心エコー図検査による加齢の心機能に対する影響の検討” J Med Ultrasonics Vol.29 No.2 p145-151 (2002)
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、従来の自動診断支援装置においては次のような問題がある。
【０００７】
すなわち、画像パターン認識技術としてはニューラルネットワーク、主成分分析法、多変量解析(回帰分析/因子分析)が主流であり、交互作用、モデルパターンの抽出、パラメータの変動範囲、主成分の選別方式などに統一性がない。従って、現在においても汎用的な解析手法が確立されておらず、信憑性／再現性に乏しいものとなっている。
【０００８】
また、例えばＭＴシステムを用いた肝疾患種別の判別方法についても、健常者の状態空間をもとにした異常パターンでの、個々の症例の距離が離れる要因パターンの分析を行っている。このため、例えば肝臓では好適な自動判定ができても、他の疾患判別に役立つかどうか(汎用性があるかどうか)不明である。
【０００９】
本発明は、上記事情を鑑みてなされたもので、ＭＴシステムを用いて疾患症例毎の尺度を作成し、当該尺度と自動計測によって取得された項目とを用いて汎用性の高い自動診断支援を実行可能な自動診断支援装置、超音波診断装置及び自動診断支援プログラムを提供することを目的としている。
【課題を解決するための手段】
【００１０】
本発明は、上記目的を達成するため、次のような手段を講じている。
【００１１】
請求項１に記載の発明は、医用画像診断装置によって取得された所定の被検体に関する医用情報を用いて、当該被検体に関する複数のパラメータを抽出すると共に、予め取得された健常者及び非健常者に関する医用情報を用いて、前記健常者に関する複数の第１のパラメータと、前記非健常者における複数の疾患のそれぞれに関する前記複数の第２のパラメータとを抽出する抽出手段と、前記複数の第１のパラメータを項目因子としてＭＴシステムを実行することで生成される前記健常者の状態空間における第１の物差しを計算し、前記非健常者の複数の第２のパラメータを項目因子としてＭＴシステムを実行することで生成される前記疾患症例毎の状態空間における複数の第２の物差しを計算し、前記被検体に関する複数のパラメータを項目因子としてＭＴシステムを実行することで、前記被検体と前記健常者及び前記複数の疾患の状態空間との間のそれぞれの距離を計算する計算手段と、前記第１の前記物差しおよび前記疾患症例毎の前記複数の第２の物差しと前記被検体の前記複数の疾患症例毎の状態空間からの前記それぞれの距離とを用いて、前記被検体の健康状態が正常であるか異常であるいかを判定し、及び健康状態が異常である場合の疾患種別を判定する判定手段と、前記判定の結果に基づいて、診断支援情報を生成する支援情報生成手段と、前記診断支援情報を所定の形態で表示する表示手段と、を具備することを特徴とする自動診断支援装置である。
請求項９に記載の発明は、所定の被検体を超音波走査し、当該被検体に関する超音波データを取得する超音波データ取得手段と、前記超音波データを含む前記被検体に関する医用情報を用いて、当該被検体に関する複数のパラメータを抽出すると共に、予め取得された健常者及び非健常者に関する医用情報を用いて、前記健常者に関する複数の第１のパラメータと、前記非健常者における複数の疾患のそれぞれに関する前記複数の第２のパラメータとを抽出する抽出手段と、前記複数の第１のパラメータを項目因子としてＭＴシステムを実行することで生成される前記健常者の状態空間における第１の物差しを計算し、前記非健常者の複数の第２のパラメータを項目因子としてＭＴシステムを実行することで生成される前記疾患症例毎の状態空間における複数の第２の物差しを計算し、前記被検体に関する複数のパラメータを項目因子としてＭＴシステムを実行することで、前記被検体と前記健常者及び前記複数の疾患の状態空間との間のそれぞれの距離を計算する計算手段と、前記第１の前記物差しおよび前記疾患症例毎の前記複数の第２の物差しと前記被検体の前記複数の疾患症例毎の状態空間からの前記それぞれの距離とを用いて、前記被検体の健康状態が正常であるか異常であるいかを判定し、及び健康状態が異常である場合の疾患種別を判定する判定手段と、前記判定の結果に基づいて、診断支援情報を生成する支援情報生成手段と、前記診断支援情報を所定の形態で表示する表示手段と、を具備することを特徴とする超音波診断装置である。
請求項１０に記載の発明は、コンピュータに、医用画像診断装置によって取得された所定の被検体に関する医用情報を用いて、当該被検体に関する複数のパラメータを抽出すると共に、予め取得された健常者及び非健常者に関する医用情報を用いて、前記健常者に関する複数の第１のパラメータと、前記非健常者における複数の疾患のそれぞれに関する前記複数の第２のパラメータとを抽出する抽出機能と、前記複数の第１のパラメータを項目因子としてＭＴシステムを実行することで生成される前記健常者の状態空間における第１の物差しを計算し、前記非健常者の複数の第２のパラメータを項目因子としてＭＴシステムを実行することで生成される前記疾患症例毎の状態空間における複数の第２の物差しを計算し、前記被検体に関する複数のパラメータを項目因子としてＭＴシステムを実行することで、前記被検体と前記健常者及び前記複数の疾患の状態空間との間のそれぞれの距離を計算する計算機能と、前記第１の前記物差しおよび前記疾患症例毎の前記複数の第２の物差しと前記被検体の前記複数の疾患症例毎の状態空間からの前記それぞれの距離とを用いて、前記被検体の健康状態が正常であるか異常であるいかを判定し、及び健康状態が異常である場合の疾患種別を判定する判定機能と、前記判定の結果に基づいて、診断支援情報を生成する支援情報生成機能と、前記診断支援情報を所定の形態で表示する表示機能と、を実現させることを特徴とする自動診断支援プログラムである。
【発明の効果】
【００１４】
以上本発明によれば、ＭＴシステムを用いて疾患症例毎の尺度を作成し、当該尺度と自動計測によって取得された項目とを用いて汎用性の高い自動診断支援を実行可能な自動診断支援装置、超音波診断装置及び自動診断支援プログラムを実現することができる。
【図面の簡単な説明】
【００１５】
【図１】図１は、本実施形態に係る自動診断支援システム１のブロック構成図を示している。
【図２】図２は、マハラノビス距離の一般的な計算の流れを示したフローチャートである。
【図３】図３は、マハラノビス距離を用いた判定処理の概念を説明するための図である。
【図４】図４は、ＭＴシステムのＳＮ比を用いた原因診断を説明するための図である。
【図５】図５は、本自動診断支援機能に従う処理（自動診断支援処理）の流れを示したフローチャートである。
【図６】図６（ａ）、（ｂ）、（ｃ）は、ステップＳ４において解析ユニット３４が実行する判定処理の概念を説明するための図である。
【図７】図７は、相対的距離を基準とする当該患者甲についての各状態空間（単位空間）からの近さを説明するための図である。
【図８】図８は、ステップＳ５において生成される診断支援情報としてのヒストグラムの一例を示した図である。
【図９】図９（ａ）、（ｂ）は、良好な結果が得られたＤＣＭのヒストグラムと、判別困難だったＡＲのヒストグラムを示した図である。
【図１０】図１０は、本ステップＳ５において生成される診断支援情報としてのレーダーチャートの他の例を示した図である。
【図１１】図１１は、本ステップＳ５において生成される診断支援情報としてのレーダーチャートの他の例を示した図である。
【図１２】図１２は、健常者と非健常者(全疾患)のマハラノビス距離（ＭＤ値）ヒストグラムを示した図である。
【図１３】図１３は、疾患別の健常者との判別を評価するための各疾患毎（ＤＣＭ、ＭＲ、ＡＳ、ＡＲ）のヒストグラムを示した図である。
【図１４】図１４は、本実施例において生成された、マハラノビス距離（ＭＤ値）に関する疾患毎のヒストグラムを示している。
【発明を実施するための形態】
【００１６】
以下、本発明の実施形態を図面に従って説明する。なお、以下の説明において、略同一の機能及び構成を有する構成要素については、同一符号を付し、重複説明は必要な場合にのみ行う。
【００１７】
図１は、本実施形態に係る自動診断支援システム１のブロック構成図を示している。同図に示すように、本自動診断支援システム１は、超音波診断装置２と、自動診断支援装置３とから構成されている。なお、本自動診断支援システム１を構成する超音波診断装置２と自動診断支援装置３とは、別体でも一体でも構わない。また、自動診断支援装置３を超音波診断装置２に内蔵させ、超音波診断装置３そのものが自動診断支援システム１となる構成であってもよい。超音波診断装置２、自動診断支援装置３のそれぞれの構成について、以下説明する。
【００１８】
（超音波診断装置）
超音波診断装置２は、超音波プローブ２０、送受信ユニット２１、Ｂモード処理ユニット２２、ドプラ処理ユニット２３、画像・波形生成ユニット２４、画像合成ユニット２５、表示ユニット２６を具備している。
【００１９】
超音波プローブ２０は、超音波送受信ユニット２１からの駆動信号に基づき超音波を発生し、被検体からの反射波を電気信号に変換する複数の圧電振動子、当該圧電振動子に設けられる整合層、当該圧電振動子から後方への超音波の伝播を防止するバッキング材等を有している。当該超音波プローブ２０から被検体Ｐに超音波が送信されると、当該送信超音波は、体内組織の音響インピーダンスの不連続面で次々と反射され、エコー信号として超音波プローブ２０に受信される。このエコー信号の振幅は、反射することになった反射することになった不連続面における音響インピーダンスの差に依存する。また、送信された超音波パルスが、移動している血流や組織で反射された場合のエコーは、ドプラ効果により移動体の超音波送信方向の速度成分を依存して、周波数偏移を受ける。
【００２０】
送受信ユニット２１は、図示しないトリガ発生回路、遅延回路およびパルサ回路等を有している。パルサ回路では、所定のレート周波数ｆｒ  Ｈｚ（周期；１／ｆｒ秒）で、送信超音波を形成するためのレートパルスが繰り返し発生される。また、遅延回路では、チャンネル毎に超音波をビーム状に集束し且つ送信指向性を決定するのに必要な遅延時間が、各レートパルスに与えられる。トリガ発生回路は、このレートパルスに基づくタイミングで、プローブ２０に駆動パルスを印加する。
【００２１】
また、送受信ユニット２１は、図示していないアンプ回路、Ａ／Ｄ変換器、加算器等を有している。アンプ回路では、プローブ２０を介して取り込まれたエコー信号をチャンネル毎に増幅する。Ａ／Ｄ変換器では、増幅されたエコー信号に対し受信指向性を決定するのに必要な遅延時間を与え、その後加算器において加算処理を行う。この加算により、エコー信号の受信指向性に応じた方向からの反射成分が強調され、受信指向性と送信指向性とにより超音波送受信の総合的なビームが形成される。
【００２２】
Ｂモード処理ユニット２２は、送受信ユニット２１からエコー信号を受け取り、対数増幅、包絡線検波処理などを施し、信号強度が輝度の明るさで表現されるＢモード画像データを生成する。
【００２３】
ドプラ処理ユニット２３は、送受信ユニット２１から受け取ったエコー信号から速度情報を周波数解析し、ドプラ効果による血流や組織、造影剤エコー成分を抽出し、平均速度、分散、パワー等の血流情報を多点について求める。得られた血流情報は画像生成ユニット２５に送られ、平均速度画像、分散画像、パワー画像、これらの組み合わせ画像としてモニター１４にカラー表示される。
【００２４】
画像・波形生成ユニット２４は、超音波スキャンの走査線信号列を、テレビなどに代表される一般的なビデオフォーマットの走査線信号列に変換し、表示画像としての超音波診断画像やドプラ波形を生成する。また、画像・波形生成ユニット２４は、画像データを格納する記憶メモリを搭載しており、３次元画像の再構成処理などを行うことが可能である。また、例えば診断の後に操作者が検査中に記録された画像を呼び出すことが可能となっている。なお、当該画像・波形生成ユニット２４に入る以前のデータは、「生データ」と呼ばれることがある。
【００２５】
画像合成ユニット２５は、画像・波形生成ユニット２４から受け取った画像を種々のパラメータの文字情報や目盛等と共に合成し、ビデオ信号としてモニター１４に出力する。３次元の再構成プログラムや本発明の画像処理プログラム等も、ここに格納されており、操作者の指示などにより、これらのプログラムは起動される。
【００２６】
表示ユニット２６は、画像合成ユニット２５からのビデオ信号に基づいて、生体内の形態学的情報や、血流情報を画像として表示する。
【００２７】
なお、本超音波診断装置２は、上記各構成要素のほか、例えば、オペレータからの各種指示、条件、関心領域（ＲＯＩ）の設定指示、種々の画質条件設定指示等を当該装置２にとりこむための各種スイッチ、ボタン、トラックボール、マウス、キーボード等からなる入力ユニットや、ネットワークを介して他の装置とデータを送受信するための送受信ユニット等を具備している。
【００２８】
（自動診断支援装置）
自動診断支援装置３は、データベースユニット３０、ＭＴシステム実行ユニット３２、解析ユニット３４、診断支援情報生成ユニット３６、表示ユニット３８を具備している。
【００２９】
データベースユニット３０は、健常者データベース３０１と非健常者データベース３０３とを有する。
【００３０】
健常者データベース３０１は、複数の健常者に関する複数の自動計測値、超音波診断装置によって取得されたＢモード画像データ、超音波診断装置によって取得されたドプラ計測データ、超音波診断装置以外の医用画像診断装置（例えば、Ｘ線コンピュータ断層撮影装置、磁気共鳴イメージング装置、Ｘ線診断装置、核医学診断装置等）によって取得された画像データ、ＥＣＧ波形等を記憶する。なお健常者とは、計測値が計測された時点において「心臓機能に関する疾患を患っていない」と医師により判断された人のことである。
【００３１】
非健常者データベース３０３は、各疾患について、複数の非健常者に関する複数の自動計測値、超音波診断装置によって取得されたＢモード画像データ、超音波診断装置によって取得されたドプラ計測データ、超音波診断装置以外の医用画像診断装置（例えば、Ｘ線コンピュータ断層撮影装置、磁気共鳴イメージング装置、Ｘ線診断装置、核医学診断装置等）によって取得された画像データ、ＥＣＧ波形等を記憶する。ここで、非健常者とは、「後述する自動診断支援機能による疾患の判別対象となっている何れかの疾患を患っている」と医師により判断された人のことである。従って、ＭＴシステムにおいては、非健常者は、健常者の状態空間に属さない。疾患の種類としては、例えば、心肥大、心筋梗塞、僧帽弁逆流、大動脈弁逆流、中隔欠損、貧脈・余脈等がある。しかしながら、疾患の種類はこれに限定する必要はなく、上記疾患以外の疾患であっても適用可能である。
【００３２】
ＭＴシステム実行ユニット３２は、健常者データベース３０１に記憶されている各種計測値（診断指標）、画像パターン、ドプラ波形からパラメータを抽出し、抽出したパラメータを項目因子として多変量の健常者の状態空間を生成し、マハラノビス距離（物差し）を計算する。また、ＭＴシステム実行ユニット３２は、非健常者データベース３０３に記憶されている疾患毎の各種計測値（診断指標）、画像パターン、ドプラ波形からパラメータを抽出し、抽出したパラメータを項目因子として、疾患毎に多変量の非健常者の状態空間を生成し、疾患毎のマハラノビス距離（物差し）を計算する。さらに、ＭＴシステム実行ユニット３２は、超音波診断装置２によって自動計測された所定患者に関する自動計測値（診断指標）、画像パターン、ドプラ波形からパラメータを抽出し、抽出したパラメータを項目因子として、当該患者に関する、各疾患に対応する各状態空間からのマハラノビス距離を計算する。
【００３３】
解析ユニット３４は、所定患者に関する各疾患の状態空間からの各マハラノビス距離と、各疾患の状態空間毎に計算された物差しとに基づいて、当該患者の健康状態が正常であるか異常であるか、及び異常である場合の疾患種別を判定する。
【００３４】
図２は、マハラノビス距離の一般的な計算の流れを示したフローチャートである。同図に示すように、まず、各計測項目ｊ（１≦ｊ≦ｎ）について、計測値ベクトルｘｉｊ＝（ｘ１ｊ,ｘ２ｊ,・・・，ｘｍｊ）の平均値と標準偏差σｊとを算出する（ステップＳａ）。次に、算出されたｘｉｊの平均値と分散値σｊとを用いて、式（１）に従って各計測値ベクトルｘｉｊを正規化し、正規化された計測値ベクトルＸｉｊ＝（Ｘ１ｊ,Ｘ２ｊ,・・・，Ｘｍｊ）を算出する（ステップＳｂ）。
【数１】
0005611546.tif 000002
【００３５】
なお、母集団のσは、マハラノビス距離の平均値が１となるように、式（２）に従って操作（正規化）される。また、上記ステップＳａ、ステップＳｂの処理は、全て（ｎ個）の計測項目について行なわれる。
【００３６】
次に、項目毎の正規化された計測値ベクトルＸｉｊを用いて、次の式（３）、式（４）に従って相関行列Ｒとその逆行列Ｒ−１が計算される（ステップＳｃ）。
【数２】
0005611546.tif 000003
【００３７】
次に、相関行列の逆行列Ｒ−１を用いて、次の式（５）、（６）、（７）、（８）に従って、マハラノビス距離（ＭＤ）の平方Ｄ２を算出する（ステップＳｄ）。
【数３】
0005611546.tif 000004
【００３８】
この様な計算により、例えば図３に示すように、健常者の状態空間と患者Ａ、患者Ｂとの各マハラノビス距離を計算し、所定の閾値と比較することにより、患者Ａ、患者Ｂのそれぞれが健常者であるか否かを判定することができる。
【００３９】
なお、ＳＮ比についても、次の式（９）〜（１４）を用いて計算することができる。
【数４】
0005611546.tif 000005
【００４０】
上記ＳＮ比を用いて、例えば図４に示すようなヒストグラムを生成し、原因診断を判定することも可能である。
【００４１】
診断支援情報生成ユニット３６は、解析ユニット３４において実行された解析処理の結果に基づいて、例えば各疾患を患っている危険性（類似度が高いものほど危険性が高い）を示すレーダーチャートＲＣ等を診断支援情報として生成する。
【００４２】
表示ユニット３８は、診断支援情報生成ユニット３６において生成された診断支援情報を、所定の形態で表示する。
【００４３】
（自動診断支援機能）
次に、本自動診断支援システム１が有する、自動診断支援機能について説明する。この機能は、超音波診断装置等の医用画像診断装置の自動計測によって得られた自動計測値（診断指標）、画像パターン、波形等を用いてパラメータを抽出し、当該パラメータをＭＴシステムの項目因子として疾患症例毎に物差しを作成し、所定の患者（すなわち、正常又は異常、或いは異常である場合の疾患種別が未知である患者（被検体））に関して、健康状態が正常であるか又は異常であるか、或いは異常である場合の疾患種別を判定し、その結果を診断支援情報として所定の形態で表示するものである。
【００４４】
図３は、本自動診断支援機能に従う処理（自動診断支援処理）の流れを示したフローチャートである。以下、各ステップにおいて実行される処理の内容について説明する。なお、本実施形態においては、自動診断支援装置３の動作は、例えば超音波診断装置２による医用情報の取得の完了に連動して開始するように制御されるものとする。
【００４５】
［自動計測・ドプラ波形・Ｂモード画像等の取得：ステップＳ１］
まず、超音波診断装置２等を用いて、所定の患者についての自動計測が実行される。この自動計測により、ＥＣＧ波形、ドプラ波形、Ｂモード画像等が取得される（ステップＳ１）。
【００４６】
［パラメータ抽出・項目因子の決定：ステップＳ２］
次に、ＭＴシステム実行ユニット３２は、ステップＳ１において取得されたＥＣＧ波形、ドプラ波形、Ｂモード画像、各種診断指標等を用いて、ＭＴシステムの項目因子とすべき各種パラメータを計算する。具体的には、ＭＴシステム実行ユニット３２は、取得された超音波画像のパターン認識データ（Ｂモード画像の二次元の相関関数等）、微分積分法等を用いた波形パターン認識データ、ＥＣＧ波形やドプラ計測に基づく診断指標を計算する（ステップＳ２）。
【００４７】
なお、この項目因子とすべきパラメータの計算は、ステップＳ１において取得された情報に限定する必要はなく、例えば、事前に他のモダリティによって取得された画像データ、過去の診断において取得されたデータ等を採用するようにしてもよい。
【００４８】
［疾患症例毎の物差し、及び当該患者に関するマハラノビス距離の計算：ステップＳ３］
ＭＴシステム実行ユニット３２は、健常者データベース３０１に記憶されている各種計測値（診断指標）、画像パターン、ドプラ波形からパラメータを抽出し、抽出したパラメータを項目因子として多変量の健常者の状態空間を生成し、マハラノビス距離（物差し）を計算する。また、ＭＴシステム実行ユニット３２は、非健常者データベース３０３に記憶されている疾患毎の各種計測値（診断指標）、画像パターン、ドプラ波形からパラメータを抽出し、抽出したパラメータを項目因子として、疾患毎に多変量の非健常者の状態空間を生成し、疾患毎のマハラノビス距離（物差し）を計算する。さらに、ＭＴシステム実行ユニット３２は、ステップＳ２において計算された各種パラメータを項目因子として、当該患者に関する、各疾患に対応する各状態空間からのマハラノビス距離、及び健常者の状態空間からのマハラノビス距離を計算する（ステップＳ３）。
【００４９】
なお、ＭＴシステム実行ユニット３２は、各疾患に関する物差し及び患者の各疾患の非健常者の状態空間からの距離を正規化することが好ましい。
【００５０】
［疾患種別の判定：ステップＳ４］
次に、解析ユニット３４は、所定患者に関する各疾患の状態空間からの各マハラノビス距離と、各疾患の状態空間毎に計算された物差しとに基づいて、当該患者の健康状態が正常であるか異常であるか、及び異常である場合の疾患種別を判定する（ステップＳ４）。
【００５１】
図６（ａ）、（ｂ）、（ｃ）は、本ステップＳ４において解析ユニット３４が実行する判定処理の概念を説明するための図である。解析ユニット３４は、例えば同図６（ａ）に示すように、健常者であるか否かを判定する場合に、当該患者甲に関する健常者の状態空間からのマハラノビス距離ａと、健常者に該当しない（非健常者である）集合αと健常者の状態空間との間の距離（物差し）Ｄ１とを用いて相対距離ａ／Ｄ１を計算する。計算の結果、相対距離ａ／Ｄ１の値が０（ゼロ）に近いほど、解析ユニット３４は、当該項目（今の場合健常社であること）に該当する可能性が高いと判定する。
【００５２】
同じく、疾患Ａに該当するか否かを判定する場合には、同図＊（ｂ）に示すように、当該患者甲に関する疾患Ａの状態空間からのマハラノビス距離ｂと、疾患Ａに該当しない集合βと疾患Ａの状態空間との間の距離（物差し）Ｄ２とを用いて、相対距離ｂ／Ｄ２を計算する。計算の結果、相対距離ｂ／Ｄ２の値が０（ゼロ）に近いほど、解析ユニット３４は、当該項目（今の場合疾患Ａであること）に該当する可能性が高いと判定する。他の疾患の判定についても同様である。
【００５３】
この様な判定処理の結果、図７に示すように、相対的距離を基準とする当該患者甲についての各状態空間（単位空間）からの近さがわかる。同図の例では、患者甲は、疾患Ａの単位空間に最も近いため、疾患Ａである可能性が高いことになる。また、患者乙は、他の管は疾患Ｂである可能性が高く、患者丙は、健常者の状態空間に最も近いため、何らの疾患も持っていない可能性が高いことになる。
【００５４】
［診断支援情報の生成・表示：ステップＳ５］
次に、診断支援情報精製ユニット３６は、ステップＳ４において実行された正常か異常かの判定結果、及び疾患種別毎の判定結果を用いて、診断支援情報としてのヒストグラム、レーダーチャート等を生成する（ステップＳ５）。
【００５５】
図８は、本ステップＳ５において生成される診断支援情報としてのヒストグラムの一例を示した図であり、甲、乙、丙の各被検体について、「正常」、「疾患Ａ」、「疾患Ｂ」、・・・・の各項目（各疾患）についての相対マハラノビス距離（相対ＭＤ値）が示されている。
【００５６】
また、より具体的な例として、図９（ａ）、（ｂ）に、良好な結果が得られたＤＣＭのヒストグラムと、判別困難だったＡＲのヒストグラムを示した。図中の横軸１〜１２は、それぞれ健常者、ＡＲ、ＡＳ−moderete、狭心症等の項目に対応した物差しであり、縦軸は距離を表している。図９（ａ）は、物差し４に相当するＤＣＭ疾患のデータをあてはめたものであり、太線は平均値である。ＤＣＭの例では物差し４にピークがあり、良好な識別ができていることがわかる。一方、自９（ｂ）は、物差し２に相当するＡＲ疾患のデータを当てはめたものであり、太字は同じく平均値である。この例では、物差し２よりも物差し５や物差し１０にもピークがあり、うまく判別できていないことがわかる。
【００５７】
図１０は、本ステップＳ５において生成される診断支援情報としてのレーダーチャート（健常者［１］と、各疾患［２］〜［１２］の距離の平均値のレーダーチャート）の一例を示した図である。外周に行くほど各単位空間に近い表現とした。ＤＳＭやＡＳ−modereteは良好な結果を示しているが、その他は、別の疾患にもピークが発生し良好ではないことがわかる。
【００５８】
また、図１１は、本ステップＳ５において生成される診断支援情報としてのレーダーチャートの他の例を示した図である。比較的良好な疾患判別ができたＡＳ−modereteとＤＣＭの平均値と適合性の良いデータの例を示した。これらの解析に用いた診断指標中、Ｅ波速度やＥ波速度と弁輪速度の比の寄与が大きく、ドプラ診断で得られる情報のみで疾患の判別ができる可能性が高いこともわかる。
【００５９】
（実施例）
次に、本自動診断支援処理の実施例について説明する。
【００６０】
自動診断支援技術の実用性を高めるため健常者と非健常者の判別以外に、心疾患の種類を判別することを試みた。同一医師の過去数年間に亘診断データをもとにMTシステムを用いてその可能性の検討を行った。
【００６１】
すなわち、健常者と非健常者の心エコー図検査の診断指標データ約３０００人分を、健常者，ＡＲ、ＡＳ−moderate、ＤＣＭ、ＨＣＭ／ＤＣＭ、ＨＨＤ、ＭＲ、ＭＳ、ＴＲ、ＯＭＩ、先天性心疾患、狭心症、分類不能の１３種類に分けた。健常者データ数はジャックナイフ法を用いて単位空間を調整したので最終的には９００となり、分類不能を除く１１種類の疾患別データは４０〜３００程度となった。診断指標として、ドプラ計測値以外に、年齢、性別、心拍数、ＭモードやＢモード像から得られた診断指標約２０種類を組合せて解析を行った。ここでは波形パターン認識データは用いていない。この様な解析から、健常者の単位空間をもとに、特定疾患の要因効果分析を通して分離性能の評価を行った。
【００６２】
図１２に、健常者と非健常者(全疾患)のマハラノビス距離（ＭＤ値）ヒストグラムを示す。健常者の空間と非健常者の空間が重複しており区別できないことがわかる。波形パターン認識データがなく、単なる計測指標だけでは健常者と非健常者(全疾患)を区別するのが困難なことがわかった。
【００６３】
図１３は、疾患別の健常者との判別を評価するための各疾患毎（ＤＣＭ、ＭＲ、ＡＳ、ＡＲ）のヒストグラムを示した図である。同図に示すように、ＤＣＲとＭＲとでは、比較的良好な分離ができた一方、ＡＳとＡＲとでは、分離が困難であることがわかる。
【００６４】
図１４は、本実施例において生成された、マハラノビス距離（ＭＤ値）に関する疾患毎のヒストグラムを示している。として示してある。同図のマハラノビス距離は、上述したように、被検者の診断指標データを、各疾患毎の物差しで正規化して得られたものである。このヒストグラムでは、距離がゼロに近ければ近いほど、それに対応した疾患パターンと一致していることとなる。この様なヒストグラムを診断支援情報として提供することで、医師等のユーザは、距離が近いものを疾患の候補としてピックアップすることができる。
【００６５】
（効果）
以上述べた構成によれば、以下の効果を得ることができる。
【００６６】
本医用画像診断支援システム或いは医用診断支援装置によれば、超音波診断装置等の医用画像診断装置の自動計測によって得られた自動計測値（診断指標）、画像パターン、波形等を用いてパラメータを抽出し、当該パラメータをＭＴシステムの項目因子として疾患症例毎に物差しを作成し、所定の患者（すなわち、正常又は異常、或いは異常である場合の疾患種別が未知である患者）に関して、健康状態が正常であるか又は異常であるか、或いは異常である場合の疾患種別を判定し、その結果を診断支援情報として所定の形態で表示する。従って、自動計測値を用いたパラメータをＭＴシステムの項目因子とした、マハラノビス距離という客観的な定量値に基づいて、多重共線性を回避しつつ高い汎用性をもって、迅速且つ簡便に正常／異常の判別及び疾患種類の判別を実行することができる。特に、パターン認識を用いる従来技術と比較した場合、本自動診断支援装置によれば、標本の項目因子が過去の例に属さない振る舞いをするものであっても、困難なく正常／異常の判別及び疾患種類の判別を実行することができる。
【００６７】
また、本医用画像診断支援システム或いは医用診断支援装置によれば、正常／異常の判別結果及び疾患種類の判別結果を、疾患種別の相対的距離を示すレーダーチャートやヒストグラム等によって表示する。従って、医師等のユーザは、疾患種別毎のピークの発生やその形状によって、迅速且つ簡便に検査結果を視認することができる。
【００６８】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。具体的な変形例としては、例えば次のようなものがある。
【００６９】
（１）本実施形態に係る各機能は、当該処理を実行するプログラムをワークステーション等のコンピュータにインストールし、これらをメモリ上で展開することによっても実現することができる。このとき、コンピュータに当該手法を実行させることのできるプログラムは、磁気ディスク（フロッピー（登録商標）ディスク、ハードディスクなど）、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤなど）、半導体メモリなどの記録媒体に格納して頒布することも可能である。
【００７０】
（２）上記実施形態においては、自動診断支援システムとして、超音波診断装置２と自動診断支援装置３とを組み合わせる場合を例とした。しかしながら、自動診断支援装置３と組み合わせる医用画像診断装置は、超音波診断装置に拘泥されない。例えば、Ｘ線コンピュータ断層撮影装置、磁気共鳴イメージング装置、Ｘ線診断装置、核医学診断装置等を採用することも可能である。
【００７１】
また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【産業上の利用可能性】
【００７２】
以上本発明によれば、ＭＴシステムを用いて疾患症例毎の尺度を作成し、当該尺度と自動計測によって取得された項目とを用いて汎用性の高い自動診断支援を実行可能な自動診断支援装置、超音波診断装置及び自動診断支援プログラムを実現することができる。
【符号の説明】
【００７３】
１…自動診断支援システム、２…超音波診断装置、３…自動診断支援装置、２０…超音波プローブ、２１…送受信ユニット、２２…Ｂモード処理ユニット、２３…ドプラ処理ユニット、２４…画像・波形生成ユニット、２５…画像合成ユニット、２６…表示ユニット、３０…データベースユニット、３２…ＭＴシステム実行ユニット、３４…解析ユニット、３６…診断支援情報生成ユニット、３８…表示ユニット
【図１】
0005611546.tif 000006
【図２】
0005611546.tif 000007
【図３】
0005611546.tif 000008
【図４】
0005611546.tif 000009
【図５】
0005611546.tif 000010
【図６】
0005611546.tif 000011
【図７】
0005611546.tif 000012
【図８】
0005611546.tif 000013
【図９】
0005611546.tif 000014
【図１０】
0005611546.tif 000015
【図１１】
0005611546.tif 000016
【図１２】
0005611546.tif 000017
【図１３】
0005611546.tif 000018
【図１４】
0005611546.tif 000019