一、实验简介
1.预测编码:
预测编码主要是利用信源相邻符合之间的相关性。根据某一模型利用以往的样本值对新样本进行预测,然后将样本的实际值与其预测值相减得到一个误差值,最后对这一误差值进行编码。
如果模型足够好,且样本序列在时间上相关性较强,则误差信号的幅度将远远小于原始信号,从而得到较大的数据压缩
2.DPCM:
DPCM是差分预测编码调制的缩写,是比较典型的预测编码系统。在DPCM系统中,需要注意的是预测器的输入是已经解码以后的样本。之所以不用原始样本来做预测,是因为在解码端无法得到原始样本,只能得到存在误差的样本。因此,在DPCM编码器中实际内嵌了一个解码器
3.预测误差:
二、实验原理
1.DPCM编码基本原理
符号 | 含义 |
---|---|
xk | 当前输入的信号样值 |
xk^ | 当前像素的预测值 |
xk’ | 当前像素的重建值 |
ek | 预测误差 |
ek^ | 量化后的预测误差 |
由以上可知,需要建立三个缓冲区,分别存原始图像、重建图像以及量化后的预测误差图像数据
2、DPCM编码系统的设计
在一个DPCM系统中,有两个因素需要设计:预测器和量化器。理想情况下,预测器
和量化器应进行联合优化。实际中,采用一种次优的设计方法:分别进行线性预测器和
量化器的优化设计。
3、量化比特数
量化比特数是指要区分所有量化级所需几位二进制数。一般来说都是2的n次方,如1、2、4、8等。量化比特数越大,画质越好。
4、PSNR和MSE
(1)MSE是“MeanSquareError”的缩写,表示原图像和量化图像之间的均方误差,指参数估计值与参数真值之差平方的期望值
(2)PSNR是“Peak Signal to Noise Ratio”的缩写,即峰值信噪比,是一种评价图像的客观标准,然而它是基于对应像素点间的误差,即基于误差敏感的图像质量评价。
- 由于并未考虑到人眼的视觉特性(人眼对空间频率较低的对比差异敏感度较高,人眼对亮度对比差异的敏感度较色度高,人眼对一个区域的感知结果会受到其周围邻近区域的影响等),因而经常出现评价结果与人的主观感觉不一致的情况。
- PSNR的单位是dB,数值越大表示失真越小
(3)PSNR和MSE计算公式如下:
5. Huffman编码
Huffman是一种无失真信源编码算法,主要目的是根据使用频率来最大化节省字符(编码)的存储空间,编码后可以得到即时的最佳码。
6. 编码效率与压缩量分析
- 压缩量依据PSNR分析
- 编码效率(压缩比和图像质量)
- 压缩比计算公式:
(1)原图像大小/DPCM+熵编码生成的预测误差图像大小
(2)原图像大小/仅进行熵编码生成的预测误差图像大小
三、实验过程
1.使用图像转换器把BMP文件转换为YUV文件。
2.读取一个256级的灰度图像YUV,提取Y的数据
3.分配缓冲区,存储数据
4.对该灰度图像进行预测,计算预测误差
5.对预测误差进行8bit、4bit、2bit、1bit均匀量化
6.对量化后的预测误差进行反量化,求重建值
7.将预测误差图像写入文件并输入Huffman编码器,得到输出码流、给出概率分布图并计算压缩比。
8.将原始图像文件输入Huffman编码器,得到输出码流、给出概率分布图并计算压缩比。
9.比较两种系统(1.DPCM+熵编码和2.仅进行熵编码)之间的编码效率(压缩比和图像质量)。压缩质量以PSNR进行计算。
四、实验代码
1.main.cpp
(1)准备工作
#define _CRT_SECURE_NO_WARNINGS
#include<stdio.h>
#include<math.h>
#include<malloc.h>
#include"dpcm.h"
int Width = 768;
int Height = 512;
int N = 8; //量化比特数
int main(int argc, char** argv)
{
FILE* srcFile;
FILE* errFile;
FILE* recFile;
FILE* fre1File;
FILE* fre2File;
FILE* fre3File;
char* srcFileName = NULL;
char* errFileName = NULL;
char* recFileName = NULL;
char* fre1FileName = NULL;
char* fre2FileName = NULL;
char* fre3FileName = NULL;
srcFileName = argv[1]; //原图yuv
errFileName = argv[2]; //预测误差图yuv
recFileName = argv[3]; //重建图yuv
fre1FileName = argv[4]; //原始图像频率分布txt
fre2FileName = argv[5]; //重建图像频率分布txt
fre3FileName = argv[6]; //预测误差图像频率分布txt
unsigned char* y, * u, * v, * re, * err;
double* fre1, * fre2, * fre3;
double MSE = 0, PSNR = 0;
//分配缓冲区
fre1 = (double*)malloc(sizeof(double) * 256);
fre2 = (double*)malloc(sizeof(double) * 256);
fre3 = (double*