(57)【要約】【課題】異常陰影候補の検出情報の表示数を減少させることにより、医用画像の読影効率を向上させる。【解決手段】画像処理装置２では、フィルタバンクを用いた検出処理１と、特徴量を用いた検出処理２とを直接的に行う処理パターン１又は２の異常陰影候補の検出処理を行う場合は全ての異常陰影候補を抽出し、上記処理１と処理２を並列的に行う処理パターン３の場合は偽陽性か真陽性かの判断が困難な候補又は視認性が低く、医師が発見しづらい候補を抽出する。そして、抽出された異常陰影候補についてのみ、その検出情報の表示出力を行う。【選択図】図１
2006346094.tif 000002
【特許請求の範囲】
【請求項１】
異なる検出アルゴリズムを複数有する異常陰影候補検出手段により、この複数の検出アルゴリズムによって医用画像から異常陰影候補を検出する検出工程と、
前記検出工程における検出過程に対して予め定められた異常陰影候補の抽出条件に従って、前記検出された異常陰影候補のうち、検出情報を示すべき異常陰影候補を抽出する抽出工程と、
前記抽出された異常陰影候補の検出情報を出力する出力工程と、
を含むことを特徴とする検出情報の出力方法。
【請求項２】
前記抽出条件は、前記複数の異なる検出アルゴリズムによる検出処理が行われる過程に応じて定められており、
前記検出工程では、前記複数の異なる検出アルゴリズムによる検出処理を組み合わせて検出を行い、
前記抽出工程では、前記検出工程における各検出アルゴリズムの検出処理が行われる過程に応じた抽出条件により、検出情報を出力すべき異常陰影候補を抽出することを特徴とする請求項１に記載の検出情報の出力方法。
【請求項３】
前記検出アルゴリズムは、検出対象の病変種に応じて複数準備されるとともに、前記抽出条件は、各病変種の検出アルゴリズムの検出処理の過程に応じて定められており、
前記検出工程では、検出対象の病変種に応じた検出アルゴリズムにより検出処理を行い、
前記抽出工程では、前記検出工程において行われた検出アルゴリズムの検出処理過程に応じた抽出条件により、検出情報を示すべき異常陰影候補を抽出することを特徴とする請求項１に記載の検出情報の出力方法。
【請求項４】
異なる検出アルゴリズムを複数有し、この複数の検出アルゴリズムにより医用画像から異常陰影候補を検出する異常陰影候補検出手段と、
検出過程に対して予め定められた異常陰影候補の抽出条件に従って、前記検出された異常陰影候補のうち、検出情報を示すべき異常陰影候補を抽出する抽出手段と、
前記抽出された異常陰影候補の検出情報を出力する出力手段と、
を備えることを特徴とする医用画像処理システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、医用画像から異常陰影候補を検出し、その検出情報を出力する出力方法及び医用画像処理システムに関する。
【背景技術】
【０００２】
医療の分野においては、ＣＴ（Computed Tomography）やＭＲＩ（Magnetic Resonance Imaging）等の撮影装置により患者を撮影した医用画像をデジタルデータに変換し、医師が診断を行う際には、表示ディスプレイ上に表示して読影を行うようになってきた。特に近年では、読影医の負担軽減や異常陰影の見落とし減少を目的として癌化部分等の異常陰影候補を検出するコンピュータ支援装置（Computer-Aided Diagnosis；以下、ＣＡＤという）が開発されてきた。
【０００３】
上記ＣＡＤにより異常陰影候補が検出されると、その検出情報として医用画像上に異常陰影候補の検出位置を指し示すマーカ（矢印や丸等）等を表示する方法が一般的である。その表示方法としては、（１）医用画像上に直接マーカを表示する方法（例えば、特許文献１参照）、（２）読影用にライフサイズ（等倍）の医用画像を主画像として表示し、その主画像の縮小画像上にマーカを表示した副画像を作成して主画像とともに表示する方法（例えば、特許文献２参照）が用いられている。
【特許文献１】特開２０００−２７６５８７号公報
【特許文献２】特開２００４−２３０００１号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
このような状況で、ＣＡＤにより検出された異常陰影候補の数が多い場合、上記（１）の方法では、読影医が観察する画像上にたくさんのマーカが表示されることとなり、医用画像自体が観察しにくくなってしまい、かえって読影作業に悪影響を及ぼすことがある。
【０００５】
一方、上記（２）の方法では、（１）の方法に比べ、主画像自体はマーカが表示されないため、医師の読影作業に与える影響が少ないが、縮小された副画像上にはやはりたくさんのマーカが表示されることとなり、見づらいという問題があった。また、主画像と副画像との位置の対応もとりにくくなる。
【０００６】
さらに、異常陰影候補の検出時には、その候補領域の画像特徴量を演算し、閾値と比較して異常陰影であるか否かを判定する方法や、画像特徴量を用いて多変量解析することにより異常陰影の候補領域を検出する方法等、様々なアルゴリズムが用いられるが、検出結果には、明らかに異常陰影であるものと疑わしいものとが混在している。
【０００７】
従来のＣＡＤでは、このように検出された異常陰影候補の全てに関する検出情報を医師に提供し、真陽性であるか偽陽性であるかの最終判断は医師によることとしている。しかし、上述したように検出数が非常に多いと、その全てについて読影医は確認しなければならず、非効率的である。特に、経験豊かな医師にとっては、明らかに異常陰影であると容易に判断できるものより、疑わしい候補を優先的に読影した方が効率的な場合がある。
【０００８】
本発明の課題は、異常陰影候補の検出情報の表示数を減少させることにより、医用画像の読影効率を向上させることである。
【課題を解決するための手段】
【０００９】
請求項１に記載の発明は、検出情報の出力方法において、
異なる検出アルゴリズムを複数有する異常陰影候補検出手段により、この複数の検出アルゴリズムによって医用画像から異常陰影候補を検出する検出工程と、
前記検出工程における検出過程に対して予め定められた異常陰影候補の抽出条件に従って、前記検出された異常陰影候補のうち、検出情報を示すべき異常陰影候補を抽出する抽出工程と、
前記抽出された異常陰影候補の検出情報を出力する出力工程と、
を含むことを特徴とする。
【００１０】
請求項２に記載の発明は、請求項１に記載の検出情報の出力方法において、
前記抽出条件は、前記複数の異なる検出アルゴリズムによる検出処理が行われる過程に応じて定められており、
前記検出工程では、前記複数の異なる検出アルゴリズムによる検出処理を組み合わせて検出を行い、
前記抽出工程では、前記検出工程における各検出アルゴリズムの検出処理が行われる過程に応じた抽出条件により、検出情報を出力すべき異常陰影候補を抽出することを特徴とする。
【００１１】
請求項３に記載の発明は、請求項１に記載の検出情報の出力方法において、
前記検出アルゴリズムは、検出対象の病変種に応じて複数準備されるとともに、前記抽出条件は、各病変種の検出アルゴリズムの検出処理の過程に応じて定められており、
前記検出工程では、検出対象の病変種に応じた検出アルゴリズムにより検出処理を行い、
前記抽出工程では、前記検出工程において行われた検出アルゴリズムの検出処理過程に応じた抽出条件により、検出情報を示すべき異常陰影候補を抽出することを特徴とする。
【００１２】
請求項４に記載の発明は、医用画像処理システムにおいて、
異なる検出アルゴリズムを複数有し、この複数の検出アルゴリズムにより医用画像から異常陰影候補を検出する異常陰影候補検出手段と、
検出過程に対して予め定められた異常陰影候補の抽出条件に従って、前記検出された異常陰影候補のうち、検出情報を示すべき異常陰影候補を抽出する抽出手段と、
前記抽出された異常陰影候補の検出情報を出力する出力手段と、
を備えることを特徴とする。
【発明の効果】
【００１３】
請求項１、２、４に記載の発明によれば、検出処理の組み合わせによって異なる検出結果に応じて、医師に検出情報を提供する異常陰影候補を制限することができる。これにより、検出情報を出力する異常陰影候補数を減少させることができ、医師が読影しやすい出力形態とすることができる。従って、読影効率を向上させることができる。
【００１４】
請求項３に記載の発明によれば、検出目的とする病変種によって異なる検出結果に応じて、医師に医師に検出情報を提供する異常陰影候補を制限することができる。これにより、検出情報を出力する異常陰影候補数を減少させることができ、医師が読影しやすい出力形態とすることができる。従って、読影効率を向上させることができる。
【発明を実施するための最良の形態】
【００１５】
〈第１実施形態〉
本実施形態では、複数の検出アルゴリズムにより異なる検出過程を経て異常陰影候補の検出を行う際に、検出された異常陰影候補からその検出過程に応じて表示出力対象の異常陰影候補を抽出し、抽出された異常陰影候補の検出情報を表示出力する例を説明する。
【００１６】
まず、構成を説明する。
図１に、本実施形態における医用画像処理システム１００のシステム構成を示す。
医用画像処理システム１００は、医用画像を撮影し、当該医用画像に画像処理を行って異常陰影候補を検出し、医用画像とともにその検出情報を医師に提供するシステムである。
【００１７】
図１に示すように、医用画像処理システム１００は、画像生成装置１、画像処理装置２、プリンタ３、画像サーバ４、ビューア５を備えて構成されている。これら各装置１〜５は、ＬＡＮ（Local Area Network）等の医療機関内で構築された通信ネットワークＮを介して相互にデータを送受信可能に接続されている。通信ネットワークＮは、ＤＩＣＯＭ（Digital Imaging and Communication in Medicine）規格が適用されている。
【００１８】
以下、各構成装置１〜５について説明する。
画像生成装置１は、人体を撮影し、その撮影画像（医用画像）のデジタルデータを生成するものであり、例えばＣＲ（Computed Radiography）、ＦＰＤ（Flat Panel Detector）、ＣＴ、ＭＲＩ、カセッテ専用の読取装置、フィルムディジタイザ等のモダリティを適用可能である。本実施形態では、画像生成装置１として、左右乳房のＸ線撮影を行う乳房専用のＣＲを適用し、乳房画像のデータが生成されるものとする。
【００１９】
なお、画像生成装置１は、上述したＤＩＣＯＭ規格に準拠した装置であり、生成した医用画像に付帯させる各種情報（例えば、医用画像が撮影された患者に関する患者情報や、撮影又は検査に関する撮影情報、検査情報等。）を外部から入力可能であるとともに、自動生成することもできる。画像生成装置１は、生成された医用画像に上記付帯情報をヘッダ情報として付加して通信ネットワークＮを介して画像処理装置２へ送信する。なお、ＤＩＣＯＭ規格に準拠していない場合には、図示しないＤＩＣＯＭ変換装置を用いて付帯情報を画像生成装置１に入力させることも可能である。
【００２０】
画像処理装置２は、画像生成装置１から供給される医用画像に対して各種画像処理を施すとともに、当該医用画像の画像解析を行って異常陰影候補の検出処理を行う。
【００２１】
プリンタ３は、画像処理装置２又は画像サーバ４から受信された医用画像のデータに基づいて、フィルム等の記録媒体に医用画像を出力する。
【００２２】
画像サーバ４は、画像ＤＢ４ａを備え、この画像ＤＢ４ａに、画像生成装置１により生成された医用画像（原画像）、画像処理装置２から受信された画像処理済みの医用画像（処理画像）を保存し、その入出力を管理する。
【００２３】
ビューア５は、診断用に医師により使用される出力手段であり、ＬＣＤ（Liquid Crystal Display）等を備えて構成されている。ビューア５は、医師の操作指示に応じて、指定された医用画像を画像サーバ４から取得して表示する。或いは、画像処理装置２で検出された異常陰影候補の検出結果を受信して表示する。
【００２４】
次に、本発明に係る画像処理装置２について詳細に説明する。
図２に、画像処理装置２の内部構成を示す。
画像処理装置２は、制御部２１、操作部２２、表示部２３、通信部２４、記憶部２５、画像処理部２６、異常陰影候補検出部２７を備えて構成されている。
【００２５】
制御部２１は、ＣＰＵ（Central Processing Unit）、ＲＡＭ（Random Access Memory）等から構成され、ＣＰＵにより記憶部２５から各種制御プログラムを読み出してＲＡＭに展開する。そして、当該プログラムに従って処理の実行を統括的に制御し、各部の動作を集中制御する。例えば、本発明に係る表示対象判断処理プログラムに従って表示対象判断処理（後述する）を実行する。この制御部２１と表示対象判断処理プログラムとの協働により、抽出手段を実現することができる。
【００２６】
操作部２２は、キーボードやマウス等を備えて構成され、キー操作やマウス操作に対応する操作信号を生成して制御部２１に出力する。
【００２７】
表示部２３は、ＬＣＤ等を備えて構成され、制御部２１の指示に従って画像処理時の操作画面や処理後の医用画像等、各種表示を行う。
【００２８】
通信部２４は、ルータやモデム等の通信用インターフェイスを備え、制御部２１の指示に従って、通信ネットワークＮ上の外部装置と通信を行う。例えば、画像生成装置１から処理対象の医用画像を受信したり、処理済みの医用画像を画像サーバ４やプリンタ３に送信する。
【００２９】
記憶部２５は、各種制御プログラムやプログラムの実行に必要なパラメータ、条件の情報或いは処理結果のデータを記憶している。
【００３０】
画像処理部２６は、画像処理プログラムに従って医用画像に対し、階調変換処理、鮮鋭度調整処理等を実行する他、乳房画像の場合には左右の乳房の画像を隣り合わせに合成する位置合わせ処理等を実行する。
【００３１】
異常陰影候補検出部２７は、医用画像を画像解析し、異常陰影の候補領域を検出する異常陰影候補検出手段である。異常陰影候補検出部２７は、複数の検出アルゴリズムを有しており、これら複数の検出アルゴリズムを用いて検出処理を実行する。本実施形態では、フィルタバンクを用いた検出アルゴリズム、画像特徴量を用いた検出アルゴリズムの２つを有している例を説明する。具体的には、記憶部２５に記憶されている各種アルゴリズムによる検出処理プログラムを読み出し、当該プログラムとの協働により検出処理を実行する。なお、検出方法の詳細については、後述する動作説明とともに説明する。
【００３２】
次に、本実施形態における医用画像処理システム１００の動作について説明する。
最初に、医用画像の生成から保存までの流れを説明する。
まず、画像生成装置１において撮影が行われ、医用画像（ここでは、乳房画像の例を説明する）が生成される。そして、生成された乳房画像の関連情報として、患者の指名、年齢、性別等の患者に関する情報や、撮影時の管電圧、乳房の圧迫厚等の撮影情報、検査日時等の検査情報、画像の読取条件等の乳房画像生成に関する詳細情報等が乳房画像に付帯される。
【００３３】
付帯情報が付帯された乳房画像は、画像生成装置１から画像処理装置２に出力される。
画像処理装置２では、通信部２４を介して入力された乳房画像が記憶部２５に記憶され、その後画像処理部２５により当該乳房画像に必要な画像処理が施される一方で、その乳房画像について異常陰影候補検出部２７により異常陰影候補の検出処理が開始される。
【００３４】
異常陰影候補の検出処理は、２つの検出アルゴリズムを組み合わせて行われるが、その処理パターンとして、２つの検出アルゴリズムの組み合わせて検出過程を異なるものとした３態様があり、何れの処理パターンを適用するかは技師又は読影医が設定可能である。
処理パターンは第１〜第３のパターンがあり、第１の処理パターンは一の検出アルゴリズムにより検出処理を行った後、その検出処理により偽陽性と判別された領域を除く画像領域に対し、他の検出アルゴリズムにより検出処理を行うものである。また、第２の処理パターンは第１の処理パターンにおける各検出アルゴリズムの処理順番を逆にした点が異なるだけのものであり、第３の処理パターンは一の検出アルゴリズムと他の検出アルゴリズムの検出処理を並行して行い、何れかの検出アルゴリズムで検出されれば全て検出結果とするものである。
【００３５】
以下、第１の処理パターンについて図３を参照して説明した後、他の処理パターン（第２の処理パターン、第３の処理パターン）について説明する。
図３は、第１の処理パターンにおける異常陰影候補の検出処理の流れを示すフローチャートである。
【００３６】
図３に示すように、まず異常陰影候補検出部２７において、記憶部２４に保存された乳房画像が読み出され（ステップＳ１）、当該乳房画像に対して、検出アルゴリズムとしてフィルタバンクを用いた異常陰影候補検出処理（以下、処理１という）が行われる（ステップＳ２）。
【００３７】
ここで、円形又は線状パターンを検出対象とするフィルタバンクの検出アルゴリズムについて説明する（電子情報通信学会論文誌  D-II Vol.J87-D-II No.1 pp.186-196参照）。
図４に、ウェーブレット解析における多重解像度の概念を取り入れたフィルタバンクを示す。原画像をS0fとし、この解像度レベルをゼロとする。HL(zj),FL(zj)は解像度レベルｊにおける低域通過フィルタ、HH(zj),FH(zj)は高域通過フィルタである。また、Sjf,Whjf,Wvjf,wdjfは、それぞれ解像度レベルｊの平滑化部分画像、水平部分画像、垂直部分画像、対角部分画像である。フィルタバンクの帯域分割フィルタバンク側では低域側成分の分割を繰り返す。帯域合成フィルタバンク側では、逆にそれらを合成する。このフィルタバンクにおいて、各フィルタは格式（１）〜（４）で与えられる。
【数１】
2006346094.tif 000003
【数２】
2006346094.tif 000004
【数３】
2006346094.tif 000005
【数４】
2006346094.tif 000006
【００３８】
ここで、HH(z)は１階差分フィルタ、HH(z)FH(z)は２階差分フィルタであるから、解像度レベルｊにおける水平部分画像、垂直部分画像は、解像度レベルｊ−１の平滑化部分画像に対してそれぞれ、垂直方向に２階差分フィルタを加えた画像、水平方向に２階差分フィルタを加えた画像に相当する。また、対角部分画像は、水平方向及び垂直方向に１階差分フィルタを加えた画像となる。すなわち、これらの部分画像は下記式（５）で示すヘッセ行列Hの各要素に相当する。
【数５】
2006346094.tif 000007
【００３９】
解像度レベルｊの水平部分画像、垂直部分画像、対角部分画像を各要素とするヘッセ行列の最小固有値を解像度レベルｊの円形パターン強調画像λminj(x,y)として用いる（下記式（６））。
【数６】
2006346094.tif 000008
【００４０】
また、ヘッセ行列の最大固有値を解像度レベルｊの円形・線状パターン強調画像λmaxj(x,y)として用いる（下記式（７））。
【数７】
2006346094.tif 000009
【００４１】
円形パターン強調画像λminj(x,y)、円形・線状パターンλmaxj(x,y)を用いることにより、異なる大きさの円形パターン及び円形・線状パターンの検出が可能となる。
【００４２】
図５を参照し、上記説明したフィルタバンクの検出アルゴリズムによる異常陰影候補の検出処理について説明する。
まず、検出対象の医用画像がフィルタバンクに入力される（ステップＳ２１）。
次に、フィルタバンクにより得られた複数の解像度の画像から、ヘッセ行列の最大固有値及び最小固有値が算出される（ステップＳ２２）。そして、異なる大きさの円形・線状パターン強調画像が作成され（ステップＳ２３）、円形成分及び線状成分が抽出される。
【００４３】
次に、強調画像に関心領域（例えば、５ｍｍ×５ｍｍ）が設定され、その領域にある円形成分・線状成分の画像特徴量が算出される（ステップＳ２４）。円形成分・線状成分の画像特徴量には、各成分の大きさ、各成分の分布等が含まれる。
そして、円形成分・線状成分の画等特徴量に基づいて、円形成分・線状成分が真陽性又は偽陽性の何れに該当するかが判別される（ステップＳ２５）。
【００４４】
判別は、マハラノビス距離等の判別分析を用いることができる。これは、予め真陽性及び偽陽性と分かっている学習データを準備しておき、これら学習データの画像特徴量を元にマハラノビス距離を求め、真陽性と偽陽性の何れに特徴が近いかを判別するものである。判別結果が得られると、図３のステップＳ３の処理へ移行する。
ステップＳ３では、処理１における異常陰影候補の検出結果として、真陽性及び偽陽性と判別された候補領域の位置、面積等が算出され、記憶部２５に記憶される。
【００４５】
次に、記憶部２５から処理１で検出対象とした医用画像と同一の画像が読み出される（ステップＳ４）。また、記憶部２５から処理１により検出された異常陰影候補の検出結果の情報が読み出され、当該読み出された検出結果（位置情報）に基づいて前記医用画像から処理１で偽陽性と判別された画像領域を除いた画像（これを前処理画像という）を生成する前処理が行われる（ステップＳ５）。具体的には、偽陽性と判別された画像領域の各画素に検出処理の処理対象から除外することを示すフラグを設定する等の処理を行う。
【００４６】
次に、生成された前処理画像を対象として、画像特徴量を用いた異常陰影候補検出処理（以下、処理２という）が行われる（ステップＳ６）。処理２は、画像に含まれる陰影領域の画像特徴量を算出し、算出された画像特徴量に基づいて陰影が真陽性か偽陽性かを判別して異常陰影候補の検出を行う検出アルゴリズムによる検出処理である。
【００４７】
図６を参照して、処理２の検出アルゴリズムによる検出処理の流れについて説明する。
まず、前処理画像から異常陰影の候補領域が特定される（ステップＳ３１）。このときの方法としては、特開平１０−９１７５８号公報に開示されているアイリスフィルタを用いた手法や、ラプラシアンフィルタを用いた手法（電気情報通信学会論文誌(D-II),Vol.J76-D-II,no.2,pp.241-249,1993）、モルフォルジーフィルタ（電気情報通信学会論文誌(D-II),Vol.J71-D-II,no.7,pp.1170-1176,1992）、ラプラシアンフィルタ（電気情報通信学会論文誌(D-II),Vol.J71-D-II,no.10,pp.1994-2001,1998）、３重リングフィルタ等、病変種に応じて種々のアルゴリズムが開発されているので、何れのものを適用してもよい。ここでは、一例として、微小石灰化クラスタ陰影の候補領域を特定する際に用いられる３重リングフィルタについて説明する。
【００４８】
微小石灰化クラスタ陰影は、乳房画像上では略円錐形状の濃度変化を有する低濃度の微小石灰化部分が集まった（クラスタ化した）陰影として現れる。このような濃度特性に基づいて医用画像に対して正方形の領域を順次設定し、この領域（これを注目領域という）毎に、微小石灰化部分を検出するフィルタとして、特有のベクトルパターンを有する３重リングフィルタを適用して異常陰影候補領域の検出を行う。なお、この注目領域のサイズは、検出目的とする病変種に応じて設定すればよい。
【００４９】
３重リングフィルタは、濃度変化が理想的な円錐形状を示す場合の濃度勾配の強度成分及び方向成分が予め決定された３つのリングフィルタから構成される。まず、ある注目画素周辺において、各リングフィルタのそれぞれの領域上の画素値から濃度勾配の強度成分及び方向成分の代表値を求める。そして、その代表値と各リングフィルタに予め決定されている濃度勾配の強度成分及び方向成分との差に基づいて、円錐形状に近い濃度変化を有する画像領域を候補領域として検出する。
なお、前処理により偽陽性の領域として削除された領域については、フィルタが適用されず、候補領域の特定は行われない。
【００５０】
このようにして、候補領域が特定されると、当該特定された候補領域の画像特徴量が算出される（ステップＳ３２）。画像特徴量には、候補領域の面積、円形度、不整形度、平均画素値、標準偏差、候補領域とそのコントラスト等が含まれる。
【００５１】
ここで、面積は特定された候補領域を構成する画素数で示される。
また、円形度は形状の複雑さを示す画像特徴量であって、候補領域の面積Ｓ、候補領域の輪郭線の長さ（周囲長）をＬとして、下記式（８）で示される。
e=4πS/L2・・・（８）
【００５２】
或いは、候補領域の重心を中心んとした、候補領域の面積Ｓと同一免疫の円と候補領域とが重なった面積（重なり面積）をUとして、下記式（９）により示すこともできる。
e=U/S・・・（９）
円形度ｅは、円に近いほど大きくなり、その値は１に近づく。
【００５３】
不整形度fは、候補領域の周囲長をL、候補領域の周囲長を平滑化した後の周囲長をL′として、下記式（１０）により示される。
f=L′/L・・・（１０）
【００５４】
次に、算出された各画像特徴量に基づいて、処理２により検出された異常陰影候補が真陽性か偽陽性かが判別される（ステップＳ３３）。判別は、例えば予め真陽性又は偽陽性の２つのグループに属する陰影画像の画像特徴量を学習データとして準備しておき、各グループの中心から、判別対象データ（算出された画像特徴量）までのマハラノビス距離を求め、真陽性又は偽陽性の何れのグループに近いかを判別することにより行うことができる。判別結果が得られると、図３に示すステップ７に移行する。
ステップＳ７では、処理２における検出結果として、異常陰影候補の位置情報、面積等の情報が記憶部２５に記憶される（ステップＳ７）。
【００５５】
次いで、制御部２１により記憶部２５から処理パターン１による検出結果の情報とともに処理パターン１に対して設定されている抽出条件の情報が読み出される。なお、処理パターン１に対して設定されている抽出条件は、全ての異常陰影候補を抽出することである。すなわち、読み出された検出結果の情報に基づいて、検出された全ての異常陰影候補が表示対象の候補として抽出され、その検出情報が表示部２３上に表示される（ステップＳ８）。検出情報とは、異常陰影候補の検出に関して医師に提供する情報をいい、例えば検出位置を指し示すマーカ画像や、検出対象とした病変種、面積等、異常陰影候補に関する事項を示す文字情報等が挙げられる。
【００５６】
以上が、第１の処理パターンにおける処理の流れであるが、第２の処理パターンは第１の処理パターンにおける処理１と処理２の順番が逆に入れ替えられただけであるので、第２の処理パターンについては、詳細な説明は省略する。第２の処理パターンにおける異常陰影候補の抽出条件も第１の処理パターンの条件と同一である。すなわち、処理２により偽陽性の領域が削除された前処理画像について処理１が実行され、処理１により真陽性と判別された異常陰影候補の全てについてその検出情報が表示されることとなる。
【００５７】
次に、第３の処理パターンに係る異常陰影候補検出処理について、図７を参照して説明する。
図７に示すように、まず異常陰影候補検出部２７において、記憶部２４に保存された乳房画像が読み出される（ステップＴ１）。次いで、この乳房画像に対し、処理１及び処理２による検出処理がそれぞれ独立して並行で行われ（ステップＴ２１、Ｔ２２）、各検出処理による検出結果がそれぞれ記憶部２５に記憶される（ステップＴ３１、Ｔ３２）。なお、各処理１、２の内容については上述した通りであるので、ここでは説明を省略する。
【００５８】
次いで、制御部２１において、処理１及び処理２による検出結果が記憶部２５から読み出され、これら検出された異常陰影候補のうち、その検出情報を表示出力すべき異常陰影候補が抽出される（ステップＴ３３）。
【００５９】
第１又は第２の処理パターンの場合、一の検出アルゴリズムで真陽性と判別された領域について他の検出アルゴリズムでもう一度検出処理を行い、他の検出アルゴリズムでもやはり真陽性であると判別された領域を最終的な検出結果として出力している。つまり、２段階で検出処理を行っているため、検出結果に含まれる異常陰影候補は真陽性である可能性は高いものと予想される。
【００６０】
一方、第３の処理パターンの場合、一の検出アルゴリズムと他の検出アルゴリズムとを別個に並行して行い、何れかの検出アルゴリズムで検出された異常陰影候補が最終的な検出結果として出力される。そのため、第１又は第２の処理パターンに比べ、その検出精度は低いと考えられ、その分検出結果には偽陽性候補が含まれている確率が高いと思われる。このような検出結果に基づいて第１の処理パターンと同様に検出された異常陰影候補の全てについてその検出情報を表示することとすると、表示部２３におけるその表示数が増え、見づらい表示となってしまう。
よって、第３の処理パターンでは、検出された異常陰影候補のうち、表示対象とする候補を抽出することにより、その検出情報の表示数を減少させることとする。
【００６１】
ここで、検出情報を表示すべき異常陰影候補は、その異常陰影候補が真陽性又は偽陽性とも明確な判別がしづらく、医師の最終判断に委ねた方がよいものとし、そのような候補を抽出するよう抽出条件が設定されて記憶部２５に記憶されているものとする。明らかに目視で真陽性と分かる陰影については特に画像処理装置２による検出情報の表示が無くとも、医師が発見することができる。そのため、医師によっては、そのような陰影についてではなく、むしろ真陽性か偽陽性かの判断が困難な陰影についての検出情報を参照したいという要望がある。よって、このような要望に応じて抽出条件を設定する。
また、その検出位置や画像上の視認性が低く、医師が発見しづらい異常陰影候補の検出情報は提供すべきであるのでこのような異常陰影候補を抽出することとしてもよい。
【００６２】
表示対象とすべき異常陰影候補の抽出は、画像特徴量を用いて行われる。
例えば、コントラストが小さい、面積が小さい場合は、真陽性か偽陽性かの判定が難しく、読影医による最終判断を必要とする場合が多い。また、視認性についても良好とは言い難く、医師が見落としやすい陰影となる。また、真陽性の異常陰影は、辺縁部の境界が不明瞭（不鮮明）となったり、辺縁部において微細なスジ状の形状を呈する、或いはスピキュラと呼ばれる歪（いびつ）な形状を有する場合があり、その辺縁の不整形度が増すが、逆に辺縁の不整形度が小さい場合には、真陽性かどうかの判断が難しくなる。また、腫瘤の場合、円形度が大きいほど、腫瘤である可能性が高くなるが、円形度が小さいと腫瘤（真陽性）か否かの判断がしづらくなる。
よって、これらの画像特徴量を判断要素として、真陽性か否かの判断が困難である、或いは視認性が低い異常陰影候補の判断を行う。
【００６３】
具体的には、候補領域において算出されたコントラスト、標準偏差、濃度平均値、曲率、面積、フラクタル次元、円形度、辺縁の不整形度等の各種画像特徴量に基づいて、多変量解析により真陽性の可能性を示す指標値が算出される。例えば、予め真陽性と分かっている陰影から算出した画像特徴量を学習データとしてニューラルネットワーク等の多変量解析による判定ロジックを構築しておき、この判定ロジックに判定対象の陰影候補から算出された各種画像特徴量を入力して当該陰影候補の真陽性の可能性を示す指標値を得る。
【００６４】
そして、この指標値と、偽陽性候補を検出するための閾値ＴＨ１、真陽性候補を検出するための閾値ＴＨ２とを比較した結果に応じて、検出された異常陰影候補において偽陽性候補及び真陽性候補を判断し、これらを除く他の異常陰影候補を表示すべき異常陰影候補として抽出する。
【００６５】
例えば、閾値ＴＨ１＝０．２５、ＴＨ２＝０．８５であるとし、ＴＨ１より低い指標値を有する候補を偽陽性候補、ＴＨ２より高い指標値を有する候補を真陽性候補と判断する条件の場合について説明する。検出時に出力された真陽性に対する指標値（これをＳとする）が０〜１の値（０に近いほど真陽性である可能性が低く、１に近いほど真陽性である可能性が高い。）で出力されるとすると、検出候補のうち、指標値ＳがＴＨ１＜Ｓ＜ＴＨ２である異常陰影候補が抽出される。
【００６６】
すなわち、抽出された候補は、検出処理によって検出された異常陰影候補の中から、偽陽性候補及び真陽性候補を除く候補、つまり真陽性か否かが疑わしく、医師の判断が必要な候補が抽出された候補となる。
さらに、画像端部に近い乳房領域に存在する異常陰影は、その端部という位置のために医師にとって見逃しやすいものとなる。よって、このような医師が発見しづらい位置にある異常陰影候補を表示対象に加えたい場合は、予め発見が困難な領域の境界として画像端からの距離を求めてこれを閾値ＴＨ３とし、検出された異常陰影候補のうち、閾値ＴＨ３を下回る（つまり、画像端に近い位置にある）異常陰影候補を表示すべき異常陰影候補として抽出することとしてもよい。
【００６７】
なお、上記偽陽性候補及び真陽性候補を判断するために用いられる閾値ＴＨ１、ＴＨ２は、読影毎に設定可能な構成とする。具体的には、操作部２２を介して閾値ＴＨ１、ＴＨ２の設定情報を入力可能な構成としておき、当該閾値ＴＨ１、ＴＨ２の設定情報を各読影医に固有の読影医ＩＤ等と対応付けて記憶部２５に格納しておく。そして、読影時に読影医ＩＤ等を入力させることにより、その読影医に対応する閾値ＴＨ１、ＴＨ２を制御部２１が読み出し、上記処理に用いる。
【００６８】
例えば、なるべく偽陽性候補は表示させたくない医師であれば、閾値ＴＨ１を０．２５から０．３に設定変更する等、医師の技量や読影スタイル等に応じて閾値ＴＨ１又はＴＨ２を設定することにより、異常陰影候補の表示数を調整することができ、各読影医のスキル等に相応しい表示とすることができる。
【００６９】
以上のようにして、表示すべき異常陰影候補が抽出されると、当該抽出された異常陰影候補の検出情報が表示部２３上に表示される（ステップＴ３４）。
【００７０】
図８に、各処理パターン１〜３による異常陰影候補の検出情報の表示例を示す。
図８（ａ）は処理パターン１又は２の表示例を示す図であり、図８（ｂ）は処理パターン３の表示例を示す図である。図８（ａ）、（ｂ）では、乳房画像上に異常陰影候補の検出情報として、その検出位置を示すマーカ画像（図中の矢印）が表示されている。
【００７１】
処理パターン１又は２では、真陽性の可能性が高い候補のみが検出結果として出力されるため、図８（ａ）に示すように異常陰影の可能性が高い陰影に絞ってマーカを表示することができる。
一方、処理パターン３では、真陽性又は偽陽性の可能性が高い候補を削除した残りの候補のみが検出結果として出力されるため、形状や濃度等からは真陽性か偽陽性か判断が困難な陰影に絞ってマーカを表示することができる。
【００７２】
以上のように、本実施形態によれば、２つのの検出処理の組み合わせで処理パターン１〜３の異なる検出過程を経る場合に、その検出過程に応じて検出された異常陰影候補からその検出情報を表示すべき候補を抽出し、抽出された異常陰影候補についてのみ、検出情報を表示する。これにより、検出過程に応じて検出情報の表示数を制限して、読影しやすい表示形態とすることができ、医師の読影効率を向上させることができる。
【００７３】
また、真陽性候補の検出率が高い処理パターン１又は２の場合は全ての異常陰影候補を抽出し、偽陽性候補の検出率が高い処理パターン３の場合は偽陽性か真陽性かの判断が困難な候補又は視認性が低く、医師が発見しづらい候補を抽出する。このように検出過程により生じる検出率の違いに応じて、その異常陰影候補の検出情報の表示数を可変させるので、表示数が減少するように調整することができる。
【００７４】
〈第２実施形態〉
第２実施形態では、検出目的とする病変種に対応した複数の検出アルゴリズムを有し、検出する病変種に応じた検出アルゴリズムを用いて異常陰影候補の検出を行う際に、用いる検出アルゴリズムによる異なる検出過程に応じて表示対象の異常陰影候補を抽出し、抽出された異常陰影候補の検出情報を表示出力する例を説明する。
【００７５】
第２実施形態における医用画像処理システムは、第１実施形態に示す医用画像処理システム１００と同一構成であり、画像処理装置２における機能及び処理動作が異なるのみである。よって、第１実施形態と同一の符号を付し、画像処理装置２における異なる機能部分及び処理動作の説明のみ行う。
【００７６】
第２実施形態において、画像処理装置２の記憶部２５には、検出目的とする病変種に応じた検出アルゴリズムによる検出処理プログラムが格納されている。本実施形態では、乳房における病変である腫瘤、微小石灰化クラスタの２種類に対応した検出アルゴリズムのプログラムが格納されていることとする。
【００７７】
次に、第２実施形態における画像処理装置２の異常陰影候補の検出処理について、図９を参照して説明する。
図９に示す異常陰影候補の検出処理では、まず異常陰影候補検出部２７により記憶部２５から処理対象の医用画像が読み出される（ステップＥ１）。次いで、検出対象の病変種を選択するための選択画面が表示部２３に表示され、腫瘤又は微小石灰化クラスタのうち、何れの病変種の異常陰影を検出対象とするかを選択操作するようガイダンスが行われる。
【００７８】
そして、この選択画面において操作部２２を介して医師により何れか一方の病変種の選択操作がなされると（ステップＥ２）、微小石灰化クラスタが選択操作された場合は（ステップＥ２；クラスタ）、ステップＥ３の処理へ移行し、腫瘤が選択操作された場合は（ステップＥ２；腫瘤）、ステップＥ６の処理へ移行する。
【００７９】
まず、微小石灰化クラスタが選択された場合について説明する。
ステップＥ３では、異常陰影候補検出部２７において微小石灰化クラスタに応じた検出処理プログラムが記憶部２５から読み出され、当該プログラムとの協働により異常陰影候補の検出処理が行われる。すなわち、微小石灰化クラスタの検出に応じた検出アルゴリズムにより検出が行われる。
【００８０】
微小石灰化クラスタの検出アルゴリズムでは、第１実施形態で説明した３重リングフィルタ等により候補領域を特定した後、特徴量を用いた検出処理（処理２）を行う検出過程を経る。各処理内容については第１実施形態で説明した通りであるので、ここでは説明を省略する。
【００８１】
そして、全ての検出過程を終了すると、微小石灰化クラスタの検出結果として、検出された異常陰影候補の位置情報、面積等の情報が記憶部２５に保存される（ステップＥ４）。次いで、制御部２１により記憶部２５から検出結果の情報が読み出されるとともに、微小石灰化クラスタの検出過程に対して予め設定されている異常陰影候補の抽出条件が読み出される。なお、微小石灰化クラスタの検出過程について設定されている抽出条件は、全ての異常陰影候補を抽出することである。すなわち、読み出された検出結果の情報に基づいて、検出された全ての異常陰影候補が表示対象の候補として抽出され、検出された全ての異常陰影候補についてその検出情報が表示部２３上に表示される（ステップＥ５）。
【００８２】
次に、腫瘤が選択された場合について説明する。
ステップＥ６では、異常陰影候補検出部２７において腫瘤に応じた検出処理プログラムが記憶部２５から読み出され、当該プログラムとの協働により異常陰影候補の検出処理が行われる。すなわち、腫瘤の検出に応じた検出アルゴリズムにより検出が行われる。
腫瘤の検出アルゴリズムでは、フィルタバンクを用いた処理１を行う検出過程を経る。
【００８３】
そして、全ての検出過程を終了すると、腫瘤の検出結果として、検出された異常陰影候補の位置情報、面積等の情報が記憶部２５に保存される（ステップＥ７）。次いで、制御部２１により記憶部２５から検出結果の情報が読み出される。そして、検出された異常陰影候補のうち、その検出情報を表示すべき候補が抽出される（ステップＥ８）。
【００８４】
微小石灰化クラスタの場合、その陰影は逆円錐構造の濃度変化を持った微小な白い点が集まって（クラスタ化して）画像上に現れる。そのため、検出された異常陰影候補の検出情報としては、微小な陰影がクラスタ化している部分を囲むフレームのマーカ画像を出力することとなる。つまり、偽陽性候補数が多い場合はその集合領域が大きくなるが、陰影は微小な点であるためその変化量は小さい。よって、異常陰影候補の検出数が多くても、フレームの大きさが変化するだけであり、それほど視認性に影響はない。
【００８５】
しかし、腫瘤の場合、マーカは矢印等、検出された腫瘤の１候補づつを指し示すようなものが用いられるため、検出候補数が多いとその分表示領域を占めることとなるため、読影の妨げとなる場合がある。よって、腫瘤の検出結果については、検出された異常陰影候補のうち、表示対象とする候補を抽出することにより、その検出情報の表示数を減少させることとする。
【００８６】
腫瘤の場合、明らかに真陽性である陰影は目視により医師が確認することが可能である。よって、第１実施形態と同様に、真陽性又は偽陽性とも明確な判別がしづらく、医師の最終判断に委ねた方がよい異常陰影候補、或いはその検出位置や画像上の視認性が低く、医師が発見しづらい異常陰影候補を制御部２１において抽出することとする。なお、抽出方法については、第１実施形態と同一であるので、ここでは説明を省略する。
【００８７】
以上のようにして、表示すべき異常陰影候補が抽出されると、当該抽出された異常陰影候補の検出情報が表示部２３上に表示される（ステップＥ８）。
【００８８】
図１０に、異常陰影候補の検出情報の表示例を示す。
図１０（ａ）は微小石灰化クラスタ候補の表示例を示す図であり、図１０（ｂ）は腫瘤候補の表示例を示す図である。微小石灰化クラスタの場合、点状に現れる微小な石灰化部分の陰影の集合を医師が特定できるように、図１０（ａ）に示すように、その候補が寄り集まる部分を囲むフレームのマーカが使用される。なお、図１０（ａ）では、便宜上微小石灰化クラスタの陰影を黒い点で示しているが、実際には低濃度な白っぽい陰影として現れる。一方、腫瘤は３ｍｍ〜５ｍｍと有る程度の大きさをもって画像上に現れるので、腫瘤候補のマーカとしては、図１０（ｂ）に示すように検出位置を指し示す矢印が用いられている。
【００８９】
微小石灰化クラスタは、偽陽性候補が多く検出された場合でも真陽性の候補の周辺で検出される可能性が高く、偽陽性候補が多くてもその石灰化部分が集合している領域がわずかに変動するのみであると考えられる。よって、検出された全ての候補について検出情報を表示したとしても、図１０（ａ）に示すようにマーカであるフレームのサイズが変化するのみであり、フレーム数が極端に増えるわけではなく、医用画像を読影するにあたって検出情報の表示が妨げになることは少ない。
【００９０】
一方、腫瘤は、１候補について１マーカが表示されるため、偽陽性候補の数が多いと、医用画像上を占める検出情報のマーカ数が多くなり、見づらい表示形態となってしまう。よって、形状や濃度等からは真陽性か偽陽性か判断が困難な陰影に絞ってマーカを表示し、その表示数を減少させることにより、読影しやすい表示形態とすることができる。
【００９１】
なお、図１０（ａ）に示す微小石灰化クラスタの表示と、図１０（ｂ）に示す腫瘤の表示とを医師の操作に応じて切り替え表示可能としてもよいし、同一乳房画像上に表示することしてもよい。
【００９２】
以上のように、本実施形態によれば、検出目的に応じて異なる検出過程を経る場合に、その検出過程に応じて検出された異常陰影候補からその検出情報を表示すべき候補を抽出し、抽出された異常陰影候補についてのみ、検出情報を表示する。これにより、検出過程に応じて検出情報の表示数を制限して、読影しやすい表示形態とすることができ、医師の読影効率を向上させることができる。
【００９３】
また、検出情報としてフレームで検出された候補の集合部分を囲む微小石灰化クラスタの場合は検出された候補の全てを抽出し、矢印で個々の候補の検出位置を指し示す腫瘤の場合は、偽陽性か真陽性かの判断が困難な候補又は視認性が低く、医師が発見しづらい候補を抽出する。このように検出過程により生じる検出率の違いに応じて、その異常陰影候補の検出情報の表示数を可変させるので、表示数が減少するように調整することができる。
【００９４】
なお、第１及び第２実施形態では、画像処理装置２において、異常陰影候補の検出、表示対象候補の抽出等、本発明に係る機能を実現することとしたが、これに限らずビューア５等、医用画像処理装置１００の各構成装置において上記機能のそれぞれを実現させることとしてもよいし、新たにそれらの機能を有する専用装置を設けることとしてもよい。
【００９５】
また、異常陰影候補の検出情報をビューア５により表示出力する例を説明したが、これに限らず、プリンタ３によりフィルム出力する場合も本発明を適用可能である。
【図面の簡単な説明】
【００９６】
【図１】本実施形態における医用画像処理システムのシステム構成を示す図である。
【図２】図１の画像処理装置の内部構成を示す図である。
【図３】第１の実施形態における異常陰影候補の検出処理（一の処理パターン）の流れを説明するフローチャートである。
【図４】フィルタバンクによる検出アルゴリズムを説明する図である。
【図５】フィルタバンクを用いた異常陰影候補の検出処理の流れを説明するフローチャートを説明する図である。
【図６】特徴量を用いた異常陰影候補の検出処理の流れを説明するフローチャートを説明する図である。
【図７】画像処理装置により実行される異常陰影候補の検出処理（他の処理パターン）の流れを説明するフローチャートである。
【図８】各処理パターンによる異常陰影候補の検出情報の表示例を示す図である。
【図９】第２の実施形態における異常陰影候補の検出処理の流れを説明するフローチャートである。
【図１０】検出目的の病変種に応じて行われた検出処理による異常陰影候補の検出情報の表示例を示す図である。
【符号の説明】
【００９７】
１００    医用画像処理システム
１    画像生成装置
２    画像処理装置
２１    制御部
２２    操作部
２３    表示部
２５    記憶部
２７    異常陰影候補検出部
３    プリンタ
４    画像サーバ
５    ビューア
【図１】
2006346094.tif 000010
【図２】
2006346094.tif 000011
【図３】
2006346094.tif 000012
【図４】
2006346094.tif 000013
【図５】
2006346094.tif 000014
【図６】
2006346094.tif 000015
【図７】
2006346094.tif 000016
【図８】
2006346094.tif 000017
【図９】
2006346094.tif 000018
【図１０】
2006346094.tif 000019