一种面向扎龙湿地遥感图像密度峰值聚类方法技术

本发明专利技术属于图像处理领域，公开了一种面向扎龙湿地遥感图像密度峰值聚类方法包含如下步骤：(1)对得到的湿地的高光谱图像进行筛选，人为选择多幅高光谱图像；(2)用中值滤波器和主成分分析法对多幅高光谱图像进行预处理得到PC1；(3)对PC1进行尺度变换；(4)进行PC1的像素点间的相似度计算；(5)用密度峰值聚类算法进行聚类处理；(6)将局部密度矩阵以及距离矩阵和改进的粒子群演化算法结合，找出最适合作为聚类中心的所有像素点；(7)以作为聚类中心的像素点为中心，根据局部密度矩阵以及截断距离对其余的像素点进行划分；(8)输出处理完毕的遥感图像。本发明专利技术降低了算法的时间复杂度，解决了聚类结果不稳定的问题，更具有普适性。

A density peak clustering method for Zhalong wetland remote sensing image

The invention belongs to the field of image processing, and discloses a clustering method for the peak density of remote sensing image of Zhalong wetland, which comprises the following steps: (1) screening the obtained hyperspectral image of wetland, selecting multiple hyperspectral images artificially; (2) preprocessing multiple hyperspectral images with median filter and principal component analysis method. Get PC1; (3) Scale transform PC1; (4) Compute the similarity between PC1 pixels; (5) Cluster processing with density peak clustering algorithm; (6) Combine local density matrix and distance matrix with improved particle swarm optimization algorithm to find all the pixels most suitable for clustering center; (7) Cluster as clustering center. The other pixels are divided according to the local density matrix and truncation distance; (8) the processed remote sensing image is output. The invention reduces the time complexity of the algorithm, solves the problem of unstable clustering results, and has more universality.

【技术实现步骤摘要】
一种面向扎龙湿地遥感图像密度峰值聚类方法
本专利技术属于图像处理领域，尤其涉及一种面向扎龙湿地遥感图像密度峰值聚类方法。
技术介绍
遥感技术诞生于上个世纪六十年代，是一门以地理学和物理学为基础，结合了计算机，信号处理，空间监测等科学的综合学科。通过卫星或者航空器上的多个传感器记录下各种地物电磁波的图片，用于环境监测，资源勘探，军事侦察等领域。随着成像光谱技术的发展与成熟，高光谱图像成为了遥感图像的主要成像方式。高光谱图像在获得地表图像信息的同时，也获得其光谱信息，第一次真正做到了光谱与图像的结合。与多光谱遥感影像相比，不仅在信息丰富程度方面有了极大的提高，在处理技术上，对该类光谱数据进行更为合理、有效的分析处理提供了可能。在提供了更为丰富的光谱信息的同时，高光谱图像也造成了“同物异谱，同谱异物”等现象，为我们后期的遥感图像处理带来了难度。遥感图像聚类作为主要对高光谱图像处理进而得到更详细的地物信息的有效手段之一，一直以来得到了学者广泛地关注。尤其是在图像分辨率与以往相比有了较大提升的情况下，一个优秀的图像处理方法的出现对于我们进一步获取地物的详细信息显得更为重要。与分类问题不同，聚类是一种无监督的学习。对于遥感图像分类，需要预先提取部分目标图像作为训练样本，在训练集中每个类别都具有属于自己的标签，在训练完成得到一个成熟的分类模型之后再对目标图像进行处理。虽然提升了精确度，但是时间复杂度也提高很多，并且由于训练集和先验知识的影响，导致算法不稳定。而聚类算法由于时间复杂度相对较低，而且没有使用训练集预先训练模型排除了人为干扰等因素，能有效减少对含量较少的地物光谱的丢失，适用于高光谱图像中不同波段间的多幅图像。遥感图像聚类的方法主要包括以下几种：1.K-means遥感图像聚类算法：K-means算法作为一种无监督的聚类算法具有简单，快速，高效等优点。但是必须要预先指定聚类数目，同时要确定初始聚类中心点。中心点确定后，通过计算各个像素点与聚类中心得距离以所有像素点距离的标准最低为准则来更新聚类中心。通过不断地迭代过程，最终完成所有像素点的类别划分。在聚类的过程中，类别数目是预先设定好的，因此会导致最终聚类结果不准确，出现高光谱图像像素点错分漏分的情况。其次，中心点的确定完全是人为设定，选取不同的中心点会造成最终结果出现较大偏差，导致算法不稳定为后续对图像信息的进一步提取带来阻碍。当遥感图像特征差异较大时，该算法必须先分为若干个类，在经过整理合并才能得到最终的结果。当图像内的波谱特征比较复杂时，该算法在聚类中心和数目的选择上需要大量的判读工作才能确定。湿地地貌特征复杂，具有城镇，草地，湖泊，农田等十几种地貌特征，此时K-means算法显得不再适用。2.模糊C均值遥感图像聚类算法：模糊C均值算法是一种结合了无监督聚类和模糊集概念的图像处理技术，在算机视觉，图像处理等领域具有广泛应用。该算法对于像素点的划分与要求每个像素点只能属于一个类簇中心的硬划分聚类算法不同，通过对模糊集的应用，构建像素点和聚类中心的隶属度矩阵来进行聚类。在某些数据点类别不清晰的数据集上可以取得比普通聚类算法更优秀的结果。但是，该算法同样也存在以下缺陷。它易受初始聚类中心和初始隶属度矩阵的影响，可能收敛到局部极小值，影响着分割的效果；在大量数据集的情况下，比如说遥感图像，模糊C-均值的反复迭代就会造成算法耗时过长。而随着传感器的不断改进，现代遥感图像的像素和分辨率也越来越高，模糊C均值算法在聚类过程中需要反复迭代运算，更加延缓了该算法的收敛速度。从以上分析来看，目前这两种现有的聚类算法均不适用于湿地的特殊的地理环境分析。因此，针对湿地的遥感图像聚类算法的研究显得十分有意义。
技术实现思路
本专利技术的目的在于公开适用于湿地的聚类结果稳定的普适性好的一种面向扎龙湿地遥感图像密度峰值聚类方法。本专利技术的目的是这样实现的：一种面向扎龙湿地遥感图像密度峰值聚类方法，包含如下步骤：步骤(1)：首先对得到的湿地的高光谱图像进行筛选，人为选择出地物清晰，地貌特征明显，分辨率较高并且拍摄时没有或者只有较少云层遮挡的多幅高光谱图像；步骤(2)：用中值滤波器和主成分分析法对多幅高光谱图像进行预处理得到PC1：采用中值滤波器处理多幅高光谱图像；对中值滤波器处理后的多幅高光谱图像中的每一张采用主成分分析法进行降维处理，在尽量保存原始信息的同时去除与分析指标无关的以及重复的变量。选取样本方差作为综合指标函数指导样本选取，从降维处理后的多幅高光谱图像中选取方差最大的作为PC1；步骤(3)：对PC1进行尺度变换：采用弱采样的方式进行尺度变换，将待处理的子图尺寸由AxB变换为A1xB1，其中A1＝A/ksize且B1＝B/ksize；ksize为尺度变换指数；步骤(4)：进行PC1的像素点间的相似度计算：以尺度变换后的PC1的左上角为原点，水平线为X轴，竖直线为Y轴建立直角坐标系；假设经过步骤(3)处理后的子图包含M行N列个像素点，建立含有n＝M*N个数据的像素数据集；单个数据包含该像素点在RGB彩色空间下R，G，B分量方向上的像素值，以及位置信息等特征信息；假设像素点p1，p2属于该像素数据集，则像素点p1,p2的相似度dsim(p1,p2)：dsim(p1,p2)＝βdrgb+(1-β)dxy；上式中，r1表示像素点p1在RGB空间下R分量方向上的像素值,g1表示像素点p1在RGB空间下G分量方向上的像素值,b1表示像素点p1在RGB空间下B分量方向上的像素值；r2表示像素点p2在RGB空间下R分量方向上的像素值,g2表示像素点p2在RGB空间下G分量方向上的像素值,b2表示像素点p2在RGB空间下B分量方向上的像素值；(x1,y1)为像素点p1的坐标，(x2,y2)为像素点p2的坐标；β为权重因子，β∈(0,1)；步骤(5)：用密度峰值聚类算法进行聚类处理：计算各个点的局部密度，创建局部密度矩阵{ρi}：图像的像素集指标集为Is＝{1,2，...M*N}，dij表示像素点i和像素点j的相似度，dc为人为设定的截断距离；计算各个像素点间与比该像素点具有更高局部密度的像素点的距离值，创建距离矩阵设表示的一个降序排列，即：ρq1≥ρq2≥...≥ρqn；定义δqi：当i≥2,当i＝1时，步骤(6)：将局部密度矩阵以及距离矩阵和改进的粒子群演化算法结合，找出最适合作为聚类中心的所有像素点：在解空间初始化整个种群，将每个像素点均当做粒子群算法中的个体；赋予每个粒子初始速度以及位置，同时在搜索过程中记录每个粒子自身到达的最优位置，以及整个种群的最优个体；目标函数设定为：找到局部密度与距离均得到较大值的像素点作为聚类中心，即使目标函数取得最小值的解；步骤(7)：以作为聚类中心的像素点为中心，根据局部密度矩阵以及截断距离对其余的像素点进行划分；步骤(8)：输出处理完毕的遥感图像。本专利技术的有益效果为：本专利技术在尽量保留图像信息的基础上降低了整个算法的时间复杂度，解决了现有聚类算法选择聚类中心受先验知识影响导致聚类结果不稳定的问题；结合粒子群演化算法，将该问题转化为寻找最优解；密度峰值算法更具有普适性；在时间复杂度允许的情况下，本专利技术很好地解决了高光谱图像“同物异谱，同谱异物”的问题。附图说本文档来自技高网...

【技术保护点】
1.一种面向扎龙湿地遥感图像密度峰值聚类方法，其特征在于：包含如下步骤：步骤(1)：首先对得到的湿地的高光谱图像进行筛选，人为选择出地物清晰，地貌特征明显，分辨率较高并且拍摄时没有或者只有较少云层遮挡的多幅高光谱图像；步骤(2)：用中值滤波器和主成分分析法对多幅高光谱图像进行预处理得到PC1；步骤(3)：对PC1进行尺度变换；步骤(4)：进行PC1的像素点间的相似度计算；步骤(5)：用密度峰值聚类算法进行聚类处理；步骤(6)：将局部密度矩阵以及距离矩阵和改进的粒子群演化算法结合，找出最适合作为聚类中心的所有像素点；步骤(7)：以作为聚类中心的像素点为中心，根据局部密度矩阵以及截断距离对其余的像素点进行划分；步骤(8)：输出处理完毕的遥感图像。

【技术特征摘要】
2018.04.25 CN 201810375568X1.一种面向扎龙湿地遥感图像密度峰值聚类方法，其特征在于：包含如下步骤：步骤(1)：首先对得到的湿地的高光谱图像进行筛选，人为选择出地物清晰，地貌特征明显，分辨率较高并且拍摄时没有或者只有较少云层遮挡的多幅高光谱图像；步骤(2)：用中值滤波器和主成分分析法对多幅高光谱图像进行预处理得到PC1；步骤(3)：对PC1进行尺度变换；步骤(4)：进行PC1的像素点间的相似度计算；步骤(5)：用密度峰值聚类算法进行聚类处理；步骤(6)：将局部密度矩阵以及距离矩阵和改进的粒子群演化算法结合，找出最适合作为聚类中心的所有像素点；步骤(7)：以作为聚类中心的像素点为中心，根据局部密度矩阵以及截断距离对其余的像素点进行划分；步骤(8)：输出处理完毕的遥感图像。2.根据权利要求1所述的一种面向扎龙湿地遥感图像密度峰值聚类方法，其特征在于：所述的步骤(2)具体为：采用中值滤波器处理多幅高光谱图像；对中值滤波器处理后的多幅高光谱图像中的每一张采用主成分分析法进行降维处理，在尽量保存原始信息的同时去除与分析指标无关的以及重复的变量；选取样本方差作为综合指标函数指导样本选取，从降维处理后的多幅高光谱图像中选取方差最大的作为PC1。3.根据权利要求1所述的一种面向扎龙湿地遥感图像密度峰值聚类方法，其特征在于：所述的步骤(3)具体为：采用弱采样的方式进行尺度变换，将待处理的子图尺寸由AxB变换为A1xB1，其中A1＝A/ksize且B1＝B/ksize；ksize为尺度变换指数。4.根据权利要求1所述的一种面向扎龙湿地遥感图像密度峰值聚类方法，其特征在于：所述的步骤(4)具体为：以尺度变换后的PC1的左上角为原点，水平线为X轴，竖直线为Y轴建立直角坐标...

【专利技术属性】
技术研发人员：董红斌，余陈，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：黑龙江,23