python高通滤波_理想高通滤波实现Python opencv示例

理想高通滤波实现 python opencv

import numpy as np

import cv2

from matplotlib import pyplot as plt

plt.rcparams[‘font.sans-serif']=[‘simhei']

plt.rcparams[‘axes.unicode_minus'] = false

i = cv2.imread(‘capture3.png')

cv2.imshow(‘original',i)

(r,g,b) = cv2.split(i)

i = cv2.merge([b,g,r])

j = np.double(cv2.cvtcolor(i,cv2.color_rgb2gray))

d1 = 30

d2 = 60

d3 = 160

fuv = np.fft.fftshift(np.fft.fft2(j))

print(‘fuv',fuv)

print(i.shape)

m,n = i.shape[0],i.shape[1]

xo = np.floor(m/2)

yo = np.floor(n/2)

h1 = np.zeros((m,n))

h2 = np.zeros((m,n))

h3 = np.zeros((m,n))

for i in range(m):

for j in range(n):

d = np.sqrt((i-xo)**2+(j-yo)**2)

if d>=d1:

h1[i,j]=1

else:

h1[i,j]=0

if d>=d2:

h2[i,j]=1

else:

h2[i,j]=0

if d>=d3:

h3[i,j]=1

else:

h3[i,j]=0

guv1 = h1*fuv

guv2 = h2*fuv

guv3 = h3*fuv

g1 = np.fft.ifftshift(guv1)

g1 = np.uint8(np.real(np.fft.ifft2(g1)))

print(‘g1',g1)

g2 = np.fft.ifftshift(guv2)

g2 = np.uint8(np.real(np.fft.ifft2(g2)))

print(‘g2',g2)

g3 = np.fft.ifftshift(guv3)

g3 = np.uint8(np.real(np.fft.ifft2(g3)))

print(‘g3',g3)

plt.subplot(2,2,1),plt.imshow(i),plt.title(‘原图像')

plt.subplot(2,2,2),plt.imshow(g1),plt.title(‘d0=30')

plt.subplot(2,2,3),plt.imshow(g2),plt.title(‘d0=60')

plt.subplot(2,2,4),plt.imshow(g3),plt.title(‘d0=160')

plt.show()

效果：

总结

以上就是这篇文章的全部内容了，希望本文的内容对大家的学习或者工作具有一定的参考学习价值，谢谢大家对萬仟网的支持。如果你想了解更多相关内容请查看下面相关链接

如您对本文有疑问或者有任何想说的，请点击进行留言回复，万千网友为您解惑！