主权项 |
一种超声病灶图像分割方法,其特征在于,包括:获取包含病灶区域的超声子图像;初步将所述包含病灶区域的超声子图像划分为前景区域、背景区域和待分割区域;利用所述包含病灶区域的超声子图像的灰度和纹理信息建立一个总能量函数以构建图模型,通过最小化所述图模型的能量函数,以获得所述待分割区域的最优的图像分割结果;所述总能量函数为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>E</mi><mo>=</mo><munder><mi>Σ</mi><mrow><mi>i</mi><mo>∈</mo><mi>S</mi></mrow></munder><munder><mi>Σ</mi><mrow><mi>j</mi><mo>∈</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>×</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>‾</mo></mover><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><mover><mi>T</mi><mo>‾</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>‾</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><msub><mover><mi>H</mi><mo>‾</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>‾</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><msub><mover><mi>L</mi><mo>‾</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局部邻域内像素灰度值高于和低于邻域灰度均值的像素的均值。 |