Python:二进制2D数组中的轮廓

时间:2017-01-27 15:38:42

标签: python numpy contour

我想以最简单的方式(没有百万次检查图像的边界)从二进制2D数组中获得n个像素的宽度进入正区域的轮廓。

示例:

img = np.array([
               [0, 0, 0, 1, 1, 1, 1, 1, 0],
               [0, 1, 1, 1, 1, 1, 1, 1, 1],
               [0, 1, 1, 1, 0, 0, 0, 0, 0],
               ])

用于呼叫,例如width = 1.如果img [i,j] == 1且任何邻居(img [i + 1,j],img [i-1,j],img [i,j-1],img [],则像素为正i,j + 1])为0。

contour1 = get_countor(img, width = 1)
contour1 = ([
               [0, 0, 0, 1, 0, 0, 0, 1, 0],
               [0, 1, 1, 0, 1, 1, 1, 1, 1],
               [0, 1, 0, 1, 0, 0, 0, 0, 0],
            ])

或致电例如width = 2.宽度= 1的所有像素都是正的,以及满足img [i,j] == 1的像素,并且其中2个索引距离(欧几里德距离)存在值为0的像素。

contour2 = get_countor(img, width = 2)
contour2 = ([
               [0, 0, 0, 1, 1, 1, 1, 1, 0],
               [0, 1, 1, 1, 1, 1, 1, 1, 1],
               [0, 1, 1, 1, 0, 0, 0, 0, 0],
            ]) 

感谢您的帮助。

3 个答案:

答案 0 :(得分:1)

不是这个问题的确切答案,而是分享一种在图像中绘制轮廓的简单方法;对于那些正在寻找它的人来说。

from PIL import Image
from PIL import ImageFilter
import numpy as np


def draw_contour(img, mask, contour_width, contour_color):
    """Draw contour on a pillow image from a numpy 2D mask."""
    contour = Image.fromarray(mask)
    contour = contour.resize(img.size)
    contour = contour.filter(ImageFilter.FIND_EDGES)
    contour = np.array(contour)

    # make sure borders are not drawn
    contour[[0, -1], :] = 0
    contour[:, [0, -1]] = 0

    # use a gaussian to define the contour width
    radius = contour_width / 10
    contour = Image.fromarray(contour)
    contour = contour.filter(ImageFilter.GaussianBlur(radius=radius))
    contour = np.array(contour) > 0
    contour = np.dstack((contour, contour, contour))

    # color the contour
    ret = np.array(img) * np.invert(contour)
    if contour_color != 'black':
        color = Image.new(img.mode, img.size, contour_color)
        ret += np.array(color) * contour

    return Image.fromarray(ret)

检查此测试输出:enter image description here

我在为PR工作时写了这个解决方案。

答案 1 :(得分:0)

import numpy as np
import pandas as pd
import random

df = pd.DataFrame([], columns=[0,1,2,3,4,5,6,7,8,9])

for i in np.arange(10):
    df.loc[len(df)] = np.random.randint(0,2,10)

df = df.astype(bool)

contour = df & ((df-df.shift(-1, axis=0).fillna(1))|(df-df.shift(1,axis=0).fillna(1))|(df-df.shift(-1,axis=1).fillna(1))|(df-df.shift(1,axis=1).fillna(1)))

输出:

DF:

enter image description here

轮廓:

enter image description here

希望这会有所帮助

答案 2 :(得分:0)

我认为你在寻找的是什么 scipy.misc.imfilter(img, "find_edges")

给定一个二进制数组img,这将产生一个0255的数组,因此你需要除以255.正如我所看到的,滤波器的宽度=通过另一次应用width = 1的过滤器获得2,所以最后你的函数看起来像

def get_countor(img, width = 1):
    for i in range(width):
        img = scipy.misc.imfilter(img, "find_edges")/255
    return img