(57)【特許請求の範囲】
【請求項１】
原画像を表すオリジナル画像信号Ｓorg 、該オリジナル画像信号Ｓorg の非鮮鋭マスク信号Ｓus、強調係数βに基づいて、下記式（１）にしたがった強調処理を施して処理済画像信号Ｓprocを得る画像処理方法において、
Ｓproc＝Ｓorg ＋β×（Ｓorg −Ｓus）              （１）
前記原画像から該原画像中の所定の大きさの輪郭を有する画像部分に対応する大きさの輪郭を有する構造要素を用いて該画像部分を抽出するモーフォロジー演算を前記オリジナル画像信号に施して、該輪郭を有する画像部分において特徴的な値を示す特定画像信号を求め、
該特定画像信号が大きいほど前記強調度合いが大きくなるような特定画像強調係数を求め、
前記原画像のうちエッジ部分において特徴的な値を示すエッジ信号を求め、
前記特定画像強調係数を前記エッジ信号が大きいほど調整後の特定画像強調係数が大きくなるように調整し、
該調整された新たな特定画像強調係数を式（１）における強調係数βとすることを特徴とする画像処理方法。
【請求項２】
前記モーフォロジー演算がオープニング処理またはクロージング処理に基づくものであることを特徴とする請求項１記載の画像処理方法。
【請求項３】
前記特定画像強調係数は、該特定画像信号の大きさに応じて予め設定された変換テーブルにしたがって求めるものであることを特徴とする請求項１または２記載の画像処理方法。
【請求項４】
前記特定画像強調係数を求めるのに先だって、ゼロ以外の値を示す特定画像信号が存在する範囲を広げる膨張処理を前記特定画像信号に施すことを特徴とする請求項１から３のうちいずれか１項に記載の画像処理方法。
【請求項５】
前記エッジ信号は、前記原画像のうち前記エッジ部分を抽出する第２のモーフォロジー演算に基づいて求めることを特徴とする請求項１から４のうちいずれか１項に記載の画像処理方法。
【請求項６】
前記第２のモーフォロジー演算は、ダイレーション処理またはエロージョン処理に基づくものであることを特徴とする請求項５記載の画像処理方法。
【請求項７】
前記エッジ信号の大きさに応じたエッジ強調係数を求め、
前記エッジ強調係数に基づいて前記特定画像強調係数を調整することを特徴とする請求項１から６のうちいずれか１項に記載の画像処理方法。
【請求項８】
前記エッジ強調係数は、前記エッジ信号の大きさに応じて予め設定された第２の変換テーブルにしたがって求めるものであることを特徴とする請求項７記載の画像処理方法。
【請求項９】
原画像を表すオリジナル画像信号Ｓorg 、該オリジナル画像信号Ｓorg の非鮮鋭マスク信号Ｓus、強調係数βに基づいて、下記式（１）にしたがった強調処理を施して処理済画像信号Ｓprocを得る画像処理装置において、
Ｓproc＝Ｓorg ＋β×（Ｓorg −Ｓus）              （１）
前記原画像から該原画像中の所定の大きさの輪郭を有する画像部分に対応する大きさの輪郭を有する構造要素を用いて該画像部分を抽出するモーフォロジー演算を前記オリジナル画像信号に施して、該輪郭を有する画像部分において特徴的な値を示す特定画像信号を求めるモーフォロジー演算手段と、
該特定画像信号が大きいほど前記強調度合いが大きくなるように特定画像強調係数が予め設定された変換テーブルと、
前記原画像のうちエッジ部分において特徴的な値を示すエッジ信号を求めるエッジ信号検出手段と、
前記特定画像強調係数を前記エッジ信号が大きいほど調整後の特定画像強調係数が大きくなるように調整する調整手段と、
該調整された新たな特定画像強調係数を前記強調係数βとして前記式（１）の強調処理を行なう強調処理手段とを備えたことを特徴とする画像処理装置。
【請求項１０】
前記モーフォロジー演算手段によるモーフォロジー演算がオープニング処理またはクロージング処理に基づくものであることを特徴とする請求項９記載の画像処理装置。
【請求項１１】
前記特定画像信号に対して、ゼロ以外の値を示す特定画像信号が存在する範囲を広げるように膨張処理を施す膨張処理手段をさらに備えたことを特徴とする請求項９または１０記載の画像処理装置。
【請求項１２】
前記エッジ信号検出手段が、前記原画像のうち前記エッジ部分を抽出する第２のモーフォロジー演算に基づいて求める第２のモーフォロジー演算手段であることを特徴とする請求項９から１１のうちいずれか１項に記載の画像処理装置。
【請求項１３】
前記第２のモーフォロジー演算手段によるモーフォロジー演算が、ダイレーション処理またはエロージョン処理に基づくものであることを特徴とする請求項１２記載の画像処理装置
【請求項１４】
前記エッジ信号の大きさに応じたエッジ強調係数が予め設定された第２の変換テーブルをさらに備え、
前記調整手段が、前記特定画像強調係数を前記エッジ強調係数に基づいて調整するものであることを特徴とする請求項９から１３のうちいずれか１項に記載の画像処理装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は画像処理方法および画像処理装置に関し、詳細には画像のうち、異常陰影等の特定の画像部分だけ抽出して強調処理する方法および装置の改良に関するものである。
【０００２】
【従来の技術】
従来より、種々の画像取得方法により得られた原画像を表す画像信号（オリジナル画像信号という）に対して、階調処理や周波数処理等の画像処理を施し、画像の観察読影性能を向上させることが行われている。特に人体を被写体とした放射線画像のような医用画像の分野においては、医師等の専門家が、得られた画像に基づいて患者の疾病や傷害の有無を的確に診断する必要があり、その画像の読影性能を向上させる画像処理は不可欠なものとなっている。
【０００３】
この画像処理のうち、いわゆる周波数強調処理としては、例えば特開昭61-169971 号に示されるように、原画像の濃度値等の画像信号Ｓorg を、
Ｓproc＝Ｓorg ＋β×（Ｓorg −Ｓus）              （１）
なる画像信号Ｓprocに変換するものが知られている。
【０００４】
ここでβは周波数強調係数、Ｓusは非鮮鋭マスク（いわゆるボケマスク）信号である。このボケマスク信号Ｓusは、２次元に配置された画素に対してオリジナル画像信号Ｓorg を中心画素とするＮ列×Ｎ行（Ｎは奇数）の画素マトリクスからなるマスク、すなわちボケマスクを設定し、
Ｓus＝（ΣΣＳorg ）／Ｎ2                         （２）
（ただし、ΣΣＳorg はボケマスク内の画素の画像信号の和）
等として求められる超低空間周波数成分である。
【０００５】
式（１）の第２項括弧内の値（Ｓorg −Ｓus）は、オリジナル画像信号Ｓorg から超低空間周波数成分であるボケマスク信号Ｓusを減算したものであるから、オリジナル画像信号Ｓorg のうちの、超低空間周波数成分を除去した比較的高い周波数成分を選択的に抽出したものを意味する。
【０００６】
この比較的高い周波数成分に周波数強調係数βを乗じたうえで、オリジナル画像信号を加算することにより、原画像のうちこの比較的高い周波数成分だけを選択的に、かつ相対的に強調・抑制（鮮鋭度の調整）を行なうことができる。
【０００７】
一方、原画像のうち異常な陰影等の特定の画像部分やエッジ部分だけを選択的に抽出する、モーフォロジー（Morphology；モフォロジーまたはモルフォロジーとも称する）のアルゴリズムに基づく処理（以下、モーフォロジー演算またはモーフォロジー処理という）が知られている。
【０００８】
このモーフォロジー処理は、特に乳癌における特徴的形態である微小石灰化像を検出するのに有効な手法として研究されているが、対象画像としてはこのようなマンモグラムにおける微小石灰化像に限るものではなく、検出しようとする特定の画像部分（異常陰影等）の大きさや形状が予めある程度分かっているものについては、いかなる画像に対しても適用することができる。
【０００９】
そしてこのモーフォロジー処理は、抽出しようとする画像部分の大きさに応じて設定される構造要素（マスク）Ｂおよびマルチスケールλとを用いて行なう処理であり、
（１）石灰化像そのものの抽出に有効であること、
（２）複雑なバックグラウンド情報に影響されにくいこと、
（３）抽出した石灰化像がひずまないこと、
などの特徴がある。
【００１０】
すなわち、この手法は一般の微分処理に比べて、石灰化像のサイズ・形状・濃度分布などの幾何学的情報をよりよく保って検出することができる。以下、このモーフォロジー処理の概要を、マンモグラムにおける微小石灰化像の検出に適用した例について説明する。
【００１１】
（モーフォロジーの基本演算）
モーフォロジー処理は一般的にはＮ次元空間における集合論として展開されるが、直感的な理解のために２次元の濃淡画像を対象として説明する。
【００１２】
濃淡画像を座標（ｘ，ｙ）の点が濃度値ｆ（ｘ，ｙ）に相当する高さをもつ空間とみなす。ここで、濃度値ｆ（ｘ，ｙ）は、濃度が低い（ＣＲＴに表示した場合には輝度が高い）程大きな画像信号値となる高輝度高信号レベルの信号とする。
【００１３】
まず、簡単のため上記２次元の濃淡画像の断面に相当する１次元の関数ｆ（ｘ）を考える。モーフォロジー演算に用いる構造要素ｇ（前述のＢに相当）は次式（３）に示すように、原点について対称な対称関数
【００１４】
【数１】
0003669455.tif 000002
【００１５】
であり、定義域内で値が０で、その定義域Ｇが下記式（４）であるとする。
【００１６】
【数２】
0003669455.tif 000003
【００１７】
このとき、モーフォロジー演算の基本形は式（５）〜（８）に示すように、非常に簡単な演算となる。
【００１８】
【数３】
0003669455.tif 000004
【００１９】
すなわち、ダイレーション（dilation）処理は、注目画素を中心とした、±ｍ（構造要素Ｂに応じて決定される値であって、図８中のマスクサイズに相当）の幅の範囲内の最大値を探索する処理であり（図８（Ａ）参照）、一方、エロージョン（erosion ）処理は、注目画素を中心とした、±ｍの幅の範囲内の最小値を探索する処理である（図８（Ｂ）参照）。このことからダイレーション処理を最大値処理、エロージョン処理を最小値処理とも称するものとする。
【００２０】
また、オープニング（opening ）処理はエロージョン処理後にダイレーション処理を行なう処理、すなわち最小値の探索の後に最大値を探索する処理であり、クロージング（closing ）処理は、ダイレーション処理後にエロージョン処理を行なう処理、すなわち最大値の探索の後に最小値を探索する処理に相当する。
【００２１】
つまりオープニング処理は、低輝度側から濃度曲線ｆ（ｘ）を滑らかにし、マスクサイズ２ｍより空間的に狭い範囲で変動する凸状の濃度変動部分（周囲部分よりも輝度が高い部分）を取り除くことに相当する（図８（Ｃ）参照）。一方、クロージング処理は、高輝度側から濃度曲線ｆ（ｘ）を滑らかにし、マスクサイズ２ｍより空間的に狭い範囲で変動する凹状の濃度変動部分（周囲部分よりも輝度が低い部分）を取り除くことに相当する（図８（Ｄ）参照）。
【００２２】
なお、構造要素ｇが原点に対して対称ではない場合の、式（５）に示すダイレーション演算をミンコフスキー（Minkowski ）和、式（６）に示すエロージョン演算をミンコフスキー差という。
【００２３】
ここで、濃度の高いもの程大きな値となる高濃度高信号レベルの信号の場合においては、濃度値ｆ（ｘ）の画像信号値が高輝度高信号レベルの場合に対して大小関係が逆転するため、高濃度高信号レベルの信号に対するダイレーション処理と高輝度高信号レベルに対するエロージョン処理（図８（Ｂ））とは一致し、高濃度高信号レベルの信号に対するエロージョン処理と高輝度高信号レベルに対するダイレーション処理（図８（Ａ））とは一致し、高濃度高信号レベルの信号に対するオープニング処理と高輝度高信号レベルに対するクロージング処理（図８（Ｄ））とは一致し、高濃度高信号レベルの信号に対するクロージング処理と高輝度高信号レベルに対するオープニング処理（図８（Ｃ））とは一致する。
【００２４】
なお、本項では高輝度高信号レベルの画像信号（輝度値）の場合について説明する。
【００２５】
（石灰化陰影検出への応用）
石灰化陰影の検出には、原画像から平滑化した画像を引き去る差分法が考えられる。単純な平滑化法では石灰化陰影と細長い形状の非石灰化陰影（乳腺、血管および乳腺支持組織等）との識別が困難であるため、東京農工大の小畑らは、多重構造要素を用いたオープニング演算に基づく下記式（９）で表されるモーフォロジーフィルターを提案している（「多重構造要素を用いたモルフォロジーフィルタによる微小石灰化像の抽出」電子情報通信学会論文誌 D-II Vol.J75-D-II No.7 P1170 〜1176  1992年７月、「モルフォロジーの基礎とそのマンモグラム処理への応用」MEDICAL IMAGING TECHNOLOGY Vol.12 No.1 January 1994 等）。
【００２６】
【数４】
0003669455.tif 000005
【００２７】
ここでＢi （ｉ＝１，２，…，Ｍ）は、直線状のＭ個（例えば図９に示すものではＭ＝４）の構造要素（これらＭ個の構造要素の全体をもって多重構造要素という）である。構造要素Ｂi を検出対象である石灰化陰影よりも大きく設定すれば、上記オープニング演算による処理で、構造要素Ｂi よりも細かな信号変化部分（空間的に狭い範囲で変動する画像部分）であって周囲よりも輝度値の大きい凸状の部分である石灰化陰影は取り除かれる。一方、細長い形状の非石灰化陰影はその長さが構造要素Ｂi よりも長く、その傾きが４つの構造要素Ｂi のいずれかに一致すればオープニング処理（式（９）の第２項の演算）をしてもそのまま残る。したがってオープニング処理によって得られた平滑化画像（石灰化陰影のみが取り除かれた画像）を原画像ｆから引き去ることで、小さな石灰化陰影のみが含まれる画像が得られる。これが式（９）の考え方である。
【００２８】
なお、前述したように、高濃度高信号レベルの信号の場合においては、石灰化陰影は周囲の画像部分よりも濃度値が低くなり、石灰化陰影は周囲部分に対して濃度値の小さい凹状の信号変化部分となるため、オープニング処理に代えてクロージング処理を適用し、式（９）に代えて式（10）を適用する。
【００２９】
【数５】
0003669455.tif 000006
【００３０】
一方、ダイレーション処理やエロージョン処理によって得られた画像を原画像から引き去ること（式（11），（12）参照）で、原画像のエッジ部（図８（Ａ）または（Ｂ）の斜線部分）だけを選択的に抽出することができる。
【００３１】
【数６】
0003669455.tif 000007
【００３２】
【発明が解決しようとする課題】
ところで上述したように、画像の読影性能を向上させるには、対象となる画像に対して画像処理を行うことが不可欠となっているが、特開平２−1078号に開示されているように、単に濃度依存による強調処理では、例えばマンモグラムにおける放射線ノイズ成分のような、画像読影の障害となる成分まで強調されるため、読影性能をむしろ低下させることになる。
【００３３】
また特開昭60−192482号、特開平２−120985号、特表平３−502975号等に開示されているように、画像信号の分散値に依存した強調処理では、局所的に濃度変化の大きい画像部分が強く強調されるため、その付近でアンダーシュート、オーバーシュートが相対的に目立ち、特にＸ線画像に関しては高濃度側でアーチファクトが発生しやすいという問題がある。
【００３４】
そこで本願出願人は、上記各モーフォロジー信号に依存した強調を行なうことにより、ノイズ成分等の画像読影に不要な成分を強調することなく、注目する特定の画像部分だけを効率よく強調処理する方法を既に提案している（特願平７-71774号、同７-204250 号等参照）。
【００３５】
この方法によれば確かに、その大きさが構造要素よりも小さい石灰化陰影だけを注目画像部分として選択的に抽出することができるため、この抽出した信号に依存した強調処理を行なうことにより、石灰化陰影だけを効率よく強調処理することができる、という従来にはない顕著な効果を得ることができる。
【００３６】
しかし、本願出願人の提案した上記方法によっても希に不自然な強調画像が発生することがある。
【００３７】
すなわち、例えば図10に示すような各オリジナル画像信号に対してモーフォロジー演算の一つであるクロージング処理を施すと、構造要素より小さな石灰化陰影▲１▼については完全に平滑化され、構造要素より大きな陰影のうち信号値の変動が殆どない平坦部分の大きさが構造要素より十分大きな陰影▲３▼についてはオリジナル画像信号を維持するため、この得られた信号をオリジナル画像信号から引き去ることによって、構造要素より小さな石灰化陰影▲１▼については０以外の信号値を発生し、構造要素より大きな陰影のうち信号値の変動が殆どない平坦部分の大きさが構造要素より十分大きな陰影▲３▼については信号値が０となって、両者を分離することが可能であるが、構造要素より大きな陰影のうち信号値の変動が殆どない平坦部分の大きさが構造要素より小さい陰影▲２▼については、その構造要素より小さい範囲で信号値が変動する部分（陰影の中心部に相当する部分）が、上記構造要素より小さな石灰化陰影▲１▼と同様に平滑化されるため、オリジナル画像信号から引き去ると、この陰影▲２▼のうち、中心部を除いた裾野に相当する部分では信号値が０となり、一方、この裾野よりも陰影内部である中心部に相当する部分は０以外の信号を発生する。この結果、上記陰影▲２▼についてはその中心部を含む、構造要素より小さい範囲だけが強調されることになる。
【００３８】
本発明は上記事情に鑑みなされたものであって、石灰化陰影等の特定の画像部分だけを選択的に強調処理する前提となる、他の画像部分からの当該特定部分の分離精度を高めて不自然な強調がなさるれるのを防止した画像処理方法および画像処理装置を提供することを目的とするものである。
【００３９】
【課題を解決するための手段】
本発明の画像処理方法は、原画像を表すオリジナル画像信号Ｓorg 、該オリジナル画像信号Ｓorg の非鮮鋭マスク信号Ｓus、強調係数βに基づいて、下記式（１）にしたがった強調処理を施して処理済画像信号Ｓprocを得る画像処理方法において、
Ｓproc＝Ｓorg ＋β×（Ｓorg −Ｓus）              （１）
前記原画像から該原画像中の所定の大きさの輪郭を有する画像部分に対応する大きさの輪郭を有する構造要素を用いて該画像部分を抽出するモーフォロジー演算を前記オリジナル画像信号に施して、該輪郭を有する画像部分において特徴的な値を示す特定画像信号を求め、
該特定画像信号が大きいほど前記強調度合いが大きくなるような特定画像強調係数を求め、
前記原画像のうちエッジ部分において特徴的な値を示すエッジ信号を求め、
前記特定画像強調係数を前記エッジ信号が大きいほど調整後の特定画像強調係数が大きくなるように調整し、
該調整された新たな特定画像強調係数を式（１）における強調係数βとすることを特徴とするものである。
【００４０】
ここで、輪郭を有する画像部分とは、例えば石灰化陰影等の異常陰影などを意味し、濃度変化が急峻な単なるエッジ部分ではなく、そのようなエッジ部分に囲まれた、閉じた画像部分を意味する。
なお、「特定画像信号が大きいほど前記強調度合いが大きくなるような特定画像強調係数」には、特定画像信号が大きくなっても一部に強調度合いが変化しない部分を含むものであって、少なくとも特定画像信号が大きくなるほど強調度合いが小さくなることはないものであればよい。
また、「前記特定画像強調係数を前記エッジ信号が大きいほど調整後の特定画像強調係数が大きくなるように調整」には、エッジ信号が大きくなっても一部に特定画像強調係数が変化しない部分を含むものであって、少なくともエッジ信号が大きくなるほど特定画像強調係数が小さくなることはないものであればよい。
【００４１】
ここで「モーフォロジー演算」とは、画像中から抽出したい画像部分に対応する大きさの輪郭を有する構造要素を用いてオープニング処理またはクロージング処理を画像信号に施すことにより、そのような画像部分を抽出する処理である（例えば、電子情報通信学会誌第７４巻第２号１６６頁から１７３頁、同第７４巻第３号２７１頁から２７９頁参照））。
なお、上記モフォロジー演算は、一般的にはＮ次元空間における集合論として展開されるが、２次元空間である画像へ適用される場合が多く、ダイレーション（dilation）処理、イロージョン（erosion ）処理、オープニング（opening ）処理およびクロージング（closing ）処理を単独で、またはこれらの処理のうち２以上を組み合わせて適用される演算処理である。
例えば、オリジナル画像信号Ｓorg に対して、強調処理しようとする画像部分の大きさに対応して予め設定した多重構造要素Ｂi （または単一の構造要素Ｂ）を用いて式（９）の第２項（多重構造要素の場合）若しくは式（７）（単一の構造要素の場合）に示すオープニング処理を施し、または式（10）の第２項（多重構造要素の場合）若しくは式（８）（単一の構造要素の場合）に示すクロージング処理を施し、オリジナル画像信号Ｓorg からこのオープニング処理したものを引き去る演算（式（13））、またはクロージング処理したものからオリジナル画像信号Ｓorgを引き去る演算（式（14））を行えばよい。
【００４２】
Ｓ′＝Ｓorg −Ｓopening                         （13）
Ｓ′＝Ｓclosing −Ｓorg                         （14）
ただし、Ｓopening はオリジナル画像信号Ｓorg に対するオープニング処理結果、Ｓclosing はオリジナル画像信号Ｓorg に対するクロージング処理結果をそれぞれ示し、Ｓ′は特定画像信号を示す。
【００４３】
なお、オリジナル画像信号Ｓorg が高濃度高信号レベルの信号であり、強調対象の画像部分が石灰化陰影である場合には、石灰化陰影は周囲よりも信号値が小さいため、最大値の探索の後に最小値を探索する処理であるクロージング処理を行なうことが必要である。
【００４４】
一方、オリジナル画像信号Ｓorg が高輝度高信号レベルの信号であり、強調対象の画像部分が石灰化陰影である場合には、石灰化陰影は周囲よりも信号値が大きいため、最小値の探索の後に最大値を探索する処理であるオープニング処理を行なうことが必要である。
【００４５】
また特定画像強調係数とは特定画像信号に応じた係数、すなわちオリジナル画像信号Ｓorg のモーフォロジー演算結果に依存した係数を意味する。そしてこの特定画像強調係数は、例えば図２のグラフに示すような特定画像信号の大きさに応じて予め設定された変換テーブルにしたがって設定するのが望ましく、その大きさを０から１の範囲とするのが後の演算処理を行なううえで便利である。
【００４６】
なお、特定画像強調係数を求めるのに先だって、ゼロ以外の値（すなわち、特徴的な値）を示す特定画像信号が存在する範囲を広げる膨張処理を特定画像信号に施すのが望ましい。
【００４７】
この膨張処理としては、ダイレーション処理、エロージョン処理の他、非鮮鋭マスク処理（ボケマスク処理）などを適用することができる。ただし、オリジナル画像信号Ｓorg が高濃度高信号レベルの信号であり、強調対象の画像部分が石灰化陰影である場合には、最小値を探索する処理であるエロージョン処理、または非鮮鋭マスク処理を適用することが必要であり、一方、オリジナル画像信号Ｓorg が高輝度高信号レベルの信号であり、強調対象の画像部分が石灰化陰影である場合には、石灰化陰影は周囲よりも信号値が大きいため、最大値を探索する処理であるダイレーション処理、または非鮮鋭マスク処理を適用することが必要である。
【００４８】
さらにこの膨張処理による膨張の程度としては、特定画像信号がゼロ以外の値（特徴的な信号）を示す画像部分を、少なくとも特定画像部分本来の大きさ（広がり）よりも大きくするように膨張するのが適切である。
【００４９】
またエッジ信号は、原画像のうち上記エッジ部分を抽出する第２のモーフォロジー演算、特に式（５），（６）と同様にダイレーション処理またはエロージョン処理を施し、このダイレーション処理したものからオリジナル画像信号Ｓorg を引き去る演算（式（15））またはオリジナル画像信号Ｓorg からエロージョン処理したものを引き去る演算（式（16））を行なっても求めることができる。ただし、必ずしもモーフォロジー演算によって求めることに限るものではなく、公知の各種のエッジ検出方法を用いることも可能である。
【００５０】
Ｓedge＝Ｓdilation−Ｓorg                       （15）
Ｓedge＝Ｓorg −Ｓerosion                       （16）
ここで、Ｓdilationはオリジナル画像信号Ｓorg に対するダイレーション処理結果、Ｓerosion はオリジナル画像信号Ｓorg に対するエロージョン処理結果をそれぞれ示し、Ｓedgeはエッジ信号を示す。
【００５１】
さらに、このエッジ信号についてもエッジ信号の大きさに応じたエッジ強調係数を求め、このエッジ強調係数に基づいて特定画像強調係数を調整するのが望ましい。
【００５２】
ここでエッジ強調係数とはエッジ信号に応じた係数、すなわち原画像のエッジ部を示すエッジ信号に依存した係数を意味する。そしてこのエッジ強調係数は、例えば図３のグラフに示すようなエッジ信号の大きさに応じて予め設定された変換テーブルにしたがって設定するのが望ましく、その大きさを０から１の範囲とするのが後の演算処理を行なううえで便利である。
【００５３】
さらにエッジ強調係数に基づいて特定画像強調係数を調整する場合は、特定画像強調係数にエッジ強調係数を乗じるのがよい。特に、エッジ強調係数が０から１の範囲の値とされている場合であって、１をエッジ部を表す値、０をエッジ部以外の部分を表す値として設定したときは、特定画像強調係数が特定画像を表す値を示したときも、エッジ強調係数がエッジ部以外の部分を表す値０を採れば、両係数の積が０となり、特定画像強調係数が表した偽画像の強調を防止することができる。
【００５４】
これは、逆の立場から見れば、特定画像強調係数が０から１の範囲の値とされている場合であって、１を石灰化陰影を表す値、０をそれ以外の部分を表す値として設定したときは、エッジ強調係数がエッジ部を表す値を示したときも、特定画像強調係数が石灰化陰影以外の部分を表す値０を採れば、両係数の乗算の結果が０となり、エッジ強調係数が表したエッジ部のうち乳腺等の石灰化陰影以外の陰影によるエッジ部の強調を防止することができる。
【００５５】
本発明の画像処理装置は、原画像を表すオリジナル画像信号Ｓorg 、該オリジナル画像信号Ｓorg の非鮮鋭マスク信号Ｓus、強調係数βに基づいて、下記式（１）にしたがった強調処理を施して処理済画像信号Ｓprocを得る画像処理装置において、
Ｓproc＝Ｓorg ＋β×（Ｓorg −Ｓus）              （１）
前記原画像から該原画像中の所定の大きさの輪郭を有する画像部分に対応する大きさの輪郭を有する構造要素を用いて該画像部分を抽出するモーフォロジー演算を前記オリジナル画像信号に施して、該輪郭を有する画像部分において特徴的な値を示す特定画像信号を求めるモーフォロジー演算手段と、
該特定画像信号が大きいほど前記強調度合いが大きくなるように特定画像強調係数が予め設定された変換テーブルと、
前記原画像のうちエッジ部分において特徴的な値を示すエッジ信号を求めるエッジ信号検出手段と、
前記特定画像強調係数を前記エッジ信号が大きいほど調整後の特定画像強調係数が大きくなるように調整する調整手段と、
該調整された新たな特定画像強調係数を前記強調係数βとして前記式（１）の強調処理を行なう強調処理手段とを備えたことを特徴とするものである。
【００５６】
上記モーフォロジー演算としては、前記本発明の画像処理方法と同様に、オープニング処理またはクロージング処理に基づくものが望ましい。すなわちオリジナル画像信号Ｓorg に対して、予め設定した多重構造要素Ｂi （または単一の構造要素Ｂ）を用いて式（９）の第２項（多重構造要素の場合）若しくは式（７）（単一の構造要素の場合）に示すオープニング処理を施し、または式（10）の第２項（多重構造要素の場合）若しくは式（８）（単一の構造要素の場合）に示すクロージング処理を施し、オリジナル画像信号Ｓorg からこのオープニング処理したものを引き去る演算（式（13））、またはクロージング処理したものからオリジナル画像信号Ｓorg を引き去る演算（式（14））を行えばよい。
【００５７】
なお、オリジナル画像信号Ｓorg が高濃度高信号レベルの信号であり、強調対象の画像部分が石灰化陰影である場合には、石灰化陰影は周囲よりも信号値が小さいため、最大値の探索の後に最小値を探索する処理であるクロージング処理を行なうことが必要であり、一方、オリジナル画像信号Ｓorg が高輝度高信号レベルの信号であり、強調対象の画像部分が石灰化陰影である場合には、石灰化陰影は周囲よりも信号値が大きいため、最小値の探索の後に最大値を探索する処理であるオープニング処理を行なうことが必要である。
【００５８】
また変換テーブルとしては、例えば図２のグラフに示すような特定画像信号の大きさに応じて予め設定されたものが望ましく、特定画像強調係数を０から１の範囲とするのが後の演算処理を行なううえで便利である。
【００５９】
さらに、特定画像信号に対して、ゼロ以外の値を示す特定画像信号が存在する範囲を広げるように膨張処理を施す膨張処理手段をさらに備えるのが後の処理のために望ましい。
【００６０】
この膨張処理としては、ダイレーション処理、エロージョン処理の他、非鮮鋭マスク処理（ボケマスク処理）などを適用することができる。ただし、オリジナル画像信号Ｓorg が高濃度高信号レベルの信号であり、強調対象の画像部分が石灰化陰影である場合には、最小値を探索する処理であるエロージョン処理、または非鮮鋭マスク処理を適用することが必要であり、一方、オリジナル画像信号Ｓorg が高輝度高信号レベルの信号であり、強調対象の画像部分が石灰化陰影である場合には、石灰化陰影は周囲よりも信号値が大きいため、最大値を探索する処理であるダイレーション処理、または非鮮鋭マスク処理を適用することが必要である。
【００６１】
またエッジ信号検出手段は、原画像のうち上記エッジ部分を抽出する第２のモーフォロジー演算、特に式（５），（６）と同様にダイレーション処理またはエロージョン処理を施し、このダイレーション処理したものからオリジナル画像信号Ｓorg を引き去る演算（式（15））またはオリジナル画像信号Ｓorg からエロージョン処理したものを引き去る演算（式（16））を行なっても求めることができる。ただし、必ずしもモーフォロジー演算によって求めることに限るものではなく、公知の各種のエッジ信号検出手段を用いることも可能である。
【００６２】
さらに、このエッジ信号についてもエッジ信号の大きさに応じたエッジ強調係数を求め、このエッジ強調係数に基づいて調整手段が特定画像強調係数を調整するのが望ましい。
【００６３】
ここでエッジ強調係数は、例えば図３のグラフに示すようなエッジ信号の大きさに応じて予め設定された変換テーブルにしたがって設定されるのが望ましく、その大きさを０から１の範囲とするのが後の演算処理を行なううえで便利である。
【００６４】
なお、上記各発明においてモーフォロジー演算を行なう場合に用いられる単一の構造要素Ｂとしては、例えば、正方形、長方形、円形、楕円形、または菱形等の上下左右対称の要素が望ましい。これは画像信号は実際には２次元状に分布するものだからである。
【００６５】
【発明の効果】
本発明の画像処理方法および装置は、原画像を表すオリジナルの画像信号Ｓorg に対してモーフォロジー演算を施すことにより、原画像から所定の大きさの輪郭を有する画像部分、例えば乳がんの特徴的構造である石灰化陰影等を検出し、一方、原画像から濃度変化が急峻なエッジ部分を検出し、さらに石灰化陰影等の画像部分を表す指標値としての特定画像強調係数を、エッジ部分を表す指標値たるエッジ信号に基づいて調整することにより、特定画像強調係数が特定画像を表す値を示したときも、エッジ強調係数がエッジ部以外の部分を表す値を採れば、両係数の関係から、特定画像強調係数が検出した偽画像についての強調を防止することができる。
【００６６】
逆に、エッジ信号がエッジ部を表す値を示したときも、特定画像強調係数が石灰化陰影以外の部分を示す値を採れば、両係数の関係から、エッジ信号が示すエッジ部のうち乳腺等の石灰化陰影以外の陰影によるエッジ部の強調を防止することができる。
【００６７】
以上により、石灰化陰影等の特定の画像部分だけを選択的に強調処理する前提となる、他の画像部分からの当該特定部分の分離精度を高めることができる。
【００６８】
【発明の実施の形態】
以下、本発明の画像処理方法を用いた画像処理装置について図面を用いて説明する。
【００６９】
図１は本実施形態の画像処理装置を示す概略ブロック図である。図示の画像処理装置は、Ｘ線画像に対して画像処理を施して、所定の輪郭を有し濃度値が周囲の画像部分よりも低い微小な石灰化陰影を選択的に強調処理する画像処理装置である。
【００７０】
この画像処理装置は、濃度信号（高濃度高信号レベル）として画像を表すオリジナル画像信号Ｓorg の、超低空間周波数に対応する非鮮鋭マスク信号Ｓusを求めるローパスフィルタ11と、オリジナル画像信号Ｓorg から非鮮鋭マスク信号Ｓusを減算して比較的高周波の成分（Ｓorg −Ｓus）を抽出する減算器17と、オリジナル画像信号Ｓorg から石灰化陰影を抽出しこの石灰化陰影を表す石灰化陰影信号に依存した第１の強調係数αcalcを求める石灰化陰影依存強調係数算出手段20と、オリジナル画像信号Ｓorg からエッジ部分を抽出しこのエッジ部分を表すエッジ信号に依存した第２の強調係数αedgeを求めるエッジ信号依存強調係数算出手段30と、これら２つの強調係数αcalc、αedgeを互いに乗じて１つの強調係数αを算出する乗算器41と、オリジナル画像信号Ｓorg 、高周波の成分（Ｓorg −Ｓus）および強調係数αを用いて下記式（17）にしたがった信号処理を施して処理済み画像信号Ｓprocを得る乗算器42および加算器43とを備えた構成である。
【００７１】
Ｓproc＝Ｓorg ＋α×（Ｓorg −Ｓus）              （17）
ここで石灰化陰影依存強調係数算出手段20はさらに詳しくは、オリジナル画像信号Ｓorg に対して、強調処理しようとする対象の石灰化陰影等の大きさに対応した多重構造要素（multi-element ）Ｂi およびスケール係数λを用いてクロージング処理（式（10）の第２項参照）を施す第１のモーフォロジー演算手段12と、このクロージング処理を施して得られた信号Ｓclosing からオリジナル画像信号Ｓorg を減算（式（14）参照）して石灰化陰影信号Ｓ′を算出する減算器18と、この減算処理の結果に対してダイレーション処理を施し石灰化陰影信号Ｓ′が０以外の値を採る画素で囲まれる画像部分を空間的に膨張させる膨張処理手段14と、この膨張処理の結果である膨張処理信号Ｓcalcと第１の強調係数αcalcとが予め対応付けられた第１のルックアップテーブル15（図２参照）とを備えた構成である。
【００７２】
なお、構造要素としては例えば、線形の３画素からなる単一の構造要素Ｂが互いに異なる４方向（45度間隔）に延びるものとして設定された多重構造要素Ｂi （図11（Ａ）参照）を用い、スケール係数λとしては１を用いるものとし、一方、膨張処理手段14によるダイレーション処理に用いられる構造要素としては例えば、縦９画素×横９画素の正方形形状の単一の構造要素Ｂを用い、スケール係数λとしては１を用いるものとする。
【００７３】
なお、多重構造要素Ｂi を構成する各構造要素Ｂの画素数としては３画素に限るものではなく、図11（Ｂ）に示すような５画素（８方向）、同図（Ｃ）に示すような７画素（８方向）、…等のように、強調しようとする石灰化陰影の大きさに応じて設定することができる。また多重構造要素Ｂi を構成する各構造要素Ｂの延びる方向数も４方向に限るものではなく、８方向や16方向等を適用することができる。ただし、各構造要素の画素数に応じて適切な方向数を設定すべきである。
【００７４】
すなわち３画素の構造要素であれば４方向が限度であり、８方向の設定は意味のないものとなるからである。したがって、３画素の場合は４方向、５画素の場合や７画素の場合は８方向とするのが妥当である。なお図11においては正方形格子状に配列された画素を対象としているため、斜め方向については厳密な意味での直線構造要素とはなっていない。
【００７５】
またエッジ信号依存強調係数算出手段30は、Ｘ線画像のエッジ部分において特徴的な信号を示すダイレーション処理を行なう第２のモーフォロジー演算手段13と、このダイレーション処理を施して得られた信号Ｓdilationからオリジナル画像信号Ｓorg を減算（式（15）参照）してエッジ信号Ｓedgeを算出する減算器19と、減算処理により得られたエッジ信号Ｓedgeに第２の強調係数αedgeが予め対応付けられた第２のルックアップテーブル16（図３参照）とを備えた構成である。
【００７６】
なおダイレーション処理に用いる構造要素としては例えば、縦５画素×横５画素の正方形形状の単一の構造要素Ｂを、スケール係数λとしては１を用いるものとする。
【００７７】
次に本実施形態の画像処理装置の作用について説明する。
【００７８】
まず前提として、ある患者の乳房を被写体としたＸ線画像が記録されたＸ線フイルムから、図示しない画像読取手段によってこのＸ線画像が光電的に読み取られてこの画像を表すオリジナル画像信号Ｓorg が生成される。
【００７９】
得られた画像信号は本実施形態の画像処理装置に入力される。
【００８０】
入力されたオリジナル画像信号Ｓorg は、ローパスフィルター11に入力され、ローパスフィルター11はこの画像信号Ｓorg に対して、例えば３列×３行の画素マトリクスからなるボケマスクを設定し、下記式（２）（Ｎ＝３に設定）にしたがって得られたボケマスク信号Ｓusを出力する。
【００８１】
Ｓus＝（ΣΣＳorg ）／Ｎ2                         （２）
（ただし、ΣΣＳorg はボケマスク内の画素の画像信号の和）
なお、ボケマスクとしては上記式（２）に示すようにマスク内の画素値の単純平均を用いるものの他、マスクの中心画素からの距離に応じてマスク内の画素値の重み付けを変化させたものを用いてもよい。
【００８２】
次いで演算素子17が、オリジナル画像信号Ｓorg からローパスフィルター11の出力であるボケマスク信号Ｓusを減算して、オリジナル画像信号Ｓorg うち高周波成分（Ｓorg −Ｓus）を出力する。
【００８３】
またこのローパスフィルタ11および演算素子17の作用と並行して、石灰化陰影依存強調係数算出手段20、エッジ信号依存強調係数算出手段30がそれぞれオリジナル画像信号Ｓorg に基づいて強調係数αcalc、αedgeを求める。
【００８４】
すなわち石灰化陰影依存強調係数算出手段20は、入力されたオリジナル画像信号Ｓorg に対して第１のモーフォロジー演算手段12が線形３画素４方向の多重構造要素Ｂi によりクロージング処理を施す。
【００８５】
オリジナル画像信号Ｓorg が図４右側最上段のグラフに示すように、乳腺構造と石灰化陰影を有する画像である場合であって、この石灰化陰影の大きさが３画素より大きい場合について適用すると、クロージング処理によって得られる画像信号Ｓclosing において、図４右側第２段目のグラフに示すように、多重構造要素Ｂi より大きい乳腺構造はオリジナル画像信号Ｓorg による形状を維持するが、石灰化陰影については、構造要素より空間的に狭い範囲の中心部は消去されるが多重構造要素Ｂi より空間的に広い裾野の部分はオリジナル画像信号Ｓorg による形状を維持する。
【００８６】
減算器18が、このクロージング処理によって得られた画像信号Ｓclosing からオリジナル画像信号Ｓorg を減算して石灰化陰影信号Ｓ′（＝Ｓclosing −Ｓorg ）を求めると、この石灰化陰影信号Ｓ′は図４右側第３段目のグラフに示すように石灰化陰影のうち多重構造要素Ｂi より空間的に狭い中心部だけがゼロ以外の特徴的な値を示すため、この信号値の相違により多重構造要素Ｂi より空間的に狭い中心部が検出される。
【００８７】
なおここでは全体の大きさが多重構造要素Ｂi より空間的に狭い範囲の石灰化陰影を図示していないが、この図示しない石灰化陰影についても同様に検出される。この検出される部分以外の部分の信号値はゼロである。
【００８８】
従来は、この石灰化陰影信号Ｓ′に依存した強調係数を式（17）に適用しており、その結果、多重構造要素Ｂi より空間的に狭い範囲に対応した画像部分についての強調係数は大きくなるためより強調度合が強いものとなるが、石灰化陰影の前述の裾野に相当する部分では強調がなされず、結局、石灰化陰影の形状が維持されないという問題があった。
【００８９】
そこで本実施形態においては、膨張処理手段14が得られた石灰化陰影信号Ｓ′に対して縦９画素×横９画素の正方形形状の単一の構造要素を用いたダイレーション処理を施して、石灰化陰影信号Ｓ′を空間的に膨張せしめる（図４右側第４段目のグラフ参照）。
【００９０】
なお、上記膨張処理により得られる膨張処理信号Ｓcalcがゼロ以外の値を採る範囲については、上記構造要素の大きさを適宜変化させることによって所望の広がりにすることができ、本実施形態においては、この構造要素を上述の大きさに設定したことにより、オリジナル画像信号Ｓorg における石灰化陰影の空間的広がりと略同一にしている。ただし、この膨張処理信号Ｓcalcがゼロ以外の値を採る範囲は必ずしもオリジナル画像信号Ｓorg における石灰化陰影の空間的広がりと同一とする必要はなく、オリジナル画像信号Ｓorg における石灰化陰影の空間的広がりと同一若しくはそれ以上の広がりであればよい。
【００９１】
このようにして得られた膨張処理信号Ｓcalcは、第１のルックアップテーブル15によって、その膨張処理信号Ｓcalcに対応した第１の強調係数αcalcに変換される（図４右側第５段目のグラフ）。第１のルックアップテーブルによれば、膨張処理信号Ｓcalcの値に応じて、第１の強調係数αcalcは０から１の範囲の値を採り、膨張処理信号Ｓcalcの値が０の場合は、第１の強調係数αcalcは０である。
【００９２】
一方、エッジ信号依存強調係数算出手段30は、入力されたオリジナル画像信号Ｓorg に対して第２のモーフォロジー演算手段13が縦５画素×横５画素の正方形形状の単一の構造要素によりダイレーション処理を施す。
【００９３】
オリジナル画像信号Ｓorg が図４左側最上段のグラフ（図４右側最上段のグラフと同じ信号）に示すように、乳腺構造と石灰化陰影を有する画像である場合について適用すると、ダイレーション処理によって得られる画像信号Ｓdilationでは、図４左側第２段目のグラフに示すように乳腺構造の一部以外はすべて平坦な信号値となる。
【００９４】
減算器19が、このダイレーション処理によって得られた画像信号Ｓdilationからオリジナル画像信号Ｓorg を減算してエッジ信号Ｓedge（＝Ｓdilation−Ｓorg ）を求めると、このエッジ信号Ｓedgeは図４左側第３段目のグラフに示すように石灰化陰影のエッジ部のみならず、乳腺構造のエッジ部も急峻に信号値が立ち上がる。そして、この立ち上がった部分以外の部分の信号値はゼロである。
【００９５】
このようにして得られたエッジ信号Ｓedgeは、第２のルックアップテーブル16によって、そのエッジ信号Ｓedgeに対応した第２の強調係数αedgeに変換される（図４左側第４段目のグラフ）。
【００９６】
なお本実施形態においては、図４右側第５段目に示す第１の強調係数αcalcがゼロ以外の値を採る範囲が、オリジナル画像信号Ｓorg における石灰化陰影と同一の空間的広がりを示しているが、これは前述したように構造要素の選択による単なる偶然であり、第１の強調係数αcalcがゼロ以外の値を採る範囲が、オリジナル画像信号Ｓorg における石灰化陰影の空間的広がりよりも広くなるように、膨張処理手段14が石灰化陰影信号Ｓ′を膨張せしめるものであればよい。
【００９７】
以上の処理により第１の強調係数αcalc、第２の強調係数αedgeが、いずれも０から１の範囲内の値を採るものとして求められ（図４右第５段、図４左第４段）、これらは乗算器41に入力されて互いに乗じられ（図４下下段）、１つの強調係数αが求められる。
【００９８】
この作用により、第１の強調係数αcalcまたは第２の強調係数αedgeのうち少なくとも一方の値が０の画像部分では、これらの積は０となり、第１の強調係数αcalcおよび第２の強調係数αedgeがいずれも０以外の値を採る画像部分でのみこれらの積は０以外の値となる。
【００９９】
すなわち従来の方法によれば、エッジ信号依存強調係数算出手段30単独では所望の強調対象である石灰化陰影のみならず強調対象ではない乳腺構造のエッジ部分においても強調係数が０以外の値を採り、また石灰化陰影強調係数算出手段20単独では構造要素より大きな石灰化陰影については、信号値の変動が殆どない平坦部分の大きさが構造要素より小さい範囲で信号値が変動する部分（陰影の中心部に相当する部分）だけが０を採っていた。
【０１００】
これに対して本実施形態においては、最終的に求められた強調係数αは、石灰化陰影のエッジ部に相当する輪郭部およびその内部（裾野部分含む）において０以外の値を採り、乳腺構造のエッジ部等の強調対象でない部分では０を採る。
【０１０１】
このような値を採る強調係数αと、オリジナル画像信号Ｓorg と、高周波成分（Ｓorg −Ｓus）とが、乗算器42および加算器43により式（17）にしたがって乗算および加算され、処理済み画像信号Ｓprocが出力される。
【０１０２】
この式（17）において高周波成分（Ｓorg −Ｓus）を強調する強調係数は上述したように、石灰化陰影のエッジ部に相当する輪郭部およびその内部においてのみ０以外の値を採るため、石灰化陰影の輪郭部およびその内部だけが選択的に強調処理されるとともに、従来の強調処理方法、装置よりも石灰化陰影の分離度を高めることができる。
【０１０３】
なお、図４下側下段のグラフに示した強調係数αには定数（Const.）が乗じられているが、これは図１の乗算器41の次段に設けられた増巾器（符番せず）の作用によるものであり、増幅器を設けない場合はConst.＝１とすればよい。
【０１０４】
以上説明した実施形態は本発明の画像処理装置の一実施形態に過ぎず、この構成に限るものではない。例えば、第１のモーフォロジー演算手段12を、線形の３画素４方向の多重構造要素を用いてクロージング処理を行なうものではなく、縦３画素×横３画素の正方形状の単一の構造要素を用いてクロージング処理を行なうモーフォロジー演算手段12′を用いてもよく（図５参照）、この場合も上記実施形態と同様の効果を得ることができる。
【０１０５】
また膨張処理手段14も、必ずしもダイレーション処理を行なうものに限らず、非鮮鋭マスク処理を行なう膨張処理手段14′を用いることもできる（図６参照）。
【０１０６】
さらに、縦３画素×横３画素の正方形状の単一の構造要素を用いてクロージング処理を行なう第１のモーフォロジー演算手段12′と、非鮮鋭マスク処理を行なう膨張処理手段14′との組み合わせを同時に用いる構成（図７参照）であってもよい。
【０１０７】
なお、上記各実施形態において、ローパスフィルター11のボケマスクサイズを縦３画素×横３画素として最終的に得られた石灰化陰影のみを強調処理した画像信号Ｓprocに対して、さらに縦３画素×横３画素以上の大きなボケマスクを用いて低周波数強調処理（通常の濃度依存強調処理）を施すことにより、乳癌における他の特徴的形態の一つである腫瘤陰影（濃度値の勾配線が陰影の中心に集中する形態の陰影）をも強調処理することができ、腫瘤陰影および石灰化陰影が強調処理された、診断性能の高い画像処理を行なうことができる。
【０１０８】
以上の各実施形態に置いて用いた第１および第２のルックアップテーブル15，16における各閾値Thld１，Thld１′，Thld２，Thld２′は経験的に求められるものであるが、各ルックアップテーブルに入力される信号のヒストグラムの極大値をThld１，Thld１′、極小値をThld２，Thld２′とするなど、ヒストグラム解析に基づいて設定することもできる。
【図面の簡単な説明】
【図１】本発明の画像処理装置の第１実施形態を示すブロック図
【図２】第１のルックアップテーブルを表す関数のグラフ
【図３】第２のルックアップテーブルを表す関数のグラフ
【図４】本発明の画像処理装置の作用を説明するための信号値を示すグラフ
【図５】本発明の画像処理装置の第２実施形態を示すブロック図
【図６】本発明の画像処理装置の第３実施形態を示すブロック図
【図７】本発明の画像処理装置の第４実施形態を示すブロック図
【図８】モーフォロジー演算の基本的な作用を説明する図
【図９】モーフォロジーフィルターにおける構造要素Ｂi （ｉ＝１，２，…，Ｍ；Ｍ＝４）を示す図
【図１０】オリジナル画像信号が示す形状の相違に応じた処理結果の相違を説明するための図
【図１１】（Ａ）図１に示した実施形態の画像処理装置で用いられる線形の３画素４方向の多重構造要素を示す図、（Ｂ）５画素８方向の多重構造要素を示す図、（Ｃ）７画素８方向の多重構造要素を示す図
【符号の説明】
11    ローパスフィルタ
12，13    モーフォロジー演算手段
12′      モーフォロジー演算手段（多重構造要素）
14    膨張処理手段（ダイレーション処理）
14′  膨張処理手段（非鮮鋭マスク処理）
15，16    ルックアップテーブル
17，18，19    減算器
20    石灰化陰影依存強調係数算出手段
30    エッジ信号依存強調係数算出手段
41，42    乗算器
43    加算器
Ｓorg   オリジナル画像信号
Ｓus    ボケマスク信号
Ｓproc  処理済画像信号
【図１】
0003669455.tif 000008
【図２】
0003669455.tif 000009
【図３】
0003669455.tif 000010
【図４】
0003669455.tif 000011
【図５】
0003669455.tif 000012
【図６】
0003669455.tif 000013
【図７】
0003669455.tif 000014
【図８】
0003669455.tif 000015
【図９】
0003669455.tif 000016
【図１０】
0003669455.tif 000017
【図１１】
0003669455.tif 000018