亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于對(duì)稱凸包和平滑輪廓的單木通透度計(jì)算*

        2020-11-30 00:37:08張衛(wèi)正張偉偉李燦林萬瀚文張秋聞金保華
        林業(yè)科學(xué) 2020年10期
        關(guān)鍵詞:單木大孔樹冠

        張衛(wèi)正 張偉偉 李燦林 萬瀚文 張秋聞 劉 巖 金保華

        (1.鄭州輕工業(yè)大學(xué)計(jì)算機(jī)與通信工程學(xué)院 鄭州 450002;2.鄭州大學(xué)國際學(xué)院 鄭州 450052)

        樹冠通透度指基于樹木冠形的、透過樹木冠層可見天空區(qū)域的系數(shù),是評(píng)估森林健康和樹木健康的重要指標(biāo)之一(Svein,2009;Spiecker,1996;宋文龍等,2015)?,F(xiàn)階段,樹冠通透度缺乏估測(cè)的標(biāo)準(zhǔn)方法,在實(shí)際作業(yè)中,觀察者通常手工繪制樹木冠層的二維形態(tài)和輪廓,估計(jì)透過冠層的可見天空區(qū)域占樹冠區(qū)域的百分比,從而實(shí)現(xiàn)等級(jí)評(píng)估(Bond,2012);也有研究人員通過參考不同通透度等級(jí)的照片(卡片),人工判讀通透度的大概程度或等級(jí)(Schomakeretal.,2007)。但這些方法均依賴樹冠形態(tài)和葉量的人工視覺判斷,受觀察者視角和主觀性影響較大,具有高度的可變性,重復(fù)性差;而且,由于每株樹木的樹冠形態(tài)和葉片分布各不相同,難以建立統(tǒng)一的判別標(biāo)準(zhǔn)。

        將數(shù)字圖像處理技術(shù)用于定量評(píng)估樹冠通透度,具有快速準(zhǔn)確、自動(dòng)化程度高等優(yōu)勢(shì),能夠克服人工判別可變性大、重復(fù)性差及易受主觀影響等局限。Zhang等(2007)采用分形維數(shù)對(duì)樹冠的空間格局進(jìn)行了研究;Córica(2009)分析了陽光透過樹木冠層后在外墻上產(chǎn)生光斑的分布模式;Barthélémy等(2007)將樹冠孔隙度指標(biāo)與樹冠通透度和樹木發(fā)育階段聯(lián)系起來,但是該指標(biāo)依賴特定環(huán)境,并需要根據(jù)樹種進(jìn)行校準(zhǔn);Clark等(2003)探討樹冠通透度估算方法,以“收縮包裹”的視角定義樹冠光滑邊界,提出了廣義的冠層區(qū)域和梯度變化概念,但沒有詳細(xì)說明如何設(shè)置閾值以及閾值對(duì)樹木分枝、嫩枝等的影響。

        本研究基于樹木數(shù)字圖像進(jìn)行單木結(jié)構(gòu)特征識(shí)別與整合,不僅對(duì)樹冠所形成內(nèi)部區(qū)域的通透度進(jìn)行計(jì)算,采用K-均值聚類算法將樹冠中足夠大的可見天空區(qū)域設(shè)定為大孔(通常是由于主要樹枝的缺失造成的),將較小的可見天空區(qū)域設(shè)定為小孔(大多是由于水分脅迫、葉子病變或過早退化造成的),并計(jì)算大孔、小孔的分布密度,而且還對(duì)樹冠外部區(qū)域的通透度進(jìn)行研究,當(dāng)可見天空區(qū)域位于樹冠輪廓外部時(shí),探索可見天空區(qū)域凹陷進(jìn)入樹冠的深度,對(duì)樹冠邊界輪廓進(jìn)行平滑處理,建立樹冠的對(duì)稱凸包,計(jì)算平滑輪廓距凸包的距離和可見天空區(qū)域的凹陷深度,通過對(duì)距離設(shè)置閾值,計(jì)算深度凹陷密度?;跇涔谥械拇罂缀托】酌芏取涔谳喞纳疃劝枷菝芏?個(gè)可量化指標(biāo)計(jì)算樹冠通透度。

        本研究的創(chuàng)新之處在于:1) 將2個(gè)原始輪廓定義為參考形狀,即用于評(píng)估樹冠深度凹陷的對(duì)稱凸包,以便于計(jì)算小孔和大孔的樹冠平滑輪廓;2) 采用K-均值聚類算法,將計(jì)算得到的平滑輪廓上各點(diǎn)到凸包最近距離的集合自動(dòng)分為2類,計(jì)算類間閾值(查找小的類中的最大值與大的類中的最小值,計(jì)算二者之和的均值作為閾值),通過閾值區(qū)分輕度和深度凹陷,并統(tǒng)計(jì)樹冠中的大孔和小孔;3) 根據(jù)對(duì)樹冠通透度的貢獻(xiàn),賦予深度凹陷密度、大孔和小孔密度3個(gè)參數(shù)不同權(quán)重?;谝陨蟿?chuàng)新之處,計(jì)算單木通透度,以期為單木的健康狀況監(jiān)測(cè)和生長狀態(tài)分析提供技術(shù)支持。

        1 圖像采集與預(yù)處理

        圖像采集與預(yù)處理流程如圖1所示。

        圖1 單木通透度計(jì)算流程Fig.1 Flowchart for calculating the single tree permeability

        1.1 圖像采集 采用微軟公司開發(fā)的便攜式平板電腦Surface Pro 4獲取單木圖像。該設(shè)備觸控顯示屏尺寸為12.3英寸,分辨率為2 736×1 824像素,采用第6代Intel酷睿處理器,系統(tǒng)內(nèi)存為16 G,具有800萬像素后置自動(dòng)對(duì)焦攝像頭,配備1 024級(jí)壓感觸控筆,電池續(xù)航時(shí)間長達(dá)9 h,機(jī)身厚度8.45 mm,質(zhì)量786 g,性能強(qiáng)大又不失輕盈,攜帶方便,便于移動(dòng)操作,滿足本研究所需圖像拍攝和圖像處理要求。

        以雪松(Cedrusdeodara)為研究對(duì)象,該樹種為常綠喬木,樹冠尖塔形,大枝平展,小枝略下垂,在我國多地栽培引種,除了可用作庭園觀賞外,也是一種重要的建筑用材。采用設(shè)備自帶的高清攝像頭獲取單木圖像,圖像中可清晰辨識(shí)枯枝、新枝等反映樹木生長狀態(tài)的細(xì)節(jié)信息。在采集圖像時(shí),確保觀察者處于合適的距離,以相機(jī)可拍攝樹木全貌且圖像顯示清晰為準(zhǔn),同時(shí)盡量保持仰視角小于60°、俯視角小于30°,水平距離大于樹高,如圖2所示。標(biāo)記觀察者獲取圖像時(shí)站立的位置、樹木種類、天氣情況、拍攝日期等信息,便于后續(xù)不同時(shí)段進(jìn)行重復(fù)觀測(cè),分析和評(píng)估單木通透度的變化。

        圖2 拍攝單木圖像示意Fig.2 Schematic diagram of taking an image of single tree

        1.2 單株樹木提取 考慮到實(shí)際采集環(huán)境,目標(biāo)樹木兩側(cè)和背景中可能有建筑、車輛、相鄰粘連交叉樹木等干擾,會(huì)對(duì)單木圖像分割和提取產(chǎn)生影響,而常用的圖像分割和提取方法難以達(dá)到理想效果,故本研究利用Surface Pro 4自帶的壓感觸控筆手工圈存圖像中冠層區(qū)域,以消除圖像中的無關(guān)信息,增強(qiáng)單木冠層的可檢測(cè)性并最大限度簡化處理數(shù)據(jù)的復(fù)雜度,從而改進(jìn)特征提取、圖像分割、匹配和識(shí)別的可靠性(劉博峰,2016)。樹冠大都包含樹枝和樹葉等詮釋樹木主要形態(tài)特征和生物學(xué)特性的組件,本研究針對(duì)樹冠進(jìn)行通透度分析,排除了樹干。

        采用Matlab中的函數(shù)get(axesHandle,‘CurrentPoint’)獲取鼠標(biāo)點(diǎn)擊處的圖像坐標(biāo),并將這些坐標(biāo)組成封閉多邊形。為了提取多邊形所包圍的感興趣區(qū)域(即單木區(qū)域),使用函數(shù)inpolygon(m1,n1,im,in)判斷圖像中各像素點(diǎn)是否在多邊形內(nèi),其中im、in為構(gòu)成多邊形邊界的頂點(diǎn)坐標(biāo),m1、n1為圖像中的像素點(diǎn)坐標(biāo)。get函數(shù)返回結(jié)果為邏輯類型的0或1,如果該點(diǎn)在多邊形內(nèi)則返回1,否則為0。運(yùn)用imwrite()函數(shù)將多邊形包圍的圖像區(qū)域存儲(chǔ),區(qū)域外的圖像背景填充為黑色,可實(shí)現(xiàn)手動(dòng)點(diǎn)擊產(chǎn)生多邊形,并圈存多邊形區(qū)域,完成單木區(qū)域提取,如圖3所示。

        圖3 手工圈存目標(biāo)樹木Fig.3 Artificial circles target treea.原始圖像Original image;b.人工圈存軌跡Manually tracing trajectory;c.存儲(chǔ)樹木區(qū)域格Storage tree area grid.

        1.3 感興趣區(qū)域的灰度化和二值化 將圈存的單木區(qū)域作為感興趣區(qū)域進(jìn)行圖像灰度化和二值化,以便于后續(xù)處理。采用平均值法(馮志新等,2013)將RGB圖像轉(zhuǎn)換為灰度圖像:

        Gray=(R+G+B)/3。

        (1)

        運(yùn)用最大類間差法(Ostu算法)確定的閾值進(jìn)行圖像二值化,該算法不受圖像亮度和對(duì)比度的影響,按圖像的灰度特性,將圖像分成背景和前景2部分(廉寧等,2014)。設(shè)灰度圖像灰度級(jí)為L,則灰度范圍為[0,L-1],運(yùn)用Ostu算法計(jì)算圖像的最佳閾值為:

        σ=Max{w0(t)×[u0(t)-u]2+w1(t)×
        [u1(t)-u]2)}。

        (2)

        式中:w0為背景比例;u0為背景均值;w1為前景比例;u1為前景均值;u為整幅圖像的均值;σ取最大值時(shí)的分割閾值t即為最佳閾值T。

        使用閾值T對(duì)原圖像f(x,y)進(jìn)行分割,將圖像分為2部分,g(x,y)即為通過閾值處理得到的二值化圖像:

        (3)

        通過對(duì)圈存區(qū)域進(jìn)行圖像灰度化和二值化,得到樹木二值化圖像,如圖4a、b所示。

        對(duì)圖3c中圈存樹冠區(qū)域的綠色通道進(jìn)行二值化,得到圖4c,并將4b與c合并,最大限度保留樹冠形態(tài),得到圖4d。

        圖4 感興趣區(qū)域二值化Fig.4 Region of interest binarizationa.圖像灰度化Grayscale image;b.圖像二值化Image binarization;c.綠色通道二值化Green channel binarization;d.合并b與c Merge b and c.

        將圖4d中的圈存區(qū)域取反,得到圖5a。采用Matlab中的imclose函數(shù)對(duì)圖5a進(jìn)行閉運(yùn)算,得到圖5b,其中閉運(yùn)算采用直徑為2的圓形結(jié)構(gòu)元素D1。閉運(yùn)算對(duì)圖像先膨脹后腐蝕,具有填充物體內(nèi)細(xì)小空洞、連接鄰近物體和平滑邊界的作用。

        采用Matlab中的bwlabel函數(shù)對(duì)圖5b進(jìn)行連通區(qū)域標(biāo)記,并保留面積最大的連通區(qū)域,得到圖5c。

        圖5 單株樹木的二值化圖像Fig.5 Binarized images of single tree a.取反Negate;b.閉運(yùn)算Closed operation;c.保留面積最大的連通區(qū)域Keep the connected area with the largest area.

        1.4 樹冠平滑輪廓 首先,采用Matlab中的imfill函數(shù)對(duì)樹冠二值化圖像進(jìn)行孔洞填充,如圖6a所示。然后,采用半徑為2的圓形結(jié)構(gòu)元素對(duì)圖6a進(jìn)行閉運(yùn)算,實(shí)現(xiàn)樹冠區(qū)域的邊界平滑,運(yùn)用bwboundaries函數(shù)提取樹冠邊界點(diǎn),得到樹冠平滑輪廓,如圖6b、c所示。

        圖6 樹冠平滑輪廓Fig.6 Smooth outline of the canopy a.孔洞填充Hole filling;b.平滑輪廓Smooth outline;c.顯示平滑輪廓Display smooth outline.

        2 圖像處理

        2.1 對(duì)稱軸確定 正常生長樹木的冠形在豎直方向基本保持平衡和對(duì)稱,其對(duì)稱軸通常是樹木主干所在的直線。確定樹冠在豎直方向的對(duì)稱軸,建立樹冠基于對(duì)稱軸的鏡像,構(gòu)建對(duì)稱樹冠凸包,從而有助于后續(xù)與通透度相關(guān)參數(shù)的計(jì)算。

        針對(duì)已經(jīng)手工圈存的單木圖像,其大小為m×n(m為行數(shù),n為列數(shù))。設(shè)定對(duì)稱軸所在的列為p,則對(duì)稱軸左側(cè)的樹木像素[1~ (p-1)列的所有樹木像素]之和與對(duì)稱軸右側(cè)的樹木像素[(p+1)~n列的所有樹木像素]之和的差值最小。在Matlab中計(jì)算出p,并在圖像中畫出對(duì)稱軸,如圖7所示。

        圖7 樹冠的對(duì)稱軸Fig.7 Symmetry axis of the canopya.樹冠的對(duì)稱軸The symmetry axis of the canopy;b.在全景圖像中顯示對(duì)稱軸Show the symmetry axis in the panoramic image.

        2.2 樹冠對(duì)稱鏡像和對(duì)稱凸包確立 計(jì)算樹冠通透度時(shí),考慮樹冠中的分枝、嫩枝、新枝、葉片以及可能缺失的枝葉,以對(duì)稱軸為中心左右兩側(cè)的樹冠分別做鏡像,然后合并成為新的樹冠,如圖8b所示,白色像素為原始樹冠,灰色像素為鏡像。

        圖8 對(duì)稱樹冠與凸包Fig.8 Symmetrical crown and convex hulla.樹冠Canopy;b.對(duì)稱樹冠Symmetrical canopy;c.對(duì)稱樹冠的Delaunay三角網(wǎng)Delaunay triangulation of symmetrical canopy; d.對(duì)稱樹冠凸包Convex hull of symmetrical canopy.

        將圖像中表示對(duì)稱樹冠的每個(gè)像素當(dāng)成一個(gè)點(diǎn),每個(gè)點(diǎn)都具有行坐標(biāo)和列坐標(biāo),所有這些點(diǎn)構(gòu)成集合X。先用DelaunayTri函數(shù)得到集合X的三角網(wǎng),后用Convex hull函數(shù)得到樹冠凸包。

        采用Matlab中的DelaunayTri函數(shù)建立約束三角網(wǎng),將集合X中的所有點(diǎn)進(jìn)行規(guī)則化處理,使得集合中每個(gè)點(diǎn)位于三角形的頂點(diǎn)且三角形構(gòu)成的網(wǎng)絡(luò)具有空?qǐng)A特性和最大最小角特性(余杰等,2010;Zhanetal.,2017)。采用Matlab中的Convex hull函數(shù)獲得構(gòu)成凸包的邊界點(diǎn)(Yangetal.,2014),依次連接凸包邊界點(diǎn)構(gòu)成樹冠凸包,如圖8d所示。

        2.3 深度凹陷密度計(jì)算 為了衡量樹冠伸展情況,采用凹陷深度表示樹冠填充凸包的緊湊程度。深度凹陷密度指樹冠平滑輪廓外部可見天空區(qū)域凹陷進(jìn)入樹冠平滑輪廓較深的區(qū)域面積與凸包面積的比值。

        將樹冠平滑輪廓和對(duì)稱樹冠凸包添加到樹冠區(qū)域,如圖9a、b所示。圖9c顯示了構(gòu)成樹冠平滑輪廓和凸包的點(diǎn),并將其中一部分(黑色矩形框)放大顯示。從平滑輪廓中黑色圓點(diǎn)處開始沿平滑輪廓順時(shí)針方向行走,計(jì)算平滑輪廓上各點(diǎn)到凸包的最短距離,并畫出距離圖,如圖10a所示。

        圖9 樹冠的對(duì)稱凸包及平滑輪廓Fig.9 Symmetrical convex hull and smooth outline of canopya.平滑輪廓和凸包Smooth outline and convex hulls;b.原始圖像上顯示平滑輪廓和凸包Show smooth outline and convex hull on the original image;c.放大顯示局部的平滑輪廓和凸包Zoom in to show local smooth outline and convex hull.

        采用K-均值聚類算法(顏佩等,2017)確定深度和輕度凹陷閾值,當(dāng)凹陷深度超過此閾值即設(shè)定為深度凹陷。對(duì)距離圖進(jìn)行K-均值聚類(K=2)處理獲得閾值Ddepression,如圖10a中紅色線段所示。

        將位于樹冠平滑輪廓以外、凸包以內(nèi)的區(qū)域作為凹陷區(qū)域,如圖10b中的白色區(qū)域所示。計(jì)算該凹陷區(qū)域內(nèi)像素點(diǎn)到凸包邊界的最短距離,如果最短距離大于Ddepression,則該點(diǎn)為深度凹陷點(diǎn),所有深度凹陷點(diǎn)構(gòu)成深度凹陷區(qū)域。遍歷凹陷區(qū)域得到所有深度凹陷點(diǎn),如圖10c中灰色區(qū)域、圖10d中紅色區(qū)域所示。

        圖10 深度凹陷區(qū)域確定Fig.10 Establishment of the deep depression areaa.平滑輪廓上各點(diǎn)到凸包的最短距離The shortest distance from each point on the smooth outline to the convex hull;b.白色區(qū)域?yàn)榘枷軹he white area is the depression;c.灰色區(qū)域?yàn)樯疃劝枷軹he gray area is the depth depression;d.紅色區(qū)域?yàn)樯疃劝枷軹he red area is the depth depression.

        2.4 大小孔檢測(cè)及密度計(jì)算 將樹冠平滑輪廓內(nèi)的區(qū)域(狹義的樹冠區(qū)域)進(jìn)行連通區(qū)域標(biāo)記,采用K-均值聚類算法(設(shè)定為2類)將標(biāo)記的連通區(qū)域分為大孔和小孔,所有標(biāo)記為小孔的連通區(qū)域像素總數(shù)除以樹冠平滑輪廓所包圍的像素總數(shù)得到小孔密度,然后計(jì)算大孔密度。

        采用Matlab中的連通區(qū)域標(biāo)記函數(shù)bwlabel對(duì)樹冠中所有可見天空區(qū)域進(jìn)行標(biāo)記,統(tǒng)計(jì)其像素?cái)?shù),并對(duì)所有連通區(qū)域的像素?cái)?shù)進(jìn)行K-均值聚類,其中K為2(樹冠中所有可見天空區(qū)域分為小孔和大孔,其中小孔對(duì)應(yīng)類別1,大孔對(duì)應(yīng)類別2)。

        樹冠中樹干、樹枝和樹葉的值設(shè)為0,小孔的值設(shè)為0.333 3,大孔的值設(shè)為0.666 6,樹冠輪廓外部像素的值設(shè)為1,如圖11所示。

        圖11 樹冠的小孔和大孔Fig.11 Small and large holes of canopy

        根據(jù)所提取的樹冠、大孔和小孔,計(jì)算得大孔密度Dmacro和小孔密度Dmicro分別為0.160 5和0.179 2。

        2.5 樹冠通透度計(jì)算 參考林業(yè)專家經(jīng)驗(yàn)和樹木生長規(guī)律,采用深度凹陷密度、大孔和小孔密度3個(gè)系數(shù)定量評(píng)估單木通透度Tc:

        (4)

        深度凹陷和大孔是因樹枝缺失造成的,考慮深度凹陷對(duì)冠形和通透度的影響高于大孔,對(duì)深度凹陷賦予更大加權(quán)。此外,式(4)中也加入了小孔對(duì)通透度的貢獻(xiàn)。

        2.6 多視圖聚合 從0°、30°、60°、90°、120°和150°共6個(gè)角度獲取6個(gè)樹冠圖像,分別計(jì)算通透度系數(shù)Tc1、Tc2、Tc3、Tc4、Tc5和Tc6,以平均值Tca作為單木通透度系數(shù),減小因視角變化引起的通透度波動(dòng),盡可能精確反映單木真實(shí)狀況。如果難以從多個(gè)視角獲取樹木圖像,則根據(jù)實(shí)際獲取數(shù)量計(jì)算平均值作為通透度系數(shù):

        (5)

        3 系統(tǒng)實(shí)現(xiàn)與分析

        3.1 系統(tǒng)實(shí)現(xiàn) 由Matlab GUI (graphical user interface,為圖形用戶界面,又稱圖形用戶接口)設(shè)計(jì)開發(fā)的系統(tǒng)界面如圖12所示,依次點(diǎn)擊打開圖像、提取樹木、確立對(duì)稱軸、凸包與平滑輪廓、深度凹陷、大小孔檢測(cè)等計(jì)算單幅圖像通透度,分別顯示在右下方的表格中;也可在提取冠層后,直接點(diǎn)擊“顯示”按鈕計(jì)算通透度。

        圖12 系統(tǒng)界面Fig.12 System interface

        本研究對(duì)其他雪松進(jìn)行圖像采集和通透度分析,如圖13所示。

        圖13 實(shí)例Fig.13 Example

        3.2 樹冠對(duì)稱鏡像和對(duì)稱凸包確立 為了對(duì)本研究所提出方法進(jìn)行試驗(yàn)驗(yàn)證和精度分析,設(shè)計(jì)一個(gè)單木驗(yàn)證模型,如圖14a所示。

        圖14 單木驗(yàn)證模型Fig.14 Single tree verification model a.本研究所設(shè)計(jì)的單木驗(yàn)證模型The single-wood verification model designed by this research institute;b.人工判讀識(shí)別驗(yàn)證模型并標(biāo)示大孔、小孔及深度凹陷區(qū)域Manual interpretation to identify the verification model and mark the large holes,small holes and deep depression areas;c.將驗(yàn)證模型顯示在液晶平面顯示器 上Display the verification model on the LCD flat panel display.

        本研究所提出方法對(duì)尺度不敏感(目標(biāo)樹等比例縮放對(duì)研究結(jié)果無影響),計(jì)算出的樹冠通透度為一比值,故采用AutoCAD 2010結(jié)合Photoshop CC 2017設(shè)計(jì)單木驗(yàn)證模型(尺寸無需單位)。組成樹冠3個(gè)等腰直角三角形的腰長為4,從上到下垂直排列;樹冠中大孔是邊長為1的正方形,小孔是腰長為1的等腰直角三角形。采用人工方法確定該模型凸包,并人工計(jì)算深度凹陷閾值,畫出深度凹陷區(qū)域,如圖14b所示。計(jì)算樹冠中小孔、大孔、樹冠區(qū)域、凸包和深度凹陷區(qū)域的面積分別為3、3、24、40和6.514 8,理想情況下的小孔密度、大孔密度、深度凹陷密度和通透度分別為0.125 0、0.125 0、0.162 9和0.264 6。

        將單木驗(yàn)證模型顯示在液晶平面顯示器上,屏幕尺寸為23.8英寸,分辨率為1 920×1 080像素,垂直和水平可視角度均為178°。采用Surface Pro 4拍攝該顯示器,獲得單木模型圖像。利用本研究所提出方法計(jì)算的小孔密度、大孔密度、深度凹陷密度和通透度分別為0.117 8、0.124 1、0.164 0和0.258 6,由此得到本研究所提出方法的精度為97.73%(0.258 6 / 0.264 6)。

        3.3 深度凹陷密度計(jì)算 采用K-均值聚類算法自動(dòng)實(shí)現(xiàn)大孔和小孔分類以及深度和輕度凹陷分類。比較2個(gè)類的相似度,滿足以下條件,則認(rèn)定2個(gè)類可以合并:

        Mhigher<1.4×Mlower。

        (6)

        式中:Mhigher為較大的類均值;Mlower為較小的類均值。

        當(dāng)進(jìn)行凹陷深度計(jì)算時(shí),如果出現(xiàn)2個(gè)類合并,則認(rèn)為樹冠只呈現(xiàn)輕度凹陷,此時(shí)深度凹陷對(duì)通透度系數(shù)的貢獻(xiàn)為零。同樣,大孔和小孔均值滿足式(6),則全部顯示為小孔,此時(shí)大孔對(duì)通透度系數(shù)的貢獻(xiàn)為零。

        如果多角度采集單木圖像時(shí)背景存在黏連、遮擋等問題,導(dǎo)致難以差距6個(gè)角度的單木圖像,則式(5)中Tca應(yīng)以實(shí)際采集的圖像數(shù)量為準(zhǔn)。采用手工圈存樹冠區(qū)域方法,未能做到全自動(dòng)分析,但可減小樹冠周圍復(fù)雜背景的干擾,且人工工作量小。如果被檢測(cè)樹木的背景為其他樹木或暗色干擾物時(shí),在人工圈存樹冠區(qū)域后,可采用He等(2011)提出的基于暗通道的去霧算法進(jìn)行圖像增強(qiáng),從而提高對(duì)冠層大孔和小孔的識(shí)別效果。

        4 結(jié)論

        為解決現(xiàn)有單木通透度評(píng)估受觀察者視角和主觀性影響較大、難以建立統(tǒng)一判別標(biāo)準(zhǔn)的問題,本研究采用Surface Pro 4獲取單木圖像,進(jìn)行單木結(jié)構(gòu)特征識(shí)別與整合,通過建立樹冠對(duì)稱凸包和樹冠平滑輪廓,不僅對(duì)樹冠所形成的內(nèi)部區(qū)域進(jìn)行研究,標(biāo)記內(nèi)部的可見天空區(qū)域,并統(tǒng)計(jì)各連通區(qū)的像素?cái)?shù),采用K-均值聚類算法將各連通區(qū)域分為大孔和小孔,而且還對(duì)樹冠輪廓外部可見天空區(qū)域凹陷進(jìn)入樹冠的深度進(jìn)行探索,采用K-均值聚類算法判別輕度和深度凹陷。通過自動(dòng)量化大孔和小孔密度、深度凹陷密度3個(gè)指標(biāo)計(jì)算單木通透度系數(shù),并對(duì)過程和結(jié)果進(jìn)行詳細(xì)分析。結(jié)果表明,本研究所提出方法的精度高達(dá)97.73%,測(cè)量速度快、人工工作量小,可為單木健康狀況監(jiān)測(cè)和生長狀態(tài)分析提供技術(shù)支持,同時(shí)該研究思路和方法也可以推廣應(yīng)用到其他樹木和作物的監(jiān)測(cè)分析,具有一定的實(shí)用價(jià)值。

        猜你喜歡
        單木大孔樹冠
        地基與無人機(jī)激光雷達(dá)結(jié)合提取單木參數(shù)
        融合LiDAR點(diǎn)云與高分影像的單木檢測(cè)方法研究
        大孔ZIF-67及其超薄衍生物的光催化CO2還原研究
        樹冠羞避是什么原理?
        榕樹
        樹冠
        文學(xué)港(2019年5期)2019-05-24 14:19:42
        無人機(jī)影像匹配點(diǎn)云單木識(shí)別算法
        遙感信息(2019年1期)2019-03-22 01:38:16
        大孔鏜刀的設(shè)計(jì)
        基于雙尺度體元覆蓋密度的TLS點(diǎn)云數(shù)據(jù)單木識(shí)別算法
        森林工程(2018年5期)2018-05-14 13:54:30
        茶樹冠層光譜特征及混合像元模擬研究
        中文字幕一区二区中文| 九九99久久精品午夜剧场免费| 成在线人免费视频播放| 国产精品日韩av一区二区三区| 一本色道久久综合狠狠躁篇| 久久精品一区二区三区av| 亚洲精品国产老熟女久久| 女同另类专区精品女同| 精品久久久久久综合日本| 欧美亚洲国产片在线播放| 自拍亚洲一区欧美另类| 亚洲熟女av一区少妇| 免费无码专区毛片高潮喷水| 婷婷色中文字幕综合在线| 亚洲免费视频网站在线| 黄片免费观看视频播放| 中文字幕日韩人妻在线视频| 无码精品黑人一区二区三区| 久久久精品中文无码字幕| 粉嫩人妻91精品视色在线看| 亚洲av永久中文无码精品综合| 无码人妻丰满熟妇精品区| 男人的av天堂狠狠操| 中文字幕日韩有码在线| 国产美女在线精品免费观看| 无码中文av有码中文av| 精品久久人妻一区二区| 十四以下岁毛片带血a级| 欧美性猛交xxxx乱大交丰满| 亚洲性色ai无码| 亚洲成人福利在线视频| 亚洲老妈激情一区二区三区 | 欧美午夜刺激影院| 狠狠久久精品中文字幕无码| 国产91大片在线观看| 成人艳情一二三区| 婷婷五月综合缴情在线视频| 熟女系列丰满熟妇av| 午夜福利影院成人影院| 欧洲熟妇色xxxx欧美老妇多毛图片| 甲状腺囊实性结节三级|