Correct edge MC for chroma
[ffmpeg.git] / libavcodec / imgconvert.c
1 /*
2  * Misc image convertion routines
3  * Copyright (c) 2001, 2002, 2003 Fabrice Bellard.
4  *
5  * This library is free software; you can redistribute it and/or
6  * modify it under the terms of the GNU Lesser General Public
7  * License as published by the Free Software Foundation; either
8  * version 2 of the License, or (at your option) any later version.
9  *
10  * This library is distributed in the hope that it will be useful,
11  * but WITHOUT ANY WARRANTY; without even the implied warranty of
12  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
13  * Lesser General Public License for more details.
14  *
15  * You should have received a copy of the GNU Lesser General Public
16  * License along with this library; if not, write to the Free Software
17  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
18  */
19
20 /**
21  * @file imgconvert.c
22  * Misc image convertion routines.
23  */
24
25 /* TODO:
26  * - write 'ffimg' program to test all the image related stuff
27  * - move all api to slice based system
28  * - integrate deinterlacing, postprocessing and scaling in the conversion process
29  */
30
31 #include "avcodec.h"
32 #include "dsputil.h"
33
34 #ifdef USE_FASTMEMCPY
35 #include "fastmemcpy.h"
36 #endif
37
38 #ifdef HAVE_MMX
39 #include "i386/mmx.h"
40 #endif
41
42 #define xglue(x, y) x ## y
43 #define glue(x, y) xglue(x, y)
44
45 #define FF_COLOR_RGB      0 /* RGB color space */
46 #define FF_COLOR_GRAY     1 /* gray color space */
47 #define FF_COLOR_YUV      2 /* YUV color space. 16 <= Y <= 235, 16 <= U, V <= 240 */
48 #define FF_COLOR_YUV_JPEG 3 /* YUV color space. 0 <= Y <= 255, 0 <= U, V <= 255 */
49
50 #define FF_PIXEL_PLANAR   0 /* each channel has one component in AVPicture */
51 #define FF_PIXEL_PACKED   1 /* only one components containing all the channels */
52 #define FF_PIXEL_PALETTE  2  /* one components containing indexes for a palette */
53
54 typedef struct PixFmtInfo {
55     const char *name;
56     uint8_t nb_channels;     /* number of channels (including alpha) */
57     uint8_t color_type;      /* color type (see FF_COLOR_xxx constants) */
58     uint8_t pixel_type;      /* pixel storage type (see FF_PIXEL_xxx constants) */
59     uint8_t is_alpha : 1;    /* true if alpha can be specified */
60     uint8_t x_chroma_shift;  /* X chroma subsampling factor is 2 ^ shift */
61     uint8_t y_chroma_shift;  /* Y chroma subsampling factor is 2 ^ shift */
62     uint8_t depth;           /* bit depth of the color components */
63 } PixFmtInfo;
64
65 /* this table gives more information about formats */
66 static const PixFmtInfo pix_fmt_info[PIX_FMT_NB] = {
67     /* YUV formats */
68     [PIX_FMT_YUV420P] = {
69         .name = "yuv420p",
70         .nb_channels = 3,
71         .color_type = FF_COLOR_YUV,
72         .pixel_type = FF_PIXEL_PLANAR,
73         .depth = 8,
74         .x_chroma_shift = 1, .y_chroma_shift = 1,
75     },
76     [PIX_FMT_YUV422P] = {
77         .name = "yuv422p",
78         .nb_channels = 3,
79         .color_type = FF_COLOR_YUV,
80         .pixel_type = FF_PIXEL_PLANAR,
81         .depth = 8,
82         .x_chroma_shift = 1, .y_chroma_shift = 0,
83     },
84     [PIX_FMT_YUV444P] = {
85         .name = "yuv444p",
86         .nb_channels = 3,
87         .color_type = FF_COLOR_YUV,
88         .pixel_type = FF_PIXEL_PLANAR,
89         .depth = 8,
90         .x_chroma_shift = 0, .y_chroma_shift = 0,
91     },
92     [PIX_FMT_YUV422] = {
93         .name = "yuv422",
94         .nb_channels = 1,
95         .color_type = FF_COLOR_YUV,
96         .pixel_type = FF_PIXEL_PACKED,
97         .depth = 8,
98         .x_chroma_shift = 1, .y_chroma_shift = 0,
99     },
100     [PIX_FMT_UYVY422] = {
101         .name = "uyvy422",
102         .nb_channels = 1,
103         .color_type = FF_COLOR_YUV,
104         .pixel_type = FF_PIXEL_PACKED,
105         .depth = 8,
106         .x_chroma_shift = 1, .y_chroma_shift = 0,
107     },
108     [PIX_FMT_YUV410P] = {
109         .name = "yuv410p",
110         .nb_channels = 3,
111         .color_type = FF_COLOR_YUV,
112         .pixel_type = FF_PIXEL_PLANAR,
113         .depth = 8,
114         .x_chroma_shift = 2, .y_chroma_shift = 2,
115     },
116     [PIX_FMT_YUV411P] = {
117         .name = "yuv411p",
118         .nb_channels = 3,
119         .color_type = FF_COLOR_YUV,
120         .pixel_type = FF_PIXEL_PLANAR,
121         .depth = 8,
122         .x_chroma_shift = 2, .y_chroma_shift = 0,
123     },
124
125     /* JPEG YUV */
126     [PIX_FMT_YUVJ420P] = {
127         .name = "yuvj420p",
128         .nb_channels = 3,
129         .color_type = FF_COLOR_YUV_JPEG,
130         .pixel_type = FF_PIXEL_PLANAR,
131         .depth = 8,
132         .x_chroma_shift = 1, .y_chroma_shift = 1,
133     },
134     [PIX_FMT_YUVJ422P] = {
135         .name = "yuvj422p",
136         .nb_channels = 3,
137         .color_type = FF_COLOR_YUV_JPEG,
138         .pixel_type = FF_PIXEL_PLANAR,
139         .depth = 8,
140         .x_chroma_shift = 1, .y_chroma_shift = 0,
141     },
142     [PIX_FMT_YUVJ444P] = {
143         .name = "yuvj444p",
144         .nb_channels = 3,
145         .color_type = FF_COLOR_YUV_JPEG,
146         .pixel_type = FF_PIXEL_PLANAR,
147         .depth = 8,
148         .x_chroma_shift = 0, .y_chroma_shift = 0,
149     },
150
151     /* RGB formats */
152     [PIX_FMT_RGB24] = {
153         .name = "rgb24",
154         .nb_channels = 3,
155         .color_type = FF_COLOR_RGB,
156         .pixel_type = FF_PIXEL_PACKED,
157         .depth = 8,
158         .x_chroma_shift = 0, .y_chroma_shift = 0,
159     },
160     [PIX_FMT_BGR24] = {
161         .name = "bgr24",
162         .nb_channels = 3,
163         .color_type = FF_COLOR_RGB,
164         .pixel_type = FF_PIXEL_PACKED,
165         .depth = 8,
166         .x_chroma_shift = 0, .y_chroma_shift = 0,
167     },
168     [PIX_FMT_RGBA32] = {
169         .name = "rgba32",
170         .nb_channels = 4, .is_alpha = 1,
171         .color_type = FF_COLOR_RGB,
172         .pixel_type = FF_PIXEL_PACKED,
173         .depth = 8,
174         .x_chroma_shift = 0, .y_chroma_shift = 0,
175     },
176     [PIX_FMT_RGB565] = {
177         .name = "rgb565",
178         .nb_channels = 3,
179         .color_type = FF_COLOR_RGB,
180         .pixel_type = FF_PIXEL_PACKED,
181         .depth = 5,
182         .x_chroma_shift = 0, .y_chroma_shift = 0,
183     },
184     [PIX_FMT_RGB555] = {
185         .name = "rgb555",
186         .nb_channels = 4, .is_alpha = 1,
187         .color_type = FF_COLOR_RGB,
188         .pixel_type = FF_PIXEL_PACKED,
189         .depth = 5,
190         .x_chroma_shift = 0, .y_chroma_shift = 0,
191     },
192
193     /* gray / mono formats */
194     [PIX_FMT_GRAY8] = {
195         .name = "gray",
196         .nb_channels = 1,
197         .color_type = FF_COLOR_GRAY,
198         .pixel_type = FF_PIXEL_PLANAR,
199         .depth = 8,
200     },
201     [PIX_FMT_MONOWHITE] = {
202         .name = "monow",
203         .nb_channels = 1,
204         .color_type = FF_COLOR_GRAY,
205         .pixel_type = FF_PIXEL_PLANAR,
206         .depth = 1,
207     },
208     [PIX_FMT_MONOBLACK] = {
209         .name = "monob",
210         .nb_channels = 1,
211         .color_type = FF_COLOR_GRAY,
212         .pixel_type = FF_PIXEL_PLANAR,
213         .depth = 1,
214     },
215
216     /* paletted formats */
217     [PIX_FMT_PAL8] = {
218         .name = "pal8",
219         .nb_channels = 4, .is_alpha = 1,
220         .color_type = FF_COLOR_RGB,
221         .pixel_type = FF_PIXEL_PALETTE,
222         .depth = 8,
223     },
224     [PIX_FMT_XVMC_MPEG2_MC] = {
225         .name = "xvmcmc",
226     },
227     [PIX_FMT_XVMC_MPEG2_IDCT] = {
228         .name = "xvmcidct",
229     },
230     [PIX_FMT_UYVY411] = {
231         .name = "uyvy411",
232         .nb_channels = 1,
233         .color_type = FF_COLOR_YUV,
234         .pixel_type = FF_PIXEL_PACKED,
235         .depth = 8,
236         .x_chroma_shift = 2, .y_chroma_shift = 0,
237     },
238 };
239
240 void avcodec_get_chroma_sub_sample(int pix_fmt, int *h_shift, int *v_shift)
241 {
242     *h_shift = pix_fmt_info[pix_fmt].x_chroma_shift;
243     *v_shift = pix_fmt_info[pix_fmt].y_chroma_shift;
244 }
245
246 const char *avcodec_get_pix_fmt_name(int pix_fmt)
247 {
248     if (pix_fmt < 0 || pix_fmt >= PIX_FMT_NB)
249         return "???";
250     else
251         return pix_fmt_info[pix_fmt].name;
252 }
253
254 enum PixelFormat avcodec_get_pix_fmt(const char* name)
255 {
256     int i;
257
258     for (i=0; i < PIX_FMT_NB; i++)
259          if (!strcmp(pix_fmt_info[i].name, name))
260              break;
261     return i;
262 }
263
264 /* Picture field are filled with 'ptr' addresses. Also return size */
265 int avpicture_fill(AVPicture *picture, uint8_t *ptr,
266                    int pix_fmt, int width, int height)
267 {
268     int size, w2, h2, size2;
269     const PixFmtInfo *pinfo;
270
271     if(avcodec_check_dimensions(NULL, width, height))
272         goto fail;
273
274     pinfo = &pix_fmt_info[pix_fmt];
275     size = width * height;
276     switch(pix_fmt) {
277     case PIX_FMT_YUV420P:
278     case PIX_FMT_YUV422P:
279     case PIX_FMT_YUV444P:
280     case PIX_FMT_YUV410P:
281     case PIX_FMT_YUV411P:
282     case PIX_FMT_YUVJ420P:
283     case PIX_FMT_YUVJ422P:
284     case PIX_FMT_YUVJ444P:
285         w2 = (width + (1 << pinfo->x_chroma_shift) - 1) >> pinfo->x_chroma_shift;
286         h2 = (height + (1 << pinfo->y_chroma_shift) - 1) >> pinfo->y_chroma_shift;
287         size2 = w2 * h2;
288         picture->data[0] = ptr;
289         picture->data[1] = picture->data[0] + size;
290         picture->data[2] = picture->data[1] + size2;
291         picture->linesize[0] = width;
292         picture->linesize[1] = w2;
293         picture->linesize[2] = w2;
294         return size + 2 * size2;
295     case PIX_FMT_RGB24:
296     case PIX_FMT_BGR24:
297         picture->data[0] = ptr;
298         picture->data[1] = NULL;
299         picture->data[2] = NULL;
300         picture->linesize[0] = width * 3;
301         return size * 3;
302     case PIX_FMT_RGBA32:
303         picture->data[0] = ptr;
304         picture->data[1] = NULL;
305         picture->data[2] = NULL;
306         picture->linesize[0] = width * 4;
307         return size * 4;
308     case PIX_FMT_RGB555:
309     case PIX_FMT_RGB565:
310     case PIX_FMT_YUV422:
311         picture->data[0] = ptr;
312         picture->data[1] = NULL;
313         picture->data[2] = NULL;
314         picture->linesize[0] = width * 2;
315         return size * 2;
316     case PIX_FMT_UYVY422:
317         picture->data[0] = ptr;
318         picture->data[1] = NULL;
319         picture->data[2] = NULL;
320         picture->linesize[0] = width * 2;
321         return size * 2;
322     case PIX_FMT_UYVY411:
323         picture->data[0] = ptr;
324         picture->data[1] = NULL;
325         picture->data[2] = NULL;
326         picture->linesize[0] = width + width/2;
327         return size + size/2;
328     case PIX_FMT_GRAY8:
329         picture->data[0] = ptr;
330         picture->data[1] = NULL;
331         picture->data[2] = NULL;
332         picture->linesize[0] = width;
333         return size;
334     case PIX_FMT_MONOWHITE:
335     case PIX_FMT_MONOBLACK:
336         picture->data[0] = ptr;
337         picture->data[1] = NULL;
338         picture->data[2] = NULL;
339         picture->linesize[0] = (width + 7) >> 3;
340         return picture->linesize[0] * height;
341     case PIX_FMT_PAL8:
342         size2 = (size + 3) & ~3;
343         picture->data[0] = ptr;
344         picture->data[1] = ptr + size2; /* palette is stored here as 256 32 bit words */
345         picture->data[2] = NULL;
346         picture->linesize[0] = width;
347         picture->linesize[1] = 4;
348         return size2 + 256 * 4;
349     default:
350 fail:
351         picture->data[0] = NULL;
352         picture->data[1] = NULL;
353         picture->data[2] = NULL;
354         picture->data[3] = NULL;
355         return -1;
356     }
357 }
358
359 int avpicture_layout(const AVPicture* src, int pix_fmt, int width, int height,
360                      unsigned char *dest, int dest_size)
361 {
362     const PixFmtInfo* pf = &pix_fmt_info[pix_fmt];
363     int i, j, w, h, data_planes;
364     const unsigned char* s;
365     int size = avpicture_get_size(pix_fmt, width, height);
366
367     if (size > dest_size || size < 0)
368         return -1;
369
370     if (pf->pixel_type == FF_PIXEL_PACKED || pf->pixel_type == FF_PIXEL_PALETTE) {
371         if (pix_fmt == PIX_FMT_YUV422 ||
372             pix_fmt == PIX_FMT_UYVY422 ||
373             pix_fmt == PIX_FMT_RGB565 ||
374             pix_fmt == PIX_FMT_RGB555)
375             w = width * 2;
376         else if (pix_fmt == PIX_FMT_UYVY411)
377           w = width + width/2;
378         else if (pix_fmt == PIX_FMT_PAL8)
379           w = width;
380         else
381           w = width * (pf->depth * pf->nb_channels / 8);
382
383         data_planes = 1;
384         h = height;
385     } else {
386         data_planes = pf->nb_channels;
387         w = (width*pf->depth + 7)/8;
388         h = height;
389     }
390
391     for (i=0; i<data_planes; i++) {
392          if (i == 1) {
393              w = width >> pf->x_chroma_shift;
394              h = height >> pf->y_chroma_shift;
395          }
396          s = src->data[i];
397          for(j=0; j<h; j++) {
398              memcpy(dest, s, w);
399              dest += w;
400              s += src->linesize[i];
401          }
402     }
403
404     if (pf->pixel_type == FF_PIXEL_PALETTE)
405         memcpy((unsigned char *)(((size_t)dest + 3) & ~3), src->data[1], 256 * 4);
406
407     return size;
408 }
409
410 int avpicture_get_size(int pix_fmt, int width, int height)
411 {
412     AVPicture dummy_pict;
413     return avpicture_fill(&dummy_pict, NULL, pix_fmt, width, height);
414 }
415
416 /**
417  * compute the loss when converting from a pixel format to another
418  */
419 int avcodec_get_pix_fmt_loss(int dst_pix_fmt, int src_pix_fmt,
420                              int has_alpha)
421 {
422     const PixFmtInfo *pf, *ps;
423     int loss;
424
425     ps = &pix_fmt_info[src_pix_fmt];
426     pf = &pix_fmt_info[dst_pix_fmt];
427
428     /* compute loss */
429     loss = 0;
430     pf = &pix_fmt_info[dst_pix_fmt];
431     if (pf->depth < ps->depth ||
432         (dst_pix_fmt == PIX_FMT_RGB555 && src_pix_fmt == PIX_FMT_RGB565))
433         loss |= FF_LOSS_DEPTH;
434     if (pf->x_chroma_shift > ps->x_chroma_shift ||
435         pf->y_chroma_shift > ps->y_chroma_shift)
436         loss |= FF_LOSS_RESOLUTION;
437     switch(pf->color_type) {
438     case FF_COLOR_RGB:
439         if (ps->color_type != FF_COLOR_RGB &&
440             ps->color_type != FF_COLOR_GRAY)
441             loss |= FF_LOSS_COLORSPACE;
442         break;
443     case FF_COLOR_GRAY:
444         if (ps->color_type != FF_COLOR_GRAY)
445             loss |= FF_LOSS_COLORSPACE;
446         break;
447     case FF_COLOR_YUV:
448         if (ps->color_type != FF_COLOR_YUV)
449             loss |= FF_LOSS_COLORSPACE;
450         break;
451     case FF_COLOR_YUV_JPEG:
452         if (ps->color_type != FF_COLOR_YUV_JPEG &&
453             ps->color_type != FF_COLOR_YUV &&
454             ps->color_type != FF_COLOR_GRAY)
455             loss |= FF_LOSS_COLORSPACE;
456         break;
457     default:
458         /* fail safe test */
459         if (ps->color_type != pf->color_type)
460             loss |= FF_LOSS_COLORSPACE;
461         break;
462     }
463     if (pf->color_type == FF_COLOR_GRAY &&
464         ps->color_type != FF_COLOR_GRAY)
465         loss |= FF_LOSS_CHROMA;
466     if (!pf->is_alpha && (ps->is_alpha && has_alpha))
467         loss |= FF_LOSS_ALPHA;
468     if (pf->pixel_type == FF_PIXEL_PALETTE &&
469         (ps->pixel_type != FF_PIXEL_PALETTE && ps->color_type != FF_COLOR_GRAY))
470         loss |= FF_LOSS_COLORQUANT;
471     return loss;
472 }
473
474 static int avg_bits_per_pixel(int pix_fmt)
475 {
476     int bits;
477     const PixFmtInfo *pf;
478
479     pf = &pix_fmt_info[pix_fmt];
480     switch(pf->pixel_type) {
481     case FF_PIXEL_PACKED:
482         switch(pix_fmt) {
483         case PIX_FMT_YUV422:
484         case PIX_FMT_UYVY422:
485         case PIX_FMT_RGB565:
486         case PIX_FMT_RGB555:
487             bits = 16;
488             break;
489         case PIX_FMT_UYVY411:
490             bits = 12;
491             break;
492         default:
493             bits = pf->depth * pf->nb_channels;
494             break;
495         }
496         break;
497     case FF_PIXEL_PLANAR:
498         if (pf->x_chroma_shift == 0 && pf->y_chroma_shift == 0) {
499             bits = pf->depth * pf->nb_channels;
500         } else {
501             bits = pf->depth + ((2 * pf->depth) >>
502                                 (pf->x_chroma_shift + pf->y_chroma_shift));
503         }
504         break;
505     case FF_PIXEL_PALETTE:
506         bits = 8;
507         break;
508     default:
509         bits = -1;
510         break;
511     }
512     return bits;
513 }
514
515 static int avcodec_find_best_pix_fmt1(int pix_fmt_mask,
516                                       int src_pix_fmt,
517                                       int has_alpha,
518                                       int loss_mask)
519 {
520     int dist, i, loss, min_dist, dst_pix_fmt;
521
522     /* find exact color match with smallest size */
523     dst_pix_fmt = -1;
524     min_dist = 0x7fffffff;
525     for(i = 0;i < PIX_FMT_NB; i++) {
526         if (pix_fmt_mask & (1 << i)) {
527             loss = avcodec_get_pix_fmt_loss(i, src_pix_fmt, has_alpha) & loss_mask;
528             if (loss == 0) {
529                 dist = avg_bits_per_pixel(i);
530                 if (dist < min_dist) {
531                     min_dist = dist;
532                     dst_pix_fmt = i;
533                 }
534             }
535         }
536     }
537     return dst_pix_fmt;
538 }
539
540 /**
541  * find best pixel format to convert to. Return -1 if none found
542  */
543 int avcodec_find_best_pix_fmt(int pix_fmt_mask, int src_pix_fmt,
544                               int has_alpha, int *loss_ptr)
545 {
546     int dst_pix_fmt, loss_mask, i;
547     static const int loss_mask_order[] = {
548         ~0, /* no loss first */
549         ~FF_LOSS_ALPHA,
550         ~FF_LOSS_RESOLUTION,
551         ~(FF_LOSS_COLORSPACE | FF_LOSS_RESOLUTION),
552         ~FF_LOSS_COLORQUANT,
553         ~FF_LOSS_DEPTH,
554         0,
555     };
556
557     /* try with successive loss */
558     i = 0;
559     for(;;) {
560         loss_mask = loss_mask_order[i++];
561         dst_pix_fmt = avcodec_find_best_pix_fmt1(pix_fmt_mask, src_pix_fmt,
562                                                  has_alpha, loss_mask);
563         if (dst_pix_fmt >= 0)
564             goto found;
565         if (loss_mask == 0)
566             break;
567     }
568     return -1;
569  found:
570     if (loss_ptr)
571         *loss_ptr = avcodec_get_pix_fmt_loss(dst_pix_fmt, src_pix_fmt, has_alpha);
572     return dst_pix_fmt;
573 }
574
575 void ff_img_copy_plane(uint8_t *dst, int dst_wrap,
576                            const uint8_t *src, int src_wrap,
577                            int width, int height)
578 {
579     if((!dst) || (!src))
580         return;
581     for(;height > 0; height--) {
582         memcpy(dst, src, width);
583         dst += dst_wrap;
584         src += src_wrap;
585     }
586 }
587
588 /**
589  * Copy image 'src' to 'dst'.
590  */
591 void img_copy(AVPicture *dst, const AVPicture *src,
592               int pix_fmt, int width, int height)
593 {
594     int bwidth, bits, i;
595     const PixFmtInfo *pf = &pix_fmt_info[pix_fmt];
596
597     pf = &pix_fmt_info[pix_fmt];
598     switch(pf->pixel_type) {
599     case FF_PIXEL_PACKED:
600         switch(pix_fmt) {
601         case PIX_FMT_YUV422:
602         case PIX_FMT_UYVY422:
603         case PIX_FMT_RGB565:
604         case PIX_FMT_RGB555:
605             bits = 16;
606             break;
607         case PIX_FMT_UYVY411:
608             bits = 12;
609             break;
610         default:
611             bits = pf->depth * pf->nb_channels;
612             break;
613         }
614         bwidth = (width * bits + 7) >> 3;
615         ff_img_copy_plane(dst->data[0], dst->linesize[0],
616                        src->data[0], src->linesize[0],
617                        bwidth, height);
618         break;
619     case FF_PIXEL_PLANAR:
620         for(i = 0; i < pf->nb_channels; i++) {
621             int w, h;
622             w = width;
623             h = height;
624             if (i == 1 || i == 2) {
625                 w >>= pf->x_chroma_shift;
626                 h >>= pf->y_chroma_shift;
627             }
628             bwidth = (w * pf->depth + 7) >> 3;
629             ff_img_copy_plane(dst->data[i], dst->linesize[i],
630                            src->data[i], src->linesize[i],
631                            bwidth, h);
632         }
633         break;
634     case FF_PIXEL_PALETTE:
635         ff_img_copy_plane(dst->data[0], dst->linesize[0],
636                        src->data[0], src->linesize[0],
637                        width, height);
638         /* copy the palette */
639         ff_img_copy_plane(dst->data[1], dst->linesize[1],
640                        src->data[1], src->linesize[1],
641                        4, 256);
642         break;
643     }
644 }
645
646 /* XXX: totally non optimized */
647
648 static void yuv422_to_yuv420p(AVPicture *dst, const AVPicture *src,
649                               int width, int height)
650 {
651     const uint8_t *p, *p1;
652     uint8_t *lum, *cr, *cb, *lum1, *cr1, *cb1;
653     int w;
654
655     p1 = src->data[0];
656     lum1 = dst->data[0];
657     cb1 = dst->data[1];
658     cr1 = dst->data[2];
659
660     for(;height >= 1; height -= 2) {
661         p = p1;
662         lum = lum1;
663         cb = cb1;
664         cr = cr1;
665         for(w = width; w >= 2; w -= 2) {
666             lum[0] = p[0];
667             cb[0] = p[1];
668             lum[1] = p[2];
669             cr[0] = p[3];
670             p += 4;
671             lum += 2;
672             cb++;
673             cr++;
674         }
675         if (w) {
676             lum[0] = p[0];
677             cb[0] = p[1];
678             cr[0] = p[3];
679             cb++;
680             cr++;
681         }
682         p1 += src->linesize[0];
683         lum1 += dst->linesize[0];
684         if (height>1) {
685             p = p1;
686             lum = lum1;
687             for(w = width; w >= 2; w -= 2) {
688                 lum[0] = p[0];
689                 lum[1] = p[2];
690                 p += 4;
691                 lum += 2;
692             }
693             if (w) {
694                 lum[0] = p[0];
695             }
696             p1 += src->linesize[0];
697             lum1 += dst->linesize[0];
698         }
699         cb1 += dst->linesize[1];
700         cr1 += dst->linesize[2];
701     }
702 }
703
704 static void uyvy422_to_yuv420p(AVPicture *dst, const AVPicture *src,
705                               int width, int height)
706 {
707     const uint8_t *p, *p1;
708     uint8_t *lum, *cr, *cb, *lum1, *cr1, *cb1;
709     int w;
710
711     p1 = src->data[0];
712
713     lum1 = dst->data[0];
714     cb1 = dst->data[1];
715     cr1 = dst->data[2];
716
717     for(;height >= 1; height -= 2) {
718         p = p1;
719         lum = lum1;
720         cb = cb1;
721         cr = cr1;
722         for(w = width; w >= 2; w -= 2) {
723             lum[0] = p[1];
724             cb[0] = p[0];
725             lum[1] = p[3];
726             cr[0] = p[2];
727             p += 4;
728             lum += 2;
729             cb++;
730             cr++;
731         }
732         if (w) {
733             lum[0] = p[1];
734             cb[0] = p[0];
735             cr[0] = p[2];
736             cb++;
737             cr++;
738         }
739         p1 += src->linesize[0];
740         lum1 += dst->linesize[0];
741         if (height>1) {
742             p = p1;
743             lum = lum1;
744             for(w = width; w >= 2; w -= 2) {
745                 lum[0] = p[1];
746                 lum[1] = p[3];
747                 p += 4;
748                 lum += 2;
749             }
750             if (w) {
751                 lum[0] = p[1];
752             }
753             p1 += src->linesize[0];
754             lum1 += dst->linesize[0];
755         }
756         cb1 += dst->linesize[1];
757         cr1 += dst->linesize[2];
758     }
759 }
760
761
762 static void uyvy422_to_yuv422p(AVPicture *dst, const AVPicture *src,
763                               int width, int height)
764 {
765     const uint8_t *p, *p1;
766     uint8_t *lum, *cr, *cb, *lum1, *cr1, *cb1;
767     int w;
768
769     p1 = src->data[0];
770     lum1 = dst->data[0];
771     cb1 = dst->data[1];
772     cr1 = dst->data[2];
773     for(;height > 0; height--) {
774         p = p1;
775         lum = lum1;
776         cb = cb1;
777         cr = cr1;
778         for(w = width; w >= 2; w -= 2) {
779             lum[0] = p[1];
780             cb[0] = p[0];
781             lum[1] = p[3];
782             cr[0] = p[2];
783             p += 4;
784             lum += 2;
785             cb++;
786             cr++;
787         }
788         p1 += src->linesize[0];
789         lum1 += dst->linesize[0];
790         cb1 += dst->linesize[1];
791         cr1 += dst->linesize[2];
792     }
793 }
794
795
796 static void yuv422_to_yuv422p(AVPicture *dst, const AVPicture *src,
797                               int width, int height)
798 {
799     const uint8_t *p, *p1;
800     uint8_t *lum, *cr, *cb, *lum1, *cr1, *cb1;
801     int w;
802
803     p1 = src->data[0];
804     lum1 = dst->data[0];
805     cb1 = dst->data[1];
806     cr1 = dst->data[2];
807     for(;height > 0; height--) {
808         p = p1;
809         lum = lum1;
810         cb = cb1;
811         cr = cr1;
812         for(w = width; w >= 2; w -= 2) {
813             lum[0] = p[0];
814             cb[0] = p[1];
815             lum[1] = p[2];
816             cr[0] = p[3];
817             p += 4;
818             lum += 2;
819             cb++;
820             cr++;
821         }
822         p1 += src->linesize[0];
823         lum1 += dst->linesize[0];
824         cb1 += dst->linesize[1];
825         cr1 += dst->linesize[2];
826     }
827 }
828
829 static void yuv422p_to_yuv422(AVPicture *dst, const AVPicture *src,
830                               int width, int height)
831 {
832     uint8_t *p, *p1;
833     const uint8_t *lum, *cr, *cb, *lum1, *cr1, *cb1;
834     int w;
835
836     p1 = dst->data[0];
837     lum1 = src->data[0];
838     cb1 = src->data[1];
839     cr1 = src->data[2];
840     for(;height > 0; height--) {
841         p = p1;
842         lum = lum1;
843         cb = cb1;
844         cr = cr1;
845         for(w = width; w >= 2; w -= 2) {
846             p[0] = lum[0];
847             p[1] = cb[0];
848             p[2] = lum[1];
849             p[3] = cr[0];
850             p += 4;
851             lum += 2;
852             cb++;
853             cr++;
854         }
855         p1 += dst->linesize[0];
856         lum1 += src->linesize[0];
857         cb1 += src->linesize[1];
858         cr1 += src->linesize[2];
859     }
860 }
861
862 static void yuv422p_to_uyvy422(AVPicture *dst, const AVPicture *src,
863                               int width, int height)
864 {
865     uint8_t *p, *p1;
866     const uint8_t *lum, *cr, *cb, *lum1, *cr1, *cb1;
867     int w;
868
869     p1 = dst->data[0];
870     lum1 = src->data[0];
871     cb1 = src->data[1];
872     cr1 = src->data[2];
873     for(;height > 0; height--) {
874         p = p1;
875         lum = lum1;
876         cb = cb1;
877         cr = cr1;
878         for(w = width; w >= 2; w -= 2) {
879             p[1] = lum[0];
880             p[0] = cb[0];
881             p[3] = lum[1];
882             p[2] = cr[0];
883             p += 4;
884             lum += 2;
885             cb++;
886             cr++;
887         }
888         p1 += dst->linesize[0];
889         lum1 += src->linesize[0];
890         cb1 += src->linesize[1];
891         cr1 += src->linesize[2];
892     }
893 }
894
895 static void uyvy411_to_yuv411p(AVPicture *dst, const AVPicture *src,
896                               int width, int height)
897 {
898     const uint8_t *p, *p1;
899     uint8_t *lum, *cr, *cb, *lum1, *cr1, *cb1;
900     int w;
901
902     p1 = src->data[0];
903     lum1 = dst->data[0];
904     cb1 = dst->data[1];
905     cr1 = dst->data[2];
906     for(;height > 0; height--) {
907         p = p1;
908         lum = lum1;
909         cb = cb1;
910         cr = cr1;
911         for(w = width; w >= 4; w -= 4) {
912             cb[0] = p[0];
913             lum[0] = p[1];
914             lum[1] = p[2];
915             cr[0] = p[3];
916             lum[2] = p[4];
917             lum[3] = p[5];
918             p += 6;
919             lum += 4;
920             cb++;
921             cr++;
922         }
923         p1 += src->linesize[0];
924         lum1 += dst->linesize[0];
925         cb1 += dst->linesize[1];
926         cr1 += dst->linesize[2];
927     }
928 }
929
930
931 static void yuv420p_to_yuv422(AVPicture *dst, const AVPicture *src,
932                               int width, int height)
933 {
934     int w, h;
935     uint8_t *line1, *line2, *linesrc = dst->data[0];
936     uint8_t *lum1, *lum2, *lumsrc = src->data[0];
937     uint8_t *cb1, *cb2 = src->data[1];
938     uint8_t *cr1, *cr2 = src->data[2];
939
940     for(h = height / 2; h--;) {
941         line1 = linesrc;
942         line2 = linesrc + dst->linesize[0];
943
944         lum1 = lumsrc;
945         lum2 = lumsrc + src->linesize[0];
946
947         cb1 = cb2;
948         cr1 = cr2;
949
950         for(w = width / 2; w--;) {
951                 *line1++ = *lum1++; *line2++ = *lum2++;
952                 *line1++ =          *line2++ = *cb1++;
953                 *line1++ = *lum1++; *line2++ = *lum2++;
954                 *line1++ =          *line2++ = *cr1++;
955         }
956
957         linesrc += dst->linesize[0] * 2;
958         lumsrc += src->linesize[0] * 2;
959         cb2 += src->linesize[1];
960         cr2 += src->linesize[2];
961     }
962 }
963
964 static void yuv420p_to_uyvy422(AVPicture *dst, const AVPicture *src,
965                               int width, int height)
966 {
967     int w, h;
968     uint8_t *line1, *line2, *linesrc = dst->data[0];
969     uint8_t *lum1, *lum2, *lumsrc = src->data[0];
970     uint8_t *cb1, *cb2 = src->data[1];
971     uint8_t *cr1, *cr2 = src->data[2];
972
973     for(h = height / 2; h--;) {
974         line1 = linesrc;
975         line2 = linesrc + dst->linesize[0];
976
977         lum1 = lumsrc;
978         lum2 = lumsrc + src->linesize[0];
979
980         cb1 = cb2;
981         cr1 = cr2;
982
983         for(w = width / 2; w--;) {
984                 *line1++ =          *line2++ = *cb1++;
985                 *line1++ = *lum1++; *line2++ = *lum2++;
986                 *line1++ =          *line2++ = *cr1++;
987                 *line1++ = *lum1++; *line2++ = *lum2++;
988         }
989
990         linesrc += dst->linesize[0] * 2;
991         lumsrc += src->linesize[0] * 2;
992         cb2 += src->linesize[1];
993         cr2 += src->linesize[2];
994     }
995 }
996
997 #define SCALEBITS 10
998 #define ONE_HALF  (1 << (SCALEBITS - 1))
999 #define FIX(x)    ((int) ((x) * (1<<SCALEBITS) + 0.5))
1000
1001 #define YUV_TO_RGB1_CCIR(cb1, cr1)\
1002 {\
1003     cb = (cb1) - 128;\
1004     cr = (cr1) - 128;\
1005     r_add = FIX(1.40200*255.0/224.0) * cr + ONE_HALF;\
1006     g_add = - FIX(0.34414*255.0/224.0) * cb - FIX(0.71414*255.0/224.0) * cr + \
1007             ONE_HALF;\
1008     b_add = FIX(1.77200*255.0/224.0) * cb + ONE_HALF;\
1009 }
1010
1011 #define YUV_TO_RGB2_CCIR(r, g, b, y1)\
1012 {\
1013     y = ((y1) - 16) * FIX(255.0/219.0);\
1014     r = cm[(y + r_add) >> SCALEBITS];\
1015     g = cm[(y + g_add) >> SCALEBITS];\
1016     b = cm[(y + b_add) >> SCALEBITS];\
1017 }
1018
1019 #define YUV_TO_RGB1(cb1, cr1)\
1020 {\
1021     cb = (cb1) - 128;\
1022     cr = (cr1) - 128;\
1023     r_add = FIX(1.40200) * cr + ONE_HALF;\
1024     g_add = - FIX(0.34414) * cb - FIX(0.71414) * cr + ONE_HALF;\
1025     b_add = FIX(1.77200) * cb + ONE_HALF;\
1026 }
1027
1028 #define YUV_TO_RGB2(r, g, b, y1)\
1029 {\
1030     y = (y1) << SCALEBITS;\
1031     r = cm[(y + r_add) >> SCALEBITS];\
1032     g = cm[(y + g_add) >> SCALEBITS];\
1033     b = cm[(y + b_add) >> SCALEBITS];\
1034 }
1035
1036 #define Y_CCIR_TO_JPEG(y)\
1037  cm[((y) * FIX(255.0/219.0) + (ONE_HALF - 16 * FIX(255.0/219.0))) >> SCALEBITS]
1038
1039 #define Y_JPEG_TO_CCIR(y)\
1040  (((y) * FIX(219.0/255.0) + (ONE_HALF + (16 << SCALEBITS))) >> SCALEBITS)
1041
1042 #define C_CCIR_TO_JPEG(y)\
1043  cm[(((y) - 128) * FIX(127.0/112.0) + (ONE_HALF + (128 << SCALEBITS))) >> SCALEBITS]
1044
1045 /* NOTE: the clamp is really necessary! */
1046 static inline int C_JPEG_TO_CCIR(int y) {
1047     y = (((y - 128) * FIX(112.0/127.0) + (ONE_HALF + (128 << SCALEBITS))) >> SCALEBITS);
1048     if (y < 16)
1049         y = 16;
1050     return y;
1051 }
1052
1053
1054 #define RGB_TO_Y(r, g, b) \
1055 ((FIX(0.29900) * (r) + FIX(0.58700) * (g) + \
1056   FIX(0.11400) * (b) + ONE_HALF) >> SCALEBITS)
1057
1058 #define RGB_TO_U(r1, g1, b1, shift)\
1059 (((- FIX(0.16874) * r1 - FIX(0.33126) * g1 +         \
1060      FIX(0.50000) * b1 + (ONE_HALF << shift) - 1) >> (SCALEBITS + shift)) + 128)
1061
1062 #define RGB_TO_V(r1, g1, b1, shift)\
1063 (((FIX(0.50000) * r1 - FIX(0.41869) * g1 -           \
1064    FIX(0.08131) * b1 + (ONE_HALF << shift) - 1) >> (SCALEBITS + shift)) + 128)
1065
1066 #define RGB_TO_Y_CCIR(r, g, b) \
1067 ((FIX(0.29900*219.0/255.0) * (r) + FIX(0.58700*219.0/255.0) * (g) + \
1068   FIX(0.11400*219.0/255.0) * (b) + (ONE_HALF + (16 << SCALEBITS))) >> SCALEBITS)
1069
1070 #define RGB_TO_U_CCIR(r1, g1, b1, shift)\
1071 (((- FIX(0.16874*224.0/255.0) * r1 - FIX(0.33126*224.0/255.0) * g1 +         \
1072      FIX(0.50000*224.0/255.0) * b1 + (ONE_HALF << shift) - 1) >> (SCALEBITS + shift)) + 128)
1073
1074 #define RGB_TO_V_CCIR(r1, g1, b1, shift)\
1075 (((FIX(0.50000*224.0/255.0) * r1 - FIX(0.41869*224.0/255.0) * g1 -           \
1076    FIX(0.08131*224.0/255.0) * b1 + (ONE_HALF << shift) - 1) >> (SCALEBITS + shift)) + 128)
1077
1078 static uint8_t y_ccir_to_jpeg[256];
1079 static uint8_t y_jpeg_to_ccir[256];
1080 static uint8_t c_ccir_to_jpeg[256];
1081 static uint8_t c_jpeg_to_ccir[256];
1082
1083 /* init various conversion tables */
1084 static void img_convert_init(void)
1085 {
1086     int i;
1087     uint8_t *cm = cropTbl + MAX_NEG_CROP;
1088
1089     for(i = 0;i < 256; i++) {
1090         y_ccir_to_jpeg[i] = Y_CCIR_TO_JPEG(i);
1091         y_jpeg_to_ccir[i] = Y_JPEG_TO_CCIR(i);
1092         c_ccir_to_jpeg[i] = C_CCIR_TO_JPEG(i);
1093         c_jpeg_to_ccir[i] = C_JPEG_TO_CCIR(i);
1094     }
1095 }
1096
1097 /* apply to each pixel the given table */
1098 static void img_apply_table(uint8_t *dst, int dst_wrap,
1099                             const uint8_t *src, int src_wrap,
1100                             int width, int height, const uint8_t *table1)
1101 {
1102     int n;
1103     const uint8_t *s;
1104     uint8_t *d;
1105     const uint8_t *table;
1106
1107     table = table1;
1108     for(;height > 0; height--) {
1109         s = src;
1110         d = dst;
1111         n = width;
1112         while (n >= 4) {
1113             d[0] = table[s[0]];
1114             d[1] = table[s[1]];
1115             d[2] = table[s[2]];
1116             d[3] = table[s[3]];
1117             d += 4;
1118             s += 4;
1119             n -= 4;
1120         }
1121         while (n > 0) {
1122             d[0] = table[s[0]];
1123             d++;
1124             s++;
1125             n--;
1126         }
1127         dst += dst_wrap;
1128         src += src_wrap;
1129     }
1130 }
1131
1132 /* XXX: use generic filter ? */
1133 /* XXX: in most cases, the sampling position is incorrect */
1134
1135 /* 4x1 -> 1x1 */
1136 static void shrink41(uint8_t *dst, int dst_wrap,
1137                      const uint8_t *src, int src_wrap,
1138                      int width, int height)
1139 {
1140     int w;
1141     const uint8_t *s;
1142     uint8_t *d;
1143
1144     for(;height > 0; height--) {
1145         s = src;
1146         d = dst;
1147         for(w = width;w > 0; w--) {
1148             d[0] = (s[0] + s[1] + s[2] + s[3] + 2) >> 2;
1149             s += 4;
1150             d++;
1151         }
1152         src += src_wrap;
1153         dst += dst_wrap;
1154     }
1155 }
1156
1157 /* 2x1 -> 1x1 */
1158 static void shrink21(uint8_t *dst, int dst_wrap,
1159                      const uint8_t *src, int src_wrap,
1160                      int width, int height)
1161 {
1162     int w;
1163     const uint8_t *s;
1164     uint8_t *d;
1165
1166     for(;height > 0; height--) {
1167         s = src;
1168         d = dst;
1169         for(w = width;w > 0; w--) {
1170             d[0] = (s[0] + s[1]) >> 1;
1171             s += 2;
1172             d++;
1173         }
1174         src += src_wrap;
1175         dst += dst_wrap;
1176     }
1177 }
1178
1179 /* 1x2 -> 1x1 */
1180 static void shrink12(uint8_t *dst, int dst_wrap,
1181                      const uint8_t *src, int src_wrap,
1182                      int width, int height)
1183 {
1184     int w;
1185     uint8_t *d;
1186     const uint8_t *s1, *s2;
1187
1188     for(;height > 0; height--) {
1189         s1 = src;
1190         s2 = s1 + src_wrap;
1191         d = dst;
1192         for(w = width;w >= 4; w-=4) {
1193             d[0] = (s1[0] + s2[0]) >> 1;
1194             d[1] = (s1[1] + s2[1]) >> 1;
1195             d[2] = (s1[2] + s2[2]) >> 1;
1196             d[3] = (s1[3] + s2[3]) >> 1;
1197             s1 += 4;
1198             s2 += 4;
1199             d += 4;
1200         }
1201         for(;w > 0; w--) {
1202             d[0] = (s1[0] + s2[0]) >> 1;
1203             s1++;
1204             s2++;
1205             d++;
1206         }
1207         src += 2 * src_wrap;
1208         dst += dst_wrap;
1209     }
1210 }
1211
1212 /* 2x2 -> 1x1 */
1213 void ff_shrink22(uint8_t *dst, int dst_wrap,
1214                      const uint8_t *src, int src_wrap,
1215                      int width, int height)
1216 {
1217     int w;
1218     const uint8_t *s1, *s2;
1219     uint8_t *d;
1220
1221     for(;height > 0; height--) {
1222         s1 = src;
1223         s2 = s1 + src_wrap;
1224         d = dst;
1225         for(w = width;w >= 4; w-=4) {
1226             d[0] = (s1[0] + s1[1] + s2[0] + s2[1] + 2) >> 2;
1227             d[1] = (s1[2] + s1[3] + s2[2] + s2[3] + 2) >> 2;
1228             d[2] = (s1[4] + s1[5] + s2[4] + s2[5] + 2) >> 2;
1229             d[3] = (s1[6] + s1[7] + s2[6] + s2[7] + 2) >> 2;
1230             s1 += 8;
1231             s2 += 8;
1232             d += 4;
1233         }
1234         for(;w > 0; w--) {
1235             d[0] = (s1[0] + s1[1] + s2[0] + s2[1] + 2) >> 2;
1236             s1 += 2;
1237             s2 += 2;
1238             d++;
1239         }
1240         src += 2 * src_wrap;
1241         dst += dst_wrap;
1242     }
1243 }
1244
1245 /* 4x4 -> 1x1 */
1246 void ff_shrink44(uint8_t *dst, int dst_wrap,
1247                      const uint8_t *src, int src_wrap,
1248                      int width, int height)
1249 {
1250     int w;
1251     const uint8_t *s1, *s2, *s3, *s4;
1252     uint8_t *d;
1253
1254     for(;height > 0; height--) {
1255         s1 = src;
1256         s2 = s1 + src_wrap;
1257         s3 = s2 + src_wrap;
1258         s4 = s3 + src_wrap;
1259         d = dst;
1260         for(w = width;w > 0; w--) {
1261             d[0] = (s1[0] + s1[1] + s1[2] + s1[3] +
1262                     s2[0] + s2[1] + s2[2] + s2[3] +
1263                     s3[0] + s3[1] + s3[2] + s3[3] +
1264                     s4[0] + s4[1] + s4[2] + s4[3] + 8) >> 4;
1265             s1 += 4;
1266             s2 += 4;
1267             s3 += 4;
1268             s4 += 4;
1269             d++;
1270         }
1271         src += 4 * src_wrap;
1272         dst += dst_wrap;
1273     }
1274 }
1275
1276 /* 8x8 -> 1x1 */
1277 void ff_shrink88(uint8_t *dst, int dst_wrap,
1278                      const uint8_t *src, int src_wrap,
1279                      int width, int height)
1280 {
1281     int w, i;
1282
1283     for(;height > 0; height--) {
1284         for(w = width;w > 0; w--) {
1285             int tmp=0;
1286             for(i=0; i<8; i++){
1287                 tmp += src[0] + src[1] + src[2] + src[3] + src[4] + src[5] + src[6] + src[7];
1288                 src += src_wrap;
1289             }
1290             *(dst++) = (tmp + 32)>>6;
1291             src += 8 - 8*src_wrap;
1292         }
1293         src += 8*src_wrap - 8*width;
1294         dst += dst_wrap - width;
1295     }
1296 }
1297
1298 static void grow21_line(uint8_t *dst, const uint8_t *src,
1299                         int width)
1300 {
1301     int w;
1302     const uint8_t *s1;
1303     uint8_t *d;
1304
1305     s1 = src;
1306     d = dst;
1307     for(w = width;w >= 4; w-=4) {
1308         d[1] = d[0] = s1[0];
1309         d[3] = d[2] = s1[1];
1310         s1 += 2;
1311         d += 4;
1312     }
1313     for(;w >= 2; w -= 2) {
1314         d[1] = d[0] = s1[0];
1315         s1 ++;
1316         d += 2;
1317     }
1318     /* only needed if width is not a multiple of two */
1319     /* XXX: veryfy that */
1320     if (w) {
1321         d[0] = s1[0];
1322     }
1323 }
1324
1325 static void grow41_line(uint8_t *dst, const uint8_t *src,
1326                         int width)
1327 {
1328     int w, v;
1329     const uint8_t *s1;
1330     uint8_t *d;
1331
1332     s1 = src;
1333     d = dst;
1334     for(w = width;w >= 4; w-=4) {
1335         v = s1[0];
1336         d[0] = v;
1337         d[1] = v;
1338         d[2] = v;
1339         d[3] = v;
1340         s1 ++;
1341         d += 4;
1342     }
1343 }
1344
1345 /* 1x1 -> 2x1 */
1346 static void grow21(uint8_t *dst, int dst_wrap,
1347                    const uint8_t *src, int src_wrap,
1348                    int width, int height)
1349 {
1350     for(;height > 0; height--) {
1351         grow21_line(dst, src, width);
1352         src += src_wrap;
1353         dst += dst_wrap;
1354     }
1355 }
1356
1357 /* 1x1 -> 2x2 */
1358 static void grow22(uint8_t *dst, int dst_wrap,
1359                    const uint8_t *src, int src_wrap,
1360                    int width, int height)
1361 {
1362     for(;height > 0; height--) {
1363         grow21_line(dst, src, width);
1364         if (height%2)
1365             src += src_wrap;
1366         dst += dst_wrap;
1367     }
1368 }
1369
1370 /* 1x1 -> 4x1 */
1371 static void grow41(uint8_t *dst, int dst_wrap,
1372                    const uint8_t *src, int src_wrap,
1373                    int width, int height)
1374 {
1375     for(;height > 0; height--) {
1376         grow41_line(dst, src, width);
1377         src += src_wrap;
1378         dst += dst_wrap;
1379     }
1380 }
1381
1382 /* 1x1 -> 4x4 */
1383 static void grow44(uint8_t *dst, int dst_wrap,
1384                    const uint8_t *src, int src_wrap,
1385                    int width, int height)
1386 {
1387     for(;height > 0; height--) {
1388         grow41_line(dst, src, width);
1389         if ((height & 3) == 1)
1390             src += src_wrap;
1391         dst += dst_wrap;
1392     }
1393 }
1394
1395 /* 1x2 -> 2x1 */
1396 static void conv411(uint8_t *dst, int dst_wrap,
1397                     const uint8_t *src, int src_wrap,
1398                     int width, int height)
1399 {
1400     int w, c;
1401     const uint8_t *s1, *s2;
1402     uint8_t *d;
1403
1404     width>>=1;
1405
1406     for(;height > 0; height--) {
1407         s1 = src;
1408         s2 = src + src_wrap;
1409         d = dst;
1410         for(w = width;w > 0; w--) {
1411             c = (s1[0] + s2[0]) >> 1;
1412             d[0] = c;
1413             d[1] = c;
1414             s1++;
1415             s2++;
1416             d += 2;
1417         }
1418         src += src_wrap * 2;
1419         dst += dst_wrap;
1420     }
1421 }
1422
1423 /* XXX: add jpeg quantize code */
1424
1425 #define TRANSP_INDEX (6*6*6)
1426
1427 /* this is maybe slow, but allows for extensions */
1428 static inline unsigned char gif_clut_index(uint8_t r, uint8_t g, uint8_t b)
1429 {
1430     return ((((r)/47)%6)*6*6+(((g)/47)%6)*6+(((b)/47)%6));
1431 }
1432
1433 static void build_rgb_palette(uint8_t *palette, int has_alpha)
1434 {
1435     uint32_t *pal;
1436     static const uint8_t pal_value[6] = { 0x00, 0x33, 0x66, 0x99, 0xcc, 0xff };
1437     int i, r, g, b;
1438
1439     pal = (uint32_t *)palette;
1440     i = 0;
1441     for(r = 0; r < 6; r++) {
1442         for(g = 0; g < 6; g++) {
1443             for(b = 0; b < 6; b++) {
1444                 pal[i++] = (0xff << 24) | (pal_value[r] << 16) |
1445                     (pal_value[g] << 8) | pal_value[b];
1446             }
1447         }
1448     }
1449     if (has_alpha)
1450         pal[i++] = 0;
1451     while (i < 256)
1452         pal[i++] = 0xff000000;
1453 }
1454
1455 /* copy bit n to bits 0 ... n - 1 */
1456 static inline unsigned int bitcopy_n(unsigned int a, int n)
1457 {
1458     int mask;
1459     mask = (1 << n) - 1;
1460     return (a & (0xff & ~mask)) | ((-((a >> n) & 1)) & mask);
1461 }
1462
1463 /* rgb555 handling */
1464
1465 #define RGB_NAME rgb555
1466
1467 #define RGB_IN(r, g, b, s)\
1468 {\
1469     unsigned int v = ((const uint16_t *)(s))[0];\
1470     r = bitcopy_n(v >> (10 - 3), 3);\
1471     g = bitcopy_n(v >> (5 - 3), 3);\
1472     b = bitcopy_n(v << 3, 3);\
1473 }
1474
1475 #define RGBA_IN(r, g, b, a, s)\
1476 {\
1477     unsigned int v = ((const uint16_t *)(s))[0];\
1478     r = bitcopy_n(v >> (10 - 3), 3);\
1479     g = bitcopy_n(v >> (5 - 3), 3);\
1480     b = bitcopy_n(v << 3, 3);\
1481     a = (-(v >> 15)) & 0xff;\
1482 }
1483
1484 #define RGBA_OUT(d, r, g, b, a)\
1485 {\
1486     ((uint16_t *)(d))[0] = ((r >> 3) << 10) | ((g >> 3) << 5) | (b >> 3) | \
1487                            ((a << 8) & 0x8000);\
1488 }
1489
1490 #define BPP 2
1491
1492 #include "imgconvert_template.h"
1493
1494 /* rgb565 handling */
1495
1496 #define RGB_NAME rgb565
1497
1498 #define RGB_IN(r, g, b, s)\
1499 {\
1500     unsigned int v = ((const uint16_t *)(s))[0];\
1501     r = bitcopy_n(v >> (11 - 3), 3);\
1502     g = bitcopy_n(v >> (5 - 2), 2);\
1503     b = bitcopy_n(v << 3, 3);\
1504 }
1505
1506 #define RGB_OUT(d, r, g, b)\
1507 {\
1508     ((uint16_t *)(d))[0] = ((r >> 3) << 11) | ((g >> 2) << 5) | (b >> 3);\
1509 }
1510
1511 #define BPP 2
1512
1513 #include "imgconvert_template.h"
1514
1515 /* bgr24 handling */
1516
1517 #define RGB_NAME bgr24
1518
1519 #define RGB_IN(r, g, b, s)\
1520 {\
1521     b = (s)[0];\
1522     g = (s)[1];\
1523     r = (s)[2];\
1524 }
1525
1526 #define RGB_OUT(d, r, g, b)\
1527 {\
1528     (d)[0] = b;\
1529     (d)[1] = g;\
1530     (d)[2] = r;\
1531 }
1532
1533 #define BPP 3
1534
1535 #include "imgconvert_template.h"
1536
1537 #undef RGB_IN
1538 #undef RGB_OUT
1539 #undef BPP
1540
1541 /* rgb24 handling */
1542
1543 #define RGB_NAME rgb24
1544 #define FMT_RGB24
1545
1546 #define RGB_IN(r, g, b, s)\
1547 {\
1548     r = (s)[0];\
1549     g = (s)[1];\
1550     b = (s)[2];\
1551 }
1552
1553 #define RGB_OUT(d, r, g, b)\
1554 {\
1555     (d)[0] = r;\
1556     (d)[1] = g;\
1557     (d)[2] = b;\
1558 }
1559
1560 #define BPP 3
1561
1562 #include "imgconvert_template.h"
1563
1564 /* rgba32 handling */
1565
1566 #define RGB_NAME rgba32
1567 #define FMT_RGBA32
1568
1569 #define RGB_IN(r, g, b, s)\
1570 {\
1571     unsigned int v = ((const uint32_t *)(s))[0];\
1572     r = (v >> 16) & 0xff;\
1573     g = (v >> 8) & 0xff;\
1574     b = v & 0xff;\
1575 }
1576
1577 #define RGBA_IN(r, g, b, a, s)\
1578 {\
1579     unsigned int v = ((const uint32_t *)(s))[0];\
1580     a = (v >> 24) & 0xff;\
1581     r = (v >> 16) & 0xff;\
1582     g = (v >> 8) & 0xff;\
1583     b = v & 0xff;\
1584 }
1585
1586 #define RGBA_OUT(d, r, g, b, a)\
1587 {\
1588     ((uint32_t *)(d))[0] = (a << 24) | (r << 16) | (g << 8) | b;\
1589 }
1590
1591 #define BPP 4
1592
1593 #include "imgconvert_template.h"
1594
1595 static void mono_to_gray(AVPicture *dst, const AVPicture *src,
1596                          int width, int height, int xor_mask)
1597 {
1598     const unsigned char *p;
1599     unsigned char *q;
1600     int v, dst_wrap, src_wrap;
1601     int y, w;
1602
1603     p = src->data[0];
1604     src_wrap = src->linesize[0] - ((width + 7) >> 3);
1605
1606     q = dst->data[0];
1607     dst_wrap = dst->linesize[0] - width;
1608     for(y=0;y<height;y++) {
1609         w = width;
1610         while (w >= 8) {
1611             v = *p++ ^ xor_mask;
1612             q[0] = -(v >> 7);
1613             q[1] = -((v >> 6) & 1);
1614             q[2] = -((v >> 5) & 1);
1615             q[3] = -((v >> 4) & 1);
1616             q[4] = -((v >> 3) & 1);
1617             q[5] = -((v >> 2) & 1);
1618             q[6] = -((v >> 1) & 1);
1619             q[7] = -((v >> 0) & 1);
1620             w -= 8;
1621             q += 8;
1622         }
1623         if (w > 0) {
1624             v = *p++ ^ xor_mask;
1625             do {
1626                 q[0] = -((v >> 7) & 1);
1627                 q++;
1628                 v <<= 1;
1629             } while (--w);
1630         }
1631         p += src_wrap;
1632         q += dst_wrap;
1633     }
1634 }
1635
1636 static void monowhite_to_gray(AVPicture *dst, const AVPicture *src,
1637                                int width, int height)
1638 {
1639     mono_to_gray(dst, src, width, height, 0xff);
1640 }
1641
1642 static void monoblack_to_gray(AVPicture *dst, const AVPicture *src,
1643                                int width, int height)
1644 {
1645     mono_to_gray(dst, src, width, height, 0x00);
1646 }
1647
1648 static void gray_to_mono(AVPicture *dst, const AVPicture *src,
1649                          int width, int height, int xor_mask)
1650 {
1651     int n;
1652     const uint8_t *s;
1653     uint8_t *d;
1654     int j, b, v, n1, src_wrap, dst_wrap, y;
1655
1656     s = src->data[0];
1657     src_wrap = src->linesize[0] - width;
1658
1659     d = dst->data[0];
1660     dst_wrap = dst->linesize[0] - ((width + 7) >> 3);
1661
1662     for(y=0;y<height;y++) {
1663         n = width;
1664         while (n >= 8) {
1665             v = 0;
1666             for(j=0;j<8;j++) {
1667                 b = s[0];
1668                 s++;
1669                 v = (v << 1) | (b >> 7);
1670             }
1671             d[0] = v ^ xor_mask;
1672             d++;
1673             n -= 8;
1674         }
1675         if (n > 0) {
1676             n1 = n;
1677             v = 0;
1678             while (n > 0) {
1679                 b = s[0];
1680                 s++;
1681                 v = (v << 1) | (b >> 7);
1682                 n--;
1683             }
1684             d[0] = (v << (8 - (n1 & 7))) ^ xor_mask;
1685             d++;
1686         }
1687         s += src_wrap;
1688         d += dst_wrap;
1689     }
1690 }
1691
1692 static void gray_to_monowhite(AVPicture *dst, const AVPicture *src,
1693                               int width, int height)
1694 {
1695     gray_to_mono(dst, src, width, height, 0xff);
1696 }
1697
1698 static void gray_to_monoblack(AVPicture *dst, const AVPicture *src,
1699                               int width, int height)
1700 {
1701     gray_to_mono(dst, src, width, height, 0x00);
1702 }
1703
1704 typedef struct ConvertEntry {
1705     void (*convert)(AVPicture *dst,
1706                     const AVPicture *src, int width, int height);
1707 } ConvertEntry;
1708
1709 /* Add each new convertion function in this table. In order to be able
1710    to convert from any format to any format, the following constraints
1711    must be satisfied:
1712
1713    - all FF_COLOR_RGB formats must convert to and from PIX_FMT_RGB24
1714
1715    - all FF_COLOR_GRAY formats must convert to and from PIX_FMT_GRAY8
1716
1717    - all FF_COLOR_RGB formats with alpha must convert to and from PIX_FMT_RGBA32
1718
1719    - PIX_FMT_YUV444P and PIX_FMT_YUVJ444P must convert to and from
1720      PIX_FMT_RGB24.
1721
1722    - PIX_FMT_422 must convert to and from PIX_FMT_422P.
1723
1724    The other conversion functions are just optimisations for common cases.
1725 */
1726 static const ConvertEntry convert_table[PIX_FMT_NB][PIX_FMT_NB] = {
1727     [PIX_FMT_YUV420P] = {
1728         [PIX_FMT_YUV422] = {
1729             .convert = yuv420p_to_yuv422,
1730         },
1731         [PIX_FMT_RGB555] = {
1732             .convert = yuv420p_to_rgb555
1733         },
1734         [PIX_FMT_RGB565] = {
1735             .convert = yuv420p_to_rgb565
1736         },
1737         [PIX_FMT_BGR24] = {
1738             .convert = yuv420p_to_bgr24
1739         },
1740         [PIX_FMT_RGB24] = {
1741             .convert = yuv420p_to_rgb24
1742         },
1743         [PIX_FMT_RGBA32] = {
1744             .convert = yuv420p_to_rgba32
1745         },
1746         [PIX_FMT_UYVY422] = {
1747             .convert = yuv420p_to_uyvy422,
1748         },
1749     },
1750     [PIX_FMT_YUV422P] = {
1751         [PIX_FMT_YUV422] = {
1752             .convert = yuv422p_to_yuv422,
1753         },
1754         [PIX_FMT_UYVY422] = {
1755             .convert = yuv422p_to_uyvy422,
1756         },
1757     },
1758     [PIX_FMT_YUV444P] = {
1759         [PIX_FMT_RGB24] = {
1760             .convert = yuv444p_to_rgb24
1761         },
1762     },
1763     [PIX_FMT_YUVJ420P] = {
1764         [PIX_FMT_RGB555] = {
1765             .convert = yuvj420p_to_rgb555
1766         },
1767         [PIX_FMT_RGB565] = {
1768             .convert = yuvj420p_to_rgb565
1769         },
1770         [PIX_FMT_BGR24] = {
1771             .convert = yuvj420p_to_bgr24
1772         },
1773         [PIX_FMT_RGB24] = {
1774             .convert = yuvj420p_to_rgb24
1775         },
1776         [PIX_FMT_RGBA32] = {
1777             .convert = yuvj420p_to_rgba32
1778         },
1779     },
1780     [PIX_FMT_YUVJ444P] = {
1781         [PIX_FMT_RGB24] = {
1782             .convert = yuvj444p_to_rgb24
1783         },
1784     },
1785     [PIX_FMT_YUV422] = {
1786         [PIX_FMT_YUV420P] = {
1787             .convert = yuv422_to_yuv420p,
1788         },
1789         [PIX_FMT_YUV422P] = {
1790             .convert = yuv422_to_yuv422p,
1791         },
1792     },
1793     [PIX_FMT_UYVY422] = {
1794         [PIX_FMT_YUV420P] = {
1795             .convert = uyvy422_to_yuv420p,
1796         },
1797         [PIX_FMT_YUV422P] = {
1798             .convert = uyvy422_to_yuv422p,
1799         },
1800     },
1801     [PIX_FMT_RGB24] = {
1802         [PIX_FMT_YUV420P] = {
1803             .convert = rgb24_to_yuv420p
1804         },
1805         [PIX_FMT_RGB565] = {
1806             .convert = rgb24_to_rgb565
1807         },
1808         [PIX_FMT_RGB555] = {
1809             .convert = rgb24_to_rgb555
1810         },
1811         [PIX_FMT_RGBA32] = {
1812             .convert = rgb24_to_rgba32
1813         },
1814         [PIX_FMT_BGR24] = {
1815             .convert = rgb24_to_bgr24
1816         },
1817         [PIX_FMT_GRAY8] = {
1818             .convert = rgb24_to_gray
1819         },
1820         [PIX_FMT_PAL8] = {
1821             .convert = rgb24_to_pal8
1822         },
1823         [PIX_FMT_YUV444P] = {
1824             .convert = rgb24_to_yuv444p
1825         },
1826         [PIX_FMT_YUVJ420P] = {
1827             .convert = rgb24_to_yuvj420p
1828         },
1829         [PIX_FMT_YUVJ444P] = {
1830             .convert = rgb24_to_yuvj444p
1831         },
1832     },
1833     [PIX_FMT_RGBA32] = {
1834         [PIX_FMT_RGB24] = {
1835             .convert = rgba32_to_rgb24
1836         },
1837         [PIX_FMT_RGB555] = {
1838             .convert = rgba32_to_rgb555
1839         },
1840         [PIX_FMT_PAL8] = {
1841             .convert = rgba32_to_pal8
1842         },
1843         [PIX_FMT_YUV420P] = {
1844             .convert = rgba32_to_yuv420p
1845         },
1846         [PIX_FMT_GRAY8] = {
1847             .convert = rgba32_to_gray
1848         },
1849     },
1850     [PIX_FMT_BGR24] = {
1851         [PIX_FMT_RGB24] = {
1852             .convert = bgr24_to_rgb24
1853         },
1854         [PIX_FMT_YUV420P] = {
1855             .convert = bgr24_to_yuv420p
1856         },
1857         [PIX_FMT_GRAY8] = {
1858             .convert = bgr24_to_gray
1859         },
1860     },
1861     [PIX_FMT_RGB555] = {
1862         [PIX_FMT_RGB24] = {
1863             .convert = rgb555_to_rgb24
1864         },
1865         [PIX_FMT_RGBA32] = {
1866             .convert = rgb555_to_rgba32
1867         },
1868         [PIX_FMT_YUV420P] = {
1869             .convert = rgb555_to_yuv420p
1870         },
1871         [PIX_FMT_GRAY8] = {
1872             .convert = rgb555_to_gray
1873         },
1874     },
1875     [PIX_FMT_RGB565] = {
1876         [PIX_FMT_RGB24] = {
1877             .convert = rgb565_to_rgb24
1878         },
1879         [PIX_FMT_YUV420P] = {
1880             .convert = rgb565_to_yuv420p
1881         },
1882         [PIX_FMT_GRAY8] = {
1883             .convert = rgb565_to_gray
1884         },
1885     },
1886     [PIX_FMT_GRAY8] = {
1887         [PIX_FMT_RGB555] = {
1888             .convert = gray_to_rgb555
1889         },
1890         [PIX_FMT_RGB565] = {
1891             .convert = gray_to_rgb565
1892         },
1893         [PIX_FMT_RGB24] = {
1894             .convert = gray_to_rgb24
1895         },
1896         [PIX_FMT_BGR24] = {
1897             .convert = gray_to_bgr24
1898         },
1899         [PIX_FMT_RGBA32] = {
1900             .convert = gray_to_rgba32
1901         },
1902         [PIX_FMT_MONOWHITE] = {
1903             .convert = gray_to_monowhite
1904         },
1905         [PIX_FMT_MONOBLACK] = {
1906             .convert = gray_to_monoblack
1907         },
1908     },
1909     [PIX_FMT_MONOWHITE] = {
1910         [PIX_FMT_GRAY8] = {
1911             .convert = monowhite_to_gray
1912         },
1913     },
1914     [PIX_FMT_MONOBLACK] = {
1915         [PIX_FMT_GRAY8] = {
1916             .convert = monoblack_to_gray
1917         },
1918     },
1919     [PIX_FMT_PAL8] = {
1920         [PIX_FMT_RGB555] = {
1921             .convert = pal8_to_rgb555
1922         },
1923         [PIX_FMT_RGB565] = {
1924             .convert = pal8_to_rgb565
1925         },
1926         [PIX_FMT_BGR24] = {
1927             .convert = pal8_to_bgr24
1928         },
1929         [PIX_FMT_RGB24] = {
1930             .convert = pal8_to_rgb24
1931         },
1932         [PIX_FMT_RGBA32] = {
1933             .convert = pal8_to_rgba32
1934         },
1935     },
1936     [PIX_FMT_UYVY411] = {
1937         [PIX_FMT_YUV411P] = {
1938             .convert = uyvy411_to_yuv411p,
1939         },
1940     },
1941
1942 };
1943
1944 int avpicture_alloc(AVPicture *picture,
1945                            int pix_fmt, int width, int height)
1946 {
1947     int size;
1948     void *ptr;
1949
1950     size = avpicture_get_size(pix_fmt, width, height);
1951     if(size<0)
1952         goto fail;
1953     ptr = av_malloc(size);
1954     if (!ptr)
1955         goto fail;
1956     avpicture_fill(picture, ptr, pix_fmt, width, height);
1957     return 0;
1958  fail:
1959     memset(picture, 0, sizeof(AVPicture));
1960     return -1;
1961 }
1962
1963 void avpicture_free(AVPicture *picture)
1964 {
1965     av_free(picture->data[0]);
1966 }
1967
1968 /* return true if yuv planar */
1969 static inline int is_yuv_planar(const PixFmtInfo *ps)
1970 {
1971     return (ps->color_type == FF_COLOR_YUV ||
1972             ps->color_type == FF_COLOR_YUV_JPEG) &&
1973         ps->pixel_type == FF_PIXEL_PLANAR;
1974 }
1975
1976 /**
1977  * Crop image top and left side
1978  */
1979 int img_crop(AVPicture *dst, const AVPicture *src,
1980               int pix_fmt, int top_band, int left_band)
1981 {
1982     int y_shift;
1983     int x_shift;
1984
1985     if (pix_fmt < 0 || pix_fmt >= PIX_FMT_NB || !is_yuv_planar(&pix_fmt_info[pix_fmt]))
1986         return -1;
1987
1988     y_shift = pix_fmt_info[pix_fmt].y_chroma_shift;
1989     x_shift = pix_fmt_info[pix_fmt].x_chroma_shift;
1990
1991     dst->data[0] = src->data[0] + (top_band * src->linesize[0]) + left_band;
1992     dst->data[1] = src->data[1] + ((top_band >> y_shift) * src->linesize[1]) + (left_band >> x_shift);
1993     dst->data[2] = src->data[2] + ((top_band >> y_shift) * src->linesize[2]) + (left_band >> x_shift);
1994
1995     dst->linesize[0] = src->linesize[0];
1996     dst->linesize[1] = src->linesize[1];
1997     dst->linesize[2] = src->linesize[2];
1998     return 0;
1999 }
2000
2001 /**
2002  * Pad image
2003  */
2004 int img_pad(AVPicture *dst, const AVPicture *src, int height, int width, int pix_fmt,
2005             int padtop, int padbottom, int padleft, int padright, int *color)
2006 {
2007     uint8_t *optr, *iptr;
2008     int y_shift;
2009     int x_shift;
2010     int yheight;
2011     int i, y;
2012
2013     if (pix_fmt < 0 || pix_fmt >= PIX_FMT_NB || !is_yuv_planar(&pix_fmt_info[pix_fmt]))
2014         return -1;
2015
2016     for (i = 0; i < 3; i++) {
2017         x_shift = i ? pix_fmt_info[pix_fmt].x_chroma_shift : 0;
2018         y_shift = i ? pix_fmt_info[pix_fmt].y_chroma_shift : 0;
2019
2020         if (padtop || padleft) {
2021             memset(dst->data[i], color[i], dst->linesize[i] * (padtop >> y_shift) + (padleft >> x_shift));
2022         }
2023
2024         if (padleft || padright || src) {
2025             if (src) { /* first line */
2026                 iptr = src->data[i];
2027                 optr = dst->data[i] + dst->linesize[i] * (padtop >> y_shift) + (padleft >> x_shift);
2028                 memcpy(optr, iptr, src->linesize[i]);
2029                 iptr += src->linesize[i];
2030             }
2031             optr = dst->data[i] + dst->linesize[i] * (padtop >> y_shift) + (dst->linesize[i] - (padright >> x_shift));
2032             yheight = (height - 1 - (padtop + padbottom)) >> y_shift;
2033             for (y = 0; y < yheight; y++) {
2034                 memset(optr, color[i], (padleft + padright) >> x_shift);
2035                 if (src) {
2036                     memcpy(optr + ((padleft + padright) >> x_shift), iptr, src->linesize[i]);
2037                     iptr += src->linesize[i];
2038                 }
2039                 optr += dst->linesize[i];
2040             }
2041         }
2042
2043         if (padbottom || padright) {
2044             optr = dst->data[i] + dst->linesize[i] * ((height - padbottom) >> y_shift) - (padright >> x_shift);
2045             memset(optr, color[i], dst->linesize[i] * (padbottom >> y_shift) + (padright >> x_shift));
2046         }
2047     }
2048     return 0;
2049 }
2050
2051 /* XXX: always use linesize. Return -1 if not supported */
2052 int img_convert(AVPicture *dst, int dst_pix_fmt,
2053                 const AVPicture *src, int src_pix_fmt,
2054                 int src_width, int src_height)
2055 {
2056     static int inited;
2057     int i, ret, dst_width, dst_height, int_pix_fmt;
2058     const PixFmtInfo *src_pix, *dst_pix;
2059     const ConvertEntry *ce;
2060     AVPicture tmp1, *tmp = &tmp1;
2061
2062     if (src_pix_fmt < 0 || src_pix_fmt >= PIX_FMT_NB ||
2063         dst_pix_fmt < 0 || dst_pix_fmt >= PIX_FMT_NB)
2064         return -1;
2065     if (src_width <= 0 || src_height <= 0)
2066         return 0;
2067
2068     if (!inited) {
2069         inited = 1;
2070         img_convert_init();
2071     }
2072
2073     dst_width = src_width;
2074     dst_height = src_height;
2075
2076     dst_pix = &pix_fmt_info[dst_pix_fmt];
2077     src_pix = &pix_fmt_info[src_pix_fmt];
2078     if (src_pix_fmt == dst_pix_fmt) {
2079         /* no conversion needed: just copy */
2080         img_copy(dst, src, dst_pix_fmt, dst_width, dst_height);
2081         return 0;
2082     }
2083
2084     ce = &convert_table[src_pix_fmt][dst_pix_fmt];
2085     if (ce->convert) {
2086         /* specific conversion routine */
2087         ce->convert(dst, src, dst_width, dst_height);
2088         return 0;
2089     }
2090
2091     /* gray to YUV */
2092     if (is_yuv_planar(dst_pix) &&
2093         src_pix_fmt == PIX_FMT_GRAY8) {
2094         int w, h, y;
2095         uint8_t *d;
2096
2097         if (dst_pix->color_type == FF_COLOR_YUV_JPEG) {
2098             ff_img_copy_plane(dst->data[0], dst->linesize[0],
2099                      src->data[0], src->linesize[0],
2100                      dst_width, dst_height);
2101         } else {
2102             img_apply_table(dst->data[0], dst->linesize[0],
2103                             src->data[0], src->linesize[0],
2104                             dst_width, dst_height,
2105                             y_jpeg_to_ccir);
2106         }
2107         /* fill U and V with 128 */
2108         w = dst_width;
2109         h = dst_height;
2110         w >>= dst_pix->x_chroma_shift;
2111         h >>= dst_pix->y_chroma_shift;
2112         for(i = 1; i <= 2; i++) {
2113             d = dst->data[i];
2114             for(y = 0; y< h; y++) {
2115                 memset(d, 128, w);
2116                 d += dst->linesize[i];
2117             }
2118         }
2119         return 0;
2120     }
2121
2122     /* YUV to gray */
2123     if (is_yuv_planar(src_pix) &&
2124         dst_pix_fmt == PIX_FMT_GRAY8) {
2125         if (src_pix->color_type == FF_COLOR_YUV_JPEG) {
2126             ff_img_copy_plane(dst->data[0], dst->linesize[0],
2127                      src->data[0], src->linesize[0],
2128                      dst_width, dst_height);
2129         } else {
2130             img_apply_table(dst->data[0], dst->linesize[0],
2131                             src->data[0], src->linesize[0],
2132                             dst_width, dst_height,
2133                             y_ccir_to_jpeg);
2134         }
2135         return 0;
2136     }
2137
2138     /* YUV to YUV planar */
2139     if (is_yuv_planar(dst_pix) && is_yuv_planar(src_pix)) {
2140         int x_shift, y_shift, w, h, xy_shift;
2141         void (*resize_func)(uint8_t *dst, int dst_wrap,
2142                             const uint8_t *src, int src_wrap,
2143                             int width, int height);
2144
2145         /* compute chroma size of the smallest dimensions */
2146         w = dst_width;
2147         h = dst_height;
2148         if (dst_pix->x_chroma_shift >= src_pix->x_chroma_shift)
2149             w >>= dst_pix->x_chroma_shift;
2150         else
2151             w >>= src_pix->x_chroma_shift;
2152         if (dst_pix->y_chroma_shift >= src_pix->y_chroma_shift)
2153             h >>= dst_pix->y_chroma_shift;
2154         else
2155             h >>= src_pix->y_chroma_shift;
2156
2157         x_shift = (dst_pix->x_chroma_shift - src_pix->x_chroma_shift);
2158         y_shift = (dst_pix->y_chroma_shift - src_pix->y_chroma_shift);
2159         xy_shift = ((x_shift & 0xf) << 4) | (y_shift & 0xf);
2160         /* there must be filters for conversion at least from and to
2161            YUV444 format */
2162         switch(xy_shift) {
2163         case 0x00:
2164             resize_func = ff_img_copy_plane;
2165             break;
2166         case 0x10:
2167             resize_func = shrink21;
2168             break;
2169         case 0x20:
2170             resize_func = shrink41;
2171             break;
2172         case 0x01:
2173             resize_func = shrink12;
2174             break;
2175         case 0x11:
2176             resize_func = ff_shrink22;
2177             break;
2178         case 0x22:
2179             resize_func = ff_shrink44;
2180             break;
2181         case 0xf0:
2182             resize_func = grow21;
2183             break;
2184         case 0xe0:
2185             resize_func = grow41;
2186             break;
2187         case 0xff:
2188             resize_func = grow22;
2189             break;
2190         case 0xee:
2191             resize_func = grow44;
2192             break;
2193         case 0xf1:
2194             resize_func = conv411;
2195             break;
2196         default:
2197             /* currently not handled */
2198             goto no_chroma_filter;
2199         }
2200
2201         ff_img_copy_plane(dst->data[0], dst->linesize[0],
2202                        src->data[0], src->linesize[0],
2203                        dst_width, dst_height);
2204
2205         for(i = 1;i <= 2; i++)
2206             resize_func(dst->data[i], dst->linesize[i],
2207                         src->data[i], src->linesize[i],
2208                         dst_width>>dst_pix->x_chroma_shift, dst_height>>dst_pix->y_chroma_shift);
2209         /* if yuv color space conversion is needed, we do it here on
2210            the destination image */
2211         if (dst_pix->color_type != src_pix->color_type) {
2212             const uint8_t *y_table, *c_table;
2213             if (dst_pix->color_type == FF_COLOR_YUV) {
2214                 y_table = y_jpeg_to_ccir;
2215                 c_table = c_jpeg_to_ccir;
2216             } else {
2217                 y_table = y_ccir_to_jpeg;
2218                 c_table = c_ccir_to_jpeg;
2219             }
2220             img_apply_table(dst->data[0], dst->linesize[0],
2221                             dst->data[0], dst->linesize[0],
2222                             dst_width, dst_height,
2223                             y_table);
2224
2225             for(i = 1;i <= 2; i++)
2226                 img_apply_table(dst->data[i], dst->linesize[i],
2227                                 dst->data[i], dst->linesize[i],
2228                                 dst_width>>dst_pix->x_chroma_shift,
2229                                 dst_height>>dst_pix->y_chroma_shift,
2230                                 c_table);
2231         }
2232         return 0;
2233     }
2234  no_chroma_filter:
2235
2236     /* try to use an intermediate format */
2237     if (src_pix_fmt == PIX_FMT_YUV422 ||
2238         dst_pix_fmt == PIX_FMT_YUV422) {
2239         /* specific case: convert to YUV422P first */
2240         int_pix_fmt = PIX_FMT_YUV422P;
2241     } else if (src_pix_fmt == PIX_FMT_UYVY422 ||
2242         dst_pix_fmt == PIX_FMT_UYVY422) {
2243         /* specific case: convert to YUV422P first */
2244         int_pix_fmt = PIX_FMT_YUV422P;
2245     } else if (src_pix_fmt == PIX_FMT_UYVY411 ||
2246         dst_pix_fmt == PIX_FMT_UYVY411) {
2247         /* specific case: convert to YUV411P first */
2248         int_pix_fmt = PIX_FMT_YUV411P;
2249     } else if ((src_pix->color_type == FF_COLOR_GRAY &&
2250                 src_pix_fmt != PIX_FMT_GRAY8) ||
2251                (dst_pix->color_type == FF_COLOR_GRAY &&
2252                 dst_pix_fmt != PIX_FMT_GRAY8)) {
2253         /* gray8 is the normalized format */
2254         int_pix_fmt = PIX_FMT_GRAY8;
2255     } else if ((is_yuv_planar(src_pix) &&
2256                 src_pix_fmt != PIX_FMT_YUV444P &&
2257                 src_pix_fmt != PIX_FMT_YUVJ444P)) {
2258         /* yuv444 is the normalized format */
2259         if (src_pix->color_type == FF_COLOR_YUV_JPEG)
2260             int_pix_fmt = PIX_FMT_YUVJ444P;
2261         else
2262             int_pix_fmt = PIX_FMT_YUV444P;
2263     } else if ((is_yuv_planar(dst_pix) &&
2264                 dst_pix_fmt != PIX_FMT_YUV444P &&
2265                 dst_pix_fmt != PIX_FMT_YUVJ444P)) {
2266         /* yuv444 is the normalized format */
2267         if (dst_pix->color_type == FF_COLOR_YUV_JPEG)
2268             int_pix_fmt = PIX_FMT_YUVJ444P;
2269         else
2270             int_pix_fmt = PIX_FMT_YUV444P;
2271     } else {
2272         /* the two formats are rgb or gray8 or yuv[j]444p */
2273         if (src_pix->is_alpha && dst_pix->is_alpha)
2274             int_pix_fmt = PIX_FMT_RGBA32;
2275         else
2276             int_pix_fmt = PIX_FMT_RGB24;
2277     }
2278     if (avpicture_alloc(tmp, int_pix_fmt, dst_width, dst_height) < 0)
2279         return -1;
2280     ret = -1;
2281     if (img_convert(tmp, int_pix_fmt,
2282                     src, src_pix_fmt, src_width, src_height) < 0)
2283         goto fail1;
2284     if (img_convert(dst, dst_pix_fmt,
2285                     tmp, int_pix_fmt, dst_width, dst_height) < 0)
2286         goto fail1;
2287     ret = 0;
2288  fail1:
2289     avpicture_free(tmp);
2290     return ret;
2291 }
2292
2293 /* NOTE: we scan all the pixels to have an exact information */
2294 static int get_alpha_info_pal8(const AVPicture *src, int width, int height)
2295 {
2296     const unsigned char *p;
2297     int src_wrap, ret, x, y;
2298     unsigned int a;
2299     uint32_t *palette = (uint32_t *)src->data[1];
2300
2301     p = src->data[0];
2302     src_wrap = src->linesize[0] - width;
2303     ret = 0;
2304     for(y=0;y<height;y++) {
2305         for(x=0;x<width;x++) {
2306             a = palette[p[0]] >> 24;
2307             if (a == 0x00) {
2308                 ret |= FF_ALPHA_TRANSP;
2309             } else if (a != 0xff) {
2310                 ret |= FF_ALPHA_SEMI_TRANSP;
2311             }
2312             p++;
2313         }
2314         p += src_wrap;
2315     }
2316     return ret;
2317 }
2318
2319 /**
2320  * Tell if an image really has transparent alpha values.
2321  * @return ored mask of FF_ALPHA_xxx constants
2322  */
2323 int img_get_alpha_info(const AVPicture *src,
2324                        int pix_fmt, int width, int height)
2325 {
2326     const PixFmtInfo *pf = &pix_fmt_info[pix_fmt];
2327     int ret;
2328
2329     pf = &pix_fmt_info[pix_fmt];
2330     /* no alpha can be represented in format */
2331     if (!pf->is_alpha)
2332         return 0;
2333     switch(pix_fmt) {
2334     case PIX_FMT_RGBA32:
2335         ret = get_alpha_info_rgba32(src, width, height);
2336         break;
2337     case PIX_FMT_RGB555:
2338         ret = get_alpha_info_rgb555(src, width, height);
2339         break;
2340     case PIX_FMT_PAL8:
2341         ret = get_alpha_info_pal8(src, width, height);
2342         break;
2343     default:
2344         /* we do not know, so everything is indicated */
2345         ret = FF_ALPHA_TRANSP | FF_ALPHA_SEMI_TRANSP;
2346         break;
2347     }
2348     return ret;
2349 }
2350
2351 #ifdef HAVE_MMX
2352 #define DEINT_INPLACE_LINE_LUM \
2353                     movd_m2r(lum_m4[0],mm0);\
2354                     movd_m2r(lum_m3[0],mm1);\
2355                     movd_m2r(lum_m2[0],mm2);\
2356                     movd_m2r(lum_m1[0],mm3);\
2357                     movd_m2r(lum[0],mm4);\
2358                     punpcklbw_r2r(mm7,mm0);\
2359                     movd_r2m(mm2,lum_m4[0]);\
2360                     punpcklbw_r2r(mm7,mm1);\
2361                     punpcklbw_r2r(mm7,mm2);\
2362                     punpcklbw_r2r(mm7,mm3);\
2363                     punpcklbw_r2r(mm7,mm4);\
2364                     paddw_r2r(mm3,mm1);\
2365                     psllw_i2r(1,mm2);\
2366                     paddw_r2r(mm4,mm0);\
2367                     psllw_i2r(2,mm1);\
2368                     paddw_r2r(mm6,mm2);\
2369                     paddw_r2r(mm2,mm1);\
2370                     psubusw_r2r(mm0,mm1);\
2371                     psrlw_i2r(3,mm1);\
2372                     packuswb_r2r(mm7,mm1);\
2373                     movd_r2m(mm1,lum_m2[0]);
2374
2375 #define DEINT_LINE_LUM \
2376                     movd_m2r(lum_m4[0],mm0);\
2377                     movd_m2r(lum_m3[0],mm1);\
2378                     movd_m2r(lum_m2[0],mm2);\
2379                     movd_m2r(lum_m1[0],mm3);\
2380                     movd_m2r(lum[0],mm4);\
2381                     punpcklbw_r2r(mm7,mm0);\
2382                     punpcklbw_r2r(mm7,mm1);\
2383                     punpcklbw_r2r(mm7,mm2);\
2384                     punpcklbw_r2r(mm7,mm3);\
2385                     punpcklbw_r2r(mm7,mm4);\
2386                     paddw_r2r(mm3,mm1);\
2387                     psllw_i2r(1,mm2);\
2388                     paddw_r2r(mm4,mm0);\
2389                     psllw_i2r(2,mm1);\
2390                     paddw_r2r(mm6,mm2);\
2391                     paddw_r2r(mm2,mm1);\
2392                     psubusw_r2r(mm0,mm1);\
2393                     psrlw_i2r(3,mm1);\
2394                     packuswb_r2r(mm7,mm1);\
2395                     movd_r2m(mm1,dst[0]);
2396 #endif
2397
2398 /* filter parameters: [-1 4 2 4 -1] // 8 */
2399 static void deinterlace_line(uint8_t *dst,
2400                              const uint8_t *lum_m4, const uint8_t *lum_m3,
2401                              const uint8_t *lum_m2, const uint8_t *lum_m1,
2402                              const uint8_t *lum,
2403                              int size)
2404 {
2405 #ifndef HAVE_MMX
2406     uint8_t *cm = cropTbl + MAX_NEG_CROP;
2407     int sum;
2408
2409     for(;size > 0;size--) {
2410         sum = -lum_m4[0];
2411         sum += lum_m3[0] << 2;
2412         sum += lum_m2[0] << 1;
2413         sum += lum_m1[0] << 2;
2414         sum += -lum[0];
2415         dst[0] = cm[(sum + 4) >> 3];
2416         lum_m4++;
2417         lum_m3++;
2418         lum_m2++;
2419         lum_m1++;
2420         lum++;
2421         dst++;
2422     }
2423 #else
2424
2425     {
2426         mmx_t rounder;
2427         rounder.uw[0]=4;
2428         rounder.uw[1]=4;
2429         rounder.uw[2]=4;
2430         rounder.uw[3]=4;
2431         pxor_r2r(mm7,mm7);
2432         movq_m2r(rounder,mm6);
2433     }
2434     for (;size > 3; size-=4) {
2435         DEINT_LINE_LUM
2436         lum_m4+=4;
2437         lum_m3+=4;
2438         lum_m2+=4;
2439         lum_m1+=4;
2440         lum+=4;
2441         dst+=4;
2442     }
2443 #endif
2444 }
2445 static void deinterlace_line_inplace(uint8_t *lum_m4, uint8_t *lum_m3, uint8_t *lum_m2, uint8_t *lum_m1, uint8_t *lum,
2446                              int size)
2447 {
2448 #ifndef HAVE_MMX
2449     uint8_t *cm = cropTbl + MAX_NEG_CROP;
2450     int sum;
2451
2452     for(;size > 0;size--) {
2453         sum = -lum_m4[0];
2454         sum += lum_m3[0] << 2;
2455         sum += lum_m2[0] << 1;
2456         lum_m4[0]=lum_m2[0];
2457         sum += lum_m1[0] << 2;
2458         sum += -lum[0];
2459         lum_m2[0] = cm[(sum + 4) >> 3];
2460         lum_m4++;
2461         lum_m3++;
2462         lum_m2++;
2463         lum_m1++;
2464         lum++;
2465     }
2466 #else
2467
2468     {
2469         mmx_t rounder;
2470         rounder.uw[0]=4;
2471         rounder.uw[1]=4;
2472         rounder.uw[2]=4;
2473         rounder.uw[3]=4;
2474         pxor_r2r(mm7,mm7);
2475         movq_m2r(rounder,mm6);
2476     }
2477     for (;size > 3; size-=4) {
2478         DEINT_INPLACE_LINE_LUM
2479         lum_m4+=4;
2480         lum_m3+=4;
2481         lum_m2+=4;
2482         lum_m1+=4;
2483         lum+=4;
2484     }
2485 #endif
2486 }
2487
2488 /* deinterlacing : 2 temporal taps, 3 spatial taps linear filter. The
2489    top field is copied as is, but the bottom field is deinterlaced
2490    against the top field. */
2491 static void deinterlace_bottom_field(uint8_t *dst, int dst_wrap,
2492                                     const uint8_t *src1, int src_wrap,
2493                                     int width, int height)
2494 {
2495     const uint8_t *src_m2, *src_m1, *src_0, *src_p1, *src_p2;
2496     int y;
2497
2498     src_m2 = src1;
2499     src_m1 = src1;
2500     src_0=&src_m1[src_wrap];
2501     src_p1=&src_0[src_wrap];
2502     src_p2=&src_p1[src_wrap];
2503     for(y=0;y<(height-2);y+=2) {
2504         memcpy(dst,src_m1,width);
2505         dst += dst_wrap;
2506         deinterlace_line(dst,src_m2,src_m1,src_0,src_p1,src_p2,width);
2507         src_m2 = src_0;
2508         src_m1 = src_p1;
2509         src_0 = src_p2;
2510         src_p1 += 2*src_wrap;
2511         src_p2 += 2*src_wrap;
2512         dst += dst_wrap;
2513     }
2514     memcpy(dst,src_m1,width);
2515     dst += dst_wrap;
2516     /* do last line */
2517     deinterlace_line(dst,src_m2,src_m1,src_0,src_0,src_0,width);
2518 }
2519
2520 static void deinterlace_bottom_field_inplace(uint8_t *src1, int src_wrap,
2521                                              int width, int height)
2522 {
2523     uint8_t *src_m1, *src_0, *src_p1, *src_p2;
2524     int y;
2525     uint8_t *buf;
2526     buf = (uint8_t*)av_malloc(width);
2527
2528     src_m1 = src1;
2529     memcpy(buf,src_m1,width);
2530     src_0=&src_m1[src_wrap];
2531     src_p1=&src_0[src_wrap];
2532     src_p2=&src_p1[src_wrap];
2533     for(y=0;y<(height-2);y+=2) {
2534         deinterlace_line_inplace(buf,src_m1,src_0,src_p1,src_p2,width);
2535         src_m1 = src_p1;
2536         src_0 = src_p2;
2537         src_p1 += 2*src_wrap;
2538         src_p2 += 2*src_wrap;
2539     }
2540     /* do last line */
2541     deinterlace_line_inplace(buf,src_m1,src_0,src_0,src_0,width);
2542     av_free(buf);
2543 }
2544
2545
2546 /* deinterlace - if not supported return -1 */
2547 int avpicture_deinterlace(AVPicture *dst, const AVPicture *src,
2548                           int pix_fmt, int width, int height)
2549 {
2550     int i;
2551
2552     if (pix_fmt != PIX_FMT_YUV420P &&
2553         pix_fmt != PIX_FMT_YUV422P &&
2554         pix_fmt != PIX_FMT_YUV444P &&
2555         pix_fmt != PIX_FMT_YUV411P)
2556         return -1;
2557     if ((width & 3) != 0 || (height & 3) != 0)
2558         return -1;
2559
2560     for(i=0;i<3;i++) {
2561         if (i == 1) {
2562             switch(pix_fmt) {
2563             case PIX_FMT_YUV420P:
2564                 width >>= 1;
2565                 height >>= 1;
2566                 break;
2567             case PIX_FMT_YUV422P:
2568                 width >>= 1;
2569                 break;
2570             case PIX_FMT_YUV411P:
2571                 width >>= 2;
2572                 break;
2573             default:
2574                 break;
2575             }
2576         }
2577         if (src == dst) {
2578             deinterlace_bottom_field_inplace(dst->data[i], dst->linesize[i],
2579                                  width, height);
2580         } else {
2581             deinterlace_bottom_field(dst->data[i],dst->linesize[i],
2582                                         src->data[i], src->linesize[i],
2583                                         width, height);
2584         }
2585     }
2586 #ifdef HAVE_MMX
2587     emms();
2588 #endif
2589     return 0;
2590 }
2591
2592 #undef FIX