《遥感原理与应用 》课程实习报告.doc

《《遥感原理与应用 》课程实习报告.doc》由会员分享，可在线阅读，更多相关《《遥感原理与应用 》课程实习报告.doc（52页珍藏版）》请在三一办公上搜索。

1、遥感原理与应用 课程实习实验一 高光谱数据分析一， 实验目的高光谱分辨率遥感用很窄而连续的光谱通道对礼物持续遥感城乡的技术。在可见光 到短波红外波段其光谱分辨率高达纳米数量级，通常具有波段多的特点，本部分实习将完成以下的内容：波谱库的操作、浏览和提取影像的反射率波谱，ENVI中感兴趣区域（ROI）的定义以及惊醒彩色合成影像的选取，其目的是为了鉴别波谱类型。二， 实验数据介绍Cuprite地区ATREM校正后的反射率数据，50个波段：lichuanjun cup95_at.intENVI相应的头文件：lichuanjun cup95_at.hdrENVI格式的JPL波谱库：jpl1.sliENV

2、I相应的头文件：jpl1.hdrENVI格式的USGS波谱库：usgs_min.sliENVI相应的头文件：rsga_min.gdr三， 实验过程在主菜单中选择File-Open Image File,选择lichuanjun cup95_at.int文件，点击Open1. 浏览影像波谱同波谱库进行比较在可用波段对话框中选择Band193(2.2008um),点击Load Band 将灰度影像加载到显示窗口中选择Tool-Profiles-Z Profile(Seqpectrum),提取表现反射率波普曲线：在ENVI菜单中选择Spectral-Spectral Library-Spectral

3、 Library Viewer按下面选择：在Spectral Library Input File对话框中点击Open- Spectral Library，从Spec_lib/jpl_sli子目录中选择jpl1.sli波谱库文件，点击“OK”，选择Select Input File区域中的jpl_sli在Spectral Library Viewer中选择Options-Edit(x,y)Scale Factors,并在Y DataMultiplier文本框中输入值1.000，点击“OK”在Spectral Library Viewer中选择下列波普名称，绘制他们的波普曲线：ALUNITE S

4、O-4ABUDDINGTONITE FELDS TS-11ACALCITE C-3DKAOLINITE WELL ORDERED PS-1A得到如下波普曲线绘制图：从绘制（plot）窗口中选择Edit-Plot Parameters,自定义波普曲线的绘制图。显示波普名称的图例，可以在绘制窗口中点击右键，从弹出的快捷菜单中选择Plot Key在绘制窗口中右键选择Stack plots 分别查看绘制的波普曲线，绘制的波普曲线如下：绘制的波普曲线：在#1Spectral Profile 窗口中点击鼠标右键选择Plot Key将显示出波普曲线的图例，该图例指出了所选像素的X、Y的像素坐标：在#1Spe

5、ctral Profile中选择Options-New Window:Blank，打开一个新的绘图窗口，然后重新放置#1Spectral Profile绘图窗口和新的绘图窗口的位置，使能够同时看到这两个窗口。在主图像窗口中选择Tool-Pixel Locator,使用Pixel Locator对话框定位到下列各点精确的像素位置上：列：590 435 494 531 502 448 260行：570 555 514 541 589 505 613在Pixel Locator对话框中输入像素的坐标列590、行570，使缩放指示矩形框移动到以这个像素为中心影像地区，然后点击Apply，将矩形框移动到

6、这个位置上同时#1Spectral Profile绘制窗口中将更新显示所选点的波谱曲线，其所对应的图例为：“X:590Y:570”。在新的绘图窗口中，点击鼠标右键选择Plot Key，打开显示了X、Y坐标位置的图例。#1Spectral Profile绘制窗口中用鼠标左键点击并按住图例“X:590Y:570”，将这个波谱曲线图例拖到新的绘图窗口中。对所列没一个像素点重复上面步骤，直到新的绘图窗口中包含了所有 7种波谱曲线。在新的绘图窗口中选择Options-Stack Plots,新的绘图窗口如下：2. 鉴别波谱曲线（1）在主菜单中选择Spectral-Spectral Analyst（2）在

7、Spectral Analyst Input Spectral Library打开并选择usgs_min.sli波谱库文件，在Edit Identify Methods weighting 对话框中点击OK（3）在主影像窗口中选择Tool-Profiles-Z Profile(Seqpectrum)，然后在#1Spectral Profile绘制窗口中点击右键，选择Plot Key（4）在主影像窗口中选择Tool- Pixel Locator,输入列502，行589，点Apply（5）在Spectral Analyst对话框中选择Options-Edit Methods weights, 在E

8、dit Identify Methods weighting对话框中为每一个Weight文本框输入0.33然后点击OK（6）在Spectral Analyst对话框中点击Apply,选择像素列502行589，则在Spectral Analyst对话框中出现下图：注意到列表中第一行显示kaolinite的波普曲线评分最高，这个相对高的分数值表明了该元素对应的地物与kaolinite最相似。（7）双击列表中的第一条波谱曲线，绘制出未知地物与波谱库中的波谱曲线，如下图：3. 定义感兴区（1）在主影像窗口中选择Overlay-Region of Interest 打开ROI Tool对话框。（2）在影

9、像中点击鼠标左键，绘制多边形，点击右键封闭该多边形，完成感兴区的定义，再次单击右键锁定感性区的位置，点击中键取消感性区。点击Stats计算感性区的数据：（3）按下图选择Cancel删除所选感性区4.二维散点图（1） 查看二维散点图在主影像窗口中选择Tools-2-D Scatter plots ,打开Scatter plot对话框在X、Y列表中选波段band193和band207，点击OK，下面的散点图显示了X、Y表现反射率值的关系（2） 散点图的密度分割（3）散点图中的跳跃像素从主影像窗口中点击左键并拖动鼠标，触发散点图中的跳跃像素。选择Options-Set Patch Size 改变散点

10、图中红色矩形光标的大小，并观察改变后的差异。（4）将散点图连接到波谱剖面廓线选择Options-Z Profile，选择一个输入文件点击OK将鼠标放到散点图中点击右键提取具有该散点位置特征的相应空间像素点的波谱曲线。（5）散点图中的感性区在 散点图中点击左键绘制感性区多边形，点击右键封闭多边形，在主窗口会以彩色编码显示出来：要移除某个类在散点图窗口中点击右键选择 Clear Class（6）散点图与波谱混合选择Options-Change Bands修改两个波段组合四， 实验体会通过这次实验，我学会了波谱库的操作、浏览和提取影像的反射率波谱，ENVI中感兴趣区域（ROI）的定义，惊醒彩色合成影

11、像的选取以及如何从感兴趣区中提取波谱信息，然后还将进行彩色合成，并使用二维散点图进行简单的分类。实验二 影像镶嵌一， 实验目的通过上机操作，掌握遥感图像镶嵌和裁切的基本方法、步骤、原理、并得到成果。影像拼接时将多付影像拼接成一幅图像，拼接成功后还要对其做进一步的平滑处理。二， 实验数据介绍基于象元的影像镶嵌： lichuanjundv06_2.img 、lichuanjundv06_3.img基于地理坐标的影像镶嵌：lichuanjun lch_01w.img、lichuanjun lch_02w.img三， 实验过程1. 基于像素的影像镶嵌File-Open Image File-lichu

12、anjundv06_2.img(和lichuanjundv06_3.img)在主菜单中选择Map-Mosaicking-Fixel Based在出现的Pixed Based Mosaic 窗口中选择Import -Import File在出现的Mosaic Input File 窗口中同时选择lichuanjun dvo6_2.img和lichuanjun dvo6_3.img在Select Mosaic Size 对话框中的Xsize输入614，Ysize输入1024，指定镶嵌影像的大小，点击“OK”在出现的Pixel Mosaic 窗口中点击lichuanjun dv06_3.img,在Y

13、O中输入513，按Enter键，这样lichuanjun dv06_3.img就出现在lichuanjun dv06_2.img下面了在Pixel Mosaic中选择File-Apply,在Mosaic Parameters中点击memory，点“OK生成镶嵌影像文件：在Pixel Mosaic中选择File-Save Template，在出现的Output osaic Template中输入输出的文件名：lichuanjun dv06.mos,点“OK”生成虚拟镶嵌影像文件：边缘羽化的处理操作：在Pixel Mosaic窗口中选择在Select Mosaic Size 对话框中的Xsize输

14、入614，Ysize输入1024，指定镶嵌影像的大小，点击“OK”在Pixel Mosaic窗口中右键单击影像的绿色轮廓框选择Edit Entry输入如下：点击“OK”对另一幅影像做相同操作选择File-Apply,输入255，点“OK”生成羽化后的影像2. 基于地理坐标的镶嵌将有地理坐标的影像嵌在一起时，常常需要进行羽化处理。专题的这一部分将展示如何使用ENVI的镶嵌工具，进行羽化处理，并创建基于地理坐标的镶嵌影像。（1）在ENVI主菜单中执行，将lichuanjun lch01w.img和lichuanjun lch02w.img文件加载到显示窗口中（2）在ENVI主菜单中，执行，开始进行

15、ENVI基于地理坐标的镶嵌。（3）执行，将weijinman_lch01w.img 和weijinman_lch02w.img文件加载到显示窗口中来。点击OK这时，我们可以看到两张图片的相对位置已经是正确的，但是两张影响有重叠部分。下面对其进行羽化处理。(4)右键单击绿色框中的任意部位，点击，出现如下的对话框。在Data Value to Ignore 文本框中输入数值0，在Feathering Distance 文本框中输入数值20，点击OK。 这时，我们可以看到原来的两张影像的拼接部分已经不是很明显了，羽化过程就此完成，基于坐标的影响镶嵌也就此完成。 四， 实验体会通过本次上机操作，我掌握

16、来遥感图像镶嵌和裁切的基本方法、步骤、原理、并得到成果。影像拼接时将多付影像拼接成一幅图像，拼接成功后还要对其做进一步的平滑处理。实验三 影像配准一， 实验目的本专题旨在介绍如何在ENVI 中对影像进行地理校正，添加地理坐标，以及如何使用ENVI 进行影像到影像的配准和影像到地图的校正。本专题介绍了使用ENVI 生成影像地图的步骤，并举例演示说明了全色影像和多光谱影像进行HSV 融合的步骤。二， 实验数据介绍lichuanjun bldr_sp.img Boulder SPOT 带地理坐标的影像子集lichuanjun bldr_sp.hdr ENVI 对应的头文件lichuanjun bld

17、r_sp.grd Boulder SPOT 地理公里网参数lichuanjun bldr_sp.ann Boulder SPOT 地图注记lichuanjun bldr_tm.img Boulder TM 没有地理坐标的影像lichuanjun bldr_tm.hdr ENVI 相应的头文件三， 实验过程A. 带地理坐标的数据和影像地图1.要打开带地理坐标的 SPOT 数据： 从 ENVI 主菜单中，选择 File Open Image File。 从列表中选择 lichuanjun bldr_sp.img 文件。 点击 OK。 当可用波段列表对话框出现后，点击 Gray Scale 单选按钮

18、，使用鼠标左键，点击相应的波段名，从对话框顶部所列波段中选中 SPOT 波段。所选择的波段名显示在 Selected Band:字段区域中。 点击 Load Band 按钮，加载这幅影像到一个新的显示窗口中。2.修改 ENVI 头文件中的地图信息 (1)在可用波段列表中，右键点击 bldr_sp.img 文件名下的 Map Info 图标，从弹出的快捷菜单中选择 Edit Map Information。 Edit Map Information对话框出现在屏幕上。这个对话框列出了在 ENVI 中添加地理坐标所用的地理信息。可以调整 ENVI 使用的 Magic Pixel（作为地图坐标系统的

19、起始点）相对应的影像坐标。因为 ENVI 可以从相应头文件信息和地图投影文件中，识别出地图投影、像元大小以及地图投影参数，所以用它能够计算出影像中任意像元的地理坐标。既可以输入地图坐标，也可以输入地理坐标（纬度/经度）。 (2) 点击 Projection/Datum文本旁边的箭头切换按钮，显示 UTM Zone 13 North地图投影的纬度/经度坐标。ENVI 在处理过程中才进行转换。 (3) 点击当前的 DMS 或者 DDEG 按钮，分别在度-分-秒（Degrees-Minutes- Seconds）和十进制的度（Decimal Degrees）之间进行切换。 (4) 点击 Cancel

20、，退出 Edit Map Information 对话框。3.光标位置/值(1)从主影像窗口菜单栏中，选择 Tools Cursor Location/Value。也可以从 ENVI 主菜单和主影像窗口菜单栏中，选择 Window Cursor Location/Value，打开这个对话框。注意到对于这幅带地理坐标的影像，对话框同时给出了像素坐标和地理坐标。 (2) 在影像中移动光标，查看特定位置的坐标值，并注意地图坐标和经纬度之间的关系。3. 叠合地图公里网(1)从主影像窗口菜单栏中，选择 Overlay Grid Lines。#1 Grid Line Parameters 对话框出现在屏幕

21、上，同时一个虚拟的边框添加到影像中，允许在影像外部显示地图公里网的标注。(2) 在这个新的对话框中，选择 File Restore Setup。 (3) 在 Enter Grid Parameters Filename 对话框中，选中 bldr_sp.grd 文件，点击 Open。先前保存过的公里网参数就会被加载到对话框中。 (4) 在#1 Grid Line Parameters 对话框中，从对话框菜单栏中，选择 Options Edit Map Grid Attributes，来查看地图参数。这将打开 Edit Map Attributes 对话框。 (5) 在 Edit Map Attr

22、ibutes 对话框中，注意公里网的间隔以及控制线条、标签、公里网交角以及矩形框（轮廓框）相应颜色和其它特征的参数。 (6) 处理完成后，点击 Cancel，关闭该对话框。 (7) 在#1 Grid Line Parameters 对话框，从菜单栏中，选择 Options Edit Geographic Grid Attributes，查看地理坐标。在 Edit Grid Attributes对话框中，再次注意地理坐标（纬度/经度）公里网的参数。处理完成后，点击 Cancel 来关闭该对话框。 (8) 在 Grid Line Parameters dialog 对话框中，点击 Apply，在影

23、像中放置公里网。 ENVI 允许同时放置像素、地图、和地理坐标公里网。4. 叠合地图注记(1) 在主影像窗口中，选择 Overlay Annotation。 (2) 在#1 Annotation: Text对话框中，选择 File Restore Annotation。打开一个标准的文件选择对话框。 (3) 在 Enter Annotation Filename 对话框中，从文件列表中选择 bldr_sp.ann 文件，点击 OK。先前保存过的地图注记被加载到影像上。(4) 按住滚动窗口的一角，并拖动鼠标，拉大该滚动窗口。重新放置改变了大小的滚动窗口，这样就可以同时看到主影像窗口。 (5) 在

24、改变了大小的滚动窗口中，使用鼠标左键，移动主影像指示矩形框，查看主影像窗口中出现的地图要素。 (6) 在#1 Annotation: Text对话框中，点击并按住 Object菜单，查看可以用来注记地图的对象。5. 输出到影像或 Postscript 文件将影像保存为 ENVI 的影像格式:(1) 在主影像窗口中，选择 File Save Image As Image File。 (2) 在 Output Display to Image File对话框中，选择 Output File Type 的下拉式按钮（默认情况下文件类型设置为 ENVI），查看可用的不同格式。 Change Graph

25、ics Overlay Selections 按钮可以打开一个同样名字的对话框，这个对话框允许添加或删除许多制图叠合选项（graphics options），包括注记和公里网。 Change Image Border Size 按钮也可以打开一个同名的对话框。这个对话框允许改变顶部、底部、左边和右边的边框宽度，如果需要，也可以改变边框的颜色。 如果带注记和公里网的彩色影像已经显示在显示窗口中，那么注记和公里网都将自动地列在制图叠合选项中。同时也可以选择其它要叠置在输出影像上的注记文件。 (3) 通过选择 Memory 或者 File 单选按钮，决定是将结果保存为一个磁盘文件，还是保存到内存中。

26、选择 Memory，点击 OK，输出影像。 (4) 新生成的影像文件自动列在可用波段列表中。在可用波段列表中，点击 Display #1 下拉式按钮，从菜单中选择 New Display，打开一个新的显示窗口。 (5) 选择 RGB Color 单选按钮，要将影像从内存中加载到显示窗口，可以连续选择 R、G 和 B（带地理坐标的 SPOT 数据）波段。 (6) 点击 Load RGB 按钮，添加注记后的影像作为一幅栅格图显示出来。将影像保存为 Postscript 文件:(1) 在主影像窗口中，选择 File Save Image As Postscript File。在 Output Dis

27、play to Postscript File 对话框中，注记和公里网都将自动地列在制图选项中。一个表述输出页的图形出现在对话框顶部靠右的地方。 (2) 在 xsize 和 ysize 参数文本框中，输入所需的输出影像大小。用鼠标左键点击对话框中代表输出页的图形，可以看到新影像的轮廓大小及其位置。 (3) 在代表输出页的图形上点击鼠标右键，把影像放置在输出页的中部。 如果想要缩放输出的地图，在 Map Scale 文本框中输入所需的地图比例，然后在代表输出页的图形上点击鼠标左键，查看结果。 如果缩放操作使影像超过了可用页的大小，那么 ENVI 会自动地创建多页 Postscript 文件。（4

28、） 如果有一个大比例的绘图仪，那么将 Page 页大小改为绘图大小，缩放的影像将会输出到一个Postscript 文件，该文件能按比例直接绘制在绘图仪上。 注意：只有在能打印彩色输出，并且能使用操作系统标准程序来打印 Postscript 输出文件的前提下，才能创建 Postscript 文件。 （5） 点击 OK，保存 postscript打印设置。如果不能打印彩色输出，点击 Cancel 取消输出操作。B.影像到影像的配准:1. 开并显示 Landsat TM 影像文件(1) 从 ENVI 主菜单中，选择 File Open Image File。 (2) 当 Enter Data Fil

29、enames 对话框出现后，选择进入 envidata 目录下的 bldr_reg 子目录，从列表中选择 lichuanjun bldr_tm.img 文件。 (3) 在文件选择对话框中，点击 Open，把 TM 影像波段加载到可用波段列表中。 (4) 在列表中选中波段 3 ，点击 No Display 按钮，并从下拉式菜单中选择 New Display。 (5)点击 Load Band 按钮，来把 TM 第3 波段的影像加载到一个新的显示窗口中。2. 显示光标位置/值(1) 从主影像窗口菜单栏中，选择 Tools Cursor Location/Value。 (2) 在主影像窗口、滚动窗口和

30、缩放窗口的 TM 影像上，移动鼠标光标。 注意到坐标是以像素单位给出的，这是因为这个影像是基于像素坐标的，它不同于上面带有地理坐标的 SPOT 影像。 (3) 选择 File Cancel，关闭 Cursor Location/Value对话框。4. 开始进行影像配准并加载地面控制点(1) 从 ENVI 主菜单栏中，选择 Map Registration Select GCPs: Image to Image。 (2) 在 Image to Image Registration对话框中，点击并选择 Display #1 （SPOT 影像），作为Base Image。点击 Display #2（

31、TM 影像），作为 Warp Image。 (3) 点击 OK，启动配准程序。通过将光标放置在两幅影像的相同地物点上，来添加单独的地面控制点。(4) 在Ground Control Points Selection对话框的Base X和Y文本框中，分别输入753和826，将SPOT影像中的光标移动到相应的点上。 (5) 使用同样的方法，在 Warp X 和 Y文本框中，分别输入 331和 433，将 TM 影像中的光标移动到相应的点上。 (6) 在两个缩放窗口中，查看光标点所处位置。如果需要，在每个缩放窗口所需位置上，点击鼠标左键，调整光标点所处的位置。 注意：在缩放窗口中支持亚像元（sub-

32、pixel）级的定位。缩放比例越大，地面控制点定位的精度就越好。 (7) 在 Ground Control Points Selection对话框中，点击 Add Point，把该地面控制点添加到列表中。点击 Show List 查看地面控制点列表。尝试选择几个地面控制点找到选择地面控制点的感觉。 注意对话框中所列的实际影像点和预测点坐标。一旦已经选择了至少 4个地面控制点以后，RMS误差就会显示出来。 (8) 在 Ground Control Points Selection 对话框中，选择 Options Clear All Points，可以清除掉所有已选择的地面控制点。 (9) 从 G

33、round Control Points Selection 对话框中，选择 File Restore GCPs from ASCII。 (10) 在 Enter Ground Control Points Filename 对话框中，选择文件 bldr_tm.pts，然后点击 OK，加载这个预先保存过的地面控制点坐标。5. 操作处理地面控制点在 Image to Image GCP List 对话框中，选择相应的地面控制点，然后在 Ground Control Points Selection 对话框中进行修改，这样可以编辑单个控制点的坐标位置。可以通过输入一个新的像素坐标，或使用对话框中的

34、方向箭头逐像素地移动坐标位置。 在 Image to Image GCP List对话框中，点击 On/Off 按钮，屏蔽掉所选择的地面控制点。这样在校正模型和 RMS 计算中都将不会考虑这个地面控制点坐标。这些地面控制点并没有被真正地删除，仅仅是被忽略掉了，可以使用 On/Off 按钮重新激活这些地面控制点。 在 Image to Image GCP List对话框中，点击 Delete 按钮，可以从列表中删除一个地面控制点。 在两个缩放窗口中调整光标位置，然后点击 Image to Image GCP List 对话框中的 Update 按钮，更新所选的地面控制点，将其修改到当前光标的所在

35、位置。 Image to Image GCP List 对话框中的 Predict按钮，允许对新的地面控制点进行预测。它以当前的校正模型为基础。 将包含 SPOT 影像的那个主影像窗口的光标放置到一个新的位置上。然后点击 Predict按钮，放置在 TM 影像上的光标就会根据校正模型，移动到预测的匹配点上去。 通过在 TM数据中，轻微地移动光标，能够对所提取的位置点进行交互式的精确定位。 在 Ground Control Points Selection 对话框中，点击 Add Point，把这个新的控制点添加到列表中。6. 校正影像(1) 从 Ground Control Points Se

36、lection 对话框中，选择 Options Warp Displayed Band。(2) 在 Registration Parameters 对话框中的 Warp Method按钮菜单中，选择 RST。在 Resampling 的按钮菜单中选择 Nearest Neighbor 重采样法。(3) 输入文件名 lichuanjun bldr_tm1.wrp，点击 OK。 (4) 重复步骤 1和步骤 2，还是使用 RST 校正法，但是要相应地选择 Bilinear 和 Cubic Convolution重采样法。 (5) 将结果分别输出到lichuanjun bldr_tm2.wrp 和li

37、chuanjun bldr_tm3.wrp 文件中。 (6) 再一次重复步骤 1 和步骤 2，这一次选择一次多项式 Polynomial 校正法，并使用 Cubic Convolution重采样法。然后再选择 Delaunay 三角网的 Triangulation 校正法，相应地使用 Cubic Convolution重采样法。 (7) 将结果分别输出到 bldr_tm4.wrp 和 bldr_tm5.wrp 文件中。7. 比较校正结果(1) 在可用波段列表中，点击原始的 TM 波段影像名 bldr_tm.img，然后从菜单栏中，选择 File Close Selected File。 (2)

38、 在随后出现的 ENVI 警告对话框中，点击 Yes关闭相应的影像文件。(3) 在可用波段列表中，选择 BLDRTM_1.WRP 文件。在 Display #下拉式按钮中选择 New Display，点击 Load Band将该文件加载到一个新的显示窗口中。 (4) 在主影像窗口中，点击鼠标右键，选择 Tools Link Link Displays。 (5) 在 Link Displays 对话框中，点击 OK，把 SPOT 影像和已添加了地理坐标的 TM 影像链接起来。 (6) 在主影像显示窗口中，点击鼠标左键，使用动态链接功能，对 SPOT 影像和 TM影像进行比较。 (7) 将 bld

39、r_tm2.wrp 和 bldr_tm3.wrp 影像加载到新的显示窗口中，使用影像动态链接功能，比较采用三种不同的重采样法（最临近法、双线性内插法和三次卷积法）所产生的效果。 注意观察，在使用最近邻法重采样的影像中的锯齿状像素，而使用双线性内插法重采样的影像看起来更加平滑，使用三次卷积法重采样的影像是最好的结果，不但有平滑效果，而且保持了影像的细节特征。 (8) 在相应的主影像窗口中，选择 File Cancel，关闭 bldr_tm1.wrp（RST 校正，最近邻法重采样）和 bldr_tm2.wrp（RST 校正，双线性内插法重采样）影像的显示窗口。 (9) 将 bldr_tm4.wrp

40、 和 bldr_tm5.wrp 影像加载到新的显示窗口中，使用影像动态链接功能，同bldr_tm3.wrp 影像（RST 校正）进行比较。C.影像到地图的配准1. 打开并显示 Landsat TM 影像文件(1) 从 ENVI 主菜单中，选择 File Open Image File。(2) 当 Enter Data Filenames对话框出现后，选择进入 envidata 目录下的 bldr_reg 子目录，从列表中选择 lichuanjun bldr_tm.img 文件。 (3) 点击 OK。 TM 影像波段被加载到可用波段列表中，同时一幅彩色影像被加载到一个新的显示窗口中。 (4) 在

41、可用波段列表中，点击 Gray Scale 按钮，选择波段 3。 (5) 点击 Load Band 按钮，把 TM 影像的第 3 波段加载到已打开的显示窗口中。2. 选择影像到地图的配准并恢复控制点坐标(1) 从 ENVI 主菜单中，选择 Map Registration Select GCPs:Image to Map。 (2) 如果打开了多个影像显示窗口，那么就在 Image to Map Registration 对话框中，点击选择包含该灰阶影像的那个显示窗口的显示号。(3) 从投影列表中选择 UTM，并在 Zone 文本框中输入 13。 (4) 设置像素大小为 30m，点击 OK，启动

42、配准程序。 (5) 在要校正的影像中，把光标移动到一个已知地图坐标的地面点上（可以从一幅地图或者 ENVI 矢量文件中见下一部分读取所需的地图坐标），来添加单个的地面控制点。 (6) Ground Control Points Selection 对话框中的 E（东向）和 N（北向）文本框中，手动地输入已知的地图坐标，然后点击 Add Point 来添加新的地面控制点。 (7) 在 Ground Control Points Selection 对话框中，选择 File Restore GCPs from ASCII，打开bldrtm_m.pts 文件。(8) 在Ground Control

43、Points Selection对话框中，点击Show List按钮。可以在Image to Map GCP List对话框中查看影像的地图坐标、实际影像点和预测点的坐标以及 RMS误差。3. 使用矢量显示的数字线划图（DLGs）来添加地图控制点(1) 从 ENVI 主菜单中，选择 File Open Vector File USGS DLG。 (2) 在文件选择对话框中，选择 bldr_rd.dlg 文件。 (3) 在 Import Optional DLG File Parameters 对话框中，选择 Memory 单选按钮，点击 OK，读入所需的数字线划图（DLG）数据。 (4) 在可

44、用矢量列表中高亮选择 ROADS AND TRAILS: BOULDER, CO 文件，点击 Load Selected按钮。 (5) 在 Load Vector 对话框中，点击 New Vector Window。把该矢量加载到一个新的矢量显示窗口中。(6) 在 Vector Window #1 窗口中，点击并拖曳鼠标左键，激活一个十字形光标。 光标处的地图坐标会在 Vector Window #1 窗口的底部列出。 (7) 在主影像显示窗口中，选择 Tools Pixel Locator，并输入 402 和418，然后点击Apply，将影像光标移动到道路交叉口相应的点上去。 注意，在缩放窗

45、口中，同样可以获取到亚像元（sub-pixel）级的定位精度。 (8) 在矢量窗口中，用鼠标左键点击并拖曳矢量光标，当十字形光标位于所需的道路交叉口时，松开鼠标左键，把矢量光标放置在道路的交叉口上，其坐标为 477593.74，4433240.0（北纬 40 度 3 分3 秒，东经-105度 15 分45 秒）。 (9) 在矢量窗口中，点击鼠标右键，并从弹出的快捷菜单中，选择Export Map Location。新的地图坐标就会出现在 Ground Control Points Selection 对话框中。 (10) 在 Ground Control Points Selection 对话框中，点击 Add Point，添加该地图坐标/影像像素对，并观察 RMS误差的变化。4.RST和三次卷积校正(1) 在 Ground Control Points Selection 对话框中，选择 Options Warp File。 (2) 在 Input Warp Image 对话框中，高亮选择文件名 bldr_tm.img，点击 OK，对 TM的 6 个波段进行校正。 (3) 在