清早起床撸啊撸
9月 日
地理与资源科学学院 级 实验名称:专题地图制作 姓 名: 成 绩: 学 号: 指导教师: 1. 实验目的和要求
利用一幅遥感影像, 运用所学知识对其进行处理, 再提取专题信息, 制成一幅专题地图。 2.主要仪器与试剂
ERDAS 软件、 Modeler output影像 3.实验步骤:
一、图像投影变换
过程:(1)启动投影变换
在图像解译模块(Image Interpreter)中启动 :ERDAS图标 面板菜单条 :Main->Image Interpreter->Utilities->Reproject Images->Reproject Images对话框(图 1)
(2)投影变换操作 ① 确定输入图像文件(Input File) : Modeler output.img; ② 定义输出图像文件 (Output File) : Modeler output sub.img ; ③ 定义投影类型(Categories ) : UTM Clark 1866 North; ④ 定义投影参数(Projection ) : UTM Zone 50 (Range
114E-120E ) ; ⑤ 定义输出图像单位(Units ) : Meters; ⑥ 确定输出统计缺省零值 :Ignore Zero in Stats; ⑦ 定义输出像元大小(Output Cell Sizes) : X: 0.5 Y: 0.5; ⑧ 选择重采样方法(Resample Method) : Nearest Neighbor; ⑨ 定义转换方法 : Polynomial Approximation(应用多项近似
拟合实现变换)
多项式最大次方(Maximum Poly Order) : 3;定义像元误差(Tolerance Pixels) : 1; (3) OK 分析:图像投影变换的目的在于将图像文件从一种地图投影类型转换到另一种投影类型, 这 种转换可以对单幅图像进行, 也可对多幅图像进行。 与图像几何校正过程中的投影变换相比, 这种直接的投影变换可以避免多项式近似值的拟合,对于大范围的地理参考是非常有意义 的。 在变换过程中, 如果在设置的最大次方内没有达到像元误差要求, 则按照下列设置执行:If Tolerance Exceeded: Continue Approximate; Rigorous Transformation
二、光谱增强
过程:彩色变换
突出显示水体信息 用 4、 3、 2作假彩色合成,具体操作如下:在 Viewer 窗口中打开 Modeler output.img图像后,在 Viewer 窗口中单击 Raster 工具条,单击 Band Combinations,
图 1Reproject Images 对
话框(设置参数之后)
出现 Set Layer Combinations for 窗口(如图 2)
在对话框中分别选择 Red 为 4, Green 为 3, Blue 为 2,单击 Apply 和 ok 。 分析:光谱增强处理是基于多波段数据对每个像元的灰度值进行变换, 达到图像增强的目的。 根据加色法的原理,选择遥感影像的某三个波段,分别赋予红、绿、蓝三种原色,就可以合 成彩色影像。 多波段影像合成时, 方案的选择十分重要, 它决定了彩色图像能否显示较丰富 的地物信息或突出某一方面的信息。
通过实验发现:合成之后的图像可以更好的反映水体污染的情况, 清洁水体的影像呈均 匀的深蓝 -蓝黑色,轻污染区的水体影像蓝中夹杂着不均匀的呈条状、块状分布的淡蓝、淡 青色调,重污染区的水体影像呈现淡蓝色 -青色。
三、图像分类
过程:监督分类
1、 定义分类模板
(1)显示需要进行分类的图像
在视窗中显示 Modeler output.img。选择 Fit to Frame,其它采用缺省设置。 (2)打开模板编辑器并调整显示字段
ERDAS 图标面板工具条 :点击 Classifier 图标 ①
Classification 菜单;② Signature Editor菜单项; ③ Signature Editor对话框 (图 3)
图 2 Set Layer Combinations for 窗口
图 3 Signature Editor对话框
对字段进行如下调整 :Signature Editor对话框菜单条 : 单击 View->Columns, 打开 View Signature Columns对话框(图 4)
①点击最上一个字段的 Column 字段向下拖拉直到最后 一个段。此时,所有字段都被选择上,并用黄色(缺省 色)标识出来。
②按住 Shift 键的同时分别点击 Red, Green, Blue三个字 段 Red, Green, Blue三个字段将从选择集中被清除 . ③点击 Apply 按钮 ④点击 Close 按钮 (3)保存分类模板
Signature Editor对话框菜单条 :单击 File->Save ①打开 Save Signature Fiel As对话框
②确定是保存所有的模板还是只保存被选中的模板 ③确定文件的目录和名字(.sig 文件) ④点击 OK 按钮
2、 评价分类模板 过程:产生报警掩膜
(1)产生分类预警掩摸(首先打开获取分类摸板后的 Signature Editor对话框,图 5)
图 5获取分类摸板后的 Signature Editor对话框
在 Signature Editor对话框 :①选择某一类或者某几类模板; ②单击 View->Image Alarm命令;
③打开 Signature Alarm对话框(图 6) ④选中 Indicate Overlap;⑤点击 Edit Parallelepiped Limits 按钮;⑥ Limits 对话框(图 7);⑦点击 Set 按钮,打开 Set Parallelepiped Limits对话 框 (图 7) ; ⑧设置计算方法 : Minimum/ Maximum; ⑨选择使用的模板 (Signature ) :Current
;
图 4 View Signature Columns对话框
图 6 Signature Alarm对话框
图 7 Limits对话框
⑩ OK ;○
11返回 Limits 对话框;○ 12Close ;○ 13返回 Signature Alarm对话框;○ 14 OK ;○ 15 Close (2)利用 Flicker 功能查看报警掩膜。 (3)删除分类报警掩膜
视窗菜单条 :View->Arrange Layers菜单,设置:①打开 Arrange Layers对话框;②右键点击 Alarm Mask图层;③弹出 Layer Options菜单:选择 Delete Layer;④ Alarm Mask图层被删 除;⑤ Apply ;⑥提示 Verify Save on Close->No;⑦ Close 3、执行监督分类
ERDAS 图标面板菜单条 :Main->Image Classification->Classification菜单 (1) Supervised Classification菜单项
(2) Supervised Classification对话框,需要确定下列参数 : ①确定输入原始文件(Input Raster File) : Modeler output.img ②定义输出分类文件(Classified File) : Modeler outputa.img ③确定分类模板文件(Input Signature File) : supervise.sig ④选择输出分类距离文件 :Distance File
⑤定义分类距离文件(Filename ) : Modeler output.img ⑥选择非参数规则(Non-Parametric Rule) :Feature Space ⑦选择叠加规则(Overlay Rule) : Parametric Rule
⑧选择未分类规则(Unclassified Rule) : Parametric Rule ⑨选择参数规则(Parametric Rule) : Maximum Likelihood ⑩不选择 Ciassify zeros ○
11OK (1)在 Attribue Options对话框上作出选择 (2) OK
(3)返回 Supervised Classification对话框 4、评价分类结果
分类叠加就是将专题分类图像与分类原始图像同时 在一个视窗中打开,将分类专题层置于上层,通过改变分
类专题的透明度及颜色等属性, 查看分类专题与原始图像
之间的关系。 对于非监督分类对果。 通过分类叠加方法来确定类别的专题特性、 并评价分类 结果。对监瞥分类结果。该方法只是查看分类结果的准确性。
分析:图像分类是基于图像像元的数据文件值, 将像元归并成有限几种类型、 等级或数据集 的过程。常规图像分类主要有两种方法:非监督分类与监督分类。
本实验采用监督分类, 监督分类比非监督分类更多地要用户来控制, 常用于对研究区域 比较了解的情况。监督分类一般有以下几个步骤 :定义分类模板、评价分类模板、进行监督 分类、评价分类结果。 ERDAS IMAGINE的监督分类是基于分类模板来进行的,而分类模 板的生成、管理、评价、和编辑等功能是由分类模板编辑器来负责的。
如果对同一个专题类型 (如水体) 采集了多个 AOI 并分别生成了模板, 可以将这些模板 合并, 以使该分类模板具多区域的综合特性。 此分类模板的检验精度比较好, 可以进行监督 分类操作。无论监督分类还是非监督分类,都是按照图像光谱特征进行聚类分析的,因此, 都带有一定的盲目性。 所以, 对获得的分类结果需要再进行一些处理工作, 才能得到最终相 对理想的分类结果。
四、制专题地图
过程:1专题地图编辑过程
图 8 Attribue Options对话框
1.1准备专题制图数据
视窗菜单条 :File->Open-Raster Layer ① File Name: Modeler output.img; ② Raster Options: Fit to Frame; ③ OK
1.2产生专题制图文件
ERDAS 图标面板菜单条 :Main->Map Composer->New Map Composition New Map Composition对话框(图 9),需要定义下列参数 :
①专题制图文件名 (New Name) : composer.map; ② (Map Width) :28; ③ (Map Height) : 20;④(Display Scale) : 1;⑤(Unit ) :centimeters;⑥(Background ) : white; OK 打开 Map Composer视窗和 Annotation 工具面板(图 10)
1.3确定专题制图范围 地图图框的绘制过程如下 :
Annotation 工具面板 ->点击 Create Map Frame图标
(1)在地图编辑视窗的图形窗口中,按住鼠标左键拖动绘制一个矩形框(Map Frame)。 (2)完成图框绘制、释放鼠标左键后,打开 Map Frame Data Source对话框(图 11) (3)点击 Viewer
(4)打开 Create Frame Instructions指示器(图 12)
在显示图像的视窗中任意位置点击左键,表示对该图像进行专题制图 ①打开 Map Frame对话框(图 13)
② Map Frame对话框中需要定义下列参数 :
③ rame Width: 24 Frame Height: 16 (Map Area Width&Height相应变化) ④ hange Scale and Frame Area, Scale 1:25000 ⑤(Map Angle) : 0
⑥(Upper Left Map Coordinates) : X: 1700000.00 / Y: 290000.00 ⑦(Upper Left Frame Coordinates) : X: 2.00 /Y:18.00) ⑧ ok
图 9 New Map Composition对话框
图 10 Map Composer视 窗 和 Annotation 工 具 面 板
图 11Map Frame Data Source对话框
图 12 Create Frame Instructions指示器框
⑨(View->Scale-> Map to Window)(图 14)
2放置图面整饰要素
2.1绘制格网线与坐标注记
Annotation 工具面板 :点击 Create Grid / Ticks图标
在位于地图编辑视窗图形窗口中的图框内点击左键 在 Set Grid / Tick Info对话框中,需要设置下列参数 : ①(Name ) : composer_grid
②(Description ) : grid, tick and neatline of composer ③选择放置地理坐标注记要素 :Geographic Wicks 复选框 ④选择放置地图图廓线要素 :Neat Line ⑤(Margin ) : 0.200 Centimeters ⑥(Unit ) : Feet ⑦(Horizontal Axis) ⑧(Length Outside) :0 ⑨(Length Inside) : 0 ⑩(Length Intersect) : 0
○
11(Starting at) : 1700000 Feet ○
12(Spacing ) : 1000 Feet ○
13选择使用完整格网线 :Use Full Grid ○
14
Apply 图 13 Create Frame Instructions对 话 框 图 14专题制图图面(充满窗口 )
○ 15 Close (如图 15)
图
15
2.2绘制地图比例尺
Annotation 工具面板 :点击 Create Scale Bar图标
在 Map Composer图形窗口中合适的位置按住左键拖动鼠标,绘制比例尺放置框 打开 Scale Bar Instructions指示器(图 16)
打开 Scale Bar Properties对话框(图 17) ①(Name ) : Scale Bar
②(Description ) : Scale Bar for Composer ③(Title ) :比例尺
④(Alignment ) : Zero ⑤(Units ) : Meters
⑥(Maximum Length) : 3 centimeters ⑦ Apply ⑧ Close (图 18) 3绘制地图图例
在 Annotation 工具面板 :点击 Create Legend图标 ①在 Map Composer视窗的图形窗口中合适的位置点
图 16 Scale Bar Instructions指示器
图 17 Scale Bar Properties对话框 图 18
图 19 Legend Properties对 话 框
击左键,定义放置图例左上角位置 ②打开 Legend Instruction指示器
, 在 Map Composer视窗的图形窗口制图框中点击左键,
指定绘制图例的依据;打开 Legend Properties对话框(图 19),需要分别设置参数
(2)标题参数 :
①(Title Content) :图例;② Underline Title 复选框;③(Title / Underline Gap) : 2 points; ④(Title / Legend Gap) : 12 points;⑤(Title Alignment) : Centered
⑥(Legend Unit) : Point;⑦(Columns Properties);⑧选中 Use Multiple Column 复选框;
⑨ (Entries per Column) : 15 ⑩ (Gap Between Columns) : 20 points; ○ 11(Gap Between Entries) :
7.5 points;○
12(Heading / First Entries Gap) : 12 points;○ 13(Text Gap) : 5 points;○ 14选中 Vertically Stack Descriptor Text 复选框 (3)色标参数(Color Patches) :
①选中 Place Patch Left of Text复选框; ②选中 Outline Color/Fill Patch 复选框; ③选中 Outline Symbol/LinelText Patch复选框; ④(Patch Width) : 30 points; ⑤(Patch Height) : 10 points; ⑥(Patch /Text Gap ) : 10 points;
⑦(Patch/ Text Alignment) : Centered;
(4)(Legend Units) : Points 图 20 (5) Apply
(6) Close (如图 20) 4绘制指北针
(1)确定指北针符号类型
Map Composer视窗菜单条 :Annotation->Styles ①打开 Styles for Composer对话框 ②选择 Symbol Styles ③选择 Other
④打开 Symbol Chooser对话框(图 21) 在 Symbol Chooser对话框中, 确定指北针类 型 :Standard->North Arrows->North Arrow 2 ①(Use Color)并选择指北针颜色; ②(Size ) : 30
③(Units ) : Paper pts; ③ Apply
④ OK ;⑤ Close
(2)放置指北针符号(Create Symbol) ① Annotation 工具面板 :点击 Create Symbol,图标
;②在 Map Composer视窗的图形窗口
中点击左键, 放置指北针; ③双击刚才放置的指北针符号; ④打开 Symbol Properties对话框 (图 22);
图 21 Symbol Chooser对话框
图 22 Symbol Properties对话框
在 Symbol Properties对话框中,确定指北针要素特性 :
①指北针要素名称 (Name ) :North Arrow; ②指北针要素描述 (Description ) :North Arrow for Composer ;③指北针符号中心位置坐标 :Center X: 25.5 , Center Y: 1.5;④选择中心位置坐标 类型与单位 :Type: Map Units : Centimeters;⑤指北针符号旋转角度及单位 :Angle: 0.00 Units: Degree ;⑥指北针符号大小尺寸(Size ) : 30;⑦选择符号尺寸类型与单位 :Type:Paper Units:Points;⑧ Apply (应用指北针符号特性定义参数);⑨ Close (图 23)
图 23 composer 对话框
5放置地图图名
(1)确定图名字体
Map Composer视窗菜单条 :Annotation->Styles
①打开 Styles for Composer对话框;②选择 Text Styles(;③选择 Other ;④打开 Text Style Chooser 对话框(图 24)
Standard 栏目,设置下列参数 :①选择图名字体 :Black Galaxy Bold;②(Size ) : 10;③ (Units ) :Paper pts;④ Apply ;⑤
OK
图 24 Text Style Chooser对话框
(2)在 Symbol Chooser对话框中,点击 Custom 标签,进入 Custom 栏目(图 25)
① (Size ) : 10 paper pts; ② Goudy-Old-Style ; ③ (Italic Andle) : 15.0; ④ (Underline Offset / Width) : 15 / 5;⑤(Shadow Offset X / Y) : 2 / 2;⑥(Fill Style);⑦ Apply ;⑧
OK
图 25 Text Style Chooser对话框
(3)放置地图图名
Annotation 工具面板 :点击 Create Text图标
①在 Map Composer视窗的图形窗口中点击左键,确定放置图名位置 ②打开 Annotation Text对话框(图 26
)
图 26 Annotation Text对话框
③在 Annotation Text对话框中输入图名字符串 ,“Land Use and Land Cover Image Map of China
④ OK (如图 25)
图 25
6保存专题制图文件
通过上述过程所生成的专题制图文件可以保存起来,以便修改和应用,具体过程如下 : Map Composer工具条 :点击 Save Composition图标 ->保存制图文件(*.Map)
Map Composer菜单条 :File->Save->Map Composition->保存制图文件(*.Map)
分析:ERADS IMAGINE的专题地图编辑器是一种所见即所得编辑器,用于产生地图质量 的图像和演示图。 地图图框确定了专题制图图面的主要内容与区域, 在此基础上, 放置图廓 线、 格网线、 坐标注记、 图名、 图例、 指北针、 比例尺等各种辅助要素, 以便图面美观实用。 格网线的样式取决于 style 对话框中的缺省参数设置,可预先设置,亦可随时修改。通 过 Create Grid / Tick功能放置的格网线、 坐标注记、 图廓线是一个自然的图形组或组合对象, 可以调整其类型,并对其某一要素进行编辑。 比例尺、图例也是一个组合要素,如果进行局 部修改,需首先解散要素组合,然后再编辑单个要素。
11
实验一erdas
实验一 计算机遥感图像处理
一、实验目的要求
1.了解计算机遥感图像处理系统和基本内容;
2.了解计算机遥感图像处理方法及过程;
3. 熟悉遥感图像波段合成和边缘提取的方法, 并掌握具体操作过程; 二、实验用片及仪器
软件准备:ERDAS 遥感图像处理系统
数据准备:利用 example 中的 lanier.img 数据及其相关的附属数据。 三、实验内容
了解计算机遥感图像处理系统
LANDSAT 的 TM 图像的真彩色合成是 3 2 1
LANDSAT 的 TM 图像的标准假彩色把红设成 4绿设成 3蓝设成 2, 就是标准假彩色, 其余的组合就叫假彩色, 在标准加彩色图像中红色 的是植物,黑色的一般是水体。
1图像波段合成
在 erdas 的 interpreter 模块中将单波段影像进行合成,生成多波段文 件,具体操作步骤为:
interpreter ->utilities->layer stack,
① 在出现的对话框中 import 框中依次选择需要合成的波段,每选择 输入一个波段用 Add 添加一次;
② output file选择导出文件路径及命名文件。
③ Data type 设为 Unsigned 8 bit;
④ Output option 设置为 Union ,选中 ignore zero stats ;
⑤进行操作
2 卷积增强处理
单击 Interpreter/Spatial Enhancement/Convolution; 打开 Convolution 对 话框 :
注意:确定输入文件(Input File)为 lanier.img;
选择卷积算子(Kernel Selection)
卷积算子文件(Kernel Library)为 default.klb;
卷积 算子类型 5×5 Edge Detect (7 ×7 low pass; 7×7 Vertical Edge detection; 7×7 Horizontal Edge detection)
文件坐标类型(Coordinate Type)为 Map;
输出数据类型(Output Data Type)为 Unsigned 8 bit. 其他默认。
单击 OK 。
3去条带处理
单击 Interpreter/Radiometric Enhancement/Destripe TM Data, 打开 Destripe TM 对话框
注意:确定输入文件(Input File)为 tm_striped.img;
输出数据类型(Output Data Type)为 Unsigned 8 bit.
输出数据统计时忽略零值,选中 Ignore Zero in Stats 复选 框 ;
边缘处理方法(Handle Edges by)为 Reflection. 其他默认
单击 OK 。
四 作业要求 :
输出标准假彩色图像,边缘提取的成果图(一个) 。
erdas实验报告
ERDAS实验报告 图像融合实验
运用HIS变换法，对Atlanta多光谱影像与全色光影像进行融合，以得到高分辨率的彩色图像。
孟晓于
2014/5/10
2014-5-10
ERDAS实验报告
图像融合实验
, 数据来源
采用Erdas中examples文件内的2000年Atlanta多光谱TM数据和高清全色Pan数据。两图为同一地区不同坐标影像，故使用前需预处理从而得到实验区域。
, 目的
多光谱TM数据分辨率较低但包含多波段色彩，而全色Pan数据只包含一层高清影像，为了得到研究区域的高清彩色影像，我们将TM和Pan数据在Erdas2014中进行融合以达到实验目的。
, 方法
在遥感领域运用较多的融合方法有主成分变换法、比值变换法、小波变换法和HIS变换法。本实验则运用HIS变换法。IHS属于色度空间变换，从多光谱彩色合成影像上分离出代表信息的明度(I)和代表光谱信息的色调(H)、饱和度(S)等3个分量，并采用相同区域的高分辨率全色波段数据代替明度(I)进行空间信息融合。
, 步骤
1. 几何校正
因原始图像空间坐标不同 ，需选取控制点进行几何校正。本实验校正方法为多项式法，选取6个控制点进行校正，其校正叠加截图如下:
1
2014-5-10
2. 叠加剪切
由校正结果可知两图像只有部分区域重合，所以建立AOI对重合区域进行剪切，以得到研究区域，截图如下:
3. 重采样
因多光谱图像分辨率较低，像元点较大，若要与全色图融合出高清影像需进行重采样来调整像元大小，以达到与高清图一致。
4. 二次剪切
因图为栅格，统一像元后，边缘区必然会有一定的扩展(如下图)，虽说扩展的范围较小，但在科研应用方面不符合要求，故须二次剪切。
5. RGB转HIS
2
2014-5-10
TM图像选取前三层再分别赋予蓝、绿、红三色，转化为HIS格式，如下图:
6. 直方图匹配
将高清图像直方图以标准图像的直方图为标准作变换,使全色光图和HIS图中I层两图像的直方图相同和近似,从而使两幅图像具有类似的色调和反差，以便作进一步的运算。
7. 图像叠加
运用Layer stack功能将全色光高清图和H、S图层进行叠加即所谓的图像融合。它将多波段图层组合到了一起，从而得到新的包含多个有助于研究者使用的多波段影像。
8. IHS转RGB
3
2014-5-10
为了得到真彩色，我们将IHS结果转化为RGB格式，如图:
, 结果分析
1. 多源遥感影像融合法之HIS变换立足于从多波段数据中选取RGB(321)3层，包括了人类视力所能感知的所有颜色，是目前运用最广的颜色系统之一，去除其他空间信息或光谱信息，并按照一定的算法或规则进行处理，将选取的3图层变换为以IHS表示。 下图为TMAtlanta在RGB格式下的6波段(图层)信息，我们选取前3层分别赋予蓝色、绿色、红色，如图所示:
将RGB3层转换为以明度，色调和饱和度为显示格式的HIS图像。由教材可知IHS只是保留全部了TM数据中的蓝色、红色和绿色波段，而删去了近红外、短波红外和热红外波段。 2. 在详细研究图像色调、饱和度与明度关系的基础上，利用图像明度与色调、饱和度有相似平滑这一性质我们将最小二乘原理应用到HIS图像融合中，该方法通过图像平滑度限定条件并用最小二乘准则根据未删除原始数据波段前的空间信息估计(预测)出融合后的色调、饱和度分量信息并运用到融合结果中，这样就避免了旧版IHS变换法中因为直接由高清全色图替换I层空间信息而造成信息丢失的误差，从而能够得到信息比较全的融合图像，这种方法成为改进后的HIS变换法，即Modified HIS Resolution Merge。下图为传统方法(左)和改进后(右)的结果对比图。由图可知改进后的方法所得到的图像较为清晰，减少了色彩偏差。
4
2014-5-10
3. IHS表示亮度(空间信息)、色调和饱和度(光谱信息)。将全色波段的高分辨率图像代替I(明度)图层，融合的图像结果既有高的空间分辨率，又有与原图像相同的色度和饱和度。本实验结果如右图所示。由图可知，该图像保留较清晰的部分为道路和居民点。由TM应用范围表可以分析:前3个图层即RGB波段在其各自应用领域都有相同的应用范围——区分人造地物类型。故可以总结出HIS变换大在保留居民点和道路上有着绝对优势。在目视解译城市道路和居住房分布等方面我们主要运用HIS变换法进行图像融合。
TM影像的主要应用范围
波段序号 波段 主要应用领域
对水体有透射能力，能够反射浅水水下特征，可区分土壤和植被、编蓝色 1 制森林类型图、区分人造地物类型
探测绿色健康植被反射率、可区分植被类型和评估作物长势，区分人绿色 2 造地物类型，对水体有一定的透射能力
3 红色 可测量植被绿色素吸收率，并依次进行植被分类，区分人造地物类型
测定生物量和作物长势，区分植被类型，绘制水体边界、探测水中生4 近红外 物的含量和土壤湿度
5 短波红外 用于探测植物含水量及土壤湿度，区别云与雪
探测地球表面不同物质的自身热辐射的主要波段，可用于热分布制6 热红外 图，岩石识别和地质探矿等方面
融
合
结
果 图
5
2014-5-10
下面是诗情画意的句子欣赏,不需要的朋友可以编辑删除!!
谢谢!!!!!
1. 染火枫林，琼壶歌月，长歌倚楼。岁岁年年，花前月下，一尊芳酒。水落红莲，唯闻玉磬，但此情依旧。
2. 玉竹曾记凤凰游，人不见，水空流。
3. 他微笑着，在岁月的流失中毁掉自己。
4. 还能不动声色饮茶，踏碎这一场，盛世烟花。
5. 红尘嚣 浮华一世转瞬空。
6
2014-5-10
6. 我不是我 你转身一走苏州里的不是我 。
7. 几段唏嘘几世悲欢 可笑我命由我不由天。
8. 经流年 梦回曲水边 看烟花绽出月圆。
9. 人生在世，恍若白驹过膝，忽然而已。然，我长活一世，却能记住你说的每一话。
10. 雾散，梦醒，我终于看见真实，那是千帆过尽的沉寂。
11. 纸张有些破旧，有些模糊。可每一笔勾勒，每一抹痕迹，似乎都记载着跨越千年万载的思念。
12. 生生的两端，我们彼此站成了岸 。
13. 缘聚缘散缘如水，背负万丈尘寰，只为一句，等待下一次相逢。
14. 握住苍老，禁锢了时空，一下子到了地老天荒
15. 人永远看不破的镜花水月，不过我指间烟云 世间千年,如我一瞬。
16. 相逢一醉是前缘，风雨散，飘然何处。
17. 虚幻大千两茫茫，一邂逅，终难忘。相逢主人留一笑，不相识，又何妨。
18. 天下风云出我辈,一入江湖岁月催;皇图霸业谈笑间,不胜人生一场醉。
19. 得即高歌失即休，多愁多恨亦悠悠，今朝有酒今朝醉，明日愁来明日愁。
20. 直道相思了无益，未妨惆怅是清狂。
21. 看那天地日月，恒静无言;青山长河，世代绵延;就像在我心中，你从未离去，也从未改变。
22. 就这样吧，从此山水不相逢。
23. 人天自两空，何相忘，何笑何惊人。
24. 既不回头,何必不忘。 既然无缘,何须誓言。 今日种种,似水无痕。 明夕何夕,君已陌路。
25. 有缘相遇,无缘相聚,天涯海角,但愿相忆。有幸相知,无幸相守,苍海明月,天长地久。
26. 相见得恨晚，相爱的太慢，进退让我两难。缘过了远分，缘过了聚散，是否回头就能够上岸
27. 天凉了,凉尽了天荒 地老了,人间的沧桑，爱哭了,这么难舍 心都空了,想放不能放。天亮了,照亮了泪光 泪干了,枕边地彷徨
28. 心微动奈何情己远.物也非,人也非,事事非,往日不可追
29. 渺渺时空,茫茫人海,与君相遇,莫失莫忘。
30. 如果换我先开口，日子是否还一样细水长流
31. 也许是前世的姻 也许是来生的缘 错在今生相见 徒增一段无果的恩怨
32. 人道海水深，不抵相思半。海水尚有涯，相思渺无畔。
7
2014-5-10
33. 醉眼看别人成双作对，
34. 无人处暗弹相思泪。
35. 终于为那一身江南烟雨覆了天下，容华谢后，不过一场，山河永寂。 36. 千秋功名，一世葬你，玲珑社稷，可笑却无君王命。 37. 凤凰台上凤凰游，负约而去，一夜苦等，从此江南江北，万里哀哭。 38. 嗟叹红颜泪、英雄殁，人世苦多。山河永寂、怎堪欢颜。 39. 风华是一指流砂，苍老是一段年华。
40. 夜雨染成天水碧。有些人不需要姿态，也能成就一场惊鸿。 41. 你要记得，紫檀未灭，我亦未去。
42. 谁在岁月里长长叹息。
43. 汉霄苍茫，牵住繁华哀伤，弯眉间，命中注定，成为过往。 44. 红尘初妆，山河无疆。 最初的面庞，碾碎梦魇无常，命格无双。 45. 江南风骨，天水成碧，天教心愿与身违。
46. 山河拱手，为君一笑 。
47. 如是颠簸生世亦无悔。
48. 荏苒岁月覆盖的过往，白驹过隙，匆匆的铸成一抹哀伤。 49. 那被岁月覆盖的花开，一切白驹过隙成为空白。
50. 褪尽风华，我依然在彼岸守护你。
51. 那些繁华哀伤终成过往，
52. 请不要失望，平凡是为了最美的荡气回肠。
53. 你的路途，从此不见我的苍老。
54. 长歌当哭，为那些无法兑现的诺言，为生命中最深的爱恋，终散作云烟。 55. 随你走在天际，看繁花满地。
56. 我自是年少，韶华倾负。
57. 你要记得，那年那月，垂柳紫陌洛城东。
58. 苍茫大地一剑尽挽破，何处繁华笙歌落。 59. 寄君一曲，不问曲终人聚散。
60. 谁将烟焚散，散了纵横的牵绊;听弦断，断那三千痴缠。61. 清风湿润，茶
烟轻扬。重温旧梦，故人已去。
1. 水滴虽小，却可以折射出太阳的光彩。
2. 梦落三千尺愁深似海，繁华遗落散满地。记忆轮回里，我举杯，在奈何桥
上满口饮尽。
8
2014-5-10
3. 人生没有轮回，就像花，人活一世，花开一季、人生如花，花似梦。
4. 生活的苦涩和美好给了我对人生的领悟，如今，千山万水走遍，我发现自己再也不愿离开文学的蓝天，再也不愿离开那个让我痴迷的文学舞台。
5. 在烟雨红尘中，轻拾季节花瓣飘落的音符，组成美妙曲符，然后，倚在时光的路口，撷一缕明媚，许自己一份唯美的怀想，与快乐、浪漫相约，闲淡清欢。
6. 未经历坎坷泥泞的艰难，哪能知道阳光大道的可贵;未经历风雪交加的黑夜，哪能体会风和日丽的可爱;未经历挫折和磨难的考验，怎能体会到胜利和成功的喜悦。挫折，想说恨你不容易
7. 燕子斜飞人家，炊烟零乱，柳絮飘飘，弥漫了山里人家。
8. 这样知解自己的生命即使是心灵空荡我也无所畏惧
9. 中秋之曰不可能岁月明如水，偶然的暗淡，恰似镜子的背后之面，有所缺憾，人生才会是积翠如云的空濛山色。
10. 在经受了失败和挫折后，我学会了坚韧;在遭受到误解和委屈时，我学会了宽容;在经历了失落和离别后，我懂得了珍惜。
11. 曾经盛开的蔷薇，虽经风吹雨打，但和着微风，还有屡屡暗香飘过。
12. 我只希望，不管三年，五年，或是十年以后。某一天，我们相遇，还能相认，你大喊一声，我想死你了。那一刻，我定会泪流满面。我们是朋友，永远的朋友。
13. 最爱的未必适合在一起，相爱是让彼此做自己。
14. 时间断想，时间不断。流逝，像是水，可弯可直，像是风，可柔可刚。
15. 如果说人生是一望无际的大海，那么挫折则是一个骤然翻起的浪花。如果说人生是湛蓝的天空，那么失意则是一朵飘浮的淡淡的白云。
16. 云层雾气，缠着几户古木人家，清新自然，如诗如画。
17. 我喜欢你，只是一个现在;我爱你，却是一整个未来。
18. 夜雨染成天水碧。有些人不需要姿态，也能成就一场惊鸿。
19. 再大的风不会永不停息，在浓得雾不会经久不散，风息雾散仍是阳光灿烂。
20. 牵着时光的衣襟，走进芳菲五月，轻轻地将春光拥入怀中，于一抹素白流韵中，弹奏一曲江南的婉韵，把盏，将似水流年浅斟轻酌。
21. 我穿越轮回而来，在奈何桥相思盈袖，凄然守候。莫落泪，纵若水落三千尺东流，云动八万里西散，我依旧会化身城碟，翩翩起舞跨过奈何桥与你相会。
22. 如果我爱你，我就会理解你，通过你的眼睛去看世界。我能理解你，是因为我能在你身上看到我自己，在我身上也看到了你。
23. 似乎风在转向，送走了缓缓袭来的味道，又将刚刚溜走的风，静静地换回来。
9
2014-5-10
24. 生活告诉我，童话只不过是小孩子幻想的游戏。
25. 人生就像穿着一件长满虱子的华丽睡袍，外表美丽，而内心却充满了干啊和恐慌。
26. 必须用另一种真实方式来代替时光里已经逝去的东西
27. 岁月，依一抹浅香于心间，看年华向晚，闻花香送暖。给时光一个浅浅的回眸;给自己一份微笑从容。沉淀，馨香;念起，温暖。
28. 人生的起起落落间，总会有一些情怀需要安静回味;总会有一些伤痛需要独自体会;总会有一段路需要一个人走;总会有一些事需要坦然面对。
29. 疏影横斜水清浅，暗香浮动月黄昏。
30. 心若没有栖息的地方，到哪里都是在流浪。
31. 今后，我会从尘世中的纷争走出，远离喧嚣，把岁月打磨成诗，让自己的文字静如睡莲，动如涟漪，无论何时都能描绘成美丽的水墨丹青。
32. 全是理智的心，恰如一柄全是锋刃的刀，它叫使用它的人手上流血。——泰戈尔
33. 我们都不擅长表达，以至于我们习惯了揣测。去肯定，去否定，反反复复，后来我们就变得敏感而脆弱。
34. 心心念念的往事、曾经深爱过的人、年少琐碎的过往，它们就像缠绕之间的一阵风，来的缱绻，去的时候让人来不及挽留。
35. 如果在乎的没有那么多，想要的没有那么多，生活便会简单得很多。
36. 在极度的喧嚣中，独自微笑独自平静是憾，落花是美的，淡淡的书香，淡淡的花香，淡淡的馨香。
37. 曾芬芳过的那片土地，幸福的花儿虽早已凋谢，只留下风雨吹打的痕迹。
38. 辗转半世红尘，缘去缘灭，空留满池伤痕。雨花迟落，霜雪纷飞，池水泛冰，已益处月的苍凉。
39. 一条古道，一匹瘦马，一个人影，被落日的余晖缓缓拉长。
40. 我们人生的大幕才刚刚拉启:刀光剑影，英雄本色;是非恩怨，儿女情常。
41. 我们要去流浪，虔诚地定格住每一寸记忆;我们要去成长，潇洒地忘却掉每一条纹路。
42. 嗅着昨日芬芳遗留的气息，寻寻觅觅，仍不见踪迹。邂逅了一场烟火，终还是那般凄凉。迷失的夜晚，点缀了无数颗孤单的星星，不知道那是否有属于我的一颗。
43. 像这样轻飘飘的日子和平平静静的心情，也算是生活中的一种享受吧。
44. 想着远方的你，绝美的笑容，只为你一个人展露，那一泓羞涩的笑容，悄悄。
10
2014-5-10
45. 细碎的声音，如羞涩的蓓蕾，夜暮花影，轻浅六月，寂寂流年，拢一阙清绝，归隐在宋词里。
46. 夕阳沉落在海水深处却不见浪花翻滚，淡淡的只留下一个让人沉思的背影。落雨是晚风中的殇，带着晨曦的翘首滑落最后的伤痕!雨尽含羞，淡抹嫣红!
47. 无影击碎了泪水，岁月在那个光年划下的痕迹原来是一刀一刀地刻在了我的心上。
48. 我收拢了梦想的翅膀，我停却了信念的脚步，却再也作不回曾经的那一天。
49. 我宁愿用尽此生，为那些尘世的硝烟尘雾，潸然泪下，为菩提落花，为世间繁华。
50. 阳光依旧在，我们穿越光影，沿着历史的足迹继续前行，创造美好生活，走向美丽明天!
51. 洋溢着春日的微笑，坚强了外表，却虚伪了内心，脆弱了，是不敢触及的。
52. 也许，就在那一刻里，梦境还在，柔情亦在。
53. 一个人的戏，自己独自导演，诠释精彩。在剧中尽情释放着自己的喜怒哀乐。笑得凄然绝美;哭得肝肠寸断。
54. 但生命中被你刻上痕迹的那些岁月无法抹去。
55. 往事悠然一笑间，不必空忧。我们一路走来，只是为了告别往事，走入下一段风景。倘若让忧伤填补了生命的空白，就真的是亵渎了生命。
56. 人生只有回不去的过去，没有过不去的当下。上帝只会给你过得去的坎，再不好过的生活，再难过的坎，咬咬牙，也就过去了。
57. 我一直以为山是水的故事，云是风的故事，你是我的故事。可是却不知道，我是不是你的故事。
58. 生命并不是一场竞赛，而是一段旅程。如果你在途中一直都试图给他人留下深刻印象，超过别人，那你就浪费了这段旅程。
59. 比如新的朋友新的感情新的思绪我想要知道的
60. 我以为我已经将爱情忘记，将你忘记。可是有一天，我听到一首歌，我的眼泪就出来了。因为这首歌，我们曾一起听过。
61. 忍花开花落，云卷云舒，品人生似棋。
62. 我离开你这一种信仰又会以怎样全新的姿势去面临更深沉的挑战
63. 人生路，路迢迢，谁道自古英雄多寂寥，若一朝，看透了，一身清风挣多少。
64. 只有夕阳站在那里。灵魂像无数的雪花飘过，光明闪烁，渐渐清醒。
65. 终于为那一身江南烟雨覆了天下，容华谢后，不过一场，山河永寂。
66. 荏苒岁月覆盖的过往，白驹过隙，匆匆的铸成一抹哀伤。
11
2014-5-10
67. 忘川水不枯，记忆不散;奈何桥不断，思卿不弃;今夕，彼岸花又放，佳期约又到，我轮回践约而来，等你归来。红尘路上，伊人在否?
68. 十年生死两茫茫，不思量，自难忘，千里孤坟，无处话凄凉，纵使相逢应不识，尘满面，鬓如霜。
69. 天空飘过一朵云，有时是晴，有时是阴。但白昼终归还是白昼。
70. 我知道回不去，但还是会想念会回忆会心疼到无法自拔。
71. 天空不曾留下鸟的痕迹，但是我已飞过、在大地上画满窗子，让所有习惯黑暗的眼睛都习惯光明。
72. 人生首先要是望远镜，看远;再就是显微镜，看细;接下来是放大镜，看透;其次是太阳镜，看淡;最后是哈哈镜，笑看生活。
73. 我不是公主，也不会有等待救赎我的王子。
12
erdas实验报告
前 言一、实习任务:练习使用 ERDAS IMAGINE 软件、撰写实习报告。二、实习时间:2011 年 9 月三、实习地点:第四教学楼五楼机房四、实习目的: 1、了解 ERDAS IMAGINE 的应用基础及应用领域。 2、掌握图像校正、拼接、投影变换、分幅裁剪、融合等预处理操作，图像增强，图像分类等遥感数字图像处理基本内容。 3、通过对 ERDAS 软件的学习和实习操作，了解遥感图像处理的基本原理、流程以及软件系统的基本构成和功能。 4、加深对所学课程原理的理解，为从事相关项目的研究和开发奠定基础。五、实习内容: 1、遥感图像处理软件概述及 ERDAS 软件基本操作。 2、遥感图像的输入/输出、波段组合及图像显示。 3、遥感图像预处理，主要包括图像校正、分幅裁剪、拼接、投影变换、融合等操作。 4、遥感图像增强。 5、遥感图像分类，分为监督分类和非监督分类。 实习一、对 j50e023013.img
和 j50e024013.img 进行几何校正 1、加载图像文件: (1)在 ERDAS 图标面板菜单条选择 Main—Start Image Viewer 命令，打开 Viewer 窗口 Viewer1。 在 (2) Viewer1 菜单条选择 File—Open—Raster Layer 命令，打开 Selectlayer to Add 窗口， 选 选择需要校正的图像 j50e023013.img。 择 Raster Options标签，选择 Display as 下拉条，选择 Gray Scale在 Layer 中输入 2:;选中 Fitto Frame 复选框，以 使添 加的 图像全 幅显 示。 单击 OK ，加载 需 要校正 的j50e023013.img 图像。 2、启动几何校正模块: 在 Viewer1 菜单条选择 Raster Geometric Correction 命令，打开选择几何校正模型(Set Geometric Model)对话框，选择多项式变换模型 Polynomial，单击 OK。在 Polynomial Model Properties 中定义多项式次方 Polynomial Order为 2，单击 Apply 按钮应用设置。单击 Close 按钮关闭当前对话框，打开 GCP ToolReference Setup 对话框。 3、采集地面控制点: (1)在 GCP Tool Reference Setup 对话框中选择 Keyboard 键盘输入，修改投影信息。单击 OK 后自动启动 Reference Map Information 组合对话框，输入 X、Y 的坐标。 (2)点击?求解误差，使误差保持在 0.5 范围内。 4、图像重采样: (1)用 Geo Correction
Tools 窗口进行重采样，选择 Image Resample 按钮。输出图像(Output File)文件名及路径，为 023013.img。 (2)选择重采样方法——最邻近采样(nearest neighbor) (3)定义输出图像范围:在 ULX、ULY、LRX、LRY 微调框中分别输入需要的数值，现在采用默认值。 (4)定义输出像元大小，X 值 30，Y 值 30，一般与数据像元大小一致。 (5)设置输出统计中忽略零值，即选中 Ignore Zero in Stats 复选框。 (6)单击 OK 按钮执行重采样。 (7)对 024013.img 的几何校正方法步骤同 023013.img 的几何校正: 实验二、对 023013.img、024013.img 图像进行裁剪 对两张图像进行裁剪运用的方法是不规则分幅裁剪。不规则分幅裁剪是指裁剪图像的边界范围是任意多边形，不能通过左上角和右下角两点的坐标确定裁剪范围，而必须事先设置一个完整的闭合的多边形区域，可以利用 AOI 工具创建裁剪多边形，然后利用分幅工具进行分割。 实现步骤如下: (1)打开要裁剪的图像 023013.img，在 Viewer 图标面板菜单条选择 AOI—Tools 菜单，打开 AOI 工具条。 (2)应用 AOI 工具绘制多边形 AOI，将多边形 AOI 保存在 023013_aoi.aoi文件中。 在 (3) ERDAS 图标面板菜单条选择 Main—Data Preparation—Subset Image命令，打开 Subset 对话框。 (4)选择处理图像文件 023013.img。 (5)输出文件名称为 023013_sub.img并设置存储路径。 (6)单击 AOI，打开 Choose AOI 对话框，选择 AOI 来源为 AOI File 中的023013_aoi.aoi 文件。 (7)输出数据类型为
Unsigned 8 BitContinues。输出统计忽略零值，选中 Ignore Zero in Stats 复选框。 (8)设置输出波段—1:7。 (9)单击 OK 按钮关闭 Subset 对话框，执行图像裁剪。 023013.img 和 024013.img 裁剪后图像如下: 实验三、023013_sub.img 与 024013_sub.img 的拼接此处采用多波段图像拼接主要步骤如下:1、启动图像拼接工具: 在 ERDAS 图标面板菜单条选择 Main—Data Preparation—Mosaic Images—Mosaic Tool 命令，打开 Mosaic Tool 对话框。2、加载拼接图像: (1)选择 Main—Data Preparation—Mosaic Images—Mosaic Tool 工具对话框，选择 Display Add Images Dialog 按钮，打开 Add Images 对话框。选择需要进行拼接的两个图像。 (2)再选择 Images Area Options 标签，进入 Images Area Options 对话框，进行拼接图像的选择。 (3)这里选择计算活动区 Compute Active Area 按钮，并单击 Set 打开Active Area Options 对话框设置如下参数: Select Search Layer:
指定哪个图层用于活动区的选择。 Background Value Range:背景值范围，即根据 from，to 设置某一光谱段或光谱值为背景值，在运行拼接过程中落入该光谱范围内的图像不参与拼接运算。 Boundary Search Type:边界搜索类型，包括 Corner 和 Edge 选项。选择Corner 时可以对 Crop Area 进行设置，将对输入影响进行裁剪。 4单击 OK，加载两幅图像。3、图像叠置组合: (1)在 Mosaic Tool 工具条选择 Set Model for Input Images 按钮，进入图像设置模式状态。Mosaic Tool 工具条会出现与该模式对应 的 调整图像跌至次序的编辑按钮 (2 )利用工具 库对图像进行上下调整，确定拼接方 案。 拼接次序调整好后，意味着重叠区将随之确定，不同的图像叠置组合，用于拼接的重叠区域也不同，查看重叠区域采用的是 Mode for Intersection 按钮4、图像匹配设置: 在 (1) Mosaic Tool 工具条选择 Display Color Corrections 按钮 ，打开色彩校正Color Corrections对话框。 (2)选中 Use Histogram Matching 按钮，单击 Set，打开 HistogramMatching直方图匹配执行影响的色调调整。 (3)匹配方法Matching Method为 Overlap Areas，即只利用叠加区直方图进行匹配。直方图类型 Histogram Type为 Band by Band，即分别从红绿蓝三个波段进行灰度的调整如果是多波段，则表示逐波段进行一对应的灰度。 4单击 OK 按钮，保存设置返回到 Color Corrections 对话框，在 ColorCorrections 窗口中再次单击 OK 按钮退出。 5在 Mosaic Tool 工具条选择 Set Mode for Intersection 按钮 ，进人设置图像关系模式的状态。 6在 Mosaic Tool 工具条选择叠加函数Set Overlap Function按钮 ( 7 ) 设 置 叠 加 方 法 Intersection Method 为 无 剪 切 线 No OutlineExists，重叠区像元灰度计算Select Function为均值Average，即叠加区各个波段的灰度值是所有覆盖该区域图像灰度的均值。 8单击 Apply 按钮应用设置，单击 Close 按钮关闭 Set Overlap Function对话框。5、运行 Mosaic 工具: 在 (1) Mosaic Tool 工具条选择输出影响模型Set MODE FOR Output Images按钮 ，进入输出模式设置状态。选择 Run the Mosaic Process to Disk 按钮 ，打开 Output File Name 对话框。 (2)输出文件名为 023013024013_ mosaic.img 选择 Output Options 标签，选中忽略零值(Stats Ignore Value)复选框。 (3)单击 OK 按钮，关闭 Run Mosaic 对话框，运行图像拼接。 023013_sub.img 与 024013_sub.img 的拼接图 023013024013_ mosaic.img 如下:实验四、用 023013024013_ mosaic.img 对 spotapan01422.img 的配准 1、加载图像文件: (1)在 ERDAS 图标面板菜单条选择 Main—Start Image Viewer 命令，打开 Viewer 窗口 Viewer1。同步骤 1 打开
一个新的 Viewer 窗口 Viewer2. 在 (2) Viewer1 菜单条选择 File—Open—Raster Layer 命令，打开 Select 选择需要校正的图 spotapan01422.img。择 Raster Optionslayer to Add 窗口， 选标签，选择 Display as 下拉条，选择 Gray Scale在 Layer 中输入 2:;选中 Fitto Frame 复选框，以使添加的图像全幅显示。单击 OK，加载需要校正的图像spotapan01422.img。 在 (3) Viewer2 菜单条选择 File—Open—Raster Layer 命令，打开 Selectlayer to Add 窗口，选择参考图像 023013024013_ mosaic.img。 2、启动几何校正模块: (1)在 Viewer1 菜单条选择 Raster Geometric Correction 命令，打开选择几何校正模型(Polynomial Model Properties)对话框。 在 (2) Polynomial Model Properties 中定义多项式次方 Polynomial Order为 2，单击 Apply 按钮应用设置。单击 Close 按钮关闭当前对话框，打开 GCP ToolReference Setup 对话框。 3、启动控制点工具: 在 即为 Existing Viewer (1) GCP Tool Reference Setup 窗口选择采点模式，按钮。单击 OK 按钮关闭该窗口，打开 Viewer Selection Instruction 指示器。 ( 2 ) 鼠 标 点 击 显 示 作 为 地 理 参 考 图 像 023013024013_ mosaic.img 的Viewer2 窗口，打开 Reference Map Information 对话框，显示参考图像的投影信息。 (3)单击 OK 按钮，控制点工具被启动，进入控制点采集状态。采集过程如下。4、GCP 具体采集过程如下: (1)在 Viewer1 移动关联方框，寻找特征明显的地物点，作为输入 GCP，并在 Viewer3 中单击确定相对应的点，GCP 数据表将记录一个输入 GCP，包括其编号、标识码、X 坐标、Y 坐标。 (2)为使 GCP1 容易识别，单击 GCP 数据列表的 Color 列 GCP1 对应的空白处，在弹出的颜色列表中选择比较醒目的颜色。 (3)在 Viewer2 移动关联方框，寻找特征明显的地物点，作为参考 GCP。 (4)为使 GCP1 容易识别，单击 GCP 数据列表的 Color 列 GCP1 对应的空白处，在弹出的颜色列表中选择比较醒目的颜色。 (5)不断重复步骤(1)—(4)，采集若干 GCP，直到满足所选定 6 个控制点。 (6)点击?求解误差，使误差保持在 0.5 范围内。 每采集一个 Input GCP，系统就自动产生一个参考控制点(Ref. GCP)通过移动 Ref. GCP 可以逐步优化校正模型。5、图像重采样: (1)用 Geo Correction Tools 窗口进行重采样，选择 Image Resample 按钮。输出图像(Output File)文件名及路径，为 01422.img。 (2)选择重采样方法——最邻近采样(nearest neighbor) (3)定义输出图像范围:在 ULX、ULY、LRX、LRY 微调框中分别输入需要的数值，现在采用默认值。 (4)定义输出像元大小，X 值 30，Y 值 30，一般与数据像元大小一致。 (5)设置输出统计中忽略零值，即选中 Ignore Zero in Stats 复选框。 (6)单击 OK 按钮执行重采样。 023013024013_ mosaic.img 对 spotapan01422.img 的配准图 01422.img 如下:最后校正成图: 实验五、用 01422.img 对 ETM.img 进行配准 1、加载图像文件: (1)在 ERDAS 图标面板菜单条选择 Main—Start Image Viewer 命令，打开 Viewer 窗口 Viewer1。同步骤 1 打开一个新的 Viewer 窗口 Viewer2. 在 (2) Viewer1 菜单条选择 File—Open—Raster Layer 命令，打开 Selectlayer to Add 窗口，选择需要校正的图像 ETM.img。选择 Raster Options 标签，选择 Display as 下拉条，选择 Gray Scale在 Layer 中输入 2:;选中 Fit toFrame 复选框，以使添加的图像全幅显示。单击 OK，加载需要校正的图像 。 在 (3) Viewer2 菜单条选择 File—Open—Raster Layer 命令，打开 Selectlayer to Add 窗口，选择参考图像 01422.img。 2、启动几
何矫正模块: (1)在 Viewer1 菜单条选择 Raster Geometric Correction 命令，打开选择几何校正模型(Polynomial Model Properties)对话框。 在 (2) Polynomial Model Properties 中定义多项式次方 Polynomial Order为 2，单击 Apply 按钮应用设置。单击 Close 按钮关闭当前对话框，打开 GCP ToolReference Setup 对话框。 3、启动控制点工具: 在 即为 Existing Viewer (1) GCP Tool Reference Setup 窗口选择采点模式，按钮。单击 OK 按钮关闭该窗口，打开 Viewer Selection Instruction
指示器。 (2)鼠标点击显示作为地理参考图像 01422.img 的 Viewer2 窗口，打开Reference Map Information 对话框，显示参考图像的投影信息。 (3)单击 OK 按钮，控制点工具被启动，进入控制点采集状态。 4、GCP 具体采集过程如下: (1)在 Viewer1 移动关联方框，寻找特征明显的地物点，作为输入 GCP，并在 Viewer3 中单击确定相对应的点，GCP 数据表将记录一个输入 GCP，包括其编号、标识码、X 坐标、Y 坐标。 (2)为使 GCP1 容易识别，单击 GCP 数据列表的 Color 列 GCP1 对应的空白处，在弹出的颜色列表中选择比较醒目的颜色。 (3)在 Viewer2 移动关联方框，寻找特征明显的地物点，作为参考 GCP。 (4)为使 GCP1 容易识别，单击 GCP 数据列表的 Color 列 GCP1 对应的空白处，在弹出的颜色列表中选择比较醒目的颜色。 (5)不断重复步骤(1)—(4)，采集若干 GCP，直到满足所选定 6 个控制点。 (6)点击?求解误差，使误差保持在 0.5 范围内。 每采集一个 Input GCP，系统就自动产生一个参考控制点(Ref. GCP)通过移动 Ref. GCP 可以逐步优化校正模型。 5、图像重采样: (1)用 Geo Correction Tools 窗口进行重采样，选择 Image Resample 按钮。输出图像(Output File)文件名及路径，为 ETM配准.img。 (2)选择重采样方法——最邻近采样(nearest neighbor) (3)定义输出图像范围:在 ULX、ULY、LRX、LRY 微调框中分别输入需要的数值，现在采用默认值。 (4)定义输出像元大小，X 值 30，Y 值 30，一般与数据像元大小一致。(5)设置输出统计中忽略零值，即选中 Ignore Zero in Stats 复选框。(6)单击 OK 按钮执行重采样。 01422.img 对 ETM.img 的配准图 ETM配准.img 如下:校正成图为: 实验六、01422.img 和 ETM配准.img 的融合 (1)在 ERDAS 图标面板菜单条选择 Main—Image Interpreter—SpatialEnhancement—Resolution Merge 菜单，打开 Resolution Merge 对话框。 确 (2) 定 高分辨率输入文件 01422.img，确定多光谱输入文件 ETM配准.im，定义输出文件 spota—ETM.img。 (3)选择融合方法Method—主成分变换(Principle Component) (4)选择重采样方法(Resample Techniques)—Bilinear Interpolation双线性插值。 (5)输出数据选择(Output Options)为 Stretch Unsigned 8 bit。 (6)输出波段选择为 Select Layers 1:7。 (7)单击 OK 按钮，关闭 Resolution Merge 对话框，执行分辨率融合。 01422.img 和 ETM配准.img 融合图像 spota—ETM.img 如下: 实验七、对 spota—ETM.img 的分类 1、启动非监督分类模块，选择输入、输出影像: (1)在 ERDAS 面板工具中选择 DATAPrep—Unsupervised Classification命令，打开非监督分类对话框。 (2)选择 Input Raster File 和 Output Cluster Layer Filename 设置被分类的图像和分类结果，并选择 Output Signature Set 产生一个模板文件。 2、设置初始参数，执行非监督分类: 在非监督分类对话框中分别设置聚类选项和处理选项: (1)选择 Classifier Unsupervised Classification 命令，打开非监督对话框，设定输入输出数据，设置聚类选项，确定初始聚类方法和分类数: ?设置聚类方法
为 Initial from Statistics确定初始分类数为 10，设置分类颜色属性为 4、3、2。 ?设置 ISODATA 的统计参数 Initializing Means Along 和 Scaling Range分别为计算初始均值所沿的轴和度量范围。 (2)设置处理选项，确定循环次数和阈值 ?在 Unsupervised Classification 面板中的 Processing options 栏选定最大循环次数为 24。 ?设置循环收敛阈值，其值是指两次分类结果相比保持不变的像元所占的最 ，目的在于避免 ISODATA 无限循环下去。单击 OK，执行非监督分大比例(0.95)类。 对 spota—ETM.img 的分类图像如下: 实习心得 通过这次的实习，使我更好的掌握了课本知识，对遥感原理和 ERDAS 软件也有了进一步的了解，一开始觉得它不是一个轻易学好的软件，觉得用处不是很大。但是经过一段时间的学习以后，发现它没有想象的那么难，而且应用范围很广，发展前景也相当可观。未来几年将会被应用与各个领域甚至融入我们的日常生活中，它会使我们的生活更加方便、快捷、省去很多不必要的麻烦。 另外在实习过程中遇到问题后积极寻找解决问题的办法，增强了自己独立思考解决问题的能力，对知识的理解和应用也更透彻。 而且通过这次学习的机会，使我们多把握一项技能，就多一次选择，在竞争激烈的社会中获得多一点的机会。 最后，实习也锻炼了我认真细心的品质，不骄不躁，坚持到底，提高了我的专业知识，使我收获颇丰。
erdas实验报告
遥感实验报告
学 院: 资源与环境工程学院 专 业 班 级: 资源环境与城乡规划管理0701班 姓 名: 陈吉 学 号: 200701122011 完 成 学 期: 2009-2010-12 指 导 教 师: 胡林
实验一 数据输入
一(实习目的:掌握TM图像数据输入的主要方法。
二(实习内容:主要包括单波段TM图像数据输入、多波段组合文件的生成。
1. 单波段二进制图像数据输入
在ERDAS图标面板工具条中，点击——打开Import/ Export对话框，如图(1)。
(1)选择Import 。选择数据输入类型(Type)为普通二进制:Generic Binary 。选择数据输入媒体(Media)为文件:File。
(2)Input File下确定输入文件路径及文件名:D盘/实验数据/实验一TM数据/Band1.dat。
(3)Output File下确定输出文件路径及文件名:D盘/实验结果/1.img
图(1)
(4)单击OK。 打开Import Generic Binary Data对话框，如图(2)所示。
(5)Data Format选择框里定义数据格式:BSQ。Data Type为Unsigned 8 Bit
(6) 数据文件行数(Row):5733。数据文件列数(Cols):6773。文件波段数量(Bands):1。(7)保存参数设置(Save Option):*.gen。单击OK退出Save Option File。
图(2)
(8)OK执行输入操作。出现进程状态条
100%后单击OK完成数据输入。
(9)重复上述过程，可依此将多波段数据全部输入，转换为 .IMG文件。 2、组合多波段数据
(1)在ERDAS图标面板工具条中，点击Interpreter>Utilities>Layer Stack出现Layer Selection and Stacking的对话框。如图(3)所示。
实验三 对比度变换
一(实验目的:了解直方图的含义， 掌握灰度变换的意义，掌握利用软件对影像进行灰度变换，掌握线性变换的原理。
二(实验内容:
1. 遥感影像的直方图
在Viewer1中选择file>open>raster layer.在select layer to add 对话框中选择D盘/实验数据/实验三数据/lanier.img。选择Raster Option选项卡No strech打勾，单击OK显示图像。打开Utility / layer info。出现对话框如下左图所示。选择Histogarm选项卡，观察影像直方图，如下右图所示
在中点击上下按钮，换一个工作波段重复, 观察影像色调与对应直方图的变化。
2.对比度拉伸:
(1)在Viewer1中选择Raster>contrast> General Contrast。弹出Contrast Adjust对话框，Method后选择Histogram Equalization ，单击Apply。直方图均衡化后图像发生变化。(如下左图)。单击BreakPoint，出现一对话框如下右图:
图(3)
(2)在Output File:*.img下选择要加载的单波段IMG图像(如D盘/实验结果/1.img)。单击Add依此选择并加载，把七个波段图像都加载。Output File:*下确定输出图像路径及名称D盘/实验结果/1234567.img 。单击OK执行并完成波段组合。
(3)在Viewer1中打开多波段组合后的图像1234567.img。如图(4)所示:
图(4)
实验二 遥感图像几何纠正
一( 实习目的 掌握几何纠正的方法，学会找控制点。
二(实习内容 采集控制点、检查点，计算变换参数，灰度重采样，纠正精度评定。
1.对TM影象进行多项式纠正。
(1)在Viewer1中选择file>open>raster layer.在select layer to add 对话框中选择D盘/实验数据/实验二几何纠正数据/tmatlanta.img。单击OK打开显示图像。单击工具条上viewer按钮对话框中选择classic viewer单击OK，创建viewer2窗口。同理在viewer2窗口中显示
/实验数据/实验二几何纠正数据/panatlanta.img的图像。 D盘
(2)在Viewer1窗口菜单中选择Raster> Geometric Correction调用几何纠正模型。在对话框中选择多项式模型Polynomial，然后单击OK。Polynomial Order 后输入2，单击apply。关掉对话框。出现GCP Tool Reference Setup对话框，直接点击OK。出现Viewer selection Instructions对话框，要求选择参考图像，所以将鼠标在Viewer2中图像上单击，出现Reference Map Information对话框。该对话框显示参考图像panatlanta.img的投影信息，点击OK，启动GCP TOOL。如图(一)所示:包括两个放大的视窗Viewer3，Viewer4，分别是两个主视窗Viewer1，Viewer2中两个关联方框的放大框。以及GCP TOOL:包含菜单和工具图标，图标下面显示两幅图像坐标数据。
图(一)
(3)在GCP TOOL中点击，进入控制点选择状态，光标显示为十字，表示可以选择同名点，在Viewer3中明显点上点左键，在GCP 数据记录表中显示该点的图像坐标。在Viewer2中移动连接框，到与Viewer3中对应的同名点位置，在GCP TOOL中点击，在Viewer4中明显点上点左键，系统在GCP 数据记录表中自动显示该点的地面坐标。在主窗口和放大窗口中系统会显示GCP #1，表示选择的第一对同名点，其颜色可以在GCP TOOL中重新设定。重复以上步骤，选取第2、3、4、5、6个控制点，第七个及以上控制点的选择:只需在Viewer1中选择一个明显点，其同名点会在Viewer2中自动显示，如果位置有误差，可以用鼠标左键调整点位，确保选择的是同名点。
全部控制点平差后的精度显示在右上角，如下图所示
因为Total小于一个像元，满足纠正精度要求，所以继续后面的重采样。 (,)为了验证控制点的精度，进行检查点采集。
首先选择GCP TOOL菜单栏上Edit/Set Point Type/Check设置点的类型为Check，然后按照前面选择控制点的方法选择检查点。选四个控制点，点击，在相应的栏内显示检查点的残差。如下图所示:
Total小于一个像元，所以原先选择的8个控制点的精度满足要求。
(5)图上控制点，检查点选取如下图所示:
(6)在中单击第一个图标显示Polynomial Model Projection对话框，选择Transformation选项卡，显示多项式中a,b的值。见下图:
至此就完成了多项式纠正的参数计算部分。退出该模块时，系统会提示对有关文件做好保存。
2.灰度重采样
选择灰度重采样模块，出现重采样对话框:在Output File下选择保存路径，并以jiuzheng为文件名保存纠正以后的图像，单击ok。出现重采样的状态条:
当Percent Down为100%时，说明重采样已经完成，单击OK。
3.实习的成果分析
(1)在Viewer1中打开jiuzheng图像，通过与Viewer2中Spot图像关联起来检查两边的位置是否对应如果两边对应，说明纠正图像满足要求。
(2)在Viewer1中的图像上点击鼠标右键，选择Geo.Link/Linlink...出现一对话框提示在要连接的窗口中点击鼠标，将两个窗口关联起来。在Viewer2中单击鼠标左键。过程如下图所示:
4(实验体会
(1) 选取控制点的时候 一定要掌握选取原则:(1)选取图像上易分辨且较精细的特征点，在本实验中有道路交叉口(四处:点1、3、5、7)，河流弯曲处(一处:点6)，明显地物(一处:点4)，建筑区轮廓转角处(两处:点2、8)。(2)控制点一定要分布均匀，不然纠正后图像变形大。
(2) 做实验要有耐心，做错了要善于总结，这样重做的时候才不会再犯错。
2)同理在弹出的Contrast Adjust对话框，Method后选择Standrad Deviations方(
法，观察影像显示结果。
(3) 重新打开实验三数据，在Viewer1中选择Raster>contrast> breakpoint观察对三个波段影像的直方图特点，如下所示:
利用分段线性变化，通过增加\删除\移动断点来调整三个波段影像的直方图,使之分布均匀\或接近正态分布。结果如图所示， 点击Apply All，影像如下图所示
3.实验总结:
掌握线性变化的原理是做此实验的关键，把直方图调整到分布均匀\或接近正态分布有点困难。
