(57)【要約】【課題】  画像データ（特に、血管内超音波画像データ）内における組織性状（例えば、繊維性組織、脂質性組織）をより正確に判別する組織性状判別方法、プログラムおよび装置を提供する。【解決手段】  画像データ（超音波画像）の注目点近傍の局所画像情報（超音波ＲＦ信号）から、前記注目点の組織性状を判別する組織性状判別方法（プログラム、装置）であって、前記注目点近傍の局所画像情報（超音波ＲＦ信号）を、予めスパースコーディング法による学習（ステップ５１，５２，５３，５４）で得られた組織性状ごとの基底関数群により分解し、基底関数ごとの係数を算出する係数算出ステップ（ステップ５５）と、前記係数算出ステップにより算出された基底関数ごとの係数に基づいて、前記注目点の組織性状を判別する組織性状判別ステップ（ステップ５６）と、を有する。【選択図】図５
2010075616.tif 000002
【特許請求の範囲】
【請求項１】
超音波画像の注目点近傍の超音波ＲＦ信号から、前記注目点の組織性状を判別する組織性状判別方法であって、
前記注目点近傍の超音波ＲＦ信号を、予めスパースコーディング法による学習で得られた組織性状ごとの基底関数群により分解し、基底関数ごとの係数を算出する係数算出ステップと、
前記係数算出ステップにより算出された基底関数ごとの係数に基づいて、前記注目点の組織性状を判別する組織性状判別ステップと、を有する
組織性状判別方法。
【請求項２】
前記超音波画像内に組織性状ごとに学習領域を設定し、前記学習領域内の複数の超音波ＲＦ信号を用いて、スパースコーディング法による学習により前記組織性状ごとの基底関数群を算出する基底関数群学習ステップ、をさらに有する
請求項１記載の組織性状判別方法。
【請求項３】
前記基底関数群学習ステップにより算出された組織性状ごとの基底関数群を用いて、前記基底関数群を算出したのと同じ組織性状の前記学習領域内の複数の超音波ＲＦ信号を分解し、基底関数ごとの係数値を累積することにより第１係数累積値群を求める第１係数累積値群算出ステップと、
前記基底関数群学習ステップにより算出された組織性状ごとの基底関数群を用いて、前記基底関数群を算出したのと異なる組織性状の前記学習領域内の複数の超音波ＲＦ信号を分解し、基底関数ごとの係数値を累積することにより第２係数累積値群を求める第２係数累積値群算出ステップと、をさらに有し、
前記組織性状判別ステップは、前記第１係数累積値群および前記第２係数累積値群を用いて前記注目点の組織性状を判別する、
請求項２記載の組織性状判別方法。
【請求項４】
前記第１係数累積値群および前記第２係数累積値群を２値化するステップ、をさらに有し、
前記組織性状判別ステップは、２値化された前記第１係数累積値群および前記第２係数累積値群を用いて前記注目点の組織性状を判別する、
請求項３記載の組織性状判別方法。
【請求項５】
前記超音波画像は、血管内超音波診断装置によるＢモード画像であり、
前記組織性状は、少なくとも繊維性組織または脂質性組織を含む、
請求項１〜４いずれか記載の組織性状判別方法。
【請求項６】
超音波画像の注目点近傍の超音波ＲＦ信号から、前記注目点の組織性状を判別する組織性状判別プログラムであって、
前記注目点近傍の超音波ＲＦ信号を、予めスパースコーディング法による学習で得られた組織性状ごとの基底関数群により分解し、基底関数ごとの係数を算出する係数算出ステップと、
前記係数算出ステップにより算出された基底関数ごとの係数に基づいて、前記注目点の組織性状を判別する組織性状判別ステップと、を有する
組織性状判別プログラム。
【請求項７】
超音波画像の注目点近傍の超音波ＲＦ信号から、前記注目点の組織性状を判別する組織性状判別装置であって、
予めスパースコーディング法による学習で得られた組織性状ごとの基底関数群を記憶している記憶手段と、
前記注目点の組織性状を判別する演算手段と、を有し、
前記演算手段は、
前記注目点近傍の超音波ＲＦ信号を、前記記憶手段に記憶されている組織性状ごとの基底関数群により分解し、基底関数ごとの係数を算出する係数算出手段と、
前記係数算出手段により算出された基底関数ごとの係数に基づいて、前記注目点の組織性状を判別する組織性状判別手段と、を有する
組織性状判別装置。
【請求項８】
画像データの注目点近傍の局所画像情報から、前記注目点の組織性状を判別する組織性状判別方法であって、
前記注目点近傍の局所画像情報を、予めスパースコーディング法による学習で得られた組織性状ごとの基底関数群により分解し、基底関数ごとの係数を算出する係数算出ステップと、
前記係数算出ステップにより算出された基底関数ごとの係数と、組織性状ごとに予め算出された第１係数累積値群および第２係数累積値群とに基づいて、前記注目点の組織性状を判別する組織性状判別ステップと、を有し、
前記組織性状ごとの基底関数群は、前記画像データ内に組織性状ごとに学習領域を設定し、前記学習領域内の局所画像情報を用いて、スパースコーディング法による学習により前記組織性状ごとの基底関数群を算出する基底関数群学習ステップにより算出され、
前記第１係数累積値群は、前記基底関数群学習ステップにより算出された組織性状ごとの基底関数群を用いて、前記基底関数群を算出したのと同じ組織性状の前記学習領域内の局所画像情報を分解し、基底関数ごとの係数値を累積することにより第１係数累積値を算出する第１係数累積値群算出ステップにより算出され、
前記第２係数累積値群は、前記基底関数群学習ステップにより算出された組織性状ごとの基底関数群を用いて、前記基底関数群を算出したのと異なる組織性状の前記学習領域内の局所画像情報を分解し、基底関数ごとの係数値を累積することにより第２係数累積値を算出する第２係数累積値群算出ステップにより算出される、
組織性状判別方法。
【請求項９】
画像データの注目点近傍の局所画像情報から、前記注目点の組織性状を判別する組織性状判別プログラムであって、
前記注目点近傍の局所画像情報を、予めスパースコーディング法による学習で得られた組織性状ごとの基底関数群により分解し、基底関数ごとの係数を算出する係数算出ステップと、
前記係数算出ステップにより算出された基底関数ごとの係数と、組織性状ごとに予め算出された第１係数累積値群および第２係数累積値群とに基づいて、前記注目点の組織性状を判別する組織性状判別ステップと、を有し、
前記組織性状ごとの基底関数群は、前記画像データ内に組織性状ごとに学習領域を設定し、前記学習領域内の局所画像情報を用いて、スパースコーディング法による学習により前記組織性状ごとの基底関数群を算出する基底関数群学習ステップにより算出され、
前記第１係数累積値群は、前記基底関数群学習ステップにより算出された組織性状ごとの基底関数群を用いて、前記基底関数群を算出したのと同じ組織性状の前記学習領域内の局所画像情報を分解し、基底関数ごとの係数値を累積することにより第１係数累積値を算出する第１係数累積値群算出ステップにより算出され、
前記第２係数累積値群は、前記基底関数群学習ステップにより算出された組織性状ごとの基底関数群を用いて、前記基底関数群を算出したのと異なる組織性状の前記学習領域内の局所画像情報を分解し、基底関数ごとの係数値を累積することにより第２係数累積値を算出する第２係数累積値群算出ステップにより算出される、
組織性状判別プログラム。
【請求項１０】
画像データの注目点近傍の局所画像情報から、前記注目点の組織性状を判別する組織性状判別装置であって、
予めスパースコーディング法による学習で得られた組織性状ごとの基底関数群と、第１係数累積値群および第２係数累積値群とを記憶している記憶手段と、
前記注目点の組織性状を判別する演算手段と、を有し、
前記演算手段は、
前記注目点近傍の局所画像情報を、前記記憶手段に記憶されている組織性状ごとの基底関数群により分解し、基底関数ごとの係数を算出する係数算出手段と、
前記係数算出手段により算出された基底関数ごとの係数と、前記記憶手段に記憶されている組織性状ごとに予め算出された第１係数累積値群および第２係数累積値群とに基づいて、前記注目点の組織性状を判別する組織性状判別手段と、を有し、
前記組織性状ごとの基底関数群は、前記画像データ内に組織性状ごとに学習領域を設定し、前記学習領域内の局所画像情報を用いて、スパースコーディング法による学習により算出されたものであり、
前記第１係数累積値群は、前記組織性状ごとの基底関数群を用いて、前記基底関数群を算出したのと同じ組織性状の前記学習領域内の局所画像情報を分解し、基底関数ごとの係数値を累積することにより算出されたものであり、
前記第２係数累積値群は、前記組織性状ごとの基底関数群を用いて、前記基底関数群を算出したのと異なる組織性状の前記学習領域内の局所画像情報を分解し、基底関数ごとの係数値を累積することにより算出されたものである、
組織性状判別装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、スパースコーディングを用いた組織性状判別に関する。特に、血管内超音波診断装置（Intravascular
Ultrasound：IVUS）による血管超音波データから、血管内の組織性状（例えば、繊維性組織、脂質性組織など）を判別する技術に関する。
【背景技術】
【０００２】
近年、急性心筋梗塞や不安定狭心症などの急性冠症候群が増加している。急性冠症候群は、血管壁に付着した脂肪成分の固まりであるプラークが破綻し、血栓が生じることによって引き起こされる。そのため、心臓に酸素や栄養を送っている冠動脈の内腔が狭くなり、血流が悪くなる病気である。
【０００３】
急性冠症候群の検査には、血管内超音波エコー法が用いられる。血管内超音波エコー法とは、血管内にカテーテルと呼ばれる細長い管を挿入し、その先端に装着された超音波探触子（プローブ）から超音波を組織に照射し、組織からの反射ＲＦ信号を解析することで、血管内部の状態を調べるという検査方法である。ＲＦ信号の振幅値を輝度値に変換することで、血管断面を画像として表示することができる。図１に血管内超音波エコー法の概略図を示す。例えば、超音波の送信周波数は４０ＭＨｚであり、受信した超音波信号（ＲＦ信号）のサンプリング周波数は４００ＭＨｚである。
【０００４】
血管内超音波エコー法には、血管内部の状態をリアルタイムに表示可能、送波信号のエネルギーがＸ線に比べて小さいので検査対象に対する負担が少ないという長所がある。また、血管内超音波エコー法によって得られるＲＦ信号は高周波であるので，血管内部構造を詳細に表示することができる。しかし、出力画像の解像度がよくないこと、減衰が激しく到達距離が短いこと、赤血球などによって拡散され、スペックルと呼ばれるノイズが生じることなどの短所もある。このような性質をもつ画像に対して、実際の医療現場では、熟練した医師が経験をもとに組織部分を目視で診断している。組織部分を正確に見極めることは、病気の早期発見や進行度合いなどを診断する上で非常に重要である。
【０００５】
冠動脈プラークの破綻は急性冠症候群の主な原因であるが、このプラークには破綻しやすいプラーク（不安定プラーク）と破綻しにくいプラーク（安定プラーク）がある。急性冠症候群の診断ではプラークの組織性状判別を行い、プラークの安定度を評価する必要がある。しかしながら、現在の血管内超音波エコー法は、血管内部の形状、プラークサイズ、冠動脈の狭窄程度の判断はできるが、プラークの組織性状判別については適当な解析方法が開発されていない。
【０００６】
以下、冠動脈プラークについて説明する。
心筋（心臓の筋肉）に血液を供給する血管を冠動脈という。冠動脈は直径３〜４mmの大きさで、内膜、中膜、外膜の三層構造から成り立っている（図２（ａ）参照）。これを総称して血管壁という。内膜と中膜は内弾性板により、中膜と外膜は外弾性板と隔てられている。弾性板は弾性線維から構成され、動脈の弾性力や強靭さを維持している。
【０００７】
狭心症や心筋梗塞などの冠症候群を引き起こす原因となる動脈硬化の一つにアテローム（粥状）硬化がある。これは、大動脈や脳動脈、冠動脈などの動脈において、図２（ｂ）のように動脈の内側にコレステロールなどの脂肪からなる粥状物質がたまることによりアテロームプラーク（粥状硬化斑）ができ、次第に肥大することで動脈の内腔を狭める。
【０００８】
アテロームプラークの脂質成分に富む部分を脂質コアと呼ぶ。図３（ａ）のように、この脂質コアが線維性組織によってしっかり囲まれ安定している動脈硬化巣を安定プラークという。この安定プラークが破綻することは少なく、狭心症の不安定化や心筋梗塞にはなりにくい。
【０００９】
一方、図３（ｂ）のように薄い線維性被膜に覆われた脂質コアが主体の場合は不安定プラークと呼ばれる。不安定プラークは図３（ｃ）のような薄い線維性被膜の破裂や亀裂により、図３（ｄ）のようにプラークが破綻して血栓を生じ、急性冠症候群を起こしやすい。
【００１０】
急性冠症候群にはプラークが破綻し心筋の血流が極端に低下した不安定狭心症と、血流の途絶により心筋の一部が壊死に陥る急性心筋梗塞がある。急性冠症候群の検査には内腔の狭窄度だけでなく、血管壁のアテロームプラークの組織性状の診断も重要となる。
【００１１】
一般的な信号解析方法として、フーリエ解析やウェーブレット解析が知られている。しかしながら、これらの信号解析方法は基底関数を固定しているため、多くの基底関数を用いなければ信号を表わせない欠点がある。また、信号から基底関数を決定する手法として、代表的なものに主成分分析（PCA）がある。PCAは、信号の２次統計量を用いて基底関数を決定するため、ある程度信号の統計的性質を考慮した基底関数を得ることが可能である。しかしながら、基底関数に直交条件を付加しているため、信号の表現には多くの基底関数が必要となる。
【００１２】
本発明に関連する先行技術文献として、特許文献１〜４がある。
特許文献１には、画像の特徴を局所的に基底関数で展開し、画像の認識を行う視覚認識システムが記載されている。しかしながら、特許文献１には、超音波ＲＦ信号を基底関数により展開する点、および、スパースコーディングを用いて組織性状ごとに基底関数群を学習することは記載も示唆もされていない。
特許文献２には、人体の心臓や胸部のモデルを基底体（基底関数）で定義して表現する数学的モデルが記載されている。しかしながら、特許文献２には、超音波ＲＦ信号を基底関数により展開する点、および、スパースコーディングを用いて組織性状ごとに基底関数群を学習することは記載も示唆もされていない。
特許文献３には、コスト関数が最小になるように基底関数群を算出し、該基底関数を用いて物体を認識するする物体認識方法が記載されている。しかしながら、特許文献３には、超音波ＲＦ信号を基底関数により展開する点、および、スパースコーディングを用いて組織性状ごとに基底関数群を学習することは記載も示唆もされていない。
特許文献４には、ＩＣメモリ試験結果の不良ビットマップを二次元ウェーブレット変換することにより、不良ビットのアドレスと不良数を求める不良解析方法および装置が記載されている。しかしながら、特許文献４には、超音波ＲＦ信号を基底関数により展開する点、および、スパースコーディングを用いて組織性状ごとに基底関数群を学習することは記載も示唆もされていない。
【特許文献１】特開２００２−３４２７５８号公報「視覚認識システム」
【特許文献２】特開２００１−２２２７０５号公報「数学的モデル並びに該モデルを用いる方法及び装置」
【特許文献３】特開２０００−３００６７号公報「物体認識方法」
【特許文献４】特開平１１−３０６７９３号公報「不良解析方法および装置」
【発明の開示】
【発明が解決しようとする課題】
【００１３】
画像データの信号解析には従来はフーリエ解析やウェーブレット解析などの信号解析が用いられてきたが、基底関数が固定であるため柔軟性に欠け、画像データから組織性状を判別するには不十分であった。本発明は、画像データ（特に、血管内超音波画像データ）内における組織性状（例えば、繊維性組織、脂質性組織）をより正確に判別する組織性状判別方法、プログラムおよび装置を提供することを目的とする。
【課題を解決するための手段】
【００１４】
上記目的を達成するため、本発明は以下の構成を有する。
超音波画像の注目点近傍の超音波ＲＦ信号から、前記注目点の組織性状を判別する組織性状判別方法であって、
前記注目点近傍の超音波ＲＦ信号を、予めスパースコーディング法による学習で得られた組織性状ごとの基底関数群により分解し、基底関数ごとの係数を算出する係数算出ステップと、
前記係数算出ステップにより算出された基底関数ごとの係数に基づいて、前記注目点の組織性状を判別する組織性状判別ステップと、を有する
組織性状判別方法。
超音波画像の注目点近傍の超音波ＲＦ信号から、前記注目点の組織性状を判別する組織性状判別プログラムであって、
前記注目点近傍の超音波ＲＦ信号を、予めスパースコーディング法による学習で得られた組織性状ごとの基底関数群により分解し、基底関数ごとの係数を算出する係数算出ステップと、
前記係数算出ステップにより算出された基底関数ごとの係数に基づいて、前記注目点の組織性状を判別する組織性状判別ステップと、を有する
組織性状判別プログラム。
超音波画像の注目点近傍の超音波ＲＦ信号から、前記注目点の組織性状を判別する組織性状判別装置であって、
予めスパースコーディング法による学習で得られた組織性状ごとの基底関数群を記憶している記憶手段と、
前記注目点の組織性状を判別する演算手段と、を有し、
前記演算手段は、
前記注目点近傍の超音波ＲＦ信号を、前記記憶手段に記憶されている組織性状ごとの基底関数群により分解し、基底関数ごとの係数を算出する係数算出手段と、
前記係数算出手段により算出された基底関数ごとの係数に基づいて、前記注目点の組織性状を判別する組織性状判別手段と、を有する
組織性状判別装置。
【００１５】
この発明の構成によれば、スパースコーディング法を用いて組織性状ごとに基底関数群を学習しておき、該基底関数群を用いて組織性状を判別するので、従来の直交固定型の基底関数群を用いた信号解析に比べて、より正確な組織性状の判別が可能である。つまり、直交固定型の基底関数群ではないので、組織性状ごとに異なる最適な基底関数群を学習により設定することができる。
また、組織性状の判別や基底関数群の学習などの信号解析に超音波ＲＦ信号を用いているため、Ｂモード化された超音波画像データに比べてより多くの情報を含んでおり、より正確な信号解析が可能である。
超音波画像としては血管内超音波画像が好ましいが、これに限定されず、通常の体外プローブによって得られた超音波画像でも良い。また、組織性状の判別としては、血管壁内の繊維性組織や脂質性組織の判別が好ましいがこれに限定されず、他の臓器や組織の判別に用いても良い。
【００１６】
また、本発明は以下の構成を有する。
画像データの注目点近傍の局所画像情報から、前記注目点の組織性状を判別する組織性状判別方法であって、
前記注目点近傍の局所画像情報を、予めスパースコーディング法による学習で得られた組織性状ごとの基底関数群により分解し、基底関数ごとの係数を算出する係数算出ステップと、
前記係数算出ステップにより算出された基底関数ごとの係数と、組織性状ごとに予め算出された第１係数累積値群および第２係数累積値群とに基づいて、前記注目点の組織性状を判別する組織性状判別ステップと、を有し、
前記組織性状ごとの基底関数群は、前記画像データ内に組織性状ごとに学習領域を設定し、前記学習領域内の局所画像情報を用いて、スパースコーディング法による学習により前記組織性状ごとの基底関数群を算出する基底関数群学習ステップにより算出され、
前記第１係数累積値群は、前記基底関数群学習ステップにより算出された組織性状ごとの基底関数群を用いて、前記基底関数群を算出したのと同じ組織性状の前記学習領域内の局所画像情報を分解し、基底関数ごとの係数値を累積することにより第１係数累積値を算出する第１係数累積値群算出ステップにより算出され、
前記第２係数累積値群は、前記基底関数群学習ステップにより算出された組織性状ごとの基底関数群を用いて、前記基底関数群を算出したのと異なる組織性状の前記学習領域内の局所画像情報を分解し、基底関数ごとの係数値を累積することにより第２係数累積値を算出する第２係数累積値群算出ステップにより算出される、
組織性状判別方法。
画像データの注目点近傍の局所画像情報から、前記注目点の組織性状を判別する組織性状判別プログラムであって、
前記注目点近傍の局所画像情報を、予めスパースコーディング法による学習で得られた組織性状ごとの基底関数群により分解し、基底関数ごとの係数を算出する係数算出ステップと、
前記係数算出ステップにより算出された基底関数ごとの係数と、組織性状ごとに予め算出された第１係数累積値群および第２係数累積値群とに基づいて、前記注目点の組織性状を判別する組織性状判別ステップと、を有し、
前記組織性状ごとの基底関数群は、前記画像データ内に組織性状ごとに学習領域を設定し、前記学習領域内の局所画像情報を用いて、スパースコーディング法による学習により前記組織性状ごとの基底関数群を算出する基底関数群学習ステップにより算出され、
前記第１係数累積値群は、前記基底関数群学習ステップにより算出された組織性状ごとの基底関数群を用いて、前記基底関数群を算出したのと同じ組織性状の前記学習領域内の局所画像情報を分解し、基底関数ごとの係数値を累積することにより第１係数累積値を算出する第１係数累積値群算出ステップにより算出され、
前記第２係数累積値群は、前記基底関数群学習ステップにより算出された組織性状ごとの基底関数群を用いて、前記基底関数群を算出したのと異なる組織性状の前記学習領域内の局所画像情報を分解し、基底関数ごとの係数値を累積することにより第２係数累積値を算出する第２係数累積値群算出ステップにより算出される、
組織性状判別プログラム。
画像データの注目点近傍の局所画像情報から、前記注目点の組織性状を判別する組織性状判別装置であって、
予めスパースコーディング法による学習で得られた組織性状ごとの基底関数群と、第１係数累積値群および第２係数累積値群とを記憶している記憶手段と、
前記注目点の組織性状を判別する演算手段と、を有し、
前記演算手段は、
前記注目点近傍の局所画像情報を、前記記憶手段に記憶されている組織性状ごとの基底関数群により分解し、基底関数ごとの係数を算出する係数算出手段と、
前記係数算出手段により算出された基底関数ごとの係数と、前記記憶手段に記憶されている組織性状ごとに予め算出された第１係数累積値群および第２係数累積値群とに基づいて、前記注目点の組織性状を判別する組織性状判別手段と、を有し、
前記組織性状ごとの基底関数群は、前記画像データ内に組織性状ごとに学習領域を設定し、前記学習領域内の局所画像情報を用いて、スパースコーディング法による学習により算出されたものであり、
前記第１係数累積値群は、前記組織性状ごとの基底関数群を用いて、前記基底関数群を算出したのと同じ組織性状の前記学習領域内の局所画像情報を分解し、基底関数ごとの係数値を累積することにより算出されたものであり、
前記第２係数累積値群は、前記組織性状ごとの基底関数群を用いて、前記基底関数群を算出したのと異なる組織性状の前記学習領域内の局所画像情報を分解し、基底関数ごとの係数値を累積することにより算出されたものである、
組織性状判別装置。
【００１７】
この発明の構成によれば、スパースコーディング法を用いて組織性状ごとに基底関数群を学習しておき、該基底関数群を用いて組織性状を判別するので、従来の直交固定型の基底関数群を用いた信号解析に比べて、より正確な組織性状の判別が可能である。つまり、直交固定型の基底関数群ではないので、組織性状ごとに異なる最適な基底関数群を学習により設定することができる。
また、組織性状ごとに係数累積値群も予め学習により求めている。特に、『第２係数累積値群』は、学習により求められた基底関数群とは異なる組織性状の学習領域内の情報に基づいて算出されている。この第２係数累積値群を用いることにより、より正確な組織性状の判別ができる。３種類以上の組織性状を判別する場合には、第２係数累積値群は異なる組織性状の種類ごとに算出される。
画像データとしては、超音波画像に限られず、レントゲン画像、ＣＴ画像、ＭＲＩ画像などの医療用画像でも良いし、より一般的な画像の認識に用いても良い。局所画像情報は、局所的な画像に関する情報であれば良く、１次元の情報でも良いし、２次元のビットマップ情報でも良い。
【００１８】
また、好ましい実施態様として以下のものがある。
前記超音波画像内に組織性状ごとに学習領域を設定し、前記学習領域内の複数の超音波ＲＦ信号を用いて、スパースコーディング法による学習により前記組織性状ごとの基底関数群を算出する基底関数群学習ステップ、をさらに有する。
【００１９】
また、好ましい実施態様として以下のものがある。
前記基底関数群学習ステップにより算出された組織性状ごとの基底関数群を用いて、前記基底関数群を算出したのと同じ組織性状の前記学習領域内の複数の超音波ＲＦ信号を分解し、基底関数ごとの係数値を累積することにより第１係数累積値群を求める第１係数累積値群算出ステップと、
前記基底関数群学習ステップにより算出された組織性状ごとの基底関数群を用いて、前記基底関数群を算出したのと異なる組織性状の前記学習領域内の複数の超音波ＲＦ信号を分解し、基底関数ごとの係数値を累積することにより第２係数累積値群を求める第２係数累積値群算出ステップと、をさらに有し、
前記組織性状判別ステップは、前記第１係数累積値群および前記第２係数累積値群を用いて前記注目点の組織性状を判別する。
組織性状ごとに係数累積値群も予め学習により求めている。特に、『第２係数累積値群』は、学習により求められた基底関数群とは異なる組織性状の学習領域内の情報に基づいて算出されている。この第２係数累積値群を用いることにより正確な組織性状の判別ができる。３種類以上の組織性状を判別する場合には、第２係数累積値群は異なる組織性状の種類ごとに算出される。
【００２０】
また、好ましい実施態様として以下のものがある。
前記第１係数累積値群および前記第２係数累積値群を２値化するステップ、をさらに有し、
前記組織性状判別ステップは、２値化された前記第１係数累積値群および前記第２係数累積値群を用いて前記注目点の組織性状を判別する。
例えば、第１係数累積値群および第２係数累積値群を、平均値の上下で２値化する。２値化することにより、演算を単純化することができ、組織性状の判別の演算が容易になる。
【００２１】
また、好ましい実施態様として以下のものがある。
前記超音波画像は、血管内超音波診断装置によるＢモード画像であり、
前記組織性状は、少なくとも繊維性組織または脂質性組織を含む。
判別する組織性状は特に限定されないが、血管内の繊維性組織および脂質性組織の判別を行うことが特に好ましい。
【発明の効果】
【００２２】
本発明の組織性状判別方法、プログラムおよび装置によれば、画像データ（特に、血管内超音波画像データ）内における組織性状（例えば、繊維性組織、脂質性組織）をより正確に判別することができる。
【発明を実施するための最良の形態】
【００２３】
以下、本発明の実施形態の一例について説明する。以下の演算やプロセスは、ＰＣ等で実行されても良いし、超音波診断装置内の演算装置により実行されても良い。
【００２４】
＜＜スパースコーディングについて＞＞
まずは、スパースコーディングについて説明する。スパースコーディングとは、重みに疎（スパース性）の条件を付加することであり、信号を少数の基底関数の線形和で効率的に表現することである。得られる基底関数は、信号の統計的性質をよく反映しているという特徴がある。
【００２５】
Ｎ次元の入力信号Ｘj（j=1,2,…,N）は基底関数φijと係数ａiの線形和
【数１】
2010075616.tif 000003
で表現されると仮定する。ただし、Ｍは基底関数の数である。再構成二乗誤差と出力信号のスパース性を考慮したコスト関数Ｅを以下の式で定義する。
【数２】
2010075616.tif 000004
ここで、第１項は誤差項、第２項はスパース項である。また、βは正の定数であり、σは入力データの標準偏差である。誤差項は、入力信号Ｘjと再構成した信号の平均二乗誤差を表し、入力信号がどれだけ復元されているかを表す。スパース項は出力信号のスパース性の程度を示すものである。Ｓ（x）はx＝０で最小になり、xが少量変化すると値が急激に増加するような関数を用いる。例えば、log（１＋x2）、−ｅ-x^2などが考えられる。また、急激に値が増加する関数ではないが|x|（絶対値x）などを用いても良い。
【００２６】
ａiとφijは学習によって求められる。基底関数の学習手順を、図４に示す。
＜ステップ４２＞新たな入力信号Ｘjを入力。
＜ステップ４３＞信号が入力されるごとに、以下の更新式でａiが更新される。
【数３】
2010075616.tif 000005
ここで、
【数４】
2010075616.tif 000006
である。なお、式（数３）は、式（数２）のＥをａiについて偏微分することで求められる。
＜ステップ４４＞ステップ４３を数回繰り返した後、一定数の信号に対して以下の更新式によりφijを更新する。
【数５】
2010075616.tif 000007
式（数５）は、式（数２）のＥをφijについて偏微分することで求められる。ここで、ηは学習係数であり、小さな正の定数である。また、<>は平均値を表す。
【００２７】
＜＜組織性状判別の全体フロー＞＞
本実施形態における全体フローについて説明する。図５に、本実施形態における組織性状の判別の全体フロー図を示す。
＜ステップ５１＞組織性状（繊維性組織、脂質性組織）ごとに基底関数群を学習するための学習領域を設定する。
＜ステップ５２＞学習領域内の複数の注目点近傍の超音波ＲＦ信号を用いて、スパースコーディングにより組織性状ごとに基底関数群を学習する。
＜ステップ５３＞ステップ５２で求められた基底関数群を用いて、学習領域内の注目点近傍の超音波ＲＦ信号を分解し、それぞれの基底関数ごとの係数を算出する。学習領域内の複数の注目点ごとに係数を算出し、それらの絶対値を累積することにより係数累積値群を算出する。
＜ステップ５４＞ステップ５３で求められた係数累積値群を平均値の上下で２値化する。
＜ステップ５５＞組織性状を判別したい注目点近傍の超音波ＲＦ信号（未知信号）を、ステップ５２で算出された基底関数群で分解し、基底関数ごとの係数を算出する。
＜ステップ５６＞ステップ５５で求められた係数を、ステップ５４で求められた２値化された係数累積値群を用いて未知信号がどの組織性状に属するかを判別する。
以下で、各プロセスについて詳細に説明する。
【００２８】
＜＜組織性状ごとの基底関数の学習プロセス＞＞
図５における、ステップ５１およびステップ５２の学習プロセスについて説明する。
図６に、注目点近傍の超音波ＲＦ信号の説明図を示す。図６（ａ）は血管超音波Ｂモード画像における注目点を示している。注目点を通る矢印は、超音波の送出方向を示している。血管超音波像なので、超音波は血管内にあるカテーテルにより送受信される。したがって、注目点近傍における超音波ＲＦ信号は図６（ｂ）のようになる。本実施形態では、この注目点近傍の超音波ＲＦ信号を、基底関数の学習や、組織性状の判別に用いる。
【００２９】
図７に繊維性組織の基底関数の学習の説明図、図８に脂質性組織の基底関数の学習の説明図を示す。図７および図８に示すように、組織性状（繊維性組織、脂質性組織）ごとの学習領域を設定する（ステップ５１）。これらの学習領域は、各種所見に基づき既に組織性状が判別されている領域を設定する。続いて、学習領域内の複数の注目点近傍の超音波ＲＦ信号を学習データＸjとし、これらの複数の学習データＸjを用いて上述のスパースコーディングを用いて基底関数φijを学習する（ステップ５２）。Ｘjは注目点近傍の超音波ＲＦ信号を時間方向に離散化したデータであり、例えば、３２点（j=1,2,…,32）のデータである。また基底関数の数は任意であるが、本実施形態では３２個の基底関数を設定している。このようにして、繊維性組織の基底関数群φfibijおよび脂質性組織の基底関数群φlipijが学習される。
【００３０】
＜＜組織性状判別のためのフィルタ作成プロセス＞＞
図５における、ステップ５３およびステップ５４のフィルタ作成プロセスについて説明する。
図９に示すように、繊維性組織の学習領域の複数の超音波ＲＦ信号を、ステップ５２で求めた繊維性組織の基底関数群φfibijおよび脂質性組織の基底関数群φlipijで分解し、基底関数ごとの係数の絶対値を累積する（ステップ５３）。繊維性組織の学習領域の超音波ＲＦ信号を繊維性組織の基底関数群φfibijで分解した各基底関数の係数の絶対値の累積値を図１０（ａ）に示す。また、繊維性組織の学習領域の超音波ＲＦ信号を脂質性組織の基底関数群φlipijで分解した各基底関数の係数の絶対値の累積値を図１０（ｂ）に示す。同様に、脂質性組織の学習領域の超音波ＲＦ信号を、繊維性組織の基底関数群φfibijで分解した各基底関数の係数の絶対値の累積値を図１０（ｃ）に、脂質性組織の基底関数群φlipijで分解した各基底関数の係数の絶対値の累積値を図１０（ｄ）に示す。
【００３１】
続いて、図１０の各係数累積値群を平均値の上下で２値化して２値フィルタを作成する（ステップ５４）。２値化した結果を図１１に示す。ここで、繊維性組織の学習領域−繊維性組織の基底関数群の組み合わせにより求められた２値フィルタをＡ1（ｉ）とする。ｉは、基底関数の番号で、例えば１〜３２の整数である。同様に、繊維性組織の学習領域−脂質性組織の基底関数群の組み合わせによる２値フィルタをＢ1（ｉ）、脂質性組織の学習領域−繊維性組織の基底関数群の組み合わせによる２値フィルタをＡ2（ｉ）、脂質性組織の学習領域−脂質性組織の基底関数群の組み合わせによる２値フィルタをＢ2（ｉ）とする。
【００３２】
＜＜判別プロセス＞＞
図５における、ステップ５５およびステップ５６の判別プロセスについて説明する。
判別したい注目点からの超音波ＲＦ信号（未知信号）を、繊維性組織の基底関数群φfibijおよび脂質性組織の基底関数群φlipijで分解し、基底関数ごとの係数を算出する（ステップ５５）。
続いて、未知信号から判別係数ＦおよびＬを算出して、判別アルゴリズムにより未知信号が繊維性組織の信号か脂質性組織の信号かを判別する（ステップ５６）。未知信号を、繊維性組織の基底関数群φfibijで分解した際の各基底関数の係数をＣ（ｉ）とし、脂質性組織の基底関数群φlipijで分解した際の各基底関数の係数をＤ（ｉ）とすると、判別係数ＦおよびＬは以下の式で求められる。
【数６】
2010075616.tif 000008
【数７】
2010075616.tif 000009
ここで、Ｍは基底関数の数であり、本実施形態においてはＭ＝３２である。
判別係数ＦおよびＬを用いて、図１２に示されたアルゴリズムに基づいて未知信号が繊維性組織の信号か脂質性組織の信号かを判別する。
【００３３】
＜＜実験結果＞＞
本実施形態における判別の分類正解率は以下の通りである。
・繊維性組織からの信号を「繊維性組織」と判定（正解）：  ８９．３６％
・繊維性組織からの信号を「脂質性組織」と判定（不正解）：  ７．８５％
・繊維性組織からの信号を「判別不能」と判定（判別不能）：  ２．７８％
・脂質性組織からの信号を「繊維性組織」と判定（不正解）：１０．５５％
・脂質性組織からの信号を「脂質性組織」と判定（正解）：  ８６．０３％
・脂質性組織からの信号を「判別不能」と判定（判別不能）：  ３．４２％
【００３４】
＜＜変形例＞＞
上述の例では、判別係数ＦおよびＬの算出に、図１１のように２値フィルタＡ1（ｉ）、Ａ2（ｉ）、Ｂ1（ｉ）、Ｂ2（ｉ）を用いたが、図１３のように重み付けされたフィルタを用いても良い。図１３のフィルタを用いた判別の分類正解率は以下の通りである。
・繊維性組織からの信号を「繊維性組織」と判定（正解）：  ９５．５８％
・繊維性組織からの信号を「脂質性組織」と判定（不正解）：  ４．２６％
・繊維性組織からの信号を「判別不能」と判定（判別不能）：  ０．１６％
・脂質性組織からの信号を「繊維性組織」と判定（不正解）：１６．５８％
・脂質性組織からの信号を「脂質性組織」と判定（正解）：  ８３．０５％
・脂質性組織からの信号を「判別不能」と判定（判別不能）：  ０．３７％
【００３５】
以上、本発明の実施形態の一例を説明したが、本発明はこれに限定されるものではなく、特許請求の範囲に記載された技術的思想の範疇において各種の変更が可能であることは言うまでもない。例えば、組織性状ごとの基底関数群を用いた判別アルゴリズムについては、様々なバリエーションが考えられる。また、スパースコーディングによる基底関数や係数累積値の学習はその都度行う必要はなく、事前に学習した組織性状ごとの基底関数や係数累積値などの情報を記憶しておき、それらの情報を用いて未知信号（組織性状を判別したい注目点近傍のＲＦ信号）の組織性状の判別を行っても良い。さらに、学習および判別に用いるデータは超音波ＲＦ信号に限られず、一般的な医療画像などの組織性状の判別に用いても良い。
【図面の簡単な説明】
【００３６】
【図１】血管内超音波エコー法の概略図。
【図２】（ａ）冠動脈の構造、（ｂ）アテロームプラークの説明図。
【図３】（ａ）安定プラーク、（ｂ）不安定プラーク、（ｃ）繊維性被膜の破裂、（ｄ）血栓形成の説明図。
【図４】スパースコーディングにおける基底関数の学習手順のフロー図。
【図５】組織性状の判別の全体フロー図。
【図６】注目点近傍の超音波ＲＦ信号の説明図。
【図７】繊維性組織の基底関数群の学習の説明図。
【図８】脂質性組織の基底関数群の学習の説明図。
【図９】係数累積値群の算出の説明図。
【図１０】係数累積値群の一例のグラフ。
【図１１】２値フィルタの一例のグラフ。
【図１２】判別アルゴリズムのフロー図。
【図１３】重みフィルタの一例のグラフ。
【図１】
2010075616.tif 000010
【図２】
2010075616.tif 000011
【図３】
2010075616.tif 000012
【図４】
2010075616.tif 000013
【図５】
2010075616.tif 000014
【図１０】
2010075616.tif 000015
【図１１】
2010075616.tif 000016
【図１２】
2010075616.tif 000017
【図１３】
2010075616.tif 000018
【図６】
2010075616.tif 000019
【図７】
2010075616.tif 000020
【図８】
2010075616.tif 000021
【図９】
2010075616.tif 000022