数字视频降分辨率算法的滤波器设计

数字视频降分辨率算法的滤波器设计-电气论文数字视频降分辨率算法的滤波器设计袁莉，初秀娟，王志国（绥化学院，黑龙江绥化152061）摘要数字视频在不通终端上应用时往往需要将分辨率降低到原始分辨率的1/2或1/4因而需要有合适的算法改善降分辨率后图像的清晰度，满足用户对显示图像质量的要求基于视频转码的需要，针对分辨率转换算法中最常用的滤波与子采样法，设计了一种用于滤波与子采样法的5阶数字滤波器，并在编译器上进行仿真验证与评估仿真结果显示，新滤波器表现出优异的峰值信噪比，并且在分块算法中依然优异总之，这里提出的方法不但复杂度适中，适应性强，而且能够高效、高质量地完成数字视频降分辨率关键词视频转码；分辨率；滤波与子采样法；滤波器中图分类号TN91334文献标识码A文章编号1004373X

（2015）20002605DesignoffilterforresolutionreductionalgorithmusedindigitalvideoYUANLi，CHUXiujuan，WANGZhiguo（SuihuaUniversity，Suihua152061，China）Abstract Whenthedigitalvideoisusedindifferentterminals，theresolutionisoftenreducedbyafactorof2or4，soasuitablealgorithmisneededtoimprovetheimagesharpnessafterresolutionreductiontomeetuserdemandsforimagequality.Inviewoftheneedofvideotranscoding，afifthorderdigitalfilterusedforfilteringandsubsamplingmethodwasdesigned，whichwasassessedandverifiedoninterpreter.ThesimulationresultsshowthatthenewfilterhashighPSNR（peaksignaltonoiseratio），andisstillexcellentintheblockalgorithm.Inshort，theproposedmethodhasmoderatecomplexityandstrongadaptability，andcanachievehighefficiencyandqualitydigitalvideoresolutionreduction.Keywords videotranscoding；resolution；filterandsample；filter各种数字视频实际应用在各种不同的终端设备上，如电脑、平板电脑、电视、手机等，他们支持不同分辨率和不同压缩标准的视频，因此经常存在对同一视频进行不同分辨率间转换的需求例如使用手持终端设备接收观看运营商提供的标准清晰度数字电视节目，往往需要将节目分辨率降低到原始分辨率的12或14因而需要有合适的算法改善降分辨率后图像的清晰度，满足人们对显示图像质量的要求对于这样的实际需求，滤波子采样算法降分辨率的算法和物理结构实现相对简单，适合在软硬件系统中集成但是此方法的滤波器设计尤为关键，需要大量的数学和信号系统相关的建模和演算设计一个性能良好的滤波子采样降分变率方法，并且进行全面系统地测试，对分辨率转换相关软硬件设计有指导意义另一方面，国内外大部分算法都是在像素域分块后在变换域（如AVS域、DCT域等）进行，保证了速度但因为分块操作引入了峰值信噪比（PSNR）的损失（称之为块效应）对块效应的分析评估可以在转换速度与精度中找到平衡点，可以有效地指导以及优化各种算法的设计和实现，为图像领域未来的各种算法设计、程序设计和硬件设计都有深远影响本文基于视频转码的需要，设计并实现了一种应用于此方法的5阶数字滤波器，并在编译器上进行仿真验证与评估，取得了令人满意的结果此滤波器在像素域8×8分块后使用，表现依然良好1滤波与子采样法滤波与子采样法降分辨率基本过程如图1所示将图像X的行与列看做数字序列，用一个系统函数对其滤波，滤去高频分量，得到图像Y，然后等间隔取出像素点组成新的图像这个过程中滤波后得到的m×n像素图像（见图1中

②）有34的像素是无用的，因为采样只取14的像素所以只要确定了采样点，则只需对采样点进行滤波即可，可以省去很多计算，过程如图2所示综上所述，滤波与子采样降分变率方法的第1步为滤波，对于k阶的滤波与子采样降分变率方法，所用到的低通滤波器系统的方程为滤波与子采样降分变率方法的第2步子采样，及对变换后的Y等间隔采样对于m×n到m/2×n/2的降分辨率，则要求每2行取出1行，每2列取出1列，组成新的图像行与列可以分开考虑成一维的序列，如果滤波器系统是的低通滤波器，输入序列，则滤波后序列是以1阶滤波子采样为例，像素域平均法是最简单形式其实现为，每2×2矩阵的4个元素取均值，得到1个新的像素值，代替原来的4个像素本质是对行和列同时做1阶低通滤波并采样在一个维度上，此低通滤波器的系统函数是，低通滤波系统的差分方程是长度为n的序列，滤波后得到的新序列长度仍为n，对其进行每2像素点取一点采样，得到的12n像素序列即为降分变率后的序列对于二维图像，把系统函数写作H1D=[1/21/2]，则二维的系统函数为在实际实现中，其实并不需要对每一个像素点进行滤波得到m×n像素图像再采样，滤波只要对采样点进行即可所以，像素域平均法降分辨率具体实现就是每2×2矩阵的4个元素取均值，得到1个新的像素值，代替原来的4个像素以2阶滤波子采样为例对连续的一行中的3个点，或对连续的3行进行变换示，为一低通数字滤波器二维具体实现，对每3×3个像素乘以如下矩阵得到一个新的像素点，即可同时完成对列和行的滤波，此时2滤波器设计

2.1设计滤波器阶数简化考虑，先低通滤波再采样，先只考虑空间一维，因为二维只是2个一维变换的简单结合首先讨论不同阶数滤波后的最佳采样结果，以1~4阶为例，这些阶数的最小偏移情形如图4所示因此，滤波之后必须选择合理的采样方法，尽可能减小图像偏移对于折半降分辨率，最好使用奇数阶的低通滤波器来滤波偶数阶滤波后无论如何采样，都会引起图像中信息的偏移举一个简单实例，图5中

①中背景是模拟图像，·表示原采样点（即数字图像像素点），×表示新采样点（降分变率后图像像素点）方框表示降分变率数字图像一个采样点（红色×）所涵盖的图片信息图5中

②为1阶（或2抽头，或像素平均值法）滤波，采样点无偏；图5中

③为2阶（或3抽头），采样点偏移；图5中

④为3阶（或4抽头），采样点无偏；图5中

⑤为4阶（或5抽头），采样点偏移可以看出，对于二维情形，依然是奇数阶的低通滤波器偏移较小考虑无偏性和算法复杂度，5阶滤波（6抽头）

2.2设计滤波器系数理想的数字低通滤波器幅频特性为矩形窗，如图6所示根据傅里叶变换，空域Sa函数，则频域为矩形窗对Sa函数采样，则可得到矩形窗的系统函数系数不同的采样疏密程度将得到不同的窗口宽度，如图7所示严格来说，可以采用FIR等方法，直接设计滤波器采用有限冲击响应数字滤波器（FIR）设计步骤如图8所示，这个系统函数在处衰减为-5dB，通带起伏1dB，阻带衰减-15dB，在5阶滤波器中具有较好的低通特性实质上，H1D=120[-23993-2]正是前面讨论到的对Sa函数的采样，如图7所示3算法评估与结果

3.1测试材料本文通过CIF格式与QCIF格式之间的分辨率转换来测试算法CIF与QCIF格式文件（扩展名“.yuv”）每一帧图像包含3个矩阵每1个像素点由8b整数（unsignedchar）记录CIF与QCIF之间做分辨率转换时，需要对Y，Cb，Cr矩阵分别做变换矩阵分别做变换对于降分辨率结果的测试方式，除了编译器下直接测试之外，也可以借助YUVViewer等软件例如，在C++下，从测试文件中读出若干帧图像，变换后写入新文件，用YUVViewer等软件查看，可以验证正确性，并做定性分析在做无偏性测试时，将有偏和无偏的输出图像写入连续的两帧，通过YUVViewer等软件播放时，可以明显感受到图像偏移

3.2评估标准通过某种参数来评价算法的结果是否合理可以采用的参数很多，比如差值的绝对值的和，但是不同样本间这个参数不够明显，不宜采用通常采用峰值信噪比来评价算法输出结果峰值性噪比（PeakSignalNoiseRatio，PSNR）定义为PSNR=10×log2552MSE其中分别指被评价图像像素点和参考标准图像像素点；Size指图像尺寸，宽×高算法输出结果与参考标准对比，计算PSNR，用于评估算法本文中的参考标准有两种，由DCT域法降分辨率法生成DCT域法降分辨率是通过DCT变换（离散余弦变换），在变换域内进行的降分变率过程DCT系数的特点是图像经DCT变换后，大部分能量都保留在DCT矩阵的低频带（左上角矩阵），去除高频系数，得到的左上角矩阵保持着图像的大部分信息再经过逆DCT变换，可以得到低分辨率图像参考标准1是全图的二维DCT变换法，如图9所示参考标准1的算法计算量很大，尤其是对整体图像做DCT变换和IDCT（InverseDCT）变换（反DCT变换）的步骤，但是结果是较为理想的，仅用于算法对比，一般不用于实际的视频降分辨率注

①二维DCT变换；

②取低频分量；

③形成仅含有低频成分的DCT系数；

④二维逆DCT变换；

⑤等间隔抽样图9参考标准1产生步骤全图的二维DCT变换法参考标准2是DCT域8×8分块算法，如图10所示实际的降分辨率算法中，为了降低算法复杂度，往往先对图像进行等大小分块，然后对分割得到的小图像分别进行处理参考标准2的算法计算量适中，但是结果较参考标准1稍差，可作为实际的视频降分辨率算法本文将其用于对像素域平均法、滤波与子采样法的评估这种分块的算法有时的确能大幅降低算法复杂度（尤其是DCT变换），然而它并不适用于全部算法，而且即使适用，分块也会引起的图像质量变差，这种由分块处理引起的图像质量下降称作块效应使用时，要综合考虑计算复杂度和块效应，将图像质量的下降控制在合理、可接受的范围内高阶的滤波子采样法不宜应用分块方法，因为会有大量的滤波抽头点在分块以外

3.3像素域单帧对比测试参加对比的算法像素平均值法（1阶滤波子采样法）H=

[11]；（5阶）滤波子采样法H=[1-52020-51]；（5阶）过滤子采样法H=[-23993-2]测试图像3个不同文件的第一帧，CIF格式参考标准1为QCIF文件Ref_DCT.yuv；参考标准2为QCIF文件Ref_DCT

8.yuv从表1可看出，宽矩形窗H=[1-52020-51]得到的降分辨率图像的PSNR并不好，甚至不如1阶滤波子采样原因是变换分辨率后保留了过多的高频分量而本文提出的H=[-23993-2]20表现出优异的PSNR

3.4块效应的影响整块滤波采样对比8×8分块滤波采样4个yuv输入文件（Bridge，Foreman，Mobile，Motherdaughter），分别连续测试前300帧测试1以整块DCT变换为参考标准，比较PSNR；测试2以8×8分块DCT变换为参考标准，比较PSNR测试数据如图11所示现在计算图11中整块与分块结果的PSNR的差值的绝对值的最大值Max|ΔPSNR|，如表2所示可以看出，在8×8分块中，PSNR的变化并不多，在合理范围之内，说明了这个5阶滤波器在分块算法中的适应性

3.5综合对比测试测试对象为CIF图像“Mobile.yuv”，一共300帧要测试的算法包括对整个图像的像素平均值法H=

[11]（8×8分块后的像素平均值法与之相同）；对整个图像的滤波子采样法H=[-23993-2]；8×8分块后的滤波子采样法H=[-23993-2]参考标准有2个整个图像的DCT变换法；8×8分块DCT变换法综合对比测试的结果见图12分析如下首先像素平均值法在不同参考标准下质量表现都是最差，其次结果也体现了算法分块方式相同则结果接近的特点本文提出的5阶滤波H1D=120[-23993-2]子采样法复杂度适中，结果令人满意此滤波方法在8×8分块中表现依然良好4结语本文基于视频转码的需要，设计并实现了一种分辨率转换的算法，着重研究了CIF格式到QCIF格式的分辨率转换，并在编译器上进行仿真验证与评估特别的，本文设计提出了一个性能良好的5阶滤波èH1D=120[-23993-2]子采样降分变率方法，对分辨率转换相关软硬件设计有指导意义参考文献

[1]RICHARDSONIE.TheH.264advancedvideocompressionstandard[M].Hoboken，NewJersey JohnWileySons，

2011.

[2]RAOKR，KIMDN，HWANGJJH.264/MPEG4advancedvideocodinginvideocodingstandards[M].Berlin，Germany SpringerNetherlands.

2014.

[3]GAOW，MAS.videocodingbasicprincipleinadvancedvideocodingsystems[M].Berlin，Germany SpringerInternationalPublishing，

2014.

[4]王峥，郑江滨.基于DCT域的任意分辨率视频可伸缩编码[J].微电子学与计算机，2012

（4）

4548.

[5]周晋.视频编码中自适应插值滤波器的算法研究[D].西安西安电子科技大学，

2011.

[6]石巍，胡海兵，邢岩.基于多分辨率控制算法的中频有源电力滤波器[J].中国电机工程学报，2011

（3）

4751.作者简介袁莉（1961—），女，教授研究方向为计算机应用。