【正文】 具體步。在ironoxide圖中，地表含鐵礦物多的地方，尤其是開采礦區(qū)、礦渣堆放區(qū)等，ironoxide指數(shù)就高，光譜值就大，在圖上呈亮白色。在NDVI圖中，植被覆蓋多的地方光譜值就大，呈灰白色，相反像水體裸地等無植被覆蓋光譜值就小，呈暗黑色。兩幅同樣行、列數(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ò)展，以突出不同波段見地物光譜的差異，提高對(duì)比度。由于不同地物反射率差值不同，兩波段亮度值相減后，差值大的被突出出來。2圖像運(yùn)算兩幅或多幅單波段影像，完成空間配準(zhǔn)后，通過一系列運(yùn)算，可以實(shí)現(xiàn)圖像增強(qiáng)，達(dá)到提取某些信息或去掉某些不必要信息的目的。在531波段合成圖像上，礦區(qū)的影像特征大致為：顏色呈淡綠色、青偏綠色，較大的礦區(qū)呈不規(guī)則多邊形狀，無植被覆蓋。數(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作假彩色合成后的圖像（見圖36），我們可以發(fā)現(xiàn)：合成之后的圖像可以更好的反映水體污染的情況，清潔水體的影像呈均勻的深藍(lán)藍(lán)黑色，輕污染區(qū)的水體影像藍(lán)中夾雜著不均勻的呈條狀、塊狀分布的淡藍(lán)、淡青色調(diào)，重污染區(qū)的水體影像呈現(xiàn)淡藍(lán)色青色。Landsat TM（ETM）的7個(gè)波段中，TMETM+2（）適用于水體污染的研究。由于原色的選擇與原來遙感波段所代表的真實(shí)顏色不同，因此生成的合成色不是地物真實(shí)的顏色，因此這種合成叫做假彩色合成。 Haze Reduction 對(duì)話框（1）ERDAS 圖標(biāo)面板菜單條.: 單擊MainImage InterpreterRadiometric Enhancement打開 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í)行去霾處理。對(duì)于全色圖像，該方法采用點(diǎn)擴(kuò)展卷積反轉(zhuǎn)（Inverse Point Spread Convolution）進(jìn)行處理，根據(jù)情況選擇55或33的卷積運(yùn)算于分別用于高頻模糊度（Highhaze）或低頻模糊度（Lowhaze）的去除。練習(xí)1 去霾處理（Haze Reduction）去霾處理的目的是降低多波段圖像（Landsat TM）或全色圖像的模糊度（霾）。8實(shí)用分析功能（Utilities）實(shí)用分析功能包括了基本的圖像處理操作。6地形分析（Topographic Analysis）地形分析功能主要是在點(diǎn)、線、面高程基礎(chǔ)上，對(duì)多種地形因素進(jìn)行分析，井對(duì)圖像進(jìn)行地形校正。3光譜增強(qiáng)（Spectral Enhancement）光譜增強(qiáng)處理是基于多波段數(shù)據(jù)對(duì)每個(gè)像元的灰度值進(jìn)行變換，達(dá)到圖像增強(qiáng)的目的。1空間增強(qiáng)（Spatial Enhancement）空問增強(qiáng)技術(shù)是利用像元自身及其周圍像元的灰度值進(jìn)行運(yùn)算，達(dá)到增強(qiáng)整個(gè)圖像之目的。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)沒有達(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í)行過程中都需要您通過各種按鍵或?qū)υ捒蚨x參數(shù)，多數(shù)解譯功能都借助模型生成器（Model Maker）建立了圖形模型算法，很容易調(diào)用或編輯。⑧ 選擇重采樣方法（Resample Method）: Nearest Neighbor。⑥ 確定輸出統(tǒng)計(jì)缺省零值:Ignore Zero in Stats。④ 定義投影參數(shù)（Projection）: UTM Zone 50 （Range 114E120E）。② 定義輸出圖像文件（Output File）: 。數(shù)據(jù)文件存放在%sampledate/project下（1）啟動(dòng)投影變換（Start Reproject）圖像投影變換功能既可以數(shù)據(jù)預(yù)處理模塊中啟動(dòng)，也可以在圖像解譯模塊中啟動(dòng)。圖標(biāo)，進(jìn)入設(shè)置圖像關(guān)系模式點(diǎn)擊Overlap Function圖標(biāo)打開Set Overlap Function對(duì)話框（） Run Mosaic對(duì)話框在Set Overlap Function對(duì)話框中，設(shè)置以下參數(shù):①設(shè)置相交關(guān)系（Intersection Method） No Cutline Exists（沒有裁切線）；②設(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打開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工具）。圖標(biāo)，進(jìn)入設(shè)置輸入圖像模式,點(diǎn)擊。調(diào)整完成后，在Mosaic Tool視窗圖形窗口點(diǎn)擊一下，退出圖像疊置組合狀態(tài)。（3）圖像疊置組合（Images Stacking）Mosaic 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)。（1） 啟動(dòng)圖像拼接工具（Start Mosaic Tools）圖像拼接工具可以通過下列兩種途徑啟動(dòng):ERDAS圖標(biāo)面板菜單條:MainData PreparationData Preparation菜單（）選擇Mosaic Images打開Mosaic Tool視窗（）ERDAS圖標(biāo)面板工具條:點(diǎn)擊Data Prep圖標(biāo)打開Data Preparation菜單（）選擇Mosaic Images打開Mosaic Tool視窗（）（2）加載Mosaic圖像（Add Images for Mosaic） Mosaic Tool視窗菜單條:EditAdd Images打開Add Images for Mosaic對(duì)話框（）Mosaic Tool視窗工具條:點(diǎn)擊Add圖標(biāo)，打開Add Images for Mosaic對(duì)話框（）在Add Images for Mosaic對(duì)話框中，需要設(shè)置以下參數(shù):①拼接圖像文件（Image FileName）: wasial ；②圖像拼接區(qū)域（Image Area Option）Compute Active Area （edge）；③Add（圖像wasial_ ）。本次實(shí)習(xí)我們通過三幅陸地資源衛(wèi)星（, , ）拼接處理，介紹衛(wèi)星圖像的拼接過程。在進(jìn)行圖像拼接時(shí)，需要確定一幅參考圖像，參考圖像將作為輸出拼接圖像的基準(zhǔn)，決定拼接圖像的對(duì)比度匹配以及輸出圖像的地圖投影、像元大小和數(shù)據(jù)類型。練習(xí)3 衛(wèi)星圖像拼接（Satellite Image Mosaic）圖像拼接處理（Mosaic Image）的說明：圖像拼接處理是要將具有地理參考的若干相鄰圖像合并成一幅圖像或一組圖像，需要拼接的輸入圖像必須含有地圖投影信息，或者說輸入圖像必須經(jīng)過幾何校正處理（Rectified）或進(jìn)行過校正標(biāo)定（Calibrated）。（13）檢驗(yàn)校正結(jié)果（Verify Rectification Result）檢驗(yàn)校正結(jié)果的基本方法是:同時(shí)在兩個(gè)視窗中打開兩幅圖像，其中一幅是校正以后的圖像，一幅是當(dāng)時(shí)的參考圖像，通過視窗地理連接（Geo Link/Unlink）功能及查詢光標(biāo)（Inquirecursor）功能進(jìn)行目視定性檢臉，具體過程如下:（a）打開兩個(gè)平鋪視窗（open and Tile two Viewer）,視窗菜單條:FileOpenRasterOption圖像文件,ERDAS圖標(biāo)面板:SessionTile Viewer平鋪視窗（b）建立視窗地理連接關(guān)系（Geo Link two Viewer）,在Viewer 1中:按住右鍵快捷菜單Geo Link /Unlink,在Viewer 2中:點(diǎn)擊左鍵一建立與Viewer 1的連接（c）通過查詢光標(biāo)進(jìn)行檢驗(yàn)（Check with Inquire Cursor）,在Viewer 1中:按住右鍵快捷菜單Inquire Cursor打開光標(biāo)查詢對(duì)話框,在Viewer 1中:移動(dòng)查詢光標(biāo)，觀測(cè)其在兩屏幕中的位置及匹配程度，并注意光標(biāo)查詢對(duì)話框中數(shù)據(jù)的變化。⑤設(shè)置輸出統(tǒng)計(jì)中忽略零值:Ignore Zero in Stats。③定義輸出圖像范圍（0output Corners）: ULX, ULY, LRX, LRY。圖標(biāo)（b）打開Image Resample（圖像重采樣）對(duì)話框（c）然后，在Image Resample對(duì)話框中，定義重采樣參數(shù):①輸出圖像文件名（Output File）:。（3） Cubic Convolution:立方卷積插值法，用立方方程和4X4窗口計(jì)算輸出像元值。ERDAS IMAGINE提供三種最常用的重采樣方法: （1） Neatest Neighbor:鄰近點(diǎn)插值法，將最鄰近像元值直接賦予輸出像元。否則，若控制點(diǎn)定義不精確，檢查點(diǎn)就無法匹配，誤差會(huì)超標(biāo)。（f）計(jì)算檢查點(diǎn)誤差:在GCP Tool工具條中點(diǎn)擊Compute Error圖標(biāo)，檢查點(diǎn)的誤差就會(huì)顯示在GCP Tool的上方，只有所有檢查點(diǎn)的誤差均小于一個(gè)像元（Pixel），才能繼續(xù)進(jìn)行合理的重采樣。依然在GCP Tool對(duì)話框狀態(tài)下:（a）在GCP Tool菜單條中確定GCP類型:EditSet Point TypeCheck（b）在GCP Tool菜單條中確定GCP匹配參數(shù)（Matching Parameter）: Edit Point Matching打開GCP Matching對(duì)話框（略） GCP Matching對(duì)話框 在GCP Matching對(duì)話框中，需要定義下列參數(shù):①匹配參數(shù)（Matching Parameters）；②最大搜索半徑（Max. Search Radius）為3；③搜索窗口大小（Search Window Size）: X:5 Y:5。）（9）采集地面檢查點(diǎn) Tool 對(duì)話框與GCP數(shù)據(jù)表以上所采集的GCP的類型為Control Point（控制點(diǎn)），用于控制計(jì)算，建立轉(zhuǎn)換模型及多項(xiàng)式方程。（j）不斷重復(fù) （a） （i），采集若干GCP,直到滿足所選定的幾何校正模型為止；而后，每采集一個(gè)Input GCP, 系統(tǒng)自動(dòng)產(chǎn)生 一個(gè)Ref. GCP, 可以逐步優(yōu)化校正模型。（b）在GCP數(shù)據(jù)表中將輸入GCP的顏色（Color）設(shè)置為比較明顯的黃色，（c）在Viewer 1中移動(dòng)關(guān)聯(lián)方框位置，尋找明顯的地物特征點(diǎn)，作為輸入GCP:（d）在GCP工具對(duì)話框中點(diǎn)擊Create GCP圖標(biāo)，并在Viewer 3中點(diǎn)擊左鍵定點(diǎn)，包括其編號(hào)、標(biāo)識(shí)碼、X坐標(biāo)、Y坐標(biāo):（e）在GCP工具對(duì)話框中點(diǎn)擊Select GCP圖標(biāo)，重新進(jìn)入GCP選擇狀態(tài)。（e）參考GCP也可以類似地保存在參考圖像中（Save Reference）或GCP文件中（（Save Reference As），便于以后調(diào)用。（d）所有的輸入GCP都可以直接保存在圖像文件中（Save Ii）put），也可以保