局部霍夫变换与形态学优化的sar图像电力线检测方法
【技术领域】
[0001]本发明属于SAR图像电力线检测技术领域，特别涉及局部霍夫变换与形态学优化的SAR图像电力线检测方法。
【背景技术】
[0002]电力线是低空飞行环境下的重要威胁要素之一，严重威胁着飞行器的飞行安全，所以如何利用微波雷达快速检测危险电力线目标，保障飞行器低空飞行安全，是十分重要的科学问题。
[0003]合成孔径雷达(Synthetic Aperture Radar，SAR)是一种全天候、全天时、作用距离远的高分辨率微波遥感成像雷达，是低空环境感知的重要手段之一，利用所获得的SAR图像进行电力线检测，对低空环境威胁要素检测具有重要的现实意义和广阔的应用前景。但电力线目标较小、反射强度弱，只有当雷达视线与电力线垂直或左右偏离两个小偏角出现布拉格谐振(Bragger)效应时才会有较强散射，而且电力线目标易于与道路、田埂等地物线段混杂，造成SAR图像电力线检测误判导致虚警率高。
【发明内容】
[0004]本发明的目的在于针对所需解决的问题，提出局部霍夫变换与形态学优化的SAR图像电力线检测方法，提高了 SAR图像电力线检测能力，可应用于微波雷达在低空环境下的静止障碍物检测与定位，提升了飞行器低空环境中威胁要素的检测能力。
[0005]为实现上述技术目的，本发明采用如下技术方案予以实现。
[0006]局部霍夫变换与形态学优化的SAR图像电力线检测方法包括以下步骤:
[0007]步骤1，获取SAR原始图像，对SAR原始图像采用Edison算子进行边缘检测，获得包含边缘信息的二值图像；所述包含边缘信息的二值图像中，白色像素部分表示边缘信息，黑色像素部分表示非边缘信息；
[0008]步骤2，在上述包含边缘信息的二值图像中划分出多个边缘带；
[0009]步骤3，在划分出多个边缘带的二值图像中，对每个边缘带采用基于极坐标变换空间的霍夫变换方法进行直线检测，得出霍夫变换检测结果；霍夫变换检测结果为包含至少一条由白色像素线段的二值图像，所述白色像素线段指由白色像素点组成的线段；
[0010]步骤4，在霍夫变换检测结果的每个边缘带内，选取长度最长的线段，如果其余任一线段与长度最长的线段的最近距离小于设定距离阈值，并且对应线段与长度最长的线段的夹角小于设定夹角阈值，则将对应线段舍去，否则，将对应线段保留；最终得出SAR图像电力线检测结果。
[0011]本发明的有益效果为:1)本发明提出基于边缘带划分的局部霍夫变换直线检测方法。该方法在有效减少全局霍夫变换伪峰的同时提高了处理实时性。2)为了抑制局部霍夫变换后图像中重叠与分叉线段对后续电力线检测的干扰，本发明中提出了联合距离和夹角阈值的形态学优化方法，滤除了重叠与分叉线段，提高了电力线检测的可靠性。
【附图说明】
[0012]图1为本发明的局部霍夫变换与形态学优化的SAR图像电力线检测方法的流程图；
[0013]图2为实测数据实验中待处理的SAR原始图像；
[0014]图3为实测数据实验按照本发明采用Edison算子进行边缘检测后得到的二值图像；
[0015]图4为实测数据实验按照本发明得出的滤波处理后的含边缘带的二值图像；
[0016]图5为实测数据实验中对图3进行全局霍夫变换直线提取后结果示意图；
[0017]图6为实测数据实验中对图3按照本发明进行步骤4的局部霍夫变换处理的二值图像；
[0018]图7为实测数据实验中按照本发明的步骤4舍弃每个边缘带中与长度最大线段距离相近且夹角较小的重复或分叉线段后的二值图像；
[0019]图8为实测数据实验中按照本发明针对图7滤除长度较短线段后的SAR图像电力线检测结果示意图。
【具体实施方式】
[0020]下面结合附图对本发明作进一步说明:
[0021]参照图1，为本发明的局部霍夫变换与形态学优化的SAR图像电力线检测方法的流程图。该局部霍夫变换与形态学优化的SAR图像电力线检测方法包括以下步骤:
[0022]步骤1，获取SAR原始图像，对SAR原始图像采用Edison算子进行边缘检测，获得包含边缘信息的二值图像；上述包含边缘信息的二值图像中，白色像素部分表示边缘信息，黑色像素部分表示非边缘信息，边缘信息主要指电力线、田埂、房屋、道路等的边缘信息。
[0023]步骤2，在上述包含边缘信息的二值图像中划分出多个边缘带，每个边缘带就是一条连续的线段、折线或曲线。
[0024]其具体子步骤为:
[0025](2.1)将上述包含边缘信息的二值图像的每个像素点赋予标志位flag，每个像素点的标志位flag的初始值为O。令k = 1，2，...，K，K表示上述包含边缘信息的二值图像中像素点的个数；当k = I时跳至子步骤(2.2)。
[0026](2.2)若上述包含边缘信息的二值图像的第k个像素点为白色像素点(根据灰度值判断，灰度值为255的像素点为白色像素点)且其标志位flag的值为0，则将上述包含边缘信息的二值图像的第k个像素点作为边缘带像素点，跳至子步骤(2.3);若上述包含边缘信息的二值图像的第k个像素点为黑色像素点(根据灰度值判断，灰度值为O的像素点为黑色像素点)或上述包含边缘信息的二值图像的第k个像素点的标志位flag的值为1，则跳至子步骤(2.5)。
[0027](2.3)记录每个边缘带像素点的坐标(即在二值图像中的行数和列数)并将其标志位flag的值更新为I，跳至子步骤(2.4)。
[0028](2.4)在每个边缘带像素点的8连通邻域中搜索标志位flag的值为O的白色像素点，如果搜索到，则将搜索到的每个像素点作为更新后的边缘带像素点，返回至子步骤(2.3)，如果未搜索到，则跳至子步骤(2.5) ο
[0029]下面对像素点的8连通邻域作出解释，像素点的8连通领域指:由上方相邻的像素点、下方相邻的像素点、左方相邻的像素点、右方相邻的像素点左上方相邻的像素点、左下方相邻的像素点、右上方相邻的像素点、右下方相邻的像素点组成的像素点的集合。显然当像素点不是边缘像素点时，像素点的8连通领域包括8个像素点；当像素点位于图像的一边但不是图像四角的像素点时，像素点的8连通领域包括5个像素点；当像素点是图像四角的像素点时，像素点的8连通领域包括3个像素点。
[0030](2.5)若k〈K，则令k的值自增1，返回至子步骤(2.2);若k = K，则说明在上述包含边缘信息的二值图像中划分出多个边缘带的过程结束。
[0031]步骤2的过程可以看做是在上述包含边缘信息的二值图像中对像素点进行逐个搜索的过程，搜索的顺序可以人为设定。在步骤2之后，上述包含边缘信息的二值图像中所有相邻的边缘像素点被划分各个不同的边缘带。
[0032]步骤3，像素点个数阈值滤波并利用局部霍夫变换方法进行直线检测。
[0033]在划分出多个边缘带的二值图像中，将像素点个数小于设定像素点数阈值的边缘带舍去(将对应的边缘带的每个像素点更新后黑色像素点)，得到滤波处理后的含边缘带的二值图像；在上述滤波处理后的含边缘带的二值图像中，对每个边缘带采用基于极坐标变换空间的霍夫变换方法进行直线检测，得出霍夫变换检测结果；霍夫变换检测结果为包含至少一条由白色像素点组成的线段的二值图像。
[0034]优选地，在步骤3中，设定像素点数阈值为30。
[0035]具体地，霍夫变换方法通过空间变换将在图像空间中的直线检测问题转变为参数空间里点的检测问题，在参数空间中进行的累加统计完成检测工作。根据参数域选取的不同，有直角坐标系中的参数空间变换和极坐标系中的参数空间变换两种，本发明采用基于极坐标变换空间的霍夫变换方法。其变换关系如下:设图像空间中一条直线L，直线上的每一点(x，y)都满足以下方程时，参数空间H为r-θ空间，
[0036]r = X cos Θ+y sin Θ (l)r为原点到直线L的距离，Θ为直线L过原点的垂线与X轴正方向的夹角。
[0037]对每个边缘带采用基于极坐标变换空间的霍夫变换方法进行直线检测的过程包括以下子步骤:
[0038](3.1)令i = 1，2，...，M，M表示待处理的二值图像(上述滤波处理后的二值图像)中边缘带的个数；当i = I时，跳至子步骤(3.2);
[0039](3.2)对第i个边缘带进行基于极坐标变换空间的霍夫变换，将第i个边缘带按照公式(I)变换到参数空间中，获得相应参数空间矩阵H，其行对应着Θ，列对应着r，其数值对应着过该点直线数的的累加值。
[0040](3.3)对参数空间矩阵H进行峰值检测，保留矩阵H中数值最大的5个元素，记录其对应的行列号，若参数空间矩阵H中元素个数不足5个，则保留矩阵H的全部元素。
[0041](3.