CN107194896B - 一种基于邻域结构的背景抑制方法和系统 - Google Patents
一种基于邻域结构的背景抑制方法和系统 Download PDFInfo
- Publication number
- CN107194896B CN107194896B CN201710415948.7A CN201710415948A CN107194896B CN 107194896 B CN107194896 B CN 107194896B CN 201710415948 A CN201710415948 A CN 201710415948A CN 107194896 B CN107194896 B CN 107194896B
- Authority
- CN
- China
- Prior art keywords
- background
- neighborhood
- pixel
- value
- point
- 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
- 230000001629 suppression Effects 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims abstract description 45
- 239000013598 vector Substances 0.000 claims abstract description 45
- 238000005259 measurement Methods 0.000 claims abstract description 21
- 230000002401 inhibitory effect Effects 0.000 claims abstract description 6
- 238000011524 similarity measure Methods 0.000 claims description 18
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 abstract description 15
- 230000000694 effects Effects 0.000 abstract description 5
- 230000005764 inhibitory process Effects 0.000 abstract description 4
- 238000001914 filtration Methods 0.000 description 29
- 238000002474 experimental method Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 7
- 238000001514 detection method Methods 0.000 description 5
- 230000000877 morphologic effect Effects 0.000 description 4
- 230000011218 segmentation Effects 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- 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)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Image Analysis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种基于邻域结构的背景抑制方法和系统,其中方法的实现包括:基于邻域像素点的灰度值得到结构向量,计算结构向量与邻域结构模板的相似性度量,将邻域结构模板分为目标模板和背景模板,利用相似性度量得到像素点为背景点的背景似然值和像素点为目标点的目标似然值;若目标似然值与背景似然值的比值大于等于阈值,则像素点为目标点,否则,像素点为背景点;以邻域像素点集合中背景点的灰度均值作为背景估计值,若像素点为背景点,直接抑制,若像素点为目标点,将目标点的灰度值减去背景估计值进行背景抑制。本发明同时具有良好的抑制效果、对不同背景类型的自适应性以及满足实时计算要求。
Description
技术领域
本发明属于图像处理领域,更具体地,涉及一种基于邻域结构的背景抑制方法和系统。
背景技术
红外技术的快速发展,使其广泛应用于军事及民用领域,在精确制导,预警监视,光学遥感,夜间导航等领域发挥重要作用。红外小目标检测、跟踪一直以来都是国内外研究热点,远距离的弱小红外小目标极易湮没在复杂背景中而难以被探测到。
目前的小目标检测主要分空域和时域两大类,时域算法由于有多帧信息,往往可以取得较好的效果,然而,在动平台探测系统下,该类算法需要进行帧间图像配准,以消除背景运动杂波,从而导致算法复杂,计算量大,增大了探测系统的复杂性。
空域背景抑制方法可分为线性方法和非线性方法两大类,大致归为三类:模板滤波,预测滤波,形态学滤波。模板滤波有传统的均值滤波,中值滤波,去最大中值滤波,高斯滤波。这类方法计算量小,参数少,但对复杂背景抑制效果能力还不足。典型的预测滤波算法有二维最小均方差滤波TDLMS,该算法利用输出信号与期望信号均方差最小达到滤波的目的。该算法中的收敛率参数需根据不同的图像特性来设置,其收敛速度常难于满足实时要求,对强杂波背景的抑制能力也不足。top-hat滤波算法是小目标检测中常用的形态学滤波方法,形态学滤波方法需要根据背景杂波及目标的特性设计合适形态学算子和结构元素,以实现背景杂波的抑制,但图像的结构往往是不同的,而且同一幅图像中不同区域的结构也往往不一致,导致该类方法对不同背景类型的适应性不强。
由此可见,现有技术存在难以同时具有良好的抑制效果、对不同背景类型的自适应性以及实时计算要求的技术问题。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种基于邻域结构的背景抑制方法和系统,由此解决现有技术存在难以同时具有良好的抑制效果、对不同背景类型的自适应性以及实时计算要求的技术问题。
为实现上述目的,按照本发明的一个方面,提供了一种基于邻域结构的背景抑制方法,包括:
(1)选取图像中的任意一个像素点,得到以像素点为中心的邻域像素点集合,将邻域像素点集合划分为前景像素点子集和背景像素点子集,基于前景像素点子集的前景像素点灰度平均值和背景像素点子集的背景像素点灰度平均值构建像素点的邻域结构模板;
(2)基于邻域像素点的灰度值得到结构向量,计算结构向量与邻域结构模板的相似性度量,将邻域结构模板分为目标模板和背景模板,利用相似性度量得到像素点为背景点的背景似然值和像素点为目标点的目标似然值;
(3)若目标似然值与背景似然值的比值大于等于阈值,则像素点为目标点,否则,像素点为背景点;
(4)以邻域像素点集合中背景点的灰度均值作为背景估计值,若像素点为背景点,直接抑制,若像素点为目标点,将目标点的灰度值减去背景估计值进行背景抑制。
进一步的,邻域结构模板包括: 其中,f表示前景像素点灰度平均值,b表示背景像素点灰度平均值。
进一步的,步骤(2)包括:
(2-1)基于图像像素点pi,j处的邻域像素点的灰度值得到结构向量:
其中,Ii-1,j表示坐标为(i-1,j)的邻域像素点的灰度值,Ii,j-1表示坐标为(i,j-1)的邻域像素点的灰度值,Ii,j表示坐标为(i,j)的邻域像素点的灰度值,Ii+1,j表示坐标为(i+1,j)的邻域像素点的灰度值,Ii+1,j+1表示坐标为(i+1,j+1)的邻域像素点的灰度值,Ii-1,j-1表示坐标为(i-1,j-1)的邻域像素点的灰度值,Ii-1,j+1表示坐标为(i-1,j+1)的邻域像素点的灰度值,Ii+1,j-1表示坐标为(i+1,j-1)的邻域像素点的灰度值;
(2-2)计算结构向量与邻域结构模板的相似性度量pkl,
其中,d表示数据维度,为邻域结构模板的核宽矩阵;
(2-3)将邻域结构模板分为目标模板和背景模板,利用相似性度量得到像素点为背景点的背景似然值Pb和像素点为目标点的目标似然值Pf,
Pl=max{pkl|k=1,2,...,8}
其中,Pl为像素点与邻域结构模板的相似性度量,
P1f表示水平垂直方向的四邻域像素点与目标模板的相似性度量,P2f表示对角线方向的四邻域像素点与目标模板的相似性度量,则有:
P1f=max{pkl|k=1,2,...,4;l=8}
P2f=max{pkl|k=5,6,...,8;l=8}
则目标似然值Pf为:Pf=(P1f+P2f)/2。
进一步的,背景估计值为
其中,card(Sbk)为Sbk中的像素点个数,Sbk表示背景点集合,pm+i,n+j表示坐标为(m+i,n+j)的邻域像素点,Ii+m,j+n表示坐标为(m+i,n+j)的邻域像素点的灰度值。
按照本发明的另一方面,提供了一种基于邻域结构的背景抑制系统,包括:
构建邻域结构模块,用于选取图像中的任意一个像素点,得到以像素点为中心的邻域像素点集合,将邻域像素点集合划分为前景像素点子集和背景像素点子集,基于前景像素点子集的前景像素点灰度平均值和背景像素点子集的背景像素点灰度平均值构建像素点的邻域结构模板;
获取似然值模块,用于基于邻域像素点的灰度值得到结构向量,计算结构向量与结构模板的相似性度量,将邻域结构模板分为目标模板和背景模板,利用相似性度量得到像素点为背景点的背景似然值和像素点为目标点的目标似然值;
判断像素点类型模块,用于判断像素点类型,若目标似然值与背景似然值的比值大于等于阈值,则像素点为目标点,否则,像素点为背景点;
进行背景抑制模块,用于进行背景抑制,以邻域像素点集合中背景点的灰度均值作为背景估计值,若像素点为背景点,直接抑制,若像素点为目标点,将目标点的灰度值减去背景估计值进行背景抑制。
进一步的,邻域结构模板包括: 其中,f表示前景像素点灰度平均值,b表示背景像素点灰度平均值。
进一步的,获取似然值模块包括:
获取结构向量模块,用于基于图像像素点pi,j处的邻域像素点的灰度值得到结构向量:
其中,Ii-1,j表示坐标为(i-1,j)的邻域像素点的灰度值,Ii,j-1表示坐标为(i,j-1)的邻域像素点的灰度值,Ii,j表示坐标为(i,j)的邻域像素点的灰度值,Ii+1,j表示坐标为(i+1,j)的邻域像素点的灰度值,Ii+1,j+1表示坐标为(i+1,j+1)的邻域像素点的灰度值,Ii-1,j-1表示坐标为(i-1,j-1)的邻域像素点的灰度值,Ii-1,j+1表示坐标为(i-1,j+1)的邻域像素点的灰度值,Ii+1,j-1表示坐标为(i+1,j-1)的邻域像素点的灰度值;
相似性度量模块,用于计算结构向量与邻域结构模板 的相似性度量pkl,
其中,d表示数据维度,为邻域结构模板的核宽矩阵;
计算似然值模块,用于将邻域结构模板分为目标模板和背景模板,利用相似性度量得到像素点为背景点的背景似然值Pb和像素点为目标点的目标似然值Pf,
Pl=max{pkl|k=1,2,...,8}
其中,Pl为像素点与邻域结构模板的相似性度量,
P1f表示水平垂直方向的四邻域像素点与目标模板的相似性度量,P2f表示对角线方向的四邻域像素点与目标模板的相似性度量,则有:
P1f=max{pkl|k=1,2,...,4;l=8}
P2f=max{pkl|k=5,6,...,8;l=8}
则目标似然值Pf为:Pf=(P1f+P2f)/2。
进一步的,背景估计值为
其中,card(Sbk)为Sbk中的像素点个数,Sbk表示背景点集合,pm+i,n+j表示坐标为(m+i,n+j)的邻域像素点,Ii+m,j+n表示坐标为(m+i,n+j)的邻域像素点的灰度值。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有益效果:
(1)本发明基于结构向量与邻域结构模板得到相似性度量,利用相似性度量得到目标似然值与背景似然值的比值,比值大于等于阈值,则像素点为目标点,将目标点的灰度值减去背景估计值进行背景抑制,否则,像素点为背景点,直接抑制;本发明同时具有良好的抑制效果、对不同背景类型的自适应性以及满足实时计算要求。本发明不仅背景抑制效果好,并且实验中无需设置参数。
(2)本发明对图像中每一个像素点构建邻域结构模板,相较于单一地对整幅图构建同样的邻域结构模板,更接近真实的邻域结构;像素点的邻域结构模板既包含了图像的灰度信息,又能描述邻域结构,反映了图像的局部背景特性;减少了数据量。
(3)优选的,基于邻域像素点的灰度值得到结构向量组,以水平垂直方向四邻域和倾斜方向组合的方式构建结构向量,既不丢失八邻域信息,又减少了计算量,同时采用结构向量旋转的方式能够避免结构向量与邻域结构模板恰好错位的情况。另外,根据目标点和背景点的特点以不同的方式求目标似然值和背景似然值,使得目标点的衡量尽可能准确,同时尽可能排除背景点。
(4)阈值的选取可以根据实际需求做出调整,具有灵活性。基于似然比,将图像中的像素点划分为背景点与目标点的方法适用于不同的图像背景,因而使得该背景抑制方法不局限于某一类特定的背景。
(5)以邻域像素点集合中背景点的灰度均值作为背景估计值,可以使背景估计值更准确,同时在背景抑制时最大限度地保留目标信息;另外,对背景点直接抑制的方式可以使背景抑制效果更佳。
附图说明
图1是本发明实施例提供的一种基于邻域结构的背景抑制方法的流程图;
图2是本发明实施例1提供的一阶四邻域结构的模板示意图;
图3是本发明实施例1提供的最大类间方差分割的三种情况示意图;
图4是本发明实施例1提供的对于两种情况的分割示意图;
图5是本发明实施例1提供的邻域示意图;
图6是本发明实施例提供的边缘抑制效果图;
图7是本发明实施例提供的不同背景下的背景抑制效果图;
图8是本发明实施例提供的不同背景抑制效果比较图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
如图1所示,一种基于邻域结构的背景抑制方法,包括:
(1)选取图像中的任意一个像素点,得到以像素点为中心的邻域像素点集合,将邻域像素点集合划分为前景像素点子集和背景像素点子集,基于前景像素点子集的前景像素点灰度平均值和背景像素点子集的背景像素点灰度平均值构建像素点的邻域结构模板;
(2)基于邻域像素点的灰度值得到结构向量,计算结构向量与结构模板的相似性度量,将邻域结构模板分为目标模板和背景模板,利用相似性度量得到像素点为背景点的背景似然值和像素点为目标点的目标似然值;
(3)若目标似然值与背景似然值的比值大于等于阈值,则像素点为目标点,否则,像素点为背景点;
(4)以邻域像素点集合中背景点的灰度均值作为背景估计值,若像素点为背景点,直接抑制,若像素点为目标点,将目标点的灰度值减去背景估计值进行背景抑制。
进一步的,步骤(1)还包括:
(1-1)选取图像中任意一个坐标为(i,j)的像素点pi,j,像素点pi,j的灰度值为Ii,j,以像素点pi,j为中心的八邻域像素点集合为S={pm+i,n+j|-1≤m,n≤1},pm+i,n+j表示坐标为(m+i,n+j)的邻域像素点,以像素点pi,j为中心的八邻域像素点对应的灰度值集合为I={Ii+m,j+n|-1≤m,n≤1},Ii+m,j+n表示坐标为(m+i,n+j)的邻域像素点的灰度值;
(1-2)将S划分为前景像素点子集Sf和背景像素点子集Sb:
(1-2-1)若Ii,j是集合I的最大值,前景像素点子集和背景像素点子集分别为:
Sf={pi,j},
Sb=S-Sf;
(1-2-2)若Ii,j是集合I的次大值,前景像素点子集和背景像素点子集分别为:
Sf={pi,j,pk,l|Ik,l=max{Ii+m,j+n,-1≤m,n≤1}},
Sb=S-Sf,
其中,pk,l为集合I的最大值对应位置为k,l所表示的像素点,Ik,l为集合I的最大值;
(1-2-3)T为灰度集I的分割阈值,设定初值T0=min{Ii+m,j+n,-1≤m,n≤1}和步长ΔT,不同的增长率Q下对应不同的分割阈值TQ:
TQ=T0+Q*ΔT,Q=0,1,2,...,[max{Ii+m,j+n,-1≤m,n≤1}/ΔT]
按阈值TQ将S划分为对应的前景像素点子集和背景像素点子集
按下列各式计算不同的增长率Q下的类间方差ICVQ:
其中,为前景像素频率,为背景像素频率,fQ为前景灰度均值,bQ为背景灰度均值,M为灰度集合的均值,式子中card(*)表示集合*中的像素点个数,表示中的像素点个数,card(S)表示S中的像素点个数,表示中的像素点个数:
根据计算得到的一系列的ICVQ,Q=0,1,2,…,取使类间方差最大的增长率K=argmax{ICVQ,Q=0,1,2,...},则:
Sb=S-Sf
(1-3)计算前景像素点子集的前景像素点灰度平均值f和背景像素点子集的背景像素点灰度平均值b:
其中,card(Sf)表示Sf中像素点的个数,card(Sb)表示Sb中像素点的个数,则图像像素点pi,j处的邻域结构模板由以下12种模板组成:
进一步的,步骤(2)还包括:
(2-1)基于图像像素点pi,j处的邻域像素点的灰度值得到结构向量:
其中,Ii-1,j表示坐标为(i-1,j)的邻域像素点的灰度值,Ii,j-1表示坐标为(i,j-1)的邻域像素点的灰度值,Ii,j表示坐标为(i,j)的邻域像素点的灰度值,Ii+1,j表示坐标为(i+1,j)的邻域像素点的灰度值,Ii+1,j+1表示坐标为(i+1,j+1)的邻域像素点的灰度值,Ii-1,j-1表示坐标为(i-1,j-1)的邻域像素点的灰度值,Ii-1,j+1表示坐标为(i-1,j+1)的邻域像素点的灰度值,Ii+1,j-1表示坐标为(i+1,j-1)的邻域像素点的灰度值,为像素点pi,j水平垂直方向的四邻域构成的向量组,为像素点pi,j处对角线方向的四邻域构成的向量组;
(2-2)计算结构向量与邻域结构模板的相似性度量pkl,
其中,数据维度d=5,为邻域结构模板下的核宽矩阵,以的构造为例,其结构与相对应,可以表示为其中,Fbin和Bbin分别为Sf和Sb对应的灰度集方差:
A、若card(Sf)≠1,则转步骤B,否则:
Fbin=ε<<Bbin
B、若card(Sb)≠1,则转步骤C,否则:
Bbin=ε<<Fbin
C、Fbin和Bbin分别为Sf,Sb的方差
其中,ε表示非常小的数值;
(2-3)将邻域结构模板分为目标模板和背景模板,其中,为目标模板,其余为背景模板;利用相似性度量得到像素点为背景点的背景似然值Pb和像素点为目标点的目标似然值Pf,
Pl=max{pkl|k=1,2,...,8}
其中,Pl为像素点与邻域结构模板的相似性度量,
P1f表示水平垂直方向的四邻域像素点与目标模板的相似性度量,P2f表示对角线方向的四邻域像素点与目标模板的相似性度量,则有:
P1f=max{pkl|k=1,2,...,4;l=8}
P2f=max{pkl|k=5,6,...,8;l=8}
则目标似然值Pf为:
Pf=(P1f+P2f)/2。
进一步的,步骤(3)的具体实现方式为:
若目标似然值与背景似然值的比值l大于等于阈值L,则像素点为目标点,否则,像素点为背景点;
l=Pf/Pb
若l≥L,则pi,j∈Stgt,否则,pi,j∈Sbk,Stgt,Sbk分别为目标点集和背景点集。
进一步的,步骤(4)的具体实现方式为:
以邻域像素点集合中背景点的灰度均值作为背景估计值,若像素点为背景点,背景抑制值为0,直接抑制,若像素点为目标点,将目标点的灰度值减去背景估计值进行背景抑制;
图像中像素点pi,j的背景抑制值Isupp为:
式中,为背景估计值,card(Sbk)为Sbk中的像素点个数:
实施例1
邻域结构模板的构建源于目标点与背景点的不同特征。目标点往往结构较为单一,主要表现为目标点的灰度值高于周围点的灰度值。对于图像中的背景点,有亮区,暗区,边界,边缘点,角点等。理想情况下,假设图像只有两个灰度值,被划分为高灰度值点与低灰度值点,以1表示高灰度值,以0表示低灰度值。如若只考虑中心点与其邻域点,则可以列出所有可能的结构,按照目标的特征找出目标模板,其余为背景模板。考虑八邻域结构,一共9个像素点,每个点的选择有高值和低值两种选择,一共2∧9种,去除邻域八个点旋转以后会重合的情况,一共128种结构。考虑四邻域结构,一共5个像素点,一共2∧5种,去除邻域四个点旋转以后会重合的情况,一共12种结构。相较八邻域结构,四邻域结构的种类更少,更易于分析。一阶四邻域结构的模板示意图如图2所示,图2中(1)、(2)所示结构与图像中暗区接近,图2中(11)、(12)与图像中亮区接近,图2中(8)所示结构为理想目标结构,图2中其余所示结构接近图像中的边缘,边界和角点。
为了得到邻域结构模板中的高灰度值和低灰度值,首先需要将像素集分割为两部分,再取两部分的灰度均值作为模板中的高灰度值和低灰度值。以较为简单的最大类间分割方法来分割像素集,但仅采用最大类间分割方法,会出现三种情况,示意图如图3所示。target表示目标点,图3(1)为正确分割的情况,图3(2)所示如若中心点为目标点,其八邻域中有一点灰度值极低,用最大类间方差将像素分为两个部分时,划分情况会出现图3(2)所示的情况,灰度值极低的像素点单独为一类,其余像素点为一类,这样很可能导致将该像素的结构归为背景结构,导致目标点丢失。图3(3)所示如若目标点在角点附近,角点较亮,目标点将被归为暗区,也会导致目标丢失。因而对这三种情况作不同的处理。
图4(1)表示图3(2)所示如若中心点为目标点,其八邻域中有一点灰度值极低时的分割情况,图4(2)表示图3(3)所示如若目标点在角点附近,角点较亮时的分割情况。
一阶结构向量组的构建用到了像素点的八邻域集合,因为八邻域比四邻域包含的信息更多,如若只考虑四邻域,则漏了很多信息,因而将八邻域拆分为两个四邻域结构,水平垂直方向的四邻域与倾斜方向的四邻域。邻域示意图如图5所示,图5(1)表示水平垂直方向四邻域,图5(2)表示倾斜方向四邻域,图5(3)表示八邻域理想目标。考虑到像素点的结构与邻域像素点的具体位置无关,只有邻域像素点的分布有关,如若邻域像素点按原来的分布同一方向旋转一定得角度,角度可取{90°,180°,270°},仍为同一个结构,因而将水平垂直方向四邻域结构向量依次旋转90°,180°,270°,得到结构向量同理,将倾斜方向四邻域结构向量依次旋转90°,180°,270°,得到结构向量
在已求得邻域结构模板和一阶结构向量组后,关键在于如何衡量像素点成为目标点、背景点的可能性。采用高斯核密度函数可以衡量像素点pi,j的8个一阶结构向量与12个邻域结构模板相似性pkl;紧接着通过取8个一阶结构向量与邻域结构模板中相似性最大的值衡量像素点与邻域结构模板的相似性度量Pl;最后用均值的统计的方式取所有背景模板的均值作为像素点为背景的背景似然值Pb。但是对于目标模板,如若像素点的结构与目标模板相似,要求水平垂直方向的结构向量与倾斜方向的结构向量同时与目标模板的相似。因而,像素水平垂直方向的结构与目标模板得相似性P1f取水平垂直方向的一阶结构向量组与目标模板中相似性最大的值;像素倾斜方向的结构与目标模板得相似性P2f取倾斜方向的一阶结构向量组与目标模板中相似性最大的值,两者的均值即为目标似然值Pf。以Pf/Pb衡量一个像素点更倾向于目标点还是背景点。设定阈值L,若目标似然值与背景似然值的比值大于等于阈值L,则像素点为目标点,否则,为背景点。
对整幅图求得每个像素点的点类型后,采用基于点划分的背景抑制策略。对于被判定为背景点的像素点,背景抑制值直接设为0;对于被判定为目标点的像素点,以邻域像素点集合中背景点的灰度均值作为背景估计值,以像素点灰度值减去背景估计值的大小作为背景抑制值。对整幅图所有像素点采取同样的背景抑制策略,最终即可得到整幅图的背景抑制结果。
以表1中的数据示例图像中某个像素点及其八邻域。
表1图像中某个像素点示例
64 | 69 | 71 |
64 | 81 | 69 |
63 | 65 | 66 |
其邻域集合S={64,69,71,64,81,69,63,65,66},示例数据中中心点的灰度值高于周围像素点的灰度值,因而采用计算公式:
Sf={pi,j}
Sb=S-Sf
实施例1中具体计算为:
Sf={81}
Sb={64,69,71,64,69,63,65,66}
计算前景像素点子集的前景像素点灰度平均值f和背景像素点子集的背景像素点灰度平均值b:
f=81
基于图像像素点pi,j处的邻域像素点的灰度值得到结构向量:
中心点的灰度值为八邻域中的最大值,Sf集合中只有一个点,Fbin可以设定为一个较小的值,使Fbin远小于Bbin。计算结果为:
Fbin=10-6
因而,12个核宽矩阵可以表示为以下形式:
紧接着计算8个结构向量与12个邻域结构模板的相似性度量pkl,求得的结果如表2所示。
表2相似性计算结果
将邻域结构模板分为目标模板和背景模板,其中,为目标模板,其余为背景模板;利用相似性度量得到像素点为背景点的背景似然值Pb和像素点为目标点的目标似然值Pf,Pl=max{pkl|k=1,2,...,8}
其中,Pl为像素点与邻域结构模板的相似性度量,
Pl=5.04826×10-11
P2=0
P3=0
P4=0
P5=0
P6=0
P7=0
P9=0
P10=0
P11=0
P12=0
P1f表示水平垂直方向的四邻域像素点与目标模板的相似性度量,P2f表示对角线方向的四邻域像素点与目标模板的相似性度量,则有:
P1f=max{pkl|k=1,2,...,4;l=8}
P2f=max{pkl|k=5,6,...,8;l=8}
则目标似然值Pf为:
P1f=max{p18,p28,p38,p48}=0.0397375
P2f=max{p58,p68,p78,p88}=0.0144970
Pf=(P1f+P2f)/2=0.0271
若目标似然值与背景似然值的比值l大于等于阈值L,则像素点为目标点,否则,像素点为背景点;
l=Pf/Pb=5908758301.66552
为了最大限度的保留疑似目标点,选择L=1,本实施例1中似然比l远远大于1,因而可以判定为目标点,即:
pi,j∈Stgt
本实施例1中,中心像素点周围的点似然比均小于1,判定为背景点,最后本实施例中中心像素点背景抑制的值计算得:
实验(1)
为了考察本发明的背景抑制效果,选用边缘较为明晰的图像,如图6(a)所示,用本发明背景抑制方法对图像处理,背景抑制结果如图6(b)所示。图6中(c),(d),(e),(f)依次为去最大中值滤波,非局部均值滤波,TDLMS滤波,Top-hat滤波的结果。
为了衡量背景抑制后的效果,以信杂比(SCR)为指标,其定义为目标It减去背景μb的强度与局部背景的标准差σb的比值。
信杂比计算结果如下表3所示。
表3信杂比计算结果
从图6(b)中可以看出本发明背景抑制方法能够抑制背景边缘,保留图像中疑似目标点结构的像素点。相较于其他方法,本发明的背景抑制效果更好,信杂比大幅提高。
实验(2)不同背景下的实验
第二个实验测试的数据取四种不同的图像背景,用提出的算法在本文下运算。实验(1)中考察了本发明背景抑制方法在边缘较为明晰情况下的背景抑制效果,然而在大片云层,边缘较为明显的情况下,图像各个像素间的相关性较好。因此选取背景时,除了选取大片云层的背景图7-B1外,还选取了另外三种背景,图7-A1小片区域较多,图7-C1云层起伏较大,图图7-D1中角点较多。图7中第2列为背景抑制的结果,图7中第三列为背景抑制后以目标点为中心大小为21×21的局部结构。求得的信杂比如下表所示。表4为该实验下信杂比计算结果,结果表明对于不同的背景图像,本发明背景抑制方法均能自适应滤波,无需进行参数调整。
表4信杂比计算结果
图7-A1 | 图7-B1 | 图7-C1 | 图7-D1 | |
SCR | 6.1754 | 6.1650 | 5.2964 | 7.4425 |
实验(3)对比实验
在实验(2)的基础上,与其他方法作对比实验。对比的方法有去最大中值滤波,非局部均值滤波,TDLMS滤波,top-hat滤波。其中去最大中值滤波采用3×3的窗口。非局部均值滤波的搜索框半径R取10,相似框半径r为1,核宽h为5。TDLMS的滤波器宽为3×3,迭代步长需要根据不同的背景调整至最优。处理后的结果如下表5所示。对比结果表明提出的方法在不同的背景下效果优于其他方法。
表5信杂比计算结果
SCR | 本发明方法 | 去最大中值 | 非局部均值 | TDLMS | Top-hat |
图7-A1 | 6.1754 | 2.4008 | 3.1541 | 2.2591 | 4.7968 |
图7-B1 | 6.1650 | 1.8120 | 3.2009 | 3.9338 | 4.3172 |
图7-C1 | 5.2964 | 1.6999 | 3.9156 | 2.5714 | 2.6245 |
图7-D1 | 7.4425 | 4.1645 | 4.6758 | 4.0682 | 5.7187 |
图8展示了用该算法处理后信杂比较低的图与其他方法作对比的结果展示。图8A,图8B,图8C,图8D,图8E依次为本发明背景抑制方法,去最大中值滤波,非局部均值滤波,TDLMS滤波,Top-hat滤波的背景抑制结果,以及背景抑制后以目标点位中心的局部结构图。
实验(1)的结果表明本发明背景抑制方法可以有效滤除背景的边缘,在图像背景相关性较强时,背景抑制后信杂比远远高于其他算法的信杂比。实验(2)选取了四种不同的背景,目标强度较弱,本发明背景抑制方法仍能自适应地抑制背景,突出小目标。实验(3)在不同背景下与其他算法做对比,从计算的信杂比中可以看出本发明背景抑制方法在背景抑制的效果上依然优于其他算法。本发明背景抑制方法不仅背景抑制效果好,并且实验中无需设置参数。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种基于邻域结构的背景抑制方法,其特征在于,包括:
(1)选取图像中的任意一个像素点,得到以像素点为中心的邻域像素点集合,将邻域像素点集合划分为前景像素点子集和背景像素点子集,基于前景像素点子集的前景像素点灰度平均值和背景像素点子集的背景像素点灰度平均值构建像素点的邻域结构模板;
(2)基于邻域像素点的灰度值得到结构向量,计算结构向量与邻域结构模板的相似性度量,将邻域结构模板分为目标模板和背景模板,利用相似性度量得到像素点为背景点的背景似然值和像素点为目标点的目标似然值;
(3)若目标似然值与背景似然值的比值大于等于阈值,则像素点为目标点,否则,像素点为背景点;
(4)以邻域像素点集合中背景点的灰度均值作为背景估计值,若像素点为背景点,直接抑制,若像素点为目标点,将目标点的灰度值减去背景估计值进行背景抑制;
所述步骤(2)包括:
(2-1)基于图像像素点pi,j处的邻域像素点的灰度值得到结构向量:
其中,Ii-1,j表示坐标为(i-1,j)的邻域像素点的灰度值,Ii,j-1表示坐标为(i,j-1)的邻域像素点的灰度值,Ii,j表示坐标为(i,j)的邻域像素点的灰度值,Ii+1,j表示坐标为(i+1,j)的邻域像素点的灰度值,Ii+1,j+1表示坐标为(i+1,j+1)的邻域像素点的灰度值,Ii-1,j-1表示坐标为(i-1,j-1)的邻域像素点的灰度值,Ii-1,j+1表示坐标为(i-1,j+1)的邻域像素点的灰度值,Ii+1,j-1表示坐标为(i+1,j-1)的邻域像素点的灰度值;
(2-2)计算结构向量与邻域结构模板的相似性度量pkl,1≤k≤8,1≤l≤12,
其中,d表示数据维度,为邻域结构模板的核宽矩阵;
(2-3)将邻域结构模板分为目标模板和背景模板,利用相似性度量得到像素点为背景点的背景似然值Pb和像素点为目标点的目标似然值Pf,
Pl=max{pkl|k=1,2,...,8}
其中,Pl为像素点与邻域结构模板的相似性度量,
P1f表示水平垂直方向的四邻域像素点与目标模板的相似性度量,P2f表示对角线方向的四邻域像素点与目标模板的相似性度量,则有:
P1f=max{pkl|k=1,2,...,4;l=8}
P2f=max{pkl|k=5,6,...,8;l=8}
则目标似然值Pf为:Pf=(P1f+P2f)/2,为目标模板。
2.如权利要求1所述的一种基于邻域结构的背景抑制方法,其特征在于,所述邻域结构模板包括: 其中,f表示前景像素点灰度平均值,b表示背景像素点灰度平均值。
3.如权利要求1所述的一种基于邻域结构的背景抑制方法,其特征在于,所述背景估计值为
其中,card(Sbk)为Sbk中的像素点个数,Sbk表示背景点集合,pm+i,n+j表示坐标为(m+i,n+j)的邻域像素点,Ii+m,j+n表示坐标为(m+i,n+j)的邻域像素点的灰度值。
4.一种基于邻域结构的背景抑制系统,其特征在于,包括:
构建邻域结构模块,用于选取图像中的任意一个像素点,得到以像素点为中心的邻域像素点集合,将邻域像素点集合划分为前景像素点子集和背景像素点子集,基于前景像素点子集的前景像素点灰度平均值和背景像素点子集的背景像素点灰度平均值构建像素点的邻域结构模板;
获取似然值模块,用于基于邻域像素点的灰度值得到结构向量,计算结构向量与邻域结构模板的相似性度量,将邻域结构模板分为目标模板和背景模板,利用相似性度量得到像素点为背景点的背景似然值和像素点为目标点的目标似然值;
判断像素点类型模块,用于判断像素点类型,若目标似然值与背景似然值的比值大于等于阈值,则像素点为目标点,否则,像素点为背景点;
进行背景抑制模块,用于进行背景抑制,以邻域像素点集合中背景点的灰度均值作为背景估计值,若像素点为背景点,直接抑制,若像素点为目标点,将目标点的灰度值减去背景估计值进行背景抑制;
所述获取似然值模块包括:
基于图像像素点pi,j处的邻域像素点的灰度值得到结构向量:
其中,Ii-1,j表示坐标为(i-1,j)的邻域像素点的灰度值,Ii,j-1表示坐标为(i,j-1)的邻域像素点的灰度值,Ii,j表示坐标为(i,j)的邻域像素点的灰度值,Ii+1,j表示坐标为(i+1,j)的邻域像素点的灰度值,Ii+1,j+1表示坐标为(i+1,j+1)的邻域像素点的灰度值,Ii-1,j-1表示坐标为(i-1,j-1)的邻域像素点的灰度值,Ii-1,j+1表示坐标为(i-1,j+1)的邻域像素点的灰度值,Ii+1,j-1表示坐标为(i+1,j-1)的邻域像素点的灰度值;
计算结构向量与邻域结构模板的相似性度量pkl,1≤k≤8,1≤l≤12,
其中,d表示数据维度,为邻域结构模板的核宽矩阵;
将邻域结构模板分为目标模板和背景模板,利用相似性度量得到像素点为背景点的背景似然值Pb和像素点为目标点的目标似然值Pf,
Pl=max{pkl|k=1,2,...,8}
其中,Pl为像素点与邻域结构模板的相似性度量,
P1f表示水平垂直方向的四邻域像素点与目标模板的相似性度量,P2f表示对角线方向的四邻域像素点与目标模板的相似性度量,则有:
P1f=max{pkl|k=1,2,...,4;l=8}
P2f=max{pkl|k=5,6,...,8;l=8}
则目标似然值Pf为:Pf=(P1f+P2f)/2,为目标模板。
5.如权利要求4所述的一种基于邻域结构的背景抑制系统,其特征在于,所述邻域结构模板包括: 其中,f表示前景像素点灰度平均值,b表示背景像素点灰度平均值。
6.如权利要求5所述的一种基于邻域结构的背景抑制系统,其特征在于,所述背景估计值为
其中,card(Sbk)为Sbk中的像素点个数,Sbk表示背景点集合,pm+i,n+j表示坐标为(m+i,n+j)的邻域像素点,Ii+m,j+n表示坐标为(m+i,n+j)的邻域像素点的灰度值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710415948.7A CN107194896B (zh) | 2017-06-05 | 2017-06-05 | 一种基于邻域结构的背景抑制方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710415948.7A CN107194896B (zh) | 2017-06-05 | 2017-06-05 | 一种基于邻域结构的背景抑制方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107194896A CN107194896A (zh) | 2017-09-22 |
CN107194896B true CN107194896B (zh) | 2019-12-17 |
Family
ID=59876397
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710415948.7A Expired - Fee Related CN107194896B (zh) | 2017-06-05 | 2017-06-05 | 一种基于邻域结构的背景抑制方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107194896B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108491796B (zh) * | 2018-03-22 | 2021-10-22 | 电子科技大学 | 一种时域周期点目标检测方法 |
CN109002777B (zh) * | 2018-06-29 | 2021-03-30 | 电子科技大学 | 一种面向复杂场景的红外小目标检测方法 |
CN109949324B (zh) * | 2019-02-01 | 2022-04-22 | 广西科技大学 | 基于非经典感受野非线性亚单元响应的轮廓检测方法 |
CN110930426B (zh) * | 2019-11-11 | 2022-09-20 | 中国科学院光电技术研究所 | 一种基于峰域形态辨识的弱小点目标提取方法 |
CN112132911B (zh) * | 2020-11-25 | 2021-03-16 | 恒银金融科技股份有限公司 | 一种用于红外数字图像检测的阈值自适应方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101957984A (zh) * | 2010-08-30 | 2011-01-26 | 西安电子科技大学 | 基于非局部萎缩因子参数估计的图像去噪方法 |
CN104537690A (zh) * | 2014-12-26 | 2015-04-22 | 中国科学院西安光学精密机械研究所 | 一种基于最大值-时间索引联合的运动点目标检测方法 |
CN105096317A (zh) * | 2015-07-03 | 2015-11-25 | 吴晓军 | 一种复杂背景中的高性能相机全自动标定方法 |
CN105242258A (zh) * | 2010-11-16 | 2016-01-13 | 高通股份有限公司 | 用于基于超声反射信号的对象位置估计的系统和方法 |
-
2017
- 2017-06-05 CN CN201710415948.7A patent/CN107194896B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101957984A (zh) * | 2010-08-30 | 2011-01-26 | 西安电子科技大学 | 基于非局部萎缩因子参数估计的图像去噪方法 |
CN105242258A (zh) * | 2010-11-16 | 2016-01-13 | 高通股份有限公司 | 用于基于超声反射信号的对象位置估计的系统和方法 |
CN104537690A (zh) * | 2014-12-26 | 2015-04-22 | 中国科学院西安光学精密机械研究所 | 一种基于最大值-时间索引联合的运动点目标检测方法 |
CN105096317A (zh) * | 2015-07-03 | 2015-11-25 | 吴晓军 | 一种复杂背景中的高性能相机全自动标定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107194896A (zh) | 2017-09-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107194896B (zh) | 一种基于邻域结构的背景抑制方法和系统 | |
US9147255B1 (en) | Rapid object detection by combining structural information from image segmentation with bio-inspired attentional mechanisms | |
CN109886218B (zh) | 基于超像素统计相异性的sar图像舰船目标检测方法 | |
CN106485740B (zh) | 一种结合稳定点和特征点的多时相sar图像配准方法 | |
CN107169972B (zh) | 一种非合作目标快速轮廓跟踪方法 | |
CN109086724B (zh) | 一种加速的人脸检测方法及存储介质 | |
CN108198201A (zh) | 一种多目标跟踪方法、终端设备及存储介质 | |
Zhu et al. | Edge-preserving guided filtering based cost aggregation for stereo matching | |
CN110390338B (zh) | 一种基于非线性引导滤波与比率梯度的sar高精度匹配方法 | |
CN112861983A (zh) | 图像匹配方法、装置、电子设备以及存储介质 | |
Zhang et al. | Fast moving pedestrian detection based on motion segmentation and new motion features | |
CN106600613B (zh) | 基于嵌入式gpu的改进lbp红外目标检测方法 | |
Wang et al. | Inshore ship detection using high-resolution synthetic aperture radar images based on maximally stable extremal region | |
Peng et al. | Infrared small-target detection based on multi-directional multi-scale high-boost response | |
CN107451595A (zh) | 基于混合算法的红外图像显著性区域检测方法 | |
Mo et al. | Sea-sky line detection in the infrared image based on the vertical grayscale distribution feature | |
Wang et al. | Road detection based on illuminant invariance and quadratic estimation | |
Ma et al. | Hierarchical Sea‐Land Segmentation for Panchromatic Remote Sensing Imagery | |
CN113205494B (zh) | 基于自适应尺度图像块加权差测量的红外小目标检测方法和系统 | |
Li et al. | Salient object detection based on meanshift filtering and fusion of colour information | |
Talai et al. | A fast edge detection using fuzzy rules | |
CN108830864A (zh) | 图像分割方法 | |
CN110363105B (zh) | 一种全极化sar图像相干斑抑制方法 | |
CN116030252A (zh) | 融合边缘信息的模糊迭代sar图像超像素分割方法 | |
Zhang et al. | A road extraction method based on high resolution remote sensing image |
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: 20191217 Termination date: 20200605 |