yading@10: /*
yading@10:  * Misc image conversion routines
yading@10:  * Copyright (c) 2001, 2002, 2003 Fabrice Bellard
yading@10:  *
yading@10:  * This file is part of FFmpeg.
yading@10:  *
yading@10:  * FFmpeg is free software; you can redistribute it and/or
yading@10:  * modify it under the terms of the GNU Lesser General Public
yading@10:  * License as published by the Free Software Foundation; either
yading@10:  * version 2.1 of the License, or (at your option) any later version.
yading@10:  *
yading@10:  * FFmpeg is distributed in the hope that it will be useful,
yading@10:  * but WITHOUT ANY WARRANTY; without even the implied warranty of
yading@10:  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
yading@10:  * Lesser General Public License for more details.
yading@10:  *
yading@10:  * You should have received a copy of the GNU Lesser General Public
yading@10:  * License along with FFmpeg; if not, write to the Free Software
yading@10:  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
yading@10:  */
yading@10: 
yading@10: /**
yading@10:  * @file
yading@10:  * misc image conversion routines
yading@10:  */
yading@10: 
yading@10: /* TODO:
yading@10:  * - write 'ffimg' program to test all the image related stuff
yading@10:  * - move all api to slice based system
yading@10:  * - integrate deinterlacing, postprocessing and scaling in the conversion process
yading@10:  */
yading@10: 
yading@10: #include "avcodec.h"
yading@10: #include "dsputil.h"
yading@10: #include "imgconvert.h"
yading@10: #include "internal.h"
yading@10: #include "libavutil/avassert.h"
yading@10: #include "libavutil/colorspace.h"
yading@10: #include "libavutil/common.h"
yading@10: #include "libavutil/pixdesc.h"
yading@10: #include "libavutil/imgutils.h"
yading@10: 
yading@10: #if HAVE_MMX_EXTERNAL
yading@10: #include "x86/dsputil_mmx.h"
yading@10: #endif
yading@10: 
yading@10: #define FF_COLOR_NA      -1
yading@10: #define FF_COLOR_RGB      0 /**< RGB color space */
yading@10: #define FF_COLOR_GRAY     1 /**< gray color space */
yading@10: #define FF_COLOR_YUV      2 /**< YUV color space. 16 <= Y <= 235, 16 <= U, V <= 240 */
yading@10: #define FF_COLOR_YUV_JPEG 3 /**< YUV color space. 0 <= Y <= 255, 0 <= U, V <= 255 */
yading@10: 
yading@10: #if HAVE_MMX_EXTERNAL
yading@10: #define deinterlace_line_inplace ff_deinterlace_line_inplace_mmx
yading@10: #define deinterlace_line         ff_deinterlace_line_mmx
yading@10: #else
yading@10: #define deinterlace_line_inplace deinterlace_line_inplace_c
yading@10: #define deinterlace_line         deinterlace_line_c
yading@10: #endif
yading@10: 
yading@10: #define pixdesc_has_alpha(pixdesc) \
yading@10:     ((pixdesc)->nb_components == 2 || (pixdesc)->nb_components == 4 || (pixdesc)->flags & PIX_FMT_PAL)
yading@10: 
yading@10: 
yading@10: void avcodec_get_chroma_sub_sample(enum AVPixelFormat pix_fmt, int *h_shift, int *v_shift)
yading@10: {
yading@10:     const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(pix_fmt);
yading@10:     av_assert0(desc);
yading@10:     *h_shift = desc->log2_chroma_w;
yading@10:     *v_shift = desc->log2_chroma_h;
yading@10: }
yading@10: 
yading@10: static int get_color_type(const AVPixFmtDescriptor *desc) {
yading@10:     if(desc->nb_components == 1 || desc->nb_components == 2)
yading@10:         return FF_COLOR_GRAY;
yading@10: 
yading@10:     if(desc->name && !strncmp(desc->name, "yuvj", 4))
yading@10:         return FF_COLOR_YUV_JPEG;
yading@10: 
yading@10:     if(desc->flags & PIX_FMT_RGB)
yading@10:         return  FF_COLOR_RGB;
yading@10: 
yading@10:     if(desc->nb_components == 0)
yading@10:         return FF_COLOR_NA;
yading@10: 
yading@10:     return FF_COLOR_YUV;
yading@10: }
yading@10: 
yading@10: static int get_pix_fmt_depth(int *min, int *max, enum AVPixelFormat pix_fmt)
yading@10: {
yading@10:     const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(pix_fmt);
yading@10:     int i;
yading@10: 
yading@10:     if (!desc || !desc->nb_components) {
yading@10:         *min = *max = 0;
yading@10:         return AVERROR(EINVAL);
yading@10:     }
yading@10: 
yading@10:     *min = INT_MAX, *max = -INT_MAX;
yading@10:     for (i = 0; i < desc->nb_components; i++) {
yading@10:         *min = FFMIN(desc->comp[i].depth_minus1+1, *min);
yading@10:         *max = FFMAX(desc->comp[i].depth_minus1+1, *max);
yading@10:     }
yading@10:     return 0;
yading@10: }
yading@10: 
yading@10: static int get_pix_fmt_score(enum AVPixelFormat dst_pix_fmt,
yading@10:                               enum AVPixelFormat src_pix_fmt,
yading@10:                               unsigned *lossp, unsigned consider)
yading@10: {
yading@10:     const AVPixFmtDescriptor *src_desc = av_pix_fmt_desc_get(src_pix_fmt);
yading@10:     const AVPixFmtDescriptor *dst_desc = av_pix_fmt_desc_get(dst_pix_fmt);
yading@10:     int src_color, dst_color;
yading@10:     int src_min_depth, src_max_depth, dst_min_depth, dst_max_depth;
yading@10:     int ret, loss, i, nb_components;
yading@10:     int score = INT_MAX;
yading@10: 
yading@10:     if (dst_pix_fmt >= AV_PIX_FMT_NB || dst_pix_fmt <= AV_PIX_FMT_NONE)
yading@10:         return ~0;
yading@10: 
yading@10:     /* compute loss */
yading@10:     *lossp = loss = 0;
yading@10: 
yading@10:     if (dst_pix_fmt == src_pix_fmt)
yading@10:         return INT_MAX;
yading@10: 
yading@10:     if ((ret = get_pix_fmt_depth(&src_min_depth, &src_max_depth, src_pix_fmt)) < 0)
yading@10:         return ret;
yading@10:     if ((ret = get_pix_fmt_depth(&dst_min_depth, &dst_max_depth, dst_pix_fmt)) < 0)
yading@10:         return ret;
yading@10: 
yading@10:     src_color = get_color_type(src_desc);
yading@10:     dst_color = get_color_type(dst_desc);
yading@10:     nb_components = FFMIN(src_desc->nb_components, dst_desc->nb_components);
yading@10: 
yading@10:     for (i = 0; i < nb_components; i++)
yading@10:         if (src_desc->comp[i].depth_minus1 > dst_desc->comp[i].depth_minus1 && (consider & FF_LOSS_DEPTH)) {
yading@10:             loss |= FF_LOSS_DEPTH;
yading@10:             score -= 65536 >> dst_desc->comp[i].depth_minus1;
yading@10:         }
yading@10: 
yading@10:     if (consider & FF_LOSS_RESOLUTION) {
yading@10:         if (dst_desc->log2_chroma_w > src_desc->log2_chroma_w) {
yading@10:             loss |= FF_LOSS_RESOLUTION;
yading@10:             score -= 256 << dst_desc->log2_chroma_w;
yading@10:         }
yading@10:         if (dst_desc->log2_chroma_h > src_desc->log2_chroma_h) {
yading@10:             loss |= FF_LOSS_RESOLUTION;
yading@10:             score -= 256 << dst_desc->log2_chroma_h;
yading@10:         }
yading@10:         // dont favor 422 over 420 if downsampling is needed, because 420 has much better support on the decoder side
yading@10:         if (dst_desc->log2_chroma_w == 1 && src_desc->log2_chroma_w == 0 &&
yading@10:             dst_desc->log2_chroma_h == 1 && src_desc->log2_chroma_h == 0 ) {
yading@10:             score += 512;
yading@10:         }
yading@10:     }
yading@10: 
yading@10:     if(consider & FF_LOSS_COLORSPACE)
yading@10:     switch(dst_color) {
yading@10:     case FF_COLOR_RGB:
yading@10:         if (src_color != FF_COLOR_RGB &&
yading@10:             src_color != FF_COLOR_GRAY)
yading@10:             loss |= FF_LOSS_COLORSPACE;
yading@10:         break;
yading@10:     case FF_COLOR_GRAY:
yading@10:         if (src_color != FF_COLOR_GRAY)
yading@10:             loss |= FF_LOSS_COLORSPACE;
yading@10:         break;
yading@10:     case FF_COLOR_YUV:
yading@10:         if (src_color != FF_COLOR_YUV)
yading@10:             loss |= FF_LOSS_COLORSPACE;
yading@10:         break;
yading@10:     case FF_COLOR_YUV_JPEG:
yading@10:         if (src_color != FF_COLOR_YUV_JPEG &&
yading@10:             src_color != FF_COLOR_YUV &&
yading@10:             src_color != FF_COLOR_GRAY)
yading@10:             loss |= FF_LOSS_COLORSPACE;
yading@10:         break;
yading@10:     default:
yading@10:         /* fail safe test */
yading@10:         if (src_color != dst_color)
yading@10:             loss |= FF_LOSS_COLORSPACE;
yading@10:         break;
yading@10:     }
yading@10:     if(loss & FF_LOSS_COLORSPACE)
yading@10:         score -= (nb_components * 65536) >> FFMIN(dst_desc->comp[0].depth_minus1, src_desc->comp[0].depth_minus1);
yading@10: 
yading@10:     if (dst_color == FF_COLOR_GRAY &&
yading@10:         src_color != FF_COLOR_GRAY && (consider & FF_LOSS_CHROMA)) {
yading@10:         loss |= FF_LOSS_CHROMA;
yading@10:         score -= 2 * 65536;
yading@10:     }
yading@10:     if (!pixdesc_has_alpha(dst_desc) && (pixdesc_has_alpha(src_desc) && (consider & FF_LOSS_ALPHA))) {
yading@10:         loss |= FF_LOSS_ALPHA;
yading@10:         score -= 65536;
yading@10:     }
yading@10:     if (dst_pix_fmt == AV_PIX_FMT_PAL8 && (consider & FF_LOSS_COLORQUANT) &&
yading@10:         (src_pix_fmt != AV_PIX_FMT_PAL8 && (src_color != FF_COLOR_GRAY || (pixdesc_has_alpha(src_desc) && (consider & FF_LOSS_ALPHA))))) {
yading@10:         loss |= FF_LOSS_COLORQUANT;
yading@10:         score -= 65536;
yading@10:     }
yading@10: 
yading@10:     *lossp = loss;
yading@10:     return score;
yading@10: }
yading@10: 
yading@10: int avcodec_get_pix_fmt_loss(enum AVPixelFormat dst_pix_fmt,
yading@10:                              enum AVPixelFormat src_pix_fmt,
yading@10:                              int has_alpha)
yading@10: {
yading@10:     int loss;
yading@10:     int ret = get_pix_fmt_score(dst_pix_fmt, src_pix_fmt, &loss, has_alpha ? ~0 : ~FF_LOSS_ALPHA);
yading@10:     if (ret < 0)
yading@10:         return ret;
yading@10:     return loss;
yading@10: }
yading@10: 
yading@10: enum AVPixelFormat avcodec_find_best_pix_fmt_of_2(enum AVPixelFormat dst_pix_fmt1, enum AVPixelFormat dst_pix_fmt2,
yading@10:                                             enum AVPixelFormat src_pix_fmt, int has_alpha, int *loss_ptr)
yading@10: {
yading@10:     enum AVPixelFormat dst_pix_fmt;
yading@10:     int loss1, loss2, loss_mask;
yading@10:     const AVPixFmtDescriptor *desc1 = av_pix_fmt_desc_get(dst_pix_fmt1);
yading@10:     const AVPixFmtDescriptor *desc2 = av_pix_fmt_desc_get(dst_pix_fmt2);
yading@10:     int score1, score2;
yading@10: 
yading@10:     loss_mask= loss_ptr?~*loss_ptr:~0; /* use loss mask if provided */
yading@10:     if(!has_alpha)
yading@10:         loss_mask &= ~FF_LOSS_ALPHA;
yading@10: 
yading@10:     dst_pix_fmt = AV_PIX_FMT_NONE;
yading@10:     score1 = get_pix_fmt_score(dst_pix_fmt1, src_pix_fmt, &loss1, loss_mask);
yading@10:     score2 = get_pix_fmt_score(dst_pix_fmt2, src_pix_fmt, &loss2, loss_mask);
yading@10: 
yading@10:     if (score1 == score2) {
yading@10:         if(av_get_padded_bits_per_pixel(desc2) != av_get_padded_bits_per_pixel(desc1)) {
yading@10:             dst_pix_fmt = av_get_padded_bits_per_pixel(desc2) < av_get_padded_bits_per_pixel(desc1) ? dst_pix_fmt2 : dst_pix_fmt1;
yading@10:         } else {
yading@10:             dst_pix_fmt = desc2->nb_components < desc1->nb_components ? dst_pix_fmt2 : dst_pix_fmt1;
yading@10:         }
yading@10:     } else {
yading@10:         dst_pix_fmt = score1 < score2 ? dst_pix_fmt2 : dst_pix_fmt1;
yading@10:     }
yading@10: 
yading@10:     if (loss_ptr)
yading@10:         *loss_ptr = avcodec_get_pix_fmt_loss(dst_pix_fmt, src_pix_fmt, has_alpha);
yading@10:     return dst_pix_fmt;
yading@10: }
yading@10: 
yading@10: #if AV_HAVE_INCOMPATIBLE_FORK_ABI
yading@10: enum AVPixelFormat avcodec_find_best_pix_fmt2(enum AVPixelFormat *pix_fmt_list,
yading@10:                                             enum AVPixelFormat src_pix_fmt,
yading@10:                                             int has_alpha, int *loss_ptr){
yading@10:     return avcodec_find_best_pix_fmt_of_list(pix_fmt_list, src_pix_fmt, has_alpha, loss_ptr);
yading@10: }
yading@10: #else
yading@10: enum AVPixelFormat avcodec_find_best_pix_fmt2(enum AVPixelFormat dst_pix_fmt1, enum AVPixelFormat dst_pix_fmt2,
yading@10:                                             enum AVPixelFormat src_pix_fmt, int has_alpha, int *loss_ptr)
yading@10: {
yading@10:     return avcodec_find_best_pix_fmt_of_2(dst_pix_fmt1, dst_pix_fmt2, src_pix_fmt, has_alpha, loss_ptr);
yading@10: }
yading@10: #endif
yading@10: 
yading@10: enum AVPixelFormat avcodec_find_best_pix_fmt_of_list(enum AVPixelFormat *pix_fmt_list,
yading@10:                                             enum AVPixelFormat src_pix_fmt,
yading@10:                                             int has_alpha, int *loss_ptr){
yading@10:     int i;
yading@10: 
yading@10:     enum AVPixelFormat best = AV_PIX_FMT_NONE;
yading@10: 
yading@10:     for(i=0; pix_fmt_list[i] != AV_PIX_FMT_NONE; i++)
yading@10:         best = avcodec_find_best_pix_fmt_of_2(best, pix_fmt_list[i], src_pix_fmt, has_alpha, loss_ptr);
yading@10: 
yading@10:     return best;
yading@10: }
yading@10: 
yading@10: /* 2x2 -> 1x1 */
yading@10: void ff_shrink22(uint8_t *dst, int dst_wrap,
yading@10:                      const uint8_t *src, int src_wrap,
yading@10:                      int width, int height)
yading@10: {
yading@10:     int w;
yading@10:     const uint8_t *s1, *s2;
yading@10:     uint8_t *d;
yading@10: 
yading@10:     for(;height > 0; height--) {
yading@10:         s1 = src;
yading@10:         s2 = s1 + src_wrap;
yading@10:         d = dst;
yading@10:         for(w = width;w >= 4; w-=4) {
yading@10:             d[0] = (s1[0] + s1[1] + s2[0] + s2[1] + 2) >> 2;
yading@10:             d[1] = (s1[2] + s1[3] + s2[2] + s2[3] + 2) >> 2;
yading@10:             d[2] = (s1[4] + s1[5] + s2[4] + s2[5] + 2) >> 2;
yading@10:             d[3] = (s1[6] + s1[7] + s2[6] + s2[7] + 2) >> 2;
yading@10:             s1 += 8;
yading@10:             s2 += 8;
yading@10:             d += 4;
yading@10:         }
yading@10:         for(;w > 0; w--) {
yading@10:             d[0] = (s1[0] + s1[1] + s2[0] + s2[1] + 2) >> 2;
yading@10:             s1 += 2;
yading@10:             s2 += 2;
yading@10:             d++;
yading@10:         }
yading@10:         src += 2 * src_wrap;
yading@10:         dst += dst_wrap;
yading@10:     }
yading@10: }
yading@10: 
yading@10: /* 4x4 -> 1x1 */
yading@10: void ff_shrink44(uint8_t *dst, int dst_wrap,
yading@10:                      const uint8_t *src, int src_wrap,
yading@10:                      int width, int height)
yading@10: {
yading@10:     int w;
yading@10:     const uint8_t *s1, *s2, *s3, *s4;
yading@10:     uint8_t *d;
yading@10: 
yading@10:     for(;height > 0; height--) {
yading@10:         s1 = src;
yading@10:         s2 = s1 + src_wrap;
yading@10:         s3 = s2 + src_wrap;
yading@10:         s4 = s3 + src_wrap;
yading@10:         d = dst;
yading@10:         for(w = width;w > 0; w--) {
yading@10:             d[0] = (s1[0] + s1[1] + s1[2] + s1[3] +
yading@10:                     s2[0] + s2[1] + s2[2] + s2[3] +
yading@10:                     s3[0] + s3[1] + s3[2] + s3[3] +
yading@10:                     s4[0] + s4[1] + s4[2] + s4[3] + 8) >> 4;
yading@10:             s1 += 4;
yading@10:             s2 += 4;
yading@10:             s3 += 4;
yading@10:             s4 += 4;
yading@10:             d++;
yading@10:         }
yading@10:         src += 4 * src_wrap;
yading@10:         dst += dst_wrap;
yading@10:     }
yading@10: }
yading@10: 
yading@10: /* 8x8 -> 1x1 */
yading@10: void ff_shrink88(uint8_t *dst, int dst_wrap,
yading@10:                      const uint8_t *src, int src_wrap,
yading@10:                      int width, int height)
yading@10: {
yading@10:     int w, i;
yading@10: 
yading@10:     for(;height > 0; height--) {
yading@10:         for(w = width;w > 0; w--) {
yading@10:             int tmp=0;
yading@10:             for(i=0; i<8; i++){
yading@10:                 tmp += src[0] + src[1] + src[2] + src[3] + src[4] + src[5] + src[6] + src[7];
yading@10:                 src += src_wrap;
yading@10:             }
yading@10:             *(dst++) = (tmp + 32)>>6;
yading@10:             src += 8 - 8*src_wrap;
yading@10:         }
yading@10:         src += 8*src_wrap - 8*width;
yading@10:         dst += dst_wrap - width;
yading@10:     }
yading@10: }
yading@10: 
yading@10: /* return true if yuv planar */
yading@10: static inline int is_yuv_planar(const AVPixFmtDescriptor *desc)
yading@10: {
yading@10:     int i;
yading@10:     int planes[4] = { 0 };
yading@10: 
yading@10:     if (     desc->flags & PIX_FMT_RGB
yading@10:         || !(desc->flags & PIX_FMT_PLANAR))
yading@10:         return 0;
yading@10: 
yading@10:     /* set the used planes */
yading@10:     for (i = 0; i < desc->nb_components; i++)
yading@10:         planes[desc->comp[i].plane] = 1;
yading@10: 
yading@10:     /* if there is an unused plane, the format is not planar */
yading@10:     for (i = 0; i < desc->nb_components; i++)
yading@10:         if (!planes[i])
yading@10:             return 0;
yading@10:     return 1;
yading@10: }
yading@10: 
yading@10: int av_picture_crop(AVPicture *dst, const AVPicture *src,
yading@10:                     enum AVPixelFormat pix_fmt, int top_band, int left_band)
yading@10: {
yading@10:     const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(pix_fmt);
yading@10:     int y_shift;
yading@10:     int x_shift;
yading@10: 
yading@10:     if (pix_fmt < 0 || pix_fmt >= AV_PIX_FMT_NB)
yading@10:         return -1;
yading@10: 
yading@10:     y_shift = desc->log2_chroma_h;
yading@10:     x_shift = desc->log2_chroma_w;
yading@10: 
yading@10:     if (is_yuv_planar(desc)) {
yading@10:     dst->data[0] = src->data[0] + (top_band * src->linesize[0]) + left_band;
yading@10:     dst->data[1] = src->data[1] + ((top_band >> y_shift) * src->linesize[1]) + (left_band >> x_shift);
yading@10:     dst->data[2] = src->data[2] + ((top_band >> y_shift) * src->linesize[2]) + (left_band >> x_shift);
yading@10:     } else{
yading@10:         if(top_band % (1<<y_shift) || left_band % (1<<x_shift))
yading@10:             return -1;
yading@10:         if(left_band) //FIXME add support for this too
yading@10:             return -1;
yading@10:         dst->data[0] = src->data[0] + (top_band * src->linesize[0]) + left_band;
yading@10:     }
yading@10: 
yading@10:     dst->linesize[0] = src->linesize[0];
yading@10:     dst->linesize[1] = src->linesize[1];
yading@10:     dst->linesize[2] = src->linesize[2];
yading@10:     return 0;
yading@10: }
yading@10: 
yading@10: int av_picture_pad(AVPicture *dst, const AVPicture *src, int height, int width,
yading@10:                    enum AVPixelFormat pix_fmt, int padtop, int padbottom, int padleft, int padright,
yading@10:             int *color)
yading@10: {
yading@10:     const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(pix_fmt);
yading@10:     uint8_t *optr;
yading@10:     int y_shift;
yading@10:     int x_shift;
yading@10:     int yheight;
yading@10:     int i, y;
yading@10: 
yading@10:     if (pix_fmt < 0 || pix_fmt >= AV_PIX_FMT_NB ||
yading@10:         !is_yuv_planar(desc)) return -1;
yading@10: 
yading@10:     for (i = 0; i < 3; i++) {
yading@10:         x_shift = i ? desc->log2_chroma_w : 0;
yading@10:         y_shift = i ? desc->log2_chroma_h : 0;
yading@10: 
yading@10:         if (padtop || padleft) {
yading@10:             memset(dst->data[i], color[i],
yading@10:                 dst->linesize[i] * (padtop >> y_shift) + (padleft >> x_shift));
yading@10:         }
yading@10: 
yading@10:         if (padleft || padright) {
yading@10:             optr = dst->data[i] + dst->linesize[i] * (padtop >> y_shift) +
yading@10:                 (dst->linesize[i] - (padright >> x_shift));
yading@10:             yheight = (height - 1 - (padtop + padbottom)) >> y_shift;
yading@10:             for (y = 0; y < yheight; y++) {
yading@10:                 memset(optr, color[i], (padleft + padright) >> x_shift);
yading@10:                 optr += dst->linesize[i];
yading@10:             }
yading@10:         }
yading@10: 
yading@10:         if (src) { /* first line */
yading@10:             uint8_t *iptr = src->data[i];
yading@10:             optr = dst->data[i] + dst->linesize[i] * (padtop >> y_shift) +
yading@10:                     (padleft >> x_shift);
yading@10:             memcpy(optr, iptr, (width - padleft - padright) >> x_shift);
yading@10:             iptr += src->linesize[i];
yading@10:             optr = dst->data[i] + dst->linesize[i] * (padtop >> y_shift) +
yading@10:                 (dst->linesize[i] - (padright >> x_shift));
yading@10:             yheight = (height - 1 - (padtop + padbottom)) >> y_shift;
yading@10:             for (y = 0; y < yheight; y++) {
yading@10:                 memset(optr, color[i], (padleft + padright) >> x_shift);
yading@10:                 memcpy(optr + ((padleft + padright) >> x_shift), iptr,
yading@10:                        (width - padleft - padright) >> x_shift);
yading@10:                 iptr += src->linesize[i];
yading@10:                 optr += dst->linesize[i];
yading@10:             }
yading@10:         }
yading@10: 
yading@10:         if (padbottom || padright) {
yading@10:             optr = dst->data[i] + dst->linesize[i] *
yading@10:                 ((height - padbottom) >> y_shift) - (padright >> x_shift);
yading@10:             memset(optr, color[i],dst->linesize[i] *
yading@10:                 (padbottom >> y_shift) + (padright >> x_shift));
yading@10:         }
yading@10:     }
yading@10:     return 0;
yading@10: }
yading@10: 
yading@10: #if FF_API_DEINTERLACE
yading@10: 
yading@10: #if !HAVE_MMX_EXTERNAL
yading@10: /* filter parameters: [-1 4 2 4 -1] // 8 */
yading@10: static void deinterlace_line_c(uint8_t *dst,
yading@10:                              const uint8_t *lum_m4, const uint8_t *lum_m3,
yading@10:                              const uint8_t *lum_m2, const uint8_t *lum_m1,
yading@10:                              const uint8_t *lum,
yading@10:                              int size)
yading@10: {
yading@10:     const uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;
yading@10:     int sum;
yading@10: 
yading@10:     for(;size > 0;size--) {
yading@10:         sum = -lum_m4[0];
yading@10:         sum += lum_m3[0] << 2;
yading@10:         sum += lum_m2[0] << 1;
yading@10:         sum += lum_m1[0] << 2;
yading@10:         sum += -lum[0];
yading@10:         dst[0] = cm[(sum + 4) >> 3];
yading@10:         lum_m4++;
yading@10:         lum_m3++;
yading@10:         lum_m2++;
yading@10:         lum_m1++;
yading@10:         lum++;
yading@10:         dst++;
yading@10:     }
yading@10: }
yading@10: 
yading@10: static void deinterlace_line_inplace_c(uint8_t *lum_m4, uint8_t *lum_m3,
yading@10:                                        uint8_t *lum_m2, uint8_t *lum_m1,
yading@10:                                        uint8_t *lum, int size)
yading@10: {
yading@10:     const uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;
yading@10:     int sum;
yading@10: 
yading@10:     for(;size > 0;size--) {
yading@10:         sum = -lum_m4[0];
yading@10:         sum += lum_m3[0] << 2;
yading@10:         sum += lum_m2[0] << 1;
yading@10:         lum_m4[0]=lum_m2[0];
yading@10:         sum += lum_m1[0] << 2;
yading@10:         sum += -lum[0];
yading@10:         lum_m2[0] = cm[(sum + 4) >> 3];
yading@10:         lum_m4++;
yading@10:         lum_m3++;
yading@10:         lum_m2++;
yading@10:         lum_m1++;
yading@10:         lum++;
yading@10:     }
yading@10: }
yading@10: #endif /* !HAVE_MMX_EXTERNAL */
yading@10: 
yading@10: /* deinterlacing : 2 temporal taps, 3 spatial taps linear filter. The
yading@10:    top field is copied as is, but the bottom field is deinterlaced
yading@10:    against the top field. */
yading@10: static void deinterlace_bottom_field(uint8_t *dst, int dst_wrap,
yading@10:                                     const uint8_t *src1, int src_wrap,
yading@10:                                     int width, int height)
yading@10: {
yading@10:     const uint8_t *src_m2, *src_m1, *src_0, *src_p1, *src_p2;
yading@10:     int y;
yading@10: 
yading@10:     src_m2 = src1;
yading@10:     src_m1 = src1;
yading@10:     src_0=&src_m1[src_wrap];
yading@10:     src_p1=&src_0[src_wrap];
yading@10:     src_p2=&src_p1[src_wrap];
yading@10:     for(y=0;y<(height-2);y+=2) {
yading@10:         memcpy(dst,src_m1,width);
yading@10:         dst += dst_wrap;
yading@10:         deinterlace_line(dst,src_m2,src_m1,src_0,src_p1,src_p2,width);
yading@10:         src_m2 = src_0;
yading@10:         src_m1 = src_p1;
yading@10:         src_0 = src_p2;
yading@10:         src_p1 += 2*src_wrap;
yading@10:         src_p2 += 2*src_wrap;
yading@10:         dst += dst_wrap;
yading@10:     }
yading@10:     memcpy(dst,src_m1,width);
yading@10:     dst += dst_wrap;
yading@10:     /* do last line */
yading@10:     deinterlace_line(dst,src_m2,src_m1,src_0,src_0,src_0,width);
yading@10: }
yading@10: 
yading@10: static void deinterlace_bottom_field_inplace(uint8_t *src1, int src_wrap,
yading@10:                                              int width, int height)
yading@10: {
yading@10:     uint8_t *src_m1, *src_0, *src_p1, *src_p2;
yading@10:     int y;
yading@10:     uint8_t *buf;
yading@10:     buf = av_malloc(width);
yading@10: 
yading@10:     src_m1 = src1;
yading@10:     memcpy(buf,src_m1,width);
yading@10:     src_0=&src_m1[src_wrap];
yading@10:     src_p1=&src_0[src_wrap];
yading@10:     src_p2=&src_p1[src_wrap];
yading@10:     for(y=0;y<(height-2);y+=2) {
yading@10:         deinterlace_line_inplace(buf,src_m1,src_0,src_p1,src_p2,width);
yading@10:         src_m1 = src_p1;
yading@10:         src_0 = src_p2;
yading@10:         src_p1 += 2*src_wrap;
yading@10:         src_p2 += 2*src_wrap;
yading@10:     }
yading@10:     /* do last line */
yading@10:     deinterlace_line_inplace(buf,src_m1,src_0,src_0,src_0,width);
yading@10:     av_free(buf);
yading@10: }
yading@10: 
yading@10: int avpicture_deinterlace(AVPicture *dst, const AVPicture *src,
yading@10:                           enum AVPixelFormat pix_fmt, int width, int height)
yading@10: {
yading@10:     int i;
yading@10: 
yading@10:     if (pix_fmt != AV_PIX_FMT_YUV420P &&
yading@10:         pix_fmt != AV_PIX_FMT_YUVJ420P &&
yading@10:         pix_fmt != AV_PIX_FMT_YUV422P &&
yading@10:         pix_fmt != AV_PIX_FMT_YUVJ422P &&
yading@10:         pix_fmt != AV_PIX_FMT_YUV444P &&
yading@10:         pix_fmt != AV_PIX_FMT_YUV411P &&
yading@10:         pix_fmt != AV_PIX_FMT_GRAY8)
yading@10:         return -1;
yading@10:     if ((width & 3) != 0 || (height & 3) != 0)
yading@10:         return -1;
yading@10: 
yading@10:     for(i=0;i<3;i++) {
yading@10:         if (i == 1) {
yading@10:             switch(pix_fmt) {
yading@10:             case AV_PIX_FMT_YUVJ420P:
yading@10:             case AV_PIX_FMT_YUV420P:
yading@10:                 width >>= 1;
yading@10:                 height >>= 1;
yading@10:                 break;
yading@10:             case AV_PIX_FMT_YUV422P:
yading@10:             case AV_PIX_FMT_YUVJ422P:
yading@10:                 width >>= 1;
yading@10:                 break;
yading@10:             case AV_PIX_FMT_YUV411P:
yading@10:                 width >>= 2;
yading@10:                 break;
yading@10:             default:
yading@10:                 break;
yading@10:             }
yading@10:             if (pix_fmt == AV_PIX_FMT_GRAY8) {
yading@10:                 break;
yading@10:             }
yading@10:         }
yading@10:         if (src == dst) {
yading@10:             deinterlace_bottom_field_inplace(dst->data[i], dst->linesize[i],
yading@10:                                  width, height);
yading@10:         } else {
yading@10:             deinterlace_bottom_field(dst->data[i],dst->linesize[i],
yading@10:                                         src->data[i], src->linesize[i],
yading@10:                                         width, height);
yading@10:         }
yading@10:     }
yading@10:     emms_c();
yading@10:     return 0;
yading@10: }
yading@10: 
yading@10: #endif /* FF_API_DEINTERLACE */
yading@10: 
yading@10: #ifdef TEST
yading@10: 
yading@10: int main(void){
yading@10:     int i;
yading@10:     int err=0;
yading@10:     int skip = 0;
yading@10: 
yading@10:     for (i=0; i<AV_PIX_FMT_NB*2; i++) {
yading@10:         AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(i);
yading@10:         if(!desc || !desc->name) {
yading@10:             skip ++;
yading@10:             continue;
yading@10:         }
yading@10:         if (skip) {
yading@10:             av_log(NULL, AV_LOG_INFO, "%3d unused pixel format values\n", skip);
yading@10:             skip = 0;
yading@10:         }
yading@10:         av_log(NULL, AV_LOG_INFO, "pix fmt %s yuv_plan:%d avg_bpp:%d colortype:%d\n", desc->name, is_yuv_planar(desc), av_get_padded_bits_per_pixel(desc), get_color_type(desc));
yading@10:         if ((!(desc->flags & PIX_FMT_ALPHA)) != (desc->nb_components != 2 && desc->nb_components != 4)) {
yading@10:             av_log(NULL, AV_LOG_ERROR, "Alpha flag mismatch\n");
yading@10:             err = 1;
yading@10:         }
yading@10:     }
yading@10:     return err;
yading@10: }
yading@10: 
yading@10: #endif