Mercurial > hg > forks > libbpg
annotate bpgenc.c @ 4:bbd61622856a
Initial implementation of TIFF format read support via libtiff.
author | Matti Hamalainen <ccr@tnsp.org> |
---|---|
date | Wed, 07 Dec 2016 10:37:38 +0200 |
parents | ffcd1967fb5c |
children | 524eae707ba4 |
rev | line source |
---|---|
0 | 1 /* |
2 * BPG encoder | |
3 * | |
4 * Copyright (c) 2014 Fabrice Bellard | |
5 * | |
6 * Permission is hereby granted, free of charge, to any person obtaining a copy | |
7 * of this software and associated documentation files (the "Software"), to deal | |
8 * in the Software without restriction, including without limitation the rights | |
9 * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell | |
10 * copies of the Software, and to permit persons to whom the Software is | |
11 * furnished to do so, subject to the following conditions: | |
12 * | |
13 * The above copyright notice and this permission notice shall be included in | |
14 * all copies or substantial portions of the Software. | |
15 * | |
16 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR | |
17 * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, | |
18 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL | |
19 * THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER | |
20 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, | |
21 * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN | |
22 * THE SOFTWARE. | |
23 */ | |
24 #include <stdlib.h> | |
25 #include <stdio.h> | |
26 #include <string.h> | |
27 #include <inttypes.h> | |
28 #include <getopt.h> | |
29 #include <math.h> | |
30 #include <assert.h> | |
31 | |
32 #include <png.h> | |
33 #include <jpeglib.h> | |
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
34 #include <tiffio.h> |
0 | 35 |
36 #include "bpgenc.h" | |
37 | |
38 typedef uint16_t PIXEL; | |
39 | |
40 static void put_ue(uint8_t **pp, uint32_t v); | |
41 | |
42 static inline int clamp_pix(int a, int pixel_max) | |
43 { | |
44 if (a < 0) | |
45 return 0; | |
46 else if (a > pixel_max) | |
47 return pixel_max; | |
48 else | |
49 return a; | |
50 } | |
51 | |
52 static inline int sub_mod_int(int a, int b, int m) | |
53 { | |
54 a -= b; | |
55 if (a < 0) | |
56 a += m; | |
57 return a; | |
58 } | |
59 | |
60 static inline int add_mod_int(int a, int b, int m) | |
61 { | |
62 a += b; | |
63 if (a >= m) | |
64 a -= m; | |
65 return a; | |
66 } | |
67 | |
68 typedef struct { | |
69 int c_shift; | |
70 int c_rnd; | |
71 int c_0_25, c_0_5, c_one; | |
72 int rgb_to_ycc[3 * 3]; | |
73 int y_one; | |
74 int y_offset; | |
75 int bit_depth; | |
76 int pixel_max; | |
77 int c_center; | |
78 } ColorConvertState; | |
79 | |
80 static void convert_init(ColorConvertState *s, int in_bit_depth, | |
81 int out_bit_depth, BPGColorSpaceEnum color_space, | |
82 int limited_range) | |
83 { | |
84 double k_r, k_b, mult, mult_y, mult_c; | |
85 int in_pixel_max, out_pixel_max, c_shift, i; | |
86 double rgb_to_ycc[3 * 3]; | |
87 | |
88 /* XXX: could use one more bit */ | |
89 c_shift = 31 - out_bit_depth; | |
90 in_pixel_max = (1 << in_bit_depth) - 1; | |
91 out_pixel_max = (1 << out_bit_depth) - 1; | |
92 mult = (double)out_pixel_max * (1 << c_shift) / (double)in_pixel_max; | |
93 // printf("mult=%f c_shift=%d\n", mult, c_shift); | |
94 if (limited_range) { | |
95 mult_y = (double)(219 << (out_bit_depth - 8)) * (1 << c_shift) / | |
96 (double)in_pixel_max; | |
97 mult_c = (double)(224 << (out_bit_depth - 8)) * (1 << c_shift) / | |
98 (double)in_pixel_max; | |
99 } else { | |
100 mult_y = mult; | |
101 mult_c = mult; | |
102 } | |
103 switch(color_space) { | |
104 case BPG_CS_YCbCr: | |
105 k_r = 0.299; | |
106 k_b = 0.114; | |
107 goto convert_ycc; | |
108 | |
109 case BPG_CS_YCbCr_BT709: | |
110 k_r = 0.2126; | |
111 k_b = 0.0722; | |
112 goto convert_ycc; | |
113 | |
114 case BPG_CS_YCbCr_BT2020: | |
115 k_r = 0.2627; | |
116 k_b = 0.0593; | |
117 convert_ycc: | |
118 rgb_to_ycc[0] = k_r; | |
119 rgb_to_ycc[1] = 1 - k_r - k_b; | |
120 rgb_to_ycc[2] = k_b; | |
121 rgb_to_ycc[3] = -0.5 * k_r / (1 - k_b); | |
122 rgb_to_ycc[4] = -0.5 * (1 - k_r - k_b) / (1 - k_b); | |
123 rgb_to_ycc[5] = 0.5; | |
124 rgb_to_ycc[6] = 0.5; | |
125 rgb_to_ycc[7] = -0.5 * (1 - k_r - k_b) / (1 - k_r); | |
126 rgb_to_ycc[8] = -0.5 * k_b / (1 - k_r); | |
127 | |
128 for(i = 0; i < 3; i++) | |
129 s->rgb_to_ycc[i] = lrint(rgb_to_ycc[i] * mult_y); | |
130 for(i = 3; i < 9; i++) | |
131 s->rgb_to_ycc[i] = lrint(rgb_to_ycc[i] * mult_c); | |
132 break; | |
133 case BPG_CS_YCgCo: | |
134 s->c_0_25 = lrint(0.25 * mult_y); | |
135 s->c_0_5 = lrint(0.5 * mult_y); | |
136 break; | |
137 default: | |
138 break; | |
139 } | |
140 | |
141 s->c_one = lrint(mult); | |
142 s->c_shift = c_shift; | |
143 s->c_rnd = (1 << (c_shift - 1)); | |
144 if (limited_range) { | |
145 s->y_offset = s->c_rnd + (16 << (c_shift + out_bit_depth - 8)); | |
146 s->y_one = lrint(mult_y); | |
147 } else { | |
148 s->y_offset = s->c_rnd; | |
149 s->y_one = s->c_one; | |
150 } | |
151 | |
152 s->bit_depth = out_bit_depth; | |
153 s->c_center = 1 << (out_bit_depth - 1); | |
154 s->pixel_max = out_pixel_max; | |
155 } | |
156 | |
157 /* 8 bit input */ | |
158 static void rgb24_to_ycc(ColorConvertState *s, | |
159 PIXEL *y_ptr, PIXEL *cb_ptr, PIXEL *cr_ptr, | |
160 const void *src1, int n, int incr) | |
161 { | |
162 const uint8_t *src = src1; | |
163 int i, r, g, b, c0, c1, c2, c3, c4, c5, c6, c7, c8, shift, rnd, center; | |
164 int pixel_max, y_offset; | |
165 | |
166 c0 = s->rgb_to_ycc[0]; | |
167 c1 = s->rgb_to_ycc[1]; | |
168 c2 = s->rgb_to_ycc[2]; | |
169 c3 = s->rgb_to_ycc[3]; | |
170 c4 = s->rgb_to_ycc[4]; | |
171 c5 = s->rgb_to_ycc[5]; | |
172 c6 = s->rgb_to_ycc[6]; | |
173 c7 = s->rgb_to_ycc[7]; | |
174 c8 = s->rgb_to_ycc[8]; | |
175 shift = s->c_shift; | |
176 rnd = s->c_rnd; | |
177 y_offset = s->y_offset; | |
178 center = s->c_center; | |
179 pixel_max = s->pixel_max; | |
180 for(i = 0; i < n; i++) { | |
181 r = src[0]; | |
182 g = src[1]; | |
183 b = src[2]; | |
184 y_ptr[i] = clamp_pix((c0 * r + c1 * g + c2 * b + | |
185 y_offset) >> shift, pixel_max); | |
186 cb_ptr[i] = clamp_pix(((c3 * r + c4 * g + c5 * b + | |
187 rnd) >> shift) + center, pixel_max); | |
188 cr_ptr[i] = clamp_pix(((c6 * r + c7 * g + c8 * b + | |
189 rnd) >> shift) + center, pixel_max); | |
190 src += incr; | |
191 } | |
192 } | |
193 | |
194 static void rgb24_to_rgb(ColorConvertState *s, | |
195 PIXEL *y_ptr, PIXEL *cb_ptr, PIXEL *cr_ptr, | |
196 const void *src1, int n, int incr) | |
197 { | |
198 const uint8_t *src = src1; | |
199 int i, r, g, b, c, shift, rnd; | |
200 | |
201 c = s->y_one; | |
202 shift = s->c_shift; | |
203 rnd = s->y_offset; | |
204 for(i = 0; i < n; i++) { | |
205 r = src[0]; | |
206 g = src[1]; | |
207 b = src[2]; | |
208 y_ptr[i] = (c * g + rnd) >> shift; | |
209 cb_ptr[i] = (c * b + rnd) >> shift; | |
210 cr_ptr[i] = (c * r + rnd) >> shift; | |
211 src += incr; | |
212 } | |
213 } | |
214 | |
215 static void rgb24_to_ycgco(ColorConvertState *s, | |
216 PIXEL *y_ptr, PIXEL *cb_ptr, PIXEL *cr_ptr, | |
217 const void *src1, int n, int incr) | |
218 { | |
219 const uint8_t *src = src1; | |
220 int i, r, g, b, t1, t2, pixel_max, c_0_5, c_0_25, rnd, shift, center; | |
221 int y_offset; | |
222 | |
223 c_0_25 = s->c_0_25; | |
224 c_0_5 = s->c_0_5; | |
225 rnd = s->c_rnd; | |
226 shift = s->c_shift; | |
227 pixel_max = s->pixel_max; | |
228 center = s->c_center; | |
229 y_offset = s->y_offset; | |
230 for(i = 0; i < n; i++) { | |
231 r = src[0]; | |
232 g = src[1]; | |
233 b = src[2]; | |
234 t1 = c_0_5 * g; | |
235 t2 = c_0_25 * (r + b); | |
236 y_ptr[i] = clamp_pix((t1 + t2 + y_offset) >> shift, pixel_max); | |
237 cb_ptr[i] = clamp_pix(((t1 - t2 + rnd) >> shift) + center, | |
238 pixel_max); | |
239 cr_ptr[i] = clamp_pix(((c_0_5 * (r - b) + | |
240 rnd) >> shift) + center, pixel_max); | |
241 src += incr; | |
242 } | |
243 } | |
244 | |
245 /* Note: used for alpha/W so no limited range */ | |
246 static void gray8_to_gray(ColorConvertState *s, | |
247 PIXEL *y_ptr, const uint8_t *src, int n, int incr) | |
248 { | |
249 int i, g, c, shift, rnd; | |
250 | |
251 c = s->c_one; | |
252 shift = s->c_shift; | |
253 rnd = s->c_rnd; | |
254 for(i = 0; i < n; i++) { | |
255 g = src[0]; | |
256 y_ptr[i] = (c * g + rnd) >> shift; | |
257 src += incr; | |
258 } | |
259 } | |
260 | |
261 static void luma8_to_gray(ColorConvertState *s, | |
262 PIXEL *y_ptr, const uint8_t *src, int n, int incr) | |
263 { | |
264 int i, g, c, shift, rnd; | |
265 | |
266 c = s->y_one; | |
267 shift = s->c_shift; | |
268 rnd = s->y_offset; | |
269 for(i = 0; i < n; i++) { | |
270 g = src[0]; | |
271 y_ptr[i] = (c * g + rnd) >> shift; | |
272 src += incr; | |
273 } | |
274 } | |
275 | |
276 /* 16 bit input */ | |
277 | |
278 static void rgb48_to_ycc(ColorConvertState *s, | |
279 PIXEL *y_ptr, PIXEL *cb_ptr, PIXEL *cr_ptr, | |
280 const void *src1, int n, int incr) | |
281 { | |
282 const uint16_t *src = src1; | |
283 int i, r, g, b, c0, c1, c2, c3, c4, c5, c6, c7, c8, shift, rnd, center; | |
284 int pixel_max, y_offset; | |
285 | |
286 c0 = s->rgb_to_ycc[0]; | |
287 c1 = s->rgb_to_ycc[1]; | |
288 c2 = s->rgb_to_ycc[2]; | |
289 c3 = s->rgb_to_ycc[3]; | |
290 c4 = s->rgb_to_ycc[4]; | |
291 c5 = s->rgb_to_ycc[5]; | |
292 c6 = s->rgb_to_ycc[6]; | |
293 c7 = s->rgb_to_ycc[7]; | |
294 c8 = s->rgb_to_ycc[8]; | |
295 shift = s->c_shift; | |
296 rnd = s->c_rnd; | |
297 y_offset = s->y_offset; | |
298 center = s->c_center; | |
299 pixel_max = s->pixel_max; | |
300 for(i = 0; i < n; i++) { | |
301 r = src[0]; | |
302 g = src[1]; | |
303 b = src[2]; | |
304 y_ptr[i] = clamp_pix((c0 * r + c1 * g + c2 * b + | |
305 y_offset) >> shift, pixel_max); | |
306 cb_ptr[i] = clamp_pix(((c3 * r + c4 * g + c5 * b + | |
307 rnd) >> shift) + center, pixel_max); | |
308 cr_ptr[i] = clamp_pix(((c6 * r + c7 * g + c8 * b + | |
309 rnd) >> shift) + center, pixel_max); | |
310 src += incr; | |
311 } | |
312 } | |
313 | |
314 static void rgb48_to_ycgco(ColorConvertState *s, | |
315 PIXEL *y_ptr, PIXEL *cb_ptr, PIXEL *cr_ptr, | |
316 const void *src1, int n, int incr) | |
317 { | |
318 const uint16_t *src = src1; | |
319 int i, r, g, b, t1, t2, pixel_max, c_0_5, c_0_25, rnd, shift, center; | |
320 int y_offset; | |
321 | |
322 c_0_25 = s->c_0_25; | |
323 c_0_5 = s->c_0_5; | |
324 rnd = s->c_rnd; | |
325 y_offset = s->y_offset; | |
326 shift = s->c_shift; | |
327 pixel_max = s->pixel_max; | |
328 center = s->c_center; | |
329 for(i = 0; i < n; i++) { | |
330 r = src[0]; | |
331 g = src[1]; | |
332 b = src[2]; | |
333 t1 = c_0_5 * g; | |
334 t2 = c_0_25 * (r + b); | |
335 y_ptr[i] = clamp_pix((t1 + t2 + y_offset) >> shift, pixel_max); | |
336 cb_ptr[i] = clamp_pix(((t1 - t2 + rnd) >> shift) + center, | |
337 pixel_max); | |
338 cr_ptr[i] = clamp_pix(((c_0_5 * (r - b) + | |
339 rnd) >> shift) + center, pixel_max); | |
340 src += incr; | |
341 } | |
342 } | |
343 | |
344 /* Note: use for alpha/W so no limited range */ | |
345 static void gray16_to_gray(ColorConvertState *s, | |
346 PIXEL *y_ptr, const uint16_t *src, int n, int incr) | |
347 { | |
348 int i, g, c, shift, rnd; | |
349 | |
350 c = s->c_one; | |
351 shift = s->c_shift; | |
352 rnd = s->c_rnd; | |
353 for(i = 0; i < n; i++) { | |
354 g = src[0]; | |
355 y_ptr[i] = (c * g + rnd) >> shift; | |
356 src += incr; | |
357 } | |
358 } | |
359 | |
360 static void luma16_to_gray(ColorConvertState *s, | |
361 PIXEL *y_ptr, const uint16_t *src, int n, int incr) | |
362 { | |
363 int i, g, c, shift, rnd; | |
364 | |
365 c = s->y_one; | |
366 shift = s->c_shift; | |
367 rnd = s->y_offset; | |
368 for(i = 0; i < n; i++) { | |
369 g = src[0]; | |
370 y_ptr[i] = (c * g + rnd) >> shift; | |
371 src += incr; | |
372 } | |
373 } | |
374 | |
375 static void rgb48_to_rgb(ColorConvertState *s, | |
376 PIXEL *y_ptr, PIXEL *cb_ptr, PIXEL *cr_ptr, | |
377 const void *src1, int n, int incr) | |
378 { | |
379 const uint16_t *src = src1; | |
380 | |
381 luma16_to_gray(s, y_ptr, src + 1, n, incr); | |
382 luma16_to_gray(s, cb_ptr, src + 2, n, incr); | |
383 luma16_to_gray(s, cr_ptr, src + 0, n, incr); | |
384 } | |
385 | |
386 typedef void RGBConvertFunc(ColorConvertState *s, | |
387 PIXEL *y_ptr, PIXEL *cb_ptr, PIXEL *cr_ptr, | |
388 const void *src, int n, int incr); | |
389 | |
390 static RGBConvertFunc *rgb_to_cs[2][BPG_CS_COUNT] = { | |
391 { | |
392 rgb24_to_ycc, | |
393 rgb24_to_rgb, | |
394 rgb24_to_ycgco, | |
395 rgb24_to_ycc, | |
396 rgb24_to_ycc, | |
397 }, | |
398 { | |
399 rgb48_to_ycc, | |
400 rgb48_to_rgb, | |
401 rgb48_to_ycgco, | |
402 rgb48_to_ycc, | |
403 rgb48_to_ycc, | |
404 } | |
405 }; | |
406 | |
407 /* val = 1.0 - val */ | |
408 static void gray_one_minus(ColorConvertState *s, PIXEL *y_ptr, int n) | |
409 { | |
410 int pixel_max = s->pixel_max; | |
411 int i; | |
412 | |
413 for(i = 0; i < n; i++) { | |
414 y_ptr[i] = pixel_max - y_ptr[i]; | |
415 } | |
416 } | |
417 | |
418 /* val = -val for chroma */ | |
419 static void gray_neg_c(ColorConvertState *s, PIXEL *y_ptr, int n) | |
420 { | |
421 int pixel_max = s->pixel_max; | |
422 int i, v; | |
423 | |
424 for(i = 0; i < n; i++) { | |
425 v = y_ptr[i]; | |
426 if (v == 0) | |
427 v = pixel_max; | |
428 else | |
429 v = pixel_max + 1 - v; | |
430 y_ptr[i] = v; | |
431 } | |
432 } | |
433 | |
434 | |
435 /* decimation */ | |
436 | |
437 /* phase = 0 */ | |
438 #define DP0TAPS2 7 | |
439 #define DP0TAPS (2 * DP0TAPS + 1) | |
440 #define DP0C0 64 | |
441 #define DP0C1 40 | |
442 #define DP0C3 (-11) | |
443 #define DP0C5 4 | |
444 #define DP0C7 (-1) | |
445 | |
446 /* phase = 0.5 */ | |
447 #define DP1TAPS2 5 | |
448 #define DP1TAPS (2 * DP1TAPS2) | |
449 #define DP1C0 57 | |
450 #define DP1C1 17 | |
451 #define DP1C2 (-8) | |
452 #define DP1C3 (-4) | |
453 #define DP1C4 2 | |
454 | |
455 #define DTAPS_MAX 7 | |
456 | |
457 /* chroma aligned with luma samples */ | |
458 static void decimate2p0_simple(PIXEL *dst, PIXEL *src, int n, int bit_depth) | |
459 { | |
460 int n2, i, pixel_max; | |
461 pixel_max = (1 << bit_depth) - 1; | |
462 n2 = (n + 1) / 2; | |
463 for(i = 0; i < n2; i++) { | |
464 dst[i] = clamp_pix(((src[-7] + src[7]) * DP0C7 + | |
465 (src[-5] + src[5]) * DP0C5 + | |
466 (src[-3] + src[3]) * DP0C3 + | |
467 (src[-1] + src[1]) * DP0C1 + | |
468 src[0] * DP0C0 + 64) >> 7, pixel_max); | |
469 src += 2; | |
470 } | |
471 } | |
472 | |
473 /* same with more precision and no saturation */ | |
474 static void decimate2p0_simple16(int16_t *dst, PIXEL *src, int n, int bit_depth) | |
475 { | |
476 int n2, i, shift, rnd; | |
477 shift = bit_depth - 7; | |
478 rnd = 1 << (shift - 1); | |
479 n2 = (n + 1) / 2; | |
480 for(i = 0; i < n2; i++) { | |
481 dst[i] = ((src[-7] + src[7]) * DP0C7 + | |
482 (src[-5] + src[5]) * DP0C5 + | |
483 (src[-3] + src[3]) * DP0C3 + | |
484 (src[-1] + src[1]) * DP0C1 + | |
485 src[0] * DP0C0 + rnd) >> shift; | |
486 src += 2; | |
487 } | |
488 } | |
489 | |
490 | |
491 /* chroma half way between luma samples */ | |
492 static void decimate2p1_simple(PIXEL *dst, PIXEL *src, int n, int bit_depth) | |
493 { | |
494 int n2, i, pixel_max; | |
495 pixel_max = (1 << bit_depth) - 1; | |
496 n2 = (n + 1) / 2; | |
497 for(i = 0; i < n2; i++) { | |
498 dst[i] = clamp_pix(((src[-4] + src[5]) * DP1C4 + | |
499 (src[-3] + src[4]) * DP1C3 + | |
500 (src[-2] + src[3]) * DP1C2 + | |
501 (src[-1] + src[2]) * DP1C1 + | |
502 (src[0] + src[1]) * DP1C0 + 64) >> 7, pixel_max); | |
503 src += 2; | |
504 } | |
505 } | |
506 | |
507 /* same with more precision and no saturation */ | |
508 static void decimate2p1_simple16(int16_t *dst, PIXEL *src, int n, int bit_depth) | |
509 { | |
510 int n2, i, shift, rnd; | |
511 shift = bit_depth - 7; | |
512 rnd = 1 << (shift - 1); | |
513 n2 = (n + 1) / 2; | |
514 for(i = 0; i < n2; i++) { | |
515 dst[i] = ((src[-4] + src[5]) * DP1C4 + | |
516 (src[-3] + src[4]) * DP1C3 + | |
517 (src[-2] + src[3]) * DP1C2 + | |
518 (src[-1] + src[2]) * DP1C1 + | |
519 (src[0] + src[1]) * DP1C0 + rnd) >> shift; | |
520 src += 2; | |
521 } | |
522 } | |
523 | |
524 static void decimate2_h(PIXEL *dst, PIXEL *src, int n, int bit_depth, int phase) | |
525 { | |
526 PIXEL *src1, v; | |
527 int d, i; | |
528 | |
529 if (phase == 0) | |
530 d = DP0TAPS2; | |
531 else | |
532 d = DP1TAPS2; | |
533 /* add edge pixels */ | |
534 src1 = malloc(sizeof(PIXEL) * (n + 2 * d)); | |
535 v = src[0]; | |
536 for(i = 0; i < d; i++) | |
537 src1[i] = v; | |
538 memcpy(src1 + d, src, n * sizeof(PIXEL)); | |
539 v = src[n - 1]; | |
540 for(i = 0; i < d; i++) | |
541 src1[d + n + i] = v; | |
542 if (phase == 0) | |
543 decimate2p0_simple(dst, src1 + d, n, bit_depth); | |
544 else | |
545 decimate2p1_simple(dst, src1 + d, n, bit_depth); | |
546 free(src1); | |
547 } | |
548 | |
549 /* src1 is a temporary buffer of length n + 2 * DTAPS */ | |
550 static void decimate2_h16(int16_t *dst, PIXEL *src, int n, PIXEL *src1, | |
551 int bit_depth, int phase) | |
552 { | |
553 PIXEL v; | |
554 int d, i; | |
555 | |
556 if (phase == 0) | |
557 d = DP0TAPS2; | |
558 else | |
559 d = DP1TAPS2; | |
560 /* add edge pixels */ | |
561 v = src[0]; | |
562 for(i = 0; i < d; i++) | |
563 src1[i] = v; | |
564 memcpy(src1 + d, src, n * sizeof(PIXEL)); | |
565 v = src[n - 1]; | |
566 for(i = 0; i < d; i++) | |
567 src1[d + n + i] = v; | |
568 if (phase == 0) | |
569 decimate2p0_simple16(dst, src1 + d, n, bit_depth); | |
570 else | |
571 decimate2p1_simple16(dst, src1 + d, n, bit_depth); | |
572 | |
573 } | |
574 | |
575 static void decimate2_v(PIXEL *dst, int16_t **src, int pos, int n, | |
576 int bit_depth) | |
577 { | |
578 int16_t *src0, *src1, *src2, *src3, *src4, *src5, *srcm1, *srcm2, *srcm3, *srcm4; | |
579 int i, shift, offset, pixel_max; | |
580 | |
581 pos = sub_mod_int(pos, 4, DP1TAPS); | |
582 srcm4 = src[pos]; | |
583 pos = add_mod_int(pos, 1, DP1TAPS); | |
584 srcm3 = src[pos]; | |
585 pos = add_mod_int(pos, 1, DP1TAPS); | |
586 srcm2 = src[pos]; | |
587 pos = add_mod_int(pos, 1, DP1TAPS); | |
588 srcm1 = src[pos]; | |
589 pos = add_mod_int(pos, 1, DP1TAPS); | |
590 src0 = src[pos]; | |
591 pos = add_mod_int(pos, 1, DP1TAPS); | |
592 src1 = src[pos]; | |
593 pos = add_mod_int(pos, 1, DP1TAPS); | |
594 src2 = src[pos]; | |
595 pos = add_mod_int(pos, 1, DP1TAPS); | |
596 src3 = src[pos]; | |
597 pos = add_mod_int(pos, 1, DP1TAPS); | |
598 src4 = src[pos]; | |
599 pos = add_mod_int(pos, 1, DP1TAPS); | |
600 src5 = src[pos]; | |
601 | |
602 shift = 21 - bit_depth; | |
603 offset = 1 << (shift - 1); | |
604 pixel_max = (1 << bit_depth) - 1; | |
605 for(i = 0; i < n; i++) { | |
606 dst[i] = clamp_pix(((srcm4[i] + src5[i]) * DP1C4 + | |
607 (srcm3[i] + src4[i]) * DP1C3 + | |
608 (srcm2[i] + src3[i]) * DP1C2 + | |
609 (srcm1[i] + src2[i]) * DP1C1 + | |
610 (src0[i] + src1[i]) * DP1C0 + offset) >> shift, pixel_max); | |
611 } | |
612 } | |
613 | |
614 /* Note: we do the horizontal decimation first to use less CPU cache */ | |
615 static void decimate2_hv(uint8_t *dst, int dst_linesize, | |
616 uint8_t *src, int src_linesize, | |
617 int w, int h, int bit_depth, int h_phase) | |
618 { | |
619 PIXEL *buf1; | |
620 int16_t *buf2[DP1TAPS]; | |
621 int w2, pos, i, y, y1, y2; | |
622 | |
623 w2 = (w + 1) / 2; | |
624 | |
625 buf1 = malloc(sizeof(PIXEL) * (w + 2 * DTAPS_MAX)); | |
626 /* init line buffer */ | |
627 for(i = 0; i < DP1TAPS; i++) { | |
628 buf2[i] = malloc(sizeof(int16_t) * w2); | |
629 y = i; | |
630 if (y > DP1TAPS2) | |
631 y -= DP1TAPS; | |
632 if (y < 0) { | |
633 /* copy from first line */ | |
634 memcpy(buf2[i], buf2[0], sizeof(int16_t) * w2); | |
635 } else if (y >= h) { | |
636 /* copy from last line (only happens for small height) */ | |
637 memcpy(buf2[i], buf2[h - 1], sizeof(int16_t) * w2); | |
638 } else { | |
639 decimate2_h16(buf2[i], (PIXEL *)(src + src_linesize * y), w, | |
640 buf1, bit_depth, h_phase); | |
641 } | |
642 } | |
643 | |
644 for(y = 0; y < h; y++) { | |
645 pos = y % DP1TAPS; | |
646 if ((y & 1) == 0) { | |
647 /* filter one line */ | |
648 y2 = y >> 1; | |
649 decimate2_v((PIXEL *)(dst + y2 * dst_linesize), buf2, | |
650 pos, w2, bit_depth); | |
651 } | |
652 /* add a new line in the buffer */ | |
653 y1 = y + DP1TAPS2 + 1; | |
654 pos = add_mod_int(pos, DP1TAPS2 + 1, DP1TAPS); | |
655 if (y1 >= h) { | |
656 /* copy last line */ | |
657 memcpy(buf2[pos], buf2[sub_mod_int(pos, 1, DP1TAPS)], | |
658 sizeof(int16_t) * w2); | |
659 } else { | |
660 /* horizontally decimate new line */ | |
661 decimate2_h16(buf2[pos], (PIXEL *)(src + src_linesize * y1), w, | |
662 buf1, bit_depth, h_phase); | |
663 } | |
664 } | |
665 | |
666 for(i = 0; i < DP1TAPS; i++) | |
667 free(buf2[i]); | |
668 free(buf1); | |
669 } | |
670 | |
671 static void get_plane_res(Image *img, int *pw, int *ph, int i) | |
672 { | |
673 if (img->format == BPG_FORMAT_420 && (i == 1 || i == 2)) { | |
674 *pw = (img->w + 1) / 2; | |
675 *ph = (img->h + 1) / 2; | |
676 } else if (img->format == BPG_FORMAT_422 && (i == 1 || i == 2)) { | |
677 *pw = (img->w + 1) / 2; | |
678 *ph = img->h; | |
679 } else { | |
680 *pw = img->w; | |
681 *ph = img->h; | |
682 } | |
683 } | |
684 | |
685 #define W_PAD 16 | |
686 | |
687 Image *image_alloc(int w, int h, BPGImageFormatEnum format, int has_alpha, | |
688 BPGColorSpaceEnum color_space, int bit_depth) | |
689 { | |
690 Image *img; | |
691 int i, linesize, w1, h1, c_count; | |
692 | |
693 img = malloc(sizeof(Image)); | |
694 memset(img, 0, sizeof(*img)); | |
695 | |
696 img->w = w; | |
697 img->h = h; | |
698 img->format = format; | |
699 img->has_alpha = has_alpha; | |
700 img->bit_depth = bit_depth; | |
701 img->color_space = color_space; | |
702 img->pixel_shift = 1; | |
703 img->c_h_phase = 1; | |
704 | |
705 if (img->format == BPG_FORMAT_GRAY) | |
706 c_count = 1; | |
707 else | |
708 c_count = 3; | |
709 if (has_alpha) | |
710 c_count++; | |
711 for(i = 0; i < c_count; i++) { | |
712 get_plane_res(img, &w1, &h1, i); | |
713 /* multiple of 16 pixels to add borders */ | |
714 w1 = (w1 + (W_PAD - 1)) & ~(W_PAD - 1); | |
715 h1 = (h1 + (W_PAD - 1)) & ~(W_PAD - 1); | |
716 | |
717 linesize = w1 << img->pixel_shift; | |
718 img->data[i] = malloc(linesize * h1); | |
719 img->linesize[i] = linesize; | |
720 } | |
721 return img; | |
722 } | |
723 | |
724 void image_free(Image *img) | |
725 { | |
726 int i, c_count; | |
727 if (img->format == BPG_FORMAT_GRAY) | |
728 c_count = 1; | |
729 else | |
730 c_count = 3; | |
731 if (img->has_alpha) | |
732 c_count++; | |
733 for(i = 0; i < c_count; i++) | |
734 free(img->data[i]); | |
735 free(img); | |
736 } | |
737 | |
738 int image_ycc444_to_ycc422(Image *img, int h_phase) | |
739 { | |
740 uint8_t *data1; | |
741 int w1, h1, bpp, linesize1, i, y; | |
742 | |
743 if (img->format != BPG_FORMAT_444 || img->pixel_shift != 1) | |
744 return -1; | |
745 bpp = 2; | |
746 w1 = (img->w + 1) / 2; | |
747 w1 = (w1 + (W_PAD - 1)) & ~(W_PAD - 1); | |
748 h1 = (img->h + (W_PAD - 1)) & ~(W_PAD - 1); | |
749 linesize1 = bpp * w1; | |
750 for(i = 1; i <= 2; i++) { | |
751 data1 = malloc(linesize1 * h1); | |
752 for(y = 0; y < img->h; y++) { | |
753 decimate2_h((PIXEL *)(data1 + y * linesize1), | |
754 (PIXEL *)(img->data[i] + y * img->linesize[i]), | |
755 img->w, img->bit_depth, h_phase); | |
756 } | |
757 free(img->data[i]); | |
758 img->data[i] = data1; | |
759 img->linesize[i] = linesize1; | |
760 } | |
761 img->format = BPG_FORMAT_422; | |
762 img->c_h_phase = h_phase; | |
763 return 0; | |
764 } | |
765 | |
766 int image_ycc444_to_ycc420(Image *img, int h_phase) | |
767 { | |
768 uint8_t *data1; | |
769 int w1, h1, bpp, linesize1, i; | |
770 | |
771 if (img->format != BPG_FORMAT_444 || img->pixel_shift != 1) | |
772 return -1; | |
773 bpp = 2; | |
774 w1 = (img->w + 1) / 2; | |
775 h1 = (img->h + 1) / 2; | |
776 w1 = (w1 + (W_PAD - 1)) & ~(W_PAD - 1); | |
777 h1 = (h1 + (W_PAD - 1)) & ~(W_PAD - 1); | |
778 linesize1 = bpp * w1; | |
779 for(i = 1; i <= 2; i++) { | |
780 data1 = malloc(linesize1 * h1); | |
781 decimate2_hv(data1, linesize1, | |
782 img->data[i], img->linesize[i], | |
783 img->w, img->h, img->bit_depth, h_phase); | |
784 free(img->data[i]); | |
785 img->data[i] = data1; | |
786 img->linesize[i] = linesize1; | |
787 } | |
788 img->format = BPG_FORMAT_420; | |
789 img->c_h_phase = h_phase; | |
790 return 0; | |
791 } | |
792 | |
793 /* duplicate right and bottom samples so that the image has a width | |
794 and height multiple of cb_size (power of two) */ | |
795 void image_pad(Image *img, int cb_size) | |
796 { | |
797 int w1, h1, x, y, c_count, c_w, c_h, c_w1, c_h1, h_shift, v_shift, c_idx; | |
798 PIXEL *ptr, v, *ptr1; | |
799 | |
800 assert(img->pixel_shift == 1); | |
801 if (cb_size <= 1) | |
802 return; | |
803 w1 = (img->w + cb_size - 1) & ~(cb_size - 1); | |
804 h1 = (img->h + cb_size - 1) & ~(cb_size - 1); | |
805 | |
806 if (img->format == BPG_FORMAT_GRAY) | |
807 c_count = 1; | |
808 else | |
809 c_count = 3; | |
810 if (img->has_alpha) | |
811 c_count++; | |
812 for(c_idx = 0; c_idx < c_count; c_idx++) { | |
813 if (img->format == BPG_FORMAT_420 && | |
814 (c_idx == 1 || c_idx == 2)) { | |
815 h_shift = 1; | |
816 v_shift = 1; | |
817 } else if (img->format == BPG_FORMAT_422 && | |
818 (c_idx == 1 || c_idx == 2)) { | |
819 h_shift = 1; | |
820 v_shift = 0; | |
821 } else { | |
822 h_shift = 0; | |
823 v_shift = 0; | |
824 } | |
825 | |
826 c_w = (img->w + h_shift) >> h_shift; | |
827 c_h = (img->h + v_shift) >> v_shift; | |
828 c_w1 = w1 >> h_shift; | |
829 c_h1 = h1 >> v_shift; | |
830 | |
831 /* pad horizontally */ | |
832 for(y = 0; y < c_h; y++) { | |
833 ptr = (PIXEL *)(img->data[c_idx] + img->linesize[c_idx] * y); | |
834 v = ptr[c_w - 1]; | |
835 for(x = c_w; x < c_w1; x++) { | |
836 ptr[x] = v; | |
837 } | |
838 } | |
839 | |
840 /* pad vertically */ | |
841 ptr1 = (PIXEL *)(img->data[c_idx] + img->linesize[c_idx] * (c_h - 1)); | |
842 for(y = c_h; y < c_h1; y++) { | |
843 ptr = (PIXEL *)(img->data[c_idx] + img->linesize[c_idx] * y); | |
844 memcpy(ptr, ptr1, c_w1 * sizeof(PIXEL)); | |
845 } | |
846 } | |
847 img->w = w1; | |
848 img->h = h1; | |
849 } | |
850 | |
851 /* convert the 16 bit components to 8 bits */ | |
852 void image_convert16to8(Image *img) | |
853 { | |
854 int w, h, stride, y, x, c_count, i; | |
855 uint8_t *plane; | |
856 | |
857 if (img->bit_depth > 8 || img->pixel_shift != 1) | |
858 return; | |
859 if (img->format == BPG_FORMAT_GRAY) | |
860 c_count = 1; | |
861 else | |
862 c_count = 3; | |
863 if (img->has_alpha) | |
864 c_count++; | |
865 for(i = 0; i < c_count; i++) { | |
866 get_plane_res(img, &w, &h, i); | |
867 stride = w; | |
868 plane = malloc(stride * h); | |
869 for(y = 0; y < h; y++) { | |
870 const uint16_t *src; | |
871 uint8_t *dst; | |
872 dst = plane + stride * y; | |
873 src = (uint16_t *)(img->data[i] + img->linesize[i] * y); | |
874 for(x = 0; x < w; x++) | |
875 dst[x] = src[x]; | |
876 } | |
877 free(img->data[i]); | |
878 img->data[i] = plane; | |
879 img->linesize[i] = stride; | |
880 } | |
881 img->pixel_shift = 0; | |
882 } | |
883 | |
884 typedef struct BPGMetaData { | |
885 uint32_t tag; | |
886 uint8_t *buf; | |
887 int buf_len; | |
888 struct BPGMetaData *next; | |
889 } BPGMetaData; | |
890 | |
891 BPGMetaData *bpg_md_alloc(uint32_t tag) | |
892 { | |
893 BPGMetaData *md; | |
894 md = malloc(sizeof(BPGMetaData)); | |
895 memset(md, 0, sizeof(*md)); | |
896 md->tag = tag; | |
897 return md; | |
898 } | |
899 | |
900 void bpg_md_free(BPGMetaData *md) | |
901 { | |
902 BPGMetaData *md_next; | |
903 | |
904 while (md != NULL) { | |
905 md_next = md->next; | |
906 free(md->buf); | |
907 free(md); | |
908 md = md_next; | |
909 } | |
910 } | |
911 | |
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
912 Image *read_tiff(BPGMetaData **pmd, const char *filename, |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
913 BPGColorSpaceEnum color_space, int out_bit_depth, |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
914 int limited_range, int premultiplied_alpha) |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
915 { |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
916 TIFF *tif; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
917 Image *img = NULL; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
918 BPGImageFormatEnum img_format; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
919 uint8_t *buf = NULL; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
920 uint32_t img_width, img_height; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
921 int img_alpha, img_bpp; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
922 size_t img_linesize; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
923 uint16_t img_depth, img_spp, img_sfmt, img_pconfig; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
924 ColorConvertState cvt; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
925 RGBConvertFunc *convert_func; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
926 |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
927 *pmd = NULL; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
928 |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
929 // Open and read TIFF header etc. |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
930 if ((tif = TIFFOpen(filename, "rb")) == NULL) |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
931 return NULL; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
932 |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
933 TIFFGetField(tif, TIFFTAG_IMAGEWIDTH, &img_width); |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
934 TIFFGetField(tif, TIFFTAG_IMAGELENGTH, &img_height); |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
935 TIFFGetField(tif, TIFFTAG_SAMPLESPERPIXEL, &img_spp); |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
936 TIFFGetField(tif, TIFFTAG_BITSPERSAMPLE, &img_depth); |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
937 TIFFGetField(tif, TIFFTAG_SAMPLEFORMAT, &img_sfmt); |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
938 TIFFGetField(tif, TIFFTAG_PLANARCONFIG, &img_pconfig); |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
939 |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
940 // Check basics |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
941 if ((img_depth != 8 && img_depth != 16) || img_spp > 4) |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
942 goto err; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
943 |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
944 img_linesize = TIFFScanlineSize(tif); |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
945 img_format = BPG_FORMAT_444; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
946 img_bpp = sizeof(uint32_t) * img_depth / 8; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
947 |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
948 if (img_spp > 3) |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
949 { |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
950 uint16_t img_esmp, *img_esmp_types; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
951 |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
952 TIFFGetField(tif, TIFFTAG_EXTRASAMPLES, &img_esmp, &img_esmp_types); |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
953 img_alpha = img_esmp == 1 && img_esmp_types[0] == EXTRASAMPLE_UNASSALPHA; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
954 } |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
955 else |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
956 img_alpha = 0; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
957 |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
958 // Allocate temporary image space |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
959 buf = (uint8_t *) _TIFFmalloc(img_linesize * img_spp); |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
960 if (buf == NULL) |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
961 goto err; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
962 |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
963 // Allocate target image |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
964 img = image_alloc(img_width, img_height, img_format, img_alpha, color_space, out_bit_depth); |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
965 if (img == NULL) |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
966 goto err; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
967 |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
968 img->limited_range = limited_range; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
969 img->premultiplied_alpha = premultiplied_alpha; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
970 |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
971 convert_init(&cvt, img_depth, out_bit_depth, color_space, limited_range); |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
972 |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
973 int idx = img_depth == 16; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
974 convert_func = rgb_to_cs[idx][color_space]; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
975 |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
976 for (int y = 0; y < img->h; y++) |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
977 { |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
978 // Read TIFF image to raster |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
979 if (!TIFFReadScanline(tif, buf, y, 0)) |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
980 goto err; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
981 |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
982 // Convert |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
983 convert_func(&cvt, |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
984 (PIXEL *)(img->data[0] + y * img->linesize[0]), |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
985 (PIXEL *)(img->data[1] + y * img->linesize[1]), |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
986 (PIXEL *)(img->data[2] + y * img->linesize[2]), |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
987 buf, img->w, img_spp); |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
988 } |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
989 |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
990 |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
991 err: |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
992 if (buf != NULL) |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
993 _TIFFfree(buf); |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
994 |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
995 TIFFClose(tif); |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
996 return img; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
997 } |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
998 |
0 | 999 Image *read_png(BPGMetaData **pmd, |
1000 FILE *f, BPGColorSpaceEnum color_space, int out_bit_depth, | |
1001 int limited_range, int premultiplied_alpha) | |
1002 { | |
1003 png_structp png_ptr; | |
1004 png_infop info_ptr; | |
1005 int bit_depth, color_type; | |
1006 Image *img; | |
1007 uint8_t **rows; | |
1008 int y, has_alpha, linesize, bpp; | |
1009 BPGImageFormatEnum format; | |
1010 ColorConvertState cvt_s, *cvt = &cvt_s; | |
1011 BPGMetaData *md, **plast_md, *first_md; | |
1012 | |
1013 png_ptr = png_create_read_struct(PNG_LIBPNG_VER_STRING, | |
1014 NULL, NULL, NULL); | |
1015 if (png_ptr == NULL) { | |
1016 return NULL; | |
1017 } | |
1018 | |
1019 info_ptr = png_create_info_struct(png_ptr); | |
1020 if (info_ptr == NULL) { | |
1021 png_destroy_read_struct(&png_ptr, NULL, NULL); | |
1022 return NULL; | |
1023 } | |
1024 | |
1025 if (setjmp(png_jmpbuf(png_ptr))) { | |
1026 png_destroy_read_struct(&png_ptr, &info_ptr, NULL); | |
1027 return NULL; | |
1028 } | |
1029 | |
1030 png_init_io(png_ptr, f); | |
1031 | |
1032 png_read_info(png_ptr, info_ptr); | |
1033 | |
1034 bit_depth = png_get_bit_depth(png_ptr, info_ptr); | |
1035 color_type = png_get_color_type(png_ptr, info_ptr); | |
1036 | |
1037 switch (color_type) { | |
1038 case PNG_COLOR_TYPE_PALETTE: | |
1039 png_set_palette_to_rgb(png_ptr); | |
1040 bit_depth = 8; | |
1041 break; | |
1042 case PNG_COLOR_TYPE_GRAY: | |
1043 case PNG_COLOR_TYPE_GRAY_ALPHA: | |
1044 if (bit_depth < 8) { | |
1045 png_set_expand_gray_1_2_4_to_8(png_ptr); | |
1046 bit_depth = 8; | |
1047 } | |
1048 break; | |
1049 } | |
1050 assert(bit_depth == 8 || bit_depth == 16); | |
1051 | |
1052 #if __BYTE_ORDER__ != __ORDER_BIG_ENDIAN__ | |
1053 if (bit_depth == 16) { | |
1054 png_set_swap(png_ptr); | |
1055 } | |
1056 #endif | |
1057 | |
1058 if (color_type == PNG_COLOR_TYPE_GRAY || | |
1059 color_type == PNG_COLOR_TYPE_GRAY_ALPHA) { | |
1060 format = BPG_FORMAT_GRAY; | |
1061 color_space = BPG_CS_YCbCr; | |
1062 } else { | |
1063 format = BPG_FORMAT_444; | |
1064 } | |
1065 | |
1066 has_alpha = (color_type == PNG_COLOR_TYPE_GRAY_ALPHA || | |
1067 color_type == PNG_COLOR_TYPE_RGB_ALPHA); | |
1068 | |
1069 if (png_get_valid(png_ptr, info_ptr, PNG_INFO_tRNS)) { | |
1070 png_set_tRNS_to_alpha(png_ptr); | |
1071 has_alpha = 1; | |
1072 } | |
1073 | |
1074 if (premultiplied_alpha) { | |
1075 png_set_alpha_mode(png_ptr, PNG_ALPHA_ASSOCIATED, PNG_GAMMA_LINEAR); | |
1076 } | |
1077 | |
1078 img = image_alloc(png_get_image_width(png_ptr, info_ptr), | |
1079 png_get_image_height(png_ptr, info_ptr), | |
1080 format, has_alpha, color_space, | |
1081 out_bit_depth); | |
1082 img->limited_range = limited_range; | |
1083 img->premultiplied_alpha = premultiplied_alpha; | |
1084 | |
1085 rows = malloc(sizeof(rows[0]) * img->h); | |
1086 if (format == BPG_FORMAT_GRAY) | |
1087 bpp = (1 + has_alpha) * (bit_depth / 8); | |
1088 else | |
1089 bpp = (3 + has_alpha) * (bit_depth / 8); | |
1090 linesize = bpp * img->w; | |
1091 for (y = 0; y < img->h; y++) { | |
1092 rows[y] = malloc(linesize); | |
1093 } | |
1094 | |
1095 png_read_image(png_ptr, rows); | |
1096 | |
1097 convert_init(cvt, bit_depth, out_bit_depth, color_space, limited_range); | |
1098 | |
1099 if (format != BPG_FORMAT_GRAY) { | |
1100 int idx; | |
1101 RGBConvertFunc *convert_func; | |
1102 | |
1103 idx = (bit_depth == 16); | |
1104 convert_func = rgb_to_cs[idx][color_space]; | |
1105 | |
1106 for (y = 0; y < img->h; y++) { | |
1107 convert_func(cvt, (PIXEL *)(img->data[0] + y * img->linesize[0]), | |
1108 (PIXEL *)(img->data[1] + y * img->linesize[1]), | |
1109 (PIXEL *)(img->data[2] + y * img->linesize[2]), | |
1110 rows[y], img->w, 3 + has_alpha); | |
1111 if (has_alpha) { | |
1112 if (idx) { | |
1113 gray16_to_gray(cvt, (PIXEL *)(img->data[3] + y * img->linesize[3]), | |
1114 (uint16_t *)rows[y] + 3, img->w, 4); | |
1115 } else { | |
1116 gray8_to_gray(cvt, (PIXEL *)(img->data[3] + y * img->linesize[3]), | |
1117 rows[y] + 3, img->w, 4); | |
1118 } | |
1119 } | |
1120 } | |
1121 } else { | |
1122 if (bit_depth == 16) { | |
1123 for (y = 0; y < img->h; y++) { | |
1124 luma16_to_gray(cvt, (PIXEL *)(img->data[0] + y * img->linesize[0]), | |
1125 (uint16_t *)rows[y], img->w, 1 + has_alpha); | |
1126 if (has_alpha) { | |
1127 gray16_to_gray(cvt, (PIXEL *)(img->data[1] + y * img->linesize[1]), | |
1128 (uint16_t *)rows[y] + 1, img->w, 2); | |
1129 } | |
1130 } | |
1131 } else { | |
1132 for (y = 0; y < img->h; y++) { | |
1133 luma8_to_gray(cvt, (PIXEL *)(img->data[0] + y * img->linesize[0]), | |
1134 rows[y], img->w, 1 + has_alpha); | |
1135 if (has_alpha) { | |
1136 gray8_to_gray(cvt, (PIXEL *)(img->data[1] + y * img->linesize[1]), | |
1137 rows[y] + 1, img->w, 2); | |
1138 } | |
1139 } | |
1140 } | |
1141 } | |
1142 | |
1143 for (y = 0; y < img->h; y++) { | |
1144 free(rows[y]); | |
1145 } | |
1146 free(rows); | |
1147 | |
1148 png_read_end(png_ptr, info_ptr); | |
1149 | |
1150 /* get the ICC profile if present */ | |
1151 first_md = NULL; | |
1152 plast_md = &first_md; | |
1153 { | |
1154 png_charp name; | |
1155 int comp_type; | |
1156 png_bytep iccp_buf; | |
1157 png_uint_32 iccp_buf_len; | |
1158 | |
1159 if (png_get_iCCP(png_ptr, info_ptr, | |
1160 &name, &comp_type, &iccp_buf, &iccp_buf_len) == | |
1161 PNG_INFO_iCCP) { | |
1162 md = bpg_md_alloc(BPG_EXTENSION_TAG_ICCP); | |
1163 md->buf_len = iccp_buf_len; | |
1164 md->buf = malloc(iccp_buf_len); | |
1165 memcpy(md->buf, iccp_buf, iccp_buf_len); | |
1166 *plast_md = md; | |
1167 plast_md = &md->next; | |
1168 } | |
1169 } | |
1170 | |
1171 png_destroy_read_struct(&png_ptr, &info_ptr, NULL); | |
1172 | |
1173 *pmd = first_md; | |
1174 return img; | |
1175 } | |
1176 | |
1177 static BPGMetaData *jpeg_get_metadata(jpeg_saved_marker_ptr first_marker) | |
1178 { | |
1179 static const char app1_exif[] = "Exif"; | |
1180 static const char app1_xmp[] = "http://ns.adobe.com/xap/1.0/"; | |
1181 static const char app2_iccp[] = "ICC_PROFILE"; | |
1182 jpeg_saved_marker_ptr marker; | |
1183 BPGMetaData *md, **plast_md, *first_md; | |
1184 int has_exif, has_xmp, l, iccp_chunk_count, i; | |
1185 jpeg_saved_marker_ptr iccp_chunks[256]; | |
1186 | |
1187 iccp_chunk_count = 0; | |
1188 has_exif = 0; | |
1189 has_xmp = 0; | |
1190 first_md = NULL; | |
1191 plast_md = &first_md; | |
1192 for (marker = first_marker; marker != NULL; marker = marker->next) { | |
1193 #if 0 | |
1194 printf("marker=APP%d len=%d\n", | |
1195 marker->marker - JPEG_APP0, marker->data_length); | |
1196 #endif | |
1197 if (!has_exif && marker->marker == JPEG_APP0 + 1 && | |
1198 marker->data_length > sizeof(app1_exif) && | |
1199 !memcmp(marker->data, app1_exif, sizeof(app1_exif))) { | |
1200 md = bpg_md_alloc(BPG_EXTENSION_TAG_EXIF); | |
1201 l = sizeof(app1_exif); | |
1202 md->buf_len = marker->data_length - l; | |
1203 md->buf = malloc(md->buf_len); | |
1204 memcpy(md->buf, marker->data + l, md->buf_len); | |
1205 *plast_md = md; | |
1206 plast_md = &md->next; | |
1207 has_exif = 1; | |
1208 } else if (!has_xmp && marker->marker == JPEG_APP0 + 1 && | |
1209 marker->data_length > sizeof(app1_xmp) && | |
1210 !memcmp(marker->data, app1_xmp, sizeof(app1_xmp)) && | |
1211 !has_xmp) { | |
1212 md = bpg_md_alloc(BPG_EXTENSION_TAG_XMP); | |
1213 l = sizeof(app1_xmp); | |
1214 md->buf_len = marker->data_length - l; | |
1215 md->buf = malloc(md->buf_len); | |
1216 memcpy(md->buf, marker->data + l, md->buf_len); | |
1217 *plast_md = md; | |
1218 plast_md = &md->next; | |
1219 has_xmp = 1; | |
1220 } else if (marker->marker == JPEG_APP0 + 2 && | |
1221 marker->data_length > (sizeof(app2_iccp) + 2) && | |
1222 !memcmp(marker->data, app2_iccp, sizeof(app2_iccp))) { | |
1223 int chunk_count, chunk_index; | |
1224 l = sizeof(app2_iccp); | |
1225 chunk_index = marker->data[l]; | |
1226 chunk_count = marker->data[l]; | |
1227 if (chunk_index == 0 || chunk_count == 0) | |
1228 continue; | |
1229 if (iccp_chunk_count == 0) { | |
1230 iccp_chunk_count = chunk_count; | |
1231 for(i = 0; i < chunk_count; i++) { | |
1232 iccp_chunks[i] = NULL; | |
1233 } | |
1234 } else { | |
1235 if (chunk_count != iccp_chunk_count) | |
1236 continue; | |
1237 } | |
1238 if (chunk_index > iccp_chunk_count) | |
1239 continue; | |
1240 iccp_chunks[chunk_index - 1] = marker; | |
1241 } | |
1242 } | |
1243 | |
1244 if (iccp_chunk_count != 0) { | |
1245 int len, hlen, idx; | |
1246 /* check that no chunk are missing */ | |
1247 len = 0; | |
1248 hlen = sizeof(app2_iccp) + 2; | |
1249 for(i = 0; i < iccp_chunk_count; i++) { | |
1250 if (!iccp_chunks[i]) | |
1251 break; | |
1252 len += iccp_chunks[i]->data_length - hlen; | |
1253 } | |
1254 if (i == iccp_chunk_count) { | |
1255 md = bpg_md_alloc(BPG_EXTENSION_TAG_ICCP); | |
1256 md->buf_len = len; | |
1257 md->buf = malloc(md->buf_len); | |
1258 idx = 0; | |
1259 for(i = 0; i < iccp_chunk_count; i++) { | |
1260 l = iccp_chunks[i]->data_length - hlen; | |
1261 memcpy(md->buf + idx, iccp_chunks[i]->data + hlen, l); | |
1262 idx += l; | |
1263 } | |
1264 assert(idx == len); | |
1265 *plast_md = md; | |
1266 plast_md = &md->next; | |
1267 } | |
1268 } | |
1269 return first_md; | |
1270 } | |
1271 | |
1272 Image *read_jpeg(BPGMetaData **pmd, FILE *f, | |
1273 int out_bit_depth) | |
1274 { | |
1275 struct jpeg_decompress_struct cinfo; | |
1276 struct jpeg_error_mgr jerr; | |
1277 int w, h, w1, i, y_h, c_h, y, v_shift, c_w, y1, idx, c_idx, h_shift; | |
1278 int h1, plane_idx[4], has_alpha, has_w_plane; | |
1279 Image *img; | |
1280 BPGImageFormatEnum format; | |
1281 BPGColorSpaceEnum color_space; | |
1282 ColorConvertState cvt_s, *cvt = &cvt_s; | |
1283 BPGMetaData *first_md = NULL; | |
1284 uint32_t comp_hv; | |
1285 | |
1286 cinfo.err = jpeg_std_error(&jerr); | |
1287 jpeg_create_decompress(&cinfo); | |
1288 | |
1289 jpeg_save_markers(&cinfo, JPEG_APP0 + 1, 65535); | |
1290 jpeg_save_markers(&cinfo, JPEG_APP0 + 2, 65535); | |
1291 | |
1292 jpeg_stdio_src(&cinfo, f); | |
1293 | |
1294 jpeg_read_header(&cinfo, TRUE); | |
1295 | |
1296 cinfo.raw_data_out = TRUE; | |
1297 cinfo.do_fancy_upsampling = TRUE; | |
1298 | |
1299 w = cinfo.image_width; | |
1300 h = cinfo.image_height; | |
1301 | |
1302 has_w_plane = 0; | |
1303 comp_hv = 0; | |
1304 if (cinfo.num_components < 1 || cinfo.num_components > 4) | |
1305 goto unsupported; | |
1306 for(i = 0; i < cinfo.num_components; i++) { | |
1307 comp_hv |= cinfo.comp_info[i].h_samp_factor << (i * 8 + 4); | |
1308 comp_hv |= cinfo.comp_info[i].v_samp_factor << (i * 8); | |
1309 } | |
1310 switch(cinfo.jpeg_color_space) { | |
1311 case JCS_GRAYSCALE: | |
1312 if (cinfo.num_components != 1 || comp_hv != 0x11) | |
1313 goto unsupported; | |
1314 format = BPG_FORMAT_GRAY; | |
1315 color_space = BPG_CS_YCbCr; | |
1316 break; | |
1317 case JCS_YCbCr: | |
1318 if (cinfo.num_components != 3) | |
1319 goto unsupported; | |
1320 switch(comp_hv) { | |
1321 case 0x111111: | |
1322 format = BPG_FORMAT_444; | |
1323 break; | |
1324 case 0x111121: | |
1325 format = BPG_FORMAT_422; | |
1326 break; | |
1327 case 0x111122: | |
1328 format = BPG_FORMAT_420; | |
1329 break; | |
1330 default: | |
1331 cinfo.raw_data_out = FALSE; | |
1332 format = BPG_FORMAT_444; | |
1333 cinfo.out_color_space = JCS_YCbCr; | |
1334 break; | |
1335 } | |
1336 color_space = BPG_CS_YCbCr; | |
1337 break; | |
1338 case JCS_RGB: | |
1339 if (cinfo.num_components != 3) | |
1340 goto unsupported; | |
1341 format = BPG_FORMAT_444; | |
1342 color_space = BPG_CS_RGB; | |
1343 cinfo.raw_data_out = FALSE; | |
1344 cinfo.out_color_space = JCS_RGB; | |
1345 break; | |
1346 case JCS_YCCK: | |
1347 if (cinfo.num_components != 4) | |
1348 goto unsupported; | |
1349 switch(comp_hv) { | |
1350 case 0x11111111: | |
1351 format = BPG_FORMAT_444; | |
1352 color_space = BPG_CS_YCbCr; | |
1353 break; | |
1354 case 0x22111121: | |
1355 format = BPG_FORMAT_422; | |
1356 color_space = BPG_CS_YCbCr; | |
1357 break; | |
1358 case 0x22111122: | |
1359 format = BPG_FORMAT_420; | |
1360 color_space = BPG_CS_YCbCr; | |
1361 break; | |
1362 default: | |
1363 cinfo.raw_data_out = FALSE; | |
1364 format = BPG_FORMAT_444; | |
1365 cinfo.out_color_space = JCS_CMYK; | |
1366 color_space = BPG_CS_RGB; | |
1367 break; | |
1368 } | |
1369 has_w_plane = 1; | |
1370 break; | |
1371 case JCS_CMYK: | |
1372 if (cinfo.num_components != 4) | |
1373 goto unsupported; | |
1374 format = BPG_FORMAT_444; | |
1375 color_space = BPG_CS_RGB; | |
1376 has_w_plane = 1; | |
1377 cinfo.raw_data_out = FALSE; | |
1378 cinfo.out_color_space = JCS_CMYK; | |
1379 break; | |
1380 default: | |
1381 unsupported: | |
1382 fprintf(stderr, "Unsupported JPEG parameters (cs=%d n_comp=%d comp_hv=%x)\n", | |
1383 cinfo.jpeg_color_space, cinfo.num_components, comp_hv); | |
1384 img = NULL; | |
1385 goto the_end; | |
1386 } | |
1387 | |
1388 v_shift = (format == BPG_FORMAT_420); | |
1389 h_shift = (format == BPG_FORMAT_422 || format == BPG_FORMAT_420); | |
1390 has_alpha = (cinfo.num_components == 4); | |
1391 img = image_alloc(w, h, format, has_alpha, color_space, out_bit_depth); | |
1392 img->has_w_plane = has_w_plane; | |
1393 | |
1394 convert_init(cvt, 8, out_bit_depth, color_space, 0); | |
1395 | |
1396 jpeg_start_decompress(&cinfo); | |
1397 | |
1398 if (color_space == BPG_CS_RGB) { | |
1399 plane_idx[0] = 2; | |
1400 plane_idx[1] = 0; | |
1401 plane_idx[2] = 1; | |
1402 } else { | |
1403 plane_idx[0] = 0; | |
1404 plane_idx[1] = 1; | |
1405 plane_idx[2] = 2; | |
1406 } | |
1407 plane_idx[3] = 3; | |
1408 | |
1409 if (cinfo.raw_data_out) { | |
1410 JSAMPROW rows[4][16]; | |
1411 JSAMPROW *plane_pointer[4]; | |
1412 | |
1413 y_h = 8 * cinfo.max_v_samp_factor; | |
1414 if (cinfo.num_components == 1) { | |
1415 c_h = 0; | |
1416 c_w = 0; | |
1417 } else { | |
1418 c_h = 8; | |
1419 if (h_shift) | |
1420 c_w = (w + 1) / 2; | |
1421 else | |
1422 c_w = w; | |
1423 } | |
1424 w1 = (w + 15) & ~15; | |
1425 for(c_idx = 0; c_idx < cinfo.num_components; c_idx++) { | |
1426 if (c_idx == 1 || c_idx == 2) { | |
1427 h1 = c_h; | |
1428 } else { | |
1429 h1 = y_h; | |
1430 } | |
1431 for(i = 0; i < h1; i++) { | |
1432 rows[c_idx][i] = malloc(w1); | |
1433 } | |
1434 plane_pointer[c_idx] = rows[c_idx]; | |
1435 } | |
1436 | |
1437 while (cinfo.output_scanline < cinfo.output_height) { | |
1438 y = cinfo.output_scanline; | |
1439 jpeg_read_raw_data(&cinfo, plane_pointer, y_h); | |
1440 | |
1441 for(c_idx = 0; c_idx < cinfo.num_components; c_idx++) { | |
1442 if (c_idx == 1 || c_idx == 2) { | |
1443 h1 = c_h; | |
1444 w1 = c_w; | |
1445 y1 = (y >> v_shift); | |
1446 } else { | |
1447 h1 = y_h; | |
1448 w1 = img->w; | |
1449 y1 = y; | |
1450 } | |
1451 idx = plane_idx[c_idx]; | |
1452 for(i = 0; i < h1; i++) { | |
1453 PIXEL *ptr; | |
1454 ptr = (PIXEL *)(img->data[idx] + | |
1455 img->linesize[idx] * (y1 + i)); | |
1456 gray8_to_gray(cvt, ptr, rows[c_idx][i], w1, 1); | |
1457 if (color_space == BPG_CS_YCbCr && has_w_plane) { | |
1458 /* negate color */ | |
1459 if (c_idx == 0) { | |
1460 gray_one_minus(cvt, ptr, w1); | |
1461 } else if (c_idx <= 2) { | |
1462 gray_neg_c(cvt, ptr, w1); | |
1463 } | |
1464 } | |
1465 } | |
1466 } | |
1467 } | |
1468 | |
1469 for(c_idx = 0; c_idx < cinfo.num_components; c_idx++) { | |
1470 if (c_idx == 1 || c_idx == 2) { | |
1471 h1 = c_h; | |
1472 } else { | |
1473 h1 = y_h; | |
1474 } | |
1475 for(i = 0; i < h1; i++) { | |
1476 free(rows[c_idx][i]); | |
1477 } | |
1478 } | |
1479 } else { | |
1480 JSAMPROW rows[1]; | |
1481 uint8_t *buf; | |
1482 int c_count; | |
1483 | |
1484 c_count = 3 + has_w_plane; | |
1485 buf = malloc(c_count * w); | |
1486 rows[0] = buf; | |
1487 while (cinfo.output_scanline < cinfo.output_height) { | |
1488 y = cinfo.output_scanline; | |
1489 jpeg_read_scanlines(&cinfo, rows, 1); | |
1490 | |
1491 for(c_idx = 0; c_idx < c_count; c_idx++) { | |
1492 idx = plane_idx[c_idx]; | |
1493 gray8_to_gray(cvt, (PIXEL *)(img->data[idx] + | |
1494 img->linesize[idx] * y), | |
1495 buf + c_idx, w, c_count); | |
1496 } | |
1497 } | |
1498 free(buf); | |
1499 } | |
1500 | |
1501 first_md = jpeg_get_metadata(cinfo.marker_list); | |
1502 | |
1503 the_end: | |
1504 jpeg_finish_decompress(&cinfo); | |
1505 | |
1506 jpeg_destroy_decompress(&cinfo); | |
1507 *pmd = first_md; | |
1508 return img; | |
1509 } | |
1510 | |
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1511 static const char PROBE_TIFF_MAGIC_LE[] = { 'I', 'I', 42, 0 }; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1512 static const char PROBE_TIFF_MAGIC_BE[] = { 'M', 'M', 0, 42 }; |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1513 #define PROBE_BUF_SIZE 8 |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1514 |
0 | 1515 Image *load_image(BPGMetaData **pmd, const char *infilename, |
1516 BPGColorSpaceEnum color_space, int bit_depth, | |
1517 int limited_range, int premultiplied_alpha) | |
1518 { | |
1519 FILE *f; | |
3
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1520 Image *img = NULL; |
0 | 1521 BPGMetaData *md; |
3
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1522 uint8_t buf[PROBE_BUF_SIZE]; |
0 | 1523 |
1524 *pmd = NULL; | |
1525 | |
3
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1526 if ((f = fopen(infilename, "rb")) == NULL) |
0 | 1527 return NULL; |
3
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1528 |
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1529 if (fread(buf, 1, PROBE_BUF_SIZE, f) != PROBE_BUF_SIZE) |
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1530 goto err; |
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1531 |
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1532 fseek(f, 0, SEEK_SET); |
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1533 |
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1534 if (png_sig_cmp(buf, 0, PROBE_BUF_SIZE) == 0) |
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1535 img = read_png(&md, f, color_space, bit_depth, limited_range, premultiplied_alpha); |
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1536 else |
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1537 if (memcmp(buf, PROBE_TIFF_MAGIC_LE, sizeof(PROBE_TIFF_MAGIC_LE)) == 0 || |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1538 memcmp(buf, PROBE_TIFF_MAGIC_BE, sizeof(PROBE_TIFF_MAGIC_BE)) == 0) |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1539 img = read_tiff(&md, infilename, color_space, bit_depth, limited_range, premultiplied_alpha); |
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1540 else |
0 | 1541 img = read_jpeg(&md, f, bit_depth); |
3
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1542 |
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1543 err: |
0 | 1544 fclose(f); |
1545 *pmd = md; | |
1546 return img; | |
1547 } | |
1548 | |
1549 void save_yuv1(Image *img, FILE *f) | |
1550 { | |
1551 int c_w, c_h, i, c_count, y; | |
1552 | |
1553 if (img->format == BPG_FORMAT_GRAY) | |
1554 c_count = 1; | |
1555 else | |
1556 c_count = 3; | |
1557 for(i = 0; i < c_count; i++) { | |
1558 get_plane_res(img, &c_w, &c_h, i); | |
1559 for(y = 0; y < c_h; y++) { | |
1560 fwrite(img->data[i] + y * img->linesize[i], | |
1561 1, c_w << img->pixel_shift, f); | |
1562 } | |
1563 } | |
1564 } | |
1565 | |
1566 void save_yuv(Image *img, const char *filename) | |
1567 { | |
1568 FILE *f; | |
1569 | |
1570 f = fopen(filename, "wb"); | |
1571 if (!f) { | |
1572 fprintf(stderr, "Could not open %s\n", filename); | |
1573 exit(1); | |
1574 } | |
1575 save_yuv1(img, f); | |
1576 fclose(f); | |
1577 } | |
1578 | |
1579 /* return the position of the end of the NAL or -1 if error */ | |
1580 static int find_nal_end(const uint8_t *buf, int buf_len) | |
1581 { | |
1582 int idx; | |
1583 | |
1584 idx = 0; | |
1585 if (buf_len >= 4 && | |
1586 buf[0] == 0 && buf[1] == 0 && buf[2] == 0 && buf[3] == 1) { | |
1587 idx = 4; | |
1588 } else if (buf_len >= 3 && | |
1589 buf[0] == 0 && buf[1] == 0 && buf[2] == 1) { | |
1590 idx = 3; | |
1591 } else { | |
1592 return -1; | |
1593 } | |
1594 /* NAL header */ | |
1595 if (idx + 2 > buf_len) | |
1596 return -1; | |
1597 /* find the last byte */ | |
1598 for(;;) { | |
1599 if (idx + 2 >= buf_len) { | |
1600 idx = buf_len; | |
1601 break; | |
1602 } | |
1603 if (buf[idx] == 0 && buf[idx + 1] == 0 && buf[idx + 2] == 1) | |
1604 break; | |
1605 if (idx + 3 < buf_len && | |
1606 buf[idx] == 0 && buf[idx + 1] == 0 && buf[idx + 2] == 0 && buf[idx + 3] == 1) | |
1607 break; | |
1608 idx++; | |
1609 } | |
1610 return idx; | |
1611 } | |
1612 | |
1613 /* return the position of the end of the NAL or -1 if error */ | |
1614 static int extract_nal(uint8_t **pnal_buf, int *pnal_len, | |
1615 const uint8_t *buf, int buf_len) | |
1616 { | |
1617 int idx, start, end, len; | |
1618 uint8_t *nal_buf; | |
1619 int nal_len; | |
1620 | |
1621 end = find_nal_end(buf, buf_len); | |
1622 if (end < 0) | |
1623 return -1; | |
1624 if (buf[2] == 1) | |
1625 start = 3; | |
1626 else | |
1627 start = 4; | |
1628 len = end - start; | |
1629 | |
1630 nal_buf = malloc(len); | |
1631 nal_len = 0; | |
1632 idx = start; | |
1633 while (idx < end) { | |
1634 if (idx + 2 < end && buf[idx] == 0 && buf[idx + 1] == 0 && buf[idx + 2] == 3) { | |
1635 nal_buf[nal_len++] = 0; | |
1636 nal_buf[nal_len++] = 0; | |
1637 idx += 3; | |
1638 } else { | |
1639 nal_buf[nal_len++] = buf[idx++]; | |
1640 } | |
1641 } | |
1642 while (idx < end) { | |
1643 nal_buf[nal_len++] = buf[idx++]; | |
1644 } | |
1645 *pnal_buf = nal_buf; | |
1646 *pnal_len = nal_len; | |
1647 return idx; | |
1648 } | |
1649 | |
1650 /* big endian variable length 7 bit encoding */ | |
1651 static void put_ue(uint8_t **pp, uint32_t v) | |
1652 { | |
1653 uint8_t *p = *pp; | |
1654 int i, j; | |
1655 | |
1656 for(i = 1; i < 5; i++) { | |
1657 if (v < (1 << (7 * i))) | |
1658 break; | |
1659 } | |
1660 for(j = i - 1; j >= 1; j--) | |
1661 *p++ = ((v >> (7 * j)) & 0x7f) | 0x80; | |
1662 *p++ = v & 0x7f; | |
1663 *pp = p; | |
1664 } | |
1665 | |
1666 typedef struct { | |
1667 const uint8_t *buf; | |
1668 int idx; | |
1669 int buf_len; | |
1670 } GetBitState; | |
1671 | |
1672 static void init_get_bits(GetBitState *s, const uint8_t *buf, int buf_len) | |
1673 { | |
1674 s->buf = buf; | |
1675 s->buf_len = buf_len; | |
1676 s->idx = 0; | |
1677 } | |
1678 | |
1679 static void skip_bits(GetBitState *s, int n) | |
1680 { | |
1681 s->idx += n; | |
1682 } | |
1683 | |
1684 /* 1 <= n <= 25. return '0' bits if past the end of the buffer. */ | |
1685 static uint32_t get_bits(GetBitState *s, int n) | |
1686 { | |
1687 const uint8_t *buf = s->buf; | |
1688 int p, i; | |
1689 uint32_t v; | |
1690 | |
1691 p = s->idx >> 3; | |
1692 if ((p + 3) < s->buf_len) { | |
1693 v = (buf[p] << 24) | (buf[p + 1] << 16) | | |
1694 (buf[p + 2] << 8) | buf[p + 3]; | |
1695 } else { | |
1696 v = 0; | |
1697 for(i = 0; i < 3; i++) { | |
1698 if ((p + i) < s->buf_len) | |
1699 v |= buf[p + i] << (24 - i * 8); | |
1700 } | |
1701 } | |
1702 v = (v >> (32 - (s->idx & 7) - n)) & ((1 << n) - 1); | |
1703 s->idx += n; | |
1704 return v; | |
1705 } | |
1706 | |
1707 /* 1 <= n <= 32 */ | |
1708 static uint32_t get_bits_long(GetBitState *s, int n) | |
1709 { | |
1710 uint32_t v; | |
1711 | |
1712 if (n <= 25) { | |
1713 v = get_bits(s, n); | |
1714 } else { | |
1715 n -= 16; | |
1716 v = get_bits(s, 16) << n; | |
1717 v |= get_bits(s, n); | |
1718 } | |
1719 return v; | |
1720 } | |
1721 | |
1722 /* at most 32 bits are supported */ | |
1723 static uint32_t get_ue_golomb(GetBitState *s) | |
1724 { | |
1725 int i; | |
1726 i = 0; | |
1727 for(;;) { | |
1728 if (get_bits(s, 1)) | |
1729 break; | |
1730 i++; | |
1731 if (i == 32) | |
1732 return 0xffffffff; | |
1733 } | |
1734 if (i == 0) | |
1735 return 0; | |
1736 else | |
1737 return ((1 << i) | get_bits_long(s, i)) - 1; | |
1738 } | |
1739 | |
1740 typedef struct { | |
1741 uint8_t *buf; | |
1742 int idx; | |
1743 } PutBitState; | |
1744 | |
1745 static void init_put_bits(PutBitState *s, uint8_t *buf) | |
1746 { | |
1747 s->buf = buf; | |
1748 s->idx = 0; | |
1749 } | |
1750 | |
1751 static void put_bit(PutBitState *s, int bit) | |
1752 { | |
1753 s->buf[s->idx >> 3] |= bit << (7 - (s->idx & 7)); | |
1754 s->idx++; | |
1755 } | |
1756 | |
1757 static void put_bits(PutBitState *s, int n, uint32_t v) | |
1758 { | |
1759 int i; | |
1760 | |
1761 for(i = 0; i < n; i++) { | |
1762 put_bit(s, (v >> (n - 1 - i)) & 1); | |
1763 } | |
1764 } | |
1765 | |
1766 static void put_ue_golomb(PutBitState *s, uint32_t v) | |
1767 { | |
1768 uint32_t a; | |
1769 int n; | |
1770 | |
1771 v++; | |
1772 n = 0; | |
1773 a = v; | |
1774 while (a != 0) { | |
1775 a >>= 1; | |
1776 n++; | |
1777 } | |
1778 if (n > 1) | |
1779 put_bits(s, n - 1, 0); | |
1780 put_bits(s, n, v); | |
1781 } | |
1782 | |
1783 typedef struct { | |
1784 uint8_t *buf; | |
1785 int size; | |
1786 int len; | |
1787 } DynBuf; | |
1788 | |
1789 static void dyn_buf_init(DynBuf *s) | |
1790 { | |
1791 s->buf = NULL; | |
1792 s->size = 0; | |
1793 s->len = 0; | |
1794 } | |
1795 | |
1796 static int dyn_buf_resize(DynBuf *s, int size) | |
1797 { | |
1798 int new_size; | |
1799 uint8_t *new_buf; | |
1800 | |
1801 if (size <= s->size) | |
1802 return 0; | |
1803 new_size = (s->size * 3) / 2; | |
1804 if (new_size < size) | |
1805 new_size = size; | |
1806 new_buf = realloc(s->buf, new_size); | |
1807 if (!new_buf) | |
1808 return -1; | |
1809 s->buf = new_buf; | |
1810 s->size = new_size; | |
1811 return 0; | |
1812 } | |
1813 | |
1814 /* suppress the VPS NAL and keep only the useful part of the SPS | |
1815 header. The decoder can rebuild a valid HEVC stream if needed. */ | |
1816 static int build_modified_sps(uint8_t **pout_buf, int *pout_buf_len, | |
1817 const uint8_t *buf, int buf_len) | |
1818 { | |
1819 int nal_unit_type, nal_len, idx, i, ret, msps_buf_len; | |
1820 int out_buf_len, out_buf_len_max; | |
1821 uint8_t *nal_buf, *msps_buf, *out_buf; | |
1822 GetBitState gb_s, *gb = &gb_s; | |
1823 PutBitState pb_s, *pb = &pb_s; | |
1824 uint8_t *p; | |
1825 | |
1826 idx = extract_nal(&nal_buf, &nal_len, buf, buf_len); | |
1827 if (idx < 0) | |
1828 return -1; | |
1829 if (nal_len < 2) { | |
1830 free(nal_buf); | |
1831 return -1; | |
1832 } | |
1833 nal_unit_type = (nal_buf[0] >> 1) & 0x3f; | |
1834 free(nal_buf); | |
1835 if (nal_unit_type != 32) { | |
1836 fprintf(stderr, "expecting VPS nal (%d)\n", nal_unit_type); | |
1837 return -1; /* expect VPS nal */ | |
1838 } | |
1839 | |
1840 ret = extract_nal(&nal_buf, &nal_len, buf + idx, buf_len); | |
1841 if (ret < 0) | |
1842 return -1; | |
1843 idx += ret; | |
1844 if (nal_len < 2) | |
1845 return -1; | |
1846 nal_unit_type = (nal_buf[0] >> 1) & 0x3f; | |
1847 if (nal_unit_type != 33) { | |
1848 fprintf(stderr, "expecting SPS nal (%d)\n", nal_unit_type); | |
1849 return -1; /* expect SPS nal */ | |
1850 } | |
1851 | |
1852 /* skip the initial part of the SPS up to and including | |
1853 log2_min_cb_size */ | |
1854 { | |
1855 int vps_id, max_sub_layers, profile_idc, sps_id; | |
1856 int chroma_format_idc, width, height, bit_depth_luma, bit_depth_chroma; | |
1857 int log2_max_poc_lsb, sublayer_ordering_info, log2_min_cb_size; | |
1858 int log2_diff_max_min_coding_block_size, log2_min_tb_size; | |
1859 int log2_diff_max_min_transform_block_size; | |
1860 int max_transform_hierarchy_depth_inter; | |
1861 int max_transform_hierarchy_depth_intra; | |
1862 int scaling_list_enable_flag, amp_enabled_flag, sao_enabled; | |
1863 int pcm_enabled_flag, nb_st_rps; | |
1864 int long_term_ref_pics_present_flag, sps_strong_intra_smoothing_enable_flag, vui_present; | |
1865 int sps_temporal_mvp_enabled_flag; | |
1866 int pcm_sample_bit_depth_luma_minus1; | |
1867 int pcm_sample_bit_depth_chroma_minus1; | |
1868 int log2_min_pcm_luma_coding_block_size_minus3; | |
1869 int log2_diff_max_min_pcm_luma_coding_block_size; | |
1870 int pcm_loop_filter_disabled_flag; | |
1871 int sps_extension_flag, sps_range_extension_flag, sps_extension_7bits; | |
1872 int sps_range_extension_flags; | |
1873 | |
1874 init_get_bits(gb, nal_buf, nal_len); | |
1875 skip_bits(gb, 16); /* nal header */ | |
1876 vps_id = get_bits(gb, 4); | |
1877 if (vps_id != 0) { | |
1878 fprintf(stderr, "VPS id 0 expected\n"); | |
1879 return -1; | |
1880 } | |
1881 max_sub_layers = get_bits(gb, 3); | |
1882 if (max_sub_layers != 0) { | |
1883 fprintf(stderr, "max_sub_layers == 0 expected\n"); | |
1884 return -1; | |
1885 } | |
1886 skip_bits(gb, 1); /* temporal_id_nesting_flag */ | |
1887 /* profile tier level */ | |
1888 skip_bits(gb, 2); /* profile_space */ | |
1889 skip_bits(gb, 1); /* tier_flag */ | |
1890 profile_idc = get_bits(gb, 5); | |
1891 for(i = 0; i < 32; i++) { | |
1892 skip_bits(gb, 1); /* profile_compatibility_flag */ | |
1893 } | |
1894 skip_bits(gb, 1); /* progressive_source_flag */ | |
1895 skip_bits(gb, 1); /* interlaced_source_flag */ | |
1896 skip_bits(gb, 1); /* non_packed_constraint_flag */ | |
1897 skip_bits(gb, 1); /* frame_only_constraint_flag */ | |
1898 skip_bits(gb, 44); /* XXX_reserved_zero_44 */ | |
1899 skip_bits(gb, 8); /* level_idc */ | |
1900 | |
1901 sps_id = get_ue_golomb(gb); | |
1902 if (sps_id != 0) { | |
1903 fprintf(stderr, "SPS id 0 expected (%d)\n", sps_id); | |
1904 return -1; | |
1905 } | |
1906 chroma_format_idc = get_ue_golomb(gb); | |
1907 if (chroma_format_idc == 3) { | |
1908 get_bits(gb, 1); /* separate_colour_plane_flag */ | |
1909 } | |
1910 width = get_ue_golomb(gb); | |
1911 height = get_ue_golomb(gb); | |
1912 /* pic conformance_flag */ | |
1913 if (get_bits(gb, 1)) { | |
1914 get_ue_golomb(gb); /* left_offset */ | |
1915 get_ue_golomb(gb); /* right_offset */ | |
1916 get_ue_golomb(gb); /* top_offset */ | |
1917 get_ue_golomb(gb); /* bottom_offset */ | |
1918 } | |
1919 bit_depth_luma = get_ue_golomb(gb) + 8; | |
1920 bit_depth_chroma = get_ue_golomb(gb) + 8; | |
1921 log2_max_poc_lsb = get_ue_golomb(gb) + 4; | |
1922 if (log2_max_poc_lsb != 8) { | |
1923 fprintf(stderr, "log2_max_poc_lsb must be 8 (%d)\n", log2_max_poc_lsb); | |
1924 return -1; | |
1925 } | |
1926 sublayer_ordering_info = get_bits(gb, 1); | |
1927 get_ue_golomb(gb); /* max_dec_pic_buffering */ | |
1928 get_ue_golomb(gb); /* num_reorder_pics */ | |
1929 get_ue_golomb(gb); /* max_latency_increase */ | |
1930 | |
1931 log2_min_cb_size = get_ue_golomb(gb) + 3; | |
1932 log2_diff_max_min_coding_block_size = get_ue_golomb(gb); | |
1933 log2_min_tb_size = get_ue_golomb(gb) + 2; | |
1934 log2_diff_max_min_transform_block_size = get_ue_golomb(gb); | |
1935 | |
1936 max_transform_hierarchy_depth_inter = get_ue_golomb(gb); | |
1937 max_transform_hierarchy_depth_intra = get_ue_golomb(gb); | |
1938 if (max_transform_hierarchy_depth_inter != max_transform_hierarchy_depth_intra) { | |
1939 fprintf(stderr, "max_transform_hierarchy_depth_inter must be the same as max_transform_hierarchy_depth_intra (%d %d)\n", max_transform_hierarchy_depth_inter, max_transform_hierarchy_depth_intra); | |
1940 return -1; | |
1941 } | |
1942 | |
1943 scaling_list_enable_flag = get_bits(gb, 1); | |
1944 if (scaling_list_enable_flag != 0) { | |
1945 fprintf(stderr, "scaling_list_enable_flag must be 0\n"); | |
1946 return -1; | |
1947 } | |
1948 amp_enabled_flag = get_bits(gb, 1); | |
1949 if (!amp_enabled_flag) { | |
1950 fprintf(stderr, "amp_enabled_flag must be set\n"); | |
1951 return -1; | |
1952 } | |
1953 sao_enabled = get_bits(gb, 1); | |
1954 pcm_enabled_flag = get_bits(gb, 1); | |
1955 if (pcm_enabled_flag) { | |
1956 pcm_sample_bit_depth_luma_minus1 = get_bits(gb, 4); | |
1957 pcm_sample_bit_depth_chroma_minus1 = get_bits(gb, 4); | |
1958 log2_min_pcm_luma_coding_block_size_minus3 = get_ue_golomb(gb); | |
1959 log2_diff_max_min_pcm_luma_coding_block_size = get_ue_golomb(gb); | |
1960 pcm_loop_filter_disabled_flag = get_bits(gb, 1); | |
1961 } | |
1962 nb_st_rps = get_ue_golomb(gb); | |
1963 if (nb_st_rps != 0) { | |
1964 fprintf(stderr, "nb_st_rps must be 0 (%d)\n", nb_st_rps); | |
1965 return -1; | |
1966 } | |
1967 long_term_ref_pics_present_flag = get_bits(gb, 1); | |
1968 if (long_term_ref_pics_present_flag) { | |
1969 fprintf(stderr, "nlong_term_ref_pics_present_flag must be 0 (%d)\n", nb_st_rps); | |
1970 return -1; | |
1971 } | |
1972 sps_temporal_mvp_enabled_flag = get_bits(gb, 1); | |
1973 if (!sps_temporal_mvp_enabled_flag) { | |
1974 fprintf(stderr, "sps_temporal_mvp_enabled_flag must be set\n"); | |
1975 return -1; | |
1976 } | |
1977 sps_strong_intra_smoothing_enable_flag = get_bits(gb, 1); | |
1978 vui_present = get_bits(gb, 1); | |
1979 if (vui_present) { | |
1980 int sar_present, sar_idx, overscan_info_present_flag; | |
1981 int video_signal_type_present_flag, chroma_loc_info_present_flag; | |
1982 int default_display_window_flag, vui_timing_info_present_flag; | |
1983 int vui_poc_proportional_to_timing_flag; | |
1984 int vui_hrd_parameters_present_flag, bitstream_restriction_flag; | |
1985 | |
1986 sar_present = get_bits(gb, 1); | |
1987 if (sar_present) { | |
1988 sar_idx = get_bits(gb, 8); | |
1989 if (sar_idx == 255) { | |
1990 skip_bits(gb, 16); /* sar_num */ | |
1991 skip_bits(gb, 16); /* sar_den */ | |
1992 } | |
1993 } | |
1994 | |
1995 overscan_info_present_flag = get_bits(gb, 1); | |
1996 if (overscan_info_present_flag) { | |
1997 skip_bits(gb, 1); /* overscan_appropriate_flag */ | |
1998 } | |
1999 | |
2000 video_signal_type_present_flag = get_bits(gb, 1); | |
2001 if (video_signal_type_present_flag) { | |
2002 fprintf(stderr, "video_signal_type_present_flag must be 0\n"); | |
2003 return -1; | |
2004 } | |
2005 chroma_loc_info_present_flag = get_bits(gb, 1); | |
2006 if (chroma_loc_info_present_flag) { | |
2007 get_ue_golomb(gb); | |
2008 get_ue_golomb(gb); | |
2009 } | |
2010 skip_bits(gb, 1); /* neutra_chroma_indication_flag */ | |
2011 skip_bits(gb, 1); | |
2012 skip_bits(gb, 1); | |
2013 default_display_window_flag = get_bits(gb, 1); | |
2014 if (default_display_window_flag) { | |
2015 fprintf(stderr, "default_display_window_flag must be 0\n"); | |
2016 return -1; | |
2017 } | |
2018 vui_timing_info_present_flag = get_bits(gb, 1); | |
2019 if (vui_timing_info_present_flag) { | |
2020 skip_bits(gb, 32); | |
2021 skip_bits(gb, 32); | |
2022 vui_poc_proportional_to_timing_flag = get_bits(gb, 1); | |
2023 if (vui_poc_proportional_to_timing_flag) { | |
2024 get_ue_golomb(gb); | |
2025 } | |
2026 vui_hrd_parameters_present_flag = get_bits(gb, 1); | |
2027 if (vui_hrd_parameters_present_flag) { | |
2028 fprintf(stderr, "vui_hrd_parameters_present_flag must be 0\n"); | |
2029 return -1; | |
2030 } | |
2031 } | |
2032 bitstream_restriction_flag = get_bits(gb, 1); | |
2033 if (bitstream_restriction_flag) { | |
2034 skip_bits(gb, 1); | |
2035 skip_bits(gb, 1); | |
2036 skip_bits(gb, 1); | |
2037 get_ue_golomb(gb); | |
2038 get_ue_golomb(gb); | |
2039 get_ue_golomb(gb); | |
2040 get_ue_golomb(gb); | |
2041 get_ue_golomb(gb); | |
2042 } | |
2043 } | |
2044 sps_extension_flag = get_bits(gb, 1); | |
2045 sps_range_extension_flag = 0; | |
2046 sps_range_extension_flags = 0; | |
2047 if (sps_extension_flag) { | |
2048 sps_range_extension_flag = get_bits(gb, 1); | |
2049 sps_extension_7bits = get_bits(gb, 7); | |
2050 if (sps_extension_7bits != 0) { | |
2051 fprintf(stderr, "sps_extension_7bits must be 0\n"); | |
2052 return -1; | |
2053 } | |
2054 if (sps_range_extension_flag) { | |
2055 sps_range_extension_flags = get_bits(gb, 9); | |
2056 if (sps_range_extension_flags & ((1 << (8 - 3)) | | |
2057 (1 << (8 - 4)) | | |
2058 (1 << (8 - 6)) | | |
2059 (1 << (8 - 8)))) { | |
2060 fprintf(stderr, "unsupported range extensions (0x%x)\n", | |
2061 sps_range_extension_flags); | |
2062 return -1; | |
2063 } | |
2064 } | |
2065 } | |
2066 | |
2067 /* build the modified SPS */ | |
2068 msps_buf = malloc(nal_len + 32); | |
2069 memset(msps_buf, 0, nal_len + 16); | |
2070 | |
2071 init_put_bits(pb, msps_buf); | |
2072 put_ue_golomb(pb, log2_min_cb_size - 3); | |
2073 put_ue_golomb(pb, log2_diff_max_min_coding_block_size); | |
2074 put_ue_golomb(pb, log2_min_tb_size - 2); | |
2075 put_ue_golomb(pb, log2_diff_max_min_transform_block_size); | |
2076 put_ue_golomb(pb, max_transform_hierarchy_depth_intra); | |
2077 put_bits(pb, 1, sao_enabled); | |
2078 put_bits(pb, 1, pcm_enabled_flag); | |
2079 if (pcm_enabled_flag) { | |
2080 put_bits(pb, 4, pcm_sample_bit_depth_luma_minus1); | |
2081 put_bits(pb, 4, pcm_sample_bit_depth_chroma_minus1); | |
2082 put_ue_golomb(pb, log2_min_pcm_luma_coding_block_size_minus3); | |
2083 put_ue_golomb(pb, log2_diff_max_min_pcm_luma_coding_block_size); | |
2084 put_bits(pb, 1, pcm_loop_filter_disabled_flag); | |
2085 } | |
2086 put_bits(pb, 1, sps_strong_intra_smoothing_enable_flag); | |
2087 put_bits(pb, 1, sps_extension_flag); | |
2088 if (sps_extension_flag) { | |
2089 put_bits(pb, 1, sps_range_extension_flag); | |
2090 put_bits(pb, 7, 0); | |
2091 if (sps_range_extension_flag) { | |
2092 put_bits(pb, 9, sps_range_extension_flags); | |
2093 } | |
2094 } | |
2095 msps_buf_len = (pb->idx + 7) >> 3; | |
2096 | |
2097 out_buf_len_max = 5 + msps_buf_len; | |
2098 out_buf = malloc(out_buf_len_max); | |
2099 | |
2100 // printf("msps_n_bits=%d\n", pb->idx); | |
2101 p = out_buf; | |
2102 put_ue(&p, msps_buf_len); /* header length */ | |
2103 | |
2104 memcpy(p, msps_buf, msps_buf_len); | |
2105 p += msps_buf_len; | |
2106 | |
2107 out_buf_len = p - out_buf; | |
2108 free(msps_buf); | |
2109 free(nal_buf); | |
2110 } | |
2111 *pout_buf = out_buf; | |
2112 *pout_buf_len = out_buf_len; | |
2113 return idx; | |
2114 } | |
2115 | |
2116 static int add_frame_duration_sei(DynBuf *out_buf, uint16_t frame_ticks) | |
2117 { | |
2118 uint8_t nal_buf[128], *q; | |
2119 int nut, nal_len; | |
2120 | |
2121 q = nal_buf; | |
2122 *q++ = 0x00; | |
2123 *q++ = 0x00; | |
2124 *q++ = 0x01; | |
2125 nut = 39; /* prefix SEI NUT */ | |
2126 *q++ = (nut << 1); | |
2127 *q++ = 1; | |
2128 *q++ = 0xff; /* payload_type = 257 */ | |
2129 *q++ = 0x02; | |
2130 *q++ = 2; /* payload_size = 2 */ | |
2131 *q++ = frame_ticks >> 8; | |
2132 *q++ = frame_ticks; | |
2133 *q++ = 0x80; /* extra '1' bit and align to byte */ | |
2134 /* Note: the 0x00 0x00 b pattern with b <= 3 cannot happen, so no | |
2135 need to escape */ | |
2136 nal_len = q - nal_buf; | |
2137 if (dyn_buf_resize(out_buf, out_buf->len + nal_len) < 0) | |
2138 return -1; | |
2139 memcpy(out_buf->buf + out_buf->len, nal_buf, nal_len); | |
2140 out_buf->len += nal_len; | |
2141 return 0; | |
2142 } | |
2143 | |
2144 static int build_modified_hevc(uint8_t **pout_buf, | |
2145 const uint8_t *cbuf, int cbuf_len, | |
2146 const uint8_t *abuf, int abuf_len, | |
2147 const uint16_t *frame_duration_tab, | |
2148 int frame_count) | |
2149 { | |
2150 DynBuf out_buf_s, *out_buf = &out_buf_s; | |
2151 uint8_t *msps; | |
2152 const uint8_t *nal_buf; | |
2153 int msps_len, cidx, aidx, is_alpha, nal_len, first_nal, start, l, frame_num; | |
2154 | |
2155 dyn_buf_init(out_buf); | |
2156 | |
2157 /* add alpha MSPS */ | |
2158 aidx = 0; /* avoids warning */ | |
2159 if (abuf) { | |
2160 aidx = build_modified_sps(&msps, &msps_len, abuf, abuf_len); | |
2161 if (aidx < 0) | |
2162 goto fail; | |
2163 if (dyn_buf_resize(out_buf, out_buf->len + msps_len) < 0) | |
2164 goto fail; | |
2165 memcpy(out_buf->buf + out_buf->len, msps, msps_len); | |
2166 out_buf->len += msps_len; | |
2167 free(msps); | |
2168 } | |
2169 | |
2170 /* add color MSPS */ | |
2171 cidx = build_modified_sps(&msps, &msps_len, cbuf, cbuf_len); | |
2172 if (cidx < 0) | |
2173 goto fail; | |
2174 if (dyn_buf_resize(out_buf, out_buf->len + msps_len) < 0) | |
2175 goto fail; | |
2176 memcpy(out_buf->buf + out_buf->len, msps, msps_len); | |
2177 out_buf->len += msps_len; | |
2178 free(msps); | |
2179 | |
2180 /* add the remaining NALs, alternating between alpha (if present) | |
2181 and color. */ | |
2182 is_alpha = (abuf != NULL); | |
2183 first_nal = 1; | |
2184 frame_num = 0; | |
2185 for(;;) { | |
2186 if (!is_alpha) { | |
2187 if (cidx >= cbuf_len) { | |
2188 if (abuf) { | |
2189 fprintf(stderr, "Incorrect number of alpha NALs\n"); | |
2190 goto fail; | |
2191 } | |
2192 break; | |
2193 } | |
2194 nal_buf = cbuf + cidx; | |
2195 nal_len = find_nal_end(nal_buf, cbuf_len - cidx); | |
2196 // printf("cidx=%d/%d nal_len=%d\n", cidx, cbuf_len, nal_len); | |
2197 if (nal_len < 0) | |
2198 goto fail; | |
2199 cidx += nal_len; | |
2200 } else { | |
2201 if (aidx >= abuf_len) | |
2202 break; | |
2203 nal_buf = abuf + aidx; | |
2204 nal_len = find_nal_end(nal_buf, abuf_len - aidx); | |
2205 // printf("aidx=%d/%d nal_len=%d\n", aidx, abuf_len, nal_len); | |
2206 if (nal_len < 0) | |
2207 goto fail; | |
2208 aidx += nal_len; | |
2209 } | |
2210 start = 3 + (nal_buf[2] == 0); | |
2211 if (!is_alpha) { | |
2212 int nut; | |
2213 /* add SEI NAL for the frame duration (animation case) */ | |
2214 nut = (nal_buf[start] >> 1) & 0x3f; | |
2215 if ((nut <= 9 || (nut >= 16 && nut <= 21)) && | |
2216 start + 2 < nal_len && (nal_buf[start + 2] & 0x80)) { | |
2217 int frame_ticks; | |
2218 assert(frame_num < frame_count); | |
2219 frame_ticks = frame_duration_tab[frame_num]; | |
2220 if (frame_ticks > 1) { | |
2221 add_frame_duration_sei(out_buf, frame_ticks); | |
2222 } | |
2223 frame_num++; | |
2224 } | |
2225 } | |
2226 if (first_nal) { | |
2227 /* skip first start code */ | |
2228 l = start; | |
2229 } else { | |
2230 l = 0; | |
2231 } | |
2232 if (dyn_buf_resize(out_buf, out_buf->len + nal_len - l) < 0) | |
2233 goto fail; | |
2234 // printf("add nal len=%d\n", nal_len - l); | |
2235 memcpy(out_buf->buf + out_buf->len, nal_buf + l, nal_len - l); | |
2236 if (is_alpha) { | |
2237 /* set nul_layer_id of alpha to '1' */ | |
2238 out_buf->buf[out_buf->len + (start - l) + 1] |= 1 << 3; | |
2239 } | |
2240 out_buf->len += nal_len - l; | |
2241 | |
2242 if (abuf) { | |
2243 is_alpha ^= 1; | |
2244 } | |
2245 first_nal = 0; | |
2246 } | |
2247 *pout_buf = out_buf->buf; | |
2248 return out_buf->len; | |
2249 fail: | |
2250 free(out_buf->buf); | |
2251 return -1; | |
2252 } | |
2253 | |
2254 typedef enum { | |
2255 #if defined(USE_X265) | |
2256 HEVC_ENCODER_X265, | |
2257 #endif | |
2258 #if defined(USE_JCTVC) | |
2259 HEVC_ENCODER_JCTVC, | |
2260 #endif | |
2261 | |
2262 HEVC_ENCODER_COUNT, | |
2263 } HEVCEncoderEnum; | |
2264 | |
2265 static char *hevc_encoder_name[HEVC_ENCODER_COUNT] = { | |
2266 #if defined(USE_X265) | |
2267 "x265", | |
2268 #endif | |
2269 #if defined(USE_JCTVC) | |
2270 "jctvc", | |
2271 #endif | |
2272 }; | |
2273 | |
2274 static HEVCEncoder *hevc_encoder_tab[HEVC_ENCODER_COUNT] = { | |
2275 #if defined(USE_X265) | |
2276 &x265_hevc_encoder, | |
2277 #endif | |
2278 #if defined(USE_JCTVC) | |
2279 &jctvc_encoder, | |
2280 #endif | |
2281 }; | |
2282 | |
2283 #define IMAGE_HEADER_MAGIC 0x425047fb | |
2284 | |
2285 #define DEFAULT_OUTFILENAME "out.bpg" | |
2286 #define DEFAULT_QP 29 | |
2287 #define DEFAULT_BIT_DEPTH 8 | |
2288 | |
2289 #ifdef RExt__HIGH_BIT_DEPTH_SUPPORT | |
2290 #define BIT_DEPTH_MAX 14 | |
2291 #else | |
2292 #define BIT_DEPTH_MAX 12 | |
2293 #endif | |
2294 #define DEFAULT_COMPRESS_LEVEL 8 | |
2295 | |
2296 | |
2297 typedef struct BPGEncoderContext BPGEncoderContext; | |
2298 | |
2299 typedef struct BPGEncoderParameters { | |
2300 int qp; /* 0 ... 51 */ | |
2301 int alpha_qp; /* -1 ... 51. -1 means same as qp */ | |
2302 int lossless; /* true if lossless compression (qp and alpha_qp are | |
2303 ignored) */ | |
2304 BPGImageFormatEnum preferred_chroma_format; | |
2305 int sei_decoded_picture_hash; /* 0, 1 */ | |
2306 int compress_level; /* 1 ... 9 */ | |
2307 int verbose; | |
2308 HEVCEncoderEnum encoder_type; | |
2309 int animated; /* 0 ... 1: if true, encode as animated image */ | |
2310 uint16_t loop_count; /* animations: number of loops. 0=infinite */ | |
2311 /* animations: the frame delay is a multiple of | |
2312 frame_delay_num/frame_delay_den seconds */ | |
2313 uint16_t frame_delay_num; | |
2314 uint16_t frame_delay_den; | |
2315 } BPGEncoderParameters; | |
2316 | |
2317 typedef int BPGEncoderWriteFunc(void *opaque, const uint8_t *buf, int buf_len); | |
2318 | |
2319 struct BPGEncoderContext { | |
2320 BPGEncoderParameters params; | |
2321 BPGMetaData *first_md; | |
2322 HEVCEncoder *encoder; | |
2323 int frame_count; | |
2324 HEVCEncoderContext *enc_ctx; | |
2325 HEVCEncoderContext *alpha_enc_ctx; | |
2326 int frame_ticks; | |
2327 uint16_t *frame_duration_tab; | |
2328 int frame_duration_tab_size; | |
2329 }; | |
2330 | |
2331 void *mallocz(size_t size) | |
2332 { | |
2333 void *ptr; | |
2334 ptr = malloc(size); | |
2335 if (!ptr) | |
2336 return NULL; | |
2337 memset(ptr, 0, size); | |
2338 return ptr; | |
2339 } | |
2340 | |
2341 BPGEncoderParameters *bpg_encoder_param_alloc(void) | |
2342 { | |
2343 BPGEncoderParameters *p; | |
2344 p = mallocz(sizeof(BPGEncoderParameters)); | |
2345 if (!p) | |
2346 return NULL; | |
2347 p->qp = DEFAULT_QP; | |
2348 p->alpha_qp = -1; | |
2349 p->preferred_chroma_format = BPG_FORMAT_420; | |
2350 p->compress_level = DEFAULT_COMPRESS_LEVEL; | |
2351 p->frame_delay_num = 1; | |
2352 p->frame_delay_den = 25; | |
2353 p->loop_count = 0; | |
2354 return p; | |
2355 } | |
2356 | |
2357 void bpg_encoder_param_free(BPGEncoderParameters *p) | |
2358 { | |
2359 free(p); | |
2360 } | |
2361 | |
2362 BPGEncoderContext *bpg_encoder_open(BPGEncoderParameters *p) | |
2363 { | |
2364 BPGEncoderContext *s; | |
2365 | |
2366 s = mallocz(sizeof(BPGEncoderContext)); | |
2367 if (!s) | |
2368 return NULL; | |
2369 s->params = *p; | |
2370 s->encoder = hevc_encoder_tab[s->params.encoder_type]; | |
2371 s->frame_ticks = 1; | |
2372 return s; | |
2373 } | |
2374 | |
2375 void bpg_encoder_set_extension_data(BPGEncoderContext *s, | |
2376 BPGMetaData *md) | |
2377 { | |
2378 s->first_md = md; | |
2379 } | |
2380 | |
2381 static int bpg_encoder_encode_trailer(BPGEncoderContext *s, | |
2382 BPGEncoderWriteFunc *write_func, | |
2383 void *opaque) | |
2384 { | |
2385 uint8_t *out_buf, *alpha_buf, *hevc_buf; | |
2386 int out_buf_len, alpha_buf_len, hevc_buf_len; | |
2387 | |
2388 out_buf_len = s->encoder->close(s->enc_ctx, &out_buf); | |
2389 if (out_buf_len < 0) { | |
2390 fprintf(stderr, "Error while encoding picture\n"); | |
2391 exit(1); | |
2392 } | |
2393 s->enc_ctx = NULL; | |
2394 | |
2395 alpha_buf = NULL; | |
2396 alpha_buf_len = 0; | |
2397 if (s->alpha_enc_ctx) { | |
2398 alpha_buf_len = s->encoder->close(s->alpha_enc_ctx, &alpha_buf); | |
2399 if (alpha_buf_len < 0) { | |
2400 fprintf(stderr, "Error while encoding picture (alpha plane)\n"); | |
2401 exit(1); | |
2402 } | |
2403 s->alpha_enc_ctx = NULL; | |
2404 } | |
2405 | |
2406 hevc_buf = NULL; | |
2407 hevc_buf_len = build_modified_hevc(&hevc_buf, out_buf, out_buf_len, | |
2408 alpha_buf, alpha_buf_len, | |
2409 s->frame_duration_tab, s->frame_count); | |
2410 if (hevc_buf_len < 0) { | |
2411 fprintf(stderr, "Error while creating HEVC data\n"); | |
2412 exit(1); | |
2413 } | |
2414 free(out_buf); | |
2415 free(alpha_buf); | |
2416 | |
2417 if (write_func(opaque, hevc_buf, hevc_buf_len) != hevc_buf_len) { | |
2418 fprintf(stderr, "Error while writing HEVC data\n"); | |
2419 exit(1); | |
2420 } | |
2421 free(hevc_buf); | |
2422 return 0; | |
2423 } | |
2424 | |
2425 int bpg_encoder_set_frame_duration(BPGEncoderContext *s, int frame_ticks) | |
2426 { | |
2427 if (frame_ticks >= 1 && frame_ticks <= 65535) { | |
2428 s->frame_ticks = frame_ticks; | |
2429 return 0; | |
2430 } else { | |
2431 return -1; | |
2432 } | |
2433 } | |
2434 | |
2435 /* Warning: currently 'img' is modified. When encoding animations, img | |
2436 = NULL indicates the end of the stream. */ | |
2437 int bpg_encoder_encode(BPGEncoderContext *s, Image *img, | |
2438 BPGEncoderWriteFunc *write_func, | |
2439 void *opaque) | |
2440 { | |
2441 const BPGEncoderParameters *p = &s->params; | |
2442 Image *img_alpha; | |
2443 HEVCEncodeParams ep_s, *ep = &ep_s; | |
2444 uint8_t *extension_buf; | |
2445 int extension_buf_len; | |
2446 int cb_size, width, height; | |
2447 | |
2448 if (p->animated && !img) { | |
2449 return bpg_encoder_encode_trailer(s, write_func, opaque); | |
2450 } | |
2451 | |
2452 /* extract the alpha plane */ | |
2453 if (img->has_alpha) { | |
2454 int c_idx; | |
2455 | |
2456 img_alpha = malloc(sizeof(Image)); | |
2457 memset(img_alpha, 0, sizeof(*img_alpha)); | |
2458 if (img->format == BPG_FORMAT_GRAY) | |
2459 c_idx = 1; | |
2460 else | |
2461 c_idx = 3; | |
2462 | |
2463 img_alpha->w = img->w; | |
2464 img_alpha->h = img->h; | |
2465 img_alpha->format = BPG_FORMAT_GRAY; | |
2466 img_alpha->has_alpha = 0; | |
2467 img_alpha->color_space = BPG_CS_YCbCr; | |
2468 img_alpha->bit_depth = img->bit_depth; | |
2469 img_alpha->pixel_shift = img->pixel_shift; | |
2470 img_alpha->data[0] = img->data[c_idx]; | |
2471 img_alpha->linesize[0] = img->linesize[c_idx]; | |
2472 | |
2473 img->data[c_idx] = NULL; | |
2474 img->has_alpha = 0; | |
2475 } else { | |
2476 img_alpha = NULL; | |
2477 } | |
2478 | |
2479 if (img->format == BPG_FORMAT_444 && img->color_space != BPG_CS_RGB) { | |
2480 if (p->preferred_chroma_format == BPG_FORMAT_420 || | |
2481 p->preferred_chroma_format == BPG_FORMAT_420_VIDEO) { | |
2482 int c_h_phase = (p->preferred_chroma_format == BPG_FORMAT_420); | |
2483 if (image_ycc444_to_ycc420(img, c_h_phase) != 0) | |
2484 goto error_convert; | |
2485 } else if (p->preferred_chroma_format == BPG_FORMAT_422 || | |
2486 p->preferred_chroma_format == BPG_FORMAT_422_VIDEO) { | |
2487 int c_h_phase = (p->preferred_chroma_format == BPG_FORMAT_422); | |
2488 if (image_ycc444_to_ycc422(img, c_h_phase) != 0) { | |
2489 error_convert: | |
2490 fprintf(stderr, "Cannot convert image\n"); | |
2491 exit(1); | |
2492 } | |
2493 } | |
2494 } | |
2495 | |
2496 cb_size = 8; /* XXX: should make it configurable. We assume the | |
2497 HEVC encoder uses the same value */ | |
2498 width = img->w; | |
2499 height = img->h; | |
2500 image_pad(img, cb_size); | |
2501 if (img_alpha) | |
2502 image_pad(img_alpha, cb_size); | |
2503 | |
2504 /* convert to the allocated pixel width to 8 bit if needed by the | |
2505 HEVC encoder */ | |
2506 if (img->bit_depth == 8) { | |
2507 image_convert16to8(img); | |
2508 if (img_alpha) | |
2509 image_convert16to8(img_alpha); | |
2510 } | |
2511 | |
2512 if (s->frame_count == 0) { | |
2513 memset(ep, 0, sizeof(*ep)); | |
2514 ep->qp = p->qp; | |
2515 ep->width = img->w; | |
2516 ep->height = img->h; | |
2517 ep->chroma_format = img->format; | |
2518 ep->bit_depth = img->bit_depth; | |
2519 ep->intra_only = !p->animated; | |
2520 ep->lossless = p->lossless; | |
2521 ep->sei_decoded_picture_hash = p->sei_decoded_picture_hash; | |
2522 ep->compress_level = p->compress_level; | |
2523 ep->verbose = p->verbose; | |
2524 | |
2525 s->enc_ctx = s->encoder->open(ep); | |
2526 if (!s->enc_ctx) { | |
2527 fprintf(stderr, "Error while opening encoder\n"); | |
2528 exit(1); | |
2529 } | |
2530 | |
2531 if (img_alpha) { | |
2532 if (p->alpha_qp < 0) | |
2533 ep->qp = p->qp; | |
2534 else | |
2535 ep->qp = p->alpha_qp; | |
2536 ep->chroma_format = 0; | |
2537 | |
2538 s->alpha_enc_ctx = s->encoder->open(ep); | |
2539 if (!s->alpha_enc_ctx) { | |
2540 fprintf(stderr, "Error while opening alpha encoder\n"); | |
2541 exit(1); | |
2542 } | |
2543 } | |
2544 | |
2545 /* prepare the extension data */ | |
2546 if (p->animated) { | |
2547 BPGMetaData *md; | |
2548 uint8_t buf[15], *q; | |
2549 | |
2550 md = bpg_md_alloc(BPG_EXTENSION_TAG_ANIM_CONTROL); | |
2551 q = buf; | |
2552 put_ue(&q, p->loop_count); | |
2553 put_ue(&q, p->frame_delay_num); | |
2554 put_ue(&q, p->frame_delay_den); | |
2555 md->buf_len = q - buf; | |
2556 md->buf = malloc(md->buf_len); | |
2557 memcpy(md->buf, buf, md->buf_len); | |
2558 md->next = s->first_md; | |
2559 s->first_md = md; | |
2560 } | |
2561 | |
2562 extension_buf = NULL; | |
2563 extension_buf_len = 0; | |
2564 if (s->first_md) { | |
2565 BPGMetaData *md1; | |
2566 int max_len; | |
2567 uint8_t *q; | |
2568 | |
2569 max_len = 0; | |
2570 for(md1 = s->first_md; md1 != NULL; md1 = md1->next) { | |
2571 max_len += md1->buf_len + 5 * 2; | |
2572 } | |
2573 extension_buf = malloc(max_len); | |
2574 q = extension_buf; | |
2575 for(md1 = s->first_md; md1 != NULL; md1 = md1->next) { | |
2576 put_ue(&q, md1->tag); | |
2577 put_ue(&q, md1->buf_len); | |
2578 memcpy(q, md1->buf, md1->buf_len); | |
2579 q += md1->buf_len; | |
2580 } | |
2581 extension_buf_len = q - extension_buf; | |
2582 | |
2583 bpg_md_free(s->first_md); | |
2584 s->first_md = NULL; | |
2585 } | |
2586 | |
2587 { | |
2588 uint8_t img_header[128], *q; | |
2589 int v, has_alpha, has_extension, alpha2_flag, alpha1_flag, format; | |
2590 | |
2591 has_alpha = (img_alpha != NULL); | |
2592 has_extension = (extension_buf_len > 0); | |
2593 | |
2594 | |
2595 if (has_alpha) { | |
2596 if (img->has_w_plane) { | |
2597 alpha1_flag = 0; | |
2598 alpha2_flag = 1; | |
2599 } else { | |
2600 alpha1_flag = 1; | |
2601 alpha2_flag = img->premultiplied_alpha; | |
2602 } | |
2603 } else { | |
2604 alpha1_flag = 0; | |
2605 alpha2_flag = 0; | |
2606 } | |
2607 | |
2608 q = img_header; | |
2609 *q++ = (IMAGE_HEADER_MAGIC >> 24) & 0xff; | |
2610 *q++ = (IMAGE_HEADER_MAGIC >> 16) & 0xff; | |
2611 *q++ = (IMAGE_HEADER_MAGIC >> 8) & 0xff; | |
2612 *q++ = (IMAGE_HEADER_MAGIC >> 0) & 0xff; | |
2613 | |
2614 if (img->c_h_phase == 0 && img->format == BPG_FORMAT_420) | |
2615 format = BPG_FORMAT_420_VIDEO; | |
2616 else if (img->c_h_phase == 0 && img->format == BPG_FORMAT_422) | |
2617 format = BPG_FORMAT_422_VIDEO; | |
2618 else | |
2619 format = img->format; | |
2620 v = (format << 5) | (alpha1_flag << 4) | (img->bit_depth - 8); | |
2621 *q++ = v; | |
2622 v = (img->color_space << 4) | (has_extension << 3) | | |
2623 (alpha2_flag << 2) | (img->limited_range << 1) | | |
2624 p->animated; | |
2625 *q++ = v; | |
2626 put_ue(&q, width); | |
2627 put_ue(&q, height); | |
2628 | |
2629 put_ue(&q, 0); /* zero length means up to the end of the file */ | |
2630 if (has_extension) { | |
2631 put_ue(&q, extension_buf_len); /* extension data length */ | |
2632 } | |
2633 | |
2634 write_func(opaque, img_header, q - img_header); | |
2635 | |
2636 if (has_extension) { | |
2637 if (write_func(opaque, extension_buf, extension_buf_len) != extension_buf_len) { | |
2638 fprintf(stderr, "Error while writing extension data\n"); | |
2639 exit(1); | |
2640 } | |
2641 free(extension_buf); | |
2642 } | |
2643 } | |
2644 } | |
2645 | |
2646 /* store the frame duration */ | |
2647 if ((s->frame_count + 1) > s->frame_duration_tab_size) { | |
2648 s->frame_duration_tab_size = (s->frame_duration_tab_size * 3) / 2; | |
2649 if (s->frame_duration_tab_size < (s->frame_count + 1)) | |
2650 s->frame_duration_tab_size = (s->frame_count + 1); | |
2651 s->frame_duration_tab = realloc(s->frame_duration_tab, | |
2652 sizeof(s->frame_duration_tab) * s->frame_duration_tab_size); | |
2653 } | |
2654 s->frame_duration_tab[s->frame_count] = s->frame_ticks; | |
2655 | |
2656 s->encoder->encode(s->enc_ctx, img); | |
2657 | |
2658 if (img_alpha) { | |
2659 s->encoder->encode(s->alpha_enc_ctx, img_alpha); | |
2660 image_free(img_alpha); | |
2661 } | |
2662 | |
2663 s->frame_count++; | |
2664 | |
2665 if (!p->animated) | |
2666 bpg_encoder_encode_trailer(s, write_func, opaque); | |
2667 | |
2668 return 0; | |
2669 } | |
2670 | |
2671 void bpg_encoder_close(BPGEncoderContext *s) | |
2672 { | |
2673 free(s->frame_duration_tab); | |
2674 bpg_md_free(s->first_md); | |
2675 free(s); | |
2676 } | |
2677 | |
2678 static int my_write_func(void *opaque, const uint8_t *buf, int buf_len) | |
2679 { | |
2680 FILE *f = opaque; | |
2681 return fwrite(buf, 1, buf_len, f); | |
2682 } | |
2683 | |
2684 static int get_filename_num(char *buf, int buf_size, const char *str, int n) | |
2685 { | |
2686 const char *p, *r; | |
2687 char *q; | |
2688 int l, c; | |
2689 | |
2690 q = buf; | |
2691 p = str; | |
2692 for(;;) { | |
2693 c = *p++; | |
2694 if (c == '\0') | |
2695 break; | |
2696 if (c == '%') { | |
2697 r = p - 1; | |
2698 l = 0; | |
2699 for(;;) { | |
2700 c = *p; | |
2701 if (c < '0' || c > '9') | |
2702 break; | |
2703 l = l * 10 + (c - '0'); | |
2704 p++; | |
2705 } | |
2706 c = *p++; | |
2707 if (c == '%') { | |
2708 goto add_char; | |
2709 } else if (c != 'd') { | |
2710 return -1; | |
2711 } | |
2712 snprintf(q, buf + buf_size - q, "%0*u", l, n); | |
2713 q += strlen(q); | |
2714 | |
2715 } else { | |
2716 add_char: | |
2717 if ((q - buf) < buf_size - 1) | |
2718 *q++ = c; | |
2719 } | |
2720 } | |
2721 *q = '\0'; | |
2722 return 0; | |
2723 } | |
2724 | |
2725 void help(int is_full) | |
2726 { | |
2727 char hevc_encoders[128]; | |
2728 int i; | |
2729 | |
2730 hevc_encoders[0] = '\0'; | |
2731 for(i = 0; i < HEVC_ENCODER_COUNT; i++) { | |
2732 if (i != 0) | |
2733 strcat(hevc_encoders, " "); | |
2734 strcat(hevc_encoders, hevc_encoder_name[i]); | |
2735 } | |
2736 | |
2737 printf("BPG Image Encoder version " CONFIG_BPG_VERSION "\n" | |
2738 "usage: bpgenc [options] infile.[jpg|png]\n" | |
2739 "\n" | |
2740 "Main options:\n" | |
2741 "-h show the full help (including the advanced options)\n" | |
2742 "-o outfile set output filename (default = %s)\n" | |
2743 "-q qp set quantizer parameter (smaller gives better quality,\n" | |
2744 " range: 0-51, default = %d)\n" | |
2745 "-f cfmt set the preferred chroma format (420, 422, 444,\n" | |
2746 " default=420)\n" | |
2747 "-c color_space set the preferred color space (ycbcr, rgb, ycgco,\n" | |
2748 " ycbcr_bt709, ycbcr_bt2020, default=ycbcr)\n" | |
2749 "-b bit_depth set the bit depth (8 to %d, default = %d)\n" | |
2750 "-lossless enable lossless mode\n" | |
2751 "-e encoder select the HEVC encoder (%s, default = %s)\n" | |
2752 "-m level select the compression level (1=fast, 9=slow, default = %d)\n" | |
2753 "\n" | |
2754 "Animation options:\n" | |
2755 "-a generate animations from a sequence of images. Use %%d or\n" | |
2756 " %%Nd (N = number of digits) in the filename to specify the\n" | |
2757 " image index, starting from 0 or 1.\n" | |
2758 "-fps N set the frame rate (default = 25)\n" | |
2759 "-loop N set the number of times the animation is played. 0 means\n" | |
2760 " infinite (default = 0)\n" | |
2761 "-delayfile file text file containing one number per image giving the\n" | |
2762 " display delay per image in centiseconds.\n" | |
2763 , DEFAULT_OUTFILENAME, DEFAULT_QP, BIT_DEPTH_MAX, DEFAULT_BIT_DEPTH, | |
2764 hevc_encoders, hevc_encoder_name[0], DEFAULT_COMPRESS_LEVEL); | |
2765 | |
2766 if (is_full) { | |
2767 printf("\nAdvanced options:\n" | |
2768 "-alphaq set quantizer parameter for the alpha channel (default = same as -q value)\n" | |
2769 "-premul store the color with premultiplied alpha\n" | |
2770 "-limitedrange encode the color data with the limited range of video\n" | |
2771 "-hash include MD5 hash in HEVC bitstream\n" | |
2772 "-keepmetadata keep the metadata (from JPEG: EXIF, ICC profile, XMP, from PNG: ICC profile)\n" | |
2773 "-v show debug messages\n" | |
2774 ); | |
2775 } | |
2776 | |
2777 exit(1); | |
2778 } | |
2779 | |
2780 struct option long_opts[] = { | |
2781 { "hash", no_argument }, | |
2782 { "keepmetadata", no_argument }, | |
2783 { "alphaq", required_argument }, | |
2784 { "lossless", no_argument }, | |
2785 { "limitedrange", no_argument }, | |
2786 { "premul", no_argument }, | |
2787 { "loop", required_argument }, | |
2788 { "fps", required_argument }, | |
2789 { "delayfile", required_argument }, | |
2790 { NULL }, | |
2791 }; | |
2792 | |
2793 int main(int argc, char **argv) | |
2794 { | |
2795 const char *infilename, *outfilename, *frame_delay_file; | |
2796 Image *img; | |
2797 FILE *f; | |
2798 int c, option_index; | |
2799 int keep_metadata; | |
2800 int bit_depth, i, limited_range, premultiplied_alpha; | |
2801 BPGColorSpaceEnum color_space; | |
2802 BPGMetaData *md; | |
2803 BPGEncoderContext *enc_ctx; | |
2804 BPGEncoderParameters *p; | |
2805 | |
2806 p = bpg_encoder_param_alloc(); | |
2807 | |
2808 outfilename = DEFAULT_OUTFILENAME; | |
2809 color_space = BPG_CS_YCbCr; | |
2810 keep_metadata = 0; | |
2811 bit_depth = DEFAULT_BIT_DEPTH; | |
2812 limited_range = 0; | |
2813 premultiplied_alpha = 0; | |
2814 frame_delay_file = NULL; | |
2815 | |
2816 for(;;) { | |
2817 c = getopt_long_only(argc, argv, "q:o:hf:c:vm:b:e:a", long_opts, &option_index); | |
2818 if (c == -1) | |
2819 break; | |
2820 switch(c) { | |
2821 case 0: | |
2822 switch(option_index) { | |
2823 case 0: | |
2824 p->sei_decoded_picture_hash = 1; | |
2825 break; | |
2826 case 1: | |
2827 keep_metadata = 1; | |
2828 break; | |
2829 case 2: | |
2830 p->alpha_qp = atoi(optarg); | |
2831 if (p->alpha_qp < 0 || p->alpha_qp > 51) { | |
2832 fprintf(stderr, "alpha_qp must be between 0 and 51\n"); | |
2833 exit(1); | |
2834 } | |
2835 break; | |
2836 case 3: | |
2837 p->lossless = 1; | |
2838 color_space = BPG_CS_RGB; | |
2839 p->preferred_chroma_format = BPG_FORMAT_444; | |
2840 bit_depth = 8; | |
2841 limited_range = 0; | |
2842 break; | |
2843 case 4: | |
2844 limited_range = 1; | |
2845 break; | |
2846 case 5: | |
2847 premultiplied_alpha = 1; | |
2848 break; | |
2849 case 6: | |
2850 p->loop_count = strtoul(optarg, NULL, 0); | |
2851 break; | |
2852 case 7: | |
2853 p->frame_delay_num = 1; | |
2854 p->frame_delay_den = strtoul(optarg, NULL, 0); | |
2855 if (p->frame_delay_den == 0) { | |
2856 fprintf(stderr, "invalid frame rate\n"); | |
2857 exit(1); | |
2858 } | |
2859 break; | |
2860 case 8: | |
2861 frame_delay_file = optarg; | |
2862 break; | |
2863 default: | |
2864 goto show_help; | |
2865 } | |
2866 break; | |
2867 case 'h': | |
2868 show_help: | |
2869 help(1); | |
2870 break; | |
2871 case 'q': | |
2872 p->qp = atoi(optarg); | |
2873 if (p->qp < 0 || p->qp > 51) { | |
2874 fprintf(stderr, "qp must be between 0 and 51\n"); | |
2875 exit(1); | |
2876 } | |
2877 break; | |
2878 case 'o': | |
2879 outfilename = optarg; | |
2880 break; | |
2881 case 'f': | |
2882 if (!strcmp(optarg, "420")) { | |
2883 p->preferred_chroma_format = BPG_FORMAT_420; | |
2884 } else if (!strcmp(optarg, "422")) { | |
2885 p->preferred_chroma_format = BPG_FORMAT_422; | |
2886 } else if (!strcmp(optarg, "444")) { | |
2887 p->preferred_chroma_format = BPG_FORMAT_444; | |
2888 } else if (!strcmp(optarg, "422_video")) { | |
2889 p->preferred_chroma_format = BPG_FORMAT_422_VIDEO; | |
2890 } else if (!strcmp(optarg, "420_video")) { | |
2891 p->preferred_chroma_format = BPG_FORMAT_420_VIDEO; | |
2892 } else { | |
2893 fprintf(stderr, "Invalid chroma format\n"); | |
2894 exit(1); | |
2895 } | |
2896 break; | |
2897 case 'c': | |
2898 if (!strcmp(optarg, "ycbcr")) { | |
2899 color_space = BPG_CS_YCbCr; | |
2900 } else if (!strcmp(optarg, "rgb")) { | |
2901 color_space = BPG_CS_RGB; | |
2902 p->preferred_chroma_format = BPG_FORMAT_444; | |
2903 } else if (!strcmp(optarg, "ycgco")) { | |
2904 color_space = BPG_CS_YCgCo; | |
2905 } else if (!strcmp(optarg, "ycbcr_bt709")) { | |
2906 color_space = BPG_CS_YCbCr_BT709; | |
2907 } else if (!strcmp(optarg, "ycbcr_bt2020")) { | |
2908 color_space = BPG_CS_YCbCr_BT2020; | |
2909 } else { | |
2910 fprintf(stderr, "Invalid color space format\n"); | |
2911 exit(1); | |
2912 } | |
2913 break; | |
2914 case 'm': | |
2915 p->compress_level = atoi(optarg); | |
2916 if (p->compress_level < 1) | |
2917 p->compress_level = 1; | |
2918 else if (p->compress_level > 9) | |
2919 p->compress_level = 9; | |
2920 break; | |
2921 case 'b': | |
2922 bit_depth = atoi(optarg); | |
2923 if (bit_depth < 8 || bit_depth > BIT_DEPTH_MAX) { | |
2924 fprintf(stderr, "Invalid bit depth (range: 8 to %d)\n", | |
2925 BIT_DEPTH_MAX); | |
2926 exit(1); | |
2927 } | |
2928 break; | |
2929 case 'v': | |
2930 p->verbose++; | |
2931 break; | |
2932 case 'e': | |
2933 for(i = 0; i < HEVC_ENCODER_COUNT; i++) { | |
2934 if (!strcmp(optarg, hevc_encoder_name[i])) | |
2935 break; | |
2936 } | |
2937 if (i == HEVC_ENCODER_COUNT) { | |
2938 fprintf(stderr, "Unsupported encoder. Available ones are:"); | |
2939 for(i = 0; i < HEVC_ENCODER_COUNT; i++) { | |
2940 fprintf(stderr, " %s", hevc_encoder_name[i]); | |
2941 } | |
2942 fprintf(stderr, "\n"); | |
2943 exit(1); | |
2944 } | |
2945 p->encoder_type = i; | |
2946 break; | |
2947 case 'a': | |
2948 p->animated = 1; | |
2949 break; | |
2950 default: | |
2951 exit(1); | |
2952 } | |
2953 } | |
2954 | |
2955 if (optind >= argc) | |
2956 help(0); | |
2957 infilename = argv[optind]; | |
2958 | |
2959 f = fopen(outfilename, "wb"); | |
2960 if (!f) { | |
2961 perror(outfilename); | |
2962 exit(1); | |
2963 } | |
2964 | |
2965 enc_ctx = bpg_encoder_open(p); | |
2966 if (!enc_ctx) { | |
2967 fprintf(stderr, "Could not open BPG encoder\n"); | |
2968 exit(1); | |
2969 } | |
2970 | |
2971 if (p->animated) { | |
2972 int frame_num, first_frame, frame_ticks; | |
2973 char filename[1024]; | |
2974 FILE *f1; | |
2975 | |
2976 if (frame_delay_file) { | |
2977 f1 = fopen(frame_delay_file, "r"); | |
2978 if (!f1) { | |
2979 fprintf(stderr, "Could not open '%s'\n", frame_delay_file); | |
2980 exit(1); | |
2981 } | |
2982 } else { | |
2983 f1 = NULL; | |
2984 } | |
2985 | |
2986 first_frame = 1; | |
2987 for(frame_num = 0; ; frame_num++) { | |
2988 if (get_filename_num(filename, sizeof(filename), infilename, frame_num) < 0) { | |
2989 fprintf(stderr, "Invalid filename syntax: '%s'\n", infilename); | |
2990 exit(1); | |
2991 } | |
2992 img = load_image(&md, filename, color_space, bit_depth, limited_range, | |
2993 premultiplied_alpha); | |
2994 if (!img) { | |
2995 if (frame_num == 0) | |
2996 continue; /* accept to start at 0 or 1 */ | |
2997 if (first_frame) { | |
2998 fprintf(stderr, "Could not read '%s'\n", filename); | |
2999 exit(1); | |
3000 } else { | |
3001 break; | |
3002 } | |
3003 } | |
3004 frame_ticks = 1; | |
3005 if (f1) { | |
3006 float fdelay; | |
3007 if (fscanf(f1, "%f", &fdelay) == 1) { | |
3008 frame_ticks = lrint(fdelay * p->frame_delay_den / (p->frame_delay_num * 100)); | |
3009 if (frame_ticks < 1) | |
3010 frame_ticks = 1; | |
3011 } | |
3012 } | |
3013 | |
3014 if (p->verbose) | |
3015 printf("Encoding '%s' ticks=%d\n", filename, frame_ticks); | |
3016 | |
3017 if (keep_metadata && first_frame) { | |
3018 bpg_encoder_set_extension_data(enc_ctx, md); | |
3019 } else { | |
3020 bpg_md_free(md); | |
3021 } | |
3022 bpg_encoder_set_frame_duration(enc_ctx, frame_ticks); | |
3023 bpg_encoder_encode(enc_ctx, img, my_write_func, f); | |
3024 image_free(img); | |
3025 | |
3026 first_frame = 0; | |
3027 } | |
3028 if (f1) | |
3029 fclose(f1); | |
3030 /* end of stream */ | |
3031 bpg_encoder_encode(enc_ctx, NULL, my_write_func, f); | |
3032 } else { | |
3033 img = load_image(&md, infilename, color_space, bit_depth, limited_range, | |
3034 premultiplied_alpha); | |
3035 if (!img) { | |
3036 fprintf(stderr, "Could not read '%s'\n", infilename); | |
3037 exit(1); | |
3038 } | |
3039 | |
3040 if (!keep_metadata && md) { | |
3041 bpg_md_free(md); | |
3042 md = NULL; | |
3043 } | |
3044 | |
3045 bpg_encoder_set_extension_data(enc_ctx, md); | |
3046 | |
3047 bpg_encoder_encode(enc_ctx, img, my_write_func, f); | |
3048 image_free(img); | |
3049 } | |
3050 | |
3051 fclose(f); | |
3052 | |
3053 bpg_encoder_close(enc_ctx); | |
3054 | |
3055 bpg_encoder_param_free(p); | |
3056 | |
3057 return 0; | |
3058 } |