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;
    }
}



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

DCT变换C语言源代码

  • 2010年08月31日 17:51
  • 11KB
  • 下载

一个C编写的快速DCT程序

  • 2010年06月10日 11:07
  • 2.25MB
  • 下载

C语言实现DCT变化

  • 2014年12月20日 18:57
  • 16.16MB
  • 下载

DCT变换编码及C语言实现

离散余弦变换(Discrete Cosine Transform,简称DCT变换)是一种与傅立叶变换紧密相关的数学运算。在傅立叶级数展开式中,如果被展开的函数是实偶函数,那么其傅立叶级数中只包含余弦...
  • yuechuxuan
  • yuechuxuan
  • 2017年07月13日 21:39
  • 608

离散余弦变换的C++实现

期末会有一个大作业是搞这个东西,自己还不是完全明白,先转载过来,有空再看吧。。。。 最近网上有朋友因为要进行图像处理,而问及离散余弦正逆变换的问题,我在网上找了一些资料,网下也找了一些参考书...
  • lwfcgz
  • lwfcgz
  • 2012年10月04日 23:19
  • 4553

离散余弦变换(C++实现)

理论部分转载自这篇blog: http://blog.csdn.net/luoweifu/article/details/8214959  该blog给出的是java代码,我用c++将其实现了。 理论...
  • Lu597203933
  • Lu597203933
  • 2015年04月17日 20:01
  • 3089

DCT变换及量化的c++实现(基于opencv矩阵运算)

由于DCT的数学原理不好描述,直接放代码了: #include #include #include #include using namespace std; using namespace cv;...
  • Qer_computerscience
  • Qer_computerscience
  • 2017年05月27日 23:17
  • 1444

离散余弦变换(含源码)

图像处理中常用的正交变换除了傅里叶变换外,还有其他一些有用的正交变换,其中离散余弦就是一种。离散余弦变换表示为DCT( Discrete Cosine Transformation),常用于图像处理和...
  • luoweifu
  • luoweifu
  • 2012年11月23日 09:38
  • 22078

DCT离散余弦变换编程

离散傅里叶变换和离散余弦变换公式如下: --------------------------------------------------------------------------------...
  • jfu22
  • jfu22
  • 2014年11月28日 14:21
  • 1258

离散余弦变换(DCT)

离散余弦变换(DCT for Discrete Cosine Transform)是与傅里叶变换相关的一种变换,它类似于离散傅里叶变换(DFT for Discrete Fourier Trans...
  • u012596983
  • u012596983
  • 2014年10月09日 09:40
  • 2216
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:DCT离散余弦变换C版本
举报原因:
原因补充:

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