CN102697521A - 一种自适应的医疗超声系统成像中的声学参数优化的方法 - Google Patents
一种自适应的医疗超声系统成像中的声学参数优化的方法 Download PDFInfo
- Publication number
- CN102697521A CN102697521A CN201210124907XA CN201210124907A CN102697521A CN 102697521 A CN102697521 A CN 102697521A CN 201210124907X A CN201210124907X A CN 201210124907XA CN 201210124907 A CN201210124907 A CN 201210124907A CN 102697521 A CN102697521 A CN 102697521A
- Authority
- CN
- China
- Prior art keywords
- sound
- cloth
- width
- image
- velocity
- 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.)
- Pending
Links
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种自适应的医疗超声系统成像参数的优化方法,设计了一种通过优化声速和组织声衰减系数来提高超声成像质量的方法,该方法通过迭代改变声速和组织声衰减系数,并使用一系列处理方法来获得最优的扫描声速和组织声衰减系数,最后把获得系数设置为该次超声成像所使用的参数,获得亮度更均匀,分辨率更好,对比度更强的超声影像。利用自适应的医疗超声系统成像参数的优化技术对诊断具有重要临床意义,适用于医学技术领域。
Description
技术领域
本发明属于医学技术领域,涉及一种超声成像技术,具体涉及一种自适应的医疗超声系统成像中的声学参数优化方法。
背景技术
在医疗超声系统成像系统中,声速和组织声衰减系数是两个非常重要的成像参数。超声波在人体不同组织中传播时,声速和组织声衰减系数都差异较大,例如,超声波在人体脂肪含量较高的脂肪传播时,声速大约为1450m/s,组织声衰减系为0.63dB/cmMHz。而在肌肉中的声速为1580m/s左右,组织声衰减系为1.8dB/cmMHz。某些肌肉组织的声速甚至可高达1640m/s,如表1所示。现在的医疗超声系统通常把1540m/s用于动态延迟聚焦波束形成,然而使用声速偏离真实声速的超声波经过了局部组织时,就会出现相位偏差,导致组织位置的偏移,点扩散函数的拓宽以及旁瓣的增加,影响图像的动态范围,对比度,降低图像的质量,影响医疗诊断的准确性。
组织 | 声衰减系数(dB/cmMHz) | 声速(m/s) |
血液 | 0.2 | 1580 |
肌肉 | 1.8 | 1580 |
肾 | 1.0 | 1560 |
肝 | 0.94 | 1550 |
脂肪 | 0.63 | 1450 |
表1
同样在医疗超声系统中通常使用固定的组织声衰减系数扫描特定的人体部位,这会导致图像的亮度不均匀。为了解决这一问题,超声仪器生产厂商往往通过时间增益补偿(TimeGain Control,TGC)滑条,让超声操作者手动给予探头接收到的回波信号一定数值的能量补偿。虽然此种交互式的增益补偿方法操作简单,且易于实现,但此方法需要操作员的参与及主观判断,不仅费时,而且必须随着深度、频率、扫描部位等成像参数的改变而再次调整。
因此,使用一种自适应的医疗超声系统成像参数的优化方法,对声速,组织声衰减系数等声学参数进行优化,具有重要的临床意义。
发明内容
本发明的目的在于克服上述技术缺陷,提供一种能提高图像质量并且省时的自适应的医疗超声系统成像中的声学参数优化方法,该方法首先在声衰减参数优化器中通过多次迭代改变组织声衰减系数得到一系列的二维灰度图像,然后使用纹理分析的方法计算出最优的组织声衰减系数,并设其为当前扫描使用的声衰减系数,从而得到亮度均匀的超声影像,为下一步声速优化打下基础。在此基础上在声速优化器中通过多次迭代改变发射声速,获得一系列的二维灰度图像。最后通过二维图像分析找到最佳的声速,并把该声速设置为当前扫描使用的系统声速。
其技术方案为:
一种自适应的医疗超声系统成像中的声学参数优化方法,包括如下步骤:
1)在声衰减优化器中使用不同的组织声衰减系数补偿同一幅二维灰度图像,声衰减系数的迭代范围为0.3dB/cmMHz~1.3dB/cmMHz,从而获得多幅使用不同声衰减系数得到的二维灰度图像;
2)以小窗口遍历每幅二维灰度图像,根据如下公式计算每个小窗口的纹理能量值,根据各个小窗口的纹理能量值计算每张灰度图像的平均纹理能量值;
能量:∑iPs(i)2∑jPd(j)2
3)在所有二维灰度图像中,找到能量值最大的那幅灰度图对应的组织声衰减系数;
4)把得到的组织声衰减系数设置为当前扫描的组织声衰减系数;
5)在声速优化器中使用不同声速对同一扫描部位发射超声波,声速的迭代范围为1450m/s~1640m/s,从而获得多幅不同声速对应的二维灰度图像;
6)图像分析:
a)以小窗口遍历每幅二维灰度图像,根据如下公式计算每个小窗口的纹理熵值,并根据各个小窗口的纹理熵值计算每幅二维灰度图像的平均纹理熵值熵:-∑iPs(i)logPs-∑jPd(j)logPd
b)选出熵值最大的Z(Z>5)幅二维图像;
c)根据如下公式计算选出的Z幅二维灰度的侧向对比度;
此处,N为从图像中沿探头侧向方向选取的长度为M的整条采样线的数目,di、μi分别代表第i条采样线所有像素点的逐点灰度差异绝对值总和与平均灰度值;
7)找出Z幅二维灰度图像的侧向对比度最大的那幅灰度图像对应的声速;
8)把得到声速设置为当前扫描的使用的系统声速。
与现有技术相比,本发明的有益效果为:
本发明一种自适应的医疗超声系统成像中的声学参数优化方法,操作简单,省时,对声速,组织声衰减系数等声学参数进行优化,进而提高超声系统成像质量,获得亮度更均匀,分辨率更好,对比度更强的超声影像,具有重要临床意义。
附图说明
图1是超声B模式成像系统框图;
图2是本发明方法步骤1)中组织声衰减优化器的的示意图;
图3是本发明方法步骤1)的示意图;
图4是本发明方法步骤2)和步骤6)的示意图;
图5是本发明方法步骤5)中声速优化器的示意图;
图6是本发明方法步骤5)的示意图;
图7是在使用本发明方法前的超声系统图像;
图8是使用本发明方法方法后的超声图像。
具体实施方式
下面结合附图与具体实施方式对本发明作进一步详细地说明。
一种自适应的医疗超声系统成像中的声学参数优化方法,包括如下步骤:
1)图1为传统的B模式超声成像标准的系统框图,它包括波束形成器、回波信号处理器,扫描转换器、图像处理器、显示这五大模块。在此基础上,我们添加一个组织声衰减系 数优化器,此模块只在进行声衰减系数优化时打开,参照图2。在图2的201这个子模块中,使用不同的组织声衰减系数补偿同一幅二维灰度图像,声衰减系数的迭代范围为0.3dB/cmMHz~1.3dB/cmMHz,每次迭代0.05dB/cmMHz。从而在进入扫描转化器前获取多幅使用不同声衰减系数得到的二维灰度图像,并保存图像和当前图像对应的组织声衰减系数。整个迭代过程如图3所示;
2)在图2的202中,使用小窗口遍历每幅二维灰度图像,并根据如下公式计算每个小窗口的纹理能量值,平均各个小窗口的纹理能量值,从而每张灰度图像的平均纹理能量值。小窗口遍历一幅二维灰度图像的过程如图4所示,小窗口的大小为3*3,9*9或16*16等等。
能量:∑iPs(i)2∑jPd(j)2
3)在图2的203中,在所有二维灰度图像中找到“能量值”最大的那幅灰度图对应的组织声衰减系数;
4)把得到的组织声衰减系数设置为当前扫描的组织声衰减系数;
5)图5为声速优化器的系统框图,首先在声速优化器的301中,使用不同声速对同一扫描部位发射超声波,声速的迭代范围为1450m/s~1640m/s,每次迭代5m/s。从而获得多幅不同声速对应的二维灰度图像,保存灰度图像和该图像对应的声速参数,整个迭代发射的过程如图6所示;
6)图像分析:
a)在图5的302中,使用小窗口遍历每幅二维灰度图像,如图4所示。根据如下公式计算每个小窗口的纹理熵值,并根据各个小窗口的纹理熵值计算每幅二维灰度图像的平均纹理熵值
熵:-ΣPs(i)logPs-ΣjPd(j)logPd
b)选出熵值最大的Z(Z>5)幅二维图像;
c)在图5的303中,根据如下公式计算选出的Z幅二维灰度的侧向对比度;
此处,N为从图像中沿探头侧向方向选取的长度为M的整条采样线的数目,di、μi分别代表第i条采样线所有像素点的逐点灰度差异绝对值总和与平均灰度值;
7)在图5的304中,找出Z幅二维灰度图像的侧向对比度最大的那幅灰度图像对应的声速;
8)把得到声速设置为当前扫描的使用的系统声速。
通过图7和图8的对比可以明显看出使用了本发明方法后图像的对比度明显提高,分辨率更好,也更细腻,聚焦效果也更好。
以上所述,仅为本发明较佳的具体实施方式,本发明的保护范围不限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可显而易见地得到的技术方案的简单变化或等效替换均落入本发明的保护范围内。
Claims (1)
1.一种自适应的医疗超声系统成像中的声学参数优化方法,其特征在于,包括如下步骤:
1)在声衰减优化器中使用不同的组织声衰减系数补偿同一幅二维灰度图像,声衰减系数的迭代范围为0.3dB/cmMHz~1.3dB/cmMHz,从而获得多幅使用不同声衰减系数得到的二维灰度图像;
2)以小窗口遍历每幅二维灰度图像,根据如下公式计算每个小窗口的纹理能量值,根据各个小窗口的纹理能量值计算每张灰度图像的平均纹理能量值;
能量:∑iPs(i)2∑jPd(j)2
3)在所有二维灰度图像中,找到能量值最大的那幅灰度图对应的组织声衰减系数;
4)把得到的组织声衰减系数设置为当前扫描的组织声衰减系数;
5)在声速优化器中使用不同声速对同一扫描部位发射超声波,声速的迭代范围为1450m/s~1640m/s,从而获得多幅不同声速对应的二维灰度图像;
6)图像分析:
a)以小窗口遍历每幅二维灰度图像,根据如下公式计算每个小窗口的纹理熵值,并根据各个小窗口的纹理熵值计算每幅二维灰度图像的平均纹理熵值
熵:-∑iPs(i)logPs-∑jPd(j)logPd
b)选出熵值最大的Z(>5)幅二维图像;
c)根据如下公式计算选出的Z幅二维灰度的侧向对比度;
此处,N为从图像中沿探头侧向方向选取的长度为M的整条采样线的数目,di、μi分别代表第i条采样线所有像素点的逐点灰度差异绝对值总和与平均灰度值;
7)找出Z幅二维灰度图像的侧向对比度最大的那幅灰度图像对应的声速;
8)把得到声速设置为当前扫描的使用的系统声速。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210124907XA CN102697521A (zh) | 2012-06-01 | 2012-06-01 | 一种自适应的医疗超声系统成像中的声学参数优化的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210124907XA CN102697521A (zh) | 2012-06-01 | 2012-06-01 | 一种自适应的医疗超声系统成像中的声学参数优化的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102697521A true CN102697521A (zh) | 2012-10-03 |
Family
ID=46890792
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210124907XA Pending CN102697521A (zh) | 2012-06-01 | 2012-06-01 | 一种自适应的医疗超声系统成像中的声学参数优化的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102697521A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103202713A (zh) * | 2013-01-31 | 2013-07-17 | 声泰特(成都)科技有限公司 | 一种医用超声基波和谐波融合的图像优化方法 |
CN104188684A (zh) * | 2014-09-15 | 2014-12-10 | 声泰特(成都)科技有限公司 | 一种自适应医疗超声成像声速优化和信号修正方法及系统 |
CN104463846A (zh) * | 2014-11-04 | 2015-03-25 | 浙江捷尚视觉科技股份有限公司 | 一种用于数字图像处理中的参数调整方法 |
US20160157831A1 (en) * | 2014-12-03 | 2016-06-09 | Samsung Electronics Co., Ltd. | Apparatus and method for supporting computer-aided diagnosis |
CN107997784A (zh) * | 2017-11-29 | 2018-05-08 | 声泰特(成都)科技有限公司 | 一种基于声速自适应修正的超声波束合成方法和系统 |
CN108135579A (zh) * | 2015-10-22 | 2018-06-08 | 株式会社日立制作所 | 超声波诊断装置以及衰减特性测量方法 |
CN110313940A (zh) * | 2019-08-01 | 2019-10-11 | 无锡海斯凯尔医学技术有限公司 | 信号衰减计算方法、装置、设备及计算机可读存储介质 |
CN112137650A (zh) * | 2019-06-28 | 2020-12-29 | 美国西门子医疗系统股份有限公司 | 具有基于脂肪分数优化的声速的超声医学成像 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101273903A (zh) * | 2007-03-28 | 2008-10-01 | 株式会社东芝 | 超声波成像装置以及超声波速度优化方法 |
JP2011229817A (ja) * | 2010-04-30 | 2011-11-17 | Hitachi Aloka Medical Ltd | 超音波診断装置 |
-
2012
- 2012-06-01 CN CN201210124907XA patent/CN102697521A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101273903A (zh) * | 2007-03-28 | 2008-10-01 | 株式会社东芝 | 超声波成像装置以及超声波速度优化方法 |
JP2011229817A (ja) * | 2010-04-30 | 2011-11-17 | Hitachi Aloka Medical Ltd | 超音波診断装置 |
Non-Patent Citations (2)
Title |
---|
李晓瑛,刘东权: "基于图像对比度分析的自适应超声声速优化", 《生物医学工程学杂志》, vol. 28, no. 6, 31 December 2011 (2011-12-31), pages 1092 * |
李晓瑛,刘东权: "基于纹理分析的自适应超声影像衰减补偿", 《四川大学学报(工程科学版)》, vol. 43, no. 3, 31 May 2011 (2011-05-31), pages 141 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103202713A (zh) * | 2013-01-31 | 2013-07-17 | 声泰特(成都)科技有限公司 | 一种医用超声基波和谐波融合的图像优化方法 |
CN104188684A (zh) * | 2014-09-15 | 2014-12-10 | 声泰特(成都)科技有限公司 | 一种自适应医疗超声成像声速优化和信号修正方法及系统 |
CN104188684B (zh) * | 2014-09-15 | 2016-08-31 | 声泰特(成都)科技有限公司 | 一种自适应医疗超声成像声速优化和信号修正方法及系统 |
CN104463846A (zh) * | 2014-11-04 | 2015-03-25 | 浙江捷尚视觉科技股份有限公司 | 一种用于数字图像处理中的参数调整方法 |
CN104463846B (zh) * | 2014-11-04 | 2017-05-17 | 浙江捷尚视觉科技股份有限公司 | 一种用于数字图像处理中的参数调整方法 |
US20160157831A1 (en) * | 2014-12-03 | 2016-06-09 | Samsung Electronics Co., Ltd. | Apparatus and method for supporting computer-aided diagnosis |
CN108135579A (zh) * | 2015-10-22 | 2018-06-08 | 株式会社日立制作所 | 超声波诊断装置以及衰减特性测量方法 |
CN108135579B (zh) * | 2015-10-22 | 2020-08-14 | 株式会社日立制作所 | 超声波诊断装置以及衰减特性测量方法 |
CN107997784A (zh) * | 2017-11-29 | 2018-05-08 | 声泰特(成都)科技有限公司 | 一种基于声速自适应修正的超声波束合成方法和系统 |
CN107997784B (zh) * | 2017-11-29 | 2021-01-12 | 声泰特(成都)科技有限公司 | 一种基于声速自适应修正的超声波束合成方法和系统 |
CN112137650A (zh) * | 2019-06-28 | 2020-12-29 | 美国西门子医疗系统股份有限公司 | 具有基于脂肪分数优化的声速的超声医学成像 |
CN110313940A (zh) * | 2019-08-01 | 2019-10-11 | 无锡海斯凯尔医学技术有限公司 | 信号衰减计算方法、装置、设备及计算机可读存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102697521A (zh) | 一种自适应的医疗超声系统成像中的声学参数优化的方法 | |
US10813595B2 (en) | Fully automated image optimization based on automated organ recognition | |
CN102499711B (zh) | 三维或四维超声图像自动优化调节方法 | |
KR100908252B1 (ko) | 영상 처리 시스템 및 방법 | |
DE102011114943B4 (de) | Bestimmen von Scherwelleninformationen in der medizinischen Ultraschallbildgebung | |
US7833159B2 (en) | Image processing system and method of enhancing the quality of an ultrasound image | |
JP4898809B2 (ja) | 超音波診断装置 | |
US20060094962A1 (en) | Aperture shading estimation techniques for reducing ultrasound multi-line image distortion | |
CN101897600A (zh) | 用于自动超声图像优化的系统和方法 | |
DE102012022010A1 (de) | Adaptive Bildoptimierung in der welleninduzierten Ultraschallbildgebung | |
EP1500371A1 (en) | Ultrasonograph | |
CN106170253B (zh) | 超声波观测装置以及超声波观测装置的工作方法 | |
US9320490B2 (en) | Ultrasound diagnostic apparatus and control method thereof | |
KR102606973B1 (ko) | 초음파 신호 데이터의 깊이 의존 감쇠를 보상하는 방법 및 시스템 | |
CN102695456B (zh) | 波束形成方法、超声波诊断装置、集成电路 | |
CN102753100A (zh) | 超声波诊断装置以及超声波图像处理装置 | |
CN105982694A (zh) | 抑制超声噪声的信号处理方法 | |
CN110840484A (zh) | 自适应匹配最优声速的超声成像方法、装置及超声设备 | |
CN107997784B (zh) | 一种基于声速自适应修正的超声波束合成方法和系统 | |
DE102020116735A1 (de) | Medizinische ultraschallbildgebung mit einer optimierten schallgeschwindigkeit auf der grundlage des fettanteils | |
DE112020001601T5 (de) | Asymmetrie für akustischen Strahlungskraftimpuls | |
CN114847983B (zh) | 一种基于预扫描的乳腺超声ct反射成像参数自适应补偿方法 | |
US20230061869A1 (en) | System and methods for beamforming sound speed selection | |
CN104188684B (zh) | 一种自适应医疗超声成像声速优化和信号修正方法及系统 | |
US8795179B2 (en) | Methods, modules, and systems for gain control in B-mode ultrasonic imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20121003 |