(57)【要約】  本発明は、ディジタル画像を補間および補正する方法および装置に関する。ソース画像は複数のフィルタ計数を有するＦＩＲフィルタを用いて目標画像にマッピングされる。フィルタ係数は画像を補間するための情報も、画像を補正するための情報も包含するので、単一の処理ステップで補間も補正も実施される。ソース画像を目標画像にマッピングする本発明による方法をハードウェア回路により簡単に実現することができるので、装置、殊にスキャナにおけるマッピングは全体のプロセス経過に遅延が生じることなく実施される。
2007538436.tif 000002
【特許請求の範囲】
【請求項１】
ディジタル画像を補間および補正する方法において、
複数のフィルタ係数を有するＦＩＲフィルタを用いてソース画像を目標画像にマッピングし、
前記フィルタ係数は前記画像を補間するための情報も、前記画像を補正するための情報も包含し、
単一のステップで前記画像の補間も前記画像の補正も実施することを特徴とする、ディジタル画像を補間および補正する方法。
【請求項２】
前記ソース画像を前記目標画像にマッピングする前に、ｋ個の離散的な位置周波数を用いて次式にしたがいｎ個のフィルタ係数を計算し、
【数１】
2007538436.tif 000015
ここでａ1〜ａnは前記フィルタ係数であり、Ｈ（ωe）はそれぞれの位置周波数ωiに関する位相要素と乗算された伝達関数である、請求項１記載の方法。
【請求項３】
複数のフィルタ係数セットを所定のシフトそれぞれに対して計算する、請求項２記載の方法。
【請求項４】
それぞれ１つのフィルタ係数セットが計算される前記所定のシフトは、０から前記ソース画像内の隣接する２つのピクセルの間隔までの全てのシフト、並びに一定の間隔の中間ステップを含む、請求項３記載の方法。
【請求項５】
半分のナイキスト周波数の領域において周波数特性が上昇している伝達関数を使用する、請求項２から４までのいずれか１項記載の方法。
【請求項６】
前記上昇がガウス型の経過を有し、かつ１５〜３０％である、請求項５記載の方法。
【請求項７】
原本を光学的に走査する対物レンズを有する光学走査装置において使用し、前記伝達関数を前記対物レンズの設定に依存して変化させる、請求項１から６までのいずれか１項記載の方法。
【請求項８】
−フィルタ係数を計算する手段と、
−請求項１から７までのいずれか１項記載の方法によりソース画像を目標画像にマッピングする手段と、
を有することを特徴とする、請求項１から７までのいずれか１項記載の方法を実施する装置。
【請求項９】
ソース画像を目標画像にマッピングする前記手段はハードウェア回路として構成されており、
−列カウンタ（１）を有し、該列カウンタ（１）の後段には列数をスケーリング係数と乗算する乗算器（３）が接続されており、
−行カウンタ（２）を有し、該行カウンタ（２）の後段にはスケーリング係数と乗算する乗算器（４）が接続されており、
−列数および行数に関する補正手段（５，６，７，８，９，１０）を有し、
−前記列数に対応付けられているＦＩＲフィルタテーブル（１３）を有し、
−前記行数に対応付けられているＦＩＲフィルタテーブル（１４）を有し、
−ローカル列カウンタ（１５）を有し、該ローカル列カウンタ（１５）の値に基づき、対応するフィルタ係数が前記ＦＩＲフィルタテーブル（１３）から読み出され、且つ該ローカル列カウンタ（１５）の値は加算器（１７）を用いて前記列数に加算され、
−ローカル行カウンタ（１６）を有し、該ローカル行カウンタ（１６）の値に基づき、フィルタ係数が前記ＦＩＲフィルタテーブル（１４）から読み出され、且つ該ローカル行カウンタ（１６）の値は加算器（１８）を用いて前記行数に加算され、
−アドレス発生器を有し、該アドレス発生器は前記ローカル列カウンタないし前記ローカル行カウンタにより補完されており、且つ補間されている列数／行数に基づき、画像メモリ（２０）に関するアドレスを形成し、該画像メモリ（２０）には前記ソース画像が記憶されており、
−２つの乗算器（２１，２２）を有し、該乗算器（２１，２２）は前記アドレス発生器によって形成されたアドレスによって設定されたピクセルの値を前記フィルタ係数と乗算し、
−アキュムレータ（２３）を有し、該アキュムレータ（２３）は前記フィルタ係数と乗算される全ての画像値を、前記ローカル行カウンタの各値に対する前記ローカル列カウンタの実行にわたり累算する、
請求項８記載の装置。
【請求項１０】
マイクロプロセッサを有し、該マイクロプロセッサにおいて請求項１から７までのいずれか１項記載の方法を実施するためのデータ処理プログラムが記憶されており、且つ実行される、請求項８記載の装置。
【請求項１１】
マイクロプロセッサを有し、前記フィルタ係数を計算するデータ処理プログラムが記憶されており、且つ実行される、請求項８または９記載の装置。
【請求項１２】
請求項８から１１までのいずれか１項記載の装置を有することを特徴とするスキャナ。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像を補間および補正する方法および装置に関する。
【０００２】
コンピュータ読み取り可能なデータファイルとして存在するディジタル画像においては、画像をオリジナルの解像度とは異なる解像度で提供することが頻繁に必要とされる。このことは殊にスキャナを用いてスキャンされた画像に該当し、この種の画像はスキャナの光学系およびセンサによって検出された解像度とは異なる解像度で出力されるべきことが多い。
【０００３】
さらには実際の光学系（対物レンズ）によるセンサへの光学的なマッピングによって、光学系、センサ、後続のアナログ的な増幅また非理想的なスキャナによって惹起される多数の欠陥が生じる。要求の高いスキャナにおいては、公知の方法を用いてこの偏差が補正される。画像を補正するためにＦＩＲフィルタ（ＦＩＲ：有限長インパルス応答）、ＩＩＲフィルタ（ＩＩＲ：無限長インパルス応答）および非線形フィルタ、例えばプルーニング（ｐｒｕｎｉｎｇ）またはメディアン（ｍｅｄｉａｎ）フィルタを使用することが公知である。非線形のフィルタ方法は、例えば原本における傷のような局所的なノイズを除去することに比べて、光学的／電子的な記録システムを補正することには適していない。
【０００４】
画像の解像度の変更は補間方法を用いて実施される。最も頻繁に使用される補間方法はＢスプラインを用いる双一次補間、四次補間および二次元補間である。
【０００５】
要求の高いスキャナにおいては、先ず検出された画像が補間により所望の解像度に変換され、続いて画像の補正が実施される。この際、例えば先ず画像の歪みを修正し、続いて解像度を変更するために、補間を複数のステップで実施することができる。
【０００６】
以下では図４〜７に基づき、従来の補間方法が画像に及ぼす作用を周波数領域で表し説明する。周波数領域においては、画像が線に沿って振動によって表される。図４は、２つのピクセル間の距離の１／４を表すｄｘ＝０．２５ピクセルのシフトの際の双一次補間の影響を示す。横座標にはナイキスト周波数（０．５）までの正規化された空間周波数がプロットされている。ナイキスト周波数は、線に沿って連続する黒と白のピクセルを表す空間周波数である。空間周波数が大きくなればなるほど、振幅の減少はますます大きくなることが見て取れる。このことは双一次補間の場合にｄｘ＝０．２５のシフトにおいては線が密接している状態を表す鮮明なコントラストが弱まるということを意味している。ｄｘ＝０．２５のシフトは例えば画像を３００ＤＰＩから２３３ＤＰＩに変換する際に頻繁に適用される。
【０００７】
図５においてはｄｘ＝０．２５のシフトに関して正規化された周波数にわたる位相経過が示されている。画像シフトが１／４ピクセルだけ実施されるべき最適な補間関数に関しては、周波数にわたる位相シフトが、原点を通過し且つナイキスト周波数において４５°の値を有する直線を生じさせなければならない。この目標関数は図５において破線で示されている。位相誤差は半分のナイキスト限界において４．５°であり、すなわちピクセル長に関して２．５％であり、それ以降はナイキスト周波数まで大幅に上昇する。
【０００８】
図６および図７はｄｘ＝０．５ピクセルについての相応のグラフを示す。図６からはｄｘ＝０．５ピクセルの補間においては振幅がさらに大幅に低減することが見て取れる。これに対して位相誤差はこの補間においては生じない（図７）。
【０００９】
四次補間ないしスプライン補間においては双一次補間よりも誤差は少ない。しかしながら全ての従来の補間方法に関して、シフトの絶対値に依存する振幅誤差および位相誤差が生じる。任意の補間に関しては画像にわたり必要とされるシフトが変化するので、これらの誤差は補間された画像においては局所的に様々に表されている。
【００１０】
US 5,301,266からは、画像が先ず周波数領域においてフィルタリングされ、その後画像が位置空間に移され、そこにおいて補間が実施される方法が公知である。
【００１１】
EP 1 041 510 A2からは、高周波成分の範囲においてノイズを有していると思われる信号を処理する方法およびシステムが公知である。この種のノイズは殊に線形の補間において惹起される。殊に複数の画像を処理する際にその種のノイズを補正するために小さい２次元ＦＩＲフィルタ、殊にインバースフィルタが使用され、これらのフィルタは補間によって惹起されるローパス特性を補正する。個々の補間ステップに対して異なるフィルタが使用され、これらのフィルタは事前に計算され、またフィルタバンクに記憶される。その種の各フィルタは、先行する補間によって惹起されている周波数特性におけるノイズを補正する。
【００１２】
ディジタル画像の解像度の変更および補正を行う公知の方法は、これらのエラーを高品質で補償することはできない。何故ならば、このために必要とされる情報を１つの補間ステップの終了後ではもはや使用することができないからである。補間誤差は領域ごとに少なくないので、従来の方法は補正不可能な重大な欠点を有する。
【００１３】
COURMONTAGEN, P.; CAVASSILAS, J.Fによる「An interpolation-filtering method for noise-corrupted Images」、Digital Signal Processing、１９９７年、第２巻、第１０９１〜１０９４頁からは、画像の補間とフィルタリングを同時に行う方法が公知である。この方法においては画像がフーリエ変換を用いて周波数領域に変換される。周波数領域においてはフィルタリングも補間もそれぞれ乗法演算によって実施されるので、両方の演算、すなわちフィルタリングと補間を同時に実施することができる。しかしながらこの演算は周波数領域において行われているので、これらの演算をその都度画像全体にわたり実施する必要がある。２５６×２５６ピクセルの画像大きさを有する画像では従来のコンピュータを用いて演算を容易に実施することができる。より大きい画像フォーマットでは、性能の良いパーソナルコンピュータを用いるだけでは、演算を短時間で実施することは不可能である。いずれにせよ周波数領域におけるこの方法を、リアルタイムで画像の補間とフィルタリングを同時に実施するために使用することはできない。この方法が基礎としているフーリエ変換の一般的な原理は、IWAINSKY A., WILHELMI W. Lexikon der Computergrafik und Bildverarbeitung, Vieweg Verlagsgesellschaft、１９９４年、第１０８頁に記載されている。
【００１４】
本発明が基礎とする課題は、上述の欠点の影響を受けず、それと同時に簡単且つ迅速に実施することができ、その結果簡潔で廉価なハードウェア回路においても実現することができる、ディジタル画像の解像度を変更と補正を実施する方法を提供することである。
【００１５】
この課題は、請求項１の特徴を有する方法および請求項８の特徴を有する装置によって解決される。本発明の有利な実施形態は従属請求項に記載されている。
【００１６】
ディジタル画像を補間および補正する本発明による方法によれば、ソース画像が複数のフィルタ係数を有するＦＩＲフィルタを用いて目標画像にマッピングされる。フィルタ係数は画像を補間するための情報も、画像を補正するための情報も包含するので、単一のステップで画像の補間も補正も実施される。
【００１７】
本発明によれば、１つのステップで解像度を変更するための補間も補正も実施する単一のフィルタが使用されるので、解像度の変更に関しても周波数補正に関しても全ての本来の情報が存在し、また計算されるので、従来の方法から公知であるような連続的な誤差は発生しない。このことは補間による誤差にも丸め誤差にも該当する。
【００１８】
ＦＩＲフィルタを分割することができるので、このＦＩＲフィルタを使用することにより位置空間における計算を画像の水平線または垂直線に沿った方向において行うことができる。このことはソフトウェア処理手段を用いる本発明による方法の非常に迅速な実施も、簡単且つ廉価なハードウェアによる実現も可能にする。
【００１９】
本発明による方法を簡単にハードウェア回路として実現することができる（請求項９）。
【００２０】
以下では本発明を図面に基づき詳細に説明する。図面において
図１は、数学的な畳み込みに応じたフィルタ計算を用いるソース画像のピクセルの目標画像へのマッピングを概略的に示し、
図２は、スキャンされたディジタル画像を補正するための例示的な伝達関数を示し、
図３は、各行に本発明によるＦＩＲフィルタに関するフィルタ係数のセットをそれぞれ１つずつ包含するテーブルを示し、フィルタはそれぞれ異なるシフトないし補間のために設けられている。
図４は、ｄｘ＝０．２５のシフトに関する従来の双一次補間における振幅の経過を示し、
図５は、ｄｘ＝０．２５のシフトに関する従来の双一次補間における位相の経過を示し、
図６は、ｄｘ＝０．５のシフトに関する従来の双一次補間における振幅の経過を示し、
図７は、ｄｘ＝０．５のシフトに関する従来の双一次補間における位相の経過を示し、
図８ａ、８ｂから１６ａ、１６ｂはそれぞれｄｘ＝０ないし０．１２５ないし０．２５ないし０．３７５ないし０．５ないし０．６２５ないし０．７５ないし０．８７５ないし１のシフトに関する振幅ないし各位相の経過を示し、
図１７は本発明による方法を実施するための装置の概略的なブロック図を示し、
図１８ａ、１８ｂは周波数特性を補正する別の伝達関数ならびにその要素を示す。
【００２１】
本発明は以下の状況を基礎とする：
画像処理においてはフィルタが因果的であってはならない。これにより振幅と位相を相互に依存せずに調節することができる。
【００２２】
本発明によれば線形で離散的なＦＩＲフィルタが使用される（「Digitale Bildverarbeitung」第４版、Bernd Jaehne - ISBN 3-540-61379-x, Kapitel 4 Nachbarschafts-operatoren）。目標画像の個々のピクセルは行毎にソース画像のピクセルから計算され、この際数学的な畳み込みに応じたフィルタ計算には、目標画像の計算すべきピクセルに相当するソース画像の所定数のピクセルが考慮される。このことは以下の式により表される：
【数１】
2007538436.tif 000003
ここで、Ｇ′mnはｎ番目の列およびｍ番目の行におけるマッピング点G′であり、Ｇmnはｎ番目の列およびｍ番目の行におけるソース画像のピクセルであり、ａiはフィルタ係数である。このことは図１に概略的に示されており、ここではソース画像の点Ｇmnが３つの係数ａiを有するフィルタを用いて目標画像における点Ｇ′mnにマッピングされる。
【００２３】
本発明の本質は、フィルタ係数ａiにより計算ステップにおいて画像の補間も周波数特性の補正も行われるようフィルタ係数ａiが選定されることである。したがって、補間ならびに周波数補正を実施するためにはただ１つのフィルタしか使用されない。
【００２４】
以下ではフィルタ係数ａiの計算を詳細に説明する。このために理想的には連続的な伝達関数
2007538436.tif 000004
を基礎とする。伝達関数
2007538436.tif 000005
はソース画像を目標画像にマッピングする際の周波数特性を補正するために使用される。この伝達関数は近似的に離散的なフーリエ変換H（ω）の伝達関数に置換される：
【数２】
2007538436.tif 000006
【００２５】
位置空間においては、ソース画像が連続関数ｆ（ｘ）によって１つの行におけるいて表される。サブピクセルシフト、すなわち隣接する２つのピクセルの間隔よりも短い距離ｄｘだけシフトする場合には、位置空間における関数はｆ（ｘ＋ｄｘ）で表される。フーリエ変換を用いる位置空間から周波数領域への移行を以下のように表すことができる：
【数３】
2007538436.tif 000007
ここでＦ（ω）は周波数領域に変換されたｆ（ｘ）の関数である。このことは、周波数領域におけるサブピクセルシフトを位相ｉωｄｘについての位相シフトとして表すことができ、且つこの位相シフトないしサブピクセルシフトがシフトされていない関数Ｆ（ω）と位相係数ｅiωdxの積として表されるということを意味している。この位相シフトないしサブピクセルシフトは補間をなす。何故ならば、このシフトによって２つのピクセル間の値が検出されるからである。したがってこれにより伝達関数に関しては、補間を表すためにこの伝達関数を位相要素と乗算するだけでよいということになる：
【数４】
2007538436.tif 000008
【００２６】
この式により本発明の実施が実質的に容易になる。何故ならば、この式は周波数領域において伝達関数Ｈ（ω）および位相シフトないし補間ｅiωdxが、相互に依存せずに検出することができる積の２つの独立したファクタであることを示すからである。
【００２７】
ＦＩＲフィルタの伝達関数ＨFIRを周波数領域においてフーリエ列を用いて以下のように表すことができる：
【数５】
2007538436.tif 000009
【００２８】
式（５）によるフーリエ行列の係数ａnは式（１）についてのフィルタ係数ａiに対応する。このことは、それ自体は自由に選択可能であり、且つ式（４）により位相要素と乗算される所定の伝達関数に関するフーリエ行列の係数を計算することによって、所定の伝達関数および所定の補間に関する式（１）のフィルタ係数を計算できるということを意味している。したがってフィルタ係数のセットを計算することができ、これらのセットを適用することによりその都度所定の周波数補正および補間が行われる。
【００２９】
したがって、（ＨFIR−Ｈdx（ω））の偏差が最小になるようにａnを決定することが目標とされる。この一致はｋに関して伝達関数Ｈ（ω）のｋ個の離散的な空間周波数に関して要求される：
【数６】
2007538436.tif 000010
【００３０】
このことはｎ＝７の係数およびｋ＝９の離散的な空間周波数を有するマトリクス式において表される：
【数７】
2007538436.tif 000011
【００３１】
解ベクトルａnを発見するために、この系を複素マトリクス要素ｅ-iωxによって解く必要がある。
【００３２】
このことは線形の最適化の別個の問題を生じさせ、この問題は「２次制約付き２次計画（Quadratically Constrained Quadratic program）」として公知である。この場合係数はチェビシェフフィルタのような特性を有する伝達関数をもたらす。
【００３３】
解決ストラテジ集はGoldfarb D.およびIyengar G.によるCORC Report TR-2002-04「Robust convex quadratically constrained programs」に記載されている。このレポートはインターネットにおいて下記のアドレスから入手することができる：
http://www.gsb.stanford.edu/facseminars/events/oit/pdfs/lyengar_Robust.pdf
以下では本発明によるＦＩＲ補間フィールドの実施例をｎ＝７の係数およびｋ＝１５の離散的な空間周波数について説明する。
【００３４】
周波数補正を例えば、周波数特性の上昇が半分のナイキスト周波数の領域において得られるように構成することができる。上昇はガウス型の経過を表す。上昇は例えば２５％である。理想的な補間フィルタはナイキスト周波数において増幅０を有する。これへの移行はカイザー窓を介して達成される。伝達関数はこれらの関数の組み合わせにより表される。
【００３５】
この種の周波数補正は、多くの光学系が比較的高い周波数、すなわち非常に密に隣接している点ないし線のコントラスト（良質なコントラスト）を比較的低い周波数、すなわち比較的大きい点ないし比較的太い線に比べて通常の場合粗悪にマッピングするという事実を考慮する。したがって増幅は半分のナイキスト周波数まで若干高められる。
【００３６】
図２にはこの伝達関数の経過（絶対値関数）が示されており、この伝達関数を用いて周波数補正が実施される。この際、横軸には１５の離散的な空間周波数がωにより正規化されてプロットされており、また縦軸には増幅率Ａがプロットされている。離散的な空間周波数は、最大値がナイキスト周波数の１５に相当するように選定されており、したがって最大値は座標系において７．５、すなわち半分のナイキスト周波数を表している。
【００３７】
別の光学的な特性を有する他の光学系においては、周波数特性を補正する伝達関数を異なるように構成することが勿論好適である。
【００３８】
本明細書において説明する伝達関数は周波数補正のみを内容とするので、補間を実施するためにこの伝達関数をさらに位相要素ｅiωdx（式（４）を参照されたい）と乗算する必要がある。このために１５の離散的な空間周波数に対して、また０×ｄ，０．１×ｄ，０．２×ｄ，．．．，０．９×ｄの１０のシフトそれぞれに対して位相要素が計算されており、ここでｄはソース画像における隣接する２つのピクセルの距離である。したがって、０．１ｄのパターンでの考えられる各シフトに対して、１５の離散的な空間周波数に関するセット位相要素が計算されており、また伝達関数と乗算されている。
【００３９】
つまり所定のシフトｄｘに対して与えられている値は式（７）の右側におけるベクトルＨdx（ω）を表す。
【００４０】
したがって１０のシフトｄｘそれぞれに対して、式（７）による方程式のセットが作成された。
【００４１】
線形の最適化のための上記の方法を用いて、１０の解ベクトルａnを求めるための１０の方程式のセットが計算された。各補間ステップに対してこれらの複数の方程式のセットを用いることにより、位置空間において作用するＦＩＲフィルタのフィルタ係数に補間を精密に組み込むことができる。
【００４２】
この際に得られるフィルタ係数ａnの値は図３の表に記載されており、ここでは各行に１つのＦＩＲフィルタに対してそれぞれ７個のフィルタ係数が所定のシフトｄｘに関して示されている。式（１）および図１によれば値ｄｘ＝０だけシフトする場合、すなわち一般的な語法ではシフトが行われない場合には、補間により、その位置がソース画像内のピクセルに対応する目標画像内のピクセルが計算され、このピクセルは係数ａ4と乗算される。図３の表から見て取れるように、係数ａ4はｄｘ＝０のシフトにおいて１．００１６９であるので、ソース画像からのこの値は目標画像の対応する値に非常に良好に移行する。しかしながらｄｘ＝０．９のシフトの際には、係数ａ5と乗算されるピクセルへの方向における、係数ａ4と乗算されるピクセルから隣接する２つのピクセルまでの距離の０．９倍のシフトに基づき、目標画像において計算すべきピクセルの位置がシフトしている。図３のテーブルに基づき、ｄｘ＝０．９ｄのシフトの場合には、係数ａ4が単に０．２３３０４２であるが、係数ａ5がこれに対して０．９９４０２６であることが見て取れる。すなわち、ソース画像内での位置が目標画像内での相応の位置のより近くに位置し、且つ係数ａ5と乗算されるピクセルが、係数ａ4と乗算される遠く離れたピクセルよりも実質的に大きく計算において考慮される。したがって係数ａ4はシフトが大きくなるに連れ段階的に低下し、それと共に係数ａ5は相応にシフトが大きくなるに連れ増加する。これにより計算されたフィルタの補間効果を具体的に跡づけすることができる。
【００４３】
周波数補正のための情報も補間のための情報も包含するフィルタ係数を前述の方法により簡単に計算することができるので、所定の用途において伝達関数を変更し、相応に異なるフィルタ係数を計算することが可能である。このことは例えば、焦点距離が異なる場合に異なる伝達関数が対応付けられているズーム対物レンズが使用される場合に有効である。したがって焦点距離を調節する度に、フィルタ係数が相応に適合された伝達関数を用いて計算される。
【００４４】
上述したように、周波数特性を補正する他の伝達関数も使用することができる。以下では非常に有利な伝達関数を図１８ａおよび１８ｂを用いて説明する。
【００４５】
スキャンの際にはスキャンすべき画像がスキャナの対物レンズの焦点内に正確に位置してないのでピントのボケが生じ、このボケは以下のボケ伝達関数ｈ１（ω）（図１８ａ）によって表される：
【数８】
2007538436.tif 000012
【００４６】
ピクセルは有限の大きさを有するので、高周波数においてはボケが生じ、このボケは以下のアパーチャ伝達関数ｈ２（ω）によって表される：
【数９】
2007538436.tif 000013
【００４７】
全体のボケはｈ１（ω）とｈ２（ω）の積ｐ（ω）によって表され、これは図１８ａに示されている。
【００４８】
周波数特性を補正する伝達関数は第１のアプローチにおいては積ｐ（ω）の逆数として計算される。しかしながら殊にωが高い場合には、この関数が著しく上昇する可能性があるので顕著なノイズが生じる虞がある。
【００４９】
最適なフィルタは上述した伝達関数ｈ１（ω）およびｈ２（ω）を以下の式に対応するウィーナーフィルタの計算に考慮することができる：
【数１０】
2007538436.tif 000014
【００５０】
周波数特性を補正するこの伝達関数は図１８ｂに示されており、また図２に示されている伝達関数と同様に本発明によるＦＩＲフィルタのフィルタ係数を生じさせるために使用することができる。図１８ａ、１８ｂにおいては伝達関数が系の最大の波数について正規化されており、この際０．５がナイキスト周波数に対応する。
【００５１】
図１７においては本発明による方法を実施するための装置がブロック図で概略的に示されている。
【００５２】
この装置は列カウンタ１および行カウンタ２を有する。列カウンタ１は目標画像における列をカウントし、行カウンタ２は目標画像における行をカウントする。列カウンタ１が目標画像における全ての列を一通りカウントし終えると、行カウンタが行についてのカウントをさらに実施する。したがって目標画像のピクセルは行毎に順次カウントされていく。したがってこれら２つのカウンタ１，２は、目標画像内のどのピクセルに対して、ソース画像からのマッピングが実施されるべきであるかを設定する。
【００５３】
列カウンタ１および行カウンタ２の後段にはそれぞれ乗算器３，４が接続されており、これらの乗算器３，４を用いて列数ないし行数が対応するスケーリング係数と乗算される。スケーリング係数はソース画像におけるステップ幅を表し、その結果目標画像内の２つのピクセルの間隔、すなわち隣接する２つの列ないし行の間隔は、ソース画像内の隣接する２つのピクセルの間隔の単位でソース画像にマッピングされる。
【００５４】
乗算器３，４の後段にはそれぞれ加算器５，６が接続されており、これらの加算器５，６はそれぞれの列数ないし行数にそれぞれ１つの開始座標を加算する。この開始座標を用いて一方では、目標画像にマッピングされるべきソース画像における画像の一部が確定される。他方では開始座標がＦＩＲフィルタの中心をソース画像内のそれぞれのピクセルに置く。以下ではこれを詳細に説明する。
【００５５】
このようにして補正された列数および行数はそれぞれ補正データ発生器７，８および加算器９，１０に供給される。補正データ発生器７，８は所定の関数またはルックアップテーブルに基づき、供給された列数および行数から補正値を計算し、この補正値は例えば所定の光学系を使用してソース画像を形成する際に惹起される歪みエラーを補正する。この補正値はそれぞれの加算器９，１０に供給され、相応の列数ないし行数と加算される。このようにして得られた列数および行数はソース画像内の座標を表し、この座標からは、列カウンタ１および行カウンタ２を用いて設定された列数および行数を用いて目標画像のピクセルを計算するための上述のＦＩＲフィルタを計算することができる。
【００５６】
加算器９の後段には調整補正を行うエレメント１１が接続されており、また分岐路には減算器１２が接続されている。調整補正を行うエレメント１１および減算器１２は本発明による装置を説明する際には差し当たり考慮されない。ＦＩＲフィルタテーブル１３に案内されるこの分岐路を介して、ＦＩＲフィルタテーブル１３には列数の分数部分が供給される。
【００５７】
加算器１０の後段には相応にＦＩＲフィルタテーブル１４が接続されており、このＦＩＲフィルタテーブル１４は行数の分数部分が供給される。列数ないし行数の整数部分は加算器１７ないし１８に供給される。加算器１７，１８とＦＩＲフィルタテーブル１３，１４の間にはそれぞれローカル列カウンタ１５ないしローカル行カウンタ１６が接続されている。所定の列数ないし行数が加算器１７，１８に加えられると、ローカル列カウンタ１５およびローカル行カウンタ１６はそれぞれ０〜ｎ−１までカウントし、この際ｎはフィルタ係数の数である。ローカル列カウンタ１５のカウント数ｉがＦＩＲフィルタテーブル１３に転送され、（ｉ＋１番目のフィルタ値）にしたがい２番目の分数部分に対応するフィルタから相応のフィルタ値が選択される。相応にして、ＦＩＲフィルタテーブル１４においてはフィルタ値が分数部分およびローカル行カウンタ１６のカウント値に依存して選択される。ローカル列カウンタ１５のカウンタ値は加算器１７において整数の列数に加えられ、これによって対応するフィルタ値に応じた列数が計算される。相応にして、ローカル行カウンタ１６のカウンタ値が、加算器１８に供給された整数の行数に加えられる。加算器１７の出力値（列数）および加算器１８の出力値（行数）は座標組を形成し、この座標組はアドレス発生器１９に供給される。アドレス発生器１９においては、これらの座標が対応するアドレスデータに変換され、このアドレスデータは画像メモリ２０内のアドレスを表す。画像メモリ２０はソース画像のデータを包含し、またアドレス発生器１９によって形成されたアドレスは画像メモリにおいて、アドレス発生器１９に供給された座標組に対応するソース画像のデータを包含する。画像メモリからは相応の値が読み出され、先ず第１の乗算器２１においてＦＩＲフィルタテーブル１４からのフィルタ値と乗算され、次いで第２の乗算器２２においてＦＩＲフィルタテーブル１３からのフィルタ値と乗算される。ピクセルを計算するために、ローカル列カウンタ１５はローカル行カウンタ１６における各ステップにつきカウントを実施する。すなわちローカル列カウンタ１５が実施される度に、ローカル行カウンタ１６のカウンタ値は１増える。この際に求められた値はアキュムレータ２３において加算され、目標画像におけるピクセルの値として出力され、この値は列カウンタ１によって設定された列数および行カウンタ２によって設定された行数によって規定されている。
【００５８】
相応の画像値の出力後には列カウンタのカウンタ値が１増え、新たな画像値が計算される。この際ローカル列カウンタ１５およびローカル行カウンタ１６はその値範囲を一度カウントし終える。
【００５９】
列カウンタが最後の列に到達すると、列カウンタは再び０にセットされ、第１の列において開始され、また行カウンタのカウンタ値が１増える。列カウンタも行カウンタもその最大値に達すると、完全な目標画像が形成される。
【００６０】
ＦＩＲフィルタは分離可能であるので、本発明によるフィルタ計算を列方向においても行方向においても同時に実施し、第１の乗算器２１および第２の乗算器２２における連続的な乗算によってオーバラップさせることができる。
【００６１】
上述の調整補正を行うエレメントを用いることにより、ピクセルを表すために８ビットしか必要とされないにもかかわらず、通常の場合３２ビットを包含するデータ語により処理する計算機の特性を十分に利用することができる。したがってデータ語の読み出しによって、１回のアクセスで４つのピクセルの値を取り出すことができる。調整補正を行うエレメントはこのことを考慮し、アドレス発生器１９において４で整除できる次に小さいアドレスが生じるように列数をセットする。この丸め込みは左シフトの作用を有する。
【００６２】
このシフトは相応の乗算装置によって再び元に戻され、この乗算装置はこのシフトを再び補償調整する。しかしながらより上質には、この整数のシフトがフィルタにおけるサブピクセルシフトと共に処理され、またいずれにせよ存在するＦＩＲフィルタテーブルを用いて補償される。このために差分器１２がアドレス発生器１９に供給された列数と加算器９から出力された列数との差を形成する。差は整数シフトの和を表し、また分数部分によって形成されるサブピクセルシフトの和を表す。したがってフィルタを用いてサブピクセルシフトのみでなく、メモリアクセスアドレスの丸め込みに対応する複数のピクセルについてのシフトも実施される。
【００６３】
図１７に基づき説明した上記の装置はソース画像を目標画像にマッピングするためのハードウェア回路である。このハードウェア回路は簡単で標準的なコンポーネントから構成されており、したがってこのハードウェア回路を廉価に製造することができる。さらにこのハードウェア回路は非常に迅速にソース画像を目標画像にマッピングすることができる。
【００６４】
もちろん本発明を、相応のマイクロプロセッサにおいて実施されるデータ処理プログラムによっても実現することができる。
【００６５】
本発明による方法はスキャナ、殊にＤＩＮのＡ４サイズの紙を毎分１００ページまで走査し、補間し、また相応の伝達関数を補正することができる高性能スキャナにおいて使用される。
【００６６】
本発明による原理を、例えばディジタルカメラのような他の適用領域にも使用することができる。
【００６７】
本発明を以下のように要約することができる：
本発明は、ディジタル画像を補間および補正する方法および装置に関する。ソース画像は複数のフィルタ計数を有するＦＩＲフィルタを用いて目標画像にマッピングされる。フィルタ係数は画像を補間するための情報も、画像を補正するための情報も包含するので、単一の処理ステップで補間も補正も実施される。
【００６８】
ソース画像を目標画像にマッピングする本発明による方法をハードウェア回路により簡単に実現することができるので、装置、殊にスキャナにおけるマッピングは全体のプロセス経過に遅延が生じることなく実施される。
【図面の簡単な説明】
【００６９】
【図１】数学的な畳み込みに応じたフィルタ計算を用いるソース画像のピクセルの目標画像へのマッピング。
【図２】スキャンされたディジタル画像を補正するための例示的な伝達関数。
【図３】各行に本発明によるＦＩＲフィルタに関するフィルタ係数のセットをそれぞれ１つずつ包含するテーブル。
【図４】ｄｘ＝０．２５のシフトに関する従来の双一次補間における振幅の経過。
【図５】ｄｘ＝０．２５のシフトに関する従来の双一次補間における位相の経過。
【図６】ｄｘ＝０．５のシフトに関する従来の双一次補間における振幅の経過。
【図７】ｄｘ＝０．５のシフトに関する従来の双一次補間における位相の経過
【図８ａ】ｄｘ＝０のシフトに関する振幅の経過。
【図８ｂ】ｄｘ＝０のシフトに関する位相の経過。
【図９ａ】ｄｘ＝０．１２５のシフトに関する振幅の経過。
【図９ｂ】ｄｘ＝０．１２５のシフトに関する位相の経過。
【図１０ａ】ｄｘ＝０．２５のシフトに関する振幅の経過。
【図１０ｂ】ｄｘ＝０．２５のシフトに関する位相の経過。
【図１１ａ】ｄｘ＝０．３７５のシフトに関する振幅の経過。
【図１１ｂ】ｄｘ＝０．３７５のシフトに関する位相の経過。
【図１２ａ】ｄｘ＝０．５のシフトに関する振幅の経過。
【図１２ｂ】ｄｘ＝０．５のシフトに関する位相の経過。
【図１３ａ】ｄｘ＝０．６２５のシフトに関する振幅の経過。
【図１３ｂ】ｄｘ＝０．６２５のシフトに関する位相の経過。
【図１４ａ】ｄｘ＝０．７５のシフトに関する振幅の経過。
【図１４ｂ】ｄｘ＝０．７５のシフトに関する位相の経過。
【図１５ａ】ｄｘ＝０．８７５のシフトに関する振幅の経過。
【図１５ｂ】ｄｘ＝０．８７５のシフトに関する位相の経過。
【図１６ａ】ｄｘ＝１のシフトに関する振幅の経過。
【図１６ｂ】ｄｘ＝１のシフトに関する位相の経過。
【図１７】本発明による方法を実施するための装置の概略的なブロック図。
【図１８】周波数特性を補正する別の伝達関数ならびにその要素。
【符号の説明】
【００７０】
１  列カウンタ、  ２  行カウンタ、  ３、４  乗算器、  ５，６，９，１０，１７，１８  加算器、  ７，８  補正データ発生器、  １１  調整補正を行うエレメント、  １２  減算器、  １３，１４  ＦＩＲフィルタテーブル、  １５  ローカル列カウンタ、  １６  ローカル行カウンタ、  １９  アドレス発生器、  ２０  画像メモリ、  ２１  第１の乗算器、  ２２  第２の乗算器、  ２３  アキュムレータ、  ２４  出力側
【図１】
2007538436.tif 000016
【図２】
2007538436.tif 000017
【図３】
2007538436.tif 000018
【図４】
2007538436.tif 000019
【図５】
2007538436.tif 000020
【図６】
2007538436.tif 000021
【図７】
2007538436.tif 000022
【図８ａ】
2007538436.tif 000023
【図８ｂ】
2007538436.tif 000024
【図９ａ】
2007538436.tif 000025
【図９ｂ】
2007538436.tif 000026
【図１０ａ】
2007538436.tif 000027
【図１０ｂ】
2007538436.tif 000028
【図１１ａ】
2007538436.tif 000029
【図１１ｂ】
2007538436.tif 000030
【図１２ａ】
2007538436.tif 000031
【図１２ｂ】
2007538436.tif 000032
【図１３ａ】
2007538436.tif 000033
【図１３ｂ】
2007538436.tif 000034
【図１４ａ】
2007538436.tif 000035
【図１４ｂ】
2007538436.tif 000036
【図１５ａ】
2007538436.tif 000037
【図１５ｂ】
2007538436.tif 000038
【図１６ａ】
2007538436.tif 000039
【図１６ｂ】
2007538436.tif 000040
【図１７】
2007538436.tif 000041
【図１８ａ】
2007538436.tif 000042
【図１８ｂ】
2007538436.tif 000043
【国際調査報告】
2007538436.tif 000044
2007538436.tif 000045
2007538436.tif 000046
2007538436.tif 000047
2007538436.tif 000048
2007538436.tif 000049