結合光達資料與航空影像重建屋頂面之研究

(1)

結合光達資料與航空影像重建屋頂

面之研究

指導教授：趙鍵哲　　　　　　　　學生姓名：鄭傑中

(2)

大綱  前言  文獻  研究理論  模擬實驗  實際資料  問題與討論  參考文獻

(3)

前言

 目前所熟知的光達系統已可快速且精確地獲得地表三維資訊。然而，僅由光達系統來獲得三維資訊可能會產生離散點解析度不足的問題，且光達資料在平面精度表現上比攝影測量差，因此，在應用上仍有其侷限。相反的，影像具有平面精度的優勢與光達所欠

缺的光譜資訊、 Scene information [Schen

k & Csatho, 2002] ，若能結合光達原有的高程資訊，將相輔相成。

(4)

文獻回顧

 資料融合：

 賴彥中 [2004] ，從資訊融合角度出發，結合光達資料

及數位化彩色空照影像進行三維建物模型之重建。

 Schenk, T[ 2002] 　從 Lidar 角度出發將 Lidar 資料結

構化，並加入影像的平面資訊使其獲得地表更完整的資訊。

(5)

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

(6)

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

(7)

研究理論

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

LIDAR Aerial imagery

Advantages _{1.高程精度較佳}

2.waveform analysis

3.It’s not perspective projection ， range image is similar such as ortho-photo

1.rich in scene information 2.平面精度高

Disadvantages 1.no scene information

2.horizontal accuracy? 1. Stereo matching (occluded areas)2. Degree of automation 3. Because of perspective

(8)

Lidar 前資料處理

 由於 Lidar 原始資料為離散型態，為了方便後續做影像處理的應用，首先將原始不規則 Lida r 資料利用 Sufer 6.0 軟體內插成規則的網格點，以高程值設定為規則網格的儲存值當作灰階的變化，如此可將此網格化資料製作成為高程影像 (Range image) ，再利用邊緣線偵測的技術萃取出 range image 的邊緣線，此即為高程變化較明顯的區域。

(9)

(10)

空間幾何原理

 平面和平面的關係  兩平面的夾角  兩平面的交一直線 L1  直線與直線的關係 ( 平行 , 垂直 , 交點 )  點與直線的關係

(11)

線條為基礎的攝影測量

 計算步驟如下：  Step1: 　利用像點觀測量組成判斷面。  Step2: 　兩判斷面進行外積獲得空間直線方向向量。  Step3: 計算三維直線和 XY, YZ ,XZ 平面的夾角，以最小者當做　　　　最佳幾何交會面，決定穿刺點。  Step4: 再計算投影中心和相片觀測值所組成的光束與三維質線的交會點，進而獲得三維點座標，因此由圖二示意圖可看出一條空間直線上會產生四個點位座標值。

(12)

(13)

模擬實驗

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

(14)

(15)

以線條為基礎的前方交會求解精度 0 5 10 15 20 25 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5x 10 4 0 -->90 [ ] 判斷面交角由度度取樣點數 Z [ :m ] 地面方向誤差變化量單位

(16)

組數角度 X (m) Y(m) Z(m) 1 0 度 0 分 29.47 秒 26600 8.8161 49525 2 0 度 4 分 54.66 秒 98.975 0.40409 319.4 3 0 度 24 分 33.18 秒 15.754 0.35842 54.577 4 0 度 49 分 5.52 秒 7.9127 0.35753 27.206 5 2 度 26 分 49.81 秒 2.7338 0.35795 9.097 6 4 度 3 分 15.04 秒 1.7034 0.35936 5.4997 7 6 度 23 分 39.32 秒 1.1282 0.36282 3.5017 8 7 度 53 分 25.69 秒 0.93829 0.36575 2.849 9 9 度 19 分 28.92 秒 0.81298 0.36929 2.4214 10 11 度 20 分 43.1 秒 0.6899 0.37559 2.0045 11 14 度 20 分 6.66 秒 0.57041 0.3882 1.6037 12 18 度 55 分 28.78 秒 0.45812 0.42029 1.2382 13 22 度 0 分 6.24 秒 0.40808 0.4598 1.0794 14 29 度 44 分 41.67 秒 0.277 0.40021 0.80679 15 31 度 36 分 27.11 秒 0.27295 0.39867 0.75634 16 36 度 01 分 38.67 秒 0.2656 0.39953 0.66386 17 41 度 38 分 0.88 秒 0.25935 0.408 0.58349 18 48 度 48 分 50.83 秒 0.25404 0.42968 0.51452 19 57 度 59 分 40.82 秒 0.24984 0.48091 0.45938 20 69 度 26 分 38.47 秒 0.24679 0.61975 0.41958 21 75 度 57 分 49.78 秒 0.24568 0.80251 0.40701

(17)

實際資料實驗

 為了讓整套作業流程更加自動化，因此當輸入房屋區塊時首先進行地形標準偏差計算，根據常理，若是表面平坦的平面屋頂，地形變化度不大，若地形標準偏差大表示該平面高程變化大，此時有可能為平面屋頂上有高塔或者為斜面屋頂，需再做進一步分析。  設定地形複雜因子為 1m

(18)

 一般對於斜屋頂的認知有以下幾點： (a)以屋脊線為中心左右對稱。 (b)絕大多數和屋脊線平行的房屋邊線都是屬於比較長的那一邊。 (c)左右兩半面的斜率相近。 (d)沒有女兒牆。 (e)山形房屋的房屋區塊形狀大部分是規則四邊形的 ( 不是矩形就是正方形 ) 。

(19)

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

(20)

門檻值的給定

 影像上 Buffer 的大小與角度值 :  原始 Lidar 掃描時的誤差  內插網格化的影響量  獲得 Lidar 線條的誤差  外方位誤差

(21)

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

(22)

平面屋頂物空間的篩選流程高程門檻距離門檻角度門檻場景分析 ( 加入圍牆厚度和圍牆高差門檻值 ) 女兒牆內外側的線條以和 Lidar 距離最接近的影像線條當作最佳解圍牆厚度門檻值設定為： 0.2~ 0.8m 　圍牆高差門檻設定為 0.7m

(23)

實驗成果

_{( 實驗一 )}

(24)

第二步房屋區塊分析  地形複雜度因子計算值為 _{(+-)0.984478< 預設} 門檻值 ₁  判定為表面大致平坦的平面屋頂

(25)

最後結果

2.5443 2.5444 2.5445 2.5446 2.5447 2.5448 2.5449 2.545 x 105 2.741 2.741 2.741 2.7411 2.7411 2.7411 2.7411 2.7411 x 106 140 145 150

(26)

(27)

第二步房屋區塊分析

地形高度的中誤差為 0.8 03201< 1

(28)

2.5434 2.5435 2.5436 2.5437 2.5438 2.5439 2.544 x 105 2.7411 2.7411 2.7411 2.7411 2.7411 2.7411 2.7412 x 106 140 145 150

(29)

(30)

第二步房屋區塊分析

 地形高度的中誤差為 1.8

96835>1

(31)

第三部將房屋切半計算左右半部

的斜率值

 左半部斜率平均值 :  -0.3378  左半部斜率中誤差 :  0.1587  右半部斜率平均值 :  0.0093  右半部斜率中誤差 :  0.0131  設定門檻值為 0.05 兩部分斜率一者為正另一為負，就常理來判不可能為山形屋，因此判定為平面房屋上有高樓 2.5494 2.5495 2.5495 2.5496 2.5496 2.5497 2.5497 2.5498 2.5498 x 105 2.7407 2.7407 2.7407 2.7407 2.7407 2.7407 2.7407 2.7407 2.7407 2.7407 x 106 155 160 165

(32)

Result

2.5494 2.5495 2.5495 2.5496 2.5496 2.5497 2.5497 2.5498 2.5498 2.5499 x 105 2.7407 2.7407 2.7407 2.7407 2.7407 2.7407 x 106 150 155 160

(33)

(34)

第二步房屋區塊分析  地形高度的中誤差為 1.5 34595>1  因此不可能為平坦的平面房屋，再做進一步分析為平面房屋有高樓，或者為斜面房屋

(35)

(36)

 左半部斜率平均值 :

 0.4141

 右半部斜率平均值 :

(37)

2.5483 2.5484 2.5485 2.5486 2.5487 2.5488 2.5489 2.549 2.5491 2.5492 2.5493 x 105 2.741 2.741 2.741 2.741 2.741 2.741 2.741 2.741 2.741 2.741 2.7411 x 106 150 155 160 2.5483 2.5484 2.5485 2.5486 2.5487 2.5488 2.5489 2.549 2.5491 2.5492 2.5493 x 105 2.741 2.741 2.741 2.741 2.741 2.741 2.741 2.741 2.741 2.741 2.7411 x 106

(38)

分成兩半開始重建

2.5483 2.5484 2.5485 2.5486 2.5487 2.5488 2.5489 2.549 2.5491 2.5492 2.5493 x 105 2.741 2.741 2.741 2.741 2.741 2.741 2.741 2.741 2.741 2.741 2.7411 x 106 150 155 160

(39)

 _{經過四組實驗的測試發現，} _{Lidar 線條的確具有輔助過濾大} 部分線條的優點。  實驗中在 scene-information 中增加兩道決策的門檻值 ( 圍牆厚度門檻，圍牆高差門檻 ) ，斜屋頂則增加線條長度的約制，可以有效的獲得正確的房屋邊線。  _{房屋重建的精度取決於影像線條，當輸入的影像線條精度很} 差時，當然重建的精度也會不佳。  上述山形房屋有一側遺失了，原因在於影像的灰階反差太小，因此利用人工量測，但可能是量測誤差太大，被篩選系統刪除了

(40)

參考文獻

 賴彥中，「結合光達資料與數位空照影像重建三維建物模型」，國立中央大學土木工程研究所碩士論文，中壢， 2004 。  陳世崇，「以屋頂面法線為約制條件協助立體像對中建物邊緣線之篩選」，國立成功大學測量工程學系碩士論文 ,1999 。  郭子審，「空照立體像對中半自動化矩形人工建物重建之研究」，國立中央大學土木工程研究所碩士論文，中壢， 2001 。  王冠華，「多重疊航照立體對半自動房屋模型重建」，國立中央大學土木工程研究所碩士論文，中壢， 2003 。  李志宏 ,2002, 應用線特徵物求解攝影測量方位參數與物型重建 , 碩士論文