JPH1098628A - ハーフトーン画像からの連続色調画像の作成方法 - Google Patents
ハーフトーン画像からの連続色調画像の作成方法Info
- Publication number
- JPH1098628A JPH1098628A JP9241963A JP24196397A JPH1098628A JP H1098628 A JPH1098628 A JP H1098628A JP 9241963 A JP9241963 A JP 9241963A JP 24196397 A JP24196397 A JP 24196397A JP H1098628 A JPH1098628 A JP H1098628A
- Authority
- JP
- Japan
- Prior art keywords
- image
- halftone
- sub
- subband
- band
- 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.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N1/00—Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
- H04N1/40—Picture signal circuits
- H04N1/40075—Descreening, i.e. converting a halftone signal into a corresponding continuous-tone signal; Rescreening, i.e. combined descreening and halftoning
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
- Color Image Communication Systems (AREA)
- Dot-Matrix Printers And Others (AREA)
- Color, Gradation (AREA)
- Control Or Security For Electrophotography (AREA)
Abstract
(57)【要約】
【課題】 ハーフトーン画像を非スクリーニングするた
めの非スクリーニング方法を提供する。 【解決手段】 ハーフトーン画像から連続色調画像を作
るための方法であって、ハーフトーン画像を複数のサブ
バンド画像に分解する工程（Ｓ１２０）と、この複数の
サブバンド画像から雑音とパターンの内の少なくとも一
つを選択的に除去する工程（Ｓ１３０）と、複数のサブ
バンド画像を連続色調画像に再構成する工程（Ｓ１４
０）とを含む。
めの非スクリーニング方法を提供する。 【解決手段】 ハーフトーン画像から連続色調画像を作
るための方法であって、ハーフトーン画像を複数のサブ
バンド画像に分解する工程（Ｓ１２０）と、この複数の
サブバンド画像から雑音とパターンの内の少なくとも一
つを選択的に除去する工程（Ｓ１３０）と、複数のサブ
バンド画像を連続色調画像に再構成する工程（Ｓ１４
０）とを含む。
Description
【０００１】
【発明の属する技術分野】本発明は、元の連続色調画像
から発生させたハーフトーン画像を連続色調の画像に逆
変換するための方法と装置に関する。特に、ウエーブレ
ット(wavelet) に基づく技術を使用し、ハーフトーン技
術およびハーフトーンの知識とは無関係にハーフトーン
画像から高品質の連続色調画像を発生できる方法と装置
に関する。
から発生させたハーフトーン画像を連続色調の画像に逆
変換するための方法と装置に関する。特に、ウエーブレ
ット(wavelet) に基づく技術を使用し、ハーフトーン技
術およびハーフトーンの知識とは無関係にハーフトーン
画像から高品質の連続色調画像を発生できる方法と装置
に関する。
【０００２】
【従来の技術及び発明が解決しようとする課題】ハーフ
トーン化は、連続色調、すなわちコントーン(contone)
またはグレースケール（無彩色スケール）の画像から２
値化した画像、すなわち２色調の画像を発生することを
指す。一般に、グレースケールの画像は、黒（または他
の単一色）と白とによる連続色調画像を指すが、コント
ーン画像は、フルカラーの画像または単一色と白とによ
る画像であることができる。何れの場合も、ハーフトー
ン画像は、種々のハーフトーン化技術の任意のものを使
用してコントーン（フルカラーまたはグレースケール）
の画像から発生される。これらの種々の技術には次のも
のが含まれる。１）集合ドット、分散ドット、確率スク
リーンなどの閾値アレーまたはディサリング（ditherin
g ）；２）エラー拡散などの適応プロセス；３）最小自
乗法、直接バイナリーサーチなどの対話型プロセス。な
お、このリストは完全なものではないことを理解された
い。
トーン化は、連続色調、すなわちコントーン(contone)
またはグレースケール（無彩色スケール）の画像から２
値化した画像、すなわち２色調の画像を発生することを
指す。一般に、グレースケールの画像は、黒（または他
の単一色）と白とによる連続色調画像を指すが、コント
ーン画像は、フルカラーの画像または単一色と白とによ
る画像であることができる。何れの場合も、ハーフトー
ン画像は、種々のハーフトーン化技術の任意のものを使
用してコントーン（フルカラーまたはグレースケール）
の画像から発生される。これらの種々の技術には次のも
のが含まれる。１）集合ドット、分散ドット、確率スク
リーンなどの閾値アレーまたはディサリング（ditherin
g ）；２）エラー拡散などの適応プロセス；３）最小自
乗法、直接バイナリーサーチなどの対話型プロセス。な
お、このリストは完全なものではないことを理解された
い。
【０００３】また、ハーフトーン化は多くのものを一つ
にする（many-to-one ）プロセスである。すなわち、特
定のハーフトーン化技術が使用されたか否かに関わら
ず、元のコントーン画像の多くの変形（variations）が
同一のハーフトーン画像になるであろう。
にする（many-to-one ）プロセスである。すなわち、特
定のハーフトーン化技術が使用されたか否かに関わら
ず、元のコントーン画像の多くの変形（variations）が
同一のハーフトーン画像になるであろう。
【０００４】非スクリーニング(unscreening) は、元の
コントーン画像のハーフトーン化を「取り消す」ための
プロセスを指す。より正確には、非スクリーニングは、
元のハーフトーン画像からコントーン画像を作り出すこ
とを指す。ハーフトーン化は多くのものを一つにするプ
ロセスであるため、非スクリーニングは、一般に一つの
ものを多くのものにする（one-to-many ）プロセスであ
る。すなわち、得られたコントーン画像の多くの変形の
任意の一つを、元のハーフトーン画像から非スクリーニ
ングによって作り出すことができる。
コントーン画像のハーフトーン化を「取り消す」ための
プロセスを指す。より正確には、非スクリーニングは、
元のハーフトーン画像からコントーン画像を作り出すこ
とを指す。ハーフトーン化は多くのものを一つにするプ
ロセスであるため、非スクリーニングは、一般に一つの
ものを多くのものにする（one-to-many ）プロセスであ
る。すなわち、得られたコントーン画像の多くの変形の
任意の一つを、元のハーフトーン画像から非スクリーニ
ングによって作り出すことができる。
【０００５】非スクリーニングの目的は、元のハーフト
ーン自体を作るのに使用された元のコントーン画像に最
も近似した非スクリーニングされたコントーン画像を発
生することである。したがって、従来の非スクリーニン
グの技術は、元のハーフトーン画像を作るのに使用した
特定のハーフトーン技術に基づいて開発されている。す
なわち、従来の各非スクリーニング技術は、特定のハー
フトーン技術を使用して作られたハーフトーン画像を非
スクリーニングするためのみに設計された。
ーン自体を作るのに使用された元のコントーン画像に最
も近似した非スクリーニングされたコントーン画像を発
生することである。したがって、従来の非スクリーニン
グの技術は、元のハーフトーン画像を作るのに使用した
特定のハーフトーン技術に基づいて開発されている。す
なわち、従来の各非スクリーニング技術は、特定のハー
フトーン技術を使用して作られたハーフトーン画像を非
スクリーニングするためのみに設計された。
【０００６】例えば、Ｐ．Ｗ．ウォン（Wong）による
「誤差拡散のための逆ハーフトーン化および核推定（In
verse Halftoning and Kernel Estimation for Error D
iffusion）」、画画像処理に関するＩＥＥＥ議事録、Ｖ
ｏｌ．４、Ｎｏ．４、４８６〜４９８頁、１９９５年４
月では、誤差拡散されたハーフトーン画像を非スクリー
ニングするための繰り返し法が三つ開示されている。こ
の三つの繰り返し法とは、１）線型フィルタリングと統
計的平滑化、２）最小平均二乗誤差（ＭＭＳＥ）投影、
および３）最大アプリオリ確率（ＭＡＰ）投影である。
これらの繰り返し法の各々は誤差拡散したハーフトーン
画像中に存在する特徴に基づくものである。非スクリー
ニングされるべきハーフトーン画像が、誤差拡散によっ
て元のコントーン画像から形成されていない場合、これ
らの方法は適切でなくなる。なぜなら、それらの方法は
誤差拡散されていないハーフトーン画像から質の劣った
コントーン画像を発生するためである。
「誤差拡散のための逆ハーフトーン化および核推定（In
verse Halftoning and Kernel Estimation for Error D
iffusion）」、画画像処理に関するＩＥＥＥ議事録、Ｖ
ｏｌ．４、Ｎｏ．４、４８６〜４９８頁、１９９５年４
月では、誤差拡散されたハーフトーン画像を非スクリー
ニングするための繰り返し法が三つ開示されている。こ
の三つの繰り返し法とは、１）線型フィルタリングと統
計的平滑化、２）最小平均二乗誤差（ＭＭＳＥ）投影、
および３）最大アプリオリ確率（ＭＡＰ）投影である。
これらの繰り返し法の各々は誤差拡散したハーフトーン
画像中に存在する特徴に基づくものである。非スクリー
ニングされるべきハーフトーン画像が、誤差拡散によっ
て元のコントーン画像から形成されていない場合、これ
らの方法は適切でなくなる。なぜなら、それらの方法は
誤差拡散されていないハーフトーン画像から質の劣った
コントーン画像を発生するためである。
【０００７】例えば、４９４頁と４９６頁のアルゴリズ
ム２および３においてウォンが開示したＭＡＰ投影非ス
クリーニング技術では、誤差拡散したハーフトーン画像
を作るのに使用した誤差拡散核(kernel)が事前に分かっ
ていなくても、誤差拡散されたハーフトーン画像を非ス
クリーニングすることができる。しかし、ＭＡＰ投影非
スクリーニング技術では、ハーフトーン画像は誤差拡散
にによって形成されたものであると仮定される。したが
って、ＭＡＰ投影非スクリーニング技術を実行するシス
テムでは、誤差拡散されたハーフトーン画像のみが供給
されるか、誤差拡散されていないハーフトーン画像を排
除するためのゲートキーパー（門番）が必要となる。さ
もないと、そのようなシステムに、例えば集合ドット技
術によって形成されたハーフトーン画像が供給された場
合、ＭＡＰ投影非スクリーニングシステムが評価し、画
像を非スクリーニングするための、ハーフトーン画像に
埋め込まれた誤差拡散核が存在しないであろう。したが
って、そのようなシステムが誤差拡散されていないハー
フトーン画像に対して動作した場合には、アルゴリズム
２のステップ５の要件を満たさない可能性が高く、した
がって、終わりのないループに入ったままになるであろ
う。
ム２および３においてウォンが開示したＭＡＰ投影非ス
クリーニング技術では、誤差拡散したハーフトーン画像
を作るのに使用した誤差拡散核(kernel)が事前に分かっ
ていなくても、誤差拡散されたハーフトーン画像を非ス
クリーニングすることができる。しかし、ＭＡＰ投影非
スクリーニング技術では、ハーフトーン画像は誤差拡散
にによって形成されたものであると仮定される。したが
って、ＭＡＰ投影非スクリーニング技術を実行するシス
テムでは、誤差拡散されたハーフトーン画像のみが供給
されるか、誤差拡散されていないハーフトーン画像を排
除するためのゲートキーパー（門番）が必要となる。さ
もないと、そのようなシステムに、例えば集合ドット技
術によって形成されたハーフトーン画像が供給された場
合、ＭＡＰ投影非スクリーニングシステムが評価し、画
像を非スクリーニングするための、ハーフトーン画像に
埋め込まれた誤差拡散核が存在しないであろう。したが
って、そのようなシステムが誤差拡散されていないハー
フトーン画像に対して動作した場合には、アルゴリズム
２のステップ５の要件を満たさない可能性が高く、した
がって、終わりのないループに入ったままになるであろ
う。
【０００８】ディサリング（dithering ）されたハーフ
トーン画像を非スクリーニングすることに関する他のシ
ステムがＣ．Ｍ．ミセリ（Miceli）らによる「逆ハーフ
トーン化（Inverse Halftoning）」、Ｊ．Ｅｌｅｃｔｒ
ｏｎ．Ｉｍａｇｉｎｇ、Ｖｏｌ．１、１４３〜１５１
頁、１９９２年４月、Ｚ．ファン（Fan ）による「デジ
タルハーフトーンからのグレーイメージの検索（Retrie
val of Gray Images From Digital Halftones ）」、Ｉ
ＳＣＡＳ会報、２４７７〜２４８０頁、１９９２年５
月、Ａ．アナロウイ（Analoui ）らによる「ハーフトー
ンからの連続色調の再構築に関する新たな結果（New Re
sults On Reconstruction of Continuous-Tone From Ha
lftone）」、ＩＣＡＳＳＰ会報、ｖｏｌ．ＩＩＩ、３１
３〜３１６頁、１９９２年３月、およびＲ．Ｓ．カーン
（Kern）らによる「線型フィルタリングと繰り返し技術
による非スクリーニング（Descreening Via Linear Fil
teringand Iterative Techniques ）」、ＳＰＩＥ会
報、ｖｏｌ．１９１３、２９９〜３０９頁、１９９３年
２月に開示されている。特に、これらの文献に開示され
た非スクリーニング技術は何れも空間領域において動作
し、インパルス除去（ミセリ）、ロジックフィルタリン
グ（ファン）、コンベックスセットへの投影（アナロウ
イ）、および連続近似（アナロウイ、カーン）の内の種
々のものを使用したローパスフィルタリングを使用して
いる。
トーン画像を非スクリーニングすることに関する他のシ
ステムがＣ．Ｍ．ミセリ（Miceli）らによる「逆ハーフ
トーン化（Inverse Halftoning）」、Ｊ．Ｅｌｅｃｔｒ
ｏｎ．Ｉｍａｇｉｎｇ、Ｖｏｌ．１、１４３〜１５１
頁、１９９２年４月、Ｚ．ファン（Fan ）による「デジ
タルハーフトーンからのグレーイメージの検索（Retrie
val of Gray Images From Digital Halftones ）」、Ｉ
ＳＣＡＳ会報、２４７７〜２４８０頁、１９９２年５
月、Ａ．アナロウイ（Analoui ）らによる「ハーフトー
ンからの連続色調の再構築に関する新たな結果（New Re
sults On Reconstruction of Continuous-Tone From Ha
lftone）」、ＩＣＡＳＳＰ会報、ｖｏｌ．ＩＩＩ、３１
３〜３１６頁、１９９２年３月、およびＲ．Ｓ．カーン
（Kern）らによる「線型フィルタリングと繰り返し技術
による非スクリーニング（Descreening Via Linear Fil
teringand Iterative Techniques ）」、ＳＰＩＥ会
報、ｖｏｌ．１９１３、２９９〜３０９頁、１９９３年
２月に開示されている。特に、これらの文献に開示され
た非スクリーニング技術は何れも空間領域において動作
し、インパルス除去（ミセリ）、ロジックフィルタリン
グ（ファン）、コンベックスセットへの投影（アナロウ
イ）、および連続近似（アナロウイ、カーン）の内の種
々のものを使用したローパスフィルタリングを使用して
いる。
【０００９】ディサリングおよび／または誤差拡散され
たハーフトーン画像のための他の非スクリーニング技術
は、Ｓ．ハイン（Hein）らによる「誤差拡散コード化画
像のハーフトーンから連続色調への変換（Halftone to
Continuous-Tone Conversionof Error-Diffusion Coded
Images）」、ＩＥＥＥ会報、画画像処理、Ｖｏｌ．
４、２０８〜２１６頁、１９９５年２月、Ｓ．Ｍ．シュ
ワイザー（Schweizer ）らによる「色調復元の対するベ
イズ的アプローチ（A Bayesian Approach to Palette R
estoration）」、ＳＰＩＥ会報、ｖｏｌ．１９１３、２
８２〜２９２頁、およびＭ．Ｙ．チン（Ting）らによる
「ハーフトーンからグレースケールへのデコーダと予測
的に切り詰めた木構造ベクトル量子化を使用した誤差拡
散した画像の圧縮（Error Diffused Image Compression
Using a Halftone-to-Grayscale Decoder and Predict
ive Pruned Tree-Structured Vector Quantizatio
n）」、ＩＥＥＥ会報、画画像処理、Ｖｏｌ．３、８５
４〜８５８、１９９４年１１月に開示されている。特
に、ハインは、非スクリーニングの問題を、平均二乗誤
差を最小にする二次計画法の問題として述べている。シ
ュワイザーでは、ディサリングされ誤差拡散されたハー
フトーン画像から連続色調画像を再構築するために統計
的アプローチが適用されている。チンにおいては、誤差
拡散されたハーフトーン画像の非スクリーニングが、ベ
クトル量子化技術に基づきルックアップテーブルを使用
して行われている。しかし、以上概説した非スクリーニ
ング技術の何れも、ハーフトーン画像に関しての情報と
無関係には動作しない。
たハーフトーン画像のための他の非スクリーニング技術
は、Ｓ．ハイン（Hein）らによる「誤差拡散コード化画
像のハーフトーンから連続色調への変換（Halftone to
Continuous-Tone Conversionof Error-Diffusion Coded
Images）」、ＩＥＥＥ会報、画画像処理、Ｖｏｌ．
４、２０８〜２１６頁、１９９５年２月、Ｓ．Ｍ．シュ
ワイザー（Schweizer ）らによる「色調復元の対するベ
イズ的アプローチ（A Bayesian Approach to Palette R
estoration）」、ＳＰＩＥ会報、ｖｏｌ．１９１３、２
８２〜２９２頁、およびＭ．Ｙ．チン（Ting）らによる
「ハーフトーンからグレースケールへのデコーダと予測
的に切り詰めた木構造ベクトル量子化を使用した誤差拡
散した画像の圧縮（Error Diffused Image Compression
Using a Halftone-to-Grayscale Decoder and Predict
ive Pruned Tree-Structured Vector Quantizatio
n）」、ＩＥＥＥ会報、画画像処理、Ｖｏｌ．３、８５
４〜８５８、１９９４年１１月に開示されている。特
に、ハインは、非スクリーニングの問題を、平均二乗誤
差を最小にする二次計画法の問題として述べている。シ
ュワイザーでは、ディサリングされ誤差拡散されたハー
フトーン画像から連続色調画像を再構築するために統計
的アプローチが適用されている。チンにおいては、誤差
拡散されたハーフトーン画像の非スクリーニングが、ベ
クトル量子化技術に基づきルックアップテーブルを使用
して行われている。しかし、以上概説した非スクリーニ
ング技術の何れも、ハーフトーン画像に関しての情報と
無関係には動作しない。
【００１０】更に、全てではないものの、以上概説した
非スクリーニングの多くのものは、二進ハーフトーン値
の空間分布に基づいて動作する。非スクリーニングの実
際的な使用の一つは、新聞、雑誌、本などの書類の印刷
されたハーフトーン画像を走査して得た電子的なハーフ
トーン画像から電子的なコントーン画像を復元すること
である。この場合、電子的なハーフトーン画像を作るた
めにハーフトーン画像が走査される。しかし、ハーフト
ーンの空間解画像度がスキャナの空間解画像度と正確に
一致している可能性は非常に低いため、スキャナは電子
的ハーフトーン画像を作る時にハーフトーン画像を空間
的にフィルタにかけてしまう。したがって、走査したハ
ーフトーン画像を非スクリーニングする時、非スクリー
ニング技術は、ハーフトーン化プロセスから得られたハ
ーフトーン値を、走査中にハーフトーン画像を空間的に
フィルタにかけることによって得られた値と区別する
か、さもなければそれを分離しなければならない。
非スクリーニングの多くのものは、二進ハーフトーン値
の空間分布に基づいて動作する。非スクリーニングの実
際的な使用の一つは、新聞、雑誌、本などの書類の印刷
されたハーフトーン画像を走査して得た電子的なハーフ
トーン画像から電子的なコントーン画像を復元すること
である。この場合、電子的なハーフトーン画像を作るた
めにハーフトーン画像が走査される。しかし、ハーフト
ーンの空間解画像度がスキャナの空間解画像度と正確に
一致している可能性は非常に低いため、スキャナは電子
的ハーフトーン画像を作る時にハーフトーン画像を空間
的にフィルタにかけてしまう。したがって、走査したハ
ーフトーン画像を非スクリーニングする時、非スクリー
ニング技術は、ハーフトーン化プロセスから得られたハ
ーフトーン値を、走査中にハーフトーン画像を空間的に
フィルタにかけることによって得られた値と区別する
か、さもなければそれを分離しなければならない。
【００１１】したがって、非スクリーニングされるべき
ハーフトーンを作ったハーフトーン化プロセスに関する
知識（ハーフトーン化プロセスの形式など）とは無関係
に動作する汎用の非スクリーニング技術に対するニーズ
がある。さらに、非スクリーニングされるべきハーフト
ーン画像が走査によって得たものであるか否かとは無関
係に動作する汎用の非スクリーニング技術に対するニー
ズもある。
ハーフトーンを作ったハーフトーン化プロセスに関する
知識（ハーフトーン化プロセスの形式など）とは無関係
に動作する汎用の非スクリーニング技術に対するニーズ
がある。さらに、非スクリーニングされるべきハーフト
ーン画像が走査によって得たものであるか否かとは無関
係に動作する汎用の非スクリーニング技術に対するニー
ズもある。
【００１２】したがって、本発明は、任意のハーフトー
ン画像を作ったハーフトーン化プロセスに関する知識
（ハーフトーン化プロセスの形式など）とは無関係に任
意のハーフトーン画像を非スクリーニングするための非
スクリーニング方法および装置を提供する。
ン画像を作ったハーフトーン化プロセスに関する知識
（ハーフトーン化プロセスの形式など）とは無関係に任
意のハーフトーン画像を非スクリーニングするための非
スクリーニング方法および装置を提供する。
【００１３】また、本発明は、ハーフトーン画像が走査
によって得たものであるか否かに関わらず任意のハーフ
トーン画像を非スクリーニングするための非スクリーニ
ング方法および装置を提供する。
によって得たものであるか否かに関わらず任意のハーフ
トーン画像を非スクリーニングするための非スクリーニ
ング方法および装置を提供する。
【００１４】このために、本発明は、ウエーブレットに
基づく汎用の非スクリーニング方法と装置を提供する。
基づく汎用の非スクリーニング方法と装置を提供する。
【００１５】さらに、本発明は、ハーフトーン雑音を除
去して、元の画像の内容の殆どを保存する一連の空間的
および周波数的なフィルタを使用してハーフトーン画像
を選択的に処理するための方法と装置を提供する。
去して、元の画像の内容の殆どを保存する一連の空間的
および周波数的なフィルタを使用してハーフトーン画像
を選択的に処理するための方法と装置を提供する。
【００１６】また、本発明は、得られたコントーン画像
を非線型的にフィルタ処理して最終的なコントーン画像
を作るための方法と装置を提供する。
を非線型的にフィルタ処理して最終的なコントーン画像
を作るための方法と装置を提供する。
【００１７】さらに、本発明は、パラメータ予測とは無
関係であり、このため走査によって得たハーフトーン画
像などの全ての形式のハーフトーン画像に対して汎用的
に使用できる方法と装置を提供する。
関係であり、このため走査によって得たハーフトーン画
像などの全ての形式のハーフトーン画像に対して汎用的
に使用できる方法と装置を提供する。
【００１８】
【課題を解決するための手段】本発明の請求項１の態様
に従うと、ハーフトーン画像から連続色調画像を作るた
めの方法であって、ハーフトーン画像を複数のサブバン
ド画像に分解する工程と、この複数のサブバンド画像か
ら雑音とパターンの内の少なくとも一つを選択的に除去
する工程と、複数のサブバンド画像を連続色調画像に再
構成する工程とを含むことを特徴とする。
に従うと、ハーフトーン画像から連続色調画像を作るた
めの方法であって、ハーフトーン画像を複数のサブバン
ド画像に分解する工程と、この複数のサブバンド画像か
ら雑音とパターンの内の少なくとも一つを選択的に除去
する工程と、複数のサブバンド画像を連続色調画像に再
構成する工程とを含むことを特徴とする。
【００１９】特に、本発明の方法は、元のハーフトーン
画像を複数のサブバンド画像に分解する工程と、この複
数のサブバンド画像の各々の縁でない領域から雑音を減
衰させる工程と、この雑音が減衰された複数のサブバン
ド画像の各々に対して方向性（oriented）フィルタ処理
を行う工程と、この方向性フィルタ処理され、雑音が減
衰された複数のサブバンド画像を初期コントーン画像に
再構成する工程と、非線型の縁保存型フィルタを使用し
て初期のコントーン画像をフィルタ処理し、最終的なコ
ントーン画像を形成する工程とを含む。
画像を複数のサブバンド画像に分解する工程と、この複
数のサブバンド画像の各々の縁でない領域から雑音を減
衰させる工程と、この雑音が減衰された複数のサブバン
ド画像の各々に対して方向性（oriented）フィルタ処理
を行う工程と、この方向性フィルタ処理され、雑音が減
衰された複数のサブバンド画像を初期コントーン画像に
再構成する工程と、非線型の縁保存型フィルタを使用し
て初期のコントーン画像をフィルタ処理し、最終的なコ
ントーン画像を形成する工程とを含む。
【００２０】また、本発明による方法の分解工程は、離
散ウエーブレット変換（discretewavelet transform ）
に基づき、元のハーフトーン画像を複数の異なる周波数
のサブバンド画像に分解する工程を含む。この分解工程
は、複数の上位レベルの異なる周波数のサブバンド画像
の一つを複数の下位レベルの異なる周波数のサブバンド
画像に分解する工程と、この上位レベルサブバンド画像
分解工程を複数回繰り返す工程とを更に含む。
散ウエーブレット変換（discretewavelet transform ）
に基づき、元のハーフトーン画像を複数の異なる周波数
のサブバンド画像に分解する工程を含む。この分解工程
は、複数の上位レベルの異なる周波数のサブバンド画像
の一つを複数の下位レベルの異なる周波数のサブバンド
画像に分解する工程と、この上位レベルサブバンド画像
分解工程を複数回繰り返す工程とを更に含む。
【００２１】雑音減衰工程は、高周波数サブバンド画像
の画素の大きさをクリップする工程を含み、さらに各レ
ベルのサブバンド画像のについて親サブバンド画像の大
きさに重み付けをする工程と、高周波数サブバンド画像
の大きさを、それらの親サブバンド画像の重みの付けら
れた大きさまでの範囲でクリップする工程とを含む。さ
らにまた、重み付けの値はレベル毎に変化する。また、
このクリップ工程は、高周波数サブバンド画像の縁に適
応される。
の画素の大きさをクリップする工程を含み、さらに各レ
ベルのサブバンド画像のについて親サブバンド画像の大
きさに重み付けをする工程と、高周波数サブバンド画像
の大きさを、それらの親サブバンド画像の重みの付けら
れた大きさまでの範囲でクリップする工程とを含む。さ
らにまた、重み付けの値はレベル毎に変化する。また、
このクリップ工程は、高周波数サブバンド画像の縁に適
応される。
【００２２】方向性フィルタ処理工程は、一次元ローパ
スフィルタを各サブバンド画像の方位に沿って適用する
工程を含む。すなわち、垂直サブバンド画像に対しては
水平ローパスフィルタを適用し、水平サブバンド画像に
対しては垂直ローパスフィルタを適用し、傾斜したサブ
バンド画像に対してはＸ形状のローパスフィルタを適用
する。また、この方向性フィルタ処理工程は、サブバン
ド画像の解画像度にフィルタ処理を適応させる工程を含
む。加えて、この方向性フィルタ処理工程は、スクリー
ンパラメータの推定値に基づいて適切なフィルタを選択
する工程を含む。
スフィルタを各サブバンド画像の方位に沿って適用する
工程を含む。すなわち、垂直サブバンド画像に対しては
水平ローパスフィルタを適用し、水平サブバンド画像に
対しては垂直ローパスフィルタを適用し、傾斜したサブ
バンド画像に対してはＸ形状のローパスフィルタを適用
する。また、この方向性フィルタ処理工程は、サブバン
ド画像の解画像度にフィルタ処理を適応させる工程を含
む。加えて、この方向性フィルタ処理工程は、スクリー
ンパラメータの推定値に基づいて適切なフィルタを選択
する工程を含む。
【００２３】初期コントーン画像をフィルタ処理する工
程は、フーバー・マルコフ（HuberMarkov）のランダム
フィールド（ＨＭＲＦ）を使用してモデル化した空間平
滑度拘束を使用して初期コントーン画像を非線型フィル
タ処理する工程を含み、またσ−フィルタに基づくフィ
ルタ処理などの縁保存法を使用して初期のコントーンを
非線型的にフィルタ処理する工程を、ＨＭＲＦ空間平滑
度拘束の代わりに、またはそれに加えて含んでいてもよ
い。
程は、フーバー・マルコフ（HuberMarkov）のランダム
フィールド（ＨＭＲＦ）を使用してモデル化した空間平
滑度拘束を使用して初期コントーン画像を非線型フィル
タ処理する工程を含み、またσ−フィルタに基づくフィ
ルタ処理などの縁保存法を使用して初期のコントーンを
非線型的にフィルタ処理する工程を、ＨＭＲＦ空間平滑
度拘束の代わりに、またはそれに加えて含んでいてもよ
い。
【００２４】上記のおよびその他の本発明の特徴と利点
は、好適な実施の形態についての以下の詳細な説明に記
載されるか、または詳細な説明から明らかである。
は、好適な実施の形態についての以下の詳細な説明に記
載されるか、または詳細な説明から明らかである。
【００２５】
【発明の実施の形態】本発明は、図１および図２に示す
ようなハーフトーン画像を、図２３および図２４に示す
ような連続色調のすなわち「コントーン」な画像に変換
することに関するものである。図１は４００×４００画
素の誤差拡散ハーフトーン画像を示す。図２は３×３集
合ドットハーフトーン化によって形成した同じ画像を示
す。
ようなハーフトーン画像を、図２３および図２４に示す
ような連続色調のすなわち「コントーン」な画像に変換
することに関するものである。図１は４００×４００画
素の誤差拡散ハーフトーン画像を示す。図２は３×３集
合ドットハーフトーン化によって形成した同じ画像を示
す。
【００２６】図３に示すように、これらの画像は、バイ
ナリープリンタなどによる印刷のためにハーフトーン画
像を発生するようにプログラムされたコンピュータなど
のハーフトーン画像発生器２００によって形成されたも
のである。汎用ハーフトーン画像非スクリーニング装置
１００へ直接入力されるであろう電子的データとして形
成され記憶される画像を電子的ハーフトーン画像と呼
ぶ。一般に、ハーフトーン画像発生器３００は任意のハ
ーフトーン画像のソースでよい。
ナリープリンタなどによる印刷のためにハーフトーン画
像を発生するようにプログラムされたコンピュータなど
のハーフトーン画像発生器２００によって形成されたも
のである。汎用ハーフトーン画像非スクリーニング装置
１００へ直接入力されるであろう電子的データとして形
成され記憶される画像を電子的ハーフトーン画像と呼
ぶ。一般に、ハーフトーン画像発生器３００は任意のハ
ーフトーン画像のソースでよい。
【００２７】これに対し、図２５に示すように、ハーフ
トーン画像を、バイナリープリンタで印刷してその後ス
キャナ３００で走査し、これによって汎用ハーフトーン
画像非スクリーニング装置１００のための入力データを
発生することができる。図２５に示す走査によって得た
ハーフトーン画像では、元の画像の解画像度は４００×
４００画素であった。この元の画像を２倍にし、４×８
集合ドットを用いてハーフトーン化した。このハーフト
ーン画像を３００ｓｐｉで印刷し、その後２００ｓｐｉ
で走査した。
トーン画像を、バイナリープリンタで印刷してその後ス
キャナ３００で走査し、これによって汎用ハーフトーン
画像非スクリーニング装置１００のための入力データを
発生することができる。図２５に示す走査によって得た
ハーフトーン画像では、元の画像の解画像度は４００×
４００画素であった。この元の画像を２倍にし、４×８
集合ドットを用いてハーフトーン化した。このハーフト
ーン画像を３００ｓｐｉで印刷し、その後２００ｓｐｉ
で走査した。
【００２８】何れの場合も、処理すべきハーフトーン画
像が元の電子的ハーフトーン画像であるか走査されたハ
ーフトーン画像であるによって、画像はハーフトーン画
像発生器２００またはスキャナ３００から入力ライン１
０３を介して汎用ハーフトーン画像非スクリーニング装
置１００のハーフトーン画像入力システム１３０に入力
される。
像が元の電子的ハーフトーン画像であるか走査されたハ
ーフトーン画像であるによって、画像はハーフトーン画
像発生器２００またはスキャナ３００から入力ライン１
０３を介して汎用ハーフトーン画像非スクリーニング装
置１００のハーフトーン画像入力システム１３０に入力
される。
【００２９】一般に、本発明の一つの実施の形態におい
ては、汎用ハーフトーン画像非スクリーニング装置１０
０は、コントローラ１１０、メモリ１２０、ハーフトー
ン画像入力システム１３０、Ｎステージハーフトーン画
像デコンポーザ１４０、サブバンド画像雑音および／ま
たはパターンフィルタシステム１５０、サブバンドから
連続色調画像を再構成するＮステージリコンポーザ１６
０、連続色調画像フィルタシステム１７０、および連続
色調画像出力システム１８０を含む。これらの要素の各
々は制御バス１９０とデータバス１９５によって相互接
続されている。汎用ハーフトーン画像非スクリーニング
装置１００は、プログラムされた汎用コンピュータ上で
実現するのが好ましい。しかし、汎用ハーフトーン画像
非スクリーニング装置１００は専用コンピュータ、プロ
グラムされたマイクロプロセッサまたはマイクロコント
ローラと周辺集積回路要素、およびＡＳＩＣまたはその
他の集積回路、ディスクリート回路要素などのハードウ
エアによる電子的または論理的回路、ＰＬＤ、ＰＬＡま
たはＰＡＬなどのプログラマブル論理デバイスなどによ
って実現することができる。一般に、図１１〜図１６に
示すフローチャートを実施できる任意のデバイスを、汎
用ハーフトーン画像非スクリーニング装置１００を実施
するために使用できる。
ては、汎用ハーフトーン画像非スクリーニング装置１０
０は、コントローラ１１０、メモリ１２０、ハーフトー
ン画像入力システム１３０、Ｎステージハーフトーン画
像デコンポーザ１４０、サブバンド画像雑音および／ま
たはパターンフィルタシステム１５０、サブバンドから
連続色調画像を再構成するＮステージリコンポーザ１６
０、連続色調画像フィルタシステム１７０、および連続
色調画像出力システム１８０を含む。これらの要素の各
々は制御バス１９０とデータバス１９５によって相互接
続されている。汎用ハーフトーン画像非スクリーニング
装置１００は、プログラムされた汎用コンピュータ上で
実現するのが好ましい。しかし、汎用ハーフトーン画像
非スクリーニング装置１００は専用コンピュータ、プロ
グラムされたマイクロプロセッサまたはマイクロコント
ローラと周辺集積回路要素、およびＡＳＩＣまたはその
他の集積回路、ディスクリート回路要素などのハードウ
エアによる電子的または論理的回路、ＰＬＤ、ＰＬＡま
たはＰＡＬなどのプログラマブル論理デバイスなどによ
って実現することができる。一般に、図１１〜図１６に
示すフローチャートを実施できる任意のデバイスを、汎
用ハーフトーン画像非スクリーニング装置１００を実施
するために使用できる。
【００３０】ハーフトーン画像入力システム１３０は、
ネットワークへのＬＡＮ接続、公衆電話回線網を介して
大型コンピュータ、マイクロコンピュータ、ミニコンピ
ュータ、またはパーソナルコンピュータに接続されたモ
デム、スキャナ３００に接続されたパラレルポートまた
はシリアルポートなどであることができる。ハーフトー
ン画像データがハーフトーン画像入力システム１３０を
介して汎用ハーフトーン画像非スクリーニング装置１０
０に入力されると、ハーフトーン画像データはメモリ１
２０に記憶される。
ネットワークへのＬＡＮ接続、公衆電話回線網を介して
大型コンピュータ、マイクロコンピュータ、ミニコンピ
ュータ、またはパーソナルコンピュータに接続されたモ
デム、スキャナ３００に接続されたパラレルポートまた
はシリアルポートなどであることができる。ハーフトー
ン画像データがハーフトーン画像入力システム１３０を
介して汎用ハーフトーン画像非スクリーニング装置１０
０に入力されると、ハーフトーン画像データはメモリ１
２０に記憶される。
【００３１】メモリ１２０は、ＲＡＭ、フラッシュメモ
リ、ＰＲＯＭ、ＥＰＲＯＭおよびＥＥＰＲＯＭなどの書
込み可能なメモリとＲＯＭなどの書込みの不可能なメモ
リの両方を含むことができる。ＲＡＭはスタティックま
たはダイナミックＲＡＭであることが好ましい。しか
し、メモリ１２０は、フロッピーディスクとディスクド
ライブ、書込み可能型光ディスクとディスクドライブ、
ハードディスクなどを使用して実現できる。
リ、ＰＲＯＭ、ＥＰＲＯＭおよびＥＥＰＲＯＭなどの書
込み可能なメモリとＲＯＭなどの書込みの不可能なメモ
リの両方を含むことができる。ＲＡＭはスタティックま
たはダイナミックＲＡＭであることが好ましい。しか
し、メモリ１２０は、フロッピーディスクとディスクド
ライブ、書込み可能型光ディスクとディスクドライブ、
ハードディスクなどを使用して実現できる。
【００３２】元のハーフトーン画像がメモリ１２０に記
憶されると、それはＮステージハーフトーン画像デコン
ポーザ１４０に出力され、そこで複数のサブバンド画像
レベルに分解される。各サブバンド画像レベルは複数の
サブバンド画像を含んでいる。元のハーフトーン画像が
Ｎステージハーフトーン画像デコンポーザ１４０によっ
てサブバンドに一旦分解されると、このサブバンド画像
はメモリ１２０に記憶される。その後、サブバンド画像
からハーフトーンエネルギーを濾去するためにサブバン
ド画像はサブバンド雑音および／またはパターンフィル
タシステム１５０に出力される。フィルタ処理されたサ
ブバンド画像はその後メモリ１２０に記憶される。
憶されると、それはＮステージハーフトーン画像デコン
ポーザ１４０に出力され、そこで複数のサブバンド画像
レベルに分解される。各サブバンド画像レベルは複数の
サブバンド画像を含んでいる。元のハーフトーン画像が
Ｎステージハーフトーン画像デコンポーザ１４０によっ
てサブバンドに一旦分解されると、このサブバンド画像
はメモリ１２０に記憶される。その後、サブバンド画像
からハーフトーンエネルギーを濾去するためにサブバン
ド画像はサブバンド雑音および／またはパターンフィル
タシステム１５０に出力される。フィルタ処理されたサ
ブバンド画像はその後メモリ１２０に記憶される。
【００３３】そして、記憶されたフィルタ処理済のサブ
バンド画像は、サブバンドから連続色調画像を再構成す
るＮステージリコンポーザ１６０に出力される。このリ
コンポーザ１６０は、サブバンド画像をコントーンの初
期単一画像に再構成する。その後、この初期コントーン
画像はメモリ１２０に記憶される。そして、初期コント
ーン画像はメモリ１２０から読み出され、連続色調画像
フィルタシステム１７０に出力される。この連続色調画
像フィルタシステム１７０は元のコントーン画像をフィ
ルタ処理し、強調されたコントーン画像を形成する。初
期コントーン画像は、空間平滑度拘束を有する非線型フ
ィルタまたはその他の縁保存フィルタリング方法を使用
して強調し、初期コントーンを平滑にするとともにフィ
ルタ化または再構成による人為的結果（artifact）を除
去することが好ましい。最終コントーン画像はその後メ
モリ１２０に記憶される。
バンド画像は、サブバンドから連続色調画像を再構成す
るＮステージリコンポーザ１６０に出力される。このリ
コンポーザ１６０は、サブバンド画像をコントーンの初
期単一画像に再構成する。その後、この初期コントーン
画像はメモリ１２０に記憶される。そして、初期コント
ーン画像はメモリ１２０から読み出され、連続色調画像
フィルタシステム１７０に出力される。この連続色調画
像フィルタシステム１７０は元のコントーン画像をフィ
ルタ処理し、強調されたコントーン画像を形成する。初
期コントーン画像は、空間平滑度拘束を有する非線型フ
ィルタまたはその他の縁保存フィルタリング方法を使用
して強調し、初期コントーンを平滑にするとともにフィ
ルタ化または再構成による人為的結果（artifact）を除
去することが好ましい。最終コントーン画像はその後メ
モリ１２０に記憶される。
【００３４】最終コントーン画像を印刷する時、最終コ
ントーン画像がメモリ１２０から連続色調画像出力シス
テム１８０に出力される。連続色調画像出力システム１
８０は、最終コントーン画像を、プリンタ４００、ハー
フトーン画像発生器２００、またはコンピュータ、デー
タ記憶装置などのその他の装置に出力ライン１０６を介
して出力する。
ントーン画像がメモリ１２０から連続色調画像出力シス
テム１８０に出力される。連続色調画像出力システム１
８０は、最終コントーン画像を、プリンタ４００、ハー
フトーン画像発生器２００、またはコンピュータ、デー
タ記憶装置などのその他の装置に出力ライン１０６を介
して出力する。
【００３５】さらに、印刷する代わりに、最終コントー
ン画像を更に処理することもできる。この追加の処理に
は、圧縮、画画像認識、異なるハーフトーンプロセスま
たは異なるハーフトーン化パラメータを用いた再スクリ
ーニング、および種々の画像操作処理を含むことができ
る。これらの種々の画像操作処理には、拡大、縮小、フ
リッピング（flipping）、回転などが含まれる。このよ
うな画画像操作処理は元のハーフトーン画像に適用でき
るが、一般にそれらは非スクリーニングされたコントー
ン画像により簡単にまたは全体的に適用できる。例え
ば、広い範囲の回転角度をコントーン画像に使用でき
る。
ン画像を更に処理することもできる。この追加の処理に
は、圧縮、画画像認識、異なるハーフトーンプロセスま
たは異なるハーフトーン化パラメータを用いた再スクリ
ーニング、および種々の画像操作処理を含むことができ
る。これらの種々の画像操作処理には、拡大、縮小、フ
リッピング（flipping）、回転などが含まれる。このよ
うな画画像操作処理は元のハーフトーン画像に適用でき
るが、一般にそれらは非スクリーニングされたコントー
ン画像により簡単にまたは全体的に適用できる。例え
ば、広い範囲の回転角度をコントーン画像に使用でき
る。
【００３６】したがって、連続色調画像出力システム１
８０は、パラレルまたはシリアルのプリンタポート、Ｌ
ＡＮ接続、モデムなどでよい。同様に、コントーン画像
を受け取る装置は、上記した機能の何れかを行える任意
の装置でよい。
８０は、パラレルまたはシリアルのプリンタポート、Ｌ
ＡＮ接続、モデムなどでよい。同様に、コントーン画像
を受け取る装置は、上記した機能の何れかを行える任意
の装置でよい。
【００３７】したがって、汎用ハーフトーン画像非スク
リーニング装置１００、連続色調画像出力システム１８
０、およびハーフトーン画像入力システム１３０は、実
際には同一の装置であり、入力ライン１０３と出力ライ
ン１０６は、実際には同一のラインである。また、ハー
フトーン画像入力システム１３０、Ｎステージハーフト
ーン画像デコンポーザ１４０、サブバンド画像雑音およ
び／またはパターンフィルタシステム１５０、サブバン
ドから連続色調画像を再構成するＮステージリコンポー
ザ１６０、または連続色調画像フィルタシステム１７０
の各々からの処理された画像出力を記憶することは絶対
に必要なものでない。各サブシステムから出力されるサ
ブバンド画像またはコントーン画像は次の下流側のサブ
システムに直接入力することができる。
リーニング装置１００、連続色調画像出力システム１８
０、およびハーフトーン画像入力システム１３０は、実
際には同一の装置であり、入力ライン１０３と出力ライ
ン１０６は、実際には同一のラインである。また、ハー
フトーン画像入力システム１３０、Ｎステージハーフト
ーン画像デコンポーザ１４０、サブバンド画像雑音およ
び／またはパターンフィルタシステム１５０、サブバン
ドから連続色調画像を再構成するＮステージリコンポー
ザ１６０、または連続色調画像フィルタシステム１７０
の各々からの処理された画像出力を記憶することは絶対
に必要なものでない。各サブシステムから出力されるサ
ブバンド画像またはコントーン画像は次の下流側のサブ
システムに直接入力することができる。
【００３８】また、図３ではサブシステムの各々が別々
のものとして示されているが、プログラム式汎用コンピ
ュータ、専用コンピュータ、マイクロプロセッサ、マイ
クロコントローラ、ＡＳＩＣ、または本発明の汎用ハー
フトーン画像非スクリーニング装置１００を実施できる
その他の装置においては、これらのサブシステムを実現
する実際の回路要素は、汎用ハーフトーン画像非スクリ
ーニング装置１００のサブシステム１１０、１３０〜１
８０の全てとはいかないまでも、それらの内の複数のも
のに対して共通である可能性が高い。
のものとして示されているが、プログラム式汎用コンピ
ュータ、専用コンピュータ、マイクロプロセッサ、マイ
クロコントローラ、ＡＳＩＣ、または本発明の汎用ハー
フトーン画像非スクリーニング装置１００を実施できる
その他の装置においては、これらのサブシステムを実現
する実際の回路要素は、汎用ハーフトーン画像非スクリ
ーニング装置１００のサブシステム１１０、１３０〜１
８０の全てとはいかないまでも、それらの内の複数のも
のに対して共通である可能性が高い。
【００３９】本発明の実施の形態では、離散ウエーブレ
ット変換を用いてハーフトーン画像を異なる周波数のサ
ブバンド画像に分解している。離散ウエーブレット変換
とその逆変換は、Ｍ．ベターリ（Vetterli）とＪ．コバ
セビック（Kovacevic）による「ウエーブレットとサブ
バンドコーディグ（Wavelets and Subband Coding
）」、Ｐｒｅｎｔｉｃｅ Ｈａｌｌ ＰＴＲ（１９９
５年）、Ｇ．ストラング（Strang）とＴ．ニューエン
（Nguyen）による「ウエーブレットとフィルタバンク
（Wavelets and Filter Banks ）」、Ｗｅｌｌｅｓｌｅ
ｙ−ＣａｍｂｒｉｄｇｅＰｒｅｓｓ（１９９６年）、お
よびＹ．Ｔ．チャン（Chan）による「ウエーブレットの
基礎（Wavelet Basics )」、Ｋｌｕｗｅｒ Ａｃａｄｅ
ｍｉｃ Ｐｕｂｌｉｓｈｅｒｓ（１９９５年）に完全に
記載されており、これらを本明細書の一部を構成するも
のとして参照する。離散ウエーブレット変換は、得られ
るサブバンド画像において非スクリーニングをより高い
効率で行えるために使用される。ハーフトーン化プロセ
スのエネルギーは主に高周波数のサブバンド画像に存在
するが、周波数が低く、解画像度が粗いサブバンド画像
は信号エネルギーを主に含んでいる。
ット変換を用いてハーフトーン画像を異なる周波数のサ
ブバンド画像に分解している。離散ウエーブレット変換
とその逆変換は、Ｍ．ベターリ（Vetterli）とＪ．コバ
セビック（Kovacevic）による「ウエーブレットとサブ
バンドコーディグ（Wavelets and Subband Coding
）」、Ｐｒｅｎｔｉｃｅ Ｈａｌｌ ＰＴＲ（１９９
５年）、Ｇ．ストラング（Strang）とＴ．ニューエン
（Nguyen）による「ウエーブレットとフィルタバンク
（Wavelets and Filter Banks ）」、Ｗｅｌｌｅｓｌｅ
ｙ−ＣａｍｂｒｉｄｇｅＰｒｅｓｓ（１９９６年）、お
よびＹ．Ｔ．チャン（Chan）による「ウエーブレットの
基礎（Wavelet Basics )」、Ｋｌｕｗｅｒ Ａｃａｄｅ
ｍｉｃ Ｐｕｂｌｉｓｈｅｒｓ（１９９５年）に完全に
記載されており、これらを本明細書の一部を構成するも
のとして参照する。離散ウエーブレット変換は、得られ
るサブバンド画像において非スクリーニングをより高い
効率で行えるために使用される。ハーフトーン化プロセ
スのエネルギーは主に高周波数のサブバンド画像に存在
するが、周波数が低く、解画像度が粗いサブバンド画像
は信号エネルギーを主に含んでいる。
【００４０】ウエーブレット分解は、周波数領域および
空間領域の何れにおいても処理を容易にする。高周波数
サブバンド画像では、信号エネルギーとハーフトーン化
エネルギーが大量に混じっている。また、高周波数サブ
バンド画像は、縁位置における重要な縁情報に対応する
信号エネルギーを含んでいる。分離したフィルタバンク
を使用することにより、クアッドツリー（quadtree）構
造のサブバンド画像に分解される。このクアッドツリー
構造のサブバンド画像の各々は異なった方位を有する。
所定のサブバンド画像レベルの各々について、３個の高
周波数サブバンド画像、（“ＬＨ”サブバンド画像、
“ＨＬ”サブバンド画像、および“ＨＨ”サブバンド画
像）が存在し、これらは水平方位、垂直方位、傾斜方位
にそれぞれ対応している。解画像度が最も粗いレベルで
は、サブバンド画像は、“ＬＬ”サブバンド画像と呼ば
れるローパス周波数サブバンド画像も含んでいる。
空間領域の何れにおいても処理を容易にする。高周波数
サブバンド画像では、信号エネルギーとハーフトーン化
エネルギーが大量に混じっている。また、高周波数サブ
バンド画像は、縁位置における重要な縁情報に対応する
信号エネルギーを含んでいる。分離したフィルタバンク
を使用することにより、クアッドツリー（quadtree）構
造のサブバンド画像に分解される。このクアッドツリー
構造のサブバンド画像の各々は異なった方位を有する。
所定のサブバンド画像レベルの各々について、３個の高
周波数サブバンド画像、（“ＬＨ”サブバンド画像、
“ＨＬ”サブバンド画像、および“ＨＨ”サブバンド画
像）が存在し、これらは水平方位、垂直方位、傾斜方位
にそれぞれ対応している。解画像度が最も粗いレベルで
は、サブバンド画像は、“ＬＬ”サブバンド画像と呼ば
れるローパス周波数サブバンド画像も含んでいる。
【００４１】異なる解画像度レベルのサブバンド画像の
間には親子関係が存在する。この関係は“ツリー”と呼
ばれる。この関係が図１０に示されている。サブバンド
画像の各係数は、次のより細い解画像度の同じ方位のサ
ブバンド画像において同じ空間位置にある４個の子供を
有する。ただし、各係数が３個の子供しか有しないロー
パスＬＬサブバンド画像は例外である。このように、レ
ベル全体を通して親子関係が存在している。一般に、子
孫の係数の大きさは、それらの親に対して増加すること
はない。すなわち、画像のパワースペクトルは、周波数
が増加するにつれて衰える。したがって、もし係数がそ
の大きさに関して重要でないことが分かった場合、その
子供も重要でない可能性が高い。
間には親子関係が存在する。この関係は“ツリー”と呼
ばれる。この関係が図１０に示されている。サブバンド
画像の各係数は、次のより細い解画像度の同じ方位のサ
ブバンド画像において同じ空間位置にある４個の子供を
有する。ただし、各係数が３個の子供しか有しないロー
パスＬＬサブバンド画像は例外である。このように、レ
ベル全体を通して親子関係が存在している。一般に、子
孫の係数の大きさは、それらの親に対して増加すること
はない。すなわち、画像のパワースペクトルは、周波数
が増加するにつれて衰える。したがって、もし係数がそ
の大きさに関して重要でないことが分かった場合、その
子供も重要でない可能性が高い。
【００４２】多くの画像については、この仮定は一般的
に有効である。ツリー拘束に違反しているということ
は、ハーフトーン化エネルギーが存在していることを示
す場合が多い。すなわち、ハーフトーンエネルギーの存
在は、画像信号のスペクトル特性と矛盾する。したがっ
て、高周波サブバンド画像の係数はクリップされ、それ
らの大きさはそれら親の係数の大きさ以下である。この
ツリーに基づくクリッピングは、解画像度が粗い高周波
サブバンド画像には適用されない。なぜなら、これらの
サブバンド画像の殆どは一般に信号情報を含んでおり、
高周波数の親サブバンド画像を有していないからであ
る。
に有効である。ツリー拘束に違反しているということ
は、ハーフトーン化エネルギーが存在していることを示
す場合が多い。すなわち、ハーフトーンエネルギーの存
在は、画像信号のスペクトル特性と矛盾する。したがっ
て、高周波サブバンド画像の係数はクリップされ、それ
らの大きさはそれら親の係数の大きさ以下である。この
ツリーに基づくクリッピングは、解画像度が粗い高周波
サブバンド画像には適用されない。なぜなら、これらの
サブバンド画像の殆どは一般に信号情報を含んでおり、
高周波数の親サブバンド画像を有していないからであ
る。
【００４３】したがって、ハーフトーンエネルギーの殆
どは、縁の回りを除いてツリーに基づくクリッピングプ
ロセスにおいて抑制されるか濾去される。各高周波サブ
バンド画像は、主にその方位に沿った情報を含んでいる
ため、上記のプロセスは、画像データから残存するハー
フトーン雑音を更に濾去するために使用できる。垂直ま
たは水平に方向づけられたサブバンド画像に対してはそ
れらの方位に沿って一次元のローパスフィルタが適用さ
れる。傾いて方向づけられたサブバンド画像に対しては
Ｘ形状のローパスフィルタが適用される。すなわち、傾
斜したサブバンド画像は、より厳しくローパスフィルタ
処理される。なぜなら、それらは一般に、含んでいる信
号情報が少なく、より多くのハーフトーン雑音を含んで
いるからである。この第２レベルのフィルタ処理は、ロ
ーパスＬＬサブバンド画像と同じ粗さの解画像度のレベ
ルのサブバンド画像を含め、全ての高周波サブバンド画
像に適用される。方向性フィルタ処理は、例えば、Ｌポ
イントランニング平均化であることが好ましい。 一般
に、２つのサブバンド画像フィルタリング工程を経た初
期の再構成されたコントーン画像は、信号情報の消失に
起因して幾分曖昧になるかまたはむらが生じている。特
に、高周波信号情報の消失により、“リンギング（ring
ing ）”アーティファクト（artifact）が生じる。した
がって、再構築された連続色調画像は、マルコフのラン
ダムフィールドによってモデル化した空間平滑度拘束を
有する非線型フィルタ処理を用いて強調処理することが
好ましい。マルコフのランダムフィールドは、フーバー
ミニマックス関数と呼ばれる特定のギブズポテンシャル
関数を使用して組み込まれる。この関数が図１８に示さ
れている。この関数の二次セグメントは閾値Ｔ以下での
不連続部分の平滑化を容易にし、直線セグメントは鋭い
縁部を保存する。この様な画像強調技術は、アーティフ
ァクトを大幅に減少させ、元の画像の細部を保存する。
σ−フィルタに基づくフィルタ処理などの他の縁保存フ
ィルタ処理方法も使用できる。σ−フィルタに基づくフ
ィルタ処理は、Ｊ．リー（Lee ）による「デジタルイメ
ージスムージングおよびσ−フィルタ（Digital Image
Smoothing and the σ-Filter ）、Ｃｏｍｐｕｔｅｒ
Ｉｍａｇｅ Ｇｒａｐｈｉｃｓ ａｎｄ Ｉｍａｇｅ
Ｐｒｏｃｅｓｓｉｎｇ、Ｖｏｌ．２４、２５５〜２６９
頁、１９８３年に記載されており、これを本明細書の一
部を構成するものしてここに援用する。実際には、ＨＭ
ＲＦに基づくフィルタ処理とσ−フィルタに基づくフィ
ルタ処理を用いて再構成されたコントーン画像を連続的
に処理することによって最良の結果が得られる。
どは、縁の回りを除いてツリーに基づくクリッピングプ
ロセスにおいて抑制されるか濾去される。各高周波サブ
バンド画像は、主にその方位に沿った情報を含んでいる
ため、上記のプロセスは、画像データから残存するハー
フトーン雑音を更に濾去するために使用できる。垂直ま
たは水平に方向づけられたサブバンド画像に対してはそ
れらの方位に沿って一次元のローパスフィルタが適用さ
れる。傾いて方向づけられたサブバンド画像に対しては
Ｘ形状のローパスフィルタが適用される。すなわち、傾
斜したサブバンド画像は、より厳しくローパスフィルタ
処理される。なぜなら、それらは一般に、含んでいる信
号情報が少なく、より多くのハーフトーン雑音を含んで
いるからである。この第２レベルのフィルタ処理は、ロ
ーパスＬＬサブバンド画像と同じ粗さの解画像度のレベ
ルのサブバンド画像を含め、全ての高周波サブバンド画
像に適用される。方向性フィルタ処理は、例えば、Ｌポ
イントランニング平均化であることが好ましい。 一般
に、２つのサブバンド画像フィルタリング工程を経た初
期の再構成されたコントーン画像は、信号情報の消失に
起因して幾分曖昧になるかまたはむらが生じている。特
に、高周波信号情報の消失により、“リンギング（ring
ing ）”アーティファクト（artifact）が生じる。した
がって、再構築された連続色調画像は、マルコフのラン
ダムフィールドによってモデル化した空間平滑度拘束を
有する非線型フィルタ処理を用いて強調処理することが
好ましい。マルコフのランダムフィールドは、フーバー
ミニマックス関数と呼ばれる特定のギブズポテンシャル
関数を使用して組み込まれる。この関数が図１８に示さ
れている。この関数の二次セグメントは閾値Ｔ以下での
不連続部分の平滑化を容易にし、直線セグメントは鋭い
縁部を保存する。この様な画像強調技術は、アーティフ
ァクトを大幅に減少させ、元の画像の細部を保存する。
σ−フィルタに基づくフィルタ処理などの他の縁保存フ
ィルタ処理方法も使用できる。σ−フィルタに基づくフ
ィルタ処理は、Ｊ．リー（Lee ）による「デジタルイメ
ージスムージングおよびσ−フィルタ（Digital Image
Smoothing and the σ-Filter ）、Ｃｏｍｐｕｔｅｒ
Ｉｍａｇｅ Ｇｒａｐｈｉｃｓ ａｎｄ Ｉｍａｇｅ
Ｐｒｏｃｅｓｓｉｎｇ、Ｖｏｌ．２４、２５５〜２６９
頁、１９８３年に記載されており、これを本明細書の一
部を構成するものしてここに援用する。実際には、ＨＭ
ＲＦに基づくフィルタ処理とσ−フィルタに基づくフィ
ルタ処理を用いて再構成されたコントーン画像を連続的
に処理することによって最良の結果が得られる。
【００４４】したがって、図４に示すように、Ｎステー
ジのハーフトーンデコンポーザ１４０は、第１ステージ
（ハーフトーンからサブバンドへの）デコンポーザ１４
１、第２ステージ（サブバンドからサブバンドへの）デ
コンポーザ１４２など第Ｎステージ（サブバンドからサ
ブバンドへの）デコンポーザまでを含んでいる。入力ハ
ーフトーン画像は、データバス１９５を介してＮステー
ジのハーフトーンデコンポーザ１４０に入力され、第１
ステージデコンポーザ１４１に到る。第１ステージデコ
ンポーザ１４１は、ローパスＬＬ1 サブバンド画像と３
個の高周波ＬＨ1 、ＨＬ1 、およびＨＨ1 サブバンド画
像を発生し、これらは分解されたサブバンド画像として
データバス１９５に出力される。
ジのハーフトーンデコンポーザ１４０は、第１ステージ
（ハーフトーンからサブバンドへの）デコンポーザ１４
１、第２ステージ（サブバンドからサブバンドへの）デ
コンポーザ１４２など第Ｎステージ（サブバンドからサ
ブバンドへの）デコンポーザまでを含んでいる。入力ハ
ーフトーン画像は、データバス１９５を介してＮステー
ジのハーフトーンデコンポーザ１４０に入力され、第１
ステージデコンポーザ１４１に到る。第１ステージデコ
ンポーザ１４１は、ローパスＬＬ1 サブバンド画像と３
個の高周波ＬＨ1 、ＨＬ1 、およびＨＨ1 サブバンド画
像を発生し、これらは分解されたサブバンド画像として
データバス１９５に出力される。
【００４５】低周波ＬＬ1 サブバンド画像は第２ステー
ジデコンポーザ１４２に入力される。この第２ステージ
デコンポーザ１４２は、低周波ＬＬ1 サブバンド画像を
低周波ＬＬ2 サブバンド画像と３個の高周波ＬＨ2 、Ｈ
Ｌ2 、およびＨＨ2 サブバンド画像に分解する。３個の
高周波サブバンド画像もデータバス１９５に出力され
る。分解プロセスが続けられ、各後続のデコンポーザス
テージは前のデコンポーザステージからのＬＬサブバン
ド画像を分解する。これが第Ｎステージデコンポーザ１
４Ｎまで行われる。低周波ＬＬN サブバンド画像は更に
分解されるのではなく、第Ｎステージデコンポーザ１４
Ｎによって３個の高周波ＬＨN 、ＨＬN 、およびＨＨN
サブバンド画像と一緒にデータバス１９５に出力され
る。
ジデコンポーザ１４２に入力される。この第２ステージ
デコンポーザ１４２は、低周波ＬＬ1 サブバンド画像を
低周波ＬＬ2 サブバンド画像と３個の高周波ＬＨ2 、Ｈ
Ｌ2 、およびＨＨ2 サブバンド画像に分解する。３個の
高周波サブバンド画像もデータバス１９５に出力され
る。分解プロセスが続けられ、各後続のデコンポーザス
テージは前のデコンポーザステージからのＬＬサブバン
ド画像を分解する。これが第Ｎステージデコンポーザ１
４Ｎまで行われる。低周波ＬＬN サブバンド画像は更に
分解されるのではなく、第Ｎステージデコンポーザ１４
Ｎによって３個の高周波ＬＨN 、ＨＬN 、およびＨＨN
サブバンド画像と一緒にデータバス１９５に出力され
る。
【００４６】同様に、図５に示すように、サブバンド画
像を連続色調画像に再構成するＮステージリコンポーザ
１６０は、第１ステージ（サブバンドからサブバンドへ
の）リコンポーザ１６１、第２ステージ（サブバンドか
らサブバンドへの）リコンポーザ１６２など第Ｎステー
ジ（サブバンドから連続色調画像への）リコンポーザ１
６Ｎまでを含んでいる。図５に示すように、サブバンド
画像は、データバス１９５を介してサブバンド画像を連
続色調画像に再構成するＮステージリコンポーザ１６０
に入力される。すなわち、高周波および低周波のＬ
ＬN 、ＬＨN 、ＨＬN 、およびＨＨN サブバンド画像は
第１ステージリコンポーザ１６１へ入力される。第１ス
テージリコンポーザ１６１は、再構成された低周波ＬＬ
(N-1) サブバンド画像を第２ステージリコンポーザ１６
２出力する。第２ステージリコンポーザ１６２は、高周
波ＬＨ(N-1) 、ＨＬ(N-1) 、ＨＨ(N-1) サブバンド画像
も入力する。第２ステージリコンポーザ１６２は、これ
らの高周波および低周波のサブバンド画像を再構成し、
再構成された低周波ＬＬN-2 サブバンド画像を出力す
る。
像を連続色調画像に再構成するＮステージリコンポーザ
１６０は、第１ステージ（サブバンドからサブバンドへ
の）リコンポーザ１６１、第２ステージ（サブバンドか
らサブバンドへの）リコンポーザ１６２など第Ｎステー
ジ（サブバンドから連続色調画像への）リコンポーザ１
６Ｎまでを含んでいる。図５に示すように、サブバンド
画像は、データバス１９５を介してサブバンド画像を連
続色調画像に再構成するＮステージリコンポーザ１６０
に入力される。すなわち、高周波および低周波のＬ
ＬN 、ＬＨN 、ＨＬN 、およびＨＨN サブバンド画像は
第１ステージリコンポーザ１６１へ入力される。第１ス
テージリコンポーザ１６１は、再構成された低周波ＬＬ
(N-1) サブバンド画像を第２ステージリコンポーザ１６
２出力する。第２ステージリコンポーザ１６２は、高周
波ＬＨ(N-1) 、ＨＬ(N-1) 、ＨＨ(N-1) サブバンド画像
も入力する。第２ステージリコンポーザ１６２は、これ
らの高周波および低周波のサブバンド画像を再構成し、
再構成された低周波ＬＬN-2 サブバンド画像を出力す
る。
【００４７】この再構成プロセスは、第Ｎステージリコ
ンポーザ１６Ｎまで、一連のリコンポーザステージにお
いて継続される。第Ｎステージリコンポーザ１６Ｎは、
第Ｎ−１ステージ（サブバンドからサブバンドへの）リ
コンポーザから出力される再構成された低周波のＬＬN
サブバンド画像を入力し、またデータバス１９５から高
周波のＬＨ1 、ＨＬ1 、およびＨＨ1 サブバンド画像を
入力し、初期連続色調画像を出力する。その後、サブバ
ンド画像を連続色調画像に再構成するＮステージリコン
ポーザ１６０は初期連続色調画像をデータバス１９５に
出力する。
ンポーザ１６Ｎまで、一連のリコンポーザステージにお
いて継続される。第Ｎステージリコンポーザ１６Ｎは、
第Ｎ−１ステージ（サブバンドからサブバンドへの）リ
コンポーザから出力される再構成された低周波のＬＬN
サブバンド画像を入力し、またデータバス１９５から高
周波のＬＨ1 、ＨＬ1 、およびＨＨ1 サブバンド画像を
入力し、初期連続色調画像を出力する。その後、サブバ
ンド画像を連続色調画像に再構成するＮステージリコン
ポーザ１６０は初期連続色調画像をデータバス１９５に
出力する。
【００４８】図６は第１ステージデコンポーザ１４１の
構成をより詳細に示す。すなわち、データバス１９５上
に入力された入力ハーフトーン画像はローパスフィルタ
１４１１とハイパスフィルタ１４１２の各々に入力され
る。ローパスフィルタ１４１１の出力はローパスフィル
タ１４１３とハイパスフィルタ１４１４の各々に入力さ
れる。同様に、ハイパスフィルタ１４１２の出力はロー
パスフィルタ１４１５とハイパスフィルタ１４１６の各
々に入力される。ローパスフィルタ１４１３の出力がロ
ーパスＬＬ1 サブバンド画像であり、ハイパスフィルタ
１４１４の出力が高周波ＬＨ1 サブバンド画像である。
同様に、ローパスフィルタ１４１５の出力が高周波ＨＬ
1 サブバンド画像であり、ハイパスフィルタ１４１６の
出力が高周波ＨＨ1 サブバンド画像である。第１ステー
ジデコンポーザ１４１と同じ構成が第２ステージ〜第Ｎ
スージデコンポーザ１４２−１４Ｎに使用できる。
構成をより詳細に示す。すなわち、データバス１９５上
に入力された入力ハーフトーン画像はローパスフィルタ
１４１１とハイパスフィルタ１４１２の各々に入力され
る。ローパスフィルタ１４１１の出力はローパスフィル
タ１４１３とハイパスフィルタ１４１４の各々に入力さ
れる。同様に、ハイパスフィルタ１４１２の出力はロー
パスフィルタ１４１５とハイパスフィルタ１４１６の各
々に入力される。ローパスフィルタ１４１３の出力がロ
ーパスＬＬ1 サブバンド画像であり、ハイパスフィルタ
１４１４の出力が高周波ＬＨ1 サブバンド画像である。
同様に、ローパスフィルタ１４１５の出力が高周波ＨＬ
1 サブバンド画像であり、ハイパスフィルタ１４１６の
出力が高周波ＨＨ1 サブバンド画像である。第１ステー
ジデコンポーザ１４１と同じ構成が第２ステージ〜第Ｎ
スージデコンポーザ１４２−１４Ｎに使用できる。
【００４９】同様に、図７は第Ｎステージリコンポーザ
１６Ｎの好適な構成を示す。図７に示すように、低周波
数ＬＬ1 サブバンド画像がローパスフィルタ１６１１に
入力される。高周波ＬＨ1 サブバンド画像がハイパスフ
ィルタ１６１２に入力される。ローパスフィルタ１６１
１とハイパスフィルタ１６１２の出力は合成され、ロー
パスフィルタ１６１５に出力される。同様に、高周波Ｈ
Ｌ1 サブバンド画像がローパスフィルタ１６１３に入力
され、高周波ＨＨ1 サブバンド画像がハイパスフィルタ
１６１４に入力される。ローパスフィルタ１６１３とハ
イパスフィルタ１６１４の出力は合成され、ハイパスフ
ィルタ１６１６に出力される。ローパスフィルタ１６１
５とハイパスフィルタ１６１６の出力は合成され、連続
色調画像として出力される。サブバンド画像を連続色調
画像に再構成するＮステージリコンポーザ１６０の第１
スージから第Ｎ−１ステージの各々は同一の構成である
ことが好ましい。これらのステージの各々では、ローパ
スフィルタ１６１５とハイパスフィルタ１６１６の出力
は、連続色調画像であり、低周波ＬＬサブバンド画像で
ある。
１６Ｎの好適な構成を示す。図７に示すように、低周波
数ＬＬ1 サブバンド画像がローパスフィルタ１６１１に
入力される。高周波ＬＨ1 サブバンド画像がハイパスフ
ィルタ１６１２に入力される。ローパスフィルタ１６１
１とハイパスフィルタ１６１２の出力は合成され、ロー
パスフィルタ１６１５に出力される。同様に、高周波Ｈ
Ｌ1 サブバンド画像がローパスフィルタ１６１３に入力
され、高周波ＨＨ1 サブバンド画像がハイパスフィルタ
１６１４に入力される。ローパスフィルタ１６１３とハ
イパスフィルタ１６１４の出力は合成され、ハイパスフ
ィルタ１６１６に出力される。ローパスフィルタ１６１
５とハイパスフィルタ１６１６の出力は合成され、連続
色調画像として出力される。サブバンド画像を連続色調
画像に再構成するＮステージリコンポーザ１６０の第１
スージから第Ｎ−１ステージの各々は同一の構成である
ことが好ましい。これらのステージの各々では、ローパ
スフィルタ１６１５とハイパスフィルタ１６１６の出力
は、連続色調画像であり、低周波ＬＬサブバンド画像で
ある。
【００５０】図８はサブバンド雑音および／またはパタ
ーンフィルタシステム１５０の好適な構成をより詳細に
示す。すなわち、Ｎステージハーフトーンデコンポーザ
１４０によって発生されたサブバンド画像はデータバス
１９５を介して非縁領域雑音減衰システム１５１に入力
される。上記したように、この非縁領域雑音減衰システ
ム１５１は好ましくはツリークリッパである。非縁領域
雑音減衰システム１５１の出力は、方向性フィルタシス
テム１５２に出力される。
ーンフィルタシステム１５０の好適な構成をより詳細に
示す。すなわち、Ｎステージハーフトーンデコンポーザ
１４０によって発生されたサブバンド画像はデータバス
１９５を介して非縁領域雑音減衰システム１５１に入力
される。上記したように、この非縁領域雑音減衰システ
ム１５１は好ましくはツリークリッパである。非縁領域
雑音減衰システム１５１の出力は、方向性フィルタシス
テム１５２に出力される。
【００５１】図８に示すように、方向性フィルタシステ
ム１５２は、サブバンド方位決定器１５３を含んでいる
ことが好ましい。このサブバンド方位決定器１５３は非
縁領域雑音減衰システム１５１からの出力を入力し、Ｎ
ステージハーフトーンデコンポーザ１４０によって発生
されたサブバンド画像の各々について方位を決定する。
水平サブバンド画像は水平サブバンドフィルタ１５４に
出力され、垂直サブバンド画像は垂直サブバンドフィル
タ１５５に出力され、傾斜サブバンド画像は傾斜サブバ
ンドフィルタ１５６に出力される。好ましくは、水平サ
ブバンドフィルタ１５４は垂直ローパスフィルタを水平
サブバンド画像に適用し、垂直サブバンドフィルタ１５
５は水平ローパスフィルタを垂直サブバンド画像に適用
する。傾斜サブバンドフィルタ１５６によりＸ形ローバ
スフィルタが傾斜サブバンド画像に対して適用されるこ
とが好ましい。水平、垂直および傾斜サブバンド画像フ
ィルタ１５４〜１５６の出力がサブバンド雑音および／
またはパターンフィルタシステム１５０によってデータ
バス１９５上に出力され、メモリ１２０またはサブバン
ド画像を連続色調画像に再構成するＮステージリコンポ
ーザ１６０に到る。
ム１５２は、サブバンド方位決定器１５３を含んでいる
ことが好ましい。このサブバンド方位決定器１５３は非
縁領域雑音減衰システム１５１からの出力を入力し、Ｎ
ステージハーフトーンデコンポーザ１４０によって発生
されたサブバンド画像の各々について方位を決定する。
水平サブバンド画像は水平サブバンドフィルタ１５４に
出力され、垂直サブバンド画像は垂直サブバンドフィル
タ１５５に出力され、傾斜サブバンド画像は傾斜サブバ
ンドフィルタ１５６に出力される。好ましくは、水平サ
ブバンドフィルタ１５４は垂直ローパスフィルタを水平
サブバンド画像に適用し、垂直サブバンドフィルタ１５
５は水平ローパスフィルタを垂直サブバンド画像に適用
する。傾斜サブバンドフィルタ１５６によりＸ形ローバ
スフィルタが傾斜サブバンド画像に対して適用されるこ
とが好ましい。水平、垂直および傾斜サブバンド画像フ
ィルタ１５４〜１５６の出力がサブバンド雑音および／
またはパターンフィルタシステム１５０によってデータ
バス１９５上に出力され、メモリ１２０またはサブバン
ド画像を連続色調画像に再構成するＮステージリコンポ
ーザ１６０に到る。
【００５２】図９は、元の画像“Ｏ”とサブバンド画像
“Ｓ”の各々との間の全体的な関係を示す。すなわち、
図９に示すように、元の画像Ｏはｍ＋１個の第１レベル
のサブバンド画像Ｓ（１，０）〜Ｓ（１，Ｍ）に分解さ
れる。そして、各レベルの第０番目のサブバンド画像、
すなわちＳ（ｎ，０）サブバンド画像が最大でｍ＋１個
の第２レベルのサブバンド画像Ｓ（２，０）〜Ｓ（２，
Ｍ）に分解される。この処理は、サブバンド画像Ｓ
（Ｎ，０）〜Ｓ（Ｎ，Ｍ）を含む第ｎ番目のレベルが形
成されるまで継続される。
“Ｓ”の各々との間の全体的な関係を示す。すなわち、
図９に示すように、元の画像Ｏはｍ＋１個の第１レベル
のサブバンド画像Ｓ（１，０）〜Ｓ（１，Ｍ）に分解さ
れる。そして、各レベルの第０番目のサブバンド画像、
すなわちＳ（ｎ，０）サブバンド画像が最大でｍ＋１個
の第２レベルのサブバンド画像Ｓ（２，０）〜Ｓ（２，
Ｍ）に分解される。この処理は、サブバンド画像Ｓ
（Ｎ，０）〜Ｓ（Ｎ，Ｍ）を含む第ｎ番目のレベルが形
成されるまで継続される。
【００５３】必要とされる処理を適用するためにサブバ
ンド画像が処理された後、連続色調画像がプロセスを逆
に行うことによって作られる。本発明のウエーブレット
に基づく非スクリーニング装置および方法では、各レベ
ルにおいて４個のサブバンド画像を発生するが、図９に
示された方法はそれに限定されていない。なぜなら、本
発明の汎用ハーフトーン画像非スクリーニング装置およ
び方法は、離散ウエーブレット変換に限定されないから
である。余弦変換（離散余弦変換を含む）、重ね変換
（lapped transforms ）、フーリエ変換、サブサンプリ
ングのないウエーブレット変換などを本発明の汎用ハー
フトーン画像非スクリーニング装置および方法に使用で
きる。したがって、本発明の最も一般的な方法が図１１
に示されている。ステップＳ１００で起動した後、図１
１〜図１５のフローチャートを実施する制御システムは
ステップＳ１１０においてハーフトーン画像を入力す
る。そして、ステップＳ１２０において、ハーフトーン
画像がサブバンド画像に分解される。次に、ステップＳ
１３０において、雑音および／またはパターンがサブバ
ンド画像から選択的に除去される。そして、ステップＳ
１４０において、サブバンド画像が初期連続色調画像に
再構成される。
ンド画像が処理された後、連続色調画像がプロセスを逆
に行うことによって作られる。本発明のウエーブレット
に基づく非スクリーニング装置および方法では、各レベ
ルにおいて４個のサブバンド画像を発生するが、図９に
示された方法はそれに限定されていない。なぜなら、本
発明の汎用ハーフトーン画像非スクリーニング装置およ
び方法は、離散ウエーブレット変換に限定されないから
である。余弦変換（離散余弦変換を含む）、重ね変換
（lapped transforms ）、フーリエ変換、サブサンプリ
ングのないウエーブレット変換などを本発明の汎用ハー
フトーン画像非スクリーニング装置および方法に使用で
きる。したがって、本発明の最も一般的な方法が図１１
に示されている。ステップＳ１００で起動した後、図１
１〜図１５のフローチャートを実施する制御システムは
ステップＳ１１０においてハーフトーン画像を入力す
る。そして、ステップＳ１２０において、ハーフトーン
画像がサブバンド画像に分解される。次に、ステップＳ
１３０において、雑音および／またはパターンがサブバ
ンド画像から選択的に除去される。そして、ステップＳ
１４０において、サブバンド画像が初期連続色調画像に
再構成される。
【００５４】次のステップＳ１５０において、制御シス
テムは、初期連続色調画像のフィルタ後処理が必要であ
るか、好ましいか、または要求されているか否かを判定
する。フィルタ後処理が行われる場合、制御はステップ
Ｓ１６０へ移行し、ここで初期連続色調画像がフィルタ
処理されて最終連続色調画像が作られる。その後、制御
はステップＳ１７０に移行する。もし、ステップＳ１５
０においてフィルタ後処理を行わないと判断した場合、
制御はステップＳ１７０へ直接ジャンプする。ステップ
Ｓ１７０でプロセスが終了する。
テムは、初期連続色調画像のフィルタ後処理が必要であ
るか、好ましいか、または要求されているか否かを判定
する。フィルタ後処理が行われる場合、制御はステップ
Ｓ１６０へ移行し、ここで初期連続色調画像がフィルタ
処理されて最終連続色調画像が作られる。その後、制御
はステップＳ１７０に移行する。もし、ステップＳ１５
０においてフィルタ後処理を行わないと判断した場合、
制御はステップＳ１７０へ直接ジャンプする。ステップ
Ｓ１７０でプロセスが終了する。
【００５５】図１２はハーフトーン画像入力ステップＳ
１１０をより詳細に示す。すなわち、ステップＳ１１２
において、制御システムは、走査したハーフトーンをス
キャナから入力しなければならないかを判定する。もし
そうである場合、制御はステップＳ１１４に移行し、こ
こでハーフトーン画像が走査され入力される。そうでな
い場合、制御はステップＳ１１６に移行し、ここでハー
フトーン画像は、ネットワークまたはホストコンピュー
タからモデムを介してメモリなどから電子的ハーフトー
ン画像として入力される。制御はステップＳ１１４また
はステップＳ１１６の何れかよりステップＳ１１８に移
行し、ステップＳ１２０に戻る。
１１０をより詳細に示す。すなわち、ステップＳ１１２
において、制御システムは、走査したハーフトーンをス
キャナから入力しなければならないかを判定する。もし
そうである場合、制御はステップＳ１１４に移行し、こ
こでハーフトーン画像が走査され入力される。そうでな
い場合、制御はステップＳ１１６に移行し、ここでハー
フトーン画像は、ネットワークまたはホストコンピュー
タからモデムを介してメモリなどから電子的ハーフトー
ン画像として入力される。制御はステップＳ１１４また
はステップＳ１１６の何れかよりステップＳ１１８に移
行し、ステップＳ１２０に戻る。
【００５６】図１３はハーフトーン画像分解ステップＳ
１２０をより詳細に示す。すなわち、図１３は、元のハ
ーフトーン画像を分解するために離散ウエーブレット変
換を使用する場合の、ステップＳ１２０の好ましい方法
を示す。ステップＳ１２０のための図１３に示すステッ
プは、ハーフトーン画像を分割するのに使用した特定の
方法に大きく依存する。
１２０をより詳細に示す。すなわち、図１３は、元のハ
ーフトーン画像を分解するために離散ウエーブレット変
換を使用する場合の、ステップＳ１２０の好ましい方法
を示す。ステップＳ１２０のための図１３に示すステッ
プは、ハーフトーン画像を分割するのに使用した特定の
方法に大きく依存する。
【００５７】ステップＳ１２１において、カウンタＮが
１にセットされ、元のハーフトーン画像はＬＬ1 、ＬＨ
1 、ＨＬ1 およびＨＨ1 サブバンド画像に分解される。
そして、ステップＳ１２２でこれらのサブバンド画像が
メモリに記憶される。次のステップＳ１２３において、
制御システムは現在の低周波ＬＬN サブバンド画像（但
し、ここではＮ＝１）を第Ｎ＋１レベルのサブバンド画
像に更に分解すべきか否かを決定する。もし、分解しな
い場合、制御はステップＳ１２７へ直接ジャンプする。
そうでない場合、制御はステップＳ１２４へ移行する。
１にセットされ、元のハーフトーン画像はＬＬ1 、ＬＨ
1 、ＨＬ1 およびＨＨ1 サブバンド画像に分解される。
そして、ステップＳ１２２でこれらのサブバンド画像が
メモリに記憶される。次のステップＳ１２３において、
制御システムは現在の低周波ＬＬN サブバンド画像（但
し、ここではＮ＝１）を第Ｎ＋１レベルのサブバンド画
像に更に分解すべきか否かを決定する。もし、分解しな
い場合、制御はステップＳ１２７へ直接ジャンプする。
そうでない場合、制御はステップＳ１２４へ移行する。
【００５８】ステップＳ１２４では、低周波ＬＬN サブ
バンド画像がメモリから入力され、Ｎの値が１だけ増加
される。その後、ステップＳ１２５において、低周波Ｌ
Ｌ(N-1) サブバンド画像となっている入力低周波サブバ
ンド画像が次のレベルのＬＬN 、ＬＨN 、ＨＬN および
ＨＨN サブバンド画像に分解される。その後、これらの
サブバンド画像はステップＳ１２６においてメモリに記
憶され、制御はステップＳ１２３へ戻る。制御システム
は、全てのサブバンド画像レベルが形成されるまでステ
ップＳ１２３〜Ｓ１２６を通るループを継続する。制御
はその後ステップＳ１２７へ移行する。
バンド画像がメモリから入力され、Ｎの値が１だけ増加
される。その後、ステップＳ１２５において、低周波Ｌ
Ｌ(N-1) サブバンド画像となっている入力低周波サブバ
ンド画像が次のレベルのＬＬN 、ＬＨN 、ＨＬN および
ＨＨN サブバンド画像に分解される。その後、これらの
サブバンド画像はステップＳ１２６においてメモリに記
憶され、制御はステップＳ１２３へ戻る。制御システム
は、全てのサブバンド画像レベルが形成されるまでステ
ップＳ１２３〜Ｓ１２６を通るループを継続する。制御
はその後ステップＳ１２７へ移行する。
【００５９】ステップＳ１２７では、記憶された第１〜
第Ｎレベルのサブバンド画像が、雑音および／またはパ
ターンの選択的除去のために出力される。その後、制御
はステップＳ１２８に移行し、ステップＳ１３０に戻
る。サブバンド画像は必ずしもメモリに記憶しなくても
よい。この場合、ステップＳ１２２とＳ１２６において
サブバンド画像をメモリに記憶する代わりに、ステップ
Ｓ１２１とＳ１２５で形成された低周波ＬＬN サブバン
ド画像を、ステップＳ１２３において追加のサブバンド
画像を形成すべきであると判断された場合にステップＳ
１２４に直接入力する。
第Ｎレベルのサブバンド画像が、雑音および／またはパ
ターンの選択的除去のために出力される。その後、制御
はステップＳ１２８に移行し、ステップＳ１３０に戻
る。サブバンド画像は必ずしもメモリに記憶しなくても
よい。この場合、ステップＳ１２２とＳ１２６において
サブバンド画像をメモリに記憶する代わりに、ステップ
Ｓ１２１とＳ１２５で形成された低周波ＬＬN サブバン
ド画像を、ステップＳ１２３において追加のサブバンド
画像を形成すべきであると判断された場合にステップＳ
１２４に直接入力する。
【００６０】図１４は、選択的に雑音および／またはパ
ターンを除去するステップＳ１３０をより詳細に示す。
すなわち、ステップＳ１３２において、サブバンド画像
の縁でない領域内の雑音が減衰される。そして、ステッ
プＳ１３４において、方向性フィルタを使用してサブバ
ンド画像がフィルタ処理される。そして、ステップＳ１
３６において、制御はステップＳ１４０に戻る。
ターンを除去するステップＳ１３０をより詳細に示す。
すなわち、ステップＳ１３２において、サブバンド画像
の縁でない領域内の雑音が減衰される。そして、ステッ
プＳ１３４において、方向性フィルタを使用してサブバ
ンド画像がフィルタ処理される。そして、ステップＳ１
３６において、制御はステップＳ１４０に戻る。
【００６１】図１５は、連続色調画像を再構成するステ
ップＳ１４０をより詳細に示す。すなわち、ステップＳ
１４１において、第ＮレベルのＬＬN 、ＬＨN 、ＨＬN
およびＨＨN サブバンド画像をメモリから入力する。そ
して、ステップＳ１４２において、新しい第Ｎ−１レベ
ルの低周波ＬＬ(N-1) サブバンド画像を形成するために
第Ｎレベルのサブバンド画像が再構成される。そして、
ステップＳ１４３において、ステップＳ１２４で増加さ
れたＮの値が１だけ減少される。その後、制御はステッ
プＳ１４４に移行する。
ップＳ１４０をより詳細に示す。すなわち、ステップＳ
１４１において、第ＮレベルのＬＬN 、ＬＨN 、ＨＬN
およびＨＨN サブバンド画像をメモリから入力する。そ
して、ステップＳ１４２において、新しい第Ｎ−１レベ
ルの低周波ＬＬ(N-1) サブバンド画像を形成するために
第Ｎレベルのサブバンド画像が再構成される。そして、
ステップＳ１４３において、ステップＳ１２４で増加さ
れたＮの値が１だけ減少される。その後、制御はステッ
プＳ１４４に移行する。
【００６２】ステップＳ１４４では、制御システムはＮ
が零に等しいか否かを判断する。もし零でない場合、制
御はステップＳ１４５に移行し、ここで新しいＮの値に
基づき、第Ｎレベルの高周波ＬＨN 、ＨＬN およびＨＨ
N サブバンド画像が入力され、再構成された低周波ＬＬ
N サブバンド画像と関連づけられる。その後、制御はス
テップＳ１４２戻り、そこで次のより高いレベルの低周
波ＬＬサブバンド画像を形成するためにサブバンド画像
が再び構成される。
が零に等しいか否かを判断する。もし零でない場合、制
御はステップＳ１４５に移行し、ここで新しいＮの値に
基づき、第Ｎレベルの高周波ＬＨN 、ＨＬN およびＨＨ
N サブバンド画像が入力され、再構成された低周波ＬＬ
N サブバンド画像と関連づけられる。その後、制御はス
テップＳ１４２戻り、そこで次のより高いレベルの低周
波ＬＬサブバンド画像を形成するためにサブバンド画像
が再び構成される。
【００６３】ステップＳ１４２〜１４５を通るこのルー
プは、ステップＳ１４４においてＮ＝０と判断されるま
で継続される。この時、制御はステップＳ１４６に移行
し、ここで“ＬＬ0 ”画像、すなわち再構成されたコン
トーン画像がメモリに出力される。その後、ステップＳ
１４７において制御はステップＳ１５０に戻る。
プは、ステップＳ１４４においてＮ＝０と判断されるま
で継続される。この時、制御はステップＳ１４６に移行
し、ここで“ＬＬ0 ”画像、すなわち再構成されたコン
トーン画像がメモリに出力される。その後、ステップＳ
１４７において制御はステップＳ１５０に戻る。
【００６４】図１６は、方向性フィルタを使用してサブ
バンド画像をフィルタ処理するステップＳ１３４のため
の一つの実施の形態を示す。すなわち、ステップＳ２０
０において、分解プロセスによって出力または発生され
たサブバンド画像の最初の一つが選択される。その後、
ステップＳ２１０において、制御システムは選択された
サブバンド画像が水平サブバンド画像であるか否かを判
定する。もし水平サブバンド画像でない場合、制御はス
テップＳ２３０へ移行する。そうでなければ、制御はス
テップＳ２２０へ移行する。ステップＳ２２０では、垂
直ローパスフィルタを使用して水平サブバンド画像がフ
ィルタ処理される。その後、制御はステップＳ２６０へ
移行する。ステップＳ２３０では、制御システムは選択
されたサブバンド画像が垂直サブバンド画像であるか否
かを判定する。もし垂直サブバンド画像である場合、制
御はステップＳ２４０へ移行し、ここで水平ローパスフ
ィルタを使用して垂直サブバンド画像がフィルタ処理さ
れる。その後、制御はステップＳ２６０へ移行する。一
方、ステップＳ２３０において、選択されたサブバンド
画像が垂直サブバンド画像でないと判断された場合に
は、それは傾斜サブバンド画像である。したがって、制
御はステップＳ２５０へ移行し、ここでＸ形ローパスフ
ィルタを使用して傾斜サブバンド画像がフィルタ処理さ
れる。その後、制御はステップＳ２６０へ移行する。
バンド画像をフィルタ処理するステップＳ１３４のため
の一つの実施の形態を示す。すなわち、ステップＳ２０
０において、分解プロセスによって出力または発生され
たサブバンド画像の最初の一つが選択される。その後、
ステップＳ２１０において、制御システムは選択された
サブバンド画像が水平サブバンド画像であるか否かを判
定する。もし水平サブバンド画像でない場合、制御はス
テップＳ２３０へ移行する。そうでなければ、制御はス
テップＳ２２０へ移行する。ステップＳ２２０では、垂
直ローパスフィルタを使用して水平サブバンド画像がフ
ィルタ処理される。その後、制御はステップＳ２６０へ
移行する。ステップＳ２３０では、制御システムは選択
されたサブバンド画像が垂直サブバンド画像であるか否
かを判定する。もし垂直サブバンド画像である場合、制
御はステップＳ２４０へ移行し、ここで水平ローパスフ
ィルタを使用して垂直サブバンド画像がフィルタ処理さ
れる。その後、制御はステップＳ２６０へ移行する。一
方、ステップＳ２３０において、選択されたサブバンド
画像が垂直サブバンド画像でないと判断された場合に
は、それは傾斜サブバンド画像である。したがって、制
御はステップＳ２５０へ移行し、ここでＸ形ローパスフ
ィルタを使用して傾斜サブバンド画像がフィルタ処理さ
れる。その後、制御はステップＳ２６０へ移行する。
【００６５】ステップＳ２６０では、フィルタ処理され
たサブバンド画像がメモリに記憶される。その後、制御
はステップＳ２７０へ移行する。ステップＳ２７０で
は、制御システムは、選択されたサブバンド画像がフィ
ルタ処理すべき最後のサブバンド画像であるか否かを判
定する。すなわち、制御システムは、全てのサブバンド
画像がフィルタ処理されたか否かを判定する。もし、全
てのサブバンド画像がフィルタ処理されていない場合、
制御はステップＳ２８０へ移行し、サブバンド画像の他
の一つが選択される。その後、制御はステップＳ２１０
へ戻る。一方、全てのサブバンド画像がフィルタ処理さ
れている場合、制御はステップＳ２９０へ移行し、ステ
ップＳ１３６へ戻る。
たサブバンド画像がメモリに記憶される。その後、制御
はステップＳ２７０へ移行する。ステップＳ２７０で
は、制御システムは、選択されたサブバンド画像がフィ
ルタ処理すべき最後のサブバンド画像であるか否かを判
定する。すなわち、制御システムは、全てのサブバンド
画像がフィルタ処理されたか否かを判定する。もし、全
てのサブバンド画像がフィルタ処理されていない場合、
制御はステップＳ２８０へ移行し、サブバンド画像の他
の一つが選択される。その後、制御はステップＳ２１０
へ戻る。一方、全てのサブバンド画像がフィルタ処理さ
れている場合、制御はステップＳ２９０へ移行し、ステ
ップＳ１３６へ戻る。
【００６６】図１７は、各サブバンド画像のそれらの親
の画像に対する関係と、親の係数と子供の係数について
の親子関係を示す。図１９は、図１１と図１３のハーフ
トーン画像分解ステップＳ１２０が、図２の集合ドット
ハーフトーン画像に適用された場合に得られるサブバン
ド画像を示す。図２０は図１９に示すサブバンド画像に
非縁領域雑音減衰ステップＳ１３２を適用した後のサブ
バンド画像を示す。すなわち、図２０に示す雑音の減衰
されたサブバンド画像では、ステップＳ１３２がツリー
クリッピングによって実施される。
の画像に対する関係と、親の係数と子供の係数について
の親子関係を示す。図１９は、図１１と図１３のハーフ
トーン画像分解ステップＳ１２０が、図２の集合ドット
ハーフトーン画像に適用された場合に得られるサブバン
ド画像を示す。図２０は図１９に示すサブバンド画像に
非縁領域雑音減衰ステップＳ１３２を適用した後のサブ
バンド画像を示す。すなわち、図２０に示す雑音の減衰
されたサブバンド画像では、ステップＳ１３２がツリー
クリッピングによって実施される。
【００６７】次に、図２１は、図１４と図１６に示す方
向性フィルタを使用したサブバンドフィルタ処理ステッ
プＳ１３４を図２０の雑音を減衰されたサブバンド画像
に適用した後のサブバンド画像を示す。図２２は、図１
１と図１５に示すサブバンド画像を連続色調画像に再構
成するステップＳ１４０に基づき、図２１に示すフィル
タ処理された画像から再構成された初期再構成コントー
ン画像を示す。
向性フィルタを使用したサブバンドフィルタ処理ステッ
プＳ１３４を図２０の雑音を減衰されたサブバンド画像
に適用した後のサブバンド画像を示す。図２２は、図１
１と図１５に示すサブバンド画像を連続色調画像に再構
成するステップＳ１４０に基づき、図２１に示すフィル
タ処理された画像から再構成された初期再構成コントー
ン画像を示す。
【００６８】図２３は、図１１の連続色調画像フィルタ
処理ステップＳ１６０を図２２に示す初期コントーン画
像に適用することによって形成された最終コントーン画
像を示す。 図２４は、図１１に概要が示されたプロセ
スを図１に示す誤差拡散されたハーフトーン画像に適用
して得られた同様の最終コントーン画像を示す。
処理ステップＳ１６０を図２２に示す初期コントーン画
像に適用することによって形成された最終コントーン画
像を示す。 図２４は、図１１に概要が示されたプロセ
スを図１に示す誤差拡散されたハーフトーン画像に適用
して得られた同様の最終コントーン画像を示す。
【００６９】図２５は、上記したように、図１と図２に
示す画像を印刷し走査した画像を示す。すなわち、図１
と図２に示す画像を、ハーフトーン画像発生器から直接
発生される電子的画像として本発明の汎用非スクリーニ
ング装置および方法に入力した。これに対し、図２５に
示す画像は、印刷されたハーフトーン画像を走査するこ
とによって本発明の汎用非スクリーニング装置および方
法に入力した。図２６は、図１１のステップを図２５に
示す走査されたハーフトーン画像に適用して発生された
コントーン画像を示す。図２３、図２４、および図２６
の比較から明らかなように、図１１の示した方法は、使
用されたハーフトーン化方法とは無関係にかつ画像が電
子的ハーフトーン画像であるか走査によって得たハーフ
トーン画像であるかに関わらず高品質のコントーン画像
を発生できる。
示す画像を印刷し走査した画像を示す。すなわち、図１
と図２に示す画像を、ハーフトーン画像発生器から直接
発生される電子的画像として本発明の汎用非スクリーニ
ング装置および方法に入力した。これに対し、図２５に
示す画像は、印刷されたハーフトーン画像を走査するこ
とによって本発明の汎用非スクリーニング装置および方
法に入力した。図２６は、図１１のステップを図２５に
示す走査されたハーフトーン画像に適用して発生された
コントーン画像を示す。図２３、図２４、および図２６
の比較から明らかなように、図１１の示した方法は、使
用されたハーフトーン化方法とは無関係にかつ画像が電
子的ハーフトーン画像であるか走査によって得たハーフ
トーン画像であるかに関わらず高品質のコントーン画像
を発生できる。
【００７０】図２３、図２４、および図２６のコントー
ンは、ハーフトーンおよび／または元のハーフトーン画
像を発生しおよび／または汎用ハーフトーン画像非スク
リーニング装置１００に入力するために使用された走査
パラメータの任意のものに関する初期知識なしに発生さ
れたものであることを理解されたい。しかし、もしその
ような情報が入手可能かまたは予測される場合、図１１
〜図１６に示した方法は、そのような情報を利用するた
めに容易に変更できる。特に、前記したＰ．Ｗ．ウォン
（Wong）による「誤差拡散のための逆ハーフトーン化お
よび核推定」（本明細書の一部を構成するものとして援
用）は、そのような情報を利用するための一つの方法を
開示している。
ンは、ハーフトーンおよび／または元のハーフトーン画
像を発生しおよび／または汎用ハーフトーン画像非スク
リーニング装置１００に入力するために使用された走査
パラメータの任意のものに関する初期知識なしに発生さ
れたものであることを理解されたい。しかし、もしその
ような情報が入手可能かまたは予測される場合、図１１
〜図１６に示した方法は、そのような情報を利用するた
めに容易に変更できる。特に、前記したＰ．Ｗ．ウォン
（Wong）による「誤差拡散のための逆ハーフトーン化お
よび核推定」（本明細書の一部を構成するものとして援
用）は、そのような情報を利用するための一つの方法を
開示している。
【００７１】上記の方向性フィルタ処理は、垂直、水
平、および傾斜方向の方向性フィルタ処理のそれぞれに
ついて下記のマスクを有するたたみ込みによって実施で
きる。垂 直 水 平 傾 斜 ０１０ ０００ １０１ ０１０ １１１ ０１０ ０１０ ０００ １０１ 本発明を上記の具体的な実施の形態に関連して記載した
が、多くの代替、変更、変種は当業者にとって明白であ
ることは明らかである。したがって、上記の本発明の実
施の形態は例示を意図するもので、限定を意図するもの
ではない。特許請求の範囲に定義された本発明の精神と
範囲から逸脱することなく種々の変更を行うことができ
る。
平、および傾斜方向の方向性フィルタ処理のそれぞれに
ついて下記のマスクを有するたたみ込みによって実施で
きる。垂 直 水 平 傾 斜 ０１０ ０００ １０１ ０１０ １１１ ０１０ ０１０ ０００ １０１ 本発明を上記の具体的な実施の形態に関連して記載した
が、多くの代替、変更、変種は当業者にとって明白であ
ることは明らかである。したがって、上記の本発明の実
施の形態は例示を意図するもので、限定を意図するもの
ではない。特許請求の範囲に定義された本発明の精神と
範囲から逸脱することなく種々の変更を行うことができ
る。
【図１】誤差拡散されたハーフトーン画像を示す写真で
ある。
ある。
【図２】３×３の集合ドットハーフトーン画像を示す写
真である。
真である。
【図３】汎用ハーフトーン画像非スクリーニング装置の
一つの実施の形態である。
一つの実施の形態である。
【図４】図３のＮステージハーフトーン画像デコンポー
ザの一つの実施の形態である。
ザの一つの実施の形態である。
【図５】サブバンドから連続色調画像を再構成する本発
明によるＮステージリコンポーザの一つの実施の形態で
ある。
明によるＮステージリコンポーザの一つの実施の形態で
ある。
【図６】図４の第１ステージのデコンポーザの一つの実
施の形態である。
施の形態である。
【図７】図５の第Ｎステージのリコンポーザの一つの実
施の形態である。
施の形態である。
【図８】図３のサブバンド雑音および／またはパターン
フィルターシステムの一つの実施の形態である。
フィルターシステムの一つの実施の形態である。
【図９】元の画像がどの様にして繰り返し的にサブバン
ド画像に分解されるのかを示す図である。
ド画像に分解されるのかを示す図である。
【図１０】隣接するレベルのサブバンド画像における係
数間の親子関係を示す図である。
数間の親子関係を示す図である。
【図１１】本発明の方法の一つの実施の形態を示すフロ
ーチャートである。
ーチャートである。
【図１２】図１１のハーフトーン画像を入力するステッ
プをより詳細に示す図である。
プをより詳細に示す図である。
【図１３】図１１のハーフトーン画像を分解するステッ
プをより詳細に示す図である。
プをより詳細に示す図である。
【図１４】図１１の選択的に雑音および／またはパター
ンを除去するステップをより詳細に示す図である。
ンを除去するステップをより詳細に示す図である。
【図１５】図１１の連続色調画像を再構成するステップ
をより詳細に示す図である。
をより詳細に示す図である。
【図１６】図１４の方向性フィルタを使用してサブバン
ド画像をフィルタ処理するステップをより詳細に示す図
である。
ド画像をフィルタ処理するステップをより詳細に示す図
である。
【図１７】ウエーブレット分解のための画像係数につい
ての親子関係を示す図である。
ての親子関係を示す図である。
【図１８】フーバーミニマックス関数を示す図である。
【図１９】図２のハーフトーン画像のウエーブレット分
解によって形成されたサブバンド画像を示す写真であ
る。
解によって形成されたサブバンド画像を示す写真であ
る。
【図２０】ツリークリッピングの後の図１９のサブバン
ド画像を示す写真である。
ド画像を示す写真である。
【図２１】方向性フィルタ処理の後の図２０のサブバン
ド画像を示す写真である。
ド画像を示す写真である。
【図２２】フィルタ処理されたサブバンド画像のウエー
ブレット再構成後の図２１の初期連続色調画像を示す写
真である。
ブレット再構成後の図２１の初期連続色調画像を示す写
真である。
【図２３】非線型フィルタ処理後の図２２の最終連続色
調画像を示す写真である。
調画像を示す写真である。
【図２４】図１の誤差拡散されたハーフトーン画像の最
終連続色調画像を示す写真である。
終連続色調画像を示す写真である。
【図２５】走査によって得られたハーフトーン画像を示
す写真である。
す写真である。
【図２６】図２５の走査によって得たハーフトーン画像
の最終連続色調画像を示す写真である。
の最終連続色調画像を示す写真である。
───────────────────────────────────────────────────── フロントページの続き (72)発明者 ジエボ ルオ アメリカ合衆国 14620 ニューヨーク州 ロチェスター ユニバーシティー パー ク 805 (72)発明者 ジガン ファン アメリカ合衆国 14580 ニューヨーク州 ウェブスター ヨークタウン ドライブ 153
Claims (1)
- 【請求項１】 ハーフトーン画像から連続色調画像を作
るための方法であって、 ハーフトーン画像を複数のサブバンド画像に分解する工
程と、 この複数のサブバンド画像から雑音とパターンの内の少
なくとも一つを選択的に除去する工程と、 複数のサブバンド画像を連続色調画像に再構成する工程
とを含むことを特徴とするハーフトーン画像からの連続
色調画像の作成方法。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US708014 | 1996-08-30 | ||
| US08/708,014 US5799112A (en) | 1996-08-30 | 1996-08-30 | Method and apparatus for wavelet-based universal halftone image unscreening |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH1098628A true JPH1098628A (ja) | 1998-04-14 |
| JP3964503B2 JP3964503B2 (ja) | 2007-08-22 |
Family
ID=24844043
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP24196397A Expired - Fee Related JP3964503B2 (ja) | 1996-08-30 | 1997-08-22 | ハーフトーン画像からの連続色調画像の作成方法 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US5799112A (ja) |
| EP (1) | EP0827333B1 (ja) |
| JP (1) | JP3964503B2 (ja) |
| DE (1) | DE69738897D1 (ja) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007116728A (ja) * | 2001-07-31 | 2007-05-10 | Ricoh Co Ltd | 圧縮された画像の画質の向上 |
| JP2007300624A (ja) * | 2006-04-28 | 2007-11-15 | Xerox Corp | 保存された２値画像データを増強する装置および方法 |
| WO2013128781A1 (ja) * | 2012-02-29 | 2013-09-06 | 独立行政法人科学技術振興機構 | 画像処理用ディジタルフィルタ、および、文字列傾斜錯視生成装置 |
| JP2014059916A (ja) * | 2012-05-14 | 2014-04-03 | Japan Science & Technology Agency | 印刷媒体、および、記録媒体 |
| CN108734680A (zh) * | 2018-05-29 | 2018-11-02 | 电子科技大学 | 一种自适应选择样本块大小的图像修复方法 |
Families Citing this family (52)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6108098A (en) * | 1995-12-28 | 2000-08-22 | Canon Kabushiki Kaisha | Image processing apparatus and method |
| JPH09214807A (ja) * | 1996-01-31 | 1997-08-15 | Canon Inc | 画像処理装置および画像処理方法 |
| US6163629A (en) * | 1998-01-30 | 2000-12-19 | Compaq Computer Corporation | Method for low complexity low memory inverse dithering |
| US6233060B1 (en) * | 1998-09-23 | 2001-05-15 | Seiko Epson Corporation | Reduction of moiré in screened images using hierarchical edge detection and adaptive-length averaging filters |
| US6768518B1 (en) * | 1998-12-16 | 2004-07-27 | Xerox Corporation | Method and apparatus for removing a checkerboard-like noise artifact from a captured composite NTSC video frame |
| US6580837B1 (en) * | 1999-12-07 | 2003-06-17 | Intel Corporation | Up-sampling decimated color plane data |
| US6973127B1 (en) * | 1999-12-23 | 2005-12-06 | Xvd Corporation | Apparatus and method for memory saving wavelet based video coding |
| US6864994B1 (en) * | 2000-01-19 | 2005-03-08 | Xerox Corporation | High-speed, high-quality descreening system and method |
| US6721458B1 (en) | 2000-04-14 | 2004-04-13 | Seiko Epson Corporation | Artifact reduction using adaptive nonlinear filters |
| US6751352B1 (en) * | 2000-05-25 | 2004-06-15 | Hewlett-Packard Development Company, L.P. | Method and apparatus for generating and decoding a visually significant barcode |
| US7107453B2 (en) * | 2000-05-25 | 2006-09-12 | Hewlett-Packard Development Company, L.P. | Authenticatable graphical bar codes |
| US6636645B1 (en) * | 2000-06-29 | 2003-10-21 | Eastman Kodak Company | Image processing method for reducing noise and blocking artifact in a digital image |
| US7187475B1 (en) * | 2000-09-28 | 2007-03-06 | Zoran Corporation | System and method for a robust de-screening filter |
| US6731821B1 (en) * | 2000-09-29 | 2004-05-04 | Hewlett-Packard Development Company, L.P. | Method for enhancing compressibility and visual quality of scanned document images |
| EP1339018B1 (en) * | 2000-11-30 | 2008-12-17 | Canon Kabushiki Kaisha | Image processing device, image processing method, recording medium and program |
| KR100574536B1 (ko) | 2000-11-30 | 2006-04-27 | 캐논 가부시끼가이샤 | 화상처리장치, 화상처리방법, 기억매체 및 프로그램 |
| US6938017B2 (en) * | 2000-12-01 | 2005-08-30 | Hewlett-Packard Development Company, L.P. | Scalable, fraud resistant graphical payment indicia |
| US6839152B2 (en) | 2000-12-06 | 2005-01-04 | Xerox Corporation | Adaptive filtering method and apparatus for descreening scanned halftoned image representations |
| US6738524B2 (en) | 2000-12-15 | 2004-05-18 | Xerox Corporation | Halftone detection in the wavelet domain |
| US7006255B2 (en) * | 2001-03-29 | 2006-02-28 | Sharp Laboratories Of America | Adaptive image filtering based on a distance transform |
| US6839450B2 (en) * | 2001-04-26 | 2005-01-04 | Hewlett-Packard Development Company, L.P. | Detecting halftone modulations embedded in an image |
| WO2002091219A1 (en) * | 2001-05-07 | 2002-11-14 | Hrl Laboratories, Llc | Method and apparatus for jointly optimizing linear signal processing filters with subband filters |
| US7187806B2 (en) * | 2001-08-23 | 2007-03-06 | Hewlett-Packard Development Company, L.P. | System and method for embedding information within a printed image using block error diffusion halftoning |
| KR100403601B1 (ko) * | 2001-12-21 | 2003-10-30 | 삼성전자주식회사 | 영상의 윤곽선 보정 장치 및 방법 |
| US7085015B2 (en) | 2002-02-01 | 2006-08-01 | Seiko Epson Corporation | Inverse halftoning for multi-level halftones |
| US7120305B2 (en) * | 2002-04-16 | 2006-10-10 | Ricoh, Co., Ltd. | Adaptive nonlinear image enlargement using wavelet transform coefficients |
| US7185388B2 (en) * | 2002-08-06 | 2007-03-06 | Harper Brush Works, Inc. | Power wave floor squeegee and handle connector |
| US6993185B2 (en) | 2002-08-30 | 2006-01-31 | Matsushita Electric Industrial Co., Ltd. | Method of texture-based color document segmentation |
| US7430334B2 (en) * | 2003-07-31 | 2008-09-30 | Hewlett Packard Development Company, L.P. | Digital imaging systems, articles of manufacture, and digital image processing methods |
| US20080159640A1 (en) * | 2003-11-06 | 2008-07-03 | Chaoqiang Liu | Document Image Encoding Decoding |
| US7653255B2 (en) | 2004-06-02 | 2010-01-26 | Adobe Systems Incorporated | Image region of interest encoding |
| US7639886B1 (en) | 2004-10-04 | 2009-12-29 | Adobe Systems Incorporated | Determining scalar quantizers for a signal based on a target distortion |
| WO2006106919A1 (ja) * | 2005-03-31 | 2006-10-12 | Nikon Corporation | 画像処理方法 |
| US7561747B2 (en) * | 2005-12-23 | 2009-07-14 | Xerox Corporation | Image data processing system and method |
| WO2008020487A1 (en) * | 2006-08-17 | 2008-02-21 | Fujitsu Limited | Image processing device, image processing program, and image processing method |
| WO2009081485A1 (ja) * | 2007-12-25 | 2009-07-02 | Fujitsu Limited | 画像処理装置、画像処理方法および画像処理プログラム |
| JP5213670B2 (ja) * | 2008-01-16 | 2013-06-19 | 三洋電機株式会社 | 撮像装置及びぶれ補正方法 |
| US8120679B2 (en) * | 2008-08-01 | 2012-02-21 | Nikon Corporation | Image processing method |
| US8077356B2 (en) * | 2008-09-30 | 2011-12-13 | Xerox Corporation | Content-aware halftone image resizing using iterative determination of energy metrics |
| US8023153B2 (en) * | 2008-09-30 | 2011-09-20 | Xerox Corporation | Content-aware halftone image resizing |
| US8098404B2 (en) * | 2008-09-30 | 2012-01-17 | Xerox Corporation | Content-aware uniform rosette color halftone image resizing using iterative determination of energy metrics |
| US8077357B2 (en) * | 2008-09-30 | 2011-12-13 | Xerox Corporation | Content-aware resizing of uniform rosette color halftone images |
| US8515182B2 (en) | 2009-02-11 | 2013-08-20 | Ecole De Technologie Superieure | Method and system for determining a quality measure for an image using multi-level decomposition of images |
| US8326046B2 (en) * | 2009-02-11 | 2012-12-04 | Ecole De Technologie Superieure | Method and system for determining structural similarity between images |
| US8515181B2 (en) | 2009-02-11 | 2013-08-20 | Ecole De Technologie Superieure | Method and system for determining a quality measure for an image using a variable number of multi-level decompositions |
| US8760572B2 (en) * | 2009-11-19 | 2014-06-24 | Siemens Aktiengesellschaft | Method for exploiting structure in sparse domain for magnetic resonance image reconstruction |
| US8983206B2 (en) | 2011-05-04 | 2015-03-17 | Ecole de Techbologie Superieure | Method and system for increasing robustness of visual quality metrics using spatial shifting |
| JP5773781B2 (ja) * | 2011-06-30 | 2015-09-02 | 株式会社東芝 | 超音波診断装置、画像処理装置及びプログラム |
| CA2899571A1 (en) * | 2013-01-30 | 2014-08-07 | Japan Science And Technology Agency | Digital filter for image processing, image processing apparatus, printing medium, recording medium, image processing method, and program |
| US9215345B2 (en) * | 2013-03-26 | 2015-12-15 | Xerox Corporation | Method and system for inverse halftoning utilizing inverse projection of predicted errors |
| CN108876741B (zh) * | 2018-06-22 | 2021-08-24 | 中国矿业大学(北京) | 一种复杂光照条件下的图像增强方法 |
| CN110677140B (zh) * | 2019-10-28 | 2021-07-02 | 北京航空航天大学 | 一种含未知输入和非高斯量测噪声的随机系统滤波器 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4630125A (en) * | 1983-06-01 | 1986-12-16 | Xerox Corporation | Unscreening of stored digital halftone images |
| US5027078A (en) * | 1989-10-10 | 1991-06-25 | Xerox Corporation | Unscreening of stored digital halftone images by logic filtering |
| US5339170A (en) * | 1991-12-31 | 1994-08-16 | Xerox Corporation | Image processing system and method employing hybrid filtering to provide improved reconstruction of continuous tone images from halftone screen-structured images |
| US5243444A (en) * | 1992-03-25 | 1993-09-07 | Xerox Corporation | Image processing system and method with improved reconstruction of continuous tone images from halftone images including those without a screen structure |
| US5239390A (en) * | 1992-06-05 | 1993-08-24 | Eastman Kodak Company | Image processing method to remove halftone screens |
| US5343309A (en) * | 1992-08-12 | 1994-08-30 | Xerox Corporation | Image processing system and method employing adaptive filtering to provide improved reconstruction of continuous tone images from halftone images including those without a screen structure |
| JP3293920B2 (ja) * | 1993-01-22 | 2002-06-17 | キヤノン株式会社 | 画像処理装置およびその方法 |
| US5471320A (en) * | 1994-05-11 | 1995-11-28 | Xerox Corporation | Stack filters for 1-to-N bit image processing in electronic printers |
-
1996
- 1996-08-30 US US08/708,014 patent/US5799112A/en not_active Expired - Lifetime
-
1997
- 1997-08-21 DE DE69738897T patent/DE69738897D1/de not_active Expired - Lifetime
- 1997-08-21 EP EP97306400A patent/EP0827333B1/en not_active Expired - Lifetime
- 1997-08-22 JP JP24196397A patent/JP3964503B2/ja not_active Expired - Fee Related
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007116728A (ja) * | 2001-07-31 | 2007-05-10 | Ricoh Co Ltd | 圧縮された画像の画質の向上 |
| JP2007300624A (ja) * | 2006-04-28 | 2007-11-15 | Xerox Corp | 保存された２値画像データを増強する装置および方法 |
| KR20140119176A (ko) * | 2012-02-29 | 2014-10-08 | 도쿠리츠교세이호진 가가쿠기쥬츠신고기코 | 화상처리용 디지털 필터 및 문자열 경사 착시 생성 장치 |
| KR20150103321A (ko) * | 2012-02-29 | 2015-09-09 | 고쿠리츠켄큐카이하츠호진 카가쿠기쥬츠신코키코 | 화상처리용 디지털 필터 및 문자열 경사 착시 생성 장치 |
| JP5385487B1 (ja) * | 2012-02-29 | 2014-01-08 | 独立行政法人科学技術振興機構 | スーパーハイブリッド画像生成装置、スーパーハイブリッド画像生成方法、印刷媒体製造方法、電子媒体製造方法、および、プログラム |
| JP2014017004A (ja) * | 2012-02-29 | 2014-01-30 | Japan Science & Technology Agency | 文字列傾斜錯視生成装置、文字列傾斜錯視生成方法、印刷媒体製造方法、電子媒体製造方法、および、プログラム |
| JP2013257901A (ja) * | 2012-02-29 | 2013-12-26 | Japan Science & Technology Agency | 画像処理用ディジタルフィルタ、画像生成装置、画像生成方法、ディジタルフィルタ作成方法、および、プログラム |
| US9721331B2 (en) | 2012-02-29 | 2017-08-01 | National Institute Of Japan Science And Technology Agency | Digital filter, and image generating, superhybrid image generating, electronic medium manufacturing, and letter-row tilt illusion generating apparatus, method and program |
| WO2013128781A1 (ja) * | 2012-02-29 | 2013-09-06 | 独立行政法人科学技術振興機構 | 画像処理用ディジタルフィルタ、および、文字列傾斜錯視生成装置 |
| CN104137144A (zh) * | 2012-02-29 | 2014-11-05 | 独立行政法人科学技术振兴机构 | 图像处理用数字滤波器及字符串倾斜错觉生成装置 |
| TWI587244B (zh) * | 2012-02-29 | 2017-06-11 | 獨立行政法人科學技術振興機構 | 影像處理用數位濾波器、影像產生裝置、超級混合影像產生裝置、影像產生方法、數位濾波器作成方法、超級混合影像產生方法、印刷媒體製造方法、電子媒體製造方法及程式 |
| JP2014059916A (ja) * | 2012-05-14 | 2014-04-03 | Japan Science & Technology Agency | 印刷媒体、および、記録媒体 |
| US9292910B2 (en) | 2012-05-14 | 2016-03-22 | National Institute Of Japan Science And Technology Agency | Image processing apparatus, image processing method, program, printing medium, and recording medium |
| JP2016139430A (ja) * | 2012-05-14 | 2016-08-04 | 国立研究開発法人科学技術振興機構 | 画像処理装置、画像処理方法、および、画像処理プログラム |
| KR20150008877A (ko) * | 2012-05-14 | 2015-01-23 | 도쿠리츠교세이호진 가가쿠기쥬츠신고기코 | 화상 처리 장치, 화상 처리 방법, 프로그램, 인쇄 매체, 및 기록 매체 |
| JP5599520B2 (ja) * | 2012-05-14 | 2014-10-01 | 独立行政法人科学技術振興機構 | 画像処理装置、画像処理方法、および、プログラム |
| CN108734680A (zh) * | 2018-05-29 | 2018-11-02 | 电子科技大学 | 一种自适应选择样本块大小的图像修复方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| EP0827333B1 (en) | 2008-08-13 |
| US5799112A (en) | 1998-08-25 |
| EP0827333A2 (en) | 1998-03-04 |
| JP3964503B2 (ja) | 2007-08-22 |
| EP0827333A3 (en) | 1999-12-29 |
| DE69738897D1 (de) | 2008-09-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3964503B2 (ja) | ハーフトーン画像からの連続色調画像の作成方法 | |
| Stevenson | Inverse halftoning via MAP estimation | |
| US5323247A (en) | Method and apparatus for halftoning and inverse halftoning and the transmission of such images | |
| US5343309A (en) | Image processing system and method employing adaptive filtering to provide improved reconstruction of continuous tone images from halftone images including those without a screen structure | |
| US7599095B2 (en) | Wavelet-based image processing path | |
| Luo et al. | A robust technique for image descreening based on the wavelet transform | |
| EP0562813B1 (en) | System and method for converting halftone images to continuous-tone images | |
| US6101285A (en) | Filter for producing continuous tone images from halftone digital images data | |
| US9215345B2 (en) | Method and system for inverse halftoning utilizing inverse projection of predicted errors | |
| CN100405814C (zh) | 以一维多阶网点分配为基础的数字半色调方法 | |
| Foi et al. | Inverse halftoning based on the anisotropic LPA-ICI deconvolution | |
| Shen et al. | A robust nonlinear filtering approach to inverse halftoning | |
| Walker et al. | Image denoising using tree-based wavelet subband correlations and shrinkage | |
| Mese et al. | Tree-structured method for LUT inverse halftoning and for image halftoning | |
| Egiazarian et al. | Adaptive denoising and lossy compression of images in transform domain | |
| Chen et al. | An adaptive inverse halftoning algorithm | |
| EP1569168B1 (en) | Finite impulse response filter method and apparatus | |
| Zhang et al. | Image inverse halftoning and descreening: a review | |
| Atkins et al. | Halftone postprocessing for improved rendition of highlights and shadows | |
| Kite | Design and quality assessment of forward and inverse error diffusion halftoning algorithms | |
| CN1189841C (zh) | 多分辨率表示的误差扩散 | |
| Luo et al. | Universal descreening technique via wavelet analysis | |
| Yang et al. | Inverse halftoning based on sparse representation with boosted dictionary | |
| Kuo et al. | Color halftone document segmentation and descreening | |
| Shen et al. | Inverse halftoning via robust nonlinear filtering |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20040803 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20061031 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070124 |
|
| 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: 20070424 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070524 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110601 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120601 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130601 Year of fee payment: 6 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |