JP2002027249A - Half tone processor, and medium having stored its procedure - Google Patents

Half tone processor, and medium having stored its procedure

Info

Publication number
JP2002027249A
JP2002027249A JP2000212464A JP2000212464A JP2002027249A JP 2002027249 A JP2002027249 A JP 2002027249A JP 2000212464 A JP2000212464 A JP 2000212464A JP 2000212464 A JP2000212464 A JP 2000212464A JP 2002027249 A JP2002027249 A JP 2002027249A
Authority
JP
Japan
Prior art keywords
image
halftone processing
block
halftone
processing device
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.)
Pending
Application number
JP2000212464A
Other languages
Japanese (ja)
Inventor
Yoshito Abe
淑人 阿部
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dai Nippon Printing Co Ltd
Original Assignee
Dai Nippon Printing Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Dai Nippon Printing Co Ltd filed Critical Dai Nippon Printing Co Ltd
Priority to JP2000212464A priority Critical patent/JP2002027249A/en
Publication of JP2002027249A publication Critical patent/JP2002027249A/en
Pending legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/40Picture signal circuits
    • H04N1/405Halftoning, i.e. converting the picture signal of a continuous-tone original into a corresponding signal showing only two levels
    • H04N1/4051Halftoning, 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/4052Halftoning, 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Image Processing (AREA)
  • Facsimile Image Signal Circuits (AREA)

Abstract

PROBLEM TO BE SOLVED: To provide a half tone processor which can get a high-quality half tone image high in reproducibility, and can finish the processing in a short time, and a medium having recorded its procedure. SOLUTION: These are a half tone processor consisting of a block dividing means, which divides a input image into blocks consisting of specified size of picture elements, an adaptation computing means which makes the picture element of each block a vector and compares this with a code word to serve as a candidate of conversion and computes the adaptation, an image generating means which generates an image in the block, based on the code word where the adaptation is maximum, an error diffusing means which diffuses the average stratum error of the block to the neighboring block of a block being processed thereafter, and a storage medium having recorded its procedure.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、コンピュータによ
り画像データを処理する技術分野に属する。特に、濃淡
画像からハーフトーン画像(擬似濃淡画像)を生成する
ハーフトーン処理に関する。
The present invention belongs to the technical field of processing image data by a computer. In particular, the present invention relates to a halftone process for generating a halftone image (pseudo grayscale image) from a grayscale image.

【0002】[0002]

【従来技術】ハーフトーン処理は、原理的に表現可能な
階調数に制限を有する再現システムにおいて、その制限
を越えて本来の濃淡画像(フルトーン画像)により近い
表現ができるようにする処理である。たとえば、オフセ
ット印刷においては、原理的にインキが印刷用紙に着肉
する網点の部分と着肉しない用紙そのものの部分との2
つの部分だけで画像を再現する。カラー印刷の場合に
は、CMYK(cyan,magenta,yellow,black)の4色の
網点と、その網点の刷り重ねの組み合わせにより、16
色だけで画像を再現する。
2. Description of the Related Art Halftone processing is a processing which, in a reproduction system having a limit on the number of gray scales that can be expressed in principle, allows expression closer to the original grayscale image (full-tone image) beyond the limit. . For example, in offset printing, in principle, two parts, ie, a halftone dot portion where the ink is deposited on the printing paper and a portion of the paper itself which is not deposited on the printing paper.
Reproduce the image with only one part. In the case of color printing, a combination of four color dots of CMYK (cyan, magenta, yellow, black) and the overprinting of the halftone dots is used.
Reproduce an image using only colors.

【0003】周知のように、オフセット印刷において
は、このような制限を越えるため網点%を変化させるこ
とが行なわれる。網点%は印刷面において網点が占める
面積の比率である。この網点は微細であるため、その微
細構造が見えない程度の遠目(およそ20cm以上)で
見ると、視覚の特性により中間色として認識され、擬似
的にフルカラーフルトーン画像が再現される。
As is well known, in offset printing, halftone dot% is changed to exceed such a limit. Halftone dot% is the ratio of the area occupied by halftone dots on the printed surface. Since these halftone dots are fine, when viewed from a distance (approximately 20 cm or more) at which the fine structure cannot be seen, the halftone dots are recognized as intermediate colors due to the visual characteristics, and a full-color full-tone image is reproduced in a pseudo manner.

【0004】[0004]

【発明が解決しようとする課題】オフセット印刷におい
てこの網点を生成する処理には、従来は、コンタクト・
スクリーンが用いられる。すなわち、コンタクト・スク
リーンと原画像のフィルムを重ねて感光フィルムに密着
露光し現像することにより網点画像のフィルムを得てい
る。コンタクト・スクリーンは周期的画像であるから、
この網点画像における網点も周期的に配列する(図7
(A)参照)。
Conventionally, the process of generating the halftone dots in offset printing includes a contact / printing process.
A screen is used. That is, a halftone dot film is obtained by superposing the contact screen and the film of the original image on the photosensitive film, and exposing the photosensitive film to exposure. Since the contact screen is a periodic image,
The dots in this dot image are also arranged periodically (FIG. 7).
(A)).

【0005】この周期的な網点をカラー印刷において刷
り重ねるとモアレを発生する原因となる。これを回避す
る必要性から、実際は各色においてスクリーン角度を変
化させている。しかし、モアレの発生を回避できても、
視覚上の障害であり印刷品質を低下させる亀甲模様(ロ
ゼッタ・パターン)等が残存し、その発生を同時に回避
することができない(図6参照)。
[0005] Overprinting these periodic halftone dots in color printing causes moiré. In order to avoid this, the screen angle is actually changed for each color. However, even if moiré can be avoided,
A glazing pattern (rosetta pattern) which is a visual obstacle and lowers the print quality remains and cannot be avoided at the same time (see FIG. 6).

【0006】この亀甲模様の発生は、網点が周期的であ
ることとともに、網点を構成する微細なドットがその周
期的部位に集中していることに原因があると認識されて
いる。そこで、網点を分散型にしたり非周期的にするこ
とが提案されている。これは、近年になってコンピュー
タ処理により任意のドットを生成することが可能となっ
たためである。
[0006] It is recognized that the generation of the turtle pattern is caused by the fact that the halftone dots are periodic and the fine dots constituting the halftone dots are concentrated in the periodic part. Therefore, it has been proposed to make the halftone dots dispersed or aperiodic. This is because arbitrary dots can be generated by computer processing in recent years.

【0007】その処理の基本は、原画像と所定サイズの
閾値行列(ディザマトリックス)との比較演算である。
この方法として、固定マスク法、組織的ディザ法、ラン
ダムディザ法、等が知られているが、全体画像において
は閾値行列が規則的に配列されることから周期性を排除
できないという問題を有する(図7(B)参照)。
The basis of the processing is a comparison operation between an original image and a threshold matrix (dither matrix) of a predetermined size.
As this method, a fixed mask method, a systematic dither method, a random dither method, and the like are known, but there is a problem that the periodicity cannot be excluded because the threshold matrix is regularly arranged in the entire image ( FIG. 7 (B)).

【0008】この周期性の認識を低下させるために、
ブロックサイズの大型化、異なる閾値行列の適用、
原画像に対するディザ値、乱数、ブルーノイズの加算処
理、等の適用が提案されているが、やはり周期性が認識
される(図8(D)参照)。
In order to reduce the recognition of the periodicity,
Larger block size, application of different threshold matrices,
Although application of a dither value, a random number, addition processing of blue noise, and the like to the original image has been proposed, periodicity is also recognized (see FIG. 8D).

【0009】閾値行列との比較演算に基づかない方法と
しては、平均誤差最小法、誤差拡散(ED)法、等が知
られている。周期性が無いという利点とともに、量子化
雑音の低減を主眼とするため平均階調の保存性や細部の
再現性に優れている。しかし、粒状感の残存、階調特性
の変化、独特のワーム状のテクスチャを発生するという
問題がある(図8(C)参照)。また、演算量の増加に
よりデータ処理の負荷が大きく相応の処理時間を必要と
するという問題がある。
As a method not based on a comparison operation with a threshold matrix, there are known an average error minimum method, an error diffusion (ED) method, and the like. In addition to the advantage that there is no periodicity, the preservation of the average gradation and the reproducibility of details are excellent because the main purpose is to reduce the quantization noise. However, there is a problem that the graininess remains, the gradation characteristic changes, and a unique worm-like texture is generated (see FIG. 8C). In addition, there is a problem that the load of data processing is large due to an increase in the amount of calculation and a corresponding processing time is required.

【0010】本発明は上記の問題を解決するためになさ
れたものである。その目的は、再現性の高い、高品質の
ハーフトーン画像を得ることができ、しかもその処理を
短時間に済ますことができるハーフトーン処理装置およ
びその処理手順を記録した媒体を提供することにある。
The present invention has been made to solve the above problems. An object of the present invention is to provide a halftone processing apparatus which can obtain a high-quality halftone image with high reproducibility and can perform the processing in a short time, and a medium in which the processing procedure is recorded. .

【0011】[0011]

【課題を解決するための手段】上記課題は下記の本発明
によって解決される。すなわち、本発明の請求項1に係
るハーフトーン処理装置は、L,M,Nを正の整数とし
て、Mビット表現の画素値を有する画像から、M>Nで
あるNビット表現の画素値を有する画像を導出するハー
フトーン処理装置であって、ブロック分割手段と、適合
度算出手段と、画像生成手段と、誤差拡散手段とを有
し、前記ブロック分割手段は、入力画像を所定サイズL
の画素からなるブロックに分割し、前記適合度算出手段
は、前記分割されたブロックの画素をL次元Mビットの
ベクトルとし、これを変換候補となるL次元Nビットの
コードワードと比較し適合度を算出し、前記画像生成手
段は適合度が最大のコードワードに基づいて前記ブロッ
クにL次元Nビットの画像を生成し、前記誤差拡散手段
は前記ブロックの平均階調誤差を以降に処理する前記ブ
ロックの近傍ブロックへ誤差拡散するようにしたもので
ある。
The above object is achieved by the present invention described below. That is, the halftone processing apparatus according to claim 1 of the present invention converts an N-bit pixel value with M> N from an image having an M-bit pixel value, where L, M, and N are positive integers. A halftone processing device for deriving an image having the image, comprising: a block dividing unit; a fitness calculating unit; an image generating unit; and an error diffusion unit.
, And the degree-of-fit calculation means converts the pixels of the divided block into an L-dimensional M-bit vector, compares the vector with an L-dimensional N-bit code word that is a conversion candidate, and The image generating means generates an L-dimensional N-bit image in the block based on the codeword having the highest fitness, and the error diffusion means processes the average gradation error of the block thereafter. The error diffusion is performed to a block adjacent to the block.

【0012】本発明によれば、ブロック分割手段により
入力画像が所定サイズLの画素からなるブロックに分割
され、適合度算出手段によりそのブロックの画素がL次
元Mビットのベクトルとされ、これを変換候補となるL
次元Nビットのコードワードと比較し適合度が算出さ
れ、画像生成手段により適合度が最大のコードワードに
基づいてそのブロックにL次元Nビットの画像が生成さ
れ、誤差拡散手段によりそのブロックの平均階調誤差を
以降に処理するそのブロックの近傍ブロックへ誤差拡散
する。すなわち、ブロックに生成される画像は適合度が
最大のコードワードに基づくものであるからハーフトー
ン画像の全体において最適化された画像を得ることがで
きる。また、誤差拡散が行なわれるからその画像は階調
再現性が優れている。したがって、再現性の高い、高品
質のハーフトーン画像が得られるハーフトーン処理装置
が提供される。
According to the present invention, the input image is divided into blocks of pixels of a predetermined size L by the block dividing means, and the pixels of the blocks are converted into L-dimensional M-bit vectors by the fitness calculating means. Candidate L
The degree of conformity is calculated by comparing with the code word of dimension N bits, an image of L dimension N bits is generated in the block based on the code word having the maximum degree of conformity by the image generation means, and the average of the block is averaged by the error diffusion means. The tone error is diffused to a block adjacent to the block to be processed thereafter. That is, since the image generated in the block is based on the codeword having the highest degree of matching, an optimized image can be obtained from the entire halftone image. Further, since the error diffusion is performed, the image has excellent gradation reproducibility. Therefore, there is provided a halftone processing device capable of obtaining a high-quality halftone image with high reproducibility.

【0013】また本発明の請求項2に係るハーフトーン
処理装置は、請求項1に係るハーフトーン処理装置にお
いて、前記誤差拡散手段が行なう前記誤差拡散は、前記
ブロックの平均階調誤差を算出する過程でその誤差を相
殺することで実行するようにしたものである。本発明に
よれば、平均階調誤差を算出する過程でその誤差が相殺
されるから、誤差拡散における処理の負荷が極めて小さ
い。
According to a second aspect of the present invention, in the halftone processing apparatus according to the first aspect, the error diffusion performed by the error diffusion means calculates an average gradation error of the block. The process is executed by canceling the error in the process. According to the present invention, since the error is canceled in the process of calculating the average gradation error, the processing load in error diffusion is extremely small.

【0014】また本発明の請求項3に係るハーフトーン
処理装置は、請求項1または2に係るハーフトーン処理
装置において、J,Kを正の整数として、前記ブロック
分割手段は、入力画像を所定サイズであるJ×K=Lの
画素からなる矩形のブロックに分割するようにしたもの
である。本発明によれば、矩形のブロックにより処理が
行なわれ基本的で最も簡明な処理あるとともに、水平
線、垂直線の再現性が優れる。
The halftone processing apparatus according to a third aspect of the present invention is the halftone processing apparatus according to the first or second aspect, wherein J and K are positive integers, and the block dividing means converts the input image into a predetermined image. It is divided into rectangular blocks composed of pixels of size J × K = L. According to the present invention, processing is performed by rectangular blocks, and the basic and simplest processing is performed, and the reproducibility of horizontal lines and vertical lines is excellent.

【0015】また本発明の請求項4に係るハーフトーン
処理装置は、請求項1〜3のいずれかに係るハーフトー
ン処理装置において、前記ブロック分割手段は、入力画
像を所定サイズLの画素からなる、六角形または三角
形、あるいはそれらを組み合わせたタイリング(tilin
g)形のブロックに分割するようにしたものである。本
発明によれば、周期性の認識が緩和され方向による再現
性が均等化される。
According to a fourth aspect of the present invention, in the halftone processing apparatus according to any one of the first to third aspects, the block dividing means converts the input image into pixels of a predetermined size L. , Hexagons or triangles, or tilings that combine them (tilin
g) is divided into blocks. According to the present invention, the recognition of the periodicity is eased, and the reproducibility in the directions is equalized.

【0016】また本発明の請求項5に係るハーフトーン
処理装置は、請求項1〜4のいずれかに係るハーフトー
ン処理装置において、適合度を算出するための最適化パ
ラメータを設定する最適化パラメータ設定手段を有する
ようにしたものである。本発明によれば、最適化パラメ
ータを設定することにより適合度の算出が適正化される
から、より再現性の高い、より高品質のハーフトーン画
像が得られる。
According to a fifth aspect of the present invention, there is provided a halftone processing apparatus according to any one of the first to fourth aspects, wherein an optimization parameter for setting an optimization parameter for calculating a degree of conformity is set. This has a setting means. According to the present invention, since the calculation of the degree of conformity is optimized by setting the optimization parameter, a higher-quality halftone image with higher reproducibility can be obtained.

【0017】また本発明の請求項6に係るハーフトーン
処理装置は、請求項1〜5のいずれかに係るハーフトー
ン処理装置において、あらかじめ生成したL次元Nビッ
トのコードワードが記載されたコードブックを記憶する
コードブック記憶手段を有するようにしたものである。
本発明によれば、コードワードはあらかじめ生成された
コードブックのものが適用される。したがって、ハーフ
トーン処理においてはコードワード演算の時間だけ処理
時間が短縮される。また、コードワードの数についてあ
らかじめ適正な減縮を行なうことができるから、さらに
処理時間が短縮される。
According to a sixth aspect of the present invention, there is provided the halftone processing apparatus according to any one of the first to fifth aspects, wherein the codebook in which an L-dimensional N-bit code word generated in advance is described. Is stored in the storage unit.
According to the present invention, a codeword from a previously generated codebook is applied. Therefore, in the halftone processing, the processing time is reduced by the time of the codeword operation. Also, since the number of codewords can be appropriately reduced in advance, the processing time is further reduced.

【0018】また本発明の請求項7に係るハーフトーン
処理装置は、請求項1〜5のいずれかに係るハーフトー
ン処理装置において、L次元Nビットのコードワードを
逐次生成するコードワード逐次生成手段を有するように
したものである。本発明によれば、ハーフトーン処理の
状態に適合したコードワードの生成を行なうことがで
き、不必要なコードワードに係わる処理が省略される。
A halftone processing apparatus according to a seventh aspect of the present invention is the halftone processing apparatus according to any one of the first to fifth aspects, wherein the codeword sequential generation means for sequentially generating an L-dimensional N-bit codeword. Is provided. According to the present invention, a codeword suitable for the state of the halftone processing can be generated, and processing relating to unnecessary codewords is omitted.

【0019】また本発明の請求項8に係るハーフトーン
処理装置は、請求項1〜7のいずれかに係るハーフトー
ン処理装置において、前記L次元Nビットのコードワー
ドは、その総数を2^L(2のL乗)未満とするように
したものである。本発明によれば、ハーフトーン処理が
高速化される。
According to an eighth aspect of the present invention, in the halftone processing apparatus according to any one of the first to seventh aspects, the total number of the L-dimensional N-bit codewords is 2 ^ L. (2 to the power of L). According to the present invention, the speed of the halftone process is increased.

【0020】また本発明の請求項9に係るハーフトーン
処理装置は、請求項1〜8のいずれかに係るハーフトー
ン処理装置において、そのハーフトーン処理装置は複数
の処理部を具備し、前記ブロック分割手段が分割したブ
ロックを前記複数の処理部の各々に割り当て、前記複数
の処理部によって並列分散処理するようにしたものであ
る。本発明によれば、ハーフトーン処理が高速化され
る。
According to a ninth aspect of the present invention, there is provided the halftone processing apparatus according to any one of the first to eighth aspects, wherein the halftone processing apparatus includes a plurality of processing units, The block divided by the dividing means is assigned to each of the plurality of processing units, and the plurality of processing units perform parallel distributed processing. According to the present invention, the speed of the halftone process is increased.

【0021】また本発明の請求項10に係るハーフトー
ン処理装置は、請求項9に係るハーフトーン処理装置に
おいて、前記適合度算出手段が行なう前記適合度の算出
を前記複数の処理部によって並列分散処理するようにし
たものである。本発明によれば、ハーフトーン処理が高
速化される。
According to a tenth aspect of the present invention, in the halftone processing apparatus according to the ninth aspect, the calculation of the fitness performed by the fitness calculating means is performed in parallel by the plurality of processing units. This is to be processed. According to the present invention, the speed of the halftone process is increased.

【0022】また本発明の請求項11に係るハーフトー
ン処理装置は、請求項1〜10のいずれかに係るハーフ
トーン処理装置において、前記適合度算出手段による適
合度の算出は、前記コードワードを全探索するのではな
く、適合度が最大のコードワードを最適化手法により探
索するようにしたものである。本発明によれば、ハーフ
トーン処理が著しく高速化される。
In the halftone processing apparatus according to claim 11 of the present invention, in the halftone processing apparatus according to any one of claims 1 to 10, the calculation of the fitness by the fitness calculating means includes the step of: Instead of performing a full search, a codeword with the highest fitness is searched by an optimization technique. According to the present invention, the halftone processing is remarkably speeded up.

【0023】また本発明の請求項12に係るハーフトー
ン処理装置は、請求項1〜11のいずれかに係るハーフ
トーン処理装置において、Pを2以上の整数として前記
画像はP色画像であって、前記ブロックの画素を(L×
P)次元Mビットのベクトルとして処理するようにした
ものである。本発明によれば、多色画像を対象とするハ
ーフトーン処理が行なわれる。
The halftone processing apparatus according to claim 12 of the present invention is the halftone processing apparatus according to any one of claims 1 to 11, wherein P is an integer of 2 or more, and the image is a P color image. , The pixels of the block are (L ×
P) It is processed as a vector of dimension M bits. According to the present invention, halftone processing is performed on a multicolor image.

【0024】また本発明の請求項13に係るハーフトー
ン処理装置は、請求項1〜11のいずれかに係るハーフ
トーン処理装置において、前記画像は多色画像であっ
て、その多色画像の各色画像ごとに前記処理を行なうよ
うにしたものである。本発明によれば、多色画像を対象
とするハーフトーン処理が高速化される。
According to a thirteenth aspect of the present invention, in the halftone processing apparatus according to any one of the first to eleventh aspects, the image is a multicolor image, and each color of the multicolor image is The above processing is performed for each image. According to the present invention, the speed of halftone processing for a multicolor image is increased.

【0025】また本発明の請求項14に係るハーフトー
ン処理装置は、請求項13に係るハーフトーン処理装置
において、前記各色画像ごとの処理を前記複数の処理部
によって前記各色画像ごとに並列分散処理するようにし
たものである。本発明によれば、多色画像を対象とする
ハーフトーン処理が著しく高速化される。
According to a fourteenth aspect of the present invention, in the halftone processing apparatus according to the thirteenth aspect, the processing for each color image is performed by the plurality of processing units in parallel distributed processing for each color image. It is something to do. According to the present invention, halftone processing for a multicolor image is remarkably speeded up.

【0026】また本発明の請求項15に係るハーフトー
ン処理手順が記録された媒体は、L,M,Nを正の整数
として、Mビット表現の画素値を有する画像から、M>
NであるNビット表現の画素値を有する画像を導出する
ハーフトーン処理手順が記録された媒体であって、前記
ハーフトーン処理手順は、ブロック分割過程と、適合度
算出過程と、画像生成過程とを有し、前記ブロック分割
過程において、入力画像を所定サイズLの画素からなる
ブロックに分割し、前記適合度算出過程において、前記
分割されたブロックの画素をL次元Mビットのベクトル
とし、これを変換候補となるL次元Nビットのコードワ
ードと比較し適合度を算出し、前記画像再生過程におい
て、適合度が最大のコードワードに基づいて前記ブロッ
クにL次元Nビットの画像を生成され、前記誤差拡散過
程において、前記ブロックの平均階調誤差を以降に処理
する前記ブロックの近傍ブロックへ誤差拡散するように
したものである。
The medium on which the halftone processing procedure according to claim 15 of the present invention has been recorded can be obtained from an image having pixel values of M bits where L, M, and N are positive integers.
A medium in which a halftone processing procedure for deriving an image having a pixel value of N bits represented by N is recorded, wherein the halftone processing procedure includes a block division step, a fitness calculation step, an image generation step, In the block dividing step, the input image is divided into blocks each having a pixel of a predetermined size L, and in the fitness calculating step, the pixels of the divided block are defined as L-dimensional M-bit vectors, An L-dimensional N-bit image is generated by comparing with an L-dimensional N-bit code word that is a conversion candidate, and in the image reproducing process, an L-dimensional N-bit image is generated in the block based on the code word having the highest degree of conformity. In the error diffusion process, the average gradation error of the block is error-diffused to a block adjacent to the block to be processed thereafter.

【0027】本発明によれば、ブロック分割過程におい
て入力画像が所定サイズLの画素からなるブロックに分
割され、適合度算出過程においてそのブロックの画素が
L次元Mビットのベクトルとされ、これを変換候補とな
るL次元Nビットのコードワードと比較し適合度が算出
され、画像生成手段において適合度が最大のコードワー
ドに基づいてそのブロックにL次元Nビットの画像が生
成される。すなわち、ブロックに生成される画像は適合
度が最大のコードワードに基づくものであるからハーフ
トーン画像の全体において最適化された画像を得ること
ができる。また、誤差拡散が行なわれるからその画像は
階調再現性が優れている。したがって、再現性の高い、
高品質のハーフトーン画像が得られるハーフトーン処理
の処理手順を記録した媒体が提供される。
According to the present invention, the input image is divided into blocks each having a predetermined size L in the block division process, and the pixels of the block are converted into L-dimensional M-bit vectors in the fitness calculation step. Comparability is calculated by comparing with a candidate L-dimensional N-bit codeword, and an L-dimensional N-bit image is generated in the block based on the codeword having the highest conformity in the image generation means. That is, since the image generated in the block is based on the codeword having the highest degree of matching, an optimized image can be obtained from the entire halftone image. Further, since the error diffusion is performed, the image has excellent gradation reproducibility. Therefore, high reproducibility,
A medium in which a processing procedure of a halftone process for obtaining a high-quality halftone image is provided.

【0028】[0028]

【発明の実施の形態】次に、本発明について実施の形態
を説明する。本発明は、単純な誤差拡散を基本技術とす
るものである。さらに、本発明は、ベクトル量子化、最
適化、等の基本技術に対して誤差拡散を加え拡張したも
のである。したがって、閾値処理、等を基本技術とする
従来の方法とは本質的に異なっている。
Next, an embodiment of the present invention will be described. The present invention is based on simple error diffusion. Further, the present invention extends the basic techniques such as vector quantization and optimization by adding error diffusion. Therefore, it is fundamentally different from the conventional method based on threshold processing or the like.

【0029】まず、誤差拡散技術について簡単に説明し
ておく。たとえば、ブロックサイズが4画素×4画素で
2値化によるハーフトーン画像を得るとする。この場合
は、単純には、平均階調は1+16=17段階で表現す
ることになる。標準的な8〜12ビット階調のフルトー
ンの原画像をハーフトーン画像化する方法としては、階
調表現能力が不足する。誤差拡散技術により、この階調
表現能力が不足を補足することができる。
First, the error diffusion technique will be briefly described. For example, assume that a halftone image is obtained by binarization with a block size of 4 pixels × 4 pixels. In this case, simply, the average gradation is expressed in 1 + 16 = 17 steps. As a method of converting a full-tone original image of a standard 8 to 12-bit gradation into a halftone image, the gradation expression ability is insufficient. With the error diffusion technique, this lack of gradation expression ability can be complemented.

【0030】誤差拡散処理について説明図を図3に示
す。原画像を8ビット階調とすると、その画素値は0〜
255の値を有する。ハーフトーン画像を1ビット階調
とすると、その画素値は0か255(0か1としてもよ
い)の値を有する。図3において、量子化器の入口に入
力した画素値は量子化器において0か255に変換され
て量子化器の出口に出力される。たとえば、127以下
を0とし128以上を255とする変換が行なわれる。
変換された画素値は量子化器の出口からさらに出力側に
出力される。この変換された画素値がハーフトーン画像
の画素値である。
FIG. 3 is a diagram for explaining the error diffusion process. If the original image has 8-bit gradation, the pixel value is 0 to
It has a value of 255. If the halftone image has a 1-bit gradation, the pixel value has a value of 0 or 255 (or 0 or 1). In FIG. 3, the pixel value input to the entrance of the quantizer is converted to 0 or 255 by the quantizer and output to the exit of the quantizer. For example, conversion is performed in which 127 or less is set to 0 and 128 or more is set to 255.
The converted pixel value is output from the outlet of the quantizer to the output side. This converted pixel value is the pixel value of the halftone image.

【0031】この量子化器における量子化により、0〜
255の値が0か255の値に丸められるのであるか
ら、もとの値と量子化後の値との間には誤差(量子化誤
差)を生じることになる。図3に示すように、量子化器
の入口に入力した画素値は、量子化器の出口に出力され
た画素値(変換された画素値)によって減算される。す
なわち、量子化誤差が演算される。この量子化誤差はフ
ィルターの入口に入力される。
By the quantization in this quantizer, 0 to 0
Since the value of 255 is rounded to 0 or 255, an error (quantization error) occurs between the original value and the value after quantization. As shown in FIG. 3, the pixel value input to the entrance of the quantizer is subtracted by the pixel value (converted pixel value) output to the exit of the quantizer. That is, a quantization error is calculated. This quantization error is input to the entrance of the filter.

【0032】フィルターは誤差をどのように拡散するか
を決定する役割を果たす。最も単純な拡散方法は、単純
遅延を行なって、次の画素に加算する方法である。その
他としては、原画像の走査は、多くは左上の画素から右
下の画素へと行なわれるラスター走査であるから、その
場合に量子化誤差を演算した画素における右方向および
下方向の複数の画素に配分する方法もある。
The filter plays a role in determining how to spread the error. The simplest spreading method is a method of performing a simple delay and adding the result to the next pixel. Otherwise, since scanning of the original image is mostly raster scanning performed from the upper left pixel to the lower right pixel, in this case, a plurality of pixels in the right direction and the lower direction in the pixel on which the quantization error has been calculated are used. There is also a way to distribute.

【0033】図3において、入力側から原画像を走査し
て画素値を順次入力する。その画素値に量子化誤差が加
算される。その量子化誤差が加算された画素値を量子化
器が入力し量子化する。その量子化された画素値は量子
化器の出口に出力される。この変換された画素値がハー
フトーン画像の画素値である。この一連の過程により、
誤差拡散が行なわれる。この結果、比較的細部の再現性
と、平均階調の再現性に優れるハーフトーン画像を得る
ことができる。
In FIG. 3, an original image is scanned from the input side and pixel values are sequentially input. A quantization error is added to the pixel value. The quantizer inputs the pixel value to which the quantization error has been added and quantizes it. The quantized pixel value is output to the exit of the quantizer. This converted pixel value is the pixel value of the halftone image. Through this series of processes,
Error diffusion is performed. As a result, it is possible to obtain a halftone image that is excellent in reproducibility of relatively fine details and reproducibility of average gradation.

【0034】上述においては、1つの数値で表現される
1つの画素ごとに量子化するスカラー量子化における誤
差拡散を説明した。ベクトル量子化における誤差拡散も
基本的に同様である。画像のベクトル量子化の過程にお
いて誤差拡散を適用することについての説明は後述する
ものとして、ここで、ベクトル量子化について簡単に説
明しておく。一般の量子化、すなわちスカラー量子化は
連続あるいは不連続なスカラー値(一次元の値)を、離
散的なスカラー値(量子化代表値)に丸めるものであ
る。たとえば、小数点以下を四捨五入して整数値を得る
操作はスカラー量子化である。この離散的なスカラー値
は等間隔である必要性はなく、整数値である必要性もな
い。
In the above, error diffusion in scalar quantization in which quantization is performed for each pixel represented by one numerical value has been described. The error diffusion in the vector quantization is basically the same. The description of applying error diffusion in the process of image vector quantization will be described later, and here, the vector quantization will be briefly described. In general quantization, that is, scalar quantization, a continuous or discontinuous scalar value (one-dimensional value) is rounded to a discrete scalar value (quantized representative value). For example, the operation of rounding off the decimal part to obtain an integer value is scalar quantization. The discrete scalar values need not be equally spaced and need not be integer values.

【0035】一方、ベクトル量子化はベクトル値(多次
元の値)に発展させたものである。ベクトル量子化は、
ベクトルによって張られる空間を、量子化代表値を母点
とするボロノイ多角形(多次元の多角形)に分割し、そ
の多角形に含まれるベクトル値はすべてその多角形の母
点である量子化代表値に丸めるものである。
On the other hand, vector quantization is developed into vector values (multidimensional values). Vector quantization is
The space spanned by the vectors is divided into Voronoi polygons (multi-dimensional polygons) with the quantized representative value as the generating point, and all vector values included in the polygon are quantized as the generating points of the polygon. Rounded to a representative value.

【0036】スカラー値における四捨五入の例と、最も
簡単な2次元のベクトル値における量子化の例を図4に
示す。図4(A)はスカラー値の場合を示し、図4
(B)は、ベクトル値の場合を示す。図4において、黒
点(黒丸)はボロノイ多角形の母点、すなわち量子化代
表値を表している。量子化は、各境界線で囲まれた領域
内のスカラー値またはベクトル値をその黒点で表される
量子化代表値に丸める処理である。
FIG. 4 shows an example of rounding in a scalar value and an example of quantization in a simplest two-dimensional vector value. FIG. 4A shows a case of a scalar value.
(B) shows the case of a vector value. In FIG. 4, black points (black circles) represent generating points of the Voronoi polygon, that is, quantization representative values. The quantization is a process of rounding a scalar value or a vector value in an area surrounded by each boundary line to a quantization representative value represented by a black point.

【0037】図4において、各境界線は隣接する黒点を
結ぶ線分の垂直2等分線であり等距離点の集合である。
図4と同等の境界線は、如何なる次元のユークリッド空
間であっても存在する。量子化代表値がベクトルのとき
には、量子化代表値を量子化代表値ベクトルとも呼ぶ。
多数ある量子化代表値ベクトルの内で、特定のベクトル
値が丸められる量子化代表値ベクトルは、図4からも明
らかなように、そのベクトル値の最近傍の量子化代表値
ベクトルである。
In FIG. 4, each boundary line is a vertical bisector of a line connecting adjacent black points and is a set of equidistant points.
A boundary line equivalent to that in FIG. 4 exists in any dimension of the Euclidean space. When the quantized representative value is a vector, the quantized representative value is also referred to as a quantized representative value vector.
Among the many quantized representative value vectors, the quantized representative value vector in which a specific vector value is rounded is the nearest quantized representative value vector of the vector value, as is clear from FIG.

【0038】すなわち、特定のベクトル値は、ユークリ
ッド距離が最小となる量子化代表値ベクトルに丸められ
る。数値的には、2つのベクトルにおける対応する要素
の差の2乗和が最小となるように量子化代表値ベクトル
が選択される。ここでは、量子化代表値ベクトルをコー
ドワード、その集合をコードブックと呼ぶ。
That is, a specific vector value is rounded to a quantized representative value vector that minimizes the Euclidean distance. Numerically, the quantized representative value vector is selected such that the sum of squares of the difference between the corresponding elements in the two vectors is minimized. Here, the quantized representative value vector is referred to as a codeword, and a set thereof is referred to as a codebook.

【0039】なお、本発明において、原画像が256値
(8ビット)の階調を有する画像で、得ようとするハー
フトーン画像が2値(1ビット)の画像であるとする。
このとき、ベクトルが2次元(2画素単位)であれば、
量子化代表値ベクトルは(0,0)、(0,255)、
(255,0)、(255,255)の4とおりとな
る。ベクトルが4次元(4画素単位)であれば、16と
おりとなる。ベクトルが4次元の場合の量子化代表値ベ
クトルを図5に示す。図5において、白矩形を“0”、
黒矩形を“255”とすると、左上より順番に平均階調
は、0,64,64,128,64,128,128,
191,64,128,128,191,128,19
1,191,255(四捨五入)であり、それぞれの階
調は、階調の順番に1、4、6、4、1とおりづつであ
る(2項係数になる)。
In the present invention, it is assumed that the original image is an image having 256-level (8-bit) gradation and the halftone image to be obtained is a binary (1-bit) image.
At this time, if the vector is two-dimensional (in units of two pixels),
Quantized representative value vectors are (0,0), (0,255),
(255, 0) and (255, 255). If the vector is four-dimensional (in units of four pixels), there are 16 types. FIG. 5 shows a quantized representative value vector when the vector is four-dimensional. In FIG. 5, the white rectangle is “0”,
Assuming that the black rectangle is “255”, the average gradation is 0, 64, 64, 128, 64, 128, 128,
191,64,128,128,191,128,19
1,191,255 (rounded off), and each gradation is 1, 4, 6, 4, 1 in the order of gradation (to be a binomial coefficient).

【0040】次に、本発明におけるハーフトーン処理に
ついて、一例を挙げて処理過程を説明する。M=8ビッ
ト(256)階調表現の画素値を有する画像(原画像)
から、N=1ビット(2)階調表現の画素値を有する画
像(ハーフトーン画像)を導出するハーフトーン処理に
ついて説明する。当然のことであるが、これから説明す
る処理における具体的な数値は、一例を示すものであ
り、本発明の具体的な実施においては、個々の場合に適
合するように説明とは異なる数値を選択することができ
る。画像は多色画像であってもよいし、画像のサイズや
ブロックのサイズ、原画像の階調数、ハーフトーン画像
の階調数が異なっていてもよい。
Next, the halftone process in the present invention will be described with reference to an example. M = 8 bits (256) Image having pixel values of gradation expression (original image)
, A halftone process for deriving an image (halftone image) having pixel values of N = 1 bit (2) gradation expression will be described. As a matter of course, specific numerical values in the processing to be described below are merely examples, and in specific embodiments of the present invention, numerical values different from the description are selected to suit individual cases. can do. The image may be a multi-color image, or the image size and block size, the number of tones of the original image, and the number of tones of the halftone image may be different.

【0041】本発明のハーフトーン処理装置における処
理過程の一例を図1に示す。まず、図1のステップS1
において、原画像をハーフトーン処理装置の記憶部に入
力する。原画像は1024×1024=1048576
画素の濃淡画像であるものとする。次に、ステップS2
において、所定のブロックのサイズを指定する。一例と
して4画素×4画素のブロックのサイズを指定する。
FIG. 1 shows an example of the processing steps in the halftone processing apparatus of the present invention. First, step S1 in FIG.
In, the original image is input to the storage unit of the halftone processing device. The original image is 1024 × 1024 = 1048576
It is assumed that the image is a grayscale image of pixels. Next, step S2
In, the size of a predetermined block is designated. As an example, a block size of 4 pixels × 4 pixels is specified.

【0042】次に、ステップS3において、原画像を4
画素×4画素のブロックに分割する。すなわちJ=K=
4、L=4×4=16である。原画像が1024×10
24画素の濃淡画像であるから、256×256=65
536のブロックに分割される。このブロックは16次
元ベクトルとして表現することができる。その6553
6個の16次元ベクトルに対してこれから処理を行な
う。原画像からブロックへの分割、および16次元ベク
トルの関係の説明図を図2に示す。図2(A)は原画像
であり、図2(B)はブロックである。
Next, in step S3, the original image is
It is divided into blocks of pixels × 4 pixels. That is, J = K =
4, L = 4 × 4 = 16. Original image is 1024 × 10
Since it is a grayscale image of 24 pixels, 256 × 256 = 65
It is divided into 536 blocks. This block can be represented as a 16-dimensional vector. Its 6553
Processing will be performed on six 16-dimensional vectors. FIG. 2 is an explanatory diagram of the division of the original image into blocks and the relationship between 16-dimensional vectors. FIG. 2A is an original image, and FIG. 2B is a block.

【0043】なお、上述における16次元ベクトルの要
素(画素値)の並び順は簡単のため平方ブロックの左上
から右下にラスタースキャンすることとし、ベクトル
(ブロック)の処理順も画像の左上から右下にラスター
スキャンすることとする。
Note that, for simplicity, the arrangement order of the elements (pixel values) of the 16-dimensional vector is raster-scanned from the upper left to the lower right of the square block, and the processing order of the vector (block) is also from the upper left to the right of the image. Raster scan below.

【0044】次に、ステップS4において、各ベクトル
ごとにその要素の算術平均を求める。すなわち、各ベク
トルの平均階調を求める。算術平均は、一般的には、非
整数となるが、簡単のため整数値に丸めることができ
る。その場合は、各ベクトルの平均階調における階調数
は、原画像のM=8ビット(256)階調表現と一致す
る。勿論、処理時間は長くなるが丸めずに処理すること
もできる。
Next, in step S4, the arithmetic average of the element is determined for each vector. That is, the average gradation of each vector is obtained. The arithmetic average is generally a non-integer, but can be rounded to an integer value for simplicity. In that case, the number of gradations in the average gradation of each vector matches the expression of M = 8 bits (256) gradation of the original image. Of course, the processing time is long but the processing can be performed without rounding.

【0045】次に、ステップS5において、適合度を判
定するための評価関数を設定する。前述したベクトル量
子化では、ユークリッド距離を最小とするコードワード
(量子化代表値ベクトル)が最も適合していると判定さ
れ、そのコードワードを選択することになる。ここで
は、ハーフトーン画像を導出することを目的としている
ため、ステップS4において得た平均階調という因子を
考慮に入れる。
Next, in step S5, an evaluation function for determining the degree of conformity is set. In the above-described vector quantization, a codeword (quantized representative value vector) that minimizes the Euclidean distance is determined to be most suitable, and that codeword is selected. Here, since the purpose is to derive a halftone image, the factor of the average gradation obtained in step S4 is taken into consideration.

【0046】たとえば、評価関数を下記の数1のように
設定する。
For example, an evaluation function is set as in the following equation 1.

【数1】 F(Gi,Ho)=α×|ΔD|+β×ΔQ/L ただし、 F(Gi,Ho) ; 評価関数 Gi ; 原画像のベクトル(原画像のブロック) Ho ; コードワード α、β ; 係数(α+β=1、0≦α、β≦1) |ΔD|; 平均階調の誤差(|x|はxの絶対値を
表す) ΔQ ; ベクトル量子化の誤差 L ; 次元数(16)
F (Gi, Ho) = α × | ΔD | + β × ΔQ / L where F (Gi, Ho); evaluation function Gi; vector of original image (block of original image) Ho; codeword α, β; coefficient (α + β = 1, 0 ≦ α, β ≦ 1) | ΔD |; average gradation error (| x | represents an absolute value of x) ΔQ; vector quantization error L; number of dimensions (16) )

【0047】次に、ステップS6において、適合度を算
出する。適合度は数1に示す評価関数によって与えられ
る。数1に示す評価関数においては、評価関数の値が小
さいほど適合度が高いことになる。したがって、適合度
を算出するということは、評価関数が最小とするような
HjをGiごとに探索するということである。順次、コ
ードワードを発生させ(2の16乗とおり)、数1に示
す評価関数Fを最小化するものを全探索する。
Next, in step S6, the degree of conformity is calculated. The fitness is given by the evaluation function shown in Equation 1. In the evaluation function shown in Expression 1, the smaller the value of the evaluation function is, the higher the fitness is. Therefore, calculating the degree of matching means searching for Hj that minimizes the evaluation function for each Gi. Codewords are sequentially generated (2 16 powers), and a search for minimizing the evaluation function F shown in Expression 1 is fully performed.

【0048】ここで、添え字oは0≦o<2^16(2
の16乗)なる数である。すなわち、コードワードの総
数である。また、添え字iは0≦i<65536なる数
である。すなわち、原画像における16次元のベクトル
(ブロック)の順序数である。また、添え字jは0≦j
<16=Lなる数である。すなわち、ベクトルの各要素
の順序を表す数である。
Here, the subscript o is 0 ≦ o <2 ^ 16 (2
16). That is, the total number of codewords. The subscript i is a number satisfying 0 ≦ i <65536. That is, it is the ordinal number of a 16-dimensional vector (block) in the original image. The subscript j is 0 ≦ j
<16 = L. That is, it is a number representing the order of each element of the vector.

【0049】次に、ステップS7において、適合度が最
大のコードワードに基づいて対応するブロックにL次元
Nビットの画像を生成する。このようにして原画像にお
ける特定の16次元のベクトル(ブロック)を処理する
ことによって、そのブロックにL次元Nビットの画像、
すなわちハーフトーン画像のそのブロックの部分が生成
される。
Next, in step S7, an L-dimensional N-bit image is generated in the corresponding block based on the codeword having the highest degree of matching. By processing a specific 16-dimensional vector (block) in the original image in this manner, an L-dimensional N-bit image,
That is, that block portion of the halftone image is generated.

【0050】次に、ステップS8において、適合度が最
大のコードワードの平均階調と、原画像における16次
元のベクトル(ブロック)の平均階調との差異、すなわ
ち平均階調誤差を演算し、以降に処理するブロックに対
して誤差拡散する。ここにおける、誤差拡散はベクトル
量子化における誤差拡散である。
Next, in step S8, the difference between the average gray level of the code word having the highest fitness and the average gray level of the 16-dimensional vector (block) in the original image, that is, the average gray level error, is calculated. Error diffusion is performed on blocks to be processed thereafter. Here, error diffusion is error diffusion in vector quantization.

【0051】ここで、画像のベクトル量子化の過程にお
いて誤差拡散を適用することについて説明する。ベクト
ル量子化においては、前述のように、原画像における特
定のブロックの部分における平均階調と、ハーフトーン
画像におけるそのブロックの部分における平均階調との
差異、すなわち誤差pが存在する。この誤差拡散は、た
とえば、以降に処理するブロックを構成する画素値にす
べて誤差pを加算することによって行なうことができ
る。または、平均階調誤差を算出するときに誤差pによ
って補正することにより行なうことができる。
Here, application of error diffusion in the process of image vector quantization will be described. In the vector quantization, as described above, there is a difference between the average tone in a specific block portion in the original image and the average tone in the block portion in the halftone image, that is, an error p. This error diffusion can be performed, for example, by adding an error p to all pixel values constituting a block to be processed thereafter. Alternatively, the calculation can be performed by correcting the average tone error with the error p.

【0052】その平均階調誤差を算出するときの補正
は、具体的には、数1に示す評価関数における平均階調
誤差である|ΔD|の値を誤差pによって補正すること
により行なうことができる。この画像のベクトル量子化
の過程において誤差拡散を適用する場合の評価関数は、
たとえば、下記の数2のようになる。
More specifically, the correction for calculating the average tone error is performed by correcting the value of | ΔD |, which is the average tone error in the evaluation function shown in Expression 1, by the error p. it can. The evaluation function when error diffusion is applied in the process of vector quantization of this image is:
For example, Equation 2 below is obtained.

【数2】 F(Gi,Ho)=α×|ΔD−p|+β×ΔQ/L ただし、 F(Gi,Ho) ; 評価関数 Gi ; 原画像のベクトル(原画像のブロック) Ho ; コードワード α、β ; 係数(α+β=1、0≦α、β≦1) |ΔD|; 平均階調の誤差(|x|はxの絶対値を
表す) p ; 以前のブロックにおける平均階調の誤差 ΔQ ; ベクトル量子化の誤差 L ; 次元数(16)
F (Gi, Ho) = α × | ΔD−p | + β × ΔQ / L where F (Gi, Ho); evaluation function Gi; vector of original image (block of original image) Ho; codeword α, β; coefficient (α + β = 1, 0 ≦ α, β ≦ 1) | ΔD |; average grayscale error (| x | represents the absolute value of x) p; average grayscale error in the previous block ΔQ: error of vector quantization L: number of dimensions (16)

【0053】次に、ステップS8において、すべてのブ
ロックに対してステップS6とステップS7の処理を行
なったか否かが判定される。すべてのブロックに対して
その処理が済んでいない場合にはステップS6に戻り以
降のステップを繰り返す。すべてのブロックに対してそ
の処理が済んでいる場合にはステップS9に進む。次
に、ステップS9において、生成したハーフトーン画像
を保存する。
Next, in step S8, it is determined whether or not the processing in steps S6 and S7 has been performed on all blocks. If the processing has not been completed for all the blocks, the process returns to step S6 and the subsequent steps are repeated. If the processing has been completed for all blocks, the process proceeds to step S9. Next, in step S9, the generated halftone image is stored.

【0054】以上、本発明におけるハーフトーン処理に
ついて図1に基づいて一例を説明した。次に、その変形
例について説明する。前述の数1、数2に示す評価関数
において、前の項|ΔD|、|ΔD−p|は平均階調の
誤差であるから大まかな濃淡の再現性を表している。ま
た、後の項ΔQ/LはRMSE(root mean square err
or)であるから微細構造の再現性を表している。したが
って、αを大きくすると濃淡は良好に再現される代わり
に微細構造の再現が犠牲となる。逆にβを大きくすると
微細構造は良好に再現される代わりに濃淡の再現が犠牲
となる。そこで、適合度の算出が最適化されるように、
指示入力によって加重を与えるためのパラメータ(助変
数)としてのαとβを任意の値に設定できるように構成
すると好適である。総合的にみて良好な再現を得ること
ができるようにするためには、たとえば、α=β=0.
5とする。
An example of the halftone process according to the present invention has been described above with reference to FIG. Next, a modified example thereof will be described. In the evaluation functions shown in Expressions 1 and 2, the previous terms | ΔD | and | ΔD-p | represent errors in average gradation, and thus represent rough reproducibility of shading. The later term ΔQ / L is RMSE (root mean square err)
or) indicates the reproducibility of the microstructure. Therefore, when α is increased, the reproduction of the fine structure is sacrificed at the expense of good reproduction of the density. Conversely, when β is increased, the fine structure is not reproduced well, but the reproduction of light and shade is sacrificed. Therefore, to optimize the calculation of fitness,
It is preferable that α and β as parameters (auxiliary variables) for giving weight by instruction input can be set to arbitrary values. In order to obtain a good overall reproduction, for example, α = β = 0.
5 is assumed.

【0055】なお、本発明における評価関数は数1に示
す評価関数に限定されない。数1に示す評価関数におい
ては、平均階調が考慮されたが、ドットゲインやドット
の連続性を考慮に入れることもできる。
Note that the evaluation function in the present invention is not limited to the evaluation function shown in Expression 1. In the evaluation function shown in Expression 1, the average gradation is considered, but dot gain and dot continuity can be taken into consideration.

【0056】また、ステップS3のブロック分割におい
て、原画像を8画素×8画素のブロックに分割した。す
なわち、ブロックの形状は正方形である。本発明におけ
るブロック形状は正方形に限定されない。J×Kの矩形
であってもよい。また、本質的には平面を充填できれば
(擬似)六角形や(擬似)三角形でもよい。さらに、そ
れらを交互に組み合わせたタイリング(tiling) でも
よい。
In the block division in step S3, the original image is divided into blocks of 8 pixels × 8 pixels. That is, the shape of the block is a square. The block shape in the present invention is not limited to a square. It may be a J × K rectangle. In addition, a (pseudo) hexagon or (pseudo) triangle may be used as long as the plane can be essentially filled. Further, tiling which combines them alternately may be used.

【0057】また、前述のステップS6の説明におい
て、一例として全探索を行なうことを説明したが、本発
明における適合度の算出は全探索に限定されない。たと
えば、遺伝アルゴリズム(GA;genetic algorithm)
などの最適化手法を適用することができる。
Further, in the above description of step S6, it has been described that a full search is performed as an example. However, the calculation of the fitness in the present invention is not limited to the full search. For example, genetic algorithm (GA)
And other optimization techniques can be applied.

【0058】また、全探索による処理時間を短縮するた
めに、複数のコードワードとの比較を多数の処理部(C
PU;central processor unit、等)を用いて並列分散
処理することにより高速化することができる。また、原
画像に部分画像を設定し、その部分画像ごとにブロック
を複数の処理部の各々に割り当て、複数の処理部を用い
て並列分散処理することにより高速化することができ
る。また、評価関数の各項、すなわち前述の数1におい
ては、α×|ΔD|の項とβ×ΔQ/Lの項を、それぞ
れ別の処理部を用いて並列分散処理することにより高速
化することができる。
In order to reduce the processing time of the full search, a comparison with a plurality of codewords is performed by a number of processing units (C
(PU; central processor unit, etc.) to increase the processing speed. Further, the speed can be increased by setting a partial image in the original image, allocating a block to each of the plurality of processing units for each of the partial images, and performing parallel distributed processing using the plurality of processing units. In addition, in the respective terms of the evaluation function, that is, in the above-described equation 1, the α × | ΔD | and β × ΔQ / L terms are speeded up by performing parallel distributed processing using different processing units. be able to.

【0059】また、評価関数の演算を条件つきで途中で
打ち切ることにより処理を高速化することができる。た
とえば、α×|ΔD|を求めた時点でそれまでに求めら
れた最小値よりも大きな値となった場合には、明らかに
最小値にはなり得ないので次の評価に移ってよい。
Further, the processing can be speeded up by terminating the operation of the evaluation function conditionally on the way. For example, when α × │ΔD│ is found to be larger than the minimum value obtained so far, it cannot be the minimum value, and the next evaluation may be performed.

【0060】また、コードワードの平均階調値をあらか
じめ演算しておいてその順序にしたがって並べておけ
ば、この操作を簡略化することができる。また、前述の
ステップS6の説明において、コードワードを逐次生成
する説明を行なったが、コードワードをあらかじめ生成
しておき、その平均階調値とともに、コードブックとし
て記憶させておき、そのコードブックを参照して処理を
行なうと処理が高速化され好適である。
This operation can be simplified if the average gradation values of the code words are calculated in advance and arranged in the order. Further, in the above description of step S6, a description has been given of the method of sequentially generating codewords. However, a codeword is generated in advance and stored as a codebook together with its average gradation value. It is preferable to perform the processing by referring to the processing because the processing speed is increased.

【0061】また、図1に示した一例においては、原画
像が単色画像であるかのように説明した。多色画像を同
様に扱う場合には、画素の配列順または多色のプレーン
(レイヤー)順にベクトルを連接させればよい。たとえ
ば、Pを2以上の整数として原画像がP色画像であると
きには、ブロックの画素を(16×P)次元のベクトル
として処理することになる。
In the example shown in FIG. 1, description has been made as if the original image was a monochromatic image. When a multicolor image is handled in the same manner, the vectors may be connected in the pixel arrangement order or the multicolor plane (layer) order. For example, when P is an integer of 2 or more and the original image is a P-color image, the pixels of the block are processed as (16 × P) -dimensional vectors.

【0062】勿論、ベクトルを連接する代わりに多色画
像における各色画像の各々に対して、前述の処理を行な
うことができる。その場合において、各色画像の各々に
対する処理を複数の処理部によって並列分散処理するこ
とにより処理を高速化することができる。
Of course, the above-described processing can be performed on each of the color images in the multi-color image instead of connecting the vectors. In such a case, the processing can be speeded up by performing the parallel processing on each of the color images by a plurality of processing units.

【0063】以上、本発明について実施の形態により説
明を行なった。そこでは、ハーフトーン処理装置におけ
る処理を中心に説明した。しかし、本発明はハーフトー
ン処理装置に限定されるものではない。当然ながら、そ
の処理手順が記録された媒体も含まれるものである。
The present invention has been described with reference to the embodiment. Here, the processing in the halftone processing apparatus has been mainly described. However, the present invention is not limited to a halftone processing device. Of course, a medium on which the processing procedure is recorded is also included.

【0064】[0064]

【発明の効果】以上のとおりであるから、本発明の請求
項1に係るハーフトーン処理装置によれば、再現性の高
い、高品質のハーフトーン画像を得ることができるハー
フトーン処理装置が提供される。また本発明の請求項2
に係るハーフトーン処理装置によれば、平均階調誤差を
算出する過程でその誤差が相殺され誤差拡散における処
理の負荷を極めて小さくすることができる。また本発明
の請求項3に係るハーフトーン処理装置によれば、基本
的で最も簡明な処理であるとともに、水平線、垂直線の
再現性が優れる。また本発明の請求項4に係るハーフト
ーン処理装置によれば、周期性の認識が緩和され方向に
よる再現性が均等化される。また本発明の請求項5に係
るハーフトーン処理装置によれば、より再現性の高い、
より高品質のハーフトーン画像を得ることができる。ま
た本発明の請求項6に係るハーフトーン処理装置によれ
ば、コードワード演算にともなう処理時間を短縮でき、
コードワード数の適正減縮によりさらに処理時間を短縮
できる。また本発明の請求項7に係るハーフトーン処理
装置によれば、状況に応じて不必要なコードワードに係
わる処理を省略することができる。また本発明の請求項
8に係るハーフトーン処理装置によれば、ハーフトーン
処理を高速化することができる。また本発明の請求項9
に係るハーフトーン処理装置によれば、ハーフトーン処
理を高速化することができる。また本発明の請求項10
に係るハーフトーン処理装置によれば、ハーフトーン処
理を著しく高速化することができる。また本発明の請求
項11に係るハーフトーン処理装置によれば、ハーフト
ーン処理を高速化することができる。また本発明の請求
項12に係るハーフトーン処理装置によれば、多色画像
を対象とするハーフトーン処理を行なうことができる。
また本発明の請求項13に係るハーフトーン処理装置に
よれば、多色画像を対象とするハーフトーン処理を高速
化することができる。また本発明の請求項14に係るハ
ーフトーン処理装置によれば、多色画像を対象とするハ
ーフトーン処理を著しく高速化することができる。また
本発明の請求項15に係るハーフトーン処理手順が記録
された媒体によれば、再現性の高い、高品質のハーフト
ーン画像を得ることができるハーフトーン処理の処理手
順を記録した媒体が提供される。
As described above, the halftone processing apparatus according to the first aspect of the present invention provides a halftone processing apparatus capable of obtaining a high-quality halftone image with high reproducibility. Is done. Claim 2 of the present invention
According to the halftone processing device according to (1), the error is canceled in the process of calculating the average gradation error, and the processing load in error diffusion can be extremely reduced. Further, according to the halftone processing apparatus according to the third aspect of the present invention, the processing is basic and simplest, and the reproducibility of horizontal lines and vertical lines is excellent. Further, according to the halftone processing device of the fourth aspect of the present invention, recognition of periodicity is eased, and reproducibility in directions is equalized. Further, according to the halftone processing device according to claim 5 of the present invention, the reproducibility is higher,
A higher quality halftone image can be obtained. Further, according to the halftone processing device according to claim 6 of the present invention, the processing time involved in the codeword operation can be reduced,
The processing time can be further reduced by appropriately reducing the number of codewords. Further, according to the halftone processing device of the seventh aspect of the present invention, it is possible to omit processing relating to unnecessary codewords depending on the situation. Further, according to the halftone processing device of the present invention, the speed of the halftone processing can be increased. Claim 9 of the present invention
According to the halftone processing device according to the above, the speed of the halftone process can be increased. Claim 10 of the present invention
According to the halftone processing device according to the first aspect, the halftone processing can be remarkably speeded up. Further, according to the halftone processing device of the present invention, the speed of the halftone processing can be increased. Further, according to the halftone processing apparatus according to the twelfth aspect of the present invention, it is possible to perform halftone processing for a multicolor image.
Further, according to the halftone processing device according to the thirteenth aspect of the present invention, it is possible to speed up halftone processing for a multicolor image. Further, according to the halftone processing apparatus according to the fourteenth aspect of the present invention, it is possible to remarkably speed up halftone processing for a multicolor image. Further, according to the medium in which the halftone processing procedure according to claim 15 of the present invention is recorded, a medium in which the processing procedure of the halftone processing capable of obtaining a high-quality halftone image with high reproducibility is provided. Is done.

【図面の簡単な説明】[Brief description of the drawings]

【図1】本発明のハーフトーン処理装置における処理過
程の一例を示す図である。
FIG. 1 is a diagram showing an example of a processing process in a halftone processing device of the present invention.

【図2】原画像からブロックへの分割、および16次元
ベクトルの関係の説明図である。
FIG. 2 is an explanatory diagram of the division of an original image into blocks and the relationship between 16-dimensional vectors.

【図3】誤差拡散処理に関する説明図である。FIG. 3 is an explanatory diagram relating to an error diffusion process.

【図4】スカラー値における四捨五入の例と、最も簡単
な2次元のベクトル値における量子化の例を示す図であ
る。
FIG. 4 is a diagram illustrating an example of rounding in a scalar value and an example of quantization in a simplest two-dimensional vector value.

【図5】ベクトルが4次元の場合の量子化代表値ベクト
ルを示す図である。
FIG. 5 is a diagram showing a quantized representative value vector when the vector is four-dimensional.

【図6】周期的な網点を刷り重ねるとき発生する亀甲模
様(ロゼッタ・パターン)を示す図である。
FIG. 6 is a diagram showing a turtle pattern (rosetta pattern) generated when a periodic halftone dot is overprinted.

【図7】ハーフトーン画像の一例(Aはコンタクトスク
リーン、Bは周期的分散網点による画像)を示す図であ
る。
FIG. 7 is a diagram showing an example of a halftone image (A is a contact screen, B is an image based on a periodically dispersed halftone dot).

【図8】ハーフトーン画像の一例(Cは誤差拡散法、D
はブルーノイズマスク法による画像)を示す図である。
FIG. 8 shows an example of a halftone image (C is an error diffusion method,
FIG. 3 is a diagram showing an image by a blue noise mask method.

Claims (15)

【特許請求の範囲】[Claims] 【請求項1】L,M,Nを正の整数として、Mビット表
現の画素値を有する画像から、M>NであるNビット表
現の画素値を有する画像を導出するハーフトーン処理装
置であって、ブロック分割手段と、適合度算出手段と、
画像生成手段と、誤差拡散手段とを有し、 前記ブロック分割手段は、入力画像を所定サイズLの画
素からなるブロックに分割し、 前記適合度算出手段は、前記分割されたブロックの画素
をL次元Mビットのベクトルとし、これを変換候補とな
るL次元Nビットのコードワードと比較し適合度を算出
し、 前記画像生成手段は適合度が最大のコードワードに基づ
いて前記ブロックにL次元Nビットの画像を生成し、 前記誤差拡散手段は、前記ブロックの平均階調誤差を以
降に処理する前記ブロックの近傍ブロックへ誤差拡散す
る、 ことを特徴とするハーフトーン処理装置。
1. A halftone processing device which derives an image having pixel values of N bits where M> N from an image having pixel values of M bits, where L, M and N are positive integers. A block dividing unit, a fitness calculating unit,
An image generating unit; and an error diffusion unit. The block dividing unit divides the input image into blocks each having a predetermined size L. The fitness calculating unit divides the pixels of the divided blocks into L. A vector of M dimensions is obtained, and this is compared with an L-dimensional N-bit code word which is a conversion candidate to calculate the degree of conformity. The image generating means assigns L-dimensional N to the block based on the code word having the maximum degree of conformity. A halftone processing apparatus, wherein a bit image is generated, and the error diffusion unit performs error diffusion on an average gradation error of the block to a block adjacent to the block to be processed thereafter.
【請求項2】請求項1記載のハーフトーン処理装置にお
いて、前記誤差拡散手段が行なう前記誤差拡散は、平均
階調誤差を算出する過程でその誤差を相殺することで実
行することを特徴とするハーフトーン処理装置。
2. A halftone processing apparatus according to claim 1, wherein said error diffusion performed by said error diffusion means is executed by canceling the average gradation error in the process of calculating the error. Halftone processing device.
【請求項3】請求項1または2記載のハーフトーン処理
装置において、J,Kを正の整数として、前記ブロック
分割手段は、入力画像を所定サイズであるJ×K=Lの
画素からなる矩形のブロックに分割することを特徴とす
るハーフトーン処理装置。
3. The halftone processing apparatus according to claim 1, wherein J and K are positive integers, and said block dividing means divides the input image into rectangles each having a predetermined size of J × K = L pixels. A halftone processing device, wherein the halftone processing device is divided into blocks.
【請求項4】請求項1〜3のいずれかに記載のハーフト
ーン処理装置において、前記ブロック分割手段は、入力
画像を所定サイズLの画素からなる、六角形または三角
形、あるいはそれらを組み合わせたタイリング(tilin
g)形のブロックに分割することを特徴とするハーフト
ーン処理装置。
4. The halftone processing device according to claim 1, wherein said block dividing means converts the input image into hexagons or triangles consisting of pixels of a predetermined size L, or a combination of these. Ring (tilin
g) A halftone processing device, which is divided into shaped blocks.
【請求項5】請求項1〜4のいずれかに4記載のハーフ
トーン処理装置において、適合度を算出するための最適
化パラメータを設定する最適化パラメータ設定手段を有
することを特徴とするハーフトーン処理装置。
5. The halftone processing device according to claim 4, further comprising an optimization parameter setting means for setting an optimization parameter for calculating a fitness. Processing equipment.
【請求項6】請求項1〜5のいずれかに記載のハーフト
ーン処理装置において、あらかじめ生成したL次元Nビ
ットのコードワードが記載されたコードブックを記憶す
るコードブック記憶手段を有することを特徴とするハー
フトーン処理装置。
6. A halftone processing apparatus according to claim 1, further comprising a codebook storing means for storing a codebook in which a previously generated L-dimensional N-bit codeword is described. And a halftone processing device.
【請求項7】請求項1〜5のいずれかに記載のハーフト
ーン処理装置において、L次元Nビットのコードワード
を逐次生成するコードワード逐次生成手段を有すること
を特徴とするハーフトーン処理装置。
7. The halftone processing apparatus according to claim 1, further comprising a codeword sequential generation means for sequentially generating an L-dimensional N-bit codeword.
【請求項8】請求項1〜7のいずれかに記載のハーフト
ーン処理装置において、前記L次元Nビットのコードワ
ードは、その総数を2^L(2のL乗)未満とすること
を特徴とするハーフトーン処理装置。
8. The halftone processing apparatus according to claim 1, wherein the total number of said L-dimensional N-bit code words is less than 2 ^ L (2 to the power of L). And a halftone processing device.
【請求項9】請求項1〜8のいずれかに記載のハーフト
ーン処理装置において、そのハーフトーン処理装置は複
数の処理部を具備し、前記ブロック分割手段が分割した
ブロックを前記複数の処理部の各々に割り当て、前記複
数の処理部によって並列分散処理することを特徴とする
ハーフトーン処理装置。
9. The halftone processing device according to claim 1, wherein said halftone processing device includes a plurality of processing units, and the block divided by said block dividing means is divided into said plurality of processing units. Wherein the plurality of processing units perform parallel distributed processing.
【請求項10】請求項9記載のハーフトーン処理装置に
おいて、前記適合度算出手段が行なう前記適合度の算出
を前記複数の処理部によって並列分散処理することを特
徴とするハーフトーン処理装置。
10. The halftone processing apparatus according to claim 9, wherein the calculation of the fitness performed by the fitness calculating means is performed in parallel by the plurality of processing units.
【請求項11】請求項1〜10のいずれかに記載のハー
フトーン処理装置において、前記適合度算出手段による
適合度の算出は、前記コードワードを全探索するのでは
なく、適合度が最大のコードワードを最適化手法により
探索することを特徴とするハーフトーン処理装置。
11. The halftone processing apparatus according to claim 1, wherein the calculation of the fitness by the fitness calculating means does not perform a full search of the codeword, but calculates the maximum fitness. A halftone processing device for searching for a codeword by an optimization technique.
【請求項12】請求項1〜11のいずれかに記載のハー
フトーン処理装置において、Pを2以上の整数として前
記画像はP色画像であって、前記ブロックの画素を(L
×P)次元Mビットのベクトルとして処理することを特
徴とするハーフトーン処理装置。
12. The halftone processing apparatus according to claim 1, wherein P is an integer of 2 or more, the image is a P-color image, and pixels of the block are (L
× P) A halftone processing device for processing as a vector of dimension M bits.
【請求項13】請求項1〜11のいずれかに記載のハー
フトーン処理装置において、前記画像は多色画像であっ
て、その多色画像の各色画像ごとに前記処理を行なうこ
とを特徴とするハーフトーン処理装置。
13. The halftone processing apparatus according to claim 1, wherein said image is a multicolor image, and said processing is performed for each color image of said multicolor image. Halftone processing device.
【請求項14】請求項13記載のハーフトーン処理装置
において、前記各色画像ごとの処理を前記複数の処理部
によって前記各色画像ごとに並列分散処理することを特
徴とするハーフトーン処理装置。
14. The halftone processing apparatus according to claim 13, wherein the processing for each color image is performed by the plurality of processing units in parallel and distributed processing for each color image.
【請求項15】L,M,Nを正の整数として、Mビット
表現の画素値を有する画像から、M>NであるNビット
表現の画素値を有する画像を導出するハーフトーン処理
手順が記録された媒体であって、前記ハーフトーン処理
手順は、ブロック分割過程と、適合度算出過程と、画像
生成過程、誤差拡散処理過程とを有し、 前記ブロック分割過程において、入力画像を所定サイズ
Lの画素からなるブロックに分割し、 前記適合度算出過程において、前記分割されたブロック
の画素をL次元Mビットのベクトルとし、これを変換候
補となるL次元Nビットのコードワードと比較し適合度
を算出し、 前記画像再生過程において、適合度が最大のコードワー
ドに基づいて前記ブロックにL次元Nビットの画像を生
成し、 前記誤差拡散過程において、前記ブロックの平均階調誤
差を以降に処理する前記ブロックの近傍ブロックへ誤差
拡散する、 ことを特徴とするハーフトーン処理手順が記録された媒
体。
15. A halftone processing procedure for deriving an image having a pixel value of an N-bit expression where M> N from an image having a pixel value of an M-bit value, where L, M, and N are positive integers, is recorded. The halftone processing procedure includes a block division process, a fitness calculation process, an image generation process, and an error diffusion process process. In the block division process, the input image is converted to a predetermined size L. In the fitness calculation step, the pixels of the divided block are defined as an L-dimensional M-bit vector, which is compared with an L-dimensional N-bit code word which is a conversion candidate. In the image reproduction process, an L-dimensional N-bit image is generated in the block based on the codeword having the highest fitness, and in the error diffusion process, A medium in which a halftone processing procedure is recorded, wherein error diffusion of an average gradation error of the block to a block adjacent to the block to be processed thereafter is performed.
JP2000212464A 2000-07-13 2000-07-13 Half tone processor, and medium having stored its procedure Pending JP2002027249A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2000212464A JP2002027249A (en) 2000-07-13 2000-07-13 Half tone processor, and medium having stored its procedure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2000212464A JP2002027249A (en) 2000-07-13 2000-07-13 Half tone processor, and medium having stored its procedure

Publications (1)

Publication Number Publication Date
JP2002027249A true JP2002027249A (en) 2002-01-25

Family

ID=18708384

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2000212464A Pending JP2002027249A (en) 2000-07-13 2000-07-13 Half tone processor, and medium having stored its procedure

Country Status (1)

Country Link
JP (1) JP2002027249A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006071877A (en) * 2004-09-01 2006-03-16 Dainippon Printing Co Ltd Method for producing screen image
JP2006295899A (en) * 2005-03-18 2006-10-26 Fuji Xerox Co Ltd Dispersion evaluation apparatus, layout determination apparatus, dispersion evaluation method, layout determination method and program therefor
CN100336077C (en) * 2003-11-21 2007-09-05 财团法人工业技术研究院 Internal storage management method of error diffusion method and its half colour tint treatment module

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100336077C (en) * 2003-11-21 2007-09-05 财团法人工业技术研究院 Internal storage management method of error diffusion method and its half colour tint treatment module
JP2006071877A (en) * 2004-09-01 2006-03-16 Dainippon Printing Co Ltd Method for producing screen image
JP4591951B2 (en) * 2004-09-01 2010-12-01 大日本印刷株式会社 Screen image generation method
JP2006295899A (en) * 2005-03-18 2006-10-26 Fuji Xerox Co Ltd Dispersion evaluation apparatus, layout determination apparatus, dispersion evaluation method, layout determination method and program therefor
JP4721112B2 (en) * 2005-03-18 2011-07-13 富士ゼロックス株式会社 Distributed evaluation apparatus, arrangement determination apparatus, and program thereof

Similar Documents

Publication Publication Date Title
US8406295B2 (en) Image processing device and image processing method
US5394250A (en) Image processing capable of handling multi-level image data without deterioration of image quality in highlight areas
US6552822B1 (en) Image processing method and apparatus
JP3571716B2 (en) Improved error variance system
JPH05308514A (en) Method for quantizing pixel value
US6330075B1 (en) Image processing method and apparatus
JP3973734B2 (en) Electronic image processing system and processing method
EP0591977A1 (en) Image processor and printing apparatus
JPH10112801A (en) Input image processing method
US7224487B2 (en) Ink reduction error diffusion
JPH11234513A (en) Pseudo gradation image processing unit
JPH04286465A (en) Picture processor
JPH09284564A (en) Screened multilevel image compression method
JP2002027249A (en) Half tone processor, and medium having stored its procedure
JPH0324673A (en) Method for processing image data
JP2002027250A (en) Half tone processor, and medium having stored its procedure
JP2003110849A (en) Multilevel halftoning method of digital image
US7149351B2 (en) Image processing apparatus and method, computer program and computer-readable storage medium
JP3728025B2 (en) Image processing method and apparatus
US6011907A (en) Method of and system for promoting predetermined patterns at a desirable orientation in an output image
US6006011A (en) Target patterns controlled error management
JP2662405B2 (en) Image processing device
JP2675792B2 (en) Image processing device
JP4027300B2 (en) Image processing method, image processing apparatus, image forming apparatus, computer program, and recording medium
JP3015049B2 (en) Image data processing method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20070427

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20080529

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20080603

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20081007