📄 loopfilter.c
字号:
blkQ = ((yQ>>2)<<2) + (xQ>>2);
blkP = ((yP>>2)<<2) + (xP>>2);
if ((p->slice_type==SP_SLICE)||(p->slice_type==SI_SLICE) )
{
Strength[idx] = (edge == 0 && (((!p->MbaffFrameFlag && (p->structure==FRAME)) ||
(p->MbaffFrameFlag && !MbP->mb_field && !MbQ->mb_field)) ||
((p->MbaffFrameFlag || (p->structure != FRAME)) && !dir))) ? 4 : 3;
}
else
{
// Start with Strength=3. or Strength=4 for Mb-edge
Strength[idx] = (edge == 0 && (((!p->MbaffFrameFlag && (p->structure==FRAME)) ||
(p->MbaffFrameFlag && !MbP->mb_field && !MbQ->mb_field)) ||
((p->MbaffFrameFlag || (p->structure!=FRAME)) && !dir))) ? 4 : 3;
if( !(MbP->mb_type==I4MB || MbP->mb_type==I16MB || MbP->mb_type==I8MB || MbP->mb_type==IPCM)
&& !(MbQ->mb_type==I4MB || MbQ->mb_type==I16MB || MbQ->mb_type==I8MB || MbQ->mb_type==IPCM) )
{
if( ((MbQ->cbp_blk & (1 << blkQ )) != 0) || ((MbP->cbp_blk & (1 << blkP)) != 0) )
Strength[idx] = 2 ;
else
{ // if no coefs, but vector difference >= 1 set Strength=1
// if this is a mixed mode edge then one set of reference pictures will be frame and the
// other will be field
if (mixedModeEdgeFlag)
{
(Strength[idx] = 1);
}
else
{
get_mb_block_pos (MbQAddr, &mb_x, &mb_y);
blk_y = (mb_y<<2) + (blkQ >> 2) ;
blk_x = (mb_x<<2) + (blkQ & 3) ;
blk_y2 = pixP.pos_y >> 2;
blk_x2 = pixP.pos_x >> 2;
{
int64 ref_p0,ref_p1,ref_q0,ref_q1;
ref_p0 = list0_refIdxArr[blk_y] [blk_x] <0 ? INT64_MIN : list0_refPicIdArr[blk_y] [blk_x];
ref_q0 = list0_refIdxArr[blk_y2][blk_x2]<0 ? INT64_MIN : list0_refPicIdArr[blk_y2][blk_x2];
ref_p1 = list1_refIdxArr[blk_y] [blk_x] <0 ? INT64_MIN : list1_refPicIdArr[blk_y] [blk_x];
ref_q1 = list1_refIdxArr[blk_y2][blk_x2]<0 ? INT64_MIN : list1_refPicIdArr[blk_y2][blk_x2];
if ( ((ref_p0==ref_q0) && (ref_p1==ref_q1)) ||
((ref_p0==ref_q1) && (ref_p1==ref_q0)))
{
Strength[idx]=0;
// L0 and L1 reference pictures of p0 are different; q0 as well
if (ref_p0 != ref_p1)
{
// compare MV for the same reference picture
if (ref_p0==ref_q0)
{
Strength[idx] = (abs( list0_mv[blk_y][blk_x][0] - list0_mv[blk_y2][blk_x2][0]) >= 4) |
(abs( list0_mv[blk_y][blk_x][1] - list0_mv[blk_y2][blk_x2][1]) >= mvlimit) |
(abs( list1_mv[blk_y][blk_x][0] - list1_mv[blk_y2][blk_x2][0]) >= 4) |
(abs( list1_mv[blk_y][blk_x][1] - list1_mv[blk_y2][blk_x2][1]) >= mvlimit);
}
else
{
Strength[idx] = (abs( list0_mv[blk_y][blk_x][0] - list1_mv[blk_y2][blk_x2][0]) >= 4) |
(abs( list0_mv[blk_y][blk_x][1] - list1_mv[blk_y2][blk_x2][1]) >= mvlimit) |
(abs( list1_mv[blk_y][blk_x][0] - list0_mv[blk_y2][blk_x2][0]) >= 4) |
(abs( list1_mv[blk_y][blk_x][1] - list0_mv[blk_y2][blk_x2][1]) >= mvlimit);
}
}
else
{ // L0 and L1 reference pictures of p0 are the same; q0 as well
Strength[idx] = ((abs( list0_mv[blk_y][blk_x][0] - list0_mv[blk_y2][blk_x2][0]) >= 4) |
(abs( list0_mv[blk_y][blk_x][1] - list0_mv[blk_y2][blk_x2][1]) >= mvlimit ) |
(abs( list1_mv[blk_y][blk_x][0] - list1_mv[blk_y2][blk_x2][0]) >= 4) |
(abs( list1_mv[blk_y][blk_x][1] - list1_mv[blk_y2][blk_x2][1]) >= mvlimit))
&&
((abs( list0_mv[blk_y][blk_x][0] - list1_mv[blk_y2][blk_x2][0]) >= 4) |
(abs( list0_mv[blk_y][blk_x][1] - list1_mv[blk_y2][blk_x2][1]) >= mvlimit) |
(abs( list1_mv[blk_y][blk_x][0] - list0_mv[blk_y2][blk_x2][0]) >= 4) |
(abs( list1_mv[blk_y][blk_x][1] - list0_mv[blk_y2][blk_x2][1]) >= mvlimit));
}
}
else
{
Strength[idx] = 1;
}
}
}
}
}
}
}
}
#define CQPOF(qp, uv) (Clip3(0, 51, qp + p->chroma_qp_offset[uv]))
/*!
*****************************************************************************************
* \brief
* Filters one edge of 16 (luma) or 8 (chroma) pel
*****************************************************************************************
*/
void EdgeLoop(imgpel** Img, byte Strength[16],struct img_par *img, int MbQAddr, int AlphaC0Offset, int BetaOffset,
int dir, int edge, int width, int yuv, int uv, StorablePicture *p)
{
int pel, ap = 0, aq = 0, Strng ;
int incP, incQ;
int C0, c0, Delta, dif, AbsDelta ;
int L2 = 0, L1, L0, R0, R1, R2 = 0, RL0, L3, R3 ;
int Alpha = 0, Beta = 0 ;
byte* ClipTab = NULL;
int small_gap;
int indexA, indexB;
int PelNum;
int StrengthIdx;
imgpel *SrcPtrP, *SrcPtrQ;
int QP;
int xP, xQ, yP, yQ;
Macroblock *MbQ, *MbP;
PixelPos pixP, pixQ;
int bitdepth_scale;
int pelnum_cr[2][4] = {{0,8,16,16}, {0,8, 8,16}}; //[dir:0=vert, 1=hor.][yuv_format]
if (!yuv)
bitdepth_scale = 1<<(img->bitdepth_luma - 8);
else
bitdepth_scale = 1<<(img->bitdepth_chroma - 8);
PelNum = yuv ? pelnum_cr[dir][p->chroma_format_idc] : 16 ;
for( pel=0 ; pel<PelNum ; pel++ )
{
xQ = dir ? pel : edge << 2;
yQ = dir ? (edge < 4 ? edge << 2 : 1) : pel;
getNeighbour(MbQAddr, xQ, yQ, 1-yuv, &pixQ);
getNeighbour(MbQAddr, xQ - (1 - dir), yQ - dir, 1-yuv, &pixP);
xP = pixP.x;
yP = pixP.y;
MbQ = &(img->mb_data[MbQAddr]);
MbP = &(img->mb_data[pixP.mb_addr]);
fieldModeFilteringFlag = MbQ->mb_field || MbP->mb_field;
StrengthIdx = (yuv&&(PelNum==8)) ? ((MbQ->mb_field && !MbP->mb_field) ? pel<<1 :((pel>>1)<<2)+(pel%2)) : pel ;
if (pixP.available || (MbQ->LFDisableIdc== 0))
{
incQ = dir ? ((fieldModeFilteringFlag && !MbQ->mb_field) ? 2 * width : width) : 1;
incP = dir ? ((fieldModeFilteringFlag && !MbP->mb_field) ? 2 * width : width) : 1;
SrcPtrQ = &(Img[pixQ.pos_y][pixQ.pos_x]);
SrcPtrP = &(Img[pixP.pos_y][pixP.pos_x]);
// Average QP of the two blocks
QP = yuv ? (QP_SCALE_CR[CQPOF(MbP->qp,uv)] + QP_SCALE_CR[CQPOF(MbQ->qp,uv)] + 1) >> 1 : (MbP->qp + MbQ->qp + 1) >> 1;
indexA = IClip(0, MAX_QP, QP + AlphaC0Offset);
indexB = IClip(0, MAX_QP, QP + BetaOffset);
Alpha =ALPHA_TABLE[indexA] * bitdepth_scale;
Beta =BETA_TABLE[indexB] * bitdepth_scale;
ClipTab=CLIP_TAB[indexA];
L0 = SrcPtrP[0] ;
R0 = SrcPtrQ[0] ;
L1 = SrcPtrP[-incP] ;
R1 = SrcPtrQ[ incQ] ;
L2 = SrcPtrP[-incP*2] ;
R2 = SrcPtrQ[ incQ*2] ;
L3 = SrcPtrP[-incP*3] ;
R3 = SrcPtrQ[ incQ*3] ;
if( (Strng = Strength[StrengthIdx]) )
{
AbsDelta = abs( Delta = R0 - L0 ) ;
if( AbsDelta < Alpha )
{
C0 = ClipTab[ Strng ] * bitdepth_scale;
if( ((abs( R0 - R1) - Beta ) & (abs(L0 - L1) - Beta )) < 0 )
{
if( !yuv)
{
aq = (abs( R0 - R2) - Beta ) < 0 ;
ap = (abs( L0 - L2) - Beta ) < 0 ;
}
RL0 = L0 + R0 ;
if(Strng == 4 ) // INTRA strong filtering
{
if( yuv) // Chroma
{
SrcPtrQ[0] = ((R1 << 1) + R0 + L1 + 2) >> 2;
SrcPtrP[0] = ((L1 << 1) + L0 + R1 + 2) >> 2;
}
else // Luma
{
small_gap = (AbsDelta < ((Alpha >> 2) + 2));
aq &= small_gap;
ap &= small_gap;
SrcPtrQ[0] = aq ? ( L1 + ((R1 + RL0) << 1) + R2 + 4) >> 3 : ((R1 << 1) + R0 + L1 + 2) >> 2 ;
SrcPtrP[0] = ap ? ( R1 + ((L1 + RL0) << 1) + L2 + 4) >> 3 : ((L1 << 1) + L0 + R1 + 2) >> 2 ;
SrcPtrQ[ incQ] = aq ? ( R2 + R0 + R1 + L0 + 2) >> 2 : R1;
SrcPtrP[-incP] = ap ? ( L2 + L1 + L0 + R0 + 2) >> 2 : L1;
SrcPtrQ[ incQ*2] = aq ? (((R3 + R2) <<1) + R2 + R1 + RL0 + 4) >> 3 : R2;
SrcPtrP[-incP*2] = ap ? (((L3 + L2) <<1) + L2 + L1 + RL0 + 4) >> 3 : L2;
}
}
else // normal filtering
{
c0 = yuv? (C0+1):(C0 + ap + aq) ;
dif = IClip( -c0, c0, ( (Delta << 2) + (L1 - R1) + 4) >> 3 ) ;
if(!yuv)
{
SrcPtrP[0] = IClip(0, img->max_imgpel_value, L0 + dif) ;
SrcPtrQ[0] = IClip(0, img->max_imgpel_value, R0 - dif) ;
}
else
{
SrcPtrP[0] = IClip(0, img->max_imgpel_value_uv, L0 + dif) ;
SrcPtrQ[0] = IClip(0, img->max_imgpel_value_uv, R0 - dif) ;
}
if( !yuv )
{
if( ap )
SrcPtrP[-incP] += IClip( -C0, C0, ( L2 + ((RL0 + 1) >> 1) - (L1<<1)) >> 1 ) ;
if( aq )
SrcPtrQ[ incQ] += IClip( -C0, C0, ( R2 + ((RL0 + 1) >> 1) - (R1<<1)) >> 1 ) ;
} ;
} ;
} ;
} ;
} ;
} ;
}
}
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -