Matlab代码运行

xiaozhubanna

求助各位大佬们，我想Matlab实现对图像进行快速傅里叶变换，振幅谱乘以以0°（垂直）或90°（水平）为中心的高斯过滤器（标准差为14°），生成只剩水平信息的图片和只剩垂直信息的图片，目前用chatGPT给编写了个代码 但我那显示运转不了，有没有人能帮看下

``` python
import cv2
import numpy as np

# 读取图像
img = cv2.imread('image.jpg', 0)

# 进行快速傅里叶变换
f = np.fft.fft2(img)
fshift = np.fft.fftshift(f)

# 将频谱分成水平和垂直方向
rows, cols = 网页链接 crow, ccol = int(rows / 2), int(cols / 2)
fshift[crow-7:crow+8, :] = 0  # 垂直方向滤波器
#fshift[:, ccol-7:ccol+8] = 0  # 水平方向滤波器

# 将频谱转换回到原始状态
f_ishift = np.fft.ifftshift(fshift)
img_back = np.fft.ifft2(f_ishift)
img_back = np.abs(img_back)

# 显示原始图像和处理后的图像
cv2.imshow('Original', img)
cv2.imshow('Filtered', img_back)
cv2.waitKey(0)
cv2.destroyAllWindows()
```

代码中，我们使用了OpenCV库来读取和显示图像。在进行FFT变换后，我们将频谱沿垂直方向的中心部分设置为0，表示我们只保留垂直方向的频率分量。如果想保留水平方向的频率分量，则可以将fshift[:, ccol-7:ccol+8]设置为0。最后，我们使用ifft2函数将频谱转换回原始图像，并使用abs函数计算像素值的绝对值来获得处理后的图像。