开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > Windows编程 > GDI/图象编程 > 查看源码
ifdct_8x8.c
资源名称：ifdct_8x8.rar　[点击查看]
上传用户：szuygs
上传日期：2009-09-28
资源大小：2k
文件大小：7k
源码类别：

GDI/图象编程

开发平台：

C/C++

ifdct_8x8.c：源码内容
  1. /**************************************************************************
  2. *          Copyright (c) 1993 Stefan Horbelt                               *
  3. *                                                                          *
  4. * Permission to use, copy, modify, and distribute this material for any    *
  5. * purpose and without fee is hereby granted, provided that the above       *
  6. * copyright notice and this permission notice appear in all copies.        *
  7. * WE MAKE NO REPRESENTATIONS ABOUT THE ACCURACY OR SUITABILITY OF THIS     *
  8. * MATERIAL FOR ANY PURPOSE.  IT IS PROVIDED "AS IS", WITHOUT ANY EXPRESS   *
  9. * OR IMPLIED WARRANTIES.                                                   *
  10. **************************************************************************/
  11. /**************************************************************************
  12. *                      *
  13. *  File    : ifdct_8x8.c                               *
  14. *  Date    : February 1993              *
  15. *  Author  : (c) Stefan Horbelt, horbelt@writeme.com      *
  16. *--------------------------------------------------------------------------*
  17. *  Inversion of the fast DCT program ffct.c of Thierry Turletti            *
  18. *  1992 INRIA Sophia Antipolis, FRANCE                                     *
  19. *--------------------------------------------------------------------------*
  20. *  Description : Perform a fast inverse DCT on a 8x8 block.                *
  21. *    by Inverted Matrices                                                  *
  22. *    The algorithm is refered to "A new Two-Dimensional Fast Cosine        *
  23. *    Transform Algorithm". IEEE Trans. on signal proc. vol 39 no. 2        *
  24. *    (S.C. Chan & K.L. Ho).                                                *
  25. *                                                                          *
  26. **************************************************************************/
  28. #include "h261.h"
  29. #define COL_MAX 360
  30. #define U_CHAR(x) ( (x)<-0.5 ? 0 : ( (x)>254.5 ? 255 : (u_char) (x+.5) ) )
  32. #define UnSurRac2  0.70710678
  33. #define Rac2Sur8   0.1767767
  34. #define UnSur8     0.125
  35. #define icf1_0      0.5097955792
  36. #define icf1_1      0.8999762231
  37. #define icf1_2      -2.562915448
  38. #define icf1_3      -0.6013448869
  39. #define icf2_0      0.5411961001
  40. #define icf2_1      -1.306562965
  41. #define Rac2    1.4142135625
  42. #define TwoRac2    2.828427125
  44. static float tmp0_0, tmp0_1, tmp0_2, tmp0_3, tmp0_4, tmp0_5, tmp0_6, tmp0_7;
  45. static float tmp1_0, tmp1_1, tmp1_2, tmp1_3, tmp1_4, tmp1_5, tmp1_6, tmp1_7;
  46. static float tmp2_0, tmp2_1, tmp2_2, tmp2_3, tmp2_4, tmp2_5, tmp2_6, tmp2_7;
  47. static float tmpblk[8][8];
  48. float *tmpblk_ptr;
  50. ifdct_8x8(image_in, block_out, lig, col)
  51. int image_in[8][8];  
  52. u_char block_out[][COL_MAX]; 
  53. int lig,col;
  54. {
  55.     register int u,v;
  57. #ifdef DEBUG
  58.     int      j,i;
  59.     fprintf(stderr,"nIFDCT:IMAGE_IN:n");
  60.     for (j=0;j<8;j++,fprintf(stderr,"n"))
  61.       for (i=0;i<8;fprintf(stderr,"%d, ",image_in[lig+j][col+i++]));
  62. #endif
  64. /*---------------------------------------------------------------*/
  66. image_in[0][0]=(image_in[0][0])<<3;
  68.       for(u=0; u<8; u++){
  69. /* 29 add; 17.5 mul (sum) 1D-IDCT*/
  71. /* this nonoptimised part is let inside the listing to make
  72.    the code mor transparent: *
  74. { register int *pt_tmpblk = &(image_in[u][0]);
  75. t0 = *pt_tmpblk++;
  76. t1 = *pt_tmpblk++;
  77. t2 = *pt_tmpblk++;
  78. t3 = *pt_tmpblk++;
  79. t4 = *pt_tmpblk++;
  80. t5 = *pt_tmpblk++;
  81. t6 = *pt_tmpblk++;
  82. t7 = *pt_tmpblk++;
  83. }
  85. tmp0_0 = t0;
  86. tmp0_1 =       t1;
  87. tmp0_2 =             t2;
  88. tmp0_3 =       t1      +  t3;
  89. tmp0_4 =                       t4;
  90. tmp0_5 =                  t3       +t5;
  91. tmp0_6 =             t2                 +t6;
  92. tmp0_7 =       t1+        t3       +t5       +t7;
  93. */
  95. /* same algorithmen as obove only optimised*/
  96. /* 5 add; 0 mul*/
  97. { register int *pt_tmpblk = &(image_in[u][0]),t3;
  98. tmp0_0 = *pt_tmpblk++;
  99. tmp0_1 = *pt_tmpblk++;
  100. tmp0_2 = *pt_tmpblk++;
  101. t3 = *pt_tmpblk++;
  102. tmp0_3 = tmp0_1 +t3;
  103. tmp0_4 = *pt_tmpblk++;
  104. tmp0_5 = t3 + *pt_tmpblk++;
  105. tmp0_6 = tmp0_2 + *pt_tmpblk++;
  106. tmp0_7 = tmp0_1 + tmp0_5 + *pt_tmpblk++;
  107. }
  109. if(u==0){/* 8 add; 10 mul*/
  111. tmp1_0 = (tmp0_0 + tmp0_4);
  112. tmp1_1 = (tmp0_0 - tmp0_4);
  114. tmp1_4 = (tmp0_1*Rac2 + tmp0_5);
  115. tmp1_5 = (tmp0_1*Rac2 - tmp0_5);
  117. tmp1_2 = (tmp0_2*Rac2 + tmp0_6)*icf2_0;
  118. tmp1_3 = (tmp0_2*Rac2 - tmp0_6)*icf2_1;
  120. tmp1_6 = (tmp0_3*Rac2 + tmp0_7)*icf2_0;
  121. tmp1_7 = (tmp0_3*Rac2 - tmp0_7)*icf2_1;
  122. }else{/* 8 add; 14 mul*/
  124. tmp1_0 = UnSurRac2*(tmp0_0 + tmp0_4);
  125. tmp1_1 = UnSurRac2*(tmp0_0 - tmp0_4);
  127. tmp1_4 = UnSurRac2*(tmp0_1*Rac2 + tmp0_5);
  128. tmp1_5 = UnSurRac2*(tmp0_1*Rac2 - tmp0_5);
  130. tmp1_2 = UnSurRac2*(tmp0_2*Rac2 + tmp0_6)*icf2_0;
  131. tmp1_3 = UnSurRac2*(tmp0_2*Rac2 - tmp0_6)*icf2_1;
  133. tmp1_6 = UnSurRac2*(tmp0_3*Rac2 + tmp0_7)*icf2_0;
  134. tmp1_7 = UnSurRac2*(tmp0_3*Rac2 - tmp0_7)*icf2_1;
  135. }
  137. /* 8 add; 4 mul*/
  138. tmp2_0 = (tmp1_0 + tmp1_2 );
  139. tmp2_2 = (tmp1_0 - tmp1_2 );
  141. tmp2_1 = (tmp1_1 + tmp1_3 );
  142. tmp2_3 = (tmp1_1 - tmp1_3 );
  144. tmp2_4 = (tmp1_4 + tmp1_6 )*icf1_0;
  145. tmp2_6 = (tmp1_4 - tmp1_6 )*icf1_2;
  147. tmp2_5 = (tmp1_5 + tmp1_7 )*icf1_1;
  148. tmp2_7 = (tmp1_5 - tmp1_7 )*icf1_3;
  150. { register float *pt_tmpblk = &(tmpblk[u][0]);
  151. /* 8 add; 0 mul*/
  152. *pt_tmpblk++ = (tmp2_0 + tmp2_4 );
  153. *pt_tmpblk++ = (tmp2_3 - tmp2_7 );
  154. *pt_tmpblk++ = (tmp2_1 + tmp2_5 );
  155. *pt_tmpblk++ = (tmp2_2 - tmp2_6 );
  156. *pt_tmpblk++ = (tmp2_2 + tmp2_6 );
  157. *pt_tmpblk++ = (tmp2_1 - tmp2_5 );
  158. *pt_tmpblk++ = (tmp2_3 + tmp2_7 );
  159. *pt_tmpblk++ = (tmp2_0 - tmp2_4 );
  160. }
  162.       }
  164. /*---------------------------------------------------------------*/
  166. for(v=0; v<8; v++){
  168.         register int colv=col+v;
  170. tmp0_0 =  tmpblk[0][v];
  171. tmp0_1 =  tmpblk[1][v];
  172. tmp0_2 =  tmpblk[2][v];
  173. tmp0_3 =  tmp0_1+ tmpblk[3][v];
  174. tmp0_4 =  tmpblk[4][v];
  175. tmp0_5 =  tmpblk[5][v]+tmpblk[3][v];
  176. tmp0_6 =  tmp0_2+tmpblk[6][v];
  177. tmp0_7 =  tmp0_3+tmpblk[5][v]+tmpblk[7][v];
  179. tmp1_0 = 0.125*(tmp0_0 + tmp0_4 *Rac2);
  180. tmp1_1 = 0.125*(tmp0_0 - tmp0_4 *Rac2);
  182. tmp1_4 = Rac2Sur8*(tmp0_1*Rac2 + tmp0_5);
  183. tmp1_5 = Rac2Sur8*(tmp0_1*Rac2 - tmp0_5);
  185. tmp1_2 = Rac2Sur8*(tmp0_2*Rac2 + tmp0_6)*icf2_0;
  186. tmp1_3 = Rac2Sur8*(tmp0_2*Rac2 - tmp0_6)*icf2_1;
  188. tmp1_6 = Rac2Sur8*(tmp0_3*Rac2 + tmp0_7)*icf2_0;
  189. tmp1_7 = Rac2Sur8*(tmp0_3*Rac2 - tmp0_7)*icf2_1;
  191. tmp2_0 = (tmp1_0 + tmp1_2 );
  192. tmp2_2 = (tmp1_0 - tmp1_2 );
  194. tmp2_1 = (tmp1_1 + tmp1_3 );
  195. tmp2_3 = (tmp1_1 - tmp1_3 );
  197. tmp2_4 = (tmp1_4 + tmp1_6 )*icf1_0;
  198. tmp2_6 = (tmp1_4 - tmp1_6 )*icf1_2;
  200. tmp2_5 = (tmp1_5 + tmp1_7 )*icf1_1;
  201. tmp2_7 = (tmp1_5 - tmp1_7 )*icf1_3;
  203. { register u_char *pt_tmpblk = &(block_out[lig][colv]);
  204. *pt_tmpblk = U_CHAR(tmp2_0 + tmp2_4 );
  205. pt_tmpblk+=COL_MAX;
  206. *pt_tmpblk = U_CHAR(tmp2_3 - tmp2_7 );
  207. pt_tmpblk+=COL_MAX;
  208. *pt_tmpblk = U_CHAR(tmp2_1 + tmp2_5 );
  209. pt_tmpblk+=COL_MAX;
  210. *pt_tmpblk = U_CHAR(tmp2_2 - tmp2_6 );
  211. pt_tmpblk+=COL_MAX;
  212. *pt_tmpblk = U_CHAR(tmp2_2 + tmp2_6 );
  213. pt_tmpblk+=COL_MAX;
  214. *pt_tmpblk = U_CHAR(tmp2_1 - tmp2_5 );
  215. pt_tmpblk+=COL_MAX;
  216. *pt_tmpblk = U_CHAR(tmp2_3 + tmp2_7 );
  217. pt_tmpblk+=COL_MAX;
  218. *pt_tmpblk = U_CHAR(tmp2_0 - tmp2_4 );
  219. }
  221.     }
  223. #ifdef DEBUG
  224. fprintf(stderr,"nIDFCT:n"); 
  225. for (j=0;j<8;j++,fprintf(stderr,"n"))
  226.   for (i=0;i<8;fprintf(stderr,"%d, ",(int)block_out[lig+j][col+i++]));    
  227. #endif
  228.   }