(57)【特許請求の範囲】
【請求項１】
画像を表すオリジナルの画像信号Ｄorg に対して、多重構造要素Ｂi およびスケール係数λを用いたモーフォロジー演算を施すことにより、前記画像信号が空間的に前記多重構造要素Ｂi より小さく変動する画像部分の画像信号値の変動量を表すモーフォロジー信号Ｄmor を抽出し、
前記オリジナルの画像信号Ｄorg の、超低空間周波数に対応する非鮮鋭マスク信号Ｄusを求め、
前記モーフォロジー信号Ｄmor に基づく関数β（Ｄmor ）および前記オリジナルの画像信号Ｄorg と前記非鮮鋭マスク信号Ｄusとの差信号（Ｄorg −Ｄus）に基づく関数ｆ（Ｄorg −Ｄus）を用いて、前記オリジナルの画像信号Ｄorg に対して、
【数１】
0003731760.tif 000018
なる演算を行って、前記画像部分のオリジナルの画像信号Ｄorgに、前記モーフォロジー信号Ｄmorの大きさに応じた強調処理を施す画像処理方法であって、
前記関数β（Ｄmor）が、単調増加関数であり、かつ、前記モーフォロジー信号Ｄmorの値が小さい領域Ｃ1および前記モーフォロジー信号Ｄmorの値が大きい領域Ｃ３における関数β（Ｄmor）の増加率が、前記モーフォロジー信号Ｄmorの値が中間の領域Ｃ２における関数β（Ｄmor）の増加率より小さい関数であることを特徴とする画像処理方法。
【請求項２】
前記モーフォロジー演算が、下記式（２）で示されるものであり、前記画像部分が、該画像部分の周囲の画像部分よりも前記オリジナルの画像信号Ｄorg の値が大きい部分であることを特徴とする請求項１記載の画像処理方法。
【数２】
0003731760.tif 000019
【請求項３】
前記モーフォロジー演算が、下記式（３）で示されるものであり、前記画像部分が、該画像部分の周囲の画像部分よりも前記オリジナルの画像信号Ｄorg の値が小さい部分であることを特徴とする請求項１記載の画像処理方法。
【数３】
0003731760.tif 000020
【請求項４】
前記モーフォロジー演算が、下記式（４）で示されるものであり、前記画像部分が、該画像部分の周囲の画像部分よりも前記オリジナルの画像信号Ｄorg の値が大きい部分であることを特徴とする請求項１記載の画像処理方法。
【数４】
0003731760.tif 000021
【請求項５】
前記モーフォロジー演算が、下記式（５）で示されるものであり、前記画像部分が、該画像部分の周囲の画像部分よりも前記オリジナルの画像信号Ｄorg の値が小さい部分であることを特徴とする請求項１記載の画像処理方法。
【数５】
0003731760.tif 000022
【請求項６】
前記モーフォロジー演算が、下記式（６）で示されるものであり、前記画像部分が、該画像部分の周囲の画像部分よりも前記オリジナルの画像信号Ｄorg の値が大きい部分であることを特徴とする請求項１記載の画像処理方法。
【数６】
0003731760.tif 000023
【請求項７】
前記モーフォロジー演算が、下記式（７）で示されるものであり、前記画像部分が、該画像部分の周囲の画像部分よりも前記オリジナルの画像信号Ｄorg の値が小さい部分であることを特徴とする請求項１記載の画像処理方法。
【数７】
0003731760.tif 000024
【請求項８】
前記モーフォロジー演算が、下記式（８）で示されるものであることを特徴とする請求項１記載の画像処理方法。
【数８】
0003731760.tif 000025
【発明の詳細な説明】
【０００１】
【産業上の利用分野】
本発明は画像処理方法に関し、詳細には画像のうち、異常陰影等の特定の画像部分だけを強調処理する画像処理方法に関するものである。
【０００２】
【従来の技術】
従来より、種々の画像取得方法により得られた画像を表す画像信号に対して、階調処理や周波数処理等の画像処理を施し、画像の観察読影性能を向上させることが行われている。特に人体を被写体とした放射線画像のような医用画像の分野においては、医師等の専門家が、得られた画像に基づいて患者の疾病や傷害の有無を的確に診断する必要があり、その画像の読影性能を向上させる画像処理は不可欠なものとなっている。
【０００３】
この画像処理のうち、いわゆる周波数強調処理としては、例えば特開昭61-169971 号に示されるように、原画像の濃度値等の画像信号（オリジナル画像信号という）Ｄorg を、
Ｄproc＝Ｄorg ＋β×（Ｄorg −Ｄus）                    （９）
なる画像信号Ｄprocに変換するものが知られている。ここでβは周波数強調係数、Ｄusは非鮮鋭マスク（いわゆるボケマスク）信号である。このボケマスク信号Ｄusは、２次元に配置された画素に対してオリジナル画像信号Ｄorg を中心画素とするＮ列×Ｎ行（Ｎは奇数）の画素マトリクスからなるマスク、すなわちボケマスクを設定し、
Ｄus＝（ΣＤorg ）／Ｎ2                                 （10）
（ただし、ΣＤorg はボケマスク内各画素の画像信号の和）
等として求められる超低空間周波数成分である。
【０００４】
式（９）の第２項括弧内の値（Ｄorg −Ｄus）は、オリジナル画像信号から超低空間周波数成分であるボケマスク信号を減算したものであるから、オリジナル画像信号のうちの、超低空間周波数成分を除去した比較的高い周波数成分を選択的に抽出することができる。この比較的高い周波数成分に周波数強調係数βを乗じたうえで、オリジナル画像信号を加算することにより、この比較的高い周波数成分を強調することができる。
【０００５】
一方、画像のうち異常な陰影等の特定の画像部分だけを選択的に抽出する、モーフォロジー（Morphology；モフォロジーまたはモルフォロジーとも称する）のアルゴリズムに基づく処理（以下、モーフォロジー演算またはモーフォロジー処理という）が知られている。このモーフォロジー処理は、特に乳癌における特徴的形態である微小石灰化像を検出するのに有効な手法として研究されているが、対象画像としては、このようなマンモグラムにおける微小石灰化像に限るものではなく、検出しようとする特定の画像部分（異常陰影等）の大きさや形状が予め分かっているものについては、いかなる画像に対しても適用することができる。
【０００６】
以下、このモーフォロジー処理を、マンモグラムにおける微小石灰化像の検出に適用した例により、モーフォロジー処理の概要について説明する。
【０００７】
（モーフォロジーの基本演算）
モーフォロジー処理は一般的にはＮ次元空間における集合論として展開されるが、直感的な理解のために２次元の濃淡画像を対象として説明する。
【０００８】
濃淡画像を座標（ｘ，ｙ）の点が濃度値ｆ（ｘ，ｙ）に相当する高さをもつ空間とみなす。ここで、濃度値ｆ（ｘ，ｙ）は、ＣＲＴに表示するための信号のように、輝度の高いもの程大きな値となる高輝度高信号レベルの信号とする。
【０００９】
まず、簡単のために、その断面に相当する１次元の関数ｆ（ｘ）を考える。モーフォロジー演算に用いる構造要素ｇは次式（11）に示すように、原点について対称な対称関数
【００１０】
【数９】
0003731760.tif 000002
【００１１】
であり、定義域内で値が０で、その定義域Ｇが下記式（12）であるとする。
【００１２】
【数１０】
0003731760.tif 000003
【００１３】
このとき、モーフォロジー演算の基本形は式（13）〜（16）に示すように、非常に簡単な演算となる。
【００１４】
【数１１】
0003731760.tif 000004
【００１５】
すなわち、dilation（ダイレーション）処理は、注目画素を中心とした、±ｍ（構造要素Ｂに応じて決定される値）の幅の中の最大値を探索する処理であり（図４（Ａ）参照）、一方、erosion （エロージョン）処理は、注目画素を中心とした、±ｍの幅の中の最小値を探索する処理である（図４（Ｂ）参照）。また、opening （またはclosing ）処理は最小値（または最大値）の探索の後に、最大値（または最小値）を探索することに相当する。すなわちopening （オープニング）処理は、低輝度側から濃度曲線ｆ（ｘ）を滑らかにし、マスクサイズ２ｍより小さい凸状の濃度変動部分（周囲部分よりも輝度が高い部分）を取り除くことに相当する（図４（Ｃ）参照）。一方、closing （クロージング）処理は、高輝度側から濃度曲線ｆ（ｘ）を滑らかにし、マスクサイズ２ｍより小さい凹状の濃度変動部分（周囲部分よりも輝度が低い部分）を取り除くことに相当する（図４（Ｄ）参照）。
【００１６】
なお、構造要素ｇが原点に対して対称ではない場合の、式（13）に示すダイレーション演算をMinkowski （ミンコフスキー）和、式（14）に示すエロージョン演算をMinkowski 差という。
【００１７】
ここで、濃度値ｆ（ｘ）が、ネガフイルムに記録するための信号のように、濃度の高いもの程大きな値となる高濃度高信号レベルの信号の場合においては、輝度と濃度との関係が逆転するため、高濃度高信号レベルの信号におけるダイレーション処理は、高輝度高信号レベルにおけるエロージョン処理（図４（Ｂ））と一致し、高濃度高信号レベルの信号におけるエロージョン処理は、高輝度高信号レベルにおけるダイレーション処理（図４（Ａ））と一致し、高濃度高信号レベルの信号におけるオープニング処理は、高輝度高信号レベルにおけるクロージング処理（図４（Ｄ））と一致し、高濃度高信号レベルの信号におけるクロージング処理は、高輝度高信号レベルにおけるオープニング処理（図４（Ｃ））と一致する。
【００１８】
なお、本項では高輝度高信号レベルの画像信号（輝度値）の場合について説明する。
【００１９】
（石灰化陰影検出への応用）
石灰化陰影の検出には、原画像から平滑化した画像を引き去る差分法が考えられる。単純な平滑化法では石灰化陰影と細長い形状の非石灰化陰影（乳腺、血管および乳腺支持組織等）との識別が困難であるため、東京農工大の小畑らは、多重構造要素を用いたオープニング演算に基づく下記式（17）で表されるモーフォロジーフィルターを提案している（「多重構造要素を用いたモルフォロジーフィルタによる微小石灰化像の抽出」電子情報通信学会論文誌 D-II Vol.J75-D-II No.7 P1170 〜1176  1992年７月等）。
【００２０】
【数１２】
0003731760.tif 000005
【００２１】
ここでＢi （ｉ＝１，２，…，Ｍ）は、例えば図５に示す直線状の４つ（この場合Ｍ＝４）の構造要素（これら４つの構造要素の全体をもって多重構造要素という）である。構造要素Ｂi を検出対象である石灰化陰影よりも大きく設定すれば、上記オープニング演算による処理で、構造要素Ｂi よりも細かな凸状の信号変化部分である石灰化陰影は取り除かれる。一方、細長い形状の非石灰化陰影はその長さが構造要素Ｂi よりも長く、その傾きが４つの構造要素Ｂi のいずれかに一致すればオープニング処理（式（17）の第２項の演算）をしてもそのまま残る。したがってオープニング処理によって得られた平滑化画像（石灰化陰影のみが取り除かれた画像）を原画像ｆから引き去ることで、小さな石灰化陰影のみが含まれる画像が得られる。これが式（17）の考え方である。
【００２２】
なお、前述したように、高濃度高信号レベルの信号の場合においては、石灰化陰影は周囲の画像部分よりも濃度値が低くなり、石灰化陰影は周囲部分に対して凹状の信号変化部分となるため、オープニング処理に代えてクロージング処理を適用し、式（17）に代えて式（18）を適用する。
【００２３】
【数１３】
0003731760.tif 000006
【００２４】
このように、モーフォロジー処理は、
（１）石灰化陰影そのものの抽出に有効であること
（２）複雑なバックグラウンド情報に影響されにくいこと
（３）抽出した石灰化陰影がひずまないこと
などの特徴がある。すなわち、この手法は一般の微分処理に比べて、石灰化陰影のサイズ・形状・濃度分布などの幾何学的情報をより良く保った検出が可能である。
【００２５】
【発明が解決しようとする課題】
ところで上述したように、画像の読影性能を向上させるには、対象となる画像に対して画像処理を行うことが不可欠となっているが、特開平２−1078号に開示されているように、単に濃度依存による強調処理では、例えばマンモグラムにおける放射線ノイズ成分のような、画像読影の障害となる成分まで強調されるため、読影性能をむしろ低下させることになる。
【００２６】
また特公昭60−192482号、特開平２−120985号、特表平３−502975号等に開示されているように、画像信号の分散値に依存した強調処理では、局所的に濃度変化の大きい画像部分が強く強調されるため、その付近でアンダーシュート、オーバーシュートが相対的に目立ち、特にＸ線画像に関しては高濃度側でアーチファクトが発生しやすいという問題がある。
【００２７】
本発明は上記事情に鑑みなされたものであって、ノイズ成分等の画像読影に不要な成分を強調することなく、注目する特定の画像部分だけを効率よく強調処理し、アーチファクトの発生を抑制した画像処理方法を提供することを目的とするものである。
【００２８】
【課題を解決するための手段】
本発明の画像処理方法は、画像を表すオリジナルの画像信号Ｄorg に対して、多重構造要素Ｂi およびスケール係数λを用いたモーフォロジー演算を施すことにより、前記画像信号が空間的に前記多重構造要素Ｂi より小さく変動する画像部分の画像信号値の変動量を表すモーフォロジー信号Ｄmor を抽出し、
前記オリジナルの画像信号Ｄorg の、超低空間周波数に対応する非鮮鋭マスク信号Ｄusを求め、
前記モーフォロジー信号Ｄmor に基づく関数β（Ｄmor ）および前記オリジナルの画像信号Ｄorg と前記非鮮鋭マスク信号Ｄusとの差信号（Ｄorg −Ｄus）に基づく関数ｆ（Ｄorg −Ｄus）を用いて、前記オリジナルの画像信号Ｄorg に対して、
【００２９】
【数１】
0003731760.tif 000007
【００３０】
なる演算を行って、前記画像部分のオリジナルの画像信号Ｄorgに、前記モーフォロジー信号Ｄmorの大きさに応じた強調処理を施す画像処理方法であって、前記関数β（Ｄmor）が、単調増加関数であり、かつ、前記モーフォロジー信号Ｄmorの値が小さい領域Ｃ1および前記モーフォロジー信号Ｄmorの値が大きい領域Ｃ３における関数β（Ｄmor）の増加率が、前記モーフォロジー信号Ｄmorの値が中間の領域Ｃ２における関数β（Ｄmor）の増加率より小さい関数であることを特徴とするものである。なお、関数β（Ｄmor ）は、図３（ａ）、（ｂ）に示すように、｜Ｄmor ｜に対する関数である。
【００３１】
また、構造要素Ｂとしては、例えば、正方形、長方形、円形、楕円形、または菱形等の上下左右対称の要素が望ましい。
【００３２】
さらに、モーフォロジー演算としては、下記式（２）〜（８）で示される種々のものを適用することができる。
【００３３】
【数２】
0003731760.tif 000008
【００３４】
【数３】
0003731760.tif 000009
【００３５】
【数４】
0003731760.tif 000010
【００３６】
【数５】
0003731760.tif 000011
【００３７】
【数６】
0003731760.tif 000012
【００３８】
【数７】
0003731760.tif 000013
【００３９】
【数８】
0003731760.tif 000014
【００４０】
すなわち、式（２）で表されるモーフォロジー演算を適用することにより、モーフォロジー信号Ｄmor として、オリジナルの画像信号Ｄorg の値が周囲の画像部分よりも大きく、かつ前記多重構造要素Ｂi より小さく変動する画像部分（例えば、高輝度高信号レベルの画像信号における石灰化陰影）を構成する画素の信号を抽出することができ、この画像部分を効果的に強調処理することができる。
【００４１】
なお、上記式（２）〜（８）において、
【００４２】
【数１４】
0003731760.tif 000015
【００４３】
また、式（３）で表されるモーフォロジー演算を適用することにより、モーフォロジー信号Ｄmor として、オリジナルの画像信号Ｄorg の値が周囲の画像部分よりも小さく、かつ前記多重構造要素Ｂi より小さく変動する画像部分（例えば、高濃度高信号レベルの画像信号における石灰化陰影）を構成する画素の信号を抽出することができ、この画像部分を効果的に強調処理することができる。
【００４４】
式（４）で表されるモーフォロジー演算を適用することにより、モーフォロジー信号Ｄmor として、オリジナルの画像信号Ｄorg の値が周囲の画像部分よりも大きくかつ前記多重構造要素Ｂi より小さく変動する画像部分や、輝度（濃度）が急激に変化するエッジ部分を構成する画素の信号を抽出することができ、このような画像部分を効果的に強調処理することができる。
【００４５】
式（５）で表されるモーフォロジー演算を適用することにより、モーフォロジー信号Ｄmor として、オリジナルの画像信号Ｄorg の値が周囲の画像部分よりも小さくかつ前記多重構造要素Ｂi より小さく変動する画像部分や、輝度（濃度）が急激に変化するエッジ部分を構成する画素の信号を抽出することができ、このような画像部分を効果的に強調処理することができる。
【００４６】
式（６）で表されるモーフォロジー演算を適用することにより、モーフォロジー信号Ｄmor として、オリジナルの画像信号Ｄorg の値が周囲の画像部分よりも大きくかつ前記多重構造要素Ｂi より小さく変動する、濃度（輝度）変化の大きい画像部分（例えばオリジナルの画像信号Ｄorg が表す画像の骨格部分）を構成する画素の信号を抽出することができ、このような画像部分（骨格部分）を効果的に強調処理することができる。式（６）を具体的に適用した例を図６に示す。原画像Ｘに対する構造要素Ｂ（半径ｒの円形構造とする）での、エロージョン処理した画像と、エロージョン処理した画像に対するオープニング処理した画像との差信号のλ（λ＝１，２，…，Ｎ）回までの和集合は、骨格部分ａおよびｂとなる。
【００４７】
式（７）で表されるモーフォロジー演算を適用することにより、モーフォロジー信号Ｄmor として、オリジナルの画像信号Ｄorg の値が周囲の画像部分よりも小さくかつ前記多重構造要素Ｂi より小さく変動する、濃度（輝度）変化の大きい画像部分（例えばオリジナルの画像信号Ｄorg が表す画像の骨格部分）を構成する画素の信号を抽出することができ、このような画像部分（骨格部分）を効果的に強調処理することができる。
【００４８】
式（６）、（７）で表されるモーフォロジー演算を一般にスケルトン処理と称し、このスケルトン処理によれば、特に骨りょう部分の画像信号に適用した場合、骨格要素だけを選択的に、効果的に強調処理することができる。
【００４９】
なお、モーフォロジー演算として式（２）〜（７）を適用する場合は、関数ｆ（Ｄorg −Ｄus）としては、ｆ（Ｄorg −Ｄus）＝Ｄorg −Ｄusとしてもよい。
【００５０】
式（８）で表されるモーフォロジー演算を適用することにより、局所的に輝度（濃度）変化の大きい画像部分を構成する画素の信号を抽出することができ、このような画像部分を効果的に強調処理することができる。
【００５１】
なお、モーフォロジー演算として式（８）を適用する場合はオーバーシュート、アンダーシュートを抑制するため、関数ｆ（Ｄorg −Ｄus）は、Ｄorg が所定の大きさ以下の範囲、あるいは所定の大きさ以上の範囲において、例えば図７に示すように、出力を抑制する関数とすることが必要である。
【００５２】
【作用および発明の効果】
本発明の画像処理方法は、画像を表すオリジナル画像信号Ｄorg に対して、多重構造要素Ｂi およびスケール係数λを用いたモーフォロジー演算を施すことにより、画像信号が空間的に多重構造要素Ｂi より小さく変動する画像部分の画像信号値の変動量を表すモーフォロジー信号Ｄmorを抽出する。
【００５３】
一方、式（１）第２項の計算処理により、オリジナル画像信号Ｄorg から、その超低空間周波数成分Ｄusを減算することにより、オリジナル画像信号Ｄorg のうちの比較的高い周波数成分（超低空間周波数成分を除いたものをいう）だけを抽出することができる。この抽出された比較的高い周波数成分には、いわゆる高周波成分である放射線ノイズも含まれる。
【００５４】
ここで、式（１）に示すように、全体画像の各画素がモーフォロジー演算で得られた石灰化陰影等の、空間的に多重構造要素Ｂi より小さく変動する特定の画像部分、を構成する画素か否かに応じた信号Ｄmor に基づく強調係数β（Ｄmor ）で高周波成分（Ｄorg −Ｄus）を強調するため、この高周波成分（Ｄorg −Ｄus）に例えば量子ノイズ等の不要な成分が含まれていても、その画素が石灰化陰影等の画像部分を構成するものでない場合は、その画素についてのβ（Ｄmor ）の値は小さいため、その画素に対する強調度合は小さいものとなる。
【００５５】
一方、画素が石灰化陰影等の画像部分を構成するものである場合は、その画素についてのβ（Ｄmor ）の値は大きいため、その画素に対する強調度合は大きいものとなる。
【００５６】
したがって、画像の高周波成分（Ｄorg −Ｄus）に放射線ノイズが含まれているか否かに拘らず、石灰化陰影等の特定の画像部分であるか否かに応じた関数β（Ｄmor ）により、特定の画像部分のオリジナルの画像信号Ｄorgに、前記モーフォロジー信号Ｄmorの大きさに応じた強調処理を施すことができる。
【００５７】
また、モーフォロジー信号Ｄmorの大きさに応じて特定の画像部分が強く強調されるため、アンダーシュート、オーバーシュートを抑制することができアーチファクトの発生を低減することができる。さらに、画像信号の分散値を計算する処理に代えて最大値または最小値を計算する処理を行うため、計算処理に掛かる演算時間を短縮することができる。
【００５８】
以下、モーフォロジー演算の一例である式（３）を具体的に説明する。
【００５９】
すなわち、高濃度高信号レベルの画像信号である濃度値Ｄorg についてのモーフォロジー演算によれば、例えば図８（１）の実線に示すような濃度値Ｄorg の分布を有する画像データに対して、同図（２）に示すような直線状の３画素の構造要素Ｂで、Ｄorg ［＋］Ｂの演算（ミンコフスキー和）を行うことにより、ある注目画素の濃度値Ｄi は、その注目画素を中心として互いに隣接する３画素（構造要素Ｂにより決定される）の中の最大値Ｄi+1 を採用したＤi ′に変換される。この演算を全画素について行うことにより、濃度値Ｄorg ′の分布を有する同図（１）の破線で示す画像データに変換される。
【００６０】
次に、このミンコフスキー和の演算で得られた濃度値（Ｄorg ［＋］Ｂ）の構造要素Ｂによるミンコフスキー差（Ｄorg ［＋］Ｂ）［−］Ｂを考えると、同図（１）の破線で示された注目画素の濃度値Ｄi ′は、その注目画素を中心として互いに隣接する３画素の中の最小値Ｄi-1 ′を採用したＤi ″（＝Ｄi ）に変換される。この演算を全画素について行うことにより、濃度値Ｄorg ″の分布を有する同図（１）の一点鎖線で示す画像データに変換される。この一点鎖線で示された画像データは、もとの実線のオリジナルの画像データに対して、構造要素Ｂよりも細かい信号値の変化の画像部分が消え、構造要素Ｂよりも大きい信号値の変化の画像部分は演算前の元の状態を保持することを示している。すなわち、以上の処理（クロージング処理）は、画像濃度の分布を高濃度側から平滑化する処理である。
【００６１】
このようにクロージング処理で得られた値（Ｄorg ［＋］Ｂ）［−］Ｂを元の画像信号Ｄorg から差し引くことにより得られた値Ｄmor は、上記クロージング処理で消された細かい信号値の変化の画像部分を表す。
【００６２】
ここで、本来、画像信号は２次元の要素である位置（ｘ，ｙ）と、３次元目の要素である信号値ｆ（ｘ，ｙ）を有するが、上記説明においては、理解の容易化のために、この２次元上に展開された画像の所定の断面に現れた、１次元状の画像信号分布曲線について説明した。
【００６３】
したがって、本発明の画像処理方法は、以上の説明を２次元画像に拡大適用したものであり、式（１）に示した多重構造要素Ｂi は、このような断面におけるモーフォロジー演算を２次元面に拡大適用する場合に、この２次元面内での向きが互いに異なる構造要素Ｂとして準備されたｉ個の構造要素Ｂの集合を意味するものである。
【００６４】
またこれらすべての多重構造要素Ｂi についてクロージング処理を行った結果、多重構造要素Ｂi のうちのいずれかと、その延びる方向が一致し、かつその大きさよりも大きく変化する画像部分については、式（３）の第２項ｍｉｎ｛（Ｄorg ［＋］λＢi ）［−］λＢi ｝の値がＤorg そのものとなるため、Ｄmor の値はゼロとなり、その部分については強調処理されない。
【００６５】
なお、スケール係数λは上記ミンコフスキー和の演算およびミンコフスキー差の演算を行う回数を意味し、回数を増加するに応じて平滑化の程度が進む。
【００６６】
このような本発明のモーフォロジー演算を利用した強調処理と、従来のボケマスク処理による強調処理との比較を、図９に示す。
【００６７】
【実施例】
以下、本発明の画像処理方法を用いた計算機支援画像診断装置について図面を用いて説明する。
【００６８】
図１は本実施例の計算機支援画像診断装置の概略構成を示すブロック図、図２はこの計算機支援画像診断装置により画像診断に供される乳房の放射線画像（マンモグラム）を示す図である。図示の計算機支援画像診断装置は、マンモグラムの全体の放射線画像（全体画像）Ｐを表す、各画素の濃度値Ｄorg の集合である画像データ（全体画像データ）Ｓを記憶する全体画像メモリ10、全体画像データＳに基づいて、もしくは一旦全体画像メモリ10に記憶された全体画像データＳに基づいて全体画像Ｐを表示するＣＲＴ等の全体画像表示手段30、全体画像メモリ10に記憶された全体画像データＳに基づいて、全体画像Ｐのうち異常陰影Ｐ1 を検出する異常陰影検出手段40、異常陰影検出手段40により異常陰影Ｐ1 が検出されたか否かを判定する判定手段50、異常陰影Ｐ1 が検出されたと判定手段50により判定された場合に、全体画像メモリ10に記憶された全体画像データＳのうち、異常陰影Ｐ1 を含む局所領域の画像Ｐ2 を表す画像データ（局所画像データ）Ｓ2 を抽出する局所領域抽出手段60、局所領域抽出手段60により抽出された局所画像データＳ2 に基づいた局所領域の画像Ｐ2 のうち異常陰影Ｐ1 が、全体画像表示手段30に表示される全体画像Ｐよりも読影性能が向上するように、異常陰影を示す画像データ（異常陰影画像データ）Ｓ1 に対して画像強調処理を行う局所画像強調手段70、およびこの画像強調処理のなされた局所画像データＳ2 に基づいて局所領域の画像Ｐ2 を表示するＣＲＴ等の局所画像表示手段90、を備えた構成である。
【００６９】
なお、異常陰影は石灰化陰影、異常陰影検出手段40はモーフォロジー演算を行うモーフォロジーフィルターであって、石灰化陰影を示す画素について、その画素についての式（３）で示した、石灰化陰影を構成する画素であるか否かを示す信号Ｄmor として出力する。
【００７０】
ただし、対象画像は、本実施例のような医用画像に限るものではなく、工業製品の検査用画像等であってもよい。例えば、内部に巣のある鋳物製品のＸ線画像について、異常陰影が当該巣の陰影であってもよい。
【００７１】
また上記局所領域とは、異常陰影である石灰化陰影を含むこの石灰化陰影近傍の領域をいうものとする。
【００７２】
なお、実施例記載中において、画像を構成する各画素の画像データを濃度信号値Ｄorg と表記し、これら画素の集合により構成される領域の画像データを画像データＳと表記するものとする。また、濃度値Ｄorg は高濃度高信号レベルの信号値である。
【００７３】
ここでモーフォロジーフィルター40は、前述した特定の画像部分の検出処理のアルゴリズムをいうが、本実施例におけるモーフォロジーフィルター40は、このアルゴリズム自体を指すのではなく、このアルゴリズムによる石灰化陰影の検出処理（式（３）で表す処理）をなす手段を意味する。
【００７４】
また、局所画像強調手段70は詳しくは、局所画像データＳ2 を構成する各画素（濃度値Ｄorg ）について、その画素を中心としたＮ列×Ｎ行（Ｎは例えば「５」などの奇数）の画素マトリクスからなるマスク（以下、単にボケマスクという）信号Ｄusを下記式（10）により算出し、その後に、濃度値Ｄorg からボケマスク信号値Ｄusを減算（Ｄorg −Ｄus）して高周波成分を算出する手段71と、
Ｄus＝（ΣＤorg ）／Ｎ2                                 （10）
（ただし、ΣＤorg はボケマスク内各画素の画像信号の和）
モーフォロジーフィルター40により抽出された石灰化陰影を示す画素についての出力Ｄmor を、図３（ａ）で表される単調に増加するβ（Ｄmor ）に変換して出力する変換テーブル72と、
この出力β（Ｄmor ）と前述の高周波成分（Ｄorg −Ｄus）とを乗じて、周囲に比べて濃度値の低い石灰化陰影に対して、より強い周波数強調を施す石灰化陰影強調手段73とを備えた構成である。
【００７５】
以下、本実施例の計算機支援画像診断装置の作用について説明する。
【００７６】
全体画像メモリ10に光磁気ディスクや画像読取装置等から、内部に石灰化部を有する乳房を含む全体画像Ｐを表す全体画像データＳが入力される。また、この全体画像データＳは直接、全体画像表示手段30にも入力され（図１のＡの経路）、もしくは一旦全体画像メモリ10に記憶されたものとして全体画像表示手段30に入力されて（図１のＢの経路）、全体画像表示手段30は、この全体画像データＳに基づいて全体画像Ｐを表示する。
【００７７】
一方、全体画像メモリ10に記憶された全体画像データＳはモーフォロジーフィルター40にも入力される。モーフォロジーフィルター40は、式（３）にしたがって、石灰化陰影Ｐ1 を示す画像データ（以下、石灰化画像データという）Ｓ1 を検出する。
【００７８】
すなわち、図２（１）に示したマンモグラムのI-I 線断面における濃度値Ｄorg の分布は同図（２）に示すものとなるが、微小石灰化陰影Ｐ1 は、その周辺の濃度値の分布が多重構造要素Ｂi の大きさに比べて細かく変化するため、式（３）のＤmor の値はゼロ以外の所定の値を有し、クロージング処理により平滑化される。一方、構造要素Ｂよりも大きい変動である血管や乳腺の画像Ｐ3 は、その濃度値の分布が多重構造要素Ｂi の大きさに比べて大きいため、Ｄmor の値がゼロとなり、クロージング処理で平滑化されない。
【００７９】
これにより、モーフォロジーフィルター40により微小石灰化陰影Ｐ1 を示す画像データＳ1 の画素（位置）が特定される。判定手段50は、モーフォロジーフィルター40によりこの微小石灰化陰影Ｐ1 を示す微小石灰化像データＳ1 が検出されたことを判定するとともに、微小石灰化像データＳ1 の画素位置を特定する位置データ（以下、石灰化画素位置データという）Ｄ1 、およびモーフォロジーフィルター40により抽出された微小石灰化陰影の濃度値の変動Ｄmor を、局所領域抽出手段60に入力する。
【００８０】
判定手段50は、モーフォロジーフィルター40によりこの石灰化陰影Ｐ1 を示す石灰化画像データが検出されなかったと判定した場合には、石灰化画像データＳ1 の画素位置を特定する石灰化画素位置データＤ1 を出力せずに処理は終了する。
【００８１】
一方、石灰化画像データが検出されたと判定した場合は、局所領域抽出手段60には全体画像メモリ10に記憶された全体画像データＳも入力され、局所領域抽出手段60は入力された全体画像データＳのうち、石灰化画素位置データＤ1 に基づいて石灰化画像データＳ1 の画素を含む近傍の画素（これらの画素の集合としての局所領域）を、予め設定された処理手順にしたがって特定したうえで、この局所領域の画像Ｐ2 を表す局所画像データＳ2 を抽出する。
【００８２】
抽出された局所画像データＳ2 および信号Ｄmor は局所画像強調手段70に入力される。
【００８３】
局所画像強調手段70に入力された局所画像データＳ2 を構成する各画素（濃度値Ｄorg ）について、まず高周波成分を算出する手段71により、ボケマスク信号Ｄusが算出されたうえで、高周波成分（Ｄorg −Ｄus）が算出される。次に、モーフォロジーフィルター40からの入力Ｄmor が変換テーブル72によりβ（Ｄmor ）に変換される。この変換テーブル72は、図３（ａ）に示すように単調増加関数である。すなわち、信号Ｄmor が大きな値を採る場合は、その画素が石灰化陰影に対応した画素であることを示す。したがって、変換テーブル72の出力β（Ｄmor ）が大きな値を出力するのは、その画素が石灰化陰影を構成する画素の場合である。
【００８４】
また、この変換テーブル72によれば、Ｄmor の小さい領域Ｃ１においてその出力β（Ｄmor ）は低く抑制されるため、Ｄmor の小さい領域Ｃ１において検出される高周波の放射線ノイズを低減させることができる。すなわち、図２（２）に示した濃度値の分布曲線上には実際には、同図（３）の拡大図に示すように、放射線ノイズが重畳しているが、式（３）の第２項のクロージング処理により、この放射線ノイズは同図（３）の破線で示すように平滑化される。したがってＤmor の値を非常に細かく変動させることとなるが、微小石灰化陰影による変動Ｄmor に比べてその変動量自体は小さいため、領域Ｃ１とＣ２との境界値を予め適当に設定することにより、この高周波ノイズが強調されるのを防止することができる。
【００８５】
一方、領域Ｃ３においては、その出力β（Ｄmor ）の変化がＤmor の変化に対して抑制されるが、これは、既にある程度のコントラストを有する画像部分について過度の強調処理がなされるのを防止するためである。過度の強調処理がなされると、その画像部分以外の画像部分のコントラストが相対的に低下するため、反って読影性能が低下するからである。
【００８６】
石灰化陰影強調手段73は、変換テーブル72からの出力β（Ｄmor ）と高周波成分を算出する手段71からの出力である高周波成分（Ｄorg −Ｄus）との積β（Ｄmor ）×（Ｄorg −Ｄus）を算出し、その積に原画像の濃度値Ｄorg を加算した式（１）で示す周波数強調処理を施す（本実施例において、式（１）のｆ（Ｄorg −Ｄus）をｆ（Ｄorg −Ｄus）＝Ｄorg −Ｄus、とする）。
【００８７】
この周波数強調処理によれば、モーフォロジーフィルター40により得られた石灰化陰影を構成する画素か否かに応じた信号Ｄmor に基づく強調係数β（Ｄmor ）で、高周波成分（Ｄorg −Ｄus）を強調するため、この高周波成分（Ｄorg −Ｄus）に例えば量子ノイズ等が含まれていても、その画素が石灰化陰影等の画像部分を構成するものでない場合（例えば、血管等の陰影の場合）は、その画素についてのβ（Ｄmor ）の値は小さいため、強調度合は小さいものとなる。一方、画素が石灰化陰影等の画像部分を構成するものである場合は、その画素についてのβ（Ｄmor ）の値は大きいため、その強調度合は大きいものとなる。
【００８８】
したがって、画像の高周波成分（Ｄorg −Ｄus）に放射線ノイズが含まれているか否かに拘らず、石灰化陰影等の特定の画像部分であるか否かに応じた関数β（Ｄmor ）により、特定の画像部分を選択的に強調処理することができる。
【００８９】
局所画像表示手段90は、この局所画像強調手段70により局所領域の画像Ｐ2 のうち石灰化陰影Ｐ1 が強調処理された画像を表示する。
【００９０】
このように全体画像のうち、石灰化陰影Ｐ1 の画像だけが別個に局所画像表示手段90に表示されるため、読影者はその表示された局所領域の画像に観察意識、診断意識を集中することができ、診断性能を向上させることができる。
【００９１】
また、全体画像表示手段30が局所画像表示手段90を兼ねる構成であってもよく、その場合であっても、表示された全体画像Ｐのうち石灰化陰影Ｐ1 の画像だけが選択的に強調されるため、オーバーシュートやアンダーシュートが抑制され、これによるアーチファクトが低減されて診断性能の向上した再生画像を得ることができる。
【００９２】
なお本実施例の画像処理方法は、式（３）で示したモーフォロジー演算についての例について説明したが、本発明の画像処理方法はこれに限らず、式（２）、（４）〜（８）に示したうちのいずれのモーフォロジー演算を用いた強調処理であってもよい。
【００９３】
ところで、モーフォロジーフィルター40の式（２）〜（８）による微小石灰化陰影の検出処理でだけでは、微小石灰化陰影と類似の陰影も検出される場合がある。すなわち、式（２）〜（８）のＤmor の値はゼロにならないような、石灰化陰影とほぼ同一の大きさの、石灰化陰影ではない画像（以下、非石灰化陰影という）も検出される場合がある。このような非石灰化陰影について強調処理をしたのでは、適切な診断に支障をきたす虞がある。
【００９４】
そこでこのような非石灰化陰影が石灰化陰影とともに検出されるのを防止して、より正確に石灰化陰影だけを検出するために、モーフォロジーフィルター40に以下に示す判定機能を付加してもよい。
【００９５】
すなわち、式（19）に示すモーフォロジー演算に基づく微分演算を行い、Ｍgradの値が大きいほど石灰化陰影である可能性が大きいため、その後に、式（２）〜（８）に代えて、式（20）に示す論理演算を行う。
【００９６】
【数１５】
0003731760.tif 000016
【００９７】
この式（20）で得られたＤmor の値がゼロの場合は、その画素は非石灰化陰影であるから、式（１）による強調処理はなされず、Ｄmor の値がゼロ以外の場合は、その画素は石灰化陰影であるから、式（１）による強調処理がなされる。なお、式（20）中のＴ１，Ｔ２は実験的に設定された閾値である。
【００９８】
また、この式（19）および（20）による判定の他に、マルチスケールのオープニング処理とクローシング処理との組合せによって、石灰化陰影と非石灰化陰影との判別の判定をすることもできる。
【００９９】
具体的には、下式（21）〜（23）によってＤmor の値を設定すればよい。
【０１００】
【数１６】
0003731760.tif 000017
【０１０１】
なお、式（23）中のＴは実験的に設定された閾値である。
【図面の簡単な説明】
【図１】計算機支援画像診断装置の概略構成を示すブロック図
【図２】図１に示した計算機支援画像診断装置により画像診断に供される乳房の放射線画像（マンモグラム）を示す図
【図３】変換テーブルを表す関数のグラフ
【図４】モーフォロジー演算の基本的な作用を説明する図
【図５】モーフォロジーフィルターにおける多重構造要素を示す図
【図６】スケルトン処理を示す説明図
【図７】オーバーシュート、アンダーシュートを抑制する関数ｆ（Ｄorg −Ｄus）を表すグラフ
【図８】モーフォロジー演算による処理を具体的に説明するための濃度分布図
【図９】従来のボケマスク処理と、本発明の画像処理とを比較した図
【符号の説明】
10    全体画像メモリ
30    全体画像表示手段
40    異常陰影検出手段
50    判定手段
60    局所領域抽出手段
70    局所画像強調手段
71    高周波成分算出手段
72    変換テーブル
73    石灰化陰影強調手段
90    局所画像表示手段
Ｄorg   オリジナル濃度値
【図１】
0003731760.tif 000026
【図２】
0003731760.tif 000027
【図３】
0003731760.tif 000028
【図４】
0003731760.tif 000029
【図５】
0003731760.tif 000030
【図６】
0003731760.tif 000031
【図７】
0003731760.tif 000032
【図８】
0003731760.tif 000033
【図９】
0003731760.tif 000034