CN102281437A - 用于图像压缩的提升结构二维离散小波变换交织扫描方法 - Google Patents

用于图像压缩的提升结构二维离散小波变换交织扫描方法 Download PDF

Info

Publication number
CN102281437A
CN102281437A CN201110148246XA CN201110148246A CN102281437A CN 102281437 A CN102281437 A CN 102281437A CN 201110148246X A CN201110148246X A CN 201110148246XA CN 201110148246 A CN201110148246 A CN 201110148246A CN 102281437 A CN102281437 A CN 102281437A
Authority
CN
China
Prior art keywords
row
dimensional discrete
discrete wavelet
formula
read
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
Application number
CN201110148246XA
Other languages
English (en)
Other versions
CN102281437B (zh
Inventor
刘文松
王健
朱恩
徐龙涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN 201110148246 priority Critical patent/CN102281437B/zh
Publication of CN102281437A publication Critical patent/CN102281437A/zh
Application granted granted Critical
Publication of CN102281437B publication Critical patent/CN102281437B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Compression Or Coding Systems Of Tv Signals (AREA)
  • Compression Of Band Width Or Redundancy In Fax (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

该发明涉及用于图像压缩的提升结构二维离散小波变换交织扫描方法，特别是应用于符合联合图像专家组（JPEG）2000图像压缩标准的、具有高效率低存储特点的、基于5/3提升结构的二维离散小波变换交织扫描方法。该发明突破了传统的光栅扫描方法所采用的先行后列的数据读入方式，按照行列交织扫描方法读入图像数据，使得5/3提升结构离散小波变换过程中处理，水平方向一维离散小波变换输出的低频系数无需存储，仅需缓存N个高频系数，便可同时进行垂直方向一维离散小波变换，同时输出四个频率子带的小波系数。相对于传统最优的设计方案，该发明将内部存储器的消耗降低了50%。

Description

用于图像压缩的提升结构二维离散小波变换交织扫描方法
技术领域
本发明涉及用于图像压缩的提升结构二维离散小波变换交织扫描方法，特别是应用于符合联合图像专家组(JPEG)2000图像压缩标准的、低内存消耗的，基于5/3提升结构的二维离散小波变换交织扫描方法。 
背景技术
目前国际最新的图像压缩标准JPEG2000（参考文献[1]）具备帧内编码、压缩效果好、支持感兴趣区域编码和渐进式传输等品质，因而在高端数字图像应用中得到了推广。在图像压缩过程中，首先采用提升结构二维离散小波变换（DWT）原始图像，之后采用位平面编码器（BPC)处理小波系数子带生成上下文和码值(CX/D)，再送入算术编码器(AC)进行压缩编码。现有的DWT已经能够满足图像编码器的实时处理要求。但随着各类高清数字图像的应用普及，所消耗的内部存储器急剧增高。如何降低DWT在计算过程中所消耗的内部存储器已经成了该领域研究热点之一。 
提升结构离散小波变换由Mallat提出，极大减小了离散小波变换的计算复杂度。JPEG2000采用的提升结构二维离散小波变换有两种，一种是9/7提升结构离散小波变换，另一种是5/3提升结构离散小波变换。无损图像压缩采用基于5/3提升结构离散小波变换，即按照光栅扫描方法，自左向右自上而下的读入图像数据进行水平方向的一维离散小波变换，即用公式1计算高频系数，用公式2计算低频系数，公式中X(m, n)表示一副 
Figure 201110148246X100002DEST_PATH_IMAGE001
行列图像中第m行第n列像素，
Figure 894614DEST_PATH_IMAGE002
表示水平方向一维离散小波变换后生成的第m行第n列个高频系数，
Figure 201110148246X100002DEST_PATH_IMAGE003
表示第m行第n列个高频系数；用公式3和公式4对高频系数
Figure 969886DEST_PATH_IMAGE002
进行垂直方向一维离散小波变换，产生甚高频系数
Figure 852392DEST_PATH_IMAGE004
和高低频系数；用公式5和公式6对低频系数
Figure 332439DEST_PATH_IMAGE003
进行垂直方向一维离散小波变换，生成低高频系数
Figure 274987DEST_PATH_IMAGE006
、甚低频系数
Figure 201110148246X100002DEST_PATH_IMAGE007
。 
                 公式1 
           
Figure 201110148246X100002DEST_PATH_IMAGE009
       公式2
               公式3
           
Figure 201110148246X100002DEST_PATH_IMAGE011
   公式4
Figure 855376DEST_PATH_IMAGE012
      公式5
           
Figure 201110148246X100002DEST_PATH_IMAGE013
     公式6
提升结构二维离散小波变换的实质是采用光栅扫描方法读入图像数据，先进行水平方向一维小波离散变换，后进行垂直方向一维离散小波变换。在进行垂直方向一维离散小波变换之前，需要存储水平方向一维离散小波变换所生成的所有H(m, n)和L(m, n)。对于N行N列的图像，其完成5/3提升结构二维离散小波变换需要内部存储器Trans_Mem为:
Figure 155776DEST_PATH_IMAGE014
                            公式7
根据公式7，随着高清高分辨率数字图像的广泛应用，5/3提升结构二维离散小波变换所需的内部存储器会随之急剧增高。如何降低5/3提升结构离散小波变换的内部存储器成为该领域的研究热点之一。
具有代表性的研究成果有：参考文献[2]提出了一种对称折叠流水线方法；参考文献[3]提出了一种双模结构，通过复用5/3提升结构离散小波变换与9/7提升结构离散小波变换来提高转换速率；参考文献[4]提出了递归金字塔结构；参考文献[5]采用了多相分解及系数折叠方法；等等。 
传统方法通过改进电路结构来降低内部存储器，但由于依然按照光栅扫描方法读入图像数据，本质上依旧是先进行水平方向的行变换，后进行垂直方向的列变换。所消耗的内部存储器并没有显著降低。传统方法对
Figure 18690DEST_PATH_IMAGE001
行
Figure 432878DEST_PATH_IMAGE001
列图像进行5/3提升结构离散小波变换时，所消耗的内部存储器最小为
Figure 201110148246X100002DEST_PATH_IMAGE015
，最多为。 
参考文献： 
[1] JPEG2000 Part 1: Final Draft, document ISO/IEC JTCI/SC29/WGI N1855.doc, International Standard (ISO/IEC FDISF DIS15444-1), Aug.2000.
[2] K. Andra, C. Chakrabarti, and T. Acharya, “A VLSI architecture for lifting-based forward and inverse wavelet transform,” IEEE Trans. on Signal Processing, vol. 50, no.4, pp. 966-977, 2002.
[3] G.-C Jung and S.-M. Park, “VLSI implement of lifting wavelet transform of JPEG2000 with efficient RPA (recursive pyramid algorithm) realization,”  IEICE Trans. on Fundamentals, vol. E88-A, no. 12, pp. 3508-3515, 2005.
[4] C.-T. Huang, P.-C. Tseng, and L.-G. Chen, “Efficient VLSI architecture of lifting-based discrete wavelet transform by systematic design method,” IEEE Int. Symposium Circuits and Systems, vol. 5, pp. 26-29, 2002.
[5] P.-C. Wu and L.-G. Chen, “An efficient architecture for two-dimensional discrete wavelet transform,” IEEE Trans. on Circuits and Systems for Video Technology, vol. 11, no.4, pp. 536-545, 2001.
[5] P.-Y. Chen, “VLSI implementation of discrete wavelet transform using the 5/3 filter,” IEICE Trans. on Information and Systems,vol. E85-D, no.12, pp. 1893-1897, 2002.
发明内容
技术问题：本发明的目的是提供一种用于图像压缩的提升结构二维离散小波变换交织扫描方法，该方法按照行列交织的顺序读入图像数据，同时进行水平方向和垂直方向的5/3提升结构一维离散小波变换，从而完成整幅图像的5/3提升结构二维离散小波变换。本发明在变换过程中，水平方向一维离散小波变换后输出的低频系数无需存储，仅需缓存N个高频系数，便可同时进行垂直方向一维离散小波变换，输出四个频率子带的小波系数。相对于传统最优的设计方案，消耗的内部存储器降低了50％。 
技术方案：本发明的用于图像压缩的提升结构二维离散小波变换交织扫描方法，按照行列交织的顺序读入N行N列图像数据，同时进行水平方向和垂直方向的5/3提升结构一维离散小波变换，其处理完成N行N列大小的图像仅需N个存储单元，其中N是自然数。下面，用5/3提升结构离散小波变换交织扫描方法对N行N列图像进行变换实例来阐明本发明的技术方案。定义N行N列图像中的第m行第n列的图像像素点为X(m, n)；定义水平方向一维离散小波变换后，第m行第n列高频系数为H（m, n），第m行第n列低频系数L(m, n)；定义垂直方向一维离散小波变换后，第m行第n列甚高频小波系数为HH(m, n)，第m行第n列高低频子带小波系数为HL(m, n)，第m行第n列低高频子带小波系数为LH(m, n)，第m行第n列甚低频子带小波系数为LL(m, n)。 
所述的按照行列交织的顺序读入N行N列的图像数据，是指从第0行第0列开始，读取第0行第0列图像数据X(0,0)、第0行第1列图像数据X(0,1)、第0行第2列图像数据X(0,2)进行水平方向一维离散小波变换，采用公式1和公式2算得第0行第0列高频系数H(0,0)、第0行第0列低频系数L(0,0)；公式1和公式2中，m和n都是小于等于N的自然数： 
           
Figure 106622DEST_PATH_IMAGE008
      公式1
                  公式2
           
Figure 281569DEST_PATH_IMAGE010
    公式3
              公式4
Figure 839775DEST_PATH_IMAGE012
      公式5
           
Figure 677281DEST_PATH_IMAGE013
     公式6
将H(0,0)存入内部存储器；换行读取第1行第0列图像数据X(1,0)、第1行第1列图像数据X(1,1)、第1行第2列图像数据X(1,2)进行水平方向一维离散小波变换，算得高频系数H(1,0)、低频系数L(1,0)并将H(1,0)存入内部存储器；按此方法逐行读入数据并进行转换直至第N行，并将相应的高频系数H(m, 0)存入内部存储器；
返回第0行，换列读取第0行第2列数据X(0,2)、第0行第3列数据X(0,3)、第0行第4列数据X(0,4)进行水平方向一维离散小波变换，计算出H(0,1)后，从内部存储器中读取H(0,0)用来计算L(0,1),之后将H(0,1)存入内部存储器并覆盖掉
Figure 201110148246X100002DEST_PATH_IMAGE017
；然后跳转到第1行读取第1行第2列数据X(1,2)、第1行第3列数据X(1,3)、第1行第4列数据X(1,4) 进行水平方向一维离散小波变换，计算出H(1,1)后，从内部存储器中读取H(1,0)用来计算L(1,1),之后将H(1,1)存入内部存储器并覆盖掉H (1,0)；按此方法逐行读取数据并进行转换直至第N行，并将原内部存储器中的H(m,0)逐个更新为H(m,1)；
返回第0行，循环进行以上操作，按照这种行列交织的方法读取全部图像数据进行计算。
所述同时进行水平方向和垂直方向的5/3提升结构离散小波变换，是水平方向一维离散小波变换所产生的高频系数和低频系数
Figure 201110148246X100002DEST_PATH_IMAGE019
是按照垂直方向逐列输出的；只需用三个寄存器缓存高频系数
Figure 683793DEST_PATH_IMAGE018
，便可用公式3和公式4进行垂直方向一维离散小波变换，生成甚高频子带小波系数HH(m,n)和高低频子带小波系数HL(m,n)；与此同时，用三个寄存器缓存高频系数，便可用公式5和6进行垂直方向一维离散小波变换，生成低高频子带小波系数LH(m,n)和甚低频子带小波系数LL(m,n)；即同时计算出四个频率子带的变换结果，从而完成5/3提升结构二维离散小波变换。 
所述的处理完成N行N列大小的图像仅需N个存储单元，是指从第0行第0列开始，读取第0行第0列图像数据X(0,0)、第0行第1列图像数据X(0,1)、第0行第2列图像数据X(0,2)进行水平方向一维离散小波变换，算得H(0,0)、L(0,0)并将H(0,0)存入内部存储器；换行读取第1行第0列图像数据X(1,0)、第1行第1列图像数据X(1,1)、第1行第2列图像数据X(1,2)进行水平方向一维离散小波变换，算得H(1,0)、L(1,0)并将H(1,0)存入内部存储器；按此方法逐行读入数据并进行转换直至第N行，并将N个高频系数H(m,0)存入内部存储器； 
返回第0行，换列读取第0行第2列数据X(0,2)、第0行第3列数据X(0,3)、第0行第4列数据X(0,4)进行水平方向一维离散小波变换，计算出H(0,1)后，从内部存储器中读取H(0,0)用来计算L(0,1),之后将H(0,1)存入内部存储器并覆盖掉；然后跳转到第1行读取第1行第2列数据X(1,2)、第1行第3列数据X(1,3)、第1行第4列数据X(1,4)进行水平方向一维离散小波变换，计算出H(1,1)后，从内部存储器中读取H(1,0)用来计算L(1,1),之后将H(1,1)存入内部存储器并覆盖掉    H (1,0)；按此方法逐行读取数据并进行转换直至第N行，并将原内部存储器中的N个H(m,0)逐个更新为N个H(m,1)；
返回第0行，循环进行以上操作，按照这种行列交织的方式读取全部图像数据进行完水平方向一维离散小波变换；水平方向一维离散小波变换所产生的高频系数
Figure 738971DEST_PATH_IMAGE002
和低频系数
Figure 201110148246X100002DEST_PATH_IMAGE021
是按照垂直方向逐列输出的，只需用三个寄存器缓存高频系数，便可用公式3和公式4进行垂直方向一维离散小波变换，生成甚高频子带小波系数HH(m,n)和高低频子带小波系数HL(m,n)；与此同时，用三个寄存器缓存低频系数，便可用公式5和公式6进行垂直方向一维离散小波变换，生成低高频子带小波系数LH(m,n)和甚低频子带小波系数LL(m,n)；即同时计算出四个频率子带的二维离散小波变换结果；在整个变换过程中，仅需N个存储单元。
有益效果：该发明突破了传统的光栅扫描方法所采用的先行后列的数据读入方式，按照行列交织扫描方法读入图像数据，使得5/3提升结构离散小波变换过程中，水平方向一维离散小波变换输出的低频系数无需存储，仅需缓存N个高频系数，便可同时进行垂直方向一维离散小波变换，同时输出四个频率子带的小波系数。该发明将内部存储器的消耗降低至，相对于传统最优的设计方案降低了50%。 
  
附图说明
下面结合附图和实施例对本发明做进一步说明。 
图1 采用传统的5/3提升结构二维离散小波变换光栅扫描方法处理6
Figure 551255DEST_PATH_IMAGE022
6图像的 
原理流程图 。                     
图2 用于图像压缩的提升结构二维离散小波变换交织扫描方法进行水平方向一
维离散小波变换时的数据读取和计算顺序。
图3用于图像压缩的提升结构二维离散小波变换交织扫描方法进行垂直方向一 
维离散小波变换时的数据读取和计算顺序。
图4用于图像压缩的提升结构二维离散小波变换交织扫描方法流程图。 
  
具体实施方式
用于图像压缩的提升结构二维离散小波变换交织扫描方法，按照行列交织的顺序读入图像数据，同时进行水平方向一维离散小波变换和垂直方向一维离散小波变换。其转换流程为： 
从第0行第0列开始，读取第0行第0列的X(0,0)、第0行第1列的X(0,1)、第0行第2列的X(0,2)进行水平方向一维离散小波变换，算得H(0,0)、L(0,0)并将H(0,0)存入内部存储器；换行读取X(1,0)、X(1,1)、X(1,2)进行水平方向一维离散小波变换，算得H(1,0)、L(1,0)并将H(1,0)存入内部存储器。按此方法逐行读入数据并进行转换直至第N行，并将相应的高频系数H(m,0)存入内部存储器。
返回第0行第2列，读取第0行第2列的X(0,2)、第0行第3列的X(0,3)、第0行第4列的X(0,4)进行水平方向一维离散小波变换，计算出H(0,1)后，从内部存储器中读取H(0,0)用来计算L(0,1),之后将H(0,1)存入内部存储器并覆盖掉
Figure 347041DEST_PATH_IMAGE017
。然后跳转到第1行读取X(1,2)、X(1,3)、X(1,4) 进行水平方向一维离散小波变换，计算出H(1,1)后，从内部存储器中读取H(1,0)用来计算L(1,1),之后将H(1,1)存入内部存储器并覆盖掉H (1,0)。按此方法逐行读取数据并进行转换直至第N行，并将原内部存储器中的H(m,0)逐个更新为H(m,1)。 
返回第0行，循环以上行列交织的读取操作和计算过程，直至全部图像数据进行完水平方向一维离散小波变换。在此过程中，水平方向一维离散小波变换所产生的高频系数
Figure 458217DEST_PATH_IMAGE002
和低频系数
Figure 468898DEST_PATH_IMAGE021
是按照垂直方向逐列输出的。只需用三个寄存器缓存高频系数
Figure 623805DEST_PATH_IMAGE002
，便可用公式3和公式4进行垂直方向一维离散小波变换，生成HH(m,n)和HL(m,n)；与此同时，用三个寄存器缓存高频系数
Figure 657620DEST_PATH_IMAGE003
，便可用公式5和公式6进行垂直方向一维离散小波变换，生成LH(m,n)和LL(m,n)。即同时计算出四个频率子带的二维离散小波变换结果。 
基于以上分析，具体的算法流程为： 
1）初始化m=0, n=0；
2) 读入原始图像数据X(m, 0)、X(m, 1)、X(m, 2)，由公式(1)、公式(2)分别计算出中间系数H(m, 0)、L(m, 0)，并将H(m, 0)存入内部存储器。将m加1；
3）重复步骤2)，至m=N，将m置零并将n加1；
4）读入原始图像数据X(m, 2n)、X(m, 2n+1)、X(m, 2n+2)，由公式(1)计算出H(m,n)从内部存储器中读入H(m, n-1)，由公式(2)计算出L(m, n)。将H(m, n)存入内部存储器，覆盖掉H(m, n-1)，将m加1；
5）重复步骤4)至m=N，将m置零并将n加1；
6）重复步骤4)和5)至n=N/2-1，计算结束；
水平方向一维离散小波变换和垂直方向一维离散小波变换是并行的计算过程。因为水平方向一维离散小波变换输出数据的顺序是垂直方向逐列输出，所以垂直方向一维离散小波变换不需要内部存储器，而只需3个寄存器缓存三个中间变换系数后，就直接计算并输出最终的变换结果。
 下面以
Figure 201110148246X100002DEST_PATH_IMAGE023
图像进行提升结构二维小波变换的实例，来说明具体实施方式： 
1）读入原始图像第一行前三个数据X(0,0)、X(0,1)、X(0,2)，计算出中间系数H(0,0)、L(0,0)，并将H(0,0)存入内部存储器；
2）读入原始图像第二行前三个数据X(1,0)、X(1,1)、X(1,2)，计算出中间系数H(1,0)、L(1,0)，并将H(1,0)存入内部存储器；
3）读入原始图像第三行前三个数据X(2,0)、X(2,1)、X(2,2)，计算出中间系数H(2,0)、L(2,0)，并将H(2,0)存入内部存储器；
4）第二级变换单元由H(0,0)、L(0,0)， H(1,0)、L(1,0)， H(2,0)、L(2,0)，计算并输出HH(0,0)、HL(0,0)、LH(0,0) 、 LL(0,0)；
5）读入原始图像第四行前三个数据X(3,0)、X(3,1)、X(3,2)，计算出中间系数H(3,0)、L(3,0)，并将H(3,0)存入内部存储器；
6）读入原始图像第五行前三个数据X(4,0)、X(4,1)、X(4,2)，计算出中间系数H(4,0)、L(4,0)，并将H(4,0)存入内部存储器；
7）第二级变换单元由H(2,0)、L(2,0)， H(3,0)、L(3,0)， H(4,0)、L(4,0)，计算并输出HH(1,0)、HL(1,0)、LH(1,0)、LL(1,0)；
8）读入原始图像第一行后三个数据X(0,2)、X(0,3)、X(0,4)，计算出中间系数H(0,1)，从内部存储器中读入H(0,0)，计算出L(0,1)。将H(0,1)存入内部存储器，覆盖掉H(0,0)；
9）读入原始图像第二行后三个数据X(1,2)、X(1,3)、X(1,4)，计算出中间系数H(1,1)，从内部存储器中读入H(1,0)，计算出L(1,1)。将H(1,1)存入内部存储器，覆盖掉H(1,0)；
10）读入原始图像第三行后三个数据X(2,2)、X(2,3)、X(2,4)，计算出中间系数H(2,1)，从内部存储器中读入H(2,0)，计算出L(2,1)。将H(2,1)存入内部存储器，覆盖掉H(2,0)；
11）第二级变换单元由H(0,1)、L(0,1)， H(1,1)、L(1,1)， H(2,1)、L(2,1)，计算并输出HH(0,1)、HL(0,1)、LH(0,1)、LL(0,1)；
12）读入原始图像第四行后三个数据X(3,2)、X(3,3)、X(3,4)，计算出中间系数H(3,1)，从内部存储器中读入H(3,0)，计算出L(3,1)。将H(3,1)存入内部存储器，覆盖掉H(3,0)；
13）读入原始图像第四行后三个数据X(4,2)、X(4,3)、X(4,4)，计算出中间系数H(4,1)，从内部存储器中读入H(4,0)，计算出L(4,1)。将H(4,1)存入内部存储器，覆盖掉H(4,0)；
14）第二级变换单元由H(2,1)、L(2,1)， H(3,1)、L(3,1)， H(4,1)、L(4,1)，计算并输出HH(1,1)、HL(1,1)、LH(1,1)、LL(1,1)，计算结束。  

Claims (4)

1.一种用于图像压缩的提升结构二维离散小波变换交织扫描方法，其特征在于按照行列交织的顺序读入N行N列图像数据，同时进行水平方向和垂直方向的5/3提升结构一维离散小波变换，其处理完成N行N列大小的图像仅需N个存储单元，其中N是自然数。
2.如权利要求1所述的用于图像压缩的提升结构二维离散小波变换交织扫描方法，其特征在于所述的按照行列交织的顺序读入N行N列的图像数据，是指从第0行第0列开始，读取第0行第0列图像数据X(0，0)、第0行第1列图像数据X(0，1)、第0行第2列图像数据X(0，2)进行水平方向一维离散小波变换，采用公式1和公式2算得第0行第0列高频系数H(0，0)、第0行第0列低频系数L(0，0)；公式1和公式2中，m和n都是小于等于N的自然数：
Figure FDA0000065978870000011
公式1
Figure FDA0000065978870000012
公式2
公式3
Figure FDA0000065978870000014
公式4
Figure FDA0000065978870000015
公式5
Figure FDA0000065978870000016
公式6
将H(0，0)存入内部存储器；换行读取第1行第0列图像数据X(1，0)、第1行第1列图像数据X(1，1)、第1行第2列图像数据X(1，2)进行水平方向一维离散小波变换，算得高频系数H(1，0)、低频系数L(1，0)并将H(1，0)存入内部存储器；按此方法逐行读入数据并进行转换直至第N行，并将相应的高频系数H(m，0)存入内部存储器；
返回第0行，换列读取第0行第2列数据X(0，2)、第0行第3列数据X(0，3)、第0行第4列数据X(0，4)进行水平方向一维离散小波变换，计算出H(0，1)后，从内部存储器中读取H(0，0)用来计算L(0，1)，之后将H(0，1)存入内部存储器并覆盖掉H(0，0)；然后跳转到第1行读取第1行第2列数据X(1，2)、第1行第3列数据X(1，3)、第1行第4列数据X(1，4)进行水平方向一维离散小波变换，计算出H(1，1)后，从内部存储器中读取H(1，0)用来计算L(1，1)，之后将H(1，1)存入内部存储器并覆盖掉H(1，0)；按此方法逐行读取数据并进行转换直至第N行，并将原内部存储器中的H(m，0)逐个更新为H(m，1)；
返回第0行，循环进行以上操作，按照这种行列交织的方法读取全部图像数据进行计算。
3.如权利要求1所述的用于图像压缩的提升结构二维离散小波变换交织扫描方法，其特征在于所述同时进行水平方向和垂直方向的5/3提升结构离散小波变换，是水平方向一维离散小波变换所产生的高频系数H(m，n)和低频系数L(m，n)是按照垂直方向逐列输出的；只需用三个寄存器缓存高频系数H(m，n)，便可用公式3和公式4进行垂直方向一维离散小波变换，生成甚高频子带小波系数HH(m，n)和高低频子带小波系数HL(m，n)；与此同时，用三个寄存器缓存高频系数L(m，n)，便可用公式5和公式6进行垂直方向一维离散小波变换，生成低高频子带小波系数LH(m，n)和甚低频子带小波系数LL(m，n)；即同时计算出四个频率子带的变换结果，从而完成5/3提升结构二维离散小波变换。
4.如权利要求1所述的用于图像压缩的提升结构二维离散小波变换交织扫描方法，其特征在于所述的处理完成N行N列大小的图像仅需N个存储单元，是指从第0行第0列开始，读取第0行第0列图像数据X(0，0)、第0行第1列图像数据X(0，1)、第0行第2列图像数据X(0，2)进行水平方向一维离散小波变换，算得H(0，0)、L(0，0)并将H(0，0)存入内部存储器；换行读取第1行第0列图像数据X(1，0)、第1行第1列图像数据X(1，1)、第1行第2列图像数据X(1，2)进行水平方向一维离散小波变换，算得H(1，0)、L(1，0)并将H(1，0)存入内部存储器；按此方法逐行读入数据并进行转换直至第N行，并将N个高频系数H(m，0)存入内部存储器；
返回第0行，换列读取第0行第2列数据X(0，2)、第0行第3列数据X(0，3)、第0行第4列数据X(0，4)进行水平方向一维离散小波变换，计算出H(0，1)后，从内部存储器中读取H(0，0)用来计算L(0，1)，之后将H(0，1)存入内部存储器并覆盖掉H(0，0)；然后跳转到第1行读取第1行第2列数据X(1，2)、第1行第3列数据X(1，3)、第1行第4列数据X(1，4)进行水平方向一维离散小波变换，计算出H(1，1)后，从内部存储器中读取H(1，0)用来计算L(1，1)，之后将H(1，1)存入内部存储器并覆盖掉H(1，0)；按此方法逐行读取数据并进行转换直至第N行，并将原内部存储器中的N个H(m，0)逐个更新为N个H(m，1)；
返回第0行，循环进行以上操作，按照这种行列交织的方式读取全部图像数据进行完水平方向一维离散小波变换；水平方向一维离散小波变换所产生的高频系数H(m，n)和低频系数L(m，n)是按照垂直方向逐列输出的，只需用三个寄存器缓存高频系数，便可用公式3和公式4进行垂直方向一维离散小波变换，生成甚高频子带小波系数HH(m，n)和高低频子带小波系数HL(m，n)；与此同时，用三个寄存器缓存低频系数，便可用公式5和公式6进行垂直方向一维离散小波变换，生成低高频子带小波系数LH(m，n)和甚低频子带小波系数LL(m，n)；即同时计算出四个频率子带的二维离散小波变换结果；在整个变换过程中，仅需N个存储单元。
CN 201110148246 2011-06-02 2011-06-02 用于图像压缩的提升结构二维离散小波变换交织扫描方法 Expired - Fee Related CN102281437B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110148246 CN102281437B (zh) 2011-06-02 2011-06-02 用于图像压缩的提升结构二维离散小波变换交织扫描方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110148246 CN102281437B (zh) 2011-06-02 2011-06-02 用于图像压缩的提升结构二维离散小波变换交织扫描方法

Publications (2)

Publication Number Publication Date
CN102281437A true CN102281437A (zh) 2011-12-14
CN102281437B CN102281437B (zh) 2013-01-09

Family

ID=45106575

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110148246 Expired - Fee Related CN102281437B (zh) 2011-06-02 2011-06-02 用于图像压缩的提升结构二维离散小波变换交织扫描方法

Country Status (1)

Country Link
CN (1) CN102281437B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102572429A (zh) * 2011-12-29 2012-07-11 东南大学 一种二维离散小波变换的硬件架构
CN102685501A (zh) * 2012-05-14 2012-09-19 西安电子科技大学 实现jpeg2000图像压缩的定点小波变换方法
CN111684484A (zh) * 2019-07-30 2020-09-18 深圳市大疆创新科技有限公司 Dwt运算装置、方法、图像处理装置和可移动平台

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6178269B1 (en) * 1998-08-06 2001-01-23 Intel Corporation Architecture for computing a two-dimensional discrete wavelet transform
US20030046322A1 (en) * 2001-06-01 2003-03-06 David Guevorkian Flowgraph representation of discrete wavelet transforms and wavelet packets for their efficient parallel implementation
CN1589021A (zh) * 2004-09-28 2005-03-02 华中科技大学 一种并行二维离散小波变换电路
CN101430737A (zh) * 2008-11-19 2009-05-13 西安电子科技大学 提升小波变换的vlsi结构设计方法
CN101534439A (zh) * 2008-03-13 2009-09-16 中国科学院声学研究所 一种低功耗并行的小波变换的vlsi结构

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6178269B1 (en) * 1998-08-06 2001-01-23 Intel Corporation Architecture for computing a two-dimensional discrete wavelet transform
US20030046322A1 (en) * 2001-06-01 2003-03-06 David Guevorkian Flowgraph representation of discrete wavelet transforms and wavelet packets for their efficient parallel implementation
CN1589021A (zh) * 2004-09-28 2005-03-02 华中科技大学 一种并行二维离散小波变换电路
CN101534439A (zh) * 2008-03-13 2009-09-16 中国科学院声学研究所 一种低功耗并行的小波变换的vlsi结构
CN101430737A (zh) * 2008-11-19 2009-05-13 西安电子科技大学 提升小波变换的vlsi结构设计方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
PO-CHENG WU ET AL.: "An efficient architecture for two-dimensional discrete wavelet transform", 《IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS FOR VIDEO TECHNOLOGY》, vol. 11, no. 4, 30 April 2001 (2001-04-30), pages 536 - 545 *
丁媛媛等: "基于提升小波变换的图像压缩编码的VLSI实现", 《吉林大学学报(工学版)》, vol. 37, no. 03, 31 May 2007 (2007-05-31), pages 675 - 680 *
陈磊等: "一种用于JPEG2000的小波变换VLSI结构设计方法", 《中国图象图形学报》, vol. 12, no. 10, 31 October 2007 (2007-10-31), pages 1730 - 1734 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102572429A (zh) * 2011-12-29 2012-07-11 东南大学 一种二维离散小波变换的硬件架构
CN102572429B (zh) * 2011-12-29 2014-07-09 东南大学 一种二维离散小波变换的硬件架构
CN102685501A (zh) * 2012-05-14 2012-09-19 西安电子科技大学 实现jpeg2000图像压缩的定点小波变换方法
CN102685501B (zh) * 2012-05-14 2014-11-12 西安电子科技大学 实现jpeg2000图像压缩的定点小波变换方法
CN111684484A (zh) * 2019-07-30 2020-09-18 深圳市大疆创新科技有限公司 Dwt运算装置、方法、图像处理装置和可移动平台
WO2021016893A1 (zh) * 2019-07-30 2021-02-04 深圳市大疆创新科技有限公司 Dwt运算装置、方法、图像处理装置和可移动平台

Also Published As

Publication number Publication date
CN102281437B (zh) 2013-01-09

Similar Documents

Publication Publication Date Title
US6671413B1 (en) Embedded and efficient low-complexity hierarchical image coder and corresponding methods therefor
CN101039430B (zh) 一种视频编码中对残差矩阵的快速扫描方法
Fang et al. Interpolation-based direction-adaptive lifting DWT and modified SPIHT for image compression in multimedia communications
Hsia et al. Improved low-complexity algorithm for 2-D integer lifting-based discrete wavelet transform using symmetric mask-based scheme
CN101631243B (zh) 一种基于小波变换的图像编和解码的方法
CN108810534B (zh) 物联网下基于方向提升小波及改进spiht的图像压缩方法
CN102833546A (zh) 基于小波子带交织最优量化的高速图像压缩方法及装置
CN102281437B (zh) 用于图像压缩的提升结构二维离散小波变换交织扫描方法
CN102223534B (zh) 用于图像压缩的全并行位平面编码方法
CN109168005A (zh) 一种提高存储效率的二维离散小波变换结构设计方法
Jin et al. Post-processing for intra coding through perceptual adversarial learning and progressive refinement
CN110740333B (zh) 基于小波模极大值重建的改进spiht图像编码和解码方法
CN102333222B (zh) 二维离散小波变换电路及应用该电路的图像压缩方法
Liu et al. An improved SPIHT algorithm for image compression in low bit rate
Modrzyk et al. A high-performance architecture of JPEG2000 encoder
CN102572423B (zh) 一种基于重要性概率平衡树的视频编码方法
Zhu et al. An improved SPIHT algorithm based on wavelet coefficient blocks for image coding
Bin et al. An improved SPIHT wavelet transform in the underwater acoustic image compression
CN103402090A (zh) 一种基于霍夫曼无损压缩算法的改进方法
Puttaraju et al. FPGA implementation of 5/3 integer dwt for image compression
CN110365990A (zh) 一种窄带准无损视频编码系统
Talukder et al. Enhancement of discrete wavelet transform (dwt) for image transmission over internet
CN109300086A (zh) 基于清晰度的图像分块方法
CN103533357B (zh) 基于fpga的多级复用小波变换器
CN103561275A (zh) 一种改进条带波变换的图像压缩方法

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: 20130109

Termination date: 20150602

EXPY Termination of patent right or utility model