1、遥感实习报告 小实验实验一、典型地物的光谱反射特征一、 实验目的熟悉ENVI软件提供的各种光谱库，针对五种典型地物：雪、植被、水体、土壤、矿物岩石，通过绘制地物的反射光谱特性曲线，说明典型地物的反射光谱特性，并分别比较属于同一大类但处于在不同状态下的地物反射光谱特性。二、 实验数据与原理数据：ENVI自带的波谱库。原理：入射到物体表面的电磁波与物体发生三中作用：反射、吸收和投射。不同地物的反射、投射和吸收能力不同，既地物的波普特性。遥感传感器能够记录地物本身发射的电磁波信息和地物反射太阳光的电磁波信息。三、 实验过程与结果启动ENVI软件，拷贝实验所需原始数据并修改文件名。在主菜单中打开Spe
2、ctral-Spectral Libraries-Spectral Library Viewer；打开Spectral Library Input File对话框，点击open-new file;打开原始文件夹，然后分别选择雪、植被、水体、土壤、矿物岩石的波普图文件，得到五种物质在不同状态下的波普曲线图图1为雪在五种不同状态下的波普图，观察图可得，不同状态的雪在波长0.5微米附近有个波峰，随着波长增加反射率逐渐降低，在可见光波段基本上是非选择性吸收体，既高反射体，但在近红外波段吸收很强。 图2为三种植物的波普反射曲线，三种植物分别为confier(针叶树)、decidous(落叶植物)、gra
3、ss(草)。观察图可看出在可见光波段0.55微米（绿光）附近的反射率较低，10%20%左右，两侧0.45微米（蓝光）和0.67微米（红光）则有两个吸收带。在近红外波段1.21.5微米件有一个反射的陡坡，至1.6微米附近有一峰值在23微米件吸收率大增，反射率大大下降图（3）绘出了纯净水、海水的波普曲线，观察图可得出水在波长1微米左右的反射率较高，近红外波段的反射率很低。 四、 实验体会不同的地物一般有不同的反射率，根据此原理可以区分不同的地物。实验二：遥感图像目视解译一、 实验目的利用徐州四环范围的Landsat-5 TM遥感影像进行目视解译。二、 实验数据与原理利用ENVI的折线注记功能，获取
4、徐州三环路的矢量图；利用ENVI的文本注记等功能，标识出遥感影像中的山、河、湖等典型地物。三、实验过程与结果File Open Image File，选择2000年的徐州TM影像图。选择overlayAnnotation ，利用注记功能在地物上添加文字和符号。在经典地物上选择合适的符号，for example Rectangle ，and put it in a right place .操作方法：以云龙湖为例，在右窗口中选择symbol，选择三角形，鼠标左键在主窗口中定位，然后调整大小。Take Yunlong as example ,select symbol and choose it
5、,click the left mouse and adjust it to a proper place ,then click the right mouse to ensure the move.公路注记：重新打开2000年的影像图，overlayAnnotation-object-polygon，选择合适的线条和color，选择windows为scroll，在scroll中进行操作。The final result:四、实验体会 添加注记比较有技巧，需掌握。实验三 遥感影像匹配一、实验目的 在ENVI平台下，以2001年徐州城区遥感图像作为基准图像，对1994年徐 州城区遥感图像进行配
6、准。二、 实验数据与原理徐州2000年与1994年遥感影像图三、 实验过程与结果打开并显示影像文件：1、 File- Open Image File 选择徐州市1994年和2000年的影像文件2、 在Available bands List中选择2000年影像作为#1，94年得影像作为#2，分别加载开始进行影像配准并加载地面控制点 1. 从 ENVI 主菜单栏中，选择 Map Registration Select GCPs: Image to Image。 2. 在 Image to Image Registration对话框中，点击并选择 Display #1 （2000年影像），作为Ba
7、se Image。点击 Display #2（1994年 影像），作为 Warp Image。然后选择同名地物点，并选择添加选择GCP时，最好选择道路的交叉处或者明显的地方，均匀选点并选择8个左右的点。并选择保存注记。最终选择了九个地物点，结果如下配准到此结束下面进行影像校正(仅对显示波段进行校正)1. 从 Ground Control Points Selection 对话框中，选择 Options Warp Displayed Band。 2. 在 Registration Parameters 对话框中的 Warp Method按钮菜单中，选择 RST。在 Resampling 的按钮菜
8、单中选择 Nearest Neighbor 重采样法。输入文件名 test3-ljs-01.wrp，点击 OK。 4. 重复步骤 1和步骤 2，还是使用 RST 校正法，但是要相应地选择 Bilinear 和 Cubic Convolution重采样法。 5. 将结果分别输出到 test-ljs-02.wrp 和 test-ljs-03.wrp 文件中。 6. 再一次重复步骤 1 和步骤 2，这一次选择一次多项式 Polynomial 校正法，并使用 Cubic Convolution重采样法。然后再选择 Delaunay 三角网的 Triangulation 校正法，相应地使用 Cubic
9、Convolution重采样法。 7. 将结果分别输出到 test-ljs-04.wrp 和 test-ljs-05.wrp 文件中。图像编号校正重采样Test3-ljs-01RSTNearest NeighbourTest3-ljs-02RSTBilinearTest3-ljs-03RSTCubic ConvolutionTest3-ljs-04PolynomialCubic ConvolutionTest3-ljs-05DelaunayCubic Convolution比较校正结果 使用动态链接来比较校正结果： 1. 在可用波段列表中，点击原始的 TM 波段影像名 bldr_tm.img，
10、然后从菜单栏中，选择 File Close Selected File。 2. 在随后出现的 ENVI 警告对话框中，点击 Yes关闭相应的影像文件。 3. 在可用波段列表中，选择test-ljs-01.WRP 文件。在 Display #下拉式按钮中选择 New Display，点击 Load Band将该文件加载到一个新的显示窗口中。 4. 在主影像窗口中，点击鼠标右键，选择 Tools Link Link Displays。 5. 在 Link Displays 对话框中，点击 OK，把 2000年 影像和校正后的 1994年影像链接起来。 6. 在主影像显示窗口中，点击鼠标左键，使用动
11、态链接功能，对 2000年 影像和校正后的 1994年影像进行比较。7. 将 test-ljs-02.wrp 和test-ljs-03.wrp 影像加载到新的显示窗口中，使用影像动态链接功能，比较采用三种不同的重采样法（最临近法、双线性内插法和三次卷积法）所产生的效果。 结果发现在使用最近邻法重采样的影像中的锯齿状像素，而使用双线性内插法重采样的影像看起来更加平滑，使用三次卷积法重采样的影像是最好的结果，不但有平滑效果，而且保持了影像的细节特征。 8. 在相应的主影像窗口中，选择 File Cancel，关闭test-ljs-01.wrp（RST 校正，最近邻法重采样）和 test-ljs-0
12、2.wrp（RST 校正，双线性内插法重采样）影像的显示窗口。 9. 将 test-ljs-04.wrp 和 test-ljs-05.wrp 影像加载到新的显示窗口中，使用影像动态链接功能，同test-ljs-03.wrp 影像（RST 校正）进行比较。 输出结果：Test3-ljs-01（RST Nearest Neighbour）与2000年进行比较Test3-ljs-02 (RST Bilinear )与2000年进行比较Test3-ljs-03(RST Cubic Convolution )与2000年影像进行比较Test3-ljs-04（Polynomial Cubic Convol
13、ution）与2000年影像进行比较Test3-ljs-05（Delaunay Cubic Convolution）与2000年影像进行比较四、 实验体会 配准时，需选择一定的地物点，如道路的交叉处，河流房屋等边缘位置。调节位置，可提高配准精度。实验四 遥感影像增强一、 实验目的 学习ENVI软件的基本操作，能够将图像进行相应变换和增强处理，在此操作中加深对理论知识的理解和掌握，完成（1）空间域滤波增强（2）频率域滤波增强（3）彩色变换，包括真彩色增强、伪彩色增强、假彩色增强（4）植被指数（5）多光谱变换二、 数据介绍实验数据为徐州2000年TM数据三、实验过程（一）空间域增强 点运算1.线性
14、对比度拉伸Linear Contrast Stretch）线性变换所用的变换函数是线性的或分段线性的，是将像元值的变动范围按线性关系扩展到指定范围，变换函数y=a*x+b.目的是为了改善图像的对比度，改变图像像元的灰度值。线性对比度拉伸是系统默认的交互式拉伸。线性拉伸的最小和最大值分别设置为 0 和 255，两者之间的所有其它值设置为中间的线性输出值具体做法如下所示：Enhance-interactive stretching, 从 Interactive Contrast Stretching 对话框内，选择 Stretch_Type Linear Contrast Stretch，要限定最
15、小和最大输入值，点击 “Apply” ，把拉伸应用于显示的数据。如图图像变化如下：2分段线性对比度拉伸（Piecewise Linear Contrast Stretch） 分段线性变换就是在一些灰度段拉伸，另一些灰度段压缩分段线性对比度拉伸可以通过使用鼠标在输入直方图中放置几个点进行交互地限定。当在点之间提供线性拉伸时，线段在点处连接起来。具体做法如下所示：选择 Stretch_Type Piecewise Linear.，要限定最小和最大输入值，点击 “Apply” ，把拉伸应用于显示的数据。如图3高斯对比度拉伸（Gaussian Contrast Stretch）系统默认的 Gaussi
16、an 拉伸是围绕DN平均值127的三个标准差的数据分布（centered at a mean DN of 127 with the data distributed over a range of 3 standard deviations）。选择 Stretch_Type Gaussian，输入拉伸的最小和最大值，点击 “Apply” ，把拉伸应用于显示的数据。4直方图均衡化对比度拉伸（Histogram Equalization Contrast Stretch）直方图均衡化的随机分布的图像直方图修改成均匀的直方图，又叫拉 平扩展，实质是对图像进行分线性拉伸，重新分配图像像元值，使一定灰度
17、 范围内的像元值数量大致相等。要自动缩放数据以与每个直方图中bin的DN数均衡：选择 Stretch_Type Eualization.，点击 “Apply” ，输入拉伸的最小和最大值，拉伸应用于显示的数据。 平方根对比度拉伸（Square Root Contrast Stretch和任意对比度拉伸和直方图匹配（Arbitrary Contrast Stretching and Histogram Matching）操作过程同上邻域运算1. 图像卷积运算卷积运算是在空间域上对图像进行领域检测的运算，采用模板，实际上是一个M*N的小图像。Functions Interactive Analysi
18、s Spatial Pixel Editor. 要为一个像元输入一个新值，点击表格中的像元值，使用空格键删除当前值，输入所需要的值，然后按回车。 要用 4 个或 8 个周围像元的平均值来替代一个像元的值，点击表格中的像元值，然后选择 Options Replace cell with 4 pixel average or Replace cell with 8 pixel average.2. 锐化（Sharpen）锐化也称高通滤波，主要是增强图像中的高频成分，突出图像的边缘信息，提高图像的细节反差，所以叫边缘增强。锐化滤波对图像显示窗口内（image、scroll 和 zoom）的数据执行高
19、通卷积。可利用三种锐化滤波器类型，每个后面添加一个不同的数据值。“Sharpen” 附近方括号内的数字是核的中心值。因此，括号中数值较高的滤波在滤波以后，有较多的原始数据保留并附加到滤波图象上。Enhancements Filter Sharpen，结果见下图3. 平滑（Smooth）平滑又叫低通滤波，目的在于消除图像中的各种干扰噪声，使图像中高频成分消退，平滑掉图像的细节，使反差降低，保存低频成分。可使用两种平滑滤波器。“Smooth 33” 滤波器使用一个大小为 33的核，“Smooth 55” 滤波器使用一个大小为 55的核。核越大，滤波后越平滑。Functions Enhancemen
20、ts Filter Smooth结果见图4. 中值（Median）中值滤波是将领域中的像素按灰度级排序，取其中间值为输出像素。可使用两种中值滤波器，核大小分别为 33 和 55。中值滤波器用核的中值来代替中央的像元值。这些滤波器有助于减少盐点和胡椒粉类型的噪音或斑点。Functions Display Enhancements Filter Median见图（二）频率域增强1. 傅立叶变换傅立叶分析是一种将图像分成空间上各种频率成分的数学方法。实际上，快速的傅立叶变换被原来将数据变换成一个复杂的强调频率分布的图像。ENVI 中 FFT 滤波（从 Filters 下拉菜单中选择）包括图像正向的
21、FFT、频率滤波器的交互式建立、滤波器的应用，以及 FFT 向原始数据空间的逆变换。当前，FFT 处理没有用到 ENVI tiling 程序，因此能被处理的图像大小受到系统可利用内存的限制。FFT 图像是“复数”数据类型，它占用了类似大小的字节图像的 8 倍内存。2. 频率域平滑去除噪声，改善图像质量，常采用低通滤波器来抑制高频部分，然后在进行傅立叶逆变换的 滤波图像，从而可达到平滑图像的目的。常用的低通滤波器有四种：理想低通滤波器、Butterworth低通滤波器、指数低通滤波器、梯形低通滤波器。3. 频率域锐化突出图像的边缘和轮廓，采用高通滤波器来阻止或削弱低频成分，让高频成分通过，然后再
22、进行傅里叶的到边缘锐化的图像。常用的高通滤波器有四种：理想高通滤波器、Butterworth高通滤波器、指数高通滤波器、梯形高通滤波器。具体做法如下所示： 选择Filters FFT Filtering Forward FFT，现Forward FFT Input File对话框时，选择要被处理的数据，需要时，可以用子集。 文件选择程序与卷积滤波中用到的一样。点击“OK”，现Forward FFT Parameters对话框，选用“File” 或 “Memory”输出。选择Filters FFT Filtering Filter Definition，选择好以后，如果显示出一幅FFT图像，应迅
23、速选择包含FFT图像的显示数（号），或指出滤波器不与一次特殊显示相关。从Filter_Type下拉菜单下方的滤波列表里，选择一种滤波类型，并设置需要的参数。用“Number of Border Pixels”参数控制用于taper滤波器（平滑滤波边缘）的像元数。零值代表没有平滑。 选择输出到“File” or “Memory”，若需要，选择一个输出文件名。 点击 “Apply” 建立 FFT 滤波器。选择Filters FFT Filtering Inverse FFT。 出现Inverse FFT Input File对话框时，选择要处理的正向的FFT数据。出现Inverse FFT Fil
24、ter File对话框，选择应用的滤波图像。 滤波图像必须是以前用上面描述的滤波定义程序生成的。 ？（三）彩色增强1. 伪彩色增强密度分割是一种单波段图像彩色变换的方法，它是把单波段的黑白遥感图像按照亮度来分层，对于每一层赋予不同的颜色，使之成为一幅彩色图像，其中每一层包含的亮度值范围可以不同。如果分层方案与地物光谱差异对应的很好，我们通过这种方法就可以区分出各种地物的类别。具体做法如下所示：Tools Color Mapping Density Slice或Overlay Density Slice所要进行灰度值范围定义的波段Density slice 框中对已经划分好的灰度值范围区间属性进
25、行编辑，删除，变颜色。只需要你点击 Edit Range 、Delete Range、Clear Range 选项框进行操作。点击Apply 应用 变化前的灰度图像 假彩色变换之后的灰度图像2. 假彩色合成根据加色法的合成原理，选择遥感影像的某三个波段，分别赋予红，绿，蓝三种原色，就可以做成彩色合成的图像，这种图像处理方法就叫做彩色合成，它能够有效的突出更多的地物信息。假彩色合成是彩色合成中的最常用的一种方法。具体做法如下所示：首先打开一幅TM图像17个波段，FileOpen Image File:点击 Available Bands List 对话框中RGB Color
26、单选按钮。在显示的波段列表中分别点击band 3 做R（红）变换、band 2 做G（绿）变换、band1做B（蓝）变换。点击 Load RGB 显示合成后的彩色图像。灰度图像 彩色图像3. 彩色变换在图像处理中通常应用的有两种彩色坐标系统。一种是由红（R）、绿（G）、蓝（B）三原色构成的彩色空间或坐标系统，这是我们之前所讲的彩色变化和增强处理的基础；另一种是彩色空间是由色调（H）、饱和度（S）和明度（I）三个变量构成的。彩色系统变换主要是指这两种坐标系统之间的变换。具体操作如下RGB 向HSV 的转换：在主菜单中点击Transforms Color Transforms
27、H：S:四、植被指数五、多光谱变换1、主成分变换主成分变换是多元统计分析中的一个重要变换，是一种很重要的特征提取方法，减少计算量，可以用作分类图像的特征图像。也可以单独应用作为一种重要的信息提取方法。主成分与原始波段间的定量关系，主要是系数矩阵的解求，主成分的原理时已讲到系数矩阵T 是原始波段协方差矩阵的特征向量矩阵的一个转置。Transforms- Principal Components-Forward PC Rotation-Compute New Statistics and Rotate出现如下Forward PCRotation Parameters对话框：该对话框参数设置如下：“
28、Stats X/Y Resize Factor” 文本框键入小于 1 的调整系数，对计算统计值的数据进行二次抽样。键入一个小于 1 的调整系数，以提高统计计算的速度。例如，在统计计算时，用一个 0.1 的调整系数将只用到十分之一的像元。点击按钮，选择基本“Covariance Matrix”或“Correlation Matrix”。计算主成分时，有代表性地要用到协方差矩阵。当波段之间数据范围差异较大时，要用到相关系数矩阵，并且需要标准化。选用 “File” 或 “Memory” 输出。在该对话框中选择 EditData Values，查看各成分的特征值的大小 主成分图像的6个波段如下图 缨帽
29、变换穗帽变换（又称KT 变换）是一种特殊的主成分分析，和主成分分析不同的是其转换系73数是固定的，因此它独立于单个图像，不同图像产生的土壤亮度和绿度可以互相转化比较。随着植被的生长，在绿度图像上信息增强，土壤亮度上的信息减弱，当植被成熟和逐渐凋落时，其在绿色度图像特征减少，在黄色度上的信息增强。这样的解释可以应用于不同的区域上的不同植被和作物，但缨帽变换无法包含一些不是绿色的植被和不同的土壤类型的信息。选择Transforms Tassled Cap,出现Tasseled Cap Transformation Input File 对话框，选择输入的文件，用ENVI标准的选择程序输入文件。
