本发明涉及高分辨SAR图像的检测技术领域，具体为一种基于超像素和随机森林的舰船目标检测方法。本发明融合了超像素分割与随机森林无监督聚类，将图像分块后进行聚类，准确检测出稀有舰船目标的同时，又较好地保留了船只的轮廓形状信息。较传统的CFAR舰船目标检测方法，本发明避开了杂波建模、参数估计、滑窗等耗时费力的操作。
1.一种基于超像素和随机森林的SAR图像舰船目标检测方法，其特征在于，包括以下步骤：步骤1、定义原图有M个像素，生成的超像素目为N；对图像的像素进行采样，采样步长为T初始化聚类中心Ck＝[lk,ak,bk,xk,yk] ，其中lk,ak,bk为像素在CIE-Lab空间的颜色值，xk,yk为像素坐标；
步骤2、将聚类中心移动到3×3邻域内的最小梯度位置；
步骤3、对每一个像素i，设置标签l(i)＝1，距离度量d(i)＝∞；
步骤4、对于聚类中心Ck的2S×2S邻域里的每一个像素i，计算Ck和像素i之间的距离D；
判断D是否小于d(i)，若是，就令d(i)＝D,l(i)＝k；若否，则重复步骤4直至完成对所有像素的判断；计算D的公式为：其中，li，ai，bi为像素i在CIE-Lab空间的颜色值，xi，yi为像素i的坐标；dlab表示颜色距离；dxy表示空间距离；m用来调整dxy的权值，取值范围为1-20；
步骤6、对每一个聚类中心，重复步骤4和步骤5，直到E小于阈值；
步骤7、得到SAR图像的超像素分割集合S＝{s1,s2,...,sN}，定义某一个像素区域为sq，计算sq在水平方向和垂直方向的灰度共生矩阵G0o和G90o；
步骤8、分别计算步骤7所得到的G0o和G90o的角二阶矩和熵，其中，角二阶矩熵 最终得到sq的纹理特征fq＝(ASM0o,ENT0o,ASM90o,ENT90o)T，q＝1,...,N；G(p,j)为灰度共生矩阵G的第p行第j列的值，k为灰度共生矩阵行数；ASM0o和ENT0o分别为G0o的角二阶矩和熵，ASM90o，ENT90o分别为G90o的角二阶矩和熵；
步骤9、对超像素集合S＝{s1,s2,...,sN}中的每一个超像素重复步骤7和步骤8，最终得到整幅图像的特征F＝{f1,f2,...,fN}；
步骤10、创建随机森林：从1到4中产生两个随机数，定义为a1,a2，设置树的最大深度depth；
步骤11、计算a1和a2维度上的平均特征差δ(a1,a2),计算公式为：步骤12、把根结点S分成左右孩子结点，分裂准则为：
其中，Sl和Sr分别为左孩子结点和右孩子结点；F(a1,v)是特征矩阵F的第a1行第v列的值，F(a2,v)同理；
步骤13、将得到的左右孩子结点视为后续分裂的根节点，继续执行步骤10到步骤12，直到达到树的最大深度depth；
步骤14、选取所包含元素最少的叶子节点，将其中每个超像素里的每一个像素灰度值置为255，其余的超像素区域置为0，由此检测显示出船只目标。
2.如权利要求1所述基于超像素和随机森林的SAR图像舰船目标检测方法，其特征在于：步骤14产生的是一幅二值图像，船只目标为白色，背景为黑色。
一种基于超像素和随机森林的SAR图像舰船目标检测方法
技术领域
[0001] 本发明涉及高分辨SAR图像的检测技术领域，特别涉及一种船只形状信息低受损的检测方法：基于超像素和随机森林的SAR图像舰船目标检测。
背景技术
[0002] 合成孔径雷达(Synthetic Aperture Radar，SAR)能全天时、全天候对地球进行实时监测。不同的地物目标有不同的散射特性，SAR图像能提供这种具有区分性的特征进行场景解译。高分辨SAR图像包含有更多的结构和形状信息和弱散射体，传统的检测方法很难检测出来。
[0003] 船舶检测作为一种重要的海事应用，得到了广泛的研究。恒虚警率检测(CFAR)是最传统的SAR图像舰船检测技术。它基于假设检验，依据背景的分布，获得自适应阈值，从而把虚警率维持在预设的水平上。其性能取决于滑动窗口结构的设计、杂波统计模型的分布和模型参数估计。但是随着模型的复杂度增加，参数估计成为一个棘手的问题。某些模型的概率密度函数甚至没有封闭的表达式，因此传统的参数估计方法，如最大似然估计(ML)，将不再适用。避免杂波建模和参数估计等难题是目标检测的一个新思路。
[0004] 把图像分块进行处理成为一种趋势。如果能给按照图像的边缘纹理进行切割，那么目标的轮廓信息就能被很好地保留，对于检测来说意义重大。随机森林是一种经典的集成学习模型，常常用来监督分类。但是分类需要大量的训练数据，而且在大多情况下，很难获得足够的数据用于训练。在无监督聚类的情况下，则不需要考虑数据量的问题，也不会出现过拟合现象。
发明内容
[0005] 针对上述存在问题，避免检测方法受杂波建模和参数估计的影响，本发明提出了一种基于超像素和随机森林的SAR图像舰船目标检测方法，以实现准确检测舰船目标并保留其轮廓的目的。
[0006] 本发明由以下步骤实现，其检测流程算法见附图1。
[0007] 步骤1、设原图有N个像素，生成的超像素数目为k。对图像的像素进行采样，采样步长为 初始化聚类中心 其中 为像素在CIE-Lab空间的颜色值， 为像素坐标。
[0008] 步骤2、将聚类中心移动到3×3邻域内的最小梯度位置。
[0009] 步骤3、对每一个像素i，设置标签l(i)＝1，距离度量d(i)＝∞。
[0010] 步骤4、对于Ck的2S×2S邻域里的每一个像素i，计算Ck和i之间的距离D；判断D是否小于d(i)。如果小于d(i)，就令d(i)＝D,l(i)＝k。计算D的公式为：
[0011]
[0012]
[0013]
[0014] m用来调整dxy的权值，一般为1-20，在算法中设置为10。
[0015] 步骤5、更新聚类中心，计算残差E
[0016] 步骤6、对每一个聚类中心，重复步骤4和步骤5，直到E小于阈值。
[0017] 步骤7、经过前六个步骤，得到SAR图像的超像素分割集合S＝{s1,s2,...,sN}。假设某一个像素区域为si，计算si在水平方向(0°)和垂直方向(90°)的灰度共生矩阵 和[0018] 步骤8、分别计算步骤7所得到的 和 的角二阶矩和熵。其中，角二阶矩熵 最终得到si的纹理特征
[0019] 步骤9、对超像素集合S＝{s1,s2,...,sN}中的每一个超像素重复步骤7和步骤8，最终得到整幅图像的特征F＝{f1,f2,...,fN}。
[0020] 步骤10、从1到4中产生两个随机数，假设为a1,a2。设置树的最大深度depth。
[0021] 步骤11、计算a1和a2维度上的平均特征差
[0022] 步骤12、把根结点S分成左右孩子结点，分裂准则为：
[0023]
[0024] 步骤13、将得到的左右孩子结点视为后续分裂的根节点，继续执行步骤10到步骤12，直到达到树的最大深度depth。
[0025] 步骤14、选取所包含元素(超像素)最少的叶子节点，将其中每个超像素里的每一个像素灰度值置为255。其余的超像素区域置为0，由此检测显示出船只目标。
[0026] 步骤1到步骤6是利用SLIC算法对SAR图像进行超像素分割的过程；步骤7到步骤9计算超像素块的纹理特征；步骤10到步骤14为创建随机森林，对超像素进行聚类的过程。本发明所利用的SLIC算法能生成紧凑、近似均匀的超像素，在运算速度，物体轮廓保持、超像素形状方面具有较高的水平，比较符合期望的分割效果。而使用随机森林进行无监督聚类，不需要考虑数据量的问题，并且综合每棵树的聚类结果，增强了检测结果的准确性。
[0027] 综上所述，较传统的SAR舰船目标检测方法，本发明避开了杂波建模、参数估计、滑窗等耗时费力的操作，将图像分块后进行聚类，准确检测出稀有舰船目标的同时，又较好地保留了船只的轮廓形状信息。
附图说明
[0028] 图1为本发明舰船检测流程；
[0029] 图2为262×292像素大小的直布罗陀海峡SAR图像船只检测图；
[0030] 图3为图2中标记船只的检测效果图；
[0031] 图4为217×304像素大小的直布罗陀海峡SAR图像船只检测图；
[0032] 图5为图4中标记船只的检测效果图。
具体实施方式
[0033] 下面以两个截取自英国直布罗陀海峡、分辨率为3米的TerraSAR图像为例，对本发明做进一步说明。
[0034] 将图2.(a)分割为1008个超像素区域，得到的划分结果如图2.(b)。不规则的划分保留了船只的轮廓信息。图2(d)是本发明的检测结果，它显示9个船只目标都被正确地检测到，并且它们的形状保存完好。然而，当采用通用的滑窗分割方法和基于灰度聚类的方法时，如图2.(c)所示，船只信息损失严重。为了查看检测结果的细节，我们在图2中标注了3个ship目标，其检测结果如图3所示。图3的第二行，是使用滑窗切割图像并直接利用像素灰度值作为聚类的特征所得到的结果，ship1,2和3仅仅是一些零散的亮块，很容易被误解为小的船舶目标。第三行显示出本发明在弱散射情况下，将船只作为整体完整地检测出来，保留了船舶结构。
[0035] 图4.(a)中有更多的船舶目标，它们尺寸的差异较大。同样，使用SLIC算法将其分割成如图4.(b)所示的1800个超像素区域。本发明的结果显示正确的检测，而没有使用超像素分割和纹理特征聚类的方法，船只信息丢失严重，如图4.(c)。我们在图4中标注了两个船只目标，如图5所示。在图5的第二行中，利用滑窗切割以及灰度作为特征检测到的ship 1是完全错误的，ship 2的轮廓被扭曲，而第三行应用本发明的结果则保留了船舶结构信息。
[0036] 基于超像素和随机森林的SAR图像舰船目标检测方法的优越性体现在对船只目标形状和轮廓上的较好保留，不需要大量的训练数据就能将船只作为一个整体检测出来，且回避了杂波建模、参数估计等棘手的问题。适用于形状结构信息丰富的高分辨SAR图像。
