自動診断支援装置、超音波診断装置及び自動診断支援プログラム
【課題】 MTシステムを用いて疾患症例毎の尺度を作成し、当該尺度と自動計測によって取得された項目とを用いて汎用性の高い自動診断支援を実行可能な自動診断支援装置、超音波診断装置及び自動診断支援プログラムを提供すること。
【解決手段】 超音波診断装置等の医用画像診断装置の自動計測によって得られた自動計測値(診断指標)、画像パターン、波形等を用いてパラメータを抽出し、当該パラメータをMTシステムの項目因子として疾患症例毎に物差しを作成し、所定の患者(すなわち、正常又は異常、或いは異常である場合の疾患種別が未知である患者)に関して、健康状態が正常であるか又は異常であるか、或いは異常である場合の疾患種別を判定し、その結果を診断支援情報として所定の形態で表示するものである。
【解決手段】 超音波診断装置等の医用画像診断装置の自動計測によって得られた自動計測値(診断指標)、画像パターン、波形等を用いてパラメータを抽出し、当該パラメータをMTシステムの項目因子として疾患症例毎に物差しを作成し、所定の患者(すなわち、正常又は異常、或いは異常である場合の疾患種別が未知である患者)に関して、健康状態が正常であるか又は異常であるか、或いは異常である場合の疾患種別を判定し、その結果を診断支援情報として所定の形態で表示するものである。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、画像を用いたコンピュータ支援診断(computer-aided diagnosis: CAD)に関するもので、特に、異なるアルゴリズムに従う複数の読影結果を用いるコンピュータ支援診断に関する。
【背景技術】
【0002】
コンピュータ支援診断とは、診断における最終判断は医師が行うことを前提としつつ、医師や技師による画像診断の負荷を低減させるため、画像パターン認識技術等を用いて、入力したデジタル画像から検出した病巣候補の位置や、コンピュータを用いて解析した定量的な結果を画像上にわかりやすく提示するものである。この技術は、例えば胸部X線写真や集団検診用CT画像による肺癌の検診、乳癌マンモグラフィの自動スクリーニングシステム、胃癌の検出システム、超音波画像診断における循環器ドプラ診断等に用いられる。
【0003】
この様なコンピュータ支援診断を用いた自動診断支援装置では、近年、MT(Mahalanobis(マハラノビス)・Taguchi(タグチ))システムと呼ばれる手法が用いられている。MTシステムとは、品質工学において多変の項目(多変量)を総合して新しい尺度を構成する手法であり、その類型として、MT(Mahalanobis Taguchi)法、T(Taguchi)法、MTA(Mahalanobis Taguchi Ajoint)法、TS(Taguchi schmitt)法等がある。このMTシステムが被検体の正常/異常の識別判断に応用され始めて約5年経ち、近年では、特殊な応用として肝疾患種類の判別にSN比の要因効果パターンによる識別が試みられている。
【0004】
なお、本願発明に関連する公知文献としては、次の様なものがある。
【先行技術文献】
【特許文献】
【0005】
【非特許文献1】“心臓超音波検査法トレーニングプロトコールの評価 ―正確性と再現性について― ” J Med Ultrasonics Vol.29 No.6 p537-544 (2002)
【非特許文献2】“我国における心エコー図のルーチン検査の現状(アンケート調査 ―第2報― ” Jpn J Med Ultrasonics Vol.32 No.3 p329-337 (2005)
【非特許文献3】“心エコー図検査による加齢の心機能に対する影響の検討” J Med Ultrasonics Vol.29 No.2 p145-151 (2002)
【発明の概要】
【発明が解決しようとする課題】
【0006】
しかしながら、従来の自動診断支援装置においては次のような問題がある。
【0007】
すなわち、画像パターン認識技術としてはニューラルネットワーク、主成分分析法、多変量解析(回帰分析/因子分析)が主流であり、交互作用、モデルパターンの抽出、パラメータの変動範囲、主成分の選別方式などに統一性がない。従って、現在においても汎用的な解析手法が確立されておらず、信憑性/再現性に乏しいものとなっている。
【0008】
また、例えばMTシステムを用いた肝疾患種別の判別方法についても、健常者の状態空間をもとにした異常パターンでの、個々の症例の距離が離れる要因パターンの分析を行っている。このため、例えば肝臓では好適な自動判定ができても、他の疾患判別に役立つかどうか(汎用性があるかどうか)不明である。
【0009】
本発明は、上記事情を鑑みてなされたもので、MTシステムを用いて疾患症例毎の尺度を作成し、当該尺度と自動計測によって取得された項目とを用いて汎用性の高い自動診断支援を実行可能な自動診断支援装置、超音波診断装置及び自動診断支援プログラムを提供することを目的としている。
【課題を解決するための手段】
【0010】
本発明は、上記目的を達成するため、次のような手段を講じている。
【0011】
請求項1に記載の発明は、医用画像診断装置によって取得された所定の被検体に関する医用情報を用いて、当該被検体に関する複数のパラメータを抽出すると共に、予め取得された非健常者に関する医用情報を用いて、非健常者に関する疾患毎の前記複数のパラメータを抽出する抽出手段と、前記非健常者に関する疾患毎の前記複数のパラメータを項目因子としてMTシステムを実行することで、各疾患症例に関する物差しを計算すると共に、前記被検体に関する複数のパラメータを項目因子としてMTシステムを実行することで、前記被検体の前記各疾患の非健常者の状態空間からの距離を計算する計算手段と、前記各疾患症例に関する物差しと前記被検体の前記各疾患の非健常者の状態空間からの距離とを用いて、前記被検体の疾患種別を判定する判定手段と、前記判定の結果に基づいて、診断支援情報を生成する支援情報生成手段と、前記診断支援情報を所定の形態で表示する表示手段と、を具備することを特徴とする自動診断支援装置である。
【0012】
請求項10に記載の発明は、所定の被検体を超音波走査し、当該被検体に関する超音波データを取得する超音波データ取得手段と、前記超音波データを含む前記被検体に関する医用情報を用いて、当該被検体に関する複数のパラメータを抽出すると共に、予め取得された非健常者に関する医用情報を用いて、非健常者に関する疾患毎の前記複数のパラメータを抽出する抽出手段と、前記非健常者に関する疾患毎の前記複数のパラメータを項目因子としてMTシステムを実行することで、各疾患症例に関する物差しを計算すると共に、前記被検体に関する複数のパラメータを項目因子としてMTシステムを実行することで、前記被検体の前記各疾患の非健常者の状態空間からの距離を計算する計算手段と、前記各疾患症例に関する物差しと前記被検体の前記各疾患の非健常者の状態空間からの距離とを用いて、前記被検体の疾患種別を判定する判定手段と、前記判定の結果に基づいて、診断支援情報を生成する支援情報生成手段と、前記診断支援情報を所定の形態で表示する表示手段と、を具備することを特徴とする超音波診断装置である。
【0013】
請求項11に記載の発明は、コンピュータに、医用画像診断装置によって取得された所定の被検体に関する医用情報を用いて、当該被検体に関する複数のパラメータを抽出させると共に、予め取得された非健常者に関する医用情報を用いて、非健常者に関する疾患毎の前記複数のパラメータを抽出させる抽出機能と、前記非健常者に関する疾患毎の前記複数のパラメータを項目因子としてMTシステムを実行させることで、各疾患症例に関する物差しを計算させると共に、前記被検体に関する複数のパラメータを項目因子としてMTシステムを実行させることで、前記被検体の前記各疾患の非健常者の状態空間からの距離を計算させる計算機能と、前記各疾患症例に関する物差しと前記被検体の前記各疾患の非健常者の状態空間からの距離とを用いて、前記被検体の疾患種別を判定させる判定機能と、前記判定の結果に基づいて、診断支援情報を生成させる支援情報生成機能と、前記診断支援情報を所定の形態で表示させる表示機能と、を実現させることを特徴とする自動診断支援プログラムである。
【発明の効果】
【0014】
以上本発明によれば、MTシステムを用いて疾患症例毎の尺度を作成し、当該尺度と自動計測によって取得された項目とを用いて汎用性の高い自動診断支援を実行可能な自動診断支援装置、超音波診断装置及び自動診断支援プログラムを実現することができる。
【図面の簡単な説明】
【0015】
【図1】図1は、本実施形態に係る自動診断支援システム1のブロック構成図を示している。
【図2】図2は、マハラノビス距離の一般的な計算の流れを示したフローチャートである。
【図3】図3は、マハラノビス距離を用いた判定処理の概念を説明するための図である。
【図4】図4は、MTシステムのSN比を用いた原因診断を説明するための図である。
【図5】図5は、本自動診断支援機能に従う処理(自動診断支援処理)の流れを示したフローチャートである。
【図6】図6(a)、(b)、(c)は、ステップS4において解析ユニット34が実行する判定処理の概念を説明するための図である。
【図7】図7は、相対的距離を基準とする当該患者甲についての各状態空間(単位空間)からの近さを説明するための図である。
【図8】図8は、ステップS5において生成される診断支援情報としてのヒストグラムの一例を示した図である。
【図9】図9(a)、(b)は、良好な結果が得られたDCMのヒストグラムと、判別困難だったARのヒストグラムを示した図である。
【図10】図10は、本ステップS5において生成される診断支援情報としてのレーダーチャートの他の例を示した図である。
【図11】図11は、本ステップS5において生成される診断支援情報としてのレーダーチャートの他の例を示した図である。
【図12】図12は、健常者と非健常者(全疾患)のマハラノビス距離(MD値)ヒストグラムを示した図である。
【図13】図13は、疾患別の健常者との判別を評価するための各疾患毎(DCM、MR、AS、AR)のヒストグラムを示した図である。
【図14】図14は、本実施例において生成された、マハラノビス距離(MD値)に関する疾患毎のヒストグラムを示している。
【発明を実施するための形態】
【0016】
以下、本発明の実施形態を図面に従って説明する。なお、以下の説明において、略同一の機能及び構成を有する構成要素については、同一符号を付し、重複説明は必要な場合にのみ行う。
【0017】
図1は、本実施形態に係る自動診断支援システム1のブロック構成図を示している。同図に示すように、本自動診断支援システム1は、超音波診断装置2と、自動診断支援装置3とから構成されている。なお、本自動診断支援システム1を構成する超音波診断装置2と自動診断支援装置3とは、別体でも一体でも構わない。また、自動診断支援装置3を超音波診断装置2に内蔵させ、超音波診断装置3そのものが自動診断支援システム1となる構成であってもよい。超音波診断装置2、自動診断支援装置3のそれぞれの構成について、以下説明する。
【0018】
(超音波診断装置)
超音波診断装置2は、超音波プローブ20、送受信ユニット21、Bモード処理ユニット22、ドプラ処理ユニット23、画像・波形生成ユニット24、画像合成ユニット25、表示ユニット26を具備している。
【0019】
超音波プローブ20は、超音波送受信ユニット21からの駆動信号に基づき超音波を発生し、被検体からの反射波を電気信号に変換する複数の圧電振動子、当該圧電振動子に設けられる整合層、当該圧電振動子から後方への超音波の伝播を防止するバッキング材等を有している。当該超音波プローブ20から被検体Pに超音波が送信されると、当該送信超音波は、体内組織の音響インピーダンスの不連続面で次々と反射され、エコー信号として超音波プローブ20に受信される。このエコー信号の振幅は、反射することになった反射することになった不連続面における音響インピーダンスの差に依存する。また、送信された超音波パルスが、移動している血流や組織で反射された場合のエコーは、ドプラ効果により移動体の超音波送信方向の速度成分を依存して、周波数偏移を受ける。
【0020】
送受信ユニット21は、図示しないトリガ発生回路、遅延回路およびパルサ回路等を有している。パルサ回路では、所定のレート周波数fr Hz(周期;1/fr秒)で、送信超音波を形成するためのレートパルスが繰り返し発生される。また、遅延回路では、チャンネル毎に超音波をビーム状に集束し且つ送信指向性を決定するのに必要な遅延時間が、各レートパルスに与えられる。トリガ発生回路は、このレートパルスに基づくタイミングで、プローブ20に駆動パルスを印加する。
【0021】
また、送受信ユニット21は、図示していないアンプ回路、A/D変換器、加算器等を有している。アンプ回路では、プローブ20を介して取り込まれたエコー信号をチャンネル毎に増幅する。A/D変換器では、増幅されたエコー信号に対し受信指向性を決定するのに必要な遅延時間を与え、その後加算器において加算処理を行う。この加算により、エコー信号の受信指向性に応じた方向からの反射成分が強調され、受信指向性と送信指向性とにより超音波送受信の総合的なビームが形成される。
【0022】
Bモード処理ユニット22は、送受信ユニット21からエコー信号を受け取り、対数増幅、包絡線検波処理などを施し、信号強度が輝度の明るさで表現されるBモード画像データを生成する。
【0023】
ドプラ処理ユニット23は、送受信ユニット21から受け取ったエコー信号から速度情報を周波数解析し、ドプラ効果による血流や組織、造影剤エコー成分を抽出し、平均速度、分散、パワー等の血流情報を多点について求める。得られた血流情報は画像生成ユニット25に送られ、平均速度画像、分散画像、パワー画像、これらの組み合わせ画像としてモニター14にカラー表示される。
【0024】
画像・波形生成ユニット24は、超音波スキャンの走査線信号列を、テレビなどに代表される一般的なビデオフォーマットの走査線信号列に変換し、表示画像としての超音波診断画像やドプラ波形を生成する。また、画像・波形生成ユニット24は、画像データを格納する記憶メモリを搭載しており、3次元画像の再構成処理などを行うことが可能である。また、例えば診断の後に操作者が検査中に記録された画像を呼び出すことが可能となっている。なお、当該画像・波形生成ユニット24に入る以前のデータは、「生データ」と呼ばれることがある。
【0025】
画像合成ユニット25は、画像・波形生成ユニット24から受け取った画像を種々のパラメータの文字情報や目盛等と共に合成し、ビデオ信号としてモニター14に出力する。3次元の再構成プログラムや本発明の画像処理プログラム等も、ここに格納されており、操作者の指示などにより、これらのプログラムは起動される。
【0026】
表示ユニット26は、画像合成ユニット25からのビデオ信号に基づいて、生体内の形態学的情報や、血流情報を画像として表示する。
【0027】
なお、本超音波診断装置2は、上記各構成要素のほか、例えば、オペレータからの各種指示、条件、関心領域(ROI)の設定指示、種々の画質条件設定指示等を当該装置2にとりこむための各種スイッチ、ボタン、トラックボール、マウス、キーボード等からなる入力ユニットや、ネットワークを介して他の装置とデータを送受信するための送受信ユニット等を具備している。
【0028】
(自動診断支援装置)
自動診断支援装置3は、データベースユニット30、MTシステム実行ユニット32、解析ユニット34、診断支援情報生成ユニット36、表示ユニット38を具備している。
【0029】
データベースユニット30は、健常者データベース301と非健常者データベース303とを有する。
【0030】
健常者データベース301は、複数の健常者に関する複数の自動計測値、超音波診断装置によって取得されたBモード画像データ、超音波診断装置によって取得されたドプラ計測データ、超音波診断装置以外の医用画像診断装置(例えば、X線コンピュータ断層撮影装置、磁気共鳴イメージング装置、X線診断装置、核医学診断装置等)によって取得された画像データ、ECG波形等を記憶する。なお健常者とは、計測値が計測された時点において「心臓機能に関する疾患を患っていない」と医師により判断された人のことである。
【0031】
非健常者データベース303は、各疾患について、複数の非健常者に関する複数の自動計測値、超音波診断装置によって取得されたBモード画像データ、超音波診断装置によって取得されたドプラ計測データ、超音波診断装置以外の医用画像診断装置(例えば、X線コンピュータ断層撮影装置、磁気共鳴イメージング装置、X線診断装置、核医学診断装置等)によって取得された画像データ、ECG波形等を記憶する。ここで、非健常者とは、「後述する自動診断支援機能による疾患の判別対象となっている何れかの疾患を患っている」と医師により判断された人のことである。従って、MTシステムにおいては、非健常者は、健常者の状態空間に属さない。疾患の種類としては、例えば、心肥大、心筋梗塞、僧帽弁逆流、大動脈弁逆流、中隔欠損、貧脈・余脈等がある。しかしながら、疾患の種類はこれに限定する必要はなく、上記疾患以外の疾患であっても適用可能である。
【0032】
MTシステム実行ユニット32は、健常者データベース301に記憶されている各種計測値(診断指標)、画像パターン、ドプラ波形からパラメータを抽出し、抽出したパラメータを項目因子として多変量の健常者の状態空間を生成し、マハラノビス距離(物差し)を計算する。また、MTシステム実行ユニット32は、非健常者データベース303に記憶されている疾患毎の各種計測値(診断指標)、画像パターン、ドプラ波形からパラメータを抽出し、抽出したパラメータを項目因子として、疾患毎に多変量の非健常者の状態空間を生成し、疾患毎のマハラノビス距離(物差し)を計算する。さらに、MTシステム実行ユニット32は、超音波診断装置2によって自動計測された所定患者に関する自動計測値(診断指標)、画像パターン、ドプラ波形からパラメータを抽出し、抽出したパラメータを項目因子として、当該患者に関する、各疾患に対応する各状態空間からのマハラノビス距離を計算する。
【0033】
解析ユニット34は、所定患者に関する各疾患の状態空間からの各マハラノビス距離と、各疾患の状態空間毎に計算された物差しとに基づいて、当該患者の健康状態が正常であるか異常であるか、及び異常である場合の疾患種別を判定する。
【0034】
図2は、マハラノビス距離の一般的な計算の流れを示したフローチャートである。同図に示すように、まず、各計測項目j(1≦j≦n)について、計測値ベクトルxij=(x1j,x2j,・・・,xmj)の平均値と標準偏差σjとを算出する(ステップSa)。次に、算出されたxijの平均値と分散値σjとを用いて、式(1)に従って各計測値ベクトルxijを正規化し、正規化された計測値ベクトルXij=(X1j,X2j,・・・,Xmj)を算出する(ステップSb)。
【数1】
【0035】
なお、母集団のσは、マハラノビス距離の平均値が1となるように、式(2)に従って操作(正規化)される。また、上記ステップSa、ステップSbの処理は、全て(n個)の計測項目について行なわれる。
【0036】
次に、項目毎の正規化された計測値ベクトルXijを用いて、次の式(3)、式(4)に従って相関行列Rとその逆行列R−1が計算される(ステップSc)。
【数2】
【0037】
次に、相関行列の逆行列R−1を用いて、次の式(5)、(6)、(7)、(8)に従って、マハラノビス距離(MD)の平方D2を算出する(ステップSd)。
【数3】
【0038】
この様な計算により、例えば図3に示すように、健常者の状態空間と患者A、患者Bとの各マハラノビス距離を計算し、所定の閾値と比較することにより、患者A、患者Bのそれぞれが健常者であるか否かを判定することができる。
【0039】
なお、SN比についても、次の式(9)〜(14)を用いて計算することができる。
【数4】
【0040】
上記SN比を用いて、例えば図4に示すようなヒストグラムを生成し、原因診断を判定することも可能である。
【0041】
診断支援情報生成ユニット36は、解析ユニット34において実行された解析処理の結果に基づいて、例えば各疾患を患っている危険性(類似度が高いものほど危険性が高い)を示すレーダーチャートRC等を診断支援情報として生成する。
【0042】
表示ユニット38は、診断支援情報生成ユニット36において生成された診断支援情報を、所定の形態で表示する。
【0043】
(自動診断支援機能)
次に、本自動診断支援システム1が有する、自動診断支援機能について説明する。この機能は、超音波診断装置等の医用画像診断装置の自動計測によって得られた自動計測値(診断指標)、画像パターン、波形等を用いてパラメータを抽出し、当該パラメータをMTシステムの項目因子として疾患症例毎に物差しを作成し、所定の患者(すなわち、正常又は異常、或いは異常である場合の疾患種別が未知である患者(被検体))に関して、健康状態が正常であるか又は異常であるか、或いは異常である場合の疾患種別を判定し、その結果を診断支援情報として所定の形態で表示するものである。
【0044】
図3は、本自動診断支援機能に従う処理(自動診断支援処理)の流れを示したフローチャートである。以下、各ステップにおいて実行される処理の内容について説明する。なお、本実施形態においては、自動診断支援装置3の動作は、例えば超音波診断装置2による医用情報の取得の完了に連動して開始するように制御されるものとする。
【0045】
[自動計測・ドプラ波形・Bモード画像等の取得:ステップS1]
まず、超音波診断装置2等を用いて、所定の患者についての自動計測が実行される。この自動計測により、ECG波形、ドプラ波形、Bモード画像等が取得される(ステップS1)。
【0046】
[パラメータ抽出・項目因子の決定:ステップS2]
次に、MTシステム実行ユニット32は、ステップS1において取得されたECG波形、ドプラ波形、Bモード画像、各種診断指標等を用いて、MTシステムの項目因子とすべき各種パラメータを計算する。具体的には、MTシステム実行ユニット32は、取得された超音波画像のパターン認識データ(Bモード画像の二次元の相関関数等)、微分積分法等を用いた波形パターン認識データ、ECG波形やドプラ計測に基づく診断指標を計算する(ステップS2)。
【0047】
なお、この項目因子とすべきパラメータの計算は、ステップS1において取得された情報に限定する必要はなく、例えば、事前に他のモダリティによって取得された画像データ、過去の診断において取得されたデータ等を採用するようにしてもよい。
【0048】
[疾患症例毎の物差し、及び当該患者に関するマハラノビス距離の計算:ステップS3]
MTシステム実行ユニット32は、健常者データベース301に記憶されている各種計測値(診断指標)、画像パターン、ドプラ波形からパラメータを抽出し、抽出したパラメータを項目因子として多変量の健常者の状態空間を生成し、マハラノビス距離(物差し)を計算する。また、MTシステム実行ユニット32は、非健常者データベース303に記憶されている疾患毎の各種計測値(診断指標)、画像パターン、ドプラ波形からパラメータを抽出し、抽出したパラメータを項目因子として、疾患毎に多変量の非健常者の状態空間を生成し、疾患毎のマハラノビス距離(物差し)を計算する。さらに、MTシステム実行ユニット32は、ステップS2において計算された各種パラメータを項目因子として、当該患者に関する、各疾患に対応する各状態空間からのマハラノビス距離、及び健常者の状態空間からのマハラノビス距離を計算する(ステップS3)。
【0049】
なお、MTシステム実行ユニット32は、各疾患に関する物差し及び患者の各疾患の非健常者の状態空間からの距離を正規化することが好ましい。
【0050】
[疾患種別の判定:ステップS4]
次に、解析ユニット34は、所定患者に関する各疾患の状態空間からの各マハラノビス距離と、各疾患の状態空間毎に計算された物差しとに基づいて、当該患者の健康状態が正常であるか異常であるか、及び異常である場合の疾患種別を判定する(ステップS4)。
【0051】
図6(a)、(b)、(c)は、本ステップS4において解析ユニット34が実行する判定処理の概念を説明するための図である。解析ユニット34は、例えば同図6(a)に示すように、健常者であるか否かを判定する場合に、当該患者甲に関する健常者の状態空間からのマハラノビス距離aと、健常者に該当しない(非健常者である)集合αと健常者の状態空間との間の距離(物差し)D1とを用いて相対距離a/D1を計算する。計算の結果、相対距離a/D1の値が0(ゼロ)に近いほど、解析ユニット34は、当該項目(今の場合健常社であること)に該当する可能性が高いと判定する。
【0052】
同じく、疾患Aに該当するか否かを判定する場合には、同図*(b)に示すように、当該患者甲に関する疾患Aの状態空間からのマハラノビス距離bと、疾患Aに該当しない集合βと疾患Aの状態空間との間の距離(物差し)D2とを用いて、相対距離b/D2を計算する。計算の結果、相対距離b/D2の値が0(ゼロ)に近いほど、解析ユニット34は、当該項目(今の場合疾患Aであること)に該当する可能性が高いと判定する。他の疾患の判定についても同様である。
【0053】
この様な判定処理の結果、図7に示すように、相対的距離を基準とする当該患者甲についての各状態空間(単位空間)からの近さがわかる。同図の例では、患者甲は、疾患Aの単位空間に最も近いため、疾患Aである可能性が高いことになる。また、患者乙は、他の管は疾患Bである可能性が高く、患者丙は、健常者の状態空間に最も近いため、何らの疾患も持っていない可能性が高いことになる。
【0054】
[診断支援情報の生成・表示:ステップS5]
次に、診断支援情報精製ユニット36は、ステップS4において実行された正常か異常かの判定結果、及び疾患種別毎の判定結果を用いて、診断支援情報としてのヒストグラム、レーダーチャート等を生成する(ステップS5)。
【0055】
図8は、本ステップS5において生成される診断支援情報としてのヒストグラムの一例を示した図であり、甲、乙、丙の各被検体について、「正常」、「疾患A」、「疾患B」、・・・・の各項目(各疾患)についての相対マハラノビス距離(相対MD値)が示されている。
【0056】
また、より具体的な例として、図9(a)、(b)に、良好な結果が得られたDCMのヒストグラムと、判別困難だったARのヒストグラムを示した。図中の横軸1〜12は、それぞれ健常者、AR、AS−moderete、狭心症等の項目に対応した物差しであり、縦軸は距離を表している。図9(a)は、物差し4に相当するDCM疾患のデータをあてはめたものであり、太線は平均値である。DCMの例では物差し4にピークがあり、良好な識別ができていることがわかる。一方、自9(b)は、物差し2に相当するAR疾患のデータを当てはめたものであり、太字は同じく平均値である。この例では、物差し2よりも物差し5や物差し10にもピークがあり、うまく判別できていないことがわかる。
【0057】
図10は、本ステップS5において生成される診断支援情報としてのレーダーチャート(健常者[1]と、各疾患[2]〜[12]の距離の平均値のレーダーチャート)の一例を示した図である。外周に行くほど各単位空間に近い表現とした。DSMやAS−modereteは良好な結果を示しているが、その他は、別の疾患にもピークが発生し良好ではないことがわかる。
【0058】
また、図11は、本ステップS5において生成される診断支援情報としてのレーダーチャートの他の例を示した図である。比較的良好な疾患判別ができたAS−modereteとDCMの平均値と適合性の良いデータの例を示した。これらの解析に用いた診断指標中、E波速度やE波速度と弁輪速度の比の寄与が大きく、ドプラ診断で得られる情報のみで疾患の判別ができる可能性が高いこともわかる。
【0059】
(実施例)
次に、本自動診断支援処理の実施例について説明する。
【0060】
自動診断支援技術の実用性を高めるため健常者と非健常者の判別以外に、心疾患の種類を判別することを試みた。同一医師の過去数年間に亘診断データをもとにMTシステムを用いてその可能性の検討を行った。
【0061】
すなわち、健常者と非健常者の心エコー図検査の診断指標データ約3000人分を、健常者,AR、AS−moderate、DCM、HCM/DCM、HHD、MR、MS、TR、OMI、先天性心疾患、狭心症、分類不能の13種類に分けた。健常者データ数はジャックナイフ法を用いて単位空間を調整したので最終的には900となり、分類不能を除く11種類の疾患別データは40〜300程度となった。診断指標として、ドプラ計測値以外に、年齢、性別、心拍数、MモードやBモード像から得られた診断指標約20種類を組合せて解析を行った。ここでは波形パターン認識データは用いていない。この様な解析から、健常者の単位空間をもとに、特定疾患の要因効果分析を通して分離性能の評価を行った。
【0062】
図12に、健常者と非健常者(全疾患)のマハラノビス距離(MD値)ヒストグラムを示す。健常者の空間と非健常者の空間が重複しており区別できないことがわかる。波形パターン認識データがなく、単なる計測指標だけでは健常者と非健常者(全疾患)を区別するのが困難なことがわかった。
【0063】
図13は、疾患別の健常者との判別を評価するための各疾患毎(DCM、MR、AS、AR)のヒストグラムを示した図である。同図に示すように、DCRとMRとでは、比較的良好な分離ができた一方、ASとARとでは、分離が困難であることがわかる。
【0064】
図14は、本実施例において生成された、マハラノビス距離(MD値)に関する疾患毎のヒストグラムを示している。として示してある。同図のマハラノビス距離は、上述したように、被検者の診断指標データを、各疾患毎の物差しで正規化して得られたものである。このヒストグラムでは、距離がゼロに近ければ近いほど、それに対応した疾患パターンと一致していることとなる。この様なヒストグラムを診断支援情報として提供することで、医師等のユーザは、距離が近いものを疾患の候補としてピックアップすることができる。
【0065】
(効果)
以上述べた構成によれば、以下の効果を得ることができる。
【0066】
本医用画像診断支援システム或いは医用診断支援装置によれば、超音波診断装置等の医用画像診断装置の自動計測によって得られた自動計測値(診断指標)、画像パターン、波形等を用いてパラメータを抽出し、当該パラメータをMTシステムの項目因子として疾患症例毎に物差しを作成し、所定の患者(すなわち、正常又は異常、或いは異常である場合の疾患種別が未知である患者)に関して、健康状態が正常であるか又は異常であるか、或いは異常である場合の疾患種別を判定し、その結果を診断支援情報として所定の形態で表示する。従って、自動計測値を用いたパラメータをMTシステムの項目因子とした、マハラノビス距離という客観的な定量値に基づいて、多重共線性を回避しつつ高い汎用性をもって、迅速且つ簡便に正常/異常の判別及び疾患種類の判別を実行することができる。特に、パターン認識を用いる従来技術と比較した場合、本自動診断支援装置によれば、標本の項目因子が過去の例に属さない振る舞いをするものであっても、困難なく正常/異常の判別及び疾患種類の判別を実行することができる。
【0067】
また、本医用画像診断支援システム或いは医用診断支援装置によれば、正常/異常の判別結果及び疾患種類の判別結果を、疾患種別の相対的距離を示すレーダーチャートやヒストグラム等によって表示する。従って、医師等のユーザは、疾患種別毎のピークの発生やその形状によって、迅速且つ簡便に検査結果を視認することができる。
【0068】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。具体的な変形例としては、例えば次のようなものがある。
【0069】
(1)本実施形態に係る各機能は、当該処理を実行するプログラムをワークステーション等のコンピュータにインストールし、これらをメモリ上で展開することによっても実現することができる。このとき、コンピュータに当該手法を実行させることのできるプログラムは、磁気ディスク(フロッピー(登録商標)ディスク、ハードディスクなど)、光ディスク(CD−ROM、DVDなど)、半導体メモリなどの記録媒体に格納して頒布することも可能である。
【0070】
(2)上記実施形態においては、自動診断支援システムとして、超音波診断装置2と自動診断支援装置3とを組み合わせる場合を例とした。しかしながら、自動診断支援装置3と組み合わせる医用画像診断装置は、超音波診断装置に拘泥されない。例えば、X線コンピュータ断層撮影装置、磁気共鳴イメージング装置、X線診断装置、核医学診断装置等を採用することも可能である。
【0071】
また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【産業上の利用可能性】
【0072】
以上本発明によれば、MTシステムを用いて疾患症例毎の尺度を作成し、当該尺度と自動計測によって取得された項目とを用いて汎用性の高い自動診断支援を実行可能な自動診断支援装置、超音波診断装置及び自動診断支援プログラムを実現することができる。
【符号の説明】
【0073】
1…自動診断支援システム、2…超音波診断装置、3…自動診断支援装置、20…超音波プローブ、21…送受信ユニット、22…Bモード処理ユニット、23…ドプラ処理ユニット、24…画像・波形生成ユニット、25…画像合成ユニット、26…表示ユニット、30…データベースユニット、32…MTシステム実行ユニット、34…解析ユニット、36…診断支援情報生成ユニット、38…表示ユニット
【技術分野】
【0001】
本発明は、画像を用いたコンピュータ支援診断(computer-aided diagnosis: CAD)に関するもので、特に、異なるアルゴリズムに従う複数の読影結果を用いるコンピュータ支援診断に関する。
【背景技術】
【0002】
コンピュータ支援診断とは、診断における最終判断は医師が行うことを前提としつつ、医師や技師による画像診断の負荷を低減させるため、画像パターン認識技術等を用いて、入力したデジタル画像から検出した病巣候補の位置や、コンピュータを用いて解析した定量的な結果を画像上にわかりやすく提示するものである。この技術は、例えば胸部X線写真や集団検診用CT画像による肺癌の検診、乳癌マンモグラフィの自動スクリーニングシステム、胃癌の検出システム、超音波画像診断における循環器ドプラ診断等に用いられる。
【0003】
この様なコンピュータ支援診断を用いた自動診断支援装置では、近年、MT(Mahalanobis(マハラノビス)・Taguchi(タグチ))システムと呼ばれる手法が用いられている。MTシステムとは、品質工学において多変の項目(多変量)を総合して新しい尺度を構成する手法であり、その類型として、MT(Mahalanobis Taguchi)法、T(Taguchi)法、MTA(Mahalanobis Taguchi Ajoint)法、TS(Taguchi schmitt)法等がある。このMTシステムが被検体の正常/異常の識別判断に応用され始めて約5年経ち、近年では、特殊な応用として肝疾患種類の判別にSN比の要因効果パターンによる識別が試みられている。
【0004】
なお、本願発明に関連する公知文献としては、次の様なものがある。
【先行技術文献】
【特許文献】
【0005】
【非特許文献1】“心臓超音波検査法トレーニングプロトコールの評価 ―正確性と再現性について― ” J Med Ultrasonics Vol.29 No.6 p537-544 (2002)
【非特許文献2】“我国における心エコー図のルーチン検査の現状(アンケート調査 ―第2報― ” Jpn J Med Ultrasonics Vol.32 No.3 p329-337 (2005)
【非特許文献3】“心エコー図検査による加齢の心機能に対する影響の検討” J Med Ultrasonics Vol.29 No.2 p145-151 (2002)
【発明の概要】
【発明が解決しようとする課題】
【0006】
しかしながら、従来の自動診断支援装置においては次のような問題がある。
【0007】
すなわち、画像パターン認識技術としてはニューラルネットワーク、主成分分析法、多変量解析(回帰分析/因子分析)が主流であり、交互作用、モデルパターンの抽出、パラメータの変動範囲、主成分の選別方式などに統一性がない。従って、現在においても汎用的な解析手法が確立されておらず、信憑性/再現性に乏しいものとなっている。
【0008】
また、例えばMTシステムを用いた肝疾患種別の判別方法についても、健常者の状態空間をもとにした異常パターンでの、個々の症例の距離が離れる要因パターンの分析を行っている。このため、例えば肝臓では好適な自動判定ができても、他の疾患判別に役立つかどうか(汎用性があるかどうか)不明である。
【0009】
本発明は、上記事情を鑑みてなされたもので、MTシステムを用いて疾患症例毎の尺度を作成し、当該尺度と自動計測によって取得された項目とを用いて汎用性の高い自動診断支援を実行可能な自動診断支援装置、超音波診断装置及び自動診断支援プログラムを提供することを目的としている。
【課題を解決するための手段】
【0010】
本発明は、上記目的を達成するため、次のような手段を講じている。
【0011】
請求項1に記載の発明は、医用画像診断装置によって取得された所定の被検体に関する医用情報を用いて、当該被検体に関する複数のパラメータを抽出すると共に、予め取得された非健常者に関する医用情報を用いて、非健常者に関する疾患毎の前記複数のパラメータを抽出する抽出手段と、前記非健常者に関する疾患毎の前記複数のパラメータを項目因子としてMTシステムを実行することで、各疾患症例に関する物差しを計算すると共に、前記被検体に関する複数のパラメータを項目因子としてMTシステムを実行することで、前記被検体の前記各疾患の非健常者の状態空間からの距離を計算する計算手段と、前記各疾患症例に関する物差しと前記被検体の前記各疾患の非健常者の状態空間からの距離とを用いて、前記被検体の疾患種別を判定する判定手段と、前記判定の結果に基づいて、診断支援情報を生成する支援情報生成手段と、前記診断支援情報を所定の形態で表示する表示手段と、を具備することを特徴とする自動診断支援装置である。
【0012】
請求項10に記載の発明は、所定の被検体を超音波走査し、当該被検体に関する超音波データを取得する超音波データ取得手段と、前記超音波データを含む前記被検体に関する医用情報を用いて、当該被検体に関する複数のパラメータを抽出すると共に、予め取得された非健常者に関する医用情報を用いて、非健常者に関する疾患毎の前記複数のパラメータを抽出する抽出手段と、前記非健常者に関する疾患毎の前記複数のパラメータを項目因子としてMTシステムを実行することで、各疾患症例に関する物差しを計算すると共に、前記被検体に関する複数のパラメータを項目因子としてMTシステムを実行することで、前記被検体の前記各疾患の非健常者の状態空間からの距離を計算する計算手段と、前記各疾患症例に関する物差しと前記被検体の前記各疾患の非健常者の状態空間からの距離とを用いて、前記被検体の疾患種別を判定する判定手段と、前記判定の結果に基づいて、診断支援情報を生成する支援情報生成手段と、前記診断支援情報を所定の形態で表示する表示手段と、を具備することを特徴とする超音波診断装置である。
【0013】
請求項11に記載の発明は、コンピュータに、医用画像診断装置によって取得された所定の被検体に関する医用情報を用いて、当該被検体に関する複数のパラメータを抽出させると共に、予め取得された非健常者に関する医用情報を用いて、非健常者に関する疾患毎の前記複数のパラメータを抽出させる抽出機能と、前記非健常者に関する疾患毎の前記複数のパラメータを項目因子としてMTシステムを実行させることで、各疾患症例に関する物差しを計算させると共に、前記被検体に関する複数のパラメータを項目因子としてMTシステムを実行させることで、前記被検体の前記各疾患の非健常者の状態空間からの距離を計算させる計算機能と、前記各疾患症例に関する物差しと前記被検体の前記各疾患の非健常者の状態空間からの距離とを用いて、前記被検体の疾患種別を判定させる判定機能と、前記判定の結果に基づいて、診断支援情報を生成させる支援情報生成機能と、前記診断支援情報を所定の形態で表示させる表示機能と、を実現させることを特徴とする自動診断支援プログラムである。
【発明の効果】
【0014】
以上本発明によれば、MTシステムを用いて疾患症例毎の尺度を作成し、当該尺度と自動計測によって取得された項目とを用いて汎用性の高い自動診断支援を実行可能な自動診断支援装置、超音波診断装置及び自動診断支援プログラムを実現することができる。
【図面の簡単な説明】
【0015】
【図1】図1は、本実施形態に係る自動診断支援システム1のブロック構成図を示している。
【図2】図2は、マハラノビス距離の一般的な計算の流れを示したフローチャートである。
【図3】図3は、マハラノビス距離を用いた判定処理の概念を説明するための図である。
【図4】図4は、MTシステムのSN比を用いた原因診断を説明するための図である。
【図5】図5は、本自動診断支援機能に従う処理(自動診断支援処理)の流れを示したフローチャートである。
【図6】図6(a)、(b)、(c)は、ステップS4において解析ユニット34が実行する判定処理の概念を説明するための図である。
【図7】図7は、相対的距離を基準とする当該患者甲についての各状態空間(単位空間)からの近さを説明するための図である。
【図8】図8は、ステップS5において生成される診断支援情報としてのヒストグラムの一例を示した図である。
【図9】図9(a)、(b)は、良好な結果が得られたDCMのヒストグラムと、判別困難だったARのヒストグラムを示した図である。
【図10】図10は、本ステップS5において生成される診断支援情報としてのレーダーチャートの他の例を示した図である。
【図11】図11は、本ステップS5において生成される診断支援情報としてのレーダーチャートの他の例を示した図である。
【図12】図12は、健常者と非健常者(全疾患)のマハラノビス距離(MD値)ヒストグラムを示した図である。
【図13】図13は、疾患別の健常者との判別を評価するための各疾患毎(DCM、MR、AS、AR)のヒストグラムを示した図である。
【図14】図14は、本実施例において生成された、マハラノビス距離(MD値)に関する疾患毎のヒストグラムを示している。
【発明を実施するための形態】
【0016】
以下、本発明の実施形態を図面に従って説明する。なお、以下の説明において、略同一の機能及び構成を有する構成要素については、同一符号を付し、重複説明は必要な場合にのみ行う。
【0017】
図1は、本実施形態に係る自動診断支援システム1のブロック構成図を示している。同図に示すように、本自動診断支援システム1は、超音波診断装置2と、自動診断支援装置3とから構成されている。なお、本自動診断支援システム1を構成する超音波診断装置2と自動診断支援装置3とは、別体でも一体でも構わない。また、自動診断支援装置3を超音波診断装置2に内蔵させ、超音波診断装置3そのものが自動診断支援システム1となる構成であってもよい。超音波診断装置2、自動診断支援装置3のそれぞれの構成について、以下説明する。
【0018】
(超音波診断装置)
超音波診断装置2は、超音波プローブ20、送受信ユニット21、Bモード処理ユニット22、ドプラ処理ユニット23、画像・波形生成ユニット24、画像合成ユニット25、表示ユニット26を具備している。
【0019】
超音波プローブ20は、超音波送受信ユニット21からの駆動信号に基づき超音波を発生し、被検体からの反射波を電気信号に変換する複数の圧電振動子、当該圧電振動子に設けられる整合層、当該圧電振動子から後方への超音波の伝播を防止するバッキング材等を有している。当該超音波プローブ20から被検体Pに超音波が送信されると、当該送信超音波は、体内組織の音響インピーダンスの不連続面で次々と反射され、エコー信号として超音波プローブ20に受信される。このエコー信号の振幅は、反射することになった反射することになった不連続面における音響インピーダンスの差に依存する。また、送信された超音波パルスが、移動している血流や組織で反射された場合のエコーは、ドプラ効果により移動体の超音波送信方向の速度成分を依存して、周波数偏移を受ける。
【0020】
送受信ユニット21は、図示しないトリガ発生回路、遅延回路およびパルサ回路等を有している。パルサ回路では、所定のレート周波数fr Hz(周期;1/fr秒)で、送信超音波を形成するためのレートパルスが繰り返し発生される。また、遅延回路では、チャンネル毎に超音波をビーム状に集束し且つ送信指向性を決定するのに必要な遅延時間が、各レートパルスに与えられる。トリガ発生回路は、このレートパルスに基づくタイミングで、プローブ20に駆動パルスを印加する。
【0021】
また、送受信ユニット21は、図示していないアンプ回路、A/D変換器、加算器等を有している。アンプ回路では、プローブ20を介して取り込まれたエコー信号をチャンネル毎に増幅する。A/D変換器では、増幅されたエコー信号に対し受信指向性を決定するのに必要な遅延時間を与え、その後加算器において加算処理を行う。この加算により、エコー信号の受信指向性に応じた方向からの反射成分が強調され、受信指向性と送信指向性とにより超音波送受信の総合的なビームが形成される。
【0022】
Bモード処理ユニット22は、送受信ユニット21からエコー信号を受け取り、対数増幅、包絡線検波処理などを施し、信号強度が輝度の明るさで表現されるBモード画像データを生成する。
【0023】
ドプラ処理ユニット23は、送受信ユニット21から受け取ったエコー信号から速度情報を周波数解析し、ドプラ効果による血流や組織、造影剤エコー成分を抽出し、平均速度、分散、パワー等の血流情報を多点について求める。得られた血流情報は画像生成ユニット25に送られ、平均速度画像、分散画像、パワー画像、これらの組み合わせ画像としてモニター14にカラー表示される。
【0024】
画像・波形生成ユニット24は、超音波スキャンの走査線信号列を、テレビなどに代表される一般的なビデオフォーマットの走査線信号列に変換し、表示画像としての超音波診断画像やドプラ波形を生成する。また、画像・波形生成ユニット24は、画像データを格納する記憶メモリを搭載しており、3次元画像の再構成処理などを行うことが可能である。また、例えば診断の後に操作者が検査中に記録された画像を呼び出すことが可能となっている。なお、当該画像・波形生成ユニット24に入る以前のデータは、「生データ」と呼ばれることがある。
【0025】
画像合成ユニット25は、画像・波形生成ユニット24から受け取った画像を種々のパラメータの文字情報や目盛等と共に合成し、ビデオ信号としてモニター14に出力する。3次元の再構成プログラムや本発明の画像処理プログラム等も、ここに格納されており、操作者の指示などにより、これらのプログラムは起動される。
【0026】
表示ユニット26は、画像合成ユニット25からのビデオ信号に基づいて、生体内の形態学的情報や、血流情報を画像として表示する。
【0027】
なお、本超音波診断装置2は、上記各構成要素のほか、例えば、オペレータからの各種指示、条件、関心領域(ROI)の設定指示、種々の画質条件設定指示等を当該装置2にとりこむための各種スイッチ、ボタン、トラックボール、マウス、キーボード等からなる入力ユニットや、ネットワークを介して他の装置とデータを送受信するための送受信ユニット等を具備している。
【0028】
(自動診断支援装置)
自動診断支援装置3は、データベースユニット30、MTシステム実行ユニット32、解析ユニット34、診断支援情報生成ユニット36、表示ユニット38を具備している。
【0029】
データベースユニット30は、健常者データベース301と非健常者データベース303とを有する。
【0030】
健常者データベース301は、複数の健常者に関する複数の自動計測値、超音波診断装置によって取得されたBモード画像データ、超音波診断装置によって取得されたドプラ計測データ、超音波診断装置以外の医用画像診断装置(例えば、X線コンピュータ断層撮影装置、磁気共鳴イメージング装置、X線診断装置、核医学診断装置等)によって取得された画像データ、ECG波形等を記憶する。なお健常者とは、計測値が計測された時点において「心臓機能に関する疾患を患っていない」と医師により判断された人のことである。
【0031】
非健常者データベース303は、各疾患について、複数の非健常者に関する複数の自動計測値、超音波診断装置によって取得されたBモード画像データ、超音波診断装置によって取得されたドプラ計測データ、超音波診断装置以外の医用画像診断装置(例えば、X線コンピュータ断層撮影装置、磁気共鳴イメージング装置、X線診断装置、核医学診断装置等)によって取得された画像データ、ECG波形等を記憶する。ここで、非健常者とは、「後述する自動診断支援機能による疾患の判別対象となっている何れかの疾患を患っている」と医師により判断された人のことである。従って、MTシステムにおいては、非健常者は、健常者の状態空間に属さない。疾患の種類としては、例えば、心肥大、心筋梗塞、僧帽弁逆流、大動脈弁逆流、中隔欠損、貧脈・余脈等がある。しかしながら、疾患の種類はこれに限定する必要はなく、上記疾患以外の疾患であっても適用可能である。
【0032】
MTシステム実行ユニット32は、健常者データベース301に記憶されている各種計測値(診断指標)、画像パターン、ドプラ波形からパラメータを抽出し、抽出したパラメータを項目因子として多変量の健常者の状態空間を生成し、マハラノビス距離(物差し)を計算する。また、MTシステム実行ユニット32は、非健常者データベース303に記憶されている疾患毎の各種計測値(診断指標)、画像パターン、ドプラ波形からパラメータを抽出し、抽出したパラメータを項目因子として、疾患毎に多変量の非健常者の状態空間を生成し、疾患毎のマハラノビス距離(物差し)を計算する。さらに、MTシステム実行ユニット32は、超音波診断装置2によって自動計測された所定患者に関する自動計測値(診断指標)、画像パターン、ドプラ波形からパラメータを抽出し、抽出したパラメータを項目因子として、当該患者に関する、各疾患に対応する各状態空間からのマハラノビス距離を計算する。
【0033】
解析ユニット34は、所定患者に関する各疾患の状態空間からの各マハラノビス距離と、各疾患の状態空間毎に計算された物差しとに基づいて、当該患者の健康状態が正常であるか異常であるか、及び異常である場合の疾患種別を判定する。
【0034】
図2は、マハラノビス距離の一般的な計算の流れを示したフローチャートである。同図に示すように、まず、各計測項目j(1≦j≦n)について、計測値ベクトルxij=(x1j,x2j,・・・,xmj)の平均値と標準偏差σjとを算出する(ステップSa)。次に、算出されたxijの平均値と分散値σjとを用いて、式(1)に従って各計測値ベクトルxijを正規化し、正規化された計測値ベクトルXij=(X1j,X2j,・・・,Xmj)を算出する(ステップSb)。
【数1】
【0035】
なお、母集団のσは、マハラノビス距離の平均値が1となるように、式(2)に従って操作(正規化)される。また、上記ステップSa、ステップSbの処理は、全て(n個)の計測項目について行なわれる。
【0036】
次に、項目毎の正規化された計測値ベクトルXijを用いて、次の式(3)、式(4)に従って相関行列Rとその逆行列R−1が計算される(ステップSc)。
【数2】
【0037】
次に、相関行列の逆行列R−1を用いて、次の式(5)、(6)、(7)、(8)に従って、マハラノビス距離(MD)の平方D2を算出する(ステップSd)。
【数3】
【0038】
この様な計算により、例えば図3に示すように、健常者の状態空間と患者A、患者Bとの各マハラノビス距離を計算し、所定の閾値と比較することにより、患者A、患者Bのそれぞれが健常者であるか否かを判定することができる。
【0039】
なお、SN比についても、次の式(9)〜(14)を用いて計算することができる。
【数4】
【0040】
上記SN比を用いて、例えば図4に示すようなヒストグラムを生成し、原因診断を判定することも可能である。
【0041】
診断支援情報生成ユニット36は、解析ユニット34において実行された解析処理の結果に基づいて、例えば各疾患を患っている危険性(類似度が高いものほど危険性が高い)を示すレーダーチャートRC等を診断支援情報として生成する。
【0042】
表示ユニット38は、診断支援情報生成ユニット36において生成された診断支援情報を、所定の形態で表示する。
【0043】
(自動診断支援機能)
次に、本自動診断支援システム1が有する、自動診断支援機能について説明する。この機能は、超音波診断装置等の医用画像診断装置の自動計測によって得られた自動計測値(診断指標)、画像パターン、波形等を用いてパラメータを抽出し、当該パラメータをMTシステムの項目因子として疾患症例毎に物差しを作成し、所定の患者(すなわち、正常又は異常、或いは異常である場合の疾患種別が未知である患者(被検体))に関して、健康状態が正常であるか又は異常であるか、或いは異常である場合の疾患種別を判定し、その結果を診断支援情報として所定の形態で表示するものである。
【0044】
図3は、本自動診断支援機能に従う処理(自動診断支援処理)の流れを示したフローチャートである。以下、各ステップにおいて実行される処理の内容について説明する。なお、本実施形態においては、自動診断支援装置3の動作は、例えば超音波診断装置2による医用情報の取得の完了に連動して開始するように制御されるものとする。
【0045】
[自動計測・ドプラ波形・Bモード画像等の取得:ステップS1]
まず、超音波診断装置2等を用いて、所定の患者についての自動計測が実行される。この自動計測により、ECG波形、ドプラ波形、Bモード画像等が取得される(ステップS1)。
【0046】
[パラメータ抽出・項目因子の決定:ステップS2]
次に、MTシステム実行ユニット32は、ステップS1において取得されたECG波形、ドプラ波形、Bモード画像、各種診断指標等を用いて、MTシステムの項目因子とすべき各種パラメータを計算する。具体的には、MTシステム実行ユニット32は、取得された超音波画像のパターン認識データ(Bモード画像の二次元の相関関数等)、微分積分法等を用いた波形パターン認識データ、ECG波形やドプラ計測に基づく診断指標を計算する(ステップS2)。
【0047】
なお、この項目因子とすべきパラメータの計算は、ステップS1において取得された情報に限定する必要はなく、例えば、事前に他のモダリティによって取得された画像データ、過去の診断において取得されたデータ等を採用するようにしてもよい。
【0048】
[疾患症例毎の物差し、及び当該患者に関するマハラノビス距離の計算:ステップS3]
MTシステム実行ユニット32は、健常者データベース301に記憶されている各種計測値(診断指標)、画像パターン、ドプラ波形からパラメータを抽出し、抽出したパラメータを項目因子として多変量の健常者の状態空間を生成し、マハラノビス距離(物差し)を計算する。また、MTシステム実行ユニット32は、非健常者データベース303に記憶されている疾患毎の各種計測値(診断指標)、画像パターン、ドプラ波形からパラメータを抽出し、抽出したパラメータを項目因子として、疾患毎に多変量の非健常者の状態空間を生成し、疾患毎のマハラノビス距離(物差し)を計算する。さらに、MTシステム実行ユニット32は、ステップS2において計算された各種パラメータを項目因子として、当該患者に関する、各疾患に対応する各状態空間からのマハラノビス距離、及び健常者の状態空間からのマハラノビス距離を計算する(ステップS3)。
【0049】
なお、MTシステム実行ユニット32は、各疾患に関する物差し及び患者の各疾患の非健常者の状態空間からの距離を正規化することが好ましい。
【0050】
[疾患種別の判定:ステップS4]
次に、解析ユニット34は、所定患者に関する各疾患の状態空間からの各マハラノビス距離と、各疾患の状態空間毎に計算された物差しとに基づいて、当該患者の健康状態が正常であるか異常であるか、及び異常である場合の疾患種別を判定する(ステップS4)。
【0051】
図6(a)、(b)、(c)は、本ステップS4において解析ユニット34が実行する判定処理の概念を説明するための図である。解析ユニット34は、例えば同図6(a)に示すように、健常者であるか否かを判定する場合に、当該患者甲に関する健常者の状態空間からのマハラノビス距離aと、健常者に該当しない(非健常者である)集合αと健常者の状態空間との間の距離(物差し)D1とを用いて相対距離a/D1を計算する。計算の結果、相対距離a/D1の値が0(ゼロ)に近いほど、解析ユニット34は、当該項目(今の場合健常社であること)に該当する可能性が高いと判定する。
【0052】
同じく、疾患Aに該当するか否かを判定する場合には、同図*(b)に示すように、当該患者甲に関する疾患Aの状態空間からのマハラノビス距離bと、疾患Aに該当しない集合βと疾患Aの状態空間との間の距離(物差し)D2とを用いて、相対距離b/D2を計算する。計算の結果、相対距離b/D2の値が0(ゼロ)に近いほど、解析ユニット34は、当該項目(今の場合疾患Aであること)に該当する可能性が高いと判定する。他の疾患の判定についても同様である。
【0053】
この様な判定処理の結果、図7に示すように、相対的距離を基準とする当該患者甲についての各状態空間(単位空間)からの近さがわかる。同図の例では、患者甲は、疾患Aの単位空間に最も近いため、疾患Aである可能性が高いことになる。また、患者乙は、他の管は疾患Bである可能性が高く、患者丙は、健常者の状態空間に最も近いため、何らの疾患も持っていない可能性が高いことになる。
【0054】
[診断支援情報の生成・表示:ステップS5]
次に、診断支援情報精製ユニット36は、ステップS4において実行された正常か異常かの判定結果、及び疾患種別毎の判定結果を用いて、診断支援情報としてのヒストグラム、レーダーチャート等を生成する(ステップS5)。
【0055】
図8は、本ステップS5において生成される診断支援情報としてのヒストグラムの一例を示した図であり、甲、乙、丙の各被検体について、「正常」、「疾患A」、「疾患B」、・・・・の各項目(各疾患)についての相対マハラノビス距離(相対MD値)が示されている。
【0056】
また、より具体的な例として、図9(a)、(b)に、良好な結果が得られたDCMのヒストグラムと、判別困難だったARのヒストグラムを示した。図中の横軸1〜12は、それぞれ健常者、AR、AS−moderete、狭心症等の項目に対応した物差しであり、縦軸は距離を表している。図9(a)は、物差し4に相当するDCM疾患のデータをあてはめたものであり、太線は平均値である。DCMの例では物差し4にピークがあり、良好な識別ができていることがわかる。一方、自9(b)は、物差し2に相当するAR疾患のデータを当てはめたものであり、太字は同じく平均値である。この例では、物差し2よりも物差し5や物差し10にもピークがあり、うまく判別できていないことがわかる。
【0057】
図10は、本ステップS5において生成される診断支援情報としてのレーダーチャート(健常者[1]と、各疾患[2]〜[12]の距離の平均値のレーダーチャート)の一例を示した図である。外周に行くほど各単位空間に近い表現とした。DSMやAS−modereteは良好な結果を示しているが、その他は、別の疾患にもピークが発生し良好ではないことがわかる。
【0058】
また、図11は、本ステップS5において生成される診断支援情報としてのレーダーチャートの他の例を示した図である。比較的良好な疾患判別ができたAS−modereteとDCMの平均値と適合性の良いデータの例を示した。これらの解析に用いた診断指標中、E波速度やE波速度と弁輪速度の比の寄与が大きく、ドプラ診断で得られる情報のみで疾患の判別ができる可能性が高いこともわかる。
【0059】
(実施例)
次に、本自動診断支援処理の実施例について説明する。
【0060】
自動診断支援技術の実用性を高めるため健常者と非健常者の判別以外に、心疾患の種類を判別することを試みた。同一医師の過去数年間に亘診断データをもとにMTシステムを用いてその可能性の検討を行った。
【0061】
すなわち、健常者と非健常者の心エコー図検査の診断指標データ約3000人分を、健常者,AR、AS−moderate、DCM、HCM/DCM、HHD、MR、MS、TR、OMI、先天性心疾患、狭心症、分類不能の13種類に分けた。健常者データ数はジャックナイフ法を用いて単位空間を調整したので最終的には900となり、分類不能を除く11種類の疾患別データは40〜300程度となった。診断指標として、ドプラ計測値以外に、年齢、性別、心拍数、MモードやBモード像から得られた診断指標約20種類を組合せて解析を行った。ここでは波形パターン認識データは用いていない。この様な解析から、健常者の単位空間をもとに、特定疾患の要因効果分析を通して分離性能の評価を行った。
【0062】
図12に、健常者と非健常者(全疾患)のマハラノビス距離(MD値)ヒストグラムを示す。健常者の空間と非健常者の空間が重複しており区別できないことがわかる。波形パターン認識データがなく、単なる計測指標だけでは健常者と非健常者(全疾患)を区別するのが困難なことがわかった。
【0063】
図13は、疾患別の健常者との判別を評価するための各疾患毎(DCM、MR、AS、AR)のヒストグラムを示した図である。同図に示すように、DCRとMRとでは、比較的良好な分離ができた一方、ASとARとでは、分離が困難であることがわかる。
【0064】
図14は、本実施例において生成された、マハラノビス距離(MD値)に関する疾患毎のヒストグラムを示している。として示してある。同図のマハラノビス距離は、上述したように、被検者の診断指標データを、各疾患毎の物差しで正規化して得られたものである。このヒストグラムでは、距離がゼロに近ければ近いほど、それに対応した疾患パターンと一致していることとなる。この様なヒストグラムを診断支援情報として提供することで、医師等のユーザは、距離が近いものを疾患の候補としてピックアップすることができる。
【0065】
(効果)
以上述べた構成によれば、以下の効果を得ることができる。
【0066】
本医用画像診断支援システム或いは医用診断支援装置によれば、超音波診断装置等の医用画像診断装置の自動計測によって得られた自動計測値(診断指標)、画像パターン、波形等を用いてパラメータを抽出し、当該パラメータをMTシステムの項目因子として疾患症例毎に物差しを作成し、所定の患者(すなわち、正常又は異常、或いは異常である場合の疾患種別が未知である患者)に関して、健康状態が正常であるか又は異常であるか、或いは異常である場合の疾患種別を判定し、その結果を診断支援情報として所定の形態で表示する。従って、自動計測値を用いたパラメータをMTシステムの項目因子とした、マハラノビス距離という客観的な定量値に基づいて、多重共線性を回避しつつ高い汎用性をもって、迅速且つ簡便に正常/異常の判別及び疾患種類の判別を実行することができる。特に、パターン認識を用いる従来技術と比較した場合、本自動診断支援装置によれば、標本の項目因子が過去の例に属さない振る舞いをするものであっても、困難なく正常/異常の判別及び疾患種類の判別を実行することができる。
【0067】
また、本医用画像診断支援システム或いは医用診断支援装置によれば、正常/異常の判別結果及び疾患種類の判別結果を、疾患種別の相対的距離を示すレーダーチャートやヒストグラム等によって表示する。従って、医師等のユーザは、疾患種別毎のピークの発生やその形状によって、迅速且つ簡便に検査結果を視認することができる。
【0068】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。具体的な変形例としては、例えば次のようなものがある。
【0069】
(1)本実施形態に係る各機能は、当該処理を実行するプログラムをワークステーション等のコンピュータにインストールし、これらをメモリ上で展開することによっても実現することができる。このとき、コンピュータに当該手法を実行させることのできるプログラムは、磁気ディスク(フロッピー(登録商標)ディスク、ハードディスクなど)、光ディスク(CD−ROM、DVDなど)、半導体メモリなどの記録媒体に格納して頒布することも可能である。
【0070】
(2)上記実施形態においては、自動診断支援システムとして、超音波診断装置2と自動診断支援装置3とを組み合わせる場合を例とした。しかしながら、自動診断支援装置3と組み合わせる医用画像診断装置は、超音波診断装置に拘泥されない。例えば、X線コンピュータ断層撮影装置、磁気共鳴イメージング装置、X線診断装置、核医学診断装置等を採用することも可能である。
【0071】
また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【産業上の利用可能性】
【0072】
以上本発明によれば、MTシステムを用いて疾患症例毎の尺度を作成し、当該尺度と自動計測によって取得された項目とを用いて汎用性の高い自動診断支援を実行可能な自動診断支援装置、超音波診断装置及び自動診断支援プログラムを実現することができる。
【符号の説明】
【0073】
1…自動診断支援システム、2…超音波診断装置、3…自動診断支援装置、20…超音波プローブ、21…送受信ユニット、22…Bモード処理ユニット、23…ドプラ処理ユニット、24…画像・波形生成ユニット、25…画像合成ユニット、26…表示ユニット、30…データベースユニット、32…MTシステム実行ユニット、34…解析ユニット、36…診断支援情報生成ユニット、38…表示ユニット
【特許請求の範囲】
【請求項1】
医用画像診断装置によって取得された所定の被検体に関する医用情報を用いて、当該被検体に関する複数のパラメータを抽出すると共に、予め取得された非健常者に関する医用情報を用いて、非健常者に関する疾患毎の前記複数のパラメータを抽出する抽出手段と、
前記非健常者に関する疾患毎の前記複数のパラメータを項目因子としてMTシステムを実行することで、各疾患症例に関する物差しを計算すると共に、前記被検体に関する複数のパラメータを項目因子としてMTシステムを実行することで、前記被検体の前記各疾患の非健常者の状態空間からの距離を計算する計算手段と、
前記各疾患症例に関する物差しと前記被検体の前記各疾患の非健常者の状態空間からの距離とを用いて、前記被検体の疾患種別を判定する判定手段と、
前記判定の結果に基づいて、診断支援情報を生成する支援情報生成手段と、
前記診断支援情報を所定の形態で表示する表示手段と、
を具備することを特徴とする自動診断支援装置。
【請求項2】
前記抽出手段は、予め取得された健常者に関する医用情報を用いて、健常者に関する前記複数のパラメータを抽出し、
前記計算手段は、前記健常者に関する前記複数のパラメータを項目因子としてMTシステムを実行することで、健常者に関する物差しを計算すると共に、前記被検体に関する複数のパラメータを項目因子としてMTシステムを実行することで、前記被検体の健常者の状態空間からの距離を計算し、
前記判定手段は、前記健常者に関する物差しと前記被検体の健常者の状態空間からの距離とを用いて、前記被検体が健常者か非健常者かを判定し、
前記支援情報生成手段は、前記被検体が健常者か非健常者かの判定結果を含む、前記診断支援情報を生成すること、
を特徴とする請求項1記載の自動診断支援装置。
【請求項3】
前記医用画像診断装置は、超音波診断装置であることを特徴とする請求項1又は2記載の自動診断支援装置。
【請求項4】
前記所定のパラメータは、Bモード画像に関するパターン認識データ、ドプラ波形による波形パターン認識データ、ドプラ自動計測による診断指標、心電図に基づく診断指標のうちの少なくとも一つを含むことを特徴とする請求項1乃至3のうちいずれか一項記載の自動診断支援装置。
【請求項5】
前記計算手段は、前記各疾患に関する物差し及び前記被検体の前記各疾患の非健常者の状態空間からの距離を正規化し、
前記判定手段は、
正規化された前記各疾患症例に関する物差し及び正規化された前記被検体の前記各疾患の非健常者の状態空間からの距離を用いて、疾患判別の相対的距離を計算し、
前記疾患判別の相対的距離に基づいて、前記被検体の疾患種別を判定すること、
を特徴とする請求項1乃至4のうちいずれか一項記載の自動診断支援装置。
【請求項6】
前記支援情報生成手段は、前記診断支援情報として、各疾患を変数するレーダーチャートを生成することを特徴とする請求項1乃至5のうちいずれか一項記載の自動診断支援装置。
【請求項7】
前記支援情報生成手段は、各疾患につき外周に行くほど当該疾患の可能性が高くなり、且つ正常と推定される範囲を参照表示する前記レーダーチャートを生成することを特徴とする請求項6記載の自動診断支援装置。
【請求項8】
前記支援情報生成手段は、前記診断支援情報として、前記各疾患の可能性を示すための定量値又はヒストグラムを生成することを特徴とする請求項1乃至7のうちいずれか一項記載の自動診断支援装置。
【請求項9】
前記医用画像診断装置による前記医用情報の取得の完了に連動して、当該装置の動作タイミングを制御する制御手段をさらに具備することを特徴とする請求項1乃至8のうちいずれか一項記載の自動診断支援装置。
【請求項10】
所定の被検体を超音波走査し、当該被検体に関する超音波データを取得する超音波データ取得手段と、
前記超音波データを含む前記被検体に関する医用情報を用いて、当該被検体に関する複数のパラメータを抽出すると共に、予め取得された非健常者に関する医用情報を用いて、非健常者に関する疾患毎の前記複数のパラメータを抽出する抽出手段と、
前記非健常者に関する疾患毎の前記複数のパラメータを項目因子としてMTシステムを実行することで、各疾患症例に関する物差しを計算すると共に、前記被検体に関する複数のパラメータを項目因子としてMTシステムを実行することで、前記被検体の前記各疾患の非健常者の状態空間からの距離を計算する計算手段と、
前記各疾患症例に関する物差しと前記被検体の前記各疾患の非健常者の状態空間からの距離とを用いて、前記被検体の疾患種別を判定する判定手段と、
前記判定の結果に基づいて、診断支援情報を生成する支援情報生成手段と、
前記診断支援情報を所定の形態で表示する表示手段と、
を具備することを特徴とする超音波診断装置。
【請求項11】
コンピュータに、
医用画像診断装置によって取得された所定の被検体に関する医用情報を用いて、当該被検体に関する複数のパラメータを抽出させると共に、予め取得された非健常者に関する医用情報を用いて、非健常者に関する疾患毎の前記複数のパラメータを抽出させる抽出機能と、
前記非健常者に関する疾患毎の前記複数のパラメータを項目因子としてMTシステムを実行させることで、各疾患症例に関する物差しを計算させると共に、前記被検体に関する複数のパラメータを項目因子としてMTシステムを実行させることで、前記被検体の前記各疾患の非健常者の状態空間からの距離を計算させる計算機能と、
前記各疾患症例に関する物差しと前記被検体の前記各疾患の非健常者の状態空間からの距離とを用いて、前記被検体の疾患種別を判定させる判定機能と、
前記判定の結果に基づいて、診断支援情報を生成させる支援情報生成機能と、
前記診断支援情報を所定の形態で表示させる表示機能と、
を実現させることを特徴とする自動診断支援プログラム。
【請求項1】
医用画像診断装置によって取得された所定の被検体に関する医用情報を用いて、当該被検体に関する複数のパラメータを抽出すると共に、予め取得された非健常者に関する医用情報を用いて、非健常者に関する疾患毎の前記複数のパラメータを抽出する抽出手段と、
前記非健常者に関する疾患毎の前記複数のパラメータを項目因子としてMTシステムを実行することで、各疾患症例に関する物差しを計算すると共に、前記被検体に関する複数のパラメータを項目因子としてMTシステムを実行することで、前記被検体の前記各疾患の非健常者の状態空間からの距離を計算する計算手段と、
前記各疾患症例に関する物差しと前記被検体の前記各疾患の非健常者の状態空間からの距離とを用いて、前記被検体の疾患種別を判定する判定手段と、
前記判定の結果に基づいて、診断支援情報を生成する支援情報生成手段と、
前記診断支援情報を所定の形態で表示する表示手段と、
を具備することを特徴とする自動診断支援装置。
【請求項2】
前記抽出手段は、予め取得された健常者に関する医用情報を用いて、健常者に関する前記複数のパラメータを抽出し、
前記計算手段は、前記健常者に関する前記複数のパラメータを項目因子としてMTシステムを実行することで、健常者に関する物差しを計算すると共に、前記被検体に関する複数のパラメータを項目因子としてMTシステムを実行することで、前記被検体の健常者の状態空間からの距離を計算し、
前記判定手段は、前記健常者に関する物差しと前記被検体の健常者の状態空間からの距離とを用いて、前記被検体が健常者か非健常者かを判定し、
前記支援情報生成手段は、前記被検体が健常者か非健常者かの判定結果を含む、前記診断支援情報を生成すること、
を特徴とする請求項1記載の自動診断支援装置。
【請求項3】
前記医用画像診断装置は、超音波診断装置であることを特徴とする請求項1又は2記載の自動診断支援装置。
【請求項4】
前記所定のパラメータは、Bモード画像に関するパターン認識データ、ドプラ波形による波形パターン認識データ、ドプラ自動計測による診断指標、心電図に基づく診断指標のうちの少なくとも一つを含むことを特徴とする請求項1乃至3のうちいずれか一項記載の自動診断支援装置。
【請求項5】
前記計算手段は、前記各疾患に関する物差し及び前記被検体の前記各疾患の非健常者の状態空間からの距離を正規化し、
前記判定手段は、
正規化された前記各疾患症例に関する物差し及び正規化された前記被検体の前記各疾患の非健常者の状態空間からの距離を用いて、疾患判別の相対的距離を計算し、
前記疾患判別の相対的距離に基づいて、前記被検体の疾患種別を判定すること、
を特徴とする請求項1乃至4のうちいずれか一項記載の自動診断支援装置。
【請求項6】
前記支援情報生成手段は、前記診断支援情報として、各疾患を変数するレーダーチャートを生成することを特徴とする請求項1乃至5のうちいずれか一項記載の自動診断支援装置。
【請求項7】
前記支援情報生成手段は、各疾患につき外周に行くほど当該疾患の可能性が高くなり、且つ正常と推定される範囲を参照表示する前記レーダーチャートを生成することを特徴とする請求項6記載の自動診断支援装置。
【請求項8】
前記支援情報生成手段は、前記診断支援情報として、前記各疾患の可能性を示すための定量値又はヒストグラムを生成することを特徴とする請求項1乃至7のうちいずれか一項記載の自動診断支援装置。
【請求項9】
前記医用画像診断装置による前記医用情報の取得の完了に連動して、当該装置の動作タイミングを制御する制御手段をさらに具備することを特徴とする請求項1乃至8のうちいずれか一項記載の自動診断支援装置。
【請求項10】
所定の被検体を超音波走査し、当該被検体に関する超音波データを取得する超音波データ取得手段と、
前記超音波データを含む前記被検体に関する医用情報を用いて、当該被検体に関する複数のパラメータを抽出すると共に、予め取得された非健常者に関する医用情報を用いて、非健常者に関する疾患毎の前記複数のパラメータを抽出する抽出手段と、
前記非健常者に関する疾患毎の前記複数のパラメータを項目因子としてMTシステムを実行することで、各疾患症例に関する物差しを計算すると共に、前記被検体に関する複数のパラメータを項目因子としてMTシステムを実行することで、前記被検体の前記各疾患の非健常者の状態空間からの距離を計算する計算手段と、
前記各疾患症例に関する物差しと前記被検体の前記各疾患の非健常者の状態空間からの距離とを用いて、前記被検体の疾患種別を判定する判定手段と、
前記判定の結果に基づいて、診断支援情報を生成する支援情報生成手段と、
前記診断支援情報を所定の形態で表示する表示手段と、
を具備することを特徴とする超音波診断装置。
【請求項11】
コンピュータに、
医用画像診断装置によって取得された所定の被検体に関する医用情報を用いて、当該被検体に関する複数のパラメータを抽出させると共に、予め取得された非健常者に関する医用情報を用いて、非健常者に関する疾患毎の前記複数のパラメータを抽出させる抽出機能と、
前記非健常者に関する疾患毎の前記複数のパラメータを項目因子としてMTシステムを実行させることで、各疾患症例に関する物差しを計算させると共に、前記被検体に関する複数のパラメータを項目因子としてMTシステムを実行させることで、前記被検体の前記各疾患の非健常者の状態空間からの距離を計算させる計算機能と、
前記各疾患症例に関する物差しと前記被検体の前記各疾患の非健常者の状態空間からの距離とを用いて、前記被検体の疾患種別を判定させる判定機能と、
前記判定の結果に基づいて、診断支援情報を生成させる支援情報生成機能と、
前記診断支援情報を所定の形態で表示させる表示機能と、
を実現させることを特徴とする自動診断支援プログラム。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【公開番号】特開2011−10775(P2011−10775A)
【公開日】平成23年1月20日(2011.1.20)
【国際特許分類】
【出願番号】特願2009−156318(P2009−156318)
【出願日】平成21年6月30日(2009.6.30)
【出願人】(000003078)株式会社東芝 (54,554)
【出願人】(594164542)東芝メディカルシステムズ株式会社 (4,066)
【Fターム(参考)】
【公開日】平成23年1月20日(2011.1.20)
【国際特許分類】
【出願日】平成21年6月30日(2009.6.30)
【出願人】(000003078)株式会社東芝 (54,554)
【出願人】(594164542)東芝メディカルシステムズ株式会社 (4,066)
【Fターム(参考)】
[ Back to top ]