JP4663349B2 - スクリーン周波数を推定する方法及び装置 - Google Patents

Description

本出願は、本出願と同じ日に提出され、その全体を参照してここに組み込む米国出願（代理人整理番号１１７５４４、１１７５２１、１１７７４５、１１７７４６、１１７７４７、１１７７４８、１１８５８４、１１８５９１及び１１８６０１）に関連する。

本出願は、すべて２００２年７月０１日に提出され、すべて共通して本件出願人に譲受され、その内容を参照してここに組み込む「文書のディジタル式網消し処理」なる米国出願番号１０／１８７，４９９号（代理人整理番号Ｄ／Ａ１２７０）、「文書のディジタル式網消し処理のための制御システム」なる米国出願番号１０／１８８，０２６号（代理人整理番号Ｄ／Ａ１２７０Ｑ）、「文書のマルチプルラスターコンテンツ（ＭＲＣ）表示のための動的閾値システム」なる米国出願番号１０／１８８，２７７号（代理人整理番号Ｄ／Ａ１２７１Ｑ）、「文書のマルチプルラスターコンテンツ（ＭＲＣ）表示のための分離システム」なる米国出願番号１０／１８８，１５７号（代理人整理番号Ｄ／Ａ１２７１Ｑ１）及び「マルチプルラスターコンテンツ（ＭＲＣ）ＴＩＦＦ及びＰＤＦのためのセグメント化技法」なる米国出願番号６０／３９３，２４４号（代理人整理番号Ｄ／Ａ２３０３Ｐ）を含む同時係属出願に関連する。

本発明は一般的には画像処理の方法及びシステムに関し、より詳しくはディジタル走査された文書を網消しする方法とシステムに関する。

ハロゲン化銀式写真撮影を例外としてほとんど全ての印刷物は、ハーフトーンスクリーンを用いて印刷されている。ハーフトーンの周波数と大きさを推定する必要性は、染料昇華式の又はハロゲン化銀式写真撮影のような少数のデバイスを例外として、ほとんど全ての印刷物がハーフトーンスクリーンを用いて印刷されている事実に基づく。このようなハーフトーンは印刷デバイスによって非常に特定的なものであり、走査されて再ハーフトーン化される場合、適切に除去されないと、可視アーティファクト及び／又は受け入れがたいモアレ（干渉縞）パターンを発生させる可能性がある。ハーフトーンの抑制はカラー文書の場合特に重要であるが、それは、このようなハーフトーンは一般的に、角度や周波数が異なるとそれに対応して少し異なったスクリーンを含む４つ以上の色分解で印刷され、また、互いに干渉して好ましくない空間アーティファクトを引き起こすからである。

オリジナルのハーフトーンスクリーンを成功裏に除去できるかどうかは、ローカル周波数を正確に推定する能力しだいである。したがって、ハーフトーンスクリーンの周波数と大きさを推定する改良型の方法と装置に対する必要性が存在する。

テキストと線画との品質とシャープさとを維持しながらも、走査された文書からハーフトーンスクリーンを消去する高度に効率的な方法とシステムを開示する。

単一の画像データ処理チャネルのみを用いて下流の画像処理、特に、画像からのハーフトーン信号の網消し処理のためのスクリーン周波数推定値を生成するスクリーン推定モジュールを開示する。

ただ１つのチャネルしか用いないとはいえ、スクリーン推定モジュールは、ハーフトーンスクリーン周波数の高品質で信頼できる推定値を生成することが可能である。コントラストが低い又はゼロに近い場合、周波数推定値は高度にフィルタリングされた画像信号に基づいていて正確な周波数測定値ではないため、単一チャネルスクリーン周波数推定値で十分である。しかしながら、コントラストが低い場合、ごくわずかなハーフトーンノイズが入力画像において測定され、その正確な周波数に関する知識は、下流の処理では必要とされない。広範囲にわたるエッジシャープ化効果を、好ましくないハーフトーンアーティファクトを悪化させることなく、この状況下での画像データに施すことができる。コントラストが高い場合、周波数推定値は画像中に発生するハーフトーン周波数の信頼できる測定値に基づく。したがって、単一チャネルスクリーン周波数推定値を、あらゆる範囲の画像コントラストにわたって用いることが可能である。

ただ１つのチャネルしか用いないため、本発明を実施することによって、コスト、電力及びデバイスのパッケージサイズをかなり節約することが可能である。本発明が特に効果的である領域は、解像度が３００〜６００ｄｐｉの範囲である。

本発明の特徴及び利点は、以下の本発明に関する詳細な説明及び添付の図面から明らかであろう。

ディジタル式に走査された文書を網消し、これにより、潜在的なハーフトーン干渉および好ましくないモアレパターンを消去又は顕著に減少させるための新しい方法およびシステムを説明する。ここで図１を参照すると、本発明の方法とシステムのブロック図が、単一チャネルスクリーン推定器モジュールＳＥＭ４０で表されている。このスクリーン推定器モジュール４０は、目的とする現行の画素における瞬間的なハーフトーン周波数を推定する役目を負っている。スクリーン推定器モジュール４０は、８ビットのソース画像Ｓｒｃ２８に対して動作して８ビットのハーフトーン周波数推定値Ｓｃｍ７２を生成する。

ハーフトーン周波数と大きさを推定する必要性は、染料昇華式又はハロゲン化銀式写真撮影などの少数のデバイスを例外としてほとんど全ての印刷物がハーフトーンスクリーンを用いて印刷されるという事実に基づいている。このようなハーフトーンは印刷デバイスに非常に固有のものであり、印刷目的で走査されたり再ハーフトーン化される場合、これを適切に除去しないと、可視アーティファクト及び／又は受け入れがたいモアレパターンを発生させる可能性がある。本件出願人の同時係属出願の一つ（代理人整理番号１１８５９１）に記載された網消しモジュール（ＤＳＣ）は、オリジナルの走査済み画像からオリジナルのハーフトーンパターンを消去（フィルタリング除去）すべく、スクリーン推定器モジュールによって生成された情報に依存している。ハーフトーンの抑制はカラー文書の場合特に重要であるが、それは、このようなハーフトーンは一般的に、角度や周波数が異なるとそれに対応して少し異なったスクリーンを含む４つ以上の色分解で印刷され、また、互いに干渉して好ましくない空間アーティファクトを引き起こすからである。

先行技術によるスクリーン推定器モジュールは、異なる感度レベルで最大で３つの周波数チャネルを用いていた。上位チャネルは、最大ソース解像度で周波数感度が最大となるようにチューニングされており、したがって、スクリーン周波数推定信号を導くために用いられる。しかしながら、このチャネルは非常に敏感であり、スクリーンが非常に弱い場合であっても、周波数の存在が報告されるのが常であった。したがって、スクリーン周波数は、スクリーンの大きさＳｃｍ７２によってさらに限定されていた。

図１の単一チャネル４０は、中程度の周波数感度となるようにチューニングされており、最大ソース解像度で動作する。スクリーン周波数信号Ｓｃｍ７２は、単一チャネルによって生成された周波数推定値を分析することによって導かれる。

単一チャネル４０は、以下に説明する最小／最大テクスチャ検出器ＭＭ３ ３２と、後述の平均化フィルタ４２とによって構成されている。単一チャネルＭＭ３ ３２ユニットは、単一チャネルの８ビットの入力ソース信号Ｓｒｃ２８に基づいて動作する。ＭＭ３最小／最大モジュール３２は、２Ｄ入力信号中の山と谷を発見するために用いられる。最小／最大検出器ユニットを以下に詳細に説明する。このユニットは、基本的に、目的とする現行の画素に中心がある３×３のウインドウのコンテンツを検査し、中心画素がそれを取り囲む８個の隣接画素より著しく大きかったり小さかったりするかを適応閾値を用いて分析する。そうであれば、中心画素は、それぞれ、山（隣接画素より大きい場合）又は谷（隣接画素より小さい場合）と見なされる。単位面積あたりの山と谷の数をカウントすることによって、ローカル周波数の測定値（尺度）が得られる。

ＭＭ３ユニットの出力３２の精度はわずか１ビットであるが、第１の後続のフィルタリングステージの前に構成係数ＤｏｔＧａｉｎでスケーリングされる。このユニットは、入力信号の１つ又は複数のカラーチャネルで動作する。しかしながら、本実施形態では、ただ１つのチャネル、すなわち、輝度チャネルが用いられる。フィルタ出力のスケーリングは、第１の後続のフィルタ正規化ステップまで、このステージの正規化係数を調整することによって延期することが可能である。

ＭＭ３最小／最大検出器３２の出力は、さまざまな平均化フィルタとサブサンプリングフィルタとを通過する。このサブサンプリングに伴うエイリアシング問題を回避するため、どの場合でも空間フィルタの範囲は、サブサンプリング率の２倍から１を減算した値である。

同様に、単一チャネルＭＭ３ ３２の出力は、２つの三角形２ＤサブサンプリングフィルタすなわちＦ３１／１６Ｆフィルタ４２とＦ３／２フィルタ４６とからなるカスケードに印加される。これらのカスケーディングされたフィルタリングユニットの出力は、各々の方向において係数３２×で（第１のフィルタでは１６×で第２のフィルタでは２×で）サブサンプリングされる。

単一チャネル４０では、解像度１／１６の信号のサンプルがＭＸ３ユニット４４を通過する。このユニットは最大３×３の演算（グレー膨張）を実行する。その出力は、それぞれ双線形補間ユニットＢＩＵ５４の入力ｂに送られる。

この単一チャネルは、空間ノイズをさらに減少させるさらなる平滑化／平均化Ｆ５ユニットステージ６４を含んでいる。Ｆ５ユニットステージ６４は５×５の三角形重み（非サブサンプリング）フィルタである。このユニットからのフィルタリング済み出力は、双線形補間ユニットＢＩＵ５４の入力に送られる。その出力もまた、Ｃ３コントラストユニット４８を通過するが、このユニットは、現行の画素に中心がある３×３ウインドウ中の最大の差を探索する。出力Ｃ３はＢＩＵユニット５４への入力ｃとなる。

単一チャネル４０によって生成された３つの信号はＢＩＵユニット５４に送られる。このユニットは、双線形補間を実行して、サブサンプリングされた入力解像度をオリジナルのソース解像度に戻す。ＢＩＵの入力ａとｃは解像度１／３２であり、入力ｂは解像度１／１６である。補間ユニットからの出力の帯域幅は、入力より実質的に高い。たとえば、上記の係数３２×の場合、補間ユニットは各入力画素に対して１０２４の出力画素を生成する。

双線形補間ユニット５４の出力は、スクリーン大きさ信号Ｓｃｍ７２の８ビットの推定値である。この推定されたスクリーン周波数信号Ｓｃｍ７２は、網消しモジュールＤＳＣ５８とセグメント化モジュールＳＥＧ５６とにエクスポートされる。信号チャネルスクリーン推定器モジュール４０のさまざまな構成要素は、以下により詳細に説明する。

図２にさまざまなフィルタユニットの一次元フィルタ応答を示し、図３〜５にさまざまなユニットの二次元フィルタ応答を示す。これらのフィルタユニットは、入力信号を平滑化すなわち平均化することによって高周波数を除去する目的で用いられる。これらフィルタユニットは各々が、方形で、分離式で、対称性のある２Ｄ ＦＩＲ（有限インパルス応答）フィルタを実行する。そのフィルタ応答は水平方向と垂直方法とで同じである。フィルタに対する入力がカラー信号である場合、同じフィルタ応答が色成分の１つずつに対して独立に適用される。１Ｄフィルタ６０の応答は、図２に示すように、整数の係数を持つ対照的な三角形の形状をしている。一回の実施毎に特定のフィルタ形状が選ばれている（しかし、他の任意のフィルタ形状も含まれる）。

一般的なフィルタ形状はＦｎ／ｋフィルタと呼ばれるが、ここでｎはフィルタサイズ（ｘ方向又はｙ方向の全体的な範囲）であり、ｋはフィルタリングされた出力に対して各方向で適用されるサブサンプリングの分量である。このサブサンプリング係数ｋは、ｋ＝１であれば省略される。本明細書において、フィルタ範囲ｎは、２Ｄフィルタ応答が有効な中心画素位置で有限のピーク値を有するように奇数（ｎ＝１、３、５、．．．）であると仮定されることに注意されたい。

１Ｄフィルタ応答と２Ｄフィルタ応答の例を図１と図２に示す。図２はＦ３とＦ１１の非正規化１Ｄフィルタ６０応答を示し、図３〜５はそれぞれＦ３ ６２、Ｆ５ ６４及びＦ７ ６６に関して結果的に得られる非正規化２Ｄ係数を示す。

フィルタは分離可能であるため、２Ｄフィルタ応答は、２つの１Ｄフィルタを水平方向と垂直方向にカスケーディングすることによって実施することが可能である。これらのフィルタは全て最高の入力データ速度で動作しているが、出力は各方向で係数ｋでサブサンプリングされる。いつもではないが、多くの場合、フィルタサイズｎとサブサンプリング係数ｋは次の関係を満足する：

これは、サブサンプリングされた面積に対する被覆オーバラップ率が５０％であることを表している。例として、最小の３×３フィルタであるＦ３ ６２の全体的な２Ｄ応答は次の通りである：

これより大きいフィルタも同様に説明される。これらのフィルタは分離可能であるので、これらは互いに直交する２つの１Ｄステップで実施するのが最良である。各フィルタ出力は、係数の和で正規化され、これにより、８ビット範囲に再び適合するようにされる。Ｆ３フィルタ６２などの一部のフィルタの場合は、その重みの合計は２の累乗の値となる。これらのフィルタは、丸めて２つだけ右方向に桁移動させる演算として簡単に実施することが可能なので、正規化ステップで除算を必要としない。たとえば、Ｆ３フィルタ６２の１Ｄ重みの合計値は１＋２＋１＝４である。この重みによる丸め除算は、２を加算して、その後で２つだけ右方向に桁移動させて実施することが可能である。

正規化結果＝（和＋２）＞＞２

一般に、丸めることが必要とされる場合、それは一般的に、除数を半分にして加算してから桁移動させることによって実施される。２の補数であるコーディング済みの二進数に対して実行される右方向桁移動は、フロアー（ｆｌｏｏｒ）（分子／２^桁移動）と等価であるため、除数の半分を加算すると、符号付の分子の場合と符号なしの分子の場合の双方で、最も近い整数に丸められる。

フィルタの重みの合計が２の累乗とならない場合、演算集約的な除算は、２つの数の比で乗算して近似させることによって回避されるが、ここで分母は選ばれた２の累乗の数である。

サブサンプリングフィルタＦ３／２Ｆ１５／８Ｆ３１／１６及びＦ６３／３２はすべて、２の累乗の１Ｄ重み、すなわち、４，６４，２５６と１０２４をそれぞれ有している。したがって、正規化は単に丸め右方向桁移動である。Ｆ５フィルタ６４は９という１Ｄ重みを有し、この場合、５７で乗算して近似させ、その後で、９だけ丸め右方向桁移動することが可能である。ここで、数ｘを５７で乗算する場合、次のような桁移動／加算／減算を用いることによって、可変乗算なしで実行可能である：

図６Ａと６Ｂを参照すると、単一チャネルで用いられるＭＸ３最大値ユニット４４は、目的とする現行の画素７４に中心を持つ３×３ウインドウ中の最大値を探索する。入力は８ビットの信号である。この最大値探索は、３×３ウインドウの９画素にわたって実行される。このグレー膨張モジュールは、この探索ウインドウの境界内で発見される最大の画素値７６から成る８ビット出力を生成する。ＭＸ３最大値アルゴリズムを図６Ｂに示す。

ここで図７Ａと７Ｂを参照すると、これらのＣ３コントラストモジュール４８は、入力でのローカルコントラストの量を測定するように設計されている。コントラストは、目的とする現行の画素７４に中心を持つウインドウ内での最大画素値と最小画素値の差として定義される。Ｃ３コントラストユニット４８は、目的とする現行の画素７４に中心を持つ３×３のウインドウのサイズを利用する。このコントラストユニットへの入力は８ビットの信号である。コントラストモジュール４８は、８ビットの白黒出力（単一チャネル）８４を生成する。Ｃ３コントラストユニット４８の動作を図７Ｂに示す。この動作は次の通りである：各画素位置に対して、３×３ウインドウのコンテンツを独立に探索して、最小画素値と最大画素値とを求める。その出力のコントラスト値は次のように定義される：

コントラスト＝最大値−最小値

符号なし８ビット入力信号の場合、最大画素値と最小画素値は常に０と２５５の間に存在するため、コントラストは［０．．．．２５５］という範囲内にあることが保証され、そのため、特殊な正規化は不必要である。

最小／最大検出モジュール３２は、入力信号中の山と谷を発見するために用いられる。単位面積あたりの山と谷の数をカウントすることによって、ローカル周波数の尺度が得られる。

ＭＭ３ユニット３２は、１成分のグレーソースに基づいて動作する。このユニットは、３×３ウインドウを利用して、現行画素がその８個の隣接する画素に対していつ極値（山又は谷）に達しているかを、以下に示すロジックにしたがって示す。最小／最大検出ユニット３２の出力は、対応するＳｒｃ画素が極値状態にある（他のカラーチャネルにも拡大可能である）ことを示す１ビット信号である。

ＭＭ３最小／最大検出構造を図８に示す。各画素に対して、それ（目的とする現行画素）を取り囲む８つの画素の外側のリングが最初に分析される。この８個の外側の画素リングをさらに、図８に示すように、各々が４画素から成る２つの集合に分割する。外側の画素リングをこのように２つの集合に区分化すると、直線のセグメントをハーフトーンとして検出する際に間違って警告を発する可能性を減少させるのに役に立つ（最も一般的に遭遇するハーフトーンはクラスタ化されたドットである可能性が高いからである）。

次のように、各集合に対して、画素値を、この集合のメンバー７８と８６の間で比較して、各集合内での最小値と最大値を独立に決定する：
Ａ_max＝集合Ａに属す全ての（ｉ，ｊ）の中の最大値（Ａｉｊ）
Ａ_min＝集合Ａに属す全ての（ｉ，ｊ）の中の最小値（Ａｉｊ）
Ｂ_max＝集合Ｂに属す全ての（ｉ，ｊ）の中の最大値（Ｂｉｊ）
Ｂ_min＝集合Ｂに属す全ての（ｉ，ｊ）の中の最小値（Ｂｉｊ）
これらから、外側リング全体と最小値の合計とが計算される。最小値合計と２つの構成パラメータとを用いて、ノイズレベルが次のように計算される：

中心画素７４の値Ｘは、次に示すように、それがどちらかの集合の最大画素値よりも「かなり」大きい場合、山にあると定義される：

同様に、中心画素７４の値Ｘは、次に示すように、それがどちらかの集合の最小画素値よりも「かなり」小さい場合、谷にあると定義される：

上の式は、３×３検出ウインドウからの出力が１に設定される２つの状況を決定するが、他の場合では全て、出力は０に設定される。

単一チャネルスクリーン推定器モジュールＳＥＭは、双線形補間ユニットＢＩＵを利用する。双線形補間ユニットは、信号を補間（アップサンプリング）してソース解像度に戻す。入力信号は各方向で３２という係数でアップサンプリングされ、そのオリジナルの解像度を回復する。各補間ユニットは、双線形補間を実行して、基本的に、各オリジナル画素に対して３２＊３２＝１０２４個の画素を生成する。この双線形補間のステップサイズは、オリジナルの画素グリッドの１／３２である。双線形補間ユニットを以下により詳細に説明する。

双線形補間ユニットＢＩＵ５４の１つの構造を図１０に示す。双線形補間ユニットは、単一チャネル４０中で生成される３つの信号９４、９６及び９８に基づいて動作する。

図１０から分かるように、双線形補間ユニット５４の各々は、２つの補間ステージ１００と１０２から成っている。第１のステージは、各方向で入力Ａ₅９４とＣ₅９８を２ｘ倍する補間ステージ１００を含んでいる。この補間ステージ１００は、簡単な双線形補間技法を用いている。入力Ａ₅９４は、Ｆ５フィルタ６４ユニットに対応している。ここで、図１０の添え字はサブサンプリングのレベルに対応している。添え字５は、その信号が１／２という係数で５回（合計で１／３２）サブサンプリングされていることを示している。Ｃ₅入力９８は、３×３コントラストユニットの出力に対応している。図１に示すように、これらの入力は双方ともが、各方向で３２×という係数で前もってサブサンプリングされている。この第１の補間ステージの出力Ａ₄とＣ₄は、補間された後、１／１６でサブサンプリングされる。これは、Ｂ₄入力９６と同じサブサンプリングレベルである。これで、Ｂ₄マイナスＡ₄差分信号１０４であるＢｍＡ₄を計算することが可能となる。ＢｍＡ₄は、Ｃ₄をＭａｇＦｉｎｅＢｌｅｎＶｓＣｏｎ関数１０６に適用することによって生成された大きさ微細混合係数ＭＦＢ₄で乗算される。ＢｍＡｘＣ₄信号は、１０８ＢｍＡ₄にＭＦＢ₄を乗算して８だけ桁移動させることによって得られた結果である。これをステージ１１０中のＡ₄に加算して、チャネルに応じてＨＩ₄又はＬＯ₄信号を作成する。次に、その結果を１６×双線形補間ユニット１１２に供給して、チャネルに応じてＬｏ出力やＨｉ出力を生成する。

上記のＭａｇＦｉｎｅＢｌｅｎＶｓＣｏｎ関数１０６は、プログラム可能関数である。１実施形態において、一般的なＭａｇＦｉｎｅＢｌｅｎＶｓＣｏｎ関数１０６は、ｙ＝（ｘ−１６）＊１２として容易に計算可能であり、この出力は次に０と１９２との間で固定される。下に示す式は、この一般的なＭａｇＦｉｎｅＢｌｅｎＶｓＣｏｎ１０６の構成値を組み込んでいる。

ＭＦＢ₄の場合に定義される上記の関数関係は、単一チャネル４０の動作を成功裏に実行させることを保証する助けとなる。この関係は、１６から始まり、傾斜１２をたどり、１９２でゼロに交差する直線を定義している。制御信号がａとｂの混合を調整し、この結果、ＭＸ３ ４４から誘導され、従って３×３画素ウインドウ内で発生する周波数の測定値であるＢの割合が、コントラストが増加するにつれ増加することが、この関係によって保証される。

コントラストが非常に低い場合、双線形補間モジュールは、高度にフィルタリングされた画像信号である入力Ａに基づいて信号を出力する。コントラストが大きい場合、この出力信号は、周波数測定値である成分Ｂに向けてより重み付けされる。したがって、出力信号Ｓｃｍはさらなる大きさ推定値なしで単独で用いることが可能であるが、それは、出力信号Ｓｃｍが小さい場合、比較的少量のハーフトーンノイズが存在し、また、その周波数の正確な測定は下流の処理では必要とされないからである。しかしながら、出力信号Ｓｃｍが大きい場合、該信号は主として周波数測定値に基づいており、したがって、画像中に存在するハーフトーン周波数の信頼性のある推定値である。

以上、特定の例示的な実施形態を詳細に説明し、添付図面に示したが、当業者であれば、本発明は上記の実施形態に限られないこと、また、さまざまな修正が、上記の本発明の図示の実施形態及び他の実施形態に対して、その広範な発明の範囲から逸脱することなく可能であることが認識されよう。したがって、本発明は開示されている特定の実施形態や装置に制限されず、むしろ、添付のクレームによって定義される本発明の範囲と精神を超えない任意の変更、適応又は修正をカバーすることが理解されよう。

単一チャネルスクリーン推定モジュールのシステムのブロック図である。 さまざまなフィルタユニットの一次元フィルタ応答を示す図である。 さまざまなフィルタユニットの二次元フィルタ応答を示す図である。 さまざまなフィルタユニットの二次元フィルタ応答を示す図である。 さまざまなフィルタユニットの二次元フィルタ応答を示す図である。 一般的な３×３の最大モジュール構造を示す図である。 一般的な３×３の最大モジュール構造を示す図である。 一般的な３×３のコントラストモジュール構造を示す図である。 一般的な３×３のコントラストモジュール構造を示す図である。 ３×３のウインドウ内の最小／最大検出構造を示す図である。 単一補間ユニットを示す図である。 単一補間ユニットを示す図である。 双線形補間ユニットの構造のブロック図である。 スクリーン周波数推定方程式を示す図である。

Claims (8)

1. ハーフトーンスクリーンを用いて形成された画像を走査して得られた画像データから前記ハーフトーンスクリーンのパターンを除去するためのスクリーン周波数を推定する方法であって、
   前記画像データに対応する第１の画像データ信号の解像度が該解像度よりも小さい第１の解像度になるように該第１の画像データ信号をサブサンプリングすることにより第２の画像データ信号を生成する第１生成ステップと、
   前記第１生成ステップで生成された第２の画像データ信号の前記第１の解像度が該第１の解像度よりも小さい第２の解像度になるように該第２の画像データ信号をサブサンプリングすることにより第３の画像データ信号を生成する第２生成ステップと、
   前記第１生成ステップで生成された第２の画像データ信号に対して、対応する画素毎に、注目画素を中心に配置した最大値探索用のウィンドウを用いて最大値を探索し、探索して得られた最大値を示す最大値信号を出力する最大値出力ステップと、
   前記第２生成ステップで生成された第３の画像データ信号に対して、対応する画素毎に、前記注目画素を中心に配置したコントラスト算出用のウィンドウを用いてコントラストを算出し、算出したコントラストを示すコントラスト信号を出力するコントラスト出力ステップと、
   前記コントラスト出力ステップによって出力されたコントラスト信号の解像度が前記第１の解像度になるように該コントラスト信号を補間することにより第１補間済み信号を生成する第３生成ステップと、
   前記第２生成ステップで生成された第３の画像データ信号の解像度が前記第１の解像度となるように該第３の画像データ信号を補間することにより第２補間済み信号を生成する第４生成ステップと、
   前記最大値出力ステップによって出力された最大値信号から、前記第４生成ステップで生成された第２補間済み信号を、対応する画素毎に減算する減算ステップと、
   前記第３生成ステップで生成された第１補間済み信号を対応する画素毎に下記（１）式に代入して得られた制御信号（ＭＦＢ_４）を、対応する画素毎に、前記減算ステップで減算されて得られた信号に乗算する乗算ステップと、
   前記乗算ステップで乗算されて得られた信号を、対応する画素毎に、前記第４生成ステップで生成された第２補間済み信号に加算する加算ステップと、
   前記加算ステップによって加算されて得られた加算済み信号の解像度が前記第１の画像データ信号の解像度になるように該加算済み信号を補間することにより前記スクリーン周波数の推定値を示す推定値信号を生成する第５生成ステップと、
   を含む方法。
   ＭＦＢ_４＝ｍａｘ（０，ｍｉｎ（１９２，（Ｃ_４−１６）＊１２））・・・・（１）
   Ｃ_４；第１補間済み信号
2. 前記画像データにおいて、前記注目画素を中心に配置した前記ウィンドウ内の前記注目画素に隣接する８個の画素を４画素からなる２つの集合に分割し、分割して得られた２つの集合の各々の最大の画素値の何れかよりも前記注目画素の画素値が所定閾値以上に大きい場合に及び前記２つの集合の各々の最小の画素値の何れかよりも前記注目画素の画素値が所定閾値以上に小さい場合にのみ前記注目画素が極値状態にあることを示す極値信号を用いて前記第１の画像データ信号を生成し、生成した前記第１の画像データ信号を前記第１生成ステップに提供する提供ステップを更に含む請求項１記載の方法。
3. 前記第２生成ステップで生成された第３の画像データ信号を、空間ノイズを減少させるための５×５の非サブサンプリングフィルタによりフィルタリングするフィルタリングステップを更に含み、
   前記第４生成ステップは、前記フィルタリングステップによってフィルタリングして得られた信号を前記第２の係数によりアップサンプリングすることにより前記第２補間済み信号を生成する請求項１または請求項２記載の方法。
4. 前記第５生成ステップで生成された推定値信号を網消しモジュールに出力する出力ステップを更に含む請求項１〜請求項３の何れか１項に記載の方法。
5. ハーフトーンスクリーンを用いて形成された画像を走査して得られた画像データから前記ハーフトーンスクリーンのパターンを除去するためのスクリーン周波数を推定する装置であって、
   前記画像データに対応する第１の画像データ信号の解像度が該解像度よりも小さい第１の解像度になるように該第１の画像データ信号をサブサンプリングすることにより第２の画像データ信号を生成する第１生成手段と、
   前記第１生成手段で生成された第２の画像データ信号の前記第１の解像度が該第１の解像度よりも小さい第２の解像度になるように該第２の画像データ信号をサブサンプリングすることにより第３の画像データ信号を生成する第２生成手段と、
   前記第１生成手段で生成された第２の画像データ信号に対して、対応する画素毎に、注目画素を中心に配置した最大値探索用のウィンドウを用いて最大値を探索し、探索して得られた最大値を示す最大値信号を出力する最大値出力手段と、
   前記第２生成手段で生成された第３の画像データ信号に対して、対応する画素毎に、前記注目画素を中心に配置したコントラスト算出用のウィンドウを用いてコントラストを算出し、算出したコントラストを示すコントラスト信号を出力するコントラスト出力手段と、
   前記コントラスト出力手段によって出力されたコントラスト信号の解像度が前記第１の解像度になるように該コントラスト信号を補間することにより第１補間済み信号を生成する第３生成手段と、
   前記第２生成手段で生成された第３の画像データ信号の解像度が前記第１の解像度となるように該第３の画像データ信号を補間することにより第２補間済み信号を生成する第４生成手段と、
   前記最大値出力手段によって出力された最大値信号から、前記第４生成手段で生成された第２補間済み信号を、対応する画素毎に減算する減算手段と、
   前記第３生成手段で生成された第１補間済み信号を対応する画素毎に（１）式に代入して得られた制御信号（ＭＦＢ_４）を、対応する画素毎に、前記減算手段で減算されて得られた信号に乗算する乗算手段と、
   前記乗算手段で乗算されて得られた信号を、対応する画素毎に、前記第４生成手段で生成された第２補間済み信号に加算する加算手段と、
   前記加算手段によって加算されて得られた加算済み信号の解像度が前記第１の画像データ信号の解像度になるように該加算済み信号を補間することにより前記スクリーン周波数の推定値を示す推定値信号を生成する第５生成手段と、
   を含む装置。
6. 前記画像データにおいて、前記注目画素を中心に配置した前記ウィンドウ内の前記注目画素に隣接する８個の画素を４画素からなる２つの集合に分割し、分割して得られた２つの集合の各々の最大の画素値の何れかよりも前記注目画素の画素値が所定閾値以上に大きい場合に及び前記２つの集合の各々の最小の画素値の何れかよりも前記注目画素の画素値が所定閾値以上に小さい場合にのみ前記注目画素が極値状態にあることを示す極値信号を用いて前記第１の画像データ信号を生成し、生成した前記第１の画像データ信号を前記第１生成手段に提供する提供手段を更に含む請求項５記載の装置。
7. 前記第２生成手段で生成された第３の画像データ信号を、空間ノイズを減少させるための５×５の非サブサンプリングフィルタによりフィルタリングするフィルタリング手段を更に含み、
   前記第４生成手段は、前記フィルタリング手段によってフィルタリングして得られた信号を前記第２の係数によりアップサンプリングすることにより前記第２補間済み信号を生成する請求項５または請求項６記載の装置。
8. 前記第５生成手段で生成された推定値信号を網消しモジュールに出力する出力手段を更に含む請求項５〜請求項７の何れか１項に記載の装置。

Images