DCT离散余弦变换C版本

原创 2015年07月08日 16:16:01
#include<iostream>
#include<math.h>
using namespace std;
#define PI 3.14159265354

#define DCT_SIZE 8

float DCT_Mat[DCT_SIZE][DCT_SIZE];
float DctMap[DCT_SIZE][DCT_SIZE]=
{
    130,130,130,129,134,133,129,130,
    130,130,130,129,134,133,130,130,
    130,130,130,129,132,132,130,130,
    129,130,130,129,130,130,129,129,
    127,128,127,129,131,129,131,130,
    127,128,127,128,127,128,132,132,
    125,126,129,129,127,129,133,132,
    127,125,128,128,126,130,131,131
};
   //需要DCT处理的图像输入口

/*float DctMap[DCT_SIZE][DCT_SIZE]=
{
    0,0,0,0,
    0,1,1,0,
    0,1,1,0,
    0,0,0,0
};*/

float DctMapTmp[DCT_SIZE][DCT_SIZE];

void InitDctMat()
{
    int i,j,t;
    for(i=0;i<DCT_SIZE;i++) DCT_Mat[0][i]=sqrt(2.0/DCT_SIZE)*sqrt(1.0/2);
    for(i=1;i<DCT_SIZE;i++)
    {
        t=0;
        for(j=0;j<DCT_SIZE;j++)
        {
            DCT_Mat[i][j]=sqrt(2.0/DCT_SIZE)*cos(PI/(2.0*DCT_SIZE)*(i+t));t+=2*i;
        }
    }
}

void DCT()
{
    float t=0;
    int i,j,k;
    for(i=0;i<DCT_SIZE;i++){
        for(j=0;j<DCT_SIZE;j++){
            t=0;
            for(k=0;k<DCT_SIZE;k++)
                t+=DCT_Mat[i][k]*DctMap[k][j];
            DctMapTmp[i][j]=t;
        }
    }
    for(i=0;i<DCT_SIZE;i++){
        for(j=0;j<DCT_SIZE;j++){
            t=0;
            for(k=0;k<DCT_SIZE;k++)
                t+=DctMapTmp[i][k]*DCT_Mat[j][k];
            DctMap[i][j]=t;
        }
    }
}

void IDCT()
{
    float t=0;
    int i,j,k;
    for(i=0;i<DCT_SIZE;i++){
        for(j=0;j<DCT_SIZE;j++){
            t=0;
            for(k=0;k<DCT_SIZE;k++)
                t+=DCT_Mat[k][i]*DctMap[k][j];
            DctMapTmp[i][j]=t;
        }
    }
    for(i=0;i<DCT_SIZE;i++){
        for(j=0;j<DCT_SIZE;j++){
            t=0;
            for(k=0;k<DCT_SIZE;k++)
                t+=DctMapTmp[i][k]*DCT_Mat[k][j];
            DctMap[i][j]=(int)(t);
        }
    }
}

int main()
{
    InitDctMat();
    DCT();
    for(int i=0;i<DCT_SIZE;i++)
    {
        for(int j=0;j<DCT_SIZE;j++)
        {
            cout<<DctMap[i][j]<<"  ";
        }
        cout<<endl;
    }
    IDCT();
    for(int i=0;i<DCT_SIZE;i++)
    {
        for(int j=0;j<DCT_SIZE;j++)
        {
            cout<<DctMap[i][j]<<"  ";
        }
        cout<<endl;
    }
}



版权声明:本文为博主原创文章,未经博主允许不得转载。 举报

相关文章推荐

OpenCV离散余弦变换原理与源码

环境配置:VS2012  OpenCV2.4.9 代码: #include "highgui.h" #include "math.h" #include "cv.h" #include "stdli...

离散余弦变换(FDCT和IDCT)

原文地址:http://hi.baidu.com/alphame/blog/item/2a4e52b7aa50a1f130add1a1.html 1.为什么要进行变换?        图像数据一般有较...

我是如何成为一名python大咖的?

人生苦短,都说必须python,那么我分享下我是如何从小白成为Python资深开发者的吧。2014年我大学刚毕业..

优化IPOL网站中基于DCT(离散余弦变换)的图像去噪算法(附源代码)。

优化IPOL网站中基于DCT(离散余弦变换)的图像去噪算法(附源代码)。   在您阅读本文前,先需要告诉你的是:即使是本文优化过的算法,DCT去噪的计算量依旧很大,请不要向这个算法提出实时运行的...

离散余弦变换(DCT)

离散余弦变换(DCT for Discrete Cosine Transform)是与傅里叶变换相关的一种变换,它类似于离散傅里叶变换(DFT for Discrete Fourier Trans...
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)