📄 transform8x8.c
字号:
img->mpr[3+ioff][3+joff] =
img->mpr[2+ioff][5+joff] =
img->mpr[1+ioff][7+joff] = (P_E + P_G + 2*P_F + 2) >> 2;
img->mpr[5+ioff][1+joff] =
img->mpr[4+ioff][3+joff] =
img->mpr[3+ioff][5+joff] =
img->mpr[2+ioff][7+joff] = (P_F + P_H + 2*P_G + 2) >> 2;
img->mpr[6+ioff][1+joff] =
img->mpr[5+ioff][3+joff] =
img->mpr[4+ioff][5+joff] =
img->mpr[3+ioff][7+joff] = (P_G + P_I + 2*P_H + 2) >> 2;
img->mpr[7+ioff][1+joff] =
img->mpr[6+ioff][3+joff] =
img->mpr[5+ioff][5+joff] =
img->mpr[4+ioff][7+joff] = (P_H + P_J + 2*P_I + 2) >> 2;
img->mpr[7+ioff][3+joff] =
img->mpr[6+ioff][5+joff] =
img->mpr[5+ioff][7+joff] = (P_I + P_K + 2*P_J + 2) >> 2;
img->mpr[7+ioff][5+joff] =
img->mpr[6+ioff][7+joff] = (P_J + P_L + 2*P_K + 2) >> 2;
img->mpr[7+ioff][7+joff] = (P_K + P_M + 2*P_L + 2) >> 2;
break;
case DIAG_DOWN_RIGHT_PRED:
if ((!block_available_up)||(!block_available_left)||(!block_available_up_left))
printf ("warning: Intra_8x8_Diagonal_Down_Right prediction mode not allowed at mb %d\n",img->current_mb_nr);
// Mode DIAG_DOWN_RIGHT_PRED
img->mpr[0+ioff][7+joff] = (P_X + P_V + 2*(P_W) + 2) >> 2;
img->mpr[0+ioff][6+joff] =
img->mpr[1+ioff][7+joff] = (P_W + P_U + 2*(P_V) + 2) >> 2;
img->mpr[0+ioff][5+joff] =
img->mpr[1+ioff][6+joff] =
img->mpr[2+ioff][7+joff] = (P_V + P_T + 2*(P_U) + 2) >> 2;
img->mpr[0+ioff][4+joff] =
img->mpr[1+ioff][5+joff] =
img->mpr[2+ioff][6+joff] =
img->mpr[3+ioff][7+joff] = (P_U + P_S + 2*(P_T) + 2) >> 2;
img->mpr[0+ioff][3+joff] =
img->mpr[1+ioff][4+joff] =
img->mpr[2+ioff][5+joff] =
img->mpr[3+ioff][6+joff] =
img->mpr[4+ioff][7+joff] = (P_T + P_R + 2*(P_S) + 2) >> 2;
img->mpr[0+ioff][2+joff] =
img->mpr[1+ioff][3+joff] =
img->mpr[2+ioff][4+joff] =
img->mpr[3+ioff][5+joff] =
img->mpr[4+ioff][6+joff] =
img->mpr[5+ioff][7+joff] = (P_S + P_Q + 2*(P_R) + 2) >> 2;
img->mpr[0+ioff][1+joff] =
img->mpr[1+ioff][2+joff] =
img->mpr[2+ioff][3+joff] =
img->mpr[3+ioff][4+joff] =
img->mpr[4+ioff][5+joff] =
img->mpr[5+ioff][6+joff] =
img->mpr[6+ioff][7+joff] = (P_R + P_Z + 2*(P_Q) + 2) >> 2;
img->mpr[0+ioff][0+joff] =
img->mpr[1+ioff][1+joff] =
img->mpr[2+ioff][2+joff] =
img->mpr[3+ioff][3+joff] =
img->mpr[4+ioff][4+joff] =
img->mpr[5+ioff][5+joff] =
img->mpr[6+ioff][6+joff] =
img->mpr[7+ioff][7+joff] = (P_Q + P_A + 2*(P_Z) + 2) >> 2;
img->mpr[1+ioff][0+joff] =
img->mpr[2+ioff][1+joff] =
img->mpr[3+ioff][2+joff] =
img->mpr[4+ioff][3+joff] =
img->mpr[5+ioff][4+joff] =
img->mpr[6+ioff][5+joff] =
img->mpr[7+ioff][6+joff] = (P_Z + P_B + 2*(P_A) + 2) >> 2;
img->mpr[2+ioff][0+joff] =
img->mpr[3+ioff][1+joff] =
img->mpr[4+ioff][2+joff] =
img->mpr[5+ioff][3+joff] =
img->mpr[6+ioff][4+joff] =
img->mpr[7+ioff][5+joff] = (P_A + P_C + 2*(P_B) + 2) >> 2;
img->mpr[3+ioff][0+joff] =
img->mpr[4+ioff][1+joff] =
img->mpr[5+ioff][2+joff] =
img->mpr[6+ioff][3+joff] =
img->mpr[7+ioff][4+joff] = (P_B + P_D + 2*(P_C) + 2) >> 2;
img->mpr[4+ioff][0+joff] =
img->mpr[5+ioff][1+joff] =
img->mpr[6+ioff][2+joff] =
img->mpr[7+ioff][3+joff] = (P_C + P_E + 2*(P_D) + 2) >> 2;
img->mpr[5+ioff][0+joff] =
img->mpr[6+ioff][1+joff] =
img->mpr[7+ioff][2+joff] = (P_D + P_F + 2*(P_E) + 2) >> 2;
img->mpr[6+ioff][0+joff] =
img->mpr[7+ioff][1+joff] = (P_E + P_G + 2*(P_F) + 2) >> 2;
img->mpr[7+ioff][0+joff] = (P_F + P_H + 2*(P_G) + 2) >> 2;
break;
case VERT_RIGHT_PRED:/* diagonal prediction -22.5 deg to horizontal plane */
if ((!block_available_up)||(!block_available_left)||(!block_available_up_left))
printf ("warning: Intra_8x8_Vertical_Right prediction mode not allowed at mb %d\n",img->current_mb_nr);
img->mpr[0+ioff][0+joff] =
img->mpr[1+ioff][2+joff] =
img->mpr[2+ioff][4+joff] =
img->mpr[3+ioff][6+joff] = (P_Z + P_A + 1) >> 1;
img->mpr[1+ioff][0+joff] =
img->mpr[2+ioff][2+joff] =
img->mpr[3+ioff][4+joff] =
img->mpr[4+ioff][6+joff] = (P_A + P_B + 1) >> 1;
img->mpr[2+ioff][0+joff] =
img->mpr[3+ioff][2+joff] =
img->mpr[4+ioff][4+joff] =
img->mpr[5+ioff][6+joff] = (P_B + P_C + 1) >> 1;
img->mpr[3+ioff][0+joff] =
img->mpr[4+ioff][2+joff] =
img->mpr[5+ioff][4+joff] =
img->mpr[6+ioff][6+joff] = (P_C + P_D + 1) >> 1;
img->mpr[4+ioff][0+joff] =
img->mpr[5+ioff][2+joff] =
img->mpr[6+ioff][4+joff] =
img->mpr[7+ioff][6+joff] = (P_D + P_E + 1) >> 1;
img->mpr[5+ioff][0+joff] =
img->mpr[6+ioff][2+joff] =
img->mpr[7+ioff][4+joff] = (P_E + P_F + 1) >> 1;
img->mpr[6+ioff][0+joff] =
img->mpr[7+ioff][2+joff] = (P_F + P_G + 1) >> 1;
img->mpr[7+ioff][0+joff] = (P_G + P_H + 1) >> 1;
img->mpr[0+ioff][1+joff] =
img->mpr[1+ioff][3+joff] =
img->mpr[2+ioff][5+joff] =
img->mpr[3+ioff][7+joff] = (P_Q + P_A + 2*P_Z + 2) >> 2;
img->mpr[1+ioff][1+joff] =
img->mpr[2+ioff][3+joff] =
img->mpr[3+ioff][5+joff] =
img->mpr[4+ioff][7+joff] = (P_Z + P_B + 2*P_A + 2) >> 2;
img->mpr[2+ioff][1+joff] =
img->mpr[3+ioff][3+joff] =
img->mpr[4+ioff][5+joff] =
img->mpr[5+ioff][7+joff] = (P_A + P_C + 2*P_B + 2) >> 2;
img->mpr[3+ioff][1+joff] =
img->mpr[4+ioff][3+joff] =
img->mpr[5+ioff][5+joff] =
img->mpr[6+ioff][7+joff] = (P_B + P_D + 2*P_C + 2) >> 2;
img->mpr[4+ioff][1+joff] =
img->mpr[5+ioff][3+joff] =
img->mpr[6+ioff][5+joff] =
img->mpr[7+ioff][7+joff] = (P_C + P_E + 2*P_D + 2) >> 2;
img->mpr[5+ioff][1+joff] =
img->mpr[6+ioff][3+joff] =
img->mpr[7+ioff][5+joff] = (P_D + P_F + 2*P_E + 2) >> 2;
img->mpr[6+ioff][1+joff] =
img->mpr[7+ioff][3+joff] = (P_E + P_G + 2*P_F + 2) >> 2;
img->mpr[7+ioff][1+joff] = (P_F + P_H + 2*P_G + 2) >> 2;
img->mpr[0+ioff][2+joff] =
img->mpr[1+ioff][4+joff] =
img->mpr[2+ioff][6+joff] = (P_R + P_Z + 2*P_Q + 2) >> 2;
img->mpr[0+ioff][3+joff] =
img->mpr[1+ioff][5+joff] =
img->mpr[2+ioff][7+joff] = (P_S + P_Q + 2*P_R + 2) >> 2;
img->mpr[0+ioff][4+joff] =
img->mpr[1+ioff][6+joff] = (P_T + P_R + 2*P_S + 2) >> 2;
img->mpr[0+ioff][5+joff] =
img->mpr[1+ioff][7+joff] = (P_U + P_S + 2*P_T + 2) >> 2;
img->mpr[0+ioff][6+joff] = (P_V + P_T + 2*P_U + 2) >> 2;
img->mpr[0+ioff][7+joff] = (P_W + P_U + 2*P_V + 2) >> 2;
break;
case HOR_DOWN_PRED:/* diagonal prediction -22.5 deg to horizontal plane */
if ((!block_available_up)||(!block_available_left)||(!block_available_up_left))
printf ("warning: Intra_8x8_Horizontal_Down prediction mode not allowed at mb %d\n",img->current_mb_nr);
img->mpr[0+ioff][0+joff] =
img->mpr[2+ioff][1+joff] =
img->mpr[4+ioff][2+joff] =
img->mpr[6+ioff][3+joff] = (P_Q + P_Z + 1) >> 1;
img->mpr[0+ioff][1+joff] =
img->mpr[2+ioff][2+joff] =
img->mpr[4+ioff][3+joff] =
img->mpr[6+ioff][4+joff] = (P_R + P_Q + 1) >> 1;
img->mpr[0+ioff][2+joff] =
img->mpr[2+ioff][3+joff] =
img->mpr[4+ioff][4+joff] =
img->mpr[6+ioff][5+joff] = (P_S + P_R + 1) >> 1;
img->mpr[0+ioff][3+joff] =
img->mpr[2+ioff][4+joff] =
img->mpr[4+ioff][5+joff] =
img->mpr[6+ioff][6+joff] = (P_T + P_S + 1) >> 1;
img->mpr[0+ioff][4+joff] =
img->mpr[2+ioff][5+joff] =
img->mpr[4+ioff][6+joff] =
img->mpr[6+ioff][7+joff] = (P_U + P_T + 1) >> 1;
img->mpr[0+ioff][5+joff] =
img->mpr[2+ioff][6+joff] =
img->mpr[4+ioff][7+joff] = (P_V + P_U + 1) >> 1;
img->mpr[0+ioff][6+joff] =
img->mpr[2+ioff][7+joff] = (P_W + P_V + 1) >> 1;
img->mpr[0+ioff][7+joff] = (P_X + P_W + 1) >> 1;
img->mpr[1+ioff][0+joff] =
img->mpr[3+ioff][1+joff] =
img->mpr[5+ioff][2+joff] =
img->mpr[7+ioff][3+joff] = (P_Q + P_A + 2*P_Z + 2) >> 2;
img->mpr[1+ioff][1+joff] =
img->mpr[3+ioff][2+joff] =
img->mpr[5+ioff][3+joff] =
img->mpr[7+ioff][4+joff] = (P_Z + P_R + 2*P_Q + 2) >> 2;
img->mpr[1+ioff][2+joff] =
img->mpr[3+ioff][3+joff] =
img->mpr[5+ioff][4+joff] =
img->mpr[7+ioff][5+joff] = (P_Q + P_S + 2*P_R + 2) >> 2;
img->mpr[1+ioff][3+joff] =
img->mpr[3+ioff][4+joff] =
img->mpr[5+ioff][5+joff] =
img->mpr[7+ioff][6+joff] = (P_R + P_T + 2*P_S + 2) >> 2;
img->mpr[1+ioff][4+joff] =
img->mpr[3+ioff][5+joff] =
img->mpr[5+ioff][6+joff] =
img->mpr[7+ioff][7+joff] = (P_S + P_U + 2*P_T + 2) >> 2;
img->mpr[1+ioff][5+joff] =
img->mpr[3+ioff][6+joff] =
img->mpr[5+ioff][7+joff] = (P_T + P_V + 2*P_U + 2) >> 2;
img->mpr[1+ioff][6+joff] =
img->mpr[3+ioff][7+joff] = (P_U + P_W + 2*P_V + 2) >> 2;
img->mpr[1+ioff][7+joff] = (P_V + P_X + 2*P_W + 2) >> 2;
img->mpr[2+ioff][0+joff] =
img->mpr[4+ioff][1+joff] =
img->mpr[6+ioff][2+joff] = (P_Z + P_B + 2*P_A + 2) >> 2;
img->mpr[3+ioff][0+joff] =
img->mpr[5+ioff][1+joff] =
img->mpr[7+ioff][2+joff] = (P_A + P_C + 2*P_B + 2) >> 2;
img->mpr[4+ioff][0+joff] =
img->mpr[6+ioff][1+joff] = (P_B + P_D + 2*P_C + 2) >> 2;
img->mpr[5+ioff][0+joff] =
img->mpr[7+ioff][1+joff] = (P_C + P_E + 2*P_D + 2) >> 2;
img->mpr[6+ioff][0+joff] = (P_D + P_F + 2*P_E + 2) >> 2;
img->mpr[7+ioff][0+joff] = (P_E + P_G + 2*P_F + 2) >> 2;
break;
case HOR_UP_PRED:/* diagonal prediction -22.5 deg to horizontal plane */
if (!block_available_left)
printf ("warning: Intra_8x8_Horizontal_Up prediction mode not allowed at mb %d\n",img->current_mb_nr);
img->mpr[0+ioff][0+joff] = (P_Q + P_R + 1) >> 1;
img->mpr[0+ioff][1+joff] =
img->mpr[2+ioff][0+joff] = (P_R + P_S + 1) >> 1;
img->mpr[0+ioff][2+joff] =
img->mpr[2+ioff][1+joff] =
img->mpr[4+ioff][0+joff] = (P_S + P_T + 1) >> 1;
img->mpr[0+ioff][3+joff] =
img->mpr[2+ioff][2+joff] =
img->mpr[4+ioff][1+joff] =
img->mpr[6+ioff][0+joff] = (P_T + P_U + 1) >> 1;
img->mpr[0+ioff][4+joff] =
img->mpr[2+ioff][3+joff] =
img->mpr[4+ioff][2+joff] =
img->mpr[6+ioff][1+joff] = (P_U + P_V + 1) >> 1;
img->mpr[0+ioff][5+joff] =
img->mpr[2+ioff][4+joff] =
img->mpr[4+ioff][3+joff] =
img->mpr[6+ioff][2+joff] = (P_V + P_W + 1) >> 1;
img->mpr[0+ioff][6+joff] =
img->mpr[2+ioff][5+joff] =
img->mpr[4+ioff][4+joff] =
img->mpr[6+ioff][3+joff] = (P_W + P_X + 1) >> 1;
img->mpr[6+ioff][4+joff] =
img->mpr[7+ioff][4+joff] =
img->mpr[4+ioff][5+joff] =
img->mpr[5+ioff][5+joff] =
img->mpr[6+ioff][5+joff] =
img->mpr[7+ioff][5+joff] =
img->mpr[2+ioff][6+joff] =
img->mpr[3+ioff][6+joff] =
img->mpr[4+ioff][6+joff] =
img->mpr[5+ioff][6+joff] =
img->mpr[6+ioff][6+joff] =
img->mpr[7+ioff][6+joff] =
img->mpr[0+ioff][7+joff] =
img->mpr[1+ioff][7+joff] =
img->mpr[2+ioff][7+joff] =
img->mpr[3+ioff][7+joff] =
img->mpr[4+ioff][7+joff] =
img->mpr[5+ioff][7+joff] =
img->mpr[6+ioff][7+joff] =
img->mpr[7+ioff][7+joff] = P_X;
img->mpr[1+ioff][6+joff] =
img->mpr[3+ioff][5+joff] =
img->mpr[5+ioff][4+joff] =
img->mpr[7+ioff][3+joff] = (P_W + 3*P_X + 2) >> 2;
img->mpr[1+ioff][5+joff] =
img->mpr[3+ioff][4+joff] =
img->mpr[5+ioff][3+joff] =
img->mpr[7+ioff][2+joff] = (P_X + P_V + 2*P_W + 2) >> 2;
img->mpr[1+ioff][4+joff] =
img->mpr[3+ioff][3+joff] =
img->mpr[5+ioff][2+joff] =
img->mpr[7+ioff][1+joff] = (P_W + P_U + 2*P_V + 2) >> 2;
img->mpr[1+ioff][3+joff] =
img->mpr[3+ioff][2+joff] =
img->mpr[5+ioff][1+joff] =
img->mpr[7+ioff][0+joff] = (P_V + P_T + 2*P_U + 2) >> 2;
img->mpr[1+ioff][2+joff] =
img->mpr[3+ioff][1+joff] =
img->mpr[5+ioff][0+joff] = (P_U + P_S + 2*P_T + 2) >> 2;
img->mpr[1+ioff][1+joff] =
img->mpr[3+ioff][0+joff] = (P_T + P_R + 2*P_S + 2) >> 2;
img->mpr[1+ioff][0+joff] = (P_S + P_Q + 2*P_R + 2) >> 2;
break;
default:
printf("Error: illegal intra_4x4 prediction mode: %d\n",predmode);
return SEARCH_SYNC;
break;
}
return DECODING_OK;
}
/*!
*************************************************************************************
* \brief
* Prefiltering for Intra8x8 prediction
*************************************************************************************
*/
void LowPassForIntra8x8Pred(int *PredPel, int block_up_left, int block_up, int block_left)
{
int i;
int LoopArray[25];
for(i = 0; i < 25; i++)
LoopArray[i] = PredPel[i] ;
if(block_up)
{
if(block_up_left)
{
LoopArray[1] = ((&P_Z)[0] + ((&P_Z)[1]<<1) + (&P_Z)[2] + 2)>>2;
}
else
LoopArray[1] = ((&P_Z)[1] + ((&P_Z)[1]<<1) + (&P_Z)[2] + 2)>>2;
for(i = 2; i <16; i++)
{
LoopArray[i] = ((&P_Z)[i-1] + ((&P_Z)[i]<<1) + (&P_Z)[i+1] + 2)>>2;
}
LoopArray[16] = (P_P + (P_P<<1) + P_O + 2)>>2;
}
if(block_up_left)
{
if(block_up && block_left)
{
LoopArray[0] = (P_Q + (P_Z<<1) + P_A +2)>>2;
}
else
{
if(block_up)
LoopArray[0] = (P_Z + (P_Z<<1) + P_A +2)>>2;
else
if(block_left)
LoopArray[0] = (P_Z + (P_Z<<1) + P_Q +2)>>2;
}
}
if(block_left)
{
if(block_up_left)
LoopArray[17] = (P_Z + (P_Q<<1) + P_R + 2)>>2;
else
LoopArray[17] = (P_Q + (P_Q<<1) + P_R + 2)>>2;
for(i = 18; i <24; i++)
{
LoopArray[i] = ((&P_Z)[i-1] + ((&P_Z)[i]<<1) + (&P_Z)[i+1] + 2)>>2;
}
LoopArray[24] = (P_W + (P_X<<1) + P_X + 2)>>2;
}
for(i = 0; i < 25; i++)
PredPel[i] = LoopArray[i];
}
/*!
***********************************************************************
* \brief
* Inverse 8x8 transformation
***********************************************************************
*/
void itrans8x8(struct img_par *img, //!< image parameters
int ioff, //!< index to 4x4 block
int joff) //!<
{
int i,j;
int m6[8][8];
Boolean lossless_qpprime = (Boolean)((img->qp + img->bitdepth_luma_qp_scale)==0 && img->lossless_qpprime_flag==1);
int residue_transform_flag = img->residue_transform_flag;
for( i=0; i<8 && !lossless_qpprime; i++)
{
int a[8], b[8];
a[0] = img->m7[ioff + 0][joff + i] + img->m7[ioff + 4][joff + i];
a[4] = img->m7[ioff + 0][joff + i] - img->m7[ioff + 4][joff + i];
a[2] = (img->m7[ioff + 2][joff + i]>>1) - img->m7[ioff + 6][joff + i];
a[6] = img->m7[ioff + 2][joff + i] + (img->m7[ioff + 6][joff + i]>>1);
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -