erdas遥感图像基本处理傻瓜式步骤
ERDAS的操作手册
ERDAS的操作手册纠正,融合,镶嵌是遥感处理中比较常见的三种处理方法。
对于初学遥感的人来说,掌握这三种方法是十分必要的。
下面,我们通过一些实例,在ERDAS 中的操作,来分别介绍这三种处理方法。
1、纠正纠正又叫几何校正,就是将图像数据投影到平面上,使其符合地图投影系统的过程;而将地图坐标赋予图像数据的过程,称为地理参考(Geo-referencing)由于所有地图投影系统都遵从于一定的地图坐标系统,所以几何校正包含了地图参考。
(1)启动在ERDAS中启动几何校正有三种方法:A、菜单方式B、图标方式C、窗口栅格操作窗口启动这种方法比较常用,启动之前在窗口中打开需要纠正的图像,然后在栅格操作菜单中启动几何校正模块。
建议使用这种启动方法,更直观简便。
(2)设置几何校正模型常用模型:功能Affine 图像仿射变换(不做投影变换)Camera 航空影像正射校正Landsat Landsat卫星影像正射校正Polynomial 多项式变换(同时做投影变换)Rubber Sheeting 非线性、非均匀变换Spot Spot卫星图像正射校正其中,多项式变换(Polynomial)在卫星图像校正过程中应用较多,在调用多项式模型时,需要确定多项式的次方数(Order),通常整景图象选择3次方。
次方数与所需的最少控制点数是相关的,最少控制点数计算公式为((t+1)*(t+2))/2,公式中t为次方数,即1次方最少需要3个控制点,2次方需要6个控制点,3次方需要10个控制点,依此类推。
(3)几何校正采点模式A、Viewer to Viewer 已经拥有需要校正图像区域的数字地图、或经过校正的图像,就可以采用Viewer to Viewer的模式。
B、File to Viewer 事先已经通过GPS测量、或摄影测量、或其它途径获得了控制点坐标,并保存为ERDAS IMAGINE的控制点格式或ASCII数据文件,就可以采用File to Viewer模式,直接在数据文件中读取控制点坐标。
erdas操作-掩膜、BCF批处理、平均光谱值
4.
需要准备掩膜文件img格式,待分析区域像素值为1,其余区域为0.该文件可用model maker生成。
掩膜输出文件中,需主要勾选忽略0值
5.
1.一般各类功能在设置好各类参数之后,都有batch的按钮
选择第三个,自动执行命令
点next,单击variables按钮,选择输入、输出文件的选择方式(一般output文件的文件名和路径可以由程序自动设定)
E
1.
在erdas的Viewer中加载shp文件。
Viewer中鼠标点击shp文件,该文件变黄色高亮。
在Viewer菜单栏中单击AOI,在AOI下拉菜单中选择Copy Selection to AOI
在Viewer菜单栏中单击File,在File下拉菜单中选择Save->AOI Layer as
保存
#
SET CELLSIZE MIN;
#
# set window for the model
#
SET WINDOW INTERSECTION;
#
# set area of interest for the model
#
SET AOI NONE;
#
# declarations
#
Integer RASTER n1_b1 FILE OLD PUBINPUT NEAREST NEIGHBOR AOI NONE "d:/temp/b1.img";
#
# function definitions
#
n6_layerstack = STACKLAYERS ( $n1_b1 , ($n2_b2 + $n4_b4 / 3) / 2 , $n3_b3) ;
ERDAS详细操作
目录1. 影像阅读2. 遥感影像分幅裁剪与拼接处理3. 影像几何校正及正射影像制作4. 影像增强1. 影像阅读1.1 设置erdas的各种默认参数1)在ERDAS IMAGINE的主菜单栏上找到sessio n→Preferences,单击出现Preferences editor对话框。
2)通过拖动Category的滚动条,可以看到右方对应出现的各个参数,同时也可以在文本编辑处修改这些参数。
3)在Category下选择Viewer,拖动滚动条查看它的各种参数。
4)查看Category的帮助信息。
点击右下方的“help”和“Category Help”,则出现以下的界面,如果有不懂的地方我们就可以通过这个帮助信息寻求答案。
1.2 显示图像1)在ERDAS主菜单上点击图标,新建一个经典窗口,如下图:2)在Viewer界面上点击File→Open →Raster Layer,在默认路径中打开lanier.im g。
3)点击Raster Options栏设置图层的红绿蓝三个波段的分配。
将原来的4 3 2 改为4 5 3后,图象的色调明显变化了。
1.3 查询像素信息1)使用查询功能选择Utility→Inquire Cursor出现下图中的对话框,通过左下方的四个三角形的符号来分别调整查询指针的上下左右的位置,圆圈表示使查询指针回到中心处,指针的移动,其中的X和Y坐标的数值也会跟着作相应的变化。
指针所指的像素的信息被显示在单元格里。
2)改变查询指针的类型选择Utility→Inquire Color,选择为黄色,则查询指针的十字框的颜色由白色变为了黄色。
3)改变查询指针的形状选择Utility→Inquire Shape,呈现的滚动条列表中选择circle.cursor,再点击Use Cursor button, 然后点击Apply。
4)量测通过这个工具可以实现在所在图层中的点,线,面,矩形,椭圆形的长度(周长)和面积。
ERDAS实验详细操作步骤
目录1. 影像阅读2. 遥感影像分幅裁剪与拼接处理3. 影像几何校正及正射影像制作4. 影像增强1. 影像阅读1.1 设置erdas的各种默认参数1)在ERDAS IMAGINE的主菜单栏上找到sessio n→Preferences,单击出现Preferences editor对话框。
2)通过拖动Category的滚动条,可以看到右方对应出现的各个参数,同时也可以在文本编辑处修改这些参数。
3)在Category下选择Viewer,拖动滚动条查看它的各种参数。
4)查看Category的帮助信息。
点击右下方的“help”和“Category Help”,则出现以下的界面,如果有不懂的地方我们就可以通过这个帮助信息寻求答案。
1.2 显示图像1)在ERDAS主菜单上点击图标,新建一个经典窗口,如下图:2)在Viewer界面上点击File→Open →Raster Layer,在默认路径中打开lanier.im g。
3)点击Raster Options栏设置图层的红绿蓝三个波段的分配。
将原来的4 3 2 改为4 5 3后,图象的色调明显变化了。
1.3 查询像素信息1)使用查询功能选择Utility→Inquire Cursor出现下图中的对话框,通过左下方的四个三角形的符号来分别调整查询指针的上下左右的位置,圆圈表示使查询指针回到中心处,指针的移动,其中的X和Y坐标的数值也会跟着作相应的变化。
指针所指的像素的信息被显示在单元格里。
选择Utility→Inquire Color,选择为黄色,则查询指针的十字框的颜色由白色变为了黄色。
选择Utility→Inquire Shape,呈现的滚动条列表中选择circle.cursor,再点击Use Cursor button, 然后点击Apply。
4)量测通过这个工具可以实现在所在图层中的点,线,面,矩形,椭圆形的长度(周长)和面积。
1.4 图层管理1)排列多个图层打开lnsoils.img,并在Raster Options栏中将锁定的ClearDisplay项取消,即不清除视窗中已经打开的图像,使得多个图层能够在一个窗口中存在和切换。
ERDASIMAGINE遥感图像处理教程
《ERDAS IMAGINE遥感图像处理教程》根据作者多年遥感应用研究和ERDAS IMAGINE软件应用经验编著而成,系统地介绍了ERDAS IMAGINE 9.3的软件功能及遥感图像处理方法。
全书分基础篇和扩展篇两部分,共25章。
基础篇涵盖了视窗操作、数据转换、几何校正、图像拼接、图像增强、图像解译、图像分类、子像元分类、矢量功能、雷达图像、虚拟GIS、空间建模、命令工具、批处理工具、图像库管理、专题制图等ERDAS IMAGINE Professional级的所有功能,以及扩展模块Subpixel、Vector、OrthoRadar、VirtualGIS等;扩展篇则主要针对ERDAS IMAGINE 9.3的新增扩展模块进行介绍,包括图像大气校正(ATCOR)、图像自动配准(AutoSync)、高级图像镶嵌(MosaicPro)、数字摄影测量(LPS)、三维立体分析(Stereo Analyst)、自动地形提取(Automatic Terrain Extraction)、面向对象信息提取(Objective)、智能变化检测(DeltaCue)、智能矢量化(Easytrace)、二次开发(EML)等十个扩展模块的功能。
《ERDAS IMAGINE遥感图像处理教程》将遥感图像处理的理论和方法与ERDAS IMAGINE软件功能融为一体,可以作为ERDAS IMAGINE软件用户的使用教程,对其他从事遥感技术应用研究的科技人员和高校师生也有参考价值。
基础篇第1章概述21.1 遥感技术基础21.1.1 遥感的基本概念21.1.2 遥感的主要特点21.1.3 遥感的常用分类31.1.4 遥感的物理基础31.2 ERDAS IMAGINE软件系统6 1.2.1 ERDAS IMAGINE概述6 1.2.2 ERDAS IMAGINE安装7 1.3 ERDAS IMAGINE图标面板11 1.3.1 菜单命令及其功能111.3.2 工具图标及其功能141.4 ERDAS IMAGINE功能体系14第2章视窗操作162.1 视窗功能概述162.1.1 视窗菜单功能172.1.2 视窗工具功能172.1.3 快捷菜单功能182.1.4 常用热键功能182.2 文件菜单操作192.2.1 图像显示操作202.2.2 图形显示操作222.3 实用菜单操作232.3.1 光标查询功能232.3.2 量测功能242.3.3 数据叠加显示252.3.4 文件信息操作272.3.5 三维图像操作292.4 显示菜单操作332.4.1 文件显示顺序332.4.2 显示比例操作332.4.3 显示变换操作342.5 AOI菜单操作342.5.1 打开AOI工具面板35 2.5.2 定义AOI显示特性35 2.5.3 定义AOI种子特征35 2.5.4 保存AOI数据层36 2.6 栅格菜单操作372.6.1 栅格工具面板功能37 2.6.2 图像对比度调整392.6.3 栅格属性编辑402.6.4 图像剖面工具432.7 矢量菜单操作452.7.1 矢量工具面板功能46 2.7.2 矢量文件生成与编辑472.7.3 改变矢量要素形状482.7.4 调整矢量要素特征482.7.5 编辑矢量属性数据492.7.6 定义要素编辑参数502.8 注记菜单操作502.8.1 创建注记文件512.8.2 设置注记要素类型522.8.3 放置注记要素522.8.4 注记要素属性编辑542.8.5 添加坐标格网55第3章数据输入/输出563.1 数据输入/输出概述563.2 二进制图像数据输入573.2.1 输入单波段数据573.2.2 组合多波段数据583.3 其他图像数据输入/输出59 3.3.1 HDF图像数据输入操作59 3.3.2 JPG图像数据输入/输出60 3.3.3 TIFF图像数据输入/输出61第4章数据预处理624.1 遥感图像处理概述62 4.1.1 遥感图像几何校正62 4.1.2 遥感图像裁剪与镶嵌63 4.1.3 数据预处理模块概述63 4.2 三维地形表面处理64 4.2.1 启动三维地形表面64 4.2.2 定义地形表面参数65 4.2.3 生成三维地形表面66 4.2.4 显示三维地形表面67 4.3 图像几何校正674.3.1 图像几何校正概述67 4.3.2 资源卫星图像校正70 4.3.3 遥感图像仿射变换76 4.3.4 航空图像正射校正78 4.4 图像裁剪处理814.4.1 图像规则裁剪814.4.2 图像不规则裁剪824.4.3 图像分块裁剪844.5 图像镶嵌处理844.5.1 图像镶嵌功能概述84 4.5.2 卫星图像镶嵌处理90 4.5.3 航空图像镶嵌处理934.6 图像投影变换954.6.1 启动投影变换954.6.2 投影变换操作964.7 其他预处理功能964.7.1 生成单值栅格图像964.7.2 重新计算图像高程974.7.3 数据发布与浏览准备97 4.7.4 产生或更新图像目录98 4.7.5 图像范围与金字塔计算99第5章图像解译1005.1 图像解译功能概述1005.1.1 图像空间增强1005.1.2 图像辐射增强1015.1.3 图像光谱增强1015.1.4 高光谱基本工具1025.1.5 高光谱高级工具1035.1.6 傅里叶变换1035.1.7 地形分析功能1045.1.8 地理信息系统分析104 5.1.9 实用分析功能1055.2 空间增强处理1065.2.1 卷积增强处理106 5.2.2 非定向边缘增强107 5.2.3 聚焦分析1085.2.4 纹理分析1095.2.5 自适应滤波1105.2.6 统计滤波1115.2.7 分辨率融合1115.2.8 改进IHS融合112 5.2.9 HPF图像融合114 5.2.10 小波变换融合115 5.2.11 删减法融合1165.2.12 Ehlers图像融合117 5.2.13 锐化增强处理118 5.3 辐射增强处理1205.3.1 查找表拉伸1205.3.2 直方图均衡化120 5.3.3 直方图匹配1215.3.4 亮度反转处理122 5.3.5 去霾处理1235.3.6 降噪处理1235.3.7 去条带处理1245.4 光谱增强处理1245.4.1 主成分变换1245.4.2 主成分逆变换1255.4.3 独立分量分析1265.4.4 去相关拉伸1275.4.5 缨帽变换1275.4.6 色彩变换1295.4.7 色彩逆变换1295.4.8 指数计算1305.4.9 自然色彩变换1315.4.10 ETM反射率变换131 5.4.11 光谱混合器1335.5 高光谱基本工具135 5.5.1 自动相对反射1355.5.2 自动对数残差1365.5.3 归一化处理1365.5.4 内部平均相对反射137 5.5.5 对数残差1375.5.6 数值调整1385.5.7 光谱均值1395.5.8 信噪比功能1395.5.9 像元均值1405.5.10 光谱剖面1415.5.11 光谱数据库1425.6 高光谱高级工具142 5.6.1 异常探测1425.6.2 目标探测1475.6.3 地物制图1495.6.4 光谱分析工程向导153 5.6.5 光谱分析工作站154 5.7 傅里叶变换1565.7.1 快速傅里叶变换156 5.7.2 傅里叶变换编辑器157 5.7.3 傅里叶图像编辑158 5.7.4 傅里叶逆变换1685.7.5 傅里叶显示变换169 5.7.6 周期噪声去除1695.7.7 同态滤波1705.8 地形分析1715.8.1 坡度分析1715.8.2 坡向分析1715.8.3 高程分带1725.8.4 地形阴影1735.8.5 彩色地势1735.8.6 地形校正1755.8.7 栅格等高线1755.8.8 点视域分析1765.8.9 路径视域分析181 5.8.10 三维浮雕1825.8.11 高程转换1835.9 地理信息系统分析184 5.9.1 邻域分析1845.9.2 周长计算1865.9.3 查找分析1865.9.4 指标分析1875.9.5 叠加分析1885.9.6 矩阵分析1895.9.7 归纳分析1905.9.8 区域特征1905.10 实用分析功能191 5.10.1 变化检测1915.10.2 函数分析1925.10.3 代数运算1925.10.4 色彩聚类1935.10.5 高级色彩聚类194 5.10.6 数值调整1955.10.7 图像掩膜1965.10.8 图像退化197 5.10.9 去除坏线197 5.10.10 投影变换198 5.10.11 聚合处理199 5.10.12 形态学计算199第6章图像分类202 6.1 图像分类简介202 6.1.1 非监督分类202 6.1.2 监督分类2036.1.3 专家系统分类206 6.2 非监督分类2086.2.1 获取初始分类209 6.2.2 调整分类结果210 6.3 监督分类2126.3.1 定义分类模板213 6.3.2 评价分类模板221 6.3.3 执行监督分类226 6.3.4 评价分类结果227 6.4 分类后处理2316.4.1 聚类统计2326.4.2 过滤分析2326.4.3 去除分析2336.4.4 分类重编码2336.5 专家分类器2346.5.1 知识工程师2356.5.2 变量编辑器2396.5.3 建立知识库2426.5.4 知识分类器248第7章子像元分类2517.1 子像元分类简介2517.1.1 子像元分类的基本特征251 7.1.2 子像元分类的基本原理252 7.1.3 子像元分类的应用领域253 7.1.4 子像元分类模块概述254 7.2 子像元分类方法2567.2.1 子像元分类流程2567.2.2 图像质量确认2587.2.3 图像预处理2597.2.4 自动环境校正2607.2.5 分类特征提取2637.2.6 分类特征组合2697.2.7 分类特征评价2717.2.8 感兴趣物质分类2747.2.9 分类后处理2767.3 子像元分类实例2777.3.1 图像预处理2777.3.2 自动环境校正2777.3.3 分类特征提取2787.3.4 感兴趣物质分类2797.3.5 查看验证文件2817.3.6 分类结果比较282第8章矢量功能2838.1 空间数据概述2838.1.1 矢量数据2838.1.2 栅格数据2848.1.3 矢量和栅格数据结构比较285 8.1.4 矢量数据和栅格数据转换286 8.2 矢量模块功能简介2898.3 矢量图层基本操作2898.3.1 显示矢量图层2898.3.2 改变矢量特性2908.3.3 改变矢量符号2918.4 要素选取与查询2988.4.1 查看选择要素属性2988.4.2 多种工具选择要素2998.4.3 判别函数选择要素3008.4.4 显示矢量图层信息3028.5 创建矢量图层3038.5.1 创建矢量图层的基本方法303 8.5.2 由ASCII文件创建点图层307 8.5.3 镶嵌多边形矢量图层3088.5.4 创建矢量图层子集3108.6 矢量图层编辑3118.6.1 编辑矢量图层的基本方法311 8.6.2 变换矢量图层3138.6.3 产生多边形Label点3148.7 建立拓扑关系3148.7.1 Build矢量图层3158.7.2 Clean矢量图层3158.8 矢量图层管理3168.8.1 重命名矢量图层3168.8.2 复制矢量图层3178.8.3 删除矢量图层3178.8.4 导出矢量图层3188.9 矢量与栅格转换3188.9.1 栅格转换矢量3188.9.2 矢量转换栅格3208.10 表格数据管理3228.10.1 INFO表管理3228.10.2 区域属性统计3288.10.3 属性转换为注记329 8.11 Shapefile文件操作331 8.11.1 重新计算高程3318.11.2 投影变换操作332第9章雷达图像处理3349.1 雷达图像处理基础334 9.1.1 雷达图像增强处理334 9.1.2 雷达图像几何校正336 9.1.3 干涉雷达DEM提取336 9.2 雷达图像模块概述337 9.3 基本雷达图像处理337 9.3.1 斑点噪声压缩3389.3.2 边缘增强处理3409.3.3 雷达图像增强3419.3.4 图像纹理分析3449.3.5 图像亮度调整3459.3.6 图像斜距调整3469.4 正射雷达图像校正3479.4.1 正射雷达图像校正概述347 9.4.2 地理编码SAR图像3489.4.3 正射校正SAR图像3529.4.4 GCP正射较正SAR图像355 9.4.5 比较OrthoRadar校正效果358 9.5 雷达像对DEM提取3599.5.1 雷达像对DEM提取概述359 9.5.2 雷达立体像对数据准备359 9.5.3 立体像对提取DEM工程360 9.6 干涉雷达DEM提取3699.6.1 干涉雷达DEM提取概述369 9.6.2 干涉雷达图像数据准备369 9.6.3 干涉雷达DEM提取工程370 9.6.4 DEM高程生成3759.7 干涉雷达变化检测3769.7.1 干涉雷达变化检测模块376 9.7.2 干涉雷达变化检测操作377第10章虚拟地理信息系统381 10.1 VirtualGIS概述38110.2 VirtualGIS视窗38210.2.1 启动VirtualGIS视窗382 10.2.2 VirtualGIS视窗功能382 10.3 VirtualGIS工程38510.3.1 创建VirtualGIS工程385 10.3.2 编辑VirtualGIS视景387 10.4 VirtualGIS分析39110.4.1 洪水淹没分析39110.4.2 矢量图形分析39410.4.3 叠加文字注记39610.4.4 叠加三维模型39810.4.5 模拟雾气分析40510.4.6 威胁性与通视性分析406 10.4.7 立体视景操作40910.4.8 叠加标识图像41010.4.9 模拟云层分析41210.5 VirtualGIS导航41410.5.1 设置导航模式41410.5.2 VirtualGIS漫游415 10.6 VirtualGIS飞行41610.6.1 定义飞行路线41710.6.2 编辑飞行路线41910.6.3 执行飞行操作42010.7 三维动画制作42010.7.1 三维飞行记录42110.7.2 三维动画工具42210.8 虚拟世界编辑器42210.8.1 虚拟世界编辑器简介422 10.8.2 创建一个虚拟世界425 10.8.3 虚拟世界的空间操作429 10.9 空间视域分析43110.9.1 视域分析数据准备431 10.9.2 生成多层视域数据432 10.9.3 虚拟世界视域分析434 10.10 设置VirtualGIS默认值436 10.10.1 默认值设置环境436 10.10.2 默认值设置选项436 10.10.3 保存默认值设置439第11章空间建模工具44011.1 空间建模工具概述44011.1.1 空间建模工具的组成440 11.1.2 图形模型的基本类型441 11.1.3 图形模型的创建过程44111.2 模型生成器功能组成442 11.2.1 模型生成器菜单命令442 11.2.2 模型生成器工具图标443 11.2.3 模型生成器工具面板444 11.3 空间建模操作过程444 11.3.1 创建图形模型44411.3.2 注释图形模型44711.3.3 生成文本程序44811.3.4 打印图形模型44911.4 条件操作函数应用450第12章图像命令工具453 12.1 图像信息管理技术453 12.1.1 图像金字塔45312.1.2 图像世界文件45312.2 图像命令工具概述454 12.3 图像命令功能操作455 12.3.1 改变栅格图像类型455 12.3.2 计算图像统计值456 12.3.3 图像金字塔操作457 12.3.4 图像地图模式操作458 12.3.5 图像地图投影操作45912.3.6 图像高程信息操作45912.3.7 图像文件常规操作461第13章批处理操作46213.1 批处理功能概述46213.2 批处理系统设置46213.3 批处理操作过程46313.3.1 单文件单命令批处理463 13.3.2 多文件单命令立即批处理465 13.3.3 多文件单命令随后批处理467 13.3.4 多文件多命令批处理469第14章图像库管理47314.1 图像库管理概述47314.2 图像库环境设置47314.3 图像库功能介绍47414.3.1 打开默认图像库47414.3.2 图像库管理功能47514.3.3 图像库图形查询476第15章地图编制47915.1 地图编制概述47915.1.1 地图编制工作流程479 15.1.2 地图编制模块概述479 15.2 地图编制操作过程48015.2.1 准备制图数据48015.2.2 创建制图文件48015.2.3 确定地图制图范围481 15.2.4 放置整饰要素48215.2.5 地图打印输出48915.3 制图文件路径编辑48915.4 系列地图编制工具49015.4.1 准备系列地图编辑文件490 15.4.2 启动系列地图编辑工具490 15.4.3 显示系列地图分幅信息491 15.4.4 系列地图输出编辑491 15.4.5 保存系列地图文件492 15.4.6 系列地图输出预览492 15.5 地图数据库工具492扩展篇第16章图像大气校正49616.1 大气校正模块概述49616.1.1 ATCOR模块主要特征49616.1.2 ATCOR模块功能组成497 16.2 太阳位置的计算49716.3 ATCOR2工作站49816.3.1 ATCOR2工程文件498 16.3.2 光谱分析模块50216.3.3 常数大气模块50716.3.4 增值产品模块51016.4 ATCOR3工作站51116.4.1 ATCOR3生成地形512 16.4.2 ATCOR3工程文件512 16.4.3 光谱分析模块51516.4.4 常数大气模块51516.4.5 增值产品模块516第17章图像自动配准51817.1 图像自动配准模块概述518 17.2 地理参考配准51917.2.1 准备图像数据51917.2.2 产生自动配准点52017.2.3 选择几何模型52317.2.4 定义投影类型52517.2.5 确定输出图像52617.3 图像边缘匹配52717.3.1 准备输入图像52717.3.2 产生自动匹配点52817.3.3 定义匹配策略52817.3.4 选择投影类型52917.3.5 确定输出图像52917.4 自动配准工程52917.4.1 保存自动配准工程文件529 17.4.2 打开自动配准工程文件530 17.5 自动配准工作站53017.5.1 自动配准工作站功能概述530 17.5.2 自动配准工作站应用流程535第18章高级图像镶嵌54318.1 高级图像镶嵌功能概述543 18.1.1 MosaicPro模块特点543 18.1.2 MosaicPro启动过程543 18.1.3 MosaicPro视窗功能544 18.2 高级图像镶嵌工作流程546 18.2.1 航空图像镶嵌54618.2.2 卫星图像镶嵌55318.2.3 图像匀光处理559第19章数字摄影测量56119.1 数字摄影测量基本原理561 19.1.1 数字摄影测量处理过程561 19.1.2 数字图像的内定向56219.1.3 图像核线数字相关56319.1.4 建立规则格网DEM 56319.1.5 图像正射校正处理56419.2 LPS工程管理器56619.2.1 LPS工程管理器功能概述566 19.2.2 LPS工程管理器视窗组成567 19.3 摄影图像摄影测量处理570 19.3.1 摄影图像处理流程57019.3.2 创建LPS工程文件57019.3.3 向LPS工程加载图像572 19.3.4 定义摄影相机几何模型573 19.3.5 定义地面控制点与检查点576 19.3.6 图像同名点自动量测583 19.3.7 执行航空三角测量58419.3.8 图像正射校正处理58719.4 数码图像摄影测量处理588 19.4.1 数码图像处理流程58919.4.2 创建LPS工程文件58919.4.3 向LPS工程加载图像59219.4.4 定义数码相机几何模型59319.4.5 自动量测图像同名点59519.4.6 执行航空三角测量59719.4.7 图像正射校正处理60019.5 扫描图像摄影测量处理60119.5.1 扫描图像处理流程60119.5.2 创建LPS工程文件60119.5.3 加载并定义第一幅图像60319.5.4 加载并定义第二幅图像61019.5.5 图像同名点自动量测61319.5.6 执行空间三角测量61419.5.7 图像正射校正处理616第20章三维立体分析61820.1 三维立体分析基本原理61820.1.1 基于立体像对的高程模型提取618 20.1.2 三维场景重建的实现方法619 20.2 三维立体分析模块概述62120.2.1 三维立体分析模块特点62220.2.2 三维立体分析模块功能62220.3 创建非定向数字立体模型62320.3.1 启动三维立体分析模块62320.3.2 加载三维立体分析图像62320.3.3 调整图像显示参数62520.3.4 保存三维立体模型627第21章自动地形提取62821.1 LPS自动地形提取概述62821.1.1 DTM及其自动提取方法62821.1.2 LPS自动地形提取功能62821.1.3 LPS自动地形提取过程62921.2 LPS自动地形提取操作63021.2.1 创建LPS工程文件63021.2.2 DTM提取参数设置63221.2.3 DTM提取选项设置63321.2.4 DTM自动提取和检查637第22章面向对象的信息提取63922.1 面向对象的信息提取简介63922.1.1 IMAGINE Objective框架设计639 22.1.2 IMAGINE Objective关键特征639 22.2 道路信息提取64022.2.1 道路信息提取模型640 22.2.2 道路信息提取过程640第23章智能变化检测64323.1 智能变化检测原理64323.1.1 图像预处理64323.1.2 变化检测方法64323.1.3 变化定量分析64823.2 智能变化检测应用特点649 23.2.1 智能变化检测技术特征649 23.2.2 智能变化检测工作特点650 23.3 智能变化检测应用操作650 23.3.1 智能变化检测向导模式651 23.3.2 智能变化检测图像显示653 23.3.3 智能变化检测场地检测656第24章智能矢量化65724.1 智能矢量化模块概述657 24.1.1 模块的关键特征65724.1.2 模块的局限性65724.2 智能矢量化模块应用658 24.2.1 模块操作快捷键65824.2.2 启动智能矢量化模块658 24.2.3 跟踪线状地物中心线659 24.2.4 跟踪面状地物边界线662 24.3 智能矢量化模块使用技巧663第25章二次开发工具66525.1 二次开发宏语言EML概述665 25.2 编写EML二次开发程序666 25.2.1 编写EML程序的过程666 25.2.2 执行EML程序的过程667 25.2.3 丰富EML程序的功能667 25.3 EML接口C程序开发包671参考文献674。
Erdas遥感图像处理
ERDAS 遥感图像处理
4.1 图像增强
空间增强处理 辐射增强处理 光谱增强处理 高光谱工具 傅里叶变换
.
23
ERDAS 遥感图像处理
图像辐射增强处理
辐射增强(Radiometric Enhance-ment)技 术是对单个像元的灰度值进行变换达到图 像增强的目的。
• 查找表拉伸
• 直方图均衡化
傅里叶变换——契形滤波
.
54
ERDAS 遥感图像处理
傅里叶变换——组合滤波
.
55
ERDAS 遥感图像处理
傅里叶变换——滤波效果
.
56
ERDAS 遥感图像处理
傅里叶变换——显示变换
傅里叶图像变为ERDAS图像
.
57
ERDAS 遥感图像处理
傅里叶变换——周期噪声去除
通过傅里叶变化消除遥感图像中诸如条带等周期性噪声
通过处理抑制一些无用的信息,增强人或
机器对某些信息的辨别能力。
.
3
ERDAS 遥感图像处理
4.1 图像增强
空间增强处理 辐射增强处理 光谱增强处理 高光谱工具 傅里叶变换
.
4
ERDAS 遥感图像处理
4.1 图像增强
空间增强处理 辐射增强处理 光谱增强处理 高光谱工具 傅里叶变换
.
49
ERDAS 遥感图像处理
滤波半径,低通滤波, 圆外的区域被处理; 高通滤波,圆内的区 域被处理。
.
50
ERDAS 遥感图像处理
傅里叶变换——低通滤波
.
51
ERDAS 遥感图像处理
傅里叶变换——傅立叶逆变换
.
52
ERDAS 遥感图像处理
ERDAS Imagine遥感图像处理方法
Hale Waihona Puke A11ERDAS Imagine遥感图像处理方法
1.3 ERDAS8.4功能体系
输入
栅格图像数据
矢量图形数据
文本属性数据
ERDAS IMAGINE遥感图像处理系统
视窗操作 空间建模 命令工具 批处理 图像库管理
数据输入输出
ERDAS Imagine遥感图像处理方法
ERDAS 8.4
遥感图像处理方法
A
1
ERDAS Imagine遥感图像处理方法
主要内容:
➢ ERDAS 8.4系统简介
➢ ERDAS 8.4视窗操作
➢ ERDAS 8.4数据处理
➢ ERDAS 8.4矢量功能
➢ ERDAS 8.4虚拟GIS
A
2
ERDAS Imagine遥感图像处理方法
A
4
ERDAS Imagine遥感图像处理方法
Essential级
完成二维/三维显示、数据输 入、排序与管理、地图配准、专题 制图以及简单的分析。
可扩充的模块:
(1)Vector模块 (2)Virtual GIS模块 (3) Developer’s Toolkit 模 块
A
5
ERDAS Imagine遥感图像处理方法
1.2 ERDAS8.4图标面板
图标
命令 IMAGINE Credits Start IMAGINE Viewer Import / Export Data Preparation Map Composer Image Interpreter Image Catalog Image Classification Spatial Modeler Radar Vector Virtual GIS A
遥感数字图像处理(ERDAS)
• 色彩变换(RGB-IHS)
– 将图像从红绿蓝彩色空间转换到以亮度、色度、 饱和度为定位参数的彩色空间,以便使图像颜 色与人眼看到的更接近。
• 指数计算
– 应用一定的数学方法,将遥感图像中不同波段 的灰度值进行各种组合运算,计算反映矿物及 植被等的常用比率和指数(植被指数,裸露指 数等)
傅立叶变换
植被指数: 水体指数: 建筑指数:
专题制图输出
根据工作需要和制图区域的地别特点,进行地图四面的整体设计,设计内 容也括图幅大小尺寸、图面布置方式、地图比例尺、图名及图例说明等; 需要淮备专题制图输出的数据层,也就是要在窗口中打开有关的图像或图 形文件; 启动地图编辑器,正式开始制作专题地图; 确走地图的内图框,同时确定输出地图所也含的实际区域范围,生成基本 的的制图输出图面内容: 在主要图面内容周围,放置图廓线、格网线、坐标注记,以及图名、图例、 比例尺、指北针等图廓外要素; 设首打印机,打印输出地图。
20 20
250/500/1000
产品
蓝
绿
红
近红外
短波红外
中波红外
热红外
全色
Landsat-7
1
1
1
1
2
1
1
Landsat-5
1
1
1
1
2
1
Spot-4
1
1
1
1
1
Spot-5
1
1
1
1
1
QuickBird
1
1
1
1
1
Modis
36
产品
多光谱数据
PAN数据
备注
全景范围
SPOT-2/4
erdas基本操作指导
Erdas基本操作指导书说明本指导书专为erdas新手编写,内容为erdas基本操作,包括打开图像的设置、波段叠加、图像融合、图像剪裁和简单的图像增强。
虽然简单,但是胜在语言通俗易懂,步骤详尽,傻瓜式操作。
作者鼓励读者任意转发,希望能够帮助到各位遥感初学者和erdas新手。
同时声明:本人保留著作权和署名权。
——L.R,2011.12.30 1.打开图像在viewer中点击File——open——raster layers,或者直接点击打开文件设置:点击Raster Option,将Clear Display 勾掉可在viewer窗口同时打开两幅图像。
勾选No Stretch 可显示图像原始色彩。
勾选Fit to Frame 可将图像调整到视窗范围内显示。
设置好后点击OK,即可加载图像。
2.波段叠加波段叠加是将几个波段合成一幅图像。
波段来源可以是单波段图像,也可以是多波段图像中的一个或几个波段。
来源图像空间分辨率可以不同(结果图像默认为分辨率最高的)、图像格式也可以不同(结果图像格式可自行设置)。
在erdas主界面上点击Interpreter——Utilities——Layer Stack。
弹出图像选择和叠加界面。
在Input File栏下选择输入图像,并在layer下拉框中选择要加入的波段(选择1、2、3…为单波段,选择ALL为所有波段),点击Add。
加载所有要合成的波段后,在Output File 栏下设置输出文件路径和名称。
在Data Type栏下方的Output下拉框中选择输出文件的数据率。
点击OK,等待运算完成即可。
3.图像融合图像融合是将分辨率较高、单波段的图像信息与分辨率较低、多波段的多光谱图像信息融合起来,使得结果图像即具有高分辨率的优点,又有丰富的多光谱信息。
操作方法:在erdas主界面上点击Interpreter——Spatial Enhancement——Resulotion Merge。
ERDAS IMAGINE遥感影像处理操作流程-2010全国土地利用变更调查监测与核查项目
北京望神州科技有限公司 北京海淀区上地三街9号嘉华大厦F座507室 邮编:100085 联系人:王鹤 Tel:(010)82781747/49‐107 手机:13439429782 邮箱:wanghe@
6
在状态栏可正要求 GCP 总体误差一般平坦区域要小于 1,山区小于 2。如 果误差较大,需要进行修改,删除点或增加新的控制点降低误差。 3、影像重采样 我们采集的点的数量、分布、精度都满足要求后,就可以进行影像的重采样 工作。在工具栏点击影像重采样 图标,打开影像重采样对话框。
RapidEye 原始数据 (5 波段)
每个波段的 NITF数据
DEM
控制资料
逐波段正射纠正 (选用NITF RPC模型)
波段合成
影像镶嵌
DOM成果
图1 使用RapidEye制作DOM的基本流程
北京望神州科技有限公司 北京海淀区上地三街9号嘉华大厦F座507室 邮编:100085 联系人:王鹤 Tel:(010)82781747/49‐107 手机:13439429782 邮箱:wanghe@
二 正射纠正 .......................................................................................................... 2 2.1 2.2 相关概念 ................................................................................................ 2 RapidEye 影像正射纠正 ...................................................................... 3 2.2.1 2.2.2 2.3 使用核心模块 .............................................................................. 3 使用 AutuSync ........................................................................... 9
ERDAS遥感图像基本处理傻瓜式步骤
ERDAS遥感图像基本处理傻瓜式步骤1、图像导入在erdas的Import/Export模块中,分别导入TM图像的第1、2、3、4、5、7波段,具体操作步骤为① 点击import模块,打开对话框②选择type类型为TIFF③ media为file;④ 然后选择输入、输出文件名路径和文件名⑤ 分别对123457波段进行导入;⑥在此之前可以选择session->preference,选择输入、输出主目录。
2、图像波段合成在erdas的interpreter模块中将单波段影像进行合成,生成多波段文件,具体操作步骤为:interpreter->utilities->layer stack,① 在出现的对话框中import框中依次选择需要合成的波段,每选择输入一个波段用Add添加一次;② output file选择导出文件路径及命名文件。
③ Data type 设为 Unsigned 8 bit;④Output option 设置为Union ,选中 ignore zero stats;⑤进行操作。
3、用shape文件进行图像切割3.1 Shape文件制作AOI文件:①在ERDAS中点击Import图标,出现Import/Export对话框②选中Imput,Type栏选择Shapefile,Media栏选择File,在Input File (*.shp)中确定要转换的shape文件,在Output File(*.arcinfo)中确定输出路径及名称,单击OK按钮,出现Import Shapefile对话框,单击Import Shapefile Now。
③注意此步骤中输出路径及输出名称均为英文字母④建立拓扑多边形⑤在Arcgis中打开ArcToolbox,Data ManagementTools—>Topology—>Build,双击Build,出现Build对话框,在Input 中填入*.arcinfo文件的路径,Feature选择Poly⑥单击OK按钮。
ERDAS入门基础教程-图像增强
遥感数字图像处理入门基础教程 ----图像增强
昆明信飞科技有限公司
2019/11/17
ERDAS IMAGINE
图像增强
遥感图像在获取的过程中由于受到大气层的散射、反射、折射或者天气的 影响,获得的图像难免会带有噪声或目视效果不好,如对比度不够、图像模糊; 有时总体效果较好,但是反需要的信息不够突出,如线状地物或地物的边缘部分; 或者,有些图像的波段较多,数据量较大,如TM影像,但各波段的信息量存在 一定的相关性,为进一步的处理造成困难。针对上述问题,需要对图像进行增强 处理。通过图像增强技术,改善图像质量,提高图像目视效果,突出所需要的信 息、压缩图像的数据量,为进一步的图像判读做好准备。
卷积增强
非定向边缘增强 聚焦分析 纹理分析 自适应滤波 统计滤波 分辨率融合 改进的HIS融合 HPF融合 小波融合 Ehlers融合 锐化处理
ERDAS IMAGINE
ERDAS IMAGINE
打开Convolution对话框后设置如下:
输出文件
输入文件
调整类型
卷积算子 类型
边缘处理 方法
输出数据 类型
输入文件
调整类型
输出数 据类型
输出文件
乘积倍数
输出数据 统计时忽 略O值
统计滤波前后效果对比图:
统计滤波前
统计滤波后
谢谢!
输入文件
调整类型
输出数 据类型
输出数据 统计时忽 略O值
输出文件 定义处理函数
纹理分析前后效果对比图:
纹理分析前
纹理分析后
ERDAS IMAGINE
六、自适应滤波
自适应滤波(Adaptive filter)是应用Wallis Adaptive Filter方法对图像的感兴趣区(AOI) 进行对比度拉伸,从而达到图像增强的目的。自适应滤波的关键问题是移动窗口(Moving Windows Size)的选择和乘积倍数大小(Multiplier)的定义,移动窗口的大小可以任意选择, 通常都定为奇数,而乘积倍数大小是为扩大图像反差或对比度。 设置如下:
ERDAS的操作手册
ERDAS的操作手册纠正,融合,镶嵌是遥感处理中比较常见的三种处理方法。
对于初学遥感的人来说,掌握这三种方法是十分必要的。
下面,我们通过一些实例,在ERDAS 中的操作,来分别介绍这三种处理方法。
1、纠正纠正又叫几何校正,就是将图像数据投影到平面上,使其符合地图投影系统的过程;而将地图坐标赋予图像数据的过程,称为地理参考(Geo-referencing)由于所有地图投影系统都遵从于一定的地图坐标系统,所以几何校正包含了地图参考。
(1)启动在ERDAS中启动几何校正有三种方法:A、菜单方式B、图标方式C、窗口栅格操作窗口启动这种方法比较常用,启动之前在窗口中打开需要纠正的图像,然后在栅格操作菜单中启动几何校正模块。
建议使用这种启动方法,更直观简便。
(2)设置几何校正模型常用模型:功能Affine 图像仿射变换(不做投影变换)Camera 航空影像正射校正Landsat Landsat卫星影像正射校正Polynomial 多项式变换(同时做投影变换)Rubber Sheeting 非线性、非均匀变换Spot Spot卫星图像正射校正其中,多项式变换(Polynomial)在卫星图像校正过程中应用较多,在调用多项式模型时,需要确定多项式的次方数(Order),通常整景图象选择3次方。
次方数与所需的最少控制点数是相关的,最少控制点数计算公式为((t+1)*(t+2))/2,公式中t为次方数,即1次方最少需要3个控制点,2次方需要6个控制点,3次方需要10个控制点,依此类推。
(3)几何校正采点模式A、Viewer to Viewer 已经拥有需要校正图像区域的数字地图、或经过校正的图像,就可以采用Viewer to Viewer的模式。
B、File to Viewer 事先已经通过GPS测量、或摄影测量、或其它途径获得了控制点坐标,并保存为ERDAS IMAGINE的控制点格式或ASCII数据文件,就可以采用File to Viewer模式,直接在数据文件中读取控制点坐标。
Erdas遥感图像处理(二)
ERDAS 遥感图像处理
图像空间增强处理——小波分辨率融合
ERDAS 遥感图像处理
4.1 图像增强
空间增强处理 辐射增强处理 光谱增强处理 高光谱工具 傅里叶变换
ERDAS 遥感图像处理
图像辐射增强处理
辐射增强(Radiometric Enhance-ment)技 术是对单个像元的灰度值进行变换达到图 像增强的目的。
ERDAS 遥感图像处理
空间增强(Spatial Enhancement)
• 利用像元自身及其周围像元的灰度值进行运算, 达到增强整个图像的目的。
• 卷积增强 • 非定向边缘增强 • 聚焦分析 • 纹理分析 • 自适应滤波 • 统计滤波 • 分辩率融合 • 锐化处理
ERDAS 遥感图像处理
图像空间增强处理——卷积处理
ERDAS遥感图像处理 (二)
刘海 测绘遥感信息工程国家重点实验室
2010-08-24
ERDAS 遥感图像处理
四、图像处理
图像增强 地形分析 实用工具 图像分类 空间建模
ERDAS 遥感图像处理
4.1 图像增强
目的
• 采用一系列技术改善图像的视觉效果, 提高图像的清晰度;
• 将图像转换成一种更适合于人或机器进 行解译和分析处理的形式。
ERDAS 遥感图像处理
图像空间增强处理——分辨率融合
ERDAS 遥感图像处理
图像空间增强处理——分辨率融合
(Resolution Merge )
对不同空间分辨率遥感图像的融合处理,使融合后的遥感图像既具 有较好空间分辨率,又具有多光谱特性,从而达到图像增强的目的.
ERDAS 遥感图像处理
图像增强不是以图像保真度为原则,而是 通过处理抑制一些无用的信息,增强人或 机器对某些信息的辨别能力。
erdas使用详细操作
erdas专题制图详细操作实验数据:1、裁剪144029(行列号)所需矢量文件2、卫星影像:LE71440292000188SGS01(数据获取来源为usgs网站:/下载所得。
)初步查看分析影像信息:首先在erdas中用tiff格式打开我们下载的影像:点击查看我们下载的影像的信息:可以看到我们下载的影像的分辨率,投影信息等。
图片格式转换:过程如下:在import\expot窗口中:这里的Type我们可以选择tiff 也可以选择geotiff,本人尝试过,选择这两种任意一种类型对输出影像几乎没有区别。
(其区别还有待查证)以此方法,将七个波段的.tiff影像图都转换成.img格式的影像。
由于band8 为全色影像,我们的实验中不需要用它,所以这里我不进行转换了。
只转换以下七个波段即可。
波段组合:将七个波段组合成.img文件设置完毕点OK。
于是我们就得到了.img格式的影像文件。
几何精校正:前面已经知道这张影像的地理信息了,为了验证其实图片投影大致的准确性,我们可以进行几何精校正进行验证。
几何精校正校正过程如下:这里我们使用多项式几何校正。
其中,多项式变换(Polynomial)在卫星图像校正过程中应用较多,在调用多项式模型时,需要确定多项式的次方数(Polynomial Order),通常整景图像选择3次方。
次方数与所需要的最少控制点数是相关的,最少控制点数计算公式为((t+1)*(t+2))/2,式中t为次方数,即1次方最少需要3个控制点,2次方最少需要6个控制点,3次方至少需要10个控制点,依次类推。
设置完毕点close。
会弹出GCP Tool Refrence Setup窗口。
由于我们要使用的是谷歌上的数据,因此我们要选择用键盘输入。
确认信息:信息正确,点OK。
在谷歌中设置:菜单栏中:工具->选项3D视图选项卡中,显示纬度\经度项中使用通用横轴墨卡托投影。
设置完后点应用和确定即可。
在谷歌中以及卫星影像中分别选取同名点,找到后,点击中的(或者点击窗口中的均可。
ERDAS遥感图像的几何校正
遥感图像的几何校正3、图像校正的具体过程第一步:显示图像文件(Display Image Files)首先,在ERDAS图标面板中点击Viewer图表两次,打开两个视窗(Viewer1/Viewer2),并将两个视窗平铺放置,操作如下:ERDAS图表面板菜单条:Session→Title Viewers然后,在Viewer1中打开需要校正的Lantsat图像:tmatlanta.img在Viewer2中打开作为地理参考的校正过的SPOT图像:panatlanta.img第二步:启动几何校正模块(Geometric Correction Tool)Viewer1菜单条:Raster→Geometric Correction→打开Set Geometric Model对话框,如图2图2 Set Geometric Model对话框→选择多项式几何校正模型:Polynomial→OK→同时自动弹出Geo Correction Tools对话框(见图3)和Polynomial Model Properties对话框(见图4)。
在Polynomial Model Properties对话框中,定义多项式模型参数以及投影参数:→定义多项式次方(Polynomial Order):2→定义投影参数:(PROJECTION):略→Apply→Close,自动弹出GCP Tool Referense Setup 对话框(见图5)图3 Geo Correction Tools对话框图4 Polynomial Properties对话框图5 GCP Tool Referense Setup 对话框第三步:启动控制点工具(Start GCP Tools)首先,在GCP Tool Referense Setup对话框(见图5)中选择采点模式:→选择视窗采集控制点的模式:Existing Viewer→OK→打开Viewer Selection Instructions指示器(见图6)图6 Viewer Selection Instructions→左键点击Viewer2中的panatlanta.img,使其作为地理参考图像(有准确的地理坐标)。
ERDAS遥感图像处理教程
《遥感数字图象处理实验指导书》实习须知实验室资源:1 硬件设备:局域网系统;高级PC 计算机 <每人一台,在规定时间使用);2 软件系统Windows 2000 或 Windows XPERDAS IMAGINE9.1或9.2ArcToolboxArcMapMicrosoft WordMicrosoft ExcelZip program3 数据资源数据源:老师提供实习目的与内容1 实习目的本实习为已具有RS的基本概念和理论基础的学生设计,目的是帮助学生在了解RS基本组成与数据结构模型的基础上,重点学习使用ERDAS IMAGINE9.1软件进行视窗操作、数据数据预处理、图像解译、图像分类和矢量功能;了解地图投影系统的使用;学习多种数据输入的方法,不同数据格式转换,数据库模式的定义等多种前后期处理工作;掌握遥感图像前后处理和解译、分类地理的技术流程和方法,数据库建设以及地理数据的编辑和管理;配合具体实例运用GIS空间分析工具。
通过系列实习过程,重点培养学习者掌握RS提取信息的基本过程和技巧,并可初步用来解决运用遥感提取信息的问题。
2 实习内容实习1:ERDAS IMAGINE 9.1系统简介与入门;包括:软件概述,视窗操作中的菜单工具条的介绍,数据的输入输出。
实习2:数据预处理;包括:图像分幅剪裁,图像几何校正,图像拼接处理,图像投影变换;实习3:图像解译;包括功能简介,辐射增强处理<去霾处理)、常用的光谱增强处理<假彩色合成与指数计算)、空间增强处理<分辨率融合)。
实习4:图像分类;包括;非监督分类和监督分类,专家分类器<在高级练习中学习)实习5:地理信息系统分析和矢量功能介绍,综合运用GIS工具解决实例提出的问题。
实习6:专题制图输出全过程学习。
3 本书有关约定实习所用原始数据存放在%sampledata目录下;实习1ERDAS IMAGINE 9.1系统简介与入门实习内容:1、了解有关ERDAS IMAGINE系统的基本概念和功能。
《ERDAS IMAGE遥感图像处理方法》操作
《ERDAS IMAGE遥感图像处理方法》操作一空间增强(Spatial Enhancement)1卷积增强处理(Convolution)功能:用一个系数矩阵将整个图像按照象元分块进行平均处理,用于改变图像的空间频率特征。
to效果:地物的轮廓和线条勾勒变清晰了。
2非定向边缘增强(Non-directional Edge)功能:应用两个非常通用的滤波器(Sobel 滤波器和Prewitt 滤波器),首先通过两个正交卷积算子(Horizontal算子和Vertical算子)分别对遥感图像进行边缘检测,然后将两个正交结果进行平均化处理。
to效果:效果明显而且强烈分别出邻区不同的部分。
3.聚焦分析(Focal Analysis)功能:使用类似卷积滤波的方法,选择一定的窗口呼函数,对输入图像文件的数值进行多种变换,应用窗口范围内的象元数值计算窗口中心象元的值,达到图像增强的目的。
to效果:深色地方变模糊,浅色地物图象得到增强,但也变得不清晰。
4.纹理分析(Texture Analysis)功能:通过二次变异等分析使图象的纹理结构更加清晰。
to效果:纹理边缘部分十分清晰。
5.自适应滤波(Adaptive Filter)功能:应用自适应滤波器对图像的感兴趣区域进行对比度拉伸处理。
to效果:颜色变浅了。
6.分辨率融合(Resolution Merge)功能:对不同空间分辨率遥感图像的融合处理,使处理后的遥感图像即具有较好的空间分辨率,又具有多光谱特征,达到图象增强的目的。
+ =效果:处理后图象既有高分辨率又有多光谱特征(彩色)。
7.锐化增强处理(Crisp Enhancement)功能:对图像进行卷积滤波处理,使整景图像的亮度得到增强而不使其专题内容发生变化。
效果:区别不大,亮度得到些许增强。
二.辐射增强(Radiometric Enhancement)1.查找表拉伸(LUT Stretch)功能:通过修改图像查找表使输出图像值发生变化。
ERDAS软件快速入门和处理影像之前的主要简单步骤
3、参数
地理空间数据云产品
(三)Landsat 4-5 TM与MSS 1、简介 发射时间:Landsat 4发射于1982年7月16日(2001年6月15日TM传感器失效,退役) Landsat 5发射于1984年3月1日(设计寿命3年,实际运行27年) 传感器类型:Landsat 4-5都配有主题成像仪(TM)和多光谱扫描仪(MSS)
2、命名;选择要素类型; 3、Click Edit > Select,选择需要的投影方式; Import > 找到需要处理的文件,直接引用投影方式。
3、在Arcgis中裁剪:掩膜提取(Extract by Mask) (1)打开ArcToolbox,Spatial Analyst Tools > Extraction > Extract by Mask
5、雾霾处理
6、解译
一、 Landsat卫星数据简介
(一)Landsat-8 1、简介
发射时间:2013 年2月11日
传感器类型:OLI陆地成像仪(Operational Land Imager)和TIRS热红外传感 器 (Thermal Infrared Sensபைடு நூலகம்r) 周期:16天,即每16天覆盖全球一次
确定输出图像保存位置:Output Cluster Layer Filename
确定分类模板文件位置:Output Signature Set Filename (2)Cluster Options分类参数设置 选择Initialize from Statistics按照图像的统计值产生自由聚类(另一种Use Signature Means按照选定的模板文件进行非监督分类); K Means:硬聚类算法(常用非监督分类方法) Isodata:迭代自组织数据分析技术
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
erdas遥感图像基本处理傻瓜式步骤2007年09月30日星期日 11:471、图像导入在erdas的Import/Export模块中,分别导入TM图像的第1、2、3、4、5、7波段,具体操作步骤为① 点击import模块,打开对话框②选择type类型为TIFF③ media为file;④ 然后选择输入、输出文件名路径和文件名⑤ 分别对123457波段进行导入;⑥在此之前可以选择session->preference,选择输入、输出主目录。
2、图像波段合成在erdas的interpreter模块中将单波段影像进行合成,生成多波段文件,具体操作步骤为:interpreter->utilities->layer stack,① 在出现的对话框中import框中依次选择需要合成的波段,每选择输入一个波段用Add添加一次;② output file选择导出文件路径及命名文件。
③ Data type 设为 Unsigned 8 bit;④Output option 设置为Union ,选中 ignore zero stats;⑤进行操作。
3、用shape文件进行图像切割3.1 Shape文件制作AOI文件:①在ERDAS中点击Import图标,出现Import/Export对话框②选中Imput,Type栏选择Shapefile,Media栏选择File,在Input File (*.shp)中确定要转换的shape文件,在Output File(*.arcinfo)中确定输出路径及名称,单击OK按钮,出现Import Shapefile对话框,单击Import Shapefile Now。
③注意此步骤中输出路径及输出名称均为英文字母④建立拓扑多边形⑤在Arcgis中打开ArcToolbox,Data ManagementTools—>Topology—>Build,双击Build,出现Build对话框,在Input 中填入*.arcinfo文件的路径,Feature选择Poly⑥单击OK按钮。
⑦在ERDAS中打开一个viewer窗口,打开arc coverage文件,新建一个aoi层(New—>AOI Layer)⑧ View—>Arrange Layers Viewer打开Arrange Layers Viewer对话框,在Vector图层上单击右键,选择Show Properties,打开Properties对话框,选中Polygon,点击Apply按钮。
⑨在View窗口中打开AOI工具栏,先选中内部区域,再点击,产生aoi,选中该aoi,在aoi下拉菜单选择“Copy Selection to AOI”,单击File—>Save—AOI Layer as,保存为aoi文件。
OK!3.2 用AOI文件进行对遥感图像切割在ERDAS图标面板工具条中单击Data Prep图标,Subset,打开Subset对话框。
在Subset对话框中需要设置下列参数:⑩输入文件名称(Input File)⑪输出文件名称(Output File)⑫单击AOI按钮确定裁剪范围⑬打开Choose AOI对话框⑭在Choose AOI对话框中确定AOI的来源为File(或Viewer)⑮如果选择了文件(File),刚进一步确定AOI文件;否则,直接进入下一步⑯输出数据类型(Output Data Type)为Unsigned 8 Bit,输出文件类型(Output Layer Type)为Themetic⑰输出象元波段(Select Layers)为1 :6(表示选择1-6六个波段)⑱输出统计忽略零值,选中Ignore Zero In Output Stats复选框⑲单击OK按钮4、图像预览在开始进行分类之前,需要先仔细查看合成后的图像,确定卫片所覆盖的地理区域及海拔范围,确定主要的地标性元素。
由于卫片原始数据的差异,同样的地表覆盖物斑块在不同时期的卫片中或不同地区的卫片中显示的颜色都可能不同,所以在预览中,还要熟悉整个范围内的地表覆盖类型和不同地物之间的变化。
①打开显示窗口,加载裁切后的6通道的图像(4,3,2)或者(4、5、3)、(7,4,2);②把图像缩小至适合窗口,浏览图像,注意河流、城镇、植被、水体、土壤的分布;③选择特定区域放大,查看各种不同的地表覆盖物类型的分布及色调变化;根据经验,在4,3,2(RGB)的波段组合下,各种地表覆盖物类型的特点如下:a.森林——森林显示出棕色、红色、褐色等一系列多变的色调。
在高海拔地区,成熟针叶林为很浓的棕色或暗红色;在中低海拔地区,森林的颜色多变,从棕色到红色到暗绿色都有,部分落叶林在冬季呈现出锈黄色;b.灌丛和草甸——相对于临近的森林斑块,灌丛和草甸呈现出明亮许多的红色到浅红色。
在高海拔地区,大片的草地在夏季可能为浅红或锈红色,而冬季则呈现青绿色;c.湖泊、河流——湖泊通常为边界清晰的黑色斑块,河流则显示为黑色或深蓝色。
在冬季,水面结冰或覆盖有雪则显示出不同深浅的紫红色;d.城镇——很明显的比较亮的灰色或青灰色斑块,通常可见有规则的灰色线条(公路)穿过;e.农田——颜色多变的绿色、灰色、淡紫色、浅红色斑块,通常沿河谷两侧不规则分布,在平原区则大片分布。
河道边的水田往往显示出富含水分的青灰色。
5、图像分类5.1 进行非监督分类步骤:第1步:启动非监督分类在ERDAS图标面板工具条中单击Classifier图标,打开Classification 对话框,单击Unsupervised Classification 按钮,打开Unsupervised Classification对话框第2步:进行非监督分类在Unsupervised Classification对话框中进行下列设置:① 确定输入文件(Input Raster File)(要进行分类的文件);②确定输出文件(Output File)(产生的分类文件),文件名定为ppprrr_YYYYMMDD_123457_unsupervised_15.img;③选择生成分类模板文件Output Signature Set,确定模板文件名称,命名同上;④确定聚类参数(Clustering Options),需要确定初始聚类方法与分类数:⑤默认选择Initialize from Statistics(按照图像的统计值产生自由聚类);⑥确定初始分类(Number of classes)为15(分为15类);⑦单击Initializing Options按钮,打开File Statistics Options对话框,设置一些统计参数,一般采用默认值;⑧单击Color Scheme Options按钮,打开Output Color Scheme Options 对话框,设置分类图像彩色属性,此处单击Approximate True Color,采用RGB 对应453波段合成。
⑨ 其他参数采用默认值。
⑩单击OK按钮(关闭Unsupervised Classification对话框,执行非监督分类)。
5.2 定义分类模板(1)步骤:Main->Image Classification->Classification->Signature Editor,打开分类模板编辑器。
在Viewer窗口下的Raster下打开Tools图标,选择多边形AOI绘制。
(2)定义模板原则①必须在分类之前就知道研究区域的森林类型、覆盖范围以及图像的叠和现象,以保证输出分类的连续性。
②当创建训练区时,对于每一个类别都有一些子类,每个子类选择的AOI区域应该不少于5个,并且每个AOI区域内象素的颜色类型一致,跳跃不能很大,即不出现杂色。
5.3 执行监督分类①依次选择:Main->Image Classification->Classification->Supervised Classification,打开监督分类对话框。
②输入原始文件③定义输出文件④确定分类模板文件⑤选择输出分类距离文件为Distance File⑥定义分类距离文件⑦选择非参数规则(Non-Parametric Rule)为Feature Space⑧选择叠加规则(Overlay Rule)为Parametric Rule⑨选择未分类规则(Unclassified Rule)为Parametric Rule⑩选择参数规则为Maximum Likelihood(即最大似然法)⑪取消选中Classify zeros复选框⑫ OK执行监督分类。
5.4 后期检查修正打开两个viewer窗口,进行链接(选择link工具)。
并可以选择aoi 的显示功能检查子类选择的正确性。
反复验证、修改模板。
5.5 重新分类重复以上步骤,重新分类,达到最佳分类结果6、图像拼接7、分类重编码将分类结果图像进行分类重编码,减少分类数量。
判断每个分类的专题属性,对相近或类似的分类通过图像重编码进行合并,并定义分类名称和颜色。
① Main->Image Interpreter->GIS Analysis->Recode②确定输入、输出文件;③ 设置新的分类编码(Setup Recode),打开Thematic Recode表格,根据需要改变“New Value”字段取值(直接输入);④单击OK;⑤单击OK8、滤波①Image Interpreter | GIS Analysis... | Eliminate...;② 输入文件为“.._ clp4.img”,输出文件为“..._ elim25.img”;③ “Minimum”选择“25”;(因为象素分辨率为28.5M,25个象素接近于2公顷;这样最小图斑为2公顷)④ “Output”选择“8bit”9.erdas配准步骤配准步骤:要求1。
在viewer模块里打开你的正确得影像。
2。
点击dataprep模块,打开他下面得image geometric correction子模块,3选择from image file,选择你要配得图像,打开4,在出现得set geometric model里选择第2个polyonial5。
在下面得对话框里有个polynomial order 一般选2就可以6,然后应用,在close那个对话框7 在下面得对话筐里,选第一个8出现一个小的对话框,不用管他,点击你第一步打开得准确得图像界面,就可以配准了9在你得准确图像里面点击一个点,然后转到你的需要配的图象里面找与他相似的点,点击,选够6个点后,在gcp tool那个界面里面将出现误差参数,可以看你得点是不是准确。