1. ffmpeg
  2. Commit 8d48db314eeb088e264dedf088b90e948b151b90
Browse code

CCITT Fax Group 3 and 4 decompression by kostya. Will be activated after a little more reviewing and fixing by me. Commited now, so as to not lose history of that.

Originally committed as revision 16323 to svn://svn.ffmpeg.org/ffmpeg/trunk

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