基于改进OGMRF-RC模型的SAR图像分类方法
A SAR image classification method based on an improved OGMRF-RC model
通讯作者: 郭拯危(1963-),女,教授,研究生导师,研究方向为SAR遥感信息处理。Email:gzw@henu.edu.cn。
责任编辑: 张仙
收稿日期: 2020-12-31 修回日期: 2021-04-23
基金资助: |
|
Received: 2020-12-31 Revised: 2021-04-23
作者简介 About authors
李 渊(1996-),男,硕士研究生,主要研究方向为SAR遥感图像处理。Email:
合成孔径雷达(synthetic aperture Radar, SAR)图像分类是遥感应用中的关键技术之一。针对对象高斯-马尔可夫随机场(object-based Gaussian-Markov random field,OGMRF)模型中区域类别标签对分类精度影响的问题,提出了区域类别模糊概率(regional category fuzzy probability, RCFP)标签场方法,使临界对象具有多种类别划分的可能性,避免唯一标签导致的错分类现象。该方法综合考虑区域特征与邻域特征,利用区域边缘信息和后验概率获得RCFP,并将其纳入特征场参数求解过程中,使特征场参数更加接近真实情况,从而提高SAR图像分类精度。以河南省开封市东部约1 400 km2的区域为研究区,采用Sentinel-1卫星SAR图像开展农田、建筑、水域3类地物的分类验证实验,与K-means,FCM,马尔可夫随机场和具有区域系数的OGMRF等方法相比较,所提出方法的总体分类精度达到94.16%,Kappa系数为0.895 7,在5种方法中效果最好。
关键词:
The classification of synthetic aperture Radar (SAR) images is one of the key technologies in the field of remote sensing applications. To address the problem that regional class labels affect the classification accuracy in the object-based Markov random field (OMRF) model, this paper proposes the concept of regional category fuzzy probability (RCFP) label field, which can effectively avoid the misclassification caused by wrong class labels by fully considering the possible categories of a single object. The RCFP of every region can be obtained using the regional edge information and posterior probability according to the features of the region and its adjacent regions. Then it is included in the calculation of feature field parameters to make the feature field parameters highly close to the real conditions of objects. The study area is located in the eastern part of Kaifeng City, Henan Province, covering an area of about 1 400 km2. Sentinel-1 SAR images were used for the classification experiment of farmlands, buildings, and water in the study area, and the performance of the improved method in this study was compared with that of the method of K-means, fuzzy C-means (FCM), MRF, and OGMRF-RC. The experimental results show that the overall accuracy (OA) and the Kappa coefficient of the proposed method are 94.16% and 0.8957 respectively, which are higher than those of other methods.
Keywords:
本文引用格式
李渊, 毋琳, 戚雯雯, 郭拯危, 李宁.
LI Yuan, WU Lin, QI Wenwen, GUO Zhengwei, LI Ning.
0 引言
MRF是一种概率图模型,分别使用特征场和标签场表示图像特征和类别。面向对象的MRF(object-based Markov random field, OMRF)模型采用超像素分割处理,按像素特征将图像划分为同质区域,作为分类的基本单元,以区域邻接图(region adjacency graph, RAG)[10] 表示各区域概率关系,并以此进行OMRF建模,将图像分类转化为区域类别最优标记问题。与传统MRF模型相比,OMRF方法用于SAR图像分类可有效利用区域对象大小、形状和边缘等信息,降低相干斑噪声对分类结果的干扰。
近年来,基于OMRF模型的图像分类方法不断涌现各种改进形式。Chen等[11] 建立了统一的MRF(unified Markov random field, UMRF)模型,将像素和区域的似然函数相乘,同时获取像素和区域信息,使得图像分类精度得到提高; Zheng等[12] 提出了一种具有区域系数的对象高斯-马尔可夫随机场模型(object-based Gaussian-Markov random field model with region coefficients, OGMRF-RC),以边缘信息和区域大小为系数,使用基于对象的线性回归方程(object-based linear regression equation, OLRE)获取相邻区域的线性关系,用于特征场建模; Xu等[13] 提出了一种极化辅助场OMRF(OMRF with polarimetric auxiliary fields, OMRF-PA)模型,采用极化指数评估OMRF模型迭代过程中的极化信息损失,结合改进的条件概率分布,有效提高了PolSAR图像分类精度。
然而,上述改进模型中,迭代计算的区域对象标签只有一个类别,这极有可能影响图像中临界区域的类型界定。本文针对该问题,提出了区域类别模糊概率(regional category fuzzy probability, RCFP)的方法,利用边缘信息和后验概率共同获取区域对象对所有可能类别的RCFP,作为OGMRF-RC模型的标签场信息,同时在特征场参数计算中考虑RCFP的影响,并将该改进模型应用于SAR图像分类中。经实验验证,与K-means,FCM和其他MRF模型等方法相比,本文方法在分类过程中以概率形式赋予各区域多类别标签,充分利用了区域对象的图像特征与空间关系信息,使得分类精度有了明显的提升。
1 OGMRF-RC模型
1.1 OLRE邻域关系
设图像I超像素分割后区域集合为
区域
式中:
若区域
式中: 特征参数
1.2 OGMRF-RC模型
OGMRF-RC模型是将图像分类问题转化为给定特征场条件下,求解区域类别标签最大后验概率的问题。根据Bayesian公式,区域类别标签
式中:
将图像分类转化为求解最小化能量的问题,则式(3)可转化为:
式中:
设待分类图像满足高斯分布,则OGMRF-RC模型的目标函数为:
式中: 特征参数
式中
2 改进的OGMRF-RC模型
图1
因而,图像分类时,除了区域对象本身的纹理特征外,将其与邻域区域对象的边缘特征作为分类参数之一,对区域类别的正确划分具有较高的参考价值。
OGMRF-RC模型中,标签场采用对象唯一类别标签的形式。对于不同类别交界的区域对象,一旦出现错误的类别标签将会直接影响特征场参数的计算和最终分类结果的准确性。针对此问题,本文提出依据区域对象自身的后验概率和邻域相关边缘信息,获得该区域对所有可能类别的RCFP,作为标签场,并参与特征场参数的迭代计算,使特征参数的更新考虑了图像的全局信息,利用图像空间特征不断优化并纠正每个区域对象的类别归属,获得更好的分类效果。
2.1 标签场模糊概率
基于边缘信息获取的区域
式中
后验概率
式中:
式中Z为归一化常数。
区域
2.2 特征场参数求解
用RCFP作为区域的类别标签,会使区域
3 改进模型流程
改进的SAR图像分类算法流程如图2所示。首先采用像素级分类方法与图像超像素分割方法对SAR图像进行处理,用于获取区域对象分割结果及初始类别标签。然后应用改进的OGMRF-RC模型对区域对象特征进行迭代计算,不断更新标签场信息,直至满足结束条件完成对SAR图像分类,获得分类结果。算法执行结束条件为: 达到最大迭代次数T,算法结束。
图2
4 实验与结果
4.1 研究区域及数据
图3
表1 实验使用的Sentinel-1卫星SAR图像参数
Tab.1
参数 | SAR图像 |
---|---|
产品类型 | SLC |
成像模式 | 干涉宽幅 |
波段 | C |
成像日期 | 2020年2月27日 |
分辨率 | 4 m×20 m |
图像尺寸 | 1 900像素×1 900像素 |
极化方式 | VV极化 |
实验数据为一景Sentinel-1卫星SAR图像,在图像覆盖区域内,选取若干农田、建筑、水域等不同类别地物的样本,进行真实地物标记,样本分布如图3(c)所示。
4.2 实验结果
为验证本文方法的有效性,选择基于像素的K-means,FCM和MRF模型,以及基于区域的OGMRF-RC模型等方法进行对比实验。实验中本文改进的OGMRF-RC方法选用MRF模型作为像素级分类方法,采用简单线性迭代聚类(simple linear iterative clustering, SLIC)[15] 方法对图像进行超像素分割,依据经验设置势能参数
图4为5种对比方法的分类结果,其中图4(a)为谷歌地图光学图像。对照真实地物分布情况,中部区域不存在水域,但图4(c)和图4(d)在该区域中存在大量农田错分为水域的情况。左上矩形区域内地物多为农田,仅有少量的建筑和水域,但图4(b)—(e)中该区域存在一些地物错分的情况。直观对比不同方法的分类结果,图4(f)所示的本文方法分类结果与真实地物图最为接近。为了更清晰地展示算法改进后的分类效果,选取了图4左上矩形区域进行放大展示,结果如图5所示。其中图5(a)为所选区域的谷歌地图光学图像,图5(b)—(f)为所选区域SAR图像的分类结果。结果显示图5(b)—(e)的分类结果中,农田和建筑存在一些错误分类,本文方法(图5(f))对选定区域的分类结果最接近地物真实情况。
图4-1
图4-1
5种对比方法的SAR图像分类结果
Fig.4-1
SAR image classification results of five methods
图4-2
图4-2
5种对比方法的SAR图像分类结果
Fig.4-2
SAR image classification results of five methods
图5
4.3 定量分析
总体精度(overall accuracy,OA)、Kappa系数和用户精度(user accuracy,UA)常用作分类精度评价指标,其计算方法分别为:
式中:
5种方法的OA和Kappa定量分析结果如图6所示,对不同类别地物的UA如表2所示。图6中,3种MRF模型分类结果的OA和Kappa系数都高于FCM方法和K-means方法,这也进一步验证了MRF模型对于SAR图像分类的适用性。而本文方法由于考虑了对象多类别属性概率,获得了最高的分类精度,OA为94.16%,Kappa系数为0.895 7。由表2的UA结果可知,5种方法对农田的分类结果差别不大; 对水域的分类结果,由于FCM和MRF方法中存在将农田错分为水域的情况,因此这两种方法的UA精度较低,而其他方法都取得了较高的精度; 对于建筑的分类结果,其他方法UA处于80%左右,而本文方法UA为89.05%。综合考虑本文方法的分类效果最好。综上所述,与其他方法相比,本文方法的分类结果,对于水域、农田和建筑物等都具有较好的效果,量化指标OA和Kappa结果同样优于对比算法,证明了本文方法的有效性。
图6
表2 5种对比方法的UA结果
Tab.2
地物类型 | K-means | FCM | MRF | OGMRF-RC | 本文方法 |
---|---|---|---|---|---|
农田 | 0.934 8 | 0.953 7 | 0.952 0 | 0.954 4 | 0.959 4 |
建筑 | 0.754 2 | 0.798 3 | 0.824 4 | 0.829 0 | 0.890 5 |
水 | 0.969 0 | 0.509 4 | 0.838 6 | 0.954 3 | 0.955 9 |
5 结论
本文提出了一种基于改进OGMRF-RC模型的SAR图像分类方法,该方法中充分考虑区域对象自身与其邻域的特征信息,提出RCFP的方法。利用区域后验概率与边缘信息共同计算RCFP,作为区域对象标签场,并纳入特征场参数计算过程,让参数的计算考虑了图像的全局信息,使SAR图像分类精度得到了有效的提高。
为验证本文方法的有效性,采用K-means,FCM,MRF,OGMRF-RC等方法开展对比实验。结果表明,本文方法对水域、农田、建筑等地物的总体分类精度最高可达0.941 6,Kappa系数最高为0.895 7,优于其他对比方法,验证了本文方法的有效性。
本文在进行特征场参数计算时,只利用了图像的灰度特征和纹理特征。在后期的研究中,可以在分类过程中增加极化等特征,从而进一步提高SAR图像的分类精度。
参考文献
基于小波域广义高斯分布的SAR图像分割算法
[J].
SAR image segmentation algorithm based on generalized gauss distribution in wavelet domain
[J].
基于三聚类中心K-Means算法的SAR船只检测方法
[J].
A ship detection method for SAR imagery based on three-centroid K-Means clustering algorithm
[J].
基于改进FCM与MRF的SAR图像分割
[J].
SAR image segmentation based on improved FCM and MRF
[J].
基于SVM+SFS策略的多时相紧致极化SAR水稻精细分类
[J].
Fine classification of rice with multi-temporal compact polarimetric SAR based on SVM+SFS strategy
[J].
Classification of very-high resolution SAR images of urban areas using copulas and texture in a hierarchical markov random field model
[J].DOI:10.1109/LGRS.2012.2193869 URL [本文引用: 1]
SAR image classification through information-theoretic textural features,MRF segmentation,and object-oriented learning vector quantization
[J].DOI:10.1109/JSTARS.4609443 URL [本文引用: 1]
A bilevel contextual MRF model for supervised classification of high spatial resolution remote sensing images
[J].
Polarimetric SAR image semantic segmentation with 3D discrete wavelet transform and Markov random field
[J].DOI:10.1109/TIP.2020.2992177 URL [本文引用: 1]
Hyperspectral image classification with transfer learning and Markov random fields
[J].DOI:10.1109/LGRS.8859 URL [本文引用: 1]
Hierarchical classification of land-cover types using RAG-based merging
[C]
Image segmentation using a unified Markov random field model
[J].
Segmentation for remote-sensing imagery using the object-based Gaussian-Markov random field model with region coefficients
[J].DOI:10.1080/01431161.2018.1563841 URL [本文引用: 1]
PolSAR image classification based on object-based Markov random field with polarimetric auxiliary label field
[J].DOI:10.1109/LGRS.8859 URL [本文引用: 1]
地理相似性:地理学的第三定律?
[J].
DOI:10.12082/dqxxkx.2020.200069
[本文引用: 1]
定律对知识的理解和传播具有极大的推动作用。本文通过对比地理学中的定律和物理学中的定律的异同之处,认为地理学的定律虽然在形式上和提出的方式上与经典的力学定律有所区别,但仍符合人们对“定律”的定义。在共识地理学定律存在的前提下,本文阐述了地理学的第一定律和第二定律中所体现出的地理学定律的特征——描述性和归纳性。基于这些对地理学定律所具有的特征的认识,本文首先介绍了Zhu等从空间推测的角度将所提出的地理学第三定律,即“地理环境越相似,地理特征越相近”,或称地理相似性定律。本文从普遍性、独立性和应用性3个方面探讨了地理相似性是否可以作为地理学的一个定律,认为地理相似性与地理学第一和第二定律相比,有着同样的普遍性,且与第一和第二定律所涉及的规律有着本质上的区别,在应用上可能解决第一和第二定律所面临的挑战,因此应可作为地理学第三定律。
Geographic similarity:Third law of geography?
[J].
/
〈 |
|
〉 |
