本发明提供一种基于X、Y并行视觉通道响应的仿生型轮廓检测方法，包括以下步骤：A、输入经灰度处理的待检测图像，构建视网膜神经节细胞与外膝体中的X型细胞和Y型细胞仿真模型；B、X型细胞的初始响应值的计算；C、Y型细胞的初始响应值的计算；D、X、Y型细胞的经典感受野响应值的计算；E、X、Y型细胞的轮廓响应值的计算；F、最终轮廓响应值的计算；G、最终轮廓响应值的计算。该方法克服现有技术轮廓识别率低的缺陷，具有仿真效果好、轮廓识别率高的特点。
1.一种基于X、Y并行视觉通道响应的仿生型轮廓检测方法，其特征在于：A、输入经灰度处理的待检测图像，并构建视网膜神经节细胞与外膝体中的X型细胞和Y型细胞仿真模型；
预设X型中心二维高斯函数、X型周边二维高斯函数、X型中心周边拮抗强度，X型周边二维高斯函数的空间标准差大于X型中心二维高斯函数，对于各像素点，分别采用X型中心二维高斯函数、X型周边二维高斯函数对其灰度值进行二维高斯滤波，分别得到各像素点的X型中心响应值和X型周边响应值，将各像素点的X型中心响应值减去其X型周边响应值与X型中心周边拮抗强度的乘积，得到各像素点的X型初始响应值；
预设Y型中心二维高斯函数、Y型周边二维高斯函数、Y型中心周边拮抗强度，Y型中心二维高斯函数的空间标准差大于X型中心二维高斯函数，Y型周边二维高斯函数的空间标准差大于Y型中心二维高斯函数；对于各像素点，分别采用Y型中心二维高斯函数进行二维高斯滤波，得到各像素点的Y型的中心响应值；预设圆形的的亚单元区域，对于各像素点，采用X型中心二维高斯函数分别对以该像素点为中心的亚单元区域内的各像素点的Y型的中心响应值进行二维高斯滤波，并将得到的全部响应值求和取均值，得到各像素点的Y型亚单元总响应值；对于各像素点，采用Y型周边二维高斯函数对其Y型亚单元总响应值进行二维高斯滤波，得到各像素点的Y型亚单元抑制响应值；对于各像素点，将其Y型亚单元总响应值减去其Y型亚单元抑制响应值与Y型中心周边拮抗强度的乘积，得到各像素点的Y型初始响应值；
D、X、Y型细胞的经典感受野响应值的计算具体为：预设分别与X型细胞、Y型细胞相对应的包含多个方向参数的二维高斯一阶偏导函数，对于各像素点，采用二维高斯一阶偏导函数分别对其X型初始响应值、Y型初始响应值进行滤波，得到各像素点各方向参数的X型经典感受野初始响应值、Y型经典感受野初始响应值；
对于各像素点，分别在其各方向参数的X型经典感受野初始响应值及Y型经典感受野初始响应值中取最大值，分别作为该像素点的X型经典感受野响应值、Y型经典感受野响应值；
预设分别与X型细胞、Y型细胞相对应的归一化高斯差分函数、非经典感受野拮抗强度，对于各像素点，采用归一化高斯差分函数分别对其X型经典感受野响应值、Y型经典感受野响应值进行滤波，得到各像素点的X型非经典感受野响应值、Y型非经典感受野响应值；对于各像素点，将其X型经典感受野响应值、Y型经典感受野响应值分别减去X型非经典感受野响应值与非经典感受野拮抗强度的乘积、Y型非经典感受野响应值与非经典感受野拮抗强度的乘积，分别得到X型轮廓响应值、Y型轮廓响应值；
预设局部区域、局部亮度对比度函数、亮度对比度阈值，所述的局部区域为方形；对于各像素点，采用局部亮度对比度函数对以其为中心的局部区域内的各像素点的灰度值进行计算，得到各像素点的局部亮度对比度；对于各像素点，对其局部亮度对比度进行归一化处理，得到各像素点的局部亮度对比度归一值；对于各像素点，若其局部亮度对比度归一值大于或等于亮度对比度阈值，则将其X型轮廓响应值作为该像素点的最终轮廓响应值，若其局部亮度对比度归一值小于亮度对比度阈值，则将其Y型轮廓响应值作为该像素点的最终轮廓响应值；
对于各像素点，对其最终轮廓响应值进行非极大值抑制和双阈值处理，得到各像素点的最终轮廓值。
2.如权利要求1所述的基于X、Y并行视觉通道响应的仿生型轮廓检测方法，其特征在于：所述的步骤B具体如下：
所述的X型初始响应值RX(x,y)＝I(x,y)*Gxc(x,y)-αx·I(x,y)*Gxf(x,y) (1)；
其中I(x,y)为各像素点的灰度值，*为卷积运算；
σx为X型细胞的空间标准差；αx为X型中心周边拮抗强度。
3.如权利要求2所述的基于X、Y并行视觉通道响应的仿生型轮廓检测方法，其特征在于：所述的步骤C具体如下：
所述的Y型的中心响应值CY(x,y)＝I(x,y)*Gyc(x,y) (2)；
其中亚单元区域S1的半径为rs＝15σx， 表示亚单元区域S1内的各像素点的Y型的中心响应值CY(x,y)，(x',y')表示亚单元区域S1内的各像素点的坐标，N表示亚单元区域S1内的像素点数目；
Y型初始响应值RY(x,y)＝NY(x,y)-αy·NY(x,y)*Gyf(x,y) (4)；
4.如权利要求3所述的基于X、Y并行视觉通道响应的仿生型轮廓检测方法，其特征在于：所述的步骤D具体如下：
X型经典感受野初始响应值CRFX(x,y,θi)＝|RX(x,y)*VX(x,y,θi)| (5)；
Y型经典感受野初始响应值CRFY(x,y,θi)＝|RY(x,y)*VY(x,y,θi)| (6)；
X型经典感受野响应值CRFX(x,y)＝max(CRFX(x,y,θi)|i＝1,2,...Nθ) (7)；
Y型经典感受野响应值CRFY(x,y)＝max(CRFY(x,y,θi)|i＝1,2,...Nθ) (8)。
5.如权利要求4所述的基于X、Y并行视觉通道响应的仿生型轮廓检测方法，其特征在于：所述的步骤E具体如下：
X型轮廓响应值ResX(x,y)＝CRFX(x,y)-β·CRFX(x,y)*wX(x,y) (9)；
Y型轮廓响应值ResY(x,y)＝CRFY(x,y)-β·CRFY(x,y)*wY(x,y) (10)；
其中||·||1为L1范数，H(x)＝max(0,x)，β为非经典感受野拮抗强度。
6.如权利要求5所述的基于X、Y并行视觉通道响应的仿生型轮廓检测方法，其特征在于：所述的步骤F具体如下：
其中s2表示以当前像素点为中心的11*11像素的方形局部区域，f(xj,yj)表示区域s2内的各像素点的灰度值，j表示区域s2内的相对坐标，j∈[-5,5]；
其中ΔCmax(x,y)为各个像素点的局部亮度对比度ΔC(x,y)中的最大值，ΔCmin(x,y)为各个像素点的局部亮度对比度ΔC(x,y)中的最小值；
一种基于X、Y并行视觉通道响应的仿生型轮廓检测方法
技术领域
[0001] 本发明涉及图像处理领域，具体涉及一种基于X、Y并行视觉通道响应的仿生型轮廓检测方法。
背景技术
[0002] 轮廓定义目标的形状，轮廓是目标识别中的重要任务之一，而从杂乱场景中获取的目标轮廓是一项重要且相当困难的任务，主要是因为轮廓周围通常存在大量纹理背景的边缘，因此这项工作主要需要排除由于纹理区域的无意义边缘，而保留目标轮廓。提高检测率的关键在于能基于上下文将局部信息优化整合成一致的全局特征。人类视觉系统具有快速和有效的从复杂场景中提取轮廓特征的能力，有效促进了以生物特性作为启发的轮廓检测算法研究的发展。由光引起的神经元信号起始于视网膜，这些信号经过通过视网膜的神经节细胞的轴突传递至中继站外膝体(LGN)，然后由外膝体神经元直接经视放射到视皮层，形成了视觉的第一通路。生理研究表明，V1层神经元具有方位选择性，且在其经典感受野(Classical Receptive Rield，CRF)外存在非经典感受野(Non-Classical Receptive Rield，NCRF)区域，虽然单独刺激该区域没有响应，但可以对CRF具有一定的调制作用。
[0003] Rodieck于1965年提出了关于同心圆拮抗式感受野的数学模型，它由一个兴奋作用强的中心机制和一个作用较弱但面积更大的抑制性周边机制所构成。Rodieck模型又称高斯差模型，difference oftwo Gaussians，DOG。1966年Enroth-Cugell和Robson观察猫GC可以按其反应的空间-时间总和性质而划分为两类：一类细胞反应的空间总和特性大体符合Rodieck模型，即它们的感受野兴奋和抑制作用是可以线性相加的，称X型细胞；另一类GC的空间总和性质是非线性的，Rodieck模型对其不甚适用，称为Y型细胞。Enroth-Cugell et al.发现神经节细胞中一些细胞(X型细胞)的感受野空间特性近似于线性的，而另一些细胞(Y型细胞)很大程度上表现出非线性的空间特征，这些特性同样存在于外膝体细胞中。目前许多受生物学启发的轮廓检测模型都只是简单地模拟视觉系统中生理特性的一部分，在轮廓提取的过程中，轮廓信息的缺失以及纹理信息不能完全被抑制这类问题依旧存在，表明轮廓检测性能还需要进一步的提高。
发明内容
[0004] 本发明旨在提供一种基于X、Y并行视觉通道响应的仿生型轮廓检测方法，该方法克服现有技术仿真效果差、轮廓识别率低的缺陷，具有仿真效果好、轮廓识别率高的特点。
[0005] 本发明的技术方案如下：一种基于X、Y并行视觉通道响应的仿生型轮廓检测方法，包括以下步骤：A、输入待检测图像，构建视网膜神经节细胞与外膝体中的X型细胞和Y型细胞仿真模型；B、对X型细胞的初始响应值进行计算；C、对Y型细胞的初始响应值进行计算；D、结合X型细胞的初始响应值、Y型细胞的初始响应值对X、Y型细胞的经典感受野响应值进行计算；E、对X、Y型细胞的轮廓响应值进行计算；F、对最终轮廓响应值的计算；G、对最终轮廓值进行计算。
[0006] 优选地，基于X、Y并行视觉通道响应的仿生型轮廓检测方法，包括以下步骤：
[0007] A、输入经灰度处理的待检测图像，构建视网膜神经节细胞与外膝体中的X型细胞和Y型细胞仿真模型；
[0008] B、X型细胞的初始响应值的计算：
[0009] 预设X型中心二维高斯函数、X型周边二维高斯函数、X型中心周边拮抗强度，X型周边二维高斯函数的空间标准差大于X型中心二维高斯函数，对于各像素点，分别采用X型中心二维高斯函数、X型周边二维高斯函数对其灰度值进行二维高斯滤波，分别得到各像素点的X型中心响应值和X型周边响应值，将各像素点的X型中心响应值减去其X型周边响应值与X型中心周边拮抗强度的乘积，得到各像素点的X型初始响应值；
[0010] C、Y型细胞的初始响应值的计算：
[0011] 预设Y型中心二维高斯函数、Y型周边二维高斯函数、Y型中心周边拮抗强度，Y型中心二维高斯函数的空间标准差大于X型中心二维高斯函数，Y型周边二维高斯函数的空间标准差大于Y型中心二维高斯函数；对于各像素点，分别采用Y型中心二维高斯函数进行二维高斯滤波，得到各像素点的Y型的中心响应值；预设圆形的的亚单元区域，对于各像素点，采用X型中心二维高斯函数分别对以该像素点为中心的亚单元区域内的各像素点的Y型的中心响应值进行二维高斯滤波，并将得到的全部响应值求和取均值，得到各像素点的Y型亚单元总响应值；对于各像素点，采用Y型周边二维高斯函数对其Y型亚单元总响应值进行二维高斯滤波，得到各像素点的Y型亚单元抑制响应值；对于各像素点，将其Y型亚单元总响应值减去其Y型亚单元抑制响应值与Y型中心周边拮抗强度的乘积，得到各像素点的Y型初始响应值；
[0012] D、X、Y型细胞的经典感受野响应值的计算：
[0013] 预设分别与X型细胞、Y型细胞相对应的包含多个方向参数的二维高斯一阶偏导函数，对于各像素点，采用二维高斯一阶偏导函数分别对其X型初始响应值、Y型初始响应值进行滤波，得到各像素点各方向参数的X型经典感受野初始响应值、Y型经典感受野初始响应值；对于各像素点，分别在其各方向参数的X型经典感受野初始响应值及Y型经典感受野初始响应值中取最大值，分别作为该像素点的X型经典感受野响应值、Y型经典感受野响应值；
[0014] E、X、Y型细胞的轮廓响应值的计算：
[0015] 预设分别与X型细胞、Y型细胞相对应的归一化高斯差分函数、非经典感受野拮抗强度，对于各像素点，采用归一化高斯差分函数分别对其X型经典感受野响应值、Y型经典感受野响应值进行滤波，得到各像素点的X型非经典感受野响应值、Y型非经典感受野响应值；对于各像素点，将其X型经典感受野响应值、Y型经典感受野响应值分别减去X型非经典感受野响应值与非经典感受野拮抗强度的乘积、Y型非经典感受野响应值与非经典感受野拮抗强度的乘积，分别得到X型轮廓响应值、Y型轮廓响应值；
[0016] F、最终轮廓响应值的计算：
[0017] 预设局部区域、局部亮度对比度函数、亮度对比度阈值，所述的局部区域为方形；对于各像素点，采用局部亮度对比度函数对以其为中心的局部区域内的各像素点的灰度值进行计算，得到各像素点的局部亮度对比度；对于各像素点，对其局部亮度对比度进行归一化处理，得到各像素点的局部亮度对比度归一值；对于各像素点，若其局部亮度对比度归一值大于或等于亮度对比度阈值，则将其X型轮廓响应值作为该像素点的最终轮廓响应值，若其局部亮度对比度归一值小于亮度对比度阈值，则将其Y型轮廓响应值作为该像素点的最终轮廓响应值；
[0018] G、最终轮廓响应值的计算：
[0019] 对于各像素点，对其最终轮廓响应值进行非极大值抑制和双阈值处理，得到各像素点的最终轮廓值。
[0020] 优选地，所述的步骤B具体如下：
[0021] 所述的X型初始响应值RX(x,y)＝I(x,y)*Gxc(x,y)-αx·I(x,y)*Gxf(x,y) (1)；
[0022] 其中I(x,y)为各像素点的灰度值，*为卷积运算；
[0023] X型中心二维高斯函数
[0024] X型周边二维高斯函数
[0025] σx为X型细胞的空间标准差；αx为X型中心周边拮抗强度。
[0026] 优选地，所述的步骤C如下：
[0027] 所述的Y型的中心响应值CY(x,y)＝I(x,y)*Gyc(x,y) (2)；
[0028] 其中Y型中心二维高斯函数
[0029] Y型细胞的空间标准差
[0030] 所述的Y型亚单元总响应值
[0031] 其中亚单元区域S1的半径为rs＝15σx， 表示亚单元区域S1内的各像素点的Y型的中心响应值CY(x,y)，(x',y')表示亚单元区域S1内的各像素点的坐标，N表示亚单元区域S1内的像素点数目；
[0032] Y型初始响应值RY(x,y)＝NY(x,y)-αy·NY(x,y)*Gyf(x,y) (4)；
[0033] 其中Y型周边二维高斯函数
[0034] αy为Y型中心周边拮抗强度。
[0035] 优选地，所述的步骤D如下：
[0036] X型经典感受野初始响应值CRFX(x,y,θi)＝|RX(x,y)*VX(x,y,θi)| (5)；
[0037] Y型经典感受野初始响应值CRFY(x,y,θi)＝|RY(x,y)*VY(x,y,θi)| (6)；
[0038] 其中对应于X型细胞的二维高斯一阶偏导函数
[0039] 对应于Y型细胞的二维高斯一阶偏导函数
[0040] Nθ为方向参数的数目；γ为表示感受野椭圆度的常数；
[0041] X型经典感受野响应值CRFX(x,y)＝max(CRFX(x,y,θi)|i＝1,2,…Nθ) (7)；
[0042] Y型经典感受野响应值CRFY(x,y)＝max(CRFY(x,y,θi)|i＝1,2,…Nθ) (8)。
[0043] 优选地，所述的步骤E如下：
[0044] X型轮廓响应值ResX(x,y)＝CRFX(x,y)-β·CRFX(x,y)*wX(x,y) (9)；
[0045] Y型轮廓响应值ResY(x,y)＝CRFY(x,y)-β·CRFY(x,y)*wY(x,y) (10)；
[0046] 对应X型细胞的归一化高斯差分函数
[0047] 对应Y型细胞的归一化高斯差分函数
[0048] 其中，
[0049]
[0050] 其中||·||1为L1范数，H(x)＝max(0,x)，β为非经典感受野拮抗强度。
[0051] 优选地，所述的步骤F如下：
[0052] 局部亮度对比度
[0053] 其中
[0054]
[0055] 其中s2表示以当前像素点为中心的11*11像素的方形局部区域，f(xj,yj)表示区域s2内的各像素点的灰度值，j表示区域s2内的相对坐标，j∈[-5,5]；
[0056] 归一化局部亮度对比度
[0057] 其中ΔCmax(x,y)为各个像素点的局部亮度对比度ΔC(x,y)中的最大值，ΔCmin(x,y)为各个像素点的局部亮度对比度ΔC(x,y)中的最小值；
[0058] 最终轮廓响应值
[0059] 本发明对分别对X型细胞及Y型细胞进行的感受响应进行模拟，通过中心及周边响应的差异计算X型细胞的初始响应；同时还基于Y型细胞感受野比X型细胞大的特点，对于Y型细胞选用更大的感受野空间标准差，并采用亚单元区域对Y型细胞进行计算，以模拟Y型细胞的整流特性；而后通过DOG模板对X、Y型细胞的初始响应进行滤波，以模拟非经典感受野的响应；最后通过引入局部亮度对比度函数，对X、Y型细胞的双通道进行整合，通过阈值判断选择X型或Y型细胞的轮廓响应，实现X、Y型细胞的双通道的合理融合，提高轮廓检测的成功率。
附图说明
[0060] 图1为本发明的轮廓检测方法流程图；
[0061] 图2为实施例1方法与文献1轮廓检测模型的检测效果对比图。
具体实施方式
[0062] 下面结合实施例具体说明本发明。
[0063] 实施例1
[0064] 如图1所示，本实施例提供的基于X、Y并行视觉通道响应的仿生型轮廓检测方法，包括以下步骤：
[0065] A、输入经灰度处理的待检测图像，构建视网膜神经节细胞与外膝体中的X型细胞和Y型细胞仿真模型；
[0066] B、X型细胞的初始响应值的计算：
[0067] 预设X型中心二维高斯函数、X型周边二维高斯函数、X型中心周边拮抗强度，X型周边二维高斯函数的空间标准差大于X型中心二维高斯函数，对于各像素点，分别采用X型中心二维高斯函数、X型周边二维高斯函数对其灰度值进行二维高斯滤波，分别得到各像素点的X型中心响应值和X型周边响应值，将各像素点的X型中心响应值减去其X型周边响应值与X型中心周边拮抗强度的乘积，得到各像素点的X型初始响应值；
[0068] 所述的步骤C具体如下：
[0069] 所述的X型初始响应值RX(x,y)＝I(x,y)*Gxc(x,y)-αx·I(x,y)*Gxf(x,y) (1)；
[0070] 其中I(x,y)为各像素点的灰度值，*为卷积运算；
[0071] X型中心二维高斯函数
[0072] X型周边二维高斯函数
[0073] σx为X型细胞的空间标准差；αx为X型中心周边拮抗强度；
[0074] C、Y型细胞的初始响应值的计算：
[0075] 预设Y型中心二维高斯函数、Y型周边二维高斯函数、Y型中心周边拮抗强度，Y型中心二维高斯函数的空间标准差大于X型中心二维高斯函数，Y型周边二维高斯函数的空间标准差大于Y型中心二维高斯函数；对于各像素点，分别采用Y型中心二维高斯函数进行二维高斯滤波，得到各像素点的Y型的中心响应值；预设圆形的的亚单元区域，对于各像素点，采用X型中心二维高斯函数分别对以该像素点为中心的亚单元区域内的各像素点的Y型的中心响应值进行二维高斯滤波，并将得到的全部响应值求和取均值，得到各像素点的Y型亚单元总响应值；对于各像素点，采用Y型周边二维高斯函数对其Y型亚单元总响应值进行二维高斯滤波，得到各像素点的Y型亚单元抑制响应值；对于各像素点，将其Y型亚单元总响应值减去其Y型亚单元抑制响应值与Y型中心周边拮抗强度的乘积，得到各像素点的Y型初始响应值；
[0076] 所述的步骤C具体如下：
[0077] 所述的Y型的中心响应值CY(x,y)＝I(x,y)*Gyc(x,y) (2)；
[0078] 其中Y型中心二维高斯函数
[0079] Y型细胞的空间标准差
[0080] 所述的Y型亚单元总响应值
[0081] 其中亚单元区域S1的半径为rs＝15σx， 表示亚单元区域S1内的各像素点的Y型的中心响应值CY(x,y)，(x',y')表示亚单元区域S1内的各像素点的坐标，N表示亚单元区域S1内的像素点数目；
[0082] Y型初始响应值RY(x,y)＝NY(x,y)-αy·NY(x,y)*Gyf(x,y) (4)；
[0083] 其中Y型周边二维高斯函数
[0084] αy为Y型中心周边拮抗强度；
[0085] D、X、Y型细胞的经典感受野响应值的计算：
[0086] 预设分别与X型细胞、Y型细胞相对应的包含多个方向参数的二维高斯一阶偏导函数，对于各像素点，采用二维高斯一阶偏导函数分别对其X型初始响应值、Y型初始响应值进行滤波，得到各像素点各方向参数的X型经典感受野初始响应值、Y型经典感受野初始响应值；对于各像素点，分别在其各方向参数的X型经典感受野初始响应值及Y型经典感受野初始响应值中取最大值，分别作为该像素点的X型经典感受野响应值、Y型经典感受野响应值；
[0087] 所示的步骤D具体如下：
[0088] X型经典感受野初始响应值CRFX(x,y,θi)＝|RX(x,y)*VX(x,y,θi)| (5)；
[0089] Y型经典感受野初始响应值CRFY(x,y,θi)＝|RY(x,y)*VY(x,y,θi)| (6)；
[0090] 其中对应X型细胞的二维高斯一阶偏导函数
[0091] 对应Y型细胞的二维高斯一阶偏导函数
[0092] Nθ为方向参数的数目；γ为表示感受野椭圆度的常数；
[0093] X型经典感受野响应值CRFX(x,y)＝max(CRFX(x,y,θi)|i＝1,2,…Nθ) (7)；
[0094] Y型经典感受野响应值CRFY(x,y)＝max(CRFY(x,y,θi)|i＝1,2,…Nθ) (8)；
[0095] E、X、Y型细胞的轮廓响应值的计算：
[0096] 预设分别与X型细胞、Y型细胞相对应的归一化高斯差分函数、非经典感受野拮抗强度，对于各像素点，采用归一化高斯差分函数分别对其X型经典感受野响应值、Y型经典感受野响应值进行滤波，得到各像素点的X型非经典感受野响应值、Y型非经典感受野响应值；对于各像素点，将其X型经典感受野响应值、Y型经典感受野响应值分别减去X型非经典感受野响应值与非经典感受野拮抗强度的乘积、Y型非经典感受野响应值与非经典感受野拮抗强度的乘积，分别得到X型轮廓响应值、Y型轮廓响应值；
[0097] 所述的步骤E具体如下：
[0098] X型轮廓响应值ResX(x,y)＝CRFX(x,y)-β·CRFX(x,y)*wX(x,y) (9)；
[0099] Y型轮廓响应值ResY(x,y)＝CRFY(x,y)-β·CRFY(x,y)*wY(x,y) (10)；
[0100] 对应X型细胞的归一化高斯差分函数
[0101] 对应Y型细胞的归一化高斯差分函数
[0102] 其中，
[0103]
[0104] 其中||·||1为L1范数，H(x)＝max(0,x)，β为非经典感受野拮抗强度；
[0105] F、最终轮廓响应值的计算：
[0106] 预设局部区域、局部亮度对比度函数、亮度对比度阈值，所述的局部区域为方形；对于各像素点，采用局部亮度对比度函数对以其为中心的局部区域内的各像素点的灰度值进行计算，得到各像素点的局部亮度对比度；对于各像素点，对其局部亮度对比度进行归一化处理，得到各像素点的局部亮度对比度归一值；对于各像素点，若其局部亮度对比度归一值大于或等于亮度对比度阈值，则将其X型轮廓响应值作为该像素点的最终轮廓响应值，若其局部亮度对比度归一值小于亮度对比度阈值，则将其Y型轮廓响应值作为该像素点的最终轮廓响应值；
[0107] 所示的步骤F具体如下：
[0108] 局部亮度对比度函数
[0109] 其中
[0110]
[0111] 其中s2表示以当前像素点为中心的11*11像素的方形局部区域，f(xj,yj)表示区域s2内的各像素点的灰度值，j表示区域s2内的相对坐标，j∈[-5,5]；
[0112] 归一化局部亮度对比度
[0113] 其中ΔCmax(x,y)为各个像素点的局部亮度对比度ΔC(x,y)中的最大值，ΔCmin(x,y)为各个像素点的局部亮度对比度ΔC(x,y)中的最小值；
[0114] 最终轮廓响应值
[0115] G、最终轮廓响应值的计算：
[0116] 对于各像素点，对其最终轮廓响应值进行非极大值抑制和双阈值处理，得到各像素点的最终轮廓值。
[0117] 下面将本实施例的轮廓检测方法与文献1提供的轮廓检测方法进行有效性对比，文献1如下：
[0118] 文献1：Yang K F,Gao S B,Guo C F,et al.Boundary detection using double-opponency and spatial sparseness constraint[J].IEEE Transactions on Image Processing,2015,24(8):2565-2578.
[0119] 为保证对比的有效性，对于本实施例采用与文献1中相同的非极大值抑制和双阈值处理进行最后的轮廓整合，其中包含的两个阈值th,tl，设置为tl＝0.5th，由阈值分位数p计算而得；
[0120] 其中性能评价指标F采用文献1中给出的如下标准：
[0121]
[0122] 其中P表示精确率，R表示召回率，性能评测指标F的取值在[0,1]之间，越接近1表示轮廓检测的效果越好，另外，定义容忍度为：在5*5的邻域内检测到的都算正确检测。
[0123] 选取伯克利分割数据集(BSDS300)随机的三幅自然图像，分别采用实施例1方案与文献1方案进行检测，对应的真实轮廓图、文献1方法检测的最优轮廓图，实施例1方法检测的最优轮廓图如图2所示；其中，文献1方法检测的最优轮廓图，实施例1方法检测的最优轮廓图中的右上角的数字为其对应的性能评价指标F的数值，表1为实施例1选用的参数值；
[0124] 表1实施例1参数组表
[0125]
[0126] 由图2可以看出，实施例1方案的轮廓检测结果优于文献1方案。
