JM8.6 色度帧内预测函数剖析（IntraChromaPrediction8x8）

最新推荐文章于 2021-12-11 22:59:32 发布
最新推荐文章于 2021-12-11 22:59:32 发布
阅读量1.1k 收藏 1
点赞数
文章标签： image up c
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/jack_incredible/article/details/6201202
版权

 

By jack_incredible

//位于macroblock.c中

void IntraChromaPrediction8x8 (int *mb_up, int *mb_left, int*mb_up_left)

{

 

  Macroblock *currMB = &img->mb_data[img->current_mb_nr];

  int     s, s0, s1, s2, s3, i, j, k;

  pel_t** image;

  int     block_x, block_y;

  int     mb_nr             = img->current_mb_nr;

  int     mb_available_up;

  int     mb_available_left[2];

  int     mb_available_up_left;

  int     ih,iv;

  int     ib,ic,iaa;

  int     uv;

  int     hline[8], vline[9];

  int     mode;

  int     best_mode = DC_PRED_8;         //just an initilaization here, should always be overwritten

  int     cost;

  int     min_cost;

  int     diff[16];

  PixelPos up;       //!< pixel position p(0,-1)

  PixelPos left[9];  //!< pixel positions p(-1, -1..8)

 

 

  for (i=0;i<9;i++)

  {

    getNeighbour(mb_nr, -1 ,  i-1 , 0, &left[i]);  //判断左上边像素加上左边8个像素（共9个像素）的有效性，若有效取得其位置信息。

  }

  

  getNeighbour(mb_nr, 0     ,  -1 , 0, &up);  //判断上边左起第一个像素的有效性

 

  mb_available_up                             = up.available;

  mb_available_up_left                        = left[0].available;

  mb_available_left[0] = mb_available_left[1] = left[1].available;

 

  if(input->UseConstrainedIntraPred)  // input->UseConstrainedIntraPred参数来自配置文件，一般情况下为0，即不执行下面块中的程序

  {

    mb_available_up      = up.available ? img->intra_block[up.mb_addr] : 0;

    for (i=1, mb_available_left[0]=1; i<5;i++)

      mb_available_left[0]  &= left[i].available ? img->intra_block[left[i].mb_addr]: 0;

    for (i=5, mb_available_left[1]=1; i<9;i++)

      mb_available_left[1]  &= left[i].available ? img->intra_block[left[i].mb_addr]: 0;

    mb_available_up_left = left[0].available ? img->intra_block[left[0].mb_addr]: 0;

  }

 

  if (mb_up)

    *mb_up = mb_available_up;

  if (mb_left)

    *mb_left = mb_available_left[0] && mb_available_left[1];

  if( mb_up_left )

    *mb_up_left = mb_available_up_left;

 

  // compute all chroma intra prediction modes for both U and V

  for (uv=0; uv<2; uv++)

  {

 image = enc_picture->imgUV[uv];

 

 // DC prediction

 for (block_y=0; block_y<8; block_y+=4)

 for (block_x=0; block_x<8; block_x+=4)

 {

 s=128;

 s0=s1=s2=s3=0;

 //===== get prediction value =====

 switch ((block_y>>1) + (block_x>>2))

 {

 case 0:  //===== TOP LEFT =====    //8x8块的左上边4x4部分进行预测

 if      (mb_available_up)       for (i=0;i<4;i++)  s0 += image[up.pos_y][up.pos_x + i];

 if      (mb_available_left[0])  for (i=1;i<5;i++)  s2 += image[left[i].pos_y][left[i].pos_x];

 if      (mb_available_up && mb_available_left[0])  s  = (s0+s2+4) >> 3;

 else if (mb_available_up)                          s  = (s0   +2) >> 2;

 else if (mb_available_left[0])                     s  = (s2   +2) >> 2;

 break;

 case 1: //===== TOP RIGHT =====   //8x8块的右上边4x4部分进行预测

 if      (mb_available_up)       for (i=4;i<8;i++)  s1 += image[up.pos_y][up.pos_x + i];

 else if (mb_available_left[0])  for (i=1;i<5;i++)  s2 += image[left[i].pos_y][left[i].pos_x];

 if      (mb_available_up)                          s  = (s1   +2) >> 2;

 else if (mb_available_left[0])                     s  = (s2   +2) >> 2;

 break;

 case 2: //===== BOTTOM LEFT =====    //8x8块的左下边4x4部分进行预测

 if      (mb_available_left[1])  for (i=5;i<9;i++)  s3 += image[left[i].pos_y][left[i].pos_x];

 else if (mb_available_up)       for (i=0;i<4;i++)  s0 += image[up.pos_y][up.pos_x + i];

 if      (mb_available_left[1])                     s  = (s3   +2) >> 2;

 else if (mb_available_up)                          s  = (s0   +2) >> 2;

 break;

 case 3: //===== BOTTOM RIGHT =====    //8x8块的右下边4x4部分进行预测

 if      (mb_available_up)       for (i=4;i<8;i++)  s1 += image[up.pos_y][up.pos_x + i];

 if      (mb_available_left[1])  for (i=5;i<9;i++)  s3 += image[left[i].pos_y][left[i].pos_x];

 if      (mb_available_up && mb_available_left[1])  s  = (s1+s3+4) >> 3;

 else if (mb_available_up)                          s  = (s1   +2) >> 2;

 else if (mb_available_left[1])                     s  = (s3   +2) >> 2;

 break;

 }           

 

 //===== prediction =====

 for (j=block_y; j<block_y+4; j++)

 for (i=block_x; i<block_x+4; i++)

 {

 img->mprr_c[uv][DC_PRED_8][i][j] = s;

 }

 }

 

 // vertical prediction

 if (mb_available_up)

 {

 for (i=0; i<8; i++)

 hline[i] = image[up.pos_y][up.pos_x + i];

 for (i=0; i<8; i++)

 for (j=0; j<8; j++)

 img->mprr_c[uv][VERT_PRED_8][i][j] = hline[i];

 }

 

 // horizontal prediction 

 if (mb_available_left[0] && mb_available_left[1])

 {

 for (i=1; i<9; i++)

 vline[i] = image[left[i].pos_y][left[i].pos_x];

 for (i=0; i<8; i++)

 for (j=0; j<8; j++)

 img->mprr_c[uv][HOR_PRED_8][i][j] = vline[j+1]; 

 }

 

 // plane prediction 

 if (mb_available_left[0] && mb_available_left[1] && mb_available_up && mb_available_up_left)

 {

 ih = 4*(hline[7] - image[left[0].pos_y][left[0].pos_x]);

 iv = 4*(vline[7+1] - image[left[0].pos_y][left[0].pos_x]);

 for (i=1;i<4;i++)

 {

 ih += i*(hline[3+i] - hline[3-i]);

 iv += i*(vline[3+i+1] - vline[3-i+1]);

 }

 ib=(17*ih+16)>>5;

 ic=(17*iv+16)>>5;

 

 iaa=16*(hline[7]+vline[7+1]);

 for (j=0; j<8; j++)

 for (i=0; i<8; i++)

 img->mprr_c[uv][PLANE_8][i][j]=max(0,min(255,(iaa+(i-3)*ib +(j-3)*ic + 16)/32));// store plane prediction

 }

  }  //end…for (uv=0; uv<2; uv++)

 

//这个地方貌似好像有点问题，觉得input->rdopt前不应该取！

if (!input->rdopt) // the rd-opt part does not work correctly (see encode_one_macroblock)

  {                       // since ipredmodes could be overwritten => encoder-decoder-mismatches

 // pick lowest cost prediction mode

 min_cost = 1<<20;

 for (i=0;i<8;i++)

 {

 getNeighbour(mb_nr, 0 ,  i , 0, &left[i]);  //获得当前块左边第一列的坐标

 }

 for (mode=DC_PRED_8; mode<=PLANE_8; mode++)

 {

 if ((mode==VERT_PRED_8 && !mb_available_up) ||

 (mode==HOR_PRED_8 && (!mb_available_left[0] || !mb_available_left[1])) ||

 (mode==PLANE_8 && (!mb_available_left[0] || !mb_available_left[1] || !mb_available_up || !mb_available_up_left)))

 continue;

 

 cost = 0;

 for (uv=0; uv<2; uv++)

 {

 image = imgUV_org[uv];

 for (block_y=0; block_y<8; block_y+=4)

 for (block_x=0; block_x<8; block_x+=4)

 {

 for (k=0,j=block_y; j<block_y+4; j++)

 for (i=block_x; i<block_x+4; i++,k++)

 {

 diff[k] = image[left[j].pos_y][left[j].pos_x+i] - img->mprr_c[uv][mode][i][j];  //求差分值

 }

 cost += SATD(diff, input->hadamard);  //SATD变换

 }

 }

 if (cost < min_cost)

 {

 best_mode = mode;

 min_cost = cost;

 }

 }

 

 currMB->c_ipred_mode = best_mode;

  }

 

}

 

 

jack_incredible
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
[JM8.6]帧内帧间模式选择
05-15
JM8.6中帧内与帧间模式的选择是其中非常重要的一部分,这一部分其实和码率控制中的代价函数计算是重叠在一起的, 在进行代价函数的计算过程中也就实现了模式的选择, 代价计算完毕, 最优模式也就选择完毕.
Intra Chroma Prediction
weixin_34146410的博客
12-07 148
  帧内预测依赖于当前宏块的相邻宏块，如果任何一个相邻宏块不可用，那么会直接影响到当前宏块的预测方式。 那么宏块怎么才谓之可用？ 满足以下几个条件的相邻宏块为不可用： 相邻宏块超出边界，即（x&lt;0 || x&gt;PicWidthInMbs），（y&lt;0 || y&gt;PicHeightInMbs） 相邻宏块与当前处理的宏块不在同一slice 如果强...
H.264编码器调试
m0_47121112的博客
12-11 557
1,实验内容 1，找出第一帧第二行第三个宏块的mode（编码模式）； 2，找出第二帧第二行第三个宏块的mode和MV（运动矢量）； 2，编码器调试 宏块（MB）类型（define.h） // Available MB modes typedef enum { PSKIP = 0, BSKIP_DIRECT = 0, P16x16 = 1, P16x8 = 2, P8x16 = 3, SMB8x8
色度分量帧内预测的代码
Lambda的专栏
01-07 844
JM模型中8x8块的色度分量帧内预测的代码 先保存下来 一边读一边注释 void IntraChromaPrediction8x8 (int *mb_up, int *mb_left, int*mb_up_left){  Macroblock *currMB = &img->mb_data[img->current_mb_nr];  int     s, s0, s1, s2, s3, i, j,
帧内预测模式提取
szu030606的专栏
10-12 1135
<br />if (input->rdopt)<br />    {<br />      int mb_available_up;<br />      int mb_available_left;<br />      int mb_available_up_left;<br />      <br />      min_rdcost = max_rdcost;<br />      <br />      // precompute all new chroma intra prediction m
JM73中encode_one_marcoblock()中关于帧内预测的函数流程===附：帧内预测的算法流程
timebomb的专栏
01-06 1551
  if (input->rdopt)       IntraChromaPrediction8x8() //8*8色度块的帧内预测，挑选最佳模式       Loop over 8*8 chroma intra prediction modes:              垂直、水平、平面预测模式不能进行的条件，跳出循环
JM8.6中帧内帧间模式的选择
08-29
在H.264视频编码标准中，JM8.6是一个特定版本的参考软件实现，它涉及到了帧内预测和帧间预测等关键技术。帧内预测是指在编码一个宏块时仅利用同一帧内的相邻像素作为参考，而帧间预测则涉及到对相邻帧的运动信息分析...
JM8.6源代码的解析
12-10
【JM8.6源代码解析】 JM8.6是一款基于Java技术开发的开源软件，其源代码提供了深入了解软件架构、设计模式以及编程实践的机会。在分析这个版本的源代码时，我们可以从中学习到许多关键的IT知识点。 1. **Java编程...
H.264测试软件JM8.6
04-07
通过分析JM8.6的源码，开发者和研究人员可以深入理解H.264的内部工作原理，包括宏块级的决策过程、熵编码的细节以及如何实现率失真优化。同时，JM8.6还提供了丰富的测试用例和参数设置，便于进行各种实验和性能评估...
JM8.6_一条余弦Cos.rar
06-11
JM8.6 源代码，设置编解码配置文件。并附详细配置步骤讲解和注释， 运行编码器，完成视频序列的H.264/AVC编码， 运行解码器，完成H.264/AVC视频流解码。 关注博主,下载YUVPlayer ...
8X8dzqmrj_LED_
09-30
Arduino LED8*8模块代码文件
Intra-prediction(4)之三种开销的计算模式
workspace
11-23 492
在Intra-prediction的最优模式选择中，我们需要计算各种模式的开销，而开销也有三种计算模式，分别为： SAD：Sum of Absolute Difference 绝对误差和，有 SAD=∑i=0,j=0Blocksize∣pij−pij′∣ SAD=\sum_{i=0,j=0}^{Blocksize}|p_{ij}-p&amp;amp;amp;amp;amp;amp;#x27;_{ij}| SAD=i=0,j=0∑Blo...
关于HM配置文件的TEncCfg类
小熊的专栏
01-07 2713
关于 class TEncCfg { protected: //==== File I/O ======== Int m_iFrameRate;//帧率 Int m_FrameSkip;//跳帧 Int m_iSourceWidth;//宽度 Int m_iSourceHeight;//高度 Int m_conf
JM、H.264中一些参数的设置和含义
学海无涯的专栏
03-16 4773
一·参数说明  这一节阐述的是encoder.cfg 中的参数对编码过程的影响  要注意的是encoder.cfg 中的参数跟input 结构体中的变量是一一对应的     StartFrame：从视频流的第几帧开始编码    FramesToBeEncoded：指明了除去 B帧后将要被编码的帧数  input->no_frames = FramesToBeEncoded    F
【H.264/AVC视频编解码技术详解】十七：帧内预测编码的预测实现方法
Workshop of Wenjie.Yin
12-16 7186
《H.264/AVC视频编解码技术详解》视频教程已经在“CSDN学院”上线，视频中详述了H.264的背景、标准协议和实现，并通过一个实战工程的形式对H.264的标准进行解析和实现，欢迎观看！“纸上得来终觉浅，绝知此事要躬行”，只有自己按照标准文档以代码的形式操作一遍，才能对视频压缩编码标准的思想和方法有足够深刻的理解和体会！链接地址：H.264/AVC视频编解码技术详解GitHub代码地址：点击这里
如何在JM8.6代码的编码端提取残差(为简便起见，仅以I4x4宏块为例)
认知 行动 坚持
11-20 9383
要在JM8.6代码中找出某一数据，不能瞎撞，不能用蛮力，而是要根据H.264的流程来找.比如，要在编码端找到像素残差，就可以按照下面的思路来找. （为简便起见，仅以I4x4宏块为例) 提前说一句：有时候对于同一个量，有好几个地方多可以提取，所以思路不要太局限，重要的是知道该量过去从何而来，将来又何去何从. 下面的思路是提取I4x4宏块中残差的思路.(啰嗦一句：如何要提取所有类型宏块
merry的共享空间 - 频道预览 - 抓虾
lukeamon的专栏
04-29 644
JM模型I帧帧内预测流程  I帧只存在帧内编码，没有帧间运动估计，不用参考其他的帧，所以I帧具有同步作用，.付出的代价就是效率稍差，不过也十分必要的。 I帧帧内编码分为亮度编码和色度编码，需要完成预测，计算RD代价，来判别宏块分块模式. I帧亮度度分块模式分为16X16,8X8,4X4三种模式，色度分块模式只有一种8X8模式，每种分块模式，又有不同的预测方式,在JM模型中，需要对这些模
x264 亮度信号8x8帧内预测模式
weixin_30516243的博客
09-20 224
1 该模式的8个预测方向与4x4帧内预测模式一样。 2该模式只有High profile及更高的Profile的才有可能使用，Baseline、Main Profile、Preset为ultrafast时不支持，命令行参数为--8x8dct与--no-8x8dct控制是否使用，x264内部参数结构x264_param_t的成员b_transform_8x8的值控制是否使用。 　　1 x264...
H.264 Intra Prediction
瑜媛的专栏
06-06 3229
对I宏块编码，不参考其所在的条带之外的任何数据。I宏块可以出现在任何类型的条带中，但是I条带中只能包含I宏块。对I宏块进行编码，只能使用帧内预测方式。对于一个包含亮度或者色度采样点的块，在采样点之间，以及块和临近块的采样点之间，存在着较高的相关性。因此，帧内预测采用临近块的采样点来预测块中的采样点的值。 1、Intra Prediction类型 对于I宏块，根据帧内预测块的尺寸，对亮度分量的预
H.264帧内预测编码matlab
最新发布
05-23
H.264是一种视频编码标准，其中的帧内预测是其中一种压缩技术。在MATLAB中实现H.264编码可以使用一些开源库和工具，例如JM、x264等。以下是一个示例代码，实现了H.264帧内预测编码： ```matlab % Read input video vidObj = VideoReader('input_video.mp4'); % Create output video object outputVideo = VideoWriter('output_video.mp4','MPEG-4'); outputVideo.FrameRate = vidObj.FrameRate; open(outputVideo); % Set H.264 encoder parameters params = struct('Profile','High','Level','4.0','BitRate',5000000,'FrameRate',vidObj.FrameRate); % Loop through each video frame while hasFrame(vidObj) % Read current frame frame = readFrame(vidObj); % Convert to YCbCr color space YCbCr = rgb2ycbcr(frame); % Split YCbCr into separate components Y = YCbCr(:,:,1); Cb = YCbCr(:,:,2); Cr = YCbCr(:,:,3); % Perform frame intra prediction on Y component Y_pred = intra_prediction(Y); % Combine YCbCr components YCbCr_pred = cat(3,Y_pred,Cb,Cr); % Convert back to RGB color space frame_pred = ycbcr2rgb(YCbCr_pred); % Write encoded frame to output video writeVideo(outputVideo,frame_pred); end % Close output video object close(outputVideo); function Y_pred = intra_prediction(Y) % Block size (8x8) block_size = 8; % Pad input image [rows,cols] = size(Y); Y_pad = padarray(Y,[block_size-1 block_size-1],'replicate','post'); % Perform intra prediction on each block for r = 1:block_size:rows for c = 1:block_size:cols % Get current block block = Y_pad(r:r+block_size-1,c:c+block_size-1); % Calculate prediction mode mode = calculate_mode(block); % Perform prediction switch mode case 0 % DC mode pred = mean(block(:)); case 1 % Horizontal mode pred = repmat(block(1,:),[block_size 1]); case 2 % Vertical mode pred = repmat(block(:,1),[1 block_size]); case 3 % Diagonal mode pred = diag(block); pred = repmat(pred(:),[1 block_size]) + repmat(pred(:)',[block_size 1]); pred = pred/(block_size+1); otherwise error('Invalid prediction mode'); end % Subtract prediction from current block Y_pred(r:r+block_size-1,c:c+block_size-1) = block - pred; end end end function mode = calculate_mode(block) % Calculate sum of absolute differences (SAD) for each mode SAD = zeros(4,1); SAD(1) = sum(abs(block(:)-mean(block(:)))); SAD(2) = sum(abs(block(1,:)-mean(block(:)))); SAD(3) = sum(abs(block(:,1)-mean(block(:)))); SAD(4) = sum(abs(diag(block)-mean(block(:)))); % Find mode with minimum SAD [~,mode] = min(SAD); mode = mode-1; end ``` 这段代码实现了一个简单的帧内预测编码器，基于直接模式选择方法（mode selection），并且使用直流、水平、垂直和对角线预测模式进行预测。在实际应用中，需要使用更复杂的编码器和预测算法，以达到更好的压缩效果。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值