CN104951765A - 基于形状先验信息和视觉对比度的遥感图像目标分割方法 - Google Patents
基于形状先验信息和视觉对比度的遥感图像目标分割方法 Download PDFInfo
- Publication number
- CN104951765A CN104951765A CN201510341210.1A CN201510341210A CN104951765A CN 104951765 A CN104951765 A CN 104951765A CN 201510341210 A CN201510341210 A CN 201510341210A CN 104951765 A CN104951765 A CN 104951765A
- Authority
- CN
- China
- Prior art keywords
- image
- contrast
- confidence map
- region
- pending image
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 93
- 230000008569 process Effects 0.000 claims abstract description 23
- 238000012549 training Methods 0.000 claims abstract description 14
- 230000011218 segmentation Effects 0.000 claims description 30
- 230000004927 fusion Effects 0.000 claims description 22
- 230000000007 visual effect Effects 0.000 claims description 13
- 238000003709 image segmentation Methods 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000013341 scale-up Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 6
- 239000000284 extract Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000010191 image analysis Methods 0.000 description 2
- 238000000513 principal component analysis Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 206010038743 Restlessness Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002902 bimodal effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000004304 visual acuity Effects 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/13—Satellite images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2415—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on parametric or probabilistic models, e.g. based on likelihood ratio or false acceptance rate versus a false rejection rate
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/26—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
- G06V10/267—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion by performing operations on regions, e.g. growing, shrinking or watersheds
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Probability & Statistics with Applications (AREA)
- Astronomy & Astrophysics (AREA)
- Remote Sensing (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于形状先验信息和视觉对比度的遥感图像目标分割方法,所述方法包括:对样本图像进行训练得到分类器;利用所述分类器对待处理图像进行处理,得到所述待处理图像的图像目标的形状先验置信图;对所述待处理图像进行分割,得到所述待处理图像的图像目标的对比度置信图;通过所述形状先验置信图和对比度置信图得到所述待处理图像中的目标。本发明通过所述形状先验置信图和对比度置信图对待处理图像进行分析,利用了分类器的分类正确性和待处理图像自身的特性,提高了待处理图像的图像目标的分割精度。
Description
技术领域
本发明属于图像处理技术领域,尤其涉及一种基于形状先验信息和视觉对比度的遥感图像目标分割方法。
背景技术
随着卫星遥感技术的发展,遥感卫星分辨率逐渐提高,世界各国已经发射了越来越多的高分辨率遥感卫星,例如:QuickBird、WorldView-2、GeoEye-1等。2014年8月发射的WorldView-3,能提供0.3米分辨率图像,飞机、舰船、车辆等目标的细节清晰可见。随着我国高分辨率对地观测***国家科技重大专项的实施,我国也将拥有更多高分辨率遥感卫星。因此,利用卫星遥感图像进行目标分析与理解越来越引起关注。其中,飞机、舰船、车辆等感兴趣目标的分割,是在目标检测的基础上,确定目标边界并提取出目标的技术,是判定目标形状与识别目标类型的重要步骤,直接影响到图像内容分析的准确性。但是,遥感图像中目标背景复杂、噪声较大,阴影、地面斑点等因素经常影响目标的准确分割,当前方法经常难以准确分割目标。
研究人员已经提出了大量的感兴趣目标分割方法,其中常用的典型目标分割方法从分割策略上讲,可以分为基于阈值、基于对比度、引入形状先验信息等方法。(1)基于阈值分割的方法准确性与图像的光谱分布密切相关,当直方图呈现双峰或近似双峰特性时,能够选取阈值,但当背景复杂或者目标光谱与背景较为接近时,很难准确分割目标。因此,此类方法适用于背景相对较简单的情况,鲁棒性和适应性较差。(2)基于对比度的方法假定目标与背景存在较大的光谱差异,通过度量图像中各个区域和图像整体光谱的差异来确定目标区域。但是该方法容易受图像中的阴影噪声干扰,因为当光照较强,目标亮度较高时,阴影则相对较显著,利用该方法进行分割容易得到阴影区域,目标区域缺失严重。(3)引入形状先验信息的方法主要是利用主成分分析、傅立叶描述子等方法学习目标形状信息,但是当前的形状先验信息建模方法对于颜色先验信息利用不充分,难以分割复杂背景中的目标。综合而言,对于遥感图像中复杂背景下的感兴趣目标分割,当前方法容易受阴影、油污等背景噪声干扰,难以取得准确的分割结果。
发明内容
有鉴于此,本发明实施例期望提供一种基于形状先验信息和视觉对比度的遥感图像目标分割方法,至少能解决现有图像目标识别精度低等技术问题。
本发明实施例的技术方案是这样实现的:
本发明实施例提供了基于形状先验信息和视觉对比度的遥感图像目标分割方法,所述方法包括:
对样本图像进行训练得到分类器;
利用所述分类器对待处理图像进行处理,得到所述待处理图像的图像目标的形状先验置信图;
对所述待处理图像进行分割,得到所述待处理图像的图像目标的对比度置信图;
通过所述形状先验置信图和对比度置信图得到所述待处理图像的图像目标。
上述方案中,所述对样本图像进行训练得到分类器包括:
从所述样本图像中提取第一图像特征;
对所述第一图像特征进行归一化处理,并通过线性支持向量机训练归一化处理后的所述第一图像特征,得到包含所述样本图像的图像目标的形状先验信息的分类器。
上述方案中,所述第一图像特征包括像素的RGB值、坐标值和直方图对比度。
上述方案中,所述利用所述分类器对待处理图像进行处理,得到所述待处理图像的图像目标的形状先验置信图包括:
从所述待处理图像中提取第二图像特征;
通过所述分类器得到所述第二图像特征属于所述待处理图像的图像目标的概率;
对所述概率进行归一化处理得到包含所述待处理图像的图像目标的形状先验置信图。
上述方案中,所述对所述待处理图像进行分割,得到所述待处理图像的图像目标的对比度置信图包括:
将所述待处理图像按设定尺度分割成图像区域;
计算所述图像区域之间的局部对比度和全局对比度;
根据所述局部对比度和全局对比度得到每个所述图像区域的对比度;
将所述对比度进行求和,得到所述待处理图像的图像目标的对比度置信图。
上述方案中,所述全局对比度的计算过程为:
其中:ri m为所述待处理图像按尺度m进行分割得到的第i个区域;gc(ri m)为区域ri m上的全局对比度;为区域ri m上的平均颜色值;区域上的平均颜色值;j为除第i个区域外的其他区域;Nm为所述待处理图像按尺度m进行分割得到的区域的总数。
上述方案中,所述对比度的计算过程为:
f(ri m)=lc(ri m)×gc(ri m)×g(x,y)
其中,f(ri m)为区域的对比度;lc(ri m)为区域的局部对比度;g(x,y)为高斯衰减函数;x和y分别是区域的像素点的横坐标和纵坐标。
上述方案中,所述高斯衰减函数的计算过程为:
其中,x0和y0分别是区域中心的横坐标和纵坐标;σx为图像三分之一的宽度;σy为图像三分之一的高度。
上述方案中,所述通过所述形状先验置信图和对比度置信图得到所述待处理图像的图像目标包括:
通过所述形状先验置信图和对比度置信图得到目标融合置信图;
将所述目标融合置信图进行二值化处理得到连通区域和背景区域;
根据所述连通区域的第一阈值和背景区域的第二阈值对所述待处理图像的像素点进行判断,得到所述待处理图像的图像目标。
上述方案中,所述目标融合置信图的计算过程为:
C=CC·exp(k·CS)
其中,C为目标融合置信图的像素值;CC为形状先验置信图的像素值;CS为对比度置信图的像素值;k为比例系数。
本发明实施例所提供的基于形状先验信息和视觉对比度的遥感图像目标分割方法,通过所述形状先验置信图和对比度置信图对待处理图像进行分析,利用了分类器的分类正确性和待处理图像自身的特性,提高了待处理图像的图像目标的识别精度。
附图说明
图1为实施例1的基于形状先验信息和视觉对比度的遥感图像目标分割方法的流程图;
图2为实施例2的基于形状先验信息和视觉对比度的遥感图像目标分割方法的示意图;
图3为实施例2的分类器的获取过程示意图;
图4为实施例2的形状先验置信图的获取过程示意图;
图5为实施例2的方法与现有方法的效果比较图。
为了能明确实现本发明的实施例的结构,在图中标注了特定的尺寸、结构和器件,但这仅为示意需要,并非意图将本发明限定在该特定尺寸、结构、器件和环境中,根据具体需要,本领域的普通技术人员可以将这些器件和环境进行调整或者修改,所进行的调整或者修改仍然包括在后附的权利要求的范围中。
具体实施方式
在以下的描述中,将描述本发明的多个不同的方面,然而,对于本领域内的普通技术人员而言,可以仅仅利用本发明的一些或者全部结构或者流程来实施本发明。为了解释的明确性而言,阐述了特定的数目、配置和顺序,但是很明显,在没有这些特定细节的情况下也可以实施本发明。在其他情况下,为了不混淆本发明,对于一些众所周知的特征将不再进行详细阐述。
实施例1
为了解决现有图像目标识别精度低等技术问题,本发明实施例提供了基于形状先验信息和视觉对比度的遥感图像目标分割方法,如图1所示,本实施例的方法包括以下步骤:
步骤S101:对样本图像进行训练得到分类器;
本实施例的样本图像是能够精确识别图像目标的遥感图像,通过对样本图像训练,就可以得到对图像目标的分类器。
步骤S102:利用所述分类器对待处理图像进行处理,得到所述待处理图像的图像目标的形状先验置信图;
通过所述分类器对待处理图像进行处理,能够得到待处理图像的图像目标的基本形状,即形状先验置信图。此处的待处理图像也是遥感图像,后续不再赘述。
步骤S103:对所述待处理图像进行分割,得到所述待处理图像的图像目标的对比度置信图;
为了进一步将待处理图像的图像目标分离出来,可以按不同尺度对待处理图像进行分割,得到对比度置信图。
步骤S104:通过所述形状先验置信图和对比度置信图得到所述待处理图像的图像目标。
通过所述形状先验置信图和对比度置信图对待处理图像进行分析,利用了分类器的分类正确性和待处理图像自身的特性,提高了待处理图像的图像目标的识别精度。
具体的,所述步骤S101具体包括:
步骤S1011:从所述样本图像中提取第一图像特征;所述第一图像特征包括像素的RGB(红RED、绿GREEN、蓝BLUE)值、坐标值和直方图对比度。
步骤S1012:对所述第一图像特征进行归一化处理,并通过线性支持向量机训练归一化处理后的所述第一图像特征,得到包含所述样本图像的图像目标的形状先验信息的分类器。
所述步骤S102具体包括:
步骤S1021:从所述待处理图像中提取第二图像特征;第二图像特征和第一图像特征一样,也是像素的RGB值、坐标值和直方图对比度。
步骤S1022:通过所述分类器得到所述第二图像特征属于所述待处理图像的图像目标的概率;
步骤S1023:对所述概率进行归一化处理得到包含所述待处理图像的图像目标的形状先验置信图。
所述步骤S103具体包括:
步骤S1031:将所述待处理图像按设定尺度分割成图像区域;
此处的按设定尺度是指按设定的大小将待处理图像进行分割,尺度大小的具体取值需要根据待处理图像的实际情况确定。
步骤S1032:计算所述图像区域之间的局部对比度和全局对比度;
局部对比度是指某区域与局部邻域间的颜色差异;全局对比度是指某区域与图像全局其它区域间的颜色差异。
步骤S1033:根据所述局部对比度和全局对比度得到每个所述图像区域的对比度;
步骤S1034:将所述对比度进行求和,得到所述待处理图像的图像目标的对比度置信图。
具体的,所述全局对比度的计算过程为:
其中:ri m为所述待处理图像按尺度m进行分割得到的第i个区域;gc(ri m)为区域ri m上的全局对比度;为区域ri m上的平均颜色值;区域上的平均颜色值;j为除第i个区域外的其他区域;Nm为所述待处理图像按尺度m进行分割得到的区域的总数。
所述对比度的计算过程为:
f(ri m)=lc(ri m)×gc(ri m)×g(x,y)
其中,f(ri m)f为区域的对比度;lc(ri m)为区域的局部对比度;g(x,y)为高斯衰减函数;x和y分别是区域的像素点的横坐标和纵坐标。
所述高斯衰减函数的计算过程为:
其中,x0和y0分别是区域中心的横坐标和纵坐标;σx为待处理图像三分之一的宽度;σy为待处理图像三分之一的高度。
所述步骤S104具体包括:
步骤S1041:通过所述形状先验置信图和对比度置信图得到目标融合置信图;
步骤S1042:将所述目标融合置信图进行二值化处理得到连通区域和背景区域;
步骤S1043:根据所述连通区域的第一阈值和背景区域的第二阈值对所述待处理图像的像素点进行判断,得到所述待处理图像的图像目标。
所述目标融合置信图的计算过程为:
C=CC·exp(k·CS)
其中,C为目标融合置信图的像素值;CS为形状先验置信图的像素值;CC为对比度置信图的像素值;k为比例系数,用于调节形状先验置信图与对比度置信图权重。其中,所有的C组合起来构成目标融合置信图;类似的,CS组合起来构成形状先验置信图;CC组合起来构成对比度置信图。为了方便起见,下述直接用C代表目标融合置信图;CS代表形状先验置信图;CC代表对比度置信图。
实施例2
本实施例对本发明进行详细说明。
本实施例的示意图如图2所示,包括以下步骤:
步骤1:对已标注的图像进行训练,逐像素提取r、g、b、x、y、hc特征。其中,r、g和b特征分别为像素的RGB值;x和y特征为像素在图像中的坐标位置;hc特征为基于直方图的对比度;对特征进行归一化,通过训练得到包含目标形状先验信息的分类器。
将样本图像逐像素的进行二值标注:目标区域的像素标注为1,背景区域的像素标注为-1。训练包含目标先验状信息的分类器时,使用线性支持向量机进行训练。图3显示了使用WorldView-2图像作为样本,训练包含目标形状先验信息分类器的过程。
步骤2:对待处理图像逐像素提取r、g、b、x、y、hc特征,利用步骤1中获得的分类器逐像素预测图像中各像素点属于目标的概率,进而得到待处理图像中目标的形状先验置信图CS。
对输入图像中的每个像素利用分类器分类后,得到范围为-1至+1的值,表示该像素属于目标区域的置信度。全图计算结束后,归一化到[0,1],即得到形状先验置信图CS。图4显示了对输入图像利用分类器逐像素预测后得到的形状先验置信图。
步骤3:对待处理图像进行多尺度的超像素分割,以ri m表示所述分割后尺度m上的第i个区域,分别计算各个区域的局部对比度、全局对比度,并结合高斯衰减函数得到各个区域的对比度计算结果f(ri m);将多尺度下各个区域对比度的计算结果进行融合,得到待处理图像中目标的对比度置信图CC。
将待处理图像使用Felzenszwalb和Huttenlocher提出的超像素分割方法分割,生成设定个尺度的分割结果。超像素分割方法可以将具有相似纹理、颜色特征的相邻像素分成图像块;在尺度m上,图像被分割为区域序列{ri m},i=1,…,Nm,Nm表示该尺度m上分割得到区域的数量。区域ri m的局部对比度计算方法为2011年Huaizu Jiang在论文“Automatic salient object segmentation basedon context and shape prior”中所提出的方法,局部对比度计算方法认为,图像中的目标区域与其相邻区域的颜色差异较大,通过计算区域与其领域间的颜色差异得到图像中的目标区域。
区域ri m的全局对比度的计算方法为:
其中,|ri m|表示区域ri m中包含的像素数;和分别表示区域ri m和的平均颜色值。
高斯衰减函数g(x,y)定义为:
其中,x0和y0分别是区域中心的横坐标和纵坐标;σx为待处理图像三分之一的宽度,σy为待处理图像三分之一的高度。这样越靠近图像中心的地方权重越大。
对单个尺度中各个区域按照下面公式进行计算,得到单尺度下的各区域对比度计算结果f(ri m):
f(ri m)=lc(ri m)×gc(ri m)×g(x,y)
其中,lc(ri m)为区域的局部对比度;gc(ri m)为区域的全局对比度。
将多个尺度上的对比度计算结果f(ri m)相加求得,就可以得到输入图像中的对比度置信图CC。
步骤4:将形状先验置信图CS和对比度置信图CC进行融合,得到集成形状先验信息和对比度信息的目标融合置信图C;将融合置信图C分割后的结果作为三分图利用GrabCut算法进行M次迭代分割,得到目标分割结果。其中,三分图就是将图像中的像素归为三类:目标区域、背景区域、可能为目标也可能为背景的未知区域;利用GrabCut算法进行迭代分割时,若图像中的目标长度大于128像素时,最优迭代次数M为3,若目标长度小于128像素时,最优迭代次数M为2。
目标融合置信图C的计算方法为:
C=CC·exp(k·CS)
通过大量实验观察,发现形状先验置信图CS更具有区分性,因此,使用指数函数来强调CS的作用,k用来调节置信图融合时CC和CS的比例,优选为4。
利用目标融合置信图C生成三分图时,将C归一化到[0,1]区间,首先对C利用5×5的模板进行中值滤波,利用粗分割阈值T(优选为0.3)进行二值化,大于T的像素点标注为1,小于T的像素点标注为0;对二值化后的图像提取包含像素数目最大的连通区域,区域外的像素点指定为背景区域;区域内的像素点依据目标融合置信图C进行判断,若目标融合置信图C中像素值小于背景阈值T1(优选为0.15),则认定该像素点为背景区域;若目标融合置信图C中像素值大于目标阈值T2(优选为0.85),则认定该像素点为目标区域;若目标融合置信图C中像素值大于背景阈值T1而小于目标阈值T2,则认定该像素点为可能的目标区域。其中,背景阈值T1和目标阈值T2都是归一化到[0,1]区间后的取值。
本发明在微机平台上进行了算法实验,实验对象为遥感图像中的飞机目标,实验过程中使用约0.5米分辨率、160×160像素的遥感图像(全色与多光谱融合图像)。样本量共500张图片,其中200张图片用于训练,300张图片用于测试。图4给出了本发明方法各步骤的处理结果,其中第一列为输入图像,第二列为形状先验置信图CS,第三列为对比度置信图CC,第四列为利用CS和CC计算得到的目标融合置信图C,第五列为目标分割结果。从图中可以看出,在先验知识指导下得到的目标置信图能够很好的反映飞机的整体形状,但是周围噪声点较多;通过对比度分析得到的目标置信图能够较好地提取与背景差异较大的区域;将两者融合,能够有效地排除噪声干扰,准确分割目标。
本发明方法与现有比较先进的目标分割方法的对比结果参见图5。对比试验中的Normalized Cuts(下述简写为NC)方法是文章“Normalized cuts and imagesegmentation”中提出的;Local Binary Fitting(下述简写为LBF)方法是文章“Implicit active contours driven by local binary fitting energy”中提出的;Contextual Hypergraph Modeling(下述简写为CHM)方法是文章“Contextualhypergraph modeling for salient object detection”中提出的;Region-based Contrast(下述简写为RC)方法是文章“Global contrast based salient region detection”中提出的;GrabCut(下述简写为GC)方法是文章“Grabcut-Interactive foregroundextraction using iterated graph cuts”中提出的。从图5中可以看出,当背景干扰因素较多时,本发明方法仍能够较准确的提取感兴趣目标。如图5前四行所示,其它方法容易受到飞机本身阴影的影响,难以准确分割目标;本发明方法在先验形状信息的指导下,可以较好地排除阴影因素的干扰;同时,如图5后两行所示,其它方法容易受到地面斑块、油污等因素的影响,本发明方法通过在对比度分析中使用高斯衰减函数,并结合先验形状信息的指导,可以较好地排除地面斑块、油污等因素的干扰。
本发明使用F测度(F-measure)方法在客观指标上与上述方法进行对比,计算方法如下:
其中,Precision为查准率,Recall为查全率。较大的F测度值说明算法分割结果较为准确。表1显示了本发明与上述方法的对比如果。从表1中可以看到,在本发明的方法F测度值是最高的,这说明本发明方法能够较准确地分割遥感图像中的飞机目标。
表1 本遥感图像目标分割方法实验结果
方法 | F测度 |
NC | 0.623±0.002 |
LBF | 0.675±0.008 |
CHM | 0.795±0.004 |
RC | 0.809±0.008 |
GC | 0.839±0.006 |
本发明方法 | 0.898±0.002 |
本发明所提供的基于形状先验信息和视觉对比度的方法主要为提高遥感图像分析与理解专门提出的。但显然,本方法也适用于医学图像分析以及普通成像设备如数码相机所拍摄图像的分析,所取得的有益效果也是相似的。
以上对本发明所提供的基于形状先验信息和视觉对比度的遥感图像目标分割方法作了详细的说明,但显然本发明的具体实现形式并不局限于此。对于本技术领域的一般技术人员来说,在不背离本发明的权利要求范围的情况下对它进行的各种显而易见的改变都在本发明的保护范围之内。
本发明所提供的基于形状先验信息和视觉对比度的遥感图像目标分割方法具有如下优点:
1.利用像素空间位置与颜色、全局对比度等底层视觉特征建立包含目标先验形状信息的分类器,可以准确地预测输入图像中目标的大致形状,保证在复杂环境中的鲁棒性;
2.多尺度的对比度信息分析可以保证区域提取的完整性,在此基础上与在先验知识指导下得到的目标区域进行融合,可以有效削弱阴影、噪声等因素的干扰,保证分割结果的准确性;
3.可以对具有规则形状的感兴趣目标进行分割,适应性较强。
在本申请所提供的几个实施例中,应该理解到,所揭露的设备和方法,可以通过其它的方式实现。以上所描述的设备实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,如:多个单元或组件可以结合,或可以集成到另一个***,或一些特征可以忽略,或不执行。另外,所显示或讨论的各组成部分相互之间的耦合、或直接耦合、或通信连接可以是通过一些接口,设备或单元的间接耦合或通信连接,可以是电性的、机械的或其它形式的。
上述作为分离部件说明的单元可以是、或也可以不是物理上分开的,作为单元显示的部件可以是、或也可以不是物理单元,即可以位于一个地方,也可以分布到多个网络单元上;可以根据实际的需要选择其中的部分或全部单元来实现本实施例方案的目的。
另外,在本发明各实施例中的各功能单元可以全部集成在一个处理模块中,也可以是各单元分别单独作为一个单元,也可以两个或两个以上单元集成在一个单元中;上述集成的单元既可以采用硬件的形式实现,也可以采用硬件加软件功能单元的形式实现。
本领域普通技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成,前述的程序可以存储于一计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包括:移动存储设备、只读存储器(Read-Only Memory,ROM)、随机存取存储器(Random Access Memory,RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (10)
1.一种基于形状先验信息和视觉对比度的遥感图像目标分割方法,其特征在于,所述方法包括:
对样本图像进行训练得到分类器;
利用所述分类器对待处理图像进行处理,得到所述待处理图像的图像目标的形状先验置信图;
对所述待处理图像进行分割,得到所述待处理图像的图像目标的对比度置信图;
通过所述形状先验置信图和对比度置信图得到所述待处理图像的图像目标。
2.根据权利要求1所述的方法,其特征在于,所述对样本图像进行训练得到分类器包括:
从所述样本图像中提取第一图像特征;
对所述第一图像特征进行归一化处理,并通过线性支持向量机训练归一化处理后的所述第一图像特征,得到包含所述样本图像的图像目标的形状先验信息的分类器。
3.根据权利要求2所述的方法,其特征在于,所述第一图像特征包括像素的RGB值、坐标值和直方图对比度。
4.根据权利要求1所述的方法,其特征在于,所述利用所述分类器对待处理图像进行处理,得到所述待处理图像的图像目标的形状先验置信图包括:
从所述待处理图像中提取第二图像特征;
通过所述分类器得到所述第二图像特征属于所述待处理图像的图像目标的概率;
对所述概率进行归一化处理得到包含所述待处理图像的图像目标的形状先验置信图。
5.根据权利要求1所述的方法,其特征在于,所述对所述待处理图像进行分割,得到所述待处理图像的图像目标的对比度置信图包括:
将所述待处理图像按设定尺度分割成图像区域;
计算所述图像区域之间的局部对比度和全局对比度;
根据所述局部对比度和全局对比度得到每个所述图像区域的对比度;
将所述对比度进行求和,得到所述待处理图像的图像目标的对比度置信图。
6.根据权利要求5所述的方法,其特征在于,所述全局对比度的计算过程为:
其中:为所述待处理图像按尺度m进行分割得到的第i个区域;为区域上的全局对比度;为区域上的平均颜色值;区域上的平均颜色值;j为除第i个区域外的其他区域;Nm为所述待处理图像按尺度m进行分割得到的区域的总数。
7.根据权利要求6所述的方法,其特征在于,所述对比度的计算过程为:
其中,为区域的对比度;为区域的局部对比度;g(x,y)为高斯衰减函数;x和y分别是区域的像素点的横坐标和纵坐标。
8.根据权利要求7所述的方法,其特征在于,所述高斯衰减函数的计算过程为:
其中,x0和y0分别是区域中心的横坐标和纵坐标;σx为待处理图像三分之一的宽度;σy为待处理图像三分之一的高度。
9.根据权利要求1所述的方法,其特征在于,所述通过所述形状先验置信图和对比度置信图得到所述待处理图像的图像目标包括:
通过所述形状先验置信图和对比度置信图得到目标融合置信图;
将所述目标融合置信图进行二值化处理得到连通区域和背景区域;
根据所述连通区域的第一阈值和背景区域的第二阈值对所述待处理图像的像素点进行判断,得到所述待处理图像的图像目标。
10.根据权利要求9所述的方法,其特征在于,所述目标融合置信图的计算过程为:
C=CC·exp(k·CS)
其中,C为目标融合置信图的像素值;CC为形状先验置信图的像素值;CS为对比度置信图的像素值;k为比例系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510341210.1A CN104951765B (zh) | 2015-06-18 | 2015-06-18 | 基于形状先验信息和视觉对比度的遥感图像目标分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510341210.1A CN104951765B (zh) | 2015-06-18 | 2015-06-18 | 基于形状先验信息和视觉对比度的遥感图像目标分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104951765A true CN104951765A (zh) | 2015-09-30 |
CN104951765B CN104951765B (zh) | 2017-12-15 |
Family
ID=54166407
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510341210.1A Expired - Fee Related CN104951765B (zh) | 2015-06-18 | 2015-06-18 | 基于形状先验信息和视觉对比度的遥感图像目标分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104951765B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106991411A (zh) * | 2017-04-17 | 2017-07-28 | 中国科学院电子学研究所 | 基于深度形状先验的遥感图像目标精细化提取方法 |
CN107527054A (zh) * | 2017-09-19 | 2017-12-29 | 西安电子科技大学 | 基于多视角融合的前景自动提取方法 |
CN109215049A (zh) * | 2018-07-18 | 2019-01-15 | 中国科学院自动化研究所 | 基于多尺度三维先验信息的屋顶分割方法、***与设备 |
CN109886133A (zh) * | 2019-01-25 | 2019-06-14 | 北京市遥感信息研究所 | 一种基于光学遥感图像的船只检测方法及*** |
CN109902555A (zh) * | 2019-01-11 | 2019-06-18 | 西安电子科技大学 | 一种基于对象的多尺度层次表达学习的变化检测方法 |
CN110807485A (zh) * | 2019-10-31 | 2020-02-18 | 电子科技大学 | 基于高分遥感影像二分类语义分割图融合成多分类语义图的方法 |
CN112528962A (zh) * | 2021-01-01 | 2021-03-19 | 生态环境部卫星环境应用中心 | 一种基于高分卫星遥感影像的牧区牛马群的监测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102722891A (zh) * | 2012-06-12 | 2012-10-10 | 大连理工大学 | 一种图像显著度检测的方法 |
CN104680529A (zh) * | 2015-02-26 | 2015-06-03 | 浙江工业大学 | 基于形状先验和相似性约束的颈动脉内外壁自动分割方法 |
-
2015
- 2015-06-18 CN CN201510341210.1A patent/CN104951765B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102722891A (zh) * | 2012-06-12 | 2012-10-10 | 大连理工大学 | 一种图像显著度检测的方法 |
CN104680529A (zh) * | 2015-02-26 | 2015-06-03 | 浙江工业大学 | 基于形状先验和相似性约束的颈动脉内外壁自动分割方法 |
Non-Patent Citations (1)
Title |
---|
窦燕 等: "结合基元对比度与边界先验的显著性区域检测", 《高技术通讯》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106991411B (zh) * | 2017-04-17 | 2019-03-22 | 中国科学院电子学研究所 | 基于深度形状先验的遥感图像目标精细化提取方法 |
CN106991411A (zh) * | 2017-04-17 | 2017-07-28 | 中国科学院电子学研究所 | 基于深度形状先验的遥感图像目标精细化提取方法 |
CN107527054B (zh) * | 2017-09-19 | 2019-12-24 | 西安电子科技大学 | 基于多视角融合的前景自动提取方法 |
CN107527054A (zh) * | 2017-09-19 | 2017-12-29 | 西安电子科技大学 | 基于多视角融合的前景自动提取方法 |
CN109215049A (zh) * | 2018-07-18 | 2019-01-15 | 中国科学院自动化研究所 | 基于多尺度三维先验信息的屋顶分割方法、***与设备 |
CN109215049B (zh) * | 2018-07-18 | 2020-11-06 | 中国科学院自动化研究所 | 基于多尺度三维先验信息的屋顶分割方法、***与设备 |
CN109902555A (zh) * | 2019-01-11 | 2019-06-18 | 西安电子科技大学 | 一种基于对象的多尺度层次表达学习的变化检测方法 |
CN109902555B (zh) * | 2019-01-11 | 2020-09-22 | 西安电子科技大学 | 一种基于对象的多尺度层次表达学习的变化检测方法 |
CN109886133A (zh) * | 2019-01-25 | 2019-06-14 | 北京市遥感信息研究所 | 一种基于光学遥感图像的船只检测方法及*** |
CN109886133B (zh) * | 2019-01-25 | 2020-11-03 | 北京市遥感信息研究所 | 一种基于光学遥感图像的船只检测方法及*** |
CN110807485A (zh) * | 2019-10-31 | 2020-02-18 | 电子科技大学 | 基于高分遥感影像二分类语义分割图融合成多分类语义图的方法 |
CN110807485B (zh) * | 2019-10-31 | 2022-05-03 | 电子科技大学 | 基于高分遥感影像二分类语义分割图融合成多分类语义图的方法 |
CN112528962A (zh) * | 2021-01-01 | 2021-03-19 | 生态环境部卫星环境应用中心 | 一种基于高分卫星遥感影像的牧区牛马群的监测方法 |
CN112528962B (zh) * | 2021-01-01 | 2021-07-20 | 生态环境部卫星环境应用中心 | 一种基于高分卫星遥感影像的牧区牛马群的监测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104951765B (zh) | 2017-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gavankar et al. | Automatic building footprint extraction from high-resolution satellite image using mathematical morphology | |
CN104951765A (zh) | 基于形状先验信息和视觉对比度的遥感图像目标分割方法 | |
CN108510467B (zh) | 基于深度可变形卷积神经网络的sar图像目标识别方法 | |
CN107767400B (zh) | 基于层次化显著性分析的遥感图像序列动目标检测方法 | |
Li et al. | Road lane detection with gabor filters | |
CN105678760A (zh) | 一种基于Canny边缘检测算法的绝缘子图像识别方法 | |
CN104408482A (zh) | 一种高分辨率sar图像目标检测方法 | |
CN104657980A (zh) | 一种改进的基于Meanshift的多通道图像分割算法 | |
Yuan et al. | Road segmentation in aerial images by exploiting road vector data | |
KR101742115B1 (ko) | 건물 인식을 위한 멀티뷰 영상에서의 인라이어 선택 및 잉여 제거 방법 | |
CN104680545A (zh) | 光学图像中存在显著目标的检测方法 | |
CN103927759A (zh) | 一种航空图像自动云检测方法 | |
Zang et al. | Traffic lane detection using fully convolutional neural network | |
KR20190059083A (ko) | 영상 분할 기반의 해양 상황 인식 장치 및 그 방법 | |
CN104050674B (zh) | 一种显著性区域检测方法及装置 | |
CN105023264A (zh) | 一种结合对象性和背景性的红外图像显著特征检测方法 | |
Liu et al. | Vehicle detection from aerial color imagery and airborne LiDAR data | |
CN109635679B (zh) | 一种实时的靶纸定位及环线识别方法 | |
Raikar et al. | Automatic building detection from satellite images using internal gray variance and digital surface model | |
Deepan et al. | Road recognition from remote sensing imagery using machine learning | |
CN104616007B (zh) | 一种基于显著性检测和颜色直方图模型的车辆识别方法 | |
CN106228553A (zh) | 高分辨率遥感图像阴影检测装置与方法 | |
Abraham et al. | Unsupervised building extraction from high resolution satellite images irrespective of rooftop structures | |
Moughal et al. | An automatic unsupervised method based on context-sensitive spectral angle mapper for change detection of remote sensing images | |
Poornima et al. | A method to align images using image segmentation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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 |
Granted publication date: 20171215 |
|
CF01 | Termination of patent right due to non-payment of annual fee |