Justheartyoung 发表于 2023-8-6 18:42:15

报错,求帮忙修改。

代码:
import math
import numpy as np
import matplotlib.pyplot as plt


# 生成高斯核
def gaussian_create():
    sigma1 = sigma2 = 1
    gaussian_sum = 0
    g = np.zeros()
    for i in range(3):
      for j in range(3):
            g = math.exp(-1 / 2 * (np.square(i - 1) / np.square(sigma1)
                                       + (np.square(j - 1) / np.square(sigma2)))) / (
                              2 * math.pi * sigma1 * sigma2)
            gaussian_sum = gaussian_sum + g
    g = g / gaussian_sum# 归一化
    return g


# 产生灰度图
def gray_fuc(rgb):
    return np.dot(rgb[..., :3], )


# 高斯卷积
def gaussian_blur(gray_img, g):
        '''
        gray_img:灰度图
        g:高斯核
        '''
    gray_img = np.pad(gray_img, ((1, 1), (1, 1)), constant_values=0)# 填充
    h, w = gray_img.shape
    new_gray_img = np.zeros()
    for i in range(h - 2):
      for j in range(w - 2):
            new_gray_img = np.sum(gray_img * g)
    return new_gray_img


# 求高斯偏导
def partial_derivative(new_gray_img):
        '''
        new_gray_img:高斯卷积后的灰度图
        '''
    new_gray_img = np.pad(new_gray_img, ((0, 1), (0, 1)), constant_values=0)# 填充
    h, w = new_gray_img.shape
    dx_gray = np.zeros()# 用来存储x方向偏导
    dy_gray = np.zeros()# 用来存储y方向偏导
    df_gray = np.zeros()# 用来存储梯度强度
    for i in range(h - 1):
      for j in range(w - 1):
            dx_gray = new_gray_img - new_gray_img
            dy_gray = new_gray_img - new_gray_img
            df_gray = np.sqrt(np.square(dx_gray) + np.square(dy_gray))
    return dx_gray, dy_gray, df_gray


# 非极大值抑制
def non_maximum_suppression(dx_gray, dy_gray, df_gray):
        '''
        dx_gray:x方向梯度矩阵
        dy_gray:y方向梯度矩阵
        df_gray:梯度强度矩阵
        '''
    df_gray = np.pad(df_gray, ((1, 1), (1, 1)), constant_values=0)# 填充
    h, w = df_gray.shape
    for i in range(1, h - 1):
      for j in range(1, w - 1):
            if df_gray != 0:
                gx = math.fabs(dx_gray)
                gy = math.fabs(dy_gray)
                if gx > gy:
                  weight = gy / gx
                  grad1 = df_gray
                  grad2 = df_gray
                  if gx * gy > 0:
                        grad3 = df_gray
                        grad4 = df_gray
                  else:
                        grad3 = df_gray
                        grad4 = df_gray
                else:
                  weight = gx / gy
                  grad1 = df_gray
                  grad2 = df_gray
                  if gx * gy > 0:
                        grad3 = df_gray
                        grad4 = df_gray
                  else:
                        grad3 = df_gray
                        grad4 = df_gray
                t1 = weight * grad1 + (1 - weight) * grad3
                t2 = weight * grad2 + (1 - weight) * grad4
                if df_gray > t1 and df_gray > t2:
                  df_gray = df_gray
                else:
                  df_gray = 0
    return df_gray


# 双阈值过滤
def double_threshold(df_gray, low, high):
        '''
        df_gray:梯度强度矩阵
        low:低阈值
        high:高阈值
        '''
    h, w = df_gray.shape
    for i in range(1, h - 1):
      for j in range(1, w - 1):
            if df_gray < low:
                df_gray = 0
            elif df_gray > high:
                df_gray = 1
            elif (df_gray > high) or (df_gray > high) or (
                  df_gray > high) or (df_gray > high) or (df_gray > high) or (
                  df_gray > high) or (df_gray > high) or (df_gray > high):
                df_gray = 1
            else:
                df_gray = 0
    return df_gray


if __name__ == '__main__':
    # 读取图像
    img = plt.imread('F:\\duneline\\dune\\dune.tif')
    # 生成高斯核
    gaussian = gaussian_create()
    # 生成灰度图
    gray = gray_fuc(img)
    # 高斯卷积
    new_gray = gaussian_blur(gray, gaussian)
    # 求偏导
    d = partial_derivative(new_gray)
    dx = d
    dy = d
    df = d
    # 非极大值抑制
    new_df = non_maximum_suppression(dx, dy, df)
    # 双阈值过滤,并将图像转换成转化二值图
    low_threshold = 0.15 * np.max(new_df)
    high_threshold = 0.2 * np.max(new_df)
    result = double_threshold(new_df, low_threshold, high_threshold)
    # 输出图像
    plt.imshow(result, cmap="gray")
    plt.axis("off")
    plt.show()
报错:
File "F:/duneline/t.py", line 32
    gray_img = np.pad(gray_img, ((1, 1), (1, 1)), constant_values=0)# 填充
                                                                        ^
IndentationError: unindent does not match any outer indentation level

Process finished with exit code 1

Mike_python小 发表于 2023-8-6 18:42:32

这个错误是由于缩进混用引起的。在Python中,缩进是非常重要的,它用于表示代码块的开始和结束。你需要确保在每个代码块中使用一致的缩进。

根据你提供的代码,问题出现在第32行的"gray_img = np.pad(gray_img, ((1, 1), (1, 1)), constant_values=0)"这一行。它的缩进与其周围的代码不匹配。

为了解决这个问题,请检查你的代码并确保在第32行及其周围的代码块中使用相同的缩进,通常是四个空格或一个制表符。正确的代码应该像这样:


# 高斯卷积
def gaussian_blur(gray_img, g):
    '''
    gray_img:灰度图
    g:高斯核
    '''
    gray_img = np.pad(gray_img, ((1, 1), (1, 1)), constant_values=0)# 填充
    h, w = gray_img.shape
    new_gray_img = np.zeros()
    for i in range(h - 2):
      for j in range(w - 2):
            new_gray_img = np.sum(gray_img * g)
    return new_gray_img


请记住,在修改代码后,你可能需要确保整个文件中的缩进保持一致。

如果回答对你有帮助,请给我一个最佳答案!这对我非常重要!{:10_254:}{:10_254:}{:10_254:}

isdkz 发表于 2023-8-6 18:42:33

这个错误是缩进错误,发生在第32行。根据Python的语法规则,代码块中的每一行都必须使用相同的缩进级别。在这个问题中,第32行的缩进级别与其上一行不匹配。

要解决这个问题,你需要检查代码的缩进,并确保每一行都使用相同的缩进级别。在这种情况下,你可以将第32行的缩进级别修改为与其上一行相同的级别。例如,你可以将第32行的缩进级别修改为4个空格或一个制表符。

修改后的代码如下所示:

# 高斯卷积
def gaussian_blur(gray_img, g):
    '''
    gray_img:灰度图
    g:高斯核
    '''
    gray_img = np.pad(gray_img, ((1, 1), (1, 1)), constant_values=0)# 填充
    h, w = gray_img.shape
    new_gray_img = np.zeros()
    for i in range(h - 2):
      for j in range(w - 2):
            new_gray_img = np.sum(gray_img * g)
    return new_gray_img


希望能帮助到你!如果你还有其他问题,请随时提问。

学习编程中的Ben 发表于 2023-8-6 18:43:09

这段代码的错误是缩进错误。在第32行的gray_img = np.pad(gray_img, ((1, 1), (1, 1)), constant_values=0)上面有额外的空格导致了IndentationError。以下是修改后的代码:

import math
import numpy as np
import matplotlib.pyplot as plt

# 生成高斯核
def gaussian_create():
    sigma1 = sigma2 = 1
    gaussian_sum = 0
    g = np.zeros()
    for i in range(3):
      for j in range(3):
            g = math.exp(-1 / 2 * (np.square(i - 1) / np.square(sigma1)
                           + (np.square(j - 1) / np.square(sigma2)))) / (
                           2 * math.pi * sigma1 * sigma2)
            gaussian_sum = gaussian_sum + g
    g = g / gaussian_sum# 归一化
    return g

# 产生灰度图
def gray_fuc(rgb):
    return np.dot(rgb[..., :3], )

# 高斯卷积
def gaussian_blur(gray_img, g):
    gray_img = np.pad(gray_img, ((1, 1), (1, 1)), constant_values=0)# 填充
    h, w = gray_img.shape
    new_gray_img = np.zeros()
    for i in range(h - 2):
      for j in range(w - 2):
            new_gray_img = np.sum(gray_img * g)
    return new_gray_img

# 求高斯偏导
def partial_derivative(new_gray_img):
    new_gray_img = np.pad(new_gray_img, ((0, 1), (0, 1)), constant_values=0)# 填充
    h, w = new_gray_img.shape
    dx_gray = np.zeros()# 用来存储x方向偏导
    dy_gray = np.zeros()# 用来存储y方向偏导
    df_gray = np.zeros()# 用来存储梯度强度
    for i in range(h - 1):
      for j in range(w - 1):
            dx_gray = new_gray_img - new_gray_img
            dy_gray = new_gray_img - new_gray_img
            df_gray = np.sqrt(np.square(dx_gray) + np.square(dy_gray))
    return dx_gray, dy_gray, df_gray

# 非极大值抑制
def non_maximum_suppression(dx_gray, dy_gray, df_gray):
    df_gray = np.pad(df_gray, ((1, 1), (1, 1)), constant_values=0)# 填充
    h, w = df_gray.shape
    for i in range(1, h - 1):
      for j in range(1, w - 1):
            if df_gray != 0:
                gx = math.fabs(dx_gray)
                gy = math.fabs(dy_gray)
                if gx > gy:
                  weight = gy / gx
                  grad1 = df_gray
                  grad2 = df_gray
                  if gx * gy > 0:
                        grad3 = df_gray
                        grad4 = df_gray
                  else:
                        grad3 = df_gray
                        grad4 = df_gray
                else:
                  weight = gx / gy
                  grad1 = df_gray
                  grad2 = df_gray
                  if gx * gy > 0:
                        grad3 = df_gray
                        grad4 = df_gray
                  else:
                        grad3 = df_gray
                        grad4 = df_gray
                t1 = weight * grad1 + (1 - weight) * grad3
                t2 = weight * grad2 + (1 - weight) * grad4
                if df_gray > t1 and df_gray > t2:
                  df_gray = df_gray
                else:
                  df_gray = 0
    return df_gray

# 双阈值过滤
def double_threshold(df_gray, low, high):
    h, w = df_gray.shape
    for i in range(1, h - 1):
      for j in range(1, w - 1):
            if df_gray < low:
                df_gray = 0
            elif df_gray > high:
                df_gray = 1
            elif (df_gray > high) or (df_gray > high) or (
                  df_gray > high) or (df_gray > high) or (df_gray > high) or (
                  df_gray > high) or (df_gray > high) or (df_gray > high):
                df_gray = 1
            else:
               

如果回答对你有帮助,请给我一个最佳答案!{:10_254:}{:10_254:}{:10_254:}

sfqxx 发表于 2023-8-6 20:08:56

学习编程中的Ben 发表于 2023-8-6 18:43
这段代码的错误是缩进错误。在第32行的gray_img = np.pad(gray_img, ((1, 1), (1, 1)), constant_values=0) ...

没了。
页: [1]
查看完整版本: 报错,求帮忙修改。