对遥感图像的分类处理综合应用

一、课程设计的研究现状、目的与意义

随着遥感技术的发展与成熟，遥感技术在各个领域的应用越来越广泛。遥感技术正朝着定量化、智能化、动态化、 络化、实用化等方向发展，随着遥感数据获取手段的加强，需要处理的遥感信息急剧增加。遥感技术正在进入一个能够快速准确地提供多种对地观测海量数据以及应用研究的新阶段，有达到一个新的高潮。
遥感图像分类在遥感图像处理中就是一个很重要的方向，根据遥感图像自身的特点使计算机对遥感图像按照一定的意义进行自动的分类处理，将使人们在面对海量的遥感数据时更方便的进行筛选以及分析应用。因而，不断有遥感科学领域研究者在图像分类的算法的借鉴，改进和创新中探索。
加深理解和巩固遥感的相关知识；培养运用ENVI软件独立分析问题、解决问题的实际工作能力；培养良好的工作习惯和科学素养，为今后工作打好基础。
更大加深对遥感相关知识的认知，以及对遥感图像处理软件的操作熟练度的掌握。无论以后是考研还是从事相关一类的工作,都为此打下坚厚的基础。

二、ENVI系统介绍

ENVI是遥感科学家用交互式数据语言开发的一套功能强大的遥感图像处理软件，能够有效的从遥感影像中提取各种目标信息。
自 ENVI5.0 版本开始，ENVI 采用了全新的软件界面，从整体上增强了用户体验，ENVI5.3 延续了ENVI5 的界面风格，对图标做了更现代化的设计。启动 ENVI5.3，如图所示，包括菜单项、工具栏、图层管理、工具箱、状态栏几个部分组成。

3.2 辐射定标

将传感器记录的电压或数字量化值(DN)转换成绝对辐射亮度值(辐射率)的过程，或者转换与地表(表观)反射率、表面(表观)温度等物理量有关的相对值的处理过程。消除传感器本身的误差，确定传感器入口的准确辐射值。对于一般的线性传感器，定标通过一个线性关系式完成数字量化与辐射亮度值的转换：
L=Gain*DN+Offset
在ENVI中打开实验数据利用ENVI自带的定标工具进行定标，获取辐射亮度或反射率。
辐射定标过程中可以仅对某一局部地区进行定标，此次课设是仅对唐海县地区进行辐射定标。

由图可以看出辐射定标后遥感图像的亮度灰度值转换为绝对亮度值，显示的地物更容易操作观察了。

3.3 大气校正

将辐射亮度或者表观反射率转换为地表实际反射率，目的是消除大气散射、吸收、反射引起的误差。基本原理：首先估计大气对地表热辐射的影响，然后把这部分大气影响从卫星传感器所观测到的热辐射总量中减去，从而得到地表热辐射强度，再把这热辐射强度转化为相应的地表温度。
利用Flaash大气校正。FLAASH校正文件的路径中心点经纬度Sceneenterocation自动获取；其对应的传感器高度以及影像数据的分辨率自动读取；设置唐海县的地面高程数据：0.001km；影像生成时的飞行过境时间（可以从元文件“_MTL.txt”中找到），大气模型参数选择，根据成像时间和纬度信息选择，参考图进行选择。

由图对比可以看出0.5-0.7um波段，植被的光谱曲线在大气校正前后发生很大的变化，主要是因为消除了大气散射，主要是消除瑞利散射的原因。

可以发现消除大气散射、吸收、反射引起的误差后图像几乎没有发生变化，但在某一局部区域会有比较清晰的反应，所以有时候可以完全忽略遥感数据的大气影响。例如，对某些分类和变化检测而言，大气校正并不是必需的。

四、图像分析

4.1 不同地物的光谱曲线

任何地物都有自身的电磁辐射规律（反射，发射，吸收电磁波），根据地物反射率与波长之间的关系而绘制的曲线。选择Display>Profiles>Spectral查看各个的地物波谱曲线。
下图为水体、植被、鱼塘、耕地的光谱曲线。

植被色素吸收决定着可见光波段的光谱反射率，细胞结构决定近红外波段的光谱反射率，而水汽吸收决定了短波红外的光谱反射率特性。由于植被叶子的细胞结构的影像，除了吸收和透射的部分形成高反射率，在近红外波段的0.8-1.0微米之间有一个反射的陡坡，0.9微米附近有一个峰值，形成植被的独有特征。

耕地表面没有明显的峰值和谷值，影响土壤的光谱反射率的是：一方面由于土壤颗粒大小本身的影响，另一方面由于土壤质地影响土壤持水能力而产生的间接影响。
由以上五种不同的光谱曲线可知不同的地物类型变化很大。地物反射光谱是指地物的反射率随入射波长而变化的规律。根据地物的反射光谱所绘制的曲线成为地物反射光谱曲线，通过地物反射光谱曲线的不同辨别地物是遥感识别地物性质的基本原理。
地物的反射光谱:不同的地物在不同波段反射率存在差异(如水体、鱼塘的光谱曲线)；相同地物光谱曲线有相似性，但是也存在差异性；地物光谱特征具有事件性和空间性(不同时间与空间光谱特征不同)。

4.2 不同波段组合生成的假彩色图案

（1）5,4,3 NIR、Red、Green——标准假彩色合成（CIR）

在这种波段组合下，植被显示为红色，植被越健康红色越亮，而且还可以区分出植被的种类，这种波段组合方式非常常用，用来监测植被、农作物和湿地，如图。

（3）5,6,4 NIR、SWIR1、Red——假彩色合成，有效区分陆地和水体

这种波段组合，深浅的橙色和绿色是陆地，深浅蓝色是水。

（5）6,5,2——假彩色合成，农作物监测

这种波段组合，对监测农作物很有效，农作物显示为高亮的绿色，裸地显示为品红色，休耕地显示为很弱的墨绿色。

由上图可知:不同波段的最小值、最大值各不同。波段一有一个转折点，转折点是最大值为7738在672处取得，标准差为781.2;波段二有一个转折点，转折点是最大值为7459在665处取得，标准差为314.7;波段三有四个转折点，最大值为76035在882处取得，标准差为359.9;波段四有三个转折点,最大值为7708在882处取得，标准差为502.2;波段五有七个转折点，转折点之间起伏较大，其中最大值为8015在129处取得，标准差为1376.4;波段六有三个转折点，最大值为12722在59处取得，标准差为1344.1。波段七有四个转折点，最大值为14950在43处取得，标准差为1167.3。

4.4 主成分分析

ENVI主成分分析(PCA)是通过使用Principal Components选项生成互不相关的输出波段，去除波段多余信息，将多波段的图像信息压缩到比原波段更有效的少数几个转化波段的方法，达到隔离噪声和减少数据集的维数的方法。 由于多波段数据经常是高度相关的，主成分变换寻找一个原点在数据均值的新的坐标系统，通过坐标轴的旋转来使数据的方差达到最大。主成分分析的目的：一是进行遥感影像数据的压缩处理。对于一幅多光谱或是高光谱影像，虽然有着众多的波段信息，这些信息一方面对探测地物的细微差别有着巨大的潜能，但另一方面，这些大量的信息又存在很大的冗余，这对后期的信息处理和利用会带来一些不便。 通过主成分分析法的处理，能很好地调整信息的重新分布而又不损失信息的总量,将大部分信息集中到前几个主分量图像上，从而起到压缩数据的目的。
二是突显影像信息的作用。虽然遥感传感器本身的谱段设计所考虑的因素之一也是地物的光谱特征，这使获取的多光谱或高光谱各波段影像数据上具有了不同地物或强或弱的信息。但在利用影像数据提取某些想要的信息时，原始的影像各波段数据上所存在的信息分布得不够突显，经常会被背景信息或干扰信息所“淹没”。而利用主成分分析法，可以在分析原始影像上地物的光谱特征基础上，选择部分原始波段影像数据进行主成分分析变换，变换后信息得到调整的主分量中，就有想要的经过信息重新聚集的图像，在这个图像数据里，所要的信息被增强突显，其他信息被抑制。
下图为以大气校正的图像为基础，进行主成分分析的结果。

上述图亮度、绿度、湿度三个分量分别反映了土壤岩石、植被及土壤和植被中的水分信息。
前三个分量合成彩色图

NDVI众所周知的一种植被指数，如图为NDVI计算结果图，一般绿色植被区的范围是0.2-0.8。

4.7.2 比值植被指数（SR）

SR指数也是众所周知的一种植被指数，指在近红外波段范围绿叶的散射与红色波段范围叶绿素吸收的比值，SRI值的范围是0-30+，如图，一般绿色植被区的范围是2-8。在LAI值很高，即植被茂密时其灵敏度会降低，如图所示。

五、遥感图像分类

利用计算机通过对遥感图像中各类地物的光谱信息和空间信息进行分析，选择特征，将图像中各个像元按照某种规则或算法划分不同的类别，然后获得遥感图像中与实际地物的对应信息，从而实现图像的分类。

5.1 监督分类

监督分类是在分类前人们已对遥感影像样本区中的类别属性有了先验知识，进而可利用这些样本类别的特征作为依据建立和训练分类器，进而完成整幅影像的类型划分，将每个像元归并到相对应的一个类别中去。换句话说，监督分类就是根据地表覆盖分类体系、方案进行遥感影像的对比分析，据此建立影像分类判别规则，最后完成整景影像的分类。

5.1.1 样本选择

2）最大似然

假设每一个波段的每一类统计都呈正态分布，计算给定像元属于某一训练样本的似然度，像元最终被归并到似然度最大的一类当中，计算结果如图。

遥感影响被分为15类，相对于监督分类，ISODATA非监督分类的精度很低，分类结果很不理想，并且存在同一类型地类被划分为多块区域的情况，不适用合并内容的情况，后续处理复杂冗长。

六、分类后处理

监督分类和非监督分类方法得到的结果一般是初步结果，难于达到最终的目的。因此，需要对初步分类结果进行一些处理，才能满足需求的分类结果。

6.1小斑块去除

应用监督分类或者非监督分类的分类结果中不可避免会产生一些面积很小的图斑。无论从专题制图的角度，还是从实际应用的角度，都有必要对这些小图斑进行剔除，利用Majority/Minority分析方法进行处理。
Majority/Minority分析采用类似于卷积滤波的方法将较大类别中的虚假像元到该类中。其中，Kernel Size 为核的大小，核越大，则处理后结果越平滑；中心像元权重为在判定在变换核中哪个类别占主体地位，中心像元权重用于设定中心像元类别将被计算多少次。下图为Majority分析结果（左）和原始分类结果（右）。

八、课程设计的总结与体会

通过此次的课程设计主要考察了对ENVI软件以及遥感知识的综合运用，更是对问题的整体分析，思考思路的方式以及对透过问题看本质的磨练。感受良多，大大加强了我的动手能力和解决实际问题的能力，使我对已学的课本知识理论联系现实的灵活应用。
通过将近两周的遥感图像处理，不仅仅是对所学习知识的运用，更是对动手能力以及思考能力的双重考核。学到了很多以前不清楚甚至盲区的知识点，深深锻炼了实际中上机操作能力以及学习能力。为以后无论是考研还是工作都打下了一定的基础。