基于變分貝葉斯推理的圖像盲反卷積算法研究
發(fā)布時(shí)間:2021-09-18 06:54
圖像是一種重要的信息載體,無論日常生活、工業(yè)生產(chǎn)還是科學(xué)研究都依賴于清晰圖像。但圖像采集過程中多種退化源的存在導(dǎo)致圖像模糊,不能為人所用。圖像復(fù)原技術(shù)逆退化過程,復(fù)原模糊圖像得到清晰圖像,能夠幫助人們恢復(fù)丟失的重要信息。本文研究了圖像復(fù)原的基礎(chǔ)理論,改進(jìn)了Kotera等人的基于變分貝葉斯推理的圖像盲反卷積算法。基于自然圖像統(tǒng)計(jì)和貝葉斯建模,對(duì)觀測模糊圖像、潛在的清晰圖像和模糊核應(yīng)用自動(dòng)相關(guān)確定先驗(yàn)?zāi)P汀W詣?dòng)相關(guān)確定先驗(yàn)?zāi)P捅举|(zhì)上是學(xué)生t-分布,可以表示為有限數(shù)量的具有共同均值的高斯分布的疊加,其精度(方差的倒數(shù))服從伽馬分布。該模型能夠與變分貝葉斯推理很好的結(jié)合;谧兎重惾~斯推理和自動(dòng)相關(guān)確定先驗(yàn)?zāi)P?本文提出了求解反卷積問題的圖像盲復(fù)原算法。圖像盲反卷積算法包括兩個(gè)階段,模糊核估計(jì)和非盲反卷積。在改進(jìn)的圖像盲反卷積算法中,考慮到模糊核的固有結(jié)構(gòu),將自動(dòng)相關(guān)確定先驗(yàn)?zāi)P投x在模糊核的梯度上。在求解反卷積問題時(shí),變分貝葉斯推理用分解形式的分布來近似后驗(yàn)概率,并通過最小化度量近似分布和后驗(yàn)分布的差異的Kullback-Leibler散度來得到對(duì)清晰圖像和模糊核的正確估計(jì)。基于變分貝葉斯推...
【文章來源】:哈爾濱工業(yè)大學(xué)黑龍江省 211工程院校 985工程院校
【文章頁數(shù)】:63 頁
【學(xué)位級(jí)別】:碩士
【部分圖文】:
高斯模型對(duì)圖像梯度
哈爾濱工業(yè)大學(xué)工學(xué)碩士學(xué)位論文-15-來源于同一張?jiān)紙D的圖像集合。圖2-6表示的是一般的圖像金字塔結(jié)構(gòu)。將模糊核估計(jì)算法在金字塔圖像集合上運(yùn)行。在最低分辨率的圖像級(jí)中,初始將模糊核設(shè)為delta函數(shù),估計(jì)出小尺寸的模糊核,再將該模糊核上采樣(up-sampling)后作為下一級(jí)模糊核盲估計(jì)的初始值。重復(fù)將由低分辨率圖像估計(jì)的模糊核上采樣后作為下一級(jí)圖像復(fù)原中模糊核盲估計(jì)的初始值,直至到原始圖一級(jí)。這樣可以避免估計(jì)的模糊核為局部最小值,保證模糊核估計(jì)的正確性,也進(jìn)一步加快收斂速度,提高運(yùn)算效率。第J級(jí),原始圖NNN/2N/2442211第J-1級(jí)第2級(jí)第1級(jí)第0級(jí)(頂點(diǎn))圖2-6Lena的圖像金字塔示意圖本課題采用圖像高斯金字塔,圖像高斯金字塔的生成分為兩個(gè)步驟:(1)進(jìn)行高斯模糊(濾波);(2)降采樣(down-sampling)。為了便于描述,這里用一維高斯濾波函數(shù)說明。一維高斯函數(shù)及其傅里葉變換可分別由公式(2-18)(2-19)表示。221()exp22xfx(2-18)22222222()expexp2exp212F(2-19)式中——高斯分布的方差;x,,——分別表示空間坐標(biāo)、圓頻率和空間頻率。傅里葉變換后,高斯函數(shù)仍為高斯函數(shù),只是幅度和方差發(fā)生了變化[52],公式(2-18)(2-19)也說明了這一點(diǎn)。進(jìn)行高斯模糊前,首先要確定高斯模糊的方差,可通過空間截止頻率得到。令空間截止頻率等于奈奎斯特頻率,那么
卣螅輝俑?莨?劍?-52)-(3-56)計(jì)算近似后驗(yàn)分布的參數(shù)。圖像盲反卷積算法的流程圖如圖3-4所示。為了正確估計(jì)尺寸較大的模糊核,必須使用圖像金字塔結(jié)構(gòu),以避免算法陷入局部最小值。本文第二章圖像預(yù)處理中已詳細(xì)介紹如何建立圖像金字塔。首先確定高斯濾波器的尺寸,本課題計(jì)算當(dāng)1()(0)20FF時(shí)的坐標(biāo)值,即傅里葉系數(shù)為最大值的1/20時(shí)的坐標(biāo)值,濾波器關(guān)于原點(diǎn)對(duì)稱,這個(gè)坐標(biāo)為高斯濾波器尺寸的一半。降采樣采用1.5:1的縮放比例,建立5層金字塔結(jié)構(gòu),采樣點(diǎn)的值則通過雙線性插值獲得,MATLAB函數(shù)interp2可實(shí)現(xiàn)。圖3-5描述了模糊核尺寸為33x33,采用1.5:1的縮放比例,每一級(jí)估計(jì)的模糊核。a)b)c)d)e)圖3-5多尺度模糊核逐級(jí)優(yōu)化示意圖模糊核尺寸:a)77;b)1010;c)1515;d)2222;e)3333本文提出的圖像盲反卷積算法同樣借助快速傅里葉變換來加速計(jì)算,所以也需要考慮邊界值問題。這里與文獻(xiàn)[32]一樣將位于圖像邊界處模糊核尺寸一半的像素的i值設(shè)為0,認(rèn)為在圖像邊界區(qū)域的模型差異最大,反卷積算法完全忽略邊界區(qū)域。超參數(shù)()a除a均設(shè)為0;()b設(shè)置為常值,其中b和b設(shè)為710,b設(shè)為410,vb設(shè)為310;3ab10[32]。本課題沒有對(duì)參數(shù)進(jìn)行優(yōu)化,對(duì)超參數(shù)的估計(jì)可參見Tzikas等的工作[31]758。
【參考文獻(xiàn)】:
期刊論文
[1]工程意義上的高斯脈沖信號(hào)不失真采樣分析[J]. 王云靜,張曉林,趙雷. 遙測遙控. 2019(01)
博士論文
[1]基于變分法的圖像復(fù)原算法研究[D]. 童蓓蕾.中國科學(xué)技術(shù)大學(xué) 2018
碩士論文
[1]自然圖像統(tǒng)計(jì)在圖像處理領(lǐng)域的應(yīng)用[D]. 張東.浙江大學(xué) 2013
本文編號(hào):3399688
【文章來源】:哈爾濱工業(yè)大學(xué)黑龍江省 211工程院校 985工程院校
【文章頁數(shù)】:63 頁
【學(xué)位級(jí)別】:碩士
【部分圖文】:
高斯模型對(duì)圖像梯度
哈爾濱工業(yè)大學(xué)工學(xué)碩士學(xué)位論文-15-來源于同一張?jiān)紙D的圖像集合。圖2-6表示的是一般的圖像金字塔結(jié)構(gòu)。將模糊核估計(jì)算法在金字塔圖像集合上運(yùn)行。在最低分辨率的圖像級(jí)中,初始將模糊核設(shè)為delta函數(shù),估計(jì)出小尺寸的模糊核,再將該模糊核上采樣(up-sampling)后作為下一級(jí)模糊核盲估計(jì)的初始值。重復(fù)將由低分辨率圖像估計(jì)的模糊核上采樣后作為下一級(jí)圖像復(fù)原中模糊核盲估計(jì)的初始值,直至到原始圖一級(jí)。這樣可以避免估計(jì)的模糊核為局部最小值,保證模糊核估計(jì)的正確性,也進(jìn)一步加快收斂速度,提高運(yùn)算效率。第J級(jí),原始圖NNN/2N/2442211第J-1級(jí)第2級(jí)第1級(jí)第0級(jí)(頂點(diǎn))圖2-6Lena的圖像金字塔示意圖本課題采用圖像高斯金字塔,圖像高斯金字塔的生成分為兩個(gè)步驟:(1)進(jìn)行高斯模糊(濾波);(2)降采樣(down-sampling)。為了便于描述,這里用一維高斯濾波函數(shù)說明。一維高斯函數(shù)及其傅里葉變換可分別由公式(2-18)(2-19)表示。221()exp22xfx(2-18)22222222()expexp2exp212F(2-19)式中——高斯分布的方差;x,,——分別表示空間坐標(biāo)、圓頻率和空間頻率。傅里葉變換后,高斯函數(shù)仍為高斯函數(shù),只是幅度和方差發(fā)生了變化[52],公式(2-18)(2-19)也說明了這一點(diǎn)。進(jìn)行高斯模糊前,首先要確定高斯模糊的方差,可通過空間截止頻率得到。令空間截止頻率等于奈奎斯特頻率,那么
卣螅輝俑?莨?劍?-52)-(3-56)計(jì)算近似后驗(yàn)分布的參數(shù)。圖像盲反卷積算法的流程圖如圖3-4所示。為了正確估計(jì)尺寸較大的模糊核,必須使用圖像金字塔結(jié)構(gòu),以避免算法陷入局部最小值。本文第二章圖像預(yù)處理中已詳細(xì)介紹如何建立圖像金字塔。首先確定高斯濾波器的尺寸,本課題計(jì)算當(dāng)1()(0)20FF時(shí)的坐標(biāo)值,即傅里葉系數(shù)為最大值的1/20時(shí)的坐標(biāo)值,濾波器關(guān)于原點(diǎn)對(duì)稱,這個(gè)坐標(biāo)為高斯濾波器尺寸的一半。降采樣采用1.5:1的縮放比例,建立5層金字塔結(jié)構(gòu),采樣點(diǎn)的值則通過雙線性插值獲得,MATLAB函數(shù)interp2可實(shí)現(xiàn)。圖3-5描述了模糊核尺寸為33x33,采用1.5:1的縮放比例,每一級(jí)估計(jì)的模糊核。a)b)c)d)e)圖3-5多尺度模糊核逐級(jí)優(yōu)化示意圖模糊核尺寸:a)77;b)1010;c)1515;d)2222;e)3333本文提出的圖像盲反卷積算法同樣借助快速傅里葉變換來加速計(jì)算,所以也需要考慮邊界值問題。這里與文獻(xiàn)[32]一樣將位于圖像邊界處模糊核尺寸一半的像素的i值設(shè)為0,認(rèn)為在圖像邊界區(qū)域的模型差異最大,反卷積算法完全忽略邊界區(qū)域。超參數(shù)()a除a均設(shè)為0;()b設(shè)置為常值,其中b和b設(shè)為710,b設(shè)為410,vb設(shè)為310;3ab10[32]。本課題沒有對(duì)參數(shù)進(jìn)行優(yōu)化,對(duì)超參數(shù)的估計(jì)可參見Tzikas等的工作[31]758。
【參考文獻(xiàn)】:
期刊論文
[1]工程意義上的高斯脈沖信號(hào)不失真采樣分析[J]. 王云靜,張曉林,趙雷. 遙測遙控. 2019(01)
博士論文
[1]基于變分法的圖像復(fù)原算法研究[D]. 童蓓蕾.中國科學(xué)技術(shù)大學(xué) 2018
碩士論文
[1]自然圖像統(tǒng)計(jì)在圖像處理領(lǐng)域的應(yīng)用[D]. 張東.浙江大學(xué) 2013
本文編號(hào):3399688
本文鏈接:http://www.sikaile.net/kejilunwen/shengwushengchang/3399688.html
最近更新
教材專著