indent
[libav.git] / libavcodec / faxcompr.c
CommitLineData
8d48db31
MN
1/*
2 * CCITT Fax Group 3 and 4 decompression
3 * Copyright (c) 2008 Konstantin Shishkov
4 *
5 * This file is part of FFmpeg.
6 *
7 * FFmpeg is free software; you can redistribute it and/or
8 * modify it under the terms of the GNU Lesser General Public
9 * License as published by the Free Software Foundation; either
10 * version 2.1 of the License, or (at your option) any later version.
11 *
12 * FFmpeg is distributed in the hope that it will be useful,
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15 * Lesser General Public License for more details.
16 *
17 * You should have received a copy of the GNU Lesser General Public
18 * License along with FFmpeg; if not, write to the Free Software
19 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
20 */
21
22/**
23 * CCITT Fax Group 3 and 4 decompression
24 * @file faxcompr.c
25 * @author Konstantin Shishkov
26 */
27#include "avcodec.h"
28#include "bitstream.h"
29#include "faxcompr.h"
30
31#define CCITT_SYMS 104
32
33static const uint16_t ccitt_syms[CCITT_SYMS] = {
34 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12,
35 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25,
36 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38,
37 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51,
38 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64,
39 128, 192, 256, 320, 384, 448, 512, 576, 640, 704, 768, 832, 896,
40 960, 1024, 1088, 1152, 1216, 1280, 1344, 1408, 1472, 1536, 1600, 1664, 1728,
41 1792, 1856, 1920, 1984, 2048, 2112, 2176, 2240, 2304, 2368, 2432, 2496, 2560
42};
43
44static const uint8_t ccitt_codes_bits[2][CCITT_SYMS] =
45{
46 {
47 0x35, 0x07, 0x07, 0x08, 0x0B, 0x0C, 0x0E, 0x0F, 0x13, 0x14, 0x07, 0x08, 0x08,
48 0x03, 0x34, 0x35, 0x2A, 0x2B, 0x27, 0x0C, 0x08, 0x17, 0x03, 0x04, 0x28, 0x2B,
49 0x13, 0x24, 0x18, 0x02, 0x03, 0x1A, 0x1B, 0x12, 0x13, 0x14, 0x15, 0x16, 0x17,
50 0x28, 0x29, 0x2A, 0x2B, 0x2C, 0x2D, 0x04, 0x05, 0x0A, 0x0B, 0x52, 0x53, 0x54,
51 0x55, 0x24, 0x25, 0x58, 0x59, 0x5A, 0x5B, 0x4A, 0x4B, 0x32, 0x33, 0x34, 0x1B,
52 0x12, 0x17, 0x37, 0x36, 0x37, 0x64, 0x65, 0x68, 0x67, 0xCC, 0xCD, 0xD2, 0xD3,
53 0xD4, 0xD5, 0xD6, 0xD7, 0xD8, 0xD9, 0xDA, 0xDB, 0x98, 0x99, 0x9A, 0x18, 0x9B,
54 0x08, 0x0C, 0x0D, 0x12, 0x13, 0x14, 0x15, 0x16, 0x17, 0x1C, 0x1D, 0x1E, 0x1F
55 },
56 {
57 0x37, 0x02, 0x03, 0x02, 0x03, 0x03, 0x02, 0x03, 0x05, 0x04, 0x04, 0x05, 0x07,
58 0x04, 0x07, 0x18, 0x17, 0x18, 0x08, 0x67, 0x68, 0x6C, 0x37, 0x28, 0x17, 0x18,
59 0xCA, 0xCB, 0xCC, 0xCD, 0x68, 0x69, 0x6A, 0x6B, 0xD2, 0xD3, 0xD4, 0xD5, 0xD6,
60 0xD7, 0x6C, 0x6D, 0xDA, 0xDB, 0x54, 0x55, 0x56, 0x57, 0x64, 0x65, 0x52, 0x53,
61 0x24, 0x37, 0x38, 0x27, 0x28, 0x58, 0x59, 0x2B, 0x2C, 0x5A, 0x66, 0x67, 0x0F,
62 0xC8, 0xC9, 0x5B, 0x33, 0x34, 0x35, 0x6C, 0x6D, 0x4A, 0x4B, 0x4C, 0x4D, 0x72,
63 0x73, 0x74, 0x75, 0x76, 0x77, 0x52, 0x53, 0x54, 0x55, 0x5A, 0x5B, 0x64, 0x65,
64 0x08, 0x0C, 0x0D, 0x12, 0x13, 0x14, 0x15, 0x16, 0x17, 0x1C, 0x1D, 0x1E, 0x1F
65 }
66};
67
68static const uint8_t ccitt_codes_lens[2][CCITT_SYMS] =
69{
70 {
71 8, 6, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 7, 7,
72 7, 7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8,
73 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8,
74 8, 8, 8, 8, 5, 5, 6, 7, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9,
75 9, 9, 9, 9, 9, 9, 9, 9, 9, 6, 9, 11, 11, 11, 12, 12, 12, 12, 12, 12,
76 12, 12, 12, 12
77 },
78 {
79 10, 3, 2, 2, 3, 4, 4, 5, 6, 6, 7, 7, 7, 8, 8, 9, 10, 10, 10, 11,
80 11, 11, 11, 11, 11, 11, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12,
81 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12,
82 12, 12, 12, 12, 10, 12, 12, 12, 12, 12, 12, 13, 13, 13, 13, 13, 13, 13, 13, 13,
83 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 11, 11, 11, 12, 12, 12, 12, 12, 12,
84 12, 12, 12, 12
85 }
86};
87
88static const uint8_t ccitt_group3_2d_bits[11] = {
89 1, 1, 2, 2, 2, 1, 3, 3, 3, 1, 1
90};
91
92static const uint8_t ccitt_group3_2d_lens[11] = {
93 4, 3, 7, 6, 3, 1, 3, 6, 7, 7, 9
94};
95
96static VLC ccitt_vlc[2], ccitt_group3_2d_vlc;
97
98av_cold void ff_ccitt_unpack_init()
99{
100 static VLC_TYPE code_table1[528][2];
101 static VLC_TYPE code_table2[648][2];
102 int i;
103 static int initialized = 0;
104
105 if(initialized)
106 return;
107 ccitt_vlc[0].table = code_table1;
108 ccitt_vlc[0].table_allocated = 528;
109 ccitt_vlc[1].table = code_table2;
110 ccitt_vlc[1].table_allocated = 648;
111 for(i = 0; i < 2; i++){
112 init_vlc_sparse(&ccitt_vlc[i], 9, CCITT_SYMS,
113 ccitt_codes_lens[i], 1, 1,
114 ccitt_codes_bits[i], 1, 1,
115 ccitt_syms, 2, 2,
116 INIT_VLC_USE_NEW_STATIC);
117 }
118 INIT_VLC_STATIC(&ccitt_group3_2d_vlc, 9, 11,
119 ccitt_group3_2d_lens, 1, 1,
120 ccitt_group3_2d_bits, 1, 1, 512);
121 initialized = 1;
122}
123
124
125static int decode_group3_1d_line(AVCodecContext *avctx, GetBitContext *gb,
11c5999e 126 unsigned int pix_left, int *runs, const int *runend)
8d48db31 127{
11c5999e
MN
128 int mode = 0;
129 unsigned int run=0;
8d48db31
MN
130 unsigned int t;
131 for(;;){
132 t = get_vlc2(gb, ccitt_vlc[mode].table, 9, 2);
133 run += t;
134 if(t < 64){
8d48db31 135 *runs++ = run;
1d798507
MN
136 if(runs >= runend){
137 av_log(avctx, AV_LOG_ERROR, "Run overrun\n");
138 return -1;
139 }
11c5999e
MN
140 if(pix_left <= run){
141 if(pix_left == run)
8d48db31 142 break;
8d48db31
MN
143 av_log(avctx, AV_LOG_ERROR, "Run went out of bounds\n");
144 return -1;
145 }
11c5999e 146 pix_left -= run;
8d48db31
MN
147 run = 0;
148 mode = !mode;
149 }else if((int)t == -1){
150 av_log(avctx, AV_LOG_ERROR, "Incorrect code\n");
151 return -1;
152 }
153 }
154 *runs++ = 0;
155 return 0;
156}
157
158static int decode_group3_2d_line(AVCodecContext *avctx, GetBitContext *gb,
26830975 159 unsigned int width, int *runs, const int *runend, const int *ref)
8d48db31 160{
26830975 161 int mode = 0, saved_run = 0, t;
8d48db31 162 int run_off = *ref++;
26830975 163 unsigned int offs=0, run= 0;
8d48db31 164
1d798507
MN
165 runend--; // for the last written 0
166
43ba8f35 167 while(offs < width){
8d48db31
MN
168 int cmode = get_vlc2(gb, ccitt_group3_2d_vlc.table, 9, 1);
169 if(cmode == -1){
170 av_log(avctx, AV_LOG_ERROR, "Incorrect mode VLC\n");
171 return -1;
172 }
8d48db31
MN
173 if(!cmode){//pass mode
174 run_off += *ref++;
175 run = run_off - offs;
2a522579 176 offs= run_off;
8d48db31 177 run_off += *ref++;
43ba8f35 178 if(offs > width){
8d48db31
MN
179 av_log(avctx, AV_LOG_ERROR, "Run went out of bounds\n");
180 return -1;
181 }
8d48db31
MN
182 saved_run += run;
183 }else if(cmode == 1){//horizontal mode
184 int k;
185 for(k = 0; k < 2; k++){
186 run = 0;
187 for(;;){
188 t = get_vlc2(gb, ccitt_vlc[mode].table, 9, 2);
189 if(t == -1){
190 av_log(avctx, AV_LOG_ERROR, "Incorrect code\n");
191 return -1;
192 }
193 run += t;
194 if(t < 64)
195 break;
196 }
197 *runs++ = run + saved_run;
1d798507
MN
198 if(runs >= runend){
199 av_log(avctx, AV_LOG_ERROR, "Run overrun\n");
200 return -1;
201 }
8d48db31 202 saved_run = 0;
43ba8f35 203 offs += run;
26830975 204 if(offs > width || run > width){
8d48db31
MN
205 av_log(avctx, AV_LOG_ERROR, "Run went out of bounds\n");
206 return -1;
207 }
8d48db31
MN
208 mode = !mode;
209 }
210 }else if(cmode == 9 || cmode == 10){
211 av_log(avctx, AV_LOG_ERROR, "Special modes are not supported (yet)\n");
212 return -1;
213 }else{//vertical mode
214 run = run_off - offs + (cmode - 5);
a56f82ea 215 run_off -= *--ref;
43ba8f35 216 offs += run;
eb0094b2 217 if(offs > width || run > width){
8d48db31
MN
218 av_log(avctx, AV_LOG_ERROR, "Run went out of bounds\n");
219 return -1;
220 }
8d48db31 221 *runs++ = run + saved_run;
1d798507
MN
222 if(runs >= runend){
223 av_log(avctx, AV_LOG_ERROR, "Run overrun\n");
224 return -1;
225 }
8d48db31
MN
226 saved_run = 0;
227 mode = !mode;
228 }
6f5c00eb
MN
229 //sync line pointers
230 while(run_off <= offs){
231 run_off += *ref++;
232 run_off += *ref++;
233 }
8d48db31
MN
234 }
235 *runs++ = saved_run;
236 *runs++ = 0;
237 return 0;
238}
239
240static void put_line(uint8_t *dst, int size, int width, const int *runs)
241{
242 PutBitContext pb;
bc0f7a9c 243 int run, mode = ~0, pix_left = width, run_idx = 0;
8d48db31
MN
244
245 init_put_bits(&pb, dst, size*8);
246 while(pix_left > 0){
247 run = runs[run_idx++];
bc0f7a9c 248 mode = ~mode;
8d48db31
MN
249 pix_left -= run;
250 for(; run > 16; run -= 16)
bc0f7a9c 251 put_sbits(&pb, 16, mode);
68d48cc5 252 if(run)
679f3fa9 253 put_sbits(&pb, run, mode);
8d48db31
MN
254 }
255}
256
257static int find_group3_syncmarker(GetBitContext *gb, int srcsize)
258{
eb12536a
MN
259 unsigned int state = -1;
260 srcsize -= get_bits_count(gb);
261 while(srcsize-- > 0){
262 state+= state + get_bits1(gb);
ec5bdf49 263 if((state & 0xFFF) == 1)
eb12536a 264 return 0;
8d48db31 265 }
eb12536a 266 return -1;
8d48db31
MN
267}
268
66909bc4 269int ff_ccitt_unpack(AVCodecContext *avctx,
8d48db31 270 const uint8_t *src, int srcsize,
66909bc4 271 uint8_t *dst, int height, int stride, enum TiffCompr compr)
8d48db31
MN
272{
273 int j;
274 GetBitContext gb;
1d798507 275 int *runs, *ref, *runend;
8d48db31 276 int ret;
6fe68c7b 277 int runsize= avctx->width + 2;
8d48db31 278
6fe68c7b
MN
279 runs = av_malloc(runsize * sizeof(runs[0]));
280 ref = av_malloc(runsize * sizeof(ref[0]));
8d48db31
MN
281 ref[0] = avctx->width;
282 ref[1] = 0;
283 ref[2] = 0;
284 init_get_bits(&gb, src, srcsize*8);
285 for(j = 0; j < height; j++){
6fe68c7b 286 runend = runs + runsize;
66909bc4 287 if(compr == TIFF_G4){
1d798507 288 ret = decode_group3_2d_line(avctx, &gb, avctx->width, runs, runend, ref);
8d48db31
MN
289 if(ret < 0){
290 av_free(runs);
291 av_free(ref);
292 return -1;
293 }
294 }else{
295 if(find_group3_syncmarker(&gb, srcsize*8) < 0)
296 break;
66909bc4 297 if(compr==TIFF_CCITT_RLE || get_bits1(&gb))
1d798507 298 ret = decode_group3_1d_line(avctx, &gb, avctx->width, runs, runend);
8d48db31 299 else
1d798507 300 ret = decode_group3_2d_line(avctx, &gb, avctx->width, runs, runend, ref);
8d48db31
MN
301 }
302 if(ret < 0){
303 put_line(dst, stride, avctx->width, ref);
304 }else{
305 put_line(dst, stride, avctx->width, runs);
306 FFSWAP(int*, runs, ref);
307 }
308 dst += stride;
309 }
310 av_free(runs);
311 av_free(ref);
312 return 0;
313}