[deb_ffmpeg.git] / ffmpeg / libavcodec / imgconvert.c

/*
 * Misc image conversion routines
 * Copyright (c) 2001, 2002, 2003 Fabrice Bellard
 *
 * This file is part of FFmpeg.
 *
 * FFmpeg is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2.1 of the License, or (at your option) any later version.
 *
 * FFmpeg is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
 * License along with FFmpeg; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
 */

/**
 * @file
 * misc image conversion routines
 */

/* TODO:
 * - write 'ffimg' program to test all the image related stuff
 * - move all api to slice based system
 * - integrate deinterlacing, postprocessing and scaling in the conversion process
 */

#include "avcodec.h"
#include "imgconvert.h"
#include "internal.h"
#include "mathops.h"
#include "libavutil/avassert.h"
#include "libavutil/colorspace.h"
#include "libavutil/common.h"
#include "libavutil/pixdesc.h"
#include "libavutil/imgutils.h"

void avcodec_get_chroma_sub_sample(enum AVPixelFormat pix_fmt, int *h_shift, int *v_shift)
{
    const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(pix_fmt);
    av_assert0(desc);
    *h_shift = desc->log2_chroma_w;
    *v_shift = desc->log2_chroma_h;
}

int avcodec_get_pix_fmt_loss(enum AVPixelFormat dst_pix_fmt,
                             enum AVPixelFormat src_pix_fmt,
                             int has_alpha)
{
    return av_get_pix_fmt_loss(dst_pix_fmt, src_pix_fmt, has_alpha);
}

enum AVPixelFormat avcodec_find_best_pix_fmt_of_2(enum AVPixelFormat dst_pix_fmt1, enum AVPixelFormat dst_pix_fmt2,
                                            enum AVPixelFormat src_pix_fmt, int has_alpha, int *loss_ptr)
{
    return av_find_best_pix_fmt_of_2(dst_pix_fmt1, dst_pix_fmt2, src_pix_fmt, has_alpha, loss_ptr);
}

#if AV_HAVE_INCOMPATIBLE_LIBAV_ABI
enum AVPixelFormat avcodec_find_best_pix_fmt2(const enum AVPixelFormat *pix_fmt_list,
                                            enum AVPixelFormat src_pix_fmt,
                                            int has_alpha, int *loss_ptr){
    return avcodec_find_best_pix_fmt_of_list(pix_fmt_list, src_pix_fmt, has_alpha, loss_ptr);
}
#else
enum AVPixelFormat avcodec_find_best_pix_fmt2(enum AVPixelFormat dst_pix_fmt1, enum AVPixelFormat dst_pix_fmt2,
                                            enum AVPixelFormat src_pix_fmt, int has_alpha, int *loss_ptr)
{
    return avcodec_find_best_pix_fmt_of_2(dst_pix_fmt1, dst_pix_fmt2, src_pix_fmt, has_alpha, loss_ptr);
}
#endif

enum AVPixelFormat avcodec_find_best_pix_fmt_of_list(const enum AVPixelFormat *pix_fmt_list,
                                            enum AVPixelFormat src_pix_fmt,
                                            int has_alpha, int *loss_ptr){
    int i;

    enum AVPixelFormat best = AV_PIX_FMT_NONE;

    for(i=0; pix_fmt_list[i] != AV_PIX_FMT_NONE; i++)
        best = avcodec_find_best_pix_fmt_of_2(best, pix_fmt_list[i], src_pix_fmt, has_alpha, loss_ptr);

    return best;
}

/* 2x2 -> 1x1 */
void ff_shrink22(uint8_t *dst, int dst_wrap,
                     const uint8_t *src, int src_wrap,
                     int width, int height)
{
    int w;
    const uint8_t *s1, *s2;
    uint8_t *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) >> 2;
            d[1] = (s1[2] + s1[3] + s2[2] + s2[3] + 2) >> 2;
            d[2] = (s1[4] + s1[5] + s2[4] + s2[5] + 2) >> 2;
            d[3] = (s1[6] + s1[7] + s2[6] + s2[7] + 2) >> 2;
            s1 += 8;
            s2 += 8;
            d += 4;
        }
        for(;w > 0; w--) {
            d[0] = (s1[0] + s1[1] + s2[0] + s2[1] + 2) >> 2;
            s1 += 2;
            s2 += 2;
            d++;
        }
        src += 2 * src_wrap;
        dst += dst_wrap;
    }
}

/* 4x4 -> 1x1 */
void ff_shrink44(uint8_t *dst, int dst_wrap,
                     const uint8_t *src, int src_wrap,
                     int width, int height)
{
    int w;
    const uint8_t *s1, *s2, *s3, *s4;
    uint8_t *d;

    for(;height > 0; height--) {
        s1 = src;
        s2 = s1 + src_wrap;
        s3 = s2 + src_wrap;
        s4 = s3 + src_wrap;
        d = dst;
        for(w = width;w > 0; w--) {
            d[0] = (s1[0] + s1[1] + s1[2] + s1[3] +
                    s2[0] + s2[1] + s2[2] + s2[3] +
                    s3[0] + s3[1] + s3[2] + s3[3] +
                    s4[0] + s4[1] + s4[2] + s4[3] + 8) >> 4;
            s1 += 4;
            s2 += 4;
            s3 += 4;
            s4 += 4;
            d++;
        }
        src += 4 * src_wrap;
        dst += dst_wrap;
    }
}

/* 8x8 -> 1x1 */
void ff_shrink88(uint8_t *dst, int dst_wrap,
                     const uint8_t *src, int src_wrap,
                     int width, int height)
{
    int w, i;

    for(;height > 0; height--) {
        for(w = width;w > 0; w--) {
            int tmp=0;
            for(i=0; i<8; i++){
                tmp += src[0] + src[1] + src[2] + src[3] + src[4] + src[5] + src[6] + src[7];
                src += src_wrap;
            }
            *(dst++) = (tmp + 32)>>6;
            src += 8 - 8*src_wrap;
        }
        src += 8*src_wrap - 8*width;
        dst += dst_wrap - width;
    }
}

/* return true if yuv planar */
static inline int is_yuv_planar(const AVPixFmtDescriptor *desc)
{
    int i;
    int planes[4] = { 0 };

    if (     desc->flags & AV_PIX_FMT_FLAG_RGB
        || !(desc->flags & AV_PIX_FMT_FLAG_PLANAR))
        return 0;

    /* set the used planes */
    for (i = 0; i < desc->nb_components; i++)
        planes[desc->comp[i].plane] = 1;

    /* if there is an unused plane, the format is not planar */
    for (i = 0; i < desc->nb_components; i++)
        if (!planes[i])
            return 0;
    return 1;
}

int av_picture_crop(AVPicture *dst, const AVPicture *src,
                    enum AVPixelFormat pix_fmt, int top_band, int left_band)
{
    const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(pix_fmt);
    int y_shift;
    int x_shift;

    if (pix_fmt < 0 || pix_fmt >= AV_PIX_FMT_NB)
        return -1;

    y_shift = desc->log2_chroma_h;
    x_shift = desc->log2_chroma_w;

    if (is_yuv_planar(desc)) {
    dst->data[0] = src->data[0] + (top_band * src->linesize[0]) + left_band;
    dst->data[1] = src->data[1] + ((top_band >> y_shift) * src->linesize[1]) + (left_band >> x_shift);
    dst->data[2] = src->data[2] + ((top_band >> y_shift) * src->linesize[2]) + (left_band >> x_shift);
    } else{
        if(top_band % (1<<y_shift) || left_band % (1<<x_shift))
            return -1;
        if(left_band) //FIXME add support for this too
            return -1;
        dst->data[0] = src->data[0] + (top_band * src->linesize[0]) + left_band;
    }

    dst->linesize[0] = src->linesize[0];
    dst->linesize[1] = src->linesize[1];
    dst->linesize[2] = src->linesize[2];
    return 0;
}

int av_picture_pad(AVPicture *dst, const AVPicture *src, int height, int width,
                   enum AVPixelFormat pix_fmt, int padtop, int padbottom, int padleft, int padright,
            int *color)
{
    const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(pix_fmt);
    uint8_t *optr;
    int y_shift;
    int x_shift;
    int yheight;
    int i, y;

    if (pix_fmt < 0 || pix_fmt >= AV_PIX_FMT_NB ||
        !is_yuv_planar(desc)) return -1;

    for (i = 0; i < 3; i++) {
        x_shift = i ? desc->log2_chroma_w : 0;
        y_shift = i ? desc->log2_chroma_h : 0;

        if (padtop || padleft) {
            memset(dst->data[i], color[i],
                dst->linesize[i] * (padtop >> y_shift) + (padleft >> x_shift));
        }

        if (padleft || padright) {
            optr = dst->data[i] + dst->linesize[i] * (padtop >> y_shift) +
                (dst->linesize[i] - (padright >> x_shift));
            yheight = (height - 1 - (padtop + padbottom)) >> y_shift;
            for (y = 0; y < yheight; y++) {
                memset(optr, color[i], (padleft + padright) >> x_shift);
                optr += dst->linesize[i];
            }
        }

        if (src) { /* first line */
            uint8_t *iptr = src->data[i];
            optr = dst->data[i] + dst->linesize[i] * (padtop >> y_shift) +
                    (padleft >> x_shift);
            memcpy(optr, iptr, (width - padleft - padright) >> x_shift);
            iptr += src->linesize[i];
            optr = dst->data[i] + dst->linesize[i] * (padtop >> y_shift) +
                (dst->linesize[i] - (padright >> x_shift));
            yheight = (height - 1 - (padtop + padbottom)) >> y_shift;
            for (y = 0; y < yheight; y++) {
                memset(optr, color[i], (padleft + padright) >> x_shift);
                memcpy(optr + ((padleft + padright) >> x_shift), iptr,
                       (width - padleft - padright) >> x_shift);
                iptr += src->linesize[i];
                optr += dst->linesize[i];
            }
        }

        if (padbottom || padright) {
            optr = dst->data[i] + dst->linesize[i] *
                ((height - padbottom) >> y_shift) - (padright >> x_shift);
            memset(optr, color[i],dst->linesize[i] *
                (padbottom >> y_shift) + (padright >> x_shift));
        }
    }
    return 0;
}

#if FF_API_DEINTERLACE

#if HAVE_MMX_EXTERNAL
#define deinterlace_line_inplace ff_deinterlace_line_inplace_mmx
#define deinterlace_line         ff_deinterlace_line_mmx
#else
#define deinterlace_line_inplace deinterlace_line_inplace_c
#define deinterlace_line         deinterlace_line_c

/* filter parameters: [-1 4 2 4 -1] // 8 */
static void deinterlace_line_c(uint8_t *dst,
                             const uint8_t *lum_m4, const uint8_t *lum_m3,
                             const uint8_t *lum_m2, const uint8_t *lum_m1,
                             const uint8_t *lum,
                             int size)
{
    const uint8_t *cm = ff_crop_tab + MAX_NEG_CROP;
    int sum;

    for(;size > 0;size--) {
        sum = -lum_m4[0];
        sum += lum_m3[0] << 2;
        sum += lum_m2[0] << 1;
        sum += lum_m1[0] << 2;
        sum += -lum[0];
        dst[0] = cm[(sum + 4) >> 3];
        lum_m4++;
        lum_m3++;
        lum_m2++;
        lum_m1++;
        lum++;
        dst++;
    }
}

static void deinterlace_line_inplace_c(uint8_t *lum_m4, uint8_t *lum_m3,
                                       uint8_t *lum_m2, uint8_t *lum_m1,
                                       uint8_t *lum, int size)
{
    const uint8_t *cm = ff_crop_tab + MAX_NEG_CROP;
    int sum;

    for(;size > 0;size--) {
        sum = -lum_m4[0];
        sum += lum_m3[0] << 2;
        sum += lum_m2[0] << 1;
        lum_m4[0]=lum_m2[0];
        sum += lum_m1[0] << 2;
        sum += -lum[0];
        lum_m2[0] = cm[(sum + 4) >> 3];
        lum_m4++;
        lum_m3++;
        lum_m2++;
        lum_m1++;
        lum++;
    }
}
#endif /* !HAVE_MMX_EXTERNAL */

/* 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_t *dst, int dst_wrap,
                                    const uint8_t *src1, int src_wrap,
                                    int width, int height)
{
    const uint8_t *src_m2, *src_m1, *src_0, *src_p1, *src_p2;
    int y;

    src_m2 = src1;
    src_m1 = src1;
    src_0=&src_m1[src_wrap];
    src_p1=&src_0[src_wrap];
    src_p2=&src_p1[src_wrap];
    for(y=0;y<(height-2);y+=2) {
        memcpy(dst,src_m1,width);
        dst += dst_wrap;
        deinterlace_line(dst,src_m2,src_m1,src_0,src_p1,src_p2,width);
        src_m2 = src_0;
        src_m1 = src_p1;
        src_0 = src_p2;
        src_p1 += 2*src_wrap;
        src_p2 += 2*src_wrap;
        dst += dst_wrap;
    }
    memcpy(dst,src_m1,width);
    dst += dst_wrap;
    /* do last line */
    deinterlace_line(dst,src_m2,src_m1,src_0,src_0,src_0,width);
}

static void deinterlace_bottom_field_inplace(uint8_t *src1, int src_wrap,
                                             int width, int height)
{
    uint8_t *src_m1, *src_0, *src_p1, *src_p2;
    int y;
    uint8_t *buf;
    buf = av_malloc(width);

    src_m1 = src1;
    memcpy(buf,src_m1,width);
    src_0=&src_m1[src_wrap];
    src_p1=&src_0[src_wrap];
    src_p2=&src_p1[src_wrap];
    for(y=0;y<(height-2);y+=2) {
        deinterlace_line_inplace(buf,src_m1,src_0,src_p1,src_p2,width);
        src_m1 = src_p1;
        src_0 = src_p2;
        src_p1 += 2*src_wrap;
        src_p2 += 2*src_wrap;
    }
    /* do last line */
    deinterlace_line_inplace(buf,src_m1,src_0,src_0,src_0,width);
    av_free(buf);
}

int avpicture_deinterlace(AVPicture *dst, const AVPicture *src,
                          enum AVPixelFormat pix_fmt, int width, int height)
{
    int i;

    if (pix_fmt != AV_PIX_FMT_YUV420P &&
        pix_fmt != AV_PIX_FMT_YUVJ420P &&
        pix_fmt != AV_PIX_FMT_YUV422P &&
        pix_fmt != AV_PIX_FMT_YUVJ422P &&
        pix_fmt != AV_PIX_FMT_YUV444P &&
        pix_fmt != AV_PIX_FMT_YUV411P &&
        pix_fmt != AV_PIX_FMT_GRAY8)
        return -1;
    if ((width & 3) != 0 || (height & 3) != 0)
        return -1;

    for(i=0;i<3;i++) {
        if (i == 1) {
            switch(pix_fmt) {
            case AV_PIX_FMT_YUVJ420P:
            case AV_PIX_FMT_YUV420P:
                width >>= 1;
                height >>= 1;
                break;
            case AV_PIX_FMT_YUV422P:
            case AV_PIX_FMT_YUVJ422P:
                width >>= 1;
                break;
            case AV_PIX_FMT_YUV411P:
                width >>= 2;
                break;
            default:
                break;
            }
            if (pix_fmt == AV_PIX_FMT_GRAY8) {
                break;
            }
        }
        if (src == dst) {
            deinterlace_bottom_field_inplace(dst->data[i], dst->linesize[i],
                                 width, height);
        } else {
            deinterlace_bottom_field(dst->data[i],dst->linesize[i],
                                        src->data[i], src->linesize[i],
                                        width, height);
        }
    }
    emms_c();
    return 0;
}

#endif /* FF_API_DEINTERLACE */

#ifdef TEST

int main(void){
    int i;
    int err=0;
    int skip = 0;

    for (i=0; i<AV_PIX_FMT_NB*2; i++) {
        const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(i);
        if(!desc || !desc->name) {
            skip ++;
            continue;
        }
        if (skip) {
            av_log(NULL, AV_LOG_INFO, "%3d unused pixel format values\n", skip);
            skip = 0;
        }
        av_log(NULL, AV_LOG_INFO, "pix fmt %s yuv_plan:%d avg_bpp:%d\n", desc->name, is_yuv_planar(desc), av_get_padded_bits_per_pixel(desc));
        if ((!(desc->flags & AV_PIX_FMT_FLAG_ALPHA)) != (desc->nb_components != 2 && desc->nb_components != 4)) {
            av_log(NULL, AV_LOG_ERROR, "Alpha flag mismatch\n");
            err = 1;
        }
    }
    return err;
}

#endif
Commit	Line	Data
2ba45a60 DM	1	/*
	2	* Misc image conversion routines
	3	* Copyright (c) 2001, 2002, 2003 Fabrice Bellard
	4	*
	5	* This file is part of FFmpeg.
	6	*
	7	* FFmpeg is free software; you can redistribute it and/or
	8	* modify it under the terms of the GNU Lesser General Public
	9	* License as published by the Free Software Foundation; either
	10	* version 2.1 of the License, or (at your option) any later version.
	11	*
	12	* FFmpeg is distributed in the hope that it will be useful,
	13	* but WITHOUT ANY WARRANTY; without even the implied warranty of
	14	* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
	15	* Lesser General Public License for more details.
	16	*
	17	* You should have received a copy of the GNU Lesser General Public
	18	* License along with FFmpeg; if not, write to the Free Software
	19	* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
	20	*/
	21
	22	/**
	23	* @file
	24	* misc image conversion routines
	25	*/
	26
	27	/* TODO:
	28	* - write 'ffimg' program to test all the image related stuff
	29	* - move all api to slice based system
	30	* - integrate deinterlacing, postprocessing and scaling in the conversion process
	31	*/
	32
	33	#include "avcodec.h"
	34	#include "imgconvert.h"
	35	#include "internal.h"
	36	#include "mathops.h"
	37	#include "libavutil/avassert.h"
	38	#include "libavutil/colorspace.h"
	39	#include "libavutil/common.h"
	40	#include "libavutil/pixdesc.h"
	41	#include "libavutil/imgutils.h"
	42
	43	void avcodec_get_chroma_sub_sample(enum AVPixelFormat pix_fmt, int h_shift, int v_shift)
	44	{
	45	const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(pix_fmt);
	46	av_assert0(desc);
	47	*h_shift = desc->log2_chroma_w;
	48	*v_shift = desc->log2_chroma_h;
	49	}
	50
	51	int avcodec_get_pix_fmt_loss(enum AVPixelFormat dst_pix_fmt,
	52	enum AVPixelFormat src_pix_fmt,
	53	int has_alpha)
	54	{
	55	return av_get_pix_fmt_loss(dst_pix_fmt, src_pix_fmt, has_alpha);
	56	}
	57
	58	enum AVPixelFormat avcodec_find_best_pix_fmt_of_2(enum AVPixelFormat dst_pix_fmt1, enum AVPixelFormat dst_pix_fmt2,
	59	enum AVPixelFormat src_pix_fmt, int has_alpha, int *loss_ptr)
	60	{
	61	return av_find_best_pix_fmt_of_2(dst_pix_fmt1, dst_pix_fmt2, src_pix_fmt, has_alpha, loss_ptr);
	62	}
	63
	64	#if AV_HAVE_INCOMPATIBLE_LIBAV_ABI
65	enum AVPixelFormat avcodec_find_best_pix_fmt2(const enum AVPixelFormat *pix_fmt_list,
66	enum AVPixelFormat src_pix_fmt,
67	int has_alpha, int *loss_ptr){
68	return avcodec_find_best_pix_fmt_of_list(pix_fmt_list, src_pix_fmt, has_alpha, loss_ptr);
69	}
70	#else
71	enum AVPixelFormat avcodec_find_best_pix_fmt2(enum AVPixelFormat dst_pix_fmt1, enum AVPixelFormat dst_pix_fmt2,
72	enum AVPixelFormat src_pix_fmt, int has_alpha, int *loss_ptr)
73	{
74	return avcodec_find_best_pix_fmt_of_2(dst_pix_fmt1, dst_pix_fmt2, src_pix_fmt, has_alpha, loss_ptr);
75	}
76	#endif
77
78	enum AVPixelFormat avcodec_find_best_pix_fmt_of_list(const enum AVPixelFormat *pix_fmt_list,
79	enum AVPixelFormat src_pix_fmt,
80	int has_alpha, int *loss_ptr){
81	int i;
82
83	enum AVPixelFormat best = AV_PIX_FMT_NONE;
84
85	for(i=0; pix_fmt_list[i] != AV_PIX_FMT_NONE; i++)
86	best = avcodec_find_best_pix_fmt_of_2(best, pix_fmt_list[i], src_pix_fmt, has_alpha, loss_ptr);
87
88	return best;
89	}
90
91	/* 2x2 -> 1x1 */
92	void ff_shrink22(uint8_t *dst, int dst_wrap,
93	const uint8_t *src, int src_wrap,
94	int width, int height)
95	{
96	int w;
97	const uint8_t s1, s2;
98	uint8_t *d;
99
100	for(;height > 0; height--) {
101	s1 = src;
102	s2 = s1 + src_wrap;
103	d = dst;
104	for(w = width;w >= 4; w-=4) {
105	d[0] = (s1[0] + s1[1] + s2[0] + s2[1] + 2) >> 2;
106	d[1] = (s1[2] + s1[3] + s2[2] + s2[3] + 2) >> 2;
107	d[2] = (s1[4] + s1[5] + s2[4] + s2[5] + 2) >> 2;
108	d[3] = (s1[6] + s1[7] + s2[6] + s2[7] + 2) >> 2;
109	s1 += 8;
110	s2 += 8;
111	d += 4;
112	}
113	for(;w > 0; w--) {
114	d[0] = (s1[0] + s1[1] + s2[0] + s2[1] + 2) >> 2;
115	s1 += 2;
116	s2 += 2;
117	d++;
118	}
119	src += 2 * src_wrap;
120	dst += dst_wrap;
121	}
122	}
123
124	/* 4x4 -> 1x1 */
125	void ff_shrink44(uint8_t *dst, int dst_wrap,
126	const uint8_t *src, int src_wrap,
127	int width, int height)
128	{
129	int w;
130	const uint8_t s1, s2, s3, s4;
131	uint8_t *d;
132
133	for(;height > 0; height--) {
134	s1 = src;
135	s2 = s1 + src_wrap;
136	s3 = s2 + src_wrap;
137	s4 = s3 + src_wrap;
138	d = dst;
139	for(w = width;w > 0; w--) {
140	d[0] = (s1[0] + s1[1] + s1[2] + s1[3] +
141	s2[0] + s2[1] + s2[2] + s2[3] +
142	s3[0] + s3[1] + s3[2] + s3[3] +
143	s4[0] + s4[1] + s4[2] + s4[3] + 8) >> 4;
144	s1 += 4;
145	s2 += 4;
146	s3 += 4;
147	s4 += 4;
148	d++;
149	}
150	src += 4 * src_wrap;
151	dst += dst_wrap;
152	}
153	}
154
155	/* 8x8 -> 1x1 */
156	void ff_shrink88(uint8_t *dst, int dst_wrap,
157	const uint8_t *src, int src_wrap,
158	int width, int height)
159	{
160	int w, i;
161
162	for(;height > 0; height--) {
163	for(w = width;w > 0; w--) {
164	int tmp=0;
165	for(i=0; i<8; i++){
166	tmp += src[0] + src[1] + src[2] + src[3] + src[4] + src[5] + src[6] + src[7];
167	src += src_wrap;
168	}
169	*(dst++) = (tmp + 32)>>6;
170	src += 8 - 8*src_wrap;
171	}
172	src += 8src_wrap - 8width;
173	dst += dst_wrap - width;
174	}
175	}
176
177	/* return true if yuv planar */
178	static inline int is_yuv_planar(const AVPixFmtDescriptor *desc)
179	{
180	int i;
181	int planes[4] = { 0 };
182
183	if ( desc->flags & AV_PIX_FMT_FLAG_RGB
184	\|\| !(desc->flags & AV_PIX_FMT_FLAG_PLANAR))
185	return 0;
186
187	/* set the used planes */
188	for (i = 0; i < desc->nb_components; i++)
189	planes[desc->comp[i].plane] = 1;
190
191	/* if there is an unused plane, the format is not planar */
192	for (i = 0; i < desc->nb_components; i++)
193	if (!planes[i])
194	return 0;
195	return 1;
196	}
197
198	int av_picture_crop(AVPicture dst, const AVPicture src,
199	enum AVPixelFormat pix_fmt, int top_band, int left_band)
200	{
201	const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(pix_fmt);
202	int y_shift;
203	int x_shift;
204
205	if (pix_fmt < 0 \|\| pix_fmt >= AV_PIX_FMT_NB)
206	return -1;
207
208	y_shift = desc->log2_chroma_h;
209	x_shift = desc->log2_chroma_w;
210
211	if (is_yuv_planar(desc)) {
212	dst->data[0] = src->data[0] + (top_band * src->linesize[0]) + left_band;
213	dst->data[1] = src->data[1] + ((top_band >> y_shift) * src->linesize[1]) + (left_band >> x_shift);
214	dst->data[2] = src->data[2] + ((top_band >> y_shift) * src->linesize[2]) + (left_band >> x_shift);
215	} else{
216	if(top_band % (1<<y_shift) \|\| left_band % (1<<x_shift))
217	return -1;
218	if(left_band) //FIXME add support for this too
219	return -1;
220	dst->data[0] = src->data[0] + (top_band * src->linesize[0]) + left_band;
221	}
222
223	dst->linesize[0] = src->linesize[0];
224	dst->linesize[1] = src->linesize[1];
225	dst->linesize[2] = src->linesize[2];
226	return 0;
227	}
228
229	int av_picture_pad(AVPicture dst, const AVPicture src, int height, int width,
230	enum AVPixelFormat pix_fmt, int padtop, int padbottom, int padleft, int padright,
231	int *color)
232	{
233	const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(pix_fmt);
234	uint8_t *optr;
235	int y_shift;
236	int x_shift;
237	int yheight;
238	int i, y;
239
240	if (pix_fmt < 0 \|\| pix_fmt >= AV_PIX_FMT_NB \|\|
241	!is_yuv_planar(desc)) return -1;
242
243	for (i = 0; i < 3; i++) {
244	x_shift = i ? desc->log2_chroma_w : 0;
245	y_shift = i ? desc->log2_chroma_h : 0;
246
247	if (padtop \|\| padleft) {
248	memset(dst->data[i], color[i],
249	dst->linesize[i] * (padtop >> y_shift) + (padleft >> x_shift));
250	}
251
252	if (padleft \|\| padright) {
253	optr = dst->data[i] + dst->linesize[i] * (padtop >> y_shift) +
254	(dst->linesize[i] - (padright >> x_shift));
255	yheight = (height - 1 - (padtop + padbottom)) >> y_shift;
256	for (y = 0; y < yheight; y++) {
257	memset(optr, color[i], (padleft + padright) >> x_shift);
258	optr += dst->linesize[i];
259	}
260	}
261
262	if (src) { /* first line */
263	uint8_t *iptr = src->data[i];
264	optr = dst->data[i] + dst->linesize[i] * (padtop >> y_shift) +
265	(padleft >> x_shift);
266	memcpy(optr, iptr, (width - padleft - padright) >> x_shift);
267	iptr += src->linesize[i];
268	optr = dst->data[i] + dst->linesize[i] * (padtop >> y_shift) +
269	(dst->linesize[i] - (padright >> x_shift));
270	yheight = (height - 1 - (padtop + padbottom)) >> y_shift;
271	for (y = 0; y < yheight; y++) {
272	memset(optr, color[i], (padleft + padright) >> x_shift);
273	memcpy(optr + ((padleft + padright) >> x_shift), iptr,
274	(width - padleft - padright) >> x_shift);
275	iptr += src->linesize[i];
276	optr += dst->linesize[i];
277	}
278	}
279
280	if (padbottom \|\| padright) {
281	optr = dst->data[i] + dst->linesize[i] *
282	((height - padbottom) >> y_shift) - (padright >> x_shift);
283	memset(optr, color[i],dst->linesize[i] *
284	(padbottom >> y_shift) + (padright >> x_shift));
285	}
286	}
287	return 0;
288	}
289
290	#if FF_API_DEINTERLACE
291
292	#if HAVE_MMX_EXTERNAL
293	#define deinterlace_line_inplace ff_deinterlace_line_inplace_mmx
294	#define deinterlace_line ff_deinterlace_line_mmx
295	#else
296	#define deinterlace_line_inplace deinterlace_line_inplace_c
297	#define deinterlace_line deinterlace_line_c
298
299	/* filter parameters: [-1 4 2 4 -1] // 8 */
300	static void deinterlace_line_c(uint8_t *dst,
301	const uint8_t lum_m4, const uint8_t lum_m3,
302	const uint8_t lum_m2, const uint8_t lum_m1,
303	const uint8_t *lum,
304	int size)
305	{
306	const uint8_t *cm = ff_crop_tab + MAX_NEG_CROP;
307	int sum;
308
309	for(;size > 0;size--) {
310	sum = -lum_m4[0];
311	sum += lum_m3[0] << 2;
312	sum += lum_m2[0] << 1;
313	sum += lum_m1[0] << 2;
314	sum += -lum[0];
315	dst[0] = cm[(sum + 4) >> 3];
316	lum_m4++;
317	lum_m3++;
318	lum_m2++;
319	lum_m1++;
320	lum++;
321	dst++;
322	}
323	}
324
325	static void deinterlace_line_inplace_c(uint8_t lum_m4, uint8_t lum_m3,
326	uint8_t lum_m2, uint8_t lum_m1,
327	uint8_t *lum, int size)
328	{
329	const uint8_t *cm = ff_crop_tab + MAX_NEG_CROP;
330	int sum;
331
332	for(;size > 0;size--) {
333	sum = -lum_m4[0];
334	sum += lum_m3[0] << 2;
335	sum += lum_m2[0] << 1;
336	lum_m4[0]=lum_m2[0];
337	sum += lum_m1[0] << 2;
338	sum += -lum[0];
339	lum_m2[0] = cm[(sum + 4) >> 3];
340	lum_m4++;
341	lum_m3++;
342	lum_m2++;
343	lum_m1++;
344	lum++;
345	}
346	}
347	#endif /* !HAVE_MMX_EXTERNAL */
348
349	/* deinterlacing : 2 temporal taps, 3 spatial taps linear filter. The
350	top field is copied as is, but the bottom field is deinterlaced
351	against the top field. */
352	static void deinterlace_bottom_field(uint8_t *dst, int dst_wrap,
353	const uint8_t *src1, int src_wrap,
354	int width, int height)
355	{
356	const uint8_t src_m2, src_m1, src_0, src_p1, *src_p2;
357	int y;
358
359	src_m2 = src1;
360	src_m1 = src1;
361	src_0=&src_m1[src_wrap];
362	src_p1=&src_0[src_wrap];
363	src_p2=&src_p1[src_wrap];
364	for(y=0;y<(height-2);y+=2) {
365	memcpy(dst,src_m1,width);
366	dst += dst_wrap;
367	deinterlace_line(dst,src_m2,src_m1,src_0,src_p1,src_p2,width);
368	src_m2 = src_0;
369	src_m1 = src_p1;
370	src_0 = src_p2;
371	src_p1 += 2*src_wrap;
372	src_p2 += 2*src_wrap;
373	dst += dst_wrap;
374	}
375	memcpy(dst,src_m1,width);
376	dst += dst_wrap;
377	/* do last line */
378	deinterlace_line(dst,src_m2,src_m1,src_0,src_0,src_0,width);
379	}
380
381	static void deinterlace_bottom_field_inplace(uint8_t *src1, int src_wrap,
382	int width, int height)
383	{
384	uint8_t src_m1, src_0, src_p1, src_p2;
385	int y;
386	uint8_t *buf;
387	buf = av_malloc(width);
388
389	src_m1 = src1;
390	memcpy(buf,src_m1,width);
391	src_0=&src_m1[src_wrap];
392	src_p1=&src_0[src_wrap];
393	src_p2=&src_p1[src_wrap];
394	for(y=0;y<(height-2);y+=2) {
395	deinterlace_line_inplace(buf,src_m1,src_0,src_p1,src_p2,width);
396	src_m1 = src_p1;
397	src_0 = src_p2;
398	src_p1 += 2*src_wrap;
399	src_p2 += 2*src_wrap;
400	}
401	/* do last line */
402	deinterlace_line_inplace(buf,src_m1,src_0,src_0,src_0,width);
403	av_free(buf);
404	}
405
406	int avpicture_deinterlace(AVPicture dst, const AVPicture src,
407	enum AVPixelFormat pix_fmt, int width, int height)
408	{
409	int i;
410
411	if (pix_fmt != AV_PIX_FMT_YUV420P &&
412	pix_fmt != AV_PIX_FMT_YUVJ420P &&
413	pix_fmt != AV_PIX_FMT_YUV422P &&
414	pix_fmt != AV_PIX_FMT_YUVJ422P &&
415	pix_fmt != AV_PIX_FMT_YUV444P &&
416	pix_fmt != AV_PIX_FMT_YUV411P &&
417	pix_fmt != AV_PIX_FMT_GRAY8)
418	return -1;
419	if ((width & 3) != 0 \|\| (height & 3) != 0)
420	return -1;
421
422	for(i=0;i<3;i++) {
423	if (i == 1) {
424	switch(pix_fmt) {
425	case AV_PIX_FMT_YUVJ420P:
426	case AV_PIX_FMT_YUV420P:
427	width >>= 1;
428	height >>= 1;
429	break;
430	case AV_PIX_FMT_YUV422P:
431	case AV_PIX_FMT_YUVJ422P:
432	width >>= 1;
433	break;
434	case AV_PIX_FMT_YUV411P:
435	width >>= 2;
436	break;
437	default:
438	break;
439	}
440	if (pix_fmt == AV_PIX_FMT_GRAY8) {
441	break;
442	}
443	}
444	if (src == dst) {
445	deinterlace_bottom_field_inplace(dst->data[i], dst->linesize[i],
446	width, height);
447	} else {
448	deinterlace_bottom_field(dst->data[i],dst->linesize[i],
449	src->data[i], src->linesize[i],
450	width, height);
451	}
452	}
453	emms_c();
454	return 0;
455	}
456
457	#endif /* FF_API_DEINTERLACE */
458
459	#ifdef TEST
460
461	int main(void){
462	int i;
463	int err=0;
464	int skip = 0;
465
466	for (i=0; i<AV_PIX_FMT_NB*2; i++) {
467	const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(i);
468	if(!desc \|\| !desc->name) {
469	skip ++;
470	continue;
471	}
472	if (skip) {
473	av_log(NULL, AV_LOG_INFO, "%3d unused pixel format values\n", skip);
474	skip = 0;
475	}
476	av_log(NULL, AV_LOG_INFO, "pix fmt %s yuv_plan:%d avg_bpp:%d\n", desc->name, is_yuv_planar(desc), av_get_padded_bits_per_pixel(desc));
477	if ((!(desc->flags & AV_PIX_FMT_FLAG_ALPHA)) != (desc->nb_components != 2 && desc->nb_components != 4)) {
478	av_log(NULL, AV_LOG_ERROR, "Alpha flag mismatch\n");
479	err = 1;
480	}
481	}
482	return err;
483	}
484
485	#endif