多级二维整数小波变换的FPGA实现研究

·PAGE34· 电子测量与仪器学报 2009年第7期 多级二维整数小波变换的FPGA实现研究 ·PAGE33·第23卷第7期 电子测量与仪器学报 Vol.23No.72009年7月 ·PAGE29·本文于2009年2月收到。*基金项目:湖南省自然科学基金(编号:07JJ3118)资助项目更多电子资料请登录赛微电子网多级二维整数小波变换的FPGA实现研究*谭会生(湖南工业大学电气与信息工程学院,株洲412008)摘要:为了满足整数小波变换实时应用的需要,研究了整数小波变换的FPGA实现问题。相对于DSP等传统实现方式,用FPGA实现整数小波变换具有处理速度快,可重新配置硬件,易于修改移植等优点。论文首先描述了二维(5,3)整数小波变换的算法,接着阐述了一种多级二维(5,3)整数小波变换的FPGA实现结构,最后给出了硬件资源消耗、最大时钟频率和功能测试结果等FPGA实现结果。为了提高系统的处理速度,降低系统的资源消耗,本设计采用了参数可配置、共享一维小波变换单元等方法进行结构优化。实验结果证明了本设计结构的有效性,逻辑功能的正确性,修改移植方便,具有良好的应用价值。关键词:图像无损压缩；多级整数小波变换；FPGA实现；共享硬件资源；参数可配置结构中图分类号:TP391文献标识码:A国家标准学科分类代码:510.4050FPGAimplementationofmultileveltwo-dimensionintegerwavelettransformTanHuisheng(SchoolofElectricalandInformationEngineering,HunanUniversityofTechnology,Zhuzhou412008,China)Abstract:Inordertomeetthedemandofreal-timeapplicationsofintegerwavelettransform(IWT),aFPGAimplementationofIWTisresearched.ComparingwithclassicalimplementationsuchasDSPetc,theFPGAimplementationofIWThastheadvantagessuchashigherprocessingspeed,reconfigurablehardware,moreeasytomodifyandtransplant.Firstofall,thispaperdescribesatwo-dimensional(5,3)IWTalgorithm,thenproposesaFPGAimplementationstructureofamultileveltwo-dimensional(5,3)IWT,andfinallypresentsitsexperimentalresultssuchashardwareresourceconsumption,maximumclockfrequencyandfunctiontestingresults.Inordertoimprovetheprocessingspeed,andreducetheresourceconsumption,theproposedstructureisoptimizedbyusingconfigurableparameter,sharingone-dimensionalwavelettransformunitetc.Theexperimentalresultsshowthatthestructureofthisdesignisavailable,thelogicfunctioniscorrect,themodificationandtransplantismoreeasy,anditisofgreatvalueinapplication.Keywords:imagelosslesscompression;multilevelintegerwavelettransform;FPGAimplementation;sharinghardwareresources;configurableparameterstructure1引言随着社会的不断发展和进步,出于资源共享和远程测控的需要,出现了基于互联网Internet的远程测控系统[1-2]。而远程测控系统一般包括共享设备端、服务器和用户端3个组成部分,用户端主要是执行现场信息的采集、反馈信息的接受及有关信息的处理等任务,其测控信息大多是视频信息,要保证其正常而有效的工作,必须解决两个关键问题:1)高速实时测控;2)视频信息的高效压缩处理。传统上的电子系统硬件实现主要有三种途径:实现图像处理通用处理器、数字信号处理器(DSP)、专用集成电路(ASIC)。但是通用处理器数据处理速度慢,不能满足高速处理的要求。DSP拥有流水线结构和优化的算法,可以加速图像处理,但不适合所有运算。ASIC功能专一,造价昂贵,适应性差。因此传统上的电子系统硬件实现途径已不能满足远程测控系统的要求。随着微电子和计算机技术的快速发展,出现了现场可编程门阵列FPGA器件,它具有结构灵活、设计周期短、密度高、体积小、性能好等优点,因此近年国内外比较流行用在FPGA芯片来实现一些复杂的电子系统[3-5],也为远程测控系统的应用和发展提供了可行的技术基础。小波变换作为一种全局变换,在时域和频域都具有良好的局部化性能,而且在应用中易于考虑人类的视觉特性,从而成为图像压缩编码的主要技术之一,并且离散小波变换已经纳入MPEG-4和JPEG2000编码标准。但小波变换需要大量的计算,因此小波变换的硬件实现成为小波实时应用的一个基础而关键的问题。JPEG2000中推荐了2种小波滤波器基于整数的(5,3)滤波器和基于实数的(9,7)滤波器,它们分别用于无损和有损压缩。基于提升方法的(5,3)整数小波变换不仅具有计算更快捷,能够在当前位置完成小波变换从而节省内存,能对任意尺寸图像进行小波变换等优点,而且可从整数变换系数中完全重建图像,因而可实现图像无损压缩[6]。对于小波变换的VLSI实现,有关学者进行了大量的研究工作,提出了资源共享、并行技术、流水线技术、折叠技术等设计优化技术及应用实例[7-12],但这些研究工作主要集中在一维和二维小波变换的单级实现上,对于多级小波变换的实现则研究很少,并且一般采用多个硬件的串级实现。本文就是基于参数可配置、资源共享和二维小波变换的分离技术来探讨基于提升方法的多级二维(5,3)整数小波正变换的FPGA实现问题,其重点是FPGA实现结构设计及其VHDL程序的设计与验证。2多级二维整数小波变换算法的描述根据多级二维金字塔小波分解的一般原理和二维整数提升模式算法的表示方法[6,8],可将多级二维(5,3)整数小波变换算法描述如公式(1)～(6)所示,式中的l,c,r分别表示各系数所属分解级数,列索引,行索引。行变换: H: (1)L: (2)列变换: HH: (3) LH: (4)HL: (5) LL: (6)3多级二维整数小波变换的结构设计3.1二维整数小波变换的总体结构设计分析公式(1)～(6)可看出,二维(5,3)整数小波变换的各步数据处理,其数据处理的本质是一样的,只是在处理过程中所处理的对象不同而已,行变换是对原始图像数据或上一级变换的低频数据,而列变换则是针对本级变换得到的高频数据和低频数据。同时,在小波变换过程中对边界数据需要进行延拓处理,JPEG2000标准中推荐采用对称的延拓算法。根据基于提升结构小波的运算特点,本文拟采用嵌入式数据延拓算法,它通过在起始与结尾端改变小波的计算过程,将数据的延拓嵌入到数据的计算过程中去,可避免进行数据延拓所需要的额外的运算和存储单元。因此可将二维(5,3)小波变换中的数据处理用一个通用的基于镜像对称延拓边界处理的一维变换算法表示公式(7)～(8)所示[6],式中k为偶数,xi表示输入数据,y2i+1、y2i分别表示变换后的高频成分H和低频成分L。 (7) (8)对于多级二维(5,3)整数小波正变换,为了减少数据处理的复杂性,可利用其变换的可分离特性,将其分开成2个独立的阶段来实现行变换和列变换,同时为了共享数据处理硬件,可设计一个行变换和列变换共享的一维(5,3)整数小波变换核,并设计一个可控制系统变换级数的参数,使系统可根据需要进行多级变换。根据以上设计思想,可得到二维(5,3)整数小波正变换的总体结构如图1所示。系统的工作原理如下:用户输入相关的参数到用户设置单元后,首先水平变换控制单元控制(5,3)小波正变换核按行完成小波的行变换,接着启动垂直方向的变换,垂直变换控制单元控制(5,3)小波正变换核按列完成小波的所有列变换,再由变换级数控制单元控制是否进行下一级的变换,若要进行下一级变换的话,则重复上述的行变换和列变换。3.2变换数据存储器的读写控制若变换的图像为N×M个像素,则需要2块能存储N×M个像素的存储器,其中一块用于存取源图像和结果图像数据,另一块则用于存取变换的临时数据。图2为小波变换过程中数据存储器的读写示意图。图中的箭头表示数据的读取和存放方向。3.3一维整数小波变换模块DWT的设计3.3.1一维整数小波变换数据处理分析根据嵌入式数据延拓(5,3)整数小波变换算法,可得到一维(5,3)小波变换数据处理原理如图3所示。图1多级二维正数小波变换的总体结构框图Fig.1Blockdiagramofmultileveltwo-dimensionintegerwavelettransform图2多级小波变换数据存储器的读写示意图Fig.2Schematicdiagramofdataread/writtenformultilevelintegerwavelettransform图3一维整数小波变换数据处理原理图Fig.3Blockdiagramofdataprocessingforone-dimensionintegerwavelettransformalgorithm3.3.2一维整数小波变换的控制根据嵌入式数据延拓(5,3)整数小波变换算法,并考虑到基于VHDL进行硬件设计的需要,可用一个状态机来控制一维(5,3)小波变换,如图4所示。该状态机不但控制小波变换的计算,还控制着变换过程中对存储器的读写操作。其中变换准备状态包括图像源数据首地址、图像数据地址步进、低频系数存放首地址、高频系数存放首地址、低频系数地址步进、高频系数地址步进和像素数目的设置。3.4系统控制模块DWT2D的设计系统控制模块DWT2D可通过一个状态机来完成其功能,该状态机如图5所示。图4一维整数小波变换的控制状态机Fig.4Controlstateplaneofone-dimensionintegerwavelettransformalgorithm图5多级二维整数小波变换控制状态机Fig.5Controlstateplaneofmultileveltwo-dimensionintegerwavelettransform4多级二维整数小波变换的实现结果根据前述的结构设计,首先进行了各模块的VHDL程序设计和顶层测试文件testdwt2d.vhd的设计,并采用ModelSimSE6.0进行仿真。测试数据是一幅64×64的Lena图像,图像的数据已存在文件testdata.txt上,在运行仿真程序时testdata被以文件的形式读到RAM中进行变换。在仿真结束后,小波变换的结果被从RAM中读出来写到文件DWTResult.txt上。使用ModelSimSE6.0进行仿真的部分波形如图6所示(限于篇幅,这里只给出了第一级行变换和列变换的仿真结果)。图7是在用ModelSimSE6.0仿真时对一幅64×64的Lena头像所做的分解实验。首先将Lena原始图像的数据存放在TestData.txt,在设置好分解的级数后,运行ModelSimSE6.0对testDWT2D.vhd进行仿真,仿真后的结果被存放在DWTResult.txt上,然后将DWTResult.txt上的数据变成图像可得图7(b)和图7(c)。(a)第一级行变换仿真结果(b)第一级列变换仿真结果图6多级二维整数小波变换部分仿真结果Fig.6Simulationresultofmultileveltwo-dimensionintegerwavelettransform(a)Lena原始图像(b)一级小波分解数据图像(c)两级小波分解数据图像图7Lena头像的小波分解数据合成的图像Fig.7WaveletdecompositiondatasynthesisimageofLenahead5结论本文在分析多级二维(5,3)整数小波变换各种数据运算本质的基础上,提出了一种基于参数可配置和共享一维小波变换核的多级二维(5,3)整数小波变换的FPGA实现结构,并进行了VHDL的程序设计与验证。实验结果表明本结构可行、有效、正确,并且易于修改和移植,具有良好的应用价值。参考文献:唐进,马树元,吴平东,等.基于Internet的共享测控系统实例[J].仪器仪表学报,2008,29(6):1319-1322.TANGJ,MASHY,WUPD,etal.Exampleofinternet-basedsharingmeasurementandcontrolsystem[J].ChineseJournalofScientificInstrument,2008,29(6):1319-1322.许丽丽.基于Internet的远程测量系统实现方案研究[J].电子测量与仪器学报,2008,22(2):52-55.XULL.Researchonsolutiontointernetbasedmeasurementsystem[J].JournalofElectronicMeasurementandInstrument,2008,22(2):52-55.TESSIERR,BURLESONW.Reconfigurablecomputingfordigitalsignalprocessing:Asurvey[J].JournalofVLSISignalProcessing,2001,28:7-27.tuningPIDControllerBasedonFPGA[J].ChineseJournalofScientificInstrument,2005,26(8):833-837.谭会生,张昌凡.EDA技术及应用(第二版)[M].西安:西安电子科技大学出版社,2004.TANHSH,ZHANGCHF.EDATechnologyandApplication(2ndEdition)[M].XI’AN:XidianUniversityPress,2004.谭会生.(5,3)整数小波正变换的FPGA实现研究[J].微计算机信息(嵌入式与SOC),2007,23(20):176-177,124.谭会生TANHSH.FPGAImplementationofIntegerWaveletForwardTransformwith(5,3)Filters[J].MicrocomputerInformation,2007,23(20):176-177,124.GRANGETTOM,MAGLIE,MATINAM,etal.Optimizationandimplementationoftheintegerwavelettransformforimagecoding[J]．IEEETrans.onImageProcessing,2002,11(6):596-604.FERRETTIM,RIZZOD.Aparallelarchitectureforthe2-Ddiscretewavelettransformwithintegerliftingscheme[J].JournalofVLSISignalProcessing,2001,28:165-185.PREMKUMARAB.AnefficientVLSIarchitectureforthecomputationof1-Ddiscretewavelettransform[J].JournalofVLSISignalProcessing2002,31:231-241.TAIPL,LIUCT,WANGJSH.Anintegratedsystolicarraydesignforvideocompression[J].JournalofVLSISignalProcessing,2003,33:157-169.LIMQUECOJC,BAYOUMIMA.AVLSIarchitectureforsepara