CN111951291B - 基于多结构形态学与fodpso混合处理红外图像边缘检测方法 - Google Patents
基于多结构形态学与fodpso混合处理红外图像边缘检测方法 Download PDFInfo
- Publication number
- CN111951291B CN111951291B CN202010810993.4A CN202010810993A CN111951291B CN 111951291 B CN111951291 B CN 111951291B CN 202010810993 A CN202010810993 A CN 202010810993A CN 111951291 B CN111951291 B CN 111951291B
- Authority
- CN
- China
- Prior art keywords
- image
- fodpso
- structural element
- value
- particle
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000012545 processing Methods 0.000 title claims abstract description 46
- 238000003708 edge detection Methods 0.000 title claims abstract description 37
- 230000000694 effects Effects 0.000 claims abstract description 26
- 230000000877 morphologic effect Effects 0.000 claims abstract description 23
- 238000000605 extraction Methods 0.000 claims abstract description 14
- 239000002245 particle Substances 0.000 claims description 101
- 230000007797 corrosion Effects 0.000 claims description 19
- 238000005260 corrosion Methods 0.000 claims description 19
- 238000001914 filtration Methods 0.000 claims description 17
- 230000004927 fusion Effects 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 5
- 241001270131 Agaricus moelleri Species 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000007547 defect Effects 0.000 abstract description 44
- 230000001131 transforming effect Effects 0.000 abstract 1
- 230000011218 segmentation Effects 0.000 description 10
- 238000001514 detection method Methods 0.000 description 9
- 238000003709 image segmentation Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 239000000463 material Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000003044 adaptive effect Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000001931 thermography Methods 0.000 description 4
- VZSRBBMJRBPUNF-UHFFFAOYSA-N 2-(2,3-dihydro-1H-inden-2-ylamino)-N-[3-oxo-3-(2,4,6,7-tetrahydrotriazolo[4,5-c]pyridin-5-yl)propyl]pyrimidine-5-carboxamide Chemical compound C1C(CC2=CC=CC=C12)NC1=NC=C(C=N1)C(=O)NCCC(N1CC2=C(CC1)NN=N2)=O VZSRBBMJRBPUNF-UHFFFAOYSA-N 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000012002 interactive response technology Methods 0.000 description 2
- 238000002156 mixing Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- QNRATNLHPGXHMA-XZHTYLCXSA-N (r)-(6-ethoxyquinolin-4-yl)-[(2s,4s,5r)-5-ethyl-1-azabicyclo[2.2.2]octan-2-yl]methanol;hydrochloride Chemical compound Cl.C([C@H]([C@H](C1)CC)C2)CN1[C@@H]2[C@H](O)C1=CC=NC2=CC=C(OCC)C=C21 QNRATNLHPGXHMA-XZHTYLCXSA-N 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000005530 etching Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 230000004660 morphological change Effects 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000007500 overflow downdraw method Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000011179 visual inspection Methods 0.000 description 1
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
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
- G06T5/30—Erosion or dilatation, e.g. thinning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- 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/155—Segmentation; Edge detection involving morphological operators
-
- 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/10048—Infrared image
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computational Linguistics (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
基于多结构形态学与FODPSO混合处理的红外图像边缘提取方法,它涉及红外图像边缘提取方法。它是要解决现有的红外图像边缘提取方法的图像周边存在无用噪声、边界不真实、对图像缺陷特征信息提取效果差的技术问题。本方法:将原始图像灰度变换,然后用不同的结构元素腐蚀并计算出权值向量;构建并联滤波器进行自适应权值多结构形态学去噪处理,得到去噪后图像;然后用FODPSO算法处理,得到分割后图像;再构建新的结构元素进行不同刻度的边缘提取,最终进行融合,从而提取到最终完整的边缘检测图像。本方法能有效滤除图像噪声,提高图像对比度,分割后图像边缘轮廓更平滑、准确,提高了对缺陷特征提取的效果,用于红外图像处理领域。
Description
技术领域
本发明涉及红外图像边缘提取方法。
背景技术
红外热成像(Infrared Thermography,IRT)作为一种新兴的无损检测(Non-destructive testing,NDT)技术,可检测材料表面的温度变化。IRT作为一种快速的非接触式大面积可视化检查技术越来越受到人们的认可,已成功地在许多领域中用作NDT技术。其工作原理在于:材料内部缺陷(如裂纹、脱粘、磨损等)导致的不连续性会造成激励能量在材料不连续处上表面产生明显温度对比。理论上,材料内部的任何不连续都会导致这种温度变化,这种变化可以通过红外热像仪来检测。针对进行缺陷试件红外热波检测时,红外热像仪在采集的红外序列热图的过程中,由于实验设备和环境自身的红外辐射产生的背景噪声、材料内部结构的不均匀以及材料表面受热不均产生的不利影响,导致采集到的序列热图信号中缺陷的信息被大量的无关噪声所覆盖,这些噪声往往会严重影响序列热图信号的分析精度及其处理分析结果,增加了后续红外图像处理的难度。因此需要对红外图像处理技术进行改进,减少红外图像中无关噪声的影响,改善对红外图像的分析效果。
数学形态学以图像的形态特征为研究对象,通过膨胀、腐蚀、开启与闭合运算来描述图像的基本特征与结构。脉冲噪声是一种常见的图像噪声,根据噪声位置灰度值与其领域灰度值的比较可以分为正、负脉冲。正脉冲噪声的位置其灰度值要大于领域的灰度值,负脉冲则相反。开运算先腐蚀后膨胀,可以过滤图像中的正脉冲噪声;闭运算先膨胀后腐蚀,可以过滤图像中的负脉冲噪声。在图像去噪的过程中,通过适当地选取结构元素的形状和维数可以提高滤波去噪的效果。
图像分割是图像处理中的一项关键技术,主要用于把图像分割成若干具有某种相同或相似特性的区域,精确定位这些区域并对其特性进行进一步分析和处理。随着图像处理技术的发展,图像分割技术已广泛地应用在诸多领域中。当前,在红外图像处理中,阈值图像分割算法对噪声敏感度高,对图像灰度差异不明显以及不同缺陷目标灰度值有重叠分割不明显,因此其更加适用于相差度很大的不同缺陷目标和背景的图像。生长区域法图像分割需要把具有某种相似性质的像素连通,进而构成最终的分割区域,这样会造成红外图像缺陷的过度分割。模糊C均值聚类算法FCM更适合存在不确定性和模糊性的缺陷图像,具有很大的局限性。遗传算法GA、粒子群算法PSO等算法在运算中很容易陷于局部最优解,其分割结果不具有全局最优解的效果。
发明内容
本发明的目的是提供一种基于多结构形态学与FODPSO混合处理红外图像边缘检测方法,达到减少图像周边大量无用噪声和不真实的边界信息,提高对图像缺陷特征信息提取的效果。
本发明的基于多结构形态学与FODPSO混合处理红外图像边缘检测方法,其算法流程图如图1所示,包括以下步骤:
(1)对原始图像In进行灰度变换,得到灰度图像f;
(2)构建不同的结构元素Ai,其中i=1、2、……、I;I=3~5;对变换后的灰度图像f进行I种结构元素腐蚀运算,得到腐蚀效果图Imei;
(3)将各结构元素Ai的权值相应地记为αi,根据其匹配图像的次数βi来计算权值αi;
(4)对I种结构元素Ai进行扩维,得到新的结构元素Bi,其中Bi的维数比Ai的维数多1~4;通过串联构建同种元素串联滤波器,串联是指将同一形状的结构元素按照维数从小到大的顺序来对图像进行滤波,得到I种串联滤波结果Imi;
(5)根据串联滤波结果Imi与灰度图像f进行差异值计算,将结果值作为权值向量
(6)将I种结构元素所构建的串联滤波器进行并联,形成并联滤波器,通过串联滤波结果Imi、权值向量进行并联加权求和计算,得到自适应权值多结构形态学并联去噪结果图像Φ(x,y);
(7)初始化分数阶粒子群FODPSO算法参数,进行粒子速度与个体位置更新;针对特定的红外缺陷图像特征,设置合适的参数,实现对去噪后图像Φ(x,y)的FODPSO算法处理,进而得到分割后图像fd;
(8)由于较大维数的结构元素会对图像缺陷的几何特征造成影响,导致形态处理结果不准确,因此对第(2)步聚的结构元素Ai进行降维,并增加一种圆盘型结构元素,构建出新的结构元素Cj,其中j=1、2、3、……、J;J=I+1;其中C1、C2、……、CI分别是A1、A2、……、AI降维得到的;CJ为圆盘型结构元素,圆盘型结构元素CJ的维数与A1、A2、……、AI的维数相等;
(9)对多结构形态学算法进行变化,计算J种结构元素在不同刻度下开闭、闭开运算下的均值图像qs(x,y),其中s=1、2、3、……、J;
(10)计算J种结构元素在不同刻度下的图像方差其中s=1、2、3、……、J;
(11)根据图像方差计算图像权值us,其中s=1、2、3、……、J;
(12)对J种结构元素Cj进行膨胀运算,获得J种结构元素在不同刻度下的边缘融合图像Gfj;
(13)将J个边缘融合图像Gf1、Gf2、Gf3、……、GfJ进行图像边缘检测结果融合,得到最终边缘提取融合图像Gf。边缘提取融合图像Gf具有更加完整的缺陷特征信息。
本发明在两次进行形态学处理时要用到结构元素,一是利用形态学算法进行图像去噪时,选取结构元素的方式是根据具体图像特征来构建适合图像的结构元素形状;二是利用形态学算法进行图像边缘检测时,依据分割后图像中对象特征,构建适合图像不同方向边缘检测的结构元素进行不同刻度运算,得到最终检测结果图像。
更进一步地,所述原始图像In为红外图像序列中一幅最优图像;
更进一步地,所述结构元素Ai为同一形状、不同维数与角度的结构元素;
更进一步地,结构元素Ai结构可为扁平型、菱型、圆盘型、八边型、线型或矩型结构元素。
结构元素的形状是影响其滤波效果的重要因素,根据相似性原则,结合输入的红外图像的缺陷特征,并保证缺陷边缘的滤除噪声效果;更优选的是线型结构元素,线型结构元素能够更好的实现对红外图像缺陷周围噪声的有效滤除;
对于无噪声或噪声影响较小的图像,结构元素的维数一般选取1~3,对于噪声影响严重的图像,结构元素的维数一般选取5~7,方向依据图像缺陷特征进行适当选取;本发明更进一步优选的A1、A2、A3、A4的维数与方向如下:
A1的维数为5,方向为-45°;
A2的维数为5,方向为45°;
A3的维数为3,方向为90°;
A4的维数为3,方向为0°;
所述结构元素A1,A2,A3,A4的半径由图像中对象的半径来确定,对图像中对象进行试验测取,更进一步地取结构元素Ai直径为4~10。
更进一步地,所述结构元素Bi是由结构元素Ai扩维得到的,为同一形状与角度、不同维数的结构元素;
对于无噪声或噪声影响较小的图像,结构元素的维数一般选取1~3,对于噪声影响严重的图像,结构元素的维数一般选取5~7,方向依据图像缺陷特征进行适当选取;本发明更进一步优选择的结构元素B1、B2、B3、B4维数与方向如下:
B1的维数为7,方向为-45°;
B2的维数为7,方向为45°;
B3的维数为5,方向为90°;
B4的维数为5,方向为0°;
结构元素Cj,j=1、2、3、4、5时的结构元素C1,C2,C3,C4,C5中,C1、C2、C3、C4为线型结构元素,C5为圆盘型结构元素;
其中C1、C2、C3、C4的维数与方向如下:
C1的维数为3,方向为90°;
C2的维数为3,方向为0°;
C3的维数为3,方向为45°;
C4的维数为3,方向为-45°;
圆盘型结构元素C5的维数为3;
所述结构元素Cj的半径由图像中对象的半径来确定,更进一步地取结构元素直径为4~10。
更进一步地,步骤(2)中所述的腐蚀运算的公式为:
式中,f(x,y)为输入图像,g(m,n)为结构元素;
更进一步地,步骤(5)中所述的权值向量的构造方式为:
更进一步地,步骤(6)中所述的整合串联结果Imi是对串联结果进行加权求和,得到去噪后图像记为Φ(x,y),即:
更进一步地,步骤(7)中所述的初始化分数阶粒子群FODPSO算法参数,进行粒子速度与个体位置更新,对去噪后图像Φ(x,y)进行FODPSO算法处理,得到分割后图像fd的算法流程图如图2所示,具体操作步骤如下:
①.初始化粒子群FODPSO算法参数,设定当前粒子群的粒子数N,设定最低粒子数POPmin,最高粒子数POPmax,设定当前粒子群数E,最低粒子群数Emin,最高粒子群数Emax;设定学习因子c1、c2;设定最大迭代次数Tmax,分数阶系数α,随机生成粒子速度Vw与位置Xx、粒子最大飞翔速度Vmax、粒子最小飞翔速度Vmin;其中,对N定义如下:
N=POPmin+(POPmax-POPmin)×rand
POPmin=0,POPmax=L
常见的nint8类型的红外图像,其灰度值为0~255;L为最大灰度值,L=256;
②.对每个粒子进行个体适应度值T(N)计算:
k代表一个具体的粒子群类,wk表示第k类的概率,μk为第k类的均值,μT表示图像灰度平均值;
③.对②中所计算的每个粒子进行适应度值评估,将种群内粒子的最优个体值Pd设置为当前位置,种群最优个体值Pg设置为初始群体中最佳粒子的位置;
④.对粒子速度进行更新,如下:
上式中t代表当前迭代次数;c1、c2为学习因子,且满足3<c1+c2<4;r1、r2为[0,1]内均匀分布的随机数;系数取值根据下面的二项式来确定:
⑤.对粒子速度进行约束:
⑥.考虑到过去粒子对当前粒子的影响,对粒子进行速度更新:
Vw(t-3)=Vw(t-2)
Vw(t-2)=Vw(t-1)
Vw(t-1)=Vw(t)
⑦.对粒子个体位置进行更新:
Xx(t+1)=Xx(t)+Vw(t+1)
⑧.对更新后粒子的适应度值进行计算,并与上次迭代最优适应度值比较,判断当前种群所有个体对应的适应度值的最小值是否小于上一次的种群最优个体适应度值;
如果是,则替换上一次种群最优个体值,接着判断是否达到最大迭代次数,如果已经达到最大迭代次数,则以当前种群所有个体对应的适应度值的最小值作为最优种群个体位置结果值,得到分割后图像fd;如果未达到最大迭代次数,返回到第②步继续执行;
如果否,则以当前种群所有个体对应的适应度值的最小值作为最优种群个体位置结果值,得到分割后图像fd。
更进一步地,步骤①中,设定当前粒子群的粒子数N=20~30,最低粒子数POPmin=1~10,最高粒子数POPmax=40~50;设定当前粒子群数E=4~5,最低粒子群数Emin=1~3,最高粒子群数Emax=5~6;设定学习因子c1=1.5~1.8,学习因子c2=1.5~1.8;设定最大迭代次数Tmax=50~150,分数阶系数α=0.5~0.8,随机生成粒子速度Vw与位置Xx,粒子最大飞翔速度Vmax=-1.5,粒子最小飞翔速度Vmin=1.5。
本发明的基于多结构形态学与FODPSO混合处理的红外图像边缘提取方法,首先对原始图像In进行灰度变换,得到灰度图像f;其次根据灰度图像噪声特点,构建不同的结构元素Ai对灰度图像f进行多结构元素腐蚀运算,得到图像腐蚀效果图Imei,通过计算匹配次数βi来确定各自权值αi,对结构元素Ai进行扩维得到结构元素Bi,通过串联得到串联滤波结果Imi,根据串联结果Imi进行差异值计算,得到权值向量其次构建并联滤波器进行自适应权值多结构形态学去噪处理,得到去噪后图像Φ(x,y);本步骤结合形态学腐蚀过程中结构元素的探测匹配原理,根据结构元素在图像中可匹配次数来进行自适应权值计算,实现权重自适应的多结构形态学去噪。然后初始化FODPSO参数,对去噪后图像Φ(x,y)进行分数阶粒子群FODPSO算法处理,得到分割后图像fd;本步骤选择基于多结构形态学与FODPSO混合处理算法,能够在整体上克服各种单一算法的不足,改善图像缺陷检测的效果。最后根据分割图像fd,对多结构形态学算法进行变化,提出基于多结构形态学算法在多刻度上的红外图像边缘检测算法,依据结构元素Ai变化得到结构元素Cj,采用对结构元素Cj进行膨胀运算对分割后图像fd进行不同刻度的边缘提取,得到不同刻度下的最终结果Gf1,Gf2,Gf3,Gf4,Gf5,将最终结果Gf1,Gf2,Gf3,Gf4,Gf5进行融合,从而提取到最终完整的边缘检测图像Gf,其具有更加完整的缺陷特征信息。本步骤基于权重自适应的多结构形态学去噪后图像,进行分数阶粒子群FODPSO图像分割,得到更加优化的红外图像。对优化后图像进行膨胀运算,通过多结构元素对图像进行边缘检测,通过加权融合的思想来整合检测到的边缘,可以在一定程度上提高边缘检测的精度,具有较好的抗干扰性和定位准确性。
本发明与现有处理方法相比较,所提出的混合算法具有如下有益效果:
(1)本发明采用自适应权值多结构形态学去噪算法对图像进行滤波去噪,能够有效滤除图像噪声,缺陷区域的信息更加清晰,提高图像对比度,便于后续图像分割操作。
(2)本发明采用FODPSO算法对图像进行分割,分割后图像的缺陷边缘轮廓更平滑、尺寸更准确,几何特征更加明显。
(3)本发明采用多结构形态学与FODPSO混合处理算法对图像进行边缘检测,与只用Canny算子进行边缘检测结果对比,经混合算法处理后的图像缺陷边缘提取准确、缺陷特征清晰明显、滤除大量的冗余边界信息、有效滤除了无用信息的干扰、提高了对缺陷特征提取的效果。
本发明可应用于红外热波无损检测领域,同时应用于红外图像处理技术领域,主要用于红外图像处理中的图像边缘检测处理。
附图说明
图1为本发明的基于多结构形态学与FODPSO混合处理的红外图像边缘提取方法流程图;
图2为本发明步骤(7)中FODPSO算法流程图;
图3为实施例1中步骤(2)中的4种结构元素A1,A2,A3,A4示意图;
图4为实施例1中步骤(4)中的A1,A2,A3,A4维数扩增后的B1,B2,B3,B4示意图;
图5为实施例1中变换后的灰度图像;
图6为实施例1中4种结构元素对图像腐蚀效果图;
图7为实施例1中自适应权值多结构形态学去噪后图像;
图8为实施例1中本算法FODPSO分割后图像;
图9为实施例1中对比PSO分割后图像;
图10为实施例1步骤(8)四中的5种结构元素C1,C2,C3,C4,C5示意图;
图11为实施例1中为多结构形态学与FODPSO混合处理算法对分割后图像进行缺陷边缘提取融合效果图;
图12为实施例1中只用Canny算子进行边缘检测效果图;
具体实施方式
为了更好地说明本发明的技术方案,下面结合附图对本发明的技术方案作进一步的说明,但并不局限于此,凡是对本发明技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,均应涵盖在本发明的保护范围内。
实施例1:本实施例的基于多结构形态学与FODPSO混合处理的红外图像边缘提取方法,按以下步骤进行:
(1)读取原始图像In,转换为灰度图像f;
(2)构建不同的结构元素Ai,其中i=1、2、3、4;其中A1,A2,A3,A4为同一形状、不同维数与角度的结构元素,有利于滤除同种噪声影响;A1,A2,A3,A4的图形如图3所示,其中A1的维数为5,方向为-45°;
A2的维数为5,方向为45°;
A3的维数为3,方向为90°;
A4的维数为3,方向为0°;
结构元素A1,A2,A3,A4的直径为4~10,在算法运行中进行取值变化,附图中并未表示出直径大小变化。
利用A1,A2,A3,A4对变换后图像f进行腐蚀运算,运算公式如下:
式中,f(x,y)为输入图像,g(m,n)为结构元素;
腐蚀运算后,得到腐蚀效果图Ime1,Ime2,Ime3,Ime4;腐蚀运算过程相当于对图像中可以匹配结构元素的位置进行探测并标记处理。
(3)将4种线型形状的结构元素A1,A2,A3,A4的权值依次记为α1,α2,α3,α4,在图像进行腐蚀的过程中,各结构元素匹配图像的次数依次为β1,β2,β3,β4,按自适应权值计算公式计算权值依次记为α1,α2,α3,α4,
(4)对4种结构元素A1,A2,A3,A4进行扩维,得到新的结构元素B1,B2,B3,B4;扩维结构元素B1,B2,B3,B4如图4所示,B1、B2、B3、B4的维数与方向如下:
B1的维数为7,方向为-45°;
B2的维数为7,方向为45°;
B3的维数为5,方向为90°;
B4的维数为5,方向为0°;
结构元素B1,B2,B3,B4的直径为4~10,在算法运行中进行取值变化,图4中并未表示出直径大小变化。
通过串联构建同种元素串联滤波器,对4种腐蚀运算结果图像Ime1,Ime2,Ime3,Ime4进行串联滤波,得到经过4种线型结构元素腐蚀的串联滤波结果Im1,Im2,Im3,Im4;
(5)将串联处理结果Im1,Im2,Im3,Im4与原始灰度图像f进行差异值计算,其结果值作为权值向量,记为
(6)将4种不同结构元素串联滤波器进行并联,形成并联滤波器,对串联滤波结果Im1,Im2,Im3,Im4、权值向量进行加权求和,以此得到并联去噪后结果图像记为Φ(x,y),即:
传统的形态学算法基本形态变化包括腐蚀、膨胀、开运算以及闭运算,运用传统的形态学算法对图像进行处理后,图像中仍然会存在大量与缺陷区域重合的短线和孤立噪声点。因此本算法对传统形态学处理算法进行了改进,改进后的多结构形态学处理算法,构建不同角度的结构元素A1,A2,A3,A4对图像进行腐蚀运算,得到腐蚀图Ime1,Ime2,Ime3,Ime4,并计算权值α1,α2,α3,α4,对结构元素A1,A2,A3,A4扩维变为新结构元素B1,B2,B3,B4,对腐蚀图Ime1,Ime2,Ime3,Ime4进行串联滤波得到串联滤波结果Im1,Im2,Im3,Im4,通过差异值计算得到权值向量最后通过并联滤波整合串联滤波结果得到去噪后图像Φ(x,y)。改进后的形态学处理算法相比于传统的形态学处理算法,其去噪效果更好。
(7)初始化分数阶粒子群FODPSO算法参数,进行粒子速度与个体位置更新,对去噪后图像Φ(x,y)进行FODPSO算法处理,得到分割后图像fd,具体如下:
①.初始化粒子群FODPSO算法参数,设定当前粒子群的粒子数N,设定最低粒子数POPmin,最高粒子数POPmax,设定当前粒子群数E,最低粒子群数Emin,最高粒子群数Emax;设定学习因子c1、c2;设定最大迭代次数Tmax,分数阶系数α,随机生成粒子速度Vw与位置Xx,粒子最大飞翔速度Vmax,粒子最小飞翔速度Vmin。其中,对N定义如下:
常见的nint8类型的红外图像,其灰度值为0~255。最大灰度值为L=256。
设定当前粒子群的粒子数为N=30,最低粒子数POPmin=10,最高粒子数POPmax=50;设定当前粒子群数E=4,最低粒子群数Emin=2,最高粒子群数Emax=6;设定学习因子c1=c2=1.8;设定最大迭代次数Tmax=150,分数阶系数α=0.6,随机生成粒子速度Vw与位置Xx,粒子最大飞翔速度Vmax=-1.5,粒子最小飞翔速度Vmin=1.5。
②.对每个粒子进行个体适应度值T(N)计算:
k代表一个具体的粒子群类,wk表示第k类的概率,μk为第k类的均值,μT表示图像灰度平均值;
③.对②中所计算的每个粒子进行适应度值评估,将种群内粒子的最优个体值Pd设置为当前位置,种群最优个体值Pg设置为初始群体中最佳粒子的位置;
④.对粒子速度进行更新,如下:
式(10)中t代表当前迭代次数;c1、c2为学习因子,且满足3<c1+c2<4;r1、r2为[0,1]内均匀分布的随机数;式中的系数取值可根据二项式(11)来计算和理解:
⑤.对粒子速度进行约束:
⑥.考虑到过去粒子对当前粒子的影响,对粒子进行速度更新:
⑦.对粒子个体位置进行更新:
Xx(t+1)=Xx(t)+Vw(t+1) (14)
⑧.对更新后粒子的适应度值进行计算,并与上次迭代最优适应度值比较,判断当前种群所有个体对应的适应度值的最小值是否小于上一次的种群最优个体适应度值,如果是,则替换上一次种群最优个体值,接着判断是否达到最大迭代次数,如果已经达到最大迭代次数,则以当前种群所有个体对应的适应度值的最小值作为最优种群个体位置结果值,得到分割后图像fd;如果未达到最大迭代次数,返回到第②步继续执行;如果否,则以当前种群所有个体对应的适应度值的最小值作为最优种群个体位置结果值,得到分割后图像fd。
(8)构建出新的结构元素:
①.为了能够尽可能匹配图像不同方向上的边缘,需要对前述的4种结构元素A1,A2,A3,A4进行变化,构建适合本处算法的5种结构元素C1,C2,C3,C4,C5来对图像进行边缘检测,如图10所示,其中C1、C2、C3、C4的维数与方向如下:
C1的维数为3,方向为90°;
C2的维数为3,方向为0°;
C3的维数为3,方向为45°;
C4的维数为3,方向为-45°;
结构元素C5为圆盘型结构元素,维数为3;
结构元素C1、C2、C3、C4、C5的直径为4~10,在算法运行中进行取值变化,附图中并未表示出直径大小变化;
同时定义5种结构元素在不同刻度上的变化为:
式中:n为刻度参数。
多结构元素在多刻度上的灰度值开运算定义为:
多结构元素在多刻度上的灰度值闭运算定义为:
则多结构元素在多刻度上的边缘检测算法定义为:
式18中所述J:J=1、2、3、4、5。
②.为了保证5种结构元素在不同刻度上所检测到的缺陷边缘能够进行图像融合,采用图像加权融合的方法来进行图像融合,对多结构元素在多刻度上的边缘融合算法定义如下:
式19中,us为多结构元素在不同刻度上边缘检测图像进行融合时的加权系数,[1,n]为刻度范围。
(9)对多结构形态学算法进行变化,计算J种结构元素在不同刻度下开闭、闭开运算下的均值图像qs(x,y),具体如下:,
一般u1=u2=…=uk=1/n,此时,多结构元素在不同刻度上所检测的边缘对最终边缘贡献一样大,但这样提取的边缘效果不是很好。为了能更好地检测最佳边缘图像,取非均值权值,基本原则是依据多结构元素在不同刻度上检测性能的不同来确定。大刻度的结构元素检测缺陷直径大的能力强,小刻度的结构元素检测缺陷直径小的能力强,因而可将大刻度的权重取得大一些,小刻度取得小一些。按照以上原则,本步骤设计如下计算不同权值的算法步骤:
运用第1小步的灰度值开运算与闭运算,分别求得5种结构元素在不同刻度下开闭、闭开运算下的均值图像:
式中,qs(x,y)表示均值图像。
(10)计算5种结构元素在不同刻度下的图像方差
(11)以各个刻度下的图像方差的比例来确定权值us:
(12)对5个结构元素分别进行膨胀,用膨胀后的5个结构元素C1,C2,C3,C4,C5在刻度n=1时,对图像进行边缘检测,将获得的5个边缘信息检测结果依据式(19),只计算n=1时的边缘融合图像,获得多结构元素在单刻度下的边缘检测结果Gf1。其中,多结构形态学算法在不同刻度下的膨胀运算定义为:
对5个结构元素分别进行膨胀,用膨胀后的5个结构元素C1,C2,C3,C4,C5在刻度n=2时,对图像进行边缘检测,将获得的5个检测结果按照步骤(12)进行图像边缘融合,获得多结构元素在刻度为n=2时的边缘检测结果Gf2。
同理,按照上述步骤对5个结构元素C1,C2,C3,C4,C5分别进行膨胀,进行刻度n=3、4、5时的边缘检测及融合,获得多结构元素在刻度n=3、4、5时的边缘融合图像:Gf3,Gf4,Gf5。
(13)将5个边缘融合图像Gf1、Gf2、Gf3、Gf4、Gf5按式(19)进行多结构元素在多刻度上的图像边缘检测结果融合,得到最终边缘提取融合图像Gf。边缘检测图像Gf具有更加完整的缺陷特征信息。对多结构形态学算法进行变化,首先对4种结构元素进行变化得到5种不同结构元素C1,C2,C3,C4,C5,然后通过采用不同刻度的结构元素进行边缘检测,再通过加权融合的思想来整合检测到的边缘,可以得到基于多结构形态学算法在多刻度上的红外图像边缘检测算法,这种算法能够有效地提高缺陷边缘检测的精度,可以更加直观的对红外图像中的缺陷进行辨析。
本实施例中,原始图像In为经脉冲相位法处理后的红外图像序列中效果最优图像。对原始图像In进行处理,首先进行灰度图转换,得到灰度图f,如图5所示。再进行自适应权值多结构形态学去噪处理,效果图Φ(x,y)如图6所示。其中(a)、(b)、(c)、(d)分别为通过4种结构元素A1,A2,A3,A4对灰度图像f进行腐蚀运算所得效果图Ime1,Ime2,Ime3,Ime4;图7为自适应权值多结构形态学去噪后图像Φ(x,y)。经过去噪后的图像噪声干扰显著降低,缺陷轮廓更加清晰,图像的对比度增加,有利于接下来的分数阶粒子群图像分割算法处理。FODPSO算法所处理的对象为图7所示的经自适应权值多结构形态学去噪后图像Φ(x,y),FODPSO分割算法处理后的分割效果图fd为图8所示,图9为PSO算法对比分割效果图。对比发现,直接采用PSO算法进行处理后的图像,缺陷周围仍然存在部分无关噪声信息,导致缺陷轮廓十分不清晰,难以进行缺陷边缘的提取。而本发明所采用的FODPSO算法处理经过多结构形态学算法处理过的去噪图像Φ(x,y),使得分割处理后图像fd的缺陷区域完整的分割出来,未产生过度分割,滤除掉了大量无关噪声,使得缺陷几何特征更加容易判断,缺陷边缘轮廓更加清晰,更有利于接下来的图像缺陷边缘提取。在经过上述的自适应权值多结构形态学对图像进行去噪,FODPSO算法对去噪后图像Φ(x,y)进行分割,多结构形态学在多刻度上的红外图像边缘检测算法对分割后图像fd进行边缘检测,效果图如图11所示,将上述混合算法处理后的图像与只用Canny算子边缘检测的图像12进行对比。从结果图12中可以看出,Canny算子虽然有较强的抑制噪声的效果,但是容易将红外图像的噪声作为边界,由Canny算子处理过后的边缘图像存在大量的冗余边界,难以有效地滤除边缘无关信息的干扰。因此,只有当红外图像噪声较小时,Canny算子边缘精度相对较高。由结果图11可知,经过多结构形态学与FODPSO混合处理算法处理后的图像,对缺陷边缘识别更加准确,缺陷轮廓更加清晰,滤除了周围边缘大量无用的边界信息,降低了噪声对缺陷边缘识别的影响,提高了缺陷边缘检测的能力。
Claims (10)
1.基于多结构形态学与FODPSO混合处理的红外图像边缘检测方法,其特征在于该方法包括以下步骤:
(1)对原始图像In进行灰度变换,得到灰度图像f;
(2)构建不同的结构元素Ai,其中i=1、2、……、I;I=3~5;对变换后的灰度图像f进行I种结构元素腐蚀运算,得到腐蚀效果图Imei;
(3)将各结构元素Ai的权值相应地记为αi,根据其匹配图像的次数βi来计算权值αi;
(4)对I种结构元素Ai进行扩维,得到新的结构元素Bi,其中Bi的维数比Ai的维数多1~4;通过串联构建同种元素串联滤波器,得到I种串联滤波结果Imi;
(5)根据串联滤波结果Imi与灰度图像f进行差异值计算,将结果值作为权值向量
(6)将I种结构元素所构建的串联滤波器进行并联,形成并联滤波器,通过串联滤波结果Imi、权值向量进行并联加权求和计算,得到自适应权值多结构形态学并联去噪结果图像Φ(x,y);
(7)初始化分数阶粒子群FODPSO算法参数,进行粒子速度与个体位置更新;对去噪后图像Φ(x,y)进行FODPSO算法处理,得到分割后图像fd;
(8)第(2)步聚的结构元素Ai进行降维,并增加一种圆盘型结构元素,构建出新的结构元素Cj,其中j=1、2、3、……、J;J=I+1;其中C1、C2、……、CI分别是A1、A2、……、AI降维得到的;CJ为圆盘型结构元素,圆盘型结构元素CJ的维数与A1、A2、……、AI的维数相等;
(9)对多结构形态学算法进行变化,计算J种结构元素在不同刻度下开闭、闭开运算下的均值图像qs(x,y),其中s=1、2、3、……、J;
(10)计算J种结构元素在不同刻度下的图像方差其中s=1、2、3、……、J;
(11)根据图像方差计算图像权值us,其中s=1、2、3、……、J;
(12)对J种结构元素Cj进行膨胀运算,获得J种结构元素在不同刻度下的边缘融合图像Gfj;
(13)将J个边缘融合图像Gf1、Gf2、Gf3、……、GfJ进行图像边缘检测结果融合,得到最终边缘提取融合图像Gf。
2.根据权利要求1所述的基于多结构形态学与FODPSO混合处理的红外图像边缘检测方法,其特征在于所述原始图像I为红外图像序列中一幅最优图像。
3.根据权利要求1或2所述的基于多结构形态学与FODPSO混合处理的红外图像边缘检测方法,其特征在于结构元素Ai结构为线型结构元素。
4.根据权利要求1或2所述的基于多结构形态学与FODPSO混合处理的红外图像边缘检测方法,其特征在于步骤(2)中I=4;且A1、A2、A3、A4的维数与方向如下:
A1的维数为5,方向为-45°;
A2的维数为5,方向为45°;
A3的维数为3,方向为90°;
A4的维数为3,方向为0°。
5.根据权利要求1或2所述的基于多结构形态学与FODPSO混合处理的红外图像边缘检测方法,其特征在于步骤(4)中所述结构元素Bi是由结构元素Ai扩维得到的,为同一形状与角度、不同维数的结构元素。
6.根据权利要求1或2所述的基于多结构形态学与FODPSO混合处理的红外图像边缘检测方法,其特征在于步骤(4)中所述的Bi由4个结构元素组成,B1、B2、B3、B4的维数与方向如下:
B1的维数为7,方向为-45°;
B2的维数为7,方向为45°;
B3的维数为5,方向为90°;
B4的维数为5,方向为0°。
7.根据权利要求1或2所述的基于多结构形态学与FODPSO混合处理的红外图像边缘检测方法,其特征在于步骤(8)中所述的结构元素Cj由5个结构元素C1,C2,C3,C4,C5组成,其中C1、C2、C3、C4为线型结构元素,C5为圆盘型结构元素;
其中C1、C2、C3、C4的维数与方向如下:
C1的维数为3,方向为90°;
C2的维数为3,方向为0°;
C3的维数为3,方向为45°;
C4的维数为3,方向为-45°;
圆盘型结构元素C5的维数为3。
8.根据权利要求1或2所述的基于多结构形态学与FODPSO混合处理的红外图像边缘检测方法,其特征在于步骤(2)中所述的腐蚀运算的公式为:
式中,f(x,y)为输入图像,g(m,n)为结构元素。
9.根据权利要求1或2所述的基于多结构形态学与FODPSO混合处理的红外图像边缘检测方法,其特征在于步骤(5)中所述的权值向量的构造方式为:
10.根据权利要求1或2所述的基于多结构形态学与FODPSO混合处理的红外图像边缘检测方法,其特征在于步骤(7)中所述的初始化分数阶粒子群FODPSO算法参数,进行粒子速度与个体位置更新,对去噪后图像Φ(x,y)进行FODPSO算法处理,得到分割后图像fd,具体操作步骤如下:
①.初始化粒子群FODPSO算法参数,选取当前粒子群的粒子数N,设定最低粒子数POPmin,最高粒子数POPmax,设定当前粒子群数E,最低粒子群数Emin,最高粒子群数Emax;设定学习因子c1、c2;设定最大迭代次数Tmax,分数阶系数α,随机生成粒子速度Vw与位置Xx,粒子最大飞翔速度Vmax,粒子最小飞翔速度Vmin;其中,对N定义如下:
N=POPmin+(POPmax-POPmin)×rand
POPmin=0,POPmax=L
上式中L为最大灰度值,L=256;
②.对每个粒子进行个体适应度值T(N)计算:
其中k代表一个具体的粒子群类,wk表示第k类的概率,μk为第k类的均值,μT表示图像灰度平均值;
③.对②中所计算的每个粒子进行适应度值评估,将种群内粒子的最优个体值Pd设置为当前位置,种群最优个体值Pg设置为初始群体中最佳粒子的位置;
④.对粒子速度进行更新,如下:
上式中t代表当前迭代次数;c1、c2为学习因子,且满足3<c1+c2<4;r1、r2为[0,1]内均匀分布的随机数;系数取值根据下面的二项式来确定:
⑤.对粒子速度进行约束:
⑥.考虑到过去粒子对当前粒子的影响,对粒子进行速度更新:
Vw(t-3)=Vw(t-2)
Vw(t-2)=Vw(t-1)
Vw(t-1)=Vw(t)
⑦.对粒子个体位置进行更新:
Xx(t+1)=Xx(t)+Vw(t+1)
⑧.对更新后粒子的适应度值进行计算,并与上次迭代最优适应度值比较,判断当前种群所有个体对应的适应度值的最小值是否小于上一次的种群最优个体适应度值,
如果是,则替换上一次种群最优个体值,接着判断是否达到最大迭代次数,如果已经达到最大迭代次数,则以当前种群所有个体对应的适应度值的最小值作为最优种群个体位置结果值,得到分割后图像fd;如果未达到最大迭代次数,返回到第②步继续执行;
如果否,则以当前种群所有个体对应的适应度值的最小值作为最优种群个体位置结果值,得到分割后图像fd。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010810993.4A CN111951291B (zh) | 2020-08-13 | 2020-08-13 | 基于多结构形态学与fodpso混合处理红外图像边缘检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010810993.4A CN111951291B (zh) | 2020-08-13 | 2020-08-13 | 基于多结构形态学与fodpso混合处理红外图像边缘检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111951291A CN111951291A (zh) | 2020-11-17 |
CN111951291B true CN111951291B (zh) | 2024-02-06 |
Family
ID=73333332
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010810993.4A Active CN111951291B (zh) | 2020-08-13 | 2020-08-13 | 基于多结构形态学与fodpso混合处理红外图像边缘检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111951291B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113393387A (zh) * | 2021-05-26 | 2021-09-14 | 房超 | 一种数学形态学自适应结构元素构造方法 |
CN117541579B (zh) * | 2024-01-08 | 2024-04-26 | 江西省兆驰光电有限公司 | 一种封装支架缺陷检测方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106296709A (zh) * | 2016-08-18 | 2017-01-04 | 北京航空航天大学 | 一种基于粒子群及模糊均值聚类的细胞图像分割方法 |
WO2018072351A1 (zh) * | 2016-10-20 | 2018-04-26 | 北京工业大学 | 一种基于粒子群优化算法对支持向量机的优化方法 |
CN109919880A (zh) * | 2019-03-18 | 2019-06-21 | 郑州轻工业学院 | 一种基于粒子群优化的红外图像增强方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8417012B2 (en) * | 2008-11-04 | 2013-04-09 | Beckman Coulter, Inc. | Non-linear histogram segmentation for particle analysis |
-
2020
- 2020-08-13 CN CN202010810993.4A patent/CN111951291B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106296709A (zh) * | 2016-08-18 | 2017-01-04 | 北京航空航天大学 | 一种基于粒子群及模糊均值聚类的细胞图像分割方法 |
WO2018072351A1 (zh) * | 2016-10-20 | 2018-04-26 | 北京工业大学 | 一种基于粒子群优化算法对支持向量机的优化方法 |
CN109919880A (zh) * | 2019-03-18 | 2019-06-21 | 郑州轻工业学院 | 一种基于粒子群优化的红外图像增强方法 |
Non-Patent Citations (1)
Title |
---|
基于多尺度多方向结构元素的形态学图像边缘检测算法;吴朔媚;韩明;王敬涛;;量子电子学报(03);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111951291A (zh) | 2020-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | Survey: State of the art in NDE data fusion techniques | |
CN108921201B (zh) | 基于特征组合与cnn的大坝缺陷识别与分类方法 | |
CN110210463A (zh) | 基于Precise ROI-Faster R-CNN的雷达目标图像检测方法 | |
US8413493B1 (en) | Method for detecting a defect on an operating turbine rotor blade | |
CN111507990A (zh) | 一种基于深度学习的隧道表面缺陷分割方法 | |
CN111951291B (zh) | 基于多结构形态学与fodpso混合处理红外图像边缘检测方法 | |
CN108596930B (zh) | 基于天牛须算法和数学形态学的射线图像焊缝提取方法 | |
CN112819775B (zh) | 一种航空航天复合材料损伤检测图像的分割强化方法 | |
Xu et al. | An infrared thermal image processing framework based on superpixel algorithm to detect cracks on metal surface | |
CN110287898A (zh) | 一种光学卫星遥感影像云检测方法 | |
CN107992818B (zh) | 一种光学遥感图像海面船只目标的检测方法 | |
CN113838054A (zh) | 基于人工智能的机械零件表面损伤检测方法 | |
CN112330538A (zh) | 一种基于特征点优化提取的损伤温度重构图像拼接方法 | |
CN110673138A (zh) | 一种基于奇异值分解和模糊c均值法的探地雷达图像处理方法 | |
CN117853722A (zh) | 一种融合超像素信息的钢材金相组织分割方法 | |
CN116823725A (zh) | 一种基于支持向量机的航空发动机叶片表面缺陷检测方法 | |
CN115861226A (zh) | 基于特征值梯度变化的深度神经网络智能识别表面缺陷方法 | |
CN116310355A (zh) | 一种复杂结构件的激光点云去噪和缺陷检测方法 | |
Fan et al. | Application of YOLOv5 neural network based on improved attention mechanism in recognition of Thangka image defects | |
CN112581548B (zh) | 一种星敏感器伪星目标滤除方法及系统 | |
CN112396580B (zh) | 一种圆形零件缺陷检测方法 | |
CN107346549B (zh) | 一种利用遥感影像多特征的多类别变化动态阈值检测方法 | |
CN113781445A (zh) | 一种多区域复杂损伤缺陷特征提取融合方法 | |
CN113706588A (zh) | 基于改进四点快速鲁棒匹配算法的环形锻件点云配准方法 | |
CN114972353B (zh) | 一种金属高压特种管件质量评估方法及系统 |
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 |