50 const uint8_t *val_table,
int nb_codes,
51 int use_static,
int is_ac)
54 uint16_t huff_code[256];
55 uint16_t huff_sym[256];
62 for (i = 0; i < 256; i++)
63 huff_sym[i] = i + 16 * is_ac;
66 huff_sym[0] = 16 * 256;
69 huff_code, 2, 2, huff_sym, 2, 2, use_static);
91 if (len > 14 && buf[12] == 1)
93 if (len > 14 && buf[12] == 2)
132 "error using external huffman table, switching back to internal\n");
175 for (i = 0; i < 64; i++) {
210 for (i = 1; i <= 16; i++) {
215 if (len < n || n > 256)
219 for (i = 0; i <
n; i++) {
230 class, index, code_max + 1);
231 if ((ret =
build_vlc(&s->
vlcs[
class][index], bits_table, val_table,
232 code_max + 1, 0,
class > 0)) < 0)
237 if ((ret =
build_vlc(&s->
vlcs[2][index], bits_table, val_table,
238 code_max + 1, 0, 0)) < 0)
260 if (bits > 16 || bits < 1) {
290 if (nb_components <= 0 ||
296 "nb_components changing in interlaced picture\n");
300 if (s->
ls && !(bits <= 8 || nb_components == 1)) {
302 "JPEG-LS that is not <= 8 "
303 "bits/component or 16-bit gray");
309 memset(h_count, 0,
sizeof(h_count));
310 memset(v_count, 0,
sizeof(v_count));
311 for (i = 0; i < nb_components; i++) {
317 if (h_count[i] > s->
h_max)
318 s->
h_max = h_count[i];
319 if (v_count[i] > s->
v_max)
320 s->
v_max = v_count[i];
326 if (!h_count[i] || !v_count[i]) {
328 "Invalid sampling factor in component %d %d:%d\n",
329 i, h_count[i], v_count[i]);
334 i, h_count[i], v_count[i],
347 || memcmp(s->
h_count, h_count,
sizeof(h_count))
348 || memcmp(s->
v_count, v_count,
sizeof(v_count))) {
353 memcpy(s->
h_count, h_count,
sizeof(h_count));
354 memcpy(s->
v_count, v_count,
sizeof(v_count));
382 if (s->
v_max == 1 && s->
h_max == 1 && s->
lossless==1 && (nb_components==3 || nb_components==4))
387 pix_fmt_id = ((unsigned)s->
h_count[0] << 28) | (s->
v_count[0] << 24) |
394 if (!(pix_fmt_id & 0xD0D0D0D0))
395 pix_fmt_id -= (pix_fmt_id & 0xF0F0F0F0) >> 1;
396 if (!(pix_fmt_id & 0x0D0D0D0D))
397 pix_fmt_id -= (pix_fmt_id & 0x0F0F0F0F) >> 1;
399 for (i = 0; i < 8; i++) {
400 int j = 6 + (i&1) - (i&6);
401 int is = (pix_fmt_id >> (4*i)) & 0xF;
402 int js = (pix_fmt_id >> (4*j)) & 0xF;
404 if (is == 1 && js != 2 && (i < 2 || i > 5))
405 js = (pix_fmt_id >> ( 8 + 4*(i&1))) & 0xF;
406 if (is == 1 && js != 2 && (i < 2 || i > 5))
407 js = (pix_fmt_id >> (16 + 4*(i&1))) & 0xF;
409 if (is == 1 && js == 2) {
415 switch (pix_fmt_id) {
503 if (pix_fmt_id == 0x14111100)
536 if (pix_fmt_id == 0x42111100) {
540 }
else if (pix_fmt_id == 0x24111100) {
568 else if (s->
bits <= 8)
587 for (i = 0; i < 4; i++)
594 if (len != (8 + (3 * nb_components)))
605 int bw = (width + s->
h_max * 8 - 1) / (s->
h_max * 8);
606 int bh = (height + s->
v_max * 8 - 1) / (s->
v_max * 8);
626 if (code < 0 || code > 16) {
628 "mjpeg_decode_dc: bad vlc: %d:%d (%p)\n",
629 0, dc_index, &s->
vlcs[0][dc_index]);
641 int dc_index,
int ac_index, int16_t *quant_matrix)
647 if (val == 0xfffff) {
651 val = val * quant_matrix[0] + s->
last_dc[component];
661 i += ((unsigned)code) >> 4;
669 int sign = (~cache) >> 31;
670 level = (
NEG_USR32(sign ^ cache,code) ^ sign) - sign;
680 block[j] = level * quant_matrix[j];
689 int component,
int dc_index,
690 int16_t *quant_matrix,
int Al)
695 if (val == 0xfffff) {
699 val = (val * quant_matrix[0] << Al) + s->
last_dc[component];
707 uint8_t *last_nnz,
int ac_index,
708 int16_t *quant_matrix,
709 int ss,
int se,
int Al,
int *EOBRUN)
720 for (i = ss; ; i++) {
724 run = ((unsigned) code) >> 4;
733 int sign = (~cache) >> 31;
734 level = (
NEG_USR32(sign ^ cache,code) ^ sign) - sign;
742 block[j] = level * quant_matrix[j] << Al;
749 block[j] = level * quant_matrix[j] << Al;
778 #define REFINE_BIT(j) { \
779 UPDATE_CACHE(re, &s->gb); \
780 sign = block[j] >> 15; \
781 block[j] += SHOW_UBITS(re, &s->gb, 1) * \
782 ((quant_matrix[j] ^ sign) - sign) << Al; \
783 LAST_SKIP_BITS(re, &s->gb, 1); \
791 av_log(s->avctx, AV_LOG_ERROR, "error count: %d\n", i); \
796 j = s->scantable.permutated[i]; \
799 else if (run-- == 0) \
806 int ac_index, int16_t *quant_matrix,
807 int ss,
int se,
int Al,
int *EOBRUN)
809 int code, i = ss, j, sign,
val,
run;
810 int last =
FFMIN(se, *last_nnz);
821 run = ((unsigned) code) >> 4;
828 block[j] = ((quant_matrix[j]^
val) - val) << Al;
836 run = ((unsigned) code) >> 4;
857 for (; i <= last; i++) {
878 for (i = 0; i < nb_components; i++)
892 for (i = 0; i < nb_components; i++)
907 int left[4], top[4], topleft[4];
908 const int linesize = s->
linesize[0];
909 const int mask = ((1 << s->
bits) - 1) << point_transform;
925 for (i = 0; i < 4; i++)
928 for (mb_y = 0; mb_y < s->
mb_height; mb_y++) {
932 ptr += linesize >> 1;
934 for (i = 0; i < 4; i++)
935 top[i] = left[i] = topleft[i] =
buffer[0][i];
937 for (mb_x = 0; mb_x < s->
mb_width; mb_x++) {
945 top[i] = left[i]= topleft[i]= 1 << (s->
bits - 1);
947 if (mb_y == resync_mb_y || mb_y == resync_mb_y+1 && mb_x < resync_mb_x || !mb_x)
948 modified_predictor = 1;
950 for (i=0;i<nb_components;i++) {
956 PREDICT(pred, topleft[i], top[i], left[i], modified_predictor);
962 left[i] =
buffer[mb_x][i] =
963 mask & (pred + (dc << point_transform));
972 for(i=0; i<nb_components; i++) {
975 for(mb_x = 0; mb_x < s->
mb_width; mb_x++) {
976 ptr[4*mb_x+3-
c] =
buffer[mb_x][i];
978 }
else if(s->
bits == 9) {
981 for(mb_x = 0; mb_x < s->
mb_width; mb_x++) {
982 ((uint16_t*)ptr)[4*mb_x+
c] =
buffer[mb_x][i];
987 for (mb_x = 0; mb_x < s->
mb_width; mb_x++) {
988 ptr[3*mb_x + 1] =
buffer[mb_x][0] - ((
buffer[mb_x][1] +
buffer[mb_x][2] - 0x200) >> 2);
989 ptr[3*mb_x + 0] =
buffer[mb_x][1] + ptr[3*mb_x + 1];
990 ptr[3*mb_x + 2] =
buffer[mb_x][2] + ptr[3*mb_x + 1];
993 for (mb_x = 0; mb_x < s->
mb_width; mb_x++) {
995 ptr[3*mb_x + 0] =
buffer[mb_x][1] + ptr[3*mb_x + 1];
996 ptr[3*mb_x + 2] =
buffer[mb_x][2] + ptr[3*mb_x + 1];
999 for(i=0; i<nb_components; i++) {
1002 for(mb_x = 0; mb_x < s->
mb_width; mb_x++) {
1003 ptr[3*mb_x+2-
c] =
buffer[mb_x][i];
1005 }
else if(s->
bits == 9) {
1008 for(mb_x = 0; mb_x < s->
mb_width; mb_x++) {
1009 ((uint16_t*)ptr)[3*mb_x+2-
c] =
buffer[mb_x][i];
1019 int point_transform,
int nb_components)
1021 int i, mb_x, mb_y,
mask;
1023 int resync_mb_y = 0;
1024 int resync_mb_x = 0;
1026 point_transform += bits - s->
bits;
1027 mask = ((1 << s->
bits) - 1) << point_transform;
1029 av_assert0(nb_components>=1 && nb_components<=4);
1031 for (mb_y = 0; mb_y < s->
mb_height; mb_y++) {
1032 for (mb_x = 0; mb_x < s->
mb_width; mb_x++) {
1039 if(!mb_x || mb_y == resync_mb_y || mb_y == resync_mb_y+1 && mb_x < resync_mb_x || s->
interlaced){
1040 int toprow = mb_y == resync_mb_y || mb_y == resync_mb_y+1 && mb_x < resync_mb_x;
1041 int leftcol = !mb_x || mb_y == resync_mb_y && mb_x == resync_mb_x;
1042 for (i = 0; i < nb_components; i++) {
1045 int n, h,
v, x,
y,
c, j, linesize;
1054 if(bits>8) linesize /= 2;
1056 for(j=0; j<
n; j++) {
1065 if(x==0 && leftcol){
1066 pred= 1 << (bits - 1);
1071 if(x==0 && leftcol){
1072 pred= ptr[-linesize];
1074 PREDICT(pred, ptr[-linesize-1], ptr[-linesize], ptr[-1], predictor);
1079 ptr += linesize >> 1;
1081 *ptr= pred + (dc << point_transform);
1083 ptr16 = (uint16_t*)(s->
picture_ptr->
data[c] + 2*(linesize * (v * mb_y + y)) + 2*(h * mb_x + x));
1085 if(x==0 && leftcol){
1086 pred= 1 << (bits - 1);
1091 if(x==0 && leftcol){
1092 pred= ptr16[-linesize];
1094 PREDICT(pred, ptr16[-linesize-1], ptr16[-linesize], ptr16[-1], predictor);
1099 ptr16 += linesize >> 1;
1101 *ptr16= pred + (dc << point_transform);
1110 for (i = 0; i < nb_components; i++) {
1113 int n, h,
v, x,
y,
c, j, linesize,
dc;
1122 if(bits>8) linesize /= 2;
1124 for (j = 0; j <
n; j++) {
1132 (linesize * (v * mb_y +
y)) +
1134 PREDICT(pred, ptr[-linesize-1], ptr[-linesize], ptr[-1], predictor);
1137 *ptr = pred + (dc << point_transform);
1139 ptr16 = (uint16_t*)(s->
picture_ptr->
data[c] + 2*(linesize * (v * mb_y + y)) + 2*(h * mb_x + x));
1140 PREDICT(pred, ptr16[-linesize-1], ptr16[-linesize], ptr16[-1], predictor);
1143 *ptr16= pred + (dc << point_transform);
1164 int linesize,
int lowres)
1169 case 1:
copy_block4(dst, src, linesize, linesize, 4);
1171 case 2:
copy_block2(dst, src, linesize, linesize, 2);
1173 case 3: *dst = *
src;
1180 int block_x, block_y;
1183 for (block_y=0; block_y<
size; block_y++)
1184 for (block_x=0; block_x<
size; block_x++)
1185 *(uint16_t*)(ptr + 2*block_x + block_y*linesize) <<= 16 - s->
bits;
1187 for (block_y=0; block_y<
size; block_y++)
1188 for (block_x=0; block_x<
size; block_x++)
1189 *(ptr + block_x + block_y*linesize) <<= 8 - s->
bits;
1194 int Al,
const uint8_t *mb_bitmask,
1195 int mb_bitmask_size,
1203 int bytes_per_pixel = 1 + (s->
bits > 8);
1215 for (i = 0; i < nb_components; i++) {
1218 reference_data[
c] = reference ? reference->
data[
c] : NULL;
1223 for (mb_y = 0; mb_y < s->
mb_height; mb_y++) {
1224 for (mb_x = 0; mb_x < s->
mb_width; mb_x++) {
1235 for (i = 0; i < nb_components; i++) {
1237 int n, h,
v, x,
y,
c, j;
1245 for (j = 0; j <
n; j++) {
1246 block_offset = (((linesize[
c] * (v * mb_y +
y) * 8) +
1247 (h * mb_x + x) * 8 * bytes_per_pixel) >> s->
avctx->
lowres);
1250 block_offset += linesize[
c] >> 1;
1251 ptr = data[
c] + block_offset;
1263 "error y=%d x=%d\n", mb_y, mb_x);
1281 "error y=%d x=%d\n", mb_y, mb_x);
1285 av_dlog(s->
avctx,
"mb: %d %d processed\n", mb_y, mb_x);
1288 (v * mb_y + y) * 8, (h * mb_x + x) * 8);
1303 int se,
int Ah,
int Al)
1312 int bytes_per_pixel = 1 + (s->
bits > 8);
1315 if (se < ss || se > 63) {
1328 data += linesize >> 1;
1332 for (mb_y = 0; mb_y < s->
mb_height; mb_y++) {
1337 for (mb_x = 0; mb_x < s->
mb_width; mb_x++,
block++, last_nnz++) {
1344 quant_matrix, ss, se, Al, &EOBRUN);
1347 quant_matrix, ss, se, Al, &EOBRUN);
1350 "error y=%d x=%d\n", mb_y, mb_x);
1368 int mb_bitmask_size,
const AVFrame *reference)
1370 int len, nb_components, i, h,
v,
predictor, point_transform;
1372 const int block_size = s->
lossless ? 1 : 8;
1373 int ilv, prev_shift;
1377 "Can not process SOS before SOF, skipping\n");
1387 "decode_sos: nb_components (%d) unsupported\n", nb_components);
1390 if (len != 6 + 2 * nb_components) {
1394 for (i = 0; i < nb_components; i++) {
1403 "decode_sos: index(%d) out of components\n", index);
1419 index = (index+2)%3;
1439 prev_shift = point_transform = 0;
1441 if (nb_components > 1) {
1445 }
else if (!s->
ls) {
1448 s->
mb_width = (s->
width + h * block_size - 1) / (h * block_size);
1457 s->
lossless ?
"lossless" :
"sequential DCT", s->
rgb ?
"RGB" :
"",
1467 for (i = 0; i < nb_components; i++)
1472 if (CONFIG_JPEGLS_DECODER && s->
ls) {
1477 point_transform, ilv)) < 0)
1486 nb_components)) < 0)
1495 point_transform)) < 0)
1499 prev_shift, point_transform,
1500 mb_bitmask, mb_bitmask_size, reference)) < 0)
1584 int t_w, t_h, v1, v2;
1596 "mjpeg: JFIF header found (version: %x.%x) SAR=%d/%d\n",
1605 if (len -10 - (t_w * t_h * 3) > 0)
1606 len -= t_w * t_h * 3;
1628 "Pegasus lossless jpeg header found\n");
1658 if (
id ==
AV_RL32(
"colr") && len > 0) {
1665 if (
id ==
AV_RL32(
"xfrm") && len > 0) {
1693 }
else if (type == 1) {
1705 if (!(flags & 0x04)) {
1715 int ret,
le, ifd_offset, bytes_read;
1772 "mjpeg: error, decode_app parser read over the end\n");
1786 for (i = 0; i < len - 2; i++)
1788 if (i > 0 && cbuf[i - 1] ==
'\n')
1797 if (!strncmp(cbuf,
"AVID", 4)) {
1799 }
else if (!strcmp(cbuf,
"CS=ITU601"))
1801 else if ((!strncmp(cbuf,
"Intel(R) JPEG Library, version 1", 32) && s->
avctx->
codec_tag) ||
1802 (!strncmp(cbuf,
"Metasoft MJPEG Codec", 20)))
1821 buf_ptr = *pbuf_ptr;
1822 while (buf_end - buf_ptr > 1) {
1825 if ((v == 0xff) && (v2 >= 0xc0) && (v2 <= 0xfe) && buf_ptr < buf_end) {
1834 av_dlog(NULL,
"find_marker skipped %d bytes\n", skipped);
1835 *pbuf_ptr = buf_ptr;
1841 const uint8_t **unescaped_buf_ptr,
1842 int *unescaped_buf_size)
1852 if (start_code ==
SOS && !s->
ls) {
1856 while (src < buf_end) {
1862 while (src < buf_end && x == 0xff)
1865 if (x >= 0xd0 && x <= 0xd7)
1872 *unescaped_buf_ptr = s->
buffer;
1873 *unescaped_buf_size = dst - s->
buffer;
1874 memset(s->
buffer + *unescaped_buf_size, 0,
1878 (buf_end - *buf_ptr) - (dst - s->
buffer));
1879 }
else if (start_code ==
SOS && s->
ls) {
1887 while (src + t < buf_end) {
1890 while ((src + t < buf_end) && x == 0xff)
1913 *unescaped_buf_ptr = dst;
1914 *unescaped_buf_size = (bit_count + 7) >> 3;
1915 memset(s->
buffer + *unescaped_buf_size, 0,
1918 *unescaped_buf_ptr = *buf_ptr;
1919 *unescaped_buf_size = buf_end - *buf_ptr;
1930 int buf_size = avpkt->
size;
1932 const uint8_t *buf_end, *buf_ptr;
1933 const uint8_t *unescaped_buf_ptr;
1935 int unescaped_buf_size;
1946 buf_end = buf + buf_size;
1947 while (buf_ptr < buf_end) {
1951 &unescaped_buf_size);
1953 if (start_code < 0) {
1955 }
else if (unescaped_buf_size > INT_MAX / 8) {
1957 "MJPEG packet 0x%x too big (%d/%d), corrupt data?\n",
1958 start_code, unescaped_buf_size, buf_size);
1962 start_code, buf_end - buf_ptr);
1976 if (start_code >= 0xd0 && start_code <= 0xd7)
1978 "restart marker: %d\n", start_code & 0x0f);
1980 else if (start_code >=
APP0 && start_code <=
APP15)
1983 else if (start_code ==
COM)
1988 if (!CONFIG_JPEGLS_DECODER &&
1989 (start_code ==
SOF48 || start_code ==
LSE)) {
1994 switch (start_code) {
2039 if (!CONFIG_JPEGLS_DECODER ||
2048 "Found EOI before any SOF, ignoring\n");
2066 int qpw = (s->
width + 15) / 16;
2069 memset(qp_table_buf->
data, qp, qpw);
2098 "mjpeg: unsupported coding type (%x)\n", start_code);
2105 "marker parser used %d bytes (%d bits)\n",
2137 for (p = 0; p<4; p++) {
2150 for (i = 0; i < h; i++) {
2151 if (is16bit) ((uint16_t*)line)[w - 1] = ((uint16_t*)line)[(w - 1) / 2];
2152 else line[w - 1] = line[(w - 1) / 2];
2153 for (index = w - 2; index > 0; index--) {
2155 ((uint16_t*)line)[
index] = (((uint16_t*)line)[index / 2] + ((uint16_t*)line)[(index + 1) / 2]) >> 1;
2157 line[
index] = (line[index / 2] + line[(index + 1) / 2]) >> 1;
2180 for (p = 0; p < 4; p++) {
2191 for (i = h - 1; i; i--) {
2194 if (src1 == src2 || i == h - 1) {
2195 memcpy(dst, src1, w);
2197 for (index = 0; index < w; index++)
2198 dst[index] = (src1[index] + src2[index]) >> 1;
2207 for (index=0; index<4; index++) {
2211 if(index && index<3){
2217 for (i=0; i<h/2; i++) {
2219 FFSWAP(
int, dst[j], dst2[j]);
2229 for (i=0; i<h; i++) {
2232 for (index=0; index<4; index++) {
2236 for (j=0; j<w; j++) {
2238 int r = dst[0][j] * k;
2239 int g = dst[1][j] * k;
2240 int b = dst[2][j] * k;
2241 dst[0][j] = g*257 >> 16;
2242 dst[1][j] = b*257 >> 16;
2243 dst[2][j] = r*257 >> 16;
2251 for (i=0; i<h; i++) {
2254 for (index=0; index<4; index++) {
2258 for (j=0; j<w; j++) {
2260 int r = (255 - dst[0][j]) * k;
2261 int g = (128 - dst[1][j]) * k;
2262 int b = (128 - dst[2][j]) * k;
2263 dst[0][j] = r*257 >> 16;
2264 dst[1][j] = (g*257 >> 16) + 128;
2265 dst[2][j] = (b*257 >> 16) + 128;
2286 return buf_ptr -
buf;
2309 for (i = 0; i < 3; i++) {
2310 for (j = 0; j < 4; j++)
2327 #if CONFIG_MJPEG_DECODER
2328 #define OFFSET(x) offsetof(MJpegDecodeContext, x)
2329 #define VD AV_OPT_FLAG_VIDEO_PARAM | AV_OPT_FLAG_DECODING_PARAM
2331 {
"extern_huff",
"Use external huffman table.",
2336 static const AVClass mjpegdec_class = {
2355 .priv_class = &mjpegdec_class,
2358 #if CONFIG_THP_DECODER