彭金英，蔣金豹，崔希民，張 建，王思佳，王明國，聶 鋮

(1.中國礦業(yè)大學(北京) 地球與測繪工程學院,北京 100083；2.北京師范大學 地理科學學部，北京 100875；3.寧夏農業(yè)技術推廣總站，寧夏 銀川 750001)

1 研究區(qū)及實驗數(shù)據(jù)

1.1 野外采樣

研究區(qū)位于寧夏回族自治區(qū)平羅縣，野外采樣在2018年4月8～10日進行，每個采樣區(qū)由5個采樣點組成，采樣點間隔約10 m，呈矩形分布，每個樣點采集0～5 cm表層土壤，并記錄采樣區(qū)中心點GPS位置信息，共采集59份土壤樣本。

1.2 實驗數(shù)據(jù)獲取

將土壤樣本帶回實驗室烘干、研磨、去除砂礫、植物根莖殘體及磚瓦片等侵入體，并過2 mm孔篩，然后將每份土樣分成兩份，分別用于光譜數(shù)據(jù)采集和土壤理化性質測試分析。

土壤光譜測定在暗室內進行，采用的光譜儀為SVC HR-1024i，其光譜范圍是350～2 500 nm，探頭視場角為8°，光源采用鹵素燈；土樣放入黑色橡膠器皿，用尺子將土面刮平，將其放置于黑色底墊上，光源照射方向與垂直方向夾角15°，光源距離土壤表面30 cm，探頭位于土壤表面垂直上方15 cm處。每個樣本在測量前均進行了白板校正，每個土樣測量4個方向(互為90°)，每個方向上測兩次光譜曲線數(shù)據(jù)，共8次，取其算術平均值作為最后的光譜反射率，如圖1所示。反射率光譜數(shù)值在5～70之間，高反射率與低反射率曲線間差異較大，但各曲線的變化趨勢大致相同。

圖1 原始光譜反射率圖

1.3 光譜數(shù)據(jù)預處理

光譜微分和連續(xù)統(tǒng)去除是光譜數(shù)據(jù)預處理的常用方法，一階導數(shù)可以扣除斜線背景影響，提高光譜分辨率；連續(xù)統(tǒng)去除將反射率歸一化為0～1，光譜的吸收特征也歸一化到一致的光譜背景上，有利于與其他光譜曲線進行特征數(shù)值的比較。

1.3.1 光譜一階微分

光譜低階微分處理對噪聲的敏感性較低，采用導數(shù)光譜技術可以消除光譜數(shù)據(jù)之間的系統(tǒng)誤差，減弱對目標光譜的影響，以便提取可識別地物的光譜吸收峰參數(shù)[13]；光譜一階微分可以消除背景噪聲，分辨重疊光譜[14]。

1.3.2 連續(xù)統(tǒng)去除

連續(xù)統(tǒng)去除(Continuum Removed, CR)可以凸顯地物的吸收特征，原始光譜經(jīng)連續(xù)統(tǒng)去除處理后，可明顯提高反演模型的預測性能[15]。

1.4 相關性分析與特征波段選擇

利用pearson相關性分析，計算各離子間的相關系數(shù)(絕對值)，計算公式見式(1)：

(1)

圖2 一階微分(左)、連續(xù)統(tǒng)去除光譜(右)與離子的相關系數(shù)

根據(jù)相關系數(shù)結果選擇相關性最強的波段為敏感波段。將光譜在1 501 nm、1 916 nm、2 097 nm和2 366 nm波段作為一階微分的敏感波段；2 035 nm、2 164 nm、2 258 nm和2 401 nm波段為連續(xù)統(tǒng)去除的敏感波段。在0.01顯著性水平下，相關系數(shù)絕對值均在0.7以上(**P0.01=0.330)。

2 模型構建方法

2.1 偏最小二乘回歸

偏最小二乘回歸(Partial Least Squares Regression，PLSR)建立在一般最小二乘的基礎上，可以在自變量存在嚴重多重相關性、樣本點個數(shù)少于變量個數(shù)的條件下進行回歸建模，可以較好地解決以往用普通多元回歸無法解決的問題，更易于辨識系統(tǒng)信息與噪聲[16]。

2.2 粒子群優(yōu)化徑向基神經(jīng)網(wǎng)絡

徑向基(Radial Basis Function，RBF)網(wǎng)絡結構簡單、訓練簡潔且學習收斂速度快，能夠逼近任意非線性函數(shù)，而粒子群(Particle Swarm Optimization，PSO)是一種群體智能的優(yōu)化算法，可實現(xiàn)個體在可解空間中的尋優(yōu)[17]。針對RBF神經(jīng)網(wǎng)絡模型結構和參數(shù)難以設置的缺點，利用PSO算法尋優(yōu)能力和收斂能力較強的優(yōu)點進行優(yōu)化，得到PSO優(yōu)化的RBF神經(jīng)網(wǎng)絡算法(PSO-RBF神經(jīng)網(wǎng)絡算法)，提高了網(wǎng)絡的泛化能力和學習能力[18]。

2.3 模型精度驗證

在建模過程中，為保證精度，從59個樣本中選取41個樣本用來建立模型，其余18個樣本用來對模型進行檢驗。模型精度采用決定系數(shù)R2、均方根誤差RMSE和預測方差比RPD這3個指標來評價。決定系數(shù)R2越大，均方根誤差RMSE越小，建模精度越高；預測方差比RPD≥2，表明模型預測能力極好，1.4

3 模型對比分析

表1 模型精度對比分析

4 結 論

(3)本次野外實驗采集樣本數(shù)量較少，定量反演模型的普遍適用性還需要進一步驗證。在后續(xù)的工作中需要采集更多具有代表性的樣點數(shù)據(jù)，以便完善模型，使其具有普適性與穩(wěn)健性，并進一步把模型應用到我國高分五號遙感影像上，從而實現(xiàn)大范圍土壤鹽堿化遙感監(jiān)測與評估。