遥感图处理实验指导
实验 图处理软件ENVI4功介绍
实验二 影配准
实验三 图融合
实验四 波段组合计算图增强
实验五 图分类
实验六 图制图
实验 图处理软件ENVI4功介绍
1 ENVI4支持文件格式
ENVI通栅格数格式包含简单二进制文件相关ASCII头文件读该文件格式允许ENVI影文件包括身嵌入头文件信息影文件
通栅格数会存储二进制字节流通常BSQ〔波段序〕BIP〔波段元交叉〕者BIL〔波段行交叉〕方式进行存储
ENVI软件支持数类型包括字节型整型符号整型长整型符号长整型浮点型双精度浮点型复数型双精度复数型64位整型符号64位整型
2 ENVI窗口显示
ENVI软件时翻开图通常三窗口影窗口〔Main Image Window〕滚动窗口〔Scroll Window〕缩放窗口〔Zoom Window〕见图11
图11 ENVI显示窗口组：影窗口滚动窗口缩放窗口
21 ENVI菜单
ENVI菜单位显示窗口顶部默认情况显示影窗口顶部提供交互式影显示分析功见图12果选显示组中包括影窗口菜单会出现滚动窗口者缩放窗口顶部ENVI菜单样该菜单中选择意菜单项选择项外显示窗口点击右键会弹出快捷菜单通菜单访问显示函数改变显示设置
图12影窗口中overlay菜单快捷菜单
22 波段列表
ENVI翻开影文件者文件中单波段波段列表〔Available Bands List〕特殊ENVI话框包括翻开文件中影波段相关图信息列表见图13波段列表彩色灰阶影加载显示窗口中直接翻开新显示窗口点击话框底部Display钮中选择显示窗口号翻开相应显示窗口然点击Gray Scale 者RGB单项选择钮列表中点击波段名选择需波段加载显示影
提示：果加载单波段影直接波段双击完成
图13波段列表话框
3 ENVI根功
翻开影文件
首先启动ENVI软件双击软件图标翻开影文件分步骤：
1 选择File→Open Image file屏幕弹出文件选择话框Enter Input Data File〞
2 选择进入实验课程提供数目录envidata中双击翻开文件夹选择can_tm文件夹文件夹中选择can_tmrimg文件然点击ok
3 波段列表中会显示刚刚选择文件can_tmrimg鼠标左键点击选择影显示方式〔Gray ScaleRGB〕选择某影波段三波段进行显示选波段标Selected Band〞显示出
4 点击Load BandLoad RGB影加载新显示窗口
图加载熟悉窗口界面图窗口滚动窗口缩放窗口进行相关操作
显示光标位置
显示鼠标光标位置值影窗口菜单栏者ENVI菜单栏中选择Window→Cursor LocationValue者影窗口中点击鼠标右键弹出快捷菜单中选择Cursor LocationValue
接着屏幕出现Cursor LocationValue话框显示出光标影窗口滚动窗口者缩放窗口中位置见图14该话框显示十字丝光标应象素屏幕颜色值实际数值
关闭话框Cursor LocationValue话框顶部菜单中选择file→cancel
旦Cursor LocationValue话框翻开隐藏显示该话框影窗口中双击鼠标左键
图14 Cursor LocationValue话框
显示影剖面廓线
交互式选择显示X轴〔水〕Y轴〔竖直〕Z轴〔波谱〕剖面廓线图剖面廓线图显示穿影横线X线Y者波谱波段Z数值
1 影窗口菜单栏中选择Tools→Profiles→X profle会翻开绘制窗口该窗口根影中选择行绘制出幅数值列号〔sample number〕间关系曲线图见图15
2 重复述程选择Y profile绘制出幅数值行号〔line number〕间关系曲线图选择Z profile绘制出相应波谱面廓线图
提示：通意影窗口中快捷菜单翻开Z轴波谱剖面廓线图
3 选择Window→Mouse Button Descriptions查鼠标键剖面廓线显示窗口中键功描述
4 放置三剖面廓线窗口便时
5 影移动移动十字丝查三幅影剖面廓线图新位置样显示新数值
6 绘制图窗口中选择File→Cancel关闭绘制窗口
图15水轴〔X〕Y轴〔竖直〕Z轴〔波谱〕剖面廓线绘制曲线图
进行快速拉伸
影窗口缩放窗口者滚动窗口中默认参数数进行快速度拉伸选择影窗口菜单栏中Enhance菜单进行种样度拉伸〔线性拉伸0－255间线性拉伸2％线性拉伸高斯拉伸均衡化拉伸方根拉伸〕
1 影窗口缩放窗口者滚动窗口中影作拉伸数源尝试进行种类型拉伸变换
2 拟窗口显示组中线性高斯均衡方根拉伸结果
显示交互式散点图
ENVI绘制出两选影数值关系图分选定两波段XY轴面坐标绘制两者散点图
1 影窗口菜单栏中选择Tools → 2D Scatter Plots接着Scatter Plot Band Choice话框会出现屏幕该话框中选择进行拟两影波段
2 选择中波段作X轴波段作Y轴然点击OK生成散点图
3 旦翻开散点图绘制窗口〔图16〕鼠标光标放影窗口中意位置住鼠标左键拖动光标时十字丝光标周围10×10范围素散点图中应点会红色突出显示出
提示：选择Window→Mouse Button Descriptions 显示鼠标键Scatter Plot显示窗口中功
4 影窗口移动鼠标光标观察产生跳跃素〔dancing pixels 〕效果
5 散点图突出显示影窗口中含应波谱值素鼠标光标放散点图窗口点击拖动鼠标中间键
6 散点图菜单栏中选择file→cancel 关闭Scatter Plot 窗口
图16 交互式散点图波段2波段3进行拟
加载幅彩色图
1果屏幕没显示波段列表话框ENVI菜单栏中选择Window →Available Bands list翻开该话框图13
注意：果没翻开影ENVI菜单中选择File→Open Image File然选择翻开影接着该影波段加载波段列表中
2 点击RGB单项选择钮然显示窗口中加载幅彩色影
3 通点击列表中波段名列表中种颜色〔红绿蓝〕选择波段点击列表中波段名时指定RGB单项选择钮会动前选择
4 三种颜色波段应时点击Display＃1拉式菜单钮中选择New Display
5 点击Load RGB钮幅彩色影加载新显示窗口中
链接两显示窗口
两显示窗口链接起进行拟两显示窗口链接起显示窗口中〔滚动缩放等〕进行操作会相链接显示窗口中产生相响应两显示窗口中链接起操作步骤：
1 影窗口菜单中选择Tools→link→Link Displays者影中点击鼠标右键弹出快捷菜单中选择Link DisplaysLink Displays话框会翻开
2 Link Displays话框中点击ok建立链接
3 现尝试中窗口进行滚动者缩放操作观察窗口反响
选择感兴趣区域〔Region of Interest〕
ENVI允许影中定义感兴趣区〔region of interestROIs〕感兴趣区提取分类信息统计生成掩膜操作
1 影窗口中选择Overlay→Region of interest者影中点击鼠标右键弹出菜单话框中选择ROI Tool接着该影显示窗口相应ROI Tool话框会出现屏幕〔图17〕
2 绘制感兴趣区边形ROI_Type中选择ROI_Type类型包括边形矩形圆形等具体绘制步骤：
影窗口中点击鼠标左键建立感兴趣边形第顶点
次点击鼠标左键序选择更边线点点击鼠标右键闭合该边形鼠标中键删新定义点者删整边形次点击鼠标右键固定边形位置
通选择ROI Tool话框顶部相应单项选择钮感兴趣区缩放窗口中定义
完成感兴趣区定义该区域会Available Regions of Interest列表中显示出时显示感兴趣区名称颜色包含素数
3 定义新感兴趣区ROI Tool话框中点击New Region定义效果见图17
图17 定义感兴趣区话框
图18 定义三感兴趣区图
4保存ROI
ROI TOOL窗口菜单：FileSave ROIs…
5合ROI
ROI TOOL窗口菜单：OptionsMerge Regions
6裁剪ROI
ROI TOOL窗口菜单：FileSubset Data via ROIs…
：菜单Basic Tools Subset Data via ROIs…
影进行注记
ENVI软件中图影中添加文符号色标条符号注记
1 幅图进行注记影菜单项选择项中选择Overlay →Annotation接着影窗口相应Annotation：Text话框会出现屏幕〔见图19〕
图19文〔Text〕模式Annotation话框
2 绘制图3D外表相似象进行注记绘制窗口菜单栏中选择Options →Annotation
注记类型：Annotation：Text话框中允许添加类型注记类型object拉菜单中进行选择注记类型文符号矩形椭圆边形折线箭头图例尺三北方图标〔Declination Diagrams〕图图例颜色标注影注记默认情况注记类型文〔Text〕话框中选项控制文注记字体〔Font〕〔Size〕颜色〔color〕位置角度
放置注记：方影窗口中添加文注记：
1 Annotation：Text话框中添加输入文〔for example Center South University of Forestry & Technology 〕
话框相应菜单参数设置中选择文字体颜色然影窗口适宜位置点击鼠标左键接着输入文会显示选点位置图110
图110 带注记影
注意：选择Window →Mouse Button Description 描述注记话框中鼠标键项功
2 鼠标左键拖动文注记圆柄窗口中放置文注记话框中改变相应区域设置值者鼠标左键拖动文符号注记改变注记属性位置
3 完成文注记设置点击鼠标右键锁定注记位置
保存恢复注记：
1 Annotation：Text话框菜单栏中选择File→Save Annotation保存影注记
2 翻开Output Annotation Filename 话框保存影注记该话框中指定保存路径保存文件名该文件扩展名ann
注意：果没影注记保存文件中关闭Annotation：Text话框时注记丧失
3 Annotation：Text话框中选择File→Restore annotation恢复保存注记文件
修改先前放置注记
先前已放置注记进行修改步骤进行修改：
1 Annotation：Text话框中选择Object→SelectionEdit
2 鼠标点击托拽出矩形框包含修改注记象
3 圆柄出现点击拖动注记圆柄修改应属性设置新注记象时样
添加网格
1影中添加公里网格信息影窗口中选择Overlay →Grid line Parameters话框
注意：影添加公里网格时影边框会动添加进
2 Grid line Parameters话框设置公里网线宽颜色公里网间隔修改公里网属性
图111 添加注记网格影
3 Edit Pixel Attributes话框中改变公里网标注网格线矩形边框交叉角颜色宽度公里网间隔完成属性设定点击Edit Pixel Attributes话框中ok钮做更改应影中
4 加公里网满意点击Grid line Parameters话框中Apply钮
保存输出影
添加标注进行图处理图进行保存影保存ENVI影文件格式者保存种通图格式然印者导入软件中直接通印机进行印输出
影保存ENVI影格式
1 影窗口菜单栏种选择File →Save Image File Output Display
To image file 话框会出现屏幕
2 选择24特彩色者8特灰阶进行输出然选择图形选项〔包括注记公里网〕边框设置果注记公里网添加彩色影显示中注记公里网会动列graphics选项中选择注记文件应输出图
3 需单项选择钮影结果输出Memory者file中果选择file输入输出文件名
4 点击ok保存影
实验二 影配准
专题旨介绍ENVI中影进行理校正添加理坐标ENVI进行影影校正
1图文件头文件修改
翻开显示SPOT数
首先翻开软件翻开软件带数库中SPOT数步骤：
1 ENVI菜单中选择file →open image file
2 enter data filename文件选择话框出现选择进入envidata目录bldr_spimg文件
3 点击ok
4 波段列表话框出现点击Grey scale单项选择钮鼠标左键点击相应波段名话框顶部列波段中选中SPOT波段选择波段名显示select Band：字段区域中
5 点击Load Band钮加载幅影新显示窗口中
修改ENVI头文件中图信息
1 波段列表中右键点击bldr_spimg文件名Map info图标弹出快捷菜单中选择Edit Map Information
Edit Map Information话框出现屏幕图21
话框列出ENVI中添加理坐标理信息调整ENVIMagic Pixel〔作图坐标系统起始点〕相应影坐标ENVI相应头文件信息图投影文件中识出图投影元图投影参数计算出影中意元理坐标输入图坐标输入理坐标〔纬度〕
2 点击ProjectionDatum文旁边箭头切换钮显示UTM Zone13 North图投影纬度度坐标ENVI处理程中进行转换
3 点击前DMS者DDEG钮分度分秒〔Degree－Minutes－Second〕十进制度〔Decimal Degrees〕间进行切换
4 点击Cancel推出Edit Map Information话框
图21Edit Map Information话框
5 修改图pixel size信息添加公里网格图标注
6 保存图影窗口中选择file →save image as →image file选择输出路径填写输出文件名称
2 影影配准
利SPOT图校正Landsat TM
翻开显示Landsat TM图
1 ENVI菜单中选择file →open image file
2 enter data filenames话框出现选择进入envidata目录bldr_reg子目录列表中选择bldr_tmimg文件
3 文件选择话框中点击open TM影波段加载波段列表中
4 列表中选择band3点击display#1钮拉式菜单中选择new display
5 点击Load Band钮TMband3波段影加载新显示窗口中
显示光标位置值
翻开显示影窗口滚动窗口者缩放窗口中光标位置信息话框步骤进行操作
1 影窗口菜单栏中选择tools →Cursor Locationvalue
2 影窗口滚动窗口缩放窗口TM影移动鼠标光标
注意坐标素单位出影基素坐标面带理坐标SPOT影
3 选择file →Cancel关闭Cursor Locationvalue话框
开始进行影配准
1 ENVI菜单栏中选择Map →Registration →Select GCPs：Image to Image
2 image to Image Registration话框中点击选择display＃1〔spot影〕作Base Image点击display＃2〔TM影〕作Warp image
图22image to Image Registration话框
3 点击ok启动动配准程序通光标放置两幅影相物点添加单独面控制点
注意：名物点需两张图般选择河流道路交叉点者拐角点交叉点拐角点选择图空间分辨率确定Quickbird空间分辨率061米Quickbird进行图校正时道路交叉点宽度应1米
4 启动Ground control Points Selection话框首先两张图找名物点现SPOT图找点鼠标移动该点样方法TM影找该点
5 两缩放窗口〔Zoom〕中查光标位置十字丝移应位置记录点XY值
图23缩放窗口〔Zoom〕话框
6 Ground Control Points Selection话框base XY文框中分输入286255SPOT影中光标位置移动相应位置点样方法Warp XY文框中分输入133264TM影中光标移动相应位置
7 Ground Control Points Selection话框中点击Add point该面控制点添加列表中点击Show list查面控制点列表见图24然样方式寻找余名物点
图24进行影配准Ground Control Points Selection话框Show list列表
注意：缩放窗口中支持亚元〔subpixel〕级定位缩放例越面控制点精度越面控制点选择通移动光标外直接名物点通鼠标进行点击完成show list话框中旦已选择4面控制点RMS误差会显示出
8 果选择控制点中某点误差删该点重新寻找新点代show list话框中鼠标点击该点选择delete删该点果选择控制满意话通Ground Control Points Selection话框选择Options →Clear All Points掉已选择面控制点
9 假设选择面控制点拟满意图进行校正前控制点信息进行保存具体方法Ground Control Points Selection话框中选择file →Save GCPs to ASCII…选择保存路径文件名控制点信息进行保存备次
图25影配准中Image to Image GCP List话框
校正影
里校正影TM影波段时校正波段影中波段
1 Ground Control Points Selection话框中选择Options →Warp Displayed Band
2 Registration Parameters话框中Warp Method钮菜单中选择RSTResampling钮菜单中选择Nearest Neighbor重采样法
图26 Registration Parameters话框
3 输入文件名bldr_tm1wrp点击ok
4 重复步骤1步骤2RST校正法相应选择BilinearCubic Convolution重采样法
5 结果分输出bldr_tm2wrpbldr_tm3wrp文件中
6 次重复步骤1步骤2次选择次项式Polynomial校正法Cubic convolution重采样法然选择Delaunay 三角网Triangulation校正法相应Cubic convolution重采样法
7 述结果分输出bldr_tm4wrpbldr_tm5wrp文件中
样方法TM影余波段进行校正述步骤区方余波段进行校正需重新找点需保存控制点信息直接调
Ground Control Points Selection话框中选择file →Restore GCPs from ASCII调保存控制点信息图进行校正校正图进行合成
8 合成步骤ENVI菜单栏中选择basic tools →lay stacking通lay stacking选择校正波段进行图层叠加叠加图进行保存保存文件名bldrtm_mimg
拟校正结果
动态链接拟校正结果：
1 波段列表中点击原始SPOT波段影名bldr_spimg然菜单中选择Grey scale显示方式显示图
2波段列表中选择bldr_tm1wrp文件Display＃拉式钮中选择New Display点击Load Band该文件加载新显示框中
3 影窗口中点击鼠标右键选择tools →Link →link displays
4 link displays话框中点击okSPOT影已添加理坐标bldr_tm1wrp影链接起
5 动态链接功带理坐标SPOT影进行拟
6 取消动态链接功选择tools →Link →unlink displays
查图坐标
1 影窗口菜单栏中选择Tools →cursor LocationValue
2 浏览带理坐标数集注意采样方法校正方法数值产生效果
3 选择File →Cancel关闭该话框
3 影图配准
局部逐步演示影图配准处理程许步骤影影配准步骤相似步骤会详细讨带理坐标 SPOT 影中获取图坐标矢量数字线划图〔DLG〕作基准数然基素坐标 Landsat TM 影进行校正匹配相应图数
31翻开显示 Landsat TM 影文件
ENVI 菜单中选择 File → Open Image File Enter Data Filenames话框出现选择进入相应文件夹目录表中lxmbldr_tmimg 文件点击 OKTM 影波段加载波段列表中时幅彩色影加载新显示窗口中
波段列表中点击 Gray Scale 钮选择波段 3点击 Load Band 钮 TM 影第 3 波段加载已翻开显示窗口中
32选择影图配准恢复控制点坐标
ENVI 菜单中选择 Map>Registration>Select GCPsImage to Map果翻开影显示窗口 Image to Map Registration 话框中点击选择包含该灰阶影显示窗口显示号投影列表中选择 UTM Zone 文框中输入 13设置素 30m点击 OK启动配准程序
校正影中光标移动图坐标面点〔幅图者 ENVI 矢量文件中[见局部]读取需图坐标〕 添加单面控制点
Ground Control Points Selection 话框中 E〔东〕 N〔北〕文框中手动输入图坐标然点击 Add Point 添加新面控制点 Ground Control Points Selection 话框中选择File>Restore GCPs from ASCII翻开
lxm PZhunpts 文件
Ground Control Points Selection话框中 点击Show List钮 Image to Map GCP List话框中查影图坐标实际影点预测点坐标 RMS误差
33矢量显示数字线划图〔DLGs〕添加图控制点：
ENVI 菜单中选择 File>Open Vector File>USGS DLG文件选择话框中选择 lxmbldr_rddlg 文件 Import Optional DLG File Parameters 话框中选择 Memory 单项选择钮点击 OK读入需数字线划图〔DLG〕数
矢量列表中高亮选择 ROADS AND TRAILS BOULDER CO 文件点击 Load Selected钮 Load Vector 话框中点击 New Vector Window该矢量加载新矢量显示窗口中
Vector Window #1 窗口中点击拖曳鼠标左键激活十字形光标光标处图坐标会 Vector Window #1 窗口底部列出
影显示窗口中选择 Tools>Pixel Locator输入 402 418然点击Apply影光标移动道路交叉口相应点注意缩放窗口中样获取亚元〔subpixel〕级定位精度
带十字形光标矢量窗口显示出图坐标（图）
矢量窗口中鼠标左键点击拖曳矢量光标十字形光标位需道路交叉口时松开鼠标左键矢量光标放置道路交叉口坐标 4775937444332400〔北纬 40 度 3 分3 秒东105度 15 分45 秒〕 矢量窗口中点击鼠标右键弹出快捷菜单中选择Export Map Location新图坐标会出现 Ground Control Points Selection 话框中 Ground Control Points Selection 话框中点击 Add Point添加该图坐标影素观察 RMS误差变化
34 RST三次卷积校正：
Ground Control Points Selection 话框中选择 Options>Warp File Input Warp Image 话框中高亮选择文件名 zhumuqian_tmimg点击 OK TM 6 波段进行校正出现 Registration Parameters dialog 话框中 Warp Method 选 RST Resampling 设置 Cubic Convolution background值改 255 output file 文框中输入输出文件名 zhumuqiantm_mimg点击 OK开始进行影图校正
35 显示结果进行评价
光标位置值〔Cursor LocationValue〕校正彩色影进行评价
波段列表中 点击Gray Scale单项选择钮 接着点击校正影波段Wrap（Band 3 zhumuqiantm_mimg） Display #拉式菜单钮中选择 New Display点击 Load Band加载幅校正 TM影注意校正影倾斜消 Landsat TM 轨道方影响原时影已带理坐标注意空间分辨率 30 米SPOT 影 10 米分辨率
实验三 图融合
数融合幅影组合单合成影处理程般高空间分辨率全色波段低空间分辨率光谱影进行融合达增强光谱影空间分辨率问题
进行数融合影文件必须含理坐标种情况空间采样会动运行果进行数融合影文件没理坐标影必须覆盖理区域相素影相方位
1 理坐标系统影融合
读取显示ER Mapper影
London TM影SPOT影二进制文件含ER Mapper头文件够动ENVIER Mapper读取程序读取
1 选择file →open external file →IP software →ER Mapper然进入lontmsp子目录选择lon_tmers
2 波段列表中单击RGB color单项选择钮次选取红绿蓝字段应波段点击Load RGB显示幅彩色TM影
3 选择file →open external file →IP software →ER Mapper然进入lontmsp子目录选择lon_spoters文件
4 波段列表中单击grey scale单项选择钮选取pseudo layer波段display拉菜单中选择new display点击load band显示幅灰阶SPOT影
调整影
1 波段列表中单击SPOT影发现空间尺寸2820×1569样方法知道TM影空间尺寸1007×560TM影素28米SPOT影素10米TM影必须28倍率调整产生SPOT影匹配10米元
2 选择Basic tools →Resize Data〔spatialSpectral〕选择lon_tmers单击okresize data parameters话框xfac文框中输入28yfac文框中输入28009选择保存路径文件名lontm点击ok调整TM影者直接samples文框中输入2820lines文框中输入1569选择保存路径文件名lon_tm点击ok调整TM影
图 31 Resize Data 话框
3 显示调整影选择Tools →link →Link Displays调整TM影SPOT全色影链接起动态链接分析拟两幅影
图32TM彩色合成影SPOT全色影
进行HIS数融合
解整处理程选择进行手工数融合第步TM彩色影转换色度饱度数值（称亮度值）〔huesaturationvalue（intensity）〕彩色空间高分辨率SPOT影换数值〔value〕波段拉伸01间满足正确数范围TM影获取色度饱度SPOT影中获取数值进行反转换转回红绿蓝彩色空间程产生出幅输出影中包含TM影中获取颜色信息SPOT影中获取空间分辨率信息
HSV正变换
1 ENVI菜单项选择择Transform →color Transforms →RGB to HSV然选择调整TM数作输入RGB影输入输出文件名点击ok执行变换
2 作灰阶影RGB彩色影显示色度饱度数值影
拉伸SPOT影换TM数值波段
1 ENVI菜单项选择择Basic Tools →Stretch Data单击lon_spot文件然点击ok
2 Data stretching话框Output Data中min文框中输入0Max文框中输入10输入输出文件名单击okSPOT影数拉伸浮点型范围010
图33 Data stretching话框
HSV反变换
1 ENVI菜单中选择Transform →Color Transforms →HSV to RGB选择变换TM影HueSaturation波段作变换HS波段
2 选择拉伸SPOT影作V波段点击okHSV to RGB Parameters话框中输入输出文件名点击ok进行反变换
图34 HSV to RGB Parameters话框
显示结果
1 波段列表中点击RGB Color单项选择钮序点击反变换RGB波段点击Load RGB钮显示幅融合TMSPOT彩色影
2 显示融合影选择Tools →Link →Link Displays融合影调整TM影SPOT全色影链接起动态分析拟影
进行动HSV变换融合
1 ENVI菜单项选择择Transform →Image Sharpening →HSV
2 果调整TM彩色影已显示窗口中Select Input RGB话框中直接选择适宜影显示窗口否Select Input RGB话框中选择Red Layer〞Green Layer〞Blue layer〞应调整TM影波段然点击ok
3 High Resolution Input File话框中选择SPOT影点击ok
4 输入输出文件名：lontmspimgHSV Sharpening Parameters话框中点击ok
结果显示链接拟
1 显示融合彩色影波段列表中选择RGB Color单项选择钮RGB波段中选择融合影中相应波段点击Load RGB
2 通影窗口菜单中选择Tools →Link →Link Displays原始TM彩色影SPOT影动融合彩色影进行链接拟
图35 TM彩色影SPOT影融合彩色影链接拟
2 带理坐标系统影融合
翻开进行融合图
1选择file →open image file 然进入envidata目录双击p123r043ETM文件夹选择p123r043_7p19991224_z49_multihdr点击翻开影出现波段例表中
2 波段列表中点击RGB Color选择R〔4〕G〔5〕B〔3〕三波段进行波段显示
3 选择file →open image file 然进入envidata目录双击p123r043ETM文件夹选择p123r043_7p19991224_z49_panhdr点击翻开影出现波段例表中
4波段列表中点击Grey scale选择p123r043_7p19991224_z49_panhdr点击Display＃1拉菜单点击New display点击Load band
Color Normalized〔Brovey〕sharpening〔标准化彩色变换融合〕
1ENVI菜单项选择择Transform →Image Sharpening →Color Normalized〔Brovey〕
2 select Input for color Bands话框中选择Available Bands List点击ok
图36 select Input for color Bands话框
3 出现波段列表话框中p123r043_7p19991224_z49_multihdr文件中选择三波段应RGB选项点击ok
4 high resolution Input file话框中选择p123r043_7p19991224_z49_panhdr文件波段点击ok
图37 high resolution Input file话框
5 Color Normalized sharpening话框中Resampling重采样方法选择Nearest Neighbor输入输出文件名保存路径点击ok会进行动融合融合结果会出现波段例表中
图38 Color Normalized sharpening话框
6 查波段融合结果波段列表话框中选择RGB color显示方式选择R〔4〕G〔5〕B〔3〕点击Load RGB显示融合图
Gramschmidt spectral sharpening〔GS变换融合〕
1ENVI菜单项选择择Transform →Image Sharpening →Gramschmidt spectral sharpening
2 select low spatial resolution multi band input file话框中选择p123r043_7p19991224_z49_multi文件点击ok
4 high resolution Input file话框中选择p123r043_7p19991224_z49_pan文件波段点击ok
5 Gramschmidt spectral sharpening话框中Resampling重采样方法选择Nearest Neighbor输入输出文件名保存路径点击ok会进行动融合融合结果会出现波段例表中
6 查波段融合结果波段列表话框中选择RGB color显示方式选择R〔4〕G〔5〕B〔3〕点击Load RGB显示融合图
图39 select low spatial resolution multi band input file话框
PC spectral sharpening（成分变换融合）
1ENVI菜单项选择择Transform →Image Sharpening →PC spectral sharpening
2 select low spatial resolution multi band input file话框中选择p123r043_7p19991224_z49_multi文件点击ok
4 high resolution Input file话框中选择p123r043_7p19991224_z49_pan文件波段点击ok
5 PC spectral sharpening话框中Resampling重采样方法选择Nearest Neighbor输入输出文件名保存路径点击ok会进行动融合融合结果会出现波段例表中
6 查波段融合结果波段列表话框中选择RGB color显示方式选择R〔4〕G〔5〕B〔3〕点击Load RGB显示融合图
注意：改变图融合重采样方法否会图显示产生影响图清晰度色彩两方面拟图融合方法图影响
实验四 波段组合计算图增强
专题中ETM数象介绍图处理程中波段组合方式波段间运算方式数拉伸增强处理程
翻开影
1选择file →open image file 然进入envidata目录双击p123r043ETM文件夹选择p123r043_7p19991224_z49_multihdr点击翻开影出现波段例表中
2 波段列表中点击RGB Color选择R〔4〕G〔5〕B〔3〕三波段进行波段显示
3 点击Load RGB幅假彩色图显示影窗口中
波段组合
1 显示影中p123r043_7p19991224_z49_multihdr文件三波段文件6波段里尝试6波段中选取3波段组合拟波段组合间图显示效果
2 波段列表中点击RGB Color选择R〔7〕G〔4〕B〔1〕三波段进行波段显示
图41波段组合间影显示效果
3 波段列表中点击Display＃1拉菜单中选择New display点击Load RGB幅真彩色图显示影窗口中
4 动态链接拟前两幅影效果影窗口中选择tools →link →link displayslink displays话框中点击ok两幅影链接起
5 波段列表中选择余波段组合方式加拟
波段运算
Band math功户提供灵活图处理工具Band math话框中实现波段间加减等运算里p123r043_7p19991224_z49_multihdr文件象进行波段运算
1 菜单栏中选择file →open image file 然进入envidata目录双击p123r043ETM文件夹选择p123r043_7p19991224_z49_multihdr点击翻开影出现波段例表中
2 菜单栏中选择Basic Tools →Band math出现band math话框
图42 Band math话框
3 enter an expression文框中输入需进行波段计算IDL〔Interactive Data Language〕表达式变量代波段名文件名变量名必须字符b〞B〞开头面着5数字字符
注意：例假设想计算三波段均值数学表达式：
〔float（b1）+float（b2）+float（b3）〕3输入文框中该表达式中三变量b1表示第变量b2表示第二变量 b3表示第三变量变量前浮点型字节防止计算时出现字节溢出错误输入效表达式点击Add to list表达式添加话框中点击ok
4 变量赋值variables to Bands Pairings话框中允许输入波段列表中波段赋enter an expression〞文框中包含变量
图43 variables to Bands Pairings话框
值原先实例中变量b1〞variables to Bands Pairings话框中标签Variables used in expression〞文框点击变量B1〞
标签Available Bands List〞列表中点击需波段数集显示变量名注意旦第波段选中具相空间尺寸波段显示波段列表中样方法B2〞B3〞等赋予值
5 输出波段运算结果三变量选定相应波段点击ok运算结果进行输出Variables to Bands Pairings话框中选择file〞Memory〞切换钮运算结果显示波段列表中
提醒：波段均值计算外尝试波段间值加减者法等系列运算波段运算般满足四条件：
1〕必须IDL语言书写波段运算表达式2〕输入波段必须具相维数3〕表达式中变量必须Bn〔bn〕命名4〕结果波段必须输入波段维数相
交互式数拉伸〔Interactive Stretching〕
翻开进行拉伸数p123r043_7p19991224_z49_multihdr影窗口选择Enhance→interactive stretching交互式度拉伸话框中显示输入直方图输出直方图说明前输入数分应拉伸两条垂线〔虚线〕说明前拉伸值值彩色影说直方图颜色选波段颜色致〔系统默认显示红色波段直方图〕交互式直方图窗口底部会列出拉伸类型直方图源浏览幅彩色图绿波段蓝波段直方图点击窗口中G〞B〞切换钮拉伸直方图变换动应图选择Options →Auto Apply般情况点击Apply〞钮变换应图通取消选中Options →Auto Apply关闭动功
图44 交互式度拉伸话框
设置拉伸参数
交互式度拉伸话框中Options〞菜单鼠标交互功拉伸设置参数时浏览交互式直方图窗口中信息
1） 更改拉伸值值
更改拉伸值值意条垂直虚线点击鼠标左键然拖放新位置者直接话框顶部文框中输入需值
提示：文框中指定值时输入需实际值〔例值37值68〕者数百分〔例值2％值98％〕设置数范围外值根输入直方图数范围进行拉伸输入数输出直方图会动更新反映输入直方图中变化应新拉伸显示数分布
2） 锁定拉伸条〔locking stretch bars〕
锁定值值拉伸条〔垂直虚线〕间距离选择Options →locking stretch bars拉伸条间距离锁定时移动两拉伸条解锁定取消选中locking stretch bars
拉伸类型菜单
交互式度拉伸话框中Stretch_type菜单中包含供选择交互式拉伸类型列表
1） 线性度拉伸
线性度拉伸系统默认交互式拉伸线性拉伸值值分设置0255两者间值设置中间线性输出值
交互式度拉伸话框中选择Stretch_type →linear contrast stretch限定输入值鼠标左键移动输入直方图中垂直线〔白色虚线〕需位置者stretch文框中输入需DN值数百分〔例值5％值95％〕重新显示原拉伸选择Options →Reset Stretch
2） 分段线性度拉伸〔piecewise linear〕
分段线性度拉伸通鼠标输入直方图中放置点进行交互限定点间局部采线性拉伸
选择Stretch_type →piecewise linear转换函数〔初始条白色直线〕绘制输入直方图输入直方图位置点击鼠标中键转换函数增加节点绘制线段端点绘制节点标记连接起
图45分段线性度拉伸话框
移动点位置标记点击鼠标左键然拖放新位置删点标记点击鼠标右键手动键入输入输出值选择Options →Edit piecewise linear结果绘制输出直方图中该直方图显示应新拉伸数分布情况点击Apply拉伸应显示数
3〕高斯度拉伸〔Gaussian〕
系统默认Gaussian拉伸均值DN127应0255正负3标准差值进行拉伸选择Stretch_type →Gaussian输入拉伸值值输出直方图条红色曲线显示选择Gaussian变换函数拉伸数分布呈白色叠加显示红色Gaussian函数
注意：手动输入需标准差选择Options →Set Gaussian Stdv点击Apply拉伸应显示数
4〕直方图均衡化度拉伸
动缩放数直方图bin中DN数相均衡选择Stretch_type →Equalization输入直方图显示未修改数分布输出直方图条红色曲线显示均衡化函数拉伸数分布呈白色叠加显示点击Apply拉伸应显示数
5〕方根度拉伸〔Square Root〕
计算输入直方图方根然应线性拉伸选择Stretch_type →Square Root输入直方图显示未修改数分布输出直方图条红色曲线显示方根函数拉伸数分布呈白色叠加显示点击Apply拉伸应显示数
6〕定义度拉伸直方图匹配〔Arbitrary〕
Arbitrary选型允许输出直方图顶部绘制形状直方图图直方图相匹配
选择Stretch_type →Arbitrary输入直方图未修改数分布通点击住拖放鼠标左键Output Histogram窗口绘制输出直方图定义直方图绿色显示点击Apply拉伸应显示数
注意：删直方图点击鼠标中键点击鼠标右键接受绘制输出直方图数统计资料匹配图中输出直方图红色显示绘制直方图匹配数函数白色曲线绘制
7〕直方图匹配
定义度拉伸功幅图直方图幅图直方图进行匹配幅图中获输入输出直方图图表顶部Input Histogram〞Output Histogram〞文标签点击鼠标左键拖放名字拖放定义输出直方图中然释放钮导入直方图红色绘制输出直方图拉伸导入直方图匹配点击Apply拉伸应显示数
8〕定义查找表拉伸〔Defining Look UpTable〕
户定义查找表输入DN值拉伸输出值恢复交互式定义
选择Stretch_type →User Defined LUT选择Options →Edit User Defined LUT出现编辑话框时包含输入DN值应拉伸输出值列表显示Edit User Defined LUT〞标签值反映前拉伸情况值点击进行编辑出现Edit Select item〞文框中时输入需值然回车键
注意：该功启动重新设置查找表初始值点击Reset〞点击ok点击Apply拉伸应显示数
图46 Edit User Defined LUT话框
图增强滤波
图增强功影窗口中enhance拉菜单中线性增强高斯增强直方图均衡增强方根增强种类型种增强类型三选择影窗口卷积窗口缩放窗口分进行增强处理
滤波通常通特定空间频率图增强空间频率通常描述亮度DN值距离方差图包括种空间频率消幅图高频信息图显滑ENVI软件中进行图滤波菜单栏中filter进行操作包括卷积滤波〔Convolutions and Morphology〕形态学滤波纹理滤波适应滤波频率滤波
卷积滤波指空间域图进行滤波形态学滤波形态根底图进行处理纹理滤波图中提取纹理信息适应滤波器抑制噪声时保存图锋利信息细节傅立叶滤波频率域图进行滤波
1） 卷积滤波
菜单中选择filter →Convolutions and Morphology出现Convolutions and Morphology Tool话框
图47 Convolutions and Morphology Tool话框
选择滤波器类型Convolutions and Morphology Tool话框中点击Convolutions拉菜单中选择滤波器类型〔高通低通拉普拉斯直通高斯高通高斯低通中值sobelRoberts户定义滤波〕系统默认滤波器类型高通类型
编辑变换核Convolutions and Morphology Tool话框中鼠标点击Kernel Size〞增减箭头时变换核尺寸2单位增减鼠标中键点击Kernel Size〞 增减箭头时变换核尺寸10单位增减鼠标右键点击Kernel Size〞 增减箭头时变换核尺寸恢复默认尺寸变换核默认设置正方形变换核变非正方形核选择Option →Square Kernel：No默认Image Add Back0
2） 纹理滤波器〔基概率统计滤波〕
Texture选项应基概率统计二阶概率统计纹理滤波许图包含区域亮度变化特征仅局限亮度值纹理指图色调作等级函数空间变化定义纹理清晰区域灰度等级相纹理区定拟接
菜单中选择filter →Texture →occurrence MeasuresTexture input file话框中选择进行滤波文件点击ok
图48 Texture input file话框
Occurrence Texture Parameter话框中选择输出文件路径核文件名余选项采默认值点击ok文件进行纹理滤波处理选纹理图计算出放置波段例表中
图49 Occurrence Texture Parameter话框
3） 适应增强型Lee滤波器
适应滤波运围绕元标准差计算新元值般说原始元值基周围效元〔符合标准差元〕计算新值代增强型Lee滤波器Lee滤波器改良根单独滤波窗口中计算出统计〔方差系数〕数进行滤波元分3类型中：相似元〔homogenous〕差异元〔heterogeneous〕指目标元〔point target〕
菜单中选择filter →Adaptive →enhance LeeEnhance Lee filter Input file话框中选择进行滤波文件点击okEnhance Lee filter Parameters话框中输入filter SizeDamping FactorHomogenous Areas等参数采默认方式参数填写完毕选择输出路径输出文件名点击ok
图410 Enhance Lee filter Parameters话框
4〕快速傅立叶变换滤波器
FFT Filtering〔Fast Fourier Transform Filtering 快速傅立叶变换滤波 〕图变换成显示空间频率成分合成输出图包括操作：频率域滤波器交互式建立滤波器应FFT原始数空间逆变换
正FFT滤波
翻开影ENVI菜单栏中选择Filters →FFT Filtering →Forward FFT出现Forward FFT Input File话框选择进行滤波文件点击ok
Forward FFT Parameters话框中选择输出file〞memory〞点击ok开始FFT计算
定义FFT滤波器
Filter Definition选项交互式直接定义滤波器通显示正变换图中绘制定义滤波器滤波类型包括：circle pass and cutuserdrawn pass and cut
ENVI菜单栏中选择Filters →FFT Filtering →Filter Definition出现Filter Definition选择话框选项中选择：
果前显示中包含幅FFT影选择包含FFT图显示号选择No Display〞指出该滤波器某特定图相关果前显示中含FFT图Samples〞Lines〞文框中键入数值限定滤波器
图411 Filter Definition话框
Filter Definition选择话框中选择Filter_type滤波器类型Circular pass〞Circular cut〞滤波器需Radius〞文框中元单位输入滤波半径
Band Pass〞Band Cut〞滤波器Inner Radius〞Outer Radius〞文框中元单位键入需值
User Defined〞User Defined Cut〞滤波器允许ENVI注记导入滤波器〔般FFT影annotation功〕
箭头钮选择Number of Border Pixels〞参数细化滤波器值0表示没滑
选择输出file〞memory〞点击apply构建FFT滤波器
反FFT变换
反FFT变换程序实际包括两步骤先应FFT滤波然FFT图反变换回原始数空间
选择Filter →FFT Filtering →Inverse FFT出现Inverse FFT Input File话框选择输入文件波段假设需选取意子集点击ok图进行反变换
实验五 图分类
专题介绍遥感影分类操作处理拟监督分类非监督分类影分类影进行相应处理处理包括聚合〔clump〕处理筛选〔sieve〕处理类〔combine〕处理精度评估
翻开显示分类影
1ENVI菜单中选择File →Open Image file选择进入envidata文件夹目录点击can_tm子目录列表中选择can_tmimg文件点击ok
2 波段列表中选择RGB Color单项选择钮然鼠标左键次点击波段4波段5波段3选波段会话框中适文框中显示出
3 点击Load RGB钮该影加载新显示窗口中影窗口影进行简单区分目视判读工作
非监督分类
ENVI菜单中选择Classification →Unsupervised →选择Kmeans者IsoData
K均值〔Kmeans〕分类法
K－均值非监督分类器聚类分析方法需首先选定需分类数机查找聚类中心位置然迭代重新配置直达优化波谱分类
1ENVI菜单中选择Classification →Unsupervised →KmeansClassification Input file话框中选择分类影分类波段数点击ok
2 Kmeans parameters话框中确定分类数6类余参数采默认方式选择文件输入名称输出路径点击ok影进行分类
3 波段列表中添加新图显示窗口Display＃1拉菜单中选择New Display
4 双击Kmeans〔can_tmimg〕分类图分类图显示新影窗口中
5 影窗口菜单中选择Tools →Link →Link Displays然话框中点击ok
图51 Kmeans parameters话框
6 鼠标左键影点击拖动动态叠加显示区域Kmeans分类结果原始彩色合成影进行拟
7 处理完成选择Tools →Link →Unlink Display关闭动态链接
IsoData〔迭代组织数分类技术〕
IsoData非监督分类法计算数空间中均匀分布类均值然距离规剩余元进行迭代聚合次迭代重新计算均值根新均值元进行分类处理程持续类元数变化少选元变化阈值者达迭代次数
1ENVI菜单中选择Classification →Unsupervised →IsoDataClassification Input file话框中选择分类影分类波段数点击ok
2 IsoData parameters话框中确定分类数6类余参数采默认方式选择文件输入名称输出路径点击ok影进行分类
3 波段列表中添加新图显示窗口Display＃1拉菜单中选择New Display
4 双击IsoData〔can_tmimg〕分类图分类图显示新影窗口中
5 影窗口菜单中选择Tools →Link →Link Displays然话框中点击ok
6 鼠标左键影点击拖动动态叠加显示区域IsoData分类结果原始彩色合成影进行拟住鼠标左键时次点击鼠标右键触发第三幅影动态叠加显示IsoData分类结果Kmeans分类结果进行拟
7 处理完成选择Tools →Link →Unlink Display关闭动态链接
监督分类法
监督分类种边训练边学分类方法进行监督分类时首先训练分类样确定分类数样进行保存然选择监督分类方法分类结束进行分类处理工作
1选择非监督分类影can_tmimg文件影窗口中选择Overlay →Region of Interest接着应显示窗口ROI Tool话框
图52 ROI Tool话框
2 ROI Tool话框中选择ROI_typepolygon然影窗口中绘制边形该边形代表新创立感兴趣区完成步请面步骤进行：
l 影窗口中点击鼠标左键建立感兴趣边形第点
l 次点击鼠标左键序选择更边线点点击鼠标右键闭合该边形鼠标中键删新定义点者删已闭合整边形次点击鼠标右键固定边形位置
l 通选ROI Controls话框顶部应单项选择钮感兴趣区缩放窗口中定义
感兴趣区域定义完会话框区域列表中显示出时显示感兴趣区域名字颜色包含素总数
3 定义新感兴趣区点击New Region钮点击Edit钮翻开Edit ROI Parameters话框中输入感兴趣名字选择感兴趣区颜色填充方式述步骤定义新感兴趣区分类影定义6类
图53 感兴趣话框训练样
4 训练样进行保存ROI Tool话框中选择file →save ROIs
5 Save ROIs to file话框中选择Select All Items点击choose选择保存路径文件名
6 点击ok保存文件
图54 Save ROIs to file话框
距离法分类〔Minimum Distance〕
距离分类法感兴趣区均值矢量计算未知元类均值矢量欧氏距离非户指定标准差距离阈值否元分类感兴趣区中接类
1ENVI菜单中选择classification →supervised →Minimum DistanceClassification Input File话框中点击选择进行监督分类文件can_tmimg点击ok
图55 Minimum Distance Parameters话框
2 Minimum Distance Parameters话框中选择Select All Items set Max stdev from Meanset Max Distance Error采默认方式选择输出路径文件名话框右边选择输出规图Output Rule Image钮鼠标点击
3点击ok进行距离监督分类计算结果出现波段列表中
4波段列表中添加新图显示窗口Display＃1拉菜单中选择New Display
5 双击min Dist〔can_tmimg〕分类图分类图显示新影窗口中
注意：监督分类方法采距离分类外采似然法行六面体等分类技术分类步骤距离法致
影分类处理
需分类影进行处理评价分类效果精度者类概括出导入图影矢量GIS中
分类统计
分类统计分类图中提取统计信息统计信息根统计信息〔值值标准差方差等〕直方图者选类中计算出均波谱
1 选择classification →post classification →class Statistics进行统计处理选择分类影然点击ok
2 接着选择原始分类影can_tmimg点击ok
3Class Selection话框选择进行统计类点击Select All Items然点击ok
4 Compute Statistics Parameters话框中选择计算统计信息点击Compute Statistics Parameters话框底部ok钮然根选择统计选项绘制图报表会出现屏幕
混淆矩阵〔Confusion Matrix〕
ENVI计算混淆矩阵功允许两幅图〔分类影真实影〕者分类影感兴趣区进行拟面真实影幅分类影者根面真实测量生成影
1 选择Classification →Post Classification →Confusion Matrix →两选项供选择Using Ground Truth Image者Using Ground Truth ROIs
2 专题中选择Using Ground Truth ROIs进行混淆矩阵计算
3 Classification Input file话框中选择分类影mindistanceimg点击ok
4 Confusion Matrix Parameters话框中点击ok钮
5 查混淆矩阵混淆影通波谱剖面廓线CursorLocation Vlaue分类影原始反射率影中感兴趣区进行拟确定误差源
图56 分类影混淆矩阵结果
聚合筛选处理〔Clump and Sieve〕
1 分类影进行筛选处理选择ENVI菜单栏中classification →Post classification →sieve classes选择分类影输入memory中然点击ok
2 输出筛选处理结果作聚合处理输入图选择classification →Post classification →clump classes存中选择先前生成筛选处理影点击ok
3 输出Memory点击Clump Parameters话框中ok钮
4 拟处理前三幅影
合类〔Combine Classes〕
1选择ENVI菜单栏中classification →Post classification →combine classes选择分类影点击ok
2 combine classification parameters话框中选择进行合两类先选择合类选择新输出类〔列unclassified 入bare类组成新bare类Input class选择unclassifiedOutput class选择bare〕点击ok生成新合类
图57 Combine Classification Parameters话框
修改类颜色
显示分类影时通修改类颜色改变特定类应颜色
1 影窗口中选择tools →Color Mapping →Class Color Mapping
2 Class Color Mapping话框中点击某类名字拖动相应颜色条者输入需颜色值改变类颜色做改变会立刻应分类影进行永久性改变话框中选择Options→ Save Changes
实验六 图制图
专题旨介绍ENVI图制图处理流程注记功添加额外信息ENVI图制图功够方便快捷交互式添加图例等素幅影绘制成图制图程般包括步骤组成首先ENVI快速制图〔QuickMap〕动生成根模板〔者恢复保存模板〕然ENVI注记功影叠合功需进行交互式制图快速制图允许设定图例输出页方位够选择影空间子集进行制图方便添加根图素图公里网例尺图图标标识图投影信息根图注记
翻开进行制图Landsat TM影
1 ENVI菜单栏中选择file →Open Image file
2 出现Enter Input Data File文件选择话框中选择进入envidata目录中ys_tmsub子目录列表中选择ysratioimg文件点击Open相应文件波段会列波段列表中根默认设置573134率波段会动加载R〞G〞B〞字段时选择RGB Color钮
3 点击话框底部Load RGB钮该影加载新显示窗口中
生成快速制图模块
1 影显示窗口菜单中选择File →QuickMap →New QuickMap翻开QuickMap Default layout话框
话框修改输出页页方位图例
图61 QuickMap Default layout话框
2 QuickMap Default layout话框中Map Scale改200000余参数采默认方式点击ok选定影子集
3 QuickMap image Selection话框中鼠标左键点击红色方框左角拖动方框选中整影
注意：选择影子集输出成图鼠标左键放置调整红色矩形框
图62 QuickMap image Selection话框
4 点击okQuickMap Parameters话框会出现屏幕
5 话框中鼠标左键点击Main Title文框键入文Yellowstone National Park Image Map
6 话框中鼠标右键点击Lower left Text文框弹出菜单栏中选择Load Projection InfoENVI头文件中加载影投影信息
7 话框中鼠标右键Lower Right Text文框键入Map Generated Using然回车〔Enter键〕次输入ENVI QuickMapCopyright2001Research System Inc次键入时回车键
图63 QuickMap Parameters话框
8 根专题目里建议ScaleGrid Lines North 前面复选框〔check box〕选中更改设置然QuickMap Parameters话框右边复选框选项进行
9 点击Declination Diagram复选框选中
10 话框底部选择Save Template输入Output Filename文件名ysratioqm点击ok快速制图结果保存快速制图模板文件模板处理相素影时进行调需显示需影选择File →QucikMap →from Previous Template恢复已保存快速制图模板
11 QuickMap Parameters话框左角点击Apply标准ENVI显示窗口中显示快速制图结果果需修改QuickMap Parameters话框中设置然点击Apply更新显示结果
12 ENVI显示窗口中查结果观察图公里网例尺指北针默认文位置
定义图版面设计中元素
ENVI提供种定制图制图选项包括添加虚边框〔Virtual borders〕文注记公里网等值线绘图插入矢量叠合分类叠合ENVI影窗口滚动窗口者缩放窗口进行定制制图操作〔果滚动窗口拖住角调整显示〕
1〕添加虚边框
默认ENVI显示包含影周围没空白空间图制图般需某图素放置影外ENVI提供虚拟边框〞功注记放置影周围创立新影虚拟边框采种方法添加影次介绍
ENVI快速制图中公里网动添加
虚拟边框会添加图影中快速制图中公里网相协调默认公里网显示出参见具体说明解修改公里网属性参数必边框动添加影两边
1 改变默认边框快速制图影显示窗口菜单中选择Overlay →Grid linesGrid line Parameters话框出现选择Options →Set Display Borders翻开Display Borders话框
2 输入100400150100点击ok
图65 设置Display Borders话框
新虚拟边框属性马会应图影中果保存公里网参数文件边框属性信息会时保存公里网参数文件中恢复公里网属性时边框属性信息会恢复
显示参数设置
Display Preference话框改变虚拟边框属性显示设置
1 快速制图影显示窗口菜单中选择File →Preference接着Display Preference话框会出现话框顶部包含刚刚介绍相类似文框
2 输入需数值边框选择需颜色
3 点击OK新边框属性会应图影中
注记功
ENVI中注记功调整虚拟边框属性
1 快速制图影显示窗口菜单中选择Overlay →Annotation
2 Annotation话框出现选择Options →Set Display Borders翻开Display Borders话框
3 输入需边框属性参数点击OK新虚拟边框属性马会应图影中果保存注记文件边框属性信息会时保存注记文件中恢复公里网属性时边框属性信息会恢复
2〕添加公里网
ENVI支持时显示素公里网理坐标公里网理坐标〔纬度度〕网公里网应图影中时100素虚拟边框〔没添加虚拟边框情况〕会动添加图影中协调公里网标签〔labels〕添加修改图影公里网参考步骤进行
1 快速制图影显示窗口菜单中选择Overlay →Grid Lines出现Grid Lines Parameters话框中默认公里网间距会显示出
2 Grid Spacing文框中输入4000公里网间距改4000米改变公里网标签属性参数选择Options →Edit Map Grid Attributes者Edit Geographic Grid Attributes修改选公里网属性样Quick Map Parameters话框中点击Additional Properties钮修改公里网属性参数
3 点击OK确定选属性
4 Grid Line Parameters话框中点击Apply新公里网应图影显示中保存公里网属性便续Grid Parameters菜单栏中选择File →Save Setup选择输出文件样保存公里网属性参数模板便加载图制图中
3〕操作处理注记
注记ENVI显示者图制图中插入〔放置〕图素种常规手段ENVI根注记功种图素放置位置
1 快速制图影显示窗口菜单中选择Overlay →Annotation翻开Annotation话框
2 Annotation话框菜单栏Object拉菜单中选择需注记素
3 选择ImageScroll者Zoom单项选择钮指定注记放置窗口
4 鼠标左键放置注记素点击鼠标右键进行确定锁定注记位置
图67 注记〔Annotation〕话框
Object拉菜单中文〔Text〕符号〔symbol〕形状注记线段箭头注记图例尺〔Scale Bar〕三北方表〔Declination〕图例注记〔Map key〕嵌入式影注记〔Image〕绘制图〔Plot〕嵌入式注记等
专题中求添加注记
①文文档：Landsat TM Dtata
Ratios 573134（RGB）
②图例注记〔Map key〕
Annotation话框菜单栏Object →Map key选择Edit Key Items添加删者修改单图例项修改成功回车键鼠标左键放置图例鼠标右键锁定图例位置
③嵌入影注记
首先确保插入影已波段列表中列出果没波段列表中话先翻开插入影ys_loctif
Annotation话框菜单栏Object →ImageAnnotation：Image话框中选择Select New Image选择插入影ys_loctif
图68 Annotation：Image话框 图610 Annotation：Declination话框
波段列表中选取插入影果需处理中选取空间子集调整影鼠标左键放置嵌入影鼠标右键锁定嵌入影
图611 ENVI注记功定制图制图
注记添加完毕注记文件进行保存Annotation话中选择File →Save Annotation注记文件进行保存
输出图制图结果
1 快速制图影显示窗口菜单中选择File →Save Image As →Image file
2 Out Display to Image file话框中Out file type类型选择BMP者JPEG选择 保存路径文件名点击ok进行保存幅完整图制图完成
文档香网(httpswwwxiangdangnet)户传
《香当网》用户分享的内容，不代表《香当网》观点或立场，请自行判断内容的真实性和可靠性！
该内容是文档的文本内容，更好的格式请下载文档
遥感图处理实验指导
实验 图处理软件ENVI4功介绍
实验二 影配准
实验三 图融合
实验四 波段组合计算图增强
实验五 图分类
实验六 图制图
实验 图处理软件ENVI4功介绍
1 ENVI4支持文件格式
ENVI通栅格数格式包含简单二进制文件相关ASCII头文件读该文件格式允许ENVI影文件包括身嵌入头文件信息影文件
通栅格数会存储二进制字节流通常BSQ〔波段序〕BIP〔波段元交叉〕者BIL〔波段行交叉〕方式进行存储
ENVI软件支持数类型包括字节型整型符号整型长整型符号长整型浮点型双精度浮点型复数型双精度复数型64位整型符号64位整型
2 ENVI窗口显示
ENVI软件时翻开图通常三窗口影窗口〔Main Image Window〕滚动窗口〔Scroll Window〕缩放窗口〔Zoom Window〕见图11
图11 ENVI显示窗口组：影窗口滚动窗口缩放窗口
21 ENVI菜单
ENVI菜单位显示窗口顶部默认情况显示影窗口顶部提供交互式影显示分析功见图12果选显示组中包括影窗口菜单会出现滚动窗口者缩放窗口顶部ENVI菜单样该菜单中选择意菜单项选择项外显示窗口点击右键会弹出快捷菜单通菜单访问显示函数改变显示设置
图12影窗口中overlay菜单快捷菜单
22 波段列表
ENVI翻开影文件者文件中单波段波段列表〔Available Bands List〕特殊ENVI话框包括翻开文件中影波段相关图信息列表见图13波段列表彩色灰阶影加载显示窗口中直接翻开新显示窗口点击话框底部Display钮中选择显示窗口号翻开相应显示窗口然点击Gray Scale 者RGB单项选择钮列表中点击波段名选择需波段加载显示影
提示：果加载单波段影直接波段双击完成
图13波段列表话框
3 ENVI根功
翻开影文件
首先启动ENVI软件双击软件图标翻开影文件分步骤：
1 选择File→Open Image file屏幕弹出文件选择话框Enter Input Data File〞
2 选择进入实验课程提供数目录envidata中双击翻开文件夹选择can_tm文件夹文件夹中选择can_tmrimg文件然点击ok
3 波段列表中会显示刚刚选择文件can_tmrimg鼠标左键点击选择影显示方式〔Gray ScaleRGB〕选择某影波段三波段进行显示选波段标Selected Band〞显示出
4 点击Load BandLoad RGB影加载新显示窗口
图加载熟悉窗口界面图窗口滚动窗口缩放窗口进行相关操作
显示光标位置
显示鼠标光标位置值影窗口菜单栏者ENVI菜单栏中选择Window→Cursor LocationValue者影窗口中点击鼠标右键弹出快捷菜单中选择Cursor LocationValue
接着屏幕出现Cursor LocationValue话框显示出光标影窗口滚动窗口者缩放窗口中位置见图14该话框显示十字丝光标应象素屏幕颜色值实际数值
关闭话框Cursor LocationValue话框顶部菜单中选择file→cancel
旦Cursor LocationValue话框翻开隐藏显示该话框影窗口中双击鼠标左键
图14 Cursor LocationValue话框
显示影剖面廓线
交互式选择显示X轴〔水〕Y轴〔竖直〕Z轴〔波谱〕剖面廓线图剖面廓线图显示穿影横线X线Y者波谱波段Z数值
1 影窗口菜单栏中选择Tools→Profiles→X profle会翻开绘制窗口该窗口根影中选择行绘制出幅数值列号〔sample number〕间关系曲线图见图15
2 重复述程选择Y profile绘制出幅数值行号〔line number〕间关系曲线图选择Z profile绘制出相应波谱面廓线图
提示：通意影窗口中快捷菜单翻开Z轴波谱剖面廓线图
3 选择Window→Mouse Button Descriptions查鼠标键剖面廓线显示窗口中键功描述
4 放置三剖面廓线窗口便时
5 影移动移动十字丝查三幅影剖面廓线图新位置样显示新数值
6 绘制图窗口中选择File→Cancel关闭绘制窗口
图15水轴〔X〕Y轴〔竖直〕Z轴〔波谱〕剖面廓线绘制曲线图
进行快速拉伸
影窗口缩放窗口者滚动窗口中默认参数数进行快速度拉伸选择影窗口菜单栏中Enhance菜单进行种样度拉伸〔线性拉伸0－255间线性拉伸2％线性拉伸高斯拉伸均衡化拉伸方根拉伸〕
1 影窗口缩放窗口者滚动窗口中影作拉伸数源尝试进行种类型拉伸变换
2 拟窗口显示组中线性高斯均衡方根拉伸结果
显示交互式散点图
ENVI绘制出两选影数值关系图分选定两波段XY轴面坐标绘制两者散点图
1 影窗口菜单栏中选择Tools → 2D Scatter Plots接着Scatter Plot Band Choice话框会出现屏幕该话框中选择进行拟两影波段
2 选择中波段作X轴波段作Y轴然点击OK生成散点图
3 旦翻开散点图绘制窗口〔图16〕鼠标光标放影窗口中意位置住鼠标左键拖动光标时十字丝光标周围10×10范围素散点图中应点会红色突出显示出
提示：选择Window→Mouse Button Descriptions 显示鼠标键Scatter Plot显示窗口中功
4 影窗口移动鼠标光标观察产生跳跃素〔dancing pixels 〕效果
5 散点图突出显示影窗口中含应波谱值素鼠标光标放散点图窗口点击拖动鼠标中间键
6 散点图菜单栏中选择file→cancel 关闭Scatter Plot 窗口
图16 交互式散点图波段2波段3进行拟
加载幅彩色图
1果屏幕没显示波段列表话框ENVI菜单栏中选择Window →Available Bands list翻开该话框图13
注意：果没翻开影ENVI菜单中选择File→Open Image File然选择翻开影接着该影波段加载波段列表中
2 点击RGB单项选择钮然显示窗口中加载幅彩色影
3 通点击列表中波段名列表中种颜色〔红绿蓝〕选择波段点击列表中波段名时指定RGB单项选择钮会动前选择
4 三种颜色波段应时点击Display＃1拉式菜单钮中选择New Display
5 点击Load RGB钮幅彩色影加载新显示窗口中
链接两显示窗口
两显示窗口链接起进行拟两显示窗口链接起显示窗口中〔滚动缩放等〕进行操作会相链接显示窗口中产生相响应两显示窗口中链接起操作步骤：
1 影窗口菜单中选择Tools→link→Link Displays者影中点击鼠标右键弹出快捷菜单中选择Link DisplaysLink Displays话框会翻开
2 Link Displays话框中点击ok建立链接
3 现尝试中窗口进行滚动者缩放操作观察窗口反响
选择感兴趣区域〔Region of Interest〕
ENVI允许影中定义感兴趣区〔region of interestROIs〕感兴趣区提取分类信息统计生成掩膜操作
1 影窗口中选择Overlay→Region of interest者影中点击鼠标右键弹出菜单话框中选择ROI Tool接着该影显示窗口相应ROI Tool话框会出现屏幕〔图17〕
2 绘制感兴趣区边形ROI_Type中选择ROI_Type类型包括边形矩形圆形等具体绘制步骤：
影窗口中点击鼠标左键建立感兴趣边形第顶点
次点击鼠标左键序选择更边线点点击鼠标右键闭合该边形鼠标中键删新定义点者删整边形次点击鼠标右键固定边形位置
通选择ROI Tool话框顶部相应单项选择钮感兴趣区缩放窗口中定义
完成感兴趣区定义该区域会Available Regions of Interest列表中显示出时显示感兴趣区名称颜色包含素数
3 定义新感兴趣区ROI Tool话框中点击New Region定义效果见图17
图17 定义感兴趣区话框
图18 定义三感兴趣区图
4保存ROI
ROI TOOL窗口菜单：FileSave ROIs…
5合ROI
ROI TOOL窗口菜单：OptionsMerge Regions
6裁剪ROI
ROI TOOL窗口菜单：FileSubset Data via ROIs…
：菜单Basic Tools Subset Data via ROIs…
影进行注记
ENVI软件中图影中添加文符号色标条符号注记
1 幅图进行注记影菜单项选择项中选择Overlay →Annotation接着影窗口相应Annotation：Text话框会出现屏幕〔见图19〕
图19文〔Text〕模式Annotation话框
2 绘制图3D外表相似象进行注记绘制窗口菜单栏中选择Options →Annotation
注记类型：Annotation：Text话框中允许添加类型注记类型object拉菜单中进行选择注记类型文符号矩形椭圆边形折线箭头图例尺三北方图标〔Declination Diagrams〕图图例颜色标注影注记默认情况注记类型文〔Text〕话框中选项控制文注记字体〔Font〕〔Size〕颜色〔color〕位置角度
放置注记：方影窗口中添加文注记：
1 Annotation：Text话框中添加输入文〔for example Center South University of Forestry & Technology 〕
话框相应菜单参数设置中选择文字体颜色然影窗口适宜位置点击鼠标左键接着输入文会显示选点位置图110
图110 带注记影
注意：选择Window →Mouse Button Description 描述注记话框中鼠标键项功
2 鼠标左键拖动文注记圆柄窗口中放置文注记话框中改变相应区域设置值者鼠标左键拖动文符号注记改变注记属性位置
3 完成文注记设置点击鼠标右键锁定注记位置
保存恢复注记：
1 Annotation：Text话框菜单栏中选择File→Save Annotation保存影注记
2 翻开Output Annotation Filename 话框保存影注记该话框中指定保存路径保存文件名该文件扩展名ann
注意：果没影注记保存文件中关闭Annotation：Text话框时注记丧失
3 Annotation：Text话框中选择File→Restore annotation恢复保存注记文件
修改先前放置注记
先前已放置注记进行修改步骤进行修改：
1 Annotation：Text话框中选择Object→SelectionEdit
2 鼠标点击托拽出矩形框包含修改注记象
3 圆柄出现点击拖动注记圆柄修改应属性设置新注记象时样
添加网格
1影中添加公里网格信息影窗口中选择Overlay →Grid line Parameters话框
注意：影添加公里网格时影边框会动添加进
2 Grid line Parameters话框设置公里网线宽颜色公里网间隔修改公里网属性
图111 添加注记网格影
3 Edit Pixel Attributes话框中改变公里网标注网格线矩形边框交叉角颜色宽度公里网间隔完成属性设定点击Edit Pixel Attributes话框中ok钮做更改应影中
4 加公里网满意点击Grid line Parameters话框中Apply钮
保存输出影
添加标注进行图处理图进行保存影保存ENVI影文件格式者保存种通图格式然印者导入软件中直接通印机进行印输出
影保存ENVI影格式
1 影窗口菜单栏种选择File →Save Image File Output Display
To image file 话框会出现屏幕
2 选择24特彩色者8特灰阶进行输出然选择图形选项〔包括注记公里网〕边框设置果注记公里网添加彩色影显示中注记公里网会动列graphics选项中选择注记文件应输出图
3 需单项选择钮影结果输出Memory者file中果选择file输入输出文件名
4 点击ok保存影
实验二 影配准
专题旨介绍ENVI中影进行理校正添加理坐标ENVI进行影影校正
1图文件头文件修改
翻开显示SPOT数
首先翻开软件翻开软件带数库中SPOT数步骤：
1 ENVI菜单中选择file →open image file
2 enter data filename文件选择话框出现选择进入envidata目录bldr_spimg文件
3 点击ok
4 波段列表话框出现点击Grey scale单项选择钮鼠标左键点击相应波段名话框顶部列波段中选中SPOT波段选择波段名显示select Band：字段区域中
5 点击Load Band钮加载幅影新显示窗口中
修改ENVI头文件中图信息
1 波段列表中右键点击bldr_spimg文件名Map info图标弹出快捷菜单中选择Edit Map Information
Edit Map Information话框出现屏幕图21
话框列出ENVI中添加理坐标理信息调整ENVIMagic Pixel〔作图坐标系统起始点〕相应影坐标ENVI相应头文件信息图投影文件中识出图投影元图投影参数计算出影中意元理坐标输入图坐标输入理坐标〔纬度〕
2 点击ProjectionDatum文旁边箭头切换钮显示UTM Zone13 North图投影纬度度坐标ENVI处理程中进行转换
3 点击前DMS者DDEG钮分度分秒〔Degree－Minutes－Second〕十进制度〔Decimal Degrees〕间进行切换
4 点击Cancel推出Edit Map Information话框
图21Edit Map Information话框
5 修改图pixel size信息添加公里网格图标注
6 保存图影窗口中选择file →save image as →image file选择输出路径填写输出文件名称
2 影影配准
利SPOT图校正Landsat TM
翻开显示Landsat TM图
1 ENVI菜单中选择file →open image file
2 enter data filenames话框出现选择进入envidata目录bldr_reg子目录列表中选择bldr_tmimg文件
3 文件选择话框中点击open TM影波段加载波段列表中
4 列表中选择band3点击display#1钮拉式菜单中选择new display
5 点击Load Band钮TMband3波段影加载新显示窗口中
显示光标位置值
翻开显示影窗口滚动窗口者缩放窗口中光标位置信息话框步骤进行操作
1 影窗口菜单栏中选择tools →Cursor Locationvalue
2 影窗口滚动窗口缩放窗口TM影移动鼠标光标
注意坐标素单位出影基素坐标面带理坐标SPOT影
3 选择file →Cancel关闭Cursor Locationvalue话框
开始进行影配准
1 ENVI菜单栏中选择Map →Registration →Select GCPs：Image to Image
2 image to Image Registration话框中点击选择display＃1〔spot影〕作Base Image点击display＃2〔TM影〕作Warp image
图22image to Image Registration话框
3 点击ok启动动配准程序通光标放置两幅影相物点添加单独面控制点
注意：名物点需两张图般选择河流道路交叉点者拐角点交叉点拐角点选择图空间分辨率确定Quickbird空间分辨率061米Quickbird进行图校正时道路交叉点宽度应1米
4 启动Ground control Points Selection话框首先两张图找名物点现SPOT图找点鼠标移动该点样方法TM影找该点
5 两缩放窗口〔Zoom〕中查光标位置十字丝移应位置记录点XY值
图23缩放窗口〔Zoom〕话框
6 Ground Control Points Selection话框base XY文框中分输入286255SPOT影中光标位置移动相应位置点样方法Warp XY文框中分输入133264TM影中光标移动相应位置
7 Ground Control Points Selection话框中点击Add point该面控制点添加列表中点击Show list查面控制点列表见图24然样方式寻找余名物点
图24进行影配准Ground Control Points Selection话框Show list列表
注意：缩放窗口中支持亚元〔subpixel〕级定位缩放例越面控制点精度越面控制点选择通移动光标外直接名物点通鼠标进行点击完成show list话框中旦已选择4面控制点RMS误差会显示出
8 果选择控制点中某点误差删该点重新寻找新点代show list话框中鼠标点击该点选择delete删该点果选择控制满意话通Ground Control Points Selection话框选择Options →Clear All Points掉已选择面控制点
9 假设选择面控制点拟满意图进行校正前控制点信息进行保存具体方法Ground Control Points Selection话框中选择file →Save GCPs to ASCII…选择保存路径文件名控制点信息进行保存备次
图25影配准中Image to Image GCP List话框
校正影
里校正影TM影波段时校正波段影中波段
1 Ground Control Points Selection话框中选择Options →Warp Displayed Band
2 Registration Parameters话框中Warp Method钮菜单中选择RSTResampling钮菜单中选择Nearest Neighbor重采样法
图26 Registration Parameters话框
3 输入文件名bldr_tm1wrp点击ok
4 重复步骤1步骤2RST校正法相应选择BilinearCubic Convolution重采样法
5 结果分输出bldr_tm2wrpbldr_tm3wrp文件中
6 次重复步骤1步骤2次选择次项式Polynomial校正法Cubic convolution重采样法然选择Delaunay 三角网Triangulation校正法相应Cubic convolution重采样法
7 述结果分输出bldr_tm4wrpbldr_tm5wrp文件中
样方法TM影余波段进行校正述步骤区方余波段进行校正需重新找点需保存控制点信息直接调
Ground Control Points Selection话框中选择file →Restore GCPs from ASCII调保存控制点信息图进行校正校正图进行合成
8 合成步骤ENVI菜单栏中选择basic tools →lay stacking通lay stacking选择校正波段进行图层叠加叠加图进行保存保存文件名bldrtm_mimg
拟校正结果
动态链接拟校正结果：
1 波段列表中点击原始SPOT波段影名bldr_spimg然菜单中选择Grey scale显示方式显示图
2波段列表中选择bldr_tm1wrp文件Display＃拉式钮中选择New Display点击Load Band该文件加载新显示框中
3 影窗口中点击鼠标右键选择tools →Link →link displays
4 link displays话框中点击okSPOT影已添加理坐标bldr_tm1wrp影链接起
5 动态链接功带理坐标SPOT影进行拟
6 取消动态链接功选择tools →Link →unlink displays
查图坐标
1 影窗口菜单栏中选择Tools →cursor LocationValue
2 浏览带理坐标数集注意采样方法校正方法数值产生效果
3 选择File →Cancel关闭该话框
3 影图配准
局部逐步演示影图配准处理程许步骤影影配准步骤相似步骤会详细讨带理坐标 SPOT 影中获取图坐标矢量数字线划图〔DLG〕作基准数然基素坐标 Landsat TM 影进行校正匹配相应图数
31翻开显示 Landsat TM 影文件
ENVI 菜单中选择 File → Open Image File Enter Data Filenames话框出现选择进入相应文件夹目录表中lxmbldr_tmimg 文件点击 OKTM 影波段加载波段列表中时幅彩色影加载新显示窗口中
波段列表中点击 Gray Scale 钮选择波段 3点击 Load Band 钮 TM 影第 3 波段加载已翻开显示窗口中
32选择影图配准恢复控制点坐标
ENVI 菜单中选择 Map>Registration>Select GCPsImage to Map果翻开影显示窗口 Image to Map Registration 话框中点击选择包含该灰阶影显示窗口显示号投影列表中选择 UTM Zone 文框中输入 13设置素 30m点击 OK启动配准程序
校正影中光标移动图坐标面点〔幅图者 ENVI 矢量文件中[见局部]读取需图坐标〕 添加单面控制点
Ground Control Points Selection 话框中 E〔东〕 N〔北〕文框中手动输入图坐标然点击 Add Point 添加新面控制点 Ground Control Points Selection 话框中选择File>Restore GCPs from ASCII翻开
lxm PZhunpts 文件
Ground Control Points Selection话框中 点击Show List钮 Image to Map GCP List话框中查影图坐标实际影点预测点坐标 RMS误差
33矢量显示数字线划图〔DLGs〕添加图控制点：
ENVI 菜单中选择 File>Open Vector File>USGS DLG文件选择话框中选择 lxmbldr_rddlg 文件 Import Optional DLG File Parameters 话框中选择 Memory 单项选择钮点击 OK读入需数字线划图〔DLG〕数
矢量列表中高亮选择 ROADS AND TRAILS BOULDER CO 文件点击 Load Selected钮 Load Vector 话框中点击 New Vector Window该矢量加载新矢量显示窗口中
Vector Window #1 窗口中点击拖曳鼠标左键激活十字形光标光标处图坐标会 Vector Window #1 窗口底部列出
影显示窗口中选择 Tools>Pixel Locator输入 402 418然点击Apply影光标移动道路交叉口相应点注意缩放窗口中样获取亚元〔subpixel〕级定位精度
带十字形光标矢量窗口显示出图坐标（图）
矢量窗口中鼠标左键点击拖曳矢量光标十字形光标位需道路交叉口时松开鼠标左键矢量光标放置道路交叉口坐标 4775937444332400〔北纬 40 度 3 分3 秒东105度 15 分45 秒〕 矢量窗口中点击鼠标右键弹出快捷菜单中选择Export Map Location新图坐标会出现 Ground Control Points Selection 话框中 Ground Control Points Selection 话框中点击 Add Point添加该图坐标影素观察 RMS误差变化
34 RST三次卷积校正：
Ground Control Points Selection 话框中选择 Options>Warp File Input Warp Image 话框中高亮选择文件名 zhumuqian_tmimg点击 OK TM 6 波段进行校正出现 Registration Parameters dialog 话框中 Warp Method 选 RST Resampling 设置 Cubic Convolution background值改 255 output file 文框中输入输出文件名 zhumuqiantm_mimg点击 OK开始进行影图校正
35 显示结果进行评价
光标位置值〔Cursor LocationValue〕校正彩色影进行评价
波段列表中 点击Gray Scale单项选择钮 接着点击校正影波段Wrap（Band 3 zhumuqiantm_mimg） Display #拉式菜单钮中选择 New Display点击 Load Band加载幅校正 TM影注意校正影倾斜消 Landsat TM 轨道方影响原时影已带理坐标注意空间分辨率 30 米SPOT 影 10 米分辨率
实验三 图融合
数融合幅影组合单合成影处理程般高空间分辨率全色波段低空间分辨率光谱影进行融合达增强光谱影空间分辨率问题
进行数融合影文件必须含理坐标种情况空间采样会动运行果进行数融合影文件没理坐标影必须覆盖理区域相素影相方位
1 理坐标系统影融合
读取显示ER Mapper影
London TM影SPOT影二进制文件含ER Mapper头文件够动ENVIER Mapper读取程序读取
1 选择file →open external file →IP software →ER Mapper然进入lontmsp子目录选择lon_tmers
2 波段列表中单击RGB color单项选择钮次选取红绿蓝字段应波段点击Load RGB显示幅彩色TM影
3 选择file →open external file →IP software →ER Mapper然进入lontmsp子目录选择lon_spoters文件
4 波段列表中单击grey scale单项选择钮选取pseudo layer波段display拉菜单中选择new display点击load band显示幅灰阶SPOT影
调整影
1 波段列表中单击SPOT影发现空间尺寸2820×1569样方法知道TM影空间尺寸1007×560TM影素28米SPOT影素10米TM影必须28倍率调整产生SPOT影匹配10米元
2 选择Basic tools →Resize Data〔spatialSpectral〕选择lon_tmers单击okresize data parameters话框xfac文框中输入28yfac文框中输入28009选择保存路径文件名lontm点击ok调整TM影者直接samples文框中输入2820lines文框中输入1569选择保存路径文件名lon_tm点击ok调整TM影
图 31 Resize Data 话框
3 显示调整影选择Tools →link →Link Displays调整TM影SPOT全色影链接起动态链接分析拟两幅影
图32TM彩色合成影SPOT全色影
进行HIS数融合
解整处理程选择进行手工数融合第步TM彩色影转换色度饱度数值（称亮度值）〔huesaturationvalue（intensity）〕彩色空间高分辨率SPOT影换数值〔value〕波段拉伸01间满足正确数范围TM影获取色度饱度SPOT影中获取数值进行反转换转回红绿蓝彩色空间程产生出幅输出影中包含TM影中获取颜色信息SPOT影中获取空间分辨率信息
HSV正变换
1 ENVI菜单项选择择Transform →color Transforms →RGB to HSV然选择调整TM数作输入RGB影输入输出文件名点击ok执行变换
2 作灰阶影RGB彩色影显示色度饱度数值影
拉伸SPOT影换TM数值波段
1 ENVI菜单项选择择Basic Tools →Stretch Data单击lon_spot文件然点击ok
2 Data stretching话框Output Data中min文框中输入0Max文框中输入10输入输出文件名单击okSPOT影数拉伸浮点型范围010
图33 Data stretching话框
HSV反变换
1 ENVI菜单中选择Transform →Color Transforms →HSV to RGB选择变换TM影HueSaturation波段作变换HS波段
2 选择拉伸SPOT影作V波段点击okHSV to RGB Parameters话框中输入输出文件名点击ok进行反变换
图34 HSV to RGB Parameters话框
显示结果
1 波段列表中点击RGB Color单项选择钮序点击反变换RGB波段点击Load RGB钮显示幅融合TMSPOT彩色影
2 显示融合影选择Tools →Link →Link Displays融合影调整TM影SPOT全色影链接起动态分析拟影
进行动HSV变换融合
1 ENVI菜单项选择择Transform →Image Sharpening →HSV
2 果调整TM彩色影已显示窗口中Select Input RGB话框中直接选择适宜影显示窗口否Select Input RGB话框中选择Red Layer〞Green Layer〞Blue layer〞应调整TM影波段然点击ok
3 High Resolution Input File话框中选择SPOT影点击ok
4 输入输出文件名：lontmspimgHSV Sharpening Parameters话框中点击ok
结果显示链接拟
1 显示融合彩色影波段列表中选择RGB Color单项选择钮RGB波段中选择融合影中相应波段点击Load RGB
2 通影窗口菜单中选择Tools →Link →Link Displays原始TM彩色影SPOT影动融合彩色影进行链接拟
图35 TM彩色影SPOT影融合彩色影链接拟
2 带理坐标系统影融合
翻开进行融合图
1选择file →open image file 然进入envidata目录双击p123r043ETM文件夹选择p123r043_7p19991224_z49_multihdr点击翻开影出现波段例表中
2 波段列表中点击RGB Color选择R〔4〕G〔5〕B〔3〕三波段进行波段显示
3 选择file →open image file 然进入envidata目录双击p123r043ETM文件夹选择p123r043_7p19991224_z49_panhdr点击翻开影出现波段例表中
4波段列表中点击Grey scale选择p123r043_7p19991224_z49_panhdr点击Display＃1拉菜单点击New display点击Load band
Color Normalized〔Brovey〕sharpening〔标准化彩色变换融合〕
1ENVI菜单项选择择Transform →Image Sharpening →Color Normalized〔Brovey〕
2 select Input for color Bands话框中选择Available Bands List点击ok
图36 select Input for color Bands话框
3 出现波段列表话框中p123r043_7p19991224_z49_multihdr文件中选择三波段应RGB选项点击ok
4 high resolution Input file话框中选择p123r043_7p19991224_z49_panhdr文件波段点击ok
图37 high resolution Input file话框
5 Color Normalized sharpening话框中Resampling重采样方法选择Nearest Neighbor输入输出文件名保存路径点击ok会进行动融合融合结果会出现波段例表中
图38 Color Normalized sharpening话框
6 查波段融合结果波段列表话框中选择RGB color显示方式选择R〔4〕G〔5〕B〔3〕点击Load RGB显示融合图
Gramschmidt spectral sharpening〔GS变换融合〕
1ENVI菜单项选择择Transform →Image Sharpening →Gramschmidt spectral sharpening
2 select low spatial resolution multi band input file话框中选择p123r043_7p19991224_z49_multi文件点击ok
4 high resolution Input file话框中选择p123r043_7p19991224_z49_pan文件波段点击ok
5 Gramschmidt spectral sharpening话框中Resampling重采样方法选择Nearest Neighbor输入输出文件名保存路径点击ok会进行动融合融合结果会出现波段例表中
6 查波段融合结果波段列表话框中选择RGB color显示方式选择R〔4〕G〔5〕B〔3〕点击Load RGB显示融合图
图39 select low spatial resolution multi band input file话框
PC spectral sharpening（成分变换融合）
1ENVI菜单项选择择Transform →Image Sharpening →PC spectral sharpening
2 select low spatial resolution multi band input file话框中选择p123r043_7p19991224_z49_multi文件点击ok
4 high resolution Input file话框中选择p123r043_7p19991224_z49_pan文件波段点击ok
5 PC spectral sharpening话框中Resampling重采样方法选择Nearest Neighbor输入输出文件名保存路径点击ok会进行动融合融合结果会出现波段例表中
6 查波段融合结果波段列表话框中选择RGB color显示方式选择R〔4〕G〔5〕B〔3〕点击Load RGB显示融合图
注意：改变图融合重采样方法否会图显示产生影响图清晰度色彩两方面拟图融合方法图影响
实验四 波段组合计算图增强
专题中ETM数象介绍图处理程中波段组合方式波段间运算方式数拉伸增强处理程
翻开影
1选择file →open image file 然进入envidata目录双击p123r043ETM文件夹选择p123r043_7p19991224_z49_multihdr点击翻开影出现波段例表中
2 波段列表中点击RGB Color选择R〔4〕G〔5〕B〔3〕三波段进行波段显示
3 点击Load RGB幅假彩色图显示影窗口中
波段组合
1 显示影中p123r043_7p19991224_z49_multihdr文件三波段文件6波段里尝试6波段中选取3波段组合拟波段组合间图显示效果
2 波段列表中点击RGB Color选择R〔7〕G〔4〕B〔1〕三波段进行波段显示
图41波段组合间影显示效果
3 波段列表中点击Display＃1拉菜单中选择New display点击Load RGB幅真彩色图显示影窗口中
4 动态链接拟前两幅影效果影窗口中选择tools →link →link displayslink displays话框中点击ok两幅影链接起
5 波段列表中选择余波段组合方式加拟
波段运算
Band math功户提供灵活图处理工具Band math话框中实现波段间加减等运算里p123r043_7p19991224_z49_multihdr文件象进行波段运算
1 菜单栏中选择file →open image file 然进入envidata目录双击p123r043ETM文件夹选择p123r043_7p19991224_z49_multihdr点击翻开影出现波段例表中
2 菜单栏中选择Basic Tools →Band math出现band math话框
图42 Band math话框
3 enter an expression文框中输入需进行波段计算IDL〔Interactive Data Language〕表达式变量代波段名文件名变量名必须字符b〞B〞开头面着5数字字符
注意：例假设想计算三波段均值数学表达式：
〔float（b1）+float（b2）+float（b3）〕3输入文框中该表达式中三变量b1表示第变量b2表示第二变量 b3表示第三变量变量前浮点型字节防止计算时出现字节溢出错误输入效表达式点击Add to list表达式添加话框中点击ok
4 变量赋值variables to Bands Pairings话框中允许输入波段列表中波段赋enter an expression〞文框中包含变量
图43 variables to Bands Pairings话框
值原先实例中变量b1〞variables to Bands Pairings话框中标签Variables used in expression〞文框点击变量B1〞
标签Available Bands List〞列表中点击需波段数集显示变量名注意旦第波段选中具相空间尺寸波段显示波段列表中样方法B2〞B3〞等赋予值
5 输出波段运算结果三变量选定相应波段点击ok运算结果进行输出Variables to Bands Pairings话框中选择file〞Memory〞切换钮运算结果显示波段列表中
提醒：波段均值计算外尝试波段间值加减者法等系列运算波段运算般满足四条件：
1〕必须IDL语言书写波段运算表达式2〕输入波段必须具相维数3〕表达式中变量必须Bn〔bn〕命名4〕结果波段必须输入波段维数相
交互式数拉伸〔Interactive Stretching〕
翻开进行拉伸数p123r043_7p19991224_z49_multihdr影窗口选择Enhance→interactive stretching交互式度拉伸话框中显示输入直方图输出直方图说明前输入数分应拉伸两条垂线〔虚线〕说明前拉伸值值彩色影说直方图颜色选波段颜色致〔系统默认显示红色波段直方图〕交互式直方图窗口底部会列出拉伸类型直方图源浏览幅彩色图绿波段蓝波段直方图点击窗口中G〞B〞切换钮拉伸直方图变换动应图选择Options →Auto Apply般情况点击Apply〞钮变换应图通取消选中Options →Auto Apply关闭动功
图44 交互式度拉伸话框
设置拉伸参数
交互式度拉伸话框中Options〞菜单鼠标交互功拉伸设置参数时浏览交互式直方图窗口中信息
1） 更改拉伸值值
更改拉伸值值意条垂直虚线点击鼠标左键然拖放新位置者直接话框顶部文框中输入需值
提示：文框中指定值时输入需实际值〔例值37值68〕者数百分〔例值2％值98％〕设置数范围外值根输入直方图数范围进行拉伸输入数输出直方图会动更新反映输入直方图中变化应新拉伸显示数分布
2） 锁定拉伸条〔locking stretch bars〕
锁定值值拉伸条〔垂直虚线〕间距离选择Options →locking stretch bars拉伸条间距离锁定时移动两拉伸条解锁定取消选中locking stretch bars
拉伸类型菜单
交互式度拉伸话框中Stretch_type菜单中包含供选择交互式拉伸类型列表
1） 线性度拉伸
线性度拉伸系统默认交互式拉伸线性拉伸值值分设置0255两者间值设置中间线性输出值
交互式度拉伸话框中选择Stretch_type →linear contrast stretch限定输入值鼠标左键移动输入直方图中垂直线〔白色虚线〕需位置者stretch文框中输入需DN值数百分〔例值5％值95％〕重新显示原拉伸选择Options →Reset Stretch
2） 分段线性度拉伸〔piecewise linear〕
分段线性度拉伸通鼠标输入直方图中放置点进行交互限定点间局部采线性拉伸
选择Stretch_type →piecewise linear转换函数〔初始条白色直线〕绘制输入直方图输入直方图位置点击鼠标中键转换函数增加节点绘制线段端点绘制节点标记连接起
图45分段线性度拉伸话框
移动点位置标记点击鼠标左键然拖放新位置删点标记点击鼠标右键手动键入输入输出值选择Options →Edit piecewise linear结果绘制输出直方图中该直方图显示应新拉伸数分布情况点击Apply拉伸应显示数
3〕高斯度拉伸〔Gaussian〕
系统默认Gaussian拉伸均值DN127应0255正负3标准差值进行拉伸选择Stretch_type →Gaussian输入拉伸值值输出直方图条红色曲线显示选择Gaussian变换函数拉伸数分布呈白色叠加显示红色Gaussian函数
注意：手动输入需标准差选择Options →Set Gaussian Stdv点击Apply拉伸应显示数
4〕直方图均衡化度拉伸
动缩放数直方图bin中DN数相均衡选择Stretch_type →Equalization输入直方图显示未修改数分布输出直方图条红色曲线显示均衡化函数拉伸数分布呈白色叠加显示点击Apply拉伸应显示数
5〕方根度拉伸〔Square Root〕
计算输入直方图方根然应线性拉伸选择Stretch_type →Square Root输入直方图显示未修改数分布输出直方图条红色曲线显示方根函数拉伸数分布呈白色叠加显示点击Apply拉伸应显示数
6〕定义度拉伸直方图匹配〔Arbitrary〕
Arbitrary选型允许输出直方图顶部绘制形状直方图图直方图相匹配
选择Stretch_type →Arbitrary输入直方图未修改数分布通点击住拖放鼠标左键Output Histogram窗口绘制输出直方图定义直方图绿色显示点击Apply拉伸应显示数
注意：删直方图点击鼠标中键点击鼠标右键接受绘制输出直方图数统计资料匹配图中输出直方图红色显示绘制直方图匹配数函数白色曲线绘制
7〕直方图匹配
定义度拉伸功幅图直方图幅图直方图进行匹配幅图中获输入输出直方图图表顶部Input Histogram〞Output Histogram〞文标签点击鼠标左键拖放名字拖放定义输出直方图中然释放钮导入直方图红色绘制输出直方图拉伸导入直方图匹配点击Apply拉伸应显示数
8〕定义查找表拉伸〔Defining Look UpTable〕
户定义查找表输入DN值拉伸输出值恢复交互式定义
选择Stretch_type →User Defined LUT选择Options →Edit User Defined LUT出现编辑话框时包含输入DN值应拉伸输出值列表显示Edit User Defined LUT〞标签值反映前拉伸情况值点击进行编辑出现Edit Select item〞文框中时输入需值然回车键
注意：该功启动重新设置查找表初始值点击Reset〞点击ok点击Apply拉伸应显示数
图46 Edit User Defined LUT话框
图增强滤波
图增强功影窗口中enhance拉菜单中线性增强高斯增强直方图均衡增强方根增强种类型种增强类型三选择影窗口卷积窗口缩放窗口分进行增强处理
滤波通常通特定空间频率图增强空间频率通常描述亮度DN值距离方差图包括种空间频率消幅图高频信息图显滑ENVI软件中进行图滤波菜单栏中filter进行操作包括卷积滤波〔Convolutions and Morphology〕形态学滤波纹理滤波适应滤波频率滤波
卷积滤波指空间域图进行滤波形态学滤波形态根底图进行处理纹理滤波图中提取纹理信息适应滤波器抑制噪声时保存图锋利信息细节傅立叶滤波频率域图进行滤波
1） 卷积滤波
菜单中选择filter →Convolutions and Morphology出现Convolutions and Morphology Tool话框
图47 Convolutions and Morphology Tool话框
选择滤波器类型Convolutions and Morphology Tool话框中点击Convolutions拉菜单中选择滤波器类型〔高通低通拉普拉斯直通高斯高通高斯低通中值sobelRoberts户定义滤波〕系统默认滤波器类型高通类型
编辑变换核Convolutions and Morphology Tool话框中鼠标点击Kernel Size〞增减箭头时变换核尺寸2单位增减鼠标中键点击Kernel Size〞 增减箭头时变换核尺寸10单位增减鼠标右键点击Kernel Size〞 增减箭头时变换核尺寸恢复默认尺寸变换核默认设置正方形变换核变非正方形核选择Option →Square Kernel：No默认Image Add Back0
2） 纹理滤波器〔基概率统计滤波〕
Texture选项应基概率统计二阶概率统计纹理滤波许图包含区域亮度变化特征仅局限亮度值纹理指图色调作等级函数空间变化定义纹理清晰区域灰度等级相纹理区定拟接
菜单中选择filter →Texture →occurrence MeasuresTexture input file话框中选择进行滤波文件点击ok
图48 Texture input file话框
Occurrence Texture Parameter话框中选择输出文件路径核文件名余选项采默认值点击ok文件进行纹理滤波处理选纹理图计算出放置波段例表中
图49 Occurrence Texture Parameter话框
3） 适应增强型Lee滤波器
适应滤波运围绕元标准差计算新元值般说原始元值基周围效元〔符合标准差元〕计算新值代增强型Lee滤波器Lee滤波器改良根单独滤波窗口中计算出统计〔方差系数〕数进行滤波元分3类型中：相似元〔homogenous〕差异元〔heterogeneous〕指目标元〔point target〕
菜单中选择filter →Adaptive →enhance LeeEnhance Lee filter Input file话框中选择进行滤波文件点击okEnhance Lee filter Parameters话框中输入filter SizeDamping FactorHomogenous Areas等参数采默认方式参数填写完毕选择输出路径输出文件名点击ok
图410 Enhance Lee filter Parameters话框
4〕快速傅立叶变换滤波器
FFT Filtering〔Fast Fourier Transform Filtering 快速傅立叶变换滤波 〕图变换成显示空间频率成分合成输出图包括操作：频率域滤波器交互式建立滤波器应FFT原始数空间逆变换
正FFT滤波
翻开影ENVI菜单栏中选择Filters →FFT Filtering →Forward FFT出现Forward FFT Input File话框选择进行滤波文件点击ok
Forward FFT Parameters话框中选择输出file〞memory〞点击ok开始FFT计算
定义FFT滤波器
Filter Definition选项交互式直接定义滤波器通显示正变换图中绘制定义滤波器滤波类型包括：circle pass and cutuserdrawn pass and cut
ENVI菜单栏中选择Filters →FFT Filtering →Filter Definition出现Filter Definition选择话框选项中选择：
果前显示中包含幅FFT影选择包含FFT图显示号选择No Display〞指出该滤波器某特定图相关果前显示中含FFT图Samples〞Lines〞文框中键入数值限定滤波器
图411 Filter Definition话框
Filter Definition选择话框中选择Filter_type滤波器类型Circular pass〞Circular cut〞滤波器需Radius〞文框中元单位输入滤波半径
Band Pass〞Band Cut〞滤波器Inner Radius〞Outer Radius〞文框中元单位键入需值
User Defined〞User Defined Cut〞滤波器允许ENVI注记导入滤波器〔般FFT影annotation功〕
箭头钮选择Number of Border Pixels〞参数细化滤波器值0表示没滑
选择输出file〞memory〞点击apply构建FFT滤波器
反FFT变换
反FFT变换程序实际包括两步骤先应FFT滤波然FFT图反变换回原始数空间
选择Filter →FFT Filtering →Inverse FFT出现Inverse FFT Input File话框选择输入文件波段假设需选取意子集点击ok图进行反变换
实验五 图分类
专题介绍遥感影分类操作处理拟监督分类非监督分类影分类影进行相应处理处理包括聚合〔clump〕处理筛选〔sieve〕处理类〔combine〕处理精度评估
翻开显示分类影
1ENVI菜单中选择File →Open Image file选择进入envidata文件夹目录点击can_tm子目录列表中选择can_tmimg文件点击ok
2 波段列表中选择RGB Color单项选择钮然鼠标左键次点击波段4波段5波段3选波段会话框中适文框中显示出
3 点击Load RGB钮该影加载新显示窗口中影窗口影进行简单区分目视判读工作
非监督分类
ENVI菜单中选择Classification →Unsupervised →选择Kmeans者IsoData
K均值〔Kmeans〕分类法
K－均值非监督分类器聚类分析方法需首先选定需分类数机查找聚类中心位置然迭代重新配置直达优化波谱分类
1ENVI菜单中选择Classification →Unsupervised →KmeansClassification Input file话框中选择分类影分类波段数点击ok
2 Kmeans parameters话框中确定分类数6类余参数采默认方式选择文件输入名称输出路径点击ok影进行分类
3 波段列表中添加新图显示窗口Display＃1拉菜单中选择New Display
4 双击Kmeans〔can_tmimg〕分类图分类图显示新影窗口中
5 影窗口菜单中选择Tools →Link →Link Displays然话框中点击ok
图51 Kmeans parameters话框
6 鼠标左键影点击拖动动态叠加显示区域Kmeans分类结果原始彩色合成影进行拟
7 处理完成选择Tools →Link →Unlink Display关闭动态链接
IsoData〔迭代组织数分类技术〕
IsoData非监督分类法计算数空间中均匀分布类均值然距离规剩余元进行迭代聚合次迭代重新计算均值根新均值元进行分类处理程持续类元数变化少选元变化阈值者达迭代次数
1ENVI菜单中选择Classification →Unsupervised →IsoDataClassification Input file话框中选择分类影分类波段数点击ok
2 IsoData parameters话框中确定分类数6类余参数采默认方式选择文件输入名称输出路径点击ok影进行分类
3 波段列表中添加新图显示窗口Display＃1拉菜单中选择New Display
4 双击IsoData〔can_tmimg〕分类图分类图显示新影窗口中
5 影窗口菜单中选择Tools →Link →Link Displays然话框中点击ok
6 鼠标左键影点击拖动动态叠加显示区域IsoData分类结果原始彩色合成影进行拟住鼠标左键时次点击鼠标右键触发第三幅影动态叠加显示IsoData分类结果Kmeans分类结果进行拟
7 处理完成选择Tools →Link →Unlink Display关闭动态链接
监督分类法
监督分类种边训练边学分类方法进行监督分类时首先训练分类样确定分类数样进行保存然选择监督分类方法分类结束进行分类处理工作
1选择非监督分类影can_tmimg文件影窗口中选择Overlay →Region of Interest接着应显示窗口ROI Tool话框
图52 ROI Tool话框
2 ROI Tool话框中选择ROI_typepolygon然影窗口中绘制边形该边形代表新创立感兴趣区完成步请面步骤进行：
l 影窗口中点击鼠标左键建立感兴趣边形第点
l 次点击鼠标左键序选择更边线点点击鼠标右键闭合该边形鼠标中键删新定义点者删已闭合整边形次点击鼠标右键固定边形位置
l 通选ROI Controls话框顶部应单项选择钮感兴趣区缩放窗口中定义
感兴趣区域定义完会话框区域列表中显示出时显示感兴趣区域名字颜色包含素总数
3 定义新感兴趣区点击New Region钮点击Edit钮翻开Edit ROI Parameters话框中输入感兴趣名字选择感兴趣区颜色填充方式述步骤定义新感兴趣区分类影定义6类
图53 感兴趣话框训练样
4 训练样进行保存ROI Tool话框中选择file →save ROIs
5 Save ROIs to file话框中选择Select All Items点击choose选择保存路径文件名
6 点击ok保存文件
图54 Save ROIs to file话框
距离法分类〔Minimum Distance〕
距离分类法感兴趣区均值矢量计算未知元类均值矢量欧氏距离非户指定标准差距离阈值否元分类感兴趣区中接类
1ENVI菜单中选择classification →supervised →Minimum DistanceClassification Input File话框中点击选择进行监督分类文件can_tmimg点击ok
图55 Minimum Distance Parameters话框
2 Minimum Distance Parameters话框中选择Select All Items set Max stdev from Meanset Max Distance Error采默认方式选择输出路径文件名话框右边选择输出规图Output Rule Image钮鼠标点击
3点击ok进行距离监督分类计算结果出现波段列表中
4波段列表中添加新图显示窗口Display＃1拉菜单中选择New Display
5 双击min Dist〔can_tmimg〕分类图分类图显示新影窗口中
注意：监督分类方法采距离分类外采似然法行六面体等分类技术分类步骤距离法致
影分类处理
需分类影进行处理评价分类效果精度者类概括出导入图影矢量GIS中
分类统计
分类统计分类图中提取统计信息统计信息根统计信息〔值值标准差方差等〕直方图者选类中计算出均波谱
1 选择classification →post classification →class Statistics进行统计处理选择分类影然点击ok
2 接着选择原始分类影can_tmimg点击ok
3Class Selection话框选择进行统计类点击Select All Items然点击ok
4 Compute Statistics Parameters话框中选择计算统计信息点击Compute Statistics Parameters话框底部ok钮然根选择统计选项绘制图报表会出现屏幕
混淆矩阵〔Confusion Matrix〕
ENVI计算混淆矩阵功允许两幅图〔分类影真实影〕者分类影感兴趣区进行拟面真实影幅分类影者根面真实测量生成影
1 选择Classification →Post Classification →Confusion Matrix →两选项供选择Using Ground Truth Image者Using Ground Truth ROIs
2 专题中选择Using Ground Truth ROIs进行混淆矩阵计算
3 Classification Input file话框中选择分类影mindistanceimg点击ok
4 Confusion Matrix Parameters话框中点击ok钮
5 查混淆矩阵混淆影通波谱剖面廓线CursorLocation Vlaue分类影原始反射率影中感兴趣区进行拟确定误差源
图56 分类影混淆矩阵结果
聚合筛选处理〔Clump and Sieve〕
1 分类影进行筛选处理选择ENVI菜单栏中classification →Post classification →sieve classes选择分类影输入memory中然点击ok
2 输出筛选处理结果作聚合处理输入图选择classification →Post classification →clump classes存中选择先前生成筛选处理影点击ok
3 输出Memory点击Clump Parameters话框中ok钮
4 拟处理前三幅影
合类〔Combine Classes〕
1选择ENVI菜单栏中classification →Post classification →combine classes选择分类影点击ok
2 combine classification parameters话框中选择进行合两类先选择合类选择新输出类〔列unclassified 入bare类组成新bare类Input class选择unclassifiedOutput class选择bare〕点击ok生成新合类
图57 Combine Classification Parameters话框
修改类颜色
显示分类影时通修改类颜色改变特定类应颜色
1 影窗口中选择tools →Color Mapping →Class Color Mapping
2 Class Color Mapping话框中点击某类名字拖动相应颜色条者输入需颜色值改变类颜色做改变会立刻应分类影进行永久性改变话框中选择Options→ Save Changes
实验六 图制图
专题旨介绍ENVI图制图处理流程注记功添加额外信息ENVI图制图功够方便快捷交互式添加图例等素幅影绘制成图制图程般包括步骤组成首先ENVI快速制图〔QuickMap〕动生成根模板〔者恢复保存模板〕然ENVI注记功影叠合功需进行交互式制图快速制图允许设定图例输出页方位够选择影空间子集进行制图方便添加根图素图公里网例尺图图标标识图投影信息根图注记
翻开进行制图Landsat TM影
1 ENVI菜单栏中选择file →Open Image file
2 出现Enter Input Data File文件选择话框中选择进入envidata目录中ys_tmsub子目录列表中选择ysratioimg文件点击Open相应文件波段会列波段列表中根默认设置573134率波段会动加载R〞G〞B〞字段时选择RGB Color钮
3 点击话框底部Load RGB钮该影加载新显示窗口中
生成快速制图模块
1 影显示窗口菜单中选择File →QuickMap →New QuickMap翻开QuickMap Default layout话框
话框修改输出页页方位图例
图61 QuickMap Default layout话框
2 QuickMap Default layout话框中Map Scale改200000余参数采默认方式点击ok选定影子集
3 QuickMap image Selection话框中鼠标左键点击红色方框左角拖动方框选中整影
注意：选择影子集输出成图鼠标左键放置调整红色矩形框
图62 QuickMap image Selection话框
4 点击okQuickMap Parameters话框会出现屏幕
5 话框中鼠标左键点击Main Title文框键入文Yellowstone National Park Image Map
6 话框中鼠标右键点击Lower left Text文框弹出菜单栏中选择Load Projection InfoENVI头文件中加载影投影信息
7 话框中鼠标右键Lower Right Text文框键入Map Generated Using然回车〔Enter键〕次输入ENVI QuickMapCopyright2001Research System Inc次键入时回车键
图63 QuickMap Parameters话框
8 根专题目里建议ScaleGrid Lines North 前面复选框〔check box〕选中更改设置然QuickMap Parameters话框右边复选框选项进行
9 点击Declination Diagram复选框选中
10 话框底部选择Save Template输入Output Filename文件名ysratioqm点击ok快速制图结果保存快速制图模板文件模板处理相素影时进行调需显示需影选择File →QucikMap →from Previous Template恢复已保存快速制图模板
11 QuickMap Parameters话框左角点击Apply标准ENVI显示窗口中显示快速制图结果果需修改QuickMap Parameters话框中设置然点击Apply更新显示结果
12 ENVI显示窗口中查结果观察图公里网例尺指北针默认文位置
定义图版面设计中元素
ENVI提供种定制图制图选项包括添加虚边框〔Virtual borders〕文注记公里网等值线绘图插入矢量叠合分类叠合ENVI影窗口滚动窗口者缩放窗口进行定制制图操作〔果滚动窗口拖住角调整显示〕
1〕添加虚边框
默认ENVI显示包含影周围没空白空间图制图般需某图素放置影外ENVI提供虚拟边框〞功注记放置影周围创立新影虚拟边框采种方法添加影次介绍
ENVI快速制图中公里网动添加
虚拟边框会添加图影中快速制图中公里网相协调默认公里网显示出参见具体说明解修改公里网属性参数必边框动添加影两边
1 改变默认边框快速制图影显示窗口菜单中选择Overlay →Grid linesGrid line Parameters话框出现选择Options →Set Display Borders翻开Display Borders话框
2 输入100400150100点击ok
图65 设置Display Borders话框
新虚拟边框属性马会应图影中果保存公里网参数文件边框属性信息会时保存公里网参数文件中恢复公里网属性时边框属性信息会恢复
显示参数设置
Display Preference话框改变虚拟边框属性显示设置
1 快速制图影显示窗口菜单中选择File →Preference接着Display Preference话框会出现话框顶部包含刚刚介绍相类似文框
2 输入需数值边框选择需颜色
3 点击OK新边框属性会应图影中
注记功
ENVI中注记功调整虚拟边框属性
1 快速制图影显示窗口菜单中选择Overlay →Annotation
2 Annotation话框出现选择Options →Set Display Borders翻开Display Borders话框
3 输入需边框属性参数点击OK新虚拟边框属性马会应图影中果保存注记文件边框属性信息会时保存注记文件中恢复公里网属性时边框属性信息会恢复
2〕添加公里网
ENVI支持时显示素公里网理坐标公里网理坐标〔纬度度〕网公里网应图影中时100素虚拟边框〔没添加虚拟边框情况〕会动添加图影中协调公里网标签〔labels〕添加修改图影公里网参考步骤进行
1 快速制图影显示窗口菜单中选择Overlay →Grid Lines出现Grid Lines Parameters话框中默认公里网间距会显示出
2 Grid Spacing文框中输入4000公里网间距改4000米改变公里网标签属性参数选择Options →Edit Map Grid Attributes者Edit Geographic Grid Attributes修改选公里网属性样Quick Map Parameters话框中点击Additional Properties钮修改公里网属性参数
3 点击OK确定选属性
4 Grid Line Parameters话框中点击Apply新公里网应图影显示中保存公里网属性便续Grid Parameters菜单栏中选择File →Save Setup选择输出文件样保存公里网属性参数模板便加载图制图中
3〕操作处理注记
注记ENVI显示者图制图中插入〔放置〕图素种常规手段ENVI根注记功种图素放置位置
1 快速制图影显示窗口菜单中选择Overlay →Annotation翻开Annotation话框
2 Annotation话框菜单栏Object拉菜单中选择需注记素
3 选择ImageScroll者Zoom单项选择钮指定注记放置窗口
4 鼠标左键放置注记素点击鼠标右键进行确定锁定注记位置
图67 注记〔Annotation〕话框
Object拉菜单中文〔Text〕符号〔symbol〕形状注记线段箭头注记图例尺〔Scale Bar〕三北方表〔Declination〕图例注记〔Map key〕嵌入式影注记〔Image〕绘制图〔Plot〕嵌入式注记等
专题中求添加注记
①文文档：Landsat TM Dtata
Ratios 573134（RGB）
②图例注记〔Map key〕
Annotation话框菜单栏Object →Map key选择Edit Key Items添加删者修改单图例项修改成功回车键鼠标左键放置图例鼠标右键锁定图例位置
③嵌入影注记
首先确保插入影已波段列表中列出果没波段列表中话先翻开插入影ys_loctif
Annotation话框菜单栏Object →ImageAnnotation：Image话框中选择Select New Image选择插入影ys_loctif
图68 Annotation：Image话框 图610 Annotation：Declination话框
波段列表中选取插入影果需处理中选取空间子集调整影鼠标左键放置嵌入影鼠标右键锁定嵌入影
图611 ENVI注记功定制图制图
注记添加完毕注记文件进行保存Annotation话中选择File →Save Annotation注记文件进行保存
输出图制图结果
1 快速制图影显示窗口菜单中选择File →Save Image As →Image file
2 Out Display to Image file话框中Out file type类型选择BMP者JPEG选择 保存路径文件名点击ok进行保存幅完整图制图完成
文档香网(httpswwwxiangdangnet)户传
《香当网》用户分享的内容，不代表《香当网》观点或立场，请自行判断内容的真实性和可靠性！
该内容是文档的文本内容，更好的格式请下载文档
