[发明专利]面向输电线路基于区域宽度信息的图像骨架提取方法

摘要

本发明公布了一种面向输电线路基于区域宽度信息的图像骨架提取方法，包括步骤如下：利用基于图像灰度梯度和灰度统计直方图的自适应预分割方法将输入的彩色图像转化为灰度图像并进行自适应图像预分割得到灰度图像IG0；分别从四个方向扫描灰度图像IG0，得到图像在各个方向的扫描线段序列；标记扫描骨架图像和扫描宽度信息图；将四个方向的扫描骨架图合成得到灰度图像IG0的综合扫描骨架图S；通过对综合扫描骨架图S的每个初始骨架点进行规格化，即得到灰度图像IG0的携带区域宽度信息的骨架图。本发明算法实现简单、计算开销小、运算速度快。

主权项

一种面向输电线路基于区域宽度信息的图像骨架提取方法，包括如下步骤：(1)利用基于图像灰度梯度和灰度统计直方图的自适应预分割方法将输入的彩色图像转化为灰度图像并进行自适应图像预分割得到灰度图像IG0：(1.1)将输入的彩色图像转换为灰度图像IG；(1.2)在水平方向逐行扫描步骤(1.1)所述的灰度图像IG，得到灰度图像IG的每个像素在水平方向的灰度梯度，形成灰度图像IG的灰度梯度图Igrad；(1.3)计算灰度梯度：对步骤(1.2)所述的灰度梯度图Igrad中属于[4，255]区间的灰度梯度值进行统计平均，得到灰度图像IG的统计平均灰度梯度Gs，其中4为视觉器官可感知的最小灰度梯度，255为灰度图像IG的最大灰度梯度；将统计平均梯度Gs与视觉器官能明显感知的灰度梯度区间Ge＝[16，32]比较，得到灰度图像IG的灰度梯度Gg：如Gs＜16，则Gg＝16；如Gs＞32，则Gg＝32；否则Gg＝Gs；(1.4)采用灰度图像IG的所有灰度等级建立其灰度直方图HG，所述灰度直方图HG表示灰度图像IG的每个灰度等级k对应的像素数量nk，其中横轴为灰度级别k；竖轴为相应灰度级别k的像素个数nk，灰度等级k的取值范围为0至255；(1.5)扫描步骤(1.4)所述的灰度直方图HG，得到灰度属性列表Ak，即针对每个灰度等级k，根据其像素数量nk与其前一个灰度等级(k-1)的像素数量nk-1的大小关系，将对应的灰度属性列表Ak标记为“上升”、“下降”或“不变”：当nk＞nk-1，则标记灰度级k为“上升”；当nk＜nk-1，则标记灰度级k为“下降”；当nk＝nk-1，则标记灰度级k为“不变”；灰度级别0和255均标记为“不变”；(1.6)扫描灰度属性列表Ak，并在灰度属性列表Ak中标记出“谷底”区域序列Bi，j(i，j∈0-255，j≥i)，“谷底”区域Bi，j(i，j∈0-255，j≥i)是指灰度属性列表Ak中属性Ai为“下降”、Aj为“上升”的区间；(1.7)找到步骤(1.4)所述的灰度直方图HG的峰值np＝MAX(nk)，确定灰度直方图HG的峰值np对应的灰度等级p所对应的图像分割灰度阈值区间Tp＝[t0，t1]，其中MAX()表示求最大值；t0＝MAX(0，p-Gg)，t1＝MIN(255，p+Gg)；在灰度属性列表Ak上，从t0向峰值对应的灰度等级p方向(从小到大)扫描，当在[t0，p]区间找到第一个“谷底”，则取该“谷底”对应的灰度等级为新的t0值，找不到“谷底”则t0不变；同样地，从t1向峰值对应的灰度等级p方向(从大到小)扫描，当在[p，t1]区间找到第一个“谷底”，则取该“谷底”对应的灰度等级为新的t1值，找不到“谷底”则t1不变；(1.8)扫描整个灰度图像IG，将灰度等级属于步骤(1.7)所述的图像分割灰度阈值区间Tp＝[t0，t1]的像素的灰度等级设置为p；(1.9)将灰度直方图HG中满足k∈Tp＝[t0，t1]的nk设为0，即nk＝0；(1.10)重复步骤(1.7)、(1.8)、(1.9)，直到灰度图像直方图的峰值为0，完成图像像素的聚类，实现图像的预分割，得到预分割后的灰度图像IG0；其特征在于还包括以下步骤：(2)分别从四个方向扫描整个灰度图像IG0，得到灰度图像IG0在各个方向的扫描线段集合{L0(r1，c1)，(r1，c2)}、{L1(r1，c1)，(r2，c2)}、{L2(r1，c1)，(r2，c1)}、{L3(r1，c1)，(r2，c2)}，其中集合元素的上标表示扫描线段的起点和终点，r1、r2为行号，取值区间为[0，ht-1]，c1、c2为列号，取值区间为[0，wid-1]，ht、wid分别为所述图像的行数和列数，即图像的高度和宽度，为大于1的自然数，集合元素的数字下标表示扫描方向d＝0、1、2、3，所述四个方向如下：方向0--水平方向、方向1--右上45度方向、方向2--垂直方向、方向3--左上45度方向，下同；(3)建立与灰度图像IG0有相同尺寸大小的扫描骨架图像Sd(d＝0、1、2、3)和综合扫描骨架图S，并将所述Sd(d＝0、1、2、3)和综合扫描骨架图S中所有像素的灰度值初始化为0，分别用于记录灰度图像IG0在扫描方向d的扫描骨架信息以及综合扫描骨架信息；建立与灰度图像IG0有相同尺寸大小的扫描宽度信息图Wd(d＝0、1、2、3)，并将扫描宽度信息图Wd(d＝0、1、2、3)中所有像素的灰度值初始化为0，分别用于记录灰度图像IG0中每个像素在所在扫描方向d的扫描线段的长度信息；(4)将步骤(2)所述的每个方向的扫描线段集合中的每条扫描线段的中点在对应扫描方向的扫描骨架图像Sd(d＝0、1、2、3)的对应像素用所述线段的长度作为该像素的灰度值标记出来，将所述扫描线段上的每个像素在对应扫描方向的扫描宽度信息图Wd(d＝0、1、2、3)上用所述线段的长度作为灰度值标记出来，分别得到灰度图像IG0在四个方向的扫描骨架图Sd(d＝0、1、2、3)和扫描宽度信息图Wd(d＝0、1、2、3)；(5)将步骤(4)得到四个方向的扫描骨架图Sd(d＝0、1、2、3)合成得到灰度图像IG0的综合扫描骨架图S；(6)通过对步骤(5)得到综合扫描骨架图S的每个初始骨架点进行规格化，即得到灰度图像IG0的携带区域宽度信息的骨架图。