离散余弦正逆变换.pdfVIP

离散余弦正逆变换 本文来自：/pathuang68/archive/2009/05/17/4193990.aspx 最近网上有朋友因为要进行图像处理，而问及离散余弦正逆变换的问题，我在网上找了 一些资料，网下也找了一些参考书，发现居然有很多错误，尤其是逆变换(IDCT)的公式都是 错的！怪不得网友做完 DCT 后，再 IDCT ，却得不到原来的出使数据了 四舍五入的误差除 外) ，实在误人甚深。写此博文，以正其谬。 说明： 下面实现的代码，严格遵循离散余弦正逆变换的公式，仅供理解离散余弦逆变换的原理 而写。并未考虑代码优化，也未考虑快速算法，快速算法见快速离散余弦变换代码实现 (FDCT) 和快速逆离散余弦变换代码实现(FIDCT) 。 离散余弦变换(DCT)公式： N −1 N −1 2i 1 u  2j1 v  F u , v = u   v  ∑ ∑ f i , j cos [ ] cos [ ] 2N 2N i =0 j =0 1  = , 当=0 时  N 2  = , 当≠0 时  N 离散余弦变化(IDCT)公式： N −1 N −1 2i 1 u  2j1 v  f i , j =∑ ∑  u  v F u , v cos [ ] cos [ ] u =0 v =0 2N 2N 1  = , 当=0 时  N 2  = , 当≠0 时  N 注：严格的离散余弦变换公式中，其中一个N 可以写为 M ，且M 和N 可以不相等。这里由 于为图像处理服务，所以都写成了N 。 代码实现： #include iostream #include math.h using namespace std; #define NUM 8 #define PI 3.1415926 short round(double a) { if (a = 0) { return (short)(a + 0.5); } else { return (short)(a - 0.5); } } // DCT - Discrete Cosine Transform void DCT(short data[NUM][NUM]) { short output[NUM][NUM]; double ALPHA, BETA; short u = 0; short v = 0; short i = 0; short j = 0; for(u = 0; u NUM; u++) { for(v = 0; v NUM; v++) { if(u == 0) { ALPHA = sqrt(1.0 / NUM); }else { ALPHA = sqrt(2.0 / NUM); } if(v == 0) { BETA = sqrt(1.0 / NUM); }else { BETA = sqrt(2.0 / NUM); } double tmp = 0.0; for(i = 0;