PDPC_0">位置决定的帧内预测组合 PDPC

1. 在 VVC 中，对于帧内预测的 Planar 模式、DC 模式和几种角度模式需要使用 PDPC (position dependent intra prediction combination) 方法进一步处理。 PDPC 用于 DC 模式、Planar 模式、小于等于水平模式(模式 18) 的角度模式、大于等于垂直模式(模式 50)且小于等于模式 80 的角度模式，即对模式 19～49 之间的角度模式不使用 PDPC 技术。如果当前块是 Bdpcm 模式或 MRL 索引大于 0，则不使用 PDPC。

2. 预测值是由参考像素和预测值再次加权求得；其中wL，wT是它们各自对应的基于水平或竖直距离的权重，如下公式：
   $pred(x^{,},y^{,})=Clip(0,(1<<BitDepth)-1,(wL.R_{-1,y}+wT.R_{x,-1}+(64-wL-wT).pred(x^{,},y^{,})+32)>>6)$

3. 如果 PDPC 用于 Planar、DC、水平、垂直模式则不需要额外的边界滤波。对于 Planar 和 DC 模式的 PDPC 操作是相同的。对于角度模式，如果是 HOR_IDX 或 VER_IDX 则左侧或上方的参考像素不会被使用。PDPC 的权值和缩放因子取决 于预测模式和块尺寸。只有块的宽和高都大于等于 4 时才可以使用 PDPC。

4. 图 6 是各种预测模式时 PDPC 参考像素定义( R_x,−1 和 R_−1,y )。预测值 pred (x’, y’)坐标为 (x’, y’)。对于对角线模式（模式 2 和模式 66），参考像素 R_x,−1 的坐标 x 为 x = x’ + y’ + 1， 参考像素 R_−1,y 的坐标 y 为 y = x’ + y’ + 1。对于其他角度模式 R_x,−1 和 R_−1,y 可能落地非整数像素位置，此时使用距其最近的整数像素值作为参考。
   

5. 修正权重wL，wT 的计算根据不同的模式有不同的计算方法，具体可以根据下表查找。这也可以从 VVenC 编码器源码中找到对应的实现，逻辑相对比较简单。
   

   • 偏移量 scale 的计算：
     • 如果角度模式50~66时，scale = Min( 2, Log2( nTbH ) − Floor( Log2( 3 * invAngle − 2 ) ) + 8 )
     • 如果角度模式为2~18时，scale = Min( 2, Log2( nTbW ) − Floor( Log2( 3 * invAngle − 2 ) ) + 8 )
     • 否则，scale = ( ( Log2( nTbW ) + Log2( nTbH ) − 2 ) >> 2 )

6. VVenC 编码器中是否使用 PDPC 技术在 IntraPrediction.cpp 文件中 initPredIntraParams 函数判断。
   

// Function for initialization of intra prediction parameters
void IntraPrediction::initPredIntraParams(const CodingUnit& cu, const CompArea area, const SPS& sps)
{const ComponentID compId = area.compID;const ChannelType chType = toChannelType(compId);const bool        useISP = NOT_INTRA_SUBPARTITIONS != cu.ispMode && isLuma( chType );const Size   cuSize    = Size( cu.blocks[compId].width, cu.blocks[compId].height );const Size   puSize    = Size( area.width, area.height );const Size&  blockSize = useISP ? cuSize : puSize;const int      dirMode = CU::getFinalIntraMode(cu, chType);const int     predMode = getWideAngle( blockSize.width, blockSize.height, dirMode );m_ipaParam.isModeVer            = predMode >= DIA_IDX;m_ipaParam.multiRefIndex        = isLuma (chType) ? cu.multiRefIdx : 0 ;m_ipaParam.refFilterFlag        = false;m_ipaParam.interpolationFlag    = false;m_ipaParam.applyPDPC            = (puSize.width >= MIN_TB_SIZEY && puSize.height >= MIN_TB_SIZEY) && m_ipaParam.multiRefIndex == 0;const int    intraPredAngleMode = (m_ipaParam.isModeVer) ? predMode - VER_IDX : -(predMode - HOR_IDX);int absAng = 0;if (dirMode > DC_IDX && dirMode < NUM_LUMA_MODE) // intraPredAngle for directional modes{static const int angTable[32]    = { 0,    1,    2,    3,    4,    6,     8,   10,   12,   14,   16,   18,   20,   23,   26,   29,   32,   35,   39,  45,  51,  57,  64,  73,  86, 102, 128, 171, 256, 341, 512, 1024 };static const int invAngTable[32] = {0,   16384, 8192, 5461, 4096, 2731, 2048, 1638, 1365, 1170, 1024, 910, 819, 712, 630, 565,512, 468,   420,  364,  321,  287,  256,  224,  191,  161,  128,  96,  64,  48,  32,  16};   // (512 * 32) / Angleconst int     absAngMode         = abs(intraPredAngleMode);const int     signAng            = intraPredAngleMode < 0 ? -1 : 1;absAng             = angTable  [absAngMode];m_ipaParam.absInvAngle           = invAngTable[absAngMode];m_ipaParam.intraPredAngle        = signAng * absAng;if (intraPredAngleMode < 0){m_ipaParam.applyPDPC = false;}else if (intraPredAngleMode > 0){const int sideSize = m_ipaParam.isModeVer ? puSize.height : puSize.width;const int maxScale = 2;m_ipaParam.angularScale = std::min(maxScale, floorLog2(sideSize) - (floorLog2(3 * m_ipaParam.absInvAngle - 2) - 8));m_ipaParam.applyPDPC &= m_ipaParam.angularScale >= 0;}}// high level conditions and DC intra predictionif( !isLuma( chType )|| useISP|| CU::isMIP( cu, chType ) //th remove this|| m_ipaParam.multiRefIndex|| DC_IDX == dirMode){}else if (cu.bdpcmM[chType]){m_ipaParam.refFilterFlag = false;}else if (dirMode == PLANAR_IDX) // Planar intra prediction{m_ipaParam.refFilterFlag = puSize.width * puSize.height > 32 ? true : false;}else if (!useISP)// HOR, VER and angular modes (MDIS){bool filterFlag = false;{const int diff = std::min<int>( abs( predMode - HOR_IDX ), abs( predMode - VER_IDX ) );const int log2Size = (Log2(puSize.width * puSize.height) >> 1);CHECK( log2Size >= MAX_INTRA_FILTER_DEPTHS, "Size not supported" );filterFlag = (diff > m_aucIntraFilter[log2Size]);}// Selelection of either ([1 2 1] / 4 ) refrence filter OR Gaussian 4-tap interpolation filterif (filterFlag){const bool isRefFilter       =  isIntegerSlope(absAng);CHECK( puSize.width * puSize.height <= 32, "DCT-IF interpolation filter is always used for 4x4, 4x8, and 8x4 luma CB" );m_ipaParam.refFilterFlag     =  isRefFilter;m_ipaParam.interpolationFlag = !isRefFilter;}}
}

7. DC、Planar 模式的 PDPC 处理在 IntraPrediction.cpp 文件中 IntraPredSampleFilter_Core函数中完成。
   

void  IntraPredSampleFilter_Core(PelBuf& dstBuf, const CPelBuf& pSrc)
{const int iWidth  = dstBuf.width;const int iHeight = dstBuf.height;const int scale = ((Log2(iWidth*iHeight) - 2) >> 2);CHECK(scale < 0 || scale > 31, "PDPC: scale < 0 || scale > 31");for (int y = 0; y < iHeight; y++){const int wT   = 32 >> std::min(31, ((y << 1) >> scale));const Pel left = pSrc.at(y + 1, 1);for (int x = 0; x < iWidth; x++){const int wL    = 32 >> std::min(31, ((x << 1) >> scale));const Pel top   = pSrc.at(x + 1, 0);const Pel val   = dstBuf.at(x, y);dstBuf.at(x, y) = val + ((wL * (left - val) + wT * (top - val) + 32) >> 6);}}
}

8. 水平、垂直模式的 PDPC 处理在 IntraPrediction.cpp 文件中 IntraHorVerPDPC_Core函数中完成。
   

void IntraHorVerPDPC_Core(Pel* pDsty,const int dstStride,Pel* refSide,const int width,const int height,int scale,const Pel* refMain, const ClpRng& clpRng)
{const Pel topLeft = refMain[0];for( int y = 0; y < height; y++ ){memcpy(pDsty,&refMain[1],width*sizeof(Pel));const Pel left    = refSide[1 + y];for (int x = 0; x < std::min(3 << scale, width); x++){const int wL  = 32 >> (2 * x >> scale);const Pel val = pDsty[x];pDsty[x]      = ClipPel(val + ((wL * (left - topLeft) + 32) >> 6), clpRng);}pDsty += dstStride;}
}

9. 其他角度模式（2～17、51～66）的 PDPC 处理在 IntraPrediction.cpp 文件中 IntraAnglePDPC_Core函数中完成。
   

void IntraAnglePDPC_Core(Pel* pDsty,const int dstStride,Pel* refSide,const int width,const int height,int scale,int invAngle)
{for (int y = 0; y<height; y++, pDsty += dstStride){int       invAngleSum = 256;for (int x = 0; x < std::min(3 << scale, width); x++){invAngleSum += invAngle;int wL   = 32 >> (2 * x >> scale);Pel left = refSide[y + (invAngleSum >> 9) + 1];pDsty[x] = pDsty[x] + ((wL * (left - pDsty[x]) + 32) >> 6);}}
}