高通滤波

import cv2
import numpy as np
import matplotlib.pyplot as plt
img=cv2.imread('C:/Users/xpp/Desktop/Lena.png',0)#原始图像
f=np.fft.fft2(img)
fshift=np.fft.fftshift(f)#傅里叶变换
rows,cols=img.shape
crow,ccol=int(rows/2),int(cols/2)#计算中心位置坐标
fshift[crow-30:crow+30,ccol-30:ccol+30]=0#高通滤波
ishift=np.fft.ifftshift(fshift)#逆傅里叶变换
iimg=np.fft.ifft2(ishift)
iimg=np.abs(iimg)
plt.subplot(121)
plt.imshow(img,cmap='gray')
plt.title('original')
plt.axis('off')
plt.subplot(122)
plt.imshow(iimg,cmap='gray')
plt.title('iimg')
plt.axis('off')
plt.show()