FFmpeg coverage

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