如何使用skimage获得hough线峰值的x，y坐标

Question

h, theta, d = transform.hough_line(outlines)
for acum, angle, dist in zip(*transform.hough_line_peaks(h, theta, d)):
    y0 = (dist - 0 * np.cos(angle)) / np.sin(angle)
    y1 = (dist - outlines.shape[1] * np.cos(angle)) / np.sin(angle)
    x0 = ...
    x1 = ...
    rr,cc,_ = draw.line_aa(x0,y0,x1,y1)

我想要的是我的轮廓形状范围之间的x0和x1值，即640,640（2D）。 我想将y0和y1缩放到outline.shape的大小。

Answer 1

使用该公式计算的y0, y1坐标对应于线与图像边缘相交的位置。 这就是它包含0和outlines.shape[1] 。

y0对应于线与列0相交的行，因此0 * cos(angle) 。

y1对应于线与图像的最后一列相交的行，即其宽度（即outlines.shape[1] ）

因此，您可以从(0, y0)到(width, y1)绘制一条线来强调检测到的线条。 使用，例如outlines[line(0, y0, width-1, y1] = 1注意我放宽width - 1因为索引从0开始， width超出范围。在你的公式中不是这种情况，因为它是从dist减去

本教程很好地说明了它的工作原理以及如何将已发现的行添加到图像中（在第一部分中）。 用您选择的图像替换未出现的X形图像，您将看到线条。 理想情况下，您的图像应该被二值化并且没有太多点，因此请尝试将其传递给边缘检测器（Canny）或骨架化过程。