鈕立明,蒙繼華,吳炳方,陳雪洋,杜 鑫,張飛飛

(1.北京天宏金晴信息技術(shù)有限公司,北京 100101;2.中國科學(xué)院遙感應(yīng)用研究所,北京 100101;3.遙感科學(xué)國家重點(diǎn)實(shí)驗(yàn)室,北京 100101)

HJ-1A星 HSI數(shù)據(jù) 2級(jí)產(chǎn)品處理流程研究

鈕立明1,蒙繼華2,3,吳炳方2,陳雪洋2,杜 鑫2,張飛飛2

(1.北京天宏金晴信息技術(shù)有限公司,北京 100101;2.中國科學(xué)院遙感應(yīng)用研究所,北京 100101;3.遙感科學(xué)國家重點(diǎn)實(shí)驗(yàn)室,北京 100101)

研究了 HJ-1A星 HSI數(shù)據(jù) 2級(jí)產(chǎn)品的數(shù)據(jù)預(yù)處理流程及相關(guān)算法,包括絕對(duì)輻亮度值轉(zhuǎn)換、條紋去除、大氣校正及幾何糾正,得到了具有精確地理位置信息的地表光譜反射率圖像;基于相同位置同期的一景 Hyperion數(shù)據(jù)標(biāo)準(zhǔn)化處理流程得到的地表反射率,進(jìn)行了 HSI數(shù)據(jù)的光譜模擬,并將模擬的地表反射率與真實(shí) HSI數(shù)據(jù)的地表反射率進(jìn)行對(duì)比分析,評(píng)價(jià)本研究所采用預(yù)處理流程的有效性。研究結(jié)果表明:應(yīng)用該流程處理的 HSI數(shù)據(jù)地表反射率與模擬 HSI的地表反射率平均相關(guān)系數(shù)為 0.947,標(biāo)準(zhǔn)差為 0.017,兩種數(shù)據(jù)的光譜反射率變化表現(xiàn)出高度一致性;各波段模擬與真實(shí) HSI地表反射率差值的平均值和標(biāo)準(zhǔn)差都接近 0值,說明本研究所發(fā)展的數(shù)據(jù)處理流程可以為 HSI數(shù)據(jù)定量分析與應(yīng)用提供參考。

HJ-1A HSI;標(biāo)準(zhǔn)化預(yù)處理流程;Hyperion;對(duì)比分析

0 引言

2008年 9月 6日發(fā)射的環(huán)境和災(zāi)害監(jiān)測(cè)預(yù)報(bào)小衛(wèi)星星座 A星 (以下簡稱 HJ-1A星)上除 CCD相機(jī)外,還搭載了高光譜成像儀 (Hyperspectral Imaging Radiometer,HSI)。HSI采用干涉成像光譜技術(shù),其光譜分辨率最高達(dá) 2.08 nm,相對(duì)于現(xiàn)有的星載成像高光譜傳感器,具有較高的光譜分辨率,可在對(duì)地觀測(cè)的同時(shí)獲取眾多連續(xù)波段的地物光譜圖像,達(dá)到從空間直接識(shí)別地球表面物體的目的[1]。

HJ-1A星 HSI數(shù)據(jù)的出現(xiàn)填補(bǔ)了目前國產(chǎn)數(shù)據(jù)在高光譜傳感器領(lǐng)域的空白,其 2級(jí)產(chǎn)品正在向用戶免費(fèi)分發(fā),具有極大的應(yīng)用前景。由于 HSI傳感器投入運(yùn)行時(shí)間不長,其數(shù)據(jù)質(zhì)量還存在一些問題:如各波段上存在不同程度的條紋噪聲;圖像的光譜特征曲線有一定不合理的抖動(dòng)等。針對(duì)這些問題,目前仍缺少有效的數(shù)據(jù)預(yù)處理流程,這極大地限制了 HSI數(shù)據(jù)的推廣應(yīng)用。因此,有必要針對(duì)數(shù)據(jù)產(chǎn)品存在的問題深入研究完整的預(yù)處理流程,得到較為真實(shí)的地表光譜反射率圖像,為進(jìn)一步的高光譜定量遙感應(yīng)用提供準(zhǔn)確的數(shù)據(jù)支持。

本文研究了 HJ-1A星 HSI數(shù)據(jù) 2級(jí)產(chǎn)品處理流程,著重闡述預(yù)處理的內(nèi)容和方法,并將處理結(jié)果與同期的一景 Hyperion數(shù)據(jù)進(jìn)行對(duì)比分析,評(píng)價(jià)HSI數(shù)據(jù)預(yù)處理的效果。

1 HSI數(shù)據(jù)產(chǎn)品與研究區(qū)概況

1.1 HSI數(shù)據(jù)產(chǎn)品

HSI數(shù)據(jù)產(chǎn)品分 0～5共 6個(gè)級(jí)別,包括原始數(shù)據(jù)產(chǎn)品、輻射校正產(chǎn)品、系統(tǒng)幾何糾正產(chǎn)品、幾何精糾正產(chǎn)品、高程校正產(chǎn)品和標(biāo)準(zhǔn)鑲嵌圖像產(chǎn)品。普通用戶通過訪問中國資源衛(wèi)星應(yīng)用中心的網(wǎng)站便可以免費(fèi)獲得 HSI數(shù)據(jù)的 2級(jí)產(chǎn)品,即系統(tǒng)幾何糾正產(chǎn)品。本研究即在 HSI2級(jí)產(chǎn)品基礎(chǔ)上進(jìn)行的。

HSI2級(jí)產(chǎn)品經(jīng)過輻射校正和系統(tǒng)幾何糾正。幾何定位精度在 1 km左右,誤差分布比較均勻,糾正后的影像映射到指定的地圖投影坐標(biāo)系下[1],其特性如表 1所示。

表 1 HSI2級(jí)產(chǎn)品特性說明Tab.1 Introduction of HSI level2 data product

1.2 研究區(qū)概況

以吉林省洮南縣為中心的 Hyperion數(shù)據(jù)與環(huán)境星 HSI數(shù)據(jù)共同覆蓋范圍為實(shí)驗(yàn)區(qū),實(shí)驗(yàn)區(qū)中心位置為東經(jīng) 121°44′31″、北緯 45°53′54″,其空間位置及遙感數(shù)據(jù)如圖 1所示。

圖 1遙感數(shù)據(jù)為 HSI第 105、70、40波段的假彩色合成圖像。實(shí)驗(yàn)區(qū)由北向南依次覆蓋了烏蘭浩特市南部、科爾沁右翼前旗東南部、洮南市西部和突泉縣東北部的小部分地區(qū)。

2 HSI數(shù)據(jù)預(yù)處理

HSI數(shù)據(jù)預(yù)處理流程如圖 2所示。

圖 2 HSI數(shù)據(jù)預(yù)處理流程Fig.2 Preprocessing flow for H IS data

HSI數(shù)據(jù)的 2級(jí)產(chǎn)品采用 HDF5格式。目前大多數(shù)遙感圖像處理軟件都無法對(duì)其進(jìn)行直接操作,因此需要先通過中國資源衛(wèi)星中心提供的軟件將其轉(zhuǎn)換成 GEOTIFF格式。該級(jí)別產(chǎn)品經(jīng)過了輻射校正及壞像元恢復(fù)等圖像質(zhì)量檢查等過程,已經(jīng)去除了圖像中的壞線,但仍有部分波段有條紋噪聲等異常存在,在對(duì)圖像數(shù)據(jù)進(jìn)行應(yīng)用之前,必須將這些異常加以校正。

2.1 絕對(duì)輻亮度值的轉(zhuǎn)換

HSI2級(jí)產(chǎn)品數(shù)據(jù)集為輻亮度產(chǎn)品,以無符號(hào)整型數(shù)據(jù)記錄,數(shù)值范圍為 0～65 535。實(shí)際上,由于地物的輻射亮度值很小,在產(chǎn)品生成時(shí)乘了一個(gè)擴(kuò)大系數(shù) 100,因此在應(yīng)用時(shí),只需要將各波段數(shù)值除以系數(shù)100,即可以得到單位為 W/(m2·sr·μm)的絕對(duì)輻射亮度值圖像。

2.2 條紋的去除

在 HSI數(shù)據(jù) 2級(jí)產(chǎn)品的部分波段上可以看到明顯的條紋,這些像元灰度值比較小,與周圍有明顯差異,條紋的存在嚴(yán)重影響了圖像的質(zhì)量和應(yīng)用。波長較短的波段條紋較清楚,與周邊差異明顯,但數(shù)量較少;波長較長的波段條紋較模糊,與周邊差異不明顯,但數(shù)量較多。圖 3為波長 479.6 nm和波長529.6 nm的原始數(shù)據(jù)產(chǎn)品,從圖中可以明顯地看到條紋。

圖 3 實(shí)驗(yàn)區(qū) HSI第 10波段 (479.6 nm,左圖)和第 30波段 (529.6 nm,右圖)原始圖像Fig.3 The original image of band 10(479.6 nm,left)and band 30(529.6 nm,right)

采用“全局去條紋”方法[2]去除影像中的垂直條紋。對(duì)于 HSI數(shù)據(jù) 2級(jí)產(chǎn)品,由于已經(jīng)過系統(tǒng)幾何糾正,圖像發(fā)生了旋轉(zhuǎn),圖像上的條紋也不再呈垂直分布,因此在進(jìn)行去條紋之前需要先將圖像進(jìn)行一定的旋轉(zhuǎn),使條紋能沿垂直方向分布。

2.2.1 圖像旋轉(zhuǎn)

快速、自動(dòng)地確定圖像旋轉(zhuǎn)角度,減少人工參與,可以在很大程度上保證處理方法的精確性和普適性,從而大大提高圖像處理效率。通過確定圖像兩個(gè)頂點(diǎn)坐標(biāo),并求得兩點(diǎn)所在直線與法線夾角的正切值,進(jìn)而確定該夾角為圖像旋轉(zhuǎn)的角度,如圖 4所示。

圖 4 圖像旋轉(zhuǎn)角度計(jì)算Fig.4 Calculation of image rotation angle

利用圖像背景值為 0這一特點(diǎn),通過求每列像元值的和,確定第一個(gè)和不等于 0的列號(hào)作為 x1,將這一列中第一個(gè)值不等于 0的行號(hào)作為 y1。用類似的方法確定 x2和 y2。通過多次試驗(yàn),確定圖像旋轉(zhuǎn)角度α的計(jì)算公式為

將各波段數(shù)據(jù)旋轉(zhuǎn)角度α后,圖像上的條紋基本上均沿垂直方向分布,這樣即可以采用“全局去條紋法”進(jìn)行條紋去除。

2.2.2 垂直條紋去除

針對(duì)條紋垂直分布的特征,采用全局去條紋法[2,6],通過像元的列平均值、標(biāo)準(zhǔn)差和波段平均值、標(biāo)準(zhǔn)差之間的差異對(duì)像元進(jìn)行分波段線性化修正,消除垂直條紋的影響,即

式中,DNi,j,k和 DN′i,j,k分別為原始和消除垂直條紋影響后的像元值。增益 gi,k和偏移量 bi,k計(jì)算公式為

式中,Stdvk和分別為第 k波段像元值的標(biāo)準(zhǔn)差和平均值;Stdvi,k和分別為第 k波段第 i列像元值的標(biāo)準(zhǔn)差和平均值。

垂直條紋去除前后效果對(duì)比如圖 5所示?？梢钥闯?原有的垂直條紋已經(jīng)被去除掉,處理后圖像的灰度值在空間上過渡平滑。

圖 5 第 30波段垂直條紋去除前 (左)后 (右)效果對(duì)比Fig.5 Comparison between before(left)and after(right)destripe in band 30

將去除條紋后的圖像按之前計(jì)算的角度進(jìn)行反向旋轉(zhuǎn),即完成圖像的條紋去除。

2.2.3 條紋去除效果評(píng)價(jià)

通過對(duì)原始圖像和去除條紋后的圖像進(jìn)行最小噪聲分離變換 (Minimum Noise Fraction,MNF)來識(shí)別和分離數(shù)據(jù)中的噪聲[7]。經(jīng)過變換,高光譜數(shù)據(jù)中的噪聲基本被分離出來了,通過對(duì)比去條紋前后MNF特征圖像,可以非常明顯地看出條紋去除的效果。圖 6比較了去條紋前后的MNF特征圖像第10波段。

圖 6 去除條紋前 (左 )后 (右 )第 10波段MNF特征圖像比較Fig.6 Comparison between before(lift)and after(right)destripe in M NF(band 10)

2.3 大氣校正

2.3.1 光譜吸收特征波段的確定

本研究采用軟件 FLAASH(Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes)進(jìn)行大氣校正。由于軟件中并沒有與 HSI數(shù)據(jù)相對(duì)應(yīng)的傳感器類型,因此需要自定義一個(gè)描述光譜吸收特征波段的文件 (Spectrograph Definition File)用以進(jìn)行大氣校正[8]。

大氣頂層的輻射亮度具有隨電磁輻射波長增加而減弱的趨勢(shì),這是因?yàn)?在傳感器端的太陽輻射是乘性因子。通過計(jì)算各波段大氣頂層反射率,可以去除太陽輻射的影響;由大氣頂層反射率曲線可以確定 820 nm附近的水汽吸收谷波段位置。將相關(guān)的波長、FWHM以及波段號(hào)等信息輸入 FLAASH模型可以計(jì)算水汽透過率。由于 HSI數(shù)據(jù)的波長范圍最長只到 950 nm,在 FLAASH中無法進(jìn)行氣溶膠反演,因此只能根據(jù)影像情況設(shè)定初始能見度為40 km。大氣頂層反射率計(jì)算公式為

式中,ρp為大氣頂層反射率;Lλ為大氣頂層輻射亮度;d為日地距離;ESUNλ為各波段大氣頂層的平均太陽輻照度;θs為太陽天頂角。

2.3.2 大氣校正

首先,在 FLAASH軟件中輸入相關(guān)參數(shù),包括圖像成像時(shí)間、地理位置、海拔、大氣模式以及氣溶膠模式,定義吸收特征波段描述文件,然后進(jìn)行大氣校正。典型地物 (植被)大氣校正前后光譜曲線對(duì)比如圖 7所示。

圖 7 HSI數(shù)據(jù)大氣校正前后植被反射光譜曲線Fig.7 Spectral curves before and after the atmospheric correction

從圖 7可以看出,經(jīng)過大氣校正,近紅外波段范圍內(nèi)由于水汽吸收而形成的反射率低谷明顯消失,植被反射光譜曲線的特征基本保持。

2.3.3 光譜曲線鋸齒的去除

經(jīng)過大氣校正后的圖像,部分像元的光譜特征曲線存在光譜抖動(dòng),這是由于在輻射定標(biāo)等數(shù)據(jù)預(yù)處理過程中不可避免地會(huì)引入一些誤差和噪音,因此需要對(duì)大氣校正后的圖像進(jìn)行處理。

通過MNF變換將數(shù)據(jù)分離成包含有用信息且連貫的MNF特征波段和以噪聲為主的MNF特征波段,將含有信息量大的前 10個(gè)波段進(jìn)行MNF反變換,得到光譜優(yōu)化后的圖像。這種方法可以有效去除光譜抖動(dòng),而且不會(huì)降低圖像的空間分辨率[6]。圖 8為光譜優(yōu)化前后 HSI數(shù)據(jù)光譜特征曲線的比較。從圖中可以看到,曲線上不正常的光譜抖動(dòng)和鋸齒基本上都被消除了。

圖 8 HSI數(shù)據(jù)光譜優(yōu)化前后植被反射光譜曲線Fig.8 Spectral curves before and after the optimization

2.4 幾何糾正

以地形圖或經(jīng)過幾何精糾正的圖像為參考,通過人工在一景圖像上選取 25～30個(gè)控制點(diǎn),利用二次多項(xiàng)式糾正模型和最鄰近像元重采樣方法對(duì)大氣校正后的 HSI數(shù)據(jù)進(jìn)行幾何糾正,誤差控制在一個(gè)像元以內(nèi)。

2.5 對(duì)比分析

Hyperion作為第一個(gè)星載民用成像光譜儀,已運(yùn)行近 7 a。經(jīng)過多年的研究,其數(shù)據(jù)的預(yù)處理流程已十分成熟和完善,得到的地表反射率也比較可靠。本研究采用相同位置且同期的一景 Hyperion數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化預(yù)處理[3,9],得到地表反射率數(shù)據(jù),并以此作為參照,對(duì)比分析 HSI數(shù)據(jù)預(yù)處理結(jié)果的可靠性,如圖 9所示。

圖 9 HSI(左)與 Hyperion(右)地表反射率數(shù)據(jù)Fig.9 HSI(lift)and Hyper ion data(right)for comparison

Hyperion數(shù)據(jù)的空間分辨率為 30m,光譜分辨率在 10 nm左右,其波段設(shè)置與 HSI數(shù)據(jù)并不相同。為了使被對(duì)比的數(shù)據(jù)與 HSI數(shù)據(jù)在空間和光譜上具有一致性,本研究將 Hyperion數(shù)據(jù)的預(yù)處理結(jié)果與HSI數(shù)據(jù)的地表反射率結(jié)果進(jìn)行了幾何配準(zhǔn),誤差控制在 1個(gè)像元以內(nèi),并降空間分辨率至 100m,使之空間尺度保持一致;將 Hyperion數(shù)據(jù)中與 HSI數(shù)據(jù)相同光譜范圍的波段提取出來,并插值成光譜分辨率為 1 nm的數(shù)據(jù),結(jié)合高斯函數(shù)模擬出 HSI光譜響應(yīng)函數(shù),得到模擬的 HSI數(shù)據(jù),作為參照。

2.5.1 地表反射率一致性分析

對(duì) 450～950 nm光譜范圍內(nèi)的 HSI與 Hyperion模擬的光譜反射率進(jìn)行逐像元相關(guān)分析,其結(jié)果如圖 10所示。

圖 10 逐像元的地表反射率相關(guān)系數(shù)及其直方圖Fig.10 Correlation coefficient of two images and its histogram

可以看出,在整個(gè)研究區(qū),HSI與 Hyperion模擬數(shù)據(jù)在整個(gè)波段范圍內(nèi)表現(xiàn)為極顯著相關(guān)。從相關(guān)系數(shù)及其空間差異上來看,兩種數(shù)據(jù)的光譜反射率變化表現(xiàn)得非常一致。

2.5.2 差值分析

對(duì) Hyperion模擬與真實(shí)環(huán)境星 HSI的地表反射率進(jìn)行差值計(jì)算,各波段差值圖的相關(guān)統(tǒng)計(jì)信息變化趨勢(shì)如圖 11所示。

圖 11 Hyperion與 HSI地表反射率差值統(tǒng)計(jì)Fig.11 Statistics for the differential image of reflection between Hyper ion and HSI

從圖 11來看,Hyperion模擬與真實(shí) HSI地表反射率差值的均值基本在 0附近,表明兩種數(shù)據(jù)中大部分區(qū)域的地表反射率大致相同;標(biāo)準(zhǔn)差均遠(yuǎn)小于其值域分布范圍,說明其分布集中、波動(dòng)程度小。各波段的差值大體服從正態(tài)分布。以第 10波段為例(圖 12),差值最值約為 0.1,最值像元主要出現(xiàn)在不同覆被類型的邊緣、破碎地區(qū)和特殊土地覆蓋區(qū)域(如水體),這主要是由于兩種數(shù)據(jù)在獲取時(shí)的空間尺度差異造成的,這種差異無法通過像元分辨率的變化得到完全消除,并在地表覆被破碎地區(qū)表現(xiàn)的更為突出。

圖 12 第 10波段逐像元的地表反射率差值及其直方圖Fig.12 Differential image and its histogram of reflection between Hyper ion and HSI

從相關(guān)性及差值統(tǒng)計(jì)分析結(jié)果來看,經(jīng)過標(biāo)準(zhǔn)化預(yù)處理得到的環(huán)境星 HSI地表反射率與 Hyperion模擬的地表反射率一致性較高,說明針對(duì) HSI的數(shù)據(jù)預(yù)處理方法是有效的。

3 結(jié)論

高光譜數(shù)據(jù)在反演植被生化參數(shù)等方面有著其他傳感器所不能比擬的優(yōu)勢(shì),是進(jìn)行精準(zhǔn)農(nóng)業(yè)遙感的重要數(shù)據(jù)源,而預(yù)處理結(jié)果的好壞直接關(guān)系著遙感分析、應(yīng)用的精度。

(1)針對(duì) HSI數(shù)據(jù) 2級(jí)產(chǎn)品,設(shè)計(jì)和提出了完整的圖像預(yù)處理流程,得到了具有精確地理位置信息的地表光譜反射率圖像。將相同位置且同期Hyperion圖像得到的地表反射率結(jié)果與之對(duì)比,認(rèn)為該流程可為 HSI數(shù)據(jù)的定量應(yīng)用提供參考。

(2)受 HSI傳感器波長范圍的限制,目前仍無法準(zhǔn)確地計(jì)算氣溶膠光學(xué)厚度及大氣水汽含量,今后還需要開展針對(duì)環(huán)境星 HSI傳感器大氣校正算法的研究,以提高大氣校正精度。

(3)由于缺乏相關(guān)的參數(shù),暫時(shí)還無法進(jìn)行 HSI傳感器 Sm ile效應(yīng)的糾正,Sm ile效應(yīng)會(huì)影響紅邊位置等高光譜參數(shù)提取的準(zhǔn)確性,今后將在這方面開展相關(guān)的研究工作。

致謝:感謝中國資源衛(wèi)星應(yīng)用中心提供了HJ-1A HSI數(shù)據(jù)。

[1] 中國資源衛(wèi)星應(yīng)用中心.環(huán)境和災(zāi)害監(jiān)測(cè)預(yù)報(bào)小衛(wèi)星星座 AB星用戶手冊(cè)[Z].北京:中國資源衛(wèi)星應(yīng)用中心,2009.

[2] 譚炳香,李增元,陳爾學(xué),等.EO-1 Hyperion高光譜數(shù)據(jù)的預(yù)處理[J].遙感信息,2005(6):36-41.

[3] 張 東,張 鷹,李 歡.海岸帶星載高光譜遙感影像預(yù)處理方法[J].海洋科學(xué)進(jìn)展,2009,27(1):92-97.

[4] 李傳榮,賈媛媛,胡 堅(jiān),等.HJ-1光學(xué)衛(wèi)星遙感應(yīng)用前景分析[J].國土資源遙感,2008(3):1-3.

[5] 李石華,角媛梅.環(huán)境與災(zāi)害監(jiān)測(cè)預(yù)報(bào)小衛(wèi)星 A星 CCD影像質(zhì)量評(píng)價(jià)[J].紅外技術(shù),2009,31(9):167-172.

[6] Datt B,McVicar TR,Van Niel TG,et al.Pre-processing EO-1 Hyperion Hyperspectral Data to Support the Application of Agricultural Indices[J].IEEE Transactions on Geoscience and Remote Sensing,2003,41(6):1246-1259.

[7] Boardman JW,Kruse FA.Automated Spectral Analysis:A Geo log-ical Example U sing AVIR IS Data[C]//Proceedings of 10 th Thematic Conference on Geologic Remote Sensing.Ann Arbor North Grapevine Mountains,Nevada,Environmental Research Institute of Michigan,1994,407-418.

[8] Research System,FLAASH,Boulder[M/OL].CO:Research System,Inc,2003.http://www.rsinc.com/envi/flaash.asp

[9] David G,Goodenough,Andrew D,et al.Processing Hyperion and ALI for Forest Classification[J].IEEE Transaction on Geoscience and Remote Sensing,2003,41(6):1321-1331.

[10]Chen CM,Hepner G F,Forster R R.Fusion of Hyperspectral and Radar Data U sing the IHS Transformation to Enhance Urban Surface Features[J].ISPRS Journal of Photogrammetry&Remote Sensing,2003,58:19-30.

(責(zé)任編輯:刁淑娟)

Research on Standard Preprocessing Flow for HJ-1A HSIL evel 2 Data Product

NIU Li-ming1,MENG Ji-hua2,3,WU Bing-fang2,CHEN Xue-yang2,DU Xin2,ZHANG Fei-fei2
(1.Beijing Wisewatch IT Company,Beijing 100101,China;2.Institute of Remote Sensing Applications,Chinese Academy of Sciences,Beijing 100101,China;3.State Key Laboratory of Remote Sensing Science,Jointly Sponsored by the Institute of Remote Sensing Applications of Chinese Academy of Sciences and Beijing Normal University,Beijing 100101,China)

To deal with the Level2 HSI data from the newly-launched HJ-1A satellite,this paper introduced in details the entire flow and relevant algorithms for data preprocessing.The introduction includes calibration,vertical stripes elimination,and atmospheric correction geometric correction.Standard spectral reflectance products with precise geo-locations were produced.Spectral reflectance data from EO-1 Hyperion of c lose dates was used to simulate the band reflectance of HJ-1A HSI.Comparisons of spectral reflectance data between simulated and actual HJ-1A HSI were made to validate the effect of the data preprocessing.The average correlation coefficient of spectral reflectance between actual and simulated HJ-1A HSI is0.947 with its standard deviation being 0.017,suggesting a high consistency.The mean and standard deviation of differential bands between real and simulated HJ-1A HSI are close to 0.The result show s that the reflectance from HJ-1A HSI is consistent with that of simulated data from Hyperion,and hence the data processing flow could provide necessary support for quantitative use of HJ-1A HSI data.

HJ-1A HSI;Standard preprocessing flow;Hyperion;Comparative analysis

蒙繼華 (1977-),男,主要從事農(nóng)作物遙感監(jiān)測(cè)及其在精準(zhǔn)農(nóng)業(yè)領(lǐng)域中的應(yīng)用研究。m engjh@irsa.ac.cn。

TP 75

A

1001-070X(2011)01-0077-06

2010-03-12;

2010-04-05

HJ-1衛(wèi)星數(shù)據(jù)應(yīng)用研究專題(編號(hào):2009A 01A 1000、2009A 01A 0900);遙感科學(xué)國家重點(diǎn)實(shí)驗(yàn)室青年人才項(xiàng)目;中國科學(xué)院青年人才前沿領(lǐng)域?qū)ｍ?xiàng)項(xiàng)目(編號(hào):08S01700CX);中國科學(xué)院知識(shí)創(chuàng)新工程重大項(xiàng)目 (編號(hào):KSCX-YW-09-01);國家青年自然科學(xué)基金項(xiàng)目(編號(hào):NSFC40801144)。

鈕立明 (1984-),男,工程師,從事農(nóng)情遙感監(jiān)測(cè)及作物參數(shù)遙感反演技術(shù)研究。