基于高光譜定量反演模型的污水綜合水質(zhì)評(píng)價(jià)

陳俊英 邢 正 張智韜 勞聰聰 栗現(xiàn)文 王海峰

(1.西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院， 陜西楊凌 712100；2.西北農(nóng)林科技大學(xué)中國(guó)旱區(qū)節(jié)水農(nóng)業(yè)研究院， 陜西楊凌 712100)

0 引言

近年來(lái)，利用高光譜遙感技術(shù)評(píng)價(jià)和監(jiān)測(cè)水資源水質(zhì)信息狀況方面的研究愈發(fā)深入[1]。而應(yīng)用高光譜技術(shù)檢測(cè)水體水質(zhì)的關(guān)鍵在于水質(zhì)的綜合評(píng)價(jià)和光譜數(shù)據(jù)與水質(zhì)參數(shù)間數(shù)學(xué)模型的建立。對(duì)于光譜數(shù)據(jù)和水質(zhì)參數(shù)間數(shù)學(xué)模型的建立，國(guó)內(nèi)外學(xué)者對(duì)影響水體質(zhì)量的幾個(gè)主要參數(shù)指標(biāo)的遙感估算進(jìn)行了大量的研究，如化學(xué)需氧量(Chemical oxygen demand, COD)、濁度、總磷、生物耗氧量(Biological oxygen demand, BOD)、總氮等。YE等[2]應(yīng)用UVE-SPA-LS-SVM的方法實(shí)現(xiàn)了對(duì)COD的建模預(yù)測(cè)；呂航等[3]利用 HJ-1A衛(wèi)星HSI高光譜遙感數(shù)據(jù)，建立了9個(gè)水質(zhì)參數(shù)與水體光譜反射率之間的估算模型；曹引等[4]建立了水體濁度的高光譜定量反演模型，為水體濁度大面積遙感監(jiān)測(cè)的業(yè)務(wù)化管理提供了技術(shù)支持；BANSOD等[5]通過(guò)高光譜的圖像數(shù)據(jù)，對(duì)恒河的水質(zhì)參數(shù)進(jìn)行了反演評(píng)價(jià)。在對(duì)于水質(zhì)單一參數(shù)的高光譜數(shù)據(jù)反演模型建立上，目前已經(jīng)達(dá)到可以定量的效果[6-9]。但影響水體質(zhì)量的水質(zhì)因子數(shù)目眾多，利用單獨(dú)的某項(xiàng)水質(zhì)參數(shù)來(lái)描述水質(zhì)的信息狀況不夠全面，因此往往會(huì)對(duì)水質(zhì)進(jìn)行綜合評(píng)價(jià)。對(duì)于水質(zhì)的綜合評(píng)價(jià)，目前的研究也較為成熟。馬小雪等[10]利用主成分分析法對(duì)溫瑞塘河流域多項(xiàng)水質(zhì)參數(shù)進(jìn)行時(shí)空分異特征分析和潛在污染源的識(shí)別；徐國(guó)賓等[11]利用模糊標(biāo)識(shí)指數(shù)對(duì)水質(zhì)達(dá)標(biāo)狀況、水質(zhì)類別和主要污染因子進(jìn)行綜合評(píng)價(jià)。這些水質(zhì)綜合評(píng)價(jià)方法通過(guò)已有的水質(zhì)參數(shù)資料能夠很好地表征水質(zhì)的信息狀況，但在實(shí)時(shí)性方面存在不足。因此，需建立高光譜結(jié)合水質(zhì)綜合評(píng)價(jià)方法反演水質(zhì)模型，充分發(fā)揮高光譜遙感的實(shí)時(shí)性、大范圍性與水質(zhì)綜合評(píng)價(jià)方法的全面性、準(zhǔn)確性的優(yōu)勢(shì)。目前對(duì)于高光譜結(jié)合水質(zhì)綜合評(píng)價(jià)方法反演水質(zhì)模型的研究較少。

本文將一組來(lái)自污水處理廠各處理工藝處的水體樣品分成兩份，一份經(jīng)由室內(nèi)理化試驗(yàn)，檢測(cè)各項(xiàng)水質(zhì)參數(shù)，并利用主成分分析對(duì)水體水質(zhì)進(jìn)行綜合評(píng)價(jià)，得到水質(zhì)綜合評(píng)價(jià)因子；同時(shí)對(duì)另一份水體樣品進(jìn)行高光譜數(shù)據(jù)的采集，將采集到的光譜數(shù)據(jù)進(jìn)行不同的預(yù)處理，采用偏最小二乘法、逐步回歸法和極限學(xué)習(xí)機(jī)法對(duì)光譜數(shù)據(jù)和水質(zhì)綜合評(píng)價(jià)因子進(jìn)行建模預(yù)測(cè)以及驗(yàn)證。比較各預(yù)處理方法及對(duì)應(yīng)的建模方法的驗(yàn)證結(jié)果，選出更適合用于水質(zhì)綜合評(píng)價(jià)高光譜反演的數(shù)據(jù)預(yù)處理及對(duì)應(yīng)建模方法，為建立高光譜結(jié)合水質(zhì)綜合評(píng)價(jià)方法反演水質(zhì)模型，實(shí)現(xiàn)對(duì)水質(zhì)信息狀況的大范圍實(shí)時(shí)監(jiān)測(cè)提供可行的路徑。

1 試驗(yàn)材料與方法

1.1 樣本采集及樣本水質(zhì)指標(biāo)的測(cè)定

試驗(yàn)用水水樣取自某生活污水處理廠，取水位置分別為生活污水處理的不同工藝處，即進(jìn)水口、厭氧池、好氧池、沉淀池、出水口，對(duì)照的水樣為純凈水。各水樣的各項(xiàng)水質(zhì)參數(shù)經(jīng)由室內(nèi)理化試驗(yàn)測(cè)定，結(jié)果見表1(部分)。

1.2 光譜測(cè)定

污水樣品采用ASD Field Spec 3型地物光譜儀測(cè)量高光譜數(shù)據(jù)。光譜儀波長(zhǎng)范圍為350～2 500 nm，采樣間隔為1.4 nm(350～1 000 nm)和2 nm(1 000～2 500 nm)，重采樣間隔為1 nm。光譜測(cè)量在暗室中進(jìn)行，光源為DH-2000型氘鹵鎢燈光源[12]。

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

本試驗(yàn)中87個(gè)樣品獲得的光譜波段為350～2 500 nm。由于試驗(yàn)條件以及其他因素的影響，測(cè)量的光譜中可能包含了一些冗余信息以及噪聲，因此需要對(duì)獲得的光譜波段進(jìn)行選擇以提高建模的準(zhǔn)確度。

由圖1可看到，在350～400 nm、2 300～2 500 nm波段范圍由于處于邊緣噪聲較大，不適用于建模。而大于2 000 nm波段，反射率很小，可利用的信息很少，難以找出不同樣本光譜圖的差別。故本研究選用400～2 000 nm的光譜波段。在建模前需要對(duì)光譜數(shù)據(jù)進(jìn)行一定的預(yù)處理以削弱由測(cè)試環(huán)境及其他干擾因素導(dǎo)致的影響，提高數(shù)據(jù)信噪比。本文采用預(yù)處理方法有Savitzky-Golay(SG)平滑、標(biāo)準(zhǔn)正態(tài)化(SNV)、去包絡(luò)線(CR)和倒數(shù)之對(duì)數(shù)(LR)預(yù)處理等。

1.3.1Savitzky-Golay平滑處理

平滑濾波是光譜分析中常用的預(yù)處理方法之

表1 主要水質(zhì)參數(shù)Tab.1 Main water quality parameters

圖1 全波段原始光譜反射率曲線Fig.1 Full-band original spectral reflectance curves

一，通常利用Savitzky-Golay方法進(jìn)行平滑濾波。Savitzky-Golay方法是一種在時(shí)域內(nèi)基于局域多項(xiàng)式最小二乘法擬合的濾波方法。其最大的特點(diǎn)在于在濾除噪聲的同時(shí)可以確保信號(hào)的形狀、寬度不變，可以提高光譜的平滑性，并降低噪聲的干擾[13]。本次處理中移動(dòng)窗口寬度為5及多項(xiàng)式次數(shù)為3。

利用Savitzky-Golay濾波法對(duì)所有光譜數(shù)據(jù)(400～2 000 nm)進(jìn)行平滑、去噪，取去噪聲后的部分光譜數(shù)據(jù)曲線如圖2a所示。

1.3.2標(biāo)準(zhǔn)正態(tài)化處理

標(biāo)準(zhǔn)正態(tài)變量變換(SNV)預(yù)處理針對(duì)一條光譜進(jìn)行處理，主要是消除光程變化、散射和顆粒大小之間的干擾[13]。計(jì)算公式為

(1)

式中Xi——i個(gè)樣品光譜的平均值(標(biāo)量)

m——波長(zhǎng)點(diǎn)數(shù)

n——校正集樣品數(shù)

采用標(biāo)準(zhǔn)正態(tài)變量變換(SNV)對(duì)所有光譜數(shù)據(jù)(400～2 000 nm)進(jìn)行處理，取處理后的部分光譜數(shù)據(jù)曲線如圖2b所示。

1.3.3去包絡(luò)線處理

去包絡(luò)線(CR)處理將光譜曲線歸一化到0～1之間，能夠突出光譜曲線的吸收和反射特征，增強(qiáng)光譜曲線各波段之間的對(duì)比性[14]。

采用去包絡(luò)線(CR)對(duì)所有光譜數(shù)據(jù)(400～2 000 nm)進(jìn)行處理，取處理后的部分光譜數(shù)據(jù)曲線如圖2c所示。

1.3.4倒數(shù)之對(duì)數(shù)處理

在高光譜研究中，常將反射率進(jìn)行倒數(shù)之對(duì)數(shù)變換,該變換形式有利于處理非線性問(wèn)題，增強(qiáng)相似光譜之間的差異，并適當(dāng)減少隨機(jī)誤差[15]。

采用倒數(shù)之對(duì)數(shù)處理方法對(duì)所有光譜數(shù)據(jù)(400～2 000 nm)進(jìn)行處理，取處理后的部分光譜數(shù)據(jù)曲線如圖2d所示。

平滑光譜反射率SG、LR在ViewSpec Pro V6.0.11軟件中處理獲得，指標(biāo)CR利用ENVI 5.1的Continuum Removed模塊處理得到。其他數(shù)據(jù)預(yù)處理通過(guò)The Unscrambler X 10.4實(shí)現(xiàn)。

1.4 模型建立與驗(yàn)證

1.4.1樣品集的劃分

樣品集的劃分采用Kennard-Stone算法(簡(jiǎn)稱K-S)。K-S算法是根據(jù)已經(jīng)被選擇的樣品計(jì)算未被選擇的樣品的最小歐氏距離，然后通過(guò)選擇經(jīng)由計(jì)算的歐氏距離最大的樣品進(jìn)入校正集，以此反復(fù)，直至選出的樣品數(shù)達(dá)到指定要求[16]。K-S算法在選擇具有代表性的樣品方面已經(jīng)被證明有著較好的效果[17]。本研究選取58個(gè)水質(zhì)樣本作為建模集，29個(gè)水質(zhì)樣本作為驗(yàn)證集，分別用于模型的建立以及精度驗(yàn)證。

圖2 不同預(yù)處理后的光譜反射率曲線Fig.2 Spectral reflectance curves after different pretreatments

1.4.2模型方法和評(píng)價(jià)指標(biāo)

采用偏最小二乘回歸法(PLSR)、逐步回歸法(SR)和極限學(xué)習(xí)機(jī)(ELM)這3種回歸方法建立高光譜遙感對(duì)水質(zhì)綜合評(píng)價(jià)的反演模型。其中PLSR在建模過(guò)程中具有降維、信息集成和波段優(yōu)選等方法極大地提高了系統(tǒng)提取主成分的能力，得到了廣泛地應(yīng)用，可較好地解決自變量之間存在的共線性問(wèn)題[18-20]。而SR是一種便捷高效的模型優(yōu)化方法，對(duì)高光譜數(shù)據(jù)的“降維”具有較好的作用[21-22]。ELM是一類基于前饋神經(jīng)網(wǎng)絡(luò)的機(jī)器學(xué)習(xí)算法，與傳統(tǒng)的前饋神經(jīng)網(wǎng)絡(luò)相比較，ELM 有著學(xué)習(xí)效率高、精度高且參數(shù)調(diào)整簡(jiǎn)單等優(yōu)點(diǎn)[23-25]。

(2)

(3)

(4)

式中yi、i——驗(yàn)證樣本的觀測(cè)值和預(yù)測(cè)值

n——驗(yàn)證樣本數(shù)

SD——樣本觀測(cè)值方差

RMSE——均方根誤差

2 結(jié)果與分析

2.1 不同污水樣本高光譜曲線特征分析

選12條較為典型的污水樣本經(jīng)過(guò)平滑(SG)、倒數(shù)之對(duì)數(shù)(LR)、去包絡(luò)線(CR)、標(biāo)準(zhǔn)正態(tài)化(SNV)4種預(yù)處理后的光譜曲線，見圖2。

由圖2a可以發(fā)現(xiàn)，12條光譜曲線波形基本相似。圖2b、2c、2d分別是水體原始光譜經(jīng)標(biāo)準(zhǔn)正態(tài)化(SNV)、去包絡(luò)線(CR)、倒數(shù)之對(duì)數(shù)(LR)3種預(yù)處理后的反射率。從圖2a、2b可以看出，在410、480 nm波段處有明顯的反射峰值，在440、760、900、1 000 nm處有明顯的吸收峰值。這與呂航等[3]的研究較為符合，是由于在410、470 nm波段處有較多的水質(zhì)參數(shù)對(duì)反射更為敏感，而在440、760、900、1 000 nm處有較多的水質(zhì)參數(shù)對(duì)光譜的吸收更為敏感。從圖2c、2d可以看到，經(jīng)由去包絡(luò)線CR和倒數(shù)之對(duì)數(shù)LR處理后，放大了1 000～2 000 nm處的光譜波段，使得光譜吸收帶變得更加明顯，可以看到在1 400 nm附近處和1 900 nm附近處也有著光譜的吸收敏感波段區(qū)，也驗(yàn)證了CR預(yù)處理能夠突出光譜曲線的吸收和反射特征，增強(qiáng)光譜曲線各波段之間的對(duì)比性，以及 LR預(yù)處理在增強(qiáng)相似光譜之間差異的優(yōu)越性。

2.2 水質(zhì)綜合評(píng)價(jià)

水質(zhì)系統(tǒng)是一個(gè)由各種水質(zhì)污染指標(biāo)變量組成的復(fù)雜系統(tǒng)，其內(nèi)里蘊(yùn)含眾多能夠影響水體質(zhì)量的水質(zhì)因子，每一種水質(zhì)因子都只從某一方面表征了水體質(zhì)量，而這些水質(zhì)因子之間又往往有著不同程度的相關(guān)性，因此在對(duì)水質(zhì)進(jìn)行評(píng)價(jià)時(shí)，將這些水質(zhì)因子都利用起來(lái)有一定的難度[27]。

主成分分析(Principal component analysis，PCA)是一種多元統(tǒng)計(jì)分析方法，其主要是利用降維的思想，把最初錯(cuò)綜復(fù)雜的多個(gè)變量歸納總結(jié)成為少數(shù)幾個(gè)綜合變量，其中每一個(gè)綜合變量都是原始變量的線性組合，各綜合變量之間不存在相關(guān)性，從而實(shí)現(xiàn)利用少數(shù)幾個(gè)綜合變量來(lái)反映原始變量的絕大部分信息，且所含的信息互不重疊[28]。

主成分分析法(PCA)分析水質(zhì)參數(shù)的基本思想是將n個(gè)水體樣本的P個(gè)污染指標(biāo)變量Xj(j=1，2，…，P)，通過(guò)對(duì)這P個(gè)污染指標(biāo)變量相關(guān)性的研究，重新組合構(gòu)造出m(m

2.3 PLSR建模及預(yù)測(cè)

偏最小二乘回歸PLSR具有主成分分析、典型相關(guān)分析和多元線性回歸等的優(yōu)點(diǎn)。運(yùn)用The Unscrambler X 10.4軟件將全波段(400～2 000 nm)4種光譜指標(biāo)(SG、LR、SNV、CR)作為自變量，以水質(zhì)綜合評(píng)價(jià)因子作為因變量，通過(guò)將均方根誤差(RMSE)和決定系數(shù)R2對(duì)主因子數(shù)作圖的方法確定最佳主因子數(shù)，建立PLSR回歸模型。建模以及驗(yàn)證結(jié)果見表4。

表2 特征值貢獻(xiàn)率Tab.2 Eigenvalue contribution rate

表3 水質(zhì)分級(jí)標(biāo)準(zhǔn)Tab.3 Water quality classification standards

表4 水質(zhì)指標(biāo)的偏最小二乘模型Tab.4 PLSR model of water quality indicators

2.4 SR建模及預(yù)測(cè)

逐步回歸是一種線性回歸模型自變量選擇方法，其基本思想是將自變量逐個(gè)引入，根據(jù)自變量對(duì)因變量的解釋程度或顯著性，將對(duì)因變量解釋程度小或者顯著性低的自變量進(jìn)行剔除，保留顯著的解釋變量，如此反復(fù)，直到既沒有顯著的解釋變量選入回歸方程，也沒有不顯著的解釋變量從回歸方程中剔除為止，完成對(duì)數(shù)據(jù)的大幅降維，得到了最優(yōu)的解釋變量集。本文運(yùn)用全波段4種光譜指標(biāo)(SG、LR、SNV、CR)作為自變量，水質(zhì)綜合評(píng)價(jià)因子為因變量。變量入選和剔除的顯著水平分別設(shè)為0.15和0.25，由“最優(yōu)”解釋變量集所建立模型的結(jié)果見表5。

表5 水質(zhì)指標(biāo)的SR模型Tab.5 SR model of water quality indicators

由表5可以看出，逐步回歸方法通過(guò)對(duì)波段的“篩選”，剔除了大量對(duì)水質(zhì)綜合評(píng)價(jià)因子不顯著的波段數(shù)據(jù)，僅保留了原數(shù)據(jù)約1%的顯著波段數(shù)據(jù)，其降維效果非常顯著。在SR建模對(duì)數(shù)據(jù)的“降維”過(guò)程中，原始光譜SG保留的波段數(shù)目最多，表明在“降維”過(guò)程中，原始光譜的波段和水質(zhì)綜合評(píng)價(jià)因子間具有較好相關(guān)性；而LR處理后的光譜數(shù)據(jù)保留的波段數(shù)目最少，表明在“降維”過(guò)程中，原始光譜的波段數(shù)值和水質(zhì)綜合評(píng)價(jià)因子間具有較差的相關(guān)性。

2.5 ELM建模及預(yù)測(cè)

極限學(xué)習(xí)機(jī)(ELM)是由HUANG等[32]提出來(lái)的求解單隱層神經(jīng)網(wǎng)絡(luò)的算法。ELM的網(wǎng)絡(luò)訓(xùn)練模型由輸入層、隱含層和輸出層組成。其中，模型的訓(xùn)練效果受隱含層的神經(jīng)元數(shù)量影響較大，且隱含層的神經(jīng)元數(shù)量需人為確定。輸入層和輸出層的神經(jīng)元數(shù)量取決于所分析問(wèn)題的自變量和因變量數(shù)量。具體推導(dǎo)過(guò)程及訓(xùn)練步驟詳見文獻(xiàn)[33]。ELM最大的特點(diǎn)是可以在保證學(xué)習(xí)精度的前提下相對(duì)于傳統(tǒng)的神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)算法速度更快。

以經(jīng)過(guò)逐步回歸SR降維后的4種光譜指標(biāo)(SG、LR、SNV、CR)的光譜數(shù)據(jù)作為自變量，水質(zhì)綜合評(píng)價(jià)因子為因變量建立SR-ELM模型。建模以及驗(yàn)證結(jié)果見表6。

表6 水質(zhì)指標(biāo)的極限學(xué)習(xí)機(jī)模型Tab.6 ELM model of water quality indicators

作為機(jī)器學(xué)習(xí)的一種算法，ELM在本組數(shù)據(jù)建模中，隱含層的神經(jīng)元數(shù)量對(duì)訓(xùn)練效果的影響較大。從表6可以看出，建模預(yù)測(cè)效果較好的SG-SR-ELM模型和SNV-SR-ELM模型的相對(duì)最佳隱含層單元數(shù)高于建模預(yù)測(cè)效果較差的LR-SR-ELM模型和CR-SR-ELM模型，反映了SG和SNV數(shù)據(jù)和水質(zhì)綜合評(píng)價(jià)因子間更具有相關(guān)性?？傮w來(lái)說(shuō)，ELM模型在本組數(shù)據(jù)中對(duì)于光譜的擬合和預(yù)測(cè)具有很好的效果。

2.6 模型對(duì)比

運(yùn)用3種不同的回歸方法對(duì)4種光譜指標(biāo)進(jìn)行建模，各個(gè)模型對(duì)污水水質(zhì)的預(yù)測(cè)效果見圖3。

圖3 水質(zhì)綜合評(píng)價(jià)反演模型Fig.3 Inversion model for comprehensive evaluation of water quality

2.7 討論

高光譜遙感在定量反演水質(zhì)參數(shù)時(shí)，由于具有光譜分辨率高和波段連續(xù)性強(qiáng)等特點(diǎn)，可以獲得更為全面廣泛的光譜波段數(shù)據(jù)，而由于光譜測(cè)量中的某些人為和自然因素的干擾，光譜數(shù)據(jù)需要進(jìn)行不同的數(shù)學(xué)變換以增強(qiáng)信噪比，從而提高光譜數(shù)據(jù)與水質(zhì)參數(shù)的相關(guān)性，進(jìn)而提高模型的預(yù)測(cè)精度[22]。

雖然高光譜遙感在實(shí)際應(yīng)用中，可以獲得更為精細(xì)的光譜信息，但因此也造成了數(shù)據(jù)和計(jì)算量的冗雜，為實(shí)現(xiàn)數(shù)據(jù)的篩選和模型的簡(jiǎn)化，本文通過(guò)逐步回歸方法對(duì)光譜數(shù)據(jù)進(jìn)行分析篩選建模。結(jié)果表明，逐步回歸方法能夠有效地對(duì)數(shù)據(jù)進(jìn)行大幅度的降維(可達(dá)到99%)，使得篩選后留下的數(shù)據(jù)既是重要的，又沒有嚴(yán)重的多重共線性。同時(shí)SR模型也有著很好的精度和預(yù)測(cè)效果，而以經(jīng)過(guò)逐步回歸篩選后的波段所建立的SR-ELM模型的精度和預(yù)測(cè)效果明顯比PLSR模型和SR模型更優(yōu)越，這為下一步的對(duì)高光譜數(shù)據(jù)通過(guò)波段篩選和數(shù)據(jù)降維以實(shí)現(xiàn)模型的簡(jiǎn)化提供了可行性。由水質(zhì)綜合評(píng)價(jià)因子和光譜數(shù)據(jù)建立的水質(zhì)綜合反演模型較由單項(xiàng)水質(zhì)參數(shù)建立的反演模型精度有所下降，沒有能夠很好地表征各不同的單項(xiàng)水質(zhì)參數(shù)對(duì)于光譜數(shù)據(jù)各波段的影響和作用，有待進(jìn)一步探索各水質(zhì)參數(shù)對(duì)于光譜不同波段數(shù)據(jù)的單獨(dú)影響和作用以及綜合的交互影響和作用。另外，由于內(nèi)陸水體光學(xué)特征的復(fù)雜性、水質(zhì)影響因子的多樣性，如何更好地將水質(zhì)信息狀況的綜合評(píng)價(jià)與高光譜技術(shù)相結(jié)合，以達(dá)到通過(guò)高光譜遙感技術(shù)實(shí)時(shí)全面地反映水體的污染程度，主要污染物的類別、來(lái)源、成因、時(shí)空分布規(guī)律以及變化趨勢(shì)，值得去進(jìn)一步深入地研究和探索，是下一步研究的方向和目標(biāo)。

3 結(jié)論

(1)PLSR模型、SR模型、SR-ELM模型均能得到具有很好精度和預(yù)測(cè)效果的模型，其中SR-ELM模型的精度高于其他兩個(gè)模型，更適用于處理本組的水質(zhì)綜合評(píng)價(jià)因子。

(3) ELM模型預(yù)測(cè)水質(zhì)綜合評(píng)價(jià)因子時(shí)具有很好的效果，為機(jī)器學(xué)習(xí)在水質(zhì)參數(shù)反演預(yù)測(cè)方面的應(yīng)用驗(yàn)證了可行性，其中SNV-SR-ELM模型為水質(zhì)的綜合評(píng)價(jià)方法和高光譜反演模型的結(jié)合提供了參考。

(4)基于標(biāo)準(zhǔn)正態(tài)化變換指標(biāo)建立的偏最小二乘模型、逐步回歸模型與極限學(xué)習(xí)機(jī)模型，其決定系數(shù)和RPD均最高，反演精度最優(yōu)，SNV為本組光譜數(shù)據(jù)的最佳預(yù)處理方法。其中SNV-SR-ELM模型決定系數(shù)為0.954，RPD為4.651，為本組數(shù)據(jù)的最佳模型。

(5)水質(zhì)綜合評(píng)價(jià)方法和高光譜技術(shù)反演水質(zhì)參數(shù)的結(jié)合具有可行性，其模型的建立可以為水體的快速監(jiān)測(cè)和綜合評(píng)價(jià)提供參考。