譚金,杜超,李彤
(中國航天科工集團七八〇一研究所,湖南長沙410205)
目前,以系留氣球和自控飛艇為代表的飛艇以其耗能少、滯空時間長、載重量大等優(yōu)點廣泛用于軍用、民用等領(lǐng)域。近年來,以具有3 km滯空高度能力的系留氣球和平流層自控飛艇浮空平臺為代表的大型無人飛艇項目更是成為研究的熱點,具有廣闊和良好的應(yīng)用前景。大型無人飛艇實際可升空高度和滯空時間與飛艇氣囊的充氣量大小密切相關(guān),因此,為保證飛艇的升空與回收安全,必須對氣囊的充氣量大小進行監(jiān)測,為放飛決策提供可靠的數(shù)據(jù)支撐。傳統(tǒng)的監(jiān)測手段都是以氣囊飽滿度定性估計,具有很大的主觀性和不確定性,所以對大型無人飛艇氣囊的體積定量測量就顯得非常迫切。
以激光掃描為代表的光學(xué)三維大尺寸測量技術(shù)與傳統(tǒng)的激光點對點的測距技術(shù)不同,激光掃描測量技術(shù)的發(fā)展為空間信息的獲取提供了全新的技術(shù)手段,由傳統(tǒng)的人工單點數(shù)據(jù)獲取轉(zhuǎn)變?yōu)檫B續(xù)自動數(shù)據(jù)獲取,提高了觀測的速度和準(zhǔn)確度,由于其融合了激光反射強度和物體色彩等光譜信息,可以真實描述目標(biāo)的整體結(jié)構(gòu)、形態(tài)特性以及光譜特征,具有測量范圍大、準(zhǔn)確度高、通用性強等特點,已成為大型飛行器、地形地貌、城市建筑三維重建等大尺寸物體幾何量測量的主要手段之一?;谌S激光測量的大型無人飛艇氣囊體積監(jiān)測系統(tǒng),是通過激光掃描獲得氣囊曲面點到激光掃描儀的距離,而后通過一系列的坐標(biāo)轉(zhuǎn)換、數(shù)據(jù)處理最終構(gòu)建氣囊的三維幾何模型,從而定量計算出氣囊的體積[1]。由于激光掃描獲取點云的速度較快,可以滿足對飛艇氣囊體積進行即時監(jiān)控的要求。
在飛艇氣囊底腹部中心位置安裝一個轉(zhuǎn)動能力不小于180°的云臺,具有180°扇區(qū)跨度掃描能力的二維激光陣列掃描儀裝在云臺上,實現(xiàn)對充氣后氣囊外形特征點的快速掃描,激光掃描系統(tǒng)結(jié)構(gòu)如圖1所示。設(shè)計總體路線是:艇載計算機對激光掃描原始數(shù)據(jù)包進行解算,轉(zhuǎn)化為三維坐標(biāo)體系,隨后通過內(nèi)插值、濾波技術(shù)重構(gòu)氣囊外形輪廓,最后通過積分獲得氣囊的體積。
如圖2所示,以掃描儀為原點O,囊體的平行切平面XbOYb為基準(zhǔn)面。云臺0時刻從零位線起,在設(shè)定的角速度ω下勻速轉(zhuǎn)動,考慮云臺零位線與掃描基準(zhǔn)線相差一個角度ψ0,則通過掃描基準(zhǔn)線的時刻為tb=ψ0/ω。從tb時刻開始采集數(shù)據(jù),每隔Δt(即每隔ωΔt的間隔角)對氣囊基準(zhǔn)面以上的斷面進行掃描,掃描儀按均分原理保留每個掃描斷面特征點到掃描原點的距離數(shù)據(jù),當(dāng)云臺工作時間達到tb+180/ω時,完成對基準(zhǔn)面XOZ以上的氣囊特征點的掃描,采集工作停止,云臺復(fù)位,等待下一個掃描采集指令。
圖2 掃描系統(tǒng)安裝位置與掃描方式
對于掃描斷面上的任意點Ai,假設(shè)云臺轉(zhuǎn)過掃描基準(zhǔn)線的角度為ψi,Ai所對應(yīng)的掃描斷面方位角為θi,到掃描原點的距離為di,則Ai點相對于坐標(biāo)系OXbYbZb的坐標(biāo)可以表示為
通過公式(1)即可把氣囊掃描斷面上的任意點的距離特征進行轉(zhuǎn)換,獲得基準(zhǔn)面以上的氣囊空間特征坐標(biāo)。若已知云臺安裝初始誤差角,則掃描點Ai的空間坐標(biāo)可按攝影測量原理中的共線方程計算:
式中:
式(2)中的矩陣分量均與坐標(biāo)系安裝初始誤差角有關(guān),通過式(2)變換可以消除系統(tǒng)安裝時引入的初始誤差角[2]。
無人飛艇氣囊體積監(jiān)測系統(tǒng)涉及的關(guān)鍵技術(shù)之一是掃描儀采集到氣囊外形特征點后,如何將已有特征點通過網(wǎng)格插值,重構(gòu)出氣囊的三維外形輪廓。本設(shè)計采用雙線性插值算法構(gòu)建三角網(wǎng)格結(jié)構(gòu),然后構(gòu)建計算網(wǎng)格,對每個計算網(wǎng)格點在三角網(wǎng)格結(jié)構(gòu)中進行搜尋插值,通過查找均分點位于哪個三角形中來構(gòu)建其高程差值,獲得網(wǎng)格點整齊均分的計算網(wǎng)格坐標(biāo),最終構(gòu)建氣囊的三維特征外形。圖3為計算網(wǎng)格點在三角網(wǎng)格中的位置。為了確保計算網(wǎng)格點高程插值的一一對應(yīng)屬性,在對氣囊外形進行三維重建時采用了區(qū)域分塊技術(shù),把氣囊分為多個部分分別進行計算,最后通過面拼接將各部分體積累加即為總體積。
圖3 平面三角網(wǎng)與待插值點
由于無人飛艇氣囊為柔性囊體,飛艇氣囊體積監(jiān)測裝置使用時無法進行剛性固定,掃描儀在掃描時的晃動振動將形成散亂點或者空洞等雜波或噪聲,需要通過對點云數(shù)據(jù)進行去噪濾波,以保證原始數(shù)據(jù)點的平滑特性[3-4]。
假設(shè)某一斷面采集了n個數(shù)據(jù)點,當(dāng)對點Si((nj)>i>j)進行濾波時,先根據(jù)不同的測量環(huán)境選定參數(shù)值j(j的取值一般為2到5之間),求出Si及兩邊相鄰的各j個點Si-j,…,Si,Si+1,…,Si+j到激光器S0的距離Di-j,…,Di,Di+1,…,Di+j;而后對距離設(shè)定權(quán)值,對Di其權(quán)值設(shè)定為2j,余下的Dk(i-j≤k≤i+j)權(quán)值設(shè)為1,根據(jù)激光測量的環(huán)境不同,濾波閾值fd應(yīng)滿足
當(dāng)距離信號低于閾值時,將距離信號值濾除;利用中值濾波,剔除奇異點;利用曲面擬合去除前端遮擋物。通過上述處理方法,可以剔除原始點云數(shù)據(jù)中的錯誤點和含有粗大誤差的點。
控制及采集系統(tǒng)由激光掃描儀、小型云臺、串口/以太網(wǎng)信號轉(zhuǎn)換器、艇載以太網(wǎng)交換機、光纖收發(fā)器和地面測控計算機組成,為避免定位誤差的累積,每個掃描周期后云臺都將復(fù)位至原始位置,小型云臺在水平面從0°轉(zhuǎn)動至200°再復(fù)位至0°的時間為一個掃描周期,掃描周期T0的值隨著云臺的水平轉(zhuǎn)速的大小而變化,云臺的水平轉(zhuǎn)速可通過地面測控計算機上的云臺控制軟件來設(shè)定??紤]到測量誤差,舍棄云臺(0°,10°)和(190°,200°)兩個不勻速運動的區(qū)間,只選取掃描周期中云臺轉(zhuǎn)速均勻的中間段(10°,190°)進行采集,采集角度范圍依然保持為180°。地面測控計算機的采集頻率根據(jù)采集周期和最小采樣角度來確定,本方案中采樣角度為1°~5°,采集頻率f與掃描周期T0之間的關(guān)系為
式中:α為采樣角度范圍;β為最小采樣角度。
采用某型飛艇氣囊對無人飛艇氣囊體積監(jiān)測裝置進行測試驗證,無人飛艇氣囊體積監(jiān)測裝置對氣囊進行激光掃描后,通過坐標(biāo)轉(zhuǎn)換、數(shù)據(jù)處理重構(gòu)出氣囊的三維特征外形如圖4所示。
經(jīng)與結(jié)構(gòu)設(shè)計工程師確認(rèn),氣囊充氣飽和后的體積設(shè)計理論值為234 m3,對氣囊連續(xù)進行6次測量,測試結(jié)果如表1所示,6次測量均值為235.08 m3,重復(fù)性為0.47 m3,實測均值與設(shè)計理論值相對誤差為0.46%,相對誤差控制在±1%以內(nèi),單次測量時間小于10 s,無人飛艇氣囊體積監(jiān)測裝置可以滿足對囊體體積即時測量的要求。
圖4 軟件重構(gòu)測量結(jié)果界面
表1 氣囊即時檢測值 m3
本文將三維激光掃描技術(shù)引入無人飛艇氣囊體積測量,通過激光掃描獲取氣囊點云數(shù)據(jù),采用雙線性插值算法構(gòu)建三角網(wǎng)格,最終重構(gòu)氣囊的三維幾何模型,實現(xiàn)了對無人飛艇氣囊體積即時監(jiān)測,實驗結(jié)果表明,體積測量相對誤差在±1%以內(nèi)。另外,針對現(xiàn)場實際使用情況還可以做以下幾方面的改進:①對于采樣點擬合表面與真實表面的誤差,以及測量過程中氣囊外形由于外界干擾所產(chǎn)生的形變誤差,可采用k-最近鄰方法構(gòu)造氣囊表面的三角形網(wǎng)絡(luò),結(jié)合設(shè)計外形以及制造過程中引入的常見外形偏差及分布的分析,獲得每個三角形的先驗曲面形式,最終重建三維表面;②對于氣囊飽和度較小時產(chǎn)生的掃描盲區(qū)所引入的擬合處理誤差,可配合艇內(nèi)攝像機來綜合評定特征外形;③通過優(yōu)化軟件算法進一步縮短單次測量的時間,提高系統(tǒng)測量準(zhǔn)確性。
[1]李欣,周佳瑋.三維激光掃描技術(shù)在船體外形測量中的試驗性研究[J].測繪信息與工程,2006(31):36-37.
[2]王昌翰,向澤君,劉潔.三維激光掃描技術(shù)在文物三維重建中的應(yīng)用研究[J].城市勘測,2010(6):67-70.
[3]嚴(yán)劍鋒,鄧喀中,邢正全.基于最小二乘擬合的三維激光掃描點云濾波[J].測繪通報,2013(5):43-46.
[4]胡向紅,陳康寧.激光掃描數(shù)據(jù)的脈沖噪聲過濾[J].計算機輔助設(shè)計與圖形學(xué)學(xué)報,2003(9):1056-1059.