上理论:
希望对大家有用。
Hilditch 细化算法的步骤为:
对图像从左向右从上向下迭代每个像素,是为一个迭代周期。在每个迭代周期中,对于每一个像素p,如果它同时满足6个条件,则标记它。在当前迭代周期结束时,则把所有标记的像素的值设为背景值。如果某次迭代周期中不存在标记点(即满足6个条件的像素),则算法结束。假设背景值为0,前景值为1,则:
6个条件为:
(I):p 为1,即p不是背景;
(2):x1,x3,x5,x7不全部为1(否则把p标记删除,图像空心了);
(3):x1~x8 中,至少有2个为1(若只有1个为1,则是线段的端点。若没有为1的,则为孤立点);
(4):p的8连通联结数为1;
联结数指在像素p的3*3邻域中,和p连接的图形分量的个数:
(5)假设x3已经标记删除,那么当x3为0时,p的8联通联结数为1;
(6)假设x5已经标记删除,那么当x5为0时,p的8联通联结数为1。
以上的理论选择网络博客:http://www.cnblogs.com/xiaotie/archive/2010/08/12/1797760.html (此博客上有C#版本的Hilditch细化算法,但是使用unsafe的代码不是太好理解的。
我上一个自己写的代码
/// <summary>
/// Hilditch细化算法
/// </summary>
/// <param name="input"></param>
/// <returns></returns>
private int[,] ThinnerHilditch(int[,] input)
{
int lWidth = input.GetLength(0);
int lHeight = input.GetLength(1);
bool IsModified = true;
int Counter = 1;
int[] nnb = new int[9];
//去掉边框像素
for (int i = 0; i < lWidth; i++)
{
input[i, 0] = 0;
input[i, lHeight - 1] = 0;
}
for (int j = 0; j < lHeight; j++)
{
input[0, j] = 0;
input[lWidth - 1, j] = 0;
}
do
{
Counter++;
IsModified = false;
int[,] nb = new int[3, 3];
for (int i = 1; i < lWidth; i++)
{
for (int j = 1; j < lHeight; j++)
{
//条件1必须为黑点
if (input[i, j] != 1)
{
continue;
}
//取3*3领域
for (int m = 0; m < 3; m++)
{
for (int n = 0; n < 3; n++)
{
nb[m, n] = input[i - 1 + m, j - 1 + n];
}
}
//复制
nnb[0] = nb[2, 1]==1?0:1;
nnb[1] = nb[2, 0]==1?0:1;
nnb[2] = nb[1, 0]==1?0:1;
nnb[3] = nb[0, 0]==1?0:1;
nnb[4] = nb[0, 1]==1?0:1;
nnb[5] = nb[0, 2]==1?0:1;
nnb[6] = nb[1, 2]==1?0:1;
nnb[7] = nb[2, 2]==1?0:1;
// 条件2:p0,p2,p4,p6 不皆为前景点
if (nnb[0] == 0 && nnb[2] == 0 && nnb[4] == 0 && nnb[6] == 0)
{
continue;
}
// 条件3: p0~p7至少两个是前景点
int iCount = 0;
for (int ii = 0; ii < 8; ii++)
{
iCount += nnb[ii];
}
if (iCount > 6) continue;
// 条件4:联结数等于1
if (DetectConnectivity(nnb) != 1)
{
continue;
}
// 条件5: 假设p2已标记删除,则令p2为背景,不改变p的联结数
if (input[i, j - 1] == -1)
{
nnb[2] = 1;
if (DetectConnectivity(nnb) != 1)
continue;
nnb[2] = 0;
}
// 条件6: 假设p4已标记删除,则令p4为背景,不改变p的联结数
if (input[i, j + 1] == -1)
{
nnb[6] = 1;
if (DetectConnectivity(nnb) != 1)
continue;
nnb[6] = 0;
}
input[i, j] = -1;
IsModified = true;
}
}
for (int i = 0; i < lWidth; i++)
{
for (int j = 0; j < lHeight; j++)
{
if (input[i, j] == -1)
{
input[i, j] = 0;
}
}
}
} while (IsModified);
return input;
}
希望对大家有用。