CN105496459B - 用于超声成像设备的自动调节方法及系统 - Google Patents

用于超声成像设备的自动调节方法及系统 Download PDF

Info

Publication number
CN105496459B
CN105496459B CN201610027605.9A CN201610027605A CN105496459B CN 105496459 B CN105496459 B CN 105496459B CN 201610027605 A CN201610027605 A CN 201610027605A CN 105496459 B CN105496459 B CN 105496459B
Authority
CN
China
Prior art keywords
sampling area
energy
index
depth
data
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.)
Active
Application number
CN201610027605.9A
Other languages
English (en)
Other versions
CN105496459A (zh
Inventor
吴方刚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
FEIYINUO TECHNOLOGY (SUZHOU) CO LTD
Feiyinuo Technology Co ltd
Original Assignee
Vinno Technology Suzhou Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Vinno Technology Suzhou Co Ltd filed Critical Vinno Technology Suzhou Co Ltd
Priority to CN201610027605.9A priority Critical patent/CN105496459B/zh
Publication of CN105496459A publication Critical patent/CN105496459A/zh
Application granted granted Critical
Publication of CN105496459B publication Critical patent/CN105496459B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00—Medical imaging apparatus involving image processing or analysis

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Radiology & Medical Imaging (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明提供的超声成像设备的自动调节方法及系统，所述方法包括：在当前扫查深度下，自动遍历每个扫查频率段对应的穿透力指数、分辨力指数、以及对比度指数；根据获取的穿透力指数、分辨力指数、以及对比度指数获取当前扫查频率段对应的图像质量指数；获取所述图像质量指数最大的扫查频率段，将其作为最优的扫查频率段；将超声成像设备调节到当前最优的扫查频率段下进行扫描。本发明的超声成像设备的自动调节方法及系统，通过每个扫查频率段对应的图像质量，自动调节当前的扫查频率段，减少工作流程，提高了超声成像设备临床诊断的方便性和使用效率，提升了超声诊断图像的质量。

Description

用于超声成像设备的自动调节方法及系统
技术领域
本发明属于超声诊断成像领域，涉及一种用于超声成像设备的自动调节方法及系统。
背景技术
医学超声成像以其使用安全、成像速度快、价格便宜、使用方便的优势在临床诊断中大量被采用，是临床诊断的重要工具之一；超声成像可以用于结构形态的观测，血流测量，也可以用于功能成像和组织定征。
超声诊断设备中的一种宽带换能器，其发展使得目前的医学成像系统支持变频成像，即超声成像系统配合探头在其带宽范围内能支持从低频到高频范围的成像；相应的，超声成像过程中，较高的扫查频率能获取高的图像分辨力，但是会牺牲图像的穿透力；反之，较低的扫查频率能获取高的图像穿透力，但是会损失图像的分辨力。
通常情况下，对于同一个超声成像系统，会设置若干扫查频率供医生选择和调节，这些频率基本覆盖探头的带宽范围，相应的，对于每个扫查频率下的图像均经过预先优化。
当医生在这种默认的扫查频率下使用设备时，若发现图像不满足需求，可凭借感官及经验值，手动反复调节扫查频率，以获取最优图像；如此，增加医生的工作量，且操作过程较为繁琐；特别是对于缺乏超声扫查经验的操作者，上述操作也会增加操作和诊断难度，大大降低工作效率。特别的，对于部分操作者，当其在系统的默认设置下通过扫查用户获取图像时，若发现图像未达到其要求，一般只会抱怨图像质量不好，通常不会主动调节设备，如此，导致诊断结果错误。
发明内容
本发明的目的在于提供一种超声成像设备的自动调节方法及系统。
为了实现上述发明目的之一，本发明一实施方式的超声成像设备的自动调节方法，所述方法包括：在当前扫查深度下，自动遍历每个扫查频率段对应的穿透力指数、分辨力指数、以及对比度指数；
根据获取的穿透力指数、分辨力指数、以及对比度指数获取当前扫查频率段对应的图像质量指数；
所述图像质量指数QIQ＝KpQp+KcQGQR，
Kp+Kc＝1，
其中，所述Qp表示穿透力指数，QR表示分辨力指数，QG表示对比度指数，Kp表示穿透力指数的加权系数，Kc表示分辨力指数和对比度指数的加权系数，0.5<Kp<1；
获取所述图像质量指数最大的扫查频率段，将其作为最优的扫查频率段；
将超声成像设备调节到当前最优的扫查频率段下进行扫描。
作为本发明一实施方式的进一步改进，所述穿透力指数的获取方法包括：
在当前扫查深度下，将所述图像等距划分为多个采样区域；
获取每个采样区域内的能量值；
根据各个采样区域的能量值大小变化获取能量分区变化突变点，并获取该能量分区变化突变点所在采样区域的穿透深度；
根据当前的扫查深度以及所述能量分区变化突变点的穿透深度获取所述穿透力指数。
作为本发明一实施方式的进一步改进，所述穿透力指数的获取方法具体包括：
分别对每个采样区域内的图像数据做二维傅里叶变换，获取每个采样区域对应的频域数据；所述图像数据包括：包络数据或图像数据；
每个采样区域对应的频域数据以Wk(u,v)表示，(u,v)表示当前采样区域的频域坐标；k表示当前采样区域，k＝1，2，…，M-1，M，u＝0,1,2,…,U-1，v＝0,1,2,…,V-1；M表示当前扫查深度上，采样区域的数量；
将每个采样区域内的原点变换到频率坐标的(U/2，V/2)处；
对每个采样区域中截止频率半径DT内的频谱能量求和获取每个采样区域对应的能量值；
从最后一个采样区域开始，依次向第一采样区域搜索，将每一个采样区域的能量值分别与预设能量阈值进行比较，若当前采样区域的能量值大于所述预设能量阈值，则将当前能量值归属的采样区域设定为能量分区变化突变点。
作为本发明一实施方式的进一步改进，所述分辨力指数的获取方法包括：
获取能量分区变化突变点分属的采样区域，将该区域以能量突变采样区域表示；
获取当前扫查深度下，第一个采样区域到所述能量突变采样区域分别对应的频域数据；
并根据第一个采样区域到所述能量突变采样区域分别对应的频域数据获取所述分辨力指数。
作为本发明一实施方式的进一步改进，所述对比度指数的获取方法包括：
获取当前扫查深度下，从第一个采样区域到所述能量突变采样区域之间的所述图像数据的直方图；
根据所述直方图获取所述对比度指数。
为了实现上述发明目的之一，本发明一实施方式的超声成像设备的自动调节系统，所述系统包括：
数据获取模块，用于在当前扫查深度下，自动遍历每个扫查频率段对应的穿透力指数、分辨力指数、以及对比度指数；
数据处理模块，用于根据获取的穿透力指数、分辨力指数、以及对比度指数获取当前扫查频率段对应的图像质量指数；
所述图像质量指数QIQ＝KpQp+KcQGQR，
Kp+Kc＝1，
其中，所述Qp表示穿透力指数，QR表示分辨力指数，QG表示对比度指数，Kp表示穿透力指数的加权系数，Kc表示分辨力指数和对比度指数的加权系数，0.5<Kp<1；
调节模块，用于获取所述图像质量指数最大的扫查频率段，将其作为最优的扫查频率段；
将超声成像设备调节到当前最优的扫查频率段下进行扫描。
作为本发明一实施方式的进一步改进，所述数据获取模块具体用于：
在当前扫查深度下，将所述图像等距划分为多个采样区域；
获取每个采样区域内的能量值；
根据各个采样区域的能量值大小变化获取能量分区变化突变点，并获取该能量分区变化突变点所在采样区域的穿透深度；
根据当前的扫查深度以及所述能量分区变化突变点的穿透深度获取所述穿透力指数。
作为本发明一实施方式的进一步改进，所述数据获取模块具体用于：
分别对每个采样区域内的图像数据做二维傅里叶变换，获取每个采样区域对应的频域数据；所述图像数据包括：包络数据或图像数据；
每个采样区域对应的频域数据以Wk(u,v)表示，(u,v)表示当前采样区域的频域坐标；k表示当前采样区域，k＝1，2，…，M-1，M，u＝0,1,2,…,U-1，v＝0,1,2,…,V-1；M表示当前扫查深度上，采样区域的数量；
将每个采样区域内的原点变换到频率坐标的(U/2，V/2)处；
对每个采样区域中截止频率半径DT内的频谱能量求和获取每个采样区域对应的能量值；
从最后一个采样区域开始，依次向第一采样区域搜索，将每一个采样区域的能量值分别与预设能量阈值进行比较，若当前采样区域的能量值大于所述预设能量阈值，则将当前能量值归属的采样区域设定为能量分区变化突变点。
作为本发明一实施方式的进一步改进，所述数据获取模块具体用于：
获取能量分区变化突变点分属的采样区域，将该区域以能量突变采样区域表示；
获取当前扫查深度下，第一个采样区域到所述能量突变采样区域分别对应的频域数据；
根据第一个采样区域到所述能量突变采样区域分别对应的频域数据获取所述分辨力指数。
作为本发明一实施方式的进一步改进，所述数据获取模块具体用于：
获取当前扫查深度下，从第一个采样区域到所述能量突变采样区域之间的所述图像数据的直方图；
根据所述直方图获取所述对比度指数。
与现有技术相比，本发明的有益效果是：本发明的超声成像设备的自动调节方法及系统，通过每个扫查频率段对应的图像质量，自动调节当前的扫查频率段，减少工作流程，提高了超声成像设备临床诊断的方便性和使用效率，提升了超声诊断图像的质量。
附图说明
图1为本发明一实施方式提供的用于超声成像设备的自动调节方法的流程示意图；
图2是本发明一实施方式提供的用于超声成像设备的自动调节系统的模块示意图；
图3是本发明一具体实施方式中，采样区间划分示意图；
图4是本发明一具体实施方式中，其中一个采样区间内像素数据的频域表达图；
图5A、5B是本发明一具体实施方式中，图像的能量曲线图；
图6是本发明一具体示例中，其中两个频率段下的图像直方图示意图。
具体实施方式
以下将结合附图所示的实施方式对本发明进行详细描述。但实施方式并不限制本发明，本领域的普通技术人员根据这些实施方式所做出的结构、方法、或功能上的变换均包含在本发明的保护范围内。
需要说明的是，为了方便描述，将当前扫查深度下，扫查区域中靠近探头表面最近的采样区域定义为第一个采样区域，将远离探头表面最远的采样区域定义为最后一个采样区域。
如图1所示，图1本发明一实施方式提供的用于超声成像设备的自动调节方法，所述方法包括：
S1、在当前扫查深度下，自动遍历每个扫查频率段对应的穿透力指数、分辨力指数、以及对比度指数；
通常情况下，对应一个图像深度，超声探头会设置若干个扫查频率段，每个扫查频率段对应的图像均做事先优化，该优化的图像通常适应大多数的用户，但对个别用户，还需要重新调整影响图像质量的各种参数。
本发明具体实施方式中，对应每个扫查频率段，均获取其对应的穿透力指数、分辨力指数、以及对比度指数，以供使用。
本发明优选实施方式中，所述穿透力指数的获取方法包括：
P1、在当前扫查深度下，将所述图像等距划分为多个采样区域；
结合图3所示，假设当前的扫查深度为D，并在当前扫查深度D下，将所述图像等距划分为M个采样区域，所述扫查区域从靠近所述超声探头表面开始，分别为采样区域1、采样区域2，…，采样区域M-1，采样区域M；
进一步的，每一个采样区域均对应一个穿透深度，从靠近所述超声探头表面开始，每个采样区域依次对应的穿透深度分别为：Z1，Z2，…，ZM-1，ZM。
进一步的，所述穿透力指数的获取方法还包括：
P2、获取每个采样区域内的能量值。
P3、根据各个采样区域的能量值大小变化获取能量分区变化突变点，并获取该能量分区变化突变点的穿透深度。
获取每个采样区域的能量值的方法有多种，本发明一优选实施方式，所述步骤P2具体包括：
M1、分别对每个采样区域内的图像数据做二维傅里叶变换，获取每个采样区域对应的频域数据。
在对图像处理过程中，如果所选择的基础图像质量不高，则在进一步的对图像进行后处理时，也不能获取好的图像质量，因此，选择对基础图像的评价较为客观。
在本发明的具体实施方式中，作为基础图像的所述图像数据包括：包络数据或图像数据，以该两种图像作为基础图像，可以较为客观和准确的评价图像质量。
本发明具体实施方式中，将原始的图像进行信号中处理抽取之后，且未做扫描变化之前的图像定义为包络数据，如此，当对包络数据进行处理过程中，可以大大减少分析的运算量，提高扫查频率自动控制的速度和效率。
将正常的信号链路上关闭除扫描变换之外的所有后处理而获取的图像称之为图像数据，采用图像数据作为获得各种参数的基础，同样可以提高扫查频率自动控制的速度和效率。
本发明一具体示例中，将图像数据以Ak(l,c)表示，
则，
其中，Wk(u,v)表示当前采样区域对应的频域数据，(u,v)表示当前采样区域的频域坐标；k表示当前采样区域，k＝1，2，…，M-1，M，u＝0,1,2,…,U-1，
v＝0,1,2,…,V-1；M表示当前扫查深度上，采样区域的数量。进一步的，所述方法还包括：
M2、将每个采样区域内的原点变换到频率坐标的(U/2，V/2)处；对每个采样区域中截止频率半径DT内的频谱能量求和获取每个采样区域对应的能量值。
结合图4所示，若D(u，v)≤DT，则能量加权系数为1，若D(u，v)＞DT，则能量加权系数为0。
以公式表示如下：
其中，H(u，v)表示能量加权系数，D(u,v)表示像素点(u,v)距离中心点的距离；DT为系统预设的阈值，其可以根据需要自行设定，在此不做详细赘述。
进一步的，每个采样区域对应的能量值以公式表示如下：
其中，P(k)表示每个采样区域对应的能量值，H(u，v)表示能量加权系数。
进一步的，所述P3具体包括：
从最后一个采样区域开始，依次向第一采样区域搜索，将每一个采样区域的能量值分别与预设能量阈值进行比较，若当前采样区域的能量值大于所述预设能量阈值，则将当前能量值归属的采样区域设定为能量分区变化突变点，并获取该能量分区变化突变点的穿透深度。
所述预设能量阈值可根据需要自行设定，也可以根据能量值的分布区间获取。
例如：采用统计算法获取预设能量阈值；或者将前面n个采样区间和后面n个采样区间的能量值的均值作为预设能量阈值，所述n为正整数，且n＜M/2。
通常情况下，信号穿透不到的地方表现为噪声，噪声水平要比组织信号弱，故，通过能量的大小变化可以确定能量分区变化突变点。
本发明一优选实施方式中，结合图5A、5B所示，根据各个采样区间的能量值绘制能量曲线图，从最后一个采样区域开始，依次向第一采样区域搜索，查找系统允许波动范围外的，能量值由小到大的突变点，将该突变点确定为能量分区变化突变点。
本发明一具体实施方式中，能量分区变化突变点的获取具有两种情况：
其中一种情况如图5A所示，该能量分区变化突变点所在的采样区间涵盖了部分组织区域和部分噪声区域，此时，能量分区变化突变点表现并不明显，该种情况下，可以将A点确定为能量分区变化突变点，也可以将B点确定为能量分区变化突变点。
其中另一中情况如图5B所示，该能量分区变化突变点所在的采样区间或者涵盖组织区域，或者涵盖噪声区域，此时，能量分区变化突变点表现十分明显，将C点确定为能量分区变化突变点。
进一步的，获取所述能量分区变化突变点的穿透深度。
为了方便描述，将所述能量分区变化突变点的穿透深度以Zbp表示，Zbp的取值范围为(Z1，Z2，…，ZM-1，ZM)。
进一步的，所述穿透力指数的获取方法还包括：
P4、根据当前的扫查深度以及所述能量分区变化突变点的穿透深度获取所述穿透力指数。
以公式表示如下：
其中，Qp表示穿透力指数，Zbp表示能量分区变化突变点的穿透深度，D表示扫查深度。
本发明中，QP反应了图像的穿透力水平，QP越大，图像的穿透力越强。
进一步的，本发明的分辨力指数的获取方法具体包括：
获取能量分区变化突变点分属的采样区域，将该区域以能量突变采样区域表示；获取当前扫查深度下，第一个采样区域到所述能量突变采样区域分别对应的频域数据；
根据第一个采样区域到所述能量突变采样区域分别对应的频域数据获取所述分辨力指数。
本发明一具体示例中，将能量突变采样区域以采样区域B表示，1<＝B<＝M；
其中，QR表示分辨力指数，x的取值范围为：1、2、3…B-1，B，Wx(u,v)表示第一个采样区域到所述能量突变采样区域中任一个采样区域对应的频域数据；H(u，v)表示分辨力加权系数，其获取方式与能量加权系数相同，在此不做详细赘述。
本发明中，QR反应的是图像穿透力范围内的高频分量的多少,高频分量越多，图像的分辨率越高；反之，图像的分辨力较差；亦即分辨力指数QR越高，图像的分辨率越高，反之越弱。
进一步的，本发明的所述对比度指数的获取方法具体包括：
获取当前扫查深度下，从第一个采样区域到所述能量突变采样区域之间的所述图像数据的直方图；根据所述直方图获取所述对比度指数。
结合图6所示，本发明一具体示例中，将每个频率段下的第一个采样区域到所述能量突变采样区域之间的所述图像数据的直方图以H(g)表示；
其中，QG表示对比度指数，g＝0,1,2，…G，G为图像数据对应的灰阶数，为直方图上反映出来的图像的平均亮度。
本发明中，QG反映直方图的分布差异，也能反映图像的对比度，QG值越小，分布较集中的直方图对应的图像对比度较差，QG值越大，图像的对比度越强。
如图6所示，直方图A对应的图像其对比度要比直方图B对应的图像的对比度要高。
进一步的，本发明一实施方式中，用于超声成像设备的自动调节方法还包括：
S2、根据获取的穿透力指数、分辨力指数、以及对比度指数获取当前扫查频率段对应的图像质量指数；
所述图像质量指数
QIQ＝KpQp+KcQGQR，
Kp+Kc＝1，
其中，所述QIQ表示图像质量指数，所述Qp表示穿透力指数，QR表示分辨力指数，QG表示对比度指数，Kp表示穿透力指数的加权系数，Kc表示分辨力指数和对比度指数的加权系数，0.5<Kp<1；
本发明一优选实施方式中，Kp优选取值参数为0.618。
进一步的，本发明一实施方式中，用于超声成像设备的自动调节方法还包括：
S3、获取所述图像质量指数最大的扫查频率段，将其作为最优的扫查频率段；将超声成像设备调节到当前最优的扫查频率段下进行扫描。
遍历每个扫查频率段对应的图像质量指数，获取图像质量指数最大时对应的扫查频率段，并将超声成像设备调节到该扫查频率段下进行扫描。
可以理解的是，超声成像设备的自动调节方法可以在设备刚启动时自动运行，也可以在超声成像设备上控制一个自动控制按钮，当通过外接设备开启自动控制按钮后，超声成像设备的自动调节方法即开始运行，如此，在运算成本消耗不高的情况下实现自动频率控制，在此不做详细赘述。
结合图2所示，本发明一实施方式中，用于超声成像设备的自动调节系统包括：数据获取模块100、数据处理模块200、调节模块300。
数据获取模块100用于在当前扫查深度下，自动遍历每个扫查频率段对应的穿透力指数、分辨力指数、以及对比度指数；
通常情况下，对应一个图像深度，超声探头会设置若干个扫查频率段，每个扫查频率段对应的图像均做事先优化，该优化的图像通常适应大多数的用户，但对个别用户，还需要重新调整影响图像质量的各种参数。
本发明具体实施方式中，对应每个扫查频率段，数据获取模块100均获取其对应的穿透力指数、分辨力指数、以及对比度指数，以供使用。
本发明优选实施方式中，获取所述穿透力指数时，数据获取模块100具体用于：在当前扫查深度下，将所述图像等距划分为多个采样区域；
结合图3所示，假设当前的扫查深度为D，并在当前扫查深度D下，数据获取模块100将所述图像等距划分为M个采样区域，所述扫查区域从靠近所述超声探头表面开始，分别为采样区域1、采样区域2，…，采样区域M-1，采样区域M；
进一步的，每一个采样区域均对应一个穿透深度，从靠近所述超声探头表面开始，每个采样区域依次对应的穿透深度分别为：Z1，Z2，…，ZM-1，ZM。
进一步的，数据获取模块100还用于：
获取每个采样区域内的能量值；根据各个采样区域的能量值大小变化获取能量分区变化突变点，并获取该能量分区变化突变点的穿透深度。
数据获取模块100获取每个采样区域的能量值的方法有多种，本发明一优选实施方式，数据获取模块100具体用于：分别对每个采样区域内的图像数据做二维傅里叶变换，获取每个采样区域对应的频域数据。
在对图像处理过程中，如果所选择的基础图像质量不高，则在进一步的对图像进行后处理时，也不能获取好的图像质量，因此，选择对基础图像的评价较为客观。
在本发明的具体实施方式中，作为基础图像的所述图像数据包括：包络数据或图像数据，以该两种图像作为基础图像，可以较为客观和准确的评价图像质量。
本发明具体实施方式中，数据获取模块100将原始的图像进行信号中处理抽取之后，且未做扫描变化之前的图像定义为包络数据，如此，当对包络数据进行处理过程中，可以大大减少分析的运算量，提高扫查频率自动控制的速度和效率。
将正常的信号链路上关闭除扫描变换之外的所有后处理而获取的图像称之为图像数据，采用图像数据作为获得各种参数的基础，同样可以提高扫查频率自动控制的速度和效率。
本发明一具体示例中，将图像数据以Ak(l,c)表示，
则，
其中，Wk(u,v)表示当前采样区域对应的频域数据，(u,v)表示当前采样区域的频域坐标；k表示当前采样区域，k＝1，2，…，M-1，M，u＝0,1,2,…,U-1，v＝0,1,2,…,V-1；M表示当前扫查深度上，采样区域的数量。
获取所述穿透力指数时，数据获取模块100还用于：将每个采样区域内的原点变换到频率坐标的(U/2，V/2)处；对每个采样区域中截止频率半径DT内的频谱能量求和获取每个采样区域对应的能量值。
结合图4所示，若D(u，v)≤DT，则能量加权系数为1，若D(u，v)＞DT，则能量加权系数为0。
以公式表示如下：
其中，H(u，v)表示能量加权系数，D(u,v)表示像素点(u,v)距离中心点的距离；DT为系统预设的阈值，其可以根据需要自行设定，在此不做详细赘述。
进一步的，每个采样区域对应的能量值以公式表示如下：
其中，P(k)表示每个采样区域对应的能量值，H(u，v)表示能量加权系数。
进一步的，获取所述穿透力指数时，数据获取模块100还用于：
从最后一个采样区域开始，依次向第一采样区域搜索，将每一个采样区域的能量值分别与预设能量阈值进行比较，若当前采样区域的能量值大于所述预设能量阈值，则将当前能量值归属的采样区域设定为能量分区变化突变点，并获取该能量分区变化突变点的穿透深度。
所述预设能量阈值可根据需要自行设定，也可以根据能量值的分布区间获取。
例如：采用统计算法获取预设能量阈值；或者将前面n个采样区间和后面n个采样区间的能量值的均值作为预设能量阈值，所述n为正整数，且n＜M/2。
通常情况下，信号穿透不到的地方表现为噪声，噪声水平要比组织信号弱，故，通过能量的大小变化可以确定能量分区变化突变点。
本发明一优选实施方式中，结合图5A、5B所示，根据各个采样区间的能量值绘制能量曲线图，从最后一个采样区域开始，依次向第一采样区域搜索，查找系统允许波动范围外的，能量值由小到大的突变点，将该突变点确定为能量分区变化突变点。
本发明一具体实施方式中，能量分区变化突变点的获取具有两种情况：
其中一种情况如图5A所示，该能量分区变化突变点所在的采样区间涵盖了部分组织区域和部分噪声区域，此时，能量分区变化突变点表现并不明显，该种情况下，可以将A点确定为能量分区变化突变点，也可以将B点确定为能量分区变化突变点。
其中另一中情况如图5B所示，该能量分区变化突变点所在的采样区间或者涵盖组织区域，或者涵盖噪声区域，此时，能量分区变化突变点表现十分明显，将C点确定为能量分区变化突变点。
进一步的，获取所述穿透力指数时，数据获取模块100还用于：获取所述能量分区变化突变点的穿透深度。
为了方便描述，将所述能量分区变化突变点的穿透深度以Zbp表示，Zbp的取值范围为(Z1，Z2，…，ZM-1，ZM)。
进一步的，获取所述穿透力指数时，数据获取模块100还用于：
根据当前的扫查深度以及所述能量分区变化突变点的穿透深度获取所述穿透力指数。
以公式表示如下：
其中，Qp表示穿透力指数，Zbp表示能量分区变化突变点的穿透深度，D表示扫查深度。
本发明中，QP反应了图像的穿透力水平，QP越大，图像的穿透力越强。
进一步的，获取所述分辨力指数时，数据获取模块100具体用于：
获取能量分区变化突变点分属的采样区域，将该区域以能量突变采样区域表示；获取当前扫查深度下，第一个采样区域到所述能量突变采样区域分别对应的频域数据；
根据第一个采样区域到所述能量突变采样区域分别对应的频域数据获取所述分辨力指数。
本发明一具体示例中，将能量突变采样区域以采样区域B表示，1<＝B<＝M；
其中，QR表示分辨力指数，x的取值范围为：1、2、3…B-1，B，Wx(u,v)表示第一个采样区域到所述能量突变采样区域中任一个采样区域对应的频域数据；H(u，v)表示分辨力加权系数，其获取方式与能量加权系数相同，在此不做详细赘述。
本发明中，QR反应的是图像穿透力范围内的高频分量的多少,高频分量越多，图像的分辨率越高；反之，图像的分辨力较差；亦即分辨力指数QR越高，图像的分辨率越高，反之越弱。
进一步的，获取所述对比度指数时，数据获取模块100具体用于：
获取当前扫查深度下，从第一个采样区域到所述能量突变采样区域之间的所述图像数据的直方图；根据所述直方图获取所述对比度指数。
结合图6所示，本发明一具体示例中，将每个频率段下的第一个采样区域到所述能量突变采样区域之间的所述图像数据的直方图以H(g)表示；
其中，QG表示对比度指数，g＝0,1,2，…G，G为当前图像数据对应的灰阶数，为直方图上反映出来的图像的平均亮度。
如图6所示，直方图A对应的图像其对比度要比直方图B对应的图像的对比度要高。
本发明中，QG反映直方图的分布差异，也能反映图像的对比度，QG值越小，分布较集中的直方图对应的图像对比度较差，QG值越大，图像的对比度越强。
进一步的，本发明一实施方式中，数据处理模块用于：根据获取的穿透力指数、分辨力指数、以及对比度指数获取当前扫查频率段对应的图像质量指数；
所述图像质量指数
QIQ＝KpQp+KcQGQR，
Kp+Kc＝1，
其中，所述QIQ表示图像质量指数，所述Qp表示穿透力指数，QR表示分辨力指数，QG表示对比度指数，Kp表示穿透力指数的加权系数，Kc表示分辨力指数和对比度指数的加权系数，0.5<Kp<1；
本发明一优选实施方式中，Kp优选取值参数为0.618。
进一步的，本发明一实施方式中，调整模块用于：获取所述图像质量指数最大的扫查频率段，将其作为最优的扫查频率段；将超声成像设备调节到当前最优的扫查频率段下进行扫描。
遍历每个扫查频率段对应的图像质量指数，获取图像质量指数最大时对应的扫查频率段，并将超声成像设备调节到该扫查频率段下进行扫描。
综上所述，本发明的超声成像设备的自动调节方法及系统，通过每个扫查频率段对应的图像质量，自动调节当前的扫查频率段，减少工作流程，提高了超声成像设备临床诊断的方便性和使用效率，提升了超声诊断图像的质量。
为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现。
通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以保存在保存介质中，如ROKK/RAKK、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，信息推送服务器，或者网络设备等)执行本申请各个实施方式或者实施方式的某些部分所述的方法。
以上所描述的装置实施方式仅仅是示意性的，其中所述作为分离部件说明的模块可以是或者也可以不是物理上分开的，作为模块显示的部件可以是或者也可以不是物理模块，即可以位于一个地方，或者也可以分布到多个网络模块上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施方式方案的目的。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。
本申请可用于众多通用或专用的计算系统环境或配置中。例如：个人计算机、信息推送服务器计算机、手持设备或便携式设备、平板型设备、多处理模块系统、基于微处理模块的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括保存设备在内的本地和远程计算机保存介质中。
应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清楚起见，本领域技术人员应当将说明书作为一个整体，各实施方式中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。
上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明，它们并非用以限制本发明的保护范围，凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。

Claims (10)

1.一种用于超声成像设备的自动调节方法，其特征在于，所述方法包括：
在当前扫查深度下，自动遍历每个扫查频率段对应的穿透力指数、分辨力指数、以及对比度指数；
根据获取的穿透力指数、分辨力指数、以及对比度指数获取当前扫查频率段对应的图像质量指数；
所述图像质量指数QIQ＝KpQp+KcQGQR，
Kp+Kc＝1，
其中，所述Qp表示穿透力指数，QR表示分辨力指数，QG表示对比度指数，Kp表示穿透力指数的加权系数，Kc表示分辨力指数和对比度指数的加权系数，0.5<Kp<1；
获取所述图像质量指数最大的扫查频率段，将其作为最优的扫查频率段；
将超声成像设备调节到当前最优的扫查频率段下进行扫描。
2.根据权利要求1所述的用于超声成像设备的自动调节方法，其特征在于，所述穿透力指数的获取方法包括：
在当前扫查深度下，将图像等距划分为多个采样区域；
获取每个采样区域内的能量值；
根据各个采样区域的能量值大小变化获取能量分区变化突变点，并获取该能量分区变化突变点所在采样区域的穿透深度；
根据当前的扫查深度以及所述能量分区变化突变点的穿透深度获取所述穿透力指数。
3.根据权利要求2所述的用于超声成像设备的自动调节方法，其特征在于，所述穿透力指数的获取方法具体包括：
分别对每个采样区域内的包络数据或图像数据做二维傅里叶变换，获取每个采样区域对应的频域数据；
每个采样区域对应的频域数据以Wk(u,v)表示，(u,v)表示当前采样区域的频域坐标；k表示当前采样区域，k＝1，2，…，M-1，M，u＝0,1,2,…,U-1，v＝0,1,2,…,V-1；M表示当前扫查深度上，采样区域的数量；
将每个采样区域内的原点变换到频率坐标的(U/2，V/2)处；
对每个采样区域中截止频率半径DT内的频谱能量求和获取每个采样区域对应的能量值；
从最后一个采样区域开始，依次向第一采样区域搜索，将每一个采样区域的能量值分别与预设能量阈值进行比较，若当前采样区域的能量值大于所述预设能量阈值，则将当前能量值归属的采样区域设定为能量分区变化突变点。
4.根据权利要求3所述的用于超声成像设备的自动调节方法，其特征在于，所述分辨力指数的获取方法包括：
获取能量分区变化突变点分属的采样区域，将该区域以能量突变采样区域表示；
获取当前扫查深度下，第一个采样区域到所述能量突变采样区域分别对应的频域数据；
并根据第一个采样区域到所述能量突变采样区域分别对应的频域数据获取所述分辨力指数。
5.根据权利要求3所述的用于超声成像设备的自动调节方法，其特征在于，所述对比度指数的获取方法包括：
获取当前扫查深度下，从第一个采样区域到能量突变采样区域之间的所述图像数据的直方图；
根据所述直方图获取所述对比度指数。
6.一种用于超声成像设备的自动调节系统，其特征在于，所述系统包括：
数据获取模块，用于在当前扫查深度下，自动遍历每个扫查频率段对应的穿透力指数、分辨力指数、以及对比度指数；
数据处理模块，用于根据获取的穿透力指数、分辨力指数、以及对比度指数获取当前扫查频率段对应的图像质量指数；
所述图像质量指数QIQ＝KpQp+KcQGQR，
Kp+Kc＝1，
其中，所述Qp表示穿透力指数，QR表示分辨力指数，QG表示对比度指数，Kp表示穿透力指数的加权系数，Kc表示分辨力指数和对比度指数的加权系数，0.5<Kp<1；
调节模块，用于获取所述图像质量指数最大的扫查频率段，将其作为最优的扫查频率段；
将超声成像设备调节到当前最优的扫查频率段下进行扫描。
7.根据权利要求6所述的用于超声成像设备的自动调节系统，其特征在于，所述数据获取模块具体用于：
在当前扫查深度下，将图像等距划分为多个采样区域；
获取每个采样区域内的能量值；
根据各个采样区域的能量值大小变化获取能量分区变化突变点，并获取该能量分区变化突变点所在采样区域的穿透深度；
根据当前的扫查深度以及所述能量分区变化突变点的穿透深度获取所述穿透力指数。
8.根据权利要求7所述的用于超声成像设备的自动调节系统，其特征在于，所述数据获取模块具体用于：
分别对每个采样区域内的包络数据或图像数据做二维傅里叶变换，获取每个采样区域对应的频域数据；
每个采样区域对应的频域数据以Wk(u,v)表示，(u,v)表示当前采样区域的频域坐标；k表示当前采样区域，k＝1，2，…，M-1，M，u＝0,1,2,…,U-1，v＝0,1,2,…,V-1；M表示当前扫查深度上，采样区域的数量；
将每个采样区域内的原点变换到频率坐标的(U/2，V/2)处；
对每个采样区域中截止频率半径DT内的频谱能量求和获取每个采样区域对应的能量值；
从最后一个采样区域开始，依次向第一采样区域搜索，将每一个采样区域的能量值分别与预设能量阈值进行比较，若当前采样区域的能量值大于所述预设能量阈值，则将当前能量值归属的采样区域设定为能量分区变化突变点。
9.根据权利要求8所述的用于超声成像设备的自动调节系统，其特征在于，所述数据获取模块具体用于：
获取能量分区变化突变点分属的采样区域，将该区域以能量突变采样区域表示；
获取当前扫查深度下，第一个采样区域到所述能量突变采样区域分别对应的频域数据；
根据第一个采样区域到所述能量突变采样区域分别对应的频域数据获取所述分辨力指数。
10.根据权利要求8所述的用于超声成像设备的自动调节系统，其特征在于，所述数据获取模块具体用于：
获取当前扫查深度下，从第一个采样区域到能量突变采样区域之间的所述图像数据的直方图；
根据所述直方图获取所述对比度指数。
CN201610027605.9A 2016-01-15 2016-01-15 用于超声成像设备的自动调节方法及系统 Active CN105496459B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610027605.9A CN105496459B (zh) 2016-01-15 2016-01-15 用于超声成像设备的自动调节方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610027605.9A CN105496459B (zh) 2016-01-15 2016-01-15 用于超声成像设备的自动调节方法及系统

Publications (2)

Publication Number Publication Date
CN105496459A CN105496459A (zh) 2016-04-20
CN105496459B true CN105496459B (zh) 2018-09-21

Family

ID=55705030

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610027605.9A Active CN105496459B (zh) 2016-01-15 2016-01-15 用于超声成像设备的自动调节方法及系统

Country Status (1)

Country Link
CN (1) CN105496459B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109996496A (zh) * 2016-08-16 2019-07-09 戈尔丹斯医疗公司 用于超声审查和成像的系统和方法
CN106815840B (zh) * 2017-01-22 2020-06-05 飞依诺科技(苏州)有限公司 一种肝部扫查图像的处理方法及装置
CN107890355B (zh) * 2017-12-19 2021-08-10 飞依诺科技(苏州)有限公司 一种调节超声波发射功率的方法和装置
CN108303471A (zh) * 2018-02-09 2018-07-20 业成科技（成都）有限公司 超声波成像系统及方法
CN111839588B (zh) * 2020-07-29 2024-06-25 深圳开立生物医疗科技股份有限公司 一种超声图像生成方法、装置、设备及存储介质

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8235905B2 (en) * 2009-05-26 2012-08-07 General Electric Company System and method for automatic ultrasound image optimization
US9668714B2 (en) * 2010-04-14 2017-06-06 Maui Imaging, Inc. Systems and methods for improving ultrasound image quality by applying weighting factors
CN101853504B (zh) * 2010-05-07 2012-04-25 厦门大学 基于视觉特性与结构相似度的图像质量评测方法
CN102354388B (zh) * 2011-09-22 2013-03-20 北京航空航天大学 一种自适应计算图像底层特征重要性权重的方法
WO2014159819A1 (en) * 2013-03-13 2014-10-02 Jinhyoung Park System and methods for producing an image from a rotational intravascular ultrasound device
CN103356241B (zh) * 2013-08-02 2015-08-19 中国十九冶集团有限公司职工医院 二维超声设备成像质量评估系统
CN103426173B (zh) * 2013-08-12 2017-05-10 浪潮电子信息产业股份有限公司 一种立体图像质量的客观评价方法
CN104616294A (zh) * 2015-01-21 2015-05-13 同济大学 一种基于视觉显著性的色调映射图像质量评价方法

Also Published As

Publication number Publication date
CN105496459A (zh) 2016-04-20

Similar Documents

Publication Publication Date Title
CN105496459B (zh) 用于超声成像设备的自动调节方法及系统
JP5433097B2 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP4405017B2 (ja) Ｂモード超音波イメージングにおける自動的な時間及び／又は横方向ゲイン補償のための方法及び装置
US9028414B2 (en) Ultrasonic observation apparatus, operation method of the same, and computer readable recording medium
CN1768709A (zh) 超声波多普勒测量装置及其控制方法
JP5079177B2 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
US20120310087A1 (en) Ultrasonic diagnosis apparatus, operation method of the same, and computer readable recording medium
EP2995259A1 (en) Ultrasound optimization method and ultrasonic medical device therefor
JP5054253B2 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP5066306B2 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP2007152111A (ja) スペクトルドプラ利得を自動調整するための方法及び装置
CN105828726A (zh) 医用诊断装置、医用诊断装置的工作方法以及医用诊断装置的工作程序
WO2012063930A1 (ja) 超音波診断装置、超音波診断装置の作動方法および超音波診断装置の作動プログラム
US20190282205A1 (en) Ultrasound imaging system and ultrasound imaging method
CN106983524A (zh) 一种反映生物组织异常的参数及其测量方法
CN105792756A (zh) 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序
WO2012011414A1 (ja) 超音波診断装置、超音波診断装置の作動方法および超音波診断装置の作動プログラム
US20230119063A1 (en) Methods and Systems for Evaluating Echo Data Contemporaneous with an Electrodiagnostic Study
US10338033B2 (en) Typing and imaging of biological and non-biological materials using quantitative ultrasound
Soylu et al. Calibrating data mismatches in deep learning-based quantitative ultrasound using setting transfer functions
CN105581812B (zh) 用于超声成像设备的自动调节方法及系统
CN106102590A (zh) 超声波诊断装置
CN112566559B (zh) 具有像素外推图像增强的超声成像系统
Wang et al. Effects of instrument settings on radiofrequency ultrasound local estimator images: A preliminary study in a gallbladder model
WO2014084278A1 (ja) 超音波診断装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder

Address after: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province

Patentee after: Feiyinuo Technology Co.,Ltd.

Address before: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province

Patentee before: Feiyinuo Technology (Suzhou) Co.,Ltd.

CP01 Change in the name or title of a patent holder
CP03 Change of name, title or address

Address after: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province

Patentee after: Feiyinuo Technology (Suzhou) Co.,Ltd.

Address before: 215123 unit 501, building C8, bio nano Park, No. 218, Xinghu street, industrial park, Suzhou, Jiangsu Province

Patentee before: VINNO TECHNOLOGY (SUZHOU) Co.,Ltd.

CP03 Change of name, title or address