CN104598922A - 基于模糊c均值的全极化sar图像分类方法 - Google Patents
基于模糊c均值的全极化sar图像分类方法 Download PDFInfo
- Publication number
- CN104598922A CN104598922A CN201510007661.1A CN201510007661A CN104598922A CN 104598922 A CN104598922 A CN 104598922A CN 201510007661 A CN201510007661 A CN 201510007661A CN 104598922 A CN104598922 A CN 104598922A
- Authority
- CN
- China
- Prior art keywords
- fuzzy
- class
- distance
- scattering
- algorithm
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2415—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on parametric or probabilistic models, e.g. based on likelihood ratio or false acceptance rate versus a false rejection rate
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Probability & Statistics with Applications (AREA)
- Radar Systems Or Details Thereof (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于模糊C均值的全极化SAR图像分类方法,针对传统的基于Cloude-Pottier分解的极化SAR图像分类方法中由于分割阈值固定导致的地类边界不准确的问题,采用Wishart距离取代传统模糊C均值聚类中的欧氏距离,并引入一个距离因子,使得分类结果与真实地表更加接近,准确度更高。本发明所达到的有益效果:本发明中基于改进的模糊C均值聚类方法采用Wishart距离取代传统的欧氏距离,更加符合全极化SAR图像中像素点的分布规律,同时还引入了一个距离因子来增加散射机制相似的地类之间的可分度。该方法有效解决了传统的基于Cloude-Pottier分解的极化SAR图像分类中由于分割阈值固定导致的地类边界不准确的问题,因此得到的分类结果与真实地表更加接近。
Description
技术领域
本发明涉及一种基于模糊C均值的全极化SAR图像分类方法,属于极化SAR数据处理技术领域。
背景技术
近些年,随着星载极化SAR***的陆续投入使用,极化SAR技术在遥感应用中占据越来越重要的位置。由于具有四个极化通道,全极化SAR数据获取的地物信息比单极化SAR更加丰富。如何有效地分析全极化SAR数据,提取全极化SAR数据中包含的散射体的特征一直是极化研究的难点。许多极化特征参数曾经被提出,如通道强度、通道强度比、通道相位差等,但这些特征参数只对特定的应用环境适用,受具体实验和需求的限制。
Cloude-Pottier算法作为一种典型的目标极化分解方法在解决以上难题中有绝对的优势。该算法不需要知道数据的概率分布状态就可以对分类结果进行合理的解译,但是由于散射机制与真实地表之间缺乏一一对应关系,所以每种地物的分类边界比较模糊。
发明内容
为解决现有技术的不足,本发明的目的在于提供一种基于模糊C均值的全极化SAR图像分类方法,针对传统的基于Cloude-Pottier分解的极化SAR图像分类方法中由于分割阈值固定导致的地类边界不准确的问题,采用Wishart距离取代传统模糊C均值聚类中的欧氏距离,并引入一个距离因子,使得分类结果与真实地表更加接近,准确度更高。
为了实现上述目标,本发明采用如下的技术方案:
基于模糊C均值的全极化SAR图像分类方法,其特征在于:包括如下步骤:
1)对原始全极化SAR数据进行预处理,采用多视处理和Refined Lee滤波方法消除SAR图像中斑点噪声的影响;
2)对预处理后的图像进行Cloude-Pottier极化分解,得到散射熵H、散射角α、平均散射强度λ三个极化参数;
3)根据平均散射强度的值将图像分为三大类,分别为:高散射强度区域、中等散射强度区域和低散射强度区域,将每一大类利用H/α平面进一步分割,共得到24小类;
4)采用层次聚类算法将得到的24类合并为n(0<n≤24)类;
5)计算每个像素点与每一类聚类中心的距离,采用FCM算法调整每一类地物的边界,直至满足目标函数最小,迭代终止并输出分类结果;所述FCM算法中的像素点到聚类中心的距离d采用Wishart距离和一个距离因子Wij重新定义:
其中<T>为像素点的相干矩阵,Vm为每一类的平均相干矩阵,即聚类中心,μij代表数据点xi隶属于类别j的概率;所述距离因子Wij表示的含义为Wij越大,每个像素与类别中心Vm的距离越小,离聚类中心近的点变得更近,与聚类中心远的点变得更远。
前述的基于模糊C均值的全极化SAR图像分类方法,其特征在于:所述步骤2)中,散射熵H、散射角α及平均散射强度λ三个极化参数采用Cloude-Pottier分解算法对预处理后的数据进行极化分解得到:
λ=λ1pl+λ2p2+λ3p3
其中λi为相干矩阵T的特征值,且∞>λ1≥λ2≥λ3>0,pi的表达式为:
前述的基于模糊C均值的全极化SAR图像分类方法,其特征在于:所述步骤3)中,根据平均散射强度将极化数据分成三大类,每类的λ边界的确定采用取中值的方法,即根据λ值的动态范围将图像分成均匀的三部分,每部分分别对应高散射强度区域、中等散射强度区域和低散射强度区域。
前述的基于模糊C均值的全极化SAR图像分类方法,其特征在于:所述步骤4)中,利用层次聚类算法合并过分割的类别,公式如下:
其中Ti与Tj代表两个类别的类别中心,Dij为Wishart距离,表示类别i与类别j之间的分离度;根据所选试验区域地物覆盖的真实情况,将图像分为n(0<n≤24)类。
前述的基于模糊C均值的全极化SAR图像分类方法,其特征在于:所述步骤5)中模糊C均值聚类算法的目标函数为:
其中m为任意大于1的实数,C为类别数目,N为图像中像素点的个数,μij代表数据点xi隶属于类别j的概率,cj为第j类的类别中心。d为像素点到聚类中心的距离。
本发明所达到的有益效果:本发明在传统Cloude-Pottier分解的基础上,利用平均散射熵将初始图像分为三大类,避免了散射机制相似但是散射强度不同的地物之间的混淆;在模糊C均值聚类的过程中,用Wishart距离取代传统的欧氏距离,与全极化SAR图像中像素点的分布更为吻合;引入了一个距离因子,在减小相同类别像素点之间距离的同时,增大不同类别像素点的距离。该方法有效解决了传统的基于Cloude-Pottier分解的极化SAR图像分类中由于分割阈值固定导致的地类边界不准确的问题,因此得到的分类结果与真实地表更加接近,对散射机制相似的地物类型的分类结果更准确。
附图说明
图1为本发明的流程图;
图2为本发明利用Cloude-Pottier分解提取的极化参数图像;
图3为利用本发明所提分类方法、以及现有的H/α/FCM、H/α/λ/FCM方法的最终结果。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
如图1所示的基于模糊C均值的全极化SAR图像分类方法,包括以下步骤:
步骤1,对原始全极化SAR数据进行预处理,分别采用多视处理及RefinedLee滤波等消除SAR图像中斑点噪声的影响,增强图像的可读性。
为了定性和定量分析本发明方法的有效性,数据采用L波段ALOS PALSAR全极化数据,图像入射角度为23.858°,空间分辨率为9.37m×3.57m(距离向×方位向)。为了增强图像的可读性及减少相干斑噪声对实验结果的影响,首先对原始数据进行预处理。其中多视处理在方位向和距离向的比例为6:1,RefinedLee滤波器的窗口大小为3×3。
需要说明的是:本发明的方法不仅适用于该实验中选取的ALOS PALSAR数据,对其他星载及机载数据同样适用,只是在预处理的时候多视处理的比例和滤波处理的窗口大小会有不同,需要针对不同的数据源进行选择。
步骤2,对预处理后的图像进行Cloude-Pottier极化分解,得到散射熵H(图2b)、散射角α(图2c)、平均散射强度λ(图2a)等三个极化参数;
利用Cloude-Pottier分解算法对预处理后的数据进行极化分解的公式如下:
λ=λ1p1+λ2p2+λ3p3
其中λi为相干矩阵T的特征值,且∞>λ1≥λ2≥λ3>0,pi的表达式为:
步骤3,根据平均散射强度的值将图像分为三大类,每类的λ边界的确定采用取中值的方法,即根据λ值的动态范围将图像分成均匀的三部分,每部分分别对应:高散射强度区域、中等散射强度区域和低散射强度区域,在每一大类中,利用H/α平面进一步分割,共得到24小类;
步骤4,为了避免上一步骤中出现过分割,根据实验中选取试验区域的地物覆盖情况,采用层次聚类算法将以上得到的24类合并为7类。需要说明的是:聚合成多少类别是由所选区域的地表覆盖类型确定的,该参数可以根据选取的试验区域进行调整。
层次聚类算法的计算公式如下:
其中Ti与Tj代表两个类别的类别中心,Dij为Wishart距离,表示类别i与类别j之间的分离度。根据对图像覆盖区域进行野外调查及与高分辨率光学卫星图像进行对比分析,将图像分为7类。
步骤5,计算每个像素点与每一类聚类中心的距离,采用改进的FCM算法调整每一类地物的边界,直至满足目标函数最小,迭代终止并输出分类结果。
采用模糊C均值聚类算法的目标函数为:
其中m为任意大于1的实数,C为类别数目,N为图像中像素点的个数,μij代表数据点xi隶属于类别j的概率,cj为第j类的类别中心。d为像素点到聚类中心的距离。本方法采用Wishart距离和一个距离因子Wij重新定义距离d,来降低传统FCM算法中欧氏距离对极化SAR计算不准确的缺陷。
其中<T>为像素点的相干矩阵,Vm为每类的平均相干矩阵,即聚类中心,μij代表数据点xi隶属于类别j的概率。Wij越大,每个像素与类别中心Vm的距离越小,该因子使得离聚类中心近的点变得更近,与聚类中心远的点变得更远。
为了定量分析本发明全极化SAR图像分类方法的效果(图3c),采用总体精度(OA)、生产者精度(PA)和用户精度(UA)评价该发明中分类结果的精度,并与H/α/FCM算法(图3a)和H/α/λ/FCM算法(图3b)进行比较(表1)。
需要说明的是:在H/α/FCM方法中对预处理后的图像在H/α平面分割,然后采用传统的模糊C均值算法进行聚类;在H/α/λ/FCM方法中,平均散射强度也参与了分类,但是聚类方法仍然采用传统的模糊C均值聚类算法。
表1 不同极化SAR图像分类方法的精度比较
从表1和图3可以看出,与H/α/FCM算法和H/α/λ/FCM算法相比具有更好的效果:利用本发明所提方法得到的分类结果中,地物的边界更清晰且连续,散射机制类似的区域得到了较好的区分,最终结果与真实地表更加接近。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (5)
1.基于模糊C均值的全极化SAR图像分类方法,其特征在于:包括如下步骤:
1)对原始全极化SAR数据进行预处理,采用多视处理和Refined Lee滤波方法消除SAR图像中斑点噪声的影响;
2)对预处理后的图像进行Cloude-Pottier极化分解,得到散射熵H、散射角α、平均散射强度λ三个极化参数;
3)根据平均散射强度的值将图像分为三大类,分别为:高散射强度区域、中等散射强度区域和低散射强度区域,将每一大类利用H/α平面进一步分割,共得到24小类;
4)采用层次聚类算法将得到的24类合并为n(0<n≤24)类;
5)计算每个像素点与每一类聚类中心的距离,采用FCM算法调整每一类地物的边界,直至满足目标函数最小,迭代终止并输出分类结果;所述FCM算法中的像素点到聚类中心的距离d采用Wishart距离和一个距离因子Wij重新定义:
其中<T>为像素点的相干矩阵,Vm为每一类的平均相干矩阵,即聚类中心,μij代表数据点xi隶属于类别j的概率;所述距离因子Wij表示的含义为Wij越大,每个像素与类别中心Vm的距离越小,离聚类中心近的点变得更近,与聚类中心远的点变得更远。
2.根据权利要求1所述的基于模糊C均值的全极化SAR图像分类方法,其特征在于:
所述步骤2)中,散射熵H、散射角α及平均散射强度λ三个极化参数采用Cloude-Pott ier分解算法对预处理后的数据进行极化分解得到:
λ=λ1p1+λ2p2+λ3p3
其中λi为相干矩阵T的特征值,且∞>λ1≥λ2≥λ3>0,pi的表达式为:
3.根据权利要求1所述的基于模糊C均值的全极化SAR图像分类方法,其特征在于:
所述步骤3)中,根据平均散射强度将极化数据分成三大类,每类的λ边界的确定采用取中值的方法,即根据λ值的动态范围将图像分成均匀的三部分,每部分分别对应高散射强度区域、中等散射强度区域和低散射强度区域。
4.根据权利要求1所述的基于模糊C均值的全极化SAR图像分类方法,其特征在于:
所述步骤4)中,利用层次聚类算法合并过分割的类别,公式如下:
其中Ti与Tj代表两个类别的类别中心,Dij为Wishart距离,表示类别i与类别j之间的分离度;根据所选试验区域地物覆盖的真实情况,将图像分为n(0<n≤24)类。
5.根据权利要求1所述的基于模糊C均值的全极化SAR图像分类方法,其特征在于:
所述步骤5)中模糊C均值聚类算法的目标函数为:
其中m为任意大于1的实数,C为类别数目,N为图像中像素点的个数,μij代表数据点xi隶属于类别j的概率,cj为第j类的类别中心。d为像素点到聚类中心的距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510007661.1A CN104598922B (zh) | 2015-01-07 | 2015-01-07 | 基于模糊c均值的全极化sar图像分类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510007661.1A CN104598922B (zh) | 2015-01-07 | 2015-01-07 | 基于模糊c均值的全极化sar图像分类方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104598922A true CN104598922A (zh) | 2015-05-06 |
CN104598922B CN104598922B (zh) | 2017-09-05 |
Family
ID=53124694
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510007661.1A Expired - Fee Related CN104598922B (zh) | 2015-01-07 | 2015-01-07 | 基于模糊c均值的全极化sar图像分类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104598922B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105403862A (zh) * | 2015-12-04 | 2016-03-16 | 西安电子科技大学 | 一种基于证据c均值的雷达信号类别分选方法 |
CN106295688A (zh) * | 2016-08-02 | 2017-01-04 | 浙江工业大学 | 一种基于稀疏均值的模糊聚类方法 |
CN106408587A (zh) * | 2016-09-19 | 2017-02-15 | 辽宁工程技术大学 | 多视sar图像分割方法及装置 |
CN106683261A (zh) * | 2016-12-22 | 2017-05-17 | 深圳怡化电脑股份有限公司 | 一种纸币的鉴伪方法和装置 |
CN108364011A (zh) * | 2018-02-05 | 2018-08-03 | 中国民航大学 | PolSAR图像多级特征提取与无监督分类方法 |
CN109740475A (zh) * | 2018-12-25 | 2019-05-10 | 杭州世平信息科技有限公司 | 一种遥感图像地面场景分类方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070269078A1 (en) * | 2001-10-15 | 2007-11-22 | Lockheed Martin Corporation | Two dimension autonomous isotropic detection technique |
CN101320431A (zh) * | 2008-06-13 | 2008-12-10 | 北京航空航天大学 | 基于局域标准差的自适应sar图像分类方法 |
CN104123563A (zh) * | 2014-07-17 | 2014-10-29 | 电子科技大学 | 一种基于Cloude分解的极化SAR图像分类方法 |
-
2015
- 2015-01-07 CN CN201510007661.1A patent/CN104598922B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070269078A1 (en) * | 2001-10-15 | 2007-11-22 | Lockheed Martin Corporation | Two dimension autonomous isotropic detection technique |
CN101320431A (zh) * | 2008-06-13 | 2008-12-10 | 北京航空航天大学 | 基于局域标准差的自适应sar图像分类方法 |
CN104123563A (zh) * | 2014-07-17 | 2014-10-29 | 电子科技大学 | 一种基于Cloude分解的极化SAR图像分类方法 |
Non-Patent Citations (2)
Title |
---|
刘利敏等: "一种利用地物散射特性进行后续类别调整的极化SAR影像分类方法", 《光谱学与光谱分析》 * |
曹芳等: "基于Cloude-Pottier目标分解和聚合的层次聚类算法的全极化SAR 数据的非监督分类算法研究", 《电子学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105403862A (zh) * | 2015-12-04 | 2016-03-16 | 西安电子科技大学 | 一种基于证据c均值的雷达信号类别分选方法 |
CN106295688A (zh) * | 2016-08-02 | 2017-01-04 | 浙江工业大学 | 一种基于稀疏均值的模糊聚类方法 |
CN106295688B (zh) * | 2016-08-02 | 2019-10-18 | 浙江工业大学 | 一种基于稀疏均值的模糊聚类方法 |
CN106408587A (zh) * | 2016-09-19 | 2017-02-15 | 辽宁工程技术大学 | 多视sar图像分割方法及装置 |
CN106408587B (zh) * | 2016-09-19 | 2018-11-09 | 辽宁工程技术大学 | 多视sar图像分割方法及装置 |
CN106683261A (zh) * | 2016-12-22 | 2017-05-17 | 深圳怡化电脑股份有限公司 | 一种纸币的鉴伪方法和装置 |
CN106683261B (zh) * | 2016-12-22 | 2019-07-26 | 深圳怡化电脑股份有限公司 | 一种纸币的鉴伪方法和装置 |
CN108364011A (zh) * | 2018-02-05 | 2018-08-03 | 中国民航大学 | PolSAR图像多级特征提取与无监督分类方法 |
CN109740475A (zh) * | 2018-12-25 | 2019-05-10 | 杭州世平信息科技有限公司 | 一种遥感图像地面场景分类方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104598922B (zh) | 2017-09-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104598922A (zh) | 基于模糊c均值的全极化sar图像分类方法 | |
CN103632363B (zh) | 基于多尺度融合的对象级高分辨率遥感影像变化检测方法 | |
Wang et al. | Unsupervised classification of fully polarimetric SAR images based on scattering power entropy and copolarized ratio | |
CN106897679A (zh) | 一种基于改进模糊c均值聚类的语义变化检测方法及*** | |
CN102982338B (zh) | 基于谱聚类的极化sar图像分类方法 | |
CN102999761B (zh) | 基于Cloude分解和K-wishart分布的极化SAR图像分类方法 | |
CN103971115A (zh) | 一种基于NDVI和PanTex指数的高分辨率遥感影像新增建设用地图斑自动提取方法 | |
Xie et al. | PolSAR image classification via D-KSVD and NSCT-domain features extraction | |
CN104166859A (zh) | 基于ssae和fsals-svm极化sar图像分类 | |
CN104867150A (zh) | 遥感影像模糊聚类的波段修正变化检测方法及*** | |
CN103578110A (zh) | 基于灰度共生矩阵的多波段高分辨率遥感影像分割方法 | |
CN103839075B (zh) | 一种基于联合稀疏表示的sar图像分类方法 | |
CN105718942A (zh) | 基于均值漂移和过采样的高光谱图像不平衡分类方法 | |
CN107742133A (zh) | 一种用于极化sar图像的分类方法 | |
CN102073867A (zh) | 一种遥感图像分类方法及装置 | |
Zhang et al. | Improved SLIC superpixel generation algorithm and its application in polarimetric SAR images classification | |
Ju et al. | A novel fully convolutional network based on marker-controlled watershed segmentation algorithm for industrial soot robot target segmentation | |
CN112989940B (zh) | 基于高分三号卫星sar影像的筏式养殖区提取方法 | |
Corr et al. | Classification of urban SAR imagery using object oriented techniques | |
Wu et al. | Fusing optical and synthetic aperture radar images based on shearlet transform to improve urban impervious surface extraction | |
CN102800101A (zh) | 一种星载红外遥感图像机场roi快速检测方法 | |
CN106934805A (zh) | 基于Gamma滤波的SAR图像超像素分割方法 | |
CN116543300A (zh) | 一种基于语义分割的云-气溶胶层次分类方法 | |
Yang et al. | Supervised land-cover classification of TerraSAR-X imagery over urban areas using extremely randomized clustering forests | |
Sun et al. | Feature space fusion classification of remote sensing image based on ant colony optimisation algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170905 Termination date: 20210107 |