Merge commit '4978850ca2cb1ec6908f5bc79cc592ca454d11e8'
[ffmpeg.git] / libavcodec / mjpegenc_common.c
1 /*
2  * lossless JPEG shared bits
3  * Copyright (c) 2000, 2001 Fabrice Bellard
4  * Copyright (c) 2003 Alex Beregszaszi
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 #include <stdint.h>
24 #include <string.h>
25
26 #include "libavutil/common.h"
27 #include "libavutil/pixdesc.h"
28 #include "libavutil/pixfmt.h"
29
30 #include "avcodec.h"
31 #include "idctdsp.h"
32 #include "jpegtables.h"
33 #include "put_bits.h"
34 #include "mjpegenc.h"
35 #include "mjpegenc_common.h"
36 #include "mjpeg.h"
37
38 /* table_class: 0 = DC coef, 1 = AC coefs */
39 static int put_huffman_table(PutBitContext *p, int table_class, int table_id,
40                              const uint8_t *bits_table, const uint8_t *value_table)
41 {
42     int n, i;
43
44     put_bits(p, 4, table_class);
45     put_bits(p, 4, table_id);
46
47     n = 0;
48     for(i=1;i<=16;i++) {
49         n += bits_table[i];
50         put_bits(p, 8, bits_table[i]);
51     }
52
53     for(i=0;i<n;i++)
54         put_bits(p, 8, value_table[i]);
55
56     return n + 17;
57 }
58
59 static void jpeg_table_header(AVCodecContext *avctx, PutBitContext *p,
60                               ScanTable *intra_scantable,
61                               uint16_t luma_intra_matrix[64],
62                               uint16_t chroma_intra_matrix[64],
63                               int hsample[3])
64 {
65     int i, j, size;
66     uint8_t *ptr;
67
68     if (avctx->codec_id != AV_CODEC_ID_LJPEG) {
69         int matrix_count = 1 + !!memcmp(luma_intra_matrix,
70                                         chroma_intra_matrix,
71                                         sizeof(luma_intra_matrix[0]) * 64);
72     /* quant matrixes */
73     put_marker(p, DQT);
74     put_bits(p, 16, 2 + matrix_count * (1 + 64));
75     put_bits(p, 4, 0); /* 8 bit precision */
76     put_bits(p, 4, 0); /* table 0 */
77     for(i=0;i<64;i++) {
78         j = intra_scantable->permutated[i];
79         put_bits(p, 8, luma_intra_matrix[j]);
80     }
81
82         if (matrix_count > 1) {
83             put_bits(p, 4, 0); /* 8 bit precision */
84             put_bits(p, 4, 1); /* table 1 */
85             for(i=0;i<64;i++) {
86                 j = intra_scantable->permutated[i];
87                 put_bits(p, 8, chroma_intra_matrix[j]);
88             }
89         }
90     }
91
92     if(avctx->active_thread_type & FF_THREAD_SLICE){
93         put_marker(p, DRI);
94         put_bits(p, 16, 4);
95         put_bits(p, 16, (avctx->width-1)/(8*hsample[0]) + 1);
96     }
97
98     /* huffman table */
99     put_marker(p, DHT);
100     flush_put_bits(p);
101     ptr = put_bits_ptr(p);
102     put_bits(p, 16, 0); /* patched later */
103     size = 2;
104     size += put_huffman_table(p, 0, 0, avpriv_mjpeg_bits_dc_luminance,
105                               avpriv_mjpeg_val_dc);
106     size += put_huffman_table(p, 0, 1, avpriv_mjpeg_bits_dc_chrominance,
107                               avpriv_mjpeg_val_dc);
108
109     size += put_huffman_table(p, 1, 0, avpriv_mjpeg_bits_ac_luminance,
110                               avpriv_mjpeg_val_ac_luminance);
111     size += put_huffman_table(p, 1, 1, avpriv_mjpeg_bits_ac_chrominance,
112                               avpriv_mjpeg_val_ac_chrominance);
113     AV_WB16(ptr, size);
114 }
115
116 static void jpeg_put_comments(AVCodecContext *avctx, PutBitContext *p)
117 {
118     int size;
119     uint8_t *ptr;
120
121     if (avctx->sample_aspect_ratio.num > 0 && avctx->sample_aspect_ratio.den > 0) {
122         /* JFIF header */
123         put_marker(p, APP0);
124         put_bits(p, 16, 16);
125         avpriv_put_string(p, "JFIF", 1); /* this puts the trailing zero-byte too */
126         put_bits(p, 16, 0x0102);         /* v 1.02 */
127         put_bits(p,  8, 0);              /* units type: 0 - aspect ratio */
128         put_bits(p, 16, avctx->sample_aspect_ratio.num);
129         put_bits(p, 16, avctx->sample_aspect_ratio.den);
130         put_bits(p, 8, 0); /* thumbnail width */
131         put_bits(p, 8, 0); /* thumbnail height */
132     }
133
134     /* comment */
135     if (!(avctx->flags & CODEC_FLAG_BITEXACT)) {
136         put_marker(p, COM);
137         flush_put_bits(p);
138         ptr = put_bits_ptr(p);
139         put_bits(p, 16, 0); /* patched later */
140         avpriv_put_string(p, LIBAVCODEC_IDENT, 1);
141         size = strlen(LIBAVCODEC_IDENT)+3;
142         AV_WB16(ptr, size);
143     }
144
145     if (avctx->pix_fmt == AV_PIX_FMT_YUV420P ||
146         avctx->pix_fmt == AV_PIX_FMT_YUV422P ||
147         avctx->pix_fmt == AV_PIX_FMT_YUV444P) {
148         put_marker(p, COM);
149         flush_put_bits(p);
150         ptr = put_bits_ptr(p);
151         put_bits(p, 16, 0); /* patched later */
152         avpriv_put_string(p, "CS=ITU601", 1);
153         size = strlen("CS=ITU601")+3;
154         AV_WB16(ptr, size);
155     }
156 }
157
158 void ff_mjpeg_init_hvsample(AVCodecContext *avctx, int hsample[3], int vsample[3])
159 {
160     int chroma_h_shift, chroma_v_shift;
161
162     av_pix_fmt_get_chroma_sub_sample(avctx->pix_fmt, &chroma_h_shift,
163                                      &chroma_v_shift);
164     if (avctx->codec->id == AV_CODEC_ID_LJPEG &&
165         (   avctx->pix_fmt == AV_PIX_FMT_BGR0
166          || avctx->pix_fmt == AV_PIX_FMT_BGRA
167          || avctx->pix_fmt == AV_PIX_FMT_BGR24)) {
168         vsample[0] = hsample[0] =
169         vsample[1] = hsample[1] =
170         vsample[2] = hsample[2] = 1;
171     } else if (avctx->pix_fmt == AV_PIX_FMT_YUV444P || avctx->pix_fmt == AV_PIX_FMT_YUVJ444P) {
172         vsample[0] = vsample[1] = vsample[2] = 2;
173         hsample[0] = hsample[1] = hsample[2] = 1;
174     } else {
175         vsample[0] = 2;
176         vsample[1] = 2 >> chroma_v_shift;
177         vsample[2] = 2 >> chroma_v_shift;
178         hsample[0] = 2;
179         hsample[1] = 2 >> chroma_h_shift;
180         hsample[2] = 2 >> chroma_h_shift;
181     }
182 }
183
184 void ff_mjpeg_encode_picture_header(AVCodecContext *avctx, PutBitContext *pb,
185                                     ScanTable *intra_scantable,
186                                     uint16_t luma_intra_matrix[64],
187                                     uint16_t chroma_intra_matrix[64])
188 {
189     const int lossless = avctx->codec_id != AV_CODEC_ID_MJPEG && avctx->codec_id != AV_CODEC_ID_AMV;
190     int hsample[3], vsample[3];
191     int i;
192     int chroma_matrix = !!memcmp(luma_intra_matrix,
193                                  chroma_intra_matrix,
194                                  sizeof(luma_intra_matrix[0])*64);
195
196     ff_mjpeg_init_hvsample(avctx, hsample, vsample);
197
198     put_marker(pb, SOI);
199
200     // hack for AMV mjpeg format
201     if(avctx->codec_id == AV_CODEC_ID_AMV) goto end;
202
203     jpeg_put_comments(avctx, pb);
204
205     jpeg_table_header(avctx, pb, intra_scantable, luma_intra_matrix, chroma_intra_matrix, hsample);
206
207     switch (avctx->codec_id) {
208     case AV_CODEC_ID_MJPEG:  put_marker(pb, SOF0 ); break;
209     case AV_CODEC_ID_LJPEG:  put_marker(pb, SOF3 ); break;
210     default: av_assert0(0);
211     }
212
213     put_bits(pb, 16, 17);
214     if (lossless && (  avctx->pix_fmt == AV_PIX_FMT_BGR0
215                     || avctx->pix_fmt == AV_PIX_FMT_BGRA
216                     || avctx->pix_fmt == AV_PIX_FMT_BGR24))
217         put_bits(pb, 8, 9); /* 9 bits/component RCT */
218     else
219         put_bits(pb, 8, 8); /* 8 bits/component */
220     put_bits(pb, 16, avctx->height);
221     put_bits(pb, 16, avctx->width);
222     put_bits(pb, 8, 3); /* 3 components */
223
224     /* Y component */
225     put_bits(pb, 8, 1); /* component number */
226     put_bits(pb, 4, hsample[0]); /* H factor */
227     put_bits(pb, 4, vsample[0]); /* V factor */
228     put_bits(pb, 8, 0); /* select matrix */
229
230     /* Cb component */
231     put_bits(pb, 8, 2); /* component number */
232     put_bits(pb, 4, hsample[1]); /* H factor */
233     put_bits(pb, 4, vsample[1]); /* V factor */
234     put_bits(pb, 8, lossless ? 0 : chroma_matrix); /* select matrix */
235
236     /* Cr component */
237     put_bits(pb, 8, 3); /* component number */
238     put_bits(pb, 4, hsample[2]); /* H factor */
239     put_bits(pb, 4, vsample[2]); /* V factor */
240     put_bits(pb, 8, lossless ? 0 : chroma_matrix); /* select matrix */
241
242     /* scan header */
243     put_marker(pb, SOS);
244     put_bits(pb, 16, 12); /* length */
245     put_bits(pb, 8, 3); /* 3 components */
246
247     /* Y component */
248     put_bits(pb, 8, 1); /* index */
249     put_bits(pb, 4, 0); /* DC huffman table index */
250     put_bits(pb, 4, 0); /* AC huffman table index */
251
252     /* Cb component */
253     put_bits(pb, 8, 2); /* index */
254     put_bits(pb, 4, 1); /* DC huffman table index */
255     put_bits(pb, 4, lossless ? 0 : 1); /* AC huffman table index */
256
257     /* Cr component */
258     put_bits(pb, 8, 3); /* index */
259     put_bits(pb, 4, 1); /* DC huffman table index */
260     put_bits(pb, 4, lossless ? 0 : 1); /* AC huffman table index */
261
262     put_bits(pb, 8, lossless ? avctx->prediction_method + 1 : 0); /* Ss (not used) */
263
264     switch (avctx->codec_id) {
265     case AV_CODEC_ID_MJPEG:  put_bits(pb, 8, 63); break; /* Se (not used) */
266     case AV_CODEC_ID_LJPEG:  put_bits(pb, 8,  0); break; /* not used */
267     default: av_assert0(0);
268     }
269
270     put_bits(pb, 8, 0); /* Ah/Al (not used) */
271
272 end:
273     if (!lossless) {
274         MpegEncContext *s = avctx->priv_data;
275         av_assert0(avctx->codec->priv_data_size == sizeof(MpegEncContext));
276
277         s->esc_pos = put_bits_count(pb) >> 3;
278         for(i=1; i<s->slice_context_count; i++)
279             s->thread_context[i]->esc_pos = 0;
280     }
281 }
282
283 void ff_mjpeg_escape_FF(PutBitContext *pb, int start)
284 {
285     int size;
286     int i, ff_count;
287     uint8_t *buf = pb->buf + start;
288     int align= (-(size_t)(buf))&3;
289     int pad = (-put_bits_count(pb))&7;
290
291     if (pad)
292         put_bits(pb, pad, (1<<pad)-1);
293
294     flush_put_bits(pb);
295     size = put_bits_count(pb) - start * 8;
296
297     av_assert1((size&7) == 0);
298     size >>= 3;
299
300     ff_count=0;
301     for(i=0; i<size && i<align; i++){
302         if(buf[i]==0xFF) ff_count++;
303     }
304     for(; i<size-15; i+=16){
305         int acc, v;
306
307         v= *(uint32_t*)(&buf[i]);
308         acc= (((v & (v>>4))&0x0F0F0F0F)+0x01010101)&0x10101010;
309         v= *(uint32_t*)(&buf[i+4]);
310         acc+=(((v & (v>>4))&0x0F0F0F0F)+0x01010101)&0x10101010;
311         v= *(uint32_t*)(&buf[i+8]);
312         acc+=(((v & (v>>4))&0x0F0F0F0F)+0x01010101)&0x10101010;
313         v= *(uint32_t*)(&buf[i+12]);
314         acc+=(((v & (v>>4))&0x0F0F0F0F)+0x01010101)&0x10101010;
315
316         acc>>=4;
317         acc+= (acc>>16);
318         acc+= (acc>>8);
319         ff_count+= acc&0xFF;
320     }
321     for(; i<size; i++){
322         if(buf[i]==0xFF) ff_count++;
323     }
324
325     if(ff_count==0) return;
326
327     flush_put_bits(pb);
328     skip_put_bytes(pb, ff_count);
329
330     for(i=size-1; ff_count; i--){
331         int v= buf[i];
332
333         if(v==0xFF){
334             buf[i+ff_count]= 0;
335             ff_count--;
336         }
337
338         buf[i+ff_count]= v;
339     }
340 }
341
342 void ff_mjpeg_encode_stuffing(MpegEncContext *s)
343 {
344     int i;
345     PutBitContext *pbc = &s->pb;
346     int mb_y = s->mb_y - !s->mb_x;
347
348     ff_mjpeg_escape_FF(pbc, s->esc_pos);
349
350     if((s->avctx->active_thread_type & FF_THREAD_SLICE) && mb_y < s->mb_height)
351         put_marker(pbc, RST0 + (mb_y&7));
352     s->esc_pos = put_bits_count(pbc) >> 3;
353
354     for(i=0; i<3; i++)
355         s->last_dc[i] = 128 << s->intra_dc_precision;
356 }
357
358 void ff_mjpeg_encode_picture_trailer(PutBitContext *pb, int header_bits)
359 {
360     av_assert1((header_bits & 7) == 0);
361
362     put_marker(pb, EOI);
363 }
364
365 void ff_mjpeg_encode_dc(PutBitContext *pb, int val,
366                         uint8_t *huff_size, uint16_t *huff_code)
367 {
368     int mant, nbits;
369
370     if (val == 0) {
371         put_bits(pb, huff_size[0], huff_code[0]);
372     } else {
373         mant = val;
374         if (val < 0) {
375             val = -val;
376             mant--;
377         }
378
379         nbits= av_log2_16bit(val) + 1;
380
381         put_bits(pb, huff_size[nbits], huff_code[nbits]);
382
383         put_sbits(pb, nbits, mant);
384     }
385 }