紅邊位置改進算法的冬小麥葉綠素含量反演

錢彬祥，黃文江，葉回春，孔維平，任　淯，邢乃琛，焦全軍

·農業信息與電氣技術·

紅邊位置改進算法的冬小麥葉綠素含量反演

錢彬祥1,2,3，黃文江2,3※，葉回春2,3，孔維平4，任淯1,3，邢乃琛1,3，焦全軍2,3

（1. 中國科學院大學，北京 100049；2. 海南省地球觀測重點實驗室，三亞 572029；3. 中國科學院空天信息創新研究院，數字地球重點實驗室，北京 100094； 4. 中國科學院空天信息創新研究院，定量遙感信息技術重點實驗室，北京 100094）

植被反射光譜的紅邊位置對葉綠素含量高度敏感，利用遙感數據建立基于紅邊位置的作物葉綠素含量反演模型，可實現大范圍作物及時的長勢監測。該研究以冬小麥為研究對象，在學習6種經典紅邊位置求解算法的基礎上，提出牛頓-切比雪夫插值法和牛頓八點插值法2種改進紅邊位置求解算法。根據不同算法的紅邊位置分布特征綜合分析了改進算法的優缺點，并在此基礎上建立基于紅邊位置的冬小麥葉綠素含量反演模型。結果表明，與傳統算法相比，2種改進算法均顯著改善了雙峰現象和紅邊位移，且基于改進算法的模型預測值與葉綠素含量實測值的決定系數>0.619，較最大一階導數法，牛頓八點插值法提高了6.321%～9.947%，牛頓-切比雪夫插值法提高了5.024%～10.480%，具有更高的精度。同時，在2種改進算法中，牛頓八點插值法具有更高的穩定性與實用性。研究結果為植被理化參數反演與農業生產應用提供理論與技術支撐。

算法；遙感；模型；冬小麥；葉綠素反演；牛頓-切比雪夫插值法；牛頓八點插值法；紅邊位置

0 引 言

葉綠素含量是一種綜合反映作物受外界環境脅迫情況、光合作用能力強弱和新陳代謝旺盛程度的非常敏感的指示劑[1-2]。研究表明，作物反射光譜的紅邊位置（Red Edge Position，REP）與葉綠素含量具有相當高的相關關系，是作物生理生化參數遙感反演的一個重要指標[3-5]。建立基于紅邊位置的作物葉綠素含量反演模型，為實現大面積作物及時的長勢監測提供了一種快速而實用的方法[6]。紅邊波段一般認為在可見光波段的660～770 nm范圍內，由于帶寬較寬，不利于其與作物理化參數相關模型的建立，所以通常用位于紅光-近紅外波段作物反射光譜的突變點來精確定位紅邊位置[7-8]。

目前，紅邊位置求解算法可分為基于微分插值和基于曲線擬合兩大類。基于曲線擬合思想的算法提出的更早，Miller等[3]提出倒高斯模型法（Red Edge Position calculated by Inverted Gaussian，REP_IG），即用一條傾斜的倒高斯曲線近似紅邊波段范圍內的作物反射光譜，該模型較為復雜，需先經等價變換才能進行下一步解算。為了簡化計算，Demetriades-Shah等[9]用紅邊波段范圍內作物反射光譜的一階導數最大值對應的波段位置確定紅邊位置。Dawson等[10]研究表明最大一階導數法（Red Edge Position calculated by Maximum First Derivative，REP_MFD）提取紅邊位置的精度取決于傳感器平臺的波段范圍和光譜分辨率，為減少紅邊位置求解精度對光譜分辨率的過度依賴，提出了拉格朗日三點內插法（Red Edge Position calculated by Lagrange，REP_LAGR）。Clevers等[11]認為拉格朗日三點內插法只適用于粗分辨率光譜，且不能解決紅邊位置求解過程中存在雙峰現象（bimodal phenomenon，BP）的問題，為此提出了線性四點內插法（Red Edge Position calculated by Linear Four-Point Insert，REP_LFPI），該方法是最簡單、易行的紅邊位置求解算法，只需要4個波段和簡單的線性插值運算。此外，Pu等[12]采用5次多項式擬合作物反射光譜，多項式一階差分最大值對應的波段位置即為紅邊位置，并利用實測光譜數據驗證了該方法的合理性。Lin等[13]開展了拉格朗日三點內插法和多項式擬合法（Red Edge Position calculated by Polynomial fitting，REP_POLY）的精度對比研究，發現拉格朗日三點內插法的精度取決于波段組合，且對噪音尤為敏感，而多項式擬合法的適用性更高、穩定性更強。Cho等[14]為減少雙峰現象帶來的不穩定性的影響，提出了線性外推法（Red Edge Position calculated by Line Extrapolate，REP_LE）。陳西亮等[4]將牛頓插值法（Red Edge Position calculated by Newton Interpolation，REP_NI）應用于紅邊位置求解過程，利用6個等距的紅邊波段反射光譜數據內插出紅邊位置，求解精度位于倒高斯模型法和線性四點內插法之間。陳西亮等[4]的研究中存在以下不足：1）選取ENVI自帶的USGS植被波譜庫中的矮松樹、山楊樹、草坪和藍云杉4種植被的4條反射光譜做分析，數據量過少，分析結果受偶然誤差和樣本選取的影響較大，且與實際生產相脫離，研究缺乏普適性；2）由于缺乏試驗數據，無法分析不同算法求解出的紅邊位置的分布特征（如雙峰現象、紅邊位移），僅利用紅邊位置的相對誤差來衡量算法的優劣缺乏科學性；3）沒有將新算法與作物生理生化參數遙感反演相結合，無法檢驗新算法的實用性。

與傳統算法相比，牛頓插值法具有明顯優勢。傳統算法普遍對傳感器提出光譜分辨率高、采樣均勻的要求，且求解的紅邊位置存在雙峰現象、紅邊位移（Displacement of Red Edge Position，REPD）和精度不足等缺點。牛頓插值法屬于微分插值方法，相較于曲線擬合方法需要的數據量更少，并且能夠完全保留原始數據，以及引入差分（Difference，DF）和差商（Difference Quotient，DFQ）的概念，可以計算非等距節點的光譜反射率。同時，與拉格朗日插值法相比較，牛頓插值法不僅克服了每增加一個節點整個計算必須重新開始的缺點，而且可以節省乘、除法運算的次數，提高了計算速度[15]。

鑒于牛頓插值法的優越性以及前人研究工作還有待進一步完善，本研究以冬小麥為試驗對象，采用數值分析技術，在學習經典算法的基礎上提出2種改進的紅邊位置求解算法：牛頓八點插值法（Red Edge Position calculated by Newton Eight-Point Interpolation，REP_NEPI）和牛頓-切比雪夫插值法（Red Edge Position calculated by Newton-Chebyshev-Node Interpolation，REP_NCNI），并進行兩者的綜合比較，擬建立更高精度的基于紅邊位置的冬小麥葉綠素含量反演模型。研究結果可以為植被理化參數遙感反演增添新的理論依據，為農業工程領域的冬小麥葉綠素含量監測、長勢評估等實際生產應用提供新的技術支撐。

1 材料與方法

1.1 研究區概況及試驗設計

研究區位于北京市昌平區小湯山國家精準農業示基地（40°11′N，116°26′E），氣候類型為暖溫帶半濕潤大陸季風氣候。本研究選取2002—2003年和2018—2019年小麥生長季內不同品種、不同肥水處理下的冬小麥連續開展了野外測量試驗，在其關鍵生育期采集常規農學數據與近地高光譜數據。具體處理方式如下：

1）小麥品種：2002—2003年選用3種不同株型、籽粒品質的冬小麥，包括株型緊湊、籽粒品質弱筋型的“京411”，株型松散、籽粒面筋質量好的“中優9507”和株型與籽粒品質中間型的“京冬8”[7-8]。2018—2019年選用株型半緊湊、籽粒面筋質量好的“輪選167”和株型緊湊、籽粒品質中間型的“京冬18”。

2）肥水處理：2002—2003年的研究區共有48個肥水處理小區，各小區面積為32.4 m×30 m，在試驗區內對不同品種的冬小麥分別進行肥水脅迫處理。4個梯度施肥處理：0N處理，返青期至成熟期不施肥；②1N處理，返青期和拔節期各追施尿素50 kg/hm2；8N處理，返青期和拔節期各追施尿素200 kg/hm2；④14N處理，返青期和拔節期各追施尿素350 kg/hm2[7-8]。4個梯度供水處理：0水處理，全生育期不灌溉；②1水處理，全生育期灌水225 m3/hm2；2水處理，全生育期灌水450 m3/hm2；3水處理，全生育期灌水675 m3/hm2[7-8]。2018—2019年的研究區共有32個肥水處理小區，各小區面積為10 m×15 m，在試驗區內正常供水，對不同品種的冬小麥進行肥脅迫處理：N0處理，返青期至成熟期不施肥；N1處理，基肥（2018年9月）97.5 kg/hm2和拔節期（2019年4月）追肥97.5 kg/hm2；N2處理，基肥（2018年9月）195 kg/hm2和拔節期（2019年4月）追肥195 kg/hm2；N3處理，基肥（2018年9月）292.5 kg/hm2和拔節期（2019年4月）追肥292.5 kg/hm2。

3）測定項目：不同生育時期的冠層光譜反射率、配套的數碼照片和包括葉綠素含量（mg/g）、葉面積指數在內的常規田間調查農學參數。

4）測定時期：包括返青期（3月25日前后）、拔節期（4月15日前后）、挑旗期（4月25日前后）、開花期（5月10日前后）、灌漿中期（5月20日前后）、灌漿后期（5月30日前后）在內的冬小麥6個關鍵生育期。通過詢問試驗區內經驗豐富的技術人員，獲取試驗區內該年度冬小麥的長勢情況、降水量和氣溫變化等信息，并就此制定具體的測量時間與試驗方案。

1.2 光譜和農學參數測定

小麥冠層光譜采用美國ASD公司生產的FieldSpec FR光譜輻射儀測定[16]，選取晴朗無云、風力較小的10:30－14:00（北京時間）進行，采樣范圍為350～2 500 nm，采樣間隔為1 nm。測定時，探頭在距離冠層約1 m處垂直向下觀測，測定前后用標準參考板校正，每個點重復測定20次，取全部測數的均值作為該點的冠層光譜反射率。同時，在目標區域附近利用打孔器獲取測試樣品，稱量后在黑暗環境下浸泡于95%乙醇溶液中24～48 h，使用UV765PC型分光光度計測定溶液在可見光440、649和665 nm處的吸光度值，并由此計算葉綠素含量（mg/g）。

1.3 紅邊位置求解算法

1.3.1 傳統的紅邊位置求解算法

傳統的紅邊位置求解算法主要有以下6種：

1）最大一階導數法

在紅邊波段范圍內，植被光譜曲線反射率的一階差分最大值對應的波長位置即為紅邊位置[17]，其計算如式（1）和式（2）所示：

2）倒高斯模型法

采用一條傾斜的倒高斯曲線近似紅邊范圍內的植被反射光譜，倒高斯曲線函數表達式[18]如式（3）和式（4）所示：

式中()為倒高斯函數；R為植被光譜反射率在紅邊波段范圍內的最大值；0、0分別為紅邊位置范圍內植被光譜反射率最小值及其所對應的波長位置；為高斯模型標準差系數，nm。

3）拉格朗日三點內插法

拉格朗日三點內插法是一種適用于粗采樣光譜的三點插值算法，具體原理是利用二次多項式擬合植被反射率光譜一階導數[16]，如式（5）和式（6）所示：

4）線性四點內插法

其假設紅邊波段范圍內的植被反射率光譜近似為一條直線，紅邊位置由4個特征點的反射率值（670、700、740和780）內插得到[19]，如式（7）和式（8）所示：

式中rep為拐點反射率。

5）多項式擬合法

多項式擬合法是利用高次多項式高精度逼近紅邊波段的植被光譜曲線，經多次試驗表明9階多項式的擬合效果最好，均方根誤差<0.15和決定系數>0.999 93，其計算如式（9）和式（10）所示：

6）線性外推法

基于一階導數光譜在遠紅色（680～700 nm）和近紅外（725～760 nm）區域的2條直線線性外推紅邊位置[20]，如式（11）所示：

式中1、1分別為遠紅色波段直線的截距與斜率；2、2分別為近紅外波段直線的截距與斜率。

1.3.2 本研究改進的紅邊位置求解算法

在學習傳統方法的基礎上，本研究運用數值分析技術將牛頓插值法應用到紅邊位置的求解過程中。在利用插值多項式進行近似逼近時，通常認為構建多項式的原始插值節點越多、多項式階數越高，多項式逼近原始數據的效果就越好，可事實并非如此。龍格在20世紀初就證明了高次插值多項式具有病態性：基于等距節點的高階插值多項式插值點數量越多，多項式在逼近區間兩端產生的振蕩現象就會越明顯，即插值結果越偏離原函數[4]，這種病態性現象被稱為龍格現象[21]。龍格現象在牛頓插值法中尤為明顯，為了有效減少龍格現象的影響，本研究提出了以下2種解決方法：一是選擇適當的插值次數；二是借助切比雪夫零點替換等距節點。

1）牛頓八點插值法

為確定適當的插值次數，本研究進行多次試驗得出：在紅邊位置范圍內采集跨度為20 nm左右的8個插值點的插值效果最好，8個插值點分別位于651、671、691、711、731、751、771和790 nm處。利用差商原理將8個插值點代入牛頓插值多項式求取各階系數[4]如式（12）所示：

2）牛頓-切比雪夫插值法

牛頓八點插值法中的節點選取比較固定，會降低算法的實用性，將切比雪夫零點應用到插值節點選取中，用零點替換等距節點，可以使得插值點的選取更具有科學性和普適性[22]，其計算如式（13）和式（14）所示：

式中T()為切比雪夫多項式；、分別為零點區間的上下限；x為零點解。

進行上述9種算法核心思想的對比，如表1所示。

表1 8種紅邊位置求解算法的特征描述

1.4 葉綠素含量反演方法與精度評價指標

研究采用紅邊位置分布特征的描述統計量對不同算法進行特性分析，這些統計量包括最值、平均值、相對誤差（Relative Error，RE）等。并運用最小二乘回歸原理（Least Square Regression，LSR）建立基于不同算法紅邊位置的冬小麥葉綠素含量反演模型，對反演模型的預測值與葉綠素含量實測值進行線性擬合，評定各個模型的精度，精度評價指標包括擬合方程的斜率、截距、決定系數（coefficient of determination，2）、均方根誤差（Root Mean Squared Error，RMSE）和標準均方根誤差（Normalized Root Mean Squared Error，NRMSE）。各指標計算如式（15）～式（18）所示：

式中x為紅邊位置的相對誤差；為各類算法所求的紅邊位置，nm；此處，、y分別為植被光譜反射率實測值的數組和元素；、ye分別為植被光譜反射率擬合值的數組和元素；?為植被光譜反射率實測值的平均值；和分別表示實測數據的編號和總數。一般認為，當NRMSE≤0.1時，認為模型的擬合優度很好；當0.10.3時，認為模型的擬合優度較差[23-25]。

挑旗期、開花期、灌漿中期是冬小麥生長周期中新陳代謝最為旺盛的3個時期，本研究選取以上3個生育期的原始試驗數據進行研究，其中，2002－2003年的樣本數據（樣本數為136）用于不同算法的特性分析與冬小麥葉綠素含量反演模型的建立，2018－2019年的樣本數據（樣本數為64）進行模型的檢驗。

2 結果與分析

2.1 不同算法的特性分析

實現8種紅邊位置求解算法，并對建模集數據進行處理，得到不同算法紅邊位置的分布特征和算法特性（表2）。同時，為了直觀地觀察不同算法的紅邊位移與雙峰現象，繪制出的紅邊位置與葉綠素含量的散點圖（圖1）。各種算法紅邊位置的最小值在694～717 nm，主要集中在694～706 nm，線性四點內插法的最小值顯著偏大（717 nm）；紅邊位置的最大值在729～736 nm，主要集中在735～736 nm，倒高斯模型法的最大值相對偏小（729 nm）；紅邊位置的平均值在719.5～727.2 nm，主要集中在725.4～727.2 nm，倒高斯模型法和線性外推法的最大值明顯偏小，分別為715.9和721.9 nm。紅邊位置的變幅是評價紅邊算法對植被反射光譜敏感程度的重要指標之一，在一定程度上，紅邊位置的變幅越大，說明該算法對植被光譜曲線越敏感；變幅較小會存在紅邊位置的飽和現象。以上算法紅邊位置的變幅在14～41 nm，主要集中在30～41 nm，線性四點內插法的變幅明顯偏小（14 nm）。考慮到最大一階導數法是嚴格遵守紅邊位置定義的算法，所以將其計算結果作為參考標準計算得到其他算法的RE在0～0.882%，主要集中在0～0.179%。其中，拉格朗日三點內插法的RE最小，可忽略不計，倒高斯模型法和線性外推法的RE較大，分別為0.882%和0.551%。

雙峰現象和紅邊位移對葉綠素含量反演模型的建立是非常不利的。根據散點圖中紅邊位置與葉綠素含量的分布特征可知，最大一階導數法、拉格朗日三點內插法、倒高斯模型法和九階多項式擬合法在進行紅邊位置求解時均出現明顯的雙峰現象，即紅邊位置出現在中心波長分別為698和730 nm的集中區。與之相比，其他4種方法在不同程度上削減了雙峰現象的不利影響，且牛頓插值法的效果最佳。同時，線性四點內插法的紅邊位置向長波方向收攏（不是紅移，是收攏，整體集中在720～730 nm）；倒高斯模型法、線性外推法存在紅邊位置向短波方向位移（藍移）。

表2 不同紅邊位置求解算法的特性分析

注：進行不同紅邊位置求解算法的特性分析所用的數據來源于2002－2003年的試驗數據（樣本數為136）。

Note: Data used to analyze the characteristics of different Red Edge Positions (REP) solving algorithms come from the experimental data from 2002 to 2003 (The number of samples are 136).

注：用于制作散點圖的數據來源于2002－2003年的試驗數據（樣本數為136）。

經以上分析，可以得到以下結論：1）最大一階導數法、拉格朗日三點內插法的紅邊位置變幅最大（41 nm），對葉綠素含量十分敏感，然而存在明顯的雙峰現象；2）倒高斯模型法的結果介于695～729 nm之間，平均值最小（719.5 nm），整體向短波方向移動（藍移），且RE最大（0.882%）；3）線性四點內插法的結果介于717～731 nm之間，整體向長波方向收攏，平均值為725.4 nm，紅邊位置變幅最小（14 nm），對葉綠素含量變化不敏感；4）線性外推法的結果變幅較好（37.1 nm），但平均值偏低（721.9 nm），整體向短波方向移動（藍移），且RE較大（0.551%）；5）九階多項式擬合法的結果整體較好，不過雙峰現象最為顯著；6）牛頓八點插值法和牛頓-切比雪夫插值法的結果最小值偏大（705和706 nm），平均值、變幅和RE較為理想，并有效克服了雙峰現象和紅邊位移。

2.2 基于紅邊位置的冬小麥葉綠素含量反演

2.2.1 基于紅邊位置的冬小麥葉綠素含量反演模型的建立

利用2002－2003年冬小麥的葉綠素a（Chlorophyll-a，Chla）、葉綠素b（Chlorophyll-b，Chlb）和葉綠素a+b（Chlorophyll-ab，Chlab）含量實測數據與紅邊位置計算結果進行相關性分析，得到基于不同算法的冬小麥葉綠素含量反演模型（表3），研究表明8種方法計算出的紅邊位置和葉綠素含量均有很好的相關性。

表3 基于不同紅邊位置求解算法的冬小麥葉綠素含量反演模型

注：建立不同紅邊位置求解算法的冬小麥葉綠素反演模型的數據來源于2002—2003年的試驗數據（樣本數為136）。

Note: Data used to build the inversion models of winter wheat chlorophyll content based on different Red Edge Position (REP) solving algorithms are the experimental data from 2002 to 2003 (The number of samples are 136).

2.2.2 基于紅邊位置的冬小麥葉綠素含量反演模型的精度評定

同時，基于幾種算法建立的Chlab預測值與測量值的擬合方程式的斜率在0.406～1.191，截距在?1.524～3.118。其中，線性四點內插法和牛頓八點插值法的斜率分別為1.052和0.839，接近于1，截距分別為?0.101和0.378，接近于0，說明基于以上2種方法建立的反演模型的預測值能夠高精度接近葉綠素含量實測值。相比之下，牛頓-切比雪夫插值法存在系統誤差，其斜率在0.498～0.628之間，與理論值不相符。

表4 基于不同紅邊位置求解算法的冬小麥葉綠素含量反演模型的精度檢驗

注：反演模型的精度檢驗數據來源于2018－2019年的試驗數據（樣本數為64）。

Note: The accuracy test data of the inversion model comes from the test data from 2018 to 2019 (The number of samples are 64).

3 討 論

3.1 改進算法的性能評述

不難發現，與傳統方法相比，改進方法在改善雙峰現象和紅邊位移方面具有顯著優勢。線性四點內插法在改善雙峰現象上的效果同樣明顯，但其紅邊位置變幅在8種算法中最小（14 nm），對葉綠素含量變化不敏感，而且存在紅邊位置整體向長波收攏，這對葉綠素反演模型的建立是不利的。考慮到擬合方程的斜率與截距，牛頓-切比雪夫插值法便暴露出其缺點，這可能與切比雪夫零點的選取區間過大導致其出現系統誤差有關，下一步可以縮小零點區間，驗證切比雪夫零點替換等距節點的科學性。同時，牛頓-切比雪夫插值法對光譜分辨率有較高的要求，很大程度上限制了其應用范圍，而牛頓八點插值法不僅繼承了牛頓插值法的優越性，而且融合了線性四點內插法對插值點需求少的特性，具有更高的精度、穩定性和實用性。綜合考慮算法特性和各項精度指標，研究認為牛頓八點插值法優于傳統算法和牛頓-切比雪夫插值法。

3.2 基于紅邊參數的冬小麥葉綠素含量反演模型的評價

8種算法擬合方程的2（0.396～0.656）普遍低于反演方程的2（0.455～0.758）的現象，是由于2002－2003和2018－2019年試驗區內的冬小麥品種、長勢與生境狀況不完全相同，對模型的適用性產生了一定的影響。下一步可以在同一實驗條件下，設置建模組與驗證組，以期提高模型的精度。

4 結 論

本研究以冬小麥為研究對象，采用最小二乘回歸分析方法，在實現6種經典紅邊位置求解算法的基礎上，將牛頓插值法應用于紅邊位置求解，提出牛頓八點插值法和牛頓-切比雪夫插值法；并根據不同算法的紅邊位置分布特征綜合分析了新舊算法的優缺點，對比了不同算法的綜合屬性信息；此外，建立了基于紅邊參數的冬小麥葉綠素含量反演模型，進行了葉綠素含量反演模型的預測值與實測值的線性擬合分析。通過以上研究得出，1）本研究提出的牛頓八點插值法和牛頓-切比雪夫插值法不僅有效地改善了紅邊位置求解時的雙峰現象和紅邊位移，具有更高的精度與穩定性，而且使用牛頓插值法建立的葉綠素含量反演模型的決定系數（coefficient of determination，2）較其他方法高，牛頓八點插值法和牛頓-切比雪夫插值法的2分別為0.728和0.751；2）基于改進算法反演模型的預測值與葉綠素含量實測值的擬合方程2較其他方法高，均>0.619；標準均方根誤差（Normalized Root Mean Squared Error，NRMSE）較其他方法小，均<0.151，即模型優度較好；3）同時，研究表明牛頓八點插值法比牛頓-切比雪夫插值法具有更高的精度、穩定性和實用性；4）關于牛頓-切比雪夫插值法中零點區間的選取和利用切比雪夫零點替換等距節點的科學性還有待進一步研究。

根據本研究對不同紅邊算法的紅邊位置分布特征分析、算法的綜合屬性信息對比以及基于紅邊參數的冬小麥葉綠素含量反演模型的精度評定，可以認為將牛頓插值法應用到冬小麥的紅邊位置計算和葉綠素含量反演的效果優于傳統方法。此外，考慮到牛頓八點插值法計算紅邊位置只需要8個特定的光譜波段，為制作簡易的作物葉綠素含量測定傳感器提供了理論基礎。

[1]Velichkova K, Krezhova D. Extraction of the red edge position from hyperspectral reflectance data for plant stress monitoring [C]//10th Jubilee International Conference of the Balkan Physical Union, 2019.

[2]Ptushenko V, Ptushenko O, Tikhonov A. Chlorophyll fluorescence induction, chlorophyll content, and chromaticity characteristics of leaves as indicators of photosynthetic apparatus senescence in arboreous plants[J]. Biochemistry, 2014, 79(3): 260-272.

[3]Miller J R, Hare E W, WU J. Quantitative characterization of the vegetation red edge reflectance 1. An inverted-Gaussian reflectance model[J]. International Journal of Remote Sensing, 1990, 11(10): 1755-1773.

[4]陳西亮，張佳華. 牛頓插值法在植被紅邊擬合中的應用[J]. 湖北農業科學，2016，55(7)：1828-1831.

Chen Xiliang, Zhang Jiahua. Application of Newton interpolation method in red edge fitting of vegetation[J]. Hubei Agricultural Science, 2016, 55(7): 1828-1831. (in Chinese with English abstract).

[5]黃林生，江靜，黃文江，等. Sentinel-2影像和BP神經網絡結合的小麥條銹病監測方法[J]. 農業工程學報，2019，35(17)：178-185.

Huang Linsheng, Jiang Jing, Huang Wenjiang, et al. Sentinel-2 imaging and BP neural network monitoring method for wheat stripe rust[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(17): 178-185. (in Chinese with English abstract).

[6]Sun Yuanheng, Qin Qiming, Ren Huazhong, et al. Red-edge band vegetation indices for leaf area index estimation from Sentinel-2/MSI imagery[J]. IEEE Transactions on Geoscience and Remote Sensing, 2020, 58(2): 826-840.

[7]Curran P J, Dungan J L, Macler B A, et al. The effect of a red leaf pigment on the relationship between red edge and chlorophyll concentration[J]. Remote Sensing of Environment, 1991, 35(1): 69-76.

[8]黃文江，王紀華，劉良云，等. 冬小麥紅邊參數變化規律及其營養診斷[J]. 遙感技術與應用，2003(4)：206-211.

Huang Wenjiang, Wang Jihua, Liu Liangyun, et al. Variation of red border parameters in winter wheat and its nutritional diagnosis[J]. Remote sensing technology and Applications, 2003(4): 206-211. (in Chinese with English abstract).

[9]Demetriades-Shah T H, Steven M D, Clark J A. High resolution derivative spectra in remote sensing[J]. Remote Sensing of Environment, 1991, 33(1): 55-64.

[10]Dawson T P, Curran P J. Technical note A new technique for interpolating the reflectance red edge position[J]. International Journal of Remote Sensing, 1998, 19(11): 2133-2139.

[11]Clevers J G P W, De Jong S M, Epema G F, et al. The use of the MERIS standard band setting for deriving the red edge index[J]. International Journal of Remote Sensing, 2002, 23(16): 3169-3184.

[12]Pu R, Peng G, Biging G S, et al. Extraction of red edge optical parameters from Hyperion data for estimation of forest leaf area index[J]. IEEE Transactions on Geoscience and Remote Sensing, 2003, 41(4):916-921.

[13]Lin L, Ustin S L, Lay M. Application of AVIRIS data in detection of oil-induced vegetation stress and cover change at Jornada, New Mexico[J]. Remote Sensing of Environment, 2005, 94(1): 1-16.

[14]Cho M A, Skidmore A K. A new technique for extracting the red edge position from hyperspectral data: The linear extrapolation method[J]. Remote Sensing of Environment, 2006, 101(2): 181-193.

[15]何天榮. 牛頓插值法計算函數值的優缺點分析[J]. 考試與評價，2016(6)：70.

He Tianrong. Analysis of advantages and disadvantages of Newton interpolation method for calculating function values[J]. Examination and Evaluation, 2016(6): 70. (in Chinese with English abstract)

[16]曹中盛，李艷大，葉春，等. 基于高光譜的雙季稻分蘗數監測模型[J]. 農業工程學報，2020，36(4)：185-192.

Cao Zhongsheng, Li Yanda, Ye Chun, et al. Monitoring model of tiller number of double cropping rice based on hyperspectrum[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(4): 185-192. (in Chinese with English abstract).

[17]Ghule A, Deshmukh R R, Gaikwad C. MFDS-m red edge position detection algorithm for discrimination between healthy and unhealthy vegetable plants[J]. Recent Trends in Image Processing and Pattern Recognition, 2019, 7(1): 372-379.

[18]高興，李斐，楊海波，等. 基于紅邊位置的馬鈴薯植株氮濃度估測方法研究[J]. 植物營養與肥料學報，2019，25(2)：134-148.

Gao Xing, Li Fei, Yang Haibo, et al. Study on nitrogen concentration estimation method of potato plants based on red edge location[J]. Journal of Plant Nutrition and Fertilizer, 2019, 25(2): 134-148. (in Chinese with English abstract).

[19]丁永軍，張晶晶，李修華，等. 基于光譜紅邊位置提取算法的番茄葉片葉綠素含量估測[J]. 農業機械學報，2016，47(3)：292-297.

Ding Yongjun, Zhang Jingjing, Li Xiuhua, et al. Estimation of chlorophyll content in tomato leaves based on spectral red edge location extraction algorithm[J]. Journal of Agricultural Machinery, 2016, 47(3): 292-297. (in Chinese with English abstract)

[20]郝瑞娟，王周鋒，王文科，等. CO2脅迫下大豆葉片紅邊位置最優算法的研究[J]. 大氣與環境光學學報，2016，11(1)：51-59.

Hao Ruijuan, Wang Zhoufeng, Wang Wenke, et al. Study on the optimal algorithm of red edge position for soybean leaves under CO2stress[J]. Journal of Atmospheric and Environmental Optics, 2016, 11(1): 51-59. (in Chinese with English abstract)

[21]Endre S, Mayers D F. An Introduction to Numerical Analysis[M]. Cambridge: Cambridge University Press, 2003.

[22]余鵬，孫學金，趙世軍. GPS定位中衛星坐標計算的切比雪夫多項式擬合法[J]. 氣象科技，2004，32(3)：198，4.

Yu Peng, Sun Xuejin, Zhao Shijun. Chebyshev polynomial fitting method for satellite coordinate calculation in GPS positioning[J]. Meteorological Science and Technology, 2004, 32(3): 198, 4. (in Chinese with English abstract)

[23]任嘉衍，劉慧敏，丁圣彥，等. 伊河流域景觀格局變化及其驅動機制[J]. 應用生態學報，2017，28(8)：2611-2620.

Ren Jiayan, Liu Huimin, Ding Shengyan, et al. Changes of landscape pattern in Yi River basin and its driving mechanism[J]. Journal of Applied Ecology, 2017, 28(8): 2611-2620. (in Chinese with English abstract).

[24]于汧卉，楊貴軍，王崇倡. 地面高光譜和PROSAIL模型的冬小麥葉綠素反演[J]. 測繪科學，2019，44(11)：96-102.

Yu Qianhui, Yang Guijun, Wang Chongchang. Ground winter wheat chlorophyll of hyperspectral and PROSAIL model inversion[J]. Mapping Science, 2019, 44(11): 96-102. (in Chinese with English abstract).

[25]吳偉斌，李佳雨，張震邦，等. 基于高光譜圖像的茶樹LAI與氮含量反演[J]. 農業工程學報，2018，34(3)：195-201.

Wu Weibin, Li Jiayu, Zhang Zhenbang, et al. Inversion of tea tree LAI and nitrogen content based on hyperspectral images[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(3): 195-201. (in Chinese with English abstract)

Inversion of winter wheat chlorophyll contents based on improved algorithms for red edge position

Qian Binxiang1,2,3, Huang Wenjiang2,3※, Ye Huichun2,3, Kong Weiping4, Ren Yu1,3, Xing Naichen1,3, Jiao Quanjun2,3

(1100049; 2.,572029,; 3.,100094,; 4.,,100094,)

Red Edge Position (REP) of vegetation spectral reflectance is highly sensitive to chlorophyll content. The inversion model of crop chlorophyll content based on REP enables timely growth monitoring of the crops on a large scale. The displacement of REP and bimodal phenomenon are ubiquitous in 6 traditional algorithms of REP. To reduce the adverse effects effectively, the Newton interpolation method was applied to calculate REP in this study. And two improved REP solving algorithms, Newton-Chebyshev-Node Interpolation (REP_NCNI) and Newton Eight-Point Interpolation (REP_NEPI) were proposed. The strengths and weaknesses of the old and improved algorithms were analyzed, according to the distribution characteristics of REP from different algorithms, and the comprehensive attribute information of the different algorithms was compared. It was found that: 1) Maximum First Derivative (REP_MFD) method and Lagrange Three-Point Interpolation (REP_LAGR) method had the largest variation of REP (41 nm), which was sensitive to chlorophyll content, however, there was an obvious bimodal phenomenon. 2) The REP calculated by an Inverted Gaussian (REP_IG) model method ranged from 695 nm to 729 nm with the lowest mean value (719.5 nm). The whole model moved towards the short-wave direction (blue shift) with the highest Relative Error (RE) (0.882%). 3) The REP calculated by the Linear Four-Point Interpolation method (REP_LFPI) were between 717 and 731 nm, with an average value of 725.4 nm. The whole result was clustered in the direction of a long wave, and the REP had the smallest variation (14 nm), which was not sensitive to the change of chlorophyll content. 4) The variation of the Linear Extrapolation (REP_LE) method was better (39 nm), but the average value was lower (721.9 nm). The whole value moved towards the short-wave direction (blue shift), and the RE was larger (0.551%). 5) The results of polynomial fitting of the ninth order (REP_POLY) were generally good, but the bimodal phenomenon was the most serious. 6) The REP_NCNI and REP_NEPI overcame the bimodal phenomenon and displacement of REP effectively with ideal mean value, amplitude, and RE. And the least square regression was adopted to establish the inversion model of chlorophyll content of winter wheat based on REP. The study revealed that compared with traditionalalgorithms, the improved algorithms exhibited the most accurate and robust performance, where the coefficient of determination of the chlorophyll content inversion model established by improved algorithms was higher than that of traditional algorithms with the coefficient of determination of 0.728 and 0.751, respectively. Moreover, in the improved algorithms, the coefficient of determination between the predicted value and the measured value was greater than 0.619, which was 10.480% higher than that of the REP_MFD method, and the standard root mean square error was less than 0.151, indicating that the goodness of the model was better. At the same time, the coefficient of determination of the inverse equation (0.455-0.758) was higher than that of the fitting equation (0.396-0.656). And the inversion model was ranked chlorophyll-a model, chlorophyll-ab model, and chlorophyll-b model, according to the coefficient of determination from large to small. Besides, the research showed that in the two improved algorithms, the REP_NEPI demonstrated the best and satisfactory performance than REP_NCNI. Considering that only 8 bands were needed to calculate REP by REP_NEPI, it provided a theoretical basis for making a simple sensor to determine the chlorophyll content of crops. The results showed that REP_NEPI was the optimal selection for the calculation of REP and the inversion of chlorophyll content of winter wheat. And this study should provide theoretical and technical support for the inversion of biophysical and biochemical parameters of vegetation and the application in agricultural production.

algorithms; remote sensing; models; winter wheat; inversion of chlorophyll content; Newton-Chebyshev-node interpolation; Newton eight-point interpolation; red edge position

錢彬祥，黃文江，葉回春，等. 紅邊位置改進算法的冬小麥葉綠素含量反演[J]. 農業工程學報，2020，36(23)：162-170.doi：10.11975/j.issn.1002-6819.2020.23.019 http://www.tcsae.org

Qian Binxiang, Huang Wenjiang, Ye Huichun, et al. Inversion of winter wheat chlorophyll contents based on improved algorithms for red edge position[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(23): 162-170. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2020.23.019 http://www.tcsae.org

2020-09-24

2020-11-23

國家重點研發計劃課題（2016YFD0300601）；國家自然科學基金青年科學基金項目（41871339，41901369）；國家高層次人才特殊支持計劃項目（黃文江）

錢彬祥，博士生，主要從事植被定量遙感及應用研究。Email：2789281754@qq.com

黃文江，博士，研究員，博士生導師，主要從事植被定量遙感及應用研究。Email：huangwj@radi.ac.cn

10.11975/j.issn.1002-6819.2020.23.019

S127; S512.1+1

A

1002-6819(2020)-23-0162-09