JP3932470B2 - 画像処理方法および装置 - Google Patents
画像処理方法および装置 Download PDFInfo
- Publication number
- JP3932470B2 JP3932470B2 JP05481998A JP5481998A JP3932470B2 JP 3932470 B2 JP3932470 B2 JP 3932470B2 JP 05481998 A JP05481998 A JP 05481998A JP 5481998 A JP5481998 A JP 5481998A JP 3932470 B2 JP3932470 B2 JP 3932470B2
- Authority
- JP
- Japan
- Prior art keywords
- processing
- image signal
- image
- closing
- opening
- 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 - Lifetime
Links
- 238000003672 processing method Methods 0.000 title claims description 27
- 238000000034 method Methods 0.000 claims description 117
- 230000000877 morphologic effect Effects 0.000 claims description 35
- 238000010586 diagram Methods 0.000 description 6
- 230000010339 dilation Effects 0.000 description 5
- 230000003628 erosive effect Effects 0.000 description 5
- 230000001629 suppression Effects 0.000 description 4
- 201000005569 Gout Diseases 0.000 description 2
- 239000003086 colorant Substances 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 102100035971 Molybdopterin molybdenumtransferase Human genes 0.000 description 1
- 101710119577 Molybdopterin molybdenumtransferase Proteins 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は画像処理方法および画像処理装置に関し、詳細には、画像の粒状（高周波ノイズ）を抑制する処理方法および装置に関するものである。
【０００２】
【従来の技術】
従来より、写真フイルムやプリントに記録された画像（カラー画像を含む）をＣＣＤ等のセンサーにより光電的に読み取って画像信号を得、これに種々の画像処理を施して、処理済画像をプリントやＣＲＴ等に再生することが行われている。このような画像処理としては、画像に含まれる所定の空間周波数成分に作用させる周波数処理や、画像の濃度などに作用する階調処理等がある。
【０００３】
さらに周波数処理としては、画像の輪郭のぼけを抑制する鮮鋭度強調処理や、感材の粒状度等に起因するノイズ（粒状）を抑制する平滑化処理等の粒状抑制処理などがある。
【０００４】
ここで粒状抑制処理の一つの手法として、近年、モルフォロジー（Morphology）演算処理が研究されている。
【０００５】
このモルフォロジー演算処理は、モフォロジーまたはモーフォロジーとも称し、一般的にはＮ次元空間における集合論として展開されるが、２次元空間である画像へ適用される場合が多い（特開平8-272961号、同9-248291号、同9-91421 号等）。ここでは、このモルフォロジー演算処理を、濃淡画像を例にして以下に簡単に説明する。
【０００６】
濃淡画像を座標（ｘ，ｙ）の点が濃度値ｆ（ｘ，ｙ）に相当する高さをもつ空間とみなし、この断面に相当する１次元の関数ｆ（ｘ）を考える。モルフォロジー演算処理に用いる構造要素ｇは次式（３）に示すように、原点について対称な対称関数
【０００７】
【数３】
【０００８】
であり、定義域内で値が０で、その定義域Ｇが下記式（４）であるとする。
【０００９】
【数４】
【００１０】
このとき、モルフォロジー演算の基本形は式（５）〜（８）に示すように、非常に簡単な演算となる。
【００１１】
【数５】
【００１２】
すなわち、ダイレーション（dilation）処理は、注目画素を中心とした、±ｍ（構造要素Ｂに応じて決定される値であって、図４中のマスクサイズに相当）の幅の範囲内の最大値を探索する処理であり（同図（Ａ）参照）、一方、イロージョン（erosion ）処理は、注目画素を中心とした、±ｍの幅の範囲内の最小値を探索する処理である（同図（Ｂ）参照）。また、オープニング（opening ）処理はイロージョン処理後にダイレーション処理を行なう処理、すなわち最小値の探索の後に最大値を探索する処理であり、クロージング（closing ）処理は、ダイレーション処理後にイロージョン処理を行なう処理、すなわち最大値の探索の後に最小値を探索する処理に相当する。
【００１３】
つまりオープニング処理は、低濃度側から濃度曲線ｆ（ｘ）を滑らかにし、マスクサイズ２ｍより空間的に狭い範囲で変動する凸状の濃度変動部分（周囲部分よりも濃度が高い部分）を抑制することに相当する（同図（Ｃ）参照）。
【００１４】
一方、クロージング処理は、高濃度側から濃度曲線ｆ（ｘ）を滑らかにし、マスクサイズ２ｍより空間的に狭い範囲で変動する凹状の濃度変動部分（周囲部分よりも濃度が低い部分）を抑制することに相当する（同図（Ｄ）参照）。
【００１５】
ここで、濃度の高いもの程大きな値となる高濃度高信号レベルの信号の場合においては、濃度値ｆ（ｘ）の画像信号値が高輝度高信号レベルの場合に対して大小関係が逆転するため、高濃度高信号レベルの信号に対するダイレーション処理と高輝度高信号レベルに対するイロージョン処理（同図（Ｂ））とは一致し、高濃度高信号レベルの信号に対するイロージョン処理と高輝度高信号レベルに対するダイレーション処理（同図（Ａ））とは一致し、高濃度高信号レベルの信号に対するオープニング処理と高輝度高信号レベルに対するクロージング処理（同図（Ｄ））とは一致し、高濃度高信号レベルの信号に対するクロージング処理と高輝度高信号レベルに対するオープニング処理（同図（Ｃ））とは一致する。
【００１６】
そして、このように原画像を表す画像信号に対して、モルフォロジー演算処理によるオープニング処理若しくはクロージング処理を施すことにより、画像から高周波ノイズ（パルス性雑音）を抑制（または除去）することができるが、上述したように、オープニング処理によれば、凸状の濃度変動部分（以下、本明細書中「正のノイズ」又は「正の粒状」という）のみ、クロージング処理によれば、凹状の濃度変動部分（以下、本明細書中「負のノイズ」又は「負の粒状」という）のみが抑制される。
【００１７】
そこで、オープニング処理により正のノイズを抑制した後に、そのオープニング処理を施して得られたオープニング処理画像信号にさらにクロージング処理を施し、またはクロージング処理により負のノイズを抑制した後に、そのクロージング処理を施して得られたクロージング処理画像信号にさらにオープニング処理を施すことで、正負両方のノイズを抑制することが提案されている（小畑「モルフォロジー」（コロナ社刊）P93,P104等）。
【００１８】
また、同一の画像信号から各別にクロージング処理画像信号とオープニング処理画像信号とを求め、これら得られたクロージング処理画像信号とオープニング処理画像信号との平均値を用いることによってノイズを抑制することも提案されている（「Detail-preserving smoothing with morphology and fuzzy reasoning 」（Journal of Electronic Imaging(Jul.1996) vol.5(3),P396-P401）。
【００１９】
さらに、上記文献にも紹介されているTsukasamoto's Fuzzy Modelsと称される方法では、クロージング処理画像信号とオープニング処理画像信号との差を、あるテーブルと比較し、ある閾値以上であれば当該部は粒状であると判定し、閾値以下であれば粒状ではなく画像の微細構造と判定することで、粒状と判定した部分を処理するものである。
【００２０】
【発明が解決しようとする課題】
しかし、上記オープニング処理およびクロージング処理を順次に行う方法は、後に行う側の処理（オープニング処理後にクロージング処理を行う場合はクロージング処理を意味し、クロージング処理後にオープニング処理を行う場合はオープニング処理を意味する）の対象画像が、原画像ではなくオープニング処理後の画像（正のノイズが抑制された画像）またはクロージング処理後の画像（負のノイズが抑制された画像）であるため、このような最初の処理後の画像にクロージング処理またはオープニング処理を施しても、元の画像（原画像）に残存する負のノイズまたは正のノイズを正確に抑制する処理とはいえない。すなわち、既にオープニング処理が施されたことによって、本来はクロージング処理で抑制されるべき画像中の負のノイズ部分を抑制できなくなり、またはこれとは反対に、既にオープニング処理が施されたことによって、本来はクロージング処理で抑制されるべきでない画像部分が負のノイズとして抑制されてしまう場合等がある。
【００２１】
また、全画素を単に、クロージング処理画像信号とオープニング処理画像信号との平均値とする処理方法（上記「Detail-preserving ……reasoning 」に記載された処理方法）は、視覚的に得られる粒状の抑制程度が小さい。
【００２２】
さらにまた、Tsukasamoto's Fuzzy Modelsは、参照テーブルを、画像ごとに準備する必要があり、しかも同一画像であっても画像上の粒状の大小に応じて相異なる参照テーブルを適用すべき点において、準備すべき多数のテーブルを決定するのは煩雑で、実用性において問題がある。
【００２３】
本発明は上記事情に鑑みなされたものであって、正および負のノイズを精度よく抑制する画像処理方法および画像処理装置を提供することを目的とするものである。
【００２４】
【課題を解決するための手段】
本発明の画像処理方法は、画像信号に対して、モルフォロジー演算によるオープニング処理とクロージング処理とを各別に施して、得られた各処理後の画像信号同士を、画素を対応させて演算処理して１つの処理済画像信号を得るものである。
【００２５】
すなわち本発明の第１の画像処理方法は、画像を構成する画素の値を規定する画像信号（輝度や色差の信号、カラー画像に基づく画像信号等も含む）ｆinに対して、モルフォロジー演算によるオープニング処理およびクロージング処理を各別に施し、
前記オープニング処理により取得されたオープニング処理画像信号ｆopn および前記クロージング処理により取得されたクロージング処理画像信号ｆcls に基づいて、該両処理画像信号の画素を対応させた下記式（１）で定義される演算処理ｇ（ｆopn ，ｆcls ）（＝fout ）により処理済画像信号fout を取得することを特徴とするものである。
【００２６】
【数１】
【００２７】
すなわち、クロージング処理画像信号ｆcls が原画像信号ｆinに等しい画素についてはオープニング処理画像信号ｆopn を処理済画像信号ｆout とし、オープニング処理画像信号ｆopn が原画像信号ｆinに等しい画素については、クロージング処理画像信号ｆcls を処理済画像信号ｆout とし、クロージング処理画像信号ｆcls またはオープニング処理画像信号ｆopn のいずれも原画像信号ｆinに等しくない画素については、クロージング処理画像信号ｆcls とオープニング処理画像信号ｆopn とを、重み付け係数w1，w2を用いた重み付け加算処理を適用するものである。
【００２８】
ここで重み付け係数w1，w2をどのように決定するかについては種々の決定方法を適用することができるが、各重み付け係数w1、w2は好ましくは 0.2から 0.8までの値、より好ましくは 0.3から 0.7までの値、さらに好ましくは 0.4から 0.6までの値、最適にはいずれも 0.5を採用するのが望ましい。
【００２９】
また、本発明の第２の画像処理方法は、画像を構成する画素の値を規定する画像信号（輝度や色差の信号、カラー画像に基づく画像信号等も含む）ｆinに対して、モルフォロジー演算によるオープニング処理およびクロージング処理を各別に施し、
前記オープニング処理により取得されたオープニング処理画像信号ｆopn および前記クロージング処理により取得されたクロージング処理画像信号ｆcls に基づいて、該両処理画像信号の画素を対応させた下記式（２）で定義される演算処理ｇ（ｆopn ，ｆcls ）（＝fout ）により処理済画像信号fout を取得することを特徴とするものである。
【００３０】
【数２】
【００３１】
すなわち、オープニング処理画像信号ｆopn が原画像信号ｆinに等しい画素については、クロージング処理画像信号ｆcls を処理済画像信号ｆout とし、クロージング処理画像信号ｆcls が原画像信号ｆinに等しい画素については、オープニング処理画像信号ｆopn を処理済画像信号ｆout とし、オープニング処理画像信号ｆopn またはクロージング処理画像信号ｆcls のいずれも原画像信号ｆinに等しくない画素については、オープニング処理画像信号ｆopn とクロージング処理画像信号ｆcls とのうち原画像信号ｆinとの差が大きい方の処理画像信号（ｆopn またはｆcls ）と、原画像信号ｆinとの差（｜ｆopn −ｆin｜または｜ｆcls −ｆin｜）、による重みづけ加算した値を処理済画像信号ｆout とする演算処理を適用するものである。
【００３２】
本発明の第１の画像処理装置は、上記本発明の第１の画像処理方法を実施するための装置であって、画像を構成する画素の値を規定する画像信号に対して、モルフォロジー演算によるオープニング処理およびクロージング処理を各別に施すモルフォロジー演算処理手段と、
前記モルフォロジー演算処理手段による前記オープニング処理により取得されたオープニング処理画像信号および前記クロージング処理により取得されたクロージング処理画像信号に基づいて、該両処理画像信号の画素を対応させた上記式（１）で定義される演算処理により処理済画像信号を取得する信号処理手段とを備えたことを特徴とするものである。
【００３３】
なお重み付け係数w1，w2の決定方法および適切な係数値については、上述した本発明の第１の画像処理方法において説明したものと同様である。
【００３４】
本発明の第２の画像処理装置は、上記本発明の第２の画像処理方法を実施するための装置であって、画像を構成する画素の値を規定する画像信号に対して、モルフォロジー演算によるオープニング処理およびクロージング処理を各別に施すモルフォロジー演算処理手段と、
前記モルフォロジー演算処理手段による前記オープニング処理により取得されたオープニング処理画像信号および前記クロージング処理により取得されたクロージング処理画像信号に基づいて、該両処理画像信号の画素を対応させた上記式（２）で定義される演算処理により処理済画像信号を取得する信号処理手段とを備えたことを特徴とするものである。
【００３５】
上記本発明の画像処理方法および画像処理装置における、オープニング処理の構造要素とクロージング処理の構造要素とは、同一の大きさおよび形状とするのが、処理の簡易化のうえで望ましい。
【００３６】
【発明の効果】
本発明の画像処理方法および装置は、画像信号に対して、モルフォロジー演算によるオープニング処理とクロージング処理とを各別に施して、オープニング処理により正のノイズが抑制された画像信号とクロージング処理により負のノイズが抑制された画像信号とをそれぞれ得、これらの得られた両画像信号を、画素を対応させて演算処理して１つの処理済画像信号を得るため、原画像に存在する正および負のノイズを精度よく抑制することができる。
【００３７】
また、この演算処理の内容を全画素について共通のものとせずに、画素ごとに異なるものとすることにより、画像中の構造物の内容等に応じて、適切にノイズを抑制することができるため、全ての画素についてオープニング処理画像信号とクロージング処理画像信号との平均値を適用する従来の処理方法（上述した「Detail-preserving ……reasoning 」に記載の方法）に比して、的確に粒状を抑制することができ、視覚的な粒状の抑制程度を向上させることができる。また、多数の参照テーブルを用いる従来の処理方法に比べて簡易的な手法で粒状を抑制することができる。
【００３８】
ここで演算処理式（１）に示すものを適用する本発明の第１の画像処理方法および装置によれば、オープニング処理により抑制されない負のノイズを表す画素（オープニング処理画像信号が原画像信号に等しい画素）については、クロージング処理画像信号を採用することにより、この負のノイズが抑制された処理済画像信号となり、クロージング処理により抑制されない正のノイズを表す画素（クロージング処理画像信号が原画像信号に等しい画素）については、オープニング処理画像信号を採用することにより、この正のノイズが抑制された処理済画像信号となり、正のノイズと負のノイズとの両方を抑制することができる。さらに、オープニング処理画像信号またはクロージング処理画像信号のいずれも原画像信号に等しくない画素については、オープニング処理画像信号とクロージング処理画像信号との加重平均値を処理済画像信号とすることにより、良好に正負のノイズを抑制することができる。このとき特に加重平均の加重割合を５対５とすれば、最適に粒状を抑制することができる。
【００３９】
また、演算処理式（２）に示すものを適用する本発明の第２の画像処理方法および装置によれば、オープニング処理により抑制されない負のノイズを表す画素（オープニング処理画像信号が原画像信号に等しい画素）については、クロージング処理画像信号を採用することにより、この負のノイズが抑制された処理済画像信号となり、クロージング処理により抑制されない正のノイズを表す画素（クロージング処理画像信号が原画像信号に等しい画素）については、オープニング処理画像信号を採用することにより、この正のノイズが抑制された処理済画像信号となり、正のノイズと負のノイズとの両方を抑制することができる。さらに、オープニング処理画像信号またはクロージング処理画像信号のいずれも原画像信号に等しくない画素については、オープニング処理画像信号ｆopn とクロージング処理画像信号ｆcls とのうち原画像信号ｆinとの差が大きい方の処理画像信号と、原画像信号ｆinとの差、による重みづけ加算した値を処理済画像信号ｆout とする処理済画像信号とすることにより、良好に正負のノイズを抑制することができる。
【００４０】
【発明の実施の形態】
以下、本発明の画像処理方法および画像処理装置の具体的な実施の形態について、図面を用いて説明する。
【００４１】
図１は本発明の画像処理方法の一実施形態の処理フローを示すフローチャート、図２は図１に示した画像処理方法を実施する画像処理装置を示す図である。
【００４２】
図示の画像処理装置は、デジタルスチルカメラにより撮影された、またはカラー写真プリントもしくはカラーフイルムから光電的に読み取られた、カラー画像を表すＲＧＢの画像信号Ｒin，Ｇin，Ｂinの入力を受けて、これらのＲＧＢ画像信号Ｒin，Ｇin，Ｂinに対して各別に、モルフォロジー演算によるオープニング処理およびクロージング処理をそれぞれ施すモルフォロジー演算処理手段10と、モルフォロジー演算処理手段10によるオープニング処理により取得されたオープニング処理画像信号Ｒopn ，Ｇopn ，Ｂopn およびクロージング処理により取得されたクロージング処理画像信号Ｒcls ，Ｇcls ，Ｂcls に基づいて、両処理画像信号の画素を対応させた画素ごとの演算処理ｇ（Ｒopn ，Ｒcls ），ｇ（Ｇopn ，Ｇcls ），ｇ（Ｂopn ，Ｂcls ）により処理済画像信号Ｒout ，Ｇout ，Ｂout を取得する信号処理手段20とを備えた構成である。
【００４３】
ここで、信号処理手段20による各演算処理ｇとしては例えば、下記式（１）に示す演算処理が適用される。
【００４４】
【数１】
【００４５】
なお、式中の記号ｆは信号値の一般的な表示であり、本実施形態のカラー画像信号を適用するにあたっては、この一般的な表示である記号ｆをそれぞれ各画像信号Ｒ、Ｇ、Ｂに置き換えて適用するものである。
【００４６】
この演算処理は、クロージング処理画像信号Ｒcls が原画像信号Ｒinに等しい画素については、オープニング処理画像信号Ｒopn を処理済画像信号Ｒout とし、オープニング処理画像信号Ｒopn が原画像信号Ｒinに等しい画素については、クロージング処理画像信号Ｒcls を処理済画像信号Ｒout とし、オープニング処理画像信号Ｒopn またはクロージング処理画像信号Ｒcls のいずれも原画像信号Ｒinに等しくない画素については、オープニング処理画像信号Ｒopn とクロージング処理画像信号Ｒcls とを、重み付け係数w1，w2を用いた重み付け加算処理を適用する演算処理である。他の色Ｇ，Ｂの画像信号についても同様である。
【００４７】
ここで重み付け係数w1，w2は、いずれも 0.5を適用するものとする。
【００４８】
モルフォロジー演算処理手段10による上記オープニング処理およびクロージング処理に適用される構造要素（マスク）は、図５に示す４つの直線構造要素からなる多重構造要素（図示のものは５画素×５画素のマスクサイズとなる）である。多重構造要素を用いたオープニング処理およびクロージング処理は、多重構造要素を構成する各構造要素を適用したオープニング処理およびクロージング処理をそれぞれ各別に施し、オープニング処理についてはそれらの各別に行ったオープニング処理のうち最大値を、クロージング処理についてはそれらの各別に行ったクロージング処理のうち最小値を、それぞれ採る処理であり、画像のように２次元状に分布する信号値を処理する場合に有用である。なお、マスクサイズ（構造要素の大きさ）は、削除しようとする粒状の大きさに応じて変更することができる。一般的な画像においては、３画素×３画素から13画素×13画素程度の大きさとなる多重構造要素を適用するのが好適であり、３画素×３画素から７画素×７画素程度の大きさとなる多重構造要素を適用するのが最適である。
【００４９】
ただし、以下においては、説明の簡単化のため、２次元画像の一断面である１次元の信号分布を例にして説明するため、多重構造要素のうち１つの構造要素を用いた説明としている。
【００５０】
次に、本実施形態の画像処理装置の作用について、原画像信号の１つである図３（１）に示すＲ画像信号Ｒinを例に説明する。
【００５１】
まず、デジタルスチルカメラ等から、カラー画像を表すＲＧＢの画像信号Ｒin，Ｇin，Ｂinがモルフォロジー演算手段10に入力される。画像信号Ｒinは上述したように図３（１）に示すような分布を有する信号であり、他の画像信号Ｇin，Ｂinも所定の分布を有している（図示せず）。
【００５２】
原画像信号Ｒinには、比較的大きな周期の濃度分布に、この濃度分布より高周波のノイズが重畳している。この高周波ノイズ（同図（１）のＡ部）を拡大したものを同図（２）に示す。
【００５３】
モルフォロジー演算手段10は、入力された原画像信号Ｒin，Ｇin，Ｂinに対して各別に、予め設定された構造要素を用いたモルフォロジー演算によるオープニング処理およびクロージング処理をそれぞれ施し、オープニング処理画像信号Ｒopn ，Ｇopn ，Ｂopn およびクロージング処理画像信号Ｒcls ，Ｇcls ，Ｂcls を算出する。図３（３）は、原画像信号Ｒinに対するオープニング処理によるオープニング処理画像信号Ｒopn の、図３（２）に相当する部分（同図（１）のＡ部）を示す図であり、同図（４）は、原画像信号Ｒinに対するクロージング処理によるクロージング処理画像信号Ｒcls の、図３（２）に相当する部分（同図（１）のＡ部）を示す図である。
【００５４】
各図から解されるように、オープニング処理画像信号Ｒopn は、原画像信号Ｒin（図３（３）中において破線で示す）のうち、凸状の濃度変動部分（正のノイズ）が抑制された信号、クロージング処理画像信号Ｒcls は、原画像信号Ｒin（図３（４）中において破線で示す）のうち、凹状の濃度変動部分（負のノイズ）が抑制された信号とされる。
【００５５】
他のオープニング処理画像信号Ｇopn ，Ｂopn もオープニング処理画像信号Ｒopn と同様に、正のノイズが抑制された信号とされ、他のクロージング処理画像信号Ｇcls ，Ｂcls もクロージング処理画像信号Ｒcls と同様に、負のノイズが抑制された信号とされる。
【００５６】
モルフォロジー演算手段10により算出された各オープニング処理画像信号Ｒopn ，Ｇopn ，Ｂopn および各クロージング処理画像信号Ｒcls ，Ｇcls ，Ｂcls は、信号処理手段20に入力される。
【００５７】
信号処理手段20は、入力された各画像信号について、各色ＲＧＢごとに、画素を対応させて、オープニング処理画像信号とクロージング処理画像信号とを上記式（１）に従った演算処理を行い、各色ごとに処理済画像信号を算出する。
【００５８】
図３（５）は、同図（２）に示したオープニング処理画像信号Ｒopn と、同図（４）に示したクロージング処理画像信号Ｒcls とを式（１）に従って演算処理して得られた処理済画像信号Ｒout を示す。
【００５９】
図から解されるように、処理済画像信号Ｒout は、原画像信号Ｒin（図３（５）中において破線で示す）のうち、正および負のいずれのノイズも精度よく抑制された信号とされる。なお、重みづけ係数w1，w2の値はいずれも0.5 としたが、この値に限るものではなく、適当な範囲で変動させることもできる。
【００６０】
他の処理済画像信号Ｇout ，Ｂout も処理済画像信号Ｒout と同様に、正および負のいずれのノイズも精度よく抑制された信号とされる。
【００６１】
このように本実施形態の画像処理装置によれば、入力された画像信号に対して、モルフォロジー演算によるオープニング処理とクロージング処理とを各別に施して、オープニング処理により正のノイズが抑制された画像信号とクロージング処理により負のノイズが抑制された画像信号とをそれぞれ得、これらの得られた両画像信号を、画素を対応させて演算処理して１つの処理済画像信号を得るため、原画像に存在する正および負のノイズを精度よく抑制することができる。
【００６２】
なお、上記画像処理装置は、信号処理手段20による各演算処理ｇとして式（１）に示す演算処理を適用した実施形態であるが、本発明の画像処理方法および画像処理装置は、式（１）による演算処理のものに限るものではなく、種々の適切な演算処理を適用することができる。
【００６３】
例えば、下記式（２）に示す演算処理を適用することもできる。
【００６４】
【数２】
【００６５】
すなわち、オープニング処理画像信号Ｒopn が原画像信号Ｒinに等しい画素については、クロージング処理画像信号Ｒcls を処理済画像信号Ｒout とし、クロージング処理画像信号Ｒcls が原画像信号Ｒinに等しい画素については、オープニング処理画像信号Ｒopn を処理済画像信号Ｒout とし、オープニング処理画像信号Ｒopn またはクロージング処理画像信号Ｒcls のいずれも原画像信号Ｒinに等しくない画素については、オープニング処理画像信号Ｒopn とクロージング処理画像信号Ｒcls とのうち原画像信号Ｒinとの差が大きい方の処理画像信号（Ｒopn またはＲcls ）と、原画像信号Ｒinとの差（｜Ｒopn −Ｒin｜または｜Ｒcls −Ｒin｜）、による重みづけ加算した値を処理済画像信号Ｒout とする演算処理である。他の色Ｇ，Ｂの画像信号についても同様である。
【００６６】
このような演算処理を適用したものについても、式（１）に示した演算処理の実施形態と同様に、良好に正負のノイズを抑制することができる。
【００６７】
また本実施形態の画像処理装置においては、カラー画像を表すＲＧＢの画像信号を処理の対象としたが、本発明の画像処理装置はそのようなカラー画像信号を対象とするものに限るものではなく、例えば輝度や色差等の各種信号を処理の対象とすることができるものである。
【図面の簡単な説明】
【図１】本発明の画像処理方法の一実施形態の処理フローを示すフローチャート
【図２】図１に示した画像処理方法を実施する画像処理装置を示す図
【図３】（１）原画像信号Ｒinを示す図、（２）（１）に示した原画像信号Ｒinのうち高周波ノイズ（Ａ部）を示す拡大図、（３）原画像信号Ｒinに基づくオープニング処理画像信号Ｒopn のうちＡ部に相当する部分を示す図、（４）原画像信号Ｒinに基づくクロージング処理画像信号Ｒcls のうちＡ部に相当する部分を示す図、（５）処理済画像信号Ｒout のうちＡ部に相当する部分を示す図
【図４】モルフォロジー演算の基本的な作用を説明する図
【図５】多重構造要素を示す図
【符号の説明】
10 モルフォロジー演算処理手段
20 信号処理手段
【発明の属する技術分野】
本発明は画像処理方法および画像処理装置に関し、詳細には、画像の粒状（高周波ノイズ）を抑制する処理方法および装置に関するものである。
【０００２】
【従来の技術】
従来より、写真フイルムやプリントに記録された画像（カラー画像を含む）をＣＣＤ等のセンサーにより光電的に読み取って画像信号を得、これに種々の画像処理を施して、処理済画像をプリントやＣＲＴ等に再生することが行われている。このような画像処理としては、画像に含まれる所定の空間周波数成分に作用させる周波数処理や、画像の濃度などに作用する階調処理等がある。
【０００３】
さらに周波数処理としては、画像の輪郭のぼけを抑制する鮮鋭度強調処理や、感材の粒状度等に起因するノイズ（粒状）を抑制する平滑化処理等の粒状抑制処理などがある。
【０００４】
ここで粒状抑制処理の一つの手法として、近年、モルフォロジー（Morphology）演算処理が研究されている。
【０００５】
このモルフォロジー演算処理は、モフォロジーまたはモーフォロジーとも称し、一般的にはＮ次元空間における集合論として展開されるが、２次元空間である画像へ適用される場合が多い（特開平8-272961号、同9-248291号、同9-91421 号等）。ここでは、このモルフォロジー演算処理を、濃淡画像を例にして以下に簡単に説明する。
【０００６】
濃淡画像を座標（ｘ，ｙ）の点が濃度値ｆ（ｘ，ｙ）に相当する高さをもつ空間とみなし、この断面に相当する１次元の関数ｆ（ｘ）を考える。モルフォロジー演算処理に用いる構造要素ｇは次式（３）に示すように、原点について対称な対称関数
【０００７】
【数３】
であり、定義域内で値が０で、その定義域Ｇが下記式（４）であるとする。
【０００９】
【数４】
このとき、モルフォロジー演算の基本形は式（５）〜（８）に示すように、非常に簡単な演算となる。
【００１１】
【数５】
すなわち、ダイレーション（dilation）処理は、注目画素を中心とした、±ｍ（構造要素Ｂに応じて決定される値であって、図４中のマスクサイズに相当）の幅の範囲内の最大値を探索する処理であり（同図（Ａ）参照）、一方、イロージョン（erosion ）処理は、注目画素を中心とした、±ｍの幅の範囲内の最小値を探索する処理である（同図（Ｂ）参照）。また、オープニング（opening ）処理はイロージョン処理後にダイレーション処理を行なう処理、すなわち最小値の探索の後に最大値を探索する処理であり、クロージング（closing ）処理は、ダイレーション処理後にイロージョン処理を行なう処理、すなわち最大値の探索の後に最小値を探索する処理に相当する。
【００１３】
つまりオープニング処理は、低濃度側から濃度曲線ｆ（ｘ）を滑らかにし、マスクサイズ２ｍより空間的に狭い範囲で変動する凸状の濃度変動部分（周囲部分よりも濃度が高い部分）を抑制することに相当する（同図（Ｃ）参照）。
【００１４】
一方、クロージング処理は、高濃度側から濃度曲線ｆ（ｘ）を滑らかにし、マスクサイズ２ｍより空間的に狭い範囲で変動する凹状の濃度変動部分（周囲部分よりも濃度が低い部分）を抑制することに相当する（同図（Ｄ）参照）。
【００１５】
ここで、濃度の高いもの程大きな値となる高濃度高信号レベルの信号の場合においては、濃度値ｆ（ｘ）の画像信号値が高輝度高信号レベルの場合に対して大小関係が逆転するため、高濃度高信号レベルの信号に対するダイレーション処理と高輝度高信号レベルに対するイロージョン処理（同図（Ｂ））とは一致し、高濃度高信号レベルの信号に対するイロージョン処理と高輝度高信号レベルに対するダイレーション処理（同図（Ａ））とは一致し、高濃度高信号レベルの信号に対するオープニング処理と高輝度高信号レベルに対するクロージング処理（同図（Ｄ））とは一致し、高濃度高信号レベルの信号に対するクロージング処理と高輝度高信号レベルに対するオープニング処理（同図（Ｃ））とは一致する。
【００１６】
そして、このように原画像を表す画像信号に対して、モルフォロジー演算処理によるオープニング処理若しくはクロージング処理を施すことにより、画像から高周波ノイズ（パルス性雑音）を抑制（または除去）することができるが、上述したように、オープニング処理によれば、凸状の濃度変動部分（以下、本明細書中「正のノイズ」又は「正の粒状」という）のみ、クロージング処理によれば、凹状の濃度変動部分（以下、本明細書中「負のノイズ」又は「負の粒状」という）のみが抑制される。
【００１７】
そこで、オープニング処理により正のノイズを抑制した後に、そのオープニング処理を施して得られたオープニング処理画像信号にさらにクロージング処理を施し、またはクロージング処理により負のノイズを抑制した後に、そのクロージング処理を施して得られたクロージング処理画像信号にさらにオープニング処理を施すことで、正負両方のノイズを抑制することが提案されている（小畑「モルフォロジー」（コロナ社刊）P93,P104等）。
【００１８】
また、同一の画像信号から各別にクロージング処理画像信号とオープニング処理画像信号とを求め、これら得られたクロージング処理画像信号とオープニング処理画像信号との平均値を用いることによってノイズを抑制することも提案されている（「Detail-preserving smoothing with morphology and fuzzy reasoning 」（Journal of Electronic Imaging(Jul.1996) vol.5(3),P396-P401）。
【００１９】
さらに、上記文献にも紹介されているTsukasamoto's Fuzzy Modelsと称される方法では、クロージング処理画像信号とオープニング処理画像信号との差を、あるテーブルと比較し、ある閾値以上であれば当該部は粒状であると判定し、閾値以下であれば粒状ではなく画像の微細構造と判定することで、粒状と判定した部分を処理するものである。
【００２０】
【発明が解決しようとする課題】
しかし、上記オープニング処理およびクロージング処理を順次に行う方法は、後に行う側の処理（オープニング処理後にクロージング処理を行う場合はクロージング処理を意味し、クロージング処理後にオープニング処理を行う場合はオープニング処理を意味する）の対象画像が、原画像ではなくオープニング処理後の画像（正のノイズが抑制された画像）またはクロージング処理後の画像（負のノイズが抑制された画像）であるため、このような最初の処理後の画像にクロージング処理またはオープニング処理を施しても、元の画像（原画像）に残存する負のノイズまたは正のノイズを正確に抑制する処理とはいえない。すなわち、既にオープニング処理が施されたことによって、本来はクロージング処理で抑制されるべき画像中の負のノイズ部分を抑制できなくなり、またはこれとは反対に、既にオープニング処理が施されたことによって、本来はクロージング処理で抑制されるべきでない画像部分が負のノイズとして抑制されてしまう場合等がある。
【００２１】
また、全画素を単に、クロージング処理画像信号とオープニング処理画像信号との平均値とする処理方法（上記「Detail-preserving ……reasoning 」に記載された処理方法）は、視覚的に得られる粒状の抑制程度が小さい。
【００２２】
さらにまた、Tsukasamoto's Fuzzy Modelsは、参照テーブルを、画像ごとに準備する必要があり、しかも同一画像であっても画像上の粒状の大小に応じて相異なる参照テーブルを適用すべき点において、準備すべき多数のテーブルを決定するのは煩雑で、実用性において問題がある。
【００２３】
本発明は上記事情に鑑みなされたものであって、正および負のノイズを精度よく抑制する画像処理方法および画像処理装置を提供することを目的とするものである。
【００２４】
【課題を解決するための手段】
本発明の画像処理方法は、画像信号に対して、モルフォロジー演算によるオープニング処理とクロージング処理とを各別に施して、得られた各処理後の画像信号同士を、画素を対応させて演算処理して１つの処理済画像信号を得るものである。
【００２５】
すなわち本発明の第１の画像処理方法は、画像を構成する画素の値を規定する画像信号（輝度や色差の信号、カラー画像に基づく画像信号等も含む）ｆinに対して、モルフォロジー演算によるオープニング処理およびクロージング処理を各別に施し、
前記オープニング処理により取得されたオープニング処理画像信号ｆopn および前記クロージング処理により取得されたクロージング処理画像信号ｆcls に基づいて、該両処理画像信号の画素を対応させた下記式（１）で定義される演算処理ｇ（ｆopn ，ｆcls ）（＝fout ）により処理済画像信号fout を取得することを特徴とするものである。
【００２６】
【数１】
すなわち、クロージング処理画像信号ｆcls が原画像信号ｆinに等しい画素についてはオープニング処理画像信号ｆopn を処理済画像信号ｆout とし、オープニング処理画像信号ｆopn が原画像信号ｆinに等しい画素については、クロージング処理画像信号ｆcls を処理済画像信号ｆout とし、クロージング処理画像信号ｆcls またはオープニング処理画像信号ｆopn のいずれも原画像信号ｆinに等しくない画素については、クロージング処理画像信号ｆcls とオープニング処理画像信号ｆopn とを、重み付け係数w1，w2を用いた重み付け加算処理を適用するものである。
【００２８】
ここで重み付け係数w1，w2をどのように決定するかについては種々の決定方法を適用することができるが、各重み付け係数w1、w2は好ましくは 0.2から 0.8までの値、より好ましくは 0.3から 0.7までの値、さらに好ましくは 0.4から 0.6までの値、最適にはいずれも 0.5を採用するのが望ましい。
【００２９】
また、本発明の第２の画像処理方法は、画像を構成する画素の値を規定する画像信号（輝度や色差の信号、カラー画像に基づく画像信号等も含む）ｆinに対して、モルフォロジー演算によるオープニング処理およびクロージング処理を各別に施し、
前記オープニング処理により取得されたオープニング処理画像信号ｆopn および前記クロージング処理により取得されたクロージング処理画像信号ｆcls に基づいて、該両処理画像信号の画素を対応させた下記式（２）で定義される演算処理ｇ（ｆopn ，ｆcls ）（＝fout ）により処理済画像信号fout を取得することを特徴とするものである。
【００３０】
【数２】
すなわち、オープニング処理画像信号ｆopn が原画像信号ｆinに等しい画素については、クロージング処理画像信号ｆcls を処理済画像信号ｆout とし、クロージング処理画像信号ｆcls が原画像信号ｆinに等しい画素については、オープニング処理画像信号ｆopn を処理済画像信号ｆout とし、オープニング処理画像信号ｆopn またはクロージング処理画像信号ｆcls のいずれも原画像信号ｆinに等しくない画素については、オープニング処理画像信号ｆopn とクロージング処理画像信号ｆcls とのうち原画像信号ｆinとの差が大きい方の処理画像信号（ｆopn またはｆcls ）と、原画像信号ｆinとの差（｜ｆopn −ｆin｜または｜ｆcls −ｆin｜）、による重みづけ加算した値を処理済画像信号ｆout とする演算処理を適用するものである。
【００３２】
本発明の第１の画像処理装置は、上記本発明の第１の画像処理方法を実施するための装置であって、画像を構成する画素の値を規定する画像信号に対して、モルフォロジー演算によるオープニング処理およびクロージング処理を各別に施すモルフォロジー演算処理手段と、
前記モルフォロジー演算処理手段による前記オープニング処理により取得されたオープニング処理画像信号および前記クロージング処理により取得されたクロージング処理画像信号に基づいて、該両処理画像信号の画素を対応させた上記式（１）で定義される演算処理により処理済画像信号を取得する信号処理手段とを備えたことを特徴とするものである。
【００３３】
なお重み付け係数w1，w2の決定方法および適切な係数値については、上述した本発明の第１の画像処理方法において説明したものと同様である。
【００３４】
本発明の第２の画像処理装置は、上記本発明の第２の画像処理方法を実施するための装置であって、画像を構成する画素の値を規定する画像信号に対して、モルフォロジー演算によるオープニング処理およびクロージング処理を各別に施すモルフォロジー演算処理手段と、
前記モルフォロジー演算処理手段による前記オープニング処理により取得されたオープニング処理画像信号および前記クロージング処理により取得されたクロージング処理画像信号に基づいて、該両処理画像信号の画素を対応させた上記式（２）で定義される演算処理により処理済画像信号を取得する信号処理手段とを備えたことを特徴とするものである。
【００３５】
上記本発明の画像処理方法および画像処理装置における、オープニング処理の構造要素とクロージング処理の構造要素とは、同一の大きさおよび形状とするのが、処理の簡易化のうえで望ましい。
【００３６】
【発明の効果】
本発明の画像処理方法および装置は、画像信号に対して、モルフォロジー演算によるオープニング処理とクロージング処理とを各別に施して、オープニング処理により正のノイズが抑制された画像信号とクロージング処理により負のノイズが抑制された画像信号とをそれぞれ得、これらの得られた両画像信号を、画素を対応させて演算処理して１つの処理済画像信号を得るため、原画像に存在する正および負のノイズを精度よく抑制することができる。
【００３７】
また、この演算処理の内容を全画素について共通のものとせずに、画素ごとに異なるものとすることにより、画像中の構造物の内容等に応じて、適切にノイズを抑制することができるため、全ての画素についてオープニング処理画像信号とクロージング処理画像信号との平均値を適用する従来の処理方法（上述した「Detail-preserving ……reasoning 」に記載の方法）に比して、的確に粒状を抑制することができ、視覚的な粒状の抑制程度を向上させることができる。また、多数の参照テーブルを用いる従来の処理方法に比べて簡易的な手法で粒状を抑制することができる。
【００３８】
ここで演算処理式（１）に示すものを適用する本発明の第１の画像処理方法および装置によれば、オープニング処理により抑制されない負のノイズを表す画素（オープニング処理画像信号が原画像信号に等しい画素）については、クロージング処理画像信号を採用することにより、この負のノイズが抑制された処理済画像信号となり、クロージング処理により抑制されない正のノイズを表す画素（クロージング処理画像信号が原画像信号に等しい画素）については、オープニング処理画像信号を採用することにより、この正のノイズが抑制された処理済画像信号となり、正のノイズと負のノイズとの両方を抑制することができる。さらに、オープニング処理画像信号またはクロージング処理画像信号のいずれも原画像信号に等しくない画素については、オープニング処理画像信号とクロージング処理画像信号との加重平均値を処理済画像信号とすることにより、良好に正負のノイズを抑制することができる。このとき特に加重平均の加重割合を５対５とすれば、最適に粒状を抑制することができる。
【００３９】
また、演算処理式（２）に示すものを適用する本発明の第２の画像処理方法および装置によれば、オープニング処理により抑制されない負のノイズを表す画素（オープニング処理画像信号が原画像信号に等しい画素）については、クロージング処理画像信号を採用することにより、この負のノイズが抑制された処理済画像信号となり、クロージング処理により抑制されない正のノイズを表す画素（クロージング処理画像信号が原画像信号に等しい画素）については、オープニング処理画像信号を採用することにより、この正のノイズが抑制された処理済画像信号となり、正のノイズと負のノイズとの両方を抑制することができる。さらに、オープニング処理画像信号またはクロージング処理画像信号のいずれも原画像信号に等しくない画素については、オープニング処理画像信号ｆopn とクロージング処理画像信号ｆcls とのうち原画像信号ｆinとの差が大きい方の処理画像信号と、原画像信号ｆinとの差、による重みづけ加算した値を処理済画像信号ｆout とする処理済画像信号とすることにより、良好に正負のノイズを抑制することができる。
【００４０】
【発明の実施の形態】
以下、本発明の画像処理方法および画像処理装置の具体的な実施の形態について、図面を用いて説明する。
【００４１】
図１は本発明の画像処理方法の一実施形態の処理フローを示すフローチャート、図２は図１に示した画像処理方法を実施する画像処理装置を示す図である。
【００４２】
図示の画像処理装置は、デジタルスチルカメラにより撮影された、またはカラー写真プリントもしくはカラーフイルムから光電的に読み取られた、カラー画像を表すＲＧＢの画像信号Ｒin，Ｇin，Ｂinの入力を受けて、これらのＲＧＢ画像信号Ｒin，Ｇin，Ｂinに対して各別に、モルフォロジー演算によるオープニング処理およびクロージング処理をそれぞれ施すモルフォロジー演算処理手段10と、モルフォロジー演算処理手段10によるオープニング処理により取得されたオープニング処理画像信号Ｒopn ，Ｇopn ，Ｂopn およびクロージング処理により取得されたクロージング処理画像信号Ｒcls ，Ｇcls ，Ｂcls に基づいて、両処理画像信号の画素を対応させた画素ごとの演算処理ｇ（Ｒopn ，Ｒcls ），ｇ（Ｇopn ，Ｇcls ），ｇ（Ｂopn ，Ｂcls ）により処理済画像信号Ｒout ，Ｇout ，Ｂout を取得する信号処理手段20とを備えた構成である。
【００４３】
ここで、信号処理手段20による各演算処理ｇとしては例えば、下記式（１）に示す演算処理が適用される。
【００４４】
【数１】
なお、式中の記号ｆは信号値の一般的な表示であり、本実施形態のカラー画像信号を適用するにあたっては、この一般的な表示である記号ｆをそれぞれ各画像信号Ｒ、Ｇ、Ｂに置き換えて適用するものである。
【００４６】
この演算処理は、クロージング処理画像信号Ｒcls が原画像信号Ｒinに等しい画素については、オープニング処理画像信号Ｒopn を処理済画像信号Ｒout とし、オープニング処理画像信号Ｒopn が原画像信号Ｒinに等しい画素については、クロージング処理画像信号Ｒcls を処理済画像信号Ｒout とし、オープニング処理画像信号Ｒopn またはクロージング処理画像信号Ｒcls のいずれも原画像信号Ｒinに等しくない画素については、オープニング処理画像信号Ｒopn とクロージング処理画像信号Ｒcls とを、重み付け係数w1，w2を用いた重み付け加算処理を適用する演算処理である。他の色Ｇ，Ｂの画像信号についても同様である。
【００４７】
ここで重み付け係数w1，w2は、いずれも 0.5を適用するものとする。
【００４８】
モルフォロジー演算処理手段10による上記オープニング処理およびクロージング処理に適用される構造要素（マスク）は、図５に示す４つの直線構造要素からなる多重構造要素（図示のものは５画素×５画素のマスクサイズとなる）である。多重構造要素を用いたオープニング処理およびクロージング処理は、多重構造要素を構成する各構造要素を適用したオープニング処理およびクロージング処理をそれぞれ各別に施し、オープニング処理についてはそれらの各別に行ったオープニング処理のうち最大値を、クロージング処理についてはそれらの各別に行ったクロージング処理のうち最小値を、それぞれ採る処理であり、画像のように２次元状に分布する信号値を処理する場合に有用である。なお、マスクサイズ（構造要素の大きさ）は、削除しようとする粒状の大きさに応じて変更することができる。一般的な画像においては、３画素×３画素から13画素×13画素程度の大きさとなる多重構造要素を適用するのが好適であり、３画素×３画素から７画素×７画素程度の大きさとなる多重構造要素を適用するのが最適である。
【００４９】
ただし、以下においては、説明の簡単化のため、２次元画像の一断面である１次元の信号分布を例にして説明するため、多重構造要素のうち１つの構造要素を用いた説明としている。
【００５０】
次に、本実施形態の画像処理装置の作用について、原画像信号の１つである図３（１）に示すＲ画像信号Ｒinを例に説明する。
【００５１】
まず、デジタルスチルカメラ等から、カラー画像を表すＲＧＢの画像信号Ｒin，Ｇin，Ｂinがモルフォロジー演算手段10に入力される。画像信号Ｒinは上述したように図３（１）に示すような分布を有する信号であり、他の画像信号Ｇin，Ｂinも所定の分布を有している（図示せず）。
【００５２】
原画像信号Ｒinには、比較的大きな周期の濃度分布に、この濃度分布より高周波のノイズが重畳している。この高周波ノイズ（同図（１）のＡ部）を拡大したものを同図（２）に示す。
【００５３】
モルフォロジー演算手段10は、入力された原画像信号Ｒin，Ｇin，Ｂinに対して各別に、予め設定された構造要素を用いたモルフォロジー演算によるオープニング処理およびクロージング処理をそれぞれ施し、オープニング処理画像信号Ｒopn ，Ｇopn ，Ｂopn およびクロージング処理画像信号Ｒcls ，Ｇcls ，Ｂcls を算出する。図３（３）は、原画像信号Ｒinに対するオープニング処理によるオープニング処理画像信号Ｒopn の、図３（２）に相当する部分（同図（１）のＡ部）を示す図であり、同図（４）は、原画像信号Ｒinに対するクロージング処理によるクロージング処理画像信号Ｒcls の、図３（２）に相当する部分（同図（１）のＡ部）を示す図である。
【００５４】
各図から解されるように、オープニング処理画像信号Ｒopn は、原画像信号Ｒin（図３（３）中において破線で示す）のうち、凸状の濃度変動部分（正のノイズ）が抑制された信号、クロージング処理画像信号Ｒcls は、原画像信号Ｒin（図３（４）中において破線で示す）のうち、凹状の濃度変動部分（負のノイズ）が抑制された信号とされる。
【００５５】
他のオープニング処理画像信号Ｇopn ，Ｂopn もオープニング処理画像信号Ｒopn と同様に、正のノイズが抑制された信号とされ、他のクロージング処理画像信号Ｇcls ，Ｂcls もクロージング処理画像信号Ｒcls と同様に、負のノイズが抑制された信号とされる。
【００５６】
モルフォロジー演算手段10により算出された各オープニング処理画像信号Ｒopn ，Ｇopn ，Ｂopn および各クロージング処理画像信号Ｒcls ，Ｇcls ，Ｂcls は、信号処理手段20に入力される。
【００５７】
信号処理手段20は、入力された各画像信号について、各色ＲＧＢごとに、画素を対応させて、オープニング処理画像信号とクロージング処理画像信号とを上記式（１）に従った演算処理を行い、各色ごとに処理済画像信号を算出する。
【００５８】
図３（５）は、同図（２）に示したオープニング処理画像信号Ｒopn と、同図（４）に示したクロージング処理画像信号Ｒcls とを式（１）に従って演算処理して得られた処理済画像信号Ｒout を示す。
【００５９】
図から解されるように、処理済画像信号Ｒout は、原画像信号Ｒin（図３（５）中において破線で示す）のうち、正および負のいずれのノイズも精度よく抑制された信号とされる。なお、重みづけ係数w1，w2の値はいずれも0.5 としたが、この値に限るものではなく、適当な範囲で変動させることもできる。
【００６０】
他の処理済画像信号Ｇout ，Ｂout も処理済画像信号Ｒout と同様に、正および負のいずれのノイズも精度よく抑制された信号とされる。
【００６１】
このように本実施形態の画像処理装置によれば、入力された画像信号に対して、モルフォロジー演算によるオープニング処理とクロージング処理とを各別に施して、オープニング処理により正のノイズが抑制された画像信号とクロージング処理により負のノイズが抑制された画像信号とをそれぞれ得、これらの得られた両画像信号を、画素を対応させて演算処理して１つの処理済画像信号を得るため、原画像に存在する正および負のノイズを精度よく抑制することができる。
【００６２】
なお、上記画像処理装置は、信号処理手段20による各演算処理ｇとして式（１）に示す演算処理を適用した実施形態であるが、本発明の画像処理方法および画像処理装置は、式（１）による演算処理のものに限るものではなく、種々の適切な演算処理を適用することができる。
【００６３】
例えば、下記式（２）に示す演算処理を適用することもできる。
【００６４】
【数２】
すなわち、オープニング処理画像信号Ｒopn が原画像信号Ｒinに等しい画素については、クロージング処理画像信号Ｒcls を処理済画像信号Ｒout とし、クロージング処理画像信号Ｒcls が原画像信号Ｒinに等しい画素については、オープニング処理画像信号Ｒopn を処理済画像信号Ｒout とし、オープニング処理画像信号Ｒopn またはクロージング処理画像信号Ｒcls のいずれも原画像信号Ｒinに等しくない画素については、オープニング処理画像信号Ｒopn とクロージング処理画像信号Ｒcls とのうち原画像信号Ｒinとの差が大きい方の処理画像信号（Ｒopn またはＲcls ）と、原画像信号Ｒinとの差（｜Ｒopn −Ｒin｜または｜Ｒcls −Ｒin｜）、による重みづけ加算した値を処理済画像信号Ｒout とする演算処理である。他の色Ｇ，Ｂの画像信号についても同様である。
【００６６】
このような演算処理を適用したものについても、式（１）に示した演算処理の実施形態と同様に、良好に正負のノイズを抑制することができる。
【００６７】
また本実施形態の画像処理装置においては、カラー画像を表すＲＧＢの画像信号を処理の対象としたが、本発明の画像処理装置はそのようなカラー画像信号を対象とするものに限るものではなく、例えば輝度や色差等の各種信号を処理の対象とすることができるものである。
【図面の簡単な説明】
【図１】本発明の画像処理方法の一実施形態の処理フローを示すフローチャート
【図２】図１に示した画像処理方法を実施する画像処理装置を示す図
【図３】（１）原画像信号Ｒinを示す図、（２）（１）に示した原画像信号Ｒinのうち高周波ノイズ（Ａ部）を示す拡大図、（３）原画像信号Ｒinに基づくオープニング処理画像信号Ｒopn のうちＡ部に相当する部分を示す図、（４）原画像信号Ｒinに基づくクロージング処理画像信号Ｒcls のうちＡ部に相当する部分を示す図、（５）処理済画像信号Ｒout のうちＡ部に相当する部分を示す図
【図４】モルフォロジー演算の基本的な作用を説明する図
【図５】多重構造要素を示す図
【符号の説明】
10 モルフォロジー演算処理手段
20 信号処理手段
Claims (6)
- 画像を構成する画素の値を規定する画像信号に対して、モルフォロジー演算によるオープニング処理およびクロージング処理を各別に施し、
前記オープニング処理により取得されたオープニング処理画像信号および前記クロージング処理により取得されたクロージング処理画像信号に基づいて、該両処理画像信号の画素を対応させた下記式（１）で定義される演算処理により処理済画像信号を取得することを特徴とする画像処理方法。
- 前記式（１）におけるw1およびw2の値がいずれも 0.5であることを特徴とする請求項１記載の画像処理方法。
- 画像を構成する画素の値を規定する画像信号に対して、モルフォロジー演算によるオープニング処理およびクロージング処理を各別に施し、
前記オープニング処理により取得されたオープニング処理画像信号および前記クロージング処理により取得されたクロージング処理画像信号に基づいて、該両処理画像信号の画素を対応させた下記式（２）で定義される演算処理により処理済画像信号を取得することを特徴とする画像処理方法。
- 画像を構成する画素の値を規定する画像信号に対して、モルフォロジー演算によるオープニング処理およびクロージング処理を各別に施すモルフォロジー演算処理手段と、
前記モルフォロジー演算処理手段による前記オープニング処理により取得されたオープニング処理画像信号および前記クロージング処理により取得されたクロージング処理画像信号に基づいて、該両処理画像信号の画素を対応させた下記式（１）で定義される演算処理により処理済画像信号を取得する信号処理手段とを備えたことを特徴とする画像処理装置。
- 前記信号処理手段による前記式（１）におけるw1およびw2の値がいずれも 0.5であることを特徴とする請求項４記載の画像処理装置。
- 画像を構成する画素の値を規定する画像信号に対して、モルフォロジー演算によるオープニング処理およびクロージング処理を各別に施すモルフォロジー演算処理手段と、
前記モルフォロジー演算処理手段による前記オープニング処理により取得されたオープニング処理画像信号および前記クロージング処理により取得されたクロージング処理画像信号に基づいて、該両処理画像信号の画素を対応させた下記式（２）で定義される演算処理により処理済画像信号を取得する信号処理手段とを備えたことを特徴とする画像処理装置。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP05481998A JP3932470B2 (ja) | 1998-03-06 | 1998-03-06 | 画像処理方法および装置 |
| US09/244,950 US6445831B1 (en) | 1998-02-10 | 1999-02-10 | Image processing method and apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP05481998A JP3932470B2 (ja) | 1998-03-06 | 1998-03-06 | 画像処理方法および装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH11259651A JPH11259651A (ja) | 1999-09-24 |
| JP3932470B2 true JP3932470B2 (ja) | 2007-06-20 |
Family
ID=12981312
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP05481998A Expired - Lifetime JP3932470B2 (ja) | 1998-02-10 | 1998-03-06 | 画像処理方法および装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3932470B2 (ja) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3730872B2 (ja) | 2000-03-31 | 2006-01-05 | 富士通株式会社 | 画像処理装置及び画像処理プログラム |
| JP3958036B2 (ja) * | 2001-12-10 | 2007-08-15 | 三菱電機株式会社 | 乱気流検出装置および乱気流検出方法 |
| JP4530834B2 (ja) * | 2003-12-26 | 2010-08-25 | 富士フイルム株式会社 | 超音波画像処理方法及び超音波画像処理装置、並びに、超音波画像処理プログラム |
| JP4652717B2 (ja) | 2004-04-26 | 2011-03-16 | 株式会社ミツトヨ | 画像処理装置及び方法並びにプログラム |
| FR2974220B1 (fr) * | 2011-04-18 | 2013-04-12 | Michelin Soc Tech | Analyse de l'image numerique de la surface interne d'un pneumatique - traitement des points de fausse mesure |
| JP7152958B2 (ja) * | 2019-01-10 | 2022-10-13 | 富士フイルムヘルスケア株式会社 | 超音波撮像装置及び画像処理方法 |
-
1998
- 1998-03-06 JP JP05481998A patent/JP3932470B2/ja not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JPH11259651A (ja) | 1999-09-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6373992B1 (en) | Method and apparatus for image processing | |
| JP4053185B2 (ja) | 画像処理方法および装置 | |
| US6667815B1 (en) | Method and apparatus for processing images | |
| US8675104B2 (en) | Noise filter for bayer pattern image data | |
| US7308152B2 (en) | Image processing method, apparatus, and program | |
| US6724942B1 (en) | Image processing method and system | |
| US20030228064A1 (en) | Multiresolution method of spatially filtering a digital image | |
| JP2004517384A (ja) | 動的な画像訂正および画像システム | |
| JP2004147333A (ja) | 選択的空間フィルタを使用するディジタル画像の階調及び空間特性の向上方法 | |
| JPH0622135A (ja) | 画像処理方法 | |
| JP2001118062A (ja) | ダイナミックレンジ自動圧縮方法 | |
| JPH05225332A (ja) | 空間的可変濾波の方法及び装置 | |
| JP3932470B2 (ja) | 画像処理方法および装置 | |
| JP2001285641A (ja) | 画像処理方法、画像処理装置および記録媒体 | |
| JP2000307868A (ja) | 画像処理方法および装置 | |
| JP3966987B2 (ja) | 画像処理方法および装置 | |
| JP2003132352A (ja) | 画像処理方法、画像処理装置、画像処理プログラム、これを記録したコンピュータ読み取り可能な記録媒体 | |
| JP4008087B2 (ja) | 画像処理方法および装置 | |
| JP4019239B2 (ja) | 画像鮮鋭化方法および画像鮮鋭化装置 | |
| JP2001008038A (ja) | 画像処理方法および画像処理装置およびマルチバンドカラー写真システム | |
| JP4053167B2 (ja) | 画像処理方法および装置 | |
| JP3804880B2 (ja) | ディジタル画像のノイズ抑制およびシャープネス強調のための画像処理方法 | |
| JP4445026B2 (ja) | 画像処理方法および装置並びにプログラム | |
| JPH11306344A (ja) | 画像処理方法および装置 | |
| JP3901837B2 (ja) | 画像処理方法および装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20040305 |
|
| A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20061128 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20070215 |
|
| 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: 20070306 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070307 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100330 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110330 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110330 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120330 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120330 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130330 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130330 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140330 Year of fee payment: 7 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| EXPY | Cancellation because of completion of term |