位置: IT常识 - 正文

OpenCV中的图像处理 —— 傅里叶变换+模板匹配(opencv如何显示图片)

编辑:rootadmin
OpenCV中的图像处理 —— 傅里叶变换+模板匹配 OpenCV中的图像处理 —— 傅里叶变换+模板匹配

推荐整理分享OpenCV中的图像处理 —— 傅里叶变换+模板匹配(opencv如何显示图片),希望有所帮助,仅作参考,欢迎阅读内容。

文章相关热门搜索词:opencv图形界面,opencv如何显示图片,opencv 图像,opencv的图像处理函数,opencv 图像,opencv 图像,opencv的图像处理函数,opencv的图像处理函数,内容如对您有帮助,希望把文章链接给更多的朋友!

现在也在逐渐深入啦,希望跟大家一起进步越来越强

目录OpenCV中的图像处理 —— 傅里叶变换+模板匹配1. 傅里叶变换1.1 Numpy实现傅里叶变换1.2 OpenCV实现傅里叶变换1.3 DFT的性能优化2. 模板匹配2.1 单对象的模板匹配2.2 多对象的模板匹配1. 傅里叶变换

关于傅里叶变换最重要的两个概念:时域与频域。以时间作为参照来观察动态世界的方法我们称其为时域分析,而频域是什么呢,它是描述信号在频率方面特性时用到的一种坐标系,频域图显示了在一个频率范围内每个给定频带内的信号量。贯穿时域与频域的方法之一就是大名鼎鼎的傅里叶分析,它可以分为傅里叶级数和傅里叶变换,傅里叶变换也就是我们这一部分要说的东西

傅里叶变换是分析线性系统的一个有力工具。 它告诉我们任何周期函数,都可以看作是不同振幅,不同相位正弦波的叠加。从数学意义上说,傅里叶变换将一个任意的周期函数分解成为无穷个正弦函数的和的形式;从物理效果上看,傅里叶变换实现了将信号从空间域到频率域的转换

在计算机视觉中傅立叶变换用于分析各种滤波器的频率特性,对于图像,使用2D离散傅里叶变换(DFT)查找频域(还有一种称为快速傅立叶变换(FFT)的快速算法)这一段文字是不是不太好理解,因为里面涉及太多比较深奥的东西了,傅里叶变换本身是比较难的一个点,在这里我就不细说了,我们只说说在计算机视觉领域我们是怎么用它的,想要深入了解的同学来看看这篇文章:深入浅出的讲解傅里叶变换(真正的通俗易懂)

对于正弦信号,如果幅度在短时间内变化比较快,则可以说它是高频信号,如果变化缓慢,则为低频信号,我们可以将相同的想法扩展到图像,图像中的振幅在哪里急剧变化?当然是在边缘点或噪声,因此,可以说边缘和噪声是图像中的高频内容

1.1 Numpy实现傅里叶变换

Numpy提供了FFT软件包来查找傅里叶变换,np.fft.fft2()为我们提供了频率转换,它将是一个复杂的数组,它的第一个参数是输入图像(灰度图像),第二个参数是可选的,它决定输出数组的大小。如果它大于输入图像的大小,则在计算FFT之前用零填充输入图像。如果小于输入图像,将裁切输入图像。如果未传递任何参数,则输出数组的大小将与输入的大小相同,但是现在获得的结果它的零频率分量(DC分量)将位于左上角,为了便于分析我们要把它居中,居中处理关系到np.fft.fftshift()函数

import cv2 as cvimport numpy as npfrom matplotlib import pyplot as pltimg = cv.imread(r'E:\image\test16.png', 0)f = np.fft.fft2(img)fshift = np.fft.fftshift(f)magnitude_spectrum = 20 * np.log(np.abs(fshift))plt.subplot(121), plt.imshow(img, cmap='gray')plt.title('Input Image'), plt.xticks([]), plt.yticks([])plt.subplot(122), plt.imshow(magnitude_spectrum, cmap='gray')plt.title('Magnitude Spectrum'), plt.xticks([]), plt.yticks([])plt.show()

我们可以看到幅度谱的中心有更多白色区域,说明图像低频内容更多。找到了幅度谱那我们是不是可以在频域中进行一些操作呢?例如高通滤波和重建图像,实质就是找到逆DFT,我们首先要用尺寸为60*60的矩形窗口遮罩抵消低频信号,然后使用np.fft.ifftshift()应用反向移位,以使DC分量再次出现在左上角。然后使用np.ifft2()函数找到逆FFT,结果同样是一个复数

import cv2 as cvimport numpy as npfrom matplotlib import pyplot as pltimg = cv.imread(r'E:\image\test15.png', 0)f = np.fft.fft2(img)fshift = np.fft.fftshift(f)magnitude_spectrum = 20 * np.log(np.abs(fshift))rows, cols = img.shapecrow, ccol = rows//2, cols//2fshift[crow - 30:crow + 31, ccol - 30:ccol + 31] = 0f_ishift = np.fft.ifftshift(fshift)img_back = np.fft.ifft2(f_ishift)img_back = np.abs(img_back)plt.subplot(131), plt.imshow(img, cmap='gray'),plt.title('Input Image'), plt.xticks([]), plt.yticks([])plt.subplot(132), plt.imshow(magnitude_spectrum, cmap='gray')plt.title('Magnitude Spectrum'), plt.xticks([]), plt.yticks([])plt.subplot(133), plt.imshow(img_back)plt.title('Result in JET'), plt.xticks([]), plt.yticks([])plt.show()OpenCV中的图像处理 —— 傅里叶变换+模板匹配(opencv如何显示图片)

1.2 OpenCV实现傅里叶变换

OpenCV为此提供了cv.dft()和cv.idft()函数。它返回与前一个相同的结果,但是有两个通道。第一个通道是结果的实部,第二个通道是结果的虚部。输入图像首先应转换为np.float32

import cv2 as cvimport numpy as npfrom matplotlib import pyplot as pltimg = cv.imread(r'E:\image\test17.png', 0)dft = cv.dft(np.float32(img), flags=cv.DFT_COMPLEX_OUTPUT)dft_shift = np.fft.fftshift(dft)magnitude_spectrum = 20 * np.log(cv.magnitude(dft_shift[:, :, 0], dft_shift[:, :, 1]))plt.subplot(121), plt.imshow(img, cmap='gray')plt.title('Input Image'), plt.xticks([]), plt.yticks([])plt.subplot(122), plt.imshow(magnitude_spectrum, cmap='gray')plt.title('Magnitude Spectrum'), plt.xticks([]), plt.yticks([])plt.show()

这一块儿代码有几个难懂的地方,没关系我们来分析一下:第一次看这段代码会有几个疑问,cv.dft()函数的参数怎么传递?cv.magnitude()函数是怎么用的?

cv.dft()函数的作用是对一维或者二维浮点数数组进行正向或反向离散傅里叶变换,其中包括4个参数,第一个即源图像,第二个参数是OutputArray类型的dst,函数调用后返回的运算结果存在这里,它的尺寸和类型取决于第三个参数flags转换标识符,它的默认值为0(参考自:opencv:dft()函数详解)

cv.magnitude()函数用来计算二维矢量的幅值,其中包括3个参数,第一个是InputArray类型的x,表示矢量的浮点型X坐标值,也就是实部,第二个参数是InputArray类型的y,表示矢量的浮点型Y坐标值,也就是虚部,第三个参数是输出的幅值

接下来我们需要做OpenCV中DFT的逆变换,上一节用了高通滤波器HPF,这一部分我们会将低通滤波器LPF应用到图像中

注意:通常,OpenCV函数cv.dft()和cv.idft()比Numpy函数更快,大约快3倍,但是Numpy函数更容易使用

我们把这一部分的代码放在后面,与DFT的性能优化放在一起更容易理解

1.3 DFT的性能优化

对于某些数组尺寸,DFT的计算性能较好,例如当数组大小为2的幂时,速度最快,对于大小为2、3和5的乘积的数组,也可以非常有效地进行处理,关于代码的性能问题,我们可以在找到DFT之前将数组的大小修改为任何最佳大小(通过填充零),对于OpenCV,我们必须手动填充零,但是对于Numpy,指定FFT计算的新大小,它将自动为您填充零

关于寻找最优大小,OpenCV为此提供了一个函数:cv.getOptimalDFTSize()

import cv2import numpy as npfrom matplotlib import pyplot as pltimg = cv2.imread(r'E:\image\test17.png', 0)rows, cols = img.shapeprint(rows, cols)# 计算DFT效率最佳的尺寸nrows = cv2.getOptimalDFTSize(rows)ncols = cv2.getOptimalDFTSize(cols)print(nrows, ncols)nimg = np.zeros((nrows, ncols))nimg[:rows, :cols] = imgimg = nimg# OpenCV计算快速傅里叶变换,输入图像应首先转换为np.float32,然后使用函数cv2.dft()和cv2.idft()。# 返回结果与Numpy相同,但有两个通道。第一个通道为有结果的实部,第二个通道为有结果的虚部。dft = cv2.dft(np.float32(img), flags=cv2.DFT_COMPLEX_OUTPUT)dft_shift = np.fft.fftshift(dft)magnitude_spectrum = 20 * np.log(cv2.magnitude(dft_shift[:, :, 0], dft_shift[:, :, 1]))plt.subplot(121), plt.imshow(img, cmap='gray')plt.title('Input Image'), plt.xticks([]), plt.yticks([])plt.subplot(122), plt.imshow(magnitude_spectrum, cmap='gray')plt.title('Magnitude Spectrum'), plt.xticks([]), plt.yticks([])plt.show()rows, cols = img.shapecrow, ccol = rows // 2, cols // 2# 首先创建一个mask,中心正方形为1,其他均为0# 如何删除图像中的高频内容,即我们将LPF应用于图像。它实际上模糊了图像。# 为此首先创建一个在低频时具有高值的掩码,即传递LF内容,在HF区域为0。mask = np.zeros((rows, cols, 2), np.uint8)mask[crow - 30:crow + 30, ccol - 30:ccol + 30] = 1# 应用掩码Mask和求逆DTFfshift = dft_shift * maskf_ishift = np.fft.ifftshift(fshift)img_back = cv2.idft(f_ishift)img_back = cv2.magnitude(img_back[:, :, 0], img_back[:, :, 1])plt.subplot(121), plt.imshow(img, cmap='gray')plt.title('Input Image'), plt.xticks([]), plt.yticks([])plt.subplot(122), plt.imshow(img_back, cmap='gray')plt.title('Magnitude Spectrum'), plt.xticks([]), plt.yticks([])plt.show()

版权

本文链接地址:https://www.jiuchutong.com/zhishi/300297.html 转载请保留说明!

上一篇:第09章_异常处理(异常处理流程为哪几个部分)

下一篇:疑问搞懂,python中文词频统计,让你真能学会(python答疑)

  • 增值税专用发票几个点
  • 息税前利润增长率与财务杠杆系数
  • 一般纳税人材料销售税率
  • 软件企业两免三减半备案条件
  • 资产负债表中的货币资金怎么算
  • 交易性金融资产入账价值怎么计算
  • 项目所在地怎么填
  • 防伪税控维护费会计分录
  • 营改增后开不了增值税发票怎么样办?
  • 视频制作费属于劳务费吗
  • 没进项发票要交多少税
  • 开增值税专票交多少税费
  • 小规模纳税人未建账处罚
  • 企业创业初期能干什么
  • 折旧费和摊销费怎么算
  • 所得税季报季末从业人数怎么填
  • 没有成立工会的企业要交工会筹备金吗
  • 怎么核算外汇业务
  • 个体户生产经营所得怎么报税
  • 商品和安装能开什么发票
  • 苹果mac怎么复制文件到u盘 免费软件
  • 企业分红需要缴纳什么税
  • 光伏电站运维费用清单
  • scanexplicit.exe - scanexplicit是什么进程 作用是什么
  • php类型判断
  • dwm.exe是啥
  • 2020年前端面试
  • PHP:mcrypt_enc_get_iv_size()的用法_Mcrypt函数
  • 工伤索赔的赔偿金额
  • php操作步骤
  • 资本公积转增注册资本的条件
  • 行政事业单位 银行
  • 工业企业制造费用占比
  • 民办非企业单位登记管理暂行条例
  • 为什么那么多人考公务员
  • 固定资产到期日
  • 库存商品的盘盈盘亏
  • 计提增值税怎么计提
  • 加计扣除内容
  • 公司捐款怎样抵扣
  • 公司出资认缴和实缴的区别
  • 怎么安装vue-cli脚手架
  • 公司车辆出售要交多少税
  • 商品流通企业的基本要素
  • 年底进项税额比年初多
  • 租车费开票属于哪个类别
  • 用友T3财务报表没有数据
  • mysql的主从
  • 购买方收到红字发票怎么做凭证
  • 电子税务局实名核验失败
  • 暂估原材料,材料成本如何冲
  • 小规模纳税人出租不动产免征增值税
  • 其他应收款审计说明怎么写
  • 机动车发票红冲多了税盘被锁
  • 小规模纳税人减按1%政策
  • 主营业务收入借方表示
  • 装修费一定要摊销吗
  • 什么是暂估入库核算
  • 福利费的概念是指
  • 开出去的发票
  • 超市会计怎么做会计分录
  • debian安装zsh
  • win7关闭窗口的方式有哪些
  • win8无法安装其它软件
  • DWHeartbeatMonitor.exe进程是什么意思 DWHeartbeatMonitor是安全的进程吗
  • win10安装不了ie
  • windows 8升级
  • win10的qq
  • win7系统不兼容开不了机
  • js倒计时结束操作
  • rhel7.6安装
  • angular form表单
  • js跨域是什么意思
  • jquery过滤选择器按照过滤规则分类包括?
  • android入门基础知识
  • python最优化
  • 研发费用加计扣除75%还是100%
  • 外埠企业如何在经营地缴税
  • 电子税务局登录不上,显示用户名不匹配
  • 员工持股平台税务稽查
  • 免责声明:网站部分图片文字素材来源于网络,如有侵权,请及时告知,我们会第一时间删除,谢谢! 邮箱:opceo@qq.com

    鄂ICP备2023003026号

    网站地图: 企业信息 工商信息 财税知识 网络常识 编程技术

    友情链接: 武汉网站建设