干涉成像光譜儀切趾函數(shù)對(duì)復(fù)原光譜的影響分析

干涉成像光譜儀切趾函數(shù)對(duì)復(fù)原光譜的影響分析
摘要：隨著搭載干涉成像光譜儀HJY20—1一A的我國(guó)環(huán)境與減災(zāi)遙感衛(wèi)星HJ一1A即將發(fā)射，我國(guó)干涉光譜成像研究 也從實(shí)驗(yàn)室開始走向?qū)嵱没�．在干涉光譜成像過程中，切趾函數(shù)處理是干涉成像光譜儀光譜復(fù)原過程中的一個(gè)重 要環(huán)節(jié)，對(duì)復(fù)原光譜的精度有著極其重要的影響．根據(jù)HJY20—1一A的參數(shù)設(shè)置，文中首先模擬了24種典型地物對(duì) 應(yīng)于HJY20—1一A和其它最大光程差設(shè)置的干涉成像光譜儀數(shù)據(jù)，在不同切趾函數(shù)作用下的復(fù)原光譜，結(jié)果表明 Hanning函數(shù)是其中最有效、最為穩(wěn)定的切趾函數(shù)，同時(shí)發(fā)現(xiàn)切趾函數(shù)的應(yīng)用雖然可以提高復(fù)原光譜的精度，但與 真實(shí)光譜仍存在一定差距，尤其對(duì)應(yīng)HJY20—1一A，復(fù)原光譜的精度更加有限．在以上分析基礎(chǔ)上，提出了基于儀器線 型函數(shù)標(biāo)準(zhǔn)化的光譜復(fù)原改進(jìn)算法，實(shí)驗(yàn)結(jié)果證實(shí)了該方法可以顯著提高復(fù)原光譜精度，尤其適用于最大光程差 較小的空間調(diào)制型干涉成像光譜儀．最后，就HJY20—1．A復(fù)原光譜對(duì)3種典型植被指數(shù)求解，進(jìn)一步證明了該方法 的有效性．

2．1 光譜復(fù)原算法 切趾函數(shù)滿足的必要條件為：在零光程差處為 極大值1，隨著光程差的增大函數(shù)值變小，在最大光 程差 處為零．試驗(yàn)中采用了矩形函數(shù)、三角函數(shù)、 Hanning函數(shù)和Blackman函數(shù)四種常用切趾函數(shù)進(jìn) 行影響分析，這四種切趾函數(shù)以及對(duì)應(yīng)的儀器線型函 數(shù)(ILS)如表1所示．由于干涉成像光譜儀探測(cè)到的 干涉信號(hào)是有限光程差范圍內(nèi)的信號(hào)，即干涉信號(hào)在 最大光程差位置被突然截?cái)�，這相當(dāng)于對(duì)干涉成像光 譜儀探測(cè)到的信號(hào)采用了矩形函數(shù)進(jìn)行處理． 切趾函數(shù)W( 對(duì)應(yīng)的儀器線型函數(shù)／LS可以 看作是通過最大光程差為 的干涉成像光譜儀的波 數(shù)為 ．的單色光，經(jīng)過切趾函數(shù)W( 作用后的復(fù) 原光譜數(shù)據(jù)．若波數(shù)為 的單色光輻亮度值為 B( 。)，則復(fù)原后的輻亮度值B ( )為：日 ( 。)= B( ) ／LS( 表示卷積)． 因此，當(dāng)波數(shù)范圍[ ]的入瞳光輻亮度 是B( )時(shí)，其經(jīng)過切趾函數(shù)W( 作用的復(fù)原光譜 輻亮度值日 ( )是： 盯m ax 日 ( )=ΣB(o-) 兒s ， (2) O’mi 不同切趾函數(shù)作用下的復(fù)原光譜精度存在差異， 而干涉成像光譜儀的最大光程差不同，復(fù)原光譜 的精度也差異很大．因此，測(cè)厚儀| 測(cè)速儀| 轉(zhuǎn)速表| 壓力表| 壓力計(jì)| 真空表| 硬度計(jì)| 探傷儀| 電子稱| 熱像儀| 頻閃儀| 測(cè)高儀| 測(cè)距儀| 金屬探測(cè)器| 試驗(yàn)機(jī)| 扭力計(jì)| 流速儀| 粗糙度儀| 流量計(jì)| 平衡儀| 本文中對(duì)干涉成像光譜 儀分別設(shè)置了四種最大光程差值，以分析不同最 大光程差時(shí)切趾函數(shù)對(duì)24種典型地物復(fù)原光譜 精度的影響．這四種最大光程差值是0．4cm、0． 1cm、0．05cm以及HJY20．1-A對(duì)應(yīng)的最大光程差， 其中HJY20．1．A的最大光程差由光譜分辨率求得 是0．0069cm．

2．2 切趾函數(shù)對(duì)復(fù)原光譜的影響分析 針對(duì)24種典型地物，根據(jù)公式(2)求得對(duì)應(yīng)不 同最大光程差、不同切趾函數(shù)的復(fù)原光譜．這些不同 條件下的24種典型地物復(fù)原光譜的相對(duì)誤差根據(jù) 公式(3)求得，相對(duì)誤差圖如圖2所示． Error( )： ． (3) L 圖2(a)～2(d)分別對(duì)應(yīng)不同最大光程差時(shí)， 四種切趾函數(shù)作用下的24種典型地物復(fù)原光譜的 相對(duì)誤差．從圖中可以看出，雖然地物類型不同，光 譜信息存在差異，但復(fù)原光譜的相對(duì)誤差結(jié)果卻幾 乎相同．這表明切趾函數(shù)作用下的復(fù)原光譜精度不 隨地物類型變化． 通過對(duì)圖2(a)～2(d)中不同切趾函數(shù)作用下 的復(fù)原光譜的比較分析發(fā)現(xiàn)，當(dāng)最大光程差一定時(shí)， 矩形函數(shù)作用下的復(fù)原光譜相對(duì)誤差最大，其它切 趾函數(shù)作用后的復(fù)原光譜精度明顯提高，尤其是當(dāng) 最大光程差L=0．1cm和L=0．4cm時(shí)，這表明切趾 函數(shù)處理是光譜復(fù)原過程中不可忽略的環(huán)節(jié)：其中， 相對(duì)于其他切趾函數(shù)，Hanning函數(shù)作用后的復(fù)原 光譜在各個(gè)最大光程差下相對(duì)誤差都較�。�此外，同 一最大光程差下，各個(gè)切趾函數(shù)作用下的復(fù)原光譜 相對(duì)誤差都隨著波長(zhǎng)增大而增大，這是隨著波長(zhǎng)變 大波長(zhǎng)分辨率變低的一種表現(xiàn)． 從圖2(a)至圖2(d)，復(fù)原光譜對(duì)應(yīng)的最大光 程差逐漸變大，比較不同最大光程差、同一切趾函數(shù) 作用下的復(fù)原光譜發(fā)現(xiàn)，隨著最大光程差增大，相對(duì) 誤差逐漸變小．圖2(a)中HJY20．1．A對(duì)應(yīng)的復(fù)原光 譜誤差很大，圖2(b)中的復(fù)原光譜精度相對(duì)圖2 (a)有了顯著提高，而當(dāng)最大光程差L=0．4cm時(shí)， 復(fù)原光譜則已經(jīng)非常貼近真實(shí)光譜． 可見，雖然切趾函數(shù)的作用可以提高復(fù)原光 譜精度，但干涉成像光譜儀的最大光程差值的影 響更大．而對(duì)于空間調(diào)制型干涉成像光譜儀，由于 受到硬件工藝限制，通常最大光程差很難設(shè)置的 很大．因此，在最大光程差較小的情況下，如何提 高復(fù)原光譜精度，對(duì)于空間調(diào)制型干涉成像光譜 儀尤其重要． 一 表1 4種切趾函數(shù)及對(duì)應(yīng)的儀器線型函數(shù) Table 1 Four diferent apodization functions and corresponding ILS functions 其中：一L≤6≤ L， =2叮r( l一 )L，=2：2叮r( l一 )L一叮r， =2叮r( l一 )L+叮r，Z4=2叮r( l一 )L一2叮r， f=2叮r(GrI一 )L+2叮r 0．4 0．3 0．2 0．1 0 — 0．1 ． 0．2 ． 0．3 — 0．4 — 0．5 4 0．4 0．3 0．2 0．1 0 一O．1 ．o．2 — 0．3 4 ．{ 三角函數(shù) ． 450 550 650 750 850 Z／nm Z／nm (a) 醐八 -Un．I0)01 、／、 ／ l l ．n⋯02 L Z／nm (c) Z／nm Z／nm (d) A／nm ／nm 圖2 四種不同最大光程差設(shè)置下的干涉成像光譜儀數(shù)據(jù)經(jīng)過四種不同切趾函數(shù)作用后的復(fù)原光譜相對(duì)誤差圖

3 光譜復(fù)原算法改進(jìn)與分析 3．1 光譜復(fù)原算法改進(jìn) 切趾函數(shù)的作用在一定程度上提高了復(fù)原光 譜精度，當(dāng)L=0．1cm時(shí)，950nm處的光譜分辨率 已經(jīng)高達(dá)0．55nm，但從圖2可以看出，復(fù)原光譜相 對(duì)誤差仍然很大．同時(shí)，HJY20．1．A經(jīng)切趾函數(shù)作 用后的復(fù)原光譜誤差仍然很大．因此如何提高復(fù) 原光譜精度對(duì)于HJY20．1．A 的應(yīng)用顯得非常 重要． 干涉成像光譜儀干涉信號(hào)的突然截?cái)嗍沟脧?fù)原 光譜能量外泄，這導(dǎo)致波數(shù) 處的復(fù)原光譜是不僅 含有該波數(shù)處的光譜信息，同時(shí)還疊加了周圍波譜 信息．因此采用基于儀器線型函數(shù)標(biāo)準(zhǔn)化的光譜復(fù) 原方法，以提高復(fù)原光譜的精度．該方法得到的復(fù)原 光譜B_std (17")為： B — std ( )= m “ ΣB(o-) min lLs Σ兒s O’min 對(duì)24種典型地物，根據(jù)公式(4)求得對(duì)應(yīng)4種 不同最大光程差、4種不同切趾函數(shù)的復(fù)原光譜．這 些不同條件下的24種典型地物復(fù)原光譜對(duì)應(yīng)的相 對(duì)誤差采用與公式(3)相同方法求得，相對(duì)誤差圖 如圖3所示． 對(duì)應(yīng)比較圖2與圖3中的各幅圖，發(fā)現(xiàn)同一最 大光程差、同一切趾函數(shù)作用下的復(fù)原光譜，經(jīng)過儀 器線型函數(shù)標(biāo)準(zhǔn)化后的復(fù)原光譜相對(duì)誤差明顯變 小，同時(shí)還削弱了復(fù)原光譜相對(duì)誤差隨波長(zhǎng)的變化 性．HJY20—1-A對(duì)應(yīng)的復(fù)原光譜采用改進(jìn)算法后精 度顯著提高． 此外，比較圖3中同一最大光程差下不同切趾 0 0 0 ． 0 ． 0 - 0 0 0 舢 0 0 0 。0 。0 。0 。O ×10。5 1 0．5 0 。0．5 一l 450 550 650 750 850 950 E／nm 圖3 四種不同最大光程差設(shè)置的24種典型地物的干涉成像光譜儀數(shù)據(jù)，經(jīng)過四種不同切趾函數(shù)作用以及儀器線型函數(shù) 標(biāo)準(zhǔn)化后的復(fù)原光譜相對(duì)誤差圖(a)L=0．0069cm (b)L=0．05cm (C)L=0．1cm(d)L=0。4cm Fig 3 Relative elTors of reconstructed spec~a calculated by the improved algorithm with diferent MPD and diferent apodization functions 函數(shù)作用后的復(fù)原光譜可見，切趾函數(shù)作用后的誤 差明顯降低，其中，Hanning函數(shù)作用后的復(fù)原光譜 在不同最大光程差時(shí)相對(duì)誤差都較小．但對(duì)于 HJY20—1一A，其它切趾函數(shù)作用后的復(fù)原光譜精度 卻稍微下降． 由上述分析可知，對(duì)干涉成像光譜儀進(jìn)行光譜 復(fù)原時(shí)，儀器線型函數(shù)標(biāo)準(zhǔn)化能有效提高復(fù)原光譜 的精度，Hanning函數(shù)作用下的復(fù)原光譜在各個(gè)最 大光程差時(shí)仍有很好的精度。

3．2 植被指數(shù)應(yīng)用分析 植被遙感是HJY20—1一A的一個(gè)重要應(yīng)用方向． 本文以三種典型植被指數(shù)為例，針對(duì)24種典型地 物，比較分析HJY20—1一A探測(cè)數(shù)據(jù)在不同切趾函數(shù) 作用時(shí)，采用改進(jìn)算法后復(fù)原光譜的植被指數(shù)反演 精度，進(jìn)而分析討論不同切趾函數(shù)在HJY20—1一A具 體應(yīng)用中的影響．這3種典型植被指數(shù)是： 歸一化植被指數(shù) ： NDVI = R76 2- R6s0 葉綠素吸收反射率指數(shù) ： C RI = R7o0 l 670a + R670+ b R670／ra—— ：——+————1一 ， 其中，a=(R7o0一R550)／150，b=R550—550a。 修正的三角植被指數(shù)2 J： MTVI2= ： ： ： ： 。f7) √(2 8o0+1) 一(6R8o0—5~／ 670—0。5) 利用3．1節(jié)中HJY20—1一A對(duì)應(yīng)4種切趾函數(shù)作 用下的復(fù)原光譜，根據(jù)公式(5)一(7)以及相對(duì)誤差 公式，求得植被指數(shù)結(jié)果以及相對(duì)誤差，其中，根據(jù) 干涉成像光譜儀HJY20—1一A的波段設(shè)置，采用 668nm、681nm、700nm、762nm替代公式中的670nm、 681nm、699nm和765nm．分析植被光譜曲線特性可 知，這種鄰近波段替代對(duì)植被指數(shù)求解影響很�。�植 室 Sample Number Sample Number (c) 圖4 不同切趾函數(shù)作用下HJY20．1一A的24種典型地物復(fù)原光譜的三種植被指數(shù)結(jié)果以及相對(duì)誤差圖．左邊是植被 指數(shù)結(jié)果，右邊為相對(duì)誤差圖(a)NDVI(b)CARI(C)MTVI2 被指數(shù)及其相對(duì)誤差結(jié)果如圖4所示，圖中橫坐標(biāo) 對(duì)應(yīng)于圖1中的24種典型地物． 由圖4可見，圖4(a)和圖4(c)中不同切趾函數(shù) 對(duì)應(yīng)的植被指數(shù)精度相當(dāng)，而圖4(b)中矩形函數(shù)對(duì) 應(yīng)的植被指數(shù)精度最高，這與圖3(a)中矩形函數(shù)對(duì) 應(yīng)的復(fù)原光譜精度最高相符．總體而言，采用改進(jìn)算 法后求得的四種不同切趾函數(shù)作用下的HJY20—1一A 的復(fù)原光譜對(duì)應(yīng)的植被指數(shù)精度較高，尤其前l(fā)6種 地物(均為典型植物)的誤差非常�。�植被指數(shù)的精度 依賴于復(fù)原光譜的精度以及相應(yīng)的植被指數(shù)公式，改 進(jìn)的算法顯著提高了HJY20-1-A復(fù)原光譜精度，進(jìn)而 提高了HJY20—1一A在植被方面的應(yīng)用潛力．

4 結(jié)論 比較分析24種典型地物對(duì)應(yīng)四種最大光程差、 四種切趾函數(shù)作用下的復(fù)原光譜，以及基于儀器線 型函數(shù)標(biāo)準(zhǔn)化的光譜復(fù)原改進(jìn)算法后復(fù)原光譜，得到以下結(jié)論：
(1)已有的復(fù)原光譜算法中，同一切趾函數(shù)作 用下的復(fù)原光譜精度不隨地物類型變化．各個(gè)切趾 函數(shù)作用下的復(fù)原光譜相對(duì)誤差總體都隨著波長(zhǎng)增 大而增大，這是隨著波長(zhǎng)變大波長(zhǎng)分辨率變低的一 種表現(xiàn)．
(2)已有的復(fù)原光譜算法中，切趾函數(shù)的作用 雖然可以提高復(fù)原光譜精度，但復(fù)原光譜精度受干 涉成像光譜儀最大光程差值的影響更大．而對(duì)于空 間調(diào)制型干涉成像光譜儀，由于受到硬件工藝限制， 通常最大光程差很難設(shè)置的很大．因此，在最大光程 差較小的情況下，如何提高復(fù)原光譜精度，對(duì)于空間 調(diào)制型干涉成像光譜儀尤其重要．
(3)使用基于儀器線型函數(shù)標(biāo)準(zhǔn)化的改進(jìn)算法 后，同一最大光程差、同一切趾函數(shù)作用下的復(fù)原 光譜精度明顯提高，同時(shí)還削弱了復(fù)原光譜相對(duì)誤