set.c
上传用户:lctgjx
上传日期:2022-06-04
资源大小:8887k
文件大小:10k
源码类别:

流媒体/Mpeg4/MP4

开发平台:

Visual C++

  1. /*****************************************************************************
  2.  * set.c: h264 encoder library
  3.  *****************************************************************************
  4.  * Copyright (C) 2005-2008 Loren Merritt <lorenm@u.washington.edu>
  5.  *
  6.  * This program is free software; you can redistribute it and/or modify
  7.  * it under the terms of the GNU General Public License as published by
  8.  * the Free Software Foundation; either version 2 of the License, or
  9.  * (at your option) any later version.
  10.  *
  11.  * This program 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
  14.  * GNU General Public License for more details.
  15.  *
  16.  * You should have received a copy of the GNU General Public License
  17.  * along with this program; if not, write to the Free Software
  18.  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02111, USA.
  19.  *****************************************************************************/
  20. #include "common.h"
  21. #define SHIFT(x,s) ((s)<0 ? (x)<<-(s) : (s)==0 ? (x) : ((x)+(1<<((s)-1)))>>(s))
  22. #define DIV(n,d) (((n) + ((d)>>1)) / (d))
  23. static const int dequant4_scale[6][3] =
  24. {
  25.     { 10, 13, 16 },
  26.     { 11, 14, 18 },
  27.     { 13, 16, 20 },
  28.     { 14, 18, 23 },
  29.     { 16, 20, 25 },
  30.     { 18, 23, 29 }
  31. };
  32. static const int quant4_scale[6][3] =
  33. {
  34.     { 13107, 8066, 5243 },
  35.     { 11916, 7490, 4660 },
  36.     { 10082, 6554, 4194 },
  37.     {  9362, 5825, 3647 },
  38.     {  8192, 5243, 3355 },
  39.     {  7282, 4559, 2893 },
  40. };
  41. static const int quant8_scan[16] =
  42. {
  43.     0,3,4,3, 3,1,5,1, 4,5,2,5, 3,1,5,1
  44. };
  45. static const int dequant8_scale[6][6] =
  46. {
  47.     { 20, 18, 32, 19, 25, 24 },
  48.     { 22, 19, 35, 21, 28, 26 },
  49.     { 26, 23, 42, 24, 33, 31 },
  50.     { 28, 25, 45, 26, 35, 33 },
  51.     { 32, 28, 51, 30, 40, 38 },
  52.     { 36, 32, 58, 34, 46, 43 },
  53. };
  54. static const int quant8_scale[6][6] =
  55. {
  56.     { 13107, 11428, 20972, 12222, 16777, 15481 },
  57.     { 11916, 10826, 19174, 11058, 14980, 14290 },
  58.     { 10082,  8943, 15978,  9675, 12710, 11985 },
  59.     {  9362,  8228, 14913,  8931, 11984, 11259 },
  60.     {  8192,  7346, 13159,  7740, 10486,  9777 },
  61.     {  7282,  6428, 11570,  6830,  9118,  8640 }
  62. };
  63. int x264_cqm_init( x264_t *h )
  64. {
  65.     int def_quant4[6][16];
  66.     int def_quant8[6][64];
  67.     int def_dequant4[6][16];
  68.     int def_dequant8[6][64];
  69.     int quant4_mf[4][6][4][4];
  70.     int quant8_mf[2][6][8][8];
  71.     int q, i, j, i_list;
  72.     int deadzone[4] = { 32 - h->param.analyse.i_luma_deadzone[1],
  73.                         32 - h->param.analyse.i_luma_deadzone[0],
  74.                         32 - 11, 32 - 21 };
  75.     int max_qp_err = -1;
  76.     int max_chroma_qp_err = -1;
  77.     for( i = 0; i < 6; i++ )
  78.     {
  79.         int size = i<4 ? 16 : 64;
  80.         for( j = (i<4 ? 0 : 4); j < i; j++ )
  81.             if( !memcmp( h->pps->scaling_list[i], h->pps->scaling_list[j], size*sizeof(uint8_t) ) )
  82.                 break;
  83.         if( j < i )
  84.         {
  85.             h->  quant4_mf[i] = h->  quant4_mf[j];
  86.             h->dequant4_mf[i] = h->dequant4_mf[j];
  87.             h->unquant4_mf[i] = h->unquant4_mf[j];
  88.         }
  89.         else
  90.         {
  91.             CHECKED_MALLOC( h->  quant4_mf[i], 52*size*sizeof(uint16_t) );
  92.             CHECKED_MALLOC( h->dequant4_mf[i],  6*size*sizeof(int) );
  93.             CHECKED_MALLOC( h->unquant4_mf[i], 52*size*sizeof(int) );
  94.         }
  95.         for( j = (i<4 ? 0 : 4); j < i; j++ )
  96.             if( deadzone[j&3] == deadzone[i&3] &&
  97.                 !memcmp( h->pps->scaling_list[i], h->pps->scaling_list[j], size*sizeof(uint8_t) ) )
  98.                 break;
  99.         if( j < i )
  100.             h->quant4_bias[i] = h->quant4_bias[j];
  101.         else
  102.             CHECKED_MALLOC( h->quant4_bias[i], 52*size*sizeof(uint16_t) );
  103.     }
  104.     for( q = 0; q < 6; q++ )
  105.     {
  106.         for( i = 0; i < 16; i++ )
  107.         {
  108.             int j = (i&1) + ((i>>2)&1);
  109.             def_dequant4[q][i] = dequant4_scale[q][j];
  110.             def_quant4[q][i]   =   quant4_scale[q][j];
  111.         }
  112.         for( i = 0; i < 64; i++ )
  113.         {
  114.             int j = quant8_scan[((i>>1)&12) | (i&3)];
  115.             def_dequant8[q][i] = dequant8_scale[q][j];
  116.             def_quant8[q][i]   =   quant8_scale[q][j];
  117.         }
  118.     }
  119.     for( q = 0; q < 6; q++ )
  120.     {
  121.         for( i_list = 0; i_list < 4; i_list++ )
  122.             for( i = 0; i < 16; i++ )
  123.             {
  124.                 h->dequant4_mf[i_list][q][0][i] = def_dequant4[q][i] * h->pps->scaling_list[i_list][i];
  125.                      quant4_mf[i_list][q][0][i] = DIV(def_quant4[q][i] * 16, h->pps->scaling_list[i_list][i]);
  126.             }
  127.         for( i_list = 0; i_list < 2; i_list++ )
  128.             for( i = 0; i < 64; i++ )
  129.             {
  130.                 h->dequant8_mf[i_list][q][0][i] = def_dequant8[q][i] * h->pps->scaling_list[4+i_list][i];
  131.                      quant8_mf[i_list][q][0][i] = DIV(def_quant8[q][i] * 16, h->pps->scaling_list[4+i_list][i]);
  132.             }
  133.     }
  134.     for( q = 0; q < 52; q++ )
  135.     {
  136.         for( i_list = 0; i_list < 4; i_list++ )
  137.             for( i = 0; i < 16; i++ )
  138.             {
  139.                 h->unquant4_mf[i_list][q][i] = (1ULL << (q/6 + 15 + 8)) / quant4_mf[i_list][q%6][0][i];
  140.                 h->  quant4_mf[i_list][q][i] = j = SHIFT(quant4_mf[i_list][q%6][0][i], q/6 - 1);
  141.                 // round to nearest, unless that would cause the deadzone to be negative
  142.                 h->quant4_bias[i_list][q][i] = X264_MIN( DIV(deadzone[i_list]<<10, j), (1<<15)/j );
  143.                 if( j > 0xffff && q > max_qp_err && (i_list == CQM_4IY || i_list == CQM_4PY) )
  144.                     max_qp_err = q;
  145.                 if( j > 0xffff && q > max_chroma_qp_err && (i_list == CQM_4IC || i_list == CQM_4PC) )
  146.                     max_chroma_qp_err = q;
  147.             }
  148.         if( h->param.analyse.b_transform_8x8 )
  149.         for( i_list = 0; i_list < 2; i_list++ )
  150.             for( i = 0; i < 64; i++ )
  151.             {
  152.                 h->unquant8_mf[i_list][q][i] = (1ULL << (q/6 + 16 + 8)) / quant8_mf[i_list][q%6][0][i];
  153.                 h->  quant8_mf[i_list][q][i] = j = SHIFT(quant8_mf[i_list][q%6][0][i], q/6);
  154.                 h->quant8_bias[i_list][q][i] = X264_MIN( DIV(deadzone[i_list]<<10, j), (1<<15)/j );
  155.                 if( j > 0xffff && q > max_qp_err )
  156.                     max_qp_err = q;
  157.             }
  158.     }
  159.     if( !h->mb.b_lossless && max_qp_err >= h->param.rc.i_qp_min )
  160.     {
  161.         x264_log( h, X264_LOG_ERROR, "Quantization overflow.  Your CQM is incompatible with QP < %d,n", max_qp_err+1 );
  162.         x264_log( h, X264_LOG_ERROR, "but min QP is set to %d.n", h->param.rc.i_qp_min );
  163.         return -1;
  164.     }
  165.     if( !h->mb.b_lossless && max_chroma_qp_err >= h->chroma_qp_table[h->param.rc.i_qp_min] )
  166.     {
  167.         x264_log( h, X264_LOG_ERROR, "Quantization overflow.  Your CQM is incompatible with QP < %d,n", max_chroma_qp_err+1 );
  168.         x264_log( h, X264_LOG_ERROR, "but min chroma QP is implied to be %d.n", h->chroma_qp_table[h->param.rc.i_qp_min] );
  169.         return -1;
  170.     }
  171.     return 0;
  172. fail:
  173.     x264_cqm_delete( h );
  174.     return -1;
  175. }
  176. void x264_cqm_delete( x264_t *h )
  177. {
  178.     int i, j;
  179.     for( i = 0; i < 6; i++ )
  180.     {
  181.         for( j = 0; j < i; j++ )
  182.             if( h->quant4_mf[i] == h->quant4_mf[j] )
  183.                 break;
  184.         if( j == i )
  185.         {
  186.             x264_free( h->  quant4_mf[i] );
  187.             x264_free( h->dequant4_mf[i] );
  188.             x264_free( h->unquant4_mf[i] );
  189.         }
  190.         for( j = 0; j < i; j++ )
  191.             if( h->quant4_bias[i] == h->quant4_bias[j] )
  192.                 break;
  193.         if( j == i )
  194.             x264_free( h->quant4_bias[i] );
  195.     }
  196. }
  197. static int x264_cqm_parse_jmlist( x264_t *h, const char *buf, const char *name,
  198.                            uint8_t *cqm, const uint8_t *jvt, int length )
  199. {
  200.     char *p;
  201.     char *nextvar;
  202.     int i;
  203.     p = strstr( buf, name );
  204.     if( !p )
  205.     {
  206.         memset( cqm, 16, length );
  207.         return 0;
  208.     }
  209.     p += strlen( name );
  210.     if( *p == 'U' || *p == 'V' )
  211.         p++;
  212.     nextvar = strstr( p, "INT" );
  213.     for( i = 0; i < length && (p = strpbrk( p, " tn," )) && (p = strpbrk( p, "0123456789" )); i++ )
  214.     {
  215.         int coef = -1;
  216.         sscanf( p, "%d", &coef );
  217.         if( i == 0 && coef == 0 )
  218.         {
  219.             memcpy( cqm, jvt, length );
  220.             return 0;
  221.         }
  222.         if( coef < 1 || coef > 255 )
  223.         {
  224.             x264_log( h, X264_LOG_ERROR, "bad coefficient in list '%s'n", name );
  225.             return -1;
  226.         }
  227.         cqm[i] = coef;
  228.     }
  229.     if( (nextvar && p > nextvar) || i != length )
  230.     {
  231.         x264_log( h, X264_LOG_ERROR, "not enough coefficients in list '%s'n", name );
  232.         return -1;
  233.     }
  234.     return 0;
  235. }
  236. int x264_cqm_parse_file( x264_t *h, const char *filename )
  237. {
  238.     char *buf, *p;
  239.     int b_error = 0;
  240.     h->param.i_cqm_preset = X264_CQM_CUSTOM;
  241.     buf = x264_slurp_file( filename );
  242.     if( !buf )
  243.     {
  244.         x264_log( h, X264_LOG_ERROR, "can't open file '%s'n", filename );
  245.         return -1;
  246.     }
  247.     while( (p = strchr( buf, '#' )) != NULL )
  248.         memset( p, ' ', strcspn( p, "n" ) );
  249.     b_error |= x264_cqm_parse_jmlist( h, buf, "INTRA4X4_LUMA",   h->param.cqm_4iy, x264_cqm_jvt4i, 16 );
  250.     b_error |= x264_cqm_parse_jmlist( h, buf, "INTRA4X4_CHROMA", h->param.cqm_4ic, x264_cqm_jvt4i, 16 );
  251.     b_error |= x264_cqm_parse_jmlist( h, buf, "INTER4X4_LUMA",   h->param.cqm_4py, x264_cqm_jvt4p, 16 );
  252.     b_error |= x264_cqm_parse_jmlist( h, buf, "INTER4X4_CHROMA", h->param.cqm_4pc, x264_cqm_jvt4p, 16 );
  253.     b_error |= x264_cqm_parse_jmlist( h, buf, "INTRA8X8_LUMA",   h->param.cqm_8iy, x264_cqm_jvt8i, 64 );
  254.     b_error |= x264_cqm_parse_jmlist( h, buf, "INTER8X8_LUMA",   h->param.cqm_8py, x264_cqm_jvt8p, 64 );
  255.     x264_free( buf );
  256.     return b_error;
  257. }