[发明专利]超声颈动脉血管膜的自动分割方法有效
申请号: | 201210216191.6 | 申请日: | 2012-06-28 |
公开(公告)号: | CN102800087A | 公开(公告)日: | 2012-11-28 |
发明(设计)人: | 丁明跃;李鹤;程洁玉;杨鑫 | 申请(专利权)人: | 华中科技大学 |
主分类号: | G06T7/00 | 分类号: | G06T7/00 |
代理公司: | 华中科技大学专利中心 42201 | 代理人: | 李智 |
地址: | 430074 湖北*** | 国省代码: | 湖北;42 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明属于计算机技术与医学图像的交叉领域,具体涉及一种超声图像颈动脉血管横截面方向上的内外轮廓自动分割方法。具体步骤如下:选取当前帧图像;人工定义基准轮廓点并插值形成基准轮廓;将基准轮廓外推形成第一感兴趣区域;自动检测颈动脉轮廓:包括在第一感兴趣区域转化得到的极坐标图像上寻找径向梯度最大值点作为初始外轮廓点、利用混合分布估计第一感兴趣区域内的灰度概率密度分布、水平集演化得到外轮廓、将外轮廓内推形成初始内轮廓、水平集演化得到内轮廓。本发明能对超声颈动脉图像进行有效的分割,减少医生工作量;同时基于本方法分割得到的颈动脉内外轮廓可用于计算颈动脉血管壁体积,为颈动脉粥样硬化病变研究提供准确有效的信息。 | ||
搜索关键词: | 超声 颈动脉 血管 自动 分割 方法 | ||
【主权项】:
1.超声颈动脉血管膜的自动分割方法,包括以下步骤:(1)若当前帧图像为颈动脉三维超声体数据的第一帧图像,则在当前帧图像上选择明显位于血管外轮廓上的像素点集作为基准点,通过基准点间插值形成闭合曲线,将其作为当前帧的基准轮廓;否则,将上一帧得到的备用基准轮廓作为当前帧的基准轮廓;(2)利用形态学膨胀法将基准轮廓外推形成第一感兴趣区域ROI1;(3)在第一感兴趣区域ROI1中检测颈动脉血管外轮廓CMAB:(3.1)在当前帧图像中提取包含第一感兴趣区域ROI1的矩形窗口图像,将其转化到极坐标图像,在极坐标图像上每隔相同角度在径向方向靠近基准轮廓的像素点中搜索梯度值最大者作为初始轮廓点,将初始轮廓点转化到直角坐标系并连线形成初始颈动脉血管外轮廓(3.2)利用混合分布估计第一感兴趣区域ROI1的灰度概率密度分布,从而得到混合分布的shape参数Kj和scale参数θj,j=1,…,M,M为混合分布中单项分布的类别数;(3.3)初始化迭代次数t=1以及符号距离函数表示第一感兴趣区域ROI1中初始颈动脉血管外轮廓以外的像素区域,x1为第一感兴趣区域ROI1中的像素点,||||为求欧式距离;(3.4)计算φ 1 t ( x 1 ) = φ 1 t - 1 ( x 1 ) - ρ 1 · δ ( φ 1 t - 1 ( x 1 ) ) · [ ln ( Σ j M q Ai ( j ) · G ( I x 1 | K j , θ j ) - ln ( Σ j M q Ae ( j ) · G ( I x 1 | K j , θ j ) ] ]]>+ λ 1 · δ ( φ 1 t - 1 ( x 1 ) ) · div ( ▿ φ 1 t - 1 ( x 1 ) | ▿ φ 1 t - 1 ( x 1 ) | ) ]]> 其中,这里0<ε<0.001,0<ρ1≤10,0<λ1≤20,qAi(j)和qAe(j)分别为前一次迭代演化轮廓的内、外区域中第j项单项分布的权重,为参数为Kj,θj的第j项单项分布函数,为像素点x1的灰度值;(3.5)计算本次迭代演化外轮廓C 1 t = { x 1 | φ 1 t ( x 1 ) = 0 } ]]> (3.6)若本次迭代演化外轮廓与前一次迭代得到的演化外轮廓相同或者t达到迭代次数上限,则本次迭代得到的演化轮廓即为颈动脉血管外轮廓CMAB,迭代结束,进入步骤(4);否则,t=t+1,返回步骤(3.4);(4)在当前帧图像中将颈动脉血管外轮廓CMAB围成的区域作为第二感兴趣区域ROI2;(5)在第二感兴趣区域ROI2中检测颈动脉血管内轮廓CLIB:(5.1)利用形态学腐蚀法将颈动脉血管外轮廓CMAB内推,将内推得到的轮廓线作为初始颈动脉血管内轮廓(5.2)初始化迭代次数t′=1以及符号距离函数表示第二感兴趣区域ROI2中初始颈动脉血管内轮廓以外的像素区域,x2为第二感兴趣区域ROI2中的像素点;(5.3)计算φ 2 t ′ ( x 2 ) = φ 2 t ′ - 1 ( x 2 ) - ρ 2 · δ ( φ 2 t ′ - 1 ( x 2 ) ) · [ ln ( Σ j M q Li ( j ) · G ( I x 3 | K j , θ j ) - ln ( Σ j M q Le ( j ) · G ( I x 3 | K j , θ j ) ] ]]>+ λ 2 · δ ( φ 2 t ′ - 1 ( x 2 ) ) · div ( ▿ φ 2 t ′ - 1 ( x 2 ) | ▿ φ 2 t ′ - 1 ( x 2 ) | ) + β 2 · δ ( φ 2 t ′ - 1 ( x 2 ) ) · B T ( x 2 ) ]]> 其中,0<ρ2≤10,0<λ2≤20,0<β2≤100,qLi(j)和qLe(j)分别为前一次演化轮廓的内、外区域中第j项单项分布的权重,为参数为Kj,θj的第j项单项分布函数,为像素点x2的灰度值;(5.4)计算本次迭代演化内轮廓C 2 t ′ = { x 2 | φ 2 t ′ ( x 2 ) = 0 } ]]> (5.5)若本次迭代演化内轮廓与前一次迭代得到的演化内轮廓相同或者达到迭代次数上限,则本次迭代得到的演化轮廓为颈动脉血管内轮廓CLIB,迭代结束,进入步骤(6);否则,t′=t′+1,返回步骤(5.3);(6)跟踪确定下一帧图像的备用基准轮廓:(6.1)在当前帧图像中,通过形态学膨胀法将颈动脉血管外轮廓CMAB外推,将外推得到的闭合轮廓线与CMAB构成的区域作为跟踪区域Ω;(6.2)在下一帧图像中搜索使得灰度差异和S=∑ΩΔI2最小的像素点集(x4,y4),其中ΔI2=[I(x3,y3)-J(x4,y4)]2,I(x3,y3)为当前帧图像的跟踪区域Ω中像素点(x3,y3)的灰度值,J(x4,y4)为下一帧图像的像素点(x4,y4)的灰度值;(6.3)令当前帧图像的跟踪区域中的点(x3,y3)与下一帧图像的像素点(x4,y4)的映射关系为x 4 y 4 = τ x τ y + 1 0 0 1 x 3 y 3 , ]]> 结合像素点(x3,y3)与(x4,y4)求解映射关系中的参数τ=(τx,τy),从而确定映射关系;(6.4)根据步骤(6.3)建立的映射关系将当前帧图像的颈动脉血管外轮廓CMAB映射为下一帧图像的备用基准轮廓。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于华中科技大学,未经华中科技大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201210216191.6/,转载请声明来源钻瓜专利网。
- 上一篇:试条、试条校正参数自动识别系统及血糖仪
- 下一篇:止动垫片