如何使用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）或骨架化過程。