JP2002141778A - ディジタルフィルタ - Google Patents
ディジタルフィルタInfo
- Publication number
- JP2002141778A JP2002141778A JP2000334603A JP2000334603A JP2002141778A JP 2002141778 A JP2002141778 A JP 2002141778A JP 2000334603 A JP2000334603 A JP 2000334603A JP 2000334603 A JP2000334603 A JP 2000334603A JP 2002141778 A JP2002141778 A JP 2002141778A
- Authority
- JP
- Japan
- Prior art keywords
- filter
- characteristic
- frequency
- digital filter
- image
- 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.)
- Pending
Links
Landscapes
- Image Processing (AREA)
Abstract
(57)【要約】
【課題】 少ないタップ数で良好な結果を与える画像処
理用ディジタルフィルタを提供する。 【解決手段】 画像信号の自己回帰モデルのパワースペ
クトルを重み関数として設計することにより、周波数軸
上の極値点における目標特性との差が周波数0の点およ
び折り返し点に近づくにつれ徐々に小さくなる特性のフ
ィルタとする。
理用ディジタルフィルタを提供する。 【解決手段】 画像信号の自己回帰モデルのパワースペ
クトルを重み関数として設計することにより、周波数軸
上の極値点における目標特性との差が周波数0の点およ
び折り返し点に近づくにつれ徐々に小さくなる特性のフ
ィルタとする。
Description
【0001】
【発明の属する技術分野】本発明は、ディジタルフィル
タ、特に画像の拡大および縮少の際に用いられる画像処
理用ディジタルフィルタに関する。
タ、特に画像の拡大および縮少の際に用いられる画像処
理用ディジタルフィルタに関する。
【0002】
【従来の技術】図1(a)は画像の一方の座標軸方向の
画素値の分布の一例を示している。また、図2(a)は
それを周波数領域で表わしたものである。図1(b)に
示すように、画素と画素の値に0を挿入することにより
画素数を2倍にすると、周波数領域では、図2(b)に
示すように、処理周波数範囲の中央、すなわち、処理周
波数範囲を0〜2πに規格化すると、ω=πの点を中心
にエイリアス成分を生じる(図では0〜1に規格化され
ているので0.5の点に生じる)。これをフィルタで除
去することにより、図2(c)に示すような原画像の周
波数スペクトル(図2(a))と類似した周波数スペク
トルが得られる。
画素値の分布の一例を示している。また、図2(a)は
それを周波数領域で表わしたものである。図1(b)に
示すように、画素と画素の値に0を挿入することにより
画素数を2倍にすると、周波数領域では、図2(b)に
示すように、処理周波数範囲の中央、すなわち、処理周
波数範囲を0〜2πに規格化すると、ω=πの点を中心
にエイリアス成分を生じる(図では0〜1に規格化され
ているので0.5の点に生じる)。これをフィルタで除
去することにより、図2(c)に示すような原画像の周
波数スペクトル(図2(a))と類似した周波数スペク
トルが得られる。
【0003】この操作をx軸方向のすべてのyの値につ
いて行ない、さらに、y軸方向について同様に行なうこ
とにより、原画像の周波数スペクトルと類似した周波数
スペクトルを有する拡大画像が得られる。また、画像縮
少の際には、画素の間引きによる折り返しノイズを除去
するため、画素の間引きの前にローパスフィルタにより
高域成分が除去される。
いて行ない、さらに、y軸方向について同様に行なうこ
とにより、原画像の周波数スペクトルと類似した周波数
スペクトルを有する拡大画像が得られる。また、画像縮
少の際には、画素の間引きによる折り返しノイズを除去
するため、画素の間引きの前にローパスフィルタにより
高域成分が除去される。
【0004】画素拡大の際のエイリアス成分を完全に除
去する理想フィルタの特性は、
去する理想フィルタの特性は、
【0005】
【数2】
【0006】である。この、エイリアス成分の除去は時
間(空間)領域で見ると、ディジタルフィルタによって
0値が値を持った信号に変えられたように見える(図1
(c))。画像の拡大という意味では画素と画素の間が
補間された結果となるので、この際に用いるディジタル
フィルタは補間フィルタと呼ばれる。
間(空間)領域で見ると、ディジタルフィルタによって
0値が値を持った信号に変えられたように見える(図1
(c))。画像の拡大という意味では画素と画素の間が
補間された結果となるので、この際に用いるディジタル
フィルタは補間フィルタと呼ばれる。
【0007】最も簡単で広く使用されている補間フィル
タは単純な線形補間演算であるが、そのインパルス応答
を図3に示す。周波数特性をみると、図4のようにな
り、このようななだらかな周波数特性のフィルタは優れ
たフィルタとは言えず、画質劣化が著しい。なお、図4
に示した特性はエイリアス成分除去のためのバンドリジ
ェクト特性ではなくローパス特性となっている。しか
し、処理周波数範囲を2倍に拡大すれば右側に折り返し
が出現するので、バンドリジェクト特性となる。本明細
書では、ローパス特性のフィルタとバンドリジェクト特
性のフィルタのいずれも、図4のような形式で表わすも
のとする。
タは単純な線形補間演算であるが、そのインパルス応答
を図3に示す。周波数特性をみると、図4のようにな
り、このようななだらかな周波数特性のフィルタは優れ
たフィルタとは言えず、画質劣化が著しい。なお、図4
に示した特性はエイリアス成分除去のためのバンドリジ
ェクト特性ではなくローパス特性となっている。しか
し、処理周波数範囲を2倍に拡大すれば右側に折り返し
が出現するので、バンドリジェクト特性となる。本明細
書では、ローパス特性のフィルタとバンドリジェクト特
性のフィルタのいずれも、図4のような形式で表わすも
のとする。
【0008】画像処理用フィルタとしては、等リプル特
性のフィルタや最大平坦特性のフィルタも使用される。
図5および図6に等リプル特性および最大平坦特性のフ
ィルタの周波数特性を示す。等リプル特性のフィルタと
は、リプルの振幅が通過域および遮断域でそれぞれ一定
値となるように設計されたフィルタである。最大平坦特
性のフィルタとは、ω=0および折り返し周波数点ω=
πでの平坦さを最大にするため微分係数が可能な限り0
になるように設計されたフィルタである。
性のフィルタや最大平坦特性のフィルタも使用される。
図5および図6に等リプル特性および最大平坦特性のフ
ィルタの周波数特性を示す。等リプル特性のフィルタと
は、リプルの振幅が通過域および遮断域でそれぞれ一定
値となるように設計されたフィルタである。最大平坦特
性のフィルタとは、ω=0および折り返し周波数点ω=
πでの平坦さを最大にするため微分係数が可能な限り0
になるように設計されたフィルタである。
【0009】
【発明が解決しようとする課題】等リプル特性のフィル
タでは低い次数で急峻な遮断特性を得られるが、「リン
ギング」(画像のエッジ付近で縞模様が生じる)が大き
いという欠点がある。最大平坦特性のフィルタでは緩や
かな遮断特性から「リンギング」が小さいと言われてい
るが、画像の「ぼけ」を小さくするためには次数を上げ
ねばならず、結局「リンギング」が生じてしまう。
タでは低い次数で急峻な遮断特性を得られるが、「リン
ギング」(画像のエッジ付近で縞模様が生じる)が大き
いという欠点がある。最大平坦特性のフィルタでは緩や
かな遮断特性から「リンギング」が小さいと言われてい
るが、画像の「ぼけ」を小さくするためには次数を上げ
ねばならず、結局「リンギング」が生じてしまう。
【0010】したがって本発明の目的は、画像の「ぼ
け」「ぎざぎざ(jaggy)」「リンギング」のどの
面でも優れたフィルタ特性を見い出し、その設計法も確
立すること。それによりディジタルカメラにおけるディ
ジタルズームやその他の画像処理において、画質を確保
した上で拡大・縮小を可能にすることにある。
け」「ぎざぎざ(jaggy)」「リンギング」のどの
面でも優れたフィルタ特性を見い出し、その設計法も確
立すること。それによりディジタルカメラにおけるディ
ジタルズームやその他の画像処理において、画質を確保
した上で拡大・縮小を可能にすることにある。
【0011】
【課題を解決するための手段】本発明によれば、周波数
軸上の極値点における目標特性との差が周波数0に近づ
くにつれ徐々に小さくなる特性を有するディジタルフィ
ルタが提供される。前記目標特性は、所定の通過域およ
び阻止域を有する特性である。また、前記差は、画像信
号の自己回帰モデルのパワースペクトルの逆数に応じて
変化することが望ましい。
軸上の極値点における目標特性との差が周波数0に近づ
くにつれ徐々に小さくなる特性を有するディジタルフィ
ルタが提供される。前記目標特性は、所定の通過域およ
び阻止域を有する特性である。また、前記差は、画像信
号の自己回帰モデルのパワースペクトルの逆数に応じて
変化することが望ましい。
【0012】人間の目は画像の直流成分(平坦な画像)
に対する感度が高く、画像信号の自己回帰モデルで重み
付けすることは、人間の視覚特性と合っており、好都合
である。目標特性に対して、直流などの低い周波数成分
の精度を上げ、重要性の低い帯域で誤差を許すことで、
低い次数で良好な画質を得る、総合的に優れたフィルタ
を設計することができる。
に対する感度が高く、画像信号の自己回帰モデルで重み
付けすることは、人間の視覚特性と合っており、好都合
である。目標特性に対して、直流などの低い周波数成分
の精度を上げ、重要性の低い帯域で誤差を許すことで、
低い次数で良好な画質を得る、総合的に優れたフィルタ
を設計することができる。
【0013】
【発明の実施の形態】本発明では、画像信号は一般的に
隣り合う成分の相関が強く、一次の自己回帰モデルで良
く近似できるという性質を利用して、一次の自己回帰モ
デルのパワースペクトルをフィルタ設計の際の重み関数
として用いてフィルタを設計する。自己回帰モデルと
は、信号x(n)がその一つ前の信号x(n−1)によ
って決まるモデルのことをいう。一次の自己回帰モデル
は、
隣り合う成分の相関が強く、一次の自己回帰モデルで良
く近似できるという性質を利用して、一次の自己回帰モ
デルのパワースペクトルをフィルタ設計の際の重み関数
として用いてフィルタを設計する。自己回帰モデルと
は、信号x(n)がその一つ前の信号x(n−1)によ
って決まるモデルのことをいう。一次の自己回帰モデル
は、
【0014】
【数3】
【0015】と表せる。ここで、w(n)は白色雑音で
あり、ρはx(n)がx(n−1)にどの程度依存して
いるかを表す数で自己相関係数と呼ばれる。ρが大きけ
れば大きいほどx(n)はx(n−1)に近い値をとっ
ていることを意味している。式(2)の両辺をフーリエ
変換すると、
あり、ρはx(n)がx(n−1)にどの程度依存して
いるかを表す数で自己相関係数と呼ばれる。ρが大きけ
れば大きいほどx(n)はx(n−1)に近い値をとっ
ていることを意味している。式(2)の両辺をフーリエ
変換すると、
【0016】
【数4】
【0017】よって、
【0018】
【数5】
【0019】したがって、
【0020】
【数6】
【0021】となり、スペクトルはρが大きければ直流
付近に集中している。フィルタの設計にその性質を生か
すための重み関数は、
付近に集中している。フィルタの設計にその性質を生か
すための重み関数は、
【0022】
【数7】
【0023】となるが、画像拡大の際の0挿入された信
号はスペクトルが1/2に圧縮されている(図2
(b))ので、補完フィルタのためには、
号はスペクトルが1/2に圧縮されている(図2
(b))ので、補完フィルタのためには、
【0024】
【数8】
【0025】を重み関数として使用する。一般に、画像
をn倍に拡大する場合には、重み関数は以下のようにな
る。
をn倍に拡大する場合には、重み関数は以下のようにな
る。
【0026】
【数9】
【0027】ρ=0.9としたときの式(3)の重み関
数W(ω)の概形を図7に示す。フィルタの設計の際に
は、式(1)の理想フィルタの特性をD(ω)とすると
き、 W(ω)|D(ω)−P(ω)| の最大値を最小にするフィルタ特性P(ω)を与えるフ
ィルタのパラメータを決定する。これは周知のReme
zのアルゴリズムを実現するプログラムに式(3)また
は式(4)の重み関数を与えることにより計算すること
ができる。
数W(ω)の概形を図7に示す。フィルタの設計の際に
は、式(1)の理想フィルタの特性をD(ω)とすると
き、 W(ω)|D(ω)−P(ω)| の最大値を最小にするフィルタ特性P(ω)を与えるフ
ィルタのパラメータを決定する。これは周知のReme
zのアルゴリズムを実現するプログラムに式(3)また
は式(4)の重み関数を与えることにより計算すること
ができる。
【0028】なお、従来技術として述べた等リプル特性
のフィルタもこれと同様にして得られるが、重み関数W
が通過域と阻止域でそれぞれ一定値δ1 とδ2 で与えら
れ、本発明のようにそれぞれの領域で変化するものでは
ない。得られたフィルタの特性の一例を図8に示す。図
からも分かるように直流付近および折り返し周波数点付
近での周波数応答が良くなっているため、チェッカーボ
ード歪みは起こりにくい。
のフィルタもこれと同様にして得られるが、重み関数W
が通過域と阻止域でそれぞれ一定値δ1 とδ2 で与えら
れ、本発明のようにそれぞれの領域で変化するものでは
ない。得られたフィルタの特性の一例を図8に示す。図
からも分かるように直流付近および折り返し周波数点付
近での周波数応答が良くなっているため、チェッカーボ
ード歪みは起こりにくい。
【0029】画像処理の評価に良く使われるLenaの
画像について、図8に示すようなフィルタを用いて拡大
処理を行なったときの被験者17名による評価の結果を
表1に示し、比較のため、等リプル特性のフィルタを用
いたときの結果を表2に、最大平坦特性のフィルタを用
いたときの結果を表3に示す。
画像について、図8に示すようなフィルタを用いて拡大
処理を行なったときの被験者17名による評価の結果を
表1に示し、比較のため、等リプル特性のフィルタを用
いたときの結果を表2に、最大平坦特性のフィルタを用
いたときの結果を表3に示す。
【0030】
【表1】
【0031】
【表2】
【0032】
【表3】
【0033】表中、“ぼ”はぼやけが大きいことを示
し、“リ”はリンギングが認められることを示し、
“○”はそのいずれでもないことを示す。この結果か
ら、ぼけ、ぎざぎざ、リンギングのどれも認識されない
ためには、等リプル特性では23タップ、最大平坦特性
では33タップが必要であったが、本発明により設計さ
れたフィルタでは、16タップでよいことがわかり、画
質を損なうこと無く3割以上の演算量削減が可能となっ
た。少ないタップ数で良い結果が得られた16タップで
遷移域幅(図8中のwsとwpの間の幅)が0.2のフ
ィルタのフィルタ係数は以下の通りである。このフィル
タの周波数特性を図9に示す。
し、“リ”はリンギングが認められることを示し、
“○”はそのいずれでもないことを示す。この結果か
ら、ぼけ、ぎざぎざ、リンギングのどれも認識されない
ためには、等リプル特性では23タップ、最大平坦特性
では33タップが必要であったが、本発明により設計さ
れたフィルタでは、16タップでよいことがわかり、画
質を損なうこと無く3割以上の演算量削減が可能となっ
た。少ないタップ数で良い結果が得られた16タップで
遷移域幅(図8中のwsとwpの間の幅)が0.2のフ
ィルタのフィルタ係数は以下の通りである。このフィル
タの周波数特性を図9に示す。
【0034】 h(0)=−0.00399707030885 h(1)=−0.01937213597755 h(2)=0.02064040261517 h(3)=0.04575287378736 h(4)=−0.05172929806458 h(5)=−0.08493069857275 h(6)=0.14094533608183 h(7)=0.45269059043937 h(8)=0.45269059043937 h(9)=0.14094533608183 h(10)=−0.08493069857275 h(11)=−0.05172929806458 h(12)=0.04575287378736 h(13)=0.02064040261517 h(14)=−0.01937213597755 h(15)=−0.00399707030885
【0035】
【発明の効果】以上説明したように本発明によれば、少
ないタップ数で良好な結果を与える画像処理用ディジタ
ルフィルタが提供される。
ないタップ数で良好な結果を与える画像処理用ディジタ
ルフィルタが提供される。
【図1】画像拡大処理における時間(空間)領域での画
素値の分布の一例を示す図である。
素値の分布の一例を示す図である。
【図2】図1の画像の周波数スペクトルを示す図であ
る。
る。
【図3】線形補間演算のインパルス応答を示す図であ
る。
る。
【図4】線形補間演算の周波数特性を示す図である。
【図5】等リプル特性のフィルタの周波数特性を示す図
である。
である。
【図6】最大平坦特性のフィルタの周波数特性を示す図
である。
である。
【図7】本発明で用いられる重み関数の一例を示す図で
ある。
ある。
【図8】本発明の一実施例に係るディジタルフィルタの
特性の一例を示す図である。
特性の一例を示す図である。
【図9】少ないタップ数で良い結果が得られたフィルタ
の周波数特性を示す図である。
の周波数特性を示す図である。
Claims (4)
- 【請求項1】 周波数軸上の極値点における目標特性と
の差が周波数0に近づくにつれ徐々に小さくなる特性を
有するディジタルフィルタ。 - 【請求項2】 前記目標特性は、所定の通過域および阻
止域を有する特性である請求項1記載のディジタルフィ
ルタ。 - 【請求項3】 前記差が、画像信号の自己回帰モデルの
パワースペクトルの逆数に応じて変化する請求項1また
は2記載のディジタルフィルタ。 - 【請求項4】 前記自己回帰モデルのパワースペクトル
は、ρを0以上1以下の実数、ωを規格化された角周波
数とするとき、式 【数1】 で表わされる請求項3記載のディジタルフィルタ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000334603A JP2002141778A (ja) | 2000-11-01 | 2000-11-01 | ディジタルフィルタ |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000334603A JP2002141778A (ja) | 2000-11-01 | 2000-11-01 | ディジタルフィルタ |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2002141778A true JP2002141778A (ja) | 2002-05-17 |
Family
ID=18810489
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2000334603A Pending JP2002141778A (ja) | 2000-11-01 | 2000-11-01 | ディジタルフィルタ |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2002141778A (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006049213A1 (ja) * | 2004-11-02 | 2006-05-11 | Fujitsu Ten Limited | 映像信号処理方法、映像信号処理装置、及び、表示装置 |
JP2008079301A (ja) * | 2006-08-23 | 2008-04-03 | Matsushita Electric Ind Co Ltd | 撮像装置 |
US7956907B2 (en) | 2006-08-23 | 2011-06-07 | Panasonic Corporation | Image capture device with low-frequency noise suppression |
-
2000
- 2000-11-01 JP JP2000334603A patent/JP2002141778A/ja active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006049213A1 (ja) * | 2004-11-02 | 2006-05-11 | Fujitsu Ten Limited | 映像信号処理方法、映像信号処理装置、及び、表示装置 |
JP2008079301A (ja) * | 2006-08-23 | 2008-04-03 | Matsushita Electric Ind Co Ltd | 撮像装置 |
US7956907B2 (en) | 2006-08-23 | 2011-06-07 | Panasonic Corporation | Image capture device with low-frequency noise suppression |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4098319B2 (ja) | 解像度変換方法および解像度変換装置 | |
JP4534594B2 (ja) | 画像処理装置、画像処理方法、画像処理方法のプログラム及び画像処理方法のプログラムを記録した記録媒体 | |
EP1379074A2 (en) | Digital descreening of image signals | |
US7245785B2 (en) | Suppression of ringing artifacts during image resizing | |
Koschan et al. | A comparison of median filter techniques for noise removal in color images | |
KR20050041849A (ko) | 적응적 다상 필터에 근거한 이미지 보간 방법 및 장치 | |
Smolka et al. | Fast modified vector median filter | |
US20040001642A1 (en) | Control system for digital de-screening of documents | |
CN110517206B (zh) | 彩色摩尔纹消除方法及装置 | |
GB2547842A (en) | Image processing device and method, image pickup device, program, and recording medium | |
JP2003529294A (ja) | 画像をエンハンスするための電子回路及び方法 | |
WO2001097171A2 (en) | Magnifying digital image using edge mapping | |
Yang et al. | A novel fractional-order signal processing based edge detection method | |
JP2002141778A (ja) | ディジタルフィルタ | |
Kao et al. | Mltistage bilateral noise filtering and edge detection for color image enhancement | |
EP1195054A1 (en) | Video-apparatus with peaking filter | |
Rieder et al. | New concepts on denoising and sharpening of video signals | |
JP2671637B2 (ja) | ディジタル画像の拡大・縮小の方法およびその装置 | |
JP4030655B2 (ja) | 階調変換方法および装置 | |
JP2004524630A (ja) | 適応コンボリューション・フィルタを使用する画像処理用の方法および装置 | |
Beyerer et al. | Preprocessing and image enhancement | |
JP2014074777A (ja) | 合焦評価値生成装置、合焦評価値生成方法、及び、合焦評価値生成プログラム | |
KR100190039B1 (ko) | 고해상도 영상 확대방법 및 장치 | |
JPH11191861A (ja) | 画像処理装置及び画像処理システム | |
JP4596496B2 (ja) | ノイズ除去装置およびノイズ除去方法 |