附录数值地形资料及卫星影像取得.docVIP

附錄一、數值地形資料及衛星影像取得 1.1 衛星影像 種類 解析度 圖幅大小 1.LANDSAT MSS 80 x 80 公尺 185 x 185 KM 2.LANDSAT TM 30 x 30 公尺 185 x 185 KM 3.SPOT XS SPOT PAN 20 x 20 公尺 10 x 10 公尺 60 x 60 60 x 60 KM KM 4.SAR 12.5 x 12.5 公尺 100 x 100 KM 1.2 影像價格： 1.LANDSAT MSS ： 1000 元 2.LANDSAT TM 5000 元 3.SPOT XS 視大小及處理等級 SPOT PAN 視大小及處理等級 1.3 影像購買處所： 國立中央大學太空及遙測研究中心代理 (03) 4227151轉7601 1.4 影像處理等級： 等級 輻 射校 正 幾何校正 幾何誤差 影像方向 解析力 Ⅰ YES NONE NONE * X:20M P:10M Ⅳ YES BULK 1500 M * X:20M P:10M Ⅴ YES 系統 1500 M 衛星軌道 X:20M P:10M Ⅵ YES 控制點 12 M 衛星軌道 X:20M P:10M Ⅶ YES 控制點 + DTM 12 M 衛星軌道 X:20M P:10M Ⅷ YES 系統 1500 M 地圖投影 X:12.5M P:6.25M Ⅸ YES 控制點 12 M 地圖投影 X:12.5M P:6.25M Ⅹ YES 控制點 + DTM 12 M 地圖投影 X:12.5M P:6.25M 1.5 數值地形： 現況：本校遙測中心 依五千分之一基本圖編號供4748幅 (其中含 121幅密圖) 價格：XYZ 值之原始影像每幅40元 影像檔 每幅20元 1.6 數值地形購買處所： 國立中央大學太空及遙測研究中心代理 (03) 4227151轉7601 附錄二 利用掃瞄地形圖影像建立數值地形模型 利用掃瞄地形圖影像建立數值地形模型的核心問題在顏色值的分配設計及影像－地理座標與高程轉換的程式設計兩項。 顏色分配 本研究顏色值分配的規則如下： B值：IB=mod(IZ,256).......................………............................(1) G值：IG=255-IZ/256........…………............………..................(2) R值：IR=255 for Z≧0 , IR=0 for Z0......…….......................(3) 其中，Z為高程，IZ為Z的整數值，IB、IG、IR分別為B、G、R波段的顏色值。假如高程值有小數位時，則 IR=iabs(IR-ifix((Z-IZ)*100))..........……….....……….............(4) 上述方程式均為FORTRAN程式的語法。 程式設計 程式主要做兩項工作：1)影像座標及地圖座標轉換，2)顏色值及高程值轉換。後者實際上即為式(1)~(4)的反轉換。前者則可以仿射轉換(Affine transformation)為之，轉換式如下： E=aX+bY+c.................................................(5) N=dX+eY+f.................................................(6) 其中，E、N為二度分帶座標，X、Y為像元座標；a, b, c, d, e, f為轉換參數，可由預設之控制點事先求得。 數化過程 數化地形圖及資料處理的過程如下： 以灰階的方式選用適當的解析度掃瞄地形圖。如為一張五千分之一的基本圖，適當的解析度約為每吋150點(150dpi)。 進入適當影像處理軟體，對於掃瞄入之地形圖影像做一色調加強處理(Color enhancement)，使得等高線在影像中顯得較清楚，以利點圖，並將此影像轉元彩色影像檔。 在等高線上點入先前所設定之顏色值(R, G, B value)，並且在點與點之間保持適當之間隔，例如在等高線彎曲的地方，則點與點之間隔需要密一點。 在影像上的四個角落點上顏色值為255,0,0的像元作為控制點。 儲存完成之影像且離開影像處理軟體。 編輯一個控制點之資料檔。此資料檔儲存控制點之座標值(二度TM座標)。 執行IMGTIZER程式，將顏色值轉為高程值，並按像元座標賦予地理座