CN102289812B - 一种基于先验形状和cv模型的目标分割方法 - Google Patents
一种基于先验形状和cv模型的目标分割方法 Download PDFInfo
- Publication number
- CN102289812B CN102289812B CN 201110247798 CN201110247798A CN102289812B CN 102289812 B CN102289812 B CN 102289812B CN 201110247798 CN201110247798 CN 201110247798 CN 201110247798 A CN201110247798 A CN 201110247798A CN 102289812 B CN102289812 B CN 102289812B
- Authority
- CN
- China
- Prior art keywords
- centerdot
- theta
- psi
- partiald
- 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
Images
Landscapes
- Image Analysis (AREA)
Abstract
一种基于先验形状以及CV模型的目标分割方法,首先选取先验形状构建符号距离函数,再根据仿射变换参数向量对先验形状符号距离函数进行仿射变换,仿射变换参数在对水平集函数进行迭代时会发生改变并最终趋于稳定,最后根据水平集迭代公式演化活动轮廓曲线,同时根据各仿射参数迭代公式求出下一次的各仿射变换参数。本发明在保持先验形状模型具有旋转、缩放和平移不变性的基础上,增加了X、Y方向拉伸以及剪切不变约束能量项,通过对先验形状的自适应变换的拓展,本发明可以实现对复杂背景下姿态变换较大的空间目标进行较好的分割。
Description
技术领域
本发明涉及图像处理技术,具体是一种基于先验形状和CV模型的目标分割方法。
背景技术
目前,在已有的各种目标分割方法中,曲线演化方法对目标的分割有相当好的结果,具体包括Snake方法、活动轮廓线方法、变形模型以及水平集方法等,参数化的Snake方法允许与模型直接交互,并且模型的表达紧凑,有利于模型的快速实现,但是难以处理模型拓扑结构的变化。基于变分水平集方法的活动轮廓模型能自然地处理演化曲线或曲面拓扑结构的变化,并能够自然地将边界信息和区域信息整合在一起。
Mumford在1989年提出通过分片光滑函数的最佳逼近解决边缘检测问题的M-S水平集模型(Mumford D,Shah J.Optimal approximation by piecewise smooth functionsand associated variational problems.Communication on Pure and Applied Mathematics,1989,42(5):577-685.),Chan和Vese在M-S模型的基础上提出了简化的CV水平集模型(Chan T F,Vese L A.Active contours without edges.IEEE Transactions on ImageProcessing,2001,10(2):266-277.),该模型的一个显著优点就是全局优化,在边界模糊处也能获得较好的分割结果,并且初始轮廓线可以放置在图像的任意位置。但该模型是基于灰度相似性来划分目标区域的,因此该模型存在三个方面的缺陷:①不能分割目标灰度与背景相似的目标,②不能有效分割纹理图像,③不能分割被遮挡、数据缺失的目标。
为了克服以上缺陷,近年来国内外许多学者重点研究了把先验形状信息和水平集结合进行图像分割的问题。Tony Chan等提出了利用先验形状知识进行图像分割处理的可形变模型(Chan T and Zhu W.Level set based shape prior segmentation.IEEEConference on Computer Vision and Pattern Recognition(CVPR),2005:1164-1170.),该模型在CV模型的基础上加入先验形状信息,能够在图像中分割出数据缺失、被遮挡或目标灰度与背景相似的目标。但该模型中的先验形状项只具有旋转、缩放和平移不变特性,对于发生剪切或在X、Y方向有不同拉伸系数的目标,上述模型的分割效果较差。
发明内容
本发明为了克服现有的先验形状模型对目标的自适应变化存在的局限性,即先验形状项只具有X、Y方向等大小缩放、平移以及旋转三个不变特性,提供一种基于先验形状和CV模型的目标分割方法,本发明在原有先验形状变分水平集模型基础上增加了X、Y方向拉伸以及剪切不变约束能量项,通过对先验形状的自适应变换的拓展,本方法可以实现对复杂背景下姿态变换较大的目标进行准确分割。
本发明通过以下技术方案实现:
一种基于先验形状和CV模型的目标分割方法,为基于变分水平集方法的图像分割模型,它在融入图像区域信息基础上,自然地融入先验形状信息,很好地处理了一些传统的几何活动轮廓模型无法解决的问题。
本发明所构建的模型对应的能量函数表示如下:
E(c1,c2,φ,ψ)=ECV(c1,c2,φ)+λEshape(φ,ψ) (1)
其中:c1和c2分别是φ对应的零水平集曲线内、外部区域的图像灰度平均值,λ为先验形状能量项对应的权重系数,CV能量项ECV(c1,c2,φ)和先验形状项Eshape(φ,ψ)分别表示如下:
Eshape(φ,ψ)=∫Ω(H(φ)-H(ψ))2dxdy (3)
其中:
XT=(Sx,Sy,θ,shx,shy,Tx,Ty),其中Sx及Sy分别为X、Y方向上的两个缩放系数,θ为角度旋转参数,shx,shy分别为X、Y方向上的两个修剪参数,Tx,Ty分别为X、Y方向上的两个平移参数;f为定义在图像域Ω上的灰度值;ψ为先验形状的水平集函数ψ0经式(4)的仿射变换得到,用于对式(3)的先验形状能量项进行约束。
对式(1)的能量项E(c1,c2,φ,ψ)进行最小化,采用变分法及梯度下降流方程构建演化模型对应的数值计算格式,令对式(1)能量项E(c1,c2,φ,ψ)求变分,可得c1,c2,仿射变换参数XT=(Sx,Sy,θ,shx,shy,Tx,Ty)以及演化的水平集函数φ对应的最速梯度下降方程为:
本发明包括如下步骤:
1、利用先验形状样本构建符号距离函数:
从序列图中选取标准姿态的目标图像作为构建先验形状模型样本,然后对选中的先验形状样本进行阈值分割,阈值的大小根据样本背景的复杂度进行手动调整,其次采用数学形态算子对阈值分割后的图像进行背景抑制并用Sobel算子提取目标的二值边缘轮廓,最后利用生成的二值轮廓线构建先验形状符号距离函数;
2、根据初始化的仿射变换参数,利用公式(4)、(5)对先验形状模型样本对应的符号距离函数进行仿射变换:
本发明拓展了原有先验形状的仿射变换项,即仿射变换从四元关系(a,b,r,θ)拓展为(Sx,Sy,θ,shx,shy,Tx,Ty),其中Sx及Sy分别为X、Y方向上的两个缩放系数,θ为角度旋转参数,shx,shy分别为X、Y方向上的两个修剪参数,Tx,Ty分别为X、Y方向上的两个平移参数;
3、利用式(6~14)对图像进行水平集全局迭代运算:
数值化计算采用正则化函数Dirac函数分别代替H(z)和δ(z),使得下降流方程(14)能够作用在所有的水平集上,从而能够对图像区域的所有目标进行自动检测,并且使能量函数E(c1,c2,φ,ψ)达到全局极小值。
本发明的有益效果是:
1、由于本发明拓展了先验形状对目标的自适应变化方式,在使用改进的先验形状水平集模型对目标进行分割时,可以分割出复杂背景下姿态变化较大且复杂的目标。
2、本发明测试的实验仿真数据是由美国AGI公司卫星仿真工具包STK软件产生的一组视频场景,利用STK8.0计算轨道姿态信息,并得到从观测平台看到的目标飞行器模拟图像信息。为使仿真数据尽量逼真,以检验目标检测、跟踪与识别方法,在生成模拟图像过程中,考虑了以下几点:(1)目标飞行器在视场中的大小必须按照两飞行器真实距离变化;(2)目标飞行器要有一定的旋转;(3)目标飞行器不应始终处于视场中央;(4)背景应有一定的变化,有星空背景,有地球背景,或星空地球背景交替变化。
下面结合附图和实施案例对本发明作详细说明。
附图说明
图1是从地球背景和星空背景下选取的5幅典型卫星姿态图片。
图2是选取的先验形状样本(图2(a)及其对应的二值图图2(b)、二值轮廓图图2(c)。
图3是基于图2(b)生成的人工合成图。
图4是图3(a)对应的初始轮廓线图以及不同模型下的分割效果图。
图5是图3(b)对应的初始轮廓线图以及不同模型下的分割效果图。
图6是图3(c)对应的初始轮廓线图以及不同模型下的分割效果图。
图7是图1中的各图对应的初始轮廓线叠加图。
图8是传统的CV模型对图1中各图的分割结果。
图9是传统的基于先验形状的CV模型对图1中各图的分割结果。
图10是本发明基于先验形状的CV模型对图1中各图的分割结果。
图11是本发明方法的流程图。
具体实施方式
以下结合附图和实施例对本发明的作详细说明:本实施例在以本发明技术方案为前提下进行实施的例子,给出了详细的实施方式和过程,但本发明的保护范围不应限于下述的实施例。
从太空对目标进行跟踪时,采集的序列图像由于拍摄的角度以及自身的姿态调整,提取的图像目标轮廓线之间基本符合仿射变换关系,本实施例选取地球背景和星空背景下的5幅典型卫星姿态图片用于检验本发明所阐述的“基于先验形状和CV模型的目标分割方法”性能,如图1所示,本实施例同时选取图2(a)作为先验形状样本,本发明基于先验形状和CV模型的目标分割方法的具体实施步骤如下:
1、利用先验形状样本构建符号距离函数:
对先验形状样本图2(a)进行阈值分割,本实施例分割的阈值选为128,然后利用形态学算子对抑制后的图像进行形态学处理并生成二值图,即2(b),其次利用Sobel算子对图2(b)进行边缘提取,得到先验形状对应的二值轮廓线,即图2(c),最后基于图2(b)及图2(c)构建先验形状样本(图2(a)对应的符号距离函数ψ与ψ0;2、本实施例初始化的仿射参数依次为:(1,1,0,0,0,0,0),即式(3)中先验形状符号距离函数ψ与ψ0相同,在对能量函数式(3)的最小化过程中,仿射参数向量(Sx,Sy,θ,shx,shy,Tx,Ty)对应的数值不断变化,从而使得先验形状能够逐渐匹配待分割的目标,在活动轮廓线稳定收敛于待分割的目标边界时,仿射参数对应的数值趋于稳定,本实施案例中对应的稳定值为(1,1,0,0,0,0,0),即先验形状不再发生变化。
3、给定水平集的迭代次数为n=1000。
4、利用式(6~14)对图3和图1进行水平集迭代运算,当活动轮廓线稳定时,获得待分割的目标边界。
图4(b)、图5(b)、图6(b)为传统CV模型的分割结果,由图可得,当被分割的目标发生遮挡时,传统的CV模型会将目标连同遮挡物一起检测出来,图4(c)、图5(c)、图6(c)为传统的基于先验形状的CV模型分割结果,由图可得,当被分割的目标相对先验形状发生剪切(图3(c))或者在X、Y缩放系数不同时(图3(b)),传统的基于先验形状的CV模型无法给出正确的分割结果,原因就在于先验形状约束项对放射变换的局限性。图4(d)、图5(d)、图6(d)为本发明改进的基于先验形状的CV模型分割结果,由图可得,由于引入了X、Y方向的拉伸及剪切不变特性,改进的图像分割模型可以很好的分割出被遮挡的各个目标。
由图8(e)可得,在目标或背景区域具有较好的同质性,且受到光照、阴影、遮挡、杂波等因数影响较小时,传统的CV方法可以获得比较满意的分割结果。否则传统的CV方法分割的效果较差,如图8(a)、8(b)、8(c)、8(d)所示。由图9可得,当被分割的目标相对先验形状的剪切系数很小并且在X、Y方向的拉伸系数相差也很小时,传统的基于先验形状的CV模型可以分割出比较满意的结果,如图9(a)、9(b)、9(c);反之传统的基于先验形状的CV模型的分割效果则比较差,如图9(d)、9(e)。相对于传统CV模型和基于先验形状的CV模型,本发明改进的变分水平集合模型可以显著地提高复杂地球以及星空背景下的卫星分割效果,如图10所示。
Claims (1)
1.一种基于先验形状和CV模型的目标分割方法,其特征在于包括下列步骤:
①利用先验形状样本构建符号距离函数:
从序列图中选取标准姿态的目标图像作为先验形状样本,然后对选中的先验形状样本进行阈值分割,阈值的大小根据样本背景的复杂度进行手动调整,其次采用数学形态算子对阈值分割后的图像进行背景抑制并用Sobel算子提取目标的边缘二值轮廓线,最后利用所述的二值轮廓线构建先验形状符号距离函数ψ0;
②设定并根据初始化的仿射变换参数,对先验形状样本的符号距离函数ψ0进行仿射变换,得到变换后的符号距离函数ψ:
仿射变换参数为Sx,Sy,θ,shx,shy,Tx,Ty,其中Sx及Sy分别为X、Y方向上的两个缩放系数,θ为角度旋转参数,shx,shy分别为X、Y方向上的修剪参数,Tx,Ty分别为X、Y方向上的平移参数;
其中:
XT=(Sx,Sy,θ,shx,shy,Tx,Ty);
③利用下列公式(3)~公式(11)对图像进行水平集全局迭代运算,当活动轮廓线稳定时,即获得待分割的目标边界:
其中f为定义在图像域Ω上的灰度值,c1和c2分别是φ对应的零水平集曲线内、外部区域的图像灰度平均值,λ为先验形状能量项对应的权重系数,ψx,ψy分别为水平集函数ψ在x和y方向上的梯度,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110247798 CN102289812B (zh) | 2011-08-26 | 2011-08-26 | 一种基于先验形状和cv模型的目标分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110247798 CN102289812B (zh) | 2011-08-26 | 2011-08-26 | 一种基于先验形状和cv模型的目标分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102289812A CN102289812A (zh) | 2011-12-21 |
CN102289812B true CN102289812B (zh) | 2013-03-27 |
Family
ID=45336202
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110247798 Expired - Fee Related CN102289812B (zh) | 2011-08-26 | 2011-08-26 | 一种基于先验形状和cv模型的目标分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102289812B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102663427B (zh) * | 2012-03-29 | 2014-03-05 | 浙江大学 | 基于先验形状稀疏凸组合的同步目标分割和识别方法 |
CN103886622B (zh) * | 2012-12-21 | 2017-10-31 | 腾讯科技(深圳)有限公司 | 自动图像区域划分的实现方法及实现装置 |
CN103093477B (zh) * | 2013-02-08 | 2016-01-20 | 河北大学 | 一种基于食道超声的主动脉瓣快速分割方法 |
CN103606145B (zh) * | 2013-10-29 | 2016-08-17 | 河北大学 | 基于帧间形状约束的gcv模型分割主动脉瓣超声图像序列方法 |
CN104657736B (zh) * | 2013-11-19 | 2017-10-27 | 中国科学院沈阳自动化研究所 | 一种基于主动轮廓的声呐图像水雷目标识别方法 |
CN103646398A (zh) * | 2013-12-04 | 2014-03-19 | 山西大学 | 一种皮肤镜病灶自动分割方法 |
CN103971377A (zh) * | 2014-05-27 | 2014-08-06 | 中国科学院遥感与数字地球研究所 | 基于先验形状水平集分割的建筑物提取方法 |
CN104732529B (zh) * | 2015-03-05 | 2017-07-28 | 北京空间机电研究所 | 一种遥感图像形状特征配准方法 |
CN106296672A (zh) * | 2016-08-03 | 2017-01-04 | 南京邮电大学 | 一种结合梯度信息的分块式局部加权的图像分割方法 |
CN106570867B (zh) * | 2016-10-18 | 2019-03-29 | 浙江大学 | 基于灰度形态学能量法的活动轮廓模型图像快速分割方法 |
CN106651889B (zh) * | 2016-11-30 | 2019-08-16 | 太原科技大学 | 一种x射线焊接接头缺陷分割方法及其分割*** |
CN107424153B (zh) * | 2017-04-18 | 2020-08-14 | 辽宁科技大学 | 基于深度学***集的人脸分割方法 |
CN107240114B (zh) * | 2017-05-27 | 2019-11-12 | 南京林业大学 | 一种基于点距离函数形状约束的半自动医学图像分割方法 |
CN107909587A (zh) * | 2017-12-26 | 2018-04-13 | 淮阴工学院 | 基于alr‑cv模型与边缘转换的图像分割算法 |
CN113344878B (zh) * | 2021-06-09 | 2022-03-18 | 北京容联易通信息技术有限公司 | 一种图像处理方法及*** |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7177471B2 (en) * | 2002-02-01 | 2007-02-13 | Siemens Corporate Research, Inc. | Shape priors for level set representations |
US7391882B2 (en) * | 2004-01-29 | 2008-06-24 | Siemens Medical Solutions Usa, Inc. | Prior knowledge, level set representations and visual grouping |
CN1870006A (zh) * | 2005-04-19 | 2006-11-29 | 西门子共同研究公司 | 用于水平集分割的形状和亮度先验的有效核密度估计 |
CN101964112B (zh) * | 2010-10-29 | 2013-04-03 | 上海交通大学 | 自适应的基于先验形状的图像分割方法 |
-
2011
- 2011-08-26 CN CN 201110247798 patent/CN102289812B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102289812A (zh) | 2011-12-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102289812B (zh) | 一种基于先验形状和cv模型的目标分割方法 | |
Das et al. | Use of salient features for the design of a multistage framework to extract roads from high-resolution multispectral satellite images | |
CN104933755A (zh) | 一种静态物体重建方法和*** | |
CN102663820B (zh) | 三维头部模型重建方法 | |
CN102722890B (zh) | 基于光流场模型的非刚性心脏图像分级配准方法 | |
CN109544677A (zh) | 基于深度图像关键帧的室内场景主结构重建方法及*** | |
CN102722697B (zh) | 一种无人飞行器视觉自主导引着陆的目标跟踪方法 | |
CN103822616A (zh) | 一种图分割与地形起伏约束相结合的遥感影像匹配方法 | |
CN101571951B (zh) | 基于邻域概率密度函数特征的水平集图像分割方法 | |
CN101980285B (zh) | 结合gmrf先验的非局部图像修复方法 | |
CN101866497A (zh) | 基于双目立体视觉的智能三维人脸重建方法及*** | |
CN101777178A (zh) | 一种图像修复方法 | |
CN110619368B (zh) | 行星表面导航特征成像匹配检测方法 | |
Guillaume et al. | Curvature tensor based triangle mesh segmentation with boundary rectification | |
CN101964112B (zh) | 自适应的基于先验形状的图像分割方法 | |
CN104751493A (zh) | 基于梯度纹理特征的稀疏跟踪方法 | |
CN104282000A (zh) | 一种基于旋转及尺度变换的图像修复方法 | |
CN104036299B (zh) | 一种基于局部纹理aam的人眼轮廓跟踪方法 | |
CN115511952A (zh) | 一种基于深度相机的装配式铺面三维点云预处理方法 | |
CN115471749A (zh) | 地外探测无监督学习的多视角多尺度目标识别方法及*** | |
CN107784284A (zh) | 人脸识别方法及*** | |
CN103914817B (zh) | 一种基于区域划分和插值的多光谱和全色图像融合方法 | |
CN102034275A (zh) | 面向大尺度变形的非刚性注册方法 | |
CN103745220A (zh) | 一种图像仿射局部不变特征的获取方法及装置 | |
Liu et al. | A fragment fracture surface segmentation method based on learning of local geometric features on margins used for automatic utensil reassembly |
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: 20130327 Termination date: 20150826 |
|
EXPY | Termination of patent right or utility model |