发明名称 一种超声病灶图像分割方法、装置及系统
摘要 本发明提出一种超声病灶图像分割方法,利用灰度信息和纹理信息建立一个总能量函数以构建图模型,然后对图模型求解以获得图像最优分割结果,从而将病灶区域从超声图像中快速分离出来。
申请公布号 CN104091331B 申请公布日期 2015.07.29
申请号 CN201410299718.5 申请日期 2014.06.27
申请人 深圳开立生物医疗科技股份有限公司 发明人 冯乃章;高梁
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 代理人
主权项 一种超声病灶图像分割方法,其特征在于,包括:获取包含病灶区域的超声子图像;初步将所述包含病灶区域的超声子图像划分为前景区域、背景区域和待分割区域;利用所述包含病灶区域的超声子图像的灰度和纹理信息建立一个总能量函数以构建图模型,通过最小化所述图模型的能量函数,以获得所述待分割区域的最优的图像分割结果;所述总能量函数为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>E</mi><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>i</mi><mo>&Element;</mo><mi>S</mi></mrow></munder><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><msub><mi>N</mi><mi>i</mi></msub></mrow></munder><msub><mi>E</mi><mi>Gray</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&times;</mo><msub><mi>E</mi><mi>Tex</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000670766010000011.GIF" wi="644" he="135" /></maths>其中,E<sub>Grag</sub>为基于灰度信息的能量函数,表示当前像素的邻域与相邻像素的邻域之间的灰度差异;E<sub>Tex</sub>为基于纹理信息的能量函数,表示当前像素的邻域与相邻像素的邻域之间的纹理差异;其中,S为所有像素的集合,i为当前像素,N<sub>i</sub>为像素i的4邻域,j为i的4邻域内的像素;灰度能量函数的定义为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>E</mi><mi>Gray</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mo>|</mo><mover><mi>T</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><mover><mi>T</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>|</mo><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000670766010000012.GIF" wi="678" he="80" /></maths>其中,<img file="FDA0000670766010000015.GIF" wi="103" he="86" />和<img file="FDA0000670766010000016.GIF" wi="108" he="75" />分别为相邻像素i和j局部邻域内的灰度均值;纹理能量函数的定义为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>E</mi><mi>Tex</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mrow><mo>(</mo><mo>|</mo><msub><mover><mi>H</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><msub><mover><mi>H</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>|</mo><mo>-</mo><mo>|</mo><msub><mover><mi>L</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><msub><mover><mi>L</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>|</mo><mo>)</mo></mrow><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000670766010000013.GIF" wi="1072" he="92" /></maths>其中,<img file="FDA0000670766010000017.GIF" wi="125" he="86" />和<img file="FDA0000670766010000018.GIF" wi="109" he="76" />分别为像素i局部邻域内像素灰度值高于和低于邻域灰度均值的像素的均值;<img file="FDA0000670766010000019.GIF" wi="142" he="80" />和<img file="FDA0000670766010000014.GIF" wi="112" he="59" />为像素j局部邻域内像素灰度值高于和低于邻域灰度均值的像素的均值。
地址 518051 广东省深圳市南山区玉泉路毅哲大厦4楼