遙感數據圖像處理實驗三、遙感圖像的幾何校正與裁剪

實驗三、遙感圖像的幾何校正與裁剪實驗內容：1.圖像分幅裁剪（SubsetImage）2.圖像幾何校正（GeometricCorrection）3.圖像拼接處理（MosaicImgaes）4.生成三維地形表面（3DSurfacing）1.圖像分幅裁剪在實際工作中，經常需要根據研究工作范圍對圖像進行分幅裁剪，按照ERDASIMAGINE8.4實現(xiàn)圖像分幅裁剪的過程，可以將圖像分幅裁剪為兩類型：規(guī)則分幅裁剪，不規(guī)則分幅裁剪。1.1規(guī)則分幅裁剪(以c:\ProgramFile\IMAGINE8.4\examples\lanier.img為例)規(guī)則分幅裁剪是指裁剪圖像的范圍是一個矩形，通過左上角和右上角兩點的坐標可以確定圖像的裁剪位置，過程如下：方法一:→ERDASIMAGINE8.4圖標面板菜單條：Main→DataPreparation(或單擊ERDASIMAGINE8.4圖標面板工具條“DataPrep”圖標)→打開DataPreparation對話框→單擊SubsetImage按鈕，打開Subset對話框在Subset對話框中需要設置下列參數：→輸入文件名（InputFile）:lanier.img→輸出文件名（OutputFile）:lanier_sub.img→坐標類型（CoordinateType）:Map→裁剪范圍（SubsetDefinition）:ULX、ULY、LRX、LRY（注：ULX，ULY是指左上角的坐標，LRX,LRY是指右上角的坐標，缺省狀態(tài)為整個圖像范圍）→輸出數據類型（OutputDataType）:Unsigned8Bit→輸出文件類型（OutputLayerType）:Continuous→輸出統(tǒng)計忽略零值:IgnoreZeroInOutputStats→輸出像元波段（SelectLayers）:2,3,4→OK(關閉Subset對話框，執(zhí)行圖像裁剪)方法二:→ERDASIMAGINE8.4圖標面板菜單條:Main→StartIMAGINEViewer(或單擊RDASIMAGINE8.4圖標面板工具條“Viewer”圖標)→打開一個二維視窗→單擊視窗工具條最左端的“打開文件”圖標→打開SelectLayerToAdd對話框在SelectLayerToAdd對話框完成以下設置：→LookIn:examples→FileName:lanier.img→Filesoftype:IMAGINEImage→雙擊OK按鈕→在二維視窗中打開lanier.img文件→單擊ERDASIMAGINE8.4圖標面板工具條“DataPrep”圖標→打開DataPreparation對話框→單擊SubsetImage按鈕→打開Subset對話框在Subset對話框中需要設置下列參數：→輸入文件名（InputFile）:lanier.img→輸出文件名（OutputFile）:lanier_sub.img→坐標類型（CoordinateType）:Map→輸出數據類型（OutputDataType）:Unsigned8Bit→輸出文件類型（OutputLayerType）:Continuous→輸出統(tǒng)計忽略零值:IgnoreZeroInOutputStats→輸出像元波段（SelectLayers）:2,3,4→單擊FromInquireBox按鈕→打開InvalidCoordinateType對話框→單擊Continue→在顯示圖像文件lanier.img視窗中單擊工具條的“+”按鈕,打開InquireCursor對話框,在視窗中移動十字光標,確定裁剪范圍左上角和右下角,讀取其坐標分別填入SubsetImage對話框的ULX,ULY中和LRX,LRY中→單擊OK按鈕(關閉Subset對話框，執(zhí)行圖像裁剪)方法三:首先在視窗中打開lanier.img文件→AOI→Tools打開AOI工具面板→單擊矩形框確定裁剪范圍→File→Save→AOILayerAs→打開SaveAOIAs對話框,輸入文件名:2→單擊OK(退出SaveAOIAs對話框)→單擊ERDASIMAGINE8.4圖標面板工具條“DataPrep”圖標→打開DataPreparation對話框→單擊SubsetImage按鈕→打開Subset對話框在Subset對話框中需要設置下列參數：→輸入文件名（InputFile）:lanier.img→輸出文件名（OutputFile）:lanier_sub.img→坐標類型（CoordinateType）:Map→輸出數據類型（OutputDataType）:Unsigned8Bit→輸出文件類型（OutputLayerType）:Continuous→輸出統(tǒng)計忽略零值:IgnoreZeroInOutputStats→輸出像元波段（SelectLayers）:2,3,4→單擊AOI按鈕→打開ChooseAOI對話框→在ChooseAOI對話框作如下設置：→AOISource:File→AOIFile:2→單擊OK(退出ChooseAOI對話框)→單擊OK(退出Subset對話框,執(zhí)行圖像裁剪)→單擊OK(退出Modeler對話框,完成圖像裁剪)1.2不規(guī)則分幅裁剪不規(guī)則分幅裁剪是指裁剪圖像的邊界范圍是個任意多邊形，無法通過左上角和右下角兩點的坐標確定圖像的裁剪位置，而必須事先生成一個完整的閉合多邊形區(qū)域，可以是一個AOI多邊形，也可以是ArcInfo的一個PolygonCoverage，針對不同的情況采用不同的裁剪過程。（一）AOI多邊形裁剪(以c:\ProgramFile\IMAGINE8.4\examples\lanier.img為例)→在視窗中打開需要裁剪的圖像→應用AOI工具繪制多邊形可以將多邊形AOI文件保存在文件中（＊.aoi）,也可以暫時不退出視窗，將圖像與AOI多邊形保留在視窗中，然后進行以下操作：→ERDASIMAGINE8.4圖標面板菜單條：Main→DataPreparation(或單擊ERDASIMAGINE8.4圖標面板工具條“DataPrep”圖標)，打開DataPreparation對話框→單擊SubsetImage按鈕，打開Subset對話框在Subset對話框中需要設置下列參數：→輸入文件名（InputFile）:lanier.img→輸出文件名（OutputFile）:lanier_sub.img→CoordinateType:Map→Output:Unsigned8bit→Oupput:Continuous→選定IgnoreZeroinOutputStats→SelectLayers:2,3,4→單擊AOI按鈕，打開ChooseAOI對話框→選定File→選擇文件：2.aoi→單擊OK(退出ChooseAOI對話框OK（退出Subset對話框→單擊Modeler對話框中的OK按鈕完成圖像裁剪（二）ArcInfo多邊形裁剪)→單擊,執(zhí)行圖像裁剪）(以c:\ProgramFile\IMAGINE8.4\examples\gy_resample.img為例)如果是按照行政區(qū)劃邊界或自然區(qū)劃邊界進行圖像的分幅裁剪，往往是首先利用ArcInfo或者ERDAS的Vector模塊繪制精確的邊界多邊形（Polygon），然后以ArcInfo的Polygon為邊界進行圖像裁剪。對于這種情況，需要調用ERDAS其它模塊的功能分兩步完成。第一步：將ArcInfo多邊形轉換為柵格圖像文件打開VectorToRaster對話框有兩種方法：方法一：ERDASIMAGINE8.4圖標面板菜單條：Main→ImageInterpreter(或單擊ERDASIMAGINE8.4圖標面板工具條“Interpreter”圖標)→打開ImageInterpreter對話框→單擊Utilities→打開Utilities對話框→選擇VectorToRaster→打開VectorToRaster對話框方法二：ERDASIMAGINE8.4圖標面板菜單條：Main→Vector(或單擊ERDASIMAGINE8.4圖標面板工具條“Vector”圖標)→打開VectorUtilities對話框→選擇VectortoRaster→打開VectortoRaster對話框本例中選擇前一種方法：在VectortoRaster對話框中設置下列參數：→InputVectorFile（輸入矢量文件名稱）：boundgy.shp→VectorType(確定矢量文件類型)：polygon→UseAttributeasValue(使用矢量屬性值)：GYBOUND_ID→OutputFile（輸出柵格文件名稱）：raster.img→DataType(柵格數據類型)：Unsigned8bit→LayerType(柵格文件類型)：Thematic→SizeDifinition(轉換范圍大小)：ULX,ULY,LRX，LRY（缺省條件下為整個圖像范圍）→Units（坐標單位）：Meters→CellSize（輸出像元大?。篨:30/Y:30→選擇SquareCell（正方形像元）→單擊OK（關閉VectorToRaster對話框，執(zhí)行矢量到柵格的轉換）矢量圖：柵格圖：第二步：通過掩模運算（Mask）實現(xiàn)圖像的不規(guī)則裁剪ERDASIMAGINE8.4圖標面板菜單條：Main→ImageInterpreter(或單擊ERDASIMAGINE8.4圖標面板工具條“Interpreter”圖標)→ImageInterpreter→Utilities→Mask→打開Mask對話框在Mask對話框中設置下列參數：→InputFile(輸入需要裁剪的圖像文件名稱):gy_resample.img→InputMaskFile(輸入掩模文件名稱)：raster.img→點擊SetupRecode按鈕設置裁剪區(qū)域內新值（NewValue）為1，區(qū)域外取0值→選擇Intersection（確定掩模區(qū)域作交集運算）→OutputFile（輸出圖像文件名稱）：mask.img→Output(輸出數據類型)：Unsigned8bit→單擊OK（關閉Mask對話框，執(zhí)行掩模運算）2.圖像幾何校正(GeometricCorrection)幾何校正就是將圖像數據投影到平面上，使其符合地圖投影系統(tǒng)的過程；而將地圖坐標系統(tǒng)賦予圖像數據的過程，稱為地理參考（Georeferencing）。由于所有地圖投影系統(tǒng)都遵包含了地理參考過程。在ERDASIMAGINE8.4中進行圖像幾何校正,通常有兩種途徑啟動幾何校正模塊：（1）數據預處理途徑從于一定的地圖坐標系統(tǒng)，所以幾何校正過程(以c:\ProgramFile\IMAGINE8.4\examples\lanier.img為例)方法一：首先在視窗中打開lanier.img文件→ERDASIMAGINE8.4圖標面板菜單條：Main→DataPreparation(或單擊ERDASIMAGINE8.4圖標面板工具條“DataPrep”圖標)，打開DataPreparation對話框→選定ImageGeometricCorrection選項→打開SetGeoCorrectionInputFile對話框→選定FromViewer→單擊SelectViewer按鈕→打開ViewerSelectionInstructions對話框→在顯示lanier.img圖像的視窗中單擊鼠標左鍵→打開SetGeometricModel對話框方法二：→ERDASIMAGINE8.4圖標面板菜單條：Main→DataPreparation(或單擊ERDASIMAGINE8.4圖標面板工具條“DataPrep”圖標)，打開DataPreparation對話框→選定ImageGeometricCorrection選項→打開SetGeoCorrectionInputFile對話框→選擇FromImageFile→在InputImageFile中選擇需要進行幾何校正的圖像：lanier.img→單擊OK→同時打開一個二維視窗（lanier.img圖像文件顯示在其中）和SetGeometricModel對話框（2）窗柵格操作途徑(以c:\ProgramFile\IMAGINE8.4\examples\lanier.img為例)首先在視窗中打開lanier.img文件視窗菜單條：Raster→GeometricCorrection→打開SetGeometricModel對話框資源衛(wèi)星圖像校正具體過程：①視窗采點模式：如果已經擁有需要校正圖像區(qū)域的數字圖像、或經過校正的圖像、或注記圖層，就可以應用視窗采點模式，直接以數字地圖、或經過校正的圖像、或注記圖層作為地理參考，在視窗中打開該圖層，從中采集控制點。（以c:\ProgramFile\IMAGINE8.4\examples下的panAtlanta.img和tmAtlanta.img為例，前者是作為地理參考的圖像，后者為需要校正的圖像）第一步：顯示圖像文件→在ERDASIMAGINE8.4圖標面板工具條點擊“Viewer”圖標兩次，打開兩個視窗→ERDASIMAGINE8.4圖標面板菜單條：Session→TileViewer然后，在Viewer#1中打開需要校正的LandsatTM圖像：在Viewer#2中打開作為地理參考的校正過的SPOT圖像：tmAtlanta.imgpanAtlanta.img第二步：啟動幾何校正模塊Viewer#1菜單條：Raster→GeometricCorrection→打開SetGeometricModel對話框→選擇多項式幾何校正計算模型：Polynormial→單擊OK按鈕→同時打開GeoCorrectionTools對話框和PolynormialModelProperties對話框在PolynormialModelProperties對話框中，定義多項式模型參數及投影參數→定義多項式次方(Ploynomial)：2，其余取缺省值→單擊Apply按鈕→單擊Close按鈕→打開GCPToolReferenceSetup對話框說明：該實例是采用視窗采點模式，作為地理參考的SPOT圖像已經含有投影信息，所以這里不需要定義投影參數。如果不是采用視窗采點模式，或者參考圖像沒有包含投影信息，則必須在這里定義投影信息，包括投影類型及其對應的投影參數。第三步：啟動控制點工具在GCPToolReferenceSetup對話框中→選擇視窗采點模式：ExistingViewer→單擊OK按鈕（關閉GCPToolReferenceSetup對話框）→打開ViewerSelectionInstruction指示器→在顯示作為地理參考圖像panAtlanta.img的Viewer#2中點擊左鍵→打開ReferenceMapInformation提示框→單擊OK按鈕（關閉ReferenceMapInformation提示框）此時整個屏幕將顯示：兩個主視窗、兩個放大窗口、兩個關聯(lián)方框（分別位于兩個視窗中，指示放大窗口與主視窗的關系）、控制點工具窗口、幾何校正工具。表明控制點工具啟動，進入控制點采集狀態(tài)第四步：采集地面控制點控制點工具窗口簡介：控制點工具窗口（GCPTool）有菜單條、工具條、控制點數據表和狀態(tài)條四個部分組成。GCPTool窗口菜單命令及其功能命令功能File:文件操作：LoadInput調用輸入控制點文件（*.gcc）SaveInput將輸入控制點保存在圖像中SaveInputAs保存輸入控制點文件（*.gcc）LoadReferenceSaveReferenceSaveReferenceASClose調用參考控制點文件（*.gcc）將參考控制點保存到圖像中保存參考控制點文件（*.gcc）關閉控制點工具View:顯示操作：ViewOnlySelectedGCP’sShowSelectedGCPinTableArrangeFramesonScreenTools視窗僅顯示所選擇的控制點在表格顯示所選擇的控制點重新排列屏幕中的組成要素調出控制點工具圖標面板重新打開放大窗口StartChipViewerEdit:編輯操作：SetPointType(Control/Check)ResetReferenceSourceReferenceProjectionPointPrediction設置采集點的類型（控制／檢查）改變參考控制點源文件改變參考文件的投影參數按照轉換方程計算下一點位置匹配兩幅圖像像元的灰度值PointMatchingHelp:聯(lián)機幫助：HelpforGCPTool關于GCP工具的聯(lián)機幫助GCPTool窗口工具圖標命令及其功能(按圖標從左到右的順序)命令功能Toggle自動GCP編輯模式開關鍵Calculate依據控制點求解幾何校正模型設置自動轉換計算開關AutomaticComputeErrorSelectGCPCreateGCPLock/UnlockFindinInputFindinReferZValue計算檢查點的誤差，更新RMS誤差激活GCP選擇工具，在視窗中選擇GCP在視窗中選擇定義GCP鎖?。尫女斍懊?，以便重復使用選擇尋找輸入圖像中的GCP選擇尋找參考文件中的GCP計算更新所選GCP的Z值Auto-ZValue自動更新所有GCP的Z值GCP窗口數據字段及其含義字段含義GCP順序號,系統(tǒng)自動產生GCP標識碼,用戶可以定義GCP當前選擇狀態(tài)提示符號輸入GCP顯示顏色輸入GCP的X坐標輸入GCP的Y坐標參考GCP顯示顏色參考GCP的X坐標參考GCP的Y坐標GCP的類型(控制點／檢查點)單個GCP的X殘差Point#PointID>ColorXInputYInputColorXReferenceYReferenceTypeXResidualYResidualRMSError單個GCP的Y殘差單個GCP的RMS誤差Contribution單個GCP的貢獻率Match兩幅圖像像元灰度值的匹配程度關于GCPTool窗口的說明:1.輸入控制點(InputGCP)是在原始文件視窗中采集的，具有原文件的坐標系統(tǒng)；而參考控制點(ReferenceGCP)是在參考文件視窗中采集的，具有已知的參考坐標系統(tǒng)，GCP工具將根據對應點的坐標值自動生成轉換模型。2.在GCP數據中，殘差(Residuals)、中誤差(RMS)、貢獻率(Contribution)及匹配程度(Match)等參數，是在編輯GCP的過程中自動計算更新的，用戶是不能任意改變的，但可以通過精確GCP位置來調整3.每個IMG文件都可以有一個GCP數據集與之相關聯(lián)，GCP數據集保存在一個柵格層數據文件中；如果IMG文件有一個GCP數據集存在的話，只要打開GCPTool，GCP點都會出現(xiàn)在視窗中4.所有的輸入控制點(InputGCP)都可以保存在圖像文件中(SaveInput)，也可以保存在控制點文件中(SaveInputAs)。如果保存在圖像文件中，調用的方法如3所述，如果是保存在控制點文件中，可以通過加載調用(LoadInput)。5.參考控制點(ReferenceInput)也可以類似地保存在參考圖像中(SaveReference)或控制點文件中(SaveReferenceAs)，可以通過加載調用(LoadReference)采集地面控制點具體過程如下：→在GCPTool窗口中點擊SelectGCP圖標，進入GCP選擇狀態(tài)→在GCP數據表中將輸入GCP的顏色（Color）設置為比較明顯的黃色→在Viewer#1中移動關聯(lián)方框位置，尋找明顯的地物特征點，作為輸入GCP→在GCPTool窗口中點擊CreateGCP圖標，并在Viewer#3中點擊左鍵定點，GCP數據表將記錄一個輸入GCP，包括其編碼、標識碼。X坐標、Y坐標→在GCPTool窗口中點擊SelectGCP圖標，重新進入GCP選擇狀態(tài)→在GCP數據表中將參考GCP的顏色（Color）設置為比較明顯的紅色→在Viewer#2中移動關聯(lián)方框位置，尋找對應的地物特→在GCPTool窗口中點擊CreateGCP圖標，并在Viewer#4中點擊左鍵定點，系統(tǒng)將自動把參考點的坐標（XReference，YReference）顯示在GCP數據表中征點，作為參考GCP→在GCPTool窗口中點擊SelectGCP圖標，重新進入GCP選擇狀態(tài)，并將光標移回到Viewer#1，準備采集另一個輸入控制點→不斷重復以上步驟，采集若干GCP，直到滿足所選定的幾何校正模型為止（多項式變換的次方數為二時，最少需要六個控制點），而后每采集一個InputGCP，系統(tǒng)就自動產生一個ReferenceGCP。第五步：采集地面檢查點以上所采集的GCP的類型(Type)均為Control(控制點),用于控制計算,建立轉換模型及多項式方程.下面所要采集的GCP的類型均是Check(檢查點),用于檢驗所建立的轉換方程的精度和適用性。在GCPTool窗口中→在GCPTool窗口菜單條中確定GCP類型：Edit→SetPointType→Check→在GCPTool窗口菜單條中確定GCP匹配參數（MatchingParameter）：Edit→PointMatching→打開GCPMatching對話框在GCPMatching對話框，需要定義下列參數：MatchingParameter（匹配參數）:→Max.SearchRadius（最大搜索半徑）：3→SearchWindow（搜索窗口大小）：X:5；Y:5ThresholdParameters（約束參數）：→CorrectionThreshold（相關閾值）:0.8→選定DiscardUnmatchedPoint（刪除不匹配的點）→單擊Close（關閉GCPMatching對話框）→確定地面檢查點：在GCPTool工具條中選擇CreateGCP圖標，并將Lock圖標打開鎖住CreateGCP功能，象選擇控制點一樣，分別在Viewer#1和Viewer#2中定義檢查點，定義完畢后點擊Unlock圖標，解除CreateGCP功能→計算檢查點誤差：在GCPTool工具條中點擊“ComputeError”圖標，檢查點的誤差值會顯示在GCPTool的上方，只有所有檢查點的誤差均小于一個像元，才能繼續(xù)進行合理的重采樣第六步：圖像重采樣重采樣過程就是依據未校正圖像像元值計算生成一幅校正圖像的過程，原圖像中所有柵格數據層都將進行重采樣。ERDASIMAGINE8.4提供三種最常用的重采樣方法：NearestNeighbor：鄰近點插值法，將最鄰近像元值直接賦予輸出像元BilinearInterpolation:雙線性插值法，用雙線性方程和2×2窗口計算輸出像元值CubicConvolution:立方卷積插值法,用立方方程和4×4窗口計算輸出像元值具體操作如下:在GeoCorrectionTools對話框中選擇“ImageResample”圖標→打開ImageResample對話框在Resample對話框中,定義重采樣參數:→OutputFile(輸出圖像文件):rectify.img→ResampleMethod(選擇重采樣方法)：NearestNeighbor→OutputCellSize(定義輸出圖像范圍)：ULX、ULY；LRX、LRY→選定IgnoreZeroinStats→RecalculateOutputDefault(設置重新計算缺省值)：SkipFactor:X:1;Y:1→單擊Ok（關閉Resample對話框）第八步：保存幾何校正模型在GeoCorrectionTools對話框中點擊Exit按鈕，退出圖像幾何校正過程。按照系統(tǒng)提示將圖像幾何校正模型保存為模式文件resample.gms.第九步：檢驗校正結果檢驗校正結果的基本方法是：同時在兩個視窗中打開兩幅圖像，其中一幅是校正后的圖像，一幅是當時的參考圖像，通過視窗地理連接（GeoLink／Unlink）功能及查詢光標功能進行目視定性檢驗，具體過程如下：→在ERDASIMAGINE8.4圖標面板具工條點擊“Viewer”圖標兩次，打開兩個視窗→ERDASIMAGINE8.4圖標面板菜單條：Session→TileViewers然后，在Viewer#1中打開校正后的圖像：rectify.img在Viewer#2中打開參考圖像：panAtlanta.img→在Viewer#1中：按住右鍵→快捷菜單→GeoLink/Unlink→在Viewer#2中：點擊左鍵→建立與Viewer#1的連接→在Viewer#1中：按住右鍵→快捷菜單→InquireCursor→打開光標查詢對話框→在Viewer#1中：移動查詢光標,觀測其在兩屏幕中的位置及匹配程度,并注意光標查詢對話框中數據的變化。如果滿意的話，關閉光標查詢窗口。②地圖采點模式：如果只有硬拷貝的地圖或坐標紙作為參考，可以采用地圖采點模式。首先在地圖上選點并量算其坐標，在地圖上尋找六個以上明顯的地物特征點，并量算起坐標（多項式變換的次方數為二時，最少需要六個控制點）,然后通過鍵盤輸入坐標數據。(以c:\ProgramFile\IMAGINE8.4\examples\nj321.img為例)首先在一個二維視窗中打開nj321.img圖像文件

視窗菜單條：Raster→GeometricCorrection→打開SetGeometricModel對話框→選擇多項式幾何校正計算模型：Polynormial(注：衛(wèi)星圖像校正過程中多采用Polynormial（多項式變換）→單擊OK按鈕→同時打開GeoCorrectionTools對話框和PolynormialModelProperties對話框在PolynormialModelProperties對話框中，定義多項式模型參數及投影參數:→定義多項式次方(Ploynomial)：2，其余取缺省值→單擊Apply按鈕→單擊Close按鈕→打開GCPToolReferenceSetup對話框→選擇KeyboardOnly（通過鍵盤輸入控制點）→單擊Ok→打開ReferenceMapInformation對話框→在MapUnits（地圖單位）中選擇：Other（說明：因為本例中采用的地圖單位為km,在該下拉菜單中沒有列出）→單擊OK→打開GCPTool對話框在顯示nj321.img圖像的視窗Viewer#1中產生一個關聯(lián)方框，同時產生一個放大窗口Viewer#2，放大窗口顯示關聯(lián)方框中圖像的放大效果。→在Viewer#1中用關聯(lián)方框查找與地圖對應的一個控制點，單擊GCPTool對話框工具條的“CreateGCP”圖標，在放大窗口Viewer#2中點擊左鍵定點，XInput和YInput值會自動產生，在XRef.和YRef.中輸入控制點在地圖上的坐標→按照類似方五個控制點的坐標→在GCPTool窗口菜單條中確定GCP類型：Edit→SetPointType→Check→當在放大窗口Viewer#2中點擊第七個控制點時，XRef.和YRef.值會自動產生。法輸入其余的左鍵定點→繼續(xù)在Viewer#1中尋找四個明顯的點，在放大窗口Viewer#2中點擊左鍵定點，這五個點作為檢查點→計算檢查點誤差：在GCPTool工具條中點擊顯示在GCPTool的上方，均小于一個像元（即<30米），才能繼續(xù)進直到滿足要求為止?！癈omputeError”圖標，檢查點的誤差值會只有所有檢查點的誤差行合理的重采樣，如果誤差大，則必須重復以上操作→單擊GeoCorrectionTools對話框的“Resample”圖標→打開Resample對話框在Resample對話框中,定義重采樣參數:→OutputFile(輸出圖像文件名):nj321resample.img→ResampleMethod(重采樣方法)：NearestNeighbor（鄰近點插值法）→選定IgnoreZeroinStats→其它設置均取系統(tǒng)默認值→單擊OK（執(zhí)行重采樣）3.圖像拼接處理（以c:\ProgramFile\IMAGINE8.4\examples下的wasia1_mss.img、wasia2_mss.img、wasia3_mss.img）圖像拼接處理是將具有地理參考的若干相鄰圖像合并成一幅圖像或一圖組像，需要拼接的輸入圖像必須含有地圖投影信息，或者說輸入圖像必須經過幾何校正處理或進行過校正標定。具體操作過程如下：→ERDASIMAGINE8.4圖標面板菜單條：Main→DataPreparation(或單擊ERDASIMAGINE8.4圖標面板工具條“DataPrep”圖標)→打開DataPreparation對話框→MosaicTool菜單→打開MosaicTool對話框MosaicTool視窗菜單條：Edit→AddImage→打開AddImageforMosaic對話框在AddImageforMosaic對話框中，需要設置下列參數：→ImageFilename（拼接圖像文件）：→ImageAreaOptions(圖像拼接區(qū)域)：ComputeActiveAreawasia1_mss.img→單擊Add（圖像wasia1_mss.img被加載到Mosaic視窗中）→按照同樣的過程加載wasia2_mss.img和wasia3_mss.img→單擊Close（關閉AddImageforMosaic對話框）→MosaicTool視窗菜單條：Edit→ImageMatching→打開MatchingOptions對話框→在MatchingOptions對話框中作如下設置：→在MatchingMethod中選擇OverlapArea→單擊OK按鈕（關閉MatchingOptions對話框）→MosaicTool視窗菜單條：Edit→SetOverlapFunction→打開SetOverlapFunction對話框在SetOverlapFunction對話框，設置以下參數：→在IntersectionType中設置NoCutlineExits（沒有裁切線）→在SelectFunction中設置Average（均值）→單擊Apply→單擊Close（關閉SetOverlapFunction對話框）→MosaicTool視窗菜單條：Process→RunMosaic→打開RunMosaic對話框在RunMosaic對話框中設置下列參數：→OutputFileName（確定輸出文件名）：wasia_mosaic.img→WitchOutputs（確定輸出圖形區(qū)域）：All→在IgnoreInputValue中設置0→在OutputBackgroundValue中設置0→選定StatsIgnoreValue，并設置為0→單擊OK（關閉RunMosaic對話框，運行圖像拼接）→MosaicTool視窗菜單條：File→Close→打開VerifySaveonClose對話框→單擊“否”按鈕(關閉MosaicTool視窗,退出Mosaic工具)4.生成三維地形表面三維地形表面工具允許在不規(guī)則空間點的基礎上產生三維地形表面，括：ASCII碼點文件、ArcInfo的Coverage點文件和線文件、ERDASIMAGINE的注記數據層以及柵格圖像文件IMG。所支持的類型包①以ASCII碼點文件為例具體操作過程如下：(數據源以c:\ProgramFile\IMAGINE8.4\examples\lnpts.dat為例)第一步：啟動三維地形表面ERDASIMAGINE8.4圖標面板菜單條：Main→DataPreparation(或單擊ERDASIMAGINE8.4圖標面板工具條“DataPrep”圖標)，→打開DataPreparation對話框→單擊CreateSurface選項→打開3DSurfacing對話框第二步：定義地形表面參數→3DSurfacing對話框菜單條：File→Read→打開ReadPoints對話框→在SourceFileType中設置：ASCIIFile→在SourceFileName中設置：inpts.dat→單擊OK（關閉ReadPoints對話框）→打開ImportOptions對話框在ImportOptions對話框中，需要設置如下參數：→FieldType（字段類型）：選擇DelimitedbySeparator(分割字符)→SeparatorCharacter(分割字符)：Comma（逗號分割）→RowTerminatorCharacter(每行結束字符)：ReturnNewLine（DOS）→NumberofToSkip(跳過行數)：0→點擊InputPreview選項，進入InputPreview欄從ImportOptions對話框（InputPreview欄）顯示的原始數據可以看到數據的記錄方式是一行一個點，每一行數據包括點號、X坐標、Y坐標、高程值四個字段，每兩個字段之間的間隔號在此處讀入時不需要，因此，必須在ImportOptions對話框的ColumnMappingZ與數據文件中字段的對應關系：→OutputColumnName：X對應InputFieldNumber:2→OutputColumnName：Y對應InputFieldNumber:3→OutputColumnName：Z對應InputFieldNumber:4→單擊Ok（關閉ImportOptions對話框，讀入數據）中確定X、Y、→返回3DSurfacing對話框第三步：生成三維地形表面→3DSurfacing對話框菜單條：Surface→Surfacing→打開Surfacing對話框在Surfacing對話框中，設置如下參數：→OutputFile（輸出文件名）：surface.img→SurfacingMethod(表面插值方法)：LinearRubberSheeting（線性插值）→OutputCorners（輸