郭玉芳,陳 浩*,巨小文,任 福
(1.自然資源部測(cè)繪標(biāo)準(zhǔn)化研究所,陜西 西安 710054;2.自然資源部陜西測(cè)繪產(chǎn)品質(zhì)量監(jiān)督檢驗(yàn)站,陜西 西安 710054;3.武漢大學(xué) 資源與環(huán)境科學(xué)學(xué)院,湖北 武漢 430079)
激光雷達(dá)技術(shù)具有高精度、高自動(dòng)化、高效率等顯著優(yōu)勢(shì),同時(shí)點(diǎn)云數(shù)據(jù)蘊(yùn)含的豐富信息為L(zhǎng)iDAR數(shù)據(jù)的使用提供了更多的可能性[1-6],但是LiDAR點(diǎn)云數(shù)據(jù)處理的算法、相關(guān)軟件以及激光雷達(dá)產(chǎn)品的規(guī)范性等方面仍制約著LiDAR 技術(shù)大范圍的普及性應(yīng)用[7-8]。我國(guó)在激光雷達(dá)技術(shù)領(lǐng)域的發(fā)展和應(yīng)用起步較晚,美國(guó)等發(fā)達(dá)國(guó)家在激光雷達(dá)領(lǐng)域的發(fā)展和應(yīng)用研究已有近二十余年的成功經(jīng)驗(yàn),在LiDAR的標(biāo)準(zhǔn)化領(lǐng)域也制定發(fā)布了相關(guān)標(biāo)準(zhǔn)。本文在梳理美國(guó)激光雷達(dá)標(biāo)準(zhǔn)化背景和發(fā)展歷程的基礎(chǔ)上,以USGS 頒布的《激光雷達(dá)基礎(chǔ)規(guī)范》最新版本為例,從LiDAR數(shù)據(jù)采集和數(shù)據(jù)處理兩方面分析了機(jī)載LiDAR標(biāo)準(zhǔn)化的關(guān)鍵技術(shù)內(nèi)容和技術(shù)指標(biāo),探討了美國(guó)激光雷達(dá)新技術(shù)標(biāo)準(zhǔn)化的特點(diǎn)與優(yōu)勢(shì),對(duì)于我國(guó)測(cè)繪地理信息新技術(shù)標(biāo)準(zhǔn)的研發(fā)和制定具有一定的借鑒和參考價(jià)值。
2003 年美國(guó)攝影測(cè)量與遙感學(xué)會(huì)制定發(fā)布了《LAS 規(guī)范》(1.0 版本),該標(biāo)準(zhǔn)為用戶提供一種開放、通用的數(shù)據(jù)格式,極大地推動(dòng)了LiDAR數(shù)據(jù)的共享和應(yīng)用[9-10]。2009 年,美國(guó)地質(zhì)調(diào)查局(USGS)啟動(dòng)了國(guó)家空間項(xiàng)目(National Geospatial Program,NGP),進(jìn)行美國(guó)全境范圍的激光雷達(dá)數(shù)據(jù)獲取。NGP計(jì)劃管理的一個(gè)重要項(xiàng)目是美國(guó)三維高程計(jì)劃(3DEP)項(xiàng)目,其旨在獲取三維空間數(shù)據(jù),目標(biāo)是到2023 年底獲取全面覆蓋美國(guó)本土、夏威夷及其屬地的LiDAR 數(shù)據(jù)和覆蓋阿拉斯加地區(qū)的干涉合成孔徑雷達(dá)(InSAR)數(shù)據(jù)[11]。3DEP 項(xiàng)目LiDAR 數(shù)據(jù)采集與全國(guó)范圍的眾多數(shù)據(jù)生產(chǎn)商合作,出現(xiàn)了交付成果不統(tǒng)一、數(shù)據(jù)集難以兼容、整體項(xiàng)目的分析處理異常困難等問題[12],因而亟需制定一個(gè)基礎(chǔ)性的標(biāo)準(zhǔn)來規(guī)范美國(guó)全境激光雷達(dá)數(shù)據(jù)采集的最低技術(shù)指標(biāo)和成果交付的統(tǒng)一格式和要求,以保障NGP整體數(shù)據(jù)的質(zhì)量和統(tǒng)一性。
LiDAR技術(shù)應(yīng)用實(shí)踐的需求推動(dòng)了LiDAR標(biāo)準(zhǔn)化的發(fā)展,USGS 于2012 年8 月正式發(fā)布了《激光雷達(dá)基礎(chǔ)規(guī)范》(以下簡(jiǎn)稱LBS 規(guī)范)1.0 版本,系統(tǒng)全面地規(guī)范了機(jī)載LiDAR從數(shù)據(jù)采集、數(shù)據(jù)處理、精度檢驗(yàn)到成果交付的具體要求,為L(zhǎng)iDAR技術(shù)的生產(chǎn)實(shí)踐和推廣應(yīng)用提供了很好的指導(dǎo)作用。該標(biāo)準(zhǔn)第一個(gè)版本發(fā)布之后迅速被美國(guó)的許多州、縣以及行業(yè)企業(yè)和數(shù)據(jù)生產(chǎn)商所接受,一些國(guó)家還將該標(biāo)準(zhǔn)作為本國(guó)制定激光雷達(dá)規(guī)范的基礎(chǔ),成為了國(guó)際地理信息領(lǐng)域影響頗為廣泛的一個(gè)LiDAR基礎(chǔ)性標(biāo)準(zhǔn)[13-14]。
測(cè)繪新技術(shù)的標(biāo)準(zhǔn)化并不是一蹴而就的,通常需要隨著技術(shù)的發(fā)展和應(yīng)用需求不斷發(fā)展完善。美國(guó)激光雷達(dá)的標(biāo)準(zhǔn)化是隨著LiDAR 技術(shù)在NGP 應(yīng)用實(shí)踐,在近10 a 間不斷進(jìn)行更新修改完善的過程。USGS 自2012年8月發(fā)布LBS規(guī)范1.0版本之后,隨后不斷進(jìn)行了標(biāo)準(zhǔn)的更新和修訂,實(shí)現(xiàn)了對(duì)標(biāo)準(zhǔn)隨技術(shù)應(yīng)用實(shí)踐和需求的同步動(dòng)態(tài)更新,LBS規(guī)范最新版本是2021年6 月發(fā)布的,各版本的發(fā)布時(shí)間和主要變化內(nèi)容具體見表1。
表1 激光雷達(dá)基礎(chǔ)規(guī)范的版本簡(jiǎn)況
數(shù)據(jù)采集和數(shù)據(jù)處理是保障LiDAR數(shù)據(jù)質(zhì)量的關(guān)鍵,也是關(guān)系LiDAR技術(shù)應(yīng)用的重要環(huán)節(jié),同時(shí)亦是標(biāo)準(zhǔn)化的重點(diǎn)和難點(diǎn)所在。LBS規(guī)范對(duì)機(jī)載LiDAR的數(shù)據(jù)采集和數(shù)據(jù)處理進(jìn)行了較為詳細(xì)的規(guī)定,以下以LBS 規(guī)范的最新版本2021 Rev.A 為例,從數(shù)據(jù)采集和數(shù)據(jù)處理兩方面對(duì)機(jī)載LiDAR標(biāo)準(zhǔn)化的關(guān)鍵技術(shù)環(huán)節(jié)和技術(shù)指標(biāo)進(jìn)行具體分析。
激光雷達(dá)的數(shù)據(jù)采集在很大程度上受項(xiàng)目區(qū)的環(huán)境及許多其他因素的影響。為了保障3DEP 項(xiàng)目全國(guó)范圍內(nèi)激光雷達(dá)采集數(shù)據(jù)的一致性,LBS 規(guī)范主要規(guī)定了必須滿足的最低采集要求,包括采集區(qū)域與天氣條件、質(zhì)量等級(jí)、采集標(biāo)準(zhǔn)等作了相關(guān)要求,具體內(nèi)容如下:
1)采集區(qū)域和天氣條件。LBS規(guī)范規(guī)定數(shù)據(jù)采集的區(qū)域應(yīng)覆蓋感興趣區(qū)域并外擴(kuò)100 m 的緩沖區(qū)域;天氣條件要求在無云霧且無雪的情況下進(jìn)行,且地面尚未發(fā)生大面積洪水或其他類型的淹沒情況,植被覆蓋區(qū)應(yīng)首選在無葉條件下進(jìn)行數(shù)據(jù)采集。在植被有葉情況下的數(shù)據(jù)采集工作在簽發(fā)任務(wù)單或合同之前應(yīng)得到USGS相關(guān)部門的批準(zhǔn)。
2)質(zhì)量等級(jí)要求。LBS 規(guī)范將質(zhì)量等級(jí)(QL)分為4 級(jí),規(guī)定從QL3 到QL0 數(shù)據(jù)質(zhì)量越來越高,并對(duì)不同質(zhì)量水平級(jí)別做了相應(yīng)的質(zhì)量要求[15],規(guī)定最低可接受額度質(zhì)量水平為QL2 級(jí)別,同時(shí)USGS 鼓勵(lì)3DEP項(xiàng)目的數(shù)據(jù)商采集提供高質(zhì)量水平的數(shù)據(jù)。
3)額定點(diǎn)云間距要求:總體額定點(diǎn)云密度(aggregate nominal pulse density, ANPD)和總體額定點(diǎn)云間距(aggregate nominal pulse spacing,ANPS)是衡量激光雷達(dá)點(diǎn)云數(shù)據(jù)采集質(zhì)量的重要參數(shù)。ANPD 是指定區(qū)域內(nèi)單位面積點(diǎn)云的典型或平均密度,通常表示為每平方米的點(diǎn)云數(shù);ANPS 是激光雷達(dá)數(shù)據(jù)集中點(diǎn)云之間的典型或平均橫向距離,通常以米為單位。LBS規(guī)范對(duì)不同等級(jí)下的ANPD和ANPS作了具體的要求,見表3。其中,在QL2等級(jí)下,ANPD應(yīng)不低于每平方米2個(gè)點(diǎn);ANPS應(yīng)不大于0.71 m;另外LBS規(guī)范特別指出特定的項(xiàng)目要依據(jù)當(dāng)?shù)鼐唧w的地形和土地覆蓋條件,可能數(shù)據(jù)采集時(shí)需要更高的點(diǎn)云密度。
機(jī)載LiDAR點(diǎn)云數(shù)據(jù)具有海量、高冗余、點(diǎn)云密度不均、非結(jié)構(gòu)化等特點(diǎn),點(diǎn)云數(shù)據(jù)的處理是當(dāng)前研究的重點(diǎn)和熱點(diǎn)問題,也是LiDAR標(biāo)準(zhǔn)化的關(guān)鍵技術(shù)環(huán)節(jié)。LBS 規(guī)范對(duì)點(diǎn)云數(shù)據(jù)處理的基礎(chǔ)性內(nèi)容如數(shù)據(jù)基準(zhǔn)、精度檢驗(yàn)、數(shù)據(jù)分類等方面的內(nèi)容進(jìn)行了具體規(guī)定,主要技術(shù)要求如下:
1)數(shù)據(jù)基準(zhǔn)。LBS規(guī)范對(duì)數(shù)據(jù)基準(zhǔn)主要規(guī)定了點(diǎn)云的文件格式、時(shí)間基準(zhǔn)、空間基準(zhǔn)、坐標(biāo)參考系等內(nèi)容。
2)精度檢驗(yàn)。LBS規(guī)范規(guī)定在對(duì)點(diǎn)云數(shù)據(jù)進(jìn)行分類和生成衍生產(chǎn)品之前,應(yīng)對(duì)點(diǎn)云數(shù)據(jù)的絕對(duì)高程精度和相對(duì)高程精度進(jìn)行檢驗(yàn),并提供檢驗(yàn)過程的詳細(xì)報(bào)告。
相對(duì)高程精度指在沒有地面測(cè)量控制點(diǎn)的情況下,激光雷達(dá)數(shù)據(jù)集內(nèi)部的幾何質(zhì)量。單點(diǎn)的精度計(jì)算公式為:
式中,Range 為每個(gè)像素中最高和最低激光雷達(dá)點(diǎn)云數(shù)據(jù)之間的差異,單位為m;Slope為指該單元相對(duì)于其8 個(gè)相鄰單元的最大坡度,以小數(shù)表示,由每個(gè)單元的最低高程計(jì)算;Cellsize為指單元格的邊長(zhǎng),單位為米,用ANPS 表示,經(jīng)四舍五入后向上取整,然后再加倍;CEILING 是一個(gè)將ANPS 向上取整的函數(shù);1.414為計(jì)算像素對(duì)角線尺寸時(shí)的系數(shù)[13]。
最終精度應(yīng)以劃定樣本區(qū)域的面Shapefile文件形式進(jìn)行報(bào)告。以每個(gè)柵格單元為樣本值,計(jì)算樣本區(qū)域Z 方向的均方根差(RMSDz)值是否超過規(guī)定采集QL等級(jí)的限制。QL2級(jí)別下,單幅內(nèi)精度和跨幅重疊一致性的RMSDz限值分別為≤0.06 m和≤0.08 m。
絕對(duì)高程精度、激光雷達(dá)數(shù)據(jù)及其衍生DEM的絕對(duì)高程精度應(yīng)根據(jù)ASPRS(2014)進(jìn)行評(píng)估和報(bào)告。劃分為植被和非植被土地覆蓋類型分別評(píng)估其絕對(duì)高程精度。評(píng)估應(yīng)說明4 個(gè)絕對(duì)精度值:點(diǎn)云數(shù)據(jù)和DEM 各 自 的NVA (nonvegetated vertical accuracy,NVA)和VVA(vegetated vertical accuracy,VVA)。采用ASPRS 方法時(shí),不同的質(zhì)量等級(jí)(QL)對(duì)數(shù)據(jù)的最低NVA 和VVA 要求不同,且需要同時(shí)滿足NVA 和VVA 的最低要求值。在QL2 質(zhì)量水平下,NVA 和VVA 在95%的置信度下的要求值分別為≤0.196 m 和≤0.30 m[13,15]。
3)數(shù)據(jù)分類。LBS規(guī)范規(guī)定了激光雷達(dá)數(shù)據(jù)必需的、最低的分類方案。在已分類的LAS 交付成果中,除非被標(biāo)記為保留點(diǎn)外,任何點(diǎn)不允許分配為0 類。所有屬于最低分類方案且未被標(biāo)記為保留點(diǎn)的點(diǎn)云都應(yīng)進(jìn)行正確分類,可以根據(jù)特別的要求或目標(biāo)添加新的分類,但需要特別注意的是,為了保證分類的連續(xù)性,整個(gè)項(xiàng)目中點(diǎn)的分類要保持一致[13]。
與歐美等發(fā)達(dá)國(guó)家相比,我國(guó)在激光雷達(dá)技術(shù)的發(fā)展和應(yīng)用研究領(lǐng)域的起步較晚,LiDAR的標(biāo)準(zhǔn)化工作也相對(duì)滯后。原國(guó)家測(cè)繪地理信息局也曾頒布了少數(shù)激光雷達(dá)的相關(guān)行業(yè)標(biāo)準(zhǔn),例如CH/T 8024-2011《機(jī)載激光雷達(dá)數(shù)據(jù)獲取技術(shù)規(guī)范》、CH/T 8023-2011《機(jī)載激光雷達(dá)數(shù)據(jù)處理技術(shù)規(guī)范》等[14]。受限于我國(guó)標(biāo)準(zhǔn)化管理機(jī)制,這些標(biāo)準(zhǔn)在實(shí)踐中未能隨LiDAR技術(shù)發(fā)展應(yīng)用的變化而及時(shí)更新,標(biāo)準(zhǔn)的適用性也非常有限。而在LiDAR數(shù)據(jù)質(zhì)量檢驗(yàn)和分級(jí)、LiDAR數(shù)據(jù)分類等關(guān)鍵技術(shù)方面,目前我國(guó)仍缺失相關(guān)的標(biāo)準(zhǔn)。
隨著我國(guó)實(shí)景三維建設(shè)的推進(jìn)以及自動(dòng)駕駛對(duì)高精度地圖市場(chǎng)需求的增長(zhǎng),LiDAR的標(biāo)準(zhǔn)化需求也愈加迫切。因此,筆者通過對(duì)近10 年USGS 發(fā)布的《激光雷達(dá)基礎(chǔ)規(guī)范》7 個(gè)版本的跟蹤和分析研究,認(rèn)為USGS 在LiDAR 新技術(shù)標(biāo)準(zhǔn)的制定、更新以及應(yīng)用等方面具有顯著特點(diǎn),其標(biāo)準(zhǔn)化的模式和思路對(duì)于我國(guó)研制測(cè)繪新技術(shù)標(biāo)準(zhǔn)具有很好的啟示作用,具體特點(diǎn)分析如下:
1)更新快速及時(shí),靈活性好。從上文的分析可以看出,USGS 在2012 年8 月發(fā)布了LBS 規(guī)范1.0 版本,之后不斷更新修訂推出更高級(jí)版本,到2021年共發(fā)布了7個(gè)版本,且2020年之后版本的標(biāo)識(shí)也發(fā)生了變化,版本標(biāo)識(shí)按年和字母組合進(jìn)行標(biāo)識(shí),意味著該標(biāo)準(zhǔn)的更新將更加及時(shí)快速,或?qū)崿F(xiàn)動(dòng)態(tài)化更新。這既是LiDAR技術(shù)迅速發(fā)展和應(yīng)用的需求,也是新技術(shù)標(biāo)準(zhǔn)自身不斷完善的需求。這種靈活、快速、動(dòng)態(tài)化的更新模式對(duì)于我國(guó)測(cè)繪新技術(shù)標(biāo)準(zhǔn)的發(fā)展具有很好的借鑒意義。
2)適用性好,保持與現(xiàn)有相關(guān)標(biāo)準(zhǔn)的協(xié)調(diào)性。LBS 規(guī)范的相關(guān)內(nèi)容保持了與現(xiàn)有標(biāo)準(zhǔn)的高度協(xié)調(diào)性,例如點(diǎn)云數(shù)據(jù)格式采用美國(guó)攝影測(cè)量與遙感學(xué)會(huì)(ASPRS) 發(fā) 布 的《LAS 規(guī) 范》 version 1.4-R15(ASPRS,2011)中所規(guī)定的點(diǎn)云數(shù)據(jù)記錄格式中較為通用的數(shù)據(jù)格式;成果的精度分級(jí)保持了與DEM相關(guān)標(biāo)準(zhǔn)的一致;交付成果的元數(shù)據(jù)保持了與美國(guó)聯(lián)邦地理數(shù)據(jù)委員會(huì)(FGDC)發(fā)布的《數(shù)字空間元數(shù)據(jù)內(nèi)容標(biāo)準(zhǔn)》(1998)的一致,同時(shí)還增加了反映LiDAR 特色的其他元數(shù)據(jù)內(nèi)容。該標(biāo)準(zhǔn)保持與現(xiàn)有相關(guān)標(biāo)準(zhǔn)的協(xié)調(diào)一致性,大大提高了該標(biāo)準(zhǔn)的市場(chǎng)適用性。
3)可操作性強(qiáng),涵蓋了從數(shù)據(jù)采集、數(shù)據(jù)處理到成果交付的全過程。LBS規(guī)范是一個(gè)全面系統(tǒng)的Li-DAR 技術(shù)規(guī)范,涵蓋了從LiDAR 飛行數(shù)據(jù)采集的要求、天氣條件、數(shù)據(jù)格式等,到數(shù)據(jù)處理的點(diǎn)云分類、定位精度檢驗(yàn)、瓦片要求等,以及最終成果交付的內(nèi)容都進(jìn)行了具體詳細(xì)的規(guī)定。與其他一些LiDAR數(shù)據(jù)處理規(guī)范僅僅側(cè)重于從LiDAR 點(diǎn)云數(shù)據(jù)生成DEM 不同的是。該標(biāo)準(zhǔn)著重強(qiáng)調(diào)了原始激光點(diǎn)云的處理過程,這樣可以保障在USGS NGP項(xiàng)目獲得完整的數(shù)據(jù)源包括元數(shù)據(jù)、特征線、質(zhì)量信息、描述性文檔和其他一些輔助數(shù)據(jù),更好地支持LiDAR技術(shù)不同使用目的應(yīng)用。
4)成果按精度分級(jí),內(nèi)容豐富,保障了LiDAR數(shù)據(jù)的廣泛應(yīng)用。USGS 在LBS 規(guī)范中規(guī)定成果按精度分為4 個(gè)等級(jí),不同等級(jí)的質(zhì)量和技術(shù)指標(biāo)不同。成果類型豐富多樣,包括分類的點(diǎn)云數(shù)據(jù)、裸地球表面數(shù)字高程模型、全部的特征線及其相關(guān)元數(shù)據(jù),這種多元化豐富的數(shù)據(jù)成果既凸顯了LiDAR 技術(shù)的優(yōu)勢(shì),兼顧了不同部門對(duì)高程數(shù)據(jù)的不同需求,同時(shí)還實(shí)現(xiàn)了一次性數(shù)據(jù)的獲取,多部門數(shù)據(jù)的共享和重復(fù)利用,推動(dòng)了地理信息資源價(jià)值的提升[14],對(duì)于我國(guó)目前開展的實(shí)景三維中國(guó)建設(shè)具有很好的借鑒意義。
激光雷達(dá)作為近年來迅速發(fā)展起來的一門新興測(cè)繪技術(shù),具有高精度、高自動(dòng)化、高效率等傳統(tǒng)測(cè)繪無法比擬的優(yōu)勢(shì),在當(dāng)前我國(guó)智慧城市建設(shè)、資源調(diào)查監(jiān)測(cè)、基礎(chǔ)測(cè)繪等領(lǐng)域的應(yīng)用越來越廣泛。當(dāng)前隨著實(shí)景三維、新型基礎(chǔ)測(cè)繪、自動(dòng)駕駛高精度地圖等研究的開展,隨之LiDAR 的標(biāo)準(zhǔn)化需求也愈加迫切。本文在梳理分析美國(guó)LiDAR標(biāo)準(zhǔn)化背景和標(biāo)準(zhǔn)化歷程的基礎(chǔ)上,從數(shù)據(jù)采集和數(shù)據(jù)處理兩方面對(duì)機(jī)載Li-DAR 標(biāo)準(zhǔn)化的關(guān)鍵技術(shù)環(huán)節(jié)和技術(shù)指標(biāo)進(jìn)行具體分析,并探討了美國(guó)激光雷達(dá)新技術(shù)的標(biāo)準(zhǔn)化特點(diǎn)和優(yōu)勢(shì),對(duì)于我國(guó)測(cè)繪地理信息新技術(shù)標(biāo)準(zhǔn)的研發(fā)和制定具有一定的參考價(jià)值。