mpeglayer1.cpp
上传用户:sun1608
上传日期:2007-02-02
资源大小:6116k
文件大小:3k
源码类别:

流媒体/Mpeg4/MP4

开发平台:

Visual C++

  1. /* MPEG/WAVE Sound library
  2.    (C) 1997 by Jung woo-jae */
  3. // Mpeglayer1.cc
  4. // It's for MPEG Layer 1
  5. #ifdef HAVE_CONFIG_H
  6. #include "config.h"
  7. #endif
  8. #include "MPEGaudio.h"
  9. #if defined(_WIN32) && defined(_MSC_VER)
  10. // disable warnings about double to float conversions
  11. #pragma warning(disable: 4244 4305)
  12. #endif
  13. // Tables for layer 1
  14. static const REAL factortable[15] = 
  15. {
  16.   0.0, 
  17.   (1.0/2.0)    * (4.0/3.0),         (1.0/4.0)     * (8.0/7.0),
  18.   (1.0/8.0)    * (16.0/15.0),       (1.0/16.0)    * (32.0/31.0),
  19.   (1.0/32.0)   * (64.0/63.0),       (1.0/64.0)    * (128.0/127.0),
  20.   (1.0/128.0)  * (256.0/255.0),     (1.0/256.0)   * (512.0/511.0),
  21.   (1.0/512.0)  * (1024.0/1023.0),   (1.0/1024.0)  * (2048.0/2047.0),
  22.   (1.0/2048.0) * (4096.0/4095.0),   (1.0/4096.0)  * (8192.0/8191.0),
  23.   (1.0/8192.0) * (16384.0/16383.0), (1.0/16384.0) * (32768.0/32767.0)
  24. };
  25. static const REAL offsettable[15] = 
  26. {
  27.   0.0, 
  28.   ((1.0/2.0)-1.0)    * (4.0/3.0),         ((1.0/4.0)-1.0)     * (8.0/7.0), 
  29.   ((1.0/8.0)-1.0)    * (16.0/15.0),       ((1.0/16.0)-1.0)    * (32.0/31.0),
  30.   ((1.0/32.0)-1.0)   * (64.0/63.0),       ((1.0/64.0)-1.0)    * (128.0/127.0),
  31.   ((1.0/128.0)-1.0)  * (256.0/255.0),     ((1.0/256.0)-1.0)   * (512.0/511.0),
  32.   ((1.0/512.0)-1.0)  * (1024.0/1023.0),   ((1.0/1024.0)-1.0)  * (2048.0/2047.0),
  33.   ((1.0/2048.0)-1.0) * (4096.0/4095.0),   ((1.0/4096.0)-1.0)  * (8192.0/8191.0),
  34.   ((1.0/8192.0)-1.0) * (16384.0/16383.0), ((1.0/16384.0)-1.0) * (32768.0/32767.0)
  35. };
  36. // Mpeg layer 1
  37. void MPEGaudio::extractlayer1(void)
  38. {
  39.   REAL fraction[MAXCHANNEL][MAXSUBBAND];
  40.   REAL scalefactor[MAXCHANNEL][MAXSUBBAND];
  41.   int bitalloc[MAXCHANNEL][MAXSUBBAND],
  42.       sample[MAXCHANNEL][MAXSUBBAND];
  43.   register int i,j;
  44.   int s=stereobound,l;
  45. // Bitalloc
  46.   for(i=0;i<s;i++)
  47.   {
  48.     bitalloc[LS][i]=getbits(4);
  49.     bitalloc[RS][i]=getbits(4);
  50.   }
  51.   for(;i<MAXSUBBAND;i++)
  52.     bitalloc[LS][i]=
  53.     bitalloc[RS][i]=getbits(4);
  54. // Scale index
  55.   if(inputstereo)
  56.     for(i=0;i<MAXSUBBAND;i++)
  57.     {
  58.       if(bitalloc[LS][i])scalefactor[LS][i]=scalefactorstable[getbits(6)];
  59.       if(bitalloc[RS][i])scalefactor[RS][i]=scalefactorstable[getbits(6)];
  60.     }
  61.   else
  62.     for(i=0;i<MAXSUBBAND;i++)
  63.       if(bitalloc[LS][i])scalefactor[LS][i]=scalefactorstable[getbits(6)];
  64.   for(l=0;l<SCALEBLOCK;l++)
  65.   {
  66.     // Sample
  67.     for(i=0;i<s;i++)
  68.     {
  69.       if((j=bitalloc[LS][i]))sample[LS][i]=getbits(j+1);
  70.       if((j=bitalloc[RS][i]))sample[RS][i]=getbits(j+1);
  71.     }
  72.     for(;i<MAXSUBBAND;i++)
  73.       if((j=bitalloc[LS][i]))sample[LS][i]=sample[RS][i]=getbits(j+1);
  74.     // Fraction
  75.     if(outputstereo)
  76.       for(i=0;i<MAXSUBBAND;i++)
  77.       {
  78. if((j=bitalloc[LS][i]))
  79.   fraction[LS][i]=(REAL(sample[LS][i])*factortable[j]+offsettable[j])
  80.   *scalefactor[LS][i];
  81. else fraction[LS][i]=0.0;
  82. if((j=bitalloc[RS][i]))
  83.   fraction[RS][i]=(REAL(sample[RS][i])*factortable[j]+offsettable[j])
  84.   *scalefactor[RS][i];
  85. else fraction[RS][i]=0.0;
  86.       }
  87.     else
  88.       for(i=0;i<MAXSUBBAND;i++)
  89. if((j=bitalloc[LS][i]))
  90.   fraction[LS][i]=(REAL(sample[LS][i])*factortable[j]+offsettable[j])
  91.   *scalefactor[LS][i];
  92. else fraction[LS][i]=0.0;
  93.     subbandsynthesis(fraction[LS],fraction[RS]);
  94.   }
  95. }