消費(fèi)級無人機(jī)傾斜攝影測量監(jiān)測棄土場的精度評估

胡晉茹,賴林楓,蘆治源,張曉峰,李 元,趙廷寧,魏廣闊

(1.交通運(yùn)輸部科學(xué)研究院,北京 100029; 2.交科院科技集團(tuán)有限公司,北京 100088;3.北京林業(yè)大學(xué)水土保持學(xué)院,北京 100083)

棄土(渣)場是生產(chǎn)建設(shè)項目產(chǎn)生的松散堆積體,對其開展監(jiān)測是預(yù)防和治理水土流失的基礎(chǔ),是國家生態(tài)文明建設(shè)的重要基礎(chǔ)支撐[1]。

近年來,無人機(jī)低空攝影測量技術(shù)快速發(fā)展,并在測繪地理信息、應(yīng)急保障、農(nóng)業(yè)、水利、電力、交通等民用領(lǐng)域得到廣泛應(yīng)用[2-4]。無人機(jī)攝影測量以其便捷靈活、快速高效的優(yōu)勢在棄土(渣)場監(jiān)測中發(fā)揮著重要作用[5]。無人機(jī)通過搭載傳感器快速而準(zhǔn)確地獲取高分辨率的影像,結(jié)合地面合理分布的像控點,運(yùn)用計算機(jī)手段提取地物觀測對象空間立體信息,如坡度、坡長、高程、位置、外觀特征等[6]。傾斜攝影測量作為一種新興的測量技術(shù)方法正在快速發(fā)展與推廣[7]。與傳統(tǒng)的豎直航空攝影測量技術(shù)不同,傾斜攝影測量不僅從豎直方向拍攝,而且通過調(diào)整鏡頭的方向以傾斜的方式進(jìn)行多角度攝影,采集更豐富的紋理信息,對地物信息表達(dá)更詳盡,在三維建模和測繪生產(chǎn)中展現(xiàn)出更加精細(xì)、準(zhǔn)確的潛力[8-10]。但搭載多鏡頭的傾斜攝影測量無人機(jī)價格昂貴,因此基于消費(fèi)級無人機(jī)的傾斜攝影測量技術(shù)應(yīng)運(yùn)而生[11-12]。消費(fèi)級無人機(jī)具有低價格、易操作、無須證照和申請空域等優(yōu)勢[13],其搭載的單鏡頭可通過云臺調(diào)整鏡頭俯仰角模擬多鏡頭無人機(jī)進(jìn)行傾斜攝影測量[14],從而構(gòu)建實景三維模型,在滿足精度的要求下,提高了效率并降低了經(jīng)濟(jì)與技術(shù)門檻。

目前針對消費(fèi)級無人機(jī)攝影測量的精度分析研究主要有:文獻(xiàn)[15]利用小型消費(fèi)級無人機(jī)進(jìn)行正射低空攝影測量,發(fā)現(xiàn)垂直方向的平均誤差是水平方向的8倍;文獻(xiàn)[16]基于消費(fèi)級無人機(jī)通過搭載傾斜相機(jī)進(jìn)行攝影測量,構(gòu)建了城市精細(xì)化的三維模型,利用檢查點進(jìn)行幾何精度評估,其中平面中誤差、高程中誤差均滿足《三維地理信息模型數(shù)據(jù)產(chǎn)品規(guī)范》(CH/T 9015—2012)中的Ⅰ級1∶500成圖比例尺測圖精度;文獻(xiàn)[17]基于消費(fèi)級無人機(jī)傾斜攝影測量構(gòu)建了精細(xì)三維模型并測算了研究區(qū)的土方量,與GNSS法對比具有可靠的精度。但是,2021年自然資源部發(fā)布了最新的《低空數(shù)字航空攝影測量內(nèi)業(yè)規(guī)范》(CH/T 3003—2021)[18],其中對攝影測量的精度有更高要求。現(xiàn)有對消費(fèi)級無人機(jī)攝影測量精度是否能滿足最新規(guī)范要求的研究較少。本文以公路建設(shè)項目產(chǎn)生的棄土場為對象,評估其在監(jiān)測過程中應(yīng)用消費(fèi)級無人機(jī)傾斜攝影測量的精度,并與傳統(tǒng)正射攝影測量方式進(jìn)行對比分析,為實際生產(chǎn)建設(shè)活動快速監(jiān)測提供參考。

1 試驗與方法

1.1 消費(fèi)級無人機(jī)傾斜攝影測量原理

基于消費(fèi)級無人機(jī)傾斜攝影測量的原理與專業(yè)傾斜航攝無人機(jī)相同[8,13-14],通過可旋轉(zhuǎn)俯仰角的云臺搭載單鏡頭模擬多鏡頭傾斜航攝儀。首先垂直于地面采集正射影像,以井字飛行方式從前、后、左、右4個方向與地面形成相同角度獲取斜片;然后對采集的影像通過尺度不變特征變換(SIFT)匹配算法進(jìn)行多視影像匹配提取特征點,采用光束法區(qū)域網(wǎng)平差嚴(yán)密解完成空中三角測量解算地面坐標(biāo),即基于共線方程與空間后方交會原理配合無人機(jī)獲取數(shù)據(jù)中的位置和姿態(tài)系統(tǒng)(position and orientation,POS)數(shù)據(jù),求出影像外方位元素和加密點物方坐標(biāo)[17,19],公式為

(1)

式中,(XA,YA,ZA)為地物點坐標(biāo);(XS,YS,ZS)為攝影中心在地面輔助坐標(biāo)系下的地物點坐標(biāo);f為焦距;λ為比例因子;ai、bi、ci(i=1,2,3)為影像的3個外方位角元素組成的9個方向余弦; (x,y)為像點的像平面坐標(biāo)。

對航測區(qū)域內(nèi)所有影像聯(lián)合為統(tǒng)一的平差區(qū)域,從而構(gòu)建特征點云,基于引入RTK測量得到的高精度地面像控點,消減POS數(shù)據(jù)誤差使得特征點云精度較高,從而提高空中三角測量解算的精度。基于空中三角測量特征點云成果進(jìn)行多視影像密集匹配生成超高密度點云,進(jìn)一步構(gòu)建三維模型。

1.2 試驗區(qū)概況

試驗區(qū)為湖南省邵陽縣白倉至新寧縣清江橋公路項目樁號K1+300(1#)、K10+600(2#)兩處棄土場,如圖1所示。原始地貌類型為低山丘陵間溝谷地貌。1#棄土場為坡地形,總占地面積為20 680 m2,場內(nèi)高程為295～310 m,總體地勢較為平緩,西南處略陡。2#棄土場為溝道形,總占地面積為14 653.3 m2,場內(nèi)高程為248～276 m,總體地勢復(fù)雜,由北至南高程逐漸降低,呈九級臺階地形分布。

圖1 研究區(qū)衛(wèi)星影像

1.3 試驗設(shè)備與軟件工具

采用DJI Phantom 4Pro消費(fèi)級四旋翼無人機(jī)獲取航圖,該無人機(jī)集飛行平臺、動力系統(tǒng)、控制系統(tǒng)、云臺相機(jī)于一體,配有遙控器與飛控軟件。其相機(jī)為FC6310單鏡頭,采用1英寸CMOS影像傳感器(1 inch=2.54 cm),有效像素為2000萬,鏡頭焦距為8.8 mm/24 mm(35 mm格式等效),云臺可控俯仰角范圍為-90°～+30°,控制精度為±0.02°,能夠在飛行中拍攝穩(wěn)定的圖像。

通過上海華測導(dǎo)航X9 RTK測量系統(tǒng)測量地面像控點的三維坐標(biāo),使用CORS模式時平面定位精度為8 mm+1×10-6D,高程定位精度為15 mm+1×10-6D,D為基站距離。使用與無人機(jī)POS數(shù)據(jù)相同的WGS-84坐標(biāo)系。

棄土場三維模型的構(gòu)建采用Metashape軟件(原Photoscan),該軟件基于最新的運(yùn)動恢復(fù)結(jié)構(gòu)——多視圖立體匹配(SfM-MVS)算法,對未經(jīng)畸變校正的影像和控制點進(jìn)行特征點提取與匹配、空中三角測量,自動化地完成相機(jī)參數(shù)解算與三維模型重建。

1.4 無人機(jī)航測方案設(shè)計

根據(jù)《低空數(shù)字航攝與處理規(guī)范》(GB/T 39612—2020)[20]中對航攝影像分辨率及航線重疊度的相關(guān)要求,綜合考慮地形條件、飛行效率,設(shè)置飛行高度相對測區(qū)最低點為100 m,航向重疊度為80%,旁向重疊度為70%。根據(jù)前人的經(jīng)驗[19],單鏡頭無人機(jī)傾斜攝影規(guī)劃5次航線,航線1云臺垂直正射采集正射影像,依次航線2、3、4、5調(diào)節(jié)云臺俯仰角為-45°,分別從東、南、西、北方向采集傾斜航片。在棄土場內(nèi)均勻布設(shè)地面標(biāo)志作為像控點,根據(jù)無人機(jī)航測的最優(yōu)控制點方案[21]選取四周和區(qū)域中心的6個像控點作為控制點參與空三解算,其他像控點作為精度評估的檢查點使用。各棄土場具體航測方案見表1。

表1 棄土場航測方案

1.5 精度評價指標(biāo)

目前,通常以像控點的誤差作為航空攝影測量的評價指標(biāo)。具體而言,基于航空攝影測量成果,識別出各像控點并提取其三維坐標(biāo),以RTK測量系統(tǒng)測量的三維坐標(biāo)為真值,對比分析得到各像控點在X、Y、Z方向的誤差,計算平面和高程的均方根誤差(RMSE)作為整體精度的評價指標(biāo)。根據(jù)《低空數(shù)字航空攝影測量內(nèi)業(yè)規(guī)范》(CH/T 3003—2021)[1]1∶500比例成圖中規(guī)定的誤差標(biāo)準(zhǔn)評價消費(fèi)級無人機(jī)航空攝影測量的精度。計算公式為

(2)

(3)

(4)

(5)

式中,Xi、Yi、Zi為第i個像控點的三維坐標(biāo);(XRTKi,YRTKi,ZRTKi)為RTK測量系統(tǒng)測量得到的第i個像控點的三維坐標(biāo);RMSEX,RMSEY、RMSEZ、RMSEH分別為X、Y、Z及高程的均方根誤差;n為像控點個數(shù)。

2 結(jié)果與分析

2.1 航測整體精度分析與評估

首先統(tǒng)計兩個棄土場兩種攝影測量方式下控制點和檢查點在X(經(jīng)度)、Y(緯度)、Z(高程)方向的誤差,生成各個方向上的誤差箱線(如圖2所示),對各攝影測量整體的精度進(jìn)行分析。由圖2(a)—(b)可知:①1#棄土場分別通過正射和傾斜攝影進(jìn)行空中三角測量時,檢查點與控制點在X、Y方向上的誤差分布范圍小且無較大差異;②在Z方向上,正射航測的檢查點比控制點誤差分布范圍顯著較寬,傾斜航測的檢查點誤差分布范圍顯著降低,略大于控制點的誤差分布范圍。這說明在1#棄土場中,傾斜航測相較于正射在平面位置的精度基本接近,而在高程上前者具有較高的精度提升。

圖2 各航測方式中控制點和檢查點在X、Y、Z方向上誤差的箱線

由圖2(c)—(d)可知:①2#棄土場在正射航測時,檢查點X方向上的誤差分布范圍在控制點X方向上的誤差分布范圍之內(nèi),而檢查點和控制點在Y、Z方向上的誤差分布范圍均較廣且具有差異。②2#棄土場在傾斜航測時,檢查點與控制點在X、Y、Z方向上的誤差分布范圍均較小且無明顯差異。這說明在2#棄土場中,正射攝影測量在平面位置和高程位置上的精度較低,而傾斜攝影測量成果精度能顯著改善。根據(jù)《低空數(shù)字航空攝影測量內(nèi)業(yè)規(guī)范》(CH/T 3003—2021)[1]1∶500比例成圖規(guī)定,以控制點合格數(shù)量和檢查點均方根誤差對各航測成果整體精度進(jìn)行評估,結(jié)果見表2。

表2 各航測方式整體精度評估

2#棄土場正射攝影測量的整體精度最低,其6個控制點中有1個不符合精度要求,檢查點平面位置的均方根誤差為11.1 cm,高程的均方根誤差為36.5 cm,超出了0.28 m的規(guī)定;其他航測成果的精度均符合規(guī)定的要求,2#棄土場傾斜攝影測量的整體精度最高,檢查點平面位置和高程的均方根誤差在5 cm內(nèi),相較于正射,檢查點平面位置均方根誤差低76.70%,高程均方根誤差低90.68%。1#棄土場正射攝影測量的檢查點平面位置均方根誤差為6.1 cm,高程均方根誤差為18.4 cm,而傾斜攝影測量的檢查點平面位置和高程的均方根誤差分布分別為4.6和7.4 cm,相較于正射,分別低25.19%和59.82%。與1#棄土場相比,2#棄土場正射航測的檢查點平面位置和高程的均方根誤差分別高82.31%和98.06%,傾斜航測的平面位置和高程的均方根誤差分別低43.21%和54.03%。

因此,傾斜攝影相比于正射攝影整體在平面和高程上具有不同程度的精度優(yōu)勢,在高程上更為明顯;2#棄土場相較于1#具有更復(fù)雜的地形條件,對正射攝影測量整體精度影響顯著。

2.2 檢查點精度分析與評估

為進(jìn)一步分析各航測成果的精度,對各成果內(nèi)檢查點的坐標(biāo)與RTK實測坐標(biāo)進(jìn)行對比,得到各航測成果檢查點的殘差,如圖3—圖5所示。

圖3 各航測成果檢查點誤差

圖4 各航測成果檢查點平面與高程誤差

圖5 各航測成果檢查點X與Y方向誤差

2#棄土場傾斜航測的各檢查點誤差集中分布在坐標(biāo)原點,1#棄土場傾斜航測次之;而2#棄土場正射航測的各檢查點誤差分散在外圍,1#棄土場正射航測的檢查點誤差分布在外圍的較少;正射攝影測量的檢查點在高程上的誤差要顯著大于平面的誤差。

根據(jù)《低空數(shù)字航空攝影測量內(nèi)業(yè)規(guī)范》(CH/T 3003—2021)[1]1∶500比例成圖規(guī)定,以平面和高程誤差的范圍和合格率對各航測方式檢查點的精度進(jìn)行評估,結(jié)果見表3。

表3 各航測方式檢查點精度評估 (%)

在平面誤差方面:各棄土場空中三角測量后檢查點平面誤差均符合規(guī)定的要求,1#棄土場通過正射攝影測量后平面誤差在10 cm以下的檢查點占94.12%,大于10 cm的占5.88%(僅1個),通過傾斜攝影測量后的檢查點平面誤差均小于10 cm;2#棄土場通過正射攝影測量的平面誤差在10 cm以下的檢查點占53.33%,大于10 cm的占46.67%,通過傾斜攝影測量后檢查點平面誤差均小于10 cm。

在高程誤差方面:1#棄土場通過正射攝影測量后高程誤差符合規(guī)定的檢查點占94.12%,小于10 cm的占52.94%,大于10 cm的占47.06%,通過傾斜攝影測量的檢查點高程誤差均符合規(guī)定的要求,小于10 cm的檢查點占70.59%,大于10 cm的占29.41%;2#棄土場通過正射攝影測量后高程誤差符合規(guī)定的檢查點僅占60.00%,小于10 cm的檢查點占33.33%,大于10 cm的占66.67%,通過傾斜攝影測量的檢查點高程誤差均符合規(guī)定的要求,小于10 cm的檢查點占93.33%,大于10 cm的占6.67%(僅1個)。

因此,各航測成果的檢查點誤差在高程位置要大于平面位置,這一差異在正射攝影測量中更為明顯;傾斜攝影測量的檢查點平面與高程精度均高于正射攝影測量,平面精度可由分米級小幅度提升至厘米級,高程精度大幅度提升至厘米級;復(fù)雜地形區(qū)域?qū)φ鋽z影測量的檢查點精度影響顯著,在高程方面更為明顯。

2.3 三維模型構(gòu)建

對各棄土場空中三角測量加密解算生成的特征點云成果通過多視圖立體匹配 (MVS)算法構(gòu)建三維模型,如圖6所示。由正射攝影測量構(gòu)建的三維模型在地形起伏區(qū)域噪聲點較多,甚至出現(xiàn)輕微扭曲、形變的情況(如圖6(a)和(c)所示),而由傾斜攝影測量構(gòu)建的三維模型則沒有出現(xiàn)上述情況(如圖6(b)和(d)所示)。對各棄土場空中三角測量提取的特征點進(jìn)行分析(如圖6(c)、(d)、(g)、(h)所示)可知,通過正射攝影三維重建提取的特征點較稀疏,1#和2#棄土場特征點密度分別為5.01和6.29個/m2;而通過傾斜攝影三維重建時特征點數(shù)量能顯著提升,1#和2#棄土場特征點密度分別為24.19和27.32個/m2。由此可知,在棄土場這種紋理單一、重復(fù)的地形條件下,通過正射空三加密解算后匹配的點云數(shù)量,不足以對棄土場的地形精確而細(xì)致的表達(dá),這導(dǎo)致三維模型出現(xiàn)較多噪聲點甚至扭曲、形變。傾斜攝影從豎直、傾斜等5個角度采集地物信息,因此可以提取足夠的特征點,使匹配的密集點云對地形進(jìn)行詳盡的表達(dá)。

圖6 各三維模型對比

3 結(jié) 論

本文量化分析了消費(fèi)級無人機(jī)在不同棄土場應(yīng)用正射和傾斜兩種攝影測量模式,以及根據(jù)攝影測量成果重建三維模型,結(jié)論如下:

(1)在丘陵地區(qū),無論是正射模式還是傾斜模式,消費(fèi)級無人機(jī)攝影測量成果的高程精度要低于平面的精度。正射攝影測量成果的精度受地形條件影響顯著,地形復(fù)雜、落差較大區(qū)域的精度不能滿足《低空數(shù)字航空攝影測量內(nèi)業(yè)規(guī)范》中1∶500比例的成圖精度要求。

(2)傾斜攝影能夠消減地形條件對測量成果的影響,相比于正射影像具有更高的精度優(yōu)勢,使平面精度由分米級小幅度提升至厘米級,高程精度大幅度提升至厘米級,完全符合《低空數(shù)字航空攝影測量內(nèi)業(yè)規(guī)范》中1∶500比例的成圖精度要求。

(3)消費(fèi)級無人機(jī)傾斜攝影測量從五個角度采集地物信息,構(gòu)建三維模型時提取特征點數(shù)量顯著提升,能對棄土場地形進(jìn)行精確而細(xì)致的表達(dá)。

本文數(shù)據(jù)采集及內(nèi)業(yè)處理耗時2 d,所用無人機(jī)系統(tǒng)價格約在萬元;消費(fèi)級無人機(jī)進(jìn)行傾斜攝影測量能節(jié)省大量成本、降低工作門檻,試驗評估了其在丘陵地區(qū)進(jìn)行1∶500比例成圖測量精度的可靠性,為公路建設(shè)項目棄土場快速監(jiān)測提供參考。