CN102930273B - 基于亮度自适应水平集的极光卵分割方法 - Google Patents
基于亮度自适应水平集的极光卵分割方法 Download PDFInfo
- Publication number
- CN102930273B CN102930273B CN201210388695.6A CN201210388695A CN102930273B CN 102930273 B CN102930273 B CN 102930273B CN 201210388695 A CN201210388695 A CN 201210388695A CN 102930273 B CN102930273 B CN 102930273B
- Authority
- CN
- China
- Prior art keywords
- mrow
- image
- aurora
- level set
- phi
- 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
- 238000000034 method Methods 0.000 title claims abstract description 59
- 239000005441 aurora Substances 0.000 claims abstract description 156
- 230000011218 segmentation Effects 0.000 claims abstract description 72
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 230000000877 morphologic effect Effects 0.000 claims description 50
- 210000004681 ovum Anatomy 0.000 claims description 46
- 102000002322 Egg Proteins Human genes 0.000 claims description 41
- 108010000912 Egg Proteins Proteins 0.000 claims description 41
- 230000003044 adaptive effect Effects 0.000 claims description 7
- 239000002131 composite material Substances 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 238000007781 pre-processing Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 230000006978 adaptation Effects 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 229920007962 Styrene Methyl Methacrylate Polymers 0.000 claims 1
- 238000013507 mapping Methods 0.000 claims 1
- 230000007547 defect Effects 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 30
- 238000002474 experimental method Methods 0.000 description 7
- 238000003709 image segmentation Methods 0.000 description 7
- 238000001514 detection method Methods 0.000 description 4
- 239000002245 particle Substances 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241001676573 Minium Species 0.000 description 1
- 244000078856 Prunus padus Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000003064 k means clustering Methods 0.000 description 1
- 210000002569 neuron Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002243 precursor Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Landscapes
- Image Analysis (AREA)
Abstract
本发明公开了一种基于亮度自适应水平集的极光卵分割方法,主要解决现有极光卵分割方法分割精度较低、鲁棒性较差和适用面较窄的缺陷。其实现步骤为:①采用形态学成分分析法对紫外极光图像进行预处理;②构建一个形态学显著性图作为极光卵的形状特征;③应用形态学显著性图的边缘曲线对水平集函数初始化;④计算亮度自适应的水平集演化速度和停止函数;⑤根据亮度自适应的水平集演化方程更新水平集函数;⑥提取更新完成后的零水平集曲线,并将其作为极光卵边界输出。本发明避免了现有分割方法出现的结果偏移和边缘泄漏现象,具有分割精度高和鲁棒性强的优点,可适用于各类紫外极光图像的分割。
Description
技术领域
本发明属于图像处理技术领域,涉及紫外波段捕获图像的目标分割,可用于从复杂的紫外极光图像中将极光卵区域分割出来。
背景技术
极光是极地区域空间光学辐射中一种重要的自然发光现象。它是由太阳风和磁尾带电粒子体与大气中的分子、原子在近地空间相互作用激发而成的。它的发生和发展与太阳活动周期、地磁活动强度密切相关。紫外波段成像探测由于可以捕捉到任何其它波段都难以观测到的极其微弱的能量散射,进而对太阳进入地球极区的能量波动进行全天候、全区域的探测,从而成为极光探测的主要手段。紫外极光图像中的极光卵区域是极光粒子以磁极为中心沉降在地球南北极地区形成的椭圆带状区域,它的赤道向边界和极向边界是重要的地球物理参数,因此,有效提取紫外极光图像中的极光卵边界对于研究极光粒子沉降规律,进而理解太阳风-磁层能量耦合过程起着十分重要的作用。
由于紫外极光图像的强噪声、低对比度以及众多的混淆现象,如日晖、彗星、宇宙射线等,所以将极光卵区域从图像中分割出来具有很大的挑战性。目前,对紫外极光图像极光卵区域进行分割的方法主要可以分为两类:无形状知识的方法和有形状知识的方法。前者不需要任何形状先验知识,仅考虑图像自身的灰度和像素间的空间关系信息。例如,Hung等人在文献“C.Hung,G.Germany,K-means anditerative selection algorithms in image segmentation,in Proc.IEEE Southeastcon,2003,(Session 1:Software Development).”中提出的基于直方图的k-均值法(Histogram-basedK-means,HKM),该方法对极光图像的灰度直方图进行k-均值聚类后,将高于灰度阈值的几个聚类作为极光卵区域;Li等人在文献“X.Li,R.Ramachandran,M.He,S.Movva,J.A.Rushing,S.J.Graves,W.Lyatsky,A.Tan,G.A.Germany,Comparingdifferent thresholding algorithms for segmenting auroras,in Proc.IEEE Conf.InformationTechnology:Coding and Computing,2004,pp.594-601.”中提出的自适应最小误差阈值法(Adaptive Minimum Error Thresholding,AMET),该方法根据磁地方时(MagneticLocal Time,MLT)将紫外极光图像划分成若干个子区域,然后对每个子区域应用最小误差阈值法进行分割处理;Germany等人在文献“G.A.Germany,G.K.Parks,H.Ranganath,R.Elsen,P.G.Richards,W.Swift,J.P.Spann,M.Brittnacher,Analysis ofauroral morphology:substorm precursor and onset on January 10,1997,GeophysicalResearch Letters 25(15)(1998)3043-3046.”中提出的脉冲耦合神经网络法(Pulse-coupled Neural Network,PCNN),该方法将图像的每个像素点看作一个分割神经元,利用相邻像素之间的灰度关系,采用神经网络法对极光卵区域进行分割。此类无形状知识的极光卵分割方法对于目标和背景灰度相差较大的极光图像分割效果较好,但是,在众多低对比度图像中会出现噪点大、误分率和漏分率较高等问题。而有形状知识的极光卵分割方法则结合了极光卵的形状特征先验知识,例如Cao等人在文献“C.Cao,T.S.Newman,G.Germany,A new shape-based auroral oval segmentationdrivenby LLS-RHT,Pattern Recognition 42(2009)607-618.”中提出的基于椭圆形状先验知识的线性最小二乘-随机霍夫变换法(Linear least-squares(LLS)-randomized Houghtransform(RHT),LLS-RHT),该方法结合极光图像灰度信息和形状特征,用一个类似椭圆的图形来拟合极光卵区域,虽然能够对形态为完整型椭圆的极光卵进行较好的分割,但是,对于形态是缺口状椭圆的极光卵,分割结果会出现很大的偏差。
水平集方法是近年来图像分割领域研究的热点之一。其中,Li等人在文献“C.Li,C.Xu,C.Gui,M.D.Fox,Level set evolution without re-initialization:A newvariational formulation,in Proc.IEEE Conf.Comput.Vis.Pattern Recog.,2005,pp.430-436.”中提出的不需要初始化的水平集方法,由于在保证了高分割准确率的同时,较高地提升了分割速度而广泛应用在医学图像分割等领域。然而,Li等人的方法对曲线的初始化是随机的,而极光图像的背景和目标灰度相差较小,若初始曲线离目标边缘较远,分割结果就会出现很大的偏离。同时,该方法中的曲线演化速度和停止函数与曲线周围的灰度值没有相关性,因而在处理极光卵分割这类弱边缘问题时,即使选用了离目标边缘很近的初始化曲线,也会出现较大的边缘泄漏现象。
综上所述,现有的极光卵分割方法和目前常用的水平集图像分割方法对紫外极光图像都不能达到满意的分割效果。
发明内容
本发明的目的在于针对上述极光卵分割方法的不足,提出一种基于亮度自适应水平集的极光卵分割方法,以提高极光卵区域分割的精度和准确率,为空间物理研究提供尽可能精确的极光卵参数信息。
实现本发明的技术思路是:通过采用形态学成分分析法对紫外极光图像进行预处理,达到降噪的目的;通过构建一个形态学显著性图,作为极光卵的形状特征,并以形态学显著性图的边缘曲线作为水平集演化的初始曲线,采用亮度自适应的水平集演化方程对曲线进行演化,使曲线最终逼近目标的边缘,从而实现极光卵的分割,具体实现步骤包括如下:
(1)输入一幅紫外极光图像I,运用形态学成分分析的方法将紫外极光图像I分解为极光形态学子层IA和噪声形态学子层IN,去掉噪声形态学子层IN,将极光形态学子层IA输出,得到预处理后的图像IP;
(2)根据磁地方时将预处理后的图像IA分为24个子区域,对每个子区域进行基于最小平方误差准则的阈值选取,根据每个子区域不同的阈值对预处理后的图像IA进行二值化操作,得到二值化后的图像IB;
(3)将二值化后的图像IB通过一组结构元素逐渐变大的复合形态学滤波器F,进行噪声点去除,得到滤波后的图像IF;
(4)寻找滤波后的图像IF中的最大连通域,找出最大连通域外边界曲线围成的图像区域,计算该图像区域内值为0的像素点个数N,设定经验阈值TA=1200,判断紫外极光图像的类型,若N大于等于TA,则该图像属于完整型椭圆极光卵图像,若N小于TA,则属于缺口状椭圆极光卵图像;
(5)根据紫外极光图像类型计算其形态学显著性图SMM,若图像属于完整型椭圆极光卵图像,则对滤波后的图像IF进行椭圆拟合,得到椭圆拟合后的图像IE,若图像属于缺口状椭圆极光卵图像,则对滤波后的图像IF不做处理,即:
(6)将图像的形态学显著性图SMM的边界作为水平集演化的初始曲线,初始化水平集函数:
其中,φ0为水平集函数的初始值,d为正的常数,Ω表示整个紫外极光图像区域,(x,y)表示紫外极光图像中第x行第y列的像素点,若该像素点属于形态学显著性图SMM以内的区域,则给φ0赋值为d,若该像素点在SMM的边界上,则给φ0赋值为0,若该像素点属于SMM以外的区域,则给φ0赋值为-d;
(7)计算像素点(x,y)及其八邻域的平均灰度值,将其归一化到[0,1]区间作为归一化的局部灰度平均值Il(x,y),计算亮度自适应的水平集演化速度v(x,y)和停止函数g(x,y):
其中,q为控制演化速度大小的常数,p为控制演化速度随Il(x,y)非线性变化程度的常数,Gλ为标准差为λ的高斯核,*表示卷积运算,表示梯度运算,||表示取模运算,m为常数;
(8)设置水平集函数更新停止条件的阈值ε=0.01,根据亮度自适应的水平集演化方程更新水平集函数:
(8a)输入t时刻的水平集函数φt,根据亮度自适应的水平集演化方程,计算水平集函数的变化值
其中,φ为水平集函数,表示水平集函数随时间的变化值,u和β为控制水平集演化曲线的形状如何变化的常数,Δ为拉普拉斯算子,div表示散度运算,δ(·)为单变量脉冲函数(Dirac函数),v和g分别为亮度自适应的水平集演化速度和停止函数;
(8b)计算更新后的水平集函数φt+1:
其中,Δt为时间步长;
(8c)判断水平集函数更新前后的数值差异是否满足水平集函数更新的停止条件:
|φt+1-φt|≤ε,
若满足,则停止水平集函数的更新,输出水平集函数,若不满足,则继续进行水平集函数的更新;
(9)提取更新完成后的水平集函数中值为0的像素点,将其作为极光卵区域的边界,输出极光卵区域。
本发明与现有的极光卵分割方法和目前常用的水平集图像分割方法相比,具有以下优点:
1)本发明通过构建极光图像的形态学显著性图,获得了良好的极光卵形状特征,将该形态学显著性图的边缘作为水平集算法的初始演化曲线,能够实现极光卵区域的精确分割,避免了目前常用水平集分割方法出现的结果偏移现象;
2)本发明结合极光图像的局部灰度信息,设计了亮度自适应的水平集演化速度和停止函数,能够实现极光卵弱边缘的良好分割,避免了目前常用水平集分割方法出现的边缘泄漏现象;
3)本发明与现有的极光卵分割方法相比,不但能够实现目标与背景灰度相差较大的极光图像的精确分割,而且对目标和背景灰度相差不大的极光图像也能获得良好的分割效果,有效降低了误分率和漏分率,并且对完整型椭圆和缺口状椭圆的极光卵图像同时适用,扩大了极光卵分割方法的适用面。
本发明的技术过程和效果可结合以下附图详细说明。
附图说明
图1是本发明的流程图;
图2是本发明中基于磁地方时区域划分的示意图;
图3是现有完整型椭圆极光卵原图像;
图4是本发明与现有五种分割方法对图3的分割效果图;
图5是现有缺口状椭圆极光卵原图像;
图6是本发明与现有五种分割方法对图5的分割效果图。
具体实施方式
参照图1,本发明基于亮度自适应水平集的极光卵分割方法主要包括三个阶段:图像预处理、极光卵形状特征提取、亮度自适应水平集演化,具体步骤如下:
步骤1:对输入的紫外极光图像进行预处理:
1.1)构建紫外极光图像数据库,即将原始紫外极光频道定义CDF数据转化为标签图像文件TIFF数据,每个图像文件为大小是200×228的灰度图;
1.2)输入一幅紫外极光图像I,并对其进行预处理:
1.2.1)构建形态学成分分析方法的稀疏表示字典ΦA和ΦN;
1.2.2)设定形态学成分分析方法的条件为:
其中,为目标形态学子层,为背景形态学子层;
1.2.3)将目标形态学子层和背景形态学子层作为变量,代入形态学成分分析方法的目标函数:
即求解在稀疏表示字典ΦA和ΦN下,目标形态学子层和背景形态学子层的最优表达,其中,|| ||1表示1范数运算;
1.2.4)输出目标形态学子层和背景形态学子层的最优表达,将目标形态学子层的最优表达作为极光形态学子层IA,将背景形态学子层的最优表达作为噪声形态学子层IN;
1.2.5)去掉噪声形态学子层IN,将极光形态学子层IA输出,得到预处理后的图像IP。
步骤2:提取紫外极光图像中极光卵的形状特征:
2.1)对预处理后的图像IP进行基于磁地方时自适应门限的二值化操作:
2.1.1)参照图2,根据磁地方时将预处理后的图像IP分为24个子区域,每个子区域跨越2小时的磁地方时,相邻的两个子区域有1小时磁地方时的重叠;
2.1.2)选取24个子区域中的一块,假设它的目标和背景服从均值为μk,标准差为σk的正态分布,其中k=1,2,当k=1时表示目标,当k=2时表示背景;
2.1.3)根据均值μk和标准差σk估计所选子区域的直方图分布He(i,t):
其中,i表示该子区域中某一点像素的灰度值,T表示分割目标和背景的灰度阈值门限,Pk为目标和背景所占该子区域面积的比例;
2.1.4)计算所选子区域实际的直方图分布H(i);
2.1.5)根据最小平方误差准则选取最优阈值门限Tm:
其中,(He(i,T)-H(i))为估计值对实际值的误差,Tm为能够使该误差的平方最小的阈值门限T的取值;
2.1.6)对所选子区域进行二值化操作,即将灰度值大于等于Tm的像素点当作目标并设定取值为1,将灰度值小于Tm的像素点当作背景并设定取值为0;
2.1.7)对24个子区域都进行二值化操作后,得到二值化后的图像IB;
2.2)将二值化后的图像IB进行多结构元素复合形态学滤波:
2.2.1)选择一组形状为圆形,大小逐步增大的结构元素{B}={B1,B2,...,Bn},Br为这组结构元素中的一个,r=1,2,…n,n为这组结构元素的个数;
2.2.2)构建多结构元素复合形态学滤波器F:
F=(((AοBj)·Bj)οBj+1)·Bj+1… j=1,...,n,
其中,A为将要通过滤波器F的图像,j为结构元素的序号,取值从1到n,ο表示形态学开操作,·表示形态学闭操作;
2.2.3)通过滤波器F对二值化后的图像IB进行处理,得到输出的图像IF:
IF=IB*F+IR*F,
其中,IR=(IB-IB*F),其为二值化后的图像IB与滤波器F卷积后的剩余部分;
2.3)寻找滤波后的图像IF中的最大连通域,找出最大连通域外边界曲线围成的图像区域,计算该图像区域内值为0的像素点个数N,设定经验阈值TA=1200,判断紫外极光图像的类型,若N大于等于TA,则该图像属于完整型椭圆极光卵图像,若N小于TA,则属于缺口状椭圆极光卵图像;
2.4)若紫外极光图像属于完整型椭圆极光卵图像,则对滤波后的图像IF进行椭圆拟合:
2.4.1)对滤波后的图像IF随机采样l个点{p(x,y)l},l=1,2,…L,L≥5,其中(x,y)为采样点的坐标值;
2.4.2)用随机采样点{p(x,y)l}拟合二次方程式:
A*x2+B*y2+C*xy+D*x+E*y+F*=0,
其中,A*为x的二次项系数,B*为y的二次项系数,C*为x和y的乘积项系数,D*为x的一次项系数,E*为y的一次项系数,F*为常数项;
2.4.3)判断拟合的二次方程式是否满足如下条件:
若满足,则进行下一步,若不满足,则重新进行步骤2.4.1);
2.4.4)将二次方程式变形为椭圆的标准方程式:
其中,(xc,yc)是椭圆中心,a和b为椭圆长短轴,θ为椭圆方向角;
2.4.5)根据椭圆中心(xc,yc)、椭圆长短轴a和b以及椭圆方向角θ,确定椭圆拟合后的图像IE;
2.5)若紫外极光图像属于缺口状椭圆极光卵图像,则对滤波后的图像IF不做处理;
2.6)将紫外极光图像的形态学显著性图SMM作为极光卵的形状特征:
步骤3:对紫外极光图像进行基于亮度自适应的水平集演化:
3.1)初始化水平集函数,即将紫外极光图像的形态学显著性图SMM的边界作为水平集演化的初始曲线,初始化水平集函数:
其中,φ0为水平集函数的初始值,d为正的常数,Ω表示整个紫外极光图像区域,(x,y)表示紫外极光图像中第x行第y列的像素点,若该像素点属于形态学显著性图SMM以内的区域,则给φ0赋值为d,若该像素点在SMM的边界上,则给φ0赋值为0,若该像素点属于SMM以外的区域,则给φ0赋值为-d;
3.2)计算亮度自适应的水平集演化速度和停止函数:
3.2.1)计算像素点(x,y)及其八邻域的平均灰度值,将其归一化到[0,1]区间作为归一化的局部灰度平均值Il(x,y),
3.2.2)计算亮度自适应的水平集演化速度v(x,y):
其中,q为控制演化速度大小的常数,p为控制演化速度随Il(x,y)非线性变化程度的常数;
3.2.3)计算亮度自适应的水平集停止函数g(x,y):
其中,Gλ为标准差为λ的高斯核,*表示卷积运算,表示梯度运算,||表示取模运算,m为常数;
3.3)对水平集函数进行更新:
3.3.1)设定水平集函数更新停止条件的阈值ε=0.01;
3.3.2)输入t时刻的水平集函数φt;
3.3.3)根据亮度自适应的水平集演化方程,计算水平集函数的变化值
其中,φ为水平集函数,表示水平集函数随时间的变化值,u和β为控制水平集演化曲线形状变化的常数,Δ为拉普拉斯算子,div表示散度运算,δ(·)为单变量脉冲函数(Dirac函数),v和g分别为亮度自适应的水平集演化速度和停止函数;
3.3.4)计算更新后的水平集函数φt+1:
其中,Δt为时间步长;
3.3.5)判断水平集函数更新前后的数值差异是否满足水平集函数更新的停止条件:
|φt+1-φt|≤ε,
若满足,则停止水平集函数的更新,输出水平集函数,若不满足,则继续进行水平集函数的更新;
3.4)提取更新完成后的水平集函数中值为0的像素点,将这些像素点作为极光卵区域的边界并输出。
本发明的效果可通过以下仿真实验进一步说明:
实验1:极光卵分割效果比较
实验条件:实验用到的完整型椭圆极光卵原图像如图3所示,其中,图3(a)为原始极光卵图像,图3(b)为极光专家对图3(a)进行手工标定的原图像;实验用到的缺口状椭圆极光卵原图像如图5所示,其中,图5(a)为原始极光卵图像,图5(b)为极光专家对图5(a)进行手工标定的原图像。
实验内容:本实验将比较本发明和现有五种分割方法对紫外极光图像极光卵区域分割的效果,现有五种分割方法为背景技术中提及的HKM、AMET、PCNN、LLS-RHT以及Li等人的水平集方法。实验分成两部分:第一部分为完整型椭圆极光卵图像的分割效果比较,第二部分为缺口状椭圆极光卵图像的分割效果比较。
1.完整型椭圆极光卵图像的分割效果比较:图4是本发明与现有五种分割方法对图3的分割效果图,其中,图4(a)为HKM对图3的分割效果图,图4(b)为AMET对图3的分割效果图,图4(c)为PCNN对图3的分割效果图,图4(d)为LLS-RHT对图3的分割效果图,图4(e)为Li等人的水平集方法对图3的分割效果图,图4(f)为本发明对图3的分割效果图。
从图4中可以看出,HKM、AMET、PCNN对完整型椭圆极光卵图像的分割结果含有较多的噪点,并且误分率和漏分率较高;LLS-RHT的分割结果大体符合完整型椭圆极光卵的形态,但是却没有表现出极光卵内边界的凹凸起伏;Li等人的水平集方法出现了边缘泄漏,分割不完整的现象;而本发明的分割结果则比较精确。
2.缺口状椭圆极光卵图像的分割效果比较:图6是本发明与现有五种分割方法对图5的分割效果图,其中,图6(a)为HKM对图5的分割效果图,图6(b)为AMET对图5的分割效果图,图6(c)为PCNN对图5的分割效果图,图6(d)为LLS-RHT对图5的分割效果图,图6(e)为Li等人的水平集方法对图5的分割效果图,图6(f)为本发明对图5的分割效果图。
从图6中可以看出,HKM、AMET、PCNN对缺口状椭圆极光卵图像的分割结果存在与完整型椭圆极光卵图像相同的问题;LLS-RHT的分割结果出现了很大的偏移现象;Li等人的水平集方法分割结果大体良好,但仍旧没有本发明的结果精确。
实验2:极光卵分割误差比较
实验条件:本实验用到了30幅完整型椭圆极光卵原图像和30幅缺口状椭圆极光卵原图像。
实验内容:比较本发明与现有五种分割方法对极光卵原图像中极光卵区域进行分割的总分割误差,标准边界由极光专家手工标定。
对极光卵原图像中极光卵区域进行分割的总分割误差SE计算如下:
SE=α*·FS+β*·MS+γ*·BDE,
其中,FS为误分率,表示将背景像素点当作极光卵的误差,MS为漏分率,表示将极光卵像素点当作背景的误差,BDE为边界偏离误差,表示分割结果的极光卵边界与标准边界之间的偏移距离总和,α*为误分率FS的权值比重,β*为漏分率MS的权值比重,γ*为边界偏离误差BDE的权值比重。
误分率FS、漏分率MS和边界偏离误差BDE根据下式进行计算:
其中,R为面积是A′的极光卵区域,该区域像素点取值为1,其余背景像素点取值为0,N′为分割结果图中极光卵边界像素点的个数,d2(i*)表示分割结果的第i*(i*=1,2,…N′)个边界像素点与标准边界像素点的距离,r(x,y)表示极光卵分割结果图中第x行第y列的像素点:
表6为本发明与现有五种分割方法对30幅完整型椭圆极光卵原图像的总分割误差比较。
表6.完整型椭圆极光卵图像分割结果的总分割误差比较
从表6可以看出,本发明的总分割误差的平均值最小,说明本发明对完整型椭圆极光卵图像的分割最精确,同时,本发明的总分割误差的均方差也最小,说明本发明对完整型椭圆极光卵图像的分割最鲁棒。
表7为本发明与现有五种分割方法对30幅缺口状椭圆极光卵原图像的总分割误差比较。
表7.缺口状椭圆极光卵图像分割结果的总分割误差比较
从表7可以看出,本发明的总分割误差的平均值最小,说明本发明对缺口状椭圆极光卵图像的分割最精确,同时,本发明的总分割误差的均方差也最小,说明本发明对缺口状椭圆极光卵图像的分割最鲁棒。
综上,本发明方法对紫外极光图像进行基于亮度自适应水平集的极光卵分割,提高了极光卵分割的准确性,增强了分割效果的鲁棒性,扩大了对不同形态极光卵图像的适用性。
Claims (4)
1.一种基于亮度自适应水平集的极光卵分割方法,包括如下步骤:
(1)输入一幅紫外极光图像I,运用形态学成分分析的方法将紫外极光图像I分解为极光形态学子层IA和噪声形态学子层IN,去掉噪声形态学子层IN,将极光形态学子层IA输出,得到预处理后的图像IP;
(2)根据磁地方时将预处理后的图像IP分为24个子区域,对每个子区域进行基于最小平方误差准则的阈值选取,根据每个子区域不同的阈值对预处理后的图像IP进行二值化操作,得到二值化后的图像IB;
(3)将二值化后的图像IB通过一组结构元素逐渐变大的复合形态学滤波器F,进行噪声点去除,得到滤波后的图像IF;
(4)寻找滤波后的图像IF中的最大连通域,找出最大连通域外边界曲线围成的图像区域,计算该图像区域内值为0的像素点个数N,设定经验阈值TA=1200,判断紫外极光图像的类型,若N大于等于TA,则该图像属于完整型椭圆极光卵图像,若N小于TA,则属于缺口状椭圆极光卵图像;
(5)根据紫外极光图像类型计算其形态学显著性图SMM,若图像属于完整型椭圆极光卵图像,则对滤波后的图像IF进行椭圆拟合,得到椭圆拟合后的图像IE,若图像属于缺口状椭圆极光卵图像,则对滤波后的图像IF不做处理,即:
(6)将紫外极光图像的形态学显著性图SMM的边界作为水平集演化的初始曲线,初始化水平集函数:
其中,φ0为水平集函数的初始值,d为正的常数,Ω表示整个紫外极光图像区域,(x,y)表示紫外极光图像中第x行第y列的像素点,若该像素点属于形态学显著性图SMM以内的区域,则给φ0赋值为d,若该像素点在SMM的边界上,则给φ0赋值为0,若该像素点属于SMM以外的区域,则给φ0赋值为-d;
(7)计算像素点(x,y)及其八邻域的平均灰度值,将其归一化到[0,1]区间作为归一化的局部灰度平均值Il(x,y),计算亮度自适应的水平集演化速度ν(x,y)和停止函数g(x,y):
其中,q为控制演化速度大小的常数,p为控制演化速度随Il(x,y)非线性变化程度的常数,Gλ为标准差为λ的高斯核,*表示卷积运算,表示梯度运算,||表示取模运算,m为常数;
(8)设置水平集函数更新停止条件的阈值ε=0.01,根据亮度自适应的水平集演化方程更新水平集函数:
(8a)输入t时刻的水平集函数φt,根据亮度自适应的水平集演化方程,计算水平集函数的变化值
其中,φ为水平集函数,表示水平集函数随时间的变化值,u和β为控制水平集演化曲线的形状如何变化的常数,Δ为拉普拉斯算子,div表示散度运算,δ(·)为单变量脉冲Dirac函数,ν和g分别为亮度自适应的水平集演化速度和停止函数;
(8b)计算更新后的水平集函数φt+1:
其中,Δt为时间步长;
(8c)判断水平集函数更新前后的数值差异是否满足水平集函数更新的停止条件:
|φt+1-φt|≤ε,
若满足,则停止水平集函数的更新,输出水平集函数,若不满足,则继续进行水平集函数的更新;
(9)提取更新完成后的水平集函数中值为0的像素点,将其作为极光卵区域的边界,输出极光卵区域。
2.根据权利要求1所述的基于亮度自适应水平集的极光卵分割方法,其中步骤(1)所述的运用形态学成分分析的方法将紫外极光图像I分解为极光形态学子层IA和噪声形态学子层IN,按如下公式进行:
其中,ΦA为极光形态学子层的稀疏表示字典,ΦN为噪声形态学子层的稀疏表示字典,为目标形态学子层,为背景形态学子层,||||1表示1范数运算,形态学成分分析法就是在目标形态学子层和背景形态学子层之和为输入图像I的条件下,寻找稀疏表示字典ΦA和ΦN下,目标形态学子层和背景形态学子层的最优表达,并将目标形态学子层的最优表达作为极光形态学子层IA,将背景形态学子层的最优表达作为噪声形态学子层IN。
3.根据权利要求1所述的基于亮度自适应水平集的极光卵分割方法,其中步骤(2)所述的对每个子区域进行基于最小平方误差准则的阈值选取,按如下步骤进行:
(3a)选取24个子区域中的一块,假设它的目标和背景服从均值为μk,标准差为σk的正态分布,其中k=1,2,当k=1时表示目标,当k=2时表示背景;
(3b)估计所选子区域的直方图分布He(i,T):
其中,i表示该子区域中某一点像素的灰度值,T表示分割目标和背景的灰度阈值门限,Pk为目标和背景所占该子区域面积的比例;
(3c)计算所选子区域实际的直方图分布H(i);
(3d)根据最小平方误差准则选取最优阈值门限Tm:
其中,(He(i,T)-H(i))为估计值对实际值的误差,则Tm为能够使该误差的平方最小的阈值门限T的取值。
4.根据权利要求1所述的基于亮度自适应水平集的极光卵分割方法,其中步骤(3)所述的将二值化后的图像IB通过一组结构元素逐渐变大的复合形态学滤波器F进行噪声点去除,按如下步骤进行:
(4a)选择一组形状为圆形,大小逐步增大的结构元素{B}={B1,B2,...,Bn},Br为这组结构元素中的一个,r=1,2,…n,n为这组结构元素的个数;
(4b)构建复合形态学滤波器F:
F=(((AоBj)·Bj)оBj+1)·Bj+1…j=1,...,n,
其中,A为将要通过滤波器F的图像,j为结构元素的序号,取值从1到n,ο表示形态学开操作,·表示形态学闭操作;
(4c)通过滤波器F对二值化后的图像IB进行处理,得到输出的图像IF:
IF=IB*F+IR*F,
其中,IR=(IB-IB*F),为二值化后的图像IB与滤波器F卷积后的剩余部分。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210388695.6A CN102930273B (zh) | 2012-10-15 | 2012-10-15 | 基于亮度自适应水平集的极光卵分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210388695.6A CN102930273B (zh) | 2012-10-15 | 2012-10-15 | 基于亮度自适应水平集的极光卵分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102930273A CN102930273A (zh) | 2013-02-13 |
CN102930273B true CN102930273B (zh) | 2015-06-17 |
Family
ID=47645070
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210388695.6A Expired - Fee Related CN102930273B (zh) | 2012-10-15 | 2012-10-15 | 基于亮度自适应水平集的极光卵分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102930273B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103632166A (zh) * | 2013-12-04 | 2014-03-12 | 西安电子科技大学 | 基于融合显著信息的潜在主题极光图像分类方法 |
CN103617425B (zh) * | 2013-12-11 | 2016-08-17 | 西安邮电大学 | 用于监测极光活动的变化轨迹的生成方法 |
CN103971120B (zh) * | 2014-04-18 | 2017-02-15 | 西安电子科技大学 | 基于空时极向局部二值模式的极光图像序列分类方法 |
CN104766085B (zh) * | 2014-12-30 | 2019-01-29 | 沈阳理工大学 | 一种多尺度图形识别方法 |
CN104680167B (zh) * | 2015-03-09 | 2018-03-20 | 西安电子科技大学 | 基于深度学习的极光卵位置确定方法 |
CN105118047A (zh) * | 2015-07-15 | 2015-12-02 | 陕西师范大学 | 基于行星际和地磁参数的极光卵边界位置预测方法 |
CN105069798A (zh) * | 2015-08-12 | 2015-11-18 | 西安电子科技大学 | 基于图割优化的局部主动轮廓图像分割方法 |
CN105321184B (zh) * | 2015-11-26 | 2018-01-16 | 北京交通大学 | 改进的基于边缘水平集的含噪图像分割方法与*** |
CN106570867B (zh) * | 2016-10-18 | 2019-03-29 | 浙江大学 | 基于灰度形态学能量法的活动轮廓模型图像快速分割方法 |
CN107045722B (zh) * | 2017-03-27 | 2019-07-30 | 西安电子科技大学 | 融合静态信息与动态信息的视频序列分割方法 |
CN110176021B (zh) * | 2019-05-21 | 2021-04-16 | 山东大学 | 结合亮度校正的显著性信息的水平集图像分割方法及*** |
CN111325696B (zh) * | 2020-03-03 | 2023-07-07 | 杭州瑞利海洋装备有限公司 | 一种基于正态分布区间估计的水声图像混响抑制方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101504767A (zh) * | 2009-03-20 | 2009-08-12 | 西安电子科技大学 | 基于水平集接力的图像分割方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8515171B2 (en) * | 2009-01-09 | 2013-08-20 | Rochester Institute Of Technology | Methods for adaptive and progressive gradient-based multi-resolution color image segmentation and systems thereof |
-
2012
- 2012-10-15 CN CN201210388695.6A patent/CN102930273B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101504767A (zh) * | 2009-03-20 | 2009-08-12 | 西安电子科技大学 | 基于水平集接力的图像分割方法 |
Non-Patent Citations (4)
Title |
---|
Comparing Different Thresholding Algorithms for Segmenting Auroras;Xiang Li et al.;《Proceedings of the International Conference on Information Technology: Coding and Computing》;20041231;1-8 * |
Distance Regularized Level Set Evolution and Its Application to Image Segmentation;Chunming Li et al.;《IEEE TRANSACTIONS ON IMAGE PROCESSING》;20101231;第19卷(第12期);3243-3254 * |
Level Set Evolution Without Re-initialization: A New Variational Formulation;Chunming Li et al.;《Proceedings of the 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition》;20051231;1-7 * |
一种基于边缘与区域信息的先验水平集图像分割方法;王斌等;《计算机学报》;20120531;第35卷(第5期);1067-1072 * |
Also Published As
Publication number | Publication date |
---|---|
CN102930273A (zh) | 2013-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102930273B (zh) | 基于亮度自适应水平集的极光卵分割方法 | |
Zhang et al. | Abnormal breast identification by nine-layer convolutional neural network with parametric rectified linear unit and rank-based stochastic pooling | |
CN109409292B (zh) | 基于精细化特征优化提取的异源图像匹配方法 | |
Tao et al. | A deep neural network modeling framework to reduce bias in satellite precipitation products | |
Chitradevi et al. | An overview on image processing techniques | |
Storie et al. | Deep learning neural networks for land use land cover mapping | |
Quackenbush | A review of techniques for extracting linear features from imagery | |
Bhalla et al. | Fusion of Infrared and Visible Images Using Fuzzy Based Siamese Convolutional Network. | |
CN110458192B (zh) | 基于视觉显著性的高光谱遥感图像分类方法及*** | |
Liu et al. | Multimorphological superpixel model for hyperspectral image classification | |
Kurmi et al. | Multifeature‐based medical image segmentation | |
Jain et al. | Automatic flood detection in SentineI-2 images using deep convolutional neural networks | |
Fu et al. | A novel spectral-spatial singular spectrum analysis technique for near real-time in situ feature extraction in hyperspectral imaging | |
Goel et al. | Automated detection, characterization, and tracking of sunspots from SoHO/MDI continuum images | |
Kwon et al. | ETVOS: An enhanced total variation optimization segmentation approach for SAR sea-ice image segmentation | |
Chen et al. | A novel infrared small target detection method based on BEMD and local inverse entropy | |
Shahbeig | Automatic and quick blood vessels extraction algorithm in retinal images | |
Qin et al. | Multilevel wavelet-SRNet for SAR target recognition | |
Haq et al. | 3D-CNNHSR: A 3-Dimensional Convolutional Neural Network for Hyperspectral Super-Resolution. | |
Ball et al. | Level set hyperspectral image classification using best band analysis | |
CN104700427B (zh) | 一种基于sar图像的道路损毁信息提取方法 | |
CN102496154B (zh) | 一种基于马尔科夫随机场的多时相遥感图像变化检测方法 | |
Zemouri et al. | Restoration based Contourlet Transform for historical document image binarization | |
Skuratov et al. | Application of a convolutional neural network and a kohonen network for accelerated detection and recognition of objects in images | |
Song et al. | HDTFF-Net: Hierarchical deep texture features fusion network for high-resolution remote sensing scene classification |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150617 Termination date: 20201015 |
|
CF01 | Termination of patent right due to non-payment of annual fee |