CN103886599B - 一种基于血管内超声图像的血管roi分割方法 - Google Patents
一种基于血管内超声图像的血管roi分割方法 Download PDFInfo
- Publication number
- CN103886599B CN103886599B CN201410115141.8A CN201410115141A CN103886599B CN 103886599 B CN103886599 B CN 103886599B CN 201410115141 A CN201410115141 A CN 201410115141A CN 103886599 B CN103886599 B CN 103886599B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- image
- msup
- 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 41
- 210000004204 blood vessel Anatomy 0.000 title claims abstract description 34
- 238000002608 intravascular ultrasound Methods 0.000 claims abstract description 52
- 230000011218 segmentation Effects 0.000 claims abstract description 43
- 239000012528 membrane Substances 0.000 claims abstract description 36
- 239000011159 matrix material Substances 0.000 claims description 29
- 241000270295 Serpentes Species 0.000 claims description 13
- 238000013507 mapping Methods 0.000 claims description 11
- 238000002604 ultrasonography Methods 0.000 claims description 8
- 210000004351 coronary vessel Anatomy 0.000 claims description 3
- 125000004122 cyclic group Chemical group 0.000 claims description 3
- 238000003708 edge detection Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000002792 vascular Effects 0.000 abstract description 5
- 238000000605 extraction Methods 0.000 abstract description 2
- 238000012800 visualization Methods 0.000 abstract description 2
- 210000004937 luminal membrane Anatomy 0.000 abstract 2
- 238000005192 partition Methods 0.000 abstract 1
- 230000006870 function Effects 0.000 description 26
- 210000002808 connective tissue Anatomy 0.000 description 8
- 238000003745 diagnosis Methods 0.000 description 5
- 238000003709 image segmentation Methods 0.000 description 4
- 208000024172 Cardiovascular disease Diseases 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000003143 atherosclerotic effect Effects 0.000 description 2
- 208000026106 cerebrovascular disease Diseases 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000003902 lesion Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 206010008118 cerebral infarction Diseases 0.000 description 1
- 208000029078 coronary artery disease Diseases 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 208000010125 myocardial infarction Diseases 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Landscapes
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Image Analysis (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明涉及一种基于血管内超声图像的血管ROI分割方法。所述方法首先分割血管的管腔区域以及管腔膜轮廓,通过定位管腔区域的中心以获得参数主动轮廓模型的初始轮廓,继而通过收敛得到血管的中外膜轮廓曲线,中外膜的提取充分利用了管腔区域信息的先验知识。最后,将中外膜轮廓曲线以内区域作为ROI,通过全局最小化活动轮廓模型实现对血管斑块的分割。本发明实现了血管ROI的管腔膜中外膜以及斑块的轮廓信息可视化,与基于统计学的IVUS图像分割方法相比,摒弃了其复杂的统计建模过程且分割结果不受IVUS图像伪影和斑块特征的影响;省掉了对IVUS图像中外膜边缘进行初始轮廓的预分割步骤,提高了分割效率。
Description
技术领域
本发明涉及计算机医学图像分析领域,特别涉及一种基于血管内超声IVUS(Intravascular Ultrasound)图像的血管感兴趣区域ROI(region of interest)分割方法。
背景技术
目前在大多数国家,心脑血管疾病已逐步成为人类死亡的最主要因素之一。冠状动脉粥样硬化病变是引起心肌梗塞和脑梗塞的主要病因,如果在早期就能够做到对硬化病变部分形态的识别以及诊断,这将对冠状动脉疾病的诊疗有着重大的意义。IVUS便是这样一种心血管疾病的超声诊断方法,血管内超声图像可以为医生显示实时的血管壁形态斑块形态进而对心血管疾病的临床诊断提供了指导意义。
然而,对于采集来的IVUS图像,尽管其成像可以显示血管壁的结构以及斑块的形态,但在实际诊疗中,IVUS图像中医生所感兴趣区域(ROI)为血管管腔区域,中外膜区域,以及管腔与中外膜之间的斑块区域。在IVUS图像的感兴趣区域中,医生往往需要目测或者根据经验确定管腔膜中外膜以及斑块轮廓的边缘,这样便不可避免地导致诊断结果不能非常客观的反映实际情况,同时也会给医生带来一定的工作困难。
因此,通过计算机对血管内超声图像进行快速准确自动地分割就显得很有必要,目前,血管内超声图像的计算机分割算法主要有三种:第一种为基于统计学的图像分割方法,文献1(G.Mendizabal-Ruiz,M.Rivera,et al.,“A probabilistic segmentationmethod for the identification of luminal borders in intrvascular ultrasoundimages”,IEEE Conference on Computer Vision and pattern Recognition,pp.1-8,2008.)通过对IVUS图像的灰度分布进行统计学建模实现血管内超声图像的分割,但是由于血管内超声图像中的伪影以及斑块等较为复杂的图像特征的存在会大大降低统计建模的准确性。第二种方法主要是通过机器学习的手段实现血管内超声图像分割,但是这种方法模型过于复杂,在实际的运用中往往会受到较多的限制,第三种是基于活动轮廓模型的算法,文献2(宫延新孙丰荣等,“基于血液斑点噪声抑制和T-Snake模型的血管内超声图像边缘提取”,中国图像图形学报,Vol.12,No.4,pp.655-660,2007)提出的方法虽然能够得到良好的分割结果,但依赖于初始轮廓线的选取,尤其对于边缘比较模糊、纹理自相似性较高的IVUS图像,其初始轮廓更不易确定,从而也在一定程度上影响了这种分割方法的准确性。
发明内容
针对现有技术中存在的上述问题,本发明提供一种基于血管内超声图像的血管ROI分割方法,采用改进的水平集模型算法结合窄带法首先分割血管的管腔区域以及管腔膜轮廓,通过定位管腔区域的中心以获得参数主动轮廓模型(Snake)的初始轮廓,继而通过收敛得到血管的中外膜轮廓曲线,中外膜的提取充分利用了管腔区域信息的先验知识。最后,将中外膜轮廓曲线以内区域作为ROI,通过全局最小化活动轮廓模型实现对血管斑块的分割。
基于血管内超声图像的血管ROI分割方法,是以血管内超声图像为基础,根据IVUS图像的特点,对血管内超声图像依次进行血管管腔膜分割,血管中外膜分割,血管斑块分割。该方法能得到管腔区域轮廓信息,中外膜边缘信息以及斑块形状轮廓信息等,在不增加附加设备的前提下,充分利用血管内超声设备提供的超声图像本身的灰度信息以及区域特性并结合各种分割算法,实现血管ROI的分割。
本发明的特征如下:
步骤1,利用血管内超声仪,匀速回拉导管,获得人体冠状动脉的血管内超声视频影像。
步骤2,将步骤1得到的血管内超声视频影像导入计算机,从视频中截取连续的血管内超声图像序列作为实验图像,图像分辨率为2x0*2y0,以下简称为超声图像;
步骤3,对上述超声图像进行超声图像管腔膜的分割。
步骤3.1,创建与待分割IVUS图像大小相同的初始管腔轮廓模板矩阵。
设待分割图像为I(x,y),创建大小同样为2x0*2y0的空矩阵Minit,确立中心坐标为O(x0,y0),令Minit(x0-a:x0+b,y0-c:y0+d)=1,即创建了长为c+d宽为a+b的矩形初始化区域模板。
步骤3.2,初始化符号距离函数映射矩阵。
设Minit为步骤3.1中设置的初始管腔轮廓模板矩阵,计算Minit中当前像素点与最近的非零像素点的欧氏距离,单位像素距离为1,设符号距离函数映射矩阵为:
φ=bwdist(Minit)-bwdist(1-Minit)-0.5
其中,bwdist为距离变换函数,默认计算二值图中当前像素点与最近的非0像素点的距离,并返回与原二值图相同大小的结果矩阵。
步骤3.3,设置窄带的宽度。
得到符号距离函数映射矩阵φ之后,设置窄带Wnarrowband的宽度为2k,窄带Wnarrowband为满足-k≤φ≤k条件的像素点组成的集合。
步骤3.4,计算窄带区域内的图像数据力Fimage,公式如下:
Fimage=[IWnarrowband(x,y)-ca]2-[IWnarrowband(x,y)-cb]2
其中ca、cb分别是待分割图像I(x,y)内外两个区域上的灰度均值。
步骤3.5,计算曲率速度项Vcurve。
对于窄带区域的点φ(x,y)(-k≤φ≤k),通过八邻域的方法分别构造φx、φy、φxx、φyy和φxy:
φx=φ(x+1,y)-φ(x-1,y)
φy=φ(x,y+1)-φ(x,y-1)
φxx=φ(x-1,y)+φ(x+1,y)-2φ(x,y)
φyy=φ(x,y-1)+φ(x,y+1)-2φ(x,y)
φxy=-0.25*φ(x-1,y-1)-0.25*φ(x+1,y+1)
+0.25*φ(x+1,y-1)+0.25*φ(x-1,y+1)
曲率速度项Vcurve为:
步骤3.6,计算速度函数,公式如下:
其中,α为光顺权重系数。
步骤3.7,对窄带Wnarrowband内符号距离函数映射矩阵φ进行迭代,迭代公式为:
式中,n为迭代次数,Δt1=0.45/max(F)为时间步长。
步骤3.8,经过n次迭代,最终符号距离函数映射矩阵φ的迭代结果即为血管管腔区域的二值矩阵,记做管腔区域轮廓图像为finit,提取管腔膜边缘即得到管腔膜分割结果。
步骤4,实现超声图像的中外膜分割。
步骤4.1,确定管腔轮廓区域图像的中心。
计算由步骤3中得到的血管管腔区域轮廓图像finit的p+q阶矩mpq,公式如下:
其中,p和q都是非负整数。
管腔轮廓区域图像的中心Ccentre(xc,yc)的坐标为:
其中,m00表示图像finit的0阶矩,m10和m01表示图像finit的一阶矩。
步骤4.2,初始化Snake初始轮廓。
设Snake初始轮廓曲线为C0(s),其坐标点为(x(s),y(s)),令s∈[0,1],θ∈[0,2π]。
步骤4.3,求Snake轮廓曲线满足的方程。
定义C(s)=(x(s),y(s))的能量:
其中,为图像能量,I(x,y)为待分割图像。
对s求导得:
设空间步长为h,轮廓线上点的坐标为Xi=(xi,yi),将导数用差分近似:
将上式化为Xi的矢量形式的线性方程组:
其中,A为由α,β,h构成的对称五对角循环矩阵,(x,y)为轮廓线上点的横坐标和纵坐标组成的矢量,fx(x,y),fy(x,y)为由Eext在Xi处对x,y的偏导数组成的矢量。
步骤4.4,通过迭代求解血管的中外膜轮廓。
将轮廓C(s)视作s和时间t的函数C(s,t),初值为:C(s,0)=C0(s),迭代求解结果为:
选取初始的Snake曲线作为方程的初解,当解收敛时,即Ci(s,t)≈0时得到最优曲线,该曲线即为血管的中外膜轮廓。
步骤5,利用全局最小化主动轮廓模型算法对血管斑块部分进行分割。
将步骤4得到的超声图像中外膜边缘作为超声图像ROI的临界线,选取步骤4中所分割的中外膜轮廓线以内的区域作为斑块分割的ROI,利用全局最小化主动轮廓模型算法对血管斑块部分进行分割。
步骤5.1,构造需要最小化的泛函。
基于全局最小化活动轮廓模型思想,为了得到血管内超声图像ROI的斑块轮廓,构造需要最小化的泛函:
E1(u=ψΩc,c1,c2,λ)=∫cgds+λ∫Ω((c1-f(x))2-(c2-f(x))2ψΩcdx
=TVg(ψΩc)+λ∫Ωr1(x,c1,c2)ψΩcdx
其中,TVg(u)是带有权重函数g(x)的函数u的总变分函数,g(x)是边缘检测函数,函数u是特征函数ψΩc。
步骤5.2,构造区域信息项,表达式为:
r1(x,c1,c2)=(fROI(x)-c1)2-(fROI(x)-c2)2
其中,c1为图像I(x)中u≥0.5区域的灰度均值,c2为图像I(x)中u<0.5区域的灰度均值,fROI(x)为待分割的血管内超声ROI图像。
步骤5.3,正则化泛函E1(u=ψΩc,c1,c2,λ)得:
步骤5.4,固定v,搜索u作为的解:
u=v-θdivp
其中,p=(p1,p2),由不动点方法解方程得到:
p0=0,且
其中,
步骤5.5,固定u,搜索v作为的解:
v=min{max{u(x)-θλr1(x,c1,c2),0},1}
步骤5.6,通过迭代完成斑块分割。
初始化u0=0,p0=0,设定最大迭代次数imax=N,根据步骤5.2、5.4和5.5分别计算pi,ui,vi,当迭代次数达到设定的次数N时停止迭代,保存当前变量vN即斑块分割结果。
与现有技术相比,本发明具有以下有益效果:
本发明提出了一种基于血管内超声图像的ROI分割方法,实现了血管ROI的管腔膜中外膜以及斑块的轮廓信息可视化,与基于统计学的IVUS图像分割方法相比,本发明摒弃了其复杂的统计建模过程且分割结果不受IVUS图像伪影和斑块特征的影响;与传统的基于活动轮廓模型的IVUS图像分割方法相比,本发明只需通过设置简单的矩形初始轮廓模板即可依次得到管腔区域轮廓从而定位Snake初始轮廓以进行中外膜的分割,省掉了对IVUS图像中外膜边缘进行初始轮廓的预分割这一步骤,提高了分割效率。本发明引入了血管斑块的分割,这是目前IVUS图像分割方法未涉及到的。
附图说明
图1是IVUS图像血管ROI分割组图:(a)为待分割的IVUS图像,(b)为血管管腔膜分割结果,(c)为二值化的管腔轮廓区域,(d)为血管中外膜分割结果,(e)为血管ROI图像,(f)为血管ROI的斑块轮廓显示;
图2是连续9帧血管内超声图像管腔膜分割效果图;
图3是连续9帧血管内超声图像中外膜分割效果图;
图4是连续9帧血管内超声图像ROI斑块轮廓图;
图5是本发明方法的流程图。
具体实施方式
本发明是采用以下技术手段实现的:
一种基于血管内超声图像的血管ROI分割方法。首先,对待分割图像结合改进的水平集模型算法和窄带法实现对血管内超声图像管腔膜的分割,然后,定位管腔区域的中心并生成Snake模型的初始轮廓曲线进行迭代收敛从而得到血管中外膜的边缘,最后,选取中外膜以内区域为斑块的ROI,结合全局最小化主动轮廓模型的算法实现对血管ROI的斑块轮廓的分割。
上述基于血管内超声图像的血管ROI分割方法,包括下述步骤:
步骤1,利用血管内超声仪,以0.5mm/s的速度匀速回拉导管,获得人体冠状动脉的血管内超声视频影像;
步骤2,将步骤(1)得到的血管内超声视频影像导入计算机,从视频中截取连续的血管内超声图像序列作为实验图像,图像分辨率为384*384,以下简称为超声图像;
步骤3,对超声图像进行血管管腔膜的分割;
设待分割图像为I(x,y),创建大小同样为384*384的空矩阵Minit,确立中心坐标为O(192,192),令Minit(192-a:192+b,192-c:192+d)=1,即创建了长为c+d宽为a+b的矩形初始化区域模板。
设Minit为初始轮廓模板矩阵,计算Minit中当前像素点与最近的非零像素点的欧氏距离(单位像素距离为1),设符号距离函数映射矩阵:
φ=bwdist(Minit)-bwdist(1-Minit)-0.5 (1)
其中,bwdist为距离变换函数,默认计算二值图中当前像素点与最近的非0像素点的距离,并返回与原二值图相同大小的结果矩阵。
设置窄带的宽度为2k,设窄带Wnarrowband为满足-k≤φ≤k条件的像素点组成的集合。设点集Pinner为满足条件(φ≤0)的像素点的集合,点集Pouter为满足条件(φ>0)的像素点的集合,并计算待分割图像I(x,y)在Pinner以及Pouter区域的灰度均值。
构造窄带区域内的图像数据力Fimage:
Fimage=[IWnarrowband(x,y)-ca]2-[IWnarrowband(x,y)-cb]2 (2)
对于窄带区域的点φ(x,y)(-k≤φ≤k),采用八邻域的方法分别构造φx,φy,φxx,φyy,φxy
φx=φ(x+1,y)-φ(x-1,y) (3)
φy=φ(x,y+1)-φ(x,y-1) (4)
φxx=φ(x-1,y)+φ(x+1,y)-2φ(x,y) (5)
φyy=φ(x,y-1)+φ(x,y+1)-2φ(x,y) (6)
φxy=-0.25*φ(x-1,y-1)-0.25*φ(x+1,y+1)+0.25*φ(x+1,y-1)+0.25*φ(x-1,y+1) (7)
通过下式计算得到曲率速度项:
计算改进的水平集速度函数F:
其中光顺权重系数α为常数,通过窄带Wnarrowband内的符号距离函数迭代公式进行符号距离函数的迭代,其中n为迭代次数,时间步长设置为Δt1=0.45/max(F)。经过n次迭代,最终迭代结果φ即为血管管腔区域的二值矩阵,提取管腔膜边缘,记该二值矩阵为finit,如图1(c)所示。
步骤4,实现超声图像中外膜的分割;
对于步骤3中所获得的血管管腔区域轮廓图像finit,假设:
mpq为图像finit的p+q阶矩,其中p和q都是非负整数,血管管腔的中心Ccentre(xc,yc)的坐标按下式计算:
通过下式得到Snake初始轮廓曲线C0(s)上的点的坐标(xsnake,ysnake):
其中,θ∈[0,2π]。
令其中s∈[0,1]。定义曲线C(s)=(x(s),y(s))的能量为:
图像能量Eext为:
其中I(x,y)为待分割图像。对上式求导得到:
设空间步长为h,轮廓线上点的坐标为Xi=(xi,yi),将导数用差分近似得:
上式可以化为如下Xi的矢量形式的线性方程组:
其中A为由α,β,h构成的对称五对角循环矩阵,(x,y)为轮廓线上点的横坐标和纵坐标组成的矢量,fx(x,y),fy(x,y)为由Eext在Xi处对x,y的偏导数组成的矢量。
将轮廓C视作s和时间t的函数C(s,t),则:
迭代求解结果为:
取初值C(s,0)=C0(s)作为方程的初解,经过N次迭代后当解收敛时即Ci(s,t)≈0时得到最优曲线,该曲线即为血管的中外膜轮廓,轮廓线上点的坐标为(xfin,yfin)。
步骤5,载入超声图像的ROI图像,如图1(e)所示利用全局最小化主动轮廓模型算法,分割血管ROI斑块轮廓。
由步骤4所获得的中外膜轮廓曲线,选取中外膜轮廓曲线内部区域为血管ROI,得到超声图像的血管ROI图像fROI。为得到血管内超声图像ROI的斑块轮廓,基于全局最小化活动轮廓模型思想构造需要最小化的泛函
其中TVg(u)是带有权重函数g(x)的函数u的总变分函数g(x)是边缘检测函数,函数u是特征函数ψΩc。
构造区域信息项:
r1(x,c1,c2)=(fROI(x)-c1)2-(fROI(x)-c2)2 (22)
其中c1为图像I(x)中u≥0.5区域灰度均值,c2为图像I(x)中u<0.5区域灰度均值。fROI(x)为待分割的血管ROI图像。
正则化E1(u=ψΩc,c1,c2,λ)得下式:
通过固定v搜索u作为的解,解得:
u=v-θdivp (24)
其中p=(p1,p2)由下式给出,
上式可由不动点方法解出:
其中p0=0,
通过固定u搜索v作为的解,解得:
v=min{max{u(x)-θλr1(x,c1,c2),0},1} (27)
在斑块分割的迭代过程中,初始化u0=0,p0=0,设置迭代次数i,根据式(26)(24)和(22)(27)分别计算pi,ui,vi,当迭代次数达到设定的次数N时,停止迭代,保存当前变量vN即斑块分割结果,如图1(f)所示。
Claims (2)
1.一种基于血管内超声图像的血管ROI分割方法,其特征在于,包括以下步骤:
步骤1,利用血管内超声仪,匀速回拉导管,获得人体冠状动脉的血管内超声视频影像;
步骤2,将步骤1得到的血管内超声视频影像导入计算机,从视频中截取连续的血管内超声图像序列作为实验图像,图像分辨率为2x0*2y0,以下简称为超声图像;
步骤3,对上述超声图像进行超声图像管腔膜的分割;
步骤3.1,创建与待分割血管内超声图像大小相同的初始管腔轮廓模板矩阵;
设待分割图像为I(x,y),创建大小同样为2x0*2y0的空矩阵Minit,确立中心坐标为O(x0,y0),令Minit(x0-a:x0+b,y0-c:y0+d)=1,即创建了长为c+d宽为a+b的矩形初始化区域模板;
步骤3.2,初始化符号距离函数映射矩阵;
设Minit为步骤3.1中设置的初始管腔轮廓模板矩阵,计算Minit中当前像素点与最近的非零像素点的欧氏距离,单位像素距离为1,设符号距离函数映射矩阵为:
φ=bwdist(Minit)-bwdist(1-Minit)-0.5
其中,bwdist为距离变换函数,默认计算二值图中当前像素点与最近的非0像素点的距离,并返回与原二值图相同大小的结果矩阵;
步骤3.3,设置窄带的宽度;
得到符号距离函数映射矩阵φ之后,设置窄带Wnarrowband的宽度为2k,窄带Wnarrowband为满足-k≤φ≤k条件的像素点组成的集合;
步骤3.4,计算窄带区域内的图像数据力Fimage,公式如下:
Fimage=[IWnarrowband(x,y)-ca]2-[IWnarrowband(x,y)-cb]2
其中ca、cb分别是待分割图像I(x,y)内外两个区域上的灰度均值;
步骤3.5,计算曲率速度项Vcurve;
对于窄带区域的点φ(x,y),-k≤φ≤k,通过八邻域的方法分别构造φx、φy、φxx、φyy和φxy:
φx=φ(x+1,y)-φ(x-1,y)
φy=φ(x,y+1)-φ(x,y-1)
φxx=φ(x-1,y)+φ(x+1,y)-2φ(x,y)
φyy=φ(x,y-1)+φ(x,y+1)-2φ(x,y)
φxy=-0.25*φ(x-1,y-1)-0.25*φ(x+1,y+1)
+0.25*φ(x+1,y-1)+0.25*φ(x-1,y+1)
曲率速度项Vcurve为:
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>c</mi>
<mi>u</mi>
<mi>r</mi>
<mi>v</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
<msubsup>
<mi>&phi;</mi>
<mi>y</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>&phi;</mi>
<mi>y</mi>
</msub>
<msub>
<mi>&phi;</mi>
<mi>x</mi>
</msub>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>y</mi>
<mi>y</mi>
</mrow>
</msub>
<msubsup>
<mi>&phi;</mi>
<mi>x</mi>
<mn>2</mn>
</msubsup>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&phi;</mi>
<mi>x</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&phi;</mi>
<mi>y</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mfrac>
<mn>3</mn>
<mn>2</mn>
</mfrac>
</msup>
</mfrac>
</mrow>
步骤3.6,计算速度函数,公式如下:
<mrow>
<mi>F</mi>
<mo>=</mo>
<mfrac>
<msub>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
</mrow>
</msub>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mo>|</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>|</mo>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>+</mo>
<mi>&alpha;</mi>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
<msubsup>
<mi>&phi;</mi>
<mi>y</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>&phi;</mi>
<mi>y</mi>
</msub>
<msub>
<mi>&phi;</mi>
<mi>x</mi>
</msub>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>y</mi>
<mi>y</mi>
</mrow>
</msub>
<msubsup>
<mi>&phi;</mi>
<mi>x</mi>
<mn>2</mn>
</msubsup>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&phi;</mi>
<mi>x</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&phi;</mi>
<mi>y</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mfrac>
<mn>3</mn>
<mn>2</mn>
</mfrac>
</msup>
</mfrac>
</mrow>
其中,α为光顺权重系数;
步骤3.7,对窄带Wnarrowband内符号距离函数映射矩阵φ进行迭代,迭代公式为:
<mrow>
<msubsup>
<mi>&phi;</mi>
<mi>W</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>&phi;</mi>
<mi>W</mi>
<mi>n</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&Delta;t</mi>
<mn>1</mn>
</msub>
<mo>&CenterDot;</mo>
<mi>F</mi>
</mrow>
式中,n为迭代次数,Δt1=0.45/max(F)为时间步长;
步骤3.8,经过N次迭代,N表示最大迭代次数,最终符号距离函数映射矩阵φ的迭代结果即为血管管腔区域的二值矩阵,记做管腔区域轮廓图像为finit,提取管腔膜边缘即得到管腔膜分割结果;
步骤4,实现超声图像的中外膜分割;
步骤4.1,确定管腔轮廓区域图像的中心;
计算由步骤3中得到的血管管腔区域轮廓图像finit的p+q阶矩mpq,公式如下:
<mrow>
<msub>
<mi>m</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mi>i</mi>
<mi>p</mi>
</msup>
<msup>
<mi>j</mi>
<mi>q</mi>
</msup>
<msub>
<mi>f</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
其中,p和q都是非负整数;
管腔轮廓区域图像的中心Ccentre(xc,yc)的坐标为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>m</mi>
<mn>10</mn>
</msub>
<msub>
<mi>m</mi>
<mn>00</mn>
</msub>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>m</mi>
<mn>01</mn>
</msub>
<msub>
<mi>m</mi>
<mn>00</mn>
</msub>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,m00表示图像finit的0阶矩,m10和m01表示图像finit的一阶矩;
步骤4.2,初始化Snake初始轮廓;
设Snake初始轮廓曲线为C0(s),其坐标点为(x(s),y(s)),令s∈[0,1],θ∈[0,2π];
步骤4.3,求Snake轮廓曲线满足的方程;
定义C(s)=(x(s),y(s))的能量:
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>s</mi>
<mi>n</mi>
<mi>a</mi>
<mi>k</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<mn>1</mn>
</munderover>
<mo>{</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>&lsqb;</mo>
<mi>&alpha;</mi>
<mo>|</mo>
<msup>
<mi>C</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>+</mo>
<mi>&beta;</mi>
<mo>|</mo>
<msup>
<mi>C</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>+</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>e</mi>
<mi>x</mi>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>}</mo>
<mi>d</mi>
<mi>s</mi>
</mrow>
其中,为图像能量,I(x,y)为待分割图像;
对s求导得:
<mrow>
<msup>
<mi>&alpha;C</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<mi>&beta;C</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>&dtri;</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>e</mi>
<mi>x</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
设空间步长为h,轮廓线上点的坐标为Xi=(xi,yi),将导数用差分近似,
<mrow>
<msubsup>
<mi>X</mi>
<mi>i</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>,</mo>
<msubsup>
<mi>X</mi>
<mi>i</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mo>-</mo>
<mn>2</mn>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mo>+</mo>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mrow>
<msup>
<mi>h</mi>
<mn>4</mn>
</msup>
</mfrac>
</mrow>
将上式化为Xi的矢量形式的线性方程组:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mi>A</mi>
<mi>x</mi>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>x</mi>
</msub>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
<mo>=</mo>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>A</mi>
<mi>y</mi>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>y</mi>
</msub>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
<mo>=</mo>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,A为由α,β,h构成的对称五对角循环矩阵,(x,y)为轮廓线上点的横坐标和纵坐标组成的矢量,fx(x,y),fy(x,y)为由Eext在Xi处对x,y的偏导数组成的矢量;
步骤4.4,通过迭代求解血管的中外膜轮廓;
将轮廓C(s)视作s和时间t的函数C(s,t),初值为:C(s,0)=C0(s),迭代求解结果为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mo>+</mo>
<mi>&gamma;</mi>
<mi>I</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>&gamma;x</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>x</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mo>+</mo>
<mi>&gamma;</mi>
<mi>I</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&gamma;y</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>y</mi>
</msub>
<mo>(</mo>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
选取初始的Snake曲线作为方程的初解,当解收敛时,即Ci(s,t)≈0时得到最优曲线,该曲线即为血管的中外膜轮廓;
步骤5,利用全局最小化主动轮廓模型算法对血管斑块部分进行分割;
将步骤4得到的超声图像中外膜边缘作为超声图像ROI的临界线,选取步骤4中所分割的中外膜轮廓线以内的区域作为斑块分割的ROI,利用全局最小化主动轮廓模型算法对血管斑块部分进行分割。
2.根据权利要求1所述的一种基于血管内超声图像的血管ROI分割方法,其特征在于,所述步骤5利用全局最小化主动轮廓模型算法对血管斑块部分进行分割的方法包括以下步骤:
步骤5.1,构造需要最小化的泛函;
基于全局最小化活动轮廓模型思想,为了得到血管内超声图像ROI的斑块轮廓,构造需要最小化的泛函:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>u</mi>
<mo>=</mo>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>&Omega;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mi>&lambda;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mo>&Integral;</mo>
<mi>c</mi>
</msub>
<mi>g</mi>
<mi>d</mi>
<mi>s</mi>
<mo>+</mo>
<mi>&lambda;</mi>
<msub>
<mo>&Integral;</mo>
<mi>&Omega;</mi>
</msub>
<mo>(</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>&Omega;</mi>
<mi>c</mi>
</mrow>
</msub>
<mi>d</mi>
<mi>x</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msub>
<mi>TV</mi>
<mi>g</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>&Omega;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&lambda;</mi>
<msub>
<mo>&Integral;</mo>
<mi>&Omega;</mi>
</msub>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>&Omega;</mi>
<mi>c</mi>
</mrow>
</msub>
<mi>d</mi>
<mi>x</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,TVg(u)是带有权重函数g(x)的函数u的总变分函数,g(x)是边缘检测函数,函数u是特征函数ψΩc;
步骤5.2,构造区域信息项,表达式为:
r1(x,c1,c2)=(fROI(x)-c1)2-(fROI(x)-c2)2
其中,c1为图像I(x)中u≥0.5区域的灰度均值,c2为图像I(x)中u<0.5区域的灰度均值,fROI(x)为待分割的血管内超声ROI图像;
步骤5.3,正则化泛函E1(u=ψΩc,c1,c2,λ)得:
<mrow>
<msubsup>
<mi>E</mi>
<mn>2</mn>
<mi>r</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>,</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mi>&lambda;</mi>
<mo>,</mo>
<mi>&alpha;</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>TV</mi>
<mi>g</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mi>&theta;</mi>
</mrow>
</mfrac>
<mo>|</mo>
<mo>|</mo>
<mi>u</mi>
<mo>-</mo>
<mi>v</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mo>&Integral;</mo>
<mi>&Omega;</mi>
</msub>
<msub>
<mi>&lambda;r</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>v</mi>
<mo>+</mo>
<mi>&alpha;</mi>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>x</mi>
</mrow>
步骤5.4,固定v,搜索u作为的解:
u=v-θdivp
其中,p=(p1,p2),由不动点方法解方程得到:
p0=0,且
其中,
步骤5.5,固定u,搜索v作为的解:
v=min{max{u(x)-θλr1(x,c1,c2),0},1}
步骤5.6,通过迭代完成斑块分割;
初始化u0=0,p0=0,设定最大迭代次数imax=N,根据步骤5.2、5.4和5.5分别计算pi,ui,vi,当迭代次数达到设定的次数N时停止迭代,保存当前变量vN即斑块分割结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410115141.8A CN103886599B (zh) | 2014-03-26 | 2014-03-26 | 一种基于血管内超声图像的血管roi分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410115141.8A CN103886599B (zh) | 2014-03-26 | 2014-03-26 | 一种基于血管内超声图像的血管roi分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103886599A CN103886599A (zh) | 2014-06-25 |
CN103886599B true CN103886599B (zh) | 2017-10-13 |
Family
ID=50955472
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410115141.8A Expired - Fee Related CN103886599B (zh) | 2014-03-26 | 2014-03-26 | 一种基于血管内超声图像的血管roi分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103886599B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104537645A (zh) * | 2014-12-15 | 2015-04-22 | 北京工业大学 | 一种基于血管内超声图像的roi标记点匹配方法 |
CN106388867A (zh) * | 2016-09-28 | 2017-02-15 | 深圳华声医疗技术有限公司 | 血管的内中膜自动识别测量方法及超声仪 |
WO2018218479A1 (en) * | 2017-05-31 | 2018-12-06 | Edan Instruments, Inc. | Systems and methods for adaptive enhancement of vascular imaging |
CN107909581B (zh) * | 2017-11-03 | 2019-01-29 | 杭州依图医疗技术有限公司 | Ct影像的肺叶段分割方法、装置、***、存储介质及设备 |
CN107909585B (zh) * | 2017-11-14 | 2020-02-18 | 华南理工大学 | 一种血管内超声影像的血管中内膜分割方法 |
CN109166132B (zh) * | 2018-07-16 | 2022-01-07 | 哈尔滨工程大学 | 一种可变初始距离符号函数的侧扫声呐图像目标识别方法 |
CN109674493B (zh) * | 2018-11-28 | 2021-08-03 | 深圳蓝韵医学影像有限公司 | 医用超声自动追踪颈动脉血管的方法、***及设备 |
CN110033466B (zh) * | 2019-04-01 | 2020-12-18 | 数坤(北京)网络科技有限公司 | 一种基于多灰度级的冠脉拉直图像分割边界确定方法 |
CN110246136B (zh) * | 2019-05-29 | 2021-07-02 | 山东大学 | 一种基于混合算法的血管内超声参数提取方法及*** |
CN112258533B (zh) * | 2020-10-26 | 2024-02-02 | 大连理工大学 | 一种超声图像中小脑蚓部的分割方法 |
CN112826535B (zh) * | 2020-12-31 | 2022-09-09 | 青岛海信医疗设备股份有限公司 | 一种超声成像中自动定位血管的方法和装置及设备 |
CN113222956B (zh) * | 2021-05-25 | 2023-09-15 | 南京大学 | 基于超声图像识别血管中斑块的方法 |
CN114648514B (zh) * | 2022-03-30 | 2022-11-29 | 中国人民解放军总医院第二医学中心 | 一种脑动脉定位提取方法、装置、电子设备及存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102800089A (zh) * | 2012-06-28 | 2012-11-28 | 华中科技大学 | 基于颈部超声图像的主颈动脉血管提取和厚度测量方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7397935B2 (en) * | 2004-05-10 | 2008-07-08 | Mediguide Ltd. | Method for segmentation of IVUS image sequences |
-
2014
- 2014-03-26 CN CN201410115141.8A patent/CN103886599B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102800089A (zh) * | 2012-06-28 | 2012-11-28 | 华中科技大学 | 基于颈部超声图像的主颈动脉血管提取和厚度测量方法 |
Non-Patent Citations (2)
Title |
---|
Intravascular Ultrasound Image Hard Plaque Recognition and Media-adventitia Border Detection;XING Dong等;《Chinese Journal of Biomedical Engineering》;20120930(第03期);第110-116页 * |
血管内超声图像的预处理;汪友生等;《国外电子测量技术》;20120215;第31卷(第2期);第66-70页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103886599A (zh) | 2014-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103886599B (zh) | 一种基于血管内超声图像的血管roi分割方法 | |
Qian et al. | An integrated method for atherosclerotic carotid plaque segmentation in ultrasound image | |
CN104504737B (zh) | 一种从肺部ct图像获得三维气管树的方法 | |
CN112258514B (zh) | 一种ct影像肺血管的分割方法 | |
JP2011514190A (ja) | 画像内の解剖学的構造のセグメンテーションおよびモデリングのための方法およびシステム | |
US7359538B2 (en) | Detection and analysis of lesions in contact with a structural boundary | |
Zhou et al. | Skin lesion segmentation using an improved snake model | |
WO2022105623A1 (zh) | 一种基于迁移学习的颅内血管病灶识别方法 | |
Shang et al. | Region competition based active contour for medical object extraction | |
Raj et al. | Automated liver tumor detection using markov random field segmentation | |
CN107392928A (zh) | 基于MRF‑Snake分割的乳腺癌诊断***及方法 | |
Chung et al. | Accurate liver vessel segmentation via active contour model with dense vessel candidates | |
CN113935976A (zh) | 一种增强ct图像的脏器内血管自动分割方法及*** | |
Chen et al. | Liver segmentation from CT images based on region growing method | |
CN113470060B (zh) | 基于ct影像的冠状动脉多角度曲面重建可视化方法 | |
Loizou et al. | Segmentation of atherosclerotic carotid plaque in ultrasound video | |
Bai et al. | Automatic whole heart segmentation based on watershed and active contour model in CT images | |
Marin-Castrillon et al. | 4D segmentation of the thoracic aorta from 4D flow MRI using deep learning | |
Santos et al. | Automatic segmentation of the lumen of the carotid artery in ultrasound B-mode images | |
Guo et al. | Retinal vessel segmentation using a finite element based binary level set method | |
Alom et al. | Automatic slice growing method based 3D reconstruction of liver with its vessels | |
Spuhler et al. | Fast and robust extraction of centerlines in 3D tubular structures using a scattered-snakelet approach | |
Makrogiannis et al. | Nonparametric intensity priors for level set segmentation of low contrast structures | |
Han et al. | Detection of colonic polyp candidates with level set-based thickness mapping over the colon wall | |
Gu et al. | Novel multistage three-dimensional medical image segmentation: Methodology and validation |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171013 Termination date: 20200326 |