1 /*
2 * Generate a synthetic YUV video sequence suitable for codec testing.
3 *
4 * copyright (c) Sebastien Bechet <s.bechet@av7.net>
5 *
6 * This file is part of FFmpeg.
7 *
8 * FFmpeg is free software; you can redistribute it and/or
9 * modify it under the terms of the GNU Lesser General Public
10 * License as published by the Free Software Foundation; either
11 * version 2.1 of the License, or (at your option) any later version.
12 *
13 * FFmpeg is distributed in the hope that it will be useful,
14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
16 * Lesser General Public License for more details.
17 *
18 * You should have received a copy of the GNU Lesser General Public
19 * License along with FFmpeg; if not, write to the Free Software
20 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
21 */
22
23 #include <stdlib.h>
24 #include <stdio.h>
25 #include <inttypes.h>
26
27 #include "utils.c"
28
29 #define FIXP (1 << 16)
30 #define MY_PI 205887 // (M_PI * FIX)
31
int_pow(int64_t a,int p)32 static int64_t int_pow(int64_t a, int p)
33 {
34 int64_t v = FIXP;
35
36 for (; p; p--) {
37 v *= a;
38 v /= FIXP;
39 }
40
41 return v;
42 }
43
int_sin(int64_t a)44 static int64_t int_sin(int64_t a)
45 {
46 if (a < 0)
47 a = MY_PI - a; // 0..inf
48 a %= 2 * MY_PI; // 0..2PI
49
50 if (a >= MY_PI * 3 / 2)
51 a -= 2 * MY_PI; // -PI / 2 .. 3PI / 2
52 if (a >= MY_PI / 2)
53 a = MY_PI - a; // -PI / 2 .. PI / 2
54
55 return a - int_pow(a, 3) / 6 + int_pow(a, 5) / 120 - int_pow(a, 7) / 5040;
56 }
57
58 static unsigned char tab_r[256 * 256];
59 static unsigned char tab_g[256 * 256];
60 static unsigned char tab_b[256 * 256];
61
62 static int h_cos[360];
63 static int h_sin[360];
64
ipol(uint8_t * src,int x,int y)65 static int ipol(uint8_t *src, int x, int y)
66 {
67 int int_x = x >> 16;
68 int int_y = y >> 16;
69 int frac_x = x & 0xFFFF;
70 int frac_y = y & 0xFFFF;
71 int s00 = src[( int_x & 255) + 256 * ( int_y & 255)];
72 int s01 = src[((int_x + 1) & 255) + 256 * ( int_y & 255)];
73 int s10 = src[( int_x & 255) + 256 * ((int_y + 1) & 255)];
74 int s11 = src[((int_x + 1) & 255) + 256 * ((int_y + 1) & 255)];
75 int s0 = (((1 << 16) - frac_x) * s00 + frac_x * s01) >> 8;
76 int s1 = (((1 << 16) - frac_x) * s10 + frac_x * s11) >> 8;
77
78 return (((1 << 16) - frac_y) * s0 + frac_y * s1) >> 24;
79 }
80
gen_image(int num,int w,int h)81 static void gen_image(int num, int w, int h)
82 {
83 const int c = h_cos[num % 360];
84 const int s = h_sin[num % 360];
85
86 const int xi = -(w / 2) * c;
87 const int yi = (w / 2) * s;
88
89 const int xj = -(h / 2) * s;
90 const int yj = -(h / 2) * c;
91 int i, j;
92
93 int x, y;
94 int xprime = xj;
95 int yprime = yj;
96
97 for (j = 0; j < h; j++) {
98 x = xprime + xi + FIXP * w / 2;
99 xprime += s;
100
101 y = yprime + yi + FIXP * h / 2;
102 yprime += c;
103
104 for (i = 0; i < w; i++) {
105 x += c;
106 y -= s;
107 put_pixel(i, j,
108 ipol(tab_r, x, y),
109 ipol(tab_g, x, y),
110 ipol(tab_b, x, y));
111 }
112 }
113 }
114
115 #define W 256
116 #define H 256
117
init_demo(const char * filename)118 static int init_demo(const char *filename)
119 {
120 int i, j;
121 int h;
122 int radian;
123 char line[3 * W];
124
125 FILE *input_file;
126
127 input_file = fopen(filename, "rb");
128 if (!input_file) {
129 perror(filename);
130 return 1;
131 }
132
133 if (fread(line, 1, 15, input_file) != 15)
134 return 1;
135 for (i = 0; i < H; i++) {
136 if (fread(line, 1, 3 * W, input_file) != 3 * W)
137 return 1;
138 for (j = 0; j < W; j++) {
139 tab_r[W * i + j] = line[3 * j ];
140 tab_g[W * i + j] = line[3 * j + 1];
141 tab_b[W * i + j] = line[3 * j + 2];
142 }
143 }
144 fclose(input_file);
145
146 /* tables sin/cos */
147 for (i = 0; i < 360; i++) {
148 radian = 2 * i * MY_PI / 360;
149 h = 2 * FIXP + int_sin(radian);
150 h_cos[i] = h * int_sin(radian + MY_PI / 2) / 2 / FIXP;
151 h_sin[i] = h * int_sin(radian) / 2 / FIXP;
152 }
153
154 return 0;
155 }
156
main(int argc,char ** argv)157 int main(int argc, char **argv)
158 {
159 int w, h, i;
160 char buf[1024];
161 int isdir = 0;
162
163 if (argc != 3) {
164 printf("usage: %s image.pnm file|dir\n"
165 "generate a test video stream\n", argv[0]);
166 return 1;
167 }
168
169 if (!freopen(argv[2], "wb", stdout))
170 isdir = 1;
171
172 w = DEFAULT_WIDTH;
173 h = DEFAULT_HEIGHT;
174
175 rgb_tab = malloc(w * h * 3);
176 wrap = w * 3;
177 width = w;
178 height = h;
179
180 if (init_demo(argv[1]))
181 return 1;
182
183 for (i = 0; i < DEFAULT_NB_PICT; i++) {
184 gen_image(i, w, h);
185 if (isdir) {
186 snprintf(buf, sizeof(buf), "%s%02d.pgm", argv[2], i);
187 pgmyuv_save(buf, w, h, rgb_tab);
188 } else {
189 pgmyuv_save(NULL, w, h, rgb_tab);
190 }
191 }
192
193 free(rgb_tab);
194 return 0;
195 }
196