You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

1697 lines
58KB

  1. /*
  2. * MJPEG decoder
  3. * Copyright (c) 2000, 2001 Fabrice Bellard
  4. * Copyright (c) 2003 Alex Beregszaszi
  5. * Copyright (c) 2003-2004 Michael Niedermayer
  6. *
  7. * Support for external huffman table, various fixes (AVID workaround),
  8. * aspecting, new decode_frame mechanism and apple mjpeg-b support
  9. * by Alex Beregszaszi
  10. *
  11. * This file is part of Libav.
  12. *
  13. * Libav is free software; you can redistribute it and/or
  14. * modify it under the terms of the GNU Lesser General Public
  15. * License as published by the Free Software Foundation; either
  16. * version 2.1 of the License, or (at your option) any later version.
  17. *
  18. * Libav is distributed in the hope that it will be useful,
  19. * but WITHOUT ANY WARRANTY; without even the implied warranty of
  20. * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
  21. * Lesser General Public License for more details.
  22. *
  23. * You should have received a copy of the GNU Lesser General Public
  24. * License along with Libav; if not, write to the Free Software
  25. * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  26. */
  27. /**
  28. * @file
  29. * MJPEG decoder.
  30. */
  31. // #define DEBUG
  32. #include <assert.h>
  33. #include "libavutil/imgutils.h"
  34. #include "libavutil/opt.h"
  35. #include "avcodec.h"
  36. #include "dsputil.h"
  37. #include "mjpeg.h"
  38. #include "mjpegdec.h"
  39. #include "jpeglsdec.h"
  40. static int build_vlc(VLC *vlc, const uint8_t *bits_table,
  41. const uint8_t *val_table, int nb_codes,
  42. int use_static, int is_ac)
  43. {
  44. uint8_t huff_size[256] = { 0 };
  45. uint16_t huff_code[256];
  46. uint16_t huff_sym[256];
  47. int i;
  48. assert(nb_codes <= 256);
  49. ff_mjpeg_build_huffman_codes(huff_size, huff_code, bits_table, val_table);
  50. for (i = 0; i < 256; i++)
  51. huff_sym[i] = i + 16 * is_ac;
  52. if (is_ac)
  53. huff_sym[0] = 16 * 256;
  54. return ff_init_vlc_sparse(vlc, 9, nb_codes, huff_size, 1, 1,
  55. huff_code, 2, 2, huff_sym, 2, 2, use_static);
  56. }
  57. static void build_basic_mjpeg_vlc(MJpegDecodeContext *s)
  58. {
  59. build_vlc(&s->vlcs[0][0], avpriv_mjpeg_bits_dc_luminance,
  60. avpriv_mjpeg_val_dc, 12, 0, 0);
  61. build_vlc(&s->vlcs[0][1], avpriv_mjpeg_bits_dc_chrominance,
  62. avpriv_mjpeg_val_dc, 12, 0, 0);
  63. build_vlc(&s->vlcs[1][0], avpriv_mjpeg_bits_ac_luminance,
  64. avpriv_mjpeg_val_ac_luminance, 251, 0, 1);
  65. build_vlc(&s->vlcs[1][1], avpriv_mjpeg_bits_ac_chrominance,
  66. avpriv_mjpeg_val_ac_chrominance, 251, 0, 1);
  67. build_vlc(&s->vlcs[2][0], avpriv_mjpeg_bits_ac_luminance,
  68. avpriv_mjpeg_val_ac_luminance, 251, 0, 0);
  69. build_vlc(&s->vlcs[2][1], avpriv_mjpeg_bits_ac_chrominance,
  70. avpriv_mjpeg_val_ac_chrominance, 251, 0, 0);
  71. }
  72. av_cold int ff_mjpeg_decode_init(AVCodecContext *avctx)
  73. {
  74. MJpegDecodeContext *s = avctx->priv_data;
  75. if (!s->picture_ptr)
  76. s->picture_ptr = &s->picture;
  77. s->avctx = avctx;
  78. ff_dsputil_init(&s->dsp, avctx);
  79. ff_init_scantable(s->dsp.idct_permutation, &s->scantable, ff_zigzag_direct);
  80. s->buffer_size = 0;
  81. s->buffer = NULL;
  82. s->start_code = -1;
  83. s->first_picture = 1;
  84. s->org_height = avctx->coded_height;
  85. avctx->chroma_sample_location = AVCHROMA_LOC_CENTER;
  86. build_basic_mjpeg_vlc(s);
  87. if (s->extern_huff) {
  88. int ret;
  89. av_log(avctx, AV_LOG_INFO, "mjpeg: using external huffman table\n");
  90. init_get_bits(&s->gb, avctx->extradata, avctx->extradata_size * 8);
  91. if ((ret = ff_mjpeg_decode_dht(s))) {
  92. av_log(avctx, AV_LOG_ERROR,
  93. "mjpeg: error using external huffman table\n");
  94. return ret;
  95. }
  96. }
  97. if (avctx->field_order == AV_FIELD_BB) { /* quicktime icefloe 019 */
  98. s->interlace_polarity = 1; /* bottom field first */
  99. av_log(avctx, AV_LOG_DEBUG, "mjpeg bottom field first\n");
  100. }
  101. if (avctx->codec->id == AV_CODEC_ID_AMV)
  102. s->flipped = 1;
  103. return 0;
  104. }
  105. /* quantize tables */
  106. int ff_mjpeg_decode_dqt(MJpegDecodeContext *s)
  107. {
  108. int len, index, i, j;
  109. len = get_bits(&s->gb, 16) - 2;
  110. while (len >= 65) {
  111. /* only 8 bit precision handled */
  112. if (get_bits(&s->gb, 4) != 0) {
  113. av_log(s->avctx, AV_LOG_ERROR, "dqt: 16bit precision\n");
  114. return -1;
  115. }
  116. index = get_bits(&s->gb, 4);
  117. if (index >= 4)
  118. return -1;
  119. av_log(s->avctx, AV_LOG_DEBUG, "index=%d\n", index);
  120. /* read quant table */
  121. for (i = 0; i < 64; i++) {
  122. j = s->scantable.permutated[i];
  123. s->quant_matrixes[index][j] = get_bits(&s->gb, 8);
  124. }
  125. // XXX FIXME finetune, and perhaps add dc too
  126. s->qscale[index] = FFMAX(s->quant_matrixes[index][s->scantable.permutated[1]],
  127. s->quant_matrixes[index][s->scantable.permutated[8]]) >> 1;
  128. av_log(s->avctx, AV_LOG_DEBUG, "qscale[%d]: %d\n",
  129. index, s->qscale[index]);
  130. len -= 65;
  131. }
  132. return 0;
  133. }
  134. /* decode huffman tables and build VLC decoders */
  135. int ff_mjpeg_decode_dht(MJpegDecodeContext *s)
  136. {
  137. int len, index, i, class, n, v, code_max;
  138. uint8_t bits_table[17];
  139. uint8_t val_table[256];
  140. int ret = 0;
  141. len = get_bits(&s->gb, 16) - 2;
  142. while (len > 0) {
  143. if (len < 17)
  144. return AVERROR_INVALIDDATA;
  145. class = get_bits(&s->gb, 4);
  146. if (class >= 2)
  147. return AVERROR_INVALIDDATA;
  148. index = get_bits(&s->gb, 4);
  149. if (index >= 4)
  150. return AVERROR_INVALIDDATA;
  151. n = 0;
  152. for (i = 1; i <= 16; i++) {
  153. bits_table[i] = get_bits(&s->gb, 8);
  154. n += bits_table[i];
  155. }
  156. len -= 17;
  157. if (len < n || n > 256)
  158. return AVERROR_INVALIDDATA;
  159. code_max = 0;
  160. for (i = 0; i < n; i++) {
  161. v = get_bits(&s->gb, 8);
  162. if (v > code_max)
  163. code_max = v;
  164. val_table[i] = v;
  165. }
  166. len -= n;
  167. /* build VLC and flush previous vlc if present */
  168. ff_free_vlc(&s->vlcs[class][index]);
  169. av_log(s->avctx, AV_LOG_DEBUG, "class=%d index=%d nb_codes=%d\n",
  170. class, index, code_max + 1);
  171. if ((ret = build_vlc(&s->vlcs[class][index], bits_table, val_table,
  172. code_max + 1, 0, class > 0)) < 0)
  173. return ret;
  174. if (class > 0) {
  175. ff_free_vlc(&s->vlcs[2][index]);
  176. if ((ret = build_vlc(&s->vlcs[2][index], bits_table, val_table,
  177. code_max + 1, 0, 0)) < 0)
  178. return ret;
  179. }
  180. }
  181. return 0;
  182. }
  183. int ff_mjpeg_decode_sof(MJpegDecodeContext *s)
  184. {
  185. int len, nb_components, i, width, height, pix_fmt_id;
  186. /* XXX: verify len field validity */
  187. len = get_bits(&s->gb, 16);
  188. s->bits = get_bits(&s->gb, 8);
  189. if (s->pegasus_rct)
  190. s->bits = 9;
  191. if (s->bits == 9 && !s->pegasus_rct)
  192. s->rct = 1; // FIXME ugly
  193. if (s->bits != 8 && !s->lossless) {
  194. av_log(s->avctx, AV_LOG_ERROR, "only 8 bits/component accepted\n");
  195. return -1;
  196. }
  197. height = get_bits(&s->gb, 16);
  198. width = get_bits(&s->gb, 16);
  199. // HACK for odd_height.mov
  200. if (s->interlaced && s->width == width && s->height == height + 1)
  201. height= s->height;
  202. av_log(s->avctx, AV_LOG_DEBUG, "sof0: picture: %dx%d\n", width, height);
  203. if (av_image_check_size(width, height, 0, s->avctx))
  204. return AVERROR_INVALIDDATA;
  205. nb_components = get_bits(&s->gb, 8);
  206. if (nb_components <= 0 ||
  207. nb_components > MAX_COMPONENTS)
  208. return -1;
  209. if (s->ls && !(s->bits <= 8 || nb_components == 1)) {
  210. av_log_missing_feature(s->avctx,
  211. "only <= 8 bits/component or "
  212. "16-bit gray accepted for JPEG-LS\n", 0);
  213. return AVERROR_PATCHWELCOME;
  214. }
  215. s->nb_components = nb_components;
  216. s->h_max = 1;
  217. s->v_max = 1;
  218. for (i = 0; i < nb_components; i++) {
  219. /* component id */
  220. s->component_id[i] = get_bits(&s->gb, 8) - 1;
  221. s->h_count[i] = get_bits(&s->gb, 4);
  222. s->v_count[i] = get_bits(&s->gb, 4);
  223. /* compute hmax and vmax (only used in interleaved case) */
  224. if (s->h_count[i] > s->h_max)
  225. s->h_max = s->h_count[i];
  226. if (s->v_count[i] > s->v_max)
  227. s->v_max = s->v_count[i];
  228. s->quant_index[i] = get_bits(&s->gb, 8);
  229. if (s->quant_index[i] >= 4)
  230. return AVERROR_INVALIDDATA;
  231. av_log(s->avctx, AV_LOG_DEBUG, "component %d %d:%d id: %d quant:%d\n",
  232. i, s->h_count[i], s->v_count[i],
  233. s->component_id[i], s->quant_index[i]);
  234. }
  235. if (s->ls && (s->h_max > 1 || s->v_max > 1)) {
  236. av_log_missing_feature(s->avctx,
  237. "Subsampling in JPEG-LS is not supported.\n", 0);
  238. return AVERROR_PATCHWELCOME;
  239. }
  240. if (s->v_max == 1 && s->h_max == 1 && s->lossless == 1)
  241. s->rgb = 1;
  242. /* if different size, realloc/alloc picture */
  243. /* XXX: also check h_count and v_count */
  244. if (width != s->width || height != s->height) {
  245. av_freep(&s->qscale_table);
  246. s->width = width;
  247. s->height = height;
  248. s->interlaced = 0;
  249. /* test interlaced mode */
  250. if (s->first_picture &&
  251. s->org_height != 0 &&
  252. s->height < ((s->org_height * 3) / 4)) {
  253. s->interlaced = 1;
  254. s->bottom_field = s->interlace_polarity;
  255. s->picture_ptr->interlaced_frame = 1;
  256. s->picture_ptr->top_field_first = !s->interlace_polarity;
  257. height *= 2;
  258. }
  259. avcodec_set_dimensions(s->avctx, width, height);
  260. s->qscale_table = av_mallocz((s->width + 15) / 16);
  261. s->first_picture = 0;
  262. }
  263. if (!(s->interlaced && (s->bottom_field == !s->interlace_polarity))) {
  264. /* XXX: not complete test ! */
  265. pix_fmt_id = (s->h_count[0] << 28) | (s->v_count[0] << 24) |
  266. (s->h_count[1] << 20) | (s->v_count[1] << 16) |
  267. (s->h_count[2] << 12) | (s->v_count[2] << 8) |
  268. (s->h_count[3] << 4) | s->v_count[3];
  269. av_log(s->avctx, AV_LOG_DEBUG, "pix fmt id %x\n", pix_fmt_id);
  270. /* NOTE we do not allocate pictures large enough for the possible
  271. * padding of h/v_count being 4 */
  272. if (!(pix_fmt_id & 0xD0D0D0D0))
  273. pix_fmt_id -= (pix_fmt_id & 0xF0F0F0F0) >> 1;
  274. if (!(pix_fmt_id & 0x0D0D0D0D))
  275. pix_fmt_id -= (pix_fmt_id & 0x0F0F0F0F) >> 1;
  276. switch (pix_fmt_id) {
  277. case 0x11111100:
  278. if (s->rgb)
  279. s->avctx->pix_fmt = PIX_FMT_BGRA;
  280. else
  281. s->avctx->pix_fmt = s->cs_itu601 ? PIX_FMT_YUV444P : PIX_FMT_YUVJ444P;
  282. assert(s->nb_components == 3);
  283. break;
  284. case 0x11000000:
  285. s->avctx->pix_fmt = PIX_FMT_GRAY8;
  286. break;
  287. case 0x12111100:
  288. s->avctx->pix_fmt = s->cs_itu601 ? PIX_FMT_YUV440P : PIX_FMT_YUVJ440P;
  289. break;
  290. case 0x21111100:
  291. s->avctx->pix_fmt = s->cs_itu601 ? PIX_FMT_YUV422P : PIX_FMT_YUVJ422P;
  292. break;
  293. case 0x22111100:
  294. s->avctx->pix_fmt = s->cs_itu601 ? PIX_FMT_YUV420P : PIX_FMT_YUVJ420P;
  295. break;
  296. default:
  297. av_log(s->avctx, AV_LOG_ERROR, "Unhandled pixel format 0x%x\n", pix_fmt_id);
  298. return AVERROR_PATCHWELCOME;
  299. }
  300. if (s->ls) {
  301. if (s->nb_components > 1)
  302. s->avctx->pix_fmt = PIX_FMT_RGB24;
  303. else if (s->bits <= 8)
  304. s->avctx->pix_fmt = PIX_FMT_GRAY8;
  305. else
  306. s->avctx->pix_fmt = PIX_FMT_GRAY16;
  307. }
  308. if (s->picture_ptr->data[0])
  309. s->avctx->release_buffer(s->avctx, s->picture_ptr);
  310. if (s->avctx->get_buffer(s->avctx, s->picture_ptr) < 0) {
  311. av_log(s->avctx, AV_LOG_ERROR, "get_buffer() failed\n");
  312. return -1;
  313. }
  314. s->picture_ptr->pict_type = AV_PICTURE_TYPE_I;
  315. s->picture_ptr->key_frame = 1;
  316. s->got_picture = 1;
  317. for (i = 0; i < 3; i++)
  318. s->linesize[i] = s->picture_ptr->linesize[i] << s->interlaced;
  319. // printf("%d %d %d %d %d %d\n",
  320. // s->width, s->height, s->linesize[0], s->linesize[1],
  321. // s->interlaced, s->avctx->height);
  322. if (len != (8 + (3 * nb_components)))
  323. av_log(s->avctx, AV_LOG_DEBUG, "decode_sof0: error, len(%d) mismatch\n", len);
  324. }
  325. /* totally blank picture as progressive JPEG will only add details to it */
  326. if (s->progressive) {
  327. int bw = (width + s->h_max * 8 - 1) / (s->h_max * 8);
  328. int bh = (height + s->v_max * 8 - 1) / (s->v_max * 8);
  329. for (i = 0; i < s->nb_components; i++) {
  330. int size = bw * bh * s->h_count[i] * s->v_count[i];
  331. av_freep(&s->blocks[i]);
  332. av_freep(&s->last_nnz[i]);
  333. s->blocks[i] = av_malloc(size * sizeof(**s->blocks));
  334. s->last_nnz[i] = av_mallocz(size * sizeof(**s->last_nnz));
  335. s->block_stride[i] = bw * s->h_count[i];
  336. }
  337. memset(s->coefs_finished, 0, sizeof(s->coefs_finished));
  338. }
  339. return 0;
  340. }
  341. static inline int mjpeg_decode_dc(MJpegDecodeContext *s, int dc_index)
  342. {
  343. int code;
  344. code = get_vlc2(&s->gb, s->vlcs[0][dc_index].table, 9, 2);
  345. if (code < 0) {
  346. av_log(s->avctx, AV_LOG_WARNING,
  347. "mjpeg_decode_dc: bad vlc: %d:%d (%p)\n",
  348. 0, dc_index, &s->vlcs[0][dc_index]);
  349. return 0xffff;
  350. }
  351. if (code)
  352. return get_xbits(&s->gb, code);
  353. else
  354. return 0;
  355. }
  356. /* decode block and dequantize */
  357. static int decode_block(MJpegDecodeContext *s, DCTELEM *block, int component,
  358. int dc_index, int ac_index, int16_t *quant_matrix)
  359. {
  360. int code, i, j, level, val;
  361. /* DC coef */
  362. val = mjpeg_decode_dc(s, dc_index);
  363. if (val == 0xffff) {
  364. av_log(s->avctx, AV_LOG_ERROR, "error dc\n");
  365. return AVERROR_INVALIDDATA;
  366. }
  367. val = val * quant_matrix[0] + s->last_dc[component];
  368. s->last_dc[component] = val;
  369. block[0] = val;
  370. /* AC coefs */
  371. i = 0;
  372. {OPEN_READER(re, &s->gb);
  373. do {
  374. UPDATE_CACHE(re, &s->gb);
  375. GET_VLC(code, re, &s->gb, s->vlcs[1][ac_index].table, 9, 2);
  376. i += ((unsigned)code) >> 4;
  377. code &= 0xf;
  378. if (code) {
  379. if (code > MIN_CACHE_BITS - 16)
  380. UPDATE_CACHE(re, &s->gb);
  381. {
  382. int cache = GET_CACHE(re, &s->gb);
  383. int sign = (~cache) >> 31;
  384. level = (NEG_USR32(sign ^ cache,code) ^ sign) - sign;
  385. }
  386. LAST_SKIP_BITS(re, &s->gb, code);
  387. if (i > 63) {
  388. av_log(s->avctx, AV_LOG_ERROR, "error count: %d\n", i);
  389. return AVERROR_INVALIDDATA;
  390. }
  391. j = s->scantable.permutated[i];
  392. block[j] = level * quant_matrix[j];
  393. }
  394. } while (i < 63);
  395. CLOSE_READER(re, &s->gb);}
  396. return 0;
  397. }
  398. static int decode_dc_progressive(MJpegDecodeContext *s, DCTELEM *block,
  399. int component, int dc_index,
  400. int16_t *quant_matrix, int Al)
  401. {
  402. int val;
  403. s->dsp.clear_block(block);
  404. val = mjpeg_decode_dc(s, dc_index);
  405. if (val == 0xffff) {
  406. av_log(s->avctx, AV_LOG_ERROR, "error dc\n");
  407. return AVERROR_INVALIDDATA;
  408. }
  409. val = (val * quant_matrix[0] << Al) + s->last_dc[component];
  410. s->last_dc[component] = val;
  411. block[0] = val;
  412. return 0;
  413. }
  414. /* decode block and dequantize - progressive JPEG version */
  415. static int decode_block_progressive(MJpegDecodeContext *s, DCTELEM *block,
  416. uint8_t *last_nnz, int ac_index,
  417. int16_t *quant_matrix,
  418. int ss, int se, int Al, int *EOBRUN)
  419. {
  420. int code, i, j, level, val, run;
  421. if (*EOBRUN) {
  422. (*EOBRUN)--;
  423. return 0;
  424. }
  425. {
  426. OPEN_READER(re, &s->gb);
  427. for (i = ss; ; i++) {
  428. UPDATE_CACHE(re, &s->gb);
  429. GET_VLC(code, re, &s->gb, s->vlcs[2][ac_index].table, 9, 2);
  430. run = ((unsigned) code) >> 4;
  431. code &= 0xF;
  432. if (code) {
  433. i += run;
  434. if (code > MIN_CACHE_BITS - 16)
  435. UPDATE_CACHE(re, &s->gb);
  436. {
  437. int cache = GET_CACHE(re, &s->gb);
  438. int sign = (~cache) >> 31;
  439. level = (NEG_USR32(sign ^ cache,code) ^ sign) - sign;
  440. }
  441. LAST_SKIP_BITS(re, &s->gb, code);
  442. if (i >= se) {
  443. if (i == se) {
  444. j = s->scantable.permutated[se];
  445. block[j] = level * quant_matrix[j] << Al;
  446. break;
  447. }
  448. av_log(s->avctx, AV_LOG_ERROR, "error count: %d\n", i);
  449. return AVERROR_INVALIDDATA;
  450. }
  451. j = s->scantable.permutated[i];
  452. block[j] = level * quant_matrix[j] << Al;
  453. } else {
  454. if (run == 0xF) {// ZRL - skip 15 coefficients
  455. i += 15;
  456. if (i >= se) {
  457. av_log(s->avctx, AV_LOG_ERROR, "ZRL overflow: %d\n", i);
  458. return AVERROR_INVALIDDATA;
  459. }
  460. } else {
  461. val = (1 << run);
  462. if (run) {
  463. UPDATE_CACHE(re, &s->gb);
  464. val += NEG_USR32(GET_CACHE(re, &s->gb), run);
  465. LAST_SKIP_BITS(re, &s->gb, run);
  466. }
  467. *EOBRUN = val - 1;
  468. break;
  469. }
  470. }
  471. }
  472. CLOSE_READER(re, &s->gb);
  473. }
  474. if (i > *last_nnz)
  475. *last_nnz = i;
  476. return 0;
  477. }
  478. #define REFINE_BIT(j) { \
  479. UPDATE_CACHE(re, &s->gb); \
  480. sign = block[j] >> 15; \
  481. block[j] += SHOW_UBITS(re, &s->gb, 1) * \
  482. ((quant_matrix[j] ^ sign) - sign) << Al; \
  483. LAST_SKIP_BITS(re, &s->gb, 1); \
  484. }
  485. #define ZERO_RUN \
  486. for (; ; i++) { \
  487. if (i > last) { \
  488. i += run; \
  489. if (i > se) { \
  490. av_log(s->avctx, AV_LOG_ERROR, "error count: %d\n", i); \
  491. return -1; \
  492. } \
  493. break; \
  494. } \
  495. j = s->scantable.permutated[i]; \
  496. if (block[j]) \
  497. REFINE_BIT(j) \
  498. else if (run-- == 0) \
  499. break; \
  500. }
  501. /* decode block and dequantize - progressive JPEG refinement pass */
  502. static int decode_block_refinement(MJpegDecodeContext *s, DCTELEM *block,
  503. uint8_t *last_nnz,
  504. int ac_index, int16_t *quant_matrix,
  505. int ss, int se, int Al, int *EOBRUN)
  506. {
  507. int code, i = ss, j, sign, val, run;
  508. int last = FFMIN(se, *last_nnz);
  509. OPEN_READER(re, &s->gb);
  510. if (*EOBRUN) {
  511. (*EOBRUN)--;
  512. } else {
  513. for (; ; i++) {
  514. UPDATE_CACHE(re, &s->gb);
  515. GET_VLC(code, re, &s->gb, s->vlcs[2][ac_index].table, 9, 2);
  516. if (code & 0xF) {
  517. run = ((unsigned) code) >> 4;
  518. UPDATE_CACHE(re, &s->gb);
  519. val = SHOW_UBITS(re, &s->gb, 1);
  520. LAST_SKIP_BITS(re, &s->gb, 1);
  521. ZERO_RUN;
  522. j = s->scantable.permutated[i];
  523. val--;
  524. block[j] = ((quant_matrix[j]^val) - val) << Al;
  525. if (i == se) {
  526. if (i > *last_nnz)
  527. *last_nnz = i;
  528. CLOSE_READER(re, &s->gb);
  529. return 0;
  530. }
  531. } else {
  532. run = ((unsigned) code) >> 4;
  533. if (run == 0xF) {
  534. ZERO_RUN;
  535. } else {
  536. val = run;
  537. run = (1 << run);
  538. if (val) {
  539. UPDATE_CACHE(re, &s->gb);
  540. run += SHOW_UBITS(re, &s->gb, val);
  541. LAST_SKIP_BITS(re, &s->gb, val);
  542. }
  543. *EOBRUN = run - 1;
  544. break;
  545. }
  546. }
  547. }
  548. if (i > *last_nnz)
  549. *last_nnz = i;
  550. }
  551. for (; i <= last; i++) {
  552. j = s->scantable.permutated[i];
  553. if (block[j])
  554. REFINE_BIT(j)
  555. }
  556. CLOSE_READER(re, &s->gb);
  557. return 0;
  558. }
  559. #undef REFINE_BIT
  560. #undef ZERO_RUN
  561. static int ljpeg_decode_rgb_scan(MJpegDecodeContext *s, int predictor,
  562. int point_transform)
  563. {
  564. int i, mb_x, mb_y;
  565. uint16_t (*buffer)[4];
  566. int left[3], top[3], topleft[3];
  567. const int linesize = s->linesize[0];
  568. const int mask = (1 << s->bits) - 1;
  569. av_fast_malloc(&s->ljpeg_buffer, &s->ljpeg_buffer_size,
  570. (unsigned)s->mb_width * 4 * sizeof(s->ljpeg_buffer[0][0]));
  571. buffer = s->ljpeg_buffer;
  572. for (i = 0; i < 3; i++)
  573. buffer[0][i] = 1 << (s->bits + point_transform - 1);
  574. for (mb_y = 0; mb_y < s->mb_height; mb_y++) {
  575. const int modified_predictor = mb_y ? predictor : 1;
  576. uint8_t *ptr = s->picture_ptr->data[0] + (linesize * mb_y);
  577. if (s->interlaced && s->bottom_field)
  578. ptr += linesize >> 1;
  579. for (i = 0; i < 3; i++)
  580. top[i] = left[i] = topleft[i] = buffer[0][i];
  581. for (mb_x = 0; mb_x < s->mb_width; mb_x++) {
  582. if (s->restart_interval && !s->restart_count)
  583. s->restart_count = s->restart_interval;
  584. for (i = 0; i < 3; i++) {
  585. int pred;
  586. topleft[i] = top[i];
  587. top[i] = buffer[mb_x][i];
  588. PREDICT(pred, topleft[i], top[i], left[i], modified_predictor);
  589. left[i] = buffer[mb_x][i] =
  590. mask & (pred + (mjpeg_decode_dc(s, s->dc_index[i]) << point_transform));
  591. }
  592. if (s->restart_interval && !--s->restart_count) {
  593. align_get_bits(&s->gb);
  594. skip_bits(&s->gb, 16); /* skip RSTn */
  595. }
  596. }
  597. if (s->rct) {
  598. for (mb_x = 0; mb_x < s->mb_width; mb_x++) {
  599. ptr[4 * mb_x + 1] = buffer[mb_x][0] - ((buffer[mb_x][1] + buffer[mb_x][2] - 0x200) >> 2);
  600. ptr[4 * mb_x + 0] = buffer[mb_x][1] + ptr[4 * mb_x + 1];
  601. ptr[4 * mb_x + 2] = buffer[mb_x][2] + ptr[4 * mb_x + 1];
  602. }
  603. } else if (s->pegasus_rct) {
  604. for (mb_x = 0; mb_x < s->mb_width; mb_x++) {
  605. ptr[4 * mb_x + 1] = buffer[mb_x][0] - ((buffer[mb_x][1] + buffer[mb_x][2]) >> 2);
  606. ptr[4 * mb_x + 0] = buffer[mb_x][1] + ptr[4 * mb_x + 1];
  607. ptr[4 * mb_x + 2] = buffer[mb_x][2] + ptr[4 * mb_x + 1];
  608. }
  609. } else {
  610. for (mb_x = 0; mb_x < s->mb_width; mb_x++) {
  611. ptr[4 * mb_x + 0] = buffer[mb_x][2];
  612. ptr[4 * mb_x + 1] = buffer[mb_x][1];
  613. ptr[4 * mb_x + 2] = buffer[mb_x][0];
  614. }
  615. }
  616. }
  617. return 0;
  618. }
  619. static int ljpeg_decode_yuv_scan(MJpegDecodeContext *s, int predictor,
  620. int point_transform)
  621. {
  622. int i, mb_x, mb_y;
  623. const int nb_components = 3;
  624. for (mb_y = 0; mb_y < s->mb_height; mb_y++) {
  625. for (mb_x = 0; mb_x < s->mb_width; mb_x++) {
  626. if (s->restart_interval && !s->restart_count)
  627. s->restart_count = s->restart_interval;
  628. if (mb_x == 0 || mb_y == 0 || s->interlaced) {
  629. for (i = 0; i < nb_components; i++) {
  630. uint8_t *ptr;
  631. int n, h, v, x, y, c, j, linesize;
  632. n = s->nb_blocks[i];
  633. c = s->comp_index[i];
  634. h = s->h_scount[i];
  635. v = s->v_scount[i];
  636. x = 0;
  637. y = 0;
  638. linesize = s->linesize[c];
  639. for (j = 0; j < n; j++) {
  640. int pred;
  641. // FIXME optimize this crap
  642. ptr = s->picture_ptr->data[c] +
  643. (linesize * (v * mb_y + y)) +
  644. (h * mb_x + x);
  645. if (y == 0 && mb_y == 0) {
  646. if (x == 0 && mb_x == 0)
  647. pred = 128 << point_transform;
  648. else
  649. pred = ptr[-1];
  650. } else {
  651. if (x == 0 && mb_x == 0)
  652. pred = ptr[-linesize];
  653. else
  654. PREDICT(pred, ptr[-linesize - 1],
  655. ptr[-linesize], ptr[-1], predictor);
  656. }
  657. if (s->interlaced && s->bottom_field)
  658. ptr += linesize >> 1;
  659. *ptr = pred + (mjpeg_decode_dc(s, s->dc_index[i]) << point_transform);
  660. if (++x == h) {
  661. x = 0;
  662. y++;
  663. }
  664. }
  665. }
  666. } else {
  667. for (i = 0; i < nb_components; i++) {
  668. uint8_t *ptr;
  669. int n, h, v, x, y, c, j, linesize;
  670. n = s->nb_blocks[i];
  671. c = s->comp_index[i];
  672. h = s->h_scount[i];
  673. v = s->v_scount[i];
  674. x = 0;
  675. y = 0;
  676. linesize = s->linesize[c];
  677. for (j = 0; j < n; j++) {
  678. int pred;
  679. // FIXME optimize this crap
  680. ptr = s->picture_ptr->data[c] +
  681. (linesize * (v * mb_y + y)) +
  682. (h * mb_x + x);
  683. PREDICT(pred, ptr[-linesize - 1],
  684. ptr[-linesize], ptr[-1], predictor);
  685. *ptr = pred + (mjpeg_decode_dc(s, s->dc_index[i]) << point_transform);
  686. if (++x == h) {
  687. x = 0;
  688. y++;
  689. }
  690. }
  691. }
  692. }
  693. if (s->restart_interval && !--s->restart_count) {
  694. align_get_bits(&s->gb);
  695. skip_bits(&s->gb, 16); /* skip RSTn */
  696. }
  697. }
  698. }
  699. return 0;
  700. }
  701. static int mjpeg_decode_scan(MJpegDecodeContext *s, int nb_components, int Ah,
  702. int Al, const uint8_t *mb_bitmask,
  703. const AVFrame *reference)
  704. {
  705. int i, mb_x, mb_y;
  706. uint8_t *data[MAX_COMPONENTS];
  707. const uint8_t *reference_data[MAX_COMPONENTS];
  708. int linesize[MAX_COMPONENTS];
  709. GetBitContext mb_bitmask_gb;
  710. if (mb_bitmask)
  711. init_get_bits(&mb_bitmask_gb, mb_bitmask, s->mb_width * s->mb_height);
  712. if (s->flipped && s->avctx->flags & CODEC_FLAG_EMU_EDGE) {
  713. av_log(s->avctx, AV_LOG_ERROR,
  714. "Can not flip image with CODEC_FLAG_EMU_EDGE set!\n");
  715. s->flipped = 0;
  716. }
  717. for (i = 0; i < nb_components; i++) {
  718. int c = s->comp_index[i];
  719. data[c] = s->picture_ptr->data[c];
  720. reference_data[c] = reference ? reference->data[c] : NULL;
  721. linesize[c] = s->linesize[c];
  722. s->coefs_finished[c] |= 1;
  723. if (s->flipped) {
  724. // picture should be flipped upside-down for this codec
  725. int offset = (linesize[c] * (s->v_scount[i] *
  726. (8 * s->mb_height - ((s->height / s->v_max) & 7)) - 1));
  727. data[c] += offset;
  728. reference_data[c] += offset;
  729. linesize[c] *= -1;
  730. }
  731. }
  732. for (mb_y = 0; mb_y < s->mb_height; mb_y++) {
  733. for (mb_x = 0; mb_x < s->mb_width; mb_x++) {
  734. const int copy_mb = mb_bitmask && !get_bits1(&mb_bitmask_gb);
  735. if (s->restart_interval && !s->restart_count)
  736. s->restart_count = s->restart_interval;
  737. if (get_bits_left(&s->gb) < 0) {
  738. av_log(s->avctx, AV_LOG_ERROR, "overread %d\n",
  739. -get_bits_left(&s->gb));
  740. return AVERROR_INVALIDDATA;
  741. }
  742. for (i = 0; i < nb_components; i++) {
  743. uint8_t *ptr;
  744. int n, h, v, x, y, c, j;
  745. int block_offset;
  746. n = s->nb_blocks[i];
  747. c = s->comp_index[i];
  748. h = s->h_scount[i];
  749. v = s->v_scount[i];
  750. x = 0;
  751. y = 0;
  752. for (j = 0; j < n; j++) {
  753. block_offset = ((linesize[c] * (v * mb_y + y) * 8) +
  754. (h * mb_x + x) * 8);
  755. if (s->interlaced && s->bottom_field)
  756. block_offset += linesize[c] >> 1;
  757. ptr = data[c] + block_offset;
  758. if (!s->progressive) {
  759. if (copy_mb)
  760. copy_block8(ptr, reference_data[c] + block_offset,
  761. linesize[c], linesize[c], 8);
  762. else {
  763. s->dsp.clear_block(s->block);
  764. if (decode_block(s, s->block, i,
  765. s->dc_index[i], s->ac_index[i],
  766. s->quant_matrixes[s->quant_index[c]]) < 0) {
  767. av_log(s->avctx, AV_LOG_ERROR,
  768. "error y=%d x=%d\n", mb_y, mb_x);
  769. return AVERROR_INVALIDDATA;
  770. }
  771. s->dsp.idct_put(ptr, linesize[c], s->block);
  772. }
  773. } else {
  774. int block_idx = s->block_stride[c] * (v * mb_y + y) +
  775. (h * mb_x + x);
  776. DCTELEM *block = s->blocks[c][block_idx];
  777. if (Ah)
  778. block[0] += get_bits1(&s->gb) *
  779. s->quant_matrixes[s->quant_index[c]][0] << Al;
  780. else if (decode_dc_progressive(s, block, i, s->dc_index[i],
  781. s->quant_matrixes[s->quant_index[c]],
  782. Al) < 0) {
  783. av_log(s->avctx, AV_LOG_ERROR,
  784. "error y=%d x=%d\n", mb_y, mb_x);
  785. return AVERROR_INVALIDDATA;
  786. }
  787. }
  788. // av_log(s->avctx, AV_LOG_DEBUG, "mb: %d %d processed\n",
  789. // mb_y, mb_x);
  790. // av_log(NULL, AV_LOG_DEBUG, "%d %d %d %d %d %d %d %d \n",
  791. // mb_x, mb_y, x, y, c, s->bottom_field,
  792. // (v * mb_y + y) * 8, (h * mb_x + x) * 8);
  793. if (++x == h) {
  794. x = 0;
  795. y++;
  796. }
  797. }
  798. }
  799. if (s->restart_interval) {
  800. s->restart_count--;
  801. i = 8 + ((-get_bits_count(&s->gb)) & 7);
  802. /* skip RSTn */
  803. if (show_bits(&s->gb, i) == (1 << i) - 1) {
  804. int pos = get_bits_count(&s->gb);
  805. align_get_bits(&s->gb);
  806. while (get_bits_left(&s->gb) >= 8 && show_bits(&s->gb, 8) == 0xFF)
  807. skip_bits(&s->gb, 8);
  808. if ((get_bits(&s->gb, 8) & 0xF8) == 0xD0) {
  809. for (i = 0; i < nb_components; i++) /* reset dc */
  810. s->last_dc[i] = 1024;
  811. } else
  812. skip_bits_long(&s->gb, pos - get_bits_count(&s->gb));
  813. }
  814. }
  815. }
  816. }
  817. return 0;
  818. }
  819. static int mjpeg_decode_scan_progressive_ac(MJpegDecodeContext *s, int ss,
  820. int se, int Ah, int Al,
  821. const uint8_t *mb_bitmask,
  822. const AVFrame *reference)
  823. {
  824. int mb_x, mb_y;
  825. int EOBRUN = 0;
  826. int c = s->comp_index[0];
  827. uint8_t *data = s->picture_ptr->data[c];
  828. const uint8_t *reference_data = reference ? reference->data[c] : NULL;
  829. int linesize = s->linesize[c];
  830. int last_scan = 0;
  831. int16_t *quant_matrix = s->quant_matrixes[s->quant_index[c]];
  832. GetBitContext mb_bitmask_gb;
  833. if (mb_bitmask)
  834. init_get_bits(&mb_bitmask_gb, mb_bitmask, s->mb_width * s->mb_height);
  835. if (!Al) {
  836. s->coefs_finished[c] |= (1LL << (se + 1)) - (1LL << ss);
  837. last_scan = !~s->coefs_finished[c];
  838. }
  839. if (s->interlaced && s->bottom_field) {
  840. int offset = linesize >> 1;
  841. data += offset;
  842. reference_data += offset;
  843. }
  844. for (mb_y = 0; mb_y < s->mb_height; mb_y++) {
  845. int block_offset = mb_y * linesize * 8;
  846. uint8_t *ptr = data + block_offset;
  847. int block_idx = mb_y * s->block_stride[c];
  848. DCTELEM (*block)[64] = &s->blocks[c][block_idx];
  849. uint8_t *last_nnz = &s->last_nnz[c][block_idx];
  850. for (mb_x = 0; mb_x < s->mb_width; mb_x++, block++, last_nnz++) {
  851. const int copy_mb = mb_bitmask && !get_bits1(&mb_bitmask_gb);
  852. if (!copy_mb) {
  853. int ret;
  854. if (Ah)
  855. ret = decode_block_refinement(s, *block, last_nnz, s->ac_index[0],
  856. quant_matrix, ss, se, Al, &EOBRUN);
  857. else
  858. ret = decode_block_progressive(s, *block, last_nnz, s->ac_index[0],
  859. quant_matrix, ss, se, Al, &EOBRUN);
  860. if (ret < 0) {
  861. av_log(s->avctx, AV_LOG_ERROR,
  862. "error y=%d x=%d\n", mb_y, mb_x);
  863. return AVERROR_INVALIDDATA;
  864. }
  865. }
  866. if (last_scan) {
  867. if (copy_mb) {
  868. copy_block8(ptr, reference_data + block_offset,
  869. linesize, linesize, 8);
  870. } else {
  871. s->dsp.idct_put(ptr, linesize, *block);
  872. ptr += 8;
  873. }
  874. }
  875. }
  876. }
  877. return 0;
  878. }
  879. int ff_mjpeg_decode_sos(MJpegDecodeContext *s, const uint8_t *mb_bitmask,
  880. const AVFrame *reference)
  881. {
  882. int len, nb_components, i, h, v, predictor, point_transform;
  883. int index, id, ret;
  884. const int block_size = s->lossless ? 1 : 8;
  885. int ilv, prev_shift;
  886. /* XXX: verify len field validity */
  887. len = get_bits(&s->gb, 16);
  888. nb_components = get_bits(&s->gb, 8);
  889. if (nb_components == 0 || nb_components > MAX_COMPONENTS) {
  890. av_log(s->avctx, AV_LOG_ERROR,
  891. "decode_sos: nb_components (%d) unsupported\n", nb_components);
  892. return AVERROR_PATCHWELCOME;
  893. }
  894. if (len != 6 + 2 * nb_components) {
  895. av_log(s->avctx, AV_LOG_ERROR, "decode_sos: invalid len (%d)\n", len);
  896. return AVERROR_INVALIDDATA;
  897. }
  898. for (i = 0; i < nb_components; i++) {
  899. id = get_bits(&s->gb, 8) - 1;
  900. av_log(s->avctx, AV_LOG_DEBUG, "component: %d\n", id);
  901. /* find component index */
  902. for (index = 0; index < s->nb_components; index++)
  903. if (id == s->component_id[index])
  904. break;
  905. if (index == s->nb_components) {
  906. av_log(s->avctx, AV_LOG_ERROR,
  907. "decode_sos: index(%d) out of components\n", index);
  908. return AVERROR_INVALIDDATA;
  909. }
  910. /* Metasoft MJPEG codec has Cb and Cr swapped */
  911. if (s->avctx->codec_tag == MKTAG('M', 'T', 'S', 'J')
  912. && nb_components == 3 && s->nb_components == 3 && i)
  913. index = 3 - i;
  914. s->comp_index[i] = index;
  915. s->nb_blocks[i] = s->h_count[index] * s->v_count[index];
  916. s->h_scount[i] = s->h_count[index];
  917. s->v_scount[i] = s->v_count[index];
  918. s->dc_index[i] = get_bits(&s->gb, 4);
  919. s->ac_index[i] = get_bits(&s->gb, 4);
  920. if (s->dc_index[i] < 0 || s->ac_index[i] < 0 ||
  921. s->dc_index[i] >= 4 || s->ac_index[i] >= 4)
  922. goto out_of_range;
  923. if (!s->vlcs[0][s->dc_index[i]].table ||
  924. !s->vlcs[1][s->ac_index[i]].table)
  925. goto out_of_range;
  926. }
  927. predictor = get_bits(&s->gb, 8); /* JPEG Ss / lossless JPEG predictor /JPEG-LS NEAR */
  928. ilv = get_bits(&s->gb, 8); /* JPEG Se / JPEG-LS ILV */
  929. prev_shift = get_bits(&s->gb, 4); /* Ah */
  930. point_transform = get_bits(&s->gb, 4); /* Al */
  931. if (nb_components > 1) {
  932. /* interleaved stream */
  933. s->mb_width = (s->width + s->h_max * block_size - 1) / (s->h_max * block_size);
  934. s->mb_height = (s->height + s->v_max * block_size - 1) / (s->v_max * block_size);
  935. } else if (!s->ls) { /* skip this for JPEG-LS */
  936. h = s->h_max / s->h_scount[0];
  937. v = s->v_max / s->v_scount[0];
  938. s->mb_width = (s->width + h * block_size - 1) / (h * block_size);
  939. s->mb_height = (s->height + v * block_size - 1) / (v * block_size);
  940. s->nb_blocks[0] = 1;
  941. s->h_scount[0] = 1;
  942. s->v_scount[0] = 1;
  943. }
  944. if (s->avctx->debug & FF_DEBUG_PICT_INFO)
  945. av_log(s->avctx, AV_LOG_DEBUG, "%s %s p:%d >>:%d ilv:%d bits:%d %s\n",
  946. s->lossless ? "lossless" : "sequential DCT", s->rgb ? "RGB" : "",
  947. predictor, point_transform, ilv, s->bits,
  948. s->pegasus_rct ? "PRCT" : (s->rct ? "RCT" : ""));
  949. /* mjpeg-b can have padding bytes between sos and image data, skip them */
  950. for (i = s->mjpb_skiptosod; i > 0; i--)
  951. skip_bits(&s->gb, 8);
  952. next_field:
  953. for (i = 0; i < nb_components; i++)
  954. s->last_dc[i] = 1024;
  955. if (s->lossless) {
  956. if (CONFIG_JPEGLS_DECODER && s->ls) {
  957. // for () {
  958. // reset_ls_coding_parameters(s, 0);
  959. if ((ret = ff_jpegls_decode_picture(s, predictor,
  960. point_transform, ilv)) < 0)
  961. return ret;
  962. } else {
  963. if (s->rgb) {
  964. if ((ret = ljpeg_decode_rgb_scan(s, predictor,
  965. point_transform)) < 0)
  966. return ret;
  967. } else {
  968. if ((ret = ljpeg_decode_yuv_scan(s, predictor,
  969. point_transform)) < 0)
  970. return ret;
  971. }
  972. }
  973. } else {
  974. if (s->progressive && predictor) {
  975. if ((ret = mjpeg_decode_scan_progressive_ac(s, predictor,
  976. ilv, prev_shift,
  977. point_transform,
  978. mb_bitmask,
  979. reference)) < 0)
  980. return ret;
  981. } else {
  982. if ((ret = mjpeg_decode_scan(s, nb_components,
  983. prev_shift, point_transform,
  984. mb_bitmask, reference)) < 0)
  985. return ret;
  986. }
  987. }
  988. if (s->interlaced &&
  989. get_bits_left(&s->gb) > 32 &&
  990. show_bits(&s->gb, 8) == 0xFF) {
  991. GetBitContext bak = s->gb;
  992. align_get_bits(&bak);
  993. if (show_bits(&bak, 16) == 0xFFD1) {
  994. av_dlog(s->avctx, "AVRn interlaced picture marker found\n");
  995. s->gb = bak;
  996. skip_bits(&s->gb, 16);
  997. s->bottom_field ^= 1;
  998. goto next_field;
  999. }
  1000. }
  1001. emms_c();
  1002. return 0;
  1003. out_of_range:
  1004. av_log(s->avctx, AV_LOG_ERROR, "decode_sos: ac/dc index out of range\n");
  1005. return AVERROR_INVALIDDATA;
  1006. }
  1007. static int mjpeg_decode_dri(MJpegDecodeContext *s)
  1008. {
  1009. if (get_bits(&s->gb, 16) != 4)
  1010. return AVERROR_INVALIDDATA;
  1011. s->restart_interval = get_bits(&s->gb, 16);
  1012. s->restart_count = 0;
  1013. av_log(s->avctx, AV_LOG_DEBUG, "restart interval: %d\n",
  1014. s->restart_interval);
  1015. return 0;
  1016. }
  1017. static int mjpeg_decode_app(MJpegDecodeContext *s)
  1018. {
  1019. int len, id, i;
  1020. len = get_bits(&s->gb, 16);
  1021. if (len < 5)
  1022. return AVERROR_INVALIDDATA;
  1023. if (8 * len > get_bits_left(&s->gb))
  1024. return AVERROR_INVALIDDATA;
  1025. id = get_bits_long(&s->gb, 32);
  1026. id = av_be2ne32(id);
  1027. len -= 6;
  1028. if (s->avctx->debug & FF_DEBUG_STARTCODE)
  1029. av_log(s->avctx, AV_LOG_DEBUG, "APPx %8X\n", id);
  1030. /* Buggy AVID, it puts EOI only at every 10th frame. */
  1031. /* Also, this fourcc is used by non-avid files too, it holds some
  1032. information, but it's always present in AVID-created files. */
  1033. if (id == AV_RL32("AVI1")) {
  1034. /* structure:
  1035. 4bytes AVI1
  1036. 1bytes polarity
  1037. 1bytes always zero
  1038. 4bytes field_size
  1039. 4bytes field_size_less_padding
  1040. */
  1041. s->buggy_avid = 1;
  1042. // if (s->first_picture)
  1043. // printf("mjpeg: workarounding buggy AVID\n");
  1044. i = get_bits(&s->gb, 8);
  1045. if (i == 2)
  1046. s->bottom_field = 1;
  1047. else if (i == 1)
  1048. s->bottom_field = 0;
  1049. #if 0
  1050. skip_bits(&s->gb, 8);
  1051. skip_bits(&s->gb, 32);
  1052. skip_bits(&s->gb, 32);
  1053. len -= 10;
  1054. #endif
  1055. // if (s->interlace_polarity)
  1056. // printf("mjpeg: interlace polarity: %d\n", s->interlace_polarity);
  1057. goto out;
  1058. }
  1059. // len -= 2;
  1060. if (id == AV_RL32("JFIF")) {
  1061. int t_w, t_h, v1, v2;
  1062. skip_bits(&s->gb, 8); /* the trailing zero-byte */
  1063. v1 = get_bits(&s->gb, 8);
  1064. v2 = get_bits(&s->gb, 8);
  1065. skip_bits(&s->gb, 8);
  1066. s->avctx->sample_aspect_ratio.num = get_bits(&s->gb, 16);
  1067. s->avctx->sample_aspect_ratio.den = get_bits(&s->gb, 16);
  1068. if (s->avctx->debug & FF_DEBUG_PICT_INFO)
  1069. av_log(s->avctx, AV_LOG_INFO,
  1070. "mjpeg: JFIF header found (version: %x.%x) SAR=%d/%d\n",
  1071. v1, v2,
  1072. s->avctx->sample_aspect_ratio.num,
  1073. s->avctx->sample_aspect_ratio.den);
  1074. t_w = get_bits(&s->gb, 8);
  1075. t_h = get_bits(&s->gb, 8);
  1076. if (t_w && t_h) {
  1077. /* skip thumbnail */
  1078. if (len -10 - (t_w * t_h * 3) > 0)
  1079. len -= t_w * t_h * 3;
  1080. }
  1081. len -= 10;
  1082. goto out;
  1083. }
  1084. if (id == AV_RL32("Adob") && (get_bits(&s->gb, 8) == 'e')) {
  1085. if (s->avctx->debug & FF_DEBUG_PICT_INFO)
  1086. av_log(s->avctx, AV_LOG_INFO, "mjpeg: Adobe header found\n");
  1087. skip_bits(&s->gb, 16); /* version */
  1088. skip_bits(&s->gb, 16); /* flags0 */
  1089. skip_bits(&s->gb, 16); /* flags1 */
  1090. skip_bits(&s->gb, 8); /* transform */
  1091. len -= 7;
  1092. goto out;
  1093. }
  1094. if (id == AV_RL32("LJIF")) {
  1095. if (s->avctx->debug & FF_DEBUG_PICT_INFO)
  1096. av_log(s->avctx, AV_LOG_INFO,
  1097. "Pegasus lossless jpeg header found\n");
  1098. skip_bits(&s->gb, 16); /* version ? */
  1099. skip_bits(&s->gb, 16); /* unknwon always 0? */
  1100. skip_bits(&s->gb, 16); /* unknwon always 0? */
  1101. skip_bits(&s->gb, 16); /* unknwon always 0? */
  1102. switch (get_bits(&s->gb, 8)) {
  1103. case 1:
  1104. s->rgb = 1;
  1105. s->pegasus_rct = 0;
  1106. break;
  1107. case 2:
  1108. s->rgb = 1;
  1109. s->pegasus_rct = 1;
  1110. break;
  1111. default:
  1112. av_log(s->avctx, AV_LOG_ERROR, "unknown colorspace\n");
  1113. }
  1114. len -= 9;
  1115. goto out;
  1116. }
  1117. /* Apple MJPEG-A */
  1118. if ((s->start_code == APP1) && (len > (0x28 - 8))) {
  1119. id = get_bits_long(&s->gb, 32);
  1120. id = av_be2ne32(id);
  1121. len -= 4;
  1122. /* Apple MJPEG-A */
  1123. if (id == AV_RL32("mjpg")) {
  1124. #if 0
  1125. skip_bits(&s->gb, 32); /* field size */
  1126. skip_bits(&s->gb, 32); /* pad field size */
  1127. skip_bits(&s->gb, 32); /* next off */
  1128. skip_bits(&s->gb, 32); /* quant off */
  1129. skip_bits(&s->gb, 32); /* huff off */
  1130. skip_bits(&s->gb, 32); /* image off */
  1131. skip_bits(&s->gb, 32); /* scan off */
  1132. skip_bits(&s->gb, 32); /* data off */
  1133. #endif
  1134. if (s->avctx->debug & FF_DEBUG_PICT_INFO)
  1135. av_log(s->avctx, AV_LOG_INFO, "mjpeg: Apple MJPEG-A header found\n");
  1136. }
  1137. }
  1138. out:
  1139. /* slow but needed for extreme adobe jpegs */
  1140. if (len < 0)
  1141. av_log(s->avctx, AV_LOG_ERROR,
  1142. "mjpeg: error, decode_app parser read over the end\n");
  1143. while (--len > 0)
  1144. skip_bits(&s->gb, 8);
  1145. return 0;
  1146. }
  1147. static int mjpeg_decode_com(MJpegDecodeContext *s)
  1148. {
  1149. int len = get_bits(&s->gb, 16);
  1150. if (len >= 2 && 8 * len - 16 <= get_bits_left(&s->gb)) {
  1151. char *cbuf = av_malloc(len - 1);
  1152. if (cbuf) {
  1153. int i;
  1154. for (i = 0; i < len - 2; i++)
  1155. cbuf[i] = get_bits(&s->gb, 8);
  1156. if (i > 0 && cbuf[i - 1] == '\n')
  1157. cbuf[i - 1] = 0;
  1158. else
  1159. cbuf[i] = 0;
  1160. if (s->avctx->debug & FF_DEBUG_PICT_INFO)
  1161. av_log(s->avctx, AV_LOG_INFO, "mjpeg comment: '%s'\n", cbuf);
  1162. /* buggy avid, it puts EOI only at every 10th frame */
  1163. if (!strcmp(cbuf, "AVID")) {
  1164. s->buggy_avid = 1;
  1165. // if (s->first_picture)
  1166. // printf("mjpeg: workarounding buggy AVID\n");
  1167. } else if (!strcmp(cbuf, "CS=ITU601"))
  1168. s->cs_itu601 = 1;
  1169. else if ((len > 20 && !strncmp(cbuf, "Intel(R) JPEG Library", 21)) ||
  1170. (len > 19 && !strncmp(cbuf, "Metasoft MJPEG Codec", 20)))
  1171. s->flipped = 1;
  1172. av_free(cbuf);
  1173. }
  1174. }
  1175. return 0;
  1176. }
  1177. /* return the 8 bit start code value and update the search
  1178. state. Return -1 if no start code found */
  1179. static int find_marker(const uint8_t **pbuf_ptr, const uint8_t *buf_end)
  1180. {
  1181. const uint8_t *buf_ptr;
  1182. unsigned int v, v2;
  1183. int val;
  1184. #ifdef DEBUG
  1185. int skipped = 0;
  1186. #endif
  1187. buf_ptr = *pbuf_ptr;
  1188. while (buf_ptr < buf_end) {
  1189. v = *buf_ptr++;
  1190. v2 = *buf_ptr;
  1191. if ((v == 0xff) && (v2 >= 0xc0) && (v2 <= 0xfe) && buf_ptr < buf_end) {
  1192. val = *buf_ptr++;
  1193. goto found;
  1194. }
  1195. #ifdef DEBUG
  1196. skipped++;
  1197. #endif
  1198. }
  1199. val = -1;
  1200. found:
  1201. av_dlog(NULL, "find_marker skipped %d bytes\n", skipped);
  1202. *pbuf_ptr = buf_ptr;
  1203. return val;
  1204. }
  1205. int ff_mjpeg_find_marker(MJpegDecodeContext *s,
  1206. const uint8_t **buf_ptr, const uint8_t *buf_end,
  1207. const uint8_t **unescaped_buf_ptr,
  1208. int *unescaped_buf_size)
  1209. {
  1210. int start_code;
  1211. start_code = find_marker(buf_ptr, buf_end);
  1212. av_fast_padded_malloc(&s->buffer, &s->buffer_size, buf_end - *buf_ptr);
  1213. if (!s->buffer)
  1214. return AVERROR(ENOMEM);
  1215. /* unescape buffer of SOS, use special treatment for JPEG-LS */
  1216. if (start_code == SOS && !s->ls) {
  1217. const uint8_t *src = *buf_ptr;
  1218. uint8_t *dst = s->buffer;
  1219. while (src < buf_end) {
  1220. uint8_t x = *(src++);
  1221. *(dst++) = x;
  1222. if (s->avctx->codec_id != AV_CODEC_ID_THP) {
  1223. if (x == 0xff) {
  1224. while (src < buf_end && x == 0xff)
  1225. x = *(src++);
  1226. if (x >= 0xd0 && x <= 0xd7)
  1227. *(dst++) = x;
  1228. else if (x)
  1229. break;
  1230. }
  1231. }
  1232. }
  1233. *unescaped_buf_ptr = s->buffer;
  1234. *unescaped_buf_size = dst - s->buffer;
  1235. av_log(s->avctx, AV_LOG_DEBUG, "escaping removed %td bytes\n",
  1236. (buf_end - *buf_ptr) - (dst - s->buffer));
  1237. } else if (start_code == SOS && s->ls) {
  1238. const uint8_t *src = *buf_ptr;
  1239. uint8_t *dst = s->buffer;
  1240. int bit_count = 0;
  1241. int t = 0, b = 0;
  1242. PutBitContext pb;
  1243. s->cur_scan++;
  1244. /* find marker */
  1245. while (src + t < buf_end) {
  1246. uint8_t x = src[t++];
  1247. if (x == 0xff) {
  1248. while ((src + t < buf_end) && x == 0xff)
  1249. x = src[t++];
  1250. if (x & 0x80) {
  1251. t -= 2;
  1252. break;
  1253. }
  1254. }
  1255. }
  1256. bit_count = t * 8;
  1257. init_put_bits(&pb, dst, t);
  1258. /* unescape bitstream */
  1259. while (b < t) {
  1260. uint8_t x = src[b++];
  1261. put_bits(&pb, 8, x);
  1262. if (x == 0xFF) {
  1263. x = src[b++];
  1264. put_bits(&pb, 7, x);
  1265. bit_count--;
  1266. }
  1267. }
  1268. flush_put_bits(&pb);
  1269. *unescaped_buf_ptr = dst;
  1270. *unescaped_buf_size = (bit_count + 7) >> 3;
  1271. } else {
  1272. *unescaped_buf_ptr = *buf_ptr;
  1273. *unescaped_buf_size = buf_end - *buf_ptr;
  1274. }
  1275. return start_code;
  1276. }
  1277. int ff_mjpeg_decode_frame(AVCodecContext *avctx, void *data, int *data_size,
  1278. AVPacket *avpkt)
  1279. {
  1280. const uint8_t *buf = avpkt->data;
  1281. int buf_size = avpkt->size;
  1282. MJpegDecodeContext *s = avctx->priv_data;
  1283. const uint8_t *buf_end, *buf_ptr;
  1284. const uint8_t *unescaped_buf_ptr;
  1285. int unescaped_buf_size;
  1286. int start_code;
  1287. int ret = 0;
  1288. AVFrame *picture = data;
  1289. s->got_picture = 0; // picture from previous image can not be reused
  1290. buf_ptr = buf;
  1291. buf_end = buf + buf_size;
  1292. while (buf_ptr < buf_end) {
  1293. /* find start next marker */
  1294. start_code = ff_mjpeg_find_marker(s, &buf_ptr, buf_end,
  1295. &unescaped_buf_ptr,
  1296. &unescaped_buf_size);
  1297. /* EOF */
  1298. if (start_code < 0) {
  1299. goto the_end;
  1300. } else if (unescaped_buf_size > (1U<<29)) {
  1301. av_log(avctx, AV_LOG_ERROR, "MJPEG packet 0x%x too big (0x%x/0x%x), corrupt data?\n",
  1302. start_code, unescaped_buf_size, buf_size);
  1303. return AVERROR_INVALIDDATA;
  1304. } else {
  1305. av_log(avctx, AV_LOG_DEBUG, "marker=%x avail_size_in_buf=%td\n",
  1306. start_code, buf_end - buf_ptr);
  1307. init_get_bits(&s->gb, unescaped_buf_ptr, unescaped_buf_size * 8);
  1308. s->start_code = start_code;
  1309. if (s->avctx->debug & FF_DEBUG_STARTCODE)
  1310. av_log(avctx, AV_LOG_DEBUG, "startcode: %X\n", start_code);
  1311. /* process markers */
  1312. if (start_code >= 0xd0 && start_code <= 0xd7)
  1313. av_log(avctx, AV_LOG_DEBUG,
  1314. "restart marker: %d\n", start_code & 0x0f);
  1315. /* APP fields */
  1316. else if (start_code >= APP0 && start_code <= APP15)
  1317. mjpeg_decode_app(s);
  1318. /* Comment */
  1319. else if (start_code == COM)
  1320. mjpeg_decode_com(s);
  1321. switch (start_code) {
  1322. case SOI:
  1323. s->restart_interval = 0;
  1324. s->restart_count = 0;
  1325. /* nothing to do on SOI */
  1326. break;
  1327. case DQT:
  1328. ff_mjpeg_decode_dqt(s);
  1329. break;
  1330. case DHT:
  1331. if ((ret = ff_mjpeg_decode_dht(s)) < 0) {
  1332. av_log(avctx, AV_LOG_ERROR, "huffman table decode error\n");
  1333. return ret;
  1334. }
  1335. break;
  1336. case SOF0:
  1337. case SOF1:
  1338. s->lossless = 0;
  1339. s->ls = 0;
  1340. s->progressive = 0;
  1341. if ((ret = ff_mjpeg_decode_sof(s)) < 0)
  1342. return ret;
  1343. break;
  1344. case SOF2:
  1345. s->lossless = 0;
  1346. s->ls = 0;
  1347. s->progressive = 1;
  1348. if ((ret = ff_mjpeg_decode_sof(s)) < 0)
  1349. return ret;
  1350. break;
  1351. case SOF3:
  1352. s->lossless = 1;
  1353. s->ls = 0;
  1354. s->progressive = 0;
  1355. if ((ret = ff_mjpeg_decode_sof(s)) < 0)
  1356. return ret;
  1357. break;
  1358. case SOF48:
  1359. s->lossless = 1;
  1360. s->ls = 1;
  1361. s->progressive = 0;
  1362. if ((ret = ff_mjpeg_decode_sof(s)) < 0)
  1363. return ret;
  1364. break;
  1365. case LSE:
  1366. if (!CONFIG_JPEGLS_DECODER ||
  1367. (ret = ff_jpegls_decode_lse(s)) < 0)
  1368. return ret;
  1369. break;
  1370. case EOI:
  1371. s->cur_scan = 0;
  1372. if ((s->buggy_avid && !s->interlaced) || s->restart_interval)
  1373. break;
  1374. eoi_parser:
  1375. if (!s->got_picture) {
  1376. av_log(avctx, AV_LOG_WARNING,
  1377. "Found EOI before any SOF, ignoring\n");
  1378. break;
  1379. }
  1380. if (s->interlaced) {
  1381. s->bottom_field ^= 1;
  1382. /* if not bottom field, do not output image yet */
  1383. if (s->bottom_field == !s->interlace_polarity)
  1384. goto not_the_end;
  1385. }
  1386. *picture = *s->picture_ptr;
  1387. *data_size = sizeof(AVFrame);
  1388. if (!s->lossless) {
  1389. picture->quality = FFMAX3(s->qscale[0],
  1390. s->qscale[1],
  1391. s->qscale[2]);
  1392. picture->qstride = 0;
  1393. picture->qscale_table = s->qscale_table;
  1394. memset(picture->qscale_table, picture->quality,
  1395. (s->width + 15) / 16);
  1396. if (avctx->debug & FF_DEBUG_QP)
  1397. av_log(avctx, AV_LOG_DEBUG,
  1398. "QP: %d\n", picture->quality);
  1399. picture->quality *= FF_QP2LAMBDA;
  1400. }
  1401. goto the_end;
  1402. case SOS:
  1403. if (!s->got_picture) {
  1404. av_log(avctx, AV_LOG_WARNING,
  1405. "Can not process SOS before SOF, skipping\n");
  1406. break;
  1407. }
  1408. if ((ret = ff_mjpeg_decode_sos(s, NULL, NULL)) < 0 &&
  1409. (avctx->err_recognition & AV_EF_EXPLODE))
  1410. return ret;
  1411. /* buggy avid puts EOI every 10-20th frame */
  1412. /* if restart period is over process EOI */
  1413. if ((s->buggy_avid && !s->interlaced) || s->restart_interval)
  1414. goto eoi_parser;
  1415. break;
  1416. case DRI:
  1417. mjpeg_decode_dri(s);
  1418. break;
  1419. case SOF5:
  1420. case SOF6:
  1421. case SOF7:
  1422. case SOF9:
  1423. case SOF10:
  1424. case SOF11:
  1425. case SOF13:
  1426. case SOF14:
  1427. case SOF15:
  1428. case JPG:
  1429. av_log(avctx, AV_LOG_ERROR,
  1430. "mjpeg: unsupported coding type (%x)\n", start_code);
  1431. break;
  1432. // default:
  1433. // printf("mjpeg: unsupported marker (%x)\n", start_code);
  1434. // break;
  1435. }
  1436. not_the_end:
  1437. /* eof process start code */
  1438. buf_ptr += (get_bits_count(&s->gb) + 7) / 8;
  1439. av_log(avctx, AV_LOG_DEBUG,
  1440. "marker parser used %d bytes (%d bits)\n",
  1441. (get_bits_count(&s->gb) + 7) / 8, get_bits_count(&s->gb));
  1442. }
  1443. }
  1444. if (s->got_picture) {
  1445. av_log(avctx, AV_LOG_WARNING, "EOI missing, emulating\n");
  1446. goto eoi_parser;
  1447. }
  1448. av_log(avctx, AV_LOG_FATAL, "No JPEG data found in image\n");
  1449. return AVERROR_INVALIDDATA;
  1450. the_end:
  1451. av_log(avctx, AV_LOG_DEBUG, "mjpeg decode frame unused %td bytes\n",
  1452. buf_end - buf_ptr);
  1453. // return buf_end - buf_ptr;
  1454. return buf_ptr - buf;
  1455. }
  1456. av_cold int ff_mjpeg_decode_end(AVCodecContext *avctx)
  1457. {
  1458. MJpegDecodeContext *s = avctx->priv_data;
  1459. int i, j;
  1460. if (s->picture_ptr && s->picture_ptr->data[0])
  1461. avctx->release_buffer(avctx, s->picture_ptr);
  1462. av_free(s->buffer);
  1463. av_free(s->qscale_table);
  1464. av_freep(&s->ljpeg_buffer);
  1465. s->ljpeg_buffer_size = 0;
  1466. for (i = 0; i < 3; i++) {
  1467. for (j = 0; j < 4; j++)
  1468. ff_free_vlc(&s->vlcs[i][j]);
  1469. }
  1470. for (i = 0; i < MAX_COMPONENTS; i++) {
  1471. av_freep(&s->blocks[i]);
  1472. av_freep(&s->last_nnz[i]);
  1473. }
  1474. return 0;
  1475. }
  1476. #define OFFSET(x) offsetof(MJpegDecodeContext, x)
  1477. #define VD AV_OPT_FLAG_VIDEO_PARAM | AV_OPT_FLAG_DECODING_PARAM
  1478. static const AVOption options[] = {
  1479. { "extern_huff", "Use external huffman table.",
  1480. OFFSET(extern_huff), AV_OPT_TYPE_INT, { .i64 = 0 }, 0, 1, VD },
  1481. { NULL },
  1482. };
  1483. static const AVClass mjpegdec_class = {
  1484. .class_name = "MJPEG decoder",
  1485. .item_name = av_default_item_name,
  1486. .option = options,
  1487. .version = LIBAVUTIL_VERSION_INT,
  1488. };
  1489. AVCodec ff_mjpeg_decoder = {
  1490. .name = "mjpeg",
  1491. .type = AVMEDIA_TYPE_VIDEO,
  1492. .id = AV_CODEC_ID_MJPEG,
  1493. .priv_data_size = sizeof(MJpegDecodeContext),
  1494. .init = ff_mjpeg_decode_init,
  1495. .close = ff_mjpeg_decode_end,
  1496. .decode = ff_mjpeg_decode_frame,
  1497. .capabilities = CODEC_CAP_DR1,
  1498. .long_name = NULL_IF_CONFIG_SMALL("MJPEG (Motion JPEG)"),
  1499. .priv_class = &mjpegdec_class,
  1500. };
  1501. AVCodec ff_thp_decoder = {
  1502. .name = "thp",
  1503. .type = AVMEDIA_TYPE_VIDEO,
  1504. .id = AV_CODEC_ID_THP,
  1505. .priv_data_size = sizeof(MJpegDecodeContext),
  1506. .init = ff_mjpeg_decode_init,
  1507. .close = ff_mjpeg_decode_end,
  1508. .decode = ff_mjpeg_decode_frame,
  1509. .capabilities = CODEC_CAP_DR1,
  1510. .long_name = NULL_IF_CONFIG_SMALL("Nintendo Gamecube THP video"),
  1511. };