DPCM 压缩系统的实现和分析

一、实验简介

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*