结合光达资料与航空影像重建屋顶面之研究

結合光達資料與航空影像重建屋頂面之研究指導教授：趙鍵哲學生姓名：鄭傑中大綱前言文獻研究理論模擬實驗實際資料問題與討論參考文獻前言目前所熟知的光達系統已可快速且精確地獲得地表三維資訊。然而，僅由光達系統來獲得三維資訊可能會產生離散點解析度不足的問題，且光達資料在平面精度表現上比攝影測量差，因此，在應用上仍有其侷限。相反的，影像具有平面精度的優勢與光達所欠缺的光譜資訊、Sceneinformation[Schenk&Csatho,2002]，若能結合光達原有的高程資訊，將相輔相成。文獻回顧資料融合：賴彥中[2004]，從資訊融合角度出發，結合光達資料及數位化彩色空照影像進行三維建物模型之重建。

Schenk,T[2002]從Lidar角度出發將Lidar資料結構化，並加入影像的平面資訊使其獲得地表更完整的資訊。影像線條過濾

陳世崇[2000]，提出一套房屋邊緣線的過濾程序，目的在將不屬於建物邊緣線的灰值邊緣線過濾掉，為利用區域成長法將建物的頂部重建為一近似的平面，用以取得建物段線的近似位置，再將其反投影回影像上進行建物邊線篩選。三維模型重建

王冠華[2003]，以半自動的方式為目標，進行航空照片中建物三維模型之重建。其方法是利用數位航空照片進行建物的輪廓萃取，首先以人工點選建物之近似位置據以進行局部的線性特徵物萃取，再搭配特徵線段匹配，即可前方交會獲得三維線段，進而重建房屋三維模型

研究理論

廣義性的資料融合是利用各種不同工具，方法，將不同來源的資料找出之間的互補性並進行整合，以產生精度較高且內容更完整的資訊。依據融合理論，在進行系統整合前必須對兩種系統進行比較，找出雙方的互補性。在此概括比較這兩種不同Sensor系統間的差異。

Table1:比較兩種系統的優缺點LIDARAerialimageryAdvantages高程精度較佳waveformanalysisIt’snotperspectiveprojection，rangeimageissimilarsuchasortho-photorichinsceneinformation平面精度高Disadvantagesnosceneinformationhorizontalaccuracy?1.Stereomatching(occludedareas)2.Degreeofautomation3.Becauseofperspectiveprojection,reliefdisplacementoccurred.Lidar前資料處理由於Lidar原始資料為離散型態，為了方便後續做影像處理的應用，首先將原始不規則Lidar資料利用Sufer6.0軟體內插成規則的網格點，以高程值設定為規則網格的儲存值當作灰階的變化，如此可將此網格化資料製作成為高程影像(Rangeimage)，再利用邊緣線偵測的技術萃取出rangeimage的邊緣線，此即為高程變化較明顯的區域。相關的統計指標空間幾何原理平面和平面的關係兩平面的夾角兩平面的交一直線L1直線與直線的關係(平行,垂直,交點)點與直線的關係線條為基礎的攝影測量

計算步驟如下：

Step1:利用像點觀測量組成判斷面

。

Step2:兩判斷面進行外積獲得空間直線方向向量。

Step3:計算三維直線和XY,YZ,XZ平面的夾角，以最小者當做最佳幾何交會面，決定穿刺點。

Step4:再計算投影中心和相片觀測值所組成的光束與三維質線的交會點，進而獲得三維點座標，因此由圖二示意圖可看出一條空間直線上會產生四個點位座標值。前方交會問題中的奇異解模擬實驗

利用模擬實驗估計像點量測誤差和外方位參數誤差造成地面點位誤差的影響量，步驟如下：在0度到90度之間由觀測量組成判斷面，並計算判斷面夾角，夾角取樣數目總共21組。利用matlabrandom函式模擬3000組像點量測誤差(0.025mm)，為了符合實際情形，額外再加入隨機的外方位參數誤差於該21組取樣資料內。由於cubic內插方式產生的曲線較平滑故利用cubic內差的函式推估其他未取樣角度的地面誤差。隨機誤差模擬圖

以線條為基礎的前方交會求解精度

組數角度X(m)Y(m)Z(m)10度0分29.47秒266008.81614952520度4分54.66秒98.9750.40409319.430度24分33.18秒15.7540.3584254.57740度49分5.52秒7.91270.3575327.20652度26分49.81秒2.73380.357959.09764度3分15.04秒1.70340.359365.499776度23分39.32秒1.12820.362823.501787度53分25.69秒0.938290.365752.84999度19分28.92秒0.812980.369292.42141011度20分43.1秒0.68990.375592.00451114度20分6.66秒0.570410.38821.60371218度55分28.78秒0.458120.420291.23821322度0分6.24秒0.408080.45981.07941429度44分41.67秒0.2770.400210.806791531度36分27.11秒0.272950.398670.756341636度01分38.67秒0.26560.399530.663861741度38分0.88秒0.259350.4080.583491848度48分50.83秒0.254040.429680.514521957度59分40.82秒0.249840.480910.459382069度26分38.47秒0.246790.619750.419582175度57分49.78秒0.245680.802510.40701實際資料實驗

為了讓整套作業流程更加自動化，因此當輸入房屋區塊時首先進行地形標準偏差計算，根據常理，若是表面平坦的平面屋頂，地形變化度不大，若地形標準偏差大表示該平面高程變化大，此時有可能為平面屋頂上有高塔或者為斜面屋頂，需再做進一步分析。設定地形複雜因子為1m一般對於斜屋頂的認知有以下幾點：以屋脊線為中心左右對稱。絕大多數和屋脊線平行的房屋邊線都是屬於比較長的那一邊。左右兩半面的斜率相近。沒有女兒牆。山形房屋的房屋區塊形狀大部分是規則四邊形的(不是矩形就是正方形)

。輸入房屋區塊近似邊線的端點座標程式自動獲得內縮的邊線內的區塊資訊第一步：房屋區塊地形複雜分析複雜地形地形平坦執行平屋頂演算策略以邊長較長的該邊將房屋區塊切半，比較兩者斜率的平均值，和中誤差，若差異甚大代表為平頂房屋上有高塔，反之則判定為斜屋頂。門檻值的給定

影像上Buffer的大小與角度值

:原始Lidar掃描時的誤差內插網格化的影響量獲得Lidar線條的誤差外方位誤差物空間門檻值的給定

此部分根據判斷面交會的精度而定，另外由於Lidar和航測之間存有約2~3公尺左右的系統誤差，因此除了系統建議外也可由使用者自行設定將物空間門檻值稍加放寬。針對平屋頂在物空間的判別分別提出了兩道門檻值，其一為高差門檻(由實驗所獲得的經驗值給定)，再者為圍牆厚度門檻針對圍牆在影像的的寬度反推回地面上。平面屋頂物空間的篩選流程

高程門檻距離門檻角度門檻場景分析(加入圍牆厚度和圍牆高差門檻值)女兒牆內外側的線條以和Lidar距離最接近的影像線條當作最佳解圍牆厚度門檻值設定為：0.2~0.8m圍牆高差門檻設定為0.7m實驗成果(實驗一)第一步輸入資訊第二步房屋區塊分析地形複雜度因子計算值為(+-)0.984478<預設門檻值1判定為表面大致平坦的平面屋頂最後結果實驗成果(實驗二)第二步房屋區塊分析

地形高度的中誤差為0.803201<1判定為平屋頂實驗成果(實驗三)第二步房屋區塊分析地形高度的中誤差為1.896835>1可能為平屋頂山形房屋第三部將房屋切半計算左右半部的斜率值左半部斜率平均值:-0.3378左半部斜率中誤差:0.1587右半部斜率平均值:0.0093右半部斜率中誤差:0.0131設定門檻值為0.05兩部分斜率一者為正另一為負，就常理來判不可能為山形屋，因此判定為平面房屋上有高樓Result實驗區四第二步房屋區塊分析地形高度的中誤差為1.534595>1因此不可能為平坦的平面房屋，再做進一步分析為平面房屋有高樓，或者為斜面房屋將屋頂分成兩半進行分析左半部斜率平均值:0.4141右半部斜率平均值:0.06分成兩半開始重建經過四組實驗的測試發現，Lidar線條的確具有輔助過濾大部分線條的優點。實驗中在scene-information中增加兩道決策的門檻值(圍牆厚度門檻，圍牆高差門檻)，斜屋頂則增加線條長度的約制，可以有效的獲得正確的房屋邊線。房屋重建的精度取決於影像線條，當輸入的影像線條精度很差時，當然重建的精度也會不佳。上述山形房屋有一側遺失了，原因在於影像的灰階反差太小，因此利用人工量測，但可能是量測誤差太大，被篩選系統刪除了

參考文獻賴彥中，「結合光達資料與數位空照影像重建三維建物模型」，國立中央大