3bu.cpp

预测编码中的运动估计3步搜索法 本程序用C语言编程实现了该算法 大家可以参考参考

		 old_block[ii_6][jj_6]=old_frame[x2+d/2+ii_6][y2-d/2+jj_6];
	 }
 temp[6]=mad(cur_block,old_block);

 for(int ii_7=0;ii_7<16;ii_7++)
	 for(int jj_7=0;jj_7<16;jj_7++)
	 {
		 old_block[ii_7][jj_7]=old_frame[x2+d/2+ii_7][y2+jj_7];
	 }
 temp[7]=mad(cur_block,old_block);

 for(int ii_8=0;ii_8<16;ii_8++)
	 for(int jj_8=0;jj_8<16;jj_8++)
	 {
		 old_block[ii_8][jj_8]=old_frame[x2+d/2+ii_8][y2+d/2+jj_8];
	 }
 temp[8]=mad(cur_block,old_block);

 int temp2[9];
 for(int i=0;i<9;i++)
 {	 temp2[i]=temp[i];}	 
 int temp_2=0;
 for(int i1=0;i1<9;i1++)
 {	 for(int j1=i1+1;j1<9;j1++)
	 {		 if(temp2[i1]>temp2[j1])
		    { temp_2=temp2[i1];
	         temp2[i1]=temp2[j1];
		     temp2[j1]=temp_2;
		     }
	 }
 }//从小到大排序
  for(int i2=0;i2<9;i2++)
 {	 if(temp2[0]==temp[i2])
	 { y=i2;}
  }
  return y;
}

int third_search()
{
	switch(y)
	{	case 0:x3=x2-d/2;y3=y2-d/2;break;
		case 1:x3=x2-d/2;y3=y2;break;
		case 2:x3=x2-d/2;y3=y2+d/2;break;
		case 3:x3=x2;y3=y2-d/2;break;
		case 4:x3=x2;y3=y2;break;
        case 5:x3=x2;y3=y2+d/2;break;
		case 6:x3=x2+d/2;y3=y2-d/2;break;
		case 7:x3=x2+d/2;y3=y2;break;
		case 8:x3=x2+d/2;y3=y2+d/2;break;
	}

for(int iii_0=0;iii_0<16;iii_0++)
	 for(int jjj_0=0;jjj_0<16;jjj_0++)
	 {
		 old_block[iii_0][jjj_0]=old_frame[x3-d/4+iii_0][y3-d/4+jjj_0];
	 }
 temp[0]=mad(cur_block,old_block);

 for(int iii_1=0;iii_1<16;iii_1++)
	 for(int jjj_1=0;jjj_1<16;jjj_1++)
	 {
		 old_block[iii_1][jjj_1]=old_frame[x3-d/4+iii_1][y3+jjj_1];
	 }
 temp[1]=mad(cur_block,old_block);

 for(int iii_2=0;iii_2<16;iii_2++)
	 for(int jjj_2=0;jjj_2<16;jjj_2++)
	 {
		 old_block[iii_2][jjj_2]=old_frame[x3-d/4+iii_2][y3+d/4+jjj_2];
	 }
 temp[2]=mad(cur_block,old_block);

 for(int iii_3=0;iii_3<16;iii_3++)
	 for(int jjj_3=0;jjj_3<16;jjj_3++)
	 {
		 old_block[iii_3][jjj_3]=old_frame[x3+iii_3][y3-d/4+jjj_3];
	 }
 temp[3]=mad(cur_block,old_block);

 for(int iii_4=0;iii_4<16;iii_4++)
	 for(int jjj_4=0;jjj_4<16;jjj_4++)
	 {
		 old_block[iii_4][jjj_4]=old_frame[x3+iii_4][y3+jjj_4];
	 }
 temp[4]=mad(cur_block,old_block);

 for(int iii_5=0;iii_5<16;iii_5++)
	 for(int jjj_5=0;jjj_5<16;jjj_5++)
	 {
		 old_block[iii_5][jjj_5]=old_frame[x3+iii_5][y3+d/4+jjj_5];
	 }
 temp[5]=mad(cur_block,old_block);

 for(int iii_6=0;iii_6<16;iii_6++)
	 for(int jjj_6=0;jjj_6<16;jjj_6++)
	 {
		 old_block[iii_6][jjj_6]=old_frame[x3+d/4+iii_6][y3-d/4+jjj_6];
	 }
 temp[6]=mad(cur_block,old_block);

 for(int iii_7=0;iii_7<16;iii_7++)
	 for(int jjj_7=0;jjj_7<16;jjj_7++)
	 {
		 old_block[iii_7][jjj_7]=old_frame[x3+d/4+iii_7][y3+jjj_7];
	 }
 temp[7]=mad(cur_block,old_block);

 for(int iii_8=0;iii_8<16;iii_8++)
	 for(int jjj_8=0;jjj_8<16;jjj_8++)
	 {
		 old_block[iii_8][jjj_8]=old_frame[x3+d/4+iii_8][y3+d/4+jjj_8];
	 }
 temp[8]=mad(cur_block,old_block);

  int temp3[9];
 for(int i=0;i<9;i++)
 {	 temp3[i]=temp[i];}	 
 int temp_3=0;
 for(int i1=0;i1<9;i1++)
 {	 for(int j1=i1+1;j1<9;j1++)
	 {		 if(temp3[i1]>temp3[j1])
		    { temp_3=temp3[i1];
	         temp3[i1]=temp3[j1];
		     temp3[j1]=temp_3;
		     }
	 }
 }//从小到大排序
  for(int i2=0;i2<9;i2++)
 {	 if(temp3[0]==temp[i2])
	 { z=i2;}
  }
  return z;
}