像素當(dāng)量法在立木材積測(cè)算中的應(yīng)用研究

摘要：為了研究像素當(dāng)量法在立木材積測(cè)算中的應(yīng)用，根據(jù)智能手機(jī)拍攝的立木像片，采用像素當(dāng)量法并運(yùn)用區(qū)分段不等長(zhǎng)的中央斷面區(qū)分求積式測(cè)算了銀杏立木材積，并分析了不同攝影方位和不同區(qū)分段數(shù)對(duì)立木材積精度的影響。結(jié)果顯示，隨著觀測(cè)方位數(shù)量的增加，立木材積相對(duì)誤差的絕對(duì)值在減小，其中，根據(jù)兩個(gè)互相垂直方位拍攝的像片測(cè)算的立木材積相對(duì)誤差小于1.5%；同時(shí)，隨著區(qū)分段數(shù)的增加，立木材積相對(duì)誤差的絕對(duì)值在減小。當(dāng)區(qū)分段在5～8 段時(shí)，立木材積相對(duì)誤差絕對(duì)值均在1.4%以內(nèi)。這就表明在用像素當(dāng)量法精確測(cè)算立木材積時(shí)，至少需要采用兩個(gè)互相垂直的方位拍攝立木像片，且區(qū)分段數(shù)定為5～8段左右為宜。

關(guān)鍵詞：像素當(dāng)量 法立木中央斷面區(qū)分求積式材積

中圖分類號(hào)：P232

Research on the Application of the Pixel Equivalent Method in the Measurement ofing Standing Volume

ZHAO HaoyanZHANG Minxia*ZHANG Lei

Nanjing Police CollegeUniversity， Nanjing， Jiangsu Province，210023 China

Abstract： In order to research onstudythe application of the pixel equivalent method in themeasuresureingof standing volume，the standing volume of Gginkgo trees biloba L.wais measured usingby the pixel equivalent method and the sectional measurementvolume calculation Gjo2tiEw5APAxucnAJrNkA==formula of central cross-sectionals with the unequal length of segments based onaccording totheimagephotos of standing trees capturtakened by smartphones， and the aeffects of the number ofdifferent photographicy orientations and the number ofdifferent sectionsegments on the accuracy of standing volume weare analyzed.The results reveal that the absolute value of the relative error of standing volume decreases with the increase of the number of photographicy orientations， andamong which， the relative error of the standing volume measured byphotos taken from the two mutually-perpendicularphotography orientations is less than 1.5%. Simultaneously， and thatthe absolute value of the relative error of standing volumedecreases with the increase of the number of sections， andthe absolute value of the relative error of standing volume is within 1.4% when the number of sectionsegments is 5-～8. Therefore，which shows thatthe photos of standing trees were taken in it isnecessary to use at least two mutually-perpendicularorientationsvertical directionsto takethe photos of standing treesand the number of sectionsegments isshould be at 5-～8 towhen accurately measuringe thestanding volume withbythe pixel equivalent method.

Key Words：The pPixel equivalent method; Standing tree; Sectional measurementThe volume calculation formula of central cross-sectionals; Volume

傳統(tǒng)測(cè)算立木材積的方法有形數(shù)法、形率法、一元材積表法和二元材積表法等。但是，傳統(tǒng)測(cè)算立木材積的方法存在勞動(dòng)強(qiáng)度大、效率低等問題。為了提高測(cè)量效率，一些學(xué)者探索采用各種先進(jìn)的工具來(lái)估測(cè)立木材積。練一寧等人[1]、郁壯等人[2]采用近景攝影測(cè)量中的共線性方程，解算出立木的胸徑和樹高，并將立木主干分為圓錐、臺(tái)柱和圓柱3種幾何體，從而計(jì)算材積；王俞明[3]根據(jù)地基激光雷達(dá)獲取的杉木的點(diǎn)云數(shù)據(jù)，構(gòu)建了樹干的Delaunay三角網(wǎng)模型并提取了林木材積，取得了較高的精度；潘黃儒等人[4]根據(jù)地基激光雷達(dá)獲取的桉樹的點(diǎn)云數(shù)據(jù)提取了材積，研究發(fā)現(xiàn)提取的材積與二元材積表計(jì)算的理論材積有很強(qiáng)的回歸關(guān)系；吳伊[5]采用隨機(jī)森林方法擬合了欒樹的一元和二元材積模型，并取得了很好的效果。雖然有學(xué)者采用像素當(dāng)量法，并運(yùn)用區(qū)分求積式估算了立木材積。但是，目前還沒有學(xué)者研究采用像素當(dāng)量法，至少測(cè)量幾個(gè)方位且劃分多少區(qū)分段數(shù)可較為精確地測(cè)算立木材積。而且以上測(cè)算材積用的相機(jī)、三維激光掃描儀等儀器價(jià)格昂貴，智能手機(jī)攝像頭具有像素高、價(jià)格低、使用廣泛等優(yōu)點(diǎn)，所以，本文探討基于手機(jī)攝影，采用像素當(dāng)量法精確測(cè)算立木材積的方法。

1研究對(duì)象與方法

1.1研究對(duì)象

為了確定采用像素當(dāng)量法，至少測(cè)量幾個(gè)方位且劃分多少區(qū)分段數(shù)可較為精確地測(cè)算立木材積，本論文選擇南京警察學(xué)院內(nèi)的銀杏（Ginkgo biloba L.）立木作為研究對(duì)象，這是因?yàn)檫@些銀杏樹的樹干比較通直、分杈少、枝葉少，便于測(cè)量不同部位的樹干直徑。校園中有27株銀杏樹，胸徑為（17.8±3.1）cm，樹高為（10.25±3.00）m。

1.2 研究方法

1.2.1調(diào)查方法

在27株銀杏樹中，選擇1株從不同方位觀測(cè)視線障礙物均最少的銀杏樹（樹高為10.65 m）作為目標(biāo)樹，在樹旁豎向放置1把長(zhǎng)度為50 cm的直尺作為參照物；隨后打開HUAWEI Mate50pro智能手機(jī)安裝的“手機(jī)測(cè)量水平儀”App，使手機(jī)保持與地面垂直；最終使用智能手機(jī)攝像頭（等效35mm膠片焦距f=27mm），采用“專業(yè)”模式，保持手機(jī)焦距不變，離目標(biāo)樹水平距離14.9m（該距離剛好能拍攝到整株目標(biāo)樹），分別從正南、東南、正東、東北、正北、西北、正西、西南8個(gè)方位，水平拍攝目標(biāo)樹。

1.2.2像素當(dāng)量法

根據(jù)直尺的實(shí)際長(zhǎng)度（50 cm）與成像像素?cái)?shù)之比計(jì)算像素當(dāng)量（單位為mm/ pixel），立木直徑的成像像素?cái)?shù)乘以像素當(dāng)量即為實(shí)際直徑。通過將手機(jī)拍攝的立木像片輸入“Adobe Photoshop CC 2018”圖像處理軟件，依次點(diǎn)擊“圖像—分析—標(biāo)尺工具”，可測(cè)量直尺成像的像素值和銀杏立木各區(qū)分段中央直徑的像素值，從而根據(jù)像素當(dāng)量計(jì)算立木各部位直徑。

1.2.3基于不等長(zhǎng)區(qū)分段的中央斷面區(qū)分求積法

通常采用等長(zhǎng)區(qū)分段的中央斷面區(qū)分求積法精確測(cè)算倒木材積，在用該方法測(cè)算立木材積時(shí)，經(jīng)常出現(xiàn)測(cè)量位置被枝葉遮擋的情況，從而影響測(cè)算立木材積的精度。因此，本文采用不等長(zhǎng)區(qū)分段的中央斷面區(qū)分求積式測(cè)算立木材積。具體原理圖如圖1所示。

根據(jù)積分的思想，將立木樹干分為n+1個(gè)區(qū)分段，其中，前n個(gè)區(qū)分段均近似看作拋物線體，最后一段梢頭近似看作圓錐體，所有區(qū)分段材積之和，即為立木材積（圖1）。劃分的區(qū)分段數(shù)越多，計(jì)算的林木材積越精確。因?yàn)殡x樹干基部越近，一般樹干直徑越大，區(qū)分段材積占整株立木材積的比重越大，所以，為了提高測(cè)量精度，離樹干基部越近，劃分的區(qū)分段數(shù)越多，反之，亦然。具體公式如下 。

v=∑_（i=1）^n?〖v_i+v^' 〗

=v_1+v_2 〖+v〗_3+……〖+v〗_n+v^'

=l_1 g_1+l_2 g_2+l_3 g_3+……〖+l_n g〗_n+1/3 g^' l^'

=π/40000（l_1 d_1^2+l_2 d_2^2+l_3 d_3^2+……+l_n d_n^2+1/3 l^' d_^'2）（1）

式中，v1，v2，v3，…vn分別為立木第1、第2、第3、…、第n個(gè)區(qū)分段的材積，v^'為梢頭材積，l1，l2，l3，…ln分別為第1、第2、第3…第n個(gè)區(qū)分段的長(zhǎng)度，g1，g2，g3，…，gn分別為第1、第2、第3、…、第n個(gè)區(qū)分段的中央段面積，g^'為立木梢頭的底面積，l^'為梢頭長(zhǎng)度，d_^'為梢頭底直徑；d1，d2，d3，…，dn分別為第1、第2、第3 、…、第n個(gè)區(qū)分段的中央直徑。

2研究結(jié)果

2.1 不同攝影方位對(duì)測(cè)算立木材積的影響

根據(jù)手機(jī)分別從正南、東南、正東、東北、正北、西北、正西、西南8個(gè)方位拍攝的立木像片，分別采用像素當(dāng)量法，將樹干劃分12個(gè)區(qū)分段，保持八個(gè)方位每個(gè)區(qū)分段中央直徑的位置高度一致，并運(yùn)用式（1）計(jì)算立木材積。

從表1可看出，根據(jù)8個(gè)方位拍攝的手機(jī)像片測(cè)算的立木材積有一定差異，正西和東北兩個(gè)方位測(cè)算的立木材積較大，正南和西南方位測(cè)算的立木材積較小，這主要是因?yàn)闃涓蓴嗝嫘螤畈灰?guī)則，正西和東北兩個(gè)方位劃分的區(qū)分段的中央直徑整體偏大，正南和西南方位的大多數(shù)區(qū)分段的中央直徑偏小所致。

以8個(gè)方位測(cè)算的立木材積的平均值作為真值，分析手機(jī)拍攝立木的方位數(shù)量對(duì)立木材積相對(duì)誤差的影響（如表2所示）。從圖2可以看出，隨著觀測(cè)方位數(shù)量的增加，立木材積相對(duì)誤差的絕對(duì)值在減小。這是因?yàn)榱⒛緲涓蓴嗝娴男螤畲蠖喑什灰?guī)則形狀，如果根據(jù)1個(gè)方位拍攝的像片，通過量取樹干各區(qū)分段的中央直徑計(jì)算的材積會(huì)有較大的誤差。所以，拍攝的方位越多，計(jì)算的立木材積越精確。但是，測(cè)算材積的工作量也會(huì)大大增加。當(dāng)采用兩個(gè)互相垂直的方位拍攝的像片計(jì)算立木材積時(shí)，材積的相對(duì)誤差在1.5%以內(nèi)，不僅滿足林業(yè)測(cè)量精度，而且測(cè)算材積的工作量也較小。

2.2 不同區(qū)分段數(shù)對(duì)測(cè)算立木材積的影響

根據(jù)正北和正東兩個(gè)方位拍攝的手機(jī)像片，分別將樹干劃分20個(gè)區(qū)分段，保持這兩個(gè)方位所有區(qū)分段中央直徑的位置高度一致，采用式（1）計(jì)算立木材積。以這兩個(gè)方位測(cè)算的平均材積作為真值，分析不同區(qū)分段數(shù)對(duì)立木材積相對(duì)誤差的影響（如圖3所示）。從圖3可以看出，隨著區(qū)分段數(shù)的增加，立木材積相對(duì)誤差的絕對(duì)值在減小。當(dāng)區(qū)分段在 1 ～5段時(shí)，立木材積相對(duì)誤差絕對(duì)值的變化幅度明顯，隨著區(qū)分段個(gè)數(shù)的增加而急劇減?。划?dāng)區(qū)分段在5～8段時(shí)，立木材積相對(duì)誤差絕對(duì)值均在1.4%以內(nèi)，且隨著區(qū)分段個(gè)數(shù)的增加，立木材積相對(duì)誤差的絕對(duì)值在緩慢減??；當(dāng)區(qū)分段在8～17段時(shí)，立木材積相對(duì)誤差的絕對(duì)值隨區(qū)分段數(shù)的增加呈波動(dòng)的趨勢(shì)，這可能是因?yàn)椴糠謪^(qū)分段的中央斷面形狀不規(guī)則，從而導(dǎo)致測(cè)算的立木材積誤差較大所致；當(dāng)區(qū)分段在18～20段時(shí)，曲線接近水平，表明隨著區(qū)分段個(gè)數(shù)的增加，立木材積相對(duì)誤差的絕對(duì)值變化幅度不明顯。

3結(jié)論和討論

本文主要探討了根據(jù)手機(jī)拍攝的立木像片，采用像素當(dāng)量法并運(yùn)用區(qū)分段不等長(zhǎng)的中央斷面區(qū)分求積式測(cè)算立木材積的方法。主要結(jié)論如下 。

第一，將立木主干視為拋物線體，并用中央斷面近似求積式計(jì)算其材積，這種方法計(jì)算的材積相對(duì)誤差很大，達(dá)到-31.4%，當(dāng)區(qū)分段數(shù)達(dá)到2個(gè)及以上，采用區(qū)分求積式計(jì)算立木材積的精度明顯提高（圖3） 。

第二，為了提高測(cè)量精度，本文提出了從樹干基部到最高點(diǎn)區(qū)分段數(shù)逐漸減少的區(qū)分求積方法。與區(qū)分段等長(zhǎng)的中央斷面區(qū)分求積式相比，本區(qū)分段不等長(zhǎng)的區(qū)分求積方法更為靈活，可以避開測(cè)量位置有枝節(jié)或被枝葉遮擋的情況，從而具有更高的材積測(cè)算精度。

第三，隨著觀測(cè)方位數(shù)量的增加，立木材積相對(duì)誤差的絕對(duì)值在減?。▓D2）。采用兩個(gè)互相垂直的方位拍攝的像片計(jì)算的立木材積的相對(duì)誤差在1.5%以內(nèi)，同時(shí)測(cè)算材積的工作量也較小。因此，在用像素當(dāng)量法精確測(cè)算立木材積時(shí)，至少需要采用兩個(gè)互相垂直的方位拍攝立木像片 。

第四，隨著區(qū)分段數(shù)的增加，立木材積相對(duì)誤差的絕對(duì)值在減?。▓D3）。當(dāng)區(qū)分段在5～ 8段時(shí)，立木材積相對(duì)誤差絕對(duì)值均在1.4%以內(nèi)。因此，測(cè)量一株高度為10m的立木時(shí)，從精度和工作量?jī)蓚€(gè)方面考慮，可以將區(qū)分段數(shù)定為5～8段左右。

參考文獻(xiàn)

[1] 練一寧，馮仲科，劉培斌，等. 智能機(jī)單片攝影測(cè)樹技術(shù)研究與試驗(yàn)[J].中國(guó)農(nóng)業(yè)科技導(dǎo)報(bào)， 2020， 22 （3）： 56-63.

[2] 郁壯，馮仲科，張標(biāo)，等.無(wú)損方法實(shí)現(xiàn)重慶地區(qū)馬尾松材積模型研建[J].中國(guó)農(nóng)業(yè)科技導(dǎo)報(bào)， 2023，25（7）：97-104.

[3] 王俞明.基于地基激光雷達(dá)的杉木參數(shù)提取與材積估測(cè)[D].長(zhǎng)沙：中南林業(yè)科技大學(xué)， 2019.

[4] 潘黃儒，孟想，張振林，等.基于TLS點(diǎn)云數(shù)據(jù)的桉樹二元材積表檢驗(yàn)[J].廣西林業(yè)科學(xué)， 2023，52（4）：516-521.

[5] 吳伊.魯中地區(qū)欒樹材積表與碳測(cè)算表的編制[D].泰安：山東農(nóng)業(yè)大學(xué)，2022.