大型散貨堆體積的快速測(cè)量

楊德山，董麗麗，梁倩倩，許文海

(大連海事大學(xué) 信息科學(xué)技術(shù)學(xué)院，遼寧 大連 116026)

?

大型散貨堆體積的快速測(cè)量

楊德山，董麗麗*，梁倩倩，許文海

(大連海事大學(xué) 信息科學(xué)技術(shù)學(xué)院，遼寧 大連 116026)

針對(duì)現(xiàn)有散貨測(cè)量系統(tǒng)對(duì)堆場(chǎng)環(huán)境適應(yīng)性差、盤點(diǎn)時(shí)間長(zhǎng)、效率低、操作復(fù)雜等不足，提出了一種散貨堆體積快速測(cè)量方法。同時(shí)，利用二維激光掃描儀、差分GPS和姿態(tài)測(cè)量系統(tǒng)設(shè)計(jì)了一種體積測(cè)量系統(tǒng)。該系統(tǒng)用激光掃描儀動(dòng)態(tài)測(cè)量堆體表面的幾何信息，用姿態(tài)測(cè)量系統(tǒng)實(shí)時(shí)測(cè)量掃描儀的空間姿態(tài)數(shù)據(jù)，用GPS測(cè)出掃描儀在測(cè)量過(guò)程中的三維位置；最后通過(guò)數(shù)據(jù)融合計(jì)算形成堆體的三維點(diǎn)云，利用點(diǎn)云獲得散貨堆體積。文中基于單條堆體輪廓點(diǎn)云特征，提出快速堆體下邊緣查找算法來(lái)去除掃描過(guò)程中地面點(diǎn)云的誤差影響；采用投影剖分法完成完整堆體點(diǎn)云計(jì)算體積。實(shí)驗(yàn)顯示，利用本文設(shè)計(jì)的測(cè)試系統(tǒng)可在30 s內(nèi)完成體積為69 m3的標(biāo)準(zhǔn)堆體測(cè)量，平均相對(duì)誤差為0.42%，重復(fù)測(cè)量誤差為0.41%。在實(shí)際散貨堆實(shí)驗(yàn)中，可在10 min內(nèi)完成大小約為31 500 m3的散貨堆測(cè)量，4個(gè)不同料堆體積測(cè)量的平均重復(fù)測(cè)量誤差為0.74%。結(jié)果表明，本方法可在保證測(cè)量精度的同時(shí)，簡(jiǎn)單、高效地測(cè)量散貨堆體積。

激光測(cè)量；體積測(cè)量；快速測(cè)量；大型散貨堆；點(diǎn)云邊緣；三角剖分

(CollegeofInformationSciencesandTechnology，DalianMaritimeUniversity，Dalian116026，China)*Correspondingauthor，E-mail：dll_lili@163.com

1　引　言

在電廠、港口、儲(chǔ)煤基地和糧倉(cāng)等大型散貨集散基地會(huì)堆存大量煤炭、礦石、糧食等散貨，快速、準(zhǔn)確獲取當(dāng)前堆場(chǎng)的儲(chǔ)量是國(guó)家和企業(yè)進(jìn)行經(jīng)濟(jì)規(guī)劃、戰(zhàn)略儲(chǔ)備計(jì)劃、生產(chǎn)安排的關(guān)鍵。

隨著傳感器和激光技術(shù)的發(fā)展，光學(xué)三維建模系統(tǒng)的應(yīng)用于越來(lái)越多[1-6]，基于三維模型可以進(jìn)行堆場(chǎng)儲(chǔ)量的測(cè)量[7-10]。熊友輝等[11]提出一種便攜式激光測(cè)量?jī)x，通過(guò)多站人工打點(diǎn)的方式進(jìn)行盤點(diǎn)，雖然該儀器存在測(cè)量周期長(zhǎng)、精度差、人為影響大等缺點(diǎn)，但仍是現(xiàn)在堆場(chǎng)測(cè)量的常用方法。張德津等[12]提出了一種基于多傳感器的固定式激光測(cè)量方法，該方法將激光掃描儀安裝在外界載體(如取料機(jī)等)上對(duì)堆場(chǎng)進(jìn)行測(cè)量。盛業(yè)華等[13]提出了一種基于三維激光掃描點(diǎn)云的多站拼接方法，該方法圍繞被測(cè)目標(biāo)布設(shè)多個(gè)測(cè)量點(diǎn)，將各測(cè)量獲得的點(diǎn)云進(jìn)行拼接進(jìn)而獲得完整的被測(cè)堆體模型。張小虎等[14]提出了激光投影線輔助下的堆場(chǎng)三維形貌攝影測(cè)量方法，將可見(jiàn)光圖像和激光掃描點(diǎn)云融合建立被測(cè)堆場(chǎng)的三維模型。王海波等[15]提出了車載式多站盤點(diǎn)系統(tǒng)?？傮w來(lái)說(shuō)，目前的堆場(chǎng)測(cè)量系統(tǒng)仍然依賴外界設(shè)備，存在測(cè)量周期長(zhǎng)、實(shí)時(shí)性差、過(guò)程復(fù)雜、人為影響大等缺點(diǎn)。

本文提出了一種散貨堆體積快速測(cè)量方法，利用二維激光掃描儀、差分全球定位系統(tǒng)(Global Position System，GPS)和姿態(tài)測(cè)量系統(tǒng)組合設(shè)計(jì)了一種手提式體積快速測(cè)量系統(tǒng)。相對(duì)基點(diǎn)測(cè)量盤煤儀如便攜式盤煤儀、車載定點(diǎn)式盤煤儀，本文方法既不需要圍繞煤堆選擇測(cè)量基點(diǎn)，也不需要在測(cè)量基點(diǎn)上對(duì)堆體打點(diǎn)測(cè)量；相對(duì)于固定式盤煤儀和投影輔助盤煤儀，本文方法不依賴于取料機(jī)等外界設(shè)備，數(shù)據(jù)處理放于盤煤儀內(nèi)部，測(cè)量完成后可直接生成點(diǎn)云并計(jì)算得到結(jié)果，不需額外的數(shù)據(jù)傳輸。測(cè)量人員僅需手提該測(cè)量系統(tǒng)圍繞被測(cè)堆體行走一周即可實(shí)時(shí)得到被測(cè)堆體的體積。實(shí)驗(yàn)證明，本系統(tǒng)在保證測(cè)量精度的同時(shí)，具有適應(yīng)能力強(qiáng)、快速高效、操作簡(jiǎn)便和成本低等特點(diǎn)，有良好的應(yīng)用價(jià)值。

2　測(cè)量原理與算法

本文提出的散貨堆測(cè)量方法的測(cè)量示意圖如圖1所示。其測(cè)量原理為：當(dāng)測(cè)量人員手提測(cè)量系統(tǒng)圍繞堆體行走時(shí)，激光掃描儀動(dòng)態(tài)測(cè)量堆體表面的幾何信息；同時(shí)姿態(tài)測(cè)量系統(tǒng)實(shí)時(shí)測(cè)量掃描儀的空間姿態(tài)數(shù)據(jù)，GPS測(cè)出掃描儀在測(cè)量過(guò)程中的三維位置；經(jīng)過(guò)數(shù)據(jù)融合計(jì)算形成堆體的三維點(diǎn)云，并利用點(diǎn)云進(jìn)行體積計(jì)算。

以開(kāi)始測(cè)量時(shí)刻GPS天線的相位中心為原點(diǎn)建立全局測(cè)量坐標(biāo)系，以水平正北為Y軸，Z軸沿重力鉛垂線向上，X軸垂直于YOZ平面指向正東。測(cè)量系統(tǒng)的坐標(biāo)系主要有激光掃描的局部坐標(biāo)系(簡(jiǎn)稱L坐標(biāo)系)，姿態(tài)局部坐標(biāo)系(簡(jiǎn)稱I坐標(biāo)系)和全局測(cè)量坐標(biāo)系(簡(jiǎn)稱M坐標(biāo)系)。

圖1　手提測(cè)量系統(tǒng)的測(cè)量示意圖

針對(duì)二維激光掃描儀，設(shè)任一腳點(diǎn)的測(cè)距為S，掃描角為θ，則該腳點(diǎn)在全局測(cè)量坐標(biāo)系中的坐標(biāo)為：

(1)

先沿Z軸順時(shí)針旋轉(zhuǎn)偏航角α，對(duì)應(yīng)的旋轉(zhuǎn)矩陣為：

(2)

再分別沿X軸和Y軸順時(shí)針旋轉(zhuǎn)橫滾角β和俯仰角γ，對(duì)應(yīng)的旋轉(zhuǎn)矩陣分別為：

(3)

(4)

由于坐標(biāo)系符合右手定則，因此繞y軸的旋轉(zhuǎn)變換矩陣與繞X軸和Z軸的變換矩陣在符號(hào)上有所不同。于是可以得到：

RGI=Ry(γ)Rx(β)Rz(α).

(5)

同理，RIL為分別沿Z軸，X軸和Y軸順時(shí)針旋轉(zhuǎn)安裝誤差角Δα，Δβ，Δγ的旋轉(zhuǎn)矩陣，即：

RIL=Ry(Δγ)Rx(Δβ)Rz(Δα).

(6)

將式(6)和式(5)帶入式(1)可將激光掃描儀獲取的堆體表面信息轉(zhuǎn)換為空間三維點(diǎn)云數(shù)據(jù)。

2.2地面點(diǎn)去除處理

在測(cè)量過(guò)程中，激光掃描儀不可避免掃描到被測(cè)堆體附近的地面，給堆體三維建模和體積計(jì)算帶來(lái)額外誤差，因此需要將地面點(diǎn)從整體點(diǎn)云中去除。由于在大型散貨集散基地的堆體是由堆料機(jī)堆砌成垛，所以會(huì)與地面形成明顯的分界，即激光掃描點(diǎn)云中俗稱的堆體下邊緣。由于二維激光掃描點(diǎn)為順序掃描，因此找到堆體下邊緣點(diǎn)即可去除地面點(diǎn)。目前，已有學(xué)者研究了許多點(diǎn)云邊緣提取算法[16]，本文針對(duì)掃描堆體的輪廓特征，提出一種快速邊緣查找算法，算法流程如下：

(1)在局部坐標(biāo)系y-z平面內(nèi)，一條激光掃描線得到的堆體輪廓點(diǎn)云如圖2所示。計(jì)算輪廓點(diǎn)云生成后的第一個(gè)點(diǎn)到最后一個(gè)點(diǎn)的連線的斜率，并求反正切，得到斜率對(duì)應(yīng)的角度為：

(7)

圖2　單條輪廓點(diǎn)云

(2)將該條掃描線對(duì)應(yīng)的輪廓點(diǎn)云在y-z平面內(nèi)順時(shí)針旋轉(zhuǎn)角度θ，即有：

(8)

(3)比較旋轉(zhuǎn)后點(diǎn)云的z坐標(biāo)，選擇其中z坐標(biāo)最小的點(diǎn)對(duì)應(yīng)的標(biāo)號(hào)k，此點(diǎn)即為地堆分界點(diǎn)。如圖3所示。

Zk=min{Zi}，i=1，2，…，n.

(9)

圖3　分界點(diǎn)提取

(4)在原始掃描點(diǎn)云中去除標(biāo)號(hào)大于k的地面點(diǎn)云。

由于雨水沖刷等原因，堆體也會(huì)存在下邊緣不明顯的特殊情況，如圖4所示。當(dāng)從掃描輪廓點(diǎn)云中無(wú)法直接找到顯著分界點(diǎn)時(shí)，本文首先根據(jù)輪廓點(diǎn)云分別擬合出堆體點(diǎn)和地面點(diǎn)直線；然后尋找兩條擬合線的交點(diǎn)m，將該點(diǎn)定義為擬合分界點(diǎn)；再計(jì)算輪廓點(diǎn)云中在兩條擬合線相交范圍內(nèi)各點(diǎn)(圖4中p1～p3點(diǎn))到擬合分界點(diǎn)m的距離；最后將距離m最近的點(diǎn)(圖4中為p1點(diǎn))取為分界點(diǎn)。通過(guò)直線擬合確定分界點(diǎn)依賴于擬合直線交點(diǎn)的準(zhǔn)確計(jì)算，而由于實(shí)際堆體表面和地面的平整度直接影響了堆體擬合直線的特性，最終體積的計(jì)算誤差會(huì)增大。因此通常情況下當(dāng)堆體分界點(diǎn)明顯時(shí)，本文采用直接旋轉(zhuǎn)求最小值的方法求取下邊緣分界點(diǎn)。

圖4　不明顯分界點(diǎn)提取方法

2.3投影Delaunay三角剖分體積算法

因?yàn)樯⒇浂训亩洋w形狀在空間上沒(méi)有層疊區(qū)域，因此可以利用投影法來(lái)計(jì)算體積。首先將上述得到的點(diǎn)云投影到xy平面上，對(duì)投影后的點(diǎn)云進(jìn)行Delaunay三角剖分，剖分為互不包含的獨(dú)立三角形，且每個(gè)三角形對(duì)應(yīng)堆體表面的一個(gè)三角面片，三角面片和投影三角形構(gòu)成一個(gè)組合體，如圖5所示。計(jì)算剖分得到的所有三角形對(duì)應(yīng)的組合體的體積并求和，即為整個(gè)散貨堆的體積V：

Vi=VABC-A′B′C′+VB-B″A″C+VA-BA″C，

(10)

(11)

由此可知，三棱柱的數(shù)量越多，體積測(cè)量越準(zhǔn)確，而激光掃描儀采集的點(diǎn)云數(shù)目決定了三棱柱的數(shù)量。

圖5　三角面片組合體體積的計(jì)算模型

2.4誤差來(lái)源及評(píng)價(jià)指標(biāo)

由點(diǎn)云生成公式可知，影響點(diǎn)云生成準(zhǔn)確性的誤差因素主要有兩類：一類是測(cè)量設(shè)備的測(cè)量誤差，主要有激光測(cè)距S和角度θ的誤差，姿態(tài)測(cè)量誤差，GPS定位誤差以及航向角測(cè)量誤差等；另一類是系統(tǒng)安裝誤差，主要包括不同設(shè)備集成安裝過(guò)程中的偏心量誤差和安置角誤差。其中，系統(tǒng)安裝誤差為靜態(tài)常量，通過(guò)系統(tǒng)標(biāo)定可消除此部分誤差的影響；而測(cè)量誤差為隨機(jī)誤差，不同的測(cè)量誤差同時(shí)存在且交叉影響，是誤差分析的難點(diǎn)。因此，需要對(duì)系統(tǒng)進(jìn)行標(biāo)定研究[17]。

激光點(diǎn)云生成的準(zhǔn)確性直接影響體積測(cè)量結(jié)果的準(zhǔn)確性。本文首先通過(guò)搭建標(biāo)準(zhǔn)理想堆體模型來(lái)驗(yàn)證系統(tǒng)的相對(duì)誤差E，相對(duì)誤差可以反映系統(tǒng)測(cè)量的準(zhǔn)確性。

(12)

式中：xm為測(cè)量值，xr為真值。

由于實(shí)際散貨堆的體積真值難以獲得，因此在驗(yàn)證系統(tǒng)的準(zhǔn)確性后，以重復(fù)測(cè)量誤差來(lái)衡量系統(tǒng)的測(cè)量可靠性。測(cè)量?jī)x器的重復(fù)性是測(cè)量?jī)x器隨機(jī)誤差分量的表現(xiàn)，通過(guò)測(cè)量?jī)x器進(jìn)行連續(xù)多次測(cè)量并和參考值進(jìn)行比較，計(jì)算實(shí)驗(yàn)標(biāo)準(zhǔn)差來(lái)表示被評(píng)定測(cè)量?jī)x器的重復(fù)性。一般用貝塞爾公式計(jì)算實(shí)驗(yàn)標(biāo)準(zhǔn)差S。以貝塞爾公式得到的標(biāo)準(zhǔn)差除以測(cè)量平均值代表系統(tǒng)的重復(fù)性誤差er。

(13)

(14)

3　標(biāo)準(zhǔn)模型實(shí)驗(yàn)

3.1測(cè)量系統(tǒng)設(shè)計(jì)

利用上述系統(tǒng)原理與算法，本文設(shè)計(jì)研發(fā)了手提式體積測(cè)量系統(tǒng)，系統(tǒng)采用SICK LMS511激光掃描儀和天寶差分GPS，利用雙天線GPS和動(dòng)態(tài)傾角儀測(cè)量姿態(tài)。雙天線GPS由兩個(gè)伸縮桿連接，避免了天線受到人體遮擋的影響。測(cè)量系統(tǒng)的實(shí)物圖如圖6所示，該系統(tǒng)長(zhǎng)45 cm，寬20 cm，高20 cm，質(zhì)量為12.5 kg，對(duì)一般被測(cè)散貨物品的光學(xué)特性無(wú)特殊要求。

圖6　手提測(cè)量系統(tǒng)實(shí)物照片

本文在系統(tǒng)研制完成后，首先對(duì)系統(tǒng)整體進(jìn)行了安裝誤差標(biāo)定，消除了安裝誤差對(duì)測(cè)量結(jié)果的影響。

3.2標(biāo)準(zhǔn)散貨模型測(cè)量

根據(jù)散貨堆的形狀特點(diǎn)，這里搭建了標(biāo)準(zhǔn)四棱臺(tái)散貨堆模型，模型底面為8 m×7 m的長(zhǎng)方形，堆高3 m，上表面為2 m×1 m的長(zhǎng)方形，標(biāo)準(zhǔn)模型的真實(shí)體積為69 m3，忽略模型在室外由于人工堆放及熱脹冷縮等原因造成的誤差。

利用本文研制的測(cè)量系統(tǒng)對(duì)模擬散貨堆進(jìn)行10次測(cè)量，生成的點(diǎn)云模型如圖8所示，10次體積測(cè)量結(jié)果如圖9所示。

圖7　標(biāo)準(zhǔn)散貨模型實(shí)物圖

圖8　標(biāo)準(zhǔn)散貨模型點(diǎn)云圖

圖9　標(biāo)準(zhǔn)堆10次測(cè)量結(jié)果

圖10　10次測(cè)量的相對(duì)誤差

根據(jù)圖9的測(cè)量結(jié)果，由式(12)和式(14)計(jì)算得到10次體積測(cè)量的平均相對(duì)誤差為0.42%，測(cè)量的重復(fù)性誤差為0.41%。結(jié)果表明，本文設(shè)計(jì)的手提式體積快速測(cè)量系統(tǒng)具有良好的測(cè)量精度。實(shí)驗(yàn)表明，該測(cè)量系統(tǒng)可在30 s內(nèi)得到被測(cè)堆體的空間三維點(diǎn)云圖和體積，在保證準(zhǔn)確度的情況下，有效提高了測(cè)量的速度與效率。

3.3散貨實(shí)物實(shí)驗(yàn)

利用自主研發(fā)的手提式體積快速測(cè)量系統(tǒng)，在秦皇島、天津和黃驊港神華儲(chǔ)煤基地進(jìn)行實(shí)地測(cè)量實(shí)驗(yàn)，選取4個(gè)體積不同的大型煤堆進(jìn)行典型測(cè)量，分別對(duì)每個(gè)煤堆進(jìn)行10次測(cè)量，以重復(fù)測(cè)量誤差來(lái)衡量測(cè)量結(jié)果的可靠性。測(cè)量結(jié)果如表1所示，所提取的堆體輪廓點(diǎn)云分界點(diǎn)如圖11所示。測(cè)量煤堆1得到的點(diǎn)云圖如圖12所示。

圖11　實(shí)際堆體輪廓點(diǎn)云及分界點(diǎn)圖

圖12　煤堆1的點(diǎn)云圖

由表1中4個(gè)煤堆的10次體積測(cè)量結(jié)果可以看出，系統(tǒng)的測(cè)量重復(fù)性誤差不會(huì)受到堆體體積的影響。計(jì)算結(jié)果顯示，系統(tǒng)的平均重復(fù)測(cè)量誤差為0.74%。其中，體積約為31 500 m3的煤堆1測(cè)量一次所需的時(shí)間小于10 min，而利用目前堆場(chǎng)使用的打點(diǎn)測(cè)量盤煤儀至少需要1 h。實(shí)踐表明，本文設(shè)計(jì)的手提式體積快速測(cè)量系統(tǒng)具有良好的重復(fù)準(zhǔn)確性，可在保證測(cè)量精度的條件下，快速獲得大型煤堆的體積。

表1　4個(gè)煤堆的體積測(cè)量結(jié)果

4　結(jié)　論

本文提出了一種大型散貨堆體積快速測(cè)量方法，解決了大型散貨堆場(chǎng)的實(shí)時(shí)快速體積測(cè)量問(wèn)題；提出了基于堆體特征輪廓的點(diǎn)云邊緣查找算法，從而得到完整的堆體點(diǎn)云模型，并利用投影Delaunay三角剖分和棱柱體積累加算法計(jì)算堆體體積。標(biāo)準(zhǔn)模型堆測(cè)量實(shí)驗(yàn)結(jié)果表明，利用本文方法設(shè)計(jì)的測(cè)量系統(tǒng)測(cè)量準(zhǔn)確，平均相對(duì)誤差為0.42%，重復(fù)測(cè)量誤差為0.41%。在黃驊港、秦皇島和天津港等神華儲(chǔ)煤基地選取多個(gè)煤堆進(jìn)行實(shí)地測(cè)量。結(jié)果表明，該方法可高效快速地測(cè)量大型煤堆的體積，在10 min內(nèi)完成31 500 m3堆體的測(cè)量，體積測(cè)量的平均重復(fù)測(cè)量誤差為0.74%。

本文設(shè)計(jì)的手提式體積快速測(cè)量系統(tǒng)已協(xié)助神華集團(tuán)完成黃驊港和天津港儲(chǔ)煤基地的年度審計(jì)盤點(diǎn)工作，實(shí)踐表明本系統(tǒng)在保證測(cè)量準(zhǔn)確度的情況下，大大提高了測(cè)量效率，節(jié)省了盤點(diǎn)的成本和時(shí)間。該系統(tǒng)可以在電廠、大型儲(chǔ)煤基地等企業(yè)進(jìn)行推廣應(yīng)用，具有良好的實(shí)用前景。

[1]邊心田，蘇顯渝，陳文強(qiáng).基于反向攝影的三維坐標(biāo)測(cè)量方法[J].中國(guó)激光，2010，37(7)：1832-1836.

BIAN X T，SU X Y， CHEN W Q.A novel tree-dimensional coordinates measurement method based on inverse photogrammetry[J].ChineseJ.Laser，2010，37(7):1832-1836.(in Chinese)

[2]戴彬，鐘若飛，胡競(jìng). 基于車載激光掃描數(shù)據(jù)的城市地物三維重建研究[J].首都師范大學(xué)學(xué)報(bào)：自然科學(xué)版，2011，32(3):89-96.

DAI B，ZHONG R F，HU J.Research on 3D reconstruction of urban features from data based on vehicle-borne laser scanning[J].JournalofCapitalNormalUniversity：NaturalScienceEdition，2011，32(3):89-96.(in Chinese)

[3]林義閩，呂乃光，婁小平，等.用于弱紋理場(chǎng)景三維重建的機(jī)器人視覺(jué)系統(tǒng)[J].光學(xué) 精密工程，2015，23(2):540-549.

LIN Y M，Lü N G，LOU X P，etal.. Robot vision system for 3D reconstruction in low texture environment[J].Opt.PrecisionEng.，2015，23(2):540-549.(in Chinese)

[4]AHN J K， LEE K Y， SIM J Y，etal..Large-scale 3D point cloud compression using adaptive radial distance prediction in hybrid coordinate domains[J].IEEEJournalofSelectedTopicsinSignalProcessing，2015，9(3)：422-434.

[5]CHEN J，WU X J，WANG M Y，etal..3D shape modeling using a self-developed hand-held 3D laser scanner and an efficient HT-ICP point cloud registration algorithm [J].Optics&LaserTechnology，2013，45： 414-423.

[6]時(shí)光，張福民，曲興華.經(jīng)濟(jì)型大尺寸激光自動(dòng)坐標(biāo)測(cè)量系統(tǒng)[J].光學(xué) 精密工程，2013，21(8):1957-1965.

SHI G，ZHANG F M，QU X H.Economical type large-scale laser automatic coordinate measuring system[J].Opt.PrecisionEng.，2013，21(8) :1957-1965.(in Chinese)

[7]周森，郭永彩，高潮，等.基于三維激光掃描的移動(dòng)大尺寸圓柱體工件長(zhǎng)度快速檢測(cè)系統(tǒng)[J].光學(xué) 精密工程，2014，22(6):1524-1530.ZHOU S，GUO Y C，GAO CH，etal..Rapid length measuring system for mobile and large scale cylinder workpieces based on 3D laser scanning[J].Opt.PrecisionEng.， 2014，22(6):1524-1530.(in Chinese)

[8]邵晴，徐濤，吉野辰萌，等.基于三維激光掃描技術(shù)的糧食儲(chǔ)量監(jiān)測(cè)系統(tǒng)的設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2015，31(20):262-267.

SHAO Q，XU T， YOSHINO T，etal..Design and experiment for grain storage monitoring system based on 3-D laser scanning technology[J].TransactionsoftheChineseSocietyofAgriculturalEngineering，2015，31(20): 262-267.(in Chinese)

[9]HYYPPA J，KELLE O，LEHIKOINEN M，etal.. A segmentation-based method to retrieve stem volume estimates from 3-D tree height models produced by laser scanners[J].IEEETransactionsonGeoscienceandRemoteSensing，2011，39(5)： 969-975.

[10]陳展鵬，雷延武，晏清洪，等.汶川震區(qū)滑坡堆積體體積三維激光掃描儀測(cè)量與計(jì)算方法[J].農(nóng)業(yè)工程學(xué)報(bào)，2013，29(8)：135-144.CHEN ZH P，LEI Y W，YAN Q H，etal..Measuring and calculation methods for landslide volume with 3-D laser scanner in Wenchuan earthquake area[J].TransactionsoftheChineseSocietyofAgriculturalEngineering， 2013，29(8)：135-144.(in Chinese)

[11]熊友輝，李培生，鄒顯宏，等.便攜式激光盤煤系統(tǒng)原理及應(yīng)用[J].中國(guó)電力，2003，36(6):48-51.XIONG Y H，LI P SH，ZOU X H，etal..Study for the method of automatic disk coal in the open coal yard [J].ElectricPower，2003，36(6):48-51.(in Chinese)

[12]張德津，李必軍，何莉.基于多傳感器集成的堆場(chǎng)激光測(cè)量技術(shù)應(yīng)用[J].中國(guó)激光，2012(2):159-164.

ZHANG D J，LI B J，HE L.Application on laser measurement for large storage yard based on multi-sensor integration[J].ChineseJ.Lasers，2012(2):159-164.(in Chinese)

[13]盛業(yè)華，張卡，張凱，等.地面三維激光掃描點(diǎn)云的多站數(shù)據(jù)無(wú)縫拼接[J].中國(guó)礦業(yè)大學(xué)學(xué)報(bào)，2010，39(2):233-237.

SHENG Y H，ZHANG K，ZHANG K，etal..Seamless multi-station merging of terrestrial laser scanned 3D point clouds[J].JournalofChinaUniversityofMining&Technology， 2010，39(2):233-237.(in Chinese)

[14]張小虎，歐建良，苑云，等.投影輪廓線輔助下的堆場(chǎng)三維形貌攝影測(cè)量研究[J].光學(xué)學(xué)報(bào)，2011(6):99-106.

ZHANG X H，OU J L，YUAN Y，etal..Stack-yard 3D measurement based on photogrammetry and projected contour aiding[J].ActaOpticaSinica，2011(6):99-106.(in Chinese)

[15]王海波，張德津，何莉.大型露天料場(chǎng)激光測(cè)量方法研究[J].中國(guó)激光，2013(5):189-196.WANG H B，ZHANG D J，HE L.Research on large open stockyard laser measurement methods[J].ChineseJ.Lasers， 2013(5):189-196.(in Chinese)

[16]譚賁，鐘若飛，李芹.車載激光掃描數(shù)據(jù)的地物分類方法[J].遙感學(xué)報(bào)，2012(1):50-66.

TAN B，ZHONG R F，LI Q.Objects classification with vehicle-borne laser scanning data[J].JournalofRemoteSensing，2012(1):50-66.(in Chinese)

[17]章大勇，吳文啟，吳美平.機(jī)載激光雷達(dá)系統(tǒng)標(biāo)定方法[J].光學(xué) 精密工程，2009，17(11):2806-2813.

ZHANG D Y，WU W Q，WU M P.Calibration technology of airborne lidar[J].Opt.PrecisionEng.，2009，17(11):2806-2813.(in Chinese)

董麗麗(1980-)，女，博士，副教授，2002年、2004年、2008年于哈爾濱工業(yè)大學(xué)分別獲得學(xué)士、碩士、博士學(xué)位，主要研究方向?yàn)楣怆娦畔⑴c光電檢測(cè)等。E-mail：dll_lili@163.com

導(dǎo)師簡(jiǎn)介：

許文海(1956-)，男，教授，博士生導(dǎo)師，1982年、1984年、1991年于哈爾濱工業(yè)大學(xué)分別獲得學(xué)士、碩士、博士學(xué)位，1993年于東京工業(yè)大學(xué)獲得博士學(xué)位，主要從事水下探測(cè)、光電檢測(cè)與光電信息等方面的研究。E-mail:whxu@dlmu.edu.cn

(版權(quán)所有未經(jīng)許可不得轉(zhuǎn)載)

Rapid volume measurement for large bulk cargo

YANG De-shan， DONG Li-li*， LIANG Qian-qian， XU Wen-hai

A rapid volume measurement method for a bulk cargo was proposed to overcome the shortcomings of traditional bulk measuring systems in poor adaptability for yard environments， longer complete time， lower efficiency and complex operation. Then a volume measuring system consisting of a 2D laser scanner， a difference GPS and an attitude measuring system was designed. For this system， the laser scanner was used to measure dynamically the surface geometry information of a bulk， the attitude measuring system was taken to measure the space attitude data of the scanner in real time and the GPS was utilized to obtain the 3D position of the scanner in measuring processing. Finally， the data fusion was used to calculate 3D point cloud of the bulk cargo and to obtain the volume of the bulk cargo. On the basis of the feature of single scanned outline， an algorithm to extract the down edge of the bulk cargo effectively was proposed to remove the measurement error from ground point cloud and a projection triangulation volume algorithm was used to calculate the volume of bulk’s point cloud. Experiment shows that the system could complete the measurement of bulk cargo with a size of 69 m3in 30 s， the average relative error is 0.42% and the repeated measurement error is 0.41%. In a reality bulk measurement， the measurement on about 31 500 m3takes less than 10 min， and the average repeated measurement error for 4 coal bulks is 0.74%. It shows that the system could measure the volume of bulk cargo quickly and easily on maintaining accurate measurement.

laser measurement; volume measurement; rapid measurement; large bulk; point cloud edge; triangulation

2016-04-05；

2016-06-03.

國(guó)家科技支撐計(jì)劃資助項(xiàng)目(No.2014BAB12B03)；國(guó)家自然科學(xué)基金資助項(xiàng)目(No.61501077)；中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目(No.3132016351)

1004-924X(2016)09-2126-08

TP23；TH821

A

10.3788/OPE.20162409.2126

楊德山(1987-)，男，河北石家莊人，2011年于江西理工大學(xué)獲得學(xué)士學(xué)位，2013年于大連海事大學(xué)獲得碩士學(xué)位，主要從事光電檢測(cè)、激光測(cè)量的研究。E-mail：67186486@qq.com