3bu.cpp

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

 for(int i_1=0;i_1<16;i_1++)
	 for(int j_1=0;j_1<16;j_1++)
	 {
		 old_block[i_1][j_1]=old_frame[16*m-d+i_1][16*n+j_1];
	 }
 temp[1]=mad(cur_block,old_block);

 for(int i_2=0;i_2<16;i_2++)
	 for(int j_2=0;j_2<16;j_2++)
	 {
		 old_block[i_2][j_2]=old_frame[16*m-d+i_2][16*n+d+j_2];
	 }
 temp[2]=mad(cur_block,old_block);

 for(int i_3=0;i_3<16;i_3++)
	 for(int j_3=0;j_3<16;j_3++)
	 {
		 old_block[i_3][j_3]=old_frame[16*m+i_3][16*n-d+j_3];
	 }
 temp[3]=mad(cur_block,old_block);

 for(int i_4=0;i_4<16;i_4++)
	 for(int j_4=0;j_4<16;j_4++)
	 {
		 old_block[i_4][j_4]=old_frame[16*m+i_4][16*n+j_4];
	 }
 temp[4]=mad(cur_block,old_block);

 for(int i_5=0;i_5<16;i_5++)
	 for(int j_5=0;j_5<16;j_5++)
	 {
		 old_block[i_5][j_5]=old_frame[16*m+i_5][16*n+d+j_5];
	 }
 temp[5]=mad(cur_block,old_block);

 for(int i_6=0;i_6<16;i_6++)
	 for(int j_6=0;j_6<16;j_6++)
	 {
		 old_block[i_6][j_6]=old_frame[16*m+d+i_6][16*n-d+j_6];
	 }
 temp[6]=mad(cur_block,old_block);

 for(int i_7=0;i_7<16;i_7++)
	 for(int j_7=0;j_7<16;j_7++)
	 {
		 old_block[i_7][j_7]=old_frame[16*m+d+i_7][16*n+j_7];
	 }
 temp[7]=mad(cur_block,old_block);

 for(int i_8=0;i_8<16;i_8++)
	 for(int j_8=0;j_8<16;j_8++)
	 {
		 old_block[i_8][j_8]=old_frame[16*m+d+i_8][16*n+d+j_8];
	 }
 temp[8]=mad(cur_block,old_block);


 int temp1[9];
 for(int i=0;i<9;i++)
 {	 temp1[i]=temp[i];}	 
 int temp_1=0;
 for(int i2=0;i2<9;i2++)
 {	 for(int j2=i2+1;j2<9;j2++)
	 {		 if(temp1[i2]>temp1[j2])
		    { temp_1=temp1[i2];
	         temp1[i2]=temp1[j2];
		     temp1[j2]=temp_1;
		     }
	 }
 }//从小到大排序
  for(int i3=0;i3<9;i3++)
 {	 if(temp1[0]==temp[i3])
	 { x=i3;}
  }
  return x;

}

int second_search()
{
	switch(x)
	{
	case 0:x2=16*m-d;y2=16*n-d;break;
	case 1:x2=16*m-d;y2=16*n;break;
	case 2:x2=16*m-d;y2=16*n+d;break;
	case 3:x2=16*m;y2=16*n-d;break;
	case 4:x2=16*m;y2=16*n;break;
	case 5:x2=16*m;y2=16*n+d;break;
	case 6:x2=16*m+d;y2=16*n-d;break;
	case 7:x2=16*m+d;y2=16*n;break;
	case 8:x2=16*m+d;y2=16*n+d;break;
	}

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

 for(int ii_1=0;ii_1<16;ii_1++)
	 for(int jj_1=0;jj_1<16;jj_1++)
	 {
		 old_block[ii_1][jj_1]=old_frame[x2-d/2+ii_1][y2+jj_1];
	 }
 temp[1]=mad(cur_block,old_block);

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

 for(int ii_3=0;ii_3<16;ii_3++)
	 for(int jj_3=0;jj_3<16;jj_3++)
	 {
		 old_block[ii_3][jj_3]=old_frame[x2+ii_3][y2-d/2+jj_3];
	 }
 temp[3]=mad(cur_block,old_block);

 for(int ii_4=0;ii_4<16;ii_4++)
	 for(int jj_4=0;jj_4<16;jj_4++)
	 {
		 old_block[ii_4][jj_4]=old_frame[x2+ii_4][y2+jj_4];
	 }
 temp[4]=mad(cur_block,old_block);

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

 for(int ii_6=0;ii_6<16;ii_6++)
	 for(int jj_6=0;jj_6<16;jj_6++)
	 {