基于先验知识修复图像裂缝和毛刺

void Label(byte [,] f, int x,int y,byte L)//递归标记算法
{
    f[x,y] = L;
    if (f[x+1,y]==255) Label(f,x+1,y,L);
    if (f[x-1,y]==255) Label(f,x-1,y,L);
    if (f[x,y+1]==255) Label(f,x,y+1,L);
    if (f[x,y-1]==255) Label(f,x,y-1,L);
    if (f[x+1,y+1]==255) Label(f,x+1,y+1,L);
    if (f[x-1,y+1]==255) Label(f,x-1,y+1,L);
    if (f[x+1,y-1]==255) Label(f,x+1,y-1,L);
    if (f[x-1,y-1]==255) Label(f,x-1,y-1,L);
}

void main()
{
    byte [,]f = LoadImg();
    if (f==null) return;

    ShowImg("f",f);

    int w = f.GetLength(0);
    int h = f.GetLength(1);

    //label the image
    byte L = 1;

    for (int y=0;y

         for (int x=0;x

     {
            if (f[x,y]==255)
            {
                Label(f,x,y,L);
                L+=1;
            }
        }
//    ShowImg("f",f);

    //Do erosion or dilation
    byte [,] g = new byte [w,h];
    for (int y=1;y

        for (int x=1;x

    {
            //腐蚀算法   
            if ((f[x,y]>=2&&f[x,y]<=4) && (f[x+1,y]>=2&&f[x+1,y]<=4) && (f[x-1,y]>=2&&f[x-1,y]<=4) && (f[x,y+1]

>=2&&f[x,y+1]<=4) && (f[x,y-1]>=2&&f[x,y-1]<=4))
                g[x,y] = 255;

             if(f[x,y]==1||f[x,y]==5) g[x,y] = 255;

            //膨胀算法
            if (f[x,y]==6 || f[x+1,y]==6 || f[x-1,y]==6
                || f[x,y+1]==6 || f[x,y-1]==6) g[x,y] = 255;               
        }
    ShowImg("g",g);

}

//基于先验知识修复图像裂缝和毛刺。
//先对图像中各个区域进行标记，由原图像可知，标记算法执行后，标记为2、3、4的区域存在毛刺，可使用腐蚀算法去除；标记为6的区域存在缝隙，可用膨胀算法去除。标记为1、5的区域正常，不施加额外算法.

//注意，循环中y是从1到h-2，x是从1到w-2，否则处理图像边界点时会出现执行错误！

//另外，对需修复区域只施加一次算法（腐蚀或膨胀），所以边界轮廓大小会发生变化。若想使边界轮廓大小不发生变化，可以同时施加腐

蚀和膨胀算法。

原图像：

修复后的图像：