1.1. 专题概述
本专题旨在介绍如何在ENVI中对影像进行地理校正,添加地理坐标,以及如何使用ENVI进行影像到影像的配准和影像到地图的校正。本专题介绍了使用ENVI生成影像地图的步骤,并举例演示说明了全色影像和多光谱影像进行HSV融合的步骤。在开始本专题内容前,我们假定用户已经熟悉了一般影像配准和重采样的概念。完成本专题内容需要1到2个小时。
♦ 本专题中使用的文件
光盘:《ENVI遥感影像处理专题与实践》附带光盘 #1 路径:envidata/bldr_reg
文件 所需的文件
bldr_sp.img bldr_sp.hdr bldr_sp.grd bldr_sp.ann bldr_tm.img bldr_tm.hdr bldr_tm.pts
描述
Boulder SPOT 带地理坐标的影像子集 ENVI对应的头文件
Boulder SPOT 地理公里网参数 Boulder SPOT 地图注记
Boulder TM 没有地理坐标的影像 ENVI相应的头文件
TM-SPOT 影像到影像配准中所用的控制点
bldrtm_m.pts TM-Map 影像到地图配准中所用的控制点 bldr_rd.dlg
Boulder 道路数字线划图(DLG)
bldrtmsp.grd 融合后的TM-SPOT影像的地图公里网 bldrtmsp.ann 融合后的TM-SPOT影像的注记
生成的文件
bldr_tm1.wrp 使用缩放平移和最近邻重采样法得到的影像到影像的配准结果 bldr_tm1.hdr ENVI 相应的头文件
bldr_tm2.wrp 使用RST和双线性内插重采样法进行的影像到影像的配准结果 bldr_tm2.hdr ENVI 相应的头文件
bldr_tm3.wrp 使用RST和三次卷积重采样法进行的影像到影像的配准结果 bldr_tm3.hdr ENVI 相应的头文件
bldr_tm4.wrp 使用一次多项式和三次卷积重采样法进行的影像到影像的配准结果
- 64 -
bldr_tm4.hdr ENVI 相应的头文件
bldr_tm5.wrp 使用Delaunay三角网和三次卷积重采样法进行的影像到影像配准结果bldr_tm5.hdr ENVI H相应的头文件 bldrtm_m.img Boulder TM影像到地图的配准结果,使用了RST和三次卷积重采样
法
bldrtm_m.hdr ENVI 相应的头文件
bldrtmsp.img Boulder TM/SPOT 使用HSV融合后的结果,分辨率为10米 bldrtmsp.hdr ENVI 相应的头文件
1.2. ENVI中带地理坐标的影像
ENVI对带地理坐标的影像提供了全面的支持,它能够对许多预定义的地图投影进行处理,这些地图投影可以采用UTM或State Plane投影方式。此外,ENVI的用户自定义地图投影功能能够创建自定义的地图投影,它允许使用6种基本投影类型,超过35种的不同椭球体以及100多种的基准数据集(Datum),来满足大多数地图投影的需要。
ENVI地图投影参数存储在一个ASCII文本文件map_proj.txt中,该文本文件能够被ENVI地图投影工具修改,或者直接被用户编辑。这个文件中的信息会被影像相应的头文件(ENVI Header files)所使用,而且ENVI允许使用已知的地图投影坐标来简单地指定相关联的Magic Pixel(地图坐标系统的起始点)。然后,选择的ENVI函数就能够使用该信息,在带地理坐标的数据空间中进行操作处理。
ENVI的影像配准和几何纠正工具允许用户将基于像素的影像定位到地理坐标上,然后对它们进行几何纠正,使其匹配基准影像的几何信息。使用全分辨率(主影像窗口)和缩放窗口来选择地面控制点(GCPs),进行影像到影像和影像到地图的配准。基准影像和未校正影像的控制点坐标都会显示出来,同时由指定的校正算法所得的误差也会显示出来。地面控制点预测功能能够使对地面控制点的选取简单化。
将使用重采样、缩放比例和平移(这三种方法通称RST),以及多项式函数(多项式系数可以从1到n),或者Delaunay三角网的方法,来对影像进行校正。所支持的重采样方法包括最近邻法(nearest-neighbor)、双线性内插法(bilinear interpolation)和三次卷积法(cubic convolution)。使用ENVI的多重动态链接显示功能对基准影像和校正后的影像进行比较,可以快速地评估配准的精度。
以下部分提供了一些ENVI自带的基于地理坐标的处理功能的例子。请参见《ENVI遥感影像处理实用手册》(ENVI User’s Guide)来获取更多更多的信息。
1.3. 带地理坐标的数据和影像地图
本专题的这一部分将使你熟悉ENVI中对带地理坐标的数据的处理,使用地图公里网和注记创建影像地图,并生成输出影像。
- 65 -
♦ 启动ENVI
启动前,请确认已正确安装ENVI。
z z
要在UNIX或Macintosh OS X中启动ENVI,请在UNIX命令行中输入envi。 要在Windows系统中启动ENVI,请双击ENVI的图标。
♦ 打开并显示SPOT数据
要打开带地理坐标的SPOT数据:
1. 从ENVI主菜单中,选择File → Open Image File。
2. 当Enter Data Filename文件选择对话框出现后,选择进入envidata目录下的bldr_reg子目
录,从列表中选择bldr_sp.img文件。 3. 点击OK。
4. 当可用波段列表对话框出现后,点击Gray Scale单选按钮,使用鼠标左键,点击相应的波段名,
从对话框顶部所列波段中选中SPOT波段。所选择的波段名显示在Selected Band:字段区域中。 5. 点击Load Band按钮,加载这幅影像到一个新的显示窗口中。
♦ 修改ENVI头文件中的地图信息
1. 在可用波段列表中,右键点击bldr_sp.img文件名下的Map Info图标,从弹出的快捷菜单中选
择Edit Map Information。
Edit Map Information对话框出现在屏幕上。
图5-1:Edit Map Information对话框
- 66 -
这个对话框列出了在ENVI中添加地理坐标所用的地理信息。可以调整ENVI使用的Magic Pixel(作为地图坐标系统的起始点)相对应的影像坐标。因为ENVI可以从相应头文件信息和地图投影文件中,识别出地图投影、像元大小以及地图投影参数,所以用它能够计算出影像中任意像元的地理坐标。既可以输入地图坐标,也可以输入地理坐标(纬度/经度)。
2. 点击Projection/Datum文本旁边的箭头切换按钮,显示UTM Zone 13 North地图投影的纬度/经度坐
标。ENVI在处理过程中才进行转换。
3. 点击当前的DMS或者DDEG按钮,分别在度-分-秒(Degrees-Minutes- Seconds)和十进制的度(Decimal Degrees)之间进行切换。
4. 点击Cancel,退出Edit Map Information对话框。
♦ 光标位置/值
要打开一个显示主影像窗口,滚动窗口,或者缩放窗口中光标位置信息的对话框,可以按如下几步进行。 1. 从主影像窗口菜单栏中,选择Tools → Cursor Location/Value。也可以从ENVI主菜单和主影像
窗口菜单栏中,选择Window → Cursor Location/Value,打开这个对话框。
图5-2:Cursor Location对话框,显示了像素值和相应的地理坐标(对带地理坐标的影像而言)
注意到对于这幅带地理坐标的影像,对话框同时给出了像素坐标和地理坐标。
2. 在影像中移动光标,查看特定位置的坐标值,并注意地图坐标和经纬度之间的关系。
♦ 叠合地图公里网
1. 从主影像窗口菜单栏中,选择Overlay → Grid Lines。#1 Grid Line Parameters对话框出现在屏
幕上,同时一个虚拟的边框添加到影像中,允许在影像外部显示地图公里网的标注。
- 67 -
图5-3:Grid Line Parameters对话框
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 Attributes对话框中,注意公里网的间隔以及控制线条、标签、公里网交角以及矩形
框(轮廓框)相应颜色和其它特征的参数。 6. 处理完成后,点击Cancel,关闭该对话框。
7. 在#1 Grid Line Parameters对话框,从菜单栏中,选择Options → Edit Geographic Grid
Attributes,查看地理坐标。在Edit Grid Attributes对话框中,再次注意地理坐标(纬度/经度)公里网的参数。处理完成后,点击Cancel来关闭该对话框。
8. 在Grid Line Parameters dialog对话框中,点击Apply,在影像中放置公里网。 ENVI允许同时放置像素、地图、和地理坐标公里网。
叠合地图注记
1. 在主影像窗口中,选择Overlay → Annotation。
2. 在#1 Annotation: Text对话框中,选择File → Restore Annotation。打开一个标准的文件选择对话
框。
- 68 -
3. 在Enter Annotation Filename对话框中,从文件列表中选择bldr_sp.ann文件,点击OK。先
前保存过的地图注记被加载到影像上。
4. 按住滚动窗口的一角,并拖动鼠标,拉大该滚动窗口。重新放置改变了大小的滚动窗口,这样就
可以同时看到主影像窗口。
5. 在改变了大小的滚动窗口中,使用鼠标左键,移动主影像指示矩形框,查看主影像窗口中出现的
地图要素。
6. 在#1 Annotation: Text对话框中,点击并按住Object菜单,查看可以用来注记地图的对象。
♦ 输出到影像或Postscript文件
ENVI对用户提供了几个选项,来保存和输出地图影像。可以将结果保存为ENVI的影像文件格式,或者保存为几种通用的图形格式(包括Postscript格式)。也可以直接打印或者导入到其它软件中。
将影像保存为ENVI的影像格式
要将结果影像保为ENVI自身的格式(比如一个RGB文件),按下面的步骤进行操作。 1. 在主影像窗口中,选择File → Save Image As → Image File。
2. 在Output Display to Image File对话框中,选择Output File Type的下拉式按钮(默认情况下文件
类型设置为ENVI),查看可用的不同格式。
Change Graphics Overlay Selections按钮可以打开一个同样名字的对话框,这个对话框允许添加或删除许多制图叠合选项(graphics options),包括注记和公里网。
Change Image Border Size按钮也可以打开一个同名的对话框。这个对话框允许改变顶部、底部、左边和右边的边框宽度,如果需要,也可以改变边框的颜色。
如果带注记和公里网的彩色影像已经显示在显示窗口中,那么注记和公里网都将自动地列在制图叠合选项中。同时也可以选择其它要叠置在输出影像上的注记文件。
3. 通过选择Memory或者File单选按钮,决定是将结果保存为一个磁盘文件,还是保存到内存中。
选择Memory,点击OK,输出影像。
4. 新生成的影像文件自动列在可用波段列表中。在可用波段列表中,点击Display #1下拉式按钮,从
菜单中选择New Display,打开一个新的显示窗口。
5. 选择RGB Color单选按钮,要将影像从内存中加载到显示窗口,可以连续选择R、G和B(带地
理坐标的SPOT数据)波段。
6. 点击Load RGB按钮,添加注记后的影像作为一幅栅格图显示出来。
将影像保存为Postscript文件
要将结果影像保存为Postscript文件,按下面的步骤进行操作。
- 69 -
1. 在主影像窗口中,选择File → Save Image As → Postscript File。在Output Display to Postscript
File对话框中,注记和公里网都将自动地列在制图选项中。一个表述输出页的图形出现在对话框顶部靠右的地方。
2. 在xsize和ysize参数文本框中,输入所需的输出影像大小。用鼠标左键点击对话框中代表输出页
的图形,可以看到新影像的轮廓大小及其位置。
3. 在代表输出页的图形上点击鼠标右键,把影像放置在输出页的中部。
z
如果想要缩放输出的地图,在Map Scale文本框中输入所需的地图比例,然后在代表输出页的图形上点击鼠标左键,查看结果。
如果缩放操作使影像超过了可用页的大小,那么ENVI会自动地创建多页Postscript文件。
图5-4:输出影像地图的Output Display to Postscript对话框
4. 如果有一个大比例的绘图仪,那么将Page页大小改为绘图大小,缩放的影像将会输出到一个
Postscript文件,该文件能按比例直接绘制在绘图仪上。
注意:只有在能打印彩色输出,并且能使用操作系统标准程序来打印Postscript输出文件的前提下,才能创建Postscript文件。
5. 点击OK,保存postscript打印设置。如果不能打印彩色输出,点击Cancel取消输出操作。
- 70 -
♦ 直接进行打印
ENVI也能够直接打印输出到你的操作系统所支持的设备上。
1. 从主影像窗口中,选择File → Print。操作系统的标准打印对话框出现后,就可以运行标准的打
印程序了。
2. 一旦已经选定了操作系统标准打印对话框中的所有参数,点击OK。ENVI将打开一个Output
Display to Printer对话框,允许设定额外的基本ENVI打印参数,这和postscript输出程序中所用的参数(见上)比较类似。将这些打印设置调整到所需的值,点击OK开始进行打印。
1.4. 影像到影像的配准
本专题的这一部分将逐步演示影像到影像的配准处理过程。带有地理坐标的SPOT影像被用作基准影像,一个基于像素坐标的Landsat TM影像将被进行校正,以匹配该SPOT影像。
♦ 打开并显示Landsat TM影像文件
1. 从ENVI主菜单中,选择File → Open Image File。
2. 当Enter Data Filenames对话框出现后,选择进入envidata目录下的bldr_reg子目录,从列表
中选择bldr_tm.img文件。
3. 在文件选择对话框中,点击Open(在UNIX操作系统下为OK),把TM影像波段加载到可用波段
列表中。
4. 在列表中选中波段3 ,点击No Display按钮,并从下拉式菜单中选择New Display。 5. 点击Load Band按钮,来把TM第3波段的影像加载到一个新的显示窗口中。
♦ 显示光标位置/值
要打开一个显示主影像窗口,滚动窗口,或者缩放窗口中光标位置信息的对话框,可以按一下步骤进行操作。
1. 从主影像窗口菜单栏中,选择Tools → Cursor Location/Value。 2. 在主影像窗口、滚动窗口和缩放窗口的TM影像上,移动鼠标光标。
注意到坐标是以像素单位给出的,这是因为这个影像是基于像素坐标的,它不同于上面带有地理坐标的SPOT影像。
3. 选择File → Cancel,关闭Cursor Location/Value对话框。
♦ 开始进行影像配准并加载地面控制点
1. 从ENVI主菜单栏中,选择Map → Registration → Select GCPs: Image to Image。
2. 在Image to Image Registration对话框中,点击并选择Display #1 (SPOT影像),作为Base Image。
点击Display #2(TM影像),作为Warp Image。
- 71 -
3. 点击OK,启动配准程序。通过将光标放置在两幅影像的相同地物点上,来添加单独的地面控制点。
图5-5:用来进行影像到影像配准的Ground Control Points Selection对话框
4. 在Ground Control Points Selection对话框的Base X和Y文本框中,分别输入753和826,将SPOT
影像中的光标移动到相应的点上。
5. 使用同样的方法,在Warp X和Y文本框中,分别输入331和433,将TM影像中的光标移动到相
应的点上。
6. 在两个缩放窗口中,查看光标点所处位置。如果需要,在每个缩放窗口所需位置上,点击鼠标左
键,调整光标点所处的位置。
注意:在缩放窗口中支持亚像元(sub-pixel)级的定位。缩放比例越大,地面控制点定位的精度就越好。
7. 在Ground Control Points Selection对话框中,点击Add Point,把该地面控制点添加到列表中。
点击Show List查看地面控制点列表。尝试选择几个地面控制点找到选择地面控制点的感觉。 注意对话框中所列的实际影像点和预测点坐标。一旦已经选择了至少4个地面控制点以后,RMS误差就会显示出来。
8. 在Ground Control Points Selection对话框中,选择Options → Clear All Points,可以清除掉所有
已选择的地面控制点。
9. 从Ground Control Points Selection对话框中,选择File → Restore GCPs from ASCII。 10. 在Enter Ground Control Points Filename对话框中,选择文件bldr_tm.pts,然后点击OK,加
载这个预先保存过的地面控制点坐标。
11. 在Image to Image GCP List对话框中,点击单独的地面控制点。查看两幅影像中相应地面控制点
的位置、实际影像点和预测点的坐标以及RMS误差。调整对话框的大小,观察Ground Control Points Selection对话框中所列的合计RMS误差(RMS Error)。
- 72 -
图5-6:影像到影像配准中所用的Image to Image GCP LIst对话框
♦ 操作处理地面控制点
下面的内容仅提供处理方法,并且只对有限的地面控制点按钮的处理功能进行操作。
z
在Image to Image GCP List对话框中,选择相应的地面控制点,然后在Ground Control Points Selection对话框中进行修改,这样可以编辑单个控制点的坐标位置。可以通过输入一个新的像素坐标,或使用对话框中的方向箭头逐像素地移动坐标位置。
z
在Image to Image GCP List对话框中,点击On/Off按钮,屏蔽掉所选择的地面控制点。这样在校正模型和RMS计算中都将不会考虑这个地面控制点坐标。这些地面控制点并没有被真正地删除,仅仅是被忽略掉了,可以使用On/Off按钮重新激活这些地面控制点。
z z
在Image to Image GCP List对话框中,点击Delete按钮,可以从列表中删除一个地面控制点。 在两个缩放窗口中调整光标位置,然后点击Image to Image GCP List对话框中的Update按钮,更新所选的地面控制点,将其修改到当前光标的所在位置。
z
Image to Image GCP List对话框中的Predict按钮,允许对新的地面控制点进行预测。它以当前的校正模型为基础。
1. 将包含SPOT影像的那个主影像窗口的光标放置到一个新的位置上。然后点击Predict按钮,放置
在TM影像上的光标就会根据校正模型,移动到预测的匹配点上去。
2. 通过在TM数据中,轻微地移动光标,能够对所提取的位置点进行交互式的精确定位。
3. 在Ground Control Points Selection对话框中,点击Add Point,把这个新的控制点添加到列表中。
♦ 校正影像
我们可以校正显示的影像波段,也可以同时校正多波段影像中的所有波段。这里我们仅对已显示的波段进行校正。
1. 从Ground Control Points Selection对话框中,选择Options → Warp Displayed Band。
- 73 -
2. 在Registration Parameters对话框中的Warp Method按钮菜单中,选择RST。在Resampling的
按钮菜单中选择Nearest Neighbor重采样法。
图5-7:Registration Parameters对话框
3. 输入文件名bldr_tm1.wrp,点击OK。
4. 重复步骤1和步骤2,还是使用RST校正法,但是要相应地选择Bilinear和Cubic Convolution重
采样法。
5. 将结果分别输出到bldr_tm2.wrp和bldr_tm3.wrp文件中。
6. 再一次重复步骤1和步骤2,这一次选择一次多项式Polynomial校正法,并使用Cubic Convolution
重采样法。然后再选择Delaunay三角网的Triangulation校正法,相应地使用Cubic Convolution重采样法。
7. 将结果分别输出到bldr_tm4.wrp和bldr_tm5.wrp文件中。
♦ 比较校正结果
使用动态链接来比较校正结果:
1. 在可用波段列表中,点击原始的TM波段影像名bldr_tm.img,然后从菜单栏中,选择File →
Close Selected File。
2. 在随后出现的ENVI警告对话框中,点击Yes关闭相应的影像文件。
- 74 -
3. 在可用波段列表中,选择BLDRTM_1.WRP文件。在Display #下拉式按钮中选择New Display,点
击Load Band将该文件加载到一个新的显示窗口中。
4. 在主影像窗口中,点击鼠标右键,选择Tools → Link → Link Displays。
5. 在Link Displays对话框中,点击OK,把SPOT影像和已添加了地理坐标的TM影像链接起来。 6. 在主影像显示窗口中,点击鼠标左键,使用动态链接功能,对SPOT影像和TM影像进行比较。 7. 将bldr_tm2.wrp和bldr_tm3.wrp影像加载到新的显示窗口中,使用影像动态链接功能,比
较采用三种不同的重采样法(最临近法、双线性内插法和三次卷积法)所产生的效果。
注意观察,在使用最近邻法重采样的影像中的锯齿状像素,而使用双线性内插法重采样的影像看起来更加平滑,使用三次卷积法重采样的影像是最好的结果,不但有平滑效果,而且保持了影像的细节特征。
8. 在相应的主影像窗口中,选择File → Cancel,关闭bldr_tm1.wrp(RST校正,最近邻法重采
样)和bldr_tm2.wrp(RST校正,双线性内插法重采样)影像的显示窗口。
9. 将bldr_tm4.wrp和bldr_tm5.wrp影像加载到新的显示窗口中,使用影像动态链接功能,同
bldr_tm3.wrp影像(RST校正)进行比较。
注意观察采用三种不同校正方法(RST、1次多项式和Delaunay三角网)对影像几何信息所产生的效果。
10. 使用动态链接功能,与带有地理坐标的SPOT影像进行比较。
♦ 查看地图坐标
要打开Cursor Location/Value对话框:
1. 从主影像窗口菜单栏中,选择Tools → Cursor Location/Value。
2. 浏览带地理坐标的数据集,注意不同的重采样法和校正法对数据值所产生的效果。 3. 选择File → Cancel,关闭该对话框。
♦ 关闭所有文件
在ENVI主菜单中,通过选择File → Close All Files,可以关闭所有的数据文件。
1.5. 影像到地图的配准
本专题的这一部分将逐步地演示影像到地图的配准处理过程。许多步骤同影像到影像的配准步骤相似,因此这些步骤将不会被详细地讨论。从带地理坐标的SPOT影像中获取的地图坐标以及一个矢量的数字线划图(DLG)都将被作为基准数据,然后对基于像素坐标的Landsat TM影像进行校正,以匹配相应的地图数据。
♦ 打开并显示Landsat TM影像文件
1. 从ENVI主菜单中,选择File → Open Image File。
- 75 -
2. 当Enter Data Filenames对话框出现后,选择进入envidata目录下的bldr_reg子目录,从列
表中选择bldr_tm.img文件。 3. 点击OK。
TM影像波段被加载到可用波段列表中,同时一幅彩色影像被加载到一个新的显示窗口中。 4. 在可用波段列表中,点击Gray Scale按钮,选择波段3。
5. 点击Load Band按钮,把TM影像的第3波段加载到已打开的显示窗口中。
♦ 选择影像到地图的配准并恢复控制点坐标
1. 从ENVI主菜单中,选择Map → Registration → Select GCPs:Image to Map。
2. 如果打开了多个影像显示窗口,那么就在Image to Map Registration对话框中,点击选择包含该
灰阶影像的那个显示窗口的显示号。
3. 从投影列表中选择UTM,并在Zone文本框中输入13。 4. 设置像素大小为30m,点击OK,启动配准程序。
5. 在要校正的影像中,把光标移动到一个已知地图坐标的地面点上(可以从一幅地图或者ENVI矢量
文件中[见下一部分]读取所需的地图坐标),来添加单个的地面控制点。
6. Ground Control Points Selection对话框中的E(东向)和N(北向)文本框中,手动地输入已知
的地图坐标,然后点击Add Point来添加新的地面控制点。
7. 在Ground Control Points Selection对话框中,选择File → Restore GCPs from ASCII,打开
bldrtm_m.pts文件。
图5-8:用来进行影像到地图配准的Ground Control Points Selection对话框
8. 在Ground Control Points Selection对话框中,点击Show List按钮。可以在Image to Map GCP List
对话框中查看影像的地图坐标、实际影像点和预测点的坐标以及RMS误差。
- 76 -
图5-9:影像到地图配准中所用的Image to Map GCP LIst对话框
♦ 使用矢量显示的数字线划图(DLGs)来添加地图控制点
1. 从ENVI主菜单中,选择File → Open Vector File → USGS DLG。 2. 在文件选择对话框中,选择bldr_rd.dlg文件。
3. 在Import Optional DLG File Parameters对话框中,选择Memory单选按钮,点击OK,读入所
需的数字线划图(DLG)数据。
4. 在可用矢量列表中高亮选择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,将影像
光标移动到道路交叉口相应的点上去。
注意,在缩放窗口中,同样可以获取到亚像元(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误差的变化。
- 77 -
图5-10:带十字形光标的矢量窗口,并且显示出了地图坐标
♦ RST和三次卷积校正
1. 在Ground Control Points Selection对话框中,选择Options →Warp File。
2. 在Input Warp Image对话框中,高亮选择文件名bldr_tm.img,点击OK,对TM的6个波段
进行校正。
3. 在出现的Registration Parameters dialog对话框中,将Warp Method选为RST,将Resampling设
置为Cubic Convolution。 4. 把background值改为255。
5. 在output file文本框中,输入输出文件名bldrtm_m.img。 6. 点击OK,开始进行影像到地图的校正。
♦ 显示结果并进行评价
使用光标位置/值(Cursor Location/Value)来对校正后的彩色影像进行评价。
1. 在可用波段列表中,点击RGB单选按钮,接着点击校正影像的波段4、波段3和波段2(作为RGB)。 2. 从Display #下拉式菜单按钮中,选择New Display。点击Load RGB,来加载这幅校正后的TM彩
色影像。
- 78 -
注意到校正影像是倾斜的,这是由于消除了Landsat TM轨道方向影响的原因。此时这个影像已经带有地理坐标,但注意它的空间分辨率是30米,而SPOT影像为10米分辨率。
z
如果需要,将SPOT影像加载到一个新的显示窗口中,比较影像的几何信息和比例。
♦ 关闭所选的文件
保留bldrtm_m.img和bldr_sp.img影像,因为在下一部分会继续使用这些影像文件。
1. 在可用波段列表中,点击其它的文件名,然后选择File → Close Selected File,来关闭这些影像。 2. 在Vector Window Parameters #1和Vector Window #1对话框中,选择File → Cancel关闭这两个
窗口。
3. 在可用波段列表中,选择File → Cancel,来关闭该对话框。
4. 在Ground Control Points Selection对话框中,选择File → Cancel,关闭对话框。如果需要,保
存地面控制点。
1.6. 对不同分辨率的带地理坐标的数据集进行HSV融合
本专题的这一部分将介绍融合操作,它将对两幅不同分辨率的带地理坐标的数据集进行融合处理。我们将使用配准过的TM彩色合成影像作为低分辨率的多光谱影像,而带地理坐标的SPOT影像作为高分辨率的影像。融合后的结果为增强了空间分辨率的彩色合成影像。
♦ 显示30米分辨率的TM彩色合成影像
1. 如果已关闭了配准后的TM影像,那么重新打开文件bldrtm_m.img。
2. 点击可用波段列表中的RGB单选按钮,将波段4、波段3和波段2(分别对应R、G和B)加载到
一个新的显示窗口中。
♦ 显示10米分辨率的SPOT影像
1. 如果已关闭了SPOT影像,那么重新打开文件bldr_sp.img。
2. 点击可用波段列表中的Gray Scale按钮,然后点击Display #按钮,从下拉菜单中选择New Display。
点击Load Band按钮,将SPOT影像加载到一个新的显示窗口中。
将SPOT影像同TM影像进行比较,注意影像中相似的几何信息,以及不同的空间范围和影像比例。
♦ 进行HSV变换融合
1. 从ENVI主菜单中,选择Transform → Image Sharpening → HSV。
2. 如果己经加载了彩色影像,那么就从Select Input RGB对话框中选择合适的显示窗口。否则,在
Select Input RGB Input Bands对话框中,选择TM影像的波段4、波段3和波段2,然后点击OK。
- 79 -
3. 这样就打开了High Resolution Input File对话框。在Select Input Band列表中选择SPOT影像,
点击OK。
4. 在HSV Sharpening Parameters对话框中,输入输出文件名bldrtmsp.img,点击OK。 一个显示处理进度的状态条出现在屏幕上。当处理完成后,新生成的影像会自动出现在可用波段列表中。
♦ 显示10米分辨率的彩色影像
1. 在可用波段列表中,选择RGB Color单选按钮,然后从列出的新生成的文件中,选择R、G和B
波段,点击 Load RGB,将融合后的彩色影像加载到一个新的显示窗口中。 将HSV融合后的彩色影像同原始TM彩色合成影像以及SPOT影像进行比较。
2. 用标准化彩色变换(Color Normalized (Brovey) Transform),试着进行同样的处理。选择Transforms
→ Image Sharpening → Color Normalized (Brovey),并输入所需的文件信息,然后点击OK。
♦ 叠合地图公里网
1. 在HSV变换融合后的主影像显示窗口中,选择Overlay → Grid Lines。
在Grid Line Parameters对话框的同时,一个虚拟的边框也会添加到影像中,允许在影像外部显示地图公里网的标注。
2. 在Grid Line Parameters对话框中,选择File → Restore Setup。在随后打开的Enter Grid
Parameters Filename对话框中,选择bldrtmsp.grd文件,点击Open(在UNIX操作系统下为OK)。以前保存过的公里网参数被加载到对话框中。 3. 点击Apply,在影像中放置公里网。
♦ 叠合影像注记
1. 从HSV变换融合后的主影像显示窗口中,选择Overlay → Annotation。
2. 在相应的Annotation: Text对话框中,选择File → Restore Annotation,在文件列表中选择
bldrtmsp.ann文件,点击Open(在UNIX操作系统下为OK)。 将以前保存过的地图注记加载到影像上。
z
按住滚动窗口的一角,并拖动鼠标,拉大该滚动窗口。
♦ 输出影像地图
要保存结果影像,可以按照70页“输出到影像或Postscript文件”中所描述的步骤进行。可以:
z z z
创建一个影像输出文件 创建一个Postscript文件
打印一个影像地图的副本(见72页的“直接进行打印”)
- 80 -
♦ 结束ENVI程序
在ENVI主菜单中选择File → Exit(在UNIX操作系统下是(Quit),在弹出的Terminate this ENVI Session对话框中选择Yes,并点击OK,退出ENVI程序。如果使用的是ENVI RT,退出ENVI会返回操作系统。
- 81 -
因篇幅问题不能全部显示,请点此查看更多更全内容