请大佬帮弟弟看看，问题出在哪里 怎么修改 谢谢

Justheartyoung

代码想实现对影像中沙丘脊线的提取，但是最终绘制的沙脊线感觉不对，有大佬知道问题出在哪里吗？可否给出改正代码！谢谢
from tifffile import imread
from skimage import filters, feature, img_as_ubyte, color
import matplotlib.pyplot as plt
import numpy as np
from skimage.transform import hough_line, hough_line_peaks

# 使用tifffile库读取TIFF格式的RGB影像
image = imread('F:\\duneline\\dune\\dune.tif')

# 将影像转换为灰度图像
gray_image = color.rgb2gray(image)

# 将灰度图像转换为8位无符号整数格式，因为后续的图像处理函数需要这种格式
edges = filters.sobel(img_as_ubyte(gray_image))

# 使用Hough变换找出可能的直线
tested_angles = np.linspace(-np.pi / 2, np.pi / 2, 360)
h, theta, d = hough_line(edges, theta=tested_angles)

# 在原图上画出沙丘脊线
fig, ax = plt.subplots()
ax.imshow(image, cmap=plt.cm.gray)

for _, angle, dist in zip(*hough_line_peaks(h, theta, d)):
    y0 = (dist - 0 * np.cos(angle)) / np.sin(angle)
    y1 = (dist - image.shape[1] * np.cos(angle)) / np.sin(angle)
    ax.plot((0, image.shape[1]), (y0, y1), '-r')

plt.show()

最终结果

edgs结果

tyh小号

根据您提供的代码，我看到您的问题可能出在将灰度图像转换为8位无符号整数格式时的处理上。

在使用 img_as_ubyte() 将灰度图像转换为8位无符号整数格式之前，您应该先对图像进行归一化操作，将像素值范围从0-1缩放到0-255。这样可以确保图像的像素值在正确的范围内，以便后续的图像处理函数能够正确处理。

下面是修改后的代码示例：

1. from tifffile import imread
   
2. from skimage import filters, feature, img_as_ubyte, color
   
3. import matplotlib.pyplot as plt
   
4. import numpy as np
   
5. from skimage.transform import hough_line, hough_line_peaks
   
6. 
   
7. # 使用tifffile库读取TIFF格式的RGB影像
   
8. image = imread('F:\\duneline\\dune\\dune.tif')
   
9. # 将影像转换为灰度图像
   
10. gray_image = color.rgb2gray(image)
    
11. 
    
12. # 将灰度图像归一化并转换为8位无符号整数格式
    
13. normalized_image = gray_image * 255
    
14. image_ubyte = normalized_image.astype(np.uint8)
    
15. 
    
16. # 使用Sobel边缘检测
    
17. edges = filters.sobel(image_ubyte)
    
18. # 使用Hough变换找出可能的直线
    
19. tested_angles = np.linspace(-np.pi / 2, np.pi / 2, 360)
    
20. h, theta, d = hough_line(edges, theta=tested_angles)
    
21. 
    
22. # 在原图上画出沙丘脊线
    
23. fig, ax = plt.subplots()
    
24. ax.imshow(image, cmap=plt.cm.gray)
    
25. for _, angle, dist in zip(*hough_line_peaks(h, theta, d)):
    
26.     y0 = (dist - 0 * np.cos(angle)) / np.sin(angle)
    
27.     y1 = (dist - image.shape[1] * np.cos(angle)) / np.sin(angle)
    
28.     ax.plot((0, image.shape[1]), (y0, y1), '-r')
    
29. 
    
30. plt.show()
    

复制代码

通过对灰度图像进行归一化，然后将像素值范围缩放到0-255，可以确保图像的像素值符合8位无符号整数的取值范围。这样修改后的代码应该能够正确地绘制沙丘脊线。希望对您有所帮助！如有其他问题，请随时提问。
如果问题已经解决，请设置最佳答案

学习编程中的Ben

根据您提供的代码，有几个问题可能导致最终绘制的沙脊线不正确：

1. 彩色图像转换为灰度图像时，使用了color.rgb2gray()函数，但是在导入语句中没有引入color模块。请将导入语句修改为from skimage import filters, img_as_ubyte, color。

2. 在边缘检测时，使用了filters.sobel()函数对灰度图像进行处理。这一步骤应该是正确的。

3. 在Hough变换的相关函数调用中，使用了zip(*hough_line_peaks(h, theta, d))。请将这部分修改为zip(h, theta, d)。

4. 在绘制沙脊线时，计算直线的起点和终点坐标的公式可能存在问题。您可以尝试使用以下公式计算：

1. y0 = (dist - 0 * np.cos(theta)) / np.sin(theta)
   
2. y1 = (dist - gray_image.shape[1] * np.cos(theta)) / np.sin(theta)
   

复制代码

请注意，我对代码进行了上述修改。以下是修改后的代码：

1. from tifffile import imread
   
2. from skimage import filters, img_as_ubyte, color
   
3. import matplotlib.pyplot as plt
   
4. import numpy as np
   
5. from skimage.transform import hough_line, hough_line_peaks
   
6. 
   
7. # 使用tifffile库读取TIFF格式的RGB影像
   
8. image = imread('F:\\duneline\\dune\\dune.tif')
   
9. 
   
10. # 将影像转换为灰度图像
    
11. gray_image = color.rgb2gray(image)
    
12. 
    
13. # 将灰度图像转换为8位无符号整数格式，因为后续的图像处理函数需要这种格式
    
14. edges = filters.sobel(img_as_ubyte(gray_image))
    
15. 
    
16. # 使用Hough变换找出可能的直线
    
17. tested_angles = np.linspace(-np.pi / 2, np.pi / 2, 360)
    
18. h, theta, d = hough_line(edges, theta=tested_angles)
    
19. 
    
20. # 在原图上画出沙丘脊线
    
21. fig, ax = plt.subplots()
    
22. ax.imshow(image, cmap=plt.cm.gray)
    
23. 
    
24. for _, angle, dist in zip(h, theta, d):
    
25.     y0 = (dist - 0 * np.cos(angle)) / np.sin(angle)
    
26.     y1 = (dist - gray_image.shape[1] * np.cos(angle)) / np.sin(angle)
    
27.     ax.plot((0, gray_image.shape[1]), (y0, y1), '-r')
    
28. 
    
29. plt.show()
    

复制代码

希望以上修改能够解决您的问题！如果您有任何其他疑问，请随时提问。

如果回答对你有帮助，请给我一个最佳答案！

isdkz

问题可能出在绘制沙丘脊线的部分。根据代码，绘制沙丘脊线的逻辑是使用Hough变换找出可能的直线，然后根据直线的角度和距离在原图上绘制直线。

在绘制直线的代码中，使用了 zip(*hough_line_peaks(h, theta, d)) 来遍历直线的角度和距离。然后根据直线的角度和距离计算直线上的两个点 (0, y0) 和 (image.shape[1], y1) ，并使用 ax.plot 函数绘制直线。

可能的问题在于直线上的两个点的计算方式不正确。根据代码，计算 y0 和 y1 的公式为：

1. y0 = (dist - 0 * np.cos(angle)) / np.sin(angle)
   
2. y1 = (dist - image.shape[1] * np.cos(angle)) / np.sin(angle)
   

复制代码

这里使用了 dist 、 angle 和 image.shape[1] 来计算直线上的点的纵坐标。其中 dist 是直线的距离， angle 是直线的角度， image.shape[1] 是图像的宽度。

根据代码逻辑， dist 表示直线与原点的距离， angle 表示直线的角度， image.shape[1] 表示图像的宽度。而直线上的点的纵坐标应该是根据直线方程计算得出的。

修正代码如下：

1. for _, angle, dist in zip(*hough_line_peaks(h, theta, d)):
   
2.     cos_angle = np.cos(angle)
   
3.     sin_angle = np.sin(angle)
   
4.     x0 = 0
   
5.     x1 = image.shape[1]
   
6.     y0 = (dist - x0 * cos_angle) / sin_angle
   
7.     y1 = (dist - x1 * cos_angle) / sin_angle
   
8.     ax.plot((x0, x1), (y0, y1), '-r')
   

复制代码

修正后的代码中，首先计算了 cos_angle 和 sin_angle ，然后根据直线方程 y = (dist - x * cos_angle) / sin_angle 计算直线上的点的纵坐标 y0 和 y1 ，其中 x0 为0， x1 为图像的宽度。最后使用 ax.plot 函数绘制直线。

请尝试使用修正后的代码，看看是否能够正确绘制沙丘脊线。

Justheartyoung

学习编程中的Ben 发表于 2023-8-4 16:22
根据您提供的代码，有几个问题可能导致最终绘制的沙脊线不正确：

1. 彩色图像转换为灰度图像时，使用了col ...

还是不行

Justheartyoung

tyh小号 发表于 2023-8-4 16:22
根据您提供的代码，我看到您的问题可能出在将灰度图像转换为8位无符号整数格式时的处理上。

在使用 img_ ...

运行结果和我原来还是一样，没有解决

学习编程中的Ben

Justheartyoung 发表于 2023-8-4 16:28
还是不行

如果您仍然遇到问题，可能是由于沙丘脊线提取的算法不适用于您的图像，或者图像中的沙丘脊线不够清晰。在这种情况下，您可以尝试以下方法来改进您的沙丘脊线提取结果：

调整边缘检测算法：尝试使用其他边缘检测算法，例如Canny边缘检测，看是否能够更好地检测到沙丘脊线。

参数调整：根据图像的特点和需要，尝试调整Hough变换的参数，包括角度的步长、距离的步长等，以及设置适当的阈值。

图像预处理：在进行边缘检测之前，可以尝试对图像进行预处理，例如去噪或平滑操作，以提高边缘检测的效果。可以使用滤波器如高斯滤波器或中值滤波器进行图像平滑处理。

其他算法：除了Hough变换，还有其他算法可用于提取直线特征，如RANSAC算法等。您可以尝试使用其他算法来提取沙丘脊线。

请注意，每个图像都有其独特的特点和挑战，因此并不存在一种通用的解决方案。根据您的具体情况，可能需要进行多次尝试和调整，以获得最佳的沙丘脊线提取结果。

如果您能提供更多关于图像特点和预期结果的信息，我可能能够为您提供更具体的建议。

Justheartyoung

Mike_python小 发表于 2023-8-4 16:22
根据您提供的代码，有几个问题可能导致绘制的沙脊线不正确：

1. 影像加载：首先，请确保实际路径'F:\\dune ...

运行了您的代码，和我结果一样，还是没解决

Justheartyoung

isdkz 发表于 2023-8-4 16:22
问题可能出在绘制沙丘脊线的部分。根据代码，绘制沙丘脊线的逻辑是使用Hough变换找出可能的直线，然后根据 ...

没有解决，还是那样

学习编程中的Ben

如果您仍然遇到问题，可能是由于沙丘脊线提取的算法不适用于您的图像，或者图像中的沙丘脊线不够清晰。在这种情况下，您可以尝试以下方法来改进您的沙丘脊线提取结果：

调整边缘检测算法：尝试使用其他边缘检测算法，例如Canny边缘检测，看是否能够更好地检测到沙丘脊线。

参数调整：根据图像的特点和需要，尝试调整Hough变换的参数，包括角度的步长、距离的步长等，以及设置适当的阈值。

图像预处理：在进行边缘检测之前，可以尝试对图像进行预处理，例如去噪或平滑操作，以提高边缘检测的效果。可以使用滤波器如高斯滤波器或中值滤波器进行图像平滑处理。

其他算法：除了Hough变换，还有其他算法可用于提取直线特征，如RANSAC算法等。您可以尝试使用其他算法来提取沙丘脊线。

请注意，每个图像都有其独特的特点和挑战，因此并不存在一种通用的解决方案。根据您的具体情况，可能需要进行多次尝试和调整，以获得最佳的沙丘脊线提取结果。

如果您能提供更多关于图像特点和预期结果的信息，我可能能够为您提供更具体的建议。


给个最佳答案呗!求你了!!!

学习编程中的Ben

请问还有哪里有问题QAQ

Mike_python小

Justheartyoung 发表于 2023-8-4 13:31
运行了您的代码，和我结果一样，还是没解决

1. from tifffile import imread
   
2. from skimage import filters, feature, img_as_ubyte, color
   
3. import matplotlib.pyplot as plt
   
4. import numpy as np
   
5. from skimage.transform import hough_line, hough_line_peaks
   
6. 
   
7. # 使用tifffile库读取TIFF格式的RGB影像
   
8. image = imread('F:\\duneline\\dune\\dune.tif')
   
9. 
   
10. # 将影像转换为灰度图像
    
11. gray_image = color.rgb2gray(image)
    
12. 
    
13. # 将灰度图像转换为8位无符号整数格式，因为后续的图像处理函数需要这种格式
    
14. edges = filters.sobel(img_as_ubyte(gray_image))
    
15. 
    
16. # 使用Hough变换找出可能的直线
    
17. tested_angles = np.linspace(-np.pi / 2, np.pi / 2, 360)
    
18. h, theta, d = hough_line(edges, theta=tested_angles)
    
19. 
    
20. # 设置阈值来选取直线
    
21. angle_threshold = 30  # 沙脊线的倾斜角度范围
    
22. dist_threshold = 50  # 直线到原点的距离阈值
    
23. 
    
24. # 在原图上画出沙脊线
    
25. fig, ax = plt.subplots()
    
26. ax.imshow(image, cmap=plt.cm.gray)
    
27. 
    
28. for _, angle, dist in zip(*hough_line_peaks(h, theta, d, threshold=0, min_distance=5)):
    
29.     if np.abs(angle) < angle_threshold and np.abs(dist) > dist_threshold:
    
30.         y0 = (dist - 0 * np.cos(angle)) / np.sin(angle)
    
31.         y1 = (dist - image.shape[1] * np.cos(angle)) / np.sin(angle)
    
32.         ax.plot((0, image.shape[1]), (y0, y1), '-r')
    
33. 
    
34. plt.show()
    

复制代码

这个代码试试？