李道民,于海波
(寧波上航測繪有限公司,浙江 寧波 315200)
架站式三維激光掃描可以精準(zhǔn)掃描室內(nèi)外場景,但不便于房頂屋脊及建筑物細(xì)部的數(shù)據(jù)采集,應(yīng)用機(jī)載LIDAR、傾斜攝影測量可以有效解決架站式三維激光掃描儀掃測不到的房頂與屋脊點(diǎn)云數(shù)據(jù),應(yīng)用手持式SLAM 掃描可以有效解決架站式三維激光掃描儀與機(jī)載LIDAR 和傾斜攝影測量掃測不到的建筑物細(xì)節(jié)部分。本文就上述多種測量模式應(yīng)用于歷史建筑數(shù)字化營造進(jìn)行闡述,介紹空地多源數(shù)據(jù)的融合方法并對其精度進(jìn)行評定,對數(shù)據(jù)融合的難點(diǎn)及解決辦法進(jìn)行闡述。
歷史建筑數(shù)字化營造研究主要包括對歷史建筑內(nèi)部、外部結(jié)構(gòu)與周邊環(huán)境進(jìn)行復(fù)原和創(chuàng)建歷史建筑三維瀏覽動(dòng)畫?;谌S激光掃描、機(jī)載LIDAR、傾斜攝影測量、手持式SLAM 掃描等多種方式可對歷史建筑進(jìn)行點(diǎn)云獲取并建模,進(jìn)而對歷史建筑內(nèi)部與外部結(jié)構(gòu)及周邊環(huán)境進(jìn)行復(fù)原。但單項(xiàng)掃描歷史建筑技術(shù)手段獲取點(diǎn)云具有一定的局限性[1],不能將歷史建筑內(nèi)部結(jié)構(gòu)、外部結(jié)構(gòu)與周邊環(huán)境點(diǎn)云同時(shí)獲取,建模成果不佳,對歷史建筑數(shù)字化營造效果不理想。基于上述多種測量方式測得歷史建筑點(diǎn)云數(shù)據(jù),融合相關(guān)數(shù)據(jù)并對其建模從而對歷史建筑進(jìn)行數(shù)字化營造。數(shù)字化營造三維成果如圖1所示。
點(diǎn)云融合利用特征點(diǎn)匹配算法實(shí)現(xiàn)多源數(shù)據(jù)精確融合,點(diǎn)云融合有3 種融合類型:兩個(gè)點(diǎn)云均為相對坐標(biāo);兩個(gè)點(diǎn)云均為絕對坐標(biāo);兩個(gè)點(diǎn)云中其一為絕對坐標(biāo),其二為相對坐標(biāo)。
兩個(gè)均為相對坐標(biāo)的點(diǎn)云融合,與三維激光掃描單體架站掃描數(shù)據(jù)配準(zhǔn)原理同理,即利用特征點(diǎn)匹配算法來實(shí)現(xiàn)兩種點(diǎn)云精確融合,融合前提為兩個(gè)點(diǎn)云目標(biāo)具有明顯相同特征點(diǎn),相同特征點(diǎn)愈多融合匹配效果愈好,反之相同特征點(diǎn)愈少,點(diǎn)云融合匹配精度愈差,甚至匹配不上。
兩個(gè)均為絕對坐標(biāo)的點(diǎn)云融合,在3 種融合類型中最為簡便,只需將兩個(gè)點(diǎn)云共同導(dǎo)入點(diǎn)云處理軟件中,無須人工或電腦特征點(diǎn)匹配算法配準(zhǔn),即可自行匹配。兩個(gè)點(diǎn)云均為絕對坐標(biāo),融合前提為兩個(gè)點(diǎn)云需轉(zhuǎn)換為相同坐標(biāo)系和高程基準(zhǔn)下。點(diǎn)云融合成果如圖2所示。
點(diǎn)云數(shù)據(jù)為絕對坐標(biāo)與相對坐標(biāo)進(jìn)行匹配融合的情況,一般將相對坐標(biāo)點(diǎn)云配準(zhǔn)至絕對坐標(biāo)點(diǎn)云上。將相對坐標(biāo)點(diǎn)云配準(zhǔn)到絕對坐標(biāo)點(diǎn)云有兩種類型: 一種相對坐標(biāo)點(diǎn)云只有平面坐標(biāo)屬性沒有高程屬性,將該相對坐標(biāo)點(diǎn)云配準(zhǔn)到絕對坐標(biāo)點(diǎn)云上; 一種相對坐標(biāo)點(diǎn)云即無平面坐標(biāo)屬性亦無高程屬性,將該相對坐標(biāo)點(diǎn)云配準(zhǔn)到絕對坐標(biāo)點(diǎn)云上。第一種,將只有平面坐標(biāo)屬性沒有高程屬性的相對坐標(biāo)點(diǎn)云配準(zhǔn)到絕對坐標(biāo)點(diǎn)云上,配準(zhǔn)方法類似于兩個(gè)絕對坐標(biāo)點(diǎn)云的配準(zhǔn)。因兩套點(diǎn)云平面坐標(biāo)對應(yīng)相同,僅高程不對應(yīng),只需將相對坐標(biāo)點(diǎn)云垂直拖動(dòng)到絕對坐標(biāo)點(diǎn)云上,拖動(dòng)相對坐標(biāo)點(diǎn)云到特征點(diǎn)或特征面重合即可。第二種,既無平面坐標(biāo)屬性亦無高程屬性的相對坐標(biāo)點(diǎn)云配準(zhǔn)到絕對坐標(biāo)點(diǎn)云上,配準(zhǔn)方法類似于兩個(gè)相對坐標(biāo)點(diǎn)云的配準(zhǔn)。即采用特征點(diǎn)匹配算法來實(shí)現(xiàn)兩個(gè)點(diǎn)云精確融合,該配準(zhǔn)方法的前提也是兩個(gè)點(diǎn)云目標(biāo)具有明顯相同特征點(diǎn),相同特征點(diǎn)愈多融合匹配效果愈好,反之相同特征點(diǎn)愈少,點(diǎn)云融合匹配精度愈差,甚至匹配不上。當(dāng)然,如果兩個(gè)點(diǎn)云因重合特征點(diǎn)少,特征點(diǎn)匹配算法無法將兩個(gè)點(diǎn)云配準(zhǔn)一體,此時(shí)可采用另一種配準(zhǔn)方法。該配準(zhǔn)方法類似于將只有平面坐標(biāo)屬性沒有高程屬性的相對坐標(biāo)點(diǎn)云配準(zhǔn)到絕對坐標(biāo)點(diǎn)云上的配準(zhǔn),即在垂直拖動(dòng)的基礎(chǔ)上,同時(shí)將被配準(zhǔn)點(diǎn)云進(jìn)行旋轉(zhuǎn),將該點(diǎn)云通過平移與旋轉(zhuǎn)最終實(shí)現(xiàn)拖動(dòng)相對坐標(biāo)點(diǎn)云與絕對坐標(biāo)點(diǎn)云的特征點(diǎn)或特征面重合。將絕對坐標(biāo)點(diǎn)云配準(zhǔn)到相對坐標(biāo)點(diǎn)云上,其配準(zhǔn)方法類似于兩個(gè)相對坐標(biāo)點(diǎn)云的配準(zhǔn),即利用特征點(diǎn)匹配算法來實(shí)現(xiàn)兩點(diǎn)云精確融合,融合前提為兩個(gè)點(diǎn)云目標(biāo)具有明顯相同特征點(diǎn)。
空地多源點(diǎn)云數(shù)據(jù)的融合精度分析采用內(nèi)符合與外符合兩種方法。內(nèi)符合通過內(nèi)業(yè)對點(diǎn)云融合成果進(jìn)行切片,在切片中對比兩種測量模式的點(diǎn)云分層效果并量取分層距離,即點(diǎn)云融合后同一實(shí)際點(diǎn)位的垂直距離。內(nèi)符合主要對點(diǎn)云融合的相對精度和高程進(jìn)行分析,外符合通過外業(yè)測量歷史建筑特征點(diǎn)獲取絕對坐標(biāo),并在內(nèi)業(yè)融合成果中對該特征點(diǎn)進(jìn)行量取對比,外符合主要對點(diǎn)云融合的絕對精度和平面坐標(biāo)進(jìn)行分析。融合點(diǎn)云分層厚度量取如圖3所示。
空地多源點(diǎn)云數(shù)據(jù)的融合精度取決于兩個(gè)點(diǎn)云配準(zhǔn)方式,前文所述幾種配準(zhǔn)方式中,利用特征點(diǎn)匹配算法配準(zhǔn)方式的配準(zhǔn)精度相較于拖動(dòng)、旋轉(zhuǎn)等配準(zhǔn)方式要高,特征點(diǎn)匹配算法配準(zhǔn)方式精度取決于兩個(gè)點(diǎn)云重合特征點(diǎn)的多少[2]。采用拖動(dòng)與旋轉(zhuǎn)方式融合的點(diǎn)云成果以人工配準(zhǔn)為準(zhǔn),故點(diǎn)云配準(zhǔn)精度較特征點(diǎn)匹配算法略低。
在不考慮單項(xiàng)掃描歷史建筑技術(shù)手段點(diǎn)云拼接精度基礎(chǔ)上,即本文中點(diǎn)云精度評比只考慮單項(xiàng)掃描歷史建筑技術(shù)手段點(diǎn)云拼接后的整體點(diǎn)云與其他掃描點(diǎn)云技術(shù)手段獲得的整體點(diǎn)云之間的拼接精度。融合點(diǎn)云特征點(diǎn)量取如圖4所示。
經(jīng)實(shí)際案例比較,云配準(zhǔn)融合的精度評比如下:兩種絕對坐標(biāo)點(diǎn)云配準(zhǔn)融合精度> 將帶有平面坐標(biāo)屬性(無高程)的相對坐標(biāo)點(diǎn)云配準(zhǔn)到絕對坐標(biāo)點(diǎn)云的融合精度> 利用特征點(diǎn)匹配算法配準(zhǔn)兩種相對坐標(biāo)點(diǎn)云的融合精度> 利用拖動(dòng)、旋轉(zhuǎn)配準(zhǔn)兩種點(diǎn)云的融合精度。
以浙江省寧波市某歷史建筑的點(diǎn)云融合為案例,結(jié)合三維激光掃描與傾斜攝影測量獲取的兩種點(diǎn)云融合進(jìn)行分析,采用的點(diǎn)云處理軟件為Trimble Realworks。
為利于不同方式配準(zhǔn)的精度對比,本次架站式三維激光掃描儀獲取的點(diǎn)云分別輸出為絕對坐標(biāo)點(diǎn)云 (絕對坐標(biāo)點(diǎn)云的平面坐標(biāo)為CGCS 2000 坐標(biāo)系,高程基準(zhǔn)為正常高基準(zhǔn))與相對坐標(biāo)點(diǎn)云,傾斜攝影測量獲取的點(diǎn)云輸出平面屬性正常(CGCS 2000 坐標(biāo)系),高程屬性分別為大地高與正常高的兩種點(diǎn)云。本次點(diǎn)云配準(zhǔn)融合分別配準(zhǔn):(1)將三維激光掃描儀獲取的絕對坐標(biāo)點(diǎn)云與傾斜攝影測量獲取的平面屬性正常(CGC S2000 坐標(biāo)系)、高程屬性為正常高的點(diǎn)云配準(zhǔn)融合;(2)將三維激光掃描儀獲取的絕對坐標(biāo)點(diǎn)云與傾斜攝影測量獲取的平面屬性正常(CGCS 2000 坐標(biāo)系)、高程屬性為大地高的點(diǎn)云配準(zhǔn)融合;(3)將三維激光掃描儀獲取的相對坐標(biāo)點(diǎn)云與傾斜攝影測量獲取的平面屬性正常 (CGCS 2000 坐標(biāo)系)、高程屬性為正常高的點(diǎn)云配準(zhǔn)融合。上述3 種配準(zhǔn)方式旨在將兩種點(diǎn)云融合為絕對坐標(biāo)點(diǎn)云(平面坐標(biāo)為CGCS 2000坐標(biāo)系,高程基準(zhǔn)為正常高)。
1)將三維激光掃描儀獲取的絕對坐標(biāo)點(diǎn)云與傾斜攝影測量獲取的平面屬性正常(CGCS 2000 坐標(biāo)系)、高程屬性為正常高的點(diǎn)云配準(zhǔn)融合。該兩個(gè)點(diǎn)云為統(tǒng)一的平面坐標(biāo)系統(tǒng)和高程基準(zhǔn)的點(diǎn)云,即兩個(gè)絕對坐標(biāo)的點(diǎn)云配準(zhǔn)融合,在點(diǎn)云處理軟件Trimble Realworks 中,獲取帶有絕對坐標(biāo)的點(diǎn)云融合成果[3]。
2)將三維激光掃描儀獲取的絕對坐標(biāo)點(diǎn)云與傾斜攝影測量獲取的平面屬性正常(CGCS 2000 坐標(biāo)系)、高程屬性為大地高的點(diǎn)云配準(zhǔn)融合。該兩個(gè)點(diǎn)云中一個(gè)是絕對坐標(biāo)點(diǎn)云,另一個(gè)為僅有統(tǒng)一的平面坐標(biāo)、高程不統(tǒng)一的相對坐標(biāo),在點(diǎn)云處理軟件Trimble Realworks 中將相對坐標(biāo)點(diǎn)云垂直拖動(dòng)到絕對坐標(biāo)點(diǎn)云上,拖動(dòng)相對坐標(biāo)點(diǎn)云到特征點(diǎn)或特征面重合,獲取帶有絕對坐標(biāo)的點(diǎn)云融合成果。
3)將三維激光掃描儀獲取的相對坐標(biāo)點(diǎn)云與傾斜攝影測量獲取的平面屬性正常(CGCS 2000 坐標(biāo)系)、高程屬性為正常高的點(diǎn)云配準(zhǔn)融合。該兩個(gè)點(diǎn)云類型為相對坐標(biāo)點(diǎn)云配準(zhǔn)到絕對坐標(biāo)點(diǎn)云上,這里采用上節(jié)介紹的兩種配準(zhǔn)方法進(jìn)行配準(zhǔn)融合。
第一種配準(zhǔn)方式,采用特征點(diǎn)匹配算法來實(shí)現(xiàn)兩個(gè)點(diǎn)云精確融合,即與兩個(gè)相對坐標(biāo)點(diǎn)云的配準(zhǔn)融合同理,在融合過程中以傾斜攝影測量獲取的絕對坐標(biāo)點(diǎn)云為基準(zhǔn)點(diǎn)云,將三維激光掃描獲取的相對坐標(biāo)點(diǎn)云利用特征點(diǎn)匹配算法配準(zhǔn)到以傾斜攝影測量獲取的絕對坐標(biāo)點(diǎn)云上,即可獲得帶有絕對坐標(biāo)屬性的融合點(diǎn)云[4]。
第二種配準(zhǔn)方式,將相對坐標(biāo)點(diǎn)云通過平移與旋轉(zhuǎn)最終實(shí)現(xiàn)拖動(dòng)相對坐標(biāo)點(diǎn)云與絕對坐標(biāo)點(diǎn)云的特征點(diǎn)或特征面重合,獲得帶有絕對坐標(biāo)的融合點(diǎn)云。三維激光融合傾斜攝影點(diǎn)云數(shù)據(jù)前、后成果圖如圖5、圖6所示。
本次案例空地多源點(diǎn)云數(shù)據(jù)的融合的精度分析采用內(nèi)符合與外符合兩種方法。內(nèi)符合通過內(nèi)業(yè)對點(diǎn)云融合成果進(jìn)行切片,切片寬度截取0.1 m,對比三維激光掃描點(diǎn)云與傾斜攝影測量點(diǎn)云重合區(qū)域,對比分層情況,并量取該歷史建筑房頂兩點(diǎn)云分層厚度。
經(jīng)量取本案例點(diǎn)云分層厚度為:三維激光掃描絕對坐標(biāo)點(diǎn)云與傾斜攝影測量絕對坐標(biāo)點(diǎn)云融合點(diǎn)云成果分層厚度為0.006 m,即兩點(diǎn)云融合后高程誤差為0.006 m;三維激光掃描絕對坐標(biāo)點(diǎn)云與傾斜攝影測量相對坐標(biāo)點(diǎn)云通過垂直拖動(dòng)傾斜攝影測量點(diǎn)云的融合點(diǎn)云分層厚度為0.0173 m,即兩點(diǎn)云融合后高程誤差為0.0137 m;三維激光掃描相對坐標(biāo)點(diǎn)云與傾斜攝影測量絕對坐標(biāo)點(diǎn)云通過特征點(diǎn)匹配算法融合點(diǎn)云的分層厚度為0.037 m,即兩點(diǎn)云融合后高程誤差為0.037 m; 三維激光掃描相對坐標(biāo)點(diǎn)云與傾斜攝影測量相對坐標(biāo)點(diǎn)云通過拖動(dòng)、旋轉(zhuǎn)三維激光掃描相對坐標(biāo)點(diǎn)云的融合點(diǎn)云分層厚度為0.045 m,即兩點(diǎn)云融合后高程誤差為0.045 m。三維激光掃描點(diǎn)云數(shù)據(jù)與傾斜攝影點(diǎn)云數(shù)據(jù)融合切片成果如圖7所示。
外符合通過外業(yè)測量歷史建筑特征點(diǎn)獲取絕對坐標(biāo)并與內(nèi)業(yè)融合成果對該特征點(diǎn)進(jìn)行量取對比,本案例外業(yè)通過全站儀對歷史建筑特征點(diǎn)進(jìn)行點(diǎn)測量,同內(nèi)業(yè)分別量取的該特征點(diǎn)傾斜攝影測量點(diǎn)云、三維激光掃描點(diǎn)云和融合點(diǎn)云成果的坐標(biāo),本案例中不同配準(zhǔn)融合方式特征點(diǎn)坐標(biāo)對比如表1所示。由表1對比可得:三維激光絕對坐標(biāo)點(diǎn)云與傾斜攝影絕對坐標(biāo)點(diǎn)云融合后,特征點(diǎn)較全站儀測量特征點(diǎn)平面坐標(biāo)誤差≤0.01 m;三維激光絕對坐標(biāo)點(diǎn)云與傾斜攝影相對坐標(biāo)點(diǎn)云通過垂直拖動(dòng)相對坐標(biāo)點(diǎn)云融合后,特征點(diǎn)較全站儀測量特征點(diǎn)平面坐標(biāo)誤差為0.02~0.03 m;三維激光相對坐標(biāo)點(diǎn)云與傾斜攝影絕對坐標(biāo)點(diǎn)云通過特征點(diǎn)匹配算法融合后,特征點(diǎn)較全站儀測量特征點(diǎn)平面坐標(biāo)誤差為0.02~0.04 m;三維激光相對坐標(biāo)點(diǎn)云與傾斜攝影絕對坐標(biāo)點(diǎn)云通過拖動(dòng)、旋轉(zhuǎn)相對坐標(biāo)點(diǎn)云融合后,特征點(diǎn)較全站儀測量特征點(diǎn)平面坐標(biāo)誤差為0.04~0.08 m。
表1 不同配準(zhǔn)融合方式特征點(diǎn)坐標(biāo)對比m
本文就基于三維激光掃描、機(jī)載LIDAR、傾斜攝影測量、手持式SLAM 掃描等多種方式測量歷史建筑進(jìn)而對歷史建筑數(shù)字化營造的應(yīng)用進(jìn)行對比闡述,介紹多種測量方式獲得數(shù)據(jù)的融合方法并對其精度評比以及融合難點(diǎn)和解決辦法進(jìn)行闡述。通過具體案例給出了多種測量方式獲得的空地多源數(shù)據(jù)的融合精度評比結(jié)果:即兩種絕對坐標(biāo)點(diǎn)云配準(zhǔn)融合精度>將帶有平面坐標(biāo)屬性沒有高程屬性的相對坐標(biāo)點(diǎn)云配準(zhǔn)到絕對坐標(biāo)點(diǎn)云的融合精度> 利用特征點(diǎn)匹配算法配準(zhǔn)兩種相對坐標(biāo)點(diǎn)云的融合精度> 利用拖動(dòng)、旋轉(zhuǎn)配準(zhǔn)兩種點(diǎn)云的融合精度。空地多源點(diǎn)云數(shù)據(jù)的融合可大大提高歷史建筑數(shù)字化營造的數(shù)據(jù)完整性,可為歷史建筑數(shù)字化營造提供豐富的數(shù)據(jù)成果,為后續(xù)歷史建筑保護(hù)和利用提供可靠的數(shù)據(jù)支持。