遥感原理与应用课程实习.docVIP

遥感原理与应用课程实习 一、实验目的练习使用ENVI软件对遥感影像进行基础高光谱分析，通过感兴趣区的波普信息学习分析目标地物。 二、实验数据介绍实验数据为ENVI4.7软件自带遥感影像数据，为及其相应头文件cap95eff.hdr（后改为和liuchunyan1.hdr），此数据子集共含有50个波段；ENVI格式的JPL波谱库jpl1.sli及其相应的头文件jpl.hdr(（后改为liuchunyan2.sli和liuchunyan2.hdr）)。 三、实验过程此次高光谱数据分析主要包含有以下几点内容波谱库操作、浏览和提取影像反射率波谱、ENVI感兴趣区的设定及进行彩色合成影像的选取，最终目的在于练习鉴别波谱类型。 （1.1）提取影像反射率波谱（1.1.1）启用ENVI4.1软件，选择主菜单FileOpen Image File，然后进入ENVI软件的Data包中，选择并打开。 它将列出50个波段的名字。 （1.1.2）点击Gray Scale单选按钮，再在Available Bands List中选中Band193波段，点击Load RGB。 然后在主影像窗口菜单下选怎ToolsProfilesZ Profile（Spectrum），提取Band193波段的反射率波谱曲线图 （1）BandBand193波段的反射率波谱曲线（1.2）影像波谱与波谱库的比较（1.2.1）通过移动主影像显示框中的缩放指示矩形框，可以在#1Spectral Profile窗口中实时发现反射率曲线的变化，此说明不同地物在相同波段下其反射率是存在一定区别的。 如以下两个位置的光谱反射率曲线之间的差异（图2.1）光谱反射率曲线（图2.2）光谱反射率曲线对于（图2.1）和（图2.2）之间注意缩放指示矩形框的位置差异。 （1.2.2）在ENVI主菜单中选择SpectralSpectral LibrarisSpectral LibraryViewer；在Spectral LibraryInput File对话框中点击Open File按钮，选择下拉中的Spectral Library，在弹出的Please Selecta File文件夹中选择jpl1.sli，点击打开，返回到Spectral LibraryInput File对话框，选中jpl1.sli，点击OK；（1.2.3）在弹出的Spectral LibraryViewer对话框中，选择OptionsEdite(x,y)Scale Factors，将Y Multiplier设置为1.000使匹配影像表现反射率范围为（1-1000），点击OK。 （1.2.4）在Spectral LibraryViewer对话框中，选择以下4个波谱，绘制它们的波谱曲线为:ALUNTTE SO-4A BUDDINGTONTTEFELDS TS-11A CALCTTEC-3D KAOLINTTEWELL ORDEREDPS-1A图 （3）波谱库的波谱曲线(1.2.5)在波谱库波谱曲线显示屏上选择EditePlot Parameters，对弹出的对话框作如下设置图（4.1）s设置X-Axis图（4.2）s设置Y-Axis在Y-Axis中，Axis Title改为Reflectance，即Y轴上表示反射率。 点击Appy，然后Cancel。 使得上述波谱库中的4个波谱曲线变为图 （5）波谱库曲线（1.2.6）设置波谱曲线图的图例。 在绘制窗口中点击鼠标右键，从弹出的快捷菜单中选择Plot Key，并通过调节X Axis的左右范围，使图例名称全部显示如下图 （6）自定义绘图参数后的波谱曲线图（1.2.7）在Band193波段波谱曲线的显示窗口中单击鼠标右键，从弹出的快捷菜单中选择Plot Key，则得到的图例显示的是全影像图中缩放指示矩形框所在的X和Y坐标图 （7）显示坐标的波谱曲线（1.2.8）在#1Spectral Profile窗口中，选择OptionsNew WindowsBlank，打开一个新的绘制窗口。 （1.2.9）在主影像显示窗口中选择ToolsPixel Locator，在弹出的对话框中作Sample590，Line570设置，点击Apply，#1Spectral Profile窗口发生相应的变化；用鼠标左键按住其绘制图上的XY图例拖到新的绘图窗口。 如此重复操作，坐标（Sample，Line）分别为(435,555)，(494,514)，(531,541)，(502,589)，(448,505)，(260,613)。 最后选择新绘制窗口中的OptionsStack Plots得到新绘图窗口为图 （8）某些ATREM影像表现反射率波谱曲线图 （9）地物波谱曲线与波谱库曲线通过目视比较上面两个分别影像图中某些波谱反射率曲线图和波谱库中波谱曲线图，可以发现，实验室所测波谱曲线同某些影像表现反射率波谱曲线之间存在着相似的形状及吸收特性。 （1.2.9）上述影像波谱曲线与波谱库中某些影像波谱曲线相比较可得出如下结论基于它们的一定相似特性，说明具有与alunite、buddingtonite、calcite和kaolinite实验所测波谱曲线相似的影像波谱曲线是由上面相应的矿物质所构成的。 这对矿物质的初步识别鉴定具有重要的应用。 （1.3）波谱曲线的鉴别ENVI中Spectral Analyst可用于波谱曲线鉴别分析，其提供了一个波谱匹配工具，根据对波谱库中的波谱曲线对影像中的波谱曲线进行评分（数值为0-1），其中数值1表示它们相当于完全匹配。 （1.3.1）在ENVI主菜单中选择SpectralSpectral Analyst，在Spectral AnalystInput Spectral Library对话框中选择liuchunyan.sli，点击OK；随后在弹出的Edite IdentifyMethods Weighting对话框中点击OK。 （1.3.2）在主影像显示窗口菜单中选择ToolsProfilesZ Profile（Spectrum）。 然后再弹出的绘制窗口中鼠标右击，选择Plot Key，显示图例。 （1.3.3）在主影像窗口中选择ToolsPixel Locator，在弹出的对话框中输入坐标（Sample，Line）为（502，589），点击Apply。 （1.3.4）在Spectral Analyst对话框中，选择OptionsEdite MethodWeights。 在弹出对话框中每个Weight值均设为0.33，点击OK。 然后点击Spectral Analyst对话框的Apply，弹出以下多条波谱曲线列表波谱曲线列表选择X502，Y589所对应的波谱曲线。 Spectral Analyst对话框中便显示出如下这便是Spectral Analyst对影像中的波谱曲线进行的评分结果，根据上表，我们发现其与KAOLINITE WELLORDER的匹配分值为0.873，与其它波谱相比具有最大值，说明该像素所对应的地物最接近于KAOLINITE WELLORDER。 （1.3.5）用鼠标双击列表中的第一条波谱曲线。 在同一绘图窗口中绘制出地物的波谱曲线和波谱库中的波谱曲线，以进行直观比较图 （10）地物的波谱曲线和波谱库中的波谱曲线看图可易知此地物与KAOLINITE WELLORDER波谱具有高度匹配。 （1.3.6）如此重复便可以鉴别多种矿物质了，随后便可以进行感兴趣区分析了。 如图 (11)波谱评分鉴定（1.4）感兴趣区定义感兴趣区（ROIS）是提取像素集合的统计信息及其平均波谱曲线。 （1.4.1）任意打开某一波段的影像图，在主影像图窗口中选择OverlayRegion ofInterest，打开#1ROI Tool对话框。 （1.4.2）用鼠标左键在主影像图中画多边形，点击一次右键结束多边形，第二次右击使多边形闭合。 （1.4.3）在#1ROI Tool对话框中点击New Region定义第二个感兴趣区，同理作三个感兴趣区。 Region# 1、#2和#3。 （1.4.4）在#1ROI Tool对话框中点击Stats按钮，便得下图感兴趣统计图 （12）感兴趣区波谱信息其中白线表示均值波谱曲线，均值波谱曲线上下两条绿色曲线表示标准差曲线，两条红线则表示最大最小的包络波谱曲线，其包含了感兴趣区所有的波谱。 (1.4.5)选择ROI StatisticsResult对话框中选择FileCancel，关闭窗口。 （1.4.6）选择#1ROI Tool对话框下的FileSave ROIS，选择路径为envidata/cup_95av.roi。 （1.5）感兴趣区波谱差异（1.5.1）打开#1ROI Tool对话框，选择FileRestore ROIS，在envidata下找到cup95_av.roi文件，并打开。 这些感兴趣区（Sample，Line）坐标分别为（502，589）、（435，555）和（260，613），这些都是通过光谱曲线鉴别而基本特定了其矿物质组成，将它们均加载到#1ROI Tool对话框中。 （1.5.2）选择ROI Tool对话框中的Off按钮，允许在主影像显示窗口中选取像素的位置。 （1.5.3）打开主影像窗口主菜单ToolsProfilesZ-Profile(Spectral)，打开#1Spectral Profiles绘制窗口图 （13）波谱信息（1.5.4）使用鼠标点击感兴趣的不同像素，移动光标的位置，在#1Spectral Profiles窗口中会显示不同的光谱曲线，且y轴的绘制范围会自动地调整大小，以保证与波谱剖面轮廓线相匹配。 图 （14）感兴趣区内的波谱比较通过两图可清晰看出每个感兴趣区中波谱曲线的差异。 （1.6）感兴趣区波谱均值（1.6.1）在ROI Tool对话框中，任意选择一个感兴趣区，例如（435，555），点击Stats，便会出现统计信息，相对应的均值波谱、标准差波谱以及包络波谱曲线。 如下图图 （15）感兴趣区内的均值波谱曲线对其三类波谱曲线进行比较容易发现，它们具有极其相似的变化规律。 （1.6.2）在（1.3）波谱曲线鉴别步骤中，已经得知（435，555）处的地物与jpl.sli波谱库的TINCALCONIT B-4A的波谱曲线最匹配。 则在ENVI主菜单下选择SpectralSpectral LibrarySpectral LibraryViewerOpenjpl_libliuchunyan1.sli，打开，回到Spectral LibraryInput File对话框中，选择liuchunyan1.sli文件，在Spectral LibraryViewer列表中找到TINCALCONIT B-4A。 并在OptionsEdite（X,Y）Scale Factors下设置Y Multiplier值为1000。 即得到TINCALCONIT B-4A的波谱曲线并在Plot Parameters下对其设置相关参数。 将它与（435，555）感兴趣区的统计波谱曲线进行比较图 （16）地物的波谱曲线与不同评分结果的波谱曲线比较观察上面两图，SpectralLibraryPlots窗口中，红线为TINCALCONIT B-4A的波谱曲线，白色为TRONA C-4A波谱曲线。 （435，555）感兴趣区的波谱曲线分别与它们相比较，发现其与红曲线的变化规律更接近。 因此，综上得出结论尽管相对于波谱库中其它波谱，得分最高者的波谱曲线已经和地物波谱曲线规律是最匹配了，但仍有一定的差异，说明波谱库中的波谱匹配鉴别地物只能比较出它们之间的差异大小，而并不能确切的确定此地物的物质组成。 （1.7）二维散点图 （1）在主影像窗口中，选择Tools2-D ScatterPlots，打开二维的scatter plot，绘制反映反射率影像的散点图。 （2）在Choose BandX列表中，选择波段193，在Choose BandY列表中，选择Band207。 点击OK。 即得显示X、Y表观反射率值的关系图 （17）二维散点图 （3）在散点图窗口菜单中选择OptionsDensity On，对其进行密度分割。 图 （18）密度分割后的二维散点图 （4）而后在散点图菜单中选择OptionsDensity Off。 四、实验体会此次高光谱分析试验，其中做了提取影像反射率波谱 （1）、影像波谱与波谱库的比较 （2）、波谱曲线的鉴别 （3）、感兴趣区定义 （4）、感兴趣区波谱差异 （5）和感兴趣区波谱均值 （6）以及二维散点图7个实验。 在此实验中收获最大的是进一步学习了遥感影像的处理分析方法，理解了遥感影像的基础处理过程和原理，可以通过实验室波借助波谱库的评分鉴定方法初步、大致鉴定地物的矿物质组成成分。 这在现实寻矿工作中能发挥重要的作用。 实验二影像镶嵌 一、实验目的当研究区域在不同的图像文件上时，需要将不同的图像文件拼合在一起形成一幅完整的包含感兴趣区域的图像，这就是图像的镶嵌。 通过镶嵌处理，可以获取更大范围的地面图像。 参与镶嵌的图像可以是不同时间同一传感器获取的，也可以是不同时间不同传感器获取的，但要求镶嵌的图像之间要有一定的重叠度，并校正到同一地图坐标系下，即镶嵌的图像要有相同的地理参考。 在镶嵌之前必须先对所有待参与的图像进行几何校正，将它们校正到相同的坐标系下，再将具有地理参考的若干相邻图像合成一幅图像或一组图像。 对输入的图像不要求一定要有相同的投影类型或等大的像元，但必须具有相同的波段数；而后再进行色匀处理。 由此，通过此次实验将理论与实践相结合，在实际操作中进一步理解教材影像镶嵌处理知识，同时练习几何校正的方法，学会对遥感图像的相关处理工作。 二、实验数据介绍主要使用了ENVI软件自带的TM数据，主要有liuchunyan06_2.img与liuchunyan06_3.img，以及liuchunyan06a.mos，liuchunayn_lch_a.mos和LIUCHUNYAN_LCH_MOS.IMG等。 三、实验过程（2.1）基于像素的影像镶嵌（2.1.1）未羽化影像镶嵌 （1）打开ENVI4.7，在主菜单工具中选择MapMosaickingPixel Based，弹出Pixel BasedMosaic对话框。 （2）在对话框上选择ImportImport Files。 在Mosaic InputtFiles对话框中，点击Open，在数据中选择liuchunyan06_2.img，打开；继续打开liuchunyan06_3.img。 （3）在Mosaic InputtFiles对话框中，按shift键，同时选中liuchunyan06_2.img和liuchunyan06_3.img两幅影像，点击OK。 （4）在弹出的对话框中，X Size和Y Size分别设为614和1024，指定镶嵌影像的大小。 （5）在Pixel Mosaic6141024窗口中，点击liuchunyan06_3.img，在底部的Y0文本框中输入513，按键盘上的Enter键，便得下图，其中liuchunyan06_3.img在liuchunyan06_2.img的下部图 （1）镶嵌影像输入显示窗口 （6）在Pixel Mosaic6141024窗口中，选择FileApply；在弹出的Mosaic Parameters对话框中输入文件名liuchunyan06.img，点击OK，即生成镶嵌影像文件图 （2）影像镶嵌图 （7）在Pixel Mosaic6141024窗口中，选择FileSave Template；在弹出的Output MosaicTemplate对话框中输入文件名liuchunyan06a.mos。 在Available选中BandsList中它，点击Load Band，即输出生成虚拟镶嵌影像图图 （3）虚拟镶嵌影像图（2.1.2）羽化镶嵌图像 （1）如同2.1中，打开liuchunyan06_2.img和liuchunyan06_3.img影像。 在Pixel Mosaic6141024窗口中选择OptionsChange MosaicSize。 弹出如下对话框，并对其作图示设置点击OK，改变输出镶嵌影像的大小。 、 （2）在Pixel Mosaic6141024窗口中，左键点击#1的红色轮廓框，将影像#1拖到镶嵌图的右下角图 （3）镶嵌影像输入显示窗口 （3）在镶嵌图中，右键点击#1的红色轮廓框，选择Edite Entry，打开Entryfilename对话框。 在Data Valueto Ignore文本框中，输入0；在Feathering Distance文本框中输入25。 对另一影像做同样处理。 （4）在Pixel Mosaic6141024窗口中选择FileApply，输入文件名为liuchunyan06bby.mos，设定Background Value为255，点击OK。 （5）在Available BandsList中选中它，点击Load Band，即得到经羽化处理的影像图图 （4）经羽化处理的影像图（2.2）基于地理坐标的影像镶嵌（2.2.1）切割线和虚拟未羽化镶嵌影像 （1）在ENVI主菜单工具中选择MapMosaickingGeoreferenced。 在弹出的Map BasedMosaic对话框中选择FileRestore Templateliuchunayn_lch_a.mos文件，在Mosaic Input Files中选中liuchunayn_lch_a.mos，点击OK图 （5）影像文件输入图 （2）在Available BandsList中选择lch_o1w.img，点击波段名，Load Band，显示该影像。 （3）在主影像中，点击右键弹出快捷菜单，选择ToggleDisPlay ScrollBars。 （4）从主影像窗口中，选择OverlayAnnotation，打开Annotation对话框。 （5）在Annotation对话框中，选择FileRestore Annotation，然后选择lch_01w.ann文件。 这将显示出一条红色的切割线（cutline），该切割线用来在镶嵌影像中混合两影像。 （2.2.2）创建输出羽化后的镶嵌影像 （1）在M APB ASEDM OSAIC对话框中，选择F ILEA PPLY。 在M OSAICP ARAMETERS对话框中，输入输出文件名LIUCHUNYAN_LCH_MOS.IMG，点击OK，创建羽化后的镶嵌影像。 并将镶嵌影像加载到一个新的显示窗口中:图 （6）羽化后的镶嵌影像（2.2.3）镶嵌时的色彩均衡: (1)在ENVI主菜单中，选择MapMosaickingGeoreferenced，开始进行ENVI基于地理坐标的镶嵌操作。 Map BasedMosaic对话框出现在屏幕上。 在Map BasedMosaic对话框中，选择ImportImport Files。 在文件选择对话框中，点击Open File按钮，进入avmosaic目录，选择lcymosaic1_equal.dat文件，点击Open。 在Mosaic InputFiles对话框中，继续打开lcymosaic_2.dat文件，点击Open。 （2）在Mosaic InputFiles对话框中，选中lcymosaic_2.dat文件，然后按住Shift键，选择lcymosaic1_equal.dat文件，点击OK。 （3）在Map BasedMosaic显示窗口底部，右键点击lcymosaic1_equal.dat文件名，选择Edit Entry。 在随后出现的对话框中，将Mosaic Display的箭头切换按钮改为RGB。 （4）对于Red选择1，Green选择2，Blue选择3。 这将使得ENVI在预览镶嵌影像时，用红色电子枪来加载波段1，绿色电子枪来加载波段2，蓝色电子枪来加载波段3。 点击OK。 此时，lcymosaic1_equal.dat影像将会以彩色显示出来。 对要镶嵌的另一幅影像，重复以上步骤。 （5）从Map BasedMosaic显示窗口中，选择FileApply。 在随后出现的Mosaic Parameters对话框的Output Filename中，输入保存路径。 点击OK。 在RGB彩色显示窗口中，显示新镶嵌的影像。 将波段1作为红色，波段2作为绿色，波段3作为蓝色。 即得影像图为图 （7）彩色镶嵌影像通常两幅影像的衔接线都会很明显。 （6）在Map BasedMosaic对话框中，右键点击lcymosaic1_equal.dat文件名，选择Edit Entry。 在随后出现的对话框的底部，将Color Balancing改为Adjust。 这表明将调节这幅影像的对比度，使其与另一幅影像相匹配。 点击OK。 （7）右键点击lcymosaic_2.dat文件，选择Edit Entry。 在随后出现的对话框中，把Color Balancing改为Fixed。 这表明将不会改变这幅影像的对比度。 仅调节另一幅影像的对比度，使其与这幅影像匹配。 点击OK。 （8）从Map BasedMosaic显示窗口中，选择FileApply。 （9）在随后出现的Mosaic Parameters对话框底部，有个Color Balanceusing选项。 保留stats fromoverlapping regions这个选择项。 因为仅仅利用重叠区域统计计算值来进行颜色平衡效果会更好。 另一个选择项，stats fromplete files，用在镶嵌影像只有一点或者没有重叠区域的情况下。 选择保存路径，点击OK。 （10）镶嵌完成时，拼接的结果会加载到可用波段列表中。 在RGB彩色显示窗口中，显示新镶嵌的影像。 将波段1作为红色，波段2作为绿色，波段3作为蓝色。 显示影像图为图 （8）色匀处理后镶嵌影像两幅影像的衔接线将不会很明显。 四、实验体会镶嵌的用途在测绘领域非常广泛，只有当相邻测图很好吻合匹配才能获取效果好且信息失真度很低的地图。 通过此次影像镶嵌实验练习，学习了邻近影像之间的匹配镶嵌，认识了一些关于镶嵌处理的方法，为今后影像处理奠定了一定基础。 实验三影像配准 一、实验目的练习进行地理校正、添加地理坐标以及使用ENVI进行影像到哦影像的配和影像到地图的校正等方法。 二、实验数据介绍SPOT数据liuchunyan_sp.img TM数据liuchunyan_tm.img矢量图liuchunyan_rd.dlg 三、实验过程影像校正方法主要有重采样、缩放比例和平移，以及多项式函数或者Delaunay三角网法等。 其中重采样方法包括最近邻法（nearest-neighbor）、双线性内插法（bilinear interpolation）和三次卷积法（cubic convolution）。 使用ENVI的多重动态链接显示功能对基准影像和校正后的影像进行比较，可以快速地评估配准精度。 （3.1）带地理坐标的影像 （1）打开ENVI软件，在主菜单中选择FileOpen ImageFile，在打开的Enter DataFilenames对话框中找到liuchunyan_sp.img文件，并打开。 （2）在Available BandsList对话框中选中此SPOT数据，点击Gray Scale，单击Load Band，即显示SPOT遥感影像图图 （1）SPOT影像图 （3）在Available BandsList对话框中，右键点击liuchunyan_sp.img文件名下的Map Info图标，在弹出的快捷菜单下选择Edit MapInformation，即弹出如下对话框此对话框列出了ENVI中添加的地理坐标信息。 可以调整ENVI使用的Magic Pixel相对应的影像坐标。 ENVI可以从相应的头文件信息和地图投影文件中识别出地图的投影、像元大小以及地图投影参数，所以可以用它计算出影像中任意像元的地理坐标。 且可以输入地图坐标或者地理坐标（经纬度）。 （4）点击Projection/Datum文本旁边的切换按钮，可显示UTM Zone13North地图投影的经纬度坐标，如上对话框中所示。 （5）点击当前的DMS或者DDEG按钮，分别在度-分-秒（Degrees-Minutes-Seconds）和十进制的度（Decimal Degrees）之间相互切换，如 （6）点击Cancel，退出Edit MapInformation对话框。 （7）在主影像显示窗口中显示SPOT影像图，选择ToolsCursor Location/Value，即出现此对话框同时给出了像素坐标和地理坐标。 当光标在影像中移动时，对话框中便会显示特定位置的坐标值。 在此可清晰观察到，图像（即像素）坐标的变化一般很大，而经纬度一般只在分秒处变化。 （3.2）叠合地图公里网 （1）在主影像显示窗口菜单栏选择OverlayGrid Lines。 弹出#1Grid LineParameters对话框，同时，在影像中中出现一个虚拟的边框。 （2）在#1Grid LineParameters对话框中选择FileRestore Setup。 在弹出的Enter GridParameters对话框中，选择liuchunyan_sp.grd文件，点击Open。 即将先前保存的公里网参数加载到对话框中。 （3）在#1Grid LineParameters对话框的菜单栏中选择OptionsEdit MapGrid Attributes，可以查看地图参数，打开Edit MapAttributes对话框。 （4）Edit MapAttributes对话框中，查看公里网的间隔以及控制线条、标签、公里网夹角和矩形框相应的颜色及其他特征参数。 处理完后点击Cancel关闭。 （5）在#1Grid LineParameters对话框的菜单栏中选择OptionsEdit GeographicGrid Atrributes，查看地理坐标。 在弹出的对话框中，注意地理坐标（经纬度）公里网的参数。 处理过后，点击Cancel关闭。 （6）在#1Grid LineParameters对话框中，点击Apply，在主影像显示窗口中同时显示公里网图 （2）加载公里网 （7）在主影像窗口中选择OverlayAnnotation。 在弹出的#1AnnotationText对话框，选择FileRestore Annotation，在Enter AnnotationFilename中打开liuchunyan_sp.ann文件，点击OK。 即将先前保存的地图注记加载到影像图 （3）地理坐标注记（3.3）影像保存（3.3.1）ENVI影像格式保存 （1）在主影像显示窗口中，选择FileSave ImageAsImageFile。 弹出Output Displayto ImageFile对话框，默认Output FileType格式为ENVI。 （2）选择保存路径，文件名lcy_sp.img点击OK，输出影像。 （3）新生成的影像文件自动出现在Available BandsList对话框中，在此对话框中，点击Display#1下拉中选择New Display，打开一个新的显示窗口。 （4）选择RGB Color单选按钮，将影像加载到显示窗口，连续选择R、G、B带地理坐标的SPOT数据波段。 点击Load RGB，则添加注记后的影像作为一幅栅格图显示出来图 （4）带地理坐标的SPOT数据（3.3.2）PostScipt文件保存 （1）在主影像窗口中，选择FileSave ImageAsPostscipt File。 在Output Displayto PostsciptFile对话框中，注记和公里网都将自动地列在制图选项中。 （2）在弹出的对话框中作相应的设置点击OK。 即将影像图保存为Postscipt文件图 （5）Postscript文件根据上图发现，原影像图被保存了两页的postscipt文件。 在ENVI中，如果缩放操作时影像超过了可用页的大小，ENVI便会自动地创建多页Postscipt文件。 （3.4）TM影像图 （1）在ENVI主菜单中选择FileOpen ImageFile，在Enter DataFilenames中选中liuchunyan_tm.img文件，打开。 （2）击中Band3，在一个新的显示窗口中显示Band3影像图。 在此主影像显示窗口菜单中选择ToolsCursor Locations/Value。 在主影像窗口中移动光标，会发现Cursor Locations/Value对话框中显示的是像素坐标图 （6）TM影像图说明此影像图是基于像素坐标的，它不同于上面带有地理坐标的SPOT影像。 （3.5）影像到影像的配准 （1）此实验把带有地理坐标的SPOT影像作为基准影像，将一个基于像素坐标的Landsat TM影像进行校正，以匹配该SPOT影像。 在ENVI主菜单栏中，选择MapRegistrationSelect GCPsImage to Image。 在弹出的Image toImage Registration对话框中，点击并选择Display#1（SPOT影像），作为Base Image；点击Display#2（TM影像），作为Warp Image。 （2）点击OK，启动配准程序。 通过将光标放置在两幅影像的相同地物点上，添加单独的地面控制点。 （3）在Ground ControlPoints Selection对话框的Base X和Y文本框中，分别输入753和826，按键盘上的Enter键将SOPT影像中的光标移到相应的点上；在Warp X和Y文本框中，分别输入331和433，按键盘上的Enter键将TM影像中的光标移到相应的点上图 （7）地面控制点选择 （4）在Ground ControlPoints Selection对话框中，点击Add Point，把该地面控制点添加到列表中，点击show List查看地面控制点列表。 （5）当地面控制点选取至少有4个时，RMS误差便显示出来了。 当RMS误差较大时，删除一些比较不符的点，重新选择控制点，或者更新控制点使选点得以比较好的分布，这样便可以一定程度的减小RMS值图 （8）地面控制点此地面控制点相对应的情况为图 （9）地面控制点信息此表中均方根误差RMS得到一定的降低。 校正的精度与地面控制点的精度、分布、数量及校正范围有关。 （6）在Ground ControlPoints Selection对话框中，选择OptionsWarp DisplayBand。 弹出Registration Parameters对话框，在Warp Method按钮下选择RST；在Resampling的下拉中选择Nearest Neighbor重采样法，文件保存为liuchunyan_tm1wrp。 （7）弹出Registration Parameters对话框，在Warp Method按钮下选择RST；在Resampling的下拉中分别选择Bilinear和Cubic Convolution重采样法，文件分别保存为liuchunyantm2wrp和liuchunyantm3wrp。 （8）再一次重复上述步骤，但Warp Method选择Polynomial，并使用Cubic Convolution重采样方法，以及选Delaunay三角网的Triangulation校正法，相应的使用Cubic Convolution重采样方法。 分别保存为liuchunyan_tm4.wrp和liuchunyan_tm5.wrp文件。 得到五个wrp文件影像图图（10.1）RSTNearest Neighbor图（10.2）RSTBilinear图（10.3）RSTCubic Convolution图（10.4）PolynomialCubic Convolution图（10.5）TriangulationCubic Convolution （9）在Available BandsList中点击原始TM波段影像名liuchunyan_tm.img，然后从菜单栏中选择FileClose SelectedFile。 （10）在Available BandsList中选择liuchunyan_tm1.wrp文件。 在Display#下拉中选择New Display，点击Load Band。 （11）在主影像窗口中点击鼠标右键，选择Link Display。 在弹出的Link Display对话框中点击OK，把SPOT影像和已添加了地理坐标的TM影像链接起来。 （12）在主影像中点击鼠标左键，使用动态链接功能，对SPOT影像和TM影像进行比较。 （13）将liuchunyan_tm2.wrp和liuchunyan_tm3.wrp影像加载到新的显示窗口中，同上法创建动态链接功能观察产生的效果图 （11）动态链接比较通过观察上面各影像图（此为校正方法相同但重采样方法不同的影像），发现在使用最近邻法重采样的影像中像素具有锯齿状，使用双线性内插法重采样的影像看起来更加平滑，使用三次卷积法采样的影像既有平滑效果，同时也保持了影像的细节特征。 （14）在相应的主影像显示窗口中，关闭liuchunyan_tm1.wrp和liuchunyan_tm2.wrp影像窗口。 （15）将liuchunyan_tm4.wrp和liuchunyan_tm5.wrp文件在不同显示窗口中打开，并使用动态链接功能，同liuchunyan_tm3.wrp进行比较图 （12）动态链接比较通过观察上面影像图（它们的校正方法不同但重采样方法相同），会发现使用1次多项式和Delaunay三角网的校正方法均使图像有些过于平滑，致使细节有所失真。 （16）通过在主影像图中打开Cursor Location/Value对话框，查看它们的地理坐标，容易发现不同的重采样和校正方法，它们在同一点处的坐标值有一些差别，说明地理校正存在一些误差，而方法的不同使其产生的误差大小也有所区别。 （3.6）影像到地图的配准将带有地理坐标的SPOT影像中获取地图坐标以及一个矢量的数字线划图（DLG）作为基准数据，对基于像素坐标的Landsat TM影像进行校正，以匹配相应的地图数据。 （1）在ENVI主菜单中选择FileOpen ImageFile，在Enter DataFilenames对话框中选择liuchunyan_tm.img文件，打开。 （2）在Available BandsList中点击Gray Scale按钮，选择Band3，点击Load band。 （3）在E