(57)【要約】【課題】現在撮影した検査対象画像と、過去に正常と診断された他人の画像データベースから、確実に異常部位を検出する。【解決手段】正常な多数の画像を収集したデータベースと、検査対象画像を、データベースに記憶されている画像から生成した疑似正常画像と対比して差分画像を作成する疑似正常差分画像生成部と、疑似正常差分画像中から異常陰影を検出する異常陰影検出部と、検出した異常陰影を表示する表示部とを備える。疑似正常画像を作成するために、検査対象画像の一部を切り出した検査画像パッチと、データベース画像の一部を切り出した類似画像パッチとを局所的に比較し、類似したものを選択することを、検査対象画像の全体を走査することによって、検査対象画像の正常であった時に類似の疑似正常画像を作り上げる。【選択図】        図１
2014221163.tif 000002
【特許請求の範囲】
【請求項１】
検査対象画像をそれに類似する画像と対比して差分画像を作成することにより病巣候補の検出を行う画像診断支援装置において、
病巣が無いことが分かっている正常な多数の画像を収集して記憶したデータベースと、
前記検査対象画像を、前記データベースに記憶されている画像から生成した疑似正常画像と対比して疑似正常差分画像を作成する疑似正常差分画像生成部と、
前記疑似正常差分画像中から異常陰影を検出する異常陰影検出部と、
検出した前記異常陰影を表示する表示部と、を備え、
前記疑似正常差分画像生成部は、前記疑似正常画像を作成するために、前記検査対象画像から画像の一部を局所的に切り出した検査画像パッチと、前記データベースに記憶されている画像の一部を切り出した類似画像パッチとを局所的に比較し、類似したものを選択することを、前記検査対象画像の全体を走査することによって、前記検査対象画像の正常であった時に類似の前記疑似正常画像を全体的に作り上げる、
ことから成る画像診断支援装置。
【請求項２】
さらに、検出した異常陰影を編集して、前記疑似正常差分画像の正の階調値に対して線形変換することにより異常陰影の強調表示を行う異常陰影編集部を備える請求項１に記載の画像診断支援装置。
【請求項３】
検査対象画像をそれに類似する画像と対比して差分画像を作成することにより病巣候補の検出を行う画像診断支援方法において、
病巣が無いことが分かっている正常な多数の画像を収集して記憶したデータベースを備え、
前記検査対象画像を、前記データベースに記憶されている画像から生成した疑似正常画像と対比して疑似正常差分画像を作成し、
前記疑似正常差分画像中から異常陰影を検出し、かつ、検出した異常陰影を表示し、
前記疑似正常画像を作成するために、前記検査対象画像から画像の一部を局所的に切り出した検査画像パッチと、前記データベースに記憶されている画像の一部を切り出した類似画像パッチとを局所的に比較し、類似したものを選択することを、前記検査対象画像の全体を走査することによって、前記検査対象画像の正常であった時に類似の前記疑似正常画像を全体的に作り上げる、
ことから成る画像診断支援方法。
【請求項４】
前記検査画像パッチと前記類似画像パッチとを局所的に比較して、類似度計算により各々のデータベース画像に対応してそれぞれ１つの類似画像パッチを決定し、この決定された類似画像パッチの内、最も類似度が高いものを一つの検査画像パッチに対応する最類似画像パッチとしてデータベース内から決定し、かつ、前記検査対象画像の全体で決定した最類似画像パッチを全て用いて再構築することにより、検査対象画像に対応して前記疑似正常画像を生成する請求項３に記載の画像診断支援方法。
【請求項５】
前記疑似正常差分画像の正の階調値に対して線形変換することにより異常陰影の強調補正を行う請求項３に記載の画像診断支援方法。
【請求項６】
前記強調補正をした疑似正常差分画像から異常陰影候補の拾い上げを行なって異常陰影初期候補を決定し、この異常陰影初期候補の特徴量を計算し、この特徴量を基に偽陽性候補を決定し、この偽陽性候補を前記異常陰影初期候補から除去することにより、前記異常陰影を決定する請求項５に記載の画像診断支援方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、検査対象画像をそれに類似する画像と対比して差分画像を作成することにより病巣候補の検出を行う画像診断支援装置及び支援方法に関する。
【背景技術】
【０００２】
肺癌などの肺の病巣を早期に診断できるようにすることは重要なことである。過去の正常時の撮影データがある場合には、経時差分処理により、現在の異常を検出できる技術が開発されている。診断を受ける患者の過去の正常な胸部写真が存在する場合には、その画像の差分を取れば比較的簡単に診断が可能であるが、正常（健康）時に胸部写真を撮っておく人は少ない。そこで他人の正常な胸部写真を利用して、患者の胸部写真との差分を取ることで診断できないかが試みられている。
【０００３】
特許文献１は、他人の正常な胸部写真を利用して、差分を取ることで診断する技術を開示する。診断対象の患者についてのディジタル胸部画像に対して、この患者とは異なる他人のディジタル胸部画像の中から類似する画像を類似画像として取得する。この両画像の差分を取ることによって、この２つの画像に共通の正常構造を取り除き、病巣だけを強調した画像を取得することができる。
しかし、現在、他人のディジタル胸部画像を利用する良い装置はまだできていない。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開2005-269214
【発明の概要】
【発明が解決しようとする課題】
【０００５】
既存手法では、検査対象の画像に類似する画像を生成する際、他人の胸部Ｘ線画像データベースの中から、年齢、性別等の属性により一次的に絞り込み、その中から、大局的にマッチングを行い、最も類似する1枚の画像を変形させて使用している。1枚の画像の変形のみでは画像の細部に亘って類似画像を構成することが困難であり、結果として、検出精度の劣化を招く。また、変形に伴うアーチファクト（偽所見）が生じやすく、このことも類似画像生成を困難にしている。
【０００６】
それ故に、本発明は、現在撮影した検査対象の１枚の画像（例えば、胸部Ｘ線画像）と、過去に正常と診断された他人の画像データベースから、確実に異常部位を検出することを目的としている。
【課題を解決するための手段】
【０００７】
本発明の画像診断支援装置は、検査対象画像をそれに類似する画像と対比して差分画像を作成することにより病巣候補の検出を行う。この画像診断支援装置は、病巣が無いことが分かっている正常な多数の画像を収集して記憶したデータベースと、検査対象画像を、データベースに記憶されている画像から生成した疑似正常画像と対比して疑似正常差分画像を作成する疑似正常差分画像生成部と、疑似正常差分画像中から異常陰影を検出する異常陰影検出部と、検出した異常陰影を表示する表示部とを備える。疑似正常差分画像生成部は、疑似正常画像を作成するために、検査対象画像から画像の一部を局所的に切り出した検査画像パッチと、データベースに記憶されている画像の一部を切り出した類似画像パッチとを局所的に比較し、類似したものを選択することを、検査対象画像の全体を走査することによって、検査対象画像の正常であった時に類似の疑似正常画像を全体的に作り上げる。
【０００８】
また、本発明の画像診断支援方法は、病巣が無いことが分かっている正常な多数の画像を収集して記憶したデータベースを備え、検査対象画像を、データベースに記憶されている画像から生成した疑似正常画像と対比して疑似正常差分画像を作成し、疑似正常差分画像中から異常陰影を検出し、かつ、検出した異常陰影を表示する。疑似正常画像を作成するために、検査対象画像から画像の一部を局所的に切り出した検査画像パッチと、データベースに記憶されている画像の一部を切り出した類似画像パッチとを局所的に比較し、類似したものを選択することを、検査対象画像の全体を走査することによって、検査対象画像の正常であった時に類似の疑似正常画像を全体的に作り上げる。
【０００９】
さらに、検出した異常陰影を編集して、疑似正常差分画像の正の階調値に対して線形変換することにより異常陰影の強調表示を行う異常陰影編集部を備える。
検査画像パッチと類似画像パッチとを局所的に比較して、類似度計算により各々のデータベース画像に対応してそれぞれ１つの類似画像パッチを決定し、この決定された類似画像パッチの内、最も類似度が高いものを一つの検査画像パッチに対応する最類似画像パッチとしてデータベース内から決定し、かつ、検査対象画像の全体で決定した最類似画像パッチを全て用いて再構築することにより、検査対象画像に対応して疑似正常画像を生成する。
【００１０】
疑似正常差分画像の正の階調値に対して線形変換することにより異常陰影の強調補正を行う。強調補正をした疑似正常差分画像から異常陰影候補の拾い上げを行なって異常陰影初期候補を決定し、この異常陰影初期候補の特徴量を計算し、この特徴量を基に偽陽性候補を決定し、この偽陽性候補を異常陰影初期候補から除去することにより、異常陰影を決定する。
【発明の効果】
【００１１】
本発明によれば、他人の正常画像から、患者の現在の正常画像を擬似的に生成し、現在の撮影データとの差分処理を行うことで、異常部位を確実に検出することができる。
【００１２】
本発明は、異常陰影の存在しない正常な胸部Ｘ線画像のみで構成されたデータベースから局所的にマッチングを行い、局所領域まで検査対象画像に類似した正常画像を生成することにより、画像の細部に亘って類似画像を構成することが可能となり、結果として、検出精度が向上する。また、変形を伴わないため画像の構造が保存され、このことも類似画像生成の精度を向上させている。
【図面の簡単な説明】
【００１３】
【図１】本発明を具体化する画像診断支援装置の概略構成を示す図である。
【図２】図１に例示した画像診断支援装置の動作を示す概略全体フロー図である。
【図３】疑似正常差分画像の生成の詳細を説明するフロー図である。
【図４】最類似画像パッチの決定について説明する図である。
【図５】中心座標の移動について説明する図である。
【図６】最類似画像パッチの合成を説明する図である。
【図７】差分画像のイメージ図である。
【図８】異常陰影候補の検出の動作の詳細を説明するフロー図である。
【図９】異常陰影初期候補の決定をさらに詳細に説明する図である。
【図１０】異常検出実験の各段階の画像を例示する写真である。
【発明を実施するための形態】
【００１４】
以下、例示に基づき本発明を説明する。図１は、本発明を具体化する画像診断支援装置の概略構成を示す図である。なお、以下、人の胸部Ｘ線画像を例として説明するが、本発明は、「胸部」に限らず、腹部（例えば、胃の画像）、頭部（脳の画像）などのあらゆる病巣候補の検出に適用が可能であり、また、人に限らず、人以外の動物を含むあらゆる生体を検査対象とすることができる。
【００１５】
図１において、正常画像データベースには、収集した多数の病巣の無いことが分かっている正常な胸部Ｘ線画像を記憶している。一方、診断対象の患者についての１枚の検査対象画像（ディジタル胸部Ｘ線画像）が画像データ記憶部に記憶される。疑似正常差分画像作成部は、この記憶された検査対象画像と、正常画像データベースの正常な胸部Ｘ線画像から生成した疑似正常画像とを用いて差分画像を作成する。作成した疑似正常差分画像は、疑似正常差分画像記憶部に記憶する。疑似正常差分画像作成部において疑似正常画像を作成するために、検査対象画像から切り出した検査画像パッチに最も類似する類似画像パッチを、データベース画像に設けた探査領域内から類似度計算により決定する。全データベース画像に対して同様に繰り返し、各々のデータベース画像に対応してそれぞれ１つの類似画像パッチを決定する。選出された類似画像パッチの内、最も類似度が高いものを、一つの検査画像パッチに対応する最類似画像パッチとして決定する。同様に、走査により次々に隣接する検査画像パッチに対応する最類似画像パッチの決定を、検査対象画像の全体で行う。
【００１６】
決定した最類似画像パッチを全て用いて全体的に統合し、再構築することにより、検査対象画像に対応する疑似正常画像を生成する。これによって、本人の正常であった時に類似の胸部画像である疑似正常画像を全体的に作り上げる。差分画像作成部は、この疑似正常画像と、本人の検査対象画像との差分画像（疑似正常差分画像）を作成して、疑似正常差分画像記憶部に記憶する。このように、過去に正常と診断された他人の胸部Ｘ線画像データベースから、検査対象の正常画像を擬似的に画像合成し、検査対象画像との差分をとる。これによって、肺がんなどの異常陰影のみが残る。異常陰影検出部は、この異常陰影を疑似正常差分画像中から検出する。異常陰影編集部は、検出した異常陰影を編集して、差分画像の正の階調値に対して線形変換することにより異常陰影の強調表示を行う。画像表示部は、検出しかつ編集した異常陰影を表示する。これによって、肺癌などの肺の病巣を正確に診断することが可能になる。
【００１７】
図２は、図１に例示した画像診断支援装置の動作を示す概略全体フロー図である。ステップＳ１において、診断対象の患者についての１枚のディジタル胸部Ｘ線画像と、正常な胸部Ｘ線画像のみを集めたデータベースから生成した疑似正常画像を用いて疑似正常差分画像を作成する（詳細は図３参照）。ステップS2において、ステップS1で作成した疑似正常差分画像中から異常陰影候補を検出する（詳細は図４，図５参照）。ステップS3において、検出した異常陰影を表示する。
【００１８】
図３は、疑似正常差分画像の生成（図２のステップS1）の詳細を説明するフロー図である。検査対象画像は、検査対象を撮影することにより、１枚の検査対象画像として、例えば、ディジタル胸部Ｘ線画像を取得する。データベース画像として、多数の正常な胸部Ｘ線画像を用意する。このデータベース画像には、肺がんなどの異常が存在していないことが医師により確認されているものを用いる。このとき、データベースに含まれる画像は年齢、性別、体型といった肺部の構造を変化させる要因について、多様性があることが望ましい。したがって、画像数は多いことが望ましいが、収集のコストと精度とのトレードオフにより決定される。
【００１９】
図３のステップS4において、ディジタル胸部Ｘ線画像から胸部領域（検査対象領域）を抽出する。胸部領域外部の階調値は０としておく。ステップS5において、データベース画像に設けた探査領域毎に探索するため（ステップS7参照）、検査対象画像に対して検査画像パッチを切り出す。但し、検査画像パッチは重複部分を持つように切り出す。
【００２０】
ステップS6において、データベースのディジタル胸部Ｘ線画像から胸部領域を抽出する。ステップS7において、ステップS5で切り出した検査画像パッチに最も類似する類似画像パッチを、データベース画像に設けた探査領域内から類似度計算により決定する。全データベース画像に対して同様に繰り返し、各々のデータベース画像に対応してそれぞれ１つの類似画像パッチを決定する。ステップS8において、ステップS7で選出された類似画像パッチの内、最も類似度が高いものを一つの検査画像パッチに対応する最類似画像パッチとして決定する。同様に、走査により次々に隣接する検査画像パッチに対応する最類似画像パッチの決定を、検査対象画像の全体で行う。
【００２１】
さらに、ステップS8の最類似画像パッチの決定について図４を参照して、また、中心座標の移動について図５を参照して説明する。図４において、検査対象画像から、上述したステップS5において、検査画像パッチが切り出される。この検査画像パッチに最も類似する類似画像パッチを、データベース画像に設けた探査領域内から決定するために、探索領域毎に、データベース画像から類似した類似画像パッチの探索を行う。探査領域は検査対象領域内に所定の大きさ（範囲）で設定するが、このとき探査領域の中心は検査画像パッチを切り出した際の中心座標と同じ座標とする。探査領域の大きさは、検査画像パッチよりも大きく（例えば、１０％程度大きく）し、かつ、大きいほど精度向上のために望ましいが、計算コストとのトレードオフとなる。
【００２２】
図４に示すように、画像の左上を原点とし、横方向に x軸、縦方向に y軸を設定する。最初は、原点を中心に検査画像パッチ A（p画素 × q画素）を切り出す。データベース画像から同じく原点を中心に s画素 × t画素の探査領域内（ただし、 p≦ sかつ q ≦ tである）で最も類似する類似画像パッチを p画素 × q画素のサイズで切り出す。類似性の尺度には、ユークリッド距離や相関係数などを用いる。全てのデータベース画像で同様に繰り返し、各々のデータベース画像に対する類似画像パッチを選出する。選出された類似画像パッチのうち、最も類似度が高いものを検査画像パッチの最類似画像パッチ A0とする。
【００２３】
次に、図５に示すx軸方向に、検査画像パッチの中心を k画素 (k < p、 k < q)移動させ、同様に検査画像パッチを切り出し、最類似画像パッチの決定を行う。検査画像パッチの中心が画像の右端まで達したら、y軸方向に k画素移動させ、検査対象画像の右下まで同様に繰り返す。 
【００２４】
類似度は相関係数Ｃによって決定する。相関係数Ｃは以下の式による。
【数１】
2014221163.tif 000003
【００２５】
I(i，j) ： 検査画像パッチの階調値、
T(i，J) ： データベース画像に設けた探査領域から切り出した任意の画像パッチの階調値、
【数２】
2014221163.tif 000004
【００２６】
図３のステップS5〜S8を、検査対象画像全体で繰り返し、ここまでの処理によって、各々の検査画像パッチのそれぞれに対する最類似画像パッチがデータベース内から決定する。
【００２７】
ステップS9において、決定した最類似画像パッチを全て用いて検査対象画像に対応する疑似正常画像を再構築する。検査画像パッチは領域が重複するように切り出されているので、最類似画像パッチも重複した領域が存在する。図６は、最類似画像パッチの合成を説明する図である。上下左右のいずれかに重複した部分は、それぞれ図６に示すように対応する画素の平均階調値を用いて再構築を行う。再構築を行うことにより、検査対象画像に対応する疑似正常画像を生成する。
【００２８】
ステップS１0において、検査対象画像とステップS9で生成した疑似正常画像の差を求め、疑似正常差分画像を作成する。図７は、差分画像のイメージ図である。一般的に肺がんなどの異常陰影は、他の正常部分よりも階調値が高くなっているので、検査対象画像と擬似正常画像の差を求めると、異常陰影が図７のように正の値として算出されると考えられる。図７の検査対象画像は、データベースに含まれる正常な画像に対して、擬似的に肺がん陰影を作成したものである。この検査対象画像において肺がん陰影は右肺中央部に位置している。検査対象画像と生成された擬似正常画像の差分を取ることによって、肺がんなど異常陰影のみが残る。その後、差分画像の正の階調値に対して線形変換することにより異常陰影の強調表示を行う。 
【００２９】
図８は、異常陰影候補の検出（図２のステップS2）の動作の詳細を説明するフロー図である。ステップS11において、上述のようにして生成した疑似正常差分画像に対して背景傾向補正を行う。背景傾向補正とは、おおまかな濃度変化を近似して得た背景傾向を全体の濃度変化から引くことによって撮影体位などによるアーチファクトを取り除き、異常陰影による差分成分を抽出するための処理である。背景傾向補正は次のように行う。まず、おおまかな濃度変化の近似のために、疑似正常差分画像が構成する局面の近似局面を、２次元多項近似を用いて導出する。次に、疑似正常差分画像から近似局面を引き、適当なオフセットを全体に付加して階調値を正の値になるようにする。
【００３０】
ステップS12において、ステップS11で補正した擬似正常差分画像から異常陰影候補の拾い上げを行い、異常陰影初期候補を決定する（さらに詳細な動作は、図９参照）。ステップS13において、補正済み擬似正常差分画像と検査対象画像を用いて異常陰影初期候補の特徴量を計算する。ここで特徴量としては次のようなものがあげられる。連結領域の重心座標（画像座標）である重心位置、連結領域の画素数である面積、実行直径、円形度、周囲長、平均値、標準偏差、分散、最大値、最小値、歪度、尖度、コントラスト、エントロピー、エネルギー等がある。
【００３１】
ステップS14において、ステップS13の特徴量を基にルールベース処理により判定を行い、偽陽性候補を決定する。ルールベース処理とは、医師の読影方法をルール化したものである。解剖学の知識を用いて、本来あるべき場所にない陰影、ないはずの場所にある印影を異常と疑う方法である。例としては、大動脈弓の陰影は除外する、左右対称な陰影は除外する、などである。この偽陽性候補をステップS12で決定した異常陰影初期候補から除去することにより、異常陰影候補を決定する。
【００３２】
図９は、異常陰影初期候補の決定（ステップS12）をさらに詳細に説明する図である。異常陰影初期候補は、周囲よりも階調値が相対的に高い島領域と定義できる。周囲より階調値が高い島領域は、画像空間をx-y平面、階調値をz軸で表現した鳥瞰図でのコブ領域となる。これを検出するために、多重閾値で等高線を引くことにより、コブ領域を検出することができる。まず、ステップS15において、ステップS11で補正した擬似正常差分画像の階調値に関するヒストグラムを算出する。ステップS16において、ステップS15で算出したヒストグラムを使用して多重閾値を導出する。多重閾値とは、累積ヒストグラムの割合に応じて設定された複数の閾値であり、例えば、累積ヒストグラムで、１０％、２０％、・・・、９０％などと閾値を決定していけば９つの閾値が設定できる。ヒストグラムの最高階調値を開始点として境界階調値を移動していき、全体面積に対する累積面積の比率が、予め％値が設定した値に一致する階調値を閾値として導出する。
【００３３】
ステップS17において、ステップS11で補正済みの擬似正常差分画像にステップS16で導出された多重閾値を適用し、2値画像を作成する。ステップS18において、ステップS1７で作成した2値画像に対しラベリングを行い、島領域を抽出する。ステップS19において、ステップS18で抽出した島領域の内、重複した島領域に対しては面積が最も小さいもののみを選択し、それ以外の島領域を除去する。これによって、上述したように、異常陰影初期候補を決定する。
【実施例】
【００３４】
以下の胸部Ｘ線画像を用いて、異常検出実験を行った。
胸部Ｘ線画像 日本放射線技術学会「標準ディジタル画像データベース」
・原画像サイズ 2048×2048px
・検査対象画像 腫瘤画像50枚
・データベース画像 非腫瘤画像87枚
実施時の画像サイズ 512×512 px、画像パッチサイズ
50×50 px、探査領域 80×80 px、移動量 5 px
【００３５】
図１０は、異常検出実験の各段階の画像を例示する写真である。（ａ）は検査対象画像であり、（ｂ）は検査対象画像（ａ）と類似するように作成した疑似正常画像である。（ｃ）は疑似正常画像（ｂ）と、検査対象画像（ａ）との差分画像（疑似正常差分画像）である。（ｄ）は生成した2値画像（図９のステップS17参照）を例示している。（ｅ）は特徴量を用いて偽陽性陰影やノイズを除去して作成した陰影検出画像を例示している。（ｆ）は検出結果の画像を例示している。
【図１】
2014221163.tif 000005
【図２】
2014221163.tif 000006
【図３】
2014221163.tif 000007
【図８】
2014221163.tif 000008
【図９】
2014221163.tif 000009
【図４】
2014221163.tif 000010
【図５】
2014221163.tif 000011
【図６】
2014221163.tif 000012
【図７】
2014221163.tif 000013
【図１０】
2014221163.tif 000014