yading@10: /*
yading@10:  * CCITT Fax Group 3 and 4 decompression
yading@10:  * Copyright (c) 2008 Konstantin Shishkov
yading@10:  *
yading@10:  * This file is part of FFmpeg.
yading@10:  *
yading@10:  * FFmpeg is free software; you can redistribute it and/or
yading@10:  * modify it under the terms of the GNU Lesser General Public
yading@10:  * License as published by the Free Software Foundation; either
yading@10:  * version 2.1 of the License, or (at your option) any later version.
yading@10:  *
yading@10:  * FFmpeg is distributed in the hope that it will be useful,
yading@10:  * but WITHOUT ANY WARRANTY; without even the implied warranty of
yading@10:  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
yading@10:  * Lesser General Public License for more details.
yading@10:  *
yading@10:  * You should have received a copy of the GNU Lesser General Public
yading@10:  * License along with FFmpeg; if not, write to the Free Software
yading@10:  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
yading@10:  */
yading@10: 
yading@10: /**
yading@10:  * @file
yading@10:  * CCITT Fax Group 3 and 4 decompression
yading@10:  * @author Konstantin Shishkov
yading@10:  */
yading@10: #include "avcodec.h"
yading@10: #include "get_bits.h"
yading@10: #include "put_bits.h"
yading@10: #include "faxcompr.h"
yading@10: 
yading@10: #define CCITT_SYMS 104
yading@10: 
yading@10: static const uint16_t ccitt_syms[CCITT_SYMS] = {
yading@10:     0,    1,    2,    3,    4,    5,    6,    7,    8,    9,   10,   11,   12,
yading@10:    13,   14,   15,   16,   17,   18,   19,   20,   21,   22,   23,   24,   25,
yading@10:    26,   27,   28,   29,   30,   31,   32,   33,   34,   35,   36,   37,   38,
yading@10:    39,   40,   41,   42,   43,   44,   45,   46,   47,   48,   49,   50,   51,
yading@10:    52,   53,   54,   55,   56,   57,   58,   59,   60,   61,   62,   63,   64,
yading@10:   128,  192,  256,  320,  384,  448,  512,  576,  640,  704,  768,  832,  896,
yading@10:   960, 1024, 1088, 1152, 1216, 1280, 1344, 1408, 1472, 1536, 1600, 1664, 1728,
yading@10:  1792, 1856, 1920, 1984, 2048, 2112, 2176, 2240, 2304, 2368, 2432, 2496, 2560
yading@10: };
yading@10: 
yading@10: static const uint8_t ccitt_codes_bits[2][CCITT_SYMS] =
yading@10: {
yading@10:   {
yading@10:     0x35, 0x07, 0x07, 0x08, 0x0B, 0x0C, 0x0E, 0x0F, 0x13, 0x14, 0x07, 0x08, 0x08,
yading@10:     0x03, 0x34, 0x35, 0x2A, 0x2B, 0x27, 0x0C, 0x08, 0x17, 0x03, 0x04, 0x28, 0x2B,
yading@10:     0x13, 0x24, 0x18, 0x02, 0x03, 0x1A, 0x1B, 0x12, 0x13, 0x14, 0x15, 0x16, 0x17,
yading@10:     0x28, 0x29, 0x2A, 0x2B, 0x2C, 0x2D, 0x04, 0x05, 0x0A, 0x0B, 0x52, 0x53, 0x54,
yading@10:     0x55, 0x24, 0x25, 0x58, 0x59, 0x5A, 0x5B, 0x4A, 0x4B, 0x32, 0x33, 0x34, 0x1B,
yading@10:     0x12, 0x17, 0x37, 0x36, 0x37, 0x64, 0x65, 0x68, 0x67, 0xCC, 0xCD, 0xD2, 0xD3,
yading@10:     0xD4, 0xD5, 0xD6, 0xD7, 0xD8, 0xD9, 0xDA, 0xDB, 0x98, 0x99, 0x9A, 0x18, 0x9B,
yading@10:     0x08, 0x0C, 0x0D, 0x12, 0x13, 0x14, 0x15, 0x16, 0x17, 0x1C, 0x1D, 0x1E, 0x1F
yading@10:   },
yading@10:   {
yading@10:     0x37, 0x02, 0x03, 0x02, 0x03, 0x03, 0x02, 0x03, 0x05, 0x04, 0x04, 0x05, 0x07,
yading@10:     0x04, 0x07, 0x18, 0x17, 0x18, 0x08, 0x67, 0x68, 0x6C, 0x37, 0x28, 0x17, 0x18,
yading@10:     0xCA, 0xCB, 0xCC, 0xCD, 0x68, 0x69, 0x6A, 0x6B, 0xD2, 0xD3, 0xD4, 0xD5, 0xD6,
yading@10:     0xD7, 0x6C, 0x6D, 0xDA, 0xDB, 0x54, 0x55, 0x56, 0x57, 0x64, 0x65, 0x52, 0x53,
yading@10:     0x24, 0x37, 0x38, 0x27, 0x28, 0x58, 0x59, 0x2B, 0x2C, 0x5A, 0x66, 0x67, 0x0F,
yading@10:     0xC8, 0xC9, 0x5B, 0x33, 0x34, 0x35, 0x6C, 0x6D, 0x4A, 0x4B, 0x4C, 0x4D, 0x72,
yading@10:     0x73, 0x74, 0x75, 0x76, 0x77, 0x52, 0x53, 0x54, 0x55, 0x5A, 0x5B, 0x64, 0x65,
yading@10:     0x08, 0x0C, 0x0D, 0x12, 0x13, 0x14, 0x15, 0x16, 0x17, 0x1C, 0x1D, 0x1E, 0x1F
yading@10:   }
yading@10: };
yading@10: 
yading@10: static const uint8_t ccitt_codes_lens[2][CCITT_SYMS] =
yading@10: {
yading@10:   {
yading@10:      8,  6,  4,  4,  4,  4,  4,  4,  5,  5,  5,  5,  6,  6,  6,  6,  6,  6,  7,  7,
yading@10:      7,  7,  7,  7,  7,  7,  7,  7,  7,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,
yading@10:      8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,
yading@10:      8,  8,  8,  8,  5,  5,  6,  7,  8,  8,  8,  8,  8,  8,  9,  9,  9,  9,  9,  9,
yading@10:      9,  9,  9,  9,  9,  9,  9,  9,  9,  6,  9, 11, 11, 11, 12, 12, 12, 12, 12, 12,
yading@10:     12, 12, 12, 12
yading@10:   },
yading@10:   {
yading@10:     10,  3,  2,  2,  3,  4,  4,  5,  6,  6,  7,  7,  7,  8,  8,  9, 10, 10, 10, 11,
yading@10:     11, 11, 11, 11, 11, 11, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12,
yading@10:     12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12,
yading@10:     12, 12, 12, 12, 10, 12, 12, 12, 12, 12, 12, 13, 13, 13, 13, 13, 13, 13, 13, 13,
yading@10:     13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 11, 11, 11, 12, 12, 12, 12, 12, 12,
yading@10:     12, 12, 12, 12
yading@10:   }
yading@10: };
yading@10: 
yading@10: static const uint8_t ccitt_group3_2d_bits[11] = {
yading@10:     1, 1, 2, 2, 2, 1, 3, 3, 3, 1, 1
yading@10: };
yading@10: 
yading@10: static const uint8_t ccitt_group3_2d_lens[11] = {
yading@10:     4, 3, 7, 6, 3, 1, 3, 6, 7, 7, 9
yading@10: };
yading@10: 
yading@10: static VLC ccitt_vlc[2], ccitt_group3_2d_vlc;
yading@10: 
yading@10: av_cold void ff_ccitt_unpack_init(void)
yading@10: {
yading@10:     static VLC_TYPE code_table1[528][2];
yading@10:     static VLC_TYPE code_table2[648][2];
yading@10:     int i;
yading@10:     static int initialized = 0;
yading@10: 
yading@10:     if(initialized)
yading@10:         return;
yading@10:     ccitt_vlc[0].table = code_table1;
yading@10:     ccitt_vlc[0].table_allocated = 528;
yading@10:     ccitt_vlc[1].table = code_table2;
yading@10:     ccitt_vlc[1].table_allocated = 648;
yading@10:     for(i = 0; i < 2; i++){
yading@10:         ff_init_vlc_sparse(&ccitt_vlc[i], 9, CCITT_SYMS,
yading@10:                            ccitt_codes_lens[i], 1, 1,
yading@10:                            ccitt_codes_bits[i], 1, 1,
yading@10:                            ccitt_syms, 2, 2,
yading@10:                            INIT_VLC_USE_NEW_STATIC);
yading@10:     }
yading@10:     INIT_VLC_STATIC(&ccitt_group3_2d_vlc, 9, 11,
yading@10:                     ccitt_group3_2d_lens, 1, 1,
yading@10:                     ccitt_group3_2d_bits, 1, 1, 512);
yading@10:     initialized = 1;
yading@10: }
yading@10: 
yading@10: 
yading@10: static int decode_group3_1d_line(AVCodecContext *avctx, GetBitContext *gb,
yading@10:                                  unsigned int pix_left, int *runs, const int *runend)
yading@10: {
yading@10:     int mode = 0;
yading@10:     unsigned int run=0;
yading@10:     unsigned int t;
yading@10:     for(;;){
yading@10:         t = get_vlc2(gb, ccitt_vlc[mode].table, 9, 2);
yading@10:         run += t;
yading@10:         if(t < 64){
yading@10:             *runs++ = run;
yading@10:             if(runs >= runend){
yading@10:                 av_log(avctx, AV_LOG_ERROR, "Run overrun\n");
yading@10:                 return -1;
yading@10:             }
yading@10:             if(pix_left <= run){
yading@10:                 if(pix_left == run)
yading@10:                     break;
yading@10:                 av_log(avctx, AV_LOG_ERROR, "Run went out of bounds\n");
yading@10:                 return -1;
yading@10:             }
yading@10:             pix_left -= run;
yading@10:             run = 0;
yading@10:             mode = !mode;
yading@10:         }else if((int)t == -1){
yading@10:             av_log(avctx, AV_LOG_ERROR, "Incorrect code\n");
yading@10:             return -1;
yading@10:         }
yading@10:     }
yading@10:     *runs++ = 0;
yading@10:     return 0;
yading@10: }
yading@10: 
yading@10: static int decode_group3_2d_line(AVCodecContext *avctx, GetBitContext *gb,
yading@10:                                  unsigned int width, int *runs, const int *runend, const int *ref)
yading@10: {
yading@10:     int mode = 0, saved_run = 0, t;
yading@10:     int run_off = *ref++;
yading@10:     unsigned int offs=0, run= 0;
yading@10: 
yading@10:     while(offs < width){
yading@10:         int cmode = get_vlc2(gb, ccitt_group3_2d_vlc.table, 9, 1);
yading@10:         if(cmode == -1){
yading@10:             av_log(avctx, AV_LOG_ERROR, "Incorrect mode VLC\n");
yading@10:             return -1;
yading@10:         }
yading@10:         if(!cmode){//pass mode
yading@10:             if(run_off < width)
yading@10:                 run_off += *ref++;
yading@10:             run = run_off - offs;
yading@10:             offs= run_off;
yading@10:             if(run_off < width)
yading@10:                 run_off += *ref++;
yading@10:             if(offs > width){
yading@10:                 av_log(avctx, AV_LOG_ERROR, "Run went out of bounds\n");
yading@10:                 return -1;
yading@10:             }
yading@10:             saved_run += run;
yading@10:         }else if(cmode == 1){//horizontal mode
yading@10:             int k;
yading@10:             for(k = 0; k < 2; k++){
yading@10:                 run = 0;
yading@10:                 for(;;){
yading@10:                     t = get_vlc2(gb, ccitt_vlc[mode].table, 9, 2);
yading@10:                     if(t == -1){
yading@10:                         av_log(avctx, AV_LOG_ERROR, "Incorrect code\n");
yading@10:                         return -1;
yading@10:                     }
yading@10:                     run += t;
yading@10:                     if(t < 64)
yading@10:                         break;
yading@10:                 }
yading@10:                 *runs++ = run + saved_run;
yading@10:                 if(runs >= runend){
yading@10:                     av_log(avctx, AV_LOG_ERROR, "Run overrun\n");
yading@10:                     return -1;
yading@10:                 }
yading@10:                 saved_run = 0;
yading@10:                 offs += run;
yading@10:                 if(offs > width || run > width){
yading@10:                     av_log(avctx, AV_LOG_ERROR, "Run went out of bounds\n");
yading@10:                     return -1;
yading@10:                 }
yading@10:                 mode = !mode;
yading@10:             }
yading@10:         }else if(cmode == 9 || cmode == 10){
yading@10:             av_log(avctx, AV_LOG_ERROR, "Special modes are not supported (yet)\n");
yading@10:             return -1;
yading@10:         }else{//vertical mode
yading@10:             run = run_off - offs + (cmode - 5);
yading@10:             run_off -= *--ref;
yading@10:             offs += run;
yading@10:             if(offs > width || run > width){
yading@10:                 av_log(avctx, AV_LOG_ERROR, "Run went out of bounds\n");
yading@10:                 return -1;
yading@10:             }
yading@10:             *runs++ = run + saved_run;
yading@10:             if(runs >= runend){
yading@10:                 av_log(avctx, AV_LOG_ERROR, "Run overrun\n");
yading@10:                 return -1;
yading@10:             }
yading@10:             saved_run = 0;
yading@10:             mode = !mode;
yading@10:         }
yading@10:         //sync line pointers
yading@10:         while(offs < width && run_off <= offs){
yading@10:             run_off += *ref++;
yading@10:             run_off += *ref++;
yading@10:         }
yading@10:     }
yading@10:     *runs++ = saved_run;
yading@10:     if (saved_run) {
yading@10:         if(runs >= runend){
yading@10:             av_log(avctx, AV_LOG_ERROR, "Run overrun\n");
yading@10:             return -1;
yading@10:         }
yading@10:         *runs++ = 0;
yading@10:     }
yading@10:     return 0;
yading@10: }
yading@10: 
yading@10: static void put_line(uint8_t *dst, int size, int width, const int *runs)
yading@10: {
yading@10:     PutBitContext pb;
yading@10:     int run, mode = ~0, pix_left = width, run_idx = 0;
yading@10: 
yading@10:     init_put_bits(&pb, dst, size*8);
yading@10:     while(pix_left > 0){
yading@10:         run = runs[run_idx++];
yading@10:         mode = ~mode;
yading@10:         pix_left -= run;
yading@10:         for(; run > 16; run -= 16)
yading@10:             put_sbits(&pb, 16, mode);
yading@10:         if(run)
yading@10:             put_sbits(&pb, run, mode);
yading@10:     }
yading@10:     flush_put_bits(&pb);
yading@10: }
yading@10: 
yading@10: static int find_group3_syncmarker(GetBitContext *gb, int srcsize)
yading@10: {
yading@10:     unsigned int state = -1;
yading@10:     srcsize -= get_bits_count(gb);
yading@10:     while(srcsize-- > 0){
yading@10:         state+= state + get_bits1(gb);
yading@10:         if((state & 0xFFF) == 1)
yading@10:             return 0;
yading@10:     }
yading@10:     return -1;
yading@10: }
yading@10: 
yading@10: int ff_ccitt_unpack(AVCodecContext *avctx,
yading@10:                     const uint8_t *src, int srcsize,
yading@10:                     uint8_t *dst, int height, int stride,
yading@10:                     enum TiffCompr compr, int opts)
yading@10: {
yading@10:     int j;
yading@10:     GetBitContext gb;
yading@10:     int *runs, *ref = NULL, *runend;
yading@10:     int ret;
yading@10:     int runsize= avctx->width + 2;
yading@10:     int err = 0;
yading@10:     int has_eol;
yading@10: 
yading@10:     runs = av_malloc(runsize * sizeof(runs[0]));
yading@10:     ref  = av_malloc(runsize * sizeof(ref[0]));
yading@10:     if (!runs || ! ref) {
yading@10:         err = AVERROR(ENOMEM);
yading@10:         goto fail;
yading@10:     }
yading@10:     ref[0] = avctx->width;
yading@10:     ref[1] = 0;
yading@10:     ref[2] = 0;
yading@10:     init_get_bits(&gb, src, srcsize*8);
yading@10:     has_eol = show_bits(&gb, 12) == 1 || show_bits(&gb, 16) == 1;
yading@10: 
yading@10:     for(j = 0; j < height; j++){
yading@10:         runend = runs + runsize;
yading@10:         if(compr == TIFF_G4){
yading@10:             ret = decode_group3_2d_line(avctx, &gb, avctx->width, runs, runend, ref);
yading@10:             if(ret < 0){
yading@10:                 err = -1;
yading@10:                 goto fail;
yading@10:             }
yading@10:         }else{
yading@10:             int g3d1 = (compr == TIFF_G3) && !(opts & 1);
yading@10:             if(compr!=TIFF_CCITT_RLE && has_eol && find_group3_syncmarker(&gb, srcsize*8) < 0)
yading@10:                 break;
yading@10:             if(compr==TIFF_CCITT_RLE || g3d1 || get_bits1(&gb))
yading@10:                 ret = decode_group3_1d_line(avctx, &gb, avctx->width, runs, runend);
yading@10:             else
yading@10:                 ret = decode_group3_2d_line(avctx, &gb, avctx->width, runs, runend, ref);
yading@10:             if(compr==TIFF_CCITT_RLE)
yading@10:                 align_get_bits(&gb);
yading@10:         }
yading@10:         if(ret < 0){
yading@10:             put_line(dst, stride, avctx->width, ref);
yading@10:         }else{
yading@10:             put_line(dst, stride, avctx->width, runs);
yading@10:             FFSWAP(int*, runs, ref);
yading@10:         }
yading@10:         dst += stride;
yading@10:     }
yading@10: fail:
yading@10:     av_free(runs);
yading@10:     av_free(ref);
yading@10:     return err;
yading@10: }