中國礦業(yè)大學(xué)《遙感原理與應(yīng)用》ENVI實(shí)習(xí)報(bào)告_第1頁
中國礦業(yè)大學(xué)《遙感原理與應(yīng)用》ENVI實(shí)習(xí)報(bào)告_第2頁
中國礦業(yè)大學(xué)《遙感原理與應(yīng)用》ENVI實(shí)習(xí)報(bào)告_第3頁
中國礦業(yè)大學(xué)《遙感原理與應(yīng)用》ENVI實(shí)習(xí)報(bào)告_第4頁
中國礦業(yè)大學(xué)《遙感原理與應(yīng)用》ENVI實(shí)習(xí)報(bào)告_第5頁
已閱讀5頁,還剩51頁未讀, 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡介

1、遙感原理與應(yīng)用實(shí)習(xí)報(bào)告學(xué)院:班級(jí):學(xué)號(hào):姓名:中國礦業(yè)大學(xué)環(huán)境與測繪學(xué)院遙感實(shí)習(xí)2010年12月26日遙感原理與應(yīng)用- -目錄1、實(shí)驗(yàn)一高光譜數(shù)據(jù)分析.22、實(shí)驗(yàn)二影像鑲嵌.153、實(shí)驗(yàn)三影像配準(zhǔn)264、實(shí)驗(yàn)四非監(jiān)督分類385、實(shí)驗(yàn)五監(jiān)督分類416、實(shí)習(xí)總結(jié).49遙感原理與應(yīng)用- -實(shí)驗(yàn)一高光譜數(shù)據(jù)分析一、實(shí)驗(yàn)?zāi)康谋緦n}旨在向用戶介紹波譜庫的概念,并描述如何從感興趣區(qū)中提取波譜信息,然后還將進(jìn)行彩色合成,并使用二維散點(diǎn)圖進(jìn)行簡單的分類。讓學(xué)生學(xué)會(huì)如何使用ENVI先進(jìn)的高光譜工具對多光譜數(shù)據(jù)進(jìn)行分析。更好地理解高光譜處理的概念及其工具。本專題將從特定礦物質(zhì)的感興趣區(qū)中提取其波譜曲線,并與波譜庫中

2、的波譜曲線進(jìn)行比較,找出顯示波譜信息的最佳RGB彩色組合。使用二維散點(diǎn)圖定位獨(dú)特的像元,探究其數(shù)據(jù)的分布特點(diǎn),然后進(jìn)行簡單的分類。二、實(shí)驗(yàn)數(shù)據(jù)介紹文件描述zhumuqian95_Cuprite地區(qū)ATREM校正后的反射率數(shù)據(jù),50個(gè)波段(整型)zhumuqian95at.hdrENVI相應(yīng)的頭文件Jpll.sliENVI格式的JPL波譜庫Jpll.hdrENVI相應(yīng)的頭文件zhumuqianmin.sliENVI格式的USGS波譜庫zhumuqian_min.hdrENVI相應(yīng)的頭文件zhumuqian95av.roi保存的感興趣區(qū)文件數(shù)據(jù):1995年的航空可見光/紅外成像光譜儀(

3、AirborneVisible/InfraredImagingSpectrometer,簡稱為AVIRIS)所采集的表觀反射率數(shù)據(jù),該數(shù)據(jù)是美國內(nèi)華達(dá)州(Nevada)Cuprite地區(qū)的表觀反射率數(shù)據(jù),它使用ATREM大氣糾正建模軟件進(jìn)行了校正。這個(gè)數(shù)據(jù)子集共包含50個(gè)波段,波譜分辨率近似為10nm寬,其波長范圍為1.992.48umo三、實(shí)驗(yàn)內(nèi)容本部分將介紹以下內(nèi)容:波譜庫操作、瀏覽和提取影像反射率波譜、ENVI中感興趣區(qū)(ROI)的定義及進(jìn)行彩色合成影像的選取,其目的是為了鑒別波譜類型。1加載AVIRIS影像數(shù)據(jù)并顯示灰階影像在ENVI主菜單中,選擇FilefOpenImageFile,

4、然后選擇進(jìn)入envi47/data目錄。選擇zhumuqian95_at.ir文件作為輸入文件名,點(diǎn)擊Open彈出可用波段列表,它將列出50個(gè)波段的名字。在可用波段列表對話框中,選fiBand193(2.2008um,點(diǎn)擊GrayScale單選按鈕,然后點(diǎn)擊LoadBand將灰度影像加載到顯示窗口中。遙感原理與應(yīng)用- -#1Band193;zhiimuQin95effJnirleDverlavEnhancelods時(shí)nw從主影像窗口菜單中選擇ToolsProfilesZProfile(Spectrum),提取表觀反射率波譜曲線。2瀏覽影像波譜并同波譜庫進(jìn)行比較在影像上移動(dòng)縮放指示矩形框,同時(shí)查

5、看#1SpectralProfile窗口中的波譜曲線,瀏覽整個(gè)影像的表觀反射率波譜曲線。在主影像窗口中,使用鼠標(biāo)左鍵點(diǎn)擊并拖動(dòng)縮放指示矩形框或者直接點(diǎn)擊鼠標(biāo)左鍵,將縮放指示矩形框移動(dòng)到以所選像素點(diǎn)為中心的區(qū)域中。將從影像中獲取的表觀反射率波譜曲線同所選波譜庫中的波譜曲線進(jìn)行比較ENVI提供了幾個(gè)不同的波譜庫,根據(jù)本專題的目的,我們將會(huì)使用JPL波譜庫(Groves等,1992)以及USGS波譜庫(Clarke等,1993)。從ENVI主菜單中選擇SpectralSpectralLibrariesSpectralLibraryViewer。在SpectralLibraryInputFile對話框

6、中,點(diǎn)擊OpenFile按鈕,從spec_lib/jpl_lib子目錄中,選擇jpll.sli波譜庫文件,點(diǎn)擊OK。選擇SelectInputFile區(qū)域中的jpll.sli,點(diǎn)擊OK。在SpectralLibraryViewer對話框中,選擇OptionsfEdit(x,y)ScaleFactors,并在YDataMultiplier文本框中,輸入值1.000,以匹配影像表觀反射率范圍(1T000),點(diǎn)擊OK。在SpectralLibraryViewer對話框中,選擇下列波譜名稱,繪制它們的波譜曲線:ALUNITES0-4ABUDDINGTONITEFELDSTST1ACALCITEC-3D

7、KAOLINITEWELLORDEREDPS-1A得到如下波譜曲線繪制圖:從繪制(plot)窗口菜單中,選擇EditPlotParameters,自定義波譜曲線的繪制圖。在PlotParameters對話框中,按下面的步驟進(jìn)行:將Charsize減少為0.50。選擇X-Axis單選按鈕,然后將Range調(diào)整為1.90到2.45。還是選中X-Axis單選按鈕,點(diǎn)擊Left/RightMargins的箭頭增量按鈕,直到達(dá)到所需X方向的頁邊距。選擇Y-Axis單選按鈕,將AxisTitle改為“Reflectance”。還是選中Y-Axis單選按鈕,點(diǎn)擊Top/BottomMargins的箭頭增量按

8、鈕,直到所需的Y方向的頁邊距。點(diǎn)擊Apply,然后再點(diǎn)擊Cancel。要顯示波譜名稱的圖例,可以在繪制窗口中點(diǎn)擊鼠標(biāo)右鍵,從彈出的快捷菜單中選擇PlotKey。將繪制窗口拖動(dòng)到所需的大小,以容納下波譜名稱。在繪制窗口中,選擇OptionsStackPlots,分別查看繪制的波譜曲線。繪制的遙感原理與應(yīng)用- -波譜曲線如下圖所示:卷SpectralUbraiyPlotsHitEditCptionsPlot/LnctionHdp1)52上JWavelength一-直匚口D.IO亠L(fēng)ln土口遙感原理與應(yīng)用- -在#1SpectralProfile繪圖窗口中點(diǎn)擊鼠標(biāo)右鍵,從彈出的快捷菜單中選擇Plot

9、Key,這將顯示出波譜曲線的圖例,該圖例指出了所選像素的X和Y的像素坐標(biāo)。在#1SpectralProfile繪圖窗口中,選擇OptionsNewWindow:Blank,打開一個(gè)新的繪圖窗口。然后重新放置#1SpectralProfile繪圖窗口和新繪圖窗口的位置,使能夠同時(shí)看到這兩個(gè)繪圖窗口。從主影像窗口菜單中,選擇ToolsPixelLocator。使用PixelLocator對話框,定位到下列各點(diǎn)精確的像素位置上:位置點(diǎn)名稱列(帶偏移)行(帶偏移)StonewallPlaya590570VarnishedTuff435555SilicaCap494514OpaliteZonewithA

10、lunite531541StronglyArgillizedZonewithKaolinite502589BuddingtoniteZone448505Calcite260613在PixelLocator對話框中,輸入像素的坐標(biāo),列(sample)590、行(line)570,使縮放指示矩形框移動(dòng)到以這個(gè)像素為中心的影像地區(qū),即StonewallPlaya地區(qū),然后點(diǎn)擊Apply,將矩形框移動(dòng)到這個(gè)位置上。同時(shí)#1SpectralProfile繪圖窗口將更新顯示所選點(diǎn)的波譜曲線,其所對應(yīng)的圖例為:“X:590Y:570”。在新的繪圖窗口中,點(diǎn)擊鼠標(biāo)右鍵,在彈出的快捷菜單中,選擇PlotKey,

11、打開顯示了X和Y坐標(biāo)位置的圖例。在#1SpectralProfile繪圖窗口中,使用鼠標(biāo)左鍵,點(diǎn)擊并按住圖例“X:590Y:570”,將這個(gè)波譜曲線圖例拖到新的繪圖窗口中。對上表所列的每一個(gè)像素點(diǎn)重復(fù)上面的步驟,直到新的繪圖窗口包含了所有的7種波譜曲線。在新的繪圖窗口中,選擇OptionsStackPlots。新的繪圖窗口如下圖所示:(SinviPlotwindow曹EpeHralLibrar/PloHIFitEdrtOptirjisPbt_Fun-fcn比IpSpecimlLibrarvFlat芒_:匸口p-!-器七sMLI5雖蛋FieEdkDptioraPk3t_RjRmHelpWowel

12、sn勺th將這些波譜曲線同已經(jīng)提取的波譜庫中的波譜曲線進(jìn)行目視比較。注意到,實(shí)驗(yàn)室所測的波譜曲線同某些影像表觀反射率波譜曲線之間存在著相似的形狀及吸收特性?;谶@些相似的特性,我們可以得到這樣的結(jié)論:與alunite、buddingtonite、calcite和kaolinite實(shí)驗(yàn)室所測波譜曲線相似的影像波譜曲線是由上面相應(yīng)的礦物質(zhì)所構(gòu)成的。從波譜庫波譜曲線的繪圖窗口中,將波譜曲線拖動(dòng)到#1SpectralProfile繪圖窗口中,進(jìn)行直接比較。3鑒別波譜曲線使用SpectralAnalystTM來鑒別波譜曲線:ENVI提供了一個(gè)波譜匹配工具,它根據(jù)波譜庫中的波譜曲線對影像中的波譜曲線進(jìn)行評

13、分。波譜分析使用多種方法產(chǎn)生一個(gè)在0到1之間的分?jǐn)?shù)值,其中分?jǐn)?shù)值1相當(dāng)于完全匹配。從ENVI主菜單中,選擇SpectralfSpectralAnalyst。點(diǎn)擊SpectralAnalystInputSpectralLibrary對話框底部的OpenSpecLib按鈕。選擇進(jìn)入usgs_min波譜庫目錄,選擇usgs_min.sli波譜庫文件,點(diǎn)擊Open。usgs_min.sli文件出現(xiàn)在SpectralAnalystInputSpectralLibrary對話框中,選中該文件,點(diǎn)擊OK。在EditIdentifyMethodsWeighting對話框中,點(diǎn)擊OK。從主影像窗口菜單中,選擇T

14、oolsProfilesZProfile(Spectrum)。然后在#1SpectralProfile繪圖窗口中,點(diǎn)擊鼠標(biāo)右鍵,從彈出的快捷菜單中,選擇PlotKey,顯示波譜曲線名稱的圖例。從主影像窗口菜單中,選擇ToolsfPixelLocator。在PixelLocator對話框中,輸入像素點(diǎn)的坐標(biāo)列502,行589,點(diǎn)擊Apply。遙感原理與應(yīng)用- -在SpectralAnalyst對話框中,選擇OptionsfEditMethodWeights。在EditIdentifyMethodsWeighting對話框中,為每一個(gè)Weight文本框輸入值0.33,然后點(diǎn)擊OK。在Spectra

15、lAnalyst對話框中,點(diǎn)擊Apply。如果在#1SpectralProfile繪圖窗口中顯示了多條波譜曲線,那么將會(huì)出現(xiàn)一個(gè)波譜曲線列表。如果出現(xiàn)了該波譜曲線列表,那么就選擇像素(列502,行589)所對應(yīng)的那條波譜曲線。SpectralAnalyst對話框?qū)⑷缦聢D所示:SpectralAnalystOptionsUnknox-m:X:502V:5S9_Libra.rySpectrumScoreSAMkaulini2.spckaolini6.spckaoliniS.:m口匚:kaolini?.spchalloys3.spc口丄i口i日.spcka.ulini3.spcka.olini4.s

16、pcha.lloys2.spckaulinil.spcKaiuliniUPlaoliniUKaoliniUKaoliniUHaHoysUPLs.olini0KaoliniUKaoliniUHa.HoysUKaoliniUTJl11仃.883:0.304.882:0.S04800-0399.000:0.00.376:0.323.875.0.697.871:0.308.3571:0.393.356:0.312.055:0.366口匚c1-/naAm-UJ-I-.t_Mu;-_.t.t_1i11oonl-A-u.b4-UM-Qt-i_H-nonono7-7-7-7-7-7-7_u_u-u-u_u_u

17、_u_u_u_unrrrSpectralAnalyst對話框,顯示出了礦物質(zhì)kaolinite的波譜曲線與像素(列502,行589)的波譜曲線最匹配(左圖)SpectralAnalyst將根據(jù)波譜庫中的波譜曲線對未知地物的波譜曲線進(jìn)行評分。上圖顯示了對像素列502行589)的波譜曲線進(jìn)行鑒別的結(jié)果。注意到列表的第一行顯示kaolinit的波譜曲線評分最高這個(gè)相對較高的分?jǐn)?shù)值表明了該像素對應(yīng)的地物與kaolinit最相似。用鼠標(biāo)雙擊列表中的第一條波譜曲線同一繪圖窗口中,繪制出未知地物的波譜曲線以及波譜庫中的波譜曲線,以進(jìn)行比較。該繪制圖如下圖所示:未知地物與波譜庫中最匹配的kaolinite的波

18、譜曲線的比較圖,表明兩者高度匹配72/JqrvelcnqtiIdentify:Knownvs.Unkinowri2AltEdI:OpbansPfetFLhicten刪使用波譜分析(SpectralAnalyst),比較波譜曲線圖,對剛才提取的影像波譜曲線進(jìn)行鑒別,查看其為哪一種礦物質(zhì)。當(dāng)已經(jīng)鑒別了幾種礦物質(zhì)后,就可以進(jìn)入下一部分了。也可以選擇將USGS波譜庫中的波譜曲線同影像波譜曲線以及JPL波譜庫中的波譜曲線進(jìn)行比較USGS波譜庫所對應(yīng)的波譜數(shù)據(jù)文件為usgs_min.sli。5定義感興趣區(qū)感興趣區(qū)(ROIs)用來提取像素集合的統(tǒng)計(jì)信息及其平均波譜曲線??梢栽谌魏物@示的影像上定義足夠多的感興

19、趣區(qū)。從主影像窗口菜單中,選擇OverlayRegionofInteres打開ROITool對話框。5.1創(chuàng)建新的感興趣區(qū):在影像中,點(diǎn)擊鼠標(biāo)左鍵。在多邊形頂點(diǎn)處,點(diǎn)擊鼠標(biāo)左鍵,繪制感興趣區(qū),或者點(diǎn)擊并按住鼠標(biāo)左鍵,移動(dòng)鼠標(biāo),連續(xù)繪制感興趣區(qū)。點(diǎn)擊鼠標(biāo)右鍵,封閉該多邊形,完成感興趣區(qū)的定義。再次點(diǎn)擊鼠標(biāo)右鍵鎖定感興趣區(qū)的位置。點(diǎn)擊Stats按鈕,計(jì)算感興趣區(qū)的統(tǒng)計(jì)信息,繪制均值波譜曲線(白色),均值波譜上下各有一條的標(biāo)準(zhǔn)差曲線(綠色),以及最小和最大值的包絡(luò)波譜曲線(紅色),其包含了感興趣區(qū)所有的波譜。遙感原理與應(yīng)用- -轡ROIStaitEtksResults:zhumuqiarr95eH.

20、mt4QDOpicns20h_QLL2DOSilcaLPlot-ClearPitt-2Min/Max/Me-an:代e甲口口#1”已d怨比pnts2J3:LcctdIQJJid!-H-JuJJ.ncion門n.n!門BalBalBalBalBalBalBnlBalEal口日tL32no21d3o1i44456Fb6&_bM222222222K59573fl99oBo9o133334BM434A-4i4_4422A51S21-J3299.289266303.41423631D.1111296315.0528573164S9206315.907143315.4142B631&5352365-b-b

21、-b-bnr-H-2222222Etde359750762343B2G1473997D2974310153034659154&60遙感原理與應(yīng)用- -在FileStatisticsReport對話框中,選擇FileCancel。然后在AvgSpectrum繪圖窗口中,也選擇FileCancel。點(diǎn)擊ROITool對話框中的Delete按鈕,刪除所選的感興趣區(qū)。5、2加載預(yù)先保存的感興趣區(qū)文件:(無感興趣區(qū)文件,用自己定義的感興趣區(qū)進(jìn)行實(shí)驗(yàn))為已知特定礦物質(zhì)預(yù)先定義的感興趣區(qū)就會(huì)列在ENVI的message對話框中,然后還將加載到ROITool對話框,如下圖所示:ROITool對話框,展示了恢復(fù)

22、出的感興趣區(qū)(左圖)選擇ROITool對話框頂部的Off單選按鈕,允許在主影像顯示窗口中選取像素的位置。從主影像窗口中,選ToolsProfilesZ-Profile(Spectrum),打開Z剖面廓線窗口(Z-Profile)。使用鼠標(biāo)右鍵點(diǎn)擊感興趣區(qū)中的某個(gè)像素,將當(dāng)前像素的位置或者光標(biāo)的位置移動(dòng)到每一個(gè)感興趣區(qū)上。點(diǎn)擊感興趣區(qū)中的不同像素,移動(dòng)光標(biāo)的位置,在SpectralProfile窗口中,顯示新的波譜剖面廓線。注意:對于每個(gè)新的感興趣區(qū),y軸的繪制范圍都會(huì)自動(dòng)地調(diào)整大小,以保證與波譜剖面廓線相匹配。通過以上操作可以查看每個(gè)感興趣區(qū)中波譜曲線的差異。5、3從感興趣區(qū)中提取均值波譜曲線

23、:在ROITool對話框中,選擇某個(gè)感興趣區(qū),點(diǎn)擊Stats按鈕,提取所需的統(tǒng)計(jì)信息,并且繪制所選感興趣區(qū)的波譜曲線。將感興趣區(qū)的均值波譜同相應(yīng)的標(biāo)準(zhǔn)差波譜(綠色,在均值波譜的上面或者下面各有一條)以及包絡(luò)波譜(紅色,也在均值波譜的上面或者下面各有一條)進(jìn)行比較,查看感興趣區(qū)波譜的變化。對每個(gè)感興趣區(qū)重復(fù)上面的步驟。如果想要從zhumuqianl.sli波譜庫中加載相應(yīng)的波譜曲線到繪圖窗口,進(jìn)行直接地比較或者鑒別,那么在加載波譜庫的波譜數(shù)據(jù)時(shí)不要忘記將Y的縮放系數(shù)改為1000。完成這些操作后,關(guān)閉所有的FileStatisticsReport對話框以及所有的繪圖窗口。在ROTTool對話框中,

24、選擇OptionsMeanforAllRegions,在同一個(gè)繪圖窗口中,繪制每個(gè)感興趣區(qū)的均值波譜曲線。5、4鑒別礦物質(zhì)設(shè)計(jì)選擇合適的彩色合成來鑒別礦物質(zhì):在可用波段列表對話框中,選擇RGB單選按鈕,并順次點(diǎn)擊Band183、Band193和Band207。點(diǎn)擊LoadRGB,將彩色合成影像加載到當(dāng)前影像顯示窗口中。在主影像窗口中,點(diǎn)擊ToolsProfilesZ-Profile(Spectrum)。注意到0在Z剖面廓線(Z-Profile)窗口中,被用來合成彩色影像的三波段的位置上,分別用紅、綠、藍(lán)三條垂直線標(biāo)出。在ROITool對話框中,選擇Off單選按鈕,使用Z剖面廓線窗口,在主影像窗

25、口中的感興趣區(qū)上或其附近瀏覽相應(yīng)的波譜曲線。通過使用鼠標(biāo)左鍵,點(diǎn)擊并拖動(dòng)波譜剖面廓線窗口中的顏色條(plotbars,將其拖動(dòng)到所需的波段位置上。注意:突出影像上特定地物的一種方法是將某個(gè)彩色條放置在某個(gè)特征吸收中心,而把另兩個(gè)彩色條放置在相對的波譜峰值上。在Z剖面廓線(ZProfile)窗口內(nèi),雙擊鼠標(biāo)左鍵,將新選擇的波段加載到顯示窗口中。查看了幾個(gè)地物點(diǎn)后,就會(huì)明白彩色合成后的顏色是如何同波譜特征相一致的。例如,alunitic區(qū)域在RGB彩色合成影像上將顯示出紅紫色(magenta),這是因?yàn)榫G色對應(yīng)的波段位于alunitic的吸收帶,這將產(chǎn)生較低的綠波段值,而紅色和藍(lán)色對應(yīng)的波段反射率

26、幾乎相等。紅色和藍(lán)色混合的結(jié)果就導(dǎo)致了含有aluite的像素顏色為紅紫色。6二維散點(diǎn)圖6、1查看二維散點(diǎn)圖:在主影像窗口中,選擇Toolsf2-DScatterPlots,打開二維的scatterpl(對話框,繪制表觀反射率影像的散點(diǎn)圖。在ChoosandX列表中,選擇波段193(band193),在ChooseBandY列表中,選擇波段207。3.點(diǎn)擊OK。左圖的散點(diǎn)圖就顯示了X、Y表觀反射率值的關(guān)系波段193和波段207的散點(diǎn)圖(ATREM表觀反射率影像)(左圖)FiledasOpbenshdp#1SuiltvrPluL遙感原理與應(yīng)用- #-遙感原理與應(yīng)用- -115n131Derd(2r

27、200ihumuqianBeffrini=引6、2散點(diǎn)圖的密度分割:從散點(diǎn)圖窗口菜單中,選擇OptionsfDensityOn,自動(dòng)對散點(diǎn)圖進(jìn)行密度分割處理。散點(diǎn)的顏色表明了用來繪制散點(diǎn)圖的兩個(gè)波段特定表觀反射率組合出現(xiàn)的頻率。紫色(Purple)代表了出現(xiàn)的頻率最低。隨著顏色從藍(lán)色、綠色、黃色,變化到紅色,頻率也依次變大。紅色代表了出現(xiàn)的頻率最高。從散點(diǎn)圖菜單中,選擇OptionsDensityOff,關(guān)閉密度分割功能。6、3散點(diǎn)圖中的跳躍像素(SCATTERPLOTDANCINGPIXELS):從主影像窗口中點(diǎn)擊鼠標(biāo)左鍵并拖動(dòng)鼠標(biāo),觸發(fā)散點(diǎn)圖中的“跳躍像素”(DancingPixels)。

28、散點(diǎn)圖中的紅色點(diǎn)就對應(yīng)于主影像窗口光標(biāo)周圍10X10區(qū)域內(nèi)的像素點(diǎn)。預(yù)測一下某些影像顏色在散點(diǎn)圖中的位置,然后進(jìn)行核實(shí)。注意散點(diǎn)圖中紅色跳躍像素構(gòu)成的子散點(diǎn)圖的形狀。選擇OptionsSetPatchSize,改變散點(diǎn)圖中紅色矩形光標(biāo)的大小,并觀察改變后的差異。6、4影像中的跳躍像素(IMAGEDANCINGPIXELS):在散點(diǎn)圖上點(diǎn)擊鼠標(biāo)中鍵并拖動(dòng)鼠標(biāo),觸發(fā)主影像窗口中的“跳躍像素”(DancingPixels)。影像中的紅色像素點(diǎn)就對應(yīng)于散點(diǎn)窗口光標(biāo)周圍10X10區(qū)域內(nèi)的散點(diǎn)。注意所選像素的空間分布及其空間一致性。選擇OptionsSetPatchSize,改變散點(diǎn)圖中矩形光標(biāo)的大小,并

29、觀察改變后的差異。6、5將散點(diǎn)圖連接到波譜剖面廓線(SCATTERPLOTSLINKEDTOASPECTRALPROFILE):從散點(diǎn)圖窗口菜單中,選擇OptionsZProfile,選擇一個(gè)輸入文件用于提取波譜剖面廓線,然后點(diǎn)擊OK。打開一個(gè)連接到該散點(diǎn)圖的空白的波譜剖面廓線繪制窗口。將鼠標(biāo)光標(biāo)放置到二維散點(diǎn)圖中,點(diǎn)擊鼠標(biāo)右鍵,提取具有該散點(diǎn)位置特性的相應(yīng)空間像素點(diǎn)的波譜曲線。比較散點(diǎn)圖上的其它點(diǎn),觀察散點(diǎn)圖上相對中心位置的點(diǎn),記錄這些點(diǎn)所屬的類別。遙感原理與應(yīng)用- -遙感原理與應(yīng)用- -#1ScatterPlol蟹ScatterPlotPrcrfte:zhuniuqurr9faba斗32日

30、n-32.02.4WavelengthFibEdtOptioreHelp0DQQGLn-r;-3025456、6散點(diǎn)圖中的感興趣區(qū)(SCATTERPLOTROIS):散點(diǎn)工具也可以用來作為一個(gè)快速的分類器。在散點(diǎn)圖中點(diǎn)擊鼠標(biāo)左鍵,選擇感興趣區(qū)(ROI)的第一個(gè)頂點(diǎn)。使用鼠標(biāo)左鍵,選擇所需的線段頂點(diǎn),在散點(diǎn)圖中,繪制感興趣區(qū)多邊形。點(diǎn)擊鼠標(biāo)右鍵封閉該多邊形。與該多邊形包含的散點(diǎn)相對應(yīng)的影像像素,在主影像窗口中會(huì)以彩色編碼后的顏色(如紅色)顯示出來。從散點(diǎn)圖窗口的Class下拉式菜單中,選擇其它的顏色。繪制其它的感興趣區(qū)多邊形,影像中對應(yīng)的像素就會(huì)以所選顏色高亮突出顯示出來。要移

31、除某個(gè)類,可以選擇OptionsClearClass。此外,通過在繪圖坐標(biāo)軸以外(下面)點(diǎn)擊鼠標(biāo)中鍵,也可以移除當(dāng)前所選的類(Class)。使用二維散點(diǎn)圖工具,由散點(diǎn)圖中的類追溯到影像上,查看某些對應(yīng)像素出現(xiàn)在影像上的位置。在散點(diǎn)圖窗口菜單欄中,選擇OptionsExportClassorExportAll,將類(Classes)轉(zhuǎn)換成感興趣區(qū),然后再作為訓(xùn)練樣區(qū),使用所有的影像波段,參與影像分類。以這種方式導(dǎo)入的感興趣區(qū)將會(huì)列出在ROITool對話框中,它對后續(xù)的監(jiān)督分類都有效。從ENVI主菜單中,選擇ClassificationCreateClassImagefromROIs,也可以將感興

32、趣區(qū)轉(zhuǎn)換為分類影像。在散點(diǎn)圖窗口中,選擇OptionsClearAll,移除散點(diǎn)圖中所選的類以及影像中相對應(yīng)的像素點(diǎn)的顏色。6、7影像中的感興趣區(qū)IMAGEROIS):散點(diǎn)圖工具也可以從影像上獲取感興趣區(qū),以作為簡單的分類器。在散點(diǎn)圖窗口菜單中,選擇OptionsImageROI。使用鼠標(biāo)左鍵來繪制線條,在主影像窗口中繪制出感興趣區(qū)多邊形。使用鼠標(biāo)右鍵封閉多邊形或者鎖定多邊形的位置。然后影像感興趣區(qū)多邊形內(nèi)的像素點(diǎn)就會(huì)映射到散點(diǎn)圖中相對應(yīng)的位置上,并用當(dāng)前所選顏色高亮顯示出來。相應(yīng)地,影像中所有與散點(diǎn)圖高亮顯示點(diǎn)相匹配的像素點(diǎn),也將反向映射到主影像窗口中,即用相同的顏色高亮顯示出來,這就仿佛是

33、繪制了散點(diǎn)圖感興趣區(qū)。這是一個(gè)用來形成兩波段的分類最簡單的方法,此外它也是一個(gè)強(qiáng)大的工具。繪制一些影像感興趣區(qū),注意影像顏色與散點(diǎn)圖特性之間的一致性。6、8散點(diǎn)圖和波譜混合(SCATTERPLOTSANDSPECTRALMIXING):在散點(diǎn)圖窗口菜單中,選擇OptionsChangeBands,修改散點(diǎn)圖繪制中所需的兩個(gè)波段組合。嘗試選擇至少一對臨近的波段組合以及其它相距較遠(yuǎn)的波譜波段組合。BandX:Band176BandY:Band215如下圖四、實(shí)驗(yàn)體會(huì)本次實(shí)驗(yàn)的主要目的在于使用ENVI先進(jìn)的高光譜工具對多光譜數(shù)據(jù)進(jìn)行分析。更好地理解高光譜處理的概念及其工具。本專題將從特定礦物質(zhì)的感興

34、趣區(qū)中提取其波譜曲線,并與波譜庫中的波譜曲線進(jìn)行比較,找出顯示波譜信息的最佳RGB彩色組合。最終通過自己的操作也達(dá)到了以上目的。實(shí)驗(yàn)過程中出現(xiàn)了對理論知識(shí)理解不透徹而不知所以的現(xiàn)象,通過閱讀相關(guān)課件慢慢熟悉,順利完成本次實(shí)驗(yàn)。實(shí)驗(yàn)二影像鑲嵌一、實(shí)習(xí)目的掌握使用ENVI影像鑲嵌的知識(shí)。影像鑲嵌是一門藝術(shù),它將把多幅影像連接合并,以生成一幅單一的合成影像。ENVI提供了交互式的方式來將沒有地理坐標(biāo)(Nogeoreferenced)的影像拼接在一起,或者自動(dòng)地拼接有地理坐標(biāo)的影像,并輸出成有地理坐標(biāo)的鑲嵌影像。鑲嵌程序提供了透明處理、直方圖匹配,以及顏色自動(dòng)平衡的功能oENVI的虛擬鑲嵌能力還允許用

35、戶有選擇性地生成并顯示鑲嵌影像,而不用生成一個(gè)很大的輸出文件。二、實(shí)驗(yàn)步驟1、基于像素(Pixel)的影像鑲嵌將以像元為基礎(chǔ)的圖像鑲嵌起來,MapMosaicking一PixelBased。操作步驟:選擇MapMosaicking一PixelBased。打開PixelBasedMosaic鑲嵌界面選擇Import一ImportFiles選擇要鑲嵌的影像在MosaicInputFiles對話框中,點(diǎn)擊OpenFile,選擇文件zhumuqian_2.img,zhumuqian_3.img。在MosaicInputFiles對話框中,按下鍵盤上的Shift(也可用Ctrl)鍵,并同時(shí)點(diǎn)擊zhumu

36、qian_2.img和zhumuqian_3.img文件名,選中這兩個(gè)文件,點(diǎn)擊0K。在SelectMosaicSize對話框的XSize中輸入614,YSize中輸入1024,指定鑲嵌影像的大小,單擊Ok。在PixelBasedMosaic對話框中,點(diǎn)擊dv06_3.img文件名。影像當(dāng)前的位置就會(huì)以像素值為單位,列在對話框底部的文本框中。(右圖)在YO文本框中,輸入值513,按下鍵盤上的Enter鍵。這樣zhumuqian_3.img影像就會(huì)放置在zhumuqian_2.img影像的下面。(下圖)PixelMosaict14x1024常1Band.1:zhuTiiiqian3.imgRed

37、#HOFileInportCptisnsHelpflfZE:iTLii1:zhuTiuqi:m_Z.lmgGreen選擇或鍵入鑲嵌尺寸的大小一OK。在PixelBasedMosaic對話框中,選擇FilefApply。當(dāng)MosaicParameters對話框出現(xiàn)后,輸入輸出文件名zhumuqian.img,點(diǎn)擊OK,生成鑲嵌影像文件。2、邊緣羽化的處理操作在PixelBasedMosaic對話框中,選擇OptionsChangeMosaicSize。在SelectMosaicSize對話框的XSize和YSize文本框中都輸入值768,點(diǎn)擊OK,改變輸出鑲嵌影像的大小。在PixelBasedM

38、osaic對話框中,左鍵點(diǎn)擊影像#2的綠色輪廓框。將影像#2拖動(dòng)到鑲嵌圖的右下角。在鑲嵌圖中,右鍵點(diǎn)擊影像#1的紅色輪廓框,選擇EditEntry,打開Entry:filename對話框。在DataValuetoIgnore文本框中,輸入值0。在FeatheringDistance文本框中,輸入值25,點(diǎn)擊OK。對另一幅影像,重復(fù)上面的兩步操作。結(jié)果如圖選擇FilefSaveTemplate,輸入輸出文件名zhumuqianB.mos。在可用波段列表中,點(diǎn)擊鑲嵌模板文件名,然后點(diǎn)擊LoadBand,顯示該鑲嵌影像。當(dāng)使用虛擬鑲嵌時(shí),不會(huì)進(jìn)行羽化處理。3、基于地理坐標(biāo)的影像鑲嵌3、1創(chuàng)建基于地理

39、坐標(biāo)的鑲嵌影像在ENVI主菜單中,選擇MapfMosaickingfGeoreferenced,開始進(jìn)行ENVI基于地理坐標(biāo)的鑲嵌操作。文件選用zhumuqian_01w.img;zhumuqian_02w.img在MapBasedMosaic對話框中選擇FilefRestoreTemplate,然后選擇文件lch_a.mos(數(shù)據(jù)缺失),打開文件,恢復(fù)基于地理坐標(biāo)并進(jìn)行羽化鑲嵌時(shí)所必需的參數(shù)。4、鑲嵌時(shí)的色彩均衡專題的這一部分將介紹使用自動(dòng)顏色平衡,來生成一幅基于地理坐標(biāo)的鑲嵌影像。對于本次操作,我們將使用Landsat7ETM影像的兩個(gè)有重疊區(qū)域的影像子集。其中一個(gè)子集已經(jīng)進(jìn)行了直方圖均衡

40、化拉伸,它的每一個(gè)波段都獨(dú)立地處理過,這樣兩個(gè)子集相應(yīng)波段的直方圖都會(huì)有完全不同的形狀。4、1創(chuàng)建未經(jīng)色彩均衡化處理的鑲嵌影像:首先,創(chuàng)建生成一個(gè)沒有進(jìn)行顏色平衡處理的鑲嵌影像。導(dǎo)入兩幅沒有羽化過的影像,這樣就可以清楚地看出影像間的銜接線。在ENVI主菜單中,選擇MapfMosaickingfGeoreferenced,開始進(jìn)行ENVI基于地理坐標(biāo)的鑲嵌操作。MapBasedMosaic對話框出現(xiàn)在屏幕上。在MapBasedMosaic對話框中,選擇ImportfImportFiles,在文件選擇對話框中,點(diǎn)擊OpenFile按鈕,選擇zhumuqian_equal.dat文件,點(diǎn)擊Open。

41、在MosaicInputFiles對話框中,繼續(xù)點(diǎn)擊OpenFile按鈕,選擇zhumuqian2.dat文件,點(diǎn)擊Open。在MosaicInputFiles對話框中,選中mosaic_2.dat文件,然后按住Shift鍵,選擇zhumuqian_equal.dat文件,點(diǎn)擊OK。兩幅影像將被自動(dòng)地放置在MosaicPreview預(yù)覽顯示窗口中正確的地理坐標(biāo)位置上。在預(yù)覽顯示窗口中,影像都按默認(rèn)的2%對比度進(jìn)行獨(dú)立地拉伸。RGB彩色鑲嵌預(yù)覽在預(yù)覽顯示中,也可以顯示多光譜RGB彩色合成后的影像。FieImportOptionsHdp#2d&AIkpzn.OpenImageFile。當(dāng)Enter

42、DataFilename文件選擇對話框出現(xiàn)后,選擇進(jìn)入envidata目錄下的zhumuqian_reg子目錄,從歹列表中選擇zhumuqian_sp.img文件。點(diǎn)擊OK。Flih.unuqian_3p.insijaarbfaroneadSFfl7jdHEls當(dāng)可用波段列表對話框出現(xiàn)后,點(diǎn)擊GrayScale單選按鈕,使用鼠標(biāo)左鍵,點(diǎn)擊相應(yīng)的波段名,從對話框頂部所列波段中選中SPOT波段。所選擇的波段名顯示在SelectedBand:字段區(qū)域中。點(diǎn)擊LoadBand按鈕,加載這幅影像到一個(gè)新的顯示窗口中。1、2修改ENVI頭文件中的地圖信息:在可用波段列表中,右鍵點(diǎn)擊zhumuqian_sp

43、.img文件名下的MapInfo圖標(biāo),從彈出的快捷菜單中選擇EditMapInformation。EditMapInformation對話框出現(xiàn)在屏幕上。這個(gè)對話框列出了在ENVI中添加地理坐標(biāo)所用的地理信息??梢哉{(diào)整ENVI使用的MagicPixel(作為地圖坐標(biāo)系統(tǒng)的起始點(diǎn))相對應(yīng)的影像坐標(biāo)。因?yàn)镋NVI可以從相應(yīng)頭文件信息和地圖投影文件中,識(shí)別出地圖投影、像元大小以及地圖投影參數(shù),所以用它能夠計(jì)算出影像中任意像元的地理坐標(biāo)。既可以輸入地圖坐標(biāo),也可以輸入地理坐標(biāo)(緯度/經(jīng)度)。點(diǎn)擊Projection/Datum文本旁邊的箭頭切換按鈕,顯示UTMZone13North地圖投影的緯度/經(jīng)度

44、坐標(biāo)。ENVI在處理過程中才進(jìn)行轉(zhuǎn)換。點(diǎn)擊當(dāng)前的DMS或者DDEG按鈕,分別在度-分-秒(Degrees-Minutes-Seconds)和十進(jìn)制的度(DecimalDegrees)之間進(jìn)行切換。點(diǎn)擊Cancel,退出EditMapInformation對話框。EditMapInformation對話框1、3光標(biāo)位置/值:要打開一個(gè)顯示主影像窗口,滾動(dòng)窗口,或者縮放窗口中光標(biāo)位置信息的對話框,可以按如下幾步進(jìn)行。從主影像窗口菜單欄中,選擇ToolsfCursorLocation/Value。也可以從ENVI主菜單和主影像窗口菜單欄中,選擇WindowCursorLocation/Value,打

45、開這個(gè)對話框。CursorLocation對話框,顯示了像素值和相應(yīng)的地理坐標(biāo)(對帶地理坐標(biāo)的影像而言)注意到對于這幅帶地理坐標(biāo)的影像,對話框同時(shí)給出了像素坐標(biāo)和地理坐標(biāo)。在影像中移動(dòng)光標(biāo),查看特定位置的坐標(biāo)值,并注意地圖坐標(biāo)和經(jīng)緯度之間的關(guān)系。1、4疊合地圖公里網(wǎng):從主影像窗口菜單欄中,選擇OverlayfGridLines。#1GridLineParameters對話框出現(xiàn)在屏幕上,同時(shí)一個(gè)虛擬的邊框添加到影像中,允許在影像外部顯示地圖公里網(wǎng)的標(biāo)注。GridLineParameters對話框在這個(gè)新的對話框中,選擇FileRestoreSetup。在EnterGridParametersF

46、ilename對話框中,選中zhumuqian_sp.grd文件,點(diǎn)擊Open。先前保存過的公里網(wǎng)參數(shù)就會(huì)被加載到對話框中。在#1GridLineParameters對話框中,從對話框菜單欄中,選擇OptionsEditMapGridAttributes,來查看地圖參數(shù)。這將打開EditMapAttributes對話框。在EditMapAttributes對話框中,注意公里網(wǎng)的間隔以及控制線條、標(biāo)簽、公里網(wǎng)交角以及矩形框(輪廓框)相應(yīng)顏色和其它特征的參數(shù)。處理完成后,點(diǎn)擊Cancel,關(guān)閉該對話框。在#1GridLineParameters對話框,從菜單欄中,選擇OptionsEditGeog

47、raphicGridAttributes,查看地理坐標(biāo)。在EditGridAttributes對話框中,再次注意地理坐標(biāo)(緯度/經(jīng)度)公里網(wǎng)的參數(shù)。處理完成后,點(diǎn)擊Cancel來關(guān)閉該對話框。在GridLineParametersdialog對話框中,點(diǎn)擊Apply,在影像中放置公里網(wǎng)。ENVI允許同時(shí)放置像素、地圖、和地理坐標(biāo)公里網(wǎng)。1、4、1疊合地圖注記:在主影像窗口中,選擇OverlayAnnotation。在#1Annotation:Text對話框中,選擇FileRestoreAnnotation。打開一個(gè)標(biāo)準(zhǔn)的文件選擇對話框。在EnterAnnotationFilename對話框中,

48、從文件列表中選擇zhumuqian_sp.ann文件,點(diǎn)擊OK。先前保存過的地圖注記被加載到影像上。按住滾動(dòng)窗口的一角,并拖動(dòng)鼠標(biāo),拉大該滾動(dòng)窗口。重新放置改變了大小的滾動(dòng)窗口,這樣就可以同時(shí)看到主影像窗口。在改變了大小的滾動(dòng)窗口中,使用鼠標(biāo)左鍵,移動(dòng)主影像指示矩形框,查看主影像窗口中出現(xiàn)的地圖要素。在#1Annotation:Text對話框中,點(diǎn)擊并按住Object菜單,查看可以用來注記地圖的對象。1、5輸出到影像或Postscript文件ENVI對用戶提供了幾個(gè)選項(xiàng),來保存和輸出地圖影像??梢詫⒔Y(jié)果保存為ENVI的影像文件格式,或者保存為幾種通用的圖形格式(包括Postscript格式)。

49、也可以直接打印或者導(dǎo)入到其它軟件中。1、5、1將影像保存為ENVI的影像格式:要將結(jié)果影像保為ENVI自身的格式(比如一個(gè)RGB文件),按下面的步驟進(jìn)行操作。在主影像窗口中,選擇FileSaveImageAsImageFile。在OutputDisplaytoImageFile對話框中,選擇OutputFileType的下拉式按鈕(默認(rèn)情況下文件類型設(shè)置為ENVI),查看可用的不同格式。ChangeGraphicsOverlaySelections按鈕可以打開一個(gè)同樣名字的對話框,這個(gè)對話框允許添加或刪除許多制圖疊合選項(xiàng)(graphicsoptions),包括注記和公里網(wǎng)。ChangeImag

50、eBorderSize按鈕也可以打開一個(gè)同名的對話框。這個(gè)對話框允許改變頂部、底部、左邊和右邊的邊框?qū)挾龋绻枰?,也可以改變邊框的顏色。如果帶注記和公里網(wǎng)的彩色影像已經(jīng)顯示在顯示窗口中,那么注記和公里網(wǎng)都將自動(dòng)地列在制圖疊合選項(xiàng)中。同時(shí)也可以選擇其它要疊置在輸出影像上的注記文件。通過選擇Memory或者File單選按鈕,決定是將結(jié)果保存為一個(gè)磁盤文件,還是保存到內(nèi)存中。選擇Memory,點(diǎn)擊OK,輸出影像。新生成的影像文件自動(dòng)列在可用波段列表中。在可用波段列表中,點(diǎn)擊Display#1下拉式按鈕,從菜單中選擇NewDisplay,打開一個(gè)新的顯示窗口。選擇RGBColor單選按鈕,要將影像從

51、內(nèi)存中加載到顯示窗口,可以連續(xù)選擇R、G和B(帶地理坐標(biāo)的SPOT數(shù)據(jù))波段。點(diǎn)擊LoadRGB按鈕,添加注記后的影像作為一幅柵格圖顯示出來。2影像到影像的配準(zhǔn)這一部分將逐步演示影像到影像的配準(zhǔn)處理過程。帶有地理坐標(biāo)的SPOT影像被用作基準(zhǔn)影像,一個(gè)基于像素坐標(biāo)的LandsatTM影像將被進(jìn)行校正,以匹配該SPOT影像。2、1打開并顯示LandsatTM影像文件:從ENVI主菜單中,選擇FileOpenImageFile。當(dāng)EnterDataFilenames對話框出現(xiàn)后,選擇進(jìn)入envidata目錄下的zhumuqian_reg子目錄,從列表中選擇zhumuqian_tm.img文件。在文件

52、選擇對話框中,點(diǎn)擊Open,把TM影像波段加載到可用波段列表中。在列表中選中波段3,點(diǎn)擊NoDisplay按鈕,并從下拉式菜單中選擇NewDisplay。點(diǎn)擊LoadBand按鈕,來把TM第3波段的影像加載到一個(gè)新的顯示窗口中。2、2顯示光標(biāo)位置/值要打開一個(gè)顯示主影像窗口,滾動(dòng)窗口,或者縮放窗口中光標(biāo)位置信息的對話框,可以按一下步驟進(jìn)行操作。從主影像窗口菜單欄中,選擇ToolsCursorLocation/Value。在主影像窗口、滾動(dòng)窗口和縮放窗口的TM影像上,移動(dòng)鼠標(biāo)光標(biāo)。注意到坐標(biāo)是以像素單位給出的,這是因?yàn)檫@個(gè)影像是基于像素坐標(biāo)的,它不同于上面帶有地理坐標(biāo)的SPOT影像。選擇File

53、Cancel,關(guān)閉CursorLocation/Value對話框。2、3開始進(jìn)行影像配準(zhǔn)并加載地面控制點(diǎn)從ENVI主菜單欄中,選擇MapRegistrationSelectGCPs:ImagetoImage。在ImagetoImageRegistration對話框中,點(diǎn)擊并選擇Display#1(SPOT影像),作為BaseImage。點(diǎn)擊Display#2(TM影像),作為WarpImage。點(diǎn)擊OK,啟動(dòng)配準(zhǔn)程序。通過將光標(biāo)放置在兩幅影像的相同地物點(diǎn)上,來添加單獨(dú)的地面控制點(diǎn)。用來進(jìn)行影像到影像配準(zhǔn)的GroundControlPointsSelection對話框在GroundControl

54、PointsSelection對話框的BaseX和Y文本框中,分別輸入753和826,將SPOT影像中的光標(biāo)移動(dòng)到相應(yīng)的點(diǎn)上。使用同樣的方法,在WarpX和Y文本框中,分別輸入331和433,將TM影像中的光標(biāo)移動(dòng)至肪目應(yīng)的點(diǎn)上。在兩個(gè)縮放窗口中,查看光標(biāo)點(diǎn)所處位置。如果需要,在每個(gè)縮放窗口所需位置上,點(diǎn)擊鼠標(biāo)左鍵,調(diào)整光標(biāo)點(diǎn)所處的位置。注意:在縮放窗口中支持亞像元(sub-pixel)級(jí)的定位??s放比例越大,地面控制點(diǎn)定位的精度就越好。在GroundControlPointsSelection對話框中,點(diǎn)擊AddPoint,把該地面控制點(diǎn)添加到列表中。點(diǎn)擊ShowList查看地面控制點(diǎn)列表。

55、嘗試選擇幾個(gè)地面控制點(diǎn)找到選擇地面控制點(diǎn)的感覺。注意對話框中所列的實(shí)際影像點(diǎn)和預(yù)測點(diǎn)坐標(biāo)。一旦已經(jīng)選擇了至少4個(gè)地面控制點(diǎn)以后,RMS誤差就會(huì)顯示出來。在GroundControlPointsSelection對話框中,選擇OptionsClearAllPoints,可以清除掉所有已選擇的地面控制點(diǎn)。從GroundControlPointsSelection對話框中,選擇FileRestoreGCPsfromASCII。在EnterGroundControlPointsFilename對話框中,選擇文件zhumuqian_tm.pts,然后點(diǎn)擊OK,加載這個(gè)預(yù)先保存過的地面控制點(diǎn)坐標(biāo)。在Ima

56、getoImageGCPList對話框中,點(diǎn)擊單獨(dú)的地面控制點(diǎn)。查看兩幅影像中相應(yīng)地面控制點(diǎn)的位置、實(shí)際影像點(diǎn)和預(yù)測點(diǎn)的坐標(biāo)以及RMS誤差。調(diào)整對話框的大小,觀察GroundControlPointsSelection對話框中所列的合計(jì)RMS誤差(RMSError)。影像到影像配準(zhǔn)中所用的ImagetoImageGCPLIst對話框2、4操作處理地面控制點(diǎn)下面的內(nèi)容僅提供處理方法,并且只對有限的地面控制點(diǎn)按鈕的處理功能進(jìn)行操作。在ImagetoImageGCPList對話框中,選擇相應(yīng)的地面控制點(diǎn),然后在GroundControlPointsSelection對話框中進(jìn)行修改,這樣可以編輯單個(gè)

57、控制點(diǎn)的坐標(biāo)位置??梢酝ㄟ^輸入一個(gè)新的像素坐標(biāo),或使用對話框中的方向箭頭逐像素地移動(dòng)坐標(biāo)位置。在ImagetoImageGCPList對話框中,點(diǎn)擊0n/0ff按鈕,屏蔽掉所選擇的地面控制點(diǎn)。這樣在校正模型和RMS計(jì)算中都將不會(huì)考慮這個(gè)地面控制點(diǎn)坐標(biāo)。這些地面控制點(diǎn)并沒有被真正地刪除,僅僅是被忽略掉了,可以使用On/Off按鈕重新激活這些地面控制點(diǎn)。在ImagetoImageGCPList對話框中,點(diǎn)擊Delete按鈕,可以從列表中刪除一個(gè)地面控制點(diǎn)。在兩個(gè)縮放窗口中調(diào)整光標(biāo)位置,然后點(diǎn)擊ImagetoImageGCPList對話框中的Update按鈕,更新所選的地面控制點(diǎn),將其修改到當(dāng)前光標(biāo)

58、的所在位置。ImagetoImageGCPList對話框中的Predict按鈕,允許對新的地面控制點(diǎn)進(jìn)行預(yù)測。它以當(dāng)前的校正模型為基礎(chǔ)。將包含SPOT影像的那個(gè)主影像窗口的光標(biāo)放置到一個(gè)新的位置上。然后點(diǎn)擊Predict按鈕,放置在TM影像上的光標(biāo)就會(huì)根據(jù)校正模型,移動(dòng)到預(yù)測的匹配點(diǎn)上去。通過在TM數(shù)據(jù)中,輕微地移動(dòng)光標(biāo),能夠?qū)λ崛〉奈恢命c(diǎn)進(jìn)行交互式的精確定位。在GroundControlPointsSelection對話框中,點(diǎn)擊AddPoint,把這個(gè)新的控制點(diǎn)添加到列表中。2、5校正影像我們可以校正顯示的影像波段,也可以同時(shí)校正多波段影像中的所有波段。這里我們僅對已顯示的波段進(jìn)行校正。

59、從GroundControlPointsSelection對話框中,選擇OptionsWarpDisplayedBand。在RegistrationParameters對話框中的WarpMethod按鈕菜單中,選擇RST。在Resampling的按鈕菜單中選擇NearestNeighbor重采樣法。輸入文件名zhumuqian_tml.wrp,點(diǎn)擊OK。重復(fù)步驟1和步驟2,還是使用RST校正法,但是要相應(yīng)地選擇Bilinear和CubicConvolution重采樣法。將結(jié)果分別輸出到zhumuqian_tm2.wrp和zhumuqian_tm3.wrp文件中。再一次重復(fù)步驟1和步驟2,這一次

60、選擇一次多項(xiàng)式Polynomial校正法,并使用CubicConvolution重采樣法。然后再選擇Delaunay三角網(wǎng)的Triangulation校正法,相應(yīng)地使用CubicConvolution重采樣法。將結(jié)果分別輸出到zhumuqian_tm4.wrp和zhumuqian_tm5.wrp文件中。肆#3Scroll(.二叵岡NearestNeighbor結(jié)果Bilinear結(jié)果CubicConvolution結(jié)果2、6比較校正結(jié)果使用動(dòng)態(tài)鏈接來比較校正結(jié)果:在可用波段列表中,點(diǎn)擊原始的TM波段影像名zhumuqian_tm.img,然后從菜單欄中,選擇FileCloseSelectedF

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲(chǔ)空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論