数据压缩（8）：DPCM 压缩系统的实现和分析

1、DPCM编解码原理

DPCM是差分预测编码调制的缩写，是比较典型的预测编码系统。在DPCM系统中，需要注意的是预测器的输入是已经解码以后的样本。之所以不用原始样本来做预测，是因为在解码端无法得到原始样本，只能得到存在误差的样本。因此，在DPCM编码器中实际内嵌了一个解码器，如编码器中虚线框中所示。
在一个DPCM系统中，有两个因素需要设计：预测器和量化器。理想情况下，预测器和量化器应进行联合优化。实际中，采用一种次优的设计方法：分别进行线性预测器和量化器的优化设计。

2、代码实现
main.cpp

#include "stdafx.h"
using namespace std;
#include "targetver.h"
#include<stdio.h>
#include<string.h>
#include<math.h>
#include<malloc.h>
#include "DPCM.h"



int main(int argc, char *argv[])
{
FILE *LinnaFile=NULL;
if((LinnaFile=fopen("D:\\大三下\\数据压缩\\DPCM编码\\linna.yuv","rb"))==NULL)
{
	printf("open Linnafile fail!");
}
FILE *PredictedFile=NULL;
if((PredictedFile=fopen("D:\\大三下\\数据压缩\\DPCM编码\\PredictedFile.yuv","wb"))==NULL)
{
	printf("open ForecastFile fail!");
}
FILE *RebuildFile=NULL;
if((RebuildFile=fopen("D:\\大三下\\数据压缩\\DPCM编码\\RebuildFile.yuv","wb"))==NULL)
{
	printf("open RebuildFile fail!");
}
int w = 256;
int h = 256;
int N = 8;
unsigned char *ybuffer, *ubuffer, *vbuffer,*yPre,*yRe;
ybuffer = (unsigned char *)malloc(sizeof(char)*(w*h));
yPre = (unsigned char *)malloc(sizeof(char)*(w*h));
yRe = (unsigned char *)malloc(sizeof(char)*(w*h));
ubuffer = (unsigned char *)malloc(sizeof(char)*(w*h)/4);
vbuffer = (unsigned char *)malloc(sizeof(char)*(w*h)/4);

fread(ybuffer,1,w*h,LinnaFile);
fread(ubuffer,1,w*h/4,LinnaFile);
fread(vbuffer,1,w*h/4,LinnaFile);

 
//DPCM
DPCM(ybuffer,yPre,yRe,N,w,h);

//计算频率
//输出原始图像频率分布
FILE *LinnaFre=NULL;
if((LinnaFre=fopen("D:\\大三下\\数据压缩\\DPCM编码\\linnaFre.txt","wb"))==NULL)
{
	printf("open Linnafre fail!");
}
double fre1[256]={0};
for (int i = 0; i < 256; i++)
{
	for (int j = 0; j < w*h; j++)
	{
		if (ybuffer[j] == i)
		fre1[i]++;
	}
	fre1[i] = fre1[i] / (w*h);
	fprintf(LinnaFre, "%lf\n",fre1[i]);
}
//输出重建图像频率分布
FILE *RebuildFre=NULL;
if((RebuildFre=fopen("D:\\大三下\\数据压缩\\DPCM编码\\RebuildFre.txt","wb"))==NULL)
{
	printf("open RebuildFre fail!");
}
double fre2[256]={0};
for (int i = 0; i < pow(2.0, N); i++)
{
	for (int j = 0; j < w*h; j++)
	{
		if (yRe[j] == i)
		fre2[i]++;
	}
	fre2[i] = fre2[i] / (w*h);
	fprintf(RebuildFre, "%lf\n", fre2[i]);
}



//计算MSE和PSNR
double MSE=0;
double PSNR=0;
for (int i = 0; i <w*h ;i++)
{
MSE = MSE + pow((ybuffer[i] - yRe[i]), 2.0);
}
 
MSE = MSE / (w*h);
PSNR = 10 * log10((255 * 255) / MSE);
printf("N: %d\n", N);
//printf("MSE is %lf\n", MSE);
printf("PSNR is %lf\n", PSNR);


//写入预测误差图像
fwrite(yPre,1,w*h,PredictedFile);
fwrite(ubuffer,1,w*h/4,PredictedFile);
fwrite(vbuffer,1,w*h/4,PredictedFile);

 
//写入预测图像
fwrite(yRe,1,w*h,RebuildFile);
fwrite(ubuffer,1,w*h/4,RebuildFile);
fwrite(vbuffer,1,w*h/4,RebuildFile);

fclose(LinnaFile);
fclose(PredictedFile);
fclose(RebuildFile);
fclose(LinnaFre);
fclose(RebuildFre);
free(ybuffer);
free(ubuffer);
free(vbuffer);
free(yPre);
free(yRe);
return 0;
}

DPCM.cpp

// DPCM.cpp : 定义控制台应用程序的入口点。
//

#include "stdafx.h"
using namespace std;
#include "targetver.h"
#include<stdio.h>
#include<math.h>
#include<malloc.h>
#include "DPCM.h"

void DPCM(unsigned char *ybuffer,unsigned char *yPre,unsigned char *yRe,int N,int w,int h)
{
char *error;
error = (char *)malloc(sizeof(char)*(w*h));//量化前误差
int *Qerror;
Qerror = (int *)malloc(sizeof(int)*(w*h));//量化后误差
double Nbit = pow(2.0, N);//量化bite数
double step = 512 / Nbit;//量化间隔

for (int i = 0; i < h; i++)
{
		//第一列像素
		error[i*w + 0]=(char)(ybuffer[i*w + 0]-128);
		Qerror[i*w + 0]=int((error[i*w + 0]+255)/step+0.5);//第一列像素预测误差值等于0
		yPre[i*w + 0]=(unsigned char)Qerror[i*w + 0];
		if (yPre[i*w + 0] > Nbit-1)
		yPre[i*w + 0] =(unsigned char)(Nbit-1);
		if (yPre[i*w + 0] < 0)
		yPre[i*w + 0]= 0;
		yRe[i*w + 0] =(yPre[i*w + 0]*step-255)+128;
		
}
for (int i = 0; i < h; i++)
{
	for (int j = 1; j < w; j++)
	{
			error[i*w + j]=(char)(ybuffer[i*w + j]-yRe[i*w + j-1]);
			Qerror[i*w + j]=int((error[i*w + j]+255)/step+0.5);//量化
			yPre[i*w + j]=(unsigned char)Qerror[i*w + j];
			if (yPre[i*w + j] > Nbit-1)
			yPre[i*w + j] =(unsigned char)(Nbit-1);
			if (yPre[i*w + j] < 0)
			yPre[i*w + j]= 0;
			yRe[i*w + j]=(yPre[i*w + j]*step-255)+yRe[i*w + j-1];//反量化
			if (yRe[i*w + j] > 255)
			yRe[i*w + j] = 255;
			if (yRe[i*w + j] < 0)
			yRe[i*w + j]= 0;
	}
} 
}

DPCM.h

void DPCM(unsigned char *ybuffer,unsigned char *yPre,unsigned char *yRe,int N,int w,int h);

3.结果
原图：

原图概率分布：

8bit量化
预测误差图像与重建图像：

Y概率分布：

压缩质量以PSNR进行计算：

4bit量化
2bit量化

可以看出量化比特数越大图像质量越好。

4.huffman编码

DPCM+熵编码：

原图	192K
重建图	96K
hffman压缩后	71K
DPCM压缩比	2
DPCM+熵编码压缩比	2.7

仅进行熵编码：

原图	192K
压缩后	83K
压缩比	2.3

DPCM+熵编码对图像压缩最大。