洪伟, 赵祥模, 王鹏, 李晓艳, 邸若海, 吕志刚, 王储

1.西安工业大学兵器科学与技术学院， 陕西西安 710021；2.西安工业大学电子信息工程学院， 陕西西安 710021；3.西安工业大学发展规划处， 陕西西安 710021

遥感(remote sensing),指一切非接触的超视距探测,一般指运用传感器对物体辐射及反射特性的探测。伴随着遥感技术不断更新迭代,遥感图像分辨率不断提升,遥感数据中常见的桥梁、飞机等目标也可在图像中清晰可见。由此可见遥感图像可以从宏观的角度反馈信息。因此准确地对遥感图像进行图像分析在军事、民生保障、气象测绘等领域十分重要[1]。

在深度卷积网络的框架下,遥感图像检测中大部分现有方法依赖于复杂的R-CNN框架,如SCRDet[2]、CAD-Net[3]、ReDet[4],该框架由两部分组成:区域建议网络(RPN)及R-CNN检测头[5]。RPN用于从水平锚框生成高质量的感兴趣区域(ROI),并利用ROI池化算子从ROI中提取准确的特征,最后利用R-CNN进行目标分类以及边界框回归。但水平ROI通常会导致定向边界框和目标之间严重错位[6],遥感图像中目标分布较为密集,水平ROI通常会包含多个实例。因此,需要具有不同角度、长宽比的符合需求的锚框,但生成锚框的同时会导致计算开销大幅增加。近期,ROI转换器的提出[7],利用水平ROI替换旋转ROI的方法,解决了锚框重复生成的问题,但仍需复杂的ROI操作。与基于R-CNN的检测器相比,单阶段检测器,如YOLO-ACN[8]、S2ANet[9]、SSD[10]、RetinaNet[11],通过回归边界框,并使用规则且密集的采样锚点直接对目标进行分类。这种体系结构具有较高的计算效率,但在精度上往往低于双阶段检测器。现阶段,遥感图像目标检测算法大多由传统双阶段目标检测算法针对遥感图像特点改进得到。传统算法虽能在常规场景的数据集上表现出良好的检测性能,但是由于遥感场景中存在目标成像小、变化方向任意、变化尺度大、分布密集、背景复杂等问题,改进后的传统算法仍难以有效提取遥感场景中目标的细节位置特征及充足的语义信息,导致检测效果无法达到预期[12]。

为提升遥感场景下目标识别能力,考虑到遥感图像目标变化方向的不确定性,本文针对遥感图像目标成像小、背景复杂、分布拥挤的问题提出基于感知延伸与锚框最适匹配的遥感图像目标检测算法,主要工作如下:①在特征融合前利用协同注意力模块(SEA),捕获特征像素间关系的同时扩展模型感知区域,实现目标与全局的关系建模;②通过分块裁剪输入图像以扩大目标在图像中所占比例,利用感知延伸特征金字塔(HQFPN)完成特征融合,HQFPN下采样过程中将感知增强卷积模块与常规卷积交替构建感知延伸下采样模块(EN-DilatedConv),减少下采样过程中低层细节位置特征的丢失,解决传统特征融合金字塔模块[13]在特征融合时信息丢失的问题;③在为真实边界框分配锚框时,采用高质量锚框匹配算法(MaxIoUAssigner-HQ),在保证召回率的同时,减少低质量锚框匹配的产生。

1 S2ANet目标检测算法

S2ANet模型的总体框架如图1所示。该模型使用ResNet50[14]模块从输入图像提取特征,并使用包含特征对齐模块和定向检测模块组成的目标检测头网络生成旋转预测边界框及类别置信度,最后使用非极大值抑制(NMS)筛选结果。FAM模块由锚优化网络(ARN)及对齐卷积层(ACL)组成。ARN是一个轻型网络,由锚分类分支及锚框回归分支组成。锚分类分支将锚细分为不同类别,锚回归分支将水平锚细化为高质量旋转锚框。为了加快模型推理速度,推理阶段不执行锚分类分支,转而通过锚框回归调整对齐卷积(AlignConv)中采样位置。

图1 S2ANet结构图

首先通过规则网络Rn={(rx,ry)}对输入特征图X采样,然后将加权的采样值相加。例如R={(-1,1),(-1,0),…,(0,1),(1,1)}表示核大小为3×3,膨胀为1,对于输出特征图Y上的每个p∈R,采样加权计算公式为

(1)

相对于普通卷积,对齐卷积增加偏置项O。

(2)

(3)

式中:k表示为卷积核大小;S表示采样步长,其中R(θ)=(cosθ,-sinθ;sinθ,cosθ)T为旋转矩阵;p位置的偏置项O可表示为

(4)

通过该方式可以将给定位置p的轴对齐卷积特征X(p)转换为基于对应锚框的任意方向卷积特征。

图2a)为普通2D卷积基于锚框生成的采样点,图2b)为经过AlignConv生成的旋转锚框,箭头表示采样点偏移。

图2 普通卷积与对齐卷积采样点

图3 HQ-S2ANet结构图

ACL将ARN生成的旋转锚解码为绝对锚框(x,w,h,θ)后与计算所得偏移量及图像提取特征输入对齐卷积中,以提取对齐特征。ODM模块利用主动旋转滤波器(ARF)缓解分类置信度和位置回归精度的不一致性,使得目标检测更加精准。ARF是一个K×K×N大小的滤波器,利用ARF旋转在卷积操作中生成的特征图具有N个方向,以实现对方向信息的编码。将ARF用于卷积层,可以获得方向敏感特征。ODM通过选择响应最强的方向通道汇集方向敏感特征,达到提取方向不变特征的目的。最后将方向敏感特征和方向不变特征分成2个子网络,分别进行边界框回归及目标分类。FAM与ODM组成的目标检测头,可以对齐具有不同方向的对象特征,从而提升目标分类的鲁棒性。

2 HQ-S2ANet目标检测算法

为解决遥感场景中目标成像小、背景复杂、分布拥挤的问题,本文提出一种改进的S2ANet遥感图像目标检测算法HQ-S2Anet,使用感知延伸特征融合模块对骨干网络提取特征进行融合;为锚框匹配真实标注框时采用MaxIoUAssigner-HQ算法,在保证召回率的同时,抑制低质量匹配的发生,以增强模型检测能力。

2.1 感知延伸特征融合模块

2.1.1 协同注意力模块(SEA)

自注意力机制(self-attention)作为注意力机制的重要组成部分,依赖特征个体相关性,通过接受的输入生成像素间关联性权重,避免外部信息对其的影响。通道注意力ECA-Net[15]在不增加网络计算量的同时使网络增强有用信息去除冗余特征。

为缓解自注意力获取通道间关系能力较弱以及ECA缺乏特征局部特征关系的问题,将自注意力对上下文信息的捕捉能力以及ECA能够加强通道间依赖关系的效果进行结合。输入特征经过1×1卷积后得到3个维度的向量Q(query),K(key),V(value)。将Q与K向量的转置矩阵相乘,经过softmax后得到权重系数A。当矩阵A与矩阵V相乘后得到经过自注意力处理的输入向量权重S后,利用ECA对其进行通道间关系处理。此时得到的矩阵较通过原始输入特征得到的矩阵,加强通道间依赖关系的同时避免了传统通道注意力因降维跨通道交互所带来的影响,由此提出了协同注意力模块SEA。如图4所示,利用协同注意力模块处理主干网络输出特征,以增强特征间相关性。

图4 协同注意力结构图

(5)

将权重系数A乘以V向量,得到加权后的输入特征评分S。

(6)

将S向量reshape为三维特征向量后,经过ECA强化通道间关系,与输入特征向量通过add操作得到SEA模块的输出特征Z。

(7)

2.1.2 感知延伸特征金字塔模块(HQFPN)

为缩短信息路径,利用低层存在的精确位置定位信息增强特征金字塔融合效果,构建自底向上路径。随着特征传递层数的增加,特征信息丢失逐渐增多、局部化响应能力不足。因此,在自底向上的路径中建立从低层到高层的捷径分支,以弥补信息丢失,同时通过传播低层细节的强响应,进一步增强整个特征层次的局部化响应能力。

式中:Conv表示通道为256,核为3×3的卷积操作;Downsample()表示下采样操作;EN-DilatedConv表示感知延伸下采样模块。

传统特征融合模块中,低层特征图细节定位信息较为丰富,高层特征图语义信息较为丰富,语义信息经过自上而下路径传递,中间经过卷积操作及信息交互过程较少,对语义信息影响较小;浅层特征需沿骨干网络传递,细节特征传递至顶层需要经过骨干网络中多次卷积、池化等操作,导致特征细节定位信息丢失较为严重,不利于浅层特征中小目标的定位。因此在原有特征融合模块基础上构建自下而上路径,避免底层特征经过骨干网络,使底层细节位置特征传递至深层时保留充足的有效信息。

图5中红色连接线为传统特征金字塔低层特征传递路径,蓝色连接线为HQFPN传递路径,构建捷径分支路径后,低层特征传递所经过的网络层结构大幅减少,有效地保留了低层特征中的细节位置信息,提升了小目标识别能力。

图5 低层特征传递示意图

遥感图像存在一定随机性,部分目标类别排列稀疏,数据量小。该类遥感图像包含的实际信息十分有限,图像中目标区域较为局限,其余区域多为无关背景,大尺寸卷积操作虽能较好地提取目标语义特征,但是对于遥感图像,扩大卷积核不仅会增加算法计算开销,也会影响算法对密集分布目标的识别能力。空洞卷积在普通卷积操作的基础上通过膨胀系数d控制卷积核稀疏采样尺度,并将卷积核稀疏采样间隔区域填充[16]。

Sw,h=Kw,h+(Kw,h-1)·(d-1)

(11)

式中:Sw,h为空洞卷积核尺寸;Kw,h为原卷积核尺寸。

故在HQFPN自底向上的下采样过程中使用空洞卷积,在不增加有效卷积单元数量的同时延伸感知范围。

图6为标准卷积与空洞卷积操作示意图,空洞卷积根据空洞率将特征间隔部分填充0,导致特征图出现网格隔影[17],一定程度弱化了特征,影响下采样后特征的显著性。为解决该问题,本文构建EN-DilatedConv模块,将空洞卷积与均值池化组合使用。对于低层特征,一部分采用空洞率为2的3×3空洞卷积与输入特征卷积后的特征组成,另一部分则由3×3大小的全局均值池化作用于输入特征图后得到的特征组成,两部分特征进行相加后,作为下采样最终的输出结果

图6 标准与空洞卷积操作示意图

O=F1+F2,O∈RH×W×C

(12)

F1=ConvD(X),F2=AvgPool(X)

(13)

式中:X为输入特征图;ConvD()表示核大小为3,步长为2,空洞系数为2,填充为2的空洞卷积;AvgPool表示全局平均池化操作。

在自下而上特征传递路径中,自最底层开始,对特征图采取下采样操作,逐层传递特征信息,与对应自上而下路径上的特征图相加融合。此外将C2的特征图经过8倍池化后,与自下而上传递路径中最顶层特征图融合,使最顶层输出特征图既包含自身丰富的高阶语义信息,又融合原始低层特征图池化得到的细节信息,达到增强特征图信息表征的能力。

最终将H3,H4,H5,C6,C7作为预测特征图,输出特征图既延伸了感知范围,又包含了低层丰富的细节定位信息及高层充足的语义信息。

2.2 锚框最适匹配方法(MaxIoUAssigner-HQ)

经过感知增强后的特征图传递到检测头中,FAM及ODM会得到生成的高质量旋转锚框与真实标注框的交并比(IOU),根据设定好的正样本及负样本阈值为所有锚框分配属性。在原网络中,S2ANet将IOU值低于负样本阈值的锚框判定为负样本,不参与真实标注框的匹配。由于航空遥感图像目标尺度小、分布拥挤,网络训练前期产生较多的低质量匹配均被忽略,导致召回率过低,不利于网络的损失计算,进而导致精度过低。改进后的S2ANet实现过程中,令低质量锚框也参与真实标注框的匹配,虽保证了召回率,但一定程度上导致低质量匹配的发生。

B为生成锚框(Anchor)个数,G为真实标注框个数。矩阵I记录网络生成锚框与对应真实标注框的IOU值,GA记录真实标注框与Anchor的最大IOU值,Gi记录对应Anchor的索引。A记录Anchor与真实标注框的最大IOU值,Ai记录对应真实标注框的索引。网络生成锚框与真实标注框对比时会对GA中的值做筛选,若小于负样本阈值,则被划分为背景,判定该标注框无符合要求锚框。该操作会使部分真实标注框未被分配锚框,导致召回率过低。为解决这一问题,MaxIoUAssigner设定正样本阈值下限,若IOU小于负样本阈值但是大于正样本阈值下限,则不将其划分为背景,并将该锚框与对应真实标注框做匹配。

由于遥感图像部分目标分布密集,单个锚框会与多个真实目标标注框发生重叠,且由于一个锚框只允许与一个真实标注框完成匹配,高质量匹配可能会被低质量锚框覆盖。若当前锚框已与最大IOU值对应的真实标注框完成匹配。i,j为真实标注框索引,m为Anchor索引。

此时若另一真实标注框的最大IOU也为与该锚框交并得到的

此时

Ai[Gi[m]]=j

(23)

匹配IOU大于该样本阈值,导致与该锚框IOU最大的真实标注框i被重写为j,该锚框被分配了IOU值较低的真实标注框,甚至可能为负样本,导致锚框匹配质量过低,影响边界框回归以及分类损失,进而影响模型精度。针对以上问题,本文提出MaxIoUAssigner-HQ高质量锚框匹配方法。

表1为MaxIoUAssigner-HQ算法伪代码,在为真实标注框分配其IOU最大锚框时通过常数因子控制,若该IOU值为当前锚框与真实标注框交并最大值,则允许匹配。若不为最大值,则与最大值做差值,当差值小于常数因子时,允许匹配发生,当大于常数因子时则拒绝匹配。通过常数因子限定匹配的发生,在保证召回率的同时,防止低质量匹配发生,避免高质量锚框被重写。

表1 MaxIoUAssigner-HQ算法伪代码

3 实验结果与分析

3.1 数据集与评价指标

为验证改进算法的有效性,设计实验在DOTAV1.0数据集上进行验证。DOTA数据集来源包含不同传感器和平台,包括 Google Earth、JL-1卫星拍摄,以及中国资源卫星数据和应用中心的GF-2卫星拍摄。数据集共计2 806幅航拍图,每张图像的像素尺寸在800×800到4 000×4 000的范围内,包含不同尺度、方向和形状的物体。类别包括:plane,baseball-diamond,bridge,ground-t-field,small-vehicle,large-vehicle,ship,tennis-court,basketball-court,storage-tank,soccer-ball-field,roundabout,harbor,swimming-pool,helicopter。由于DOTA测试集未公开相应类别标签文件,本次实验采用官方给定训练集进行训练,不设置验证集,划分官方验证集作为测试集,验证集与测试集无相同图片。采用mAP作为评价指标。

3.2 实验环境及参数设置

实验环境为Ubuntu18.04,Pytorch版本为1.8.1,服务器硬件配置为NVIDIA RTX3070(8GB显存)显卡,Inter core i5处理器。

本文采用的DOTA数据集中图像分辨率集中在800×800到4 000×4 000的范围内,若直接将大分辨率图像输入网络或导致参数量过大,图形计算单元无法负荷。经过实验,将输入分辨率改为2 048×2 048后,网络无法在当前硬件环境下运行。若分辨率过小,将导致可学习特征减少,网络无法获取足够的空间定位信息及语义信息,进而导致网络检测性能下降。故本文对数据集中图片进行裁剪操作,将大分辨图片裁剪为若干尺寸为1 024×1 024的子图。实验时输入图像尺寸为1 024×1 024,由于硬件环境制约,设置batch-size为2,初始学习率设为0.002 5,动量设为0.9,使用随机梯度下降(SGD)优化器。考虑到数据集规模及模型复杂度,为防止模型出现过拟合现象,将训练Epochs设为20,训练过程中,第一个Epoch使用Warmup算法进行预热;输出预测结果后,使用NMS算法进行处理。

3.3 结果分析

图7展示了S2ANet与HQ-S2ANet在DOTA遥感目标验证集的部分检测结果图,图中上半部分为原始图片,下半部分为框选区域局部放大图。图中可看出目标排列密集、尺寸小,变化范围大、变化方向任意的特点。S2ANet会出现漏检或者误检的问题。S2ANet对图7a)左侧尺寸较小的卡车不够敏感,在图7c)中会忽略右侧汽车目标且出现了误检及重复检测的问题,而HQ-S2ANet避免了上述问题的发生。

图7 DOTA数据集上部分检测结果

为验证所提算法的性能,在划分测试集上进行实验,并对比了文中模型以及其他先进模型。结果如表2所示。类别分别为PL:plane,BD:baseball-diamond,BR:bridge,GTF:ground-t-field,SV:small-vehicle,LV:large-vehicle,SH:ship,TC:tennis-court,BC:basketball-court,ST:storage-tank,SBF:soccer-ball-field,RA:roundabout,HA:harbor,SP:swimming-pool,HC:helicopter。

表2 HQ-S2ANet与其他先进模型在DOTA数据集上的性能比较 %

在DOTA数据集上,HQ-S2ANet的mAP为71.3%,与原算法S2ANet精度相比提高3.1%,与Rotated-atss算法相比提高3.3%。Rotated-retinanet精度为65.70%。Rotated-CFA精度为67.30%,该算法虽能较好地检测PL,BR这类分布较为稀疏且目标大小适中的目标,但对于密集目标,如SV,其精度较低。改进算法对SV,HC这类分布密集小目标的AP明显较高,表明其具备解决遥感图像检测难点的能力,能从复杂背景中提取特征信息较少且语义信息较弱的目标。

图8为HQ-S2ANet与S2ANet特征融合后的热力图结果,图8a)～8c)为S2ANet输出结果, 图8d～8f)为HQ-S2ANet输出结果。图8b),8e)为FPN与HQFPN最底层特征图特征可视化,包含丰富细节定位信息。图8e)中经过SEA加权后的特征对小目标给予更高关注度,与图8b)相比,HQ-S2ANet关注区域更广、对目标关注程度更高。图8f)为图8e)经过HQFPN下采样模块细化特征后的热力图。较图8c)经过传统下采样处理,图8f)扩大感受野的同时捕获邻近目标间关系,目标边缘特征更加明确。

图8 S2ANet与HQ-S2ANet特征可视化

图9 锚框阈值与mAP及部分类别召回率关系

3.4 消融实验

将HQFPN,MaxIoUAssigner-HQ、协同注意力模块依次嵌入到S2ANet模型,为验证模块有效性,结果如表3所示。

表3 DOTA数据集上的消融实验研究

嵌入各个模块后,均能获得比原始网络更高的精度。其中嵌入HQFPN模块后,由于特征图经过增强,通过大感受野的方式融合了低层细节位置信息以及高层语义信息,mAP提升0.7%。再进一步优化锚框匹配策略,为抑制锚框低质量匹配的发生,嵌入MaxIoUAssigner-HQ模块后,mAP提升0.8%,证明该模块在保证召回率的同时提升了网络识别精度。为捕获特征中的上下文信息并获取通道间依赖关系,将SEA模块加入网络。特征图经过注意力加权后,更加注重有效信息的提取,mAP提升1.6%,证明网络完成局部间信息交互的同时捕获了通道间信息依赖,能够显著提升目标检测精度。

3.4.1 注意力模块对比实验

该实验在DOTA数据集上对2种注意力融合的先后顺序进行了测试,第一种方式为利用通道注意力对自注意力中Q,K,V向量进行加权后计算关联度矩阵,第二种方式为利用通道注意力对自注意力机制得到的向量评分矩阵进行加权处理后与向量S相加。由于目标种类过多,表中只取尺寸较小及分布较为密集的2类目标,如表4所示。

表4 不同注意力融合方式的HQ-S2ANet性能比较 %

若过早地使用通道注意力对自注意力机制中Q,K,V矩阵进行加权,由于通道间权重经过softmax处理,值域为0～1之间,原始矩阵在与其相乘后,一定程度上弱化了特征。自注意力机制在后续处理过程中,多次运用到Q,K,V进行计算,导致特征的弱化影响被不同程度放大了。故使用Self-ECA构造协同注意力结构,对网络的预测精度影响较小。

3.4.2 下采样方式对比实验

表5对比了在DOTA数据集上3种下采样方式,方式1在低层使用空洞系数为2的空洞卷积,高层使用普通卷积;方式2下采样全过程使用空洞系数为2的空洞卷积;方式3在低层使用空洞系数为2的空洞卷积及均值池化对特征进行提取后将2种特征图相加传入下一层,高层则使用普通卷积。

表5 不同下采样方式的HQ-S2ANet性能比较 %

根据结果分析,下采样方式1与方式2均因空洞卷积的使用,导致了网格隔影的产生,在一定程度上弱化了特征,这对以小目标为主的遥感图像检测影响较大。方式3利用均值池化与空洞卷积相加的方式,在获得大感受野的同时弥补了空洞卷积对特征的弱化,从而提升了网络性能。

3.4.3 锚框阈值对比实验

表6为HQ-S2ANet在控制变量的基础上,在DOTA数据集中对比不同阈值下锚框匹配的检测效果。实验采用的模型为HQ-S2ANet,将阈值设定为0.4,0.6,0.8进行实验。从实验结果可以看出,当阈值为0.8时,部分高质量匹配仍会被负样本匹配重写,由于未能有效抑制锚框低质量匹配发生,导致召回率虽优于阈值为0.4的对照组,但模型精度仍不理想。阈值为0.4时,部分正样本锚框被忽略,导致召回率过低影响网络精度。根据结果判断,最优阈值处于0.4～0.8区间,故设定阈值为0.5,0.7,进一步实验判断。阈值设定为0.5,0.7时,模型精度高于0.4,0.8且低于0.6,故判断最优阈值为0.6。基于模型精度的考虑,设定阈值为0.6。

表6 不同锚框阈值的HQ-S2ANet性能比较 %

3.4.4 算法复杂度分析

表7表明,基于S2ANet与改进的HQ-S2ANet算法相比,浮点型运算量仅增加6.08G,参数量增加2.61M,识别精度更高,较S2ANet提升3.1%。与当前遥感图像目标检测领域主流算法Rotated-retinanet及Rotated-reppoints相比,虽参数量有所提升,但精度提升明显。由此证明:本文所提改进算法HQ-S2ANet在计算开销增加较少的前提下, 大幅提升了算法识别能力。

表7 模型算法复杂度与mAP统计

4 结 语

遥感图像的目标检测在多个领域具有重要研究意义,为解决遥感图像目标成像小、分布拥挤、背景复杂的问题,提出一种基于感知延伸与锚框最适匹配的遥感图像目标检测算法HQ-S2ANet。在特征融合下采样过程中将感知延伸卷积模块与常规卷积交替堆叠构建HQFPN模块,减少下采样过程中低层细节位置信息的丢失以增强网络对小目标的识别能力;为解决遥感目标图像分布拥挤的问题,在为真实边界框分配锚框时,提出高质量锚框匹配算法,在保证召回率的同时,减少低质量锚框匹配的产生;同时设计了协同注意力模块,以较低的计算开销增强了模型对目标与全局特征的建模能力,进而提升模型识别能力。通过在DOTA数据集上的实验,以及与先进模型的对比表明,所提算法可以提升在遥感图像中的目标检测能力,验证了所提模型的有效性。锚框匹配算法中,限定阈值对不同数据可适应性不足,下一步考虑网络通过学习得到最适阈值,在不影响精度的情况下提升方法的鲁棒性。