[发明专利]一种基于视通路多感受野朝向性关联的轮廓检测方法

摘要

本发明公开了一种基于视通路多感受野朝向性关联的轮廓检测方法，本发明从视通路中各层级的感受野空间尺度特性以及同一层级感受野的相互关联性出发，模拟视通路中视网膜神经元的经典感受野特性和LGN(外侧膝状体)细胞的非经典感受野特性，同时有朝向性地将多个LGN细胞的矩形感受野进行关联，用于模拟V1(初级视皮层)简单细胞的方向选择特性，提出一种基于视通路多感受野朝向性关联的轮廓检测方法，充分发挥上述视觉特性在轮廓检测中的重要作用。

主权项

1.一种基于视通路多感受野朝向性关联的轮廓检测方法，其特征在于，该方法具体包括以下步骤：步骤(1)分别构建on、off型双高斯差全局模板对输入的灰度图像IG(x,y)进行空间滤波，即把函数值作为空间分布权值，通过二维卷积计算图像像素点所对应的空间线性总和；随后对滤波结果进行非负半波整流，即对滤波结果中的负数置零，非负数保持不变，进一步得到RF1的响应其中，μ表示感受野的兴奋型或抑制型极性，具体取值为on或off，下同；σ_r是RF1抑制区域高斯函数的标准差；步骤(2)建立与灰度图像相同尺寸的神经元网络GN(x,y)；网络中每个神经元的突触权值在空间距离上呈高斯分布，神经元采用LIF简化模型；将步骤(1)得到的RF1响应以一一对应的方式输到神经元网络GN(x,y)，统计一个响应周期内各神经元的脉冲发放次数，计算脉冲发放频率，作为视网膜对轮廓信息的编码结果步骤(3)由于LGN细胞与视网膜神经元相比，具有更大的感受野，即通常RF1偏好于捕捉图像的局部细节，而RF2更多地关注图像的全局特征；因此，将步骤(2)的编码结果进行局部信息融合：对于某个特定点，以它为中心构建局部窗口，计算窗口覆盖区域的平均值，得到全局特征图并将其作为RF2中心区域的输入；如式(1)所示；其中，ws是局部方窗的窗长；step为移动步长，设置与窗长相等；此时全局特征图的大小与原图相比，将有所缩小；考虑到LGN细胞nCRF大外周区对感受野中心区域的调节作用，构建三高斯差函数的二次半波整流模型TDoG_HR(x,y)，如式(2)所示；式中，Temp(x,y)是第一次半波整流结果，σ_sur、σ_inh与σ_exc分别是RF2大外周区域、周边抑制区域以及中心兴奋区域的高斯函数标准差，根据式(3)获取RF2大外周区对中心兴奋区域输入的调节量并将全局特征图与负指数化后的调节量以权值α相减，用双线性插值方法把相减结果调整为原灰度图大小，即为RF2的响应如式(4)所示；其中resize表示双线性插值放大运算，δ是负指数化的强度参数，取0.2；nCRF表示非经典感受野，默认权值α＝0.6；步骤(4)构建N个on型与N个off型，共2N个尺寸为w×w的LGN细胞矩形感受野RF3，其中1≤N≤round(min(m,n)/w)，m、n分别表示输入灰度图像的宽和高，round为取整运算，min表示取较小值；所构建的2N个RF3在空间位置上相邻分布，并且呈现特定的朝向性关联：首先，在空间上任意确定一个位置点作为2N个RF3的综合重心，记作comCenter，并围绕comCenter建立一个三元素的元胞集合SetRF3来表示2N个RF3的相对位置，如式(5)所示；式中，μ_i为第i个RF3的极性，r_i为第i个RF3中心位置到comCenter的距离，是第i个RF3中心位置与comCenter的连线相对于水平方向的夹角；显然极坐标确定了2N个RF3的相对位置与整体朝向性；其次，每个RF3分别按照其在朝向性关联中的相对位置，从RF2响应中获取具体输出值，并记作如式(6)所示；其中，表示极坐标转换到直角坐标后每个RF3的输出值相对于RF2响应的位置偏移量；步骤(5)基于视通路中感受野的层级关联特性——V1简单细胞由若干相邻的LGN细胞同心圆式感受野的输出相连而成，根据步骤(4)中2N个RF3的朝向性关联输出，建立V1方位敏感性感受野RF4的多感受野朝向性关联模型，用于实现V1简单细胞的方向选择特性：将步骤(4)确定的综合重心comCenter作为RF4的中心，同时计算2N个RF3输出的几何平均值作为RF4的方向选择结果osθ(x,y；σr),如式(7)所示；其中，θ表示RF4的偏好方向，由集合SetRF3确定；V1表示初级视皮层；步骤(6)考虑到不同偏好方向的RF4结构相似而输入来源不同，说明对于其他偏好方向的RF4，同样具有2N个RF3输入，只是这些RF3采取了其他朝向的关联，因此，通过调整集合SetRF3中每个元胞的极坐标方向进一步获取其他RF4的方向选择结果如式(8)所示；其中，θ_k表示其他偏好方向，k表示方向索引，k＝1,...,Orens，Orens为方向个数，在空间中均匀分布；为朝向性关联时的角度增量；步骤(7)根据RF4的方向选择结果针对每个(x,y)，选择当前位置在RF4全部方向选择结果中的最大值，重构得到初级视皮层轮廓响应RV1(x,y；σ_r)，同时记录最大值所在的方向，得到最优方向矩阵BO(x,y)；步骤(8)根据最优方向矩阵BO(x,y)，对初级视皮层轮廓响应RV1(x,y；σr)做非极大值抑制处理，得到轮廓响应的细化结果thin(x,y；σr)；步骤(9)对细化结果做滞后阈值处理，得到轮廓响应二值图BC(x,y；σr)；Retina为视网膜层，LGN为外侧膝状体层，V1为初级视皮层，对感受野作如下标记：Retina层的经典感受野——记为RF1；LGN层的非经典感受野——记为RF2；LGN层的矩形感受野——记为RF3；V1方位敏感性感受野——记为RF4。