遙感圖像處理實驗課程

2023/2/131實驗一認識ERDASIMAGING軟件ERDASIMAGING軟件概述ERDASIMAGING是美國ERDAS公司開發(fā)的專業(yè)遙感圖像處理與地理信息系統(tǒng)軟件，以模塊化的方式提供給用戶，可使用戶根據(jù)自己的應用要求、資金情況合理地選擇不同的功能模塊以及不同的組合，對系統(tǒng)進行精簡，充分利用軟硬件資源，并最大限度的滿足用戶的專業(yè)應用要求。ERDASIMAGING面對不同需求的用戶，系統(tǒng)的擴展功能采樣開放的體系結(jié)構(gòu)，以IMAGINGEssentials、IMAGINGAdvantage、IMAGINGProfessional的形式為用戶提供了低、中、高三檔產(chǎn)品構(gòu)架，并有豐富的功能擴展供用戶選擇，使產(chǎn)品模塊的組合具有極大的靈活性。第1頁/共161頁第一頁，共162頁。2023/2/132IMAGINGEssentials級這是一個包括制圖和可視化核心功能的價格較低的圖像工具軟件，無論是獨立地從事工作還是處在企業(yè)協(xié)同計算的環(huán)境下，都可以借助IMAGINGEssentials完成二維/三維顯示、數(shù)據(jù)輸入、排序與管理、地圖配準、專題制圖以及簡單的分析。可以集成使用多種數(shù)據(jù)類型，并在保持相同的易于使用和易于剪裁的界面下，升級到其他級別的ERDASIMAGING軟件?？蓴U充的模塊包括：Vector模塊：直接采用GIS業(yè)界領袖ESRI的ArcInfo數(shù)據(jù)結(jié)構(gòu)Coverage，可以建立、顯示、編輯和查詢Covages，完成拓撲關系的建立和修改，實現(xiàn)矢量圖形和柵格圖像的雙向轉(zhuǎn)換等。第2頁/共161頁第二頁，共162頁。2023/2/133VirtualGIS模塊：功能強大的三維可視化分析工具，可以完成實時的三維飛行模擬、建立虛擬世界、進行空間視域分析、矢量與柵格的三維疊加、空間GIS分析等。Developer’Toolkit模塊：ERDASIMAGING的C語言開發(fā)工具包，包含了數(shù)百個函數(shù)，是ERDASIMAGING客戶化的基礎。IMAGINGAdvantage級IMAGINGAdvantage級建立在IMAGINGEssentials級基礎之上，增加了更豐富的柵格圖像GIS分析和單張航片正設校正等軟件。第3頁/共161頁第三頁，共162頁。2023/2/134IMAGINGAdvantage為用戶提供了靈活可靠的用于柵格分析、正設校正、地形編輯及圖像拼接工具。IMAGINGAdvantage級是一個完整的圖像地理信息系統(tǒng)?？蓴U充的模塊包括：第4頁/共161頁第四頁，共162頁。2023/2/135OrthoMAX模塊：全功能、高性能的數(shù)字航測軟件，依據(jù)立體像對進行正設校正、自動DEM提取、立體地形顯示及浮動光標方式的DEM交互編輯等。Radar模塊：完成雷達圖像的基本處理，包括亮度調(diào)整、斑點噪聲消除、紋理分析、邊緣提取等功能。第5頁/共161頁第五頁，共162頁。2023/2/136OrthoBase模塊：區(qū)域數(shù)字攝影測量模塊，用于航空影像的空中測量和正設校正。OrthoRadar模塊：可對Radarsat、ERS雷達圖像進行地理編碼、正設校正等處理。StereoSARDEM模塊：采用類似于立體測量的方法，從雷達圖像數(shù)據(jù)中提取DEM。IFSARDEM模塊：采用干涉方法，以像對為基礎從雷達圖像數(shù)據(jù)中提取DEM。ATCOR模塊：用于大氣因子校正和霧曦消除。第6頁/共161頁第六頁，共162頁。2023/2/137ERDASIMAGING圖標面板ERDASIMAGING圖標面板包括菜單條和工具條兩部分，提供了啟動ERDASIMAGING軟件模塊的全部菜單。命令菜單及其功能菜單菜單功能Session(綜合菜單)完成系統(tǒng)設置、面板布局、日志管理，啟動命令工具、批處理過程、實用功能、聯(lián)機幫助等Main主菜單啟動ERDAS圖標面板中包括的所有功能模塊Tools工具菜單完成文本編輯、矢量及柵格數(shù)據(jù)屬性編輯、圖形圖像文件坐標轉(zhuǎn)換、注記記字體管理、三維動畫制作Utilities實用菜單完成多種柵格數(shù)據(jù)格式的設置與轉(zhuǎn)換、圖像的比較Help幫助菜單啟動關于圖標面板的聯(lián)機幫助、ERDASIMAGING聯(lián)機文檔查看、動態(tài)連接庫瀏覽等表1.1ERDASIMAGING圖標面板菜單條第7頁/共161頁第七頁，共162頁。2023/2/138命令功能Preference面向單個用戶或全體用戶，設置多數(shù)功能模塊的系統(tǒng)默認值Configuration為ERDASIMAGING配置各種外圍設備，如打印機、磁帶機SessionLog查看ERDASIMAGING提示、命令及運行過程中的實時記錄ActiveProcessList查看與取消ERDASIMAGING系統(tǒng)當前正在運行的處理操作Commands啟動命令工具，進入命令菜單狀態(tài)，通過命令執(zhí)行處理操作EnterLogMessage向系統(tǒng)綜合日志輸入文本信息StartBatchCommands啟動或退出批處理工具，打開批處理向?qū)?，記錄批處理命令OpenBatchFile打開批處理命令文件ViewBatchQueue打開批處理進程對話框，查看、編輯、刪除批處理隊列FlipIcon確定圖標面板的水平或垂直顯示狀態(tài)TileViewers平鋪排列兩個以上已經(jīng)打開的窗口CloseAllViewers關閉當前打開的所有窗口表1.2Session菜單命令及其功能第8頁/共161頁第八頁，共162頁。2023/2/139命令功能Main進入主菜單，啟動圖標面板中包括的所有模塊Tools進入工具菜單，顯示和編輯文本及圖像文件Utilities進入實用菜單，執(zhí)行ERDAS的常用功能Help打開ERDASIMAGING聯(lián)機幫助文檔Properties打開ERDAS系統(tǒng)特性對話框，查看和配置序列號與模塊ExitIMAGING退出ERDASIMAGING軟件環(huán)境第9頁/共161頁第九頁，共162頁。2023/2/1310命令功能IMAGINGCredits查閱ERDAS信用卡StartIMAGINGViewer啟動ERDASIMAGING窗口Import/Export啟動ERDASIMAGING數(shù)據(jù)輸入/輸出模塊DataPreparation啟動ERDASIMAGING數(shù)據(jù)預處理模塊MapComposer啟動ERDASIMAGING專題制圖模塊ImageInterpreter啟動ERDASIMAGING圖像解譯模塊ImageCatalog啟動ERDASIMAGING圖像庫管理模塊ImageClassification啟動ERDASIMAGING圖像分類模塊SpatialModeler啟動ERDASIMAGING空間建模工具Vector啟動ERDASIMAGING矢量功能模塊Radar啟動ERDASIMAGING雷達圖像處理模塊VirtualGIS啟動ERDASIMAGING虛擬GIS模塊OrthoBase啟動ERDASIMAGING正射影像校正模塊SubpixelClassifier啟動ERDASIMAGING子像元分類模塊表1.3Main命令及其功能第10頁/共161頁第十頁，共162頁。2023/2/1311命令功能EditTextFiles編輯ASCII碼文本文件EditRasterAttributes編輯柵格文件屬性數(shù)據(jù)ViewBinaryData查看二進制文件的內(nèi)容ViewHFAFileStructure查看ERDASIMAGING層次文件結(jié)構(gòu)AnnotationInformation查看注記文件信息，包括元素數(shù)量與投影參數(shù)ImageInformation獲取ERDASIMAGING柵格圖像文件的所有信息VectorInformation獲取ERDASIMAGING矢量圖形文件的所有信息ImageCommandsTool打開圖像命令對話框，進入ERDAS命令操作環(huán)境CoordinateCalculator將坐標系統(tǒng)從一種橢圓體或參數(shù)轉(zhuǎn)變?yōu)榱硗庖环NCreate/DisplayMovieSequences產(chǎn)生和顯示一系列圖像畫面形成的動畫Create/DisplayViewerSequences產(chǎn)生和顯示一系列窗口畫面組成的動畫ImageDrape以DEM為基礎的三維圖像顯示與操作表1.4Tools菜單命令及其功能第11頁/共161頁第十一頁，共162頁。2023/2/1312命令功能JPEGCompressImage應用JPEG壓縮技術對柵格圖像進行壓縮，以便于保存DecompressJPEGImage將應用JPEG壓縮技術所產(chǎn)生的柵格圖像進行解壓縮ConvertPixelstoASCII將柵格圖像文件數(shù)據(jù)轉(zhuǎn)換為ASCII碼文件ConvertASCIItoPixels以ASCII嗎文件為基礎產(chǎn)生柵格圖像文件ConvertImagetoAnnotation將柵格圖像文件轉(zhuǎn)換成IMAGING的多邊形注記數(shù)據(jù)ConvertAnnotationtoRaster將IMAGING的多邊形注記數(shù)據(jù)轉(zhuǎn)換成柵格圖像文件CreateFontTables以特定的字體生成柵格符號圖FonttoSymbol將特定的字體轉(zhuǎn)換為地圖符號CompareImages打開圖像比較對話框，比較兩幅圖像之間的某種屬性ReconfigureRasterFormats重新配置系統(tǒng)中柵格圖像數(shù)據(jù)格式ReconfigureVectorFormats重新配置系統(tǒng)中矢量圖形數(shù)據(jù)格式ReconfigureResampleMethods重新設置系統(tǒng)中圖像重采樣方法ReconfigureGeometricMethods重新設置系統(tǒng)中圖像幾何校正方法表1.5Utilities菜單命令及其功能第12頁/共161頁第十二頁，共162頁。2023/2/1313命令功能HelpforICONPanel顯示ERDASIMAGING圖標面板的聯(lián)機幫助IMAGINGonLineDocumentation進入聯(lián)機幫助目錄，查看IMAGING聯(lián)機文檔IMAGINGVersion查看正在進行的ERDASIMAGING版本IMAGINGDLLInformation查看IMAGING動態(tài)連接庫的類型與常數(shù)信息表1.6Help菜單命令及其功能第13頁/共161頁第十三頁，共162頁。2023/2/1314工具圖標及其功能與IMAGINGProfessional級對應的圖標面板工具條中的圖標有10個，除了IMAGINGProfessional級功能及其應用外，還將涉及4個重要的擴展模塊：Subpixel模塊、OrthoBase模塊、Vector模塊和VirtualGIS模塊，共有14個圖標。第14頁/共161頁第十四頁，共162頁。2023/2/1315Viewer視圖，幾乎所有的操作都在視圖中進行解譯模塊：進行圖像各種增強處理輸入輸出模塊：各種數(shù)據(jù)格式的轉(zhuǎn)換數(shù)據(jù)預處理模塊：圖像分幅、幾何校正第15頁/共161頁第十五頁，共162頁。2023/2/1316分類操作：監(jiān)督分類、非監(jiān)督分類虛擬GIS模塊：三維可視化矢量功能模塊第16頁/共161頁第十六頁，共162頁。2023/2/1317ERDASIMAGING功能體系第17頁/共161頁第十七頁，共162頁。2023/2/1318實習二、視窗操作一、窗口功能簡介二、File菜單操作三、Utility菜單操作四、View菜單操作五、AOI菜單操作六、Raster菜單操作七、Vector菜單操作第18頁/共161頁第十八頁，共162頁。2023/2/1319一、窗口功能簡介矢量視窗柵格視窗第19頁/共161頁第十九頁，共162頁。2023/2/1320二、File菜單操作圖像顯示操作第一步：啟動程序（StartProgram）第二步：確定文件（DetermineFile）第三步：設置參數(shù)（RasterOptions）第四步：打開圖像（OpenRasterLayer）第20頁/共161頁第二十頁，共162頁。2023/2/1321第21頁/共161頁第二十一頁，共162頁。2023/2/1322第22頁/共161頁第二十二頁，共162頁。2023/2/1323應用不同的波段組合顯示效果第23頁/共161頁第二十三頁，共162頁。2023/2/1324圖形顯示操作第一步：啟動程序（StartProgram）第二步：確定文件（DetermineFile）第三步：設置參數(shù)（VectorOptions）第四步：打開圖形（OpenVectorLayer）第24頁/共161頁第二十四頁，共162頁。2023/2/1325第25頁/共161頁第二十五頁，共162頁。2023/2/1326第26頁/共161頁第二十六頁，共162頁。2023/2/1327應用不同的符號顯示效果第27頁/共161頁第二十七頁，共162頁。2023/2/1328矢量視窗柵格視窗第28頁/共161頁第二十八頁，共162頁。2023/2/1329Utility第29頁/共161頁第二十九頁，共162頁。2023/2/1330命令功能InquireCursor啟動光標查詢功能InquireBox啟動查詢光標框InquireColor設置查詢光標顏色InquireShape設置查詢光標形狀InquireHome查詢光標復位于視窗中心Measure啟動系統(tǒng)量測功能SelectorProperties改變當前查詢光標的特性Blend混合顯示上下兩層圖像Swipe卷簾顯示上下兩層圖像Flicker閃爍顯示上下兩層圖像LayerInfo顯示視窗上層文件信息VirtualGIS啟動虛擬地理信息系統(tǒng)功能三、Utility菜單操作表2.1實用菜單命令與功能第30頁/共161頁第三十頁，共162頁。2023/2/1331在打開InquireCursor窗口的同時，窗口中出現(xiàn)十字查詢光標，窗口與十字光標是同步關聯(lián)的，在窗口中任意移動十字光標，窗口中的信息動態(tài)更新。窗口中顯示出了十字光標當前所在位置像元的縱橫坐標—地圖坐標Map、文件坐標File、地理坐標Lat/Lon、圖紙坐標Paper；3個波段的顏色、灰度值；LUT表值；像元大小直方圖等信息。利用對話框左下方的上下左右4個移動功能鍵或一定十字光標的位置，窗口中的像元信息實時發(fā)生變化。第31頁/共161頁第三十一頁，共162頁。2023/2/1332查詢光標的位置與顏色查詢光標第32頁/共161頁第三十二頁，共162頁。2023/2/1333查詢光標的位置與顏色第33頁/共161頁第三十三頁，共162頁。2023/2/1334查詢框第34頁/共161頁第三十四頁，共162頁。2023/2/1335在MeasurementToolforViewer窗口中執(zhí)行下列操作：單擊MeasurePosition圖標，在窗口中單擊確定任意一點，量測其位置。單擊MeasureLengths圖標，在窗口中連續(xù)單擊確定一條折線，量測其長度。單擊MeasureAreas圖標，在窗口中連續(xù)單擊確定一多邊形，量測周長與面積單擊MeasureRectangular圖標，在窗口中按左鍵拖動確定一矩形，量測周長與面積。單擊MeasureEllipses圖標，在窗口中按左鍵拖動確定一橢圓，量測其周長與面積。量測結(jié)果分別以1、2、3、4、5的順序編號顯示在MeasurementToolforViewer窗口。測量功能第35頁/共161頁第三十五頁，共162頁。2023/2/1336測量點的位置，線的長度，多邊形的面積第36頁/共161頁第三十六頁，共162頁。2023/2/1337數(shù)據(jù)疊加顯示數(shù)據(jù)疊加顯示是針對具有相同地理參考系統(tǒng)(地圖投影和坐標系統(tǒng))的兩個文件進行操作的，所以，在進行數(shù)據(jù)疊加操作前，首先根據(jù)需要在一個窗口中打開兩個文件。疊加數(shù)據(jù)準備：在同一窗口中打開兩個文件lanier.img和insoils.img(或eldodem.img和eldoatm.img)，在第二個文件打開時，在RasterOptions選項卡中，取消ClearDisplay復選框。疊加顯示數(shù)據(jù)：ERDAS系統(tǒng)提供的數(shù)據(jù)疊加顯示工具由3個，分別是混合顯示工具(BlendTool)、卷簾顯示工具(SwipeTool)和閃爍顯示工具(FlickerTool)。第37頁/共161頁第三十七頁，共162頁。2023/2/1338混合顯示（BlendTool）

通過控制上層圖像顯示的透明度大小，使得上下兩層圖像混合顯示第38頁/共161頁第三十八頁，共162頁。2023/2/1339卷簾顯示（SwipeTool）水平卷簾和垂直卷簾顯示應用過渡線將上層數(shù)據(jù)文件分為不透明和透明兩個部分，移動過渡線同時顯示上下兩層數(shù)據(jù)文件，查看其相互關系第39頁/共161頁第三十九頁，共162頁。2023/2/1340閃爍工具（FlickerTool）自動比較上下兩層圖像的屬性差異及其關系，經(jīng)常應用該操作的典型實例是分類專題圖像與原始圖像之間的比較

分類專題圖像和原始圖像第40頁/共161頁第四十頁，共162頁。2023/2/1341圖像直方圖文件信息操作波段1的文件信息文件通常信息第41頁/共161頁第四十一頁，共162頁。2023/2/1342波段5的文件信息文件信息操作第42頁/共161頁第四十二頁，共162頁。2023/2/1343不同波段的直方圖波段直方圖第43頁/共161頁第四十三頁，共162頁。2023/2/1344波段像元值波段像元值第44頁/共161頁第四十四頁，共162頁。2023/2/1345四、View菜單操作文件顯示順序（ArrangeLayer）

第45頁/共161頁第四十五頁，共162頁。2023/2/1346顯示比例操作（DisplayScale）用于調(diào)整文件顯示比例及其與視窗的對比關系→ImagetoWindow：按照視窗大小調(diào)整文件顯示比例→WindowtoImage：按照文件大小調(diào)整視窗顯示比例→Extent：顯示文件整體范圍→ScaleTool：顯示文件整體范圍第46頁/共161頁第四十六頁，共162頁。2023/2/1347顯示變換操作（Rotate/Flip/Stretch）顯示變換操作實質(zhì)上是對圖像進行仿射變換，不過只是顯示而已，并不是對文件數(shù)據(jù)進行操作，變換操作可以同時進行縮放、平移、拉伸、旋轉(zhuǎn)、鏡面等線性變換第47頁/共161頁第四十七頁，共162頁。2023/2/1348五、AOI菜單操作AOI是用戶感興趣區(qū)域（AreaofInterest）的縮寫，AOI菜單包含了AOI工具（AOITools）等14項命令，分別應用于完成與AOI有關的文件操作，確定了一個AOI之后，可以使相關的ERDASIMAGING的處理針對AOI內(nèi)的像元進行。AOI區(qū)域可以保存為一個文件，便于在以后的多種場合調(diào)用，AOI區(qū)域經(jīng)常應用于圖像分類模板文件的定義以及裁剪框的定義。第48頁/共161頁第四十八頁，共162頁。2023/2/1349打開AOI工具面板（AOIToolPalette）AOI工具面板中包含了所有的AOI菜單操作命令。AOI工具面板大致可以分為三個功能區(qū)，前兩排圖標是產(chǎn)生AOI與選擇AOI功能區(qū)、中間兩排是編輯AOI功能區(qū)、而最后兩排則是定義AOI屬性功能區(qū)，掌握AOI工具面板中的命令功能，對于在圖像處理工作中正確使用AOI功能，發(fā)揮AOI的作用是非常有意義的。第49頁/共161頁第四十九頁，共162頁。2023/2/1350

命令

功能Tools打開AOI工具面板Undo取消編輯操作（可多次取消）Cut剪切AOI區(qū)域DeleteRasterMasker刪除柵格掩模Group建立AOI要素組合Ungroup取消AOI要素組合Reshape改變AOI要素形狀InvertPolygon選擇AOI區(qū)域以外的要素ElementPropertiesAOI要素特性StylesAOI顯示特性SeedPropertiesAOI種子特征CopySelectiontoAOI向AOI復制選擇要素Link建立視窗AOI連接TabletInputNewConfigurationCurrentConfiguration數(shù)字化儀輸入建立新的數(shù)字化儀配置調(diào)用已配置的數(shù)字化儀表2.2AOI菜單命令與功能第50頁/共161頁第五十頁，共162頁。2023/2/1351定義AOI顯示特性定義AOI種子特征表2.3AOI種子特征參數(shù)及含義參數(shù)含義Neighborhood：種子增長模式：4NeighborhoodMode4個相鄰像元增長模式8NeighborhoodMode8個相鄰像元增長模式GeographicConstraints：種子增長的地理約束Area（Pixels/Hectares/Acres）面積約束（像元個數(shù)、面積）Distance（Pixels/Meters/Feet）距離約束（像元個數(shù)、距離）SpectralEuclideanDistance光譜歐式距離GrowatInquire以查詢光標為種子增長SetConstraintAOI以AOI區(qū)域為約束條件Options：選擇項定義：IncludeIslandPolygon允許島狀多邊形存在UpdateRegionMean重新計算AOI區(qū)域均值BufferRegionBoundary對AOI區(qū)域進行Buffer第51頁/共161頁第五十一頁，共162頁。2023/2/1352六、柵格菜單操作柵格菜單操作是針對柵格圖像而言的，是視窗菜單操作命令中的重要內(nèi)容，柵格菜單包含了柵格工具（RasterTools）等20項命令，可以對柵格圖像數(shù)據(jù)進行多種變換處理或?qū)傩栽O置。第52頁/共161頁第五十二頁，共162頁。2023/2/1353波段組合第53頁/共161頁第五十三頁，共162頁。2023/2/1354圖像對比度調(diào)整（AdjustImageContrast）局部線性拉伸（PiecewiseLinerStretch）直方圖斷點操作（HistogramBreakpoints）柵格屬性編輯（RasterAttribute）改變分類圖斑顏色（ChangeClassColor）設置圖像的透明度（MakeLayersTransparent）增加柵格屬性表字段（AddAttributeColumn）生成分類統(tǒng)計報告（GenerateStatisticalReport）第54頁/共161頁第五十四頁，共162頁。2023/2/1355圖像剖面工具（ProfileTool）第55頁/共161頁第五十五頁，共162頁。2023/2/1356光譜剖面曲線（SpectralProfile）第56頁/共161頁第五十六頁，共162頁。2023/2/1357空間剖面曲線（SpatialProfile）第57頁/共161頁第五十七頁，共162頁。2023/2/1358三維空間剖面（SurfaceProfile）第58頁/共161頁第五十八頁，共162頁。2023/2/1359七、Vector菜單操作EnableEditing：開始編輯矢量文件Tools：矢量工具欄第59頁/共161頁第五十九頁，共162頁。2023/2/1360顯示矢量屬性特征第60頁/共161頁第六十頁，共162頁。2023/2/1361修改顯示符號第61頁/共161頁第六十一頁，共162頁。2023/2/1362第62頁/共161頁第六十二頁，共162頁。2023/2/1363第63頁/共161頁第六十三頁，共162頁。2023/2/1364第64頁/共161頁第六十四頁，共162頁。2023/2/1365實習三：輸入輸出ERDASIMAGIMG的數(shù)據(jù)輸入輸出功能（Import/Export），允許用戶輸入多種格式的數(shù)據(jù)供ERDAS使用，同時允許將IMAGING的文件轉(zhuǎn)換成多種數(shù)據(jù)格式。目前，IMAGING8.4可以輸入的數(shù)據(jù)格式達70多種，可以輸出的數(shù)據(jù)格式近30種，幾乎包括常用或常見的柵格數(shù)據(jù)和矢量數(shù)據(jù)格式，具體的數(shù)據(jù)格式都羅列在IMAGING輸入輸出對話框中。數(shù)據(jù)輸入輸出的一般操作過程如下：ERDAS圖標面板菜單條：Main→Import/Export→數(shù)據(jù)輸入輸出對話框ERDAS圖標面板菜單條：點擊Import/Export圖標→數(shù)據(jù)輸入輸出對話框第65頁/共161頁第六十五頁，共162頁。2023/2/1366表3.1常用輸入輸出數(shù)據(jù)格式數(shù)據(jù)輸入格式數(shù)據(jù)輸出格式ArcInfoCoverageE00ArcInfoCoverageE00ArcInfoGridE00ArcInfoGridE00ERDASGISERDASGISERDASLANERDASLANShapefileShapefileGenericBinaryGenericBinaryGeoTIFFGeoTIFFTIFFTIFFJPEGJPEG第66頁/共161頁第六十六頁，共162頁。2023/2/1367數(shù)據(jù)輸入：二進制數(shù)據(jù)轉(zhuǎn)換成IMG數(shù)據(jù)

第67頁/共161頁第六十七頁，共162頁。2023/2/1368二進制數(shù)據(jù)轉(zhuǎn)換成單波段IMG數(shù)據(jù)第68頁/共161頁第六十八頁，共162頁。2023/2/1369在數(shù)據(jù)輸入輸出對話框中，通常需要設置下列參數(shù)信息：→確定是輸入數(shù)據(jù)（Import）還是輸出數(shù)據(jù)（Export）→在列表中選擇輸入或輸出數(shù)據(jù)的類型（Type在）→列表中選擇輸入數(shù)據(jù)或輸出數(shù)據(jù)的媒體（Media、CD-ROM、Tape、File）→確定輸入數(shù)據(jù)路徑和文件名（InputFile：*.*）→確定輸出數(shù)據(jù)文件路徑和文件名（OutputFile：*.*）→OK（進入下一級參數(shù)的設置，隨數(shù)據(jù)類型而不同）

第69頁/共161頁第六十九頁，共162頁。2023/2/1370第70頁/共161頁第七十頁，共162頁。2023/2/1371第71頁/共161頁第七十一頁，共162頁。2023/2/1372第72頁/共161頁第七十二頁，共162頁。2023/2/1373將多個單波段圖像組合成多波段圖像Interpreter>Utilities>LayerStack第73頁/共161頁第七十三頁，共162頁。2023/2/1374第74頁/共161頁第七十四頁，共162頁。2023/2/1375二進制數(shù)據(jù)轉(zhuǎn)換成多波段IMG數(shù)據(jù)第75頁/共161頁第七十五頁，共162頁。2023/2/1376第76頁/共161頁第七十六頁，共162頁。2023/2/1377第77頁/共161頁第七十七頁，共162頁。2023/2/1378第78頁/共161頁第七十八頁，共162頁。2023/2/1379TIFF圖像數(shù)據(jù)輸入輸出

TIFF圖像數(shù)據(jù)是非常通用的圖像文件格式，用戶只要在打開圖像文件時將文件類型制定為TIFF格式就可以直接在視窗中顯示TIFF圖像。如果要在圖像解譯器（Interperter）或其他模塊下對圖像做進一步處理操作，依然需要將TIFF文件轉(zhuǎn)換為IMG文件，這種轉(zhuǎn)換非常簡單，只要在打開TIFF的視窗中將TIFF文件另存為IMG文件就可以了。同樣，如果ERDASIMAGING的IMG文件需要轉(zhuǎn)換為GeoTIFF文件，只要在打開IMG圖像文件的視窗中將IMG另存為TIFF文件就可以了。

第79頁/共161頁第七十九頁，共162頁。2023/2/1380輸出JPEG圖像數(shù)據(jù)JPEG圖像數(shù)據(jù)是一種通用的圖像文件格式，ERDAS可以將自己的IMG圖像文件輸出成JPEG圖像文件，供其他圖像處理系統(tǒng)或辦公軟件使用，轉(zhuǎn)換過程如下：ERDAS圖標面板菜單條：Main→Import/Export→數(shù)據(jù)輸入輸出對話框ERDAS圖標面板菜單條：點擊Import/Export圖標→數(shù)據(jù)輸入輸出對話框第80頁/共161頁第八十頁，共162頁。2023/2/1381作業(yè)將examples內(nèi)的dmtm.img轉(zhuǎn)換成TIFF（GeoTIFF）、JPEG，比較二者的不同將nanj*.tif組合成一個多波段圖像（4波段），并記錄轉(zhuǎn)換過程第81頁/共161頁第八十一頁，共162頁。2023/2/1382實習四、數(shù)據(jù)預處理一、預處理簡介二、三維地形表面三、圖像分幅裁剪四、圖像幾何校正五、圖像仿射變換六、航空影像正射校正七、圖像拼接處理六、圖像投影變換第82頁/共161頁第八十二頁，共162頁。2023/2/1383三維地形表面幾何校正圖像拼接投影變換DataPreparation菜單第83頁/共161頁第八十三頁，共162頁。2023/2/1384一、預處理簡介ERDASIMAGING數(shù)據(jù)預處理模塊是由一組實用的圖像數(shù)據(jù)處理工具構(gòu)成，包括三維地形表面（CreateSurface）、圖像分幅裁剪（SubsetImage）、圖像幾何校正（ImageGeometricCorrection）、圖像拼接處理（MosaicImages）、以及圖像投影變換（ReprojectionImages）等，主要是根據(jù)工作區(qū)域的地理特征和專題信息提取的客觀需要，對數(shù)據(jù)輸入模塊中獲取的IMG圖像文件進行范圍調(diào)整、誤差校正、坐標轉(zhuǎn)換等處理，以便進一步開展圖像解譯、專題分類等分析研究。數(shù)據(jù)預處理模塊簡稱DataPreparation或DataPrep，可以通過兩種途徑啟動：ERDAS圖標面板菜單條：Main→DataPreparation→DataPreparation菜單ERDAS圖標面板工具條：點擊DataPrep圖標→DataPreparation菜單第84頁/共161頁第八十四頁，共162頁。2023/2/1385二、三維地形表面三維地形表面工具允許用戶在不規(guī)則空間點的基礎上產(chǎn)生三維地形表面，所支持的輸入數(shù)據(jù)類型包括：ASCII碼點文件、ArcInfo的Coverage點文件和線文件，ERDASIMAGING的注記數(shù)據(jù)層，以及柵格圖像文件IMG。所有輸入數(shù)據(jù)必須具有X、Y、Z值，三維地形表面工具所應用的TIN插值方法，所輸出的是一個連續(xù)的柵格圖像文件。每一個已知的空間點在輸出的地形表面上保持Z值不變，而沒有Z值的空間點，其輸出表面的Z值是基于TIN其周圍的已知點插值計算獲得的。第85頁/共161頁第八十五頁，共162頁。2023/2/13861.啟動三維地形表面（SurfacingCreateSurface）2.定義地形表面參數(shù)（DefiningSurfaceParameters）ReadPoint對話框第86頁/共161頁第八十六頁，共162頁。2023/2/1387在ReadPoint對話框中需要定義下列參數(shù)：→數(shù)據(jù)源文件類型（SourceFileType）：ASCIIFile→數(shù)據(jù)源文件名稱（SourceFileName）：Inpts.dat→OK→打開ImportOptions對話框在ImportOptions欄目，需要定義下列參數(shù)：→選擇字段類型（FieldType）：DelimitedbySeparation（分割字符）→選擇分割字符（SeparationCharactor）：Comma（逗號分割）→每行結(jié)束字符（RowTerminatorCharactor）：ReturnNewLine（DOS）→確定跳過幾行（NumberofRowstoSkip）：0（從頭讀）→點擊InputPreview標簽，進入InputPreview欄目第87頁/共161頁第八十七頁，共162頁。2023/2/1388ImportOptions對話框（FieldDefinition欄目）第88頁/共161頁第八十八頁，共162頁。2023/2/1389從ImportOptions對話框顯示的原始數(shù)據(jù)可知，數(shù)據(jù)文件中的數(shù)據(jù)記錄方式是一行一個點，每一行數(shù)據(jù)包括點號、X坐標、Y坐標、Z坐標（高程值）四個字段，其中點號在此處讀入數(shù)據(jù)時不需要，因此，必須在ImportOptions對話框的ColumMapping中確定X、Y、Z與數(shù)據(jù)文件中字段的對應關系：→OutputColumnName：X對應InputFieldNumber：2→OutputColumnName：Y對應InputFieldNumber：3→OutputColumnName：Z對應InputFieldNumber：4→OK（關閉ImportOptions對話框，讀入數(shù)據(jù)）→數(shù)據(jù)讀到3DSurfacingDataCellArray中第89頁/共161頁第八十九頁，共162頁。2023/2/13903DSurfacing對話框（讀入數(shù)據(jù)之后）第90頁/共161頁第九十頁，共162頁。2023/2/13913.生成三維地形表面（Creating3DSurface）3DSurfacing對話框菜單條：Surface→Surfacing→打開Surfacing對話框3DSurfacing對話框工具條：點擊PerformSurfacing圖標→打開Surfacing對話框Surfacing對話框第91頁/共161頁第九十一頁，共162頁。2023/2/1392在Surfacing對話框中需要設置下列參數(shù)：→輸出文件名稱（OutputFile）：testsurface.img→表面插值方法（SurfacingMethod）：LinearRubberSheeting→輸出文件范圍（OutputCorners）：自動讀取ULX、ULY、LRX、LRY→輸出像元大?。∣utputCellSize）：X：30/Y：30→輸出像元形狀：SquareCells→輸出圖像背景值（BackgroundValue）：0→輸出統(tǒng)計忽略值（IgnoreZeroinOutputStats）：0→輸出數(shù)據(jù)類型（OutputDataType）：Unsigned16bit→OK（關閉Surfacing對話框，執(zhí)行地形表面過程）4.顯示三維地形表面（Display3DSurface）由三維地形表面過程生成的圖像文件，就是ERDAS系統(tǒng)的DEM文件，DEM文件常常用于三維圖像（ImageDrape）的顯示或虛擬地理信息系統(tǒng)（VirtualGIS）操作。所以，可以在Viewer視窗中打開曲面所生成的地形表面文件testsurface.img文件，顯示其二維平面效果或通過ImageInfo查看其定量信息，也可以在其上疊加對應的具有相同投影系統(tǒng)的圖像文件，顯示其三維立體效果。第92頁/共161頁第九十二頁，共162頁。2023/2/1393三、圖像分幅裁剪（SubsetImage）1.規(guī)則分幅裁剪（RectangleSubsetImage）規(guī)則分幅是指裁剪圖像的邊界范圍是一個矩形，通過左上角和右下角兩點的坐標，就可以確定圖像的裁剪位置，整個裁剪過程比較簡單。ERDAS圖標面板菜單條：Main→DataPreparation→DataPreparation菜單→選擇SubsetImage→開SubsetImage對話框ERDAS圖標面板工具條：點擊DataPrep圖標→打開DataPreparation菜單→選擇SubsetImage→打開SubsetImage對話框第93頁/共161頁第九十三頁，共162頁。2023/2/1394SubsetImage對話框。第94頁/共161頁第九十四頁，共162頁。2023/2/1395在SubsetImage對話框中需要設置下列參數(shù)：→輸入文件名稱（InputFile）：Lanier.img→輸出文件名稱（OutputFile）：Lanier_sub.img→坐標類型（CoordinateType）：File→裁剪范圍（SubsetDefinition）：輸入ULX、ULY、LRX、LRY→輸出數(shù)據(jù)類型（OutputDataType）：Unsigned8bit→輸出統(tǒng)計忽略零值：IgnoreZeroinOutputStats→輸出像元波段（SelectLayers）：1:7（表示選擇從第一波段到第七波段）

4,3,2（表示選擇4，3，2四個波段）→OK（關閉SubsetImage對話框，執(zhí)行圖像裁剪）第95頁/共161頁第九十五頁，共162頁。2023/2/13962.不規(guī)則分幅裁剪（PolygonSubsetImage）不規(guī)則分幅裁剪是指裁剪圖像的邊界范圍是個任意多邊形，無法通過左上角和右下角兩點的坐標確定圖像的裁剪位置，而必須事先生成一個完整閉合多邊形區(qū)域，可以是一個AOI多邊形，也可以是ArcInfo的一個PolygonCoverage。針對不同的情況采用不同裁剪過程。⑴AOI多邊形裁剪（PolygonAOISubsetImage）首先在視窗中打開需要裁剪的圖像，并應用AOI工具繪制多邊形AOI，可以將多邊形AOI保存在文件中（*.aoi），也可以暫時不退出視窗，將圖像與AOI多邊形保留在視窗中。DataPreparation菜單→選擇SubsetImage→開SubsetImage對話框點擊DataPrep圖標→打開DataPreparation菜單→選擇SubsetImage→打開SubsetImage對話框第96頁/共161頁第九十六頁，共162頁。2023/2/1397在SubsetImage對話框中需要設置下列參數(shù)：→輸入文件名稱（InputFile）：Lanier.img→輸出文件名稱（OutputFile）：Lanier_sub.img→應用AOI確定裁剪范圍：點擊AOI按鈕→打開選擇AOI（ChooseAOI）對話框ChooseAOI對話框第97頁/共161頁第九十七頁，共162頁。2023/2/1398→在ChooseAOI對話框中確定AOI的來源（AOISource）：File或Viewer→輸出數(shù)據(jù)類型（OutputDataType）：Unsigned8bit→輸出統(tǒng)計忽略零值：IgnoreZeroinOutputStats→輸出像元波段（SelectLayers）：1:7（表示選擇從第一波段到第七波段）4,3,2（表示選擇4，3，2四個波段）→OK（關閉SubsetImage對話框，執(zhí)行圖像裁剪）第98頁/共161頁第九十八頁，共162頁。2023/2/1399⑵ArcInfo多邊形裁剪（PolygonCoverageSubsetImage）如果按照行政區(qū)劃邊界或自然區(qū)劃邊界進行圖像的分幅裁剪，往往是首先利用ArcInfo或ERDAS的Vector模塊繪制精確的邊界多邊形（Polygon），然后以ArcInfo的Polygon為邊界條件進行圖像裁剪。對于這種情況，需要調(diào)用ERDAS其他模塊的功能分兩步完成。第一步：將ArcInfo多邊形轉(zhuǎn)換成柵格圖像文件ERDAS圖標面板菜單條：Main→ImageInterpreter→Utilities→VectortoRaster→打開VectortoRaster對話框ERDAS圖標面板工具條：Main→Vector→VectortoRaster→打開VectortoRaster對話框ERDAS圖標面板工具條：點擊Interpreter圖標→Utilities→VectortoRaster→打開VectortoRaster對話框ERDAS圖標面板工具條：點擊→Vector圖標→VectortoRaster→打開VectortoRaster對話框第99頁/共161頁第九十九頁，共162頁。2023/2/13100VectortoRaster對話框（a）VectortoRaster對話框（b）第100頁/共161頁第一百頁，共162頁。2023/2/13101在VectortoRaster對話框中需要設置下列參數(shù)：→確定矢量文件名稱（InputVectorFile）：zone88→確定輸出文件名稱（OutputRasterFile）：raster.img→確定矢量文件類型（VectorType）：Polygon→轉(zhuǎn)換范圍大?。⊿ubsetDefinition）：ULX、ULY、LRX、LRY→輸出像元大?。≒ixelSize）：30→柵格數(shù)據(jù)類型（OutputPixelDataType）：Unsigned8bit→使用矢量屬性值（AnItemasPixelValue）：ZONE88-ID→柵格文件類型（OutputImageType）：Thematic→是否創(chuàng)建圖層金字塔（CreatePyramidLayers）→輸出統(tǒng)計忽略零值：IgnoreZeroinOutputStatsOK（關閉VectortoRaster對話框，執(zhí)行失柵轉(zhuǎn)換）第101頁/共161頁第一百零一頁，共162頁。2023/2/13102第二步：通過掩模運算（Mask）實現(xiàn)圖像不規(guī)則裁剪ERDAS圖標面板菜單條：Main→ImageInterpreter→Utilities→Mask→打開Mask對話框ERDAS圖標面板工具條：點擊Interpreter圖標→Utilities→Mask→打開Mask對話框Mask對話框第102頁/共161頁第一百零二頁，共162頁。2023/2/13103在Mask對話框中需要設置下列參數(shù)：→確定輸入文件名稱（InputVectorFile）：Lanier.img→確定掩模文件名稱（InputMaskFile）：raster.img→點擊SetRecode設置裁剪區(qū)域內(nèi)新值（NewValue）為1，區(qū)域外取0值→確定掩模區(qū)域做交集運算：Intersection→輸出圖像文件名稱（OutputFile）：mask.img→輸出數(shù)據(jù)類型（OutputDataType）：Unsigned8bit→Ok（關閉Mask對話框，執(zhí)行掩模運算）→實現(xiàn)圖像不規(guī)則裁剪。第103頁/共161頁第一百零三頁，共162頁。2023/2/13104四、圖像幾何校正（GeometricCorrection）為什么要進行圖像幾何校正？在校正中應該注意哪些問題？校正的具體步驟有哪些？第104頁/共161頁第一百零四頁，共162頁。2023/2/13105概念：幾何校正就是將圖像數(shù)據(jù)投影到平面上，使其符合（Conform）地圖投影系統(tǒng)的過程。地理參考：將地圖坐標系統(tǒng)賦予圖像數(shù)據(jù)的過程，稱為地理參考（Geo-referencing）。由于所有地圖投影系統(tǒng)都遵從于一定的地圖坐標系統(tǒng)，所以幾何校正過程包含了地理參考過程。遙感圖像的幾何精校正：利用控制點進行的幾何校正，它是利用一種數(shù)學模型來近似描述遙感圖像的幾何畸變過程，并利用畸變的遙感圖像與標準地圖之間的一些對應點求得這個幾何畸變模型，然后利用此模型進行幾何畸變的校正，這種校正不考慮畸變的具體形成原因，只考慮如何利用畸變模型來校正遙感圖像。第105頁/共161頁第一百零五頁，共162頁。2023/2/13106圖像幾何校正概述⑴圖像幾何校正途徑①數(shù)據(jù)預處理途徑（StartfromDataPreparation）Main→DataPreparation→ImageGeometricCorrection→打開SetGeo-CorrectionInputFile對話框點擊DataPrep圖標→ImageGeometricCorrection→打開SetGeo-CorrectionInputFile對話框SetGeo-CorrectionInputFile對話框第106頁/共161頁第一百零六頁，共162頁。2023/2/13107在SetGeo-CorrectionInputFile對話框中，需要確定校正圖像，有兩種選擇情況其一：首先確定來自視窗（FromViewer），然后選擇顯示圖像視窗（SelectViewer）打開SetGeometricModel對話框選擇幾何校正計算模型（SelectGeometricModel）→OK打開校正模型參數(shù)與投影參數(shù)設置對話框定義校正模型參數(shù)與投影參數(shù)→Apply→CloseSetGeometricModel對話框第107頁/共161頁第一百零七頁，共162頁。2023/2/13108→打開GCPToolReferenceSetup對話框→確定采點模式，采點校正…GCPToolReferenceSetup對話框第108頁/共161頁第一百零八頁，共162頁。2023/2/13109其二：首先確定來自文件（FromImageFile）然后選擇輸入圖像（InputImageFile）打開SetGeometricModel對話框選擇幾何校正計算模型(SelectGeometricModel)→OK打開校正模型參數(shù)與投影參數(shù)設置對話框定義校正模型參數(shù)與投影參數(shù)→Apply→Close打開GCPToolReferenceSetup對話框確定采點模式，采點校正…對于第一種情況，必須事先在一個視窗中打開需要幾何校正的圖像，否則無法進行選擇。第109頁/共161頁第一百零九頁，共162頁。2023/2/13110（2）視窗柵格操作途徑（StartFromViewerRaster）這種途徑是首先在一個視窗中打開需要校正的圖像，然后在柵格操作菜單中啟動幾何校正模塊，具體過程如下：視窗菜單條：Raster→GeometricCorrection選擇幾何校正計算模型（SelectGeometricModel）→OK打開校正模型參數(shù)與投影參數(shù)設置對話框定義校正模型參數(shù)與投影參數(shù)Apply→Close打開GCPToolReferenceSetup對話框確定采點模式，采點校正…幾何正模塊啟動的兩種途徑相比，視窗柵格操作途徑更為直觀簡便。第110頁/共161頁第一百一十頁，共162頁。2023/2/13111其中，多項式變換在衛(wèi)星圖像校正過程中應用較多，在調(diào)用多項式模型時，需要確定多項式的次方數(shù)（Order），通常整景圖像選擇3次方，次方數(shù)與所需要的最少控制點數(shù)是相關的，最少控制點計算公式為(t+1)(t+2)/2，式中t為次方數(shù)，即一次方程需要3個控制點，二次方程需要6個控制點，三次方程需要10個控制點，依此類推。幾何校正計算模型（GeometricCorrectionModel）SPOT衛(wèi)星圖像正射校正SPOTLandsat衛(wèi)星圖像正射校正Landsat航空影像正射校正Camera非線性、非均勻變換RubberSheeting投影變換（轉(zhuǎn)換調(diào)用多項式變換）Reproject多項式變換（同時做投影變換）Polyomial圖像仿射變換（不做投影變換）Affine

功能

模型表4.1幾何校正計算模型與功能第111頁/共161頁第一百一十一頁，共162頁。2023/2/13112表4.2幾何校正采點模式與含義地圖采點模式：在當前數(shù)字化儀上采點在新配置數(shù)字化儀上采點通過鍵盤輸入控制點MaptoViewer：DigitizingTablet（Current）DigitizingTablet（New）KeyboardOnly文件采點模式：在控制點文件中讀點ASCII碼文件中讀點FiletoViewer：GCPFile（*.gcc）ASCIIFile視窗采點模式：在已經(jīng)打開的視窗中采點在新打開的圖像視窗中采點在新打開的矢量視窗中采點在新打開的注記視窗中采點ViewertoViewer：ExistingViewerImageLayer（NewViewer）VectorLayer（NewViewer）AnnotationLayer（NewViewer）含義模式第112頁/共161頁第一百一十二頁，共162頁。2023/2/13113圖像校正的一般流程（GeneralRectificationWorkflow）圖像校正的一般流程第113頁/共161頁第一百一十三頁，共162頁。2023/2/13114圖像校正的具體過程（ConcreteRectificationProcess）第一步：顯示圖像文件（DisplayImageFiles）第二步：啟動幾何校正模塊（GeometricCorrectionTool）PolynomialModelProperties對話框SetGeometricModel對話框第114頁/共161頁第一百一十四頁，共162頁。2023/2/13115GeoCorrectionTools對話框第115頁/共161頁第一百一十五頁，共162頁。2023/2/13116第三步：啟動控制點工具（StartGCPTools）ViewerSelectionInstruction指示器第116頁/共161頁第一百一十六頁，共162頁。2023/2/13117ReferenceMapInformation提示框第117頁/共161頁第一百一十七頁，共162頁。2023/2/13118第四步：采集地面控制點（GroundControlPoint）a.在GCP工具對話框中點擊SelectGCP圖標，進入GCP選擇狀態(tài)；b.在GCP數(shù)據(jù)表中將輸入GCP的顏色（Color）設置為比較明顯的黃色；c.在Viewer#1中移動關聯(lián)方框位置，尋找明顯的地物特征點，作為輸入GCP；d.在GCP工具對話框中點擊CreateGCP圖標，并在Viewer#3中點擊左鍵定點，GCP數(shù)據(jù)表將記錄一個輸入GCP，包括其編號、標識碼、X坐標、Y坐標；e.在GCP工具對話框中點擊SelectGCP圖標，重新進入GCP選擇狀態(tài)；f.在GCP數(shù)據(jù)表中將參考GCP的顏色設置為比較明顯的紅色；第118頁/共161頁第一百一十八頁，共162頁。2023/2/13119地面控制點在圖像上有明顯的、清晰的定位識別標志，如道路交叉點、河流叉口、建筑邊界、農(nóng)田界線等;地面控制點上的地物不隨時間而變化，以保證當兩幅不同時段的圖像或地圖幾何校正時，可以同時識別出來；在相對高差大的地區(qū)盡可能選取高程相似的地段；控制點應均勻分布在整幅圖像內(nèi)，且要有一定的數(shù)量保證；地面控制點的數(shù)量、分布和精度直接影響幾何校正的效果；控制點的精度和選取的難易程度與圖像的質(zhì)量、地物的特征及圖像的空間分辨率密切相關。第119頁/共161頁第一百一十九頁，共162頁。2023/2/13120g.在Viewer#2中移動關聯(lián)方框位置，尋找對應的地物特征點，作為參考GCP；h.在GCP工具對話框中點擊CreateGCP圖標，并在Viewer#4中點擊左鍵定點，系統(tǒng)將自動把參考點的坐標（XReference、YReference）顯示在GCP數(shù)據(jù)表中；i.在GCP工具對話框中點擊SelectGCP圖標，重新進入GCP選擇狀態(tài)，并將光標移回到Viewer#1中，準備采集另一個控制點；j.不斷重復上述步驟（a～i），采集若干GCP，直到滿足所選定的幾何校正模型為止；而后，每采集一個InputGCP，系統(tǒng)就自動產(chǎn)生一個Ref.GCP，通過移動Ref.GCP可以逐步優(yōu)化校正模型，采集GCP后，GCP數(shù)據(jù)表如下圖所示。第120頁/共161頁第一百二十頁，共162頁。2023/2/13121GCP工具對話框GCPTool對話框第121頁/共161頁第一百二十一頁，共162頁。2023/2/13122第五步：采集地面檢查點（GroundCheckPoint）Check（檢查點）：用于檢驗所建立的轉(zhuǎn)換方程的精度和實用性依然在GCP工具對話框狀態(tài)下：在GCP工具對話框菜單條中確定GCP類型：Edit→SetPoint→Type→Check在GCPTool菜單條中確定GCP匹配參數(shù)（MatchingParameter）Edit→PointMatching→打開GCPMatching對話框GCPMatching對話框第122頁/共161頁第一百二十二頁，共162頁。2023/2/13123第六步：計算轉(zhuǎn)換模型（ComputeTransformation）PolynomialModelProperties對話框第123頁/共161頁第一百二十三頁，共162頁。2023/2/13124第七步：圖像重采樣（ResampletheImage）NearestNeighbor：最鄰近點插值法，將最鄰近像元值直接賦予輸出像元。優(yōu)點是算法非常簡單且保持原光譜信息不變；缺點是幾何精度較差，使校正后的圖像亮度具有不連續(xù)性，表現(xiàn)為原來光滑的邊界出現(xiàn)鋸齒狀。BilinearInterpolation：雙線性插值法，用雙線性方程和2×2窗口計算輸出像元值，優(yōu)點是計算較為簡單，且具有一定的亮度采樣精度以及幾何上比較精確，從而使得校正后的圖像亮度連續(xù)；缺點是由于亮度值內(nèi)插，原來的光譜信息發(fā)生了變化，而且這種方法具有低通濾波的性質(zhì)，從而易造成高頻成分(如線條、邊緣等)的損失，使圖像變得模糊。CubicConvolution：立體卷積插值法，用立方方程和4×4窗口計算輸出像元值，優(yōu)點是不僅圖像亮度連續(xù)以及幾何上較精確，而且還能較好的保留高頻部分;其缺點是計算量大。第124頁/共161頁第一百二十四頁，共162頁。2023/2/13125ImageResample對話框第125頁/共161頁第一百二十五頁，共162頁。2023/2/13126第八步：保存幾何校正模式（SaveRectificationMode）在GeoCorrectionTools對話框中點擊Exit按鈕，退出圖像幾何校正過程，按照系統(tǒng)提示選擇保存圖像幾何校正模式，并定義模式文件（*.gms），以便下次直接使用。第九步：檢驗校正結(jié)果（VerityRectificationResult）檢驗校正結(jié)果的基本方法是：同時在視窗中打開兩幅圖像，其中一幅是校正以后的圖像，一幅是當時的參考圖像，通過視窗地理連接（GeoLink/Unlink）功能及查詢光標（InquireCursor）功能進行目視定性檢驗。第126頁/共161頁第一百二十六頁，共162頁。2023/2/13127五、圖像仿射變換在實際工作中，經(jīng)常需要對圖像進行一次線性轉(zhuǎn)換，諸如旋轉(zhuǎn)、位移、翻轉(zhuǎn)、拉伸等，以使圖像的北方向真正向北，這些操作都屬于圖像仿射變換的范疇，具體過程如下：第一步：在視窗中顯示需要的圖像第二步：在視窗中啟動圖像幾何校正工具第三步：選擇幾何校正模型在SetGeometricModel對話框中選擇Affine（仿射變換）→Ok→打開AffineModel對話框第127頁/共161頁第一百二十七頁，共162頁。2023/2/13128AffineModelProperties對話框第128頁/共161頁第一百二十八頁，共162頁。2023/2/13129第四步：定義圖像仿射變換參數(shù)在AffineModelProperties對話框中，定義下列參數(shù)：線性變換選擇項（LinearAdjustmentOptions）比例因子（Scale）：X1:1、Y:1:1、Both:1位移因子（Offset）：X:200、Y:200旋轉(zhuǎn)角度（RotateAngle）：30確定正向選擇方向（PositiveRotationDirection）逆時針方向為正翻轉(zhuǎn)方向選擇（ReflectOptions）：（四種任選）Apply→Close第129頁/共161頁第一百二十九頁，共162頁。2023/2/13130第五步：圖像重采樣在GeoCorrectionTools對話框中選擇ImageResample圖標打開ImageResample（圖像重采樣）對話框在ImageResample對話框中，定義重采樣參數(shù)：輸出圖像文件名（OutputFile）：affine.img選擇重采樣方法（ResampleMethod）：NearestNeighbor定義輸出圖像范圍（OutputCorner）：ULX、ULY、LRX、LRY定義輸出像元大?。∣utputCellSize）：X：30Y：30設置輸出統(tǒng)計中忽略零值：IgnoreZeroinStats設置重采樣計算輸出缺省值（RecalculateOutputDefault）：SkipFactor5OK（關閉ImageResample對話框，啟動重采樣進程）第130頁/共161頁第一百三十頁，共162頁。2023/2/13131第六步：重采樣圖像與原圖像對比同時打開兩個視窗（Viewer#1和Viewer#2），并平鋪排列在Viewer#1中打開經(jīng)過仿射變換以后的圖像，并在RasterOption中將OrientImagetomapSystem選擇項設置為Off狀態(tài)在Viewer#2中打開仿射變換以前的原始圖像將兩個視窗建立地理關聯(lián)，顯示，查看，對比變化。第131頁/共161頁第一百三十一頁，共162頁。2023/2/13132六、航空影像正射校正航空影像正射校正（Ort