二值圖像的細化

細化算法的參考http://blog.sina.com.cn/s/blog_6f611c3001017y5m.html

直接上案例，對手寫數字進行細化：

activity:

//二值圖像細化

    public void remove(View v) {

        Bitmap bim=BitmapFactory.decodeResource(getResources(),R.drawable.i);

        int w = bim.getWidth();

        int h = bim.getHeight();

        int[] neighbor = new int[8];

        int[] imgBuf = new int[w * h];

        // memset(mark,0,w*h);

        bim.getPixels(imgBuf, 0, w, 0, 0, w, h);

        int[] mark = new int[w * h];

        boolean loop = true;

        int x, y, k;

        int markNum = 0;

        while (loop) {

            loop = false;

            // 第一步

            markNum = 0;

            for (y = 1; y < h - 1; y++) {

                for (x = 1; x < w - 1; x++) {

                    // 條件1：p必須是前景點

                    if (imgBuf[y * w + x] == Color.WHITE)

                        continue;

                    neighbor[0] = (imgBuf[y * w + x + 1] >> 8 & 0xff);

                    neighbor[1] = (imgBuf[(y - 1) * w + x + 1] >> 16 & 0xff);

                    neighbor[2] = (imgBuf[(y - 1) * w + x] >> 16 & 0xff);

                    neighbor[3] = (imgBuf[(y - 1) * w + x - 1] >> 16 & 0xff);

                    neighbor[4] = (imgBuf[y * w + x - 1] >> 16 & 0xff);

                    neighbor[5] = (imgBuf[(y + 1) * w + x - 1] >> 16 & 0xff);

                    neighbor[6] = (imgBuf[(y + 1) * w + x] >> 16 & 0xff);

                    neighbor[7] = (imgBuf[(y + 1) * w + x + 1] >> 16 & 0xff);

                    // 條件2：2<=N(p）<=6

                    int np = (neighbor[0] + neighbor[1] + neighbor[2]

                            + neighbor[3] + neighbor[4] + neighbor[5]

                            + neighbor[6] + neighbor[7]) / 255;

                    if ((8 - np) < 2 || (8 - np) > 6)

                        continue;

                    // 條件3：S(p）=1

                    int sp = 0;

                    for (int i = 1; i < 8; i++) {

                        if (neighbor[i] - neighbor[i - 1] == 255)

                            sp++;

                    }

                    if (neighbor[0] - neighbor[7] == 255)

                        sp++;

                    if (sp != 1)

                        continue;

                    // 條件4：p2*p0*p6=0

                    if ((neighbor[2] | neighbor[0] | neighbor[6]) == 0)

                        continue;

                    // 條件5：p0*p6*p4=0

                    if ((neighbor[0] | neighbor[6] | neighbor[4]) == 0)

                        continue;

                    // 标記删除

                    mark[w * y + x] = 1;

                    markNum++;

                    loop = true;

                }

            }

            // 将标記删除的點置為背景色

            if (markNum > 0) {

                for (y = 0; y < h; y++) {

                    for (x = 0; x < w; x++) {

                        k = y * w + x;

                        if (mark[k] == 1) {

                            imgBuf[k] = Color.WHITE;

                        }

                    }

                }

            }

            // 第二步

            markNum = 0;

            for (y = 1; y < h - 1; y++) {

                for (x = 1; x < w - 1; x++) {

                    // 條件1：p必須是前景點

                    if (imgBuf[y * w + x] == Color.WHITE)

                        continue;

                    neighbor[0] = (imgBuf[y * w + x + 1] >> 16 & 0xff);

                    neighbor[1] = (imgBuf[(y - 1) * w + x + 1] >> 16 & 0xff);

                    neighbor[2] = (imgBuf[(y - 1) * w + x] >> 16 & 0xff);

                    neighbor[3] = (imgBuf[(y - 1) * w + x - 1] >> 16 & 0xff);

                    neighbor[4] = (imgBuf[y * w + x - 1] >> 16 & 0xff);

                    neighbor[5] = (imgBuf[(y + 1) * w + x - 1] >> 16 & 0xff);

                    neighbor[6] = (imgBuf[(y + 1) * w + x] >> 16 & 0xff);

                    neighbor[7] = (imgBuf[(y + 1) * w + x + 1] >> 16 & 0xff);

                    // 條件2：<=N(p)<=6

                    int np = (neighbor[0] + neighbor[1] + neighbor[2]

                            + neighbor[3] + neighbor[4] + neighbor[5]

                            + neighbor[6] + neighbor[7]) / 255;

                    if ((8 - np) < 2 || (8 - np) > 6)

                        continue;

                    // 條件3：S(p)=1

                    int sp = 0;

                    for (int i = 1; i < 8; i++) {

                        if (neighbor[i] - neighbor[i - 1] == 255)

                            sp++;

                    }

                    if (neighbor[0] - neighbor[7] == 255)

                        sp++;

                    if (sp != 1)

                        continue;

                    // 條件4：p2*p0*p4==0

                    if ((neighbor[2] | neighbor[0] | neighbor[4]) == 0)

                        continue;

                    // 條件5：p2*p6*p4==0

                    if ((neighbor[2] | neighbor[6] | neighbor[4]) == 0)

                        continue;

                    // 标記删除

                    mark[w * y + x] = 1;

                    markNum++;

                    loop = true;

                }

            }

            // 将标記删除的點置為背景色

            for (y = 0; y < h; y++) {

                for (x = 0; x < w; x++) {

                    k = y * w + x;

                    if (mark[k] == 1) {

                        imgBuf[k] = Color.WHITE;

                    }

                }

            }

        }

        Bitmap bitmap = Bitmap.createBitmap(w, h, Config.RGB_565);

        bitmap.setPixels(imgBuf, 0, w, 0, 0, w, h);

        destinationImage.setImageBitmap(bitmap);

    }

最後得到的圖像效果

 Zhang并行快速細化算法的實作，不僅要達到細化的目的，并且保證圖像輪廓的連通性。