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