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+熵编码对图像压缩最大。