本发明涉及海洋数据共享应用领域，具体而言，尤其涉及一种船载雷达图像海上目标空间数据的共享方法。
背景技术：
德国、英国、加拿大、俄罗斯和中国的公司不断尝试船载雷达图像与电子海图数据的共享应用。例如，1994年加拿大近海ecdis-3000系统实现了电子海图与船载雷达的图像叠加，保证了大雾天气下的安全航行。由于技术水平低、硬件设施差、电子海图信息量大、雷达信号处理量大等因素，当时的数据融合效率较低。后来，挪威的hamilton公司和英国的transas公司在这一领域不断取得突破，提高了电子海图与船载雷达图像数据融合的效率。目前，一些产品可以实现舰载雷达与电子海图的快速信息集成，如美国的dkartnavigator。然而，由于商业竞争，他们的技术还没有公开。此外，transas和sperry还开发了成熟的雷达信号处理卡，可直接安装在通用pc的pci插槽中，当雷达信号处理卡与雷达收发天线单元连接时，可直接处理雷达视频信号，记录和跟踪运动目标，控制收发器等。这样，雷达视频信号就可以通过二次开发转化为图像并保存下来。但是上述保存的船载雷达视频图像数据并不能共享应用。
技术实现要素：
根据上述提出的技术问题，而提供一种船载雷达图像海上目标空间数据的共享方法。本发明将船载雷达原始图像中的目标空间数据导出，共享应用到地理信息系统相关平台。
本发明采用的技术手段如下：
一种船载雷达图像海上目标空间数据的共享方法，包括如下步骤：
s1、提取船载雷达图像中的海上目标；
s2、提取船载雷达图像中的海上目标的轮廓；
s3、海上目标轮廓边界点空间数据输出。
进一步地，所述步骤s1具体为：
s11、坐标系转换，将原始船载雷达图像从极坐标系转换至笛卡尔坐标系；
s12、在转换坐标系后的图像中设定捕获区；
s13、采用灰度阈值法与连续像元面积阈值法，提取图像中的海上目标。
进一步地，所述步骤s2具体为：
s21、坐标系转换，将提取的海上目标图像从笛卡尔坐标系转换至极坐标系；
s22、根据船首向，对识别的目标进行图像旋转，提取目标轮廓的平面直角坐标数据。
进一步地，所述步骤s3具体为：
s31、通过投影变换，将本船的gps地理坐标数据变换成平面直角坐标数据；
s32、根据目标轮廓边界点的极坐标数据(r,θ)与本船的平面直角坐标数据(x本船，y本船)，计算得到目标轮廓边界点的平面直角坐标数据(x，y)；其计算公式如下：
s33、通过投影变换，逐步将目标轮廓每个边界点的平面直角坐标变换至以经纬度为单位的地理坐标，输出目标轮廓的地理坐标数据。
进一步地，所述步骤s13具体为：
s131、先采用灰度阈值法，对图像进行二值化，初步提取海上目标；
s132、再采用连续像元面积阈值法，提取图像中准确的海上目标。
较现有技术相比，本发明具有以下优点：
1、本发明提供的船载雷达图像中的海上目标，通过在船载雷达数据采集的原始笛卡尔坐标系下，实现了图像中高亮海上目标的快速提取。
2、本发明提供的船载雷达图像中的海上目标轮廓的经纬度信息，根据船舶gps信息、船首向信息，通过空间位置计算与投影变换，提取了船载雷达图像中的海上目标轮廓的经纬度信息，可以使其快速融入电子海图或者其他地理信息平台。
基于上述理由本发明可在海洋数据共享应用等领域广泛推广。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图做以简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明方法流程图。
图2为本发明极坐标系下的船载雷达图像。
图3为本发明极坐标系与笛卡尔坐标系的转换图。
图4为本发明笛卡尔坐标系的船载雷达图像。
图5为本发明捕获区设定图。
图6为本发明采用灰度阈值法后的图像二值化示意图。
图7为本发明提取的图像中的海上目标示意图。
图8为本发明极坐标系下识别的目标示意图。
图9为本发明识别目标进行船首向方位旋转图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
如图1所示，本发明提供了一种船载雷达图像海上目标空间数据的共享方法，包括如下步骤：
s1、提取船载雷达图像中的海上目标；
进一步地作为本发明优选的实施方式，所述步骤s1具体为：
s11、坐标系转换，将原始船载雷达图像从极坐标系转换至笛卡尔坐标系；如图2所示，为在极坐标系下船载雷达图像；如图3所示，为极坐标系与笛卡尔坐标系的转换图；如图4所示，为坐标系转换后，在笛卡尔坐标系下的船载雷达图像；
s12、如图5所示，在转换坐标系后的图像中设定捕获区；
s13、采用灰度阈值法与连续像元面积阈值法，提取图像中的海上目标。
进一步地作为本发明优选的实施方式，所述步骤s13具体为：
s131、先采用灰度阈值法，对图像进行二值化，如图6所示，本实施例在matlab图像处理软件中设置bw＝im2bw(img,120/255)，初步提取海上目标；
s132、再采用连续像元面积阈值法，提取图像中准确的海上目标。本实施例在matlab图像处理软件中设置bw＝ismember(l,find([s.area]>20))，如图7所示，为提取的图像中的海上目标。
s2、提取船载雷达图像中的海上目标的轮廓；
进一步地作为本发明优选的实施方式，所述步骤s2具体为：
s21、坐标系转换，将提取的海上目标图像从笛卡尔坐标系转换至极坐标系；如图8所示，为在极坐标系下识别的目标；
s22、如图9所示，根据船首向，对识别的目标进行图像旋转，提取目标轮廓的平面直角坐标数据。
s3、海上目标轮廓边界点空间数据输出。
进一步地作为本发明优选的实施方式，所述步骤s3具体为：
s31、通过投影变换，将本船的gps地理坐标数据变换成平面直角坐标数据；
s32、根据目标轮廓边界点的极坐标数据(r,θ)与本船的平面直角坐标数据(x本船，y本船)，计算得到目标轮廓边界点的平面直角坐标数据(x，y)；其计算公式如下：
s33、通过投影变换，逐步将目标轮廓每个边界点的平面直角坐标变换至以经纬度为单位的地理坐标，输出目标轮廓的地理坐标数据。如表1所示：
表1海上目标边界点的地理坐标采样
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的范围。
1.一种船载雷达图像海上目标空间数据的共享方法，其特征在于，包括如下步骤：
s1、提取船载雷达图像中的海上目标；
s2、提取船载雷达图像中的海上目标的轮廓；
s3、海上目标轮廓边界点空间数据输出。
2.根据权利要求1所述的船载雷达图像海上目标空间数据的共享方法，其特征在于，所述步骤s1具体为：
s11、坐标系转换，将原始船载雷达图像从极坐标系转换至笛卡尔坐标系；
s12、在转换坐标系后的图像中设定捕获区；
s13、采用灰度阈值法与连续像元面积阈值法，提取图像中的海上目标。
3.根据权利要求1所述的船载雷达图像海上目标空间数据的共享方法，其特征在于，所述步骤s2具体为：
s21、坐标系转换，将提取的海上目标图像从笛卡尔坐标系转换至极坐标系；
s22、根据船首向，对识别的目标进行图像旋转，提取目标轮廓的平面直角坐标数据。
4.根据权利要求1所述的船载雷达图像海上目标空间数据的共享方法，其特征在于，所述步骤s3具体为：
s31、通过投影变换，将本船的gps地理坐标数据变换成平面直角坐标数据；
s32、根据目标轮廓边界点的极坐标数据(r,θ)与本船的平面直角坐标数据(x本船，y本船)，计算得到目标轮廓边界点的平面直角坐标数据(x，y)；其计算公式如下：
s33、通过投影变换，逐步将目标轮廓每个边界点的平面直角坐标变换至以经纬度为单位的地理坐标，输出目标轮廓的地理坐标数据。
5.根据权利要求2所述的船载雷达图像海上目标空间数据的共享方法，其特征在于，所述步骤s13具体为：
s131、先采用灰度阈值法，对图像进行二值化，初步提取海上目标；
s132、再采用连续像元面积阈值法，提取图像中准确的海上目标。
技术总结