(57)【要約】【課題】撮影時期の異なる２つのディジタル胸部画像に経時的差分処理を適用して作成された経時的差分画像を入力画像として画像処理を行って、結節状陰影の候補を検出すると共に、その候補と鑑別診断結果を経時的差分画像上に表示させることで診断結果の精度を向上するようにした診断支援方法とその装置を提供する。【解決手段】撮影時期の異なる２つのディジタル胸部画像１に経時的差分処理を適用して作成された経時的差分画像２を入力画像として画像処理を行って、結節状陰影３の候補を検出するようにした。【選択図】　　　図２
2004048587.tif 000002
【特許請求の範囲】
【請求項１】
撮影時期の異なる２つのディジタル胸部画像に経時的差分処理を適用して作成された経時的差分画像を入力画像として画像処理を行って、結節状陰影の候補を検出することを特徴とするコンピュータ支援診断方法。
【請求項２】
前記画像処理は前記経時的差分画像上の肺野と縦隔を合わせた胸部領域のみに限定して行うことを特徴とする請求項１に記載のコンピュータ支援診断方法。
【請求項３】
前記経時的差分画像に背景傾向補正処理をした後で、請求項１又は２に記載のコンピュータ支援診断方法を適用することを特徴とするコンピュータ支援診断方法。
【請求項４】
前記経時的差分画像を使用して導出した特徴量を使用することを特徴とする請求項１乃至３のいずれか１つに記載のコンピュータ支援診断方法。
【請求項５】
前記ディジタル胸部画像及び前記経時的差分画像の両方を使用して導出した特徴量を使用することを特徴とする請求項１乃至３のいずれか１つに記載のコンピュータ支援診断方法。
【請求項６】
前記特徴量として円形度、面積を含むデータを使用することを特徴とする請求項４又は５に記載のコンピュータ支援診断方法を適用することを特徴とするコンピュータ支援診断方法。
【請求項７】
真の結節状陰影かどうかの判別のために、ルールベース処理を適用することを特徴とする請求項１乃至６のいずれか１つに記載のコンピュータ支援診断方法。
【請求項８】
結節状陰影の候補の抽出に、多数の閾値を切換えながら閾値処理を実行する多重閾値処理を適用することを特徴とする請求項１乃至７のいずれか１つに記載のコンピュータ支援診断方法。
【請求項９】
請求項１乃至８のいずれか１つに記載のコンピュータ支援診断方法を用いて導出された結節状陰影の候補、又は手動で指定された結節状陰影の候補に対して、良性腫瘍、原発性肺がん、肺転移がん等の疾病区分を識別する鑑別診断処理を実施することを特徴とするコンピュータ支援診断方法。
【請求項１０】
請求項１乃至９のいずれか１つに記載のコンピュータ診断支援方法を用いて導出された結節状陰影の候補及び鑑別診断結果を、前記経時的差分画像上又は前記ディジタル胸部画像上に強調表示することを特徴とするコンピュータ支援診断方法。
【請求項１１】
請求項１乃至９のいずれか１つに記載のコンピュータ診断支援方法を用いて導出された結節状陰影の候補及び鑑別診断結果を、前記経時的差分画像上又は前記ディジタル胸部画像上に計測データとして表示することを特徴とするコンピュータ支援診断方法。
【請求項１２】
ディジタル胸部画像を保存する画像データ記憶部と、該画像データ記憶部に接続された画像データ蓄積メモリと、該画像データ蓄積メモリに接続された画像表示処理部と、該画像表示処理部に接続され画像表示装置とを有するコンピュータ支援診断装置であって、
前記画像データ蓄積メモリには、前記画像データ記憶部に保存されている撮影時期の異なる２つのディジタル胸部画像を読み出して経時的差分画像を作成する経時的差分画像作成部と、該経時的差分画像作成部で作成した経時的差分画像と前記ディジタル胸部画像とに基づいて結節状陰影の候補を検出する結節状陰影検出部と、外部からの要求に応じて前記結節状陰影の編集処理を実行する結節状陰影編集部とが接続されており、前記結節状陰影の候補と鑑別診断結果とを前記経時的差分画像又は前記ディジタル胸部画像に重ねて前記画像表示装置に表示するようにしたことを特徴とするコンピュータ支援診断装置。
【請求項１３】
前記画像データ蓄積メモリには、前記経時的差分画像、ディジタル胸部画像及び結節状陰影の候補とに基づいて良性腫瘍、原発性肺がん、肺転移がん等の疾病を識別する鑑別診断処理を実行する鑑別診断処理部を有することを特徴とする請求項１２に記載のコンピュータ支援診断装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
この発明は、撮影時期の異なる２つのディジタル胸部画像に経時的差分処理を適用して作成された経時的差分画像を入力画像として画像処理を行って、結節状陰影の候補を検出すると共に、その候補と鑑別診断結果を記憶装置に保存し、ディスプレイに表示する方法及びその装置に関するものであり、特に、胸部正面画像を用いる医用診断に好適なものである。
【０００２】
【従来の技術】
現在、日本国内においては、様々な集団検診が実施されてきているが、そのうちには、胸部Ｘ線検査、胃部Ｘ線造影検査等の画像診断を必要とするものが含まれている。
【０００３】
特に、肺がんは近年急増しており、早期発見することが必須とされていて、集団検診による画像診断の必要性が特に大きい。しかし、従来のように、胸部単純Ｘ線画像を主に用いた画像診断では、肺がんの見落としの多さが指摘され、その改善が問題となっている。
【０００４】
そこで、集団検診における有病正診率の改善を目的として、ＣＡＤシステム（コンピュータ支援診断システム）が開発されてきているが、その一手法として、同一患者の撮影日時の異なる２つのディジタル胸部正面画像に対して、一方の画像から他方の画像を差し引くことによって、２つの画像に共通の正常構造を取り除き、新しい病巣、又は病巣の変化分だけを強調した経時的差分画像を作ることができる、いわゆる経時的差分法アルゴリズムが開発されてきた。
【０００５】
その１つとして、時間的に連続した２つのディジタル胸部画像を入力として胸部画像間の変化を検出する技術に関しては、特開平７−３７０７４号公報に開示されている。
【０００６】
また、ディジタル化された原胸部画像の信号強調画像と信号抑制画像との差分をとった差分画像を生成し、この差分画像と原胸部画像とを使用して小結節の候補を特定する技術が特開平９−１８５７１４号公報に開示されている。
【０００７】
【発明が解決しようとする課題】
しかしながら、以上のような従来技術にあっては、オリジナルのディジタル胸部画像をそのまま入力として使用するが、オリジナルのディジタル胸部画像上には人体の正常構造が重なり合って写っており、臓器の陰影と結節状陰影との区別が付きにくく、結節状陰影の候補を検出する処理を実行した場合に、誤った結節状陰影を検出してしまう可能性が大きいといった問題点があった。
【０００８】
また、オリジナルのディジタル胸部画像上には経時的変化を表す情報が存在せず、通常、時間の経過で増大又は縮小する結節状陰影の情報を活用することができなかった。
【０００９】
さらに、オリジナルのディジタル胸部画像上に結節状陰影の候補を表示するのは、正常構造と結節状陰影の区別が自明ではなく、読影医にとって見やすい情報になっていなかった。
【００１０】
そこで、この発明は、以上のような従来技術の問題点を解消するために考えられたものであり、撮影時期の異なる２つのディジタル胸部画像に経時的差分処理を適用して作成された経時的差分画像を入力画像として画像処理を行って、結節状陰影の候補を検出すると共に、その候補と鑑別診断結果を経時的差分画像上に表示させることで診断結果の精度を向上するようにした診断支援方法とその装置を提供することを課題としている。
【００１１】
【課題を解決するための手段】
かかる課題を解決するために、請求項１に記載の発明は、撮影時期の異なる２つのディジタル胸部画像に経時的差分処理を適用して作成された経時的差分画像を入力画像として画像処理を行って、結節状陰影の候補を検出することを特徴としている。
【００１２】
請求項２に記載の発明は、請求項１の構成に加えて、前記画像処理は前記経時的差分画像上の肺野と縦隔を合わせた胸部領域のみに限定して行うことを特徴としている。
【００１３】
請求項３に記載の発明は、前記経時的差分画像に背景傾向補正処理をした後で、請求項１又は２に記載のコンピュータ支援診断方法を適用することを特徴としている。
【００１４】
請求項４に記載の発明は、請求項１乃至３のいずれか１つの構成に加えて、前記経時的差分画像を使用して導出した特徴量を使用することを特徴としている。
【００１５】
ここで、特徴量とは画像の特徴を表す指標であって、形状、位置に基づく特徴量と画素値に基づく特徴量とがある。形状、位置に基づく特徴量には、重心位置、面積、実効直径、円形度、周囲長等があり、画素値に基づく特徴量には、平均値、標準偏差、分散、最大値、最小値、歪度、尖度、コントラスト、エネルギー、エントロピー等がある。
【００１６】
請求項５に記載の発明は、請求項１乃至３のいずれか１つの構成に加えて、前記ディジタル胸部画像及び前記経時的差分画像の両方を使用して導出した特徴量を使用することを特徴としている。
【００１７】
請求項６に記載の発明は、請求項４又は５の構成に加えて、前記特徴量として円形度、面積を含むデータを使用することを特徴としている。
【００１８】
請求項７に記載の発明は、請求項１乃至６のいずれか１つの構成に加えて、真の結節状陰影かどうかの判別のために、ルールベース処理を適用することを特徴としている。
【００１９】
ここで、ルールベース処理とは、人工知能の分野ではエキスパートシステムによく利用されている方法のことである。
【００２０】
例えば、この発明にあっては、結節状陰影のデータ及び結節状陰影と間違いやすいデータを多く収集してきて、それらのデータ群の間に存在するルールを解釈して、ｉｆ−ｔｈｅｎルールの形式で与えられたクラス間の決定境界により判別する処理のことをいう。
【００２１】
請求項８に記載の発明は、請求項１乃至７のいずれか１つの構成に加えて、結節状陰影の候補の抽出に、多数の閾値を切換えながら閾値処理を実行する多重閾値処理を適用することを特徴としている。
【００２２】
請求項９に記載の発明は、請求項１乃至８のいずれか１つに記載のコンピュータ支援診断方法を用いて導出された結節状陰影の候補、又は手動で指定された結節状陰影の候補に対して、良性腫瘍、原発性肺がん、肺転移がん等の疾病区分を識別する鑑別診断処理を実施することを特徴としている。
【００２３】
請求項１０に記載の発明は、請求項１乃至９のいずれか１つに記載のコンピュータ支援診断方法を用いて導出された結節状陰影の候補及び鑑別診断結果を、前記経時的差分画像上又は前記ディジタル胸部画像上に強調表示することを特徴としている。
【００２４】
請求項１１に記載の発明は、請求項１乃至９のいずれか１つに記載のコンピュータ支援診断方法を用いて導出された結節状陰影の候補及び鑑別診断結果を、前記経時的差分画像上又は前記ディジタル胸部画像上に計測データとして表示することを特徴としている。
【００２５】
請求項１２に記載の発明は、ディジタル胸部画像を保存する画像データ記憶部と、該画像データ記憶部に接続された画像データ蓄積メモリと、該画像データ蓄積メモリに接続された画像表示処理部と、該画像表示処理部に接続され画像表示装置とを有するコンピュータ支援診断装置であって、前記画像データ蓄積メモリには、前記画像データ記憶部に保存されている撮影時期の異なる２つのディジタル胸部画像を読み出して経時的差分画像を作成する経時的差分画像作成部と、該経時的差分画像作成部で作成した経時的差分画像と前記ディジタル胸部画像とに基づいて結節状陰影の候補を検出する結節状陰影検出部と、外部からの要求に応じて前記結節状陰影の編集処理を実行する結節状陰影編集部とが接続されており、前記結節状陰影の候補と鑑別診断結果とを前記経時的差分画像又は前記ディジタル胸部画像に重ねて前記画像表示装置に表示するようにしたことを特徴としている。
【００２６】
請求項１３に記載の発明は、前記画像データ蓄積メモリには、前記経時的差分画像、ディジタル胸部画像及び結節状陰影の候補とに基づいて良性腫瘍、原発性肺がん、肺転移がん等の疾病を識別する鑑別診断処理を実行する鑑別診断処理部を有することを特徴としている。
【００２７】
【発明の実施の形態】
この発明の実施の形態に係るコンピュータ支援診断方法とその装置について、図面に従って説明する。
【００２８】
図１は、この発明の実施の形態に係るコンピュータ支援診断装置の構成を示したブロック図である。
【００２９】
この発明の実施の形態に係るコンピュータ支援診断装置Ｍは、ディジタル胸部画像１を入力し、経時的差分画像２と結節状陰影３とを入出力する画像データ記憶部４と、この画像データ記憶部４と接続された画像データ蓄積メモリ５とを有している。画像データ蓄積メモリ５では、ディジタル胸部画像１を入力画像とし、このディジタル胸部画像１を使用して経時的差分画像作成部６で作成された経時的差分画像２に基づいて、結節状陰影検出部７、鑑別診断処理部８、結節状陰影編集部９の処理を行い、その結果として結節状陰影３を出力する。この出力された結節状陰影３は経時的差分画像２又はディジタル胸部画像１と一緒にして、画像表示装置１１に表示できる画像データに変換する画像表示処理部１０を経て画像表示装置１１に表示する。これにより、結節状陰影３の候補及び鑑別診断結果を、経時的差分画像２又はディジタル胸部画像１に重ねて表示することが可能となる。
【００３０】
ディジタル胸部画像１は、この発明において入力画像となる胸部正面のディジタル画像であり、ファイル、データベース等の形式で格納されている。
【００３１】
経時的差分画像２は、撮影時期の異なる２つのディジタル胸部画像１に経時的差分処理を適用して作成された経時的差分画像であって、この発明において入力画像となる共に出力画像の一部となり、ファイル、データベース等の形式で格納されている。
【００３２】
結節状陰影３は、この発明において出力画像となる結節状陰影のデータであり、ファイル、データベース等の形式で格納されている。
【００３３】
画像データ記憶部４は、ディジタル胸部画像１を指定された要求に基づいて読み込んで、画像データ蓄積メモリ５に格納する処理を行う部分である。また、経時的差分画像作成部６で作成された経時的差分画像２、及び結節状陰影検出部７、結節状陰影編集部９、鑑別診断処理部８で作成された結節状陰影３を画像データ蓄積メモリ５から取り出し、指定された要求に基づいて、ファイル、データベース等の何らかの形式で保存する部分である。
【００３４】
画像データ蓄積メモリ５は、画像データ記憶部４で読み込まれたディジタル胸部画像１と、経時的差分画像作成部６で作成された経時的差分画像２、結節状陰影検出部７、結節状陰影編集部９及び鑑別診断処理部８で作成された結節状陰影３を一時的に蓄積しておく部分である。
【００３５】
経時的差分画像作成部６は、画像データ蓄積メモリ５から入力された撮影時期の異なる時系列に並んだ２つのディジタル胸部画像１を基にして、１つの経時的差分画像２を作成して画像データ蓄積メモリ５に出力する部分である。
【００３６】
結節状陰影検出部７は、画像データ蓄積メモリ５から入力された経時的差分画像２、ディジタル胸部画像１を基にして、結節状陰影３の候補を検出して画像データ蓄積メモリ５に出力する部分である。
【００３７】
鑑別診断処理部８は、画像データ蓄積メモリ５から入力された結節状陰影３、経時的差分画像２及びディジタル胸部画像１を基にして、鑑別診断処理を実行して、鑑別診断結果を結節状陰影３の部分データとして画像データ蓄積メモリ５に出力する部分である。
【００３８】
結節状陰影編集部９は、画像表示処理部１０を経て画像表示装置１１に表示されたディジタル胸部画像１又は経時的差分画像２の上で、指定された要求に基づいて、結節状陰影３を編集、入力、削除する処理を行う部分。
【００３９】
画像表示処理部１０は、画像データ蓄積メモリ５に格納されているディジタル胸部画像１又は経時的差分画像２、及び結節状陰影３を、指定された要求に基づいて画像表示装置１１に表示する処理を行う部分である。
【００４０】
画像表示装置１１は、ＣＲＴ、プラズマディスプレイ、液晶ディスプレイ等の表示装置であり、医用診断に使用するため、走査線１０００本系以上の高精細表示装置であることが望ましい。
【００４１】
以下、この発明の実施形態に係るコンピュータ診断支援装置が実現する処理について、詳細に説明する。
【００４２】
まず、結節状陰影検出の詳細について説明する。
【００４３】
図２は、結節状陰影検出の手順を示したフローチャートである。
【００４４】
まず、経時的差分画像２の作成において、画像変形をしていない方のディジタル胸部画像１ａを用いて、胸部領域１２を抽出し、輪郭情報として出力する（ステップＳ１）。
【００４５】
胸部領域１２の抽出の具体的処理としては、ディジタル胸部画像１ａ上から、胸郭の外部輪郭境界（左右と上側の肋骨で象られる輪郭）を検出し、輪郭がない下部の部分を直線で結んで、一つの閉じた輪郭にする。胸部領域１２として、肺野と縦隔を合せた領域を抽出することができれば十分であり、前記処理はその目的を達している。
【００４６】
結節状陰影３は、胸部領域１２の内部にのみ存在するとみなすことができる。以降の処理においては、ステップＳ１で導出された胸部領域１２を利用することにより、計算効率を向上させることができる。
【００４７】
次に、経時的差分画像２に対して、背景傾向補正処理（ステップＳ２）を適用して、新しい経時的差分画像としての補正済み経時的差分画像１３を作成する。
【００４８】
ステップＳ２の背景傾向補正処理は、具体的には次の３段階を経て実行する。
【００４９】
まず、経時的差分画像２は、ステップＳ１で作成した胸部領域１２の外部を全て固定値にする。指定する固定値は、経時的差分画像２の胸部領域１２の内部の画素平均値とする。この処理により、胸部領域１２の外部の画素分布の影響が除去される。
【００５０】
次に、経時的差分画像２が構成する曲面の近似曲面を、２次元多項式近似を用いて導出する。
【００５１】
最後に、経時的差分画像２から近似曲面を差し引き、適当なオフセットを全体に付加して画素値を正の値になるようにして、画像を出力する。この３番目の処理を実際に適用するのは胸部領域１２の内部のみとし、胸部領域１２の外部については画素値０を与える。
【００５２】
ステップＳ２で導出された補正済み経時的差分画像１３に対して、多重閾値処理を適用することにより、初期結節状陰影候補１４を拾い上げる（ステップＳ３）。初期結節状陰影候補１４は、補正済み経時的差分画像１３上の連結領域として与えられる。
【００５３】
次に、補正済み経時的差分画像１３及びディジタル胸部画像１ａを用いて、ステップＳ３で導出された初期結節状陰影候補１４のそれぞれに関する特徴量１５を計算する（ステップＳ４）。
【００５４】
ここで、計算される特徴量１５の具体的な例としては、次のようなものがあげられる。
【００５５】
形状、位置に基づく特徴量１５としては、連結領域の重心座標（画像座標）である重心位置、連結領域の画素数である面積、面積が一致する円の直径である実効直径、重心位置において実効直径の円を描いたときに円内に含まれる陰影候補の面積の割合である円形度、連結領域の周囲長である周囲長等がある。画素値に基づく特徴量１５としては、平均値、標準偏差、分散、最大値、最小値、歪度、尖度等の一般統計量やコントラスト、エネルギー、エントロピー等がある。
【００５６】
直接の結節状陰影検出対象である経時的差分画像２のみならず、オリジナルのディジタル胸部画像１ａから求めた特徴量１５を使用することにより、人体の正常構造など、経時的差分画像２上では現れない特徴を生かした処理を行うことが可能となる。
【００５７】
次に、ステップＳ４で導出された特徴量１５を使用して、ルールベース処理により判定を行い、偽陽性候補を決定する（ステップＳ５）。偽陽性候補を初期結節状陰影候補１４から除去することにより、最終的な結節状陰影候補１６が得られる。
【００５８】
なお、経時的差分画像２は、原理的に、経時的変化分以外は、肋骨等の正常構造を含めて消えた画像になっているはずであるが、撮影体位のずれ、呼吸位のずれ等を原因とする偽画像成分が残っていることが多い。ステップＳ２の背景傾向補正処理によって、これらの偽画像成分のうち直流成分を取り除くことができる。その結果、多重閾値処理により結節状陰影候補１６を拾い上げるときに、画像上の位置の違いによる影響が小さくなり、安定して処理を実行できる。
【００５９】
以下、結節状陰影候補拾い上げの詳細について説明する。
【００６０】
図３は、結節状陰影候補拾い上げの手順を示したフローチャートである。　結節状陰影候補１６を拾い上げるために、多数の閾値を切換えながら閾値処理を実行する多重閾値処理を適用する。
【００６１】
まず、補正済み経時的差分画像１３の画素値に関するヒストグラム１７を作成する（ステップＳ６）。補正済み経時的差分画像１３上では、画素値０の領域は胸部領域１２の外部を表すので、ヒストグラム作成の際には使用しない。
【００６２】
ステップＳ６で導出されたヒストグラム１７を使用して、ヒストグラム面積に対する％比率で定義して多重閾値１８を導出する（ステップＳ７）。
【００６３】
多重閾値１８を導出するための具体的な処理としては、画素値の小さい方から順にヒストグラム面積を累積していって、予め％値で設定された値に一致する画素値を閾値として導出する。通常、％値は複数設定されているので、得られる閾値も複数であり、この発明ではこれを多重閾値と呼んでいる。
【００６４】
ステップＳ７で導出された多重閾値１８をそれぞれ、補正済み経時的差分画像１３に適用して、２値画像１９を作成する（ステップＳ８）。
【００６５】
ステップＳ８で導出された２値画像１９のそれぞれに対して、閾値以下を表す連結領域をラベリングして島領域２０を全て導出する（ステップＳ９）。
【００６６】
ステップＳ９で導出された島領域２０には、互いに包含された状態にある重複した島領域２０が存在する。この重複分を除去した（ステップＳ１０）後の島領域を全て、初期結節状陰影候補１４とする。
【００６７】
ステップＳ１０の島領域除去の具体的処理としては、異なる％値の島領域を組合せて、互いに包含された状態にあるかどうかを検知して、％値の大きい方の島領域（すなわち、面積の大きい方の島領域）を陰影候補から除去する。これにより、％値の小さい、すなわち画素値のできるだけ小さい部分で拾い上げられた結節状陰影３の候補を出力することができる。
【００６８】
次に、閾値導出の考え方について説明する。
【００６９】
補正済み経時的差分画像１３上では、結節状陰影候補１６は、周囲に比べて相対的に暗い（画素値の小さい）領域として表示されている。したがって、２値画像１９を作成したときに島領域２０として現れる結節状陰影候補１６は、閾値よりも小さい領域の方である。結節状陰影候補１６を洩れなく拾い上げるためには、画素値の小さい方から順に閾値を多段階に設定していって、それぞれ島領域を抽出していくものとする。
【００７０】
具体的には、画素ヒストグラムの最小画素値を開始点として境界画素値を移動していって、全体面積に対する累積面積の比率が指定した％値になる画素値を、閾値として使用する。指定する境界画素値は、３％、６％、９％等と、％値で予め設定しておく。
【００７１】
図４は、多重閾値導出の具体例を示したグラフである。
【００７２】
多重閾値は、％比率で設定された境界画素値を基に、画素値に関するヒストグラム１７から導出される。
【００７３】
以下、鑑別診断処理の手順について詳細に説明する。
【００７４】
図５は、鑑別診断処理の手順を示したフローチャートである。
【００７５】
鑑別診断処理は図５（ａ）の学習処理と、図５（ｂ）の識別処理との２つの処理で構成される。
【００７６】
学習処理の最初には、補正済み経時的差分画像１３及びディジタル胸部画像１を用いて、入力された結節状陰影３に対する特徴量１５ｇを計算する（ステップＳ１１）。
【００７７】
ここで、直接の結節状陰影検出対象である経時的差分画像２のみならず、オリジナルのディジタル胸部画像１ａを用いることにより、人体の正常構造等、経時的差分画像２上では現れない特徴を生かした処理を行うことができるといった利点が考えられる。
【００７８】
次に、ステップＳ１１で作成された特徴量１５ｇを入力として、多クラス識別器を学習して多クラス識別器２１を出力する（ステップＳ１２）して、学習処理は終了となる。出力された多クラス識別器２１は、次の識別処理において使用される。
【００７９】
識別処理では、まず、学習処理における特徴量計算（ステップＳ１１）と同様に特徴量１５ｓを計算する（ステップＳ１３）。
【００８０】
次に、ステップＳ１３で作成された特徴量１５ｓを入力及びステップＳ１２の多クラス識別器の学習により出力された多クラス識別器２１を入力として多クラス識別器を実行し（ステップＳ１４）、鑑別診断結果２２を出力する。鑑別診断処理の実行装置には、最終的には多クラス識別器２１を組込んで、識別処理を動作させるようにするとよい。
【００８１】
学習処理は、その前段階として実行させておく必要がある。多クラス識別器の構造をどんなものにするにせよ、識別対象の結節状陰影３の特徴をよく表しているデータを相当数用意しておき、ジャックナイフ法等の適切な手順に従って学習させる。また、それぞれの識別器の特徴を活かして、適用方法を工夫する。
【００８２】
多クラス識別器２１としては、例えば、線形判別の多段実行による多クラス分類や人工ニューラルネットワーク等のような識別器を使用すればよい。
【００８３】
鑑別診断の識別対象としては、良性腫瘍、原発性肺がん、肺転移がん等の疾病を識別し、鑑別診断結果２２として出力するものとする。
【００８４】
以下、画像表示処理の詳細について説明する。
【００８５】
画像表示処理においては、基本的には、ディジタル胸部画像１、経時的差分画像２を画像表示装置１１上に表示して読影を行う。この発明の実施の形態に係るコンピュータ支援診断装置Ｍに基づいて作成された結節状陰影３の候補及び鑑別診断結果２２は、経時的差分画像２上又はディジタル胸部画像１上に重畳して表示して、読影を支援する。
【００８６】
図６に、経時的差分画像をそのまま表示した例の概略図を示した。図７に、
結節状陰影の強調表示した例の概略図を示した。
【００８７】
この発明の実施の形態に係るコンピュータ支援診断装置Ｍによれば、結節状陰影検出処理により作成された結節状陰影３の候補を使用して、図６に示したような、経時的差分画像２上に淡く写っている結節状陰影３を、図７に示したように強調表示することができる。強調表示することにより、読影者は、通常の読影手順の延長線上で、自動検出されたかどうかは意識せずに、結節状陰影３を自然と認識できるようになる。
【００８８】
また、結節状陰影検出処理により作成された結節状陰影３の候補を使用して、経時的差分画像２上に淡く写っている結節状陰影３の場所を、計測データ表示の一環として、指し示すようにすることができる。この計測データ表示により、読影者は明確に結節状陰影３の候補を認識することができる。
【００８９】
鑑別診断結果については、結節状陰影３の候補を色又はマークによって区分けして表示する。
【００９０】
計測データ表示においては、例えば、良性腫瘍は緑色、原発性肺がんは黄色、肺転移がんは赤色というように交差線分を表示することで、区分の識別を明示することができる。また、＋印、×印、□印というように、計測データ表示のマークを変えることもできる。
【００９１】
図８は、計測データ表示の交差線分表示の例を示した概略図であり、結節状陰影候補の中心に交差線分を表示している。図９は、拡大表示の例を示した概略図であり、結節状陰影候補を含む領域を拡大表示している。
【００９２】
【実施例】
次に、この発明の実施の形態に係るコンピュータ支援診断装置Ｍを適用した実施例について示す。
【００９３】
まず、結節状陰影検出処理の実施例である症例の一覧を、表１に示す。
【００９４】
【表１】
2004048587.tif 000003
【００９５】
症例データベースとしては、産業医科大学病院の放射線部で選択された９０例を使用した。これらの症例は全て異常例であり、不明瞭なものが８４％を占めていて、比較的難易度が高いといえる（難易度の判定は１０人の読影者によって行われた際の平均値によった）。
【００９６】
各症例は、撮像日時の異なる２つのディジタル胸部正面画像として用意され、経時的差分画像作成処理を実行して経時的差分画像を作成した後に、結節状陰影検出処理を実行した。
【００９７】
結節状陰影は時間的に新しい画像の方に存在し、主観的難易度はこの画像を読影して決定されたものである。
【００９８】
表２は、実施例（全９０症例）に対する検出性能を示している。表３は、陰影種類別の検出感度の内訳を示している。
【００９９】
【表２】
2004048587.tif 000004
【０１００】
【表３】
2004048587.tif 000005
【０１０１】
表２及び３に示した実施例における検出性能は、予め用意された適用症例に対して設定調整されたものである。比較的難易度の高い症例を使用しているため、同一撮像条件のディジタル胸部画像であれば、安定した性能を発揮することがわかる。
【０１０２】
なお、この発明は、この実施例をもとに、さらに多くの改変を適用してから実施することが可能であることは言うまでもない。
【０１０３】
表４は、結節状陰影処理の設定パラメータの一覧表である。
【０１０４】
この実施例では、結節状陰影検出処理全体を２回ずつ実行するような仕組みにしている。それぞれの処理実行パスをパス１、パス２としている。最終結果としては、両方の実行結果を統合して出力するものとした。
【０１０５】
これにより、処理を２回実行することにより、異なる特性を持つ結節状陰影が検出できることがわかる。
【０１０６】
【表４】
2004048587.tif 000006
【０１０７】
実施例（症例）のディジタル胸部画像は１０ビットの画素値（０から１０２３までの値）を持っており、偽陽性候補判定ルールにおける標準偏差の値は、その事実に基づいて設定された。
【０１０８】
結節状陰影拾い上げ時の実行結果として、図１０は結節状陰影（陽性）の分布の例を示した円形度と面積との関係のグラフであり、図１１は偽陽性候補の分布の例を示した円形度と面積との関係のグラフである。
【０１０９】
ルールベース処理の設定について、その偽陽性候補の判定のために実行するルールベース処理の閾値の例を示したグラフを図１２と図１３に示した。図１２は面積と円形度との関係を示し、図１３に円形度と標準偏差との関係を示している。
【０１１０】
ルールベース処理の検出性能（処理別の内訳）を、拾い上げられた結節状陰影（全８１症例）に対して計測した。その結果は表５のとおりであり、検出感度は処理実行の２回目では１００％を達成した。
【０１１１】
【表５】
2004048587.tif 000007
【０１１２】
【発明の効果】
以上説明したように、請求項１に記載の発明によれば、撮影時期の異なる２つのディジタル胸部画像に経時的差分処理を適用して作成された経時的差分画像を入力画像として画像処理を行って、結節状陰影の候補を検出するので、時間経過で増大又は縮小する結節状陰影の情報を活用することが可能となり、かつ結節状陰影の候補が強調して表示されることとなるため、読影者の読影にかかる負担を軽減でき、読影時間も短縮できる。
【０１１３】
請求項２に記載の発明によれば、画像処理は経時的差分画像上の肺野と縦隔を合わせた胸部領域のみに限定して行うので、無駄な画像処理を除外することができ、請求項１の効果に加えて、画像処理にかかる時間をより短縮することが可能となる。
【０１１４】
請求項３に記載の発明によれば、経時的差分画像に背景傾向補正処理をした後で、請求項１又は２に記載のコンピュータ支援診断方法を適用するので、偽画像成分のうち直流成分を取り除くことができ、多重閾値処理により結節状陰影候補を拾い上げるときに、画像上の位置の違いによる影響が小さくなり、安定した画像処理が実行されるため、請求項１又は２の効果に加えて、より検出感度を向上させることが可能となる。
【０１１５】
請求項４に記載の発明によれば、経時的差分画像を使用して導出した特徴量を使用するので、請求項１乃至３のいずれか１つに記載の効果に加えて、より検出感度を向上させることが可能となる。
【０１１６】
請求項５に記載の発明によれば、ディジタル胸部画像及び前記経時的差分画像の両方を使用して導出した特徴量を使用するので、オリジナルのディジタル胸部画像に含まれている人体の正常構造の特徴を生かした画像処理が実行されるため、請求項１乃至３のいずれか１つに記載の効果に加えて、より検出感度を向上させることが可能となる。
【０１１７】
請求項６に記載の発明によれば、特徴量として円形度、面積を含むデータを使用するので、請求項４又は５に記載の効果に加えて、放射状をした結節状陰影がより効率的に画像処理できる。
【０１１８】
請求項７に記載の発明によれば、真の結節状陰影かどうかの判別のためにルールベース処理を適用するので、過去の経験則が精度よく反映された結果を提供することができるため、請求項１乃至６のいずれか１つに記載の効果に加えて、比較的容易により検出感度を向上させることが可能となる。
【０１１９】
請求項８に記載の発明によれば、結節状陰影の候補の抽出に、多数の閾値を切換えながら閾値処理を実行する多重閾値処理を適用するので、画素値の小さいものから大きなものまで的確に捉えることができるため、請求項１乃至７のいずれか１つに記載の効果に加えて、結節状陰影の小さなものから大きなものまで的確に捉えることができるため、より検出感度を向上させることが可能となる。
【０１２０】
請求項９に記載の発明によれば、請求項１乃至８のいずれか１つに記載のコンピュータ支援診断方法を用いて導出された結節状陰影の候補、又は手動で指定された結節状陰影の候補に対して、良性腫瘍、原発性肺がん、肺転移がん等の疾病区分を識別する鑑別診断処理を実施するので、具体的な疾病区分の情報が得られることとなり、請求項１乃至８のいずれか１つの効果に加えて、病変判別の精度が向上する。
【０１２１】
請求項１０に記載の発明によれば、請求項１乃至９のいずれか１つに記載のコンピュータ支援診断方法を用いて導出された結節状陰影の候補及び鑑別診断結果を、経時的差分画像上又はディジタル胸部画像上に強調表示するので、請求項１乃至９のいずれか１つの効果に加えて、より一層病変判別の精度が向上する。
【０１２２】
請求項１１に記載の発明によれば、請求項１乃至９のいずれか１つに記載のコンピュータ支援診断方法を用いて導出された結節状陰影の候補及び鑑別診断結果を、経時的差分画像上又はディジタル胸部画像上に計測データとして表示するので、疾病区分の情報を参考にする際の目安が与えられるため、請求項１乃至９のいずれか１つの効果に加えて、読影者の負担が一層軽減される。
【０１２３】
請求項１２に記載の発明によれば、画像データ蓄積メモリには、画像データ記憶部に保存されている撮影時期の異なる２つのディジタル胸部画像を読み出して経時的差分画像を作成する経時的差分画像作成部と、該経時的差分画像作成部で作成した経時的差分画像とディジタル胸部画像とに基づいて結節状陰影の候補を検出する結節状陰影検出部と、外部からの要求に応じて結節状陰影の編集処理を実行する結節状陰影編集部とが接続されており、結節状陰影の候補と鑑別診断結果とを経時的差分画像又はディジタル胸部画像に重ねて画像表示装置に表示するようにしたので、時間経過で増大又は縮小する結節状陰影の情報を活用することが可能となり、かつ結節状陰影の候補が強調して表示されることとなるため、読影者の読影にかかる負担を軽減でき、読影時間も短縮できる。
【０１２４】
請求項１３に記載の発明によれば、画像データ蓄積メモリには、経時的差分画像、ディジタル胸部画像及び結節状陰影の候補とに基づいて良性腫瘍、原発性肺がん、肺転移がん等の疾病を識別する鑑別診断処理を実行する鑑別診断処理部を有するので、結節状陰影の候補、又は手動で指定された結節状陰影の候補に対して、良性腫瘍、原発性肺がん、肺転移がん等の疾病区分を識別する鑑別診断処理を実施できるため、具体的な疾病区分の情報が得られることとなり、請求項１２の効果に加えて、病変判別の精度が一層向上する。
【図面の簡単な説明】
【図１】この発明の実施の形態に係るコンピュータ支援診断装置の構成を示したブロック図である。
【図２】同コンピュータ支援診断装置を使用した結節状陰影検出の手順を示したフローチャートである。
【図３】同コンピュータ支援診断装置を使用した結節状陰影候補拾い上げの手順を示したフローチャートである。
【図４】同コンピュータ支援診断装置を使用した多重閾値導出の具体例を示したグラフである。
【図５】同コンピュータ支援診断装置を使用した鑑別診断処理の手順を示したフローチャートである。
【図６】同コンピュータ支援診断装置の経時的差分画像をそのまま表示した例の概略図である。
【図７】同コンピュータ支援診断装置の結節状陰影の強調表示した例の概略図である。
【図８】同コンピュータ支援診断装置の計測データ表示の交差線分表示の例を示した概略図である。
【図９】同コンピュータ支援診断装置の計測データ表示の拡大表示の例を示した概略図である。
【図１０】同コンピュータ支援診断装置の結節状陰影（陽性）の分布の例を示した円形度と面積との関係のグラフである。
【図１１】同コンピュータ支援診断装置の偽陽性候補の分布の例を示した円形度と面積との関係のグラフである。
【図１２】同コンピュータ支援診断装置のルールベース処理の閾値の例を示した面積と円形度との関係のグラフである。
【図１３】同コンピュータ支援診断装置のルールベース処理の閾値の例を示した円形度と標準偏差との関係のグラフである。
【符号の説明】
Ｍ　コンピュータ支援診断装置
１　ディジタル胸部画像
２　経時的差分画像
３　結節状陰影
４　画像データ記憶部
５　画像データ蓄積メモリ
６　経時的差分画像作成部
７　結節状陰影検出部
８　鑑別診断処理部
９　結節状陰影編集部
１０　画像表示処理部
１１　画像表示装置
Ｓ１　胸部領域の抽出
Ｓ２　背景傾向補正
Ｓ３　結節状陰影候補の拾い上げ
Ｓ４　特徴量計算
Ｓ５　偽陽性候補の判定
【図１】
2004048587.tif 000008
【図２】
2004048587.tif 000009
【図３】
2004048587.tif 000010
【図４】
2004048587.tif 000011
【図５】
2004048587.tif 000012
【図６】
2004048587.tif 000013
【図７】
2004048587.tif 000014
【図８】
2004048587.tif 000015
【図９】
2004048587.tif 000016
【図１０】
2004048587.tif 000017
【図１１】
2004048587.tif 000018
【図１２】
2004048587.tif 000019
【図１３】
2004048587.tif 000020