劉 婧，王 威，李 驥，楊蔚蔚

（長沙理工大學計算機與通信學院，湖南 長沙 410114）

基于對偶樹復小波變換的模糊圖像質(zhì)量評價*

劉 婧，王 威，李 驥，楊蔚蔚

（長沙理工大學計算機與通信學院，湖南 長沙 410114）

小波域和結構相似度SSIM的質(zhì)量評價方法已經(jīng)成為圖像處理領域的研究熱點，然而都存在一定的缺陷：傳統(tǒng)的離散小波變換缺乏平移不變性，其方向選擇性也十分有限；對于嚴重模糊的圖像，SSIM評價結果并不十分準確。基于此，提出了一種適應于模糊圖像質(zhì)量評價的新算法。該算法用對偶樹復小波變換DTCWT將圖像進行分解來獲取復小波系數(shù)，然后對所得到的六個方向的高頻子帶系數(shù)分別進行平均梯度幅度值的結構相似度 MGSIM測量，最后將所得到的全部 MGSIM的均值作為最終的原始模糊圖像的模糊值。仿真實驗驗證了本方法比結構相似度更吻合人眼的視覺效果，與主觀評價方法具有很好的一致性，并且在各方面的性能都優(yōu)于目前有關文獻的方法。

小波域；對偶樹復小波變換；SSIM；平均梯度幅度值的結構相似度（MGSIM）

1 引言

在圖像處理過程中，如何對圖像處理結果進行質(zhì)量評價是一個值得重點關注的問題。模糊性是圖像最常見的失真類型之一，導致圖像失真的原因有多種，如光學系統(tǒng)的像差、相機與物體間的相對機械運動等。由于模糊圖像的廣泛存在，所以對模糊圖像質(zhì)量評價算法的研究具有非常重要的意義。當前經(jīng)常使用的評價方法可以分為主客觀這兩種常見的類型。其中可信度最高的評價方法是數(shù)字圖像的主觀評價DMOS（Different Mean Objective Score）方法［1］，主觀評價方法通常會進行大量重復實驗，不僅耗時多、費用高，并且無法建立對應的數(shù)學模型，脫離人工不能得到結果，因此不能應用到實際的視頻和圖像處理系統(tǒng)當中［2］。所以，圖像質(zhì)量的客觀評價方法目前已成為研究的熱點［3］。依據(jù)對原始圖像的參照程度，客觀圖像質(zhì)量評價可分為三類：全參考質(zhì)量評價、半（弱）參考質(zhì)量評價以及無參考質(zhì)量評價［4］。在許多環(huán)境下，參考圖像無法獲取，因此需要深入研究無參考質(zhì)量評價方法。

Feichtenhofer C H等人［5］根據(jù)人類視覺系統(tǒng)特性提取模糊圖像邊緣的梯度特征來建立 PSI （Perceptual Sharpness Index）指標，通過PSI指標直接計算圖像質(zhì)量；Ferzli R等人［6］采用恰可察覺模型JNB（Just Noticeable Blur）對模糊圖像進行質(zhì)量評價；Marziliano P等人［7］利用圖像中的邊緣平均寬度作為圖像模糊測度，方法簡單易行，可以取得更好的評估效果，Moorthy K A等人［8］提出先用分類算法求出失真概率，再單獨計算失真強度并加權和等。以上都是從一定角度出發(fā)，構造了模糊圖像質(zhì)量評價指標，取得了一定的成果。但是，評價結果存在一定的缺陷，還需要進一步提高和改善，例如文獻［5］的方法要求一次性先定義大批參數(shù)來進行失真概率的學習。

結合上述分析，本文提出了一種基于對偶樹復小波變換和平均梯度幅度值結構相似度 MGSIM （Mean Gradient magnitude based structure SIM-milarity）相結合的模糊圖像質(zhì)量評估方法，先將圖像進行對偶樹復小波變換分解，針對六個方向的高頻子帶系數(shù)進行平均梯度幅度值的結構相似度測量，然后將所得到的全部 MGSIM的均值作為最終的原模糊圖像的模糊值。

2 結構相似度計算

2.1 由模糊圖像得到參考圖像

原始無失真圖像在經(jīng)過高斯低通濾波器后可以得到一幅模糊圖像，所得到的圖像模糊程度是由低通濾波器的帶寬來決定的。經(jīng)過對原始無失真圖像、模糊圖像和濾波器使圖像模糊的研究，發(fā)現(xiàn)人類對模糊值變化較小的圖像之間的差異難以察覺。

自然圖像a在經(jīng)過11×11并且σ=1.5的高斯低通濾波器以后得到模糊圖像b，同樣的原理，把模糊圖像b也經(jīng)過相同高斯濾波器之后可以得到再模糊圖像c。全參考圖像質(zhì)量評價中，必須要有清晰圖像a作為參考，才能去評價模糊圖像b的好壞。根據(jù)前面分析，如果沒有清晰圖像作為參考圖像，那么可以用再模糊圖像c作為模糊圖像b的參考圖像。因此，對模糊圖像b的無參考圖像質(zhì)量評價問題，就轉化為以模糊圖像c作為參考圖像的全參考圖像質(zhì)量評價問題。與模糊圖像相比，清晰的圖像包含有更詳細的信息，即高頻分量，所以可以通過高斯低通濾波器構造模糊圖像c。

2.2 SSIM計算

圖像與圖像之間的像素差異即圖像的高度結構性，用人眼直接觀察時容易被忽略?？紤]到多方面問題Zhou Wang等人［9］提出了一種新的結構相似度SSIM（Structural SIMilarity）圖像質(zhì)量評價方法。

SSIM包含了自然圖像和失真圖像分別在亮度、對比度以及結構信息各方面的對比，這三部分均可由單獨的函數(shù)來表示，通常使用的計算方法如下：

綜合以上計算結果，可以得到模型的評判方法［10］：

其中，x為自然圖像，y為質(zhì)降圖像，ux和uy為自然圖像與質(zhì)降圖像中局部像素塊的平均亮度值，σx和σy則表示為兩幅圖像中局部像素塊的量度標準差，σxy為兩幅圖像中對應像素塊的亮度相關系數(shù)。參數(shù)α、β、γ均大于零，以調(diào)整三個部分所占的比重，一般都取1；C1、C2、C3用于增加計算結構的穩(wěn)定性，C1、C2、C3為較小的值。MSSIM（Mean Structural SIMilarity）算法將一幅圖像劃分成M個圖像塊，每一個圖像子塊的大小統(tǒng)一為N×N。由此整幅圖像的MSSIM定義為：

2.3 MGSIM計算

雖然SSIM可以用來作為一種檢測圖像質(zhì)量好壞的指標，但是對于嚴重模糊圖像評估結果SSIM可靠 性 不 高［11，12］。人 類 對 圖 像 的 邊緣 信 息部分極度敏感，梯度能夠很好地顯示對比度和細節(jié)等不易察覺的圖像紋理特征的變化，能夠更好地評估圖像的質(zhì)量好壞，因此可以把梯度看成圖像的重要結構信息，進行基于梯度的結構失真的圖像質(zhì)量評價方法［13］。

梯度幅值可以采用以下公式計算：

在圖像質(zhì)量評價中我們研究的是數(shù)字圖像，對于數(shù)字圖像一般可以用差分運算代替微分運算。公式（6）可借用差分公式來近似代替：

計算梯度的方法有很多，可用八鄰域絕對值和來代表圖像的梯度幅度值。圖像x和y的梯度相似度g為：

其中Gx（i，j）、Gy（i，j）分別表示圖像x和圖像y在 （i，j）點的梯度幅度值。

將式（4）中的s（x，y）用梯度相似度g（x，y）代替，就可以得到GSIM。那么整個圖像的梯度結構相似度可以由每一個子圖塊的GSIM的均值來求得：

設兩幅圖像x、y，xi、yi表示圖像的子塊，i∈［1，M］，整個圖像分塊數(shù)為M，則：

MGSIM著重于圖像之間細節(jié)上的相似度，因此更符合人眼視覺系統(tǒng)。

3 結合MGSIM和DTCWT無參考模糊圖像評價算法

自離散小波變換DWT（Discrete Wavelet Transform）出現(xiàn)以來，它作為信號和圖像處理的工具已經(jīng)得到了一定的肯定，但是傳統(tǒng)的離散小波變換存在缺少平移不變性和有限的方向選擇性兩個方面的弊端。

為了進一步克服離散小波變換所帶來的弊端，Kingsbury N G在1998年提出對偶樹復小波變換DTCWT（（Dual Tree Complex Wavelet Transform）［14］這一新的小波變換，它是以離散小波變換為基礎的一種新型小波變換模型，是二叉樹算法的一種，不僅保留了DWT和CWT（Complex Wavelet Transform）的優(yōu)點，而且還解決了其他小波變換的缺點，并且還具有更多比較好的優(yōu)點，它有平移不變性和多方向性。

對偶樹復小波變換是兩個實值離散小波變換的和，它會得到±15°、±45°和±75°六個方向的子帶，但是傳統(tǒng)的離散小波變換只能夠產(chǎn)生三個方向子帶，分別是0°、45°和90°。一維的對偶樹復小波變換分解公式如下：

其中，ζ表示一維向量，w（ζ）是尺度函數(shù)，φ（ζ）是擴散母小波，N是自然數(shù)集合，j和l分別是平移和擴散指數(shù)，Sj0，l是尺度系數(shù)，cj，l為復小波系數(shù)。

與一維對偶樹復小波變換類似，二維對偶樹復小波變換把二維圖像Zα，β進行分解，分解如下：

Figure 1 Principle diagram of the dual tree complex wavelet transform圖1 對偶樹復小波變換原理圖

其中，θ∈Θ=｛±15°，±45°，±75°｝是復小波函數(shù)的方向，這樣就可以把Zα，β分解成一個復數(shù)的低通子帶和六個復數(shù)值高通子帶。

圖1是DTCWT的詳細圖［15］，該圖給出了兩個互相平行的小波樹，即樹A和樹B。濾波器H0b和H1b與濾波器H0a和H1a之間的延遲恰好是一個采樣值的間隔，那么就可以確保Tree B中的第一層的向下采樣取到Tree B中因隔點采樣運算而舍棄的不能保留的采樣值。

在實現(xiàn)對偶樹復小波變換的正變換后，還要實現(xiàn)其逆變換。為了實現(xiàn)其逆變換，我們在每棵樹中都使用了雙正交濾波器，這也是因為該濾波器的完全重構特性，同時為了達到整個系統(tǒng)近似的平移不變性，必須對最后兩棵樹所得到的結果求平均。經(jīng)過二維對偶樹復小波變換后，得到一個低頻子帶和六個方向的高頻子帶，方向分別是±15°、±45°和±75°，這么多的方向可以更好地描述圖像紋理、邊界等特征，可更好地用于圖像后期處理。

4 實驗結果與分析

4.1 實驗流程

基于上述分析，模糊指標計算的具體步驟如下：

（1）為原圖像構造參考圖像。將原始模糊圖像經(jīng)過大小為11×11，σ=1.5的高斯低通濾波器得到再次模糊圖像，作為參考圖像。經(jīng)過這一步就把原有的無參考圖像質(zhì)量評價轉換成了全參考圖像質(zhì)量評價。

（2）提取對偶樹復小波系數(shù)。將原始模糊圖像和再模糊圖像分別進行對偶樹復小波變換得到± 15°、±45°和±75°的六個方向的子頻帶。

（3）MGSIM計算。把兩幅圖像得到的第i個系數(shù)分別進行平均梯度幅度值的結構相似度計算MGSIM（xi，yi）。

（4）對每一個頻帶的MGSIM求取平均值，獲得整幅圖像的平均梯度幅度值的結構相似度。最后將所得到的全部MGSIM的平均值等同于最后的評價圖像的模糊值blur：

其中，N表示子頻帶個數(shù)；MGSIM表示每個子頻帶的平均梯度幅度值的結構相似度。MGSIM越大，所求的blur值越大，表示原始圖模糊圖像的質(zhì)量越差，反之則質(zhì)量越好。

本文方法的流程圖如圖2所示。

Figure 2 System flowchart圖2 系統(tǒng)流程圖

4.2 實驗結果分析

實驗選用LIVE實驗室提供的blur圖片庫來驗證本文方法的有效性（軟硬件平臺：PC；OS：Windows XP；CPU：2.71 GHz；RAM：2 GB；Software：Matlab 7.0）。選用兩個客觀參量作為評價指標，即經(jīng)過非線性回歸之后的皮爾遜相關系數(shù)CC（Correlation Coefficent）和斯皮爾曼等級相關系數(shù)SROCC（Spearman Rank Order Correlation Coefficient）。SROCC與CC是0到1的值，它們值越接近1，代表其性能越好。

為了對小波變換與對偶樹復小波變換做深入比較，用Lena圖進行了實驗。圖3為小波變換結果，圖3a～圖3c分別是90°、0°和高頻結果圖。圖4為對偶樹復小波變換結果圖，圖4a～圖4f分別表示在方向—15°、—45°、—75°、75°、45°、15°的圖像。對比兩個變換后的結果，我們可以看出具有更多方向圖像的對偶樹復小波變換可以更好地描述圖像在各個方向上的紋理特征。

圖5a表示原圖像，圖5b～圖5g分別是高頻上方向為—15°、—45°、—75°、75°、45°、15°的變換結果。其中圖5h是低頻圖像，因為低頻圖像很容易受到噪聲、光照等因素的干擾與影響，所以它一般不用于圖像質(zhì)量評價中。

為了進一步評價本文方法的性能，本文算法與目前已有算法［5～7］的對比如表1所示。

Figure 3 Results after the wavelet transform圖3 小波變換后結果

Figure 4 Results 1 after the dual tree complex wavelet transform圖4 對偶樹復小波變換后結果1

Table 1 Per formance comparison of different evaluation algorithms表1 不同評價算法的性能比較

由表1可見，與JNB算法、分類算法中的方法比較，本文算法的結果和主觀評價方法的一致性、準確性都有比較明顯的提升，相對于分類算法在CC評價值上提升了將近2個百分點，在SROCC上提升了5個百分點。把CC和SROCC這兩個重要質(zhì)量衡量指標結合起來思考，發(fā)現(xiàn)本文算法在一定程度上能夠獲得比較好的性能。

總之，由表1中的比較可以看出，本文中提出的方法預測的圖像客觀質(zhì)量評價與圖像主觀評價更加一致，更加符合人眼視覺系統(tǒng)，而且計算復雜度低。因此，本文所提出的新方法可以用于對無參考模糊圖像進行質(zhì)量評價。

Figure 5 Results 2 after the dual tree complex wavelet transform圖5 對偶樹復小波變換后結果2

5 結束語

針對模糊失真圖像，本文提出了一種基于對偶樹復小波變換和平均梯度幅度值結構相似度相結合的模糊圖像質(zhì)量評估方法，通過高斯低通濾波將原始圖像再模糊，對原始圖像和模糊圖像經(jīng)過對偶樹復小波變換和MGSIM等步驟，最終得到原始圖像模糊值。實驗結果表明，本方法是一種較準確可靠的無參考模糊圖像評價方法。怎樣把這種構造參考圖像的評價方法應用到其他各種失真類型的無參考圖像質(zhì)量評價中將是以后工作的研究重心。

［1］ Sheikh H R，Sabir M F，Bovik A C.A statistical evaluation of recent full reference in image quality assessment algorithms［J］.IEEE Transactions on Image Processing，2006，15（11）：3440-3451.

［2］ Wang Hong-jiao.Research and implementation of the image quality assessment method［D］.Xi'an：Xi'an University of E-lectronic Science and Techology，2014.（in Chinese）

［3］ Lu W，Gao X B，Li X L，et al An image quality assessment metric based Contourlet［C］∥Proc of IEEE International Conference on Image Processing，2008：1172-1175.

［4］ Wang Zheng-you，Li Zhen-xing，Lin Wei-si，et al.Improved image quality assessment model incorporating H VS and FSIM［J］.Chinese Journal of Scientific，2013，33（7）：1606-1612.（in Chinese）

［5］ Feichtenhofer C H，Schallauer F H.A perceptual sharpness metric based on local edge gradient analysis［J］.IEEE Signal Processing Letters，2013，20（4）：379-382.

［6］ Ferzli R，Karam L J.A no-reference object image sharpness metric based on just-noticeable blur［J］.IEEE Transactions on Image Processing，2009，18（4）：717-728.

［7］ Marziliano P，Dufaux F.Winkler S，et al.Perceptual blur and ringing metrics：Application to JPEG 2000［J］.Signal Processing：Image Communication，2004，19（2）：163-172.

［8］ Moorthy K A，Bovik A C.A two-step framework for constructing blind image quality indices［J］.IEEE Signal Processing Letters，2010，17（5）：513-516.

［9］ Zhou Wang，Conrad Bovik.Image quality assessment：From error visibility to structural similarity［J］.IEEE Transactions on Image Processing，2004，13（4）：600-612.

［10］ Li Ling-ling，Wang Hui，Yan Cheng-xin.Study on regional segmentation image fusion technique based on improved structure similarity［J］.Application Research of Computers，2014（6）：1914-1917.（in Chinese）

［11］ Liu An-min，Lin Wei-si，Narwaria Manish.Image quality assessment based on gradient similarity［J］.IEEE Transactions on Image Processing，2012，47（17）：962-963.

［12］ Gu Shan-bo，Shao Feng，Jiang Gang-yi，et al.An objective quality assessment metric for stereoscopic images based on just noticeable distortion［J］.Journal of Optoelectronics Laser，2012，23（5）：999-1004.（in Chinese）

［13］ Yang Chun-ling，Kuang Kai-zhi，et al.Gradient based structural similarity for image quality assessment［J］.Journal of South University of Technology，2006，34（9）：22-25.（in Chinese）

［14］ Kingsbury N G.Image processing with complex wavelets ［J］.Philosophical Transactions：Mathematical，Physical and Engineering Science（S1364-503x），1999，357（1760）：2543-2560.

［15］ Kingsbury N.Image processing with complex wavelets［EB/ OL］.［2000-11-08］.www-sigproc.eng.cam.ac.uk/publication/ngk/ngk99a.ps.gz.

附中文參考文獻：

［2］ 王洪嬌.圖像質(zhì)量評價方法的研究及實現(xiàn)［D］.西安：西安電子科技大學，2014.

［4］ 王正友，李振興，林維斯，等.結合 HVS和相似特征的圖像質(zhì)量評價方法［J］.儀器儀表學報，2013，33，（7）：1606-1612.

［10］ 李玲玲，王蕙，閆成新.基于改進結構相似度的區(qū)域劃分圖像融合技術研究［J］，計算機應用研究，2014（6）：1914-1917.

［12］ 顧珊波，邵楓，蔣剛毅，等.失真的立體圖像質(zhì)量客觀評價方法［J］.光電子·激光，2012，23（5）：999-1004.

［13］ 楊春玲，曠開智，陳冠豪，等.基于梯度的結構相似度的圖像質(zhì)量評價［J］.華南理工大學學報，2006，34（9）：22-25.

劉婧（1991 ），女，湖南益陽人，碩士生，研究方向為智能信息處理。E-mail：123410932@qq.com

LIU Jing，born in 1991，MS candidate，her research interest includes intelligent signal processing.

王威（1974），男，山東青島人，博士，教授，研究方向為智能信息處理。E-mail：wangwei@csust.edu.cn

WANG Wei，born in 1974，PhD，professor，his research interest includes intelligent signal processing.

李驥（1981），男，湖南婁底人，博士，講師，研究方向為信號處理和電磁散射。E-mail：hangliji@163.com

LI Ji，born in 1981，PhD，lecturer，his research interests include signal processing，and electromagnetic scattering.

楊蔚蔚（1988 ），女，湖南常德人，碩士生，研究方向為智能信息處理。E-mail：1732648392@qq.com

YANG Wei-wei，born in 1988，MS candidate，her research interest includes intelligent signal processing.

Blur image quality assessment based on DTCWT

LIU Jing，WANG Wei，LI Ji，YANG Wei-wei
（School of Computer and Communication Engineering，Changsha University of Science&Technology，Changsha 410114，China）

Wavelet domain and structural similarity（SSIM）quality assessment method have become hotspots in the field of image processing，however，both of them have some flaws：the traditional discrete wavelet transform lacks of translational invariance and its direction selectivity is also highly limited；for severe blurred images，the results of the SSIM are not very accurate.Therefore，we propose a new algorithm for blur image quality evaluation.This algorithm uses dual tree complex wavelet transform（DTCWT）image decomposition to obtain the complex wavelet coefficients and the high frequency sub band coefficients of all the six directions.Then the structural similarity of the average gradient amplitude is measured.Finally all the mean gradient-magnitude-based structural similarity（MGSIM）average is calculated as the final fuzzy values of the original blur images.Experimental results show that the proposed method fits the visual characteristics of human eyes better in contrast with the structural similarity method，and well matches the results of subjective evaluation methods.The assessment results are better than the current literature in terms of overall performance.

wavelet domain；dual-tree complex wavelet transform；SSIM；mean gradient magnitude based structural similarity

TP391.4

A

10.3969/j.issn.1007-130X.2015.08.023

1007-130X（2015）08-1573-06

2014-09-05；

2014-12-16

博士后基金資助項目（2013M542467）；國防973基金資助項目

通信地址：410114湖南省長沙市長沙理工大學計算機與通信學院

Address：School of Computer and Communication Engineering，Changsha University of Science&Technology，Changsha 410114，Hunan，P.R.China