CN102215406A - 一种基于dct变换的分段编码信号的快速解码方法 - Google Patents
一种基于dct变换的分段编码信号的快速解码方法 Download PDFInfo
- Publication number
- CN102215406A CN102215406A CN 201110198330 CN201110198330A CN102215406A CN 102215406 A CN102215406 A CN 102215406A CN 201110198330 CN201110198330 CN 201110198330 CN 201110198330 A CN201110198330 A CN 201110198330A CN 102215406 A CN102215406 A CN 102215406A
- Authority
- CN
- China
- Prior art keywords
- dct
- cos
- idct
- theta
- sin
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Compression, Expansion, Code Conversion, And Decoders (AREA)
- Compression Or Coding Systems Of Tv Signals (AREA)
Abstract
本发明公开了一种基于DCT变换的分段编码信号的快速解码方法,属于信号处理技术领域。本发明将长度为N/5的信号序列{a m },{b m },{c m },{d m },{e m },(m=0,1,…,N/5–1)的DCT域系数{A i },{B i },{C i },{D i },{E i },(i=0,1,…,N/5–1)转换为长度为N的原始编码信号序列{x n },(n=0,1,…,N–1)的DCT域系数{X k }(k=0,1,…,N–1),其中{X k }的计算分成{X 5 i },{X 5i+1},{X 5i+2},{X 5i+3},{X 5i+4},(i=0,1,…,N/5–1)五个部分分别进行计算,减少了DCT变换次数,从而降低了解码过程的计算复杂度。相比现有技术,本发明方法具有较低的复杂度,解码实时性更好。
Description
技术领域
本发明涉及一种解码方法,尤其涉及一种基于DCT变换的解码方法,用于对分段编码信号进行快速解码,属于信号处理技术领域。
背景技术
编解码是数字信号处理技术中极其重要的部分,编码是指将一个输入信号转换为代码,这种代码是被优化过的以利于传输或储存,解码则是编码的反向过程。编解码过程通常由编解码装置完成。通常的信号编码过程通常包括时域正变换、量化、熵编码这几个过程,解码过程包括反熵编码、反量化以及频域反变换。
离散余弦变换(Discrete Cosine Transform,DCT)是数字信号处理中一种很重要的数学工具,它与统计最佳的K-L变换非常近似,从而被广泛应用于信号编解码中。
输入序列{xn},n=0,1,...,N-1的正向DCT定义为
反向DCT(IDCT)变换定义为
其中N是序列长度并且
DCT具有许多成熟的快速算法,现有的快速算法对于长度为N=2l,l≥2,的DCT需要的复杂度为:
其中“M”代表乘法数,“A”代表加法数。
在现有基于DCT变换的编码方法中,需要发送的信号{xn}长度通常比较长,所以需要对信号进行分段编码发送,现有的文献中考虑了将信号划分为两个长度为N/2,以及三个长度为N/3的情况。但由于实际信号的长度是多变的,为了能够处理更多的情况,有时甚至会将原始将信号序列{xn}等分成五段{am},{bm},{cm},{dm},{em}的情形,即am=xm,bm=xm+N/5,cm=xm+2N/5,dm=xm+3N/5,em=xm+4N/5,m=0,1,...,N/5-1。首先将{am},{bm},{cm},{dm},{em}分别经过DCT变换得到其相应的DCT域系数{Ai},{Bi},{Ci},{Di},{Ei},i=0,1,...,N/5-1,然后对这五组系数进行量化、熵编码等处理后得到系数{A″i},{B″i},{C″i},{D″i},{E″i}传送至接收端。在解码时,首先对接收到的系数{A″i},{B″i},{C″i},{D″i},{E″i}分别进行反熵编码和反量化处理得到恢复的系数{Ai},{Bi},{Ci},{Di},{Ei},关键问题是如何通过{Ai},{Bi},{Ci},{Di},{Ei}计算出{Xk}(其中{Xk}是{xn}的长度为N的DCT的系数)?因为信号的编解码对实时性的要求相当高,所以在保证质量的情况下,要求复杂度越低越好。对于这种分段编码的信号,现有解码方法是先将输入的长度为N/5的DCT域系数{Ai},{Bi},{Ci},{Di},{Ei}分别通过IDCT反变换回时域得到原来的时域信号{am},{bm},{cm},{dm},{em},然后将这五个序列串联组合成{xn},再计算长度为N的序列{xn}的DCT的系数{Xk}。由此可以知道,传统的解码方法需要计算五个长度为N/5的IDCT和一个长度为N的DCT,具有较高的计算复杂度,从而在一定程度上影响了解码的实时性。
发明内容
本发明所要解决的技术问题在于现有基于DCT变换的解码方法对于分段编码的信号解码时计算复杂度高、实时性差的缺陷,提供一种基于DCT变换的分段编码信号的快速解码方法,该方法具有更低的计算复杂度和更好的实时性。
本发明具体采用以下技术方案:
一种基于DCT变换的分段编码信号的快速解码方法,所述分段编码信号是通过将长度为N的原始信号序列等分成五段长度为N/5的信号序列,然后分别对这五段信号序列进行DCT变换得到其相应的DCT域系数,最后对这五组DCT域系数分别进行量化、熵编码处理得到,所述快速解码方法包括以下步骤:
步骤1、对分段编码信号进行反熵编码、反量化处理,得到恢复的五组DCT域系数;
步骤2、设步骤1中得到的恢复的五组DCT域系数为{Ai},{Bi},{Ci},{Di},{Ei},(i=0,1,...,N/5-1),按照以下公式分别计算{X5i},{X5i+1},{X5i+2},{X5i+3},{X5i+4},(i=0,1,...,N/5-1):
其中
B′i=(-1)iBi,D′i=(-1)iDi,Fi=(-1)iBi+Ci,Gi=(-1)iBi-Ci,Hi=(-1)iDi+Ei,
Ji=(-1)iDi-Ei,Li=Ai+(Fi-Hi)cos(2π/5)-Ii/2,Mi=Gi+2Jicos(2π/5),
Ri=Ai-(Fi-Hi)cos(2π/5)-Fi/2,Si=2Gicos(2π/5)-Ji,
i=0,1,...,N/5-1,m=0,1,...,N/5-1,DCTN/5(.)和IDCTN/5(.)分别表示对括号中的信号序列作长度为N/5的正向和反向DCT变换;
步骤3、将步骤2得到的五个序列{X5i},{X5i+1},{X5i+2},{X5i+3},{X5i+4}中的元素依次串联组合得到序列{Xk}(k=0,1,...,N-1),{Xk}即为长度为N的原始信号序列的DCT域系数。
进一步地,{X5i},{X5i+1},{X5i+2},{X5i+3},{X5i+4},(i=0,1,...,N/5-1)采用如下方法得到:
按照下式构造并且计算中间量Ti,l,(i=1,2,...,N/5-1,l=0,1,2,3,4),
根据下式递推计算{X5i},{X5i+1},{X5i+2},{X5i+3},{X5i+4},(i=0,1,...,N/5-1),
X5i+l=Ti,l-X5i-l,(i=1,2,...,N/5-1,l=1,2,3,4)
其中初始值为Xl=T0,l/2,(l=1,2,3,4)。
相比现有技术,本发明方法的计算复杂度较低,解码的实时性更好。传统方法需要计算五个长度为N/5的IDCT和一个长度为N的DCT,其中后者可以通过已有的快速算法将其转化为五个长度为N/5的DCT,具有较高的计算复杂度,从而在一定程度上影响了解码的实时性。而本发明方法只需要计算四个长度为N/5的IDCT和四个长度为N/5的DCT。假设一个长度为N的DCT利用五个长度为N/5的DCT来实现的话,本发明方法节省了一个长度为N/5的IDCT和一个长度为N/5的DCT,因此具有更好的解码实时性。
附图说明
图1为现有方法进行分段编码的流程示意图;
图2为现有方法进行分段解码的流程示意图;
图3为本发明的快速解码方法的信号流图,其中图a为通过长度为N/5的信号{am},{bm},{cm},{dm},{em}的DCT域系数{Ai},{B′i},{Ci},{D′i},{Ei}得到中间值Ti,l,(i=1,2,...,N/5-1,l=0,1,2,3,4),图b为通过中间值Ti,l得到系数{X5i},{X5i+1},{X5i+2},{X5i+3},{X5i+4},(i=0,1,...,N/5-1)。
具体实施方式
下面结合附图对本发明的技术方案进行详细说明:
图1显示了传统的分段编码的流程,首先将要发送的信号{xn}等分成五段{am},{bm},{cm},{dm},{em},也就是am=xm,bm=xm+N/5,cm=xm+2N/5,dm=xm+3N/5,em=xm+4N/5,m=0,1,...,N/5-1,并分别对以上五个序列分别进行DCT变换得到其相应的DCT域系数{Ai},{Bi},{Ci},{Di},{Ei},i=0,1,...,N/5-1,然后对这五组系数进行量化、熵编码等处理,得到系数{A″i},{B″i},{C″i},{D″i},{E″i},将之传送至接收端或储存在介质中。
图2显示了传统方法进行分段解码的流程,首先对接收到的系数{A″i},{B″i},{C″i},{D″i},{E″i}分别进行反熵编码和反量化等处理得到恢复的系数{Ai},{Bi},{Ci},{Di},{Ei},并将{Ai},{Bi},{Ci},{Di},{Ei}分别通过IDCT反变换回时域,得到原来的时域信号{am},{bm},{cm},{dm},{em},然后将这五个序列串联组合成{xn},再计算长度为N的序列{xn}的DCT的系数{Xk}。采用传统方法时,解码过程的计算复杂度为
图3给出了用本发明方法进行N点信号解码的具体实现流图,其中输入是长度为N/5的信号{am},{bm},{cm},{dm},{em}的DCT域系数{Ai},{B′i},{Ci},{D′i},{Ei};注意:B′i=(-1)iBi,D′i=(-1)iDi。输出是长度为N的信号{xn}的DCT域系数{Xk},k=0,1,...,N-1。图中{Xk}通过{X5i},{X5i+1},{X5i+2},{X5i+3},{X5i+4},i=0,1,...,N/5-1,五个部分进行表达。
采用本发明方法进行解码时,将一个长度为N的DCT的系数{Xk}分为{X5i},{X5i+1},{X5i+2},{X5i+3},{X5i+4},i=0,1,...,N/5-1,五个部分分别进行计算。这五个部分的系数可以通过已知的DCT域系数{Ai},{Bi},{Ci},{Di},{Ei},i=0,1,...,N/5-1,计算如下,主要分为两个步骤:
步骤一:构造并且计算中间量Ti,l,(i=1,2,...,N/5-1,l=0,1,2,3,4)
我们令 则构造的中间量Ti,l可以通过四个长度为N/5的DCT和四个长度为N/5的IDCT计算如下:
其中,
am=xm,bm=xm+N/5,cm=xm+2N/5,dm=xm+3N/5,
em=xm+4N/5,b′m=bN/5-1-m,d′m=dN/5-1-m,θm=(2m+1)π/(2N),
fm=b′m+cm,gm=b′m-cm,hm=d′m+em,jm=d′m-em,
B′i=(-1)iBi,D′i=(-1)iDi,Fi=(-1)iBi+Ci,Gi=(-1)iBi-Ci,Ii=(-1)iDi+Ei,Ji=(-1)iDi-Ei,
Li=Ai+(Fi-Ii)cos(2π/5)-Ii/2,Mi=Gi+2Jicos(2π/5),
Ri=Ai-(Fi-Ii)cos(2π/5)-Fi/2,Si=2Gicos(2π/5)-Ji,
i=0,1,...,N/5-1,DCTN/5(.)和IDCTN/5(.)分别表示对括号中的信号序列作长度为N/5的正向和反向DCT变换。步骤一具体的实现过程如图3(a)所示。
步骤二:通过中间量Ti,l,(i=0,1,...,N/5-1,l=0,1,2,3,4)按照递推方式计算出{X5i},{X5i+ 1},{X5i+2},{X5i+3},{X5i+4}
具体采用如下递推的方法得到:
X5i+l=Ti,l-X5i-l,i=1,2,...,N/5-1,l=1,2,3,4,
其中初始值为Xl=T0,l/2。需要注意的是X5(i+1)-l=X5i+(5-l),即,
X5(i+1)-1=X5i+4,X5(i+1)-2=X5i+3,X5(i+1)-3=X5i+2,X5(i+1)-4=X5i+1.
我们以长度N=10进行说明,则步骤二的递推可以通过如下方式得到:
(1)计算初值Xl=T0,l/2,即
X1=T0,1/2,X2=T0,2/2,X3=T0,3/2,X4=T0,4/2.
(2)通过递推公式X5i+l=Ti,l-X5i-l,(i=1,2,...,N/5-1,l=1,2,3,4)计算如下:
X6=T1,1-X4,X7=T1,2-X3,X8=T1,3-X2,X9=T1,4-X1.
从上式可以看出,我们计算X6(即X5i+1,i=1)的值除了需要步骤一已经计算出来的T1,1外还需要X4(即X5i+4,i=0)的值,计算X9(即X5i+4,i=1)的值除了需要步骤一已经计算出来的T1,4外还需要X1(即X5i+1,i=0)的值,因此{X5i+1}和{X5i+4},i=1,2,...,N/5-1的值需要联合在一起进行计算。同样的,我们计算X7(即X5i+2,i=1)的值除了需要步骤一已经计算出来的T1,2外还需要X3(即X5i+3,i=0)的值,计算X8(即X5i+3,i=1)的值除了需要步骤一已经计算出来的T1,3外还需要X2(即X5i+2,i=0)的值,因此{X5i+2}和{X5i+3},i=1,2,...,N/5-1的值也需要联合在一起进行计算。
图3(b)给出了通过Ti,l,(i=1,2,...,N/5-1,l=1,2,3,4)的值计算{X5i},{X5i+1},{X5i+2},{X5i+3},{X5i+4}的示意图。
采用本发明方法解码时,其计算复杂度为:
下表1显示了采用本发明方法与采用传统方法解码时的计算复杂度对比。
表1
从表1可以看出,本发明的解码方法比传统方法更加有效。当序列长度N从20增加到160时,本发明方法比传统的方法节省了3%到15%的计算复杂度。
Claims (2)
1.一种基于DCT变换的分段编码信号的快速解码方法,所述分段编码信号是通过将长度为N的原始信号序列等分成五段长度为N/5的信号序列,然后分别对这五段信号序列进行DCT变换得到其相应的DCT域系数,最后对这五组DCT域系数分别进行量化、熵编码处理得到,其特征在于,所述快速解码方法包括以下步骤:
步骤1、对分段编码信号进行反熵编码、反量化处理,得到恢复的五组DCT域系数;
步骤2、设步骤1中得到的恢复的五组DCT域系数为{Ai},{Bi},{Bi},{Di},{Ei},(i=0,1,...,N/5-1),按照以下公式分别计算{X5i},{X5i+1},{X5i+2},{X5i+3},{X5i+4},(i=0,1,...,N/5-1):
其中
B′i=(-1)iBi,D′i=(-1)iDi,Fi=(-1)iBi+Ci,Gi=(-1)iBi-Ci,Hi=(-1)iDi+Ei,
Ji=(-1)iDi-Ei,Li=Ai+(Fi-Hi)cos(2π/5)-Ii/2,Mi=Gi+2Jicos(2π/5),
Ri=Ai-(Fi-Hi)cos(2π/5)-Fi/2,Si=2Gicos(2π/5)-Ji,
i=0,1,...,N/5-1,m=0,1,...,N/5-1,DCTN/5(.)和IDCTN/5(.)分别表示对括号中的信号序列作长度为N/5的正向和反向DCT变换;
步骤3、将步骤2得到的五个序列{X5i},{X5i+1},{X5i+2},{X5i+3},{X5i+4}中的元素依次串联组合得到序列{Xk},{Xk}即为长度为N的原始信号序列的DCT域系数。
2.如权利要求1所述基于DCT变换的分段编码信号的快速解码方法,其特征在于,{X5i},{X5i+1},{X5i+2},{X5i+3},{X5i+4},(i=0,1,...,N/5-1)采用如下方法得到:按照下式构造并且计算中间量Ti,l,(i=1,2,...,N/5-1,l=0,1,2,3,4),
根据下式递推计算{X5i},{X5i+1},{X5i+2},{X5i+3},{X5i+4},(i=0,1,...,N/5-1),
X5i+l=Ti,l-X5i-l,(i=1,2,...,N/5-1,l=1,2,3,4)
其中初始值为Xl=T0,l/2,(l=1,2,3,4)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110198330 CN102215406B (zh) | 2011-07-15 | 2011-07-15 | 一种基于dct变换的分段编码信号的快速解码方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110198330 CN102215406B (zh) | 2011-07-15 | 2011-07-15 | 一种基于dct变换的分段编码信号的快速解码方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102215406A true CN102215406A (zh) | 2011-10-12 |
CN102215406B CN102215406B (zh) | 2012-10-24 |
Family
ID=44746500
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110198330 Expired - Fee Related CN102215406B (zh) | 2011-07-15 | 2011-07-15 | 一种基于dct变换的分段编码信号的快速解码方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102215406B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001086505A (ja) * | 1999-09-10 | 2001-03-30 | Ddi Corp | 符号化動画像データの復号装置および記録媒体 |
US20020041630A1 (en) * | 2000-09-29 | 2002-04-11 | Yoshiharu Uetani | Compression-encoded data decoding apparatus |
CN1499846A (zh) * | 2002-11-04 | 2004-05-26 | ���ǵ�����ʽ���� | 基于能量的自适应dct/idct装置及其控制方法 |
CN101686390A (zh) * | 2008-09-24 | 2010-03-31 | 安凯(广州)软件技术有限公司 | 用于嵌入式系统的渐进模式jpeg图像分段解码方法 |
-
2011
- 2011-07-15 CN CN 201110198330 patent/CN102215406B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001086505A (ja) * | 1999-09-10 | 2001-03-30 | Ddi Corp | 符号化動画像データの復号装置および記録媒体 |
US20020041630A1 (en) * | 2000-09-29 | 2002-04-11 | Yoshiharu Uetani | Compression-encoded data decoding apparatus |
CN1499846A (zh) * | 2002-11-04 | 2004-05-26 | ���ǵ�����ʽ���� | 基于能量的自适应dct/idct装置及其控制方法 |
CN101686390A (zh) * | 2008-09-24 | 2010-03-31 | 安凯(广州)软件技术有限公司 | 用于嵌入式系统的渐进模式jpeg图像分段解码方法 |
Non-Patent Citations (2)
Title |
---|
《计算机应用》 20080630 陈琳等 多信息分段电子编码方法和系统 第28卷, * |
《计算机应用与软件》 20060630 王俊杰等 一种DCT域的音频信息隐藏方法 第23卷, 第6期 * |
Also Published As
Publication number | Publication date |
---|---|
CN102215406B (zh) | 2012-10-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107832837B (zh) | 一种基于压缩感知原理的卷积神经网络压缩方法及解压缩方法 | |
CN101742313B (zh) | 基于压缩感知技术的分布式信源编码的方法 | |
CN110648675B (zh) | 产生hoa信号的混合的空间/系数域表示的方法和设备 | |
CN102138326B (zh) | 反量化图像的方法和设备以及解码图像的方法和设备 | |
CN102298137B (zh) | 自适应数据率星载sar原始数据压缩方法 | |
JP2005189886A5 (zh) | ||
TW200843517A (en) | Code amount estimating method and apparatus, and program and storage medium therefor | |
CN101308655B (zh) | 一种音频编解码方法与装置 | |
CN102036075B (zh) | 一种图像及数字视频编码及解码方法 | |
CN102857760B (zh) | 一种无反馈的优化码率的分布式视频编解码方法及其系统 | |
CN102572427B (zh) | 基于压缩感知的多描述图像编码和解码方法 | |
WO2011137841A1 (zh) | 压缩编码方法及装置、解码解压方法及装置、以及通信系统 | |
Siddeq et al. | Applied minimized matrix size algorithm on the transformed images by DCT and DWT used for image compression | |
Ahmed et al. | Audio compression using transforms and high order entropy encoding | |
Li et al. | SPIHT algorithm combined with Huffman encoding | |
CN102572426B (zh) | 一种数据处理的方法和装置 | |
CN101960515B (zh) | 用于不同滤波器组域之间进行变换的方法和设备 | |
CN101944235B (zh) | 基于分数傅立叶变换的图像压缩方法 | |
CN102982807B (zh) | 用于对语音信号lpc系数进行多级矢量量化的方法和系统 | |
CN102215406B (zh) | 一种基于dct变换的分段编码信号的快速解码方法 | |
CN101308656A (zh) | 音频暂态信号的编解码方法 | |
CN102103859B (zh) | 一种数字音频编码、解码方法及装置 | |
CN102281443A (zh) | 基于优化的层式离散余弦变换的压缩感知图像处理方法 | |
CN103533351B (zh) | 一种多量化表的图像压缩方法 | |
CN102163976B (zh) | 一种基于gdft-ii变换的快速解码方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121024 Termination date: 20140715 |
|
EXPY | Termination of patent right or utility model |