首页 > 专利 > 广西科技大学 > 模拟人眼微动的图像轮廓检测方法专利详情

模拟人眼微动的图像轮廓检测方法 0 0

有效专利查看PDF

申请进展

基本信息

申请人信息

代理人信息

摘要

法律状态

权利要求

说明书

专利申请流程有哪些步骤？

申请

申请号：指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。

申请日：提出专利申请之日。

2017-04-11

申请公布

申请公布指发明专利申请经初步审查合格后，自申请日（或优先权日）起18个月期满时的公布或根据申请人的请求提前进行的公布。

申请公布号：专利申请过程中，在尚未取得专利授权之前，国家专利局《专利公报》公开专利时的编号。

申请公布日：申请公开的日期，即在专利公报上予以公开的日期。

2017-09-12

授权

授权指对发明专利申请经实质审查没有发现驳回理由，授予发明专利权；或对实用新型或外观设计专利申请经初步审查没有发现驳回理由，授予实用新型专利权或外观设计专利权。

2020-01-31

预估到期

发明专利权的期限为二十年，实用新型专利权期限为十年，外观设计专利权期限为十五年，均自申请日起计算。专利届满后法律终止保护。

2037-04-11

基本信息

有效性	有效专利	专利类型	发明专利
申请号	CN201710230521.X	申请日	2017-04-11
公开/公告号	CN107067408B	公开/公告日	2020-01-31
授权日	2020-01-31	预估到期日	2037-04-11
申请年	2017年	公开/公告年	2020年
缴费截止日
分类号	G06T7/13 、G06K9/46	主分类号	G06T7/13
是否联合申请	独立申请	文献类型号	B
独权数量	1	从权数量	5
权利要求数量	6	非专利引证数量	1
引用专利数量	8	被引证专利数量	0
非专利引证	1、Chuan Lin 等.Improved contourdetection model with spatial summationproperties based on nonclassicalreceptive field《.Journal of ElectronicImaging》.2016,第25卷(第4期),Cosmin Grigorescu 等.Improved ContourDetection by Non-classical ReceptiveField Inhibition《.Biologically MotivatedComputer Vision》.2002,第2525卷林川等.考虑微动机制与感受野特性的轮廓检测模型《.计算机工程与应用》.2016,第52卷(第24期),李康群等.基于视通路多感受野朝向性关联的轮廓检测方法《.中国生物医学工程学报》.2017,第36卷(第1期),;
引用专利	CN106033606A、CN106338733A、CN104484667A、CN106033608A、CN106156779A、CN106127740A、US2015062000A1、JP4579191B2	被引证专利
专利权维持	5	专利申请国编码	CN
专利事件		事务标签	公开、实质审查、授权

申请人信息

申请人	广西科技大学	第一申请人	广西科技大学
专利权人	广西科技大学	当前专利权人	广西科技大学
发明人	林川、曹以隽	第一发明人	林川
地址	广西壮族自治区柳州市东环路268号	邮编	545006
申请人数量	1	发明人数量	2
申请人所在省	广西壮族自治区	申请人所在市	广西壮族自治区柳州市

代理人信息

代理机构

专利代理机构是经省专利管理局审核，国家知识产权局批准设立，可以接受委托人的委托，在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。

长沙正奇专利事务所有限责任公司

代理人

专利代理师是代理他人进行专利申请和办理其他专利事务，取得一定资格的人。

周晟

摘要

本发明提供一种模拟人眼微动的图像轮廓检测方法，包括以下步骤：A、输入待检测图像，预设整体抑制参数、抑制系数及Gabor滤波器组，滤波得到各像素点经典感受野刺激响应；B、对经典感受野刺激响应进行截断处理；C、利用DoG模板与临时中心区，计算得到各像素点的归一化的权重函数；计算各像素点各方向下的非经典感受野刺激响应粗算值，并求标准差；D、计算得到各像素点的标准差权重，进而计算得到各像素点的非经典感受野刺激响应终算值；E、计算得到各像素点的轮廓识别值，构成待检测图像的轮廓识别图像。该方法克服现有技术的缺陷，在抑制强纹理的同时保留弱的边缘，提高轮廓识别的成功率。

摘要附图
说明书附图：图1
说明书附图：图2
说明书附图：图3

法律状态

序号	法律状态公告日	法律状态	法律状态信息
1	2020-01-31	授权
2	2017-09-12	实质审查的生效	IPC(主分类): G06T 7/13 专利申请号: 201710230521.X 申请日: 2017.04.11
3	2017-08-18	公开

权利要求

权利要求书是申请文件最核心的部分，是申请人向国家申请保护他的发明创造及划定保护范围的文件。

1.模拟人眼微动的图像轮廓检测方法，其特征在于包括以下步骤：
A、输入经灰度处理的待检测图像，预设整体抑制参数与抑制系数，预设沿圆周均匀分布的多个方向参数的Gabor滤波器组，对待检测图像中的各像素点分别按照各方向参数进行Gabor滤波，获得各像素点的各方向的Gabor能量值；对于各像素点，选取其各方向的Gabor能量值中的最大值，作为该像素点经典感受野刺激响应；
B、对于各像素点，将其经典感受野刺激响应进行截断处理，得到各像素点截断后的经典感受野刺激响应；
C、利用高斯差分函数DoG模板，构建一组临时中心区，各个临时中心区相对于视野中心区具有不同偏移角度；对于各像素点，将其临时中心区响应与DoG模板进行整合与归一化，得到一组归一化的权重函数；
对于各像素点，在不同偏移角度下，将归一化的权重函数与DoG模板内的截断后的经典感受野刺激响应作乘积后求和，得到各像素点在各偏移角度下的非经典感受野刺激响应粗算值；对各像素点在各偏移角度下的非经典感受野刺激响应粗算值求标准差；
D、对于各像素点，结合各偏移角度下的非经典感受野刺激响应粗算值的标准差及整体抑制参数计算得到标准差权重；将标准差权重与各偏移角度下的非经典感受野刺激响应粗算值的最小值进行乘积得到该像素点的非经典感受野刺激响应终算值；
E、对于各像素点，将其经典感受野刺激响应与非经典感受野刺激响应终算值结合抑制系数计算得到该像素点的综合刺激响应，即为该像素点的轮廓识别值，将待检测图像全部像素点的轮廓识别值进行非极大值抑制和值二值化之后即得到待检测图像的轮廓识别图像。

2.如权利要求1所述的模拟人眼微动的图像轮廓检测方法，其特征在于：
所述的步骤A中经典感受野刺激响应的计算具体如下：
所述的Gabor滤波器组的二维Gabor函数表达式如下：
其中 γ为一个表示椭圆形感受野长短轴比例的常
数，参数λ为波长，σ为Gabor函数的标准差以及DoG模板中心区的带宽，1/λ为余弦函数的空间频率，是相角参数，θ为Gabor滤波的角度参数；
I(x,y)为待检测图像，*为卷积运算符；
Gabor能量值计算如下：
其中θi为Gabor滤波的某一角度，Nθ为Gabor滤波的角度的个数；
E(x,y；σ)为像素点(x,y)的各角度Gabor滤波能量值的最大值，即为像素点(x,y)的经典感受野刺激响应。

3.如权利要求2所述的模拟人眼微动的图像轮廓检测方法，其特征在于：
所述的步骤B中截断后的经典感受野刺激响应的计算过程如下：
利用上限比例PH∈(0,1)和下限比例PL∈(0,1)对E(x,y；σ)进行截断：
将各像素点的E(x,y；σ)从小到大进行选取，选取出PH对应百分比数目的E(x,y；σ)，其中的最大值设为QH，作为上限分位数；
将各像素点的E(x,y；σ)从小到大进行选取，选取出PL对应百分比数目的E(x,y；σ)，其中的最大值设为QL，作为下限分位数；
截断后的经典感受野刺激响应：

4.如权利要求3所述的模拟人眼微动的图像轮廓检测方法，其特征在于：
所述的步骤C中的DoG模板的表达式：
其中k为控制DoG模板大小的参数；
所述的临时中心区响应的表达式如下：
其中d表示视野中心到临时中心区的距离，表示临时中心区的偏移角度；
所述的各像素点的整合与归一化过程如下：
通过归一化的权重函数表达式进行，表达式如下：
其中w(x,y；σ,φ)＝wm(x,y；d,φ)·DoG(x,y；σ,k)，||·||1为(L1)范数正则化，H(X)为取正值的函数；
所述的各像素点各偏移角度下的非经典感受野刺激响应粗算值的计算过程如下：
其中：Inhe(x,y；σ,φj)为各像素点各偏移角度下的非经典感受野刺激响应粗算值；
-3kσφj表示多个偏移角度；
所述的各像素点各偏移角度下的非经典感受野刺激响应粗算值的平均值及标准差的计算过程如下：
其中STDinh(x,y)为各像素点各方向下的非经典感受野刺激响应粗算值的标准差，Aveinh(x,y)为各像素点各方向下的非经典感受野刺激响应粗算值的平均值。

5.如权利要求4所述的模拟人眼微动的图像轮廓检测方法，其特征在于：
所述的步骤D中的标准差权重的计算过程如下：
其中wstd(x,y；σ)为标准差权重，fos为整体抑制参数；
非经典感受野刺激响应终算值的计算过程如下：
Inh(x,y；σ)＝Inhm(x,y；σ)·wstd(x,y；σ)   (15)；
Inhm(x,y；σ)＝min{Inhe(x,y；σ,φj)|j＝1,2,...,Nφ}   (16)；
其中Inhm(x,y；σ)为Inhe(x,y；σ,φj)的最小值。

6.如权利要求5所述的模拟人眼微动的图像轮廓检测方法，其特征在于：
所述的步骤E中的综合刺激响应的计算过程如下：
R(x,y)＝H(E(x,y；σ)-αInh(x,y；σ))  (17)；
其中R(x,y)为像素点的综合刺激响应，α为抑制系数。

说明书

技术领域

[0001] 本发明涉及图像处理领域，具体涉及一种模拟人眼微动的图像轮廓检测方法。

背景技术

[0002] 早期的生理学和神经科学在视觉信息处理的研究大都是基于中心有限区域的经典感受野。然而，后期更多的研究表明关在更大范围内的光刺激能调制经典感受野响应，这个外周区域被叫做非经典感受野。这个调制功能使得神经元能够整合更大范围内的信息并传递给随后的视觉过程。相比较于神经节细胞和外膝体细胞的抑制，初级视觉皮层的非经典感受野有着更加复杂的特性。有文献表示非经典感受野可大概分为四种模式：(1)全范围抑制；(2)全范围促进；(3)两边抑制，两端促进；(4)两边促进，两端抑制。另外，非经典感受野有着独立的方向选择性。当经典感受野和非经典感受野接受不同类型的方向，亮度，空间频率，空间相位，和运动速度刺激时，就会产生更强的响应，当经典感受野和非经典感受野接受相同类型的方向，亮度，空间频率，空间相位，和运动速度刺激时，就会产生更弱的响应。

[0003] 然而，这些感受野的性质大多都是基于麻醉动物实验。在这种情况下，人眼的运动就被忽视了。但实际上人眼的运动对视觉信息处理的脑机制很有帮助。具体的，人眼运动能分为稳定注视，运动注视和固定的眼动。其中固定的眼动又包括震颤，漂移和微跳动。眼动影响着视网膜图像和其后续的视觉系统，包括外膝体，初级视觉皮层和高级视觉皮层，因此在进行轮廓检测时，人眼微动是需要重点考虑的因素，急需一种能够人眼微动的轮廓检测方法。

发明内容

[0004] 本发明旨在提供一种模拟人眼微动的图像轮廓检测方法，该方法克服现有技术未考虑人眼微动机制的缺陷，在很好地抑制强纹理的同时保留弱的边缘，提高轮廓识别的成功率。

[0005] 本发明的技术方案如下：模拟人眼微动的图像轮廓检测方法，包括以下步骤：

[0006] A、输入经灰度处理的待检测图像，预设整体抑制参数与抑制系数，预设沿圆周均匀分布的多个方向参数的Gabor滤波器组，对待检测图像中的各像素点分别按照各方向参数进行Gabor滤波，获得各像素点的各方向的Gabor能量值；对于各像素点，选取其各方向的Gabor能量值中的最大值，作为该像素点经典感受野刺激响应；

[0007] B、对于各像素点，将其经典感受野刺激响应进行截断处理，得到各像素点截断后的经典感受野刺激响应；

[0008] C、利用高斯差分函数DoG模板，构建一组临时中心区，各个临时中心区相对于视野中心区具有不同偏移角度；对于各像素点，将其临时中心区响应与DoG模板进行整合与归一化，得到一组归一化的权重函数；

[0009] 对于各像素点，在不同偏移角度下，将归一化的权重函数与DoG模板内的截断后的经典感受野刺激响应作乘积后求和，得到各像素点在各偏移角度下的非经典感受野刺激响应粗算值；对各像素点在各偏移角度下的非经典感受野刺激响应粗算值求标准差；

[0010] D、对于各像素点，结合各偏移角度下的非经典感受野刺激响应粗算值的标准差及整体抑制参数计算得到标准差权重；将标准差权重与各偏移角度下的非经典感受野刺激响应粗算值的最小值进行乘积得到该像素点的非经典感受野刺激响应终算值；

[0011] E、对于各像素点，将其经典感受野刺激响应与非经典感受野刺激响应终算值结合抑制系数计算得到该像素点的综合刺激响应，即为该像素点的轮廓识别值，将待检测图像全部像素点的轮廓识别值进行非极大值抑制和值二值化之后即得到待检测图像的轮廓识别图像。

[0012] 优选地，所述的步骤A中经典感受野刺激响应的计算具体如下：

[0013] 所述的Gabor滤波器组的二维Gabor函数表达式如下：

[0014]

[0015] 其中 γ为一个表示椭圆形感受野长短轴比例的常数，参数λ为波长，σ为Gabor函数的标准差以及DoG模板中心区的带宽，1/λ为余弦函数的空间频率，是相角参数，θ为Gabor滤波的角度参数；

[0016]

[0017] I(x,y)为待检测图像，*为卷积运算符；

[0018] Gabor能量值计算如下：

[0019]

[0020]

[0021] 其中θi为Gabor滤波的某一角度，Nθ为Gabor滤波的角度的个数；

[0022]

[0023] E(x,y；σ)E(x,y)为像素点(x,y)的各角度Gabor滤波能量值的最大值，即为像素点(x,y)的经典感受野刺激响应。

[0024] 优选地，所述的步骤B中截断后的经典感受野刺激响应的计算过程如下：

[0025] 利用上限比例PH∈(0,1)和下限比例PL∈(0,1)对E(x,y；σ)进行截断：

[0026] 将各像素点的E(x,y；σ)从小到大进行选取，选取出PH对应百分比数目的E(x,y；σ)，其中的最大值设为QH，作为上限分位数；

[0027] 将各像素点的E(x,y；σ)从小到大进行选取，选取出PL对应百分比数目的E(x,y；σ)，其中的最大值设为QL，作为下限分位数；

[0028] 截断后的经典感受野刺激响应：

[0029]

[0030] 优选地，所述的步骤C中的DoG模板的表达式：

[0031]

[0032] 其中k为控制DoG模板大小的参数；

[0033] 所述的临时中心区响应的表达式如下：

[0034]

[0035] 其中d表示视野中心到临时中心区的距离，表示临时中心区的偏移角度；

[0036] 所述的各像素点的整合与归一化过程如下：

[0037] 通过归一化的权重函数表达式进行，表达式如下：

[0038]

[0039] 其中w(x,y；d,φ)＝wm(x,y；d,φ)·DoG(x,y；σ,k)，||·||1为(L1)范数正则化，H(X)为取正值的函数；

[0040] 所述的各像素点各偏移角度下的非经典感受野刺激响应粗算值的计算过程如下：

[0041]

[0042] 其中：Inhe(x,y；σ,φj)为各像素点各偏移角度下的非经典感受野刺激响应粗算值；

[0043] -3kσ

[0044] φj表示多个偏移角度；

[0045] 所述的各像素点各偏移角度下的非经典感受野刺激响应粗算值的平均值及标准差的计算过程如下：

[0046]

[0047]

[0048] 其中STDinh(x,y)为各像素点各方向下的非经典感受野刺激响应粗算值的标准差，Aveinh(x,y)为各像素点各方向下的非经典感受野刺激响应粗算值的平均值。

[0049] 优选地，所述的步骤D中的标准差权重的计算过程如下：

[0050]

[0051]

[0052] 其中wstd(x,y；σ)为标准差权重，fos为整体抑制参数；

[0053] 非经典感受野刺激响应终算值的计算过程如下：

[0054] Inh(x,y；σ)＝Inhm(x,y；σ)·wstd(x,y；σ) (15)；

[0055] Inhm(x,y；σ)＝min{Inhe(x,y；σ,φj)|j＝1,2,...,Nφ} (16)；

[0056] 其中Inhm(x,y；σ)为Inhe(x,y；σ,φj)的最小值。

[0057] 优选地，所述的步骤E中的综合刺激响应的计算过程如下：

[0058] R(x,y)＝H(E(x,y；σ)-αInh(x,y；σ)) (17)；

[0059] 其中R(x,y)为像素点的综合刺激响应，α为抑制系数。

[0060] 本发明方法中设置临时的中心区域对人眼微动产生的偏移进行模拟，假设固定的眼动偏移会导致一个临时的中心区域，此处我们仅仅假设抑制作用仅仅发生于短距离或长距离的神经元连接，即临时的中心区并不会抑制临时的中心区域的响应，通过模拟人眼微动的临时的中心区域保证模拟的真实性以及轮廓检测的准确性；

[0061] 并且，通过多通道滤波器去模拟人眼微动在周边抑制中的的影响，多通道特征的选择是为了模拟人眼微动的的无方向性，提高模拟的真实性以及轮廓检测的准确性；同时，根据实验可以发现同性质的纹理会使得标准差变小，因此，在算法选择时设置当标准差小的时候，抑制程度就会增强；采用标准差的方法对不同的滤波结果融合，可以抑制显著的纹理，将显著的纹理去除掉，以降低误检测的概率；

[0062] 再有，由于过大的Gabor能量值会导致不精确的抑制响应，可能一些不较弱的边缘会被其周边较强的边缘的响应抑制，因此采用能量截断的方法增强较弱的边缘，以降低漏检测的概率，提高检测质量。

实施方案

[0066] 下面结合附图和实施例具体说明本发明。

[0067] 实施例1

[0068] 如图1所示，本实施例提供的模拟人眼微动的图像轮廓检测方法包括以下步骤：

[0069] A、输入经灰度处理的待检测图像，预设整体抑制参数与抑制系数，预设沿圆周均匀分布的多个方向参数的Gabor滤波器组，对待检测图像中的各像素点分别按照各方向参数进行Gabor滤波，获得各像素点的各方向的Gabor能量值；对于各像素点，选取其各方向的Gabor能量值中的最大值，作为该像素点经典感受野刺激响应；

[0070] 所述的步骤A中经典感受野刺激响应的计算具体如下：

[0071] 所述的Gabor滤波器组的二维Gabor函数表达式如下：

[0072]

[0073] 其中 γ为一个表示椭圆形感受野长短轴比例的常数，参数λ为波长，σ为Gabor函数的标准差以及DoG模板中心区的带宽，1/λ为余弦函数的空间频率，σ/λ为空间频率的带宽，是相角参数，θ为Gabor滤波的角度参数；

[0074]

[0075] I(x,y)为待检测图像，*为卷积运算符；

[0076] Gabor能量值计算如下：

[0077]

[0078]

[0079] 其中θi为Gabor滤波的某一角度，Nθ为Gabor滤波的角度的个数；

[0080]

[0081] E(x,y；σ)E(x,y)为像素点(x,y)的各角度Gabor滤波能量值的最大值，即为像素点(x,y)的经典感受野刺激响应；

[0082] B、对于各像素点，将其经典感受野刺激响应进行截断处理，得到各像素点截断后的经典感受野刺激响应；

[0083] 所述的步骤B中截断后的经典感受野刺激响应的计算过程如下：

[0084] 利用上限比例PH∈(0,1)和下限比例PL∈(0,1)对E(x,y；σ)进行截断：

[0085] 将各像素点的E(x,y；σ)从小到大进行选取，选取出PH对应百分比数目的E(x,y；σ)，其中的最大值设为QH，作为上限分位数；本实施例PH＝0.8，即将各像素点的E(x,y；σ)从小到大选取80％数目的值，将其中的最大值设为QH；

[0086] 将各像素点的E(x,y；σ)从小到大进行选取，选取出PL对应百分比数目的E(x,y；σ)，其中的最大值设为QL；本实施例PL＝0.1，即将各像素点的E(x,y；σ)从小到大选取10％数目的值，将其中的最大值设为QL；

[0087] 截断后的经典感受野刺激响应：

[0088]

[0089] C、利用高斯差分函数DoG模板，构建一组临时中心区，各个临时中心区相对于视野中心区具有不同偏移角度；对于各像素点，将其临时中心区响应与DoG模板进行整合与归一化，得到一组归一化的权重函数；

[0090] 对于各像素点，在不同偏移角度下，将归一化的权重函数与DoG模板内的截断后的经典感受野刺激响应作乘积后求和，得到各像素点在各偏移角度下的非经典感受野刺激响应粗算值；对各像素点在各偏移角度下的非经典感受野刺激响应粗算值求标准差；

[0091] 所述的步骤C中的DoG模板的表达式：

[0092]

[0093] 其中k为控制DoG模板大小的参数；

[0094] 所述的临时中心区响应的表达式如下：

[0095]

[0096] 其中d表示视野中心到临时中心区的距离，表示临时中心区的偏移角度；

[0097] 所述的各像素点的整合与归一化过程如下：

[0098] 通过归一化的权重函数表达式进行，表达式如下：

[0099]

[0100] 其中w(x,y；d,φ)＝wm(x,y；d,φ)·DoG(x,y；σ,k)，||·||1为(L1)范数正则化，H(X)为取正值的函数；

[0101] 所述的各像素点各偏移角度下的非经典感受野刺激响应粗算值的计算过程如下：

[0102]

[0103] 其中：Inhe(x,y；σ,φj)为各像素点各偏移角度下的非经典感受野刺激响应粗算值；

[0104] -3kσ

[0105] φj表示多个偏移角度；

[0106] 所述的各像素点各偏移角度下的非经典感受野刺激响应粗算值的平均值及标准差的计算过程如下：

[0107]

[0108]

[0109] 其中STDinh(x,y)为各像素点各方向下的非经典感受野刺激响应粗算值的标准差，Aveinh(x,y)为各像素点各方向下的非经典感受野刺激响应粗算值的平均值；

[0110] D、对于各像素点，结合各偏移角度下的非经典感受野刺激响应粗算值的标准差及整体抑制参数计算得到标准差权重；将标准差权重与各偏移角度下的非经典感受野刺激响应粗算值的最小值进行乘积得到该像素点的非经典感受野刺激响应终算值；

[0111] 优选地，所述的步骤D中的标准差权重的计算过程如下：

[0112]

[0113]

[0114] 其中wstd(x,y；σ)为标准差权重，fos为整体抑制参数；

[0115] 非经典感受野刺激响应终算值的计算过程如下：

[0116] Inh(x,y；σ)＝Inhm(x,y；σ)·wstd(x,y；σ) (15)；

[0117] Inhm(x,y；σ)＝min{Inhe(x,y；σ,φj)|j＝1,2,...,Nφ} (16)；

[0118] 其中Inhm(x,y；σ)为Inhe(x,y；σ,φj)的最小值；

[0119] E、对于各像素点，将其经典感受野刺激响应与非经典感受野刺激响应终算值结合抑制系数计算得到该像素点的综合刺激响应，即为该像素点的轮廓识别值，将待检测图像全部像素点的轮廓识别进行非极大值抑制和值二值化之后即得到待检测图像的轮廓识别图像；

[0120] 所述的步骤E中的综合刺激响应的计算过程如下：

[0121] R(x,y)＝H(E(x,y；σ)-αInh(x,y；σ)) (17)；

[0122] 其中R(x,y)为像素点的综合刺激响应，α为抑制系数。

[0123] 下面将本实施例的轮廓检测方法与文献1提供的轮廓检测各项同性模型以及各项异性模型进行有效性对比，其中选用文献1中的各向同性模型以及各项异性模型进行有效性对比，文献1如下：

[0124] 文献1：Grigorescu C,Petkov N,Westenberg M.Contour detection based on nonclassical receptive field inhibition[J].IEEE Transactions on Image Processing,2003,12(7):729-739.

[0125] 为保证对比的有效性，对于本实施例采用与文献1中相同的非极大值抑制方法进行后续的轮廓整合，其中包含的两个阈值th,tl设置为tl＝0.5th，由阈值分位数p计算而得；

[0126] 其中性能评价指标P采用文献1中给出的如下标准：

[0127]

[0128]

[0129]

[0130] 式中nTP、nFP、nFN分别表示检测得到的正确轮廓、错误轮廓以及遗漏的轮廓的数目，nGT表示真实轮廓的数目，efn表示错误检测参数，efp表示遗漏检测参数；评测指标P取值在[0,1]之间，越接近1表示轮廓检测的效果越好，另外，定义容忍度为：在5*5的邻域内检测到的都算正确检测。

[0131] 选取毛刷、大象、犀牛3副经典图像进行有效性对比，分别采用文献1中的各项同性模型、各项异性模型以及实施例1方法对上述3幅图进行轮廓检测，其中实施例1方法选用的参数组如表1所示，

[0132] 表1 实施例1参数组表

[0133]

[0134] 文献1中的各项同性模型、各项异性模型采用如下80组参数：α＝{1.0，1.2}，σ＝{1.4，1.6，1.8，2.0，2.2，2.4，2.6，2.8}，p＝{0.5，0.4，0.3，0.2，0.1}；

[0135] 选取各项同性模型、各项异性模型以及实施例1方法中效果最好的一组参数进行对比，轮廓提取效果图对比如图2所示，由图2可看出，从轮廓提取的效果上看，实施例1方法均优于文献1中的各项同性模型、各项异性模型；图3为临时中心区的示意图，其中虚线椭圆部分为临时中心区；

[0136] 其中表2为实施例1结果图对应的部分参数表，其余参数参照表1中的数据；表3、4分别为各项同性模型、各项异性模型结果图对应的参数表，表5为实施例1方法与其他对比模型的识别效果对比表，进一步证明，实施例1方法均优于文献1中的各项同性模型、各项异性模型。

[0137] 表2 实施例1结果图对应的部分参数表

[0138]

[0139] 表3 各项同性模型结果图对应的参数表

[0140]

[0141] 表4 各项异性模型结果图对应的参数表

[0142]

[0143] 表5 实验结果对比图

[0144]

附图说明

[0063] 图1为本发明的轮廓检测方法流程图

[0064] 图2为实施例1方法与文件1轮廓检测模型的检测效果对比图

[0065] 图3为实施例1方法与文件1轮廓检测模型的检测参数对比图

1模拟人眼微动的图像轮廓检测方法