JP3574711B2 - Image data binarization method and apparatus - Google Patents
Image data binarization method and apparatus Download PDFInfo
- Publication number
- JP3574711B2 JP3574711B2 JP09623096A JP9623096A JP3574711B2 JP 3574711 B2 JP3574711 B2 JP 3574711B2 JP 09623096 A JP09623096 A JP 09623096A JP 9623096 A JP9623096 A JP 9623096A JP 3574711 B2 JP3574711 B2 JP 3574711B2
- Authority
- JP
- Japan
- Prior art keywords
- image data
- pixel
- error
- binarization
- target pixel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
Description
【0001】
【発明の属する技術分野】
本発明は、多値の画像データを2値の画像データに変換する画像データ2値化方法、及び装置に関するものである。
【0002】
【従来の技術】
従来、インクジェットプリンタ等の2値記録装置によって、疑似的に多階調画像を記録するに際しては、多値の画像データを2値化するために誤差拡散方法が多用されていた。
【0003】
図14はそのような誤差拡散方法を採用した従来の画像データ2値化装置を示すブロック図である。この画像データ2値化装置は、図14に示すように、誤差メモリ120、加算器122、比較器124、減算器126、ビット変換器128、及び誤差分配回路130を備えている。
【0004】
画像信号のダイナミックレンジは0〜255とし、図中太線部は符号ビット、整数部ビット、小数部ビットによって成る多ビットのデータ線であることを示している。
【0005】
ここで、誤差拡散方法とは、或る画素について画像データを2値化した際に、2値化前の画像データと2値化後の画像データとの誤差を、所定の割合で上記画素の近傍の画素にそれぞれ分配して拡散する方法である。図15はこのような誤差拡散方法の原理を説明するための説明図である。図15において、230は画像を示しており、231の矢印は画像230に対して2値化処理する際の処理順序を示しており、また、232は画像230中の画素を示している。
【0006】
2値化対象画素(以下注目画素という)Pについての2値化前の画像データの値をXとし、2値化済みの画素から注目画素Pに分配された誤差の累計をΣeP とすれば、注目画素に対する補正データX’は次式により表される。
X’=X+ΣeP ………(1)
【0007】
注目画素Pに対する補正データX’を閾値(例えば128)と比較することにより注目画素Pの2値化が行われるが、2値化後の画像データがY(0又は255)であれば、注目画素Pの2値化に伴ってこの画素から発生する誤差eは次式で表される。
e=X’−Y ………(2)
【0008】
この誤差eの値を注目画素Pの近傍の画素にそれぞれ分配して拡散する。即ち、誤差eの値を所定の割合(誤差分配率)で分割して、近傍の画素についての2値化前の画像データにそれぞれ加算する。従来の拡散処理の例としては、誤差eの値を注目画素Pに対して右隣の画素A、右下の画素B、真下の画素C、及び左下の画素Dの4つの画素へ分配し、誤差分配率として、画素Aに対し7/16、画素Bに対し1/16、画素Cに対し5/16、画素Dに対し3/16を適用する例が知られている。
【0009】
画素A:eA=e×7/16
画素B:eB=e×1/16 ………(3)
画素C:eC=e×5/16
画素D:eD=e×3/16
【0010】
このようして、近傍の画素に拡散された誤差は、画素毎に、誤差データとして誤差メモリ120に蓄積されることになる。
【0011】
では、図14に示す画像データ2値化装置について説明する。2値化対象の注目画素がPである場合、加算器122には画素Pについての多値の画像データが入力される。また、誤差メモリ120からは、注目画素Pに蓄積された誤差データΣePが読み出され、加算器122に入力される。加算器122では、入力された多値の画像データに誤差データを加算し注目画素の画像データを補正する。
【0012】
比較器124では、加算器122からの補正済画像データを入力し、別に入力される基準値と比較して、それら値の大小関係に応じて1ビットのデータ、即ち、2値の画像データを出力する。これにより、加算器122から出力された画像データは2値化される。また、ビット変換器128では、比較器124からの1ビットの2値の画像データを、8ビットの2値の画像データ(0又は255)に変換する。
【0013】
減算器126では、加算器122からの補正済画像データ(2値化前の画像データ)とビット変換器128からの画像データ(2値化後の画像データ)とを入力し、2値化前の画像データから2値化後の画像データを減算して、注目画素Pの2値化に伴って発生する誤差eを導き出す。
【0014】
誤差分配回路130では、減算器126によって得られた誤差eを、図15において述べたように注目画素Pの近傍の画素A,B,C,Dにそれぞれ分配して拡散する。即ち、誤差分配回路130では、例えば式(2)で示したように、誤差eに対して画素A,B,C,Dに対応する係数をそれぞれ掛けた後、各演算結果eA,eB,eC,eDを、誤差メモリ120内の画素A,B,C,Dについての誤差データΣeA,ΣeB,ΣeC,ΣeDにそれぞれ加算して、各誤差データΣeA,ΣeB,ΣeC,ΣeDの値を更新する。
【0015】
以上のように、図14に示す画像データ2値化装置によれば、多値の画像データを2値化する際に誤差拡散方法を利用することによって、2値化しても、画像を疑似的に多階調にて表現することが可能となる。
【0016】
【発明が解決しようとする課題】
しかしながら、上記した従来の画像データ2値化装置においては、比較的低い階調値の画像データが入力された場合、2値化された画像データとして、値“1”となる画素が画像中において点線状に配置され、また、比較的高い階調値の画像データが入力された場合には、2値化された画像データとして、値“0”となる画素が画像中において点線状に配置されてしまい、著しく画質を損ねてしまうという問題があった。
【0017】
今、図14に示す画像データ2値化装置をプリンタ等に採用した場合、画像データの階調値は記録画像における濃度に対応することになる。従って、比較的低い階調値の画像データが入力される場合は、記録結果として比較的低い濃度による記録画像が得られ、比較的高い階調値の画像データが入力される場合には、比較的高い濃度による記録画像が得られる。
【0018】
図16及び図17はそれぞれ図14の画像データ2値化装置を用いた場合のプリンタによる記録結果を示す説明図である。図16は低濃度領域(0%〜15%)についての記録結果を示し、図17は高濃度領域(85%〜100%)についての記録結果を示している。図16及び図17では、画像の横方向において、濃度が上記した範囲で連続的に変化している。なお、2値化処理の順序は、図15に示したように、画像の左上隅を始点として左から右へ、左から右へと横方向に1ラインずつ走査することにより行なわれる。
上記のことは後述の2値化処理においても同様である。
【0019】
また、図16に示した低濃度領域のうち、特に濃度1%についての記録結果を図18に示し、図17に示した高濃度領域のうち、特に濃度99%についての記録結果を図19に示す。即ち、図18及び図19においては、画像全体を均一の濃度としている。図20は図15の画像データ2値化装置を用いた場合のプリンタによる中間濃度領域(15%〜85%)についての記録結果を示す説明図である。図20では画像の横方向において濃度が上記の範囲で連続的に変化している。
【0020】
図16または図18に示すように、0%〜3%の比較的低い濃度においては、部分的に、黒点(即ち、値“1”の画素)が右上から左下に向かって点線状に連続して配置されている。また、図17または図19に示すように、97%〜100%の比較的高い濃度においては、部分的に、白点(即ち、値“0”の画素)が右上から左下に向かって点線状に連続して配置されている。従って、何れの場合も、黒点または白点が均等に散在しておらず、そのため、著しく画質を損ねてしまっている。
【0021】
図20に示される中間濃度領域(15%〜85%)においては、一部に顕著なテクスチャが認められるが、このようなテクスチャは、2値化手段における閾値に乱数を重畳させることによって除去可能である。
【0022】
一方、低濃度領域あるいは高濃度領域に発生するテクスチャは、閾値に乱数を重畳させても除去することはできない。
【0023】
このような低濃度領域あるいは高濃度領域に発生するテクスチャは、注目画素の2値化に伴って発生する誤差をより多くの2値化前の画素に分配することにより改善されるが、このような方法はハードウェア規模が大きくなったり、あるいは2値化処理時間が長くなる等の問題を有していた。
【0024】
本発明は、上記した従来技術の問題点に鑑みなされたものであり、その目的は、より改善された、画像データの2値化を行なうことができる画像データ2値化装置を提供することにある。
【0025】
【課題を解決するための手段およびその作用・効果】
上記した目的の少なくとも一部を達成するために、
本発明の請求項1に記載された方法は、
多値の画像データを2値化する画像データ2値化方法であって、
(a)注目画素の画像データに、前記注目画素に蓄積された誤差と、前記注目画素の近傍に位置する2値化前の非注目画素に蓄積された誤差とを加算し、該加算後の画像データを前記注目画素の修正データとする工程と、
(b)前記注目画素および前記非注目画素に蓄積された誤差を更新する工程と、
(c)前記補正データと閾値に基づき前記注目画素の画像データを2値化する工程と、
(d)前記補正データと前記2値化の結果に基づき前記注目画素の2値化に伴って発生する誤差を導き出す工程と、
(e)前記導き出された誤差を2値化前の画素に分配する工程と、
(f)分配された前記誤差を画素毎に誤差データとして蓄積する工程と、
によって画像データを2値化することを要旨とする。
【0026】
請求項1記載の方法によれば、注目画素の画像データに、2値化済みの画素から注目画素に蓄積された誤差と、注目画素の近傍に位置する2値化前の非注目画素に蓄積された誤差とを加算することによって、注目画素の画像データが補正される。
【0027】
前記非注目画素に蓄積された誤差とは、注目画素の近傍に位置する2値化済みの画素から発生した誤差が分配され蓄積したものである。
【0028】
従って、注目画素は注目画素近傍に位置する数多の2値化済み画素から発生した誤差に基づき補正されることになり、このため、2値化後の画像品質が著しく改善される。
【0029】
前記請求項2記載の方法によれば、非注目画素に蓄積された誤差を注目画素の画像データに加算する代わりに、前記非注目画素に注目画素の画像データに応じて決定される係数を乗じた値を注目画素の画像データに加算して補正データを得ている。
【0030】
従って、注目画素の画像データにどの程度、2値化前の画素に蓄積された誤差を加算するかを注目画素の画像データに応じて決定することが可能であり、2値化後の画像品質の更なる改善あるいは2値化処理時間の短縮等の効果が得られる。
【0031】
請求項3に記載の発明は、
多値の画像データを2値化する画像データ2値化装置であって、
注目画素の画像データに、前記注目画素に蓄積された誤差と、前記注目画素の近傍に位置する2値化前の非注目画素に蓄積された誤差とを加算することにより前記注目画素に対する補正データを導出する補正データ導出手段と、
前記注目画素および前記非注目画素に蓄積された誤差を更新する更新手段と、
前記補正データと閾値に基づき前記注目画素の画像データを2値化する2値化手段と、
前記補正データと前記2値化の結果に基づき前記注目画素の2値化に伴って発生する誤差を2値化前の画素に分配する誤差分配手段と、
前記分配された誤差を画素毎に誤差データとして蓄積する誤差蓄積手段とを備えたことを要旨としている。
【0032】
請求項3に記載の発明によっても、注目画素は注目画素近傍に位置する数多の2値化済み画素から発生した誤差に基づき補正されることになり、このため、2値化後の画像品質が著しく改善する。
【0033】
【発明の実施の形態】
以下、本発明の実施の形態を実施例に基づいて説明する。図1は本発明の第1の実施例における誤差拡散の方法を説明するための説明図である。図1(a)において、220は画像を示しており、221の矢印は画像220に対する2値化処理の順序を示しており、また、222は画像220中の2値化対象の注目画素を示している。
【0034】
図1(b)は本実施例における2値化対象の注目画素Pに対する補正データを導出する方法を説明するための説明図である。
【0035】
本発明においては、注目画素の多値データは注目画素に蓄積された誤差ΣeP と注目画素近傍に位置する2値化前の画素に蓄積された誤差を注目画素に加算することにより、注目画素Pに対する補正データを導出する。
【0036】
本実施例においては、注目画素近傍に位置する2値化前の画素として画素Eを選び画素Eに蓄積された誤差ΣeE の1/8を注目画素の画像データに加算すると共に、注目画素に蓄積された誤差はその全部を注目画素の画像データに加算する。
【0037】
注目画素Pについての2値化前の画像データをXとし、2値化済みの画素から注目画素Pに分配された誤差の累計をΣeP とすれば、注目画素に対する補正データX’’は次式に従って導出される。
X’’=X+ΣeP+1/8ΣeE ………(4)
【0038】
注目画素Pの画像データの補正が完了すれば、誤差メモリ20内の画素P,Eに対する誤差値は次式に従って更新される。
ΣeP →0
ΣeE →7/8ΣeE ………(5)
【0039】
注目画素Pに対する補正データX’’に閾値(例えば128)と比較することにより注目画素Pの2値化が行われるが、2値化後の画像データがY(0又は255)であれば、注目画素Pの2値化に伴ってこの画素から発生する誤差eは次式で表される。
e=X’’−Y ………(6)
【0040】
導出された誤差は図1(c)に示されるように注目画素近傍の2値化前の画素A,B,C,Dに所定の比率で配分される。
【0041】
このようにして近傍画素に配分された誤差は、画素毎に誤差データとして誤差メモリに蓄積される。
【0042】
以上で画素Pに対する2値化処理は完了し、次画素Aに対する2値化処理が開始する。
【0043】
図2は本発明の第1の実施例としての画像データ2値化装置を示すブロック図である。この画像データ2値化装置は、図2に示すように、誤差メモリ20、データ補正回路22、比較器24、減算器26、ビット変換器28、誤差分配回路30、誤差更新回路32を備えている。
【0044】
誤差メモリ20では、図2に示すように画像の2ライン分しか記憶容量がないが、注目画素に対する補正が完了すれば、誤差メモリ20における注目画素に対する誤差を0に更新することにより、この2ライン分の記憶領域をトグルで利用することが可能となり、1枚の画像の各画素についての誤差データを順番に蓄積することができる。
【0045】
図2において、2値化対象の注目画素がPである場合、データ補正回路22には画素Pについての多値の画像データが入力される。この多値の画像データは例えば8ビットのデータであり、データの値(即ち、諧調値)としては“0”〜“255”の値を採り得る。
【0046】
誤差メモリ20からは、注目画素Pに対して、2値化済みの近傍画素から分配され蓄積された誤差データΣeP と、画素Eに対して、2値化済みの近傍画素から分配され蓄積された誤差データΣeE が読み出され、データ補正回路22に入力される。
【0047】
なお、誤差データΣeP と誤差データΣeE とが読み出されると、誤差更新回路32によって誤差メモリ20内に蓄積されていた誤差データΣeP は0に、誤差データΣeE は7/8ΣeE に更新される。
【0048】
データ補正回路22では、注目画素Pの画像データに、注目画素Pに蓄積された誤差ΣeP と画素Eに蓄積された誤差ΣeE の一部、即ち、1/8ΣeE を加算し、注目画素の画像データを補正する。
【0049】
比較器24では、データ補正回路22からの補正済画像データを入力し、別に入力される基準値と比較して、それら値の大小関係に応じて1ビットのデータ、即ち、2値の画像データを出力する。ここで、基準値としては“0”〜“255”の中間の値である“128”を採用する。即ち、比較器24は、データ補正回路22からの画像データの値が“128”以上の場合には“1”を出力し、“128”よりも小さい場合には“0”を出力する。こうして、データ補正回路22より出力された画像データは2値化される。
【0050】
ビット変換器28では、比較器24からの1ビットの2値の画像データを入力し、8ビットの2値の画像データに変換して出力する。即ち、比較器24からの画像データの値が“0”の場合は、8ビットの画像データとして値“0”のデータを出力し、比較器24からの画像データの値が“1”の場合は、8ビットの画像データとして値“255”のデータを出力する。なお、このようなビット変換器28は、例えば、比較器24の出力から分岐された1ビットの画像データ線を8本に分岐して、8ビットの画像データ線として減算器26の入力に接続することによって容易に実現することができる。
【0051】
減算器26では、データ補正回路22からの画像データ(2値化前の画像データ)とビット変換器28からの画像データ(2値化後の画像データ)とを入力し、2値化前の画像データから2値化後の画像データを減算して、注目画素の2値化によって発生した誤差eを導き出す。
【0052】
導出された誤差eは、誤差分配回路30において、式(3)に従ってeA ,eB ,eC ,eD と分割され、注目画素の近傍に位置する画素A,B,C,Dに分配される。即ち、誤差メモリ20内の画素A,B,C,Dに対する誤差値ΣeA ,ΣeB ,ΣeC ,ΣeD に、eA ,eB ,eC ,eD を加算することにより、ΣeA ,ΣeB ,ΣeC ,ΣeD の値を更新する。
【0053】
このようにして、注目画素Pについて一連の処理が終了したら、2値化処理が画素Pから次の画素Aに移り、今度は、その画素Aを2値化対象の注目画素として同様の処理が繰り返される。
【0054】
図3は本実施例による低濃度領域(0%〜15%)についてのプリンタによる記録結果を示す説明図であり、従来例の図16と対応するものである。
【0055】
図4は本実施例による高濃度領域(85%〜100%)についてのプリンタによる記録結果を示す説明図であり、従来例の図17と対応するものである。
【0056】
図5は本実施例による低濃度領域(1%一定)についてのプリンタによる記録結果を示す説明図であり、従来例の図18と対応するものである。
【0057】
図6は本実施例による高濃度領域(99%一定)についてのプリンタによる記録結果を示す説明図であり、従来例の図19と対応するものである。
何れも黒点、白点は均等に配置され画質は著しく改善されている。
【0058】
図7は本実施例による中間濃度領域(15%〜85%)についてのプリンタによる記録結果を示す説明図であり、従来例の図20と対応するものである。
この領域においては、従来方法と比較して顕著な差異は認められない。
【0059】
前記実施例においては、注目画素の画像データを補正するにあたっては、注目画素の画像データの値とは無関係に画素Eに蓄積された誤差ΣeE に対する一定比率(1/8)の誤差量を注目画素の画像データに分配しているが、この分配率を注目画素の画像データに対応した好ましい値に設定することにより、2値化後の画像品質の更なる改善、あるいは2値化処理時間の短縮等の効果が得られる。
【0060】
図8は本発明の第2の実施例における2値化対象の注目画素Pに対する補正データを導出する方法を説明するための説明図である。
【0061】
本実施例においては、注目画素Pに蓄積された誤差ΣeP 、画素Eに蓄積された誤差ΣeE 、及び画素Fに蓄積された誤差ΣeF を注目画素の画像データに分配することにより、注目画素の補正データを導出する。画素E,Fに蓄積された誤差ΣeE ,ΣeF の画素Pの画像データに対する分配率は一定でなく、図8(b),(c)に示されるように、画素Pの画像データの値に応じて、それぞれWE ,WF と決定される。
【0062】
このようにして決定された分配率WE ,WF に基づき、注目画素Pの画像データXに対する補正データX’’’ は次式に従って導出される。
X’’’=X+ΣeP+WE・ΣeE+WFΣeF ………(7)
【0063】
このようにして注目画素Pの画像データの補正が完了すれば、誤差メモリ内の画素P,E,Fに対する誤差値は次式に従って更新される。
ΣeP →0
ΣeE →(1−WE )・ΣeE ………(8)
ΣeF →(1−WF )・ΣeF
【0064】
以降の動作は前記の第1の実施例と全く同じである。
【0065】
図9は本発明の第2の実施例による画像データ2値化装置のブロック図である。
データ補正回路22には画素Pについての多値画像データXが入力される。データ補正回路22は、先ず、誤差メモリ20からΣeP ,ΣeE ,ΣeF を読み出し、次に、注目画素Pの画像データXの値に対応する分配率WE ,WF を決定し、式(7)に従って注目画素Pの画像データに対する補正データX’’’ を導出する。
【0066】
注目画素Pの画像データXに対応する分配率WE ,WF を決定する方法としては、例えば、画像データXをメモリにアドレスとして入力し、このアドレスに対応する記憶内容としてWE ,WF を格納しておけばよい。あるいは、画素Pの画像データXと誤差値ΣeE ,ΣeF をメモリにアドレスとして入力し、このアドレスに対応する記憶内容として画像データX+WE・ΣeE +WF・ΣeF を格納しておけば、補正データX’’’ は迅速に導出される。
【0067】
このようにして注目画素Pの画像データの補正が完了すれば、誤差メモリ20内の誤差データΣeP ,ΣeE ,ΣeF は、式(8)に従って誤差更新回路32によって更新される。
【0068】
更新データ(1−WE )・ΣeE ,(1−WF )・ΣeF を導出する方法としては、例えば、注目画素Pの画像データXと誤差値ΣeE 、あるいは注目画素Pの画像データXと誤差値ΣeF をアドレスとしてメモリに入力し、このアドレスに対応する記憶内容として、(1−WE)・ΣeE あるいは(1−WF )・ΣeF
を格納しておけばよい。
以降の動作は第1の実施例と全く同じである。
【0069】
図10は本実施例による低濃度領域(15%〜0%)についてのプリンタによる記録結果を示す説明図であり、従来例の図16と対応するものである。
【0070】
図11は本実施例による高濃度領域(85%〜100%)についてのプリンタによる記録結果を示す説明図であり、従来例の図17と対応するものである。
【0071】
図12は本実施例による低濃度領域(1%一定)についてのプリンタによる記録結果を示す説明図であり、従来例の図18と対応するものである。
【0072】
図13は本実施例による高濃度領域(99%一定)についてのプリンタによる記録結果を示す説明図であり、従来例の図19と対応するものである。
何れも黒点、白点は均等に配置され画質は著しく改善されている。
【0073】
本実施例によれば、低濃度領域(1%)及び高濃度領域(99%)において、第1の実施例と比べ若干の改善が認められる。
【0074】
また、本実施例においては、濃度領域10%〜90%の部分に対する処理は従来方法と何等変わらないので、2値化処理に要する時間も従来方法とさして変わらない。
【図面の簡単な説明】
【図1】本発明の第1の実施例における誤差拡散の方法を説明するための説明図である。
【図2】本発明の第1の実施例としての画像データ2値化装置を示すブロック図である。
【図3】本発明の第1の実施例による低濃度領域(0%〜15%)についての記録結果を示す説明図である。
【図4】本発明の第1の実施例による高濃度領域(85%〜100%)についての記録結果を示す説明図である。
【図5】本発明の第1の実施例による濃度1%についての記録結果を示す説明図である。
【図6】本発明の第1の実施例による濃度99%についての記録結果を示す説明図である。
【図7】本発明の第1の実施例による中間濃度領域(15%〜85%)についての記録結果を示す説明図である。
【図8】本発明の第2の実施例における誤差拡散の方法を説明するための説明図である。
【図9】本発明の第2の実施例としての画像データ2値化装置を示すブロック図である。
【図10】本発明の第2の実施例による低濃度領域(0%〜15%)についての記録結果を示す説明図である。
【図11】本発明の第2の実施例による高濃度領域(85%〜100%)についての記録結果を示す説明図である。
【図12】本発明の第2の実施例による濃度1%についての記録結果を示す説明図である。
【図13】本発明の第2の実施例による濃度99%についての記録結果を示す説明図である。
【図14】従来の画像データ2値化装置を示すブロック図である。
【図15】誤差拡散方法の原理を説明するための説明図である。
【図16】従来の画像データ2値化装置を用いた場合のプリンタによる低濃度領域(0%〜15%)についての記録結果を示す説明図である。
【図17】従来の画像データ2値化装置を用いた場合のプリンタによる高濃度領域(85%〜100%)についての記録結果を示す説明図である。
【図18】従来の画像データ2値化装置を用いた場合のプリンタによる濃度1%についての記録結果を示す説明図である。
【図19】従来の画像データ2値化装置を用いた場合のプリンタによる濃度99%についての記録結果を示す説明図である。
【図20】従来の画像データ2値化装置を用いた場合のプリンタによる中間濃度領域(15%〜85%)についての記録結果を示す説明図である。
【符号の説明】
20…誤差メモリ
22…データ補正回路
24…比較器
26…減算器
28…ビット変換器
30…誤差分配回路
32…誤差更新回路
120…誤差メモリ
122…加算器
124…比較器
126…減算器
128…ビット変換器
130…誤差分配回路
220,230…画像
221,231…走査線
222,232…画素
A,B,C,D,E,F…画素
P…注目画素
P100〜P111…画素
P200〜P211…画素
eA,eB,eC,eD,eE,eF…分配誤差
e…誤差
ΣeA,ΣeB,ΣeC,ΣeD,ΣeE,ΣeF…誤差データ
ΣeP…誤差データ[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention relates to an image data binarization method and apparatus for converting multi-valued image data into binary image data.
[0002]
[Prior art]
2. Description of the Related Art Conventionally, when a multi-level image is pseudo-recorded by a binary recording apparatus such as an ink jet printer, an error diffusion method has been frequently used to binarize multi-level image data.
[0003]
FIG. 14 is a block diagram showing a conventional image data binarization device employing such an error diffusion method. As shown in FIG. 14, this image data binarization device includes an
[0004]
The dynamic range of the image signal is 0 to 255, and the bold line in the figure indicates a multi-bit data line including a sign bit, an integer part bit, and a decimal part bit.
[0005]
Here, the error diffusion method means that when image data is binarized for a certain pixel, an error between the image data before binarization and the image data after binarization is determined by a predetermined ratio of the pixel. This is a method of distributing and diffusing each pixel to neighboring pixels. FIG. 15 is an explanatory diagram for explaining the principle of such an error diffusion method. In FIG. 15,
[0006]
If the value of the image data before binarization for the binarization target pixel (hereinafter referred to as the pixel of interest) P is X, and the total error distributed from the binarized pixel to the pixel of interest P is Σe P The correction data X ′ for the target pixel is expressed by the following equation.
X ′ = X + Σe P (1)
[0007]
The target pixel P is binarized by comparing the correction data X ′ for the target pixel P with a threshold value (for example, 128). If the binarized image data is Y (0 or 255), the binarization is performed. The error e generated from the pixel P due to the binarization of the pixel P is expressed by the following equation.
e = X'-Y (2)
[0008]
The value of the error e is distributed and diffused to pixels near the target pixel P. That is, the value of the error e is divided at a predetermined ratio (error distribution ratio) and added to the image data of the neighboring pixels before binarization. As an example of the conventional diffusion process, the value of the error e is distributed to four pixels: a pixel A on the right, a pixel B on the lower right, a pixel C on the lower right, and a pixel D on the lower left with respect to the pixel of interest P; As an error distribution ratio, an example in which 7/16 is applied to the pixel A, 1/16 to the pixel B, 5/16 to the pixel C, and 3/16 to the pixel D is known.
[0009]
Pixel A: e A = e × 7/16
Pixel B: e B = e × 1/16 (3)
Pixel C: e C = e × 5/16
Pixel D: e D = e × 3/16
[0010]
In this way, the error diffused to the neighboring pixels is accumulated in the
[0011]
Now, the image data binarization device shown in FIG. 14 will be described. When the target pixel to be binarized is P, multivalued image data for the pixel P is input to the
[0012]
The
[0013]
The
[0014]
The
[0015]
As described above, according to the image data binarization device shown in FIG. 14, when the multi-valued image data is binarized, the error diffusion method is used, so that even if the image is binarized, the image is pseudo-generated. Can be expressed in multiple gradations.
[0016]
[Problems to be solved by the invention]
However, in the above-described conventional image data binarization device, when image data having a relatively low gradation value is input, a pixel having a value “1” is binarized image data in the image. When image data having a relatively high gradation value is input in the form of a dotted line, pixels having a value “0” are arranged in the image in a dotted line as binarized image data. There is a problem that image quality is significantly impaired.
[0017]
Now, when the image data binarization device shown in FIG. 14 is employed in a printer or the like, the gradation value of the image data corresponds to the density in the recorded image. Therefore, when image data with a relatively low tone value is input, a print image with a relatively low density is obtained as a printing result, and when image data with a relatively high tone value is input, a comparative image is obtained. A recorded image with an extremely high density can be obtained.
[0018]
FIGS. 16 and 17 are explanatory views showing recording results by the printer when the image data binarization device of FIG. 14 is used. FIG. 16 shows a recording result for a low density area (0% to 15%), and FIG. 17 shows a recording result for a high density area (85% to 100%). 16 and 17, in the horizontal direction of the image, the density continuously changes in the above-described range. As shown in FIG. 15, the binarization process is performed by scanning one line horizontally from left to right and from left to right starting from the upper left corner of the image.
The same applies to the binarization processing described later.
[0019]
Also, FIG. 18 shows a recording result of the low-density area shown in FIG. 16, particularly for the density of 1%, and FIG. 19 shows a recording result of the high-density area shown in FIG. 17, particularly for the density of 99%. Show. That is, in FIGS. 18 and 19, the entire image has a uniform density. FIG. 20 is an explanatory diagram showing a printing result of an intermediate density area (15% to 85%) by a printer when the image data binarization apparatus of FIG. 15 is used. In FIG. 20, the density continuously changes in the above range in the horizontal direction of the image.
[0020]
As shown in FIG. 16 or FIG. 18, at a relatively low density of 0% to 3%, a black point (that is, a pixel having a value “1”) partially continues in a dotted line from upper right to lower left. Is arranged. Also, as shown in FIG. 17 or FIG. 19, at a relatively high density of 97% to 100%, a white point (that is, a pixel having a value “0”) partially becomes a dotted line from upper right to lower left. Are arranged continuously. Therefore, in each case, the black points or white points are not evenly scattered, thereby significantly deteriorating the image quality.
[0021]
In the intermediate density region (15% to 85%) shown in FIG. 20, a remarkable texture is partially observed, but such a texture can be removed by superimposing a random number on a threshold value in the binarization unit. It is.
[0022]
On the other hand, a texture generated in a low-density area or a high-density area cannot be removed even if a random number is superimposed on a threshold.
[0023]
The texture generated in such a low-density region or a high-density region is improved by distributing the error generated due to the binarization of the target pixel to more pixels before binarization. These methods have problems such as an increase in hardware scale and an increase in binarization processing time.
[0024]
SUMMARY OF THE INVENTION The present invention has been made in consideration of the above-described problems of the related art, and has as its object to provide an improved image data binarization device capable of performing binarization of image data. is there.
[0025]
[Means for Solving the Problems and Their Functions and Effects]
To achieve at least some of the above objectives,
The method according to
An image data binarization method for binarizing multi-valued image data,
(A) adding, to the image data of the target pixel, the error accumulated in the target pixel and the error accumulated in the non-target pixel before binarization located in the vicinity of the target pixel; Making image data correction data of the pixel of interest;
(B) updating the error accumulated in the target pixel and the non-target pixel;
(C) binarizing the image data of the pixel of interest based on the correction data and a threshold;
(D) deriving an error generated with the binarization of the target pixel based on the correction data and the binarization result;
(E) distributing the derived error to pixels before binarization;
(F) accumulating the distributed error as error data for each pixel;
The gist is to binarize the image data.
[0026]
According to the method of
[0027]
The error accumulated in the non-target pixel is an error generated and distributed from a binarized pixel located near the target pixel.
[0028]
Therefore, the target pixel is corrected based on errors generated from a number of binarized pixels located in the vicinity of the target pixel, so that the image quality after binarization is significantly improved.
[0029]
According to the method of claim 2, instead of adding the error accumulated in the non-target pixel to the image data of the target pixel, the non-target pixel is multiplied by a coefficient determined according to the image data of the target pixel. The obtained value is added to the image data of the target pixel to obtain correction data.
[0030]
Therefore, it is possible to determine to what extent the error accumulated in the pixel before binarization is added to the image data of the target pixel in accordance with the image data of the target pixel. And the effect of shortening the binarization processing time can be obtained.
[0031]
The invention according to
An image data binarization device for binarizing multi-valued image data,
Correction data for the pixel of interest is obtained by adding the error accumulated in the pixel of interest to the image data of the pixel of interest and the error accumulated in the non-target pixel before binarization located in the vicinity of the pixel of interest. Correction data deriving means for deriving
Updating means for updating the error accumulated in the target pixel and the non-target pixel,
Binarizing means for binarizing the image data of the pixel of interest based on the correction data and a threshold value;
Error distributing means for distributing an error generated due to binarization of the target pixel based on the correction data and the binarization result to pixels before binarization;
And an error accumulating means for accumulating the distributed error as error data for each pixel.
[0032]
According to the third aspect of the present invention, the target pixel is corrected based on errors generated from a large number of binarized pixels located near the target pixel. Therefore, the image quality after the binarization is corrected. Is significantly improved.
[0033]
BEST MODE FOR CARRYING OUT THE INVENTION
Hereinafter, embodiments of the present invention will be described based on examples. FIG. 1 is an explanatory diagram for explaining an error diffusion method according to the first embodiment of the present invention. In FIG. 1A,
[0034]
FIG. 1B is an explanatory diagram for describing a method for deriving correction data for the target pixel P to be binarized in the present embodiment.
[0035]
In the present invention, the multi-value data of the target pixel is obtained by adding the error Δe P accumulated in the target pixel and the error accumulated in the pixel before binarization located near the target pixel to the target pixel. The correction data for P is derived.
[0036]
In the present embodiment, the pixel E is selected as a pre-binarized pixel located in the vicinity of the pixel of interest, and 1/8 of the error Δe E accumulated in the pixel E is added to the image data of the pixel of interest. All of the accumulated errors are added to the image data of the target pixel.
[0037]
The image data before binarization for the pixel of interest P and X, if the total of distributed errors from binarized pixels to the pixel of interest P and Sigma] e P, the correction data X with respect to the pixel of interest '' the following It is derived according to the formula.
X '' = X + Σe P + / Σe E (4)
[0038]
When the correction of the image data of the target pixel P is completed, the error value for the pixels P and E in the
Σe P → 0
Σe E → 7 / 8Σe E ……… (5)
[0039]
The binarization of the pixel of interest P is performed by comparing the correction data X ″ for the pixel of interest P with a threshold value (eg, 128). If the image data after binarization is Y (0 or 255), An error e generated from the pixel of interest P due to binarization of the pixel is expressed by the following equation.
e = X ''-Y (6)
[0040]
The derived error is distributed at a predetermined ratio to the pixels A, B, C, and D before binarization in the vicinity of the target pixel as shown in FIG.
[0041]
The error distributed to the neighboring pixels in this manner is stored in the error memory as error data for each pixel.
[0042]
Thus, the binarization process for the pixel P is completed, and the binarization process for the next pixel A starts.
[0043]
FIG. 2 is a block diagram showing an image data binarization device as a first embodiment of the present invention. As shown in FIG. 2, the image data binarization device includes an
[0044]
The
[0045]
In FIG. 2, when the target pixel to be binarized is P, multivalued image data of the pixel P is input to the
[0046]
The
[0047]
Incidentally, when the error data Sigma] e P and error data Sigma] e E is read, the error data Sigma]
[0048]
The
[0049]
The
[0050]
The
[0051]
The
[0052]
The derived error e is divided by the
[0053]
In this way, when a series of processing is completed for the target pixel P, the binarization processing shifts from the pixel P to the next pixel A, and the same processing is performed using the pixel A as the target pixel for binarization. Repeated.
[0054]
FIG. 3 is an explanatory diagram showing a printing result by a printer in a low density area (0% to 15%) according to the present embodiment, and corresponds to FIG. 16 of the conventional example.
[0055]
FIG. 4 is an explanatory diagram showing a printing result of a high density area (85% to 100%) by the printer according to the present embodiment, and corresponds to FIG. 17 of the conventional example.
[0056]
FIG. 5 is an explanatory diagram showing a printing result of a low density area (1% constant) by a printer according to the present embodiment, and corresponds to FIG. 18 of the conventional example.
[0057]
FIG. 6 is an explanatory diagram showing a printing result of a high density area (99% constant) by the printer according to the present embodiment, and corresponds to FIG. 19 of the conventional example.
In each case, the black point and the white point are arranged evenly, and the image quality is remarkably improved.
[0058]
FIG. 7 is an explanatory diagram showing a printing result of the intermediate density area (15% to 85%) by the printer according to the present embodiment, and corresponds to FIG. 20 of the conventional example.
In this region, no remarkable difference is observed as compared with the conventional method.
[0059]
In the above-described embodiment, when correcting the image data of the target pixel, the error amount of a fixed ratio (に 対 す る) to the error Δe E accumulated in the pixel E regardless of the value of the image data of the target pixel. Although the image data is distributed to the pixel image data, by setting this distribution ratio to a preferable value corresponding to the image data of the target pixel, the image quality after the binarization is further improved, or the binarization processing time is reduced. Effects such as shortening can be obtained.
[0060]
FIG. 8 is an explanatory diagram for explaining a method for deriving correction data for a target pixel P to be binarized according to the second embodiment of the present invention.
[0061]
In this embodiment, the error Δe P accumulated in the pixel of interest P , the error Δe E accumulated in the pixel E , and the error Δe F accumulated in the pixel F are distributed to the image data of the pixel of interest. The correction data of the pixel is derived. The distribution ratios of the errors Δe E and Δe F accumulated in the pixels E and F to the image data of the pixel P are not constant, and as shown in FIGS. 8B and 8C, the values of the image data of the pixel P depending on, it is respectively W E, and W F determined.
[0062]
Thus it is determined the distribution ratio W E, based on the W F, the correction data X with respect to the image data X of the pixel of interest P '''is derived according to the following equation.
X '''= X + Σe P + W E · Σe E + W F Σe F ......... (7)
[0063]
When the correction of the image data of the target pixel P is completed in this way, the error values for the pixels P, E, and F in the error memory are updated according to the following equation.
Σe P → 0
Σe E → (1-W E ) · Σe E ......... (8)
Σe F → (1-W F ) · Σe F
[0064]
The subsequent operation is exactly the same as in the first embodiment.
[0065]
FIG. 9 is a block diagram of an image data binarization device according to a second embodiment of the present invention.
The multi-level image data X for the pixel P is input to the
[0066]
Distribution ratio W E corresponding to the image data X of the pixel of interest P, as a method of determining the W F, for example, receives image data X as an address in the memory, W E as memory contents corresponding to the address, W F Should be stored. Alternatively, the image data X and the error value Sigma] e E of the pixel P, enter the Sigma] e F as an address in the memory, if storing the image data X + W E · Σe E + W F · Σe F as memory contents corresponding to the address , Correction data X ′ ″ is quickly derived.
[0067]
When the correction of the image data of the target pixel P is completed in this way, the error data Σe P , Σe E , and Σe F in the
[0068]
Update data (1-W E) · Σe E, (1-W F) · Σe as a method of deriving F, for example, the image data X and the error value Sigma] e E of the target pixel P or image data of the pixel of interest P, X and error values Sigma] e F input to the memory as an address, as a memory content corresponding to this address, (1-W E) · Σe E or (1-W F) · Σe F
Should be stored.
The subsequent operation is exactly the same as in the first embodiment.
[0069]
FIG. 10 is an explanatory diagram showing a printing result by a printer in a low density area (15% to 0%) according to the present embodiment, and corresponds to FIG. 16 of the conventional example.
[0070]
FIG. 11 is an explanatory diagram showing a printing result of a high density area (85% to 100%) by the printer according to the present embodiment, and corresponds to FIG. 17 of the conventional example.
[0071]
FIG. 12 is an explanatory diagram showing a printing result of a low density area (1% constant) by the printer according to the present embodiment, and corresponds to FIG. 18 of the conventional example.
[0072]
FIG. 13 is an explanatory diagram showing a printing result of a high density area (99% constant) by the printer according to the present embodiment, and corresponds to FIG. 19 of the conventional example.
In each case, the black point and the white point are arranged evenly, and the image quality is remarkably improved.
[0073]
According to the present embodiment, in the low-density region (1%) and the high-density region (99%), a slight improvement compared to the first embodiment is recognized.
[0074]
Further, in the present embodiment, the processing for the portion of the density region of 10% to 90% is not different from the conventional method, so that the time required for the binarization processing is not different from the conventional method.
[Brief description of the drawings]
FIG. 1 is an explanatory diagram for explaining an error diffusion method according to a first embodiment of the present invention.
FIG. 2 is a block diagram showing an image data binarization device as a first embodiment of the present invention.
FIG. 3 is an explanatory diagram showing a recording result in a low density area (0% to 15%) according to the first embodiment of the present invention.
FIG. 4 is an explanatory diagram showing a recording result in a high density area (85% to 100%) according to the first embodiment of the present invention.
FIG. 5 is an explanatory diagram showing a recording result for a density of 1% according to the first embodiment of the present invention.
FIG. 6 is an explanatory diagram showing recording results for a density of 99% according to the first embodiment of the present invention.
FIG. 7 is an explanatory diagram showing a recording result in an intermediate density area (15% to 85%) according to the first embodiment of the present invention.
FIG. 8 is an explanatory diagram for explaining an error diffusion method according to a second embodiment of the present invention.
FIG. 9 is a block diagram illustrating an image data binarization device according to a second embodiment of the present invention.
FIG. 10 is an explanatory diagram showing a recording result in a low density area (0% to 15%) according to the second embodiment of the present invention.
FIG. 11 is an explanatory diagram showing a recording result in a high density area (85% to 100%) according to the second embodiment of the present invention.
FIG. 12 is an explanatory diagram showing a recording result for a density of 1% according to a second embodiment of the present invention.
FIG. 13 is an explanatory diagram showing a recording result for a density of 99% according to the second embodiment of the present invention.
FIG. 14 is a block diagram showing a conventional image data binarization device.
FIG. 15 is an explanatory diagram for explaining the principle of the error diffusion method.
FIG. 16 is an explanatory diagram showing a printing result in a low density area (0% to 15%) by a printer when a conventional image data binarization device is used.
FIG. 17 is an explanatory diagram showing a printing result in a high density area (85% to 100%) by a printer when a conventional image data binarization device is used.
FIG. 18 is an explanatory diagram showing a printing result for a density of 1% by a printer when a conventional image data binarization device is used.
FIG. 19 is an explanatory diagram showing a printing result of a printer at a density of 99% when a conventional image data binarization device is used.
FIG. 20 is an explanatory diagram showing a printing result of an intermediate density area (15% to 85%) by a printer when a conventional image data binarization device is used.
[Explanation of symbols]
Claims (3)
(a)注目画素の画像データに、前記注目画素に蓄積された誤差と、前記注目画素の近傍に位置する2値化前の非注目画素に蓄積された誤差とを加算し、該加算後の画像データを前記注目画素の修正データとする工程と、
(b)前記注目画素および前記非注目画素に蓄積された誤差を更新する工程と、
(c)前記補正データと閾値に基づき前記注目画素の画像データを2値化する工程と、
(d)前記補正データと前記2値化の結果に基づき前記注目画素の2値化に伴って発生する誤差を導き出す工程と、
(e)前記導き出された誤差を2値化前の画素に分配する工程と、
(f)分配された前記誤差を画素毎に誤差データとして蓄積する工程と、
によって成る画像データ2値化方法。An image data binarization method for binarizing multi-valued image data,
(A) adding, to the image data of the target pixel, the error accumulated in the target pixel and the error accumulated in the non-target pixel before binarization located in the vicinity of the target pixel; Making image data correction data of the pixel of interest;
(B) updating the error accumulated in the target pixel and the non-target pixel;
(C) binarizing the image data of the pixel of interest based on the correction data and a threshold;
(D) deriving an error generated with the binarization of the target pixel based on the correction data and the binarization result;
(E) distributing the derived error to pixels before binarization;
(F) accumulating the distributed error as error data for each pixel;
Image data binarization method comprising:
前記工程(a)においては、前記非注目画素に蓄積された誤差に前記係数を乗じた値を、前記非注目画素に蓄積された誤差の代わりに前記注目画素の画像データに加算することを特徴とする請求項1記載の画像データ2値化方法。A step of multiplying the error accumulated in the non-target pixel by a coefficient of 0 to 1 determined according to image data of the target pixel,
In the step (a), a value obtained by multiplying the error accumulated in the non-target pixel by the coefficient is added to the image data of the target pixel instead of the error accumulated in the non-target pixel. 2. The image data binarization method according to claim 1, wherein:
注目画素の画像データに、前記注目画素に蓄積された誤差と、前記注目画素の近傍に位置する2値化前の非注目画素に蓄積された誤差とを加算することにより前記注目画素に対する補正データを導出する補正データ導出手段と、
前記注目画素および前記注非目画素に蓄積された誤差を更新する更新手段と、
前記補正データと閾値に基づき前記注目画素の画像データを2値化する2値化手段と、
前記補正データと前記2値化の結果に基づき前記注目画素の2値化に伴って発生する誤差を2値化前の画素に分配する誤差分配手段と、
前記分配された誤差を画素毎に誤差データとして蓄積する誤差蓄積手段とを備えたことを特徴とする画像データ2値化装置。An image data binarization device for binarizing multi-valued image data,
Correction data for the pixel of interest is obtained by adding the error accumulated in the pixel of interest to the image data of the pixel of interest and the error accumulated in the non-target pixel before binarization located in the vicinity of the pixel of interest. Correction data deriving means for deriving
Updating means for updating the error accumulated in the noted pixel and the noted pixel,
Binarizing means for binarizing the image data of the pixel of interest based on the correction data and a threshold value;
Error distributing means for distributing an error generated due to binarization of the target pixel based on the correction data and the binarization result to pixels before binarization;
An image data binarizing device comprising: an error accumulating means for accumulating the distributed error as error data for each pixel.
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP09623096A JP3574711B2 (en) | 1996-03-25 | 1996-03-25 | Image data binarization method and apparatus |
US08/787,587 US5911009A (en) | 1996-01-25 | 1997-01-22 | Method and apparatus for binary coding of image data including adding error accumulated for a target pixel and a pixel in the vicinity to be later coded |
DE69729965T DE69729965D1 (en) | 1996-01-25 | 1997-01-24 | Method and device for binary coding of image data |
EP97101145A EP0786741B1 (en) | 1996-01-25 | 1997-01-24 | Method and apparatus for binary coding of image data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP09623096A JP3574711B2 (en) | 1996-03-25 | 1996-03-25 | Image data binarization method and apparatus |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH09261465A JPH09261465A (en) | 1997-10-03 |
JP3574711B2 true JP3574711B2 (en) | 2004-10-06 |
Family
ID=14159434
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP09623096A Expired - Fee Related JP3574711B2 (en) | 1996-01-25 | 1996-03-25 | Image data binarization method and apparatus |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3574711B2 (en) |
-
1996
- 1996-03-25 JP JP09623096A patent/JP3574711B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JPH09261465A (en) | 1997-10-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1073258B1 (en) | Error diffusion pattern shifting reduction through programmable threshold perturbation | |
JP3313399B2 (en) | Image processing device | |
JP4121631B2 (en) | Image data processing system and image data processing method | |
JPH10271331A (en) | Image processing method and device therefor | |
US5805738A (en) | Image processing apparatus and method | |
JPH11187264A (en) | Method and device for processing image | |
US5204760A (en) | System and method for converting continuous half tone image into pseudo half tone image | |
JP4068181B2 (en) | Multilevel gradation pixel value level reduction method and system | |
EP0786741B1 (en) | Method and apparatus for binary coding of image data | |
JP3574711B2 (en) | Image data binarization method and apparatus | |
JPH0738767A (en) | Image binarizing processor | |
JP3225099B2 (en) | Image processing device | |
JP3124589B2 (en) | Image processing device | |
JP2000270210A (en) | Picture processor | |
JP3124604B2 (en) | Image processing device | |
JP3157870B2 (en) | Image processing method | |
JP2738865B2 (en) | Image processing device | |
JPH09214760A (en) | Picture data binarizing device | |
JP3549475B2 (en) | Image processing method | |
JP3059457B2 (en) | Image processing device | |
JPH09214759A (en) | Method for binarizing picture data and its device | |
JPH0197066A (en) | Picture processing method | |
JP4185720B2 (en) | Image processing apparatus and image processing method | |
JPH10257303A (en) | Image-processing method and device therefor | |
JPH09205545A (en) | Method and device for binarizing image data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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: 20040622 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20040705 |
|
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: 20080709 Year of fee payment: 4 |
|
LAPS | Cancellation because of no payment of annual fees |