文档章节

DCT, IDCT变换--C语言实现

t
 tigerBin
发布于 2017/06/14 12:35
字数 498
阅读 72
收藏 0

变换矩阵,DCT = IDCT'

static const double MtxDCT[8][8] = 
{
    {0.3536,    0.3536,    0.3536,    0.3536,    0.3536,    0.3536,    0.3536,    0.3536},
    {0.4904,    0.4157,    0.2778,    0.0975,   -0.0975,   -0.2778,   -0.4157,   -0.4904},
    {0.4619,    0.1913,   -0.1913,   -0.4619,   -0.4619,   -0.1913,    0.1913,    0.4619},
    {0.4157,   -0.0975,   -0.4904,   -0.2778,    0.2778,    0.4904,    0.0975,   -0.4157},
    {0.3536,   -0.3536,   -0.3536,    0.3536,    0.3536,   -0.3536,   -0.3536,    0.3536},
    {0.2778,   -0.4904,    0.0975,    0.4157,   -0.4157,   -0.0975,    0.4904,   -0.2778},
    {0.1913,   -0.4619,    0.4619,   -0.1913,   -0.1913,    0.4619,   -0.4619,    0.1913},
    {0.0975,   -0.2778,    0.4157,   -0.4904,    0.4904,   -0.4157,    0.2778,   -0.0975}
};


static const double MtxIDCT[8][8] = 
{
    {0.3536,    0.4904,    0.4619,    0.4157,    0.3536,    0.2778,    0.1913,    0.0975},
    {0.3536,    0.4157,    0.1913,   -0.0975,   -0.3536,   -0.4904,   -0.4619,   -0.2778},
    {0.3536,    0.2778,   -0.1913,   -0.4904,   -0.3536,    0.0975,    0.4619,    0.4157},
    {0.3536,    0.0975,   -0.4619,   -0.2778,    0.3536,    0.4157,   -0.1913,   -0.4904},
    {0.3536,   -0.0975,   -0.4619,    0.2778,    0.3536,   -0.4157,   -0.1913,    0.4904},
    {0.3536,   -0.2778,   -0.1913,    0.4904,   -0.3536,   -0.0975,    0.4619,   -0.4157},
    {0.3536,   -0.4157,    0.1913,    0.0975,   -0.3536,    0.4904,   -0.4619,    0.2778},
    {0.3536,   -0.4904,    0.4619,   -0.4157,    0.3536,   -0.2778,    0.1913,   -0.0975}
};

变换方程: F(u,v) = G * f * G'

逆变换方程: f(x,y) = G' * F * G

其中 G = DCT,  G' = IDCT ,  G 与 G' 互为转置

计算:

class MCU
{
public:
	int *operator[] (int row)
	{
		return data[row];
	}


	int data[8][8];
};


MCU MtxMulI2D(MCU &left, const double right[8][8])
{
	MCU dctBuf;
	for (int i = 0; i < 8; i++)
	{
		for (int j = 0; j < 8; j++)
		{
			double tempVal = 0.0;
			for (int k = 0; k < 8; k++)
			{
				tempVal += left[i][k] * right[k][j];
			}
			dctBuf[i][j] = round(tempVal);
		}
	}
	return dctBuf;
}

MCU MtxMulD2I(const double left[8][8], MCU &right)
{
	MCU dctBuf;
	for (int i = 0; i < 8; i++)
	{
		for (int j = 0; j < 8; j++)
		{
			double tempVal = 0.0;
			for (int k = 0; k < 8; k++)
			{
				tempVal += left[i][k] * right[k][j];
			}
			dctBuf[i][j] = round(tempVal);
		}
	}
	return dctBuf;
}

void DCT(MCU &block)
{
	block = MtxMulD2I(MtxDCT, block);
	block = MtxMulI2D(block, MtxIDCT);
}


void IDCT(MCU &block)
{
	block = MtxMulD2I(MtxIDCT, block);
	block = MtxMulI2D(block, MtxDCT);
}

测试:


int data2[8][8] = 
{
 {144,146,149,152,154,156,156,156},
 {148,150,152,154,156,156,156,156},
 {155,156,157,158,158,157,156,155},
 {160,161,161,162,161,159,157,155},
 {163,163,164,163,162,160,158,156},
 {163,163,164,164,162,160,158,157},
 {160,161,162,162,162,161,159,158},
 {158,159,161,161,162,161,159,158}
};

int data1[8][8] = 
{
 {139, 144, 149, 153, 155, 155, 155, 155}, 
 {144, 151, 153, 156, 159, 156, 156, 156}, 
 {150, 155, 160, 163, 158, 156, 156, 156}, 
 {159, 161, 162, 160, 160, 159, 159, 159}, 
 {159, 160, 161, 162, 162, 155, 155, 155}, 
 {161, 161, 161, 161, 160, 157, 157, 157}, 
 {162, 162, 161, 163, 162, 157, 157, 157}, 
 {162, 162, 161, 161, 163, 158, 158, 158} 
};


void Show(MCU &block)
{
	for (int i=0;i<8;i++)
	{
		for(int j=0;j<8;j++)
		{
			printf("%4d ", block[i][j]);
		}
		printf("\n");
	}
}




int main(int argc, char *argv[])
{
	MCU block;
	for (int i=0;i<8;i++)
	{
		for(int j=0;j<8;j++)
		{
			block[i][j] = data1[i][j];
		}
	}

	DCT(block);   // 离散余弦变换
	IDCT(block);  // 逆离散余弦变换
	Show(block);
	return 0;
}

最后计算结果与原始数据相差不大

© 著作权归作者所有

共有 人打赏支持
t
粉丝 0
博文 65
码字总数 17440
作品 0
武汉
程序员
x264源代码简单分析:宏块编码(Encode)部分

===================================================== H.264源代码分析文章列表: 【编码 - x264】 x264源代码简单分析:概述 x264源代码简单分析:x264命令行工具(x264.exe) x264源代码...

leixiaohua1020
2015/05/24
0
0
从零开始手敲次世代游戏引擎(JPEG特别篇)-1

上一篇我们导入了Aili这个相对比较复杂的模型,并发现了一些问题,进行了修正。 但是目前我们所渲染的只是一个没有任何材质的模型(也常常被称为白模,或者灰模)。为了能够输出类似Blender...

陈文礼
2017/12/01
0
0
x264源代码简单分析:编码器主干部分-1

===================================================== H.264源代码分析文章列表: 【编码 - x264】 x264源代码简单分析:概述 x264源代码简单分析:x264命令行工具(x264.exe) x264源代码...

leixiaohua1020
2015/05/11
0
0
【DCT】OPENCV python 离散余弦变换问题

因为项目需要,目前在研究《基于离散余弦变换和区域生长的白粉虱图像分割算法》,使用的OPENCV python。 论文地址:%e5%9f%ba%e4%ba%8e%e7%a6%bb%e6%95%a3%e4%bd%99%e5%bc%a6%e5%8f%98%e6%8d...

majisong
2016/06/02
1K
0
OpenCV - Operations on Arrays 对数组(矩阵)的一些操作

Function (函数名) Use (函数用处) add 矩阵加法,A+B的更高级形式,支持mask scaleAdd 矩阵加法,一个带有缩放因子dst(I) = scale * src1(I) + src2(I) addWeighted 矩阵加法,两个带有缩放...

苍海一粟
2012/11/20
0
2

没有更多内容

加载失败,请刷新页面

加载更多

TypeScript基础入门之JSX(二)

转发 TypeScript基础入门之JSX(二) 属性类型检查 键入检查属性的第一步是确定元素属性类型。 内在元素和基于价值的元素之间略有不同。 对于内部元素,它是JSX.IntrinsicElements上的属性类型...

durban
今天
1
0
AVA中CAS-ABA的问题解决方案AtomicStampedReference

了解CAS(Compare-And-Swap) CAS即对比交换,它在保证数据原子性的前提下尽可能的减少了锁的使用,很多编程语言或者系统实现上都大量的使用了CAS。 JAVA中CAS的实现 JAVA中的cas主要使用的是...

码代码的小司机
今天
2
0
Android JNI开发系列(十三) JNI异常处理

JNI 异常处理 JNI异常与JAVA处理异常的区别 JAVA 有异常处理机制,而JNI没有 如果JAVA中异常没有捕获,后面的代码不会执行,JNI会执行 JAVA编译时的异常,是在方法显示的声明了某一个异常,编...

蔡小鹏
今天
2
0
简单介绍Java 的JAR包、EAR包、WAR包区别

WAR包 WAR(Web Archive file)网络应用程序文件,是与平台无关的文件格式,它允许将许多文件组合成一个压缩文件。War专用于Web方面。大部分的JAVA WEB工程,都是打成WAR包进行发布的。 War是...

Linux就该这么学
今天
3
0
Qt那些事0.0.7

在帮助文档(Overview - QML and C++ Integration)中随缘遇到一张图,是关于C++对象与QML整合介绍的,值得标记下来,虽然大部分功能也有所涉猎,但是还是留个记号,万一哪天我失忆了还想写Q...

Ev4n
今天
4
0

没有更多内容

加载失败,请刷新页面

加载更多

返回顶部
顶部