曾靜靜，張曉剛，王剛

(1.青島市勘察測繪研究院，山東 青島 266032； 2.青島市西海岸基礎(chǔ)地理信息中心有限公司，山東 青島 266555；3.海陸地理信息集成與應(yīng)用國家地方聯(lián)合工程研究中心(青島)，山東 青島 266032)

1 引 言

近年來，我國社會經(jīng)濟(jì)快速發(fā)展，城市地表變化迅速且復(fù)雜，對基礎(chǔ)地理信息數(shù)據(jù)的現(xiàn)勢性要求也越來越高。而作為基礎(chǔ)地理信息數(shù)據(jù)的重要組成部分，地形圖及數(shù)據(jù)庫成果的快速準(zhǔn)確更新對城市管理具有重要意義。更新的關(guān)鍵之一在于對變化區(qū)域的檢測，準(zhǔn)確的變化檢測結(jié)果是合理選擇更新方法和更新模式的重要依據(jù)，根據(jù)整體變化檢測結(jié)果設(shè)計(jì)合理的更新方案，可以有效提高更新效率，確保成果質(zhì)量。

城市區(qū)域地表變化的檢測方法主要有：全野外巡圖、影像人工解譯對比、影像自動檢測變化以及利用已有資料等。其中，全野外巡圖方法準(zhǔn)確度高，但效率低、投入的人力和時間成本高；利用兩期影像目視解譯、人工勾繪變化區(qū)域的方法準(zhǔn)確度高，但是效率低；利用兩期影像自動檢測變化方法誤判率高；利用已有資料方法易產(chǎn)生變化丟漏，只能作為一種補(bǔ)充手段。以上幾種方法均有各種弊端，無法滿足快速、準(zhǔn)確識別變化區(qū)域的要求。

LiDAR技術(shù)融合了激光測量、飛行器姿態(tài)控制、GNSS等多種前沿技術(shù)[1，2]，能快速獲取地表形態(tài)，數(shù)據(jù)精度高、處理成本低，是一種新型的三維信息數(shù)據(jù)獲取手段。目前，利用LiDAR數(shù)據(jù)進(jìn)行地表變化檢測已有一定的研究成果：張良提出了基于多層局部ICP匹配的地形三維變化探測算法，通過差異信息熵判斷窗口變化，在最優(yōu)窗口下匹配，實(shí)現(xiàn)了基于多時相LiDAR數(shù)據(jù)的地表變化檢測[3]。鞏翼龍等人提出了一種利用多層次規(guī)則分類算法解決點(diǎn)云與影像建筑物變化檢測難題，形成了基于多源數(shù)據(jù)的建筑物變化檢測技術(shù)方法[4]。奚以成針對城市地表變化特點(diǎn)，重點(diǎn)研究了改進(jìn)型ICP算法的匹配、基于坡度熵的地表抽稀算法，并在此基礎(chǔ)上提出了變化檢測的策略[5]?？偟膩砜?，前人研究主要集中于多期點(diǎn)云數(shù)據(jù)的配準(zhǔn)算法，而針對城市區(qū)域地表變化檢測的實(shí)際應(yīng)用研究及生產(chǎn)總結(jié)較少。本文主要探索使用成熟的商業(yè)LiDAR處理軟件對城市區(qū)域進(jìn)行變化檢測，根據(jù)生產(chǎn)經(jīng)驗(yàn)總結(jié)出一套完整的數(shù)據(jù)處理流程，對實(shí)際工程應(yīng)用具有一定的指導(dǎo)意義。

2 研究區(qū)域和數(shù)據(jù)

區(qū)域1：研究區(qū)域總面積2 km2，屬于人工島嶼一部分(圖1)。2017年～2019年整個島嶼屬于建筑施工階段。該島嶼是人工填海而成，地形平坦，建筑物較多、植被較少，地面全部是硬化地表和鋪裝路面。

圖1 區(qū)域1范圍圖

區(qū)域2：總面積為5 km2，包含平地和丘陵地兩種地形，屬于城鄉(xiāng)接合部(圖2)。建筑物包括：農(nóng)村房屋、棚房、城市居民樓房、商場、工業(yè)廠房；道路包括城市道路、等級路；橋梁包括一般橋梁、立交橋、高架橋；測區(qū)內(nèi)有三條地面河流；植被以闊葉林和針葉林為主。實(shí)驗(yàn)區(qū)地形、地物豐富，具有代表性。

圖2 區(qū)域2范圍圖

本研究所用到的兩期LiDAR點(diǎn)云數(shù)據(jù)如表1所示：

研究區(qū)域數(shù)據(jù)參數(shù)列表 表1

3 技術(shù)方法

采用兩期機(jī)載LiDAR點(diǎn)云數(shù)據(jù)進(jìn)行快速變化檢測的主要流程是：

(1)數(shù)據(jù)預(yù)處理：分別對兩期點(diǎn)云數(shù)據(jù)進(jìn)行數(shù)據(jù)預(yù)處理，主要包括解算點(diǎn)云數(shù)據(jù)、剔除粗差點(diǎn)、坐標(biāo)轉(zhuǎn)換匹配、計(jì)算數(shù)據(jù)相對精度等。其中，相對精度的作用是確定變化閾值。在計(jì)算相對精度之前首先需根據(jù)研究區(qū)域的面積判斷是否需要將整個區(qū)域分成多個區(qū)塊，并分別計(jì)算各個區(qū)塊的相對精度，以防止區(qū)域過大導(dǎo)致整體相對精度與局部最優(yōu)相對精度差異較大，影響后續(xù)變化分析。

利用兩期影像數(shù)據(jù)，通過目視解譯法均勻地選取未變化的平坦的鋪裝路面以及硬化地表作為未變化區(qū)域，計(jì)算兩期數(shù)據(jù)中相對應(yīng)的未變化區(qū)域的高程差異均值作為該兩期數(shù)據(jù)的相對高程精度。經(jīng)過計(jì)算，區(qū)域1相對精度為 6 cm，區(qū)域2相對精度為 10 cm。

基于建筑物房屋角點(diǎn)、地物特征點(diǎn)等進(jìn)行局部平面位置匹配，采用局部四參數(shù)糾正方案，進(jìn)一步提高原始數(shù)據(jù)的相對一致性精度，使得后續(xù)變化分析更加準(zhǔn)確。

(2)數(shù)據(jù)分類：可將地物點(diǎn)云分為道路、水域、房屋、植被四大類。其中部分房屋比如尖頂房屋和屋頂結(jié)構(gòu)復(fù)雜的房屋，不同期的點(diǎn)云數(shù)據(jù)柵格化后高程差異較大，如圖3所示：A、B、C點(diǎn)和D、E、F點(diǎn)為不同時期的點(diǎn)云數(shù)據(jù)，雖然都在同一個柵格中，但很明顯用不同時期的點(diǎn)云數(shù)據(jù)得到的高程均值差異是比較大的，因此房屋的檢測閾值選取原則是數(shù)值盡量大，但要小于測區(qū)內(nèi)最低房屋的高度值。地面相對比較平緩，檢測閾值選取應(yīng)盡量小，因房屋和地面點(diǎn)所用的檢測閾值不同，需要將地面和房屋點(diǎn)區(qū)分出來；道路和河流都屬于地面點(diǎn)層，應(yīng)用同一檢測閾值；植被無法根據(jù)高度判斷變化，只可以判斷是否新增和消失。數(shù)據(jù)分類將點(diǎn)云數(shù)據(jù)分為：地面層、建筑物層、植被層。

(h)A+hB+hC)/3>(h)D+hE+hF)/3

(1)

注：h是點(diǎn)的高程值。

圖3 尖頂建筑物不同時期點(diǎn)云數(shù)據(jù)高程均值差值圖

首先對原始點(diǎn)云數(shù)據(jù)構(gòu)建不規(guī)則三角網(wǎng)，設(shè)置高度閾值，通過不斷的迭代計(jì)算各個點(diǎn)的高差與閾值的關(guān)系，判斷該點(diǎn)是否為地面點(diǎn)，從而將地面點(diǎn)與非地面點(diǎn)分開[6]。然后在非地面點(diǎn)云中通過判斷點(diǎn)云的平整度以及高度識別出建筑物屋頂區(qū)域。

(3)點(diǎn)云柵格化：分別對兩期建筑物、地面、植被點(diǎn)進(jìn)行柵格化處理，在柵格化過程中，兩期數(shù)據(jù)采用的網(wǎng)格大小需要相同，即保證兩期數(shù)據(jù)生成的柵格圖像的分辨率一致[7]。網(wǎng)格大小需根據(jù)點(diǎn)云密度選取，網(wǎng)格閾值Z由有下面的公式得出：

Z=(z1+z2+z3+…+zi)/i

(2)

其中Z表示每個格網(wǎng)點(diǎn)高程值，zi表示每個點(diǎn)的高程值，i表示每個柵格內(nèi)點(diǎn)的個數(shù)。

(4)變化檢測：對建筑物柵格圖像進(jìn)行聚類約束差值計(jì)算，首先按照空間鄰域高程相近原則進(jìn)行聚類平滑，降低孤立噪聲干擾；然后基于聚類后的柵格進(jìn)行比較，計(jì)算待分析格網(wǎng)到上一期同名格網(wǎng)的高程差值，得到提取變化圖斑后計(jì)算變化率，最后對結(jié)果進(jìn)行精度驗(yàn)證。

4 變化結(jié)果分析與評價

4.1 實(shí)驗(yàn)結(jié)果

房屋的檢測閾值選取原則是數(shù)值盡量大，避免尖頂房屋或者房屋細(xì)部結(jié)構(gòu)對結(jié)果造成的影響，但要小于測區(qū)內(nèi)最低房屋的高度值。本次實(shí)驗(yàn)區(qū)最低房屋的高度 2.3 m，因此選擇 2 m為建筑物檢測閾值，分類統(tǒng)計(jì)結(jié)果如圖4、表2所示：

圖4 建筑物變化檢測結(jié)果

建筑物變化結(jié)果統(tǒng)計(jì) 表2

對兩個不同時期地面圖像進(jìn)行計(jì)算并按照閾值進(jìn)行分割，主要分為地面升高區(qū)域、地面降低區(qū)域、未變化區(qū)域，結(jié)果如圖5、表3所示：

圖5 地形變化檢測結(jié)果

表3

4.2 分析與評價

(1)總體分析

測區(qū)1是填海造陸的人工島嶼，2017年～2019年間測區(qū)變化大。基于本文提出方法得出的變化檢測結(jié)果可知，整個測區(qū)大部分區(qū)域呈高度增加趨勢，與實(shí)際相符。根據(jù)建筑物變化檢測結(jié)果可以得出近年來該區(qū)域內(nèi)建筑物變化明顯，其中消失的建筑較多，多為建設(shè)工地中臨時搭建的棚房。

測區(qū)2是屬于城鄉(xiāng)接合部，2017年～2019年總體變化較小，變化區(qū)域多集中在新建小區(qū)、新建廠房多，拆遷村莊。經(jīng)過圖像對比驗(yàn)證，實(shí)驗(yàn)提取的變化圖斑，符合實(shí)際變化情況。

(2)準(zhǔn)確性分析

建筑物分析結(jié)果：基于兩期影像利用人工判讀的方法驗(yàn)證建筑物的變化檢測結(jié)果，所有變化情況全部符合實(shí)際情況，利用該方法不但能夠判斷出是否變化，而且還能進(jìn)一步分析出建筑物的新增、擴(kuò)建、拆除情況。其中，點(diǎn)云分類結(jié)果的正確性對變化檢測結(jié)果有影響，所以選擇合適的分類參數(shù)，保證分類的正確性對于建筑物變化檢測至關(guān)重要，實(shí)驗(yàn)區(qū)建筑物的正確率達(dá)到87%。

圖6 (a)建筑物旁存在安全網(wǎng) (b)建筑物安全網(wǎng)拆除

在研究區(qū)域中存在第一期數(shù)據(jù)中建筑物已經(jīng)封頂，安全網(wǎng)沒有拆除，而在第二期數(shù)據(jù)中該建筑物安全網(wǎng)拆除的情況，如圖6所示。此時由于建筑物高程未變化，本文方法將該建筑歸為未變化建筑物，而人工統(tǒng)計(jì)方法更加依賴影像紋理信息會將該建筑歸為變化建筑物。

地形分析結(jié)果：基于兩期影像利用人工判讀的方法驗(yàn)證地形的變化檢測結(jié)果，兩個實(shí)驗(yàn)區(qū)的地形變化情況全部符合實(shí)際情況。試驗(yàn)區(qū)準(zhǔn)確率達(dá)到76.5%，可以很好地檢測出地形變化情況，且目前濾波算法比較成熟，準(zhǔn)確率和流程自動化程度高。

水域、植被分析結(jié)果：不同時期河流、湖泊水位差異和不同季節(jié)植被長勢差異會對變化分析結(jié)果產(chǎn)生較大的影響，因此直接利用該方法求水域和植被的變化情況，準(zhǔn)確率低。

(3)效率分析

本文提出的方法簡單、自動化程度高，可以借助Terrasolid、LiDAR360等成熟軟件對點(diǎn)云進(jìn)行批處理，不需要過多的人工干預(yù)，降低了生產(chǎn)成本。通過本文提到的方法，節(jié)約了38%的人工工作量，具體分析如表4所示：

表4

以建筑物和地形變化檢測結(jié)果為依據(jù)選擇基礎(chǔ)地形圖和數(shù)據(jù)庫更新方式、更新機(jī)制、更新頻次，對于城市的基礎(chǔ)地理信息數(shù)據(jù)定期更新具有重要意義[8]。該方法已經(jīng)成功應(yīng)用于地形圖更新項(xiàng)目中。

5 結(jié) 論

本文利用兩期點(diǎn)云數(shù)據(jù)檢測城市地表變化，得到城市區(qū)域地表變化圖斑和變化率。所用方法準(zhǔn)確率較高、自動化程度高，能減少人工判讀的工作量，提高了地表變化檢測的效率，滿足了地形圖、數(shù)據(jù)庫更新前期了解區(qū)域變化情況的需求。但該方法也存在局限性，對植被、水域的變化檢測準(zhǔn)確率低于50%，不建議用該方法對植被和水域進(jìn)行變化檢測，因此如何提高變化檢測的全面性仍是后續(xù)研究的重點(diǎn)。