Chapter5

Sensor-depentdent Preprocessing/感測器相關前處理功能

5-1. Calibration: Optical Calibration

(1) 用Mapla.bat介面點選「Optical calibration」，目的做光學率定，從感測器DN值推算地表真正輻射值。 1

(2) Input image (.tif) 2

(3) Output image 3

(4) 開啟dim檔檢視影像參數(建議以UltraEdit開啟)，包含衛星、太陽角度參數，以及拍攝時間參數可計算地球與太陽實際距離。 4 5

(5) 參數設定 6

(6) Gains Biases 7

開啟dim檔有Gains Biases參數，注意需要換成倒數，並將參數匯成txt檔。 8 13

(7) Solar illuminations 9

Solar illuminations的來源可到下列網站找尋，並將參數匯成txt檔。 https://books.google.com.tw/books?id=Sne9CgAAQBAJ&pg=PA238&lpg=PA238&dq=formosat-2+irradiance&source=bl&ots=ZKG7Iiwb6z&sig=M9W-AKUg2QL_TvGMEDBpeagh7Jk&hl=zh-TW&sa=X&ved=0ahUKEwif8ZOei5_OAhWItI8KHXmsCoQ4ChDoAQgbMAA#v=onepage&q=formosat-2%20irradiance&f=false 16

15

(8) 再按Execute，便可完成影像率定。 11

(9) 用Monteverdi.bat開啟calibration image 17

5-2. Geometry: Orthorectification and map projections

1. Data Source：

(1) Formosat-2原始影像 太空中心衛星資料庫快速查詢與即時服務系統(https://esoknspodb.nchc.org.tw) (2) Formosat-2參考正射影像 太空中心衛星資料庫快速查詢與即時服務系統(https://esoknspodb.nchc.org.tw) 台灣地區的拼接正射影像，建議進入太空衛星資料庫
(http://esoknspodb.nchc.org.tw/nspo_esok/index.php) 如果沒有帳號密碼的話需要先申請 a.進入"衛星資料查詢"→僅將福衛二號選項勾選→按放大鏡圖示進入查詢 b.右側選單的產品類型只將融合(i)勾選→按下最底下的篩選 c.出現2012-2007年間的2m解析度融合正射產品 (3) DEM: SRTM 下載

方法1.(建議使用此方法)
a. 用Mapla.bat介面點選”Download or list SRTM tiles to a set of images”，可直接下載影像所在區域的DEM。

3

b. Input image (.tif) 4

c. Output directory，再按Execute便可完成DEM下載。 5

方法2. https://dds.cr.usgs.gov/srtm/version2_1/SRTM3/ 選擇區域，再選擇經緯度下載影像所在區域的DEM。

1

2

(4) 控制點來源:國土測繪中心的通用版電子地圖(用IE流覽器開啟http://emap.nlsc.gov.tw/gis103/)

2.影像正射化

(1) 用Mapla.bat介面點選”Read image information”，目的是輸入原始影像(tif檔)，輸出geom檔。 6

(2) 輸入完後出現以下畫視窗。 7

(3) Input image file Chap5-1輻射改正後影像(tif檔)。 8

(4) Output geom檔，再按Execute執行。 9

(5) 用Mapla.bat介面點選”Homologous Points Extraction”，為改正影像的geom檔，利用SIFT或SURF方法找匹配點。

10

(6) Input image file 1。 11

(7) Input image file 2 (reference image)。 12

(8) 設定參數如下圖。 13

(9) Input DEM。 14

(10) Input Geoid，使用geoid model :EGM96 download。 (http://hg.orfeo-toolbox.org/OTBData/raw-file/404aa6e4b3e0/Input/DEM/egm96.grd) 15

(11) Output file，建立一檔名。 16

(12) 勾選紅色框欄處，再按Execute，產生經緯度坐標系統。 17

(13) Output.txt檔。 19

(14) 用Mapla.bat介面點選”Refine Sensor Model”，目地用前步驟產生的匹配點的經緯度坐標來更新geom檔。

20

(15) 輸入完後出現以下畫視窗 21

(16) Input geom file (.geom)。 22

(17) Output geom file (.geom)。 23

(18) Iutput file containing tie points (output_ForRefine)，輸入匹配點的經緯度坐標。 24

(19) 選擇投影坐標系統WGS84。 25

(20) DEM directory、Geoid File步驟同(9)(10)，再按Execute。 27

(21)用Mapla.bat介面點選” Ortho-rectification”，目的做正射改正。 27

(22) Input Image (.tif)。 28

(23) Output Image(輸入檔名.tif) 。 29

(24) 投影坐標系統EPSG code，EPSG code:3826。 各影像有所不同，此影像以台灣為主。 30

(25) DEM directory、Geoid File步驟同(9)(10)，再按Execute。 31

(26) 成功產生全色態正射影像(.tif)，再以Monteverdi開啟影像 + Ref_Ortho。 32 33 34

results matching ""

    No results matching ""