【正文】 Overlap Function對話框（）Mosaic Tool視窗工具條:點擊Set Intersection Mod。圖標，進入設(shè)置圖像關(guān)系模式點擊Overlap Function圖標打開Set Overlap Function對話框（） Run Mosaic對話框在Set Overlap Function對話框中，設(shè)置以下參數(shù):①設(shè)置相交關(guān)系（Intersection Method） No Cutline Exists（沒有裁切線）；②設(shè)置重疊區(qū)像元灰度計算（Select Function） : Average（均值）；③Apply（保存設(shè)置）Close（關(guān)閉Matching Options對話框）（6）運行Mosaic工具（Run the Mosaic）Mosaic Tool視窗菜單條:Process Run Mosaic打開Run Mosaic對話框（）在Run Mosaic對話框中，設(shè)置下列參數(shù):①確定輸出文件名（Output File Name） ；②確定輸出圖像區(qū)域（Witch Outputs） : All；③忽略輸入圖像值（Ignore Input Value）: 0；④輸出圖像背景值（Output Background Value）:0；⑤忽略輸出統(tǒng)計值（（Stars Ignore Value）: 0；⑥OK（關(guān)閉Run Mosaic對話框，運行圖像拼接）（7）退出Mosaic工具（Exit the Mosaic）:Mosaic Tool視窗菜單條:FileClose，系統(tǒng)將會提示是否保存Mosaic設(shè)置（8）No（關(guān)閉Mosaic Tool視窗，退出Mosaic工具）。練習4 圖像投影變換（Reproject Images）圖像投影變換的目的在于將圖像文件從一種地圖投影類型轉(zhuǎn)換到另一種投影類型，這種轉(zhuǎn)換可以對單幅圖像進行，也可以通過批處理向?qū)В˙atch Wizard），這種直接的投影變換可以避免多項式近似值的擬合，對于大范圍的地理參考是非常有意義的。數(shù)據(jù)文件存放在%sampledate/project下（1）啟動投影變換（Start Reproject）圖像投影變換功能既可以數(shù)據(jù)預(yù)處理模塊中啟動，也可以在圖像解譯模塊中啟動。（a）在數(shù)據(jù)預(yù)處理模塊（Data Preparation）中可以通過兩種途徑啟動:ERDAS圖標面板菜單條:MainData PreparationReproject Images，Reproject Images對話框（） Images對話框（設(shè)置參數(shù)之后）ERDAS圖標面板工具條:點擊Daterep圖標Reproject ImagesRePject Images對話框（）（b）在圖像解譯模塊（Image Interpreter）中也可以通過兩種途徑啟動:ERDAS圖標面板菜單條:MainImage InterpreterUtilitiesReproject ImagesReproject Images對話框（）ERDAS圖標面板工具條:點擊DataPrep圖標UtilitiesReproject ImagesReproject Images對話框（）（2）投影變換操作（Reproject Operation）在Reproject Images對話框中必須設(shè)置下列參數(shù)，方可進行投影變換:① 確定輸入圖像文件（Input File）: 。② 定義輸出圖像文件（Output File）: 。③ 定義投影類型（Categories）: UTM Clark 1866 North。④ 定義投影參數(shù)（Projection）: UTM Zone 50 （Range 114E120E）。⑤ 定義輸出圖像單位（Units）: Meters/ Feet /Degrees。⑥ 確定輸出統(tǒng)計缺省零值:Ignore Zero in Stats。⑦ 定義輸出像元大?。∣utput Cell Sizes）: X: Y: 。⑧ 選擇重采樣方法（Resample Method）: Nearest Neighbor。⑨ 定義轉(zhuǎn)換方法:Rigorous Transformation（嚴格按照投影數(shù)學模型進行變換）。Polynomial Approximation（應(yīng)用多項近似擬合實現(xiàn)變換）如果選擇Polynomial Approximation轉(zhuǎn)換方法，還需設(shè)置下列參數(shù):多項式最大次方（Maximum Poly Order）: 3；定義像元誤差（Tolerance Pixels）: 1；如果在設(shè)置的最大次方內(nèi)沒有達到像元誤差要求，則按照下列設(shè)置執(zhí)行：If Tolerance Exceeded: Continue Approximate（依然應(yīng)用多項式模型轉(zhuǎn)換）；Rigorous Transformation（嚴格按照投影模型轉(zhuǎn)換）（3）OK（關(guān)閉Reproject Images對話框，執(zhí)行投影變換）實習3　圖像解譯實習內(nèi)容: 去霾處理光譜增強圖像空間增強課時安排：4課時圖像解譯功能簡介（Introduction of Image Interpreter）ERADS IMAGINE的圖像解譯器（Image Interpreter）包含了50多個用于遙感圖像處理的功能模塊，這些功能模塊在執(zhí)行過程中都需要您通過各種按鍵或?qū)υ捒蚨x參數(shù)，多數(shù)解譯功能都借助模型生成器（Model Maker）建立了圖形模型算法，很容易調(diào)用或編輯。圖像解譯器又稱Image Interpreter或Interpreter，可以通過兩種途徑啟動:ERDAS圖標面板菜單條:MainImage IrnerpreterImage Interpreter菜單（）ERDAS圖標面板工具條:點擊Interpreter圖標Image Interpreter菜單（）Image Interpreter面板，ERDAS圖像解譯模塊包含了8個方面的功能，依次是遙感圖像的空間增強（Spatial Enhancement）、輻射增強（Radiometric Enhancemen）、光譜增強（SpectralEnhancement）、高光譜工具（Hyper Spectral Tools）、傅立葉變換（Fourier Analysis）、地形分析（Topographic Analysis）、地理信息系統(tǒng)分析（GIS Analysis）、以及其它實用功能（Utilities），每一項功能菜單中又包含若干具體的遙感圖像處理功能。1空間增強（Spatial Enhancement）空問增強技術(shù)是利用像元自身及其周圍像元的灰度值進行運算，達到增強整個圖像之目的。2輻射增強（Radiometric Enhancement）輻射增強處理是對單個像元的灰度恒進行變換運到圖像增強的目的。3光譜增強（Spectral Enhancement）光譜增強處理是基于多波段數(shù)據(jù)對每個像元的灰度值進行變換，達到圖像增強的目的。4高光譜工具（ Hyper spectral Tools）高光譜工具是通過補償大氣對光譜的混淆來增強圖像．5傅立葉變換（ Fourier Analysis）傅立葉變換是首先把遙感圖像從空間域轉(zhuǎn)換到頻率域，然后在頻率域上對圖像進行濾波處理，減少或消除周期性噪聲，再把圖像從頻率域轉(zhuǎn)換到空間域，達到增強圖像的目的。6地形分析（Topographic Analysis）地形分析功能主要是在點、線、面高程基礎(chǔ)上，對多種地形因素進行分析，井對圖像進行地形校正。7地理信息系統(tǒng)分析（GIS Analysis）地理信息系統(tǒng)分析功能主要是對圖像進行各種空間分析，涉及像元之間或?qū)ｎ}分類之間的空間關(guān)系處理，使處理后的圖像更好地表達主要的專題信息。8實用分析功能（Utilities）實用分析功能包括了基本的圖像處理操作。其中常用的圖象掩膜和投影變換已經(jīng)練習過。練習1 去霾處理（Haze Reduction）去霾處理的目的是降低多波段圖像（Landsat TM）或全色圖像的模糊度（霾）。對于多波段圖像（Landsat TM），該方法實質(zhì)上是基于纓穗變換方法（Tasseled Cap Transformation），首先對圖像逆行主成份變換，找出與模糊度相關(guān)的成份井剔除，然后再進行主成份逆變換回到RGB彩色空間，達到去霾之目的。對于全色圖像，該方法采用點擴展卷積反轉(zhuǎn)（Inverse Point Spread Convolution）進行處理，根據(jù)情況選擇55或33的卷積運算于分別用于高頻模糊度（Highhaze）或低頻模糊度（Lowhaze）的去除。數(shù)據(jù)文件存放在%sampledate/ haze下。 Haze Reduction 對話框（1）ERDAS 圖標面板菜單條.: 單擊MainImage InterpreterRadiometric Enhancement打開 Haze Reduction 對話框（）,在對話框中，需要設(shè)置下列參數(shù)：①確定輸入文件（Input File）：；②定義輸出文件（Ouput File）：；③定義坐標類型為（Coordinate Type）：Map；④處理范圍確定（Subset Definition）：在ULX／Y、LRX／Y微調(diào)框中輸入需要的數(shù)值（缺省狀態(tài)為整個圖像范圍，可以應(yīng)用InquireBox定義子區(qū)）；⑤處理方法（Method）：選擇Landsat 5 TM（或Landsat 4 TM）（2）單擊OK按鈕，關(guān)閉 Haze Reduction對話框，執(zhí)行去霾處理。去霾處理（HazeReduction）： Haze處理前的TM圖像練習2 光譜增強1彩色變換（1）突出顯示水體信息根據(jù)加色法的原理，選擇遙感影像的某三個波段，分別賦予紅、綠、藍三種原色，就可以合成彩色影像。由于原色的選擇與原來遙感波段所代表的真實顏色不同，因此生成的合成色不是地物真實的顏色，因此這種合成叫做假彩色合成。多波段影像合成時，方案的選擇十分重要，它決定了彩色圖像能否顯示較豐富的地物信息或突出某一方面的信息。Landsat TM（ETM）的7個波段中，TMETM+2（）適用于水體污染的研究。不少研究表明，水體的蘭譜段和綠譜段的輻射比值，是溶解有機質(zhì)和浮游生物出現(xiàn)情況的一種度量。數(shù)據(jù)文件存放在%sampledate/ spectral下 Set Layer Combinations for 窗口用2作假彩色合成，具體操作如下：，在Viewer窗口中單擊Raster工具條，單擊Band Combinations,出現(xiàn)Set Layer Combinations for ， 1986年TM432合成圖在對話框中分別選擇Red為4，Green為3，Blue為2，單擊Apply和ok后，退出對話框，審視2作假彩色合成后的圖像（見圖36），我們可以發(fā)現(xiàn)：合成之后的圖像可以更好的反映水體污染的情況，清潔水體的影像呈均勻的深藍藍黑色，輕污染區(qū)的水體影像藍中夾雜著不均勻的呈條狀、塊狀分布的淡藍、淡青色調(diào)，重污染區(qū)的水體影像呈現(xiàn)淡藍色青色。（2）突出礦點和居民點信息銅綠山銅鐵礦 1986年TM743合成圖銅山口鎮(zhèn)TMETM+1（）能反映巖礦石中鐵離子疊加吸收譜帶，為褐鐵礦、鐵帽特征識別譜帶，并且礦物在紅光波段TMETM+3（）的反射值較大，RGB為1或1波段組合（見圖37）有利于突出礦區(qū)信息。在531波段合成圖像上，礦區(qū)的影像特征大致為：顏色呈淡綠色、青偏綠色，較大的礦區(qū)呈不規(guī)則多邊形狀，無植被覆蓋。而據(jù)Barry（1983）研究，用3或4作假彩色合成有利于居民地的提取，743波段合成圖像上（見圖38），居民地影像顏色呈紫色。2圖像運算兩幅或多幅單波段影像，完成空間配準后，通過一系列運算，可以實現(xiàn)圖像增強，達到提取某些信息或去掉某些不必要信息的目的。兩幅同樣行、列數(shù)的圖像，對應(yīng)像元的亮度值相減就是差值計算，即：fD（x，y）= f1（x，y）f2（x，y）差值計算應(yīng)用于兩個波段時，相減后的值反映了同一地物光譜反射率之間的差。由于不同地物反射率差值不同，兩波段亮度值相減后，差值大的被突出出來。例如，當紅外波段減紅波段時，植被的反射率差異很大，相減后的差值就大，而土壤和水在這兩個波段反射率差值就很小，因此相減后的圖像可以把植被信息突出出來。兩幅同樣行、列數(shù)的圖像，對應(yīng)像元的亮度值相除（除數(shù)不為0）就是比值計算，即：fR（x，y）= f1（x，y）/f2（x，y）圖311 ETM圖像重采樣得到和SPOT5全色波段一樣的像元對話框 1986年NDVI圖 1986年Iron Oxide圖銅綠山銅鐵礦比值運算可以檢測波段的斜率信息并加以擴展，以突出不同波段見地物光譜的差異，提高對比度。該運算常用于突出遙感影像中的植被特征、提取植被類型或估算植被生物量，這種算法稱為植被指數(shù)（VI，Vegetation Index），如歸一化差異植被指數(shù)NDVI（見圖39），（TM4TM3）/（TM4+TM3），是最廣泛使用的一種VI。在NDVI圖中，植被覆蓋多的地方光譜值就大，呈灰白色，相反像水體裸地等無植被覆蓋光譜值就小，呈暗黑色。類似，鐵礦指數(shù)（即ironoxide指數(shù)，實質(zhì)為紅光波段（band3）與藍光波段（band1）反射值的比值，見圖310），體現(xiàn)了含鐵氧化物的礦石band3/band1比值較大的特點。在ironoxide圖中，地表含鐵礦物多的地方，尤其是開采礦區(qū)、礦渣堆放區(qū)等，ironoxide指數(shù)就高，光譜值就大，在圖上呈亮白色。接著通過指數(shù)運算和波段組合得到礦區(qū)指數(shù)圖（R（5/7）、G（5/4）、B（3/1）），由于TM（ETM）的波段3/1和5/4能體現(xiàn)含鐵和鐵的氧化物，而5/7表示含粘土礦物的特點，可以大致提取礦區(qū)范圍（）