基于偏振反射模型和隨機森林回歸的葉片氮含量反演

張子晗,晏 磊,2,劉思遠(yuǎn),付 瑜,姜凱文,楊 彬,劉綏華,張飛舟*

1. 北京大學(xué)地球與空間科學(xué)學(xué)院， 遙感與地理信息系統(tǒng)研究所空間信息集成與3S工程應(yīng)用北京市重點實驗室,北京 100871 2. 桂林航空工業(yè)學(xué)院， 廣西高校無人機遙測重點實驗室,廣西 桂林 541004 3. 湖南大學(xué)電氣與信息工程學(xué)院,湖南 長沙 410082 4. 貴州師范大學(xué)地理與環(huán)境科學(xué)學(xué)院,貴州 貴陽 550001

引 言

在科研人員利用高光譜遙感數(shù)據(jù)進(jìn)行植被內(nèi)部重要物質(zhì)含量反演之后的30多年里， 高光譜遙感技術(shù)已廣泛應(yīng)用于研究與植被相關(guān)的多種重要的生物化學(xué)物質(zhì)， 例如： 水分、 纖維素、 蛋白質(zhì)、 葉綠素、 氮元素、 硫元素和磷元素等。 葉片氮元素含量(leaf nitrogen concentration,LNC)與植物光合呼吸作用和其他的相關(guān)生物過程密切相關(guān)， 其含量僅占葉片總質(zhì)量的0.2%～6.4%， 但卻能極大程度上影響全球氣候變化過程， 并在廣域碳氮循環(huán)過程中起到主導(dǎo)作用， 對于解構(gòu)生態(tài)環(huán)境規(guī)律有重大科學(xué)意義。 氮元素在葉片內(nèi)部的主要依托為蛋白質(zhì)和葉綠素， 這兩種物質(zhì)在可見光近紅外(400～2 500 nm)波段有獨特的吸收特性， 因而可以利用高光譜數(shù)據(jù)進(jìn)行氮元素含量反演。

利用高光譜數(shù)據(jù)進(jìn)行氮含量反演的算法可以大致分為三類： 第一類方法建立在葉片反射率與葉片氮含量之間統(tǒng)計關(guān)系的基礎(chǔ)上， 以逐次線性擬合回歸(stepwise multiple linear regression,SMLR)和偏最小二乘回歸(partial least squares regression,PLSR)為代表[1]， 在高光譜波段數(shù)較多時， 這類算法的回歸計算比較耗時， 精度相比于后兩類方法更低， 但由于算法原理簡單明晰， 仍被廣泛使用； 第二類算法建立在高光譜植被指數(shù)的基礎(chǔ)上， 通過計算高光譜植被指數(shù)來反演葉片氮含量[2-3]， 這類算法通過提取紅邊等深層次光譜特征來構(gòu)建植被指數(shù)， 對于高光譜數(shù)據(jù)的挖掘仍舊不夠深入， 但由于計算簡單， 常用于在農(nóng)學(xué)相關(guān)的應(yīng)用領(lǐng)域粗略估計葉片氮含量； 第三類算法建立在機器學(xué)習(xí)的基礎(chǔ)上， 以支持向量機(support vector machine,SVM)、 人工神經(jīng)網(wǎng)絡(luò)(artificial neural network,ANN)和隨機森林(random forest,RF)為代表[4-5]， 這類算法屬于遙感領(lǐng)域與計算機領(lǐng)域?qū)W科交叉的范疇， 具有極大的精度優(yōu)勢。 不同的算法選擇會產(chǎn)生算法固有誤差， 影響氮含量反演精度， 為了盡可能減小算法誤差， 同時盡可能利用高光譜數(shù)據(jù)蘊含的信息， 本研究中測試了多種氮含量反演方法， 最終以反演精度為指標(biāo)， 選擇了第三類算法中的RF算法作為本研究中葉片氮含量反演使用的主體框架。

上述三類氮含量反演算法的輸入?yún)?shù)皆為高光譜反射率， 但從定量遙感分析的角度來看， 光譜測量得到的反射率包括兩個部分的貢獻(xiàn)： 一部分來源于入射光在葉片內(nèi)部的多次散射過程， 另一部分來源于入射光在葉片表面的鏡面反射過程[6]。 在多次散射過程中， 入射光大部分透過葉片表面進(jìn)入葉片內(nèi)部， 與葉片內(nèi)部結(jié)構(gòu)中的葉綠素和蛋白質(zhì)等含氮生物化學(xué)物質(zhì)發(fā)生相互作用， 在葉片內(nèi)部多次散射后， 攜帶氮元素信息離開葉片， 這部分來自于葉片內(nèi)部多次散射的漫反射率所攜帶的光譜信息是氮含量反演密切關(guān)注的。 而在鏡面反射過程中， 入射光部分在葉表蠟質(zhì)層直接發(fā)生鏡面反射， 這部分反射能量不包含與葉片內(nèi)部葉綠素和蛋白質(zhì)等含氮生物化學(xué)物質(zhì)相關(guān)的信息， 在氮含量反演中屬于高光譜反射率自帶的數(shù)據(jù)誤差， 如果使用包含鏡面反射部分的高光譜反射率直接進(jìn)行氮含量反演會產(chǎn)生機理誤差， 影響氮含量反演精度。 依據(jù)菲涅爾原理， 葉片反射光的鏡面反射組分是部分偏振的， 而葉片反射光的漫反射組分由于在葉片內(nèi)部多次散射導(dǎo)致的消光作用是非偏振的。 因此， 通過偏振反射建模方法獲得雙向偏振分布函數(shù)(bidirectional polarization distribution function,BPDF)， 可模擬計算偏振反射率， 進(jìn)而扣除測量獲取的葉片高光譜反射率中與氮含量反演無關(guān)的鏡面反射組分， 從源頭減少機理誤差， 提高氮含量反演精度。

綜上所述， 影響葉片氮含量反演精度的誤差來源有兩個： 算法誤差和機理誤差。 選擇RF算法進(jìn)行葉片氮含量反演算法可以盡可能減小算法誤差， 構(gòu)建BPDF可以去除大部分機理誤差， 最終實現(xiàn)基于BPDF和RF的葉片氮含量精確反演。

1 實驗部分

1.1 氮含量反演數(shù)據(jù)

本研究中用于氮含量反演的數(shù)據(jù)為主數(shù)據(jù)集。 葉片氮元素反演的研究地點為德國巴伐利亞國家森林公園[7](49°3′19″N,13°12′9″E)。 根據(jù)優(yōu)勢樹種的類型， 選擇使用國際地圈生物圈計劃(international geosphere-biosphere programme,IGBP)土地覆蓋分類體系， 優(yōu)勢樹種與IGBP類別的對應(yīng)情況見表1。

表1 巴伐利亞國家森林公園內(nèi)優(yōu)勢樹種對應(yīng)的IGBP類別Table 1 IGBP classes of dominant vegetation inBavarian Forest National Park

由于本研究區(qū)域內(nèi)的優(yōu)勢樹種中闊葉樹木均為落葉闊葉品種， 針葉樹木都為常綠針葉品種， 所以在研究區(qū)域中被簡單標(biāo)注為3種類型： 落葉闊葉林(IGBP04)、 常綠針葉林(IGBP01)和針闊混交林(IGBP05)， 其點位空間分布如圖1所示。

圖1 葉片氮元素反演點位空間分布Fig.1 Spatial distribution of foliar nitrogen content retrieval plot

共有26個點位作為氮元素反演研究區(qū)(8片落葉闊葉林， 8片常綠針葉林， 10片針闊混交林)。 其中， 每個采樣區(qū)在實地占地為30 m×30 m， 點位范圍用分層隨機抽樣的方法選擇8～9棵樹木采樣， 單顆樹木至少采集20片葉片， 妥善保存好后在實驗室用凱式定氮法測量氮元素含量。 高光譜數(shù)據(jù)的獲取時間為2013年7月22日， 利用HySpex傳感器系統(tǒng)在3 000 m高度獲取高光譜影像(HySpex傳感器系統(tǒng)包括2臺高光譜成像光譜儀， 一臺覆蓋400～1 000 nm的光譜范圍， 另一臺覆蓋1 000～2 500 nm的光譜范圍)， 高光譜數(shù)據(jù)共有418個窄波段， 410～2 495 nm內(nèi)， 每5 nm一個通道。 氮含量反演使用的原始反射率是經(jīng)過幾何校正和大氣校正處理后的26個樣區(qū)內(nèi)的418個波段的雙向反射因子數(shù)據(jù)(bidirectional reflectance factor,BRF)， 除了BRF數(shù)據(jù)之外無人機系統(tǒng)還記錄了獲取數(shù)據(jù)時的觀測幾何(太陽天頂角、 觀測天頂角、 相對方位角)。

1.2 偏振反射建模數(shù)據(jù)

偏振反射建模數(shù)據(jù)為輔數(shù)據(jù)集。 本研究通偏振反射建模獲取BPDF， 進(jìn)而估計偏振反射率。 用于偏振反射建模的BRDF-BPDF數(shù)據(jù)庫[8]由搭載在PARASOL衛(wèi)星上的POLDER傳感器獲取。 本研究中為了建立合適的BPDF模型所使用的數(shù)據(jù)來自于POLDER經(jīng)過校正的產(chǎn)品， 包含了數(shù)據(jù)獲取時的經(jīng)緯度和觀測幾何、 地物的IGBP類別以及地物類型同質(zhì)性、 樣區(qū)的歸一化植被指數(shù)(normalized difference vegetation index,NDVI)、 氣溶膠情況、 6個波段(490,565,670,765,865和1 020 nm)的反射率和865 nm波段偏振反射率。

該BRDF-BPDF數(shù)據(jù)庫從POLDER的連續(xù)7年的數(shù)據(jù)中挑選出質(zhì)量最好的1年， 即2008年， 同時， 在POLDER數(shù)據(jù)中挑選了質(zhì)量最優(yōu)秀(盡可能滿足對全球范圍的覆蓋、 數(shù)據(jù)有效性高、 樣區(qū)IGBP類別勻質(zhì)性高)的點位。 最終數(shù)據(jù)庫的時間跨度和點位選擇定為： 包括2008年12個月份的數(shù)據(jù)， 每個月份采取50個均勻分布在全球的點位， 在每個點位上按照不同的觀測幾何選擇高質(zhì)量POLDER/PARASOL體系的數(shù)據(jù)。 本研究使用該數(shù)據(jù)庫中2008年7月份歐洲范圍內(nèi)的數(shù)據(jù)。

1.3 數(shù)據(jù)一致性說明

使用兩組相對獨立的數(shù)據(jù)集， 主數(shù)據(jù)集用于氮含量反演， 輔數(shù)據(jù)集用于建立偏振反射模型， 使用主數(shù)據(jù)集和輔數(shù)據(jù)集的流程圖如圖2所示。 從數(shù)據(jù)獲取的時間來看， 考慮到偏振反射建模具有較強的季節(jié)周期性特征， 為了與巴伐利亞國家森林公園的氮含量反演高光譜數(shù)據(jù)的獲取時間契合， 在進(jìn)行偏振反射建模時僅使用BRDF-BPDF數(shù)據(jù)庫數(shù)據(jù)中與機載數(shù)據(jù)同月份的數(shù)據(jù)， 以保證時間一致性。 從數(shù)據(jù)獲取的地理位置來看， 由于氮含量反演數(shù)據(jù)的研究位置在德國， 同時考慮到在BRDF-BPDF數(shù)據(jù)庫數(shù)據(jù)中獲取的點位數(shù)量不能太少， 最終根據(jù)BRDF-BPDF數(shù)據(jù)庫數(shù)據(jù)的經(jīng)緯度， 篩選出歐洲范圍內(nèi)的數(shù)據(jù)進(jìn)行偏振反射建模， 以保證地理位置一致性。 從兩組數(shù)據(jù)集的空間分辨率來看， POLDER數(shù)據(jù)的分辨率為6.5 km， 無人機高光譜數(shù)據(jù)的分辨率為1.65 m(400～1 000 nm)和3.3 m(1 000～2 500 nm)， 但由于BRDF-BPDF數(shù)據(jù)庫中收錄的POLDER數(shù)據(jù)勻質(zhì)性大于75%， 可認(rèn)為在數(shù)據(jù)匹配時受空間尺度的影響較小。

圖2 數(shù)據(jù)使用流程圖Fig.2 Flowchart of data utilization

1.3 方法

1.3.1 偏振反射建模

研究中獲取的高光譜數(shù)據(jù)為418個窄波段的BRF。 根據(jù)菲涅爾定律， 地表反射是部分偏振的， 其中的偏振部分主要為線偏振， 也是由觀測幾何(太陽天頂角、 觀測天頂角、 相對方位角)決定的， 可用BPDF表征其多角度偏振反射特性。 有文獻(xiàn)表明， 在可見光近紅外波段， BPDF可視為觀測幾何穩(wěn)定時的光譜不變量， 在不同波段的偏振反射率為定值[9]。 BPDF的決定因素是地物類型和觀測幾何， 為了去除地表反射率的偏振部分， 需要利用BRDF-BPDF數(shù)據(jù)庫中的7月份歐洲范圍的數(shù)據(jù)， 建立3個分別適應(yīng)IGBP01， IGBP04和IGBP05這3種不同地表IGBP類別的BPDF模型以模擬估計偏振反射率(polarization bidirectional reflectance factor,PBRF)， 進(jìn)而從BRF中去扣除具有部分偏振特性的鏡面反射組分——偏振反射率PBRF， 最終獲得漫反射反射率(diffuse bidirectional reflectance factor,DBRF)， 以消除機理誤差的影響。

目前精度較高的BPDF模型幾乎全部為半經(jīng)驗?zāi)Ｐ停?純粹的經(jīng)驗?zāi)Ｐ碗m然能夠在局部較好擬合偏振反射率， 但由于缺乏物理機理， 模型參數(shù)嚴(yán)重依賴訓(xùn)練數(shù)據(jù)， 且過擬合明顯。 半經(jīng)驗?zāi)Ｐ屯ㄟ^較少的經(jīng)驗參數(shù)對包含物理意義的模型進(jìn)行校正調(diào)整， 其結(jié)果是更具有普適性。 本研究選擇了5個應(yīng)用最廣、 精度最高的BPDF模型用于模擬IGBP01， IGBP04和IGBP05這3種不同地表IGBP類別下的偏振反射率PBRF， 分別為： Nadal模型[10]、 Waquet模型[11]、 Maignan模型[12]、 Litvinov模型[13]和Diner模型[14]。

為了使公式簡潔明了， 在BPDF模型的公式表達(dá)中使用一些簡記。 其中， 太陽光方向與天頂方向(其夾角稱為太陽天頂角θS， 觀測方向與天頂方向的夾角稱為觀測天頂角θV， 太陽光方向在地表的投影與觀測方向在地表的投影之間的夾角稱為相對方位角φ， 鏡面反射過程中的法線與天頂方向的夾角稱為半天頂角θh， 入射方向與反射方向之間的夾角的一半為入射角αI， 折射角為αT， 入射方向與反射方向之間方向的變化為方向散射角γ， 折射率為N。 部分角度的余弦值簡記為

μS=cosθS

(1)

μV=cosθV

(2)

μI=cosαI

(3)

μT=cosαT

(4)

由菲涅爾公式Fp(γ,N)得到的變量之間的相關(guān)關(guān)系見式(5)

5種模型的偏振反射建模公式表達(dá)見表2。

表2 5種偏振反射模型的公式表達(dá)Table 2 Formulas of 5 BPDF models

利用輔數(shù)據(jù)BRDF-BPDF數(shù)據(jù)庫中的數(shù)據(jù)分別建立IGBP01，IGBP04和IGBP05這3種地表覆蓋類型下的5種BPDF模型， 分別擬合得到自由參數(shù)。 通過十折交叉驗證的數(shù)據(jù)分組策略， 以均方根誤差(rootmeansquareerror,RMSE)和可決系數(shù)(coefficientofdetermination,RSQ)為指標(biāo)， 挑選3種地表覆蓋類型下的最優(yōu)BPDF模型來模擬偏振反射率PBRF。

1.3.2 氮含量反演

氮含量反演研究基于德國巴伐利亞國家森林公園的實測氮含量主數(shù)據(jù)、 無人機高光譜主數(shù)據(jù)和偏振反射建模獲得模擬仿真結(jié)果， 研究主要包括三個步驟： 首先， 通過BPDF模擬計算偏振反射率PBRF， 進(jìn)而計算漫反射率DBRF， 獲取氮含量反演的算法輸入對照組(原始反射率BRF和去偏漫反射率DBRF)； 然后， 設(shè)置RF算法參數(shù)， 利用整體散射系數(shù)和去偏后的漫反射散射系數(shù)作為算法輸入對照組進(jìn)行葉片氮含量反演； 最后， 橫向?qū)Ρ萊F算法和其他算法的表現(xiàn)， 并對偏振反射率在氮含量反演中的影響進(jìn)行定量評估。

1.3.3 精度評估

偏振反射建模使用POLDER/PARASOL體系的數(shù)據(jù)庫中數(shù)據(jù)量較大， 采用十折交叉驗證的思路進(jìn)行分組后驗證精度， 重復(fù)10次十折交叉驗證， 取均方根誤差RMSE和可決系數(shù)RSQ的均值作為BPDF模型精度的度量。 氮含量反演使用26個點位的高光譜數(shù)據(jù)受到數(shù)據(jù)量的限制， 采用留一法交叉驗證的思路進(jìn)行數(shù)據(jù)分組后驗證精度， 同樣以均方根誤差RMSE和可決系數(shù)RSQ作為葉片氮含量反演精度的度量， 見式(6)和式(7)。

2 結(jié)果與討論

2.1 BPDF模型選擇結(jié)果與分析

首先， 按照氮含量反演研究時間、 區(qū)域進(jìn)行輔數(shù)據(jù)集BRDF-BPDF數(shù)據(jù)庫的篩選， 選擇了POLDER數(shù)據(jù)庫中2008年7月采集的歐洲范圍內(nèi)的數(shù)據(jù)； 然后， 利用十折交叉驗證的數(shù)據(jù)劃分思路將數(shù)據(jù)庫中的數(shù)據(jù)隨機劃分為訓(xùn)練集和驗證集， 利用訓(xùn)練集數(shù)據(jù)擬合獲得5個BPDF模型的自由參數(shù)； 獲得參數(shù)組后利用中位數(shù)方法求得單樣區(qū)最優(yōu)模型參數(shù)組， 根據(jù)該模型參數(shù)組預(yù)測偏振反射率PBRF； 將預(yù)測得到的偏振反射率與包含于數(shù)據(jù)庫中的實測865 nm偏振反射率數(shù)據(jù)進(jìn)行對比， 最后利用均方根誤差RMSE和可決系數(shù)RSQ來評價BPDF模型的效果， 5種模型的評價結(jié)果見表3。

表3 IGBP01/04/05地表覆蓋類型下5種BPDF模型精度評價Table 3 Accuracy assessment of 5 BPDF models on IGBP01/04/05

最終， 以可決系數(shù)RSQ為主要評判指標(biāo)， 均方根誤差RMSE為輔助評判指標(biāo)， 在RSQ盡可能大的情況下選擇RMSE盡量小的BPDF模型。 從整體來看， 在同一地表覆蓋類型下不同BPDF模型的模擬精度差異不大(保留3位小數(shù)后， 均方根誤差RMSE差異極小， 可決系數(shù)RSQ差異相對明顯)， 而在不同地表覆蓋類型下， IGBP04地表覆蓋類型的BPDF模型表現(xiàn)顯著優(yōu)于IGBP01和IGBP05， 說明闊葉樹木的偏振效應(yīng)相較于針葉樹木而言更易模擬， 推測可能與起偏介質(zhì)的面積大小有關(guān)。 模型選擇結(jié)果為： 在IGBP01地表覆蓋類型下最優(yōu)BPDF模型為Nadal模型， 在IGBP04地表覆蓋類型下最優(yōu)BPDF模型為Litvinov模型， 在IGBP05地表覆蓋類型下最優(yōu)BPDF模型為Maignan模型。 由于不同地表覆蓋類型的起偏特性不同， 最優(yōu)BPDF模型并不統(tǒng)一， 分別利用選出的3種BPDF模型模擬對應(yīng)的3種地表覆蓋類型的偏振反射率。

2.2 RF模型參數(shù)設(shè)置與葉片氮反演精度橫向?qū)Ρ确治?/h3>

在偏振反射建模的基礎(chǔ)上利用原始反射率BRF和漫反射率DBRF作為RF算法輸入進(jìn)行葉片氮含量反演。 RF的核心參數(shù)為決策樹數(shù)量， 本研究分析比對了去偏前后決策樹數(shù)量為5， 10， 50和100的RF算法的表現(xiàn)， 見表4。 研究結(jié)果說明， 在決策樹數(shù)量一定時， 使用去偏數(shù)據(jù)進(jìn)行葉片氮含量反演普遍能獲得精度提升且在RF決策樹數(shù)量設(shè)定為10時獲得最優(yōu)氮含量預(yù)測模型。 利用RF算法進(jìn)行氮含量反演能獲取的最高精度為： 可決系數(shù)達(dá)到0.803， 均方根誤差達(dá)到0.252。

表4 不同決策樹數(shù)量的RF算法精度評價Table 4 RF algorithm accuracy assessment ofdifferent estimator numbers

2.3 偏振反射率占比統(tǒng)計分析

在26個采樣點位， 偏振反射率PBRF在整體反射率BRF中的平均占比為3.580%～11.649%， 各個點位內(nèi)的偏振反射率占比箱形圖見圖3， 箱形圖形象地展示了偏振反射率PBRF占比的算數(shù)均值、 最小值、 下四分位數(shù)、 中位數(shù)、 上四分位數(shù)和最大值。 該分布表明偏振反射率PBRF在整體反射率BRF中平均占比為3.350%， 最高占比16.519%， 說明偏振反射率的影響是不可忽略的， 機理誤差足以影響反演精度， 在高光譜葉片氮含量反演中有去除偏振反射率PBRF的必要性。

圖3 偏振反射率在整體反射率中占比的分布圖中箱體內(nèi)的×表示算數(shù)均值， 箱體內(nèi)的黑色橫線表示中位數(shù)Fig.3 Ratio of PBRF to BRFIn this figure,the × in the box represents the mean value,and the transverse line in the box represents the median value

2.4 葉片氮反演精度橫向?qū)Ρ确治?/h3>

橫向?qū)Ρ攘薖LSR算法、 主成分回歸(principal component regression,PCR)算法、 支持向量回歸(support vector regression,SVR)算法、 K-近鄰(K-nearest neighbor,KNN)回歸算法和RF回歸算法在氮含量反演中的表現(xiàn)， 利用留一法進(jìn)行交叉驗證， 精度評價結(jié)果見表5。 其中， 由于偏最小二乘算法應(yīng)用最為廣泛[1]， 將其作為基線算法進(jìn)行整體精度評估。 研究結(jié)果說明， 在利用上述方法進(jìn)行氮含量反演的過程中， 去除偏振后葉片氮含量反演的精度普遍都有提升(KNN算法在選擇的最鄰近元數(shù)目K足夠大時， 最鄰近元的選擇在去除偏振前后相同， 因而模型預(yù)測結(jié)果相同， 無精度提升)。 以均方根誤差RMSE為核心指標(biāo)評估提升幅度， 平均提升幅度為4.244%， 其中RF算法去除偏振后葉片氮含量反演精度提升幅度最大， 為13.103%， 且精度最優(yōu)(達(dá)到了上述5種反演算法中的最高可決系數(shù)0.803和最低均方根誤差0.252)， 充分說明了RF算法利用高光譜數(shù)據(jù)進(jìn)行葉片氮含量反演的潛力。 相較于作為基線算法的PLSR(未進(jìn)行光譜去偏)， 通過偏振反射建模去除大部分光譜機理誤差并通過選擇隨機森林模型減小算法誤差后， 葉片氮含量反演精度整體提高32.440%， 證明了去除機理誤差和減少算法誤差的必要性， 實現(xiàn)了高精度葉片氮含量反演。

表5 5種葉片氮含量反演算法精度評價Table 5 Accuracy assessment of 5 LNC retrieval algorithms

3 結(jié) 論

研究旨在提高葉片氮含量反演精度， 從機理誤差和算法誤差兩個誤差源入手， 通過偏振反射率建模去除了大部分機理誤差， 通過選擇精度最優(yōu)的隨機森林算法并調(diào)整參數(shù)盡可能降低了算法誤差， 最終實現(xiàn)了葉片氮含量精確反演， 相比于基線方法整體精度提升為32.440%(可決系數(shù)RSQ達(dá)到0.803， 均方根誤差RMSE達(dá)到0.252)。 同時， 在討論中定量評估了偏振反射率PBRF在整體反射率BRF中的占比分布并對比了去除偏振前后多種氮含量反演算法的精度變化， 體現(xiàn)了隨機森林算法在高光譜信息挖掘中的優(yōu)勢， 也證明了去除不攜帶氮含量信息的偏振反射率的必要性。

致謝：本研究使用的德國巴伐利亞國家森林公園葉片氮含量實測數(shù)據(jù)和無人機高光譜數(shù)據(jù)由威斯康星大學(xué)麥迪遜分校的王智慧博士提供， 在此表示衷心的感謝！