FFmpeg coverage

Directory: ../../../ffmpeg/
File: src/libavcodec/jpeg2000dec.c
Date: 2026-10-19 01:23:10
Exec Total Coverage
Lines: 1193 1776 67.2%
Functions: 44 48 91.7%
Branches: 787 1350 58.3%
Line Branch Exec Source
1 /*
2 * JPEG 2000 image decoder
3 * Copyright (c) 2007 Kamil Nowosad
4 * Copyright (c) 2013 Nicolas Bertrand <nicoinattendu@gmail.com>
5 *
6 * This file is part of FFmpeg.
7 *
8 * FFmpeg is free software; you can redistribute it and/or
9 * modify it under the terms of the GNU Lesser General Public
10 * License as published by the Free Software Foundation; either
11 * version 2.1 of the License, or (at your option) any later version.
12 *
13 * FFmpeg is distributed in the hope that it will be useful,
14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
16 * Lesser General Public License for more details.
17 *
18 * You should have received a copy of the GNU Lesser General Public
19 * License along with FFmpeg; if not, write to the Free Software
20 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
21 */
22
23 /**
24 * @file
25 * JPEG 2000 image decoder
26 */
27
28 #include <inttypes.h>
29 #include <math.h>
30
31 #include "libavutil/attributes.h"
32 #include "libavutil/avassert.h"
33 #include "libavutil/common.h"
34 #include "libavutil/imgutils.h"
35 #include "libavutil/mem.h"
36 #include "libavutil/opt.h"
37 #include "libavutil/pixdesc.h"
38 #include "avcodec.h"
39 #include "bytestream.h"
40 #include "codec_internal.h"
41 #include "decode.h"
42 #include "internal.h"
43 #include "thread.h"
44 #include "jpeg2000.h"
45 #include "jpeg2000dsp.h"
46 #include "profiles.h"
47 #include "jpeg2000dec.h"
48 #include "jpeg2000htdec.h"
49
50 #define JP2_SIG_TYPE 0x6A502020
51 #define JP2_SIG_VALUE 0x0D0A870A
52 #define JP2_CODESTREAM 0x6A703263
53 #define JP2_HEADER 0x6A703268
54
55 #define HAD_COC 0x01
56 #define HAD_QCC 0x02
57
58 // Values of flag for placeholder passes
59 enum HT_PLHD_STATUS {
60 HT_PLHD_OFF,
61 HT_PLHD_ON
62 };
63
64 #define HT_MIXED 0x80 // bit 7 of SPcod/SPcoc
65
66
67 /* get_bits functions for JPEG2000 packet bitstream
68 * It is a get_bit function with a bit-stuffing routine. If the value of the
69 * byte is 0xFF, the next byte includes an extra zero bit stuffed into the MSB.
70 * cf. ISO-15444-1:2002 / B.10.1 Bit-stuffing routine */
71 8840179 static int get_bits(Jpeg2000DecoderContext *s, int n)
72 {
73 8840179 int res = 0;
74
75
2/2
✓ Branch 0 taken 17558593 times.
✓ Branch 1 taken 8840179 times.
 26398772 while (--n >= 0) {
76 17558593 res <<= 1;
77
2/2
✓ Branch 0 taken 2154406 times.
✓ Branch 1 taken 15404187 times.
 17558593 if (s->bit_index == 0) {
78
2/2
✓ Branch 1 taken 2093872 times.
✓ Branch 2 taken 60534 times.
 2154406 s->bit_index = 7 + (bytestream2_get_byte(&s->g) != 0xFFu);
79 }
80 17558593 s->bit_index--;
81 17558593 res |= (bytestream2_peek_byte(&s->g) >> s->bit_index) & 1;
82 }
83 8840179 return res;
84 }
85
86 87192 static void jpeg2000_flush(Jpeg2000DecoderContext *s)
87 {
88
2/2
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 87190 times.
 87192 if (bytestream2_get_byte(&s->g) == 0xff)
89 2 bytestream2_skip(&s->g, 1);
90 87192 s->bit_index = 8;
91 87192 }
92
93 /* decode the value stored in node */
94 1787992 static int tag_tree_decode(Jpeg2000DecoderContext *s, Jpeg2000TgtNode *node,
95 int threshold)
96 {
97 Jpeg2000TgtNode *stack[30];
98 1787992 int sp = -1, curval = 0;
99
100
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1787992 times.
 1787992 if (!node) {
101 ✗ av_log(s->avctx, AV_LOG_ERROR, "missing node\n");
102 ✗ return AVERROR_INVALIDDATA;
103 }
104
105
4/4
✓ Branch 0 taken 4250491 times.
✓ Branch 1 taken 385413 times.
✓ Branch 2 taken 2847912 times.
✓ Branch 3 taken 1402579 times.
 4635904 while (node && !node->vis) {
106 2847912 stack[++sp] = node;
107 2847912 node = node->parent;
108 }
109
110
2/2
✓ Branch 0 taken 1402579 times.
✓ Branch 1 taken 385413 times.
 1787992 if (node)
111 1402579 curval = node->val;
112 else
113 385413 curval = stack[sp]->val;
114
115
4/4
✓ Branch 0 taken 3943501 times.
✓ Branch 1 taken 161076 times.
✓ Branch 2 taken 2316585 times.
✓ Branch 3 taken 1626916 times.
 4104577 while (curval < threshold && sp >= 0) {
116
2/2
✓ Branch 0 taken 31252 times.
✓ Branch 1 taken 2285333 times.
 2316585 if (curval < stack[sp]->val)
117 31252 curval = stack[sp]->val;
118
2/2
✓ Branch 0 taken 3146479 times.
✓ Branch 1 taken 66521 times.
 3213000 while (curval < threshold) {
119 int ret;
120
2/2
✓ Branch 1 taken 2250064 times.
✓ Branch 2 taken 896415 times.
 3146479 if ((ret = get_bits(s, 1)) > 0) {
121 2250064 stack[sp]->vis++;
122 2250064 break;
123
1/2
✓ Branch 0 taken 896415 times.
✗ Branch 1 not taken.
 896415 } else if (!ret)
124 896415 curval++;
125 else
126 ✗ return ret;
127 }
128 2316585 stack[sp]->val = curval;
129 2316585 sp--;
130 }
131 1787992 return curval;
132 }
133
134 1309 static int pix_fmt_match(enum AVPixelFormat pix_fmt, int components,
135 int bpc, uint32_t log2_chroma_wh, int pal8)
136 {
137 1309 int match = 1;
138 1309 const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(pix_fmt);
139
140 av_assert2(desc);
141
142
2/2
✓ Branch 0 taken 319 times.
✓ Branch 1 taken 990 times.
 1309 if (desc->nb_components != components) {
143 319 return 0;
144 }
145
146
3/5
✓ Branch 0 taken 279 times.
✓ Branch 1 taken 680 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 31 times.
✗ Branch 4 not taken.
 990 switch (components) {
147 279 case 4:
148
2/2
✓ Branch 0 taken 242 times.
✓ Branch 1 taken 37 times.
 558 match = match && desc->comp[3].depth >= bpc &&
149
3/4
✓ Branch 0 taken 279 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 213 times.
✓ Branch 3 taken 29 times.
 771 (log2_chroma_wh >> 14 & 3) == 0 &&
150
1/2
✓ Branch 0 taken 213 times.
✗ Branch 1 not taken.
 213 (log2_chroma_wh >> 12 & 3) == 0;
151 av_fallthrough;
152 959 case 3:
153
2/2
✓ Branch 0 taken 886 times.
✓ Branch 1 taken 7 times.
 1852 match = match && desc->comp[2].depth >= bpc &&
154
4/4
✓ Branch 0 taken 893 times.
✓ Branch 1 taken 66 times.
✓ Branch 2 taken 833 times.
✓ Branch 3 taken 53 times.
 2685 (log2_chroma_wh >> 10 & 3) == desc->log2_chroma_w &&
155
2/2
✓ Branch 0 taken 831 times.
✓ Branch 1 taken 2 times.
 833 (log2_chroma_wh >> 8 & 3) == desc->log2_chroma_h;
156 av_fallthrough;
157 959 case 2:
158
1/2
✓ Branch 0 taken 831 times.
✗ Branch 1 not taken.
 1790 match = match && desc->comp[1].depth >= bpc &&
159
3/4
✓ Branch 0 taken 831 times.
✓ Branch 1 taken 128 times.
✓ Branch 2 taken 831 times.
✗ Branch 3 not taken.
 2621 (log2_chroma_wh >> 6 & 3) == desc->log2_chroma_w &&
160
1/2
✓ Branch 0 taken 831 times.
✗ Branch 1 not taken.
 831 (log2_chroma_wh >> 4 & 3) == desc->log2_chroma_h;
161 av_fallthrough;
162
163 990 case 1:
164
1/2
✓ Branch 0 taken 862 times.
✗ Branch 1 not taken.
 1852 match = match && desc->comp[0].depth >= bpc &&
165
2/2
✓ Branch 0 taken 853 times.
✓ Branch 1 taken 9 times.
 862 (log2_chroma_wh >> 2 & 3) == 0 &&
166
3/4
✓ Branch 0 taken 862 times.
✓ Branch 1 taken 128 times.
✓ Branch 2 taken 853 times.
✗ Branch 3 not taken.
 2705 (log2_chroma_wh & 3) == 0 &&
167
2/2
✓ Branch 0 taken 845 times.
✓ Branch 1 taken 8 times.
 853 (desc->flags & AV_PIX_FMT_FLAG_PAL) == pal8 * AV_PIX_FMT_FLAG_PAL;
168 }
169 990 return match;
170 }
171
172 // pix_fmts with lower bpp have to be listed before
173 // similar pix_fmts with higher bpp.
174 #define RGB_PIXEL_FORMATS AV_PIX_FMT_PAL8,AV_PIX_FMT_RGB24,AV_PIX_FMT_RGBA,AV_PIX_FMT_RGB48,AV_PIX_FMT_RGBA64
175 #define GRAY_PIXEL_FORMATS AV_PIX_FMT_GRAY8,AV_PIX_FMT_GRAY8A,AV_PIX_FMT_GRAY16,AV_PIX_FMT_YA16
176 #define YUV_PIXEL_FORMATS AV_PIX_FMT_YUV410P,AV_PIX_FMT_YUV411P,AV_PIX_FMT_YUVA420P, \
177 AV_PIX_FMT_YUV420P,AV_PIX_FMT_YUV422P,AV_PIX_FMT_YUVA422P, \
178 AV_PIX_FMT_YUV440P,AV_PIX_FMT_YUV444P,AV_PIX_FMT_YUVA444P, \
179 AV_PIX_FMT_YUV420P9,AV_PIX_FMT_YUV422P9,AV_PIX_FMT_YUV444P9, \
180 AV_PIX_FMT_YUVA420P9,AV_PIX_FMT_YUVA422P9,AV_PIX_FMT_YUVA444P9, \
181 AV_PIX_FMT_YUV420P10,AV_PIX_FMT_YUV422P10,AV_PIX_FMT_YUV444P10, \
182 AV_PIX_FMT_YUVA420P10,AV_PIX_FMT_YUVA422P10,AV_PIX_FMT_YUVA444P10, \
183 AV_PIX_FMT_YUV420P12,AV_PIX_FMT_YUV422P12,AV_PIX_FMT_YUV444P12, \
184 AV_PIX_FMT_YUV420P14,AV_PIX_FMT_YUV422P14,AV_PIX_FMT_YUV444P14, \
185 AV_PIX_FMT_YUV420P16,AV_PIX_FMT_YUV422P16,AV_PIX_FMT_YUV444P16, \
186 AV_PIX_FMT_YUVA420P16,AV_PIX_FMT_YUVA422P16,AV_PIX_FMT_YUVA444P16
187 #define XYZ_PIXEL_FORMATS AV_PIX_FMT_XYZ12
188
189 static const enum AVPixelFormat rgb_pix_fmts[] = {RGB_PIXEL_FORMATS};
190 static const enum AVPixelFormat gray_pix_fmts[] = {GRAY_PIXEL_FORMATS};
191 static const enum AVPixelFormat yuv_pix_fmts[] = {YUV_PIXEL_FORMATS};
192 static const enum AVPixelFormat xyz_pix_fmts[] = {XYZ_PIXEL_FORMATS,
193 YUV_PIXEL_FORMATS};
194 static const enum AVPixelFormat all_pix_fmts[] = {RGB_PIXEL_FORMATS,
195 GRAY_PIXEL_FORMATS,
196 YUV_PIXEL_FORMATS,
197 XYZ_PIXEL_FORMATS};
198
199 /* marker segments */
200 /* get sizes and offsets of image, tiles; number of components */
201 851 static int get_siz(Jpeg2000DecoderContext *s)
202 {
203 int i;
204 int ncomponents;
205 851 uint32_t log2_chroma_wh = 0;
206 851 const enum AVPixelFormat *possible_fmts = NULL;
207 851 int possible_fmts_nb = 0;
208 int ret;
209 int o_dimx, o_dimy; //original image dimensions.
210 int dimx, dimy;
211
212
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 851 times.
 851 if (bytestream2_get_bytes_left(&s->g) < 36) {
213 ✗ av_log(s->avctx, AV_LOG_ERROR, "Insufficient space for SIZ\n");
214 ✗ return AVERROR_INVALIDDATA;
215 }
216
217 851 s->avctx->profile = bytestream2_get_be16u(&s->g); // Rsiz
218 851 s->width = bytestream2_get_be32u(&s->g); // Width
219 851 s->height = bytestream2_get_be32u(&s->g); // Height
220 851 s->image_offset_x = bytestream2_get_be32u(&s->g); // X0Siz
221 851 s->image_offset_y = bytestream2_get_be32u(&s->g); // Y0Siz
222 851 s->tile_width = bytestream2_get_be32u(&s->g); // XTSiz
223 851 s->tile_height = bytestream2_get_be32u(&s->g); // YTSiz
224 851 s->tile_offset_x = bytestream2_get_be32u(&s->g); // XT0Siz
225 851 s->tile_offset_y = bytestream2_get_be32u(&s->g); // YT0Siz
226 851 ncomponents = bytestream2_get_be16u(&s->g); // CSiz
227
228
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 851 times.
 851 if (av_image_check_size2(s->width, s->height, s->avctx->max_pixels, AV_PIX_FMT_NONE, 0, s->avctx)) {
229 ✗ avpriv_request_sample(s->avctx, "Large Dimensions");
230 ✗ return AVERROR_PATCHWELCOME;
231 }
232
233
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 851 times.
 851 if (ncomponents <= 0) {
234 ✗ av_log(s->avctx, AV_LOG_ERROR, "Invalid number of components: %d\n",
235 s->ncomponents);
236 ✗ return AVERROR_INVALIDDATA;
237 }
238
239
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 851 times.
 851 if (ncomponents > 4) {
240 ✗ avpriv_request_sample(s->avctx, "Support for %d components",
241 ncomponents);
242 ✗ return AVERROR_PATCHWELCOME;
243 }
244
245
2/4
✓ Branch 0 taken 851 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 851 times.
✗ Branch 3 not taken.
 851 if (s->tile_offset_x < 0 || s->tile_offset_y < 0 ||
246
1/2
✓ Branch 0 taken 851 times.
✗ Branch 1 not taken.
 851 s->image_offset_x < s->tile_offset_x ||
247
1/2
✓ Branch 0 taken 851 times.
✗ Branch 1 not taken.
 851 s->image_offset_y < s->tile_offset_y ||
248
1/2
✓ Branch 0 taken 851 times.
✗ Branch 1 not taken.
 851 s->tile_width + (int64_t)s->tile_offset_x <= s->image_offset_x ||
249
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 851 times.
 851 s->tile_height + (int64_t)s->tile_offset_y <= s->image_offset_y
250 ) {
251 ✗ av_log(s->avctx, AV_LOG_ERROR, "Tile offsets are invalid\n");
252 ✗ return AVERROR_INVALIDDATA;
253 }
254
255
2/4
✓ Branch 0 taken 851 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 851 times.
 851 if (s->image_offset_x >= s->width || s->image_offset_y >= s->height) {
256 ✗ av_log(s->avctx, AV_LOG_ERROR, "image offsets outside image");
257 ✗ return AVERROR_INVALIDDATA;
258 }
259
260
4/6
✓ Branch 0 taken 1 times.
✓ Branch 1 taken 850 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
 851 if (s->reduction_factor && (s->image_offset_x || s->image_offset_y) ){
261 ✗ av_log(s->avctx, AV_LOG_ERROR, "reduction factor with image offsets is not fully implemented");
262 ✗ return AVERROR_PATCHWELCOME;
263 }
264
265 851 s->ncomponents = ncomponents;
266
267
2/4
✓ Branch 0 taken 851 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 851 times.
 851 if (s->tile_width <= 0 || s->tile_height <= 0) {
268 ✗ av_log(s->avctx, AV_LOG_ERROR, "Invalid tile dimension %dx%d.\n",
269 s->tile_width, s->tile_height);
270 ✗ return AVERROR_INVALIDDATA;
271 }
272
273
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 851 times.
 851 if (bytestream2_get_bytes_left(&s->g) < 3 * s->ncomponents) {
274 ✗ av_log(s->avctx, AV_LOG_ERROR, "Insufficient space for %d components in SIZ\n", s->ncomponents);
275 ✗ return AVERROR_INVALIDDATA;
276 }
277
278
2/2
✓ Branch 0 taken 851 times.
✓ Branch 1 taken 851 times.
 1702 for (i = 0; i < s->ncomponents; i++) {
279
1/2
✓ Branch 0 taken 851 times.
✗ Branch 1 not taken.
 851 if (s->cdef[i] < 0) {
280
2/2
✓ Branch 0 taken 2724 times.
✓ Branch 1 taken 851 times.
 3575 for (i = 0; i < s->ncomponents; i++) {
281 2724 s->cdef[i] = i + 1;
282 }
283
2/2
✓ Branch 0 taken 207 times.
✓ Branch 1 taken 644 times.
 851 if ((s->ncomponents & 1) == 0)
284 207 s->cdef[s->ncomponents-1] = 0;
285 }
286 }
287 // after here we no longer have to consider negative cdef
288
289 851 int cdef_used = 0;
290
2/2
✓ Branch 0 taken 2724 times.
✓ Branch 1 taken 851 times.
 3575 for (i = 0; i < s->ncomponents; i++)
291 2724 cdef_used |= 1<<s->cdef[i];
292
293 // Check that the channels we have are what we expect for the number of components
294
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 851 times.
 851 if (cdef_used != ((int[]){0,2,3,14,15})[s->ncomponents])
295 ✗ return AVERROR_INVALIDDATA;
296
297
2/2
✓ Branch 0 taken 2724 times.
✓ Branch 1 taken 851 times.
 3575 for (i = 0; i < s->ncomponents; i++) { // Ssiz_i XRsiz_i, YRsiz_i
298 2724 uint8_t x = bytestream2_get_byteu(&s->g);
299 2724 s->cbps[i] = (x & 0x7f) + 1;
300 2724 s->precision = FFMAX(s->cbps[i], s->precision);
301 2724 s->sgnd[i] = !!(x & 0x80);
302 2724 s->cdx[i] = bytestream2_get_byteu(&s->g);
303 2724 s->cdy[i] = bytestream2_get_byteu(&s->g);
304
3/6
✓ Branch 0 taken 2724 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2724 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 2724 times.
✗ Branch 5 not taken.
 2724 if ( !s->cdx[i] || s->cdx[i] == 3 || s->cdx[i] > 4
305
3/6
✓ Branch 0 taken 2724 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2724 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2724 times.
 2724 || !s->cdy[i] || s->cdy[i] == 3 || s->cdy[i] > 4) {
306 ✗ av_log(s->avctx, AV_LOG_ERROR, "Invalid sample separation %d/%d\n", s->cdx[i], s->cdy[i]);
307 ✗ return AVERROR_INVALIDDATA;
308 }
309
2/2
✓ Branch 0 taken 2517 times.
✓ Branch 1 taken 207 times.
 2724 int i_remapped = s->cdef[i] ? s->cdef[i]-1 : (s->ncomponents-1);
310
311 2724 log2_chroma_wh |= s->cdy[i] >> 1 << i_remapped * 4 | s->cdx[i] >> 1 << i_remapped * 4 + 2;
312 }
313
314 851 s->numXtiles = ff_jpeg2000_ceildiv(s->width - s->tile_offset_x, s->tile_width);
315 851 s->numYtiles = ff_jpeg2000_ceildiv(s->height - s->tile_offset_y, s->tile_height);
316
317 // There must be at least a SOT and SOD per tile, their minimum size is 14
318
1/2
✓ Branch 0 taken 851 times.
✗ Branch 1 not taken.
 851 if (s->numXtiles * (uint64_t)s->numYtiles > INT_MAX/sizeof(*s->tile) ||
319
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 851 times.
 851 s->numXtiles * s->numYtiles * 14LL > bytestream2_size(&s->g)
320 ) {
321 ✗ s->numXtiles = s->numYtiles = 0;
322 ✗ return AVERROR(EINVAL);
323 }
324
325 851 s->tile = av_calloc(s->numXtiles * s->numYtiles, sizeof(*s->tile));
326
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 851 times.
 851 if (!s->tile) {
327 ✗ s->numXtiles = s->numYtiles = 0;
328 ✗ return AVERROR(ENOMEM);
329 }
330
331
2/2
✓ Branch 0 taken 3216 times.
✓ Branch 1 taken 851 times.
 4067 for (i = 0; i < s->numXtiles * s->numYtiles; i++) {
332 3216 Jpeg2000Tile *tile = s->tile + i;
333
334 3216 tile->comp = av_mallocz(s->ncomponents * sizeof(*tile->comp));
335
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 3216 times.
 3216 if (!tile->comp)
336 ✗ return AVERROR(ENOMEM);
337 }
338
339 /* compute image size with reduction factor */
340 851 o_dimx = ff_jpeg2000_ceildivpow2(s->width - s->image_offset_x,
341 s->reduction_factor);
342 851 o_dimy = ff_jpeg2000_ceildivpow2(s->height - s->image_offset_y,
343 s->reduction_factor);
344 851 dimx = ff_jpeg2000_ceildiv(o_dimx, s->cdx[0]);
345 851 dimy = ff_jpeg2000_ceildiv(o_dimy, s->cdy[0]);
346
2/2
✓ Branch 0 taken 1873 times.
✓ Branch 1 taken 851 times.
 2724 for (i = 1; i < s->ncomponents; i++) {
347
2/2
✓ Branch 1 taken 1869 times.
✓ Branch 2 taken 4 times.
 1873 dimx = FFMAX(dimx, ff_jpeg2000_ceildiv(o_dimx, s->cdx[i]));
348
2/2
✓ Branch 1 taken 1869 times.
✓ Branch 2 taken 4 times.
 1873 dimy = FFMAX(dimy, ff_jpeg2000_ceildiv(o_dimy, s->cdy[i]));
349 }
350
351 851 ret = ff_set_dimensions(s->avctx, dimx << s->avctx->lowres, dimy << s->avctx->lowres);
352
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 851 times.
 851 if (ret < 0)
353 ✗ return ret;
354
355
2/2
✓ Branch 0 taken 848 times.
✓ Branch 1 taken 3 times.
 851 if (s->avctx->profile == AV_PROFILE_JPEG2000_DCINEMA_2K ||
356
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 848 times.
 848 s->avctx->profile == AV_PROFILE_JPEG2000_DCINEMA_4K) {
357 3 possible_fmts = xyz_pix_fmts;
358 3 possible_fmts_nb = FF_ARRAY_ELEMS(xyz_pix_fmts);
359 } else {
360
3/4
✓ Branch 0 taken 613 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 204 times.
✓ Branch 3 taken 31 times.
 848 switch (s->colour_space) {
361 613 case 16:
362 613 possible_fmts = rgb_pix_fmts;
363 613 possible_fmts_nb = FF_ARRAY_ELEMS(rgb_pix_fmts);
364 613 break;
365 ✗ case 17:
366 ✗ possible_fmts = gray_pix_fmts;
367 ✗ possible_fmts_nb = FF_ARRAY_ELEMS(gray_pix_fmts);
368 ✗ break;
369 204 case 18:
370 204 possible_fmts = yuv_pix_fmts;
371 204 possible_fmts_nb = FF_ARRAY_ELEMS(yuv_pix_fmts);
372 204 break;
373 31 default:
374 31 possible_fmts = all_pix_fmts;
375 31 possible_fmts_nb = FF_ARRAY_ELEMS(all_pix_fmts);
376 31 break;
377 }
378 }
379
2/2
✓ Branch 0 taken 817 times.
✓ Branch 1 taken 34 times.
 851 if ( s->avctx->pix_fmt != AV_PIX_FMT_NONE
380
2/2
✓ Branch 1 taken 3 times.
✓ Branch 2 taken 814 times.
 817 && !pix_fmt_match(s->avctx->pix_fmt, ncomponents, s->precision, log2_chroma_wh, s->pal8))
381 3 s->avctx->pix_fmt = AV_PIX_FMT_NONE;
382
2/2
✓ Branch 0 taken 37 times.
✓ Branch 1 taken 814 times.
 851 if (s->avctx->pix_fmt == AV_PIX_FMT_NONE)
383
2/2
✓ Branch 0 taken 492 times.
✓ Branch 1 taken 6 times.
 498 for (i = 0; i < possible_fmts_nb; ++i) {
384
2/2
✓ Branch 1 taken 31 times.
✓ Branch 2 taken 461 times.
 492 if (pix_fmt_match(possible_fmts[i], ncomponents, s->precision, log2_chroma_wh, s->pal8)) {
385 31 s->avctx->pix_fmt = possible_fmts[i];
386 31 break;
387 }
388 }
389
390
2/2
✓ Branch 0 taken 6 times.
✓ Branch 1 taken 845 times.
 851 if (i == possible_fmts_nb) {
391
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
 6 if (ncomponents == 4 &&
392
2/4
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
 2 s->cdy[0] == 1 && s->cdx[0] == 1 &&
393
2/4
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
 2 s->cdy[1] == 1 && s->cdx[1] == 1 &&
394
2/4
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
 2 s->cdy[2] == s->cdy[3] && s->cdx[2] == s->cdx[3]) {
395
4/8
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
 2 if (s->precision == 8 && s->cdy[2] == 2 && s->cdx[2] == 2 && !s->pal8) {
396 2 s->avctx->pix_fmt = AV_PIX_FMT_YUVA420P;
397 2 s->cdef[0] = 0;
398 2 s->cdef[1] = 1;
399 2 s->cdef[2] = 2;
400 2 s->cdef[3] = 3;
401 2 i = 0;
402 }
403
3/4
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
 4 } else if (ncomponents == 3 && s->precision == 8 &&
404
2/4
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
 2 s->cdx[0] == s->cdx[1] && s->cdx[0] == s->cdx[2] &&
405
2/4
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
 2 s->cdy[0] == s->cdy[1] && s->cdy[0] == s->cdy[2]) {
406 2 s->avctx->pix_fmt = AV_PIX_FMT_RGB24;
407 2 i = 0;
408
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
 2 } else if (ncomponents == 2 && s->precision == 8 &&
409 ✗ s->cdx[0] == s->cdx[1] && s->cdy[0] == s->cdy[1]) {
410 ✗ s->avctx->pix_fmt = AV_PIX_FMT_YA8;
411 ✗ i = 0;
412
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
 2 } else if (ncomponents == 2 && s->precision == 16 &&
413 ✗ s->cdx[0] == s->cdx[1] && s->cdy[0] == s->cdy[1]) {
414 ✗ s->avctx->pix_fmt = AV_PIX_FMT_YA16;
415 ✗ i = 0;
416
2/4
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
 2 } else if (ncomponents == 1 && s->precision == 8) {
417 2 s->avctx->pix_fmt = AV_PIX_FMT_GRAY8;
418 2 i = 0;
419 ✗ } else if (ncomponents == 1 && s->precision == 12) {
420 ✗ s->avctx->pix_fmt = AV_PIX_FMT_GRAY16LE;
421 ✗ i = 0;
422 }
423 }
424
425
426
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 851 times.
 851 if (i == possible_fmts_nb) {
427 ✗ av_log(s->avctx, AV_LOG_ERROR,
428 "Unknown pix_fmt, profile: %d, colour_space: %d, "
429 "components: %d, precision: %d\n"
430 "cdx[0]: %d, cdy[0]: %d\n"
431 "cdx[1]: %d, cdy[1]: %d\n"
432 "cdx[2]: %d, cdy[2]: %d\n"
433 "cdx[3]: %d, cdy[3]: %d\n",
434 ✗ s->avctx->profile, s->colour_space, ncomponents, s->precision,
435 s->cdx[0],
436 s->cdy[0],
437 ncomponents > 1 ? s->cdx[1] : 0,
438 ncomponents > 1 ? s->cdy[1] : 0,
439 ncomponents > 2 ? s->cdx[2] : 0,
440 ncomponents > 2 ? s->cdy[2] : 0,
441 ncomponents > 3 ? s->cdx[3] : 0,
442 ncomponents > 3 ? s->cdy[3] : 0);
443 ✗ return AVERROR_PATCHWELCOME;
444 }
445 851 s->avctx->bits_per_raw_sample = s->precision;
446 851 return 0;
447 }
448 /* get extended capabilities (CAP) marker segment */
449 2 static int get_cap(Jpeg2000DecoderContext *s, Jpeg2000CodingStyle *c)
450 {
451 uint32_t Pcap;
452 2 uint16_t Ccap_i[32] = { 0 };
453 uint16_t Ccap_15;
454 uint8_t P;
455
456
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
 2 if (bytestream2_get_bytes_left(&s->g) < 6) {
457 ✗ av_log(s->avctx, AV_LOG_ERROR, "Underflow while parsing the CAP marker\n");
458 ✗ return AVERROR_INVALIDDATA;
459 }
460
461 2 Pcap = bytestream2_get_be32u(&s->g);
462 2 s->isHT = (Pcap >> (31 - (15 - 1))) & 1;
463
2/2
✓ Branch 0 taken 64 times.
✓ Branch 1 taken 2 times.
 66 for (int i = 0; i < 32; i++) {
464
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 62 times.
 64 if ((Pcap >> (31 - i)) & 1)
465 2 Ccap_i[i] = bytestream2_get_be16u(&s->g);
466 }
467 2 Ccap_15 = Ccap_i[14];
468
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 2 if (s->isHT == 1) {
469 2 av_log(s->avctx, AV_LOG_INFO, "This codestream uses the HT block coder.\n");
470 // Bits 14-15
471
1/4
✗ Branch 0 not taken.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
 2 switch ((Ccap_15 >> 14) & 0x3) {
472 ✗ case 0x3:
473 ✗ s->Ccap15_b14_15 = HTJ2K_MIXED;
474 ✗ break;
475 ✗ case 0x1:
476 ✗ s->Ccap15_b14_15 = HTJ2K_HTDECLARED;
477 ✗ break;
478 2 case 0x0:
479 2 s->Ccap15_b14_15 = HTJ2K_HTONLY;
480 2 break;
481 ✗ default:
482 ✗ av_log(s->avctx, AV_LOG_ERROR, "Unknown CCap value.\n");
483 ✗ return AVERROR(EINVAL);
484 break;
485 }
486 // Bit 13
487
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
 2 if ((Ccap_15 >> 13) & 1) {
488 ✗ av_log(s->avctx, AV_LOG_ERROR, "MULTIHT set is not supported.\n");
489 ✗ return AVERROR_PATCHWELCOME;
490 }
491 // Bit 12
492 2 s->Ccap15_b12 = (Ccap_15 >> 12) & 1;
493 // Bit 11
494 2 s->Ccap15_b11 = (Ccap_15 >> 11) & 1;
495 // Bit 5
496 2 s->Ccap15_b05 = (Ccap_15 >> 5) & 1;
497 // Bit 0-4
498 2 P = Ccap_15 & 0x1F;
499
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
 2 if (!P)
500 ✗ s->HT_B = 8;
501
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 2 else if (P < 20)
502 2 s->HT_B = P + 8;
503 ✗ else if (P < 31)
504 ✗ s->HT_B = 4 * (P - 19) + 27;
505 else
506 ✗ s->HT_B = 74;
507
508
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
 2 if (s->HT_B > 31) {
509 ✗ av_log(s->avctx, AV_LOG_ERROR, "Codestream exceeds available precision (B > 31).\n");
510 ✗ return AVERROR_PATCHWELCOME;
511 }
512 }
513 2 return 0;
514 }
515
516 /* get common part for COD and COC segments */
517 869 static int get_cox(Jpeg2000DecoderContext *s, Jpeg2000CodingStyle *c)
518 {
519 uint8_t byte;
520
521
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 869 times.
 869 if (bytestream2_get_bytes_left(&s->g) < 5) {
522 ✗ av_log(s->avctx, AV_LOG_ERROR, "Insufficient space for COX\n");
523 ✗ return AVERROR_INVALIDDATA;
524 }
525
526 /* nreslevels = number of resolution levels
527 = number of decomposition level +1 */
528 869 c->nreslevels = bytestream2_get_byteu(&s->g) + 1;
529
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 869 times.
 869 if (c->nreslevels >= JPEG2000_MAX_RESLEVELS) {
530 ✗ av_log(s->avctx, AV_LOG_ERROR, "nreslevels %d is invalid\n", c->nreslevels);
531 ✗ return AVERROR_INVALIDDATA;
532 }
533
534
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 869 times.
 869 if (c->nreslevels <= s->reduction_factor) {
535 /* we are forced to update reduction_factor as its requested value is
536 not compatible with this bitstream, and as we might have used it
537 already in setup earlier we have to fail this frame until
538 reinitialization is implemented */
539 ✗ av_log(s->avctx, AV_LOG_ERROR, "reduction_factor too large for this bitstream, max is %d\n", c->nreslevels - 1);
540 ✗ s->reduction_factor = c->nreslevels - 1;
541 ✗ return AVERROR(EINVAL);
542 }
543
544 /* compute number of resolution levels to decode */
545 869 c->nreslevels2decode = c->nreslevels - s->reduction_factor;
546
547 869 c->log2_cblk_width = (bytestream2_get_byteu(&s->g) & 15) + 2; // cblk width
548 869 c->log2_cblk_height = (bytestream2_get_byteu(&s->g) & 15) + 2; // cblk height
549
550
2/4
✓ Branch 0 taken 869 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 869 times.
✗ Branch 3 not taken.
 869 if (c->log2_cblk_width > 10 || c->log2_cblk_height > 10 ||
551
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 869 times.
 869 c->log2_cblk_width + c->log2_cblk_height > 12) {
552 ✗ av_log(s->avctx, AV_LOG_ERROR, "cblk size invalid\n");
553 ✗ return AVERROR_INVALIDDATA;
554 }
555
556 869 c->cblk_style = bytestream2_get_byteu(&s->g);
557
2/2
✓ Branch 0 taken 12 times.
✓ Branch 1 taken 857 times.
 869 if (c->cblk_style != 0) { // cblk style
558
3/4
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 10 times.
 12 if (c->cblk_style & JPEG2000_CTSY_HTJ2K_M || c->cblk_style & JPEG2000_CTSY_HTJ2K_F) {
559 2 av_log(s->avctx,AV_LOG_TRACE,"High Throughput jpeg 2000 codestream.\n");
560 } else {
561 10 av_log(s->avctx, AV_LOG_WARNING, "extra cblk styles %X\n", c->cblk_style);
562
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
 10 if (c->cblk_style & JPEG2000_CBLK_BYPASS)
563 ✗ av_log(s->avctx, AV_LOG_WARNING, "Selective arithmetic coding bypass\n");
564 }
565 }
566 869 c->transform = bytestream2_get_byteu(&s->g); // DWT transformation type
567 /* set integer 9/7 DWT in case of BITEXACT flag */
568
4/4
✓ Branch 0 taken 867 times.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 223 times.
✓ Branch 3 taken 644 times.
 869 if ((s->avctx->flags & AV_CODEC_FLAG_BITEXACT) && (c->transform == FF_DWT97))
569 223 c->transform = FF_DWT97_INT;
570
571
2/2
✓ Branch 0 taken 16 times.
✓ Branch 1 taken 853 times.
 869 if (c->csty & JPEG2000_CSTY_PREC) {
572 int i;
573
2/2
✓ Branch 0 taken 84 times.
✓ Branch 1 taken 16 times.
 100 for (i = 0; i < c->nreslevels; i++) {
574 84 byte = bytestream2_get_byte(&s->g);
575 84 c->log2_prec_widths[i] = byte & 0x0F; // precinct PPx
576 84 c->log2_prec_heights[i] = (byte >> 4) & 0x0F; // precinct PPy
577
2/2
✓ Branch 0 taken 68 times.
✓ Branch 1 taken 16 times.
 84 if (i)
578
2/4
✓ Branch 0 taken 68 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 68 times.
 68 if (c->log2_prec_widths[i] == 0 || c->log2_prec_heights[i] == 0) {
579 ✗ av_log(s->avctx, AV_LOG_ERROR, "PPx %d PPy %d invalid\n",
580 ✗ c->log2_prec_widths[i], c->log2_prec_heights[i]);
581 ✗ c->log2_prec_widths[i] = c->log2_prec_heights[i] = 1;
582 ✗ return AVERROR_INVALIDDATA;
583 }
584 }
585 } else {
586 853 memset(c->log2_prec_widths , 15, sizeof(c->log2_prec_widths ));
587 853 memset(c->log2_prec_heights, 15, sizeof(c->log2_prec_heights));
588 }
589 869 return 0;
590 }
591
592 /* get coding parameters for a particular tile or whole image*/
593 851 static int get_cod(Jpeg2000DecoderContext *s, Jpeg2000CodingStyle *c,
594 const uint8_t *properties)
595 {
596 Jpeg2000CodingStyle tmp;
597 int compno, ret;
598
599
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 851 times.
 851 if (bytestream2_get_bytes_left(&s->g) < 5) {
600 ✗ av_log(s->avctx, AV_LOG_ERROR, "Insufficient space for COD\n");
601 ✗ return AVERROR_INVALIDDATA;
602 }
603
604 851 tmp.csty = bytestream2_get_byteu(&s->g);
605
606 // get progression order
607 851 tmp.prog_order = bytestream2_get_byteu(&s->g);
608
609 851 tmp.nlayers = bytestream2_get_be16u(&s->g);
610 851 tmp.mct = bytestream2_get_byteu(&s->g); // multiple component transformation
611
612
3/4
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 841 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 10 times.
 851 if (tmp.mct && s->ncomponents < 3) {
613 ✗ av_log(s->avctx, AV_LOG_ERROR,
614 "MCT %"PRIu8" with too few components (%d)\n",
615 ✗ tmp.mct, s->ncomponents);
616 ✗ return AVERROR_INVALIDDATA;
617 }
618
619
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 851 times.
 851 if ((ret = get_cox(s, &tmp)) < 0)
620 ✗ return ret;
621 851 tmp.init = 1;
622
2/2
✓ Branch 0 taken 2724 times.
✓ Branch 1 taken 851 times.
 3575 for (compno = 0; compno < s->ncomponents; compno++)
623
1/2
✓ Branch 0 taken 2724 times.
✗ Branch 1 not taken.
 2724 if (!(properties[compno] & HAD_COC))
624 2724 memcpy(c + compno, &tmp, sizeof(tmp));
625 851 return 0;
626 }
627
628 /* Get coding parameters for a component in the whole image or a
629 * particular tile. */
630 18 static int get_coc(Jpeg2000DecoderContext *s, Jpeg2000CodingStyle *c,
631 uint8_t *properties)
632 {
633 int compno, ret;
634 uint8_t has_eph, has_sop;
635
636
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 18 times.
 18 if (bytestream2_get_bytes_left(&s->g) < 2) {
637 ✗ av_log(s->avctx, AV_LOG_ERROR, "Insufficient space for COC\n");
638 ✗ return AVERROR_INVALIDDATA;
639 }
640
641 18 compno = bytestream2_get_byteu(&s->g);
642
643
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 18 times.
 18 if (compno >= s->ncomponents) {
644 ✗ av_log(s->avctx, AV_LOG_ERROR,
645 "Invalid compno %d. There are %d components in the image.\n",
646 compno, s->ncomponents);
647 ✗ return AVERROR_INVALIDDATA;
648 }
649
650 18 c += compno;
651 18 has_eph = c->csty & JPEG2000_CSTY_EPH;
652 18 has_sop = c->csty & JPEG2000_CSTY_SOP;
653 18 c->csty = bytestream2_get_byteu(&s->g);
654 18 c->csty |= has_eph; //do not override eph present bits from COD
655 18 c->csty |= has_sop; //do not override sop present bits from COD
656
657
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 18 times.
 18 if ((ret = get_cox(s, c)) < 0)
658 ✗ return ret;
659
660 18 properties[compno] |= HAD_COC;
661 18 c->init = 1;
662 18 return 0;
663 }
664
665 4 static int get_rgn(Jpeg2000DecoderContext *s, int n)
666 {
667 uint16_t compno;
668
1/2
✓ Branch 0 taken 4 times.
✗ Branch 1 not taken.
 4 compno = (s->ncomponents < 257)? bytestream2_get_byte(&s->g):
669 ✗ bytestream2_get_be16u(&s->g);
670
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 4 times.
 4 if (bytestream2_get_byte(&s->g)) {
671 ✗ av_log(s->avctx, AV_LOG_ERROR, "Invalid RGN header.\n");
672 ✗ return AVERROR_INVALIDDATA; // SRgn field value is 0
673 }
674 // SPrgn field
675 // Currently compno cannot be greater than 4.
676 // However, future implementation should support compno up to 65536
677
1/2
✓ Branch 0 taken 4 times.
✗ Branch 1 not taken.
 4 if (compno < s->ncomponents) {
678 int v;
679
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 4 times.
 4 if (s->curtileno == -1) {
680 ✗ v = bytestream2_get_byte(&s->g);
681 ✗ if (v > 30)
682 ✗ return AVERROR_PATCHWELCOME;
683 ✗ s->roi_shift[compno] = v;
684 } else {
685
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 4 times.
 4 if (s->tile[s->curtileno].tp_idx != 0)
686 ✗ return AVERROR_INVALIDDATA; // marker occurs only in first tile part of tile
687 4 v = bytestream2_get_byte(&s->g);
688
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 4 times.
 4 if (v > 30)
689 ✗ return AVERROR_PATCHWELCOME;
690 4 s->tile[s->curtileno].comp[compno].roi_shift = v;
691 }
692 4 return 0;
693 }
694 ✗ return AVERROR_INVALIDDATA;
695 }
696
697 /* Get common part for QCD and QCC segments. */
698 876 static int get_qcx(Jpeg2000DecoderContext *s, int n, Jpeg2000QuantStyle *q)
699 {
700 int i, x;
701
702
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 876 times.
 876 if (bytestream2_get_bytes_left(&s->g) < 1)
703 ✗ return AVERROR_INVALIDDATA;
704
705 876 x = bytestream2_get_byteu(&s->g); // Sqcd
706
707 876 q->nguardbits = x >> 5;
708 876 q->quantsty = x & 0x1f;
709
710
2/2
✓ Branch 0 taken 643 times.
✓ Branch 1 taken 233 times.
 876 if (q->quantsty == JPEG2000_QSTY_NONE) {
711 643 n -= 3;
712
2/4
✓ Branch 1 taken 643 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 643 times.
 643 if (bytestream2_get_bytes_left(&s->g) < n ||
713 n > JPEG2000_MAX_DECLEVELS*3)
714 ✗ return AVERROR_INVALIDDATA;
715
2/2
✓ Branch 0 taken 12004 times.
✓ Branch 1 taken 643 times.
 12647 for (i = 0; i < n; i++)
716 12004 q->expn[i] = bytestream2_get_byteu(&s->g) >> 3;
717
2/2
✓ Branch 0 taken 6 times.
✓ Branch 1 taken 227 times.
 233 } else if (q->quantsty == JPEG2000_QSTY_SI) {
718
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 6 times.
 6 if (bytestream2_get_bytes_left(&s->g) < 2)
719 ✗ return AVERROR_INVALIDDATA;
720 6 x = bytestream2_get_be16u(&s->g);
721 6 q->expn[0] = x >> 11;
722 6 q->mant[0] = x & 0x7ff;
723
2/2
✓ Branch 0 taken 588 times.
✓ Branch 1 taken 6 times.
 594 for (i = 1; i < JPEG2000_MAX_DECLEVELS * 3; i++) {
724 588 int curexpn = FFMAX(0, q->expn[0] - (i - 1) / 3);
725 588 q->expn[i] = curexpn;
726 588 q->mant[i] = q->mant[0];
727 }
728 } else {
729 227 n = (n - 3) >> 1;
730
2/4
✓ Branch 1 taken 227 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 227 times.
 227 if (bytestream2_get_bytes_left(&s->g) < 2 * n ||
731 n > JPEG2000_MAX_DECLEVELS*3)
732 ✗ return AVERROR_INVALIDDATA;
733
2/2
✓ Branch 0 taken 4268 times.
✓ Branch 1 taken 227 times.
 4495 for (i = 0; i < n; i++) {
734 4268 x = bytestream2_get_be16u(&s->g);
735 4268 q->expn[i] = x >> 11;
736 4268 q->mant[i] = x & 0x7ff;
737 }
738 }
739 876 return 0;
740 }
741
742 /* Get quantization parameters for a particular tile or a whole image. */
743 851 static int get_qcd(Jpeg2000DecoderContext *s, int n, Jpeg2000QuantStyle *q,
744 const uint8_t *properties)
745 {
746 Jpeg2000QuantStyle tmp;
747 int compno, ret;
748
749 851 memset(&tmp, 0, sizeof(tmp));
750
751
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 851 times.
 851 if ((ret = get_qcx(s, n, &tmp)) < 0)
752 ✗ return ret;
753
2/2
✓ Branch 0 taken 2724 times.
✓ Branch 1 taken 851 times.
 3575 for (compno = 0; compno < s->ncomponents; compno++)
754
1/2
✓ Branch 0 taken 2724 times.
✗ Branch 1 not taken.
 2724 if (!(properties[compno] & HAD_QCC))
755 2724 memcpy(q + compno, &tmp, sizeof(tmp));
756 851 return 0;
757 }
758
759 /* Get quantization parameters for a component in the whole image
760 * on in a particular tile. */
761 25 static int get_qcc(Jpeg2000DecoderContext *s, int n, Jpeg2000QuantStyle *q,
762 uint8_t *properties)
763 {
764 int compno;
765
766
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 25 times.
 25 if (bytestream2_get_bytes_left(&s->g) < 1)
767 ✗ return AVERROR_INVALIDDATA;
768
769 25 compno = bytestream2_get_byteu(&s->g);
770
771
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 25 times.
 25 if (compno >= s->ncomponents) {
772 ✗ av_log(s->avctx, AV_LOG_ERROR,
773 "Invalid compno %d. There are %d components in the image.\n",
774 compno, s->ncomponents);
775 ✗ return AVERROR_INVALIDDATA;
776 }
777
778 25 properties[compno] |= HAD_QCC;
779 25 return get_qcx(s, n - 1, q + compno);
780 }
781
782 8 static int get_poc(Jpeg2000DecoderContext *s, int size, Jpeg2000POC *p)
783 {
784 int i;
785
1/2
✓ Branch 0 taken 8 times.
✗ Branch 1 not taken.
 8 int elem_size = s->ncomponents <= 257 ? 7 : 9;
786 8 Jpeg2000POC tmp = {{{0}}};
787
788
2/4
✓ Branch 1 taken 8 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 8 times.
 8 if (bytestream2_get_bytes_left(&s->g) < 5 || size < 2 + elem_size) {
789 ✗ av_log(s->avctx, AV_LOG_ERROR, "Insufficient space for POC\n");
790 ✗ return AVERROR_INVALIDDATA;
791 }
792
793
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 8 times.
 8 if (elem_size > 7) {
794 ✗ avpriv_request_sample(s->avctx, "Fat POC not supported");
795 ✗ return AVERROR_PATCHWELCOME;
796 }
797
798 8 tmp.nb_poc = (size - 2) / elem_size;
799
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 8 times.
 8 if (tmp.nb_poc > MAX_POCS) {
800 ✗ avpriv_request_sample(s->avctx, "Too many POCs (%d)", tmp.nb_poc);
801 ✗ return AVERROR_PATCHWELCOME;
802 }
803
804
2/2
✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
 16 for (i = 0; i<tmp.nb_poc; i++) {
805 8 Jpeg2000POCEntry *e = &tmp.poc[i];
806 8 e->RSpoc = bytestream2_get_byteu(&s->g);
807 8 e->CSpoc = bytestream2_get_byteu(&s->g);
808 8 e->LYEpoc = bytestream2_get_be16u(&s->g);
809 8 e->REpoc = bytestream2_get_byteu(&s->g);
810 8 e->CEpoc = bytestream2_get_byteu(&s->g);
811 8 e->Ppoc = bytestream2_get_byteu(&s->g);
812
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 8 times.
 8 if (!e->CEpoc)
813 ✗ e->CEpoc = 256;
814
2/2
✓ Branch 0 taken 4 times.
✓ Branch 1 taken 4 times.
 8 if (e->CEpoc > s->ncomponents)
815 4 e->CEpoc = s->ncomponents;
816
2/4
✓ Branch 0 taken 8 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 8 times.
✗ Branch 3 not taken.
 8 if ( e->RSpoc >= e->REpoc || e->REpoc > 33
817
2/4
✓ Branch 0 taken 8 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 8 times.
✗ Branch 3 not taken.
 8 || e->CSpoc >= e->CEpoc || e->CEpoc > s->ncomponents
818
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 8 times.
 8 || !e->LYEpoc) {
819 ✗ av_log(s->avctx, AV_LOG_ERROR, "POC Entry %d is invalid (%d, %d, %d, %d, %d, %d)\n", i,
820 ✗ e->RSpoc, e->CSpoc, e->LYEpoc, e->REpoc, e->CEpoc, e->Ppoc
821 );
822 ✗ return AVERROR_INVALIDDATA;
823 }
824 }
825
826
3/4
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
 8 if (!p->nb_poc || p->is_default) {
827 6 *p = tmp;
828 } else {
829
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
 2 if (p->nb_poc + tmp.nb_poc > MAX_POCS) {
830 ✗ av_log(s->avctx, AV_LOG_ERROR, "Insufficient space for POC\n");
831 ✗ return AVERROR_INVALIDDATA;
832 }
833 2 memcpy(p->poc + p->nb_poc, tmp.poc, tmp.nb_poc * sizeof(tmp.poc[0]));
834 2 p->nb_poc += tmp.nb_poc;
835 }
836
837 8 p->is_default = 0;
838
839 8 return 0;
840 }
841
842
843 /* Get start of tile segment. */
844 3242 static int get_sot(Jpeg2000DecoderContext *s, int n)
845 {
846 Jpeg2000TilePart *tp;
847 uint16_t Isot;
848 uint32_t Psot;
849 unsigned TPsot;
850
851
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 3242 times.
 3242 if (bytestream2_get_bytes_left(&s->g) < 8)
852 ✗ return AVERROR_INVALIDDATA;
853
854 3242 s->curtileno = 0;
855 3242 Isot = bytestream2_get_be16u(&s->g); // Isot
856
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 3242 times.
 3242 if (Isot >= s->numXtiles * s->numYtiles)
857 ✗ return AVERROR_INVALIDDATA;
858
859 3242 s->curtileno = Isot;
860 3242 Psot = bytestream2_get_be32u(&s->g); // Psot
861 3242 TPsot = bytestream2_get_byteu(&s->g); // TPsot
862
863 /* Read TNSot but not used */
864 3242 bytestream2_get_byteu(&s->g); // TNsot
865
866
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 3242 times.
 3242 if (!Psot)
867 ✗ Psot = bytestream2_get_bytes_left(&s->g) - 2 + n + 2;
868
869
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 3242 times.
 3242 if (Psot > bytestream2_get_bytes_left(&s->g) - 2 + n + 2) {
870 ✗ av_log(s->avctx, AV_LOG_ERROR, "Psot %"PRIu32" too big\n", Psot);
871 ✗ return AVERROR_INVALIDDATA;
872 }
873
874
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 3242 times.
 3242 if (TPsot >= FF_ARRAY_ELEMS(s->tile[Isot].tile_part)) {
875 ✗ avpriv_request_sample(s->avctx, "Too many tile parts");
876 ✗ return AVERROR_PATCHWELCOME;
877 }
878
879 3242 s->tile[Isot].tp_idx = TPsot;
880 3242 tp = s->tile[Isot].tile_part + TPsot;
881 3242 tp->tile_index = Isot;
882 3242 tp->tp_end = s->g.buffer + Psot - n - 2;
883
884
2/2
✓ Branch 0 taken 3216 times.
✓ Branch 1 taken 26 times.
 3242 if (!TPsot) {
885 3216 Jpeg2000Tile *tile = s->tile + s->curtileno;
886
887 /* copy defaults */
888 3216 memcpy(tile->codsty, s->codsty, s->ncomponents * sizeof(Jpeg2000CodingStyle));
889 3216 memcpy(tile->qntsty, s->qntsty, s->ncomponents * sizeof(Jpeg2000QuantStyle));
890 3216 memcpy(&tile->poc , &s->poc , sizeof(tile->poc));
891 3216 tile->poc.is_default = 1;
892 }
893
894 3242 return 0;
895 }
896
897 4 static int read_crg(Jpeg2000DecoderContext *s, int n)
898 {
899
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 4 times.
 4 if (s->ncomponents*4 != n - 2) {
900 ✗ av_log(s->avctx, AV_LOG_ERROR, "Invalid CRG marker.\n");
901 ✗ return AVERROR_INVALIDDATA;
902 }
903 4 bytestream2_skip(&s->g, n - 2);
904 4 return 0;
905 }
906
907 2 static int read_cpf(Jpeg2000DecoderContext *s, int n)
908 {
909
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
 2 if (bytestream2_get_bytes_left(&s->g) < (n - 2))
910 ✗ return AVERROR_INVALIDDATA;
911 2 bytestream2_skip(&s->g, n - 2);
912 2 return 0;
913 }
914
915 /* Tile-part lengths: see ISO 15444-1:2002, section A.7.1
916 * Used to know the number of tile parts and lengths.
917 * There may be multiple TLMs in the header.
918 * TODO: The function is not used for tile-parts management, nor anywhere else.
919 * It can be useful to allocate memory for tile parts, before managing the SOT
920 * markers. Parsing the TLM header is needed to increment the input header
921 * buffer.
922 * This marker is mandatory for DCI. */
923 11 static int get_tlm(Jpeg2000DecoderContext *s, int n)
924 {
925 uint8_t Stlm, ST, SP, tile_tlm, i;
926 11 bytestream2_get_byte(&s->g); /* Ztlm: skipped */
927 11 Stlm = bytestream2_get_byte(&s->g);
928
929 // too complex ? ST = ((Stlm >> 4) & 0x01) + ((Stlm >> 4) & 0x02);
930 11 ST = (Stlm >> 4) & 0x03;
931
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 11 times.
 11 if (ST == 0x03) {
932 ✗ av_log(s->avctx, AV_LOG_ERROR, "TLM marker contains invalid ST value.\n");
933 ✗ return AVERROR_INVALIDDATA;
934 }
935
936 11 SP = (Stlm >> 6) & 0x01;
937 11 tile_tlm = (n - 4) / ((SP + 1) * 2 + ST);
938
2/2
✓ Branch 0 taken 38 times.
✓ Branch 1 taken 11 times.
 49 for (i = 0; i < tile_tlm; i++) {
939
2/4
✗ Branch 0 not taken.
✓ Branch 1 taken 17 times.
✓ Branch 2 taken 21 times.
✗ Branch 3 not taken.
 38 switch (ST) {
940 ✗ case 0:
941 ✗ break;
942 17 case 1:
943 17 bytestream2_get_byte(&s->g);
944 17 break;
945 21 case 2:
946 21 bytestream2_get_be16(&s->g);
947 21 break;
948 }
949
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 38 times.
 38 if (SP == 0) {
950 ✗ bytestream2_get_be16(&s->g);
951 } else {
952 38 bytestream2_get_be32(&s->g);
953 }
954 }
955 11 return 0;
956 }
957
958 514 static int get_plt(Jpeg2000DecoderContext *s, int n)
959 {
960 int i;
961 int v;
962
963 514 av_log(s->avctx, AV_LOG_DEBUG,
964 514 "PLT marker at pos 0x%X\n", bytestream2_tell(&s->g) - 4);
965
966
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 514 times.
 514 if (n < 4)
967 ✗ return AVERROR_INVALIDDATA;
968
969 514 /*Zplt =*/ bytestream2_get_byte(&s->g);
970
971
2/2
✓ Branch 0 taken 73452 times.
✓ Branch 1 taken 514 times.
 73966 for (i = 0; i < n - 3; i++) {
972 73452 v = bytestream2_get_byte(&s->g);
973 }
974
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 514 times.
 514 if (v & 0x80)
975 ✗ return AVERROR_INVALIDDATA;
976
977 514 return 0;
978 }
979
980 ✗ static int get_ppm(Jpeg2000DecoderContext *s, int n)
981 {
982 void *new;
983
984 ✗ if (n < 3) {
985 ✗ av_log(s->avctx, AV_LOG_ERROR, "Invalid length for PPM data.\n");
986 ✗ return AVERROR_INVALIDDATA;
987 }
988 ✗ bytestream2_get_byte(&s->g); //Zppm is skipped and not used
989 ✗ new = av_realloc(s->packed_headers,
990 ✗ s->packed_headers_size + n - 3);
991 ✗ if (new) {
992 ✗ s->packed_headers = new;
993 } else
994 ✗ return AVERROR(ENOMEM);
995 ✗ s->has_ppm = 1;
996 ✗ memset(&s->packed_headers_stream, 0, sizeof(s->packed_headers_stream));
997 ✗ bytestream2_get_bufferu(&s->g, s->packed_headers + s->packed_headers_size,
998 ✗ n - 3);
999 ✗ s->packed_headers_size += n - 3;
1000
1001 ✗ return 0;
1002 }
1003
1004 ✗ static int get_ppt(Jpeg2000DecoderContext *s, int n)
1005 {
1006 Jpeg2000Tile *tile;
1007 void *new;
1008
1009 ✗ if (n < 3) {
1010 ✗ av_log(s->avctx, AV_LOG_ERROR, "Invalid length for PPT data.\n");
1011 ✗ return AVERROR_INVALIDDATA;
1012 }
1013 ✗ if (s->curtileno < 0)
1014 ✗ return AVERROR_INVALIDDATA;
1015
1016 ✗ tile = &s->tile[s->curtileno];
1017 ✗ if (tile->tp_idx != 0) {
1018 ✗ av_log(s->avctx, AV_LOG_ERROR,
1019 "PPT marker can occur only on first tile part of a tile.\n");
1020 ✗ return AVERROR_INVALIDDATA;
1021 }
1022
1023 ✗ tile->has_ppt = 1; // this tile has a ppt marker
1024 ✗ bytestream2_get_byte(&s->g); // Zppt is skipped and not used
1025 ✗ new = av_realloc(tile->packed_headers,
1026 ✗ tile->packed_headers_size + n - 3);
1027 ✗ if (new) {
1028 ✗ tile->packed_headers = new;
1029 } else
1030 ✗ return AVERROR(ENOMEM);
1031 ✗ memset(&tile->packed_headers_stream, 0, sizeof(tile->packed_headers_stream));
1032 ✗ bytestream2_get_bufferu(&s->g, tile->packed_headers + tile->packed_headers_size, n - 3);
1033 ✗ tile->packed_headers_size += n - 3;
1034
1035 ✗ return 0;
1036 }
1037
1038 2881 static int init_tile(Jpeg2000DecoderContext *s, int tileno)
1039 {
1040 int compno;
1041 2881 int tilex = tileno % s->numXtiles;
1042 2881 int tiley = tileno / s->numXtiles;
1043 2881 Jpeg2000Tile *tile = s->tile + tileno;
1044
1045
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2881 times.
 2881 if (!tile->comp)
1046 ✗ return AVERROR(ENOMEM);
1047
1048 2881 tile->coord[0][0] = av_clip(tilex * (int64_t)s->tile_width + s->tile_offset_x, s->image_offset_x, s->width);
1049 2881 tile->coord[0][1] = av_clip((tilex + 1) * (int64_t)s->tile_width + s->tile_offset_x, s->image_offset_x, s->width);
1050 2881 tile->coord[1][0] = av_clip(tiley * (int64_t)s->tile_height + s->tile_offset_y, s->image_offset_y, s->height);
1051 2881 tile->coord[1][1] = av_clip((tiley + 1) * (int64_t)s->tile_height + s->tile_offset_y, s->image_offset_y, s->height);
1052
1053
2/2
✓ Branch 0 taken 9264 times.
✓ Branch 1 taken 2881 times.
 12145 for (compno = 0; compno < s->ncomponents; compno++) {
1054 9264 Jpeg2000Component *comp = tile->comp + compno;
1055 9264 Jpeg2000CodingStyle *codsty = tile->codsty + compno;
1056 9264 Jpeg2000QuantStyle *qntsty = tile->qntsty + compno;
1057 int ret; // global bandno
1058
1059 9264 comp->coord_o[0][0] = tile->coord[0][0];
1060 9264 comp->coord_o[0][1] = tile->coord[0][1];
1061 9264 comp->coord_o[1][0] = tile->coord[1][0];
1062 9264 comp->coord_o[1][1] = tile->coord[1][1];
1063
1064 9264 comp->coord_o[0][0] = ff_jpeg2000_ceildiv(comp->coord_o[0][0], s->cdx[compno]);
1065 9264 comp->coord_o[0][1] = ff_jpeg2000_ceildiv(comp->coord_o[0][1], s->cdx[compno]);
1066 9264 comp->coord_o[1][0] = ff_jpeg2000_ceildiv(comp->coord_o[1][0], s->cdy[compno]);
1067 9264 comp->coord_o[1][1] = ff_jpeg2000_ceildiv(comp->coord_o[1][1], s->cdy[compno]);
1068
1069 9264 comp->coord[0][0] = ff_jpeg2000_ceildivpow2(comp->coord_o[0][0], s->reduction_factor);
1070 9264 comp->coord[0][1] = ff_jpeg2000_ceildivpow2(comp->coord_o[0][1], s->reduction_factor);
1071 9264 comp->coord[1][0] = ff_jpeg2000_ceildivpow2(comp->coord_o[1][0], s->reduction_factor);
1072 9264 comp->coord[1][1] = ff_jpeg2000_ceildivpow2(comp->coord_o[1][1], s->reduction_factor);
1073
1074
2/2
✓ Branch 0 taken 9262 times.
✓ Branch 1 taken 2 times.
 9264 if (!comp->roi_shift)
1075 9262 comp->roi_shift = s->roi_shift[compno];
1076
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 9264 times.
 9264 if (!codsty->init)
1077 ✗ return AVERROR_INVALIDDATA;
1078
4/6
✓ Branch 0 taken 1 times.
✓ Branch 1 taken 9263 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
 9264 if (s->isHT && (!s->Ccap15_b05) && (!codsty->transform)) {
1079 ✗ av_log(s->avctx, AV_LOG_ERROR, "Transformation = 0 (lossy DWT) is found in HTREV HT set\n");
1080 ✗ return AVERROR_INVALIDDATA;
1081 }
1082
4/6
✓ Branch 0 taken 1 times.
✓ Branch 1 taken 9263 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
 9264 if (s->isHT && s->Ccap15_b14_15 != (codsty->cblk_style >> 6) && s->Ccap15_b14_15 != HTJ2K_HTONLY) {
1083 ✗ av_log(s->avctx, AV_LOG_ERROR, "SPcod/SPcoc value does not match bit 14-15 values of Ccap15\n");
1084 ✗ return AVERROR_INVALIDDATA;
1085 }
1086
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 9264 times.
 9264 if (ret = ff_jpeg2000_init_component(comp, codsty, qntsty,
1087 9264 s->cbps[compno], s->cdx[compno],
1088 s->cdy[compno], s->avctx))
1089 ✗ return ret;
1090 }
1091 2881 return 0;
1092 }
1093
1094 /* Read the number of coding passes. */
1095 890998 static int getnpasses(Jpeg2000DecoderContext *s)
1096 {
1097 int num;
1098
2/2
✓ Branch 1 taken 56247 times.
✓ Branch 2 taken 834751 times.
 890998 if (!get_bits(s, 1))
1099 56247 return 1;
1100
2/2
✓ Branch 1 taken 29425 times.
✓ Branch 2 taken 805326 times.
 834751 if (!get_bits(s, 1))
1101 29425 return 2;
1102
2/2
✓ Branch 1 taken 157732 times.
✓ Branch 2 taken 647594 times.
 805326 if ((num = get_bits(s, 2)) != 3)
1103
1/2
✓ Branch 0 taken 157732 times.
✗ Branch 1 not taken.
 157732 return num < 0 ? num : 3 + num;
1104
2/2
✓ Branch 1 taken 604660 times.
✓ Branch 2 taken 42934 times.
 647594 if ((num = get_bits(s, 5)) != 31)
1105
1/2
✓ Branch 0 taken 604660 times.
✗ Branch 1 not taken.
 604660 return num < 0 ? num : 6 + num;
1106 42934 num = get_bits(s, 7);
1107
1/2
✓ Branch 0 taken 42934 times.
✗ Branch 1 not taken.
 42934 return num < 0 ? num : 37 + num;
1108 }
1109
1110 890998 static int getlblockinc(Jpeg2000DecoderContext *s)
1111 {
1112 890998 int res = 0, ret;
1113
2/2
✓ Branch 1 taken 512784 times.
✓ Branch 2 taken 890998 times.
 1403782 while (ret = get_bits(s, 1)) {
1114
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 512784 times.
 512784 if (ret < 0)
1115 ✗ return ret;
1116 512784 res++;
1117 }
1118 890998 return res;
1119 }
1120
1121 ✗ static inline void select_header(Jpeg2000DecoderContext *s, const Jpeg2000Tile *tile,
1122 int *tp_index)
1123 {
1124 ✗ s->g = tile->tile_part[*tp_index].header_tpg;
1125 ✗ if (bytestream2_get_bytes_left(&s->g) == 0 && s->bit_index == 8) {
1126 ✗ av_log(s->avctx, AV_LOG_WARNING, "Packet header bytes in PPM marker segment is too short.\n");
1127 ✗ if (*tp_index < FF_ARRAY_ELEMS(tile->tile_part) - 1) {
1128 ✗ s->g = tile->tile_part[++(*tp_index)].tpg;
1129 }
1130 }
1131 ✗ }
1132
1133 87192 static inline void select_stream(Jpeg2000DecoderContext *s, const Jpeg2000Tile *tile,
1134 int *tp_index, const Jpeg2000CodingStyle *codsty)
1135 {
1136 int32_t is_endof_tp;
1137
1138 87192 s->g = tile->tile_part[*tp_index].tpg;
1139
3/4
✓ Branch 1 taken 9 times.
✓ Branch 2 taken 87183 times.
✓ Branch 3 taken 9 times.
✗ Branch 4 not taken.
 87192 is_endof_tp = bytestream2_get_bytes_left(&s->g) == 0 && s->bit_index == 8;
1140 // Following while loop is necessary because a tilepart may include only SOD marker.
1141 // Such a tilepart has neither packet header nor compressed data.
1142
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 87192 times.
 87202 while (is_endof_tp) {
1143
1/2
✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
 10 if (*tp_index < FF_ARRAY_ELEMS(tile->tile_part) - 1) {
1144 10 s->g = tile->tile_part[++(*tp_index)].tpg;
1145
3/4
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 9 times.
✓ Branch 3 taken 1 times.
✗ Branch 4 not taken.
 10 is_endof_tp = bytestream2_get_bytes_left(&s->g) == 0 && s->bit_index == 8;
1146 } else {
1147 ✗ is_endof_tp = 0;
1148 }
1149 }
1150
2/2
✓ Branch 0 taken 25452 times.
✓ Branch 1 taken 61740 times.
 87192 if (codsty->csty & JPEG2000_CSTY_SOP) {
1151
1/2
✓ Branch 1 taken 25452 times.
✗ Branch 2 not taken.
 25452 if (bytestream2_peek_be32(&s->g) == JPEG2000_SOP_FIXED_BYTES)
1152 25452 bytestream2_skip(&s->g, JPEG2000_SOP_BYTE_LENGTH);
1153 else
1154 ✗ av_log(s->avctx, AV_LOG_ERROR, "SOP marker not found. instead %X\n", bytestream2_peek_be32(&s->g));
1155 }
1156 87192 }
1157
1158 87264 static int jpeg2000_decode_packet(Jpeg2000DecoderContext *s, Jpeg2000Tile *tile, int *tp_index,
1159 const Jpeg2000CodingStyle *codsty,
1160 Jpeg2000ResLevel *rlevel, int precno,
1161 int layno, const uint8_t *expn, int numgbits)
1162 {
1163 int bandno, cblkno, ret, nb_code_blocks;
1164 int cwsno;
1165
1166
2/2
✓ Branch 0 taken 72 times.
✓ Branch 1 taken 87192 times.
 87264 if (layno < rlevel->band[0].prec[precno].decoded_layers)
1167 72 return 0;
1168 87192 rlevel->band[0].prec[precno].decoded_layers = layno + 1;
1169 // Select stream to read from
1170
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 87192 times.
 87192 if (s->has_ppm)
1171 ✗ select_header(s, tile, tp_index);
1172
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 87192 times.
 87192 else if (tile->has_ppt)
1173 ✗ s->g = tile->packed_headers_stream;
1174 else
1175 87192 select_stream(s, tile, tp_index, codsty);
1176
1177
2/2
✓ Branch 1 taken 3915 times.
✓ Branch 2 taken 83277 times.
 87192 if (!(ret = get_bits(s, 1))) {
1178 3915 jpeg2000_flush(s);
1179 3915 goto skip_data;
1180
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 83277 times.
 83277 } else if (ret < 0)
1181 ✗ return ret;
1182
1183
2/2
✓ Branch 0 taken 220079 times.
✓ Branch 1 taken 83277 times.
 303356 for (bandno = 0; bandno < rlevel->nbands; bandno++) {
1184 220079 Jpeg2000Band *band = rlevel->band + bandno;
1185 220079 Jpeg2000Prec *prec = band->prec + precno;
1186
1187
2/2
✓ Branch 0 taken 220077 times.
✓ Branch 1 taken 2 times.
 220079 if (band->coord[0][0] == band->coord[0][1] ||
1188
2/2
✓ Branch 0 taken 7190 times.
✓ Branch 1 taken 212887 times.
 220077 band->coord[1][0] == band->coord[1][1])
1189 7192 continue;
1190 212887 nb_code_blocks = prec->nb_codeblocks_height *
1191 212887 prec->nb_codeblocks_width;
1192
2/2
✓ Branch 0 taken 1001564 times.
✓ Branch 1 taken 212887 times.
 1214451 for (cblkno = 0; cblkno < nb_code_blocks; cblkno++) {
1193 1001564 Jpeg2000Cblk *cblk = prec->cblk + cblkno;
1194 int incl, newpasses, llen;
1195 void *tmp;
1196
1197
2/2
✓ Branch 0 taken 913967 times.
✓ Branch 1 taken 87597 times.
 1001564 if (!cblk->incl) {
1198 913967 incl = 0;
1199 913967 cblk->modes = codsty->cblk_style;
1200
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 913957 times.
 913967 if (cblk->modes >= JPEG2000_CTSY_HTJ2K_F)
1201 10 cblk->ht_plhd = HT_PLHD_ON;
1202
2/2
✓ Branch 0 taken 60567 times.
✓ Branch 1 taken 853400 times.
 913967 if (layno > 0)
1203 60567 incl = tag_tree_decode(s, prec->cblkincl + cblkno, 0 + 1) == 0;
1204 913967 incl = tag_tree_decode(s, prec->cblkincl + cblkno, layno + 1) == layno;
1205
1206
2/2
✓ Branch 0 taken 813458 times.
✓ Branch 1 taken 100509 times.
 913967 if (incl) {
1207 813458 int zbp = tag_tree_decode(s, prec->zerobits + cblkno, 100);
1208 813458 int v = expn[bandno] + numgbits - 1 - (zbp - tile->comp->roi_shift);
1209
2/4
✓ Branch 0 taken 813458 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 813458 times.
 813458 if (v < 0 || v > 30) {
1210 ✗ av_log(s->avctx, AV_LOG_ERROR,
1211 "nonzerobits %d invalid or unsupported\n", v);
1212 ✗ return AVERROR_INVALIDDATA;
1213 }
1214 813458 cblk->incl = 1;
1215 813458 cblk->nonzerobits = v;
1216 813458 cblk->zbp = zbp;
1217 813458 cblk->lblock = 3;
1218 }
1219 } else {
1220 87597 incl = get_bits(s, 1);
1221 }
1222
1223
2/2
✓ Branch 0 taken 890998 times.
✓ Branch 1 taken 110566 times.
 1001564 if (incl) {
1224 890998 uint8_t bypass_term_threshold = 0;
1225 890998 uint8_t bits_to_read = 0;
1226 890998 uint32_t segment_bytes = 0;
1227 890998 int32_t segment_passes = 0;
1228 890998 uint8_t next_segment_passes = 0;
1229 int32_t href_passes, pass_bound;
1230 890998 uint32_t tmp_length = 0;
1231 int32_t newpasses_copy, npasses_copy;
1232
1233
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 890998 times.
 890998 if ((newpasses = getnpasses(s)) <= 0)
1234 ✗ return newpasses;
1235
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 890998 times.
 890998 if (cblk->npasses + newpasses >= JPEG2000_MAX_PASSES) {
1236 ✗ avpriv_request_sample(s->avctx, "Too many passes");
1237 ✗ return AVERROR_PATCHWELCOME;
1238 }
1239
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 890998 times.
 890998 if ((llen = getlblockinc(s)) < 0)
1240 ✗ return llen;
1241
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 890998 times.
 890998 if (cblk->lblock + llen + av_log2(newpasses) > 16) {
1242 ✗ avpriv_request_sample(s->avctx,
1243 "Block with length beyond 16 bits");
1244 ✗ return AVERROR_PATCHWELCOME;
1245 }
1246 890998 cblk->nb_lengthinc = 0;
1247 890998 cblk->nb_terminationsinc = 0;
1248 890998 av_free(cblk->lengthinc);
1249 890998 cblk->lengthinc = av_calloc(newpasses, sizeof(*cblk->lengthinc));
1250
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 890998 times.
 890998 if (!cblk->lengthinc)
1251 ✗ return AVERROR(ENOMEM);
1252 890998 tmp = av_realloc_array(cblk->data_start, cblk->nb_terminations + newpasses + 1,
1253 sizeof(*cblk->data_start));
1254
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 890998 times.
 890998 if (!tmp)
1255 ✗ return AVERROR(ENOMEM);
1256 890998 cblk->data_start = tmp;
1257 890998 cblk->lblock += llen;
1258
1259 // Count number of necessary terminations for non HT code block
1260 890998 newpasses_copy = newpasses;
1261 890998 npasses_copy = cblk->npasses;
1262
2/2
✓ Branch 0 taken 890988 times.
✓ Branch 1 taken 10 times.
 890998 if (!(cblk->modes & JPEG2000_CTSY_HTJ2K_F)) {
1263 do {
1264 893516 int newpasses1 = 0;
1265
1266
2/2
✓ Branch 0 taken 12162772 times.
✓ Branch 1 taken 889354 times.
 13052126 while (newpasses1 < newpasses_copy) {
1267 12162772 newpasses1++;
1268
2/2
✓ Branch 1 taken 4162 times.
✓ Branch 2 taken 12158610 times.
 12162772 if (needs_termination(codsty->cblk_style, npasses_copy + newpasses1 - 1)) {
1269 4162 cblk->nb_terminationsinc++;
1270 4162 break;
1271 }
1272 }
1273 893516 npasses_copy += newpasses1;
1274 893516 newpasses_copy -= newpasses1;
1275
2/2
✓ Branch 0 taken 2528 times.
✓ Branch 1 taken 890988 times.
 893516 } while (newpasses_copy);
1276 }
1277
1278
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 890988 times.
 890998 if (cblk->ht_plhd) {
1279 10 href_passes = (cblk->npasses + newpasses - 1) % 3;
1280 10 segment_passes = newpasses - href_passes;
1281 10 pass_bound = 2;
1282 10 bits_to_read = cblk->lblock;
1283
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
 10 if (segment_passes < 1) {
1284 // No possible HT Cleanup pass here; may have placeholder passes
1285 // or an original J2K block bit-stream (in MIXED mode).
1286 ✗ segment_passes = newpasses;
1287 ✗ while (pass_bound <= segment_passes) {
1288 ✗ bits_to_read++;
1289 ✗ pass_bound += pass_bound;
1290 }
1291 ✗ segment_bytes = get_bits(s, bits_to_read);
1292 ✗ if (segment_bytes) {
1293 ✗ if (cblk->modes & HT_MIXED) {
1294 ✗ cblk->ht_plhd = HT_PLHD_OFF;
1295 ✗ cblk->modes &= (uint8_t) (~(JPEG2000_CTSY_HTJ2K_F));
1296 }
1297 else {
1298 ✗ av_log(s->avctx, AV_LOG_WARNING, "Length information for a HT-codeblock is invalid\n");
1299 }
1300 }
1301 } else {
1302
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
 10 while (pass_bound <= segment_passes) {
1303 ✗ bits_to_read++;
1304 ✗ pass_bound += pass_bound;
1305 }
1306 10 segment_bytes = get_bits(s, bits_to_read);
1307
1/2
✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
 10 if (segment_bytes) {
1308 // No more placeholder passes
1309
1/2
✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
 10 if (!(cblk->modes & HT_MIXED)) {
1310 // Must be the first HT Cleanup pass
1311
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
 10 if (segment_bytes < 2)
1312 ✗ av_log(s->avctx, AV_LOG_WARNING, "Length information for a HT-codeblock is invalid\n");
1313 10 next_segment_passes = 2;
1314 10 cblk->ht_plhd = HT_PLHD_OFF;
1315 // Write length information for HT CleanUp segment
1316 10 cblk->pass_lengths[0] = segment_bytes;
1317 ✗ } else if (cblk->lblock > 3 && segment_bytes > 1
1318 ✗ && (segment_bytes >> (bits_to_read - 1)) == 0) {
1319 // Must be the first HT Cleanup pass, since length MSB is 0
1320 ✗ next_segment_passes = 2;
1321 ✗ cblk->ht_plhd = HT_PLHD_OFF;
1322 // Write length information for HT CleanUp segment
1323 ✗ cblk->pass_lengths[0] = segment_bytes;
1324 } else {
1325 // Must have an original (non-HT) block coding pass
1326 ✗ cblk->modes &= (uint8_t) (~(JPEG2000_CTSY_HTJ2K_F));
1327 ✗ cblk->ht_plhd = HT_PLHD_OFF;
1328 ✗ segment_passes = newpasses;
1329 ✗ while (pass_bound <= segment_passes) {
1330 ✗ bits_to_read++;
1331 ✗ pass_bound += pass_bound;
1332 ✗ segment_bytes <<= 1;
1333 ✗ segment_bytes += get_bits(s, 1);
1334 }
1335 }
1336 } else {
1337 // Probably parsing placeholder passes, but we need to read an
1338 // extra length bit to verify this, since prior to the first
1339 // HT Cleanup pass, the number of length bits read for a
1340 // contributing code-block is dependent on the number of passes
1341 // being included, as if it were a non-HT code-block.
1342 ✗ segment_passes = newpasses;
1343 ✗ if (pass_bound <= segment_passes) {
1344 while (1) {
1345 ✗ bits_to_read++;
1346 ✗ pass_bound += pass_bound;
1347 ✗ segment_bytes <<= 1;
1348 ✗ segment_bytes += get_bits(s, 1);
1349 ✗ if (pass_bound > segment_passes)
1350 ✗ break;
1351 }
1352 ✗ if (segment_bytes) {
1353 ✗ if (cblk->modes & HT_MIXED) {
1354 ✗ cblk->modes &= (uint8_t) (~(JPEG2000_CTSY_HTJ2K_F));
1355 ✗ cblk->ht_plhd = HT_PLHD_OFF;
1356 } else {
1357 ✗ av_log(s->avctx, AV_LOG_WARNING, "Length information for a HT-codeblock is invalid\n");
1358 }
1359 }
1360 }
1361 }
1362 }
1363
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 890988 times.
 890988 } else if (cblk->modes & JPEG2000_CTSY_HTJ2K_F) {
1364 // Quality layer commences with a non-initial HT coding pass
1365 ✗ if(bits_to_read != 0)
1366 ✗ av_log(s->avctx, AV_LOG_WARNING, "Length information for a HT-codeblock is invalid\n");
1367 ✗ segment_passes = cblk->npasses % 3;
1368 ✗ if (segment_passes == 0) {
1369 // newpasses is a HT Cleanup pass; next segment has refinement passes
1370 ✗ segment_passes = 1;
1371 ✗ next_segment_passes = 2;
1372 ✗ if (segment_bytes == 1)
1373 ✗ av_log(s->avctx, AV_LOG_WARNING, "Length information for a HT-codeblock is invalid\n");
1374 } else {
1375 // newpasses = 1 means npasses is HT SigProp; 2 means newpasses is
1376 // HT MagRef pass
1377 ✗ segment_passes = newpasses > 1 ? 3 - segment_passes : 1;
1378 ✗ next_segment_passes = 1;
1379 ✗ bits_to_read = av_log2(segment_passes);
1380 }
1381 ✗ bits_to_read = (uint8_t) (bits_to_read + cblk->lblock);
1382 ✗ segment_bytes = get_bits(s, bits_to_read);
1383 // Write length information for HT Refinement segment
1384 ✗ cblk->pass_lengths[1] += segment_bytes;
1385
2/2
✓ Branch 0 taken 889354 times.
✓ Branch 1 taken 1634 times.
 890988 } else if (!(cblk->modes & (JPEG2000_CBLK_TERMALL | JPEG2000_CBLK_BYPASS))) {
1386 // Common case for non-HT code-blocks; we have only one segment
1387 889354 bits_to_read = (uint8_t) cblk->lblock + av_log2((uint8_t) newpasses);
1388 889354 segment_bytes = get_bits(s, bits_to_read);
1389 889354 segment_passes = newpasses;
1390
1/2
✓ Branch 0 taken 1634 times.
✗ Branch 1 not taken.
 1634 } else if (cblk->modes & JPEG2000_CBLK_TERMALL) {
1391 // RESTART MODE
1392 1634 bits_to_read = cblk->lblock;
1393 1634 segment_bytes = get_bits(s, bits_to_read);
1394 1634 segment_passes = 1;
1395 1634 next_segment_passes = 1;
1396 } else {
1397 // BYPASS MODE
1398 ✗ bypass_term_threshold = 10;
1399 ✗ if(bits_to_read != 0)
1400 ✗ av_log(s->avctx, AV_LOG_WARNING, "Length information for a codeblock is invalid\n");
1401 ✗ if (cblk->npasses < bypass_term_threshold) {
1402 // May have from 1 to 10 uninterrupted passes before 1st RAW SigProp
1403 ✗ segment_passes = bypass_term_threshold - cblk->npasses;
1404 ✗ if (segment_passes > newpasses)
1405 ✗ segment_passes = newpasses;
1406 ✗ while ((2 << bits_to_read) <= segment_passes)
1407 ✗ bits_to_read++;
1408 ✗ next_segment_passes = 2;
1409 ✗ } else if ((cblk->npasses - bypass_term_threshold) % 3 < 2) {
1410 // 0 means newpasses is a RAW SigProp; 1 means newpasses is a RAW MagRef pass
1411 ✗ segment_passes = newpasses > 1 ? 2 - (cblk->npasses - bypass_term_threshold) % 3 : 1;
1412 ✗ bits_to_read = av_log2(segment_passes);
1413 ✗ next_segment_passes = 1;
1414 } else {
1415 // newpasses is an isolated Cleanup pass that precedes a RAW SigProp pass
1416 ✗ segment_passes = 1;
1417 ✗ next_segment_passes = 2;
1418 }
1419 ✗ bits_to_read = (uint8_t) (bits_to_read + cblk->lblock);
1420 ✗ segment_bytes = get_bits(s, bits_to_read);
1421 }
1422 // Update cblk->npasses and write length information
1423 890998 cblk->npasses = (uint8_t) (cblk->npasses + segment_passes);
1424 890998 cblk->lengthinc[cblk->nb_lengthinc++] = segment_bytes;
1425
1426
3/4
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 890988 times.
✓ Branch 2 taken 10 times.
✗ Branch 3 not taken.
 890998 if ((cblk->modes & JPEG2000_CTSY_HTJ2K_F) && cblk->ht_plhd == HT_PLHD_OFF) {
1427 10 newpasses -= (uint8_t) segment_passes;
1428
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
 10 while (newpasses > 0) {
1429 ✗ segment_passes = newpasses > 1 ? next_segment_passes : 1;
1430 ✗ next_segment_passes = (uint8_t) (3 - next_segment_passes);
1431 ✗ bits_to_read = (uint8_t) (cblk->lblock + av_log2(segment_passes));
1432 ✗ segment_bytes = get_bits(s, bits_to_read);
1433 ✗ newpasses -= (uint8_t) (segment_passes);
1434 // This is a FAST Refinement pass
1435 // Write length information for HT Refinement segment
1436 ✗ cblk->pass_lengths[1] += segment_bytes;
1437 // Update cblk->npasses and write length information
1438 ✗ cblk->npasses = (uint8_t) (cblk->npasses + segment_passes);
1439 ✗ cblk->lengthinc[cblk->nb_lengthinc++] = segment_bytes;
1440 }
1441 } else {
1442 890988 newpasses -= (uint8_t) (segment_passes);
1443
2/2
✓ Branch 0 taken 2528 times.
✓ Branch 1 taken 890988 times.
 893516 while (newpasses > 0) {
1444
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2528 times.
 2528 if (bypass_term_threshold != 0) {
1445 ✗ segment_passes = newpasses > 1 ? next_segment_passes : 1;
1446 ✗ next_segment_passes = (uint8_t) (3 - next_segment_passes);
1447 ✗ bits_to_read = (uint8_t) (cblk->lblock + av_log2(segment_passes));
1448 } else {
1449
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2528 times.
 2528 if ((cblk->modes & JPEG2000_CBLK_TERMALL) == 0)
1450 ✗ av_log(s->avctx, AV_LOG_WARNING, "Corrupted packet header is found.\n");
1451 2528 segment_passes = 1;
1452 2528 bits_to_read = cblk->lblock;
1453 }
1454 2528 segment_bytes = get_bits(s, bits_to_read);
1455 2528 newpasses -= (uint8_t) (segment_passes);
1456
1457 // Update cblk->npasses and write length information
1458 2528 cblk->npasses = (uint8_t) (cblk->npasses + segment_passes);
1459 2528 cblk->lengthinc[cblk->nb_lengthinc++] = segment_bytes;
1460 }
1461 }
1462
1463
2/2
✓ Branch 0 taken 893526 times.
✓ Branch 1 taken 890998 times.
 1784524 for (int i = 0; i < cblk->nb_lengthinc; ++i)
1464 893526 tmp_length = (tmp_length < cblk->lengthinc[i]) ? cblk->lengthinc[i] : tmp_length;
1465
1466
2/2
✓ Branch 0 taken 822218 times.
✓ Branch 1 taken 68780 times.
 890998 if (tmp_length > cblk->data_allocated) {
1467 822218 size_t new_size = FFMAX(2 * cblk->data_allocated, tmp_length);
1468 822218 void *new = av_realloc(cblk->data, new_size);
1469
1/2
✓ Branch 0 taken 822218 times.
✗ Branch 1 not taken.
 822218 if (new) {
1470 822218 cblk->data = new;
1471 822218 cblk->data_allocated = new_size;
1472 }
1473 }
1474
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 890998 times.
 890998 if (tmp_length > cblk->data_allocated) {
1475 ✗ avpriv_request_sample(s->avctx,
1476 "Block with lengthinc greater than %zu",
1477 cblk->data_allocated);
1478 ✗ return AVERROR_PATCHWELCOME;
1479 }
1480 } else {
1481 // This codeblock has no contribution to the current packet
1482 110566 continue;
1483 }
1484 }
1485 }
1486 83277 jpeg2000_flush(s);
1487
1488
2/2
✓ Branch 0 taken 25321 times.
✓ Branch 1 taken 57956 times.
 83277 if (codsty->csty & JPEG2000_CSTY_EPH) {
1489
1/2
✓ Branch 1 taken 25321 times.
✗ Branch 2 not taken.
 25321 if (bytestream2_peek_be16(&s->g) == JPEG2000_EPH)
1490 25321 bytestream2_skip(&s->g, 2);
1491 else
1492 ✗ av_log(s->avctx, AV_LOG_ERROR, "EPH marker not found. instead %X\n", bytestream2_peek_be32(&s->g));
1493 }
1494
1495 // Save state of stream
1496
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 83277 times.
 83277 if (s->has_ppm) {
1497 ✗ tile->tile_part[*tp_index].header_tpg = s->g;
1498 ✗ select_stream(s, tile, tp_index, codsty);
1499
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 83277 times.
 83277 } else if (tile->has_ppt) {
1500 ✗ tile->packed_headers_stream = s->g;
1501 ✗ select_stream(s, tile, tp_index, codsty);
1502 }
1503
2/2
✓ Branch 0 taken 220079 times.
✓ Branch 1 taken 83277 times.
 303356 for (bandno = 0; bandno < rlevel->nbands; bandno++) {
1504 220079 Jpeg2000Band *band = rlevel->band + bandno;
1505 220079 Jpeg2000Prec *prec = band->prec + precno;
1506
1507 220079 nb_code_blocks = prec->nb_codeblocks_height * prec->nb_codeblocks_width;
1508
2/2
✓ Branch 0 taken 1008754 times.
✓ Branch 1 taken 220079 times.
 1228833 for (cblkno = 0; cblkno < nb_code_blocks; cblkno++) {
1509 1008754 Jpeg2000Cblk *cblk = prec->cblk + cblkno;
1510
4/4
✓ Branch 0 taken 1007120 times.
✓ Branch 1 taken 1634 times.
✓ Branch 2 taken 117756 times.
✓ Branch 3 taken 889364 times.
 1008754 if (!cblk->nb_terminationsinc && !cblk->lengthinc)
1511 117756 continue;
1512
2/2
✓ Branch 0 taken 893526 times.
✓ Branch 1 taken 890998 times.
 1784524 for (cwsno = 0; cwsno < cblk->nb_lengthinc; cwsno ++) {
1513
2/2
✓ Branch 0 taken 849862 times.
✓ Branch 1 taken 43664 times.
 893526 if (cblk->data_allocated < cblk->length + cblk->lengthinc[cwsno] + 4) {
1514 849862 size_t new_size = FFMAX(2*cblk->data_allocated, cblk->length + cblk->lengthinc[cwsno] + 4);
1515 849862 void *new = av_realloc(cblk->data, new_size);
1516
1/2
✓ Branch 0 taken 849862 times.
✗ Branch 1 not taken.
 849862 if (new) {
1517 849862 cblk->data = new;
1518 849862 cblk->data_allocated = new_size;
1519 }
1520 }
1521
1/2
✓ Branch 1 taken 893526 times.
✗ Branch 2 not taken.
 893526 if ( bytestream2_get_bytes_left(&s->g) < cblk->lengthinc[cwsno]
1522
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 893526 times.
 893526 || cblk->data_allocated < cblk->length + cblk->lengthinc[cwsno] + 4
1523 ) {
1524 ✗ av_log(s->avctx, AV_LOG_ERROR,
1525 "Block length %"PRIu16" or lengthinc %d is too large, left %d\n",
1526 ✗ cblk->length, cblk->lengthinc[cwsno], bytestream2_get_bytes_left(&s->g));
1527 ✗ return AVERROR_INVALIDDATA;
1528 }
1529
1530 893526 bytestream2_get_bufferu(&s->g, cblk->data + cblk->length, cblk->lengthinc[cwsno]);
1531 893526 cblk->length += cblk->lengthinc[cwsno];
1532 893526 memset(cblk->data + cblk->length, 0, 4);
1533 893526 cblk->lengthinc[cwsno] = 0;
1534
2/2
✓ Branch 0 taken 4162 times.
✓ Branch 1 taken 889364 times.
 893526 if (cblk->nb_terminationsinc) {
1535 4162 cblk->nb_terminationsinc--;
1536 4162 cblk->nb_terminations++;
1537 4162 cblk->data[cblk->length++] = 0xFF;
1538 4162 cblk->data[cblk->length++] = 0xFF;
1539 4162 cblk->data_start[cblk->nb_terminations] = cblk->length;
1540 }
1541 }
1542 890998 av_freep(&cblk->lengthinc);
1543 890998 cblk->nb_lengthinc = 0;
1544 }
1545 }
1546 // Save state of stream
1547 83277 tile->tile_part[*tp_index].tpg = s->g;
1548 83277 return 0;
1549
1550 3915 skip_data:
1551
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 3915 times.
 3915 if (codsty->csty & JPEG2000_CSTY_EPH) {
1552 ✗ if (bytestream2_peek_be16(&s->g) == JPEG2000_EPH)
1553 ✗ bytestream2_skip(&s->g, 2);
1554 else
1555 ✗ av_log(s->avctx, AV_LOG_ERROR, "EPH marker not found. instead %X\n", bytestream2_peek_be32(&s->g));
1556 }
1557
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 3915 times.
 3915 if (s->has_ppm) {
1558 ✗ tile->tile_part[*tp_index].header_tpg = s->g;
1559 ✗ select_stream(s, tile, tp_index, codsty);
1560
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 3915 times.
 3915 } else if (tile->has_ppt) {
1561 ✗ tile->packed_headers_stream = s->g;
1562 ✗ select_stream(s, tile, tp_index, codsty);
1563 }
1564 3915 tile->tile_part[*tp_index].tpg = s->g;
1565 3915 return 0;
1566 }
1567
1568 2882 static int jpeg2000_decode_packets_po_iteration(Jpeg2000DecoderContext *s, Jpeg2000Tile *tile,
1569 int RSpoc, int CSpoc,
1570 int LYEpoc, int REpoc, int CEpoc,
1571 int Ppoc, int *tp_index)
1572 {
1573 2882 int ret = 0;
1574 int layno, reslevelno, compno, precno, ok_reslevel;
1575 int x, y;
1576 int step_x, step_y;
1577
1578
4/6
✓ Branch 0 taken 259 times.
✓ Branch 1 taken 2619 times.
✓ Branch 2 taken 3 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
 2882 switch (Ppoc) {
1579 259 case JPEG2000_PGOD_RLCP:
1580 259 av_log(s->avctx, AV_LOG_DEBUG, "Progression order RLCP\n");
1581 259 ok_reslevel = 1;
1582
3/4
✓ Branch 0 taken 1298 times.
✓ Branch 1 taken 259 times.
✓ Branch 2 taken 1298 times.
✗ Branch 3 not taken.
 1557 for (reslevelno = RSpoc; ok_reslevel && reslevelno < REpoc; reslevelno++) {
1583 1298 ok_reslevel = 0;
1584
2/2
✓ Branch 0 taken 10385 times.
✓ Branch 1 taken 1298 times.
 11683 for (layno = 0; layno < LYEpoc; layno++) {
1585
2/2
✓ Branch 0 taken 31105 times.
✓ Branch 1 taken 10385 times.
 41490 for (compno = CSpoc; compno < CEpoc; compno++) {
1586 31105 Jpeg2000CodingStyle *codsty = tile->codsty + compno;
1587 31105 Jpeg2000QuantStyle *qntsty = tile->qntsty + compno;
1588
2/2
✓ Branch 0 taken 24920 times.
✓ Branch 1 taken 6185 times.
 31105 if (reslevelno < codsty->nreslevels) {
1589 24920 Jpeg2000ResLevel *rlevel = tile->comp[compno].reslevel +
1590 reslevelno;
1591 24920 ok_reslevel = 1;
1592
2/2
✓ Branch 0 taken 26420 times.
✓ Branch 1 taken 24920 times.
 51340 for (precno = 0; precno < rlevel->num_precincts_x * rlevel->num_precincts_y; precno++)
1593
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 26420 times.
 26420 if ((ret = jpeg2000_decode_packet(s, tile, tp_index,
1594 codsty, rlevel,
1595 precno, layno,
1596 26420 qntsty->expn + (reslevelno ? 3 * (reslevelno - 1) + 1 : 0),
1597
2/2
✓ Branch 0 taken 20235 times.
✓ Branch 1 taken 6185 times.
 26420 qntsty->nguardbits)) < 0)
1598 ✗ return ret;
1599 }
1600 }
1601 }
1602 }
1603 259 break;
1604
1605 2619 case JPEG2000_PGOD_LRCP:
1606 2619 av_log(s->avctx, AV_LOG_DEBUG, "Progression order LRCP\n");
1607
2/2
✓ Branch 0 taken 2698 times.
✓ Branch 1 taken 2619 times.
 5317 for (layno = 0; layno < LYEpoc; layno++) {
1608 2698 ok_reslevel = 1;
1609
4/4
✓ Branch 0 taken 21155 times.
✓ Branch 1 taken 2690 times.
✓ Branch 2 taken 21147 times.
✓ Branch 3 taken 8 times.
 23845 for (reslevelno = RSpoc; ok_reslevel && reslevelno < REpoc; reslevelno++) {
1610 21147 ok_reslevel = 0;
1611
2/2
✓ Branch 0 taken 68169 times.
✓ Branch 1 taken 21147 times.
 89316 for (compno = CSpoc; compno < CEpoc; compno++) {
1612 68169 Jpeg2000CodingStyle *codsty = tile->codsty + compno;
1613 68169 Jpeg2000QuantStyle *qntsty = tile->qntsty + compno;
1614
2/2
✓ Branch 0 taken 59595 times.
✓ Branch 1 taken 8574 times.
 68169 if (reslevelno < codsty->nreslevels) {
1615 59595 Jpeg2000ResLevel *rlevel = tile->comp[compno].reslevel +
1616 reslevelno;
1617 59595 ok_reslevel = 1;
1618
2/2
✓ Branch 0 taken 59595 times.
✓ Branch 1 taken 59595 times.
 119190 for (precno = 0; precno < rlevel->num_precincts_x * rlevel->num_precincts_y; precno++)
1619
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 59595 times.
 59595 if ((ret = jpeg2000_decode_packet(s, tile, tp_index,
1620 codsty, rlevel,
1621 precno, layno,
1622 59595 qntsty->expn + (reslevelno ? 3 * (reslevelno - 1) + 1 : 0),
1623
2/2
✓ Branch 0 taken 50997 times.
✓ Branch 1 taken 8598 times.
 59595 qntsty->nguardbits)) < 0)
1624 ✗ return ret;
1625 }
1626 }
1627 }
1628 }
1629 2619 break;
1630
1631 3 case JPEG2000_PGOD_CPRL:
1632 3 av_log(s->avctx, AV_LOG_DEBUG, "Progression order CPRL\n");
1633
2/2
✓ Branch 0 taken 9 times.
✓ Branch 1 taken 3 times.
 12 for (compno = CSpoc; compno < CEpoc; compno++) {
1634 9 Jpeg2000Component *comp = tile->comp + compno;
1635 9 Jpeg2000CodingStyle *codsty = tile->codsty + compno;
1636 9 Jpeg2000QuantStyle *qntsty = tile->qntsty + compno;
1637 9 step_x = 32;
1638 9 step_y = 32;
1639
1640
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 9 times.
 9 if (RSpoc >= FFMIN(codsty->nreslevels, REpoc))
1641 ✗ continue;
1642
1643
2/2
✓ Branch 0 taken 60 times.
✓ Branch 1 taken 9 times.
 69 for (reslevelno = RSpoc; reslevelno < FFMIN(codsty->nreslevels, REpoc); reslevelno++) {
1644 60 uint8_t reducedresno = codsty->nreslevels - 1 -reslevelno; // ==> N_L - r
1645 60 Jpeg2000ResLevel *rlevel = comp->reslevel + reslevelno;
1646 60 step_x = FFMIN(step_x, rlevel->log2_prec_width + reducedresno);
1647 60 step_y = FFMIN(step_y, rlevel->log2_prec_height + reducedresno);
1648 }
1649
2/4
✓ Branch 0 taken 9 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 9 times.
 9 if (step_x >= 31 || step_y >= 31){
1650 ✗ avpriv_request_sample(s->avctx, "CPRL with large step");
1651 ✗ return AVERROR_PATCHWELCOME;
1652 }
1653 9 step_x = 1<<step_x;
1654 9 step_y = 1<<step_y;
1655
1656
2/2
✓ Branch 0 taken 33 times.
✓ Branch 1 taken 9 times.
 42 for (y = tile->coord[1][0]; y < tile->coord[1][1]; y = (y/step_y + 1)*step_y) {
1657
2/2
✓ Branch 0 taken 243 times.
✓ Branch 1 taken 33 times.
 276 for (x = tile->coord[0][0]; x < tile->coord[0][1]; x = (x/step_x + 1)*step_x) {
1658
2/2
✓ Branch 0 taken 1464 times.
✓ Branch 1 taken 243 times.
 1707 for (reslevelno = RSpoc; reslevelno < FFMIN(codsty->nreslevels, REpoc); reslevelno++) {
1659 unsigned prcx, prcy;
1660 1464 uint8_t reducedresno = codsty->nreslevels - 1 -reslevelno; // ==> N_L - r
1661 1464 Jpeg2000ResLevel *rlevel = comp->reslevel + reslevelno;
1662 1464 int xc = x / s->cdx[compno];
1663 1464 int yc = y / s->cdy[compno];
1664
1665
3/4
✓ Branch 0 taken 816 times.
✓ Branch 1 taken 648 times.
✓ Branch 2 taken 816 times.
✗ Branch 3 not taken.
 1464 if (yc % (1LL << (rlevel->log2_prec_height + reducedresno)) && y != tile->coord[1][0]) //FIXME this is a subset of the check
1666 816 continue;
1667
1668
3/4
✓ Branch 0 taken 270 times.
✓ Branch 1 taken 378 times.
✓ Branch 2 taken 270 times.
✗ Branch 3 not taken.
 648 if (xc % (1LL << (rlevel->log2_prec_width + reducedresno)) && x != tile->coord[0][0]) //FIXME this is a subset of the check
1669 270 continue;
1670
1671 // check if a precinct exists
1672 378 prcx = ff_jpeg2000_ceildivpow2(xc, reducedresno) >> rlevel->log2_prec_width;
1673 378 prcy = ff_jpeg2000_ceildivpow2(yc, reducedresno) >> rlevel->log2_prec_height;
1674 378 prcx -= ff_jpeg2000_ceildivpow2(comp->coord_o[0][0], reducedresno) >> rlevel->log2_prec_width;
1675 378 prcy -= ff_jpeg2000_ceildivpow2(comp->coord_o[1][0], reducedresno) >> rlevel->log2_prec_height;
1676
1677 378 precno = prcx + rlevel->num_precincts_x * prcy;
1678
1679
2/4
✓ Branch 0 taken 378 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 378 times.
 378 if (prcx >= rlevel->num_precincts_x || prcy >= rlevel->num_precincts_y) {
1680 ✗ av_log(s->avctx, AV_LOG_WARNING, "prc %d %d outside limits %d %d\n",
1681 prcx, prcy, rlevel->num_precincts_x, rlevel->num_precincts_y);
1682 ✗ continue;
1683 }
1684
1685
2/2
✓ Branch 0 taken 1074 times.
✓ Branch 1 taken 378 times.
 1452 for (layno = 0; layno < LYEpoc; layno++) {
1686
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 1074 times.
 1074 if ((ret = jpeg2000_decode_packet(s, tile, tp_index, codsty, rlevel,
1687 precno, layno,
1688 1074 qntsty->expn + (reslevelno ? 3 * (reslevelno - 1) + 1 : 0),
1689
2/2
✓ Branch 0 taken 978 times.
✓ Branch 1 taken 96 times.
 1074 qntsty->nguardbits)) < 0)
1690 ✗ return ret;
1691 }
1692 }
1693 }
1694 }
1695 }
1696 3 break;
1697
1698 ✗ case JPEG2000_PGOD_RPCL:
1699 ✗ av_log(s->avctx, AV_LOG_WARNING, "Progression order RPCL\n");
1700 ✗ ok_reslevel = 1;
1701 ✗ for (reslevelno = RSpoc; ok_reslevel && reslevelno < REpoc; reslevelno++) {
1702 ✗ ok_reslevel = 0;
1703 ✗ step_x = 30;
1704 ✗ step_y = 30;
1705 ✗ for (compno = CSpoc; compno < CEpoc; compno++) {
1706 ✗ Jpeg2000Component *comp = tile->comp + compno;
1707 ✗ Jpeg2000CodingStyle *codsty = tile->codsty + compno;
1708
1709 ✗ if (reslevelno < codsty->nreslevels) {
1710 ✗ uint8_t reducedresno = codsty->nreslevels - 1 -reslevelno; // ==> N_L - r
1711 ✗ Jpeg2000ResLevel *rlevel = comp->reslevel + reslevelno;
1712 ✗ step_x = FFMIN(step_x, rlevel->log2_prec_width + reducedresno);
1713 ✗ step_y = FFMIN(step_y, rlevel->log2_prec_height + reducedresno);
1714 }
1715 }
1716 ✗ step_x = 1<<step_x;
1717 ✗ step_y = 1<<step_y;
1718
1719 ✗ for (y = tile->coord[1][0]; y < tile->coord[1][1]; y = (y/step_y + 1)*step_y) {
1720 ✗ for (x = tile->coord[0][0]; x < tile->coord[0][1]; x = (x/step_x + 1)*step_x) {
1721 ✗ for (compno = CSpoc; compno < CEpoc; compno++) {
1722 ✗ Jpeg2000Component *comp = tile->comp + compno;
1723 ✗ Jpeg2000CodingStyle *codsty = tile->codsty + compno;
1724 ✗ Jpeg2000QuantStyle *qntsty = tile->qntsty + compno;
1725 ✗ uint8_t reducedresno = codsty->nreslevels - 1 -reslevelno; // ==> N_L - r
1726 ✗ Jpeg2000ResLevel *rlevel = comp->reslevel + reslevelno;
1727 unsigned prcx, prcy;
1728 int trx0, try0;
1729
1730 ✗ if (!s->cdx[compno] || !s->cdy[compno])
1731 ✗ return AVERROR_INVALIDDATA;
1732
1733 ✗ if (reslevelno >= codsty->nreslevels)
1734 ✗ continue;
1735
1736 ✗ trx0 = ff_jpeg2000_ceildiv(tile->coord[0][0], (int64_t)s->cdx[compno] << reducedresno);
1737 ✗ try0 = ff_jpeg2000_ceildiv(tile->coord[1][0], (int64_t)s->cdy[compno] << reducedresno);
1738
1739 ✗ if (!(y % ((uint64_t)s->cdy[compno] << (rlevel->log2_prec_height + reducedresno)) == 0 ||
1740 ✗ (y == tile->coord[1][0] && ((int64_t)try0 << reducedresno) % (1ULL << (reducedresno + rlevel->log2_prec_height)))))
1741 ✗ continue;
1742
1743 ✗ if (!(x % ((uint64_t)s->cdx[compno] << (rlevel->log2_prec_width + reducedresno)) == 0 ||
1744 ✗ (x == tile->coord[0][0] && ((int64_t)trx0 << reducedresno) % (1ULL << (reducedresno + rlevel->log2_prec_width)))))
1745 ✗ continue;
1746
1747 // check if a precinct exists
1748 ✗ prcx = ff_jpeg2000_ceildiv(x, (int64_t)s->cdx[compno] << reducedresno) >> rlevel->log2_prec_width;
1749 ✗ prcy = ff_jpeg2000_ceildiv(y, (int64_t)s->cdy[compno] << reducedresno) >> rlevel->log2_prec_height;
1750 ✗ prcx -= ff_jpeg2000_ceildivpow2(comp->coord_o[0][0], reducedresno) >> rlevel->log2_prec_width;
1751 ✗ prcy -= ff_jpeg2000_ceildivpow2(comp->coord_o[1][0], reducedresno) >> rlevel->log2_prec_height;
1752
1753 ✗ precno = prcx + rlevel->num_precincts_x * prcy;
1754
1755 ✗ ok_reslevel = 1;
1756 ✗ if (prcx >= rlevel->num_precincts_x || prcy >= rlevel->num_precincts_y) {
1757 ✗ av_log(s->avctx, AV_LOG_WARNING, "prc %d %d outside limits %d %d\n",
1758 prcx, prcy, rlevel->num_precincts_x, rlevel->num_precincts_y);
1759 ✗ continue;
1760 }
1761
1762 ✗ for (layno = 0; layno < LYEpoc; layno++) {
1763 ✗ if ((ret = jpeg2000_decode_packet(s, tile, tp_index,
1764 codsty, rlevel,
1765 precno, layno,
1766 ✗ qntsty->expn + (reslevelno ? 3 * (reslevelno - 1) + 1 : 0),
1767 ✗ qntsty->nguardbits)) < 0)
1768 ✗ return ret;
1769 }
1770 }
1771 }
1772 }
1773 }
1774 ✗ break;
1775
1776 1 case JPEG2000_PGOD_PCRL:
1777 1 av_log(s->avctx, AV_LOG_WARNING, "Progression order PCRL\n");
1778 1 step_x = 32;
1779 1 step_y = 32;
1780
2/2
✓ Branch 0 taken 4 times.
✓ Branch 1 taken 1 times.
 5 for (compno = CSpoc; compno < CEpoc; compno++) {
1781 4 Jpeg2000Component *comp = tile->comp + compno;
1782 4 Jpeg2000CodingStyle *codsty = tile->codsty + compno;
1783
1784
2/2
✓ Branch 0 taken 25 times.
✓ Branch 1 taken 4 times.
 29 for (reslevelno = RSpoc; reslevelno < FFMIN(codsty->nreslevels, REpoc); reslevelno++) {
1785 25 uint8_t reducedresno = codsty->nreslevels - 1 -reslevelno; // ==> N_L - r
1786 25 Jpeg2000ResLevel *rlevel = comp->reslevel + reslevelno;
1787 25 step_x = FFMIN(step_x, rlevel->log2_prec_width + reducedresno);
1788 25 step_y = FFMIN(step_y, rlevel->log2_prec_height + reducedresno);
1789 }
1790 }
1791
2/4
✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 1 times.
 1 if (step_x >= 31 || step_y >= 31){
1792 ✗ avpriv_request_sample(s->avctx, "PCRL with large step");
1793 ✗ return AVERROR_PATCHWELCOME;
1794 }
1795 1 step_x = 1<<step_x;
1796 1 step_y = 1<<step_y;
1797
1798
2/2
✓ Branch 0 taken 1 times.
✓ Branch 1 taken 1 times.
 2 for (y = tile->coord[1][0]; y < tile->coord[1][1]; y = (y/step_y + 1)*step_y) {
1799
2/2
✓ Branch 0 taken 1 times.
✓ Branch 1 taken 1 times.
 2 for (x = tile->coord[0][0]; x < tile->coord[0][1]; x = (x/step_x + 1)*step_x) {
1800
2/2
✓ Branch 0 taken 4 times.
✓ Branch 1 taken 1 times.
 5 for (compno = CSpoc; compno < CEpoc; compno++) {
1801 4 Jpeg2000Component *comp = tile->comp + compno;
1802 4 Jpeg2000CodingStyle *codsty = tile->codsty + compno;
1803 4 Jpeg2000QuantStyle *qntsty = tile->qntsty + compno;
1804
1805
2/4
✓ Branch 0 taken 4 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 4 times.
 4 if (!s->cdx[compno] || !s->cdy[compno])
1806 ✗ return AVERROR_INVALIDDATA;
1807
1808
2/2
✓ Branch 0 taken 25 times.
✓ Branch 1 taken 4 times.
 29 for (reslevelno = RSpoc; reslevelno < FFMIN(codsty->nreslevels, REpoc); reslevelno++) {
1809 unsigned prcx, prcy;
1810 25 uint8_t reducedresno = codsty->nreslevels - 1 -reslevelno; // ==> N_L - r
1811 25 Jpeg2000ResLevel *rlevel = comp->reslevel + reslevelno;
1812 int trx0, try0;
1813
1814 25 trx0 = ff_jpeg2000_ceildiv(tile->coord[0][0], (int64_t)s->cdx[compno] << reducedresno);
1815 25 try0 = ff_jpeg2000_ceildiv(tile->coord[1][0], (int64_t)s->cdy[compno] << reducedresno);
1816
1817
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 25 times.
 25 if (!(y % ((uint64_t)s->cdy[compno] << (rlevel->log2_prec_height + reducedresno)) == 0 ||
1818 ✗ (y == tile->coord[1][0] && ((int64_t)try0 << reducedresno) % (1ULL << (reducedresno + rlevel->log2_prec_height)))))
1819 ✗ continue;
1820
1821
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 25 times.
 25 if (!(x % ((uint64_t)s->cdx[compno] << (rlevel->log2_prec_width + reducedresno)) == 0 ||
1822 ✗ (x == tile->coord[0][0] && ((int64_t)trx0 << reducedresno) % (1ULL << (reducedresno + rlevel->log2_prec_width)))))
1823 ✗ continue;
1824
1825 // check if a precinct exists
1826 25 prcx = ff_jpeg2000_ceildiv(x, (int64_t)s->cdx[compno] << reducedresno) >> rlevel->log2_prec_width;
1827 25 prcy = ff_jpeg2000_ceildiv(y, (int64_t)s->cdy[compno] << reducedresno) >> rlevel->log2_prec_height;
1828 25 prcx -= ff_jpeg2000_ceildivpow2(comp->coord_o[0][0], reducedresno) >> rlevel->log2_prec_width;
1829 25 prcy -= ff_jpeg2000_ceildivpow2(comp->coord_o[1][0], reducedresno) >> rlevel->log2_prec_height;
1830
1831 25 precno = prcx + rlevel->num_precincts_x * prcy;
1832
1833
2/4
✓ Branch 0 taken 25 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 25 times.
 25 if (prcx >= rlevel->num_precincts_x || prcy >= rlevel->num_precincts_y) {
1834 ✗ av_log(s->avctx, AV_LOG_WARNING, "prc %d %d outside limits %d %d\n",
1835 prcx, prcy, rlevel->num_precincts_x, rlevel->num_precincts_y);
1836 ✗ continue;
1837 }
1838
1839
2/2
✓ Branch 0 taken 175 times.
✓ Branch 1 taken 25 times.
 200 for (layno = 0; layno < LYEpoc; layno++) {
1840
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 175 times.
 175 if ((ret = jpeg2000_decode_packet(s, tile, tp_index, codsty, rlevel,
1841 precno, layno,
1842 175 qntsty->expn + (reslevelno ? 3 * (reslevelno - 1) + 1 : 0),
1843
2/2
✓ Branch 0 taken 147 times.
✓ Branch 1 taken 28 times.
 175 qntsty->nguardbits)) < 0)
1844 ✗ return ret;
1845 }
1846 }
1847 }
1848 }
1849 }
1850 1 break;
1851
1852 ✗ default:
1853 ✗ break;
1854 }
1855
1856 2882 return ret;
1857 }
1858
1859 2881 static int jpeg2000_decode_packets(Jpeg2000DecoderContext *s, Jpeg2000Tile *tile)
1860 {
1861 2881 int ret = AVERROR_BUG;
1862 int i;
1863 2881 int tp_index = 0;
1864
1865 2881 s->bit_index = 8;
1866
2/2
✓ Branch 0 taken 9 times.
✓ Branch 1 taken 2872 times.
 2881 if (tile->poc.nb_poc) {
1867
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 9 times.
 19 for (i=0; i<tile->poc.nb_poc; i++) {
1868 10 Jpeg2000POCEntry *e = &tile->poc.poc[i];
1869 30 ret = jpeg2000_decode_packets_po_iteration(s, tile,
1870 10 e->RSpoc, e->CSpoc,
1871 10 FFMIN(e->LYEpoc, tile->codsty[0].nlayers),
1872 10 e->REpoc,
1873 10 FFMIN(e->CEpoc, s->ncomponents),
1874
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
 10 e->Ppoc, &tp_index
1875 );
1876
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
 10 if (ret < 0)
1877 ✗ return ret;
1878 }
1879 } else {
1880 2872 ret = jpeg2000_decode_packets_po_iteration(s, tile,
1881 0, 0,
1882 2872 tile->codsty[0].nlayers,
1883 33,
1884 s->ncomponents,
1885 2872 tile->codsty[0].prog_order,
1886 &tp_index
1887 );
1888 }
1889 /* EOC marker reached */
1890 2881 bytestream2_skip(&s->g, 2);
1891
1892 2881 return ret;
1893 }
1894
1895 /* TIER-1 routines */
1896 4025580 static void decode_sigpass(Jpeg2000T1Context *t1, int width, int height,
1897 int bpno, int bandno,
1898 int vert_causal_ctx_csty_symbol)
1899 {
1900 4025580 int mask = (3u << bpno)>>1, y0, x, y;
1901
1902
2/2
✓ Branch 0 taken 14247108 times.
✓ Branch 1 taken 4025580 times.
 18272688 for (y0 = 0; y0 < height; y0 += 4)
1903
2/2
✓ Branch 0 taken 235531408 times.
✓ Branch 1 taken 14247108 times.
 249778516 for (x = 0; x < width; x++)
1904
4/4
✓ Branch 0 taken 1115006649 times.
✓ Branch 1 taken 57650602 times.
✓ Branch 2 taken 937125843 times.
✓ Branch 3 taken 177880806 times.
 1172657251 for (y = y0; y < height && y < y0 + 4; y++) {
1905 937125843 int flags_mask = -1;
1906
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 937125843 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
 937125843 if (vert_causal_ctx_csty_symbol && y == y0 + 3)
1907 ✗ flags_mask &= ~(JPEG2000_T1_SIG_S | JPEG2000_T1_SIG_SW | JPEG2000_T1_SIG_SE | JPEG2000_T1_SGN_S);
1908
2/2
✓ Branch 0 taken 572729325 times.
✓ Branch 1 taken 364396518 times.
 937125843 if ((t1->flags[(y+1) * t1->stride + x+1] & JPEG2000_T1_SIG_NB & flags_mask)
1909
2/2
✓ Branch 0 taken 255814734 times.
✓ Branch 1 taken 316914591 times.
 572729325 && !(t1->flags[(y+1) * t1->stride + x+1] & (JPEG2000_T1_SIG | JPEG2000_T1_VIS))) {
1910
2/2
✓ Branch 2 taken 80968267 times.
✓ Branch 3 taken 174846467 times.
 255814734 if (ff_mqc_decode(&t1->mqc, t1->mqc.cx_states + ff_jpeg2000_getsigctxno(t1->flags[(y+1) * t1->stride + x+1] & flags_mask, bandno))) {
1911 80968267 int xorbit, ctxno = ff_jpeg2000_getsgnctxno(t1->flags[(y+1) * t1->stride + x+1] & flags_mask, &xorbit);
1912
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 80968267 times.
 80968267 if (t1->mqc.raw) {
1913 ✗ t1->data[(y) * t1->stride + x] |= (uint32_t)(ff_mqc_decode(&t1->mqc, t1->mqc.cx_states + ctxno)) << 31;
1914 ✗ t1->data[(y) * t1->stride + x] |= mask;
1915 } else {
1916 80968267 t1->data[(y) * t1->stride + x] |= (uint32_t)(ff_mqc_decode(&t1->mqc, t1->mqc.cx_states + ctxno) ^ xorbit) << 31;
1917 80968267 t1->data[(y) * t1->stride + x] |= mask;
1918 }
1919 80968267 ff_jpeg2000_set_significance(t1, x, y,
1920 80968267 t1->data[(y) * t1->stride + x] & INT32_MIN);
1921 }
1922 255814734 t1->flags[(y + 1) * t1->stride + x + 1] |= JPEG2000_T1_VIS;
1923 }
1924 }
1925 4025580 }
1926
1927 3688064 static void decode_refpass(Jpeg2000T1Context *t1, int width, int height,
1928 int bpno, int vert_causal_ctx_csty_symbol)
1929 {
1930 int phalf;
1931 int y0, x, y;
1932
1933 3688064 phalf = 1 << bpno;
1934
1935
2/2
✓ Branch 0 taken 13057007 times.
✓ Branch 1 taken 3688064 times.
 16745071 for (y0 = 0; y0 < height; y0 += 4)
1936
2/2
✓ Branch 0 taken 217413142 times.
✓ Branch 1 taken 13057007 times.
 230470149 for (x = 0; x < width; x++)
1937
4/4
✓ Branch 0 taken 1029719586 times.
✓ Branch 1 taken 52786672 times.
✓ Branch 2 taken 865093116 times.
✓ Branch 3 taken 164626470 times.
 1082506258 for (y = y0; y < height && y < y0 + 4; y++)
1938
2/2
✓ Branch 0 taken 279193434 times.
✓ Branch 1 taken 585899682 times.
 865093116 if ((t1->flags[(y + 1) * t1->stride + x + 1] & (JPEG2000_T1_SIG | JPEG2000_T1_VIS)) == JPEG2000_T1_SIG) {
1939 ✗ int flags_mask = (vert_causal_ctx_csty_symbol && y == y0 + 3) ?
1940
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 279193434 times.
 279193434 ~(JPEG2000_T1_SIG_S | JPEG2000_T1_SIG_SW | JPEG2000_T1_SIG_SE | JPEG2000_T1_SGN_S) : -1;
1941 279193434 int ctxno = ff_jpeg2000_getrefctxno(t1->flags[(y + 1) * t1->stride + x + 1] & flags_mask);
1942 279193434 t1->data[(y) * t1->stride + x] |= phalf >> 1;
1943
2/2
✓ Branch 1 taken 104556358 times.
✓ Branch 2 taken 174637076 times.
 279193434 if (ff_mqc_decode(&t1->mqc, t1->mqc.cx_states + ctxno))
1944 104556358 t1->data[(y) * t1->stride + x] |= phalf;
1945 else {
1946 174637076 t1->data[(y) * t1->stride + x] &= ~(phalf);
1947
1948 }
1949 279193434 t1->flags[(y + 1) * t1->stride + x + 1] |= JPEG2000_T1_REF;
1950 }
1951 3688064 }
1952
1953 4399286 static void decode_clnpass(const Jpeg2000DecoderContext *s, Jpeg2000T1Context *t1,
1954 int width, int height, int bpno, int bandno,
1955 int seg_symbols, int vert_causal_ctx_csty_symbol)
1956 {
1957 4399286 int mask = (3u << bpno)>>1, y0, x, y, runlen, dec;
1958
1959
2/2
✓ Branch 0 taken 15675400 times.
✓ Branch 1 taken 4399286 times.
 20074686 for (y0 = 0; y0 < height; y0 += 4) {
1960
2/2
✓ Branch 0 taken 260469372 times.
✓ Branch 1 taken 15675400 times.
 276144772 for (x = 0; x < width; x++) {
1961 260469372 int flags_mask = -1;
1962
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 260469372 times.
 260469372 if (vert_causal_ctx_csty_symbol)
1963 ✗ flags_mask &= ~(JPEG2000_T1_SIG_S | JPEG2000_T1_SIG_SW | JPEG2000_T1_SIG_SE | JPEG2000_T1_SGN_S);
1964
2/2
✓ Branch 0 taken 258378523 times.
✓ Branch 1 taken 2090849 times.
 260469372 if (y0 + 3 < height &&
1965
2/2
✓ Branch 0 taken 116152811 times.
✓ Branch 1 taken 142225712 times.
 258378523 !((t1->flags[(y0 + 1) * t1->stride + x + 1] & (JPEG2000_T1_SIG_NB | JPEG2000_T1_VIS | JPEG2000_T1_SIG)) ||
1966
2/2
✓ Branch 0 taken 107696334 times.
✓ Branch 1 taken 8456477 times.
 116152811 (t1->flags[(y0 + 2) * t1->stride + x + 1] & (JPEG2000_T1_SIG_NB | JPEG2000_T1_VIS | JPEG2000_T1_SIG)) ||
1967
2/2
✓ Branch 0 taken 100366742 times.
✓ Branch 1 taken 7329592 times.
 107696334 (t1->flags[(y0 + 3) * t1->stride + x + 1] & (JPEG2000_T1_SIG_NB | JPEG2000_T1_VIS | JPEG2000_T1_SIG)) ||
1968
2/2
✓ Branch 0 taken 96872806 times.
✓ Branch 1 taken 3493936 times.
 100366742 (t1->flags[(y0 + 4) * t1->stride + x + 1] & (JPEG2000_T1_SIG_NB | JPEG2000_T1_VIS | JPEG2000_T1_SIG) & flags_mask))) {
1969
2/2
✓ Branch 1 taken 90566361 times.
✓ Branch 2 taken 6306445 times.
 96872806 if (!ff_mqc_decode(&t1->mqc, t1->mqc.cx_states + MQC_CX_RL))
1970 90566361 continue;
1971 6306445 runlen = ff_mqc_decode(&t1->mqc,
1972 6306445 t1->mqc.cx_states + MQC_CX_UNI);
1973 12612890 runlen = (runlen << 1) | ff_mqc_decode(&t1->mqc,
1974 6306445 t1->mqc.cx_states +
1975 MQC_CX_UNI);
1976 6306445 dec = 1;
1977 } else {
1978 163596566 runlen = 0;
1979 163596566 dec = 0;
1980 }
1981
1982
4/4
✓ Branch 0 taken 666468227 times.
✓ Branch 1 taken 167812162 times.
✓ Branch 2 taken 664377378 times.
✓ Branch 3 taken 2090849 times.
 834280389 for (y = y0 + runlen; y < y0 + 4 && y < height; y++) {
1983 664377378 int flags_mask = -1;
1984
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 664377378 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
 664377378 if (vert_causal_ctx_csty_symbol && y == y0 + 3)
1985 ✗ flags_mask &= ~(JPEG2000_T1_SIG_S | JPEG2000_T1_SIG_SW | JPEG2000_T1_SIG_SE | JPEG2000_T1_SGN_S);
1986
2/2
✓ Branch 0 taken 658070933 times.
✓ Branch 1 taken 6306445 times.
 664377378 if (!dec) {
1987
2/2
✓ Branch 0 taken 153966202 times.
✓ Branch 1 taken 504104731 times.
 658070933 if (!(t1->flags[(y+1) * t1->stride + x+1] & (JPEG2000_T1_SIG | JPEG2000_T1_VIS))) {
1988 153966202 dec = ff_mqc_decode(&t1->mqc, t1->mqc.cx_states + ff_jpeg2000_getsigctxno(t1->flags[(y+1) * t1->stride + x+1] & flags_mask,
1989 bandno));
1990 }
1991 }
1992
2/2
✓ Branch 0 taken 24012194 times.
✓ Branch 1 taken 640365184 times.
 664377378 if (dec) {
1993 int xorbit;
1994 24012194 int ctxno = ff_jpeg2000_getsgnctxno(t1->flags[(y + 1) * t1->stride + x + 1] & flags_mask,
1995 &xorbit);
1996 24012194 t1->data[(y) * t1->stride + x] |= (uint32_t)(ff_mqc_decode(&t1->mqc, t1->mqc.cx_states + ctxno) ^ xorbit) << 31;
1997 24012194 t1->data[(y) * t1->stride + x] |= mask;
1998 24012194 ff_jpeg2000_set_significance(t1, x, y, t1->data[(y) * t1->stride + x] & INT32_MIN);
1999 }
2000 664377378 dec = 0;
2001 664377378 t1->flags[(y + 1) * t1->stride + x + 1] &= ~JPEG2000_T1_VIS;
2002 }
2003 }
2004 }
2005
2/2
✓ Branch 0 taken 104 times.
✓ Branch 1 taken 4399182 times.
 4399286 if (seg_symbols) {
2006 int val;
2007 104 val = ff_mqc_decode(&t1->mqc, t1->mqc.cx_states + MQC_CX_UNI);
2008 104 val = (val << 1) + ff_mqc_decode(&t1->mqc, t1->mqc.cx_states + MQC_CX_UNI);
2009 104 val = (val << 1) + ff_mqc_decode(&t1->mqc, t1->mqc.cx_states + MQC_CX_UNI);
2010 104 val = (val << 1) + ff_mqc_decode(&t1->mqc, t1->mqc.cx_states + MQC_CX_UNI);
2011
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 104 times.
 104 if (val != 0xa)
2012 ✗ av_log(s->avctx, AV_LOG_ERROR,
2013 "Segmentation symbol value incorrect\n");
2014 }
2015 4399286 }
2016
2017 919837 static int decode_cblk(const Jpeg2000DecoderContext *s, Jpeg2000CodingStyle *codsty,
2018 Jpeg2000T1Context *t1, Jpeg2000Cblk *cblk,
2019 int width, int height, int bandpos, uint8_t roi_shift, const int M_b)
2020 {
2021 919837 int passno = cblk->npasses, pass_t = 2, bpno = cblk->nonzerobits - 1 + 31 - M_b - 1 - roi_shift;
2022 919837 int pass_cnt = 0;
2023 919837 int vert_causal_ctx_csty_symbol = codsty->cblk_style & JPEG2000_CBLK_VSC;
2024 919837 int term_cnt = 0;
2025 int coder_type;
2026
2027
2/4
✓ Branch 0 taken 919837 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 919837 times.
 919837 av_assert0(width <= 1024U && height <= 1024U);
2028
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 919837 times.
 919837 av_assert0(width*height <= 4096);
2029
2030 919837 memset(t1->data, 0, t1->stride * height * sizeof(*t1->data));
2031
2032 /* If code-block contains no compressed data: nothing to do. */
2033
2/2
✓ Branch 0 taken 108213 times.
✓ Branch 1 taken 811624 times.
 919837 if (!cblk->length)
2034 108213 return 0;
2035
2036 811624 memset(t1->flags, 0, t1->stride * (height + 2) * sizeof(*t1->flags));
2037
2038 811624 cblk->data[cblk->length] = 0xff;
2039 811624 cblk->data[cblk->length+1] = 0xff;
2040 811624 ff_mqc_initdec(&t1->mqc, cblk->data, 0, 1);
2041
2042
2/2
✓ Branch 0 taken 12112930 times.
✓ Branch 1 taken 811624 times.
 12924554 while (passno--) {
2043
2/4
✓ Branch 0 taken 12112930 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 12112930 times.
 12112930 if (bpno < -1 || bpno > 29) {
2044 ✗ av_log(s->avctx, AV_LOG_ERROR, "bpno (%d) became invalid\n", bpno);
2045 ✗ return AVERROR_INVALIDDATA;
2046 }
2047
3/4
✓ Branch 0 taken 4025580 times.
✓ Branch 1 taken 3688064 times.
✓ Branch 2 taken 4399286 times.
✗ Branch 3 not taken.
 12112930 switch(pass_t) {
2048 4025580 case 0:
2049 4025580 decode_sigpass(t1, width, height, bpno + 1, bandpos,
2050 vert_causal_ctx_csty_symbol);
2051 4025580 break;
2052 3688064 case 1:
2053 3688064 decode_refpass(t1, width, height, bpno + 1, vert_causal_ctx_csty_symbol);
2054 3688064 break;
2055 4399286 case 2:
2056 av_assert2(!t1->mqc.raw);
2057 4399286 decode_clnpass(s, t1, width, height, bpno + 1, bandpos,
2058 4399286 codsty->cblk_style & JPEG2000_CBLK_SEGSYM,
2059 vert_causal_ctx_csty_symbol);
2060 4399286 break;
2061 }
2062
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 12112930 times.
 12112930 if (codsty->cblk_style & JPEG2000_CBLK_RESET) // XXX no testcase for just this
2063 ✗ ff_mqc_init_contexts(&t1->mqc);
2064
2065
4/4
✓ Branch 0 taken 11301306 times.
✓ Branch 1 taken 811624 times.
✓ Branch 3 taken 3861 times.
✓ Branch 4 taken 11297445 times.
 12112930 if (passno && (coder_type = needs_termination(codsty->cblk_style, pass_cnt))) {
2066
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 3861 times.
 3861 if (term_cnt >= cblk->nb_terminations) {
2067 ✗ av_log(s->avctx, AV_LOG_ERROR, "Missing needed termination \n");
2068 ✗ return AVERROR_INVALIDDATA;
2069 }
2070
2/4
✓ Branch 0 taken 3861 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 3861 times.
 3861 if (FFABS(cblk->data + cblk->data_start[term_cnt + 1] - 2 - t1->mqc.bp) > 0) {
2071 ✗ av_log(s->avctx, AV_LOG_WARNING, "Mid mismatch %td in pass %d of %d\n",
2072 ✗ cblk->data + cblk->data_start[term_cnt + 1] - 2 - t1->mqc.bp,
2073 ✗ pass_cnt, cblk->npasses);
2074 }
2075
2076 3861 ff_mqc_initdec(&t1->mqc, cblk->data + cblk->data_start[++term_cnt], coder_type == 2, 0);
2077 }
2078
2079 12112930 pass_t++;
2080
2/2
✓ Branch 0 taken 4399286 times.
✓ Branch 1 taken 7713644 times.
 12112930 if (pass_t == 3) {
2081 4399286 bpno--;
2082 4399286 pass_t = 0;
2083 }
2084 12112930 pass_cnt ++;
2085 }
2086
2087
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 811624 times.
 811624 if (cblk->data + cblk->length - 2 > t1->mqc.bp) {
2088 ✗ av_log(s->avctx, AV_LOG_WARNING, "End mismatch %td\n",
2089 ✗ cblk->data + cblk->length - 2 - t1->mqc.bp);
2090 }
2091
2092
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 811624 times.
 811624 if (cblk->data + cblk->length < t1->mqc.bp) {
2093 ✗ av_log(s->avctx, AV_LOG_WARNING, "Synthetic End of Stream Marker Read.\n");
2094 }
2095
2096 /* Reconstruct the sample values */
2097
2/2
✓ Branch 0 taken 11964899 times.
✓ Branch 1 taken 811624 times.
 12776523 for (int y = 0; y < height; y++) {
2098
2/2
✓ Branch 0 taken 197511925 times.
✓ Branch 1 taken 11964899 times.
 209476824 for (int x = 0; x < width; x++) {
2099 int32_t sign, n, val;
2100 197511925 const uint32_t mask = (int64_t)UINT32_MAX >> (M_b + 1); // bit mask for ROI detection
2101
2102 197511925 n = x + (y * t1->stride);
2103 197511925 val = t1->data[n];
2104 197511925 sign = val & INT32_MIN;
2105 197511925 val &= INT32_MAX;
2106 /* ROI shift, if necessary */
2107
4/4
✓ Branch 0 taken 32768 times.
✓ Branch 1 taken 197479157 times.
✓ Branch 2 taken 30464 times.
✓ Branch 3 taken 2304 times.
 197511925 if (roi_shift && (((uint32_t)val & ~mask) == 0))
2108 30464 val = (uint32_t)val << roi_shift;
2109 197511925 t1->data[n] = val | sign; /* NOTE: Binary point for reconstruction value is located in 31 - M_b */
2110 }
2111 }
2112 811624 return 1;
2113 }
2114
2115 /* TODO: Verify dequantization for lossless case
2116 * comp->data can be float or int
2117 * band->stepsize can be float or int
2118 * depending on the type of DWT transformation.
2119 * see ISO/IEC 15444-1:2002 A.6.1 */
2120
2121 /* Float dequantization of a codeblock.*/
2122 ✗ static void dequantization_float(int x, int y, Jpeg2000Cblk *cblk,
2123 Jpeg2000Component *comp,
2124 Jpeg2000T1Context *t1, Jpeg2000Band *band, const int M_b)
2125 {
2126 int i, j;
2127 ✗ int w = cblk->coord[0][1] - cblk->coord[0][0];
2128 ✗ const int downshift = 31 - M_b;
2129 ✗ float fscale = band->f_stepsize;
2130 ✗ fscale /= (float)(1LL << downshift);
2131 ✗ for (j = 0; j < (cblk->coord[1][1] - cblk->coord[1][0]); ++j) {
2132 ✗ float *datap = &comp->f_data[(comp->coord[0][1] - comp->coord[0][0]) * (y + j) + x];
2133 ✗ int *src = t1->data + j*t1->stride;
2134 ✗ for (i = 0; i < w; ++i) {
2135 ✗ int val = src[i];
2136 ✗ if (val < 0) // Convert sign-magnitude to two's complement
2137 ✗ val = -(val & INT32_MAX);
2138 ✗ datap[i] = (float)val * fscale;
2139 }
2140 }
2141 ✗ }
2142
2143 /* Integer dequantization of a codeblock.*/
2144 606367 static void dequantization_int(int x, int y, Jpeg2000Cblk *cblk,
2145 Jpeg2000Component *comp,
2146 Jpeg2000T1Context *t1, Jpeg2000Band *band, const int M_b)
2147 {
2148 int i, j;
2149
1/2
✓ Branch 0 taken 606367 times.
✗ Branch 1 not taken.
 606367 const int downshift = FFMIN(31 - M_b, 31);
2150 606367 int w = cblk->coord[0][1] - cblk->coord[0][0];
2151
2/2
✓ Branch 0 taken 8882373 times.
✓ Branch 1 taken 606367 times.
 9488740 for (j = 0; j < (cblk->coord[1][1] - cblk->coord[1][0]); ++j) {
2152 8882373 int32_t *datap = &comp->i_data[(comp->coord[0][1] - comp->coord[0][0]) * (y + j) + x];
2153 8882373 int *src = t1->data + j*t1->stride;
2154
1/2
✓ Branch 0 taken 8882373 times.
✗ Branch 1 not taken.
 8882373 if (band->i_stepsize == 32768) {
2155
2/2
✓ Branch 0 taken 145656327 times.
✓ Branch 1 taken 8882373 times.
 154538700 for (i = 0; i < w; ++i) {
2156 145656327 int val = src[i];
2157
2/2
✓ Branch 0 taken 42982330 times.
✓ Branch 1 taken 102673997 times.
 145656327 if (val < 0) // Convert sign-magnitude to two's complement
2158 42982330 val = -((val & INT32_MAX) >> downshift);
2159 else
2160 102673997 val >>= downshift;
2161 145656327 datap[i] = val;
2162 }
2163 } else {
2164 // This should be VERY uncommon
2165 ✗ for (i = 0; i < w; ++i) {
2166 ✗ int val = src[i];
2167 ✗ if (val < 0) // Convert sign-magnitude to two's complement
2168 ✗ val = -((val & INT32_MAX) >> downshift);
2169 else
2170 ✗ val >>= downshift;
2171 ✗ datap[i] = (val * (int64_t)band->i_stepsize) / 65536;
2172 }
2173 }
2174 }
2175 606367 }
2176
2177 205267 static void dequantization_int_97(int x, int y, Jpeg2000Cblk *cblk,
2178 Jpeg2000Component *comp,
2179 Jpeg2000T1Context *t1, Jpeg2000Band *band, const int M_b)
2180 {
2181 int i, j;
2182 205267 int w = cblk->coord[0][1] - cblk->coord[0][0];
2183 205267 float fscale = band->f_stepsize;
2184 205267 const int downshift = 31 - M_b;
2185 205267 const int PRESCALE = 6; // At least 6 is required to pass the conformance tests in ISO/IEC 15444-4
2186 int scale;
2187
2188 205267 fscale /= (float)(1LL << downshift);
2189 205267 fscale *= (float)(1 << PRESCALE);
2190 205267 fscale *= (float)(1 << (16 + I_PRESHIFT));
2191 205267 scale = (int)(fscale + 0.5);
2192 205267 band->i_stepsize = scale;
2193
2/2
✓ Branch 0 taken 3082878 times.
✓ Branch 1 taken 205267 times.
 3288145 for (j = 0; j < (cblk->coord[1][1] - cblk->coord[1][0]); ++j) {
2194 3082878 int32_t *datap = &comp->i_data[(comp->coord[0][1] - comp->coord[0][0]) * (y + j) + x];
2195 3082878 int *src = t1->data + j*t1->stride;
2196
2/2
✓ Branch 0 taken 51871982 times.
✓ Branch 1 taken 3082878 times.
 54954860 for (i = 0; i < w; ++i) {
2197 51871982 int val = src[i];
2198
2/2
✓ Branch 0 taken 9002949 times.
✓ Branch 1 taken 42869033 times.
 51871982 if (val < 0) // Convert sign-magnitude to two's complement
2199 9002949 val = -(val & INT32_MAX);
2200 // Shifting down to prevent overflow in dequantization
2201 51871982 val = (val + (1LL << (PRESCALE - 1))) >> PRESCALE;
2202
2/2
✓ Branch 0 taken 9121136 times.
✓ Branch 1 taken 42750846 times.
 51871982 datap[i] = RSHIFT(val * (int64_t)band->i_stepsize, 16);
2203 }
2204 }
2205 205267 }
2206
2207 8 static inline void mct_decode(const Jpeg2000DecoderContext *s, Jpeg2000Tile *tile)
2208 {
2209 8 int i, csize = 1;
2210 void *src[3];
2211
2212
2/2
✓ Branch 0 taken 16 times.
✓ Branch 1 taken 8 times.
 24 for (i = 1; i < 3; i++) {
2213
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 16 times.
 16 if (tile->codsty[0].transform != tile->codsty[i].transform) {
2214 ✗ av_log(s->avctx, AV_LOG_ERROR, "Transforms mismatch, MCT not supported\n");
2215 ✗ return;
2216 }
2217
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 16 times.
 16 if (memcmp(tile->comp[0].coord, tile->comp[i].coord, sizeof(tile->comp[0].coord))) {
2218 ✗ av_log(s->avctx, AV_LOG_ERROR, "Coords mismatch, MCT not supported\n");
2219 ✗ return;
2220 }
2221 }
2222
2223
2/2
✓ Branch 0 taken 24 times.
✓ Branch 1 taken 8 times.
 32 for (i = 0; i < 3; i++)
2224
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 24 times.
 24 if (tile->codsty[0].transform == FF_DWT97)
2225 ✗ src[i] = tile->comp[i].f_data;
2226 else
2227 24 src[i] = tile->comp[i].i_data;
2228
2229
2/2
✓ Branch 0 taken 16 times.
✓ Branch 1 taken 8 times.
 24 for (i = 0; i < 2; i++)
2230 16 csize *= tile->comp[0].coord[i][1] - tile->comp[0].coord[i][0];
2231
2232 8 s->dsp.mct_decode[tile->codsty[0].transform](src[0], src[1], src[2], csize);
2233 }
2234
2235
2236 2881 static inline int tile_codeblocks(const Jpeg2000DecoderContext *s, Jpeg2000Tile *tile)
2237 {
2238 Jpeg2000T1Context t1;
2239
2240 int compno, reslevelno, bandno;
2241
2242 /* Loop on tile components */
2243
2/2
✓ Branch 0 taken 9264 times.
✓ Branch 1 taken 2881 times.
 12145 for (compno = 0; compno < s->ncomponents; compno++) {
2244 9264 Jpeg2000Component *comp = tile->comp + compno;
2245 9264 Jpeg2000CodingStyle *codsty = tile->codsty + compno;
2246 9264 Jpeg2000QuantStyle *quantsty = tile->qntsty + compno;
2247
2248 9264 int coded = 0;
2249 9264 int subbandno = 0;
2250
2251 9264 t1.stride = (1<<codsty->log2_cblk_width) + 2;
2252
2253 /* Loop on resolution levels */
2254
2/2
✓ Branch 0 taken 62434 times.
✓ Branch 1 taken 9264 times.
 71698 for (reslevelno = 0; reslevelno < codsty->nreslevels2decode; reslevelno++) {
2255 62434 Jpeg2000ResLevel *rlevel = comp->reslevel + reslevelno;
2256 /* Loop on bands */
2257
2/2
✓ Branch 0 taken 168774 times.
✓ Branch 1 taken 62434 times.
 231208 for (bandno = 0; bandno < rlevel->nbands; bandno++, subbandno++) {
2258 int nb_precincts, precno;
2259 168774 Jpeg2000Band *band = rlevel->band + bandno;
2260 168774 int cblkno = 0, bandpos;
2261 /* See Rec. ITU-T T.800, Equation E-2 */
2262 168774 int M_b = quantsty->expn[subbandno] + quantsty->nguardbits - 1;
2263
2264 168774 bandpos = bandno + (reslevelno > 0);
2265
2266
2/2
✓ Branch 0 taken 168772 times.
✓ Branch 1 taken 2 times.
 168774 if (band->coord[0][0] == band->coord[0][1] ||
2267
2/2
✓ Branch 0 taken 7800 times.
✓ Branch 1 taken 160972 times.
 168772 band->coord[1][0] == band->coord[1][1])
2268 7802 continue;
2269
2270
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 160972 times.
 160972 if (M_b > 31) {
2271 ✗ avpriv_request_sample(s->avctx, "M_b (%d) > 31", M_b);
2272 ✗ return AVERROR_PATCHWELCOME;
2273 }
2274
2275 160972 nb_precincts = rlevel->num_precincts_x * rlevel->num_precincts_y;
2276 /* Loop on precincts */
2277
2/2
✓ Branch 0 taken 162151 times.
✓ Branch 1 taken 160972 times.
 323123 for (precno = 0; precno < nb_precincts; precno++) {
2278 162151 Jpeg2000Prec *prec = band->prec + precno;
2279
2280 /* Loop on codeblocks */
2281 162151 for (cblkno = 0;
2282
2/2
✓ Branch 0 taken 919847 times.
✓ Branch 1 taken 162151 times.
 1081998 cblkno < prec->nb_codeblocks_width * prec->nb_codeblocks_height;
2283 919847 cblkno++) {
2284 int x, y, ret;
2285
2286 919847 Jpeg2000Cblk *cblk = prec->cblk + cblkno;
2287
2288
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 919837 times.
 919847 if (cblk->modes & JPEG2000_CTSY_HTJ2K_F)
2289 10 ret = ff_jpeg2000_decode_htj2k(s, codsty, &t1, cblk,
2290 10 cblk->coord[0][1] - cblk->coord[0][0],
2291 10 cblk->coord[1][1] - cblk->coord[1][0],
2292 10 M_b, comp->roi_shift);
2293 else
2294 919837 ret = decode_cblk(s, codsty, &t1, cblk,
2295 919837 cblk->coord[0][1] - cblk->coord[0][0],
2296 919837 cblk->coord[1][1] - cblk->coord[1][0],
2297 919837 bandpos, comp->roi_shift, M_b);
2298
2299
2/2
✓ Branch 0 taken 811634 times.
✓ Branch 1 taken 108213 times.
 919847 if (ret)
2300 811634 coded = 1;
2301 else
2302 108213 continue;
2303 811634 x = cblk->coord[0][0] - band->coord[0][0];
2304 811634 y = cblk->coord[1][0] - band->coord[1][0];
2305
2306
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 811634 times.
 811634 if (codsty->transform == FF_DWT97)
2307 ✗ dequantization_float(x, y, cblk, comp, &t1, band, M_b);
2308
2/2
✓ Branch 0 taken 205267 times.
✓ Branch 1 taken 606367 times.
 811634 else if (codsty->transform == FF_DWT97_INT)
2309 205267 dequantization_int_97(x, y, cblk, comp, &t1, band, M_b);
2310 else
2311 606367 dequantization_int(x, y, cblk, comp, &t1, band, M_b);
2312 } /* end cblk */
2313 } /*end prec */
2314 } /* end band */
2315 } /* end reslevel */
2316
2317 /* inverse DWT */
2318
1/2
✓ Branch 0 taken 9264 times.
✗ Branch 1 not taken.
 9264 if (coded)
2319
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 9264 times.
 9264 ff_dwt_decode(&comp->dwt, codsty->transform == FF_DWT97 ? (void*)comp->f_data : (void*)comp->i_data);
2320
2321 } /*end comp */
2322 2881 return 0;
2323 }
2324
2325 #define WRITE_FRAME(D, PIXEL) \
2326 static inline void write_frame_ ## D(const Jpeg2000DecoderContext * s, Jpeg2000Tile * tile, \
2327 AVFrame * picture, int precision) \
2328 { \
2329 const AVPixFmtDescriptor *pixdesc = av_pix_fmt_desc_get(s->avctx->pix_fmt); \
2330 int planar = !!(pixdesc->flags & AV_PIX_FMT_FLAG_PLANAR); \
2331 int pixelsize = planar ? 1 : pixdesc->nb_components; \
2332 \
2333 int compno; \
2334 int x, y; \
2335 \
2336 for (compno = 0; compno < s->ncomponents; compno++) { \
2337 Jpeg2000Component *comp = tile->comp + compno; \
2338 Jpeg2000CodingStyle *codsty = tile->codsty + compno; \
2339 PIXEL *line; \
2340 float *datap = comp->f_data; \
2341 int32_t *i_datap = comp->i_data; \
2342 int cbps = s->cbps[compno]; \
2343 int w = tile->comp[compno].coord[0][1] - \
2344 ff_jpeg2000_ceildiv(s->image_offset_x, s->cdx[compno]); \
2345 int h = tile->comp[compno].coord[1][1] - \
2346 ff_jpeg2000_ceildiv(s->image_offset_y, s->cdy[compno]); \
2347 int plane = 0; \
2348 ptrdiff_t dstoffset = 0; \
2349 \
2350 if (planar) \
2351 plane = s->cdef[compno] ? s->cdef[compno]-1 : (s->ncomponents-1); \
2352 else \
2353 dstoffset = s->cdef[compno] ? s->cdef[compno] - 1 : compno; \
2354 \
2355 y = tile->comp[compno].coord[1][0] - \
2356 ff_jpeg2000_ceildiv(s->image_offset_y, s->cdy[compno]); \
2357 line = (PIXEL *)picture->data[plane] + y * (picture->linesize[plane] / sizeof(PIXEL));\
2358 for (; y < h; y++) { \
2359 PIXEL *dst; \
2360 \
2361 x = tile->comp[compno].coord[0][0] - \
2362 ff_jpeg2000_ceildiv(s->image_offset_x, s->cdx[compno]); \
2363 dst = line + x * pixelsize + dstoffset; \
2364 \
2365 if (codsty->transform == FF_DWT97) { \
2366 for (; x < w; x++) { \
2367 int val = lrintf(*datap) + (1 << (cbps - 1)); \
2368 /* DC level shift and clip see ISO 15444-1:2002 G.1.2 */ \
2369 val = av_clip(val, 0, (1 << cbps) - 1); \
2370 *dst = val << (precision - cbps); \
2371 datap++; \
2372 dst += pixelsize; \
2373 } \
2374 } else { \
2375 for (; x < w; x++) { \
2376 int val = *i_datap + (1 << (cbps - 1)); \
2377 /* DC level shift and clip see ISO 15444-1:2002 G.1.2 */ \
2378 val = av_clip(val, 0, (1 << cbps) - 1); \
2379 *dst = val << (precision - cbps); \
2380 i_datap++; \
2381 dst += pixelsize; \
2382 } \
2383 } \
2384 line += picture->linesize[plane] / sizeof(PIXEL); \
2385 } \
2386 } \
2387 \
2388 }
2389
2390
14/18
✓ Branch 1 taken 1321 times.
✓ Branch 2 taken 1 times.
✓ Branch 5 taken 4 times.
✓ Branch 6 taken 3933 times.
✓ Branch 7 taken 3 times.
✓ Branch 8 taken 1 times.
✓ Branch 9 taken 3933 times.
✗ Branch 10 not taken.
✗ Branch 13 not taken.
✓ Branch 14 taken 535220 times.
✗ Branch 15 not taken.
✗ Branch 16 not taken.
✓ Branch 17 taken 95336791 times.
✓ Branch 18 taken 535220 times.
✓ Branch 19 taken 535220 times.
✓ Branch 20 taken 3937 times.
✓ Branch 21 taken 3937 times.
✓ Branch 22 taken 1322 times.
 95877270 WRITE_FRAME(8, uint8_t)
2391
14/18
✓ Branch 1 taken 909 times.
✓ Branch 2 taken 650 times.
✓ Branch 5 taken 2600 times.
✓ Branch 6 taken 2727 times.
✓ Branch 7 taken 1950 times.
✓ Branch 8 taken 650 times.
✓ Branch 9 taken 2727 times.
✗ Branch 10 not taken.
✗ Branch 13 not taken.
✓ Branch 14 taken 726092 times.
✗ Branch 15 not taken.
✗ Branch 16 not taken.
✓ Branch 17 taken 133058168 times.
✓ Branch 18 taken 726092 times.
✓ Branch 19 taken 726092 times.
✓ Branch 20 taken 5327 times.
✓ Branch 21 taken 5327 times.
✓ Branch 22 taken 1559 times.
 133791146 WRITE_FRAME(16, uint16_t)
2392
2393 #undef WRITE_FRAME
2394
2395 2881 static int jpeg2000_decode_tile(AVCodecContext *avctx, void *td,
2396 int jobnr, int threadnr)
2397 {
2398 2881 const Jpeg2000DecoderContext *s = avctx->priv_data;
2399 2881 AVFrame *picture = td;
2400 2881 Jpeg2000Tile *tile = s->tile + jobnr;
2401
2402 2881 int ret = tile_codeblocks(s, tile);
2403
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2881 times.
 2881 if (ret < 0)
2404 ✗ return ret;
2405
2406 /* inverse MCT transformation */
2407
2/2
✓ Branch 0 taken 8 times.
✓ Branch 1 taken 2873 times.
 2881 if (tile->codsty[0].mct)
2408 8 mct_decode(s, tile);
2409
2410
2/2
✓ Branch 0 taken 1322 times.
✓ Branch 1 taken 1559 times.
 2881 if (s->precision <= 8) {
2411 1322 write_frame_8(s, tile, picture, 8);
2412 } else {
2413 4675 int precision = picture->format == AV_PIX_FMT_XYZ12 ||
2414
2/2
✓ Branch 0 taken 650 times.
✓ Branch 1 taken 907 times.
 1557 picture->format == AV_PIX_FMT_RGB48 ||
2415
1/2
✓ Branch 0 taken 650 times.
✗ Branch 1 not taken.
 650 picture->format == AV_PIX_FMT_RGBA64 ||
2416
3/4
✓ Branch 0 taken 1557 times.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 650 times.
✗ Branch 3 not taken.
 3116 picture->format == AV_PIX_FMT_GRAY16 ? 16 : s->precision;
2417
2418 1559 write_frame_16(s, tile, picture, precision);
2419 }
2420
2421 2881 return 0;
2422 }
2423
2424 851 static void jpeg2000_dec_cleanup(Jpeg2000DecoderContext *s)
2425 {
2426 int tileno, compno;
2427
2/2
✓ Branch 0 taken 3216 times.
✓ Branch 1 taken 851 times.
 4067 for (tileno = 0; tileno < s->numXtiles * s->numYtiles; tileno++) {
2428
1/2
✓ Branch 0 taken 3216 times.
✗ Branch 1 not taken.
 3216 if (s->tile[tileno].comp) {
2429
2/2
✓ Branch 0 taken 10255 times.
✓ Branch 1 taken 3216 times.
 13471 for (compno = 0; compno < s->ncomponents; compno++) {
2430 10255 Jpeg2000Component *comp = s->tile[tileno].comp + compno;
2431 10255 Jpeg2000CodingStyle *codsty = s->tile[tileno].codsty + compno;
2432
2433 10255 ff_jpeg2000_cleanup(comp, codsty);
2434 }
2435 3216 av_freep(&s->tile[tileno].comp);
2436 3216 av_freep(&s->tile[tileno].packed_headers);
2437 3216 s->tile[tileno].packed_headers_size = 0;
2438 }
2439 }
2440 851 av_freep(&s->packed_headers);
2441 851 s->packed_headers_size = 0;
2442 851 memset(&s->packed_headers_stream, 0, sizeof(s->packed_headers_stream));
2443 851 av_freep(&s->tile);
2444 851 memset(s->codsty, 0, sizeof(s->codsty));
2445 851 memset(s->qntsty, 0, sizeof(s->qntsty));
2446 851 memset(s->properties, 0, sizeof(s->properties));
2447 851 memset(&s->poc , 0, sizeof(s->poc));
2448 851 memset(s->roi_shift, 0, sizeof(s->roi_shift));
2449 851 s->numXtiles = s->numYtiles = 0;
2450 851 s->ncomponents = 0;
2451 851 s->has_ppm = 0;
2452 851 s->isHT = 0;
2453 851 s->precision = 0;
2454 851 s->colour_space = 0;
2455 851 s->pal8 = 0;
2456 851 }
2457
2458 851 static int jpeg2000_read_main_headers(Jpeg2000DecoderContext *s)
2459 {
2460 851 Jpeg2000CodingStyle *codsty = s->codsty;
2461 851 Jpeg2000QuantStyle *qntsty = s->qntsty;
2462 851 Jpeg2000POC *poc = &s->poc;
2463 851 uint8_t *properties = s->properties;
2464 851 uint8_t in_tile_headers = 0;
2465
2466 9671 for (;;) {
2467 10522 int len, ret = 0;
2468 uint16_t marker;
2469 int oldpos;
2470
2471
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 10522 times.
 10522 if (bytestream2_get_bytes_left(&s->g) < 2) {
2472 ✗ av_log(s->avctx, AV_LOG_ERROR, "Missing EOC\n");
2473 ✗ break;
2474 }
2475
2476 10522 marker = bytestream2_get_be16u(&s->g);
2477 10522 oldpos = bytestream2_tell(&s->g);
2478
3/4
✓ Branch 0 taken 10522 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 10520 times.
 10522 if (marker >= 0xFF30 && marker <= 0xFF3F)
2479 2 continue;
2480
2/2
✓ Branch 0 taken 3242 times.
✓ Branch 1 taken 7278 times.
 10520 if (marker == JPEG2000_SOD) {
2481 Jpeg2000Tile *tile;
2482 Jpeg2000TilePart *tp;
2483
2484
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 3242 times.
 3242 if (!s->tile) {
2485 ✗ av_log(s->avctx, AV_LOG_ERROR, "Missing SIZ\n");
2486 ✗ return AVERROR_INVALIDDATA;
2487 }
2488
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 3242 times.
 3242 if (s->curtileno < 0) {
2489 ✗ av_log(s->avctx, AV_LOG_ERROR, "Missing SOT\n");
2490 ✗ return AVERROR_INVALIDDATA;
2491 }
2492
2493 3242 tile = s->tile + s->curtileno;
2494 3242 tp = tile->tile_part + tile->tp_idx;
2495
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 3242 times.
 3242 if (tp->tp_end < s->g.buffer) {
2496 ✗ av_log(s->avctx, AV_LOG_ERROR, "Invalid tpend\n");
2497 ✗ return AVERROR_INVALIDDATA;
2498 }
2499
2500
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 3242 times.
 3242 if (s->has_ppm) {
2501 ✗ uint32_t tp_header_size = bytestream2_get_be32(&s->packed_headers_stream);
2502 ✗ if (bytestream2_get_bytes_left(&s->packed_headers_stream) < tp_header_size)
2503 ✗ return AVERROR_INVALIDDATA;
2504 ✗ bytestream2_init(&tp->header_tpg, s->packed_headers_stream.buffer, tp_header_size);
2505 ✗ bytestream2_skip(&s->packed_headers_stream, tp_header_size);
2506 }
2507
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 3242 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
 3242 if (tile->has_ppt && tile->tp_idx == 0) {
2508 ✗ bytestream2_init(&tile->packed_headers_stream, tile->packed_headers, tile->packed_headers_size);
2509 }
2510
2511 3242 bytestream2_init(&tp->tpg, s->g.buffer, tp->tp_end - s->g.buffer);
2512 3242 bytestream2_skip(&s->g, tp->tp_end - s->g.buffer);
2513
2514 3242 continue;
2515 }
2516
2/2
✓ Branch 0 taken 851 times.
✓ Branch 1 taken 6427 times.
 7278 if (marker == JPEG2000_EOC)
2517 851 break;
2518
2519 6427 len = bytestream2_get_be16(&s->g);
2520
2/4
✓ Branch 0 taken 6427 times.
✗ Branch 1 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 6427 times.
 6427 if (len < 2 || bytestream2_get_bytes_left(&s->g) < len - 2) {
2521 ✗ if (s->avctx->strict_std_compliance >= FF_COMPLIANCE_STRICT) {
2522 ✗ av_log(s->avctx, AV_LOG_ERROR, "Invalid len %d left=%d\n", len, bytestream2_get_bytes_left(&s->g));
2523 ✗ return AVERROR_INVALIDDATA;
2524 }
2525 ✗ av_log(s->avctx, AV_LOG_WARNING, "Missing EOC Marker.\n");
2526 ✗ break;
2527 }
2528
2529
14/17
✓ Branch 0 taken 851 times.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 18 times.
✓ Branch 3 taken 851 times.
✓ Branch 4 taken 4 times.
✓ Branch 5 taken 25 times.
✓ Branch 6 taken 851 times.
✓ Branch 7 taken 8 times.
✓ Branch 8 taken 3242 times.
✓ Branch 9 taken 44 times.
✓ Branch 10 taken 4 times.
✓ Branch 11 taken 11 times.
✓ Branch 12 taken 514 times.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✗ Branch 16 not taken.
 6427 switch (marker) {
2530 851 case JPEG2000_SIZ:
2531
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 851 times.
 851 if (s->ncomponents) {
2532 ✗ av_log(s->avctx, AV_LOG_ERROR, "Duplicate SIZ\n");
2533 ✗ return AVERROR_INVALIDDATA;
2534 }
2535 851 ret = get_siz(s);
2536
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 851 times.
 851 if (!s->tile)
2537 ✗ s->numXtiles = s->numYtiles = 0;
2538 851 break;
2539 2 case JPEG2000_CAP:
2540
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
 2 if (!s->ncomponents) {
2541 ✗ av_log(s->avctx, AV_LOG_ERROR, "CAP marker segment shall come after SIZ\n");
2542 ✗ return AVERROR_INVALIDDATA;
2543 }
2544 2 ret = get_cap(s, codsty);
2545 2 break;
2546 18 case JPEG2000_COC:
2547
1/6
✗ Branch 0 not taken.
✓ Branch 1 taken 18 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
 18 if (in_tile_headers == 1 && s->isHT && (!s->Ccap15_b11)) {
2548 ✗ av_log(s->avctx, AV_LOG_ERROR,
2549 "COC marker found in a tile header but the codestream belongs to the HOMOGENEOUS set\n");
2550 ✗ return AVERROR_INVALIDDATA;
2551 }
2552 18 ret = get_coc(s, codsty, properties);
2553 18 break;
2554 851 case JPEG2000_COD:
2555
1/6
✗ Branch 0 not taken.
✓ Branch 1 taken 851 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
 851 if (in_tile_headers == 1 && s->isHT && (!s->Ccap15_b11)) {
2556 ✗ av_log(s->avctx, AV_LOG_ERROR,
2557 "COD marker found in a tile header but the codestream belongs to the HOMOGENEOUS set\n");
2558 ✗ return AVERROR_INVALIDDATA;
2559 }
2560 851 ret = get_cod(s, codsty, properties);
2561 851 break;
2562 4 case JPEG2000_RGN:
2563
2/6
✓ Branch 0 taken 4 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 4 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
 4 if (in_tile_headers == 1 && s->isHT && (!s->Ccap15_b11)) {
2564 ✗ av_log(s->avctx, AV_LOG_ERROR,
2565 "RGN marker found in a tile header but the codestream belongs to the HOMOGENEOUS set\n");
2566 ✗ return AVERROR_INVALIDDATA;
2567 }
2568 4 ret = get_rgn(s, len);
2569
2/4
✓ Branch 0 taken 4 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 4 times.
 4 if ((!s->Ccap15_b12) && s->isHT) {
2570 ✗ av_log(s->avctx, AV_LOG_ERROR, "RGN marker found but the codestream belongs to the RGNFREE set\n");
2571 ✗ return AVERROR_INVALIDDATA;
2572 }
2573 4 break;
2574 25 case JPEG2000_QCC:
2575
1/6
✗ Branch 0 not taken.
✓ Branch 1 taken 25 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
 25 if (in_tile_headers == 1 && s->isHT && (!s->Ccap15_b11)) {
2576 ✗ av_log(s->avctx, AV_LOG_ERROR,
2577 "QCC marker found in a tile header but the codestream belongs to the HOMOGENEOUS set\n");
2578 ✗ return AVERROR_INVALIDDATA;
2579 }
2580 25 ret = get_qcc(s, len, qntsty, properties);
2581 25 break;
2582 851 case JPEG2000_QCD:
2583
1/6
✗ Branch 0 not taken.
✓ Branch 1 taken 851 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
 851 if (in_tile_headers == 1 && s->isHT && (!s->Ccap15_b11)) {
2584 ✗ av_log(s->avctx, AV_LOG_ERROR,
2585 "QCD marker found in a tile header but the codestream belongs to the HOMOGENEOUS set\n");
2586 ✗ return AVERROR_INVALIDDATA;
2587 }
2588 851 ret = get_qcd(s, len, qntsty, properties);
2589 851 break;
2590 8 case JPEG2000_POC:
2591
3/6
✓ Branch 0 taken 4 times.
✓ Branch 1 taken 4 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 4 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
 8 if (in_tile_headers == 1 && s->isHT && (!s->Ccap15_b11)) {
2592 ✗ av_log(s->avctx, AV_LOG_ERROR,
2593 "POC marker found in a tile header but the codestream belongs to the HOMOGENEOUS set\n");
2594 ✗ return AVERROR_INVALIDDATA;
2595 }
2596 8 ret = get_poc(s, len, poc);
2597 8 break;
2598 3242 case JPEG2000_SOT:
2599
2/2
✓ Branch 0 taken 851 times.
✓ Branch 1 taken 2391 times.
 3242 if (!in_tile_headers) {
2600 851 in_tile_headers = 1;
2601
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 851 times.
 851 if (s->has_ppm) {
2602 ✗ bytestream2_init(&s->packed_headers_stream, s->packed_headers, s->packed_headers_size);
2603 }
2604 }
2605
1/2
✓ Branch 1 taken 3242 times.
✗ Branch 2 not taken.
 3242 if (!(ret = get_sot(s, len))) {
2606 av_assert1(s->curtileno >= 0);
2607 3242 codsty = s->tile[s->curtileno].codsty;
2608 3242 qntsty = s->tile[s->curtileno].qntsty;
2609 3242 poc = &s->tile[s->curtileno].poc;
2610 3242 properties = s->tile[s->curtileno].properties;
2611 }
2612 3242 break;
2613 44 case JPEG2000_PLM:
2614 // the PLM marker is ignored
2615 case JPEG2000_COM:
2616 // the comment is ignored
2617 44 bytestream2_skip(&s->g, len - 2);
2618 44 break;
2619 4 case JPEG2000_CRG:
2620 4 ret = read_crg(s, len);
2621 4 break;
2622 11 case JPEG2000_TLM:
2623 // Tile-part lengths
2624 11 ret = get_tlm(s, len);
2625 11 break;
2626 514 case JPEG2000_PLT:
2627 // Packet length, tile-part header
2628 514 ret = get_plt(s, len);
2629 514 break;
2630 ✗ case JPEG2000_PPM:
2631 // Packed headers, main header
2632 ✗ if (in_tile_headers) {
2633 ✗ av_log(s->avctx, AV_LOG_ERROR, "PPM Marker can only be in Main header\n");
2634 ✗ return AVERROR_INVALIDDATA;
2635 }
2636 ✗ ret = get_ppm(s, len);
2637 ✗ break;
2638 ✗ case JPEG2000_PPT:
2639 // Packed headers, tile-part header
2640 ✗ if (s->has_ppm) {
2641 ✗ av_log(s->avctx, AV_LOG_ERROR,
2642 "Cannot have both PPT and PPM marker.\n");
2643 ✗ return AVERROR_INVALIDDATA;
2644 }
2645 ✗ if ((!s->Ccap15_b11) && s->isHT) {
2646 ✗ av_log(s->avctx, AV_LOG_ERROR, "PPT marker found but the codestream belongs to the HOMOGENEOUS set\n");
2647 ✗ return AVERROR_INVALIDDATA;
2648 }
2649 ✗ ret = get_ppt(s, len);
2650 ✗ break;
2651 2 case JPEG2000_CPF:
2652 // Corresponding profile marker
2653 2 ret = read_cpf(s, len);
2654 2 break;
2655 ✗ default:
2656 ✗ av_log(s->avctx, AV_LOG_ERROR,
2657 "unsupported marker 0x%.4"PRIX16" at pos 0x%X\n",
2658 ✗ marker, bytestream2_tell(&s->g) - 4);
2659 ✗ bytestream2_skip(&s->g, len - 2);
2660 ✗ break;
2661 }
2662
2/4
✓ Branch 1 taken 6427 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 6427 times.
 6427 if (bytestream2_tell(&s->g) - oldpos != len || ret) {
2663 ✗ av_log(s->avctx, AV_LOG_ERROR,
2664 "error during processing marker segment %.4"PRIx16"\n",
2665 marker);
2666 ✗ return ret ? ret : -1;
2667 }
2668 }
2669 851 return 0;
2670 }
2671
2672 /* Read bit stream packets --> T2 operation. */
2673 817 static int jpeg2000_read_bitstream_packets(Jpeg2000DecoderContext *s)
2674 {
2675 817 int ret = 0;
2676 int tileno;
2677
2678
2/2
✓ Branch 0 taken 2881 times.
✓ Branch 1 taken 817 times.
 3698 for (tileno = 0; tileno < s->numXtiles * s->numYtiles; tileno++) {
2679 2881 Jpeg2000Tile *tile = s->tile + tileno;
2680
2681
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 2881 times.
 2881 if ((ret = init_tile(s, tileno)) < 0)
2682 ✗ return ret;
2683
2684
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 2881 times.
 2881 if ((ret = jpeg2000_decode_packets(s, tile)) < 0)
2685 ✗ return ret;
2686 }
2687
2688 817 return 0;
2689 }
2690
2691 817 static int jp2_find_codestream(Jpeg2000DecoderContext *s)
2692 {
2693 uint32_t atom_size, atom, atom_end;
2694 817 int search_range = 10;
2695
2696 817 while (search_range
2697
2/4
✓ Branch 0 taken 2452 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2452 times.
✗ Branch 3 not taken.
 4904 &&
2698 2452 bytestream2_get_bytes_left(&s->g) >= 8) {
2699 2452 atom_size = bytestream2_get_be32u(&s->g);
2700 2452 atom = bytestream2_get_be32u(&s->g);
2701
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2452 times.
 2452 if (atom_size == 1) {
2702 ✗ if (bytestream2_get_be32u(&s->g)) {
2703 ✗ avpriv_request_sample(s->avctx, "Huge atom");
2704 ✗ return 0;
2705 }
2706 ✗ atom_size = bytestream2_get_be32u(&s->g);
2707 ✗ if (atom_size < 16 || (int64_t)bytestream2_tell(&s->g) + atom_size - 16 > INT_MAX)
2708 ✗ return AVERROR_INVALIDDATA;
2709 ✗ atom_end = bytestream2_tell(&s->g) + atom_size - 16;
2710 } else {
2711
3/4
✓ Branch 0 taken 2451 times.
✓ Branch 1 taken 1 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 2451 times.
 2452 if (atom_size < 8 || (int64_t)bytestream2_tell(&s->g) + atom_size - 8 > INT_MAX)
2712 1 return AVERROR_INVALIDDATA;
2713 2451 atom_end = bytestream2_tell(&s->g) + atom_size - 8;
2714 }
2715
2716
2/2
✓ Branch 0 taken 816 times.
✓ Branch 1 taken 1635 times.
 2451 if (atom == JP2_CODESTREAM)
2717 816 return 1;
2718
2719
2/4
✓ Branch 1 taken 1635 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 1635 times.
 1635 if (bytestream2_get_bytes_left(&s->g) < atom_size || atom_end < atom_size)
2720 ✗ return 0;
2721
2722
3/4
✓ Branch 0 taken 817 times.
✓ Branch 1 taken 818 times.
✓ Branch 2 taken 817 times.
✗ Branch 3 not taken.
 1635 if (atom == JP2_HEADER &&
2723 817 atom_size >= 16) {
2724 uint32_t atom2_size, atom2, atom2_end;
2725 do {
2726
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 1634 times.
 1634 if (bytestream2_get_bytes_left(&s->g) < 8)
2727 ✗ break;
2728 1634 atom2_size = bytestream2_get_be32u(&s->g);
2729 1634 atom2 = bytestream2_get_be32u(&s->g);
2730 1634 atom2_end = bytestream2_tell(&s->g) + atom2_size - 8;
2731
3/6
✓ Branch 0 taken 1634 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1634 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1634 times.
✗ Branch 5 not taken.
 1634 if (atom2_size < 8 || atom2_end > atom_end || atom2_end < atom2_size)
2732 break;
2733 1634 atom2_size -= 8;
2734
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1634 times.
 1634 if (atom2 == JP2_CODESTREAM) {
2735 ✗ return 1;
2736
3/4
✓ Branch 0 taken 817 times.
✓ Branch 1 taken 817 times.
✓ Branch 2 taken 817 times.
✗ Branch 3 not taken.
 2451 } else if (atom2 == MKBETAG('c','o','l','r') && atom2_size >= 7) {
2737 817 int method = bytestream2_get_byteu(&s->g);
2738 817 bytestream2_skipu(&s->g, 2);
2739
1/2
✓ Branch 0 taken 817 times.
✗ Branch 1 not taken.
 817 if (method == 1) {
2740 817 s->colour_space = bytestream2_get_be32u(&s->g);
2741 }
2742
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 817 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
 817 } else if (atom2 == MKBETAG('p','c','l','r') && atom2_size >= 6) {
2743 int i, size, colour_count, colour_channels, colour_depth[3];
2744 ✗ colour_count = bytestream2_get_be16u(&s->g);
2745 ✗ colour_channels = bytestream2_get_byteu(&s->g);
2746 // FIXME: Do not ignore channel_sign
2747 ✗ colour_depth[0] = (bytestream2_get_byteu(&s->g) & 0x7f) + 1;
2748 ✗ colour_depth[1] = (bytestream2_get_byteu(&s->g) & 0x7f) + 1;
2749 ✗ colour_depth[2] = (bytestream2_get_byteu(&s->g) & 0x7f) + 1;
2750 ✗ size = (colour_depth[0] + 7 >> 3) * colour_count +
2751 ✗ (colour_depth[1] + 7 >> 3) * colour_count +
2752 ✗ (colour_depth[2] + 7 >> 3) * colour_count;
2753 ✗ if (colour_count > AVPALETTE_COUNT ||
2754 ✗ colour_channels != 3 ||
2755 ✗ colour_depth[0] > 16 ||
2756 ✗ colour_depth[1] > 16 ||
2757 ✗ colour_depth[2] > 16 ||
2758 ✗ atom2_size < size) {
2759 ✗ avpriv_request_sample(s->avctx, "Unknown palette");
2760 ✗ bytestream2_seek(&s->g, atom2_end, SEEK_SET);
2761 ✗ continue;
2762 }
2763 ✗ s->pal8 = 1;
2764 ✗ for (i = 0; i < colour_count; i++) {
2765 uint32_t r, g, b;
2766 ✗ if (colour_depth[0] <= 8) {
2767 ✗ r = bytestream2_get_byteu(&s->g) << 8 - colour_depth[0];
2768 ✗ r |= r >> colour_depth[0];
2769 } else {
2770 ✗ r = bytestream2_get_be16u(&s->g) >> colour_depth[0] - 8;
2771 }
2772 ✗ if (colour_depth[1] <= 8) {
2773 ✗ g = bytestream2_get_byteu(&s->g) << 8 - colour_depth[1];
2774 ✗ g |= g >> colour_depth[1];
2775 } else {
2776 ✗ g = bytestream2_get_be16u(&s->g) >> colour_depth[1] - 8;
2777 }
2778 ✗ if (colour_depth[2] <= 8) {
2779 ✗ b = bytestream2_get_byteu(&s->g) << 8 - colour_depth[2];
2780 ✗ b |= b >> colour_depth[2];
2781 } else {
2782 ✗ b = bytestream2_get_be16u(&s->g) >> colour_depth[2] - 8;
2783 }
2784 ✗ s->palette[i] = 0xffu << 24 | r << 16 | g << 8 | b;
2785 }
2786
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 817 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
 817 } else if (atom2 == MKBETAG('c','d','e','f') && atom2_size >= 2) {
2787 ✗ int n = bytestream2_get_be16u(&s->g);
2788 ✗ for (; n>0; n--) {
2789 ✗ int cn = bytestream2_get_be16(&s->g);
2790 ✗ av_unused int typ = bytestream2_get_be16(&s->g);
2791 ✗ int asoc = bytestream2_get_be16(&s->g);
2792 ✗ if (cn < 4 && asoc < 4)
2793 ✗ s->cdef[cn] = asoc;
2794 }
2795
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 817 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
 817 } else if (atom2 == MKBETAG('r','e','s',' ') && atom2_size >= 18) {
2796 int64_t vnum, vden, hnum, hden, vexp, hexp;
2797 uint32_t resx;
2798 ✗ bytestream2_skip(&s->g, 4);
2799 ✗ resx = bytestream2_get_be32u(&s->g);
2800 ✗ if (resx != MKBETAG('r','e','s','c') && resx != MKBETAG('r','e','s','d')) {
2801 ✗ bytestream2_seek(&s->g, atom2_end, SEEK_SET);
2802 ✗ continue;
2803 }
2804 ✗ vnum = bytestream2_get_be16u(&s->g);
2805 ✗ vden = bytestream2_get_be16u(&s->g);
2806 ✗ hnum = bytestream2_get_be16u(&s->g);
2807 ✗ hden = bytestream2_get_be16u(&s->g);
2808 ✗ vexp = bytestream2_get_byteu(&s->g);
2809 ✗ hexp = bytestream2_get_byteu(&s->g);
2810 ✗ if (!vnum || !vden || !hnum || !hden) {
2811 ✗ bytestream2_seek(&s->g, atom2_end, SEEK_SET);
2812 ✗ av_log(s->avctx, AV_LOG_WARNING, "RES box invalid\n");
2813 ✗ continue;
2814 }
2815 ✗ if (vexp > hexp) {
2816 ✗ vexp -= hexp;
2817 ✗ hexp = 0;
2818 } else {
2819 ✗ hexp -= vexp;
2820 ✗ vexp = 0;
2821 }
2822 ✗ if ( INT64_MAX / (hnum * vden) > pow(10, hexp)
2823 ✗ && INT64_MAX / (vnum * hden) > pow(10, vexp))
2824 ✗ av_reduce(&s->sar.den, &s->sar.num,
2825 ✗ hnum * vden * pow(10, hexp),
2826 ✗ vnum * hden * pow(10, vexp),
2827 INT32_MAX);
2828 }
2829 1634 bytestream2_seek(&s->g, atom2_end, SEEK_SET);
2830
2/2
✓ Branch 0 taken 817 times.
✓ Branch 1 taken 817 times.
 1634 } while (atom_end - atom2_end >= 8);
2831 } else {
2832 818 search_range--;
2833 }
2834 1635 bytestream2_seek(&s->g, atom_end, SEEK_SET);
2835 }
2836
2837 ✗ return 0;
2838 }
2839
2840 67 static av_cold int jpeg2000_decode_init(AVCodecContext *avctx)
2841 {
2842 67 Jpeg2000DecoderContext *s = avctx->priv_data;
2843
2844
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 67 times.
 67 if (avctx->lowres)
2845 ✗ av_log(avctx, AV_LOG_WARNING, "lowres is overridden by reduction_factor but set anyway\n");
2846
3/4
✓ Branch 0 taken 66 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 66 times.
✗ Branch 3 not taken.
 67 if (!s->reduction_factor && avctx->lowres < JPEG2000_MAX_RESLEVELS) {
2847 66 s->reduction_factor = avctx->lowres;
2848 }
2849
3/4
✓ Branch 0 taken 1 times.
✓ Branch 1 taken 66 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 1 times.
 67 if (avctx->lowres != s->reduction_factor && avctx->lowres)
2850 ✗ return AVERROR(EINVAL);
2851
2852 67 ff_jpeg2000dsp_init(&s->dsp);
2853 67 ff_jpeg2000_init_tier1_luts();
2854
2855 67 return 0;
2856 }
2857
2858 851 static int jpeg2000_decode_frame(AVCodecContext *avctx, AVFrame *picture,
2859 const AVPacket *avpkt)
2860 {
2861 851 Jpeg2000DecoderContext *s = avctx->priv_data;
2862 int ret;
2863
2864 851 s->avctx = avctx;
2865 851 bytestream2_init(&s->g, avpkt->data, avpkt->size);
2866 851 s->curtileno = -1;
2867 851 memset(s->cdef, -1, sizeof(s->cdef));
2868
2869
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 851 times.
 851 if (bytestream2_get_bytes_left(&s->g) < 2) {
2870 ✗ ret = AVERROR_INVALIDDATA;
2871 ✗ goto end;
2872 }
2873
2874 // check if the image is in jp2 format
2875
3/4
✓ Branch 1 taken 851 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 817 times.
✓ Branch 4 taken 34 times.
 1702 if (bytestream2_get_bytes_left(&s->g) >= 12 &&
2876
1/2
✓ Branch 1 taken 817 times.
✗ Branch 2 not taken.
 1668 (bytestream2_get_be32u(&s->g) == 12) &&
2877
1/2
✓ Branch 1 taken 817 times.
✗ Branch 2 not taken.
 1634 (bytestream2_get_be32u(&s->g) == JP2_SIG_TYPE) &&
2878 817 (bytestream2_get_be32u(&s->g) == JP2_SIG_VALUE)) {
2879
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 817 times.
 817 if (!jp2_find_codestream(s)) {
2880 ✗ av_log(avctx, AV_LOG_ERROR,
2881 "Could not find Jpeg2000 codestream atom.\n");
2882 ✗ ret = AVERROR_INVALIDDATA;
2883 ✗ goto end;
2884 }
2885 } else {
2886 34 bytestream2_seek(&s->g, 0, SEEK_SET);
2887 }
2888
2889
2/4
✓ Branch 1 taken 851 times.
✗ Branch 2 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 851 times.
 851 while (bytestream2_get_bytes_left(&s->g) >= 3 && bytestream2_peek_be16(&s->g) != JPEG2000_SOC)
2890 ✗ bytestream2_skip(&s->g, 1);
2891
2892
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 851 times.
 851 if (bytestream2_get_be16u(&s->g) != JPEG2000_SOC) {
2893 ✗ av_log(avctx, AV_LOG_ERROR, "SOC marker not present\n");
2894 ✗ ret = AVERROR_INVALIDDATA;
2895 ✗ goto end;
2896 }
2897
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 851 times.
 851 if (ret = jpeg2000_read_main_headers(s))
2898 ✗ goto end;
2899
2900
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 851 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
 851 if (s->sar.num && s->sar.den)
2901 ✗ avctx->sample_aspect_ratio = s->sar;
2902 851 s->sar.num = s->sar.den = 0;
2903
2904
2/2
✓ Branch 0 taken 34 times.
✓ Branch 1 taken 817 times.
 851 if (avctx->skip_frame >= AVDISCARD_ALL) {
2905 34 jpeg2000_dec_cleanup(s);
2906 34 return 0;
2907 }
2908
2909 /* get picture buffer */
2910
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 817 times.
 817 if ((ret = ff_thread_get_buffer(avctx, picture, 0)) < 0)
2911 ✗ goto end;
2912
2913
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 817 times.
 817 if (ret = jpeg2000_read_bitstream_packets(s))
2914 ✗ goto end;
2915
2916
4/4
✓ Branch 0 taken 2225 times.
✓ Branch 1 taken 612 times.
✓ Branch 2 taken 2020 times.
✓ Branch 3 taken 205 times.
 2837 for (int x = 0; x < s->ncomponents && s->codsty[x].transform == FF_DWT53;)
2917
2/2
✓ Branch 0 taken 612 times.
✓ Branch 1 taken 1408 times.
 2020 if (++x == s->ncomponents)
2918 612 picture->flags |= AV_FRAME_FLAG_LOSSLESS;
2919
2920 817 avctx->execute2(avctx, jpeg2000_decode_tile, picture, NULL,
2921 817 s->numXtiles * s->numYtiles);
2922
2923 817 jpeg2000_dec_cleanup(s);
2924
2925
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 817 times.
 817 if (s->avctx->pix_fmt == AV_PIX_FMT_PAL8)
2926 ✗ memcpy(picture->data[1], s->palette, 256 * sizeof(uint32_t));
2927
2928 817 return 0;
2929
2930 ✗ end:
2931 /* Do not leave a partially initialized frame on the output. */
2932 ✗ av_frame_unref(picture);
2933 ✗ jpeg2000_dec_cleanup(s);
2934 ✗ return ret;
2935 }
2936
2937 1734 static int jpeg2000_receive_frame(AVCodecContext *avctx, AVFrame *frame)
2938 {
2939 1734 AVPacket *const avpkt = avctx->internal->in_pkt;
2940 int ret;
2941
2942 /* One codestream per packet. Loop so that skipped frames (produce no
2943 * output) transparently pull the next packet instead of returning to the
2944 * caller empty-handed. */
2945 do {
2946 1768 ret = ff_decode_get_packet(avctx, avpkt);
2947
2/2
✓ Branch 0 taken 917 times.
✓ Branch 1 taken 851 times.
 1768 if (ret < 0)
2948 917 return ret;
2949
2950 851 ret = jpeg2000_decode_frame(avctx, frame, avpkt);
2951 851 av_packet_unref(avpkt);
2952
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 851 times.
 851 if (ret < 0)
2953 ✗ return ret;
2954
2/2
✓ Branch 0 taken 34 times.
✓ Branch 1 taken 817 times.
 851 } while (!frame->buf[0]);
2955
2956 817 return 0;
2957 }
2958
2959 #define OFFSET(x) offsetof(Jpeg2000DecoderContext, x)
2960 #define VD AV_OPT_FLAG_VIDEO_PARAM | AV_OPT_FLAG_DECODING_PARAM
2961
2962 static const AVOption options[] = {
2963 { "lowres", "Lower the decoding resolution by a power of two",
2964 OFFSET(reduction_factor), AV_OPT_TYPE_INT, { .i64 = 0 }, 0, JPEG2000_MAX_RESLEVELS - 1, VD },
2965 { NULL },
2966 };
2967
2968 static const AVClass jpeg2000_class = {
2969 .class_name = "jpeg2000",
2970 .item_name = av_default_item_name,
2971 .option = options,
2972 .version = LIBAVUTIL_VERSION_INT,
2973 };
2974
2975 const FFCodec ff_jpeg2000_decoder = {
2976 .p.name = "jpeg2000",
2977 CODEC_LONG_NAME("JPEG 2000"),
2978 .p.type = AVMEDIA_TYPE_VIDEO,
2979 .p.id = AV_CODEC_ID_JPEG2000,
2980 .p.capabilities = AV_CODEC_CAP_SLICE_THREADS | AV_CODEC_CAP_FRAME_THREADS | AV_CODEC_CAP_DR1,
2981 .priv_data_size = sizeof(Jpeg2000DecoderContext),
2982 .init = jpeg2000_decode_init,
2983 FF_CODEC_RECEIVE_FRAME_CB(jpeg2000_receive_frame),
2984 .p.priv_class = &jpeg2000_class,
2985 .p.max_lowres = 5,
2986 .p.profiles = NULL_IF_CONFIG_SMALL(ff_jpeg2000_profiles),
2987 .caps_internal = FF_CODEC_CAP_SKIP_FRAME_FILL_PARAM,
2988 };
2989