專利名稱：基于三字典塊匹配的圖像復(fù)原方法
技術(shù)領(lǐng)域：
本發(fā)明屬于圖像處理技術(shù)領(lǐng)域，具體地說(shuō)是一種對(duì)模糊圖像進(jìn)行復(fù)原的方法，該方法可用于對(duì)各種已知模糊類型的模糊圖像進(jìn)行復(fù)原。
背景技術(shù)：
圖像復(fù)原是指去除或減輕在獲取數(shù)字圖像過(guò)程中發(fā)生的圖像質(zhì)量下降的現(xiàn)象，它是圖像處理中重要而又富有挑戰(zhàn)性的研究?jī)?nèi)容。對(duì)于圖像復(fù)原問(wèn)題，研究者已經(jīng)提出了很多方法。傳統(tǒng)的復(fù)原方法有逆濾波，維納濾波，卡爾曼濾波和廣義逆的奇異值分解法等，這些方法已經(jīng)被廣泛地應(yīng)用于圖像復(fù)原上，但是這些方法要求模糊圖像具有較高的信噪比， 如逆濾波的方法僅適用于高信噪比的圖像，這一點(diǎn)限制了傳統(tǒng)的復(fù)原方法在實(shí)際中的應(yīng)用。這些方法的另一個(gè)缺點(diǎn)就是在復(fù)原時(shí)，圖像邊緣不能很好地恢復(fù)，同時(shí)又丟失了一些細(xì) T1 fn 息。上述經(jīng)典的復(fù)原方法不但效果差，而且在實(shí)際應(yīng)用中不能很好的實(shí)現(xiàn)。因此，目前國(guó)際上提出了一些改進(jìn)上述缺點(diǎn)的圖像復(fù)原方法。如，I. Daubechies等人提出基于小波的閾值迭代法，參見(jiàn)文獻(xiàn)〈〈An iterative thresholding algorithm for linear inverse problems with a sparsity constraint》，Commun. Pure Appl. Math. ,2004,Vol. 57,No. 11, PP. 1413-1457。這種方法將兩次迭代所得的復(fù)原結(jié)果的差值作為下一次迭代結(jié)果的補(bǔ)償，是一種有效的復(fù)原方法。但是，這種方法是在小波域進(jìn)行噪聲抑制，容易產(chǎn)生振鈴效應(yīng)，且不能銳化圖像邊緣。此后，J. Bioucas-Dias等人將閾值迭代法進(jìn)行了改進(jìn)，參見(jiàn)文獻(xiàn)《Anew TwIST :two_step iterative shrinkage/thresholding algorithms for image restoration》，IEEETrans. Image Process.，2007，Vol. 16，No. 12，pp. 2992-3004。該方法的收斂速度比一般的閾值迭代法有所提高，同時(shí)，J. Bioucas-Dias等人在他們的代碼示例中，將噪聲系數(shù)轉(zhuǎn)換到全變分域中進(jìn)行抑制，去除了振鈴效應(yīng)，但是這種方法在圖像的平滑區(qū)域容易產(chǎn)生階梯效應(yīng)，且不能很好的恢復(fù)圖像高頻細(xì)節(jié)。

發(fā)明內(nèi)容
本發(fā)明的目的在于克服上述已有技術(shù)的不足，提出一種基于三字典塊匹配的圖像復(fù)原方法，以在圖像復(fù)原時(shí)，能夠去除圖像梯度效應(yīng)，銳化圖像邊緣，恢復(fù)圖像高頻細(xì)節(jié)信息，提高模糊圖像的恢復(fù)質(zhì)量。實(shí)現(xiàn)本發(fā)明目的的技術(shù)方案為通過(guò)將待處理的模糊圖像分塊，分別找到各圖像塊在模糊字典中最匹配的圖像塊，然后利用三個(gè)字典間的對(duì)應(yīng)關(guān)系恢復(fù)出低頻結(jié)果圖和高頻結(jié)果圖，再將恢復(fù)出的低頻結(jié)果圖和高頻結(jié)果圖相加，得到最終的復(fù)原結(jié)果圖。其具體步驟包括(1)輸入一幅待處理模糊圖像\，初始化模糊圖像)(b的低頻結(jié)果圖& = 0，低頻結(jié)果圖&的大小與模糊圖像\大小相同；初始化模糊圖像\的高頻結(jié)果圖\ = 0，高頻結(jié)果圖\的大小與模糊圖像\相同；(2)將待處理模糊圖像)(b進(jìn)行大小為5X5的分塊，分塊過(guò)程中相鄰塊之間重疊4 個(gè)像素，得到圖像塊集合P = {p(i) |i = 1,2,…，G}，初始化i = 1 ；(3)設(shè)構(gòu)造字典所用的清晰樣本圖像的個(gè)數(shù)M = 5，模糊字典Db、清晰字典D。和高頻字典Dh分別含有的圖像塊數(shù)目N = 300000，聚類數(shù)目η = 600 ；(4)對(duì)M幅清晰樣本圖像進(jìn)行預(yù)處理并采樣，分別得到模糊字典Db，清晰字典Dc和高頻字典Dh (4a)對(duì)M幅清晰樣本圖像用與待處理模糊圖像)(b相同的模糊核進(jìn)行模糊化，得到 M幅模糊樣本圖像，對(duì)這M幅模糊樣本圖像分別按照從左到右，從上到下的順序提取5X5的圖像塊，把得到的N個(gè)圖像塊排成一行，得到模糊字典Db ；(4b)對(duì)M幅清晰樣本圖像分別按照從左到右，從上到下的順序提取5X5的圖像塊，把得到的N個(gè)圖像塊排成一行，得到清晰字典D。；(4c)對(duì)M幅清晰樣本圖像進(jìn)行高斯高通濾波，得到M幅高頻圖像，對(duì)這M幅高頻圖像分別按照從左到右，從上到下的順序提取5X5的圖像塊，把得到的N個(gè)圖像塊排成一行， 得到高頻字典Dh;(5)將模糊字典Db中的所有圖像塊用k-means方法聚成η類，得到聚類中心集C ={C(k) |k= 1,2, ...，n}；(6)提取模糊圖像)(b中的第i個(gè)圖像塊P(i)，在聚類中心集C= {C(k)|k= 1， 2，…，η}中搜索與圖像塊P(i)最匹配的聚類中心C(j)，C(j)為聚類中心集C中的第j個(gè)圖像塊；(7)搜索聚類中心C(j)子類中與圖像塊P(i)最匹配的前5個(gè)圖像塊P' b(t)，t =1,2,…，5，通過(guò)模糊字典Db、清晰字典D。和高頻字典Dh這三個(gè)字典圖像塊之間的位置對(duì)應(yīng)關(guān)系，分別找到清晰字典D。中與P' b(t)相對(duì)應(yīng)的5個(gè)清晰圖像塊P'。(t)和高頻字典Dh中與P' b(t)相對(duì)應(yīng)的5個(gè)高頻圖像塊P' h(t)；(8)通過(guò)非局部權(quán)重公式求出與圖像塊P(i)最匹配的5個(gè)圖像塊P' b(t)對(duì)于圖像塊p(i)的權(quán)重系數(shù)W(t)，對(duì)清晰字典D。中的5個(gè)清晰圖像塊P'。(t)進(jìn)行非局部協(xié)作處理，得到低頻結(jié)果圖像塊&(1)；對(duì)高頻字典Dh中的5個(gè)高頻圖像塊P' h(t)進(jìn)行非局部協(xié)作處理，得到高頻結(jié)果圖像塊I\(i)，將PJi)和I\(i)放到低頻結(jié)果圖&和高頻結(jié)果圖\中對(duì)應(yīng)位置，同時(shí)將圖像塊P (i)的標(biāo)記修改為i = i+Ι ；(9)判斷圖像塊標(biāo)記i是否大于待處理模糊圖像\中所含的圖像塊總數(shù)G，如果i >G，則執(zhí)行步驟(10)；否則，返回步驟(6)，直到滿足條件i > G為止；(10)分別對(duì)低頻結(jié)果圖&和高頻結(jié)果圖\中的像素重疊部分求平均，得到求平均后的低頻結(jié)果圖足與求平均后的高頻結(jié)果圖足，最終的復(fù)原結(jié)果圖X為足與足相加之和， 即;T = Iz+足。本發(fā)明與現(xiàn)有的技術(shù)相比具有以下優(yōu)點(diǎn)1、本發(fā)明利用模糊字典Db、清晰字典D。和高頻字典Dh之間的塊匹配技術(shù)，與基于小波或全變分的圖像復(fù)原方法相比，解決了復(fù)原結(jié)果圖容易產(chǎn)生階梯效應(yīng)的問(wèn)題；2、本發(fā)明采用非局部協(xié)作處理方法重構(gòu)復(fù)原圖像，能夠銳化圖像邊緣，且能夠恢復(fù)更多的圖像細(xì)節(jié)。


圖1是本發(fā)明的流程圖2是本發(fā)明建立的三個(gè)字典的子流程圖3是本發(fā)明中用三個(gè)字典塊匹配方法進(jìn)行圖像復(fù)原的示意圖4是本發(fā)明在仿真實(shí)驗(yàn)中建立三個(gè)字典時(shí)用到的5幅清晰樣本圖像；
圖5是本發(fā)明仿真實(shí)驗(yàn)中使用的Cameraman清晰圖像；
圖6是本發(fā)明仿真實(shí)驗(yàn)中使用的Cameraman模糊圖像；
圖7是本發(fā)明在仿真實(shí)驗(yàn)中得到的Cameraman復(fù)原圖8是用現(xiàn)有的基于小波的閾值迭代法在仿真實(shí)驗(yàn)中得到的Cameraman復(fù)原 圖。圖9是用現(xiàn)有的基于全變分的閾值迭代法在仿真實(shí)驗(yàn)中得到的Cameraman復(fù)原
具體實(shí)施例方式參照?qǐng)D1，發(fā)明的具體實(shí)現(xiàn)步驟如下步驟1，輸入一幅待處理模糊圖像\，初始化模糊圖像)(b的低頻結(jié)果圖&為空矩陣，即& = 0，初始化模糊圖像)(b的高頻結(jié)果圖\為空矩陣，即\ = 0，該低頻結(jié)果圖&的矩陣大小和高頻結(jié)果圖\的矩陣大小均與待處理模糊圖像\的大小相同。步驟2，對(duì)待處理模糊圖像)(b進(jìn)行大小為5X5的分塊，分塊過(guò)程中相鄰塊之間重疊4個(gè)像素，得到圖像塊集合P= {P(i)|i = 1,2,…，G}，初始化i = 1。步驟3，設(shè)構(gòu)造字典所用的清晰樣本圖像的個(gè)數(shù)M = 5，模糊字典Db、清晰字典Dc和高頻字典Dh分別含有的圖像塊總數(shù)目N = 300000，聚類數(shù)目η = 600。步驟4，對(duì)M幅清晰樣本圖像進(jìn)行預(yù)處理并采樣，分別得到模糊字典Db，清晰字典 Dc和高頻字典Dh。參照?qǐng)D2，本步驟的具體實(shí)現(xiàn)如下(4a)對(duì)M幅清晰樣本圖像用與待處理模糊圖像)(b相同的模糊核進(jìn)行模糊化，得到 M幅模糊樣本圖像，對(duì)這M幅模糊樣本圖像分別按照從左到右，從上到下的順序提取5X5的圖像塊，把得到的N個(gè)圖像塊排成一行，得到模糊字典Db ；(4b)對(duì)M幅清晰樣本圖像分別按照從左到右，從上到下的順序提取5X5的圖像塊，把得到的N個(gè)圖像塊排成一行，得到清晰字典D。；(4c)對(duì)M幅清晰樣本圖像進(jìn)行高斯高通濾波，得到M幅高頻圖像，對(duì)這M幅高頻圖像分別按照從左到右，從上到下的順序提取5X5的圖像塊，把得到的N個(gè)圖像塊排成一行， 得到高頻字典Dh。步驟5，將模糊字典Db中的所有圖像塊，利用由Tapas Kanungo等人在文獻(xiàn)《An efficient k-means clustering algorithm analysis and implementation)), IEEE transactions on pattern analysis and machine intelligence,2002, Vol.24, No. 7, pp. 881-892中提出的k-means方法聚成η類，得到聚類中心集C = {C(k) k= 1，2，…，η}。步驟6，提取待處理模糊圖像)(b中的第i個(gè)圖像塊P(i)，在聚類中心集C = {C(k) |k= 1,2,…，η}中搜索與圖像塊P(i)最匹配的聚類中心，其搜索過(guò)程通過(guò)如下匹配度公式實(shí)現(xiàn)，即找出NORMl (k)，k = 1,2,…，η中的最小值NORMl (j)，則NORMl (j)對(duì)應(yīng)的聚類中心C(j)就是與圖像塊P(i)最匹配的聚類中心，匹配度公式為
權(quán)利要求
1.一種基于三字典塊匹配的圖像復(fù)原方法，包括如下步驟(1)輸入一幅待處理模糊圖像\，初始化模糊圖像\的低頻結(jié)果圖&= 0，低頻結(jié)果圖 X1的大小與模糊圖像\大小相同；初始化模糊圖像\的高頻結(jié)果圖\ = 0，高頻結(jié)果圖\ 的大小與模糊圖像\相同；(2)將待處理模糊圖像)(b進(jìn)行大小為5X5的分塊，分塊過(guò)程中相鄰塊之間重疊4個(gè)像素，得到圖像塊集合P = {P⑴Ii = 1,2,…，G}，初始化i = 1 ；(3)設(shè)構(gòu)造字典所用的清晰樣本圖像的個(gè)數(shù)M= 5，模糊字典Db、清晰字典D。和高頻字典Dh分別含有的圖像塊數(shù)目N = 300000，聚類數(shù)目η = 600 ；(4)對(duì)M幅清晰樣本圖像進(jìn)行預(yù)處理并采樣，分別得到模糊字典Db，清晰字典D。和高頻字典Dh:(4a)對(duì)M幅清晰樣本圖像用與待處理模糊圖像)(b相同的模糊核進(jìn)行模糊化，得到M幅模糊樣本圖像，對(duì)這M幅模糊樣本圖像分別按照從左到右，從上到下的順序提取5 X 5的圖像塊，把得到的N個(gè)圖像塊排成一行，得到模糊字典Db ；(4b)對(duì)M幅清晰樣本圖像分別按照從左到右，從上到下的順序提取5X5的圖像塊，把得到的N個(gè)圖像塊排成一行，得到清晰字典D。；(4c)對(duì)M幅清晰樣本圖像進(jìn)行高斯高通濾波，得到M幅高頻圖像，對(duì)這M幅高頻圖像分別按照從左到右，從上到下的順序提取5 X 5的圖像塊，把得到的N個(gè)圖像塊排成一行，得到高頻字典Dh ；(5)將模糊字典Db中的所有圖像塊用k-means方法聚成η類，得到聚類中心集C= {C(k) |k= 1,2, ...，n}；(6)提取模糊圖像)(b中的第i個(gè)圖像塊P(i)，在聚類中心集C={C(k)|k= 1,2,…， η}中搜索與圖像塊P(i)最匹配的聚類中心C(j)，C(j)為聚類中心集C中的第j個(gè)聚類中心；(7)搜索聚類中心C(j)子類中與圖像塊P(i)最匹配的前5個(gè)圖像塊P'b(t)，t = 1， 2，…，5，通過(guò)模糊字典Db、清晰字典D。和高頻字典Dh這三個(gè)字典圖像塊之間的位置對(duì)應(yīng)關(guān)系，分別找到清晰字典D。中與P' b(t)相對(duì)應(yīng)的5個(gè)清晰圖像塊P'。(t)和高頻字典Dh中與P' b(t)相對(duì)應(yīng)的5個(gè)高頻圖像塊P' h(t)；(8)通過(guò)非局部權(quán)重公式求出與圖像塊P(i)最匹配的5個(gè)圖像塊P'b(t)對(duì)于圖像塊 P(i)的權(quán)重系數(shù)w(t)，對(duì)清晰字典D。中的5個(gè)清晰圖像塊P'。(t)進(jìn)行非局部協(xié)作處理， 得到低頻結(jié)果圖像塊P1G)；對(duì)高頻字典Dh中的5個(gè)高頻圖像塊P' h(t)進(jìn)行非局部協(xié)作處理，得到高頻結(jié)果圖像塊I\(i)，將Mi)和Ph(i)放到低頻結(jié)果圖X1和高頻結(jié)果圖\中對(duì)應(yīng)位置，同時(shí)將圖像塊P(i)的標(biāo)記修改為i = i+Ι ；(9)判斷圖像塊標(biāo)記i是否大于待處理模糊圖像)(b中所含的圖像塊總數(shù)G，如果i> G，則執(zhí)行步驟(10)；否則，返回步驟(6)，直到滿足條件i > G為止；(10)分別對(duì)低頻結(jié)果圖&和高頻結(jié)果圖\中的像素重疊部分求平均，得到求平均后的低頻結(jié)果圖足與求平均后的高頻結(jié)果圖足，最終的復(fù)原圖X為足與足相加之和，即Jf = X； + Jf Α O
2.根據(jù)權(quán)利要求1所述的基于三字典塊匹配的圖像復(fù)原方法，其中步驟(6)所述的在聚類中心集C= {C(k)|k=l，2，…，η}中搜索與圖像塊P(i)最匹配的聚類中心，其搜索過(guò)程通過(guò)如下匹配度公式實(shí)現(xiàn)，即找出NORMl (k)，k = 1，2，…，η中的最小值NORMl (j)，則 NORMl (j)對(duì)應(yīng)的聚類中心C(j)就是與圖像塊P(i)最匹配的聚類中心，匹配度公式為 NORMl (k) = I P⑴-C(k) I 2其中，NORMl (k)為匹配度估計(jì)函數(shù)，C(k)為聚類中心集C= {C(k) |k= 1,2,…，η}中的第k個(gè)聚類中心；
3.根據(jù)權(quán)利要求1所述的基于三字典塊匹配的圖像復(fù)原方法，其中步驟(7)所述的在聚類中心C(j)所對(duì)應(yīng)的子類中搜索與P(i)相匹配的前5個(gè)圖像塊P' b(t)，t = 1，2，···， 5，其搜索過(guò)程通過(guò)如下匹配度公式實(shí)現(xiàn)，即找出N0RM2 (s)，s= 1,2,…，N中最小的前5個(gè)值，根據(jù)標(biāo)記s找出相對(duì)應(yīng)的Db(S)中的5個(gè)圖像塊，得到與P(i)最相似的5個(gè)圖像塊，匹配度公式為N0RM2 (s) = I IP (i)-Db (s) | 2其中，N0RM2(s)為匹配度估計(jì)函數(shù)，圖像塊Db(S)是第j個(gè)聚類中心C(j)對(duì)應(yīng)子類中的圖像塊，也是模糊字典Db中的第s個(gè)圖像塊，S= 1，2，…，N。找出N0RM2(s)中的最小的5個(gè)值，將對(duì)應(yīng)的Db(S)中的5個(gè)圖像塊賦值給圖像塊P' b(t)，t=l，2，…，5;
4.根據(jù)權(quán)利要求1所述的基于三字典塊匹配的圖像復(fù)原方法，其中步驟(8)所述的通過(guò)非局部權(quán)重公式求出與圖像塊P(i)最匹配的5個(gè)圖像塊P' b(t)對(duì)于圖像塊P(i)的權(quán)重系數(shù)w (t)，其計(jì)算公式為 心c(j)的子類中與圖像塊P(i)最匹配的前5個(gè)圖像塊；
5.根據(jù)權(quán)利要求1所述的基于三字典塊匹配的圖像復(fù)原方法，其中步驟(8)所述的對(duì)清晰字典D。中的5個(gè)清晰圖像塊P'。(t)進(jìn)行非局部協(xié)作處理，得到低頻結(jié)果圖像塊 P1 (i)，其計(jì)算公式為 其中，W(t)為權(quán)重系數(shù)，P'。(t)為清晰字典D。中與P' b(t)相對(duì)應(yīng)的5個(gè)清晰圖像塊；
6.根據(jù)權(quán)利要求1所述的基于三字典塊匹配的圖像復(fù)原方法，其中步驟(8)所述的對(duì)高頻字典Dh中的5個(gè)清晰圖像塊P' h(t)進(jìn)行非局部協(xié)作處理，得到高頻結(jié)果圖像塊 Ph(i)，其計(jì)算公式為 其中，w(t)為權(quán)重系數(shù)，P' h(t)為高頻字典Dh中與P' b(t)相對(duì)應(yīng)的5個(gè)高頻圖像塊。
全文摘要
本發(fā)明公開(kāi)了一種基于三字典塊匹配的圖像復(fù)原方法，主要解決現(xiàn)有技術(shù)在圖像復(fù)原時(shí)，無(wú)法銳化圖像邊緣且造成部分圖像高頻細(xì)節(jié)丟失的問(wèn)題。本發(fā)明的技術(shù)方案為首先輸入一幅待處理模糊圖像，并構(gòu)造模糊字典Db，清晰字典Dc和高頻字典Dh；然后通過(guò)將待處理的模糊圖像分塊，分別找到各圖像塊在模糊字典Db中最匹配的圖像塊；再利用三個(gè)字典間的一一對(duì)應(yīng)關(guān)系恢復(fù)出低頻結(jié)果圖和高頻結(jié)果圖；最后將恢復(fù)出的低頻結(jié)果圖和高頻結(jié)果圖相加，得到最終的復(fù)原結(jié)果圖。本發(fā)明在圖像復(fù)原時(shí)，能夠去除圖像梯度效應(yīng)，銳化圖像邊緣，恢復(fù)圖像高頻細(xì)節(jié)信息，提高模糊圖像的恢復(fù)質(zhì)量?？捎糜趯?duì)各種已知模糊類型的模糊圖像進(jìn)行復(fù)原。
文檔編號(hào)G06T5/00GK102254305SQ20111018024
公開(kāi)日2011年11月23日 申請(qǐng)日期2011年6月30日 優(yōu)先權(quán)日2011年6月30日
發(fā)明者侯彪, 劉忠偉, 季佩媛, 張小華, 焦李成, 王爽, 鐘樺 申請(qǐng)人:西安電子科技大學(xué)