You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

304 lines
10KB

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