趙向苗，高潤宏*，牛若愷，盧玉婷，高健斌

(1.內(nèi)蒙古農(nóng)業(yè)大學(xué) 林學(xué)院，內(nèi)蒙古 呼和浩特 010019；2.西鄂爾多斯國家級自然保護(hù)區(qū)管理局，內(nèi)蒙古 鄂爾多斯 017000)

沙冬青(Ammopiptanthusmongolicus)又名蒙古黃花木，屬于亞洲中部荒漠地區(qū)唯一的常綠闊葉灌木，是古地中海第三孑遺植物，作為阿拉善地區(qū)的主要建群種，對荒漠區(qū)生態(tài)穩(wěn)定起著不可替代的作用。沙冬青已被列入第1批珍稀瀕危保護(hù)植物名錄中[1]。到目前為止，沙冬青人工種植不斷擴(kuò)大；但由于人為破壞、城鎮(zhèn)礦區(qū)建設(shè)等干擾及受自身生境條件的影響，天然分布面積在不斷減少，處于更加瀕危狀態(tài)。因此，通過高光譜遙感技術(shù)掌握沙冬青的空間分布，對其生長狀態(tài)以及保護(hù)利用具有重要的生態(tài)意義。

高光譜技術(shù)作為高光譜遙感技術(shù)中不可分割的一部分，憑借其在“非破壞”的前提下，無損快速的獲取植物光譜信息來監(jiān)測植物生長狀態(tài)，為高光譜遙感監(jiān)測技術(shù)提供理論基礎(chǔ)。植被高光譜遙感在生態(tài)領(lǐng)域也具有重要作用，它涉及的應(yīng)用研究領(lǐng)域較為廣泛。主要集中植被形態(tài)及活力探測研究，土地覆蓋的監(jiān)測以及關(guān)于冠層水文狀態(tài)和生物化學(xué)參數(shù)的估計等方面。在研究植被變化的過程中，傳統(tǒng)的野外調(diào)查方法大部分都具有破壞性，不僅費(fèi)時費(fèi)力，而且受到野外條件限制；特別是對于一些難以到達(dá)的地方，例如高山、沙漠等，高光譜技術(shù)的出現(xiàn)使其快速、無損監(jiān)測植物生長具有更加現(xiàn)實的意義和優(yōu)勢[2]。由于綠色植物的形態(tài)特征、化學(xué)組成成分以及對物候差異的適應(yīng)等不同，因而其吸收和反射的電磁波也不相同，在獲得的光譜曲線中表現(xiàn)為不同的光譜特征。關(guān)于通過植物的光譜特征來反映其生長狀態(tài)，國內(nèi)外的學(xué)者已經(jīng)開展了大量調(diào)查和分析研究[3-7]。

國內(nèi)關(guān)于沙冬青的研究進(jìn)展成果較多，且研究分析主要以其生理特征[8-11]、種群擴(kuò)散[12-15]、抗逆性基因及其機(jī)理[16-19]等為基礎(chǔ)進(jìn)行研究，關(guān)于其高光譜遙感監(jiān)測研究較少。劉美等[20]通過對沙冬青高光譜特征提取分析，探究冬季識別新疆沙冬青群落主要物種的可行性以及2種沙冬青的光譜特征差異，取得了較好的成果。但由于氣候及野外條件影響，未能分析沙冬青葉片理化參量與光譜曲線之間的相關(guān)性?；谇叭搜芯炕A(chǔ)，本試驗以西鄂爾多斯自然保護(hù)區(qū)廣泛分布的沙冬青為研究對象，獲取不同齡級沙冬青冠層光譜，運(yùn)用馬氏距離法選取其原始光譜、一階微分光譜、二階微分光譜特征差異顯著波段進(jìn)行分析，通過分析選擇出識別不同齡級沙冬青光譜差異顯著的波段；將冠層原始及其微分光譜與葉片葉綠素含量進(jìn)行分析，通過相關(guān)性的大小來實現(xiàn)從宏觀尺度理解沙冬青的生長狀況與葉綠素含量間的關(guān)系。

1 研究區(qū)概況

本研究選擇西鄂爾多斯自然保護(hù)區(qū)為試驗點(diǎn)。該區(qū)域地理位置106°40′-107°44′E，40°13′-26°10′N，平均海拔1 091.6 m；土壤類型主要為棕鈣土；年平均氣溫9.3℃，年日照時數(shù)3 000 h，年平均降雨量162.4 mm，蒸發(fā)量高達(dá)3 481.1 mm[21]。經(jīng)調(diào)查得出自然保護(hù)區(qū)內(nèi)共有植物335種。其中，屬于國家級重點(diǎn)保護(hù)植物有沙冬青(Ammopiptanthusmongolicus)、半日花(Helianthemumsongaricum)、四合木(Tetraenamongolica)、蒙古扁桃(Prunusmongolica)、綿刺(Potaniniamongolica)、革包菊(Tugarinoviamongolica)、胡楊(Populuseuphratica)等7種。

2 材料與方法

2.1 材料與數(shù)據(jù)獲取

2.1.1 沙冬青冠層光譜數(shù)據(jù)測定 本研究通過室外測量結(jié)合室內(nèi)實驗來完成，光譜儀使用ASD便攜式地物光譜儀(ASD FieldSpec FR)，波長范圍350～2 500 nm，采樣間隔為1.4 nm，光譜分辨率為3 nm，視場角為25°。測量對象選擇西鄂爾多斯自然保護(hù)區(qū)野外自然分布的沙冬青。測定時每隔30 min對白板進(jìn)行校正。野外測量時間以2020年7月25日和9月25日為主。以野外測量操作及規(guī)范為基礎(chǔ)，選擇當(dāng)?shù)貢r間11:00-13:00天氣晴朗、無云，風(fēng)力<3級以及受野外自然條件干擾較小的條件下測量。在對不同齡級(冠幅：Ⅰ(幼齡)≤75 cm，75 cm<Ⅱ(中齡)≤225 cm，300 cm<Ⅲ(老齡))沙冬青冠層光譜測定時，每次每個齡級隨機(jī)選擇10株進(jìn)行測定，儀器每隔5秒記錄1次數(shù)據(jù)，采集過程中每個光譜值重復(fù)采集10次，將10組數(shù)據(jù)平均值作為樣品最終反射光譜值。

2.1.2 沙冬青葉片葉綠素含量測定 沙冬青葉片采集在光譜測量之后進(jìn)行，根據(jù)光譜測定植株，選取其測定區(qū)域內(nèi)齡級、大小一致的葉片于黑暗低溫條件下迅速帶回。將葉片擦拭干凈，去除掉其主脈快速剪碎混勻，稱取0.2 g(可視樣品葉綠素含量高低而增減用量)剪碎樣品置于80%無水乙醇10 mL溶液中，置于常溫避光環(huán)境24 h直至葉片全部褪綠后，用80%乙醇頂容至25 mL，此溶液即為葉綠素提取液。最后，通過紫外分光光度計分別測定葉綠素提取液在663、645 nm和440 nm波長處的吸光度(A)，以重復(fù)3次的平均值作為測定的最終結(jié)果。通過以下公式計算色素濃度:

Chl-a=12.7A663-2.69A645

(1)

Chl-b=22.9A645-4.68A663

(2)

Chl-ab=Ca+Cb=20.21A645+8.02A663

(3)

Car=4.7A400-0.27×(Chl-a+Chl-b)

(4)

Chl-t=濃度(mg/L)×提取液體積(mL)/樣品鮮重(g)×1000

(5)

式中,Chl-a,Chl-b,Chl-ab,Car分別表示色素提取液中葉綠素a、葉綠素b、葉綠素和類胡蘿卜素的濃度，mg/L；Chl-t表示單位面積的色素含量，mg/g；A664，A645及A440分別表示提取液在663、645 nm及440 nm處的吸光度。

2.2 研究方法

2.2.1 光譜一階微分 光譜曲線一階微分可以消除儀器本身以及野外條件影響造成光譜數(shù)據(jù)之間的系統(tǒng)誤差，通過光譜微分技術(shù)有助于去除部分受背景、噪聲影響雜亂無章的線性光譜對目標(biāo)光譜的影響，更好地提取目標(biāo)光譜信息[21]。

(6)

式中,λi表示每個波段的波長；R′(λi)為波長λ的一階導(dǎo)數(shù)。

2.2.2 光譜二階微分

(7)

式中，R′(λj)為波長j和波長j+2之間的二階微分光譜；R(λj)為第j波長原始光譜反射率；R(λj+1)為第j+1波長原始光譜反射率；R(λj+2)為波長j+2原始光譜反射率。

2.2.3 馬氏距離法

(8)

式中，MD(K)為不同光譜曲線在k波段的馬氏距離；Xk為不同光譜曲線在k波段的差值矩陣；μ為零點(diǎn)向量矩陣；∑為協(xié)方差矩陣[6]。不同齡級之間的馬氏距離在SPSS中計算。

3 結(jié)果與分析

3.1 不同齡級沙冬青冠層光譜特征

3.1.1 不同齡級沙冬青冠層原始光譜特征分析 對觀測到的沙冬青冠層光譜數(shù)據(jù)進(jìn)行分析，得出不同齡級沙冬青光譜變化曲線特征(圖1)。

圖1 不同齡級沙冬青光譜特征分析

由圖1可知，不同齡級沙冬青在整個夏季生長期內(nèi)光譜曲線整體上較接近，且不同齡級光譜曲線存在較為顯著的差異。在350～650 nm的可見光范圍內(nèi)，光譜反射曲線屬于平穩(wěn)上升狀態(tài)，在達(dá)到第1個峰值后隨即下降，該波段的數(shù)值與其它波段相比較低，表明該波段吸收的紫外光和藍(lán)紫光較多，則反射和透射較少。但由于此波段范圍內(nèi)葉綠素出現(xiàn)綠色強(qiáng)反射，產(chǎn)生“藍(lán)邊”和“黃邊”峰值，因此觀察到植物葉片為綠色。此波段幼齡沙冬青反射率明顯高于其他2個齡級，說明幼齡沙冬青體內(nèi)葉綠素含量較多。在750～950 nm波長范圍的“紅邊”區(qū)域內(nèi)，光譜反射率隨波長的增加而上升，曲線較陡幾乎處于直線上升狀態(tài)，其斜率的大小與植物單位面積葉綠素的含量有關(guān)。此波段反射率是獲取“紅邊”特征的重要信息，在該區(qū)域幼齡沙冬青反射率高于老、中齡，這是由于幼齡沙冬青葉綠素含量較多且呈水膠溶狀態(tài)，對其形成強(qiáng)烈的紅外反射造成光譜吸收特征的變化。在900～1 200 nm范圍“紅肩”區(qū)域內(nèi)，不同齡級沙冬青葉片由于受葉片內(nèi)部結(jié)構(gòu)和水分吸收，以及氧吸收的阻礙等因素影響，則表現(xiàn)出2個微弱的低谷。綜合來看，夏季沙冬青處于生長期，由于葉色差異較大，則幼齡與老齡、中齡的冠層反射率相差較大。

3.1.2 不同齡級沙冬青冠層微分光譜特征分析 不同齡級沙冬青冠層反射光譜曲線由于其葉片內(nèi)部結(jié)構(gòu)差異、土壤條件差異以及外界環(huán)境的不同均會產(chǎn)生不同的特征，通過一、二階微分處理后，可以盡可能的消除由于儀器自身以及外界環(huán)境條件對光譜曲線產(chǎn)生的影響，有利于更好地反映內(nèi)在特性。從圖2、圖3可看出，3個齡級沙冬青冠層光譜曲線在700～750 nm范圍的“紅邊”區(qū)域產(chǎn)生了明顯的小波峰，由此說明在此波段沙冬青內(nèi)部組織結(jié)構(gòu)細(xì)胞產(chǎn)生了明顯的變化。通過比較發(fā)現(xiàn)幼齡沙冬青峰值起伏變化較明顯，在719 nm波長處出現(xiàn)了極大峰值，這主要是由于幼齡沙冬青長勢的原因造成的。

圖2 沙冬青一階微分光譜特征

圖3 沙冬青二階微分光譜特征

3.2 馬氏距離法對光譜波段選擇

3.2.1 馬氏距離法對原始光譜波段選擇 不同齡級沙冬青冠層原始光譜曲線如圖1所示，通過對比得出其曲線總體呈現(xiàn)一致趨勢，但在某些波段依舊存在顯著差異，反映光譜特征信息差異主要集中在這些區(qū)域。圖4表示的是馬氏距離法對原始光譜差異波段選擇，馬氏距離值越大則表示差異越明顯，越容易進(jìn)行區(qū)分；進(jìn)行波段選擇時需遵循兩個基本原則：1)為有效剔除大量差異較小波段，馬氏距離值必須高于各個波段的平均值；2)為更好地與高光譜數(shù)據(jù)結(jié)合，進(jìn)行光譜識別時馬氏距離值高的波段必須連續(xù)出現(xiàn)10 nm以上；根據(jù)計算，原始光譜各波段馬氏距離平均值為2.9，在差異高的波段的馬氏距離大于平均值的基礎(chǔ)上，得出不同齡級原始光譜中差異顯著的波段有：350～376、337～432、550～567、734～860、1 125～1 200 nm；350～376、337～432 nm波長范圍位于紫外光和藍(lán)紫光區(qū)域，此波段為植物光譜的“藍(lán)邊”區(qū)域；550～567 nm波長范圍為葉綠素的典型吸收帶；734～860 nm波段屬于沙冬青“紅邊”位置；1 125～1 200 m波段屬于葉片水分強(qiáng)烈吸收帶附近；由此可見，運(yùn)用馬氏距離法對特征波段的選取包含了植物光譜基本特征。通過馬氏距離法對差異顯著波段反射率值每10 nm求平均(以下處理選擇相同的波段范圍)，共獲得45個光譜反射率平均值用于不同齡級識別。

圖4 馬氏距離法對原始光譜差異波段選擇

3.2.2 馬氏距離對一階微分光譜波段選擇結(jié)果 圖5所示為馬氏距離法對一階微分光譜差異波段選擇。根據(jù)計算，一階光譜各波段馬氏距離平均值為11，在差異高的波段的馬氏距離大于平均值的基礎(chǔ)上，得出不同齡級一階光譜中差異顯著的波段有685～699、700～749、756～769、929～939、1 119～1 171 nm，685～699、700～749、756～769 nm波段沙冬青光譜反射率隨著波長增加而急劇升高，光譜曲線形態(tài)較陡，接近“紅邊”位置；929～939 nm波段位于水和氧的吸收帶；1 119～1 171 nm位于植物水分強(qiáng)烈吸收帶附近，通過馬氏距離法對差異顯著波段反射率值每10 nm求平均(以下處理選擇相同的波段范圍)，共獲得14個光譜反射率平均值用于不同齡級識別。

圖5 馬氏距離法對一階微分光譜差異波段選擇

3.2.3 馬氏距離對二階微分光譜波段選擇結(jié)果 圖6所示為馬氏距離法對二階微分光譜差異波段選擇，根據(jù)計算，二階光譜各波段馬氏距離平均值為7.9，在差異高的波段的馬氏距離大于平均值的基礎(chǔ)上，得出不同齡級二階光譜中差異顯著的波段有714～724、755～764、928～941、961～970、990～1 000、1 119～1 130 nm；其中714～724、755～764 nm接近“紅邊”位置；928～941、961～970、990～1 000、1 119～1 130 nm均位于植物水分強(qiáng)烈吸收帶附近，通過馬氏距離法對差異顯著波段反射率值每10 nm求平均，共獲得7個光譜反射率平均值用于不同齡級識別。

圖6 馬氏距離法對二階微分光譜差異波段選擇

3.3 沙冬青冠層光譜與葉綠素含量相關(guān)性分析

3.3.1 不同齡級沙冬青冠層原始光譜與葉綠素含量相關(guān)分析 圖7表示夏季不同齡級沙冬青冠層原始光譜曲線與葉片葉綠素含量相關(guān)系數(shù)曲線圖。通過相關(guān)文獻(xiàn)查閱以及相關(guān)系數(shù)檢驗表驗證得知，在P=0.001水平上的相關(guān)性臨界值為0.408；表示在P=0.001水平上相關(guān)系數(shù)的絕對值>0.408時即呈現(xiàn)極顯著相關(guān)，以下同上。

通過圖7可知，3個齡級(老、中、幼)沙冬青葉片葉綠素含量與冠層光譜曲線相關(guān)系數(shù)基本呈現(xiàn)一致性，三者的相關(guān)系數(shù)均顯示出先降低而后增加的趨勢，總體上呈現(xiàn)負(fù)相關(guān)性。其中，在740～1 000 nm波長范圍內(nèi)中齡和幼齡的相關(guān)系數(shù)產(chǎn)生了較大的波動性，相關(guān)系數(shù)絕對值均<0.408，相關(guān)性不顯著。而當(dāng)波長在<740 nm的可見光區(qū)域，3個齡級沙冬青葉片葉綠素含量與冠層光譜曲線呈現(xiàn)出極顯著的相關(guān)性。當(dāng)波長處于425～740 nm范圍內(nèi)，在P=0.001水平下，老齡沙冬青葉綠素含量與冠層光譜反射率相關(guān)系數(shù)絕對值>0.5，呈現(xiàn)出極顯著相關(guān)性。其中，當(dāng)波長處于701 nm處時呈現(xiàn)出最好的相關(guān)性，相關(guān)系數(shù)為-0.569，在波長處于600 nm的相關(guān)性也極顯著，相關(guān)系數(shù)達(dá)到-0.565。中齡沙冬青葉綠素含量與冠層光譜反射率呈顯著性，在波長處于716 nm處時呈現(xiàn)最好的相關(guān)性，相關(guān)系數(shù)達(dá)到-0.640，在波長處于600 nm處的相關(guān)性也呈現(xiàn)極為顯著性，相關(guān)系數(shù)為-0.419。幼齡沙冬青在波長處于675 nm處呈現(xiàn)出最好的相關(guān)性，相關(guān)系數(shù)為-0.531，當(dāng)波長處于401 nm波段時的相關(guān)性也很大，相關(guān)系數(shù)達(dá)到-0.486。由此可知，老齡的敏感波長為701 nm和600 nm；中齡的敏感波長為716 nm和600 nm；幼齡的敏感波長位于675 nm和401 nm；三者的相關(guān)性顯著波段均為<740 nm的可見光范圍內(nèi)，由此可知可見光波段屬于沙冬青葉綠素光譜反應(yīng)的敏感波段。

圖7 沙冬青葉綠素含量與冠層光譜相關(guān)系數(shù)

3.3.2 不同齡級沙冬青冠層一階微分光譜與葉綠素含量相關(guān)系數(shù) 對不同齡級沙冬青葉綠素總量與冠層一階微分光譜進(jìn)行相關(guān)分析，其相關(guān)系數(shù)如圖8所示。

與原始光譜相關(guān)系數(shù)曲線相比較，經(jīng)過導(dǎo)數(shù)處理后的一階微分光譜曲線與葉片葉綠素含量的關(guān)系呈現(xiàn)出更高的相關(guān)性，且相關(guān)性顯著波段數(shù)呈現(xiàn)出明顯增多的趨勢(圖8a)。老齡沙冬青555～677 nm、中齡沙冬青656～669 nm和幼齡沙冬青698～750、679～687、620～635 nm波段區(qū)域內(nèi)，葉片葉綠素含量與一階微分光譜曲線呈現(xiàn)極顯著正相關(guān)關(guān)系。其中老齡在569 nm、中齡在658 nm及幼齡在683 nm最大相關(guān)系數(shù)分別為0.585、0.558、0.695。老齡在500～548、683～750 nm波長范圍內(nèi)，葉片葉綠素含量與一階微分光譜曲線呈現(xiàn)顯著負(fù)相關(guān)，在512 nm和683 nm處達(dá)到最大相關(guān)系數(shù)-0.618和-0.656。中齡沙冬青在531～550、670～700 nm波長范圍內(nèi)，葉片葉綠素含量與一階微分光譜曲線呈現(xiàn)出顯著負(fù)相關(guān)關(guān)系，在532 nm和693 nm處達(dá)到最大相關(guān)系數(shù)-0.560和-0.427。幼齡沙冬青在357～379、407～462、518～584、590～605 nm波長范圍內(nèi)，葉綠素含量與一階微分光譜曲線呈現(xiàn)出顯著負(fù)相關(guān)關(guān)系，在374 nm和597 nm處達(dá)到最大相關(guān)系數(shù)-0.687和-0.646。

圖8 沙冬青葉綠素含量與一階微分光譜相關(guān)系數(shù)

4 結(jié)論與討論

不同齡級沙冬青原始反射光譜曲線變化趨勢基本一致，總體符合植被光譜特征的規(guī)律性，且幼齡沙冬青反射率差異較明顯。沙冬青光譜隨著生長期的變化也在變化；整個生長期內(nèi)，老、中和幼齡的光譜曲線差異較大，在可見光(350～650 nm)范圍，植物光譜特征差異主要受各種色素含量的影響，在這其中葉綠素含量的多少呈現(xiàn)出最為重要的作用；此波段幼齡反射率明顯高于其他2個齡級，說明幼齡沙冬青體內(nèi)葉綠素含量較多，出現(xiàn)綠色強(qiáng)反射；而750～850 nm的“紅邊”區(qū)域，在該區(qū)域幼齡沙冬青反射率高于老、中齡，這是由于幼齡沙冬青葉綠素含量較多且呈水膠溶狀態(tài)，對其形成強(qiáng)烈的紅外反射造成光譜吸收特征的變化。近紅外波段范圍內(nèi)，植物光譜特征差異主要受植物葉片水分吸收的影響，此波段通常稱為水分吸收波段，此波段可能受葉片內(nèi)部結(jié)構(gòu)和水的吸收的影響，則出現(xiàn)兩個微弱的低谷。

通過對原始光譜進(jìn)行一二階微分變換結(jié)果表明:幼齡沙冬青對綠光的反射能力較強(qiáng)；相比較于幼齡，老、中齡對綠光的反射則呈現(xiàn)出較弱的狀態(tài)。幼齡沙冬青在“紅邊”波長范圍內(nèi)光譜反射率上升速度較快，在719 nm處出現(xiàn)極大峰值，而老、中齡上升速度較慢。這主要是由于幼齡沙冬青長勢的影響所致。

通過分析得出馬氏距離法對不同齡級沙冬青光譜差異顯著波段的選擇具有顯著性，且所選擇波段均為光譜特征波段。利用馬氏距離法對原始光譜和一二階微分處理，通過對比得出其曲線總體呈現(xiàn)一致趨勢，但在某些波段依舊存在顯著差異，反映光譜特征信息差異主要集中在這些區(qū)域。在原始光譜曲線中通過馬氏距離法選擇出區(qū)分不同齡級的45個波段平均值中，11個值屬于可見光波長范圍內(nèi)，34個存在于近紅外波長范圍。一階微分光譜中選取的14個波段中，2個存在于可見光波長范圍，12個存在于于近紅外波段。二階微分光譜曲線中選擇的7個波段全部存在于近紅外波段。由此，通過分析可得基于馬氏距離法識別不同齡級沙冬青顯著差異的波段大多處于近紅外波段范圍。

通過冠層光譜與葉綠素相關(guān)性分析發(fā)現(xiàn)3個齡級與葉綠素的敏感波長均為<740 nm的可見光波長范圍內(nèi)，老齡在701 nm波長處相關(guān)性最好，相關(guān)系數(shù)為-0.569，中齡在波長716 nm處相關(guān)性最好，相關(guān)系數(shù)為-0.640，幼齡波長在675nm處相關(guān)性最好，相關(guān)系數(shù)為-0.531。由此說明，可見光波段屬于沙冬青葉綠素光譜反應(yīng)的敏感波。為進(jìn)一步找出不同齡級葉片光譜曲線對葉綠素含量的敏感波段，對沙冬青葉片原始光譜曲線進(jìn)行不同變換(一階微分)處理，微分光譜與葉綠素相關(guān)性分析，得出變換后的敏感波段分別為：老齡555～677 nm、中齡656～669 nm、幼齡698～750、679～687、620～635nm波長范圍內(nèi)，沙冬青葉片葉綠素含量與冠層一階微分光譜的相關(guān)性達(dá)到極顯著正相關(guān)。與原始光譜相關(guān)系數(shù)曲線相比較，經(jīng)過導(dǎo)數(shù)處理后的一階微分光譜曲線與葉片葉綠素含量的關(guān)系呈現(xiàn)出更高的相關(guān)性，且相關(guān)性顯著波段數(shù)呈現(xiàn)出明顯增多的趨勢。

原始光譜、光譜一階微分、二階微分，馬氏距離法特征波段選取，光譜曲線特征與葉片葉綠素含量相關(guān)性分析等均為區(qū)分不同齡級沙冬青結(jié)構(gòu)、生長狀態(tài)以及動態(tài)監(jiān)測提供了理論依據(jù)[22-24]。本研究僅涉及夏季生長季典型荒漠植物沙冬青光譜特征，對于不同植物在不同季節(jié)的光譜特征以及同種植物在不同環(huán)境中光譜特征差異等問題有待進(jìn)一步研究。