?? mc_prediction.c
字號:
cur_lineY_m1 = cur_lineY_m2 + shift_x;
cur_lineY = cur_lineY_m1 + shift_x;
cur_lineY_p1 = cur_lineY + shift_x;
cur_lineY_p2 = cur_lineY_p1 + shift_x;
cur_lineY_p3 = cur_lineY_p2 + shift_x;
for (i = 0; i < hor_block_size; i++)
{
result = (*(cur_lineY_m2++) + *(cur_lineY_p3++)) * COEF[0];
result += (*(cur_lineY_m1++) + *(cur_lineY_p2++)) * COEF[1];
result += (*(cur_lineY ++) + *(cur_lineY_p1++)) * COEF[2];
block[j][i] = (block[j][i] + iClip1(img->max_imgpel_value, ((result + 16) >> 5)) +1 ) >>1;
}
cur_lineY_m2 = cur_lineY_m1 - hor_block_size ;
}
}
}
}
else // unsafe positions
{
if (dx == 0 && dy == 0)
{ /* fullpel position */
for (j = 0; j < ver_block_size; j++)
{
cur_lineY = cur_imgY[iClip3(0, maxold_y, y_pos + j)];
for (i = 0; i < hor_block_size; i++)
{
block[j][i] = cur_lineY[iClip3(0, maxold_x, x_pos + i )];
}
}
}
else
{ /* other positions */
if (dy == 0)
{ /* No vertical interpolation */
tmp_pos = x_pos - 2;
for (i = 0; i < hor_block_size; i++)
{
ipos_m2 = iClip3(0, maxold_x, tmp_pos++);
ipos_m1 = iClip3(0, maxold_x, tmp_pos++);
ipos = iClip3(0, maxold_x, tmp_pos++);
ipos_p1 = iClip3(0, maxold_x, tmp_pos++);
ipos_p2 = iClip3(0, maxold_x, tmp_pos++);
ipos_p3 = iClip3(0, maxold_x, tmp_pos );
tmp_pos -= 4;
for (j = 0; j < ver_block_size; j++)
{
cur_lineY = cur_imgY[iClip3(0,maxold_y,y_pos+j)];
result = (cur_lineY[ipos_m2] + cur_lineY[ipos_p3]) * COEF[0];
result += (cur_lineY[ipos_m1] + cur_lineY[ipos_p2]) * COEF[1];
result += (cur_lineY[ipos ] + cur_lineY[ipos_p1]) * COEF[2];
block[j][i] = iClip1(img->max_imgpel_value, ((result + 16)>>5));
}
}
if ((dx&1) == 1)
{
for (j = 0; j < ver_block_size; j++)
{
cur_lineY = cur_imgY[iClip3(0,maxold_y,y_pos+j)];
for (i = 0; i < hor_block_size; i++)
{
block[j][i] = (block[j][i] + cur_lineY[iClip3(0, maxold_x, x_pos + i + (dx >> 1))] + 1 )>>1;
}
}
}
}
else if (dx == 0)
{ /* No horizontal interpolation */
tmp_pos = y_pos - 2;
for (j = 0; j < ver_block_size; j++)
{
jpos_m2 = iClip3(0, maxold_y, tmp_pos++);
jpos_m1 = iClip3(0, maxold_y, tmp_pos++);
jpos = iClip3(0, maxold_y, tmp_pos++);
jpos_p1 = iClip3(0, maxold_y, tmp_pos++);
jpos_p2 = iClip3(0, maxold_y, tmp_pos++);
jpos_p3 = iClip3(0, maxold_y, tmp_pos );
tmp_pos -= 4;
for (i = 0; i < hor_block_size; i++)
{
pres_x = iClip3(0,maxold_x,x_pos+i);
result = (cur_imgY[jpos_m2][pres_x] + cur_imgY[jpos_p3][pres_x]) * COEF[0];
result += (cur_imgY[jpos_m1][pres_x] + cur_imgY[jpos_p2][pres_x]) * COEF[1];
result += (cur_imgY[jpos ][pres_x] + cur_imgY[jpos_p1][pres_x]) * COEF[2];
block[j][i] = iClip1(img->max_imgpel_value, ((result+16)>>5));
}
}
if ((dy&1) == 1)
{
for (j = 0; j < ver_block_size; j++)
{
cur_lineY = cur_imgY[iClip3(0,maxold_y,y_pos+j+(dy>>1))];
for (i = 0; i < hor_block_size; i++)
{
block[j][i] = (block[j][i] + cur_lineY[iClip3(0, maxold_x, x_pos + i)] + 1 )>>1;
}
}
}
}
else if (dx == 2)
{ /* Vertical & horizontal interpolation */
tmp_pos = x_pos - 2;
for (i = 0; i < hor_block_size; i++)
{
ipos_m2 = iClip3(0, maxold_x, tmp_pos++);
ipos_m1 = iClip3(0, maxold_x, tmp_pos++);
ipos = iClip3(0, maxold_x, tmp_pos++);
ipos_p1 = iClip3(0, maxold_x, tmp_pos++);
ipos_p2 = iClip3(0, maxold_x, tmp_pos++);
ipos_p3 = iClip3(0, maxold_x, tmp_pos );
tmp_pos -= 4;
for (j = 0; j < ver_block_size + 5; j++)
{
cur_lineY = cur_imgY[iClip3(0,maxold_y,y_pos + j - 2)];
tmp_res[j][i] = (cur_lineY[ipos_m2] + cur_lineY[ipos_p3]) * COEF[0];
tmp_res[j][i] += (cur_lineY[ipos_m1] + cur_lineY[ipos_p2]) * COEF[1];
tmp_res[j][i] += (cur_lineY[ipos ] + cur_lineY[ipos_p1]) * COEF[2];
}
}
for (j = 0; j < ver_block_size; j++)
{
jpos_m2 = j ;
jpos_m1 = j + 1;
jpos = j + 2;
jpos_p1 = j + 3;
jpos_p2 = j + 4;
jpos_p3 = j + 5;
for (i = 0; i < hor_block_size; i++)
{
result = (tmp_res[jpos_m2][i] + tmp_res[jpos_p3][i]) * COEF[0];
result += (tmp_res[jpos_m1][i] + tmp_res[jpos_p2][i]) * COEF[1];
result += (tmp_res[jpos ][i] + tmp_res[jpos_p1][i]) * COEF[2];
block[j][i] = iClip1(img->max_imgpel_value, ((result+512)>>10));
}
}
if ((dy&1) == 1)
{
for (j = 0; j < ver_block_size; j++)
{
pres_y = j+2+(dy>>1);
for (i = 0; i < hor_block_size; i++)
{
block[j][i] = (block[j][i] + iClip1(img->max_imgpel_value, ((tmp_res[pres_y][i] + 16) >> 5)) + 1 )>>1;
}
}
}
}
else if (dy == 2)
{ /* Horizontal & vertical interpolation */
tmp_pos = y_pos - 2;
for (j = 0; j < ver_block_size; j++)
{
jpos_m2 = iClip3(0, maxold_y, tmp_pos++);
jpos_m1 = iClip3(0, maxold_y, tmp_pos++);
jpos = iClip3(0, maxold_y, tmp_pos++);
jpos_p1 = iClip3(0, maxold_y, tmp_pos++);
jpos_p2 = iClip3(0, maxold_y, tmp_pos++);
jpos_p3 = iClip3(0, maxold_y, tmp_pos );
tmp_pos -= 4;
for (i = 0; i < hor_block_size + 5; i++)
{
pres_x = iClip3(0,maxold_x, x_pos + i - 2);
tmp_res[j][i] = (cur_imgY[jpos_m2][pres_x] + cur_imgY[jpos_p3][pres_x])*COEF[0];
tmp_res[j][i] += (cur_imgY[jpos_m1][pres_x] + cur_imgY[jpos_p2][pres_x])*COEF[1];
tmp_res[j][i] += (cur_imgY[jpos ][pres_x] + cur_imgY[jpos_p1][pres_x])*COEF[2];
}
}
for (j = 0; j < ver_block_size; j++)
{
for (i = 0; i < hor_block_size; i++)
{
result = (tmp_res[j][i ] + tmp_res[j][i + 5]) * COEF[0];
result += (tmp_res[j][i + 1] + tmp_res[j][i + 4]) * COEF[1];
result += (tmp_res[j][i + 2] + tmp_res[j][i + 3]) * COEF[2];
block[j][i] = iClip1(img->max_imgpel_value, ((result+512)>>10));
}
}
if ((dx&1) == 1)
{
for (j = 0; j < ver_block_size; j++)
{
for (i = 0; i < hor_block_size; i++)
{
block[j][i] = (block[j][i] + iClip1(img->max_imgpel_value, ((tmp_res[j][i + 2 + (dx>>1)] + 16)>>5))+1)>>1;
}
}
}
}
else
{ /* Diagonal interpolation */
tmp_pos = x_pos - 2;
for (i = 0; i < hor_block_size; i++)
{
ipos_m2 = iClip3(0, maxold_x, tmp_pos++);
ipos_m1 = iClip3(0, maxold_x, tmp_pos++);
ipos = iClip3(0, maxold_x, tmp_pos++);
ipos_p1 = iClip3(0, maxold_x, tmp_pos++);
ipos_p2 = iClip3(0, maxold_x, tmp_pos++);
ipos_p3 = iClip3(0, maxold_x, tmp_pos );
tmp_pos -= 4;
for (j = 0; j < ver_block_size; j++)
{
cur_lineY = cur_imgY[iClip3(0,maxold_y,(dy == 1 ? y_pos+j : y_pos+j+1))];
result = (cur_lineY[ipos_m2] + cur_lineY[ipos_p3]) * COEF[0];
result += (cur_lineY[ipos_m1] + cur_lineY[ipos_p2]) * COEF[1];
result += (cur_lineY[ipos ] + cur_lineY[ipos_p1]) * COEF[2];
block[j][i] = iClip1(img->max_imgpel_value, ((result+16)>>5));
}
}
tmp_pos = y_pos - 2;
for (j = 0; j < ver_block_size; j++)
{
jpos_m2 = iClip3(0, maxold_y, tmp_pos++);
jpos_m1 = iClip3(0, maxold_y, tmp_pos++);
jpos = iClip3(0, maxold_y, tmp_pos++);
jpos_p1 = iClip3(0, maxold_y, tmp_pos++);
jpos_p2 = iClip3(0, maxold_y, tmp_pos++);
jpos_p3 = iClip3(0, maxold_y, tmp_pos );
tmp_pos -= 4;
for (i = 0; i < hor_block_size; i++)
{
pres_x = dx == 1 ? x_pos+i : x_pos+i+1;
pres_x = iClip3(0, maxold_x, pres_x);
result = (cur_imgY[jpos_m2][pres_x] + cur_imgY[jpos_p3][pres_x]) * COEF[0];
result += (cur_imgY[jpos_m1][pres_x] + cur_imgY[jpos_p2][pres_x]) * COEF[1];
result += (cur_imgY[jpos ][pres_x] + cur_imgY[jpos_p1][pres_x]) * COEF[2];
block[j][i] = (block[j][i] + iClip1(img->max_imgpel_value, ((result+16)>>5)) +1 ) >>1;
}
}
}
}
}
}
void get_block_chroma(int uv, int ref_frame, StorablePicture **list, int x_pos, int y_pos, int hor_block_size, int ver_block_size, struct img_par *img, imgpel block[MB_BLOCK_SIZE][MB_BLOCK_SIZE])
{
int subpel_x = img->mb_cr_size_x == 8 ? 7 : 3;
int subpel_y = img->mb_cr_size_y == 8 ? 7 : 3;
int shiftpel_x = img->mb_cr_size_x == 8 ? 3 : 2;
int shiftpel_y = img->mb_cr_size_y == 8 ? 3 : 2;
int total_scale = shiftpel_x + shiftpel_y;
int dx = (x_pos & subpel_x);
int dy = (y_pos & subpel_y);
int dxcur = (subpel_x + 1 - dx);
int dycur = (subpel_y + 1 - dy);
int i, j;
int maxold_x = dec_picture->size_x_cr_m1;
int maxold_y = (dec_picture->mb_field[img->current_mb_nr]) ? (dec_picture->size_y_cr >> 1) - 1 : dec_picture->size_y_cr_m1;
int result;
StorablePicture *curr_ref = list[ref_frame];
static imgpel **cur_img;
static imgpel *cur_line, *cur_line_p1;
int tmp_pos;
static int ipos, ipos_p1;
if (curr_ref == no_reference_picture && img->framepoc < img->recovery_poc)
{
printf("list[ref_frame] is equal to 'no reference picture' before RAP\n");
/* fill the block with sample value 128 */
for (j = 0; j < ver_block_size; j++)
for (i = 0; i < hor_block_size; i++)
block[j][i] = 128;
return;
}
cur_img = curr_ref->imgUV[uv];
x_pos = x_pos >> shiftpel_x;
y_pos = y_pos >> shiftpel_y;
if ((y_pos >= 0) && (y_pos < maxold_y - ver_block_size) && (x_pos >= 0) && (x_pos < maxold_x - hor_block_size))
{
if (dx == 0 && dy == 0)
{ /* fullpel position */
for (j = 0; j < ver_block_size; j++)
{
memcpy(&(block[j][0]), &(cur_img[ y_pos + j ][x_pos]), hor_block_size * sizeof(imgpel));
}
}
else if (dx == 0)
{
for (j = 0; j < ver_block_size; j++)
{
cur_line = cur_img[y_pos + j];
cur_line_p1 = cur_img[y_pos + j + 1];
tmp_pos = x_pos;
for (i = 0; i < hor_block_size; i++)
{
ipos = tmp_pos++;
result = (
dxcur * dycur * cur_line [ipos ] +
dxcur * dy * cur_line_p1[ipos ]);
block[j][i] = iClip1(img->max_imgpel_value_uv, rshift_rnd_sf(result, total_scale));
}
}
}
else if (dy == 0)
{
for (j = 0; j < ver_block_size; j++)
{
cur_line = cur_img[y_pos + j];
tmp_pos = x_pos;
for (i = 0; i < hor_block_size; i++)
{
ipos = tmp_pos++;
ipos_p1 = tmp_pos;
result = (
dxcur * dycur * cur_line [ipos ] +
dx * dycur * cur_line [ipos_p1]);
block[j][i] = iClip1(img->max_imgpel_value_uv, rshift_rnd_sf(result, total_scale));
}
}
}
else
{ /* other positions */
for (j = 0; j < ver_block_size; j++)
{
cur_line = cur_img[y_pos + j];
cur_line_p1 = cur_img[y_pos + j + 1];
tmp_pos = x_pos;
for (i = 0; i < hor_block_size; i++)
{
ipos = tmp_pos++;
ipos_p1 = tmp_pos;
result = (
dxcur * dycur * cur_line [ipos ] +
dx * dycur * cur_line [ipos_p1] +
dxcur * dy * cur_line_p1[ipos ] +
dx * dy * cur_line_p1[ipos_p1]);
block[j][i] = iClip1(img->max_imgpel_value_uv, rshift_rnd_sf(result, total_scale));
}
}
}
}
else // unsafe positions
{
if (dx == 0 && dy == 0)
{ /* fullpel position */
for (j = 0; j < ver_block_size; j++)
{
cur_line = cur_img[iClip3(0, maxold_y, y_pos + j)];
for (i = 0; i < hor_block_size; i++)
{
block[j][i] = cur_line[iClip3(0, maxold_x, x_pos + i )];
}
}
}
else if (dx == 0)
{
for (j = 0; j < ver_block_size; j++)
{
cur_line = cur_img[iClip3(0, maxold_y, y_pos + j)];
cur_line_p1 = cur_img[iClip3(0, maxold_y, y_pos + j + 1)];
tmp_pos = x_pos;
for (i = 0; i < hor_block_size; i++)
{
ipos = iClip3(0, maxold_x, tmp_pos++);
result = (
dxcur * dycur * cur_line [ipos ] +
dxcur * dy * cur_line_p1[ipos ]);
block[j][i] = iClip1(img->max_imgpel_value_uv, rshift_rnd_sf(result, total_scale));
}
}
}
else if (dy == 0)
{
for (j = 0; j < ver_block_size; j++)
{
cur_line = cur_img[iClip3(0, maxold_y, y_pos + j)];
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -