imgutils: Fix a typo in avcodec_get_pix_fmt_loss
[libav.git] / libavcodec / imgconvert.c
CommitLineData
de6d9b64 1/*
f1ea5c2a 2 * Misc image conversion routines
406792e7 3 * Copyright (c) 2001, 2002, 2003 Fabrice Bellard
de6d9b64 4 *
2912e87a 5 * This file is part of Libav.
b78e7197 6 *
2912e87a 7 * Libav is free software; you can redistribute it and/or
ff4ec49e
FB
8 * modify it under the terms of the GNU Lesser General Public
9 * License as published by the Free Software Foundation; either
b78e7197 10 * version 2.1 of the License, or (at your option) any later version.
de6d9b64 11 *
2912e87a 12 * Libav is distributed in the hope that it will be useful,
de6d9b64 13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
ff4ec49e
FB
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15 * Lesser General Public License for more details.
de6d9b64 16 *
ff4ec49e 17 * You should have received a copy of the GNU Lesser General Public
2912e87a 18 * License along with Libav; if not, write to the Free Software
5509bffa 19 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
de6d9b64 20 */
983e3246
MN
21
22/**
ba87f080 23 * @file
f1ea5c2a 24 * misc image conversion routines
983e3246
MN
25 */
26
c50c0bc8
FB
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 */
983e3246 32
de6d9b64 33#include "avcodec.h"
35685a3c 34#include "imgconvert.h"
4443c0e9 35#include "internal.h"
f9279ee7 36#include "mathops.h"
2b4abbd6 37#include "libavutil/colorspace.h"
1d9c2dc8 38#include "libavutil/common.h"
8e861e1b 39#include "libavutil/pixdesc.h"
737eb597 40#include "libavutil/imgutils.h"
de6d9b64 41
716d413c 42void avcodec_get_chroma_sub_sample(enum AVPixelFormat pix_fmt, int *h_shift, int *v_shift)
524c6b63 43{
50ba57e0
AK
44 const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(pix_fmt);
45 *h_shift = desc->log2_chroma_w;
46 *v_shift = desc->log2_chroma_h;
524c6b63
FB
47}
48
ae3822bc
LB
49static int is_gray(const AVPixFmtDescriptor *desc)
50{
e6c4ac7b 51 return desc->nb_components - (desc->flags & AV_PIX_FMT_FLAG_ALPHA) == 1;
ae3822bc 52}
2a877875 53
ae3822bc
LB
54int avcodec_get_pix_fmt_loss(enum AVPixelFormat dst_pix_fmt,
55 enum AVPixelFormat src_pix_fmt,
b6147995
FB
56 int has_alpha)
57{
50ba57e0
AK
58 const AVPixFmtDescriptor *src_desc = av_pix_fmt_desc_get(src_pix_fmt);
59 const AVPixFmtDescriptor *dst_desc = av_pix_fmt_desc_get(dst_pix_fmt);
ae3822bc
LB
60 int loss, i, nb_components = FFMIN(src_desc->nb_components,
61 dst_desc->nb_components);
b6147995
FB
62
63 /* compute loss */
64 loss = 0;
ae3822bc
LB
65
66 if (dst_pix_fmt == src_pix_fmt)
67 return 0;
68
69 for (i = 0; i < nb_components; i++)
70 if (src_desc->comp[i].depth_minus1 > dst_desc->comp[i].depth_minus1)
71 loss |= FF_LOSS_DEPTH;
72
14b903f3
SS
73 if (dst_desc->log2_chroma_w > src_desc->log2_chroma_w ||
74 dst_desc->log2_chroma_h > src_desc->log2_chroma_h)
b6147995 75 loss |= FF_LOSS_RESOLUTION;
ae3822bc 76
e6c4ac7b 77 if ((src_desc->flags & AV_PIX_FMT_FLAG_RGB) != (dst_desc->flags & AV_PIX_FMT_FLAG_RGB))
ae3822bc
LB
78 loss |= FF_LOSS_COLORSPACE;
79
e6c4ac7b 80 if (has_alpha && !(dst_desc->flags & AV_PIX_FMT_FLAG_ALPHA) &&
47b447aa 81 (src_desc->flags & AV_PIX_FMT_FLAG_ALPHA))
b6147995 82 loss |= FF_LOSS_ALPHA;
b6147995 83
ae3822bc
LB
84 if (dst_pix_fmt == AV_PIX_FMT_PAL8 && !is_gray(src_desc))
85 return loss | FF_LOSS_COLORQUANT;
b6147995 86
ae3822bc
LB
87 if (src_desc->nb_components > dst_desc->nb_components)
88 if (is_gray(dst_desc))
89 loss |= FF_LOSS_CHROMA;
90
91 return loss;
b6147995
FB
92}
93
716d413c
AK
94static enum AVPixelFormat avcodec_find_best_pix_fmt1(enum AVPixelFormat *pix_fmt_list,
95 enum AVPixelFormat src_pix_fmt,
b6147995
FB
96 int has_alpha,
97 int loss_mask)
98{
a61ec8e7 99 int dist, i, loss, min_dist;
716d413c 100 enum AVPixelFormat dst_pix_fmt;
b6147995
FB
101
102 /* find exact color match with smallest size */
716d413c 103 dst_pix_fmt = AV_PIX_FMT_NONE;
b6147995 104 min_dist = 0x7fffffff;
ac29054f 105 i = 0;
716d413c
AK
106 while (pix_fmt_list[i] != AV_PIX_FMT_NONE) {
107 enum AVPixelFormat pix_fmt = pix_fmt_list[i];
ac29054f 108
716d413c 109 if (i > AV_PIX_FMT_NB) {
ac29054f
JG
110 av_log(NULL, AV_LOG_ERROR, "Pixel format list longer than expected, "
111 "it is either not properly terminated or contains duplicates\n");
716d413c 112 return AV_PIX_FMT_NONE;
ac29054f
JG
113 }
114
115 loss = avcodec_get_pix_fmt_loss(pix_fmt, src_pix_fmt, has_alpha) & loss_mask;
116 if (loss == 0) {
ae3822bc 117 dist = av_get_bits_per_pixel(av_pix_fmt_desc_get(pix_fmt));
ac29054f
JG
118 if (dist < min_dist) {
119 min_dist = dist;
120 dst_pix_fmt = pix_fmt;
b6147995
FB
121 }
122 }
ac29054f 123 i++;
b6147995
FB
124 }
125 return dst_pix_fmt;
126}
127
716d413c
AK
128enum AVPixelFormat avcodec_find_best_pix_fmt2(enum AVPixelFormat *pix_fmt_list,
129 enum AVPixelFormat src_pix_fmt,
ac29054f
JG
130 int has_alpha, int *loss_ptr)
131{
716d413c 132 enum AVPixelFormat dst_pix_fmt;
a61ec8e7 133 int loss_mask, i;
b6147995
FB
134 static const int loss_mask_order[] = {
135 ~0, /* no loss first */
136 ~FF_LOSS_ALPHA,
137 ~FF_LOSS_RESOLUTION,
138 ~(FF_LOSS_COLORSPACE | FF_LOSS_RESOLUTION),
139 ~FF_LOSS_COLORQUANT,
140 ~FF_LOSS_DEPTH,
141 0,
142 };
143
144 /* try with successive loss */
145 i = 0;
146 for(;;) {
147 loss_mask = loss_mask_order[i++];
ac29054f 148 dst_pix_fmt = avcodec_find_best_pix_fmt1(pix_fmt_list, src_pix_fmt,
b6147995
FB
149 has_alpha, loss_mask);
150 if (dst_pix_fmt >= 0)
151 goto found;
152 if (loss_mask == 0)
153 break;
154 }
716d413c 155 return AV_PIX_FMT_NONE;
b6147995
FB
156 found:
157 if (loss_ptr)
158 *loss_ptr = avcodec_get_pix_fmt_loss(dst_pix_fmt, src_pix_fmt, has_alpha);
159 return dst_pix_fmt;
160}
161
85c242d8 162/* 2x2 -> 1x1 */
54009d42 163void ff_shrink22(uint8_t *dst, int dst_wrap,
e352ff08 164 const uint8_t *src, int src_wrap,
85c242d8
FB
165 int width, int height)
166{
167 int w;
e352ff08
FB
168 const uint8_t *s1, *s2;
169 uint8_t *d;
85c242d8
FB
170
171 for(;height > 0; height--) {
172 s1 = src;
173 s2 = s1 + src_wrap;
174 d = dst;
175 for(w = width;w >= 4; w-=4) {
0a9ad8d1
FB
176 d[0] = (s1[0] + s1[1] + s2[0] + s2[1] + 2) >> 2;
177 d[1] = (s1[2] + s1[3] + s2[2] + s2[3] + 2) >> 2;
178 d[2] = (s1[4] + s1[5] + s2[4] + s2[5] + 2) >> 2;
179 d[3] = (s1[6] + s1[7] + s2[6] + s2[7] + 2) >> 2;
85c242d8
FB
180 s1 += 8;
181 s2 += 8;
182 d += 4;
183 }
184 for(;w > 0; w--) {
0a9ad8d1 185 d[0] = (s1[0] + s1[1] + s2[0] + s2[1] + 2) >> 2;
85c242d8
FB
186 s1 += 2;
187 s2 += 2;
188 d++;
189 }
190 src += 2 * src_wrap;
191 dst += dst_wrap;
192 }
193}
194
e352ff08 195/* 4x4 -> 1x1 */
54009d42 196void ff_shrink44(uint8_t *dst, int dst_wrap,
e352ff08 197 const uint8_t *src, int src_wrap,
6742d95d
FR
198 int width, int height)
199{
200 int w;
e352ff08
FB
201 const uint8_t *s1, *s2, *s3, *s4;
202 uint8_t *d;
6742d95d
FR
203
204 for(;height > 0; height--) {
205 s1 = src;
e352ff08
FB
206 s2 = s1 + src_wrap;
207 s3 = s2 + src_wrap;
208 s4 = s3 + src_wrap;
6742d95d 209 d = dst;
e352ff08
FB
210 for(w = width;w > 0; w--) {
211 d[0] = (s1[0] + s1[1] + s1[2] + s1[3] +
212 s2[0] + s2[1] + s2[2] + s2[3] +
213 s3[0] + s3[1] + s3[2] + s3[3] +
214 s4[0] + s4[1] + s4[2] + s4[3] + 8) >> 4;
215 s1 += 4;
216 s2 += 4;
217 s3 += 4;
218 s4 += 4;
6742d95d
FR
219 d++;
220 }
e352ff08
FB
221 src += 4 * src_wrap;
222 dst += dst_wrap;
223 }
224}
225
54009d42
MN
226/* 8x8 -> 1x1 */
227void ff_shrink88(uint8_t *dst, int dst_wrap,
228 const uint8_t *src, int src_wrap,
229 int width, int height)
230{
231 int w, i;
232
233 for(;height > 0; height--) {
234 for(w = width;w > 0; w--) {
235 int tmp=0;
236 for(i=0; i<8; i++){
237 tmp += src[0] + src[1] + src[2] + src[3] + src[4] + src[5] + src[6] + src[7];
238 src += src_wrap;
239 }
240 *(dst++) = (tmp + 32)>>6;
241 src += 8 - 8*src_wrap;
242 }
243 src += 8*src_wrap - 8*width;
244 dst += dst_wrap - width;
245 }
246}
247
c50c0bc8 248/* return true if yuv planar */
ae3822bc 249static inline int is_yuv_planar(const AVPixFmtDescriptor *desc)
c50c0bc8 250{
e6c4ac7b
AK
251 return (!(desc->flags & AV_PIX_FMT_FLAG_RGB) &&
252 (desc->flags & AV_PIX_FMT_FLAG_PLANAR));
c50c0bc8
FB
253}
254
636d6a4a 255int av_picture_crop(AVPicture *dst, const AVPicture *src,
716d413c 256 enum AVPixelFormat pix_fmt, int top_band, int left_band)
f2651e7a 257{
50ba57e0 258 const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(pix_fmt);
f2651e7a
BC
259 int y_shift;
260 int x_shift;
261
ae3822bc 262 if (pix_fmt < 0 || pix_fmt >= AV_PIX_FMT_NB || !is_yuv_planar(desc))
f2651e7a
BC
263 return -1;
264
50ba57e0
AK
265 y_shift = desc->log2_chroma_h;
266 x_shift = desc->log2_chroma_w;
f2651e7a
BC
267
268 dst->data[0] = src->data[0] + (top_band * src->linesize[0]) + left_band;
269 dst->data[1] = src->data[1] + ((top_band >> y_shift) * src->linesize[1]) + (left_band >> x_shift);
270 dst->data[2] = src->data[2] + ((top_band >> y_shift) * src->linesize[2]) + (left_band >> x_shift);
271
272 dst->linesize[0] = src->linesize[0];
273 dst->linesize[1] = src->linesize[1];
274 dst->linesize[2] = src->linesize[2];
275 return 0;
276}
277
636d6a4a 278int av_picture_pad(AVPicture *dst, const AVPicture *src, int height, int width,
716d413c 279 enum AVPixelFormat pix_fmt, int padtop, int padbottom, int padleft, int padright,
6845801f 280 int *color)
5341c209 281{
50ba57e0 282 const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(pix_fmt);
79acfb0e 283 uint8_t *optr;
5341c209
LA
284 int y_shift;
285 int x_shift;
286 int yheight;
287 int i, y;
288
716d413c 289 if (pix_fmt < 0 || pix_fmt >= AV_PIX_FMT_NB ||
ae3822bc 290 !is_yuv_planar(desc)) return -1;
5341c209
LA
291
292 for (i = 0; i < 3; i++) {
50ba57e0
AK
293 x_shift = i ? desc->log2_chroma_w : 0;
294 y_shift = i ? desc->log2_chroma_h : 0;
5341c209
LA
295
296 if (padtop || padleft) {
6845801f
LB
297 memset(dst->data[i], color[i],
298 dst->linesize[i] * (padtop >> y_shift) + (padleft >> x_shift));
5341c209
LA
299 }
300
79acfb0e
LB
301 if (padleft || padright) {
302 optr = dst->data[i] + dst->linesize[i] * (padtop >> y_shift) +
303 (dst->linesize[i] - (padright >> x_shift));
304 yheight = (height - 1 - (padtop + padbottom)) >> y_shift;
305 for (y = 0; y < yheight; y++) {
306 memset(optr, color[i], (padleft + padright) >> x_shift);
307 optr += dst->linesize[i];
5341c209 308 }
79acfb0e
LB
309 }
310
311 if (src) { /* first line */
312 uint8_t *iptr = src->data[i];
313 optr = dst->data[i] + dst->linesize[i] * (padtop >> y_shift) +
314 (padleft >> x_shift);
1be97a21 315 memcpy(optr, iptr, (width - padleft - padright) >> x_shift);
79acfb0e 316 iptr += src->linesize[i];
6845801f
LB
317 optr = dst->data[i] + dst->linesize[i] * (padtop >> y_shift) +
318 (dst->linesize[i] - (padright >> x_shift));
5341c209
LA
319 yheight = (height - 1 - (padtop + padbottom)) >> y_shift;
320 for (y = 0; y < yheight; y++) {
321 memset(optr, color[i], (padleft + padright) >> x_shift);
79acfb0e 322 memcpy(optr + ((padleft + padright) >> x_shift), iptr,
1be97a21 323 (width - padleft - padright) >> x_shift);
79acfb0e 324 iptr += src->linesize[i];
5341c209
LA
325 optr += dst->linesize[i];
326 }
327 }
328
329 if (padbottom || padright) {
6845801f
LB
330 optr = dst->data[i] + dst->linesize[i] *
331 ((height - padbottom) >> y_shift) - (padright >> x_shift);
332 memset(optr, color[i],dst->linesize[i] *
333 (padbottom >> y_shift) + (padright >> x_shift));
5341c209
LA
334 }
335 }
336 return 0;
337}
338
54b298fe
RB
339#if FF_API_DEINTERLACE
340
78670fbf
DB
341#if HAVE_MMX_EXTERNAL
342#define deinterlace_line_inplace ff_deinterlace_line_inplace_mmx
343#define deinterlace_line ff_deinterlace_line_mmx
344#else
345#define deinterlace_line_inplace deinterlace_line_inplace_c
346#define deinterlace_line deinterlace_line_c
347
85c242d8 348/* filter parameters: [-1 4 2 4 -1] // 8 */
de4bc44a 349static void deinterlace_line_c(uint8_t *dst,
bb270c08
DB
350 const uint8_t *lum_m4, const uint8_t *lum_m3,
351 const uint8_t *lum_m2, const uint8_t *lum_m1,
352 const uint8_t *lum,
353 int size)
85c242d8 354{
05563cca 355 const uint8_t *cm = ff_crop_tab + MAX_NEG_CROP;
85c242d8 356 int sum;
85c242d8
FB
357
358 for(;size > 0;size--) {
5981f4e6
F
359 sum = -lum_m4[0];
360 sum += lum_m3[0] << 2;
361 sum += lum_m2[0] << 1;
362 sum += lum_m1[0] << 2;
363 sum += -lum[0];
85c242d8 364 dst[0] = cm[(sum + 4) >> 3];
5981f4e6
F
365 lum_m4++;
366 lum_m3++;
367 lum_m2++;
368 lum_m1++;
369 lum++;
85c242d8 370 dst++;
85c242d8 371 }
5981f4e6 372}
de4bc44a
VS
373
374static void deinterlace_line_inplace_c(uint8_t *lum_m4, uint8_t *lum_m3,
375 uint8_t *lum_m2, uint8_t *lum_m1,
376 uint8_t *lum, int size)
5981f4e6 377{
05563cca 378 const uint8_t *cm = ff_crop_tab + MAX_NEG_CROP;
5981f4e6
F
379 int sum;
380
381 for(;size > 0;size--) {
382 sum = -lum_m4[0];
383 sum += lum_m3[0] << 2;
384 sum += lum_m2[0] << 1;
385 lum_m4[0]=lum_m2[0];
386 sum += lum_m1[0] << 2;
387 sum += -lum[0];
388 lum_m2[0] = cm[(sum + 4) >> 3];
389 lum_m4++;
390 lum_m3++;
391 lum_m2++;
392 lum_m1++;
393 lum++;
394 }
85c242d8 395}
17337f54 396#endif /* !HAVE_MMX_EXTERNAL */
85c242d8
FB
397
398/* deinterlacing : 2 temporal taps, 3 spatial taps linear filter. The
399 top field is copied as is, but the bottom field is deinterlaced
400 against the top field. */
0c1a9eda 401static void deinterlace_bottom_field(uint8_t *dst, int dst_wrap,
da64ecc3 402 const uint8_t *src1, int src_wrap,
5981f4e6 403 int width, int height)
85c242d8 404{
da64ecc3 405 const uint8_t *src_m2, *src_m1, *src_0, *src_p1, *src_p2;
5981f4e6
F
406 int y;
407
408 src_m2 = src1;
409 src_m1 = src1;
410 src_0=&src_m1[src_wrap];
411 src_p1=&src_0[src_wrap];
412 src_p2=&src_p1[src_wrap];
413 for(y=0;y<(height-2);y+=2) {
414 memcpy(dst,src_m1,width);
85c242d8 415 dst += dst_wrap;
5981f4e6
F
416 deinterlace_line(dst,src_m2,src_m1,src_0,src_p1,src_p2,width);
417 src_m2 = src_0;
418 src_m1 = src_p1;
419 src_0 = src_p2;
420 src_p1 += 2*src_wrap;
421 src_p2 += 2*src_wrap;
85c242d8 422 dst += dst_wrap;
85c242d8 423 }
5981f4e6
F
424 memcpy(dst,src_m1,width);
425 dst += dst_wrap;
426 /* do last line */
427 deinterlace_line(dst,src_m2,src_m1,src_0,src_0,src_0,width);
428}
429
210461c0
VG
430static int deinterlace_bottom_field_inplace(uint8_t *src1, int src_wrap,
431 int width, int height)
5981f4e6 432{
0c1a9eda 433 uint8_t *src_m1, *src_0, *src_p1, *src_p2;
5981f4e6 434 int y;
0c1a9eda 435 uint8_t *buf;
44b0edda 436 buf = av_malloc(width);
210461c0
VG
437 if (!buf)
438 return AVERROR(ENOMEM);
5981f4e6
F
439
440 src_m1 = src1;
441 memcpy(buf,src_m1,width);
442 src_0=&src_m1[src_wrap];
443 src_p1=&src_0[src_wrap];
444 src_p2=&src_p1[src_wrap];
445 for(y=0;y<(height-2);y+=2) {
446 deinterlace_line_inplace(buf,src_m1,src_0,src_p1,src_p2,width);
447 src_m1 = src_p1;
448 src_0 = src_p2;
449 src_p1 += 2*src_wrap;
450 src_p2 += 2*src_wrap;
451 }
452 /* do last line */
453 deinterlace_line_inplace(buf,src_m1,src_0,src_0,src_0,width);
6000abfa 454 av_free(buf);
210461c0 455 return 0;
85c242d8
FB
456}
457
da64ecc3 458int avpicture_deinterlace(AVPicture *dst, const AVPicture *src,
716d413c 459 enum AVPixelFormat pix_fmt, int width, int height)
de6d9b64 460{
210461c0 461 int i, ret;
85c242d8 462
716d413c
AK
463 if (pix_fmt != AV_PIX_FMT_YUV420P &&
464 pix_fmt != AV_PIX_FMT_YUVJ420P &&
465 pix_fmt != AV_PIX_FMT_YUV422P &&
466 pix_fmt != AV_PIX_FMT_YUVJ422P &&
467 pix_fmt != AV_PIX_FMT_YUV444P &&
468 pix_fmt != AV_PIX_FMT_YUV411P &&
469 pix_fmt != AV_PIX_FMT_GRAY8)
85c242d8 470 return -1;
5981f4e6 471 if ((width & 3) != 0 || (height & 3) != 0)
85c242d8 472 return -1;
5981f4e6 473
85c242d8
FB
474 for(i=0;i<3;i++) {
475 if (i == 1) {
476 switch(pix_fmt) {
716d413c
AK
477 case AV_PIX_FMT_YUVJ420P:
478 case AV_PIX_FMT_YUV420P:
85c242d8
FB
479 width >>= 1;
480 height >>= 1;
481 break;
716d413c
AK
482 case AV_PIX_FMT_YUV422P:
483 case AV_PIX_FMT_YUVJ422P:
85c242d8
FB
484 width >>= 1;
485 break;
716d413c 486 case AV_PIX_FMT_YUV411P:
47017dd8
RS
487 width >>= 2;
488 break;
85c242d8
FB
489 default:
490 break;
491 }
716d413c 492 if (pix_fmt == AV_PIX_FMT_GRAY8) {
2a7feb18
AG
493 break;
494 }
85c242d8 495 }
5981f4e6 496 if (src == dst) {
210461c0
VG
497 ret = deinterlace_bottom_field_inplace(dst->data[i],
498 dst->linesize[i],
499 width, height);
500 if (ret < 0)
501 return ret;
5981f4e6
F
502 } else {
503 deinterlace_bottom_field(dst->data[i],dst->linesize[i],
504 src->data[i], src->linesize[i],
505 width, height);
506 }
de6d9b64 507 }
55ffe9df 508 emms_c();
85c242d8 509 return 0;
de6d9b64 510}
54b298fe
RB
511
512#endif /* FF_API_DEINTERLACE */