WO2003041394A1 - Image compression apparatus, image compression program, and image compression method - Google Patents

Image compression apparatus, image compression program, and image compression method Download PDF

Info

Publication number
WO2003041394A1
WO2003041394A1 PCT/JP2002/011353 JP0211353W WO03041394A1 WO 2003041394 A1 WO2003041394 A1 WO 2003041394A1 JP 0211353 W JP0211353 W JP 0211353W WO 03041394 A1 WO03041394 A1 WO 03041394A1
Authority
WO
WIPO (PCT)
Prior art keywords
image compression
image
tile
compression
image data
Prior art date
Application number
PCT/JP2002/011353
Other languages
French (fr)
Japanese (ja)
Inventor
Hideyasu Kuniba
Original Assignee
Nikon Corporation
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 Nikon Corporation filed Critical Nikon Corporation
Priority to JP2003543305A priority Critical patent/JP4065978B2/en
Publication of WO2003041394A1 publication Critical patent/WO2003041394A1/en

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N19/00Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
    • H04N19/60Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using transform coding
    • H04N19/63Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using transform coding using sub-band based transform, e.g. wavelets
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N19/00Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
    • H04N19/10Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding
    • H04N19/102Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding characterised by the element, parameter or selection affected or controlled by the adaptive coding

Definitions

  • the present invention relates to an image compression apparatus, an image compression program, and an image compression method.
  • the present invention relates to an image compression device that divides image data into a plurality of tiles and performs image compression for each tile.
  • the present invention relates to an image compression program for causing a computer to function as an image compression device.
  • the present invention relates to an image compression method for compressing image data for each tile. Description of the background art
  • the screen size (the number of vertical and horizontal pixels) of the image data fluctuates, the screen size of each tile can be kept substantially constant by changing the number of tiles. In this case, there is an advantage that it is possible to flexibly cope with various screen sizes without significantly changing the image compression processing in units of tiles.
  • tile-based image compression is an important technology in image compression technology.
  • Patent Document 1 is known as a conventional technique for image compression in pixel block units, which are division ranges of orthogonal transform.
  • an activity indicating the degree to which a spatial frequency component is included is calculated for each pixel block, and the amount of compression code assigned to each pixel block is adjusted according to the activity.
  • Patent Document 1 Japanese Unexamined Patent Publication No. Hei 2 _ 2 6 2 786 (US Pat. No. 5,085,184)
  • Patent Document 1 described above, the larger the activity of a pixel block, the greater the amount of compression code to be allocated. Therefore, in the flat part where the tile boundary is most conspicuous (the activity is almost zero), the amount of compression code allocated is small. As a result, in Patent Document 1, it is difficult to positively improve the tile boundary.
  • Patent Document 1 only the compression code amount is adjusted in units of pixel blocks, and the frequency distribution of the compression code amount (particularly, the code amount allocated to low-frequency components) is not changed. Therefore, at very high image compression ratios where tile boundaries are a real problem, it seems difficult to make fine adjustments such as suppressing image information degradation in the mid-high frequency range while improving tile boundaries. Disclosure of the invention
  • an object of the present invention is to provide an image compression technique capable of appropriately suppressing tile boundaries from being inconspicuous at the time of image expansion.
  • An image compression apparatus is an image compression apparatus that divides image data into a plurality of tiles and performs image compression for each tile, and includes a statistical processing unit, a parameter setting unit, and an image compression unit.
  • the statistical processing unit processes the image data and displays a message “if a gradation shift occurs in the image data. And the statistical properties that correlate with the "conspicuousness of the gradation shift in the case.”
  • the parameter setting unit predicts the conspicuousness of the tile boundary at the time of image expansion based on this statistical property, and sets the “image compression parameter that affects the appearance of compression distortion” according to the conspicuousness of the tile boundary. Adjust the settings.
  • the image compression unit compresses the image data using the image compression parameters adjusted and set as described above.
  • the statistical property is a value indicating the flatness of the gradation of the image data.
  • the parameter setting unit adjusts and sets the image compression parameter in such a direction as to reduce the compression distortion of the low-frequency spatial frequency component as the degree of flatness increases.
  • the statistical property is a ratio of a gradation flat tile (tile that can be regarded as a flat gradation) to a plurality of tiles.
  • the parameter setting unit adjusts and sets the image compression parameter in a direction to reduce the compression distortion of the low-frequency spatial frequency component as the ratio of the gradation flat tile increases.
  • the statistical property is a variance obtained for a transform coefficient obtained by transforming image data into a spatial frequency domain.
  • the parameter setting unit adjusts and sets the image compression parameter in such a direction as to reduce the compression distortion of the low-frequency spatial frequency component as the variance of the transform coefficient is smaller.
  • the statistical property is a variance obtained for a transform coefficient obtained by performing a wavelet transform on the image data.
  • the parameter setting unit of the present invention uniformly sets the image compression parameters for all of the plurality of tiles.
  • the statistical processing unit of the present invention obtains a statistical property for each tile.
  • the parameter setting unit predicts the conspicuousness of the tile boundaries in tile units based on the statistical properties obtained in tile units, and reduces the compression distortion of low-frequency spatial frequency components as the tiles are more conspicuous. Adjust and set the image compression parameters.
  • the parameter setting unit also adjusts and sets the image compression parameter for the tile adjacent to the tile in which the tile boundary is conspicuous so as to reduce the compression distortion of the low-frequency spatial frequency component.
  • an image compression program according to the present invention causes a computer to function as the above-described statistical processing unit, parameter setting unit, and image compression unit.
  • the image compression method of the present invention is an image compression method in which image data is divided into a plurality of tiles, and image compression is performed for each tile.
  • a statistical processing step for obtaining a statistical property that correlates with "conspicuousness in the event of a shift"; and predicting the conspicuousness of the tile boundary at the time of image expansion based on the statistical property.
  • Parameter setting step to adjust and set the "image compression parameter that affects the amount of compression distortion" according to the conspicuousness of the image, and image compression to compress the image data using the adjusted image compression parameter And steps.
  • FIG. 1 is a diagram illustrating a configuration of an electronic camera 11 according to the present embodiment.
  • FIG. 2 is a flowchart illustrating an image compression process according to the first embodiment.
  • FIG. 3 is a diagram showing a state of subband division by wavelet transform.
  • FIG. 4 is a flowchart illustrating an image compression process according to the second embodiment.
  • FIG. 5 is a flowchart illustrating an image compression process according to the third embodiment.
  • Fig. 6 is a diagram showing the standard deviation (square root of variance) of the transform coefficients for each subband, picking up the characteristic tiles of the test image.
  • FIG. 7 is a diagram illustrating an example of the weighting factor Ws for each subband used for calculating the evaluation parameter k.
  • FIG. 8 is a diagram showing characteristic parameters k of the test image picked up.
  • FIG. 9 is a map in which the evaluation parameters k are arranged according to the tile division.
  • FIG. 1 is a diagram illustrating a configuration of an electronic camera 11 according to the present embodiment.
  • a photographing lens 12 is attached to an electronic camera 11.
  • the light receiving surface of the image sensor 13 is arranged in the image space of the photographing lens 12.
  • the output of the image sensor 13 is processed through the AD conversion unit 14 and the signal processing unit 15 and then recorded in the buffer memory 16.
  • An image processing unit 17 and an image compression unit 18 are connected to the data bus of the buffer memory 16 respectively.
  • the image compression section 18 has the following configuration requirements 1 to 5.
  • the compressed file generated in the image compression section 18 is provided to the recording section 32.
  • the recording unit 32 records and stores the compressed file on the memory card 33.
  • an operation member 38 including a release button, a setting button, and the like is provided on a housing of the electronic camera 11. Output signals of these operation members 38 are supplied to a microprocessor 39.
  • the microprocessor 39 is a microprocessor for controlling the system of the electronic camera 11 and is connected to the image compression section 18 and the image processing section 17 via signal lines.
  • the parameter setting unit described in the claims corresponds to the image compression unit 18 and the microprocessor 39.
  • the image compression unit described in the claims corresponds to the image compression unit 18.
  • the analog image data generated by the image sensor 13 is digitized via the AD converter 14 and then supplied to the signal processor 15.
  • the signal processing unit 15 performs signal processing such as black level correction and gamma correction on the image data.
  • the image data after the signal processing is temporarily recorded in the buffer memory 16.
  • the image processing unit 17 reads out the image data in the buffer memory 16 in units of processing and performs image processing such as color interpolation processing and color coordinate conversion.
  • the image compression section 18 reads out the image data after the image processing from the buffer memory 16 and performs the image compression processing.
  • FIG. 2 is a flowchart illustrating this image compression processing.
  • Step S 1 The tile processing section 21 in the image compression section 18 divides the image data in the buffer memory 16 into tile units and sequentially reads the data.
  • Step S2 The wavelet transform unit 22 performs wavelet transform on each tile and divides the tile into a plurality of subbands as shown in FIG.
  • Step S3 The wavelet transformer 22 finds the standard deviation (square root of variance) of the wavelet transform coefficient for each subband for each tile. In general, the smaller these standard deviations are overall, the higher the flatness of the gradation (including the color tone) of the image data. Therefore, even if a slight gradation shift occurs in the image data, it can be determined that the gradation shift is conspicuous.
  • the wavelet transform unit 22 transmits the standard deviation thus obtained to the microprocessor 39.
  • Step S4 The microprocessor 39 selects an image according to the standard deviation of the conversion coefficient. Estimate the conspicuousness of discontinuous compression distortion generated at tile boundaries during image expansion. The microprocessor 39 adjusts and sets the visual weight (visual importance) of each subband according to the prediction situation.
  • the microprocessor 39 obtains the ratio of the gradation flat tile whose standard deviation is equal to or smaller than the threshold value (here, equal to or smaller than 3) in the 3 LL band (lowest band) of each tile, and determines the ratio as r 1. I do. As the value of 1 is larger, the gradation of the image data is globally flat, and it can be determined that the gradation deviation is more conspicuous.
  • the microprocessor 39 calculates the ratio of the flat tone tiles whose standard deviation is equal to or less than the threshold value (here, 3 or less) in the 3 LH, 3HL, and 3 HH bands of each tile, and calculates the ratio as r 2 And This is the ratio of tiles that can be regarded as gradation flat in the intermediate spatial frequency range. In addition to rl, it can be determined that the larger the value of r2, the flatter the gradation of the image data and the more noticeable the gradation deviation. Further, the microprocessor 39 obtains the ratio of the gradation flat tile whose standard deviation is equal to or smaller than the threshold value (here, equal to or smaller than 3) among the 2 LH, 2 HL, and 2HH bands of each tile. I do.
  • the threshold value here, 3 or less
  • This r 3 is also the ratio of tiles that can be regarded as having a flat gradation in the intermediate spatial frequency range. In addition to r 1 and r 2, it can be determined that the larger the value of 3 is, the flatter the gradation of the image data is, and the more the gradation deviation is more conspicuous.
  • the microprocessor 39 calculates the ratio of the gradation flat tile whose standard deviation is equal to or less than the threshold value (here, 3 or less) in the 1 LH, 1 HL, and 1 HH bands of each tile, and calculates the ratio as r And 4.
  • This r4 is the ratio of tiles that can be regarded as gradation flat in the high spatial frequency range. In addition to r l to r 3, it can be determined that as the value of r 4 is larger, the gradation of the image data is flatter like silk, and the gradation deviation is more conspicuous.
  • This parameter k is the same as r1 to r4 described above. These are integrated evaluation parameters. Therefore, it can be determined that the larger the evaluation parameter k is, the flatter the gradation of the image data is from the low frequency band to the high frequency band, and that the gradation deviation is conspicuous.
  • image data having a large evaluation parameter k for example, a picture such as a blue sky in fine weather is assumed.
  • image data is expected to be particularly conspicuous because the tile boundaries are raised because of the monotonous pattern.
  • image data with a small evaluation parameter k is assumed to be a complex picture rich in detail.
  • discontinuous compression distortions generated at tile boundaries are buried in complex patterns and are expected to be inconspicuous.
  • the microprocessor 39 determines the visual weight of each subband by using the following equation in consideration of the prediction situation based on the evaluation parameter k.
  • the microprocessor 39 transmits the selected visual weight to the quantization unit 23.
  • Step S5 The quantization unit 23 optimizes the quantization step width of the transform coefficient so as to minimize visual image quality deterioration due to quantization distortion under an appropriate coding rate constraint condition.
  • the method of optimizing the quantization step width will be described in principle.
  • the transform coefficients for each subband show a probability distribution very similar to the Laplace distribution. Therefore, the probability distribution of the conversion coefficient is approximated by the following Laplace distribution equation, and the variable ⁇ that specifies the Laplace distribution is obtained.
  • f (y) ⁇ -e- ⁇ ⁇ ' ⁇ [1]
  • the quantized transform coefficient takes the k-th quantized value.
  • the probability pk is Becomes
  • the coding rate R (Q) after quantization is ideally equal to the entropy, and
  • the quantization step width of the i-th subband is Q i
  • the coding rate at that time is R i (Q i)
  • the quantization distortion is D i (R i).
  • the quantization step width Q i is determined.
  • the quantization unit 23 refers to this data table based on the variable ⁇ and the visual weight, The optimal quantization step width can be determined immediately.
  • Step S 6 The quantization unit 23 quantizes the transform coefficient of each subband according to the quantization step width Q i thus optimized.
  • the quantization step width of the lower sub-band is set to be narrower as the evaluation parameter k is larger.
  • the lower the compression distortion (quantization distortion in this case) the lower the image data with a large evaluation parameter k.
  • Step S7 The bit modeling section 24 reads a mask image of the selected area preset in the R ⁇ I setting section 27.
  • the bit modeling section 24 It is determined whether or not the transform coefficient is included in the selected area for each subband based on the mask image of.
  • the bit modeling unit 24 shifts the quantized transform coefficient located in the selected area by S bits.
  • the S bits are the number of bits determined according to the known Max Shift method.
  • Step S8 The bit modeling unit 24 divides the quantized transform coefficient into bit planes.
  • the bit modeling unit 24 executes three kinds of encoding passes specified in JPEG 2000 in order from the bit plane of the most significant digit.
  • Step S 9 The arithmetic encoding unit 25 acquires encoded data from the bit modeling unit 24 for each processing unit. Arithmetic coding section 25 performs arithmetic coding on the coded data by using an MQ coder 1 which is a binary arithmetic code.
  • Step S10 The bitstream generation unit 26 rearranges the encoded data according to the SNR progressive and other priorities, and generates a bitstream.
  • Step S11 The bitstream generation unit 26 cuts off this bitstream with a predetermined code amount and outputs it as a compressed file.
  • the recording unit 32 records and stores the compressed file on the memory card 33.
  • the electronic camera 11 predicts the conspicuousness of the tile boundary using the standard deviation of the wavelet transform coefficient.
  • the quantization step width is adjusted and set so that the image distortion of the tile boundary is less conspicuous and the compression distortion in the low frequency band is suppressed lower.
  • FIG. 4 is a diagram illustrating an image compression process according to the second embodiment. Hereinafter, the operation of the second embodiment will be described along the step numbers shown in FIG.
  • Step S21 The microprocessor 39 extracts the image data to be compressed by dividing it into a plurality of divided areas (for example, tiles), and standard deviation of pixel values (values of luminance, color, color difference, etc.) for each divided area. Ask for.
  • Step S22 The ratio of the "segmented region where the standard deviation is equal to or less than a predetermined threshold" to the entire extracted divided region is determined, and the ratio is defined as the flatness R of the image data. It can be determined that the larger the flatness R is, the higher the flatness of the gradation (including the color tone) of the image data is and the more noticeable the gradation deviation is.
  • Step S23 It is predicted that the larger the flatness R is, the more noticeable the discontinuous compression distortion that occurs at the tile boundary during image expansion. Therefore, the microprocessor 39 relatively increases the visual weight of the low-frequency subband as the value of the flatness R increases. The microprocessor 39 transmits the visual weight adjusted and set as described above to the quantization unit 23.
  • Step S24 The tile processing section 21 divides the image data in the buffer memory 16 into tile units and sequentially reads them.
  • Step S25 The wavelet transform unit 22 independently performs wavelet transform on each tile and divides it into subbands as shown in FIG.
  • steps S26 to S32 are the same as steps S5 to S11 described in the first embodiment, and a description thereof will be omitted.
  • the conspicuousness of the tile boundary is predicted based on the standard deviation of the pixel values.
  • the prediction operation can be performed without waiting for the wavelet change.
  • the processing of predicting the conspicuousness of the tile boundaries and the wavelet transform can be processed in parallel. Suitable for improving efficiency.
  • FIG. 5 is a diagram illustrating an image compression process according to the third embodiment. Hereinafter, the operation of the third embodiment will be described along the step numbers shown in FIG.
  • Steps S41 to S42 Same as steps S1 to S2 described in the first embodiment.
  • the variance is calculated for the entire area inside the tile.
  • the variance may be calculated only for the inner circumference of the tile related to the tile boundary (that is, the central area of the tile not related to the tile boundary may be excluded from the variance calculation).
  • the variance may be calculated by dividing the tile into the vicinity of the top side, near the bottom side, near the left side, and near the right side. .
  • FIG. 6 is a diagram in which characteristic image regions of a test image (person image) are picked up and the standard deviation of the ⁇ ⁇ -Apelet transform coefficient is shown for each sub-band.
  • the tile corresponding to the “background” has a uniformly low standard deviation from the high-frequency subband to the low-frequency subband because the entire image is blurred, and the tile boundary is It can be judged that it is very noticeable.
  • the tiles of “Knitted Clothes”, “Eyes”, and “Mouths” It can be determined that the standard deviation is relatively high over the area subband, and the tile boundaries are relatively inconspicuous.
  • the “hair” tile has a high standard deviation value from the low-frequency subband to the high-frequency subband. Therefore, it can be determined that the tile boundaries are hardly noticeable because they are buried in a complicated pattern.
  • the wavelet transform unit 22 transmits the standard deviation thus obtained to the microprocessor 39.
  • Step S44 The microprocessor 39 obtains an evaluation parameter k which indicates the conspicuousness of the tile boundary by a numerical value in order to perform the evaluation judgment as shown in FIG. 6 by calculation.
  • the evaluation parameter k here is a value calculated for each tile, using the following equation, different from the first embodiment. ...]
  • the weight coefficient W s in the expression [10] is a value indicating the ratio of each subband S contributing to the conspicuousness of the tile boundary, and may be determined based on a subjective evaluation experiment of an image or the like.
  • FIG. 7 is a table showing an example of the weight coefficient Ws determined in this way.
  • a in [10] is a proportionality constant. '
  • the evaluation parameter k is determined using the [10] expression, but the present invention is not limited to this evaluation expression.
  • An evaluation formula that calculates an evaluation parameter that makes the tile boundary less noticeable as the value of the variance is larger may be used.
  • Fig. 8 is a table showing the values of the evaluation parameter k by picking up characteristic image areas of the test image (person image). It can be predicted that the larger the value of the evaluation parameter k is, the more noticeable the tile boundary is.
  • Step S45 Subsequently, the microprocessor 39 adjusts the value of the value parameter k between adjacent tiles.
  • a local area including the upper, lower, left, and right neighboring tiles is first set for each tile to be processed by the microprocessor.
  • the microprocessor 39 smoothes the evaluation parameter k (such as a simple average and a weighted average) in the local region.
  • Figure 9 [A] is a map in which the evaluation parameters k of the human image are arranged according to the tile division.
  • Fig. 9 [B] is a map showing the evaluation parameter k after the smoothing process.
  • a calculation operation such as a MAX value or a median value may be used.
  • the evaluation parameter k may be adjusted between actually adjacent sides.
  • the adjustment of the evaluation parameter k may be performed between the left side of the tile located on the right and the right side of the tile located on the left.
  • the visual weight T 3LL of the low-frequency subband is set to be relatively large.
  • the microprocessor 39 transmits the visual weight thus selected to the quantization unit 23.
  • Step S47 The quantization unit 23 performs a quantization step of the transform coefficient according to the visual weight set for each tile so as to minimize the quantization distortion and the visual quality degradation. Is optimized for each tile.
  • the procedure of this optimization method is the same as that in step S5 of the first embodiment, and thus the duplicate description is omitted here.
  • Step S48 The quantizing section 23 quantizes the transform coefficient for each tile using the quantization step width of each subband optimized for each tile.
  • the quantization step width of the lower sub-band is set to be narrower as the evaluation parameter k is larger.
  • tiles with a larger evaluation parameter k have lower compression distortion (quantization distortion in this case) in the low frequency band.
  • Steps S49 to S53 Same as steps S7 to S11 described in the first embodiment.
  • Step S54 The microprocessor 39 adds data relating to the quantization step width determined for each tile to the compressed data to generate a compressed file in preparation for using data at the time of image expansion.
  • the recording unit 32 records and stores the generated compressed file on the memory card 33.
  • tiles in which tile boundaries are conspicuous are appropriately selected and low-frequency compression distortion (that is, occurrence of tile boundaries) which is a problem is suppressed in tile units.
  • low-frequency compression distortion that is, occurrence of tile boundaries
  • tile units for example, tiles such as blue sky
  • it is possible to appropriately select tiles whose tile boundaries are not visually conspicuous for example, tiles with complicated patterns
  • moderately low-frequency compression distortion for tiles with complicated patterns, the image information of the mid- and high-frequency components is compressed as accurately as possible to the extent that tile boundaries are allowed to occur, and image information deterioration due to compression is prevented as much as possible. It becomes possible.
  • the frequency distribution of the compression code amount is adjusted for each tile based on the conspicuousness of the tile boundary.
  • tile boundaries are a real problem
  • low-frequency compression distortion is suppressed also around tiles with noticeable tile boundaries.
  • the conversion coefficient and the standard deviation (square root of the variance) of the pixel value are obtained to evaluate the conspicuousness of the gradation shift.
  • the present invention is not limited to this.
  • any statistical property of the present invention may be used as long as the quantity correlates with “conspicuousness of gradation deviation”.
  • statistical properties indicating gradation flatness, gradation uniformity, gradation monotony, and the like may be used.
  • variable of the present invention may be any measure that indicates the degree of dispersion of values, and is not limited to the statistical calculation formula expressed by the formula [9].
  • the dynamic range of the pixel value is a value that fluctuates greatly even if there is a large gradation change even momentarily. Therefore, the dynamic range of the pixel value is not a value indicating the flatness of the gradation, and is not an amount correlating with “conspicuousness of gradation deviation”.
  • the quantization step size of the transform coefficient is adjusted and set as the image compression parameter.
  • the present invention is not limited to this.
  • any image compression parameter of the present invention may be used as long as it affects the amount of appearance of compression distortion at a tile boundary.
  • the image compression parameter may be a parameter for adjusting the lower digit truncation of the bit plane, a parameter for adjusting the bit stream truncation, or a parameter for adjusting a selected region of ROI coding.
  • the image compression parameters are uniformly changed for all the plurality of tiles. In this case, since the image compression is uniformly applied to the entire plurality of tiles, there is little possibility that the compression distortion will appear on the screen unevenly.
  • the present invention is not limited to this.
  • the third embodiment As shown in the figure, the image compression parameters may be adjusted and set for each individual tile. Further, the image compression parameters may be individually adjusted and set at the tile boundary and inside the tile.
  • the image compression device of the present invention may be realized on a computer.
  • an image compression program (corresponding to claim 9) may be created by converting the image compression procedure in the above-described embodiment into a program code.
  • the present invention may be embodied in various other forms without departing from the spirit or main characteristics thereof. Therefore, the above-described embodiment is merely an example in every aspect, and should not be construed as limiting.
  • the scope of the present invention is defined by the scope of the claims, and is not limited by the text of the specification. Further, all modifications and changes belonging to the equivalent scope of the claims are within the scope of the present invention.
  • a statistical property correlating with “visual conspicuousness of gradation shift” is obtained from image data, and in accordance with the conspicuous prominence of a tile boundary predicted by the statistical property, Adjust and set image compression parameters that affect the appearance of compression distortion.

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Compression Or Coding Systems Of Tv Signals (AREA)
  • Compression Of Band Width Or Redundancy In Fax (AREA)

Abstract

An image compression apparatus for dividing image data into a plurality of tiles and performing image compression for each of the tiles. The apparatus includes a statistical processing unit, a parameter setting unit, and an image compression unit. The statistical processing unit processes image data to obtain a statistical feature correlating with “noticeability of a gradation shift when a gradation shift is caused in the image data”. According to this statistical feature, the parameter setting unit predicts noticeability of the tile boundary upon image decompression and adjusts “an image compression parameter affecting the amount of compression distortion”. The image compression unit performs image compression of image data by using the image compression parameter thus adjusted.

Description

明 細 書  Specification
画像圧縮装置、 画像圧縮プログラム、 および画像圧縮方法 発明の分野  FIELD OF THE INVENTION The present invention relates to an image compression apparatus, an image compression program, and an image compression method.
本発明は、 画像データを複数のタイルに分割し、 各タイルごとに画像圧縮する 画像圧縮装置に関する。  The present invention relates to an image compression device that divides image data into a plurality of tiles and performs image compression for each tile.
本発明は、 コンピュータを画像圧縮装置として機能させるための画像圧縮プロ グラムに関する。  The present invention relates to an image compression program for causing a computer to function as an image compression device.
本発明は、 画像データを各タイルごとに画像圧縮する画像圧縮方法に関する。 背景技術の説明  The present invention relates to an image compression method for compressing image data for each tile. Description of the background art
従来、 J P E G 2 0 0 0のように、 画像データをタイルと呼ばれる小領域に区 分し、 これらタイルごとに画像圧縮を実施する画像圧縮方式が知られている。 このようなタイル単位の画像圧縮では、 大量の画像データを一度に扱う必要が なく、 ワークメモリを節約できるという利点がある。  Conventionally, there has been known an image compression method in which image data is divided into small areas called tiles and image compression is performed for each of the tiles, as in JPEG200. This type of tile-based image compression has the advantage of saving work memory without having to deal with a large amount of image data at once.
また、 画像データの画面サイズ (縦横画素数) が変動しても、 タイルの分割数 を変化させることにより、各タイルの画面サイズを概ね一定に保つことができる。 この場合、 タイル単位の画像圧縮処理をさほど変えずに、 多様な画面サイズに柔 軟に対応できるという利点がある。  Even if the screen size (the number of vertical and horizontal pixels) of the image data fluctuates, the screen size of each tile can be kept substantially constant by changing the number of tiles. In this case, there is an advantage that it is possible to flexibly cope with various screen sizes without significantly changing the image compression processing in units of tiles.
さらに、 タイル単位に画像圧縮処理を分離実行できるので、 画像圧縮の並列処 理ゃマルチスレツド化などを容易に実現できるといった利点もある。  Furthermore, since image compression processing can be separated and executed for each tile, there is an advantage that parallel processing of image compression and multi-threading can be easily realized.
以上のような利点により、 タイル単位の画像圧縮は、 画像圧縮技術において重 要な技術となっている。  Due to the above advantages, tile-based image compression is an important technology in image compression technology.
なお、 タイル単位の画像圧縮とは異なるが、 直交変換の区分範囲である画素ブ ロック単位に画像圧縮する従来技術として、 下記の特許文献 1が知られている。 この特許文献 1では、 空間周波数成分が含まれる程度を示すァクティビティを画 素ブロックごとに算出し、 そのァクティビティに応じて各画素プロックに割り当 てる圧縮符号量を調整している。  It should be noted that, although different from image compression in tile units, Patent Document 1 below is known as a conventional technique for image compression in pixel block units, which are division ranges of orthogonal transform. In Patent Document 1, an activity indicating the degree to which a spatial frequency component is included is calculated for each pixel block, and the amount of compression code assigned to each pixel block is adjusted according to the activity.
(特許文献 1 ) 特開平 2 _ 2 6 2 7 8 6号公報 (U S 5 0 5 1 8 4 0 ) (Patent Document 1) Japanese Unexamined Patent Publication No. Hei 2 _ 2 6 2 786 (US Pat. No. 5,085,184)
ところで、 上述したタイル単位の画像圧縮では、 隣接するタイルにおいて個別 に画像圧縮が行われる。 そのため、 圧縮歪みもタイル毎に独立に発生し、 タイル 境界では圧縮歪みが不連続となる。 その結果、 伸張後の画像には、 タイル間の境 界が境目 (以下、 この境目を 『タイル境界』 という) となって見えるなどの問題 があった。  Meanwhile, in the above-described image compression on a tile-by-tile basis, image compression is individually performed on adjacent tiles. Therefore, the compression distortion also occurs independently for each tile, and the compression distortion is discontinuous at the tile boundary. As a result, the decompressed image had problems such as boundaries between tiles appearing as boundaries (hereinafter, these boundaries are referred to as “tile boundaries”).
このようなタイル境界の問題を改善するには、 可逆圧縮に近い画像圧縮を行つ て圧縮歪みの発生自体を抑制すればよい。 しかしながら、 必要以上に圧縮歪みを 抑制すれば、 画像圧縮の符号化効率が低くなり、 圧縮データ量が徒らに増えるな どの問題が生じてしまう。  In order to improve such a problem of the tile boundary, it is only necessary to perform image compression close to lossless compression to suppress the generation of compression distortion itself. However, if the compression distortion is suppressed more than necessary, the encoding efficiency of image compression will be reduced and the amount of compressed data will increase unnecessarily.
また、 上述した特許文献 1では、 アクティビティの大きな画素ブロックほど、 割り当てる圧縮符号量を増加させる。 そのため、 タイル境界が最も目立つ平坦部 (アクティビティは殆どゼロ) では、 割り当てる圧縮符号量が少なくなる。 その 結果、 特許文献 1では、 タイル境界を積極的に改善することは難しい。  Further, in Patent Document 1 described above, the larger the activity of a pixel block, the greater the amount of compression code to be allocated. Therefore, in the flat part where the tile boundary is most conspicuous (the activity is almost zero), the amount of compression code allocated is small. As a result, in Patent Document 1, it is difficult to positively improve the tile boundary.
さらに、特許文献 1では、画素プロックの単位に圧縮符号量を調整するのみで、 圧縮符号量の周波数配分 (特に、 低周波成分に割り当てる符号量) については変 更されない。 そのため、 タイル境界が現実問題となる極めて高い画像圧縮率にお いて、 タイル境界を改善しつつ、 中高域周波数の画像情報劣化を抑えるといった 微妙な調整は困難であつたと思われる。 発明の開示  Further, in Patent Document 1, only the compression code amount is adjusted in units of pixel blocks, and the frequency distribution of the compression code amount (particularly, the code amount allocated to low-frequency components) is not changed. Therefore, at very high image compression ratios where tile boundaries are a real problem, it seems difficult to make fine adjustments such as suppressing image information degradation in the mid-high frequency range while improving tile boundaries. Disclosure of the invention
そこで、 本発明では、 上述した問題点に鑑みて、 画像伸張時にタイル境界が目 立たないよう適度に抑制することが可能な画像圧縮技術を提供することを目的と する。  In view of the above-mentioned problems, an object of the present invention is to provide an image compression technique capable of appropriately suppressing tile boundaries from being inconspicuous at the time of image expansion.
以下、 本発明について説明する。  Hereinafter, the present invention will be described.
本発明の画像圧縮装置は、 画像データを複数のタイルに区分し、 各タイルごと に画像圧縮を施す画像圧縮装置であって、 統計処理部、 パラメータ設定部、 およ ぴ画像圧縮部を備える。  An image compression apparatus according to the present invention is an image compression apparatus that divides image data into a plurality of tiles and performs image compression for each tile, and includes a statistical processing unit, a parameter setting unit, and an image compression unit.
この統計処理部は、 画像データを処理して 「画像データに階調ズレを生じた場 合の階調ズレの目立ちやすさ」 と相関する統計的性質を求める。 The statistical processing unit processes the image data and displays a message “if a gradation shift occurs in the image data. And the statistical properties that correlate with the "conspicuousness of the gradation shift in the case."
パラメータ設定部は、 この統計的性質に基づいて、 画像伸張時におけるタイル 境界の目立ちやすさを予測し、 このタイル境界の目立ちやすさに応じて 「圧縮歪 みの出現量を左右する画像圧縮パラメータ」 を調整設定する。  The parameter setting unit predicts the conspicuousness of the tile boundary at the time of image expansion based on this statistical property, and sets the “image compression parameter that affects the appearance of compression distortion” according to the conspicuousness of the tile boundary. Adjust the settings.
画像圧縮部は、 このように調整設定された画像圧縮パラメータを用いて、 画像 データを画像圧縮する。  The image compression unit compresses the image data using the image compression parameters adjusted and set as described above.
さらに好ましくは、 本発明において、 統計的性質を、 画像データの階調の平坦 さを示す値にする。. パラメータ設定部は、 この平坦さの度合いが高いほど、 低域 空間周波数成分の圧縮歪みを軽減する方向に画像圧縮パラメータを調整設定する。 また好ましくは、 本発明において、 統計的性質を、 複数のタイルに占める階調 平坦タイル(平坦な階調とみなせるタイル) の割合とする。 パラメータ設定部は、 この階調平坦タイルの割合が多いほど、 低域空間周波数成分の圧縮歪みを軽減す る方向に画像圧縮パラメータを調整設定する。  More preferably, in the present invention, the statistical property is a value indicating the flatness of the gradation of the image data. The parameter setting unit adjusts and sets the image compression parameter in such a direction as to reduce the compression distortion of the low-frequency spatial frequency component as the degree of flatness increases. Also preferably, in the present invention, the statistical property is a ratio of a gradation flat tile (tile that can be regarded as a flat gradation) to a plurality of tiles. The parameter setting unit adjusts and sets the image compression parameter in a direction to reduce the compression distortion of the low-frequency spatial frequency component as the ratio of the gradation flat tile increases.
さらに好ましくは、 本発明において、 統計的性質を、 画像データを空間周波数 領域に変換した変換係数について求めた分散とする。 パラメータ設定部は、 この 変換係数の分散が小さいほど、 低域空間周波数成分の圧縮歪みを軽減する方向に 画像圧縮パラメータを調整設定する。  More preferably, in the present invention, the statistical property is a variance obtained for a transform coefficient obtained by transforming image data into a spatial frequency domain. The parameter setting unit adjusts and sets the image compression parameter in such a direction as to reduce the compression distortion of the low-frequency spatial frequency component as the variance of the transform coefficient is smaller.
また好ましくは、 本発明において、 統計的性質を、 画像データをウエーブレツ ト変換して得た変換係数について求めた分散とする。  Also preferably, in the present invention, the statistical property is a variance obtained for a transform coefficient obtained by performing a wavelet transform on the image data.
さらに好ましくは、本発明のパラメータ設定部が、複数のタイル全体について、 画像圧縮パラメータを一律に設定する。  More preferably, the parameter setting unit of the present invention uniformly sets the image compression parameters for all of the plurality of tiles.
また好ましくは、 本発明の統計処理部が、 タイル単位に統計的性質を求める。 パラメータ設定部は、 タイル単位に求めた統計的性質に基づいて、 タイル境界の 目立ちやすさをタイル単位に予測し、 タイル境界の目立ちやすいタイルほど、 低 域空間周波数成分の圧縮歪みを軽減する方向に画像圧縮パラメータを調整設定す る。  Also preferably, the statistical processing unit of the present invention obtains a statistical property for each tile. The parameter setting unit predicts the conspicuousness of the tile boundaries in tile units based on the statistical properties obtained in tile units, and reduces the compression distortion of low-frequency spatial frequency components as the tiles are more conspicuous. Adjust and set the image compression parameters.
さらに、 好ましくは、 パラメータ設定部が、 タイル境界が目立ちやすいタイル に隣接するタイルについても、 低域空間周波数成分の圧縮歪みを軽減する方向に 画像圧縮パラメータを調整設定する。 一方、 本発明の画像圧縮プログラムは、 コンピュータを、 上述した統計処理部、 パラメータ設定部、 および画像圧縮部として機能させることを特徴とする。 Further, preferably, the parameter setting unit also adjusts and sets the image compression parameter for the tile adjacent to the tile in which the tile boundary is conspicuous so as to reduce the compression distortion of the low-frequency spatial frequency component. On the other hand, an image compression program according to the present invention causes a computer to function as the above-described statistical processing unit, parameter setting unit, and image compression unit.
また、 本発明の画像圧縮方法は、 画像データを複数のタイルに区分し、 各タイ ルごとに画像圧縮を施す画像圧縮方法であって、 画像データを処理して 「画像デ 一タに階調ズレを生じた場合の目立ちやすさ」 と相関する統計的性質を求める統 計処理ステップと、 この統計的性質に基づいて画像伸張時におけるタイル境界の 目立ちやすさを予測し、 予測されたタイル境界の目立ちやすさに応じて 「圧縮歪 みの出現量を左右する画像圧縮パラメータ」 を調整設定するパラメータ設定ステ ップと、 調整設定された画像圧縮パラメータを用いて画像データを画像圧縮する 画像圧縮ステップとを有する。 図面の簡単な説明  Further, the image compression method of the present invention is an image compression method in which image data is divided into a plurality of tiles, and image compression is performed for each tile. A statistical processing step for obtaining a statistical property that correlates with "conspicuousness in the event of a shift"; and predicting the conspicuousness of the tile boundary at the time of image expansion based on the statistical property. Parameter setting step to adjust and set the "image compression parameter that affects the amount of compression distortion" according to the conspicuousness of the image, and image compression to compress the image data using the adjusted image compression parameter And steps. BRIEF DESCRIPTION OF THE FIGURES
なお、 本発明における上述した目的およびそれ以外の目的は、 以下の説明と添 付図面とによつて容易に確認することができる。  The above objects and other objects of the present invention can be easily confirmed by the following description and the accompanying drawings.
図 1は、 本実施形態における電子カメラ 1 1の構成を示す図である。  FIG. 1 is a diagram illustrating a configuration of an electronic camera 11 according to the present embodiment.
図 2は、 第 1の実施形態における画像圧縮処理を説明する流れ図である。  FIG. 2 is a flowchart illustrating an image compression process according to the first embodiment.
図 3は、 ウェーブレツト変換によるサブバンド分割の様子を示す図である。 図 4は、 第 2の実施形態における画像圧縮処理を説明する流れ図である。  FIG. 3 is a diagram showing a state of subband division by wavelet transform. FIG. 4 is a flowchart illustrating an image compression process according to the second embodiment.
図 5は、 第 3の実施形態における画像圧縮処理を説明する流れ図である。  FIG. 5 is a flowchart illustrating an image compression process according to the third embodiment.
図 6は、 テス ト画像の特徴的なタイルをピックアップして、 変換係数の標準偏 差 (分散の平方根) をサブバンド別に示した図である。  Fig. 6 is a diagram showing the standard deviation (square root of variance) of the transform coefficients for each subband, picking up the characteristic tiles of the test image.
図 7は、 評価パラメータ kの計算に使用するサブバンド別の重み係数 W sの一 例を示す図である。  FIG. 7 is a diagram illustrating an example of the weighting factor Ws for each subband used for calculating the evaluation parameter k.
図 8は、 テスト画像の特徴的なタイルをピックアップして、 評価パラメータ k を示した図である。  FIG. 8 is a diagram showing characteristic parameters k of the test image picked up.
図 9は、 評価パラメータ kをタイル区分に従って配列したマップである。 発明を実施するための最良の形態  FIG. 9 is a map in which the evaluation parameters k are arranged according to the tile division. BEST MODE FOR CARRYING OUT THE INVENTION
以下、 図面に基づいて本発明にかかる実施形態を説明する。 《第 1の実施形態》 Hereinafter, embodiments of the present invention will be described with reference to the drawings. << 1st Embodiment >>
[電子カメラの構成説明]  [Description of electronic camera configuration]
図 1は、 本実施形態における電子カメラ 1 1の構成を示す図である。  FIG. 1 is a diagram illustrating a configuration of an electronic camera 11 according to the present embodiment.
図 1において、 電子カメラ 1 1には、 撮影レンズ 1 2が装着される。 この撮影 レンズ 1 2の像空間には、 撮像素子 1 3の受光面が配置される。 この撮像素子 1 3の出力は、 A D変換部 1 4および信号処理部 1 5を介して処理された後、 ノ ッ ファメモリ 1 6に記録される。  In FIG. 1, a photographing lens 12 is attached to an electronic camera 11. The light receiving surface of the image sensor 13 is arranged in the image space of the photographing lens 12. The output of the image sensor 13 is processed through the AD conversion unit 14 and the signal processing unit 15 and then recorded in the buffer memory 16.
このバッファメモリ 1 6のデータバスには、 画像処理部 1 7および画像圧縮部 1 8がそれぞれ接続される。  An image processing unit 17 and an image compression unit 18 are connected to the data bus of the buffer memory 16 respectively.
この画像圧縮部 1 8は、 下記のような構成要件①〜⑦を備える。  The image compression section 18 has the following configuration requirements 1 to 5.
①タイル処理部 2 1  ①Tile processing unit 2 1
②ウエーブレット変換部 2 2  ②Wavelet converter 2 2
③量子化部 2 3  ③ Quantization unit 2 3
④ビットモデリング部 2 4  ④Bit modeling part 2 4
⑤算術符号化部 2 5 ⑤ Arithmetic encoder 2 5
⑥ビットストリーム生成部 2 6  ⑥ Bit stream generator 2 6
⑦ R O I設定部 2 7  ⑦ ROI setting section 2 7
このような画像圧縮部 1 8において生成された圧縮ファイルは、 記録部 3 2に 与えられる。 記録部 3 2は、 この圧縮ファイルを、 メモリカード 3 3に記録保存 する。  The compressed file generated in the image compression section 18 is provided to the recording section 32. The recording unit 32 records and stores the compressed file on the memory card 33.
さらに、 電子カメラ 1 1の筐体には、 レリーズ釦ゃ設定釦などからなる操作部 材 3 8が設けられる。 これらの操作部材 3 8の出力信号は、 マイクロプロセッサ 3 9に与えられる。 このマイクロプロセッサ 3 9は、 電子カメラ 1 1をシステム 制御するためのマイクロプロセッサであり、 画像圧縮部 1 8および画像処理部 1 7などと信号線を介して接続される。  Further, an operation member 38 including a release button, a setting button, and the like is provided on a housing of the electronic camera 11. Output signals of these operation members 38 are supplied to a microprocessor 39. The microprocessor 39 is a microprocessor for controlling the system of the electronic camera 11 and is connected to the image compression section 18 and the image processing section 17 via signal lines.
[発明との対応関係]  [Correspondence with invention]
ここで、 請求項の記載事項と、 上述した電子カメ.ラ 1 1との対応関係について 説明する。 なお、 ここでの対応関係は、 参考のために一解釈を例示するものであ り、 本発明を徒らに限定するものではない。 請求項記載の統計処理部は、 画像圧縮部 1 8およびマイクロプロセッサ 3 9に 対応する。 Here, the correspondence between the items described in the claims and the above-described electronic camera 11 will be described. It should be noted that the correspondence relationship here is merely an example for reference, and does not limit the present invention. The statistical processing unit described in the claims corresponds to the image compression unit 18 and the microprocessor 39.
請求項記載のパラメータ設定部は、 画像圧縮部 1 8およびマイクロプロセッサ 3 9に対応する。  The parameter setting unit described in the claims corresponds to the image compression unit 18 and the microprocessor 39.
請求項記載の画像圧縮部は、 画像圧縮部 1 8に対応する。  The image compression unit described in the claims corresponds to the image compression unit 18.
[第 1の実施形態の動作説明]  [Explanation of Operation of First Embodiment]
以下、 第 1の実施形態の動作について説明する。  Hereinafter, the operation of the first embodiment will be described.
まず、 撮像素子 1 3で生成されたアナログの画像データは、 A D変換部 1 4を 介してデジタル化された後、 信号処理部 1 5に与えられる。 信号処理部 1 5は、 この画像データに、 黒レベル補正、 およびガンマ補正などの信号処理を施す。 信 号処理後の画像データは、 バッファメモリ 1 6に一時記録される。 画像処理部 1 7は、 このバッファメモリ 1 6内の画像データを処理単位ずつ読み出し、 色補間 処理、 およぴ色座標変換などの画像処理を施す。  First, the analog image data generated by the image sensor 13 is digitized via the AD converter 14 and then supplied to the signal processor 15. The signal processing unit 15 performs signal processing such as black level correction and gamma correction on the image data. The image data after the signal processing is temporarily recorded in the buffer memory 16. The image processing unit 17 reads out the image data in the buffer memory 16 in units of processing and performs image processing such as color interpolation processing and color coordinate conversion.
画像圧縮部 1 8は、 バッファメモリ 1 6内から、 画像処理後の画像データを読 み出して、 画像圧縮処理を実施する。 図 2は、 この画像圧縮処理を説明する流れ 図である。  The image compression section 18 reads out the image data after the image processing from the buffer memory 16 and performs the image compression processing. FIG. 2 is a flowchart illustrating this image compression processing.
以下、 この図 2に示すステップ番号に沿って、 本発明の特徴である画像圧縮処 理を説明する。  Hereinafter, the image compression processing which is a feature of the present invention will be described along the step numbers shown in FIG.
ステップ S 1 : 画像圧縮部 1 8内のタイル処理部 2 1は、 バッファメモリ 1 6 内の画像データをタイル単位に分割して逐次読み出す。 Step S 1: The tile processing section 21 in the image compression section 18 divides the image data in the buffer memory 16 into tile units and sequentially reads the data.
ステップ S 2 : ウェーブレツト変換部 2 2は、 各タイルをウェーブレツト変換 し、 図 3に示すような複数のサブバンドに分ける。 Step S2: The wavelet transform unit 22 performs wavelet transform on each tile and divides the tile into a plurality of subbands as shown in FIG.
ステップ S 3 : ウエーブレット変換部 2 2は、 各タイルについて、 ウェーブレ ッ ト変換係数の標準偏差 (分散の平方根) をサブバンド毎に求める。 一般に、 こ れらの標準偏差が全体的に小さいほど、 画像データの階調 (色調も含む) の平坦 さは高い。 そのため、 画像データに僅かな階調ズレが生じても、 その階調ズレが 目立ちやすいと判断できる。 ウェーブレッ ト変換部 2 2は、 このように求めた標 準偏差を、 マイクロプロセッサ 3 9に伝達する。 Step S3: The wavelet transformer 22 finds the standard deviation (square root of variance) of the wavelet transform coefficient for each subband for each tile. In general, the smaller these standard deviations are overall, the higher the flatness of the gradation (including the color tone) of the image data. Therefore, even if a slight gradation shift occurs in the image data, it can be determined that the gradation shift is conspicuous. The wavelet transform unit 22 transmits the standard deviation thus obtained to the microprocessor 39.
ステップ S 4 : マイクロプロセッサ 3 9は、 変換係数の標準偏差に応じて、 画 像伸張に際してタイル境界に生じる不連続な圧縮歪みの目立ちやすさを予測する。 マイクロプロセッサ 3 9は、その予測状況に応じて、各サブバンドの視覚重み(視 覚的な重要度) を調整設定する。 Step S4: The microprocessor 39 selects an image according to the standard deviation of the conversion coefficient. Estimate the conspicuousness of discontinuous compression distortion generated at tile boundaries during image expansion. The microprocessor 39 adjusts and sets the visual weight (visual importance) of each subband according to the prediction situation.
以下、 この調整設定の具体例について説明する。  Hereinafter, a specific example of this adjustment setting will be described.
まず、 マイクロプロセッサ 39は、 各タイルの 3 L Lバンド (最低域バンド) の中で、 標準偏差が閾値以下 (ここでは 3以下) となる階調平坦タイルの割合を 求め、 その割合を r 1とする。 この 1の値が大きいほど、 画像データの階調は 大局的に平坦であり'、 階調ズレが目立ちやすいと判断できる。  First, the microprocessor 39 obtains the ratio of the gradation flat tile whose standard deviation is equal to or smaller than the threshold value (here, equal to or smaller than 3) in the 3 LL band (lowest band) of each tile, and determines the ratio as r 1. I do. As the value of 1 is larger, the gradation of the image data is globally flat, and it can be determined that the gradation deviation is more conspicuous.
さらに、 マイクロプロセッサ 3 9は、 各タイルの 3 LH, 3HL, 3 HHバン ドの中で標準偏差が閾値以下 (ここでは 3以下) となる階調平坦タイルの割合を 求め、 その割合を r 2とする。 この ι· 2は、 中間の空間周波数域において階調平 坦とみなせるタイルの割合である。 r lに加えて、 この r 2の値が大きいほど、 画像データの階調は更に平坦であり、階調ズレは更に目立ちやす 、と判断できる。 また、 マイクロプロセッサ 39は、 各タイルの 2 LH, 2 HL, 2HHパンド の中で標準偏差が閾値以下 (ここでは 3以下) となる階調平坦タイルの割合を求 め、 その割合を r 3とする。 この r 3も、 中間の空間周波数域において階調平坦 とみなせるタイルの割合である。 r 1および r 2に加えて、 こので 3の値が大き いほど、 画像データの階調は一層と平坦になり、 階調ズレは一段と目立ちやすい と判断できる。  Further, the microprocessor 39 calculates the ratio of the flat tone tiles whose standard deviation is equal to or less than the threshold value (here, 3 or less) in the 3 LH, 3HL, and 3 HH bands of each tile, and calculates the ratio as r 2 And This is the ratio of tiles that can be regarded as gradation flat in the intermediate spatial frequency range. In addition to rl, it can be determined that the larger the value of r2, the flatter the gradation of the image data and the more noticeable the gradation deviation. Further, the microprocessor 39 obtains the ratio of the gradation flat tile whose standard deviation is equal to or smaller than the threshold value (here, equal to or smaller than 3) among the 2 LH, 2 HL, and 2HH bands of each tile. I do. This r 3 is also the ratio of tiles that can be regarded as having a flat gradation in the intermediate spatial frequency range. In addition to r 1 and r 2, it can be determined that the larger the value of 3 is, the flatter the gradation of the image data is, and the more the gradation deviation is more conspicuous.
さらに、 マイクロプロセッサ 3 9は、 各タイルの 1 LH, 1 HL, 1 HHバン ドの中で標準偏差が閾値以下 (ここでは 3以下) となる階調平坦タイルの割合を 求め、 その割合を r 4とする。 この r 4は、 高域の空間周波数域において階調平 坦とみなせるタイルの割合である。 r l〜r 3に加えて、 この r 4の値が大きい ほど、 画像デ タの階調はシルクのように平坦であり、 階調ズレは非常に目立ち やすいと判断できる。  Further, the microprocessor 39 calculates the ratio of the gradation flat tile whose standard deviation is equal to or less than the threshold value (here, 3 or less) in the 1 LH, 1 HL, and 1 HH bands of each tile, and calculates the ratio as r And 4. This r4 is the ratio of tiles that can be regarded as gradation flat in the high spatial frequency range. In addition to r l to r 3, it can be determined that as the value of r 4 is larger, the gradation of the image data is flatter like silk, and the gradation deviation is more conspicuous.
このように求めた、 階調平坦タイルの割合 r 1〜4を用いて、 マイクロプロセ ッサ 3 9は、  Using the ratio r 1 to r 4 of the gradation flat tile obtained in this way, the microprocessor 39
パラメータ k = 8 (2X r 1+0.5 X r 2 + 0.25X r 3 + 0.25 X r 4) Parameter k = 8 (2Xr1 + 0.5Xr2 + 0.25Xr3 + 0.25Xr4)
を算出し、 パラメータ kを求める。 このパラメータ kは、 上述した r 1〜 r 4を 統合した評価パラメータである。 したがって、 この評価パラメータ kが大きいほ ど、 画像データの階調は、 低域から高域の空間周波数にかけて平坦であり、 階調 ズレが目立ちやすいと判断できる。 Is calculated, and the parameter k is obtained. This parameter k is the same as r1 to r4 described above. These are integrated evaluation parameters. Therefore, it can be determined that the larger the evaluation parameter k is, the flatter the gradation of the image data is from the low frequency band to the high frequency band, and that the gradation deviation is conspicuous.
このような評価パラメータ kの大きな画像データとしては、 例えば晴天時の青 空などの絵柄が想定される。 このような画像データは、 単調な絵柄のため、 タイ ル境界が浮き上がって特に目立ちやすいと予測される。  As such image data having a large evaluation parameter k, for example, a picture such as a blue sky in fine weather is assumed. Such image data is expected to be particularly conspicuous because the tile boundaries are raised because of the monotonous pattern.
一方、 評価パラメータ kの小さな画像データとしては、 ディテールを豊富に含 むような複雑な絵柄が想定される。 このような画像データは、 タイル境界に発生 する不連続な圧縮歪みが、 複雑な絵柄に埋もれて目立ちにくいと予測される。 マイクロプロセッサ 3 9は、 このような評価パラメータ kによる予測状況を加 味し、 次式を用いて各サブバンドの視覚重みを決定する。  On the other hand, image data with a small evaluation parameter k is assumed to be a complex picture rich in detail. In such image data, discontinuous compression distortions generated at tile boundaries are buried in complex patterns and are expected to be inconspicuous. The microprocessor 39 determines the visual weight of each subband by using the following equation in consideration of the prediction situation based on the evaluation parameter k.
y 3LL= 1.0 y 3LL = 1.0
γ 3LH= y 3HL= exp(-l/18 X k) γ 3LH = y 3HL = exp (-l / 18 X k)
7 3HH=exp(-2/18 Xk)  7 3HH = exp (-2/18 Xk)
y 2LH= y 2HL=exp(-4/18 X k) y 2LH = y 2HL = exp (-4/18 X k)
y 2HH= exp(-8/18 Xk) y 2HH = exp (-8/18 Xk)
y 1LH= y lHL= exp(-9/18 X k) y 1LH = y lHL = exp (-9/18 X k)
y lHH=exp(-18/18 Xk) y lHH = exp (-18/18 Xk)
上式では、 評価パラメータ kが大きくなるに従って、 低域サブバンドの視覚重 み γ 3Ι が相対的に大きく設定される。  In the above equation, as the evaluation parameter k increases, the visual weight γ 3Ι of the low band sub-band is set relatively large.
これは、 発明者が多数の画像圧縮実験に基づいて明らかにした 「低域の圧縮歪 みが発生するほど、 タイル境界が目立ちやすくなる」 という知見に基づく処置で ある。  This is a measure based on the findings that the inventor clarified based on a number of image compression experiments that "the more low-frequency compression distortion occurs, the more noticeable the tile boundary becomes."
マイクロプロセッサ 3 9は、 このように選択された視覚重みを、 量子化部 2 3 に伝達する。  The microprocessor 39 transmits the selected visual weight to the quantization unit 23.
ステップ S 5 : 量子化部 2 3は、 適当な符号化レートの拘束条件下で、 量子化 歪みによる視覚的な画質劣化を最小化するよう、 変換係数の量子化ステップ幅を 最適化する。 Step S5: The quantization unit 23 optimizes the quantization step width of the transform coefficient so as to minimize visual image quality deterioration due to quantization distortion under an appropriate coding rate constraint condition.
以下、 この量子化ステップ幅の最適化手法について、 原理的に説明する。 通常、 サブバンド毎の変換係数は、 ラプラス分布によく似た確率分布を示す。 そこで、 変換係数の確率分布を下記のラプラス分布式で近似し、 ラプラス分布を 特定する変数 αを求める。 f(y) = ^-e-^ ·'·[1] このような変換係数を、 量子化ステップ幅 Qで量子化した場合、 量子化後の変 換係数が k番目の量子化値を取る確率 p kは、
Figure imgf000011_0001
となる。 このとき、 量子化後の符号化レート R (Q) は、 理想的にはエントロピ 一に等しく、
Hereinafter, the method of optimizing the quantization step width will be described in principle. Normally, the transform coefficients for each subband show a probability distribution very similar to the Laplace distribution. Therefore, the probability distribution of the conversion coefficient is approximated by the following Laplace distribution equation, and the variable α that specifies the Laplace distribution is obtained. f (y) = ^ -e- ^ · '· [1] When such a transform coefficient is quantized with the quantization step width Q, the quantized transform coefficient takes the k-th quantized value. The probability pk is
Figure imgf000011_0001
Becomes At this time, the coding rate R (Q) after quantization is ideally equal to the entropy, and
R(Q) = -∑PklogPk R (Q) = -∑P k logP k
•[3]• [3]
Figure imgf000011_0002
となる。
Figure imgf000011_0002
Becomes
一方、 量子化歪み D (Q) は、 2乗誤差を用いて、  On the other hand, the quantization distortion D (Q) is calculated using the square error
D(Q)-∑L_ 2)(y-kQ†f(y)dy D (Q) -∑L_ 2) (y-kQ † f (y) dy
1/2)
Figure imgf000011_0003
1/2)
Figure imgf000011_0003
と表される。 It is expressed as
ここで、 i番目のサブバンドの量子化ステップ幅を Q i とし、 そのときの符号 化レートを R i (Q i ) とし、 量子化歪みを D i (R i ) と表す。  Here, the quantization step width of the i-th subband is Q i, the coding rate at that time is R i (Q i), and the quantization distortion is D i (R i).
ここでは、 サブバンドの符号化レートの総和を目標レートにするという拘束条 件のもとで、 量子化歪みによる視覚的な画質劣化を最小化する。 そこで、 ラグラ ンジェの未定乗数; を用いて、 次の関数 Jを導入する。  Here, the visual quality degradation due to quantization distortion is minimized under the constraint that the sum of the subband coding rates be the target rate. Therefore, the following function J is introduced using Lagrange's undetermined multiplier;
J =∑rA(Ri + ^∑ Qi) '··[5] この関数 Jの式中において、 γ iは、 サブバンド i に属する変換係数の視覚重 みに該当し、 上述したステップ S 4で既に求めた値群 γ 3Ι 〜γ 1HHである。 次に、 この関数 Jを R iで偏微分して、 関数 Jの第 1項が停留値となる条件式 を求めると、 , W 16] となる。 この式を変形することにより、J = ∑rA (Ri + ^ ∑ Qi) '· [5] In the expression of this function J, γ i corresponds to the visual weight of the transform coefficient belonging to subband i, and is a group of values γ 3Ι to γ 1HH already obtained in step S4 described above. Next, this function J is partially differentiated by R i to obtain a conditional expression in which the first term of the function J is a stationary value. By transforming this equation,
m M
Figure imgf000012_0001
Figure imgf000012_0001
が得られる。 Is obtained.
この [ 7 ] 式に適当な; Lの値を与えて解くことにより、 サブパンド毎の Q iが それぞれ求まる。 これらの Q iに基づいて、
Figure imgf000012_0002
を計算し、 トータルの符号化レート R λを求める。
By solving this equation [7], giving the value of L, Q i for each sub-band is obtained. Based on these Q i,
Figure imgf000012_0002
To calculate the total coding rate R λ.
この符号化レート Rえが目標符号化レートに一致するように、 未定乗数; Lを微 調整することにより、 量子化ステップ幅 Q iが確定する。  By fine-tuning the undetermined multiplier; L so that the coding rate R matches the target coding rate, the quantization step width Q i is determined.
(なお、 上述した最適化の演算結果をデータテーブルとして予めまとめておいて もよい。 この場合、 量子化部 2 3は、 このデータテーブルを、 変数 αおよび視覚 重みに基づいて参照することにより、 最適な量子化ステップ幅を即座に決定する ことができる。 )  (Note that the operation results of the above-described optimization may be summarized in advance as a data table. In this case, the quantization unit 23 refers to this data table based on the variable α and the visual weight, The optimal quantization step width can be determined immediately.)
ステップ S 6 : 量子化部 2 3は、 このように最適化された量子化ステップ幅 Q iに従って、 各サブバンドの変換係数を量子化する。 Step S 6: The quantization unit 23 quantizes the transform coefficient of each subband according to the quantization step width Q i thus optimized.
ところで、 上述した量子化ステップ幅の最適化では、 評価パラメータ kが大き いほど、 低域サブバンドの量子化ステップ幅は狭く設定される。 その結果、 評価 パラメータ kの大きな画像データほど、 低域の圧縮歪み (ここでは量子化歪み) が小さくなる。  By the way, in the above-described optimization of the quantization step width, the quantization step width of the lower sub-band is set to be narrower as the evaluation parameter k is larger. As a result, the lower the compression distortion (quantization distortion in this case), the lower the image data with a large evaluation parameter k.
ステップ S 7 : ビッ トモデリング部 2 4は、 R〇 I設定部 2 7に予め設定され る選択領域のマスク画像を読み出す。 ビッ トモデリング部 2 4は、 この選択領域 のマスク画像に基づいて、 サブバンド毎に変換係数が選択領域に含まれるか否か を判定する。 ビットモデリング部 2 4は、 選択領域内に位置する量子化後の変換 係数を、 Sビット分だけシフトする。 なお、 この Sビットは、 公知のマックスシ フト法に従って定められるビット数である。 Step S7: The bit modeling section 24 reads a mask image of the selected area preset in the R〇I setting section 27. The bit modeling section 24 It is determined whether or not the transform coefficient is included in the selected area for each subband based on the mask image of. The bit modeling unit 24 shifts the quantized transform coefficient located in the selected area by S bits. The S bits are the number of bits determined according to the known Max Shift method.
ステップ S 8 : ビットモデリング部 2 4は、 量子化後の変換係数をビットプレ ーンに分割する。 ビットモデリング部 2 4は、 最上位桁のビットプレーンから'順 に、 J P E G 2 0 0 0に規定される 3通りの符号化パスを実行する。 Step S8: The bit modeling unit 24 divides the quantized transform coefficient into bit planes. The bit modeling unit 24 executes three kinds of encoding passes specified in JPEG 2000 in order from the bit plane of the most significant digit.
ステップ S 9 : 算術符号化部 2 5は、 ビットモデリング部 2 4から符号化デー タを処理単位ごとに取得する。算術符号化部 2 5は、 この符号化データに対して、 2値算術符号である M Qコーダ一を用いて、 算術符号化を施す。 Step S 9: The arithmetic encoding unit 25 acquires encoded data from the bit modeling unit 24 for each processing unit. Arithmetic coding section 25 performs arithmetic coding on the coded data by using an MQ coder 1 which is a binary arithmetic code.
ステップ S 1 0 : ビットストリ一ム生成部 2 6は、 符号化データを S N Rプロ グレツシブその他の優先順位に従って並べ替え、 ビットス トリームを生成する。 ステップ S 1 1 : ビットストリーム生成部 2 6は、 このビットストリームを所 定の符号量で打ち切り、 圧縮ファイルとして出力する。 記録部 3 2は、 この圧縮 ファイルをメモリカード 3 3に記録保存する。 Step S10: The bitstream generation unit 26 rearranges the encoded data according to the SNR progressive and other priorities, and generates a bitstream. Step S11: The bitstream generation unit 26 cuts off this bitstream with a predetermined code amount and outputs it as a compressed file. The recording unit 32 records and stores the compressed file on the memory card 33.
上述した一連の動作により、 電子カメラ 1 1の撮像動作および画像圧縮処理が 完了する。  With the series of operations described above, the imaging operation of the electronic camera 11 and the image compression process are completed.
[第 1の実施形態の効果など]  [Effects of the first embodiment, etc.]
以上説明したように、 電子カメラ 1 1は、 ウェーブレット変換係数の標準偏差 を用いて、 タイル境界の目立ちやすさを予測する。 このとき、 タイル境界の目立 ちゃすい画像データほど、 低域の圧縮歪みを低く抑えるように量子化ステップ幅 を調整設定する。  As described above, the electronic camera 11 predicts the conspicuousness of the tile boundary using the standard deviation of the wavelet transform coefficient. At this time, the quantization step width is adjusted and set so that the image distortion of the tile boundary is less conspicuous and the compression distortion in the low frequency band is suppressed lower.
その結果、 タイル境界が視覚的に目立ちやすい画像データ (例えば青空などの 絵柄) を適切に選別し、 問題となる低域の圧縮歪みを予め抑えることが可能にな る。  As a result, it is possible to appropriately select image data (for example, a picture such as a blue sky) in which the tile boundaries are visually conspicuous, and to suppress in advance the problematic low-frequency compression distortion.
また逆に、 タイル境界が視覚的に目立ちにくい画像データ (例えば細かく複雑 な絵柄) を適切に選別し、 適度に低域の圧縮歪みを許容することが可能になる。 次に、 別の実施形態について説明する。  Conversely, it is possible to appropriately select image data (for example, fine and complicated patterns) in which tile boundaries are not visually noticeable, and to allow moderately low-frequency compression distortion. Next, another embodiment will be described.
《第 2の実施形態》 第 2の実施形態の構成は、 第 1の実施形態 (図 1 ) と同一であるため、 ここで の説明を省略する。 << 2nd Embodiment >> Since the configuration of the second embodiment is the same as that of the first embodiment (FIG. 1), the description is omitted here.
図 4は、 第 2の実施形態における画像圧縮処理を示す図である。 以下、 図 4に 示すステップ番号に沿って、 第 2の実施形態の動作を説明する。  FIG. 4 is a diagram illustrating an image compression process according to the second embodiment. Hereinafter, the operation of the second embodiment will be described along the step numbers shown in FIG.
ステップ S 2 1 : マイクロプロセッサ 3 9は、 圧縮対象の画像データを複数の 分割領域 (例えばタイル) に分けて抽出し、 分割領域ごとに画素値 (輝度、 色、 色差などの値) の標準偏差を求める。 Step S21: The microprocessor 39 extracts the image data to be compressed by dividing it into a plurality of divided areas (for example, tiles), and standard deviation of pixel values (values of luminance, color, color difference, etc.) for each divided area. Ask for.
ステップ S 2 2 : 抽出した分割領域全体に占める 『標準偏差が所定の閾値以下 となった分割領域』 の割合を求め、 その割合を画像データの平坦度 Rとする。 こ の平坦度 Rが大きいほど、 画像データの階調 (色調も含む) の平坦さが高く、 階 調ズレが目立ちやすいと判断できる。 Step S22: The ratio of the "segmented region where the standard deviation is equal to or less than a predetermined threshold" to the entire extracted divided region is determined, and the ratio is defined as the flatness R of the image data. It can be determined that the larger the flatness R is, the higher the flatness of the gradation (including the color tone) of the image data is and the more noticeable the gradation deviation is.
ステップ S 2 3 : この平坦度 Rが大きレ、ほど、 画像伸張に際してタイル境界に 生じる不連続な圧縮歪みが目立ちやすくなると予測される。 そこで、 マイクロプ 口セッサ 3 9は、 平坦度 Rの値が大きいほど、 低域サブバンドの視覚重みを相対 的に増大させる。 マイクロプロセッサ 3 9は、 このように調整設定された視覚重 みを量子化部 2 3に伝達する。 Step S23: It is predicted that the larger the flatness R is, the more noticeable the discontinuous compression distortion that occurs at the tile boundary during image expansion. Therefore, the microprocessor 39 relatively increases the visual weight of the low-frequency subband as the value of the flatness R increases. The microprocessor 39 transmits the visual weight adjusted and set as described above to the quantization unit 23.
ステップ S 2 4 : タイル処理部 2 1は、 ノ ッファメモリ 1 6内の画像データを タィル単位に分割して逐次読み出す。 Step S24: The tile processing section 21 divides the image data in the buffer memory 16 into tile units and sequentially reads them.
ステップ S 2 5 : ウェーブレット変換部 2 2は、 各タイルを独立にウェーブレ ット変換し、 図 3に示すようなサブバンドに分ける。 Step S25: The wavelet transform unit 22 independently performs wavelet transform on each tile and divides it into subbands as shown in FIG.
なお、 以降の動作 (ステップ S 2 6〜3 2 ) は、 第 1の実施形態で説明したス テツプ S 5〜S 1 1と同一動作のため、 ここでの説明を省略する。  Note that the subsequent operations (steps S26 to S32) are the same as steps S5 to S11 described in the first embodiment, and a description thereof will be omitted.
[第 2の実施形態の効果など]  [Effects of the second embodiment, etc.]
以上説明した動作により、 第 2の実施形態においても、 第 1の実施形態と同様 の効果を得ることができる。  By the operation described above, the same effects as those of the first embodiment can be obtained in the second embodiment.
特に、 第 2の実施形態では、 画素値の標準偏差に基づいて、 タイル境界の目立 ちゃすさを予測している。 その結果、 第 1の実施形態とは異なり、 ウェーブレツ ト変换を待たずに予測動作を実施することができる。 その結果、 タイル境界の目 立ちやすさの予測動作と、 ウェーブレット変換などを並列に処理できるなど、 処 理効率の向上に適している。 In particular, in the second embodiment, the conspicuousness of the tile boundary is predicted based on the standard deviation of the pixel values. As a result, unlike the first embodiment, the prediction operation can be performed without waiting for the wavelet change. As a result, the processing of predicting the conspicuousness of the tile boundaries and the wavelet transform can be processed in parallel. Suitable for improving efficiency.
以下、 別の実施形態について説明する。  Hereinafter, another embodiment will be described.
《第 3の実施形態》  << Third embodiment >>
第 3の実施形態の構成は、 第 1の実施形態 (図 1 ) と同一であるため、 ここで の説明を省略する。  Since the configuration of the third embodiment is the same as that of the first embodiment (FIG. 1), the description is omitted here.
図 5は、 第 3の実施形態における画像圧縮処理を示す図である。 以下、 図 5に 示すステップ番号に沿って、 第 3の実施形態の動作を説明する。  FIG. 5 is a diagram illustrating an image compression process according to the third embodiment. Hereinafter, the operation of the third embodiment will be described along the step numbers shown in FIG.
ステップ S 4 1 〜 S 4 2 : 第 1の実施形態で説明したステップ S 1 〜 S 2と同 じ。 Steps S41 to S42: Same as steps S1 to S2 described in the first embodiment.
ステップ S 4 3 : ウェーブレツト変換部 2 2は、 各タイルについて、 ウエーブ レッ ト変換係数の分散をサブバンド毎に求める。 下式は、 ウェーブレッ ト変換係 数の分散 (標準偏差の二乗) を算出するための式である。 び = ^∑[ ^ )- ]2 · ' · [9] Step S43: The wavelet transform unit 22 finds the variance of the wavelet transform coefficient for each subband for each tile. The following equation is used to calculate the variance (standard deviation squared) of the wavelet transform coefficient. And = ^ ∑ [^)-] 2 · '· [9]
,ヌ  , Nu
s '' 一 . s fucking one
2  Two
as ノハン r Sの分散 a s nohan r variance of S
Cs {x, y):サブバンド Sのタイル内位置 (JC, )の変換係数 C s (x, y): Conversion coefficient of position (JC,) in tile of subband S
Cj:サブバンド Sの変換係数の平均値 Cj: Average value of conversion coefficient of subband S
N:サブバンド Sの変換係数の個数 N: Number of transform coefficients of subband S
なお、 ここではタイル内部の全域について分散を算出している。 しかしながら、 タイル境界と関係するタイル内周に限って分散を算出してもよい (つまり、 タイ ル境界と関係しないタイルの中央域は分散演算から除いてもよい) 。 また、 タイ ルの各辺ごとにタイル境界の目立ちやすさを区別して予測するため、 タイルの上 辺近傍、 下辺近傍、左辺近傍、 および右辺近傍に分けて、 分散を算出してもよレ、。 図 6は、 テス ト画像 (人物画像) の特徴的な画像領域をピックアップして、 ゥ エープレツト変換係数の標準偏差をサブバンド別に示した図である。  Here, the variance is calculated for the entire area inside the tile. However, the variance may be calculated only for the inner circumference of the tile related to the tile boundary (that is, the central area of the tile not related to the tile boundary may be excluded from the variance calculation). In addition, in order to distinguish and predict the prominence of the tile boundary for each side of the tile, the variance may be calculated by dividing the tile into the vicinity of the top side, near the bottom side, near the left side, and near the right side. . FIG. 6 is a diagram in which characteristic image regions of a test image (person image) are picked up and the standard deviation of the エ ー -Apelet transform coefficient is shown for each sub-band.
このテス ト画像の場合、 『背景』 に該当するタイルは、 画像全体がぼかされて いるため、 高域サブバンドから低域サブバンドにかけて標準偏差の値が一様に低 く、 タイル境界が非常に目立ちやすいと判断できる。  In the case of this test image, the tile corresponding to the “background” has a uniformly low standard deviation from the high-frequency subband to the low-frequency subband because the entire image is blurred, and the tile boundary is It can be judged that it is very noticeable.
一方、 『ニットの服』 , 『目』 , 『口元』 のタイルは、 低域サブバンドから中 域サブバンドにかけて標準偏差の値が比較的高く、 タイル境界は比較的目立ちに くいと判断できる。 On the other hand, the tiles of “Knitted Clothes”, “Eyes”, and “Mouths” It can be determined that the standard deviation is relatively high over the area subband, and the tile boundaries are relatively inconspicuous.
また、 『髪の毛』 のタイルは、 低域サブバンドから高域サブバンドにかけて標 準偏差の値が全体的に高い。 そのため、 複雑な絵柄に埋もれて、 タイル境界は殆 ど目立たないと判断できる。  The “hair” tile has a high standard deviation value from the low-frequency subband to the high-frequency subband. Therefore, it can be determined that the tile boundaries are hardly noticeable because they are buried in a complicated pattern.
ウェーブレット変換部 2 2は、 このように求めた標準偏差を、 マイクロプロセ ッサ 3 9に伝達する。  The wavelet transform unit 22 transmits the standard deviation thus obtained to the microprocessor 39.
ステップ S 4 4 : マイクロプロセッサ 3 9は、 図 6で示したような評価判断を 計算で行うため、 タイル境界の目立ちやすさを数値で示す評価パラメータ kを求 める。 ここでの評価パラメータ kは、 第 1の実施形態とは異なり、 下式を用いて、 タィル別に算出された値である。 … ]Step S44: The microprocessor 39 obtains an evaluation parameter k which indicates the conspicuousness of the tile boundary by a numerical value in order to perform the evaluation judgment as shown in FIG. 6 by calculation. The evaluation parameter k here is a value calculated for each tile, using the following equation, different from the first embodiment. …]
Figure imgf000016_0001
Figure imgf000016_0001
[ 1 0 ] 式中の重み係数 W sは、 各サブバンド Sがタイル境界の目立ちやすさ に寄与する割合を示す値であり、 画像の主観評価実験などに基づいて決定すれば よい。 図 7は、 このように決定した重み係数 W sの一例を示す表である。 また、 [ 1 0 ] 式中の Aは比例定数である。 '  The weight coefficient W s in the expression [10] is a value indicating the ratio of each subband S contributing to the conspicuousness of the tile boundary, and may be determined based on a subjective evaluation experiment of an image or the like. FIG. 7 is a table showing an example of the weight coefficient Ws determined in this way. A in [10] is a proportionality constant. '
なお、 ステップ S 4 3において、 タイルを上辺、 下辺、 右辺、 および左辺に区 分して分散を計算している場合は、 これらの区分ごとに評価パラメータ kを求め ることが好ましい。  When the variance is calculated by dividing the tile into the upper side, the lower side, the right side, and the left side in step S43, it is preferable to obtain the evaluation parameter k for each of these sections.
ちなみに、 ここでは [ 1 0 ] 式を用いて評価パラメータ kを決定しているが、 本発明はこの評価式に限定されない。 分散の値が大きいほど、 タイル境界が目立 たないとの評価パラメータを算出する評価式であればよい。  Incidentally, here, the evaluation parameter k is determined using the [10] expression, but the present invention is not limited to this evaluation expression. An evaluation formula that calculates an evaluation parameter that makes the tile boundary less noticeable as the value of the variance is larger may be used.
図 8は、 テスト画像 (人物画像) の特徴的な画像領域をピックアップして、 評 価パラメータ kの値を示した表である。この評価パラメータ kの値が大きいほど、 タイル境界が目立ちやすいと予測できる。  Fig. 8 is a table showing the values of the evaluation parameter k by picking up characteristic image areas of the test image (person image). It can be predicted that the larger the value of the evaluation parameter k is, the more noticeable the tile boundary is.
ステップ S 4 5 : 続いて、 マイクロプロセッサ 3 9は、 隣接タイル間において 価パラメータ kの値を調整する。 ここでは、 マイクロプロセッサ 3 9力 処理 対象のタイルごとに上下左右の隣接タイルを含めた局所領域をまず設定する。 次 に、マイクロプロセッサ 39は、その局所領域内で評価パラメータ kの平滑化(単 純平均や重み付き平均など) を行う。 Step S45: Subsequently, the microprocessor 39 adjusts the value of the value parameter k between adjacent tiles. Here, a local area including the upper, lower, left, and right neighboring tiles is first set for each tile to be processed by the microprocessor. Next In addition, the microprocessor 39 smoothes the evaluation parameter k (such as a simple average and a weighted average) in the local region.
図 9 [A] は、 人物画像の評価パラメータ kをタイル区分に合わせて配列した マップである。 図 9 [B] は、 平滑化処理後の評価パラメータ kを示したマップ である。  Figure 9 [A] is a map in which the evaluation parameters k of the human image are arranged according to the tile division. Fig. 9 [B] is a map showing the evaluation parameter k after the smoothing process.
なお、 ここでの平滑化演算の代わりに、 MAX値, メディアン値などの算出演 算を使用してもよい。  Note that instead of the smoothing operation here, a calculation operation such as a MAX value or a median value may be used.
また、 ステップ S 4 3において、 タイルを上辺、 下辺、 右辺、 および左辺に区 分して標準偏差を計算している場合は、 実際に隣接する辺の間で評価パラメータ kの調整を行えばよい。 (例えば、 右に位置するタイルの左辺と、 左に位置する タイルの右辺との間で、 評価パラメータ kの調整を行えばよい。 )  When the standard deviation is calculated by dividing the tile into the upper side, the lower side, the right side, and the left side in step S43, the evaluation parameter k may be adjusted between actually adjacent sides. . (For example, the adjustment of the evaluation parameter k may be performed between the left side of the tile located on the right and the right side of the tile located on the left.)
このような評価パラメータ kの調整により、 タイル境界が目立ちやすいタィル または辺が存在した場合、 そのタイルまたは辺に接するタイルについてもタイル 境界が目立ちやすいとの判断がなされる。  By adjusting the evaluation parameter k as described above, if there is a tile or a side where the tile boundary is conspicuous, it is determined that the tile boundary is also conspicuous for the tile adjacent to the tile or the side.
ステップ S 46 : マイクロプロセッサ 39は、 このような評価パラメータ に よる予測状況を加味し、 次式を用いて各サブバンドの視覚重みを決定する。 y3LL=1.0 Step S46: The microprocessor 39 determines the visual weight of each sub-band using the following equation, taking into account the prediction situation based on such evaluation parameters. y3LL = 1.0
y 3LH= y3HL=exp(-l/18Xk) y 3LH = y3HL = exp (-l / 18Xk)
y 3HH=exp(-2/18Xk) y 3HH = exp (-2 / 18Xk)
y 2LH= γ 2HL=exp(-4/18 X k) y 2LH = γ 2HL = exp (-4/18 X k)
y2HH=exp(-8/18Xk) y2HH = exp (-8 / 18Xk)
γ 1LH= y lHL=exp(-9/18 X k) γ 1LH = y lHL = exp (-9/18 X k)
γ lHH=exp(-18/18Xk) γ lHH = exp (-18 / 18Xk)
上式では、 評価パラメータ kが大きくなるに従って、 低域サブバンドの視覚重 み T 3LL が相対的に大きく設定される。  In the above equation, as the evaluation parameter k increases, the visual weight T 3LL of the low-frequency subband is set to be relatively large.
マイクロプロセッサ 3 9は、 このように選択された視覚重みを、 量子化部 23 に伝達する。  The microprocessor 39 transmits the visual weight thus selected to the quantization unit 23.
ステップ S 47 : 量子化部 2 3は、 タイル別に設定された視覚重みに従って、 量子化歪みを視覚的な画質劣化を最小化するように、 変換係数の量子化ステップ をタイルごとに最適化する。 この最適化手法の手順については、 第 1の実施形態 のステップ S 5と同一であるため、 ここでの重複説明を省略する。 Step S47: The quantization unit 23 performs a quantization step of the transform coefficient according to the visual weight set for each tile so as to minimize the quantization distortion and the visual quality degradation. Is optimized for each tile. The procedure of this optimization method is the same as that in step S5 of the first embodiment, and thus the duplicate description is omitted here.
ステップ S 4 8 : 量子化部 2 3は、 タイル別に最適化された各サブバンドの量 子化ステップ幅を用いて、 タイル別に変換係数を量子化する。 Step S48: The quantizing section 23 quantizes the transform coefficient for each tile using the quantization step width of each subband optimized for each tile.
ところで、 上述した量子化ステップ幅の最適化では、 評価パラメータ kが大き いほど、 低域サブバンドの量子化ステップ幅は狭く設定される。 その結果、 評価 パラメータ kの大きなタイルほど、 低域の圧縮歪み (ここでは量子化歪み) が小 さくなる。  By the way, in the above-described optimization of the quantization step width, the quantization step width of the lower sub-band is set to be narrower as the evaluation parameter k is larger. As a result, tiles with a larger evaluation parameter k have lower compression distortion (quantization distortion in this case) in the low frequency band.
ステップ S 4 9〜S 5 3 : 第 1の実施形態で説明したステップ S 7〜S 1 1と 同じ。 Steps S49 to S53: Same as steps S7 to S11 described in the first embodiment.
ステップ S 5 4 : マイクロプロセッサ 3 9は、 画像伸張時のデータ使用に備え て、タイル別に決定した量子化ステップ幅に関するデータを圧縮データに付加し、 圧縮ファイルを生成する。 記録部 3 2は、 生成された圧縮ファイルをメモリカー ド 3 3に記録保存する。 Step S54: The microprocessor 39 adds data relating to the quantization step width determined for each tile to the compressed data to generate a compressed file in preparation for using data at the time of image expansion. The recording unit 32 records and stores the generated compressed file on the memory card 33.
上述した一連の動作により、 電子力メラ 1 1の撮像動作および画像圧縮処理が 完了する。  With the series of operations described above, the imaging operation of the electronic force camera 11 and the image compression processing are completed.
[第 3の実施形態の効果など]  [Effects of Third Embodiment]
以上説明した動作により、 第 3の実施形態においても、 第 1の実施形態と同様 の効果を得ることができる。  With the operation described above, the same effects as in the first embodiment can be obtained in the third embodiment.
特に、 第 3の実施形態では、 タイル境界の目立ちやすいタイル (例えば青空な どのタイル) を適切に選別し、 問題となる低域の圧縮歪み (すなわちタイル境界 の発生) をタイル単位に抑えることが可能になる。 また逆に、 タイル境界が視覚 的に目立ちにくいタイル (例えば絵柄の複雑なタイル) を適切に選別し、 適度に 低域の圧縮歪みを許容することが可能になる。 その結果、 絵柄の複雑なタイルに ついては、 タイル境界の発生を許容する分だけ、 中域〜高域周波数成分の画像情 報をなるベく正確に圧縮し、 圧縮による画像情報劣化をできる限り防ぐことが可 能になる。  In particular, in the third embodiment, tiles in which tile boundaries are conspicuous (for example, tiles such as blue sky) are appropriately selected and low-frequency compression distortion (that is, occurrence of tile boundaries) which is a problem is suppressed in tile units. Will be possible. Conversely, it is possible to appropriately select tiles whose tile boundaries are not visually conspicuous (for example, tiles with complicated patterns) and to allow moderately low-frequency compression distortion. As a result, for tiles with complicated patterns, the image information of the mid- and high-frequency components is compressed as accurately as possible to the extent that tile boundaries are allowed to occur, and image information deterioration due to compression is prevented as much as possible. It becomes possible.
このように第 3の実施形態では、 タイル境界の目立ちやすさを基準にして、 圧 縮符号量の周波数配分をタイル別に調整する。 その結果、 タイル境界が現実問題 となる極めて高い圧縮率において、 タイル境界を改善しつつ、 画像情報劣化を極 力抑えるといつた調整を実現することが可能になる。 As described above, in the third embodiment, the frequency distribution of the compression code amount is adjusted for each tile based on the conspicuousness of the tile boundary. As a result, tile boundaries are a real problem At extremely high compression ratios, it is possible to achieve the desired adjustment by minimizing image information degradation while improving tile boundaries.
さらに、 第 3の実施形態では、 タイル境界の目立つタイルの周囲についても、 低域の圧縮歪みを抑えている。 その結果、 平坦タイルの周囲に低域歪みが発生し て、 結果としてタイル境界が目立ってしまうといった事態を改善することが可能 になる。  Furthermore, in the third embodiment, low-frequency compression distortion is suppressed also around tiles with noticeable tile boundaries. As a result, it is possible to improve a situation in which low-frequency distortion occurs around a flat tile, and as a result, a tile boundary becomes conspicuous.
[実施形態の補足事項]  [Supplementary information of the embodiment]
なお、 上述した実施形態では、 変換係数や画素値の標準偏差 (分散の平方根) を求めて、 階調ズレの目立ちやすさを評価している。 しかしながら、 本発明はこ れに限定されるものではない。 一般に、 本発明の統計的性質としては、 『階調ズ レの目立ちやすさ』 と相関する量であればなんでもよい。例えば、階調の平坦さ、 階調の一様さ、 階調の単調さなどを示す統計的性質でもよい。  In the above-described embodiment, the conversion coefficient and the standard deviation (square root of the variance) of the pixel value are obtained to evaluate the conspicuousness of the gradation shift. However, the present invention is not limited to this. In general, any statistical property of the present invention may be used as long as the quantity correlates with “conspicuousness of gradation deviation”. For example, statistical properties indicating gradation flatness, gradation uniformity, gradation monotony, and the like may be used.
また、 本発明の 『分散』 は、 値の分散具合を示す尺度であればよく、 [ 9 ] 式 で示した統計学的な算出式に限定されるものではない。  In addition, the “variance” of the present invention may be any measure that indicates the degree of dispersion of values, and is not limited to the statistical calculation formula expressed by the formula [9].
ちなみに、 画素値のダイナミックレンジは、 瞬間的にでも大きな階調変化があ れば大きく変動してしまう値である。 そのため、 画素値のダイナミックレンジに ついては、 階調の平坦さを示す値ではなく、 『階調ズレの目立ちやすさ』 と相関 する量でもない。  By the way, the dynamic range of the pixel value is a value that fluctuates greatly even if there is a large gradation change even momentarily. Therefore, the dynamic range of the pixel value is not a value indicating the flatness of the gradation, and is not an amount correlating with “conspicuousness of gradation deviation”.
また、 上述した実施形態では、 画像圧縮パラメータとして、 変換係数の量子化 ステップ幅を調整設定している。 しかしながら、 本発明はこれに限定されるもの ではない。 一般に、 本発明の画像圧縮パラメータは、 タイル境界における圧縮歪 みの出現量を左右するパラメータであればなんでもよい。 例えば、 画像圧縮パラ メータとしては、 ビットプレーンの下位桁打ち切りを調整するパラメータや、 ビ ットス トリームの打ち切りを調整するパラメータや、 R O I符号化の選択領域を 調整するパラメータなどでもよい。  In the above-described embodiment, the quantization step size of the transform coefficient is adjusted and set as the image compression parameter. However, the present invention is not limited to this. In general, any image compression parameter of the present invention may be used as long as it affects the amount of appearance of compression distortion at a tile boundary. For example, the image compression parameter may be a parameter for adjusting the lower digit truncation of the bit plane, a parameter for adjusting the bit stream truncation, or a parameter for adjusting a selected region of ROI coding.
なお、 第 1および第 2の実施形態では、 複数のタイル全体に対して画像圧縮パ ラメータを一律に変更している。 この場合、 複数のタイル全体に均質に画像圧縮 を施すため、 画面上において圧縮歪みの現れ方にムラを生じるなどのおそれが少 ない。 ただし、 本発明はこれに限定されるものではない。 例えば、 第 3の実施形 態に示したように、 個々のタイルごとに画像圧縮パラメータを調整設定してもよ レ、。 また、 タイル境界とタイル内側とにおいて画像圧縮パラメータを個別に調整 設定してもよい。 In the first and second embodiments, the image compression parameters are uniformly changed for all the plurality of tiles. In this case, since the image compression is uniformly applied to the entire plurality of tiles, there is little possibility that the compression distortion will appear on the screen unevenly. However, the present invention is not limited to this. For example, the third embodiment As shown in the figure, the image compression parameters may be adjusted and set for each individual tile. Further, the image compression parameters may be individually adjusted and set at the tile boundary and inside the tile.
また、 上述した実施形態では、 電子カメラ 1 1において本発明の画像圧縮を実 施する場合について説明した。 しかしながら、 本発明はこれに限定されるもので はない。 例えば、 コンピュータ上において本発明の画像圧縮装置を実現してもよ い。 この場合、 上述した実施形態における画像圧縮の手順をプログラムコード化 することにより、 画像圧縮プログラム (請求項 9に対応) を作成すればよい。 なお、 本発明は、 その精神または主要な特徴から逸脱することなく、 他のいろ いろな形で実施することができる。 そのため、 前述の実施例はあらゆる点で単な る例示に過ぎず、 限定的に解釈してはならない。 本発明の範囲は、 特許請求の範 囲によって示すものであって、 明細書本文には、 なんら拘束されない。 さらに、 特許請求の範囲の均等範囲に属する変形や変更は、 すべて本宪明の範囲内のもの である。 産業上の利用の可能性  Further, in the above-described embodiment, the case where the image compression of the present invention is performed in the electronic camera 11 has been described. However, the present invention is not limited to this. For example, the image compression device of the present invention may be realized on a computer. In this case, an image compression program (corresponding to claim 9) may be created by converting the image compression procedure in the above-described embodiment into a program code. Note that the present invention may be embodied in various other forms without departing from the spirit or main characteristics thereof. Therefore, the above-described embodiment is merely an example in every aspect, and should not be construed as limiting. The scope of the present invention is defined by the scope of the claims, and is not limited by the text of the specification. Further, all modifications and changes belonging to the equivalent scope of the claims are within the scope of the present invention. Industrial applicability
本発明の画像圧縮では、 「階調ズレの視覚的な目立ちやすさ」 と相関する統計 的性質を画像データから求め、 この統計的性質により予測されるタイル境界の目 立ちやすさに応じて、 圧縮歪みの出現量を左右する画像圧縮パラメータを調整設 定する。  In the image compression of the present invention, a statistical property correlating with “visual conspicuousness of gradation shift” is obtained from image data, and in accordance with the conspicuous prominence of a tile boundary predicted by the statistical property, Adjust and set image compression parameters that affect the appearance of compression distortion.
したがって、 画像伸張時にタイル境界が目立ちそうな画像データに対しては、 圧縮歪みを予め抑えるなどの適切な対処が容易になる。  Therefore, for image data in which tile boundaries are likely to be conspicuous at the time of image expansion, appropriate measures such as suppressing compression distortion in advance are facilitated.

Claims

請求の範囲 The scope of the claims
( 1 ) 画像データを複数のタイルに区分し、 前記各タイルごとに画像圧縮を施 す画像圧縮装置であって、  (1) An image compression device that divides image data into a plurality of tiles and performs image compression for each of the tiles,
前記画像データを処理して、 「前記画像データに階調ズレを生じた場合の前記 階調ズレの目立ちやすさ」 と相関する統計的性質を求める統計処理部と、 前記統計的性質に基づいて画像伸張時におけるタイル境界の目立ちやすさを予 測し、 前記タイル境界の目立ちやすさに応じて 「圧縮歪みの出現量を左右する画 像圧縮パラメータ」 を調整設定するパラメータ設定部と、  A statistical processing unit that processes the image data and obtains a statistical property that is correlated with “conspicuousness of the tone shift when tone shift occurs in the image data”; A parameter setting unit for predicting the conspicuousness of the tile boundaries at the time of image expansion and adjusting and setting `` image compression parameters that affect the appearance of compression distortion '' according to the conspicuousness of the tile boundaries;
調整設定された前記画像圧縮パラメータを用いて、 前記画像データを画像圧縮 する画像圧縮部と  An image compression unit configured to perform image compression on the image data using the adjusted image compression parameter;
を備えたことを特徴とする画像圧縮装置。  An image compression device comprising:
( 2 ) 請求項 1に記載の画像圧縮装置において、  (2) In the image compression device according to claim 1,
前記統計的性質は、 前記画像データの階調の平坦さを示す値であり、 前記パラメータ設定部は、 前記平坦さの度合いが高いほど、 低域空間周波数成 分の圧縮歪みを軽減する方向に前記画像圧縮パラメータを調整設定する  The statistical property is a value indicating the flatness of the gradation of the image data.The higher the degree of flatness, the smaller the compression distortion of the low-frequency spatial frequency component. Adjust and set the image compression parameters
ことを特徴とする画像圧縮装置。  An image compression apparatus characterized by the above-mentioned.
( 3 ) 請求項 1に記載の画像圧縮装置において、  (3) In the image compression device according to claim 1,
前記統計的性質は、 前記複数のタイルに占める階調平坦タイル (平坦な階調と みなせるタイル) の割合であり、  The statistical property is a ratio of a gradation flat tile (tile that can be regarded as a flat gradation) to the plurality of tiles,
前記パラメータ設定部は、 前記階調平坦タイルの割合が多いほど、 低域空間周 波数成分の圧縮歪みを軽減する方向に前記画像圧縮パラメータを調整設定する ことを特徴とする画像圧縮装置。  The image compression device, wherein the parameter setting unit adjusts and sets the image compression parameter in a direction to reduce a compression distortion of a low-frequency spatial frequency component as the ratio of the gradation flat tile increases.
( 4 ) 請求項 1に記載の画像圧縮装置において、  (4) In the image compression device according to claim 1,
前記統計的性質は、 前記画像データを空間周波数領域に変換した変換係数につ いて求めた分散であり、  The statistical property is a variance obtained for a transform coefficient obtained by converting the image data into a spatial frequency domain.
前記パラメータ設定部は、 前記変換係数の分散が小さいほど、 低域空間周波数 成分の圧縮歪みを軽減する方向に前記画像圧縮パラメータを調整設定する ことを特徴とする画像圧縮装置。  The image compression device, wherein the parameter setting unit adjusts and sets the image compression parameter in a direction to reduce compression distortion of a low-frequency spatial frequency component as the variance of the transform coefficient is smaller.
( 5 ) 請求項 4に記載の画像圧縮装置において、 前記統計的性質は、 前記画像データをウェーブレツト変換した変換係数につい て求めた分散である (5) In the image compression device according to claim 4, The statistical property is a variance obtained for a transform coefficient obtained by performing a wavelet transform on the image data.
ことを特徴とする画像圧縮装置。  An image compression apparatus characterized by the above-mentioned.
( 6 ) 請求項 1ないし請求項 5のいずれか 1項に記載の画像圧縮装置において、 前記パラメータ設定部は、 前記複数のタイル全体について、 前記画像圧縮パラ メータを一律に設定する  (6) The image compression device according to any one of Claims 1 to 5, wherein the parameter setting unit uniformly sets the image compression parameters for the entire plurality of tiles.
ことを特徴とする画像圧縮装置。  An image compression apparatus characterized by the above-mentioned.
( 7 ) 請求項 1、 2、 4、 5のいずれか 1項に記載の画像圧縮装置において、 前記統計処理部は、 前記タイル単位に前記統計的性質を求め、  (7) The image compression device according to any one of claims 1, 2, 4, and 5, wherein the statistical processing unit obtains the statistical property for each tile.
前記パラメータ設定部は、 前記タイル単位の前記統計的性質に基づいて、 前記 タイル境界の目立ちやすさを前記タイル単位に予測し、 タイル境界の目立ちやす いタイルほど、 低域空間周波数成分の圧縮歪みを軽減する方向に前記画像圧縮パ ラメータを調整設定する  The parameter setting unit predicts the conspicuousness of the tile boundary in the tile unit based on the statistical property of the tile unit, and the compression distortion of the low-frequency spatial frequency component increases as the tile is more conspicuous. Adjust and set the image compression parameter in the direction to reduce
ことを特徴とする画像圧縮装置。  An image compression apparatus characterized by the above-mentioned.
( 8 ) 請求項 7に記載の画像圧縮装置において、  (8) In the image compression device according to claim 7,
前記パラメータ設定部は、 前記タィル境界が目立ちやすいタイルに隣接するタ ィルについても、 低域空間周波数成分の圧縮歪みを軽減する方向に前記画像圧縮 パラメータを調整設定する  The parameter setting unit adjusts and sets the image compression parameter in a direction to reduce compression distortion of a low-frequency spatial frequency component also for a tile adjacent to a tile in which the tile boundary is conspicuous.
ことを特徵とする画像圧縮装置。  An image compression apparatus characterized in that:
( 9 ) コンピュータを、 請求項 1ないし請求項 8のいずれか 1項に記載の前記 統計処理部、 前記パラメータ設定部、 および前記画像圧縮部として機能させるた めの画像圧縮プログラム。  (9) An image compression program for causing a computer to function as the statistical processing unit, the parameter setting unit, and the image compression unit according to any one of claims 1 to 8.
( 1 0 ) 画像データを複数のタイルに区分し、 前記各タイルごとに画像圧縮を 施す画像圧縮方法であって、  (10) An image compression method for dividing image data into a plurality of tiles and performing image compression on each of the tiles,
前記画像データを処理して、 「前記画像データに階調ズレを生じた場合の目立 ちゃすさ」 と相関する統計的性質を求める統計処理ステツプと、  A statistical processing step of processing the image data to determine a statistical property correlating with “inconspicuousness when a gradation shift occurs in the image data”;
前記統計的性質に基づいて画像伸張時におけるタイル境界の目立ちやすさを予 測し、 予測された前記タイル境界の目立ちやすさに応じて圧縮歪みの出現量を左 右する画像圧縮パラメータ」 を調整設定するパラメータ設定ステツプと. 調整設定された前記画像圧縮パラメータを用いて、 前記画像データを画像圧縮 する画像圧縮ステップと Predict the conspicuousness of the tile boundary at the time of image expansion based on the statistical properties, and adjust the `` image compression parameter that controls the appearance of compression distortion according to the predicted conspicuousness of the tile boundary '' Parameter setting step to be set. An image compression step of compressing the image data using the adjusted image compression parameter;
を有することを特徴とする画像圧縮方法。  An image compression method comprising:
PCT/JP2002/011353 2001-11-05 2002-10-31 Image compression apparatus, image compression program, and image compression method WO2003041394A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2003543305A JP4065978B2 (en) 2001-11-05 2002-10-31 Image compression apparatus, image compression program, and image compression method

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2001-339612 2001-11-05
JP2001339612 2001-11-05

Publications (1)

Publication Number Publication Date
WO2003041394A1 true WO2003041394A1 (en) 2003-05-15

Family

ID=19153938

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2002/011353 WO2003041394A1 (en) 2001-11-05 2002-10-31 Image compression apparatus, image compression program, and image compression method

Country Status (2)

Country Link
JP (1) JP4065978B2 (en)
WO (1) WO2003041394A1 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009503999A (en) * 2005-07-26 2009-01-29 ソニー株式会社 Encoding device, encoding method, decoding device, decoding method, and program thereof
JP2009278163A (en) * 2008-05-12 2009-11-26 Ricoh Co Ltd Image processing apparatus, image processing method, computer program and information recording medium
CN111046893A (en) * 2018-10-12 2020-04-21 富士通株式会社 Image similarity determining method and device, and image processing method and device

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0563990A (en) * 1991-08-30 1993-03-12 Fuji Xerox Co Ltd Encoder of image signal
JPH06217286A (en) * 1993-01-18 1994-08-05 Sanyo Electric Co Ltd Image encoding device
JPH06334873A (en) * 1993-05-19 1994-12-02 Mitsubishi Electric Corp Coding system for still picture
JPH0774965A (en) * 1993-06-15 1995-03-17 Ricoh Co Ltd Picture processing method
JPH07107481A (en) * 1993-09-30 1995-04-21 Sony Corp Method and device for encoding picture signal
JPH0969948A (en) * 1995-09-01 1997-03-11 Mita Ind Co Ltd Facsimile transmission control method
JPH10126777A (en) * 1996-10-21 1998-05-15 Fuji Photo Film Co Ltd Image compression method and its equipment
JP2000059558A (en) * 1998-08-05 2000-02-25 Ricoh Co Ltd Network scanner device

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0563990A (en) * 1991-08-30 1993-03-12 Fuji Xerox Co Ltd Encoder of image signal
JPH06217286A (en) * 1993-01-18 1994-08-05 Sanyo Electric Co Ltd Image encoding device
JPH06334873A (en) * 1993-05-19 1994-12-02 Mitsubishi Electric Corp Coding system for still picture
JPH0774965A (en) * 1993-06-15 1995-03-17 Ricoh Co Ltd Picture processing method
JPH07107481A (en) * 1993-09-30 1995-04-21 Sony Corp Method and device for encoding picture signal
JPH0969948A (en) * 1995-09-01 1997-03-11 Mita Ind Co Ltd Facsimile transmission control method
JPH10126777A (en) * 1996-10-21 1998-05-15 Fuji Photo Film Co Ltd Image compression method and its equipment
JP2000059558A (en) * 1998-08-05 2000-02-25 Ricoh Co Ltd Network scanner device

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009503999A (en) * 2005-07-26 2009-01-29 ソニー株式会社 Encoding device, encoding method, decoding device, decoding method, and program thereof
JP2009278163A (en) * 2008-05-12 2009-11-26 Ricoh Co Ltd Image processing apparatus, image processing method, computer program and information recording medium
US8355584B2 (en) 2008-05-12 2013-01-15 Ricoh Company, Limited Image processing apparatus, image processing method, and computer program product
CN111046893A (en) * 2018-10-12 2020-04-21 富士通株式会社 Image similarity determining method and device, and image processing method and device
CN111046893B (en) * 2018-10-12 2024-02-02 富士通株式会社 Image similarity determining method and device, image processing method and device

Also Published As

Publication number Publication date
JP4065978B2 (en) 2008-03-26
JPWO2003041394A1 (en) 2005-03-03

Similar Documents

Publication Publication Date Title
US7916961B2 (en) Compression encoder, compression encoding method and program
JP4594688B2 (en) Image encoding processing method, image decoding processing method, moving image compression processing method, moving image expansion processing method, image encoding processing program, image encoding device, image decoding device, image encoding / decoding system, extended image compression / decompression Processing system
JP2001510311A (en) Object-based rate control apparatus and method in coding scheme
US6327392B1 (en) Method of visual progressive coding
WO1999004359A1 (en) Apparatus and method for macroblock based rate control in a coding system
JP2003523652A (en) Adaptive block size DCT image compression based on change
WO1997017797A2 (en) Apparatus and method for quadtree based variable block size motion estimation
JP2003521178A (en) Quality-based image compression
JP2002335527A (en) Moving picture frame encoding method and adaptive quantization method and data compression rate deciding method and encoder
JP4045913B2 (en) Image coding apparatus, image coding method, and image processing apparatus
US8050506B2 (en) Image enhancement device
WO2012079237A1 (en) Wavelet coefficient quantization method applying human -eye visual model in image compression
JPH10336662A (en) Video coding method optimizing image quality of middle of picture, system and computer program product
JP2005530423A (en) Method and system for optimizing image sharpness during encoding
JP3743389B2 (en) Image compression apparatus and image compression program
AU762996B2 (en) An image compression system and method of determining quantisation parameters therefor
WO2003041394A1 (en) Image compression apparatus, image compression program, and image compression method
CA2656207A1 (en) Image processing method and apparatus, image processing program, and storage medium which stores the program
JP6792370B2 (en) Image processing equipment, its control method, and computer programs
EP1322116A2 (en) Image compression device, method for image compression, and electronic camera
CN104581157B (en) JPEG2000 encoding region of interest methods based on pre-cut
JP3792067B2 (en) Visual progressive coding method for images
AU2004225364B2 (en) Image compressing apparatus that achieves desired code amount
JP6813991B2 (en) Image coding device and its control method and program
JP4034236B2 (en) Specific transform coefficient priority coding apparatus, specific transform coefficient priority coding decoding system, and specific transform coefficient priority coding method

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NO NZ OM PH PL PT RO RU SD SE SG SI SK SL TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): GH GM KE LS MW MZ SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR IE IT LU MC NL PT SE SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 2003543305

Country of ref document: JP

122 Ep: pct application non-entry in european phase