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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > 图形图像处理 > 查看源码
Recon.cpp
资源名称：263.rar　[点击查看]
上传用户：szklck
上传日期：2007-01-22
资源大小：925k
文件大小：19k
源码类别：

图形图像处理

开发平台：

Visual C++

Recon.cpp：源码内容
  1. // Recon.cpp: implementation of the CRecon class.
  2. //
  3. //////////////////////////////////////////////////////////////////////
  5. #include "stdafx.h"
  6. #include "视频编解码器.h"
  7. #include "Recon.h"
  9. #ifdef _DEBUG
  10. #undef THIS_FILE
  11. static char THIS_FILE[]=__FILE__;
  12. #define new DEBUG_NEW
  13. #endif
  14. #define sign(a)         ((a) < 0 ? -1 : 1)
  16. //////////////////////////////////////////////////////////////////////
  17. // Construction/Destruction
  18. //////////////////////////////////////////////////////////////////////
  19. extern  int MV[2][5][MBR_MAX+1][MBC_MAX+2];
  20. extern int mv_outside_frame,adv_pred_mode,pb_frame,trd,trb;
  21. extern int modemap[MBR_MAX+1][MBC_MAX+2];
  22. extern int mb_width,coded_picture_width;
  23. extern  unsigned char *edgeframe[3],*newframe[3], *bframe[3],*oldrefframe[3];
  24. CRecon::CRecon()
  25. {
  27. }
  29. CRecon::~CRecon()
  30. {
  32. }
  34. void CRecon::reconstruct(int bx, int by, int P, int bdx, int bdy)
  35. {
  36.   int w,h,lx,lx2,dx,dy,xp,yp,comp,sum;
  37.   int x,y,mode,xvec,yvec;
  38.   unsigned char *src[3];
  40.   x = bx/16+1; 
  41.   y = by/16+1;
  42.   lx = coded_picture_width;
  44.   if (mv_outside_frame) //
  45.   {
  46.     lx2 = coded_picture_width + 64;
  47.     src[0] = edgeframe[0];
  48.     src[1] = edgeframe[1];
  49.     src[2] = edgeframe[2];
  50.   }
  51.   else 
  52.   {
  53.     lx2 = coded_picture_width;
  54.     src[0] = oldrefframe[0];
  55.     src[1] = oldrefframe[1];
  56.     src[2] = oldrefframe[2];
  57.   }
  59.   mode = modemap[y][x];
  61.   if (P) // P 预测模式
  62.   {
  63.     /* P prediction */
  64.     if (adv_pred_mode) 
  65. {
  66.       w = 8; h = 8;
  67.       /* Y*/
  68.       for (comp = 0; comp < 4; comp++) 
  69.   {
  70.       xp = bx + ((comp&1)<<3);
  71.       yp = by + ((comp&2)<<2);
  72.       recon_comp_obmc(src[0],newframe[0], lx,lx2,comp,w,h,xp,yp);
  73.       }
  74.       if (mode == MODE_INTER4V)
  75.   {
  76.     sum = MV[0][1][y][x]+MV[0][2][y][x]+ MV[0][3][y][x]+MV[0][4][y][x];
  77.     dx = sign(sum)*(roundtab[abs(sum)%16] + (abs(sum)/16)*2);
  78.     sum = MV[1][1][y][x]+MV[1][2][y][x]+ MV[1][3][y][x]+MV[1][4][y][x];
  79.     dy = sign(sum)*(roundtab[abs(sum)%16] + (abs(sum)/16)*2);
  80.       }
  81.       else 
  82.   {
  83.     dx = MV[0][0][y][x];
  84.     dy = MV[1][0][y][x];
  85.     // chroma rounding 
  86.     dx = ( dx % 4 == 0 ? dx >> 1 : (dx>>1)|1 );
  87.     dy = ( dy % 4 == 0 ? dy >> 1 : (dy>>1)|1 );
  88.       }
  89.       lx>>=1;bx>>=1; 
  90.   lx2>>=1; by>>=1;
  91.       // Chroma 
  92.       recon_comp(src[1],newframe[1], lx,lx2,w,h,bx,by,dx,dy,1);
  93.       recon_comp(src[2],newframe[2], lx,lx2,w,h,bx,by,dx,dy,2);
  94.     }
  95.     else 
  96. { // normal prediction mode 
  97.       // P prediction 
  98.       w = 16; h = 16;
  99.       dx = MV[0][0][y][x];
  100.       dy = MV[1][0][y][x];
  101.       
  102.       // Y image
  103.       recon_comp(src[0],newframe[0], lx,lx2,w,h,bx,by,dx,dy,0);
  104.       //src[0]->newframe[0]  
  105.       lx>>=1; w>>=1; 
  106.   bx>>=1; lx2>>=1; 
  107.       h>>=1; by>>=1;  
  108.       // chroma rounding 
  109.       dx = ( dx % 4 == 0 ? dx >> 1 : (dx>>1)|1 );
  110.       dy = ( dy % 4 == 0 ? dy >> 1 : (dy>>1)|1 );
  112.       // Chroma image
  113.       recon_comp(src[1],newframe[1], lx,lx2,w,h,bx,by,dx,dy,1);
  114.       recon_comp(src[2],newframe[2], lx,lx2,w,h,bx,by,dx,dy,2);
  115.   //src[]->newframe[]
  116.     }
  117.   }//end of(P)
  118.   else 
  119.   {
  120.     // B forward prediction 
  121.     if (adv_pred_mode) 
  122. {
  123.       if (mode == MODE_INTER4V) 
  124.   {
  125.      w = 8; h = 8;
  126.      // Y Image
  127.      xvec = yvec = 0;
  128.      for (comp = 0; comp < 4; comp++) 
  129.  {
  130.         xvec += dx = (trb)*MV[0][comp+1][y][x]/trd + bdx;
  131.         yvec += dy = (trb)*MV[1][comp+1][y][x]/trd + bdy;
  132.         xp = bx + ((comp&1)<<3);
  133.         yp = by + ((comp&2)<<2);
  134.         recon_comp(src[0],bframe[0], lx,lx2,w,h,xp,yp,dx,dy,0);
  135.  } 
  136.      // chroma rounding (table 16/H.263) 
  137.      dx = sign(xvec)*(roundtab[abs(xvec)%16] + (abs(xvec)/16)*2);
  138.      dy = sign(yvec)*(roundtab[abs(yvec)%16] + (abs(yvec)/16)*2);
  140.      lx>>=1; bx>>=1; 
  141.  lx2>>=1;by>>=1;
  142.      // Chroma  U V image 
  143.      recon_comp(src[1],bframe[1], lx,lx2,w,h,bx,by,dx,dy,1);
  144.      recon_comp(src[2],bframe[2], lx,lx2,w,h,bx,by,dx,dy,2);
  145.       } 
  146.       else
  147.   {  // adv_pred_mode but 16x16 vector 
  148.      w = 16; h = 16;
  149.      dx = (trb)*MV[0][0][y][x]/trd + bdx;
  150.      dy = (trb)*MV[1][0][y][x]/trd + bdy;
  151.      // Y image
  152.      recon_comp(src[0],bframe[0], lx,lx2,w,h,bx,by,dx,dy,0);
  154.      lx>>=1; w>>=1; bx>>=1; lx2>>=1;
  155.      h>>=1; by>>=1;  
  157.      xvec = 4*dx;
  158.      yvec = 4*dy;
  160.      // chroma rounding (table 16/H.263) 
  161.      dx = sign(xvec)*(roundtab[abs(xvec)%16] + (abs(xvec)/16)*2);
  162.      dy = sign(yvec)*(roundtab[abs(yvec)%16] + (abs(yvec)/16)*2);
  164.      // Chroma  U V image
  165.      recon_comp(src[1],bframe[1], lx,lx2,w,h,bx,by,dx,dy,1);
  166.      recon_comp(src[2],bframe[2], lx,lx2,w,h,bx,by,dx,dy,2);
  167.       }
  168.     }//end of(adv_pred_mode)
  169.     else 
  170. { // normal B forward prediction 
  172.       w = 16; h = 16;
  173.       dx = (trb)*MV[0][0][y][x]/trd + bdx;
  174.       dy = (trb)*MV[1][0][y][x]/trd + bdy;
  175.       /* Y */
  176.       recon_comp(src[0],bframe[0], lx,lx2,w,h,bx,by,dx,dy,0);
  178.       lx>>=1; w>>=1; bx>>=1; lx2>>=1;
  179.       h>>=1; by>>=1;  
  181.       xvec = 4*dx;
  182.       yvec = 4*dy;
  184.       // chroma rounding (table 16/H.263)
  185.       dx = sign(xvec)*(roundtab[abs(xvec)%16] + (abs(xvec)/16)*2);
  186.       dy = sign(yvec)*(roundtab[abs(yvec)%16] + (abs(yvec)/16)*2);
  188.       // Chroma U V image
  189.       recon_comp(src[1],bframe[1], lx,lx2,w,h,bx,by,dx,dy,1);
  190.       recon_comp(src[2],bframe[2], lx,lx2,w,h,bx,by,dx,dy,2);
  191.     }
  192.   }
  193. }
  195. void CRecon::recon_comp(unsigned char *src, unsigned char *dst, 
  196.     int lx, int lx2, int w, int h, int x, int y, int dx, int dy, int flag)
  197. {
  198.  int xint, xh, yint, yh;
  199.   unsigned char *s, *d;
  201.   xint = dx>>1;
  202.   xh = dx & 1;
  203.   yint = dy>>1;
  204.   yh = dy & 1;
  206.   /* origins */
  207.   s = src + lx2*(y+yint) + x + xint;
  208.   d = dst + lx*y + x;
  210.   if (!xh && !yh)
  211.     if (w!=8)
  212.       rec(s,d,lx,lx2,h);
  213.     else
  214.       recc(s,d,lx,lx2,h);
  215.   else if (!xh && yh)
  216.     if (w!=8)
  217.       recv(s,d,lx,lx2,h);
  218.     else 
  219.       recvc(s,d,lx,lx2,h);
  220.   else if (xh && !yh)
  221.     if (w!=8)
  222.       rech(s,d,lx,lx2,h);
  223.     else
  224.       rechc(s,d,lx,lx2,h);
  225.   else /* if (xh && yh) */
  226.     if (w!=8)
  227.       rec4(s,d,lx,lx2,h);
  228.     else
  229.       rec4c(s,d,lx,lx2,h);
  231. }
  233. void CRecon::recon_comp_obmc(unsigned char *src, unsigned char *dst, int lx, int lx2, int comp, int w, int h, int x, int y)
  234. {
  235. int j,k;
  236.   int xmb,ymb;
  237.   int c8,t8,l8,r8;
  238.   int xit,xib,xir,xil;
  239.   int yit,yib,yir,yil;
  240.   int vect,vecb,vecr,vecl;
  241.   int nx[5],ny[5],xint[5],yint[5],xh[5],yh[5];
  242.   int p[64],*pd;
  243.   unsigned char *d,*s[5];
  245.   xmb = (x>>4)+1;
  246.   ymb = (y>>4)+1;
  248.   c8 = (modemap[ymb][xmb] == MODE_INTER4V ? 1 : 0);
  249.   t8 = (modemap[ymb-1][xmb] == MODE_INTER4V ? 1 : 0);
  250.   l8 = (modemap[ymb][xmb-1] == MODE_INTER4V ? 1 : 0);
  251.   r8 = (modemap[ymb][xmb+1] == MODE_INTER4V ? 1 : 0);
  253.   switch (comp+1) {
  254.   case 1:
  255.     vect = (t8 ? 3 : 0) ; yit = ymb-1 ; xit = xmb;
  256.     vecb = (c8 ? 3 : 0) ; yib = ymb   ; xib = xmb;
  257.     vecl = (l8 ? 2 : 0) ; yil = ymb   ; xil = xmb-1;
  258.     vecr = (c8 ? 2 : 0) ; yir = ymb   ; xir = xmb;
  259.     if (ymb == 1) {
  260.       yit = ymb;
  261.       vect = (c8 ? 1 : 0);
  262.     }
  263.     if (xmb == 1) {
  264.       xil = xmb;
  265.       vecl = (c8 ? 1 : 0);
  266.     }
  267.     break;
  268.   case 2:
  269.     vect = (t8 ? 4 : 0) ; yit = ymb-1 ; xit = xmb;
  270.     vecb = (c8 ? 4 : 0) ; yib = ymb   ; xib = xmb;
  271.     vecl = (c8 ? 1 : 0) ; yil = ymb   ; xil = xmb;
  272.     vecr = (r8 ? 1 : 0) ; yir = ymb   ; xir = xmb+1;
  273.     if (ymb == 1) {
  274.       yit = ymb;
  275.       vect = (c8 ? 2 : 0);
  276.     }
  277.     if (xmb == mb_width) {
  278.       xir = xmb;
  279.       vecr = (c8 ? 2 : 0);
  280.     }
  281.     break;
  282.   case 3:
  283.     vect = (c8 ? 1 : 0) ; yit = ymb  ; xit = xmb;
  284.     vecb = (c8 ? 3 : 0) ; yib = ymb  ; xib = xmb;
  285.     vecl = (l8 ? 4 : 0) ; yil = ymb  ; xil = xmb-1;
  286.     vecr = (c8 ? 4 : 0) ; yir = ymb  ; xir = xmb;
  287.     if (xmb == 1) {
  288.       xil = xmb;
  289.       vecl = (c8 ? 3 : 0);
  290.     }
  291.     break;
  292.   case 4:
  293.     vect = (c8 ? 2 : 0) ; yit = ymb  ; xit = xmb;
  294.     vecb = (c8 ? 4 : 0) ; yib = ymb  ; xib = xmb;
  295.     vecl = (c8 ? 3 : 0) ; yil = ymb  ; xil = xmb;
  296.     vecr = (r8 ? 3 : 0) ; yir = ymb  ; xir = xmb+1;
  297.     if (xmb == mb_width) {
  298.       xir = xmb;
  299.       vecr = (c8 ? 4 : 0);
  300.     }
  301.     break;
  302.   default:
  303.     fprintf(stderr,"Illegal block number in recon_comp_obmc (recon.c)n");
  304.     exit(1);
  305.     break;
  306.   }
  308.   nx[0] = MV[0][c8 ? comp + 1 : 0][ymb][xmb];
  309.   ny[0] = MV[1][c8 ? comp + 1 : 0][ymb][xmb];
  310.   
  311.   nx[1] = MV[0][vect][yit][xit]; ny[1] = MV[1][vect][yit][xit];
  312.   nx[2] = MV[0][vecb][yib][xib]; ny[2] = MV[1][vecb][yib][xib];
  313.   nx[3] = MV[0][vecr][yir][xir]; ny[3] = MV[1][vecr][yir][xir];
  314.   nx[4] = MV[0][vecl][yil][xil]; ny[4] = MV[1][vecl][yil][xil];
  316.   for (k = 0; k < 5; k++) {
  317.     xint[k] = nx[k]>>1;
  318.     xh[k] = nx[k] & 1;
  319.     yint[k] = ny[k]>>1;
  320.     yh[k] = ny[k] & 1;
  321.     s[k] = src + lx2*(y+yint[k]) + x + xint[k];
  322.   }
  323.   
  324.   d = dst + lx*y + x;
  326.   if (!xh[0] && !yh[0])
  327.     reco(s[0],&p[0],8,lx2,0,0,0,8,0,8);
  328.   else if (!xh[0] && yh[0])
  329.     recvo(s[0],&p[0],8,lx2,0,0,0,8,0,8);
  330.   else if (xh[0] && !yh[0])
  331.     recho(s[0],&p[0],8,lx2,0,0,0,8,0,8);
  332.   else /* if (xh[] && yh[]) */
  333.     rec4o(s[0],&p[0],8,lx2,0,0,0,8,0,8);
  335.   if (!xh[1] && !yh[1])
  336.     reco(s[1],&p[0],8,lx2,1,1,0,8,0,4);
  337.   else if (!xh[1] && yh[1])
  338.     recvo(s[1],&p[0],8,lx2,1,1,0,8,0,4);
  339.   else if (xh[1] && !yh[1])
  340.     recho(s[1],&p[0],8,lx2,1,1,0,8,0,4);
  341.   else /* if (xh[] && yh[]) */
  342.     rec4o(s[1],&p[0],8,lx2,1,1,0,8,0,4);
  344.   if (!xh[2] && !yh[2])
  345.     reco(s[2]+(lx2<<2),&p[32],8,lx2,1,2,0,8,4,8);
  346.   else if (!xh[2] && yh[2])
  347.     recvo(s[2]+(lx2<<2),&p[32],8,lx2,1,2,0,8,4,8);
  348.   else if (xh[2] && !yh[2])
  349.     recho(s[2]+(lx2<<2),&p[32],8,lx2,1,2,0,8,4,8);
  350.   else /* if (xh[] && yh[]) */
  351.     rec4o(s[2]+(lx2<<2),&p[32],8,lx2,1,2,0,8,4,8);
  353.   if (!xh[3] && !yh[3])
  354.     reco(s[3],&p[0],8,lx2,1,3,4,8,0,8);
  355.   else if (!xh[3] && yh[3])
  356.     recvo(s[3],&p[0],8,lx2,1,3,4,8,0,8);
  357.   else if (xh[3] && !yh[3])
  358.     recho(s[3],&p[0],8,lx2,1,3,4,8,0,8);
  359.   else /* if (xh[] && yh[]) */
  360.     rec4o(s[3],&p[0],8,lx2,1,3,4,8,0,8);
  362.   if (!xh[4] && !yh[4])
  363.     reco(s[4],&p[0],8,lx2,1,4,0,4,0,8);
  364.   else if (!xh[4] && yh[4])
  365.     recvo(s[4],&p[0],8,lx2,1,4,0,4,0,8);
  366.   else if (xh[4] && !yh[4])
  367.     recho(s[4],&p[0],8,lx2,1,4,0,4,0,8);
  368.   else /* if (xh[] && yh[]) */
  369.     rec4o(s[4],&p[0],8,lx2,1,4,0,4,0,8);
  371.   pd = &p[0];
  372.   for (j = 0; j < 8; j++) {
  373.     d[0] = (pd[0] + 4 )>>3;
  374.     d[1] = (pd[1] + 4 )>>3;
  375.     d[2] = (pd[2] + 4 )>>3;
  376.     d[3] = (pd[3] + 4 )>>3;
  377.     d[4] = (pd[4] + 4 )>>3;
  378.     d[5] = (pd[5] + 4 )>>3;
  379.     d[6] = (pd[6] + 4 )>>3;
  380.     d[7] = (pd[7] + 4 )>>3;
  381.     d += lx;
  382.     pd += 8;
  383.   }
  385. }
  387. void CRecon::rec(unsigned char *s, unsigned char *d, int lx, int lx2, int h)
  388. {
  389.  int j;
  391.   for (j=0; j<h; j++)
  392.   {
  393.     d[0] = s[0];
  394.     d[1] = s[1];
  395.     d[2] = s[2];
  396.     d[3] = s[3];
  397.     d[4] = s[4];
  398.     d[5] = s[5];
  399.     d[6] = s[6];
  400.     d[7] = s[7];
  401.     d[8] = s[8];
  402.     d[9] = s[9];
  403.     d[10] = s[10];
  404.     d[11] = s[11];
  405.     d[12] = s[12];
  406.     d[13] = s[13];
  407.     d[14] = s[14];
  408.     d[15] = s[15];
  409.     s+= lx2;
  410.     d+= lx;
  411.   }
  412. }
  414. void CRecon::recc(unsigned char *s, unsigned char *d, int lx, int lx2, int h)
  415. {
  416.  int j;
  418.   for (j=0; j<h; j++)
  419.   {
  420.     d[0] = s[0];
  421.     d[1] = s[1];
  422.     d[2] = s[2];
  423.     d[3] = s[3];
  424.     d[4] = s[4];
  425.     d[5] = s[5];
  426.     d[6] = s[6];
  427.     d[7] = s[7];
  428.     s+= lx2;
  429.     d+= lx;
  430.   }
  432. }
  434. void CRecon::reco(unsigned char *s, int *d, int lx, int lx2, int addflag, int c, int xa, int xb, int ya, int yb)
  435. {
  436.   int i,j;
  437.   int *om;
  439.   om = &OM[c][ya][0];
  441.   if (!addflag) {
  442.     for (j = ya; j < yb; j++) {
  443.       for (i = xa; i < xb; i++) {
  444. d[i] = s[i] * om[i];
  445.       }
  446.       s+= lx2;
  447.       d+= lx;
  448.       om+= 8;
  449.     }
  450.   }
  451.   else {
  452.     for (j = ya; j < yb; j++) {
  453.       for (i = xa; i < xb; i++) {
  454. d[i] += s[i] * om[i];
  455.       }
  456.       s+= lx2;
  457.       d+= lx;
  458.       om+= 8;
  459.     }
  460.   }
  462. }
  464. void CRecon::rech(unsigned char *s, unsigned char *d, int lx, int lx2, int h)
  465. {
  466.   unsigned char *dp,*sp;
  467.   int j;
  468.   unsigned int s1,s2;
  470.   sp = s;
  471.   dp = d;
  472.   for (j=0; j<h; j++)
  473.   {
  474.     s1=sp[0];
  475.     dp[0] = (unsigned int)(s1+(s2=sp[1])+1)>>1;
  476.     dp[1] = (unsigned int)(s2+(s1=sp[2])+1)>>1;
  477.     dp[2] = (unsigned int)(s1+(s2=sp[3])+1)>>1;
  478.     dp[3] = (unsigned int)(s2+(s1=sp[4])+1)>>1;
  479.     dp[4] = (unsigned int)(s1+(s2=sp[5])+1)>>1;
  480.     dp[5] = (unsigned int)(s2+(s1=sp[6])+1)>>1;
  481.     dp[6] = (unsigned int)(s1+(s2=sp[7])+1)>>1;
  482.     dp[7] = (unsigned int)(s2+(s1=sp[8])+1)>>1;
  483.     dp[8] = (unsigned int)(s1+(s2=sp[9])+1)>>1;
  484.     dp[9] = (unsigned int)(s2+(s1=sp[10])+1)>>1;
  485.     dp[10] = (unsigned int)(s1+(s2=sp[11])+1)>>1;
  486.     dp[11] = (unsigned int)(s2+(s1=sp[12])+1)>>1;
  487.     dp[12] = (unsigned int)(s1+(s2=sp[13])+1)>>1;
  488.     dp[13] = (unsigned int)(s2+(s1=sp[14])+1)>>1;
  489.     dp[14] = (unsigned int)(s1+(s2=sp[15])+1)>>1;
  490.     dp[15] = (unsigned int)(s2+sp[16]+1)>>1;
  491.     sp+= lx2;
  492.     dp+= lx;
  493.   }
  495. }
  497. void CRecon::rechc(unsigned char *s, unsigned char *d, int lx, int lx2, int h)
  498. {
  499.   unsigned char *dp,*sp;
  500.   int j;
  501.   unsigned int s1,s2;
  503.   sp = s;
  504.   dp = d;
  505.   for (j=0; j<h; j++)
  506.   {
  507.     s1=sp[0];
  508.     dp[0] = (unsigned int)(s1+(s2=sp[1])+1)>>1;
  509.     dp[1] = (unsigned int)(s2+(s1=sp[2])+1)>>1;
  510.     dp[2] = (unsigned int)(s1+(s2=sp[3])+1)>>1;
  511.     dp[3] = (unsigned int)(s2+(s1=sp[4])+1)>>1;
  512.     dp[4] = (unsigned int)(s1+(s2=sp[5])+1)>>1;
  513.     dp[5] = (unsigned int)(s2+(s1=sp[6])+1)>>1;
  514.     dp[6] = (unsigned int)(s1+(s2=sp[7])+1)>>1;
  515.     dp[7] = (unsigned int)(s2+sp[8]+1)>>1;
  516.     sp+= lx2;
  517.     dp+= lx;
  518.   }
  520. }
  522. void CRecon::recho(unsigned char *s, int *d, int lx, int lx2, int addflag, int c, int xa, int xb, int ya, int yb)
  523. {
  524.  int *dp,*om;
  525.   unsigned char *sp;
  526.   int i,j;
  528.   sp = s;
  529.   dp = d;
  530.   om = &OM[c][ya][0];
  532.   if (!addflag) {
  533.     for (j = ya; j < yb; j++) {
  534.       for (i = xa; i < xb; i++) {
  535. dp[i] = (((unsigned int)(sp[i] + sp[i+1]+1))>>1)*om[i];
  536.       }
  537.       sp+= lx2;
  538.       dp+= lx;
  539.       om+= 8;
  540.     }
  541.   }
  542.   else {
  543.     for (j = ya; j < yb; j++) {
  544.       for (i = xa; i < xb; i++) {
  545. dp[i] += (((unsigned int)(sp[i] + sp[i+1]+1))>>1)*OM[c][j][i];
  546.       }
  547.       sp+= lx2;
  548.       dp+= lx;
  549.       om+= 8;
  550.     }
  551.   }
  553. }
  555. void CRecon::recv(unsigned char *s, unsigned char *d, int lx, int lx2, int h)
  556. {
  557.   unsigned char *dp,*sp,*sp2;
  558.   int j;
  560.   sp = s;
  561.   sp2 = s+lx2;
  562.   dp = d;
  563.   for (j=0; j<h; j++)
  564.   {
  565.     dp[0] = (unsigned int)(sp[0]+sp2[0]+1)>>1;
  566.     dp[1] = (unsigned int)(sp[1]+sp2[1]+1)>>1;
  567.     dp[2] = (unsigned int)(sp[2]+sp2[2]+1)>>1;
  568.     dp[3] = (unsigned int)(sp[3]+sp2[3]+1)>>1;
  569.     dp[4] = (unsigned int)(sp[4]+sp2[4]+1)>>1;
  570.     dp[5] = (unsigned int)(sp[5]+sp2[5]+1)>>1;
  571.     dp[6] = (unsigned int)(sp[6]+sp2[6]+1)>>1;
  572.     dp[7] = (unsigned int)(sp[7]+sp2[7]+1)>>1;
  573.     dp[8] = (unsigned int)(sp[8]+sp2[8]+1)>>1;
  574.     dp[9] = (unsigned int)(sp[9]+sp2[9]+1)>>1;
  575.     dp[10] = (unsigned int)(sp[10]+sp2[10]+1)>>1;
  576.     dp[11] = (unsigned int)(sp[11]+sp2[11]+1)>>1;
  577.     dp[12] = (unsigned int)(sp[12]+sp2[12]+1)>>1;
  578.     dp[13] = (unsigned int)(sp[13]+sp2[13]+1)>>1;
  579.     dp[14] = (unsigned int)(sp[14]+sp2[14]+1)>>1;
  580.     dp[15] = (unsigned int)(sp[15]+sp2[15]+1)>>1;
  581.     sp+= lx2;
  582.     sp2+= lx2;
  583.     dp+= lx;
  584.   }
  586. }
  588. void CRecon::recvc(unsigned char *s, unsigned char *d, int lx, int lx2, int h)
  589. {
  590.  unsigned char *dp,*sp,*sp2;
  591.   int j;
  593.   sp = s;
  594.   sp2 = s+lx2;
  595.   dp = d;
  597.   for (j=0; j<h; j++)
  598.   {
  599.     dp[0] = (unsigned int)(sp[0]+sp2[0]+1)>>1;
  600.     dp[1] = (unsigned int)(sp[1]+sp2[1]+1)>>1;
  601.     dp[2] = (unsigned int)(sp[2]+sp2[2]+1)>>1;
  602.     dp[3] = (unsigned int)(sp[3]+sp2[3]+1)>>1;
  603.     dp[4] = (unsigned int)(sp[4]+sp2[4]+1)>>1;
  604.     dp[5] = (unsigned int)(sp[5]+sp2[5]+1)>>1;
  605.     dp[6] = (unsigned int)(sp[6]+sp2[6]+1)>>1;
  606.     dp[7] = (unsigned int)(sp[7]+sp2[7]+1)>>1;
  607.     sp+= lx2;
  608.     sp2+= lx2;
  609.     dp+= lx;
  610.   }
  611. }
  613. void CRecon::recvo(unsigned char *s, int *d, int lx, int lx2, int addflag, int c, int xa, int xb, int ya, int yb)
  614. {
  615.  int *dp,*om;
  616.   unsigned char *sp,*sp2;
  617.   int i,j;
  619.   sp = s;
  620.   sp2 = s+lx2;
  621.   dp = d;
  622.   om = &OM[c][ya][0];
  624.   if (!addflag) {
  625.     for (j = ya; j < yb; j++) {
  626.       for (i = xa; i < xb; i++) {
  627. dp[i] = (((unsigned int)(sp[i] + sp2[i]+1))>>1)*om[i];
  628.       }
  629.       sp+= lx2;
  630.       sp2+= lx2;
  631.       dp+= lx;
  632.       om+= 8;
  633.     }
  634.   }
  635.   else {
  636.     for (j = ya; j < yb; j++) {
  637.       for (i = xa; i < xb; i++) {
  638. dp[i] += (((unsigned int)(sp[i] + sp2[i]+1))>>1)*om[i];
  639.       }
  640.       sp+= lx2;
  641.       sp2+= lx2;
  642.       dp+= lx;
  643.       om+= 8;
  644.     }
  645.   }
  647. }
  649. void CRecon::rec4(unsigned char *s, unsigned char *d, int lx, int lx2, int h)
  650. {
  651.   unsigned char *dp,*sp,*sp2;
  652.   int j;
  653.   unsigned int s1,s2,s3,s4;
  655.   sp = s;
  656.   sp2 = s+lx2;
  657.   dp = d;
  658.   for (j=0; j<h; j++)
  659.   {
  660.     s1=sp[0]; s3=sp2[0];
  661.     dp[0] = (unsigned int)(s1+(s2=sp[1])+s3+(s4=sp2[1])+2)>>2;
  662.     dp[1] = (unsigned int)(s2+(s1=sp[2])+s4+(s3=sp2[2])+2)>>2;
  663.     dp[2] = (unsigned int)(s1+(s2=sp[3])+s3+(s4=sp2[3])+2)>>2;
  664.     dp[3] = (unsigned int)(s2+(s1=sp[4])+s4+(s3=sp2[4])+2)>>2;
  665.     dp[4] = (unsigned int)(s1+(s2=sp[5])+s3+(s4=sp2[5])+2)>>2;
  666.     dp[5] = (unsigned int)(s2+(s1=sp[6])+s4+(s3=sp2[6])+2)>>2;
  667.     dp[6] = (unsigned int)(s1+(s2=sp[7])+s3+(s4=sp2[7])+2)>>2;
  668.     dp[7] = (unsigned int)(s2+(s1=sp[8])+s4+(s3=sp2[8])+2)>>2;
  669.     dp[8] = (unsigned int)(s1+(s2=sp[9])+s3+(s4=sp2[9])+2)>>2;
  670.     dp[9] = (unsigned int)(s2+(s1=sp[10])+s4+(s3=sp2[10])+2)>>2;
  671.     dp[10] = (unsigned int)(s1+(s2=sp[11])+s3+(s4=sp2[11])+2)>>2;
  672.     dp[11] = (unsigned int)(s2+(s1=sp[12])+s4+(s3=sp2[12])+2)>>2;
  673.     dp[12] = (unsigned int)(s1+(s2=sp[13])+s3+(s4=sp2[13])+2)>>2;
  674.     dp[13] = (unsigned int)(s2+(s1=sp[14])+s4+(s3=sp2[14])+2)>>2;
  675.     dp[14] = (unsigned int)(s1+(s2=sp[15])+s3+(s4=sp2[15])+2)>>2;
  676.     dp[15] = (unsigned int)(s2+sp[16]+s4+sp2[16]+2)>>2;
  677.     sp+= lx2;
  678.     sp2+= lx2;
  679.     dp+= lx;
  680.   }
  682. }
  684. void CRecon::rec4c(unsigned char *s, unsigned char *d, int lx, int lx2, int h)
  685. {
  686.  unsigned char *dp,*sp,*sp2;
  687.   int j;
  688.   unsigned int s1,s2,s3,s4;
  690.   sp = s;
  691.   sp2 = s+lx2;
  692.   dp = d;
  693.   for (j=0; j<h; j++)
  694.   {
  695.     s1=sp[0]; s3=sp2[0];
  696.     dp[0] = (unsigned int)(s1+(s2=sp[1])+s3+(s4=sp2[1])+2)>>2;
  697.     dp[1] = (unsigned int)(s2+(s1=sp[2])+s4+(s3=sp2[2])+2)>>2;
  698.     dp[2] = (unsigned int)(s1+(s2=sp[3])+s3+(s4=sp2[3])+2)>>2;
  699.     dp[3] = (unsigned int)(s2+(s1=sp[4])+s4+(s3=sp2[4])+2)>>2;
  700.     dp[4] = (unsigned int)(s1+(s2=sp[5])+s3+(s4=sp2[5])+2)>>2;
  701.     dp[5] = (unsigned int)(s2+(s1=sp[6])+s4+(s3=sp2[6])+2)>>2;
  702.     dp[6] = (unsigned int)(s1+(s2=sp[7])+s3+(s4=sp2[7])+2)>>2;
  703.     dp[7] = (unsigned int)(s2+sp[8]+s4+sp2[8]+2)>>2;
  704.     sp+= lx2;
  705.     sp2+= lx2;
  706.     dp+= lx;
  707.   }
  709. }
  711. void CRecon::rec4o(unsigned char *s, int *d, int lx, int lx2, int addflag, int c, int xa, int xb, int ya, int yb)
  712. {
  713.   int *dp,*om;
  714.   unsigned char *sp,*sp2;
  715.   int i,j;
  717.   sp = s;
  718.   sp2 = s+lx2;
  719.   dp = d;
  720.   om = &OM[c][ya][0];
  722.   if (!addflag) {
  723.     for (j = ya; j < yb; j++) {
  724.       for (i = xa; i < xb; i++) {
  725. dp[i] = (((unsigned int)(sp[i]+sp2[i]+sp[i+1]+sp2[i+1]+2))>>2)*om[i];
  726.       }
  727.       sp+= lx2;
  728.       sp2+= lx2;
  729.       dp+= lx;
  730.       om+= 8;
  731.     }
  732.   }
  733.   else {
  734.     for (j = ya; j < yb; j++) {
  735.       for (i = xa; i < xb; i++) {
  736. dp[i] += (((unsigned int)(sp[i]+sp2[i]+sp[i+1]+sp2[i+1]+2))>>2)*om[i];
  737.       }
  738.       sp+= lx2;
  739.       sp2+= lx2;
  740.       dp+= lx;
  741.       om += 8;
  742.     }
  743.   }
  745. }