CN113491535A - 超声波诊断装置、信号处理装置以及信号处理程序 - Google Patents
超声波诊断装置、信号处理装置以及信号处理程序 Download PDFInfo
- Publication number
- CN113491535A CN113491535A CN202110190790.4A CN202110190790A CN113491535A CN 113491535 A CN113491535 A CN 113491535A CN 202110190790 A CN202110190790 A CN 202110190790A CN 113491535 A CN113491535 A CN 113491535A
- Authority
- CN
- China
- Prior art keywords
- signal
- roi
- dimensional
- reception
- representative
- 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.)
- Granted
Links
- 238000012545 processing Methods 0.000 title claims abstract description 35
- 238000000034 method Methods 0.000 claims description 23
- 238000004458 analytical method Methods 0.000 claims description 18
- 239000000523 sample Substances 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000002604 ultrasonography Methods 0.000 claims description 8
- 230000007423 decrease Effects 0.000 claims description 5
- 238000003384 imaging method Methods 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 description 19
- 238000010586 diagram Methods 0.000 description 11
- 238000012937 correction Methods 0.000 description 8
- 210000001519 tissue Anatomy 0.000 description 8
- 238000007689 inspection Methods 0.000 description 5
- 238000003745 diagnosis Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 201000010099 disease Diseases 0.000 description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 208000004930 Fatty Liver Diseases 0.000 description 2
- 206010019708 Hepatic steatosis Diseases 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 208000010706 fatty liver disease Diseases 0.000 description 2
- 238000001727 in vivo Methods 0.000 description 2
- 210000004185 liver Anatomy 0.000 description 2
- 208000019423 liver disease Diseases 0.000 description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 231100000240 steatosis hepatitis Toxicity 0.000 description 2
- 206010012426 Dermal cyst Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 206010033675 panniculitis Diseases 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 210000004304 subcutaneous tissue Anatomy 0.000 description 1
Images
Classifications
-
- 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
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
- A61B8/467—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means
- A61B8/469—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means for selection of a region of interest
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
本发明提供超声波诊断装置、信号处理装置以及信号处理程序。不管方向如何都能精度良好地算出超声波的衰减率。在将接收波束形成后的接收信号按照接收扫描线的顺序进行排列而形成的二维信号空间设定多个二维的ROI。使用ROI内的信号求取代表值,将求得的代表值作为设定于ROI内的代表点的处理后信号值。使用该处理后信号算出衰减率等。
Description
技术领域
本发明涉及超声波诊断装置,涉及被检体内的声衰减特性的测量手法。
背景技术
以超声波、MRI(Magnetic Resonance Imaging,磁共振成像)装置、X射线CT(Computed Tomography,计算机断层扫描)装置为代表的医疗用的图像显示装置作为将不能目视的生物体内的信息以数值或图像的形态进行提示的装置而被广泛利用。尤其地,利用超声波来显示图像的超声波摄像装置与其他装置比较具备更高的时间分辨率,例如具有能将跳动下的心脏没有洇渗地图像化的性能。
超声波诊断装置向检查对象发送超声波,接收来自散射体的反射信号(以下称作接收RF(Radio Frequency,射频)信号),从而构成图像。基本上通过根据收发所需的时间和音速来测量到散射体的距离,构成基于超声波能量的强度的亮度的空间分布,由此生成超声波图像(以下称作B像)。
另外,提出开发超声波诊断装置通过使用接收信号定量测量血流速度、组织的硬度、衰减特性等声的特性来进行疾患的诊断的技术。
例如在专利文献1中公开了求取生物体组织的衰减特性的手法。在专利文献1中,首先从探头对对象物发送使焦距不同的第1以及第2会聚的超声波波束,分别得到关于对象物的给定的测量点的第1以及第2超声波波束所引起的接收信号。接下来,衰减特性算出部使用关于给定的测量点得到的第1超声波波束所引起的接收信号和第2超声波波束所引起的接收信号的超声波能量的深度方向的变化,来求取对象物的衰减特性。
另外,在专利文献2中提出如下方法:该方法是一种为了减低超声波图像中所含的各种噪声、散斑,强调与被检体组织相关的信息,进行图像的滤波处理的技术,在肌肉组织结构那样在局部区域来看的情况下组织结构向特定的方向对齐时,调整滤波器的参数,以使该结构清晰。即,算出超声波图像中的局部区域的边缘的方向、所述局部区域中的边缘的大小、和局部区域中的边缘的方向一致性,设定滤波器的参数。
现有技术文献
专利文献
专利文献1:JP专利第6457107号公报
专利文献2:JP特开2012-75882号公报
作为生物体组织的声特性之一,有在组织内部传播的声能的衰减特性。所谓衰减特性,是超声波在媒介中传播时超声波能量对应于其传播距离而通过吸收、散射的现象进行衰减的特性。肝疾患之一的脂肪肝已知与正常肝比较声能的衰减变大。因此,期待通过测量衰减率而有助于疾患的诊断。
超声波的信号强度随着传播距离变大而平滑地减少,但实际中包含随机的信号强度的变化(散斑)。衰减率的算出由于是根据每单位距离的信号强度的降低率来进行算出,因此该散斑会引起测量值的偏差,有会使衰减率的算出精度降低的问题。
另外,在超声波诊断装置中,安装有对应于深度而在接收信号上乘上增益的时间增益控制(Time Gain Control(TGC))功能,进行调整,使得将超声波能量的强度分布表征为亮度分布的B像的亮度从浅部到深部都一样。但实际的接收RF信号的深度方向的衰减有不与深度相应的情况。例如在将发送超声波会聚而照射的情况下,接收RF信号在发送超声波的焦点附近强度变大。
另外,实际的接收RF信号的衰减不仅在深度方向上,还在深度方向和方位方向上二维地产生。例如在设想肝疾患的诊断的情况下,在肝实质和皮下组织、囊肿、横隔膜中,由于音速、密度、衰减率不同,因此会在结构物的边界引起超声波反射、折射、散射,给深度方向和方位方向的声压轮廓(profile)二维地带来影响。
因此,在仅深度方向的TGC功能中,有如下情况:不能使深度方向的B像的亮度分布一样,依然在B像上难以视觉辨识生物体组织的形态,诊断变得困难。
发明内容
本发明的目的在于,不管方向如何都能精度良好地算出超声波的衰减率。
为了达成上述目的,本发明的超声波诊断装置具有:接收波束形成器,其通过从超声波探头接受信号并进行接收波束形成,来针对设定于摄像区域的多个接收扫描线分别生成接收信号;和信号处理部,其对接收信号进行处理。信号处理部具有:二维ROI设定部,其在将多个接收信号按照接收扫描线的顺序进行排列的二维信号空间设定多个二维的ROI;和代表值算出部,其使用ROI内的信号来求取代表值,将求得的代表值作为设定于ROI内的代表点的处理后信号值。
发明的效果
根据本发明的超声波诊断装置,由于设定多个二维的ROI,将各自的代表值作为处理后信号值,因此不管方向如何都能精度良好地算出超声波的衰减率。
附图说明
图1是本发明的实施方式1的超声波诊断装置的框图。
图2是表示实施方式1的超声波诊断装置的动作的流程图。
图3是在实施方式1示出向检查对象的收发的说明图。
图4是表示设定为将实施方式1的接收RF信号进行排列的二维信号空间的ROI的说明图。
图5是表示实施方式1的超声波诊断装置的动作的流程图。
图6是表示在实施方式1中判断ROI适合性的处理的说明图。
图7是表示实施方式1的超声波诊断装置的动作的流程图。
图8是在实施方式1中示出接收RF信号的二维信号空间中的有效ROI和无效ROI的说明图。
图9是在实施方式1中示出接收RF信号的二维信号空间中的有效区域(有效ROI)和无效区域(无效ROI)的分布的说明图。
图10是表示实施方式1的超声波诊断装置的动作的流程图。
图11是在实施方式1中示出代表值适合性判定部的处理的说明图。
图12是表示实施方式2的超声波诊断装置的动作的流程图。
图13是实施方式2的显示画面例。
图14是表示实施方式3的超声波诊断装置的动作的流程图。
图15是实施方式3的显示画面例。
图16是实施方式3的显示画面例。
附图标记的说明
1:超声波收发装置
10:探头
13:外部输入设备
16:显示部
20:收发控制部
21:发送波束形成器
22:接收波束形成器
21:发送波束形成器
30:信号处理部
31:二维ROI设定部
32:ROI适合性分析部
33:代表值算出部
34:代表值适合性判定部
35:插补部
36:图像生成部
37:衰减率算出部
100:检查对象
具体实施方式
以下按照附图来说明本发明的实施方式。
<<实施方式1>>
<超声波诊断装置的整体结构>
在图1示出实施方式的超声波诊断装置1的一个结构例的框图。本实施方式的超声波收发装置1具备收发控制部20和信号处理部30。另外,在收发控制部20连接探头10。在信号处理部30连接外部输入设备13和显示部16。
收发控制部20具备发送波束形成器21和接收波束形成器22。发送波束形成器21生成发送信号,对构成探头10的排列的振子输出。探头10的多个振子将发送信号变换成超声波,向检查对象100发送。在检查对象100中反射、散射的超声波当中到达探头10的超声波被探头10的各振子变换成电信号。接收波束形成器22接受各振子的电信号,关于设定于摄像区域的接收扫描线上的多个点使其延迟以使焦点对准,并进行加法运算(调相加法运算),由此生成(接收波束形成)接收信号(以下称作接收RF信号)。
信号处理部30具备二维ROI设定部31、ROI适合性分析部32、代表值算出部33、代表值适合性判定部34、插补部35、图像生成部36和衰减率算出部37而构成。
二维ROI设定部31形成将多个接收RF信号按接收扫描线的顺序进行排列的二维信号空间,对二维信号空间设定多个二维的ROI。
代表值算出部33使用ROI内的信号来求取代表值,将求得的代表值作为设定于ROI内的代表点的处理后信号值。
如此地,将二维的ROI在接收RF信号的二维信号空间中设定多个,求取各自的代表值,将其设为ROI的代表值的处理后信号值,因此能在二维空间中减低ROI内的接收信号中所含的散斑。由此,不管方向如何都能精度良好地算出超声波的衰减率。
期望上述二维ROI设定部31将多个ROI关于沿着二维信号空间的接收扫描线的方向以及接收扫描线的排列方向的至少一方,设定在相互错开的位置。
这时,二维ROI设定部31也可以基于接收RF信号适应性地变更设定的超声波的尺寸、形状。
ROI适合性分析部32判定ROI内的信号强度的分布是否适合预先确定的信号强度分布,将适合的ROI设为有效ROI,将不适合的ROI设为无效ROI。
代表值适合性判定部34判定关于代表值算出部33算出的多个ROI的代表值的分布是否适合预先确定的条件。代表值适合性判定部34并不将不适合条件的代表值采用为处理后信号值。例如,代表值适合性判定部34在多个ROI的代表值的信号强度的分布沿着声波传播方向慢慢减少的情况下,判定为适合条件。
插补部35在二维信号空间中使用多个ROI的代表值,通过插补来求取ROI与ROI之间的二维信号空间的信号值。
图像生成部36使用处理后信号来生成B模式像。图像生成部36可以设为使用处理后信号来将B模式像的亮度值最佳化的结构。
衰减率算出部37使用处理后信号来求取深度方向的衰减率。
信号处理部30由具备CPU(Central Processing Unit,中央处理器)、GPU(Graphics Processing Unit,图形处理器)等处理器、和存储器的计算机等构成,通过CPU将存放于存储器的程序读入并执行,来通过软件实现上述的二维ROI设定部31、ROI适合性分析部32、代表值算出部33、代表值适合性判定部34、插补部35、图像生成部36和衰减率算出部37的功能。另外,信号处理部30还能将其一部分或全部通过硬件实现。例如使用ASIC(Application Specific Integrated Circuit,特定用途集成电路)那样的定制IC、FPGA(Field-Programmable Gate Array,现场可编程门阵列)那样的可编程IC来构成信号处理部7,进行电路设计,以使得实现信号处理部7的各部的功能即可。
在图2示出实施方式的超声波诊断装置的一结构例的流程图。按照流程图来说明整体的动作。
首先,从探头10如图3那样收发超声波。具体地,将由发送波束形成器21生成的发送信号输出到探头10,从探头10发送作为超声波的波束的发送超声波301,探头10接收检查对象100内部的散射体303散射来的反射超声波302。接收波束形成器22通过对来自探头10的信号进行调相加法运算,来生成沿着设定的接收扫描线的接收RF信号。接收RF信号是在调相加法运算处理后进一步以常用对数进行压缩而得到的信号。
另外,在接收RF信号中包含被称作散斑的随机的信号的强度变化。这是由于从媒介内的散射体303散射的超声波随机相互干涉而产生的,根据媒介的声特性、媒介内的散射体的形状以及密度而变化。另外,设想检查对象的媒介有衰减。因此,若超声波在该媒介内部传播,则对应于传播距离而超声波能量减少。
在步骤201,对二维ROI设定部31从接收波束形成器22输入沿着接收扫描线的接收RF信号。
二维ROI设定部31通过将多个接收RF信号按照接收扫描线的顺序在横轴(方位)方向(探头的振子的排列的方向)上进行排列,来形成二维信号空间。因此,该二维信号空间的横轴成为方位方向即探头10的振子所进行排列的方向,纵轴成为相当于深度方向的时间。
接下来,在步骤202,使用发送的超声波的声压轮廓来修正接收信号,得到修正后信号。由此,对发送的超声波的声压的强弱所引起的接收RF信号的信号强度的分布进行修正。声压轮廓使用通过均质的模型(phantom)或仿真预先进行测量、算出并保存在超声波诊断装置内的存储器的声压轮廓。作为修正方法,将接收RF信号设为M(i,t),将声压轮廓设为R(i,t)、将修正后的信号设为C(i,t),从而表征为
【数学式1】
另外,i表示探头的振子的编号,t表示时间。
接下来,在步骤203,二维ROI设定部31决定设定的多个ROI的规格以及位置。在此,二维ROI设定部31将多个ROI关于沿着二维信号空间的接收扫描线的方向以及接收扫描线的排列方向的至少一方,设定在相互错开的位置。另外,二维ROI设定部31基于接收RF信号来适应性地变更设定的ROI的尺寸、形状。
使用图4来说明ROI规格、和位置的设定方法的示例。ROI401是如下那样的ROI:在接收RF信号中,深度方向的尺寸是采样数h,方位方向的尺寸是采样数w的尺寸,将其代表点(例如中心)设定为来到(xn,zn)的位置。ROI402是与ROI401尺寸相同但将中心位置设定为在深度方向上以采样数dh而不同、在方位方向上以采样数dw而不同的位置。因而,相对于ROI401的代表点(xn,zn),ROI402的代表点的位置成为(xn+dw,zn+dhs)。
使用图5来说明二维ROI设定部31所进行的ROI401、402的深度方向的尺寸(采样数h)和方位方向的尺寸(采样数w)的决定方法、以及ROI401与ROI402的位置错开量(深度方向的采样数dh、方位方向的采样数dws)的决定方法的一例。
二维ROI设定部31在图5的步骤501,在步骤201,接受将接收RF信号进行排列而形成的二维信号空间。接下来,在步骤502进行二维傅立叶变换,求取接收RF信号的空间频率的分布,二维ROI设定部31基于该空间频率分布来求取散斑的方位方向、深度方向的尺寸的分布。
接下来在步骤503,二维ROI设定部31基于散斑的空间上的尺寸分布来决定ROI的尺寸。作为具体的决定方法,例如使ROI的深度方向的采样数h、方位方向的采样数w以及深度方向的移动采样数dh、方位方向的移动采样数dw为散斑的空间尺寸的平均的任意倍。
接下来,在图2的步骤204,ROI适合性分析部32判断二维ROI设定部31所设定的ROI是否使用在步骤205以后的处理中、适合性。在图6示出适合性的判断的方法的示意图。首先,在二维ROI设定部31所设定的ROI是ROI601、ROI602的情况下,ROI适合性分析部32分别算出ROI内的各点的接收RF信号的信号强度(亮度)分布603、信号强度(亮度)分布604。然后,ROI适合性分析部32将各个信号强度分布与理想数据605的差异定量化。另外,理想数据使用预先用均质模型、仿真等算出并保存在超声波诊断装置的存储器的数据。
在图6中,ROI601的信号强度(亮度)分布603与理想分布605非常相似。另一方面,包含结构物606的ROI602的信号强度(亮度)分布604从理想分布605偏离。ROI适合性分析部32根据这样从理想分布607的偏离度来判定ROI的适合性的有无。
关于由ROI适合性分析部32判断为有适合性的ROI,在接下来的步骤205,代表值算出部33进行代表值的算出。
使用图7的流程来具体说明上述步骤204、205的ROI适合性分析部32所进行的偏离度的算出、和代表值算出部33所进行的代表值的算出的处理流程。
在图7的步骤701,二维ROI设定部31基于图4所示那样前一个设定的ROI的位置来决定ROI的设定位置,并进行设定。
在步骤702,ROI适合性分析部32算出ROI中的信号强度(亮度)直方图Rmeasured。
在步骤703,ROI适合性分析部32将ROI中的信号强度(亮度直方图Rmeasured与理想的亮度直方图Rmodel的偏离度S算出为
【数学式2】
在步骤704,ROI适合性分析部32通过预先设定的阈值与在式(2)算出的偏离度S的比较来判定适合性。即,ROI适合性分析部32在偏离度S比预先设定的阈值小的情况下判定为有适合性。
在步骤705,ROI适合性分析部32对判定为有适合性的ROI,根据ROI中的值算出ROI的代表点(xn,zn)的值μ,作为有效ROI。作为代表点的值,使用ROI中的平均值、中央值、最大值、最小值这样的值。另外,在本实施方式中,将最初的代表点的位置设为ROI的中央的点,但并不限于中央的点,也可以设为重心、任意的边的中点、四角的点的预先确定的位置的点。
另外,ROI适合性分析部32关于判断为没有适合性的ROI,在步骤706将代表点的值设为常数τ,作为无效ROI。另外,将τ决定成0等在接收RF信号上不出现的常数。
在步骤207,ROI适合性分析部32将算出的代表值保存在存储器(未图示)。
最后在步骤208(步骤707),ROI适合性分析部32直到完成全部区域的分析为止都重复上述处理。
图8是表示图7的步骤701到707的处理中途的样子的图。可知判定为适合的有效ROI803和判定为不适合的无效ROI804位于接收RF信号的二维信号空间805中。对各个ROI设定代表点802。
图9是表示从图7的步骤701完成了707的处理时的样子的图。若处理完成,则成为无效ROI的代表点的集合体即无效区域的903和有效ROI的代表点的集合体即有效区域904在接收RF信号805中分布的形式。
接下来,代表值适合性判定部34判定关于代表值算出部33算出的多个ROI的代表值的分布是否适合预先确定的条件。
首先,在图2的步骤208,代表值适合性判定部34进行代表值的选出。图10是代表值的选出的流程图。
首先在步骤1001,代表值适合性判定部34将代表点间用线段连起来,在步骤1002,作成图11所示的被3点包围的区域1102。
接下来,在步骤1003,代表值适合性判定部34判定在该被3点包围的区域1102内是否不含图9所示的无效区域(无效ROI)903。
最后,在步骤1004,代表值适合性判定部34使用3点的代表点的值来算出通过被3点包围的区域中的超声波波束方向的代表值的梯度。
图示步骤1004的代表值适合性判定部34的操作的是图11。在图11中,代表值适合性判定部34设定通过被k1、k2、k3这3点包围的区域1102的内部的波束向量1101,算出该波束向量1101的方向的代表值的梯度。
接下来在步骤1005,若波束向量1101的方向的梯度为负,则代表值适合性判定部34将代表值设为合适而选出。
在此,作为代表值适合性判定部34在波束向量1101的方向的梯度为负时设为合适而选出的理由,设想为随着在检查对象100中传播而超声波衰减,必定沿着超声波的波束方向即声波的传播方向,能量减少,即代表值的梯度成为负。
其中,由于由ROI适合性分析部32本来就判断为亮度分布不合适的无效区域(无效ROI)能量不一定减少,因此在步骤1003预先除去。
在图2、步骤209,对步骤208中选出的代表点间进行二维的插补处理。由此,能使代表点的分布更加平滑,能使图像生成部36中生成的图像更加平滑。另外,由于即使是空开间隔配置多个ROI的情况,也能算出ROI间的信号值,因此能使处理速度提升。另外,插补方法使用二维的线性插补、二维的样条插补、二维的B样条插补方法等公知的方法。
根据以上,能在二维空间生成在深度方向以及方位方向上除去散斑的信号(代表值以及插补处理值)。
最后在图2的步骤210,插补部35进行图像的生成。关于图像的生成方法,在实施方式2、实施方式3叙述。
<<实施方式2>>
图12表示图2的步骤210中的图像生成部36所进行的B像的图像的生成处理的流程图的示例。
首先,在步骤1201,对图像生成部36输入从图2的步骤209输出的插补处理后的信号。
接下来,在步骤1202,图像生成部36将声压轮廓读入。另外,该声压轮廓与图2的步骤202的声压轮廓修正中所使用的声压轮廓相同。
接下来,在步骤1203,图像生成部36算出用于进行B像的描绘区域整体的增益修正的二维修正增益。若将插补处理后的信号设为S(I,t),将声压轮廓设为P(I,t),则二维修正增益G(I,t)表征为
【数学式3】
G(i、t)=S(i、t)·P(i、t)·B …(3)。
另外,B是用户设定的所期望的B像的亮度。
接下来在步骤1204将接收RF信号读入。
接下来在步骤1205,将在步骤1205读入的接收RF信号通过二维修正增益进行修正,作成修正后的信号。若将修正后的信号设为C(i,t),将接收RF信号设为M(i,t),则表征为
【数学式4】
C(i、t)=M(i、t)-G(i、t) …(4)。
在步骤1206,图像生成部36通过使用式(4)中作成的修正后的信号C(i,t)生成B像,能得到用二维修正增益修正过的B像。图像生成部36将生成的图像显示在显示部16。
图13是在实施方式2作成的修正B像的显示例。将增益分布图1302和修正B像1303进行排列显示,在增益分布图1302中示出在图2的步骤208选出的代表点1202和在图2的步骤204分析的无效区域(无效ROI)1201。B像单体自不必说,用户能如图13所示那样通过与增益分布图(map)进行比较来判断是否适当地进行了修正。
<<实施方式3>>
图14表示图2的步骤210中的图像生成部36所进行的衰减率分布图的图像的生成处理的流程图的示例。衰减率是每单位距离、单位频率的超声波能量的衰减量。
在实施方式3中,首先在步骤1401,对图像生成部36输入从图2的步骤209输出的插补处理后的信号。
接下来,在步骤1402,图像生成部36通过在插补处理后信号的时间轴乘上音速的值(生物体内约1540m/s)来换算成距离。
接下来,在步骤1403,图像生成部36算出波束方向的梯度。该波束方向的梯度表示与超声波的传播相伴的超声波能量的衰减量本身。
进而在步骤1404,图像生成部36根据在步骤1402算出的波束方向的梯度来算出衰减率。若将波束方向的梯度设为D[dB/cm],将发送中所用的频率设为f[MHz],则衰减率α[dB/MHz/cm]表征为
【数学式5】
图像生成部36对用户经由外部输入设备13预先设定的衰减率分布图显示区域内的全区域进行上述处理,作成衰减率分布图。图像生成部36将作成的衰减率分布图显示于显示部16。
图15是衰减率分布图的显示部16中的显示形态的示例。在超声波图像1500上叠加衰减率分布图显示区域1503,显示衰减率分布图。在该显示形态中,用户能将衰减率分布图和代表点1502、无效区域1501进行比较。因而,用户能根据代表点、无效区域的被选出的方法来判断是否适当地进行了衰减率的测量。
图16是衰减率分布图的显示部16中的显示形态的另一个示例。在超声波图像1600上叠加衰减率分布图显示区域1603。衰减率的大小由与以衰减率的显示范围1605确定的色条对应的颜色示出。在该显示形态中,用户能从B像1602与衰减率分布图的比较来判断是否适当检测了无效值1604、适当进行了衰减率的测量。
如上述那样,根据本实施方式1的超声波诊断装置,能在二维空间生成在深度方向以及方位方向上除去散斑的信号(代表值以及插补处理值)。因而,可以使用该信号如实施方式2那样生成B像,也可以如实施方式3那样算出衰减率。
B像由于被精度良好地修正了超声波的传播所引起的衰减,因此成为高精细的图像。因而,用户能看着该高精细的B像来进行诊断。
另外,通过如实施方式3那样算出衰减率并进行显示,在衰减率大到给定值以上的情况下,能有助于脂肪肝等在衰减率表现出特征的疾患的判定。
Claims (12)
1.一种超声波诊断装置,其特征在于,具有:
接收波束形成器,其通过从超声波探头接受信号并进行接收波束形成,来针对设定于摄像区域的多个接收扫描线分别生成接收信号;和
信号处理部,其对所述接收信号进行处理,
所述信号处理部具有:
二维ROI设定部,其在将多个所述接收信号按照所述接收扫描线的顺序进行排列的二维信号空间设定多个二维的ROI;和
代表值算出部,其使用所述ROI内的信号来求取代表值,将求得的代表值作为设定于所述ROI内的代表点的处理后信号值。
2.根据权利要求1所述的超声波诊断装置,其特征在于,
所述二维ROI设定部将多个所述ROI针对所述二维信号空间的沿着所述接收扫描线的方向以及所述接收扫描线的排列方向的至少一方,设定于相互错开的位置。
3.根据权利要求1所述的超声波诊断装置,其特征在于,
所述二维ROI设定部基于所述接收信号来适应性地变更所述ROI的尺寸以及形状的至少一方。
4.根据权利要求1所述的超声波诊断装置,其特征在于,
所述信号处理部具备:ROI适合性分析部,其判定所述ROI内的信号强度的分布是否适合预先确定的信号强度分布,将适合的ROI作为有效ROI,将不适合的ROI作为无效ROI。
5.根据权利要求1所述的超声波诊断装置,其特征在于,
所述信号处理部具备:代表值适合性判定部,其判定所述代表值算出部算出的针对多个所述ROI的所述代表值的分布是否适合预先确定的条件,
所述代表值适合性判定部并不将不适合所述条件的代表值采用为所述处理后信号值。
6.根据权利要求5所述的超声波诊断装置,其特征在于,
所述代表值适合性判定部在所述多个ROI的所述代表值的信号强度的分布沿着声波传播方向慢慢减少的情况下,判定为适合所述条件。
7.根据权利要求1所述的超声波诊断装置,其特征在于,
所述信号处理部具有:图像生成部,其使用所述处理后信号值来生成B模式像。
8.根据权利要求1所述的超声波诊断装置,其特征在于,
所述信号处理部具有:图像生成部,其使用所述处理后信号值来将B模式像的亮度值最佳化。
9.根据权利要求1所述的超声波诊断装置,其特征在于,
所述信号处理部具有:衰减率算出部,其使用所述处理后信号值来求取深度方向的衰减率。
10.根据权利要求1所述的超声波诊断装置,其特征在于,
所述信号处理部具有:插补部,其在所述二维信号空间中使用多个所述ROI的所述代表值,通过插补来求取所述ROI与ROI之间的所述二维信号空间的信号值。
11.一种信号处理装置,其特征在于,具有:
二维ROI设定部,其在将接收波束形成后的接收信号按照接收扫描线的顺序进行排列的二维信号空间设定多个二维的ROI;和
代表值算出部,其使用所述ROI内的信号来求取代表值,将求得的代表值作为设定于所述ROI内的代表点的处理后信号值。
12.一种信号处理程序,其特征在于,使计算机作为如下单元发挥功能:
二维ROI设定单元,其在将接收波束形成后的接收信号按照接收扫描线的顺序进行排列而形成的二维信号空间设定多个二维的ROI;和
代表值算出单元,其使用所述ROI内的信号来求取代表值,将求得的代表值作为设定于所述ROI内的代表点的处理后信号值。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020-066134 | 2020-04-01 | ||
JP2020066134A JP7432426B2 (ja) | 2020-04-01 | 2020-04-01 | 超音波診断装置、信号処理装置および信号処理プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113491535A true CN113491535A (zh) | 2021-10-12 |
CN113491535B CN113491535B (zh) | 2024-04-26 |
Family
ID=77995008
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110190790.4A Active CN113491535B (zh) | 2020-04-01 | 2021-02-19 | 超声波诊断装置、信号处理装置以及记录介质 |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP7432426B2 (zh) |
CN (1) | CN113491535B (zh) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101219063A (zh) * | 2007-01-12 | 2008-07-16 | 深圳迈瑞生物医疗电子股份有限公司 | 基于二维分析的b图像均衡方法和系统结构 |
JP2010082230A (ja) * | 2008-09-30 | 2010-04-15 | Fujifilm Corp | 超音波信号処理装置及び方法 |
CN104622509A (zh) * | 2013-11-08 | 2015-05-20 | 日立阿洛卡医疗株式会社 | 超声波诊断装置以及弹性评价方法 |
JP2015126955A (ja) * | 2015-04-08 | 2015-07-09 | 日立アロカメディカル株式会社 | 超音波診断装置、及び超音波診断方法 |
US20150342568A1 (en) * | 2014-05-29 | 2015-12-03 | Ge Medical Systems Global Technology Company, Llc | Ultrasonic diagnostic apparatus |
WO2016157624A1 (ja) * | 2015-03-31 | 2016-10-06 | オリンパス株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
CN106068099A (zh) * | 2014-12-22 | 2016-11-02 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 |
WO2017068892A1 (ja) * | 2015-10-22 | 2017-04-27 | 株式会社日立製作所 | 超音波診断装置、および、減衰特性計測方法 |
JP2017113145A (ja) * | 2015-12-22 | 2017-06-29 | オリンパス株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
JP2017158917A (ja) * | 2016-03-11 | 2017-09-14 | 東芝メディカルシステムズ株式会社 | 医用診断装置、および医用解析プログラム |
US20180025492A1 (en) * | 2016-07-22 | 2018-01-25 | Toshiba Medical Systems Corporation | Analyzing apparatus and analyzing method |
US20190328364A1 (en) * | 2018-04-30 | 2019-10-31 | Esaote Spa | Method and ultrasound system for shear wave elasticity imaging |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017110756A1 (ja) | 2015-12-24 | 2017-06-29 | オリンパス株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
-
2020
- 2020-04-01 JP JP2020066134A patent/JP7432426B2/ja active Active
-
2021
- 2021-02-19 CN CN202110190790.4A patent/CN113491535B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101219063A (zh) * | 2007-01-12 | 2008-07-16 | 深圳迈瑞生物医疗电子股份有限公司 | 基于二维分析的b图像均衡方法和系统结构 |
JP2010082230A (ja) * | 2008-09-30 | 2010-04-15 | Fujifilm Corp | 超音波信号処理装置及び方法 |
CN104622509A (zh) * | 2013-11-08 | 2015-05-20 | 日立阿洛卡医疗株式会社 | 超声波诊断装置以及弹性评价方法 |
US20150342568A1 (en) * | 2014-05-29 | 2015-12-03 | Ge Medical Systems Global Technology Company, Llc | Ultrasonic diagnostic apparatus |
CN106068099A (zh) * | 2014-12-22 | 2016-11-02 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 |
WO2016157624A1 (ja) * | 2015-03-31 | 2016-10-06 | オリンパス株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
JP2015126955A (ja) * | 2015-04-08 | 2015-07-09 | 日立アロカメディカル株式会社 | 超音波診断装置、及び超音波診断方法 |
WO2017068892A1 (ja) * | 2015-10-22 | 2017-04-27 | 株式会社日立製作所 | 超音波診断装置、および、減衰特性計測方法 |
JP2017113145A (ja) * | 2015-12-22 | 2017-06-29 | オリンパス株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
JP2017158917A (ja) * | 2016-03-11 | 2017-09-14 | 東芝メディカルシステムズ株式会社 | 医用診断装置、および医用解析プログラム |
US20180025492A1 (en) * | 2016-07-22 | 2018-01-25 | Toshiba Medical Systems Corporation | Analyzing apparatus and analyzing method |
US20190328364A1 (en) * | 2018-04-30 | 2019-10-31 | Esaote Spa | Method and ultrasound system for shear wave elasticity imaging |
Also Published As
Publication number | Publication date |
---|---|
JP7432426B2 (ja) | 2024-02-16 |
CN113491535B (zh) | 2024-04-26 |
JP2021159510A (ja) | 2021-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7985182B2 (en) | Ultrasonic diagnostic apparatus and ultrasonic image acquiring method | |
JP5730979B2 (ja) | 超音波診断装置、及び弾性評価方法 | |
US7887485B2 (en) | Ultrasonic image boundary extracting method, ultrasonic image boundary extracting apparatus, and ultrasonic imaging apparatus | |
US11523774B2 (en) | Tissue property estimation with ultrasound medical imaging | |
KR102223048B1 (ko) | 정량적 초음파 이미징을 위한 관심 구역 배치 | |
JP6216736B2 (ja) | 超音波診断装置、及び超音波診断方法 | |
WO2014162966A1 (ja) | 超音波診断装置、及び弾性評価方法 | |
EP1685799B1 (en) | Ultrasonic diagnostic apparatus and ultrasonic image acquiring method | |
US20060241456A1 (en) | Ultrasonic imaging apparatus and ultrasonic imaging method | |
CN102415902B (zh) | 超声波诊断装置以及超声波图像处理装置 | |
EP2208465A1 (en) | Ultrasonographic device, operation method thereof, and ultrasonogram diagnosis program | |
EP2423872A1 (en) | Ultrasonic image processing apparatus | |
US8300909B2 (en) | Ultrasonographic device and ultrasonographic method | |
US20230404537A1 (en) | Ultrasound medical imaging with optimized speed of sound based on fat fraction | |
CN114176639A (zh) | 用于介质的超声表征的方法和系统 | |
CN108135570B (zh) | 超声成像设备和超声成像设备的控制方法 | |
EP2157442B1 (en) | Formation of an elastic image in an ultrasound system | |
US9107631B2 (en) | Ultrasonic imaging apparatus and a method for generating an ultrasonic image | |
KR20210116268A (ko) | 초음파 의료 이미징을 사용한 간 질환 활성도 추정 | |
CN113491535B (zh) | 超声波诊断装置、信号处理装置以及记录介质 | |
JP5823184B2 (ja) | 超音波診断装置、医用画像処理装置および医用画像処理プログラム | |
WO2023274763A1 (en) | Ultrasonic quantification of acoustic attenuation coefficient in the presence of elevation aperture blockage | |
CN114173673A (zh) | 使用图像数据进行的超声系统声学输出控制 | |
US20240156441A1 (en) | Tissue property estimation with ultrasound medical imaging | |
EP3685753A1 (en) | Ultrasound imaging apparatus and method of controlling the same |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
TA01 | Transfer of patent application right | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20211224 Address after: Chiba County, Japan Applicant after: Fujifilm medical health Co.,Ltd. Address before: Tokyo, Japan Applicant before: Hitachi, Ltd. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |