JP2942460B2 - 画像処理装置 - Google Patents
画像処理装置Info
- Publication number
- JP2942460B2 JP2942460B2 JP6148816A JP14881694A JP2942460B2 JP 2942460 B2 JP2942460 B2 JP 2942460B2 JP 6148816 A JP6148816 A JP 6148816A JP 14881694 A JP14881694 A JP 14881694A JP 2942460 B2 JP2942460 B2 JP 2942460B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- luminance
- pixel
- average value
- image processing
- 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
Links
- 238000004364 calculation method Methods 0.000 claims description 22
- 238000000034 method Methods 0.000 claims description 10
- 235000019557 luminance Nutrition 0.000 description 36
- 238000009499 grossing Methods 0.000 description 19
- 238000010586 diagram Methods 0.000 description 7
- 238000003672 processing method Methods 0.000 description 6
- 238000007796 conventional method Methods 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
Landscapes
- Image Processing (AREA)
Description
【０００１】
【産業上の利用分野】本発明は、画像のエッジのボケを
抑えてノイズを除去する画像処理方法及び画像処理装置
に関するものである。
抑えてノイズを除去する画像処理方法及び画像処理装置
に関するものである。
【０００２】
【従来の技術】従来の画像のノイズを除去するのに使わ
れている画像処理方法について以下に説明する。 ［１］平滑化方法 原画像の着目画素の輝度を、それを中心とするウィンド
ウの中の画素の輝度の平均で置き換える。いま、原画像
をf(i,j)、ウィンドウを一辺（２ｎ＋１）の正方形と
し、出力画像をg(i,j)、平均のための重み係数をh(k,l)
とすれば、
れている画像処理方法について以下に説明する。 ［１］平滑化方法 原画像の着目画素の輝度を、それを中心とするウィンド
ウの中の画素の輝度の平均で置き換える。いま、原画像
をf(i,j)、ウィンドウを一辺（２ｎ＋１）の正方形と
し、出力画像をg(i,j)、平均のための重み係数をh(k,l)
とすれば、
【０００３】
【数１】
【０００４】となる。
【０００５】重み係数h(k,l)は、次式を満たす。
【０００６】
【数２】
【０００７】図５は一辺３（ｎ＝１）のときの例であ
る。○は画素を、その中の数字は輝度を表す。原画像に
おいてウィンドウ内の中心の画素が着目画素で、それに
対応する出力画素の輝度はウィンドウ内の原画素の輝度
に重み係数を掛けたものの和を四捨五入したものになっ
ている。
る。○は画素を、その中の数字は輝度を表す。原画像に
おいてウィンドウ内の中心の画素が着目画素で、それに
対応する出力画素の輝度はウィンドウ内の原画素の輝度
に重み係数を掛けたものの和を四捨五入したものになっ
ている。
【０００８】［２］エッジ保存平滑化方法 原画像f(i,j)において、図６のような９つのウィンドウ
内の画素の分散を求め、この分散が最も小さいウィンド
ウ内の原画素の平均値を出力画素g(i,j)とする。 この
フィルタのアルゴリズムは次のような考えに基づいてい
る。エッジがウィンドウ内に含まれているとすると、分
散は大きくなり、逆に最も分散の小さいウィンドウには
エッジが含まれていないと考えられる。そこで、最も分
散の小さいウィンドウ内の画素を平滑化することによ
り、エッジのぼけを防ぐ。
内の画素の分散を求め、この分散が最も小さいウィンド
ウ内の原画素の平均値を出力画素g(i,j)とする。 この
フィルタのアルゴリズムは次のような考えに基づいてい
る。エッジがウィンドウ内に含まれているとすると、分
散は大きくなり、逆に最も分散の小さいウィンドウには
エッジが含まれていないと考えられる。そこで、最も分
散の小さいウィンドウ内の画素を平滑化することによ
り、エッジのぼけを防ぐ。
【０００９】
【発明が解決しようとする課題】従来のノイズ除去のた
めの画像処理方法は、中間調画像の性質を生かしていな
いため、画像情報として重要なエッジの処理が十分では
なかった。このことを数式と実験例を用いて説明する。
中間調画像の性質を見るために、図７にＳＩＤＢＡ（標
準画像データベース、Standard Image Data-Base）のgi
rl（512×512ドット、一画素８ビット）の横軸をある直
線上の画素位置（上から351番目の画素列）、縦軸に輝
度をとったグラフを示す。図７を見ると判るように１画
素の単位で輝度が変化している。しかし、画像情報とし
て重要なのは数画素単位の大きな変化であり、細かい変
化は返って目障りなもので、ノイズとみなされる。そこ
で、画像の大きな変化を残し、細かい変化を取り除くこ
とを考える。
めの画像処理方法は、中間調画像の性質を生かしていな
いため、画像情報として重要なエッジの処理が十分では
なかった。このことを数式と実験例を用いて説明する。
中間調画像の性質を見るために、図７にＳＩＤＢＡ（標
準画像データベース、Standard Image Data-Base）のgi
rl（512×512ドット、一画素８ビット）の横軸をある直
線上の画素位置（上から351番目の画素列）、縦軸に輝
度をとったグラフを示す。図７を見ると判るように１画
素の単位で輝度が変化している。しかし、画像情報とし
て重要なのは数画素単位の大きな変化であり、細かい変
化は返って目障りなもので、ノイズとみなされる。そこ
で、画像の大きな変化を残し、細かい変化を取り除くこ
とを考える。
【００１０】着目画素を原点とした座標軸をrとし、画
素の位置をri、輝度をfiとすれば、次式が成り立つ。 fi = a+bri+c(ri) (3) 図８はrとfの関係のグラフである。画像の階調数をＲと
している。ここで、c(r)を次式が成り立つようにしても
一般性は失わない。 E［c(ri)］= 0 (4) ただし、E［x］はxの平均を表す。輝度変化を(3)式のよ
うに表すと、c(r)をノイズと見なせ、理想的にはノイズ
をとった輝度は次式のようになる。 fi = a+bri (5) 以上の議論に従えば、画像処理によりノイズを除去する
ことはa+briの推定値を求めることに帰着される。統計
学において、aとbの値の推定値α、βを回帰係数と言
い、 f=α+βr (6) をfのrへの回帰直線と言う。回帰係数の求めかたにはい
ろいろあるが、c(ri)が正規分布に従うとしたとき、最
も信頼性のある求めかたは最小２乗法と言われる次式を
最小にするα、βを求めることである。 Se = Σc(ri)2 = Σ[f-α-βri]2 (7) 従って、
素の位置をri、輝度をfiとすれば、次式が成り立つ。 fi = a+bri+c(ri) (3) 図８はrとfの関係のグラフである。画像の階調数をＲと
している。ここで、c(r)を次式が成り立つようにしても
一般性は失わない。 E［c(ri)］= 0 (4) ただし、E［x］はxの平均を表す。輝度変化を(3)式のよ
うに表すと、c(r)をノイズと見なせ、理想的にはノイズ
をとった輝度は次式のようになる。 fi = a+bri (5) 以上の議論に従えば、画像処理によりノイズを除去する
ことはa+briの推定値を求めることに帰着される。統計
学において、aとbの値の推定値α、βを回帰係数と言
い、 f=α+βr (6) をfのrへの回帰直線と言う。回帰係数の求めかたにはい
ろいろあるが、c(ri)が正規分布に従うとしたとき、最
も信頼性のある求めかたは最小２乗法と言われる次式を
最小にするα、βを求めることである。 Se = Σc(ri)2 = Σ[f-α-βri]2 (7) 従って、
【００１１】
【数８】
【００１２】
【数９】
【００１３】より、 β = Σ(fi-fc)・(ri-rc)/Σ(ri-rc)2 (10) α = fc-β・rc (11) ここで、 fc = E[fi] (12) rc = E[ri] (13) 以上より、着目画素r=0でのノイズを除去した輝度は f0 = E[fi] - E[ri]Σ(fi-fc)・(ri-rc)/Σ(ri-rc)2 (14) となる。
【００１４】ところが、従来の技術である平滑化、エッ
ジ保存平滑化方法においては、出力画素f0は原画像の平
均値であるので、 f0 = E[fi] (15) 上式と(11)、(12)式より、 α = fc (16) β = 0 (17) となる。つまり輝度の傾きを考慮にいれていない。その
ためにノイズ除去と同時に、エッジなどの輝度の大きな
変化を保存できない。このことを図９の一次元画像を使
って示す。図９において、（ａ）は回帰直線を用いた方
法、（ｂ）は平滑化、（ｃ）はエッジ保存平滑化を表
す。エッジは輝度の傾斜が変化するところであるので、
（ａ）の様にそこにちょうど着目画素がきても、２つの
回帰直線によって原画像が近似されるので、出力画素の
輝度を（α1＋α2）／２とすれば、エッジが保存でき
る。これに対して（ｂ）の平滑化では、着目画素の輝度
α3は原画像より小さくなっており、このことはエッジ
がぼけたことに対応する。（ｃ）のエッジ保存平滑化で
は、ウィンドウが複数あり場所によって使われるウィン
ドウが不連続に変わるので、出力画像が安定しない。
（ｃ）では分散の小さい方のα5が出力画素の輝度にな
る。では、実際の画像に従来の方法を適用したらどうな
るかをみてみる。
ジ保存平滑化方法においては、出力画素f0は原画像の平
均値であるので、 f0 = E[fi] (15) 上式と(11)、(12)式より、 α = fc (16) β = 0 (17) となる。つまり輝度の傾きを考慮にいれていない。その
ためにノイズ除去と同時に、エッジなどの輝度の大きな
変化を保存できない。このことを図９の一次元画像を使
って示す。図９において、（ａ）は回帰直線を用いた方
法、（ｂ）は平滑化、（ｃ）はエッジ保存平滑化を表
す。エッジは輝度の傾斜が変化するところであるので、
（ａ）の様にそこにちょうど着目画素がきても、２つの
回帰直線によって原画像が近似されるので、出力画素の
輝度を（α1＋α2）／２とすれば、エッジが保存でき
る。これに対して（ｂ）の平滑化では、着目画素の輝度
α3は原画像より小さくなっており、このことはエッジ
がぼけたことに対応する。（ｃ）のエッジ保存平滑化で
は、ウィンドウが複数あり場所によって使われるウィン
ドウが不連続に変わるので、出力画像が安定しない。
（ｃ）では分散の小さい方のα5が出力画素の輝度にな
る。では、実際の画像に従来の方法を適用したらどうな
るかをみてみる。
【００１５】ＳＩＤＢＡの標準画像girl（512×512ドッ
ト、一画素８ビット）を原画像として、平滑化（３×３
ウィンドウ、重み係数(a)）とエッジ保存平滑化を行っ
た。
ト、一画素８ビット）を原画像として、平滑化（３×３
ウィンドウ、重み係数(a)）とエッジ保存平滑化を行っ
た。
【００１６】画像の性質を示すために、図１０、図１１
にそれぞれ平滑化よる画像、エッジ保存平滑化による画
像の横軸をある直線上の画素位置（上から351番目の画
素列）、縦軸に輝度をとったものを示す。図を見てわか
るように、平滑化による画像は原画像に比べてノイズが
小さくなっているが、画素位置に対する輝度の変化がな
だらかになっており、エッジがぼけている。また、エッ
ジ保存平滑化による画像はノイズが小さくなっている
が、エッジが逆に強調され、画素位置に対する輝度変化
が階段状になっており、原画像と性質の異なる画像にな
ってしまっている。本発明は以上の点を鑑みなされたも
のであり、画像の輝度の傾きを考慮に入れ、回帰直線を
用いることにより、エッジをぼかさず、強調せず、原画
像の性質を保って、ノイズを除去することができる画像
処理方法および画像処理装置を提供することを目的とす
る。
にそれぞれ平滑化よる画像、エッジ保存平滑化による画
像の横軸をある直線上の画素位置（上から351番目の画
素列）、縦軸に輝度をとったものを示す。図を見てわか
るように、平滑化による画像は原画像に比べてノイズが
小さくなっているが、画素位置に対する輝度の変化がな
だらかになっており、エッジがぼけている。また、エッ
ジ保存平滑化による画像はノイズが小さくなっている
が、エッジが逆に強調され、画素位置に対する輝度変化
が階段状になっており、原画像と性質の異なる画像にな
ってしまっている。本発明は以上の点を鑑みなされたも
のであり、画像の輝度の傾きを考慮に入れ、回帰直線を
用いることにより、エッジをぼかさず、強調せず、原画
像の性質を保って、ノイズを除去することができる画像
処理方法および画像処理装置を提供することを目的とす
る。
【００１７】
【課題を解決するための手段】本発明の画像処理装置
は、画像のエッジをぼかさず、強調せずにノイズを除去
する装置であって、原画像の画像データを格納するため
の記憶手段と、該記憶手段から抽出した着目画素を原点
とした該着目画素を含む直線上の画素を用いて輝度の回
帰係数を求める回帰係数演算手段と、該回帰係数演算手
段から出力された複数の回帰係数の平均値を求める平均
値演算手段とを備え、該平均値演算手段において、前記
直線上の画素から求められた輝度の回帰係数の平均値に
基づいて画像のノイズ除去処理を行うものである。
は、画像のエッジをぼかさず、強調せずにノイズを除去
する装置であって、原画像の画像データを格納するため
の記憶手段と、該記憶手段から抽出した着目画素を原点
とした該着目画素を含む直線上の画素を用いて輝度の回
帰係数を求める回帰係数演算手段と、該回帰係数演算手
段から出力された複数の回帰係数の平均値を求める平均
値演算手段とを備え、該平均値演算手段において、前記
直線上の画素から求められた輝度の回帰係数の平均値に
基づいて画像のノイズ除去処理を行うものである。
【００１８】
【作用】上記画像処理方法では、画像の輝度の回帰分析
に基づいて輝度の回帰係数が求められ、この回帰係数に
基づいて画像のノイズが除去される。また、上記画像処
理装置では、回帰係数演算手段が原画像から回帰係数を
求め、これに基づいて画像処理が行われて画像のノイズ
が除去される。また、平均値演算手段は前記回帰係数演
算手段により求められた複数の回帰係数の平均値を求
め、これを出力し、これに基づいて画像のノイズが除去
される。上記のように、本発明の画像処理方法及び画像
処理装置では、エッジをぼかさず、強調せず、原画像の
性質を保って、ノイズを除去する。
に基づいて輝度の回帰係数が求められ、この回帰係数に
基づいて画像のノイズが除去される。また、上記画像処
理装置では、回帰係数演算手段が原画像から回帰係数を
求め、これに基づいて画像処理が行われて画像のノイズ
が除去される。また、平均値演算手段は前記回帰係数演
算手段により求められた複数の回帰係数の平均値を求
め、これを出力し、これに基づいて画像のノイズが除去
される。上記のように、本発明の画像処理方法及び画像
処理装置では、エッジをぼかさず、強調せず、原画像の
性質を保って、ノイズを除去する。
【００１９】
【実施例】以下、図面を参照して本発明の一実施例を詳
細に説明する。図１に本発明による画像処理装置の一実
施例を示す。図１に示す本発明の一実施例である画像処
理装置は、ライン単位で入力される画像データを格納す
るためのラインバッファ１、ラインバッファ１に接続さ
れており、ラインバッファ１より抽出した着目画素を原
点とした直線rk（k=1,2,‥8）上の標本画素rk=rk,i(i=
0,1,‥‥m-1、mは正の整数)から回帰係数を求める回帰
係数演算部２のk（k=1,2,‥8）と、回帰係数演算部２の
k（k=1,2,‥8）より出力された複数の回帰係数の平均値
を求める平均値演算部３、平均値演算部３に接続されて
おり、平均値演算部３から出力される画素を格納するラ
インバッファ４から構成される。
細に説明する。図１に本発明による画像処理装置の一実
施例を示す。図１に示す本発明の一実施例である画像処
理装置は、ライン単位で入力される画像データを格納す
るためのラインバッファ１、ラインバッファ１に接続さ
れており、ラインバッファ１より抽出した着目画素を原
点とした直線rk（k=1,2,‥8）上の標本画素rk=rk,i(i=
0,1,‥‥m-1、mは正の整数)から回帰係数を求める回帰
係数演算部２のk（k=1,2,‥8）と、回帰係数演算部２の
k（k=1,2,‥8）より出力された複数の回帰係数の平均値
を求める平均値演算部３、平均値演算部３に接続されて
おり、平均値演算部３から出力される画素を格納するラ
インバッファ４から構成される。
【００２０】次に、この画像処理装置の各部の動作を説
明する。先ず、回帰係数演算部２のk（k=1,2,‥8）につ
いて説明する。 ＜回帰係数演算部２のk（k=1,2,‥8）＞図２の様に画像
のｘｙ平面上の点(x,y)を原点とした８方向の座標軸rk
（k=1,2,‥8）をとる。rk軸とx軸とのなす角度θkは θk = kπ／４ (18) とする。回帰係数演算部はラインバッファ１に格納され
ている画像データから、８方向ともそれぞれｍ個の標本
画素rk=rk,i(i=0,1,‥‥m-1、mは正の整数)の輝度fk,i
から回帰直線を求め、８つの回帰係数αk（k=1,2,‥8）
を出力する。ここで、回帰直線fk=αk+βkrkの傾きβk
と切片αkは次の式で与えられる。
明する。先ず、回帰係数演算部２のk（k=1,2,‥8）につ
いて説明する。 ＜回帰係数演算部２のk（k=1,2,‥8）＞図２の様に画像
のｘｙ平面上の点(x,y)を原点とした８方向の座標軸rk
（k=1,2,‥8）をとる。rk軸とx軸とのなす角度θkは θk = kπ／４ (18) とする。回帰係数演算部はラインバッファ１に格納され
ている画像データから、８方向ともそれぞれｍ個の標本
画素rk=rk,i(i=0,1,‥‥m-1、mは正の整数)の輝度fk,i
から回帰直線を求め、８つの回帰係数αk（k=1,2,‥8）
を出力する。ここで、回帰直線fk=αk+βkrkの傾きβk
と切片αkは次の式で与えられる。
【００２１】
【数１９】
【００２２】 αk=fk,c-βk・rk,c (20) ここで、
【００２３】
【数２１】
【００２４】
【数２２】
【００２５】ただし、位置の単位は勝手に決められるの
で、標本画素の位置rk=rk,iを rk,i=i (i=0,1,‥‥m-1) (23) とすれば、標本画素数ｍが決まれば画素位置rk,iに関す
る演算は一意的に決まる。そこで、(24)式の様にCkを定
義する。
で、標本画素の位置rk=rk,iを rk,i=i (i=0,1,‥‥m-1) (23) とすれば、標本画素数ｍが決まれば画素位置rk,iに関す
る演算は一意的に決まる。そこで、(24)式の様にCkを定
義する。
【００２６】
【数２４】
【００２７】すると(20)式は次の様に変形できる。
【００２８】
【数２５】
【００２９】従って、αkは輝度fk,i(i=0,1,‥,m-1)の
みの関数となるので、回帰係数演算部は入力がfk,i(i=
0,1,‥,m-1)のみで構成できる。(25)式をもとに構成し
た回帰係数演算部の一例を図３に示す。図３の回帰係数
演算部２のkは、ラインバッファ１に接続されており、
ラインバッファ１より抽出した着目画素を原点とした直
線rk上の標本画素rk=rk,i(i=0,1,‥m-1)の輝度fk,i(i=
0,1,‥m-1)の平均を求める平均値演算部５、ラインバッ
ファ１と平均値演算部５に接続されており輝度fk,i(i=
0,1,‥m-1)と平均値演算部５から出力された平均輝度
fk,cの差を計算し出力する減算部６、減算部６とメモリ
８に接続されておりメモリに格納されている重み係数Ci
と輝度と平均輝度との差fk,i-fk,cの内積を計算し出力
する内積演算部７、内積演算部７に接続され重み係数Ci
を格納しており内積演算部に供給するメモリ８、平均値
演算部５と内積演算部７に接続されており輝度から内積
を引いた値を回帰係数αiとして出力する減算部９から
構成される。
みの関数となるので、回帰係数演算部は入力がfk,i(i=
0,1,‥,m-1)のみで構成できる。(25)式をもとに構成し
た回帰係数演算部の一例を図３に示す。図３の回帰係数
演算部２のkは、ラインバッファ１に接続されており、
ラインバッファ１より抽出した着目画素を原点とした直
線rk上の標本画素rk=rk,i(i=0,1,‥m-1)の輝度fk,i(i=
0,1,‥m-1)の平均を求める平均値演算部５、ラインバッ
ファ１と平均値演算部５に接続されており輝度fk,i(i=
0,1,‥m-1)と平均値演算部５から出力された平均輝度
fk,cの差を計算し出力する減算部６、減算部６とメモリ
８に接続されておりメモリに格納されている重み係数Ci
と輝度と平均輝度との差fk,i-fk,cの内積を計算し出力
する内積演算部７、内積演算部７に接続され重み係数Ci
を格納しており内積演算部に供給するメモリ８、平均値
演算部５と内積演算部７に接続されており輝度から内積
を引いた値を回帰係数αiとして出力する減算部９から
構成される。
【００３０】次に、平均値演算部３について説明する。 ＜平均値演算部３＞平均値演算部３は回帰係数演算部２
のk（k=1,2,‥8）から出力されたαk（k=1,2,‥8）を使
って出力画素g(x,y)を計算し、出力してラインバッファ
４に蓄える。式で表すと、
のk（k=1,2,‥8）から出力されたαk（k=1,2,‥8）を使
って出力画素g(x,y)を計算し、出力してラインバッファ
４に蓄える。式で表すと、
【００３１】
【数２６】
【００３２】となる。本実施例の画像処理装置による処
理画像（標本画素数m=4のとき）の輝度のグラフを図４
に示す。平滑化やエッジ保存平滑化による画像にくらべ
て、エッジなど、原画像の性質を保存し、かつノイズが
除去されていることが判る。
理画像（標本画素数m=4のとき）の輝度のグラフを図４
に示す。平滑化やエッジ保存平滑化による画像にくらべ
て、エッジなど、原画像の性質を保存し、かつノイズが
除去されていることが判る。
【００３３】
【発明の効果】本発明の画像処理装置は、原画像の画像
データを格納するための記憶手段と、該記憶手段から抽
出した着目画素を原点とした該着目画素を含む直線上の
画素を用いて輝度の回帰係数を求める回帰係数演算手段
と、該回帰係数演算手段から出力された複数の回帰係数
の平均値を求める平均値演算手段とを備え、該平均値演
算手段によって求められた輝度の回帰係数の平均値に基
づいて画像のノイズを除去するものであるため、画像処
理の際、エッジがぼかされず、しかもエッジが強調され
ることなく、また、高い精度のノイズ除去が行われる。
データを格納するための記憶手段と、該記憶手段から抽
出した着目画素を原点とした該着目画素を含む直線上の
画素を用いて輝度の回帰係数を求める回帰係数演算手段
と、該回帰係数演算手段から出力された複数の回帰係数
の平均値を求める平均値演算手段とを備え、該平均値演
算手段によって求められた輝度の回帰係数の平均値に基
づいて画像のノイズを除去するものであるため、画像処
理の際、エッジがぼかされず、しかもエッジが強調され
ることなく、また、高い精度のノイズ除去が行われる。
【図１】本発明の一実施例である画像処理装置のブロッ
ク図である。
ク図である。
【図２】図１の画像処理装置の動作を説明するための座
標軸rk上の着目画素と標本画素の関係を示した説明図で
ある。
標軸rk上の着目画素と標本画素の関係を示した説明図で
ある。
【図３】図１に示した画像処理装置における回帰係数演
算部のブロック図である。
算部のブロック図である。
【図４】本発明による画像処理画像の輝度と画素位置の
関係を示したグラフである。
関係を示したグラフである。
【図５】従来の平滑化方法を説明するための説明図であ
る。
る。
【図６】従来のエッジ保存平滑化方法を説明するための
ウィンドウの説明図である。
ウィンドウの説明図である。
【図７】原画像の画素位置と輝度の関係を示したグラフ
である。
である。
【図８】画像を直線で近似したときの数式とグラフの対
応を関係を示した図である。
応を関係を示した図である。
【図９】本発明と従来の方法の原理の数式とグラフの対
応を示した図である。
応を示した図である。
【図１０】従来の平滑化方法による画像の輝度と画素位
置の関係を示したグラフである。
置の関係を示したグラフである。
【図１１】従来のエッジ保存平滑化方法による画像の輝
度と画素の位置の関係を示したグラフである。
度と画素の位置の関係を示したグラフである。
１、４ ラインバッファ ２の１〜２の８ 回帰係数演算部 ３、５ 平均値演算部 ６、９ 減算部 ７ 内積演算部 ８ メモリ
Claims (1)
- 【請求項１】 画像のノイズを除去する画像処理装置で
あって、原画像の画像データを格納するための記憶手段
と、該記憶手段から抽出した着目画素を原点とした該着
目画素を含む直線上の画素を用いて輝度の回帰係数を求
める回帰係数演算手段と、該回帰係数演算手段から出力
された複数の回帰係数の平均値を求める平均値演算手段
とを備え、 該平均値演算手段において、前記直線上の画素から求め
られた輝度の回帰係数の平均値に基づいて画像のノイズ
除去処理を行うことを特徴とした画像処理装置 。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP6148816A JP2942460B2 (ja) | 1994-06-30 | 1994-06-30 | 画像処理装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP6148816A JP2942460B2 (ja) | 1994-06-30 | 1994-06-30 | 画像処理装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0816770A JPH0816770A (ja) | 1996-01-19 |
| JP2942460B2 true JP2942460B2 (ja) | 1999-08-30 |
Family
ID=15461364
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP6148816A Expired - Fee Related JP2942460B2 (ja) | 1994-06-30 | 1994-06-30 | 画像処理装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2942460B2 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4311794B2 (ja) * | 1999-01-29 | 2009-08-12 | オリンパス株式会社 | 画像処理装置及び画像処理プログラムを格納した記録媒体 |
| JP5245715B2 (ja) * | 2008-09-30 | 2013-07-24 | Ｎｋワークス株式会社 | 写真画像処理方法、写真画像処理プログラム、及び写真画像処理装置 |
-
1994
- 1994-06-30 JP JP6148816A patent/JP2942460B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0816770A (ja) | 1996-01-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7257273B2 (en) | Hierarchical scheme for blur detection in digital image using wavelet transform | |
| US6707952B1 (en) | Method for removing ringing artifacts from locations near dominant edges of an image reconstructed after compression | |
| WO2016206087A1 (zh) | 一种低照度图像处理方法和装置 | |
| JP4727720B2 (ja) | 画像処理方法および画像処理装置 | |
| US6801672B1 (en) | Removing noise from a color image using wavelets | |
| US20080253678A1 (en) | Denoise method on image pyramid | |
| EP1909227B1 (en) | Method of and apparatus for minimizing ringing artifacts in an input image | |
| US8965141B2 (en) | Image filtering based on structural information | |
| CN109377450B (zh) | 一种边缘保护的去噪方法 | |
| US6803954B1 (en) | Filtering control method for improving image quality of bi-linear interpolated image | |
| EP2164040A1 (en) | System and method for high quality image and video upscaling | |
| Scollar et al. | Image enhancement using the median and the interquartile distance | |
| WO2016189901A1 (ja) | 画像処理装置、画像処理方法、プログラム、これを記録した記録媒体、映像撮影装置、及び映像記録再生装置 | |
| US6574374B1 (en) | Enhancing image compression performance by morphological processing | |
| US6289133B1 (en) | Image processing method and apparatus | |
| JP2001078025A (ja) | 画像の輝度補正方法 | |
| JP2005025696A (ja) | 画像処理装置、画像処理方法及びプログラム及び記録媒体 | |
| CN111383190B (zh) | 图像处理设备和方法 | |
| CN113344801A (zh) | 一种应用于燃气计量设施环境下的图像增强方法、系统、终端及存储介质 | |
| Patil et al. | Bilateral filter for image denoising | |
| JP2942460B2 (ja) | 画像処理装置 | |
| US6728415B2 (en) | Method and apparatus for image processing using adaptive convolution filters | |
| US20040184669A1 (en) | Method for removing ringing artifacts from locations near dominant edges of an image reconstructed after compression | |
| Tigga et al. | Image deblurring with impulse noise using alternating direction method of multipliers and Lucy-Richardson method | |
| Immerkær | Use of blur-space for deblurring and edge-preserving noise smoothing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090618 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100618 Year of fee payment: 11 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100618 Year of fee payment: 11 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110618 Year of fee payment: 12 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120618 Year of fee payment: 13 |
|
| LAPS | Cancellation because of no payment of annual fees |