diff --git a/libavcodec/imgconvert.c b/libavcodec/imgconvert.c index 469677b73b..403336812d 100644 --- a/libavcodec/imgconvert.c +++ b/libavcodec/imgconvert.c @@ -19,7 +19,9 @@ #include #include #include + #include "avcodec.h" +#include "dsputil.h" #ifdef USE_FASTMEMCPY #include "fastmemcpy.h" @@ -182,36 +184,366 @@ static void bgr24_to_yuv420p(UINT8 *lum, UINT8 *cb, UINT8 *cr, } } -int img_convert_to_yuv420(UINT8 *img_out, UINT8 *img, +/* XXX: use generic filter ? */ +/* 1x2 -> 1x1 */ +static void shrink2(UINT8 *dst, int dst_wrap, + UINT8 *src, int src_wrap, + int width, int height) +{ + int w; + UINT8 *s1, *s2, *d; + + for(;height > 0; height--) { + s1 = src; + s2 = s1 + src_wrap; + d = dst; + for(w = width;w >= 4; w-=4) { + d[0] = (s1[0] + s2[0]) >> 1; + d[1] = (s1[1] + s2[1]) >> 1; + d[2] = (s1[2] + s2[2]) >> 1; + d[3] = (s1[3] + s2[3]) >> 1; + s1 += 4; + s2 += 4; + d += 4; + } + for(;w > 0; w--) { + d[0] = (s1[0] + s2[0]) >> 1; + s1++; + s2++; + d++; + } + src += 2 * src_wrap; + dst += dst_wrap; + } +} + +/* 2x2 -> 1x1 */ +static void shrink22(UINT8 *dst, int dst_wrap, + UINT8 *src, int src_wrap, + int width, int height) +{ + int w; + UINT8 *s1, *s2, *d; + + for(;height > 0; height--) { + s1 = src; + s2 = s1 + src_wrap; + d = dst; + for(w = width;w >= 4; w-=4) { + d[0] = (s1[0] + s1[1] + s2[0] + s2[1] + 2) >> 1; + d[1] = (s1[2] + s1[3] + s2[2] + s2[3] + 2) >> 1; + d[2] = (s1[4] + s1[5] + s2[4] + s2[5] + 2) >> 1; + d[3] = (s1[6] + s1[7] + s2[6] + s2[7] + 2) >> 1; + s1 += 8; + s2 += 8; + d += 4; + } + for(;w > 0; w--) { + d[0] = (s1[0] + s1[1] + s2[0] + s2[1] + 2) >> 1; + s1 += 2; + s2 += 2; + d++; + } + src += 2 * src_wrap; + dst += dst_wrap; + } +} + +static void img_copy(UINT8 *dst, int dst_wrap, + UINT8 *src, int src_wrap, + int width, int height) +{ + for(;height > 0; height--) { + memcpy(dst, src, width); + dst += dst_wrap; + src += src_wrap; + } +} + +#define SCALE_BITS 10 + +#define C_Y (76309 >> (16 - SCALE_BITS)) +#define C_RV (117504 >> (16 - SCALE_BITS)) +#define C_BU (138453 >> (16 - SCALE_BITS)) +#define C_GU (13954 >> (16 - SCALE_BITS)) +#define C_GV (34903 >> (16 - SCALE_BITS)) + +#define RGBOUT(r, g, b, y1)\ +{\ + y = (y1 - 16) * C_Y;\ + r = cm[(y + r_add) >> SCALE_BITS];\ + g = cm[(y + g_add) >> SCALE_BITS];\ + b = cm[(y + b_add) >> SCALE_BITS];\ +} + +/* XXX: no chroma interpolating is done */ +static void yuv420p_to_rgb24(AVPicture *dst, AVPicture *src, + int width, int height) +{ + UINT8 *y1_ptr, *y2_ptr, *cb_ptr, *cr_ptr, *d, *d1, *d2; + int w, y, cb, cr, r_add, g_add, b_add, width2; + UINT8 *cm = cropTbl + MAX_NEG_CROP; + + d = dst->data[0]; + y1_ptr = src->data[0]; + cb_ptr = src->data[1]; + cr_ptr = src->data[2]; + width2 = width >> 1; + for(;height > 0; height -= 2) { + d1 = d; + d2 = d + dst->linesize[0]; + y2_ptr = y1_ptr + src->linesize[0]; + for(w = width2; w > 0; w --) { + cb = cb_ptr[0] - 128; + cr = cr_ptr[0] - 128; + r_add = C_RV * cr + (1 << (SCALE_BITS - 1)); + g_add = - C_GU * cb - C_GV * cr + (1 << (SCALE_BITS - 1)); + b_add = C_BU * cb + (1 << (SCALE_BITS - 1)); + + /* output 4 pixels */ + RGBOUT(d1[0], d1[1], d1[2], y1_ptr[0]); + RGBOUT(d1[3], d1[4], d1[5], y1_ptr[1]); + RGBOUT(d2[0], d2[1], d2[2], y2_ptr[0]); + RGBOUT(d2[3], d2[4], d2[5], y2_ptr[1]); + + d1 += 6; + d2 += 6; + y1_ptr += 2; + y2_ptr += 2; + cb_ptr++; + cr_ptr++; + } + d += 2 * dst->linesize[0]; + y1_ptr += 2 * src->linesize[0] - width; + cb_ptr += src->linesize[1] - width2; + cr_ptr += src->linesize[2] - width2; + } +} + +/* XXX: no chroma interpolating is done */ +static void yuv422p_to_rgb24(AVPicture *dst, AVPicture *src, + int width, int height) +{ + UINT8 *y1_ptr, *cb_ptr, *cr_ptr, *d, *d1; + int w, y, cb, cr, r_add, g_add, b_add, width2; + UINT8 *cm = cropTbl + MAX_NEG_CROP; + + d = dst->data[0]; + y1_ptr = src->data[0]; + cb_ptr = src->data[1]; + cr_ptr = src->data[2]; + width2 = width >> 1; + for(;height > 0; height --) { + d1 = d; + for(w = width2; w > 0; w --) { + cb = cb_ptr[0] - 128; + cr = cr_ptr[0] - 128; + r_add = C_RV * cr + (1 << (SCALE_BITS - 1)); + g_add = - C_GU * cb - C_GV * cr + (1 << (SCALE_BITS - 1)); + b_add = C_BU * cb + (1 << (SCALE_BITS - 1)); + + /* output 2 pixels */ + RGBOUT(d1[0], d1[1], d1[2], y1_ptr[0]); + RGBOUT(d1[3], d1[4], d1[5], y1_ptr[1]); + + d1 += 6; + y1_ptr += 2; + cb_ptr++; + cr_ptr++; + } + d += dst->linesize[0]; + y1_ptr += src->linesize[0] - width; + cb_ptr += src->linesize[1] - width2; + cr_ptr += src->linesize[2] - width2; + } +} + +/* XXX: always use linesize. Return -1 if not supported */ +int img_convert(AVPicture *dst, int dst_pix_fmt, + AVPicture *src, int pix_fmt, + int width, int height) +{ + int i; + + if (dst_pix_fmt == pix_fmt) { + switch(pix_fmt) { + case PIX_FMT_YUV420P: + for(i=0;i<3;i++) { + if (i == 1) { + width >>= 1; + height >>= 1; + } + img_copy(dst->data[i], dst->linesize[i], + src->data[i], src->linesize[i], + width, height); + } + break; + default: + return -1; + } + } else if (dst_pix_fmt == PIX_FMT_YUV420P) { + + switch(pix_fmt) { + case PIX_FMT_YUV420P: + for(i=0;i<3;i++) { + img_copy(dst->data[i], dst->linesize[i], + src->data[i], src->linesize[i], + width, height); + } + break; + case PIX_FMT_YUV422P: + img_copy(dst->data[0], dst->linesize[0], + src->data[0], src->linesize[0], + width, height); + width >>= 1; + height >>= 1; + for(i=1;i<3;i++) { + shrink2(dst->data[i], dst->linesize[i], + src->data[i], src->linesize[i], + width, height); + } + break; + case PIX_FMT_YUV444P: + img_copy(dst->data[0], dst->linesize[0], + src->data[0], src->linesize[0], + width, height); + width >>= 1; + height >>= 1; + for(i=1;i<3;i++) { + shrink22(dst->data[i], dst->linesize[i], + src->data[i], src->linesize[i], + width, height); + } + break; + case PIX_FMT_YUV422: + yuv422_to_yuv420p(dst->data[0], dst->data[1], dst->data[2], + src->data[0], width, height); + break; + case PIX_FMT_RGB24: + rgb24_to_yuv420p(dst->data[0], dst->data[1], dst->data[2], + src->data[0], width, height); + break; + case PIX_FMT_BGR24: + bgr24_to_yuv420p(dst->data[0], dst->data[1], dst->data[2], + src->data[0], width, height); + break; + default: + return -1; + } + } else if (dst_pix_fmt == PIX_FMT_RGB24) { + switch(pix_fmt) { + case PIX_FMT_YUV420P: + yuv420p_to_rgb24(dst, src, width, height); + break; + case PIX_FMT_YUV422P: + yuv422p_to_rgb24(dst, src, width, height); + break; + default: + return -1; + } + } else { + return -1; + } + return 0; +} + +/* filter parameters: [-1 4 2 4 -1] // 8 */ +static void deinterlace_line(UINT8 *dst, UINT8 *src, int src_wrap, + int size) +{ + UINT8 *cm = cropTbl + MAX_NEG_CROP; + int sum; + UINT8 *s; + + for(;size > 0;size--) { + s = src; + sum = -s[0]; + s += src_wrap; + sum += s[0] << 2; + s += src_wrap; + sum += s[0] << 1; + s += src_wrap; + sum += s[0] << 2; + s += src_wrap; + sum += -s[0]; + dst[0] = cm[(sum + 4) >> 3]; + dst++; + src++; + } +} + +/* deinterlacing : 2 temporal taps, 3 spatial taps linear filter. The + top field is copied as is, but the bottom field is deinterlaced + against the top field. */ +static void deinterlace_bottom_field(UINT8 *dst, int dst_wrap, + UINT8 *src1, int src_wrap, + int width, int height) +{ + UINT8 *src, *ptr; + int y, y1, i; + UINT8 buf[5 * width]; + + src = src1; + for(y=0;y= 0 && (y1 + 4) < height) { + /* fast case : no edges */ + deinterlace_line(dst, src, src_wrap, width); + } else { + /* in order to use the same function, we use an intermediate buffer */ + ptr = buf; + for(i=0;i<5;i++) { + if (y1 < 0) + memcpy(ptr, src1, width); + else if (y1 >= height) + memcpy(ptr, src1 + (height - 1) * src_wrap, width); + else + memcpy(ptr, src1 + y1 * src_wrap, width); + y1++; + ptr += width; + } + deinterlace_line(dst, buf, width, width); + } + dst += dst_wrap; + src += (2 + 1) * src_wrap; + } +} + + +/* deinterlace, return -1 if format not handled */ +int avpicture_deinterlace(AVPicture *dst, AVPicture *src, int pix_fmt, int width, int height) { - UINT8 *pict; - int size, size_out; - UINT8 *picture[3]; + int i; - pict = img_out; - size = width * height; - size_out = (size * 3) / 2; - picture[0] = pict; - picture[1] = pict + size; - picture[2] = picture[1] + (size / 4); - - switch(pix_fmt) { - case PIX_FMT_YUV420P: - memcpy(pict, img, size_out); - break; - case PIX_FMT_YUV422: - yuv422_to_yuv420p(picture[0], picture[1], picture[2], - img, width, height); - break; - case PIX_FMT_RGB24: - rgb24_to_yuv420p(picture[0], picture[1], picture[2], - img, width, height); - break; - case PIX_FMT_BGR24: - bgr24_to_yuv420p(picture[0], picture[1], picture[2], - img, width, height); - break; + if (pix_fmt != PIX_FMT_YUV420P && + pix_fmt != PIX_FMT_YUV422P && + pix_fmt != PIX_FMT_YUV444P) + return -1; + if ((width & 1) != 0 || (height & 3) != 0) + return -1; + + for(i=0;i<3;i++) { + if (i == 1) { + switch(pix_fmt) { + case PIX_FMT_YUV420P: + width >>= 1; + height >>= 1; + break; + case PIX_FMT_YUV422P: + width >>= 1; + break; + default: + break; + } + } + deinterlace_bottom_field(dst->data[i], dst->linesize[i], + src->data[i], src->linesize[i], + width, height); } - return size_out; + return 0; }