數(shù)字視頻降分辨率算法的濾波器設(shè)計(jì)

1、    數(shù)字視頻降分辨率算法的濾波器設(shè)計(jì)    袁莉+初秀娟+王志國(guó)摘 要： 數(shù)字視頻在不通終端上應(yīng)用時(shí)往往需要將分辨率降低到原始分辨率的1/2或1/4。因而需要有合適的算法改善降分辨率后圖像的清晰度，滿足用戶對(duì)顯示圖像質(zhì)量的要求?；谝曨l轉(zhuǎn)碼的需要，針對(duì)分辨率轉(zhuǎn)換算法中最常用的濾波與子采樣法，設(shè)計(jì)了一種用于濾波與子采樣法的5階數(shù)字濾波器，并在編譯器上進(jìn)行仿真驗(yàn)證與評(píng)估。仿真結(jié)果顯示，新濾波器表現(xiàn)出優(yōu)異的峰值信噪比，并且在分塊算法中依然優(yōu)異?？傊?，這里提出的方法不但復(fù)雜度適中，適應(yīng)性強(qiáng)，而且能夠高效、高質(zhì)量地完成數(shù)字視頻降分辨率。關(guān)鍵詞： 視頻轉(zhuǎn)碼； 分

2、辨率； 濾波與子采樣法； 濾波器： tn913?34 ： a ： 1004?373x（2015）20?0026?05design of filter for resolution reduction algorithm used in digital videoyuan li， chu xiujuan， wang zhiguo（suihua university， suihua 152061， china）abstract： when the digital video is used in different terminals， the resolution is often reduce

3、d by a factor of 2 or 4， so a suitable algorithm is needed to improve the image sharpness after resolution reduction to meet user demands for image quality. in view of the need of video transcoding， a fifth?order digital filter used for filtering and sub?sampling method was designed， which was asses

4、sed and verified on interpreter. the simulation results show that the new filter has high psnr （peak signal to noise ratio）， and is still excellent in the block algorithm. in short， the proposed method has moderate complexity and strong adaptability， and can achieve high efficiency and quality digit

5、al video resolution reduction.keywords： video transcoding； resolution； filter and sample； filter各種數(shù)字視頻實(shí)際應(yīng)用在各種不同的終端設(shè)備上，如電腦、平板電腦、電視、手機(jī)等，他們支持不同分辨率和不同壓縮標(biāo)準(zhǔn)的視頻，因此經(jīng)常存在對(duì)同一視頻進(jìn)行不同分辨率間轉(zhuǎn)換的需求。例如使用手持終端設(shè)備接收觀看運(yùn)營(yíng)商提供的標(biāo)準(zhǔn)清晰度數(shù)字電視節(jié)目，往往需要將節(jié)目分辨率降低到原始分辨率的12或14。因而需要有合適的算法改善降分辨率后圖像的清晰度，滿足人們對(duì)顯示圖像質(zhì)量的要求。對(duì)于這樣的實(shí)際需求，濾波子采樣算法降分辨率的算法和

6、物理結(jié)構(gòu)實(shí)現(xiàn)相對(duì)簡(jiǎn)單，適合在軟硬件系統(tǒng)中集成。但是此方法的濾波器設(shè)計(jì)尤為關(guān)鍵，需要大量的數(shù)學(xué)和信號(hào)系統(tǒng)相關(guān)的建模和演算。設(shè)計(jì)一個(gè)性能良好的濾波子采樣降分變率方法，并且進(jìn)行全面系統(tǒng)地測(cè)試，對(duì)分辨率轉(zhuǎn)換相關(guān)軟硬件設(shè)計(jì)有指導(dǎo)意義。另一方面，國(guó)內(nèi)外大部分算法都是在像素域分塊后在變換域（如avs域、dct域等）進(jìn)行，保證了速度但因?yàn)榉謮K操作引入了峰值信噪比（psnr）的損失（稱之為塊效應(yīng)）。對(duì)塊效應(yīng)的分析評(píng)估可以在轉(zhuǎn)換速度與精度中找到平衡點(diǎn)，可以有效地指導(dǎo)以及優(yōu)化各種算法的設(shè)計(jì)和實(shí)現(xiàn)，為圖像領(lǐng)域未來(lái)的各種算法設(shè)計(jì)、程序設(shè)計(jì)和硬件設(shè)計(jì)都有深遠(yuǎn)影響。本文基于視頻轉(zhuǎn)碼的需要，設(shè)計(jì)并實(shí)現(xiàn)了一種應(yīng)用于此方法的5階

7、數(shù)字濾波器，并在編譯器上進(jìn)行仿真驗(yàn)證與評(píng)估，取得了令人滿意的結(jié)果。此濾波器在像素域8×8分塊后使用，表現(xiàn)依然良好。1 濾波與子采樣法濾波與子采樣法降分辨率基本過(guò)程如圖1所示。將圖像x的行與列看做數(shù)字序列，用一個(gè)系統(tǒng)函數(shù)對(duì)其濾波，濾去高頻分量，得到圖像y，然后等間隔取出像素點(diǎn)組成新的圖像。這個(gè)過(guò)程中濾波后得到的m×n像素圖像（見(jiàn)圖1中）有34的像素是無(wú)用的，因?yàn)椴蓸又蝗?4的像素。所以只要確定了采樣點(diǎn)，則只需對(duì)采樣點(diǎn)進(jìn)行濾波即可，可以省去很多計(jì)算，過(guò)程如圖2所示。注：原始m×n像素圖像；濾波后的m×n像素圖像；等間隔采樣；合并成新的m2×n2圖像

8、圖1 濾波與子采樣法降分辨率（一）注：原始m×n像素圖像；濾波與等間隔采樣同時(shí)進(jìn)行；合并成新的m2×n2像素圖像圖2 濾波與子采樣法降分辨率（二）綜上所述，濾波與子采樣降分變率方法的第1步為濾波，對(duì)于k階的濾波與子采樣降分變率方法，所用到的低通濾波器系統(tǒng)的方程為：y（n）=h0x（n）+h1x（n-1） + h2x（n-2） + +hkx（n-k）濾波與子采樣降分變率方法的第2步子采樣，及對(duì)變換后的y等間隔采樣。對(duì)于m×n到m2×n2的降分辨率，則要求每2行取出1行，每2列取出1列，組成新的圖像。行與列可以分開(kāi)考慮成一維的序列，如果濾波器系統(tǒng)是y（n）=

9、h0x（n）+h1x（n-1）+h2x（n-2）的低通濾波器，輸入序列x=x1，x2，xm，則濾波后序列是：y=h2+h1+h0x1h2+h1x1+h0x2h2x1+h1x2+h0x3+h2xm-2+h1xm-1+h0xm。等間隔采樣，結(jié)果可以是如下兩種： y= h2+h1+h0x1 h2x1+h1x2+h0x3 h2x3+h1x4+h0x5；y=h2+h1x1+h0x2 h2x2+h1x3+h0x4 h2x4+h1x5+h01x6以1階濾波子采樣為例，像素域平均法是最簡(jiǎn)單形式。其實(shí)現(xiàn)為，每2×2矩陣的4個(gè)元素取均值，得到1個(gè)新的像素值，代替原來(lái)的4個(gè)像素。本質(zhì)是對(duì)行和列同時(shí)做1階低

10、通濾波并采樣。在一個(gè)維度上，此低通濾波器的系統(tǒng)函數(shù)是：h（z）=12+12z，低通濾波系統(tǒng)的差分方程是：y（n）=12x（n）+12x（n-1）。長(zhǎng)度為n的序列，濾波后得到的新序列長(zhǎng)度仍為n，對(duì)其進(jìn)行每2像素點(diǎn)取一點(diǎn)采樣，得到的12n像素序列即為降分變率后的序列。對(duì)于二維圖像，把系統(tǒng)函數(shù)h（z）=12+12z寫作h1d=1212，則二維的系統(tǒng)函數(shù)為：h=h1d×h1d=14141414。在實(shí)際實(shí)現(xiàn)中，其實(shí)并不需要對(duì)每一個(gè)像素點(diǎn)進(jìn)行濾波得到m×n像素圖像再采樣，濾波只要對(duì)采樣點(diǎn)進(jìn)行即可。所以，像素域平均法降分辨率具體實(shí)現(xiàn)就是每2×2矩陣的4個(gè)元素取均值，得到1個(gè)新的

11、像素值，代替原來(lái)的4個(gè)像素。以2階濾波子采樣為例。對(duì)連續(xù)的一行中的3個(gè)點(diǎn)，或?qū)B續(xù)的3行進(jìn)行變換：y（n）=h0x（n）+h1x（n-1）+h2x（n-2）。此濾波器系統(tǒng)函數(shù)：h（z）=h0+h11z+h21z2或?qū)懽鱤1d=h0h1。例如，取h1d=121412，此系統(tǒng)的頻率響應(yīng)如圖3所示，為一低通數(shù)字濾波器。二維具體實(shí)現(xiàn)，對(duì)每3×3個(gè)像素乘以如下矩陣得到一個(gè)新的像素點(diǎn)，即可同時(shí)完成對(duì)列和行的濾波，此時(shí)：h2d=h0h0h1h0h2h0h0h1h2h1h2h1h0h2h1h2h2h2如果：h1d=121412可知：h2d=1418141811618141814圖3 h1d=12 1

12、4 12的低通濾波器2 濾波器設(shè)計(jì)2.1 設(shè)計(jì)濾波器階數(shù)簡(jiǎn)化考慮，先低通濾波再采樣，先只考慮空間一維，因?yàn)槎S只是2個(gè)一維變換的簡(jiǎn)單結(jié)合。首先討論不同階數(shù)濾波后的最佳采樣結(jié)果，以14階為例，這些階數(shù)的最小偏移情形如圖4所示。圖4 不同階數(shù)無(wú)偏性最好的濾波采樣方式因此，濾波之后必須選擇合理的采樣方法，盡可能減小圖像偏移。對(duì)于折半降分辨率，最好使用奇數(shù)階的低通濾波器來(lái)濾波。偶數(shù)階濾波后無(wú)論如何采樣，都會(huì)引起圖像中信息的偏移。舉一個(gè)簡(jiǎn)單實(shí)例，圖5中中背景是模擬圖像，·表示原采樣點(diǎn)（即數(shù)字圖像像素點(diǎn)），×表示新采樣點(diǎn)（降分變率后圖像像素點(diǎn)）。方框表示降分變率數(shù)字圖像一個(gè)采樣點(diǎn)（紅色

13、×）所涵蓋的圖片信息。圖5中為1階（或2抽頭，或像素平均值法）濾波，采樣點(diǎn)無(wú)偏；圖5中為2階（或3抽頭），采樣點(diǎn)偏移；圖5中為3階（或4抽頭），采樣點(diǎn)無(wú)偏；圖5中為4階（或5抽頭），采樣點(diǎn)偏移?？梢钥闯?，對(duì)于二維情形，依然是奇數(shù)階的低通濾波器偏移較小?？紤]無(wú)偏性和算法復(fù)雜度，5階濾波（6抽頭）：y（n）=h0x（n）+h1x（n-1）+h2x（n-2）+h3x（n-3）+h4x（n-4）+h5x（n-5），無(wú)偏移，濾波器特性較為理想。圖5 不同濾波采樣法，模擬圖像-數(shù)字圖像-降分變率數(shù)字圖像關(guān)系2.2 設(shè)計(jì)濾波器系數(shù)理想的數(shù)字低通濾波器幅頻特性為矩形窗，如圖6所示。圖6 理想的數(shù)字低

14、通濾波器幅頻特性根據(jù)傅里葉變換，空域sa函數(shù)，則頻域?yàn)榫匦未?。?duì)sa函數(shù)采樣，則可得到矩形窗的系統(tǒng)函數(shù)系數(shù)。不同的采樣疏密程度將得到不同的窗口寬度，如圖7所示。嚴(yán)格來(lái)說(shuō)，可以采用fir等方法，直接設(shè)計(jì)濾波器。采用有限沖擊響應(yīng)數(shù)字濾波器（fir）設(shè)計(jì)步驟：（1） 確定濾波器長(zhǎng)度 n=6，群時(shí)延=12n-1=2.5；（2） 取截止頻率c=2；（3） 取矩形窗函數(shù)，過(guò)渡帶寬度0.9×2n，則由6=n=0.9×2通帶邊-阻帶邊得：通帶<0.35，阻帶>0.65 ；（4） hdn=sin2n-n-，n=r6n是矩形窗函數(shù)，群時(shí)延=2.5，所以hn=hdn×r6n

15、=sin2n-2.5n-2.5，n=0，1，2，5；（5） 取n=0，1，2，.，5 得到hn，系數(shù)歸一化，近似取整，得到：h1d=h0h1h2h3h4h5 =-0.090 0 0.150 1 0.450 2 0.450 2 0.150 1 -0.090 0 -110320920920320-110 =120-23993-2如圖8所示，這個(gè)系統(tǒng)函數(shù)在c=2處衰減為-5 db，通帶起伏1 db，阻帶衰減 <-15 db，在5階濾波器中具有較好的低通特性。實(shí)質(zhì)上，h1d=120-23993-2正是前面討論到的對(duì)sa函數(shù)的采樣，如圖7所示。圖7 采樣疏密程度與窗口寬度的關(guān)系圖8 h1d=120

16、-23993-2的幅頻特性 3 算法評(píng)估與結(jié)果3.1 測(cè)試材料本文通過(guò)cif格式與qcif格式之間的分辨率轉(zhuǎn)換來(lái)測(cè)試算法。cif與qcif格式文件（擴(kuò)展名“.yuv”）每一幀圖像包含3個(gè)矩陣y，cb，cr1。每1個(gè)像素點(diǎn)由8 b整數(shù)（unsigned char）記錄。cif與qcif之間做分辨率轉(zhuǎn)換時(shí)，需要對(duì)y，cb，cr矩陣分別做變換。對(duì)于降分辨率結(jié)果的測(cè)試方式，除了編譯器下直接測(cè)試之外，也可以借助yuv viewer等軟件。例如，在c+下，從測(cè)試文件中讀出若干幀圖像，變換后寫入新文件，用yuv viewer等軟件查看，可以驗(yàn)證正確性，并做定性分析。在做無(wú)偏性測(cè)試時(shí)，將有偏和無(wú)偏的輸出圖像寫

17、入連續(xù)的兩幀，通過(guò)yuv viewer等軟件播放時(shí)，可以明顯感受到圖像偏移。3.2 評(píng)估標(biāo)準(zhǔn)通過(guò)某種參數(shù)來(lái)評(píng)價(jià)算法的結(jié)果是否合理?？梢圆捎玫膮?shù)很多，比如差值的絕對(duì)值的和，但是不同樣本間這個(gè)參數(shù)不夠明顯，不宜采用。通常采用峰值信噪比來(lái)評(píng)價(jià)算法輸出結(jié)果。峰值性噪比（peak signal noise ratio，psnr）定義為psnr=10×log（2552mse）。其中：mse=（pixel-pixelr）2size；pixel，pixelr分別指被評(píng)價(jià)圖像像素點(diǎn)和參考標(biāo)準(zhǔn)圖像像素點(diǎn)；size指圖像尺寸，寬×高。算法輸出結(jié)果與參考標(biāo)準(zhǔn)對(duì)比，計(jì)算psnr，用于評(píng)估算法。本文

18、中的參考標(biāo)準(zhǔn)有兩種，由dct域法降分辨率法生成。dct域法降分辨率是通過(guò)dct變換（離散余弦變換），在變換域內(nèi)進(jìn)行的降分變率過(guò)程。dct系數(shù)的特點(diǎn)是圖像經(jīng)dct變換后，大部分能量都保留在dct矩陣的低頻帶（左上角矩陣），去除高頻系數(shù)，得到的左上角矩陣保持著圖像的大部分信息。再經(jīng)過(guò)逆dct變換，可以得到低分辨率圖像。參考標(biāo)準(zhǔn)1是全圖的二維dct變換法，如圖9所示。參考標(biāo)準(zhǔn)1的算法計(jì)算量很大，尤其是對(duì)整體圖像做dct變換和idct（inverse dct）變換（反dct變換）的步驟，但是結(jié)果是較為理想的，僅用于算法對(duì)比，一般不用于實(shí)際的視頻降分辨率。注：二維dct變換；取低頻分量；形成僅含有低頻成

19、分的dct系數(shù)；二維逆dct變換；等間隔抽樣圖9 參考標(biāo)準(zhǔn)1產(chǎn)生步驟：全圖的二維dct變換法參考標(biāo)準(zhǔn)2是dct域8×8分塊算法，如圖10所示。實(shí)際的降分辨率算法中，為了降低算法復(fù)雜度，往往先對(duì)圖像進(jìn)行等大小分塊，然后對(duì)分割得到的小圖像分別進(jìn)行處理。參考標(biāo)準(zhǔn)2的算法計(jì)算量適中，但是結(jié)果較參考標(biāo)準(zhǔn)1稍差，可作為實(shí)際的視頻降分辨率算法。本文將其用于對(duì)像素域平均法、濾波與子采樣法的評(píng)估。這種分塊的算法有時(shí)的確能大幅降低算法復(fù)雜度（尤其是dct變換），然而它并不適用于全部算法，而且即使適用，分塊也會(huì)引起的圖像質(zhì)量變差，這種由分塊處理引起的圖像質(zhì)量下降稱作塊效應(yīng)。使用時(shí)，要綜合考慮計(jì)算復(fù)雜度和塊

20、效應(yīng)，將圖像質(zhì)量的下降控制在合理、可接受的范圍內(nèi)。高階的濾波子采樣法不宜應(yīng)用分塊方法，因?yàn)闀?huì)有大量的濾波抽頭點(diǎn)在分塊以外。注：分割圖像；分塊二維dct變換；分塊取低頻分量，形成僅含有低頻成分的dct系數(shù)；分塊二維逆dct變換；分塊等間隔抽樣；組合圖像圖10 參考標(biāo)準(zhǔn)2產(chǎn)生步驟：dct域8×8分塊算法3.3 像素域單幀對(duì)比測(cè)試參加對(duì)比的算法：像素平均值法（1階濾波子采樣法）h=1 1；（5階）濾波子采樣法 h=1 -5 20 20 -5 1；（5階）過(guò)濾子采樣法h=-2 3 9 9 3 -2。測(cè)試圖像：3個(gè)不同文件的第一幀，cif格式。參考標(biāo)準(zhǔn)1為qcif文件ref_dct.yuv；參

21、考標(biāo)準(zhǔn)2為qcif文件ref_dct8.yuv。從表1可看出，寬矩形窗h=1 -5 20 20 -5 1得到的降分辨率圖像的psnr并不好，甚至不如1階濾波子采樣。原因是變換分辨率后保留了過(guò)多的高頻分量。而本文提出的h=-2 3 9 9 3 -220表現(xiàn)出優(yōu)異的psnr。表1 像素域單幀對(duì)比測(cè)試3.4 塊效應(yīng)的影響整塊濾波采樣對(duì)比 8×8分塊濾波采樣：4個(gè)yuv輸入文件（bridge，foreman，mobile，mother?daughter），分別連續(xù)測(cè)試前300幀。測(cè)試1：以整塊dct變換為參考標(biāo)準(zhǔn)，比較psnr；測(cè)試2：以8×8分塊dct變換為參考標(biāo)準(zhǔn)，比較psnr

22、。測(cè)試數(shù)據(jù)如圖11所示?，F(xiàn)在計(jì)算圖11中整塊與分塊結(jié)果的psnr的差值的絕對(duì)值的最大值maxpsnr，如表2所示?？梢钥闯?，在8×8分塊中，psnr的變化并不多，在合理范圍之內(nèi)，說(shuō)明了這個(gè)5階濾波器在分塊算法中的適應(yīng)性。3.5 綜合對(duì)比測(cè)試測(cè)試對(duì)象為cif圖像“mobile.yuv”，一共300幀。要測(cè)試的算法包括：對(duì)整個(gè)圖像的像素平均值法h=1 1 （8×8分塊后的像素平均值法與之相同）；對(duì)整個(gè)圖像的濾波子采樣法h=-2 3 9 9 3 -2；8×8分塊后的濾波子采樣法h=-2 3 9 9 3 -2。參考標(biāo)準(zhǔn)有2個(gè)：整個(gè)圖像的dct變換法；8×8分塊dct變換法。綜合對(duì)比測(cè)試的結(jié)果見(jiàn)圖12。表2 maxpsnr圖12 多種算法對(duì)比測(cè)試數(shù)據(jù)分析如下：首先像素平均值法在不同參考標(biāo)準(zhǔn)下質(zhì)量表現(xiàn)都是最差，其次結(jié)果也體現(xiàn)了算法分塊方式相同則結(jié)果接近的特點(diǎn)。本文提出的5階濾波h1d=120-23993