0

我有使用轮廓算法提取特征的图像（我正在进行天体物理源提取）。这种方法产生了一个“特征图”，每个像素都用一个整数“标记”（通常每张图有约 1000 个独特的特征）。

我想将每个单独的特征显示为它自己的轮廓。

我可以做到这一点的一种方法是：

for ii in range(labelmask.max()):
    contour(labelmask,levels=[ii-0.5])

但是，这非常慢，尤其是对于大图像。有没有更好（更快）的方法？

PS 一个小测试表明skimage 的 find-contours并不快。

根据@tcaswell 的评论，我需要解释为什么contour(labels, levels=np.unique(levels)+0.5))或类似的东西不起作用：

1. Matplotlib spaces each subsequent contour "inward" by a linewidth to avoid overlapping contour lines.  This is not the behavior desired for a labelmask.
2. The lowest-level contours encompass the highest-level contours
3. As a result of the above, the highest-level contours will be surrounded by a miniature version of whatever colormap you're using and will have extra-thick contours compared to the lowest-level contours.
4

1 回答 1

1

很抱歉回答我自己的问题......不耐烦（和好运）让我变得更好。

关键是使用 matplotlib 的低级 C 例程：

I = imshow(data)
E = I.get_extent()
x,y = np.meshgrid(np.linspace(E[0],E[1],labels.shape[1]), np.linspace(E[2],E[3],labels.shape[0]))

for ii in np.unique(labels):
    if ii == 0: continue
    tracer = matplotlib._cntr.Cntr(x,y,labels*(labels==ii))
    T = tracer.trace(0.5)
    contour_xcoords,contour_ycoords = T[0].T
    # to plot them:
    plot(contour_xcoords, contour_ycoords)

请注意，这labels*(labels==ii)会将每个标签的轮廓放在稍微不同的位置；labels==ii如果您想要在相邻标签之间重叠轮廓，请将其更改为。

于 2013-03-02T19:12:28.410 回答