C#圖像細化：Hilditch細化算法

http://blog.csdn.net/aofengdaxia/article/details/7439992

上理論：

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;
        }