CN105513041A - 一种大幅面遥感图像海陆分割的方法及系统 - Google Patents

Abstract

本发明适用于遥感图像处理技术领域，提供了一种大幅面遥感图像海陆分割的方法，步骤包括：步骤A，对所述遥感图像进行海陆填充，获取海陆粗分割结果图；步骤B，结合所述海陆粗分割结果图对所述遥感图像进行分块提取，获取海陆精细分割结果图；步骤C，将所述海陆粗分割结果图与所述海陆精细分割结果图进行融合，获取海陆分割结果图。本发明实现了一种自动的大幅面遥感图像的海陆分割方法，较好地解决了现有的海陆分割技术中，在近岸、近海多岛屿区域、地貌特征复杂的陆海区域等复杂背景下的海陆分割精度较低、分割效果几乎失效的问题。

Description

一种大幅面遥感图像海陆分割的方法及系统

技术领域

本发明属于遥感图像处理技术领域，尤其涉及一种大幅面遥感图像海陆分割的方法及系统。

背景技术

近年来，基于遥感图像的海上目标检测、识别越来越引起研究人员的关注。然而，大量的研究结果表明，针对遥感图像的海洋目标检测算法在遥感图像中包含有陆地区域或海岛区域的情况下，直接检测的结果总是在较低的检测概率和较高的虚警概率内进行波动。为了让海洋目标检测算法仅仅作用于海洋遥感区域，进而提高算法检测率，降低虚警、漏检，必须在对海洋目标进行检测前对大幅面遥感图像进行陆地、海岛区域屏蔽或移除，此即为对遥感图像进行海陆分割处理的目的。

现有技术提出的针对光学遥感图像进行海陆分割的方法主要有以下几种：

一、利用人工手动标注的方式实现对光学遥感图像的海陆分离。该方法十分耗时，计算效率低，不能满足大数据量的应用需求。

二、利用基于阈值分割的海陆分割方法。基于阈值分割方法主要有OTSU阈值法、基于最大熵等方法。这类分割方法一般会根据图像的直方图分布选取阈值对图像进行分割，阈值选取的好坏直接决定了海陆分割的效果。对于图像的直方图呈现双峰分布或近似双峰分布的理想情况，阈值选取较为简单，但当图像灰度分布较为复杂，如近海区域灰度接近或低于海洋区域灰度时，不容易自动选取阈值，因此该类算法有一定局限性。

三、利用基于纹理的方法，特别是将纹理特征与其他特征融合在一起进行海陆分割的方法。例如，将灰度图像和纹理特征图融合后进行海陆分离；将灰度图像与LBP特征图融合后利用阈值来实现海陆分离等。这类海陆分离方法一般对局部遥感图像的海陆分割效果较佳，但当图像背景复杂时，或将这类方法应用至大幅面遥感图像时，这类方法的分割效果有待改进。

四、利用简单统计模型的海陆分割方法。例如，利用贝叶斯统计模型的海陆分割方法，以及利用阈值分割和高斯统计模型相结合的海陆分割方法。但经调研发现，上述算法在应用过程中都需人工手动辅助，操作较为繁琐；另外，这类算法也不适合直接应用于大幅面遥感图像中。

上述这几类传统的海陆分割算法一般应用在小尺度范围内对遥感图像进行海陆分割，而对于大幅面遥感图像特别是背景复杂的遥感图像，这些算法的分割结果较差，这在近岸、近海多岛屿区域以及有云雾背景下的海洋区域表现的尤为明显。

综上所述，现有技术提出的针对光学遥感图像进行海陆分离的方法存在局部有效、分割精度低的问题。

发明内容

本发明所要解决的技术问题在于提供一种大幅面遥感图像海陆分割的方法及系统，旨在解决现有技术提出的针对光学遥感图像进行海陆分离的方法存在局部有效、分割精度低的问题。

本发明是这样实现的，一种大幅面遥感图像海陆分割的方法，步骤包括：

步骤A，对所述遥感图像进行海陆填充，获取海陆粗分割结果图；

步骤B，结合所述海陆粗分割结果图对所述遥感图像进行分块提取，获取海陆精细分割结果图；

步骤C，将所述海陆粗分割结果图与所述海陆精细分割结果图进行融合，获取海陆分割结果图。

本发明与现有技术相比，有益效果在于：本发明实现了一种自动的大幅面遥感图像的海陆分割方法，较好地解决了现有的海陆分割技术中，在近岸、近海多岛屿区域、地貌特征复杂的陆海区域等复杂背景下的海陆分割精度较低、分割效果几乎失效的问题。本发明可以有效处理复杂背景情形下的大幅面遥感图像的海陆分割问题，分割结果更为精准，对基于海陆分离的相关技术研究如海洋目标检测影响较小，适用于对遥感图像进行批量处理，实用性更强。

附图说明

图1是本发明实施例提供的一种大幅面遥感图像海陆分割的方法的流程图，

图2是本发明实施例提供的大幅面遥感图像的海陆分割方法的具体流程图；

图3a是本发明实施例中为粗分割步骤中的遥感图像添加海岸线的流程图；

图3b是本发明实施例中需用到的一卫星遥感多光谱原始示例图；

图3c是本发明实施例中，在遥感示例图中添加海岸线后的结果图；

图4a是本发明中粗分割步骤中的遥感图像海陆填充流程图；

图4b是本发明实施例中，遥感示例图中海陆粗分割后的结果图；

图4c是本发明实施例中得到的局部海陆粗分割结果图；

图5是本发明实施例中海陆精细分割步骤中的自动提取训练集子图像和测试集子图像流程图；

图6是本发明实施例中获取位于有效遥感影像区域的海陆粗分割结果实例图；

图7是本发明实施例中从原始遥感图像中分别提取全陆地区域、全海洋区域、以及需进行精细海陆分割处理的海岸线附近区域和海岛区域流程图；

图8是本发明实施例中自动提取相关区域子块的结果示例图；

图9a是图10中提取的部分全陆地训练集子图像；

图9b是图10中提取的部分全海洋训练集子图像；

图10a是本发明实施例中全陆地区域训练集对应的特征集的概率密度曲线；

图10b是本发明实施例中全海洋区域训练集对应的特征集的概率密度区曲线；

图11a是本发明实施例中全陆地训练集特征概率分布评估结果；

图11b是本发明实施例中全海洋训练集特征概率分布评估结果；

图12是本发明实施例中对测试集海陆精细分割结果进行形态学处理步骤流程图；

图13是本发明实施例中对所述海陆粗分割结果图与所述海陆精细分割结果图进行融合的流程图；

图14为本发明实施例提供的一种大幅面遥感图像海陆分割的系统的结构示意图。

图15为本发明实施例提供的一种大幅面遥感图像海陆分割的系统的详细结构示意图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

为了解决现有技术存在的问题，本发明提出的大幅面遥感图像的海陆分割方法主要包括基于地理信息的遥感图像海陆粗分割和基于统计学习的海陆精细分割两个部分。该方法首先利用遥感图像的经纬度信息，在全球一致化多层次高分辨率海岸线数据库(GSHHS,AGlobalSelf-consistentHierarchicalHigh-resolutionShorelineDatabase)的基础上填充大部分陆地区域，实现大幅面遥感图像的海陆粗分割，从而避免了陆地区域背景过于复杂、地貌信息过于丰富、不易分析的问题。因海岸线数据库精度、人工填海造陆、卫星图像拍摄角度等因素会造成部分遥感图像陆地区域，特别是海岸线附近区域和海岛区域未被粗分割填充，因此需对这类区域进行精细海陆分割。根据海陆粗分割结果，提取需进一步精细分割处理的区域子图像作为测试集，提取全陆地区域和全海洋区域子图像作为训练集，对训练集图像建立基于局部信息熵特征的自适应概率统计学习模型，利用该模型对测试集子图像进行基于局部信息熵的海陆精细分割。将基于地理信息的海陆粗分割结果和海陆精细分割结果融合在一起，并对融合结果进行形态学处理，最终得到完整的大幅面遥感图像海陆分割图像。

如图1所示，为本发明实施例提供的一种大幅面遥感图像海陆分割的方法，步骤包括：

A，对所述遥感图像进行海陆填充，获取海陆粗分割结果图；

B，结合所述海陆粗分割结果图对所述遥感图像进行分块提取，获取海陆精细分割结果图；

C，将所述海陆粗分割结果图与所述海陆精细分割结果图进行融合，获取海陆分割结果图。

以下结合图2对本发明实施例提供的一种大幅面遥感图像海陆分割的方法进行详细的阐述：

A：进行海陆粗分割

遥感图像的海陆粗分割过程是基于地理信息实现的，即借助于遥感图像4个顶点对应的经纬度信息以及全球海岸线数据库GSHHS实现的。GSHHS数据库根据用户的不同需求提供不同分辨率的海岸线数据文件，包括全分辨率数据文件gshhs_f.b、高分辨率数据文件gshhs_h.b、中等分辨率数据文件gshhs_i.b、低分辨率数据文件gshhs_l.b、粗分辨率数据文件gshhs_c.b。本节的海陆粗分割方法使用的海岸线数据文件为全分辨率数据文件gshhs_f.b，即通过对gshhs_f.b文件中的原始经纬度数据进行组织、变换来实现为大幅面光学遥感图像添加海岸线、海陆填充等功能。本实施例的海陆粗分割过程包括S1和S2两个步骤：

S1：利用GSHHS数据库为遥感图像添加海岸线。

其中，步骤S1为遥感图像添加海岸线的流程，如图3a所示，具体步骤如下：

S11：输入某一遥感图像，读取与该遥感图像相对应的地理信息文件(.xml)，获得遥感图像4个顶点A、B、C、D对应的经纬度坐标信息(lat1,lon1),(lat2,lon2),(lat3,lon3),(lat4,lon4)。

S12：根据遥感图像4个顶点对应的经纬度坐标信息，可确定一个四边形ABCD的经纬度区域，从gshhs_f.b海岸线数据文件中读取位于该经纬度区域范围内的所有离散海岸线经纬度数据，即(latt1,lonn1),(latt2,lonn2),(latt3,lonn3),…,(lattk,lonnk),…。

S13：根据遥感图像分辨率的大小获取遥感图像4个顶点A、B、C、D对应的顶点像素坐标，设遥感图像分辨率大小为MⅹN，则遥感图像4个顶点对应的顶点像素坐标分别为(1,1),(1,N),(M,1),(M,N)；计算遥感图像的顶点像素坐标与经纬度坐标之间的线性关系，根据所述线性关系获取所述遥感图像中每个像素点对应的经纬度数据。

S14：将位于四边形经纬度区域范围内的离散海岸线经纬度数据(latt1,lonn1),(latt2,lonn2),(latt3,lonn3),…,(lattk,lonnk),…插值映射为对应的遥感图像像素坐标，得到离散海岸线像素坐标。

S15：在遥感图像中，将S14步骤中获取的离散海岸线像素坐标连成线，形成海岸线，生成海岸线添加图。

图3b为一CBERS2B卫星多光谱图像，图3c为利用上述方法在图中添加海岸线后的海岸线添加图。

S2：对所述海岸线添加图进行海陆填充，得到海陆粗分割结果图。

步骤S2对海岸线添加图进行海陆填充的流程如图4a所示，具体步骤包括：

S21：读取海岸线添加图四个顶点的经纬度坐标，根据经纬度坐标确定一个四边形经纬度区域。

S22：依次读取GSHHS数据库中的经纬度数据段，其中，在该数据中每一个经纬度数据段为一封闭陆地或海岛经纬度数据。

S23：判断该经纬度数据段与海岸线添加图顶点对应的四边形经纬度范围是否有重合区域；

S24：如若该经纬度数据与海岸线添加图的四边形经纬度范围有重合区域，保留重合区域范围的经纬度数据；

S25：判断步骤S24获取的经纬度数据是否为一封闭数据段，如果该经纬度数据为一非封闭数据段数据，需将该非封闭数据段数据重新构造成一个封闭区域，并用数据段分割标志(NaN,NaN)将该非封闭数据段数据与其他数据段数据隔开；

S26：若判断步骤S24获取的经纬度数据为一封闭数据段，则将该经纬度数据均变换为图像的像素坐标，并将位于这些像素坐标围成的曲线范围内的像素点标记成陆地区域或海岛区域；

S27：利用区域填充算法，将S26像素点范围内的区域填充成一个完整的陆地块或海岛块；

S28：判断GSHHS数据库中是否还有其他封闭数据段，若存在重复步骤S22至步骤S27直至海陆填充完毕，得到海陆粗分割结果图。

图4b为遥感图像图3c进行海陆填充后的结果图，即海陆粗分割结果图。

B，海陆精细分割

虽然基于地理信息的海陆填充算法可以有效填充遥感图像中的大部分陆地区域，但由图4b也可以看出，填充的陆地区域有部分是遥感图像无效区域，如图像左上角的区域，该区域并不是遥感图像有效区域，但仍被该算法判别为陆地部分。另外，由于人工开发导致海岸线变迁、数据库中海岸线数据精准程度、遥感卫星图像拍摄角度等因素的影响，上述基于地理信息的海陆粗分割过程的陆地填充效果在海岸线附近区域、海岛周边区域效果并不理想。图4c为从图4b中截取的四个海岸线附近区域、海岛周边区域，由该图可以看出，部分陆地区域并没有被海陆填充算法划分为陆地部分，部分识别出来的海岛区域与遥感图像中实际海岛区域出现位置偏离的现象。因此，在海岸线边缘区域以及海岛区域需进行进一步的精细海陆分割处理。

本发明中，对于海岸线边缘区域以及海岛区域的海陆精细分割方法是基于全陆地区域和全海洋区域统计学习模型实现的，包括S3-S7五个步骤：

S3：在上述海陆粗分割基础上，对遥感图像进行分块，自动提取需要进一步精细分割处理的海岸线边缘区域块以及海岛区域块，同时提取一些全陆地区域块和全海洋区域块；以全陆地区域块、全海洋区域块作为训练集，以海岸线边缘区域块、海岛区域块作为测试集。

自动提取训练集子图像和测试集子图像区域是基于原始遥感图像信息和海陆粗分割后的图像相关信息来完成的。如图5所示，具体提取步骤如下：

S31：确定遥感图像的有效影像区域，形成遥感影像的区域。图3b中的倾斜四边形影像区域即为有效遥感影像区域S_p1p2p3p4(p1,p2,p3,p4对应四边形的4个顶点)，该四边形区域外围的黑色区域为无效影像图像区域，可以通过阈值化处理和形态学操作得到遥感图像的有效影像区域；

S32：利用所述有效影像区域确定所述海陆粗分割结果图中的有效陆地区域。图4b中的纯白色区域S_redlandmark为海陆粗分割算法识别出来的陆地区域，而该陆地区域与步骤S31中的有效影像区域相重叠的部分即为海陆粗分割结果图中的有效陆地区域，如图6中的纯白色区域S_truelandmark所示；

S33：根据所述有效陆地区域对所述遥感图像中进行分别提取全陆地区域、全海洋区域、需进行精细海陆分割处理的海岸线附近区域和海岛区域。具体来讲，如图7所示,该步骤又包括以下5个步骤：

S331：对遥感图像进行分块，获取若干图像子块。如将图像均分为512ⅹ512大小的图像子块。

S332：依次读取每一个图像子块；

S333：根据全陆地区域判别标准，判断其是否为全陆地区域子块，如属于全陆地区域子块，则提取该图像子块放入全陆地训练集文件夹中，并跳转至步骤S332继续执行。如不属于全陆地区域子块，进入步骤S334。

S334：根据全海洋区域判别标准，判断其是否为全海洋区域子块，如属于全海洋区域子块，则提取该图像子块放入全海洋训练集文件夹中，并跳转至步骤S332继续执行。如不属于全海洋区域子块，进入步骤S335。

S335：根据海岸线附近区域和海岛区域判别标准，判断其是否为海岸线附近区域或海岛区域子块，如属于该类，则提取该图像子块放入测试集文件夹中，并跳转至步骤S332继续执行。如不属于该类别，不进行其他操作，直接进入步骤S332。

本发明实施例中，选取的图像子块大小均为512ⅹ512。设图像子块为S_abcd，则有：

(1)对全陆地区域子块的判别：若该图像子块均位于海陆粗分割识别出来的遥感图像有效陆地区域中，即图6中的纯白色区域，也即

$S_{\mathrm{abcd}}{\⋐S}_{\mathrm{truelandmask}}---\left(1\right)$ 则该图像子块为全陆地区域子块；

(2)对全海洋区域子块的判别：若该图像子块中不包含任何海陆粗分割识别出来的遥感图像有效陆地区域(即图6中的纯白色区域)的部分，并且该图像子块处于遥感图像中的有效影像区域S_p1p2p3p4中，也即

则该图像子块为全海洋区域子块；

(3)对海岸线边缘区域子块和海岛区域子块的判别：

(3a)若该图像子块一部分为纯白色区域，剩余的部分为非纯白色区域，即

则该图像子块为海岸线附近区域子块或海岛区域子块的候选子块；

(3b)位于有效影像区域S_p1p2p3p4边界上符合条件a的图像子块，如果该图像子块中的有效区域已经全部由海陆粗分割识别出来的区域完全覆盖，则排除该图像子块；

(3c)如果该图像子块识别出的海岛面积，即纯白色区域的面积不足图像子块面积的5％，排除该图像子块。

图8为本发明实施例提供的一个自动提取相关区域子块的结果图，其中粗黑线矩形框表示该图像子块为全陆地区域子块，白色线框表示该图像子块为全海洋区域子块，黑色线框表示该图像子块为需进一步进行海陆精细分割的区域。

S4：对所述训练集子图像和所述测试集子图像进行计算，分别获取这些图像中每个像素点对应的局部信息熵。其中，所述的某个像素所对应的局部信息熵是指，选取以该像素为中心坐标的一个邻域，计算出的该邻域的局部信息熵；

信息熵是对信息概率分布的一种度量。对于灰度图像，如果一幅图像包含G个灰度级，并且各个灰度级出现的概率分别为P₁,P₂,…P_G，则图像的信息熵定义为：

$E=-\underset{k=1}{\overset{G}{\Σ}}{P}_k{\log }_2\left(P_k\right)---\left(4\right)$

在一般灰度图像中G＝255。

图像的局部信息熵是指像素在局部(如mⅹm)邻域中的信息熵，即

$E\left({\mathrm{local}}_{m\×m}\right)=-\underset{k=1}{\overset{G}{\Σ}}{P}_k{\log }_2\left(P_k\right)---\left(5\right)$

其中，local_mⅹm表示局部邻域图像，

$P_k=\frac{n_k}{m\×m}---\left(6\right)$

n_k表示灰度值为k的像素点在局部邻域中的个数。其中，本发明中所述的某个像素所对应的局部信息熵是指，选取以该像素为中心坐标的mⅹm的邻域，计算出的该邻域的局部信息熵。在后面的叙述中不再赘述。

在本发明实施例的实施过程中，对训练集和测试集子图像提取的局部信息熵特征是在像素的9ⅹ9(即m＝9)局部邻域中计算的。

S5：以所有全陆地区域块的子图像计算出的所有局部信息熵作为陆地区域特征集，以所有全海洋区域块的子图像计算出的所有局部信息熵为海洋区域特征集，分别对陆地区域特征集、海洋区域特征集进行概率分布统计，并对其概率分布进行学习，计算出和该概率分布最吻合的分布类型及相关参数，确定最优的特征概率分布模型并获得特征概率分布函数。

图9为从图8中提取的部分全陆地区域和全海洋区域，全陆地区域训练集和全海洋区域训练集分别对应的特征集的概率密度曲线如图10所示。

根据训练集特征的统计分布和概率密度曲线对训练集特征进行模型概率密度函数及参数进行估计，自动选取最适合的分布模型。根据特征数据的分布特点，建立不同类型的分布模型，本发明中使用的概率分布模型有贝塔分布、Birnbaunm-Saunders分布、指数分布、Extremevalue分布、伽马分布、广义极值分布、广义Pareto分布、逆高斯分布、Logistic分布、Log-logistic分布、对数正态分布、Nakagami分布、正态分布、瑞利分布、Rician分布、t分布、韦伯分布等，依据统计学习中的贝叶斯信息准则(BIC准则)选择与训练集特征概率分布最吻合的最佳分布模型。BIC的计算方法如下：

$\mathrm{BIC}=-2\ln \hat{L}+k\ln \left(n\right)$

其中

$\hat{L}=p\left(x\right|\widehat{\mathrm{\θ}},M)$

在上式中，x为待估计数据，n为待估计数据的个数，M为待评估的分布模型，θ为该分布模型对应的参数，则表示该模型对应似然函数的最大值，表示模型似然函数取最大值时的参数取值。

若某一分布计算出来的BIC值越低，表示该分布模型与实际特征集概率分布数据拟合效果越好。

假设对图8所对应的训练集特征进行测试，图11为本发明实施例提供的两类训练集特征概率分布评估结果。其中，全陆地区域特征集服从韦伯分布时，其BIC值最小，因此该训练集特征的最佳概率分布为韦伯分布；而全海洋区域特征集服从广义极值分布时，其BIC值最小，因此该训练集特征的最佳概率分布为广义极值分布。

因此，全陆地区域训练集特征概率分布函数为——韦伯分布，韦伯分布函为：

$f(x;\mathrm{\&lambda;},k)=\left\{\begin{array}{cc}\frac{k}{\mathrm{\&lambda;}}{\left(\frac{x}{\mathrm{\&lambda;}}\right)}^{k-1}{e}^{-{(x/\mathrm{\&lambda;})}^k}& x\&GreaterEqual;0\\ 0& x<0\end{array}\right.---\left(7\right)$

其中x为随机变量，这里表示局部信息熵特征向量；λ>0为比例参数；k>0为形状参数。经参数估计，上述全陆地训练集特征对应的韦伯分布中，参数λ＝5.2459，k＝17.3460，函数曲线对应图11a中的黑色实线。

全海洋区域训练集特征概率分布函数为——广义极值分布，广义极值分布函数表示为：

$f(x;\mathrm{\&mu;},\mathrm{\&sigma;},\mathrm{\&xi;})=\frac{1}{\mathrm{\&sigma;}}t{\left(x\right)}^{\mathrm{\&xi;}+1}{e}^{-t\left(x\right)}---\left(8\right)$

其中，

$t\left(x\right)=\{\begin{array}{cc}{(1+(\frac{x-\mathrm{\&mu;}}{\mathrm{\&sigma;}}\left)\mathrm{\&xi;}\right)}^{-1/\mathrm{\&xi;}}& \mathrm{\&xi;}\&NotEqual;0\\ e^{-(x-\mathrm{\&mu;})/\mathrm{\&sigma;}}& \mathrm{\&xi;}=0\end{array};$

μ为位置参数；σ＞0为比例参数；ξ∈R为形状参数。经参数估计，上述全海洋训练集特征对应的广义极值分布中，μ＝2.3473,σ＝0.2369,ξ＝0.0569，函数曲线对应图11b中的黑色实线。

对大量不同的遥感图像进行测试发现，所有的测试图像数据中提取的全陆地区域、全海洋区域局部信息熵特征概率统计分布均满足上述16种分布中的一种或几种，均可以利用上述方法求出对应的最佳概率分布函数估计。

S6：根据所述特征概率分布函数对所述测试集子图像进行基于局部信息熵的阈值分割，得到阈值分割图。

由上述，假设估计出的全陆地区域训练集的特征概率分布函数如公式(7)所示，代入上述估计出的参数后具体为：

$f_1\left(x\right)=\frac{17.3460}{5.2459}{\left(\frac{x}{5.2479}\right)}^{17.3460-1}{e}^{-{(x/5.2479)}^{17.3460}}---\left(9\right)$

并假设估计出的全海洋区域训练集的特征概率分布函数如公式(8)所示，代入上述估计出的参数后具体为：：

$f_2\left(x\right)=\frac{1}{0.2369}{(1+(\frac{x-2.3473}{0.2369}\left)0.0569\right)}^{-1/0.569(0.569+1)}{e}^{-{(1+(\frac{x-2.3473}{0.2369}\left)0.0569\right)}^{-1/0.0569}}---\left(10\right)$

对于测试集中的某一子图像，提取每一个像素的局部信息熵特征，设A(m,n)为某一像素点，A_x(m,n)为该像素点对应的局部信息熵特征，则有判别函数

$\left\{\begin{array}{cc}{f}_1\left(A_x\right(m,n))\&GreaterEqual;f_2\left(A_x\right(m,n)),& A\&Element;land\\ f_1\left(A_x\right(m,n))<f_2\left(A_x\right(m,n)),& A\&Element;ocean\end{array}\right.---\left(11\right)$

因为全陆地曲线函数和全海洋曲线函数之间分离度较大，为了降低时间开销，可以采用基于贝叶斯原理的阈值判别方式代替上述判别函数。

由

f₁(x)＝f₂(x)(12)解得两个函数曲线相交时的交点x作为阈值T,即T＝x，在上述两个具体函数中，T＝x＝3.9149，则有判别函数

$\left\{\begin{array}{cc}{A}_x(m,n)\&GreaterEqual;T,& A\&Element;land\\ A_x(m,n)<T,& A\&Element;ocean\end{array}\right.---\left(13\right)$

以此阈值对从图8中提取的所有测试集图像进行基于局部信息熵的阈值分割，即海陆精细分离，得到阈值分割图。

S7：对所述阈值分割图进行形态学处理，得到海陆精细分割结果图。

本步骤主要包括两个步骤，如图12所示：

S71：对于阈值分割图进行空洞和细缝填充，获得初级海陆精细分割结果图。在本步骤中，对于阈值分割图，在陆地区域可能会存在因为植被、湖水、河流等因素造成的孔洞或细缝，孔洞或细缝的存在会对海陆分割结果造成误判，可以利用形态学中的闭运算和孔洞填充的方法对阈值分割图进行连接和填补，消除孔洞与细缝。

S72：去除所述初级海陆精细分割结果图中面积较小的连通域，得到海陆精细分割结果图。在本步骤中，诸如舰船目标的近岸海面一般目标较小，为了尽量降低将在海岸线附近的海面目标判别为陆地区域，在对图像进行孔洞填充后还需去除面积较小的连通域，这些面积较小的连通域作为疑似候选海面目标需要得到保留。在本步骤中，连通域面积阈值大小的选取需根据足够的先验知识和舰船等目标外形的实际参数来进行确定。

C，粗分割与精细分割结果融合

该过程对应图2中的S8。如图13所示，分为3个步骤：

S81：将所述海陆精细分割结果图影射融合至所述海陆粗分割结果图中的相应位置，得到二值化掩模图像。在本步骤中，在粗分割的基础上，将测试集图像即海岸线边缘区域、海岛区域图像精细分割的结果映射至粗分割结果图像中相应位置，用这些图像子区域精细分割的结果代替粗分割的结果，得到二值化掩模图像。

S82：对所述二值化掩模图像进行形态学处理，得到海陆分割图像。基于统计学习的海陆精细分割后的二值化掩模图像虽然在输出前已经用形态学的相关操作对孔洞、细缝进行了填充，但由于部分测试集海陆分割结果图在上、下、左、右边界处的孔洞并没有完全被填充，在将精细化分割后的结果与粗分割二值化处理后的结果进行融合后不可避免地会在融合的边缘区域会产生孔洞，因此在融合后仍需对融合后的二值化图像进行形态学操作，填充其中的孔洞。

S83：将所述海陆分割图像影射至所述遥感图像，得到海陆分割结果图，所述海陆分割结果图包括：

大幅面遥感图像陆地区域——S_im·B_im；

大幅面遥感图像海洋区域——S_im·(1-B_im)；

其中S_im为原始遥感图像，B_im为二值化的海陆分割掩模图像。

本发明中，涉及到形态学中的相关操作，如闭运算以及孔洞填充等，以及必要的图像区域填充方法等，对图像处理和计算机视觉领域的技术人员是已知内容，此处不再赘述。

如图14所示，为本发明实施例提供的一种大幅面遥感图像海陆分割的系统，包括：

海陆粗分割模块10，用于对所述遥感图像进行海陆填充，获取海陆粗分割结果图；

海陆精细分割模块11，用于利用所述海陆粗分割结果图对所述遥感图像进行分块提取，获取海陆精细分割结果图；

海陆分割结果获取模块12，用于将所述海陆粗分割结果图与所述海陆精细分割结果图进行融合，获取海陆分割结果图。

如图15所示，海陆粗分割模块10具体包括：

海岸线添加模块101，用于利用GSHHS数据为所述遥感图像添加海岸线，生成海岸线添加图，具体的，海岸线添加模块具体用于：

a11，输入某一遥感图像，读取与所述遥感图像相对应的地理信息文件，获取所述遥感图像4个顶点对应的经纬度坐标；

a12,根据所述遥感图像4个顶点对应的经纬度坐标信息，确定一个四边形的经纬度区域；从GSHHS数据库中读取位于所述经纬度区域范围内的所有离散海岸线经纬度数据；

a13，根据所述遥感图像的分辨率大小获取4个顶点A、B、C、D对应的顶点像素坐标，并计算所述顶点像素坐标与所述经纬度坐标之间的线性关系；

a14，将位于所述经纬度区域范围内的离散海岸线经纬度数据映射为对应的遥感图像像素坐标，得到离散海岸线像素坐标；

a15，在遥感图像中，将所述离散海岸线像素坐标连成线，形成海岸线，生成海岸线添加图。

海陆填充模块102，对所述海岸线添加图进行海陆填充，获取海陆粗分割结果图。具体的，海陆填充模块102具体用于：

a21，读取所述经纬度坐标，根据所述经纬度坐标确定一个四边形经纬度区域；

a22，依次读取GSHHS数据库中的经纬度数据段，每一个经纬度数据段为一封闭陆地或海岛经纬度数据；

a23，判断所述经纬度数据段与所述四边形经纬度区域是否有重合区域；

a24，若所述经纬度数据段与所述四边形经纬度范围有重合区域，则保留重合区域范围的经纬度数据；

a25，判断步骤a24获取的经纬度数据是否为一封闭数据段；若该经纬度数据为一非封闭数据段数据，则将所述非封闭数据段数据重新构造成一个封闭区域，并将所述非封闭数据段数据与其他数据段数据隔开；

a26，若判断步骤a24获取的经纬度数据为一封闭数据段，则将所述经纬度数据变换为所述海岸线添加图的像素坐标，并将位于这些像素坐标围成的曲线范围内的像素点标记成陆地区域或海岛区域；

a27，利用区域填充算法，将步骤a26中获取的区域内部填充成一个完整的陆地块或海岛块；

a28，判断GSHHS数据库中是否仍存在其他封闭数据段，若存在重复步骤a22至步骤a27直至海陆填充完毕，得到海陆粗分割结果图。

海陆精细分割单元11具体包括：

图像分块模块111,用于对所述遥感图像进行分块，得到全陆地区域块、全海洋区域块、海岸线边缘区域块和海岛区域块；以所述全陆地区域块和所述全海洋区域块为训练集，以所述海岸线边缘区域块和所述海岛区域块为测试集；分别对所述训练集和所述测试集提取子图像，获取训练集子图像和测试集子图像。具体的，图像分块模块111具体用于：

b11，确定所述遥感图像的有效影像区域；

b12，利用所述有效影像区域确定所述海陆粗分割结果图中的有效陆地区域；

b13，根据所述有效陆地区域对所述遥感图像进行分块提取，得到所述全陆地区域块、全海洋区域块、海岸线边缘区域块和海岛区域块。在本步骤中，图像分块模块111还用于：

b131,对所述遥感图像进行分块，获取若干图像子块；

b132，依次读取每一图像子块；

b133，判断所述图像子块是否为全陆地区域子块；

若所述图像子块为全陆地区域子块，则将所述图像子块放入全陆地训练集文件夹中，并继续执行步骤b132；

若所述图像子块非全陆地区域子块，则继续进行步骤b134；

b134，判断所述图像子块是否为全海洋区域子块；

若所述图像子块为全海洋区域子块，则将所述图像子块放入全海洋训练集文件夹中，并继续执行步骤b132；

若所述图像子块非全海洋区域子块，则继续进行步骤b135；

b135，判断所述图像子块是否为海岸线边缘区域子块或海岛区域子块；

若所述图像子块为海岸线边缘区域子块或海岛区域子块，则将所述图像子块放入海岸线边缘区域或海岛区域训练集文件夹中，并继续执行步骤b132；

若所述图像子块非海岸线边缘区域子块或海岛区域子块，则继续进行步骤b132。

局部信息熵获取模块112，用于对所述训练集子图像和所述测试集子图像进行计算，分别获取每一个像素点对应的局部信息熵；

分布函数获取模块113，用于以根据全陆地区域块的子图像计算出的局部信息熵作为陆地区域特征集，以根据全海洋区域块的子图像计算出的局部信息熵为海洋区域特征集，分别对所述陆地区域特征集和海洋区域特征集进行概率分布统计及学习，计算出与所述概率分布最吻合的分布类型及相关参数，确定最优特征概率分布模型并获得相应的特征概率分布函数；

阈值分割模块114，用于根据所述特征概率分布函数对所述测试集子图像进行基于局部信息熵的阈值分割，得到阈值分割图；

精细分割获取模块115，用于对所述阈值分割图进行形态学处理，得到海陆精细分割结果图。具体的，精细分割获取模块115具体用于：

b51，对于阈值分割图进行空洞和细缝填充，获得初级海陆精细分割结果图；

b52，去除所述初级海陆精细分割结果图中面积较小的连通域，得到海陆精细分割结果图。

海陆分割结果获取单元12包括：

影射融合模块121，将所述海陆精细分割结果图影射融合至所述海陆粗分割结果图中的相应位置，得到二值化掩模图像；

处理模块122，对所述二值化掩模图像进行形态学处理，得到海陆分割图像；

图像获得模块123，将所述海陆分割图像影射至所述遥感图像，得到海陆分割结果图；所述海陆分割结果图包括大幅面遥感图像陆地区域和大幅面遥感图像海洋区域。

本发明提出的大幅面遥感图像的海陆分割方法，首先利用遥感图像的经纬度信息，在全球海岸线数据库(GSHHS)的基础上填充大部分陆地区域，实现大幅面遥感图像的海陆粗分割；根据海陆粗分割结果，提取需进一步精细分割处理的区域子图像作为测试集，提取全陆地区域和全海洋区域子图像作为训练集，对训练集图像建立基于局部信息熵特征的自适应概率统计学习模型，确定最优的概率分布模型，利用该模型对测试集子图像进行基于局部信息熵的贝叶斯海陆精细分割。将基于地理信息的海陆粗分割结果和海陆精细分割结果融合在一起，并对融合结果进行形态学处理，最终得到完整的大幅面遥感图像海陆分割图像。相对于现有的海陆分割方法，该方法可以有效处理复杂背景情形下的大幅面遥感图像的海陆分割问题，分割结果更为精准，对基于海陆分离的相关技术研究如海洋目标检测影响较小，适用于对遥感图像进行批量处理，实用性更强。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来控制相关的硬件完成，所述的程序可以在存储于一计算机可读取存储介质中，所述的存储介质，如ROM/RAM、磁盘、光盘等。

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (18)

1.一种大幅面遥感图像海陆分割的方法，其特征在于，所述方法包括下述步骤：

步骤A，对所述遥感图像进行海陆填充，获取海陆粗分割结果图；

步骤B，结合所述海陆粗分割结果图对所述遥感图像进行分块提取，获取海陆精细分割结果图；

步骤C，将所述海陆粗分割结果图与所述海陆精细分割结果图进行融合，获取海陆分割结果图。

2.如权利要求1所述的方法，其特征在于，所述步骤A包括：

步骤A1，利用GSHHS数据库为所述遥感图像添加海岸线，生成海岸线添加图；

步骤A2，对所述海岸线添加图进行海陆填充，获取海陆粗分割结果图。

3.如权利要求2所述的方法，其特征在于，所述步骤A1包括：

步骤A11，输入某一遥感图像，读取与所述遥感图像相对应的地理信息文件，获取所述遥感图像4个顶点对应的经纬度坐标；

步骤A12,根据所述遥感图像4个顶点对应的经纬度坐标信息，确定一个四边形的经纬度区域；从所述GSHHS数据库中读取位于所述经纬度区域范围内的所有离散海岸线经纬度数据；

步骤A13，根据所述遥感图像的分辨率大小获取4个顶点对应的顶点像素坐标，并计算每个顶点的像素坐标与各自经纬度坐标之间的线性关系，根据所述线性关系获取所述遥感图像中每个像素点对应的经纬度数据；

步骤A14，将位于所述经纬度区域内的离散海岸线经纬度数据映射为对应的遥感图像的像素坐标，得到离散海岸线像素坐标；

步骤A15，在所述遥感图像中，将所述离散海岸线像素坐标连成线，形成海岸线，生成海岸线添加图。

4.如权利要求3所述的方法，其特征在于，所述步骤A2包括：

步骤A21，读取所述经纬度坐标，根据所述经纬度坐标确定一个四边形经纬度区域；

步骤A22，依次读取GSHHS数据库中的经纬度数据段，每一个经纬度数据段为一封闭陆地或海岛经纬度数据；

步骤A23，判断所述经纬度数据段与所述四边形经纬度区域是否有重合区域；

步骤A24，若所述经纬度数据段与所述四边形经纬度范围有重合区域，则保留重合区域范围的经纬度数据；

步骤A25，判断步骤A24获取的经纬度数据是否为一封闭数据段；若该经纬度数据为一非封闭数据段数据，则将所述非封闭数据段数据重新构造成一个封闭区域，并将所述非封闭数据段数据与其他数据段数据隔开；

步骤A26，若判断步骤A24获取的经纬度数据为一封闭数据段，则将所述经纬度数据变换为所述海岸线添加图的像素坐标，并将位于这些像素坐标围成的曲线范围内的像素点标记成陆地区域或海岛区域；

步骤A27，利用区域填充算法，将步骤A26中获取的区域填充成一个完整的陆地块或海岛块；

步骤A28，判断GSHHS数据库中是否仍存在其他封闭数据段，若存在重复步骤A22至步骤A27直至海陆填充完毕，得到海陆粗分割结果图。

5.如权利要求1所述的方法，其特征在于，所述步骤B包括：

步骤B1，对所述遥感图像进行分块，得到全陆地区域块、全海洋区域块、海岸线边缘区域块和海岛区域块；以所述全陆地区域块和所述全海洋区域块为训练集，以所述海岸线边缘区域块和所述海岛区域块为测试集；

分别对所述训练集和所述测试集提取子图像，获取训练集子图像和测试集子图像；

步骤B2，对所述训练集子图像和所述测试集子图像进行计算，分别获取这些图像中每个像素点对应的局部信息熵；

步骤B3，以根据全陆地区域块的子图像计算出的局部信息熵作为陆地区域特征集，以根据全海洋区域块的子图像计算出的局部信息熵为海洋区域特征集，分别对所述陆地区域特征集和海洋区域特征集进行概率分布统计及学习，计算出与所述概率分布最吻合的分布类型及相关参数，确定最优特征概率分布模型并获得相应的特征概率分布函数；

步骤B4，根据所述特征概率分布函数对所述测试集子图像进行基于局部信息熵的阈值分割，得到阈值分割图；

步骤B5，对所述阈值分割图进行形态学处理，得到海陆精细分割结果图。

6.如权利要求5所述的方法，其特征在于，所述步骤B1包括：

步骤B11，确定所述遥感图像的有效影像区域；

步骤B12，利用所述有效影像区域确定所述海陆粗分割结果图中的有效陆地区域；

步骤B13，根据所述有效陆地区域对所述遥感图像中进行分块提取，得到所述全陆地区域块、全海洋区域块、海岸线边缘区域块和海岛区域块。

7.如权利要求6所述的方法，其特征在于，所述步骤B13具体包括：

步骤B131,对所述遥感图像进行分块，获取若干图像子块；

步骤B132，依次读取每一图像子块；

步骤B133，判断所述图像子块是否为全陆地区域子块；

若所述图像子块为全陆地区域子块，则将所述图像子块放入全陆地训练集文件夹中，并继续执行步骤B132；

若所述图像子块非全陆地区域子块，则继续进行步骤B134；

步骤B134，判断所述图像子块是否为全海洋区域子块；

若所述图像子块为全海洋区域子块，则将所述图像子块放入全海洋训练集文件夹中，并继续执行步骤B132；

若所述图像子块非全海洋区域子块，则继续进行步骤B135；

步骤B135，判断所述图像子块是否为海岸线边缘区域子块或海岛区域子块；

若所述图像子块为海岸线边缘区域子块或海岛区域子块，则将所述图像子块放入海岸线边缘区域或海岛区域测试集文件夹中，并继续执行步骤B132；

若所述图像子块非海岸线边缘区域子块或海岛区域子块，则继续进行步骤B132。

8.如权利要求5所述的方法，其特征在于，所述步骤B5包括：

步骤B51，对于阈值分割图进行空洞和细缝填充，获得初级海陆精细分割结果图；

步骤B52，去除所述初级海陆精细分割结果图中面积较小的连通域，得到海陆精细分割结果图。

9.如权利要求1所述的方法，其特征在于，所述步骤C包括：

步骤C1，将所述海陆精细分割结果图影射融合至所述海陆粗分割结果图中的相应位置，得到二值化掩模图像；

步骤C2，对所述二值化掩模图像进行形态学处理，得到海陆分割图像；

步骤C3，将所述海陆分割图像影射至所述遥感图像，得到海陆分割结果图；所述海陆分割结果图包括大幅面遥感图像陆地区域和大幅面遥感图像海洋区域。

10.一种大幅面遥感图像海陆分割的系统，其特征在于，所述系统包括：

海陆粗分割单元，用于对所述遥感图像进行海陆填充，获取海陆粗分割结果图；

海陆精细分割单元，用于利用所述海陆粗分割结果图对所述遥感图像进行分块提取，获取海陆精细分割结果图；

海陆分割结果获取单元，用于将所述海陆粗分割结果图与所述海陆精细分割结果图进行融合，获取海陆分割结果图。

11.如权利要求10所述的系统，其特征在于，所述海陆粗分割单元包括：

海岸线添加模块，用于利用GSHHS数据为所述遥感图像添加海岸线，生成海岸线添加图；

海陆填充模块，对所述海岸线添加图进行海陆填充，获取海陆粗分割结果图。

12.如权利要求11所述的系统，其特征在于，所述海岸线添加模块具体用于：

首先，输入某一遥感图像，读取与所述遥感图像相对应的地理信息文件，获取所述遥感图像4个顶点对应的经纬度坐标；

接着,根据所述遥感图像4个顶点对应的经纬度坐标信息，确定一个四边形的经纬度区域；从GSHHS数据库中读取位于所述经纬度区域范围内的所有离散海岸线经纬度数据；

接着，根据所述遥感图像的分辨率大小获取4个顶点对应的顶点像素坐标，并计算每个顶点的像素坐标与各自经纬度坐标之间的线性关系，根据所述线性关系获取所述遥感图像中每个像素点对应的经纬度数据；

接着，将位于所述经纬度区域范围内的离散海岸线经纬度数据映射为对应的遥感图像像素坐标，得到离散海岸线像素坐标；

最后，在遥感图像中，将所述离散海岸线像素坐标连成线，形成海岸线，生成海岸线添加图。

13.如权利要求11所述的系统，其特征在于，所述海陆填充模块具体用于：

a21，读取所述经纬度坐标，根据所述经纬度坐标确定一个四边形经纬度区域；

a22，依次读取GSHHS数据库中的经纬度数据段，每一个经纬度数据段为一封闭陆地或海岛经纬度数据；

a23，判断所述经纬度数据段与所述四边形经纬度区域是否有重合区域；

a24，若所述经纬度数据段与所述四边形经纬度范围有重合区域，则保留重合区域范围的经纬度数据；

a25，判断步骤a24获取的经纬度数据是否为一封闭数据段；若该经纬度数据为一非封闭数据段数据，则将所述非封闭数据段数据重新构造成一个封闭区域，并将所述非封闭数据段数据与其他数据段数据隔开；

a26，若判断步骤a24获取的经纬度数据为一封闭数据段，则将所述经纬度数据变换为所述海岸线添加图的像素坐标，并将位于这些像素坐标围成的曲线范围内的像素点标记成陆地区域或海岛区域；

a27，利用区域填充算法，将步骤a26中获取的区域填充成一个完整的陆地块或海岛块；

a28，判断GSHHS数据库中是否仍存在其他封闭数据段，若存在重复步骤a22至步骤a27直至海陆填充完毕，得到海陆粗分割结果图。

14.如权利要求10所述的系统，其特征在于，所述海陆精细分割单元包括：

图像分块模块,用于对所述遥感图像进行分块，得到全陆地区域块、全海洋区域块、海岸线边缘区域块和海岛区域块；以所述全陆地区域块和所述全海洋区域块为训练集，以所述海岸线边缘区域块和所述海岛区域块为测试集；

分别对所述训练集和所述测试集提取子图像，获取训练集子图像和测试集子图像；

局部信息熵获取模块，用于对所述训练集子图像和所述测试集子图像进行计算，分别获取这些图像中每个像素点对应的局部信息熵；

分布函数获取模块，用于以根据全陆地区域块的子图像计算出的每个像素点的局部信息熵作为陆地区域特征集，以根据全海洋区域块的子图像计算出的每个像素点的局部信息熵为海洋区域特征集，分别对所述陆地区域特征集和海洋区域特征集进行概率分布统计及学习，计算出与所述概率分布最吻合的分布类型及相关参数，确定最优特征概率分布模型并获得相应的特征概率分布函数；

阈值分割模块，用于根据所述特征概率分布函数对所述测试集子图像进行基于局部信息熵的阈值分割，得到阈值分割图；

精细分割获取模块，用于对所述阈值分割图进行形态学处理，得到海陆精细分割结果图。

15.如权利要求14所述的系统，其特征在于，所述图像分块模块具体用于：

b11，确定所述遥感图像的有效影像区域；

b12，利用所述有效影像区域确定所述海陆粗分割结果图中的有效陆地区域；

b13，根据所述有效陆地区域对所述遥感图像中进行分块提取，得到所述全陆地区域块、全海洋区域块、海岸线边缘区域块和海岛区域块。

16.如权利要求15所述的系统，其特征在于，在所述b13中，所述图像分块模块还用于：

b131,对所述遥感图像进行分块，获取若干图像子块；

b132，依次读取每一图像子块；

b133，判断所述图像子块是否为全陆地区域子块；

若所述图像子块为全陆地区域子块，则将所述图像子块放入全陆地训练集文件夹中，并继续执行步骤b132；

若所述图像子块非全陆地区域子块，则继续进行步骤b134；

b134，判断所述图像子块是否为全海洋区域子块；

若所述图像子块为全海洋区域子块，则将所述图像子块放入全海洋训练集文件夹中，并继续执行步骤b132；

若所述图像子块非全海洋区域子块，则继续进行步骤b135；

b135，判断所述图像子块是否为海岸线边缘区域子块或海岛区域子块；

若所述图像子块为海岸线边缘区域子块或海岛区域子块，则将所述图像子块放入海岸线边缘区域或海岛区域测试集文件夹中，并继续执行步骤b132；

若所述图像子块非海岸线边缘区域子块或海岛区域子块，则继续进行步骤b132。

17.如权利要求14所述的系统，其特征在于，所述精细分割获取模块具体用于：

b51，对于阈值分割图进行空洞和细缝填充，获得初级海陆精细分割结果图；

b52，去除所述初级海陆精细分割结果图中面积较小的连通域，得到海陆精细分割结果图。

18.如权利要求10所述的系统，其特征在于，所述海陆分割结果获取单元包括：

影射融合模块，将所述海陆精细分割结果图影射融合至所述海陆粗分割结果图中的相应位置，得到二值化掩模图像；

处理模块，对所述二值化掩模图像进行形态学处理，得到海陆分割图像；

图像获得模块，将所述海陆分割图像影射至所述遥感图像，得到海陆分割结果图；所述海陆分割结果图包括大幅面遥感图像陆地区域和大幅面遥感图像海洋区域。