JPH114438A - Image compressor and camera - Google Patents

Image compressor and camera

Info

Publication number
JPH114438A
JPH114438A JP9172918A JP17291897A JPH114438A JP H114438 A JPH114438 A JP H114438A JP 9172918 A JP9172918 A JP 9172918A JP 17291897 A JP17291897 A JP 17291897A JP H114438 A JPH114438 A JP H114438A
Authority
JP
Japan
Prior art keywords
quantization
quantization table
data
image
error
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP9172918A
Other languages
Japanese (ja)
Other versions
JP3645690B2 (en
Inventor
Nobusato 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.)
Pentax Corp
Original Assignee
Asahi Kogaku Kogyo 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 Asahi Kogaku Kogyo Co Ltd filed Critical Asahi Kogaku Kogyo Co Ltd
Priority to JP17291897A priority Critical patent/JP3645690B2/en
Publication of JPH114438A publication Critical patent/JPH114438A/en
Application granted granted Critical
Publication of JP3645690B2 publication Critical patent/JP3645690B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Television Signal Processing For Recording (AREA)
  • Compression Or Coding Systems Of Tv Signals (AREA)

Abstract

PROBLEM TO BE SOLVED: To compress an image by predicting a quantization error, so as to reduce an image deterioration due to quantization and inverse quantization. SOLUTION: This image compressor 10 is provided with a quantization table generating circuit 20. The quantization table generating circuit 20 is provided with a statistical quantity storage section 21, a quantization error predict section 22, and a quantization table calculation section 24. A DCT processing circuit 12 applies discrete cosine transform(DCT) to image data, that is, a luminance signal and color difference signals Cb, Cr to convert the data into a DCT transform coefficient. The statistic quantity storage section 21 calculates sampling number data which is a statistical quantity with a close relation to a quantization error from the DCT transform coefficient for each spatial frequency, stores sampling number data of a prescribed number, and updates the data for each image input. The quantization error predict section 22 uses the sampling number data to conduct quantization error prediction. The quantization table calculation section 24 calculates quantization tables Qya, Qca, based on the quantization error predict result by the quantization error predict section 22, and outputs them to a quantization-processing circuit 14.

Description

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

【0001】[0001]

【発明の属する技術分野】本発明は、カラー静止画像を
JPEGアルゴリズムに準拠して情報圧縮する画像圧縮
装置およびカメラに関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an image compression apparatus and a camera for compressing information of a color still image according to the JPEG algorithm.

【0002】[0002]

【従来の技術】高解像度画像を符号化(圧縮)して、情
報の授受を通信伝送路を介して行なう標準化アルゴリズ
ムがJPEG(Joint Photographic Expert Group )か
ら勧告されている。このJPEGアルゴリズムのプロセ
スでは、原画像データを2次元離散コサイン変換(以
下、2次元DCTという)によって空間周波数軸上の成
分に分解し、この空間周波数軸上で表された各データを
量子化し、量子化した各データを符号化することによ
り、大幅な情報圧縮が行なわれている。また圧縮された
画像データは復号化、逆量子化することにより伸張され
る。
2. Description of the Related Art JPEG (Joint Photographic Expert Group) has recommended a standardized algorithm for encoding (compressing) a high-resolution image and transmitting and receiving information via a communication transmission line. In the JPEG algorithm process, original image data is decomposed into components on a spatial frequency axis by a two-dimensional discrete cosine transform (hereinafter, referred to as a two-dimensional DCT), and each data represented on the spatial frequency axis is quantized. Significant information compression is performed by encoding each quantized data. The compressed image data is expanded by decoding and inverse quantization.

【0003】JPEGでは、量子化あるいは逆量子化の
ために、所定の量子化テーブルを推奨している。このD
CTと量子化による画像圧縮は、広範な被写体において
効率の良い画像圧縮であるが、原画像データと伸張され
た再生画像データとの間に誤差を伴う非可逆方式であ
る。
[0003] JPEG recommends a predetermined quantization table for quantization or inverse quantization. This D
Image compression by CT and quantization is an efficient image compression for a wide range of subjects, but is an irreversible method involving an error between the original image data and the expanded reproduced image data.

【0004】[0004]

【発明が解決しようとする課題】画像圧縮の量子化ある
いは逆量子化において、被写体の多様性にも関わらず、
単一の量子化テーブルを使用するため、被写体によって
は再生画像データの劣化が大きくなることが問題であ
る。
In the quantization or inverse quantization of image compression, regardless of the variety of subjects,
Since a single quantization table is used, there is a problem in that reproduced image data greatly deteriorates depending on a subject.

【0005】本発明は、この様な問題点に鑑み、量子化
誤差を予測し、量子化と逆量子化による画像劣化が少な
い画像圧縮を容易にする画像圧縮装置およびカメラを提
供することが目的である。
SUMMARY OF THE INVENTION In view of the above problems, an object of the present invention is to provide an image compression apparatus and a camera which predict a quantization error and facilitates image compression with less image deterioration due to quantization and inverse quantization. It is.

【0006】[0006]

【課題を解決するための手段】本発明による画像圧縮装
置は、撮影光学系から得られた原画像データに対応した
直交変換係数を、量子化テーブルを用いて量子化し、量
子化直交変換係数を求める量子化手段と、直交変換係数
に基づいて、量子化テーブルを更新する量子化テーブル
作成手段とを備えたことを特徴としている。
An image compression apparatus according to the present invention quantizes an orthogonal transform coefficient corresponding to original image data obtained from a photographing optical system using a quantization table, and quantizes the quantized orthogonal transform coefficient. It is characterized by comprising a quantization means to be obtained and a quantization table creation means for updating the quantization table based on the orthogonal transform coefficient.

【0007】画像圧縮装置において、好ましくは、量子
化テーブル作成手段が、量子化直交変換係数を逆量子化
して求められる逆量子化直交変換係数と直交変換係数と
の誤差の絶対値が2よりも大きい直交変換係数の標本数
を、直交変換係数毎にカウントし、各画像における標本
数のデータ(標本数データ)を、所定画像枚数分だけ格
納する格納部を備える。
In the image compression apparatus, preferably, the quantization table creating means has an absolute value of an error between the inversely quantized orthogonal transform coefficient and the orthogonal transform coefficient obtained by inversely quantizing the quantized orthogonal transform coefficient, which is larger than two. A storage unit is provided which counts the number of samples of a large orthogonal transform coefficient for each orthogonal transform coefficient, and stores data of the number of samples (sample number data) in each image for a predetermined number of images.

【0008】画像圧縮装置において、好ましくは、格納
部において、新たな画像の標本数データが入力される毎
に、格納部に格納された最も古い画像の標本数データが
消去されることにより、所定枚数の標本数データが更新
される。
In the image compression apparatus, preferably, the sample number data of the oldest image stored in the storage unit is deleted every time the sample number data of a new image is input to the storage unit. The number of samples data is updated.

【0009】画像圧縮装置において、好ましくは、量子
化テーブル作成手段が、所定枚数の標本数データに基づ
いて、誤差を予測する量子化誤差予測部と、量子化誤差
予測部から得られた予測誤差値に基づいて量子化テーブ
ルを算出し、量子化手段に出力する量子化テーブル算出
部とを備える。
In the image compression apparatus, preferably, the quantization table creating means includes a quantization error prediction section for predicting an error based on a predetermined number of sample data, and a prediction error obtained from the quantization error prediction section. A quantization table calculation unit that calculates a quantization table based on the value and outputs the quantization table to a quantization unit.

【0010】画像圧縮装置において、好ましくは、量子
化テーブル作成手段が、新規に入力された標本数データ
を含めた所定枚数の標本数データにより、量子化テーブ
ルを作成する。また好ましくは、量子化テーブル作成手
段が、新規に入力された標本数データより前に入力され
た所定枚数の標本数データにより、量子化テーブルを作
成する。
In the image compression apparatus, preferably, the quantization table creating means creates a quantization table based on a predetermined number of sample number data including the newly input sample number data. Preferably, the quantization table creating means creates the quantization table based on a predetermined number of sample number data input before the newly input sample number data.

【0011】また本発明によるカメラは、撮影光学系か
ら得られた原画像データに対応した直交変換係数を、量
子化テーブルを用いて量子化し、量子化直交変換係数を
求める量子化手段と、直交変換係数に基づいて、量子化
テーブルを更新する量子化テーブル作成手段とを備え、
量子化テーブル作成手段が、1画像分の画像データに対
応し、量子化直交変換係数を逆量子化して求められる逆
量子化直交変換係数と直交変換係数との誤差の絶対値が
2よりも大きい直交変換係数の標本数を、直交変換係数
毎にカウントし、各画像における標本数のデータ(標本
数データ)を、所定画像枚数分だけ格納可能な標本数デ
ータ格納手段と、標本数データ格納手段に格納されてい
る所定枚数の標本数データに基づいて、誤差を予測する
量子化誤差予測部と、量子化誤差予測部から得られた予
測誤差値に基づいて量子化テーブルを算出し、量子化手
段に出力する量子化テーブル算出部とを備え、さらに量
子化テーブル作成手段によって作成された量子化テーブ
ルを用いて画像圧縮された圧縮画像データを記録媒体に
記録する記録手段と、作成された量子化テーブルを格納
する量子化テーブル格納手段とを備える。
Further, the camera according to the present invention quantizes orthogonal transform coefficients corresponding to original image data obtained from the photographing optical system by using a quantization table to obtain a quantized orthogonal transform coefficient. A quantization table creating unit that updates the quantization table based on the transform coefficient,
The quantization table creation means corresponds to the image data of one image, and the absolute value of the error between the inversely quantized orthogonal transform coefficient and the orthogonal transform coefficient obtained by inversely quantizing the quantized orthogonal transform coefficient is greater than 2. A sample number data storage unit that counts the number of samples of the orthogonal transform coefficient for each orthogonal transform coefficient and stores data of the number of samples in each image (sample number data) for a predetermined number of images; A quantization error prediction unit that predicts an error based on a predetermined number of sample number data stored in the quantization error prediction unit, and calculates a quantization table based on a prediction error value obtained from the quantization error prediction unit. Recording means for recording on a recording medium compressed image data compressed using the quantization table created by the quantization table creating means , And a quantization table storage means for storing the quantization tables created.

【0012】カメラにおいて、好ましくは、記録手段に
よって、圧縮画像データが記録媒体の第1の領域に記録
され、量子化に用いられた量子化テーブルが記録媒体の
第2の領域に記録される。
[0012] In the camera, preferably, the recording means records the compressed image data in a first area of the recording medium, and a quantization table used for quantization is recorded in a second area of the recording medium.

【0013】カメラにおいて、好ましくは、量子化テー
ブル作成手段が、テーブル格納手段に格納されている量
子化テーブルを、量子化に用いるとともに、記録媒体へ
の圧縮画像データの記録が完了した後に、この圧縮画像
データと対応する標本数データに基づいた量子化テーブ
ルを作成し、作成された量子化テーブルを量子化テーブ
ル格納手段に格納することにより、量子化テーブルを更
新する。
[0013] In the camera, preferably, the quantization table creating means uses the quantization table stored in the table storage means for quantization, and after the recording of the compressed image data on the recording medium is completed, A quantization table is created based on the compressed image data and the corresponding sample number data, and the quantization table is updated by storing the created quantization table in the quantization table storage means.

【0014】[0014]

【発明の実施の形態】以下、本発明による画像圧縮装置
およびカメラの実施形態について添付図面を参照して説
明する。
DESCRIPTION OF THE PREFERRED EMBODIMENTS Embodiments of an image compression apparatus and a camera according to the present invention will be described below with reference to the accompanying drawings.

【0015】図1には本発明の実施形態であるカメラの
概略構成が示される。被写体Sの撮影光学系11を介し
て得られる被写体像(静止画)は、色分解光学系13に
よってレッド(R)、グリーン(G)およびブルー
(B)の画像に分解され、例えばCCDから成る3枚の
撮像素子15に結像される。即ち、各撮像素子15には
それぞれR、G、Bの画像が形成される。これらR、
G、Bの画像信号は信号処理部17において所定の処理
を施され、輝度データY、色差データCb、Crに変換
されて、画像メモリ19に入力される。画像メモリ19
は輝度データYおよび色差データCb、Crをそれぞれ
格納するために、相互に独立したメモリ領域に分割され
ており、各メモリ領域は1画像分の記憶容量を有してい
る。輝度データY、色差データCb、Crは、画像圧縮
装置10に対する入力データ(原画像データ)である。
FIG. 1 shows a schematic configuration of a camera according to an embodiment of the present invention. A subject image (still image) of the subject S obtained through the photographing optical system 11 is separated into red (R), green (G), and blue (B) images by the color separation optical system 13, and is composed of, for example, a CCD. An image is formed on three image sensors 15. That is, R, G, and B images are formed on the respective image sensors 15. These R,
The G and B image signals are subjected to predetermined processing in a signal processing unit 17, converted into luminance data Y and color difference data Cb and Cr, and input to an image memory 19. Image memory 19
Are divided into mutually independent memory areas for storing the luminance data Y and the color difference data Cb and Cr, respectively, and each memory area has a storage capacity for one image. The luminance data Y and the color difference data Cb and Cr are input data (original image data) to the image compression device 10.

【0016】輝度データY、あるいは色差データCb、
Crは、DCT処理回路12において2次元DCTを施
され、空間周波数毎にDCT係数に変換される。これら
のDCT係数に基づいて、量子化テーブル作成回路20
により輝度データ用量子化テーブルQya、および色差
データ用量子化テーブルQcaが作成され、量子化処理
回路14に出力される。量子化テーブル作成回路20に
は統計量格納部21、量子化誤差予測部22、量子化テ
ーブル算出部24が設けられているが、これらの構成お
よび動作については後で詳述する。
The luminance data Y or the color difference data Cb,
Cr is subjected to two-dimensional DCT in the DCT processing circuit 12, and is converted into DCT coefficients for each spatial frequency. Based on these DCT coefficients, the quantization table creation circuit 20
Thus, a quantization table Qya for luminance data and a quantization table Qca for color difference data are created, and are output to the quantization processing circuit 14. The quantization table creation circuit 20 is provided with a statistic storage unit 21, a quantization error prediction unit 22, and a quantization table calculation unit 24, and the configuration and operation thereof will be described later in detail.

【0017】輝度データY、あるいは色差データCb、
CrのDCT係数は、量子化処理回路14において、量
子化テーブル作成回路20により作成された量子化テー
ブルQya、Qcaを用いて量子化され、量子化DCT
係数に変換される。輝度データY、あるいは色差データ
Cb、Crの量子化DCT係数は、図示しない符号化処
理回路によってJPEGアルゴリズムに従って符号化
(圧縮)され、圧縮画像データとして記録媒体Mの圧縮
画像データ記録領域M1に記録される。量子化処理回路
14において用いられた量子化テーブルQya、Qca
は、記録媒体Mのテーブル記録領域M2に記録される。
記録媒体Mは、例えば着脱自在なICカードである。
The luminance data Y or the color difference data Cb,
The DCT coefficient of Cr is quantized by the quantization processing circuit 14 using the quantization tables Qya and Qca created by the quantization table creation circuit 20, and the quantized DCT
Converted to coefficients. The quantized DCT coefficients of the luminance data Y or the color difference data Cb and Cr are encoded (compressed) by an encoding processing circuit (not shown) according to the JPEG algorithm, and are recorded as compressed image data in the compressed image data recording area M1 of the recording medium M. Is done. The quantization tables Qya and Qca used in the quantization processing circuit 14
Is recorded in the table recording area M2 of the recording medium M.
The recording medium M is, for example, a removable IC card.

【0018】図2は、JPEGの推奨する従来の量子化
テーブルを示す。図2(a)は輝度値Y用の量子化テー
ブルQyであり、図2(b)は色差Cb、Cr用の量子
化テーブルQcである。
FIG. 2 shows a conventional quantization table recommended by JPEG. FIG. 2A shows a quantization table Qy for the luminance value Y, and FIG. 2B shows a quantization table Qc for the color differences Cb and Cr.

【0019】図3を参照して、量子化と逆量子化につい
て説明する。図3には、一例として、8×8画素の原画
像データ、すなわち輝度データYのDCT係数Fvuと、
量子化DCT係数Rvuと、逆量子化DCT係数F’vu
および量子化テーブルQyが示される。添字v,uは、
64個のDCT係数を8×8のマトリクスの形式で表示
したときの縦および横方向の位置をそれぞれ示し、添字
vは上から0,1,2,...7であり、添字uは左か
ら0,1,2,...7である。
Referring to FIG. 3, quantization and inverse quantization will be described. FIG. 3 shows, as an example, original image data of 8 × 8 pixels, that is, a DCT coefficient F vu of luminance data Y,
A quantized DCT coefficient R vu and an inversely quantized DCT coefficient F ′ vu ,
And a quantization table Qy. The subscripts v and u are
The vertical and horizontal positions when the 64 DCT coefficients are displayed in the form of an 8 × 8 matrix are shown, and the subscript v is 0, 1, 2,. . . 7 and the subscript u is 0, 1, 2,. . . 7

【0020】原画像データは、DCT処理回路12にお
ける2次元DCTによって、8×8=64個のDCT係
数Fvuに変換される。2次元DCTは公知であるため、
ここでは詳述しない。
The original image data is converted into 8 × 8 = 64 DCT coefficients F vu by two-dimensional DCT in the DCT processing circuit 12. Since two-dimensional DCT is known,
It will not be described in detail here.

【0021】64個のDCT係数のうち、位置(0,
0)にあるDCT成分F00はDC(直流)成分であり、
残り63個のDCT係数FvuはAC(交流)成分であ
る。AC成分は、係数F01もしくは係数F10から係数F
77に向かって、より高い空間周波数成分が8×8画素ブ
ロックの原画像データ中にどのくらいあるかを示してい
る。DC成分は8×8画素のブロック全体の画素値の平
均値(直流成分)を表している。すなわち、各DCT係
数Fvuはそれぞれ所定の空間周波数に対応している。
Of the 64 DCT coefficients, the position (0,
0) is a DC (direct current) component,
The remaining 63 DCT coefficients F vu are AC (alternating current) components. The AC component is calculated from the coefficient F01 or the coefficient F10 to the coefficient F
Towards 77, it shows how much higher spatial frequency components are in the original image data of the 8 × 8 pixel block. The DC component represents the average value (DC component) of the pixel values of the entire block of 8 × 8 pixels. That is, each DCT coefficient F vu corresponds to a predetermined spatial frequency.

【0022】量子化テーブルQvuを用いてDCT係数F
vuを量子化する式は(1)式により定義される。この式
におけるround は、最も近い整数への近似を意味する。
すなわち、DCT係数Fvuおよび量子化テーブルQvu
各要素同士の割算と、四捨五入とによって、量子化DC
T係数Rvuが求められる。
Using a quantization table Q vu , a DCT coefficient F
The equation for quantizing vu is defined by equation (1). Round in this equation means approximation to the nearest integer.
That is, the division of each element of the DCT coefficient F vu and the quantization table Q vu and the rounding-off result in the quantization DC
The T coefficient R vu is determined.

【0023】 Rvu=round (Fvu/Qvu) ・・・(1) {0≦u,v ≦7}R vu = round (F vu / Q vu ) (1) {0 ≦ u, v ≦ 7}

【0024】このように、量子化処理回路14において
求められた量子化DCT係数Rvuは、例えばJPEGに
準拠したハフマン符号化を用いて、DC成分、AC成分
毎に符号化され、圧縮画像データとして記録媒体Mに記
録される。ハフマン符号化については、従来公知のため
詳細な説明は省略する。
As described above, the quantized DCT coefficient R vu obtained by the quantization processing circuit 14 is encoded for each DC component and AC component using, for example, Huffman coding conforming to JPEG, and compressed image data Is recorded on the recording medium M. Since the Huffman coding is conventionally known, a detailed description is omitted.

【0025】符号化された圧縮画像信号を伸張して、画
面に表示するためには、復号化、逆量子化、2次元DC
Tの逆変換(以下、2次元IDCTという)の処理が必
要である。この復号化は、ハフマン符号化とは逆の作用
であり、従来公知であるため詳述しない。復号化によっ
て得られた量子化DCT係数は、量子化に用いた量子化
テーブルQy、Qcを用いてそれぞれ逆量子化され、逆
量子化DCT係数F’vuに変換される。これらの逆量子
化DCT係数F’vuは2次元DCTの逆変換(以下、2
次元IDCTという)を施され、それぞれ輝度データ
Y’、色差データCb’、Cr’に変換される。2次元
IDCTについても公知であるのでここでは詳述しな
い。
In order to decompress the coded compressed image signal and display it on the screen, decoding, inverse quantization, two-dimensional DC
Processing of the inverse transformation of T (hereinafter referred to as two-dimensional IDCT) is required. This decoding is the reverse of the operation of Huffman coding, and is well known in the art and will not be described in detail. The quantized DCT coefficients obtained by the decoding are inversely quantized using the quantization tables Qy and Qc used for quantization, and are converted into inversely quantized DCT coefficients F ′ vu . These inversely quantized DCT coefficients F ′ vu are inversely transformed by two-dimensional DCT (hereinafter, 2
(Referred to as dimension IDCT), and converted into luminance data Y ′ and color difference data Cb ′ and Cr ′, respectively. Since the two-dimensional IDCT is also known, it will not be described in detail here.

【0026】例えば、図3に示すDCT係数F00(=2
61)を、量子化係数q00(=16)を用いて(1)式
により求めると、量子化DCT係数R00は16になり、
この量子化DCT係数R00(=16)にq00(=16)
を掛け合わせると、逆量子化DCT係数F’00(=16 ×
16=256)が求められる。なお、DCT係数F00(=26
1)と、逆量子化DCT係数F’00(=256)との差
分値i00(256-261=-5)を係数値誤差とする。このよう
に量子化では除算の余りを丸めるため、圧縮画像データ
の伸張(逆量子化)において係数値誤差iを生じる。こ
れが量子化誤差の原因である。
For example, the DCT coefficient F00 (= 2
61) is obtained by Expression (1) using the quantization coefficient q00 (= 16), the quantization DCT coefficient R00 becomes 16, and
Q00 (= 16) is added to the quantized DCT coefficient R00 (= 16).
, The inverse quantized DCT coefficient F'00 (= 16 ×
16 = 256). The DCT coefficient F00 (= 26
The difference value i00 (256-261 = -5) between 1) and the inversely quantized DCT coefficient F'00 (= 256) is defined as a coefficient value error. In this way, in the quantization, the remainder of the division is rounded, so that a coefficient value error i occurs in the expansion (inverse quantization) of the compressed image data. This is the cause of the quantization error.

【0027】図4は、DCT係数Fvuと、量子化および
逆量子化した後の逆量子化DCT係数F’vuとの各係数
毎の係数値誤差ivuを示す図である。本実施形態は、こ
のように量子化、あるいは逆量子化によって生じる係数
値誤差ivuを予測し、各量子化係数qvuを算出して量子
化テーブルQy、Qcを作成することにより、量子化誤
差を小さくする構成を備えている。
[0027] FIG. 4 is a diagram illustrating the DCT coefficients F vu, a coefficient value error i vu for each coefficient of the inverse quantized DCT coefficients F 'vu after quantization and inverse quantization. The present embodiment predicts the coefficient value error i vu caused by the quantization or inverse quantization in this way, calculates each quantization coefficient q vu , and creates the quantization tables Qy and Qc to perform the quantization. A configuration is provided to reduce the error.

【0028】図1および図5を参照して、量子化テーブ
ル作成回路20における量子化テーブル作成処理を説明
する。例えば8×8画素のブロック数がMである画像デ
ータ1枚分が入力される(以下、最新画像データとい
う)と、まず図示しない信号処理回路においてブロック
毎に輝度データYと色差データCb、Crに変換され、
DCT処理回路12においてDCTが施されて、DCT
係数に変換され、統計量格納部21に出力される。
Referring to FIGS. 1 and 5, the quantization table creation processing in the quantization table creation circuit 20 will be described. For example, when one piece of image data of 8 × 8 pixels having M blocks is input (hereinafter referred to as the latest image data), first, in a signal processing circuit (not shown), luminance data Y and chrominance data Cb, Cr Is converted to
The DCT is performed in the DCT processing circuit 12, and the DCT is performed.
It is converted into a coefficient and output to the statistic storage unit 21.

【0029】統計量格納部21において、最新画像デー
タのDCT係数に関し、v×u×M個のDCT係数か
ら、予測誤差値の統計量が算出される。統計量格納部2
1は、最新画像データの直前に得られた所定枚数(N
枚)の画像データにおける各周波数の標本数のデータ
(以下標本数データという)を記録しており、新たに標
本数データが入力されると、一番古い標本数データが消
去され、最新の標本数データを含む新たなN枚分の標本
数データに書き換えられる。
In the statistic storage unit 21, the statistic of the prediction error value is calculated from the v × u × M DCT coefficients with respect to the DCT coefficient of the latest image data. Statistics storage 2
1 is a predetermined number (N) obtained immediately before the latest image data.
The number of samples at each frequency in the image data (hereinafter referred to as sample number data) is recorded. When new sample number data is input, the oldest sample number data is deleted and the latest sample number data is deleted. The data is rewritten to new sample number data of N sheets including numerical data.

【0030】量子化誤差予測部22では、統計量格納部
21から得られたN枚分の標本数データから量子化誤差
を予測する。量子化テーブル算出部24では、量子化誤
差予測部22による量子化誤差予測結果に基づいて、量
子化テーブルQya、Qcaが算出され、量子化処理回
路14に出力される。
The quantization error prediction section 22 predicts a quantization error from the sample number data of N pieces obtained from the statistic storage section 21. The quantization table calculator 24 calculates the quantization tables Qya and Qca based on the result of the quantization error prediction by the quantization error predictor 22, and outputs it to the quantization processing circuit 14.

【0031】量子化誤差予測部22における量子化誤差
予測について詳述する。量子化テーブルQyにおける量
子化係数をqvuとすると、DCT係数Fvuと逆量子化D
CT係数F’vuとの係数値誤差ivuの範囲は、次の
(2)式により表される。
The quantization error prediction in the quantization error prediction section 22 will be described in detail. If the quantization coefficient in the quantization table Qy is q vu , the DCT coefficient F vu and the inverse quantization D
The range of the coefficient value error i vu from the CT coefficient F ′ vu is represented by the following equation (2).

【0032】 −qvu/2≦ivu≦qvu/2 ・・・(2)−q vu / 2 ≦ ivu ≦ q vu / 2 (2)

【0033】例えば量子化係数qvu=6であれば、係数
値誤差ivuは−3、−2、−1、0、1、2、3の何れ
かに相当する。係数値誤差ivuの範囲は、量子化係数q
vuが大きくなるにつれ拡大する。
For example, if the quantization coefficient q vu = 6, the coefficient value error i vu corresponds to any of -3, -2, -1, 0, 1, 2, and 3. The range of the coefficient value error i vu is the quantization coefficient q
Expands as vu grows.

【0034】次にM×Nブロックの係数値誤差ivuが、
(2)式に示した範囲内で平均的に分布すると予測し、
(3)式によって係数値誤差ivuの予測値rmsvuを求
める。なお(3)式において、iは係数値誤差ivu、q
は量子化係数qvu、rmsは各空間周波数における係数
値誤差ivuの二乗平均値の平方根、即ちrms値を示
す。例えばq=6であれば、rms=約1.7795で
ある。
Next, the coefficient value error i vu of the M × N block is
It is predicted that the distribution is averaged within the range shown in the equation (2),
The predicted value rms vu of the coefficient value error i vu is obtained by the equation (3). In the equation (3), i is a coefficient value error i vu , q
Indicates a quantization coefficient q vu , and rms indicates a square root of a root mean square value of a coefficient value error i vu at each spatial frequency, that is, an rms value. For example, if q = 6, rms = about 1.7795.

【0035】[0035]

【数1】 (Equation 1)

【0036】図6に、量子化係数qとrms値との対応
表を示す。本実施形態では、輝度データYは8ビット
(256段階)で表され、量子化係数qは1から255
まで変化することとする。この表によく示されるよう
に、量子化係数qが大きくなると、係数値誤差iの範囲
も大きくなるため、rms値は拡大する。
FIG. 6 shows a correspondence table between the quantization coefficient q and the rms value. In the present embodiment, the luminance data Y is represented by 8 bits (256 steps), and the quantization coefficient q is 1 to 255.
To change. As is well shown in this table, as the quantization coefficient q increases, the range of the coefficient value error i also increases, so that the rms value increases.

【0037】本実施形態では、M×N個の全ブロックに
おいて所定の条件を満たす、例えば各空間周波数におけ
るDCT係数Fvuの絶対値が2より大きいDCT係数F
vuの数を、標本数Evuとする(0≦Evu≦M×N)。本
実施形態ではDCT係数Fvuの絶対値が2より大きい場
合、量子化係数qvuの値が大きくなるにつれ係数値誤差
vuが大きくなる、また標本数Evuが多くなるにつれ、
画像全体の量子化誤差が拡大することに注目して、予測
誤差合計値Gvuを次の(4)式で表している。
In the present embodiment, the DCT coefficient F which satisfies a predetermined condition in all M × N blocks, for example, the absolute value of the DCT coefficient F vu at each spatial frequency is larger than 2
The number of vu is set to the number of samples E vu (0 ≦ E vu ≦ M × N). In this embodiment, when the absolute value of the DCT coefficient F vu is greater than 2, the coefficient value error i vu increases as the value of the quantization coefficient q vu increases, and as the number of samples E vu increases,
Paying attention to the fact that the quantization error of the entire image increases, the prediction error total value G vu is expressed by the following equation (4).

【0038】 Gvu=rmsvu×Evu ・・・(4)G vu = rms vu × E vu (4)

【0039】予測誤差合計値Gvuが各空間周波数につい
て均一であれば、画像全体の量子化誤差は減少し、画像
伸張時において復元性の高い再生画像データが得られ
る。従って本実施形態では、空間周波数の低い方、特に
DC成分に画像情報が集中する性質を利用して、まずD
C成分であるDCT係数F00の予測誤差合計値G00が、
予め決定された所定の量子化係数q00から(3)式およ
び(4)式により求められる。他のDCT係数F
vu(v,u≠0)、即ち63個のAC成分の予測誤差合
計値Gvuは、DC成分の予測誤差合計値G00と同じ値と
して、(5)式により対応する予測誤差rms’vuが算
出され、量子化テーブル算出部24に出力される。
If the prediction error total value G vu is uniform for each spatial frequency, the quantization error of the entire image is reduced, and reproduced image data with high resilience at the time of image expansion is obtained. Therefore, in the present embodiment, first, using the lower spatial frequency, in particular, the property that image information concentrates on the DC component,
The total prediction error G00 of the DCT coefficient F00, which is the C component, is
It is obtained from the predetermined quantization coefficient q00 determined in advance by the equations (3) and (4). Other DCT coefficients F
vu (v, u ≠ 0), that is, the prediction error total value G vu of the 63 AC components is the same as the prediction error total value G 00 of the DC component, and the corresponding prediction error rms ′ vu is calculated by the equation (5). It is calculated and output to the quantization table calculation unit 24.

【0040】 rms’vu=G00/Evu ・・・(5) ただし u,v≠0Rms' vu = G 00 / E vu (5) where u, v ≠ 0

【0041】次に、量子化テーブル算出部24におけ
る、量子化テーブル作成処理を説明する。前述の量子化
誤差予測により得られた予測誤差rms’vuは、(6)
式に代入され、(6)式を満たすq’に最も近い整数で
ある予測量子化係数q’vu(v,u≠0)が求められ
る。
Next, the quantization table creation processing in the quantization table calculation unit 24 will be described. The prediction error rms' vu obtained by the above-described quantization error prediction is (6)
The prediction quantization coefficient q ′ vu (v, u ≠ 0), which is an integer closest to q ′ that satisfies the expression (6), is substituted into the expression, and is obtained.

【0042】[0042]

【数2】 (Equation 2)

【0043】例えば、480 ×720 画素の画像を10枚処
理する場合、画像1枚当たりのブロック数Mは5400
であり、総ブロック数M×Nは5400×10=54000である。
このときのDC成分に対応する量子化係数q00を6と
し、DCT係数F00の絶対値が2より大きい標本数E00
を53940 とすると、(3)式および(4)式により予測
誤差合計値G00=1.7795×53940 =95986.23が求められ
る。AC成分、例えばDCT係数F77の絶対値が2より
大きい標本数E77を31520 とすると、(5)式により予
測量子化係数rms’77=95986.23÷31520 =約3.0
452が求められ、(6)式を満たす予測量子化係数
q’77=11が求められる。
For example, when processing 10 images of 480 × 720 pixels, the number M of blocks per image is 5,400.
And the total number of blocks M × N is 5400 × 10 = 54000.
At this time, the quantization coefficient q00 corresponding to the DC component is set to 6, and the sample number E00 in which the absolute value of the DCT coefficient F00 is larger than 2
Is 53940, the total prediction error value G00 = 1.7995 × 53940 = 95986.23 is obtained from the equations (3) and (4). Assuming that the AC component, for example, the sample number E77 in which the absolute value of the DCT coefficient F77 is larger than 2, is 31520, the prediction quantization coefficient rms'77 = 95986.23 ÷ 31520 = approximately 3.0 by equation (5).
452 is obtained, and a predicted quantization coefficient q′77 = 11 that satisfies the expression (6) is obtained.

【0044】以上のように求められた各予測量子化係数
q’vuは、圧縮比に応じてスケーリングされ、8×8の
マトリクスQyaの形で量子化処理回路14に出力され
る。画像圧縮では、量子化係数qを大きくすると高圧縮
ができるが、画質は劣化する。逆に量子化係数qを小さ
くすると画質は向上するが、圧縮後のデータ量が大きく
なる。画質は、量子化係数qに所定の圧縮比をかけるこ
とによりコントロールされる。
Each predicted quantization coefficient q ′ vu obtained as described above is scaled according to the compression ratio, and output to the quantization processing circuit 14 in the form of an 8 × 8 matrix Qya. In image compression, high compression can be achieved by increasing the quantization coefficient q, but image quality deteriorates. Conversely, reducing the quantization coefficient q improves image quality, but increases the amount of data after compression. The image quality is controlled by multiplying the quantization coefficient q by a predetermined compression ratio.

【0045】図7は作成された量子化テーブルの一例を
示す。図7(a)は輝度値Y用の量子化テーブルQya
であり、図7(b)は色差Cb、Cr用の量子化テーブ
ルQcaである。図2と比較すると、特にAC成分にお
いてそれぞれ対応する量子化係数に違いが見られる。
FIG. 7 shows an example of the created quantization table. FIG. 7A shows a quantization table Qya for the luminance value Y.
FIG. 7B shows a quantization table Qca for the color differences Cb and Cr. Compared to FIG. 2, there is a difference particularly in the corresponding quantization coefficient in the AC component.

【0046】図8、図9のフローチャートを参照して、
統計量格納部21における統計処理を説明する。図8お
よび図9において、BLOCKSはブロック数Mを示す
変数であり、BAは処理する総ブロック数(M×N)と
等しい。FvuはDCT係数を示し、添字v、uは、それ
ぞれ0から7まで変化する。Evuは標本数を示す変数で
ある。
Referring to the flowcharts of FIGS. 8 and 9,
The statistical processing in the statistic storage unit 21 will be described. 8 and 9, BLOCKS is a variable indicating the number M of blocks, and BA is equal to the total number of blocks to be processed (M × N). F vu indicates a DCT coefficient, and the subscripts v and u change from 0 to 7, respectively. E vu is a variable indicating the number of samples.

【0047】まず図8を参照する。ステップS102で
は、変数BLOCKSと変数Evuとの初期値が、それぞ
れ0に設定される。ステップS104では変数vの初期
値が0、ステップS106では変数uの初期値が0に設
定される。ステップS108ではDCT係数Fvuの絶対
値が2よりも大きいか否かが判定される。DCT係数F
vuの絶対値が2よりも大きければ、ステップS110に
おいて標本数Evuが1インクリメントされステップS1
12に進む。DCT係数Fvuの絶対値が2よりも大きく
なければ、ステップS112に進み、変数uが1インク
リメントされる。
First, reference is made to FIG. In step S102, the initial values of the variables BLOCKS and E vu are set to 0, respectively. In step S104, the initial value of the variable v is set to 0, and in step S106, the initial value of the variable u is set to 0. In step S108, it is determined whether the absolute value of the DCT coefficient F vu is larger than 2. DCT coefficient F
If the absolute value of vu is larger than 2, the number of samples E vu is incremented by one in step S110, and the value in step S1 is increased.
Proceed to 12. If the absolute value of the DCT coefficient F vu is not larger than 2, the process proceeds to step S112, where the variable u is incremented by one.

【0048】ステップS114では変数uが8であるか
否かが判定され、uが8でない、即ち7以下であれば、
ステップS108から再実行される。u=8であればス
テップS116に進み、変数vが1インクリメントされ
る。同様に、ステップS118では変数vが8であるか
否かが判定され、vが8でない、即ち7以下であれば、
ステップS106から再実行される。v=8であればス
テップS120に進む。
In step S114, it is determined whether or not the variable u is 8, and if u is not 8, ie, 7 or less,
The processing is executed again from step S108. If u = 8, the process proceeds to step S116, where the variable v is incremented by one. Similarly, in step S118, it is determined whether or not the variable v is 8, and if v is not 8, ie, 7 or less,
The processing is executed again from step S106. If v = 8, the process proceeds to step S120.

【0049】ステップS120では、変数BLOCKS
が1インクリメントされ、ステップS122において変
数BLOCKSが総ブロック数BAであるか否かが判定
される。BLOCKSがBAでなければ、ステップS1
04から再実行され、BLOCKS=BAであればステ
ップS124に進む。
In step S120, the variable BLOCKS
Is incremented by one, and it is determined in step S122 whether the variable BLOCKS is the total block number BA. If BLOCKS is not BA, step S1
04, and if BLOCKS = BA, the process proceeds to step S124.

【0050】このように、8×8のDCT係数のマトリ
クスFvuにおいて、まず一番上の横一列の標本数Evu
左から順にカウントされ、順に次の横一列の標本数Evu
がカウントされる。即ち図8の処理では、64個の各周
波数成分において、DCT係数Fvuの絶対値が2よりも
大きい標本数Evuがそれぞれカウントされる。
As described above, in the matrix F vu of the DCT coefficients of 8 × 8, the sample number E vu of the top horizontal row is counted sequentially from the left, and the sample number E vu of the next horizontal row is sequentially counted.
Is counted. That is, in the process of FIG. 8, the number of samples E vu whose absolute value of the DCT coefficient F vu is larger than 2 is counted in each of the 64 frequency components.

【0051】次に図9を参照する。ステップS124で
は変数vの初期値が0、ステップS126では変数uの
初期値が0、ステップS128では画像数を表す変数i
の初期値が0に設定される。ステップS130では、E
vu[i+1] はEvu[i] に変換される。例えばE[1] がE
[0] になり、E[6] はE[5] になる。ステップS134
では変数iが画像数より1枚少ない数(N−1)と等し
いか否かが判定される。変数iが(N−1)と等しくな
ければステップS130が再実行され、変数iが(N−
1)と等しければステップS136に進む。ステップS
136では、Evu[N-1] がEvuに設定され、ステップS
138に進む。
Next, reference is made to FIG. In step S124, the initial value of the variable v is 0, in step S126, the initial value of the variable u is 0, and in step S128, the variable i representing the number of images
Is set to 0. In step S130, E
vu [i + 1] is converted to E vu [i]. For example, if E [1] is E
[0] and E [6] becomes E [5]. Step S134
In, it is determined whether or not the variable i is equal to the number (N-1) one less than the number of images. If the variable i is not equal to (N-1), step S130 is executed again, and the variable i becomes (N-
If it is equal to 1), the process proceeds to step S136. Step S
At 136, E vu [N-1] is set to E vu , and step S
Proceed to 138.

【0052】ステップS138では、変数uが1インク
リメントされ、ステップS140では変数uが8である
か否かが判定される。uが8でない、即ち7以下であれ
ば、ステップS128から再実行される。u=8であれ
ばステップS142に進み、変数vが1インクリメント
される。同様に、ステップS144では変数vが8であ
るか否かが判定され、vが8でない、即ち7以下であれ
ば、ステップS126から再実行される。v=8であれ
ば処理は終了する。
In step S138, the variable u is incremented by one, and in step S140, it is determined whether or not the variable u is 8. If u is not 8, that is, 7 or less, the process is repeated from step S128. If u = 8, the process proceeds to step S142, and the variable v is incremented by one. Similarly, in step S144, it is determined whether or not the variable v is 8, and if v is not 8, that is, 7 or less, the process is repeated from step S126. If v = 8, the process ends.

【0053】例えば、N=10枚分の標本数データ(E
[0] 、E[1] 、...E[9] ;括弧内の数値が大きいほ
ど最新の画像データであることを示す)が、統計量格納
部21に格納されているとする。最新標本数であるEが
統計量格納部21に入力されると、格納されていた標本
数データE[i] の添字i(1≦i≦N−1)がそれぞれ
1繰り下がり(ステップS130)、標本数データはE
[0] 、E[1] 、...E[8] となる。そして最新標本数
データEが、E[9] に更新される(ステップS13
6)。更新されたE[0] 〜E[8] 、およびE[9] は、再
び統計量格納部21に格納されるとともに、量子化誤差
予測部22に出力される。このように図9の処理では、
統計量格納部21に格納されているN枚分の標本数デー
タEを、最新の撮影画像のN枚の標本数データEで常に
更新している。
For example, N = 10 sample number data (E
[0], E [1],. . . E [9]; the larger the numerical value in parentheses, the more recent the image data is.) Is stored in the statistic storage unit 21. When the latest sample number E is input to the statistic storage unit 21, the subscript i (1 ≦ i ≦ N-1) of the stored sample number data E [i] is lowered by one (step S130). , The sample size data is E
[0], E [1],. . . E [8]. Then, the latest sample number data E is updated to E [9] (step S13).
6). The updated E [0] to E [8] and E [9] are stored again in the statistic storage unit 21 and output to the quantization error prediction unit 22. Thus, in the processing of FIG.
The N sample number data E stored in the statistic storage unit 21 is constantly updated with the N sample number data E of the latest captured image.

【0054】図10、図11のフローチャートを参照し
て、量子化誤差予測部22における量子化誤差予測処理
を説明する。図10および図11において、Evuは標本
数を示す変数である。
The quantization error prediction processing in the quantization error prediction section 22 will be described with reference to the flowcharts of FIGS. 10 and 11, E vu is a variable indicating the number of samples.

【0055】まず図10を参照する。ステップS202
では、変数Evuの初期値が0に設定される。ステップS
204では変数vの初期値が0、ステップS206では
変数uの初期値が0に設定される。ステップS208で
は変数iの初期値が0に設定される。
First, reference is made to FIG. Step S202
Then, the initial value of the variable E vu is set to 0. Step S
At 204, the initial value of the variable v is set to 0, and at step S206, the initial value of the variable u is set to 0. In step S208, the initial value of the variable i is set to 0.

【0056】ステップS210では、各空間周波数にお
ける標本数Evu[i] の加算が行なわれ、ステップS21
2において変数iが1インクリメントされ、ステップS
214に進む。ステップS214では変数iがNである
か否かが判定され、iがNでなければステップS210
から再実行される。変数iがNであると判定されるとス
テップS216に進む。即ち、ステップS208からス
テップS214までの処理によって、各空間周波数にお
ける標本数EvuのN個の総和が求められる。
In step S210, the number of samples E vu [i] at each spatial frequency is added, and step S21 is performed.
In step 2, the variable i is incremented by one,
Proceed to 214. In step S214, it is determined whether or not the variable i is N. If i is not N, step S210
Is re-executed from If it is determined that the variable i is N, the process proceeds to step S216. That is, through the processing from step S208 to step S214, N total sums of the number of samples E vu at each spatial frequency are obtained.

【0057】ステップS216において、変数uが1イ
ンクリメントされ、ステップS218では変数uが8で
あるか否かが判定され、uが8でない、即ち7以下であ
れば、ステップS208から再実行される。u=8であ
ればステップS220に進み、変数vが1インクリメン
トされる。同様に、ステップS222では変数vが8で
あるか否かが判定され、vが8でない、即ち7以下であ
れば、ステップS206から再実行される。v=8であ
ればステップS302に進む。
In step S216, the variable u is incremented by one. In step S218, it is determined whether or not the variable u is 8, and if u is not 8, ie, 7 or less, the process is re-executed from step S208. If u = 8, the process proceeds to step S220, and the variable v is incremented by one. Similarly, in step S222, it is determined whether or not the variable v is 8, and if v is not 8, ie, 7 or less, the process is re-executed from step S206. If v = 8, the process proceeds to step S302.

【0058】このように図10の処理では、図8のフロ
ーチャートにおいてカウントされた各画像iの標本数E
vu[i] の総和Evuを各空間周波数毎に求めている。
As described above, in the processing of FIG. 10, the number of samples E of each image i counted in the flowchart of FIG.
The sum Evu of vu [i] is obtained for each spatial frequency.

【0059】次に図11を参照する。ステップS302
では変数vの初期値が0、ステップS304では変数u
の初期値が0に設定される。ステップS306では変数
v、uが共に0であるか否かが判定される。変数v、u
が共に0であればステップS308に進み、変数v、u
が共に0でなければ、ステップS310に進む。
Next, reference is made to FIG. Step S302
Then, the initial value of the variable v is 0, and in step S304, the variable u
Is set to 0. In step S306, it is determined whether both the variables v and u are 0. Variables v, u
If both are 0, the process proceeds to step S308, where the variables v and u
Are not 0, the process proceeds to step S310.

【0060】ステップS308では、所定の量子化係数
q00に基づいて、(3)式を用いてq00のrms値、即
ちrms00が求められる。このrms00と標本数E00と
の乗算により、予測誤差合計値G00が求められ、ステッ
プS312に進む。
In step S308, the rms value of q00, that is, rms00, is obtained using equation (3) based on the predetermined quantization coefficient q00. By multiplying this rms00 by the number of samples E00, the total prediction error G00 is obtained, and the process proceeds to step S312.

【0061】ステップS310では、予測誤差合計値G
00と標本数Evuを(5)式に代入して、rms値である
rms’vuが求められる。その後、rms’vuを(6)
式に代入することにより予測量子化係数q’vuが求めら
れ、ステップS312に進む。
In step S310, the prediction error total value G
By substituting 00 and the number of samples E vu into the equation (5), the rms value rms ′ vu is obtained. Then, rms' vu (6)
The predicted quantization coefficient q ′ vu is obtained by substituting into the equation, and the process proceeds to step S312.

【0062】ステップS312では、変数uが1インク
リメントされ、ステップS314では変数uが8である
か否かが判定される。uが8でない、即ち7以下であれ
ば、ステップS306から再実行される。u=8であれ
ばステップS316に進み、変数vが1インクリメント
される。同様に、ステップS318では変数vが8であ
るか否かが判定され、vが8でない、即ち7以下であれ
ば、ステップS304から再実行される。v=8であれ
ば処理は終了する。
In step S312, the variable u is incremented by one, and in step S314, it is determined whether or not the variable u is 8. If u is not 8, that is, 7 or less, the process is re-executed from step S306. If u = 8, the process proceeds to step S316, where the variable v is incremented by one. Similarly, in step S318, it is determined whether or not the variable v is 8, and if v is not 8, ie, 7 or less, the process is re-executed from step S304. If v = 8, the process ends.

【0063】このように、まず量子化係数q00によりD
C成分の予測誤差合計値G00を求め、DC成分のrms
00を決定する(ステップS308)。そして、AC成分
の各rms値rms’vuを、各総和標本数Evuと予測誤
差合計値G00とから求め、予測量子化係数q’vu(v,
u≠0)を算出する(ステップS310)。即ち図11
の処理では、DC成分の予測誤差合計値G00を基に、A
C成分の各周波数成分に対応した63個の予測量子化係
数q’vuが決定される。
As described above, first, D is calculated by the quantization coefficient q00.
The total prediction error G00 of the C component is obtained, and the rms of the DC component is obtained.
00 is determined (step S308). Then, each rms value rms ' vu of the AC component is obtained from each total sample number E vu and the prediction error total value G00, and the prediction quantization coefficient q' vu (v,
u ≠ 0) is calculated (step S310). That is, FIG.
In the processing of A, based on the total prediction error G00 of the DC component, A
63 predicted quantization coefficients q ′ vu corresponding to each frequency component of the C component are determined.

【0064】第1実施形態には量子化テーブル作成回路
20が設けられ、この量子化テーブル作成回路20は、
原画像データのDCT係数から、量子化誤差に関係の深
い統計量である標本数データを各空間周波数毎に算出
し、所定枚数標本数データを格納し、画像入力毎に更新
する統計量格納部21を備えている。また、この標本数
データを用いて量子化誤差を予測する量子化誤差予測部
22と、量子化係数を量子化誤差に基づいて求め、量子
化テーブルを出力する量子化テーブル算出部24とを備
えている。従って、所定枚数の最新の画像に最適な単一
の量子化テーブルQyaを作成できる。従って同じ圧縮
画像データ量で量子化誤差の少ない再生画像データを得
ることができる。
In the first embodiment, a quantization table creation circuit 20 is provided.
A statistic storage unit that calculates, for each spatial frequency, sample number data, which is a statistic closely related to the quantization error, from DCT coefficients of the original image data, stores a predetermined number of sample number data, and updates each time an image is input 21. The image processing apparatus further includes a quantization error prediction unit 22 that predicts a quantization error using the sample number data, and a quantization table calculation unit 24 that obtains a quantization coefficient based on the quantization error and outputs a quantization table. ing. Accordingly, a single quantization table Qya optimal for a predetermined number of latest images can be created. Therefore, reproduced image data with a small quantization error can be obtained with the same amount of compressed image data.

【0065】図12は本発明の第2実施形態を示す図で
ある。第1実施形態と同一の構成には第1実施形態と同
符号を付しており、ここでは説明を省略する。第2実施
形態では、量子化テーブル作成回路20によって作成さ
れた量子化テーブルをQ’y(輝度データ用)、Q’c
(色差データ用)とする。
FIG. 12 is a view showing a second embodiment of the present invention. The same components as those of the first embodiment are denoted by the same reference numerals as those of the first embodiment, and description thereof is omitted here. In the second embodiment, the quantization tables created by the quantization table creation circuit 20 are Q′y (for luminance data), Q′c
(For color difference data).

【0066】量子化テーブル格納部30は、例えばデフ
ォオルトの量子化テーブルQy、Qcが書き込まれたメ
モリである。DCT処理回路12からにより得られるD
CT係数Fvuは、量子化処理回路14に入力されると共
に、量子化テーブル作成回路20にも入力される。量子
化処理回路14ではデフォルトの量子化テーブルQy、
Qcを用いて量子化が行なわれる。量子化された画像デ
ータは符号化され、圧縮画像データとして記録媒体Mの
圧縮画像データ記録領域M1に記録される。デフォルト
の量子化テーブルQy、Qcは、記録媒体Mのテーブル
記録領域M2に記録される。
The quantization table storage section 30 is a memory in which, for example, default quantization tables Qy and Qc are written. D obtained from the DCT processing circuit 12
The CT coefficient F vu is input to the quantization processing circuit 14 and also to the quantization table creation circuit 20. In the quantization processing circuit 14, default quantization tables Qy,
The quantization is performed using Qc. The quantized image data is encoded and recorded in the compressed image data recording area M1 of the recording medium M as compressed image data. The default quantization tables Qy and Qc are recorded in the table recording area M2 of the recording medium M.

【0067】次に量子化テーブル作成回路20では、新
たに入力されたDCT係数による量子化誤差と、統計量
格納部21に記録された過去の量子化誤差とに基づい
て、量子化テーブルQ’y、Q’cが作成され、量子化
テーブル格納部30に出力される。量子化テーブル格納
部30は、量子化テーブル作成回路から得られた量子化
テーブルQ’y、Q’cを、デフォルトの量子化テーブ
ルQy、Qcとして書き換える。
Next, the quantization table creation circuit 20 calculates the quantization table Q ′ based on the quantization error due to the newly input DCT coefficient and the past quantization error recorded in the statistic storage unit 21. y and Q′c are created and output to the quantization table storage unit 30. The quantization table storage unit 30 rewrites the quantization tables Q′y and Q′c obtained from the quantization table creation circuit as default quantization tables Qy and Qc.

【0068】従って、量子化に用いられる量子化テーブ
ルQy、Qcは、以前に入力された画像データに基づい
て作成されたものである。第2実施形態においても、第
1実施形態と同様、各画像のDCT係数から量子化誤差
を予測し、誤差の統計量を更新することにより量子化テ
ーブルQy、Qcを作成する構成であり、同じ圧縮画像
データ量で量子化誤差の少ない再生画像データを得るこ
とができる。さらに第2実施形態では、量子化テーブル
Qy、Qcは量子化テーブル格納部30から読み込むだ
けで良いので、量子化テーブルの算出結果を待たずに量
子化でき、処理時間の短縮も図られる。
Therefore, the quantization tables Qy and Qc used for the quantization are created based on the image data previously input. In the second embodiment, as in the first embodiment, the quantization errors are predicted from the DCT coefficients of the respective images, and the quantization tables Qy and Qc are created by updating the statistics of the errors. Reproduced image data with a small quantization error can be obtained with the amount of compressed image data. Further, in the second embodiment, since the quantization tables Qy and Qc need only be read from the quantization table storage unit 30, quantization can be performed without waiting for the calculation result of the quantization table, and the processing time can be reduced.

【0069】[0069]

【発明の効果】本発明によると、量子化誤差を予測し、
量子化と逆量子化による画像劣化が少ない画像圧縮を容
易にする画像圧縮装置、およびカメラを提供することが
できる。
According to the present invention, the quantization error is predicted,
It is possible to provide an image compression device and a camera that facilitate image compression with little image degradation due to quantization and inverse quantization.

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

【図1】本発明によるカメラの第1実施形態を示すブロ
ック図である。
FIG. 1 is a block diagram showing a first embodiment of a camera according to the present invention.

【図2】JPEG推奨の輝度データ用量子化テーブルと
色差データ用量子化テーブルを示す図である。
FIG. 2 is a diagram showing a quantization table for luminance data and a quantization table for color difference data recommended by JPEG.

【図3】輝度データYのDCT係数Fvuと、量子化DC
T係数Rvuと、逆量子化DCT係数F’vu、および量子
化テーブルQyを示す図である。
FIG. 3 shows a DCT coefficient F vu of luminance data Y and a quantized DC
It is a figure which shows the T coefficient Rvu , the inverse quantization DCT coefficient F'vu , and the quantization table Qy.

【図4】DCT係数Fvuと、逆量子化DCT係数F’vu
との各係数毎の係数値誤差ivuを示す図である。
FIG. 4 shows a DCT coefficient F vu and an inversely quantized DCT coefficient F ′ vu
FIG. 9 is a diagram showing a coefficient value error i vu for each coefficient.

【図5】量子化テーブル作成回路におけるフローを示す
ブロック図である。
FIG. 5 is a block diagram illustrating a flow in a quantization table creation circuit.

【図6】量子化係数qとrms値との対応を示す表であ
る。
FIG. 6 is a table showing a correspondence between a quantization coefficient q and an rms value.

【図7】作成された輝度データ用量子化テーブルQya
と色差データ用量子化テーブルQcaとを示す図であ
る。
FIG. 7 shows a created quantization table for luminance data Qya.
FIG. 7 is a diagram showing a color difference data quantization table Qca.

【図8】量子化誤差統計処理のフローチャートの前半部
分を示す図である。
FIG. 8 is a diagram illustrating a first half of a flowchart of a quantization error statistical process;

【図9】量子化誤差統計処理のフローチャートの後半部
分を示す図である。
FIG. 9 is a diagram illustrating the latter half of the flowchart of the quantization error statistical process.

【図10】量子化誤差予測処理のフローチャートの前半
部分を示す図である。
FIG. 10 is a diagram illustrating a first half of a flowchart of a quantization error prediction process.

【図11】量子化誤差予測処理のフローチャートの後半
部分を示す図である。
FIG. 11 is a diagram illustrating the latter half of the flowchart of the quantization error prediction process.

【図12】本発明によるカメラの第2実施形態を示すブ
ロック図である。
FIG. 12 is a block diagram showing a second embodiment of the camera according to the present invention.

【符号の説明】[Explanation of symbols]

10 圧縮画像装置 12 DCT処理回路 14 量子化処理回路 20 量子化テーブル作成回路 21 統計量格納部 22 量子化誤差予測部 24 量子化テーブル算出部 DESCRIPTION OF SYMBOLS 10 Compressed image apparatus 12 DCT processing circuit 14 Quantization processing circuit 20 Quantization table preparation circuit 21 Statistic amount storage part 22 Quantization error prediction part 24 Quantization table calculation part

─────────────────────────────────────────────────────
────────────────────────────────────────────────── ───

【手続補正書】[Procedure amendment]

【提出日】平成9年11月19日[Submission date] November 19, 1997

【手続補正1】[Procedure amendment 1]

【補正対象書類名】明細書[Document name to be amended] Statement

【補正対象項目名】0022[Correction target item name] 0022

【補正方法】変更[Correction method] Change

【補正内容】[Correction contents]

【0022】量子化テーブルQyを用いてDCT係数F
vuを量子化する式は(1)式により定義される。この式
におけるround は、最も近い整数への近似を意味する。
すなわち、DCT係数Fvuおよび量子化テーブルQy
各要素同士の割算と、四捨五入とによって、量子化DC
T係数Rvuが求められる。
Using the quantization table Qy , the DCT coefficient F
The equation for quantizing vu is defined by equation (1). Round in this equation means approximation to the nearest integer.
In other words, the DCT coefficient F vu and the division of each element of the quantization table Qy , and the rounding-off, make the quantization DC
The T coefficient R vu is determined.

【手続補正2】[Procedure amendment 2]

【補正対象書類名】明細書[Document name to be amended] Statement

【補正対象項目名】0023[Correction target item name] 0023

【補正方法】変更[Correction method] Change

【補正内容】[Correction contents]

【0023】 Rvu=round (Fvuvu ) ・・・(1) {0≦u,v ≦7}R vu = round (F vu / q vu ) (1) {0 ≦ u, v ≦ 7}

【手続補正3】[Procedure amendment 3]

【補正対象書類名】明細書[Document name to be amended] Statement

【補正対象項目名】0025[Correction target item name] 0025

【補正方法】変更[Correction method] Change

【補正内容】[Correction contents]

【0025】符号化された圧縮画像信号を伸張して、画
面に表示するためには、復号化、逆量子化、2次元DC
Tの逆変換(以下、2次元IDCTという)の処理が必
要である。この復号化は、ハフマン符号化とは逆の作用
であり、従来公知であるため詳述しない。復号化によっ
て得られた量子化DCT係数は、量子化に用いた量子化
テーブルQy、Qcを用いてそれぞれ逆量子化され、逆
量子化DCT係数F’vuに変換される。これらの逆量子
化DCT係数F’vu2次元DCTの逆変換である2次
元IDCTを施され、それぞれ輝度データY’、色差デ
ータCb’、Cr’に変換される。2次元IDCTにつ
いても公知であるのでここでは詳述しない。
In order to decompress the coded compressed image signal and display it on the screen, decoding, inverse quantization, two-dimensional DC
Processing of the inverse transformation of T (hereinafter referred to as two-dimensional IDCT) is required. This decoding is the reverse of the operation of Huffman coding, and is well known in the art and will not be described in detail. The quantized DCT coefficients obtained by the decoding are inversely quantized using the quantization tables Qy and Qc used for quantization, and are converted into inversely quantized DCT coefficients F ′ vu . These inversely quantized DCT coefficients F ′ vu are second-order inverse transforms of two-dimensional DCT.
The original IDCT is performed, and the data is converted into luminance data Y ′ and chrominance data Cb ′ and Cr ′, respectively. Since the two-dimensional IDCT is also known, it will not be described in detail here.

【手続補正4】[Procedure amendment 4]

【補正対象書類名】明細書[Document name to be amended] Statement

【補正対象項目名】0042[Correction target item name] 0042

【補正方法】変更[Correction method] Change

【補正内容】[Correction contents]

【0042】[0042]

【数2】 (Equation 2)

【手続補正5】[Procedure amendment 5]

【補正対象書類名】明細書[Document name to be amended] Statement

【補正対象項目名】符号の説明[Correction target item name] Explanation of sign

【補正方法】変更[Correction method] Change

【補正内容】[Correction contents]

【符号の説明】 10 画像圧縮装置 12 DCT処理回路 14 量子化処理回路 20 量子化テーブル作成回路 21 統計量格納部 22 量子化誤差予測部 24 量子化テーブル算出部[Description of Signs] 10 Image compression apparatus 12 DCT processing circuit 14 Quantization processing circuit 20 Quantization table creation circuit 21 Statistic amount storage unit 22 Quantization error prediction unit 24 Quantization table calculation unit

Claims (9)

【特許請求の範囲】[Claims] 【請求項1】 撮影光学系から得られた原画像データに
対応した直交変換係数を、量子化テーブルを用いて量子
化し、量子化直交変換係数を求める量子化手段と、 前記直交変換係数に基づいて、前記量子化テーブルを更
新する量子化テーブル作成手段とを備えたことを特徴と
する画像圧縮装置。
1. A quantization means for quantizing an orthogonal transform coefficient corresponding to original image data obtained from an imaging optical system using a quantization table to obtain a quantized orthogonal transform coefficient, based on the orthogonal transform coefficient. And a quantization table creating means for updating the quantization table.
【請求項2】 前記量子化テーブル作成手段が、 前記量子化直交変換係数を逆量子化して求められる逆量
子化直交変換係数と前記直交変換係数との誤差の絶対値
が2よりも大きい前記直交変換係数の標本数を、前記直
交変換係数毎にカウントし、各画像における前記標本数
のデータ(標本数データ)を、所定画像枚数分だけ格納
する格納部を備えたことを特徴とする請求項1に記載の
画像圧縮装置。
2. The method according to claim 1, wherein said quantization table creating means is configured to: inversely quantize said quantized orthogonal transform coefficient and obtain an orthogonal quantization coefficient whose absolute value is greater than 2; A storage unit for counting the number of samples of transform coefficients for each of the orthogonal transform coefficients, and storing data of the number of samples (sample number data) in each image for a predetermined number of images. 2. The image compression device according to 1.
【請求項3】 前記格納部において、新たな画像の前記
標本数データが入力される毎に、前記格納部に格納され
た最も古い画像の前記標本数データが消去されることに
より、前記所定枚数の前記標本数データが更新されるこ
とを特徴とする請求項2に記載の画像圧縮装置。
3. The storage unit deletes the sample number data of the oldest image stored in the storage unit every time the sample number data of a new image is input, thereby obtaining the predetermined number of samples. The image compression apparatus according to claim 2, wherein the sample number data is updated.
【請求項4】 前記量子化テーブル作成手段が、 所定枚数の前記標本数データに基づいて、前記誤差を予
測する量子化誤差予測部と、 前記量子化誤差予測部から得られた予測誤差値に基づい
て量子化テーブルを算出し、前記量子化手段に出力する
量子化テーブル算出部とを備えることを特徴とする請求
項3に記載の画像圧縮装置。
4. The quantization table creation means includes: a quantization error prediction unit that predicts the error based on a predetermined number of sample number data; and a prediction error value obtained from the quantization error prediction unit. The image compression apparatus according to claim 3, further comprising: a quantization table calculation unit that calculates a quantization table based on the quantization table and outputs the quantization table to the quantization unit.
【請求項5】 前記量子化テーブル作成手段が、新規に
入力された前記標本数データを含めた前記所定枚数の前
記標本数データにより、前記量子化テーブルを作成する
ことを特徴とする請求項4に記載の画像圧縮装置。
5. The quantization table creating means creates the quantization table from the predetermined number of sample number data including the newly input sample number data. An image compression device according to claim 1.
【請求項6】 前記量子化テーブル作成手段が、新規に
入力された前記標本数データより前に入力された前記所
定枚数の前記標本数データにより、前記量子化テーブル
を作成することを特徴とする請求項4に記載の画像圧縮
装置。
6. The quantization table creating means creates the quantization table based on the predetermined number of sample number data input before the newly input sample number data. The image compression device according to claim 4.
【請求項7】 撮影光学系から得られた原画像データに
対応した直交変換係数を、量子化テーブルを用いて量子
化し、量子化直交変換係数を求める量子化手段と、 前記直交変換係数に基づいて、前記量子化テーブルを更
新する量子化テーブル作成手段とを備え、 前記量子化テーブル作成手段が、1画像分の画像データ
に対応し、前記量子化直交変換係数を逆量子化して求め
られる逆量子化直交変換係数と前記直交変換係数との誤
差の絶対値が2よりも大きい前記直交変換係数の標本数
を、前記直交変換係数毎にカウントし、各画像における
前記標本数のデータ(標本数データ)を、所定画像枚数
分だけ格納可能な標本数データ格納手段と、 前記標本数データ格納手段に格納されている所定枚数の
前記標本数データに基づいて、前記誤差を予測する量子
化誤差予測部と、 前記量子化誤差予測部から得られた予測誤差値に基づい
て量子化テーブルを算出し、前記量子化手段に出力する
量子化テーブル算出部とを備え、 さらに前記量子化テーブル作成手段によって作成された
前記量子化テーブルを用いて画像圧縮された圧縮画像デ
ータを記録媒体に記録する記録手段と、 作成された前記量子化テーブルを格納する量子化テーブ
ル格納手段とを備えることを特徴とするカメラ。
7. A quantization means for quantizing an orthogonal transform coefficient corresponding to original image data obtained from an imaging optical system by using a quantization table to obtain a quantized orthogonal transform coefficient, based on the orthogonal transform coefficient. And a quantization table creating means for updating the quantization table, wherein the quantization table creating means corresponds to image data of one image, and is obtained by inversely quantizing the quantized orthogonal transform coefficients. The number of samples of the orthogonal transform coefficient whose absolute value of the error between the quantized orthogonal transform coefficient and the orthogonal transform coefficient is larger than 2 is counted for each of the orthogonal transform coefficients, and data of the sample number in each image (sample number) Data) for storing a predetermined number of images for a predetermined number of images, and predicting the error based on a predetermined number of the sample number data stored in the sample number data storage means. And a quantization table calculation unit that calculates a quantization table based on the prediction error value obtained from the quantization error prediction unit and outputs the quantization table to the quantization unit. Recording means for recording, on a recording medium, compressed image data image-compressed using the quantization table created by the quantization table creation means, and quantization table storage means for storing the created quantization table. A camera characterized in that:
【請求項8】 前記記録手段によって、圧縮画像データ
が前記記録媒体の第1の領域に記録され、量子化に用い
られた量子化テーブルが前記記録媒体の第2の領域に記
録されることを特徴とする請求項7に記載のカメラ。
8. The recording device according to claim 1, wherein the compressed image data is recorded in a first area of the recording medium, and a quantization table used for quantization is recorded in a second area of the recording medium. The camera according to claim 7, characterized in that:
【請求項9】 前記量子化テーブル作成手段が、 前記テーブル格納手段に格納されている前記量子化テー
ブルを、量子化に用いるとともに、 前記記録媒体への前記圧縮画像データの記録が完了した
後に、この圧縮画像データと対応する前記標本数データ
に基づいた前記量子化テーブルを作成し、 作成された前記量子化テーブルを前記量子化テーブル格
納手段に格納することにより、前記量子化テーブルを更
新することを特徴とする請求項8に記載のカメラ。
9. The quantization table creation unit uses the quantization table stored in the table storage unit for quantization, and after the recording of the compressed image data on the recording medium is completed, Updating the quantization table by creating the quantization table based on the sample number data corresponding to the compressed image data and storing the created quantization table in the quantization table storage means. 9. The camera according to claim 8, wherein:
JP17291897A 1997-06-13 1997-06-13 Image compression apparatus and camera Expired - Fee Related JP3645690B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP17291897A JP3645690B2 (en) 1997-06-13 1997-06-13 Image compression apparatus and camera

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP17291897A JP3645690B2 (en) 1997-06-13 1997-06-13 Image compression apparatus and camera

Publications (2)

Publication Number Publication Date
JPH114438A true JPH114438A (en) 1999-01-06
JP3645690B2 JP3645690B2 (en) 2005-05-11

Family

ID=15950773

Family Applications (1)

Application Number Title Priority Date Filing Date
JP17291897A Expired - Fee Related JP3645690B2 (en) 1997-06-13 1997-06-13 Image compression apparatus and camera

Country Status (1)

Country Link
JP (1) JP3645690B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4899053A (en) * 1987-10-21 1990-02-06 Criticare Systems, Inc. Solid state non-dispersive IR analyzer using electrical current-modulated microsources

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4899053A (en) * 1987-10-21 1990-02-06 Criticare Systems, Inc. Solid state non-dispersive IR analyzer using electrical current-modulated microsources

Also Published As

Publication number Publication date
JP3645690B2 (en) 2005-05-11

Similar Documents

Publication Publication Date Title
US8170355B2 (en) Image encoding/decoding method and apparatus
US20120099787A1 (en) Image encoding/decoding method and apparatus
US20020196983A1 (en) Image encoding apparatus, image encoding method, and recording medium in which image encoding program is recorded
WO2003105070A1 (en) Spatial prediction based intra coding
US20140198990A1 (en) Method and apparatus for encoding and decoding image
JP3579544B2 (en) Image compression / decompression device
US20080193028A1 (en) Method of high quality digital image compression
JP2000244935A (en) Method for compressing picture data
US6396955B1 (en) Image compression and expansion device
JPH1075448A (en) Image compressor and image expander
JP3222780B2 (en) Image compression device
JPH07327232A (en) Picture compression device
JP3375539B2 (en) Image compression device and image decompression device
JP3645690B2 (en) Image compression apparatus and camera
US6445870B1 (en) Image compression and expansion device
US20020173952A1 (en) Coding
JP3645689B2 (en) Image compression apparatus and quantization table creation apparatus
HUE033524T2 (en) Method for compression of data
US6226445B1 (en) Image compression and expansion device
EP1345449A1 (en) Digital image compression and decompression using block coding
JP3559314B2 (en) Image compression device
JP3709106B2 (en) Image compression and decompression device
US6208754B1 (en) Image compression and expansion device using pixel offset
JPH06113291A (en) Picture coder and decoder
JP3655734B2 (en) Image compression apparatus and image display apparatus

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20040224

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20040426

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20040927

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20041007

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: 20050121

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20050204

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: 20090210

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100210

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110210

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120210

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120210

Year of fee payment: 7

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120210

Year of fee payment: 7

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120210

Year of fee payment: 7

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120210

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130210

Year of fee payment: 8

LAPS Cancellation because of no payment of annual fees