自选一幅灰度图像,对其模拟大气湍流造成的模糊;然后用几何均值滤波器对其复原。

import numpy as np
from numpy.fft import fft2, fftshift, ifft2, ifftshift
import matplotlib.pyplot as plt
from skimage import io, util
import skimage.filters.rank as sfr
from scipy.signal import correlate2d

f = io.imread('aerial_view.tif')
F = fftshift(fft2(f))
[r,c] = F.shape;
u = np.arange(r).reshape((-1,1)) - np.ceil(r/2)
v = np.arange(c)-np.ceil(c/2)
k = 0.0025
H = np.exp(-k*(u*u+v*v)**(5/6)); 
G = F*H
g = np.uint8(np.real(ifft2(ifftshift(G))));
plt.figure(figsize=(12,6))
plt.rcParams['font.sans-serif'] = ['SimHei']
plt.subplot(121), io.imshow(f), plt.axis('off'), plt.title('鸟瞰图')
plt.subplot(122), io.imshow(g), plt.axis('off'), plt.title('严重湍流的影响')
plt.show()
# 几何平均滤波
g_log = np.log(np.double(g)+1);
g2 = np.exp(correlate2d(g_log, np.ones((3,3)),mode='same'))**(1/9)-1;
g2 = np.uint8(g2);
plt.figure(figsize=(12,12
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值