[转载]ERDAS遥感图像基本处理+ENVI求相关系数矩阵
(2016-04-29 11:38:40)
标签:
转载 |
分类: 遥感 |
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添加一次;
②
③ Data type 设为 Unsigned 8 bit;
④Output option 设置为Union ,选中 ignore
zero
⑤进行操作。
3、 用shape文件进行图像切割
3.1 Shape文件制作AOI文件:
①
②
③
④
⑤
⑥
⑦
⑧
⑨
3.2 用AOI文件进行对遥感图像切割
在ERDAS图标面板工具条中单击Data
Prep图标,Subset,打开Subset对话框。在Subset对话框中需要设置下列参数:
⑩
⑪
⑫
⑬
⑭
⑮
⑯
⑰
⑱
⑲
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)(要进行分类的文件);
②
③
④
⑤
⑥
⑦
⑧
⑨ 其他参数采用默认值。
⑩
5.2 定义分类模板
(1)步骤:
Main->Image
Classification->Classification->Signature
Editor,打开分类模板编辑器。在Viewer窗口下的Raster下打开Tools图标,选择多边形AOI绘制。
(2)定义模板原则
①必须在分类之前就知道研究区域的森林类型、覆盖范围以及图像的叠和现象,以保证输出分类的连续性。
②当创建训练区时,对于每一个类别都有一些子类,每个子类选择的AOI区域应该不少于5个,并且每个AOI区域内象素的颜色类型一致,跳跃不能很大,即不出现杂色。
5.3 执行监督分类
①
②
③
④
⑤
⑥
⑦
⑧
⑨
⑩
⑪
⑫
5.4 后期检查修正
5.5 重新分类
重复以上步骤,重新分类,达到最佳分类结果
6、图像拼接
7、分类重编码
将分类结果图像进行分类重编码,减少分类数量。判断每个分类的专题属性,对相近或类似的分类通过图像重编码进行合并,并定义分类名称和颜色。
①
②
③ 设置新的分类编码(Setup Recode),打开Thematic Recode表格,根据需要改变“New
Value”字段取值(直接输入);
④
⑤
8、 滤波
①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个polyonial
5。在下面得对话框里有个polynomial order 一般选2就可以
6, 然后应用,在close那个对话框
7 在下面得对话筐里,选第一个
8出现一个小的对话框,不用管他,点击你第一步打开得准确得图像界面,就可以配准了
9在你得准确图像里面点击一个点,然后转到你的需要配的图象里面找与他相似的点,点击,选够6个点后,在gcp
tool那个界面里面将出现误差参数,可以看你得点是不是准确。这6个点是控制点,所以你选择时应该分散开,尽量在全图范围里找。
10。然后选择其他点越多越好
11,如果你不想那样麻烦的找,你可以在任何一幅图像上点击鼠标右键,然后选择geo。link/unlink选项,那样子的话,你只需要在一个图象里点点,机器自动给你在另一个里面找,
不过头六个点不会出现误差参数,从第七个点开始会出现误差参数,如果误差太大,你必须在修改。
12,在你配完后,在geo correction tools里面点击第3个图标就ok
13,如果你配准到中间要休息,需要保存gcp tools对话框,
下次在匹配直接打开就行。
14。配准完后得总误差也就是那个rms误差必须在0。5个像元以里。
工作完毕,你的图就ok
=================================================================================================另外:如何利用Erdas的Make Modeler计算TM的最佳波段组合(其实用ENVI是很好作)
1。先利用Erdas的spatial modeler->make
model计算出个波段之间的相关系数矩阵。这一步的关键是要在矩阵输出时候选择输出到文件(即保存为扩展名为*.mtx的文件),根据你的遥感数据的实际波段数定制矩阵的维数,比如6*6矩阵将记录六个波段间的相关系数;
2。收集各波段的相关信息。对于OIF指数来说主要是每个波段的标准差(sta
dev),如果其他的算法还需要其他的信息,可以通过layer info分别记录;
3。根据OIF指数的计算公式,计算OIF值。这个过程可能有些麻烦,因为如果波段比较多的话,人工计算相当繁琐。
如何在ENVI中求出相关系数矩阵
很简单地,在统计功能下面就能实现,注意要在对话框中把你需要的参数都要打勾选上
有最小值,最大值,直方图,相关系数,协方差等等
ENVI中,选择 Basic Tools > Statistics
> Compute Statistics
然后选择 “Calculate Covariance Statistics” 复选框, 通过选中 “Covariance Image”
和/或 “Text Report”
复选框,选择你想把协方差、相关和特征向量矩阵输出到图像,或一个文本记录,还是输出到两者。
erdas中也可以实现:
modeler——model maker
Create a raster object—— Create a function definition——Create a
matrix object (矩阵选择output形式,意在把相关性存成一个矩阵)
双击raster,定义为你要处理的影像;
双击function,定义为 CORRELATION ( <raster>
) ,用你选择的影像名字替代括号中的<raster>;
双击matrix,定义输出文件。
ENVI里的密度分割,二值化,灰度图像赋彩色分类等操作。打开图像,在Image的Tool菜单下拉的Color map二级菜单里有密度分割density slice等,在此进行操作即可。
在erdas里植被盖度图怎么划分等级?打开图像时,Select Layer To Add窗口中点击Raster
Options卡片,Display as 后面下拉框选择Pseudo Color,打开图像后,点击菜单Raster->Attributes,弹出窗口,就能看到颜色了,单击你所看到的颜色框,就可以随意设置所需要的颜色了。另外,可以在左侧
Row一栏里同时选择多行。
=============================================================================================