发明名称 一种自适应的医疗超声系统成像中的声学参数优化的方法
摘要 本发明公开了一种自适应的医疗超声系统成像参数的优化方法,设计了一种通过优化声速和组织声衰减系数来提高超声成像质量的方法,该方法通过迭代改变声速和组织声衰减系数,并使用一系列处理方法来获得最优的扫描声速和组织声衰减系数,最后把获得系数设置为该次超声成像所使用的参数,获得亮度更均匀,分辨率更好,对比度更强的超声影像。利用自适应的医疗超声系统成像参数的优化技术对诊断具有重要临床意义,适用于医学技术领域。
申请公布号 CN102697521A 申请公布日期 2012.10.03
申请号 CN201210124907.X 申请日期 2012.06.01
申请人 声泰特(成都)科技有限公司 发明人 刘东权;高君
分类号 A61B8/00(2006.01)I 主分类号 A61B8/00(2006.01)I
代理机构 代理人
主权项 一种自适应的医疗超声系统成像中的声学参数优化方法,其特征在于,包括如下步骤:1)在声衰减优化器中使用不同的组织声衰减系数补偿同一幅二维灰度图像,声衰减系数的迭代范围为0.3dB/cmMHz~1.3dB/cmMHz,从而获得多幅使用不同声衰减系数得到的二维灰度图像;2)以小窗口遍历每幅二维灰度图像,根据如下公式计算每个小窗口的纹理能量值,根据各个小窗口的纹理能量值计算每张灰度图像的平均纹理能量值; <mrow> <msub> <mi>P</mi> <mi>s</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>;</mo> <msub> <mi>d</mi> <mi>x</mi> </msub> <mo>,</mo> <msub> <mi>d</mi> <mi>y</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mi>Card</mi> <mo>{</mo> <mrow> <mo>(</mo> <mi>k</mi> <mo>,</mo> <mi>l</mi> <mo>)</mo> </mrow> <mo>|</mo> <msub> <mi>y</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>l</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>y</mi> <mrow> <mi>k</mi> <mo>+</mo> <msub> <mi>d</mi> <mi>x</mi> </msub> <mo>,</mo> <mi>l</mi> <mo>+</mo> <msub> <mi>d</mi> <mi>y</mi> </msub> </mrow> </msub> <mo>=</mo> <mi>i</mi> <mo>}</mo> </mrow> <mrow> <msub> <mi>P</mi> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>j</mi> <mo>;</mo> <msub> <mi>d</mi> <mi>x</mi> </msub> <mo>,</mo> <msub> <mi>d</mi> <mi>y</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mi>Card</mi> <mo>{</mo> <mrow> <mo>(</mo> <mi>k</mi> <mo>,</mo> <mi>l</mi> <mo>)</mo> </mrow> <mo>|</mo> <msub> <mi>y</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>l</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>y</mi> <mrow> <mi>k</mi> <mo>+</mo> <msub> <mi>d</mi> <mi>x</mi> </msub> <mo>,</mo> <mi>l</mi> <mo>+</mo> <msub> <mi>d</mi> <mi>y</mi> </msub> </mrow> </msub> <mo>=</mo> <mi>j</mi> <mo>}</mo> </mrow>能量:∑iPs(i)2∑jPd(j)23)在所有二维灰度图像中,找到能量值最大的那幅灰度图对应的组织声衰减系数;4)把得到的组织声衰减系数设置为当前扫描的组织声衰减系数;5)在声速优化器中使用不同声速对同一扫描部位发射超声波,声速的迭代范围为1450m/s~1640m/s,从而获得多幅不同声速对应的二维灰度图像;6)图像分析:a)以小窗口遍历每幅二维灰度图像,根据如下公式计算每个小窗口的纹理熵值,并根据各个小窗口的纹理熵值计算每幅二维灰度图像的平均纹理熵值熵:‑∑iPs(i)logPs‑∑jPd(j)logPdb)选出熵值最大的Z(>5)幅二维图像;c)根据如下公式计算选出的Z幅二维灰度的侧向对比度; <mrow> <msub> <mi>d</mi> <mi>i</mi> </msub> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>M</mi> <mo>-</mo> <mn>2</mn> </mrow> </munderover> <mo>|</mo> <msub> <mi>I</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>I</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>|</mo> </mrow> <mrow> <msub> <mi>&mu;</mi> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mi>M</mi> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>M</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>I</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> </mrow> <mrow> <mi>Contrast</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mfrac> <msub> <mi>d</mi> <mi>i</mi> </msub> <msub> <mi>&mu;</mi> <mi>i</mi> </msub> </mfrac> </mrow>此处,N为从图像中沿探头侧向方向选取的长度为M的整条采样线的数目,di、μi分别 代表第i条采样线所有像素点的逐点灰度差异绝对值总和与平均灰度值;7)找出Z幅二维灰度图像的侧向对比度最大的那幅灰度图像对应的声速;8)把得到声速设置为当前扫描的使用的系统声速。
地址 610041 四川省成都市高新区创业路16号火炬大厦A601