遙感信息技術(shù)實(shí)習(xí)指導(dǎo)書(shū)

1、實(shí)習(xí)2. ERDAS IMAGINE 8.6系統(tǒng)的數(shù)據(jù)預(yù)處理本章主要學(xué)習(xí)以下內(nèi)容：4-三維地形表面圖像分幅裁減4-圖像幾何校正(此部份內(nèi)容較多，供有興趣的同學(xué)選做)4-圖像拼接處理I 圖像投影變換/、三維地形表面三維地形表面工具允許您在不規(guī)則空間的基礎(chǔ)上產(chǎn)生三維地形表面啟動(dòng)三維地形表面(Starting Create Surface)ERDAS 圖標(biāo)面板菜單條：MainData PreparationData preparation 菜單選擇 Create Surface打開(kāi) 3D Surfacing 對(duì)話(huà)框定義地形表面參數(shù)(Defining Surface Parameters)第一步：3D

2、 Surfacing 對(duì)話(huà)框菜單條：FileRead打開(kāi) Read Points第二步：在Read Points對(duì)話(huà)框中需要定義下列參數(shù)：-數(shù)據(jù)源文件類(lèi)型(Source File Type): ASCII File-數(shù)據(jù)源文件名稱(chēng)(Source File Name):用戶(hù)文件夾Inpts.dat-OK 進(jìn)入 Import Option 對(duì)話(huà)框其中參數(shù)定義如圖所示：然后點(diǎn)擊OK從Import Options對(duì)話(huà)框顯示的原始數(shù)據(jù)可知，數(shù)據(jù)文件中的數(shù)據(jù)記錄方式是一行 一個(gè)點(diǎn)，每一行數(shù)據(jù)包括點(diǎn)號(hào)M坐標(biāo)、Y坐標(biāo)、高程值四個(gè)字段，其中點(diǎn)號(hào)在此外讀 入數(shù)據(jù)時(shí)不需要，因而在Column Mapping中確定X

3、、Y、Z與數(shù)據(jù)文件中字段的對(duì)應(yīng)關(guān) 系為2、3、4。I生成三維地形表面(C reating 3D Surface3D Surfacing 對(duì)話(huà)框菜單條：SurfaceSurfacing打開(kāi) Surfacing 對(duì)話(huà)框其中參數(shù)如圖所示：將輸出文件放在個(gè)人工作目錄下保存，并截圖保存在作業(yè)中說(shuō)明；2、圖像分幅裁減在實(shí)際工作中，經(jīng)常需要根據(jù)研究工作范圍對(duì)圖像進(jìn)行分幅裁剪，按照ERDAS實(shí) 現(xiàn)圖像分幅裁剪的過(guò)程，可以將圖像分幅裁剪分為兩種類(lèi)型:規(guī)則分幅裁剪，不規(guī)則分幅 裁剪。L 規(guī)則分幅裁剪(Rectangle Subset Image)規(guī)則分幅裁剪是指裁剪圖像的邊界范圍是一個(gè)矩形，通討左下角和右下角兩點(diǎn)

4、的坐 標(biāo)，就可以確定圖像的裁剪位置，整個(gè)裁剪過(guò)程比較簡(jiǎn)單。ERDAS 圖標(biāo)面板菜單條:Main 一 Data Preparation 一 Data Preparation 菜單 選擇 Subset Image 一打開(kāi) Subset Image 對(duì)話(huà)框圖 2.1 Subset Image 對(duì)話(huà)框選擇 Subset Image 一打開(kāi) Subset Image 對(duì)話(huà)框(圖 2.1)在Subset Image對(duì)話(huà)框中需要設(shè)置下列參數(shù)：輸入文件名稱(chēng)(Input File):Lanier.img輸出文件名稱(chēng)(Output File):Lanier sub.img坐標(biāo)類(lèi)型(Coordinate Type)

5、:File裁剪范圍(Subset Definition)輸入 ULX, ULY LRX, LRY (用戶(hù)自行決定)輸出數(shù)據(jù)類(lèi)型(Output Data Type” ):Unsigned 8 Bit輸出文件類(lèi)型(Oouput Layer Type) : Continuous輸出統(tǒng)計(jì)忽略零值:Ignore Zero In Output Stats輸出像元波段(Select Layers):2:5(表示選擇2, 3, 4, 5四個(gè)波段)一 OK(關(guān)閉Subset Image對(duì)話(huà)框，執(zhí)行圖像裁剪)說(shuō)明:在上述圖像裁剪過(guò)程中，裁剪范圍是通過(guò)直接輸入左上角點(diǎn)坐標(biāo)(ULX, ULY) 和右下角點(diǎn)坐標(biāo)(LRX

6、; LRY)定義的。此外，還可以通過(guò)兩種方式定義裁剪范圍:其一是應(yīng) 用查詢(xún)框(Inquire Box)，具體過(guò)程是首先在打開(kāi)被裁剪圖像的視窗中放置查詢(xún)框，然后在 Subsetimage對(duì)話(huà)框中選擇From Inquire Box功能:其二是應(yīng)用感興趣區(qū)域(AOI)，具體過(guò) 程是首先在打開(kāi)被裁剪圖像的視窗中繪畫(huà)矩形AOI,然后在Subset Image對(duì)話(huà)框中選擇 AOI功能，打開(kāi)AOI對(duì)話(huà)框，并確定AOI區(qū)域來(lái)自圖像視窗即可。(有興趣的同學(xué)可以 試試)不規(guī)則分幅裁剪(Polygon Subset Image)在實(shí)際工作中，我們經(jīng)常遇到的問(wèn)題是不規(guī)則的邊界剪裁問(wèn)題。不規(guī)則分幅裁剪是 指裁剪圖像的

7、邊界范圍是個(gè)任意多邊形，無(wú)法通過(guò)左上角和右下角兩點(diǎn)的坐標(biāo)確定圖像 的裁剪位置，而必須事先生成一個(gè)完整的閉合多邊形區(qū)域，可以是一個(gè)AOI(AREA OF INTREST)一多邊形，也可以是ArcInfo的一個(gè)Polygon Coverage,針對(duì)不同的情況采用不 同裁剪過(guò)程.1、AOI 多邊形裁剪(Polygon AOI Subset Image )首先在視窗中打開(kāi)需要裁剪的圖像Lanier.img，并應(yīng)用AOI工具繪制多邊形AOI可以將多邊形AOI保存在文件中(*.aoi )，保存名為：subset.aoi也可以暫時(shí)不退出視窗， 將圖像與AOI多邊形保留在視窗中，然后：ERDAS 圖標(biāo)面板菜單

8、條:MainData Preparation 一 Data Preparation 菜單，選擇 Subset Image 一打開(kāi) Subset Image 對(duì)話(huà)框ERDAS圖標(biāo)面板工具條:點(diǎn)擊Data Prep圖標(biāo)一打開(kāi)Data Preparation菜單選擇Subset Image 一打開(kāi)Subset Image對(duì)話(huà)框在Subset Image對(duì)話(huà)框中需要設(shè)置下 列參數(shù)輸入文件名稱(chēng)(Input File): Lanier.img輸出文件名稱(chēng)(Output File): subLanier.img應(yīng)用AOI確定裁剪范圍:點(diǎn)擊AOI按鈕打開(kāi)選擇AOI (Choose AOI)對(duì)話(huà)框在 Choos

9、e AOI 對(duì)話(huà)框中確定 AOI 的來(lái)源(AOI Source) : File(或 Viewer)如果選擇了文件(File)，則進(jìn)一步確定AOI文件，否則，直接進(jìn)入下一步圖2.4設(shè)置Subset Image參數(shù)對(duì)話(huà)框輸出數(shù)據(jù)類(lèi)型(Output Data Type):Unsigned 8 Bit輸出像元波段(Select Layers):1:6(表示選擇1,2, 3, 4, 5,6六個(gè)波段)OK(關(guān)閉Subset Image對(duì)話(huà)框，執(zhí)行圖像裁剪)將裁剪后的圖像保存在作業(yè)中說(shuō)明ArcInfo 多邊形裁剪(Polygon Coverage Subset Image)如果是按照行政區(qū)劃邊界或自然區(qū)劃邊

10、界進(jìn)行圖像的分幅裁剪，往往是首先利用 ArcInfo或ERDAS的Vector模塊繪制精確的邊界多邊形(Polygon )，然后以.ArcInfo的 Polygon為邊界條件進(jìn)行圖像裁剪。對(duì)于這種情況，需要調(diào)用ERDAS其它模塊的功能分 兩步完成。第一步:將ArcInfo多邊形轉(zhuǎn)換成柵格圖像文件(有4種方式可以打開(kāi))ERDAS 圖標(biāo)面板菜單條:Main 一 Image Interprete: 一 UtilitiesVxtor to Raster打開(kāi)Victor to Raster對(duì)話(huà)框ERDAS，圖標(biāo)面板菜單條:Main 一 Vector 一 Victor to Raster打開(kāi)Victor

11、to Raster對(duì)話(huà)框ERDAS 圖標(biāo)面板工具條:點(diǎn)擊 Interpreter 圖標(biāo)一 Utilities 一 Victor to Raster 打開(kāi)Victor to Raster對(duì)話(huà)框ERDAS圖標(biāo)面板工具條:點(diǎn)擊Vector圖標(biāo)一 Vector to Raster 打開(kāi)Vector to Raster對(duì)話(huà)框在Vector to Raster對(duì)話(huà)框中需要設(shè)置下列參數(shù)：輸入矢量文件名稱(chēng)(Input Vector File):zone88確定矢量文件類(lèi)型(Vector Type) : Polygon輸出柵格文件名稱(chēng)(Output Image File):raster.img柵格數(shù)據(jù)類(lèi)型(Da

12、ta Type):Unsigned 8 bit柵格文件類(lèi)型(Layer Type) : Continuous轉(zhuǎn)換范圍大小(Size Definition) : ULX, ULY LRX, LRY坐標(biāo)單位(Units) :feet輸出像元大小(Cell Size):X:80/Y:80選擇正方形像元:Squire CellOK(關(guān)閉Vector to Raster對(duì)話(huà)框，執(zhí)行矢柵轉(zhuǎn)換)第二步:通過(guò)掩膜運(yùn)算(Mask)實(shí)現(xiàn)圖像不規(guī)則裁剪ERDAS 圖標(biāo)面板菜單條:Main 一 Image Interpreter 一 Utilities 一 Mask打開(kāi) Mask對(duì)話(huà)框ERDAS圖標(biāo)面板工具條:點(diǎn)擊I

13、nterpreter圖標(biāo)一 UtilitiesMask打開(kāi)Mask對(duì)話(huà)框在Mask對(duì)話(huà)框中需要設(shè)置下列參數(shù)：輸入圖像文件名稱(chēng)(Input File) :germtm .img輸入掩膜文件名稱(chēng)(Input Mask File):raster.img點(diǎn)擊Setup Recode設(shè)里裁剪區(qū)域內(nèi)新值(New Value)為1，區(qū)域外取0值確定掩膜區(qū)域做交集運(yùn)算:Intersection輸出圖像文件名稱(chēng)(Output File):mask.img輸出數(shù)據(jù)類(lèi)型(Output Data Type):Unsigned 8 bitOK(關(guān)閉Mask對(duì)話(huà)框，執(zhí)行掩膜運(yùn)算)打開(kāi)多邊形剪切后的圖像mask.img將多

14、邊形剪切后的圖像mask.img截圖保存在作業(yè)中說(shuō)明3、圖像幾何校正(Geometric Correction)幾何校正就是將圖像數(shù)據(jù)投影到平面上，使其符合(Conform)地圖投影系統(tǒng)的過(guò)程; 而將地圖坐標(biāo)系統(tǒng)賦予圖像數(shù)據(jù)的過(guò)程，稱(chēng)為地理參考(Geo-referencing)。由于所有地 圖投影系統(tǒng)都遵從于一定的地圖坐標(biāo)系統(tǒng)，所以幾何校正過(guò)程包含了地理參考過(guò)程。圖像幾何校正途徑(Geometric Correction Process)在ERDAS IMAGINE系統(tǒng)中進(jìn)行圖像幾何校正，通常有兩種途徑啟動(dòng)幾何校正模塊。(1)數(shù)據(jù)預(yù)處理途徑(Start from Data Preparatio

15、n)ERDAS 圖標(biāo)面板菜單條:Main 一 Data Preparation 一 Image Geometric Correction打開(kāi) Set Geo-Correction Input File 對(duì)話(huà)框ERDAS 圖標(biāo)面板工具條:點(diǎn)擊 Data Prep 圖標(biāo)-Image Geometric Correction打開(kāi) Set Geo-Correction Input File 對(duì)話(huà)框上 幾何校正計(jì)算模型（Geometric Correction Model）如圖所示，ERDAS提供的圖像幾何校正計(jì)算模型有10種，具體功能如表所列:幾何校正計(jì)算模型與功能模型功能Affine圖像仿射變換（不

16、做投影變換）Camera航空影像正射校正Landsat衛(wèi)星圖像正射校正Polynomial多項(xiàng)式變換（同時(shí)做投影變換）Rubber Sheeting非線(xiàn)性、非均勻變換SpotSpot衛(wèi)星圖像正射校正DPPDB運(yùn)用合理的多項(xiàng)式系數(shù)描述在圖像獲取時(shí)刻，圖像和地球表面的關(guān)系IKONOS運(yùn)用合理的多項(xiàng)式系數(shù)描述在圖像獲取時(shí)刻，IKONOS圖像和地球表面的關(guān)系NITF RPC運(yùn)用合理的多項(xiàng)式系數(shù)描述在圖像獲取時(shí)刻，NITF格式的圖像和地球表面的關(guān)系QuickBird RPC運(yùn)用合理的多項(xiàng)式系數(shù)描述在圖像獲取時(shí)刻，QuickBird圖像和地球表面的關(guān)系注：dppdb Digital Point Posit

17、ioning Data Base以資源衛(wèi)星圖像校正為例介紹的是以已經(jīng)具有地理參考的SPOT圖像為基礎(chǔ)，進(jìn)行Landsat TM圖像校正過(guò)程。 第一步：顯示圖像文件在首先，在ERDAS圖標(biāo)面板中點(diǎn)擊Viewer圖標(biāo)兩次，打開(kāi)兩個(gè)視窗Viewer #I/ viewer #2 ），并將兩個(gè)視窗平鋪放置，操作過(guò)程如下：ERDAS 圖標(biāo)面板菜單條:Session Tile Viewers然后，在Viewer #仲打開(kāi)需要校正的Landsat TM圖像:tmAtlanta.img在Viewer #2中打開(kāi)作為地理參考的校正過(guò)的SPOT圖像:panAtlanta.img:第二步:啟動(dòng)幾何校正模塊(Geome

18、tric Correction Tool)Viewer #1 菜單條:Raster Geometric Correction 打開(kāi) Set Geometric Model 對(duì)話(huà)框選擇多項(xiàng)式幾何校正計(jì)算模型:Polynomial - OK同時(shí)打開(kāi) Polynomial Model Properties 對(duì)話(huà)框和 Geo Correction Tools 對(duì)話(huà)框在Polynomial Model Properties對(duì)話(huà)框中，定義多項(xiàng)式模型參數(shù)及投影參數(shù)：定義多項(xiàng)式次方(Polynomial Order) : 2定義投影參數(shù)(Projection略)一 Apply 一 Close一打開(kāi) GCP T

19、ool Reference Setup對(duì)話(huà)框說(shuō)明:該實(shí)例是采用視窗采點(diǎn)模式，作為地理參考的SPOT圖像已經(jīng)含有投影信息，所以， 這里不需要定義投影參數(shù)。如果不是采用視窗采點(diǎn)模式，或者參考圖像沒(méi)有包含投影信 息，則必須在這里定義投影信息，包括投影類(lèi)型及其對(duì)應(yīng)的投影參數(shù)。第三步:啟動(dòng)控制點(diǎn)工具(Start GCP Tools )首先在GCP Toots Reference Setup對(duì)話(huà)框中選擇采點(diǎn)模式:一選擇視窗采點(diǎn)模式:Existing Viewer一 OK(關(guān)閉 GCP Tools Reference Setup 對(duì)話(huà)框)一打開(kāi) Viewer Selection Instiuctions 指

20、示器圖 2.12 Viewer Selection Instiuctions 指示器一在顯示作為地理參考圖像panAtlanta.img的Viewer #2中點(diǎn)擊左鍵一打開(kāi)Reference Map Information提示框(圖2.13)(顯示參考圖像的投影信息)一 OK(關(guān)閉 Reference Map Information 提示框)一整個(gè)屏幕將自動(dòng)變化為如圖所示的狀態(tài):其中包含兩個(gè)主視窗、兩個(gè)放大窗口、兩個(gè)關(guān)聯(lián)方框(分別位于兩個(gè)視窗中，指示放大窗口與主視窗的關(guān)系)、控制點(diǎn)工具對(duì)話(huà) 框、幾何校正工具等。表明控制點(diǎn)工具被啟動(dòng)，進(jìn)入控制點(diǎn)采集狀態(tài)。第四步:采集地面控制點(diǎn)(Ground Co

21、ntrol Point)控制點(diǎn)工具對(duì)話(huà)框簡(jiǎn)介(Introduction to GCP Tool)在正式開(kāi)始采集控制點(diǎn)之前，首先對(duì)控制點(diǎn)工具對(duì)話(huà)框進(jìn)行說(shuō)明 GCP工具對(duì)話(huà)框(GCP Tool)由菜單條、工具條、控制點(diǎn)數(shù)據(jù)表(GCP Cell Array)及狀態(tài)條(Status Bar)4個(gè)部分組成，菜單條中菜單命令及其功能詳見(jiàn)軟件菜單說(shuō)明。關(guān)于GCP工具對(duì)話(huà)框，還需要說(shuō)明幾點(diǎn)：0)輸入控制點(diǎn)(Input GCP)的是在原始文件視窗中采集的，具有原文件的坐標(biāo)系統(tǒng)；而 參考控制點(diǎn)(Reference GCP)是在參考文件視窗中采集的，具有己知的參考坐標(biāo)系統(tǒng)，GCP工具將根據(jù)對(duì)應(yīng)點(diǎn)的坐標(biāo)值自動(dòng)生成轉(zhuǎn)換

22、模型。在GCP數(shù)據(jù)表中，殘差(Residuals)、中誤差(RMS )、貢獻(xiàn)率(Contribution)及匹配程 度(Match)等參數(shù)，是在編輯GCP的過(guò)程中自動(dòng)計(jì)算更新的，用戶(hù)是不可以任意改變的， 但可以通過(guò)精確GCP位置來(lái)調(diào)整。每個(gè)IMG文件都可以有一個(gè)GCP數(shù)據(jù)集與之相關(guān)聯(lián)，GCP數(shù)據(jù)集保存在一個(gè)柵格 層數(shù)據(jù)文件中;如果IMG文件有一個(gè)GCP數(shù)據(jù)集存在的話(huà)，只要打開(kāi)GCP工具，GCP 點(diǎn)就會(huì)出現(xiàn)在視窗中。所有的輸入GCP都可以直接保存在圖像文件中(Save Ii)put)，也可以保存在控制點(diǎn) 文件中(Save Input As)。如果是保存在文件中，調(diào)用的方法如(c)所述，如果是保存

23、在GCP 文件中，可以通過(guò)加載調(diào)用(Load Input)。參考GCP也可以類(lèi)似地保存在參考圖像中(Save Reference)或 GCP文件中(SaveReference As)，便于以后調(diào)用。(2) GCP 的具體采集過(guò)程(Steps of GCP Selection)在圖像幾何校正過(guò)程中，采集控制點(diǎn)是一項(xiàng)非常重要和相當(dāng)繁重的工作，具體過(guò)程 如下：在GCP工具對(duì)話(huà)框中點(diǎn)擊Select GCP圖標(biāo)，進(jìn)入GCP選擇狀態(tài)；在GCP數(shù)據(jù)表中將輸入GCP的顏色(Color)設(shè)置為比較明顯的黃色，-在Viewer #1中移動(dòng)關(guān)聯(lián)方框位置，尋找明顯的地物特征點(diǎn)，作為輸入GCP:在GCP工具對(duì)話(huà)框中點(diǎn)擊

24、Create GCP圖標(biāo)，并在Viewer #3中點(diǎn)擊左鍵定點(diǎn)，GCP.數(shù)據(jù)表將記錄一個(gè)輸入GCP，包括其編號(hào)、標(biāo)識(shí)碼、X坐標(biāo)、Y坐標(biāo)：在GCP工具對(duì)話(huà)框中點(diǎn)擊Select GCP圖標(biāo)，重新進(jìn)入GCP選擇狀態(tài)；在GCP數(shù)據(jù)表中將參考GCP的顏色(Color)設(shè)置為比較明顯的紅色：在Viewer #2中移動(dòng)關(guān)聯(lián)方框位置，尋找對(duì)應(yīng)的地物特征點(diǎn)，作為參考GCP:在GCP工具對(duì)話(huà)框中點(diǎn)擊Create GCP圖標(biāo)，并在Viewer #4中點(diǎn)擊左鍵定點(diǎn)，系 統(tǒng)將自動(dòng)把參考點(diǎn)的坐標(biāo)X Reference, Y Reference 顯示在GCP表中，在GCP工具對(duì)話(huà)框中點(diǎn)擊Select GCP圖標(biāo)，重新進(jìn)入G

25、CP選擇狀態(tài):并將光標(biāo)移回 到Viewer #1，準(zhǔn)備采集另一個(gè)輸入控制點(diǎn)。(J)不斷重復(fù)(a) - (i)，采集若干GCP,直到滿(mǎn)足所選定的幾何校正模型為止；而后，每采 集一個(gè)Input GCP,系統(tǒng)自動(dòng)產(chǎn)生一個(gè)Ref. GCP,通過(guò)移動(dòng)Ref.GCP可以逐步優(yōu)化校 正模型。采集GCP以后，GCP數(shù)據(jù)表如圖2.15而示(實(shí)際數(shù)據(jù)以所采集的GSP的點(diǎn) 位不同而不同。)GCP Tool對(duì)話(huà)框與GCP數(shù)據(jù)表第5步采集地面檢查點(diǎn)以上所采集的GCP的類(lèi)型為Control Point(控制點(diǎn))，用于控制計(jì)算，建立轉(zhuǎn)換模型及多 項(xiàng)式方程。下面所要采集的GCP的類(lèi)型均是Check(檢查點(diǎn))，用于檢驗(yàn)所建立的

26、轉(zhuǎn)換方 程的精度和實(shí)用性。依然在GCP Tool對(duì)話(huà)框狀態(tài)下：在 GCP Tool菜單條中確定 GCP類(lèi)型:Edit Set Point TypeCheck在GCP Tool菜單條中確定GCP匹配參數(shù)(Matching Parameter):Edit - Point Matching 一打開(kāi) GCP Matching 對(duì)話(huà)框(略)在GCP Matching對(duì)話(huà)框中，需要定義下列參數(shù)：匹配參數(shù)(Matching Parameters):最大搜索半徑(Max. Search Radius):3搜索窗口大小(Search Window Size): X:5 Y5(4)約束參數(shù)(Threshold P

27、arameters ):設(shè)相關(guān)閾值(Correlation Threshold):0.8.刪除不匹配的點(diǎn)(Discard Unmatched Point):Active在匹配所有/選擇點(diǎn)(Match All / Selected Point)選項(xiàng)組中設(shè)從輸入到參考(Reference from Input)或從參考到輸入(Input from Reference)(6)Close(關(guān)閉 GCP Matching 對(duì)話(huà)框)GCP Matching 對(duì)話(huà)框(c)確定地面檢查點(diǎn):在GCP Tool工具條中選擇Create GCP圖標(biāo)，并將Lock圖標(biāo)打開(kāi)， 鎖住Create GCP功能，如同選擇控制

28、點(diǎn)一樣，分別在Viewer #l和Viewer #2中定義5個(gè) 檢查點(diǎn)，定義完畢后點(diǎn)擊Unlock圖標(biāo)，解除Create GCP功能。(d)計(jì)算檢查點(diǎn)誤差:在GCP Tool工具條中點(diǎn)擊Compute Error圖標(biāo)，檢查點(diǎn)的誤差就會(huì)顯 示在GCP Tool的上方，只有所有檢查點(diǎn)的誤差均小于一個(gè)像元(Pixel)，才能繼續(xù)進(jìn)行合 理的重采樣。一般來(lái)說(shuō)，如果你的控制點(diǎn)(GCP)定位選擇比較準(zhǔn)確的話(huà)，檢查點(diǎn)匹配會(huì)比 較好，誤差會(huì)在限差范圍內(nèi)。否則，若控制點(diǎn)定義不精確，檢查點(diǎn)就無(wú)法匹配，誤差會(huì) 超標(biāo)。第六步:計(jì)算轉(zhuǎn)換模型(Compute Transformation):在控制點(diǎn)采集過(guò)程中，一般是設(shè)置

29、為自動(dòng)轉(zhuǎn)換計(jì)算模式，所以，隨著控制點(diǎn)采集過(guò)程 的完成，轉(zhuǎn)換模型就自動(dòng)計(jì)算生成，下面是轉(zhuǎn)換模型的查閱過(guò)程：_在 Geo-Correction Tools 對(duì)話(huà)框中點(diǎn)擊 Display Model Properties 圖標(biāo)口打開(kāi)Polynomial Model Properties(多項(xiàng)式模型參數(shù))對(duì)話(huà)框(圖4.22)在多項(xiàng)式模型參數(shù)對(duì)話(huà)框中查閱模型參數(shù)，并記錄轉(zhuǎn)換模型Close (關(guān)閉模型特性對(duì)話(huà)框，進(jìn)入圖像重采樣階段)第七步:圖像重采樣(Resample the Image):圖像重采樣簡(jiǎn)介(Introduction to image Resample )重采樣過(guò)程就是依據(jù)未校正圖像像元值計(jì)

30、算生成一幅校正圖像的過(guò)程，原圖像中所 有柵格數(shù)據(jù)層都將進(jìn)行重采樣。ERDAS IMAGINE提供三種最常用的重采樣方法：(a) Neatest Neighbor:鄰近點(diǎn)插值法，將最鄰近像元值直接賦予輸出像元；(b ) Bilinear Interpolation:雙線(xiàn)性插值法，用雙線(xiàn)性方程和2X2窗口計(jì)算輸出像元值;(c) Cubic Convolution:立方卷積插值法，用立方方程和4X4窗口計(jì)算輸出像元值；圖像重采樣過(guò)程(Process of Image Resample )首先，在Geo Correction Too、對(duì)話(huà)框中選擇Image Resampl。圖標(biāo) 打開(kāi)Image Res

31、ample(圖像重采樣)對(duì)話(huà)框(圖423)然后，在Image Resample對(duì)話(huà)框中，定義重采樣參數(shù)：一輸出圖像文件名(Output File):rectify.img一選擇重采樣方法(Resample Method): Nearest Neighbor一定義輸出圖像范圍(0output Corners): ULX, ULY LRX, LRY一定義輸出像元大小(Output Cell Size): X:30. Y30一設(shè)置輸出統(tǒng)計(jì)中忽略零值:Ignore Zero in Stats.一設(shè)置重新計(jì)算輸出缺省值(Recalculate Output Default): Skip Factor.

32、10一 OK(關(guān)閉Image Resample。對(duì)話(huà)框，啟動(dòng)重采樣進(jìn)程)第八步:保存幾何校正模式(Save Rectification Mode)在Geo Correction Tools對(duì)話(huà)框中點(diǎn)擊Exit按鈕，退出圖像幾何校正過(guò)程，按照系統(tǒng) 提示選擇保存圖像幾何校正模式，并定義模式文件(*.gm )，以便下次直接使用。第九步:檢驗(yàn)校正結(jié)果(Verify Rectification Result)檢驗(yàn)校正結(jié)果的基本方法是:同時(shí)在兩個(gè)視窗中打開(kāi)兩幅圖像，其中一幅是校正以后 的圖像，一幅是當(dāng)時(shí)的參考圖像，通過(guò)視窗地理連接(Geo Link/Unlink)功能及查詢(xún)光標(biāo) (Inquirecurso

33、r)功能進(jìn)行目視定性檢臉，具體過(guò)程如下：打開(kāi)兩個(gè)平鋪視窗(open and Tile two Viewer)視窗菜單條:File 一 OpenRasterOption 一圖像文件ERDAS 圖標(biāo)面板:SessionTile Viewer 一平鋪視窗建立視窗地理連接關(guān)系(Geo Link two Viewer)在Viewer #1中:按住右鍵一快捷菜單一 Geo Link /Unlink在Viewer #2中:點(diǎn)擊左鍵一建立與Viewer #1的連接通過(guò)查詢(xún)光標(biāo)進(jìn)行檢驗(yàn)(Check with Inquire Cursor)在Viewer #1中:按住右鍵一快捷菜單一 Inquire Cursor

34、 一打開(kāi)光標(biāo)查詢(xún)對(duì)話(huà)框在Viewer #1中:移動(dòng)查詢(xún)光標(biāo)，觀測(cè)其在兩屏幕中的位置及匹配程度，并注意光標(biāo) 查詢(xún)對(duì)話(huà)框中數(shù)據(jù)的變化。如果滿(mǎn)意的話(huà)，關(guān)閉光標(biāo)查詢(xún)對(duì)話(huà)框。4、衛(wèi)星圖像拼接 (Satellite Image Mosaic)圖像拼接處理(Mosaic Image)的說(shuō)明：圖像拼接處理是要將具有地理參考的若干相鄰圖像合并成一幅圖像或一組圖像，需 要拼接的輸入圖像必須含有地圖投影信息，或者說(shuō)輸入圖像必須經(jīng)過(guò)幾何校正處理 (Rectified)或進(jìn)行過(guò)校正標(biāo)定(Calibrated)。雖然所有的輸入圖像可以具有不同的投影類(lèi) 型、不同的像元大小，但必須具有相同的波段數(shù)。在進(jìn)行圖像拼接時(shí)，需要確定

35、一幅參 考圖像，參考圖像將作為輸出拼接圖像的基準(zhǔn)，決定拼接圖像的對(duì)比度匹配以及輸出圖 像的地圖投影、像元大小和數(shù)據(jù)類(lèi)型。Mosaic Tool視窗由菜單條(Menu Bar)、工具條(Tool Bar)、圖形窗口(Graphic View)、 和狀態(tài)條(Status Bar)四部分組成。本次實(shí)習(xí) 我們通過(guò)三幅陸地資源衛(wèi)星(wasia_mss.img, wasia2_mss.img, wasia3_tm.img)接處理，介紹衛(wèi)星圖像的拼接過(guò)程。啟動(dòng)圖像拼接工具(Start Mosaic Tools)圖像拼接工具可以通過(guò)下列兩種途徑啟動(dòng)：ERDAS 圖標(biāo)面板菜單條:Main 一 Data Prep

36、aration-Data Preparation 菜單一選擇 Mosaic Images 一打開(kāi) Mosaic Tool 視窗ERDAS圖標(biāo)面板工具條:點(diǎn)擊Data Prep圖標(biāo)一打開(kāi)Data Preparation菜單一選擇 Mosaic Images 一打開(kāi) Mosaic Tool視窗2.2 加載 Mosaic 圖像(Add Images for Mosaic)Mosaic Tool 視窗菜單條:Edit Add Images一打開(kāi) Add Images for Mosaic 對(duì)話(huà)框Mosaic Tool視窗工具條:點(diǎn)擊Add圖標(biāo)白*一打開(kāi) Add Images for Mosaic 對(duì)話(huà)

37、框Add Images for Mosaic對(duì)話(huà)框(選擇文件后)在Add Images for Mosaic對(duì)話(huà)框中，需要設(shè)置以下參數(shù)： 拼接圖像文件(Image FileName): wasial mss.img圖像拼接區(qū)域(Image Area Option)- Compute Active Area (edge)Add(圖像wasial_ mss.img被加載到 Mosaic視窗中) 同樣的過(guò)程加載 wasia2_mss.img 和 wasia3_tm.imgClose(關(guān)閉 Add Images for Mosaic 對(duì)話(huà)框)圖像疊置組合(Images Stacking)Mosaic

38、Too:視窗工具條:點(diǎn)擊Set，Input Mode圖標(biāo)口，進(jìn)入設(shè)置輸入圖像模式的狀態(tài)，Mosaic Tool視窗工具條中會(huì)出現(xiàn)與該模式對(duì)應(yīng)的調(diào)整圖像疊置次序的編輯圖標(biāo)。充分利 用系統(tǒng)所提供的編輯工具，根據(jù)需要進(jìn)行上下層調(diào)整，這些調(diào)整工具包括：Send Image to Top:將選擇圖像置于最上層Send Image Up One:將選擇圖像上移一層Send Image to Bottom:將選擇圖像置于最下層Send Image Down One:將選擇圖像下移一層Reverse Imge Order:將選擇圖像次序顛倒調(diào)整完成后，在Mosaic Tool視窗圖形窗口點(diǎn)擊一下，退出圖像疊置

39、組合狀態(tài)。圖像匹配設(shè)置(Image Matching)Mosaic Tool視窗菜單條:Edit Image Matching -打開(kāi) Matching Options 對(duì)話(huà)框Mosaic Tool，視窗工具條：點(diǎn)擊Set Input Mode。圖標(biāo)口，進(jìn)入設(shè)置輸入圖像模式， 點(diǎn)擊;Image Matchin:圖標(biāo)甌，打開(kāi)Matching Options對(duì)話(huà)框Set Overlap Function 對(duì)話(huà)框Matching Options 對(duì)話(huà)框 在Matching Options對(duì)話(huà)框中做如下設(shè)置： 設(shè)置匹配方法（Matching Method）: Overlap Areas（重疊區(qū)域匹配）

40、 。電保存設(shè)置，關(guān)閉Matching Options對(duì)話(huà)框）Set Overlap Function 對(duì)話(huà)框Mosaic Tool 視窗菜單條:Edit Set Overlap Function 打開(kāi) Set Overlap Function 對(duì)話(huà)框，進(jìn)入設(shè)置圖像關(guān)系模式Mosaic Too;視窗工具條:點(diǎn)擊Set Intersection Mod。圖標(biāo)回，進(jìn)入設(shè)置圖像關(guān)系模式點(diǎn)擊 Overlap Function 圖標(biāo) 1打開(kāi) Set Overlap Function 對(duì)話(huà)框在Set Overlap Function對(duì)話(huà)框中，設(shè)置以下參數(shù)：一設(shè)置相交關(guān)系（Intersection Metho

41、d）- No Cutline Exists （沒(méi)有裁切線(xiàn)）一設(shè)置重疊區(qū)像元灰度計(jì)算（Select Function） : Average （均值）一 Apply（保存設(shè)置）一 Close（關(guān)閉 Matching Options 對(duì)話(huà)框）運(yùn)行 Mosaic 工具（Run the Mosaic）Mosaic Tool 視窗菜單條:Process - Run Mosaic 一打開(kāi) Run Mosaic 對(duì)話(huà)框Run Mosaic對(duì)話(huà)框在Run Mosaic對(duì)話(huà)框中，設(shè)置下列參數(shù)：確定輸出文件名(Output File Name)- wasia_mosaic.img確定輸出圖像區(qū)域(Witch Outputs) : All忽略輸入圖像值(Ignore Input Value): 0輸出圖像背景值(Output Background Value):0忽略輸出統(tǒng)計(jì)值(Stars Ignore Value): 0OK(關(guān)閉Run Mosaic對(duì)話(huà)框，運(yùn)行圖像拼接)I 退出 Mosaic 工具(Exit the Mosaic):Mosaic Tool視窗菜單條:File 一 Close系統(tǒng)提示是否保存Mosaic設(shè)置？ No(關(guān)閉Mosaic Tool視窗，退出Mosaic工具)。將拼接后圖像截圖保存至作業(yè)說(shuō)明5、圖像投影變換(Reproject Ima