CN104376570A - 一种遥感二值图像目标区域分割的方法 - Google Patents

Abstract

本发明公开了一种遥感二值图像目标区域分割的方法，采用映射表、辅助表“双表”实时记录和修正等价标号，完成图像的最终连通域标记，很好地解决了标记冲突的问题。比传统算法在处理效率上有显著提高，具有较好的应用价值。

Description

一种遥感二值图像目标区域分割的方法

技术领域

本发明涉及一种图像分割技术，特指一种遥感二值图像目标区域分割的方法。 

背景技术

计算机自动解译遥感影像时，地物目标识别的结果常常会以二值图的形式输出以待后续处理。如环境卫星影像提取湖泊水体，当工程任务要考察其中特定的湖泊时（如计算单个湖泊面积），就需要区分出二值图中各个不同的水域，这就必须对生成二值图进行连通区域标记。二值图连通域标记是指将图中满足四邻域或八邻域连通规则的像素识别为同一个目标，然后用唯一的标号来标记连通区域内的像素点。二值图像连通域分割方法，按照处理对象的不同主要分为基于像素和基于线段两大类。像素点扫描方式有区域增长法、顺序扫描法、递归标记法等。线段扫描算法主要有基于游程的标记算法以及多种改进算法，此外还有基于模块的标记算法。各种算法的主要不同点在于对冲突标记对的处理方式，最终的标记要通过算法合并或者回溯扫描，将等价标号标记为同一连通域标号。优化算法的目的在于正确解决好冲突标记的同时，提高标记算法的速度和性能。 

遥感二值图像尺寸大，单个连通区域像素多，基于像素的算法容易导致堆栈溢出或者前景像素多次反复遍历，效率不高。而基于线段的传统标记算法相比基于像素的方法效率要高，利用邻接表来记录冲突标记，一旦遇到等价对就记录在邻接表中，在数据量较多且形状复杂的遥感二值图像处理上，容易使邻接表存储大量的等价对信息，既浪费存储空间也不利于算法合并处理。 

发明内容

本发明针对现有技术存在的上述不足，提供一种遥感二值图像目标区域分割的方法。 

为实现上述目的，本发明采取的技术方案是：一种遥感二值图像目标区域分割的方法，其特征在于，采用以下算法： 

功能：二值图像初步标号，建立标号映射表

输入：遥感二值图像Fij

输出：中间图像Fij’，标号映射表tab_ys

（1）定义映射表tab_ys ，定义辅助表tab_fz ，总标号label = 0 ，映射表和辅助表各元素初始化为0。

（2）从上到下逐行扫描图像，若在第i 行中发现1 值段Fi,m~Fi,n ，为其申请结点，进行（3）；否则返回（2）进行下一行扫描。 

（3）据公式: (L1.begin≤L2.end + 1)∩(L1.end≥L2.begin - 1)，与第i - 1 行中线段结点逐个判断连通性，若都不连通（适合i = 1），则label = label + 1 赋值给F’i,m~F’i,n和当前线段结点,tab_ys[label] = label ，返回（2）；否则有k(k ≥1) 个线段结点与之连通，线段标号为label(1)~label(k) ，用label(1) 赋值给当前线段结点和F’i,m~F’i,n ，进行（4）。 

（4）若k =1，返回（2）；否则逐个检查tab_ys[label(x)](2≤x≤k) 与tab_ys[label(1)] 的关系，若不相等，则执行correct(tab_ys[label(1)],tab_ys[label(x)]) 。全部检查完毕返回（2）。 

根据初步标号结果和映射表，执行Fij' = tab_ysFij' ，即可完成图像的最终连通域标记。 

上述所述的一种遥感二值图像目标区域分割的方法，其特征在于，所述的辅助表负责记录各个连通域对应的全部标号，如连通域j对应的全部标号为{a，b，c} ，则tab_fz[ j] = a ，tab_fz[a] = b ，tab_fz[b] = c ，tab_fz[c] = 0 ，从而形成了一个标号链表。 

上述所述的一种遥感二值图像目标区域分割的方法，其特征在于，所述的correct(tab_ys[label(1)],tab_ys[label(x)])为： 

void correct（int label1，int labelx）{

int tmp=labelx；

while（tmp）{

tab_ys[tmp]=label1；

tmp=tab_fz[tmp]；}

tmp=label1；

while（tab_fz[tmp]）tmp=tab_fz[tmp]；

tab_fz[tmp]=labelx；}

具体实施方式：

下面通过具体实施例，对本发明的技术方案作进一步具体的说明，但是本发明并不限于这些实施例。

如下，表1为二值图像示意图表，表2为中间图像示意图表。 

双表的赋值是在逐行扫描图像，对目标线段进行标号的同时完成的，建立双表的过程如下： 

（1）若当前扫描线段为孤立线段（与上一行无邻接线段），则赋值新标号label 给中间图像的当前位置，如表2中第一行标号为1 的线段，第二行标号为2、3 的线段。执行tab_ys[label] = label 。

（2）当前线段与上一行唯一一条线段连通，则直接复制上一行的标号赋值给当前线段。如表2中第三行标号为1、2、3 的线段。 

（3）当前线段与上一行多于一条线段连通，则将最左边连通线段标号赋值给当前线段，然后找到当前线段标号在映射表中对应的连通域A1 ，第L(L≥2) 次与其连通线段的标号对应连通域为A2 。若A1 和A2 相同，则不作处理；否则，A1 和A2 应实现合并，先将A2 在辅助表中对应的全部标号找出，分别以这些标号作下标，修改其在映射表中对应连通域为A1 ，然后再在辅助表中找到A1 对应的标号链表末端元素，将末端元素更改为A2 。具体两个不同连通区域label1 和label2 连通合并，可用correct 函数来表述。 

void correct（int label1，int label2）{ 

int tmp=label2；

while（tmp）{

tab_ys[tmp]=label1；

tmp=tab_fz[tmp]；}

tmp=label1；

while（tab_fz[tmp]）tmp=tab_fz[tmp]；

tab_fz[tmp]=label2；}

表1 二值图像示意图表

表2 中间图像示意图表

表2 中间图像示意图

按照上述规则遍历完表2 后，辅助表和映射表中的数值如表1 所示。表2 中9 个不同标号在映射表中已分割为3 个不同的连通域，分别为5、3 和9。在辅助表中，可以发现连通域5 对应的所有标号构成了一个标号链表{4，1，2，6，7，8}。

表1 最终映射表和辅助表 

对表2进行第二次遍历，更改各个标号为映射表中对应的实际连通域，则完成了图像目标的准确分割，表3为分割结果。

表3 最终图像示意图表 

Claims (3)

1.一种遥感二值图像目标区域分割的方法，其特征在于，采用以下算法：

功能：二值图像初步标号，建立标号映射表

输入：遥感二值图像Fij

输出：中间图像Fij’，标号映射表tab_ys

（1）定义映射表tab_ys ，定义辅助表tab_fz ，总标号label = 0 ，映射表和辅助表各元素初始化为0；

（2）从上到下逐行扫描图像，若在第i 行中发现1 值段Fi,m~Fi,n ，为其申请结点，进行（3）；否则返回（2）进行下一行扫描；

（3）据公式: (L1.begin≤L2.end + 1)∩(L1.end≥L2.begin - 1)，与第i - 1 行中线段结点逐个判断连通性，若都不连通（适合i = 1），则label = label + 1 赋值给F’i,m~F’i,n和当前线段结点,tab_ys[label] = label ，返回（2）；否则有k(k ≥1) 个线段结点与之连通，线段标号为label(1)~label(k) ，用label(1) 赋值给当前线段结点和F’i,m~F’i,n ，进行（4）；

（4）若k =1，返回（2）；否则逐个检查tab_ys[label(x)](2≤x≤k) 与tab_ys[label(1)] 的关系，若不相等，则执行correct(tab_ys[label(1)],tab_ys[label(x)]) ；

全部检查完毕返回（2）；

根据初步标号结果和映射表，执行Fij' = tab_ysFij' ，即可完成图像的最终连通域标记。

2.根据权利要求1所述的一种遥感二值图像目标区域分割的方法，其特征在于，所述的辅助表负责记录各个连通域对应的全部标号，如连通域j对应的全部标号为{a，b，c} ，则tab_fz[ j] = a ，tab_fz[a] = b ，tab_fz[b] = c ，tab_fz[c] = 0 ，从而形成了一个标号链表。

3.根据权利要求1所述的一种遥感二值图像目标区域分割的方法，其特征在于，所述的correct(tab_ys[label(1)],tab_ys[label(x)])为：

void correct（int label1，int labelx）{

int tmp=labelx；

while（tmp）{

tab_ys[tmp]=label1；

tmp=tab_fz[tmp]；}

tmp=label1；

while（tab_fz[tmp]）tmp=tab_fz[tmp]；

tab_fz[tmp]=labelx；}。