摘 要：針對(duì)含節(jié)子木材力學(xué)特性不確定，不易判斷其是否可用的現(xiàn)狀，提出一種通過(guò)檢測(cè)含節(jié)子木材的抗彎性能來(lái)對(duì)含有節(jié)子的木材是否可用進(jìn)行評(píng)判的方法。選取在東北地區(qū)占到總森林面積15%～20%的常見樹木蒙古櫟為試驗(yàn)對(duì)象，首先采用目標(biāo)檢測(cè)算法對(duì)木材表面含節(jié)子區(qū)域進(jìn)行識(shí)別；然后對(duì)識(shí)別的區(qū)域進(jìn)行光譜提取，并構(gòu)建定量預(yù)測(cè)模型；最后通過(guò)深度學(xué)習(xí)對(duì)含節(jié)子木材的力學(xué)性能進(jìn)行分析。試驗(yàn)結(jié)果表明，提出的基于連續(xù)投影的支持向量機(jī)算法（Successive Projections Algorithm-Support Vector Machine，SPA-SVM）預(yù)測(cè)模型對(duì)木材抗彎性能具有優(yōu)秀的預(yù)測(cè)能力，其試驗(yàn)結(jié)果指標(biāo)決定系數(shù)R2=0.96，均方根誤差RMSE=0.58，相對(duì)分析誤差RPD=5.09。該預(yù)測(cè)模型能非常準(zhǔn)確地對(duì)含節(jié)子木材的抗彎性能進(jìn)行預(yù)測(cè)，預(yù)測(cè)結(jié)果與真實(shí)數(shù)值誤差較小，符合試驗(yàn)要求標(biāo)準(zhǔn)，預(yù)測(cè)結(jié)果可以為木材是否使用提供依據(jù)。

關(guān)鍵詞：近紅外光譜； 抗彎性能； 力學(xué)樣本制備； 光譜數(shù)據(jù)處理； 特征值提取； 支持向量機(jī)； 預(yù)測(cè)模型

中圖分類號(hào)：S784 文獻(xiàn)標(biāo)識(shí)碼：A DOI：10.7525/j.issn.1006-8023.2024.04.018

Research on the Bending Resistance of Wood Containing Knots Based on Near-infrared Spectroscopy Technology

SONG Chuang1， SUN Liping2， WANG Pengkun1， YANG Yang3*

（1.The 723 Institute of CSSC， Yangzhou 225000， China；2.College of Mechanical and Electrical Engineering， Northeast Forestry University， Harbin 150040， China；3.Harbin University， Harbin 150076， China）

Abstract： In view of the uncertain mechanical properties of knot-bearing wood and the difficulty of judging whether it is usable， this article proposes a method to evaluate the usability of wood containing knots by detecting the bending properties of wood containing knots. The common tree Quercus mongolica， which accounts for 15%-20% of the total forest area in Northeast China， is selected as the experimental object. Firstly， the object detection algorithm is used to identify the areas containing knots on the suface of wood， followed by spectral extraction of the identified area and the construction of a quantitative prediction model. Finally， the mechanical properties of wood containing knots are analyzed through deep learning. The experimental results indicate that the SPA-SVM prediction model proposed in this article has excellent predictive ability for the bending properties of wood， with experimental results indicators of R2=0.96， RMSE=0.58， and RPD=5.09. The prediction model proposed in this article can accurately predict the bending properties of wood containing knots. The predicted results have a small error with the actual values， which meets the experimental requirements and standards. The predicted results can provide a basis for whether the wood can be used.

Keywords： Near infrared spectroscopy； bending properties； machine learning； predictive model

0 引言

國(guó)家林業(yè)和草原局發(fā)布的第九次森林資源調(diào)查報(bào)告表明，中國(guó)的森林總面積約為2.20億hm2，森林總面積和木材蓄積量分別位于世界第5位和第 6位［1］，但與之相對(duì)的是我國(guó)人均森林面積僅為世界人均的1/4。國(guó)內(nèi)木材供給缺口的增大導(dǎo)致木材進(jìn)口數(shù)量近年持續(xù)增加，通過(guò)對(duì)海關(guān)2020年的數(shù)據(jù)進(jìn)行分析可以得出，我國(guó)木材僅2020一年進(jìn)口量就高達(dá)1.075億m3（原木體積），40%以上的木材沒(méi)有供應(yīng)保障。盡管如此，我國(guó)木材的平均利用率只有65%，造成了極大的浪費(fèi)，而且在被浪費(fèi)的木材中有60%左右是含有節(jié)子的木材［2］。在木材加工中，通常被廢棄的含有節(jié)子的木材有2種原因：除去因?yàn)槊烙^選用無(wú)瑕疵的木材，更重要的是擔(dān)心含節(jié)子木材力學(xué)性能不能滿足產(chǎn)品要求從而導(dǎo)致產(chǎn)品質(zhì)量低下。為了提高整體的利用率，要提高含有節(jié)子木材的利用，但在以往的研究中對(duì)于含有節(jié)子的木材的研究較少。

Horvath等［3］利用近紅外光譜技術(shù)建立了針對(duì)山楊木樹種抗彎強(qiáng)度和彈性模量的偏最小二乘預(yù)測(cè)模型；Schimleck等［4］應(yīng)用BP神經(jīng)網(wǎng)絡(luò)，通過(guò)近紅外光譜數(shù)據(jù)對(duì)無(wú)瑕疵蒙達(dá)利松的彈性模量（Modulus of elasticity，MOE）進(jìn)行了預(yù)測(cè)；Kothiyal等［5］在不對(duì)木材含水率進(jìn)行統(tǒng)一處理的情況下，對(duì)細(xì)葉桉的近紅外光譜數(shù)據(jù)與力學(xué)性能之間的關(guān)系進(jìn)行了研究，發(fā)現(xiàn)在此條件下隨機(jī)森林模型的測(cè)定系數(shù)較高；Xu等［6］研究發(fā)現(xiàn)了木材力學(xué)性能試件中水分含量的不同會(huì)影響靜曲強(qiáng)度（Modulus of Rapture，MOR）預(yù)測(cè)模型的精確度；Liang等［7］通過(guò)基于遺傳算法優(yōu)化的偏最小二乘法進(jìn)行光譜優(yōu)選，構(gòu)建了蒙古櫟的近紅外光譜值與其MOR對(duì)應(yīng)關(guān)系的模型。

對(duì)木材力學(xué)性能評(píng)估的常規(guī)方法是從待測(cè)試樣中隨機(jī)取樣進(jìn)行力學(xué)破壞性試驗(yàn)，雖然該方法的測(cè)定結(jié)果較為準(zhǔn)確，但破壞性試驗(yàn)后的試樣因發(fā)生不可逆形變因而無(wú)法繼續(xù)使，容易導(dǎo)致資源的浪費(fèi)［8］。本研究將通過(guò)目標(biāo)檢測(cè)與近紅外光譜技術(shù)結(jié)合的方法，對(duì)東北地區(qū)的一種次生林樹種蒙古櫟進(jìn)行測(cè)試，采用力學(xué)性能最重要的指標(biāo)MOR作為評(píng)判指標(biāo)［9］。

1 材料與方法

1.1 研究方法

首先，對(duì)試樣進(jìn)行采集多光譜圖像數(shù)據(jù)并測(cè)定其木材力學(xué)性能的真實(shí)值；其次，根據(jù)試樣表面的圖像標(biāo)定出節(jié)子面積與所占位置，并以此數(shù)據(jù)為依據(jù)求出節(jié)子部分和無(wú)瑕疵部分的慣性矩與整個(gè)試樣橫截面的慣性矩比，該比值作為隨機(jī)采樣后節(jié)子部分和無(wú)瑕疵部分光譜數(shù)據(jù)處理時(shí)的加權(quán)比，以此得到代表試樣的加權(quán)平均光譜值，接下來(lái)剔除異常試樣，劃分訓(xùn)練集和測(cè)試集的試樣；然后，對(duì)原始光譜數(shù)據(jù)進(jìn)行預(yù)處理和特征提?。蛔詈?，通過(guò)處理好的訓(xùn)練集數(shù)據(jù)構(gòu)建預(yù)測(cè)模型。

1.2 試樣制備與數(shù)據(jù)采集

1.2.1 含節(jié)子木材力學(xué)試樣的制備

本試驗(yàn)試樣所選擇的樣品木材為蒙古櫟（Quercus mongolica），是東北地區(qū)常見的樹種，采樣地點(diǎn)為黑龍江省五常市林業(yè)局的沖河林場(chǎng)。在試驗(yàn)區(qū)域內(nèi)按照海拔由高到低的順序隨機(jī)砍伐12株樣木，并標(biāo)出每株樣木的生長(zhǎng)朝向，在距地面約1.3 m處向上切割木材。通過(guò)切割獲得長(zhǎng)度為1 m的第1段木材，然后以2 m的間隔切割獲得長(zhǎng)度為1 m的第2段木材，重復(fù)循環(huán)，直到接近樹冠。

根據(jù)國(guó)家《無(wú)疵小試樣木材物理力學(xué)性質(zhì)試驗(yàn)方法第2部分：取樣方法和一般要求》（GB/T 1927.2—2021）要求，將木材切片制成粗條，選取粗條中含節(jié)子的部分，制備尺寸為300 mm×70 mm×70 mm的含節(jié)子木材力學(xué)試樣共計(jì)150個(gè)。對(duì)這些試樣按照編號(hào)QX1—QX150的順序進(jìn)行標(biāo)記，之后將所有試樣放入干燥箱中干燥，以確保試樣水分含量均為12%左右，最后將每個(gè)試樣分別放入單個(gè)密封袋中，以防止存放時(shí)試驗(yàn)室內(nèi)的水蒸氣對(duì)試樣造成影響。

1.2.2 可見/近紅外多光譜圖像采集

為了保證試驗(yàn)的準(zhǔn)確性，數(shù)據(jù)在室內(nèi)溫度保持20 ℃左右，室內(nèi)濕度保持在50%左右的試驗(yàn)室進(jìn)行采集。使用芬蘭SPECIM公司所生產(chǎn)的FX-10結(jié)合ENVI自帶的采點(diǎn)功能，在節(jié)子部分與無(wú)瑕部分各隨機(jī)采集20個(gè)點(diǎn)的光譜數(shù)據(jù)，將采集得到的ASCII碼轉(zhuǎn)出為txt文件方便接下來(lái)的處理。

1.2.3 含節(jié)子試樣力學(xué)真值測(cè)量

試驗(yàn)按照《無(wú)疵小試樣木材物理力學(xué)性質(zhì)試驗(yàn)方法第9部分：抗彎強(qiáng)度測(cè)定》（GB/T 1927.9—2021）《無(wú)疵小試樣木材物理力學(xué)性質(zhì)試驗(yàn)方法第10部分：抗彎彈性模量測(cè)定》（GB/T 1927.10—2021）中規(guī)定的試驗(yàn)流程和規(guī)范設(shè)計(jì)，通過(guò)使用萬(wàn)能力學(xué)試驗(yàn)機(jī)的三點(diǎn)彎曲法對(duì)含節(jié)子木材試樣進(jìn)行力學(xué)破壞性試驗(yàn)，以此來(lái)測(cè)量各試樣的靜曲強(qiáng)度（MOR）和彈性模量（MOE），具體試驗(yàn)示意圖如圖1所示。圖1中：P表示最大荷載，N；b與h分別表示試樣寬度與試樣高度，mm；L表示兩支座間跨度距離，mm。

試樣含水率為12%時(shí)，按式（1）計(jì)算MOR，按式（2）計(jì)算MOE。

σ=PL/（bh^2 ）。 （1）

E=（23PL^3）/108vbh。 （2）

式中：σ表示抗彎強(qiáng)度，MPa；E表示彈性模量，MPa；v表示試樣在經(jīng)受上、下限荷載時(shí)的變形值，mm。

1.2.4 木材節(jié)子面積測(cè)量

使用目標(biāo)檢測(cè)的方法來(lái)獲取試樣表面節(jié)子的面積，節(jié)子區(qū)域在木材上可以視作小目標(biāo)，所以可以使用小目標(biāo)檢測(cè)法進(jìn)行對(duì)節(jié)子的識(shí)別［10］。常見的小目標(biāo)檢測(cè)算法主要為一階段算法，本研究采用一階段算法中最具代表性的算法—YOLO V5s作為檢測(cè)節(jié)子位置與大小的算法并將得到的節(jié)子橫截面積用于下一步的慣性矩計(jì)算中。

1.2.5 慣性矩比值計(jì)算

在之前的相關(guān)研究中，Lisbeth等［11］針對(duì)含節(jié)子木材研究發(fā)現(xiàn)，當(dāng)選取密度和最大凹陷處距離相同的含節(jié)子木材時(shí)，用木節(jié)的慣性矩與該試樣整體的慣性矩對(duì)其足尺M(jìn)OR進(jìn)行多元回歸分析，木節(jié)慣性矩與其對(duì)應(yīng)試樣橫截面慣性矩的比值、試樣的平均密度和試樣整體的力學(xué)性能有一定的比例關(guān)系，具體如式（3）所示。

f_f=aI_k/I_g+bρ^1.25+c。 （3）

式中：f_f為足尺抗彎強(qiáng)度；I_k/I_g為木節(jié)的慣性矩（Ik）與整個(gè)橫截面的慣性矩（Ig）比；ρ為木材平均密度；a、b與c均為擬合參數(shù)。

本試驗(yàn)所選用的試樣均為同一林場(chǎng)中采伐的蒙古櫟，因此可將各試樣的平均密度看作相同值，由此可推知木節(jié)慣性矩與試樣橫截面慣性矩的比值和該試樣的MOR值有近似線性關(guān)系。同時(shí)因?yàn)楸驹囼?yàn)是通過(guò)對(duì)試樣進(jìn)行縱向物理破壞來(lái)測(cè)得各試樣對(duì)應(yīng)MOR值，因此選用式（4）來(lái)計(jì)算節(jié)子與試樣橫截面對(duì)X軸慣性矩的值。

I_x=∫_A?y^2 dA。 （4）

式中：A表示所求慣性矩對(duì)應(yīng)區(qū)域的面積；y表示試樣的高度［12］。

舉例：某塊試樣經(jīng)過(guò)計(jì)算后可以得到節(jié)子區(qū)域的慣性矩與試樣橫截面的慣性矩的比值Ik/Ig為0.119 1，即在光譜值加權(quán)平均計(jì)算時(shí)，節(jié)子部分所對(duì)應(yīng)的光譜數(shù)據(jù)占所有采集的光譜數(shù)據(jù)的比值為11.91%。

在光譜值采集階段，每塊試樣均隨機(jī)獲取了 20個(gè)節(jié)子點(diǎn)的光譜值與20個(gè)無(wú)瑕區(qū)域的光譜值，將Ik/Ig作為光譜數(shù)據(jù)計(jì)算的加權(quán)比，以此把每一份木材試樣對(duì)應(yīng)的平均光譜值求出作為原始光譜數(shù)據(jù)。

1.3 數(shù)據(jù)處理與預(yù)測(cè)模型建立

1.3.1 光譜數(shù)據(jù)的處理

多光譜鏡頭采集并加權(quán)平均處理的光譜數(shù)據(jù)不能直接用于算法訓(xùn)練，需要進(jìn)行預(yù)處理，才能提取特征值，進(jìn)行訓(xùn)練。預(yù)處理需要分為3步：去噪、去異常、求導(dǎo)。

首先，使用SG（Savitzky-Golay）卷積平滑去除光譜數(shù)據(jù)中的高頻噪聲，保證光譜數(shù)據(jù)圖像的平滑［13］。經(jīng)過(guò)擬合后的光譜數(shù)據(jù)與未處理的原始光譜數(shù)據(jù)，通過(guò)式（5）可以求得其之間的誤差（ε）。

ε=∑_（-m）^m?（X ?_λ^（（j） ）-X_λ^（（j） ） ）^2 。 （5）

式中：-m到m為窗口寬度；X_λ^（（j） ）為在對(duì)應(yīng)窗口寬度下的反射率；X ?_λ^（（j） ）為所有窗口下反射率的加權(quán)平均值。通過(guò)最小二乘法可以求出當(dāng)ε取得最小值時(shí)所對(duì)應(yīng)的各項(xiàng)加權(quán)系數(shù)。

然后，將去噪后的光譜數(shù)據(jù)通過(guò)使用馬氏距離法（Mahalanobis Distance，MD，式中記為MD）將異常的光譜值去掉［14］，其公式為

M_D=√（（x_i-x_j ） C_x （x_i-x_j ）^T ）。 （6）

C_x=（（X_c ）^T （X_c ））/（（n-1） ）。 （7）

式中：X_c表示中心化后的光譜矩陣；C_x表示協(xié)方差矩陣；x_i與x_j表示數(shù)據(jù)集中任意2條光譜數(shù)據(jù)；n表示數(shù)據(jù)集中所含試樣的個(gè)數(shù)。

最后，通過(guò)導(dǎo)數(shù)處理法來(lái)去除基線漂移對(duì)光譜數(shù)據(jù)的影響，并提供相較原始數(shù)據(jù)更加明顯的反射峰［15］。由于本試驗(yàn)的數(shù)據(jù)是離散的，所以可以直接通過(guò)差分法求出各個(gè)試樣對(duì)應(yīng)的1stDer光譜數(shù)據(jù)，公式為

dx/dλ=（x_（j+1）-x_j）/?λ。 （8）

式中：x表示反射率；j表示波長(zhǎng)點(diǎn)數(shù)所對(duì)應(yīng)的序號(hào)；?λ表示波長(zhǎng)間距。

1.3.2 特征值提取

本試驗(yàn)采用連續(xù)投影算法（SPA）對(duì)處理過(guò)后的光譜數(shù)據(jù)進(jìn)行特征值提?。?6］，因?yàn)镾PA算法通過(guò)提取光譜數(shù)據(jù)中最小共線性和數(shù)據(jù)重復(fù)性最小的特征數(shù)據(jù)，減少光譜信息的重疊，以此來(lái)達(dá)到最大程度地用幾個(gè)光譜特征數(shù)據(jù)代替被測(cè)樣品的完整光譜信息的目的［4，17-18］。

1.3.3 評(píng)價(jià)指標(biāo)

對(duì)比不同模型，需要使用相同的衡量標(biāo)準(zhǔn)以確保有明確的數(shù)據(jù)說(shuō)明模型預(yù)測(cè)性能，所以本研究采用多個(gè)衡量指標(biāo)來(lái)判斷各模型的差異，分別為決定系數(shù)（R2）、均方根誤差（RMSE）與相對(duì)分析誤差（RPD）。當(dāng)RPD>2時(shí)，表明模型有極好的預(yù)測(cè)能力；當(dāng)RPD=1.4時(shí)，表明模型可對(duì)試驗(yàn)進(jìn)行粗略的估計(jì)；當(dāng)RPD<1.4時(shí)，表明模型無(wú)法對(duì)試樣進(jìn)行估計(jì)。

1.3.4 預(yù)測(cè)模型

支持向量機(jī)（Support Vector Machine，SVM）的基本模型是在特征空間中定義的劃分試樣的線性分類器，其學(xué)習(xí)策略是最大化劃分間隔。采用SVM建模的具體算法流程如下。

首先，在試樣空間中設(shè)有一個(gè)能將所有試樣準(zhǔn)確劃分開的超平面，任一試樣（x_i，y_i）到該超平面的函數(shù)距離γ ?與幾何距離γ如式（9）與式（10）所示

γ ?=y_i （ω^T X_i+d）。 （9）

γ=γ ?/‖ω‖ 。 （10）

式中：Xi為每個(gè)樣本都包含多個(gè)特征并由此構(gòu)成特征空間，X_i=[x_1，…，x_n ]∈χ；y_i用來(lái)區(qū)分該試樣為正相關(guān)還是負(fù)相關(guān)；ω為垂直于超平面的法向量；d為偏差，規(guī)定ωTx+d=0。

然后，根據(jù)上述得出的距離，可得出將訓(xùn)練數(shù)據(jù)集劃分為幾何距離最大的分類超平面。

最后，將優(yōu)化問(wèn)題拉格朗日化，并對(duì)拉格朗日化的問(wèn)題求最小值，可得到ω，將其帶入SVM的決策式可得式（11）

f（x）=ω^T x+d=∑_（i=1）^n?〖α_i y_i 〗 （X_i，x）+d。 （11）

式中：α_i為拉格朗日乘子，可以把約束條件融合到目標(biāo)函數(shù)里去，從而只用一個(gè)函數(shù)表達(dá)式便可將每個(gè)樣本準(zhǔn)確分類。

使用SVM算法建立預(yù)測(cè)模型時(shí)，需要先確定所使用的核函數(shù)以及相應(yīng)參數(shù)，常見的核函數(shù)包含多項(xiàng)式核函數(shù)、徑向基函數(shù)核、拉普拉斯核函數(shù)以及Sigmoid核函數(shù)。本試驗(yàn)采用SVM的徑向基核函數(shù)（RBF）。

在具體的試驗(yàn)過(guò)程中，利用K-S（Kolmogorov-Smirnov）劃分法將150份試樣數(shù)據(jù)按2∶1分類，其中，校正集樣本100個(gè)，預(yù)測(cè)集樣本50個(gè)，同時(shí)使用了隨機(jī)森林（RF），反向傳播神經(jīng)網(wǎng)絡(luò)（BP神經(jīng)網(wǎng)絡(luò)）2種算法作為對(duì)比試驗(yàn)組。

2 結(jié)果與分析

2.1 含節(jié)子木材力學(xué)性能的真實(shí)值

萬(wàn)能力學(xué)試驗(yàn)機(jī)最終以Excel圖表的格式輸出荷載與形變數(shù)據(jù)的關(guān)系曲線，該曲線的名稱為“試驗(yàn)力-位移曲線”，具體結(jié)果如圖2所示。

由圖2可以看出，當(dāng)試樣所受到的壓力荷載在300～4 400 N時(shí)，試樣所受荷載壓力與試樣表面位移可近似視為正比關(guān)系；當(dāng)試樣所受到的壓力荷載達(dá)到5 700 N時(shí)，試樣發(fā)生斷裂，力學(xué)破壞性試驗(yàn)結(jié)束，此時(shí)的荷載壓力經(jīng)過(guò)計(jì)算即可得到該試樣所對(duì)應(yīng)的MOR和MOE的真實(shí)值，其結(jié)果見表1。

在之前的研究中，Schimleck等［19］通過(guò)采集巴西紅木的橫截面光譜，分別對(duì)MOR和MOE進(jìn)行了預(yù)測(cè)，結(jié)果表明MOR比MOE在近紅外光譜的預(yù)測(cè)精度方面更加精確［4］，因此本研究選擇MOR作為模型訓(xùn)練所需預(yù)測(cè)的數(shù)據(jù)。

2.2 光譜數(shù)據(jù)處理結(jié)果及分析

圖3為原始光譜圖像。將采集到的150份試樣的原始光譜數(shù)據(jù)經(jīng)過(guò)去噪、去異常、求導(dǎo)，結(jié)果如圖4和圖5所示。

原始光譜圖像在經(jīng)過(guò)以上一系列處理后，由圖6可以更直觀地看到，這些處理方法有效地突出了各光譜數(shù)據(jù)的反射峰峰值，并去除了光譜圖像其他區(qū)域數(shù)值對(duì)光譜圖像特征點(diǎn)的影響。處理后的光譜圖像能夠使接下來(lái)的特征提取部分快速高效地進(jìn)行［19］，并保證力學(xué)性能預(yù)測(cè)模型的精確度。

通過(guò)SPA對(duì)進(jìn)行預(yù)處理后的光譜數(shù)據(jù)進(jìn)行特征波長(zhǎng)的選擇與提取，可得到以下結(jié)果。

根據(jù)變量數(shù)目的選擇與RMSE之間的關(guān)系可以得知，選擇18個(gè)變量時(shí)可以保證選擇最少的特征波長(zhǎng)，同時(shí)保留最多的光譜信息［20］，在本試驗(yàn)中這 18個(gè)特征波長(zhǎng)分別是938、1 049、1 235、1 131、1 282、1 299、1 330、1 357、1 386、1 429、1 468、1 506、 1 591、1 639、1 661、1 708、1 747、1 770 nm。

2.3 基于機(jī)器學(xué)習(xí)的含節(jié)子木材力學(xué)性能預(yù)測(cè)

將特征提取后的光譜數(shù)據(jù)按照2∶1的比例作為預(yù)測(cè)模型的訓(xùn)練集與測(cè)試集，并在訓(xùn)練之前，將10%誤差較大的預(yù)測(cè)集與訓(xùn)練集數(shù)據(jù)進(jìn)行了排除，以保證模型最大程度地不受測(cè)試數(shù)據(jù)誤差的影響。采用BP神經(jīng)網(wǎng)絡(luò)、SVM支持向量機(jī)與隨機(jī)森林RF模型進(jìn)行含節(jié)子木材MOR值的無(wú)損預(yù)測(cè)。其訓(xùn)練結(jié)果如圖7—圖9所示，橫坐標(biāo)為預(yù)測(cè)試樣編號(hào)，縱坐標(biāo)為預(yù)測(cè)結(jié)果MOR值。

SPA-BP、SPA-RF和SPA-SVM這3種模型通過(guò)相同的訓(xùn)練集進(jìn)行訓(xùn)練，使用相同的測(cè)試集進(jìn)行結(jié)果測(cè)試，其試驗(yàn)結(jié)果見表2。規(guī)定R2與RPD較大的同時(shí)RMSE較小的模型為最優(yōu)模型。

從測(cè)試集的預(yù)測(cè)結(jié)果與實(shí)際數(shù)據(jù)對(duì)比來(lái)看，SPA-BP模型訓(xùn)練評(píng)價(jià)指標(biāo)R2=0.85，RMSE=1.20，RPD=2.52，說(shuō)明預(yù)測(cè)結(jié)果與實(shí)際數(shù)據(jù)擬合程度一般，該模型能在對(duì)木材MOR值進(jìn)行預(yù)測(cè)方面性能較差；SPA-RF模型訓(xùn)練評(píng)價(jià)指標(biāo)R2=0.75，RMSE=1.50，RPD=1.99，預(yù)測(cè)結(jié)果與實(shí)際數(shù)據(jù)擬合程度較差，說(shuō)明該模型無(wú)法應(yīng)用于對(duì)木材MOR值進(jìn)行預(yù)測(cè)；SPA-SVM訓(xùn)練評(píng)價(jià)指標(biāo)R2=0.95，RMSE=0.59，RPD=4.99，預(yù)測(cè)結(jié)果與實(shí)際數(shù)據(jù)的匹配程度非常好，表明該模型能夠很好地預(yù)測(cè)含節(jié)子木材的MOR值。

3 結(jié)論

本研究采用目標(biāo)檢測(cè)算法對(duì)多光譜圖像中的數(shù)字圖像部分進(jìn)行節(jié)子區(qū)域識(shí)別，輸出節(jié)子區(qū)域面積與位置后，通過(guò)計(jì)算得到節(jié)子區(qū)域和無(wú)瑕疵區(qū)域同整個(gè)木材試樣橫截面的慣性矩比，之后在節(jié)子區(qū)域與無(wú)瑕疵區(qū)域各隨機(jī)采集20個(gè)點(diǎn)的光譜數(shù)據(jù)，以上面求出的慣性矩比例作為加權(quán)系數(shù)，將每個(gè)試樣選擇的光譜數(shù)據(jù)求出加權(quán)平均值，作為代表該試樣的光譜值，再用SPA對(duì)預(yù)處理后的光譜數(shù)據(jù)進(jìn)行特征提取，將其輸入到SVM預(yù)測(cè)模型中，實(shí)現(xiàn)對(duì)含節(jié)子木材抗彎強(qiáng)度的無(wú)損預(yù)測(cè)。從預(yù)測(cè)結(jié)果與真實(shí)值的比較來(lái)看，此方法可以達(dá)到不通過(guò)損壞木板來(lái)評(píng)判木材力學(xué)性能的目的。

【參 考 文 獻(xiàn)】

［1］ 張盼盼.我國(guó)森林資源發(fā)展?fàn)顩r與趨勢(shì)分析［J］.綠色科技，2019（17）：204-206.

ZHANG P P.Analysis of the development status and trends of forest resources in china［J］.Green Technology，2019（17）：204-206.

［2］ 初偉杰.探析木材檢驗(yàn)工作的重要性［J］.農(nóng)村實(shí)用科技信息，2015（7）：32.

CHU W J.Exploring the importance of wood inspection work［J］.Practical Scientific and Technological Information in Rural Areas，2015（7）：32.

［3］ HORVATH L，PESZLEN I，PERALTA P，et al.Use of transmittance near-infrared spectroscopy to predict the mechanical properties of 1- and 2-year-old transgenic aspen［J］.Wood Science and Technology，2011，45（2）：303-314.

［4］ SCHIMLECK L R，DE MATOS，JIM OILVEIRA，et al.Non-destructive estimation of pernambuco（Caesalpinia echinata） clear wood properties using near infrared spectroscopy［J］.Journal of Near Infrared Spectorscopy，2011，19（5）：411-419.

［5］ KOTHIYAL V，RATURI A.Estimating mechanical properties and specific gravity for five-year-old Eucalyptus tereticornis having broad moisture content range by NIR spectroscopy［J］.Holzforschung，2011，65（5）：757-762.

［6］ XU Q H，XU M H，NI Y H，et al.Progresses and perspectives of near infrared spectroscopy analytical technology［J］.Canadian Journal of Forest Research，2011，41（2）：352-358.

［7］ LIANG H，CAO J，TU W J，et al.Nondestructive determination of the compressive strength of wood using near-infrared spectroscopy［J］.Bioresources，2016，11（3）：7205-7213．

［8］ 褚小立，陳瀑，李敬巖，等.近紅外光譜分析技術(shù)的最新進(jìn)展與展望［J］.分析測(cè)試學(xué)報(bào)，2020，39（10）：1181-1188.

CHU X L，CHEN P，LI J Y，et al.The latest progress and prospects of near-infrared spectroscopy analysis technology［J］.Journal of Analysis and Testing，2020，39（10）：1181-1188.

［9］ FUJIMOTO T，KOBORI H，TSUCHIKAWA S.Prediction of wood density independently of moisture conditions using near infrared spectroscopy［J］.Journal of Near Infrared Spectorscopy，2012，20（3）：353-359.

［10］ 謝富，朱定局.深度學(xué)習(xí)目標(biāo)檢測(cè)方法綜述［J］.計(jì)算機(jī)系統(tǒng)應(yīng)用，2022，31（2）：1-12.DOI：10.15888/j.cnki.csa.008303.

XIE F，ZHU D J.Overview of deep learning object detection methods［J］.Application of Computer Systems，2022，31（2）：1-12.DOI：10.15888/j.cnki.csa.008303.

［11］ LISBETH G T.Determination of dry matter content and basic density of norway spruce by near infrared reflectance and transmittance spectroscopy［J］.Journal of Near Infrared Spectroscopy，1994，2（3）.

［12］ 李明月，陳東升，王雪玉，等.節(jié)子的性質(zhì)及其對(duì)木材材性的影響［J］.西北林學(xué)院學(xué)報(bào)，2020，35（4）：197-204.

LI M Y，CHEN D S，WANG X Y，et al.Nature of the knot and its impact on wood properties［J］.Journal of Northwest Forestry University，2020，35（4）：197-204.

［13］ YUN S，HAN D，OH S J，et al.CutMix： regularization strategy to train strong classifiers with localizable features［C］.2019 IEEE/CVF International 36FFH9EizlS/bT/0U/M44w==Conference on Computer Vision（ICCV），Seoul，Korea，2019：33-42.

［14］ 張靈帥，王衛(wèi)東，谷運(yùn)紅，等.近紅外光譜的主成分分析-馬氏距離聚類判別用于卷煙的真?zhèn)舞b別［J］.光譜學(xué)與光譜分析，2011，31（5）：1254-1257.

ZHANG S L，WANG W D，GU Y H，et al.Principal component analysis of near infrared spectroscopy-Markov distance clustering discrimination for authenticity identification of cigarettes［J］.Spectroscopy and Spectral Analysis，2011，31（5）：1254-1257.

［15］ CHEN J H，YU H L，JIANG D P，et al.A novel NIRS modelling method with OPLS-SPA and MIX-PLS for timber evaluation［J］.Journal of Forestry Research，2022，33（1）：369-376.

［16］ YANG Y，WANG F H，JIANG D P.PSO-SVR model for ecology property of indoor decoration surface materials［J］.Journal of Northeast Forestry University，2019，47（4）：81-85．

［17］ SCHIMLECK L R，EVANS R，MATHESON A C，et al.Estimation of Pinus radiata D.Don clear wood properties by near-infrared spectroscopy［J］.Journal of Wood Science，2002，48（2）：132-137

［18］ 肖雄亮，陳長(zhǎng)明.乙醇濃度預(yù)測(cè)的多元線性回歸模型建立及驗(yàn)證［J］.紅外技術(shù)，2021，43（12）：1228-1233.

XIAO X L，CHEN C M.Establishment and validation of a multiple linear regression model for predicting ethanol concentration［J］.Infrared Technology，2021，43（12）：1228-1233.

［19］ SCHIMLECK L R，EVANS R，MATHESON A C.Estimation of Pinus radiata D.Don clear wood properties by near-infrared spectroscopy［J］.Journal of Wood Science，2002，48（2）：132-137.

［20］ WU S Y.Information system establishment and prediction analysis for SVM model of RBF kernel function［J］.Journal of Physics：Conference Series，2021，1982（1）：81-93.