基于和軟件進(jìn)行ads80數(shù)據(jù)生產(chǎn)

1、基于Xpro和ERDAS軟件進(jìn)行ADS80數(shù)據(jù)生產(chǎn)流程北京天圖科技目錄第 1 章處理流程1第 2 章數(shù)據(jù)準(zhǔn)備與處理2數(shù)據(jù)路徑修改22.1輔助數(shù)據(jù)整理22.2第 3 章數(shù)據(jù)預(yù)處理4第 4 章新建工程7第 5 章空三11自動(dòng)生產(chǎn)同名點(diǎn)115.1量測(cè)點(diǎn)125.2執(zhí)行空三135.3第 6 章輸出 L1 級(jí)數(shù)據(jù)18第 7 章生成DEM22自動(dòng)匹配點(diǎn)云格式高程數(shù)據(jù)227.1將點(diǎn)云轉(zhuǎn)化為柵格 DEM247.2第 8 章輸出 L2 級(jí)數(shù)據(jù)26第 9 章影像降位、鑲嵌與標(biāo)準(zhǔn)分幅29影像降位299.1影像鑲嵌與標(biāo)準(zhǔn)分幅299.2I第1章 處理流程整個(gè)處理按如程進(jìn)行：ADS 80 RAW數(shù)據(jù)數(shù)據(jù)預(yù)處理自動(dòng)點(diǎn)匹配與點(diǎn)

2、量測(cè)工程創(chuàng)建生產(chǎn)L1級(jí)數(shù)據(jù)與測(cè)圖軟件對(duì)接進(jìn)行數(shù)字測(cè)圖點(diǎn)云數(shù)據(jù)拼接轉(zhuǎn)換匹配點(diǎn)云格式DSM生成DEM降位、調(diào)色、鑲嵌、分幅生產(chǎn)標(biāo)準(zhǔn)分幅DOM生產(chǎn)L2級(jí)數(shù)據(jù)其中點(diǎn)云數(shù)據(jù)拼接轉(zhuǎn)換、拼接生成 DEM，L2 級(jí)數(shù)據(jù)降位、調(diào)色、鑲嵌、分幅使用 ERDASIMAGINE，其余部分使用Xpro。L1 級(jí)成果大部分國內(nèi)及國外測(cè)圖軟件可以直接對(duì)接使用。1空三計(jì)算第2章 數(shù)據(jù)準(zhǔn)備與處理2.1 數(shù)據(jù)路徑修改ADS80 數(shù)據(jù)包含很多內(nèi)容，處理前應(yīng)把數(shù)據(jù)整理好，放到指定的目錄。一般情況下由于拷貝移動(dòng)等，使得 ADS80 數(shù)據(jù)處理時(shí)和數(shù)據(jù)時(shí)的文件路徑不同，需要修改 sup，info 中的路徑。用文本編輯器打開所有相關(guān)文件，批

3、量修改 IMAGE_ FILE NAME 的存放路徑，使其指向正確的數(shù)據(jù)路徑。2.2 輔助數(shù)據(jù)整理1. 相機(jī)文件： 將相機(jī)文件放到一個(gè)固定路徑后， 需在C:Program FilesLeicaGeosystemsXProXpro.ini 中將該路徑定義好，本例將相機(jī)文件夾放在 C:ADScalibrations目錄下，則對(duì) Xpro.ini 的修改如下圖所示。22.資料：點(diǎn)坐標(biāo)文件，點(diǎn)之記等，點(diǎn)坐標(biāo)文件以“IDXYZ”形式的文本文件即可。3.航線信息（非必須）：航線軌跡圖，飛行方向及架次等。3第3章 數(shù)據(jù)預(yù)處理1、點(diǎn)擊開始菜單Leica GeosystemsXProLeica XPro，啟動(dòng)

4、XPro 程序打開 XPro 軟件2、點(diǎn)擊數(shù)據(jù)預(yù)處理圖標(biāo)，打開數(shù)據(jù)預(yù)處理界面3、點(diǎn)擊 Browse圖標(biāo)，選擇保存 ADS80 數(shù)據(jù)目錄（包含 gps-imu、raw 的那一級(jí)目錄），點(diǎn)擊 Open 按鈕開始驗(yàn)證文件，注意查看提示，如果出現(xiàn) error 要根據(jù)報(bào)錯(cuò)信息檢查數(shù)據(jù)準(zhǔn)備階段是否有問題。44、檢查完畢后，在 Data Preparation框上勾選中 Georeferencing with 選項(xiàng)框，并勾選上 Pyramids，Statistics框5、在 Data Preparation框上，在Take下面勾選中需要進(jìn)行數(shù)據(jù)預(yù)處理的那些ADS80 條帶，在 All Bands 右側(cè)的圖標(biāo)

5、點(diǎn)擊，選擇要處理的波段。5如果需要設(shè)定footprint 的路徑，可以在shapefile 處選擇。6、設(shè)置完成后，在Data PreparationYes框上，點(diǎn)擊 Process 按鈕，在彈出的框點(diǎn)擊在彈出框選擇Run Locally，點(diǎn)擊 Run 執(zhí)行，對(duì)ADS80 數(shù)據(jù)進(jìn)行預(yù)處理。6第4章 新建工程1)在 XPro 界面上，點(diǎn)擊新建工程圖標(biāo)，打開新建工程框2)點(diǎn)擊File/ New Block，打開 New Block M框3)在 New Block M入工程名稱框上點(diǎn)擊 Browse 按鈕選擇新建工程的路徑，在 Name 后面輸4)點(diǎn)擊Next 按鈕，設(shè)置需要新建工程的投影信息75)

6、在 Coordinate System 后面選擇 Projection System6)在 Choose Projection System框上選擇與所做的 ADS80 數(shù)據(jù)對(duì)應(yīng)的投影，點(diǎn)擊OK，關(guān)閉 Choose Projection System框7)在 Block Preparation框上點(diǎn)擊Block/ Add Take 按鈕，加載 stp 文件88)9)選擇需要處理的.stp 文件，單擊 open 加載選中的數(shù)據(jù)如果有點(diǎn)，點(diǎn)擊 Block > Edit Ground Control.點(diǎn)擊左上角 Import Control Points 導(dǎo)入點(diǎn)。根據(jù)點(diǎn)實(shí)際投影坐標(biāo)設(shè)定相關(guān)信息

7、9導(dǎo)入之后，可根據(jù)每個(gè)點(diǎn)情況修改 Type 列，修改完畢點(diǎn)擊Save and Close10) 點(diǎn)擊File/Save Block 按鈕，保存新建的工程10第5章 空三5.1自動(dòng)生產(chǎn)同名點(diǎn)1)，打開 Triangulation點(diǎn)擊空三圖標(biāo)框。2)在 Triangulation框上點(diǎn)擊File/Open Block 按鈕，打開 Block Selection框3)在 Block Selection框上點(diǎn)擊 Browse選擇需要加載的工程文件114) 在 Triangulation框上點(diǎn)擊Run APM 圖標(biāo)，打開 Run APM框5) 在 Run APM框上點(diǎn)擊 OK 按鈕，關(guān)閉Run APM框

8、，自動(dòng)生成同名點(diǎn)。5.2量測(cè)點(diǎn)1) 在空三界面中可以看到點(diǎn)以三角形顯示， 在需要量測(cè)的點(diǎn)上點(diǎn)擊右鍵選擇Measurement，會(huì)彈出與該點(diǎn)相關(guān)的所有影像122)在 Measurement 窗口處根據(jù)需要使用相關(guān)工具調(diào)節(jié)縮放和顏色，使用刺點(diǎn)圖標(biāo)添加控制點(diǎn)3)點(diǎn)擊主界面的Window > Point List 可以打開Point List for GCP，查看點(diǎn)列表，可以看到哪些點(diǎn)已經(jīng)刺過，也可以通過這個(gè)列表選擇某個(gè)點(diǎn)，點(diǎn)擊下圖的圖標(biāo)直接跳轉(zhuǎn)。5.3執(zhí)行空三1) 在 Triangulation框上，點(diǎn)擊 Create scenario 按鈕，打開 New Scenario Name 對(duì)13話

9、框2) 在 New Scenario Name框上輸入新的名字，點(diǎn)擊 OK 按鈕，關(guān)閉 New ScenarioName框。3) 空三參數(shù)設(shè)定，一般調(diào)節(jié)基本參數(shù)就能完成計(jì)算-基本參數(shù): sigma0 可使用默認(rèn) 3.5m；如果有點(diǎn)參與運(yùn)算，GPS 參數(shù)可設(shè)定為Datum，若是多架次則設(shè)為 Datum per Session，在極少數(shù)情況下如果有明顯的系統(tǒng)誤差，可嘗試設(shè)為 Drift； IMU，若是多架次則設(shè)為 Session，在極少數(shù)情況下如果有明顯的系統(tǒng)誤差，設(shè)為 Take，當(dāng)有至少兩個(gè)從相反方向飛行的重疊條帶時(shí)，可設(shè)為 Sensor。4)以上設(shè)置完畢后，在 Triangulation框上，

10、點(diǎn)擊Run Adjustment 按鈕，運(yùn)行空三5)空三完畢可以看到如下圖的整體誤差信息，一般紅域?yàn)樾枰M(jìn)行修正的區(qū)域，可以通過調(diào)整連接點(diǎn)或是匹配點(diǎn)的方式進(jìn)行修改。146)在 Triangulation框上，選中粗差較大的點(diǎn)，單擊右鍵選擇 Measurement，在該界面的下方會(huì)顯示該點(diǎn)所在影像的多個(gè)視窗，調(diào)整位置偏差較大的影像的點(diǎn)，也可以選中粗差較大的點(diǎn)，單擊右鍵選擇Remove Point 按鈕，刪除粗差較大的點(diǎn)。7)也可通過點(diǎn)列表瀏覽，查看紅色的點(diǎn)，進(jìn)行重新量測(cè)或調(diào)整。調(diào)節(jié)完可以再次執(zhí)行空三，查看誤差。158)若點(diǎn)調(diào)節(jié)完后，誤差仍難以滿足要求，可以嘗試使用 Triangulation-A

11、nalysis-APM loop功能，在基本參數(shù)設(shè)定界面，勾選 Densify Points Threshold，設(shè)定循環(huán)次數(shù)，并根據(jù)需要設(shè)定高級(jí)參數(shù)中的 Densify Points Threshold 閾值。執(zhí)行空三即可。9)若空三誤差已經(jīng)滿足要求，點(diǎn)擊 Write to Block 按鈕，將空三后的結(jié)果寫入工程中。10) 若需要導(dǎo)出點(diǎn)，點(diǎn)擊Scenario 下的 Export Ground Points16兩個(gè)文件會(huì)被生成：adjustedGroundControlPoints.txt、adjustedGroundTiePoints.txt.17第6章 輸出 L1 級(jí)數(shù)據(jù)1)點(diǎn)擊Prod

12、uct Generation 圖標(biāo)，打開 Product Generation框，點(diǎn)擊File/OpenBlock 按鈕，選擇空三后的工程文件2)在 Product Generation框上點(diǎn)擊 Generate 圖標(biāo)，打開 Generate Product框3)在 Generate Product框上選中 L1_RGB，點(diǎn)擊 Next 按鈕184)在 Generate Product框上選中Atmospheric + BRDF 選項(xiàng)，點(diǎn)擊Next 按鈕5)6)在 Generate Product在 Generate Product框上選擇 Standard 選項(xiàng)，點(diǎn)擊Next 按鈕框上點(diǎn)擊O

13、utput Root 后面的按鈕，設(shè)置輸出文件的路徑，在 Perfix 后面設(shè)置輸出影像的前綴。197)在 Generate Product框上，點(diǎn)擊 Finish 按鈕，關(guān)閉 Generate Product框8)在框上點(diǎn)擊 Execute jobs 圖標(biāo)，打開 Process Locally or Remotely框9)在 Process Locally or Remotely框上點(diǎn)擊Run Locally 按鈕，運(yùn)行程序，此成果輸出后即可得到ADS80 的 1 級(jí)RGB 數(shù)據(jù)。10) 如果需要 1 級(jí) PAN 數(shù)據(jù)，則需在Edit 點(diǎn)擊 Products2011) 在 Product T

14、ype Editor 中新建 L1_PAN，參考 L1_RGB 設(shè)定相關(guān)參數(shù)，保存。在前面數(shù)據(jù)生產(chǎn)中選擇新建的 L1_PAN 數(shù)據(jù)類型即可獲得 1 級(jí) PAN 數(shù)據(jù)。21第7章 生成 DEM7.1自動(dòng)匹配點(diǎn)云格式高程數(shù)據(jù)1)點(diǎn)擊Product Generation 圖標(biāo)，打開 Product Generation框，點(diǎn)擊File/OpenBlock 按鈕，選擇空三后的工程文件2)在 Product Generation框上點(diǎn)擊 Generate 圖標(biāo)，打開 Generate Product框223)在 Generate Product框上選中 DSM_RGB，點(diǎn)擊 Next 按鈕4)在 Gen

15、erate Product框上選中Atmospheric + BRDF 選項(xiàng)，點(diǎn)擊Next 按鈕5)6)在 Generate Product在 Generate Product框上選擇Standard 選項(xiàng)，點(diǎn)擊Next 按鈕框上點(diǎn)擊Output Root 后面的按鈕，設(shè)置輸出文件的路徑，在 Perfix 后面設(shè)置輸出影像的前綴。237)在 Generate Product框上，點(diǎn)擊 Finish 按鈕，關(guān)閉 Generate Product框8)在框上點(diǎn)擊 Execute jobs 圖標(biāo)，打開 Process Locally or Remotely框9)在 Process Locally o

16、r Remotely出后即可得到點(diǎn)云格式的高程數(shù)據(jù)?？蛏宵c(diǎn)擊Run Locally 按鈕，運(yùn)行程序，此成果輸7.2將點(diǎn)云轉(zhuǎn)化為柵格 DEM1) DEM 提取完成后，需要使用 ERDAS IMAGINE 將點(diǎn)云轉(zhuǎn)換為柵格格式 DEM，才能繼續(xù)進(jìn)行。242)點(diǎn)擊 ERDAS IMAGINE 主界面 Terrain下的 TPT(Terrain Prep Tool)進(jìn)行合并，點(diǎn)擊 TPT 中圖標(biāo)加載匹配好的 las 文件3)點(diǎn)擊下進(jìn)入 Surface DTM 界面，設(shè)置 Rasterization 選項(xiàng)下柵格 DEM 的文件名、路徑和分辨率，點(diǎn)擊 OK 執(zhí)行。4)輸出的結(jié)果就是柵格格式的 DEM。25

17、第8章 輸出 L2 級(jí)數(shù)據(jù)1)點(diǎn)擊Product Generation 圖標(biāo)，打開 Product Generation框，點(diǎn)擊File/OpenBlock 按鈕，選擇空三后的工程文件2)在 Product Generation框上點(diǎn)擊 Generate 圖標(biāo)，打開 Generate Product框3)在 Generate Product框上選中 L2_RGB，點(diǎn)擊 Next 按鈕264)在 Generate Product框上選中Atmospheric + BRDF 選項(xiàng)，點(diǎn)擊Next 按鈕5)6)在 Generate Product在 Generate Product框上選擇Standa

18、rd 選項(xiàng)，點(diǎn)擊Next 按鈕框上設(shè)定 GSD 地面分辨率；點(diǎn)擊 Elevation Source 中 Selecta DEM，選擇上一步轉(zhuǎn)換好的柵格 DEM；點(diǎn)擊 Output Root 后面的按鈕，設(shè)置輸出文件的路徑，在 Perfix 后面設(shè)置輸出影像的前綴。277)在 Generate Product框上，點(diǎn)擊 Finish 按鈕，關(guān)閉 Generate Product框8)在框上點(diǎn)擊 Execute jobs 圖標(biāo)，打開 Process Locally or Remotely框9)在 Process Locally or Remotely框上點(diǎn)擊Run Locally 按鈕，運(yùn)行程序，

19、此成果輸出后即可得到ADS80 的 2 級(jí)正射的 RGB 數(shù)據(jù)。28第9章 影像降位、鑲嵌與標(biāo)準(zhǔn)分幅9.1影像降位1)由于默認(rèn)輸出的為 16 位的數(shù)據(jù)，為了生產(chǎn)需要將數(shù)據(jù)進(jìn)行降位處理2)點(diǎn)擊 ERDAS IMAGINE 主界面 Raster 選項(xiàng)卡下 Radiometric中的 Rescale 進(jìn)行降位，設(shè)置輸入輸出數(shù)據(jù)，算法為 Minimum-要批處理，點(diǎn)擊Batch 進(jìn)行。um，輸出統(tǒng)計(jì)忽略 0，點(diǎn)擊 OK 即可，如果需9.2影像鑲嵌與標(biāo)準(zhǔn)分幅1)點(diǎn)擊 ERDAS IMAGINE 主界面 Toolbox 下 MosaicPro，選擇降好位的數(shù)據(jù)，點(diǎn)擊 Use Entire Image 按鈕（如果影像背景區(qū)域較大可選擇 Compute Act