雷宇斌,劉 寧,郭云開,劉 磊,李丹娜

(1.湖南省第二測繪院，湖南 長沙 410000;2.中國能源建設(shè)集團(tuán)湖南省電力設(shè)計院有限公司，湖南 長沙 410007;3.長沙理工大學(xué)，湖南 長沙 410076)

為了充分挖掘光譜信息，需要進(jìn)行光譜變換，提高光譜對于重金屬的響應(yīng)能力。目前常用的土壤光譜變換方法主要有一階微分、二階微分、倒數(shù)對數(shù)、包絡(luò)線去除、多元散射校正、標(biāo)準(zhǔn)正態(tài)化、小波分析與傅里葉變換等[2-7]。

由于目前土壤高光譜反演模型往往忽略了重金屬的變異程度，而本文研究區(qū)的重金屬Cd元素的變異程度高達(dá)96%，傳統(tǒng)未嵌入地理位置參數(shù)的模型難以對高度變異的Cd元素進(jìn)行準(zhǔn)確建模與預(yù)測，因此本文在光譜組合變換研究的基礎(chǔ)上選取了能有效揭示土壤變異性空間關(guān)系的GWR模型來進(jìn)行重金屬Cd元素的反演研究。

1 數(shù)據(jù)的獲取與處理

1.1 研究區(qū)概況

選取湖南瀏陽市某硫鐵礦區(qū)周邊耕地為研究對象，該地區(qū)占地總面積234.5 km2，區(qū)域地貌以丘陵和低山坡為主，境內(nèi)氣候隸屬于亞熱帶季風(fēng)性濕潤氣候，熱量充足，降水豐沛，光照較足，多年平均氣溫17 ℃。 該區(qū)域作為老工業(yè)基地，礦產(chǎn)資源豐富，己探明有磷礦石、菊花石、海泡石等礦物20多種。當(dāng)?shù)氐V產(chǎn)資源開發(fā)的同時加劇了耕地土壤生態(tài)環(huán)境的破壞，為了人民生活、生產(chǎn)幸福安全，需要研究礦區(qū)周圍耕地重金屬污染程度，因此開展該區(qū)域耕地重金屬污染評價與反演工作。

1.2 土壤與光譜數(shù)據(jù)采集

本次研究實驗采用AvaFiled-3高精度光譜地物波譜儀采集土壤光譜數(shù)據(jù)，采集時間為2017年5月，測量前先進(jìn)行白板校正，每個土壤樣點(diǎn)采集數(shù)條光譜曲線，取平均值后作為該樣點(diǎn)的實際光譜反射數(shù)據(jù)。土壤取樣根據(jù)土壤質(zhì)地、作物覆蓋類型在試驗區(qū)使用T型取土器采集5個土樣點(diǎn)，采樣點(diǎn)位根據(jù)實地情況呈蛇形分布，所采集的土樣使用密封裝置帶回實驗室，并對每個樣點(diǎn)進(jìn)行相應(yīng)行編號與簡單的信息記錄。實驗時采用GPS定位儀器獲取每個采樣點(diǎn)的地理坐標(biāo)，實驗最終測得40個樣方點(diǎn)，樣點(diǎn)圖如圖1所示。

圖1 采樣點(diǎn)位置圖

1.3 土壤重金屬理化性質(zhì)分析

重金屬含量的測定依據(jù)《土壤環(huán)境監(jiān)測技術(shù)規(guī)范》的相關(guān)標(biāo)準(zhǔn)進(jìn)行。研究測定了Pb、Cu、Cd 3種重金屬，結(jié)果如表1所示。表中三級標(biāo)準(zhǔn)均為國家土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)(GB15618-1995)所規(guī)定的土壤自然背景值與三級標(biāo)準(zhǔn)。

土壤科學(xué)中可用變異系數(shù)的數(shù)值來對土壤性質(zhì)的變異程度進(jìn)行相應(yīng)等級劃分，變異系數(shù)值在0～15%之間為小變異，16%～50%之間為中等變異，大于50%的則劃分為高度變異[8]，表中數(shù)據(jù)顯示Cu、Pb的變異系數(shù)分別為50%與21%，Cd的變異系數(shù)高達(dá)96%，3種重金屬元素中Cu元素的變異程度較低，Pb變異程度一般，Cd的變異最高，說明該區(qū)域重金屬元素在空間上的變化程度相差較大。該區(qū)域?qū)儆诘V區(qū)，適用于國家三級土壤質(zhì)量標(biāo)準(zhǔn)，從表1中可以看出，Pb，Cu元素含量平均值均未超出標(biāo)準(zhǔn)，但Cd元素嚴(yán)重超標(biāo)，其含量為國家三級土壤質(zhì)量標(biāo)準(zhǔn)的3.5倍，判定該區(qū)域受Cd污染脅迫嚴(yán)重。因此，后續(xù)研究主要針對重金屬Cd元素進(jìn)行光譜變換與建模反演研究。

表1 重金屬含量統(tǒng)計特征

2 光譜變換處理與分析

2.1 光譜變換

光譜變換處理的目的在于減弱光譜由于實驗儀器設(shè)備、測量條件以及各種隨機(jī)因素的影響所產(chǎn)生的噪聲以提升建模精度。Kemper等人的研究發(fā)現(xiàn)，以5 nm、10 nm為重采樣間隔對光譜進(jìn)行重采樣處理均可以在去除光譜噪聲、減少冗余信息的同時較好的保持光譜數(shù)據(jù)特征[9]。本研究為了考慮更多的波段信息，最終在平滑之后選了5 nm的間隔來進(jìn)行重采樣處理。由于本次室內(nèi)實驗數(shù)據(jù)在2 350 nm處附近存在嚴(yán)重噪聲，因此對此波段附近光譜數(shù)據(jù)做剔除處理，剔除前、后效果見圖2。

圖2 剔除前后土壤原始光譜曲線

光譜變換處理的方法較多，由于測量條件的不一所帶來的噪聲與誤差尤為復(fù)雜，沒有一種適用于各物質(zhì)的、各成分分析的變換處理方法，因此在具體分析過程中需盡可能嘗試多種數(shù)學(xué)變換方法以找到適用于某特定研究的最佳處理方法。本研究在對原始光譜數(shù)據(jù)進(jìn)行常規(guī)變換處理方法的同時進(jìn)行了組合變化處理，研究主要選取了一、二階微分與開平方與開立方方法進(jìn)行了結(jié)合，因為微分變換可以有效地消除背景噪聲的影響并且能有效地提取出與研究變量相關(guān)的特征波段，喬星星[10]等的研究發(fā)現(xiàn)涉及微分處理后的光譜變換，尤其是進(jìn)行開方、對數(shù)后再進(jìn)行微分變換處理能最有效的提高土壤氮元素的光譜之間的相關(guān)性。本研究以此為參照在土壤重金屬領(lǐng)域遙感反演領(lǐng)域進(jìn)行嘗試，具體光譜變換處理方法如表2所示。

A.恭喜你，你是一個絕佳的反派人選。雖然未必有人相信你的辦事能力，但你非常有潛力，會有一鳴驚人的一天。你是一個可怕的敵人，有能力做反派主角。

2.2 相關(guān)性分析

相關(guān)系數(shù)的算式為

(1)

式中:R為相關(guān)系數(shù)，R的取值范圍為[-1,1]，|R|越趨近1代表變量X與Y的相關(guān)性越高。

表2 光譜數(shù)據(jù)變換處理方法

相關(guān)性分析是指對兩個或多個具備相關(guān)性的變量元素進(jìn)行分析，從而衡量兩個變量因素的相關(guān)密切程度。相關(guān)性分析在光譜敏感波段的選取方面應(yīng)用廣泛，利用SPSS軟件來分析重金屬元素Cd含量與光譜數(shù)據(jù)的相關(guān)性。將重金屬Cd元素與各預(yù)處理變換后的光譜數(shù)據(jù)進(jìn)行相關(guān)性分析的結(jié)果如圖3所示，各光譜變換處理的最大相關(guān)性數(shù)值如表3所示。

圖3 Cd元素與各光譜處理相關(guān)性曲線圖

重金屬最大相關(guān)性YSFDSDLOG-RECCREOSREOCREOSR-FD Cd-0.2760.4410.4310.2730.321-0.275-0.2740.425EOSR-SDEOCR-FDEOCR-SDREC-FDREC-SDLOG-FDLOG-SD0.4240.4160.421-0.4930.4110.4230.414

分析圖3、表3中數(shù)據(jù)可知：倒數(shù)對數(shù)、開平方、開立方、連續(xù)統(tǒng)去除處理的光譜數(shù)據(jù)與Cd元素的相關(guān)性相比于原始數(shù)據(jù)的相關(guān)性并無改善；一、二階微分方法處理后光譜數(shù)據(jù)與Cd元素的相關(guān)性相比原始光譜數(shù)據(jù)與重金屬Cd元素0.276的相關(guān)性有所提高，相關(guān)性絕對值最大在0.41左右，且二階微分相關(guān)性圖像的整體穩(wěn)定較好；組合變換預(yù)處理方法中只有倒數(shù)一階微分的最大相關(guān)性絕對值接近0.5，其余組合變換如根號二階微分、開立方二階微分變換的最大相關(guān)性絕對值相比常規(guī)一、二階微分處理并無有效提升，且相關(guān)性曲線的整體穩(wěn)定性并不優(yōu)于二階微分相關(guān)性曲線。因此本研究最終選取二階微分、倒數(shù)一階微分處理方法中相關(guān)性較高的波段進(jìn)行進(jìn)一步建模分析并與原始光譜數(shù)據(jù)進(jìn)行比較，備選的特征波段如表4所示。

表4 Cd元素特征波段選取

3 地理加權(quán)回歸建模(GWR)

以往的各種土壤重金屬回歸建模過程中均未考慮到回歸參數(shù)與空間位置的關(guān)系，但在實際研究過程中，數(shù)據(jù)建模的回歸參數(shù)是隨著地理位置而變化的，這時傳統(tǒng)未考慮地理位置的全局回歸模型參數(shù)只是研究區(qū)的平均值，不能反映整個研究區(qū)回歸參數(shù)的真實空間特征。為了解決這一問題，F(xiàn)ortheringham等于1996年提出地理加權(quán)回歸模型(Geographieally Weighted Regression Model,GWR)，其主要改進(jìn)是將研究數(shù)據(jù)的空間位置嵌入到了模型回歸參數(shù)之中，并利用局部加權(quán)最小二乘的方法逐點(diǎn)對各參數(shù)進(jìn)行估計，其中權(quán)值是回歸點(diǎn)所在的地理空間位置到其他各觀測點(diǎn)的地理空間位置之間的距離函數(shù)[11]。

3.1 模型原理

GWR模型是對普通線性回歸模型的擴(kuò)展，將數(shù)據(jù)的地理位置嵌入到回歸參數(shù)之中，即：

yi=β0(ui+vi)+∑βk(ui,vi)xik+εi,

i=1,2,3,…，n.

(2)

式中：ui,vi為第i個采樣點(diǎn)的坐標(biāo)；βk(ui,vi)為第i個采樣點(diǎn)上的第k個回歸參數(shù)；εi是第i個樣點(diǎn)的隨機(jī)誤差。

GWR模型當(dāng)中每個樣本點(diǎn)的回歸權(quán)重值隨著地理位置的變化而變化，權(quán)重函數(shù)用于描述位置因素對建模影響因素的大小[12]。由于big-square法在函數(shù)帶寬范圍之外的權(quán)重為0，更符合實際研究中會對回歸參數(shù)幾乎沒有影響的數(shù)據(jù)刪掉的情形，因此本研究選取big-square函數(shù)作為空間加權(quán)函數(shù)進(jìn)行GWR建模分析，big-square加權(quán)函數(shù)的最優(yōu)帶寬以AIC法則確定：

(3)

由于GWR模型是基于數(shù)據(jù)的空間結(jié)構(gòu)所建立的，不同土壤采樣點(diǎn)的重金屬對光譜的影響參數(shù)是不同的，表現(xiàn)出一定的空間異質(zhì)性，而本研究區(qū)中重金屬Cd含量的變異系數(shù)高達(dá)96%，屬于高度變異，因此研究建立GWR模型進(jìn)行探討分析，其中28個樣點(diǎn)進(jìn)行建模，12個樣點(diǎn)進(jìn)行預(yù)測，地理加權(quán)回歸建?；贕WR4.0軟件實現(xiàn)。

3.2 建模結(jié)果分析

研究采用決定系數(shù)(R2)、均方根誤差(RMSE)、以及相對分析誤差(RPD)3種指標(biāo)對模型的精度進(jìn)行分析。決定系數(shù)也稱擬合優(yōu)度，表示自變量的平方和在因變量平方和中所占的比例，代表模型的所有解釋變量對因變量的聯(lián)合影響程度； MRE主要用于衡量與預(yù)測值同實測值之間的偏離程度，它是將預(yù)測值與實測值偏差的平方與樣本數(shù)比值的平方根再除以樣本的實測平均值，能有效地反映模型的預(yù)測誤差；RPD的定義為實測值標(biāo)準(zhǔn)差與均方根誤差之間的比值，其結(jié)果能夠較好的評價所建立模型的適用泛化能力，模型的R2>0.5并且RPD>2時表明模型的預(yù)測能力較好；1.4

R2=SSR/SST，

(4)

(5)

RPD=STD/RMSE.

(6)

GWR建模結(jié)果顯示(見表5)：①REC-SD光譜數(shù)據(jù)的建模效果最為理想，模型的R2高達(dá)0.82，MRE為29%，表明模型的預(yù)測精度達(dá)到71%，RPD為2.1，表明模型的預(yù)測能力較為優(yōu)良;②SD光譜數(shù)據(jù)的建模結(jié)果中,有MRE為41%，表明模型的預(yù)測精度達(dá)到59%，預(yù)測精度一般，0.67的R2與1.5的RPD表明該模型對新數(shù)據(jù)的泛化能力較好;③YS光譜數(shù)據(jù)的建模效果較差，0.9的RPD與63%的MRE均表明該模型不可用于對樣本點(diǎn)進(jìn)行預(yù)測。

進(jìn)一步分析可知，各光譜變換方式下Cd元素的GWR建模效果差異較大，變換后的光譜數(shù)據(jù)建模效果均優(yōu)于原始光譜數(shù)據(jù)，且以倒數(shù)一階微分變換的建模效果最優(yōu)，說明對于原始光譜數(shù)據(jù)進(jìn)行變換處理，尤其是倒數(shù)一階微分變換處理能夠有效地提高Cd元素GWR模型的預(yù)測精度。將重金屬Cd元素REC-SD光譜數(shù)據(jù)建模的預(yù)測值與實測值繪制的散點(diǎn)圖中顯示(見圖4)，兩者之間具有明顯的線性趨勢，預(yù)測值與實測值趨近于中心Y=X直線，表明樣本點(diǎn)的預(yù)測結(jié)果比較理想。

表5 GWR建模方案與精度評價表

圖4 預(yù)測值與實測值散點(diǎn)圖

4 結(jié)束語

選取瀏陽某礦區(qū)周邊耕地土壤為研究對象，在對研究區(qū)污染狀況進(jìn)行統(tǒng)計分析的基礎(chǔ)上選取受重金屬Cd元素嚴(yán)重脅迫的土壤光譜數(shù)據(jù)進(jìn)行處理與分析。研究結(jié)果顯示，光譜組合變換方式中的REC-SD處理能有效提升重金屬Cd元素含量與光譜數(shù)據(jù)的相關(guān)性，其建立的GWR模型精度均優(yōu)于YS與SD光譜數(shù)據(jù)模型，研究結(jié)果表明光譜組合變換處理的可行性與優(yōu)越性以及高度變異Cd元素GWR建模反演的適用性，為土壤高光譜數(shù)據(jù)的變換處理提供了一定參考價值，也為高度變異元素的反演建模建立了新的思路。但本研究僅采用了較為常見的10種數(shù)學(xué)組合變換，雖取得了一定效果但數(shù)學(xué)組合變換方法不勝枚舉，且受條件限制只研究了Cu、Pb、Cd 3種重金屬元素，其余元素變異對光譜變換的影響亟需進(jìn)一步研究，今后可深入研究以尋求效果更好且尤為普適的變換方法。