CN108648206B - 一种Robert边缘检测膜计算系统及方法 - Google Patents
一种Robert边缘检测膜计算系统及方法 Download PDFInfo
- Publication number
- CN108648206B CN108648206B CN201810429215.3A CN201810429215A CN108648206B CN 108648206 B CN108648206 B CN 108648206B CN 201810429215 A CN201810429215 A CN 201810429215A CN 108648206 B CN108648206 B CN 108648206B
- Authority
- CN
- China
- Prior art keywords
- value
- var
- enzyme
- membrane
- executed
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000003708 edge detection Methods 0.000 title claims abstract description 23
- 238000000034 method Methods 0.000 title abstract description 16
- 239000012528 membrane Substances 0.000 claims abstract description 47
- 102000004190 Enzymes Human genes 0.000 claims abstract description 43
- 108090000790 Enzymes Proteins 0.000 claims abstract description 43
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000004519 manufacturing process Methods 0.000 claims description 4
- 239000002473 artificial blood Substances 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000004913 activation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052731 fluorine Inorganic materials 0.000 description 1
- 125000001153 fluoro group Chemical group F* 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/181—Segmentation; Edge detection involving edge growing; involving edge linking
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明公开了一种Robert边缘检测膜计算系统及方法,该系统包括3个膜,采用膜计算系统结构,设计了基于膜计算系统的Robert边缘检测执行框架,利用膜计算系统在线性时间内能产生指数增长空间的特性,不管图像分辨率如何增加,所需的计算复杂度为O(1)。其核心优势在于充分利用酶数值膜系统的并行性、酶变量对算法逻辑的可控性、酶数值膜系统以数值变量为处理对象。
Description
技术领域
本发明属于Robert边缘检测技术领域,具体地说,涉及一种Robert边缘检测膜计算系统及方法。
背景技术
在现有PC串行计算模式下,对一幅n×n的图像,用Robert边缘检测算子检测边缘时,计算复杂度随图像分辨率的增加呈平方阶次增加。
Robert边缘检测原理如下:以水平与垂直方向梯度算子计算对角线方向相邻像素之差,以近似代替梯度幅度,通过梯度幅度和设定的阈值大小确定图像边缘。
水平与垂直方向梯度算子分别如下。
对图像I,假定大小为n×n,边缘检测后的图象为edg,传统的Robert边缘检测过程如下:
(1)先利用水平梯度算子计算水平梯度fx;
(2)利用垂直梯度算子计算垂直梯度fy;
(4)设定阈值θ,若g(i,j)>θ,edg(i,j)=1,反之edg(i,j)=0;
算法执行流程如下:
该算法中所有语句的频度之和(即算法的时间耗费)为:T(n)=6n2+n。算法的时间复杂度为:O(n2);随着数据量的增加,计算量以平方阶次增长,当图像分辨率较大时,计算将会非常耗时。本发明即设计一种新的计算模式,时间复杂度和图像分辨率无关,始终为常数阶O(1)。
发明内容
本发明目的在于提供一Robert边缘检测膜计算系统及方法。该系统采用酶数值膜系统结构,设计了针对Robert边缘检测的膜计算系统,给出了相应定义、系统结构图以及系统运行过程。利用膜计算系统在线性时间内能产生指数增长空间的特性,所设计的Robert边缘检测执行方法不管图像分辨率如何增加,所需的计算时间复杂度为常数阶。
其具体技术方案为:
一种Robert边缘检测膜计算系统,度为3,具体为:
П=(m,H,μ,(Var1,E1,Pr1,Var1(0)),(Var2,E2,Pr2,Var2(0)),(Var3,E3,Pr3,Var3(0)))
其中:
(1)m为П的度,代表酶数值膜系统中膜的层次数量,本方法中m=3。
(2)H为膜标号,本系统采用3个膜,标号分别为1,2,3。
(3)μ=[[[]3]2]1代表酶数值膜的结构。
(4)Var1,Var2,Var3代表3个膜内的所有对象变量的集合,其对应的初始值为 Var1(0),Var2(0),Var3(0)。Var1={xij,edgij,ed1,ed2,ed3,CT,EH},其中,xij表示图像I中,坐标为(i,j) 元素的灰度值,对应的初始值分别为:xij=vij表示图像的灰度值。edgij表示边缘检测结束后,坐标为(i,j)元素的灰度值,取值范围为0,1;初始值为edgij=0。ed1=0;ed2=1;ed3=-256。 CT数值变量,初始值CT=1。EH为酶变量,初始值为空,EH=[],表示图像I中是否所有像素位置处的边缘被计算出,若所有值被计算,其值不为空,反之为空。Var2={fxij,fyij,gij}, fxij、fyij、gij分别表示图像I中坐标(i,j)处的水平、垂直与总体梯度值,对应的初始值均为空。diffij表示图像I中坐标(i,j)处总体梯度与设定阈值的差,对应的初始值为空。Eij酶变量,表示对应的diffij是否被计算出,若被计算出,其值不为空,反之为空;酶变量,表示图像I中坐标(i,j)处的边缘是否被计算出,若被计算,其值不为空,反之为空。
(6)Pri,表示膜i内的规则集,规则由两部分构成,一部分是值产生规则,一部分是值分配规则。
举例如下:
Pri=(F(x1,x2,...xN),e,c1|v1+c2|v2+...+cn|vn)
这里,(x1,x2,...xN)是系统对象的变量值,F(x1,x2,...xN)是值产生函数;c1|v1+c2|v2+...+cn|vn是值分配协议。整个规则的含义是:如果酶变量e的值大于(x1,x2,...xN)中的最小值,则Pri这条规则激活。对象变量通过F(x1,x2,...xN)函数运算后产生一个数值d,d按分配协议分配给变量v1,v2,...vn;具体为:
在本发明中,Pr1、Pr2、Pr3具体为:
1)Pr1表示膜1-Main中的规则,共有1条规则,Pr1,main,其含义为:
Pr1,main:EC(EH→)ET:表示当EH的值大于EC时,将EC的值赋给ET。
2)Pr2表示膜2-ComputeGradient中的规则,共有2(n-1)2+n2+2n+2(n-1)条规则,具体含义为:
3)Pr3表示膜3-ComputeEdge中的规则,共有4n2+1条,具体含义为:
一种本发明所述Robert边缘检测膜计算系统的计算方法,包括以下步骤:
所述3个膜中,执行顺序如下:
(3)gij计算出以后,规则执行条件被满足,因此被执行,计算图像上每像素的梯度与设定阈值的绝对差,并将差赋给diffij和酶Eij,由于是1∶1的等比例分配,差diffij和酶Eij的值是相同的。若Eij>0,表示像素(i,j)是边缘,反之不是边缘。
(6)在膜Main中,Pr1,main被执行,ET被赋值为1,整个系统停止运行。
与现有技术相比,本发明的有益效果为:
本发明设计的3个膜,对象数:8×n2+5,其中包含酶对象:2×n2+2,规则数: 4×n2+2×(n-1)2+2n+2(n-1);具体如表1所示。
表1
3个膜,在酶的激活条件下,每条规则同时执行,因此不管图像分辨率有多大,都执行 6步。具体为:
Pr1,main在EH不为0的条件下执行,第6步;
因此,算法执行的时间复杂度为O(1)。
附图说明
图1是本发明膜计算系统结构图;
图2是一个Robert边缘检测膜计算系统实际示例,图2(a-e)给出了系统执行过程中各个变量的值,其中,图2(a)原图,图2(b)fx,图2(c)fy,图2(d)g,图2(e)edg。
具体实施方式
下面结合附图和具体实施例对本发明的技术方案作进一步详细地说明。
一种Robert边缘检测膜系统,度为3,具体为:
П=(m,H,μ,(Var1,E1,Pr1,Var1(0)),(Var2,E2,Pr2,Var2(0)),(Var3,E3,Pr3,Var3(0)))
其中:
(1)m为П的度,代表酶数值膜系统中膜的层次数量,本方法中m=3。
(2)H为膜标号,本系统采用3个膜,标号分别为1,2,3。
(3)μ=[[[]3]2]1代表酶数值膜的结构。
(4)Var1,Var2,Var3代表3个膜内的所有对象变量的集合,其对应的初始值为 Var1(0),Var2(0),Var3(0)。Var1={xij,edgij,ed1,ed2,ed3,CT,EH},其中,xij表示图像I中,坐标为(i,j) 元素的灰度值,对应的初始值分别为:xij=vij表示图像的灰度值。edgij表示边缘检测结束后,坐标为(i,j)元素的灰度值,取值范围为0,1;初始值为edgij=0。ed1=0;ed2=1;ed3=-256。 CT数值变量,初始值CT=1。EH为酶变量,初始值为空,EH=[],表示图像I中是否所有像素位置处的边缘被计算出,若所有值被计算,其值不为空,反之为空。Var2={fxij,fyij,gij}, fxij、fyij、gij分别表示图像I中坐标(i,j)处的水平、垂直与总体梯度值,对应的初始值均为空。diffij表示图像I中坐标(i,j)处总体梯度与设定阈值的差,对应的初始值为空。Eij酶变量,表示对应的diffij是否被计算出,若被计算出,其值不为空,反之为空;酶变量,表示图像I中坐标(i,j)处的边缘是否被计算出,若被计算,其值不为空,反之为空。
(6)Pri,表示膜i内的规则集,规则由两部分构成,一部分是值产生规则,一部分是值分配规则。
举例如下:
Pri=(F(x1,x2,...xN),e,c1|v1+c2|v2+...+cn|vn)
这里,(x1,x2,...xN)是系统对象的变量值,F(x1,x2,...xN)是值产生函数;c1|v1+c2|v2+...+cn|vn是值分配协议。整个规则的含义是:如果酶变量e的值大于(x1,x2,...xN)中的最小值,则Pri这条规则激活。对象变量通过F(x1,x2,...xN)函数运算后产生一个数值d,d按分配协议分配给变量v1,v2,...vn;具体为:
在本发明中,Pr1、Pr2、Pr3具体为:
1)Pr1表示膜Main中的规则,共有1条规则,Pr1,main,其含义为:
Pr1,main:EC(EH→)ET:表示当EH的值大于EC时,将EC的值赋给ET。
2)Pr2表示膜ComputeGradient中的规则,共有2(n-1)2+n2条规则,具体含义为:
3)Pr3表示膜ComputeEdge中的规则,共有4n2+1条,具体含义为:
如图1所示为所设计的膜系统结构,其运行过程如下:所述3个膜中,
(3)gij计算出以后,规则执行条件被满足,因此被执行,计算图像上每像素的梯度与设定阈值的绝对差,并将差赋给diffij和酶Eij,由于是1∶1的等比例分配,差diffij和酶Eij的值是相同的。若Eij>0,表示像素(i,j)是边缘,反之不是边缘。
(6)在膜main中,Pr1,main被执行,ET被赋值为1,整个系统停止运行。
如图2所示,现以分辨率为5*5pixel的图像为例,进一步说明本发明的执行过程与中间结果。原始图像如图2(a)所示。
(4)酶Eij实际反映了梯度和阈值的差,其结果可能为正,可能为负。对于规则当Eij为正时,其值大于ed1,被执行,ed2的值被分配给edgij。当Eij为负时,其值大于ed3,被执行,0被分配给edgij。和被执行后,得到图2(e)结果。同时酶被产生。
(6)EH产生后,膜Main中的规则Pr1,main被满足,执行后,控制系统运行的酶变量ET被赋值为1,整个系统停止运行,输出edg。
以上所述,仅为本发明较佳的具体实施方式,本发明的保护范围不限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可显而易见地得到的技术方案的简单变化或等效替换均落入本发明的保护范围内。
Claims (2)
1.一种Robert边缘检测膜系统,其特征在于,度为3,包括3个膜,具体为:
∏=(m,H,μ,(Var1,E1,Pr1,Var1(0)),(Var2,E2,Pr2,Var2(0)),(Var3,E3,Pr3,Var3(0)))
其中:(1)m为∏的度,代表酶数值膜系统中膜的层次数量,m=3;
(2)H为膜标号,本系统采用3个膜,标号分别为1,2,3;
(3)μ=[[[ ]3]2]1代表酶数值膜的结构;
(4)Var1,Var2,Var3代表3个膜内的所有对象变量的集合,其对应的初始值为Var1(0),Var2(0),Var3(0);Var1={xij,edgij,ed1,ed2,ed3,CT,EH},其中,xij表示图像I中,坐标为(i,j)元素的灰度值,对应的初始值分别为:xij=vij表示图像的灰度值;edgij表示边缘检测结束后,坐标为(i,j)元素的灰度值,取值范围为0,1;初始值为edgij=0;ed1=0;ed2=1;ed3=-256;CT数值变量,初始值CT=1;EH为酶变量,初始值为空,EH=[],表示图像I中是否所有像素位置处的边缘被计算出,若所有值被计算,其值不为空,反之为空;Var2={fxij,fyij,gij),fxij、fyij、gij分别表示图像I中坐标(i,j)处的水平、垂直与总体梯度值,对应的初始值均为空;diffij表示图像I中坐标(i,j)处总体梯度与设定阈值的差,对应的初始值为空;Eij酶变量,表示对应的diffij是否被计算出,若被计算出,其值不为空,反之为空;酶变量,表示图像I中坐标(i,j)处的边缘是否被计算出,若被计算,其值不为空,反之为空;
(6)Pri,表示膜i内的规则集,规则由两部分构成,一部分是值产生规则,一部分是值分配规则。
2.如权利要求 1所述一种Robert边缘检测膜系统,其特征在于,所述3个膜中,膜Main中的变量为图像每个元素,其初始值为图像灰度值;
(3)gij计算出以后,规则执行条件被满足,因此被执行,计算图像上每像素的梯度与设定阈值的绝对差,并将差赋给diffij和酶Eij,由于是1∶1的等比例分配,差diffij和酶Eij的值是相同的;若Eij>0,表示像素(i,j)是边缘,反之不是边缘;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810429215.3A CN108648206B (zh) | 2018-04-28 | 2018-04-28 | 一种Robert边缘检测膜计算系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810429215.3A CN108648206B (zh) | 2018-04-28 | 2018-04-28 | 一种Robert边缘检测膜计算系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108648206A CN108648206A (zh) | 2018-10-12 |
CN108648206B true CN108648206B (zh) | 2022-09-16 |
Family
ID=63749681
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810429215.3A Expired - Fee Related CN108648206B (zh) | 2018-04-28 | 2018-04-28 | 一种Robert边缘检测膜计算系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108648206B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103218605A (zh) * | 2013-04-09 | 2013-07-24 | 东南大学 | 一种基于积分投影与边缘检测的快速人眼定位方法 |
CN105493105A (zh) * | 2013-07-24 | 2016-04-13 | 意大利电信股份公司 | 关键点识别 |
CN107729644A (zh) * | 2017-06-23 | 2018-02-23 | 安徽理工大学 | 一种利用膜计算提高时间有效性实现最小支配集的方法 |
CN107807669A (zh) * | 2017-12-01 | 2018-03-16 | 西华大学 | 基于酶数值膜结构的控制方法、装置及可读取存储介质 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090124506A1 (en) * | 2007-11-02 | 2009-05-14 | Technion Research & Development Foundation Ltd. | Membrane computing |
-
2018
- 2018-04-28 CN CN201810429215.3A patent/CN108648206B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103218605A (zh) * | 2013-04-09 | 2013-07-24 | 东南大学 | 一种基于积分投影与边缘检测的快速人眼定位方法 |
CN105493105A (zh) * | 2013-07-24 | 2016-04-13 | 意大利电信股份公司 | 关键点识别 |
CN107729644A (zh) * | 2017-06-23 | 2018-02-23 | 安徽理工大学 | 一种利用膜计算提高时间有效性实现最小支配集的方法 |
CN107807669A (zh) * | 2017-12-01 | 2018-03-16 | 西华大学 | 基于酶数值膜结构的控制方法、装置及可读取存储介质 |
Non-Patent Citations (5)
Title |
---|
Area Extraction of beads in Membrane filter using Image Segmentation Techniques;Neeti Taneja 等;《International Journal of Computer Science and Information Technologies》;20141231;第5卷(第6期);7799-7803 * |
Segmenting images with gradient-based edge detection using Membrane Computing;Daniel Díaz-Pernil 等;《Pattern Recognition Letters》;20130630;第34卷(第8期);846-855 * |
Universal enzymatic numerical P systems with small number of enzymatic variables;Zhiqiang Zhang 等;《Sci. China Inf. Sci》;20171120;38-49 * |
膜计算在图像处理领域应用研究综述;袁建英 等;《安徽大学学报(自然科学版)》;20180330;第42卷(第3期);24-31 * |
酶数值膜系统仿真器设计与实现;王争一;《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》;20170115(第01期);I137-50 * |
Also Published As
Publication number | Publication date |
---|---|
CN108648206A (zh) | 2018-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Dhingra et al. | The proximal augmented Lagrangian method for nonsmooth composite optimization | |
CN101271572B (zh) | 基于免疫克隆选择聚类的图像分割方法 | |
CN109165735B (zh) | 基于生成对抗网络与自适应比例生成样本图片的方法 | |
CN103745482B (zh) | 一种基于蝙蝠算法优化模糊熵的双阈值图像分割方法 | |
CN108460783B (zh) | 一种脑部核磁共振图像组织分割方法 | |
WO2023201772A1 (zh) | 基于迭代域内适应和自训练的跨域遥感图像语义分割方法 | |
CN112836820B (zh) | 用于图像分类任务的深度卷积网络训练方法、装置及系统 | |
CN113052755A (zh) | 一种基于深度学习的高分辨率图像智能化抠图方法 | |
Pérez-Patricio et al. | An FPGA stereo matching unit based on fuzzy logic | |
CN112988815A (zh) | 一种大规模高维高速流数据在线异常检测的方法及系统 | |
CN110880773A (zh) | 基于数据驱动与物理模型驱动结合的电网调频控制方法 | |
Hao et al. | Canny edge detection enhancement by general auto-regression model and bi-dimensional maximum conditional entropy | |
CN113392854A (zh) | 一种图像纹理特征提取分类方法 | |
CN104036503B (zh) | 基于空间位置信息的图像分割方法 | |
CN108648206B (zh) | 一种Robert边缘检测膜计算系统及方法 | |
US20200092529A1 (en) | Image processing apparatus and 2d image generation program | |
CN111460974B (zh) | 一种基于优化的散乱点云数据全局特征提取方法 | |
CN114581454B (zh) | 基于背景差分法的量子图像分割方法、装置及存储介质 | |
CN110457155A (zh) | 一种样本类别标签的修正方法、装置及电子设备 | |
CN109389565A (zh) | 一种保持边缘和平滑度可调节的图像平滑方法 | |
Qiao et al. | Boosting federated learning convergence with prototype regularization | |
CN114863114A (zh) | 一种基于改进的混沌粒子群的图像阈值分割方法 | |
Liu et al. | Effective 3D boundary learning via a nonlocal deformable network | |
Ma | Inpainting technology in art pictures based on edge and texture feature | |
CN109190655B (zh) | 一种基于酶数值膜系统的ncc图像匹配方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220916 |