一种用于拉普拉斯信源的自适应变长编码方法

1.本发明属于编码技术领域，具体涉及一种用于拉普拉斯信源的自适应变长编码方法。


背景技术：

2.拉普拉斯分布信源广泛存在于编码中，如量化残差、小波变换高频系数都服从拉普拉斯分布。
[0003][0004]
任何均值为μ的拉普拉斯信源都可以通过移位来使该信源关于0对称分布。
[0005][0006]
一般来说，对信源进行量化在传输或存储能节约相应的资源消耗。当使用量化步长qstep，得到0均值拉普拉斯信源量化系数并映射到非负整数后的概率分布如下所示：
[0007][0008]
其中
[0009][0010]
对于符合上式分布的信源，目前没有针对该概率分布的除哈夫曼编码以外的高效率变长编码方式。
[0011]
拉普拉斯信源量化系数的概率分布f(c)在c＝0时与c≠0时不属于同一解析式，故考虑其f(c)的等价分布如下
[0012][0013]
其中k为任意的正整数。规定上式中c＝-0时，c 1＝ 0。
[0014]
现有的变长编码方法包括：哈夫曼编码、一元编码、哥伦布编码、指数哥伦布编码、哥伦布rice编码。
[0015]
文献“huffman d a.a method for the construction of minimum-redundancy codes[j].proceedings of the ire,1952,40(9):1098-1101.”提出的哈夫曼编码：哈夫曼编码是任何信源的最佳变长编码，其核心思想即为给概率大的符号分配较短码字，给概率小的符号分配较长码字，从而使得平均码长最短。哈夫曼编码的构造时先挑选出概率最小的两个符号，合并这两个符号为一个集合，以这两个符号的概率和代表集合的概率，在集合
内通过符号“0”和“1”来区分这两个符号，重复该操作直至只有所有的符号合并到一个集合。图1为哈夫曼编码树构造过程。
[0016]
哈夫曼编码为最佳变长码，但是出于其编码和解码的复杂性，以及哈夫曼树需要额外的比特进行存储和传输，实际的编码应用中很少使用哈夫曼编码。故以下讨论均不涉及哈夫曼编码。
[0017]
一元编码：对于符号集合x＝{x1,x2,x3,
…
}，符号x1对应的码字为x1个“0”紧接一个1。一元编码也称为逗号编码。
[0018]
文献“golomb s.run-length encodings(corresp.)[j].ieee transactions on information theory,1966,12(3):399-401.”提出的哥伦布编码是对几何分布信源分布参数时的补充，几何分布如下所示
[0019]
f(x)＝(1-ρ)ρ
x
[0020]
哥伦布编码由前缀部分和后缀部分组成。考虑下式，哥伦布编码即是用前缀部分表示xq，使用后缀部分表示xr，其中m为参数。表1为哥伦布编码不同m对应的码字：
[0021]
x＝m
·
xq xr[0022]
表1哥伦布编码不同m对应的码字
[0023][0024][0025]
文献“taubman d,marcellin m.jpeg2000 image compression fundamentals,standards and practice:image compression fundamentals,standards and practice[m].springer science&business media,2012”提出的一种内容自适应的哥伦布编码的参数m的计算方式为：
[0026]
[0027][0028]
其中e[x]为信源的期望。
[0029]
文献“teuhola j.a compression method for clustered bit-vectors[j].information processing letters,1978,7(6):308-311.”提出的指数哥伦布编码是变长编码的一种，具有很好的结构性。指数哥伦布编码由前缀和后缀组成，前缀和后缀都依赖于指数哥伦布编码的阶数k，用于表示非负整数n的k阶指数哥伦布码可由下列步骤生成：
[0030]
1.将数字n由二进制形式写出，去掉最低的k个比特位，之后加1；
[0031]
2.计算留下的比特数，将此数减1，即是需要增加的前缀0的个数。
[0032]
3.将步骤1中去掉的最低k个比特位不会比特串尾部。
[0033]
表2 k阶指数哥伦布编码的码字
[0034][0035][0036]
文献“gallager r,van voorhis d.optimal source codes for geometrically distributed integer alphabets(corresp.)[j].ieee transactions on information theory,1975,21(2):228-230.”提出的哥伦布rice编码为哥伦布编码的简化，即当m为2的整数次幂的哥伦布编码。
[0037]
现有的变长编码方式，虽然能够完成对拉普拉斯信源的编码，但其都不能完全符合拉普拉斯信源的特点，故除哈夫曼编码外其编码效率都不如本文提出的变长编码方法。
虽然哈夫曼编码为最佳变长编码，但是其编码和解码都更加复杂，并且哈夫曼树需要更多的比特来传输和存储。


技术实现要素：

[0038]
本发明提供了一种对于拉普拉斯信源的内容自适应的变长编码方式，本发明从拉普拉斯信源分布出发，构造了一种具有自适应特性且符合信源分布特点的变长码字，优于现有的除哈夫曼编码以外的变长编码方法。
[0039]
本发明的技术方案如下：
[0040]
一种用于拉普拉斯信源的自适应变长编码方法，其特征在于，包括以下步骤：
[0041]
s1、输入待编码的拉普拉斯信源，定义符号集合为x＝{x1,x2,x3,
…
}；
[0042]
s2、估计拉普拉斯信源的参数
[0043][0044]
其中median{}表示获得集合x＝{x1,x2,x3,
…
}的中位数；
[0045]
若则将信源集合x的所有符号减去使得信源关于0对称；若则不对信源集合x操作。得到信源符号集合x
′
＝{x
′1,x
′2,x
′3,
…
}，计算信源x
′
的分布参数
[0046][0047]
其中n为符号集合x
′
的数目；
[0048]
s3、对符号集合x
′
进行量化，量化步长为qstep，得到量化系数集合c＝{c1,c2,c3,
…
}，量化的公式为：
[0049]ci
＝x
′i/qstep
[0050]
s4、将属于整数集合的量化系数集合c映射到正整数，得到编码系数集合c
′
＝{c
′1,c
′2,c
′3,
…
}
[0051][0052]
其中i＝1,2,3
…
n；
[0053]
s5、计算编码自适应参数m：
[0054][0055]
其中，表示下取整；
[0056]
s6、根据编码自适应参数m计算码字后缀参数ka和kb：
[0057][0058]
ka＝kb 1
[0059]
s7、对编码系数c
′i进行处理：当m≠2
p-1,p∈n

，且当c
′i≠0时，将c
′i加1，其中n

为正整数集合；
[0060]
s8、计算c
i,q
和c
i,r
，c
i,q
和c
i,r
：
[0061]ci,q
＝c
′i/m
[0062]ci,r
＝c
′i％m
[0063]
其中“％”表示取余运算，“/”表示除法运算，其中i＝1,2,3
…
n；
[0064]
s9、计算编码系数c
′i对应的变长编码码字，码字由前缀部分和后缀部分组成：
[0065]
前缀部分：由c
i,q
个0和一个1组成；
[0066]
后缀部分：由suffixcount个二进制位表示的数值suffix，其中
[0067][0068]
其中，当c
′i＝0时，将suffixcount减1
[0069][0070]
本发明的有益效果：
[0071]
本发明根据信源的概率分布推导得出，相对于其他变长编码方式，本发明构造的码字集合更加适合信源服从的拉普拉斯概率分布，具有更高的编码效率。
[0072]
本发明中的编码自适应参数m由信源分布参数b和量化步长qstep决定，本发明能够根据不同信源分布来得到最优变长编码码字集合，达到节省比特提高编码效率的目的。
附图说明
[0073]
图1是哈夫曼编码树构造过程；
[0074]
图2是本发明提出的对于拉普拉斯信源的自适应变长编码的构造方法框图；
[0075]
图3是本发明方法中计算码字的算法框图。
具体实施方式
[0076]
下面结合附图，详细说明本发明的技术方案。
[0077]
本发明提出的自适应变长编码方法的思路为：因为符号集合{c}的分布呈现和几何分布的相似性，直接对符号集合{c}应用哥伦布编码求得哥伦布编码中的参数m，最后当m≠2
p-1,p∈n

时合并c＝-0和c＝ 0时的码字用以表示c＝0，既可得到本发明提出的变长编码方法。本发明实际上是对哥伦布编码在拉普拉斯信源上的改进，最终得到的码字集合区别于哥伦布编码码字集合，并且，相对于哥伦布编码，新得到的变长编码码字具有更高的压缩效率。
[0078]
以下为本发明的逻辑推导：
[0079]
考虑量化系数c为下式所示，将对c的编码由对商部分cq的前缀以及对余数部分cr后缀组成。其中cq的概率密度函数由表示，cr的概率密度函数由表示。
[0080]
c＝m
·cq
cr[0081]
通过改变m的取值，使得一元编码为cq的最优码字。一元编码为cq的最优码字的充分条件即为下式：
[0082][0083]
求解上述不等式即可得到m需要满足的关系为：
[0084][0085]
将代入上式，即可得到m的计算公式：
[0086][0087]
其中b为拉普拉斯分布参数，一般根据样本进行参数估计得到的替代。根据m得到，写出符号集合c＝{-0, 0,1,2,3,
…
}的哥伦布码字集合为表1所示：
[0088]
表1哥伦布码字集合
[0089][0090][0091]
当m≠2
p-1,p∈n

时合并c＝-0和c＝ 0时的码字以表示符号c＝0，则可得到表2所示的码字集合，即本发明提出的变长编码码字集合。
[0092]
表2变长编码码字集合
1,p∈n

，且当c
′i≠0时，将c
′i加1，其中n

为正整数集合。
[0109]
计算c
i,q
和c
i,r
，c
i,q
和c
i,r
均为自然数。
[0110]ci,q
＝c
′i/m
[0111]ci,r
＝c
′i％m
[0112]
其中“％”表示取余运算，“/”表示除法运算，其中i＝1,2,3
…
n。
[0113]
h、计算编码系数c
′i对应的变长编码码字，码字由前缀部分和后缀部分组成
[0114]
前缀部分：由c
i,q
个0紧接一个1组成；
[0115]
后缀部分：由suffixcount个二进制位表示的数值suffix，其中
[0116][0117]
其中，当c
′i＝0时，将suffixcount减1。
[0118][0119]
步骤b中的参数的估计方式为：
[0120][0121]
其中median{}表示获得集合x＝{x1,x2,x3,
…
}的中位数。参数的估计方式为：
[0122][0123]
其中n为符号集合x
′
的数目。
[0124]
步骤i得到的是编码系数集合c
′
对应的码字集合，表二为不同m时编码系数集合对应的码字集合。
[0125]
下面结合仿真结论证明本发明的有效性：
[0126]
将本发明公开的自适应变长编码方法应用在suvc增强层编码时，本发明相对于0阶指数格伦布获得了平均5.13％(y-psnr)、5.11％(u-psnr)、5.14％(v-psnr)的bd-rate节省，如表5所示：
[0127]
表5
[0128][0129]