JP2009514273A - マルチビット深度の画像装置でｆｍ／ａｍ混在網点画像を処理する方法 - Google Patents
マルチビット深度の画像装置でｆｍ／ａｍ混在網点画像を処理する方法 Download PDFInfo
- Publication number
- JP2009514273A JP2009514273A JP2008536906A JP2008536906A JP2009514273A JP 2009514273 A JP2009514273 A JP 2009514273A JP 2008536906 A JP2008536906 A JP 2008536906A JP 2008536906 A JP2008536906 A JP 2008536906A JP 2009514273 A JP2009514273 A JP 2009514273A
- Authority
- JP
- Japan
- Prior art keywords
- pixel
- bit
- output
- value
- diffusion
- 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
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000009792 diffusion process Methods 0.000 claims abstract description 50
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 9
- 230000009977 dual effect Effects 0.000 claims abstract description 7
- 238000012216 screening Methods 0.000 claims description 25
- 238000003384 imaging method Methods 0.000 claims description 13
- 238000009826 distribution Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 230000007246 mechanism Effects 0.000 claims description 4
- 230000002457 bidirectional effect Effects 0.000 claims description 3
- 238000012896 Statistical algorithm Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 8
- 238000005516 engineering process Methods 0.000 description 8
- 230000007480 spreading Effects 0.000 description 2
- 238000003892 spreading Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000001771 impaired effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
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/40087—Multi-toning, i.e. converting a continuous-tone signal for reproduction with more than two discrete brightnesses or optical densities, e.g. dots of grey and black inks on white paper
-
- 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/405—Halftoning, i.e. converting the picture signal of a continuous-tone original into a corresponding signal showing only two levels
- H04N1/4051—Halftoning, i.e. converting the picture signal of a continuous-tone original into a corresponding signal showing only two levels producing a dispersed dots halftone pattern, the dots having substantially the same size
- H04N1/4052—Halftoning, i.e. converting the picture signal of a continuous-tone original into a corresponding signal showing only two levels producing a dispersed dots halftone pattern, the dots having substantially the same size by error diffusion, i.e. transferring the binarising error to neighbouring dot decisions
- H04N1/4053—Halftoning, i.e. converting the picture signal of a continuous-tone original into a corresponding signal showing only two levels producing a dispersed dots halftone pattern, the dots having substantially the same size by error diffusion, i.e. transferring the binarising error to neighbouring dot decisions with threshold modulated relative to input image data or vice versa
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Discrete Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
- Color, Gradation (AREA)
Abstract
【解決手段】本発明は、マルチビット深度の画像装置でFM/AM 混在網点画像を処理する方法を開示し、画像のハードコピーの分野でのハーフトーンドットを生成する技術に属する。先行技術では、出力装置のための誤差拡散法に基づいたマルチビット画像深度での混在格子点を制御する方法の実現は、誤差拡散による影響を回避できないので、マルチビット画像深度での混在格子点の出力は、装置の必要条件を達成することができない。本方法は、ネッティング(netting) の際のオリジナル誤差拡散のデュアルフィードバック混在ネッティング(netting) 法に基づいた、マルチビット混在格子点の動的制御アルゴリズムを適用する。本発明は、マルチビット画像深度の出力装置から質が高く豊富な階調でのマルチビット網点画像を出力することができ、オリジナル１ビット深度の装置によって出力された混在格子点の縁での鋸歯状現象を効果的に解決することができ、低解像度で高解像度の効果を有するFM/AM 混在格子点の連続階調出力を保障する。
Description
本発明は、画像の複写のためのハードコピーの分野でハーフトーンドットを作成する方法に関し、特にマルチビット深度の画像装置でFM-AM 混在網点画像を処理する方法に関する。
画像の複写のためのハードコピーは、一般的にプリンタ、及び改良された印刷及び製版装置のためのスクリーニング及び製版技術に関する。ハードコピーで画像を複写するために用いられるスクリーニング技術は、デジタル画像ハーフトーン技術とも呼ばれる。デジタル画像ハーフトーン技術は振幅変調(AM)スクリーニングと周波数変調(FM)スクリーニングとを備える。振幅変調スクリーニング技術は、集められたドットの配列ディザリングとも呼ばれ、作成された網点画像の色のついた点が、ドットと呼ばれるカラー領域のクラスタを形成すべく、ペアで幾何学的に集められることを特徴とする。この技術は、原画像のグレーレベルを表わすためにドットのサイズを制御するので、ドットは振幅変調ドットと呼ばれる。
FM及びAMスクリーンに基づいた混在スクリーニング技術について、本出願者は、中国特許出願「デュアルフィードバックに基づく誤差拡散を用いた周波数変調スクリーニングのための方法」（2005年９月14日に公開された出願番号No.200510068127.8 ）を出願しており、開示されたFM-AM 混在ドットの技術は、主にシングルビットの画像深度を備えた装置に基づいている。実際の出力中では、シングルビットの装置は、ドットを出力するためにわずか２階調(白:255; 黒:0)しか備えていないので、鋸歯状現象が、シングルビット深度の画像装置によって出力されたFM-AM 混在ドットの縁、特に文字の縁に一般的に存在する。更に、混在ドットに関するAMの特徴の理由から、ドットのサイズによる画像の鮮明化を損なうため、シングルビット深度の画像装置によって出力された画像に深刻な粒状感及び質感を与え、全体画像での滑らかさ及び出力画像の質に不利な影響を与える。
シングルビット深度の画像装置の欠点を克服するために、マルチビット深度の画像装置が市場に出ている。この装置の出力階調は、画像深度に応じて２の累乗まで増加する。一般的な２ビット又は４ビットの出力装置は夫々４又は16階調を表わすことができる。マルチビット深度の画像装置の最終目的は、低解像度での画像の質の問題を克服すべく、低解像度で高解像度の出力効果を得ることである。更に、マルチビット深度の画像装置に組み込まれたPWM 階調オフセット技術と組み合わせて、ハーフトーンドットが安定して出力され、装置によって出力されたハーフトーンドットの最適化された質を得ることである。
ハードコピーの分野では、マルチビット深度の画像装置に適したハーフトーンスクリーニング技術は、主にAMスクリーニング技術である。AMドットが本来規則的な分布を有し、ドットのサイズを制御できるので、シングルビット深度での従来のAMスクリーニング技術に関して、マルチビット深度の画像装置に適したマルチビット深度での対応するスクリーニング技術を発展させることは容易である。
誤差拡散に基づいた「デュアルフィードバックに基づく誤差拡散を用いた周波数変調スクリーニングのための方法」の出願では、リアルタイムの動的ハーフトーンスクリーニング、ドットの分布の不均一性、及び様々な位置でドットの形状を制御するためのランダムな変更により、マルチビット深度でのドット生成技術を導き出すのは困難である。生成のための効果的な方法は存在せず、従来の方法は改善される必要がある。
米国特許出願公開第２００４／１３０７５３号明細書
先行技術の欠点を克服するために、本発明は、マルチビット深度の画像装置でFM-AM 混在網点画像を処理する方法を提供する。前記方法は、誤差拡散に基づきマルチビット画像深度でFM-AM 混在ドットを生成する技術に関する。マルチビットのドットを生成する技術は、従来のFM-AM 混在ドットがマルチビット画像深度の特徴を有するように、混在ドットを生成する従来の技術に基づいて導き出される。従って、マルチビットFM-AM 混在ドットが、マルチビット深度の画像装置によって出力でき、マルチビット深度の画像装置の利点を有効に利用して、混在ドットの出力の質を改善することができる。
上記の目的を達成するために、本発明は、マルチビット深度の画像装置でFM/AM 混在網点画像を処理する方法を提供し、前記方法は、以下を備える。
ステップ1) 装置のビット深度n を考慮して間隔[0,255]を2n-1階調に平均的に分割する
[0, R1], (R1, R2], … , (Ri-1, Ri], … (R2 n -2, 255]
但し、ビット出力のドット・マトリクスの対応する範囲は、
(0, Out1), (Out1, Out2), … , (Outi-1, Outi), … (Out2 n -2, 11…1)
であり、Outiはn ビット深度の２進表示であり、
各階調の中心点の閾値Miは、階調のための閾値比較パラメータとしてサンプリングされる
ステップ2) 間隔[0,255]にn ビット画像深度での出力確率閾値を設定する
但し、ステップ2)は以下を含む
2n-2個の整数L1,L2,…,Li,…,L2 n -2をn ビット画像深度での出力確率閾値として設定すること及び
装置の必要条件に基づき確率閾値を調整すること
ステップ3) デュアルフィードバックの誤差拡散を用いたFM-AM 混在スクリーニング処理に基づき、2n-1階調(Ri-1,Ri) でドットを夫々処理し、詳細な処理として以下を含む
ステップ(a) 原画像での注目画素の最後の入力値g"(m,t) に対して閾値比較の操作T を実行し、その後、操作の結果を網点画像のための注目画素の対応する値b(m,t)に変換する
但し、ステップ(a) は、原画像を走査するとき、双方向走査の処理を用い、ある行が左から右に走査されるとき、続いて次の行は右から左に走査される
ステップ(b) 画素の値b(m,t)を、同一画素の入力値g'(m,t) と比較して、b(m,t)とg'(m,t) との差を得る、但し、差は誤差値e(m,t)であり、入力値g'(m,t) は閾値を得るために用いられる
ステップ(c) 誤差値e(m,t)を、誤差拡散フィルタe により予め設定された重み付け分布係数と乗算して、その後、乗算結果を注目画素の周囲の未処理の画素に拡散する、但し、原画像の対応する画素の入力値g'(m,t) を得るために、注目画素の周囲の未処理の画素への各拡散結果を、原画像の対応する画素の原入力値g(m,t)に加重的に加算する
但し、誤差拡散フィルタe は以下に示す拡散原理と重み付け分布係数とを用いる
ステップ1) 装置のビット深度n を考慮して間隔[0,255]を2n-1階調に平均的に分割する
[0, R1], (R1, R2], … , (Ri-1, Ri], … (R2 n -2, 255]
但し、ビット出力のドット・マトリクスの対応する範囲は、
(0, Out1), (Out1, Out2), … , (Outi-1, Outi), … (Out2 n -2, 11…1)
であり、Outiはn ビット深度の２進表示であり、
各階調の中心点の閾値Miは、階調のための閾値比較パラメータとしてサンプリングされる
ステップ2) 間隔[0,255]にn ビット画像深度での出力確率閾値を設定する
但し、ステップ2)は以下を含む
2n-2個の整数L1,L2,…,Li,…,L2 n -2をn ビット画像深度での出力確率閾値として設定すること及び
装置の必要条件に基づき確率閾値を調整すること
ステップ3) デュアルフィードバックの誤差拡散を用いたFM-AM 混在スクリーニング処理に基づき、2n-1階調(Ri-1,Ri) でドットを夫々処理し、詳細な処理として以下を含む
ステップ(a) 原画像での注目画素の最後の入力値g"(m,t) に対して閾値比較の操作T を実行し、その後、操作の結果を網点画像のための注目画素の対応する値b(m,t)に変換する
但し、ステップ(a) は、原画像を走査するとき、双方向走査の処理を用い、ある行が左から右に走査されるとき、続いて次の行は右から左に走査される
ステップ(b) 画素の値b(m,t)を、同一画素の入力値g'(m,t) と比較して、b(m,t)とg'(m,t) との差を得る、但し、差は誤差値e(m,t)であり、入力値g'(m,t) は閾値を得るために用いられる
ステップ(c) 誤差値e(m,t)を、誤差拡散フィルタe により予め設定された重み付け分布係数と乗算して、その後、乗算結果を注目画素の周囲の未処理の画素に拡散する、但し、原画像の対応する画素の入力値g'(m,t) を得るために、注目画素の周囲の未処理の画素への各拡散結果を、原画像の対応する画素の原入力値g(m,t)に加重的に加算する
但し、誤差拡散フィルタe は以下に示す拡散原理と重み付け分布係数とを用いる
但し、**は注目画素の位置を表わし、他の位置での各算術比は、注目画素に対する拡散重み付け係数を表わし、拡散重み付け係数は[0,1] に属し、以下の式を満たす
2×d1+4×d2+2×d3+2×d4+2×d5∈[0, 1]
ステップ(d) 注目画素を囲む対応する未処理の画素に処理結果を夫々拡散し、対応する画素の最後の入力値g"(m,t) を得るために、原画像の対応する画素の原入力値g(m,t)に、拡散された各処理結果を加重的に加算する、但し、処理結果は、拡散フィルタw を用いて注目画素の出力値b(m,t)に乗算操作を実行し、ディザリングアルゴリズムで乗算操作の結果を処理することにより得られ、ステップ(d) はステップ(b) 及び(c) と並行して実行される
但し、拡散フィルタw の拡散モードは以下に示すように設定される
2×d1+4×d2+2×d3+2×d4+2×d5∈[0, 1]
ステップ(d) 注目画素を囲む対応する未処理の画素に処理結果を夫々拡散し、対応する画素の最後の入力値g"(m,t) を得るために、原画像の対応する画素の原入力値g(m,t)に、拡散された各処理結果を加重的に加算する、但し、処理結果は、拡散フィルタw を用いて注目画素の出力値b(m,t)に乗算操作を実行し、ディザリングアルゴリズムで乗算操作の結果を処理することにより得られ、ステップ(d) はステップ(b) 及び(c) と並行して実行される
但し、拡散フィルタw の拡散モードは以下に示すように設定される
但し、走査方向は左から右であり、**は注目画素の位置を表わし、他の位置での各パラメータは、注目画素に対する拡散重み付け係数を表わし、パラメータは[0,1] に属し、以下の式を満たす
wsum=(w0+w1+w2+w3)∈[0, 1]
ステップ(d) での拡散フィルタw のためのディザリングアルゴリズムは、以下に示す通りである
wsum=(w0+w1+w2+w3)∈[0, 1]
ステップ(d) での拡散フィルタw のためのディザリングアルゴリズムは、以下に示す通りである
但し、fRand はディザリングを微調整するためのパラメータであり、R(m,t)は注目ドットを走査するためのランダム値でのパラメータであり、R_MAX はランダムパラメータR(i)の最大値であり、cDither は、ディザリングの振幅を調整するためのパラメータであり、振幅変調の動作を決定し、dw0 〜dw3 はディザリング後の異なる方向での拡散フィルタw の拡散重み付け係数である
ステップ(e) ステップ(a)-(d) を、すべての画素の原入力値g(m,t)が処理されるまで繰り返す
ステップ4) 階調(Ri-1,Ri) での混在スクリーニングが達成されたとき、ドットの形状を制御するために、注目ドットの形状を制御するための確率閾値Fi及び蓄積された値ShapeCurを考慮して、隣接した出力グレーレベルでの公知の動的階調可変出力機構及び動的統計アルゴリズムを用いることによって、出力ドット・マトリクスデータを動的に算出する
但し、ステップ4)での動的階調可変出力機構は、以下を含む
ShapeCurを考慮して注目ドットの偽似乱数値を生成する
Fi=random(ShapeCur) 式１
但し、偽似乱数関数randomは、コンパイル環境で自動的に生成され、Fi∈[0,255] であり、出力ドット・マトリクスデータは動的に算出される
ステップ(e) ステップ(a)-(d) を、すべての画素の原入力値g(m,t)が処理されるまで繰り返す
ステップ4) 階調(Ri-1,Ri) での混在スクリーニングが達成されたとき、ドットの形状を制御するために、注目ドットの形状を制御するための確率閾値Fi及び蓄積された値ShapeCurを考慮して、隣接した出力グレーレベルでの公知の動的階調可変出力機構及び動的統計アルゴリズムを用いることによって、出力ドット・マトリクスデータを動的に算出する
但し、ステップ4)での動的階調可変出力機構は、以下を含む
ShapeCurを考慮して注目ドットの偽似乱数値を生成する
Fi=random(ShapeCur) 式１
但し、偽似乱数関数randomは、コンパイル環境で自動的に生成され、Fi∈[0,255] であり、出力ドット・マトリクスデータは動的に算出される
本発明は以下の効果を有する。本発明の方法によれば、マルチビット混在ハーフトーンドットが、デュアルフィードバックに基づいた誤差拡散を用いて従来のFM-AM 混在スクリーニング方法に基づき、装置の画像深度を考慮して生成されることが可能であり、マルチビットの画像装置を有効に利用することができる。更に、高解像度の効果を有するFM-AM 混在ハーフトーンスクリーニングを、低解像度で得ることができる。本発明の方法は、ドットの縁での鋸歯状現象を最適化し、文字の縁での滑らかさの欠如に関する問題を解決し、ドットの出力階調を改善し、階調の連続した滑らかさを保証する。
本発明を、添付図面及び実施形態を参照して以下に詳細に説明する。
本発明は、マルチビット深度の画像装置でのFM-AM 混在網点画像を処理する方法を提供し、前記方法は、（図１に示すように）マルチビット深度でドットを制御するための動的アルゴリズムを用いる。本実施形態では、２ビット深度の画像装置(n=2)が用いられる。 詳細な解決法を以下に説明する。
ステップ１では、（図２に示すように）間隔[0-255] を（装置のビット深度n が２であり、2n-1=3であるため）３階調に平均的に分割する。
[0, 84], (84, 171] 及び (171, 255]
ビット出力のドット・マトリクスの対応する範囲は、
(11, 10), (10, 01) 及び (01, 00)
である。
本実施形態での実際の出力では、黒一色は０と設定され、白一色は255 と設定される。
３階調の中心点の閾値Mi(M1=42,M2=127 及びM3=212) は、各階調のための閾値比較パラメータとしてサンプリングされる。
[0, 84], (84, 171] 及び (171, 255]
ビット出力のドット・マトリクスの対応する範囲は、
(11, 10), (10, 01) 及び (01, 00)
である。
本実施形態での実際の出力では、黒一色は０と設定され、白一色は255 と設定される。
３階調の中心点の閾値Mi(M1=42,M2=127 及びM3=212) は、各階調のための閾値比較パラメータとしてサンプリングされる。
ステップ２では、間隔[0,255] にn ビット画像深度での出力確率閾値として（2n-2=2であるため）２個の整数F1=8及びF2=24 を設定する。11、10、01及び00の出力確率は、２ビット混在ドットの出力効果を得るように閾値によって左右される。
ステップ３では、デュアルフィードバックの誤差拡散を用いた公知のFM-AM 混在スクリーニング法に基づき、３階調でドットを夫々処理する。詳細な処理を以下に説明する。
(a) 原画像を走査し、入力する。閾値比較の操作T は原画像での注目画素の画素値g"(m,t) に対して実行される。操作の結果は、網点画像のために注目画素の対応する値b(m,t)に変換される。
ランダム分布でのドットの走査方向及び周波数によって引き起こされる（一般的に「モワレ」と呼ばれる）現象を回避するために、本実施形態は原画素のデータを走査するために双方向走査方法を用いる。
(b) 注目画素の値b(m,t)を、注目画素の入力値g'(m,t) と比較して、差（つまりb(m,t)とg'(m,t) との誤差値e(m,t)）を得る。
(c) 誤差値e(m,t)を、拡散フィルタe によりある重み付け分布係数と乗算する。その後、乗算結果を注目画素の周囲の未処理の画素に拡散する。対応する画素の入力値g'(m,t) を得るために、注目画素の周囲の未処理の画素への各拡散結果を、原画像の対応する画素の原入力値g(m,t)に加重的に加算する。誤差拡散フィルタは以下に示す拡散原理と重み付け分布係数とを用いる。
(a) 原画像を走査し、入力する。閾値比較の操作T は原画像での注目画素の画素値g"(m,t) に対して実行される。操作の結果は、網点画像のために注目画素の対応する値b(m,t)に変換される。
ランダム分布でのドットの走査方向及び周波数によって引き起こされる（一般的に「モワレ」と呼ばれる）現象を回避するために、本実施形態は原画素のデータを走査するために双方向走査方法を用いる。
(b) 注目画素の値b(m,t)を、注目画素の入力値g'(m,t) と比較して、差（つまりb(m,t)とg'(m,t) との誤差値e(m,t)）を得る。
(c) 誤差値e(m,t)を、拡散フィルタe によりある重み付け分布係数と乗算する。その後、乗算結果を注目画素の周囲の未処理の画素に拡散する。対応する画素の入力値g'(m,t) を得るために、注目画素の周囲の未処理の画素への各拡散結果を、原画像の対応する画素の原入力値g(m,t)に加重的に加算する。誤差拡散フィルタは以下に示す拡散原理と重み付け分布係数とを用いる。
但し、**は注目画素の位置を表わし、他の位置での各算術比は、注目画素に対する拡散重み付け係数を表わす。
本実施形態では、係数を以下のように設定する。
d1=1/44, d2=2/44, d3=5/44, d4=4/44, d5=8/44
ステップ(b) 及び(c) は、第１拡散フィードバック操作を達成し、誤差拡散の基本原理を実施する。
(d) このステップはステップ(b) 及び(c) と並行して実行される。処理結果を、注目画素を囲む対応する未処理の画素に夫々拡散し、対応する画素の最後の入力値g"(m,t) を得るために、ステップ(b) 及び(c) での拡散された各処理結果及び誤差拡散の各パラメータを、原画像の対応する画素の原入力値g(m,t)に加重的に加算する、但し、処理結果は、第２拡散フィルタw を用いて注目画素の出力値b(m,t)に対する乗算操作を実行し、ディザリングアルゴリズムで乗算操作の結果を処理することにより得られる。第２拡散フィルタw の拡散モードは以下のように設定される。
本実施形態では、係数を以下のように設定する。
d1=1/44, d2=2/44, d3=5/44, d4=4/44, d5=8/44
ステップ(b) 及び(c) は、第１拡散フィードバック操作を達成し、誤差拡散の基本原理を実施する。
(d) このステップはステップ(b) 及び(c) と並行して実行される。処理結果を、注目画素を囲む対応する未処理の画素に夫々拡散し、対応する画素の最後の入力値g"(m,t) を得るために、ステップ(b) 及び(c) での拡散された各処理結果及び誤差拡散の各パラメータを、原画像の対応する画素の原入力値g(m,t)に加重的に加算する、但し、処理結果は、第２拡散フィルタw を用いて注目画素の出力値b(m,t)に対する乗算操作を実行し、ディザリングアルゴリズムで乗算操作の結果を処理することにより得られる。第２拡散フィルタw の拡散モードは以下のように設定される。
但し、走査方向は左から右であり、**は注目画素の位置を表わし、他の位置での各パラメータは、注目画素に対する拡散重み付け係数を表わす。パラメータは[0,1] に属し、以下の式を満たす。
wsum=(w0+w1+w2+w3)∈[0, 1]
本実施形態でのディザリングアルゴリズムは以下の式を用いる。
wsum=(w0+w1+w2+w3)∈[0, 1]
本実施形態でのディザリングアルゴリズムは以下の式を用いる。
この式では、fRand はディザリングを微調整するためのパラメータであり、R(m,t)は注目ドットを走査するためのランダム値でのパラメータであり、R_MAX はランダムパラメータR(i)の最大値であり、cDither は、ディザリングの振幅を調整するためのパラメータであり、振幅変調の動作を決定し、dw0 〜dw3 はディザリング後の異なる方向でのフィルタw の拡散重み付け係数である。
このステップ(d) は、第２拡散フィードバック操作を達成し、周波数変調スクリーンの振幅変調特性を実施する。
本実施形態では、係数は以下のように設定される。
w0=w2=0.175, w1=w3=0.025, その結果 wsum=0.4
cDither=0.2
スクリーニングの間、周波数変調ハーフトーンドットのサイズは、wsumを調整することにより変更され、周波数変調ハーフトーンドットの形状はw0〜w3の値を調整することにより制御される。
(e) ステップ(a)-(d) を、すべての画素の原入力値g(m,n)が処理されるまで繰り返す。
このステップ(d) は、第２拡散フィードバック操作を達成し、周波数変調スクリーンの振幅変調特性を実施する。
本実施形態では、係数は以下のように設定される。
w0=w2=0.175, w1=w3=0.025, その結果 wsum=0.4
cDither=0.2
スクリーニングの間、周波数変調ハーフトーンドットのサイズは、wsumを調整することにより変更され、周波数変調ハーフトーンドットの形状はw0〜w3の値を調整することにより制御される。
(e) ステップ(a)-(d) を、すべての画素の原入力値g(m,n)が処理されるまで繰り返す。
ステップ４では、３階調での混在スクリーニングが達成されるとき、動的階調可変出力機構が、注目ドットの形状を制御するための確率閾値(L1=8,L2=24)及び蓄積された値ShapeCurを考慮して用いられる。このステップのアルゴリズムは、以下の式を用いる（注目ドットの偽似乱数値はShapeCurと式１とに基づいて生成される）。
F=random(ShapeCur)
但し、偽似乱数関数randomは、コンパイル環境で自動的に生成することができる。(F∈[0,255])
F=random(ShapeCur)
但し、偽似乱数関数randomは、コンパイル環境で自動的に生成することができる。(F∈[0,255])
出力ドット・マトリクスデータは、式２に基づいて動的に算出される。ゾーン[0,84]での処理は図３に示され、ゾーン(84,171]での処理は図４に示され、ゾーン(171,255] での処理は図５に示される。
２ビット深度の画像装置を用いて、２ビットFM-AM 混在スクリーニングが本実施形態での上記のステップを実行することにより達成することができる。図６はスクリーニングの効果を示す。
本発明は上述した実施形態に制限されない。本発明での技術的解決法に応じて当業者によって得られた他の実施形態は、本発明の技術改良の範囲内である。
Claims (8)
- マルチビット深度の画像装置でFM/AM 混在網点画像を処理する方法において、
ステップ1) 前記装置のビット深度n を考慮して間隔[0,255]を2n-1階調に平均的に分割する
[0, R1], (R1, R2], … , (Ri-1, Ri], … (R2 n -2, 255]
但し、i は正の整数であり2n-2未満であり、ビット出力のドット・マトリクスの対応する範囲は、
(0, Out1), (Out1, Out2), … , (Outi-1, Outi), … (Out2 n -2, 11…1)
であり、Outiはn ビット深度の２進表示であり、
各階調の中心点の閾値Miは、階調のための閾値比較パラメータとしてサンプリングされる
ステップ2) 間隔[0,255]にn ビット画像深度での出力確率閾値を設定する
ステップ3) デュアルフィードバックの誤差拡散を用いたFM-AM 混在スクリーニング処理に基づき、2n-1階調(Ri-1,Ri) でドットを夫々処理する
ステップ4) 階調(Ri-1,Ri) での混在スクリーニングが達成されたとき、注目ドットの形状を制御するための確率閾値Fi及び蓄積された値ShapeCurを考慮して、隣接した出力グレーレベルでの動的統計アルゴリズムを用いることによって、出力ドット・マトリクスデータを動的に算出する
を備えることを特徴とする方法。 - 前記ステップ2)は、
2n-2個の整数L1,L2,…,Li,…,L2 n -2をn ビット画像深度での出力確率閾値として設定すること及び
前記装置の必要条件に基づき前記確率閾値を調整すること
を含むことを特徴とする請求項１に記載の方法。 - 前記ステップ3)は、
ステップ(a) 原画像での注目画素の最後の入力値g"(m,t) に対して閾値比較の操作T を実行し、その後、前記操作の結果を網点画像のための注目画素の対応する値b(m,t)に変換する
ステップ(b) 注目画素の値b(m,t)を、注目画素の入力値g'(m,t) と比較して、b(m,t)とg'(m,t) との差を得る、但し、前記差は誤差値e(m,t)であり、前記入力値g'(m,t) は閾値を得るために用いられる
ステップ(c) 誤差値e(m,t)を、拡散フィルタe により予め設定された重み付け分布係数と乗算して、その後、前記乗算結果を注目画素の周囲の未処理の画素に拡散する、但し、原画像の対応する画素の入力値g'(m,t) を得るために、注目画素の周囲の未処理の画素への各拡散結果を、原画像の対応する画素の原入力値g(m,t)に加重的に加算する
ステップ(d) 注目画素を囲む対応する未処理の画素に処理結果を夫々拡散し、対応する画素の最後の入力値g"(m,t) を得るために、原画像の対応する画素の原入力値g(m,t)に、拡散された各処理結果を加重的に加算する、但し、前記処理結果は、拡散フィルタw を用いて注目画素の出力値b(m,t)に乗算操作を実行し、ディザリングアルゴリズムで前記乗算操作の結果を処理することにより得られ、前記ステップ(d) は前記ステップ(b) 及び(c) と並行して実行される
ステップ(e) 前記ステップ(a)-(d) を、すべての画素の原入力値g(m,t)が処理されるまで繰り返す
を含むことを特徴とする請求項１に記載の方法。 - 前記ステップ(a) は、原画像を走査するとき、双方向走査の処理を用い、ある行が左から右に走査されるとき、続いて次の行は右から左に走査されることを特徴とする請求項３に記載の方法。
- 前記誤差拡散フィルタe は、以下に示す拡散原理と重み付け分布モードとを用いることを特徴とする請求項４に記載の方法。
但し、**は注目画素の位置を表わし、他の位置での各算術比は、注目画素に対する拡散重み付け係数を表わし、拡散重み付け係数は[0,1] に属し、以下の式を満たす
2×d1+4×d2+2×d3+2×d4+2×d5∈[0, 1] - 前記拡散フィルタw の拡散モードは以下に示すように設定されることを特徴とする請求項５に記載の方法。
但し、走査方向は左から右であり、**は注目画素の位置を表わし、他の位置での各パラメータは、注目画素に対する拡散重み付け係数を表わし、パラメータは[0,1] に属し、以下の式を満たす
wsum=(w0+w1+w2+w3)∈[0, 1] - 前記ステップ(d) での拡散フィルタw のためのディザリングアルゴリズムは、以下に示す通りであることを特徴とする請求項６に記載の方法。
但し、fRand はディザリングを微調整するためのパラメータであり、R(m,t)は注目ドットを走査するためのランダム値でのパラメータであり、R_MAX はランダムパラメータR(i)の最大値であり、cDither は、ディザリングの振幅を調整するためのパラメータであり、dw0 〜dw3 はディザリング後の異なる方向での拡散フィルタw の拡散重み付け係数である - 前記ステップ4)での動的階調可変出力機構は、
ShapeCurを考慮して注目ドットの偽似乱数値を生成することを含むことを特徴とする請求項１、２、３又は７のいずれかに記載の方法。
Fi=random(ShapeCur) 式１
但し、偽似乱数関数randomは、コンパイル環境で自動的に生成され、Fi∈[0,255] であり、
出力ドット・マトリクスデータは動的に算出される
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2005101166359A CN100377566C (zh) | 2005-10-26 | 2005-10-26 | 多位成像深度设备上的调频调幅混合半色调图像处理方法 |
| PCT/CN2006/000856 WO2007048289A1 (fr) | 2005-10-26 | 2006-04-29 | Methode de traitement d'une image en simili melange fm/am effectuee sur un appareil presentant une profondeur d'imagerie de plusieurs octets |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2009514273A true JP2009514273A (ja) | 2009-04-02 |
| JP4499175B2 JP4499175B2 (ja) | 2010-07-07 |
Family
ID=36760769
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2008536906A Expired - Fee Related JP4499175B2 (ja) | 2005-10-26 | 2006-04-29 | マルチビット深度の画像装置でｆｍ／ａｍ混在網点画像を処理する方法 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US8035860B2 (ja) |
| EP (1) | EP1942655A4 (ja) |
| JP (1) | JP4499175B2 (ja) |
| CN (1) | CN100377566C (ja) |
| WO (1) | WO2007048289A1 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2014533919A (ja) * | 2011-11-30 | 2014-12-15 | ペキン ユニバーシティ ファウンダー グループ カンパニー リミテッドPeking University Founder Group Co., Ltd. | マルチビット深度ハーフトーンａｍ網点を生成する方法と装置 |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100377567C (zh) * | 2005-10-26 | 2008-03-26 | 北京北大方正电子有限公司 | 在多位成像深度设备上的调频调幅混合网点网型控制方法 |
| CN100393099C (zh) * | 2006-08-31 | 2008-06-04 | 北京北大方正电子有限公司 | 调频调幅混合网点层次连续调控制方法 |
| CN101404714B (zh) * | 2008-10-17 | 2010-09-22 | 北大方正集团有限公司 | 一种图像处理方法及装置 |
| US9179039B2 (en) | 2014-02-12 | 2015-11-03 | Xerox Corporation | Methods and systems for processing low resolution images via error diffusion |
| CN107133936B (zh) * | 2017-05-09 | 2019-11-01 | 上海理工大学 | 数字半色调方法 |
| CN109884345A (zh) * | 2019-02-25 | 2019-06-14 | 燕山大学 | Afm轻敲模式下表征成像质量的方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004172798A (ja) * | 2002-11-19 | 2004-06-17 | Sony Corp | 画像処理装置、画像処理方法 |
| CN1668063A (zh) * | 2005-04-27 | 2005-09-14 | 北京北大方正电子有限公司 | 一种基于双反馈的误差扩散调频挂网方法 |
Family Cites Families (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5832122A (en) * | 1995-03-24 | 1998-11-03 | Fuji Photo Film Co., Ltd. | Method of processing image data |
| CN1062079C (zh) * | 1995-06-21 | 2001-02-14 | 时代集团公司 | 调幅调频挂网方法 |
| JPH11146189A (ja) | 1997-11-07 | 1999-05-28 | Dainippon Printing Co Ltd | デジタル階調画像の印刷方法および階調画像が表現された印刷物 |
| JP3884149B2 (ja) * | 1997-11-14 | 2007-02-21 | 株式会社アマダエンジニアリングセンター | 折曲げ加工装置における曲げ角度制御方法および装置 |
| JP3921678B2 (ja) | 1998-02-24 | 2007-05-30 | ソニー株式会社 | 画像処理方法および装置 |
| US6690837B1 (en) * | 1998-11-03 | 2004-02-10 | Agfa-Gevaert | Screening method for overlapping sub-images |
| US6271936B1 (en) * | 1998-12-11 | 2001-08-07 | Eastman Kodak Company | Combining error diffusion, dithering and over-modulation for smooth multilevel printing |
| US6433891B1 (en) * | 1998-12-14 | 2002-08-13 | Oak Technology, Inc. | Stochastic screening method with dot pattern regularity control and dot growth |
| JP2000287086A (ja) * | 1999-03-31 | 2000-10-13 | Minolta Co Ltd | 画像処理装置 |
| DE10022225B4 (de) * | 2000-05-08 | 2015-12-17 | Heidelberger Druckmaschinen Ag | Mehrstufen-Rasterung mit Rastern beliebiger Rasterweiten und Rasterwinkel |
| US20040130753A1 (en) * | 2003-01-06 | 2004-07-08 | Crounse Kenneth R. | Halftone method and system using hybrid AM/FM screening for highlight/shadow tonal regions |
| CN100481869C (zh) * | 2003-08-07 | 2009-04-22 | 宇东光学公司 | 自动检测图文的半色调处理方法 |
| TW200627897A (en) * | 2005-01-28 | 2006-08-01 | Univ Nat Taiwan Normal | Method for forming digital watermark with hybrid dots |
| CN100377567C (zh) * | 2005-10-26 | 2008-03-26 | 北京北大方正电子有限公司 | 在多位成像深度设备上的调频调幅混合网点网型控制方法 |
-
2005
- 2005-10-26 CN CNB2005101166359A patent/CN100377566C/zh not_active Expired - Fee Related
-
2006
- 2006-04-29 EP EP06722424A patent/EP1942655A4/en not_active Ceased
- 2006-04-29 WO PCT/CN2006/000856 patent/WO2007048289A1/zh active Application Filing
- 2006-04-29 US US12/091,290 patent/US8035860B2/en not_active Expired - Fee Related
- 2006-04-29 JP JP2008536906A patent/JP4499175B2/ja not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004172798A (ja) * | 2002-11-19 | 2004-06-17 | Sony Corp | 画像処理装置、画像処理方法 |
| CN1668063A (zh) * | 2005-04-27 | 2005-09-14 | 北京北大方正电子有限公司 | 一种基于双反馈的误差扩散调频挂网方法 |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2014533919A (ja) * | 2011-11-30 | 2014-12-15 | ペキン ユニバーシティ ファウンダー グループ カンパニー リミテッドPeking University Founder Group Co., Ltd. | マルチビット深度ハーフトーンａｍ網点を生成する方法と装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN1774029A (zh) | 2006-05-17 |
| EP1942655A4 (en) | 2008-12-03 |
| JP4499175B2 (ja) | 2010-07-07 |
| US20080278765A1 (en) | 2008-11-13 |
| EP1942655A1 (en) | 2008-07-09 |
| WO2007048289A1 (fr) | 2007-05-03 |
| CN100377566C (zh) | 2008-03-26 |
| US8035860B2 (en) | 2011-10-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4480782B2 (ja) | デュアルフィードバックに基づく誤差拡散を用いた周波数変調スクリーンのための方法 | |
| Li et al. | Tone-dependent error diffusion | |
| EP0838940B1 (en) | Multiple spatial channel processing | |
| JP4499175B2 (ja) | マルチビット深度の画像装置でｆｍ／ａｍ混在網点画像を処理する方法 | |
| JP4499176B2 (ja) | マルチビット深度の画像装置でｆｍ／ａｍ混在格子点及び格子点形状を制御する方法 | |
| JP4107059B2 (ja) | 画像処理方法および画像処理装置 | |
| US6335989B1 (en) | Halftone printing using donut filters | |
| WO2008028409A1 (fr) | Procédé de commande de l'ajustement continu d'une hiérarchie de point net combiné, à modulation d'amplitude/modulation de fréquence | |
| JP4149269B2 (ja) | 中間ドットを用いた誤差拡散方法及びシステム | |
| WO2006135810A2 (en) | Method and system for designing multilevel halftone screens using spot functions | |
| KR101098713B1 (ko) | 디더 마스크 생성 방법과 이를 포함하는 이미지 하프토닝 방법, 디더 마스크 생성 장치 및 컴퓨터 판독가능 매체 | |
| JPH0785272A (ja) | 周波数変調ハーフトーン画像および作成方法 | |
| US9167130B2 (en) | AM-FM halftone screen design | |
| US6697169B1 (en) | Enhanced error diffusion using a correlated chromatic noise pattern | |
| US8139267B2 (en) | Method, apparatus and program for creating a threshold value matrix | |
| US8867100B2 (en) | Image quantization for digital printing | |
| Mao et al. | 4-row serpentine tone dependent fast error diffusion | |
| Wu et al. | Supercell dither masks with constrained blue noise interpolation | |
| JP2018160883A (ja) | 画像処理装置、画像形成装置、画像処理方法及び画像処理プログラム | |
| CN114332253A (zh) | 基于图像阶调的半色调混合加网方法 | |
| JP2004112804A (ja) | サブドット位相変調を伴うｆｍスクリーニング | |
| He et al. | Clustered AM/FM Halftoning Algorithm | |
| JP2013514042A (ja) | 誤差拡散ドット密度の制御 | |
| JP5464906B2 (ja) | 画像処理装置及びその制御方法 | |
| Wang et al. | Research on Bilevel Rendition of Gray Pictures on Jacquard |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20091022 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20091027 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20100127 |
|
| A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20100203 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100216 |
|
| 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: 20100316 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20100414 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130423 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 Ref document number: 4499175 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130423 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140423 Year of fee payment: 4 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| LAPS | Cancellation because of no payment of annual fees |