libavcodec/faxcompr.c

   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  * @file
  24  * CCITT Fax Group 3 and 4 decompression
  25  * @author Konstantin Shishkov
  26  */
  27 #include "avcodec.h"
  28 #include "get_bits.h"
  29 #include "put_bits.h"
  30 #include "faxcompr.h"
  31
  32 #define CCITT_SYMS 104
  33
  34 static const uint16_t ccitt_syms[CCITT_SYMS] = {
  35     0,    1,    2,    3,    4,    5,    6,    7,    8,    9,   10,   11,   12,
  36    13,   14,   15,   16,   17,   18,   19,   20,   21,   22,   23,   24,   25,
  37    26,   27,   28,   29,   30,   31,   32,   33,   34,   35,   36,   37,   38,
  38    39,   40,   41,   42,   43,   44,   45,   46,   47,   48,   49,   50,   51,
  39    52,   53,   54,   55,   56,   57,   58,   59,   60,   61,   62,   63,   64,
  40   128,  192,  256,  320,  384,  448,  512,  576,  640,  704,  768,  832,  896,
  41   960, 1024, 1088, 1152, 1216, 1280, 1344, 1408, 1472, 1536, 1600, 1664, 1728,
  42  1792, 1856, 1920, 1984, 2048, 2112, 2176, 2240, 2304, 2368, 2432, 2496, 2560
  43 };
  44
  45 static const uint8_t ccitt_codes_bits[2][CCITT_SYMS] =
  46 {
  47   {
  48     0x35, 0x07, 0x07, 0x08, 0x0B, 0x0C, 0x0E, 0x0F, 0x13, 0x14, 0x07, 0x08, 0x08,
  49     0x03, 0x34, 0x35, 0x2A, 0x2B, 0x27, 0x0C, 0x08, 0x17, 0x03, 0x04, 0x28, 0x2B,
  50     0x13, 0x24, 0x18, 0x02, 0x03, 0x1A, 0x1B, 0x12, 0x13, 0x14, 0x15, 0x16, 0x17,
  51     0x28, 0x29, 0x2A, 0x2B, 0x2C, 0x2D, 0x04, 0x05, 0x0A, 0x0B, 0x52, 0x53, 0x54,
  52     0x55, 0x24, 0x25, 0x58, 0x59, 0x5A, 0x5B, 0x4A, 0x4B, 0x32, 0x33, 0x34, 0x1B,
  53     0x12, 0x17, 0x37, 0x36, 0x37, 0x64, 0x65, 0x68, 0x67, 0xCC, 0xCD, 0xD2, 0xD3,
  54     0xD4, 0xD5, 0xD6, 0xD7, 0xD8, 0xD9, 0xDA, 0xDB, 0x98, 0x99, 0x9A, 0x18, 0x9B,
  55     0x08, 0x0C, 0x0D, 0x12, 0x13, 0x14, 0x15, 0x16, 0x17, 0x1C, 0x1D, 0x1E, 0x1F
  56   },
  57   {
  58     0x37, 0x02, 0x03, 0x02, 0x03, 0x03, 0x02, 0x03, 0x05, 0x04, 0x04, 0x05, 0x07,
  59     0x04, 0x07, 0x18, 0x17, 0x18, 0x08, 0x67, 0x68, 0x6C, 0x37, 0x28, 0x17, 0x18,
  60     0xCA, 0xCB, 0xCC, 0xCD, 0x68, 0x69, 0x6A, 0x6B, 0xD2, 0xD3, 0xD4, 0xD5, 0xD6,
  61     0xD7, 0x6C, 0x6D, 0xDA, 0xDB, 0x54, 0x55, 0x56, 0x57, 0x64, 0x65, 0x52, 0x53,
  62     0x24, 0x37, 0x38, 0x27, 0x28, 0x58, 0x59, 0x2B, 0x2C, 0x5A, 0x66, 0x67, 0x0F,
  63     0xC8, 0xC9, 0x5B, 0x33, 0x34, 0x35, 0x6C, 0x6D, 0x4A, 0x4B, 0x4C, 0x4D, 0x72,
  64     0x73, 0x74, 0x75, 0x76, 0x77, 0x52, 0x53, 0x54, 0x55, 0x5A, 0x5B, 0x64, 0x65,
  65     0x08, 0x0C, 0x0D, 0x12, 0x13, 0x14, 0x15, 0x16, 0x17, 0x1C, 0x1D, 0x1E, 0x1F
  66   }
  67 };
  68
  69 static const uint8_t ccitt_codes_lens[2][CCITT_SYMS] =
  70 {
  71   {
  72      8,  6,  4,  4,  4,  4,  4,  4,  5,  5,  5,  5,  6,  6,  6,  6,  6,  6,  7,  7,
  73      7,  7,  7,  7,  7,  7,  7,  7,  7,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,
  74      8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,
  75      8,  8,  8,  8,  5,  5,  6,  7,  8,  8,  8,  8,  8,  8,  9,  9,  9,  9,  9,  9,
  76      9,  9,  9,  9,  9,  9,  9,  9,  9,  6,  9, 11, 11, 11, 12, 12, 12, 12, 12, 12,
  77     12, 12, 12, 12
  78   },
  79   {
  80     10,  3,  2,  2,  3,  4,  4,  5,  6,  6,  7,  7,  7,  8,  8,  9, 10, 10, 10, 11,
  81     11, 11, 11, 11, 11, 11, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12,
  82     12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12,
  83     12, 12, 12, 12, 10, 12, 12, 12, 12, 12, 12, 13, 13, 13, 13, 13, 13, 13, 13, 13,
  84     13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 11, 11, 11, 12, 12, 12, 12, 12, 12,
  85     12, 12, 12, 12
  86   }
  87 };
  88
  89 static const uint8_t ccitt_group3_2d_bits[11] = {
  90     1, 1, 2, 2, 2, 1, 3, 3, 3, 1, 1
  91 };
  92
  93 static const uint8_t ccitt_group3_2d_lens[11] = {
  94     4, 3, 7, 6, 3, 1, 3, 6, 7, 7, 9
  95 };
  96
  97 static VLC ccitt_vlc[2], ccitt_group3_2d_vlc;
  98
  99 av_cold void ff_ccitt_unpack_init(void)
 100 {
 101     static VLC_TYPE code_table1[528][2];
 102     static VLC_TYPE code_table2[648][2];
 103     int i;
 104     static int initialized = 0;
 105
 106     if (initialized)
 107         return;
 108     ccitt_vlc[0].table = code_table1;
 109     ccitt_vlc[0].table_allocated = 528;
 110     ccitt_vlc[1].table = code_table2;
 111     ccitt_vlc[1].table_allocated = 648;
 112     for (i = 0; i < 2; i++) {
 113         ff_init_vlc_sparse(&ccitt_vlc[i], 9, CCITT_SYMS,
 114                            ccitt_codes_lens[i], 1, 1,
 115                            ccitt_codes_bits[i], 1, 1,
 116                            ccitt_syms, 2, 2,
 117                            INIT_VLC_USE_NEW_STATIC);
 118     }
 119     INIT_VLC_STATIC(&ccitt_group3_2d_vlc, 9, 11,
 120                     ccitt_group3_2d_lens, 1, 1,
 121                     ccitt_group3_2d_bits, 1, 1, 512);
 122     initialized = 1;
 123 }
 124
 125
 126 static int decode_group3_1d_line(AVCodecContext *avctx, GetBitContext *gb,
 127                                  unsigned int pix_left, int *runs,
 128                                  const int *runend)
 129 {
 130     int mode         = 0;
 131     unsigned int run = 0;
 132     unsigned int t;
 133     for (;;) {
 134         t    = get_vlc2(gb, ccitt_vlc[mode].table, 9, 2);
 135         run += t;
 136         if (t < 64) {
 137             *runs++ = run;
 138             if (runs >= runend) {
 139                 av_log(avctx, AV_LOG_ERROR, "Run overrun\n");
 140                 return AVERROR_INVALIDDATA;
 141             }
 142             if (pix_left <= run) {
 143                 if (pix_left == run)
 144                     break;
 145                 av_log(avctx, AV_LOG_ERROR, "Run went out of bounds\n");
 146                 return AVERROR_INVALIDDATA;
 147             }
 148             pix_left -= run;
 149             run       = 0;
 150             mode      = !mode;
 151         } else if ((int)t == -1) {
 152             av_log(avctx, AV_LOG_ERROR, "Incorrect code\n");
 153             return AVERROR_INVALIDDATA;
 154         }
 155     }
 156     *runs++ = 0;
 157     return 0;
 158 }
 159
 160 static int decode_group3_2d_line(AVCodecContext *avctx, GetBitContext *gb,
 161                                  unsigned int width, int *runs,
 162                                  const int *runend, const int *ref)
 163 {
 164     int mode          = 0, saved_run = 0, t;
 165     int run_off       = *ref++;
 166     unsigned int offs = 0, run = 0;
 167
 168     while (offs < width) {
 169         int cmode = get_vlc2(gb, ccitt_group3_2d_vlc.table, 9, 1);
 170         if (cmode == -1) {
 171             av_log(avctx, AV_LOG_ERROR, "Incorrect mode VLC\n");
 172             return AVERROR_INVALIDDATA;
 173         }
 174         if (!cmode) { //pass mode
 175             if (run_off < width)
 176                 run_off += *ref++;
 177             run      = run_off - offs;
 178             offs     = run_off;
 179             if (run_off < width)
 180                 run_off += *ref++;
 181             if (offs > width) {
 182                 av_log(avctx, AV_LOG_ERROR, "Run went out of bounds\n");
 183                 return AVERROR_INVALIDDATA;
 184             }
 185             saved_run += run;
 186         } else if (cmode == 1) { //horizontal mode
 187             int k;
 188             for (k = 0; k < 2; k++) {
 189                 run = 0;
 190                 for (;;) {
 191                     t = get_vlc2(gb, ccitt_vlc[mode].table, 9, 2);
 192                     if (t == -1) {
 193                         av_log(avctx, AV_LOG_ERROR, "Incorrect code\n");
 194                         return AVERROR_INVALIDDATA;
 195                     }
 196                     run += t;
 197                     if (t < 64)
 198                         break;
 199                 }
 200                 *runs++ = run + saved_run;
 201                 if (runs >= runend) {
 202                     av_log(avctx, AV_LOG_ERROR, "Run overrun\n");
 203                     return AVERROR_INVALIDDATA;
 204                 }
 205                 saved_run = 0;
 206                 offs     += run;
 207                 if (offs > width || run > width) {
 208                     av_log(avctx, AV_LOG_ERROR, "Run went out of bounds\n");
 209                     return AVERROR_INVALIDDATA;
 210                 }
 211                 mode = !mode;
 212             }
 213         } else if (cmode == 9 || cmode == 10) {
 214             avpriv_report_missing_feature(avctx, "Special modes support");
 215             return AVERROR_PATCHWELCOME;
 216         } else { //vertical mode
 217             run      = run_off - offs + (cmode - 5);
 218             run_off -= *--ref;
 219             offs    += run;
 220             if (offs > width || run > width) {
 221                 av_log(avctx, AV_LOG_ERROR, "Run went out of bounds\n");
 222                 return AVERROR_INVALIDDATA;
 223             }
 224             *runs++ = run + saved_run;
 225             if (runs >= runend) {
 226                 av_log(avctx, AV_LOG_ERROR, "Run overrun\n");
 227                 return AVERROR_INVALIDDATA;
 228             }
 229             saved_run = 0;
 230             mode      = !mode;
 231         }
 232         //sync line pointers
 233         while (offs < width && run_off <= offs) {
 234             run_off += *ref++;
 235             run_off += *ref++;
 236         }
 237     }
 238     *runs++ = saved_run;
 239     if (saved_run) {
 240         if (runs >= runend) {
 241             av_log(avctx, AV_LOG_ERROR, "Run overrun\n");
 242             return -1;
 243         }
 244         *runs++ = 0;
 245     }
 246     return 0;
 247 }
 248
 249 static void put_line(uint8_t *dst, int size, int width, const int *runs)
 250 {
 251     PutBitContext pb;
 252     int run, mode = ~0, pix_left = width, run_idx = 0;
 253
 254     init_put_bits(&pb, dst, size);
 255     while (pix_left > 0) {
 256         run       = runs[run_idx++];
 257         mode      = ~mode;
 258         pix_left -= run;
 259         for (; run > 16; run -= 16)
 260             put_sbits(&pb, 16, mode);
 261         if (run)
 262             put_sbits(&pb, run, mode);
 263     }
 264     flush_put_bits(&pb);
 265 }
 266
 267 static int find_group3_syncmarker(GetBitContext *gb, int srcsize)
 268 {
 269     unsigned int state = -1;
 270     srcsize -= get_bits_count(gb);
 271     while (srcsize-- > 0) {
 272         state += state + get_bits1(gb);
 273         if ((state & 0xFFF) == 1)
 274             return 0;
 275     }
 276     return -1;
 277 }
 278
 279 int ff_ccitt_unpack(AVCodecContext *avctx, const uint8_t *src, int srcsize,
 280                     uint8_t *dst, int height, int stride,
 281                     enum TiffCompr compr, int opts)
 282 {
 283     int j;
 284     GetBitContext gb;
 285     int *runs, *ref = NULL, *runend;
 286     int ret;
 287     int runsize = avctx->width + 2;
 288     int has_eol;
 289
 290     runs = av_malloc(runsize * sizeof(runs[0]));
 291     ref  = av_malloc(runsize * sizeof(ref[0]));
 292     if (!runs || !ref) {
 293         ret = AVERROR(ENOMEM);
 294         goto fail;
 295     }
 296     ref[0] = avctx->width;
 297     ref[1] = 0;
 298     ref[2] = 0;
 299     init_get_bits(&gb, src, srcsize * 8);
 300     has_eol = show_bits(&gb, 12) == 1 || show_bits(&gb, 16) == 1;
 301
 302     for (j = 0; j < height; j++) {
 303         runend = runs + runsize;
 304         if (compr == TIFF_G4) {
 305             ret = decode_group3_2d_line(avctx, &gb, avctx->width, runs, runend,
 306                                         ref);
 307             if (ret < 0)
 308                 goto fail;
 309         } else {
 310             int g3d1 = (compr == TIFF_G3) && !(opts & 1);
 311             if (compr != TIFF_CCITT_RLE &&
 312                 has_eol &&
 313                 find_group3_syncmarker(&gb, srcsize * 8) < 0)
 314                 break;
 315             if (compr == TIFF_CCITT_RLE || g3d1 || get_bits1(&gb))
 316                 ret = decode_group3_1d_line(avctx, &gb, avctx->width, runs,
 317                                             runend);
 318             else
 319                 ret = decode_group3_2d_line(avctx, &gb, avctx->width, runs,
 320                                             runend, ref);
 321             if (compr == TIFF_CCITT_RLE)
 322                 align_get_bits(&gb);
 323         }
 324         if (avctx->err_recognition & AV_EF_EXPLODE && ret < 0)
 325             goto fail;
 326
 327         if (ret < 0) {
 328             put_line(dst, stride, avctx->width, ref);
 329         } else {
 330             put_line(dst, stride, avctx->width, runs);
 331             FFSWAP(int *, runs, ref);
 332         }
 333         dst += stride;
 334     }
 335     ret = 0;
 336 fail:
 337     av_free(runs);
 338     av_free(ref);
 339     return ret;
 340 }