魏 丹,李 渊,黄 丹

1(河南大学 计算机与信息工程学院,开封 475004)

2(河南大学 实验室与设备管理处,开封 475004)

近年来,随着极化SAR 系统的不断发展,可用的极化SAR 数据逐渐增多,分辨率也在不断提高.与单极化SAR 相比,极化SAR 进行的是全极化测量,可以获取更丰富的目标信息,因此极化SAR 的应用也越来越广泛[1].采用极化SAR 图像进行地物分类是极化SAR 技术在实际应用中的一个重要发展方向,也是SAR图像解译的重要研究内容.

通常地表地物类型是草地、森林、农田、城区、水域和裸地等,采用极化SAR 图像对这些地物类型的分类研究可以在民用领域和军用领域为人们提供帮助.在民用领域可以对农作物的生长、城市的变迁等进行监测,也可以测量土壤水分和对矿产资源进行勘探;在军用领域,可以识别战略目标,为之后的战略部署等方面提供理论依据.由以上可知,极化SAR 图像已经成为地理环境监测的重要组成部分.目前在区域乃至全球的地表地物分类中,采用极化SAR 图像进行地物分类表现出巨大的应用前景和潜力.

极化SAR 系统按照载体的不同可以分为机载SAR[2-8]和星载SAR[9-15],部分机载和星载极化SAR系统如表1和表2所示.表2中,高分三号(GF-3)卫星[15]是中国首颗民用高分辨率星载极化SAR 成像卫星,于2016年8月10日发射升空,空间分辨率从1 米到500 米,幅宽从10 公里到650 公里,有聚束和超精细条带等12 种成像模式.随着极化SAR 数据的增多,如何有效地利用其数据特征实现更精确地分类成为了当今研究的热门问题.

表1 典型机载极化SAR 系统

表2 典型星载极化SAR 系统

本文结构安排如下:第1 节对近几年极化SAR 图像地物分类方法及其应用进行总结和论述,第2 节对极化SAR 图像地物分类的发展趋势进行分析,第3 节对本文进行总结.

1 极化SAR 图像地物分类的方法

当雷达发射的电磁波接触地物时,不同地物对应的雷达回波信息也不同,所以极化SAR 图像中每个像元所展示的信息也不同.对极化SAR 图像进行地物分类的实质是将图像的所有像元按其性质分为若干个类别的过程,从而达到对地物进行分类的目的.

极化SAR 图像进行地物分类的基本步骤为:预处理→特征提取→分类处理.预处理通常包括斑点滤波、几何校正和去取向角等.特征提取一般为纹理特征提取和极化特征提取.纹理特征的提取方法有局部傅立叶变换、曲波变换、Gabor 小波和灰度共生矩阵等[16].极化特征可以通过目标分解获得,其中比较著名的分解有Pauli 分解[17]、Cloude 分解[17]、Krogager 分解[18]、Freeman 分解[19]、SSCM 分解[20]、Yamaguchi 四分量分解[21,22]、Touzi 分解[23]和六分量分解[24]等.

对极化SAR 图像进行地物分类处理的方法多种多样,根据所用特征信息的不同,可以将分类方法归结为基于极化目标散射特性的分类、基于统计分析的分类、基于目标散射特性和统计方法相结合的分类.根据分类过程中是否需要人工标记样本,可以分为监督分类、无监督分类和半监督分类3 大类.根据处理对象的不同可以分为基于像素级和基于对象级的分类.根据采用的技术方法的不同可以将极化SAR 图像地物分类方法分为6 种,分别是主动轮廓模型法、马尔可夫随机场模型法、模糊理论法、支持向量机法、神经网络法和融合算法,如图1所示.

图1 极化SAR 图像地物分类方法

本节从技术方法的角度出发,对近年来不同极化SAR 图像地物分类方法进行总结和论述.

1.1 基于主动轮廓模型的分类

1988年,Kass 等提出了一种主动轮廓模型(Active Contour Model,ACM),该模型主要原理是通过初始化曲线,构造能量方程式,最小化曲线的内外部能量,使曲线逐渐接近待检测物体的边缘,分割出目标.ACM方法中,按照不同的轮廓曲线表达式可分为参数主动轮廓模型(Parametric ACM,PACM)和基于水平集的几何主动轮廓模型(Geometric ACM,GACM);按照曲线的演化方式又可分为基于边界的主动轮廓模型(Edgebased ACM,EACM) 和基于区域的主动轮廓模型(Region-Based ACM,RACM)[25].ACM 方法有很强的抗噪声能力,所以也被应用于极化SAR 图像分类中[26-30].但是该方法的缺点是计算量大,效率低.

2014年,清华大学的杨健团队采用多时相水平集主动轮廓模型对AIRSAR 获取的荷兰弗莱福兰极化数据进行分类处理[30].图2为分类结果,其中图2(a)为C 波段、L 波段和P 波段的波段合成图,图2(b) 为ACM 分类结果.

图2 基于ACM 的分类结果[30]

1.2 基于MRF 模型的分类

由于相干斑噪声的存在,基于单个像素所提供的信息对极化SAR 图像进行地物分类,通常不能得到令人满意的结果.为了获得可靠和满意的结果,分类最好是基于集群的统计度量,而不是基于单个像素的度量.马尔可夫随机场(Markov Random Field,MRF)模型是一种概率图形模型,它可以统计像素之间的空间上下文信息,是处理图像分类问题的强有力工具.在MRF模型中,图像特征被提取为特征场,图像中所有像素类别标记的集合为标记场,在已知特征场的前提下,分类结果通过最大标记场概率获得.MRF 因其固有的抗噪声能力和减少类内差异的影响而被广泛采用.针对极化SAR 图像地物分类问题,各国学者提出了一系列基于MRF 模型的分类方法[31-35].

2014年,Doulgeris 以Radarsat-2 在C 波段采集的美国旧金山海湾极化SAR 图像为实验数据,采用MRF 模型提取上下文信息,对地物进行分类[35].图3为MRF 分类结果,其中图3(a)为Pauli 图像,图3(b)为分类图像,从分类图像中可以看出有取向角的建筑物和森林植被可以被区分开来,提高了分类精度.

图3 基于MRF 的分类结果[35]

1.3 基于模糊理论的分类

传统的极化SAR 图像地物分类方法是为每个像元指定一个特定的类,但是极化SAR 图像中的一个像元可能包含多个地物类别,对这类像元的分类会影响最终的分类精度,采用模糊理论可以较好的解决此问题.模糊理论思想是用0-1 之间的值表示每个像元属于每个类的程度,这些值被称为隶属度,每个像元属于所有类的隶属度之和为1,最终用隶属度的大小确定像元所属的分类.经过多年发展,模糊理论已经广泛用于极化SAR 图像地物分类[36-40].本节将介绍FCM (Fuzzy C-Means,FCM)和最近邻方法在极化SAR 图像地物分类中的应用.

1.3.1 基于FCM 的分类

FCM 算法的思想是把n个数据集合A={a1,a2,···,an}分成C个模糊组,对每个数据集可用值在0～1 之间的隶属度来确定其属于各个组的程度,通过迭代使得非相似性指标的价值函数达到最小,FCM 流程图如图4所示.

2009年,杨然等用模糊分类方法下的FCM 算法与Yamaguchi 分解相结合,提出了一种新的分类模型[39].该模型把四个散射分量组成一组归一化的特征矢量,进行FCM 聚类分析,并且用日本机载L 波段PiSAR数据验证了该算法具有较高的分类精度和较好的视觉效果.图5是FCM 分类结果,其中图5(a)为Niigata 地区的极化合成图,A 区为森林,B 区为人工建筑,C 区为河流,D 区和E 区都是耕地.图5(b)是FCM 分类结果.

图4 FCM 流程图

图5 基于FCM 的分类结果[39]

1.3.2 基于最近邻的分类

最近邻分类的原理是通过确定分类类别,选择正确的不同类别样本,计算测试样本与所有样本之间的距离,并将距离转换为0-1 之间隶属度,如果隶属度大于设定的隶属度阈值,则测试样本属于与它距离最近的样本所属的类别.

2018年,陆翔等采用极化目标分解得到的极化散射特性和用灰度共生矩阵提取的纹理信息作为特征,应用模糊分类方法下的最近邻分类器对TerraSARX 在X 波段采集的云南西双版纳勐腊县极化SAR 图像进行地物分类,图6为分类结果[40].其中图6(a)是勐腊县光学影像,图6(b)为地物分类结果图,图6(c)为不同地物的颜色表示.

图6 基于最近邻的分类结果[40]

1.4 基于支持向量机的分类

支持向量机(Support Vector Machines,SVM)是一种性能很强的分类器,它的基本思想为:建立一个分类超平面作为决策曲面,使得不同样本的隔离边缘最大.如图7所示,圆圈和正方形分别代表了两种类型的样本,H 为分类超平面,H1 与H2 为通过各类样本中最边缘的点,且与H 相互平行、距离相等的平面,它们之间的距离是分类间隔.最优分类超平面就是以最大间隔将两类样本分开的超平面.1999年,Vapnik 团队第一次将SVM 用到图像的分类当中[41].经过多年的发展,SVM 也广泛用于极化SAR 图像的分类中[42-47].

2019年,张腊梅提出了一种改进的SVM 算法,该算法结合了Pinball 损失函数,并将多成分目标分解,H/A/ α分解以及灰度共生矩阵方法分别提取到的12 种极化特征和纹理特征作为SVM 的训练样本[42].这里的12 种特征分别是奇次散射功率Ps、偶次散射功率Pd 和体散射功率Pv、极化散射熵H、极化各向异性度A、平均散射角 α、能量ENY、熵ENT、对比度CON、反差矩IDM、相关度COR 及和平均SAV.并对EMISAR 在Foulum 地区获取的L 波段极化SAR数据进行分类处理,把地物分为建筑物、森林、裸地、细径作物、阔叶作物5 种类型.图8为分类结果,其中图8(a)为Pauli 图像,图8(b)为改进SVM 分类图像,图8(c)为不同地物的颜色表示.

图7 最优分类超平面示意图

图8 基于SVM 的分类结果[42]

1.5 基于神经网络的分类

传统分类算法在进行极化SAR 图像地物分类时,需要人工提取分类特征,而神经网络是一种以人脑为模型的机器学习算法,可以自动的提取分类特征,在特征学习时不需要过多的人工干预,通过对大量样本进行内部特征分析,可以得到更好的分类效果.1991年,Pottier 等第一次利用人工多层感知神经网络对极化SAR 图像进行了分类研究.经过多年的发展,神经网络算法已经广泛用于极化SAR 图像的分类当中[48-59].

2018年,张腊梅等将传统卷积神经网络改进为三维卷积神经网络(3D CNN)[57].并对UAVSAR 系统采集的美国旧金山海湾某港口极化SAR 图像进行地物分类,该地区地物种类丰富,共有11 种地物类型,包括浅海、海洋、森林、人造目标、道路、植被、平坦裸地、粗糙裸地及3 种农田地区,分类结果如图9所示,其中图9(a)是Pauli 图像,图9(b)是三维卷积神经网络分类图像.

图9 基于3D CNN 的分类结果[57]

2019年,刘旭等提出了采用极化卷积神经网络(Polarimetric Convolutional Network,PCN) 对极化SAR 图像进行地物分类的方法[58].该方法采用极化散射编码方式来处理极化散射矩阵,该编码模式可以完全保持散射矩阵的极化信息,避免了将矩阵分解成一维向量,并获得一个二维矩阵.针对这种编码模式,刘旭等人设计了一种基于卷积神经网络的算法,并将其命名为PCN,用于对极化SAR 图像进行分类处理.图10展示了RadarSat-2 在C 波段获取的荷兰弗莱福兰极化SAR 图像进行分类处理的结果.其中图10(a)是荷兰弗莱福兰Pauli 分解图像,图10(b)是PCN 分类结果,图10(c)为不同地物的颜色表示.结果表明该方法取得了较好的分类效果,具有巨大的极化SAR 数据分类潜力.

1.6 基于融合算法的分类

上述各种算法都有其固有的特点,在极化SAR 图像地物分类中也都取得了有效的研究成果,推动了极化SAR 图像地物分类的发展,基于上述分类方法的不同特点,研究者们提出将两种或多种算法进行融合,以提高分类精度[60-66].

图10 基于PCN 的分类结果[58]

2018年,Xie 等提出了一种新的极化SAR 图像分类模型,该模型首先结合Wishart 分布与自动编码器(Auto-Encoder,AE)神经网络提出了WAE (Wishart-Auto-Encoder)分类模型,之后又结合K-means 聚类算法提出了Clustering-WAE 分类模型以提高分类能力[62].并对RadarSat-2 采集的西安地区极化SAR 图像进行地物分类处理,结果如图11所示.其中图11(a)是西安地区Pauli 分解图,图11(b)～图11(e)分别是K-means算法、Wishart 算法、AE 算法、WAE 算法分类结果,图11(f)是Clustering-WAE 融合算法分类结果,图11(g)为不同地物的颜色表示.分类结果表示,Clustering-WAE 融合算法分类效果要优于其他算法.

2018年,Guo 等将模糊理论应用到超像素提取中,采用FCM 算法并对ESAR 获取的德国奥伯菲芬-霍芬地区极化SAR 图像进行模糊超像素提取,之后采用SVM 对图像进行地物分类[63],结果如图12所示.其中图12(a)是Pauli 图像;图12(b)是地面真实图;图12(c)超像素提取结果图;图12(d)为分类结果与地面真实地物比较图,图中黑色部分表示分类结果与地面真实地物不一致;图12(e)是不同地物的颜色表示.

1.7 各方法优缺点

上述各方法在极化SAR 图像地物分类中都有其优缺点,具体如表3所示.

图11 西安地区地物分类图[62]

这些方法中,主动轮廓模型因考虑了轮廓内外部信息,可以获得较为平滑的分割轮廓,但是其计算量大,且易受噪声影响,该方法适用于较均匀,场景较小的图像.MRF 模型在图像分类时考虑了像素点的空间上下文信息,所以其抗噪声能力强,但是计算量较大,也只适用于场景较小的图像.模糊理论可以解决图像中的模糊像素问题,但只适用于分辨率较低的图像.神经网络整体分类精度高,对噪声有很强的鲁棒性,但是需要大量样本集进行训练,所以适用于有大量样本集的场景.相比于神经网络,支持向量机适用于小样本集的场景,但是该方法中的一些参数需要凭实验和经验来进行设置.融合算法可以综合多种算法的优点,但是其算法复杂度较高,根据其算法的不同,不同的融合算法也适用于不同的场景.

图12 奥伯菲芬-霍芬地区地物分类图[63]

表3 分类方法优缺点

2 极化SAR 图像地物分类发展趋势

在极化SAR 图像分类中,基于SVM、MRF 模型、模糊理论、神经网络等的分类方法仍是主流方法,也取得了不错的分类效果,但是仍然存在一些困难,如计算量大、分类特征不多、不适用于复杂环境等.而且在进行分类时,分类的效果会与相干斑噪声、雷达的工作波段以及地物的类型等因素有关.随着极化SAR图像地物分类技术的发展,更多研究人员将深度学习技术,多波段融合技术引入到极化SAR 图像地物分类当中,并着力于发现新的分类特征.这些不断涌现的新思路、新方法指示着未来极化SAR 图像地物分类的发展趋势,大致可以归为以下几点.

2.1 新的分类特征

在极化SAR 图像地物分类中,特征提取是其中重要的步骤,常用的分类特征有纹理特征和极化特征.基于这些特征,研究者们结合着分类器对极化SAR 图像进行地物分类,近年来,一些新的分类特征不断涌现,以提高分类性能.

2014年,Uhlmann 等利用伪彩图提取强大的颜色特征,并首次将其用于极化SAR 图像地物分类,结果表明,颜色特征具有较好的鲁棒性,受图像的噪声等因素的影响非常小[67].之后,研究者们也逐渐将颜色特征引入到极化SAR 图像的分类当中.如沈璐等于2017年提取颜色特征和极化特征组成特征向量,采用SVM 分类器对极化SAR 图像进行分类[68].除颜色特征外,Liu 等于2019年利用极化散射编码提取稀疏散射编码特征,这种特征可以保持数据的空间结构[58].

以上研究者们所提出的新的分类特征在实际应用中都取得了非常显著的分类效果,所以新的分类特征的发现与引入将是极化SAR 图像地物分类中的主要发展趋势.

2.2 深度学习

采用传统的神经网络对极化SAR 图像进行地物分类虽然能提取出分类特征,但是没有考虑到图像本身的二维空间信息,而深度学习使用多层神经网络对图像进行特征学习,可以提取出更高级的特征,充分挖掘数据信息,很大程度上提高了分类精度,分类结果也表现更加稳定.深度学习也是目前极化SAR 图像地物分类中,最常用和分类效果最好的方法.

目前,有大量研究者们采用深度学习方法对极化SAR 图像进行地物分类.如2016年,Liu 等采用深度置信网络(Deep Belief Network,DBN)对西安地区极化SAR 图像进行地物分类[69].2016年,Jiao 等结合深度堆叠网络(Deep Stacking Network,DSN)和Wishart距离提出了W-DSN 模型,并对荷兰弗莱福兰极化SAR图像进行分类[60].2019年,Liu 等采用生成对抗网络(Generative Adversarial Network,GAN)对西安地区极化SAR 图像进行地物分类[70].这些深度学习方法在极化SAR 图像地物分类过程中都取得了非常精确的分类效果.

由上述可知,深度学习方法已经成为目前极化SAR 图像地物分类的主流方法.可是深度学习依靠大量样本的学习,目前公开的极化SAR 图像地物样本数据集少之又少,研究者们应首先考虑如何做出一个具有代表性的极化SAR 图像地物样本数据集.

2.3 多波段融合

极化SAR 图像成像信息不仅与地物类型有关,还与雷达发射波的波长有关,不同波段的极化SAR 其传输特性、目标回波的后向散射特性等不同,相同的物体也有可能呈现不同的表现形式[71].在地物分类方面,不同波段的极化SAR 数据,在实际测量中有着不同的优势,如高频波段极化SAR 图像对地物的纹理描述清晰,有利于区分纹理差异较大的地物,如城区和裸地等.低频波段极化SAR 图像对地物湿度和植被情况比较敏感,有利于区分具有不同含水量的地物,如旱地和沼泽地等.

综上所述,不同波段的极化SAR 系统在实际测量中都有其特点和优势,因此研究者们利用信息融合技术将多波段极化SAR 图像进行波段融合,获取更丰富的目标信息以获得更好的地物分类效果.如2014年,杨健团队将C 波段、L 波段和P 波段极化SAR 图像进行融合,之后采用水平集方法对地物进行分类,结果表明相比于单波段极化SAR 数据,采用组合多波段数据进行地物分类的效果更好[30].因此,将同一地区多波段数据进行融合也是未来极化SAR 图像地物分类中的重要发展方向.

但是,由于目前绝大多数星载极化SAR 系统的工作波段只有一个,所以获取同一地区的星载极化SAR多波段数据就变得尤为困难,就极化SAR 系统的发展趋势而言,多波段星载极化SAR 系统是未来的主要发展趋势.

3 结束语

本文首先介绍了采用极化SAR 图像进行地物分类的好处和重要性,总结了近年来在国内外公开发表的基于极化SAR 图像进行地物分类的研究成果,并对其研究成果进行了简要概述.基于这些方法的大量研究,许多创新方法也在不断涌现,如引入新的分类特征;将深度学习技术引入到分类当中;利用信息融合技术将多波段数据进行波段融合,这些方法都极大推动了极化SAR 图像分类的发展.

纵观极化SAR 图像地物分类史,新的研究成果在不断增加地物分类的准确度,也在不断增强着地物分类的实用性.但是仍然存在一些问题,如怎样选取最适合的分类特征、同一地区多波段数据不多、如何更好地去除相干斑噪声等,因此要加强对极化SAR 技术的研究,获取更深层次的成果.相信在不久的将来,极化SAR 图像分类方法研究会越来越完善,分类效果越来越精确.