遥感影像监督分类与非监督分类及相关代码实现

遥感影像监督分类与非监督分类及相关代码实现

ID:14809406

大小:227.15 KB

页数:43页

时间:2018-07-30

遥感影像监督分类与非监督分类及相关代码实现_第1页
遥感影像监督分类与非监督分类及相关代码实现_第2页
遥感影像监督分类与非监督分类及相关代码实现_第3页
遥感影像监督分类与非监督分类及相关代码实现_第4页
遥感影像监督分类与非监督分类及相关代码实现_第5页
资源描述:

《遥感影像监督分类与非监督分类及相关代码实现》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、遥感影像监督分类与非监督分类摘要:遥感影像的分类方法按照是否有先验类别可分为监督分类和非监督分类,这两种分类方法有着本质的区别但也有存在一定的联系。本文从分类原理和分类方法等不同角度分别介绍了监督分类和非监督分类方法,并对两种方法的分类结果进行了对比和分析。关键词:遥感;监督分类;非监督分类;ISODATA算法;贝叶斯分类算法。1.数据来源本文使用的数据是华盛顿广场上空由卫星拍摄的高光谱遥感影像。该幅影像使用的传感器系统覆盖0.4到2.4m的可见光到近红外的210个波段。由于0.9-1.4米光谱对应的区域上空大气透光性很差,因此将这个区域内的波段从影像中删除,最后

2、得到191个波段。该数据集有1208个扫描行,每行307个像元,容量近似150MB。为了清晰地反映影像中地物的特征,本文从191个波段中选择了3个波段,分别是120、140和160。下面两幅图分别是全波段影像和三波段影像: 图1全波段影像图2三波段影像1.遥感影像的监督分类2.1监督分类的原理监督分类(supervisedclassification)又称训练场地法,是以建立统计识别函数为理论基础,依据典型样本训练方法进行分类的技术。即根据已知训练区提供的样本,通过选择特征参数,求出特征参数作为决策规则,建立判别函数以对各待分类影像进行的图像分类,是模式识别的一种

3、方法。要求训练区域具有典型性和代表性。判别准则若满足分类精度要求,则此准则成立;反之,需重新建立分类的决策规则,直至满足分类精度要求为止。常用算法有:最大释然分类法、最小距离分类法、马氏距离分类法、平行六面体分类法、K-NN分类法。本文选用最大释然分类法对遥感影像进行监督分类。最大似然判别法.也称为贝叶斯(Bayes)分类,是基于图像统计的监督分类法,也是典型的和应用最广的监督分类方法.它建立在Bayes准则的基础上,偏重于集群分布的统计特性,分类原理是假定训练样本数据在光谱空间的分布是服从高斯正态分布规律的,做出样本的概率密度等值线,确定分类,然后通过计算标本(

4、像元)属于各组(类)的概率,将标本归属于概率最大的一组.用最大似然法分类,具体分为三步:首先确定各类的训练样本,再根据训练样本计算各类的统计特征值,建立分类判别函数,最后逐点扫描影像各像元,将像元特征向量代入判别函数,求出其属于各类的概率,将待判断像元归属于最大判别函数值的一组。2.2最大释然分类算法的调试#include#include#include#include#includeconstfloatPI=3.1415;floattrain1[60][3],

5、train2[60][3],train3[60][3],train4[60][3];//存各类训练样本floatm1[3],m2[3],m3[3],m4[3],c1[3][3],c2[3][3],c3[3][3],c4[3][3];//各类均值向量及协方差矩阵floatp=0.25;//先验概率floattest[240][11];//检验样本各列分别存放:分类前类别,行号,列号,一波段灰度,二波段灰度,三波段灰度,分类后类别floata1,a2,a3,a4;//公式中常数项;floatqq,qq1,qq2,qq3,qq4;//分类精度floatk;//kappa

6、值intq[4][4];//混淆矩阵floata[4][4];//精度指标矩阵各列分别为生产者精度漏分误差//用户精度错分误差intmain(){voidgetdata1(char*filename);voidgetdata2(char*filename);voidcompute(floattrain[60][3],floatm[3],floatc[3][3],float*a);floatidentify(floatm[3],floatc[3][3],floata,floatband1,floatband2,floatband3);intclassify(floa

7、tp1,floatp2,floatp3,floatp4);voidtesting();voidaccuracy(floatp[4][4]);voidoutput(char*filename);voidoutput1(char*filename);getdata1("05training.txt");//存放的是训练样本三个波段的灰度值getdata2("05testing.txt");//检验样本行号列号类别代码三个波段的灰度值compute(train1,m1,c1,&a1);//计算各类均值向量及协方差矩阵compute(train2,m2,c2,&a2);c

8、omput

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。