本发明属于视觉神经通路技术领域,具体涉及一种基于视觉神经通路的边缘检测方法。
背景技术:
作为数字图像处理和计算机视觉的基础性工作,图像中的边缘被定义为像素及其邻域在特定方向上的阶跃性变化(颜色变化或灰度变化),边缘检测则是通过特定算法,标注出这些变化明显的像素,从而为后续的图像处理提供关键信息。
图像的边缘检测大幅度的减少了图像包含的数据量,剔除了冗余的、不相关的信息,保留图像重要的轮廓性属性,从而为目标识别和图像理解等较高层次的计算机视觉应用提供有效信息。因此,如何准确检测图像中的边缘,成为图像处理及相关研究领域的重要研究问题。
边缘检测的研究具有悠久的历史,迄今有许多边缘检测算法被提出,如一阶微分算子、二阶微分算子、log算子、dog算子、canny算子等,以及基于模糊理论、数学形态学及神经网络等新兴技术理论的方法。但是这些方法大多忽视了一个基本问题,即边缘检测之所以重要,主要依赖于人类视觉理论的假设,该假设认为人类对物体的感知在很大程度上依赖于边缘的。因此,只有符合人眼视觉机制的边缘检测算法,才可以检测出人眼认为准确的边缘。
技术实现要素:
本发明的目的是提供一种基于视觉神经通路的边缘检测方法,解决了现有技术中存缺少符合人眼视觉机制的边缘检测算法问题。
本发明所采用的技术方案是,一种基于视觉神经通路的边缘检测方法,具体按照以下步骤实施:
步骤1、建立视网膜的感受机制模型,模拟视网膜对图片信息的响应;
步骤2、建立侧膝体的感受机制模型,模拟侧膝体对图片信息的响应;
步骤3、建立大脑皮层v1区的感受机制模型,模拟脑皮层v1区对边缘轮廓的检测;
步骤4、结合步骤1~3建立的模型组合生成基于视觉神经通路的边缘检测模型,模拟视觉神经通路的边缘检测输出。
本发明的特点还在于:
步骤1具体按照以下步骤实施:
步骤1.1、建立光源信号对电信号的转换模型,即在视网膜神经外丛状层中水平细胞和双极细胞中进行光源信号对电信号的转换,具体如式(1)所示:
式中,ωi代表以xi为中心的双极细胞的感受野,k(x,t)是单个双极细胞在t时刻的时空传播模型,f(x,t)是t时刻感光细胞接受到的光信号刺激,式(1)表示外丛状层电信号的产生完全线性依赖于以xi为中心的双极细胞感受野ωi及t时刻以前的输入信号在感受野中的滞留;
步骤1.2、单个双极细胞时空传播模型k(x,t)简化为所有双极细胞时空传播模型k1(x,t),假设双极细胞数量极多且所有双极细胞均遵循同样的时空传播函数k(x,t),且时空传播函数k(x,t)不受限于ωi空间内,则式(1)可简化为:
λ(x,t)=∫t′∈r∫x′∈rk1(x-x′,t-t′)f(x′,t′)dx′dt′=(k1*f)(x,t)(2)
步骤1.3、光信号刺激f(x,t)具有时变性,并且中枢和周围反应具有不用的时间进程,故建立时空不可分离的中心-环绕型感受野模型,具体如式(3)~(5)所示:
k(x,t)=c(x,t)-s(x,t)(3)
式中,ωc是感受野中心区的常数系数,ωs是周围区的常数系数,g(x)是高斯滤波器,σc和σs是用于调节感受野半径的高斯函数标准差;
光源信号对电信号的转换模型、双极细胞时空传播模型k1(x,t)及时空不可分离的中心-环绕型感受野模型组合为视网膜的感受机制模型。
步骤1.3中v(t)是低通滤波器函数,具体如式(6)所示:
式中,δ0是单位冲击函数;
步骤2具体按照以下步骤实施:
步骤2.1、建立侧膝体的感受机制模型,侧膝体的感受机制模型包括经典感受野模型与非经典感受野模型,典感受野模型与非经典感受野模型均为加权的高斯函数,侧膝体模型具体如式(8)所示:
b(x,t)=h(x,t)*a(x,t)(8)
式中,h(x,t)表示外膝体细胞在t时刻对视网膜输出信号a(x,t)的变换;
外膝体细胞在t时刻对对视网膜输出信号a(x,t)的变换h(x,t),具体如(9)所示:
式中,rc表示经典感受野中心区域,rs表示经典感受野周围区域,rns表示非经典感受野的半径,rs≥2rc,rns≥4rc,
步骤2.2、简化侧膝体的感受机制模型,对于静态图像,假设任意时刻内的侧膝体变换是不变的,获得简化侧膝体的感受机制模型,则式(8)可简化为:
b(x)=h(x)*a(x)(10)
简化侧膝体的感受机制模型即为侧膝体的感受机制模型。
步骤3具体按照以下步骤实施:
步骤3.1、大脑皮层v1区的感受机制模型选用感受野合并模型;
步骤3.2、感受野合并模型视为简单细胞由多个子单元组成,每个子单元接受侧膝体细胞的输出作为输入,并对其附近的侧膝体细胞响应线性加权求和计算出每个子单元的响应,如式(12)~(14)所示:
式中,δi,σi,ρi,
步骤3.3、基于感受野合并模型中简单细胞对特定方向的响应,定义为所有子单元响应的几何加权平均,具体式(15)所示,即只有当简单细胞的所有子单元都被激活时,该简单细胞才被激活;
ωi是由子单元位置参数决定的系数;
步骤3.4、感受野合并模型的不同方向选择性通过改变子单元位置的参数来得到,使用具有不同方向的感受野合并模型对图像进行处理,即可以获得各个方向的边缘,并使用最大化操作方法对各个方向的简单细胞响应进行合并,获得边缘响应集合,具体如式(16)所示,针对得到简单细胞的边缘响应集合,使用canny算子中的非极大值抑制;
式中,ψ(c)是一系列不同方向角度对应的子单元位置参数集合,感受野合并模型即为大脑皮层v1区的感受机制模型。
步骤4具体按照以下步骤实施:
步骤4.1、根据步骤1~3建立的视网膜的感受机制模型、侧膝体的感受机制模型及大脑皮层v1区的感受机制模型,将视网膜的感受机制模型、侧膝体的感受机制模型及大脑皮层v1区的感受机制模型进行顺序组合生成基于视觉神经通路的边缘检测模型;
步骤4.2、向基于视觉神经通路的边缘检测模型输入感受机制模型、侧膝体的感受机制模型及大脑皮层v1区的感受机制模型的主要参数;
步骤4.3、基于视觉神经通路的边缘检测模型使用python和opencv进行计算机程序仿真模拟,输入一张彩色图片,输出边缘图像。
本发明的有益效果是:
1、针对整个视觉神经通路建模,而不仅关注视觉神经通路中的某一部分。
2、在最新生物学研究成果上,对视网膜、侧膝体细胞响应进行建模。
3、可以得到更为符合人眼视觉感受的边缘检测图。
附图说明
图1是本发明一种基于视觉神经通路的边缘检测方法的流程图;
图2是本发明一种基于视觉神经通路的边缘检测方法的视网膜细胞结构分布图;
图3是本发明一种基于视觉神经通路的边缘检测方法的侧膝体非经典感受野分布图;
图4是生物学研究发现的侧膝体细胞反应曲线图;
图5是本发明一种基于视觉神经通路的边缘检测方法的侧膝体的感受机制模型的反应曲线图;
图6是本发明一种基于视觉神经通路的边缘检测方法的简单细胞的子单元分布图;
图7是本发明实施例1的输入图片;
图8是本发明实施例1的边缘输出图像。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明一种基于视觉神经通路的边缘检测方法,如图1所示,具体按照以下步骤实施:
步骤1、建立视网膜的感受机制模型,模拟视网膜对图片信息的响应;
步骤2、建立侧膝体的感受机制模型,模拟侧膝体对图片信息的响应;
步骤3、建立大脑皮层v1区的感受机制模型,模拟脑皮层v1区对边缘轮廓的检测;
步骤4、结合步骤1~3建立的模型组合生成基于视觉神经通路的边缘检测模型,模拟视觉神经通路的边缘检测输出。
通过生物学及神经科学对视网膜的研究,使用数学工具对人眼视网膜细胞响应进行建模,并在建模过程中考虑水平细胞及双极细胞的时延特性。
通过生物学及神经科学对侧膝体的研究,使用数学工具对下丘脑中的侧膝体细胞响应进行建模,并在建模过程中考虑侧膝体非经典感受野对信号的调制作用。
通过对视网膜及侧膝体的建模,使用该模型对图片进行处理,即可得到模拟视网膜和侧膝体对图片信息的响应,并将其作为输入,传输到基于大脑皮层v1区响应建立的模型中,由于脑皮层v1区主要负责的就是边缘轮廓的检测,通过对v1区建模得到的输出,即为模拟视网膜通路得到的边缘检测输出。
步骤1具体按照以下步骤实施:
根据生物学及神经科学对眼的研究,对光信号在眼中的传输进行建模。首先,人眼接收光源信号,光源信号在视网膜中经过一系列变换以适应人脑的需求。而视网膜神经细胞主要可以分为三类,外丛状层、内丛状层及神经节层。外丛状层接受光源信号并将其线性转换为电信号,内丛状层完成电信号的非线性矫正,神经节层则用于产生一系列基于电信号幅度的脉冲信号。
步骤1.1、建立光源信号对电信号的转换模型,即在视网膜神经外丛状层中水平细胞和双极细胞中进行光源信号对电信号的转换,具体如式(1)所示:
式中,ωi代表以xi为中心的双极细胞的感受野,k(x,t)是单个双极细胞在t时刻的时空传播模型,f(x,t)是t时刻感光细胞接受到的光信号刺激,式(1)表示外丛状层电信号的产生完全线性依赖于以xi为中心的双极细胞感受野ωi及t时刻以前的输入信号在感受野中的滞留;
步骤1.2、单个双极细胞时空传播模型k(x,t)简化为所有双极细胞时空传播模型k1(x,t),假设双极细胞数量极多且所有双极细胞均遵循同样的时空传播函数k(x,t),且时空传播函数k(x,t)不受限于ωi空间内,则式(1)可简化为:
λ(x,t)=∫t′∈r∫x′∈rk1(x-x′,t-t′)f(x′,t′)dx′dt′=(k1*f)(x,t)(2)
步骤1.3、光信号刺激f(x,t)具有时变性,并且中枢和周围反应具有不用的时间进程,故建立时空不可分离的中心-环绕型感受野模型,具体如式(3)~(5)所示:
k(x,t)=c(x,t)-s(x,t)(3)
式中,ωc是感受野中心区的常数系数,ωs是周围区的常数系数,g(x)是高斯滤波器,σc和σs是用于调节感受野半径的高斯函数标准差;
步骤1.3中v(t)是低通滤波器函数,具体如式(6)所示:
式中,δ0是单位冲击函数;
光源信号对电信号的转换模型、双极细胞时空传播模型k1(x,t)及时空不可分离的中心-环绕型感受野模型组合为视网膜的感受机制模型。
对视网膜的感受机制模型的时间信息解释如下:感光细胞接收到光源信号后,于τg时刻先传送给水平细胞,同时,一小部分感光细胞在τc时刻将光源信号直接传动给双极细胞,这两部分信号加权后用于激活双极细胞的感受野;此外,τg时刻水平细胞接收到的信号经过整合交互后,再经过τs的延时传送给双极细胞,这部分信号对双极细胞起到抑制作用。综上,双极细胞就接收到两种极性相反的信号,在考虑到时间因素的基础上,使用高斯函数对感受野的中心-环绕型特性进行建模,即可得到时空不可分离的中心-环绕型感受野模型。由于图像在任意时刻是静态不变的,光信号刺激f(x,t)可直接表示为。因此直接用k(x,t)和f(x)进行卷积即可得到视网膜输出的电信号a(x,t)。
步骤2具体按照以下步骤实施:
与视网膜神经节细胞的感受野相似,侧膝体细胞同样具有感受野机制,且其经典的感受野也被认为是中心-环绕型的模型。但是,根据生物学对哺乳动物的研究发现,在侧膝体经典感受野的周围,还存在一个区域,该区域不会产生响应,但是会对位于经典感受野的响应产生调节作用,即在经典感受野外再添加一个刺激,细胞的响应会产生改变。
基于相关生物学及生理学研究提出的反应曲线如图4所示,本方案提出侧膝体的感受机制模型以包括的经典感受野模型与非经典感受野模型。侧膝体的感受机制模型假设非经典感受野是位于经典感受野周围的另一个同心圆,即经典感受野外存在一个大外环表示非经典感受野,如图3所示。
步骤2.1、建立侧膝体的感受机制模型,侧膝体的感受机制模型包括经典感受野模型与非经典感受野模型,典感受野模型与非经典感受野模型均为加权的高斯函数,侧膝体模型具体如式(8)所示:
b(x,t)=h(x,t)*a(x,t)(8)
式中,h(x,t)表示外膝体细胞在t时刻对对视网膜输出信号a(x,t)的变换;
外膝体细胞在t时刻对对视网膜输出信号a(x,t)的变换h(x,t),具体如(9)所示:
式中,rc表示经典感受野中心区域,rs表示经典感受野周围区域,rns表示非经典感受野的半径,rs≥2rc,rns≥4rc,
步骤2.2、简化侧膝体的感受机制模型,对于静态图像,假设任意时刻内的侧膝体变换是不变的,获得简化侧膝体的感受机制模型,则式(8)可简化为:
b(x)=h(x)*a(x)(10)
简化侧膝体的感受机制模型即为侧膝体的感受机制模型。
本方案的侧膝体的感受机制模型如图5所示,使用均方根,如式(11),统计本发明建模结果与图4中反应曲线的误差为0.017,即从数据上验证本发明对lgn非经典感受野所建模型的有效性。
h(xi)为本模型每个取样点xi对应的幅值,r(xi)为图4中反应曲线每个取样点对应的幅值。
步骤3具体按照以下步骤实施:
生物学相关研究发现,哺乳动物的大脑皮层v1区细胞不会对光点刺激产生响应,但是对特定方向的光栅十分敏感。除此之外,不同细胞的感受野特性也会有所区别,可被分为简单细胞,复杂细胞,和超复杂细胞等。其中,简单细胞会对特定方向和宽度的条状刺激做出强烈反应,因此,许多研究认为,简单细胞的感受野从结构上可以认为是多个侧膝体细胞按照特定的方向聚合而成。根据这一结构提出的感受野合并模型在本发明中进行应用,作为对整个视通路中大脑皮层v1区的模型。
步骤3.1、大脑皮层v1区的感受机制模型选用感受野合并模型;
步骤3.2、感受野合并模型视为简单细胞由多个子单元组成,其排布如图3所示,每个子单元接受侧膝体细胞的输出作为输入,并对其附近的侧膝体细胞响应线性加权求和计算出每个子单元的响应,如式(12)~(14)所示:
式中,δi,σi,ρi,
步骤3.3、基于感受野合并模型中简单细胞对特定方向的响应,定义为所有子单元响应的几何加权平均,具体式(15)所示,即只有当简单细胞的所有子单元都被激活时,该简单细胞才被激活;
ωi是由子单元位置参数决定的系数;
步骤3.4、感受野合并模型的不同方向选择性通过改变子单元位置的参数来得到,使用具有不同方向的感受野合并模型对图像进行处理,即可以获得各个方向的边缘,并使用最大化操作方法对各个方向的简单细胞响应进行合并,获得边缘响应集合,具体如式(16)所示,针对得到简单细胞的边缘响应集合,使用canny算子中的非极大值抑制;
式中,ψ(c)是一系列不同方向角度对应的子单元位置参数集合,感受野合并模型即为大脑皮层v1区的感受机制模型。
步骤4具体按照以下步骤实施:
步骤4.1、根据步骤1~3建立的视网膜的感受机制模型、侧膝体的感受机制模型及大脑皮层v1区的感受机制模型,将视网膜的感受机制模型、侧膝体的感受机制模型及大脑皮层v1区的感受机制模型进行顺序组合生成基于视觉神经通路的边缘检测模型;
步骤4.2、向基于视觉神经通路的边缘检测模型输入感受机制模型、侧膝体的感受机制模型及大脑皮层v1区的感受机制模型的主要参数,如表1所示;
步骤4.3、基于视觉神经通路的边缘检测模型使用python和opencv进行计算机程序仿真模拟,输入一张彩色图片,如图7所示,输出边缘图像,如图8所示。
表1
本发明的优点在于,基于视觉神经通路的边缘检测方法分别对视网膜、侧膝体及大脑皮层v1区建模,将人眼视觉机制融入边缘检测中,可以得到更为符合人眼观察特征的边缘检测效果。另一个优点在于,本发明考虑视网膜双极细胞及水平细胞信号传输的延时性,考虑侧膝体细胞的非经典感受野对经典感受野的调制作用,考虑大脑皮层v1区细胞对特定方向宽度的光源选择性,对现有模型进行更为符合生物特性的建模。