1、实验过程 启动 ENVI 并加载AVIRIS 影像数据启动前，请确认已正确安装 ENVI。 要在 UNIX 或Macintosh OS X 中启动ENVI，请在UNIX 命令行中输入envi。 要在 Windows 系统中启动ENVI，请双击ENVI 的图标。当程序成功地加载并执行后，ENVI 的主菜单将会出现在屏幕上。1. 在 ENVI 主菜单中，选择File Open Image File2. 选择 cup95_at.int 文件作为输入文件名，点击Open 弹出可用波段列表，它将列出50 个波段的名字。 显示灰阶影像1. 在可用波段列表对话框中，选择 Band 193
2、（2.2008um）2. 点击 Gray Scale 单选按钮，然后点击Load Band。将灰度影像加载到显示窗口中。3. 从主影像窗口菜单中选择 Tools Profiles Z Profile （Spectrum），提取表观反射率波谱曲线。 浏览影像波谱并同波谱库进行比较1. 在影像上移动缩放指示矩形框，同时查看#1 Spectral Profile 窗口中的波谱曲线，浏览整个影像的表观反射率波谱曲线。2. 在主影像窗口中，使用鼠标左键点击并拖动缩放指示矩形框或者直接点击鼠标左键，将缩放指示矩形框移动到以所选像素点为中心的区域中。3. 将从影像中获取的表观反射率波谱曲线同所选波谱库中的波
3、谱曲线进行比较。ENVI 提供了几个不同的波谱库，根据本专题的目的，我们将会使用JPL 波谱库（Groves 等，1992）以及USGS 波谱库（Clarke 等，1993）。4. 从 ENVI 主菜单中选择Spectral Spectral Libraries Spectral Library Viewer。5. 在 Spectral Library Input File 对话框中，点击Open File 按钮，从spec_lib/jpl_lib 子目录中，选择jpl1.sli 波谱库文件，点击OK。6. 选择 Select Input File 区域中的jpl1.sli，点击OK。7. 在
4、 Spectral Library Viewer 对话框中，选择Options Edit （x, y） Scale Factors，并在Y DataMultiplier 文本框中，输入值1000，以匹配影像表观反射率范围（1-1000），点击OK。8. 在 Spectral Library Viewer 对话框中，选择下列波谱名称，绘制它们的波谱曲线：选中ALUNITE SO-4A 、BUDDINGTONITE FELDS TS-11A、CALCITE C-3D、KAOLINITE WELL ORDERED PS-1A后波谱曲线绘制图如下：9. 从绘制（plot）窗口菜单中，选择Edit Pl
5、ot Parameters，自定义波谱曲线的绘制图。在PlotParameters 对话框中，按下面的步骤进行： 将 Charsize 减少为0.50。 选择 X-Axis 单选按钮，然后将Range 调整为1.90 到2.45。 还是选中 X-Axis 单选按钮，点击Left/Right Margins 的箭头增量按钮，直到达到所需X 方向的页边距。 选择 Y-Axis 单选按钮，将Axis Title 改为“Reflectance”。 还是选中 Y-Axis 单选按钮，点击Top/Bottom Margins 的箭头增量按钮，直到所需的Y 方向的页边距。 点击 Apply，然后再点击Can
6、cel。10. 要显示波谱名称的图例，可以在绘制窗口中点击鼠标右键，从弹出的快捷菜单中选择Plot Key。将绘制窗口拖动到所需的大小，以容纳下波谱名称。11. 在绘制窗口中，选择 Options Stack Plots，分别查看绘制的波谱曲线。绘制的波谱曲线如下图所示：12. 在#1 Spectral Profile 绘图窗口中点击鼠标右键，从弹出的快捷菜单中选择Plot Key，这将显示出波谱曲线的图例，该图例指出了所选像素的X 和Y 的像素坐标。13. 在#1 Spectral Profile 绘图窗口中，选择Options New Window: Blank，打开一个新的绘图窗口。然后
7、重新放置#1 Spectral Profile 绘图窗口和新绘图窗口的位置，使能够同时看到这两个绘图窗口。14. 从主影像窗口菜单中，选择Tools Pixel Locator。使用Pixel Locator 对话框，定位到下列各点精确的像素位置上：15. 在 Pixel Locator 对话框中，输入像素的坐标，列（sample）590、行（line）570，使缩放指示矩形框移动到以这个像素为中心的影像地区，即Stonewall Playa 地区，然后点击Apply，将矩形框移动到这个位置上。同时#1 Spectral Profile 绘图窗口将更新显示所选点的波谱曲线，其所对应的图例为：“
8、X：590 Y：570”。16. 在新的绘图窗口中，点击鼠标右键，在弹出的快捷菜单中，选择 Plot Key，打开显示了X 和Y 坐标位置的图例。17. 在#1 Spectral Profile 绘图窗口中，使用鼠标左键，点击并按住图例“X：570”，将这个波谱曲线图例拖到新的绘图窗口中。18. 对上表所列的每一个像素点重复上面的步骤，直到新的绘图窗口包含了所有的 7 种波谱曲线。19. 在新的绘图窗口中，选择 Options Stack Plots。新的绘图窗口如下图所示：将波谱曲线和波谱库波谱曲线比较，效果图如下： 鉴别波谱曲线（使用Spectral Analyst TM）：使用 Spec
9、tral AnalystTM 来鉴别波谱曲线：ENVI 提供了一个波谱匹配工具，它根据波谱库中的波谱曲线对影像中的波谱曲线进行评分。波谱分析使用多种方法产生一个在0 到1 之间的分数值，其中分数值1相当于完全匹配。1. 从 ENVI 主菜单中，选择Spectral Spectral Analyst。2. 点击 Spectral Analyst Input Spectral Library 对话框底部的Open Spec Lib 按钮。3. 选择进入 usgs_min 波谱库目录，选择usgs_min.sli 波谱库文件，点击Open。4. usgs_min.sli 文件出现在Spectral
10、Analyst Input Spectral Library 对话框中，选中该文件，点击OK。5. 在 Edit Identify Methods Weighting 对话框中，点击OK。6. 从主影像窗口菜单中，选择Tools Profiles Z Profile （Spectrum）。然后在#1 Spectral Profile绘图窗口中，点击鼠标右键，从弹出的快捷菜单中，选择Plot Key，显示波谱曲线名称的图例。7. 从主影像窗口菜单中，选择 Tools Pixel Locator。8. 在 Pixel Locator 对话框中，输入像素点的坐标列502，行589，点击Apply。9
11、. 在 Spectral Analyst 对话框中，选择Options Edit Method Weights。10. 在 Edit Identify Methods Weighting 对话框中，为每一个Weight 文本框输入值0.33，然后点击OK。不同的匹配方法在ENVI 遥感影像处理用户手册（ENVI Users Guide）中进行描述。11. 在 Spectral Analyst 对话框中，点击Apply。如果在#1 Spectral Profile 绘图窗口中显示了多条波谱曲线，那么将会出现一个波谱曲线列表。如果出现了该波谱曲线列表，那么就选择像素（列502，行589）所对应的那
12、条波谱曲线。Spectral Analyst 对话框将如下图所示：下面显示出了矿物质kaolinite 的波谱曲线与像素（列502，行589）的波谱曲线最匹配双击列表中的第一条曲线，在同一绘图窗口中，绘制出未知地物的波谱曲线及波谱库中的波谱曲线，进行比较。该绘制图如下图所示：可以得出，未知地物与波谱库中最匹配的kaolinite的波谱曲线两者高度匹配。 定义感兴趣区创建新的感兴趣区1. 在影像中，点击鼠标左键。2. 在多边形顶点处，点击鼠标左键，绘制感兴趣区，或者点击并按住鼠标左键，移动鼠标，连续绘制感兴趣区。3. 点击鼠标右键，封闭该多边形，完成感兴趣区的定义。再次点击鼠标右键锁定感兴趣区的
13、位置。4. 点击 Stats 按钮，计算感兴趣区的统计信息，绘制均值波谱曲线（白色），均值波谱上下各有一条的标准差曲线（绿色），以及最小和最大值的包络波谱曲线（红色），其包含了感兴趣区所有的波谱。计算感兴趣区的统计信息：加载预先保存的感兴趣区文件：从感兴趣区中提取均值波谱曲线：在同一个窗口中，绘制每个感兴趣区的均值波谱曲线。并将每条波谱带分别显示出来，进行比较： 鉴别矿物质设计选择合适的彩色合成来鉴别矿物质：1. 在可用波段列表对话框中，选择 RGB 单选按钮，并顺次点击Band 183、Band 193 和Band 207。2. 点击 Load RGB，将彩色合成影像加载到当前影像显示窗口中
14、。3. 在主影像窗口中，点击Tools Profiles Z-Profile （Spectrum）。注意到在Z 剖面廓线（Z-Profile）窗口中，被用来合成彩色影像的三波段的位置上，分别用红、绿、蓝三条垂直线标出。4. 在 ROI Tool 对话框中，选择Off 单选按钮，使用Z 剖面廓线窗口，在主影像窗口中的感兴趣区上或其附近浏览相应的波谱曲线。5. 通过使用鼠标左键，点击并拖动波谱剖面廓线窗口中的颜色条（plot bars），将其拖动到所需的波段位置上。在可用波段列表对话框中，选择对话框，并用183、193、207波段组合显示：在Z剖面廓线窗口中，被用来合成彩色影像的三波段位置上，分别
15、用RGB三条垂线标出：突出影像上特定地物的一种方法是将某个彩色条放置在某个特征吸收中心，而把另两个彩色条放置在相对的波谱峰值上。根据先前显示的均值波谱的波谱特性，选择RGB波段对应的位置。实验结果实验结果，在实验过程中已经阐述出来，这里就不再累赘叙述。实验体会经过这次实验让我了对地物光谱反射率，有了深刻的了解，知道了假定在 RGB 影像上已知特定像素的颜色，基本可以预测它所对应的波谱曲线是什么样子的。也可以根据训练样区的波谱特征，解释它们的颜色。测试特定的 RGB 波段选择，识别出某种矿物质。实验二 、使用envi进行影像镶嵌影像镶嵌是一门艺术，它将把多幅影像连接合并，以生成一幅单一的合成影像
16、。ENVI 提供了交互式的方式来将没有地理坐标（no-georeferenced）的影像拼接在一起，或者自动地拼接有地理坐标的影像，并输出成有地理坐标的镶嵌影像。镶嵌程序提供了透明处理、直方图匹配，以及颜色自动平衡的功能。在遥感影像处理与分析中，预处理是最初的也是最基本的影像操作。图像校正是从具有畸变的图像中消除畸变的处理过程，消除集合畸变的称为几何校正。几何校正是从具有几何畸变的图像中消除畸变的过程。其任务是定量地确定图像上的像元坐标与目标物的地理坐标的对应关系。校正的最终目的是确定校正后图像的行列数，然后找到新图像中每一像元的亮度值。Pixel-Based Mosaicking（基于像素的
17、影像镶嵌）KJLdv06_2.img AVIRIS 02 景影像KJLdv06_2.hdr ENVI 相应的头文件KJLdv06_3.img AVIRIS 03 景影像KJLdv06_3.hdr ENVI 相应的头文件KJLdv06a.mos AVIRIS 拼接影像镶嵌模板文件KJLdv06b.mos 羽化后的AVIRIS 影像镶嵌模板文件KJLdv06_fea.img 羽化后的镶嵌影像KJLdv06_fea.hdr ENVI 相应的头文件Georeferenced Mosaicking（基于地理坐标的影像镶嵌）KJLlch_01w.img 直方图匹配校正后的影像KJLlch_01w.hdr
18、ENVI 相应的头文件KJLlch_01w.ann 切割线的注记文件KJLlch_02w.img 直方图匹配校正后的影像KJLlch_02w.hdr ENVI 相应的头文件KJLlch_a.mos 带地理坐标的影像镶嵌模板文件KJLlch_mos1.img 带地理坐标的影像镶嵌结果KJLlch_mos1.hdr ENVI 相应的头文件置基于像素的影像在 ENVI 主菜单中，选择Map Mosaicking Pixel Based，开始进行ENVI 基于像素的镶嵌操作。Pixel Based Mosaic 对话框出现在屏幕上。1. 从 Pixel Based Mosaic 对话框中，选择Impo
19、rt Import Files。2. 在 Mosaic Input Files 对话框中，点击Open File，选择进入avmosaic 目录，选择文件KJLdv06_2.img。3. 在 Mosaic Input Files 对话框中，再一次点击Open File，选择KJLdv06_3.img 文件。4. 在 Mosaic Input Files 对话框中，按下键盘上的Shift 键，并同时点击KJLdv06_2.img 和KJLdv06_3.img文件名，选中这两个文件，点击OK。5. 在 Select Mosaic Size 对话框的X Size 中输入614，Y Size 中输入1
20、024，指定镶嵌影像的大小。6. 在 Pixel Based Mosaic 对话框中，点击KJLdv06_3.img 文件名。影像当前的位置就会以像素值为单位，列在对话框底部的文本框中。7. 在 YO 文本框中，输入值513，按下键盘上的Enter 键。这样KJLdv06_3.img 影像就会放置在KJLdv06_2.img 影像的下面。 8. 在 Pixel Based Mosaic 对话框中，选择File Apply。当Mosaic Parameters 对话框出现后，输入输出文件名KJLdv06.img，点击OK，生成镶嵌影像文件。 要生成虚拟镶嵌影像，而不是新的镶嵌影像文件，在 Pix
21、el Based Mosaic 对话框中，选择File Save Template。当Output Mosaic Template 对话框出现后，输入输出的文件名KJLdv06a.mos。9. 点击可用波段列表中的 KLJdv06a.mos 波段名，然后点击Load Band，显示镶嵌后的影像。第二部分展示了在拼接的镶嵌影像中布置两幅影像位置的方法，既可以输入XO 和YO 的值，也可以在对话框中把影像拖曳到所需位置。这个例子还包括了边缘羽化的处理操作。1. 在 Pixel Based Mosaic 对话框中，选择OptionsChange Mosaic Size。在Select Mosaic
22、Size 对话框的X Size 和Y Size 文本框中都输入值768，点击OK，改变输出镶嵌影像的大小。2. 在 Pixel Based Mosaic 对话框中，左键点击影像2 的绿色轮廓框。将影像2 拖动到镶嵌图的右下角。3. 在镶嵌图中，右键点击影像1 的红色轮廓框，选择Edit Entry，打开Entry:filename 对话框。4. 在 Data Value to Ignore 文本框中，输入值0。在Feathering Distance 文本框中，输入值25，点击OK。5. 对另一幅影像，重复上面的两步操作。6. 选择 File Save Template，输入输出文件名KJLd
23、v06b.mos。在可用波段列表中，点击镶嵌模板文件名，然后点击Load Band，显示该镶嵌影像。当使用虚拟镶嵌时，不会进行羽化处理。7. 现在真实地创建输出文件，同时对镶嵌影像进行羽化处理。在 Pixel Based Mosaic 对话框中，选择File Apply，点击OK。镶嵌后的影像如图：基于地理坐标的影像镶嵌将有地理坐标的影像镶嵌在一起时，常常需要进行羽化处理。专题的这一部分将展示如何使用ENVI的镶嵌工具，进行羽化处理，并创建基于地理坐标的镶嵌影像。创建基于地理坐标的镶嵌影像1. 在 ENVI 主菜单中，选择Map Mosaicking Georeferenced，开始进行ENV
24、I 基于地理坐标的镶嵌操作。2. 在 Map Based Mosaic 对话框中选择Import然后选择文件KJLlch_01w.img，KJLlch_02w.img打开文件，设置基于地理坐标并进行羽化镶嵌时所必需的参数。查看顶部影像，切割线和虚拟未羽化的镶嵌影像1. 在可用波段列表中，选择KJLlch_01w.img 文件。点击波段名，然后点击Load Band 按钮，以灰阶的方式显示该影像。2. 在主影像窗口中，点击右键弹出快捷菜单，选择Toggle Display Scroll Bars，开启滚动条。点击水平滚动条，滚动影像，直到显示出影像合适的部分。3. 从主影像窗口中，选择 Over
25、lay Annotation，打开Annotation 对话框。4. 在 Annotation 对话框中，选择File Restore Annotation，然后选择KJLlch_01w.ann 文件。这将显示出一条红色的切割线（cutline），该切割线用来在镶嵌影像中混合两影像。5. 从可用波段列表中，将KJLlch_02w.img 加载到一个新的显示窗口中。查看该影像上切割线的特性。选择并将KJLlch_01w.img显示创建输出羽化后的镶嵌影像1. 在 Map Based Mosaic 对话框中，选择File Apply。在Mosaic Parameters 对话框中，输入输出文件名K
26、JLlch_mos.img，点击OK，创建羽化后的镶嵌影像。2. 关闭包含单独的倾斜影像的两个影像显示窗口，并将镶嵌影像加载到一个新的显示窗口中。3. 使用影像动态链接功能，对羽化后的镶嵌影像和没有羽化的镶嵌影像进行比较。在下图中，左边的影像为倾斜的直方图匹配后的影像，并且已选定了切割线。右边的影像为使用切割线羽化后的镶嵌影像。镶嵌结果如下：KJLbldr_tm1.wrp 使用缩放平移和最近邻重采样法得到的影像到影像的配准结果。KJLbldr_tm1.hdr ENVI 相应的头文件。KJLbldr_tm2.wrp 使用RST 和双线性内插重采样法进行的影像到影像的配准结果。KJLbldr_tm
27、2.hdr ENVI 相应的头文件。KJLbldr_tm3.wrp 使用RST 和三次卷积重采样法进行的影像到影像的配准结果。KJLbldr_tm3.hdr ENVI 相应的头文件。KJLbldr_tm4.wrp 使用一次多项式和三次卷积重采样法进行的影像到影像的配准结果。KJLbldr_tm4.hdr ENVI 相应的头文件。KJLbldr_tm5.wrp 使用Delaunay 三角网和三次卷积重采样法进行的影像到影像配准结果KJLbldr_tm5.hdr ENVI H 相应的头文件。KJLbldrtm_m.imgBoulder TM 影像到地图的配准结果，使用了RST 和三次卷积重采样法。
28、KJLbldrtm_m.hdr ENVI 相应的头文件。实验三、影像地理坐标定位和配准本实验旨在练习如何在ENVI 中对影像进行地理校正，添加地理坐标，以及如何使用ENVI 进行影像到影像的配准和影像到地图的校正。ENVI 的影像配准和几何纠正工具允许用户将基于像素的影像定位到地理坐标上，然后对它们进行几何纠正，使其匹配基准影像的几何信息。使用全分辨率（主影像窗口）和缩放窗口来选择地面控制点（GCPs），进行影像到影像和影像到地图的配准。基准影像和未校正影像的控制点坐标都会显示出来，同时由指定的校正算法所得的误差也会显示出来。地面控制点预测功能能够使对地面控制点的选取简单化。将使用重采样、缩放比例和平移（这三种方法通称RST），以及多项式函数（多项式系数可以从1 到n），或者Delaunay 三角网的方法，来对影像进行校正。所支持的重采样方法包括最近邻法（nearest-neighbor）、双线性内插法（bilinear interpol
