基于地基激光雷達(dá)的葉傾角分布升尺度方法研究

蘇　偉　展郡鴿　李　靜　馬鴻元　吳代英　張　蕊

(1.中國農(nóng)業(yè)大學(xué)信息與電氣工程學(xué)院， 北京 100083; 2.環(huán)境保護(hù)部衛(wèi)星環(huán)境應(yīng)用中心， 北京 100094)

?

基于地基激光雷達(dá)的葉傾角分布升尺度方法研究

蘇偉1展郡鴿1李靜2馬鴻元1吳代英1張蕊1

(1.中國農(nóng)業(yè)大學(xué)信息與電氣工程學(xué)院， 北京 100083; 2.環(huán)境保護(hù)部衛(wèi)星環(huán)境應(yīng)用中心， 北京 100094)

地基激光雷達(dá)因其具有穿透力強(qiáng)，能夠提取植被冠層三維結(jié)構(gòu)信息的優(yōu)勢，是提取植被葉傾角分布(Leaf angle distribution, LAD)的理想數(shù)據(jù)源，因此將地基激光雷達(dá)數(shù)據(jù)與遙感影像結(jié)合獲取大尺度葉傾角分布結(jié)果頗具潛力。以河北省保定市北部4個(gè)縣為研究區(qū)，利用10個(gè)玉米樣地的地基激光雷達(dá)數(shù)據(jù)提取葉傾角分布結(jié)果，使用主成分正變換提取玉米實(shí)測葉傾角分布數(shù)據(jù)中信息量最大的前3個(gè)主成分，再利用神經(jīng)網(wǎng)絡(luò)模型對所提取的主成分與Landsat8反射率數(shù)據(jù)結(jié)合建立關(guān)系模型，然后將訓(xùn)練好的模型應(yīng)用于整個(gè)研究區(qū)進(jìn)行升尺度轉(zhuǎn)換，最后通過主成分逆變換，得到升尺度后平均葉傾角(Mean tilt angle, MTA)結(jié)果。對升尺度后LAD與實(shí)測LAD及升尺度后MTA與實(shí)測MTA進(jìn)行交叉驗(yàn)證，結(jié)果表明，升尺度MTA與實(shí)測MTA的驗(yàn)證精度(R2)為0.786 2，均方根誤差(RMSE)為3.04°。該結(jié)果表明，使用提取主成分方法建立光譜數(shù)據(jù)與葉傾角分布的關(guān)系模型從而達(dá)到升尺度轉(zhuǎn)換的目的具有可行性，模擬精度較高，且誤差較小。

地基激光雷達(dá); 葉傾角分布; 主成分分析; BP神經(jīng)網(wǎng)絡(luò); 升尺度

引言

地基激光雷達(dá)(Terrestrial laser scanning, TLS)是一種新型的主動(dòng)遙感技術(shù)，是傳統(tǒng)雷達(dá)與現(xiàn)代激光技術(shù)結(jié)合的產(chǎn)物，通過位置、距離、角度等觀測數(shù)據(jù)直接獲取植被表面點(diǎn)的三維坐標(biāo)，在高精度三維冠層信息的實(shí)時(shí)獲取方面有極大的優(yōu)勢[1]。近年來，諸多學(xué)者利用激光雷達(dá)數(shù)據(jù)提取植被冠層結(jié)構(gòu)信息，在提取單木參數(shù)、森林參數(shù)以及植被葉面積等方面取得了一定成果[2-7]。葉傾角分布(Leaf angle distribution, LAD)是影響植被冠層能量傳輸和輻射的重要植被結(jié)構(gòu)參數(shù)之一[8],并且決定了入射光的有效合成輻射在葉片中的分布情況，進(jìn)而直接影響植物的生產(chǎn)能力[9]。傳統(tǒng)的LAD直接測量方法是利用量角器與葉片表面直接接觸進(jìn)行量取[10]，費(fèi)時(shí)費(fèi)力。而常用的平均葉傾角(Mean tilt angle, MTA)間接測量方式有很大的隨機(jī)性[11-12]。近年來，也出現(xiàn)了使用攝影法對葉傾角進(jìn)行測量[13-16]，用于闊葉林和擴(kuò)展型草本植被的測量，但該方法是一種只能在局部和固定數(shù)量位置上使用的方法。

也有學(xué)者利用多角度光譜數(shù)據(jù)，從作物的結(jié)構(gòu)敏感性參數(shù)[17]和模型轉(zhuǎn)換來識(shí)別作物L(fēng)AD[18]。然而，由于缺少LAD的實(shí)測數(shù)據(jù)， LAD的估算結(jié)果精度沒有被直接驗(yàn)證。ZOU等[19]利用高光譜分辨率的成像光譜數(shù)據(jù)，結(jié)合PROSAIL模型分別使用藍(lán)光波長(479 nm)、紅光波長(663 nm)組合關(guān)系及紅邊波長(748 nm)對葉傾角建立非線性關(guān)系，經(jīng)過對攝影法所測量的6種植物平均葉傾角的驗(yàn)證，兩種方法的驗(yàn)證精度(決定系數(shù)R2)分別為0.34和0.60，均方根誤差(RMSE)分別為18.7°和11.4°，盡管誤差較大，但該方法為使用光譜數(shù)據(jù)計(jì)算植被平均葉傾角提供了一種新的思路。然而，該方法受限于成像光譜儀的視場范圍，只能得到區(qū)域上空間分布的反演結(jié)果，且高光譜數(shù)據(jù)并非是一種能簡單獲取且經(jīng)濟(jì)的遙感數(shù)據(jù)源，因此，目前對于直接建立常用光學(xué)遙感影像的光譜信息與植被葉傾角關(guān)系研究較少。而地基激光雷達(dá)具有提取詳細(xì)三維信息的優(yōu)勢，能夠比較精確地提取植被葉傾角，因此，探索一種結(jié)合植被冠層三維信息與常用遙感數(shù)據(jù)結(jié)合并反演大尺度葉傾角的方法頗具潛力。

本文基于地基激光雷達(dá)掃描葉傾角分布數(shù)據(jù)，將葉傾角分布作為一個(gè)角度分布集合數(shù)據(jù)，并對該集合進(jìn)行簡單劃分，得到10個(gè)葉傾角分布變量，通過主成分分析法，對葉傾角分布變量進(jìn)行綜合信息提取，將提取的綜合信息即主成分，利用神經(jīng)網(wǎng)絡(luò)模型建立其與光譜數(shù)據(jù)的關(guān)系，從而得到升尺度光譜-綜合信息關(guān)系模型，利用主成分逆變換法計(jì)算葉傾角分布結(jié)果及平均葉傾角，并使用交叉驗(yàn)證的方法對平均葉傾角分布升尺度結(jié)果進(jìn)行驗(yàn)證。

1　材料與方法

1.1研究區(qū)域概況

研究區(qū)位于河北省保定市北部，覆蓋涿州市、高碑店市、定興縣、易縣4個(gè)縣市(115°29′44.985″～116°14′17.019″E, 38°56′43.625″～39°35′53.829″N)，由于氣候溫和，四季分明，水熱條件適宜農(nóng)作物的種植，該地大田玉米廣泛種植春、夏玉米。研究中實(shí)測數(shù)據(jù)采集時(shí)間為2015年8月25—27日，正值夏玉米灌漿期，實(shí)測數(shù)據(jù)采樣點(diǎn)均勻分布于研究區(qū)范圍內(nèi)，每個(gè)采樣樣本數(shù)據(jù)的玉米點(diǎn)云掃描面積約為10 m×2 m，共獲取10個(gè)地基激光雷達(dá)玉米掃描數(shù)據(jù)樣本，圖1為研究區(qū)Landsat8 OLI遙感影像及樣本空間分布示意圖。

圖1　研究區(qū)及采樣點(diǎn)分布示意圖(Landsat8 OLI遙感影像, R:655 nm； G:561 nm； B:483 nm)Fig.1　Study area and LAD sample points distribution (Landsat8 OLI remote sensing image, R:655 nm; G:561 nm; B:483 nm)

1.2TLS玉米掃描數(shù)據(jù)獲取及預(yù)處理

使用FARO Focus 3D X330三維激光雷達(dá)掃描儀采集TLS玉米掃描數(shù)據(jù)，對目標(biāo)采用2～3點(diǎn)的掃描方式，在數(shù)據(jù)內(nèi)業(yè)處理中，利用同名點(diǎn)對點(diǎn)云數(shù)據(jù)進(jìn)行拼接配準(zhǔn)。盡管實(shí)驗(yàn)進(jìn)行時(shí)盡量避免了有風(fēng)天氣，但仍會(huì)受到微風(fēng)的干擾而導(dǎo)致噪點(diǎn)產(chǎn)生，因此，首先對拼接后的玉米點(diǎn)云數(shù)據(jù)進(jìn)行去噪處理，然后根據(jù)實(shí)際需求，對玉米的莖、穂等非葉片組織進(jìn)行剔除，得到玉米葉片點(diǎn)云數(shù)據(jù)。圖2為玉米點(diǎn)云數(shù)據(jù)去噪、去莖、去穂的結(jié)果圖。

圖2　預(yù)處理后激光雷達(dá)玉米掃描點(diǎn)云示意圖Fig.2　Results of preprocessing maize TLS cloud points

1.3Landsat8 OLI遙感影像獲取及預(yù)處理

為了使實(shí)測數(shù)據(jù)采集時(shí)間與衛(wèi)星過境時(shí)間相匹配，選取日期與實(shí)驗(yàn)時(shí)間最近的Landsat8 OLI影像作為實(shí)驗(yàn)數(shù)據(jù)，本文所采用的Landsat8數(shù)據(jù)采集時(shí)間為2015年8月22日。Landsat8過境周期為16 d，成像幅寬為185 km×185 km。Landsat8有2個(gè)主要載荷: 運(yùn)營性陸地成像儀(Operational land imager，OLI) 和熱紅外傳感器(Thermal infrared sensor， TIRS)。OLI傳感器有9個(gè)波段，其中包括一個(gè)15 m的全色波段，其余波段的空間分辨率為30 m；TIRS有2個(gè)熱紅外波段，空間分辨率為100 m。根據(jù)ZOU等[16]使用PROSAIL模型對葉傾角波段敏感性分析的結(jié)果，可知藍(lán)、綠、紅、近紅外波段與葉傾角變化具有較高的相關(guān)性，因此采用OLI傳感器采集的藍(lán)、綠、紅、近紅外4個(gè)波段對葉傾角進(jìn)行模型構(gòu)建。首先使用ENVI 5.1軟件對Landsat8 OLI數(shù)據(jù)進(jìn)行輻射定標(biāo)，再使用FLASHH模型進(jìn)行大氣校正得到地表反射率，并利用研究區(qū)邊界對遙感影像進(jìn)行剪裁，提取出研究區(qū)。

1.4研究方法

1.4.1基于TLS數(shù)據(jù)的玉米葉傾角分布提取

本研究計(jì)算葉傾角的方法為提取單點(diǎn)方向上的法向量[20]。整個(gè)冠層的激光雷達(dá)點(diǎn)云為一個(gè)點(diǎn)集合P=(p1,p2，…,pi)，在每一個(gè)單點(diǎn)上搜索鄰域內(nèi)的點(diǎn)，組成一個(gè)平面，求取該平面的法向量，通過法向量得到葉傾角。主要步驟如下：

(1)對于任何一個(gè)單點(diǎn)pi，其鄰域內(nèi)k鄰近點(diǎn)組成點(diǎn)集Pi=(pi1,pi2,…,pik)，利用點(diǎn)集Pi構(gòu)建平面Ti，ni為該平面的單位法向量，pi為向量ni的起點(diǎn)，計(jì)算式為

(1)

(2)通過主成分分析法求取單位法向量。計(jì)算點(diǎn)集P=(p1k,p2k,…,pik)的半正定協(xié)方差矩陣M，通過求解矩陣M的特征值λ1、λ2、λ3及其對應(yīng)的特征向量e1、e2、e3，找到最小特征對應(yīng)的特征向量，該向量即為單位法向量[2,21]。M計(jì)算式為

(2)

式中(pi-p)——列向量

(3)統(tǒng)一所有點(diǎn)的法向量方向。對于任何一個(gè)單點(diǎn)，其法向量的方向具有隨機(jī)性，可能指向夾角為180°的2個(gè)相反方向。使用的方法是在某一區(qū)域內(nèi)隨機(jī)選取一個(gè)單位法向量構(gòu)建黎曼面，使得此區(qū)域內(nèi)生成樹最小。

(4)通過步驟(1)～(3)計(jì)算出來的每個(gè)點(diǎn)法向量傾角范圍為[-90°,90°]，對負(fù)值葉傾角取絕對值處理，得到葉傾角范圍為[0°,90°]，對葉傾角頻率分布以5°為間隔進(jìn)行統(tǒng)計(jì)，得出葉傾角分布情況。

圖3即為10個(gè)實(shí)測點(diǎn)葉傾角[0°，90°]頻率分布結(jié)果。

圖3　實(shí)測點(diǎn)葉傾角頻率分布結(jié)果Fig.3　Measured leaf area frequency distribution results

1.4.2葉傾角分布主成分提取

主成分分析(Principal component analysis，PCA)是根據(jù)變量之間的相互關(guān)系，盡可能不丟失信息地利用若干個(gè)綜合性指標(biāo)表示多個(gè)變量的方法，在處理多元樣本數(shù)據(jù)時(shí)，常遇到多元變量之間存在相關(guān)關(guān)系的情形，使得數(shù)據(jù)的分析復(fù)雜化，使用主成分分析可以把多元變量轉(zhuǎn)換為少數(shù)幾個(gè)獨(dú)立的綜合變量[22]。葉傾角分布作為一組連續(xù)概率分布數(shù)據(jù)，包含了實(shí)測樣地玉米葉傾角的頻率分布信息，為了探索葉傾角分布與光譜數(shù)據(jù)的關(guān)系，驗(yàn)證了近紅外波段與基于葉傾角分布提取的平均值擬合模型，結(jié)果表明，光譜信息與葉傾角分布平均值之間無明顯相關(guān)關(guān)系(圖4)，因此直接使用光譜數(shù)據(jù)模擬葉傾角方法是不可行的。但在輻射傳輸理論中，葉傾角的分布與冠層的光分布模式密切相關(guān)[23]。因此，需要選擇其他方法來提取葉傾角分布信息并建立其與光譜數(shù)據(jù)之間的關(guān)系。

圖4　基于近紅外波段模擬MTA結(jié)果Fig.4　Simulated MTA result based on infrared band

綜上，將葉傾角分布所包含的信息作為數(shù)據(jù)集合，不同角度范圍所對應(yīng)的頻率作為數(shù)據(jù)集合中的分量，可以通過主成分分析法對該數(shù)據(jù)中所包含信息進(jìn)行綜合提取，建立綜合信息與玉米冠層反射率的關(guān)系模型。首先，將葉傾角頻率分布結(jié)果以5°為間隔劃分為18個(gè)葉傾角頻率分布變量，對18個(gè)變量進(jìn)行主成分分析，但根據(jù)實(shí)測玉米葉傾角分布結(jié)果來看，葉傾角主要分布于[40°,90°]，[0°,40°]葉傾角變量所包含的信息極少，所以將主成分分析對象改為以5°為間隔、范圍為[40°,90°]的10個(gè)變量。通過主成分分析結(jié)果可知(表1)，前3個(gè)主成分包含葉傾角分布信息中95%以上的信息量。

表1　葉傾角分布成分主成分分析結(jié)果Tab.1　Results of LAD using PCA

1.4.3基于BP神經(jīng)網(wǎng)絡(luò)模型的LAD升尺度轉(zhuǎn)換

BP(Back propagation)神經(jīng)網(wǎng)絡(luò)是一種按照誤差逆向傳播算法訓(xùn)練的多層前饋網(wǎng)絡(luò)，是目前應(yīng)用最廣泛的神經(jīng)網(wǎng)絡(luò)模型之一[24]。BP神經(jīng)網(wǎng)絡(luò)的基本結(jié)構(gòu)一般包含3層:輸入層、隱含層(也稱中間層)和輸出層;i表示輸入層的神經(jīng)元數(shù)目,j表示隱含層的神經(jīng)元數(shù)目,k表示輸出層的神經(jīng)元數(shù)目。每一層上的神經(jīng)元都通過權(quán)值與相鄰層上的各神經(jīng)元實(shí)現(xiàn)全連接。

本文在Matlab編程環(huán)境下建立LAD的升尺度模型，基于主成分分析所提取成分以及實(shí)測點(diǎn)對應(yīng)的冠層反射率，設(shè)定訓(xùn)練數(shù)據(jù)的輸入層參數(shù)為實(shí)測點(diǎn)玉米冠層對應(yīng)藍(lán)、綠、紅、近紅外4個(gè)波段反射率，輸出層設(shè)置為第1～10主成分，以神經(jīng)網(wǎng)絡(luò)模型作為訓(xùn)練模型將2組數(shù)據(jù)建立關(guān)系模型。訓(xùn)練結(jié)果顯示，貢獻(xiàn)率較大的成分，其訓(xùn)練精度(決定系數(shù)R2)達(dá)到0.85～0.99，而對于信息貢獻(xiàn)率較小的成分，其訓(xùn)練精度較低，僅為0.1～0.3。將基于10個(gè)實(shí)測葉傾角分布數(shù)據(jù)訓(xùn)練好的升尺度模型應(yīng)用于整個(gè)研究區(qū)，從而實(shí)現(xiàn)地面離散實(shí)測結(jié)果到空間連續(xù)分布反演結(jié)果的升尺度轉(zhuǎn)換。

2　結(jié)果驗(yàn)證與分析

2.1模擬-實(shí)測葉傾角分布的交叉驗(yàn)證結(jié)果

為了對模擬結(jié)果進(jìn)行驗(yàn)證，利用留一法交叉驗(yàn)證[25]對10個(gè)實(shí)測平均葉傾角數(shù)據(jù)進(jìn)行模擬結(jié)果驗(yàn)證。留一法交叉驗(yàn)證方法是：從10個(gè)實(shí)測數(shù)據(jù)中選擇1個(gè)數(shù)據(jù)作為驗(yàn)證數(shù)據(jù); 然后使用剩下的實(shí)測數(shù)據(jù)構(gòu)建模型，并用最先被排除的實(shí)測值來驗(yàn)證該模型的精度，如此重復(fù) 10次。利用神經(jīng)網(wǎng)絡(luò)模型建立了Landsat8 OLI影像波段反射率與主成分之間的模型，為了驗(yàn)證模型精度，基于正變換的協(xié)方差矩陣以及特征向量，對模擬的主成分結(jié)果進(jìn)行主成分逆變換，得到升尺度后每一像素所對應(yīng)的葉傾角分布來實(shí)現(xiàn)模擬結(jié)果與實(shí)測數(shù)據(jù)的驗(yàn)證，并將模擬的葉傾角分布結(jié)果與實(shí)測葉傾角分布曲線疊加，其對比結(jié)果如圖5所示。

圖5　升尺度和實(shí)測葉傾角分布結(jié)果對比Fig.5　Contrast result of upscaled and measured LAD

由升尺度葉傾角分布結(jié)果(圖5)可知，升尺度葉傾角分布結(jié)果與實(shí)測葉傾角分布結(jié)果一致，頻率分布最大即葉傾角分布的最大值范圍在[65°,75°]，但不同實(shí)測點(diǎn)會(huì)有前移或后移的趨勢，因?qū)崪y點(diǎn)選取是在研究區(qū)范圍內(nèi)隨機(jī)選擇，玉米品種、種植方式會(huì)導(dǎo)致各實(shí)測點(diǎn)之間存在差異，但總體角度分布偏差較小，較為一致。10個(gè)實(shí)測點(diǎn)的擬合結(jié)果在[40°,60°]范圍內(nèi)，實(shí)測葉傾角分布與升尺度葉傾角分布結(jié)果良好，說明該角度分布范圍內(nèi)所提取的主成分信息量充足，玉米冠層在[40°,60°]葉傾角分布情況下，冠層對光譜反射率呈現(xiàn)出較為穩(wěn)定的反射信息，這與大田玉米均一的冠層結(jié)構(gòu)較為吻合。

2.2升尺度平均葉傾角反演結(jié)果

為了對升尺度后葉傾角分布結(jié)果與實(shí)測葉傾角分布結(jié)果進(jìn)行定量化驗(yàn)證，對上述2種結(jié)果計(jì)算平均葉傾角，升尺度平均葉傾角的相對誤差分布為0.25°～6.62°，絕對誤差為0.40%～11.45%(表2)，盡管不同品種玉米和不同種植方式會(huì)影響葉傾角的分布狀況，但玉米作為冠層結(jié)構(gòu)均一的大田作物，理論上不應(yīng)出現(xiàn)類似于實(shí)測點(diǎn)4的極大誤差情況，究其原因如下：在玉米激光點(diǎn)云預(yù)處理時(shí)，會(huì)根據(jù)玉米葉片的掃描情況，盡量完整地保留完整玉米葉片形態(tài)信息，但受到風(fēng)和玉米葉片之間遮擋的影響，會(huì)導(dǎo)致玉米將葉傾角較小的葉尖剔除，而其他實(shí)測數(shù)據(jù)中這種數(shù)據(jù)處理誤差較小，從而導(dǎo)致了實(shí)測點(diǎn)4與其他點(diǎn)平均葉傾角差異較大的情況。但總體來說，實(shí)測MTA與升尺度MTA的驗(yàn)證精度(決定系數(shù)R2)達(dá)到了0.786 2，均方根誤差(RMSE)為3.04°，該結(jié)果表明，升尺度模型精度比較理想。

表2　實(shí)測與升尺度MTA誤差結(jié)果Tab.2　Results of difference between measured MTA and upscaled MTA

圖6　升尺度與實(shí)測MTA結(jié)果對比Fig.6　Contrast result of upscaled and measured MTA

2.3平均葉傾角上推空間分布結(jié)果

為了使葉傾角分布升尺度數(shù)據(jù)能夠被定量化地應(yīng)用于物理模型(如PROSAIL模型)研究中，使用主成分逆變換法，求取研究區(qū)內(nèi)每一像素上的葉傾角分布對應(yīng)的葉傾角平均值，葉傾角的空間分布結(jié)果如圖7所示。從圖7可以看出，平均葉傾角集中分布于[55°,70°]，與該生育期玉米的平均葉傾角實(shí)測結(jié)果一致，玉米葉傾角分布由北向南逐漸增大，也體現(xiàn)出作物生長的物候特征，南部成熟較早，因此平均葉傾角較大。

圖7　玉米冠層平均葉傾角的升尺度轉(zhuǎn)換結(jié)果Fig.7　Upscaled MTA results of maize canopy

3　結(jié)論

(1)通過實(shí)驗(yàn)驗(yàn)證，直接建立Landsat8 OLI數(shù)據(jù)與平均葉傾角之間關(guān)系的方法是不可行的，但葉傾角與植被冠層反射率之間的密切關(guān)系是客觀存在的，因此需要將葉傾角信息通過其他方式進(jìn)行表達(dá)，找出與冠層反射率呈強(qiáng)相關(guān)性的表達(dá)方式，因此將葉傾角分布看作連續(xù)分布的概率數(shù)據(jù)集合，并對該數(shù)據(jù)集合進(jìn)行劃分，形成合適數(shù)量的變量，并利用主成分分析法，提取變量的綜合信息，進(jìn)而代替葉傾角分布數(shù)據(jù)原始的數(shù)據(jù)表達(dá)方式，使其能夠建立起與冠層反射率之間的關(guān)系，結(jié)果表明該方法可行。

(2)利用神經(jīng)網(wǎng)絡(luò)模型將地面離散的實(shí)測葉傾角分布數(shù)據(jù)進(jìn)行尺度上推反演，從而得到誤差小、效果穩(wěn)定的全區(qū)域葉傾角分布尺度上推結(jié)果。

(3)考慮研究中對葉傾角分布數(shù)據(jù)的實(shí)際應(yīng)用，通過主成分逆變換法，將通過神經(jīng)網(wǎng)絡(luò)模型模擬得到的表達(dá)葉傾角信息的主成分?jǐn)?shù)據(jù)，轉(zhuǎn)換為每一像素上的葉傾角分布，進(jìn)而將葉傾角分布定量化反演得到每一像素(30 m×30 m)平均葉傾角。通過驗(yàn)證結(jié)果可知，平均葉傾角的驗(yàn)證精度較高，誤差較小。

1AHMED O S, FRANKLIN S E, WULDER M A. Characterizing stand-level forest canopy cover and height using Landsat time series, samples of airborne LiDAR, and the random forest algorithm[J]. ISPRS Journal of Photogrammetry & Remote Sensing, 2015, 101: 89-101.

2ZHENG G, MOSKAL L M. Spatial variability of terrestrial laser scanning based leaf area index [J]. International Journal of Applied Earth Observation & Geoinformation, 2012, 19(10): 226-237.

3XU Guangcai, PANG Yong, LI Zengyuan. Classifying land cover based on calibrated full-waveform airborne light detection and ranging data[J]. Chinese Optics Letters, 2013, 11(8): 87-92.

4趙陽,余新曉. 三維激光技術(shù)無損測量林分結(jié)構(gòu)因子適用性研究[J]. 激光與光電子學(xué)進(jìn)展, 2013, 50(7): 35-41.ZHAO Yang,YU Xinxiao. Applicability of 3D laser technology in noninvasive measurement of stand structure in forestry [J]. Laser & Optoelectronics Progress, 2013, 50(7): 35-41.(in Chinese)

5鄭莎莎,董品亮,王成,等.基于LiDAR點(diǎn)云數(shù)據(jù)的樹冠空隙度指數(shù)分析[J].國土資源遙感,2014,26(4):103-110.

ZHENG Shasha, DONG Pinliang, WANG Cheng, et al. Lacunarity analysis of LiDAR point clouds for tree crowns [J].Remote Sensing for Land and Resources,2014,26(4):103-110.(in Chinese)

6蘇偉,郭皓,趙冬玲,等. 基于優(yōu)化PROSAIL葉傾角分布函數(shù)的玉米LAI反演方法[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2016, 47(3):234-241，271.

SU Wei, GUO Hao, ZHAO Dongling, et al. Leaf area index retrivel for maize canopy using optimized leaf angle distribution function of PROSAIL model[J]. Transactions of the Chinese Society for Agricultural Machinery,2016,47(3):234-241，271.(in Chinese)

7蘇偉,展郡鴿,張明政,等. 基于機(jī)載LiDAR數(shù)據(jù)的農(nóng)作物葉面積指數(shù)估算方法研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2016, 47(3):272-277.

SU Wei, ZHAN Junge, ZHANG Mingzheng, et al. Estimation method of crop leaf area index based on airborne LiDAR data[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(3):272-277.(in Chinese)

8ROSS J. The radiation regime and architecture of plant stands [M]. Dordrecht: Springer Science & Business Media, 2012 .

9GUTSCHICK V P. Joining leaf photosynthesis models and canopyphoton-transport models [M]∥Vanderbilt V C, Grant L.Photon-Vegetation Interactions, Springer Berlin Heidelberg, 1991: 501-535.

10CAMPBELL G S, NORMAN J M. An introduction to environmental biophysics [M]. 2nd ed. New York: Springer, 1998.

11WELLES J M, NORMAN J M. Instrument for indirect measurement of canopy architecture [J]. Agronomy Journal, 1991, 83(5): 818-825.

12WELLES J M, COHEN S. Canopy structure measurement by gap fraction analysis using commercial instrumentation[J]. Journal of Experimental Botany, 1996, 47(9):1335-1342.

14PISEK J, RYU Y, ALIKAS K. Estimating leaf inclination and G-function from leveled digital camera photography in broadleaf canopies [J]. Trees, 2014, 25(5): 919-924.

15PISEK J, SONNENTAG O, RICHARDSON A D. Is the spherical leaf inclination angle distribution a valid assumption for temperate and boreal broad leaf tree species? [J]. Agricultural & Forest Meteorology, 2013, 169: 186-194.

16HUANG W J, NIU Z, WANG J H. Identifying crop leaf angle distribution based on two-temporal and bidirectional canopy reflectance [J]. IEEE Transactions on Geoscience and Remote Sensing, 2006,44(12): 3601-3609.

17JACQUEMOUD S, VERHOEF W, BARET F. PROSPECT + SAIL models: a review of use for vegetation characterization[J]. Remote Sensing of Environment, 2009,113: S56-S66.

18BACOUR C, JACQUEMOUD S, LEROY M. Reliability of the estimation of vegetation characteristics by inversion of three canopy reflectance models on airborne POLDER data[J]. Agronomie, 2002, 22(6): 555-565.

LIU Ting. Invershion of leaf area index of maize at the key growth stages based on PROSAIL model and multi-source data [D]. Beijing: China Agricultural University, 2015.(in Chinese)

21馬利霞，鄭光，何維.葉方向三維空間分布的地面激光雷達(dá)反演與分析[J].遙感學(xué)報(bào), 2015,19(4):609-617.

MA Lixia，ZHENG Guang，HE Wei.Retrieving three-dimensional spatial distribution of leaf orientation using terrestrial LiDAR data[J]. Journal of Remote Sensing, 2015, 19(4): 609-617.(in Chinese)

22密長林,馬愛功,張曉東. 主成分分析在遙感影像數(shù)據(jù)中的實(shí)例應(yīng)用[J].山東國土資源,2013,29(7):69-71,76.

MI Changlin, MA Aigong, ZHANG Xiaodong. Practical application of principal component analysis in remote sensing image processing [J]. Shandong Land and Resources, 2013, 29(7): 69-71,76.(in Chinese)

23張軍,王一鳴,趙燕東. 基于橢球函數(shù)的棉花葉傾角分布動(dòng)態(tài)模擬[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2009,40(4):157-160,183.

ZHANG Jun, WANG Yiming, ZHAO Yandong. Dynamic simulation of leaf inclination angle distribution based on ellipsoidal function[J]. Transactions of the Chinese Society for Agricultural Machinery, 2009, 40(4):157-160,183.(in Chinese)

24姜方方,基于地理要素趨勢面的LAI尺度轉(zhuǎn)換方法研究[D].北京：中國農(nóng)業(yè)大學(xué),2015.

JIANG Fangfang. Research on the method of scale transformation for LAI based on geographic trend surface [D]. Beijing: China Agricultural University, 2015.(in Chinese)

25STONE M. Cross-validatory choice and assessment of statistical predictions [J]. Journal of the Royal Statistical Society: Series B (Methodological), 1974, 36(2):111-147.

Upscaling Leaf Angle Distribution Using Terrestrial Laser Scanning Technique

Su Wei1Zhan Junge1Li Jing2Ma Hongyuan1Wu Daiying1Zhang Rui1

(1.CollegeofInformationandElectricalEngineering,ChinaAgriculturalUniversity,Beijing100083,China2.SatelliteEnvironmentCenter,MinistryofEnvironmentalProtection,Beijing100094,China)

Leaf angle distribution (LAD) can be used to describe the canopy structure of vegetation completely, such as crops, trees and grass. It’s one of the important parameters to quantitative description of vegetation canopy structure. At the present, there are few studies used the spectral data to inverse LAD, and results of the most existing studies of mean leaf tilt angle and leaf angle distribution were the locational inversion. Therefore, this study set the study site in five counties of Baoding City, Hebei Province, using terrestrial laser scanning (TLS) to acquire the leaf angle distribution data of maize. Combining the Landsat8 remote sensing data, firstly, the principle component analysis was taken to extract the principle information of measured leaf angle distribution of maize. Secondly, the back propagation artificial neural network was taken to model the relationship of principal information and spectral data. Then, the model was used in the whole study area to accomplish the upscaling transform. Finally, the upscaled mean tilt angel (MTA) was calculated based on the predicted LAD by principal component inverse transformation, in order to quantitate the leaf angle data. The cross validation result showed that the accuracy (R2) between upscaled MTA and measured MTA was 0.786 2, and the mean square root error (RMSE) was 3.04°. Consequently, it shows that this method can realize the aim of LAD upscaling.

terrestrial laser scanning; leaf angle distribution; principal components analysis; back propagation neural network; upscaling

10.6041/j.issn.1000-1298.2016.09.026

2016-01-28

2016-03-16

國家自然科學(xué)基金項(xiàng)目(41371327)

蘇偉(1979—)，女，副教授，博士生導(dǎo)師，主要從事農(nóng)業(yè)遙感應(yīng)用研究，E-mail： suwei@cau.edu.cn

李靜(1978—)，女，高級(jí)工程師，主要從事生態(tài)遙感監(jiān)測與應(yīng)用研究，E-mail: li_jingly@163.com

TN959.3； Q948.1

A

1000-1298(2016)09-0180-06