李小紅
(山西省測繪地理信息院第一測繪院,山西 太原 030002)
隨著城市的快速發(fā)展,能快速、高精度地提取建筑物的基本情況,對城市規(guī)劃和建設(shè)有極其重要的意義,對航空攝影測量產(chǎn)品的快速更新和地理國情數(shù)據(jù)建庫的更新以及對三維建模也有非常重要的意義。1∶10000DLG 房屋層生成現(xiàn)有的方法有三種:(1) 利用符合時(shí)相的大比例矢量圖縮編,此方法作業(yè)效率高,但可利用的大比例尺矢量數(shù)據(jù)成果較少;(2)利用航空影像立體測圖,勞動強(qiáng)度大,作業(yè)成本高;(3)利用正射影像提取建筑物輪廓線,在城區(qū)樓層較高區(qū)域建筑物投影差較大,精度不能保證,還有部分建筑物被遮擋及陰影等很多因素干擾, 從高分辨率衛(wèi)星影像和航空影像中提取建筑物輪廓會有好多被遮擋的建筑物的輪廓提不到的情況。
LiDAR(Light Detection And Ranging)是激光探測及測距系統(tǒng)的簡稱。LiDAR 系統(tǒng)是一種利用機(jī)載激光掃描(Airborne Laser Scanning,ALS)技術(shù)快速采集地表點(diǎn)的攝影測量裝置,可自動生成精確的數(shù)字表面模型 DSM(Digital Surface Model)。ADS100 是徠卡公司的一款采用線陣CCD 推掃式數(shù)字成像技術(shù),集成了全球衛(wèi)星定位系統(tǒng)(GPS)和高精度慣性測量系統(tǒng)(IMU),其成果可以直接用于測繪生產(chǎn)作業(yè),最大限度地減少了外業(yè)控制測量工作,提高了攝影測量的工作效率和成果質(zhì)量。
像素工廠與傳統(tǒng)的攝影測量軟件相比較, 其中一個(gè)優(yōu)勢就是高效自動大面積計(jì)算數(shù)字表面模型(DSM)。超高精度數(shù)字地面模型柵格大小能達(dá)到0.04m的空間DSM(依賴原始影像分辨率),在計(jì)算獲得的高程影像上噪聲通常比較少,適用于光學(xué)衛(wèi)星、雷達(dá)衛(wèi)星以及航拍影像,不同于傳統(tǒng)的基于相關(guān)算法的技術(shù),做完空三之后,像素工廠生成立體相對,根據(jù)前后視及下視的多重相關(guān),在并行節(jié)點(diǎn)上快速生成高分辨率的DSM。DAS100 航攝儀可以獲取高分辨率的影像,利用像素工廠生產(chǎn)高精度的DSM。
近年來,很多學(xué)者對利用點(diǎn)云數(shù)據(jù)提取建筑物輪廓進(jìn)行了大量的研究。如武漢大學(xué)的李必軍等提出的利用建筑物的幾何特征,從激光掃描數(shù)據(jù)中進(jìn)行建筑物特征提取研究;南京大學(xué)的程亮等提出在LiDAR 輔助下利用超高分辨率影像提取建筑物輪廓方法等。利用機(jī)載激光雷達(dá)獲取的點(diǎn)云數(shù)據(jù)或航空攝影測量處理計(jì)算生成的DSM 數(shù)據(jù)包括了地形和地貌的三維數(shù)據(jù),本文提出了通過DSM 數(shù)據(jù)提取建筑物的輪廓線的方法。
因?yàn)镈SM 是在DEM 的基礎(chǔ)上包括了建筑物、構(gòu)筑物、樹木、植被等要素,所以要提取建筑物的輪廓等信息,就需要首先提取出地物信息,通過機(jī)載激光雷達(dá)或航空攝影測量獲取的LAS 格式的DSM,對DSM 進(jìn)行濾波處理生成DEM,利用DSM 與DEM 的差值模型影像,再根據(jù)建筑物幾何形狀基本規(guī)則,折角一般接近直角,面積比較大,房頂高程基本一致等特性。根據(jù)高程閾值和面積閾值從地物中過濾掉植被等其他地物點(diǎn),生成建筑物輪廓線,將輪廓線轉(zhuǎn)點(diǎn),通過房子的一些特征,平行、垂直進(jìn)行建筑物輪廓線細(xì)化。流程(如圖1 所示):
圖1 技術(shù)流程
比周圍點(diǎn)的平均高程明顯高的點(diǎn)為空中點(diǎn),這些點(diǎn)通常是噪點(diǎn),低點(diǎn)通常指的是一個(gè)或一組比周圍點(diǎn)明顯低的點(diǎn),一般在地面點(diǎn)分類之前需進(jìn)行低點(diǎn)分類,而且需要分類多次才能全部找到低于地面的錯誤點(diǎn)。需要去除錯誤點(diǎn)和激光掃描儀的系統(tǒng)誤差可能造成不必要的高位或低位的粗差點(diǎn)云。首先,利用濾波軟件剔除這些粗差點(diǎn)云。剔除DSM 中的壞值,即極高、極低粗差點(diǎn)。然后,基于地形坡度學(xué)的分類方法,該方法根據(jù)假設(shè)法,假設(shè)兩個(gè)相鄰的點(diǎn)的高差不可能由地形陡坡引起,這些高程點(diǎn)更可能是非地面點(diǎn),兩點(diǎn)的高差隨著兩點(diǎn)距離縮小,較高的點(diǎn)為地面點(diǎn)的可能性降低。通過迭代法建立三角網(wǎng)的方式過濾地面點(diǎn),例如假設(shè)max building size = 70m,它則是假設(shè)在一個(gè)70m×70m 的區(qū)域內(nèi),至少有一個(gè)點(diǎn)為地面點(diǎn)。通過該方法找出一些初始的地面點(diǎn)構(gòu)成三角網(wǎng),再根據(jù)限制的條件(地形的角度、迭代角的大小、迭代距離),不斷地向三角網(wǎng)中加入LAS 點(diǎn)進(jìn)行迭代,進(jìn)一步細(xì)化三角網(wǎng),最終得到的三角網(wǎng)頂點(diǎn)則為地面點(diǎn)。根據(jù)地形類別的不同,選擇不同的參數(shù),一般城市地區(qū),Max building size 要設(shè)大點(diǎn),在山區(qū)或農(nóng)村,Max building size 要設(shè)小點(diǎn)。通過濾波后剩下的非地面點(diǎn)中還包含房屋、車輛、樹木、花草等地物。
由于地形的高低起伏太大會給提取建筑物輪廓帶來麻煩,灰度值較大的區(qū)域不一定是高建筑物,可能是其所在的地形較高。同時(shí),地形的起伏會給建筑物輪廓的提取算法增加了干擾,容易引起錯誤。因此,我們必須首先消除DSM 地形因子,提取數(shù)字地面模型(DEM),然后獲得規(guī)格化后的DSM(即nDSM),然后開始構(gòu)建提取。由于機(jī)載激光雷達(dá)或者航空攝影測量獲取的DSM 包括了地形和地貌,要提取建筑物的信息,建筑物屬于地物,首先要去掉地貌,用DSM 和DEM 做差值計(jì)算,生成地物影像。
利用邊緣分割的算法,對地物影像進(jìn)行分隔,這種算法計(jì)算很快,通過不同閾值控制,從而產(chǎn)生從粗到細(xì)的分隔。選擇高的閾值將會分隔出很少塊,選擇低的閾值將會分出很多塊,要選擇合理的閾值進(jìn)行分隔,盡量分隔出好的邊界。地物影像除了建筑物,還有植被、圍墻、道路等地物。利用建筑物比較高、面積比較大、屋頂高程一致的特性,提取建筑物影像。
(1)根據(jù)建筑物高度一般都高于3m,所以高程閾值h 設(shè)為3m,這時(shí)候把低矮植被、汽車等過濾掉,從地物影像里提取h>3m 的地物影像;余下的地物影像主要包括高植被和建筑物。
(2)檢測、標(biāo)記并計(jì)算連通地物影像的面積。
(3)根據(jù)建筑物的面積一般都大于6×6m2,設(shè)面積閾值S>36m2,可以過濾掉孤立的高植被、電線桿、路燈等地物。
(4)利用第三步提取的地物影像,生成二值圖像,如式(1):
式(1)p(xi,yj)表示點(diǎn)(xi,yj)高程值,一幅二值圖像按照從左到右、從上到下的順序搜索,找到的第一個(gè)g(x,y)=1 的點(diǎn)作為種子點(diǎn), 將種子點(diǎn)的8 鄰域范圍內(nèi)的高度與種子點(diǎn)的高度做差值,得出它與其連通地物影像每個(gè)像元的高差值Δhi,若所有的Δhi<0.5m,則該連通區(qū)域?yàn)榻ㄖ铩?/p>
(5)剩下的建筑物面邊緣雜亂,面中仍存在由數(shù)據(jù)處理中產(chǎn)生的小面積數(shù)據(jù)空白,建筑物面還存在不規(guī)則缺口和凸出塊的情況,為了處理小面積的空白和不規(guī)則的缺口和凸出塊,利用數(shù)學(xué)形態(tài)學(xué)中的閉運(yùn)算操作,閉運(yùn)算的原理是先膨脹后腐蝕的聯(lián)合運(yùn)算,閉運(yùn)算第一步需要先估算出空白數(shù)據(jù)的最大面積,然后用該面積除以單元格網(wǎng)的面積,取整數(shù)之后可得數(shù)學(xué)形態(tài)學(xué)閉運(yùn)算需要窗口B 的大小。數(shù)學(xué)形態(tài)學(xué)閉運(yùn)算可以填充比窗口B 小的缺口、孔洞和連同小間隔的間斷以及平滑邊界的功能,閉運(yùn)算可以部分彌補(bǔ)建筑物面域的空洞。閉運(yùn)算是利用結(jié)構(gòu)元素B 對初始提取的建筑物面A 進(jìn)行膨脹,然后接著再進(jìn)行一次腐蝕的結(jié)果,其公式如式(2):
式(2)中,A 為提取的建筑物面的集合;B 為窗口大小;·為閉運(yùn)算;為膨脹運(yùn)算;θ 為腐蝕運(yùn)算。
根據(jù)建筑物影像是二值圖,利用ArcGIS 矢量化建筑物邊線,先建一個(gè)面矢量文件和一個(gè)線矢量文件,設(shè)置坐標(biāo)系與影像坐標(biāo)系一致,開始編輯線矢量,在菜單欄上點(diǎn)擊“自定義”下的“工具條”,勾選上ArcScan 工具條,在ArcScan 工具條上點(diǎn)擊“矢量化”,選擇矢量化設(shè)置,在彈出的“矢量化設(shè)置”設(shè)置樣式為“輪廓”, 在ArcScan 工具條上點(diǎn)擊“矢量化”,選 “生成要素”,按照默認(rèn)設(shè)置生成建筑物矢量。生成的建筑物輪廓線有好多條線段組成。根據(jù)實(shí)際情況建筑物的邊數(shù)不會太多,建筑屋頂輪廓線存在平行、直角等特點(diǎn),需要規(guī)則化建筑屋頂邊線,修整過程包含兩步:(1)對所得結(jié)果中近似平行的一些直線進(jìn)行合并。(2)根據(jù)建筑物一般形狀是規(guī)則的幾何形狀,折角接近直角,把相交線正交化。
本文使用一種格網(wǎng)間距為0.2m-0.3m 的LiDAR點(diǎn)云數(shù)據(jù)和一種地面分辨率優(yōu)于0.2mADS100 航空影像數(shù)據(jù)。兩種數(shù)據(jù)是使用不同機(jī)載設(shè)備、采集相同區(qū)域得到的不同密度的點(diǎn)云數(shù)據(jù)。選取的該試驗(yàn)區(qū)覆蓋面積約為6.7km2,且地勢比較平坦,適合于采用本文的方法進(jìn)行實(shí)驗(yàn)。LiDAR 點(diǎn)云數(shù)據(jù)的密度是每平方米16點(diǎn), ADS100 航空攝影數(shù)據(jù)處理獲取的分辨率為0.4m的DSM 及LAS 格式的點(diǎn)密度是每平方米2.5 點(diǎn)。
首先將像素工廠生成的分辨率為0.4m 的DSM 數(shù)據(jù)轉(zhuǎn)換成LAS 格式,基于TerraScan 軟件剔除粗差點(diǎn)(及LAS 中的極低點(diǎn)和極高點(diǎn));然后,采用三角網(wǎng)迭代濾波法區(qū)分出地面點(diǎn)和非地面點(diǎn),利用二次樣條函數(shù)法從分離出來的地面點(diǎn)層內(nèi)插生成的格網(wǎng)間距為0.3m 和 0.4m 的 DEM,(如圖 2 所示);也利用二次樣條函數(shù)法將濾除低點(diǎn)的LAS 點(diǎn)內(nèi)插生成的格網(wǎng)間距為0.3m 和 0.4m 的 DSM,(如圖 3 所示)。從圖 3(a)和 3(b)中可以看出建筑物比較多,除了建筑物還有圍墻、車輛和植被等地物。再利用DSM1 減去DEM1 得到差值影像nDSM1(如圖 4(a)所示);利用 DSM2 減去 DEM2 得到差值影像nDSM2(如圖4(b)所示)。根據(jù)建筑物高度一般都高于3m,所以高程閾值h 設(shè)為3m,這時(shí)候把低矮植被、汽車等過濾掉, 從地物影像里提取h>3m 的地物影像;余下的地物影像主要包括高植被和建筑物,(如圖5 所示)。圖6 為采用面積閾值獲得,去掉低矮植物、圍墻、獨(dú)立樹的地物影像圖。圖7 為根據(jù)初步的建筑物影像圖生成二值圖,黑色值0 為非建筑物,白色值1 為建筑物。圖8 為提取的建筑物輪廓線。圖9 為對平行線近似合并及對建筑物線進(jìn)行正交化處理的結(jié)果。
為了分析提取建筑物輪廓線的精度,采用了兩種方法,一種方法是通過建立立體相對,利用ERDAS 軟件里的測圖模塊(LPS)采集建筑物矢量,將自動提取的建筑物邊界線和立體采集建筑物矢量進(jìn)行套合,結(jié)果顯示兩套矢量數(shù)據(jù)基本套合;另一種方法是利用提取的建筑物輪廓線和正射影像進(jìn)行套合,結(jié)果顯示也基本套合(如圖10 所示)。但是,由于高度閾值、面積閾值設(shè)置過大會出現(xiàn)個(gè)別漏提和高度閾值、面積閾值設(shè)置過小多提的情況。
本文提出了從激光雷達(dá)或從ADS 數(shù)據(jù)獲取的LAS 格式的數(shù)據(jù),通過濾波分離出地面點(diǎn)和非地面點(diǎn),做差值影像,分離出非地面點(diǎn),通過高度閾值、面積閾值和建筑物頂接近平行,分割出建筑物,自動提取建筑物的粗略邊線并進(jìn)行細(xì)化得到建筑物邊界線,能夠滿足1∶10000 DLG 生產(chǎn)的需要。
本文的方法仍然存在一些不足,實(shí)驗(yàn)中發(fā)現(xiàn)了分割建筑物時(shí),有少數(shù)建筑物和鄰近非建筑物分割成同一個(gè)建筑物面的情況,所以建筑物規(guī)則化處理算法有待更好地解決。今后將進(jìn)一步對利用LAS 數(shù)據(jù)結(jié)合正射影像提取建筑物輪廓線的方法和基于分類法從真正射影像中提取建筑物輪廓線的方法進(jìn)行研究。