(57)【特許請求の範囲】
【請求項１】
病変部の有無が未知である検査画像をシフト不変特徴量でベクトル記述する特徴量算出部と、
前記検査画像の前記シフト不変特徴量を、病変部を含まない複数の正常構造画像のシフト不変特徴量を記述した複数のベクトルから算出される基底ベクトルに係数を乗じた線形和に変換する基底表現部と、
前記検査画像の前記シフト不変特徴量の係数と、前記複数の正常構造画像に含まれる一の正常構造画像から算出されるシフト不変特徴量を前記基底ベクトルに係数を乗じた線形和に変換した際の係数との差が、判定閾値よりも大きい場合に、前記検査画像が病変部を含むと判定する病変判定部と、
前記病変判定部による判定結果を出力する出力部と
を備える診断支援装置。
【請求項２】
前記特徴量算出部は、前記検査画像の着目画素から、複数のシフト不変特徴量を算出し、算出した各シフト不変特徴量をベクトルの要素とする検査画像特徴量ベクトルを算出し、
前記基底表現部は、各々が正常構造の医用画像である複数の正常構造画像の各画素から算出される各シフト不変特徴量をベクトルの要素とする複数の画像特徴量ベクトルの基底ベクトルの線形和で、前記検査画像特徴量ベクトルを記述した際の係数を算出し、算出した前記係数をベクトルの要素とする検査係数ベクトルを算出し、
前記診断支援装置は、さらに、前記画像特徴量ベクトルを前記基底ベクトルの線形和で記述した際の係数をベクトルの要素とする複数の正常係数ベクトルの中から、前記検査係数ベクトルに最も類似する正常係数ベクトルを最近傍ベクトルとして検出する最近傍ベクトル検出部を備え、
前記病変判定部は、前記検査係数ベクトルと前記最近傍ベクトルとの間の距離と前記判定閾値とを比較し、前記距離が前記判定閾値より大きい場合には前記着目画素が病変部の画素であると判定し、前記距離が前記判定閾値より小さい場合には前記着目画素が正常部の画素であると判定する
請求項１記載の診断支援装置。
【請求項３】
前記シフト不変特徴量は、ウェーブレット係数、ＨＬＡＣ特徴量、ＳＩＦＴ特徴量またはＨＯＧ特徴量を含む
請求項２記載の診断支援装置。
【請求項４】
前記医用画像は、放射線画像、超音波画像または病理標本画像を含む
請求項２または３に記載の診断支援装置。
【請求項５】
前記最近傍ベクトル検出部は、前記着目画素から所定距離範囲内に位置する前記複数の正常構造画像の画素に対応する複数の正常係数ベクトルの中から、前記検査係数ベクトルに最も類似する正常係数ベクトルを前記最近傍ベクトルとして検出する
請求項２〜４のいずれか１項に記載の診断支援装置。
【請求項６】
前記判定閾値は、病変部の画素から算出される各シフト不変特徴量をベクトルの要素とする病変画像特徴量ベクトルを、前記基底ベクトルの線形和で記述した際の係数をベクトルの要素とする病変係数ベクトルと、前記複数の正常係数ベクトルのうち前記病変係数ベクトルと最も類似するものとの間の距離の平均値または中央値である
請求項２〜５のいずれか１項に記載の診断支援装置。
【請求項７】
前記病変判定部は、さらに、前記着目画素と、前記最近傍ベクトルの算出の基となった画素との間の距離が検出エラー閾値よりも大きい場合には、検出エラーと判定する
請求項２〜６のいずれか１項に記載の診断支援装置。
【請求項８】
前記特徴量算出部は、前記検査画像の着目画素から、複数のシフト不変特徴量を算出し、算出した各シフト不変特徴量をベクトルの要素とする検査画像特徴量ベクトルを算出し、
前記基底表現部は、各々が正常構造の医用画像である複数の正常構造画像の各画素から算出される各シフト不変特徴量をベクトルの要素とする複数の画像特徴量ベクトルの基底ベクトルの線形和で、前記検査画像特徴量ベクトルを記述した際の係数を算出し、算出した前記係数をベクトルの要素とする検査係数ベクトルを算出し、
前記診断支援装置は、さらに、前記画像特徴量ベクトルを前記基底ベクトルの線形和で記述した際の係数をベクトルの要素とする複数の正常係数ベクトルの中から、前記検査係数ベクトルに最も類似するものから所定個数の正常係数ベクトルを検出する近傍ベクトル検出部を備え、
前記病変判定部は、前記近傍ベクトル検出部が検出した前記所定個数の正常係数ベクトルの各々について、前記検査係数ベクトルと当該正常係数ベクトルとの間の距離と前記判定閾値とを比較し、前記距離が前記判定閾値より大きい場合には前記着目画素が病変部の画素であると判定し、前記距離が前記判定閾値より小さい場合には前記着目画素が正常部の画素であると判定する
請求項１記載の診断支援装置。
【請求項９】
病変部の有無が未知である検査画像をシフト不変特徴量でベクトル記述し、
前記検査画像の前記シフト不変特徴量を、病変部を含まない複数の正常構造画像のシフト不変特徴量を記述した複数のベクトルから算出される基底ベクトルに係数を乗じた線形和に変換し、
前記検査画像の前記シフト不変特徴量の係数と、前記複数の正常構造画像に含まれる一の正常構造画像から算出されるシフト不変特徴量を前記基底ベクトルに係数を乗じた線形和に変換した際の係数との差が、判定閾値よりも大きい場合に、前記検査画像が病変部を含むと判定し、
前記病変部の判定結果を出力する
診断支援方法。
【請求項１０】
請求項９に記載の診断支援方法をコンピュータに実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、医師による画像診断を支援する診断支援装置および診断支援方法に関する。
【背景技術】
【０００２】
従来、医師による画像診断を支援するために、正常構造の医用画像である正常構造画像と、検査対象の医用画像である検査画像とを比較して、病変部を判定する装置が提案されている（例えば、特許文献１〜３参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開平７−３７０７４号公報
【特許文献２】特開２００２−３２７３５号公報
【特許文献３】特開２００４−４１６９４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、従来の装置では、正常構造画像と検査画像との位置あわせを行わなければならず、正確に位置あわせを行うことができなければ、病変部の判定精度が低いという課題がある。【０００５】
本発明は従来の課題を解決するためになされたものであり、病変部の判定精度が高い診断支援装置を提供する。
【課題を解決するための手段】
【０００６】
本発明の一態様に係る診断支援装置は、病変部の有無が未知である検査画像をシフト不変特徴量でベクトル記述する特徴量算出部と、前記検査画像の前記シフト不変特徴量を、病変部を含まない複数の正常構造画像のシフト不変特徴量を記述した複数のベクトルから算出される基底ベクトルに係数を乗じた線形和に変換する基底表現部と、前記検査画像の前記シフト不変特徴量の係数と、前記複数の正常構造画像に含まれる一の正常構造画像から算出されるシフト不変特徴量を前記基底ベクトルに係数を乗じた線形和に変換した際の係数との差が、判定閾値よりも大きい場合に、前記検査画像が病変部を含むと判定する病変判定部と、前記病変判定部による判定結果を出力する出力部とを備える。【０００７】
この構成によると、検査画像と正常構造画像との比較を、シフト不変特徴量を基底表現した際の係数間で行っている。このため、正常構造画像と検査画像の位置合わせを行うことなく、病変の有無を判定できる。位置合わせ工程が不要となることにより、ランドマークの設定方法の違いや設定位置のばらつきによる病変部の判定精度の低下が無くなり、病変部の判定精度が高い診断支援装置を提供することができる。【０００８】
なお、これらの包括的または具体的な態様は、システム、方法、集積回路、コンピュータプログラムまたはコンピュータ読み取り可能なＣＤ−ＲＯＭなどの記録媒体で実現されてもよく、システム、方法、集積回路、コンピュータプログラムおよび記録媒体の任意な組み合わせで実現されてもよい。
【発明の効果】
【０００９】
本発明によると、病変部の判定精度が高い診断支援装置を提供することができる。
【図面の簡単な説明】
【００１０】
【図１】図１は、特許文献１におけるランドマーク設定を説明する図である。
【図２】図２は、本発明の実施の形態１に係る診断支援装置の機能的な構成を示すブロック図である。
【図３】図３は、ウェーブレット変換によるシフト不変特徴量の算出例を示す図である。
【図４】図４は、ハール型のマザーウェーブレットを用いたウェーブレット係数の算出例を示す図である。
【図５】図５は、正常構造の基底ベクトルと平均ベクトルを算出する方法を説明する図である。
【図６】図６は、病変有無の判定閾値の決定方法の一例を示す図である。
【図７】図７は、正常係数ベクトルと病変係数ベクトルとの間の距離のヒストグラムを示す図である。
【図８】図８は、診断支援装置の実行する処理のフローチャートである。
【図９】図９は、実施の形態２に係る診断支援装置の機能的な構成を示すブロック図である。
【図１０】図１０は、類似症例検索サーバーによる類似症例データの検索方法の一例を説明するための図である。
【図１１】図１１は、表示部による判定結果および類似症例データの表示例を示す図である。
【図１２】図１２は、実施の形態３に係る診断支援装置の機能的な構成を示すブロック図である。
【図１３】図１３は、実施の形態４に係る診断支援装置の機能的な構成を示すブロック図である。
【図１４】図１４は、表示部による判定結果の表示例を示す図である。
【図１５】図１５は、診断支援装置を実現するコンピュータシステムのハードウェア構成を示すブロック図である。
【図１６】図１６は、本発明に係る診断支援装置の必須の構成要素を示す図である。
【発明を実施するための形態】
【００１１】
（本発明の基礎となった知見）
本発明者は、「背景技術」の欄において記載した、従来の装置に関し、以下の問題が生じることを見出した。
【００１２】
医療業務は診断業務と治療業務の２つに大きく分かれる。
【００１３】
診断業務は病名と病状の把握が目的であり、病名および病状に応じて治療方針が決定される。診断業務は問診または触診から始まり、病気の疑いがあれば画像診断でより詳細な検査が行われる。画像診断では放射線または超音波などを用いた非破壊検査により体内の状態を診る。また内視鏡などによって臓器表面が観察される。病名および病状を確定するには、病変部から検体を摘出して細胞レベルで状態観察を行う病理診断が行われる。摘出した検体は顕微鏡で観察できる厚みにスライスされ、標本が作製される。作製された標本は顕微鏡越しにデジタルカメラまたはスキャナで撮影され、デジタル画像として保存、参照される。
【００１４】
以上のような医用画像がデジタル化されたことで、医用画像と、コンピュータによるデータ処理との整合性が高まり、医師や技師の診断業務をＩＴシステムで支援できる機会が増えてきた。ＣＡＤ（Ｃｏｍｐｕｔｅｒ  Ａｉｄｅｄ  Ｄｅｔｅｃｔｉｏｎ）はその一例であり、病変部の検出にコンピュータを活用する方法である。
【００１５】
診断医は通常、病変部を含まない正常構造のみから成る医用画像を記憶している。検査対象の医用画像である検査画像が提示されると、診断医は正常構造の医用画像である正常構造画像を思い浮かべて検査画像との比較を行う。診断医は、正常構造画像との違いを発見すると、違う箇所が病変部であると判断する。コンピュータ処理において、２つのデータの差分を計算することは基本的な機能であり、正常構造画像と検査画像との比較から病変部を検出する画像診断は、コンピュータ処理が得意とするところである。
【００１６】
しかし、差分計算を行うには正常構造画像と検査画像との位置合わせを行っておく必要がある。正常構造画像を作るには、通常、過去の医用画像が利用される。検査画像を診断医が確認した結果、病変部が見つからなければ、この画像は正常構造画像となる。患者が過去に撮影した検査画像に病変部が存在しなければ、この検査画像を正常構造画像として利用することができる。しかし、撮影条件の違いまたは体型の変化などの様々な要因から、同一患者であっても正常構造画像と検査画像との位置が食い違うことが通常である。また、検査画像を初めて撮影する場合には比較対象とする画像が存在せず、正常構造画像を得ることはできない。この場合には、他の患者の正常構造画像を利用することになるが、他の患者との体型の違いから正常構造画像と検査画像との位置合わせが必要になる。
【００１７】
位置合わせは、具体的には、回転、平行移動、拡大・縮小などの幾何変換により実現される。例えば、正常構造画像と検査画像とにそれぞれ対応点を複数設定し、各画像の対応点が一致するように、一方の画像にアフィン変換と拡大・縮小処理とを施す。また、たとえば特許文献１および２が開示するように、局所領域ごとのテンプレートマッチングによる非線形変換方法、または大局的位置合わせと局所的位置合わせを組み合わせた変換方法により位置合わせが行われる。
【００１８】
ところで、正常構造画像は過去に撮影され、病変がないことを確認した画像から作成される。その理由は、第一に、初めて撮影された患者の検査画像には比較対象となる医用画像が存在しないためであり、上述したとおりである。第二に、医学的知識は過去症例の積み上げで構築される傾向にあり、病変を含まない正常構造に関しても、過去症例から作成した方が、医学的利用価値が高くなる可能性が高いからである。医学的知識は、日進月歩で向上しており、過去症例の解釈が改善されることがしばしば起こる。ＩＴシステムに登録した医学的知識は、常に更新が必要であり、正常構造画像も例外ではない。
【００１９】
そこで、複数の患者の正常構造画像を集めて、これらを網羅的に表現できるような汎用性の高い正常構造画像を作ることが望ましい。その具体的実現方法として、たとえば特許文献３が開示するように、正常構造画像を平均形状と固有形状の線形和で表現する。すなわち、正常構造を表わす形状ベクトルｘは以下の式１で表現される。
【００２０】
【数１】
0005607839.tif 000002
【００２１】
ここで、ｘａｖｅは平均形状ベクトル、Ｐｓは形状の固有ベクトル、ｂｓは形状係数のセットを表わす。
【００２２】
式１を計算するには、平均形状ベクトルｘａｖｅと形状の固有ベクトルＰｓが必要であり、画像情報をベクトル化するために、図１に示すようなランドマークＭを画像上に設定する。図中の複数の黒点がランドマークＭである。ランドマークＭのｘｙ座標がベクトルの要素となり、画像情報がベクトル化される。検査画像Ｐ１、Ｐ２、Ｐ３に示すように、複数の正常構造画像で個別にランドマークを設定し、形状ベクトルを定義し、これらから平均ベクトルと固有ベクトルを算出する。なお、検査画像についても同様に式１で表すことができる。
【００２３】
以上説明したベクトルを用いて、正常構造画像と検査画像の位置合わせを行い、正常構造画像と検査画像の差分から病変部を検出して、画像診断の支援を行う。
【００２４】
しかしながら、従来の技術では、ランドマークの設定作業が煩雑であり、診断業務の効率低下につながる。また、ランドマークの設定方法を変えれば、ランドマークのｘｙ座標が変わりベクトル要素が変わるため、結果的に平均形状ベクトルおよび固有形状ベクトルがランドマークの設定方法の変更前と異なってくる。特許文献３は、ランドマークの設定方法について開示していないため、同じ手法を利用したとしても、ランドマークの設定方法の違いや、設定位置のばらつきで、様々な正常構造画像（正常構造を表す形状ベクトル）が生成されてしまう。医学的知見は過去症例からの蓄積で構築されており、ある一つの症例がランドマークの設定方法に応じて多様に定義されるのは再利用性の観点で課題がある。このように一つの症例が多様に定義されると、病変部の判定精度が低下するという課題がある。つまり、定義によって同じ検査画像が病変画像と判断されたり正常画像と判断されたりしてしまう。
【００２５】
このような問題を解決するために、本発明の一態様に係る診断支援装置は、病変部の有無が未知である検査画像をシフト不変特徴量でベクトル記述する特徴量算出部と、前記検査画像の前記シフト不変特徴量を、病変部を含まない複数の正常構造画像のシフト不変特徴量を記述した複数のベクトルから算出される基底ベクトルに係数を乗じた線形和に変換する基底表現部と、前記検査画像の前記シフト不変特徴量の係数と、前記複数の正常構造画像に含まれる一の正常構造画像から算出されるシフト不変特徴量を前記基底ベクトルに係数を乗じた線形和に変換した際の係数との差が、判定閾値よりも大きい場合に、前記検査画像が病変部を含むと判定する病変判定部と、前記病変判定部による判定結果を出力する出力部とを備える。
【００２６】
この構成によると、検査画像と正常構造画像との比較を、シフト不変特徴量を基底表現した際の係数間で行っている。このため、正常構造画像と検査画像の位置合わせを行うことなく、病変の有無を判定できる。位置合わせ工程が不要となることにより、ランドマークの設定方法の違いや設定位置のばらつきによる病変部の判定精度の低下が無くなり、病変部の判定精度が高い診断支援装置を提供することができる。
【００２７】
例えば、前記特徴量算出部は、前記検査画像の着目画素から、複数のシフト不変特徴量を算出し、算出した各シフト不変特徴量をベクトルの要素とする検査画像特徴量ベクトルを算出し、前記基底表現部は、各々が正常構造の医用画像である複数の正常構造画像の各画素から算出される各シフト不変特徴量をベクトルの要素とする複数の画像特徴量ベクトルの基底ベクトルの線形和で、前記検査画像特徴量ベクトルを記述した際の係数を算出し、算出した前記係数をベクトルの要素とする検査係数ベクトルを算出し、前記診断支援装置は、さらに、前記画像特徴量ベクトルを前記基底ベクトルの線形和で記述した際の係数をベクトルの要素とする複数の正常係数ベクトルの中から、前記検査係数ベクトルに最も類似する正常係数ベクトルを最近傍ベクトルとして検出する最近傍ベクトル検出部を備え、前記病変判定部は、前記検査係数ベクトルと前記最近傍ベクトルとの間の距離と前記判定閾値とを比較し、前記距離が前記判定閾値より大きい場合には前記着目画素が病変部の画素であると判定し、前記距離が前記判定閾値より小さい場合には前記着目画素が正常部の画素であると判定してもよい。
【００２８】
この構成によると、シフト不変特徴量を用い、検査画像と正常構造画像との比較を基底係数ベクトル間で行っている。このため、正常構造画像と検査画像の位置合わせを行うことなく、病変の有無を判定できる。位置合わせ工程が不要となることにより、ランドマークの設定方法の違いや設定位置のばらつきによる病変部の判定精度の低下が無くなり、病変部の判定精度が高い診断支援装置を提供することができる。
【００２９】
また、前記シフト不変特徴量は、ウェーブレット係数、ＨＬＡＣ（Ｈｉｇｈｅｒ  ｏｒｄｅｒ  Ｌｏｃａｌ  ＡｕｔｏＣｏｒｒｅｌａｔｉｏｎ）特徴量、ＳＩＦＴ（Ｓｃａｌｅ−Ｉｎｖａｒｉａｎｔ  Ｆｅａｔｕｒｅ  Ｔｒａｎｓｆｏｒｍ）特徴量またはＨＯＧ（Ｈｉｓｔｏｇｒａｍｓ  ｏｆ  Ｏｒｉｅｎｔｅｄ  Ｇｒａｄｉｅｎｔｓ）特徴量を含んでいても良い。
【００３０】
この構成によると、画素位置に依存せず画像情報をベクトル化でき、ランドマークなどを用いる位置合わせ工程が不要になる。
【００３１】
また、前記医用画像は、放射線画像、超音波画像または病理標本画像を含んでいても良い。
【００３２】
この構成によると、病名や病状の診断に必要なあらゆる医用画像が利用でき、診断支援の効率と質を向上できる。
【００３３】
また、前記最近傍ベクトル検出部は、前記着目画素から所定距離範囲内に位置する前記複数の正常構造画像の画素に対応する複数の正常係数ベクトルの中から、前記検査係数ベクトルに最も類似する正常係数ベクトルを前記最近傍ベクトルとして検出しても良い。
【００３４】
この構成によると、正常係数ベクトルの参照対象を制限することができるため、最近傍ベクトルの検出時間の短縮につながる。
【００３５】
また、前記判定閾値は、病変部の画素から算出される各シフト不変特徴量をベクトルの要素とする病変画像特徴量ベクトルを、前記基底ベクトルの線形和で記述した際の係数をベクトルの要素とする病変係数ベクトルと、前記複数の正常係数ベクトルのうち前記病変係数ベクトルと最も類似するものとの間の距離の平均値または中央値であっても良い。
【００３６】
この構成によると、過去の診断結果に裏付けされた医学的精度の高い判定ができる。
【００３７】
また、前記病変判定部は、さらに、前記着目画素と、前記最近傍ベクトルの算出の基となった画素との間の距離が検出エラー閾値よりも大きい場合には、検出エラーと判定しても良い。
【００３８】
上記距離が大きすぎる場合には、着目画素が含まれる臓器とは異なる臓器の画素から最近傍ベクトルが検出された可能性が高いため、正確な病変有無の判定を行うことができない。このため、検出エラー閾値を設けることにより、病変有無の判定精度を高めることができる。
【００３９】
また、前記特徴量算出部は、前記検査画像の着目画素から、複数のシフト不変特徴量を算出し、算出した各シフト不変特徴量をベクトルの要素とする検査画像特徴量ベクトルを算出し、前記基底表現部は、各々が正常構造の医用画像である複数の正常構造画像の各画素から算出される各シフト不変特徴量をベクトルの要素とする複数の画像特徴量ベクトルの基底ベクトルの線形和で、前記検査画像特徴量ベクトルを記述した際の係数を算出し、算出した前記係数をベクトルの要素とする検査係数ベクトルを算出し、前記診断支援装置は、さらに、前記画像特徴量ベクトルを前記基底ベクトルの線形和で記述した際の係数をベクトルの要素とする複数の正常係数ベクトルの中から、前記検査係数ベクトルに最も類似するものから所定個数の正常係数ベクトルを検出する近傍ベクトル検出部を備え、前記病変判定部は、前記近傍ベクトル検出部が検出した前記所定個数の正常係数ベクトルの各々について、前記検査係数ベクトルと当該正常係数ベクトルとの間の距離と前記判定閾値とを比較し、前記距離が前記判定閾値より大きい場合には前記着目画素が病変部の画素であると判定し、前記距離が前記判定閾値より小さい場合には前記着目画素が正常部の画素であると判定してもよい。
【００４０】
この構成によると、検査画像に最も類似する正常構造画像から得られる判定結果のみならず、２番目以降に類似する正常構造画像から得られる判定結果も出力することができる。
【００４１】
なお、これらの包括的または具体的な態様は、システム、方法、集積回路、コンピュータプログラムまたはコンピュータ読み取り可能なＣＤ−ＲＯＭなどの記録媒体で実現されてもよく、システム、方法、集積回路、コンピュータプログラムまたは記録媒体の任意な組み合わせで実現されてもよい。
【００４２】
以下、実施の形態について、図面を参照しながら具体的に説明する。
【００４３】
なお、以下で説明する実施の形態は、いずれも包括的または具体的な例を示すものである。以下の実施の形態で示される数値、形状、材料、構成要素、構成要素の配置位置及び接続形態、ステップ、ステップの順序などは、一例であり、本発明を限定する主旨ではない。また、以下の実施の形態における構成要素のうち、最上位概念を示す独立請求項に記載されていない構成要素については、任意の構成要素として説明される。
【００４４】
（実施の形態１）
本実施の形態では、対象の位置に依存しないシフト不変特徴量をベクトルの要素として、画像情報をベクトル表現することで、位置合わせのための対応点設定が不要である診断支援装置について説明する。
【００４５】
図２は、本発明の実施の形態１に係る診断支援装置の機能的な構成を示すブロック図である。
【００４６】
診断支援装置１００は、入力部１０１、特徴量算出部１０２、基底表現部１０３、最近傍ベクトル検出部１０４、正常係数ベクトル格納部１０５、病変判定部１０６、出力部１１１、表示部１０７を有し、正常構造画像と検査画像の異なる箇所である病変部を検出する。
【００４７】
入力部１０１は、診断に必要となる医用画像である検査画像を受け付ける。本実施の形態では医用画像のタイプに制限を設けず、任意の医用画像を対象とする。たとえば、ＣＴ（Ｃｏｍｐｕｔｅｒ  Ｔｏｍｏｇｒａｐｈｙ）画像、ＭＲ（Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ）画像、ＰＥＴ（Ｐｏｓｉｔｒｏｎ  Ｅｍｉｓｓｉｏｎ  Ｔｏｍｏｇｒａｐｈｙ）画像、超音波画像、病理画像、放射線画像の一種である胸部単純Ｘ線画像などを対象とする。
【００４８】
特徴量算出部１０２は、検査画像の着目画素から、複数のシフト不変特徴量を算出し、算出した各シフト不変特徴量をベクトルの要素とする検査画像特徴量ベクトルを算出する。つまり、特徴量算出部１０２は、シフト不変特徴量を使って検査画像の画像情報をベクトル化し、検査画像特徴量ベクトルｆｐを出力する。特徴量算出部１０２は、検査画像の画素ごとに検査画像特徴量ベクトルｆｐを出力する。
【００４９】
ここで、画素値群からシフト不変特徴量を算出する方法を、ウェーブレット変換を例に取って説明する。図３は、ウェーブレット変換によるシフト不変特徴量の算出例を示す図である。特徴量算出部１０２は、ウェーブレット変換により、検査画像の、スケールｔの多重解像度表現を行う。スケール１では、隣接画素との輝度差が算出され、スケール２へ移行する際に複数画素での平滑化が施される。スケール２でも、隣接画素との輝度差が算出されるが、スケール２の各画素は、スケール１の複数画素が平滑化されたものであり、スケール２の方が低い周波数成分となる。したがって、スケール１からスケールｔ（ｔは２以上の整数）まで計算を進めることで、高周波成分から徐々に低周波数成分へ移行しながら、ウェーブレット係数Ｖ、Ｈ、Ｄが各スケールで算出される。特徴量算出部１０２は、画素ごとに、各ステージで算出されたウェーブレット係数Ｖ、Ｈ、Ｄと、スケールｔの画像から算出される輝度平均値Ｌと、からなる空間周波数ベクトルＦを生成する。つまり、空間周波数ベクトルＦの次元数は（３ｔ＋１）次元になる。
【００５０】
ハール型のマザーウェーブレットを用いた場合、図４（ａ）に示すように、Ｖは処理対象の画素である着目画素３０の右隣接画素３１との輝度差分値、Ｈは着目画素３０の下隣接画素３２との輝度差分値、Ｄは右斜下隣接画素３３との輝度差分値、Ｌは着目画素３０、右隣接画素３１、下隣接画素３２、右斜下隣接画素３３からなる４つの画素の輝度平均値である。図４（ａ）はスケール１に相当し、図４（ｂ）はスケール２に相当する。スケール２の検査画像は、スケール１の検査画像の４つの画素の輝度平均値を１画素とする画像である。つまり、スケール２では、スケール１における４画素の輝度平均値である出力Ｌが輝度差分値の計算対象であるブロックの輝度値になる。スケール２における出力Ｖはブロック３４と右隣接ブロック３５との輝度差分値である。スケール２における出力Ｈはブロック３４と下隣接ブロック３６との輝度差分値である。スケール２における出力Ｄはブロック３４と右斜下隣接ブロック３７との輝度差分値である。また、スケール２における出力Ｌはブロック３４から右斜下隣接ブロック３７の４つのブロックの輝度平均値である。
【００５１】
以上より、ウェーブレット変換を用いた場合、検査画像特徴量ベクトルｆｐは図３の空間周波数ベクトルＦとして算出される。
【００５２】
なお、本実施の形態ではシフト不変特徴量として、ウェーブレット係数を用いるが、シフト不変特徴量はこれに限定されるものではなく、任意のシフト不変特徴量を用いることができる。たとえば、シフト不変特徴量として、ＳＩＦＴ特徴量、ＨＬＡＣ特徴量、またはＨＯＧ特徴量などを用いることができる。
【００５３】
基底表現部１０３は、検査画像特徴量ベクトルｆｐを基底ベクトル表現に変換し、検査係数ベクトルαｐを出力する。つまり、基底表現部１０３は、複数の正常構造画像の各画素から算出される各シフト不変特徴量をベクトルの要素とする複数の画像特徴量ベクトルの基底ベクトルの線形和で、検査画像特徴量ベクトルｆｐを記述した際の係数を算出し、算出した係数をベクトルの要素とする検査係数ベクトルαｐを算出する。基底表現部１０３は、検査画像特徴量ベクトルｆｐを、以下の式２によって、検査係数ベクトルαｐへ変換する。
【００５４】
【数２】
0005607839.tif 000003
【００５５】
ここで、行列Ｂは後述する正常構造基底ベクトル行列、ベクトルｇは後述する正常構造平均ベクトルを表わす。式２は式１を形状係数のセットｂｓについて解いたものである。式１と式２の対応関係は以下のとおりである。
【００５６】
形状ベクトルｘ  ⇔  検査画像特徴量ベクトルｆｐ
平均形状ベクトルｘａｖｅ  ⇔  正常構造平均ベクトルｇ
形状の固有ベクトルＰｓ  ⇔  正常構造基底ベクトル行列Ｂ
形状係数のセットｂｓ  ⇔  検査係数ベクトルαｐ
【００５７】
正常構造基底ベクトル行列Ｂと正常構造平均ベクトルｇは、図５に示すように、多数の正常構造画像から得られた画像特徴量ベクトルをもとに算出される。例えば、正常構造画像の幅がＷ、高さがＨとすると、１枚の正常構造画像から（Ｗ×Ｈ）本の画像特徴量ベクトルが算出される。正常構造画像の枚数をＱ枚とすると、Ｑ枚の正常構造画像から（Ｗ×Ｈ×Ｑ）本の画像特徴量ベクトルが得られる。画像特徴量ベクトルの次元数はｎ次元とする。
【００５８】
正常構造平均ベクトルｇは、画像特徴量ベクトルの要素ごとに平均値を算出することで求められる。
【００５９】
正常構造基底ベクトル行列Ｂは、主成分分析によって以下の式３の連立方程式の解である固有ベクトルｂ１、ｂ２、・・・、ｂｎとして算出される。
【００６０】
【数３】
0005607839.tif 000004
【００６１】
ここで、行列Ｓは分散共分散行列であり、以下の式４で与えられる。
【００６２】
【数４】
0005607839.tif 000005
【００６３】
ｓｉは画像特徴量ベクトルのｉ次元要素の分散である。上述したように、画像特徴量ベクトルは（Ｗ×Ｈ×Ｑ）本だけ得られ、よって、画像特徴量ベクトルのｉ次元要素は（Ｗ×Ｈ×Ｑ）個だけ存在する。ｓｉはその分散である。
【００６４】
また、固有値λは以下の式５で与えられる。
【００６５】
【数５】
0005607839.tif 000006
【００６６】
固有値λはｎ個求められ、大きい値から順番にλ１、λ２、…、λｎとする。
【００６７】
以上より、基底表現部１０３は、検査画像の画像特徴量ベクトル（検査画像特徴量ベクトルｆｐ）を基底表現に変換し、検査係数ベクトルαｐを生成する。これにより、検査画像と正常構造画像との比較が可能になる。基底表現部１０３の出力である検査係数ベクトルαｐは、最近傍ベクトル検出部１０４と病変判定部１０６へ入力され、正常構造画像との比較に用いられる。
【００６８】
なお、式２の計算に必要となる正常構造基底ベクトル行列Ｂと正常構造平均ベクトルｇは、予め計算しておき、基底表現部１０３に格納しておく。異常所見が認められない画像が正常構造画像となるが、画像の一部に病変部があっても画像特徴量を算出する領域が病変部を含まなければ、正常構造基底ベクトル行列Ｂや正常構造平均ベクトルｇの算出のサンプルになり得る。図４の場合、ブロック３４から右斜下隣接ブロック３７が画像特徴量の算出に関与する。この４つのブロック内に病変部が含まれなければ、正常構造基底ベクトル行列Ｂおよび正常構造平均ベクトルｇの算出のサンプルになる。そこで、病変部を線で囲むことで、正常構造と病変部は分離できる。また、病変部を囲んだ線をまたがないように画像特徴量を求める領域を制御すれば、正常構造基底ベクトル行列Ｂおよび正常構造平均ベクトルｇを正しく算出できる。病変部を囲む線の入力方法は限定されるものではなく、たとえば医師が診断の際に、マウス等で線を引くようにしても大きな負担にならない。あるいは、所見に部位を表す医学用語が記載されている場合がある。医学用語と医用画像の位置関係を示すデータをあらかじめ用意しておけば、人手を介さずに病変部の位置を特定することができ、正常構造と病変部を分離することができる。
【００６９】
さらに、「病変部である」と診断された過去症例を参照し、病変領域の画像特徴量ベクトルをあらかじめ計算する。前記病変領域の画像特徴量ベクトルと、正常構造を取得したい対象画像の画像特徴量ベクトルとの差が小さい場合は、病変部である可能性が高い。このため、２つのベクトル間の距離が所定の閾値よりも小さい場合は、対象画像を正常構造画像に加えない。
【００７０】
最近傍ベクトル検出部１０４は、正常係数ベクトル格納部１０５に格納された正常係数ベクトルαの中から、検査係数ベクトルαｐに最も類似するものを検出する。正常係数ベクトルαは、正常構造画像の画像特徴量ベクトルを、当該画像特徴量ベクトルの基底ベクトルの線形和で記述した際の係数をベクトルの要素とする。
【００７１】
正常係数ベクトルαは、以下の式６で算出する。
【００７２】
【数６】
0005607839.tif 000007
【００７３】
ここで、ベクトルｆは正常画像特徴量ベクトルである。正常画像特徴量ベクトルｆは、検査画像の代わりに正常構造画像を用いることにより、検査画像特徴量ベクトルｆｐと同様に計算される。式６の構成は式２と同一であり、式２における検査画像特徴量ベクトルｆｐの代わりに正常画像特徴量ベクトルｆを与えると、正常係数ベクトルαが算出される。
【００７４】
正常係数ベクトル格納部１０５は、式６の計算結果である正常係数ベクトルαを正常構造画像の枚数分だけ格納している。
【００７５】
最近傍ベクトル検出部１０４は、検査係数ベクトルαｐと各正常係数ベクトルαとの間の距離を計算し、最も距離が短い最近傍の正常係数ベクトルαを検出する。最近傍ベクトル検出部１０４は、最近傍の正常係数ベクトルαが検出されたら、これをベクトル間の距離とともに病変判定部１０６へ送る。
【００７６】
病変判定部１０６は、検査画像の画素ごとに、検査係数ベクトルαｐと正常係数ベクトルαとの間の距離に基づいて、病変部の有り、無しを判定する。つまり、病変判定部１０６は、検査係数ベクトルαｐと正常係数ベクトルαとの間の距離と判定閾値とを比較する。病変判定部１０６は、算出した距離が判定閾値より大きい場合には、検査係数ベクトルαｐが算出された着目画素が病変部の画素であると判定し、距離が判定閾値より小さい場合には着目画素が正常部の画素であると判定する。なお、病変有無の判定閾値は過去症例から算出され、第１閾値データベース１０９に予め格納されている。
【００７７】
図６は、病変有無の判定閾値の決定方法の一例を示す図である。判定閾値は、特徴量算出部１０２と、基底表現部１０３と、最近傍ベクトル検出部１０４と、正常係数ベクトル格納部１０５と、ベクトル距離算出部５０１と、検出エラー閾値算出部５０２とを備える判定閾値決定装置が決定する。
【００７８】
まず、病変部の存在が確認された病変画像が過去症例から取得され、病名ごとに分類される。次に、特徴量算出部１０２が、たとえば、病名Ｄに分類された病変画像Ｉｄの画像座標ごと（画素ごと）に複数の画像特徴量を算出し、算出した各画像特徴量をベクトルの要素とする病変画像特徴量ベクトルｆｄ（ｆｄ，１〜ｆｄ，ｎａ）を生成する。病変画像特徴量ベクトルｆｄは画素ごとに生成される。基底表現部１０３は、病変画像特徴量ベクトルｆｄを、式２のベクトルｆｐに代入し、ベクトルαｐに変換する。このベクトルαｐを病変係数ベクトルαｄとする。次に、最近傍ベクトル検出部１０４は、正常係数ベクトル格納部１０５に格納されている正常係数ベクトルαの中から、病変係数ベクトルαｄに最も類似する正常係数ベクトルαを検出する。例えば、病変係数ベクトルαｄとの間の距離が最も小さい正常係数ベクトルαが検出される。ベクトル距離算出部５０１は、検出された最近傍の正常係数ベクトルαと病変係数ベクトルαｄとの間の距離を算出し、算出した距離を病変有無の判定閾値とする。ただし、同一病名に分類される病変画像は多数存在するため、各病変部の各画素から算出される判定閾値の、平均値または中央値などを判定閾値とするのが妥当である。
【００７９】
診断支援装置１００の病変判定部１０６は、検出エラーを防ぐために、さらに、検査係数ベクトルαｐおよび最近傍として検出された正常係数ベクトルαの各々について、算出の基となった画素の位置を確認する。図６に示すように、判定閾値決定装置の最近傍ベクトル検出部１０４は、最近傍の正常係数ベクトルαが検出された際、その画像座標を最近傍ベクトル画像座標として出力し、第１閾値データベース１０９に格納する。画像座標は、ＣＴ画像であれば、スライス位置を表すｚ座標とスライス面内の位置を表すｘ、ｙ座標から成る。胸部単純Ｘ線画像はスライスを行わないため、たとえばｚ座標を１に固定すればいい。
【００８０】
本実施の形態では、シフト不変特徴量の一致をもって間接的に位置合わせを行っている。このため、複数の正常構造画像または検査画像を直接変形させて対応点が同じ画像座標となるような変形処理は行っていない。したがって、図１に示すように、正常構造画像や検査画像の対応点は異なる画像座標を有する場合がほとんどである。しかし、その画像座標の違いはある程度の大きさに収まるはずであり、大きくかけ離れた位置に対応点が存在することはない。たとえば、胸部単純Ｘ線画像の場合、撮影方法には一定のガイドラインがあり、その範囲の中で体の向きが少し異なる、という状況で撮影される。ＣＴおよびＭＲＩ（Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  Ｉｍａｇｉｎｇ）では患者がベッドに横たわり、患者と撮影装置の位置関係が固定化される。このため、ＣＴ画像およびＭＲ画像は、胸部単純Ｘ線画像よりも撮影条件による画像座標の違いは小さくなる。一方、患者の体型は様々であるため、異なる患者間で正常構造画像と検査画像を比べた場合、対応点の画像座標は一致しない。しかし、肺と心臓の位置関係、または肺と肩甲骨の位置関係など、臓器および骨などの相対的位置関係は患者間で同一であり、対応点の画像座標が著しく異なることはない。
【００８１】
以上の考え方から、病変判定部１０６は、検査係数ベクトルαｐの算出の基となった画素と最近傍として検出された正常係数ベクトルαの算出の基となった画素との間の距離が、あらかじめ設定した検出エラー閾値よりも大きい場合には、検出エラーと判定する。病変判定部１０６は、検出エラーと判定した場合に検出エラー通知信号１０８を最近傍ベクトル検出部１０４へ出力することにより、最近傍ベクトル検出部１０４に最近傍ベクトルの再検出を指示する。検出エラー通知信号１０８を受信した最近傍ベクトル検出部１０４は、他の最近傍ベクトルを検出する。
【００８２】
検出エラーの基準になる検出エラー閾値は、第２閾値データベース１１０に格納されている。検出エラー閾値の決め方は限定されるものではないが、エラー閾値は、たとえば病名ごとに設定される。図６に示すように、過去症例から病変部の画像座標は入手できる。このため、検出エラー閾値算出部５０２は、病名Ｄの全データの画像座標の分布のうち、最も離れている２つの画像座標間の距離を算出し、これを検出エラー閾値とする。あるいは、検出エラー閾値算出部５０２は、画像座標の分布の重心から最も離れている画像座標との間の距離を検出エラー閾値とする。また、検出エラー閾値算出部５０２は、患者の各臓器の位置を大まかに画像座標（ｘ，ｙ，ｚ）に従い分類し、ある臓器から別の臓器までの距離を求めて、これをエラー検出閾値としてもいい。
【００８３】
加えて、最近傍ベクトル検出部１０４は、最近傍ベクトルを検出する際に、処理対象の画素（対象画素）から所定距離範囲内に位置する画素に対応する複数の正常係数ベクトルαの中から、最近傍ベクトルを検出するようにしても良い。これにより、最近傍ベクトルの検出時間を短くすることができる。それと同時に検索エラーの発生を抑えることもできる。すなわち、最近傍ベクトルの参照範囲を限定することで、正常構造画像と検査画像の対応点の画像座標が臓器の大きさを超えるほど異ならないようにすることができる。このため、最近傍ベクトル検出部１０４は、検査係数ベクトルαｐの画像座標を中心に所定距離範囲内の画素の正常係数ベクトルαを参照すれば十分である。参照範囲の設定方法は限定されるものではないが、たとえば、第２閾値データベース１１０に格納されている検出エラー閾値の最大値を上記所定距離として参照範囲を設定しても良い。第２閾値データベース１１０には病名ごとに検出エラー閾値が格納されており、それらは病名ごとに大きさが異なる。しかし、検出エラー閾値の最大値を上記所定距離とすることであらゆる病名に対応することができる。このため、参照範囲の限定を行なっても病変部の見落としは発生しない。なお、臨床情報または既往症など、画像情報以外の患者情報から病変部検出に関する条件を絞り込むことができれば、病変部の検出時間をできるだけ短くするために、正常係数ベクトルαの参照範囲をさらに制限できる。
【００８４】
図６では、病変画像を病名ごとに分類したが、病変画像に付随する所見などを利用して病状ごとに病変画像を分類することができれば、病状判定用の判定閾値を算出できる。たとえば、腫瘍の病変画像を良性腫瘍の病変画像と悪性腫瘍の病変画像に二分し、それぞれの病変画像について判定閾値を求める。この２つの判定閾値を用いて病変判定部１０６が判定を行うことにより、腫瘍が良性であるか、悪性であるかの病状判定ができる。
【００８５】
なお、これまでの説明では、病変判定部１０６は、病変部が有るか、無いかを二択で判定したが、病変存在の可能性を示す値を算出することも可能である。図７は、病名Ｄに分類された正常係数ベクトルαと病変係数ベクトルαｄとの間の距離のヒストグラムを表す。距離の平均値および中央値ともに１．０となり、病変部の有無を二択で判定する場合は、判定閾値が１．０となる。頻度の下に括弧書きで記された百分率は累積頻度の全体に対する割合である。例えば、病変判定部１０６は、正常係数ベクトルαと病変係数ベクトルαｄとの間の距離が１．０であれば５７％の可能性で病変部が存在すると判定することができる。１．６は最大距離であり、病変判定部１０６は、距離が１．６であれば１００％病変部が存在すると判定することができる。
【００８６】
出力部１１１は、病変判定部１０６による判定結果を出力する。
【００８７】
表示部１０７は、ディスプレイ装置などにより構成され、出力部１１１より判定結果を受け、病変部が有りと判定された場合、病変部の画像座標の画素値を特定の色（たとえば赤や黄色など）で置き換えて、病変部の存在とその位置を画像として表示する。
【００８８】
図８は、診断支援装置１００の実行する処理のフローチャートである。
【００８９】
診断支援装置１００の処理に先立って、正常係数ベクトル格納部１０５に正常係数ベクトルαが格納される。つまり、複数の正常構造画像から取得したシフト不変特徴量ベクトル群から基底ベクトルが算出される。正常構造画像をこの基底ベクトルの線形和で記述した際の係数をベクトルの要素とする正常係数ベクトルαが作成される。作成された正常係数ベクトルαが正常係数ベクトル格納部１０５に格納される。
【００９０】
ステップＳ７１において、入力部１０１は、検査画像を受け付ける。
【００９１】
ステップＳ７２において、特徴量算出部１０２は、入力部１０１が受け付けた検査画像の画素ごとに、複数の画像特徴量を算出する。特徴量算出部１０２は、検査画像の画素ごとに、各画像特徴量をベクトルの要素とする検査画像特徴量ベクトルｆｐを作成し、基底表現部１０３に出力する。
【００９２】
ステップＳ７３において、基底表現部１０３は、検査画像の画素ごとに、検査画像特徴量ベクトルｆｐを検査係数ベクトルαｐへ変換する。
【００９３】
ステップＳ７４において、最近傍ベクトル検出部１０４は、正常係数ベクトル格納部１０５に格納された正常係数ベクトルαを参照して、検査画像の画素ごとに、検査係数ベクトルαｐに最も近い正常係数ベクトルαを検出する。
【００９４】
ステップＳ７５において、病変判定部１０６は、検査画像の画素ごとに、検査係数ベクトルαｐと最近傍の正常係数ベクトルαとの間の距離を算出し、算出した距離と病変の有無を判断する判定閾値との比較を行う。
【００９５】
検査係数ベクトルαｐと最近傍の正常係数ベクトルαとの間の距離が判定閾値より大きい場合は、ステップＳ７６へ進み、その処理対象画素については「病変部あり」と判定し、判手結果を出力部１１１に出力する。表示部１０７は、出力部１１１より判定結果を受け、処理対象画素の画素位置に画素値＝０（黒）を与えて、病変部を明示した差分画像を表示する。
【００９６】
一方、検査係数ベクトルαｐと最近傍の正常係数ベクトルαとの間の距離が判定閾値以下の場合は、ステップＳ７７へ進み、その処理対象画素については「病変部なし」と判定し、判定結果を出力部１１１に出力する。表示部１０７は、出力部１１１より判定結果を受け、処理対象画素の画素位置に画素値＝１（白）を与えて差分画像を表示する。
【００９７】
以上の構成により、診断支援装置１００は、正常構造画像と検査画像の違いから病変部を検出することができる。本実施の形態によると、シフト不変特徴量を利用することにより位置合わせ工程が不要となり、診断業務の効率と質が向上する。
【００９８】
（実施の形態２）
本実施の形態では、実施の形態１と同様に、画像情報のベクトル化をシフト不変特徴量を用いて実現し、位置合わせのための対応点設定が不要である診断支援装置について説明する。本実施の形態では、特に、類似症例検索サーバーと連携することで所見情報を参照し、病変部の類似症例を提示することができる診断支援装置について説明する。
【００９９】
図９は、本発明の実施の形態２に係る診断支援装置８００の機能的な構成を示すブロック図である。なお、図２と同様の構成要素については同じ符号を付し、詳細な説明を繰り返さない。
【０１００】
診断支援装置８００は、入力部１０１と、特徴量算出部１０２と、基底表現部１０３と、最近傍ベクトル検出部１０４と、正常係数ベクトル格納部１０５と、病変判定部８０１と、出力部８１１と、表示部８０２とを備える。
【０１０１】
病変判定部８０１は、実施の形態１の病変判定部１０６と同様に、検査画像の画素ごとに、検査係数ベクトルαｐと正常係数ベクトルαとの間の距離に基づいて、病変部の有り、無しを判定する。
【０１０２】
また、病変判定部８０１は、類似症例検索サーバー８０３に接続され、症例ＤＢ（データベース）部８０４に保存された診断済の症例データを類似症例検索サーバー８０３を介して参照する。過去の症例データは、医師が正式に下した診断結果である。診断支援装置８００は、医師または技師の診断を支援する位置付けにあり、診断の最終決定を行うことはない。過去の診断結果（過去症例データ）は診断支援装置８００から見れば「正解」であって、絶対的な参照情報である。一方、医師等も過去症例データを参照して新たな診断を下している。そこで、診断支援装置８００が病変の有無を提示するだけでなく、どのような病気であるかと、加えてその病状をも提示できれば、診断業務の効率と質が向上する。
【０１０３】
以上の考え方から、病変判定部８０１は、検査画像に病変ありと判定した場合には、類似症例検索サーバー８０３を介して症例ＤＢ部８０４を参照し、病変に類似する過去症例データを検出する。
【０１０４】
出力部８１１は、病変判定部８０１による判定結果と、病変判定部８０１が検出した過去症例データとを出力する。
【０１０５】
表示部８０２は、ディスプレイ装置などにより構成され、出力部８１１より判定結果と類似症例データとを受け、実施の形態１の表示部１０７と同様に、病変部の存在とその位置を画像として表示する。加えて、表示部８０２は、この病変部に類似する過去症例データを同時に表示する。
【０１０６】
図１０は、類似症例検索サーバー８０３による類似症例データの検索方法の一例を説明するための図である。
【０１０７】
病変判定部８０１は、検査画像９０２と正常構造画像９０３の違いを判定する。判定した結果は、差分画像９０１として表現される。ここで、差分画像９０１は、病変部を画素値＝０（黒）で、正常部を画素値＝１（白）で示した画像である。病変部９０４を拡大すると、分布パターン９０５が得られる。分布パターン９０５は、病変部の画素と、正常部の画素とから構成される。病変判定部８０１による病変部か否かの判定は画素単位で実行されるが、その結果を画像全体でまとめて見ると、複数の画素からなるパターンと捉えることができる。病名または病状の判断と分布パターン９０５との間には相関関係がある可能性が高い。このため、類似症例検索サーバー８０３は、分布パターン９０５を手がかりに、症例ＤＢ部８０４の中から、類似する分布パターン９０５を含む類似症例データを抽出することができる。
【０１０８】
図１０の場合、類似症例検索サーバー８０３は、症例ＤＢ部８０４に記憶されている症例データの中から、分布パターン９０５に最も類似する分布パターン９０６を含む画像９０７と、その画像９０７に付随する病名、病状、所見および臨床データなどを検出する。例えば、画像９０７に付随する病名はＡであり、病状はＸである。
【０１０９】
図１１は、表示部８０２による判定結果および類似症例データの表示例を示す図である。表示部８０２は、画像１２０１と、画像１２０３と、病名・病状１２０４と、所見１２０５とを表示する。画像１２０１は、検査画像９０２の病変部を線１２０２で囲んだ画像である。画像１２０３は、類似症例検索サーバー８０３が検出した、図１０に示した画像９０７と同じである。なお、症例ＤＢ部８０４には、画像１２０３のような症例画像が病名および病状ごとに格納されており、この分類に従って、病名はＡであり、病状はＸであることが表示される。各症例データには所見も付属しており、画像１２０３に対する所見１２０５が表示される。
【０１１０】
なお、本実施の形態では、分布パターン同士が類似しているか否かの判断方法および類似する分布パターンの探索方法については限定しない。例えば、分布パターン同士でパターンマッチングを行うことにより、分布パターン同士が類似しているか否かを判断すると共に類似する分布パターンの画像上での位置を特定しても良い。また、パターンマッチングの代わりに、分布パターンから画像特徴量を算出し、分布パターン同士で画像特徴量を比較するようにしても良い。
【０１１１】
以上の構成により、診断支援装置８００は、正常構造画像と検査画像の違いから病変部を検出することができる。本実施の形態によると、過去症例データを参照することにより、病変部と類似する症例データをユーザに提示することができる。このため、医師による病名および病状診断の業務効率および質を向上することができる。
【０１１２】
（実施の形態３）
本実施の形態では、実施の形態１および２と同様に、画像情報のベクトル化、をシフト不変特徴量を用いて実現し、位置合わせのための対応点設定が不要である診断支援装置について説明する。本実施の形態では、特に、入力装置、表示装置、電子カルテサーバーおよび診断支援サーバーを同一箇所に設置する必要がない。
【０１１３】
図１２は、本発明の実施の形態３に係る診断支援装置１０００の機能的な構成を示すブロック図である。なお、図２または図９と同様の構成要素については同じ符号を付し、詳細な説明を繰り返さない。
【０１１４】
診断支援装置１０００は、特徴量算出部１０２と、基底表現部１００４と、最近傍ベクトル検出部１０４と、正常係数ベクトル格納部１０５と、病変判定部８０１と、出力部８１１とを含む。
【０１１５】
診断支援装置１０００の特徴量算出部１０２は、ネットワークを介して入力装置１００１から検査画像を取得する。また、病変判定部８０１と基底表現部１００４は電子カルテサーバー１００２とネットワークを介して接続され、電子カルテサーバー１００２の症例ＤＢ部１００３に格納されている過去症例データを参照する。病変判定部８０１の出力である検査画像と正常構造画像との差分画像および類似症例データは、ネットワークを介して表示装置１００５へ送られる。
【０１１６】
入力装置１００１、電子カルテサーバー１００２、または表示装置１００５は、標準的なインターフェース（たとえば、ＤＩＣＯＭ（Ｄｉｇｉｔａｌ  Ｉｍａｇｉｎｇ  ａｎｄ  ＣＯｍｍｕｎｉｃａｔｉｏｎ  ｉｎ  Ｍｅｄｉｃｉｎｅ）やＨＬ７（Ｈｅａｌｔｈ  Ｌｅｖｅｌ  Ｓｅｖｅｎ）など）を介して、診断支援装置１０００と接続されている。これにより、診断支援装置１０００は、様々なモダリティから検査画像を受け付け、過去症例データを参照することにより、病変を検出することができる。これにより、診断支援装置１０００は、医師または技師による病変の検出を支援することができる。
【０１１７】
基底表現部１００４は、症例ＤＢ部１００３から正常構造画像を取得する。画像検査は病院内で毎日実施されるため、正常構造画像は毎日追加される。より多くの正常構造画像を取得して正常構造基底ベクトル行列Ｂと正常構造平均ベクトルｇを作成することにより、様々な体型、または撮影条件の違いに対応することができる。これにより、診断支援装置１０００の利用価値を高めることができる。このため、正常構造基底ベクトル行列Ｂと正常構造平均ベクトルｇは随時更新を行うことが望ましい。病変判定部８０１も電子カルテサーバー１００２とつなぐことで、症例ＤＢ部１００３の症例データ更新を反映でき、最新の症例データを含めた過去症例データの中から類似症例データを検出することができる。このため、診断支援装置１０００は、より的確な診断支援を実行することができる。
【０１１８】
以上の構成により、診断支援装置１０００は、入力装置１００１、表示装置１００５、電子カルテサーバー１００２、および診断支援サーバーを同一箇所に設置することなく、診断支援を行うことができる。そこで、たとえば、技師が検査画像を撮像し、遠隔地にいる画像診断医が表示装置１００５に表示された病変部の画像および病変部に類似する過去症例データを見ながら、診断を行うことができる。
【０１１９】
（実施の形態４）
本実施の形態では、実施の形態１〜３と同様に、画像情報のベクトル化を、シフト不変特徴量を用いて実現し、位置合わせのための対応点設定が不要である診断支援装置について説明する。本実施の形態では、特に、病変部である可能性を複数の提示方法で示し、より多面的な診断支援が可能である。
【０１２０】
図１３は、本発明の実施の形態４に係る診断支援装置の機能的な構成を示すブロック図である。なお、図２と同様の構成要素については同じ符号を付し、詳細な説明を繰り返さない。
【０１２１】
診断支援装置１３００は、入力部１０１と、特徴量算出部１０２と、基底表現部１０３と、近傍ベクトル検出部１３０１と、正常係数ベクトル格納部１０５と、病変判定部１３０２と、出力部１３０３と、表示部１０７とを含む。
【０１２２】
近傍ベクトル検出部１３０１は、正常係数ベクトル格納部１０５に格納された正常係数ベクトルαの中から、検査係数ベクトルαｐに類似するものを複数検出する。たとえば、検査係数ベクトルαｐとの間の距離が短いものから順に３本の正常係数ベクトルαを検出する。
【０１２３】
病変判定部１３０２は、近傍ベクトル検出部１３０１が検出した複数（例えば、３本）の正常係数ベクトルαを受け取り、正常係数ベクトルαごとに、検査係数ベクトルαｐと正常係数ベクトルαとの間の距離と判定閾値とを比較する。病変判定部１３０２は、比較結果に基づいて、実施の形態１の病変判定部１０６と同様に、検査係数ベクトルαｐが算出された着目画素が病変部の画素であるか否かを判定する。
【０１２４】
出力部１３０３は、病変判定部１０６による判定結果を、正常係数ベクトルαごとに分類して、表示部１０７に出力する。
【０１２５】
図１４は、表示部１０７に表示された判定結果の一例を示す図である。図１４に示すように、表示部１０７には、複数（例えば、３種類）の判定結果が表示される。例えば、検査画像である画像１２０３に対応する判定結果を示す画像１４０１〜１４０３が表示される。画像１２０３と画像１４０１〜１４０３とは同じ検査画像であるが、画像１４０１〜１４０３には、病変部１４０４〜１４０６が線で囲われて示されている。つまり、画像１４０１には、検査係数ベクトルαｐと検査係数ベクトルαｐに最も近い正常係数ベクトルαとの間の距離に基づいて判定された病変部１４０４が示されている。また、画像１４０２には、検査係数ベクトルαｐと検査係数ベクトルαｐに２番目に近い正常係数ベクトルαとの間の距離に基づいて判定された病変部１４０５が示されている。さらに、画像１４０３には、検査係数ベクトルαｐと検査係数ベクトルαｐに３番目に近い正常係数ベクトルαとの間の距離に基づいて判定された病変部１４０６が示されている。検査係数ベクトルαｐと正常係数ベクトルαとの間の距離が大きくなるほど病変部と判定される可能性が高くなる。このため、画像１４０１から画像１４０３の中で、画像１４０３の病変部１４０６の面積が最も広い。つまり、画像１４０３では、画像１４０１または画像１４０２に比べて緩く病変部が検出されていることに相当する。逆に、画像１４０１では、画像１４０２または画像１４０３に比べて厳しく病変部が検出されていることに相当する。
【０１２６】
病気や病状の医学的解釈は、日進月歩で進んでいる。特に、複数の医師および技師が参加して開催されるカンファレンスは、日々増える症例を元に、新たな診断方法や診断基準に関して議論が成される場となっている。このような検討を行う場において、図１４に示すような複数の病変判定基準に基づく判定結果を示すことは、新たな解釈を得るために貴重なデータとなる。
【０１２７】
以上説明したように、上述の実施の形態１〜４によると、画像情報のベクトル化をシフト不変特徴量で実現するため、位置合わせのための対応点設定が不要である。したがって、正常構造画像を安定に、効率良く作成できる。また、画像診断医の記憶に基づいた正常構造画像と検査画像の比較がデータ化され、画像診断医の診断を客観的に支援できる。また、診断結果はインフォームド・コンセントや医療教育、基礎医学など、様々な場面で活用できる。
【０１２８】
また、効率よく、かつ正確に病変部を検出することができる。画像診断医は、日々更新される過去症例データを参考にしながら新しい検査を行うため、計算機による病変検出は画像診断医の作業効率を高め、医療ワークフローを円滑にまわすことができる。病名および病状の確定は治療方針の決定に大きな影響を与えるため、診断支援装置による診断支援は医療全体の効率化、質的向上に貢献できる。
【０１２９】
以上、本発明の実施の形態に係る診断支援装置について説明したが、本発明は、これらの実施の形態に限定されるものではない。
【０１３０】
例えば、診断支援装置は、コンピュータにより実現することも可能である。
【０１３１】
図１５は、診断支援装置１００を実現するコンピュータシステムのハードウェア構成を示すブロック図である。診断支援装置８００および診断支援装置１０００も診断支援装置１００と同様にコンピュータシステムにより実現可能である。
【０１３２】
診断支援装置１００は、コンピュータ６０と、コンピュータ６０に指示を与えるためのキーボード６６およびマウス３８と、コンピュータ６０の演算結果等の情報を提示するためのディスプレイ６２と、コンピュータ６０で実行されるプログラムを読み取るためのＣＤ−ＲＯＭ（Ｃｏｍｐａｃｔ  Ｄｉｓｃ−Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）装置４０および通信モデム（図示せず）とを含む。
【０１３３】
診断支援装置１００が行う処理であるプログラムは、コンピュータで読取可能な記録媒体であるＣＤ−ＲＯＭ４２に記憶され、ＣＤ−ＲＯＭ装置４０で読み取られる。または、コンピュータネットワークを通じて通信モデム５２で読み取られる。
【０１３４】
コンピュータ６０は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）４４と、ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）４６と、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）４８と、ハードディスク５０と、通信モデム５２と、バス５４とを含む。
【０１３５】
ＣＰＵ４４は、ＣＤ−ＲＯＭ装置４０または通信モデム５２を介して読み取られたプログラムを実行する。ＲＯＭ４６は、コンピュータ６０の動作に必要なプログラムやデータを記憶する。ＲＡＭ４８は、プログラム実行時のパラメータなどのデータを記憶する。ハードディスク５０は、プログラムやデータなどを記憶する。通信モデム５２は、コンピュータネットワークを介して他のコンピュータとの通信を行う。バス５４は、ＣＰＵ４４、ＲＯＭ４６、ＲＡＭ４８、ハードディスク５０、通信モデム５２、ディスプレイ６２、キーボード６６、マウス３８およびＣＤ−ＲＯＭ装置４０を相互に接続する。
【０１３６】
さらに、上記の各装置を構成する構成要素の一部または全部は、１個のシステムＬＳＩ（Ｌａｒｇｅ  Ｓｃａｌｅ  Ｉｎｔｅｇｒａｔｉｏｎ：大規模集積回路）から構成されているとしても良い。システムＬＳＩは、複数の構成部を１個のチップ上に集積して製造された超多機能ＬＳＩであり、具体的には、マイクロプロセッサ、ＲＯＭ、ＲＡＭなどを含んで構成されるコンピュータシステムである。ＲＡＭには、コンピュータプログラムが記憶されている。マイクロプロセッサが、コンピュータプログラムに従って動作することにより、システムＬＳＩは、その機能を達成する。
【０１３７】
さらにまた、上記の各装置を構成する構成要素の一部または全部は、各装置に脱着可能なＩＣカードまたは単体のモジュールから構成されているとしても良い。ＩＣカードまたはモジュールは、マイクロプロセッサ、ＲＯＭ、ＲＡＭなどから構成されるコンピュータシステムである。ＩＣカードまたはモジュールは、上記の超多機能ＬＳＩを含むとしても良い。マイクロプロセッサが、コンピュータプログラムに従って動作することにより、ＩＣカードまたはモジュールは、その機能を達成する。このＩＣカードまたはこのモジュールは、耐タンパ性を有するとしても良い。
【０１３８】
また、本発明は、上記に示す方法であるとしても良い。また、これらの方法をコンピュータにより実現するコンピュータプログラムであるとしても良いし、前記コンピュータプログラムからなるデジタル信号であるとしても良い。
【０１３９】
さらに、本発明は、上記コンピュータプログラムまたは上記デジタル信号をコンピュータ読み取り可能な非一時的な記録媒体、例えば、フレキシブルディスク、ハードディスク、ＣＤ−ＲＯＭ、ＭＯ、ＤＶＤ、ＤＶＤ−ＲＯＭ、ＤＶＤ−ＲＡＭ、ＢＤ（Ｂｌｕ−ｒａｙ  Ｄｉｓｃ（登録商標））、半導体メモリなどに記録したものとしても良い。また、これらの非一時的な記録媒体に記録されている上記デジタル信号であるとしても良い。
【０１４０】
また、本発明は、上記コンピュータプログラムまたは上記デジタル信号を、電気通信回線、無線または有線通信回線、インターネットを代表とするネットワーク、データ放送等を経由して伝送するものとしても良い。
【０１４１】
また、本発明は、マイクロプロセッサとメモリを備えたコンピュータシステムであって、上記メモリは、上記コンピュータプログラムを記憶しており、上記マイクロプロセッサは、上記コンピュータプログラムに従って動作するとしても良い。
【０１４２】
また、上記プログラムまたは上記デジタル信号を上記非一時的な記録媒体に記録して移送することにより、または上記プログラムまたは上記デジタル信号を上記ネットワーク等を経由して移送することにより、独立した他のコンピュータシステムにより実施するとしても良い。
【０１４３】
さらに、上記実施の形態および上記変形例をそれぞれ組み合わせるとしても良い。
【０１４４】
なお、本発明に係る診断支援装置の必須の構成要素は、図１６に示すように、特徴量算出部１０２と、基底表現部１０３と、病変判定部１０６と、出力部１１１とである。各構成要素は、実施の形態１で説明したとおりである。これら以外の構成要素は備えられているのが望ましいが、本発明を実現する上では必ずしも必要はない。
【０１４５】
今回開示された実施の形態はすべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は上記した説明ではなくて請求の範囲によって示され、請求の範囲と均等の意味および範囲内でのすべての変更が含まれることが意図される。
【産業上の利用可能性】
【０１４６】
本発明は、医用画像から病変部を特定することにより、医師による画像診断を支援する診断支援装置等に適用可能である。
【符号の説明】
【０１４７】
３０  着目画素
３１  右隣接画素
３２  下隣接画素
３３  右斜下隣接画素
３４  ブロック
３５  右隣接ブロック
３６  下隣接ブロック
３７  右斜下隣接ブロック
３８  マウス
４０  ＣＤ−ＲＯＭ装置
４２  ＣＤ−ＲＯＭ
４４  ＣＰＵ
４６  ＲＯＭ
４８  ＲＡＭ
５０  ハードディスク
５２  通信モデム
５４  バス
６０  コンピュータ
６２  ディスプレイ
６６  キーボード
１００、８００、１３００  診断支援装置
１０１  入力部
１０２  特徴量算出部
１０３、１００４  基底表現部
１０４  最近傍ベクトル検出部
１０５  正常係数ベクトル格納部
１０６、８０１、１３０２  病変判定部
１０７、８０２  表示部
１０８  検出エラー通知信号
１０９  第１閾値データベース
１１０  第２閾値データベース
１１１、８１１、１３０３  出力部
５０１  ベクトル距離算出部
５０２  検出エラー閾値算出部
８０３  類似症例検索サーバー
８０４、１００３  症例ＤＢ部
９０１  差分画像
９０２  検査画像
９０３  正常構造画像
９０４、１２０２、１４０４〜１４０６  病変部
９０５、９０６  分布パターン
９０７、１２０１、１２０３、１４０１〜１４０３  画像
１００１  入力装置
１００２  電子カルテサーバー
１００５  表示装置
１２０４  病名・病状
１２０５  所見
１３０１  近傍ベクトル検出部
【図１】
0005607839.tif 000008
【図２】
0005607839.tif 000009
【図３】
0005607839.tif 000010
【図４】
0005607839.tif 000011
【図５】
0005607839.tif 000012
【図６】
0005607839.tif 000013
【図７】
0005607839.tif 000014
【図８】
0005607839.tif 000015
【図９】
0005607839.tif 000016
【図１０】
0005607839.tif 000017
【図１１】
0005607839.tif 000018
【図１２】
0005607839.tif 000019
【図１３】
0005607839.tif 000020
【図１４】
0005607839.tif 000021
【図１５】
0005607839.tif 000022
【図１６】
0005607839.tif 000023