霍紅元，周 萍，倪卓婭

(中國(guó)地質(zhì)大學(xué)(北京)地球科學(xué)與資源學(xué)院，北京 100083)

不同載荷指標(biāo)的模擬HyMap數(shù)據(jù)質(zhì)量評(píng)價(jià)

霍紅元，周 萍，倪卓婭

(中國(guó)地質(zhì)大學(xué)(北京)地球科學(xué)與資源學(xué)院，北京 100083)

從圖像預(yù)處理和地學(xué)應(yīng)用兩個(gè)角度對(duì)不同載荷指標(biāo)的模擬HyMap數(shù)據(jù)質(zhì)量進(jìn)行評(píng)價(jià)。在圖像預(yù)處理角度上，利用圖像的均方差異常、直方圖異常、數(shù)據(jù)相關(guān)異常等3個(gè)參量分別對(duì)模擬HyMap圖像質(zhì)量加以分析和評(píng)價(jià);在模擬數(shù)據(jù)的地學(xué)應(yīng)用角度上，利用模擬的HyMap數(shù)據(jù)進(jìn)行礦化蝕變信息提取和礦物填圖。通過對(duì)提取信息的種類、數(shù)量及礦物填圖精度的分析，獲得了對(duì)不同載荷指標(biāo)和不同尺度模擬HyMap數(shù)據(jù)質(zhì)量評(píng)價(jià)方法的許多重要認(rèn)識(shí)。

HyMap數(shù)據(jù)質(zhì)量評(píng)價(jià);圖像預(yù)處理;信息提取;礦物填圖;載荷指標(biāo)

0 引言

影響遙感圖像質(zhì)量的因素有很多，圖像評(píng)價(jià)參數(shù)因子也有多種，不同評(píng)價(jià)參數(shù)的原理不盡相同，在圖像質(zhì)量評(píng)價(jià)中所起的作用也必然不同，一般的評(píng)價(jià)因子有信噪比、灰度直方圖、圖像信息熵等［1－2］。本文著重分析影響模擬HyMap數(shù)據(jù)質(zhì)量的一個(gè)重要因素——噪聲。噪聲分析是遙感圖像質(zhì)量評(píng)價(jià)的重要一環(huán)。噪聲在理論上定義為“不可預(yù)測(cè)，只能用概率統(tǒng)計(jì)方法來認(rèn)識(shí)的隨機(jī)誤差”。因此，可將圖像噪聲看成是一個(gè)多維隨機(jī)過程，對(duì)噪聲的描述完全可以借用隨機(jī)過程及其概率分布函數(shù)和概率密度函數(shù)的方法［3］。但在很多情況下，這種描述方法很復(fù)雜，甚至不可能，而且在實(shí)際應(yīng)用中往往沒有必要，故通常用其數(shù)字特征(即均值、方差、相關(guān)函數(shù)等)來表征［4－5］。本文采用的模擬圖像質(zhì)量評(píng)價(jià)指標(biāo)為均方差異常、直方圖異常和數(shù)據(jù)相關(guān)異常;在HyMap模擬數(shù)據(jù)應(yīng)用試驗(yàn)中，根據(jù)所識(shí)別的蝕變礦物信息種類和數(shù)量及礦物填圖的精度等應(yīng)用結(jié)果對(duì)模擬數(shù)據(jù)質(zhì)量進(jìn)行分析和評(píng)價(jià)。

1 研究區(qū)與模擬數(shù)據(jù)

研究區(qū)位于新疆東天山土墩地區(qū)，地理坐標(biāo)為E94°4'～ 94°11'，N42°9'～ 42°13'。大地構(gòu)造位置處于準(zhǔn)噶爾板塊東南緣活動(dòng)帶，出露地層可分為結(jié)晶基底與蓋層。成礦帶經(jīng)歷多期構(gòu)造變形，巖漿活動(dòng)強(qiáng)烈，皺褶與斷裂系統(tǒng)極為發(fā)育，產(chǎn)出多種類型和規(guī)模不等的礦床，主要礦種為 Cu、Mo、Ni、Fe等多金屬。區(qū)內(nèi)的變質(zhì)作用主要為與構(gòu)造破碎及熱液蝕變有關(guān)的動(dòng)力變質(zhì)和熱液變質(zhì)作用。動(dòng)力變質(zhì)作用主要伴隨斷裂構(gòu)造形成，沿構(gòu)造線呈帶狀分布;熱液變質(zhì)作用主要產(chǎn)生部位為斷裂構(gòu)造帶和層間裂隙，其規(guī)模和強(qiáng)度取決于構(gòu)造的規(guī)模和性質(zhì)。蝕變作用主要形成黃鐵絹云巖化、絹云母化、高嶺土化、硅化、碳酸鹽化、綠泥石化和鉀化，同時(shí)伴有黃鐵礦化、方鉛礦化、孔雀石化及金礦化。

星載模擬HyMap數(shù)據(jù)是以澳大利亞HyVista公司生產(chǎn)的HyMap航空成像光譜數(shù)據(jù)為基礎(chǔ)、利用本文開發(fā)的模擬軟件對(duì)原始HyMap圖像加入一定比例的噪聲而模擬生成的，其不同載荷指標(biāo)與尺度如下:空間分辨率分別為10、15、50、75和100 m;光譜分辨率分別為5、15、20、25和30 nm;信噪比分別為50～100、250～300、350～400和 450～500(VNIR—SWIR)。

2 從圖像預(yù)處理角度評(píng)價(jià)模擬數(shù)據(jù)質(zhì)量

在各個(gè)不同的圖像預(yù)處理階段，經(jīng)處理的HyMap圖像數(shù)據(jù)具有不同的特征，對(duì)其進(jìn)行質(zhì)量評(píng)價(jià)時(shí)，需要選擇合適的評(píng)價(jià)參數(shù)［6－7］。HyMap圖像預(yù)處理包括飛機(jī)姿態(tài)的糾正，沿飛行和掃描方向的幾何糾正，圖像的輻射校正，圖像的配準(zhǔn)、均勻性處理和去噪聲等。

2.1 基于不同指標(biāo)與尺度的模擬HyMap圖像質(zhì)量分析

2.1.1 基于均方差異常的模擬圖像質(zhì)量分析

通常，圖像均方差大，則圖像亮度變化范圍大，圖像信息較豐富。但不能直接以均方差大小為標(biāo)準(zhǔn)來評(píng)判圖像質(zhì)量好壞。由于HyMap CCD鏡頭的故障，某個(gè)波段的圖像可能會(huì)有極高亮度值的像元(圖像亮斑)，則該波段圖像均方差比鄰近波段要大;如果高亮度的像元較少，也不會(huì)對(duì)整個(gè)波段圖像的均方差產(chǎn)生大的影響。而當(dāng)某波段圖像的均方差比相鄰波段的都小時(shí)，該波段的圖像質(zhì)量也不一定比其他波段差。如果對(duì)HyMap圖像的質(zhì)量進(jìn)行目視檢查，即可通過均方差圖像對(duì)其質(zhì)量做出大概的評(píng)價(jià)(圖1)。

圖1 HyMap圖像均方差圖解Fig.1 Diagrams of mean squared difference of HyMap image

(1)不同空間分辨率模擬圖像分析。在圖1(a)中，不同空間分辨率模擬圖像均方差的大小關(guān)系為S11＞S12＞S13＞S14＞S15(這里，原始 HyMap圖像的均方差記為S11，空間分辨率分別為15、45、60、75 m的模擬HyMap圖像的均方差分別記為S12、S13、S14、S15)?？梢钥闯?，15 m分辨率的模擬圖像質(zhì)量要明顯好于45、60、75 m分辨率的模擬圖像質(zhì)量;45 m分辨率模擬圖像的均方差稍高于60 m分辨率的模擬圖像的均方差，這2種類型模擬圖像的質(zhì)量相差不多;而分辨率為75 m的模擬圖像的質(zhì)量最差。

(2)不同光譜分辨率模擬圖像分析。記原始HyMap圖像的均方差為S21，光譜分辨率分別為5、15、20、25和30 nm的模擬HyMap圖像的均方差分別記為 S22、S23、S24、S25、S26，從圖 1(b)中可以看出，S21＞S22＞S23＞S24＞S25＞S26，表明不同光譜分辨率模擬圖像的質(zhì)量隨光譜分辨率的降低而降低。

(3)不同信噪比模擬圖像分析。從圖1(c)可以看出，不同信噪比的模擬HyMap圖像的均方差曲線基本重合。針對(duì)不同信噪比的模擬圖像，均方差異常度量指標(biāo)的可區(qū)分度不高，這可能是用模擬軟件模擬圖像時(shí)所增加的噪聲所致。

圖1(d)對(duì)原始HyMap圖像和不同空間分辨率、不同光譜分辨率、不同信噪比的模擬HyMap圖像均方差曲線進(jìn)行綜合對(duì)比，可以看出，S31＞S32≥S33＞S34(記原始HyMap圖像、不同光譜分辨率圖像、不同信噪比圖像、不同空間分辨率圖像的均方差分別為 S31、S32、S33、S34)。

從整體上看，原始HyMap圖像質(zhì)量要高于模擬圖像，不同光譜分辨率和不同信噪比的模擬HyMap圖像質(zhì)量相差無幾，都略高于空間分辨率低于30 m的模擬HyMap圖像;但對(duì)于15 m分辨率的模擬HyMap圖像而言，其均方差曲線值要高于其余載荷指標(biāo)模擬圖像的值，故15 m分辨率的模擬HyMap圖像質(zhì)量要好于其他模擬圖像。不同空間分辨率模擬圖像的可區(qū)分度較好，不同光譜分辨率模擬圖像的可區(qū)分度次之，不同信噪比模擬圖像的可區(qū)分性較差?？梢哉J(rèn)為，不同信噪比和不同光譜分辨率的模擬圖像隨著信噪比與光譜分辨率的變化，其圖像質(zhì)量變化不大;不同空間分辨率的模擬圖像隨著空間分辨率的變化，其圖像質(zhì)量變化較大。因此設(shè)計(jì)傳感器時(shí)，在光譜分辨率和信噪比一定的情況下，應(yīng)盡可能地提高傳感器的空間分辨率，這樣會(huì)使圖像質(zhì)量明顯提高。

2.1.2 基于直方圖異常的模擬圖像質(zhì)量分析

直方圖異常主要用于評(píng)價(jià)數(shù)據(jù)的動(dòng)態(tài)范圍，數(shù)據(jù)的動(dòng)態(tài)范圍越大，地物的層次越豐富、可分性越高。在對(duì)遙感圖像進(jìn)行預(yù)處理時(shí)，首先應(yīng)對(duì)圖像建立灰度直方圖以便了解圖像整體亮度分布情況，然后可通過對(duì)直方圖做均衡化和歸一化處理，對(duì)圖像質(zhì)量進(jìn)行調(diào)整?；叶戎狈綀D描述圖像灰度統(tǒng)計(jì)特性，故可以從直方圖上分析圖像噪聲情況，進(jìn)而評(píng)價(jià)模擬HyMap圖像的質(zhì)量。通常，圖像的直方圖呈近似正態(tài)分布。如果圖像受到灰度值很大、頻率很小的噪聲的影響，其直方圖將偏離正態(tài)分布，而被擠壓成很窄的形狀，表明有效信息受到抑制。為了方便做直方圖對(duì)比分析，對(duì)不同載荷指標(biāo)模擬圖像均取第一波段在450 nm處的直方圖曲線(圖2)進(jìn)行分析研究。

圖2 模擬圖像450 nm處直方圖曲線對(duì)比圖Fig.2 Comparison of histogram cures of simulated image at 450 nm

(1)不同空間分辨率模擬圖像直方圖異常分析。從圖2(a)中的對(duì)比分析可以看出，15 m分辨率模擬圖像的直方圖曲線的橫坐標(biāo)值即灰度級(jí)范圍比45、60和75 m分辨率模擬圖像的灰度級(jí)范圍大，且同一灰度級(jí)出現(xiàn)的頻率明顯偏高，而75 m分辨率模擬圖像的灰度級(jí)范圍最窄。說明隨著空間分辨率的提高，模擬圖像質(zhì)量也相應(yīng)提高。

(2)不同光譜分辨率模擬圖像直方圖異常分析。圖2(b)中隨著圖像亮度值范圍逐漸向高值方向移動(dòng)，5、15、20、25 nm 光譜分辨率圖像的直方圖曲線移動(dòng)步長(zhǎng)稍小且移動(dòng)步長(zhǎng)較均勻，而30 nm光譜分辨率圖像的直方圖曲線移動(dòng)步長(zhǎng)較其余光譜分辨率圖像直方圖曲線的移動(dòng)步長(zhǎng)大;縱軸所代表的不同灰度級(jí)值所出現(xiàn)的頻率的峰值基本穩(wěn)定且均衡，故主要差別在于灰度級(jí)閾值范圍的不同。

(3)不同信噪比模擬圖像直方圖異常分析。從圖2(c)中可看出，隨著信噪比增大，模擬圖像直方圖曲線的閾值范圍逐漸向右移動(dòng)，且移動(dòng)步長(zhǎng)不均衡，信噪比為450～500的模擬圖像直方圖曲線的移動(dòng)步長(zhǎng)較其他3種信噪比模擬圖像直方圖曲線的移動(dòng)步長(zhǎng)大;閾值范圍跨度較寬，信噪比為450～500的模擬圖像的均方差要比其余3種信噪比模擬圖像的均方差大;圖像質(zhì)量隨著信噪比增高而不斷提高，其中信噪比為450～500的模擬圖像的質(zhì)量最好，信噪比為250～300和350～400的模擬圖像幾乎重合，這3種信噪比模擬圖像質(zhì)量都要好于信噪比為50～100的模擬圖像質(zhì)量。

2.1.3 基于數(shù)據(jù)相關(guān)異常的模擬圖像質(zhì)量分析

HyMap圖像的相關(guān)性表現(xiàn)為空間相關(guān)性和譜間相關(guān)性，這里只討論譜間相關(guān)性。從光譜學(xué)角度來看，各波段之間的光譜反射率具有較強(qiáng)的相關(guān)性，且相關(guān)性的大小在很大程度上取決于光譜分辨率。同景圖像的不同波段圖像所涉及的地面目標(biāo)相同，并具有相同的空間拓?fù)浣Y(jié)構(gòu)?；谶@兩點(diǎn)因素，每個(gè)波段的光譜圖像在同一空間位置的像元具有相似性。對(duì)于不同波段間隔的HyMap圖像，各個(gè)波段包含的地物信息不一，各波段之間信息的重疊與分異程度亦不同，各類型圖像數(shù)據(jù)間的相關(guān)性也不同。

圖3 不同空間分辨率模擬圖像相關(guān)圖Fig.3 Correlation of simulated image with different spatial resolution

(1)不同空間分辨率模擬圖像相關(guān)異常。從不同空間分辨率模擬圖像相關(guān)圖(圖3)中可以看到，在整體上空間分辨率為15、45、60、70 m的模擬圖像相差無幾，但不能因此認(rèn)為這4種模擬圖像質(zhì)量基本相同。因?yàn)閺脑砩蟻砜矗@些不同空間分辨率的高光譜圖像是在相同條件下利用原始HyMap數(shù)據(jù)和模擬軟件模擬而成，故在這種情況下，空間分辨率高的模擬圖像包含更豐富的空間紋理信息，其質(zhì)量理應(yīng)比空間分辨率低的模擬圖像的質(zhì)量好。

(2)不同光譜分辨率模擬圖像相關(guān)異常。從不同光譜分辨率模擬圖像相關(guān)圖(圖4)的對(duì)比中可以看到，在各模擬圖像相關(guān)圖上都有一些值為零的黑色條帶，一般情況下這些零值的模擬圖像都是壞波段;另外，波段間的相關(guān)系數(shù)較大，相關(guān)性高，尤其是光譜分辨率為5 nm的模擬圖像的波段相關(guān)性比其余載荷指標(biāo)模擬圖像的波段相關(guān)性要高(且接近于1)，數(shù)據(jù)冗余現(xiàn)象極為明顯。光譜分辨率為15、20、25、30 nm的圖像的相關(guān)系數(shù)依次減小，所以從理論上來說，其質(zhì)量應(yīng)逐漸提高。但圖4中光譜分辨率為15 nm的模擬圖像的波段相關(guān)性卻比光譜分辨率為20、25 nm的模擬圖像的波段相關(guān)性低，造成這一異?，F(xiàn)象的可能原因是原始HyMap數(shù)據(jù)在模擬過程中所加噪聲不均一所致。

圖4 不同光譜分辨率模擬圖像相關(guān)圖Fig.4 Correlation of simulated image with different spectral resolution

(3)不同信噪比模擬圖像相關(guān)異常。從不同信噪比模擬圖像的波段相關(guān)圖(圖5)和相關(guān)系數(shù)矩陣來看，信噪比高的模擬圖像的波段相關(guān)性較低，且信噪比在250以上的模擬圖像的波段相關(guān)性要比信噪比在100以下的模擬圖像的波段相關(guān)性低。

圖5 不同信噪比的模擬圖像相關(guān)圖Fig.5 Correlation of simulated image with different SNR

綜上所述，空間分辨率高于30 m、信噪比在250以上和光譜分辨率低于15 nm的HyMap模擬圖像的整體質(zhì)量較好。

3 從地學(xué)應(yīng)用角度評(píng)價(jià)模擬數(shù)據(jù)質(zhì)量

3.1 高光譜遙感礦物識(shí)別概況

從巖礦信息提取的角度分析，國(guó)內(nèi)外發(fā)展的光譜識(shí)別方法從本質(zhì)上可以歸納為3大類型:①以重建光譜與標(biāo)準(zhǔn)光譜相似性度量為基礎(chǔ)的光譜匹配方法，如歐式距離及基于歐式距離的相似指數(shù)方法、光譜角填圖(Spectral Angle Mapper，SAM)技術(shù)、光譜相關(guān)填圖(Spectral Correlation Mapper，SCM)方法［8］、光 譜 信 息 散 度 (Spectral Information Divergence，SID)［9－10］和交叉相關(guān)光譜匹配(Cross Correlation Spectral Matching，CCSM)等;②基于診斷性光譜吸收譜帶特征參量的局部光譜識(shí)別方法，如光譜特征擬合(Spectral Feature Fitting，SFF)［11］、光譜吸收指數(shù)(Spectral Absorption Index，SAI)［4］、吸收譜帶定位分析(Analysis of Absorption Band Positioning，AABP)等;③以礦物學(xué)和礦物光譜知識(shí)為基礎(chǔ)的智能識(shí)別方法。

3.2 蝕變信息提取與礦物填圖

3.2.1 技術(shù)流程

蝕變信息提取與礦物填圖技術(shù)流程如圖6所示。

圖6 蝕變信息提取與礦物填圖技術(shù)流程圖Fig.6 Flow chart of the method of alteration information extraction and mineral mapping

3.2.2 信息提取、填圖結(jié)果和模擬參數(shù)關(guān)系分析

在礦化蝕變信息提取過程中，本文采用了SAM填圖和混合調(diào)制匹配濾波(Mixture－Tuned Matched Filtering，MTMF)技術(shù)相結(jié)合的方法，并在2010年5月進(jìn)行了野外驗(yàn)證和采樣，正確率在80%左右，效果較好。

(1)不同空間分辨率數(shù)據(jù)的提取結(jié)果。在礦物信息提取過程中，隨著空間分辨率的降低，端元獲取的難度變大，端元獲取的精度隨之降低。在進(jìn)行MTMF填圖時(shí)，隨著空間分辨率的降低，匹配度不斷變小(范圍是0～1.0)，不可行性 (Infeasability)值隨之升高，使獲取最佳匹配目標(biāo)的難度加大。從總體上來說，空間分辨率的降低，使識(shí)別端元種類明顯減少，空間分布細(xì)節(jié)減弱;而空間分辨率的提高，會(huì)使遙感影像更清晰、紋理特征更鮮明，混合像元效應(yīng)降低，蝕變信息提取的結(jié)果更好。從圖7中可以看出，空間分辨率為10 m的圖像比其余幾種空間分辨率較低的圖像所提取的蝕變礦物信息及礦物填圖的精度和準(zhǔn)確度都要高;研究認(rèn)為，要得到包含8種以上蝕變礦物類型的填圖結(jié)果，數(shù)據(jù)的空間分辨率需要高于30 m。

圖7 不同空間分辨率數(shù)據(jù)的蝕變信息提取結(jié)果Fig.7 Extraction of alteration information from different spatial resolution data

(2)不同信噪比數(shù)據(jù)提取結(jié)果。從不同信噪比數(shù)據(jù)的提取結(jié)果(圖8)可以看出，隨著信噪比的提高，識(shí)別端元種類略微增多，空間分布細(xì)節(jié)更加精細(xì)，所提取礦物種類有明顯差別，信噪比高的優(yōu)越性體現(xiàn)較好。對(duì)比分析認(rèn)為，要得到包含6種以上蝕變礦物類型的填圖結(jié)果，數(shù)據(jù)的信噪比需要高于150～200(VNIR/SWIR);要得到礦物類型更多、空間分布更精細(xì)的填圖結(jié)果，則數(shù)據(jù)的信噪比需要高于350～400;另外，從信噪比在350～400與450～500的數(shù)據(jù)所提取礦物的效果比從其余兩種信噪比的數(shù)據(jù)所得的結(jié)果好，所提取蝕變信息的精度較高。

圖8 不同信噪比數(shù)據(jù)的蝕變信息提取結(jié)果Fig.8 Extraction of alteration information from different SNR data

(3)不同光譜分辨率數(shù)據(jù)的提取結(jié)果。從不同光譜分辨率數(shù)據(jù)的提取結(jié)果(圖9)中可以看出，隨著光譜分辨率的提高，圖像所含信息顯著增強(qiáng);但同時(shí)噪聲也伴隨著圖像信息的增強(qiáng)而明顯增多，對(duì)識(shí)別蝕變礦物信息影響很嚴(yán)重，有時(shí)甚至掩蓋了某些蝕變礦物的特征光譜，導(dǎo)致蝕變礦物無法識(shí)別或加大了鑒別的難度。另外，隨著光譜分辨率的降低，識(shí)別端元種類減少得不明顯，但空間分布細(xì)節(jié)有所減弱，部分礦物出現(xiàn)混淆，應(yīng)用效果不僅受光譜分辨率的影響，而且與所選擇的數(shù)據(jù)處理方法、野外調(diào)查與驗(yàn)證工作等也有較大程度的關(guān)系。為滿足地質(zhì)調(diào)查礦物識(shí)別與填圖的應(yīng)用需求，對(duì)不同光譜分辨率的載荷性能指標(biāo)需求初步分析認(rèn)為，光譜分辨率在10～20 nm時(shí)最好。

圖9 不同光譜分辨率數(shù)據(jù)的蝕變信息提取結(jié)果Fig.9 Extraction of alteration information from different spectral resolution data

4 結(jié)論

(1)許多研究表明，很多地表物質(zhì)的吸收峰深度(banddepth)一半處為 20 ～40 nm［12］。結(jié)合本文研究，為滿足地質(zhì)調(diào)查工作中礦化蝕變信息提取與填圖的應(yīng)用需求，經(jīng)初步分析確定的星載高光譜傳感器的載荷性能指標(biāo)為:光譜分辨率10～20 nm，空間分辨率10～30 m，信噪比高于200。

(2)總的來說，光譜分辨率、空間分辨率和信噪比之間是相互制約的。要保證高光譜分辨率數(shù)據(jù)的信息，信噪比就要做出讓步(相對(duì)來說，信噪比級(jí)別可能會(huì)降低);反之，要提高信噪比，必然要加寬波段的帶寬，致使光譜分辨率降低?？臻g分辨率與信噪比的相互關(guān)系亦然——大像元肯定會(huì)有高的信噪比，且噪聲干擾相對(duì)較小。

(3)不同載荷指標(biāo)的模擬HyMap數(shù)據(jù)質(zhì)量評(píng)價(jià)還存在以下問題與不足:①數(shù)據(jù)源較單一且受地域環(huán)境限制;②從圖像預(yù)處理角度進(jìn)行數(shù)據(jù)質(zhì)量評(píng)價(jià)的指標(biāo)比較單一，且評(píng)價(jià)過程中有一定的主觀因素。

［1］馬德敏.高光譜圖像質(zhì)量評(píng)價(jià)［J］.紅外，2004(增刊):18－23.

［2］汪孔橋，Kangas J A.數(shù)字圖像的質(zhì)量評(píng)價(jià)［J］.測(cè)控技術(shù)，2000，19(5):14 －16.

［3］王長(zhǎng)耀，牛 錚，唐華俊，等，對(duì)地觀測(cè)技術(shù)與精準(zhǔn)農(nóng)業(yè)［M］.北京:科學(xué)出版社，2001.

［4］王晉年，鄭蘭芬，童慶禧.成像光譜圖像光譜吸收鑒別模型與礦物填圖研究［J］.環(huán)境遙感，1996，11(1):20 －31.

［5］Jensen J R.遙感數(shù)字影像處理導(dǎo)論［M］.陳曉玲，龔 威，李平湘，等譯.北京:機(jī)械工業(yè)出版社，2007.

［6］王 昱.數(shù)字遙感影像構(gòu)象質(zhì)量評(píng)價(jià)方法初探［J］.遙感信息，2000(4):32－33.

［7］熊興華.數(shù)字影像質(zhì)量評(píng)價(jià)方法評(píng)述［J］.測(cè)繪科學(xué)，2004，29(1):68－71.

［8］Carvalho O A J，Menese P R.Spectral Correlation Mapper(SCM):An Improvement on the Spectral Angle Mapper(SAM)［C］//Ninth JPL Airborne Earth Science Workshop.USA:JPL publication，2000，18:65 －74.

［9］Chang C I，Plaza A.A Fast Iterative Algorithm for Implementation of Pixel Purity Index［J］.IEEE Geoscience and Remote Sensing Letters，2006，3(1):63 －67

［10］Chang C I.An Information － theoretic Approach to Spectral variability，Similarity，and Discrimination for Hyperspectral Image Analysis［J］.IEEE Information Theory，2000，46(5):1927 －1932.

［11］Clark R N，King T V V，Klejwa M，et al.High Spectral Resolution Reflectance Spectroscopy of Minerals［J］.Journal of Geophysical Research，1990，95(10):12653 －12680.

［12］Hunt G R.Near－ infrared(1.3 ～2.4 μm)Spectral of Alteration Minerals － potential for Use in Remote Sensing［J］.IEEE Transactions on Geoscience and Remote Sensing，1979，42(3):501 －513.

Evaluation of the Quality of HyMap Data Simulated with Different Payload Indexes

HUO Hong－yuan，ZHOU Ping，NI Zhuo－ya
(School of Earth Sciences and Resurces，China University of Geosicences(Beijing)，Beijing 100083，China)

In this paper，the quality of HyMap data simulated with different payload indexes was evaluated in two aspects，namely，pre － processing and geoscience applications of the HyMap image data.In the aspect of simulated data pre － processing，three parameters were used to evaluate the quality of the simulated HyMap image data，i.e.，average variance abnormality，histogram abnormality and correlation abnormality.In the aspect of geoscience application of the simulated HyMap data，the simulated HyMap image data were used to extract the alteration information of mineralization and mineral mapping，and the quality of simulated HyMap data with different payload indexes or with the same payload indexes but different scales was analyzed and evaluated according to the number and kinds of the information extracted and the precision of mineral mapping.Some good results have been achieved，and several valuable advices and suggestions are put forward for the study of HyMap sensor.

Qulity evalution of HyMap data;Image pre－processing;Information extraction;Mineral mapping;Payload index

TP 79

A

1001－070X(2011)04－0052－06

2011－03－30;

2011－05－19

霍紅元(1985－)，男，碩士研究生，主要從事高光譜地學(xué)應(yīng)用方面的研究。

(責(zé)任編輯:劉心季)