[發明專利]基于Tetrolet變換的多通道衛星云圖融合方法

有效

• 申請號: 201310729319.3
• 申請日: 2013-12-25
• 公開（公告）號: CN103700075A
• 公開（公告）日: 2014-04-02
• 發明（設計）人: 張長江;陳源
• 申請（專利權）人: 浙江師范大學
• 主分類號: G06T5/00
• 分類號: G06T5/00;G06T5/50
• 代理公司: 杭州華鼎知識產權代理事務所(普通合伙) 33217
• 代理人: 胡根良

1.基于Tetrolet變換的多通道衛星云圖融合方法，其特征在于：首先對待融合的多通道衛星云圖分別進行圖像直方圖均衡化處理，再分別進行Tetrolet變換，得到低頻系數、高頻系數及其相對應的拼板覆蓋值；然后，在Tetrolet域低頻部分，利用拉普拉斯金字塔進行再次分解，對其頂層取均值，其他層取灰度絕對值大的部分后重構；在Tetrolet域高頻部分，取每個圖像拼板中高頻系數標準差較大的部分，拼板覆蓋值也取其相對應的值，最后，通過Tetrolet逆變換得到最終的融合圖像。

2.如權利要求1所述的基于Tetrolet變換的多通道衛星云圖融合方法，其特征在于：包括以下步驟：

步驟1對配準后的源圖像A和B分別進行直方圖均衡化處理，得到圖像A'和B'；

步驟2對圖像A'和B'進行Tetrolet變換，分解層數為M，得到高頻系數TH_A和TH_B、低頻系數TL_A和TL_B及相應的拼板覆蓋值TC_A和TC_B；

步驟3分別對低頻系數TL_A和TL_B做拉普拉斯金字塔分解，分解層數為N，得到分解圖像LA和LB，第n（0＜n≤N）層子圖分別為LA_n和LB_n；

步驟4對拉普拉斯金字塔頂層子圖LA_N和LB_N用均值法進行融合，得到融合結果LF_N為

LFN(i,j)=LAN(i,j)+LBN(i,j)2]]>

其中，0＜i≤CL_N，0＜j≤RL_N，CL_N是分解子圖第N層圖像的行數，RL_N是分解子圖第N層圖像的列數；

步驟5對拉普拉斯金字塔其他層子圖LA_n和LB_n用灰度絕對值取大的融合規則融合，則融合結果LF_n為

LFn(i,j)=LAn(i,j),|LAn(i,j)|≥|LBn(i,j)|LBn(i,j),|LAn(i,j)|<|LBn(i,j)|;]]>

步驟6將融合后得到的拉普拉斯金字塔LF進行重構，得到低頻部分的融合結果TL_F；

步驟7在Tetrolet變換域高頻系數部分，對每個分解子塊按照標準差取大的融合規則融合，記第m（0＜m≤M）層p行q列的高頻系數分別為(TH_A)_m(p,q)和(TH_B)_m(p,q)，其大小為12×1的矩陣，其中0＜p≤CH_m，0＜q≤RH_m，CH_m是分解子圖第m層高頻系數的行數，RH_m是分解子圖第m層高頻系數的列數，則

(THF)m(p,q)=(THA)m(p,q),std((THA)m(p,q))≥std((THB)m(p,q))(THB)m(p,q),std((THA)m(p,q))<std((THB)m(p,q))]]>

其中std是標準差，表示為

std=Σi=1KΣj=1G(x(i,j)-x‾)2K×G-1]]>

這里，是x的均值，1≤i≤K，1≤j≤G。在Tetrolet域的高頻系數矩陣中K＝12，G＝1；

步驟8Tetrolet變換中的拼板覆蓋值的選取根據高頻部分的選擇而定，即

(TCF)m(p,q)=(TCA)m(p,q),THm(p,q)=(THA)m(p,q)(TCB)m(p,q),THm(p,q)=(THB)m(p,q)]]>

其中，(TC_A)_m(p,q)和(TC_B)_m(p,q)分別表示圖像A和B經Tetrolet變換后第m層p行q列的拼板覆蓋值；

步驟9對融合處理后的Tetrolet系數值進行Tetrolet逆變換，得到最終的融合圖像F。