2 * Copyright (C) 2003 Michael Niedermayer <michaelni@gmx.at>
4 * This file is part of Libav.
6 * Libav is free software; you can redistribute it and/or
7 * modify it under the terms of the GNU Lesser General Public
8 * License as published by the Free Software Foundation; either
9 * version 2.1 of the License, or (at your option) any later version.
11 * Libav is distributed in the hope that it will be useful,
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 * Lesser General Public License for more details.
16 * You should have received a copy of the GNU Lesser General Public
17 * License along with Libav; if not, write to the Free Software
18 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
27 #undef HAVE_AV_CONFIG_H
28 #include "libavutil/imgutils.h"
29 #include "libavutil/mem.h"
30 #include "libavutil/avutil.h"
31 #include "libavutil/crc.h"
32 #include "libavutil/pixdesc.h"
33 #include "libavutil/lfg.h"
36 /* HACK Duplicated from swscale_internal.h.
37 * Should be removed when a cleaner pixel format system exists. */
40 || (x)==PIX_FMT_GRAY16BE \
41 || (x)==PIX_FMT_GRAY16LE \
43 #define hasChroma(x) (!( \
45 || (x)==PIX_FMT_MONOBLACK \
46 || (x)==PIX_FMT_MONOWHITE \
48 #define isALPHA(x) ( \
50 || (x)==PIX_FMT_BGR32_1 \
51 || (x)==PIX_FMT_RGB32 \
52 || (x)==PIX_FMT_RGB32_1 \
53 || (x)==PIX_FMT_YUVA420P \
56 static uint64_t getSSD(uint8_t *src1, uint8_t *src2, int stride1, int stride2, int w, int h)
61 //printf("%d %d\n", w, h);
65 int d= src1[x + y*stride1] - src2[x + y*stride2];
67 //printf("%d", abs(src1[x + y*stride1] - src2[x + y*stride2])/26 );
82 // test by ref -> src -> dst -> out & compare out against ref
84 static int doTest(uint8_t *ref[4], int refStride[4], int w, int h,
85 enum PixelFormat srcFormat, enum PixelFormat dstFormat,
86 int srcW, int srcH, int dstW, int dstH, int flags,
89 static enum PixelFormat cur_srcFormat;
90 static int cur_srcW, cur_srcH;
91 static uint8_t *src[4];
92 static int srcStride[4];
93 uint8_t *dst[4] = {0};
94 uint8_t *out[4] = {0};
97 uint64_t ssdY, ssdU=0, ssdV=0, ssdA=0;
98 struct SwsContext *dstContext = NULL, *outContext = NULL;
102 if (cur_srcFormat != srcFormat || cur_srcW != srcW || cur_srcH != srcH) {
103 struct SwsContext *srcContext = NULL;
106 for (p = 0; p < 4; p++)
109 av_image_fill_linesizes(srcStride, srcFormat, srcW);
110 for (p = 0; p < 4; p++) {
112 src[p] = av_mallocz(srcStride[p]*srcH+16);
113 if (srcStride[p] && !src[p]) {
120 srcContext = sws_getContext(w, h, PIX_FMT_YUVA420P, srcW, srcH,
121 srcFormat, SWS_BILINEAR, NULL, NULL, NULL);
123 fprintf(stderr, "Failed to get %s ---> %s\n",
124 av_pix_fmt_descriptors[PIX_FMT_YUVA420P].name,
125 av_pix_fmt_descriptors[srcFormat].name);
130 sws_scale(srcContext, ref, refStride, 0, h, src, srcStride);
131 sws_freeContext(srcContext);
133 cur_srcFormat = srcFormat;
138 av_image_fill_linesizes(dstStride, dstFormat, dstW);
139 for (i=0; i<4; i++) {
140 /* Image buffers passed into libswscale can be allocated any way you
141 * prefer, as long as they're aligned enough for the architecture, and
142 * they're freed appropriately (such as using av_free for buffers
143 * allocated with av_malloc). */
144 /* An extra 16 bytes is being allocated because some scalers may write
147 dst[i]= av_mallocz(dstStride[i]*dstH+16);
148 if (dstStride[i] && !dst[i]) {
156 dstContext= sws_getContext(srcW, srcH, srcFormat, dstW, dstH, dstFormat, flags, NULL, NULL, NULL);
158 fprintf(stderr, "Failed to get %s ---> %s\n",
159 av_pix_fmt_descriptors[srcFormat].name,
160 av_pix_fmt_descriptors[dstFormat].name);
165 // printf("test %X %X %X -> %X %X %X\n", (int)ref[0], (int)ref[1], (int)ref[2],
166 // (int)src[0], (int)src[1], (int)src[2]);
168 printf(" %s %dx%d -> %s %3dx%3d flags=%2d",
169 av_pix_fmt_descriptors[srcFormat].name, srcW, srcH,
170 av_pix_fmt_descriptors[dstFormat].name, dstW, dstH,
174 sws_scale(dstContext, src, srcStride, 0, srcH, dst, dstStride);
176 for (i = 0; i < 4 && dstStride[i]; i++) {
177 crc = av_crc(av_crc_get_table(AV_CRC_32_IEEE), crc, dst[i], dstStride[i] * dstH);
180 if (r && crc == r->crc) {
186 for (i=0; i<4; i++) {
188 out[i]= av_mallocz(refStride[i]*h);
189 if (refStride[i] && !out[i]) {
196 outContext= sws_getContext(dstW, dstH, dstFormat, w, h, PIX_FMT_YUVA420P, SWS_BILINEAR, NULL, NULL, NULL);
198 fprintf(stderr, "Failed to get %s ---> %s\n",
199 av_pix_fmt_descriptors[dstFormat].name,
200 av_pix_fmt_descriptors[PIX_FMT_YUVA420P].name);
205 sws_scale(outContext, dst, dstStride, 0, dstH, out, refStride);
207 ssdY= getSSD(ref[0], out[0], refStride[0], refStride[0], w, h);
208 if (hasChroma(srcFormat) && hasChroma(dstFormat)) {
209 //FIXME check that output is really gray
210 ssdU= getSSD(ref[1], out[1], refStride[1], refStride[1], (w+1)>>1, (h+1)>>1);
211 ssdV= getSSD(ref[2], out[2], refStride[2], refStride[2], (w+1)>>1, (h+1)>>1);
213 if (isALPHA(srcFormat) && isALPHA(dstFormat))
214 ssdA= getSSD(ref[3], out[3], refStride[3], refStride[3], w, h);
221 sws_freeContext(outContext);
223 for (i=0; i<4; i++) {
229 printf(" CRC=%08x SSD=%5"PRId64",%5"PRId64",%5"PRId64",%5"PRId64"\n",
230 crc, ssdY, ssdU, ssdV, ssdA);
234 sws_freeContext(dstContext);
236 for (i=0; i<4; i++) {
244 static void selfTest(uint8_t *ref[4], int refStride[4], int w, int h,
245 enum PixelFormat srcFormat_in,
246 enum PixelFormat dstFormat_in)
248 const int flags[] = { SWS_FAST_BILINEAR,
249 SWS_BILINEAR, SWS_BICUBIC,
250 SWS_X , SWS_POINT , SWS_AREA, 0 };
253 const int dstW[] = { srcW - srcW/3, srcW, srcW + srcW/3, 0 };
254 const int dstH[] = { srcH - srcH/3, srcH, srcH + srcH/3, 0 };
255 enum PixelFormat srcFormat, dstFormat;
257 for (srcFormat = srcFormat_in != PIX_FMT_NONE ? srcFormat_in : 0;
258 srcFormat < PIX_FMT_NB; srcFormat++) {
259 if (!sws_isSupportedInput(srcFormat) || !sws_isSupportedOutput(srcFormat))
262 for (dstFormat = dstFormat_in != PIX_FMT_NONE ? dstFormat_in : 0;
263 dstFormat < PIX_FMT_NB; dstFormat++) {
267 if (!sws_isSupportedInput(dstFormat) || !sws_isSupportedOutput(dstFormat))
271 av_pix_fmt_descriptors[srcFormat].name,
272 av_pix_fmt_descriptors[dstFormat].name);
275 for (k = 0; flags[k] && !res; k++) {
276 for (i = 0; dstW[i] && !res; i++)
277 for (j = 0; dstH[j] && !res; j++)
278 res = doTest(ref, refStride, w, h,
279 srcFormat, dstFormat,
280 srcW, srcH, dstW[i], dstH[j], flags[k],
283 if (dstFormat_in != PIX_FMT_NONE)
286 if (srcFormat_in != PIX_FMT_NONE)
291 static int fileTest(uint8_t *ref[4], int refStride[4], int w, int h, FILE *fp,
292 enum PixelFormat srcFormat_in,
293 enum PixelFormat dstFormat_in)
297 while (fgets(buf, sizeof(buf), fp)) {
299 enum PixelFormat srcFormat;
302 enum PixelFormat dstFormat;
308 ret = sscanf(buf, " %12s %dx%d -> %12s %dx%d flags=%d CRC=%x"
309 " SSD=%"PRId64", %"PRId64", %"PRId64", %"PRId64"\n",
310 srcStr, &srcW, &srcH, dstStr, &dstW, &dstH,
311 &flags, &r.crc, &r.ssdY, &r.ssdU, &r.ssdV, &r.ssdA);
313 srcStr[0] = dstStr[0] = 0;
314 ret = sscanf(buf, "%12s -> %12s\n", srcStr, dstStr);
317 srcFormat = av_get_pix_fmt(srcStr);
318 dstFormat = av_get_pix_fmt(dstStr);
320 if (srcFormat == PIX_FMT_NONE || dstFormat == PIX_FMT_NONE) {
321 fprintf(stderr, "malformed input file\n");
324 if ((srcFormat_in != PIX_FMT_NONE && srcFormat_in != srcFormat) ||
325 (dstFormat_in != PIX_FMT_NONE && dstFormat_in != dstFormat))
332 doTest(ref, refStride, w, h,
333 srcFormat, dstFormat,
334 srcW, srcH, dstW, dstH, flags,
344 int main(int argc, char **argv)
346 enum PixelFormat srcFormat = PIX_FMT_NONE;
347 enum PixelFormat dstFormat = PIX_FMT_NONE;
348 uint8_t *rgb_data = av_malloc (W*H*4);
349 uint8_t *rgb_src[3]= {rgb_data, NULL, NULL};
350 int rgb_stride[3]={4*W, 0, 0};
351 uint8_t *data = av_malloc (4*W*H);
352 uint8_t *src[4]= {data, data+W*H, data+W*H*2, data+W*H*3};
353 int stride[4]={W, W, W, W};
355 struct SwsContext *sws;
360 if (!rgb_data || !data)
363 sws= sws_getContext(W/12, H/12, PIX_FMT_RGB32, W, H, PIX_FMT_YUVA420P, SWS_BILINEAR, NULL, NULL, NULL);
365 av_lfg_init(&rand, 1);
367 for (y=0; y<H; y++) {
368 for (x=0; x<W*4; x++) {
369 rgb_data[ x + y*4*W]= av_lfg_get(&rand);
372 sws_scale(sws, rgb_src, rgb_stride, 0, H, src, stride);
373 sws_freeContext(sws);
376 for (i = 1; i < argc; i += 2) {
377 if (argv[i][0] != '-' || i+1 == argc)
379 if (!strcmp(argv[i], "-ref")) {
380 FILE *fp = fopen(argv[i+1], "r");
382 fprintf(stderr, "could not open '%s'\n", argv[i+1]);
385 res = fileTest(src, stride, W, H, fp, srcFormat, dstFormat);
388 } else if (!strcmp(argv[i], "-src")) {
389 srcFormat = av_get_pix_fmt(argv[i+1]);
390 if (srcFormat == PIX_FMT_NONE) {
391 fprintf(stderr, "invalid pixel format %s\n", argv[i+1]);
394 } else if (!strcmp(argv[i], "-dst")) {
395 dstFormat = av_get_pix_fmt(argv[i+1]);
396 if (dstFormat == PIX_FMT_NONE) {
397 fprintf(stderr, "invalid pixel format %s\n", argv[i+1]);
402 fprintf(stderr, "bad option or argument missing (%s)\n", argv[i]);
407 selfTest(src, stride, W, H, srcFormat, dstFormat);