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 libavcodec/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(void) 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, 126 unsigned int pix_left, int *runs, const int *runend) 127{ 128 int mode = 0; 129 unsigned int run=0; 130 unsigned int t; 131 for(;;){ 132 t = get_vlc2(gb, ccitt_vlc[mode].table, 9, 2); 133 run += t; 134 if(t < 64){ 135 *runs++ = run; 136 if(runs >= runend){ 137 av_log(avctx, AV_LOG_ERROR, "Run overrun\n"); 138 return -1; 139 } 140 if(pix_left <= run){ 141 if(pix_left == run) 142 break; 143 av_log(avctx, AV_LOG_ERROR, "Run went out of bounds\n"); 144 return -1; 145 } 146 pix_left -= run; 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, 159 unsigned int width, int *runs, const int *runend, const int *ref) 160{ 161 int mode = 0, saved_run = 0, t; 162 int run_off = *ref++; 163 unsigned int offs=0, run= 0; 164 165 runend--; // for the last written 0 166 167 while(offs < width){ 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 } 173 if(!cmode){//pass mode 174 run_off += *ref++; 175 run = run_off - offs; 176 offs= run_off; 177 run_off += *ref++; 178 if(offs > width){ 179 av_log(avctx, AV_LOG_ERROR, "Run went out of bounds\n"); 180 return -1; 181 } 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; 198 if(runs >= runend){ 199 av_log(avctx, AV_LOG_ERROR, "Run overrun\n"); 200 return -1; 201 } 202 saved_run = 0; 203 offs += run; 204 if(offs > width || run > width){ 205 av_log(avctx, AV_LOG_ERROR, "Run went out of bounds\n"); 206 return -1; 207 } 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); 215 run_off -= *--ref; 216 offs += run; 217 if(offs > width || run > width){ 218 av_log(avctx, AV_LOG_ERROR, "Run went out of bounds\n"); 219 return -1; 220 } 221 *runs++ = run + saved_run; 222 if(runs >= runend){ 223 av_log(avctx, AV_LOG_ERROR, "Run overrun\n"); 224 return -1; 225 } 226 saved_run = 0; 227 mode = !mode; 228 } 229 //sync line pointers 230 while(run_off <= offs){ 231 run_off += *ref++; 232 run_off += *ref++; 233 } 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; 243 int run, mode = ~0, pix_left = width, run_idx = 0; 244 245 init_put_bits(&pb, dst, size*8); 246 while(pix_left > 0){ 247 run = runs[run_idx++]; 248 mode = ~mode; 249 pix_left -= run; 250 for(; run > 16; run -= 16) 251 put_sbits(&pb, 16, mode); 252 if(run) 253 put_sbits(&pb, run, mode); 254 } 255} 256 257static int find_group3_syncmarker(GetBitContext *gb, int srcsize) 258{ 259 unsigned int state = -1; 260 srcsize -= get_bits_count(gb); 261 while(srcsize-- > 0){ 262 state+= state + get_bits1(gb); 263 if((state & 0xFFF) == 1) 264 return 0; 265 } 266 return -1; 267} 268 269int ff_ccitt_unpack(AVCodecContext *avctx, 270 const uint8_t *src, int srcsize, 271 uint8_t *dst, int height, int stride, enum TiffCompr compr) 272{ 273 int j; 274 GetBitContext gb; 275 int *runs, *ref, *runend; 276 int ret; 277 int runsize= avctx->width + 2; 278 279 runs = av_malloc(runsize * sizeof(runs[0])); 280 ref = av_malloc(runsize * sizeof(ref[0])); 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++){ 286 runend = runs + runsize; 287 if(compr == TIFF_G4){ 288 ret = decode_group3_2d_line(avctx, &gb, avctx->width, runs, runend, ref); 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; 297 if(compr==TIFF_CCITT_RLE || get_bits1(&gb)) 298 ret = decode_group3_1d_line(avctx, &gb, avctx->width, runs, runend); 299 else 300 ret = decode_group3_2d_line(avctx, &gb, avctx->width, runs, runend, ref); 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} 314