【正文】 息，因此可以分離出大部分光譜信息，有利 于圖像光譜信息的保持。同時 , 圖像從 RGB 轉(zhuǎn)換到 IHS這一彩色空間變換有效地分離了代表空間亮度信息的 I 分量和代表光譜彩色信息的 H 和 S分量。 算法融合流程如下圖： I H S 彩 色 空 間 變 換 提 取 I 分 量多 光 譜 段 a多 光 譜 段 b多 光 譜 段 c高 分 辨 率 圖 像以 I 分 量 為 標 準 進 行 直 方 圖 匹 配以 直 方 圖 匹 配 后 的 全 色 圖 像 替 換 原 圖 像 的 I 分 量I H S 反 變 換融 合 結(jié) 果 圖 32 基于 HIS變換融合流程圖 基于 IHS 變換的融合方法特點是 : 算法簡單易于實現(xiàn) , 圖像的高頻細節(jié)信息保留較好但光譜信息損失較大。 小波變換融合法 “小波 ”就是小區(qū)域、長度有限、均值為 0 的波形。 小波的定義及特點 （ 1） 連續(xù)小波 設 為 )(t? 一平方可積函數(shù)，即 )()( 2 RLt ?? ，若其傅里葉變換 )(?? 滿足條件 ??? ?Rd???)(2 (21) 則稱 )(t? 為一個基本小波或小波母函數(shù)。減小 小 波變換系數(shù)兀余度的做法是 將小波 基函數(shù) Rtt ???? ? ταα τατα ,0),()( 21, ??的α，τ限定 在一些離散點上取值，一種最通常的離散方法就是將尺度按冪級數(shù)進行離 散，即取 mm 0αα ? (m 為整數(shù)， 0a ≠ 1，一般取 0a =2)。2(21)22(21)(, ??? 基于 MATLAB 的圖像融合算法 為了簡化起見，把 t 軸用 sT 歸一化 ， 于是上式就變?yōu)?)2(2)( 2, ntt mmnm ?? ?? ?? 。它是計算機 視覺中常用的圖像處理方法，提供了不同尺度下分析函數(shù)的一種手段。上圖是對 0V 空間的三級分解示意圖。設 G 和 H分別為理想的低通和高通濾波器，用這對濾波器對原始信號進行多分辨率分解，則信號經(jīng) G 和 H濾波后兩支路輸出必定正交 (因為頻帶不重疊 )，且兩支路輸出的帶寬均減半，因此采樣率可以減半而不會引起信息丟失 (帶通信號的采樣率決定于其帶寬，而不受限于其頻率的上限 )。按照 二維 Mallat 算法，在尺度 j1 上有如下的 Mallat 分解公式 ： 1111????????jnmjjnmjjnmjjnmjcAGGcDcAGHcVcAHGcHcAHHcA 相應的重構(gòu)公式如下 : jmnjmnjmnjmnj cDGGcVHGcHGHcAHHcA ????????? ???? 1 其中 jjjj cDcVcHcA , 分別對 應于圖像 jCA 的低頻成分、水平方向上的高頻成分、垂直方向上的高頻成分、對角方向上的高頻成 分 。 原圖像H H1L H1H L1H H2H L2L L2L H2L H1H H1L L1H L 1HL 畢業(yè)設計 圖 37 圖像的二級小波分解示意圖 其中 L 為圖像的低頻部分，集中了其主要能量， )3,2,1(, 321 ?jHHH jjj 分別表小第 j層水平、垂直與對角方向的高頻分量， 它們 都是圖像的細節(jié)部分。它也是一種圖像的多分辨率、多尺度分解。 設 A,B 為兩幅原始圖像， F 為融合圖像。這種方法既簡化了計算，同時又具有很好的融合結(jié)果。其中 K、 α、 β為加權(quán)因子。對于不同的圖像，適當調(diào)整 K、 α、 β,可以消減模糊邊緣并確保在消減時不過度地喪失邊緣信息。同一融合算法，對同一圖像，觀察者感興趣的部分不同，則認為的效果不同 。二是增加融合圖像中的信息量或信息的精度及可靠性，為人的決策提供更豐富、更準確、更可靠的圖像信息。盡管主觀 評價通過大量的統(tǒng)計，可以獲得比較準確的判斷結(jié) 果，這個結(jié)果主要是針對人的視覺感受做出的，但是這種評價方式最大的不足是整個過程非常煩瑣，如果沒有進行大量的統(tǒng)計，所獲得的結(jié)論可能不準確。 基于光譜特征的評價 （ 1） 偏差指數(shù) 偏差指數(shù)是用來比較融合圖像和低分辨多光譜圖像偏離程度的。客觀評價可分為兩類 :一類是用來評判融合圖像所保持的光譜特征情況的評價指標，主要有偏差指數(shù)，光譜扭曲程度和相關系數(shù)等 。觀察者根據(jù)一些事先規(guī)定的評價尺度或自己的經(jīng)驗，對被評價的圖像提出質(zhì)量判決。因此，探究圖像融合效果的客觀、定量評價方法十分必要。 畢業(yè)設計 第四章 圖像融合效果評價 圖像融合效果的評價問題是一項重要 而有意義的工作。因子 K 調(diào)節(jié)兩幅圖像的戰(zhàn)優(yōu)比例，使兩幅亮度不同的圖像達到均衡。在 本融合實驗中 圖像融合過程中，融合規(guī)則及融合算子的選擇 如下 。 對各分解層進行小波反變換，上一層的小波反變換的結(jié)果就是下一層 小波金子塔圖像序列的低頻子圖像，依次類推，最終所得到的圖像就是融合圖像。通過在這種融合圖像的小波表示 逆變換即可得到融合圖像。對一幅灰度圖像進行 N層的小波分解，形成 3N+1個不同頻帶的數(shù)據(jù)，其中有 3N 個包含細節(jié)信息的高頻帶和一個包含近似分量的低頻帶。 將圖像 做一次小波分解，即將圖像分解成低頻近似分量、水平高頻分量、垂直高頻分量和對角高頻分量，二維 圖像 數(shù)據(jù)經(jīng) 三 次小波分解的塔形框架如圖 33所示，圖中 下標 表示小波分解的層數(shù)。然后再對低頻輸出分量再一次分別應用 G 和 H進行濾波，再降 2 采樣，如此重復濾波就實現(xiàn)了對原始信號的多分辨率分解，如圖 32 所示，圖中↓ 2表示降 2采樣 。由式 (31)可得各子空間之間的以下特性 : 121122110 . . .. . . WWWWVWWVWVV NNN ???????????? ? 這個結(jié)果說明，分辨率為 02 =1 的多分辨率分析子空間 0V 可以用有限個子空間來逼近。多分辨率分析可以形象地表示為一組嵌套的多分辨率子空間，小波變換是一種多分辨率分析的強有力工具。當 A=B 時，則稱該框架為緊框架， 也就是說離散小波序列 ? ?Zkjkj ?,?構(gòu)成一個框架 。為了不丟失信息，我們要求采樣間隔 τ 滿足 Nyqulst 采樣定理，即采樣頻率大于等于該尺 度下頻率通帶的二倍。 畢業(yè)設計 將小 波母函數(shù) )(t? 進行伸縮平移，設其伸縮因子 (尺度因子 )為 α ，平移因子為 τ，平移伸縮后的函數(shù)為 )(, tτα? ，則有 Rtt ???? ? ταα τατα ,0),()( 21, ?? 稱 )(, tτα? 為連續(xù)小波基函數(shù)。與 Fourier 變換相比，小波變換是時間 (空間 )頻率的局部化分析，它通過伸縮平移運算對信號 (函數(shù) )逐步進行多尺度細化，最終達到高頻處時間細分，低頻處頻率細分，能自動適應時頻信號分析的要求，從而可聚焦到信號的任意細節(jié)，解決了 Fourier 變換的困難 問題，成為繼 Fourier 變換以來在科學方法上的重大突破。比如土地利用、植被監(jiān)測等項目更注重遙感圖像的色彩信息 , 所以普通的 IHS 方法并不能很好的滿足要求。在此介紹 Harrison和 Jupp于 1990年提出的變換模型，它是目前在多光譜圖像融合領域中最常見的變換模型。通常的色彩顯示是通過 RGB( 紅綠藍 )信號的亮度值所確定的 , RGB 彩色坐標系統(tǒng)中 R、 G、 B呈非線性關系 ,使調(diào)整色調(diào)的定量操作較為困難。而色調(diào)和飽和度又統(tǒng)稱為色度，反映顏色信息。 第一主成分分量低 空 間 分 辨 率 多 光 譜 圖 像高 空間 分辨 率全 色圖 像P C A 反 變 換以 第 一 主 成 分 為 標 準 直 方 圖 匹 配 全 色 圖P C A 正 變 換第二主成分分量第N主成分分量融 合 圖 像 圖 31 基于 PCA 變換的圖像融合方法 PCA 融合算法的優(yōu)點在于，它適用于多光譜圖像的所有波段 (IHS 變換只能用 3 個波段 )，但其不足之處在于，由于 PCA 融合算法中只是用高分辨率圖像簡單替換低分辨率圖像的第一主成分，故低分辨率圖像第一主成分分量會損失一部分反映光譜特性的信息，使畢業(yè)設計 得融合后圖像的光譜畸變嚴重。 ④將變換矩陣 T 代入 Y=TX，將得到 KL 變換的具體表達式 ： XUXuuuuuuuuuY Tmmmmmm??????????????.....................212221212111 式中 Y 矩陣的行向量 ???JY 為第 j 個 主分量。 一般圖像的線性變換可以用下面的式子表示 ： TXY? 式中 X 為待變換圖像數(shù)據(jù)矩陣， Y 為變換后的數(shù)據(jù)矩陣， T 為實現(xiàn)這一線性變換的變換矩陣。能否在各個變量之間相關關系研究的基礎上，用較少的新變量代替原來較多的變量，而且使這些較少的新變量盡可能多地保留原來的變量所反應的信息，主分量分析就是實現(xiàn)這個目標的 一種強有力的方法，它是把原來多個變量化為少數(shù)幾個綜合指標的一種統(tǒng)計分析方法，從數(shù)學角度來看，這是一種降維處理技術(shù)，用較少的幾個綜合指標來代替原來較多的變量指標，而且使這些較少的綜合指標既能盡量地反映原來較多指標所反映的信息，同時它們之間又彼此獨立。但在多數(shù)應用場合，這些簡單的圖像融合方法的局限性是顯而易見的，無法獲得滿意的融合效畢業(yè)設計 果。當融合圖像的灰度值差異很大時，就會出現(xiàn)明顯的拼接痕跡，不利于人眼識別和后續(xù)的目標識別過程 。 ω1 表示加權(quán)系數(shù) 1； ω2表示加權(quán)系數(shù) 2； 通常 ω 1+ω 2=1。 加權(quán)平均方法將源圖像對應像素的灰度值進行加權(quán)平均，生成新的圖像，它是最直接的融合方法。限制方法是確定一個閾值，僅僅選取 R 值大于這個閾值的點作為角點。 該算法易受噪聲的影響。這也就說明該點是一個角點。至少需要 6 對（推薦 12 對）匹配點。 Matlab 中提供有 4 次冪的實現(xiàn)，分別至少需要 6， 10， 10 對匹配點。 如果輸入圖像呈現(xiàn)傾斜，翹起現(xiàn)象，選此法。 2．仿射（ affine）：將平行線轉(zhuǎn)換成平行線。 1．線性正投影（ linear conformal）：最簡單。 利用 Matlab Image Processing Toolbox 中的圖像配準工具實現(xiàn)線性正投影、仿射、投影、多項式、分段線性、局部加權(quán)平均配準的過程。 按自動化程度可以分為人工、半自動和全自動三種類型 。在特征匹配前，首先要從待匹配的多幅圖像中提取出灰度變化明顯的點、線、區(qū)域等特征，組成特征集。該方法不需要對圖像作特征提取，直接利用全部可用的圖像灰度，提高了配準算法的估計精度和魯棒性，但是由于匹配點周圍區(qū)域的點的灰度都參與計算，因此，其 計算量大，速度慢。圖像配準算法就是設法建立兩幅圖像之間的對應關系 ,確定相應幾何變換參數(shù) ,對兩幅圖像中的一幅進行幾何變換的方法。窗口下移一行，再從左到右進行。 ② 不處理原始圖像 4個邊緣的像素，使處理后 4個邊緣的像素灰度保持原始灰度，或人為地賦予特殊灰度。矩陣中心的元素值占的比例越小，平滑作用越強，會使圖像越模糊； ),( yxh 的作用域越大，平滑作用越強，亦導致圖像越模糊。 鄰域平均法 鄰域平均法按 ),(),(),( yxhyxfyxg ?? 直接對圖像作低通濾波處理。常用的空間濾波去噪方法主要有均值濾波、中值濾波、維納濾波以及空間域低通濾波法，其中，維納濾波是一種自適應的濾波方法，比一般的濾波器具有更好的選擇性，可以更好地保護圖像邊緣，尤其對加性高斯白噪聲的抑制效果很好。但相同的是都會導致圖像失真，使圖像模糊，甚至掩蓋圖像的重要特征，給后續(xù)的圖像分析和處理帶來不利影響。針對圖像某一部分或整幅圖像曝光不足使用灰度變換，其目的是增強圖像灰度對比度。 輸 入原 始數(shù) 字圖 像準備工作建 立糾 正變 換函 數(shù)確 定輸 出影 視范 圍逐 個 像素 的 幾何 位 置變 換像 素亮 度值 重采 樣輸 出糾 正后 的圖 像 圖 21 幾何校正處理過程 （ 2）灰度校正 根據(jù)圖像不同失真情況以及所需的不同圖像特征可以采用不同的修正方法。 圖像的校正 圖像校正主要分為兩類：幾何校正和灰度校正。涉及圖像融合 前的預處理，濾波去噪和圖像配準，采取手動配準獲取精準的配準圖像，在空域內(nèi)的像素級圖像融合采用簡單直接的加權(quán)平均法，在頻域內(nèi)的小波變換融合，最后對融合結(jié)果進行評價。 到 80 年代末，人們才開始將圖像融合技術(shù)應用于一般圖像處理 (可見光圖像、紅外圖像等 )。 90 年代開始，圖像融合技術(shù)開始成為遙感圖像處理和分析中的研究熱點之一。由于