袁 婕，張 飛,3※，葛翔宇，郭婉臻，鄧來飛

地理加權(quán)回歸模型結(jié)合高光譜反演鹽生植物葉片鹽離子含量

袁 婕1,2，張 飛1,2,3※，葛翔宇1,2，郭婉臻1,2，鄧來飛1,2

（1. 新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院，烏魯木齊 830046；2. 新疆大學(xué)綠洲生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室，烏魯木齊 830046；3. 新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院智慧城市與環(huán)境建模自治區(qū)普通高校重點(diǎn)實(shí)驗(yàn)室，烏魯木齊 830046）

快速、無損地估算鹽生植物葉片鹽離子含量在植物生長(zhǎng)監(jiān)測(cè)、耐鹽植物篩選和土壤鹽漬化監(jiān)測(cè)等方面有實(shí)用價(jià)值。該研究以新疆艾比湖保護(hù)區(qū)內(nèi)鹽生植物為研究對(duì)象，通過分析植物葉片鹽離子（K+、Na+、Ca2+、Mg2+）含量與冠層高光譜數(shù)據(jù)的光譜變換和二維植被指數(shù)(比值型植被指數(shù)（ratio vegetation index，RVI）、差值型植被指數(shù)（difference vegetation index，DVI）、歸一化型植被指數(shù)（normalized difference vegetation index，NDVI）)的相關(guān)性選取特征波段，構(gòu)建基于地理加權(quán)回歸模型（geographically weighted regression，GWR）的葉片鹽離子含量估算模型，并與BP神經(jīng)網(wǎng)絡(luò)模型（back propagation neural network）進(jìn)行對(duì)比，研究基于GWR模型估算干旱區(qū)鹽生植物葉片鹽離子的可行性。結(jié)果表明，選取特征波段集中表現(xiàn)在紅及短波紅外波段：K+含量在反射率倒數(shù)的對(duì)數(shù)選取的紅光區(qū)域內(nèi)波段使用GWR估算效果最佳；Na+的特征波段在光譜變換下集中于短波紅外區(qū)域，二維植被指數(shù)集中在近紅外、短波近紅外及黃、橙、紅區(qū)域，各種波段選取下GWR對(duì)Na+的含量估算均有較好效果，但反射率對(duì)數(shù)的一階估算效果最好；Ca2+含量在反射率平方根的一階微分下選取的短波紅外波段通過GWR模型估算效果最好；Mg2+含量在DVI選取的位于紅光區(qū)域特征波段估算效果最佳，但使用GWR模型對(duì)Mg2+的估算精度不及BP模型。分析基于GWR鹽離子模型估算模型發(fā)現(xiàn)，含量較高的離子估算效果更好，K+、Na+的模型精度優(yōu)于Ca2+、Mg2+。在使用GWR模型估算植物葉片鹽離子含量時(shí)，特征波段均指向紅及短波紅外波段，符合植被光譜機(jī)理的響應(yīng)。

干旱；葉片；高光譜；GWR模型；鹽生植物；鹽離子

0 引 言

植被是生態(tài)系統(tǒng)中的重要組成部分，研究區(qū)域植被覆蓋變化對(duì)衡量區(qū)域生態(tài)系統(tǒng)健康、合理利用植物資源和城市規(guī)劃具有重要意義[1-2]。鹽生植物對(duì)土壤中的鹽分具有一定的吸收作用，生長(zhǎng)在鹽漬土環(huán)境下下的鹽生植物對(duì)降低土壤鹽離子含量具有很好的效果[3]。新疆是中國(guó)鹽漬化土壤分布分布范圍最廣、面積最大的區(qū)域，鹽漬化土壤占全區(qū)土地總面積的7%。土壤鹽漬化是該區(qū)域重大環(huán)境風(fēng)險(xiǎn)嚴(yán)重影響生態(tài)安全和人類發(fā)展。而鹽生植物對(duì)鹽漬化土壤有很好的改良效果，可作為作為自然調(diào)節(jié)器。鹽生植物可以有效降低土壤含鹽量，疏松土壤，能夠緩解鹽漬化災(zāi)害，恢復(fù)退化的生態(tài)系統(tǒng)[4-6]。鹽生植物葉片各鹽離子含量是反映植物生理狀況和土壤鹽漬化程度的重要指標(biāo)和篩選耐鹽植物的重要參數(shù)。然而，目前檢測(cè)鹽生植物葉片各鹽離子含量需要破壞性取樣和復(fù)雜的化學(xué)分析，開發(fā)快速、無損的估算鹽生植物葉片各鹽離子含量的技術(shù)具有多方面的實(shí)用價(jià)值。

高光譜數(shù)據(jù)近年來得到廣泛應(yīng)用，尤其在定量估算地物參量中發(fā)揮出其潛力[7]。高光譜數(shù)據(jù)具較高的光譜分辨率，所攜帶的精細(xì)光譜信息將植物的遙感知識(shí)從宏觀監(jiān)測(cè)帶到生理生化過程的微觀識(shí)別[8]。植物高光譜在植物生理指標(biāo)葉綠素與冠層葉綠素密度、全氮、粗纖維、植物干物質(zhì)量等方面均有研究[9-12]。Wang等[13]通過植物光譜反射率得到葉片尺度的蒸騰量，并找出其敏感波段在2 435、2 440、2 445和2 470 nm處，這對(duì)偏最小二乘回歸法（partial least squares regression, PLSR）預(yù)測(cè)葉片蒸騰作用至關(guān)重要。Sampson等[14-15]將高光譜用于植物的病蟲害及脅迫檢測(cè)，并在疾病發(fā)生初期得到很好的檢測(cè)效果；Guo[16]研究表明反演濕地植物氮含量時(shí)，植被類型差異對(duì)反演模型影響不顯著。然而，干旱區(qū)具有其區(qū)域的特點(diǎn)，該區(qū)域內(nèi)鹽生植物的光譜信息對(duì)各鹽離子敏感程度尚不明確?？v觀近年研究，光譜的預(yù)測(cè)多采用多元線性回歸、BP神經(jīng)網(wǎng)絡(luò)、偏最小二乘、隨機(jī)森林法等[17-19]，并得到了較好的預(yù)測(cè)效果，但這些預(yù)測(cè)方法也具有一定的局限性：默認(rèn)每個(gè)樣點(diǎn)的環(huán)境因素對(duì)光譜反射率影響是相同的，即模型系數(shù)相同[20-21]。因此在考慮地理要素時(shí)，應(yīng)當(dāng)充分考慮要素的空間異質(zhì)性，每個(gè)樣點(diǎn)代入不同的系數(shù)即加入空間坐標(biāo)進(jìn)行分析。

地理加權(quán)回歸（geographically weighted regression, GWR）模型是針對(duì)不同空間子集受空間變化影響的自變量與響應(yīng)變量之間的關(guān)系構(gòu)建模型，被廣泛應(yīng)用于具有空間非平穩(wěn)性特征的空間數(shù)據(jù)領(lǐng)域。目前已在土壤屬性的空間預(yù)測(cè)中有顯著效果，然而在受空間關(guān)系約束的植被光譜方面尚待挖掘?；诖?，本文充分考慮空間因素，基于GWR模型定量估算鹽離子含量；并闡明葉片鹽離子與相關(guān)光譜參量間的關(guān)系，優(yōu)選出特征波段，旨在為后續(xù)低空遙感系統(tǒng)在葉片鹽離子估算應(yīng)用提供理論依據(jù)及技術(shù)支持。

1 研區(qū)概況

新疆艾比湖濕地國(guó)家級(jí)自然保護(hù)區(qū)范圍位于82°35′47″～83°53′21″E，44°31′05″～45°09′35″N（圖1），東西長(zhǎng)102.63 km，南北寬72.3 km[22]。艾比湖保護(hù)區(qū)景觀多樣，鹽堿化土壤面積大，鹽生植物種類多樣，常見的鹽生植物有白刺、駱駝刺、堿蓬、鹽節(jié)木、鹽爪爪、鹽穗木、鹽穗木等[23-24]。艾比湖國(guó)家自然保護(hù)區(qū)濕地的地理位置和生態(tài)位置十分重要，屬于典型的溫帶干旱區(qū)濕地，是天山北坡綠洲與沙漠化共軛演化的中心，對(duì)地區(qū)調(diào)節(jié)氣候、維護(hù)區(qū)域生態(tài)平衡具有重要意義[25]。

圖1 研究區(qū)及采樣點(diǎn)示意圖

2 數(shù)據(jù)采集與處理

2.1 植物光譜采集與預(yù)處理

2017年7月，在研究區(qū)進(jìn)行光譜采集。共設(shè)30個(gè)樣點(diǎn)，每個(gè)樣點(diǎn)選取2～3個(gè)優(yōu)勢(shì)種（主要包括花花柴、梭梭、檉柳等）進(jìn)行光譜及葉片采集。采集到的植被葉片放入牛皮紙袋中，立即用冰保存，以確保新鮮。使用ASD（Analytical Spectral Devices）公司生產(chǎn)FieldSpec3光譜儀進(jìn)行光譜采集，光譜采樣間隔為1.4 nm（采樣范圍350～1 000 nm）和2 nm（采樣范圍1 000～2 500 nm），重采樣間隔1 nm。測(cè)前使用白板進(jìn)行定標(biāo)校正，在晴朗無風(fēng)的正午進(jìn)行采集，光譜采集時(shí)間為12:00—15:00陽光幾乎直射的時(shí)間段。為盡量降低背景物影響，光譜儀探頭垂直放置在植被冠層上方約5 cm處，向下對(duì)準(zhǔn)被測(cè)植物，如植物較為稀疏，則使其盡量聚集以確保充滿整個(gè)視場(chǎng)。每點(diǎn)測(cè)量10個(gè)反射光譜，取平均值作為該點(diǎn)的原始光譜反射率。為減小野外噪聲對(duì)光譜數(shù)據(jù)的影響，對(duì)測(cè)得光譜曲線進(jìn)行均值、去噪及平滑處理[26]。利用光譜處理軟件ASD View Spec Pro對(duì)采集到的光譜曲線進(jìn)行分析。將校正后的光譜和平均值作為光譜樣本集，進(jìn)行光譜變換和統(tǒng)計(jì)分析。

2.2 葉片含鹽量計(jì)算

在野外每個(gè)樣點(diǎn)采集2～3種典型植物冠層葉片若干，剪下后裝入牛皮紙袋帶回。植物葉片樣在105 ℃殺青30 min后，70 ℃烘至恒量，再將樣品坩堝放在調(diào)溫電爐里進(jìn)行加熱，待樣品冒煙后，再燒15 min左右。將裝有樣品的坩堝移入高溫電爐中，半開坩堝蓋，由室溫升至400 ℃，保持30 min，再升至550 ℃，燒至灰分近于白色為止，冷卻后稱至恒質(zhì)量（準(zhǔn)確至0.01 g）。取一定量的灰分經(jīng)粉碎過篩后用HNO3-HCLO4溶液定容，用北京普析TAS-986型原子吸收分光光度計(jì)（分辨率0.3 nm）測(cè)定K+、Na+、Ca2+和Mg2+的含量[27]。

2.3 光譜變換與特征波段選取

光譜降維主要以波段優(yōu)選為主，選取標(biāo)準(zhǔn)應(yīng)具有相關(guān)性好、噪聲弱、保留完整、信息量大的的特點(diǎn)[28]。由于原始光譜為野外采集獲得，光譜噪聲多集中于近紅外區(qū)域，如圖2a所示，故統(tǒng)一刪除了近紅外區(qū)域1 356～ 1 479、1 801～1 999、2 341～2 500 nm 3個(gè)區(qū)域的波段，見圖2b。波段選取采用2種不同方法：1）為了突顯特征光譜的有效信息，對(duì)原始光譜進(jìn)行14種光譜變換：一階微分、二階微分、求平方根、平方根的一階微分、平方根的二階微分、求對(duì)數(shù)、對(duì)數(shù)一階、對(duì)數(shù)二階、求倒數(shù)、倒數(shù)一階、倒數(shù)二階、求倒數(shù)對(duì)數(shù)、倒數(shù)對(duì)數(shù)一階、倒數(shù)對(duì)數(shù)二階。將植物葉片鹽離子含量與變換后的光譜建立相關(guān)性，選出相關(guān)性最大的波段作為建模的光譜參數(shù)；2）考慮到光譜波段的協(xié)同響應(yīng)，利用更多有效的光譜信息，以消除土壤的噪聲干擾。

圖2 原始光譜和去干擾噪聲后光譜

選取敏感波段組合構(gòu)建3種二維植被指數(shù)[29]：歸一化型植被指數(shù)（normalized difference vegetation index，NDVI）、差值型植被指數(shù)（difference vegetation index，DVI）、比值型植被指數(shù)（ratio vegetation index，RVI）。

式中R和R為波段和獲取的光譜反射率。

2.4 模型方法

GWR對(duì)普通最小二乘回歸模型進(jìn)行了空間擴(kuò)展，將數(shù)據(jù)的地理位置嵌入到回歸參數(shù)之中，使得參數(shù)可以進(jìn)行局部估計(jì)，擴(kuò)展后的模型如下：

式中y為樣點(diǎn)的因變量；x為第個(gè)樣點(diǎn)上的第個(gè)變量（共個(gè)）的觀測(cè)值；為樣點(diǎn)總數(shù)；（μ，v）為樣點(diǎn)的地理空間坐標(biāo)；0為回歸常數(shù)項(xiàng)，β為第個(gè)回歸系數(shù)；?為誤差項(xiàng)。如果β在空間保持不變，則模型就變?yōu)槿帜Ｐ?，系?shù)估算采用加權(quán)最小二乘實(shí)現(xiàn)，用矩陣表示為

式中W為由已知點(diǎn)估計(jì)待測(cè)點(diǎn)時(shí)的權(quán)重，d為估算點(diǎn)與樣點(diǎn)間的歐氏距離，為帶寬其中帶寬由最小 AIC信息準(zhǔn)則確定[30-31]。

本文的GWR回歸過程在GWR4.0軟件支持下完成， GWR模型的預(yù)測(cè)能力與應(yīng)用較為廣泛應(yīng)用的BP神經(jīng)網(wǎng)絡(luò)模型（back propagation neural network）進(jìn)行了對(duì)比。模型精度評(píng)價(jià)采用決定系數(shù)（2）和均方根誤差（root mean squared error，RMSE）進(jìn)行評(píng)價(jià)。

3 結(jié)果與分析

3.1 葉片鹽離子含量統(tǒng)計(jì)分析

K、Na、Ca、Mg是植物生長(zhǎng)所需的大量元素，均為金屬元素并以離子的形式存在。葉片鹽離子含量如圖3所示，從均值來看，在艾比湖保護(hù)區(qū)采集的干旱區(qū)鹽生植物中：Na+含量最高，K+次之，而Ca2+與Mg2+含量最低。在后續(xù)建模中以3∶1比例選取建模集與驗(yàn)證集。

圖3 鹽離子含量箱圖

3.2 一維光譜相關(guān)性特征波段選擇

將原光譜及14種數(shù)學(xué)變換后的光譜分別與鉀、鈉、鈣、鎂4種鹽離子進(jìn)行相關(guān)性分析結(jié)果如圖4所示。

注：R為反射率，下同。樣本數(shù)64。Note：R is reflectance. Same as below. Sample size is 64.

優(yōu)選出的4種離子建模所需的特征波段如表1所示，K+的敏感波段集中在在400～700 nm的光合有效輻射區(qū)域（photosynthetically available radiation，PAR），且集中在不經(jīng)微分變換的紅光和黃光區(qū)域；Na+的敏感波段集中在949～1 355 nm近紅外區(qū)域；Ca2+敏感波段集中在665～672和919～1 283 nm可見光紅光及近紅外區(qū)域；Mg2+敏感波段主要集中在384、651～669 nm，主要為可見光紅光區(qū)域，與紫光區(qū)域波段也有一定相關(guān)性，但相關(guān)性整體偏小。原始光譜與K+、Na+的相關(guān)性較高，達(dá)到顯著水平，光譜變換使Ca2+、Mg2+含量與光譜的相關(guān)性增加，從而能夠按照標(biāo)準(zhǔn)選出建模波段，光譜變換對(duì)鹽離子含量與光譜相關(guān)性有提高作用。

3.3 二維植被指數(shù)特征波段選擇

二維相關(guān)系數(shù)圖能夠?qū)}離子含量和光譜指數(shù)之間的相關(guān)性進(jìn)行可視化表達(dá)。建立鹽生植物葉片4種鹽離子含量值分別與實(shí)測(cè)光譜反射率與RVI、DVI、NDVI的決定系數(shù)圖（圖5），依不同顏色選取相應(yīng)的敏感波段。

注：NDVI為歸一化型植被指數(shù)；DVI為差值型植被指數(shù)；RVI為比值型植被指數(shù)；下同。

表1 鹽離子建模波段及光譜變換形式

從圖5可以看出，構(gòu)建指數(shù)選取敏感波段時(shí)，同一種離子的不同指數(shù)構(gòu)建方法下二維相關(guān)圖具有相似性，即同一鹽離子含量在3種植物指數(shù)下的敏感波段集中區(qū)域相似。其中，Na+與3種指數(shù)的相關(guān)性整體大于其他3個(gè)離子，敏感波段主要集中在：1 480～1 800 nm，：500～750 nm分別位于近紅外與可見光黃、橙、紅及紅外區(qū)域；K+、Ca2+和Mg2+在選取敏感波段時(shí)，相關(guān)性較小，整體效果不明顯，K+在RVI、DVI集中在：1 650～1 800 nm，：1 650～1 800 nm的近紅外區(qū)域，NDVI集中在：650～700 nm，：650～700 nm的紅光區(qū)域；Ca2+在3種指數(shù)下均集中在：650～700 nm，：650～700 nm的紅光區(qū)域；Mg2+主要集中在：350～400、650～700 nm，：350～400、650～700 nm的可見光紫和紅光區(qū)域。

優(yōu)選出的用于建立基于二維指數(shù)的模型的敏感波段組合見表2，因?yàn)镃a2+在3種指數(shù)計(jì)算下選取的特征波段相同（684，664），為避免冗余，在3種指數(shù)計(jì)算中選取RVI指數(shù)建模。Mg2+在波段選取時(shí)，在RVI和NDVI 2種指數(shù)計(jì)算下選取的敏感波段相同，為（684，661）。后續(xù)建模時(shí)，只進(jìn)行1次建模。Na+與K+含量估算使用表中特征波段進(jìn)行建模。

3.4 GWR建模對(duì)比分析

一維相關(guān)性特征選取的策略和二維植被指數(shù)特征選取策略都具有較高的相關(guān)性，但對(duì)定量估算模型的影響不相同，為尋找更好的估算模型，通過各處理下的特征波段構(gòu)建2種鹽離子估算模型，得到如表3所示的結(jié)果。

表2 用于建模的二維指數(shù)的特征波段

表3 GWR模型與BP對(duì)比及精度檢驗(yàn)

注：建模樣本48個(gè)。驗(yàn)證樣本為16個(gè)。—，負(fù)值，不記錄數(shù)據(jù)。

Note：Samples for model establishment and validation are 48 and 16.—, negative value and not recorded.

圖6 基于建模集的不同方法預(yù)測(cè)值與實(shí)測(cè)值

圖7 基于GWR最優(yōu)模型預(yù)測(cè)值與實(shí)測(cè)值驗(yàn)證

圖8 基于BP最優(yōu)模型預(yù)測(cè)值與實(shí)測(cè)值驗(yàn)證

4 討 論

2種特征波段選取方案選取波段比較，選取的相關(guān)波段具有一定的重合性，K+特征波段在光譜變換下集中在紅光及黃光區(qū)域，在波段選取時(shí)集中在短波紅外及紅光區(qū)域，2種波段選取中重合部分為紅光區(qū)域；Na+敏感波段選取中，光譜變換下集中在短波紅外區(qū)域，指數(shù)選取下集中在紅外及可見光黃、橙、紅區(qū)域，重合部分為短波紅外區(qū)域；Ca2+在光譜變換下集中在紅光及短波紅外區(qū)域，在指數(shù)選取中集中在紅外區(qū)域，重合部分為紅光區(qū)域；Mg2+在光譜變換及指數(shù)選取下均集中在紅光及紫光區(qū)域。重合部分為紅光及短波紅外區(qū)域說明使用紅及短波紅外波段進(jìn)行離子估算的適用性及精度較好。

基于GWR模型對(duì)干旱區(qū)的鹽生植物鹽離子估算具有較優(yōu)的結(jié)果。GWR模型對(duì)有Na+估算效果最為顯著，這也是由于本研究區(qū)屬于干旱區(qū)的典型鹽漬化災(zāi)害區(qū)域，依據(jù)王勇輝等[33]的研究成果，本研究區(qū)鹽漬化土壤主要以中的陽離子含量以K+、Na+為主，植物主要吸收土壤中Na+為主的鹽分，而土壤鹽漬化存在隨機(jī)性和空間異質(zhì)性，所引起的鹽分脅迫也存在相似的屬性。故GWR模型對(duì)植物鹽離子含量的預(yù)測(cè)效果提高的作用大小取決于植物鹽離子與各變量間相關(guān)關(guān)系的空間非平穩(wěn)性程度[21]。本研究的鹽離子估算結(jié)果與Pandey等[34]的研究結(jié)果具有一定的相同點(diǎn)：K+估算精度令人滿意，Ca2+、Mg2+2種離子的估算結(jié)果具有潛在的高精度。BP神經(jīng)網(wǎng)絡(luò)在本文估算中表現(xiàn)并不好，或因機(jī)器學(xué)習(xí)需要大的樣本量進(jìn)行學(xué)習(xí)、訓(xùn)練才更易達(dá)到較好效果，而本文樣本數(shù)量較少，若加大樣本量，機(jī)器學(xué)習(xí)的精度應(yīng)所提高。依據(jù)本文所研究?jī)?nèi)容及樣本點(diǎn)設(shè)置，選取最優(yōu)模型時(shí)，不僅需具有較高的精度還需要考慮訓(xùn)練的成本。在區(qū)域范圍的小樣本條件下，GWR能夠取得高精度值得進(jìn)行推廣。

本文基于GWR方法對(duì)植物鹽離子含量得到較好的估算結(jié)果，尤其Na+預(yù)測(cè)最優(yōu)，說明地理數(shù)據(jù)的非平穩(wěn)性在其中起到關(guān)鍵作用。然而光譜數(shù)據(jù)僅僅提供地表參量的信號(hào)，從結(jié)果中可見一些光譜變換沒有起到挖掘信息的作用，仍然需要引入深度挖掘技術(shù)進(jìn)一步完善。此外，考慮到模型的局限性，易受到空間樣本量和離子含量的影響，可嘗試應(yīng)用遷移學(xué)習(xí)調(diào)試，使模型適用于不同季節(jié)、不同地區(qū)，以增加其普適性，進(jìn)一步驗(yàn)證，為GWR模型在植被鹽離子的光譜估算中提供更全面的角度。

5 結(jié) 論

通過對(duì)艾比湖保護(hù)區(qū)內(nèi)植物鹽離子的檢測(cè)，相較于其他3種鹽離子，Na+在鹽生植物中含量最高，K+次之，Ca2+、Mg2+含量相近。通過實(shí)測(cè)光譜14種波段變換及3種植物指數(shù)的構(gòu)建選取特征波段進(jìn)行建模，結(jié)果如下：

1）在建模過程中，Na+的含量原始光譜呈負(fù)相關(guān)，經(jīng)光譜變換后，選取特征波段集中在949～1 355 nm短波紅外區(qū)域，建模精度最好2均大于0.7。在構(gòu)建指數(shù)時(shí)，波段集中在1 480～1 800 nm，500～750 nm分別位于短波紅外與可見光黃、橙、紅區(qū)域。模型驗(yàn)證精度除差值型植被指數(shù)（difference vegetation index，DVI）外，擬合度均達(dá)到0.7以上。使用光譜數(shù)據(jù)估算Na+含量具有較高精度并有多種建模方法。

2）K+含量估算選取波段在400～700 nm的光合活躍區(qū)，且集中在不經(jīng)微分變換的紅光和黃光區(qū)域及在構(gòu)建植被指數(shù)后選取的近紅外區(qū)域。除通過反射率倒數(shù)一階、DVI、比值型植被指數(shù)（ratio vegetation index，RVI）方法選取波段外，建模效果較好。其中l(wèi)g(1/)選取特征波段效果最好。Ca2+、Mg2+含量較低，且與原始光譜的相關(guān)性較低，經(jīng)變換后的光譜能有效提高相關(guān)性，并在估算建模時(shí)Ca2+使用光譜變換選取特征波段建模效果較植被指數(shù)選取波段建模效果更優(yōu)，建模最優(yōu)為平方根的一階微分；Mg2+估算時(shí)，只有DVI選取的波段表現(xiàn)較好。

總體而言，在地理加權(quán)回歸（geographically weighted regression, GWR）模型下鹽離子模型估算精度整體比BP神經(jīng)網(wǎng)絡(luò)模型（back propagation neural network）高。GWR估算模型按精度評(píng)價(jià)排序?yàn)镵+、Na+＞Ca2+、Mg2+，對(duì)含量較高的離子估算效果更好。4種離子在不同波段選取方法下，最優(yōu)估算模型的特征波段均集中于紅及短波紅外波段?？煽紤]進(jìn)一步推廣紅及短波紅外波段對(duì)鹽生植物葉片鹽離子含量估算研究。

[1] 田華麗，夏藝，梁士楚，等. 桂林漓江濕地植被種類組成及其區(qū)系成分[J]. 濕地科學(xué)，2015，13(1)：103－110. Tian Huali, Xia Yi, Liang Shichu, et al. Species composition and flora of wetland vegetation in Lijiang River, Guilin[J]. Wetland Science, 2015, 13(1): 103－110. (in Chinese with English abstract)

[2] 韓濤，王大為. 2000—2014年石羊河流域植被覆蓋變化研究[J]. 中國(guó)農(nóng)學(xué)通報(bào)，2017，33(13)：66－74. Han Tao, Wang Dawei. Change of vegetation coverage in Shiyang River Basin 2000-2014[J]. Chinese Agricultural Science Bulletin, 2017, 33(13): 66－74. (in Chinese with English abstract)

[3] 趙可夫，范海，江行玉，等. 鹽生植物在鹽漬土壤改良中的作用[J]. 應(yīng)用與環(huán)境生物學(xué)報(bào)，2002，8(1)：31－35. Zhao Kefu, Fan Hai, Jiang Xingyu, et al. Improvement and utilization of saline soil by planting halophytes[J]. Chin J Appl Environ Biol, 2002, 8(1): 31－35. (in Chinese with English abstract)

[4] 趙丹，余冬梅，胡夏嵩，等. 青海柴達(dá)木盆地大柴旦鹽湖區(qū)鹽生植物降鹽效應(yīng)[J]. 鹽湖研究，2016，24(4)：37－47. Zhao Dan, Yü Dongmei, Hu Xiasong, et al. The bio- desalinization effect by halophytes in Da Qaidam salt lake area[J]. Journal of Salt Lake Research. 2016, 24(4): 37－47. (in Chinese with English abstract)

[5] 黃建，田長(zhǎng)彥，卞衛(wèi)國(guó)，等. 陳春秀.4種鹽生植物生長(zhǎng)對(duì)土壤石油污染的響應(yīng)[J]. 干旱區(qū)研究，2014，31(1)：100－104. Huang Jian, Tian Changyan, Bian Weiguo, et al. Response of growth of four halophyte species in oil-contaminated soil[J]. Arid Zone Research, 2014, 31(1): 100－104. (in Chinese with English abstract)

[6] 郗金標(biāo)，張福鎖，毛達(dá)如，等. 新疆鹽生植物群落物種多樣性及其分布[J]. 林業(yè)科學(xué)，2006，42(10)：6－12. Xi Jinbiao, Zhang Fusuo, Mao Daru, et al. Species diversity and distribution of halophytic vegetation in Xinjiang[J]. Scientia Silvae Sinicae, 2006, 42(10): 6－12. (in Chinese with English abstract)

[7] Yue J, Feng H, Jin X, et al. A comparison of crop parameters estimation using images from UAV-mounted snapshot hyperspectral sensor and high-definition digital camera[J]. Remote Sensing, 2018, 10(7): 1138.

[8] 劉煒，常慶瑞，郭曼，等. 夏玉米可見/近紅外光小波主成分提取與氮素含量神經(jīng)網(wǎng)絡(luò)檢測(cè)[J]. 紅外與毫米波學(xué)報(bào)，2011，30(1)：48－54. Liu Wei, Chang Qingrui, Guo Man, et al. Detection of leaf nitrogen content of summer corn using visible/near infrared spectra[J]. Journal of Infrared and Millimeter Waves, 2011, 30(1): 48－54. (in Chinese with English abstract)

[9] 劉克，趙文吉，郭逍宇，等. 基于地面實(shí)測(cè)光譜的濕地植物全氮含量估算研究[J]. 光譜學(xué)與光譜分析，2012，32(2)：465－471. Liu Ke, Zhao Wenji, Guo Xiaoyu, et al. Estimating Total nitrogen content in wetland vegetation based on measured reflectance spectra[J]. Spectroscopy and Spectral Analysis, 2012, 32(2): 465－471. (in Chinese with English abstract)

[10] Broge N H, Leblanc E. Comparing prediction power and stability of broadband and hyperspectral vegetation indices for estimation of green leaf area index and canopy chlorophyll density[J]. Remote Sensing of Environment, 2001, 76(2): 156－172.

[11] 易秋香，黃敬峰，王秀珍. 玉米粗纖維含量高光譜估算模型研究[J]. 紅外與毫米波報(bào)，2007，26(5)：393－395. Yi Qiuxiang, Huang Jingfeng, Wang Xiuzhen. Hyperspectral estimation models for crude fibre concentration of cone[J]. Journal of Infrared and Millimeter Waves, 2007, 26(5): 393－395. (in Chinese with English abstract)

[12] 劉冰峰，李軍，賀佳，等. 基于高光譜植被指數(shù)的夏玉米地上干物質(zhì)量估算模型研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)，2016，47(3)：254－262. Liu Bingfeng, Li Jun, He Jia, et al. Estimation models of above-ground dry matter accumulation of summer maize based on hyperspectral remote sensing vegetation indexes[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(3): 254－262. (in Chinese with English abstract)

[13] Wang Q, Jin J. Leaf transpiration of drought tolerant plant can be captured by hyperspectral reflectance using PLSR analysis[J]. iForest-Biogeosciences and Forestry, 2015, 9: e1－e8.

[14] Sampson P H, Zarcotejada P J, Mohammed G H, et al. Hyperspectral remote sensing of forest condition: estimating chlorophyll content in tolerant hardwoods.[J]. Forest Science, 2003, 49(3): 381－391.

[15] Lowe A, Harrison N, French A P. Hyperspectral image analysis techniques for the detection and classification of the early onset of plant disease and stress[J]. Plant Methods, 2017, 13(1): 80.

[16] Guo C, Guo X. Estimating leaf chlorophyll and nitrogen content of wetland emergent plants using hyperspectral data in the visible domain[J]. Spectroscopy Letters, 2016, 49(3): 180－187.

[17] 彭杰，向紅英，王家強(qiáng)，等. 葉面降塵的高光譜定量遙感模型[J]. 紅外與毫米波學(xué)報(bào)，2013，32(4)：313－319. Peng Jie, Xiang Hongying, Wang Jiaqiang, et al. Quantitative model of foliar dustfall content using Hyperspectral remote sensing[J]. Journal of Infrared and Millimeter Waves, 2013, 32(4): 313－319. (in Chinese with English abstract)

[18] 田喜，何紹蘭，呂強(qiáng)，等. 高光譜圖像信息的柑橘葉片光合色素含量分析技術(shù)研究[J]. 光譜學(xué)與光譜分析，2014，34(9)：2506－2512. Tian Xi, He Shaolan, Lü Qiang, et al. Determination of photosynthetic pigments in citrus leaves based on hyperspectral images data[J]. Spectroscopy and Spectral Analysis, 2014, 34(9): 2506－2512. (in Chinese with English abstract)

[19] 徐明星，吳紹華，周生路，等. 重金屬含量的高光譜建模反演:考古土壤中的應(yīng)用[J]. 紅外與毫米波學(xué)報(bào)，2011，30(2)：109－114.Xu Mingxing, Wu Shaohua, Zhou Shenglu, et al. Hyperspectral reflectance models for retrieving heavy metal content: application in the archaeological soil[J]. Journal of Infrared and Millimeter Waves, 2011, 30(2): 109－114. (in Chinese with English abstract)

[20] Fortheringham A S, Chanrlton M, Brunsdon C. The geographically of parameter space: An investigation of spatial nonstationarity[J]. International Journal of Geographical Information Systems, 1996, 10: 605－627.

[21] 江振藍(lán)，楊玉盛，沙晉明. GWR模型在土壤重金屬高光譜預(yù)測(cè)中的應(yīng)用[J]. 地理學(xué)報(bào)，2017，72(3)：533－544. Jiang Zhenlan, Yang Yusheng, Sha Jinming. Application of GWR model in hyperspectral prediction of soil heavy metals[J]. Acta Geographica Sinica, 2017, 72(3): 533－544. (in Chinese with English abstract)

[22] 王璐, 丁建麗. 艾比湖保護(hù)區(qū)荒漠植被時(shí)空過程變化及其植被指數(shù)影響因素分析[J]. 草業(yè)學(xué)報(bào)，2015，24(5)：4－11. Wang Lu, Ding Jianli. Vegetation index feature change and its influencing factor and spatial-temporal process analysis of desert grassland in the Ebinur Lake Nature Reserve, Xinjiang[J]. Acta Prataculturae Sinica, 2015, 24(5): 4－11. (in Chinese with English abstract)

[23] 龔雪偉，呂光輝，馬玉，等. 艾比湖流域2種典型荒漠鹽生植物冠下土與葉片的生態(tài)化學(xué)計(jì)量特征[J]. 林業(yè)科學(xué)，2017，53(4)：28－36. Gong Xuewei, Lü Guanghui, Ma Yu, et al. Ecological stoichiometry characteristics in the soil under crown an leaves of two desert halophytes with soil salinity gradients in Ebinur lake basin[J]. Scientia Silvae Sinicae, 2017, 53(4): 28－36. (in Chinese with English abstract)

[24] Wang X, Zhang F, Kung X, et al. New methods for improving the remote sensing estimation of soil organic matter content(SOMC) in the Ebinur Lake Wetland National Nature Reserve (ELWNNR)in northwest China[J]. Remote Sensing of Environment, 2018, 218(1): 104－118.

[25] 井云清，張飛，張?jiān)?，? 4個(gè)時(shí)期艾比湖濕地國(guó)家級(jí)自然保護(hù)區(qū)植被覆蓋度變化[J]. 濕地科學(xué)，2016，14(6)： 895－900. Jing Yunqing, Zhang Fei, Zhang Yue, et al. Changes of vegetation coverage in Aibihu wetland national nature reserve for 4 periods[J]. Wetland Science,2016, 14(6): 895－900. (in Chinese with English abstract)

[26] 張賢龍，張飛，張海威，等. 基于光譜變換的高光譜指數(shù)土壤鹽分反演模型優(yōu)選[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2018，34(1)：110－117.Zhang Xianlong, Zhang Fei, Zhang Haiwei, et al. Optimization of soil salt inversion model based on spectral transformation from hyperspectral index[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(1): 110－117. (in Chinese with English abstract)

[27] 李哲，張飛，馮海寬，等. 基于波段組合的植被葉片鹽離子估算研究[J]. 光學(xué)學(xué)報(bào)，2017，37(11)：317－331. Li Zhe, Zhang Fei, Feng Haikuan, et al. Research on the estimation of salt Ions of vegetation leaves based on band combination[J]. Acta Optica Sinica, 2017, 37(11): 317－331.(in Chinese with English abstract)

[28] 吳一全，周楊，盛東慧，等. 基于子空間中主成分最優(yōu)線性預(yù)測(cè)的高光譜波段選擇[J]. 紅外與毫米波學(xué)報(bào)，2018，37(1)：119－128. Wu Yiquan, Zhou Yang, Sheng Donghui, et al. Band selection of hyperspectral image based on optimal linear prediction of principal components in subspace[J]. Journal of Infrared and Millimeter Waves, 2018, 37(1): 119－128. (in Chinese with English abstract)

[29] Ge X, Wang J, Ding J, et al. Combining UAV-based hyperspectral imagery and machine learning algorithms for soil moisture content monitoring[J]. Peer J, 2019, 7: e6926.

[30] 江振藍(lán)，楊玉盛，沙晉明. 福州市土壤鉻含量高光譜預(yù)測(cè)的GWR模型研究[J]. 生態(tài)學(xué)報(bào)，2017，37(23)：533－544. Jiang Zhenlan, Yang Yusheng, Sha Jinming. Application of GWR model in hyperspectral prediction of soil heavy metals[J]. Acta Geographica Sinica, 2017, 37(23): 533－544. (in Chinese with English abstract)

[31] Mcmillen D P. Geographically weighted regression: The analysis of spatially varying relationships[J]. American Journal of Agricultural Economics, 2004, 86(2): 554－556.

[32] Steidle Neto A J, Lopes D C, Pinto F A C, et al. Vis/NIR spectroscopy and chemometrics for non-destructive estimationof water and chlorophyll status in sunflower leaves[J]. Biosystems Engineering, 2017, 155: 124－133.

[33] 王勇輝，馬蓓. 艾比湖主要補(bǔ)給河流下游河岸帶土壤鹽分特征[J]. 干旱區(qū)研究，2013，30(2)：196－202. Wang Yonghui, Ma Bei. Soil Salinity in the riparian zone in lower reaches of the main feeding rivers of the Ebinur lake[J]. Arid Zone Research, 2013, 30(2): 196－202. (in Chinese with English abstract)

[34] Pandey P, Ge Y, Stoerger V, et al. High throughput in vivo analysis of plant leaf chemical properties using hyperspectral imaging[J]. Frontiers in Plant Science, 2017, 8: 1348.

Leaf salt ion content estimation of halophyte plants based on geographically weighted regression model combined with hyperspectral data

Yuan Jie1,2, Zhang Fei1,2,3※, Ge Xiangyu1,2, Guo Wanzhen1,2, Deng Laifei1,2

(1.830046,; 2.830046,; 3.,,,830046,)

Rapid and non-destructive estimation of leaf salt ion concentrations in halophytes can provide valuable information for plant growth monitoring, selection of salt-tolerant plants and soil salinity monitoring. In this study, the canopy reflectance (350-2 500 nm) and the leaf salt ion (K+, Na+, Ca2+, Mg2+) concentration in the halophytes were measured in the Ebinur Lake Protection Zones, Xinjiang, China. Data collected includes hyperspectral data and leaf salt ion data, and the relationships between the leaf ion concentrations and the selected spectral indices were analyzed.K+sensitive wave bands on the photosynthetic effective radiation area of the 400- 700 nm (photosynthetically available radiation, PAR), and focused on the red and yellow areas without differential transform; The sensitive bands of Na+are concentrated in the near infrared region of 949- 1 355 nm. Ca2+sensitive bands were concentrated in the visible red and near-infrared regions of 665-672 and 919-1 283 nm. Mg2+sensitive bands were mainly concentrated in 384, 651- 669 nm, mainly in the visible red light region. There was a certain correlation with the ultraviolet region band, but the correlation was generally small. The correlation between the original spectrum and K+and Na+was relatively high, reaching a significant level. Spectral transformation increased the correlation between the contents of Ca2+and Mg2+and the spectrum, so that modeling bands could be selected according to the standard. Spectral transformation could improve the correlation between the content of salt ions and the spectrum.There were 64 samples in total, and the proportion of samples used for modeling and verification was 3:1.2and root mean squared error (RMSE) were used as accuracy evaluation criteria. A Geographically Weighted Regression (GWR) model and a back propagation (BP) model were constructed for estimating leaf salt ion concentrations with the spectral transform and the spectral indices as ratio vegetation index (RVI), difference vegetation index (DVI) and normalized difference vegetation index, and achieved a promising accuracy. The GWR estimation was the best in the bands in the red light region selected by the reciprocal logarithm of reciprocal of reflectance. The characteristic bands of Na+were concentrated in the short-wave infrared region under the spectral transformation, and the two-dimensional vegetation index was concentrated in the near-infrared region, short-wave near-infrared region, yellow, orange and red region. The short-wave infrared band selected under first order of square root for Ca2+content had the best estimation effect through GWR model. Mg2+content was best estimated in the characteristic bands in the red light region selected by DVI, but the GWR model was not as accurate as BP model in estimating Mg2+content. Based on the GWR salt ion model, the estimation of ions with higher content was better, and the accuracy of K+and Na+ models was better than that of Ca2+and Mg2+. When the GWR model was used to estimate the salt ion content in plant leaves, the characteristic bands all pointed to red and short-wave infrared bands. The model based on logarithms of reciprocal of reflectance and GWR for estimated K+produced the superior performance (2=0.930, RMSE=0.018 mg/kg). The optimal GWR model with the highest2and lowest RMSE was estimation model on Na+(2=0.984, RMSE=0.041 mg/kg) via processing. For the estimation model on Ca2+, the model produced reasonable outcome using first order of square root of reflectance-GWR strategy. Moreover, compared with BP model, the GWR model had insufficient estimation for Mg2+whereas DVI scheme contributed to improve accuracy of the BP estimated model. By comparison, the GWR model yielded better results in higher-content ion models. In conclusion, our study showed GWR model was effective for estimating leaf salt ions through vegetation spectral information. Sensitive bands for salt ions were prominent in the red bands and short-wave infrared bands, which were consistent with the response of vegetation spectral mechanism.

drought; leaf; hyperspectra; GWR model; halophyte; saline ions

10.11975/j.issn.1002-6819.2019.10.015

V232.4

A

1002-6819(2019)-10-0115-10

2019-01-30

2019-04-10

國(guó)家自然科學(xué)基金-新疆本地優(yōu)秀青年培養(yǎng)專項(xiàng)（U1503302）；新疆維吾爾自治區(qū)自然科學(xué)基金項(xiàng)目（2016D01C029）

袁 婕，主要從事干旱區(qū)植物遙感應(yīng)用研究。Email：yuanjie_0516@163.com

張 飛，博士，教授，主要從事干旱區(qū)資源與環(huán)境遙感應(yīng)用研究。Email：zhangfei3s@163.com

袁 婕，張 飛，葛翔宇，郭婉臻，鄧來飛. 地理加權(quán)回歸模型結(jié)合高光譜反演鹽生植物葉片鹽離子含量[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2019，35(10)：115－124. doi：10.11975/j.issn.1002-6819.2019.10.015 http://www.tcsae.org

Yuan Jie, Zhang Fei, Ge Xiangyu, Guo Wanzhen, Deng Laifei.Leaf salt ion content estimation of halophyte plants based on geographically weighted regression model combined with hyperspectral data[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(10): 115－124. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2019.10.015 http://www.tcsae.org