Mercurial > hg > forks > libbpg
annotate bpgenc.c @ 12:9b00d4206f99
Implement support for one type of grayscale TIFF.
| author | Matti Hamalainen <ccr@tnsp.org> |
|---|---|
| date | Thu, 08 Dec 2016 22:41:01 +0200 |
| parents | e70eb0e6acd5 |
| children | df9f2fee4285 |
| 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 | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
80 static void convert_init(ColorConvertState *s, int in_bit_depth, |
| 0 | 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) { | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
95 mult_y = (double)(219 << (out_bit_depth - 8)) * (1 << c_shift) / |
| 0 | 96 (double)in_pixel_max; |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
97 mult_c = (double)(224 << (out_bit_depth - 8)) * (1 << c_shift) / |
| 0 | 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; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
108 |
| 0 | 109 case BPG_CS_YCbCr_BT709: |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
110 k_r = 0.2126; |
| 0 | 111 k_b = 0.0722; |
| 112 goto convert_ycc; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
113 |
| 0 | 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); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
127 |
| 0 | 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); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
186 cb_ptr[i] = clamp_pix(((c3 * r + c4 * g + c5 * b + |
| 0 | 187 rnd) >> shift) + center, pixel_max); |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
188 cr_ptr[i] = clamp_pix(((c6 * r + c7 * g + c8 * b + |
| 0 | 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); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
237 cb_ptr[i] = clamp_pix(((t1 - t2 + rnd) >> shift) + center, |
| 0 | 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 | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
278 static void rgb48_to_ycc(ColorConvertState *s, |
| 0 | 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); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
306 cb_ptr[i] = clamp_pix(((c3 * r + c4 * g + c5 * b + |
| 0 | 307 rnd) >> shift) + center, pixel_max); |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
308 cr_ptr[i] = clamp_pix(((c6 * r + c7 * g + c8 * b + |
| 0 | 309 rnd) >> shift) + center, pixel_max); |
| 310 src += incr; | |
| 311 } | |
| 312 } | |
| 313 | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
314 static void rgb48_to_ycgco(ColorConvertState *s, |
| 0 | 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); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
336 cb_ptr[i] = clamp_pix(((t1 - t2 + rnd) >> shift) + center, |
| 0 | 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 */ | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
345 static void gray16_to_gray(ColorConvertState *s, |
| 0 | 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 | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
360 static void luma16_to_gray(ColorConvertState *s, |
| 0 | 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 | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
375 static void rgb48_to_rgb(ColorConvertState *s, |
| 0 | 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 | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
386 typedef void RGBConvertFunc(ColorConvertState *s, |
| 0 | 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 }; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
406 |
| 0 | 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++) { | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
464 dst[i] = clamp_pix(((src[-7] + src[7]) * DP0C7 + |
|
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
465 (src[-5] + src[5]) * DP0C5 + |
|
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
466 (src[-3] + src[3]) * DP0C3 + |
|
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
467 (src[-1] + src[1]) * DP0C1 + |
| 0 | 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++) { | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
481 dst[i] = ((src[-7] + src[7]) * DP0C7 + |
|
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
482 (src[-5] + src[5]) * DP0C5 + |
|
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
483 (src[-3] + src[3]) * DP0C3 + |
|
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
484 (src[-1] + src[1]) * DP0C1 + |
| 0 | 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++) { | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
498 dst[i] = clamp_pix(((src[-4] + src[5]) * DP1C4 + |
|
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
499 (src[-3] + src[4]) * DP1C3 + |
|
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
500 (src[-2] + src[3]) * DP1C2 + |
|
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
501 (src[-1] + src[2]) * DP1C1 + |
| 0 | 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++) { | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
515 dst[i] = ((src[-4] + src[5]) * DP1C4 + |
|
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
516 (src[-3] + src[4]) * DP1C3 + |
|
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
517 (src[-2] + src[3]) * DP1C2 + |
|
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
518 (src[-1] + src[2]) * DP1C1 + |
| 0 | 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 | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
529 if (phase == 0) |
| 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 | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
556 if (phase == 0) |
| 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); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
572 |
| 0 | 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]; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
601 |
| 0 | 602 shift = 21 - bit_depth; |
| 603 offset = 1 << (shift - 1); | |
| 604 pixel_max = (1 << bit_depth) - 1; | |
| 605 for(i = 0; i < n; i++) { | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
606 dst[i] = clamp_pix(((srcm4[i] + src5[i]) * DP1C4 + |
|
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
607 (srcm3[i] + src4[i]) * DP1C3 + |
|
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
608 (srcm2[i] + src3[i]) * DP1C2 + |
|
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
609 (srcm1[i] + src2[i]) * DP1C1 + |
| 0 | 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, | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
616 uint8_t *src, int src_linesize, |
| 0 | 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; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
622 |
| 0 | 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)); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
695 |
| 0 | 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); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
716 |
| 0 | 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); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
805 |
| 0 | 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++) { | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
813 if (img->format == BPG_FORMAT_420 && |
| 0 | 814 (c_idx == 1 || c_idx == 2)) { |
| 815 h_shift = 1; | |
| 816 v_shift = 1; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
817 } else if (img->format == BPG_FORMAT_422 && |
| 0 | 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 | |
|
9
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
912 void insert_md(BPGMetaData **list, BPGMetaData *node) |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
913 { |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
914 if (*list == NULL) |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
915 { |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
916 *list = node; |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
917 node->next = NULL; |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
918 } |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
919 else |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
920 { |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
921 node->next = *list; |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
922 *list = node; |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
923 } |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
924 } |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
925 |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
926 int add_md_contents(BPGMetaData **list, const int tag, const size_t len, const void *buf) |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
927 { |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
928 BPGMetaData *md; |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
929 |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
930 if ((md = bpg_md_alloc(tag)) == NULL) |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
931 goto err; |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
932 |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
933 if ((md->buf = malloc(len)) == NULL) |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
934 goto err; |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
935 |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
936 md->buf_len = len; |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
937 memcpy(md->buf, buf, len); |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
938 insert_md(list, md); |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
939 return 0; |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
940 |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
941 err: |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
942 bpg_md_free(md); |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
943 return -1; |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
944 } |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
945 |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
946 |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
947 #define TIFF_DEBUG |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
948 |
|
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
949 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
|
950 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
|
951 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
|
952 { |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
953 TIFF *tif; |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
954 Image *img = NULL; |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
955 BPGImageFormatEnum img_format; |
|
9
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
956 uint8_t *buf = NULL, *tmp_buf = NULL; |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
957 uint32_t img_width, img_height, tmp_len; |
|
10
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
958 int img_alpha, err_spp = 0; |
|
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
959 size_t img_linesize; |
|
11
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
960 uint16_t img_depth, img_spp, img_pconfig, img_pmetric, |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
961 *buf2 = NULL, *img_colormap_r, *img_colormap_g, *img_colormap_b; |
|
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
962 ColorConvertState cvt; |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
963 RGBConvertFunc *convert_func; |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
964 |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
965 // 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
|
966 if ((tif = TIFFOpen(filename, "rb")) == NULL) |
| 7 | 967 { |
| 968 fprintf(stderr, "Could not open TIFF image '%s'\n", filename); | |
|
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
969 return NULL; |
| 7 | 970 } |
|
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
971 |
| 7 | 972 if (!TIFFGetField(tif, TIFFTAG_IMAGEWIDTH, &img_width) || |
| 973 !TIFFGetField(tif, TIFFTAG_IMAGELENGTH, &img_height) || | |
| 974 !TIFFGetField(tif, TIFFTAG_SAMPLESPERPIXEL, &img_spp) || | |
| 975 !TIFFGetField(tif, TIFFTAG_BITSPERSAMPLE, &img_depth)) | |
| 976 { | |
| 977 fprintf(stderr, "TIFF file lacks basic fields!\n"); | |
| 978 goto err; | |
| 979 } | |
| 980 | |
|
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
981 TIFFGetField(tif, TIFFTAG_PLANARCONFIG, &img_pconfig); |
|
10
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
982 TIFFGetField(tif, TIFFTAG_PHOTOMETRIC, &img_pmetric); |
|
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
983 |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
984 // Check basics |
|
10
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
985 switch (img_pmetric) |
|
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
986 { |
|
12
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
987 case PHOTOMETRIC_MINISBLACK: |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
988 img_format = BPG_FORMAT_GRAY; |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
989 color_space = BPG_CS_YCbCr; |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
990 if (img_spp != 1) |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
991 err_spp = 1; |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
992 break; |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
993 |
|
10
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
994 case PHOTOMETRIC_RGB: |
|
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
995 img_format = BPG_FORMAT_444; |
|
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
996 if (img_spp < 3) |
|
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
997 err_spp = 1; |
|
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
998 break; |
|
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
999 |
|
11
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1000 case PHOTOMETRIC_PALETTE: |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1001 img_format = BPG_FORMAT_444; |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1002 if (img_spp != 1) |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1003 err_spp = 1; |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1004 |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1005 if (!TIFFGetField(tif, TIFFTAG_COLORMAP, &img_colormap_r, &img_colormap_g, &img_colormap_b)) |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1006 { |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1007 fprintf(stderr, "TIFF with palette is missing colormap entry!\n"); |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1008 goto err; |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1009 } |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1010 break; |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1011 |
|
10
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
1012 default: |
|
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
1013 fprintf(stderr, "TIFF file has unsupported photometric interpretation %d!\n", img_pmetric); |
|
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
1014 goto err; |
|
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
1015 } |
|
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
1016 |
|
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
1017 if (err_spp) |
|
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
1018 { |
|
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
1019 fprintf(stderr, "TIFF file has unsupported samples per pixel %d for format %d.\n", img_spp, img_pmetric); |
|
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
1020 goto err; |
|
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
1021 } |
|
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
1022 |
|
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
1023 if (img_depth != 8 && img_depth != 16) |
|
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
1024 { |
|
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
1025 fprintf(stderr, "TIFF file has unsupported depth %d (not 8 or 16)!\n", img_depth); |
|
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1026 goto err; |
|
10
d04d9e3a77d0
Add checking TIFF photometric interpretation tag data.
Matti Hamalainen <ccr@tnsp.org>
parents:
9
diff
changeset
|
1027 } |
|
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1028 |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1029 img_linesize = TIFFScanlineSize(tif); |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1030 |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1031 if (img_spp > 3) |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1032 { |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1033 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
|
1034 |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1035 TIFFGetField(tif, TIFFTAG_EXTRASAMPLES, &img_esmp, &img_esmp_types); |
|
6
690bf78c1ce9
Initial implementation of alpha channel support for TIFF images.
Matti Hamalainen <ccr@tnsp.org>
parents:
5
diff
changeset
|
1036 img_alpha = img_esmp == 1 && img_esmp_types[0] == EXTRASAMPLE_ASSOCALPHA; |
|
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1037 } |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1038 else |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1039 img_alpha = 0; |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1040 |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1041 // Allocate temporary image space |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1042 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
|
1043 if (buf == NULL) |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1044 goto err; |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1045 |
|
11
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1046 if (img_pmetric == PHOTOMETRIC_PALETTE) |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1047 { |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1048 // For indexed palette images we need a secondary RGB conversion buffer |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1049 buf2 = (uint16_t *) malloc(img_width * 3 * sizeof(uint16_t)); |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1050 if (buf2 == NULL) |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1051 goto err; |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1052 |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1053 // Force 16bit depth and samples per pixel 3 as the |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1054 // TIFF colormap/palette has 16bit R, G, B entries. |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1055 img_spp = 3; |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1056 img_depth = 16; |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1057 } |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1058 |
|
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1059 // Allocate target image |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1060 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
|
1061 if (img == NULL) |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1062 goto err; |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1063 |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1064 img->limited_range = limited_range; |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1065 img->premultiplied_alpha = premultiplied_alpha; |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1066 |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1067 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
|
1068 |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1069 int idx = img_depth == 16; |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1070 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
|
1071 |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1072 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
|
1073 { |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1074 // Read TIFF image to raster |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1075 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
|
1076 goto err; |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1077 |
|
11
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1078 // Convert the scanline data |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1079 switch (img_pmetric) |
|
6
690bf78c1ce9
Initial implementation of alpha channel support for TIFF images.
Matti Hamalainen <ccr@tnsp.org>
parents:
5
diff
changeset
|
1080 { |
|
11
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1081 case PHOTOMETRIC_PALETTE: |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1082 for (int x = 0, offs = 0; x < img->w; x++) |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1083 { |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1084 buf2[offs++] = img_colormap_r[buf[x]]; |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1085 buf2[offs++] = img_colormap_g[buf[x]]; |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1086 buf2[offs++] = img_colormap_b[buf[x]]; |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1087 } |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1088 convert_func(&cvt, |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1089 (PIXEL *)(img->data[0] + y * img->linesize[0]), |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1090 (PIXEL *)(img->data[1] + y * img->linesize[1]), |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1091 (PIXEL *)(img->data[2] + y * img->linesize[2]), |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1092 buf2, img->w, img_spp); |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1093 break; |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1094 |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1095 case PHOTOMETRIC_RGB: |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1096 convert_func(&cvt, |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1097 (PIXEL *)(img->data[0] + y * img->linesize[0]), |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1098 (PIXEL *)(img->data[1] + y * img->linesize[1]), |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1099 (PIXEL *)(img->data[2] + y * img->linesize[2]), |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1100 buf, img->w, img_spp); |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1101 |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1102 if (img_alpha) |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1103 { |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1104 if (idx) |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1105 { |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1106 gray16_to_gray(&cvt, |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1107 (PIXEL *)(img->data[3] + y * img->linesize[3]), |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1108 ((uint16_t *) buf) + 3, img->w, 4); |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1109 } |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1110 else |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1111 { |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1112 gray8_to_gray(&cvt, |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1113 (PIXEL *)(img->data[3] + y * img->linesize[3]), |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1114 ((uint8_t *) buf) + 3, img->w, 4); |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1115 } |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1116 } |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1117 break; |
|
12
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1118 |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1119 case PHOTOMETRIC_MINISBLACK: |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1120 if (img_depth == 16) |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1121 { |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1122 luma16_to_gray(&cvt, |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1123 (PIXEL *)(img->data[0] + y * img->linesize[0]), |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1124 (uint16_t *) buf, img->w, img_spp); |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1125 |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1126 if (img_alpha) |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1127 { |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1128 gray16_to_gray(&cvt, |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1129 (PIXEL *)(img->data[1] + y * img->linesize[1]), |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1130 ((uint16_t *) buf) + 1, img->w, 2); |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1131 } |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1132 } |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1133 else |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1134 { |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1135 luma8_to_gray(&cvt, |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1136 (PIXEL *)(img->data[0] + y * img->linesize[0]), |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1137 (uint8_t *) buf, img->w, img_spp); |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1138 |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1139 if (img_alpha) |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1140 { |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1141 gray8_to_gray(&cvt, |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1142 (PIXEL *)(img->data[1] + y * img->linesize[1]), |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1143 ((uint8_t *) buf) + 1, img->w, 2); |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1144 } |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1145 } |
|
9b00d4206f99
Implement support for one type of grayscale TIFF.
Matti Hamalainen <ccr@tnsp.org>
parents:
11
diff
changeset
|
1146 break; |
|
6
690bf78c1ce9
Initial implementation of alpha channel support for TIFF images.
Matti Hamalainen <ccr@tnsp.org>
parents:
5
diff
changeset
|
1147 } |
|
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1148 } |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1149 |
|
9
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
1150 // Get meta data |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
1151 if (TIFFGetField(tif, TIFFTAG_ICCPROFILE, &tmp_len, &tmp_buf)) |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
1152 add_md_contents(pmd, BPG_EXTENSION_TAG_ICCP, tmp_len, tmp_buf); |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
1153 |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
1154 if (TIFFGetField(tif, TIFFTAG_XMLPACKET, &tmp_len, &tmp_buf)) |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
1155 add_md_contents(pmd, BPG_EXTENSION_TAG_XMP, tmp_len, tmp_buf); |
|
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1156 |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1157 err: |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1158 if (buf != NULL) |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1159 _TIFFfree(buf); |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1160 |
|
11
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1161 if (buf2 != NULL) |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1162 free(buf2); |
|
e70eb0e6acd5
Add support for indexed palette TIFFs, at least some of them.
Matti Hamalainen <ccr@tnsp.org>
parents:
10
diff
changeset
|
1163 |
|
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1164 TIFFClose(tif); |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1165 return img; |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1166 } |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1167 |
| 0 | 1168 Image *read_png(BPGMetaData **pmd, |
| 1169 FILE *f, BPGColorSpaceEnum color_space, int out_bit_depth, | |
| 1170 int limited_range, int premultiplied_alpha) | |
| 1171 { | |
| 1172 png_structp png_ptr; | |
| 1173 png_infop info_ptr; | |
| 1174 int bit_depth, color_type; | |
| 1175 Image *img; | |
| 1176 uint8_t **rows; | |
| 1177 int y, has_alpha, linesize, bpp; | |
| 1178 BPGImageFormatEnum format; | |
| 1179 ColorConvertState cvt_s, *cvt = &cvt_s; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1180 |
| 0 | 1181 png_ptr = png_create_read_struct(PNG_LIBPNG_VER_STRING, |
| 1182 NULL, NULL, NULL); | |
| 1183 if (png_ptr == NULL) { | |
| 1184 return NULL; | |
| 1185 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1186 |
| 0 | 1187 info_ptr = png_create_info_struct(png_ptr); |
| 1188 if (info_ptr == NULL) { | |
| 1189 png_destroy_read_struct(&png_ptr, NULL, NULL); | |
| 1190 return NULL; | |
| 1191 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1192 |
| 0 | 1193 if (setjmp(png_jmpbuf(png_ptr))) { |
| 1194 png_destroy_read_struct(&png_ptr, &info_ptr, NULL); | |
| 1195 return NULL; | |
| 1196 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1197 |
| 0 | 1198 png_init_io(png_ptr, f); |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1199 |
| 0 | 1200 png_read_info(png_ptr, info_ptr); |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1201 |
| 0 | 1202 bit_depth = png_get_bit_depth(png_ptr, info_ptr); |
| 1203 color_type = png_get_color_type(png_ptr, info_ptr); | |
| 1204 | |
| 1205 switch (color_type) { | |
| 1206 case PNG_COLOR_TYPE_PALETTE: | |
| 1207 png_set_palette_to_rgb(png_ptr); | |
| 1208 bit_depth = 8; | |
| 1209 break; | |
| 1210 case PNG_COLOR_TYPE_GRAY: | |
| 1211 case PNG_COLOR_TYPE_GRAY_ALPHA: | |
| 1212 if (bit_depth < 8) { | |
| 1213 png_set_expand_gray_1_2_4_to_8(png_ptr); | |
| 1214 bit_depth = 8; | |
| 1215 } | |
| 1216 break; | |
| 1217 } | |
| 1218 assert(bit_depth == 8 || bit_depth == 16); | |
| 1219 | |
| 1220 #if __BYTE_ORDER__ != __ORDER_BIG_ENDIAN__ | |
| 1221 if (bit_depth == 16) { | |
| 1222 png_set_swap(png_ptr); | |
| 1223 } | |
| 1224 #endif | |
| 1225 | |
| 1226 if (color_type == PNG_COLOR_TYPE_GRAY || | |
| 1227 color_type == PNG_COLOR_TYPE_GRAY_ALPHA) { | |
| 1228 format = BPG_FORMAT_GRAY; | |
| 1229 color_space = BPG_CS_YCbCr; | |
| 1230 } else { | |
| 1231 format = BPG_FORMAT_444; | |
| 1232 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1233 |
| 0 | 1234 has_alpha = (color_type == PNG_COLOR_TYPE_GRAY_ALPHA || |
| 1235 color_type == PNG_COLOR_TYPE_RGB_ALPHA); | |
| 1236 | |
| 1237 if (png_get_valid(png_ptr, info_ptr, PNG_INFO_tRNS)) { | |
| 1238 png_set_tRNS_to_alpha(png_ptr); | |
| 1239 has_alpha = 1; | |
| 1240 } | |
| 1241 | |
| 1242 if (premultiplied_alpha) { | |
| 1243 png_set_alpha_mode(png_ptr, PNG_ALPHA_ASSOCIATED, PNG_GAMMA_LINEAR); | |
| 1244 } | |
| 1245 | |
| 1246 img = image_alloc(png_get_image_width(png_ptr, info_ptr), | |
| 1247 png_get_image_height(png_ptr, info_ptr), | |
| 1248 format, has_alpha, color_space, | |
| 1249 out_bit_depth); | |
| 1250 img->limited_range = limited_range; | |
| 1251 img->premultiplied_alpha = premultiplied_alpha; | |
| 1252 | |
| 1253 rows = malloc(sizeof(rows[0]) * img->h); | |
| 1254 if (format == BPG_FORMAT_GRAY) | |
| 1255 bpp = (1 + has_alpha) * (bit_depth / 8); | |
| 1256 else | |
| 1257 bpp = (3 + has_alpha) * (bit_depth / 8); | |
| 1258 linesize = bpp * img->w; | |
| 1259 for (y = 0; y < img->h; y++) { | |
| 1260 rows[y] = malloc(linesize); | |
| 1261 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1262 |
| 0 | 1263 png_read_image(png_ptr, rows); |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1264 |
| 0 | 1265 convert_init(cvt, bit_depth, out_bit_depth, color_space, limited_range); |
| 1266 | |
| 1267 if (format != BPG_FORMAT_GRAY) { | |
| 1268 int idx; | |
| 1269 RGBConvertFunc *convert_func; | |
| 1270 | |
| 1271 idx = (bit_depth == 16); | |
| 1272 convert_func = rgb_to_cs[idx][color_space]; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1273 |
| 0 | 1274 for (y = 0; y < img->h; y++) { |
| 1275 convert_func(cvt, (PIXEL *)(img->data[0] + y * img->linesize[0]), | |
| 1276 (PIXEL *)(img->data[1] + y * img->linesize[1]), | |
| 1277 (PIXEL *)(img->data[2] + y * img->linesize[2]), | |
| 1278 rows[y], img->w, 3 + has_alpha); | |
| 1279 if (has_alpha) { | |
| 1280 if (idx) { | |
| 1281 gray16_to_gray(cvt, (PIXEL *)(img->data[3] + y * img->linesize[3]), | |
| 1282 (uint16_t *)rows[y] + 3, img->w, 4); | |
| 1283 } else { | |
| 1284 gray8_to_gray(cvt, (PIXEL *)(img->data[3] + y * img->linesize[3]), | |
| 1285 rows[y] + 3, img->w, 4); | |
| 1286 } | |
| 1287 } | |
| 1288 } | |
| 1289 } else { | |
| 1290 if (bit_depth == 16) { | |
| 1291 for (y = 0; y < img->h; y++) { | |
| 1292 luma16_to_gray(cvt, (PIXEL *)(img->data[0] + y * img->linesize[0]), | |
| 1293 (uint16_t *)rows[y], img->w, 1 + has_alpha); | |
| 1294 if (has_alpha) { | |
| 1295 gray16_to_gray(cvt, (PIXEL *)(img->data[1] + y * img->linesize[1]), | |
| 1296 (uint16_t *)rows[y] + 1, img->w, 2); | |
| 1297 } | |
| 1298 } | |
| 1299 } else { | |
| 1300 for (y = 0; y < img->h; y++) { | |
| 1301 luma8_to_gray(cvt, (PIXEL *)(img->data[0] + y * img->linesize[0]), | |
| 1302 rows[y], img->w, 1 + has_alpha); | |
| 1303 if (has_alpha) { | |
| 1304 gray8_to_gray(cvt, (PIXEL *)(img->data[1] + y * img->linesize[1]), | |
| 1305 rows[y] + 1, img->w, 2); | |
| 1306 } | |
| 1307 } | |
| 1308 } | |
| 1309 } | |
| 1310 | |
| 1311 for (y = 0; y < img->h; y++) { | |
| 1312 free(rows[y]); | |
| 1313 } | |
| 1314 free(rows); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1315 |
| 0 | 1316 png_read_end(png_ptr, info_ptr); |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1317 |
| 0 | 1318 /* get the ICC profile if present */ |
| 1319 { | |
| 1320 png_charp name; | |
| 1321 int comp_type; | |
| 1322 png_bytep iccp_buf; | |
| 1323 png_uint_32 iccp_buf_len; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1324 |
| 0 | 1325 if (png_get_iCCP(png_ptr, info_ptr, |
|
9
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
1326 &name, &comp_type, &iccp_buf, &iccp_buf_len) == PNG_INFO_iCCP) |
|
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
1327 add_md_contents(pmd, BPG_EXTENSION_TAG_ICCP, iccp_buf_len, iccp_buf); |
| 0 | 1328 } |
| 1329 | |
| 1330 png_destroy_read_struct(&png_ptr, &info_ptr, NULL); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1331 |
| 0 | 1332 return img; |
| 1333 } | |
| 1334 | |
| 1335 static BPGMetaData *jpeg_get_metadata(jpeg_saved_marker_ptr first_marker) | |
| 1336 { | |
| 1337 static const char app1_exif[] = "Exif"; | |
| 1338 static const char app1_xmp[] = "http://ns.adobe.com/xap/1.0/"; | |
| 1339 static const char app2_iccp[] = "ICC_PROFILE"; | |
| 1340 jpeg_saved_marker_ptr marker; | |
| 1341 BPGMetaData *md, **plast_md, *first_md; | |
| 1342 int has_exif, has_xmp, l, iccp_chunk_count, i; | |
| 1343 jpeg_saved_marker_ptr iccp_chunks[256]; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1344 |
| 0 | 1345 iccp_chunk_count = 0; |
| 1346 has_exif = 0; | |
| 1347 has_xmp = 0; | |
| 1348 first_md = NULL; | |
| 1349 plast_md = &first_md; | |
| 1350 for (marker = first_marker; marker != NULL; marker = marker->next) { | |
| 1351 #if 0 | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1352 printf("marker=APP%d len=%d\n", |
| 0 | 1353 marker->marker - JPEG_APP0, marker->data_length); |
| 1354 #endif | |
| 1355 if (!has_exif && marker->marker == JPEG_APP0 + 1 && | |
| 1356 marker->data_length > sizeof(app1_exif) && | |
| 1357 !memcmp(marker->data, app1_exif, sizeof(app1_exif))) { | |
| 1358 md = bpg_md_alloc(BPG_EXTENSION_TAG_EXIF); | |
| 1359 l = sizeof(app1_exif); | |
| 1360 md->buf_len = marker->data_length - l; | |
| 1361 md->buf = malloc(md->buf_len); | |
| 1362 memcpy(md->buf, marker->data + l, md->buf_len); | |
| 1363 *plast_md = md; | |
| 1364 plast_md = &md->next; | |
| 1365 has_exif = 1; | |
| 1366 } else if (!has_xmp && marker->marker == JPEG_APP0 + 1 && | |
| 1367 marker->data_length > sizeof(app1_xmp) && | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1368 !memcmp(marker->data, app1_xmp, sizeof(app1_xmp)) && |
| 0 | 1369 !has_xmp) { |
| 1370 md = bpg_md_alloc(BPG_EXTENSION_TAG_XMP); | |
| 1371 l = sizeof(app1_xmp); | |
| 1372 md->buf_len = marker->data_length - l; | |
| 1373 md->buf = malloc(md->buf_len); | |
| 1374 memcpy(md->buf, marker->data + l, md->buf_len); | |
| 1375 *plast_md = md; | |
| 1376 plast_md = &md->next; | |
| 1377 has_xmp = 1; | |
| 1378 } else if (marker->marker == JPEG_APP0 + 2 && | |
| 1379 marker->data_length > (sizeof(app2_iccp) + 2) && | |
| 1380 !memcmp(marker->data, app2_iccp, sizeof(app2_iccp))) { | |
| 1381 int chunk_count, chunk_index; | |
| 1382 l = sizeof(app2_iccp); | |
| 1383 chunk_index = marker->data[l]; | |
| 1384 chunk_count = marker->data[l]; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1385 if (chunk_index == 0 || chunk_count == 0) |
| 0 | 1386 continue; |
| 1387 if (iccp_chunk_count == 0) { | |
| 1388 iccp_chunk_count = chunk_count; | |
| 1389 for(i = 0; i < chunk_count; i++) { | |
| 1390 iccp_chunks[i] = NULL; | |
| 1391 } | |
| 1392 } else { | |
| 1393 if (chunk_count != iccp_chunk_count) | |
| 1394 continue; | |
| 1395 } | |
| 1396 if (chunk_index > iccp_chunk_count) | |
| 1397 continue; | |
| 1398 iccp_chunks[chunk_index - 1] = marker; | |
| 1399 } | |
| 1400 } | |
| 1401 | |
| 1402 if (iccp_chunk_count != 0) { | |
| 1403 int len, hlen, idx; | |
| 1404 /* check that no chunk are missing */ | |
| 1405 len = 0; | |
| 1406 hlen = sizeof(app2_iccp) + 2; | |
| 1407 for(i = 0; i < iccp_chunk_count; i++) { | |
| 1408 if (!iccp_chunks[i]) | |
| 1409 break; | |
| 1410 len += iccp_chunks[i]->data_length - hlen; | |
| 1411 } | |
| 1412 if (i == iccp_chunk_count) { | |
| 1413 md = bpg_md_alloc(BPG_EXTENSION_TAG_ICCP); | |
| 1414 md->buf_len = len; | |
| 1415 md->buf = malloc(md->buf_len); | |
| 1416 idx = 0; | |
| 1417 for(i = 0; i < iccp_chunk_count; i++) { | |
| 1418 l = iccp_chunks[i]->data_length - hlen; | |
| 1419 memcpy(md->buf + idx, iccp_chunks[i]->data + hlen, l); | |
| 1420 idx += l; | |
| 1421 } | |
| 1422 assert(idx == len); | |
| 1423 *plast_md = md; | |
| 1424 plast_md = &md->next; | |
| 1425 } | |
| 1426 } | |
| 1427 return first_md; | |
| 1428 } | |
| 1429 | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1430 Image *read_jpeg(BPGMetaData **pmd, FILE *f, |
| 0 | 1431 int out_bit_depth) |
| 1432 { | |
| 1433 struct jpeg_decompress_struct cinfo; | |
| 1434 struct jpeg_error_mgr jerr; | |
| 1435 int w, h, w1, i, y_h, c_h, y, v_shift, c_w, y1, idx, c_idx, h_shift; | |
| 1436 int h1, plane_idx[4], has_alpha, has_w_plane; | |
| 1437 Image *img; | |
| 1438 BPGImageFormatEnum format; | |
| 1439 BPGColorSpaceEnum color_space; | |
| 1440 ColorConvertState cvt_s, *cvt = &cvt_s; | |
| 1441 uint32_t comp_hv; | |
| 1442 | |
| 1443 cinfo.err = jpeg_std_error(&jerr); | |
| 1444 jpeg_create_decompress(&cinfo); | |
| 1445 | |
| 1446 jpeg_save_markers(&cinfo, JPEG_APP0 + 1, 65535); | |
| 1447 jpeg_save_markers(&cinfo, JPEG_APP0 + 2, 65535); | |
| 1448 | |
| 1449 jpeg_stdio_src(&cinfo, f); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1450 |
| 0 | 1451 jpeg_read_header(&cinfo, TRUE); |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1452 |
| 0 | 1453 cinfo.raw_data_out = TRUE; |
| 1454 cinfo.do_fancy_upsampling = TRUE; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1455 |
| 0 | 1456 w = cinfo.image_width; |
| 1457 h = cinfo.image_height; | |
| 1458 | |
| 1459 has_w_plane = 0; | |
| 1460 comp_hv = 0; | |
| 1461 if (cinfo.num_components < 1 || cinfo.num_components > 4) | |
| 1462 goto unsupported; | |
| 1463 for(i = 0; i < cinfo.num_components; i++) { | |
| 1464 comp_hv |= cinfo.comp_info[i].h_samp_factor << (i * 8 + 4); | |
| 1465 comp_hv |= cinfo.comp_info[i].v_samp_factor << (i * 8); | |
| 1466 } | |
| 1467 switch(cinfo.jpeg_color_space) { | |
| 1468 case JCS_GRAYSCALE: | |
| 1469 if (cinfo.num_components != 1 || comp_hv != 0x11) | |
| 1470 goto unsupported; | |
| 1471 format = BPG_FORMAT_GRAY; | |
| 1472 color_space = BPG_CS_YCbCr; | |
| 1473 break; | |
| 1474 case JCS_YCbCr: | |
| 1475 if (cinfo.num_components != 3) | |
| 1476 goto unsupported; | |
| 1477 switch(comp_hv) { | |
| 1478 case 0x111111: | |
| 1479 format = BPG_FORMAT_444; | |
| 1480 break; | |
| 1481 case 0x111121: | |
| 1482 format = BPG_FORMAT_422; | |
| 1483 break; | |
| 1484 case 0x111122: | |
| 1485 format = BPG_FORMAT_420; | |
| 1486 break; | |
| 1487 default: | |
| 1488 cinfo.raw_data_out = FALSE; | |
| 1489 format = BPG_FORMAT_444; | |
| 1490 cinfo.out_color_space = JCS_YCbCr; | |
| 1491 break; | |
| 1492 } | |
| 1493 color_space = BPG_CS_YCbCr; | |
| 1494 break; | |
| 1495 case JCS_RGB: | |
| 1496 if (cinfo.num_components != 3) | |
| 1497 goto unsupported; | |
| 1498 format = BPG_FORMAT_444; | |
| 1499 color_space = BPG_CS_RGB; | |
| 1500 cinfo.raw_data_out = FALSE; | |
| 1501 cinfo.out_color_space = JCS_RGB; | |
| 1502 break; | |
| 1503 case JCS_YCCK: | |
| 1504 if (cinfo.num_components != 4) | |
| 1505 goto unsupported; | |
| 1506 switch(comp_hv) { | |
| 1507 case 0x11111111: | |
| 1508 format = BPG_FORMAT_444; | |
| 1509 color_space = BPG_CS_YCbCr; | |
| 1510 break; | |
| 1511 case 0x22111121: | |
| 1512 format = BPG_FORMAT_422; | |
| 1513 color_space = BPG_CS_YCbCr; | |
| 1514 break; | |
| 1515 case 0x22111122: | |
| 1516 format = BPG_FORMAT_420; | |
| 1517 color_space = BPG_CS_YCbCr; | |
| 1518 break; | |
| 1519 default: | |
| 1520 cinfo.raw_data_out = FALSE; | |
| 1521 format = BPG_FORMAT_444; | |
| 1522 cinfo.out_color_space = JCS_CMYK; | |
| 1523 color_space = BPG_CS_RGB; | |
| 1524 break; | |
| 1525 } | |
| 1526 has_w_plane = 1; | |
| 1527 break; | |
| 1528 case JCS_CMYK: | |
| 1529 if (cinfo.num_components != 4) | |
| 1530 goto unsupported; | |
| 1531 format = BPG_FORMAT_444; | |
| 1532 color_space = BPG_CS_RGB; | |
| 1533 has_w_plane = 1; | |
| 1534 cinfo.raw_data_out = FALSE; | |
| 1535 cinfo.out_color_space = JCS_CMYK; | |
| 1536 break; | |
| 1537 default: | |
| 1538 unsupported: | |
| 1539 fprintf(stderr, "Unsupported JPEG parameters (cs=%d n_comp=%d comp_hv=%x)\n", | |
| 1540 cinfo.jpeg_color_space, cinfo.num_components, comp_hv); | |
| 1541 img = NULL; | |
| 1542 goto the_end; | |
| 1543 } | |
| 1544 | |
| 1545 v_shift = (format == BPG_FORMAT_420); | |
| 1546 h_shift = (format == BPG_FORMAT_422 || format == BPG_FORMAT_420); | |
| 1547 has_alpha = (cinfo.num_components == 4); | |
| 1548 img = image_alloc(w, h, format, has_alpha, color_space, out_bit_depth); | |
| 1549 img->has_w_plane = has_w_plane; | |
| 1550 | |
| 1551 convert_init(cvt, 8, out_bit_depth, color_space, 0); | |
| 1552 | |
| 1553 jpeg_start_decompress(&cinfo); | |
| 1554 | |
| 1555 if (color_space == BPG_CS_RGB) { | |
| 1556 plane_idx[0] = 2; | |
| 1557 plane_idx[1] = 0; | |
| 1558 plane_idx[2] = 1; | |
| 1559 } else { | |
| 1560 plane_idx[0] = 0; | |
| 1561 plane_idx[1] = 1; | |
| 1562 plane_idx[2] = 2; | |
| 1563 } | |
| 1564 plane_idx[3] = 3; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1565 |
| 0 | 1566 if (cinfo.raw_data_out) { |
| 1567 JSAMPROW rows[4][16]; | |
| 1568 JSAMPROW *plane_pointer[4]; | |
| 1569 | |
| 1570 y_h = 8 * cinfo.max_v_samp_factor; | |
| 1571 if (cinfo.num_components == 1) { | |
| 1572 c_h = 0; | |
| 1573 c_w = 0; | |
| 1574 } else { | |
| 1575 c_h = 8; | |
| 1576 if (h_shift) | |
| 1577 c_w = (w + 1) / 2; | |
| 1578 else | |
| 1579 c_w = w; | |
| 1580 } | |
| 1581 w1 = (w + 15) & ~15; | |
| 1582 for(c_idx = 0; c_idx < cinfo.num_components; c_idx++) { | |
| 1583 if (c_idx == 1 || c_idx == 2) { | |
| 1584 h1 = c_h; | |
| 1585 } else { | |
| 1586 h1 = y_h; | |
| 1587 } | |
| 1588 for(i = 0; i < h1; i++) { | |
| 1589 rows[c_idx][i] = malloc(w1); | |
| 1590 } | |
| 1591 plane_pointer[c_idx] = rows[c_idx]; | |
| 1592 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1593 |
| 0 | 1594 while (cinfo.output_scanline < cinfo.output_height) { |
| 1595 y = cinfo.output_scanline; | |
| 1596 jpeg_read_raw_data(&cinfo, plane_pointer, y_h); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1597 |
| 0 | 1598 for(c_idx = 0; c_idx < cinfo.num_components; c_idx++) { |
| 1599 if (c_idx == 1 || c_idx == 2) { | |
| 1600 h1 = c_h; | |
| 1601 w1 = c_w; | |
| 1602 y1 = (y >> v_shift); | |
| 1603 } else { | |
| 1604 h1 = y_h; | |
| 1605 w1 = img->w; | |
| 1606 y1 = y; | |
| 1607 } | |
| 1608 idx = plane_idx[c_idx]; | |
| 1609 for(i = 0; i < h1; i++) { | |
| 1610 PIXEL *ptr; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1611 ptr = (PIXEL *)(img->data[idx] + |
| 0 | 1612 img->linesize[idx] * (y1 + i)); |
| 1613 gray8_to_gray(cvt, ptr, rows[c_idx][i], w1, 1); | |
| 1614 if (color_space == BPG_CS_YCbCr && has_w_plane) { | |
| 1615 /* negate color */ | |
| 1616 if (c_idx == 0) { | |
| 1617 gray_one_minus(cvt, ptr, w1); | |
| 1618 } else if (c_idx <= 2) { | |
| 1619 gray_neg_c(cvt, ptr, w1); | |
| 1620 } | |
| 1621 } | |
| 1622 } | |
| 1623 } | |
| 1624 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1625 |
| 0 | 1626 for(c_idx = 0; c_idx < cinfo.num_components; c_idx++) { |
| 1627 if (c_idx == 1 || c_idx == 2) { | |
| 1628 h1 = c_h; | |
| 1629 } else { | |
| 1630 h1 = y_h; | |
| 1631 } | |
| 1632 for(i = 0; i < h1; i++) { | |
| 1633 free(rows[c_idx][i]); | |
| 1634 } | |
| 1635 } | |
| 1636 } else { | |
| 1637 JSAMPROW rows[1]; | |
| 1638 uint8_t *buf; | |
| 1639 int c_count; | |
| 1640 | |
| 1641 c_count = 3 + has_w_plane; | |
| 1642 buf = malloc(c_count * w); | |
| 1643 rows[0] = buf; | |
| 1644 while (cinfo.output_scanline < cinfo.output_height) { | |
| 1645 y = cinfo.output_scanline; | |
| 1646 jpeg_read_scanlines(&cinfo, rows, 1); | |
| 1647 | |
| 1648 for(c_idx = 0; c_idx < c_count; c_idx++) { | |
| 1649 idx = plane_idx[c_idx]; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1650 gray8_to_gray(cvt, (PIXEL *)(img->data[idx] + |
| 0 | 1651 img->linesize[idx] * y), |
| 1652 buf + c_idx, w, c_count); | |
| 1653 } | |
| 1654 } | |
| 1655 free(buf); | |
| 1656 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1657 |
|
9
bfd933ef38b4
Implement some metadata support in TIFF and clean up meta data handling otherwise.
Matti Hamalainen <ccr@tnsp.org>
parents:
8
diff
changeset
|
1658 *pmd = jpeg_get_metadata(cinfo.marker_list); |
| 0 | 1659 |
| 1660 the_end: | |
| 1661 jpeg_finish_decompress(&cinfo); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1662 |
| 0 | 1663 jpeg_destroy_decompress(&cinfo); |
| 1664 return img; | |
| 1665 } | |
| 1666 | |
|
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1667 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
|
1668 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
|
1669 #define PROBE_BUF_SIZE 8 |
|
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1670 |
| 0 | 1671 Image *load_image(BPGMetaData **pmd, const char *infilename, |
| 1672 BPGColorSpaceEnum color_space, int bit_depth, | |
| 1673 int limited_range, int premultiplied_alpha) | |
| 1674 { | |
| 1675 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
|
1676 Image *img = NULL; |
|
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1677 uint8_t buf[PROBE_BUF_SIZE]; |
| 0 | 1678 |
| 1679 *pmd = NULL; | |
| 1680 | |
|
3
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1681 if ((f = fopen(infilename, "rb")) == NULL) |
| 0 | 1682 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
|
1683 |
|
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1684 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
|
1685 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
|
1686 |
|
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1687 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
|
1688 |
|
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1689 if (png_sig_cmp(buf, 0, PROBE_BUF_SIZE) == 0) |
|
8
5a1eec3c43cc
Get rid of indirect metadata variable sillyness in load_image().
Matti Hamalainen <ccr@tnsp.org>
parents:
7
diff
changeset
|
1690 img = read_png(pmd, f, color_space, bit_depth, limited_range, premultiplied_alpha); |
|
3
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1691 else |
|
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1692 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
|
1693 memcmp(buf, PROBE_TIFF_MAGIC_BE, sizeof(PROBE_TIFF_MAGIC_BE)) == 0) |
|
8
5a1eec3c43cc
Get rid of indirect metadata variable sillyness in load_image().
Matti Hamalainen <ccr@tnsp.org>
parents:
7
diff
changeset
|
1694 img = read_tiff(pmd, infilename, color_space, bit_depth, limited_range, premultiplied_alpha); |
|
4
bbd61622856a
Initial implementation of TIFF format read support via libtiff.
Matti Hamalainen <ccr@tnsp.org>
parents:
3
diff
changeset
|
1695 else |
|
8
5a1eec3c43cc
Get rid of indirect metadata variable sillyness in load_image().
Matti Hamalainen <ccr@tnsp.org>
parents:
7
diff
changeset
|
1696 img = read_jpeg(pmd, 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
|
1697 |
|
ffcd1967fb5c
Clean up load_image() probing code a bit in preparation for TIFF support.
Matti Hamalainen <ccr@tnsp.org>
parents:
0
diff
changeset
|
1698 err: |
| 0 | 1699 fclose(f); |
| 1700 return img; | |
| 1701 } | |
| 1702 | |
| 1703 void save_yuv1(Image *img, FILE *f) | |
| 1704 { | |
| 1705 int c_w, c_h, i, c_count, y; | |
| 1706 | |
| 1707 if (img->format == BPG_FORMAT_GRAY) | |
| 1708 c_count = 1; | |
| 1709 else | |
| 1710 c_count = 3; | |
| 1711 for(i = 0; i < c_count; i++) { | |
| 1712 get_plane_res(img, &c_w, &c_h, i); | |
| 1713 for(y = 0; y < c_h; y++) { | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1714 fwrite(img->data[i] + y * img->linesize[i], |
| 0 | 1715 1, c_w << img->pixel_shift, f); |
| 1716 } | |
| 1717 } | |
| 1718 } | |
| 1719 | |
| 1720 void save_yuv(Image *img, const char *filename) | |
| 1721 { | |
| 1722 FILE *f; | |
| 1723 | |
| 1724 f = fopen(filename, "wb"); | |
| 1725 if (!f) { | |
| 1726 fprintf(stderr, "Could not open %s\n", filename); | |
| 1727 exit(1); | |
| 1728 } | |
| 1729 save_yuv1(img, f); | |
| 1730 fclose(f); | |
| 1731 } | |
| 1732 | |
| 1733 /* return the position of the end of the NAL or -1 if error */ | |
| 1734 static int find_nal_end(const uint8_t *buf, int buf_len) | |
| 1735 { | |
| 1736 int idx; | |
| 1737 | |
| 1738 idx = 0; | |
| 1739 if (buf_len >= 4 && | |
| 1740 buf[0] == 0 && buf[1] == 0 && buf[2] == 0 && buf[3] == 1) { | |
| 1741 idx = 4; | |
| 1742 } else if (buf_len >= 3 && | |
| 1743 buf[0] == 0 && buf[1] == 0 && buf[2] == 1) { | |
| 1744 idx = 3; | |
| 1745 } else { | |
| 1746 return -1; | |
| 1747 } | |
| 1748 /* NAL header */ | |
| 1749 if (idx + 2 > buf_len) | |
| 1750 return -1; | |
| 1751 /* find the last byte */ | |
| 1752 for(;;) { | |
| 1753 if (idx + 2 >= buf_len) { | |
| 1754 idx = buf_len; | |
| 1755 break; | |
| 1756 } | |
| 1757 if (buf[idx] == 0 && buf[idx + 1] == 0 && buf[idx + 2] == 1) | |
| 1758 break; | |
| 1759 if (idx + 3 < buf_len && | |
| 1760 buf[idx] == 0 && buf[idx + 1] == 0 && buf[idx + 2] == 0 && buf[idx + 3] == 1) | |
| 1761 break; | |
| 1762 idx++; | |
| 1763 } | |
| 1764 return idx; | |
| 1765 } | |
| 1766 | |
| 1767 /* return the position of the end of the NAL or -1 if error */ | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1768 static int extract_nal(uint8_t **pnal_buf, int *pnal_len, |
| 0 | 1769 const uint8_t *buf, int buf_len) |
| 1770 { | |
| 1771 int idx, start, end, len; | |
| 1772 uint8_t *nal_buf; | |
| 1773 int nal_len; | |
| 1774 | |
| 1775 end = find_nal_end(buf, buf_len); | |
| 1776 if (end < 0) | |
| 1777 return -1; | |
| 1778 if (buf[2] == 1) | |
| 1779 start = 3; | |
| 1780 else | |
| 1781 start = 4; | |
| 1782 len = end - start; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1783 |
| 0 | 1784 nal_buf = malloc(len); |
| 1785 nal_len = 0; | |
| 1786 idx = start; | |
| 1787 while (idx < end) { | |
| 1788 if (idx + 2 < end && buf[idx] == 0 && buf[idx + 1] == 0 && buf[idx + 2] == 3) { | |
| 1789 nal_buf[nal_len++] = 0; | |
| 1790 nal_buf[nal_len++] = 0; | |
| 1791 idx += 3; | |
| 1792 } else { | |
| 1793 nal_buf[nal_len++] = buf[idx++]; | |
| 1794 } | |
| 1795 } | |
| 1796 while (idx < end) { | |
| 1797 nal_buf[nal_len++] = buf[idx++]; | |
| 1798 } | |
| 1799 *pnal_buf = nal_buf; | |
| 1800 *pnal_len = nal_len; | |
| 1801 return idx; | |
| 1802 } | |
| 1803 | |
| 1804 /* big endian variable length 7 bit encoding */ | |
| 1805 static void put_ue(uint8_t **pp, uint32_t v) | |
| 1806 { | |
| 1807 uint8_t *p = *pp; | |
| 1808 int i, j; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1809 |
| 0 | 1810 for(i = 1; i < 5; i++) { |
| 1811 if (v < (1 << (7 * i))) | |
| 1812 break; | |
| 1813 } | |
| 1814 for(j = i - 1; j >= 1; j--) | |
| 1815 *p++ = ((v >> (7 * j)) & 0x7f) | 0x80; | |
| 1816 *p++ = v & 0x7f; | |
| 1817 *pp = p; | |
| 1818 } | |
| 1819 | |
| 1820 typedef struct { | |
| 1821 const uint8_t *buf; | |
| 1822 int idx; | |
| 1823 int buf_len; | |
| 1824 } GetBitState; | |
| 1825 | |
| 1826 static void init_get_bits(GetBitState *s, const uint8_t *buf, int buf_len) | |
| 1827 { | |
| 1828 s->buf = buf; | |
| 1829 s->buf_len = buf_len; | |
| 1830 s->idx = 0; | |
| 1831 } | |
| 1832 | |
| 1833 static void skip_bits(GetBitState *s, int n) | |
| 1834 { | |
| 1835 s->idx += n; | |
| 1836 } | |
| 1837 | |
| 1838 /* 1 <= n <= 25. return '0' bits if past the end of the buffer. */ | |
| 1839 static uint32_t get_bits(GetBitState *s, int n) | |
| 1840 { | |
| 1841 const uint8_t *buf = s->buf; | |
| 1842 int p, i; | |
| 1843 uint32_t v; | |
| 1844 | |
| 1845 p = s->idx >> 3; | |
| 1846 if ((p + 3) < s->buf_len) { | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1847 v = (buf[p] << 24) | (buf[p + 1] << 16) | |
| 0 | 1848 (buf[p + 2] << 8) | buf[p + 3]; |
| 1849 } else { | |
| 1850 v = 0; | |
| 1851 for(i = 0; i < 3; i++) { | |
| 1852 if ((p + i) < s->buf_len) | |
| 1853 v |= buf[p + i] << (24 - i * 8); | |
| 1854 } | |
| 1855 } | |
| 1856 v = (v >> (32 - (s->idx & 7) - n)) & ((1 << n) - 1); | |
| 1857 s->idx += n; | |
| 1858 return v; | |
| 1859 } | |
| 1860 | |
| 1861 /* 1 <= n <= 32 */ | |
| 1862 static uint32_t get_bits_long(GetBitState *s, int n) | |
| 1863 { | |
| 1864 uint32_t v; | |
| 1865 | |
| 1866 if (n <= 25) { | |
| 1867 v = get_bits(s, n); | |
| 1868 } else { | |
| 1869 n -= 16; | |
| 1870 v = get_bits(s, 16) << n; | |
| 1871 v |= get_bits(s, n); | |
| 1872 } | |
| 1873 return v; | |
| 1874 } | |
| 1875 | |
| 1876 /* at most 32 bits are supported */ | |
| 1877 static uint32_t get_ue_golomb(GetBitState *s) | |
| 1878 { | |
| 1879 int i; | |
| 1880 i = 0; | |
| 1881 for(;;) { | |
| 1882 if (get_bits(s, 1)) | |
| 1883 break; | |
| 1884 i++; | |
| 1885 if (i == 32) | |
| 1886 return 0xffffffff; | |
| 1887 } | |
| 1888 if (i == 0) | |
| 1889 return 0; | |
| 1890 else | |
| 1891 return ((1 << i) | get_bits_long(s, i)) - 1; | |
| 1892 } | |
| 1893 | |
| 1894 typedef struct { | |
| 1895 uint8_t *buf; | |
| 1896 int idx; | |
| 1897 } PutBitState; | |
| 1898 | |
| 1899 static void init_put_bits(PutBitState *s, uint8_t *buf) | |
| 1900 { | |
| 1901 s->buf = buf; | |
| 1902 s->idx = 0; | |
| 1903 } | |
| 1904 | |
| 1905 static void put_bit(PutBitState *s, int bit) | |
| 1906 { | |
| 1907 s->buf[s->idx >> 3] |= bit << (7 - (s->idx & 7)); | |
| 1908 s->idx++; | |
| 1909 } | |
| 1910 | |
| 1911 static void put_bits(PutBitState *s, int n, uint32_t v) | |
| 1912 { | |
| 1913 int i; | |
| 1914 | |
| 1915 for(i = 0; i < n; i++) { | |
| 1916 put_bit(s, (v >> (n - 1 - i)) & 1); | |
| 1917 } | |
| 1918 } | |
| 1919 | |
| 1920 static void put_ue_golomb(PutBitState *s, uint32_t v) | |
| 1921 { | |
| 1922 uint32_t a; | |
| 1923 int n; | |
| 1924 | |
| 1925 v++; | |
| 1926 n = 0; | |
| 1927 a = v; | |
| 1928 while (a != 0) { | |
| 1929 a >>= 1; | |
| 1930 n++; | |
| 1931 } | |
| 1932 if (n > 1) | |
| 1933 put_bits(s, n - 1, 0); | |
| 1934 put_bits(s, n, v); | |
| 1935 } | |
| 1936 | |
| 1937 typedef struct { | |
| 1938 uint8_t *buf; | |
| 1939 int size; | |
| 1940 int len; | |
| 1941 } DynBuf; | |
| 1942 | |
| 1943 static void dyn_buf_init(DynBuf *s) | |
| 1944 { | |
| 1945 s->buf = NULL; | |
| 1946 s->size = 0; | |
| 1947 s->len = 0; | |
| 1948 } | |
| 1949 | |
| 1950 static int dyn_buf_resize(DynBuf *s, int size) | |
| 1951 { | |
| 1952 int new_size; | |
| 1953 uint8_t *new_buf; | |
| 1954 | |
| 1955 if (size <= s->size) | |
| 1956 return 0; | |
| 1957 new_size = (s->size * 3) / 2; | |
| 1958 if (new_size < size) | |
| 1959 new_size = size; | |
| 1960 new_buf = realloc(s->buf, new_size); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
1961 if (!new_buf) |
| 0 | 1962 return -1; |
| 1963 s->buf = new_buf; | |
| 1964 s->size = new_size; | |
| 1965 return 0; | |
| 1966 } | |
| 1967 | |
| 1968 /* suppress the VPS NAL and keep only the useful part of the SPS | |
| 1969 header. The decoder can rebuild a valid HEVC stream if needed. */ | |
| 1970 static int build_modified_sps(uint8_t **pout_buf, int *pout_buf_len, | |
| 1971 const uint8_t *buf, int buf_len) | |
| 1972 { | |
| 1973 int nal_unit_type, nal_len, idx, i, ret, msps_buf_len; | |
| 1974 int out_buf_len, out_buf_len_max; | |
| 1975 uint8_t *nal_buf, *msps_buf, *out_buf; | |
| 1976 GetBitState gb_s, *gb = &gb_s; | |
| 1977 PutBitState pb_s, *pb = &pb_s; | |
| 1978 uint8_t *p; | |
| 1979 | |
| 1980 idx = extract_nal(&nal_buf, &nal_len, buf, buf_len); | |
| 1981 if (idx < 0) | |
| 1982 return -1; | |
| 1983 if (nal_len < 2) { | |
| 1984 free(nal_buf); | |
| 1985 return -1; | |
| 1986 } | |
| 1987 nal_unit_type = (nal_buf[0] >> 1) & 0x3f; | |
| 1988 free(nal_buf); | |
| 1989 if (nal_unit_type != 32) { | |
| 1990 fprintf(stderr, "expecting VPS nal (%d)\n", nal_unit_type); | |
| 1991 return -1; /* expect VPS nal */ | |
| 1992 } | |
| 1993 | |
| 1994 ret = extract_nal(&nal_buf, &nal_len, buf + idx, buf_len); | |
| 1995 if (ret < 0) | |
| 1996 return -1; | |
| 1997 idx += ret; | |
| 1998 if (nal_len < 2) | |
| 1999 return -1; | |
| 2000 nal_unit_type = (nal_buf[0] >> 1) & 0x3f; | |
| 2001 if (nal_unit_type != 33) { | |
| 2002 fprintf(stderr, "expecting SPS nal (%d)\n", nal_unit_type); | |
| 2003 return -1; /* expect SPS nal */ | |
| 2004 } | |
| 2005 | |
| 2006 /* skip the initial part of the SPS up to and including | |
| 2007 log2_min_cb_size */ | |
| 2008 { | |
| 2009 int vps_id, max_sub_layers, profile_idc, sps_id; | |
| 2010 int chroma_format_idc, width, height, bit_depth_luma, bit_depth_chroma; | |
| 2011 int log2_max_poc_lsb, sublayer_ordering_info, log2_min_cb_size; | |
| 2012 int log2_diff_max_min_coding_block_size, log2_min_tb_size; | |
| 2013 int log2_diff_max_min_transform_block_size; | |
| 2014 int max_transform_hierarchy_depth_inter; | |
| 2015 int max_transform_hierarchy_depth_intra; | |
| 2016 int scaling_list_enable_flag, amp_enabled_flag, sao_enabled; | |
| 2017 int pcm_enabled_flag, nb_st_rps; | |
| 2018 int long_term_ref_pics_present_flag, sps_strong_intra_smoothing_enable_flag, vui_present; | |
| 2019 int sps_temporal_mvp_enabled_flag; | |
| 2020 int pcm_sample_bit_depth_luma_minus1; | |
| 2021 int pcm_sample_bit_depth_chroma_minus1; | |
| 2022 int log2_min_pcm_luma_coding_block_size_minus3; | |
| 2023 int log2_diff_max_min_pcm_luma_coding_block_size; | |
| 2024 int pcm_loop_filter_disabled_flag; | |
| 2025 int sps_extension_flag, sps_range_extension_flag, sps_extension_7bits; | |
| 2026 int sps_range_extension_flags; | |
| 2027 | |
| 2028 init_get_bits(gb, nal_buf, nal_len); | |
| 2029 skip_bits(gb, 16); /* nal header */ | |
| 2030 vps_id = get_bits(gb, 4); | |
| 2031 if (vps_id != 0) { | |
| 2032 fprintf(stderr, "VPS id 0 expected\n"); | |
| 2033 return -1; | |
| 2034 } | |
| 2035 max_sub_layers = get_bits(gb, 3); | |
| 2036 if (max_sub_layers != 0) { | |
| 2037 fprintf(stderr, "max_sub_layers == 0 expected\n"); | |
| 2038 return -1; | |
| 2039 } | |
| 2040 skip_bits(gb, 1); /* temporal_id_nesting_flag */ | |
| 2041 /* profile tier level */ | |
| 2042 skip_bits(gb, 2); /* profile_space */ | |
| 2043 skip_bits(gb, 1); /* tier_flag */ | |
| 2044 profile_idc = get_bits(gb, 5); | |
| 2045 for(i = 0; i < 32; i++) { | |
| 2046 skip_bits(gb, 1); /* profile_compatibility_flag */ | |
| 2047 } | |
| 2048 skip_bits(gb, 1); /* progressive_source_flag */ | |
| 2049 skip_bits(gb, 1); /* interlaced_source_flag */ | |
| 2050 skip_bits(gb, 1); /* non_packed_constraint_flag */ | |
| 2051 skip_bits(gb, 1); /* frame_only_constraint_flag */ | |
| 2052 skip_bits(gb, 44); /* XXX_reserved_zero_44 */ | |
| 2053 skip_bits(gb, 8); /* level_idc */ | |
| 2054 | |
| 2055 sps_id = get_ue_golomb(gb); | |
| 2056 if (sps_id != 0) { | |
| 2057 fprintf(stderr, "SPS id 0 expected (%d)\n", sps_id); | |
| 2058 return -1; | |
| 2059 } | |
| 2060 chroma_format_idc = get_ue_golomb(gb); | |
| 2061 if (chroma_format_idc == 3) { | |
| 2062 get_bits(gb, 1); /* separate_colour_plane_flag */ | |
| 2063 } | |
| 2064 width = get_ue_golomb(gb); | |
| 2065 height = get_ue_golomb(gb); | |
| 2066 /* pic conformance_flag */ | |
| 2067 if (get_bits(gb, 1)) { | |
| 2068 get_ue_golomb(gb); /* left_offset */ | |
| 2069 get_ue_golomb(gb); /* right_offset */ | |
| 2070 get_ue_golomb(gb); /* top_offset */ | |
| 2071 get_ue_golomb(gb); /* bottom_offset */ | |
| 2072 } | |
| 2073 bit_depth_luma = get_ue_golomb(gb) + 8; | |
| 2074 bit_depth_chroma = get_ue_golomb(gb) + 8; | |
| 2075 log2_max_poc_lsb = get_ue_golomb(gb) + 4; | |
| 2076 if (log2_max_poc_lsb != 8) { | |
| 2077 fprintf(stderr, "log2_max_poc_lsb must be 8 (%d)\n", log2_max_poc_lsb); | |
| 2078 return -1; | |
| 2079 } | |
| 2080 sublayer_ordering_info = get_bits(gb, 1); | |
| 2081 get_ue_golomb(gb); /* max_dec_pic_buffering */ | |
| 2082 get_ue_golomb(gb); /* num_reorder_pics */ | |
| 2083 get_ue_golomb(gb); /* max_latency_increase */ | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2084 |
| 0 | 2085 log2_min_cb_size = get_ue_golomb(gb) + 3; |
| 2086 log2_diff_max_min_coding_block_size = get_ue_golomb(gb); | |
| 2087 log2_min_tb_size = get_ue_golomb(gb) + 2; | |
| 2088 log2_diff_max_min_transform_block_size = get_ue_golomb(gb); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2089 |
| 0 | 2090 max_transform_hierarchy_depth_inter = get_ue_golomb(gb); |
| 2091 max_transform_hierarchy_depth_intra = get_ue_golomb(gb); | |
| 2092 if (max_transform_hierarchy_depth_inter != max_transform_hierarchy_depth_intra) { | |
| 2093 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); | |
| 2094 return -1; | |
| 2095 } | |
| 2096 | |
| 2097 scaling_list_enable_flag = get_bits(gb, 1); | |
| 2098 if (scaling_list_enable_flag != 0) { | |
| 2099 fprintf(stderr, "scaling_list_enable_flag must be 0\n"); | |
| 2100 return -1; | |
| 2101 } | |
| 2102 amp_enabled_flag = get_bits(gb, 1); | |
| 2103 if (!amp_enabled_flag) { | |
| 2104 fprintf(stderr, "amp_enabled_flag must be set\n"); | |
| 2105 return -1; | |
| 2106 } | |
| 2107 sao_enabled = get_bits(gb, 1); | |
| 2108 pcm_enabled_flag = get_bits(gb, 1); | |
| 2109 if (pcm_enabled_flag) { | |
| 2110 pcm_sample_bit_depth_luma_minus1 = get_bits(gb, 4); | |
| 2111 pcm_sample_bit_depth_chroma_minus1 = get_bits(gb, 4); | |
| 2112 log2_min_pcm_luma_coding_block_size_minus3 = get_ue_golomb(gb); | |
| 2113 log2_diff_max_min_pcm_luma_coding_block_size = get_ue_golomb(gb); | |
| 2114 pcm_loop_filter_disabled_flag = get_bits(gb, 1); | |
| 2115 } | |
| 2116 nb_st_rps = get_ue_golomb(gb); | |
| 2117 if (nb_st_rps != 0) { | |
| 2118 fprintf(stderr, "nb_st_rps must be 0 (%d)\n", nb_st_rps); | |
| 2119 return -1; | |
| 2120 } | |
| 2121 long_term_ref_pics_present_flag = get_bits(gb, 1); | |
| 2122 if (long_term_ref_pics_present_flag) { | |
| 2123 fprintf(stderr, "nlong_term_ref_pics_present_flag must be 0 (%d)\n", nb_st_rps); | |
| 2124 return -1; | |
| 2125 } | |
| 2126 sps_temporal_mvp_enabled_flag = get_bits(gb, 1); | |
| 2127 if (!sps_temporal_mvp_enabled_flag) { | |
| 2128 fprintf(stderr, "sps_temporal_mvp_enabled_flag must be set\n"); | |
| 2129 return -1; | |
| 2130 } | |
| 2131 sps_strong_intra_smoothing_enable_flag = get_bits(gb, 1); | |
| 2132 vui_present = get_bits(gb, 1); | |
| 2133 if (vui_present) { | |
| 2134 int sar_present, sar_idx, overscan_info_present_flag; | |
| 2135 int video_signal_type_present_flag, chroma_loc_info_present_flag; | |
| 2136 int default_display_window_flag, vui_timing_info_present_flag; | |
| 2137 int vui_poc_proportional_to_timing_flag; | |
| 2138 int vui_hrd_parameters_present_flag, bitstream_restriction_flag; | |
| 2139 | |
| 2140 sar_present = get_bits(gb, 1); | |
| 2141 if (sar_present) { | |
| 2142 sar_idx = get_bits(gb, 8); | |
| 2143 if (sar_idx == 255) { | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2144 skip_bits(gb, 16); /* sar_num */ |
|
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2145 skip_bits(gb, 16); /* sar_den */ |
| 0 | 2146 } |
| 2147 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2148 |
| 0 | 2149 overscan_info_present_flag = get_bits(gb, 1); |
| 2150 if (overscan_info_present_flag) { | |
| 2151 skip_bits(gb, 1); /* overscan_appropriate_flag */ | |
| 2152 } | |
| 2153 | |
| 2154 video_signal_type_present_flag = get_bits(gb, 1); | |
| 2155 if (video_signal_type_present_flag) { | |
| 2156 fprintf(stderr, "video_signal_type_present_flag must be 0\n"); | |
| 2157 return -1; | |
| 2158 } | |
| 2159 chroma_loc_info_present_flag = get_bits(gb, 1); | |
| 2160 if (chroma_loc_info_present_flag) { | |
| 2161 get_ue_golomb(gb); | |
| 2162 get_ue_golomb(gb); | |
| 2163 } | |
| 2164 skip_bits(gb, 1); /* neutra_chroma_indication_flag */ | |
| 2165 skip_bits(gb, 1); | |
| 2166 skip_bits(gb, 1); | |
| 2167 default_display_window_flag = get_bits(gb, 1); | |
| 2168 if (default_display_window_flag) { | |
| 2169 fprintf(stderr, "default_display_window_flag must be 0\n"); | |
| 2170 return -1; | |
| 2171 } | |
| 2172 vui_timing_info_present_flag = get_bits(gb, 1); | |
| 2173 if (vui_timing_info_present_flag) { | |
| 2174 skip_bits(gb, 32); | |
| 2175 skip_bits(gb, 32); | |
| 2176 vui_poc_proportional_to_timing_flag = get_bits(gb, 1); | |
| 2177 if (vui_poc_proportional_to_timing_flag) { | |
| 2178 get_ue_golomb(gb); | |
| 2179 } | |
| 2180 vui_hrd_parameters_present_flag = get_bits(gb, 1); | |
| 2181 if (vui_hrd_parameters_present_flag) { | |
| 2182 fprintf(stderr, "vui_hrd_parameters_present_flag must be 0\n"); | |
| 2183 return -1; | |
| 2184 } | |
| 2185 } | |
| 2186 bitstream_restriction_flag = get_bits(gb, 1); | |
| 2187 if (bitstream_restriction_flag) { | |
| 2188 skip_bits(gb, 1); | |
| 2189 skip_bits(gb, 1); | |
| 2190 skip_bits(gb, 1); | |
| 2191 get_ue_golomb(gb); | |
| 2192 get_ue_golomb(gb); | |
| 2193 get_ue_golomb(gb); | |
| 2194 get_ue_golomb(gb); | |
| 2195 get_ue_golomb(gb); | |
| 2196 } | |
| 2197 } | |
| 2198 sps_extension_flag = get_bits(gb, 1); | |
| 2199 sps_range_extension_flag = 0; | |
| 2200 sps_range_extension_flags = 0; | |
| 2201 if (sps_extension_flag) { | |
| 2202 sps_range_extension_flag = get_bits(gb, 1); | |
| 2203 sps_extension_7bits = get_bits(gb, 7); | |
| 2204 if (sps_extension_7bits != 0) { | |
| 2205 fprintf(stderr, "sps_extension_7bits must be 0\n"); | |
| 2206 return -1; | |
| 2207 } | |
| 2208 if (sps_range_extension_flag) { | |
| 2209 sps_range_extension_flags = get_bits(gb, 9); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2210 if (sps_range_extension_flags & ((1 << (8 - 3)) | |
|
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2211 (1 << (8 - 4)) | |
|
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2212 (1 << (8 - 6)) | |
| 0 | 2213 (1 << (8 - 8)))) { |
| 2214 fprintf(stderr, "unsupported range extensions (0x%x)\n", | |
| 2215 sps_range_extension_flags); | |
| 2216 return -1; | |
| 2217 } | |
| 2218 } | |
| 2219 } | |
| 2220 | |
| 2221 /* build the modified SPS */ | |
| 2222 msps_buf = malloc(nal_len + 32); | |
| 2223 memset(msps_buf, 0, nal_len + 16); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2224 |
| 0 | 2225 init_put_bits(pb, msps_buf); |
| 2226 put_ue_golomb(pb, log2_min_cb_size - 3); | |
| 2227 put_ue_golomb(pb, log2_diff_max_min_coding_block_size); | |
| 2228 put_ue_golomb(pb, log2_min_tb_size - 2); | |
| 2229 put_ue_golomb(pb, log2_diff_max_min_transform_block_size); | |
| 2230 put_ue_golomb(pb, max_transform_hierarchy_depth_intra); | |
| 2231 put_bits(pb, 1, sao_enabled); | |
| 2232 put_bits(pb, 1, pcm_enabled_flag); | |
| 2233 if (pcm_enabled_flag) { | |
| 2234 put_bits(pb, 4, pcm_sample_bit_depth_luma_minus1); | |
| 2235 put_bits(pb, 4, pcm_sample_bit_depth_chroma_minus1); | |
| 2236 put_ue_golomb(pb, log2_min_pcm_luma_coding_block_size_minus3); | |
| 2237 put_ue_golomb(pb, log2_diff_max_min_pcm_luma_coding_block_size); | |
| 2238 put_bits(pb, 1, pcm_loop_filter_disabled_flag); | |
| 2239 } | |
| 2240 put_bits(pb, 1, sps_strong_intra_smoothing_enable_flag); | |
| 2241 put_bits(pb, 1, sps_extension_flag); | |
| 2242 if (sps_extension_flag) { | |
| 2243 put_bits(pb, 1, sps_range_extension_flag); | |
| 2244 put_bits(pb, 7, 0); | |
| 2245 if (sps_range_extension_flag) { | |
| 2246 put_bits(pb, 9, sps_range_extension_flags); | |
| 2247 } | |
| 2248 } | |
| 2249 msps_buf_len = (pb->idx + 7) >> 3; | |
| 2250 | |
| 2251 out_buf_len_max = 5 + msps_buf_len; | |
| 2252 out_buf = malloc(out_buf_len_max); | |
| 2253 | |
| 2254 // printf("msps_n_bits=%d\n", pb->idx); | |
| 2255 p = out_buf; | |
| 2256 put_ue(&p, msps_buf_len); /* header length */ | |
| 2257 | |
| 2258 memcpy(p, msps_buf, msps_buf_len); | |
| 2259 p += msps_buf_len; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2260 |
| 0 | 2261 out_buf_len = p - out_buf; |
| 2262 free(msps_buf); | |
| 2263 free(nal_buf); | |
| 2264 } | |
| 2265 *pout_buf = out_buf; | |
| 2266 *pout_buf_len = out_buf_len; | |
| 2267 return idx; | |
| 2268 } | |
| 2269 | |
| 2270 static int add_frame_duration_sei(DynBuf *out_buf, uint16_t frame_ticks) | |
| 2271 { | |
| 2272 uint8_t nal_buf[128], *q; | |
| 2273 int nut, nal_len; | |
| 2274 | |
| 2275 q = nal_buf; | |
| 2276 *q++ = 0x00; | |
| 2277 *q++ = 0x00; | |
| 2278 *q++ = 0x01; | |
| 2279 nut = 39; /* prefix SEI NUT */ | |
| 2280 *q++ = (nut << 1); | |
| 2281 *q++ = 1; | |
| 2282 *q++ = 0xff; /* payload_type = 257 */ | |
| 2283 *q++ = 0x02; | |
| 2284 *q++ = 2; /* payload_size = 2 */ | |
| 2285 *q++ = frame_ticks >> 8; | |
| 2286 *q++ = frame_ticks; | |
| 2287 *q++ = 0x80; /* extra '1' bit and align to byte */ | |
| 2288 /* Note: the 0x00 0x00 b pattern with b <= 3 cannot happen, so no | |
| 2289 need to escape */ | |
| 2290 nal_len = q - nal_buf; | |
| 2291 if (dyn_buf_resize(out_buf, out_buf->len + nal_len) < 0) | |
| 2292 return -1; | |
| 2293 memcpy(out_buf->buf + out_buf->len, nal_buf, nal_len); | |
| 2294 out_buf->len += nal_len; | |
| 2295 return 0; | |
| 2296 } | |
| 2297 | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2298 static int build_modified_hevc(uint8_t **pout_buf, |
| 0 | 2299 const uint8_t *cbuf, int cbuf_len, |
| 2300 const uint8_t *abuf, int abuf_len, | |
| 2301 const uint16_t *frame_duration_tab, | |
| 2302 int frame_count) | |
| 2303 { | |
| 2304 DynBuf out_buf_s, *out_buf = &out_buf_s; | |
| 2305 uint8_t *msps; | |
| 2306 const uint8_t *nal_buf; | |
| 2307 int msps_len, cidx, aidx, is_alpha, nal_len, first_nal, start, l, frame_num; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2308 |
| 0 | 2309 dyn_buf_init(out_buf); |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2310 |
| 0 | 2311 /* add alpha MSPS */ |
| 2312 aidx = 0; /* avoids warning */ | |
| 2313 if (abuf) { | |
| 2314 aidx = build_modified_sps(&msps, &msps_len, abuf, abuf_len); | |
| 2315 if (aidx < 0) | |
| 2316 goto fail; | |
| 2317 if (dyn_buf_resize(out_buf, out_buf->len + msps_len) < 0) | |
| 2318 goto fail; | |
| 2319 memcpy(out_buf->buf + out_buf->len, msps, msps_len); | |
| 2320 out_buf->len += msps_len; | |
| 2321 free(msps); | |
| 2322 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2323 |
| 0 | 2324 /* add color MSPS */ |
| 2325 cidx = build_modified_sps(&msps, &msps_len, cbuf, cbuf_len); | |
| 2326 if (cidx < 0) | |
| 2327 goto fail; | |
| 2328 if (dyn_buf_resize(out_buf, out_buf->len + msps_len) < 0) | |
| 2329 goto fail; | |
| 2330 memcpy(out_buf->buf + out_buf->len, msps, msps_len); | |
| 2331 out_buf->len += msps_len; | |
| 2332 free(msps); | |
| 2333 | |
| 2334 /* add the remaining NALs, alternating between alpha (if present) | |
| 2335 and color. */ | |
| 2336 is_alpha = (abuf != NULL); | |
| 2337 first_nal = 1; | |
| 2338 frame_num = 0; | |
| 2339 for(;;) { | |
| 2340 if (!is_alpha) { | |
| 2341 if (cidx >= cbuf_len) { | |
| 2342 if (abuf) { | |
| 2343 fprintf(stderr, "Incorrect number of alpha NALs\n"); | |
| 2344 goto fail; | |
| 2345 } | |
| 2346 break; | |
| 2347 } | |
| 2348 nal_buf = cbuf + cidx; | |
| 2349 nal_len = find_nal_end(nal_buf, cbuf_len - cidx); | |
| 2350 // printf("cidx=%d/%d nal_len=%d\n", cidx, cbuf_len, nal_len); | |
| 2351 if (nal_len < 0) | |
| 2352 goto fail; | |
| 2353 cidx += nal_len; | |
| 2354 } else { | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2355 if (aidx >= abuf_len) |
| 0 | 2356 break; |
| 2357 nal_buf = abuf + aidx; | |
| 2358 nal_len = find_nal_end(nal_buf, abuf_len - aidx); | |
| 2359 // printf("aidx=%d/%d nal_len=%d\n", aidx, abuf_len, nal_len); | |
| 2360 if (nal_len < 0) | |
| 2361 goto fail; | |
| 2362 aidx += nal_len; | |
| 2363 } | |
| 2364 start = 3 + (nal_buf[2] == 0); | |
| 2365 if (!is_alpha) { | |
| 2366 int nut; | |
| 2367 /* add SEI NAL for the frame duration (animation case) */ | |
| 2368 nut = (nal_buf[start] >> 1) & 0x3f; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2369 if ((nut <= 9 || (nut >= 16 && nut <= 21)) && |
| 0 | 2370 start + 2 < nal_len && (nal_buf[start + 2] & 0x80)) { |
| 2371 int frame_ticks; | |
| 2372 assert(frame_num < frame_count); | |
| 2373 frame_ticks = frame_duration_tab[frame_num]; | |
| 2374 if (frame_ticks > 1) { | |
| 2375 add_frame_duration_sei(out_buf, frame_ticks); | |
| 2376 } | |
| 2377 frame_num++; | |
| 2378 } | |
| 2379 } | |
| 2380 if (first_nal) { | |
| 2381 /* skip first start code */ | |
| 2382 l = start; | |
| 2383 } else { | |
| 2384 l = 0; | |
| 2385 } | |
| 2386 if (dyn_buf_resize(out_buf, out_buf->len + nal_len - l) < 0) | |
| 2387 goto fail; | |
| 2388 // printf("add nal len=%d\n", nal_len - l); | |
| 2389 memcpy(out_buf->buf + out_buf->len, nal_buf + l, nal_len - l); | |
| 2390 if (is_alpha) { | |
| 2391 /* set nul_layer_id of alpha to '1' */ | |
| 2392 out_buf->buf[out_buf->len + (start - l) + 1] |= 1 << 3; | |
| 2393 } | |
| 2394 out_buf->len += nal_len - l; | |
| 2395 | |
| 2396 if (abuf) { | |
| 2397 is_alpha ^= 1; | |
| 2398 } | |
| 2399 first_nal = 0; | |
| 2400 } | |
| 2401 *pout_buf = out_buf->buf; | |
| 2402 return out_buf->len; | |
| 2403 fail: | |
| 2404 free(out_buf->buf); | |
| 2405 return -1; | |
| 2406 } | |
| 2407 | |
| 2408 typedef enum { | |
| 2409 #if defined(USE_X265) | |
| 2410 HEVC_ENCODER_X265, | |
| 2411 #endif | |
| 2412 #if defined(USE_JCTVC) | |
| 2413 HEVC_ENCODER_JCTVC, | |
| 2414 #endif | |
| 2415 | |
| 2416 HEVC_ENCODER_COUNT, | |
| 2417 } HEVCEncoderEnum; | |
| 2418 | |
| 2419 static char *hevc_encoder_name[HEVC_ENCODER_COUNT] = { | |
| 2420 #if defined(USE_X265) | |
| 2421 "x265", | |
| 2422 #endif | |
| 2423 #if defined(USE_JCTVC) | |
| 2424 "jctvc", | |
| 2425 #endif | |
| 2426 }; | |
| 2427 | |
| 2428 static HEVCEncoder *hevc_encoder_tab[HEVC_ENCODER_COUNT] = { | |
| 2429 #if defined(USE_X265) | |
| 2430 &x265_hevc_encoder, | |
| 2431 #endif | |
| 2432 #if defined(USE_JCTVC) | |
| 2433 &jctvc_encoder, | |
| 2434 #endif | |
| 2435 }; | |
| 2436 | |
| 2437 #define IMAGE_HEADER_MAGIC 0x425047fb | |
| 2438 | |
| 2439 #define DEFAULT_OUTFILENAME "out.bpg" | |
| 2440 #define DEFAULT_QP 29 | |
| 2441 #define DEFAULT_BIT_DEPTH 8 | |
| 2442 | |
| 2443 #ifdef RExt__HIGH_BIT_DEPTH_SUPPORT | |
| 2444 #define BIT_DEPTH_MAX 14 | |
| 2445 #else | |
| 2446 #define BIT_DEPTH_MAX 12 | |
| 2447 #endif | |
| 2448 #define DEFAULT_COMPRESS_LEVEL 8 | |
| 2449 | |
| 2450 | |
| 2451 typedef struct BPGEncoderContext BPGEncoderContext; | |
| 2452 | |
| 2453 typedef struct BPGEncoderParameters { | |
| 2454 int qp; /* 0 ... 51 */ | |
| 2455 int alpha_qp; /* -1 ... 51. -1 means same as qp */ | |
| 2456 int lossless; /* true if lossless compression (qp and alpha_qp are | |
| 2457 ignored) */ | |
| 2458 BPGImageFormatEnum preferred_chroma_format; | |
| 2459 int sei_decoded_picture_hash; /* 0, 1 */ | |
| 2460 int compress_level; /* 1 ... 9 */ | |
| 2461 int verbose; | |
| 2462 HEVCEncoderEnum encoder_type; | |
| 2463 int animated; /* 0 ... 1: if true, encode as animated image */ | |
| 2464 uint16_t loop_count; /* animations: number of loops. 0=infinite */ | |
| 2465 /* animations: the frame delay is a multiple of | |
| 2466 frame_delay_num/frame_delay_den seconds */ | |
| 2467 uint16_t frame_delay_num; | |
| 2468 uint16_t frame_delay_den; | |
| 2469 } BPGEncoderParameters; | |
| 2470 | |
| 2471 typedef int BPGEncoderWriteFunc(void *opaque, const uint8_t *buf, int buf_len); | |
| 2472 | |
| 2473 struct BPGEncoderContext { | |
| 2474 BPGEncoderParameters params; | |
| 2475 BPGMetaData *first_md; | |
| 2476 HEVCEncoder *encoder; | |
| 2477 int frame_count; | |
| 2478 HEVCEncoderContext *enc_ctx; | |
| 2479 HEVCEncoderContext *alpha_enc_ctx; | |
| 2480 int frame_ticks; | |
| 2481 uint16_t *frame_duration_tab; | |
| 2482 int frame_duration_tab_size; | |
| 2483 }; | |
| 2484 | |
| 2485 void *mallocz(size_t size) | |
| 2486 { | |
| 2487 void *ptr; | |
| 2488 ptr = malloc(size); | |
| 2489 if (!ptr) | |
| 2490 return NULL; | |
| 2491 memset(ptr, 0, size); | |
| 2492 return ptr; | |
| 2493 } | |
| 2494 | |
| 2495 BPGEncoderParameters *bpg_encoder_param_alloc(void) | |
| 2496 { | |
| 2497 BPGEncoderParameters *p; | |
| 2498 p = mallocz(sizeof(BPGEncoderParameters)); | |
| 2499 if (!p) | |
| 2500 return NULL; | |
| 2501 p->qp = DEFAULT_QP; | |
| 2502 p->alpha_qp = -1; | |
| 2503 p->preferred_chroma_format = BPG_FORMAT_420; | |
| 2504 p->compress_level = DEFAULT_COMPRESS_LEVEL; | |
| 2505 p->frame_delay_num = 1; | |
| 2506 p->frame_delay_den = 25; | |
| 2507 p->loop_count = 0; | |
| 2508 return p; | |
| 2509 } | |
| 2510 | |
| 2511 void bpg_encoder_param_free(BPGEncoderParameters *p) | |
| 2512 { | |
| 2513 free(p); | |
| 2514 } | |
| 2515 | |
| 2516 BPGEncoderContext *bpg_encoder_open(BPGEncoderParameters *p) | |
| 2517 { | |
| 2518 BPGEncoderContext *s; | |
| 2519 | |
| 2520 s = mallocz(sizeof(BPGEncoderContext)); | |
| 2521 if (!s) | |
| 2522 return NULL; | |
| 2523 s->params = *p; | |
| 2524 s->encoder = hevc_encoder_tab[s->params.encoder_type]; | |
| 2525 s->frame_ticks = 1; | |
| 2526 return s; | |
| 2527 } | |
| 2528 | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2529 void bpg_encoder_set_extension_data(BPGEncoderContext *s, |
| 0 | 2530 BPGMetaData *md) |
| 2531 { | |
| 2532 s->first_md = md; | |
| 2533 } | |
| 2534 | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2535 static int bpg_encoder_encode_trailer(BPGEncoderContext *s, |
| 0 | 2536 BPGEncoderWriteFunc *write_func, |
| 2537 void *opaque) | |
| 2538 { | |
| 2539 uint8_t *out_buf, *alpha_buf, *hevc_buf; | |
| 2540 int out_buf_len, alpha_buf_len, hevc_buf_len; | |
| 2541 | |
| 2542 out_buf_len = s->encoder->close(s->enc_ctx, &out_buf); | |
| 2543 if (out_buf_len < 0) { | |
| 2544 fprintf(stderr, "Error while encoding picture\n"); | |
| 2545 exit(1); | |
| 2546 } | |
| 2547 s->enc_ctx = NULL; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2548 |
| 0 | 2549 alpha_buf = NULL; |
| 2550 alpha_buf_len = 0; | |
| 2551 if (s->alpha_enc_ctx) { | |
| 2552 alpha_buf_len = s->encoder->close(s->alpha_enc_ctx, &alpha_buf); | |
| 2553 if (alpha_buf_len < 0) { | |
| 2554 fprintf(stderr, "Error while encoding picture (alpha plane)\n"); | |
| 2555 exit(1); | |
| 2556 } | |
| 2557 s->alpha_enc_ctx = NULL; | |
| 2558 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2559 |
| 0 | 2560 hevc_buf = NULL; |
| 2561 hevc_buf_len = build_modified_hevc(&hevc_buf, out_buf, out_buf_len, | |
| 2562 alpha_buf, alpha_buf_len, | |
| 2563 s->frame_duration_tab, s->frame_count); | |
| 2564 if (hevc_buf_len < 0) { | |
| 2565 fprintf(stderr, "Error while creating HEVC data\n"); | |
| 2566 exit(1); | |
| 2567 } | |
| 2568 free(out_buf); | |
| 2569 free(alpha_buf); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2570 |
| 0 | 2571 if (write_func(opaque, hevc_buf, hevc_buf_len) != hevc_buf_len) { |
| 2572 fprintf(stderr, "Error while writing HEVC data\n"); | |
| 2573 exit(1); | |
| 2574 } | |
| 2575 free(hevc_buf); | |
| 2576 return 0; | |
| 2577 } | |
| 2578 | |
| 2579 int bpg_encoder_set_frame_duration(BPGEncoderContext *s, int frame_ticks) | |
| 2580 { | |
| 2581 if (frame_ticks >= 1 && frame_ticks <= 65535) { | |
| 2582 s->frame_ticks = frame_ticks; | |
| 2583 return 0; | |
| 2584 } else { | |
| 2585 return -1; | |
| 2586 } | |
| 2587 } | |
| 2588 | |
| 2589 /* Warning: currently 'img' is modified. When encoding animations, img | |
| 2590 = NULL indicates the end of the stream. */ | |
| 2591 int bpg_encoder_encode(BPGEncoderContext *s, Image *img, | |
| 2592 BPGEncoderWriteFunc *write_func, | |
| 2593 void *opaque) | |
| 2594 { | |
| 2595 const BPGEncoderParameters *p = &s->params; | |
| 2596 Image *img_alpha; | |
| 2597 HEVCEncodeParams ep_s, *ep = &ep_s; | |
| 2598 uint8_t *extension_buf; | |
| 2599 int extension_buf_len; | |
| 2600 int cb_size, width, height; | |
| 2601 | |
| 2602 if (p->animated && !img) { | |
| 2603 return bpg_encoder_encode_trailer(s, write_func, opaque); | |
| 2604 } | |
| 2605 | |
| 2606 /* extract the alpha plane */ | |
| 2607 if (img->has_alpha) { | |
| 2608 int c_idx; | |
| 2609 | |
| 2610 img_alpha = malloc(sizeof(Image)); | |
| 2611 memset(img_alpha, 0, sizeof(*img_alpha)); | |
| 2612 if (img->format == BPG_FORMAT_GRAY) | |
| 2613 c_idx = 1; | |
| 2614 else | |
| 2615 c_idx = 3; | |
| 2616 | |
| 2617 img_alpha->w = img->w; | |
| 2618 img_alpha->h = img->h; | |
| 2619 img_alpha->format = BPG_FORMAT_GRAY; | |
| 2620 img_alpha->has_alpha = 0; | |
| 2621 img_alpha->color_space = BPG_CS_YCbCr; | |
| 2622 img_alpha->bit_depth = img->bit_depth; | |
| 2623 img_alpha->pixel_shift = img->pixel_shift; | |
| 2624 img_alpha->data[0] = img->data[c_idx]; | |
| 2625 img_alpha->linesize[0] = img->linesize[c_idx]; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2626 |
| 0 | 2627 img->data[c_idx] = NULL; |
| 2628 img->has_alpha = 0; | |
| 2629 } else { | |
| 2630 img_alpha = NULL; | |
| 2631 } | |
| 2632 | |
| 2633 if (img->format == BPG_FORMAT_444 && img->color_space != BPG_CS_RGB) { | |
| 2634 if (p->preferred_chroma_format == BPG_FORMAT_420 || | |
| 2635 p->preferred_chroma_format == BPG_FORMAT_420_VIDEO) { | |
| 2636 int c_h_phase = (p->preferred_chroma_format == BPG_FORMAT_420); | |
| 2637 if (image_ycc444_to_ycc420(img, c_h_phase) != 0) | |
| 2638 goto error_convert; | |
| 2639 } else if (p->preferred_chroma_format == BPG_FORMAT_422 || | |
| 2640 p->preferred_chroma_format == BPG_FORMAT_422_VIDEO) { | |
| 2641 int c_h_phase = (p->preferred_chroma_format == BPG_FORMAT_422); | |
| 2642 if (image_ycc444_to_ycc422(img, c_h_phase) != 0) { | |
| 2643 error_convert: | |
| 2644 fprintf(stderr, "Cannot convert image\n"); | |
| 2645 exit(1); | |
| 2646 } | |
| 2647 } | |
| 2648 } | |
| 2649 | |
| 2650 cb_size = 8; /* XXX: should make it configurable. We assume the | |
| 2651 HEVC encoder uses the same value */ | |
| 2652 width = img->w; | |
| 2653 height = img->h; | |
| 2654 image_pad(img, cb_size); | |
| 2655 if (img_alpha) | |
| 2656 image_pad(img_alpha, cb_size); | |
| 2657 | |
| 2658 /* convert to the allocated pixel width to 8 bit if needed by the | |
| 2659 HEVC encoder */ | |
| 2660 if (img->bit_depth == 8) { | |
| 2661 image_convert16to8(img); | |
| 2662 if (img_alpha) | |
| 2663 image_convert16to8(img_alpha); | |
| 2664 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2665 |
| 0 | 2666 if (s->frame_count == 0) { |
| 2667 memset(ep, 0, sizeof(*ep)); | |
| 2668 ep->qp = p->qp; | |
| 2669 ep->width = img->w; | |
| 2670 ep->height = img->h; | |
| 2671 ep->chroma_format = img->format; | |
| 2672 ep->bit_depth = img->bit_depth; | |
| 2673 ep->intra_only = !p->animated; | |
| 2674 ep->lossless = p->lossless; | |
| 2675 ep->sei_decoded_picture_hash = p->sei_decoded_picture_hash; | |
| 2676 ep->compress_level = p->compress_level; | |
| 2677 ep->verbose = p->verbose; | |
| 2678 | |
| 2679 s->enc_ctx = s->encoder->open(ep); | |
| 2680 if (!s->enc_ctx) { | |
| 2681 fprintf(stderr, "Error while opening encoder\n"); | |
| 2682 exit(1); | |
| 2683 } | |
| 2684 | |
| 2685 if (img_alpha) { | |
| 2686 if (p->alpha_qp < 0) | |
| 2687 ep->qp = p->qp; | |
| 2688 else | |
| 2689 ep->qp = p->alpha_qp; | |
| 2690 ep->chroma_format = 0; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2691 |
| 0 | 2692 s->alpha_enc_ctx = s->encoder->open(ep); |
| 2693 if (!s->alpha_enc_ctx) { | |
| 2694 fprintf(stderr, "Error while opening alpha encoder\n"); | |
| 2695 exit(1); | |
| 2696 } | |
| 2697 } | |
| 2698 | |
| 2699 /* prepare the extension data */ | |
| 2700 if (p->animated) { | |
| 2701 BPGMetaData *md; | |
| 2702 uint8_t buf[15], *q; | |
| 2703 | |
| 2704 md = bpg_md_alloc(BPG_EXTENSION_TAG_ANIM_CONTROL); | |
| 2705 q = buf; | |
| 2706 put_ue(&q, p->loop_count); | |
| 2707 put_ue(&q, p->frame_delay_num); | |
| 2708 put_ue(&q, p->frame_delay_den); | |
| 2709 md->buf_len = q - buf; | |
| 2710 md->buf = malloc(md->buf_len); | |
| 2711 memcpy(md->buf, buf, md->buf_len); | |
| 2712 md->next = s->first_md; | |
| 2713 s->first_md = md; | |
| 2714 } | |
| 2715 | |
| 2716 extension_buf = NULL; | |
| 2717 extension_buf_len = 0; | |
| 2718 if (s->first_md) { | |
| 2719 BPGMetaData *md1; | |
| 2720 int max_len; | |
| 2721 uint8_t *q; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2722 |
| 0 | 2723 max_len = 0; |
| 2724 for(md1 = s->first_md; md1 != NULL; md1 = md1->next) { | |
| 2725 max_len += md1->buf_len + 5 * 2; | |
| 2726 } | |
| 2727 extension_buf = malloc(max_len); | |
| 2728 q = extension_buf; | |
| 2729 for(md1 = s->first_md; md1 != NULL; md1 = md1->next) { | |
| 2730 put_ue(&q, md1->tag); | |
| 2731 put_ue(&q, md1->buf_len); | |
| 2732 memcpy(q, md1->buf, md1->buf_len); | |
| 2733 q += md1->buf_len; | |
| 2734 } | |
| 2735 extension_buf_len = q - extension_buf; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2736 |
| 0 | 2737 bpg_md_free(s->first_md); |
| 2738 s->first_md = NULL; | |
| 2739 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2740 |
| 0 | 2741 { |
| 2742 uint8_t img_header[128], *q; | |
| 2743 int v, has_alpha, has_extension, alpha2_flag, alpha1_flag, format; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2744 |
| 0 | 2745 has_alpha = (img_alpha != NULL); |
| 2746 has_extension = (extension_buf_len > 0); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2747 |
|
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2748 |
| 0 | 2749 if (has_alpha) { |
| 2750 if (img->has_w_plane) { | |
| 2751 alpha1_flag = 0; | |
| 2752 alpha2_flag = 1; | |
| 2753 } else { | |
| 2754 alpha1_flag = 1; | |
| 2755 alpha2_flag = img->premultiplied_alpha; | |
| 2756 } | |
| 2757 } else { | |
| 2758 alpha1_flag = 0; | |
| 2759 alpha2_flag = 0; | |
| 2760 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2761 |
| 0 | 2762 q = img_header; |
| 2763 *q++ = (IMAGE_HEADER_MAGIC >> 24) & 0xff; | |
| 2764 *q++ = (IMAGE_HEADER_MAGIC >> 16) & 0xff; | |
| 2765 *q++ = (IMAGE_HEADER_MAGIC >> 8) & 0xff; | |
| 2766 *q++ = (IMAGE_HEADER_MAGIC >> 0) & 0xff; | |
| 2767 | |
| 2768 if (img->c_h_phase == 0 && img->format == BPG_FORMAT_420) | |
| 2769 format = BPG_FORMAT_420_VIDEO; | |
| 2770 else if (img->c_h_phase == 0 && img->format == BPG_FORMAT_422) | |
| 2771 format = BPG_FORMAT_422_VIDEO; | |
| 2772 else | |
| 2773 format = img->format; | |
| 2774 v = (format << 5) | (alpha1_flag << 4) | (img->bit_depth - 8); | |
| 2775 *q++ = v; | |
| 2776 v = (img->color_space << 4) | (has_extension << 3) | | |
| 2777 (alpha2_flag << 2) | (img->limited_range << 1) | | |
| 2778 p->animated; | |
| 2779 *q++ = v; | |
| 2780 put_ue(&q, width); | |
| 2781 put_ue(&q, height); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2782 |
| 0 | 2783 put_ue(&q, 0); /* zero length means up to the end of the file */ |
| 2784 if (has_extension) { | |
| 2785 put_ue(&q, extension_buf_len); /* extension data length */ | |
| 2786 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2787 |
| 0 | 2788 write_func(opaque, img_header, q - img_header); |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2789 |
| 0 | 2790 if (has_extension) { |
| 2791 if (write_func(opaque, extension_buf, extension_buf_len) != extension_buf_len) { | |
| 2792 fprintf(stderr, "Error while writing extension data\n"); | |
| 2793 exit(1); | |
| 2794 } | |
| 2795 free(extension_buf); | |
| 2796 } | |
| 2797 } | |
| 2798 } | |
| 2799 | |
| 2800 /* store the frame duration */ | |
| 2801 if ((s->frame_count + 1) > s->frame_duration_tab_size) { | |
| 2802 s->frame_duration_tab_size = (s->frame_duration_tab_size * 3) / 2; | |
| 2803 if (s->frame_duration_tab_size < (s->frame_count + 1)) | |
| 2804 s->frame_duration_tab_size = (s->frame_count + 1); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2805 s->frame_duration_tab = realloc(s->frame_duration_tab, |
| 0 | 2806 sizeof(s->frame_duration_tab) * s->frame_duration_tab_size); |
| 2807 } | |
| 2808 s->frame_duration_tab[s->frame_count] = s->frame_ticks; | |
| 2809 | |
| 2810 s->encoder->encode(s->enc_ctx, img); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2811 |
| 0 | 2812 if (img_alpha) { |
| 2813 s->encoder->encode(s->alpha_enc_ctx, img_alpha); | |
| 2814 image_free(img_alpha); | |
| 2815 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2816 |
| 0 | 2817 s->frame_count++; |
| 2818 | |
| 2819 if (!p->animated) | |
| 2820 bpg_encoder_encode_trailer(s, write_func, opaque); | |
| 2821 | |
| 2822 return 0; | |
| 2823 } | |
| 2824 | |
| 2825 void bpg_encoder_close(BPGEncoderContext *s) | |
| 2826 { | |
| 2827 free(s->frame_duration_tab); | |
| 2828 bpg_md_free(s->first_md); | |
| 2829 free(s); | |
| 2830 } | |
| 2831 | |
| 2832 static int my_write_func(void *opaque, const uint8_t *buf, int buf_len) | |
| 2833 { | |
| 2834 FILE *f = opaque; | |
| 2835 return fwrite(buf, 1, buf_len, f); | |
| 2836 } | |
| 2837 | |
| 2838 static int get_filename_num(char *buf, int buf_size, const char *str, int n) | |
| 2839 { | |
| 2840 const char *p, *r; | |
| 2841 char *q; | |
| 2842 int l, c; | |
| 2843 | |
| 2844 q = buf; | |
| 2845 p = str; | |
| 2846 for(;;) { | |
| 2847 c = *p++; | |
| 2848 if (c == '\0') | |
| 2849 break; | |
| 2850 if (c == '%') { | |
| 2851 r = p - 1; | |
| 2852 l = 0; | |
| 2853 for(;;) { | |
| 2854 c = *p; | |
| 2855 if (c < '0' || c > '9') | |
| 2856 break; | |
| 2857 l = l * 10 + (c - '0'); | |
| 2858 p++; | |
| 2859 } | |
| 2860 c = *p++; | |
| 2861 if (c == '%') { | |
| 2862 goto add_char; | |
| 2863 } else if (c != 'd') { | |
| 2864 return -1; | |
| 2865 } | |
| 2866 snprintf(q, buf + buf_size - q, "%0*u", l, n); | |
| 2867 q += strlen(q); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2868 |
| 0 | 2869 } else { |
| 2870 add_char: | |
| 2871 if ((q - buf) < buf_size - 1) | |
| 2872 *q++ = c; | |
| 2873 } | |
| 2874 } | |
| 2875 *q = '\0'; | |
| 2876 return 0; | |
| 2877 } | |
| 2878 | |
| 2879 void help(int is_full) | |
| 2880 { | |
| 2881 char hevc_encoders[128]; | |
| 2882 int i; | |
| 2883 | |
| 2884 hevc_encoders[0] = '\0'; | |
| 2885 for(i = 0; i < HEVC_ENCODER_COUNT; i++) { | |
| 2886 if (i != 0) | |
| 2887 strcat(hevc_encoders, " "); | |
| 2888 strcat(hevc_encoders, hevc_encoder_name[i]); | |
| 2889 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2890 |
| 0 | 2891 printf("BPG Image Encoder version " CONFIG_BPG_VERSION "\n" |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2892 "usage: bpgenc [options] infile.[jpg|png|tiff]\n" |
| 0 | 2893 "\n" |
| 2894 "Main options:\n" | |
| 2895 "-h show the full help (including the advanced options)\n" | |
| 2896 "-o outfile set output filename (default = %s)\n" | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2897 "-q qp set quantizer parameter (smaller gives better quality,\n" |
| 0 | 2898 " range: 0-51, default = %d)\n" |
| 2899 "-f cfmt set the preferred chroma format (420, 422, 444,\n" | |
| 2900 " default=420)\n" | |
| 2901 "-c color_space set the preferred color space (ycbcr, rgb, ycgco,\n" | |
| 2902 " ycbcr_bt709, ycbcr_bt2020, default=ycbcr)\n" | |
| 2903 "-b bit_depth set the bit depth (8 to %d, default = %d)\n" | |
| 2904 "-lossless enable lossless mode\n" | |
| 2905 "-e encoder select the HEVC encoder (%s, default = %s)\n" | |
| 2906 "-m level select the compression level (1=fast, 9=slow, default = %d)\n" | |
| 2907 "\n" | |
| 2908 "Animation options:\n" | |
| 2909 "-a generate animations from a sequence of images. Use %%d or\n" | |
| 2910 " %%Nd (N = number of digits) in the filename to specify the\n" | |
| 2911 " image index, starting from 0 or 1.\n" | |
| 2912 "-fps N set the frame rate (default = 25)\n" | |
| 2913 "-loop N set the number of times the animation is played. 0 means\n" | |
| 2914 " infinite (default = 0)\n" | |
| 2915 "-delayfile file text file containing one number per image giving the\n" | |
| 2916 " display delay per image in centiseconds.\n" | |
| 2917 , DEFAULT_OUTFILENAME, DEFAULT_QP, BIT_DEPTH_MAX, DEFAULT_BIT_DEPTH, | |
| 2918 hevc_encoders, hevc_encoder_name[0], DEFAULT_COMPRESS_LEVEL); | |
| 2919 | |
| 2920 if (is_full) { | |
| 2921 printf("\nAdvanced options:\n" | |
| 2922 "-alphaq set quantizer parameter for the alpha channel (default = same as -q value)\n" | |
| 2923 "-premul store the color with premultiplied alpha\n" | |
| 2924 "-limitedrange encode the color data with the limited range of video\n" | |
| 2925 "-hash include MD5 hash in HEVC bitstream\n" | |
| 2926 "-keepmetadata keep the metadata (from JPEG: EXIF, ICC profile, XMP, from PNG: ICC profile)\n" | |
| 2927 "-v show debug messages\n" | |
| 2928 ); | |
| 2929 } | |
| 2930 | |
| 2931 exit(1); | |
| 2932 } | |
| 2933 | |
| 2934 struct option long_opts[] = { | |
| 2935 { "hash", no_argument }, | |
| 2936 { "keepmetadata", no_argument }, | |
| 2937 { "alphaq", required_argument }, | |
| 2938 { "lossless", no_argument }, | |
| 2939 { "limitedrange", no_argument }, | |
| 2940 { "premul", no_argument }, | |
| 2941 { "loop", required_argument }, | |
| 2942 { "fps", required_argument }, | |
| 2943 { "delayfile", required_argument }, | |
| 2944 { NULL }, | |
| 2945 }; | |
| 2946 | |
| 2947 int main(int argc, char **argv) | |
| 2948 { | |
| 2949 const char *infilename, *outfilename, *frame_delay_file; | |
| 2950 Image *img; | |
| 2951 FILE *f; | |
| 2952 int c, option_index; | |
| 2953 int keep_metadata; | |
| 2954 int bit_depth, i, limited_range, premultiplied_alpha; | |
| 2955 BPGColorSpaceEnum color_space; | |
| 2956 BPGMetaData *md; | |
| 2957 BPGEncoderContext *enc_ctx; | |
| 2958 BPGEncoderParameters *p; | |
| 2959 | |
| 2960 p = bpg_encoder_param_alloc(); | |
| 2961 | |
| 2962 outfilename = DEFAULT_OUTFILENAME; | |
| 2963 color_space = BPG_CS_YCbCr; | |
| 2964 keep_metadata = 0; | |
| 2965 bit_depth = DEFAULT_BIT_DEPTH; | |
| 2966 limited_range = 0; | |
| 2967 premultiplied_alpha = 0; | |
| 2968 frame_delay_file = NULL; | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
2969 |
| 0 | 2970 for(;;) { |
| 2971 c = getopt_long_only(argc, argv, "q:o:hf:c:vm:b:e:a", long_opts, &option_index); | |
| 2972 if (c == -1) | |
| 2973 break; | |
| 2974 switch(c) { | |
| 2975 case 0: | |
| 2976 switch(option_index) { | |
| 2977 case 0: | |
| 2978 p->sei_decoded_picture_hash = 1; | |
| 2979 break; | |
| 2980 case 1: | |
| 2981 keep_metadata = 1; | |
| 2982 break; | |
| 2983 case 2: | |
| 2984 p->alpha_qp = atoi(optarg); | |
| 2985 if (p->alpha_qp < 0 || p->alpha_qp > 51) { | |
| 2986 fprintf(stderr, "alpha_qp must be between 0 and 51\n"); | |
| 2987 exit(1); | |
| 2988 } | |
| 2989 break; | |
| 2990 case 3: | |
| 2991 p->lossless = 1; | |
| 2992 color_space = BPG_CS_RGB; | |
| 2993 p->preferred_chroma_format = BPG_FORMAT_444; | |
| 2994 bit_depth = 8; | |
| 2995 limited_range = 0; | |
| 2996 break; | |
| 2997 case 4: | |
| 2998 limited_range = 1; | |
| 2999 break; | |
| 3000 case 5: | |
| 3001 premultiplied_alpha = 1; | |
| 3002 break; | |
| 3003 case 6: | |
| 3004 p->loop_count = strtoul(optarg, NULL, 0); | |
| 3005 break; | |
| 3006 case 7: | |
| 3007 p->frame_delay_num = 1; | |
| 3008 p->frame_delay_den = strtoul(optarg, NULL, 0); | |
| 3009 if (p->frame_delay_den == 0) { | |
| 3010 fprintf(stderr, "invalid frame rate\n"); | |
| 3011 exit(1); | |
| 3012 } | |
| 3013 break; | |
| 3014 case 8: | |
| 3015 frame_delay_file = optarg; | |
| 3016 break; | |
| 3017 default: | |
| 3018 goto show_help; | |
| 3019 } | |
| 3020 break; | |
| 3021 case 'h': | |
| 3022 show_help: | |
| 3023 help(1); | |
| 3024 break; | |
| 3025 case 'q': | |
| 3026 p->qp = atoi(optarg); | |
| 3027 if (p->qp < 0 || p->qp > 51) { | |
| 3028 fprintf(stderr, "qp must be between 0 and 51\n"); | |
| 3029 exit(1); | |
| 3030 } | |
| 3031 break; | |
| 3032 case 'o': | |
| 3033 outfilename = optarg; | |
| 3034 break; | |
| 3035 case 'f': | |
| 3036 if (!strcmp(optarg, "420")) { | |
| 3037 p->preferred_chroma_format = BPG_FORMAT_420; | |
| 3038 } else if (!strcmp(optarg, "422")) { | |
| 3039 p->preferred_chroma_format = BPG_FORMAT_422; | |
| 3040 } else if (!strcmp(optarg, "444")) { | |
| 3041 p->preferred_chroma_format = BPG_FORMAT_444; | |
| 3042 } else if (!strcmp(optarg, "422_video")) { | |
| 3043 p->preferred_chroma_format = BPG_FORMAT_422_VIDEO; | |
| 3044 } else if (!strcmp(optarg, "420_video")) { | |
| 3045 p->preferred_chroma_format = BPG_FORMAT_420_VIDEO; | |
| 3046 } else { | |
| 3047 fprintf(stderr, "Invalid chroma format\n"); | |
| 3048 exit(1); | |
| 3049 } | |
| 3050 break; | |
| 3051 case 'c': | |
| 3052 if (!strcmp(optarg, "ycbcr")) { | |
| 3053 color_space = BPG_CS_YCbCr; | |
| 3054 } else if (!strcmp(optarg, "rgb")) { | |
| 3055 color_space = BPG_CS_RGB; | |
| 3056 p->preferred_chroma_format = BPG_FORMAT_444; | |
| 3057 } else if (!strcmp(optarg, "ycgco")) { | |
| 3058 color_space = BPG_CS_YCgCo; | |
| 3059 } else if (!strcmp(optarg, "ycbcr_bt709")) { | |
| 3060 color_space = BPG_CS_YCbCr_BT709; | |
| 3061 } else if (!strcmp(optarg, "ycbcr_bt2020")) { | |
| 3062 color_space = BPG_CS_YCbCr_BT2020; | |
| 3063 } else { | |
| 3064 fprintf(stderr, "Invalid color space format\n"); | |
| 3065 exit(1); | |
| 3066 } | |
| 3067 break; | |
| 3068 case 'm': | |
| 3069 p->compress_level = atoi(optarg); | |
| 3070 if (p->compress_level < 1) | |
| 3071 p->compress_level = 1; | |
| 3072 else if (p->compress_level > 9) | |
| 3073 p->compress_level = 9; | |
| 3074 break; | |
| 3075 case 'b': | |
| 3076 bit_depth = atoi(optarg); | |
| 3077 if (bit_depth < 8 || bit_depth > BIT_DEPTH_MAX) { | |
| 3078 fprintf(stderr, "Invalid bit depth (range: 8 to %d)\n", | |
| 3079 BIT_DEPTH_MAX); | |
| 3080 exit(1); | |
| 3081 } | |
| 3082 break; | |
| 3083 case 'v': | |
| 3084 p->verbose++; | |
| 3085 break; | |
| 3086 case 'e': | |
| 3087 for(i = 0; i < HEVC_ENCODER_COUNT; i++) { | |
| 3088 if (!strcmp(optarg, hevc_encoder_name[i])) | |
| 3089 break; | |
| 3090 } | |
| 3091 if (i == HEVC_ENCODER_COUNT) { | |
| 3092 fprintf(stderr, "Unsupported encoder. Available ones are:"); | |
| 3093 for(i = 0; i < HEVC_ENCODER_COUNT; i++) { | |
| 3094 fprintf(stderr, " %s", hevc_encoder_name[i]); | |
| 3095 } | |
| 3096 fprintf(stderr, "\n"); | |
| 3097 exit(1); | |
| 3098 } | |
| 3099 p->encoder_type = i; | |
| 3100 break; | |
| 3101 case 'a': | |
| 3102 p->animated = 1; | |
| 3103 break; | |
| 3104 default: | |
| 3105 exit(1); | |
| 3106 } | |
| 3107 } | |
| 3108 | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
3109 if (optind >= argc) |
| 0 | 3110 help(0); |
| 3111 infilename = argv[optind]; | |
| 3112 | |
| 3113 f = fopen(outfilename, "wb"); | |
| 3114 if (!f) { | |
| 3115 perror(outfilename); | |
| 3116 exit(1); | |
| 3117 } | |
| 3118 | |
| 3119 enc_ctx = bpg_encoder_open(p); | |
| 3120 if (!enc_ctx) { | |
| 3121 fprintf(stderr, "Could not open BPG encoder\n"); | |
| 3122 exit(1); | |
| 3123 } | |
| 3124 | |
| 3125 if (p->animated) { | |
| 3126 int frame_num, first_frame, frame_ticks; | |
| 3127 char filename[1024]; | |
| 3128 FILE *f1; | |
| 3129 | |
| 3130 if (frame_delay_file) { | |
| 3131 f1 = fopen(frame_delay_file, "r"); | |
| 3132 if (!f1) { | |
| 3133 fprintf(stderr, "Could not open '%s'\n", frame_delay_file); | |
| 3134 exit(1); | |
| 3135 } | |
| 3136 } else { | |
| 3137 f1 = NULL; | |
| 3138 } | |
| 3139 | |
| 3140 first_frame = 1; | |
| 3141 for(frame_num = 0; ; frame_num++) { | |
| 3142 if (get_filename_num(filename, sizeof(filename), infilename, frame_num) < 0) { | |
| 3143 fprintf(stderr, "Invalid filename syntax: '%s'\n", infilename); | |
| 3144 exit(1); | |
| 3145 } | |
| 3146 img = load_image(&md, filename, color_space, bit_depth, limited_range, | |
| 3147 premultiplied_alpha); | |
| 3148 if (!img) { | |
| 3149 if (frame_num == 0) | |
| 3150 continue; /* accept to start at 0 or 1 */ | |
| 3151 if (first_frame) { | |
| 3152 fprintf(stderr, "Could not read '%s'\n", filename); | |
| 3153 exit(1); | |
| 3154 } else { | |
| 3155 break; | |
| 3156 } | |
| 3157 } | |
| 3158 frame_ticks = 1; | |
| 3159 if (f1) { | |
| 3160 float fdelay; | |
| 3161 if (fscanf(f1, "%f", &fdelay) == 1) { | |
| 3162 frame_ticks = lrint(fdelay * p->frame_delay_den / (p->frame_delay_num * 100)); | |
| 3163 if (frame_ticks < 1) | |
| 3164 frame_ticks = 1; | |
| 3165 } | |
| 3166 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
3167 |
| 0 | 3168 if (p->verbose) |
| 3169 printf("Encoding '%s' ticks=%d\n", filename, frame_ticks); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
3170 |
| 0 | 3171 if (keep_metadata && first_frame) { |
| 3172 bpg_encoder_set_extension_data(enc_ctx, md); | |
| 3173 } else { | |
| 3174 bpg_md_free(md); | |
| 3175 } | |
| 3176 bpg_encoder_set_frame_duration(enc_ctx, frame_ticks); | |
| 3177 bpg_encoder_encode(enc_ctx, img, my_write_func, f); | |
| 3178 image_free(img); | |
| 3179 | |
| 3180 first_frame = 0; | |
| 3181 } | |
| 3182 if (f1) | |
| 3183 fclose(f1); | |
| 3184 /* end of stream */ | |
| 3185 bpg_encoder_encode(enc_ctx, NULL, my_write_func, f); | |
| 3186 } else { | |
| 3187 img = load_image(&md, infilename, color_space, bit_depth, limited_range, | |
| 3188 premultiplied_alpha); | |
| 3189 if (!img) { | |
| 3190 fprintf(stderr, "Could not read '%s'\n", infilename); | |
| 3191 exit(1); | |
| 3192 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
3193 |
| 0 | 3194 if (!keep_metadata && md) { |
| 3195 bpg_md_free(md); | |
| 3196 md = NULL; | |
| 3197 } | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
3198 |
| 0 | 3199 bpg_encoder_set_extension_data(enc_ctx, md); |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
3200 |
| 0 | 3201 bpg_encoder_encode(enc_ctx, img, my_write_func, f); |
| 3202 image_free(img); | |
| 3203 } | |
| 3204 | |
| 3205 fclose(f); | |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
3206 |
| 0 | 3207 bpg_encoder_close(enc_ctx); |
|
5
524eae707ba4
Cosmetics: Remove trailing whitespace.
Matti Hamalainen <ccr@tnsp.org>
parents:
4
diff
changeset
|
3208 |
| 0 | 3209 bpg_encoder_param_free(p); |
| 3210 | |
| 3211 return 0; | |
| 3212 } |