顧及礦物粒徑的月表虹灣地區(qū)礦物填圖

董曉瑩， 林偉華， 劉福江， 張 琪， 常 遠(yuǎn)

(1.中國(guó)地質(zhì)大學(xué)(武漢)信息工程學(xué)院，武漢 430074； 2.吉林省交通規(guī)劃設(shè)計(jì)院，長(zhǎng)春 130021)

0 引言

月球表面礦物組成的測(cè)定是理解月球成因和地質(zhì)演化的基礎(chǔ)，這些成分信息可以用來(lái)研究月球地殼形成時(shí)的巖漿模型、月球的地殼結(jié)構(gòu)、玄武巖火山、火山口/盆地結(jié)構(gòu)及噴出物的流溢，以及理解月壤的混合機(jī)理[1]。可見(jiàn)光/近紅外反射光譜(0.4～2.5 μm)是研究月壤礦物組成的重要數(shù)據(jù)之一，大多數(shù)月球礦物在該波段范圍有比較顯著的吸收特征(中心位置： 750 nm和950 nm)。隨著月表高光譜數(shù)據(jù)探測(cè)技術(shù)的發(fā)展，對(duì)月表礦物的分析逐漸趨向于定量反演[2]，高光譜數(shù)據(jù)憑借其具有連續(xù)光譜、高光譜分辨率和高空間分辨率成為當(dāng)前研究的主要數(shù)據(jù)。印度探月衛(wèi)星Chandrayaan-1獲取的月球礦物制圖儀(moon mineralogy mapper，M3)高光譜數(shù)據(jù)可用于較精確地定量分析月球礦物的光譜特征和含量分布[2]。

由于地物混合的多樣性以及高光譜數(shù)據(jù)空間分辨率的限制，混合光譜現(xiàn)象在遙感影像中普遍存在。因此，混合光譜分解[3-4]成為利用高光譜數(shù)據(jù)進(jìn)行礦物定量反演研究的熱點(diǎn)和難點(diǎn)。Johnson等[5]利用半經(jīng)驗(yàn)法進(jìn)行非線性解混反演得到了月表橄欖石、輝石和鈦鐵礦等的含量； Mustard 等[6]利用Hapke模型進(jìn)行非線性解混反演月表礦物豐度； Li 等[1]根據(jù)Hapke模型采用多端元光譜分解的方法得到了月表礦物豐度圖； 閆柏琨等[7-8]基于光譜分解利用Clementine UV/VIS/NIR數(shù)據(jù)對(duì)全月進(jìn)行了礦物填圖； Combe等[9]基于混合光譜分析利用M3對(duì)礦物表面混合特性進(jìn)行了研究。

Hapke模型是一種常用的研究礦物反射特征的傳輸模型，可以很好地模擬礦物的混合光譜特征，具有嚴(yán)格的物理意義。在任何Hapke模型計(jì)算中礦物粒徑都是必須明確的參數(shù)之一，大量的實(shí)驗(yàn)室結(jié)果均證明了礦物粒徑對(duì)礦物的反射特性有顯著的影響。Li 等[10]基于Hapke模型利用月壤特征集(lunar soil characterization consortium，LSCC)數(shù)據(jù)定量研究了礦物粒徑對(duì)熔融玻璃、單斜輝石、斜方輝石、斜長(zhǎng)石、橄欖石、鈦鐵礦和火山玻璃等月表礦物光譜的影響以及亞微觀金屬鐵(submicroscopic metallic Fe，SMFe)的含量； Nelson等[11]比較了Hapke模型中礦物粒徑不同計(jì)算方法結(jié)果的異同； 張淵智等[12]利用Hapke模型和Mie理論模擬了不同粒徑下月表橄欖石的二向性特征，發(fā)現(xiàn)光譜反射率隨著粒徑的增加而降低； 吳昀昭等[13]利用Hapke模型考慮礦物粒徑與密度研究了太空風(fēng)化對(duì)其反射光譜的影響。

本文基于Hapke模型，利用線性分解模型在顧及礦物有效粒徑大小和密度對(duì)反射光譜的影響，用反射實(shí)驗(yàn)室(reflectence experiment laboratory，Relab)光譜庫(kù)數(shù)據(jù)建立模型，采用M3數(shù)據(jù)對(duì)月表虹灣地區(qū)典型礦物的含量及其分布填圖。

1 數(shù)據(jù)源

Relab光譜庫(kù)是美國(guó)Brown大學(xué)在美國(guó)航天航空局(NASA)的資助下建立的面向公眾的免費(fèi)數(shù)據(jù)庫(kù)。本文采用的5種礦物端元光譜分別為L(zhǎng)S-CMP-009(單斜輝石)、LS-CMP-012(斜方輝石)、LR-CMP-014(橄欖石)、LS-CMP-011(斜長(zhǎng)石)和PI-CMP-006(鈦鐵礦)，使用的光譜范圍是430～2 600 nm。針對(duì)月壤，該數(shù)據(jù)庫(kù)將礦物顆粒分為<10 μm，10～20 μm，20～45 μm和>45 μm 4個(gè)等級(jí)，本文使用的是10～20 μm的數(shù)據(jù)，并依據(jù)M3影像光譜范圍和光譜分辨率進(jìn)行重采樣。M3數(shù)據(jù)技術(shù)性能指標(biāo)詳見(jiàn)表1。不同的觀測(cè)條件下，光學(xué)周期不同，空間分辨率也有明顯的差別。在軌道高度為100 km時(shí)，包含3個(gè)光學(xué)周期： OP1A，OP1B和OP2A； 當(dāng)軌道高度為200 km時(shí)，包含有OP2B和OP2C光學(xué)周期。其中，OP1B和OP2C可覆蓋虹灣地區(qū)及月表采樣點(diǎn)。

表1 M3主要技術(shù)和性能指標(biāo)Tab.1 Main technical and performance indicates of M3

2 礦物反演模型建立

2.1 礦物單次散射反照率

月表不是理想的朗伯面，平行太陽(yáng)光經(jīng)月表物質(zhì)反射后的光強(qiáng)分布不是各向同性的，受到太陽(yáng)光入射角、反射光出射角以及兩者的夾角(即相位角)的影響，具有明顯的二向性反射特性[14]，其波譜特征與目標(biāo)物體的組成成分、物體內(nèi)部結(jié)構(gòu)和物體表面結(jié)構(gòu)關(guān)系密切[15]。月表二向性反射率是進(jìn)行月表組分識(shí)別和礦物含量反演的重要依據(jù)[16]，利用月表二向性特征、根據(jù)實(shí)測(cè)反射率可以得到月表礦物含量。然而，研究發(fā)現(xiàn)礦物光譜的二向性一般不呈現(xiàn)線性特征，不能直接進(jìn)行線性混合，但是礦物的單次散射反照率(single-scattering albedo, SSA)是各向同性的，可直接進(jìn)行線性混合[17]。Hapke模型可以實(shí)現(xiàn)二向性反射率與同向性的單次散射反照率之間的轉(zhuǎn)換，具有嚴(yán)格的物理意義。假設(shè)礦物為朗伯反射體，由Hapke模型[18]可知，礦物光譜的二向性反射率r(μ0,μ,g)的公式為

(1)

(2)

P(g)=1+bcos(g)+c[1.5cos2(g)-0.5] ，

(3)

式中：μ0為入射角余弦值；μ為出射角余弦值；g為相位角；ω為礦物SSA，即僅考慮一次散射時(shí)，散射能量與顆粒散射和吸收的能量之和的比值；B(g)為后向散射函數(shù)，定義了當(dāng)相位角減小時(shí)顆粒粗糙表面亮度值的增加量，當(dāng)相位角大于15°時(shí)，可認(rèn)為所有的散射都是同向性的，B(g)可忽略；H(μ)為同向性散射函數(shù)；P(g)為平均單粒子相函數(shù)，即顆粒散射能量隨著g變化的規(guī)律；b和c為常數(shù)。

Relab光譜庫(kù)中數(shù)據(jù)測(cè)量時(shí)，入射角i=30°，出射角e=0°，相位角g=30°，故B(g)=0，P(g)=1?？珊?jiǎn)化Hapke模型如式(4)，從而得到礦物SSA，即

(4)

2.2 礦物混合反射率

基于Hapke模型得到的礦物SSA，將端元礦物的SSA進(jìn)行隨機(jī)線性混合[19]，混合時(shí)根據(jù)已有的研究成果[1]，對(duì)5種成分礦物的含量進(jìn)行約束： 斜長(zhǎng)石( 0%～100%)、單斜輝石(0%～70%)、斜方輝石(0%～50%)、橄欖石(0%～50%)和鈦鐵礦(0%～20%)，并考慮顆粒的平均有效粒徑和密度，然后再利用Hapke模型計(jì)算得到混合反射率數(shù)據(jù)。

顆粒的平均有效粒徑和密度ρ等不同，礦物的反射率也會(huì)發(fā)生變化，因此這里計(jì)算混合反射率時(shí)，將與ρ考慮進(jìn)礦物混合公式[16]，即

(5)

式中：ω′為混合礦物SSA；ωi，Mi，ρi和di分別為混合物中第i種組分的SSA、質(zhì)量分?jǐn)?shù)、密度和平均有效粒徑。其中，顆粒的采用area-weighted mean size來(lái)表示，即

≈DLln(DU/DL)，

(6)

式中：DU為最大粒徑；DL為最小有效粒徑。文中用到的5種礦物信息詳見(jiàn)表2。

表2 Relab光譜庫(kù)端元礦物光學(xué)常數(shù)信息Tab.2 Optical constant information of endmembermineralsin Relab spectral library (nm)

得到混合反射率后需要對(duì)其進(jìn)行去連續(xù)統(tǒng)處理。該方法是一種用于分離光譜吸收特征與背景的光譜分析方法，可有效消除地形和照度等對(duì)光譜特征造成的影響，從而去除端元光譜與影像光譜之間由于照度不同所導(dǎo)致的光譜變異。

2.3 全約束線性分解

遙感探測(cè)器記錄的瞬時(shí)視場(chǎng)角對(duì)應(yīng)的地面范圍內(nèi)目標(biāo)的輻射能量總和被稱為遙感像元。通常一個(gè)遙感像元內(nèi)包含多種不同類型的地物，被稱為混合像元。目前，混合像元分解[3-4]是解決混合像元最為有效的方法，即將記錄混合像元的影像光譜分解為組成光譜(端元)與其相應(yīng)豐度。針對(duì)混合像元分解的模型，目前主要分為線性光譜混合模型與非線性光譜混合模型2類。其中線性光譜混合模型因其簡(jiǎn)單、效率高、物理意義明確而得到更廣泛的應(yīng)用。本文即采用線性光譜混合模型。

假設(shè)像元內(nèi)相同的地物具有相同的光譜特征并且忽略多次散射過(guò)程，混合像元反射率與端元反射率以及端元豐度之間的關(guān)系可以表示為[4]

(7)

2.4 反演模型建立

充分考慮2個(gè)約束條件，利用線性分解得到的各端元礦物的含量記為“礦物分解含量”； SSA隨機(jī)混合時(shí)設(shè)定的各端元礦物的含量記為“真實(shí)含量”，根據(jù)統(tǒng)計(jì)關(guān)系建立5種端元礦物(單斜輝石、斜方輝石、斜長(zhǎng)石、橄欖石和鈦鐵礦)的豐度反演模型[21-22]，如圖1所示。

(a) 單斜輝石 (b) 斜方輝石 (c) 斜長(zhǎng)石

(d)橄欖石 (e) 鈦鐵礦

圖1不同端元礦物分解含量與真實(shí)含量的統(tǒng)計(jì)關(guān)系

Fig.1Statisticalrelationsbetweentheunmixingabundanceandtherealabundanceofdifferentendmemberminerals

由圖1可以看到，單斜輝石、斜方輝石、斜長(zhǎng)石、橄欖石和鈦鐵礦這5種礦物的“分解含量”與“真實(shí)含量”之間的線性相關(guān)系數(shù)分別是0.98，0.98，0.83，0.91，0.50。單斜輝石、斜方輝石、斜長(zhǎng)石和橄欖石的線性相關(guān)性均比較高，鈦鐵礦相對(duì)來(lái)說(shuō)較低，這與其在月球中本身含量較低以及不透明礦物的特性有關(guān)。與之前的研究[23]相比，將礦物粒徑對(duì)其反射特性的影響考慮進(jìn)來(lái)建立的模型，相關(guān)性均有比較明顯的提高。尤其是單斜輝石和斜方輝石的模型相關(guān)性基本接近于1，說(shuō)明本實(shí)驗(yàn)?zāi)Ｐ涂梢暂^好地反演礦物含量，可以用來(lái)進(jìn)行月表礦物豐度填圖。

2.5 模型驗(yàn)證

為了驗(yàn)證反演模型的精度，將該方法應(yīng)用到Apollo樣品數(shù)據(jù)上，得到礦物分解含量，并與樣品的實(shí)測(cè)礦物含量進(jìn)行對(duì)比。其中，輝石是指單斜輝石與斜方輝石含量之和，結(jié)果如表3和圖2所示。

表3 Apollo 采樣點(diǎn)礦物反演結(jié)果與實(shí)測(cè)結(jié)果對(duì)比Tab.3 Comparison between the inversed abundance and the measured abundance of differentendmember minerals in Apollo (%)

(a) 輝石 (b) 斜長(zhǎng)石

(c) 橄欖石 (d) 鈦鐵礦

圖2Apollo采樣點(diǎn)礦物反演結(jié)果與實(shí)測(cè)結(jié)果統(tǒng)計(jì)關(guān)系

Fig.2Statisticalrelationsbetweentheinversedabundanceandthemeasuredabundanceofdifferent

endmembermineralsinApollo

從表3和圖2可以看出，相比實(shí)測(cè)含量，斜長(zhǎng)石、橄欖石和鈦鐵礦的含量均被高估，而反演輝石的含量略低于實(shí)驗(yàn)室測(cè)量得到的真實(shí)值。造成這一現(xiàn)象的原因可能有以下幾點(diǎn)：

1)端元選擇不夠或端元特征不夠明顯會(huì)造成礦物反演結(jié)果有偏差[3,24]，本文選取了月表豐度相對(duì)較高并且有明顯反射特征的5種端元光譜來(lái)進(jìn)行實(shí)驗(yàn)。從Apollo帶回來(lái)的樣品中可以發(fā)現(xiàn)，熔融玻璃與火山玻璃也是月表的主要礦物，但是由于其是風(fēng)化產(chǎn)物，光譜反射特征已基本消失，因此沒(méi)有對(duì)其進(jìn)行反演。

2)斜長(zhǎng)石的含量被高估，且反演得到的斜長(zhǎng)石含量與實(shí)測(cè)中斜長(zhǎng)石與熔融玻璃含量之和比較接近，可能與斜長(zhǎng)石的反射特性有關(guān)。斜長(zhǎng)石整體反射率較高，且經(jīng)過(guò)月表風(fēng)化作用，吸收特征被弱化，而熔融玻璃作為風(fēng)化產(chǎn)物，其基本沒(méi)有反射特征，反射光譜曲線變得平滑，兩者比較接近，因此造成了斜長(zhǎng)石含量被高估。

3)鈦鐵礦的反演結(jié)果被高估可能與鈦鐵礦自身屬性有關(guān)，鈦鐵礦作為不透明礦物，在光譜分解時(shí)其含量通常會(huì)被高估[25]。

4)輝石的含量低于實(shí)測(cè)含量值，這可能與斜長(zhǎng)石等含量被高估有關(guān)，輝石的含量與長(zhǎng)石含量存在高度的負(fù)相關(guān)[24]，而輝石與橄欖石在1 000 nm附近有一個(gè)共同的吸收特征。

總體來(lái)講，輝石、斜長(zhǎng)石、橄欖石和鈦鐵礦這4種礦物的“分解含量”與“實(shí)測(cè)含量”均有較好的相關(guān)性，這說(shuō)明本實(shí)驗(yàn)?zāi)Ｐ湍茌^好地描述這4種礦物的分布，可以用于礦物填圖。

2.6 虹灣地區(qū)礦物填圖

本文采用M3L2級(jí)數(shù)據(jù)，已經(jīng)過(guò)光度校正和幾何校正等預(yù)處理。為了進(jìn)一步削弱噪聲的影響，利用Savitzky-Golay濾波法對(duì)圖像進(jìn)行平滑處理。Savitzky-Golay 濾波是廣泛用于數(shù)據(jù)流平滑除噪的一種方法，可在確保光譜曲線的形狀和寬度不變并保持光譜特征的情況下去除光譜噪聲，達(dá)到平滑去噪的效果[26]。處理流程如圖3所示。

圖3M3數(shù)據(jù)處理流程

Fig.3ProcessflowchartofM3data

利用上述方法得到5種端元礦物的分解含量后，采用該反演模型，得到礦物在虹灣地區(qū)的填圖結(jié)果，如圖4所示。得到5種礦物在虹灣地區(qū)含量均值分別為： 7.3 wt%，0.7 wt%，86.7 wt%，13.2 wt%和9.8 wt%。

(a) 單斜輝石 (b) 斜方輝石 (c) 斜長(zhǎng)石

(d) 橄欖石 (e) 鈦鐵礦

圖4虹灣地區(qū)礦物分布

Fig.4MineralabundancedistributioninSinusIridum

由圖4可以看出，①單斜輝石主要分布在虹灣內(nèi)部，一些撞擊坑中有較高分布，高地幾乎沒(méi)有分布； ②斜方輝石在與虹灣相接的西北部高地有少量分布，而虹灣內(nèi)部含量幾乎為0； ③斜長(zhǎng)石廣泛分布于虹灣及其附近區(qū)域且豐度很高，尤其在高地區(qū)域，虹灣內(nèi)部豐度相對(duì)較低且沿雨海延伸方向逐漸減少，到與虹灣相連的雨海中僅有少量分布； ④橄欖石主要分布在與虹灣相連的雨海中，在虹灣內(nèi)部?jī)H有零星分布，此外在高地的一些區(qū)域也發(fā)現(xiàn)有少量橄欖石分布； ⑤鈦鐵礦主要分布在與虹灣相接的雨海及虹灣內(nèi)部，在高地僅有少量分布。

將本實(shí)驗(yàn)結(jié)果與之前未考慮礦物粒徑影響得到的結(jié)果[23,26]相比，輝石的分布基本不變，只是本文得到的單斜輝石含量較低，且高地分布基本為0； 與虹灣相連的雨海內(nèi)部有零星分布這一研究結(jié)果與李嬋等[26]根據(jù)礦物吸收特征利用IIM影像在該地區(qū)得到的輝石分布相符； 本文反演得到了少量的斜方輝石，且主要分布在高地，在虹灣內(nèi)部及雨?；緵](méi)有分布，這一結(jié)果與Lucey和Staid等的研究相符[24,27]； 斜長(zhǎng)石的分布基本不變； 本文得到的橄欖石主要分布在與虹灣相連的雨海中，在虹灣內(nèi)部?jī)H有零星分布，且在高地的一些區(qū)域也發(fā)現(xiàn)有少量橄欖石分布，這一結(jié)果與Lucey之前的研究一致[24]； 本文得到的鈦鐵礦在虹灣內(nèi)部有較廣泛分布，而在高地僅散落于一些區(qū)域上，與李嬋等[26]利用根據(jù)礦物吸收特征在該地區(qū)TiO2的填圖結(jié)果基本一致，而鈦鐵礦的分布主要取決于TiO2的分布[24]。由此可知，本文將礦物粒徑考慮進(jìn)反演模型提高了模型反演的精度，結(jié)果更具可信性。

3 結(jié)論與展望

礦物粒徑的差異會(huì)引起反射光譜的變化，對(duì)輝石及鈦鐵礦的影響尤為明顯[12]。本文顧及粒徑等對(duì)礦物反射光譜造成的影響，利用全約束線性分解的方法，得到了虹灣地區(qū)5種主要礦物豐度分布。實(shí)驗(yàn)結(jié)果表明，將粒徑等影響考慮進(jìn)模型提高了模型的反演精度。反演得到的5種礦物在虹灣地區(qū)的分布為： 單斜輝石是虹灣地區(qū)輝石的主要存在形式，主要分布虹灣內(nèi)部，一些小撞擊坑中含量較高，在高地基本沒(méi)有分布； 斜方輝石在月表含量較少，主要分布在高地，在虹灣內(nèi)部和雨?；緵](méi)有分布； 斜長(zhǎng)石主要分布在高地，在虹灣內(nèi)部也有分布，相對(duì)較少，且與輝石的分布呈現(xiàn)負(fù)相關(guān)； 橄欖石主要分布與虹灣相連的雨海中，在虹灣內(nèi)部及高地有少量分布； 鈦鐵礦主要分布在雨海及虹灣內(nèi)部，在高地僅有零星分布，含量較低。

但是，本文僅選取了單斜輝石、斜方輝石、斜長(zhǎng)石、橄欖石和鈦鐵礦5種典型礦物，且每種礦物僅選取了一個(gè)端元光譜，這可能導(dǎo)致結(jié)果存在誤差； 且月球由于沒(méi)有大氣層，太空風(fēng)化現(xiàn)象十分嚴(yán)重，而并未將太空風(fēng)化效應(yīng)[28]考慮進(jìn)模型，也會(huì)使結(jié)果產(chǎn)生誤差。此外，王景然[29]研究發(fā)現(xiàn)，地形起伏雖不會(huì)引起礦物分布的變化，但是會(huì)影響其豐度。因此，為了進(jìn)一步提高反演結(jié)果的準(zhǔn)確性，下一步將進(jìn)一步優(yōu)化程序并考慮將礦物的地形、粒徑和風(fēng)化程度對(duì)礦物光譜的影響加入模型建立光譜庫(kù)，同一種礦物將對(duì)應(yīng)多個(gè)端元光譜進(jìn)行礦物反演。

志謝： 感謝NASA以及Relab實(shí)驗(yàn)室提供本文所用實(shí)驗(yàn)數(shù)據(jù)！

參考文獻(xiàn)(References)：

[1] Li L,Lucey P G.Use of multiple endmember spectral mixture analysis and radiative transfer model to derive lunar mineral abundance maps[C]//Proceedings of 40th Lunar and Planetary Science Conference.Woodlands:Lunar and Planetary Science,2009.

[2] Pieters C M,Boardman J,Buratti B,et al.The moon mineralogy mapper(M3)on Chandrayaan-1[J].Current Science,2009,96(4):500-505.

[3] 李二森.高光譜遙感圖像混合像元分解的理論與算法研究[D].鄭州:解放軍信息工程大學(xué),2011.

Li E S.Research on Theory and Algorithms of Mixed Pixel Decomposition for Hyperspectral Remote Sensing Image[D].Zhengzhou:PLA Information Engineering University,2011.

[4] Keshava N,Mustard J F.Spectral unmixing[J].IEEE Signal Processing Magazine,2002,19(1):44-57.

[5] Johnson P E,Smith M O,Taylor-George S,et al.A semiempirical method for analysis of the reflectance spectra of binary mineral mixtures[J].Journal of Geophysical Research,1983,88(B4):3557-3561.

[6] Mustard J F,Pieters C M.Photometric phase functions of common geologic minerals and applications to quantitative analysis of mineral mixture reflectance spectra[J].Journal of Geophysical Research,1989,94(B10):13619-13634.

[7] 閆柏琨,甘甫平,王潤(rùn)生,等.基于光譜分解的Clementine UV/VIS/NIR數(shù)據(jù)月表礦物填圖[J].國(guó)土資源遙感,2009,21(4):19-24.doi:10.6046/gtzyyg.2009.04.04.

Yan B K,Gan F P,Wang R S,et al.Mineral mapping of the lunar surface using Clementine UV/VIS/NIR data based on unmixing of spectral[J].Remote Sensing for Land and Resources,2009,21(4):19-24.doi:10.6046/gtzyyg.2009.04.04.

[8] Yan B K,Wang R S,Gan F P,et al.Minerals mapping of the lunar surface with Clementine UVVIS/NIR data based on spectra unmixing method and Hapke model[J].Icarus,2010,208(1):11-19.

[9] Combe J P,McCord T B,Hayne P O,et al.Mapping of lunar volatiles with moon mineralogy mapper spectra:A challenge due to thermal emission[C]//EPSC-DPS Joint Meeting 2011.Nantes,France:COPERNICUS,2011.

[10] Li S,Li L.Radiative transfer modeling for quantifying lunar surface minerals,particle size,and submicroscopic metallic Fe[J].Journal of Geophysical Research,2011,116(E9):E09001.

[11] Nelson M L,Clark R N.Comparison of different measures of the grain size in Hapke modeling[C]//Bulletin of the American Astronomical Society,1990,22:1033.

[12] 張淵智,安 璐,黃朝君.基于太空風(fēng)化的尖晶石二向性反射特性研究[J].深空探測(cè)學(xué)報(bào),2014,1(3):210-213.

Zhang Y Z,An L,Huang Z J.The study of bidirectional reflectance feature of the spinel based on the space weathering[J].Journal of Deep Space Exploration,2014,1(3):210-213.

[13] 吳昀昭,鄭永春,鄒永廖,等.基于非線性混合模型研究太空風(fēng)化對(duì)月壤光譜的影響[J].空間科學(xué)學(xué)報(bào),2010,30(2):154-159.

Wu Y Z,Zheng Y C,Zou Y L,et al.Research of the optical effects of space weathering on lunar regolith based on the nonlinear mixing model[J].Chinese Journal of Space Science,2010,30(2):154-159.

[14] 歐陽(yáng)自遠(yuǎn).月球科學(xué)概論[M].北京:中國(guó)宇航出版社,2005:38.

Ouyang Z Y.Introduction to Lunar Science[M].Beijing: China Astronautic Publishing House,2005:38.

[15] 路 鵬,陳圣波,崔騰飛,等.月球表面礦物二向性反射特性實(shí)驗(yàn)研究[J].巖石學(xué)報(bào),2016,32(1):107-112.

Lu P,Chen S B,Cui T F,et al.Experimental study on bidirectional reflectance characteristics of minerals on lunar surface[J].Acta Petrologica Sinica,2016,32(1):107-112.

[16] 崔騰飛.月表礦物二向性反射模型研究[D].長(zhǎng)春:吉林大學(xué),2012.

Cui T F.Study on Bidirectional Reflectance Model of Minerals on Lunar Surface[D].Changchun:Jilin University,2012.

[17] Hapke B.Theory of Reflectance and Emittance Spectroscopy[M].New York:Cambridge University Press,1993.

[18] Hapke B.Space weathering from Mercury to the asteroid belt[J].Journal of Geophysical Research,2001,106(E5):10039-10073.

[19] 閆柏琨,李建忠,甘甫平,等.一種月壤主要礦物組分含量反演的光譜解混方法[J].光譜學(xué)與光譜分析,2012,32(12):3335-3340.

Yan B K,Li J Z,Gan F P,et al.A spectral unmixing method of estimating main minerals abundance of lunar soils[J].Spectroscopy and Spectral Analysis,2012,32(12):3335-3340.

[20] Heinz D C,Chang C I.Fully constrained least squares linear spectral mixture analysis method for material quantification in hyperspectral imagery[J].IEEE Transactions on Geoscience and Remote Sensing,2001,39(3):529-545.

[21] 燕守勛,張 兵,趙永超,等.高光譜遙感巖礦識(shí)別填圖的技術(shù)流程與主要技術(shù)方法綜述[J].遙感技術(shù)與應(yīng)用,2004,19(1):52-63.

Yan S X,Zhang B,Zhao Y C,et al.Summarizing the technical flow and main approaches for discrimination and mapping of rocks and minerals using hyperspectral remote sensing[J].Remote Sensing Technology and Application,2004,19(1):52-63.

[22] 王振超.月壤光譜特性分析與月表礦物信息定量反演[D].北京:中國(guó)地質(zhì)大學(xué)(北京),2011.

Wang Z C.Analysis of Spectral Characteristics of Lunar Soil and Quantitative Inversion of Minerals Information[D].Beijing: China University of Geoscience(Beijing),2011.

[23] 張 琪,劉福江,李 嬋,等.全約束線性分解的月表虹灣地區(qū)礦物含量反演[J].國(guó)土資源遙感,2016,28(1):7-14.doi:10.6046/gtzyyg.2016.01.02.

Zhang Q,Liu F J,Li C,et al.Fully constrained linear-unmixing for inversion of lunar mineral abundance in Sinus Iridum[J].Remote Sensing for Land and Resources,2016,28(1):7-14.doi:10.6046/gtzyyg.2016.01.02.

[24] Lucey P G.Mineral maps of the Moon[J].Geophysical Research Letters,2004,31(8):L08701.

[25] 薛 彬,楊建峰,趙葆常.月球表面主要礦物反射光譜特性研究[J].地球物理學(xué)進(jìn)展,2004,19(3):717-720.

Xue B,Yang J F,Zhao B C.The study of spectral feature of major minerals on the lunar surface[J].Progress in Geophysics,2004,19(3):717-720.

[26] 李 嬋,劉福江,鄭小坡,等.月表虹灣地區(qū)輝石及橄欖石含量反演[J].中國(guó)科學(xué)(物理學(xué)力學(xué)天文學(xué)),2013,43(11):1387-1394.

Li C,Liu F J,Zheng X P,et al.Lunar pyroxene and olivine abundance analysis of Sinus Iridum[J].Science China Physics,Mechanics and Astronomy,2013,43(11):1387-1394.

[27] Staid M I,Pieters C M.Mineralogy of the last lunar basalts:Results from Clementine[J].Journal of Geophysical Research,2001,106(E11):27887-27900.

[28] 付曉輝,鄒永廖,鄭永春,等.月球表面太空風(fēng)化作用及其效應(yīng)[J].空間科學(xué)學(xué)報(bào),2011,31(6):705-715.

Fu X H,Zou Y L,Zheng Y C,et al.Space weathering processes and effects on the Moon[J].Chinese Journal of Space Science,2011,31(6):705-715.

[29] 王景然.基于地形校正月表礦物含量反演研究[D].長(zhǎng)春:吉林大學(xué),2012.

Wang J R.Topographic Correction Based Retrieval of Lunar Mineral Abundance[D].Changchun:Jilin University,2012.