devbox
编程谜题解决者
这个屌丝很懒，什么也没留下！
关注作者
热门标签
热门文章
当前位置:   article > 正文

yoloe 超过yolox和yolov5

作者：编程谜题解决者 | 2024-01-31 20:45:43

赞

踩

yoloe

号称超过yolox和yolov5的PP-YoLoE

开源地址：

GitHub - PaddlePaddle/PaddleDetection: Object Detection toolkit based on PaddlePaddle. It supports object detection, instance segmentation, multiple object tracking and real-time multi-person keypoint detection.

使用笔记：

yoloe 目标检测使用笔记_AI视觉网奇的博客-CSDN博客

PP-YoLoE | PP-YoLov2全面升级Anchor-Free，速度精度完美超越YoLoX和YoLov5

PP-YOLOE是基于PP-YOLOv2的卓越的单阶段Anchor-free模型，超越了多种流行的yolo模型。PP-YOLOE有一系列的模型，即s/m/l/x，可以通过width multiplier和depth multiplier配置。PP-YOLOE避免使用诸如deformable convolution或者matrix nms之类的特殊算子，以使其能轻松地部署在多种多样的硬件上。

1模型架构

PP-YOLOE由以下方法组成：

  • 可扩展的backbone和neck

  • Task Alignment Learning

  • Efficient Task-aligned head with DFL和VFL

  • SiLU激活函数

1.1、Backbone

PP-YOLOE的Backbone主要是使用RepVGG模块以及CSP的模型思想对ResNet及逆行的改进，同时也使用了SiLU激活函数、Effitive SE Attention等模块，下面我们一一道来。

1、RepVGG

RepVGG，这个网络就是在VGG的基础上面进行改进，主要的思路包括：

  1. 在VGG网络的Block块中加入了Identity和残差分支，相当于把ResNet网络中的精华应用 到VGG网络中；

  2. 模型推理阶段，通过Op融合策略将所有的网络层都转换为3×3卷积，便于网络的部署和加速。

上图展示了模型推理阶段的重参数化过程，其实就是一个OP融合和OP替换的过程。图A从结构化的角度展示了整个重参数化流程， 图B从模型参数的角度展示了整个重参数化流程。整个重参数化步骤如下所示：

步骤1：首先通过式3将残差块中的卷积层和BN层进行融合，该操作在很多深度学习框架的推理阶段都会执行。图中的蓝色框中执行3×3卷积+BN层的融合，图中的黑色矩形框中执行1×1卷积+BN层的融合，图中的黄色矩形框中执行3×3卷积(卷积核设置为全1)+BN层的融合。其中表示转换前的卷积层参数，表示BN层的均值，表示BN层的方差，和分别表示BN层的尺度因子和偏移因子，W’和b’分别表示融合之后的卷积的权重和偏置。

步骤2：将融合后的卷积层转换为3×3卷积，即将具体不同卷积核的卷积均转换为具有3×3大小的卷积核的卷积。由于整个残差块中可能包含1×1卷积分支和Identity两种分支，如图中的黑框和黄框所示。对于1×1卷积分支而言，整个转换过程就是利用3×3卷积核替换1×1卷积核，具体的细节如图中的紫框所示，即将1×1卷积核中的数值移动到3×3卷积核的中心点即可；对于Identity分支而言，该分支并没有改变输入的特征映射的数值，那么可以设置一个3×3的卷积核，将所有的9个位置处的权重值都设置为1，那么它与输入的特征映射相乘之后，保持了原来的数值，具体的细节如图中的褐色框所示。

步骤3：合并残差分支中的3×3卷积。即将所有分支的权重W和偏置B叠加起来，从而获得一个融合之后的3×3卷积层。

为什么要用VGG式模型？

除了相信简单就是美以外，VGG式极简模型至少还有5大现实的优势：

  1. 3×3卷积非常快。在GPU上，3×3卷积的计算密度（理论运算量除以所用时间）可达1×1和5×5卷积的4倍。

  2. 单路架构非常快，因为并行度高。同样的计算量，“大而整”的运算效率远超“小而碎”的运算。

  3. 单路架构省内存。例如，ResNet的shortcut虽然不占计算量，却增加了一倍的显存占用。

  4. 单路架构灵活性更好，容易改变各层的宽度（如剪枝）。

  5. RepVGG主体部分只有一种算子：3×3卷积接ReLU。在设计专用芯片时，给定芯片尺寸或造价可以集成海量的3×3卷积+ReLU计算单元来达到很高的效率。

下图表示RepVGG推理融合后的ONNX输出，可以看出简化了很多。

  1. class RepVggBlock(nn.Layer):
  2.     def __init__(self, ch_in, ch_out, act='relu'):
  3.         super(RepVggBlock, self).__init__()
  4.         self.ch_in = ch_in
  5.         self.ch_out = ch_out
  6.         self.conv1 = ConvBNLayer(ch_in, ch_out, 3, stride=1, padding=1, act=None)
  7.         self.conv2 = ConvBNLayer(ch_in, ch_out, 1, stride=1, padding=0, act=None)
  8.         self.act = get_act_fn(act) if act is None or isinstance(act, (str, dict)) else act
  9.  
  10.     def forward(self, x):
  11.         if hasattr(self, 'conv'):
  12.             y = self.conv(x)
  13.         else:
  14.             y = self.conv1(x) + self.conv2(x)
  15.         y = self.act(y)
  16.         return y
  17.  
  18.     def convert_to_deploy(self):
  19.         if not hasattr(self, 'conv'):
  20.             self.conv = nn.Conv2D(in_channels=self.ch_in, out_channels=self.ch_out, kernel_size=3, stride=1, padding=1, groups=1)
  21.         kernel, bias = self.get_equivalent_kernel_bias()
  22.         self.conv.weight.set_value(kernel)
  23.         self.conv.bias.set_value(bias)
  24.  
  25.     def get_equivalent_kernel_bias(self):
  26.         # 融合推理
  27.         kernel3x3, bias3x3 = self._fuse_bn_tensor(self.conv1)
  28.         kernel1x1, bias1x1 = self._fuse_bn_tensor(self.conv2)
  29.         return kernel3x3 + self._pad_1x1_to_3x3_tensor(kernel1x1), bias3x3 + bias1x1
  30.  
  31.     def _pad_1x1_to_3x3_tensor(self, kernel1x1):
  32.         if kernel1x1 is None:
  33.             return 0
  34.         else:
  35.             return nn.functional.pad(kernel1x1, [1, 1, 1, 1])
  36.  
  37.     def _fuse_bn_tensor(self, branch):
  38.         if branch is None:
  39.             return 0, 0
  40.         kernel = branch.conv.weight
  41.         running_mean = branch.bn._mean
  42.         running_var = branch.bn._variance
  43.         gamma = branch.bn.weight
  44.         beta = branch.bn.bias
  45.         eps = branch.bn._epsilon
  46.         std = (running_var + eps).sqrt()
  47.         t = (gamma / std).reshape((-1, 1, 1, 1))
  48.         return kernel * t, beta - running_mean * gamma / std

2、Swish激活函数

从代码和公式来看，Swish包含了SiLU，换句话说SiLU是Swish的一种特例。

所以画图基本上都使用了SiLU代替Swish，因为YOLOE中的Swish的，也就是SiLU激活函数。

β是个常数或可训练的参数。Swish 具备无上界有下界、平滑、非单调的特性。Swish 在深层模型上的效果优于 ReLU。

例如，仅仅使用 Swish 单元替换 ReLU 就能把 Mobile NASNetA 在 ImageNet 上的 top-1 分类准确率提高 0.9%，Inception-ResNet-v 的分类准确率提高 0.6%。

  1. class ConvBNLayer(nn.Layer):
  2.     def __init__(self, ch_in, ch_out, filter_size=3, stride=1, groups=1, padding=0, act=None):
  3.         super(ConvBNLayer, self).__init__()
  4.         self.conv = nn.Conv2D(in_channels=ch_in, out_channels=ch_out, kernel_size=filter_size, stride=stride, padding=padding, groups=groups, bias_attr=False)
  5.         self.bn = nn.BatchNorm2D(ch_out, weight_attr=ParamAttr(regularizer=L2Decay(0.0)), bias_attr=ParamAttr(regularizer=L2Decay(0.0)))
  6.         self.act = get_act_fn(act) if act is None or isinstance(act, (str, dict)) else act
  7.  
  8.     def forward(self, x):
  9.         x = self.conv(x)
  10.         x = self.bn(x)
  11.         x = self.act(x)
  12.  
  13.         return x

3、Effective SE Attention

该模块主要是来自于《CenterMask：Real-Time Anchor-Free Instance Segmentation》中的eSE模块；

在输出的内部添加了一个channel上的attention模块eSE。原始的SE模块中使用2个FC去进行channel权重映射，但是为了减少计算量通常会将FC中的channel给剪裁一些（小于输入的channel），这就引入了一些信息的损失，为此文章直接将2个FC替换为了1个FC。

  1. class EffectiveSELayer(nn.Layer):
  2.     def __init__(self, channels, act='hardsigmoid'):
  3.         super(EffectiveSELayer, self).__init__()
  4.         self.fc = nn.Conv2D(channels, channels, kernel_size=1, padding=0)
  5.         self.act = get_act_fn(act) if act is None or isinstance(act, (str, dict)) else act
  6.  
  7.     def forward(self, x):
  8.         x_se = x.mean((2, 3), keepdim=True)
  9.         x_se = self.fc(x_se)
  10.         return x * self.act(x_se)

4、CSPNet结构

CSPNet的主要思想还是Partial Dense Block，设计Partial Dense Block的目的是：

  • 增加梯度路径：通过分裂合并策略，可以使梯度路径的数目翻倍。由于采用了跨阶段的策略，可以减轻使用显式特征映射复制进行连接的缺点；

  • 平衡各层的计算：通常情况下，DenseNet底层的信道数远远大于增长率。由于部分dense block中的dense layer操作所涉及的底层信道只占原始信道的一半，因此可以有效地解决近一半的计算瓶颈；

  • 减少内存流量：假设dense block在DenseNet中的基本特征映射大小为，增长率为d，并且共m层。然后，该dense block的CIO为，部分dense block的CIO为。虽然m和d通常比c小得多，但部分dense block最多可以节省网络内存流量的一半。

  1. class CSPResStage(nn.Layer):
  2.     def __init__(self, block_fn, ch_in, ch_out, n, stride, act='relu', attn='eca'):
  3.         super(CSPResStage, self).__init__()
  4.         ch_mid = (ch_in + ch_out) // 2
  5.         if stride == 2:
  6.             self.conv_down = ConvBNLayer(ch_in, ch_mid, 3, stride=2, padding=1, act=act)
  7.         else:
  8.             self.conv_down = None
  9.         self.conv1 = ConvBNLayer(ch_mid, ch_mid // 2, 1, act=act)
  10.         self.conv2 = ConvBNLayer(ch_mid, ch_mid // 2, 1, act=act)
  11.         self.blocks = nn.Sequential(* [block_fn(ch_mid // 2, ch_mid // 2, act=act, shortcut=True) for i in range(n)])
  12.         if attn:
  13.             self.attn = EffectiveSELayer(ch_mid, act='hardsigmoid')
  14.         else:
  15.             self.attn = None
  16.  
  17.         self.conv3 = ConvBNLayer(ch_mid, ch_out, 1, act=act)
  18.  
  19.     def forward(self, x):
  20.         if self.conv_down is not None:
  21.             x = self.conv_down(x)
  22.         y1 = self.conv1(x)
  23.         y2 = self.blocks(self.conv2(x))
  24.         y = paddle.concat([y1, y2], axis=1)
  25.         if self.attn is not None:
  26.             y = self.attn(y)
  27.         y = self.conv3(y)
  28.         return y

5、SPP结构

SPP-Net全名为Spatial Pyramid Pooling（空间金字塔池化结构），2015年由微软研究院的何恺明提出,主要解决2个问题：

  1. 有效避免了R-CNN算法对图像区域剪裁、缩放操作导致的图像物体剪裁不全以及形状扭曲等问题。

  2. 解决了卷积神经网络对图像重复特征提取的问题，大大提高了产生候选框的速度，且节省了计算成本。

SPP 显著特点

  1. 不管输入尺寸是怎样，SPP 可以产生固定大小的输出

  2. 使用多个窗口(pooling window)

  3. SPP 可以使用同一图像不同尺寸(scale)作为输入, 得到同样长度的池化特征。

其它特点

  1. 由于对输入图像的不同纵横比和不同尺寸，SPP同样可以处理，所以提高了图像的尺度不变(scale-invariance)和降低了过拟合(over-fitting)

  2. 实验表明训练图像尺寸的多样性比单一尺寸的训练图像更容易使得网络收敛(convergence)

  3. SPP 对于特定的CNN网络设计和结构是独立的。(也就是说，只要把SPP放在最后一层卷积层后面，对网络的结构是没有影响的， 它只是替换了原来的pooling层)

  4. 不仅可以用于图像分类而且可以用来目标检测

通过spp模块实现局部特征和全局特征（所以空间金字塔池化结构的最大的池化核要尽可能的接近等于需要池化的featherMap的大小）的featherMap级别的融合，丰富最终特征图的表达能力，从而提高MAP。

  1. class SPP(nn.Layer):
  2.     def __init__(self, ch_in, ch_out, k, pool_size, act='swish', data_format='NCHW'):
  3.         super(SPP, self).__init__()
  4.         self.pool = []
  5.         self.data_format = data_format
  6.         for i, size in enumerate(pool_size):
  7.             pool = self.add_sublayer('pool{}'.format(i), nn.MaxPool2D(kernel_size=size, stride=1, padding=size // 2, data_format=data_format, ceil_mode=False))
  8.             self.pool.append(pool)
  9.         self.conv = ConvBNLayer(ch_in, ch_out, k, padding=k // 2, act=act)
  10.  
  11.     def forward(self, x):
  12.         outs = [x]
  13.         for pool in self.pool:
  14.             outs.append(pool(x))
  15.         if self.data_format == 'NCHW':
  16.             y = paddle.concat(outs, axis=1)
  17.         else:
  18.             y = paddle.concat(outs, axis=-1)
  19.  
  20.         y = self.conv(y)
  21.         return y
  22.  
  23.  
  24. class CSPStage(nn.Layer):
  25.     def __init__(self, block_fn, ch_in, ch_out, n, act='swish', spp=False):
  26.         super(CSPStage, self).__init__()
  27.  
  28.         ch_mid = int(ch_out // 2)
  29.         self.conv1 = ConvBNLayer(ch_in, ch_mid, 1, act=act)
  30.         self.conv2 = ConvBNLayer(ch_in, ch_mid, 1, act=act)
  31.         self.convs = nn.Sequential()
  32.         next_ch_in = ch_mid
  33.         for i in range(n):
  34.             self.convs.add_sublayer(str(i), eval(block_fn)(next_ch_in, ch_mid, act=act, shortcut=False))
  35.             if i == (n - 1) // 2 and spp:
  36.                 self.convs.add_sublayer('spp', SPP(ch_mid * 4, ch_mid, 1, [5, 9, 13], act=act))
  37.             next_ch_in = ch_mid
  38.         self.conv3 = ConvBNLayer(ch_mid * 2, ch_out, 1, act=act)
  39.  
  40.     def forward(self, x):
  41.         y1 = self.conv1(x)
  42.         y2 = self.conv2(x)
  43.         y2 = self.convs(y2)
  44.         y = paddle.concat([y1, y2], axis=1)
  45.         y = self.conv3(y)
  46.         return y

1.2、Neck

yoloe的neck结构采用的依旧是FPN+PAN结构模式，将Neck部分用立体图画出来，更直观的看下两部分之间是如何通过FPN结构融合的。

如图所示，FPN是自顶向下的，将高层特征通过上采样和低层特征做融合得到进行预测的特征图。

和FPN层不同，yoloe在FPN层的后面还添加了一个自底向上的特征金字塔。FPN是自顶向下，将高层的强语义特征传递下来，对整个金字塔进行增强，不过只增强了语义信息，对定位信息没有传递，而本文就是针对这一点，在FPN的后面添加一个自底向上的金字塔。这样的操作是对FPN的补充，将低层的强定位特征传递上去。

  1. class CustomCSPPAN(nn.Layer):
  2.     __shared__ = ['norm_type', 'data_format', 'width_mult', 'depth_mult', 'trt']
  3.  
  4.     def __init__(self, in_channels=[256, 512, 1024], out_channels=[1024, 512, 256], norm_type='bn', act='leaky',
  5.                  stage_fn='CSPStage', block_fn='BasicBlock', stage_num=1, block_num=3, drop_block=False,
  6.                  block_size=3, keep_prob=0.9, spp=False, data_format='NCHW', width_mult=1.0,
  7.                  depth_mult=1.0, trt=False):
  8.  
  9.         super(CustomCSPPAN, self).__init__()
  10.         out_channels = [max(round(c * width_mult), 1) for c in out_channels]
  11.         block_num = max(round(block_num * depth_mult), 1)
  12.         act = get_act_fn(act, trt=trt) if act is None or isinstance(act, (str, dict)) else act
  13.         self.num_blocks = len(in_channels)
  14.         self.data_format = data_format
  15.         self._out_channels = out_channels
  16.         in_channels = in_channels[::-1]
  17.         fpn_stages = []
  18.         fpn_routes = []
  19.         for i, (ch_in, ch_out) in enumerate(zip(in_channels, out_channels)):
  20.             if i > 0:
  21.                 ch_in += ch_pre // 2
  22.  
  23.             stage = nn.Sequential()
  24.             for j in range(stage_num):
  25.                 stage.add_sublayer(str(j), eval(stage_fn)(block_fn, ch_in if j == 0 else ch_out, ch_out, block_num, act=act, spp=(spp and i == 0)))
  26.  
  27.             if drop_block:
  28.                 stage.add_sublayer('drop', DropBlock(block_size, keep_prob))
  29.  
  30.             fpn_stages.append(stage)
  31.  
  32.             if i < self.num_blocks - 1:
  33.                 fpn_routes.append(ConvBNLayer(ch_in=ch_out, ch_out=ch_out // 2, filter_size=1, stride=1, padding=0, act=act))
  34.  
  35.             ch_pre = ch_out
  36.  
  37.         self.fpn_stages = nn.LayerList(fpn_stages)
  38.         self.fpn_routes = nn.LayerList(fpn_routes)
  39.  
  40.         pan_stages = []
  41.         pan_routes = []
  42.         for i in reversed(range(self.num_blocks - 1)):
  43.             pan_routes.append(ConvBNLayer(ch_in=out_channels[i + 1], ch_out=out_channels[i + 1], filter_size=3, stride=2, padding=1, act=act))
  44.             ch_in = out_channels[i] + out_channels[i + 1]
  45.             ch_out = out_channels[i]
  46.             stage = nn.Sequential()
  47.             for j in range(stage_num):
  48.                 stage.add_sublayer(str(j), eval(stage_fn)(block_fn, ch_in if j == 0 else ch_out, ch_out, block_num, act=act, spp=False))
  49.             if drop_block:
  50.                 stage.add_sublayer('drop', DropBlock(block_size, keep_prob))
  51.  
  52.             pan_stages.append(stage)
  53.  
  54.         self.pan_stages = nn.LayerList(pan_stages[::-1])
  55.         self.pan_routes = nn.LayerList(pan_routes[::-1])
  56.  
  57.     def forward(self, blocks, for_mot=False):
  58.         blocks = blocks[::-1]
  59.         fpn_feats = []
  60.  
  61.         for i, block in enumerate(blocks):
  62.             if i > 0:
  63.                 block = paddle.concat([route, block], axis=1)
  64.             route = self.fpn_stages[i](block)
  65.             fpn_feats.append(route)
  66.  
  67.             if i < self.num_blocks - 1:
  68.                 route = self.fpn_routes[i](route)
  69.                 route = F.interpolate(route, scale_factor=2., data_format=self.data_format)
  70.  
  71.         pan_feats = [fpn_feats[-1], ]
  72.         route = fpn_feats[-1]
  73.         for i in reversed(range(self.num_blocks - 1)):
  74.             block = fpn_feats[i]
  75.             route = self.pan_routes[i](route)
  76.             block = paddle.concat([route, block], axis=1)
  77.             route = self.pan_stages[i](block)
  78.             pan_feats.append(route)
  79.  
  80.         return pan_feats[::-1]

1.3、Head

对于PP-YOLOE的head部分，其依旧是TOOD的head，也就是T-Head，主要是包括了Cls Head和Loc Head。具体来说，T-head首先在FPN特征基础上进行分类与定位预测；然后TAL基于所提任务对齐测度计算任务对齐信息；最后T-head根据从TAL传回的信息自动调整分类概率与定位预测。

由于2个任务的预测都是基于这个交互特征来完成的，但是2个任务对于特征的需求肯定是不一样的，因为作者设计了一个layer attention来为每个任务单独的调整一下特征，这个部分的结构也很简单，可以理解为是一个channel-wise的注意力机制。这样的话就得到了对于每个任务单独的特征,然后再利用这些特征生成所需要的类别或者定位的特征图。

  1. class PPYOLOEHead(nn.Layer):
  2.     __shared__ = ['num_classes', 'trt', 'exclude_nms']
  3.     __inject__ = ['static_assigner', 'assigner', 'nms']
  4.  
  5.     def __init__(self,
  6.                  in_channels=[1024, 512, 256],
  7.                  num_classes=80,
  8.                  act='swish',
  9.                  fpn_strides=(32, 16, 8),
  10.                  grid_cell_scale=5.0,
  11.                  grid_cell_offset=0.5,
  12.                  reg_max=16,
  13.                  static_assigner_epoch=4,
  14.                  use_varifocal_loss=True,
  15.                  static_assigner='ATSSAssigner',
  16.                  assigner='TaskAlignedAssigner',
  17.                  nms='MultiClassNMS',
  18.                  eval_input_size=[],
  19.                  loss_weight={'class': 1.0, 'iou': 2.5, 'dfl': 0.5,},
  20.                  trt=False,
  21.                  exclude_nms=False):
  22.         super(PPYOLOEHead, self).__init__()
  23.         assert len(in_channels) > 0, "len(in_channels) should > 0"
  24.         self.in_channels = in_channels
  25.         self.num_classes = num_classes
  26.         self.fpn_strides = fpn_strides
  27.         self.grid_cell_scale = grid_cell_scale
  28.         self.grid_cell_offset = grid_cell_offset
  29.         self.reg_max = reg_max
  30.         self.iou_loss = GIoULoss()
  31.         self.loss_weight = loss_weight
  32.         self.use_varifocal_loss = use_varifocal_loss
  33.         self.eval_input_size = eval_input_size
  34.  
  35.         self.static_assigner_epoch = static_assigner_epoch
  36.         self.static_assigner = static_assigner
  37.         self.assigner = assigner
  38.         self.nms = nms
  39.         self.exclude_nms = exclude_nms
  40.         # stem
  41.         self.stem_cls = nn.LayerList()
  42.         self.stem_reg = nn.LayerList()
  43.         act = get_act_fn(act, trt=trt) if act is None or isinstance(act, (str, dict)) else act
  44.         for in_c in self.in_channels:
  45.             self.stem_cls.append(ESEAttn(in_c, act=act))
  46.             self.stem_reg.append(ESEAttn(in_c, act=act))
  47.         # pred head
  48.         self.pred_cls = nn.LayerList()
  49.         self.pred_reg = nn.LayerList()
  50.         for in_c in self.in_channels:
  51.             self.pred_cls.append(nn.Conv2D(in_c, self.num_classes, 3, padding=1))
  52.             self.pred_reg.append(nn.Conv2D(in_c, 4 * (self.reg_max + 1), 3, padding=1))
  53.         # projection conv
  54.         self.proj_conv = nn.Conv2D(self.reg_max + 1, 1, 1, bias_attr=False)
  55.         self._init_weights()
  56.  
  57.     @classmethod
  58.     def from_config(cls, cfg, input_shape):
  59.         return {'in_channels': [i.channels for i in input_shape], }
  60.  
  61.     def forward_train(self, feats, targets):
  62.         anchors, anchor_points, num_anchors_list, stride_tensor = generate_anchors_for_grid_cell(feats, self.fpn_strides, self.grid_cell_scale, self.grid_cell_offset)
  63.  
  64.         cls_score_list, reg_distri_list = [], []
  65.         for i, feat in enumerate(feats):
  66.             avg_feat = F.adaptive_avg_pool2d(feat, (1, 1))
  67.             cls_logit = self.pred_cls[i](self.stem_cls[i](feat, avg_feat) + feat)
  68.             reg_distri = self.pred_reg[i](self.stem_reg[i](feat, avg_feat))
  69.             # cls and reg
  70.             cls_score = F.sigmoid(cls_logit)
  71.             cls_score_list.append(cls_score.flatten(2).transpose([0, 2, 1]))
  72.             reg_distri_list.append(reg_distri.flatten(2).transpose([0, 2, 1]))
  73.         cls_score_list = paddle.concat(cls_score_list, axis=1)
  74.         reg_distri_list = paddle.concat(reg_distri_list, axis=1)
  75.  
  76.         return self.get_loss([cls_score_list, reg_distri_list, anchors, anchor_points, num_anchors_list, stride_tensor], targets)
  77.  
  78.     def forward_eval(self, feats):
  79.         if self.eval_input_size:
  80.             anchor_points, stride_tensor = self.anchor_points, self.stride_tensor
  81.         else:
  82.             anchor_points, stride_tensor = self._generate_anchors(feats)
  83.         cls_score_list, reg_dist_list = [], []
  84.         for i, feat in enumerate(feats):
  85.             b, _, h, w = feat.shape
  86.             l = h * w
  87.             avg_feat = F.adaptive_avg_pool2d(feat, (1, 1))
  88.             cls_logit = self.pred_cls[i](self.stem_cls[i](feat, avg_feat) + feat)
  89.             reg_dist = self.pred_reg[i](self.stem_reg[i](feat, avg_feat))
  90.             reg_dist = reg_dist.reshape([-1, 4, self.reg_max + 1, l]).transpose([0, 2, 1, 3])
  91.             reg_dist = self.proj_conv(F.softmax(reg_dist, axis=1))
  92.             # cls and reg
  93.             cls_score = F.sigmoid(cls_logit)
  94.             cls_score_list.append(cls_score.reshape([b, self.num_classes, l]))
  95.             reg_dist_list.append(reg_dist.reshape([b, 4, l]))
  96.  
  97.         cls_score_list = paddle.concat(cls_score_list, axis=-1)
  98.         reg_dist_list = paddle.concat(reg_dist_list, axis=-1)
  99.  
  100.         return cls_score_list, reg_dist_list, anchor_points, stride_tensor
  101.  
  102.     def forward(self, feats, targets=None):
  103.         assert len(feats) == len(self.fpn_strides), "The size of feats is not equal to size of fpn_strides"
  104.  
  105.         if self.training:
  106.             return self.forward_train(feats, targets)
  107.         else:
  108.             return self.forward_eval(feats)

2样本匹配

2.1 ATSS Assigner思想

ATSS论文指出One-Stage Anchor-Based和Center-Based Anchor-Free检测算法间的差异主要来自于正负样本的选择，基于此提出ATSS(Adaptive Training Sample Selection)方法，该方法能够自动根据GT的相关统计特征选择合适的Anchor Box作为正样本，在不带来额外计算量和参数的情况下，能够大幅提升模型的性能。

ATSS选取正样本的方法如下：其简要流程为：

  1. 计算每个 gt bbox 和多尺度输出层的所有 anchor 之间的 IoU

  2. 计算每个 gt bbox 中心坐标和多尺度输出层的所有 anchor 中心坐标的 l2 距离

  3. 遍历每个输出层，遍历每个 gt bbox，找出当前层中 topk (超参，默认是 9 )个最小 l2 距离的 anchor 。假设一共有 l 个输出层，那么对于任何一个 gt bbox，都会挑选出 topk×l 个候选位置

  4. 对于每个 gt bbox，计算所有候选位置 IoU 的均值和标准差，两者相加得到该 gt bbox 的自适应阈值

  5. 遍历每个 gt bbox，选择出候选位置中 IoU 大于阈值的位置，该位置认为是正样本，负责预测该 gt bbox

  6. 如果 topk 参数设置过大，可能会导致某些正样本位置不在 gt bbox 内部，故需要过滤掉这部分正样本，设置为背景样本

1、ATSS主要2大特性：

  1. 保证了所有的正样本Anchor都是在Ground Truth的周围。

  2. 最主要是根据不同层的特性对不同层的正样本的阈值进行了微调。

2、ATSS的贡献

  • 指出Anchor-Base检测器和Anchor-Free检测器之间的本质区别实际上是如何定义正训练样本和负训练样本；

  • 提出自适应训练样本选择，以根据目标的统计特征自动选择正负样本；

  • 证明了在图像上的每个位置上平铺多个Anchor来提升检测的性能是没效果的；

  1. class ATSSAssigner(nn.Layer):
  2.     """Bridging the Gap Between Anchor-based and Anchor-free Detection
  3.      via Adaptive Training Sample Selection
  4.     """
  5.     __shared__ = ['num_classes']
  6.  
  7.     def __init__(self, topk=9, num_classes=80, force_gt_matching=False, eps=1e-9):
  8.         super(ATSSAssigner, self).__init__()
  9.         self.topk = topk
  10.         self.num_classes = num_classes
  11.         self.force_gt_matching = force_gt_matching
  12.         self.eps = eps
  13.  
  14.     def _gather_topk_pyramid(self, gt2anchor_distances, num_anchors_list, pad_gt_mask):
  15.         pad_gt_mask = pad_gt_mask.tile([1, 1, self.topk]).astype(paddle.bool)
  16.         gt2anchor_distances_list = paddle.split(gt2anchor_distances, num_anchors_list, axis=-1)
  17.         num_anchors_index = np.cumsum(num_anchors_list).tolist()
  18.         num_anchors_index = [0, ] + num_anchors_index[:-1]
  19.         is_in_topk_list = []
  20.         topk_idxs_list = []
  21.         for distances, anchors_index in zip(gt2anchor_distances_list, num_anchors_index):
  22.             num_anchors = distances.shape[-1]
  23.             topk_metrics, topk_idxs = paddle.topk(distances, self.topk, axis=-1, largest=False)
  24.             topk_idxs_list.append(topk_idxs + anchors_index)
  25.             topk_idxs = paddle.where(pad_gt_mask, topk_idxs, paddle.zeros_like(topk_idxs))
  26.             is_in_topk = F.one_hot(topk_idxs, num_anchors).sum(axis=-2)
  27.             is_in_topk = paddle.where(is_in_topk > 1, paddle.zeros_like(is_in_topk), is_in_topk)
  28.             is_in_topk_list.append(is_in_topk.astype(gt2anchor_distances.dtype))
  29.         is_in_topk_list = paddle.concat(is_in_topk_list, axis=-1)
  30.         topk_idxs_list = paddle.concat(topk_idxs_list, axis=-1)
  31.         return is_in_topk_list, topk_idxs_list
  32.  
  33.     @paddle.no_grad()
  34.     def forward(self, anchor_bboxes, num_anchors_list, gt_labels, gt_bboxes, pad_gt_mask, bg_index, gt_scores=None, pred_bboxes=None):
  35.         """
  36.         ATSS匹配步骤如下：
  37.         1. 计算所有预测bbox与GT之间的IoU
  38.         2. 计算所有预测bbox与GT之间的距离
  39.         3. 在每个pyramid level上，对于每个gt，选择k个中心距离gt中心最近的bbox，总共选择k*l个bbox作为每个gt的候选框
  40.         4. 获取这些候选框对应的iou，计算mean和std，设 mean + std为 iou 阈值
  41.         5. 选择iou大于或等于阈值的样本为正样本
  42.         6. 将正样本的中心限制在gt内
  43.         7. 如果Anchor框被分配到多个gts，则选择具有最高的IoU的那个。
  44.         Args:
  45.             anchor_bboxes (Tensor, float32): pre-defined anchors, shape(L, 4),
  46.                     "xmin, xmax, ymin, ymax" format
  47.             num_anchors_list (List): num of anchors in each level
  48.             gt_labels (Tensor, int64|int32): Label of gt_bboxes, shape(B, n, 1)
  49.             gt_bboxes (Tensor, float32): Ground truth bboxes, shape(B, n, 4)
  50.             pad_gt_mask (Tensor, float32): 1 means bbox, 0 means no bbox, shape(B, n, 1)
  51.             bg_index (int): background index
  52.             gt_scores (Tensor|None, float32) Score of gt_bboxes,
  53.                     shape(B, n, 1), if None, then it will initialize with one_hot label
  54.             pred_bboxes (Tensor, float32, optional): predicted bounding boxes, shape(B, L, 4)
  55.         Returns:
  56.             assigned_labels (Tensor): (B, L)
  57.             assigned_bboxes (Tensor): (B, L, 4)
  58.             assigned_scores (Tensor): (B, L, C), if pred_bboxes is not None, then output ious
  59.         """
  60.         assert gt_labels.ndim == gt_bboxes.ndim and gt_bboxes.ndim == 3
  61.  
  62.         num_anchors, _ = anchor_bboxes.shape
  63.         batch_size, num_max_boxes, _ = gt_bboxes.shape
  64.  
  65.         # 1. 计算所有预测bbox与GT之间的IoU, [B, n, L]
  66.         ious = iou_similarity(gt_bboxes.reshape([-1, 4]), anchor_bboxes)
  67.         ious = ious.reshape([batch_size, -1, num_anchors])
  68.  
  69.         # 2. 计算所有预测bbox与GT之间的距离, [B, n, L]
  70.         gt_centers = bbox_center(gt_bboxes.reshape([-1, 4])).unsqueeze(1)
  71.         anchor_centers = bbox_center(anchor_bboxes)
  72.         gt2anchor_distances = (gt_centers - anchor_centers.unsqueeze(0)).norm(2, axis=-1).reshape([batch_size, -1, num_anchors])
  73.  
  74.         # 3. 在每个pyramid level上，对于每个gt，选择k个中心距离gt中心最近的bbox，总共选择k*l个bbox作为每个gt的候选框
  75.         # based on the center distance, [B, n, L]
  76.         is_in_topk, topk_idxs = self._gather_topk_pyramid(gt2anchor_distances, num_anchors_list, pad_gt_mask)
  77.  
  78.         # 4. 获取这些候选框对应的iou，计算mean和std，设 mean + std为 iou 阈值
  79.         iou_candidates = ious * is_in_topk
  80.         iou_threshold = paddle.index_sample(iou_candidates.flatten(stop_axis=-2), topk_idxs.flatten(stop_axis=-2))
  81.         iou_threshold = iou_threshold.reshape([batch_size, num_max_boxes, -1])
  82.         iou_threshold = iou_threshold.mean(axis=-1, keepdim=True) + iou_threshold.std(axis=-1, keepdim=True)
  83.         is_in_topk = paddle.where(iou_candidates > iou_threshold.tile([1, 1, num_anchors]), is_in_topk, paddle.zeros_like(is_in_topk))
  84.  
  85.         # 6. 将正样本的中心限制在gt内, [B, n, L]
  86.         is_in_gts = check_points_inside_bboxes(anchor_centers, gt_bboxes)
  87.  
  88.         # 选择正样本, [B, n, L]
  89.         mask_positive = is_in_topk * is_in_gts * pad_gt_mask
  90.  
  91.         # 7. 如果Anchor框被分配到多个gts，则选择具有最高的IoU的那个。
  92.         mask_positive_sum = mask_positive.sum(axis=-2)
  93.         if mask_positive_sum.max() > 1:
  94.             mask_multiple_gts = (mask_positive_sum.unsqueeze(1) > 1).tile([1, num_max_boxes, 1])
  95.             is_max_iou = compute_max_iou_anchor(ious)
  96.             mask_positive = paddle.where(mask_multiple_gts, is_max_iou, mask_positive)
  97.             mask_positive_sum = mask_positive.sum(axis=-2)
  98.         # 8. 确认每个gt_bbox 都匹配到了 anchor
  99.         if self.force_gt_matching:
  100.             is_max_iou = compute_max_iou_gt(ious) * pad_gt_mask
  101.             mask_max_iou = (is_max_iou.sum(-2, keepdim=True) == 1).tile([1, num_max_boxes, 1])
  102.             mask_positive = paddle.where(mask_max_iou, is_max_iou, mask_positive)
  103.             mask_positive_sum = mask_positive.sum(axis=-2)
  104.         assigned_gt_index = mask_positive.argmax(axis=-2)
  105.  
  106.         # 匹配目标
  107.         batch_ind = paddle.arange(end=batch_size, dtype=gt_labels.dtype).unsqueeze(-1)
  108.         assigned_gt_index = assigned_gt_index + batch_ind * num_max_boxes
  109.         assigned_labels = paddle.gather(gt_labels.flatten(), assigned_gt_index.flatten(), axis=0)
  110.         assigned_labels = assigned_labels.reshape([batch_size, num_anchors])
  111.         assigned_labels = paddle.where(mask_positive_sum > 0, assigned_labels, paddle.full_like(assigned_labels, bg_index))
  112.  
  113.         assigned_bboxes = paddle.gather(gt_bboxes.reshape([-1, 4]), assigned_gt_index.flatten(), axis=0)
  114.         assigned_bboxes = assigned_bboxes.reshape([batch_size, num_anchors, 4])
  115.  
  116.         assigned_scores = F.one_hot(assigned_labels, self.num_classes)
  117.         if pred_bboxes is not None:
  118.             # assigned iou
  119.             ious = batch_iou_similarity(gt_bboxes, pred_bboxes) * mask_positive
  120.             ious = ious.max(axis=-2).unsqueeze(-1)
  121.             assigned_scores *= ious
  122.         elif gt_scores is not None:
  123.             gather_scores = paddle.gather(gt_scores.flatten(), assigned_gt_index.flatten(), axis=0)
  124.             gather_scores = gather_scores.reshape([batch_size, num_anchors])
  125.             gather_scores = paddle.where(mask_positive_sum > 0, gather_scores, paddle.zeros_like(gather_scores))
  126.             assigned_scores *= gather_scores.unsqueeze(-1)
  127.  
  128.         return assigned_labels, assigned_bboxes, assigned_scores

2.2、Task-aligned Assigner思想（TOOD）

TOOD提出了Task Alignment Learning (TAL) 来显式的把2个任务的最优Anchor拉近。这是通过设计一个样本分配策略和任务对齐loss来实现的。样本分配计算每个Anchor的任务对齐度，同时任务对齐loss可以逐步将分类和定位的最佳Anchor统一起来。

类似于近期提出的One-Stage检测器，所提TOOD采用了类似的架构:Backbone-FPN-Head。考虑到效率与简单性，类似ATSS， TOOD在每个位置放置一个Anchor。

正如所讨论的，由于分类与定位任务的发散性，现有One-Stage检测器存在任务不对齐(Task Mis-Alignment)约束问题。本文提出通过显式方式采用T-head+TAL对2个任务进行对齐，见上图。T-head与TAL通过协同工作方式改善2个任务的对齐问题；

TOOD选取样本的方法具体来说:

  • 首先，T-head在FPN特征基础上进行分类与定位预测;

  • 然后，TAL基于所提任务对齐测度计算任务对齐信息;

  • 最后，T-head根据从TAL传回的信息自动调整分类概率与定位预测。

1、Task-Aligned Head

本文所提T-Head见下图，它具有非常简单的结构：特征提取器+TAP。

为增强分类与定位之间的相互作用，作者通过特征提取器学习任务交互（Task-Interactive）特征，如中蓝色框部分。这种设计不仅有助于任务交互，同时可以为2个任务提供多级多尺度特征。

假设X表示FPN特征，特征提取器采用N个连续卷积计算任务交互特征：

因此，通过特征提取器可以得到丰富的多尺度特征并用于送入到后续2个TAP模块中进行分类与定位对齐。

2、Task-Aligned Sample Assignment

为与NMS搭配，训练样例的Anchor分配需要满足以下规则：

  • 正常对齐的Anchor应当可以预测高分类得分，同时具有精确定位；

  • 不对齐的Anchor应当具有低分类得分，并在NMS阶段被抑制。

基于上述两个规则，作者设计了一种新的Anchor对齐度量以显式度量Anchor层面的对齐度。该对齐度量将集成到样本分配与损失函数中以动态提炼每个Anchor的预测。

Anchor Alignment metric 考虑到分类得分与IoU表征了预测质量，我们采用2者的高阶组合度量任务对齐度，公式定义如下：

其中，s与u分别表示分类得分与IoU值，而用于控制两者的影响。因此，t在联合优化中起着非常重要的作用，它激励网络动态的聚焦于高质量的Anchor上。

Training sample assignment 正如已有研究表明，训练样例分配对于检测器的训练非常重要。

为提升两个任务的对齐性，TOOD聚焦于任务对齐Anchor，采用一种简单的分配规则选择训练样本：对每个实例，选择m个具有最大t值的Anchor作为正样例，选择其余的Anchor作为负样例。然后，通过新的损失函数(针对分类与定位的对齐而设计的损失函数)任务进行训练。

  1. class TaskAlignedAssigner(nn.Layer):
  2.     def __init__(self, topk=13, alpha=1.0, beta=6.0, eps=1e-9):
  3.         super(TaskAlignedAssigner, self).__init__()
  4.         self.topk = topk
  5.         self.alpha = alpha
  6.         self.beta = beta
  7.         self.eps = eps
  9.     @paddle.no_grad()
  10.     def forward(self, pred_scores, pred_bboxes, anchor_points, num_anchors_list, gt_labels, gt_bboxes, pad_gt_mask, bg_index, gt_scores=None):
  11.         """
  12.         Task-Aligned Assigner计算步骤如下：
  13.         1. 计算所有 bbox与 gt 之间的对齐度
  14.         2. 选择 top-k bbox 作为每个 gt 的候选项
  15.         3. 将正样品的中心限制在 gt 内(因为Anchor-Free检测器只能预测大于0的距离)
  16.         4. 如果一个Anchor被分配给多个gt，将选择IoU最高的那个。
  17.         Args:
  18.             pred_scores (Tensor, float32): predicted class probability, shape(B, L, C)
  19.             pred_bboxes (Tensor, float32): predicted bounding boxes, shape(B, L, 4)
  20.             anchor_points (Tensor, float32): pre-defined anchors, shape(L, 2), "cxcy" format
  21.             num_anchors_list (List): num of anchors in each level, shape(L)
  22.             gt_labels (Tensor, int64|int32): Label of gt_bboxes, shape(B, n, 1)
  23.             gt_bboxes (Tensor, float32): Ground truth bboxes, shape(B, n, 4)
  24.             pad_gt_mask (Tensor, float32): 1 means bbox, 0 means no bbox, shape(B, n, 1)
  25.             bg_index (int): background index
  26.             gt_scores (Tensor|None, float32) Score of gt_bboxes, shape(B, n, 1)
  27.         Returns:
  28.             assigned_labels (Tensor): (B, L)
  29.             assigned_bboxes (Tensor): (B, L, 4)
  30.             assigned_scores (Tensor): (B, L, C)
  31.         """
  32.         assert pred_scores.ndim == pred_bboxes.ndim
  33.         assert gt_labels.ndim == gt_bboxes.ndim and gt_bboxes.ndim == 3
  35.         batch_size, num_anchors, num_classes = pred_scores.shape
  36.         _, num_max_boxes, _ = gt_bboxes.shape
  38.         # 计算GT与预测box之间的iou, [B, n, L]
  39.         ious = iou_similarity(gt_bboxes, pred_bboxes)
  40.         # 获取预测bboxes class score
  41.         pred_scores = pred_scores.transpose([0, 2, 1])
  42.         batch_ind = paddle.arange(end=batch_size, dtype=gt_labels.dtype).unsqueeze(-1)
  43.         gt_labels_ind = paddle.stack([batch_ind.tile([1, num_max_boxes]), gt_labels.squeeze(-1)], axis=-1)
  44.         bbox_cls_scores = paddle.gather_nd(pred_scores, gt_labels_ind)
  45.         # 计算bbox与 gt 之间的对齐度, [B, n, L]
  46.         alignment_metrics = bbox_cls_scores.pow(self.alpha) * ious.pow(self.beta)
  48.         # check the positive sample's center in gt, [B, n, L]
  49.         is_in_gts = check_points_inside_bboxes(anchor_points, gt_bboxes)
  51.         # 选择 top-k 预测 bbox 作为每个 gt 的候选项
  52.         is_in_topk = gather_topk_anchors(alignment_metrics * is_in_gts, self.topk, topk_mask=pad_gt_mask.tile([1, 1, self.topk]).astype(paddle.bool))
  54.         # select positive sample, [B, n, L]
  55.         mask_positive = is_in_topk * is_in_gts * pad_gt_mask
  57.         # 如果一个Anchor被分配给多个gt，将选择IoU最高的那个, [B, n, L]
  58.         mask_positive_sum = mask_positive.sum(axis=-2)
  59.         if mask_positive_sum.max() > 1:
  60.             mask_multiple_gts = (mask_positive_sum.unsqueeze(1) > 1).tile([1, num_max_boxes, 1])
  61.             is_max_iou = compute_max_iou_anchor(ious)
  62.             mask_positive = paddle.where(mask_multiple_gts, is_max_iou, mask_positive)
  63.             mask_positive_sum = mask_positive.sum(axis=-2)
  64.         assigned_gt_index = mask_positive.argmax(axis=-2)
  65.  
  66.         # assigned target
  67.         assigned_gt_index = assigned_gt_index + batch_ind * num_max_boxes
  68.         assigned_labels = paddle.gather(gt_labels.flatten(), assigned_gt_index.flatten(), axis=0)
  69.         assigned_labels = assigned_labels.reshape([batch_size, num_anchors])
  70.         assigned_labels = paddle.where(mask_positive_sum > 0, assigned_labels, paddle.full_like(assigned_labels, bg_index))
  71.  
  72.         assigned_bboxes = paddle.gather(gt_bboxes.reshape([-1, 4]), assigned_gt_index.flatten(), axis=0)
  73.         assigned_bboxes = assigned_bboxes.reshape([batch_size, num_anchors, 4])
  74.  
  75.         assigned_scores = F.one_hot(assigned_labels, num_classes)
  76.         # rescale alignment metrics
  77.         alignment_metrics *= mask_positive
  78.         max_metrics_per_instance = alignment_metrics.max(axis=-1, keepdim=True)
  79.         max_ious_per_instance = (ious * mask_positive).max(axis=-1, keepdim=True)
  80.         alignment_metrics = alignment_metrics / (max_metrics_per_instance + self.eps) * max_ious_per_instance
  81.         alignment_metrics = alignment_metrics.max(-2).unsqueeze(-1)
  82.         assigned_scores = assigned_scores * alignment_metrics
  83.  
  84.         return assigned_labels, assigned_bboxes, assigned_scores

3损失函数

3.1、分类损失varifocal loss

Focal loss定义：

其中a是前景背景的损失权重，p的y次是不同样本的权重，难分样本的损失权重会增大。当训练一个密集的物体检测器使连续的IACS回归时，本文从focal loss中借鉴了样本加权思想来解决类不平衡问题。但是，与focal loss同等对待正负样本的损失不同，而varifocal loss选择不对称地对待它们。varifocal loss定义如下：

其中p是预测的IACS得分，q是目标IoU分数。对于训练中的正样本，将q设置为生成的bbox和gt box之间的IoU（gt IoU），而对于训练中的负样本，所有类别的训练目标q均为0。

备注：Varifocal Loss会预测Iou-aware Cls_score（IACS）与分类两个得分，通过p的y次来有效降低负样本损失的权重，正样本选择不降低权重。此外，通过q（Iou感知得分）来对Iou高的正样本损失加大权重，相当于将训练重点放在高质量的样本上面。

  1. @staticmethod
  2. def _varifocal_loss(pred_score, gt_score, label, alpha=0.75, gamma=2.0):
  3.     weight = alpha * pred_score.pow(gamma) * (1 - label) + gt_score * label
  4.     loss = F.binary_cross_entropy(pred_score, gt_score, weight=weight, reduction='sum')
  5.     return loss

3.2、回归损失

1、GIoULoss

GIOU的计算很简单，对于两个bounding box A，B。我们可以算出其最小凸集（包围A、B的最小包围框）C。有了最小凸集，就可以计算GIOU：

计算方法很简单，从公式可以看出，GIOU有几个特点：

  • GIOU是IOU的下界，且取值范围为(-1, 1]。当两个框不重合时，IOU始终为0，不论A、B相隔多远，但是对于GIOU来说，A，B不重合度越高（离的越远），GIOU越趋近于-1。关于这点，下面再详细解释一下。

  • GIOU其实就是在IOU的基础上减掉了一个东西，这个减掉的东西，让我门避免了两个bbox不重合时Loss为0的情况。至于减掉的东西怎么去直观的理解，似乎不好找到一个很好的解释？

  • 可导：这一点需要强调下，由于max，min，分段函数(比如ReLU)这些都是可导的，所以用1-GIOU作为Loss是可导的。

当IOU=0时：

显然， A∪B值不变，最大化GIOU就是要最小化C，最小化C就会促成2个框不断靠近，而不是像最小化IOU那样loss为0。

YOLO V3涨了2个点，Faster RCNN，MaskRCNN这种涨点少了些。主要原因在于Faster RCNN，MaskRCNN本身的Anchor很多，出现完全无重合的情况比较少，这样GIOU和IOU Loss就无明显差别。所以提升不是太明显。

在TOOD中，bbox(Bouding box)通过对齐的anchor(具有更大的分类得分、更精确的定位)预测得到，这样的bbox通常经过NMS后仍可以得以保留。此外，t可以在训练阶段通过对损失加权选择高质量的bbox。因此，采用t度量bbox的质量，同时结合GIoU Loss定义了TOOD的Reg Loss如下:

  1. @register
  2. @serializable
  3. class GIoULoss(object):
  4.     """
  5.     Generalized Intersection over Union, see https://arxiv.org/abs/1902.09630
  6.     Args:
  7.         loss_weight (float): giou loss weight, default as 1
  8.         eps (float): epsilon to avoid divide by zero, default as 1e-10
  9.         reduction (string): Options are "none", "mean" and "sum". default as none
  10.     """
  11.  
  12.     def __init__(self, loss_weight=1., eps=1e-10, reduction='none'):
  13.         self.loss_weight = loss_weight
  14.         self.eps = eps
  15.         assert reduction in ('none', 'mean', 'sum')
  16.         self.reduction = reduction
  17.  
  18.     def bbox_overlap(self, box1, box2, eps=1e-10):
  19.         """calculate the iou of box1 and box2
  20.         Args:
  21.             box1 (Tensor): box1 with the shape (..., 4)
  22.             box2 (Tensor): box1 with the shape (..., 4)
  23.             eps (float): epsilon to avoid divide by zero
  24.         Return:
  25.             iou (Tensor): iou of box1 and box2
  26.             overlap (Tensor): overlap of box1 and box2
  27.             union (Tensor): union of box1 and box2
  28.         """
  29.         x1, y1, x2, y2 = box1
  30.         x1g, y1g, x2g, y2g = box2
  31.  
  32.         xkis1 = paddle.maximum(x1, x1g)
  33.         ykis1 = paddle.maximum(y1, y1g)
  34.         xkis2 = paddle.minimum(x2, x2g)
  35.         ykis2 = paddle.minimum(y2, y2g)
  36.         w_inter = (xkis2 - xkis1).clip(0)
  37.         h_inter = (ykis2 - ykis1).clip(0)
  38.         overlap = w_inter * h_inter
  39.  
  40.         area1 = (x2 - x1) * (y2 - y1)
  41.         area2 = (x2g - x1g) * (y2g - y1g)
  42.         union = area1 + area2 - overlap + eps
  43.         iou = overlap / union
  44.  
  45.         return iou, overlap, union
  46.  
  47.     def __call__(self, pbox, gbox, iou_weight=1., loc_reweight=None):
  48.         x1, y1, x2, y2 = paddle.split(pbox, num_or_sections=4, axis=-1)
  49.         x1g, y1g, x2g, y2g = paddle.split(gbox, num_or_sections=4, axis=-1)
  50.         box1 = [x1, y1, x2, y2]
  51.         box2 = [x1g, y1g, x2g, y2g]
  52.         iou, overlap, union = self.bbox_overlap(box1, box2, self.eps)
  53.         xc1 = paddle.minimum(x1, x1g)
  54.         yc1 = paddle.minimum(y1, y1g)
  55.         xc2 = paddle.maximum(x2, x2g)
  56.         yc2 = paddle.maximum(y2, y2g)
  57.  
  58.         area_c = (xc2 - xc1) * (yc2 - yc1) + self.eps
  59.         miou = iou - ((area_c - union) / area_c)
  60.         if loc_reweight is not None:
  61.             loc_reweight = paddle.reshape(loc_reweight, shape=(-1, 1))
  62.             loc_thresh = 0.9
  63.             giou = 1 - (1 - loc_thresh) * miou - loc_thresh * miou * loc_reweight
  64.         else:
  65.             giou = 1 - miou
  66.         if self.reduction == 'none':
  67.             loss = giou
  68.         elif self.reduction == 'sum':
  69.             loss = paddle.sum(giou * iou_weight)
  70.         else:
  71.             loss = paddle.mean(giou * iou_weight)
  72.         return loss * self.loss_weight

2、L1 loss

均绝对误差（Mean Absolute Error,MAE) 是指模型预测值f(x)和真实值y之间距离的均值，其公式如下：

忽略下标i ，设n=1，以f(x)−y为横轴，MAE的值为纵轴，得到函数的图形如下：

MAE曲线连续，但是在y−f(x)=0处不可导。而且 MAE 大部分情况下梯度都是相等的，这意味着即使对于小的损失值，其梯度也是大的。这不利于函数的收敛和模型的学习。但是，无论对于什么样的输入值，都有着稳定的梯度，不会导致梯度爆炸问题，具有较为稳健性的解。

相比于MSE，MAE有个优点就是，对于离群点不那么敏感。因为MAE计算的是误差y−f(x)的绝对值，对于任意大小的差值，其惩罚都是固定的。

loss_l1 = F.l1_loss(pred_bboxes_pos, assigned_bboxes_pos)

3、DF Loss

对于任意分布来建模框的表示，它可以用积分形式嵌入到任意已有的和框回归相关的损失函数上，例如最近比较流行的GIoU Loss。这个实际上也就够了，不过涨点不是很明显，我们又仔细分析了一下，发现如果分布过于任意，网络学习的效率可能会不高，原因是一个积分目标可能对应了无穷多种分布模式。如下图所示：

考虑到真实的分布通常不会距离标注的位置太远，所以我们又额外加了个loss，希望网络能够快速地聚焦到标注位置附近的数值，使得他们概率尽可能大。基于此，我们取了个名字叫Distribution Focal Loss (DFL)：

其形式上与QFL的右半部分很类似，含义是以类似交叉熵的形式去优化与标签y最接近的一左一右两个位置的概率，从而让网络快速地聚焦到目标位置的邻近区域的分布中去。

QFL和DFL的作用是正交的，他们的增益互不影响

  1. def _df_loss(self, pred_dist, target):
  2.     target_left = paddle.cast(target, 'int64')
  3.     target_right = target_left + 1
  4.     weight_left = target_right.astype('float32') - target
  5.     weight_right = 1 - weight_left
  6.     loss_left = F.cross_entropy(pred_dist, target_left, reduction='none') * weight_left
  7.     loss_right = F.cross_entropy(pred_dist, target_right, reduction='none') * weight_right
  8.     return (loss_left + loss_right).mean(-1, keepdim=True)

3.3、总损失

其中，a表示分类损失的权重系数，b表示回归损失的权重系数，c表示DFL损失的权重系数。

  1.     def get_loss(self, head_outs, gt_meta):
  2.         pred_scores, pred_distri, anchors, anchor_points, num_anchors_list, stride_tensor = head_outs
  3.  
  4.         anchor_points_s = anchor_points / stride_tensor
  5.         pred_bboxes = self._bbox_decode(anchor_points_s, pred_distri)
  6.  
  7.         gt_labels = gt_meta['gt_class']
  8.         gt_bboxes = gt_meta['gt_bbox']
  9.         pad_gt_mask = gt_meta['pad_gt_mask']
  10.         # Epoch小于100使用ATSS匹配
  11.         if gt_meta['epoch_id'] < self.static_assigner_epoch:
  12.             assigned_labels, assigned_bboxes, assigned_scores = \
  13.                 self.static_assigner(
  14.                     anchors,
  15.                     num_anchors_list,
  16.                     gt_labels,
  17.                     gt_bboxes,
  18.                     pad_gt_mask,
  19.                     bg_index=self.num_classes,
  20.                     pred_bboxes=pred_bboxes.detach() * stride_tensor)
  21.             alpha_l = 0.25
  22.         else:
  23.              # Epoch大于100使用TAL匹配
  24.             assigned_labels, assigned_bboxes, assigned_scores = \
  25.                 self.assigner(
  26.                 pred_scores.detach(),
  27.                 pred_bboxes.detach() * stride_tensor,
  28.                 anchor_points,
  29.                 num_anchors_list,
  30.                 gt_labels,
  31.                 gt_bboxes,
  32.                 pad_gt_mask,
  33.                 bg_index=self.num_classes)
  34.             alpha_l = -1
  35.         # rescale bbox
  36.         assigned_bboxes /= stride_tensor
  37.         # cls loss
  38.         if self.use_varifocal_loss:
  39.             one_hot_label = F.one_hot(assigned_labels, self.num_classes)
  40.             loss_cls = self._varifocal_loss(pred_scores, assigned_scores, one_hot_label)
  41.         else:
  42.             loss_cls = self._focal_loss(pred_scores, assigned_scores, alpha_l)
  43.  
  44.         assigned_scores_sum = assigned_scores.sum()
  45.         if paddle_distributed_is_initialized():
  46.             paddle.distributed.all_reduce(assigned_scores_sum)
  47.             assigned_scores_sum = paddle.clip(assigned_scores_sum / paddle.distributed.get_world_size(), min=1)
  48.         loss_cls /= assigned_scores_sum
  49.  
  50.         loss_l1, loss_iou, loss_dfl = self._bbox_loss(pred_distri, pred_bboxes, anchor_points_s, assigned_labels, assigned_bboxes, assigned_scores, assigned_scores_sum)
  51.         loss = self.loss_weight['class'] * loss_cls + self.loss_weight['iou'] * loss_iou + self.loss_weight['dfl'] * loss_dfl
  52.         out_dict = {
  53.             'loss': loss,
  54.             'loss_cls': loss_cls,
  55.             'loss_iou': loss_iou,
  56.             'loss_dfl': loss_dfl,
  57.             'loss_l1': loss_l1,
  58.         }
  59.         return out_dict

4模型推理与部署

4.1、模型推理

  1. # inference single image
  2. CUDA_VISIBLE_DEVICES=0 python tools/infer.py -c configs/ppyoloe/ppyoloe_crn_l_300e_coco.yml -o weights=https://paddledet.bj.bcebos.com/models/ppyoloe_crn_l_300e_coco.pdparams --infer_img=demo/000000014439_640x640.jpg
  3.  
  4. # inference all images in the directory
  5. CUDA_VISIBLE_DEVICES=0 python tools/infer.py -c configs/ppyoloe/ppyoloe_crn_l_300e_coco.yml -o weights=https://paddledet.bj.bcebos.com/models/ppyoloe_crn_l_300e_coco.pdparams --infer_dir=demo

4.2、导出ONNX

  1. # export inference model
  2. python tools/export_model.py configs/ppyoloe/ppyoloe_crn_l_300e_coco.yml --output_dir=output_inference -o weights=https://paddledet.bj.bcebos.com/models/ppyoloe_crn_l_300e_coco.pdparams
  3.  
  4. # install paddle2onnx
  5. pip install paddle2onnx
  6.  
  7. # convert to onnx
  8. paddle2onnx --model_dir output_inference/ppyoloe_crn_l_300e_coco --model_filename model.pdmodel --params_filename model.pdiparams --opset_version 11 --save_file ppyoloe_crn_l_300e_coco.onnx
  9.

4。3、导出TensorRT Engine

python tools/export_model.py configs/ppyoloe/ppyoloe_crn_l_300e_coco.yml -o weights=https://paddledet.bj.bcebos.com/models/ppyoloe_crn_l_300e_coco.pdparams -o trt=True

还需要测试 

python tools/export_model.py configs/ppyoloe/ppyoloe_crn_s_300e_coco.yml -o weights=ppyoloe_crn_l_300e_coco.pdparams -o trt=True

5参考

[1].https://github.com/PaddlePaddle/PaddleDetection

声明：本文内容由网友自发贡献，不代表【wpsshop博客】立场，版权归原作者所有，本站不承担相应法律责任。如您发现有侵权的内容，请联系我们。转载请注明出处：https://www.wpsshop.cn/article/detail/52146
推荐阅读
相关标签

Copyright © 2003-2013 www.wpsshop.cn 版权所有，并保留所有权利。

  

闽ICP备14008679号