JP5027872B2 - 異常塊候補を検出する装置、方法、ならびにコンピュータ読取り可能媒体 - Google Patents

異常塊候補を検出する装置、方法、ならびにコンピュータ読取り可能媒体 Download PDF

Info

Publication number
JP5027872B2
JP5027872B2 JP2009506831A JP2009506831A JP5027872B2 JP 5027872 B2 JP5027872 B2 JP 5027872B2 JP 2009506831 A JP2009506831 A JP 2009506831A JP 2009506831 A JP2009506831 A JP 2009506831A JP 5027872 B2 JP5027872 B2 JP 5027872B2
Authority
JP
Japan
Prior art keywords
iris
pixel
map
contrast
filter
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP2009506831A
Other languages
English (en)
Other versions
JP2010500049A (ja
Inventor
ニエ ヤオ
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujifilm Corp
Original Assignee
Fujifilm Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fujifilm Corp filed Critical Fujifilm Corp
Publication of JP2010500049A publication Critical patent/JP2010500049A/ja
Application granted granted Critical
Publication of JP5027872B2 publication Critical patent/JP5027872B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/502—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of breast, i.e. mammography
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30068—Mammography; Breast

Landscapes

  • Engineering & Computer Science (AREA)
  • Quality & Reliability (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Description

本発明は、入力画像内の異常塊候補（ａｂｎｏｒｍａｌ ｍａｓｓ ｃａｎｄｉｄａｔｅｓ）を検出するための装置、方法、およびコンピュータ読取り可能媒体に関する。入力画像内の癌の可能性のある塊は、１つのタイプの異常塊候補である。特に、本発明は、ある種の特徴を有する入力画像内の塊を識別することによって異常塊候補を検出することに関する。
医療分野では、画像に埋め込まれた異常塊候補を自動的に検出し、検出された異常塊候補の画質を向上させ、画質を向上された異常塊候補を含む可視画像を表示する、コンピュータ支援診断（ＣＡＤ）システムが知られている。医師は、ＣＡＤシステムによって検出された異常塊候補を含む可視画像を見て、画像に含まれる異常塊候補が、癌の塊などの患部を表す真の異常塊か否かについて最終的な判断を下す。
たとえばモフォロジカル・フィルタ技術などの異常塊候補を検出する技術が知られている。モフォロジカル・フィルタ技術を使用して、胸部画像にモフォロジカル・フィルタによる画像処理が行われ、モフォロジカル・フィルタの出力値に閾値処理が行われ、微細石灰化の塊（ｍｉｃｒｏｃａｌｃｉｆｉｃａｔｉｏｎ ｍａｓｓ）（異常塊の１つの形態）の候補が自動的に検出される。
減算処理を利用する技術も知られている。減算処理を使用して、入力された医療画像に対応する正常な構造の画像が人工的に形成され、入力された医療画像と正常な構造の画像との差を表す差分画像が形成され、差分画像内の所定値に少なくとも等しい画素値を有する塊が異常塊候補として検出される。
いかなる自動化されたシステムにおいても、正確さが向上することが望ましい。
本発明の一の実施形態による入力画像から異常塊候補を検出する方法は、入力画像に基づいて入力勾配ベクトル・マップを生成するステップと、入力勾配ベクトル・マップに基づいてアイリス・コントラスト（ｉｒｉｓ ｃｏｎｔｒａｓｔ）マップを生成するステップと、アイリス・コントラスト・マップに基づいてアイリス・リング・フィルタ（ｉｒｉｓ ｒｉｎｇ ｆｉｌｔｅｒ）応答マップを生成するステップと、異常塊候補として、画素のアイリス・コントラスト値およびアイリス・リング・フィルタ応答値の両方が、それぞれ最小アイリス・コントラスト閾値以上であり、かつ最小アイリス・リング・フィルタ応答閾値以上になっている、入力画像の画素の位置を出力するステップと、を含む。
入力勾配ベクトル・マップは、入力画像の画素のベクトル値のマップであってもよい。各画素に関するベクトル値は、画素の近傍小領域での入力画像内の画素の変化の方向および大きさを表していてもよい。
アイリス・コントラスト・マップは、入力画像の画素のアイリス・コントラスト値のマップであってもよい。各画素に関するアイリス・コントラスト値は、入力勾配ベクトル・マップ内のアイリス・コントラスト・フィルタに対応する画素の応答値を表していてもよい。
アイリス・リング・フィルタ応答マップは、入力画像の画素のアイリス・リング・フィルタ応答値のマップであってもよい。各画素に関するアイリス・リング・フィルタ応答値は、アイリス・コントラスト・マップ内のアイリス・リング・フィルタに対応する画素の応答値を表していてもよい。
本発明の一の実施形態による異常塊候補の検出装置は、入力画像に基づいて入力勾配ベクトル・マップを生成するように構成された入力勾配ベクトル・マップ生成手段と、入力勾配ベクトル・マップに基づいてアイリス・コントラスト・マップを生成するように構成されたアイリス・コントラスト・マップ生成手段と、アイリス・コントラスト・マップに基づいてアイリス・リング・フィルタ応答マップを生成するように構成されたアイリス・リング・フィルタ応答マップ生成手段と、異常塊候補として、画素のアイリス・コントラスト値およびアイリス・リング・フィルタ応答値の両方が、それぞれ最小アイリス・コントラスト閾値以上であり、かつ最小アイリス・リング・フィルタ応答閾値以上になっている、入力画像の画素の位置を出力するように構成された異常塊候補出力手段と、を備える。
本発明の一の実施形態によるコンピュータ読取り可能な媒体は、入力画像から異常塊候補を検出するための、コンピュータ上で実行可能なプログラムを含む。プログラムは、入力画像に基づいて入力勾配ベクトル・マップを生成するステップと、入力勾配ベクトル・マップに基づいてアイリス・コントラスト・マップを生成するステップと、アイリス・コントラスト・マップに基づいてアイリス・リング・フィルタ応答マップを生成するステップと、異常塊候補として、画素のアイリス・コントラスト値およびアイリス・リング・フィルタ応答値の両方が、それぞれ最小アイリス・コントラスト閾値以上であり、かつ最小アイリス・リング・フィルタ応答閾値以上になっている、入力画像の画素の位置を出力するステップと、を含む。
本発明のこれらの、およびその他の実施形態は、異常塊候補を検出する正確さを向上させる。
本発明の１つの実施形態に従って異常塊候補を検出する方法を示す。 本発明の１つの実施形態に従って異常塊候補を検出する別の方法を示す。 本発明の１つの実施形態に従って入力画像の画素のベクトル値を決定する例を示す。 本発明の１つの実施形態に従って入力画像の画素のベクトル値を決定する例を示す。 本発明の１つの実施形態に従って入力勾配ベクトル・マップの成分マップを生成する方法を示す。 本発明の１つの実施形態に従って入力勾配ベクトル・マップを生成するために使用されるアイリス・コントラスト・フィルタの例示のマスク領域を示す。 本発明の１つの実施形態による、アイリス・コントラスト・フィルタのライン・セグメントと入力画像の画素の勾配との間の例示の角度関係を示す。 本発明の１つの実施形態に従ってアイリス・フィルタ応答マップを生成する方法を示す。 本発明の１つの実施形態に従ってアイリス・フィルタ応答マップを生成する別の方法を示す。 本発明の１つの実施形態に従ってアイリス・コントラスト・マップを生成する別のプロセスを示す。 本発明の１つの実施形態に従って、画素のベクトルの大きさを調整する方法を示す。 本発明の１つの実施形態に従って、識別された異常塊候補を除去する方法を示す。 本発明の１つの実施形態に従って、異常塊の識別プロセスの一貫性を向上させる方法を示す。 本発明の１つの実施形態に従って、入力画像から異常塊候補を検出する装置を示す。
本発明の特徴および利点は、当分野の技術者には図面を参照して以下の説明から明らかになる。
単純化し、かつ例示する目的で、本発明の原理がその例示の実施形態を主に参照することによって説明される。しかし、当業者は、同じ原理が多くのタイプの画像分析システムおよび方法に等しく適用可能であることを容易に理解するであろう。
大まかには、異常塊候補を検出する方法１０は、図１に示されるような２つの概略のステップを含む。最初に、異常塊候補が、乳房Ｘ線撮影画像などの入力画像から識別される（ステップ１０１）。異常塊の１つの例には癌腫瘍がある。異常塊候補が識別された後に、異常塊である可能性が最も大きい候補を除去するために除去プロセスを行うことができる（ステップ１０３）。
図２は、図１のステップ１０１および１０３を実行するための例示の方法２０を示す。図示されるように、入力勾配ベクトル・マップが乳房Ｘ線撮影画像などの入力画像から生成される（ステップ２０１）。入力勾配ベクトル・マップは、入力画像の画素のベクトル値のマップとして表現できる。言い換えれば、入力画像の各画素に対して、入力勾配ベクトル・マップ内に画素の対応するベクトル値がある。
各画素に対して、ベクトル値は、画素の近傍小領域での画素の方向および大きさの両方の変化を表す。ベクトル値は、ベクトルの指向角度（ｏｒｉｅｎｔａｔｉｏｎ ａｎｇｌｅ）φ、およびベクトルの大きさＡの２つの成分を含み、
Figure 0005027872
のように定義できる。
ただし、ΔｘおよびΔｙはそれぞれ、画素の近傍小領域でのｘおよびｙ方向の画素の変化を表す。
図３Ａおよび３Ｂは、画素のベクトル値を決定するための例を示す。入力画像内の各画素ｊに関して、図３Ａに示されるような寸法ｎ×ｎの正方形のマスク・フィルタの中心を、画素ｊに合わせることができる。この例では、ｎは５であり、Δｘは（ｆ１５＋ｆ２５＋ｆ３５＋ｆ４５＋ｆ５５）−（ｆ１１＋ｆ２１＋ｆ３１＋ｆ４１＋ｆ５１）と定義され、Δｙは（ｆ１１＋ｆ１２＋ｆ１３＋ｆ１４＋ｆ１５）−（ｆ５１＋ｆ５２＋ｆ５３＋ｆ５４＋ｆ５５）と定義される。次いで、ベクトルの指向角度（または単に角度）φ、および大きさが、それぞれ式（１）および（２）に従って計算できる。計算されたベクトル値が図３Ｂに示される。
マスク・フィルタの寸法は、図３Ａおよび３Ｂに示されるような例に限定されない。実際に、マスク・フィルタの形も限定されない。その他の形状および寸法のマスク・フィルタが、ベクトル値を決定するために利用できる。
さらに、ベクトルを決定することは、フィルタの使用に厳密には限定されない。ベクトル値が画素の十分に局所化された領域（たとえば近傍小領域）で画素の変化を表すことのみが求められる。本発明の範囲は、説明されたようにベクトル値を生成するために使用できる任意の方法およびシステムを完全に包含する。
上述したように、ベクトル値は角度φおよび大きさＡの２つの成分を含む。したがって、入力勾配ベクトル・マップを生成することは、図４に示されるように２つの成分マップを生成することによって遂行できる。図４では、入力勾配ベクトル・マップを生成するステップ（図２のステップ２０１）は、入力画像から入力勾配角度マップを生成するステップ（ステップ４０１）、および同じく入力画像から入力勾配マグニチュード・マップを生成するステップ（ステップ４０３）によって遂行できる。入力勾配角度マップおよび入力勾配マグニチュード・マップは、それぞれ入力画像の画素の成分の角度および大きさの値を含む。入力勾配角度マップおよび入力勾配マグニチュード・マップが生成される順番は、図４に示されるような実施形態に限定されない。
再び図２を参照すると、入力画像の入力勾配ベクトル・マップが生成された後に、アイリス・コントラスト・マップが入力勾配ベクトル・マップから生成できる（ステップ２０３）。アイリス・コントラスト・マップは、入力画像の画素のアイリス・コントラスト値のマップとして表現できる。画素のアイリス・コントラスト値は、入力勾配ベクトル・マップ内のアイリス・コントラスト・フィルタに対応する画素の応答として表現できる。
図５は、入力勾配ベクトル・マップの画素に適用できるマスク領域を有する例示のアイリス・コントラスト・フィルタを示す。マスク領域は、リング（ハッシュ化された部分）によって特徴付けることができる。アイリス・コントラスト・フィルタのパラメータには、ｒ（リングの内側半径）、ｄ（リングの幅）、Ｒ（リングの外側半径）、Ｉ（リングの外側の限界）、およびＭ（フィルタ処理で検討される角度の数）が含まれる。リングの内側半径ｒは、適応性を有することができ、リングの幅ｄは固定できる。１つの例として、リングの幅ｄは５ピクセルに固定でき、Ｉは２５ピクセルに固定できる。この例では、ｒは、０から２０の範囲の値（０からＩ−ｄ）を取ることができる。また、Ｍは１６に定めることができる。
入力画像内の各画素（ｘ、ｙ）に関するアイリス・コントラスト・マップを生成するために、画素（ｘ、ｙ）のアイリス・コントラスト値Ｃ（ｘ、ｙ）は、
Figure 0005027872
ただし、
Figure 0005027872
のようにアイリス・コントラスト・フィルタを適用することによって定義できる。
式（３）および（４）では、Ａｉｊは入力勾配ベクトル・マップからの画素ｉ，ｊのベクトルの大きさであり、θｉｊは、入力勾配ベクトル・マップからの画素ｉ，ｊのベクトルの方向φｉｊと、画素ｉ，ｊをアイリス・コントラスト・フィルタの中心に連結するライン・セグメントとの間の角度である。あるいは、大きさＡｉｊおよび角度θｉｊの値は、個々の成分の入力勾配マグニチュード・マップおよび入力勾配角度マップから決定できる（図４を参照されたい）。
上記に示されたように、角度θｉｊは、画素ｉ，ｊのベクトルの方向φｉｊと、画素ｉ，ｊをアイリス・コントラスト・フィルタの中心に連結するライン・セグメントとの間の角度として定義される。図５では、そのような１６本のライン・セグメントｊ（Ｍ＝１６）が示される。図６は、式（４）で使用される角度θｉｊと、特定のライン・セグメントｊ上の画素ｉ，ｊのベクトルの方向（配向角度）φｉｊとの間の関係を示す。
また、上記に示されるように、Ｍはフィルタ処理で検討される角度の数を表す。たとえば、Ｍが１６に定められる場合、式（４）を使用して計算されたそれぞれに対して、対象とする画素（ｘ、ｙ）からの半径ｉおよび角度０、±π／１６、±π／８、±３π／１６、±π／４、±５π／１６、±３π／８、±７π／１６、±π／２での画素のみが計算で考慮される。しかし、半径ｉでの全ての画素が、Ｃｉを決定する際に考慮できることが企図される。この例では、Ｍは、対象とする画素から半径ｉにある画素の総数を単に表す。Ｍを特定の数に定めると、計算がより速くなるが、正確さが犠牲になる可能性がある。数Ｍは、正確さの犠牲が、特定の適用に関して無視できる場合に選択できる。
図２を再び参照すると、入力画像のアイリス・コントラスト・マップが生成された後に、アイリス・フィルタ応答マップがアイリス・コントラスト・マップから生成できる（ステップ２０５）。アイリス・フィルタ応答マップは、入力画像の画素のアイリス・フィルタ応答値のマップとして表現できる。画素に関するアイリス・フィルタ応答値は、アイリス・コントラスト・マップ内のアイリス・リング・フィルタに対応する画素の応答として表現できる。
図７Ａは、アイリス・フィルタ応答マップを生成するための例示の過程を示す（図２のステップ２０５）。図示されるように、アイリス・コントラスト勾配角度マップがアイリス・コントラスト・マップから生成できる（ステップ７０１）。アイリス・コントラスト勾配マップは、図３Ａおよび３Ｂに示されるように入力勾配ベクトル・マップを生成するための過程と同様である。しかし、入力画像の代わりに、アイリス・コントラスト・マップがこの過程で入力として使用される。また、勾配角度値（すなわちコントラスト角度）のみが決定される必要がある（式（１）を参照されたい）。
簡単にするため、入力勾配ベクトル・マップを生成するために使用されるマスク・フィルタが、アイリス・コントラスト勾配角度マップを生成するためにも使用できる。しかし、これは厳密には必要でなく、すなわち完全に異なるマスク・フィルタが使用できる。
図７Ａを再び参照すると、画素のアイリス・フィルタ応答値を決定するために、アイリス・フィルタ応答マップが、アイリス・リング・フィルタをアイリス・コントラスト勾配角度マップに適用することによって生成できる（ステップ７０３）。アイリス・リング・フィルタは、図５に示されるアイリス・コントラスト・フィルタと同様なものであることができるが、全く同じである必要はない。言い換えれば、アイリス・リング・フィルタのパラメータｒ、ｄ、Ｒ、Ｉ、およびＭの値は、アイリス・コントラスト・フィルタの対応するパラメータと同じである必要はない。
画素（ｘ、ｙ）のアイリス・フィルタ応答値Ｄ（ｘ、ｙ）は、アイリス・リング・フィルタを
Figure 0005027872
ただし、
Figure 0005027872
のように適用することによって決定できる。
式（５）および（６）では、θｉｊはアイリス・コントラスト勾配角度マップからの画素ｉ，ｊの方向と、画素ｉ，ｊをアイリス・リング・フィルタの中心に連結するライン・セグメントとの間の角度である（図５の説明を参照されたい）。また、Ｍはフィルタ処理で検討される角度の数である。式（３）および（４）に関係する状態と同様に、Ｍは、いくつかの角度のみが検討されることを示すために、任意の特定の値を取ることができる。同じく、全ての角度が同様に検討できる。固定された数の角度を検討するか、全ての角度を検討するかの選択は、計算上の効率および複雑さの問題である。
図２を再び参照すると、アイリス・フィルタ応答マップが生成された後に、異常塊候補位置の位置が出力できる（ステップ２０７）。異常塊候補は、以下の出力条件の両方が真である位置であるように決定できる。第１に、位置のアイリス・コントラスト値（すなわち画素位置）は、最小アイリス・コントラスト閾値以上である。第２に、位置のアイリス・リング・フィルタ応答値が、最小アイリス・リング・フィルタ応答閾値以上である。第１および第２の条件の組合せにより、方法の全体の正確さが向上する。
図２に示される方法は、より効率的にすることができる。たとえば、図７Ｂはアイリス・フィルタ応答マップを生成する別の過程を示す。図示されるように、アイリス・コントラスト勾配角度マップが図７Ａと同じく、アイリス・コントラスト・マップから生成できる（ステップ７０１）。その後、アイリス・コントラスト値が最小アイリス・コントラスト閾値以上である画素（または画素位置）が配置できる（ステップ７１１）。次いで、ステップ７１３で、アイリス・リング・フィルタ応答値が、ステップ７１１で配置された画素に対してのみ決定できる（式（５）および（６）が適用される）。このようにして、第１の出力条件を満たす画素位置のみが評価されるので、過程が速度を上げる、すなわちより効率的にすることができる。
図２に示される方法も代替を含む。上記に示したように、画素（ｘ、ｙ）のアイリス・フィルタ応答値Ｄ（ｘ、ｙ）は、アイリス・リング・フィルタを式（５）および（６）に適用することによって決定できる。しかし、代替として、画素（ｘ、ｙ）の応答値Ｄ（ｘ、ｙ）は、ハーフリング・アイリス・フィルタを
Figure 0005027872
ただし、
Figure 0005027872
のように適用することによって決定できる。
式（７）および（８）では、θｉｊはアイリス・コントラスト勾配角度マップからの画素ｉ，ｊの方向と、画素ｉ，ｊをハーフリング・アイリス・リング・フィルタの中心に連結するライン・セグメントとの間の角度である（図５の説明を参照されたい）。式（７）および（８）は、方法をより効率的にするために最小アイリス・コントラスト閾値以上のアイリス・コントラスト値を有する画素位置のみに適用できることに留意されたい。
別の代替として、方法の正確さは、アイリス・コントラスト・マップを生成する前に、入力勾配ベクトル・マップの大きさを調整することによって向上できる。図８は、アイリス・コントラスト・マップを生成する別の過程を示す。図８では、入力勾配ベクトル・マップのベクトル値の大きさが調整できる（ステップ８０１）。次いで、アイリス・コントラスト・マップが、調整された入力勾配ベクトル・マップから生成できる（ステップ８０３）。
特に、大きさは図９に示されるように調整できる。入力勾配ベクトル・マップの各画素（ｘ、ｙ）に対して、マスク・フィルタの中心を画素に合わせることができる（ステップ９０１）。マスク・フィルタは、図３Ａに示されるフィルタと同様の正方形のフィルタであることができる。しかし、この過程に関するマスク・フィルタの寸法および形状はそれほど限定されない。次いで、最小の大きさＡｍｉｎがフィルタ内の画素に対して確認できる（ステップ９０３）。画素（ｘ、ｙ）の調整された大きさＡｏｕｔは、Ａｘｙ−Ａｍｉｎであるように計算でき、ただし、Ａｘｙは、調整する前の画素ベクトルの大きさである（ステップ９０５）。
そのような調整を行う理由は、スキン・ライン付近の領域では、勾配ベクトルが同様に評価された大きさを有するという観察に基づいている。したがって、ＡｘｙとＡｍｉｎの間の差は小さく、それにより低く評価されたＡｏｕｔが、領域内の偽の勾配ベクトルを抑制するようになる。胸部の中央領域では、胸部組織または癌の塊の画素、および脂肪の領域の画素は、両方とも任意の局所的な領域にある可能性がある。通常は、Ａｍｉｎは脂肪領域の画素に一致し、非常に小さい。したがって、ＡｏｕｔはほぼＡｘｙと同じである。したがって、対象とする領域の所望の勾配情報が保存される。
入力勾配ベクトル・マップの成分マップが使用される場合、入力勾配マグニチュード・マップのみを調整することが求められる可能性がある。
上記の説明では、アイリス・コントラスト・マップ（図２のステップ２０３）は、入力勾配ベクトル・マップの画素に図５のアイリス・コントラスト・フィルタおよび式（３）および（４）を適用することによって生成できる。代替として、アイリス・コントラスト・マップは、下記に示されるように同様のアイリス・コントラスト・フィルタを使用して式
Figure 0005027872
（ただし、Ｃｉ＝ｍｅｄｉａｎ｛Ａｉｊｃｏｓθｉｊ，ｊ＝１，２，．．．，Ｍ｝ （１０））
を適用することによって生成できる。
式（９）および（１０）では、Ａｉｊは入力勾配ベクトル・マップからの画素ｉ，ｊのベクトルの大きさであり、θｉｊは、入力勾配ベクトル・マップからの画素ｉ，ｊのベクトルの方向と、画素ｉ，ｊをアイリス・コントラスト・フィルタの中心に連結するライン・セグメントとの間の角度である。あるいは、大きさＡｉｊおよび角度θｉｊの値は、個々の成分の入力勾配マグニチュード・マップおよび入力勾配角度マップから決定できる。
式（９）および（１０）を適用することから生じるアイリス・コントラスト・マップは、入力の画素のメジアン・アイリス・コントラスト値のアイリス・コントラスト・マップとして、または単にメジアン・アイリス・コントラスト・マップとして表現できる。メジアン・アイリス・コントラスト値は、第１の出力条件が満たされるか否かを決定するために、最小メジアン・アイリス・コントラスト閾値と比較できる。
アイリス・コントラスト出力およびメジアン・アイリス・コントラスト出力の両方が、入力画像、この例では入力勾配ベクトル・マップの円形の縁部に対して感度があり、すなわち高い応答性を有する。対象となる異常塊候補は、一般的に円形を有するものとして特徴付けられ、したがって、円形の縁部に対する感受性が有利であり、システムの正確さがアイリス・コントラスト出力またはメジアン・アイリス・コントラスト出力を利用することによって向上できる。
しかし、メジアン・アイリス・コントラストは、アイリス・コントラスト出力より優れた１つの利点を有する。対象でない通常の塊は、一般的に非円形である。しかし、勾配の大きさが大きい場合、アイリス・コントラスト出力は非円形の塊に対して高い応答を有することができる。一方で、メジアン・アイリス・コントラスト出力は、非円形の塊に対してアイリス・コントラスト出力よりも低い応答を有する。したがって、メジアン・アイリス・コントラスト・マップは、全体的にアイリス・コントラスト・マップよりもノイズが少ない（またはより正確である）。
メジアン・アイリス・コントラスト・マップを生成することの改善として、式（１０）でＣｉを決定するのに使用されるベクトルの大きさＡｉｊは
Figure 0005027872
のように調整できる。
マグニチュード調整式（１１）は、非常に強力な勾配が癌の塊から生じる可能性が少ないという観察に基づいている。値ＡｍａｘおよびＢは、異常塊候補の偽陽性の識別が起こることを最小限に抑えるために、強い勾配を抑制するように選択できる。パラメータＡｍａｘは、典型的な異常塊の縁部が超えない勾配マグニチュード値を表すことができる。バラメータＢは、ガウス関数の広がり（ｓｐｒｅａｄ）を表すことができ、実験的に決定できる。
大きさＡｉｊが式（１１）に従って調整される場合、式（１１）に従って大きさを調整する前に、図８および９に示されるように大きさが調整されることが好ましい。
図２を再び参照すると、アイリス・フィルタ応答マップを生成するステップ（ステップ２０５）では、メジアン・アイリス・コントラスト・マップがアイリス・コントラスト・マップの代わりに使用できる。より具体的には、アイリス・コントラスト勾配マップがメジアン・アイリス・コントラスト・マップから生成でき（図７Ａのステップ７０１）、アイリス・フィルタ応答値がアイリス・コントラスト勾配マップから決定できる（ステップ７０３）。メジアン・アイリス・コントラスト・マップは、図７Ｂに示される向上した効率の過程にも適用される。さらに、ハーフリング・アイリス・リング・フィルタ（式（７）および（８）を参照されたい）がアイリス・リング・フィルタ応答マップを生成するためにメジアン・アイリス・コントラスト・マップに適用できる。
アイリス・リング・フィルタ応答マップが生成された後に、第１および第２の出力条件の両方を満たす、入力画像の１つまたは複数の画素位置が異常塊候補として識別できる（図２のステップ２０７）。正確さおよび有用性は、識別された位置の除去によって向上できる。
図１０は、含むことができるいくつかの除去過程を示す。異常塊候補の可能性として識別される各画素位置に対して、別の画素位置によって異常として既に識別された塊の一部分である場合、それは除去（排除）できる（ステップ１００１）。たとえば、対象となる異常塊候補は、特定の最小寸法になっているものとして特徴付けることができる。したがって、２つの識別された位置が互いの所定の距離の中にある場合、２つの位置が同じ異常塊の部分になっている可能性がある。有用性は、同じ塊を重複して識別する可能性を低減することによって向上できる。
また、エンドユーザは、所定の数の最も可能性のある候補のみに興味を持つ可能性がある。第１および第２の出力条件の両方を満たす全ての位置の中で、異常である位置の可能性に基づいて位置が順序付けられる。尤度は、アイリス・コントラスト値から単独で、アイリス・リング・フィルタ応答値から単独で、または両方の組合せから決定できる。両方が使用される場合、各値は順序を決定するために重み付けすることができる。位置の順序が決定された後に、最も上にある所定の数の位置が、異常塊候補として出力できる（ステップ１００３）。
第１の出力条件を満たす位置があるが、同じ位置が第２の出力条件を満たさない状態がある可能性がある。この場合には、第１の出力条件を満たす１つまたは複数の位置が異常塊候補として識別できるように、第２の出力条件、すなわち最小アイリス・リング・フィルタ応答閾値を緩和できる。
さらに、異常塊を有する可能性が他の領域よりも少ない身体の領域がある。そのような領域の例には、胸壁、肩、スキン・ライン、および胸筋が含まれる。画素位置がこれらの所定の領域のいずれかの中に収まる場合、位置が除去できる（ステップ１００５）。
必ずしも全ての除去ステップ１００１、１００３、および１００５が行われる必要がないことが企図される。また、除去ステップの順序は、図１０に示される順序に限定されない。ステップは、同時または逐次に行うことができ、あるいは全く行われなくてもよい。
別の向上として、最小アイリス・コントラスト閾値が入力画像の特徴に従って設定できる。異なる入力画像が異なる特徴を有することが一般に認識されている。したがって、アイリス・コントラスト値の範囲は、ある入力画像と隣の入力画像とは異なる可能性がある。最小コントラスト閾値が全ての入力画像に関して固定された場合、異常塊候補の識別不足または識別過剰が生じるおそれがある。
異常塊候補の識別の一貫性を向上させる１つの方式は、最小アイリス・コントラスト閾値を入力画像から生成されたアイリス・コントラスト・マップに基づいて調整することである可能性がある。図１１は、この向上された過程を示す。図示されるように、アイリス・コントラスト・マップを生成した後に（ステップ２０３）、およびアイリス・フィルタ応答マップを生成する前に（ステップ２０５）、最小アイリス・コントラスト・レベルが１つまたは複数の入力画像の特徴に基づいて設定できる（ステップ１１０１）ことを除いて、図１１は図２と同一である。
たとえば、所与の画像内で、最上にあるわずかな割合のアイリス・コントラスト値のみが、真に異常である可能性があるとすることができる。たとえば、最も高い５％のアイリス・コントラスト値のみを対象とするこができる。この環境では、最小アイリス・コントラスト・レベルを第９５のパーセンタイル・レベルに設定できる。
別の例として、所定の数の最も可能性のある候補のみが識別できる。次いで、最小アイリス・コントラスト閾値は、アイリス・コントラスト・マップでの所定の数の位置のみが、レベル以上のアイリス・コントラスト値を有するようなレベルに設定できる。
最小アイリス・コントラスト閾値を設定することによって、異常識別の結果においての一貫性が入力画像の範囲全体にわたって達成できる。一貫性をさらに向上させるために、同様の検討が最小アイリス・リング・フィルタ応答閾値を設定するために行われる。
図１２は、入力画像から異常塊候補を検出する装置１２００を示す。装置１２００は、入力勾配ベクトル・マップ生成デバイス１２１０、アイリス・コントラスト・マップ生成デバイス１２２０、アイリス・フィルタ応答マップ生成デバイス１２３０、最小アイリス・コントラスト閾値設定デバイス１２４０、および異常塊候補出力デバイス１２５０を備えることができる。
入力勾配ベクトル・マップ生成デバイス１２１０は、入力画像から入力勾配ベクトル・マップを生成するように構成できる。入力勾配ベクトル・マップ生成デバイス１２１０は、入力勾配ベクトル・マップの成分マップを生成するように構成された、入力勾配角度マップ生成デバイス１２１１、および入力勾配マグニチュード・マップ生成デバイス１２１３を備えることができる。
アイリス・コントラスト・マップ生成デバイス１２２０は、入力勾配ベクトル・マップ生成デバイス１２１０によって生成された入力勾配ベクトル・マップに基づいてアイリス・コントラスト・マップを生成するように構成できる。アイリス・コントラスト・マップ生成デバイス１２２０は、アイリス・コントラスト・フィルタリング・デバイス１２２１、および／またはメジアン・アイリス・コントラスト・フィルタリング・デバイス１２２３を備えることができる。アイリス・コントラスト・フィルタリング・デバイス１２２１は、式（３）および（４）に基づいて画素のアイリス・コントラスト応答値を決定するように構成でき、メジアン・アイリス・コントラスト・フィルタリング・デバイス１２２３は、式（９）および（１０）に基づいて画素のメジアン・アイリス・コントラスト応答値を決定するように構成できる。アイリス・コントラスト・マップ生成デバイス１２２０は、図８および９に示され、ならびに／または式（１１）による過程の結果に従って入力勾配ベクトル・マップの大きさを調整するように構成された、勾配マグニチュード調整デバイス１２２５をさらに備えることができる。
アイリス・フィルタ応答マップ生成デバイス１２３０は、アイリス・コントラスト・マップ生成デバイス１２２０によって生成された入力コントラスト・マップに基づいてアイリス・フィルタ応答マップを生成するように構成できる。アイリス・フィルタ応答マップ生成デバイス１２３０は、アイリス・リング・フィルタリング・デバイス１２３１、および／またはハーフリング・アイリス・フィルタリング・デバイス１２３３を備えることができる。アイリス・リング・フィルタリング・デバイス１２３１は、式（５）および（６）に基づいて画素のアイリス・リング・フィルタ応答値を決定するように構成でき、ハーフリング・アイリス・フィルタリング・デバイス１２３３は式（７）および（８）に基づいて画素のハーフリング・アイリス・リング・フィルタ応答値を決定するように構成できる。アイリス・フィルタ応答マップ生成デバイス１２３５は、アイリス・フィルタ応答マップを生成するために、アイリス・リング・フィルタリング・デバイス１２３１および／またはハーフリング・アイリス・フィルタリング・デバイス１２３３に対する入力として使用できるアイリス・コントラスト勾配角度マップを生成するように構成されたアイリス・コントラスト勾配角度マップ生成デバイスをさらに備えることができる。
最小アイリス・コントラスト閾値設定デバイス１２４０は、図１１のステップ１１０１の結果に従ってアイリス・コントラスト・マップ生成デバイス１２２０によって生成される入力コントラスト・マップに基づいて最小アイリス・コントラスト閾値を設定するように構成できる。
異常塊候補出力デバイス１２５０は、ステップ２０７に従って異常塊候補として第１および第２の出力条件を満たす画素の位置を出力するように構成できる（図２および１１を参照されたい）。異常塊候補出力デバイス１２５０は、ステップ１００１、１００３、および１００５のうちの１つまたは複数の結果に従って候補の位置を除去するように構成された除去デバイス１２５１を備えることができる（図１０を参照されたい）。
異常塊候補を識別する方法は、実行された場合にコンピュータ・デバイスが異常塊候補を識別できるように、計算可能−読取り可能媒体（ｃｏｍｐｕｔａｂｌｅ−ｒｅａｄａｂｌｅ ｍｅｄｉｕｍ）でのプログラムとして記録できる。
本発明をその例示の実施形態を参照して説明してきたが、当業者は、本発明の真の趣旨および範囲から逸脱することなく本発明の説明された実施形態に様々な変形を行うことができる。本明細書に使用される用語および説明は、例示として明記されるに過ぎず、限定するものではない。特に、本発明の方法は例として説明してきたが、方法のステップを例示されたものとは異なる順序で、または同時に行うことができる。当業者は、これらの、およびその他の変形形態が、以下の特許請求の範囲に定義された本発明の趣旨および範囲、ならびにその等価物の中で可能であることを理解するであろう。

Claims (26)

  1. 入力画像に基づいて入力勾配ベクトル・マップを生成するステップであって、前記入力勾配ベクトル・マップは前記入力画像の画素のベクトル値のマップであり、各画素に関する前記ベクトル値が前記画素の近傍小領域で前記入力画像内の前記画素の変化の方向および大きさを表すステップと、
    前記入力勾配ベクトル・マップに基づいてアイリス・コントラスト・マップを生成するステップであって、前記アイリス・コントラスト・マップは前記入力画像の前記画素のアイリス・コントラスト値のマップであり、各画素に関する前記アイリス・コントラスト値が前記入力勾配ベクトル・マップ内のアイリス・コントラスト・フィルタに対応する画素の応答値を表すステップと、
    前記アイリス・コントラスト・マップに基づいてアイリス・リング・フィルタ応答マップを生成するステップであって、前記アイリス・リング・フィルタ応答マップは前記入力画像の前記画素のアイリス・リング・フィルタ応答値のマップであり、各画素に関する前記アイリス・リング・フィルタ応答値が、前記アイリス・コントラスト・マップ内のアイリス・リング・フィルタに対応する画素の応答値を表すステップと、
    異常塊候補として、前記画素の前記アイリス・コントラスト値および前記アイリス・リング・フィルタ応答値の両方が、それぞれ最小アイリス・コントラスト閾値以上であり、かつ最小アイリス・リング・フィルタ応答閾値以上になっている、前記入力画像の画素の位置を出力するステップとを含む、
    入力画像から異常塊候補を検出する方法。
  2. 前記入力勾配ベクトル・マップを生成するステップが、
    前記入力画像に基づいて入力勾配角度マップを生成するステップであって、前記入力勾配角度マップが前記入力画像の前記画素の角度のマップであり、各画素に関する前記角度が前記画素の近傍小領域での前記入力画像内の前記画素の変化の方向を表すステップと、
    前記入力画像に基づいて入力勾配マグニチュード・マップを生成するステップであって、前記入力勾配マグニチュード・マップが前記入力画像の前記画素のスカラー値のマップであり、各画素に関する前記スカラー値が前記画素の近傍小領域での前記入力画像内の前記画素の変化の大きさを表すステップと、を含む請求項１に記載の方法。
  3. 前記アイリス・コントラスト・マップを生成するステップが、
    前記入力画像の各画素位置に対して
    Figure 0005027872
    Figure 0005027872
    のように前記アイリス・コントラスト・フィルタを適用することによって応答Ｃ（ｘ，ｙ）を出力するステップを含み、
    Ｃ（ｘ，ｙ）が画素位置の前記アイリス・コントラスト値であり、Ａｉｊが前記入力勾配ベクトル・マップからの画素ｉ，ｊのベクトルの大きさであり、θｉｊが前記入力勾配ベクトル・マップからの前記画素ｉ，ｊのベクトルの方向と、前記画素ｉ，ｊを前記アイリス・コントラスト・フィルタの中心に連結するライン・セグメントとの間の角度であり、ｒが前記アイリス・コントラスト・フィルタの内側半径であり、ｄが前記アイリス・コントラスト・フィルタのリングの幅であり、ＲがＲ＝ｒ＋ｄであるように前記アイリス・コントラスト・フィルタの外側半径であり、ＩがＲの上限であり、Ｍが方向の数であり、
    前記アイリス・コントラスト・フィルタの中心が前記画素位置（ｘ，ｙ）に合わせられ、ｒが適応性を有し、ｄが固定される、請求項１に記載の方法。
  4. 前記アイリス・リング・フィルタ応答マップを生成するステップが、
    前記アイリス・コントラスト・マップに基づいてアイリス・コントラスト勾配角度マップを生成するステップであって、前記アイリス・コントラスト勾配角度マップが前記入力画像の前記画素のコントラスト角度のマップであり、各画素に関する前記コントラスト角度が、前記対応する画素の近傍小領域での前記アイリス・コントラスト・マップ内の対応する画素の変化の方向を表すステップと、
    前記アイリス・リング・フィルタを各画素位置に関して
    Figure 0005027872
    もしくは
    Figure 0005027872
    のように適用し、
    ハーフリング・アイリス・リング・フィルタを各画素位置に関して
    Figure 0005027872
    のように適用するステップのうちの１つを介して応答Ｄ（ｘ，ｙ）を出力するステップと、を含み、
    Ｄ（ｘ，ｙ）が画素位置の前記アイリス・フィルタ応答値であり、θｉｊが前記アイリス・コントラスト勾配角度マップからの画素ｉ，ｊのコントラスト角度と、前記画素ｉ，ｊを前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタの中心に連結するライン・セグメントとの間の角度であり、ｒが前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタの内側半径であり、ｄが前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタのリングの幅であり、ＲがＲ＝ｒ＋ｄであるように前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタの外側半径であり、ＩがＲの上限であり、Ｍが方向の数であり、
    前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタの中心が前記画素位置（ｘ，ｙ）に合わせられ、
    ｒが適応性を有し、ｄが固定される、請求項３に記載の方法。
  5. 前記アイリス・リング・フィルタ応答マップを生成するステップがさらに、
    そのアイリス・コントラスト値のそれぞれが、前記最小アイリス・コントラスト閾値以上である１つまたは複数の画素位置を選択するステップと、
    前記応答Ｄ（ｘ，ｙ）を前記選択された画素位置に対してのみ出力するステップと、を含む、請求項４に記載の方法。
  6. 前記アイリス・コントラスト・マップを生成するステップがさらに、
    前記アイリス・コントラスト・フィルタを適用する前に、前記入力勾配ベクトル・マップを調整するステップであって、
    前記入力勾配ベクトル・マップを調整するステップが、前記入力勾配ベクトル・マップの各画素（ｘ，ｙ）に関して、
    所定の寸法のマスク・フィルタの中心を前記画素（ｘ，ｙ）に合わせるステップと、
    Ａｍｉｎを決定するステップであって、Ａｍｉｎが前記入力勾配ベクトル・マップからの前記マスク・フィルタ内での前記画素の最小の大きさであるステップと、
    調整された大きさＡｏｕｔを前記画素に関してＡｏｕｔ＝Ａｘｙ−Ａｍｉｎであるように出力するステップであって、Ａｘｙが前記入力勾配ベクトル・マップからの前記画素（ｘ，ｙ）の大きさであるステップと、を含む、請求項３に記載の方法。
  7. 前記アイリス・コントラスト・マップを生成するステップが、前記アイリス・コントラスト・フィルタを適用することによって応答Ｃ（ｘ，ｙ）を前記入力画像の各画素の位置に対して
    Figure 0005027872
     Ｃｉ＝ｍｅｄｉａｎ｛Ａｉｊｃｏｓθｉｊ，ｊ＝１，２，．．．，Ｍ｝
    のように出力するステップを含み、
    Ｃ（ｘ，ｙ）が画素位置の前記アイリス・コントラスト値であり、Ａｉｊが前記入力勾配ベクトル・マップからの画素ｉ，ｊのベクトルの大きさであり、θｉｊが前記入力勾配ベクトル・マップからの前記画素ｉ，ｊのベクトルの方向と、前記画素ｉ，ｊを前記アイリス・コントラスト・フィルタの中心に連結するライン・セグメントとの間の角度であり、ｒが前記アイリス・コントラスト・フィルタの内側半径であり、ｄが前記アイリス・コントラスト・フィルタのリングの幅であり、ＲがＲ＝ｒ＋ｄであるように前記アイリス・コントラスト・フィルタの外側半径であり、ＩがＲの上限であり、Ｍが方向の数であり、
    前記アイリス・コントラスト・フィルタの中心が前記画素位置（ｘ，ｙ）に合わせられ、
    ｒが適応性を有し、ｄが固定される、請求項１に記載の方法。
  8. 前記アイリス・コントラスト・マップを生成するステップがさらに、
    出力Ｃ（ｘ，ｙ）に前記アイリス・コントラスト・フィルタを適用する前に、
    Figure 0005027872
    であるように、前記入力勾配ベクトル・マップから前記画素ｉ，ｊの前記ベクトルの大きさＡｉｊを調整するステップを含み、
    Ａｍａｘが所定の最大マグニチュード値であり、Ｂが所定の除数である、請求項７に記載の方法。
  9. 前記アイリス・リング・フィルタ応答マップを生成するステップが、
    前記アイリス・コントラスト・マップに基づいてアイリス・コントラスト勾配角度マップを生成するステップであって、前記アイリス・コントラスト勾配角度マップが前記入力画像の前記画素のコントラスト角度のマップであり、各画素に関する前記コントラスト角度が、前記対応する画素の近傍小領域での前記アイリス・コントラスト・マップ内の対応する画素の変化の方向を表すステップと、
    前記アイリス・リング・フィルタを各画素位置に関して
    Figure 0005027872
    もしくは
    Figure 0005027872
    のように適用し、または
    ハーフリング・アイリス・リング・フィルタを各画素位置に関して
    Figure 0005027872
    のように適用するステップのうちの１つを介して応答Ｄ（ｘ，ｙ）を出力するステップとを含み、
    Ｄ（ｘ，ｙ）が画素位置の前記アイリス・フィルタ応答値であり、θｉｊが前記アイリス・コントラスト勾配角度マップからの画素ｉ，ｊの前記コントラスト角度と、前記画素ｉ，ｊを前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタの中心に連結するライン・セグメントとの間の角度であり、ｒが前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタの内側半径であり、ｄが前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタのリングの幅であり、ＲがＲ＝ｒ＋ｄであるように前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタの外側半径であり、ＩがＲの上限であり、Ｍが方向の数であり、
    前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタの中心が前記画素位置（ｘ，ｙ）に合わせられ、
    ｒが適応性を有し、ｄが固定される、請求項７に記載の方法。
  10. 前記アイリス・リング・フィルタ応答マップを生成するステップがさらに、
    そのアイリス・コントラスト値のそれぞれが、前記最小アイリス・コントラスト閾値以上である１つまたは複数の画素位置を選択するステップと、
    前記応答Ｄ（ｘ，ｙ）を前記選択された画素位置に対してのみ出力するステップとを含む、請求項９に記載の方法。
  11. 前記アイリス・コントラスト・マップを生成するステップがさらに、
    前記アイリス・コントラスト・フィルタを適用する前に前記入力勾配ベクトル・マップを調整するステップを含み、
    前記入力勾配ベクトル・マップを調整するステップは、前記入力勾配ベクトル・マップの各画素（ｘ，ｙ）に関して、
    所定の寸法のマスク・フィルタの中心を前記画素（ｘ，ｙ）に合わせるステップと、
    Ａｍｉｎを決定するステップであって、Ａｍｉｎが前記入力勾配ベクトル・マップからの前記マスク・フィルタ内での前記画素の最小の大きさであるステップと、
    調整された大きさＡｏｕｔを前記画素に関してＡｏｕｔ＝Ａｘｙ−Ａｍｉｎであるように出力するステップであって、Ａｘｙが前記入力勾配ベクトル・マップからの前記画素（ｘ，ｙ）の大きさであるステップと、
    を含む、請求項７に記載の方法。
  12. 前記異常塊候補として前記位置を出力するステップが、
    前記位置が、前記入力画像の別の位置によって異常であるとして既に識別された塊の一部分であるかどうか決定するステップと、
    前記位置が前記入力画像の前記別の位置によって既に異常として識別された前記塊の一部分でないことが決定された場合に、前記位置を前記異常塊候補として出力するステップとを含む、請求項１に記載の方法。
  13. 前記位置が既に異常であるとして識別された前記塊の一部分であるかどうか決定するステップが、
    前記位置が前記別の位置からの最小閾値の距離内であるかどうか決定するステップと、
    前記位置が前記別の位置からの前記最小閾値の距離内にあることが決定された場合に、前記位置が前記既に識別された塊の一部分であることを決定するステップと、を含む、請求項１２に記載の方法。
  14. 前記位置が複数の異常塊候補のうちの１つであり、
    前記位置を前記異常塊候補として出力するステップが、
    前記アイリス・コントラスト・マップ値および前記位置に対応する前記アイリス・リング・フィルタ応答値の一方または両方に基づいて前記複数の前記異常塊候補の間の前記位置の順序を決定するステップと、
    前記位置の前記順序が候補の所定の最大数以内である場合に、前記位置を前記異常塊候補として出力するステップと、
    を含む、請求項１に記載の方法。
  15. 前記異常塊候補として前記位置を出力するステップが、
    前記位置が所定の領域内にあるかどうか決定するステップであって、前記所定の領域が胸壁領域、肩領域、スキン・ライン領域、および胸筋領域のうちの少なくとも１つであるステップと、
    前記位置が前記所定の領域内にないことが決定された場合に、前記位置を前記異常塊候補として出力するステップと、
    を含む、請求項１に記載の方法。
  16. 前記入力画像の所定の数の位置のみ、または所定の割合の位置のみが前記コントラスト値以上のアイリス・コントラスト値を有するように、前記アイリス・コントラスト・マップからコントラスト値レベルを決定するステップと、
    前記コントラスト値を前記最小アイリス・コントラスト閾値として設定するステップとをさらに含む、請求項１に記載の方法。
  17. 入力画像に基づいて入力勾配ベクトル・マップを生成するように構成された入力勾配ベクトル・マップ生成手段であって、前記入力勾配ベクトル・マップが前記入力画像の画素のベクトル値のマップであり、各画素に関する前記ベクトル値が前記画素の近傍小領域で前記入力画像内の前記画素の変化の方向および大きさを表す入力勾配ベクトル・マップ生成手段と、
    前記入力勾配ベクトル・マップに基づいてアイリス・コントラスト・マップを生成するように構成されたアイリス・コントラスト・マップ生成手段であって、前記アイリス・コントラスト・マップが前記入力画像の前記画素のアイリス・コントラスト値のマップであり、各画素に関する前記アイリス・コントラスト値が前記入力勾配ベクトル・マップ内のアイリス・コントラスト・フィルタに対応する画素の応答値を表すアイリス・コントラスト・マップ生成手段と、
    前記アイリス・コントラスト・マップに基づいてアイリス・リング・フィルタ応答マップを生成するように構成されたアイリス・リング・フィルタ応答マップ生成手段であって、前記アイリス・リング・フィルタ応答マップが前記入力画像の前記画素のアイリス・リング・フィルタ応答値のマップであり、各画素に関する前記アイリス・リング・フィルタ応答値が前記アイリス・コントラスト・マップ内のアイリス・リング・フィルタに対応する画素の応答値を表すアイリス・リング・フィルタ応答マップ生成手段と、
    前記異常塊候補として、前記画素の前記アイリス・コントラスト値および前記アイリス・リング・フィルタ応答値の両方が、それぞれ最小アイリス・コントラスト閾値以上であり、かつ最小アイリス・リング・フィルタ応答閾値以上になっている、前記入力画像の画素の位置を出力するように構成された異常塊候補出力手段とを備える、異常塊候補の検出装置。
  18. 前記アイリス・コントラスト・マップ生成手段が、
    前記アイリス・コントラスト・フィルタを適用することによって応答Ｃ（ｘ，ｙ）を前記入力画像の各画素位置に対して
    Figure 0005027872
    のように出力するように構成されたアイリス・コントラスト・フィルタリング手段であって、
    Ｃ（ｘ，ｙ）が画素位置の前記アイリス・コントラスト値であり、Ａｉｊが前記入力勾配ベクトル・マップからの画素ｉ，ｊのベクトルの大きさであり、θｉｊが前記入力勾配ベクトル・マップからの前記画素ｉ，ｊのベクトルの方向と、前記画素ｉ，ｊを前記アイリス・コントラスト・フィルタの中心に連結するライン・セグメントとの間の角度であり、ｒが前記アイリス・コントラスト・フィルタの内側半径であり、ｄが前記アイリス・コントラスト・フィルタのリングの幅であり、ＲがＲ＝ｒ＋ｄであるように前記アイリス・コントラスト・フィルタの外側半径であり、ＩがＲの上限であり、Ｍが方向の数であり、
    前記アイリス・コントラスト・フィルタの中心が前記画素位置（ｘ，ｙ）に合わせられ、
    ｒが適応性を有し、ｄが固定される、請求項１７に記載の装置。
  19. 前記アイリス・リング・フィルタ応答マップ生成手段が、
    前記アイリス・コントラスト・マップに基づいてアイリス・コントラスト勾配角度マップを生成するように構成されたアイリス・コントラスト勾配角度マップ生成手段であって、前記アイリス・コントラスト勾配角度マップが前記入力画像の前記画素のコントラスト角度のマップであり、各画素に関する前記コントラスト角度が、前記対応する画素の近傍小領域での前記アイリス・コントラスト・マップ内の対応する画素の変化の方向を表すアイリス・コントラスト勾配角度マップ生成手段と、
    前記アイリス・フィルタリング手段、またはハーフリング・アイリス・リング・フィルタリング手段、もしくは両方であって、前記アイリス・フィルタリング手段が各画素位置に関してアイリス・リング・フィルタを適用することによって応答Ｄ（ｘ，ｙ）を
    Figure 0005027872
    であるように出力するように構成され、
    前記ハーフリング・アイリス・リング・フィルタリング手段が各画素位置に関してハーフリング・アイリス・リング・フィルタを適用することによって応答Ｄ（ｘ，ｙ）を
    Figure 0005027872
    であるように出力するように構成され、
    Ｄ（ｘ，ｙ）が画素位置の前記アイリス・フィルタ応答値であり、θｉｊが前記アイリス・コントラスト勾配角度マップからの画素ｉ，ｊの前記コントラスト角度と、前記画素ｉ，ｊを前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタの中心に連結するライン・セグメントとの間の角度であり、ｒが前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタの内側半径であり、ｄが前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタのリングの幅であり、ＲがＲ＝ｒ＋ｄであるように前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタの外側半径であり、ＩがＲの上限であり、Ｍが方向の数であり、
    前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタの中心が前記画素位置（ｘ，ｙ）に合わせられ、
    ｒが適応性を有し、ｄが固定される、前記アイリス・フィルタリング手段、またはハーフリング・アイリス・リング・フィルタリング手段、もしくは両方、を備える、請求項１８に記載の装置。
  20. 前記アイリス・コントラスト・マップ生成手段がさらに、
    前記入力勾配ベクトル・マップ生成手段から前記入力勾配ベクトル・マップを調整するように構成された勾配マグニチュード調整手段を備え、
    前記調整された入力勾配ベクトル・マップが前記アイリス・コントラスト・フィルタリング手段に提供され、
    前記勾配マグニチュード調整手段が、前記入力勾配ベクトル・マップの各画素（ｘ，ｙ）に関して、
    所定の寸法のマスク・フィルタの中心を前記画素（ｘ，ｙ）に合わせるステップと、
    Ａｍｉｎを決定するステップであって、Ａｍｉｎが前記入力勾配ベクトル・マップからの前記マスク・フィルタ内での前記画素の最小の大きさであるステップと、
    調整された大きさＡｏｕｔを前記画素に関してＡｏｕｔ＝Ａｘｙ−Ａｍｉｎであるように出力するステップであって、Ａｘｙが前記入力勾配ベクトル・マップからの前記画素（ｘ，ｙ）の大きさであるステップとを行うことによって前記入力勾配ベクトル・マップを調整するように構成される、請求項１８に記載の装置。
  21. 前記アイリス・コントラスト・マップ生成手段が、
    前記アイリス・コントラスト・フィルタを適用することによって応答Ｃ（ｘ，ｙ）を前記入力画像の各画素位置に対して
    Figure 0005027872
     Ｃｉ＝ｍｅｄｉａｎ｛Ａｉｊｃｏｓθｉｊ，ｊ＝１，２，．．．，Ｍ｝
    であるように出力するように構成されたメジアン・アイリス・コントラスト・フィルタリング手段であって、
    Ｃ（ｘ，ｙ）が画素位置の前記アイリス・コントラスト値であり、Ａｉｊが前記入力勾配ベクトル・マップからの画素ｉ，ｊのベクトルの大きさであり、θｉｊが前記入力勾配ベクトル・マップからの前記画素ｉ，ｊのベクトルの方向と、前記画素ｉ，ｊを前記アイリス・コントラスト・フィルタの中心に連結するライン・セグメントとの間の角度であり、ｒが前記アイリス・コントラスト・フィルタの内側半径であり、ｄが前記アイリス・コントラスト・フィルタのリングの幅であり、ＲがＲ＝ｒ＋ｄであるように前記アイリス・コントラスト・フィルタの外側半径であり、ＩがＲの上限であり、Ｍが方向の数であり、
    前記アイリス・コントラスト・フィルタの中心が前記画素位置（ｘ，ｙ）に合わせられ、
    ｒが適応性を有し、ｄが固定される、請求項１７に記載の装置。
  22. 前記アイリス・コントラスト・マップ生成手段がさらに、
    前記入力勾配ベクトル・マップ生成手段から前記入力勾配ベクトル・マップを調整するように構成された勾配マグニチュード調整手段を備え、
    前記調整された入力勾配ベクトル・マップが前記メジアン・アイリス・コントラスト・フィルタリング手段に提供され、
    前記勾配マグニチュード調整手段が、前記入力勾配ベクトル・マップの各画素ｉ，ｊに関して、前記入力勾配ベクトル・マップから前記画素ｉ，ｊの前記ベクトルの大きさＡｉｊを調整するステップを
    Figure 0005027872
    であり、ただし、Ａｍａｘが所定の最大マグニチュード値であり、Ｂが所定の除数であるように、前記アイリス・コントラスト・フィルタを出力Ｃ（ｘ，ｙ）を適用する前に行うことによって前記入力勾配ベクトル・マップを調整するように構成された、請求項２１に記載の装置。
  23. 前記アイリス・フィルタ応答マップ生成手段が、
    前記アイリス・コントラスト・マップに基づいてアイリス・コントラスト勾配角度マップを生成するように構成された、アイリス・コントラスト勾配角度マップ生成手段であって、前記アイリス・コントラスト勾配角度マップが前記入力画像の前記画素のコントラスト角度のマップであり、各画素に関する前記コントラスト角度が、前記対応する画素の近傍小領域での前記アイリス・コントラスト・マップ内の対応する画素の変化の方向を表すステップと、
    アイリス・リング・フィルタリング手段、またはハーフリング・アイリス・リング・フィルタリング手段、もしくは両方であって、前記アイリス・フィルタリング手段が各画素位置に関してアイリス・リング・フィルタを適用することによって応答Ｄ（ｘ，ｙ）を
    Figure 0005027872
    であるように出力するように構成され、
    前記ハーフリング・アイリス・リング・フィルタリング手段が各画素位置に関してハーフリング・アイリス・リング・フィルタを適用することによって前記応答Ｄ（ｘ，ｙ）を
    Figure 0005027872
    であるように出力するように構成され、
    Ｄ（ｘ，ｙ）が画素位置の前記アイリス・フィルタ応答値であり、θｉｊが前記アイリス・コントラスト勾配角度マップからの画素ｉ，ｊの前記コントラスト角度と、前記画素ｉ，ｊを前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタの中心に連結するライン・セグメントとの間の角度であり、ｒが前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタの内側半径であり、ｄが前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタのリングの幅であり、ＲがＲ＝ｒ＋ｄであるように前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタの外側半径であり、ＩがＲの上限であり、Ｍが方向の数であり、
    前記アイリス・リング・フィルタまたは前記ハーフリング・アイリス・リング・フィルタの中心が前記画素位置（ｘ，ｙ）に合わせられ、
    ｒが適応性を有し、ｄが固定される、請求項２１に記載の装置。
  24. 前記アイリス・コントラスト・マップ生成手段がさらに、
    前記入力勾配ベクトル・マップ生成手段から前記入力勾配ベクトル・マップを調整するように構成された勾配マグニチュード調整手段を備え、
    前記調整された入力勾配ベクトル・マップが前記アイリス・コントラスト・フィルタリング手段に提供され、
    前記調整された入力勾配ベクトル・マップが前記アイリス・コントラスト・フィルタリング手段に提供され、前記勾配マグニチュード調整手段が、前記入力勾配ベクトル・マップの各画素（ｘ，ｙ）に関して、
    所定の寸法のマスク・フィルタの中心を前記画素（ｘ，ｙ）に合わせるステップと、
    Ａｍｉｎを決定するステップであって、Ａｍｉｎが前記入力勾配ベクトル・マップからの前記マスク・フィルタ内での前記画素の最小の大きさであるステップと、
    調整された大きさＡｏｕｔを前記画素に関してＡｏｕｔ＝Ａｘｙ−Ａｍｉｎであるように出力するステップであって、Ａｘｙが前記入力勾配ベクトル・マップからの前記画素（ｘ，ｙ）の大きさであるステップとを行うことによって、前記入力勾配ベクトル・マップを調整するように構成される、請求項２１に記載の装置。
  25. 前記異常塊候補出力手段が、前記異常塊候補として、前記画素の前記アイリス・コントラスト応答値が、前記最小アイリス・リング・フィルタ応答閾値以上になっているアイリス・リング・フィルタ応答値を有する位置が全く存在しない場合に、前記最小アイリス・コントラスト閾値以上である前記入力画像の画素の前記位置を出力するように構成された、請求項１７に記載の装置。
  26. 入力画像から異常塊候補を検出するコンピュータで実行可能なプログラムが記録されるコンピュータ読取り可能な媒体であって、前記プログラムが、
    前記入力画像に基づいて入力勾配ベクトル・マップを生成するステップであって、前記入力勾配ベクトル・マップが前記入力画像の画素のベクトル値のマップであり、各画素に関する前記ベクトル値が前記画素の近傍小領域で前記入力画像内の前記画素の変化の方向および大きさを表すステップと、
    前記入力勾配ベクトル・マップに基づいてアイリス・コントラスト・マップを生成するステップであって、前記アイリス・コントラスト・マップが前記入力画像の前記画素のアイリス・コントラスト値のマップであり、各画素に関する前記アイリス・コントラスト値が前記入力勾配ベクトル・マップ内のアイリス・コントラスト・フィルタに対応する画素の応答値を表すステップと、
    前記アイリス・コントラスト・マップに基づいてアイリス・リング・フィルタ応答マップを生成するステップであって、前記アイリス・リング・フィルタ応答マップが前記入力画像の前記画素のアイリス・リング・フィルタ応答値のマップであり、各画素に関する前記アイリス・リング・フィルタ応答値が、前記アイリス・コントラスト・マップ内のアイリス・リング・フィルタに対応する画素の応答値を表すステップと、
    前記異常塊候補として、前記画素の前記アイリス・コントラスト値および前記アイリス・リング・フィルタ応答値の両方が、それぞれ最小アイリス・コントラスト閾値以上であり、かつ最小アイリス・リング・フィルタ応答閾値以上になっている、前記入力画像の画素の位置を出力するステップとを含む、コンピュータ読取り可能な媒体。
JP2009506831A 2006-08-16 2007-08-15 異常塊候補を検出する装置、方法、ならびにコンピュータ読取り可能媒体 Expired - Fee Related JP5027872B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US11/504,630 2006-08-16
US11/504,630 US20080044070A1 (en) 2006-08-16 2006-08-16 Abnormal mass candidate detecting apparatus, method and computer-readable medium
PCT/JP2007/066217 WO2008020648A1 (en) 2006-08-16 2007-08-15 Abnormal mass candidate detecting apparatus, method and computer-readable medium

Publications (2)

Publication Number Publication Date
JP2010500049A JP2010500049A (ja) 2010-01-07
JP5027872B2 true JP5027872B2 (ja) 2012-09-19

Family

ID=39082168

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009506831A Expired - Fee Related JP5027872B2 (ja) 2006-08-16 2007-08-15 異常塊候補を検出する装置、方法、ならびにコンピュータ読取り可能媒体

Country Status (3)

Country Link
US (1) US20080044070A1 (ja)
JP (1) JP5027872B2 (ja)
WO (1) WO2008020648A1 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4498149B2 (ja) * 2005-01-17 2010-07-07 キヤノン株式会社 画像読み取り装置
US8280119B2 (en) * 2008-12-05 2012-10-02 Honeywell International Inc. Iris recognition system using quality metrics
CN106097327B (zh) * 2016-06-06 2018-11-02 宁波大学 结合流形特征及双目特性的立体图像质量客观评价方法
CN111667467B (zh) * 2020-05-28 2021-01-26 江苏大学附属医院 基于聚类算法的下肢血管钙化指数多参数累积计算方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5627907A (en) * 1994-12-01 1997-05-06 University Of Pittsburgh Computerized detection of masses and microcalcifications in digital mammograms
JPH08263641A (ja) * 1995-01-23 1996-10-11 Fuji Photo Film Co Ltd アイリスフィルターの演算方法
JP2004248817A (ja) * 2003-02-19 2004-09-09 Fuji Photo Film Co Ltd 異常陰影検出方法および装置並びにプログラム
JP4453321B2 (ja) * 2003-09-26 2010-04-21 コニカミノルタエムジー株式会社 医用画像処理装置及びプログラム
JP2005198890A (ja) * 2004-01-16 2005-07-28 Fuji Photo Film Co Ltd 異常陰影判定方法、異常陰影判定装置およびそのプログラム
JP2005253708A (ja) * 2004-03-12 2005-09-22 Fuji Photo Film Co Ltd 特徴量選択装置、異常陰影判別装置およびプログラム
JP4647360B2 (ja) * 2004-04-05 2011-03-09 富士フイルム株式会社 差分画像作成装置、差分画像作成方法、及び、そのプログラム
JP2006130212A (ja) * 2004-11-09 2006-05-25 Fuji Photo Film Co Ltd 異常陰影候補検出方法、および装置、ならびにプログラム

Also Published As

Publication number Publication date
WO2008020648A1 (en) 2008-02-21
US20080044070A1 (en) 2008-02-21
JP2010500049A (ja) 2010-01-07

Similar Documents

Publication Publication Date Title
CN105957063B (zh) 基于多尺度加权相似性测度的ct图像肝脏分割方法及系统
EP3128919B1 (en) Lung segmentation and bone suppression techniques for radiographic images
Figueiredo et al. Automated lesion detectors in retinal fundus images
Kulkarni et al. Classification of lung cancer stages on CT scan images using image processing
JP6273291B2 (ja) 画像処理装置および方法
Bandyopadhyay Pre-processing of mammogram images
JP2004135868A (ja) 異常陰影候補検出処理システム
Hogeweg et al. Suppression of translucent elongated structures: applications in chest radiography
JP5027872B2 (ja) 異常塊候補を検出する装置、方法、ならびにコンピュータ読取り可能媒体
Sethi et al. Segmentation of cancerous regions in liver using an edge-based and phase congruent region enhancement method
Yang et al. An automated method for accurate vessel segmentation
EP3258834A1 (en) Chest radiograph (cxr) image analysis
EP3472805B1 (en) A method and apparatus for mapping at least part of a structure in an image of at least part of a body of a subject
Xu et al. An automatic tongue detection and segmentation framework for computer-aided tongue image analysis
JP2006346094A (ja) 検出情報の出力方法及び医用画像処理システム
US10579895B2 (en) Information processing apparatus, extraction method, and medium
JP4124406B2 (ja) 異常陰影検出装置
Bhatia et al. Proposed algorithm to blotch grey matter from tumored and non tumored brain MRI images
JP2001299740A (ja) 異常陰影検出処理システム
Lim et al. A robust segmentation framework for spine trauma diagnosis
Kumar et al. Neural Network Based Approach for Detection of Abnormal Regions of Lung Cancer in X-Ray Image
JP5432483B2 (ja) 異常陰影検出装置およびプログラム
JP2005246032A (ja) 異常陰影検出方法および装置並びにプログラム
JP2004187998A (ja) 異常陰影形状抽出装置
Apostolopoulos et al. Computer aided classification of mammographic tissue using shapelets and support vector machines

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20100219

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20120524

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20120622

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20150629

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees