1研究區(qū)概況

黃河三角洲位于渤海灣南岸和萊州灣西岸，主要分布在中國(guó)山東省東營(yíng)市和濱州市境內(nèi)，屬暖溫帶?半?濕?潤(rùn)?大?陸?性?季?風(fēng)?氣?候?，年?平?均?氣?溫?11.7~12.6?℃，日照充足，雨熱同期。該地區(qū)土壤以濱海鹽土和濱海潮土為主，養(yǎng)分含量低，且鹽漬化嚴(yán)重，不利于喬木生長(zhǎng)。自?20?世紀(jì)?80?年代以來(lái)，刺槐作為主要的造林樹(shù)種在該地區(qū)廣泛種植，逐漸形成以孤島、軍馬場(chǎng)、黃河故道和大汶流自然保護(hù)區(qū)?4?個(gè)林場(chǎng)為主體的種植區(qū)。本研究以孤島林場(chǎng)為研究區(qū)（圖1）。

圖?1?研究區(qū)及樣方分布（以?Li?DAR?數(shù)據(jù)為背景）

2材料與方法

2.1?數(shù)據(jù)，

2.1.1?野外生態(tài)調(diào)查數(shù)據(jù)

本研究于?2017?年?6?月、2019?年?6?月和?2019?年?10?月共獲取了?30?個(gè)野外調(diào)查樣方。根據(jù)兩級(jí)采樣設(shè)計(jì)原則，在?30?m×30?m?的樣方的?4?個(gè)角和中心布置了?5?個(gè)?10?m×10?m?的次樣方。利用?GPS?記錄每個(gè)樣方的位置。根據(jù)美國(guó)林業(yè)局樹(shù)冠狀況分類(lèi)指南對(duì)樣方進(jìn)行健康分級(jí)，具體方法如下：在?5?個(gè)次樣方中選擇?1?個(gè)標(biāo)準(zhǔn)木材評(píng)估?5?個(gè)樹(shù)冠活力指標(biāo)，分別為活冠比、冠幅、郁閉度、林冠枯梢率、樹(shù)葉透光度（表?1），取?5?個(gè)次樣方指標(biāo)的平均值作為樣方的指標(biāo)值，最后根據(jù)森林健康等級(jí)判斷規(guī)則（表?2）得到樣方的健康等級(jí)。30個(gè)樣方被分為?3?類(lèi)，包括?14?個(gè)健康樣方、10?個(gè)亞健康樣方和?6?個(gè)不健康樣方。

?

2.1.2?無(wú)人機(jī)?Li?DAR?數(shù)據(jù)

本研究于?2017?年?6?月采集了黃河三角洲孤島林場(chǎng)區(qū)域的無(wú)人機(jī)?Li?DAR數(shù)?據(jù)?。無(wú)?人?機(jī)?上?搭?載激?光?掃?描?儀?、IMU（Inertial?measurement?unit）和?GPS（Global?posi?tion?system），飛行高度為?120?m，飛行速度為?4.8?m/s，巡航半徑為?2?km。激光波長(zhǎng)采用近紅外波段，光斑直徑約為?50?mm，獲取的平均點(diǎn)密度為?40?個(gè)/m2，采用?WGS84?坐標(biāo)系和?UTM?投影。

2.2?Li?DAR?數(shù)據(jù)預(yù)處理

首先采用基于統(tǒng)計(jì)的點(diǎn)云去噪方法（Statisticaloutlier?removal）對(duì)原始的無(wú)人機(jī)?Li?DAR?數(shù)據(jù)進(jìn)行去?噪?。然?后?通?過(guò)?改?進(jìn)?的?漸?進(jìn)?三?角?網(wǎng)?加?密?算?法（IPTD）對(duì)去噪后的?Li?DAR?點(diǎn)云進(jìn)行濾波，將點(diǎn)云分為地面點(diǎn)和非地面點(diǎn)。利用地面點(diǎn)插值生成數(shù)字地形模型?DTM，并用?DTM?對(duì)點(diǎn)云數(shù)據(jù)進(jìn)行歸一化消除地形的影響。再利用歸一化后的點(diǎn)云生成數(shù)字表面模型?DSM，減去?DTM?得到冠層高度模型?CHM。本研究為了獲取詳細(xì)的冠層變化，將?CHM?的空間分辨率設(shè)置為?1?m2。

2.3?Li?DAR?特征變量提取

2.3.1?垂直方向的變量?

本研究的垂直變量分為?2個(gè)部分，包括基于歸一化?Li?DAR?點(diǎn)云中提取的特征變量和基于冠層垂直剖面提取的特征變量?；邳c(diǎn)云提取的變量包括?3?個(gè)。高度百分位數(shù)（H99）是所有回波點(diǎn)按高度排序后，99%?的點(diǎn)云所在的高度，可以從整體上反映樣方內(nèi)樹(shù)木的高度。垂直復(fù)雜度（VCI）可以量化點(diǎn)云的垂直分布，其數(shù)值區(qū)間為?0~1，0?表示點(diǎn)云在垂直方向上分布不均勻，1?表示分布均勻，公式為：

式中，n?為垂直方向上劃分的高度層數(shù)，Pi為每個(gè)高程層中點(diǎn)云占點(diǎn)云總數(shù)的比例（%）。葉面積密度的變異系數(shù)（LADcv）可以反映冠層內(nèi)葉面積指數(shù)隨高度變化的變異情況，公式為：

式中，H?為垂直方向劃分的高度層數(shù)，GF?為間隙率，k?為消光系數(shù)，nab為第?i?個(gè)高度層以上的點(diǎn)云數(shù)量，nbe為第?i?個(gè)高度層以下的點(diǎn)云數(shù)量，ntot為點(diǎn)云總數(shù)。冠層垂直剖面可以反映包括枝葉在內(nèi)的整個(gè)樹(shù)冠在垂直方向上的點(diǎn)云密度分布情況，由某一高度層的點(diǎn)云與點(diǎn)云總數(shù)的比值表示。本研究采用Weibull?密度函數(shù)擬合冠層垂直剖面，以獲取尺度參數(shù)?α（weibull_α）和形狀參數(shù)?β（weibull_β），公式為：

式中，f（x）是點(diǎn)云密度，x是高度百分位數(shù)。

2.3.2?水平方向的變量?

冠層高度模型?CHM?描述了樹(shù)冠的表面形態(tài)。樹(shù)木枯梢會(huì)增加激光脈沖的穿透率，并在?CHM?的水平方向產(chǎn)生不同大小的不規(guī)則斑塊，因此可以通過(guò)?CHM?的標(biāo)準(zhǔn)差（CHMstd）反映冠層的稀疏程度以及落葉情況。郁閉度（CC）是反映林分質(zhì)量的重要參數(shù)，是冠層垂直投影與林地面積之比。使用?Li?DAR?估測(cè)郁閉度時(shí)，常用首次回波點(diǎn)與所有回波點(diǎn)的比值表示。提取的變量及其具體描述見(jiàn)表?3。

表?3?Li?DAR?變量及描述

2.4?森林健康指標(biāo)構(gòu)建

本研究采用?Person’s?相關(guān)系數(shù)度量?Li?DAR?變量之間的關(guān)系，從而確定能夠全面且系統(tǒng)地反映森林健康狀況的關(guān)鍵?Li?DAR?變量。先通過(guò)層次聚類(lèi)的類(lèi)平均法對(duì)?Li?DAR?變量進(jìn)行聚類(lèi)，把相關(guān)性相似的變量歸為一類(lèi)；再根據(jù)類(lèi)別內(nèi)相關(guān)系數(shù)最大（能代表所在類(lèi)的所有指標(biāo)）、類(lèi)別間相關(guān)系數(shù)最?。ㄅc其他類(lèi)的指標(biāo)相關(guān)性最低）的原則，在每一類(lèi)中只選取一個(gè)變量作為關(guān)鍵變量構(gòu)成森林健康指標(biāo)。為了探究森林健康指標(biāo)在解釋森林健康狀況時(shí)的機(jī)理，采用無(wú)監(jiān)督的層次聚類(lèi)法將每個(gè)指標(biāo)分為?3?類(lèi)，分別對(duì)應(yīng)刺槐健康等級(jí)中的健康、亞健康、不健康。最后根據(jù)每個(gè)指標(biāo)的聚類(lèi)結(jié)果，以類(lèi)邊界處的中值作為閾值，構(gòu)建森林健康指標(biāo)的分級(jí)閾值。與美國(guó)林業(yè)局樹(shù)冠狀況分類(lèi)指南類(lèi)似，本研究基于森林健康指標(biāo)制定了森林健康等級(jí)的判斷規(guī)則，即識(shí)別?1?個(gè)樣方是否滿足健康或不健康的規(guī)則，若不符合這?2?個(gè)等級(jí)則被視為亞健康。據(jù)此對(duì)研究區(qū)的刺槐林進(jìn)行森林健康等級(jí)判斷，采用總精度和Kappa?系數(shù)對(duì)判斷結(jié)果進(jìn)行評(píng)價(jià)。

3結(jié)果與分析

3.1?Li?DAR?變量分析

根據(jù)Li?DAR變量的Person’s?相關(guān)系數(shù)（圖?2）將其分為?4?類(lèi)。LADcv和?weibull_α?分別單獨(dú)作為一類(lèi)（與其他變量相關(guān)性低），CHMstd和?H99為一類(lèi)（r?=0.84），weibull_β、CC、VCI?為一類(lèi)，相關(guān)系數(shù)為?0.64~0.94。隨后根據(jù)選取原則，分別在后?2?個(gè)類(lèi)別中選取H99和?VCI。最終得到由?LADcv、weibull_α、H99和?VCI?4個(gè)變量組成的森林健康指標(biāo)。

圖?2?Li?DAR?變量的相關(guān)矩陣圖

?

3.3?刺槐人工林的健康判斷與評(píng)價(jià)

基于上述森林健康指標(biāo)對(duì)研究區(qū)的?30?個(gè)樣方進(jìn)行健康等級(jí)判斷，通過(guò)與野外調(diào)查結(jié)果比較，可知分類(lèi)的總精度為?86.7%，Kappa?系數(shù)為?0.79，對(duì)判斷結(jié)果的評(píng)價(jià)如表?6?所示。由表?6?可知，健康等級(jí)之間的混淆主要出現(xiàn)在健康樣方和亞健康樣方之間，健康樣方的生產(chǎn)者精度為?85.7%，說(shuō)明?14.3%?的健康樣方被錯(cuò)分為亞健康樣方，用戶精度為?92.3%，說(shuō)明7.7%?的亞健康樣方被錯(cuò)分為健康樣方；亞健康樣方的生產(chǎn)者精度為?80.0%，有?20.0%?的亞健康樣方被錯(cuò)分為健康或不健康樣方；用戶精度為?80.0%，說(shuō)明有?20.0%?的健康樣方被錯(cuò)分為亞健康樣方；重度枯梢刺槐樣方的生產(chǎn)者精度為?100.0%，而用戶精度為85.7%，說(shuō)明不健康樣方均?被?正確識(shí)別?，但是有14.3%?的亞健康樣方被錯(cuò)分為不健康。因此森林健康指標(biāo)識(shí)別健康和不健康樣方有較大優(yōu)勢(shì)，但對(duì)亞健康樣方的識(shí)別精度有待提高。

表?6?基于森林健康指標(biāo)的刺槐森林健康等級(jí)精度評(píng)價(jià)

4討論與結(jié)論?

4.1?討??論

黃河三角洲刺槐林的生長(zhǎng)主要受長(zhǎng)期土壤鹽漬化的脅迫，整體生長(zhǎng)狀況不佳，不同退化程度導(dǎo)致刺槐林的樹(shù)高、郁閉度情況存在差異。本研究利用?Li?DAR?數(shù)據(jù)從樹(shù)冠結(jié)構(gòu)角度對(duì)其健康狀況進(jìn)行判別，獲得了較理想的效果，主要是由于發(fā)生退化的刺槐多表現(xiàn)為落葉或樹(shù)木低矮，很少表現(xiàn)為樹(shù)冠枯黃變色。但是對(duì)于那些樹(shù)葉變黃但還沒(méi)有掉落的樹(shù)木，僅從結(jié)構(gòu)角度判斷仍存在缺陷，此時(shí)利用光譜信息可以檢測(cè)到這種變化。葉片變黃是因?yàn)槿~綠素含量減少，會(huì)在?525~565?nm（綠色）、690~730?nm（紅邊）和?730~760?nm（近紅外）區(qū)域發(fā)生顯著變化。森林健康是植被結(jié)構(gòu)和葉片化學(xué)特征的函數(shù)，植被結(jié)構(gòu)指樹(shù)冠的變化，而葉片化學(xué)特征是指各種色素的含量。因此，在今后的研究中可以把能描述樹(shù)冠三維結(jié)構(gòu)的?Li?DAR?技術(shù)與能反映樹(shù)葉生化特征的高光學(xué)遙感技術(shù)結(jié)合起來(lái)，將有助于提高對(duì)森林健康等級(jí)的識(shí)別精度。

4.2?結(jié)??論

本研究以黃河三角洲孤島林場(chǎng)的刺槐林為研究對(duì)象，利用無(wú)人機(jī)?Li?DAR?數(shù)據(jù)提取表征植被垂直和水平變化的?Li?DAR?特征變量，并通過(guò)聚類(lèi)分析構(gòu)建森林健康指標(biāo)，得到以下結(jié)論：①森林健康指標(biāo)包括LADcv、weibull_α、H99和?VCI，分別從樹(shù)葉茂密程度、樹(shù)冠形狀、樹(shù)高、植被垂直分布?4?個(gè)方面對(duì)森林健康狀況進(jìn)行評(píng)價(jià)；②森林健康指標(biāo)能有效識(shí)別出健康和不健康的刺槐林，對(duì)亞健康刺槐林的識(shí)別精度有待提高，總體上的分類(lèi)精度較理想（總精度為?86.7%，Kappa?系數(shù)為?0.79），說(shuō)明激光雷達(dá)技術(shù)可以在判斷森林健康狀況方面發(fā)揮重要作用。今后的研究中可與高光譜信息結(jié)合起來(lái)，構(gòu)建更為完善的森林健康指標(biāo)。

萊森光學(xué)（深圳）有限公司是一家提供光機(jī)電一體化集成解決方案的高科技公司，我們專(zhuān)注于光譜傳感和光電應(yīng)用系統(tǒng)的研發(fā)、生產(chǎn)和銷(xiāo)售。