一种基于雷达影像与电子海图信息的岸线配准方法
【专利摘要】本发明公开了一种基于雷达影像与电子海图信息的岸线配准方法，包括以下步骤：对雷达影像进行处理，获取其中的河/海岸线，并转换为全局坐标系下的像素位置坐标；根据S57海图格式标准，从电子海图中提取代表河/海岸线的线条信息，并转换为大地坐标系下的经纬度坐标；通过灰色关联处理手段，针对雷达影像和电子海图所提取出来的岸线坐标进行关联匹配，比较不同位置及尺度下的数据关联度，在这两种不同来源的数据之间实现粗配准；通过ICP配准处理手段，针对雷达影像和电子海图所提取出来的岸线坐标进行精配准。本发明提供的岸线识别方法识别精确度高，为其他目标识别划定范围做铺垫，有效降低雷达目标识别的错误率。
【专利说明】一种基于雷达影像与电子海图信息的岸线配准方法
【技术领域】
[0001] 本发明涉及地理信息【技术领域】，尤其涉及一种基于雷达影像与电子海图信息的岸 线配准方法。
【背景技术】
[0002] 在船舶现代化监管中，海事雷达是一种非常重要的监管手段。由于其工作环境复 杂，岸线周围错综复杂的建筑物、行驶车辆、风浪、水面等往往会对海事雷达产生较大的干 扰，出现杂波，有必要针对雷达图像进行处理，生成岸线的基本轮廓，进而为岸线识别，以及 水上目标识别做好铺垫，以便降低雷达目标识别的错误率。此外，对雷达影像和电子海图的 叠加显示，能够将电子海图的地理信息优势和雷达装备的实时探测能力结合在一起，从而 更加准确地获悉周围水域的障碍物信息，大幅提高船舶的避碰能力。因此，通过图像处理技 术，从雷达影像中提取出岸线数据，并与电子海图进行配准，实现雷达影像和电子海图的叠 加显示，在保障水上交通安全方面具有重要意义。
【发明内容】
[0003] 本发明要解决的技术问题在于针对现有技术中的缺陷，提供一种基于雷达影像与 电子海图信息的岸线配准方法。
[0004] 本发明解决其技术问题所采用的技术方案是：一种基于雷达影像与电子海图信息 的岸线配准方法，包括以下步骤：
[0005] 1)对雷达影像进行处理，获取其中的河/海岸线，并转换为全局坐标系下的像素 位置坐标；
[0006] 2)根据S57海图格式标准，从电子海图中提取代表河/海岸线的线条信息，并转换 为大地坐标系下的经纬度坐标；
[0007] 3)通过灰色关联处理手段，针对雷达影像和电子海图所提取出来的岸线坐标进行 关联匹配，比较不同位置及尺度下的数据关联度，在这两种不同来源的数据之间实现粗配 准；
[0008] 4)通过ICP配准处理手段，针对雷达影像和电子海图所提取出来的岸线坐标进行 精配准。
[0009] 按上述方案，所述步骤1)中雷达影像的图像处理使用包括连通域判别、边缘提 取、Gabor小波处理图像处理技术。
[0010] 按上述方案，所述步骤（1)中基于雷达影像的河/海岸线信息提取方法，包括以下 步骤：
[0011] (I. 1)对雷达影像进行灰度变换，将彩色图像转化为8位灰度图像；
[0012] (1. 2)对雷达影像进行连续帧间差分处理，提取雷达检测中的背景影像；定义第t 帧雷达检测影像为：U，其中i、j代表第i行，第j列的像素点；则连续三帧雷达影像可 以分别表示为
【权利要求】
1. 一种基于雷达影像与电子海图信息的岸线配准方法，其特征在于，包括以下步骤： 1) 对雷达影像进行处理，获取其中的河/海岸线，并转换为全局坐标系下的像素位置 坐标； 2) 根据S57海图格式标准，从电子海图中提取代表河/海岸线的线条信息，并转换为大 地坐标系下的经纬度坐标； 3) 通过灰色关联处理手段，针对雷达影像和电子海图所提取出来的岸线坐标进行关联 匹配，比较不同位置及尺度下的数据关联度，在这两种不同来源的数据之间实现粗配准； 4) 通过ICP配准处理手段，针对雷达影像和电子海图所提取出来的岸线坐标进行精配 准。
2. 根据权利要求1所述的岸线配准方法，其特征在于，所述步骤1)中雷达影像的图像 处理使用包括连通域判别、边缘提取、Gabor小波处理图像处理技术。
3. 根据权利要求1所述的岸线配准方法，其特征在于，所述步骤（1)中基于雷达影像的 河/海岸线信息提取方法，包括以下步骤： (I.1)对雷达影像进行灰度变换，将彩色图像转化为8位灰度图像； (1.2)对雷达影像进行连续帧间差分处理，提取雷达检测中的背景影像；定义第t帧雷 达检测影像为：其中i、j代表第i行，第j列的像素点；则连续三帧雷达影像可以分 别表示为：WjM7H;同时定义第t巾贞的背景图像为：(?+ (1. 3)将这三帧连续视频图像分为两组，第t-Ι帧和第t帧为一组，第t帧及第t+1 帧为另一组，分别对这两组视频的每一个象素点进行帧间差减，并将其绝对值保存在
(1.4)将{/?/}中各像素的值分别与一个预先给定的阈值Ttl进行比较，如果 某一点i，j上对应的作「/}，{1?1}的值均大于Ttl，则可以判断这一点在当前连续三帧图 像中都是运动的；反之，该点即为背景图像中的一部分，如公式（2)所示，对前30秒 的雷达检测影像反复进行背景提取处理，即可构建出完整的初始背景影像；
(1. 5)统计背景影像的灰度分布直方图，设定阈值T1，针对目标信息的灰度特性进行二 值化；如公式⑶所示：
(1.6)在得到二值化的背景图象后，利用区域生长法将灰度值为255的像素点按邻 域连通性原则合成一个区域，对连通区域像素值大小的统计和排序，如果区域面积过小或 者长宽比不符合岸线的一般特征，则认为是干扰信号而予以删除，最终选取出包含了岸 线有效信息的连通区域；其中，连通区域的选取采用四领域法，即目标像素P及其周围的 (x+l，y), (x，y-l), (x-l，y), (x，y+l)呈现联通状态则认为是同一连通域； (1. 7)采用拉普拉斯算子对包含岸线有效信息的连通域进行边缘提取，获取可能的岸 线边缘；拉普拉斯算子的定义如公式（4)、（5)所示：
(1.8)对背景边缘图像采用Hough变换获取岸线的方向信息，Hough变换的极坐标定义 为：xcosΘ+ysinΘ=P，通过对所有边缘像素点的极坐标（P，Θ)进行统计，选取最大的 Θ作为岸线的方向角信息，并以其为依据确定Gabor变换的角度，代入Gabor小波中实现对 岸线点的提取；二维Gabor小波变换以及傅里叶变换的定义如公式（6)、（7)、（8)所示：
4.根据权利要求1所述的岸线配准方法，其特征在于，所述步骤（3)中针对来自雷达影 像和电子海图的岸线信息进行灰色关联分析，确定几何形态最为接近的岸线位置，实现两 种岸线数据间的粗配准，具体包括以下步骤： (3. 1)针对雷达影像和电子海图的特征，确定灰色关联分析中的参考数据列和比较数 据列：由于雷达影像为位图形式，而电子海图为矢量图形式，因此将雷达影像中提取的岸线 坐标设定为参考数据列，定义为：R=Ir(k) |k=1，2, ···，!!，将电子海图中的岸线信息设定 为比较数据列，分别定义为=Ei =IeiGO|k=1，2, ···，!!}。之后比较R和若干AEi之间的 几何形态相似程度，判断其关联程度； (3. 2)雷达影像中的岸线信息为位图坐标，电子海图中的岸线信息为经纬度信息，数据 量纲不同，需根据公式（9)进行数据的无量纲化处理如下：
(3.3)根据公式（10)和（11)，计算比较数据列R与参考数据列Ei中对应数据的关联 系数如下：
(3. 4)取所有关联系数的平均值为数据列之间的关联度，并以此表示比较数据列R与 参考数据列Ei之间的关联程度，如公式（12)如示：
(3. 5)将关联度ri按大小排序，取关联度最大的比较数据列R与参考数据列Ei为同一 条岸线的粗配准结果。
5.根据权利要求1所述的岸线配准方法，其特征在于，所述步骤（4)中针对来自雷达影 像和电子海图的岸线信息采用了改进的ICP点云配准方法，实现了两种岸线数据间的精配 准，具体包括以下步骤： (4. 1)将粗配准后的雷达影像岸线数据存储在点集P中，定义为：{P|PieR3，i=1，2,…，N}，Pi为岸线上某点的坐标，定义为：(xpi，ypi);同时，将对应的电子海图岸线数据 存储在点集Q中，定义为：（QlqieR3,i= 1，2,…，N}，qi为岸线上某点的坐标，定义为： (Xpyqi);作为同一段岸线的非同源数据，p、Q两个点集之间应符合刚体变换，假定该刚体 变换是由一个2*2的旋转矩阵R和一个2*2的平移矩阵T实现的，如公式（13)所示：
(4. 2)根据最优匹配原则设计数据配准的目标函数，如公式（14)所示，点集R、T应使 该目标函数的值最小：
(4. 3)计算点集P中的每一个点在Q点集中的对应最近点； (4. 4)计算上述对应点对平均距离最小的刚体变换，求得旋转矩阵R和平移矩阵T; (4. 5)对Q使用上一步的平移矩阵T和旋转矩阵Q，得到新的变换点集，并替换获得新 的点集Q; (4.6)如果新的点集Q与参考点集P之间满足公式（14)中目标函数的要求，即平均距 离小于给定的阀值，则停止迭代，完成精配准；否则，以新的点集Q与参考点集P继续迭代下 去，直到达到目标函数的要求。
【文档编号】G06T7/00GK104463877SQ201410764435
【公开日】2015年3月25日 申请日期:2014年12月12日 优先权日:2014年12月12日
【发明者】谢磊, 刘建坤, 马枫, 初秀民, 曹臻, 王树武 申请人:武汉理工大学