JPH0723229A - Picture data compression processing method - Google Patents
Picture data compression processing methodInfo
- Publication number
- JPH0723229A JPH0723229A JP14755493A JP14755493A JPH0723229A JP H0723229 A JPH0723229 A JP H0723229A JP 14755493 A JP14755493 A JP 14755493A JP 14755493 A JP14755493 A JP 14755493A JP H0723229 A JPH0723229 A JP H0723229A
- Authority
- JP
- Japan
- Prior art keywords
- image data
- wavelet transform
- transform coefficient
- signal
- scanning direction
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Compression Or Coding Systems Of Tv Signals (AREA)
- Compression Of Band Width Or Redundancy In Fax (AREA)
- Compression, Expansion, Code Conversion, And Decoders (AREA)
- Image Processing (AREA)
Abstract
Description
【0001】[0001]
【産業上の利用分野】本発明は画像データの圧縮処理方
法、特に詳細にはウェーブレット変換を用いて高いデー
タ圧縮率が得られるようにした画像データの圧縮処理方
法に関するものである。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an image data compression processing method, and more particularly, to an image data compression processing method capable of obtaining a high data compression rate by using wavelet transform.
【0002】[0002]
【従来の技術】例えばTV信号等、中間調画像を担持す
る画像信号は膨大な情報量を有しているので、その伝送
には広帯域の伝送路が必要である。そこで従来より、こ
のような画像信号は冗長性が大きいことに着目し、この
冗長性を抑圧することによって画像データを圧縮する試
みが種々なされている。また最近では、例えば光ディス
クや磁気ディスク等に中間調画像を記録することが広く
行われており、この場合には記録媒体に効率良く画像信
号を記録することを目的として画像データ圧縮が広く適
用されている。2. Description of the Related Art Since an image signal carrying a halftone image such as a TV signal has an enormous amount of information, a wide band transmission line is required for its transmission. Therefore, conventionally, attention has been paid to the fact that such an image signal has large redundancy, and various attempts have been made to compress the image data by suppressing this redundancy. Further, recently, for example, recording of a halftone image on an optical disk or a magnetic disk has been widely performed, and in this case, image data compression is widely applied for the purpose of efficiently recording an image signal on a recording medium. ing.
【0003】このような画像データの圧縮方法の一つと
して、従来から、画像データを格納,伝送等する場合
に、該画像データに予測符号化による圧縮処理を施して
データ量を圧縮減少せしめた上で格納,伝送等を行い、
画像再生の際はその圧縮された画像データ(圧縮画像デ
ータ)に復号処理を施して伸長し、その伸長された画像
データ(伸長画像データ)に基づいて可視像を再生する
ような方法が採用されている。As one of such image data compression methods, conventionally, when image data is stored or transmitted, the image data is compressed by predictive coding to reduce the amount of data. Stored, transmitted, etc. above,
At the time of image reproduction, a method is adopted in which the compressed image data (compressed image data) is subjected to decoding processing and decompressed, and a visible image is reproduced based on the decompressed image data (decompressed image data). Has been done.
【0004】また、画像データ圧縮方法の一つとして、
ベクトル量子化を利用する方法が知られている。この方
法は、2次元画像データを標本数K個のブロックに分割
し、予めK個のベクトル要素を規定して作成した相異な
る複数のベクトルからなるコードブックの中で、上記ブ
ロックの各々内の画像データの組と最小歪にて対応する
ベクトルをそれぞれ選択し、この選択されたベクトルを
示す情報を各ブロックと対応させて符号化するようにし
たものである。As one of image data compression methods,
A method utilizing vector quantization is known. This method divides the two-dimensional image data into blocks of K samples and predefines K vector elements to create a codebook composed of a plurality of different vectors. The vector corresponding to the set of image data and the minimum distortion is selected, and the information indicating the selected vector is encoded in association with each block.
【0005】上述のようなブロック内の画像データは互
いに高い相関性を有しているので、各ブロック内の画像
データを、比較的少数だけ用意したベクトルのうちの1
つを用いてかなり正確に示すことが可能となる。したが
って、画像データの伝送あるいは記録は、実際のデータ
の代わりにこのベクトルを示す符号を伝送あるいは記憶
することによってなし得るから、データ圧縮が実現され
るのである。例えば256 レベル(=8bit)の濃度ス
ケールの中間調画像における64画素についての画像デー
タ量は、8×64=512bitとなるが、この64画素を
1ブロックとして該ブロック内の各画像データを64要素
からなるベクトルで表わし、このようなベクトルを256
通り用意したコードブックを作成するものとすれば、1
ブロック当りのデータ量はベクトル識別のためのデータ
量すなわち8bitとなり、結局データ量を8/(8×
64)=1/64に圧縮可能となる。Since the image data in the blocks as described above have a high correlation with each other, one of the vectors in which a relatively small number of image data in each block are prepared is prepared.
It is possible to give a fairly accurate representation by using one. Therefore, transmission or recording of image data can be performed by transmitting or storing a code indicating this vector instead of actual data, so that data compression is realized. For example, the amount of image data for 64 pixels in a halftone image of a 256-level (= 8 bit) density scale is 8 × 64 = 512 bits, and these 64 pixels are regarded as one block, and each image data in the block has 64 elements. Is represented by a vector consisting of
Assuming that you create a codebook prepared as above, 1
The data amount per block is the data amount for vector identification, that is, 8 bits, and the data amount is 8 / (8 ×
64) = 1/64 can be compressed.
【0006】以上のようにして画像データを圧縮して記
録あるいは伝送した後、ベクトル識別情報が示すベクト
ルのベクトル要素を各ブロック毎の再構成データとし、
この再構成データを用いれば原画像が再現される。After the image data is compressed and recorded or transmitted as described above, the vector element of the vector indicated by the vector identification information is reconstructed data for each block,
The original image is reproduced by using the reconstructed data.
【0007】また、上述した予測符号化によるデータ圧
縮を行う場合の圧縮率を向上させる方法の1つとして、
予測符号化処理と共に画像データのビット分解能(濃度
分解能)を低下させる、すなわち画像データをより粗く
量子化する量子化処理を行うことが考えられる。Further, as one of the methods for improving the compression rate in the case of performing the data compression by the above-mentioned predictive coding,
It is conceivable to reduce the bit resolution (density resolution) of the image data together with the predictive coding process, that is, perform the quantization process of coarsely quantizing the image data.
【0008】そこで、本願出願人により、上述した予測
符号化による方法と量子化による方法とを組み合わせた
補間符号化による画像データ圧縮方法が提案されている
(特開昭62−247676号公報)。この方法は、画像データ
を適当な間隔でサンプリングした主データと該主データ
以外の補間データとに区分し、補間データは上記主デー
タに基づいて内挿予測符号化処理、すなわち補間データ
を主データに基づいて内挿予測し、予測誤差に対してハ
フマン符号化等の可変長符号化(値により符号長が変わ
るような信号への変換)を行うことにより画像データを
圧縮するものである。Therefore, the applicant of the present application has proposed an image data compression method by interpolation coding, which is a combination of the method by predictive coding and the method by quantization described above (Japanese Patent Application Laid-Open No. 62-247676). In this method, image data is divided into main data sampled at appropriate intervals and interpolation data other than the main data, and the interpolation data is interpolation prediction coding processing based on the main data, that is, interpolation data is the main data. The image data is compressed by performing interpolative prediction based on the above, and performing variable length coding (conversion into a signal whose code length changes depending on the value) such as Huffman coding with respect to the prediction error.
【0009】また、画像データを圧縮するにあたっては
当然圧縮率は高い方が望ましい。しかしながら、上記補
間符号化において大きな圧縮率の向上を望むことは技術
的に困難であり、従ってより大きな圧縮率を達成するた
め、空間分解能を小さくする画像データ数減少処理を上
記補間符号化と組合わせることが考えられる。Further, when compressing image data, it is naturally desirable that the compression rate is high. However, it is technically difficult to expect a large improvement in the compression rate in the above-mentioned interpolation coding, and therefore, in order to achieve a higher compression rate, the image data number reduction processing for reducing the spatial resolution is combined with the above-mentioned interpolation coding. It is possible to match.
【0010】そこで本願出願人により、上述した補間符
号化と画像データ数減少処理とを組み合わせ、より高画
質を維持しつつより高い圧縮率を達成し得る画像データ
圧縮方法が提案されている(特開平2-280462号公報)。Therefore, the applicant of the present application has proposed an image data compression method capable of achieving a higher compression rate while maintaining higher image quality by combining the above-described interpolation coding and the image data number reduction process (special feature). Kaihei 2-280462).
【0011】一方、上述した画像データを処理するため
の方法としてウェーブレット変換なる方法が提案されて
いる。On the other hand, a method called wavelet transform has been proposed as a method for processing the above-mentioned image data.
【0012】ここで、ウェーブレット変換について説明
する。Here, the wavelet transform will be described.
【0013】ウェーブレット変換は、周波数解析の方法
として近年開発されたものであり、ステレオのパターン
マッチング、データ圧縮等に応用がなされているもので
ある(OLIVIER RIOUL and MARTIN VETTERLI;Wavelets a
nd Signal Processing,IEEESP MAGAZINE,P.14-38,OCTOB
ER 1991、Stephane Mallat;Zero-Crossings of a Wavel
et Transform,IEEE TRANSACTIONS ON INFORMATION THEO
RY,VOL.37,NO.4,P.1019-1033,JULY 1991 )。The wavelet transform has been recently developed as a method of frequency analysis and has been applied to stereo pattern matching, data compression, etc. (OLIVIER RIOUL and MARTIN VETTERLI; Wavelets a
nd Signal Processing, IEEESP MAGAZINE, P.14-38, OCTOB
ER 1991, Stephane Mallat; Zero-Crossings of a Wavel
et Transform, IEEE TRANSACTIONS ON INFORMATION THEO
RY, VOL.37, NO.4, P.1019-1033, JULY 1991).
【0014】このウェーブレット変換は、図7に示すよ
うな関数hを基底関数として、This wavelet transform uses a function h as shown in FIG. 7 as a basis function.
【0015】[0015]
【数1】 [Equation 1]
【0016】なる式において信号を複数の周波数帯域毎
の周波数信号に変換するため、フーリエ変換のような偽
振動の問題が発生しない。すなわち、関数hの周期およ
び縮率を変化させ、原信号を移動させることによりフィ
ルタリング処理を行えば、細かな周波数から粗い周波数
までの所望とする周波数に適合した周波数信号を作成す
ることができる。例えば、図8に示すように、信号Sorg
をウェーブレット変換し、各周波数帯域毎に逆ウェーブ
レット変換した信号と、図9に示すように信号Sorg を
フーリエ変換し、各周波数帯域毎に逆フーリエ変換した
信号で見てみると、ウェーブレット変換はフーリエ変換
と比べて原信号Sorg の振動と対応した周波数帯域の周
波数信号を得ることができる。すなわち、フーリエ変換
において原信号Sorg の部分Bと対応する周波数帯域7
の部分B′には振動が発生しているのに対し、ウェーブ
レット変換では原信号Sorg の部分Aと対応する周波数
帯域W7の部分A′には原信号と同様に振動は発生して
いないものとなる。Since the signal is converted into a frequency signal for each of a plurality of frequency bands in the equation, the problem of false vibration unlike the Fourier transform does not occur. That is, if the filtering process is performed by changing the period and reduction ratio of the function h and moving the original signal, it is possible to create a frequency signal adapted to a desired frequency from a fine frequency to a coarse frequency. For example, as shown in FIG. 8, the signal Sorg
Is wavelet-transformed and the inverse wavelet transform is performed for each frequency band, and the signal Sorg is Fourier-transformed as shown in FIG. 9 and the signal is inverse Fourier-transformed for each frequency band. Compared with the conversion, it is possible to obtain the frequency signal in the frequency band corresponding to the vibration of the original signal Sorg. That is, in the Fourier transform, the frequency band 7 corresponding to the part B of the original signal Sorg
While the vibration is generated in the portion B'of, the wavelet transform does not generate the vibration in the portion A'of the frequency band W7 corresponding to the portion A of the original signal Sorg like the original signal. Become.
【0017】また、このウェーブレット変換を用いて、
前述した画像データの圧縮を行う方法が提案されている
(Marc Antonini et al.,Image Coding Using Wavelet
Transform,IEEE TRANSACTIONS ON IMAGE PROCESSING
,VOL.1 ,NO.2,p205-220,APRIL 1992)。Further, using this wavelet transform,
A method for compressing the above-mentioned image data has been proposed (Marc Antonini et al., Image Coding Using Wavelet.
Transform, IEEE TRANSACTIONS ON IMAGE PROCESSING
, VOL.1, NO.2, p205-220, APRIL 1992).
【0018】この方法は、画像を表す原画像データにウ
ェーブレット変換を施して、原画像データを複数の周波
数帯域を有する画像データに変換し、この各画像データ
に対してノイズ成分を多く担持する高周波数帯域の画像
データにはビット数を少なく、主要被写体を担持する低
周波数帯域の画像データにはビット数を多く割り当てて
前述したベクトル量子化を施すことにより、原画像デー
タの圧縮を行うものである。この方法によれば、原画像
データの圧縮率を向上させることができ、また、圧縮さ
れた画像データに対して逆ウェーブレット変換を施すこ
とにより、原画像を完全に復元することができる。According to this method, original image data representing an image is subjected to wavelet transform to convert the original image data into image data having a plurality of frequency bands, and each image data has a high noise component. The image data in the frequency band has a small number of bits, and the image data in the low frequency band that carries the main subject has a large number of bits and is subjected to the aforementioned vector quantization to compress the original image data. is there. According to this method, the compression rate of the original image data can be improved, and the original image can be completely restored by performing the inverse wavelet transform on the compressed image data.
【0019】[0019]
【発明が解決しようとする課題】しかしながら、上述し
たウェーブレット変換を用いて画像データを圧縮する方
法においては、さらに圧縮率を向上させようとすると、
原画像の画質が劣化してしまうおそれがあり、画像の高
圧縮率化には限界があった。However, in the method of compressing image data using the above-mentioned wavelet transform, if the compression rate is further improved,
There is a possibility that the image quality of the original image may be deteriorated, and there is a limit to increase the compression rate of the image.
【0020】本発明は上記事情に鑑み、原画像の画質を
劣化させることなく高い圧縮率により画像データの圧縮
を行うことができる画像データ圧縮処理方法を提供する
ことを目的とするものである。In view of the above circumstances, it is an object of the present invention to provide an image data compression processing method capable of compressing image data at a high compression rate without degrading the image quality of the original image.
【0021】[0021]
【課題を解決するための手段】本発明による画像データ
圧縮処理方法は、所定の被写体を含む画像を表す原画像
データに圧縮処理を施す画像データ圧縮処理方法におい
て、前記原画像データにウェーブレット変換を施すこと
により、前記原画像データを複数の周波数帯域毎の画像
データに分解し、該複数の周波数帯域毎の画像データの
うち、所定の周波数帯域の画像データにおいて該画像デ
ータの非連結成分の値を0とする処理を行い、該0とす
る処理が行われた全ての前記画像データおよび前記0と
する処理が行われなかった全ての前記画像データを符号
化することを特徴とするものである。An image data compression processing method according to the present invention is an image data compression processing method for compressing original image data representing an image including a predetermined subject, wherein the original image data is wavelet transformed. By performing the above, the original image data is decomposed into image data for each of a plurality of frequency bands, and among the image data for each of the plurality of frequency bands, the value of the unconnected component of the image data in the image data of a predetermined frequency band. Is performed, and all of the image data that has been subjected to the processing of 0 and all the image data that has not been subjected to the processing of 0 is encoded. .
【0022】ここで、非連結成分とは、孤立点すなわ
ち、値が0となる画素の中に1つだけ値を有する画素が
ある場合のその値を有する画素のことをいう。Here, the non-connected component means an isolated point, that is, a pixel having a value of 0 when there is a pixel having only one value.
【0023】さらに、本発明による画像データの再構成
方法は、本発明による画像データ圧縮処理方法により得
られた圧縮画像データを再構成するための方法であり、
前記符号化された前記各画像データを逆符号化し、該逆
符号化された前記各画像データに逆ウェーブレット変換
を施すことにより前記原画像データを再構成することを
特徴とするものである。Further, the image data reconstruction method according to the present invention is a method for reconstructing compressed image data obtained by the image data compression processing method according to the present invention.
The original image data is reconstructed by inversely encoding each of the encoded image data and applying an inverse wavelet transform to each of the inversely encoded image data.
【0024】[0024]
【作用および発明の効果】画像データをウェーブレット
変換することにより得られた複数の周波数帯域毎の画像
データにおいて、ノイズ成分の画像データを担持する画
素の位置は、孤立点、すなわち非連結成分となるもので
あり、これを除去しても画質にはそれほど影響を及ぼさ
ないものである。また、画像データを圧縮する際には、
画像データの値が0である点が多いほど符号化を行う際
のビット数を減少させることができるため圧縮率を向上
させることができる。本願発明はこの点に着目してなさ
れたものである。In the image data for each of a plurality of frequency bands obtained by wavelet transforming the image data, the position of the pixel carrying the image data of the noise component becomes an isolated point, that is, a non-connected component. However, the removal of this does not significantly affect the image quality. Also, when compressing image data,
The more the image data value is 0, the more the number of bits for encoding can be reduced, so that the compression rate can be improved. The present invention has been made paying attention to this point.
【0025】すなわち、本発明による画像データ圧縮処
理方法は、ウェーブレット変換により複数の周波数帯域
に分解された画像データのうち、所定の周波数帯域にお
いて、画像データの非連結成分を0とする処理を行うよ
うにしたものである。That is, the image data compression processing method according to the present invention performs processing for setting the non-connected component of the image data to 0 in a predetermined frequency band among the image data decomposed into a plurality of frequency bands by the wavelet transform. It was done like this.
【0026】このように、非連結成分を0とする処理を
行うことにより、画像のノイズ等の不要な情報を除去す
ることができ、これらの不要な情報を符号化する際の必
要なビット数の低減を図ることができ、あるいは符号化
する必要がなくなる。これにより、画像データを圧縮す
る際の圧縮率を向上させることができる。As described above, by performing the processing for setting the non-connected components to 0, unnecessary information such as image noise can be removed, and the number of bits required for encoding these unnecessary information. Can be reduced or need not be encoded. This can improve the compression rate when compressing the image data.
【0027】また、非連結成分を0とする処理を行うよ
うにしたため、有用な情報を担持する画素の画像データ
の値を0とする処理を行うことがなく、有用な情報が欠
落することもない。Further, since the processing for setting the non-connected component to 0 is performed, the processing for setting the value of the image data of the pixel carrying useful information to 0 is not performed, and useful information may be lost. Absent.
【0028】さらに、画像のノイズ等の不要な情報を除
去することができるため、再現された画像の画質を向上
させることができる。Further, since unnecessary information such as image noise can be removed, the quality of the reproduced image can be improved.
【0029】また、本発明による画像データの再構成方
法は、本発明による画像データ圧縮処理方法により圧縮
された画像データを逆符号化すなわち復号化し、この逆
復号化された画像データに対して逆ウェーブレット変換
を施すようにしたため、画像の各部分のうち重要な部分
の画質を維持しつつ原画像を再生することができる。The image data reconstruction method according to the present invention is such that the image data compressed by the image data compression processing method according to the present invention is inversely encoded, that is, decoded, and the inversely decoded image data is inversely encoded. Since the wavelet transform is applied, it is possible to reproduce the original image while maintaining the image quality of an important part of each part of the image.
【0030】[0030]
【実施例】以下図面を参照して本発明の実施例について
説明する。Embodiments of the present invention will be described below with reference to the drawings.
【0031】図1は本発明による画像データ圧縮処理方
法の第1実施例の基本的概念を表す図である。図1に示
すように本発明による画像データ圧縮処理方法の実施例
は、原画像を表す原画像データ1に対してウェーブレッ
ト変換2を施し、複数の周波数帯域毎の画像データ3を
得る。次いで、各画像データ3のうち所定の周波数帯域
よりも高い周波数帯域の画像データ3において、この画
像データの非連結成分の値を0とする処理4を各周波数
帯域毎に行う。そして、この0とする処理4がなされた
全ての画像データ3および0とする処理4が行われなか
った全ての画像データ3に対して符号化5を行うもので
ある。FIG. 1 is a diagram showing the basic concept of a first embodiment of an image data compression processing method according to the present invention. As shown in FIG. 1, in the embodiment of the image data compression processing method according to the present invention, original image data 1 representing an original image is subjected to wavelet transform 2 to obtain image data 3 for each of a plurality of frequency bands. Next, in the image data 3 in the frequency band higher than the predetermined frequency band among the image data 3, the process 4 for setting the value of the non-connected component of the image data to 0 is performed for each frequency band. Then, the coding 5 is performed on all the image data 3 for which the processing 4 for setting 0 is performed and for all the image data 3 for which the processing 4 for setting 0 is not performed.
【0032】以下本発明の実施例の詳細について説明す
る。The details of the embodiments of the present invention will be described below.
【0033】本実施例は、例えば特開昭55-12492号公報
や特開昭56-11395号等に記録されている蓄積性蛍光体シ
ートを利用した放射線画像情報記録再生システムにおい
て、蓄積性蛍光体シートに記録された人体の放射線画像
をレーザビーム走査によりデジタル画像データとして読
み取ったものを対象としている。なお、放射線画像の読
み取りは、図2に示すように、蓄積性蛍光体シート10に
対して主走査方向(横方向)にレーザビームを走査させ
ながらシート10を副走査方向(縦方向)に移動させてシ
ート10を2次元走査することにより行われたものであ
る。This embodiment is a radiation image information recording / reproducing system using a stimulable phosphor sheet recorded in JP-A-55-12492 and JP-A-56-11395. It is intended to read a radiation image of a human body recorded on a body sheet as digital image data by laser beam scanning. As shown in FIG. 2, the radiation image is read by moving the sheet 10 in the sub-scanning direction (vertical direction) while scanning the stimulable phosphor sheet 10 with a laser beam in the main scanning direction (horizontal direction). Then, the sheet 10 is two-dimensionally scanned.
【0034】次いで、原画像データに対してウェーブレ
ット変換がなされる。Next, wavelet transform is performed on the original image data.
【0035】図3は、原画像データSorg に対するウェ
ーブレット変換の詳細を表す図である。FIG. 3 is a diagram showing details of the wavelet transform for the original image data Sorg.
【0036】なお、本実施例においては、ウェーブレッ
ト変換の各係数が直交する直交ウェーブレット変換を行
うものであり、前述したMarc Antonini らの文献に記載
されているものである。In the present embodiment, the orthogonal wavelet transform in which each coefficient of the wavelet transform is orthogonal is performed, which is described in the above-mentioned document of Marc Antonini et al.
【0037】図3に示すように、原画像データSorg の
主走査方向に基本ウェーブレット関数より求められる関
数gと関数hとによりフィルタリング処理を行う。すな
わち、このような関数g,hによる主走査方向に並ぶ画
素の一列毎のフィルタリング処理を副走査方向に一画素
ずつズラしながら行い、原画像データSorg の主走査方
向のウェーブレット変換係数信号Wg0,Wh0を求めるも
のである。As shown in FIG. 3, the filtering process is performed in the main scanning direction of the original image data Sorg by the function g and the function h obtained from the basic wavelet function. That is, the filtering process for each column of pixels arranged in the main scanning direction by such functions g and h is performed while shifting each pixel in the sub scanning direction to obtain the wavelet transform coefficient signal Wg0 of the original image data Sorg in the main scanning direction. This is for obtaining Wh0.
【0038】ここで、関数g,hは基本ウェーブレット
関数より一意に求められるものであり、例えば、関数h
は、以下の表1に示すものとなる。なお、表1において
関数h′は、ウェーブレット変換がなされた画像データ
に逆ウェーブレット変換を行う際に用いる関数を表すも
のである。また以下の式(2) に示すように関数gは関数
h′から求められ、逆ウェーブレット変換を行うための
関数g′は関数hから求められる。Here, the functions g and h are uniquely obtained from the basic wavelet function. For example, the function h
Is shown in Table 1 below. It should be noted that in Table 1, the function h'represents a function used when performing inverse wavelet transform on image data that has been subjected to wavelet transform. Further, as shown in the following expression (2), the function g is obtained from the function h ', and the function g'for performing the inverse wavelet transform is obtained from the function h.
【0039】[0039]
【表1】 [Table 1]
【0040】 g′=(−1)n h g =(−1)n h′ …(2) このようにして、ウェーブレット変換係数信号Wg0、W
h0が求められると、ウェーブレット変換係数信号Wg0、
Wh0について、主走査方向の画素を1画素おきに間引
き、主走査方向の画素数を1/2 にする。ついで、この画
素が間引かれたウェーブレット変換係数信号Wg0、Wh0
それぞれの副走査方向に関数g,hによりフィルタリン
グ処理を行い、ウェーブレット変換係数信号WW0 ,W
V0 ,VW0 およびVV0 を得る。G ′ = (− 1) n h g = (− 1) n h ′ (2) In this way, the wavelet transform coefficient signals Wg 0, W
When h0 is obtained, the wavelet transform coefficient signal Wg0,
For Wh0, pixels in the main scanning direction are thinned out every other pixel, and the number of pixels in the main scanning direction is halved. Then, the wavelet transform coefficient signals Wg0 and Wh0 from which this pixel is thinned out.
Filtering processing is performed in each sub-scanning direction by the functions g and h, and wavelet transform coefficient signals WW 0 and W
Obtain V 0 , VW 0 and VV 0 .
【0041】次いでウェーブレット変換係数信号W
W0 ,WV0 ,VW0 およびVV0 について、副走査方
向の画素を1画素おきに間引くことを行い、副走査方向
の画素数を1/2 とする処理を行う。これにより、各ウェ
ーブレット変換係数信号VV0 ,WV0 ,VW0 ,WW
0 の画素数は原画像データSorg の画素数の1/4 とな
る。次いで、ウェーブレット変換係数信号VV0 の主走
査方向に関数g,hによりフィルタリング処理を行う。Next, the wavelet transform coefficient signal W
With respect to W 0 , WV 0 , VW 0 and VV 0 , the pixels in the sub-scanning direction are thinned out every other pixel, and the number of pixels in the sub-scanning direction is halved. As a result, each wavelet transform coefficient signal VV 0 , WV 0 , VW 0 , WW
The number of pixels of 0 is 1/4 of the number of pixels of the original image data Sorg. Next, filtering processing is performed by the functions g and h in the main scanning direction of the wavelet transform coefficient signal VV 0 .
【0042】すなわち、関数g0 ,h0 により主走査方
向に並ぶ画素の一列毎のフィルタリング処理を副走査方
向に一画素づつズラながら行い、ウェーブレット変換係
数信号VV0 の主走査方向のウェーブレット変換係数信
号Wg1およびWh1を求めるものである。That is, the filtering processing for each column of pixels arranged in the main scanning direction by the functions g 0 and h 0 is performed while shifting the pixels one pixel at a time in the sub-scanning direction to obtain the wavelet transform coefficient of the wavelet transform coefficient signal VV 0 in the main scanning direction. The signals Wg1 and Wh1 are obtained.
【0043】ここでウェーブレット変換係数信号VV0
は主副両方向について画素数が原画像データの1/2 とな
っているため、画像の周波数帯域は原画像データと比較
して半分となっている。したがって、ウェーブレット変
換係数信号VV0 を関数g,hでフィルタリング処理を
施すことにより、原画像データの周波数成分のうちウェ
ーブレット変換係数信号VV0 が表す周波数成分よりも
低周波数成分を表すウェーブレット変換係数信号Wg1,
Wh1が求められる。Here, the wavelet transform coefficient signal VV 0
Since the number of pixels in both main and sub directions is half that of the original image data, the frequency band of the image is half that of the original image data. Therefore, by filtering the wavelet transform coefficient signal VV 0 with the functions g and h, the wavelet transform coefficient signal representing a frequency component lower than the frequency component represented by the wavelet transform coefficient signal VV 0 among the frequency components of the original image data. Wg1,
Wh1 is required.
【0044】このようにして、ウェーブレット変換係数
信号Wg1,Wh1が求められると、ウェーブレット変換係
数信号Wg1,Wh1について、主走査方向の画素を1画素
おきに間引き、主走査方向の画素数をさらに1/2 とす
る。次いでウェーブレット変換係数信号Wg1、Wh1それ
ぞれの副走査方向に関数g,hによりフィルタリング処
理を行い、ウェーブレット変換係数信号WW1 ,W
V1 ,VW1 およびVV1 を得る。When the wavelet transform coefficient signals Wg1 and Wh1 are obtained in this way, the wavelet transform coefficient signals Wg1 and Wh1 are thinned out every other pixel in the main scanning direction, and the number of pixels in the main scanning direction is further reduced to 1. / 2 Then the wavelet transform factor signals Wg1, Wh1 function g in each of the sub-scanning direction, performs a filtering process by h, the wavelet transform factor signals WW 1, W
Obtain V 1 , VW 1 and VV 1 .
【0045】次いでウェーブレット変換係数信号W
W1 ,WV1 ,VW1 ,VV1 について、副走査方向の
画素を1画素おきに間引き、副走査方向の画素数を1/2
とする処理を行う。これにより、各ウェーブレット変換
係数信号VV1 ,WV1 ,VW1,WW1 の画素数は原
画像データSorg の画素数の1/16となる。Next, the wavelet transform coefficient signal W
For W 1 , WV 1 , VW 1 , and VV 1 , pixels in the sub-scanning direction are thinned out every other pixel, and the number of pixels in the sub-scanning direction is halved.
And perform the process. As a result, the number of pixels of each wavelet transform coefficient signal VV 1 , WV 1 , VW 1 , WW 1 becomes 1/16 of the number of pixels of the original image data Sorg.
【0046】以下、上述したのと同様にして、画素が間
引かれたウェーブレット変換係数信号VV1 の主走査方
向に関数g,hによりフィルタリング処理を行い、さら
に得られたウェーブレット変換係数信号の主走査方向の
画素を間引き、この画素を間引いたウェーブレット変換
係数信号について、副走査方向に関数g,hによりフィ
ルタリング処理を行い、ウェーブレット変換係数信号W
W2 ,WV2 ,VW2,VV2 を得る。Thereafter, in the same manner as described above, filtering processing is performed by the functions g and h in the main scanning direction of the wavelet transform coefficient signal VV 1 from which pixels are thinned, and the main wavelet transform coefficient signal obtained is further processed. Pixels in the scanning direction are thinned out, and the wavelet transform coefficient signal obtained by thinning out the pixels is filtered in the sub-scanning direction by the functions g and h to obtain the wavelet transform coefficient signal W.
W 2 , WV 2 , VW 2 and VV 2 are obtained.
【0047】このようなウェーブレット変換をN回繰り
返すことによりウェーブレット変換係数信号WW0 〜W
WN ,WV0 〜WVN ,VW0 〜VWN ,およびVVN
を得る。ここで、N回目のウェーブレット変換により得
られるウェーブレット変換係数信号WWN ,WVN ,V
WN ,VVN は、原画像データと比較して主副両方向の
画素数が各々(1/2)N となっているため、各ウェーブレ
ット変換係数信号はNが大きいほど周波数帯域が低く、
原画像データの周波数成分のうち低周波成分を表すデー
タとなる。By repeating such wavelet transform N times, wavelet transform coefficient signals WW 0 to W
W N, WV 0 ~WV N, VW 0 ~VW N, and VV N
To get Here, the wavelet transform factor signals obtained by the wavelet transform of the N-th WW N, WV N, V
Since W N and VV N each have (1/2) N pixels in both the main and sub directions as compared with the original image data, the larger the N of each wavelet transform coefficient signal, the lower the frequency band,
It becomes the data representing the low frequency component of the frequency component of the original image data.
【0048】したがって、ウェーブレット変換係数信号
WWi (i=0〜N、以下同様)は、原画像データSor
g の主副両方向の周波数の変化を表すものであり、iが
大きいほど低周波信号となる。またウェーブレット変換
係数信号WVi は画像信号Sorg の主走査方向の周波数
の変化を表すものであり、iが大きいほど低周波信号と
なる。さらにウェーブレット変換係数信号VWi は画像
信号Sorg の副走査方向の周波数の変化を表すものであ
り、iが大きいほど低周波信号となる。Therefore, the wavelet transform coefficient signal WW i (i = 0 to N, the same applies hereinafter) is used as the original image data Sor.
It represents a change in the frequency of g in both main and sub directions, and the larger i is, the lower the frequency of the signal becomes. The wavelet transform coefficient signal WV i represents the change in the frequency of the image signal Sorg in the main scanning direction, and the larger i is, the lower the frequency signal becomes. Furthermore, the wavelet transform coefficient signal VW i represents the change in frequency of the image signal Sorg in the sub-scanning direction, and the larger i is, the lower the frequency signal becomes.
【0049】ここで、図4にウェーブレット変換係数信
号を複数の周波数帯域毎に表す図を示す。なお、図4に
おいては便宜上3回目のウェーブレット変換を行った状
態までを表すものとする。なお、図4においてウェーブ
レット変換係数信号WW3 は原画像を主副各方向につい
て(1/2)3 に縮小したものとなっている。FIG. 4 is a diagram showing the wavelet transform coefficient signal for each of a plurality of frequency bands. Note that, in FIG. 4, for convenience, the state up to the third wavelet transform is shown. In FIG. 4, the wavelet transform coefficient signal WW 3 is the original image reduced to (1/2) 3 in each of the main and sub directions.
【0050】次いで、各ウェーブレット変換係数信号W
Vi ,VWi ,WWi ,VVi の係数信号中の所定の周
波数帯域の係数信号において非連結成分の値を0とする
処理を行う。Next, each wavelet transform coefficient signal W
A process of setting the value of the non-connected component to 0 in the coefficient signal of a predetermined frequency band among the coefficient signals of V i , VW i , WW i , and VV i is performed.
【0051】例えば、図5に示すある周波数帯域の係数
信号において、画素w(i,j)について、その周囲に
ある4点近傍の画素a(w(i ,j-1)),b(w(i+1
,j)),c(w(i ,j+1)),d(w(i-1 ,j))の
値が0であった場合に、画素w(i,j)の値を0とす
る処理を行うものである。For example, in the coefficient signal of a certain frequency band shown in FIG. 5, for a pixel w (i, j), pixels a (w (i, j-1)) and b (w in the vicinity of four points around it. (I + 1
, J)), c (w (i, j + 1)), d (w (i-1, j)) has a value of 0, the value of the pixel w (i, j) is set to 0. The processing is performed.
【0052】この処理を、前述した係数信号中の必要な
情報が欠落していない程度の低い周波数帯域よりも高い
周波数帯域の全てのウェーブレット変換係数信号W
Vi ,VWi ,WWi ,VVi について行う。This processing is performed for all wavelet transform coefficient signals W in the frequency band higher than the low frequency band in which necessary information in the coefficient signal is not missing.
V i, VW i, WW i , carried out on VV i.
【0053】次いでこのような処理がなされたウェーブ
レット変換係数信号WVi ,VWi,WWi ,VVi に
ついて前述したハフマン符号化、予測符号化等の符号化
を行うことにより圧縮処理がなされる。Then, the wavelet transform coefficient signals WV i , VW i , WW i , and VV i thus processed are subjected to the above-described Huffman coding, predictive coding, etc. to perform compression processing.
【0054】この際ウェーブレット変換係数信号中のノ
イズ等と見られる非連結成分の値は0とされているた
め、ノイズ等の分だけ量子化する際のビット数の低減を
図ることができる。At this time, the value of the non-connected component, which is considered to be noise or the like, in the wavelet transform coefficient signal is set to 0, so that it is possible to reduce the number of bits for quantization by the amount of noise or the like.
【0055】このため、各ウェーブレット変換係数信号
WVi ,VWi ,WWi ,VVi のそれぞれについてノ
イズ等とみられる不要な情報を除去し、符号化する際の
ビット数を低減するとともに、必要な情報を有する部分
は、通常の符号化を行う際のビット数によりある程度の
画質を維持しつつ圧縮率を向上させることができる。Therefore, unnecessary information, such as noise, is removed from each of the wavelet transform coefficient signals WV i , VW i , WW i , and VV i to reduce the number of bits at the time of encoding, and necessary. In the portion having information, it is possible to improve the compression rate while maintaining a certain level of image quality depending on the number of bits used for normal encoding.
【0056】このように符号化がなされて圧縮された原
画像データSorg は例えば光ディスク等の記録媒体に格
納され、保存、移送等がなされる。なお、符号化につい
ては、特開昭61-90531号公報、特開昭62-247676 号公
報、特開昭63-41714号公報、特開平1-101341号公報等に
その説明が記載されているため詳細な説明はここでは省
略する。The original image data Sorg encoded and compressed in this way is stored in a recording medium such as an optical disk, and stored and transported. The encoding is described in JP-A-61-90531, JP-A-62-247676, JP-A-63-41714 and JP-A-1-01341. Therefore, detailed description is omitted here.
【0057】次に圧縮されたデータを再構成する方法に
ついて説明する。Next, a method of reconstructing compressed data will be described.
【0058】まず、圧縮された原画像データに対し、復
号化を行うことにより、各ウェーブレット変換係数信号
WVi ,VWi ,WWi ,VVi を得る。この復号化に
ついては、前述した特開昭61-90531号公報等にその詳細
が記載されているため、その詳細な説明はここでは省略
する。First, the compressed original image data is decoded to obtain the respective wavelet transform coefficient signals WV i , VW i , WW i , VV i . The details of this decryption are described in the above-mentioned Japanese Patent Laid-Open No. 61-90531, and the detailed description thereof is omitted here.
【0059】次いで、復号化がなされることにより得ら
れたウェーブレット変換係数信号WVi ,VWi ,WW
i ,VVi について逆ウェーブレット変換を施す。Next, the wavelet transform coefficient signals WV i , VW i , WW obtained by decoding are performed.
i, performs inverse wavelet transform for VV i.
【0060】図6は、逆ウェーブレット変換の詳細を表
す図である。FIG. 6 is a diagram showing details of the inverse wavelet transform.
【0061】図6に示すように、まず各ウェーブレット
変換係数信号VVN ,VWN ,WVN ,WWN について
副走査方向に並ぶ画素間に1画素分の間隔をあける処理
を行う(図では×2と表示)。次いでこの間隔があけら
れたウェーブレット変換係数信号VVN を副走査方向に
前述した関数hとは異なる関数h′により、ウェーブレ
ット変換係数信号VWN を副走査方向に前述した関数g
とは異なる関数g′によりフィルタリング処理を行う。
すなわち、関数g′,h′によるウェーブレット変換係
数信号VVN ,VWN の副走査方向に並ぶ一列の画素毎
のフィルタリング処理を主走査方向に一画素ずつズラし
ながら行い、ウェーブレット変換係数信号VVN ,VW
N の逆ウェーブレット変換係数信号を得、これを2倍し
て加算することにより逆ウェーブレット変換係数信号W
hN′を得る。[0061] As shown in FIG. 6, the wavelet transform factor signals VV N, VW N, WV N , the WW N performs processing spacing of one pixel among pixels arranged in the sub-scanning direction (× in FIG. 2). Followed by different function h 'is the function h described above wavelet transform factor signals VV N this spaced in the sub-function g mentioned above the wavelet transform factor signal VW N in the sub-scanning direction
Filtering processing is performed by a function g ′ different from.
That is, the function g ', h' wavelet transform factor signals by VV N, filtering processing for each pixel of a row aligned in the sub-scanning direction VW N in the main scanning direction is performed while Shifts pixel by pixel, the wavelet transform factor signals VV N , VW
The inverse wavelet transform coefficient signal W is obtained by doubling and adding N inverse wavelet transform coefficient signals.
get hN ′.
【0062】このようにウェーブレット変換を行う関数
と逆ウェーブレット変換を行う関数とを異なるものとし
ているのは、以下のような理由からである。ウェーブレ
ット変換と逆ウェーブレット変換で同一の関数となる、
すなわち、直交する関数を設計することは難しく、直交
性、連続性、関数の短さ、対称性のいずれかの条件を緩
める必要がある。そこで、直交性の条件を緩めることに
より他の条件を満たす関数を選択したものである。The function for performing the wavelet transform and the function for performing the inverse wavelet transform are different from each other for the following reason. Wavelet transform and inverse wavelet transform have the same function,
That is, it is difficult to design orthogonal functions, and it is necessary to loosen any of the conditions of orthogonality, continuity, function shortness, and symmetry. Therefore, a function that satisfies other conditions is selected by relaxing the orthogonality condition.
【0063】以上より、本実施例ではウェーブレット変
換を行う関数h,gと逆ウェーブレット変換を行う関数
h′,g′とを双直交の異なるものとしている。したが
って、ウェーブレット変換係数信号VVi ,VWi ,W
Vi ,WWi を関数h′,g′で逆ウェーブレット変換
することにより、原画像データを完全に復元できること
となる。As described above, in the present embodiment, the functions h and g for performing the wavelet transform and the functions h'and g'for performing the inverse wavelet transform are made bi-orthogonal different from each other. Therefore, the wavelet transform coefficient signals VV i , VW i , W
By inverse wavelet transforming V i and WW i with the functions h ′ and g ′, the original image data can be completely restored.
【0064】一方、これと並列して、ウェーブレット変
換係数信号WVN を副走査方向に関数h′により、ウェ
ーブレット変換係数信号WWN を副走査方向に関数g′
によりフィルタリング処理を行い、ウェーブレット変換
係数信号WVN ,WWN の逆ウェーブレット変換係数信
号を得、これを2倍して加算することにより逆ウェーブ
レット変換係数信号WgN′を得る。On the other hand, in parallel with this, the wavelet transform coefficient signal WV N in the sub-scanning direction is given by the function h ', and the wavelet transform coefficient signal WW N is given in the sub-scanning direction by the function g'.
The performs filtering process, to obtain an inverse wavelet transform factor signals of the wavelet transform factor signals WV N, WW N, obtain the inverse wavelet transform factor signals WGN 'by adding to 2 times this.
【0065】次いで、逆ウェーブレット変換係数信号W
hN′,WgN′について主走査方向に並ぶ画素間に1画素
分の間隔をあける処理を行う。その後逆ウェーブレット
変換係数信号WhN′を主走査方向に関数h′により、逆
ウェーブレット変換係数信号WgN′を主走査方向に関数
g′によりフィルタリング処理し、ウェーブレット変換
係数信号WhN′,WgN′の逆ウェーブレット変換係数信
号を得、これを2倍して加算することにより逆ウェーブ
レット変換係数信号VVN-1 ′を得る。Next, the inverse wavelet transform coefficient signal W
With respect to hN 'and WgN', a process of spacing one pixel between pixels arranged in the main scanning direction is performed. Then, the inverse wavelet transform coefficient signal WhN 'is filtered in the main scanning direction by the function h', and the inverse wavelet transform coefficient signal WgN 'is filtered in the main scanning direction by the function g', and the inverse wavelet of the wavelet transform coefficient signals WhN 'and WgN' is obtained. An inverse wavelet transform coefficient signal VVN -1 'is obtained by obtaining a transform coefficient signal, doubling it, and adding it.
【0066】次いでこの逆ウェーブレット変換係数信号
VVN-1 ′、ウェーブレット変換係数信号VWN-1 ,W
VN-1 ,WWN-1 について副走査方向に並ぶ画素間に1
画素分の間隔をあける処理を行う。その後この逆ウェー
ブレット変換係数信号VVN-1 ′を副走査方向に前述し
た関数h′により、ウェーブレット変換係数信号VW
N-1 を副走査方向に前述した関数g′によりフィルタリ
ング処理を行う。すなわち、関数g′,h′によるウェ
ーブレット変換係数信号VVN-1 ′,VWN-1 の副走査
方向に並ぶ一列の画素毎のフィルタリング処理を主走査
方向に一画素ずつズラしながら行い、ウェーブレット変
換係数信号VVN-1 ′,VWN-1 の逆ウェーブレット変
換係数信号を得、これを2倍して加算することにより逆
ウェーブレット変換係数信号WhN-1′を得る。Next, the inverse wavelet transform coefficient signal VVN -1 'and the wavelet transform coefficient signals VWN -1 , W
1 between pixels arranged in the sub - scanning direction for V N-1 and WW N- 1
Performs a process for spacing pixels. After that, the inverse wavelet transform coefficient signal VVN -1 'is processed in the sub-scanning direction by the above-mentioned function h'.
N-1 is filtered in the sub - scanning direction by the function g'described above. That is, the filtering processing is performed for each row of pixels of the wavelet transform coefficient signals VVN -1 ', VWN -1 by the functions g', h ', which are arranged in the sub - scanning direction, while shifting each pixel in the main-scanning direction. transform factor signals VV N-1 to obtain a ', to give the inverse wavelet transform factor signal VW N-1, the inverse wavelet transform factor signals WhN-1 by adding to 2 times this'.
【0067】一方、これと並列して、ウェーブレット変
換係数信号WVN-1 を副走査方向に関数h′により、ウ
ェーブレット変換係数信号WWN-1 を副走査方向に関数
g′によりフィルタリング処理を行い、ウェーブレット
変換係数信号WVN-1 ,WWN-1 の逆ウェーブレット変
換係数信号を得、これを2倍して加算することにより逆
ウェーブレット変換係数信号WgN-1′を得る。On the other hand, in parallel with this, the wavelet transform coefficient signal WV N-1 is filtered in the sub - scanning direction by the function h ', and the wavelet transform coefficient signal WW N-1 is filtered in the sub - scanning direction by the function g'. , to obtain a wavelet transform factor signals WV N-1, the inverse wavelet transform factor signal WW N-1, to obtain an inverse wavelet transform factor signals WGN-1 'by adding to 2 times this.
【0068】次いで、逆ウェーブレット変換係数信号W
hN-1′,WgN-1′について主走査方向に並ぶ画素間に1
画素分の間隔をあける処理を行う。その後逆ウェーブレ
ット変換係数信号WhN-1′を主走査方向に関数h′によ
り、逆ウェーブレット変換係数信号WgN-1′を主走査方
向に関数g0 ′によりフィルタリング処理し、ウェーブ
レット変換係数信号WhN-1′,WgN-1′の逆ウェーブレ
ット変換係数信号を得、これを2倍して加算することに
より逆ウェーブレット変換係数信号VVN-2 ′を得る。Next, the inverse wavelet transform coefficient signal W
1 between pixels lined up in the main scanning direction for hN-1 'and WgN-1'
Performs a process for spacing pixels. Subsequent inverse wavelet transform factor signals WhN-1 'function in the main scanning direction h', 'functions g 0 in the main scanning direction' inverse wavelet transform factor signals WGN-1 filtering process by the wavelet transform factor signals WhN-1 ', WGN-1' give the inverse wavelet transform factor signals of, obtaining the inverse wavelet transform factor signal VV N-2 'by adding to 2 times this.
【0069】以下、順次逆ウェーブレット変換係数信号
VVi ′(i=−1〜N)を作成し、最終的に逆ウェー
ブレット変換係数信号VV-1′を得る。この最終的な逆
ウェーブレット変換係数信号VV-1′が原画像データS
org を表す画像データとなる。Thereafter, the inverse wavelet transform coefficient signal VV i ′ (i = −1 to N) is sequentially created, and finally the inverse wavelet transform coefficient signal VV −1 ′ is obtained. This final inverse wavelet transform coefficient signal VV −1 ′ is the original image data S
Image data representing org.
【0070】このようにして得られたウェーブレット変
換係数信号VV-1′は図示しない画像再生装置に送られ
て、放射線画像の再生に供せられる。The wavelet transform coefficient signal VV -1 ′ thus obtained is sent to an image reproducing device (not shown) and is used for reproducing a radiation image.
【0071】この再生装置は、CRT等のディスプレイ
手段でもよいし、感光フイルムに光走査記録を行う記録
装置であってもよい。This reproducing device may be a display device such as a CRT or a recording device for performing optical scanning recording on the photosensitive film.
【0072】ここで、上述した本発明による画像データ
圧縮処理方法実施例において、非連結成分を0とする処
理として、図5における画素w(i,j)の4点近傍の
画素a,b,c,dの値が0の場合に画素w(i,j)
を0とする処理を行っているが、8点近傍の画素a,
b,c,d,e,f,g,hの値が0の場合に画素w
(i,j)を0とするようにしてもよい。さらに、ウェ
ーブレット変換係数信号WWi は原画像データSorg の
主副両方向の周波数の変化を示すものであり、ウェーブ
レット変換係数信号WVi は原画像データSorg の主走
査方向の、ウェーブレット変換係数信号VWi は副走査
方向の周波数の変化をそれぞれ示すものであることか
ら、ウェーブレット変換係数信号WWi については画素
w(i,j)の4点近傍の画素a,b,c,dを、ウェ
ーブレット変換係数信号WVi については副走査方向の
2点近傍の画素b,d、ウェーブレット変換係数信号V
Wi については、主走査方向の2点近傍の画素a,cを
用い、係数信号によって異方性をもたせて非連結成分を
0とする処理を行うようにしてもよい。Here, in the above-described embodiment of the image data compression processing method according to the present invention, as processing for setting the non-connected component to 0, the pixels a, b, and 4 in the vicinity of four points of the pixel w (i, j) in FIG. Pixel w (i, j) when the values of c and d are 0
Is set to 0, but pixels a,
Pixel w when the values of b, c, d, e, f, g, and h are 0
(I, j) may be set to 0. Further, the wavelet transform coefficient signal WW i indicates the change in frequency of the original image data Sorg in both main and sub directions, and the wavelet transform coefficient signal WV i is the wavelet transform coefficient signal VW i in the main scanning direction of the original image data Sorg. Represents the change in frequency in the sub-scanning direction, therefore, for the wavelet transform coefficient signal WW i , the pixels a, b, c, d near the four points of the pixel w (i, j) are Regarding the signal WV i , the pixels b and d near two points in the sub-scanning direction and the wavelet transform coefficient signal V
As for W i , pixels a and c near two points in the main scanning direction may be used, and anisotropy may be imparted by the coefficient signal to make the non-connected component zero.
【0073】また、この非連結成分を0とする処理を収
束するまで行うようにしてもよく、さらに、非連結成分
を0とする処理であれば、上述した処理に限らず、いか
なる処理を行うようにしてもよいものである。Further, the processing for setting the non-connected component to 0 may be performed until it converges. Further, as long as the processing for setting the non-connected component to 0, any processing is performed without being limited to the above processing. You can do it.
【0074】なお、上述した本発明による画像データ圧
縮処理方法の実施例においては、ウェーブレット変換を
行うための関数hとして表1に示すものを用いたが、こ
れに限定されるものではなく以下に示す表2、表3に示
すものを用いてもよい。In the above-described embodiment of the image data compression processing method according to the present invention, the function h for performing wavelet transformation shown in Table 1 is used, but the function h is not limited to this. You may use what is shown in Table 2 and Table 3 shown.
【0075】[0075]
【表2】 [Table 2]
【0076】[0076]
【表3】 [Table 3]
【0077】また、これ以外にもウェーブレット変換を
行うことのできる関数であれば、いかなる関数を用いて
もよく、例えば双直交ではなく対称ではないが直交する
ものを用いてもよい。In addition to this, any function may be used as long as it is a function capable of performing wavelet transformation, and for example, a function which is not symmetric but orthogonal but not symmetrical may be used.
【0078】さらに、表1,2および3に示すようにn
=0の軸に関して左右対称な関数のみではなく、n=0
の軸に関して左右非対称な関数を用いてウェーブレット
変換を行うようにしてもよいものである。このように左
右非対称な関数を用いてウェーブレット変換を行った場
合は、ウェーブレット変換を行った関数をn=0の軸に
関して左右を反転させた関数を用いて逆ウェーブレット
変換を行うものである。すなわち、左右非対称な関数
g,hについて、逆ウェーブレット変換を行なう関数
g′,h′は、 g[n]=g′[−n] h[n]=h′[−n] …(3) 但し、[−n]は左右反転を表す。Further, as shown in Tables 1, 2, and 3, n
Not only the function symmetrical about the axis of = 0, but n = 0
The wavelet transform may be performed using a function that is asymmetric with respect to the axis of. When the wavelet transform is performed using the left-right asymmetric function as described above, the inverse wavelet transform is performed by using the function obtained by inverting the left-right function of the wavelet-transformed function with respect to the axis of n = 0. That is, with respect to the left-right asymmetric functions g and h, the functions g ′ and h ′ that perform the inverse wavelet transform are: g [n] = g ′ [− n] h [n] = h ′ [− n] (3) However, [-n] represents left-right inversion.
【0079】となる。It becomes
【0080】さらに、上述した実施例においては、放射
線画像を表す原画像データを圧縮処理する実施例につい
て説明したが、本発明による画像の圧縮処理方法は、通
常の画像についても適用できるものである。Further, in the above-described embodiment, the embodiment in which the original image data representing the radiation image is compressed is described, but the image compression processing method according to the present invention can be applied to a normal image. .
【0081】例えば、主要被写体として人物等が記録さ
れた35mmネガフイルムの画像を圧縮する実施例について
説明すると、まずこのネガフイルムをデジタルスキャナ
ーで読み取り、この画像を表す画像データを得、この画
像データについて前述したような関数g,hによりフィ
ルタリング処理することによりウェーブレット変換を行
う。次いでウェーブレット変換を行うことにより得られ
たウェーブレット変換係数信号に対し、上述した本発明
の実施例と同様に、非連結成分を0とする処理を行う。For example, an example of compressing an image of a 35 mm negative film in which a person or the like is recorded as a main subject will be described. First, the negative film is read by a digital scanner to obtain image data representing this image, and the image data The wavelet transform is performed by performing the filtering process with the functions g and h as described above. Then, the wavelet transform coefficient signal obtained by performing the wavelet transform is subjected to processing for setting the non-connected component to 0, as in the above-described embodiment of the present invention.
【0082】次いで、この0とする処理がなされたウェ
ーブレット変換係数信号に対して符号化を行うことによ
り画像データを圧縮する。Next, the image data is compressed by encoding the wavelet transform coefficient signal which has been processed to be 0.
【0083】また、この圧縮された画像データを前述し
た実施例と同様に復号化し、さらに逆ウェーブレット変
換を施すことにより、原画像データを再構成することが
できる。The original image data can be reconstructed by decoding the compressed image data in the same manner as in the above-mentioned embodiment and further applying the inverse wavelet transform.
【0084】このように、圧縮処理を行うことにより、
0とされたデータの分だけ符号化する際のビット数を低
減でき、データの圧縮率を向上させることができるもの
である。By performing the compression processing in this way,
The number of bits for encoding can be reduced by the amount of data set to 0, and the data compression rate can be improved.
【0085】また、上述した実施例において、ウェーブ
レット変換として双直交ウェーブレット変換を用いた例
を示したが、これに限定されるものではなく、直交のウ
ェーブレット変換によって複数の周波数帯域に分割され
た係数信号を用いてもよい。Further, in the above-mentioned embodiment, the example in which the biorthogonal wavelet transform is used as the wavelet transform has been shown, but the present invention is not limited to this, and the coefficients divided into a plurality of frequency bands by the orthogonal wavelet transform are shown. A signal may be used.
【0086】さらに、上述した実施例においては所定の
周波数帯域のウェーブレット変換係数信号に基づいて非
連結成分を0とする処理を行うようにしているが、とく
にこれに限定されるものではなく、全ての周波数帯域の
データの非連結成分を0とする処理を行うようにしても
よいものである。すなわち、上述した本発明の実施例に
おいては、最も低い周波数帯域から最も高い周波数帯域
までの全ての周波数帯域のデータに対し、非連結成分を
0とする処理を行うようにしてもよいものである。Further, in the above-described embodiment, the processing for setting the non-connected component to 0 is performed based on the wavelet transform coefficient signal in the predetermined frequency band, but the processing is not particularly limited to this and all It is also possible to perform processing for setting the non-connected component of the data in the frequency band of 0 to 0. That is, in the above-described embodiment of the present invention, the processing of setting the non-connected component to 0 may be performed on the data in all the frequency bands from the lowest frequency band to the highest frequency band. .
【図1】本発明による画像データ圧縮処理方法の実施例
の基本的概念を表す図FIG. 1 is a diagram showing a basic concept of an embodiment of an image data compression processing method according to the present invention.
【図2】本発明に用いられる画像データの読み取り方式
を表す図FIG. 2 is a diagram showing a method of reading image data used in the present invention.
【図3】ウェーブレット変換の詳細を表す図FIG. 3 is a diagram showing details of a wavelet transform.
【図4】ウェーブレット変換係数信号を表す図FIG. 4 is a diagram showing a wavelet transform coefficient signal.
【図5】本発明による画像データ圧縮処理方法の実施例
を説明するための図FIG. 5 is a diagram for explaining an embodiment of an image data compression processing method according to the present invention.
【図6】逆ウェーブレット変換の詳細を表す図FIG. 6 is a diagram showing details of inverse wavelet transform.
【図7】ウェーブレット変換に用いられる基本ウェーブ
レット関数を表す図FIG. 7 is a diagram showing a basic wavelet function used for wavelet transform.
【図8】ウェーブレット変換を説明するための図FIG. 8 is a diagram for explaining a wavelet transform.
【図9】フーリエ変換を説明するための図FIG. 9 is a diagram for explaining a Fourier transform.
10 蓄積性蛍光体シート h,h′,g,g′ ウェーブレット変換を行うため
の関数 VVi ,VWi ,WVi ,WWi (i=1〜n) ウェーブレット変換係数信号10 stimulable phosphor sheet h, h ', g, g ' function VV i for performing a wavelet transform, VW i, WV i, WW i (i = 1~n) wavelet transform factor signals
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.6 識別記号 庁内整理番号 FI 技術表示箇所 H04N 7/30 ─────────────────────────────────────────────────── ─── Continuation of the front page (51) Int.Cl. 6 Identification code Office reference number FI technical display location H04N 7/30
Claims (2)
ータに圧縮処理を施す画像データ圧縮処理方法におい
て、 前記原画像データにウェーブレット変換を施すことによ
り、前記原画像データを複数の周波数帯域毎の画像デー
タに分解し、 該複数の周波数帯域毎の画像データのうち、所定の周波
数帯域の画像データにおいて該画像データの非連結成分
の値を0とする処理を行い、 該0とする処理が行われた全ての前記画像データおよび
前記0とする処理が行われなかった全ての前記画像デー
タを符号化することを特徴とする画像データ圧縮処理方
法。1. An image data compression processing method for performing compression processing on original image data representing an image including a predetermined subject, wherein the original image data is subjected to wavelet transformation to obtain the original image data for each of a plurality of frequency bands. Of the image data for each of the plurality of frequency bands, the process of setting the value of the non-connected component of the image data to 0 in the image data of the predetermined frequency band, and the process of setting the value to 0 are performed. An image data compression processing method, characterized in that all the image data that has been executed and all the image data that has not been processed to be 0 are encoded.
符号化し、 該逆符号化された前記各画像データに逆ウェーブレット
変換を施すことにより、請求項1記載の画像データ圧縮
処理方法により圧縮された前記原画像データを再構成す
ることを特徴とする画像データの再構成方法。2. The image data compression processing method according to claim 1, wherein the encoded image data is inversely encoded, and inversely wavelet transform is applied to the inversely encoded image data. A method for reconstructing image data, comprising reconstructing the original image data that has been created.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP14755493A JP3282136B2 (en) | 1993-06-18 | 1993-06-18 | Image data compression processing method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP14755493A JP3282136B2 (en) | 1993-06-18 | 1993-06-18 | Image data compression processing method |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH0723229A true JPH0723229A (en) | 1995-01-24 |
JP3282136B2 JP3282136B2 (en) | 2002-05-13 |
Family
ID=15432964
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP14755493A Expired - Fee Related JP3282136B2 (en) | 1993-06-18 | 1993-06-18 | Image data compression processing method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3282136B2 (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6173086B1 (en) | 1996-07-18 | 2001-01-09 | Fuji Photo Film Co., Ltd. | Image processing method |
US6801339B1 (en) | 1998-03-26 | 2004-10-05 | Fuji Photo Film Co., Ltd. | Image processing method and apparatus |
US6868187B2 (en) | 2000-05-12 | 2005-03-15 | Fuji Photo Film Co., Ltd. | Image transformation method and apparatus, and storage medium |
-
1993
- 1993-06-18 JP JP14755493A patent/JP3282136B2/en not_active Expired - Fee Related
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6173086B1 (en) | 1996-07-18 | 2001-01-09 | Fuji Photo Film Co., Ltd. | Image processing method |
US6801339B1 (en) | 1998-03-26 | 2004-10-05 | Fuji Photo Film Co., Ltd. | Image processing method and apparatus |
US6868187B2 (en) | 2000-05-12 | 2005-03-15 | Fuji Photo Film Co., Ltd. | Image transformation method and apparatus, and storage medium |
US7289678B2 (en) | 2000-05-12 | 2007-10-30 | Fujifilm Corporation | Image transformation method and apparatus, and storage medium |
Also Published As
Publication number | Publication date |
---|---|
JP3282136B2 (en) | 2002-05-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3499254B2 (en) | Image data compression processing method | |
JP4367880B2 (en) | Image processing apparatus and method, and storage medium | |
KR101303375B1 (en) | Method and device for encoding blocks of values | |
JPH0779405A (en) | Electronic still camera and picture reproducing device | |
KR100561587B1 (en) | Method and apparatus for three-dimensional wavelet transform | |
EP1705924A1 (en) | Moving picture encoding method and device, and moving picture decoding method and device | |
JPH08294119A (en) | Image coder/decoder | |
JP3282138B2 (en) | Image data compression processing method and image data reconstruction method | |
JPH04229382A (en) | Method and device for resolution conversion of digital image data | |
JP3282136B2 (en) | Image data compression processing method | |
JP3260008B2 (en) | Image data compression processing method | |
JPH09214967A (en) | Image data compression processing method | |
JP3282134B2 (en) | Image data compression processing method | |
JPH09307897A (en) | Image data compression processing method | |
JP3271098B2 (en) | Digital image signal decoding apparatus and method | |
JP3260016B2 (en) | Image data compression processing method and image data reconstruction method | |
JP3305480B2 (en) | Image encoding / decoding device | |
JP2980218B2 (en) | Image information encoding device and image information decoding device | |
US5900911A (en) | Signal encoding method separating a quantized signal into a binary signal and a sign and intensity signal | |
JPH09312845A (en) | Image data compression processing method | |
JPH10336643A (en) | Image coder and image decoder | |
JPH09307898A (en) | Image data compression processing method | |
JPH1188183A (en) | Wavelet converter, its method, wavelet inverse converter, its method, image coder, its method, image decoder and its method | |
JP3392949B2 (en) | Image compression device and image decompression device | |
JP4174254B2 (en) | Moving picture encoding apparatus, moving picture decoding apparatus, and methods thereof |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20020205 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080301 Year of fee payment: 6 |
|
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: 20080301 Year of fee payment: 6 |
|
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: 20080301 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090301 Year of fee payment: 7 |
|
LAPS | Cancellation because of no payment of annual fees |