基于非下采樣Shearlet變換與矢量C-V模型的圖像分割方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及圖像處理領(lǐng)域，尤其是一種基于非下采樣Shearlet變換與矢量C-V模 型的圖像分割方法。
【背景技術(shù)】
[0002] 圖像分割的重要性與實用性使它受到越來越多的關(guān)注，與此同時，大量的分割算 法被提出，但是到目前為止還沒有一種簡單的方法可以廣泛地適用于所有的圖像，換言之， 針對某一類圖像發(fā)展起來的分割算法并不一定適用于其他類的圖像。Mumfor-Shah(M-S)分 割技術(shù)是一種基于區(qū)域的圖像分割方法。Chan和Vese通過加入與輪廓區(qū)域特性有關(guān)的項簡 化了 Μ-S方法中的能量泛函，提出了著名的Chan-VeSe(C-V)方法，并且將C-V模型擴(kuò)展到了 矢量情況，矢量情況下的C-V模型繼承了C-V模型的優(yōu)點(diǎn)，可以檢測在標(biāo)量表示下檢測不到 或者是多譜段的目標(biāo)，矢量C-V模型與C-V模型一樣還不能準(zhǔn)確獲得非同質(zhì)圖像的準(zhǔn)確分割 結(jié)果?；谛〔ㄓ虻膱D像處理方法可以充分考慮到圖像的細(xì)節(jié)信息，所以小波變換在圖像 處理的成功運(yùn)用引起了學(xué)者們的極大興趣，提出了很多基于小波域的圖像分割方法。有學(xué) 者嘗試結(jié)合多分辨率表示和C-V模型實現(xiàn)圖像分割，獲得了不錯的分割效果，這些方法既可 以保留圖像的細(xì)節(jié)信息，又可利用C-V模型的整體分割性能。利用小波變換和C-V模型提出 的一種對物體裂痕進(jìn)行快速分割的算法，此方法通過小波極大模值粗略地定位區(qū)域，不但 降低了 C-V模型產(chǎn)生的大量數(shù)據(jù)，且能提供初始輪廓曲線，加快了 C-V模型的收斂速度。但 是，迄今為止還沒有關(guān)于結(jié)合非下采樣的剪切波變換和矢量C-V模型的分割方法的相關(guān)報 道。

【發(fā)明內(nèi)容】

[0003] 本發(fā)明是為了解決現(xiàn)有技術(shù)所存在的上述技術(shù)問題，提供一種基于非下采樣 Shearlet變換與矢量C-V模型的圖像分割方法。
[0004] 本發(fā)明的技術(shù)解決方案是:一種基于非下采樣Shearlet變換與矢量C-ν模型的圖 像分割方法，其特征在于按照如下步驟進(jìn)行： 定義C-ν能量泛函形式如下，
(1) 其中Η是Heaviside函數(shù)，定義如下：
它的分布導(dǎo)數(shù)為
對于固定的C-和f極小化關(guān)于#的函數(shù)料
通過使用最速下降流方法可以得到控制平均曲率和誤差項的Eu 1 er-Lagrange方程參； 矢量C-V模型有如下定義：
其中^^是〇中圖像的等I層信道，
是演化曲線，
和
是近似圖像強(qiáng)度的未知的常數(shù)向量，^ > 〇和是每一層信道的參 數(shù)； 其水平集函數(shù)如下：
衡 為了計算_#:的Euler-Lagrange方程，將士和氣.規(guī)范為： 和
極小化關(guān)于#的能量泛函得到的Euler-Lagrange方程為：
_ 步驟1.對給定的圖像I用NSST進(jìn)行η層分解，首先是用非下采樣的Laplacian金字塔 算法實現(xiàn)圖像的多尺度分解，再利用剪切矩陣獲得方向，具體地： 步驟11.利用Laplacian金字塔分解算法，把分解成一個低通濾波后的圖像f，和 一個高通濾波后的圖像@; 步驟12.在偽極格上計算:|^，生成矩陣麵g:; 步驟13.對矩陣戶^進(jìn)行帶通濾波； 步驟14.直接重新定義笛卡爾抽樣值，向二維快速Fourier變換或應(yīng)用逆采用上一步 中的偽極離散Fourier變換，得到剪切系數(shù)； 步驟2.將用NSST獲得的不同方向的低頻和高頻分量，送到矢量C-V模型中； 步驟3.在整個圖像上初始化一系列均勻分布演化曲線，設(shè)每層多分辨表示的圖像的 正則化參數(shù) 步驟4.極小化由步驟2得到的能量泛函； 步驟5.輸出最后的分割結(jié)果。
[0005] 本發(fā)明NSST可以提供一種靈活的針對圖像的多尺度和多方向性的分解，且能夠保 留圖像的細(xì)節(jié)信息，矢量C-V模型對前一步的圖像多分辨率表示的分割可以獲得圖像分割 的整體效果。與現(xiàn)有技術(shù)相比，本發(fā)明具有以下優(yōu)點(diǎn)：第一，具有較強(qiáng)的方向性。本發(fā)明利用 了 Shearlet變換，具有較強(qiáng)的方向性，能很好地捕獲圖像線面奇異的幾何結(jié)構(gòu)，充分發(fā)揮Ον 模型的分割性能; 第二，具有平移不變性。本發(fā)明采用非下采樣的 Shearlet 隱馬爾可夫樹 模型，通過對圖像非下采樣的Shearlet系數(shù)的概率統(tǒng)計特性的分析，結(jié)合隱馬爾可夫樹模 型，解決了圖像混疊的問題。第三，細(xì)節(jié)信息的保留。本發(fā)明結(jié)合了多尺度變換與矢量c-ν模 型，既可以檢測在標(biāo)量情況下檢測不到的多譜段目標(biāo)，又能保留圖像的細(xì)節(jié)信息。
【附圖說明】
[0006] 圖1是利用本發(fā)明的算法采用不同的NSST分解層數(shù)對同一幅遙感圖像的分割結(jié) 果。
[0007] 圖2是利用不同方法對不同類型圖像的分割結(jié)果比較，其中本發(fā)明算法采用2層 NSST分解，每種分辨率的高頻分量均分解為4個方向。
【具體實施方式】
[0008] 按照如下步驟進(jìn)行： Shearlet變換在多分辨率表示過程中采用了下采樣的操作算子，因而會丟失空間支 撐，仿造非下采樣輪廓波變換(Nonsubsampled Contour let Transform，NSCT)的構(gòu)造方法， 用非下采樣的Laplacian金字塔算法替換Laplacian金字塔算法，構(gòu)造了非下采樣剪切波變 換(Nonsubsampled Shearlet Transform,NSST)〇
[0009] 定義c-v能量泛函形式如下，
(1) 其中Η是Heaviside函數(shù)，定義如下：
它的分布導(dǎo)數(shù)爻
對于固定的e-和e~極小化關(guān)于#的函數(shù)
通過使用最速下降流方法可以得到控制平均曲率和誤差項的Eu 1 er-Lagrange方程參； 矢量C-V模型有如下定義：
其中是Ω中圖像的等?層信道，
是演化曲線，
和
是近似圖像強(qiáng)度的未知的常數(shù)向量，
是每一層信道的參 數(shù)； 其水平集函數(shù)如下：
II 為了計算的Euler-Lagrange方程，將5和-規(guī)范為： 和
極小化關(guān)于#的能量泛函得到的Euler-Lagrange方程為：
步驟1.對給定的圖像I用NSST進(jìn)行η層分解(不同層數(shù)的分解可產(chǎn)生不同的分割結(jié) 果;較小的η值可以保留原始圖像的細(xì)節(jié)，但易受噪聲干擾;較大的η值可以減少噪聲干擾， 但會丟失圖像的細(xì)節(jié)，結(jié)果如圖1所示;建議進(jìn)行2層分解），首先是用非下采樣的Laplacian 金字塔算法實現(xiàn)圖像的多尺度分解，再利用剪切矩陣獲得方向，具體地： 步驟11.利用Laplacian金字塔分解算法，把If4分解成一個低通濾波后的圖像$，和 一個高通濾波后的圖像f ; 步驟12.在偽極格上計算|:，生成矩陣_; 步驟13.對矩陣1?5進(jìn)行帶通濾波； 步驟14.直接重新定義笛卡爾抽樣值，向二維快速Fourier變換或應(yīng)用逆采用上一步 中的偽極離散Fourier變換，得到剪切系數(shù)； 步驟2.將用NSST獲得的不同方向的低頻和高頻分量，送到矢量C-V模型中； 步驟3.在整個圖像上初始化一系列均勻分布演化曲線，設(shè)每層多分辨表示的圖像的 正則化參數(shù)：V"二 步驟4.極小化由步驟2得到的能量泛函； 步驟5.輸出最后的分割結(jié)果。
[0010]將本發(fā)明實施例與典型的矢量c-ν方法和基于非下采樣小波變換的c-ν方法的分 割效果進(jìn)行對比，結(jié)果如圖2所示，可以看出，本發(fā)明方法的分割結(jié)果比其它兩種典型算法 有更好的整體性與視覺效果。
【主權(quán)項】
1. 一種基于非下采樣Shear let變換與矢量C-ν模型的圖像分割方法，其特征在于按照 如下步驟進(jìn)行： 定義C-V能量泛函形式如下，?：Β 其中Η是Heaviside函數(shù)，定義如下：它的分布導(dǎo)數(shù)為_轉(zhuǎn)=; 對于固定的?Γ和【極小化關(guān)于-的函數(shù)洞通過使用最速下降流方法可以得到控制平均曲率和誤差項的Eu 1 er-Lagrange方程; 矢量C-V模型有如下定義：CS) 其中^是e中圖像的等I層信道，是演化曲線，? 和 ?7二:iC.、*、〇是近似圖像強(qiáng)度的未知的常數(shù)向量，零，是每一層信道的 X. ·.；! ：·.. ' -·ν；；：·'· 參數(shù)； 其水平集函數(shù)如下：為了計算學(xué)的Euler-Lagrange方程，將見.和€.規(guī)范為：和極小化關(guān)于#的能量泛函得到的Eu 1 er-Lagrange方程為：步驟1.對給定的圖像J用NSST進(jìn)行η層分解，首先是用非下采樣的Laplacian金字塔 算法實現(xiàn)圖像的多尺度分解，再利用剪切矩陣獲得方向，具體地： 步驟11.利用Laplacian金字塔分解算法，把分解成一個低通濾波后的圖像:^，和 一個高通濾波后的圖像 步驟12.在偽極格上計算:^，生成矩陣1^'; 步驟13.對矩陣進(jìn)行帶通濾波； 步驟14.直接重新定義笛卡爾抽樣值，向二維快速Fourier變換或應(yīng)用逆采用上一步 中的偽極離散Fourier變換，得到剪切系數(shù)； 步驟2.將用NSST獲得的不同方向的低頻和高頻分量，送到矢量C-V模型中； 步驟3.在整個圖像上初始化一系列均勻分布演化曲線，設(shè)每層多分辨表示的圖像的 正則化參數(shù)= 1; 步驟4.極小化由步驟2得到的能量泛函； 步驟5.輸出最后的分割結(jié)果。
【專利摘要】本發(fā)明提出了一種基于非下采樣Shearlet變換與矢量C-V模型的圖像分割方法，屬于圖像處理領(lǐng)域。該方法對給定的圖像用非下采樣剪切波進(jìn)行多層分解，獲得不同方向的低頻和高頻分量，送到矢量C-V模型中，在整個圖像上初始化一系列均勻分布演化曲線，極小化得到的能量泛函，從而得到分割結(jié)果。本發(fā)明分割效果顯著，具有較高的客觀評價質(zhì)量和很好的視覺效果。
【IPC分類】G06T7/00
【公開號】CN105608705
【申請?zhí)枴緾N201610022719
【發(fā)明人】王相海, 方玲玲, 蘇欣, 朱毅歡
【申請人】遼寧師范大學(xué)
【公開日】2016年5月25日
【申請日】2016年1月14日