1/* 2 * Generate a synthetic YUV video sequence suitable for codec testing. 3 * NOTE: No floats are used to guarantee bitexact output. 4 * 5 * Copyright (c) 2002 Fabrice Bellard 6 * 7 * This file is part of FFmpeg. 8 * 9 * FFmpeg is free software; you can redistribute it and/or 10 * modify it under the terms of the GNU Lesser General Public 11 * License as published by the Free Software Foundation; either 12 * version 2.1 of the License, or (at your option) any later version. 13 * 14 * FFmpeg is distributed in the hope that it will be useful, 15 * but WITHOUT ANY WARRANTY; without even the implied warranty of 16 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 17 * Lesser General Public License for more details. 18 * 19 * You should have received a copy of the GNU Lesser General Public 20 * License along with FFmpeg; if not, write to the Free Software 21 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA 22 */ 23 24#include <stdlib.h> 25#include <stdint.h> 26#include <stdio.h> 27 28#include "utils.c" 29 30static unsigned int myrnd(unsigned int *seed_ptr, int n) 31{ 32 unsigned int seed, val; 33 34 seed = *seed_ptr; 35 seed = (seed * 314159) + 1; 36 if (n == 256) { 37 val = seed >> 24; 38 } else { 39 val = seed % n; 40 } 41 *seed_ptr = seed; 42 return val; 43} 44 45#define NOISE_X 10 46#define NOISE_Y 30 47#define NOISE_W 26 48 49#define FRAC_BITS 8 50#define FRAC_ONE (1 << FRAC_BITS) 51 52/* cosine approximate with 1-x^2 */ 53static int int_cos(int a) 54{ 55 int v, neg; 56 a = a & (FRAC_ONE - 1); 57 if (a >= (FRAC_ONE / 2)) 58 a = FRAC_ONE - a; 59 neg = 0; 60 if (a > (FRAC_ONE / 4)) { 61 neg = -1; 62 a = (FRAC_ONE / 2) - a; 63 } 64 v = FRAC_ONE - ((a * a) >> 4); 65 v = (v ^ neg) - neg; 66 return v; 67} 68 69#define NB_OBJS 10 70 71typedef struct VObj { 72 int x, y, w, h; 73 int r, g, b; 74} VObj; 75 76static VObj objs[NB_OBJS]; 77 78static unsigned int seed = 1; 79 80static void gen_image(int num, int w, int h) 81{ 82 int r, g, b, x, y, i, dx, dy, x1, y1; 83 unsigned int seed1; 84 85 if (num == 0) { 86 for (i = 0; i < NB_OBJS; i++) { 87 objs[i].x = myrnd(&seed, w); 88 objs[i].y = myrnd(&seed, h); 89 objs[i].w = myrnd(&seed, w / 4) + 10; 90 objs[i].h = myrnd(&seed, h / 4) + 10; 91 objs[i].r = myrnd(&seed, 256); 92 objs[i].g = myrnd(&seed, 256); 93 objs[i].b = myrnd(&seed, 256); 94 } 95 } 96 97 /* first a moving background with gradients */ 98 /* test motion estimation */ 99 dx = int_cos(num * FRAC_ONE / 50) * 35; 100 dy = int_cos(num * FRAC_ONE / 50 + FRAC_ONE / 10) * 30; 101 for (y = 0; y < h; y++) { 102 for (x = 0; x < w; x++) { 103 x1 = (x << FRAC_BITS) + dx; 104 y1 = (y << FRAC_BITS) + dy; 105 r = ((y1 * 7) >> FRAC_BITS) & 0xff; 106 g = (((x1 + y1) * 9) >> FRAC_BITS) & 0xff; 107 b = ((x1 * 5) >> FRAC_BITS) & 0xff; 108 put_pixel(x, y, r, g, b); 109 } 110 } 111 112 /* then some noise with very high intensity to test saturation */ 113 seed1 = num; 114 for (y = 0; y < NOISE_W; y++) { 115 for (x = 0; x < NOISE_W; x++) { 116 r = myrnd(&seed1, 256); 117 g = myrnd(&seed1, 256); 118 b = myrnd(&seed1, 256); 119 put_pixel(x + NOISE_X, y + NOISE_Y, r, g, b); 120 } 121 } 122 123 /* then moving objects */ 124 for (i = 0; i < NB_OBJS; i++) { 125 VObj *p = &objs[i]; 126 seed1 = i; 127 for (y = 0; y < p->h; y++) { 128 for (x = 0; x < p->w; x++) { 129 r = p->r; 130 g = p->g; 131 b = p->b; 132 /* add a per object noise */ 133 r += myrnd(&seed1, 50); 134 g += myrnd(&seed1, 50); 135 b += myrnd(&seed1, 50); 136 put_pixel(x + p->x, y + p->y, r, g, b); 137 } 138 } 139 p->x += myrnd(&seed, 21) - 10; 140 p->y += myrnd(&seed, 21) - 10; 141 } 142} 143 144void print_help(const char* name) 145{ 146 printf("usage: %s file|dir [w=%i] [h=%i]\n" 147 "generate a test video stream\n", 148 name, DEFAULT_WIDTH, DEFAULT_HEIGHT); 149 exit(1); 150} 151 152int main(int argc, char **argv) 153{ 154 int w, h, i; 155 char buf[1024]; 156 int isdir = 0; 157 158 if (argc < 2 || argc > 4) { 159 print_help(argv[0]); 160 } 161 162 if (!freopen(argv[1], "wb", stdout)) 163 isdir = 1; 164 165 w = DEFAULT_WIDTH; 166 if(argc > 2) { 167 w = atoi(argv[2]); 168 if (w < 1) print_help(argv[0]); 169 } 170 h = DEFAULT_HEIGHT; 171 if(argc > 3) { 172 h = atoi(argv[3]); 173 if (h < 1) print_help(argv[0]); 174 } 175 176 rgb_tab = malloc(w * h * 3); 177 wrap = w * 3; 178 width = w; 179 height = h; 180 181 for (i = 0; i < DEFAULT_NB_PICT; i++) { 182 gen_image(i, w, h); 183 if (isdir) { 184 snprintf(buf, sizeof(buf), "%s%02d.pgm", argv[1], i); 185 pgmyuv_save(buf, w, h, rgb_tab); 186 } else { 187 pgmyuv_save(NULL, w, h, rgb_tab); 188 } 189 } 190 191 free(rgb_tab); 192 return 0; 193} 194