本篇是该系列的第五篇,承接上篇huffman解码,介绍接下来的两个步骤——反量化和逆zigzag变换,即IDCT前的两个步骤。
需要说明的是,这两个步骤可以颠倒,本人的实现是,先反量化,再逆ZigZag变换。
其实,这两步不需要太多说明,无非是查表对数据进行scale和数据重排,为了完整性,还是介绍一下吧。
1. 反量化
先拿到两个表,一个是解析文件头得到的量化表,另外一个是huffman解码得到的重建的8x8的block表。如下:
如何操作?对应位置相乘,即得到的反量化表:
2. 逆ZigZag变换
反量化表dequantization table是由如下ZigZag方式扫描得到的(忽略数字,关注箭头扫描方向):
因此,需要对该表进行恢复,得到:
如何操作?无非是查表,可以参考下面的代码实现:
1 static const uint8_t zigzag[64] = 2 { 3 0, 1, 5, 6, 14, 15, 27, 28, 4 2, 4, 7, 13, 16, 26, 29, 42, 5 3, 8, 12, 17, 25, 30, 41, 43, 6 9, 11, 18, 24, 31, 40, 44, 53, 7 10, 19, 23, 32, 39, 45, 52, 54, 8 20, 22, 33, 38, 46, 51, 55, 60, 9 21, 34, 37, 47, 50, 56, 59, 61, 10 35, 36, 48, 49, 57, 58, 62, 63 11 }; 12 13 // IZigZag 14 int JpegReZigZag(struct ABitReader *abr, struct jpegParam *param, int *dst_block, const int *src_block, bool dump) 15 { 16 for (int i=0; i<64; i++) { 17 dst_block[i] = src_block[zigzag[i]]; 18 } 19 20 if (dump) { 21 puts("----after rezigzag----"); 22 for (int i=0; i<64; i++) { 23 printf("%4d ", dst_block[i]); 24 if (i % 8 == 7) 25 puts(""); 26 } 27 puts(""); 28 } 29 30 return 0; 31 }