JP3583831B2 - 信号処理装置 - Google Patents
信号処理装置 Download PDFInfo
- Publication number
- JP3583831B2 JP3583831B2 JP14277595A JP14277595A JP3583831B2 JP 3583831 B2 JP3583831 B2 JP 3583831B2 JP 14277595 A JP14277595 A JP 14277595A JP 14277595 A JP14277595 A JP 14277595A JP 3583831 B2 JP3583831 B2 JP 3583831B2
- Authority
- JP
- Japan
- Prior art keywords
- sampling
- data
- signal
- digital data
- equation
- 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
Images
Landscapes
- Image Processing (AREA)
- Editing Of Facsimile Originals (AREA)
- Facsimile Image Signal Circuits (AREA)
Description
【産業上の利用分野】
本発明は、標本化によって生じる折り返し成分を打ち消し、原信号が持つ高周波成分を復元し、また、デジタル信号を高解像度化する信号処理装置に関する。
【0002】
【従来の技術】
デジタル信号は、連続信号である原信号を標本化して作られる。標本化定理によれば、原信号が標本化周波数の半分の周波数(以下、ナイキスト周波数)以下に帯域制限されていれば、標本値から原信号を完全に復元することができ、また帯域制限されていなければ、完全な復元が不可能であることが知られている。
【0003】
デジタル信号を高解像度化する場合、高解像度化された出力信号の標本化周波数は、入力信号のそれよりも高いので、出力信号は、入力信号の表現できる帯域よりも広い帯域を持つ。つまり、出力信号の方が入力より多くの情報を表現できる。しかし、原信号に関する情報を得る手段を一つの入力デジタル信号だけに限定すれば、得られる情報は原信号(入力デジタル信号)のナイキスト周波数以下の帯域だけであり、それ以上の高周波成分を復元することは不可能である。
【0004】
この事実は、具体的にはキュービックスプライン補間や、周波数空間法による高解像度化では、高周波成分を含まない、ぼけた信号しか得られないという問題となって現れる。なお、キュービックスプライン補間とは、非格子点の濃度を、その周囲の16個の格子点における濃度を用いて、3次式によって補間するもので、その詳細は、例えば「コンピュータ画像処理入門」(総研出版)第4章に記載されている。また、周波数空間法としては、画像をDCT(2次元離散コサイン変換)によって空間領域から周波数領域に変換し、高周波成分として0を追加して、逆DCTによって空間領域に逆変換する画像の拡大方式がある(特開平2−76472号公報)。
【0005】
そこで、原信号に関するより多くの情報を総合して高解像度信号を得る方法として、一つの原信号を、標本化の位置(タイミング)を変えて同一の標本化間隔で標本化した複数組のデータを利用することを考える。図2は、同一の原信号について標本化位置を変えて同一の標本化間隔で標本化した、3組の標本化信号(1,2,3)を示す。
【0006】
標本化位置を標本化間隔の1/2だけずらして2回標本化すれば、1/2の標本化間隔で1回標本化したものと同等になることは明らかである。
【0007】
【発明が解決しようとする課題】
ところで、上記したようなデータを生成するために、例えばデジタルカメラを使い、複数回撮影することを考える。カメラの内部で1/2画素分、センサを移動させて2回撮影したとしても、丁度1/2画素だけずれたデータを得るためには、
(1)カメラ自身を完全に固定しなければならない
(2)撮影対象も完全に静止していなければならない
など、条件が厳しく、適用できる範囲が狭い。
【0008】
逆に言えば、カメラや撮影対象が移動する場合、移動量が不明であるが、とにかく異なる位置で標本化したデータが得られるので、このようなデータを統合することにより、原信号(撮影対象)の高周波成分の情報を的確に取り出すことができれば、その適用範囲は非常に広いと言える。
【0009】
本発明は、上記した事情を考慮してなされたもので、
本発明の目的は、標本化位置の異なる複数組のデジタルデータから、ナイキスト周波数以上の原信号の高周波成分を復元し、ぼけの少ない高解像度データを得るようにした信号処理装置を提供することにある。
【0010】
【課題を解決するための手段】
前記目的を達成するために、請求項1記載の発明では、標本化位置を変えて同一の標本化間隔で標本化された複数組のデジタルデータに対して、標本化周波数の半分の周波数以上の高周波成分を透過するフィルタ手段と、該フィルタ処理された複数組の信号をそれぞれ高解像度化する手段と、前記各デジタルデータの標本化位置に応じた重みを算出する手段と、対応する位置における各高解像度化データに該重みをつけて加重和をとり出力する手段とを備えたことを特徴としている。
【0011】
請求項2記載の発明では、前記標本化位置に応じた重みwkは、
【0012】
【数6】
【0013】
の関係を満たすことを特徴としている。
【0014】
ここで、
wkはk番目のデジタルデータのための重み
Δkはk番目のデジタルデータの標本化位置
2πi中のiは虚数単位。
【0015】
請求項3記載の発明では、前記標本化位置に応じた重みwkは、
【0016】
【数7】
【0017】
の関係を満たし、かつ
【0018】
【数8】
【0019】
を最小化することを特徴としている。
【0020】
請求項4記載の発明では、前記標本化位置に応じた重みwkは、
【0021】
【数9】
【0022】
の関係を満たし、かつ
【0023】
【数10】
【0024】
なる量の増加関数からなる評価関数を最小化することを特徴としている。
【0025】
請求項5記載の発明では、前記各デジタルデータの標本化位置の差を、与えられたデジタルデータから推定することを特徴としている。
【0026】
【作用】
周波数空間上で〔−2π,2π〕に帯域制限されている原信号を、3組の異なる標本化位置(その標本化間隔は同一)で標本化する。それをLPF処理すると、基本波は同一で、高調波は標本化位置のずれに従って、その位相がずれている。位相の異なるこれらの信号に、適当な重みをつけて平均することにより、高調波を打ち消し、基本波だけを残すことができる。
【0027】
【実施例】
以下、本発明の一実施例を図面を用いて具体的に説明する。
まず、本発明の原理について説明する。
〈折り返し成分復元の必要性〉
標本化定理より、理想的に帯域制限された原信号ならば、1組の標本信号から完全に復元されるので、標本化位置を変えて標本化しても得られる情報は増えない。つまり、異なる位置で標本化することで情報を増やすためには、標本化される原信号が理想的に帯域制限されておらず、折り返し歪が含まれていることが必要である。
【0028】
折り返された成分は、元々原信号が持っていたナイキスト周波数以上の成分であり、折り返しによって、より低い帯域と重なってしまったために歪となっているが、逆にその重なりを分離することができれば、原信号の高周波成分という重要な情報を得ることができる。
【0029】
本発明は、このような折り返しを含む複数の信号から、折り返し成分を取り除き、同時に原信号の高周波成分を復元するものである。
【0030】
〈解析〉
まず、異なる複数の位置で標本化した信号がどのような性質を持つかを解析する。ここでは、簡単のため一次元の原信号(連続信号)を例にして説明する。
【0031】
原信号をf(x)、それを間隔1で標本化した標本化信号をd(x)とする。また、フーリエ変換、および畳み込みを式(1)、(2)、(3)のように定義する。
【0032】
【数11】
【0033】
このとき、式(4)、(5)に示す性質を持つことが知られている。
【0034】
【数12】
【0035】
標本化は、原信号と、以下のようなδパルス列(s(x))との積をとる操作と考えられる。
【0036】
【数13】
【0037】
ここで、S(ω)はs(x)のフーリエ変換で、周波数空間での表現である。以下、大文字の英字が周波数空間での表現であり、小文字の英字が実空間での表現である。
【0038】
そうすると、標本化したデジタル信号d(x)およびそのフーリエ変換D(ω)は、式(7)、(8)で表される。
【0039】
【数14】
【0040】
図3は、実空間での標本化の過程を説明する図であり、(a)はf(x)、(b)はs(x)、(c)はf(x)s(x)である。ここで、f(x)が画像信号であるとき、s(x)が主走査信号となって、画像信号がサンプルされ(s(x)の周期が画素間隔となる)、その振幅値(f(x)s(x))が量子化される。
【0041】
図4は、実空間をフーリエ変換した周波数空間における標本化の過程を説明する図であり、(a)はf(x)のスペクトル、(b)はs(x)のスペクトル、(c)はd(x)のスペクトルであり、原信号f(x)のフーリエ変換F(ω)を2π間隔で繰り返し並べた形になっている(標本化周波数=2π、ナイキスト周波数=π)。以下、式(8)のk=0の項を基本波、k=mの項をm次高調波と呼ぶ。
【0042】
次に、同じf(x)について、標本化位置をΔだけずらして同一の標本化間隔で標本化する場合を考える。この標本化した信号をdΔ(x)とする。これは、式(6)のδパルス列s(x)を、Δだけずらしたものs(x−Δ)に変えることに相当する。位置をずらしたパルス列のフーリエ変換は、前述した式(4)の性質から
【0043】
【数15】
【0044】
となる。
【0045】
これは、式(6)にexp(−iωΔ)を掛けた形である。つまり、S(ω)をωに比例して位相を変化させた、「ねじれた」形になる(図6)。図5は、標本化位置をずらした場合の実空間における標本化の過程を示し、図6は、その周波数空間における標本化の過程を示す。S(w)は複素ベクトル(Reが実軸、Imが虚軸)であり、そのベクトルの長さは変わらずに、左右にいくに従って、ベクトルの位相がずれていく(図6(b)では、このずれを、破線の円周上でベクトルが回転している様子で表している)。
【0046】
そうすると、このずれた位置で標本化した信号(dΔ(x)=f(x)s(x−Δ))のフーリエ変換は、
【0047】
【数16】
【0048】
となり、やはり「ねじれた」形になる。図6(c)において、2πの位置のベクトルは、F(ω)を2πだけ位相をずらし(F(ω−2π))、実軸から2πΔだけ回転(exp(−2πΔ))したものであり、−2πの位置のベクトルは、F(ω)を−2πだけ位相をずらし(F(ω+2π))、実軸から2πΔだけ回転(exp(−2πΔ))したものである。
【0049】
このように、折り返し歪を含み、標本化位置の異なる複数組のデジタルデータを周波数空間で見ると、基本波は標本化位置のずれ(Δ)によらず原信号と同一で、k次高調波は2πkΔだけ位相がずれていることが分かる。この位相のずれを使って、以下説明するように、折り返された成分だけを取り除き、本来の高周波成分を復元することができる。
【0050】
〈高周波成分の復元〉
以上の解析結果を使い、一つの原信号を異なる位置で標本化した複数の組の標本化信号から、折り返し成分を取り除き、原信号の高周波成分を復元する原理を説明する。
【0051】
ここでは、原信号f(x)が、周波数空間上で、ナイキスト周波数の2倍、つまり〔−2π,2π〕に帯域制限されているものと仮定する。この場合、図7(a)に示すように、周波数空間では、隣接する信号分布の裾が重なり、折り返し成分の重なりは2重である。つまり、ある周波数に存在する成分は、基本波、高調波の内、たかだか2つを足したものとなる。
【0052】
また、ここでは標本化信号の標本化位置は既知であると仮定し、3つの位置Δ0、Δ1、Δ2で標本化した標本化信号をそれぞれd0(x)、d1(x)、d2(x)とし、またそれぞれのフーリエ変換をD0(ω)、D1(ω)、D2(ω)とする。
【0053】
まず、一つの標本化信号d0(x)にローパスフィルタLPFをかけ、〔−2π,2π〕の帯域だけを取り出す。この信号をe0(x)とする。そのフーリエ変換(E0(ω))は次のようになる。
【0054】
【数17】
【0055】
これは、基本波と±1次高調波の半分づつが重なった信号である(図7(a)、図8)。
【0056】
同様に、標本化信号d1(x)、d2(x)もそれぞれLPF処理し、e1(x)、e2(x)を求める。これら3つの信号は、式(10)より、基本波は同一で、高調波は標本化位置(Δ0,Δ1,Δ2)に従い、その位相がずれている(図9)。
【0057】
図9は、周波数空間上で、標本化位置のずれによる位相のずれを説明する図である。3つの基本波の位相は実軸にあり、信号E0(ω)の1次高調波の位相は、実軸から2πΔ0だけ位相がずれていて、信号E1(ω)の1次高調波の位相は、実軸から2πΔ1だけ位相がずれていて、信号E2(ω)の1次高調波の位相は、実軸から2πΔ2だけ位相がずれている。
【0058】
次に、これら3つの信号を、適当な重みをつけて加重和をとる。3つの信号を重みwkで加重和をとった信号、
【0059】
【数18】
【0060】
は、周波数空間で表現すれば、式(12)より、
【0061】
【数19】
【0062】
となる。よって、
【0063】
【数20】
【0064】
なる関係を満たす重み(wk)を使えば、式(15)より、高調波成分は打ち消され、従って折り返し成分が打ち消され、同時に式(14)より基本波(F(ω))はそのまま残り、結果として原信号(F(ω))を復元することができる。図7(c)は、式14、15によって打ち消された高調波成分を示し、図7(d)は、復元された基本波を示す。式(14)、(15)は、重みwkについての連立1次方程式なので、行列演算によって簡単に解を求めることができる。
【0065】
ここでは、原信号は〔−2π,2π〕に帯域制限されているものと仮定した。そのため折り返しは2重となり、一つの高調波成分を打ち消すため式(15)の条件が一つ必要であった。より広い帯域の信号を扱う場合、打ち消すべき折り返し成分が増えるため、式(15)に相当する条件式が増え、打ち消すために必要なデジタル信号の数が増えることになる。例えば、原信号が〔−(M+1)π,(M+1)π〕の帯域を持つ場合、原信号を復元するためには、M個の折り返し成分を打ち消す必要があるため、式(15)に相当する打ち消し条件式が2M個となる。重み総和条件と併せて2M+1個の連立方程式を解くためには、2M+1個のデジタル信号が必要となる。
【0066】
〈雑音への対応〉
上記した説明では、各信号には雑音は含まれないものと考え、原信号の完全な復元を行うものであった。しかし、現実には信号に必ず雑音が含まれるので、雑音に強い信号処理方法が求められる。
【0067】
ここでは、まず前述したような加重和による信号の雑音の現れ方を解析し、それに対する2つの対処方法を説明する。
【0068】
一般に、複数の信号が、互いに無相関で分散σ2の雑音を含むとすると、重みwkを使ったそれらの加重平均信号の雑音は、
【0069】
【数21】
【0070】
なる分散を持つ。つまり、加重平均信号の雑音は、それぞれの信号のΣ(wk)2倍になると考えられる。
【0071】
式(14)、(15)で求められる重み(wk)は、雑音を考慮していないので、標本化位置の組み合わせによっては、Σ(wk)2 が非常に大きくなる場合があり、得られる加重平均信号は折り返し成分が打ち消されると同時に雑音が強調される可能性がある。そこで、雑音の低減も考慮した重みを考える。
【0072】
〈「折り返し成分を打ち消し」かつ「雑音最小化」のための重み〉
3つの入力信号を使う場合は、折り返し成分を打ち消すような加重(wk)は、式(15)の条件から一意に決まる。一方、4つ以上の入力信号を使えば、自由度が増えるので、式(15)を満たした上で、残る自由度を使って式(16)を最小化する重みをとることにより、折り返し成分を打ち消した上で、更に雑音の低減を図ることができる。
【0073】
つまり、式(14)、(15)を満たす重みの代わりに、
【0074】
【数22】
【0075】
を最小化する重みを使う。
【0076】
これはwkについて、1次式の条件付き、2乗和の最小化であるので、公知の計算法(Lagrangeの未定乗数法)により、連立1次方程式を解くことで簡単に解が得られる。
【0077】
〈「折り返し成分および雑音」の最小化のための重み〉
雑音成分との兼ね合いを考えると、折り返し成分も一種の雑音であり、必ずしも0にする必要はない。折り返し成分および雑音成分をそれぞれ低減する重みも有効である。
【0078】
これは、式(15)を束縛条件とはせず、次式のような、その絶対値の2乗値
【0079】
【数23】
【0080】
を式(16)と同等の雑音評価値と考え、それぞれの雑音評価値に対して増加関数となる評価関数をとり、それを最小化する重みを採用することによって実現できる。
【0081】
例えば、Σwk=1の条件のもとで、次式のような、式(18)と式(16)の加重和を評価関数とし、
【0082】
【数24】
【0083】
これを最小化する重みを使う。この場合、α≒0では折り返し成分の打ち消しが優先となり、α≒1では基本波の雑音低減が優先となる。式(19)も、wkについての2乗和の最小化であり、やはり、連立1次方程式を解くことで簡単に解が得られる。この場合は、2組以上の入力データが必要である。
【0084】
〈標本化位置の推定〉
上記した説明では、それぞれの入力標本化信号について、標本化位置は既知であるとして考えてきた。しかし、前述したようなデジタルカメラを動かしながら撮影したデータを扱う場合などにおいては、標本化位置を予め知ることができないので、与えられたデジタルデータからそれぞれの標本化位置を推定する必要がある。その推定方法を以下、説明する。
【0085】
2つのデジタルデータの標本化位置のずれは、移動する対象物の速度とも考えることができる。そこで、速度vで平行移動する1次元の連続関数f(x,t)を例に考える。変化が微小ならば、Tayler展開の1次近似として、
【0086】
【数25】
【0087】
が各点で成り立つ。この式の∂f/∂xおよび∂f/∂tから、各点で変化速度vを求めることができる。
【0088】
しかし、上記した式では、
(1)1次近似からの外れ
(2)雑音の影響
などのため精度が悪いので、次式のような、式(20)の左辺の2乗をある領域内で積分した値を評価関数とし、この評価関数が最小値をとるvを求めることで、正確な値が得られる。
【0089】
【数26】
【0090】
このようにして、2つのデータ間の相対位置が分かれば、1つのデータを基準として、他の全てのデータとの相対位置を求めることにより、全てのデータの相対位置を推定することができる。
【0091】
〈実施例〉
図1は、本発明の実施例のブロック構成図である。本実施例では、前述した本発明の原理に基づき、同一の対象を異なる複数の標本化位置で標本化したデジタルデータを入力し、標本化される前の原信号の高周波成分を復元した、高解像度のデジタルデータを得るための信号処理方法を説明する。以下、各部分について、その機能を説明する。
【0092】
(入力デジタルデータ)
本発明で使用する入力デジタルデータ1は、同一の対象を、位置をずらして同一の標本化間隔で標本化した3組のデジタルデータである。入力デジタルデータとしては、例えば以下のようにして作成されたものが使用される。
【0093】
(a)イメージスキャナによって、一枚の原稿を、少しずつ原稿の設置位置をずらして、3回スキャンして得られる3枚のデジタル画像データ。
(b)デジタルカメラの位置、または方向を少しずつずらしながら、ほぼ静止した対象を3回撮影することによって得られる3枚のデジタル画像データ。
(c)ビデオカメラの位置、または方向を少しずつずらしながら、ほぼ静止した対象を撮影し、そのビデオデータの連続する3コマ分をデジタル化したデータ。
【0094】
一般的に、上記したスキャナの如きデジタル画像入力装置の内部では、図10に示すような信号処理が行われている。図10は、デジタル画像入力装置の信号処理を説明する図である。すなわち、
1.濃度分布orig(x)を持つ原稿に対し、
2.光学系やセンサの開口特性などによる伝達関数(点像広がり関数、psf(x))が作用し、
3.その後、式(6)のサンプルパルス列(s(x))を掛けることにより、標本化される。
【0095】
本発明の原理で説明した復元すべき原信号(f(x))とは、量子化される直前の信号である。この場合は、原稿濃度分布に点像広がり関数が作用した信号、orig(x)*psf(x)がその原信号に相当する。
【0096】
さて、このようなデジタル画像入力装置を使用した場合の原信号の特性について説明する。一般にデジタル画像入力装置においては、上記した点像広がり関数が適切なLPF特性を持つように設計することによって、折り返し歪を防止するための帯域制限を実施している。しかし、実際の帯域制限は、
理想的なLPFの特性(周波数空間で矩形関数)は実現が困難である、
エッジ再現性のため高周波成分も保存したい、
などの理由から不完全であり、得られるデジタルデータは、通常ある程度の折り返し歪を含んでいる。このような不完全に帯域制限された原信号(f(x))は、原理の説明で述べたような、〔−2π,2π〕に帯域制限した信号と考えられ、本発明の信号処理によって原信号の高周波成分を復元することができる。図11は、帯域制限フィルタの周波数特性を示し、21は〔−π,π〕で1、それ以外では0となる理想的な帯域制限フィルタの特性であり、22は実際のpsfによる帯域制限フィルタの特性であり、23は本発明の帯域制限フィルタの特性である。
【0097】
上記した入力デジタルデータの例は、2次元のデジタルデータであるが、以下では簡単のため、1次元のデータを処理する方法を説明する。以下の説明から2次元データ処理への拡張も容易である。また、2次元の画像データでも、その標本化位置のずれを横方向に限定できれば、各列毎に独立した1次元のデータとして扱うこともできる。
【0098】
本実施例では、3組のデジタルデータを使用する。以下、3組の入力デジタルデータをそれぞれin0,in1,in2とする。また、第i組の入力デジタルデータ(ini)の、第j番目のサンプル(画素)を、ini〔j〕とする。
【0099】
(基準データバッファ)
基準データバッファ2は、標本化位置の基準の位置とするため、第0組のデータ(in0)を格納するバッファである。後段の位置差推定部では、2組のデジタルデータを入力し、それらの標本化位置の差を推定する。第0組のデータ(in0)を基準データとし、in0とin1、およびin0とin2の標本化位置差を推定することにより、第0組のデータの標本化位置を基準として、3組のデータの相対的な標本化位置を推定することができる。
【0100】
(位置差推定部)
位置差推定部3は、それぞれの入力データ(in1,in2)について、基準データバッファに格納された基準データ(in0)との、標本化位置の差を推定することにより、3組の入力データの相対的な標本化位置を推定する手段である。
【0101】
第1組のデータの標本化位置をΔ1、第2組のデータの標本化位置をΔ2として出力する(基準となる第0組のデータの標本化位置は0である)。このような推定は、式(21)の微分、積分をデジタルデータの差分、加算に置き換えることにより実現できる。
【0102】
すなわち、in0とiniの位置差(Δi)を推定するには、式(21)を
【0103】
【数27】
【0104】
と置き換え、また、対象の速度(v)ではなく、標本化位置の移動量を測るので、符号を変えて、
【0105】
【数28】
【0106】
を最小化するΔiとして求められる。これは線形最小2乗法により解くことができる。
【0107】
(加重計算部)
加重計算部4は、位置差推定部3で得られた標本化位置(0,Δ1,Δ2)に基づいて、後段の積和部で使用される加重を求める手段である。第0,1,2組の入力デジタルデータに対する加重として、それぞれw0,w1,w2を出力する。具体的には、式(14)、(15)を解く演算を行う。
【0108】
また、他の方法として、4組以上の入力デジタルデータを使い、式(14)、(15)の束縛条件のもとで、式(17)を最小化する重みを求めることもできる。さらに、他の方法として、2組以上の入力デジタルデータを使い、式(14)の束縛条件のもとで、式(19)を最小化する重みを求めることもできる。
【0109】
(解像度倍率指定部)
本実施例の最終的な出力デジタルデータは、入力デジタルデータよりも画素数の多い高解像度なデータである。解像度倍率指定部5は、出力デジタルデータの解像度を、入力デジタルデータの解像度の倍率として指定する手段である。ここでは、解像度を10倍にするものとする。つまり、出力デジタルデータは入力デジタルデータに比べ、サンプル数が10倍、標本化周期が1/10となる。
【0110】
解像度の倍率は、利用者が指定するようにしてもよい。ただし、本実施例では、ナイキスト周波数の2倍の周波数成分までを復元するので、出力データに折り返し歪を含めないためには、倍率は2倍以上としなければならない。
【0111】
(広帯域LPF処理および高解像度化処理部)
この部分は、本発明の原理で説明した「〔−2π,2π〕の透過帯域を持つLPFにより、基本波全部および1次高調波の半分を足した信号を取り出す」部分に相当する。
【0112】
入力デジタルデータは、ini〔j〕というデータ列であるが、これを次式のように、位置(x)に対する連続関数(デルタパルス列)と考えたものが式(8)の標本化信号d(x)にあたる。
【0113】
【数29】
【0114】
この信号(di(x))に、式(11)の広帯域LPFを掛ける。
【0115】
理想的なLPF処理は、周波数空間では、
【0116】
【数30】
【0117】
なる矩形関数との乗算であり、実空間では、sin(2x)/2x(sinc関数)とのコンボリューションである。
【0118】
di(x)をLPF処理した信号をei(x)とすると、
【0119】
【数31】
【0120】
となる。
【0121】
さらに、式(23)より、入力デジタルデータ(ini〔j〕)を使って書き換えると、
【0122】
【数32】
【0123】
となる。
【0124】
このLPF処理によって得られる信号(ei(x))は連続信号であり、これを標本化した信号を、このブロック6の出力デジタル信号とする。そのため、連続信号の内、欲しい(標本化される)位置(x)についてだけ式(24)の計算をすればよい。
【0125】
このようにして「広帯域LPF処理および高解像度化処理部6」で得られるデジタルデータをmidi〔j〕とする。ここでは新たな標本化間隔を1/10としたので、midi〔j〕は結局、
【0126】
【数33】
【0127】
となり、予めsin(2(j/10−k))/(2(j/10−k))なるテーブル(実際には有限の範囲内でそれを近似したもの)を用意しておけば、入力デジタル信号iniの積和によって計算できる。
【0128】
LPF処理された信号ei(x)は、原信号のナイキスト周波数よりも高い周波数成分を含む信号である。そしてそれを有効に表現するためには、つまり原信号の高周波成分を復元するためには、LPF処理された信号ei(x)に対するサンプリング間隔は、入力信号のサンプリング間隔より狭い間隔でなければならない。信号ei(x)に対するサンプリング間隔が、入力信号のサンプリング間隔より狭くなければ、再び折り返し歪が発生し、意味がない。そのため、ここでは1/10としている。
【0129】
(中間データバッファ)
中間データバッファ7は、高解像度化処理された中間デジタルデータ(midi〔j〕)を保存するバッファである。
【0130】
(積和部)
この積和部7での加重和計算は、本発明の原理で説明した「式(13)の加重和をとり、周波数空間で3本のベクトルで表される高調波成分を打ち消す」部分に相当する。本実施例の構成では、実空間で処理する。実空間でもやはり各中間データに、それぞれ加重wiをかけた後、加算すればよい。
【0131】
ここで、注意しなければならないことは、原理の説明では、原信号の位置は固定されていて、サンプリング位置の方が移動していることである。一方、中間デジタルデータは、各入力デジタルデータの標本化位置のずれに係らず、入力デジタルデータだけから計算されるので、mid0〔0〕(第0組の中間データの0番目のサンプル)は、in0〔0〕(第0組の入力データの0番目のサンプル)と同じ位置というだけで、mid1〔0〕(第1組の中間データの0番目のサンプル)と同じ位置にあたる訳ではない。
【0132】
このように、入力デジタル信号の各サンプルは、それぞれ異なる位置に対応しているので、同じ位置、つまり図12の縦方向に並んだ中間データ同志を足し合わせるためには、加算するときに各中間データバッファの1番目のサンプル同志ではなく、位置をずらせて足し合わせなければならない。
【0133】
図12は、中間データ間の位置関係を説明する図である。入力データd0(x)は、サンプリング間隔(=1)でサンプルされる。すなわち、図の太い矢印で示すように、サンプル点はj=0、次のサンプル点はj=10である。入力データd1(x)は、位置差推定部で推定されたΔ1だけずれた位置から、同様に、サンプリング間隔(=1)でサンプルされ、そのサンプル点はj=0、次はj=10である。入力データd2(x)は、位置差推定部で推定されたΔ2だけずれた位置から、同様に、サンプリング間隔(=1)でサンプルされ、そのサンプル点はj=0、次はj=10である。
【0134】
eo(x),e1(x),e2(x)は、それぞれdi(x)をLPF処理した連続信号であり、ei(x)を、サンプリング間隔(=1/10)でサンプル(ei(x)上の各点がサンプル点)した値が、それぞれmidi〔j〕となる。中間データバッファmid0には、例えばサンプル点j=0を開始点として、サンプリング間隔(=1/10)でサンプルしたデータが、図13に示すように格納される。同様に、中間データバッファmid1,2には、サンプル点j=0を開始点として、サンプリング間隔(=1/10)でサンプルしたデータが格納される。図13では、位置差を付けて格納しているが、実際には各バッファにはサンプル点j=0から順に、その値が格納される。
【0135】
例えば、「位置差推定部」で、Δ1=1/10、Δ2=3/10と推定されたとする。ここでは、中間データは入力データの10倍の密度であるので、この位置差は、積和部8で10倍に変換される。従って、第1番目の中間データmid1と第0番目の中間データmid0は、1サンプル分のずれとなり、第2番目の中間データmid1と第0番目の中間データmid0は、3サンプル分のずれとなる。
【0136】
そして、図12の縦方向に並んだ中間データ同志を足し合わせるとき、各データはサンプル点jを基点としているので、例えばサンプル点j=0で中間データ同志を足し合わせる場合は、第0番目の中間データはmid0〔0〕であり、第1番目の中間データは、1サンプル前である(−1)番目のサンプルのmid0〔−1〕となり、第2番目の中間データは、3サンプル前である(−3)番目のサンプルのmid0〔−3〕となる。
【0137】
従って、積和部8は、図13に示すように位置をずらして、中間データバッファ7から、各データmid0〔j〕、mid0〔j−1〕、mid0〔j−3〕を読み出し、加重計算部からの重みw0,w1,w2を掛けて足し合わせる。以下、同様にして、j+1、j+2...の各サンプル点で出力を算出する。
【0138】
このように、推定された標本化位置(Δi)と解像度倍率に従って、中間データバッファから位置をずらして各サンプル点のデータを取り出し、加重和をとる。出力データ(out〔j〕)は、基準データ(in0)と同じ位置とすると、最終的に出力されるデータ(out〔j〕)は、
【0139】
【数34】
【0140】
となる。
【0141】
上記した実施例では、簡単のため1次元の信号について説明したが、多次元信号についても同様に説明される。ただし2次元信号の場合、折り返しが縦横2方向から来るので、原信号の帯域が縦横共に2倍に広がると、3つの折り返し成分が重なり、それらを打ち消すためには2M+1=7個のデジタルデータを必要とする。
【0142】
また、上記した実施例は、複数のデジタルデータを使い、1つの高解像度化データを得る方法であるが、動画像を入力として、例えば、フレーム1からフレーム9の9枚のデータを使ってフレーム5を高解像度化し、続いてフレーム2からフレーム10を使ってフレーム6を高解像度化し、と続けることにより、動画像から動画像への高解像度化を実現することも可能となる。
【0143】
【発明の効果】
以上、説明したように、本発明によれば、標本化間隔が同一で標本化位置を変えた複数組のデジタルデータを用いることにより、標本化による折り返し歪を打ち消し、原信号の高周波成分を復元することができる。これにより、ぼけ、歪の少ない高解像度のデータを得ることが可能となる。
【0144】
本発明によれば、標本化間隔が同一で標本化位置を変えた複数組のデジタルデータをフィルタ処理し、フィルタ処理後の信号を高解像度化処理し、標本化位置に応じた重みで積和演算を行っているので、単純な演算によって折り返し成分を打ち消して、原信号の高周波成分を復元することができる。
【0145】
本発明によれば、雑音の少ないデータに対して、高精度に原信号の高周波成分を復元することができる。
【0146】
本発明によれば、雑音を含むデータに対しても、折り返し成分と雑音を同時に低減することができる。
【0147】
本発明によれば、標本化位置が未知の場合に、入力されたデジタルデータから標本化位置を推定しているので、未知の標本化位置を持つデジタルデータに対しても、原信号の高周波成分を復元することができる。
【図面の簡単な説明】
【図1】本発明の実施例のブロック構成図である。
【図2】同一の原信号について標本化位置を変えて標本化した、3組の標本化信号を示す。
【図3】(a)、(b)、(c)は、実空間での標本化の過程を説明する図である。
【図4】(a)、(b)、(c)は、周波数空間における標本化の過程を説明する図である。
【図5】(a)、(b)、(c)は、標本化位置をずらした場合の実空間における標本化の過程を示す。
【図6】(a)、(b)、(c)は、標本化位置をずらした場合の周波数空間における標本化の過程を示す。
【図7】(a)、(b)、(c)、(d)は、本発明による高周波成分の復元を説明する図である。
【図8】LPF処理された信号を示す図である。
【図9】周波数空間上で、標本化位置のずれによる位相のずれを説明する図である。
【図10】デジタル画像入力装置の信号処理を説明する図である。
【図11】帯域制限フィルタの周波数特性を示
【図12】中間データ間の位置関係を説明する図である。
【図13】中間データバッファからのデータの取り出し方を説明する図である。
【符号の説明】
1 入力デジタルデータ
2 基準データバッファ
3 位置差推定部
4 加重計算部
5 出力解像度指定部
6 広帯域LPF処理および高解像度化処理部
7 中間データバッファ
8 積和部
9 出力デジタルデータ
8 出力デジタルデータ
Claims (5)
- 標本化位置を変えて同一の標本化間隔で標本化された複数組のデジタルデータに対して、標本化周波数の半分の周波数以上の高周波成分を透過するフィルタ手段と、該フィルタ処理された複数組の信号をそれぞれ高解像度化する手段と、前記各デジタルデータの標本化位置に応じた重みを算出する手段と、対応する位置における各高解像度化データに該重みをつけて加重和をとり出力する手段とを備えたことを特徴とする信号処理装置。
- 前記各デジタルデータの標本化位置の差を、与えられたデジタルデータから推定することを特徴とする請求項1記載の信号処理装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP14277595A JP3583831B2 (ja) | 1995-06-09 | 1995-06-09 | 信号処理装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP14277595A JP3583831B2 (ja) | 1995-06-09 | 1995-06-09 | 信号処理装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH08336046A JPH08336046A (ja) | 1996-12-17 |
JP3583831B2 true JP3583831B2 (ja) | 2004-11-04 |
Family
ID=15323308
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP14277595A Expired - Lifetime JP3583831B2 (ja) | 1995-06-09 | 1995-06-09 | 信号処理装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3583831B2 (ja) |
Families Citing this family (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4617825B2 (ja) * | 2003-10-22 | 2011-01-26 | ソニー株式会社 | データ処理装置およびデータ処理方法、プログラム、並びに記録媒体 |
JP2008243037A (ja) * | 2007-03-28 | 2008-10-09 | National Univ Corp Shizuoka Univ | 画像処理装置、画像処理方法および画像処理プログラム |
JP4988460B2 (ja) * | 2007-07-05 | 2012-08-01 | 株式会社日立製作所 | 画像信号処理装置および画像信号処理方法 |
EP2012532A3 (en) | 2007-07-05 | 2012-02-15 | Hitachi Ltd. | Video displaying apparatus, video signal processing apparatus and video signal processing method |
JP4793339B2 (ja) | 2007-07-09 | 2011-10-12 | ソニー株式会社 | 動画像変換装置、動画像復元装置、および方法、並びにコンピュータ・プログラム |
JP4876048B2 (ja) | 2007-09-21 | 2012-02-15 | 株式会社日立製作所 | 映像送受信方法、受信装置、映像蓄積装置 |
JP5250232B2 (ja) * | 2007-10-04 | 2013-07-31 | 株式会社日立製作所 | 画像信号処理装置、画像信号処理方法及び画像表示装置 |
JP5250233B2 (ja) * | 2007-10-04 | 2013-07-31 | 株式会社日立製作所 | 画像信号処理装置、画像信号処理方法及び画像表示装置 |
US8089557B2 (en) | 2007-10-04 | 2012-01-03 | Hitachi, Ltd. | Video signal processing apparatus, video signal processing method and video display apparatus |
JP5291916B2 (ja) * | 2007-11-06 | 2013-09-18 | 株式会社日立製作所 | 動画表示装置 |
JP5416899B2 (ja) * | 2007-11-22 | 2014-02-12 | 株式会社日立製作所 | 映像表示装置 |
JP4508279B2 (ja) * | 2008-07-17 | 2010-07-21 | ソニー株式会社 | 画像処理装置、画像処理方法、及び、プログラム |
JP5369526B2 (ja) * | 2008-07-28 | 2013-12-18 | 株式会社日立製作所 | 画像信号処理装置、表示装置、録画再生装置、画像信号処理方法 |
US8212927B2 (en) | 2008-07-28 | 2012-07-03 | Hitachi, Ltd. | Image signal processing apparatus, image signal processing method and video display apparatus |
KR102090499B1 (ko) | 2008-12-19 | 2020-03-18 | 가부시키가이샤 한도오따이 에네루기 켄큐쇼 | 표시장치의 구동방법 |
JP5146335B2 (ja) * | 2009-01-22 | 2013-02-20 | ソニー株式会社 | 画像処理装置および方法、並びにプログラム |
KR102104986B1 (ko) | 2009-02-06 | 2020-04-27 | 가부시키가이샤 한도오따이 에네루기 켄큐쇼 | 표시 장치 구동 방법 |
JP5448981B2 (ja) | 2009-04-08 | 2014-03-19 | 株式会社半導体エネルギー研究所 | 液晶表示装置の駆動方法 |
JP5251825B2 (ja) | 2009-10-22 | 2013-07-31 | セイコーエプソン株式会社 | 画像表示装置および画像表示方法 |
US20110157470A1 (en) * | 2009-12-28 | 2011-06-30 | Sadao Tsuruga | Receiver, receiving method and output control method |
US8665367B2 (en) | 2010-01-29 | 2014-03-04 | Sharp Laboratories Of America, Inc. | Video resolution enhancement technique |
JP2011186777A (ja) * | 2010-03-09 | 2011-09-22 | Toshiba Corp | バーコード読取装置およびバーコード読取方法 |
JP5644626B2 (ja) * | 2011-03-28 | 2014-12-24 | 株式会社Jvcケンウッド | 信号処理装置、信号処理方法、および信号処理プログラム |
JP5696560B2 (ja) * | 2011-03-28 | 2015-04-08 | 株式会社Jvcケンウッド | 信号処理装置、信号処理方法、および信号処理プログラム |
-
1995
- 1995-06-09 JP JP14277595A patent/JP3583831B2/ja not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
JPH08336046A (ja) | 1996-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3583831B2 (ja) | 信号処理装置 | |
KR100648308B1 (ko) | 해상도 변환방법 및 장치 | |
Ramponi | Warped distance for space-variant linear image interpolation | |
US8098961B2 (en) | Resolution enhancement apparatus and method | |
Katkovnik et al. | A spatially adaptive nonparametric regression image deblurring | |
Sekko et al. | A deconvolution technique using optimal Wiener filtering and regularization | |
JP4582666B2 (ja) | 画像を補間および補正する方法および装置ならびに該装置を有するスキャナ | |
JP3773563B2 (ja) | 2次元デジタル信号処理装置 | |
Jency Rubia et al. | Digital image restoration using modified richardson-lucy deconvolution algorithm | |
Reichenbach et al. | Restoration and reconstruction of AVHRR images | |
JP4004561B2 (ja) | 画像処理方法および装置 | |
JP3234232B2 (ja) | 超空間的可変アポダイゼーション(スーパーsva) | |
JP6532151B2 (ja) | 超解像装置およびプログラム | |
Saito et al. | Super-resolution interpolation with a quasi blur-hypothesis | |
Choi et al. | Color image interpolation in the DCT domain using a wavelet-based differential value | |
Boggione et al. | Restoration of Landsat-7 images | |
Schowengerdt et al. | Topics in the two-dimensional sampling and reconstruction of images | |
Cheikh et al. | Color image interpolation using vector rational filters | |
Jung et al. | Kernel design for real-time denoising implementation in low-resolution images | |
El-Khamy et al. | Adaptive least squares acquisition of high resolution images | |
Rathod et al. | Satellite image resolution enhancement using SWT and DWT with SWT | |
Wood et al. | Restoration for sampled imaging systems | |
Gurjar et al. | Study of various resampling techniques for high-resolution remote sensing imagery | |
JP6532148B2 (ja) | 超解像装置およびプログラム | |
Ruikar et al. | Image Denoising using Tri Nonlinear and Nearest Neighbour Interpolation with Wavelet Transform |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20040518 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20040705 |
|
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: 20040727 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20040730 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080806 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080806 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090806 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090806 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100806 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100806 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110806 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110806 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120806 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120806 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130806 Year of fee payment: 9 |
|
EXPY | Cancellation because of completion of term |