CN103327337A - Classification quantization coding method based on bi-orthogonal lapped transform - Google Patents
Classification quantization coding method based on bi-orthogonal lapped transform Download PDFInfo
- Publication number
- CN103327337A CN103327337A CN2013102709431A CN201310270943A CN103327337A CN 103327337 A CN103327337 A CN 103327337A CN 2013102709431 A CN2013102709431 A CN 2013102709431A CN 201310270943 A CN201310270943 A CN 201310270943A CN 103327337 A CN103327337 A CN 103327337A
- Authority
- CN
- China
- Prior art keywords
- sigma
- image
- coefficient
- coding
- classification
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Compression Or Coding Systems Of Tv Signals (AREA)
- Compression, Expansion, Code Conversion, And Decoders (AREA)
Abstract
The invention provides a classification quantization coding method based on bi-orthogonal lapped transform, and belongs to the technical field of remote sensing image data transmission. According to the classification quantization coding method, the step of coding parameter designing is completed after image sequences are tested, and then the step of realizing an image coding algorithm to input images is conducted. The images are classified into different types according to coding characteristics through the method of classification training, an appropriate quantization method is selected for each type, therefore, the problem that image coding performances are different due to the fact that on the premise of quality-fixed coding, different types of remote sensing images have different requirements for coding performance parameters is solved, and the purpose of objective-quality-fixed coding of different types of images is achieved. The classification quantization coding method is applied to remote sensing satellites and greatly improves the downloading efficiency of remote sensing satellite image data.
Description
Technical field
The invention belongs to the remote sensing image data transmission technique field, be specifically related to a kind of classification quantitative coding method based on the biorthogonal lapped transform.
Background technology
Because the needs of observation, satellite sensor need to be transmitted back to ground by the image that observation system is captured.Along with user's demand, the resolution of image will be more and more higher, and this will cause the sharp increase of image data amount, and the data link channel capacity of at present space communication is limited.For so that ground can receive high-quality image, solve the contradiction of the transmission of remote sensing images big data quantity and limited interchannel by the Remote Sensing Image Coding method, be necessary.
Focus in the present Remote Sensing Image Coding method based on the embedded encoded method of discrete cosine transform and wavelet transformation.But in conventional method, the lifting of coding efficiency tends to be accompanied by the increase of complexity, and this brings more requirement with regard to the hardware designs of giving coded system.And in satellite system, hardware device all has a lot of restrictions on operational capability, internal memory and power consumption.Therefore, conventional method is difficult to satisfy the application demand of hardware system Real-time Collection transmission.Image coding technique based on the biorthogonal lapped transform provides a kind of effective method for realizing low complex degree, high-performance image coding, and it is adopted by JPEG XR standard.
In the remote sensing earth observation, widely different between the atural object in traditional coding method, carries out coding transmission with all cartographic features according to identical mode often.Therefore, the image of texture-rich tends to have more distortion than the simple image of texture.And for the user, tend to wish that various types of images all can satisfy quality requirement, and namely having similar coding distortion, this has just proposed new demand to method for encoding images.
The people such as Li have set up a kind of forecast model for describing the JPEG2000 coding quality by analysing in depth contacting between image liveness and the Image Coding performance.(referring to document: Ling Li and Zhen-Song Wang, Compression Quality Prediction Model for JPEG2000, IEEE Transactions on Image Processing, 2010) further can find, dissimilar image has different coding characteristics, and this feature can be weighed by the image liveness.
Summary of the invention
Problem for the background technology existence, the present invention proposes a kind of classification quantitative coding method based on the biorthogonal lapped transform, but efficient solution remote sensing images never of the same type are in the Image Coding performance difference problem of deciding under the quality coded prerequisite coding efficiency parameter request difference to be caused, thereby have reached the purpose that dissimilar image is decided quality coded.
For solving the problems of the technologies described above, the present invention adopts following technical scheme:
A kind of classification quantitative coding method based on the biorthogonal lapped transform comprises the steps,
Step 1, finish the design of coding parameter by the test pattern sequence;
Step 2, the coding parameter that utilizes step 1 to obtain are encoded to input picture.
Described step 1 specifically may further comprise the steps,
Step 1.1, test pattern is carried out the biorthogonal lapped transform;
Transform method adopts the transform method of JPEG XR standard, is specially: at first test pattern is divided into size and is 16 * 16 macro block, 4 * 4 piece carries out the biorthogonal lapped transform first time as the unit in each macro block; After the conversion, upper left angular data is the DC coefficient in each 4 * 4, and remaining 15 is the HP coefficient; The piece of DC coefficient compositions 4 * 4 all in 16 * 16 macro blocks is carried out the biorthogonal lapped transform second time, finally obtain a DC coefficient and 15 LP coefficients; Data after the final conversion of each macro block of 16 * 16 comprise 1 DC coefficient, 15 LP coefficients and 240 HP coefficients;
Step 1.2, the DC coefficient to behind each macroblock transform, LP coefficient and HP coefficient quantize, and quantitative formula is:
Wherein, A
IjCoefficient behind the representation transformation, B
IjCoefficient after representative quantizes, Q
StepRepresent quantization step, initial quantization step is set to 1, round () and represents to round up computing;
Step 1.3, to B
IjCarry out inverse quantization, inverse quantization formula is:
Ci
j=Bi
j×Q
step
Wherein, C
IjRepresent the dequantized coefficients of macro block, i, j represent respectively C
IjThe row, column coordinate;
Step 1.4, according to step 1.1 couple C
IjCarry out the inverse transformation of biorthogonal lapped transform, obtain reconstructed image;
The Y-PSNR (Peak Signal to Noise Ratio, PSNR) of step 1.5, calculating test pattern and reconstructed image, computing formula is:
Wherein, M, N be length and the width of representative image respectively, I (i, j) and
Represent respectively the pixel size of original image and reconstructed image, i, j represent respectively the row, column coordinate of original image and reconstructed image;
Step 1.6, to establish expection image objective quality be T
QIf, PSNR-T
Q>0.5, Q then
Step=Q
Step+ 1, then repeating step 1.2 is to step 1.6; If PSNR-T
Q<-0.5, Q then
Step=Q
Step-1, then repeating step 1.2 is to step 1.6; Other situations are finishing iteration computational process then;
Step 1.7, determine quantization step by iterative process, according to quantization step image is divided into dissimilarly, and set the quantization step of each type;
Step 1.8, computed image liveness comprise IAMD1, IAMD2, IAME1 and IAME2, and with it characteristic of division as image, the classification type in construction feature vector and the integrating step 1.7 is finished the training of SVMs;
Specific formula for calculation is as follows:
X (i, j) represents the pixel size of original image, and i, j represent respectively the row, column coordinate of original image.
Described step 2 specifically may further comprise the steps,
Step 2.1, the utilization method identical with step 1.1 are carried out the biorthogonal lapped transform to input picture;
The liveness of the method computed image that step 2.2, utilization are identical with step 1.8;
Step 2.3, combining image liveness make up the characteristic vector that is used for classification, and the classification based training result of integrating step 1.8 determines the type of input picture by SVMs, further determine quantization step;
Step 2.4, quantize according to the coefficient of quantization step after to input picture biorthogonal lapped transform;
Step 2.5, the coefficient after quantizing is carried out the entropy coding.
Compared with prior art, but the present invention's efficient solution remote sensing images never of the same type in the Image Coding performance difference problem of deciding under the quality coded prerequisite coding efficiency parameter request difference to be caused, thereby reached the purpose that dissimilar image is decided objective quality coding.Be used for remote sensing satellite and use, can greatly improve the lower transfer efficiency of remote sensing satellite image data.
Description of drawings
Fig. 1 is flow chart of the present invention;
Fig. 2 (a) is the typical test pattern in the category-A after the classification;
Fig. 2 (b) is the typical test pattern in the category-B after the classification;
Fig. 2 (c) is the typical test pattern in the C class after the classification;
Fig. 3 (a) is input picture Road;
Fig. 3 (b) is input picture Farmland;
Fig. 3 (c) is input picture Meadow;
Fig. 3 (d) is input picture City;
Fig. 4 (a) is decompressed image Road;
Fig. 4 (b) is decompressed image Farmland;
Fig. 4 (c) is decompressed image Meadow;
Fig. 4 (d) is decompressed image City;
Embodiment
Embodiment 1:
(1) set up the test pattern storehouse, picture number is 15 width of cloth;
(2) setting the objective image quality PSNR that expects is 35dB;
(3) calculate the image liveness value of test pattern, and by behind the biorthogonal lapped transform through after quantizing inverse transformation, near the quantization parameter value in the time of 35dB, its result is as shown in table 1 respectively;
Table 1
(4) resemble 1,2,3,4,5,11,12,13,14,15 according to large young pathbreaker's resolution chart of quantization parameter and be defined as category-A, test image 6,15 is defined as category-B, and test image 7,8,9,10 is divided into the C class;
(5) the image liveness value of different test patterns is used for the training of SVMs as characteristic value;
(6) quantization step of category-A test image being got average is 47, and it is that to get average be 66 for the quantization step of 59, C class testing image that the quantization step of category-B test image is got average;
(7) input picture comprises dissimilar image Road, Farmland, Meadow, City four width of cloth, and its liveness calculated value is as shown in table 2 respectively;
Table 2
Image | IAMD1 | IAMD2 | IAME1 | IAME2 |
Road | 21.7111 | 27.7085 | 33.7382 | 38.3088 |
Farmland | 7.9463 | 13.0992 | 12.4334 | 19.5948 |
Meadow | 8.8276 | 12.1177 | 13.4544 | 16.9056 |
City | 36.7683 | 48.1544 | 52.3366 | 62.4131 |
(8) according to the training result of SVMs in image liveness value and the integrating step 5, by SVMs input picture is classified.Classification results is: Road and City are divided into category-A; Meadow is divided into category-B; Farmland is divided into the C class;
(9) input picture is carried out biorthogonal lapped transform and quantification one by one, quantization step selects respectively 47,66,59 and 47.
(10) data after quantizing are finished the entropy coding.Entropy coding method has adopted the adaptive arithmetic code method.
(11) assess for the validity to coding method, the decoding data to behind the entropy coding obtains decoded reconstructed image.Weigh the difference of input picture and reconstructed image by PSNR, its result is respectively 35.45,35.06,35.02 and 34.53, can find out that worst error is no more than 1dB, has reached the purpose that the objective quality coding is decided in expection.
Claims (3)
1. classification quantitative coding method based on the biorthogonal lapped transform is characterized in that: comprises the steps,
Step 1, finish the design of coding parameter by the test pattern sequence;
Step 2, the coding parameter that utilizes step 1 to obtain are encoded to input picture.
2. a kind of classification quantitative coding method based on the biorthogonal lapped transform as claimed in claim 1, it is characterized in that: described step 1 specifically may further comprise the steps,
Step 1.1, test pattern is carried out the biorthogonal lapped transform;
Transform method adopts the transform method of JPEG XR standard, is specially: at first test pattern is divided into size and is 16 * 16 macro block, 4 * 4 piece carries out the biorthogonal lapped transform first time as the unit in each macro block; After the conversion, upper left angular data is the DC coefficient in each 4 * 4, and remaining 15 is the HP coefficient; The piece of DC coefficient compositions 4 * 4 all in 16 * 16 macro blocks is carried out the biorthogonal lapped transform second time, finally obtain a DC coefficient and 15 LP coefficients; Data after the final conversion of each macro block of 16 * 16 comprise 1 DC coefficient, 15 LP coefficients and 240 HP coefficients;
Step 1.2, the DC coefficient to behind each macroblock transform, LP coefficient and HP coefficient quantize, and quantitative formula is:
Wherein, A
IjCoefficient behind the representation transformation, B
IjCoefficient after representative quantizes, Q
StepRepresent quantization step, initial quantization step is set to 1, round () and represents to round up computing;
Step 1.3, to B
IjCarry out inverse quantization, inverse quantization formula is:
C
ij=B
ij×Q
step
Wherein, C
IjRepresent the dequantized coefficients of macro block, i, j represent respectively C
IjThe row, column coordinate;
Step 1.4, according to step 1.1 couple C
IjCarry out the inverse transformation of biorthogonal lapped transform, obtain reconstructed image;
The Y-PSNR (Peak Signal to Noise Ratio, PSNR) of step 1.5, calculating test pattern and reconstructed image, computing formula is:
Wherein, M, N be length and the width of representative image respectively, I (i, j) and
Represent respectively the pixel size of original image and reconstructed image, i, j represent respectively the row, column coordinate of original image and reconstructed image;
Step 1.6, to establish expection image objective quality be T
QIf, PSNR-T
Q>0.5, Q then
Step=Q
Step+ 1, then repeating step 1.2 is to step 1.6; If PSNR-T
Q<-0.5, Q then
Step=Q
Step-1, then repeating step 1.2 is to step 1.6; Other situations are finishing iteration computational process then;
Step 1.7, determine quantization step by iterative process, according to quantization step image is divided into dissimilarly, and set the quantization step of each type;
Step 1.8, computed image liveness comprise IAMD1, IAMD2, IAME1 and IAME2, and with it characteristic of division as image, the classification type in construction feature vector and the integrating step 1.7 is finished the training of SVMs;
Specific formula for calculation is as follows:
Wherein, x (i, j) represents the pixel size of original image, and i, j represent respectively the row, column coordinate of original image.
3. a kind of classification quantitative coding method based on the biorthogonal lapped transform as claimed in claim 1 or 2, it is characterized in that: described step 2 specifically may further comprise the steps,
Step 2.1, the utilization method identical with step 1.1 are carried out the biorthogonal lapped transform to input picture;
The liveness of the method computed image that step 2.2, utilization are identical with step 1.8;
Step 2.3, combining image liveness make up the characteristic vector that is used for classification, and the classification based training result of integrating step 1.8 determines the type of input picture by SVMs, further determine quantization step;
Step 2.4, quantize according to the coefficient of quantization step after to input picture biorthogonal lapped transform;
Step 2.5, the coefficient after quantizing is carried out the entropy coding.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310270943.1A CN103327337B (en) | 2013-06-28 | 2013-06-28 | A kind of classification quantitative coding method based on biorthogonal lapped transform |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310270943.1A CN103327337B (en) | 2013-06-28 | 2013-06-28 | A kind of classification quantitative coding method based on biorthogonal lapped transform |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103327337A true CN103327337A (en) | 2013-09-25 |
CN103327337B CN103327337B (en) | 2015-12-23 |
Family
ID=49195834
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310270943.1A Expired - Fee Related CN103327337B (en) | 2013-06-28 | 2013-06-28 | A kind of classification quantitative coding method based on biorthogonal lapped transform |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103327337B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104036293A (en) * | 2014-06-13 | 2014-09-10 | 武汉大学 | Rapid binary encoding based high resolution remote sensing image scene classification method |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1917562A (en) * | 2002-03-27 | 2007-02-21 | 微软公司 | System and method for progressively transforming and coding digital data |
CN103164713A (en) * | 2011-12-12 | 2013-06-19 | 阿里巴巴集团控股有限公司 | Image classification method and device |
-
2013
- 2013-06-28 CN CN201310270943.1A patent/CN103327337B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1917562A (en) * | 2002-03-27 | 2007-02-21 | 微软公司 | System and method for progressively transforming and coding digital data |
CN103164713A (en) * | 2011-12-12 | 2013-06-19 | 阿里巴巴集团控股有限公司 | Image classification method and device |
Non-Patent Citations (1)
Title |
---|
彭洲: "基于可逆整数时间域重叠变换的分类矢量量化图像编码", 《电子与信息学报》, vol. 33, no. 11, 30 November 2011 (2011-11-30), pages 2547 - 2552 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104036293A (en) * | 2014-06-13 | 2014-09-10 | 武汉大学 | Rapid binary encoding based high resolution remote sensing image scene classification method |
CN104036293B (en) * | 2014-06-13 | 2017-02-22 | 武汉大学 | Rapid binary encoding based high resolution remote sensing image scene classification method |
Also Published As
Publication number | Publication date |
---|---|
CN103327337B (en) | 2015-12-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100490494C (en) | Decoding apparatus, dequantizing method, distribution determining method, and program thereof | |
CN106385584B (en) | The adaptively sampled coding method of distributed video compressed sensing based on spatial correlation | |
CN101631243B (en) | Image encoding/decoding method based on wavelet transformation | |
CN101606320A (en) | The distortion of quantized data is estimated | |
CN102123278A (en) | Signal source encoding method based on distributed compressive sensing technology | |
CN108495132B (en) | The big multiplying power compression method of remote sensing image based on lightweight depth convolutional network | |
CN104427349A (en) | Bayer image compression method | |
JP5022471B2 (en) | Encoding method of wavelet image and corresponding decoding method | |
CN102281446A (en) | Visual-perception-characteristic-based quantification method in distributed video coding | |
CN102868907A (en) | Objective evaluation method for quality of segmental reference video | |
CN113079378B (en) | Image processing method and device and electronic equipment | |
CN105847800A (en) | Image compression method based on all phase discrete sine double orthogonal transformation and system | |
CN106056638B (en) | A kind of low complexity compression perceptual image coding method of adaptive texture contrast | |
CN108810534B (en) | Image compression method based on direction lifting wavelet and improved SPIHT under Internet of things | |
CN102857760B (en) | Feedback-free code rate optimization distributed video encoding and decoding method and system | |
Unde et al. | Rate–distortion analysis of structured sensing matrices for block compressive sensing of images | |
CN102592472A (en) | Parking-space detection method based on JPEG (joint photographic experts group) compressed image | |
CN104581158A (en) | Quantization table and image compression processing methods and devices, terminal and image searching system | |
CN102724495A (en) | Wyner-Ziv frame quantification method based on rate distortion | |
CN103327337B (en) | A kind of classification quantitative coding method based on biorthogonal lapped transform | |
CN106559668B (en) | A kind of low code rate image compression method based on intelligent quantization technology | |
CN101668204A (en) | Immune clone image compression method | |
CN104320661A (en) | Image coding quality predicting method based on difference entropy and structural similarity | |
CN105872536B (en) | A kind of method for compressing image based on dual coding pattern | |
CN115294222A (en) | Image encoding method, image processing method, terminal, and medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151223 Termination date: 20180628 |
|
CF01 | Termination of patent right due to non-payment of annual fee |