1. 程式人生 > >影象融合(三)-- 拉普拉斯金字塔

影象融合(三)-- 拉普拉斯金字塔

2、拉普拉斯金字塔融合

  影象金字塔方法的原理是:將參加融合的的每幅影象分解為多尺度的金字塔影象序列,將低解析度的影象在上層,高解析度的影象在下層,上層影象的大小為前一層影象大小的1/4。層數為0,1,2……N。將所有影象的金字塔在相應層上以一定的規則融合,就可得到合成金字塔,再將該合成金字塔按照金字塔生成的逆過程進行重構,得到融合金字塔。這個總的思路就是一下所有基於金字塔融合的演算法過程,不同點就在於分解構造的金字塔不同,每層的融合規則不一樣,重構的方法不同而已。金字塔方法最先實現了這種思想,之後小波方法進一步完善和發展了這種多尺度融和的思想。

2.1、原理闡述

(1)高斯金字塔

  高斯金字塔是最基本的影象塔。首先將原影象作為最底層影象G0(高斯金字塔的第0層),利用高斯核(5*5)對其進行卷積,然後對卷積後的影象進行下采樣(去除偶數行和列)得到上一層影象G1,將此影象作為輸入,重複卷積和下采樣操作得到更上一層影象,反覆迭代多次,形成一個金字塔形的影象資料結構,即高斯金字塔。

高斯金字塔的構建過程為:假設高斯金字塔的第L層影象為Gl:

式中N為高斯金字塔頂層層號,Rl和Cl分別為高斯金字塔第l層的行數和列數W(m,n)是一個二維可分離的5*5視窗函式,表示式為:

  由G0,G1,,,GN,就構成了一個高斯金字塔,其中G0為高斯金字塔的底層(與原影象相同)GN為金字塔的頂層。由此可見高斯金字塔的當前層影象就是對其前一層影象首先進行高斯低通濾波,然後再進行隔行和隔列的降2取樣而生成的。前一層影象大小依次為當前層影象大小的4倍。

Opencv中使用pyrdown函式就可以獲得高斯金字塔。

(2)拉普拉斯金字塔

  在高斯金字塔的運算過程中,影象經過卷積和下采樣操作會丟失部分高頻細節資訊。為描述這些高頻資訊,人們定義了拉普拉斯金字塔(Laplacian Pyramid, LP)。用高斯金字塔的每一層影象減去其上一層影象上取樣並高斯卷積之後的預測影象,得到一系列的差值影象即為 LP 分解影象。

將Gl內插方法得到放大影象*Gl,使*Gl的尺寸與*Gl-1的尺寸相同,即放大運算元Expand

該式子實現兩個步驟:在偶數行和列插入0,然後使用下采樣中的高斯核進行濾波處理,得到和l-1層一樣大小的影象。

 

  N為拉普拉斯金字塔頂層的層號LPl是拉普拉斯金字塔分解的第L層影象。由LP0,LP1、LP2…LPN構成的金字塔即為拉普拉斯金字塔。它的每一層L0影象是高斯金字塔本層G0影象與其高一層影象G1經內插放大後圖像*G1的差,此過程相當於帶通濾波,因此拉普拉斯金字塔又稱為帶通金字塔分解。

  內插方法:opencv中有實現的函式pyrup。可以得到*G1。然後在兩個函式作差,相減就可以得到拉普拉斯金字塔。

  求得每個影象的拉普拉斯金字塔後需要對相應層次的影象進行融合,具體的融合規則有,取大、取小,等等。

(3)重構

     對融合後的拉普拉斯金字塔,從其頂層開始逐層從上至下按下式進行遞推,可以恢復其對應的高斯金字塔,並最終可得到原影象G0。就是從最高層開始使用內插的方法。

2.2、融合應用

  影象拉普拉斯金字塔分解的目的是將源影象分別分解到不同的空間頻帶上,融合過程是在各空間頻率層上分別進行的,這樣就可以針對不同分解層的不同頻帶上的特徵與細節,採用不同的融合運算元以達到突出特定頻帶上特徵與細節的目的。即有可能將來自不同影象的特徵與細節融合在一起。

(1)頂層處理

  設LAl和LBl分別為源影象A,B經過拉普拉斯金字塔分解後得到的第l層影象,融合後的結果為LFl。當l=N時,LAN和LBN分別為源影象A,B經過拉普拉斯金字塔分解後得到的頂層影象。對於頂層影象的融合,首先計算以其各個畫素為中心的區域大小為M*N(M、N取奇數且M >= 3、N >= 3)的區域平均梯度:

其中,Ix與Iy分別為畫素f(x,y)在x與y方向上的一階差分,定義如下:

  因此對於頂層影象中的每一個畫素LAN(i, j)和LBN(i, j)都可以得到與之相對應的區域平均梯度GA(i, j)和GB(i, j)。由於平均梯度反映了影象中的微小細節反差和紋理變化特徵,同時也反映出影象的清晰度。一般來說平均梯度越大,影象層次也豐富,則影象越清晰。因此頂層影象的融合結果為:

(2)各層次處理

  當0

 

則其他層次影象的融合結果為:

在得到金字塔各個層次的融合影象LF1、LF2、LFN後。通過前面的重構,便可得到最終的融合影象。

 第二種融合規則:

    採用最高層係數取平均,其餘各層係數絕對值取大的融合策略進行融合。融合後圖像的係數(灰度值)越接近較清晰影象的灰度值就說明融合效果好。

--------------------- 本文來自 易大飛 的CSDN 部落格 ,全文地址請點選:https://blog.csdn.net/stf1065716904/article/details/73335512?utm_source=copy