遙感數(shù)據(jù)圖像處理實(shí)驗(yàn)三、遙感圖像的幾何校正與裁剪

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

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