Commit | Line | Data |
---|---|---|
2ba45a60 DM |
1 | /* |
2 | * JPEG-LS encoder | |
3 | * Copyright (c) 2003 Michael Niedermayer | |
4 | * Copyright (c) 2006 Konstantin Shishkov | |
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-LS encoder. | |
26 | */ | |
27 | ||
28 | #include "avcodec.h" | |
29 | #include "get_bits.h" | |
30 | #include "put_bits.h" | |
31 | #include "golomb.h" | |
32 | #include "internal.h" | |
33 | #include "mathops.h" | |
34 | #include "mjpeg.h" | |
35 | #include "jpegls.h" | |
36 | ||
37 | /** | |
38 | * Encode error from regular symbol | |
39 | */ | |
40 | static inline void ls_encode_regular(JLSState *state, PutBitContext *pb, int Q, | |
41 | int err) | |
42 | { | |
43 | int k; | |
44 | int val; | |
45 | int map; | |
46 | ||
47 | for (k = 0; (state->N[Q] << k) < state->A[Q]; k++) | |
48 | ; | |
49 | ||
50 | map = !state->near && !k && (2 * state->B[Q] <= -state->N[Q]); | |
51 | ||
52 | if (err < 0) | |
53 | err += state->range; | |
54 | if (err >= (state->range + 1 >> 1)) { | |
55 | err -= state->range; | |
56 | val = 2 * FFABS(err) - 1 - map; | |
57 | } else | |
58 | val = 2 * err + map; | |
59 | ||
60 | set_ur_golomb_jpegls(pb, val, k, state->limit, state->qbpp); | |
61 | ||
62 | ff_jpegls_update_state_regular(state, Q, err); | |
63 | } | |
64 | ||
65 | /** | |
66 | * Encode error from run termination | |
67 | */ | |
68 | static inline void ls_encode_runterm(JLSState *state, PutBitContext *pb, | |
69 | int RItype, int err, int limit_add) | |
70 | { | |
71 | int k; | |
72 | int val, map; | |
73 | int Q = 365 + RItype; | |
74 | int temp; | |
75 | ||
76 | temp = state->A[Q]; | |
77 | if (RItype) | |
78 | temp += state->N[Q] >> 1; | |
79 | for (k = 0; (state->N[Q] << k) < temp; k++) | |
80 | ; | |
81 | map = 0; | |
82 | if (!k && err && (2 * state->B[Q] < state->N[Q])) | |
83 | map = 1; | |
84 | ||
85 | if (err < 0) | |
86 | val = -(2 * err) - 1 - RItype + map; | |
87 | else | |
88 | val = 2 * err - RItype - map; | |
89 | set_ur_golomb_jpegls(pb, val, k, state->limit - limit_add - 1, state->qbpp); | |
90 | ||
91 | if (err < 0) | |
92 | state->B[Q]++; | |
93 | state->A[Q] += (val + 1 - RItype) >> 1; | |
94 | ||
95 | ff_jpegls_downscale_state(state, Q); | |
96 | } | |
97 | ||
98 | /** | |
99 | * Encode run value as specified by JPEG-LS standard | |
100 | */ | |
101 | static inline void ls_encode_run(JLSState *state, PutBitContext *pb, int run, | |
102 | int comp, int trail) | |
103 | { | |
104 | while (run >= (1 << ff_log2_run[state->run_index[comp]])) { | |
105 | put_bits(pb, 1, 1); | |
106 | run -= 1 << ff_log2_run[state->run_index[comp]]; | |
107 | if (state->run_index[comp] < 31) | |
108 | state->run_index[comp]++; | |
109 | } | |
110 | /* if hit EOL, encode another full run, else encode aborted run */ | |
111 | if (!trail && run) { | |
112 | put_bits(pb, 1, 1); | |
113 | } else if (trail) { | |
114 | put_bits(pb, 1, 0); | |
115 | if (ff_log2_run[state->run_index[comp]]) | |
116 | put_bits(pb, ff_log2_run[state->run_index[comp]], run); | |
117 | } | |
118 | } | |
119 | ||
120 | /** | |
121 | * Encode one line of image | |
122 | */ | |
123 | static inline void ls_encode_line(JLSState *state, PutBitContext *pb, | |
124 | void *last, void *cur, int last2, int w, | |
125 | int stride, int comp, int bits) | |
126 | { | |
127 | int x = 0; | |
128 | int Ra, Rb, Rc, Rd; | |
129 | int D0, D1, D2; | |
130 | ||
131 | while (x < w) { | |
132 | int err, pred, sign; | |
133 | ||
134 | /* compute gradients */ | |
135 | Ra = x ? R(cur, x - stride) : R(last, x); | |
136 | Rb = R(last, x); | |
137 | Rc = x ? R(last, x - stride) : last2; | |
138 | Rd = (x >= w - stride) ? R(last, x) : R(last, x + stride); | |
139 | D0 = Rd - Rb; | |
140 | D1 = Rb - Rc; | |
141 | D2 = Rc - Ra; | |
142 | ||
143 | /* run mode */ | |
144 | if ((FFABS(D0) <= state->near) && | |
145 | (FFABS(D1) <= state->near) && | |
146 | (FFABS(D2) <= state->near)) { | |
147 | int RUNval, RItype, run; | |
148 | ||
149 | run = 0; | |
150 | RUNval = Ra; | |
151 | while (x < w && (FFABS(R(cur, x) - RUNval) <= state->near)) { | |
152 | run++; | |
153 | W(cur, x, Ra); | |
154 | x += stride; | |
155 | } | |
156 | ls_encode_run(state, pb, run, comp, x < w); | |
157 | if (x >= w) | |
158 | return; | |
159 | Rb = R(last, x); | |
160 | RItype = FFABS(Ra - Rb) <= state->near; | |
161 | pred = RItype ? Ra : Rb; | |
162 | err = R(cur, x) - pred; | |
163 | ||
164 | if (!RItype && Ra > Rb) | |
165 | err = -err; | |
166 | ||
167 | if (state->near) { | |
168 | if (err > 0) | |
169 | err = (state->near + err) / state->twonear; | |
170 | else | |
171 | err = -(state->near - err) / state->twonear; | |
172 | ||
173 | if (RItype || (Rb >= Ra)) | |
174 | Ra = av_clip(pred + err * state->twonear, 0, state->maxval); | |
175 | else | |
176 | Ra = av_clip(pred - err * state->twonear, 0, state->maxval); | |
177 | W(cur, x, Ra); | |
178 | } | |
179 | if (err < 0) | |
180 | err += state->range; | |
181 | if (err >= state->range + 1 >> 1) | |
182 | err -= state->range; | |
183 | ||
184 | ls_encode_runterm(state, pb, RItype, err, | |
185 | ff_log2_run[state->run_index[comp]]); | |
186 | ||
187 | if (state->run_index[comp] > 0) | |
188 | state->run_index[comp]--; | |
189 | } else { /* regular mode */ | |
190 | int context; | |
191 | ||
192 | context = ff_jpegls_quantize(state, D0) * 81 + | |
193 | ff_jpegls_quantize(state, D1) * 9 + | |
194 | ff_jpegls_quantize(state, D2); | |
195 | pred = mid_pred(Ra, Ra + Rb - Rc, Rb); | |
196 | ||
197 | if (context < 0) { | |
198 | context = -context; | |
199 | sign = 1; | |
200 | pred = av_clip(pred - state->C[context], 0, state->maxval); | |
201 | err = pred - R(cur, x); | |
202 | } else { | |
203 | sign = 0; | |
204 | pred = av_clip(pred + state->C[context], 0, state->maxval); | |
205 | err = R(cur, x) - pred; | |
206 | } | |
207 | ||
208 | if (state->near) { | |
209 | if (err > 0) | |
210 | err = (state->near + err) / state->twonear; | |
211 | else | |
212 | err = -(state->near - err) / state->twonear; | |
213 | if (!sign) | |
214 | Ra = av_clip(pred + err * state->twonear, 0, state->maxval); | |
215 | else | |
216 | Ra = av_clip(pred - err * state->twonear, 0, state->maxval); | |
217 | W(cur, x, Ra); | |
218 | } | |
219 | ||
220 | ls_encode_regular(state, pb, context, err); | |
221 | } | |
222 | x += stride; | |
223 | } | |
224 | } | |
225 | ||
226 | static void ls_store_lse(JLSState *state, PutBitContext *pb) | |
227 | { | |
228 | /* Test if we have default params and don't need to store LSE */ | |
229 | JLSState state2 = { 0 }; | |
230 | state2.bpp = state->bpp; | |
231 | state2.near = state->near; | |
232 | ff_jpegls_reset_coding_parameters(&state2, 1); | |
233 | if (state->T1 == state2.T1 && | |
234 | state->T2 == state2.T2 && | |
235 | state->T3 == state2.T3 && | |
236 | state->reset == state2.reset) | |
237 | return; | |
238 | /* store LSE type 1 */ | |
239 | put_marker(pb, LSE); | |
240 | put_bits(pb, 16, 13); | |
241 | put_bits(pb, 8, 1); | |
242 | put_bits(pb, 16, state->maxval); | |
243 | put_bits(pb, 16, state->T1); | |
244 | put_bits(pb, 16, state->T2); | |
245 | put_bits(pb, 16, state->T3); | |
246 | put_bits(pb, 16, state->reset); | |
247 | } | |
248 | ||
249 | static int encode_picture_ls(AVCodecContext *avctx, AVPacket *pkt, | |
250 | const AVFrame *pict, int *got_packet) | |
251 | { | |
252 | const AVFrame *const p = pict; | |
253 | const int near = avctx->prediction_method; | |
254 | PutBitContext pb, pb2; | |
255 | GetBitContext gb; | |
256 | uint8_t *buf2, *zero, *cur, *last; | |
257 | JLSState *state; | |
258 | int i, size, ret; | |
259 | int comps; | |
260 | ||
261 | if (avctx->pix_fmt == AV_PIX_FMT_GRAY8 || | |
262 | avctx->pix_fmt == AV_PIX_FMT_GRAY16) | |
263 | comps = 1; | |
264 | else | |
265 | comps = 3; | |
266 | ||
267 | if ((ret = ff_alloc_packet2(avctx, pkt, avctx->width *avctx->height * comps * 4 + | |
268 | FF_MIN_BUFFER_SIZE)) < 0) | |
269 | return ret; | |
270 | ||
271 | buf2 = av_malloc(pkt->size); | |
272 | ||
273 | init_put_bits(&pb, pkt->data, pkt->size); | |
274 | init_put_bits(&pb2, buf2, pkt->size); | |
275 | ||
276 | /* write our own JPEG header, can't use mjpeg_picture_header */ | |
277 | put_marker(&pb, SOI); | |
278 | put_marker(&pb, SOF48); | |
279 | put_bits(&pb, 16, 8 + comps * 3); // header size depends on components | |
280 | put_bits(&pb, 8, (avctx->pix_fmt == AV_PIX_FMT_GRAY16) ? 16 : 8); // bpp | |
281 | put_bits(&pb, 16, avctx->height); | |
282 | put_bits(&pb, 16, avctx->width); | |
283 | put_bits(&pb, 8, comps); // components | |
284 | for (i = 1; i <= comps; i++) { | |
285 | put_bits(&pb, 8, i); // component ID | |
286 | put_bits(&pb, 8, 0x11); // subsampling: none | |
287 | put_bits(&pb, 8, 0); // Tiq, used by JPEG-LS ext | |
288 | } | |
289 | ||
290 | put_marker(&pb, SOS); | |
291 | put_bits(&pb, 16, 6 + comps * 2); | |
292 | put_bits(&pb, 8, comps); | |
293 | for (i = 1; i <= comps; i++) { | |
294 | put_bits(&pb, 8, i); // component ID | |
295 | put_bits(&pb, 8, 0); // mapping index: none | |
296 | } | |
297 | put_bits(&pb, 8, near); | |
298 | put_bits(&pb, 8, (comps > 1) ? 1 : 0); // interleaving: 0 - plane, 1 - line | |
299 | put_bits(&pb, 8, 0); // point transform: none | |
300 | ||
301 | state = av_mallocz(sizeof(JLSState)); | |
302 | /* initialize JPEG-LS state from JPEG parameters */ | |
303 | state->near = near; | |
304 | state->bpp = (avctx->pix_fmt == AV_PIX_FMT_GRAY16) ? 16 : 8; | |
305 | ff_jpegls_reset_coding_parameters(state, 0); | |
306 | ff_jpegls_init_state(state); | |
307 | ||
308 | ls_store_lse(state, &pb); | |
309 | ||
310 | zero = av_mallocz(FFABS(p->linesize[0])); | |
311 | if (!zero) { | |
312 | av_free(state); | |
313 | return AVERROR(ENOMEM); | |
314 | } | |
315 | last = zero; | |
316 | cur = p->data[0]; | |
317 | if (avctx->pix_fmt == AV_PIX_FMT_GRAY8) { | |
318 | int t = 0; | |
319 | ||
320 | for (i = 0; i < avctx->height; i++) { | |
321 | ls_encode_line(state, &pb2, last, cur, t, avctx->width, 1, 0, 8); | |
322 | t = last[0]; | |
323 | last = cur; | |
324 | cur += p->linesize[0]; | |
325 | } | |
326 | } else if (avctx->pix_fmt == AV_PIX_FMT_GRAY16) { | |
327 | int t = 0; | |
328 | ||
329 | for (i = 0; i < avctx->height; i++) { | |
330 | ls_encode_line(state, &pb2, last, cur, t, avctx->width, 1, 0, 16); | |
331 | t = *((uint16_t *)last); | |
332 | last = cur; | |
333 | cur += p->linesize[0]; | |
334 | } | |
335 | } else if (avctx->pix_fmt == AV_PIX_FMT_RGB24) { | |
336 | int j, width; | |
337 | int Rc[3] = { 0, 0, 0 }; | |
338 | ||
339 | width = avctx->width * 3; | |
340 | for (i = 0; i < avctx->height; i++) { | |
341 | for (j = 0; j < 3; j++) { | |
342 | ls_encode_line(state, &pb2, last + j, cur + j, Rc[j], | |
343 | width, 3, j, 8); | |
344 | Rc[j] = last[j]; | |
345 | } | |
346 | last = cur; | |
347 | cur += p->linesize[0]; | |
348 | } | |
349 | } else if (avctx->pix_fmt == AV_PIX_FMT_BGR24) { | |
350 | int j, width; | |
351 | int Rc[3] = { 0, 0, 0 }; | |
352 | ||
353 | width = avctx->width * 3; | |
354 | for (i = 0; i < avctx->height; i++) { | |
355 | for (j = 2; j >= 0; j--) { | |
356 | ls_encode_line(state, &pb2, last + j, cur + j, Rc[j], | |
357 | width, 3, j, 8); | |
358 | Rc[j] = last[j]; | |
359 | } | |
360 | last = cur; | |
361 | cur += p->linesize[0]; | |
362 | } | |
363 | } | |
364 | ||
365 | av_freep(&zero); | |
366 | av_freep(&state); | |
367 | ||
368 | /* the specification says that after doing 0xff escaping unused bits in | |
369 | * the last byte must be set to 0, so just append 7 "optional" zero-bits | |
370 | * to avoid special-casing. */ | |
371 | put_bits(&pb2, 7, 0); | |
372 | size = put_bits_count(&pb2); | |
373 | flush_put_bits(&pb2); | |
374 | /* do escape coding */ | |
375 | init_get_bits(&gb, buf2, size); | |
376 | size -= 7; | |
377 | while (get_bits_count(&gb) < size) { | |
378 | int v; | |
379 | v = get_bits(&gb, 8); | |
380 | put_bits(&pb, 8, v); | |
381 | if (v == 0xFF) { | |
382 | v = get_bits(&gb, 7); | |
383 | put_bits(&pb, 8, v); | |
384 | } | |
385 | } | |
386 | avpriv_align_put_bits(&pb); | |
387 | av_free(buf2); | |
388 | ||
389 | /* End of image */ | |
390 | put_marker(&pb, EOI); | |
391 | flush_put_bits(&pb); | |
392 | ||
393 | emms_c(); | |
394 | ||
395 | pkt->size = put_bits_count(&pb) >> 3; | |
396 | pkt->flags |= AV_PKT_FLAG_KEY; | |
397 | *got_packet = 1; | |
398 | return 0; | |
399 | } | |
400 | ||
401 | static av_cold int encode_close(AVCodecContext *avctx) | |
402 | { | |
403 | av_frame_free(&avctx->coded_frame); | |
404 | return 0; | |
405 | } | |
406 | ||
407 | static av_cold int encode_init_ls(AVCodecContext *ctx) | |
408 | { | |
409 | ctx->coded_frame = av_frame_alloc(); | |
410 | if (!ctx->coded_frame) | |
411 | return AVERROR(ENOMEM); | |
412 | ||
413 | ctx->coded_frame->pict_type = AV_PICTURE_TYPE_I; | |
414 | ctx->coded_frame->key_frame = 1; | |
415 | ||
416 | if (ctx->pix_fmt != AV_PIX_FMT_GRAY8 && | |
417 | ctx->pix_fmt != AV_PIX_FMT_GRAY16 && | |
418 | ctx->pix_fmt != AV_PIX_FMT_RGB24 && | |
419 | ctx->pix_fmt != AV_PIX_FMT_BGR24) { | |
420 | av_log(ctx, AV_LOG_ERROR, | |
421 | "Only grayscale and RGB24/BGR24 images are supported\n"); | |
422 | return -1; | |
423 | } | |
424 | return 0; | |
425 | } | |
426 | ||
427 | AVCodec ff_jpegls_encoder = { | |
428 | .name = "jpegls", | |
429 | .long_name = NULL_IF_CONFIG_SMALL("JPEG-LS"), | |
430 | .type = AVMEDIA_TYPE_VIDEO, | |
431 | .id = AV_CODEC_ID_JPEGLS, | |
432 | .init = encode_init_ls, | |
433 | .close = encode_close, | |
434 | .capabilities = CODEC_CAP_FRAME_THREADS | CODEC_CAP_INTRA_ONLY, | |
435 | .encode2 = encode_picture_ls, | |
436 | .pix_fmts = (const enum AVPixelFormat[]) { | |
437 | AV_PIX_FMT_BGR24, AV_PIX_FMT_RGB24, | |
438 | AV_PIX_FMT_GRAY8, AV_PIX_FMT_GRAY16, | |
439 | AV_PIX_FMT_NONE | |
440 | }, | |
441 | }; |