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


Posted in Python onJanuary 30, 2019

理想高通滤波实现 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()

效果:

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

总结

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

Python 相关文章推荐
python实现apahce网站日志分析示例
Apr 02 Python
Python中基础的socket编程实战攻略
Jun 01 Python
Python模拟脉冲星伪信号频率实例代码
Jan 03 Python
Python使用Tkinter实现机器人走迷宫
Jan 22 Python
Python实现的txt文件去重功能示例
Jul 07 Python
删除DataFrame中值全为NaN或者包含有NaN的列或行方法
Nov 06 Python
Python列表切片操作实例总结
Feb 19 Python
全面了解django的缓存机制及使用方法
Jul 22 Python
python实现输入的数据在地图上生成热力图效果
Dec 06 Python
python七种方法判断字符串是否包含子串
Aug 18 Python
Pycharm的Available Packages为空的解决方法
Sep 18 Python
python中sys模块的介绍与实例
Apr 17 Python
对DataFrame数据中的重复行,利用groupby累加合并的方法详解
Jan 30 #Python
WIn10+Anaconda环境下安装PyTorch(避坑指南)
Jan 30 #Python
对dataframe数据之间求补集的实例详解
Jan 30 #Python
如何在Django中添加没有微秒的 DateTimeField 属性详解
Jan 30 #Python
Pandas统计重复的列里面的值方法
Jan 30 #Python
自学python的建议和周期预算
Jan 30 #Python
python DataFrame 取差集实例
Jan 30 #Python
You might like
php中判断文件存在是用file_exists还是is_file的整理
2012/09/12 PHP
php 注册时输入信息验证器的实现详解
2013/07/05 PHP
php正则匹配文章中的远程图片地址并下载图片至本地
2015/09/29 PHP
超强多功能php绿色集成环境详解
2017/01/25 PHP
thinkPHP5框架auth权限控制类与用法示例
2018/06/12 PHP
各种页面定时跳转(倒计时跳转)代码总结
2013/10/24 Javascript
用jQuery与JSONP轻松解决跨域访问的问题
2014/02/04 Javascript
JavaScript sup方法入门实例(把字符串显示为上标)
2014/10/20 Javascript
js实现点击图片改变页面背景图的方法
2015/02/28 Javascript
JavaScript把数组作为堆栈使用的方法
2015/03/20 Javascript
JavaScript 对象深入学习总结(经典)
2015/09/29 Javascript
jQuery页面加载初始化的3种方法(推荐)
2016/06/02 Javascript
Vue.js表单控件实践
2016/10/27 Javascript
H5手机端多文件上传预览插件
2017/04/21 Javascript
详解AngularJs路由之Ui-router-resolve(预加载)
2017/06/13 Javascript
js移动端事件基础及常用事件库详解
2017/08/15 Javascript
Angular2.0/4.0 使用Echarts图表的示例代码
2017/12/07 Javascript
10行原生JS实现文字无缝滚动(超简单)
2018/01/02 Javascript
JS实现数组的增删改查操作示例
2018/08/29 Javascript
vue移动端监听滚动条高度的实现方法
2018/09/03 Javascript
vue template中slot-scope/scope的使用方法
2018/09/06 Javascript
微信小程序移动拖拽视图-movable-view实例详解
2019/08/17 Javascript
python使用 request 发送表单数据操作示例
2019/09/25 Python
python 使用事件对象asyncio.Event来同步协程的操作
2020/05/04 Python
python 基于opencv去除图片阴影
2021/01/26 Python
美国网上鞋子零售商:Dr. Scholl’s Shoes
2017/11/17 全球购物
同步和异步有何异同,在什么情况下分别使用他们?举例说明
2014/02/27 面试题
幼儿教师自我鉴定
2013/11/02 职场文书
如何写好建议书
2014/03/13 职场文书
毕业实习自我鉴定范文2014
2014/09/26 职场文书
优秀班主任推荐材料
2014/12/17 职场文书
医生辞职信范文
2015/03/02 职场文书
2016年幼儿园万圣节活动总结
2016/04/05 职场文书
导游词之重庆渣滓洞
2020/01/08 职场文书
python实现网络五子棋
2021/04/11 Python
世界十大动漫制作公司排行榜,迪士尼上榜,第二是美国代表性文化符
2022/03/18 欧美动漫