ts,ps,mpeg2 decoder and analysis

哥倫布編碼的解析過程

golomb 編碼主要是針對整數(shù)進(jìn)行編碼，減少整數(shù)使用的空間位數(shù)
在golomb編碼時，有一個可以變化的參數(shù)m

本人的算法如下，很簡便的！
 

int My_get_se_golomb_31()

{

     int b = 0;

     int nValue = 0;

     b = My_get_ue_golomb_31();

     nValue = b / 2;

     if (b > 1 && b%2 ==0)

     {

         nValue = nValue * (-1);

     }

     return nValue;

}

int  My_get_ue_golomb_31()

{

     int b = 0;

     int i = 0;

     do

     {

         i ++;

         if (i <32)

              b = Bitstream_get(i,FALSE);

         else

         {

              outlog("SPS: nreserved should is zero.");

              break;

         }       

     } while(!b);

     Bitstream_get(i);

     b = (1 <<(i - 1)) -1 + Bitstream_get(i-1);    

     return b;

}

http://ludajun.blog.sohu.com/

 例如：取m = 1, 對整數(shù)x = ４進(jìn)行編碼, 算法如下
b = 2的ｍ次方
q = int((x-1)/ b)
r = x - qb - 1

由此計算出 q = 1, r = 1
二進(jìn)制編碼為q 個　１，　１　個０，　然后是r的二進(jìn)制編碼
所以編碼為：101

解碼思路：
先算出所給整數(shù)的位數(shù)n
然后從高位--〉低位找到第一個０所在的位置i，這樣就能 q = n - i, r = ｘ的第０位到第ｉ－１位所表示的數(shù)字（由低位－－〉高位）
最終解碼結(jié)果為qb + 1 + r

python 實現(xiàn)的　算法：
from math import log

m = 1

def compress(x):
    q = (x - 1) >> m
    r = x - (q << m) - 1
    result = ((((1 << q) - 1) << 1) << m) | r
    return result

def uncompress(x):
    if x == 0:
        return 1
    #計算x的位數(shù)
    n = int(log(x, 2))
    if n == 0:
        return x + 1
    for i in range(n, -1, -1):
        if (x & (1 << i)) == 0:
            q = n - i
            r = x & (1 << i - 1)
            return r + 1 + (q << m)