技术摘要：
本发明提供了一种遥感图像海陆自动分割方法，包括如下步骤：1、对待分割遥感图像进行线性增强处理，改善图像对比度；2、广义径向高斯核函数改进FLICM聚类方法，得到陆地初始区域，结合Canny算子提取海岸线粗略轮廓；3、根据粗分割结果所得区域及轮廓作为加权CV模型初始  全部
背景技术：
通过海陆分割可以实现海洋与陆地分离，抑制陆地区域复杂背景的干扰，有效降 低海面目标检测的难度。海陆分割作为海面目标检测与识别的前提与基础，分割结果的优 劣直接影响后续处理的准确性，在民用、军事侦察、海洋监控等领域具有重要意义。另外，通 过海陆分割获得的海岸线信息在潮汐观测、地图绘制、海洋资源管理等诸多应用场景起着 至关重要的作用。因此，亟需一种有效的海陆分割方法。 传统海陆分割普遍是基于海图、剖面监测、地形图等方法，耗费大量时间、人力，且 动态跟踪性差。随着遥感图像被人们所关注，基于遥感图像的海陆分割方法逐渐成为研究 热点。目前基于遥感图像的海陆分割方法主要分为基于区域信息及基于边界轮廓等。其中 模糊C均值(FCM)聚类方法作为一种有效的基于区域信息的图像分割方法，被不断改进。近 年提出的模糊局部信息C均值(FLICM)聚类方法在像素模糊隶属度判定过程中综合考虑了 像素局部的空间和灰度信息，展现出较强的鲁棒性，但未考虑使用核函数映射方式，且其仅 依据图像区域灰度对图像进行分割，未考虑图像全局信息，分割精度有待提高。另一方面， 基于无边缘主动轮廓模型Chan-Vese(CV)模型的方法利用图像全局信息建立能量函数，以 此控制曲线运动，分割精度较高但计算效率较低，且仅利用区域的均值信息确定边界，对初 始条件敏感。
技术实现要素：
发明目的：针对现有遥感图像海陆分割领域中存在的问题，进一步提高遥感图像 海陆分割的精度和速度，本发明公开了一种基于改进FLICM和加权CV的遥感图像海陆自动 分割方法。针对遥感图像对比度不理想，主观视觉效果较差的问题，首先对遥感图像线性拉 伸，突出感兴趣的信息；然后采用改进FLICM聚类方法对遥感图像进行海陆粗分割，引入广 义径向高斯核函数使像素点在高维特征空间中拥有更优的线性可聚性；为了避免CV模型对 初始条件敏感这个缺陷，再分别将粗分割得到的陆地区域和海岸线粗略轮廓作为加权CV模 型的初始迭代条件，实现海陆精细分割；最后，综合考虑陆地与海洋的不同特性，利用形态 学手段及连通域准则在进一步保持海岸线完整性的同时消除散杂点，得到最终的海陆分割 结果。 技术方案：本发明公开了一种基于改进FLICM和加权CV的遥感图像海陆自动分割 方法，包括如下步骤： 包括如下步骤： 步骤1，对待分割遥感图像进行线性变换增强处理； 步骤2，引入广义径向高斯核函数对FLICM聚类方法进行改进，实现遥感图像海陆 粗分割，结合Canny算法提取海岸线粗略轮廓； 5 CN 111582198 A 说　明　书 2/7 页 步骤3，计算得到遥感图像海陆精细分割结果； 步骤4，利用形态学手段及连通域准则对海陆精细分割结果进一步修正，最终完成 海陆分割。 步骤1包括：利用如下线性变换公式对待分割遥感图像进行拉伸处理： 式中，ge为增强后图像灰度级，g为待分割遥感图像的原始图像灰度级，gmax、gmin分 别为原始图像灰度级的最大值与最小值。 步骤2包括： 步骤2-1，初始化核模糊局部信息C均值聚类参数：设定待聚类图像像素的个数n、 聚类数目c、隶属度指数m、最大迭代次数TM、停止迭代条件ε＞0，初始化迭代次数t＝0、原始 隶属度uik和聚类中心{v1,v2,...,vc}，vc表示第c个聚类中心； 步骤2-2，引入广义径向高斯核函数，利用下式计算每个样本xk的原始隶属度uik： 其中，i、j的取值范围为1～c；Kugrb为广义径向高斯核函数： 式中，b、r、ρ为可调参数，b＞0,b∈z，z为整数集合；r≥1，P为核参数，D为数据的维 数；x、y分别表示图像上像素点的横坐标和纵坐标； 步骤2-3，利用下式计算并更新每个样本xk的新的模糊隶属度 式中，xk隶属于第i类区域的隶属度为uik；xk隶属于第j类区域的隶属度为ujk； p、q为决定原有隶属度函数和空间函数相关性的指数参数； 空间函数hik表示由局部空间像素决定的xk属于第i类区域的可能性； 空间函数hjk表示由局部空间像素决定的xk属于第j类区域的可能性； uit为像素xt属于第i类局部区域的隶属度，变量t属于Nk，Nk为以xk 为中心的局部区域； 步骤2-4，聚类中心计算与更新： 6 CN 111582198 A 说　明　书 3/7 页 利用下式计算并更新t 1时刻的聚类中心 ： v′i表示t时刻的的聚类中心，u′ik表示更新后的隶属度； 步骤2-5，如果 或t＝TM，TM是设定的阈值，则停止迭代更新，否则t ＝t 1，返回步骤2-2； 步骤2-6，依据最大模糊隶属度判定准则，如果uji＞ujk则将xj划分为第i类区域，对 增强后遥感图像进行聚类分割，得到遥感图像海陆粗分割结果，uji表示隶属于第i类区域的 隶属度，k＝1,2,...,c；i≠k； 步骤2-7，依据步骤2-6粗分割后得到的遥感图像海陆粗分割结果，采用Canny算法 提取海岸线粗略轮廓。 步骤3包括： 步骤3-1，根据步骤2-7所得的海岸线粗略轮廓，初始化零水平集函数 如下式所示： 步骤3-2，步骤2-6所得的海陆粗分割结果表示为改进核模糊局部信息聚类结果g (x,y)，将g(x,y)作为加权CV模型的输入，并根据下式改进偏微分方程进行计算，使能量函 数最小： 表示理想阶跃函数，当进行数值运算时，用下式参与运算： 式中，ε表示趋于0的小正数(一般取值为10-6)； μ、v、λ1、λ2为常数，C1、C2分别为引入自适应权值δ1和δ2的目标拟合中心和背景拟合 中心，定义如下： 7 CN 111582198 A 说　明　书 4/7 页 其中参数 参数 通过迭代实现自适应 调节，海岸线粗略轮廓经过CV模型迭代分割后得到最佳轮廓； 步骤3-3，对于得到的最佳轮廓，通过水平集函数更新及边界演化，实现海陆精细 分割。 步骤3-3中，选择符号距离函数作为水平集函数，则水平集函数变为： 其中，所述轮廓为步骤3-2得到的最佳轮廓，d为符号距离函数，表示高维空间内点 (x，y)到零水平集的距离。 步骤4包括：采用连通域判定准则，取海陆精细分割结果中最大连通域为海洋，通 过形态学开运算方法消除海洋内部可能存在的杂散点及噪声，再通过形态学闭运算方法去 除目标内部孔洞，最终完成海陆自动分割。 有益效果：与现有技术相比，本发明公开的基于改进FLICM和加权CV的遥感图像海 陆自动分割方法的优点在于： (1)有机结合现有遥感图像海陆分割中基于区域信息及基于边界轮廓两类主流算 法，充分考虑图像局部性及全局性，取长补短，提高遥感图像海陆分割的精度和速度。 (2)通过图像增强方法对遥感图像预处理，改善图像整体对比度，提高图像使用价 值，为后续分割提供了良好的数据条件。 (3)引入广义径向高斯核函数改进FLICM聚类方法，将待分割图像的样本特征非线 性映射高维特征空间进行放大，有利于特征的分辨及提取。且广义径向高斯核函数规避了 单一高斯核可调参数唯一的弊端，通过多参数提升学习及泛化能力，以此改善FLICM方法聚 类性能，实现遥感图像海陆粗分割。 (4)将遥感图像核模糊局部信息C均值聚类结果作为加权CV模型的初始区域及轮 廓，解决传统CV模型初始条件未知、收敛时间长的问题，提高方法速度及智能化程度。同时 构建自适应权值加权平均策略替代传统的算术平均计算像素点对拟合中心的贡献值，充分 考虑差异性，使分割结果更加精准，从而实现遥感图像海陆精细分割。 (5)依据不同区域特性采用形态学及连通域准则对精细分割结果修正，进一步保 证分割方法的性能。 附图说明 下面结合附图和