时域下的一维Gabor滤波器:

可以将Gabor滤波器看作是两个输出两个相位的滤波器,分别在实数域和虚数域上。

实数域上滤波器为:

虚数域上滤波器为:

傅里叶变换为频域:

上述两个滤波器对频率敏感,为了获得一个对相位不敏感,而且对正弦输入无正面响应的响应,因此,通过计算两个滤波器输出的和的平方根来得到这样一个滤波器。

在频率域上,针对一个特定频率的响应能量只是傅里叶变换的能量。       

上述式子为gaussian 函数,中心为fo,宽度比例为a。

频率带宽和峰值响应

峰值响应在频率为fo处,

半个能量级的频率为

因此

所以带宽为约等于 2*0.5a = a,所以a可以被解释为半能量滤波器的带宽。

直流电响应的消除

通过调节fo 和a可以获得较大的直流响应。为了排除直流响应通过减去一个低通高斯滤波器的输出操作。

傅里叶变换后:

通过一系列推导,得到当b=a时,满足条件,解得:

空间二维滤波

函数:,其中s(x,y)为载波(正弦函数),w_r(x,y)为二维高斯形状的滤波函数,如包络线。

载波

(u0,v0)定义空间频率,P定义正弦曲线的相位

可以将载波函数看作实轴和虚轴上的两个函数

u0,v定义正弦曲线在笛卡尔坐标系下的空间频率。空间频率也可以在极坐标下描述为能量F0和方向w0

ie

载波通过上述变换可以得到

Gaussian 包络线

其中(x0,y0)为此函数的峰值,a,b为gaussian 尺度大小调节参数,下标r表示循环操作

形状如下:

gaussian 函数主要由一下九个参数定义:

每一个Gabor函数都有两个正交的函数,分别为函数的实轴和虚轴部分。

形状如下:

空域上的gabor函数为:

极坐标下的gabor函数为:

傅里叶变换以后的Gabor函数显示:

半能级描述

由于其大小等于峰值能级的一半,而峰值能级(极值)可以由以下获得,峰值(极值点)为(u,v) = (u0,v0),得到峰值能级为K/ab,

因此,只需找到一系列满足K/2ab的(u,v)的值。

半能级频率和方向带宽

神经元半能级频率及方向带宽可以用于针对半能级响应的测量,假设优选神经元的空间频率u0,v0,

在极坐标下空间频率可以被表述为F0,w0.

我们通过用正弦图像的方向,不同空间频率能级来探测神经元,进而发现半能级频率带宽。

通过增加F,直到能级的响应是在(F0,w0)处能级的一半。定义此时为F_max,

通过较小F,直到能级的响应是在(F0,w0)处能级的一半。定义此时为F_min

定义半能级频率宽度为:

或者在倍频下为:

半能级方向带宽和半能级频率带宽相似,只是这里操作的是由F变为w

得到如下结果:

有效传播和均方根传播

一函数均方根的长度,宽度,和区域定义在他们的一阶矩和二阶矩

一个复杂函数g(x,y)的矩阵通过将这个函数转为概率密度(必须为正且和为1)然后计算其一阶矩和二阶矩(一阶矩就是随机变量的期望,二阶矩就是随机变量平方的期望)

通过g(x,y),我们建立如下概率密度

g(x,y)绝对值的平方为信号的能量级的平方,总是正的

Z保证累加和为1.0

统计度量的位置和大小:

均方根的长度和宽度定义为一个旋转型的f(x,y)的,因此上述的旋转分布的协方差为0 。

让Xr Yr代表 协方差为0的旋转变量,均方根的长度和宽度为

在频率域上通过傅里叶变换,得到的定义如下:

在空域和频率域上的均方根区域定义如下:

有效长度,宽度,区域只是乘上一个sqrt(2)*pi

在任何二维函数的限定条件下,有如下不等式:

用Gabor函数可以很简单的证明这不确定关系的下限。对于在空间域上给定的一个区域,提供了在频率域上最大可能的分辨率,反之亦然。

Gabor函数的均方根宽度和长度为:

将Gabor函数的概率密度简单的考虑为Gaussian

方差为

由不确定关系可以得到
                             

作为简单细胞感受野模型的Gabor函数(此章含有大量生物上由Gabor思想启发的设计文献,忽略此章的阅读)

7Gabor的空域滤波

考虑大量带有Gabor核函数及除了位置参数(x,y)的其他参数的感受野神经元。让这些神经元呈中心凹处均匀分布。

中心凹处的至少包含两个正交的神经元。我们可以模型化将一系列神经元的操作作为卷积操作(假设在中心凹处一个连续均匀分布的滤波器)。由于空域卷积是频率域的结果。Gabor的功能像是对中心凹处图像进行带通滤波。峰频由载波函数的空间频率控制(U0,V0)半能级区由高斯包络线的旋转角度,尺度参数a和b控制。

8能量滤波一个正交对是一系列两个具有相同振幅响应但相位响应移动90度的两个线性操作。严格的讲,sin,和cos的Gabor操作不是整觉对,因为cos相位Gabor具有一些直流响应。然而sin Gabor却没有。然而,可以让正交Gabor对看起来非常像sin/cos对。因此,sin和cos Gabor对被普遍当作是正交对。

能量机械学是累加正交对输出的平方的系统。对于随机正弦输入 ,能量机械有非调制响应。

9~end 只大体浏览未重点阅读

论文阅读笔记一:Tutorial on Gabor Filters的更多相关文章

  1. 论文阅读笔记 - YARN : Architecture of Next Generation Apache Hadoop MapReduceFramework

    作者:刘旭晖 Raymond 转载请注明出处 Email:colorant at 163.com BLOG:http://blog.csdn.net/colorant/ 更多论文阅读笔记 http:/ ...

  2. 论文阅读笔记 - Mesos: A Platform for Fine-Grained ResourceSharing in the Data Center

    作者:刘旭晖 Raymond 转载请注明出处 Email:colorant at 163.com BLOG:http://blog.csdn.net/colorant/ 更多论文阅读笔记 http:/ ...

  3. 论文阅读笔记 Word Embeddings A Survey

    论文阅读笔记 Word Embeddings A Survey 收获 Word Embedding 的定义 dense, distributed, fixed-length word vectors, ...

  4. 论文阅读笔记 Improved Word Representation Learning with Sememes

    论文阅读笔记 Improved Word Representation Learning with Sememes 一句话概括本文工作 使用词汇资源--知网--来提升词嵌入的表征能力,并提出了三种基于 ...

  5. [置顶] 人工智能(深度学习)加速芯片论文阅读笔记 (已添加ISSCC17,FPGA17...ISCA17...)

    这是一个导读,可以快速找到我记录的关于人工智能(深度学习)加速芯片论文阅读笔记. ISSCC 2017 Session14 Deep Learning Processors: ISSCC 2017关于 ...

  6. Nature/Science 论文阅读笔记

    Nature/Science 论文阅读笔记 Unsupervised word embeddings capture latent knowledge from materials science l ...

  7. 论文阅读笔记(二十一)【CVPR2017】:Deep Spatial-Temporal Fusion Network for Video-Based Person Re-Identification

    Introduction (1)Motivation: 当前CNN无法提取图像序列的关系特征:RNN较为忽视视频序列前期的帧信息,也缺乏对于步态等具体信息的提取:Siamese损失和Triplet损失 ...

  8. 论文阅读笔记(十八)【ITIP2019】:Dynamic Graph Co-Matching for Unsupervised Video-Based Person Re-Identification

    论文阅读笔记(十七)ICCV2017的扩刊(会议论文[传送门]) 改进部分: (1)惩罚函数:原本由两部分组成的惩罚函数,改为只包含 Sequence Cost 函数: (2)对重新权重改进: ① P ...

  9. [论文阅读笔记] GEMSEC,Graph Embedding with Self Clustering

    [论文阅读笔记] GEMSEC: Graph Embedding with Self Clustering 本文结构 解决问题 主要贡献 算法原理 参考文献 (1) 解决问题 已经有一些工作在使用学习 ...

  10. [论文阅读笔记] metapath2vec: Scalable Representation Learning for Heterogeneous Networks

    [论文阅读笔记] metapath2vec: Scalable Representation Learning for Heterogeneous Networks 本文结构 解决问题 主要贡献 算法 ...

随机推荐

  1. 基于ip的虚拟主机配置——在一台服务器上绑定多个 IP 地址

    进入/etc/sysconfig/network-scripts,修改ifcfg-ens33文件 输入 ip addr 查看ip 引用:https://blog.csdn.net/u013887008 ...

  2. 获取AWR报告

    1.进入sqlplus [oracle@localhost admin]$ sqlplus / as sysdba SQL Production :: Copyright (c) , , Oracle ...

  3. CF1095E Almost Regular Bracket Sequence

    题目地址:CF1095E Almost Regular Bracket Sequence 真的是尬,Div.3都没AK,难受QWQ 就死在这道水题上(水题都切不了,我太菜了) 看了题解,发现题解有错, ...

  4. 1.Python_字符串_常用办法总结

    明确:对字符串的操作方法都不会改变原来字符串的值. 1.去掉空格和特殊符号 name.strip() 去掉空格和换行符 name.strip("xx") 去掉某个字符串 name. ...

  5. opencv处理验证码python代码

    # -*- coding: utf-8 -*- # @Time : 2019-02-11 09:39 # @Author : cxa # @File : bgr2gry.py # @Software: ...

  6. unity制作背景

  7. Java的if判断对象为null时,null放在比较运算符的左边还是右边较好?

    如java中:if(name == null)和if(null == name)有什么讲究吗? 答:在java里面,它们是一样的.但是通常写为null == name.这其实是在C语言里面引申出来的. ...

  8. $Django content_type组件 缓存组件

    1  content_type组件(只能方便插入添加) 需求:课程,学位课(不同的课程字段不一样),价格策略 #免费课 class Free_classes (models.Model): id = ...

  9. python操作三大主流数据库(2)python操作mysql②python对mysql进行简单的增删改查

    python操作mysql②python对mysql进行简单的增删改查 1.设计mysql的数据库和表 id:新闻的唯一标示 title:新闻的标题 content:新闻的内容 created_at: ...

  10. python连接kafka生产者,消费者脚本

    # -*- coding: utf-8 -*- ''''' 使用kafka-Python 1.3.3模块 # pip install kafka==1.3.5 # pip install kafka- ...