CN103476344B - 被检体信息获取装置 - Google Patents
被检体信息获取装置 Download PDFInfo
- Publication number
- CN103476344B CN103476344B CN201280017094.2A CN201280017094A CN103476344B CN 103476344 B CN103476344 B CN 103476344B CN 201280017094 A CN201280017094 A CN 201280017094A CN 103476344 B CN103476344 B CN 103476344B
- Authority
- CN
- China
- Prior art keywords
- subject
- reference layer
- strain
- measurement reference
- pressure
- 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
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
- A61B8/085—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures for locating body or organic structures, e.g. tumours, calculi, blood vessels, nodules
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0858—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
-
- 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
- A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from 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/58—Testing, adjusting or calibrating the diagnostic device
- A61B8/585—Automatic set-up of the device
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52042—Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Radar, Positioning & Navigation (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Remote Sensing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Vascular Medicine (AREA)
- Gynecology & Obstetrics (AREA)
- Databases & Information Systems (AREA)
- Primary Health Care (AREA)
- Data Mining & Analysis (AREA)
- Epidemiology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
被检体中的粘弹性特性通过简单的方法被成像。在超声波发射/接收探测器与被检体之间包含弹性模量和粘度系数已知的粘弹性测量参照层,并且,从根据施加到被检体的随时间改变的压力在粘弹性测量参照层和被检体中产生的应变的随时间的变化和粘弹性测量参照层的已知的弹性模量和粘度系数,计算被检体内的弹性模量和粘度系数的分布。
Description
技术领域
本发明涉及被检体信息获取装置,特别是涉及用于将被检体的粘弹性特性成像的被检体信息获取装置。
背景技术
超声诊断装置作为被检体信息获取装置是已知的。超声诊断装置以非侵害的方式将被检体内的信息成像,使得超声诊断装置被广泛用于医疗领域中。
在超声诊断装置中,常规上使用超声断面图像和多普勒超声图像,在超声断面图像中,生物体内的结构从由组织的反射性差异产生的反射回声被成像,而在多普勒超声图像中,血液流动速度等通过利用由血液流动产生的超声多普勒效应被成像。
此外,近年来,通过超声波测量组织的硬度,并且,该硬度被用于诊断组织。这是由于组织的硬度与病理状态密切相关。例如,已知诸如乳腺癌和甲状腺癌的硬化性癌比正常组织和良性肿瘤硬。
近年来,据报道,癌与良性肿瘤之间的粘度不同,使得对于组织诊断,不仅需要组织硬度测量,而且需要包含粘度特性评价的粘弹性特性评价。因此,需要通过使用超声波测量生物体内的粘弹性特性的分布、将该分布形成为图像、组合该图像与组织结构的常规图像,并且使用组合的图像以诊断组织内的癌等。
作为用于计算被检体的粘度的超声诊断装置,PTL1公开了如下这样的装置,该装置利用发射和接收超声波的其中在换能器的表面上设置压力传感器的超声探测器、通过超声波测量应变分布、通过压力传感器测量向被检体施加的压力分布,并且从应变分布和压力分布计算弹性和粘度的值。
作为用于测量压力的方法,PTL2公开了其中弹性模量已知的压力测量变形部分被夹在超声探测器和生物组织之间以测量压力的超声诊断装置。超声诊断装置通过超声波测量压力测量变形部分的变形,从弹性模量与应变之间的关系计算施加到压力测量变形部分的压力(应力),从被检体中的应力和应变分布计算被检体中的弹性模量分布,并且显示弹性模量分布。
但是,在PTL1中描述的配置中,超声探测器的结构复杂。而且,一般地,在通过使用超声波测量组织中的应变分布时,施加的压力非常小,使得通常的压力传感器难以测量压力。并且,如果压力传感器紧接在超声探测器下面被***,那么发射和接收超声波的效率和灵敏度劣化,使得存在不能执行深部测量的问题。
在PTL2中,虽然装置的配置简单,但是仅显示被检体中的弹性分布,并且,没有公开用于测量被检体中的粘度分布的方法。因此,需要具有简单的配置并且可计算被检体的粘度的装置。
引文列表
专利文献1
PTL1日本专利公开No.2007-222605
PTL2日本专利公开No.2005-66041
发明内容
可解决以上问题的本发明是一种被检体信息获取装置,该装置用于向被检体发射弹性波、接收从被检体内反射的弹性波并且获取被检体内的信息。被检体信息获取装置包括:转换设备,其接收弹性波并将弹性波转换成电信号;参照层,其设置在转换设备与被检体之间,并且该参照层的粘度系数(ηc)和弹性模量(Ec)是已知;测量单元,其通过使用该电信号测量当向被检体和参照层施加压力时产生的被检体的应变(ε)和参照层的应变(εc);和计算单元,其通过使用参照层的粘度系数(ηc)和弹性模量(Ec)、被检体的应变(ε)和参照层的应变(εc)计算被检体的粘度系数(η)。
本发明的有利效果
本发明可提供诸如超声诊断装置的如下这样的被检体信息获取装置,该被检体信息获取装置即使具有简单的装置配置仍可测量高灵敏性和高分辨率的应变分布、从基于应变分布的应力分布的测量计算杨氏模量和粘度系数的分布、并且将生物体内的粘弹性特性的分布成像。
附图说明
图1是示出作为本发明的被检体信息获取装置的超声诊断装置的框图。
图2A和图2B是示出本发明的粘弹性测量的原理的示图。
图3A~3C是用于解释根据本发明的滞后回线(hysteresis loop)的面积的示图。
图4是用于解释根据本发明的粘弹性分布测量的示图。
具体实施方式
以下,将参照附图详细描述本实施例的实施例。
图1是示出作为本实施例的被检体信息获取装置的超声诊断装置的框图。连同使用超声波获得的被检体的诊断区域的断层图像一起,超声诊断装置显示基于生物组织的粘弹性分布、特别是除了基于弹性模量分布以外还基于粘度系数分布的图像。特别地,在图1所示的例子中,通过由使用后面描述的手持型探测器的操作员利用探测器周期性地按压被检体,测量粘弹性测量参照层和被检体的应变。通过使用粘弹性测量参照层的应变、被检体的应变以及粘弹性测量参照层的已知的弹性模量和粘度系数计算被检体的粘弹性分布,并且显示计算结果。
在图1中,附图标记1表示探测器,该探测器包含转换元件阵列2和粘弹性测量参照层3,该转换元件阵列2是兼作向被检体4发射作为弹性波的超声波的单元和接收作为由被检体4反射的弹性波的超声波并将超声波转换成电信号的单元两者的转换设备,该粘弹性测量参照层3设置在转换元件阵列2与被检体4之间并且其粘度系数(ηc)和弹性模量(Ec)是已知的。附图标记4表示被检体,附图标记5表示向探测器1发射控制信号的发射单元,附图标记6表示从包含于探测器1中的转换设备接收电信号的接收单元。附图标记7是作为如下测量单元的应变分布检测单元,该测量单元通过使用来自包含于探测器1中的转换设备的电信号来测量当向被检体4和粘弹性测量参照层3施加压力时产生的被检体4的应变和粘弹性测量参照层3的应变。附图标记8是作为如下计算单元的粘弹性分布计算单元,该计算单元通过使用粘弹性测量参照层3的测量应变、被检体4的测量应变以及粘弹性测量参照层3的已知的粘度系数和弹性模量计算被检体4的粘度系数。此外,作为实施例,图1的配置包括成像单元9、显示单元10和控制单元11。探测器1还包括控制计算粘弹性分布的操作的控制开关12。以下,将描述通过上述的超声诊断装置计算被检体4的粘度系数的原理。
在本发明的实施例中,探测器1通过手动操作或机械驱动向被检体周期性地施加压力,同时,探测器1通过其弹性模量和粘度系数已知的粘弹性测量参照层3向被检体发射以及从被检体接收作为弹性波的超声波,使得粘弹性测量参照层3中的应变分布和被检体4中的应变分布被测量。
当组织没有粘性时,当向被检体4周期性地施加压力时在被检体4的组织中周期性地产生的应变振幅值由下式表示。
[数学式1]
这里,E是组织的弹性常数(杨氏模量),σ0是周期性地施加的应力(压力)的振幅。因此,当组织不具有粘性时,可从应变振幅值获得组织的弹性。相反,当粘弹性测量参照层3没有粘性时,可从应变量获得应力。
另一方面,例如,当组织具有粘性时,例如,如果施加关于时间具有三角波形的周期压力,那么在组织中产生的周期应变的振幅值由下式表示。
[数学式2]
这里,τ是组织的粘度系数与杨氏模量的比,T是被周期性地施加的压力的周期。如上所述,从应变的测量获得的值是杨氏模量与粘度系数的比,使得不能分别获得杨氏模量与粘度系数。因此,必需在测量通过周期性地按压被检体4而导致的在时间序列上随时间改变的应变和应力的变化的同时,通过使用应变和应力的测量值计算被检体4的粘弹性。此时,由于粘弹性测量参照层3中的粘度,粘弹性测量参照层3的应变的周期性振幅由下式代表。
[数学式3]
可以理解,当周期非常大(T→∞)时,即,当通过花费足够的时间执行测量时,(式2)和(式3)两者都大致逼近(式1)。(因此,如果压力的周期被设为足够长,那么至少可通过使用(式1)测量杨氏模量)。但是,对于正常的组织和粘弹性测量参照层3,杨氏模量与粘度系数的比τ和τC是从几秒到几十秒的时间的。如果使用充分大于该时间的周期,那么测量周期变为诸如几分钟的长时间,使得难以对于实际的测量使用(式1)。(并且,在这种情况下,不能测量粘度系数)。
下面,将描述通过使用在PTL2中描述的技术暂时计算被检体的粘弹性的情况以及基于本实施例的被检体4的粘弹性的计算,同时比较这两个计算。图2A和图2B示出基于PTL2测量被检体的粘弹性的情况和基于本实施例的粘弹性测量的情况的原理。图2A是基于PTL2的被检体的粘弹性测量的动态模型。图2B是基于本发明的被检体4的粘弹性测量的动态模型。在图2A中,粘弹性测量参照层的粘度不被考虑,使得通过如(式1)所示将测量的应变乘以粘弹性测量参照层的已知的杨氏模量Ec,获得施加到被检体的应力σ。如果如PTL2所描述的,当在经过诸如几分钟到大于十几分钟的长时间之后时间变化收敛时仅测量被检体的弹性,那么通过使用粘弹性测量参照层的测量应变获得的应力和通过使用该应力获得的被检体的杨氏模量是根据(式1)的值。但是,由于测量时间太长,实际上难以通过使用手持探测器的超声诊断装置执行这种长时间的测量。此外,作为本质的问题,为了获得被检体的粘度,必须测量应变和应力随时间的变化。具体而言,在经过长时间之后的应变和应力中,不仅粘弹性测量参照层的粘弹性的影响而且被检体的粘弹性的影响消失,使得必须以短的时间间隔测量随时间的变化。在这种情况下,被检体和粘弹性测量参照层各自分别基于它们的粘弹性随时间变化,使得必须根据图2B所示的动态模型执行测量。特别地,当使用具有优异的超声传播特性以及具有为了抑制边界处的反射而与被检体相同的声学特性的部件作为粘弹性测量参照层时,如果粘弹性测量参照层的粘度被忽略,那么在实际的测量中出现大的误差。因此,当使用由实际部件形成的粘弹性测量参照层时,难以仅从粘弹性测量参照层的杨氏模量独立地获得应力。
另一方面,如本实施例中那样,当粘弹性测量参照层3的杨氏模量EC和杨氏模量与粘度系数的比τC已知时,通过使用(式3)从粘弹性测量参照层3的周期应变的振幅εC获得应力振幅σ0。具体地,在粘弹性测量参照层3的位置上测量粘弹性测量参照层3的应变。通过测量应变的最大值和最小值获得应变振幅εC,并且从(式3)获得应力振幅σ0。
当施加更一般的压力时,向组织施加的应力随时间的变化由下式表示。
[数学式4]
这里,σ(t)和εC(t)是粘弹性测量参照层3的应力和应变的随时间的变化。EC和ηC是粘弹性测量参照层3的已知的杨氏模量和粘度系数。
因此,能够通过使用(式4)从在至少一个压力周期上测量的粘弹性测量参照层3的应变获得应力的变化,并且从应力的变化的振幅获得应力振幅σ0。
在本实施例的例子中,在时间序列上测量被检体4和粘弹性测量参照层3中的应变分布。由此,比较被检体4中的特定组织的应变量和粘弹性测量参照层3中的按压该组织的区域的应变量的随时间的变化,并且计算这两个应变量之间的时间偏移(延迟时间)。此时,这两个应变量之间的延迟时间可由下式表示。
[数学式5]
这里,τ是被检体4中的组织的粘度系数与杨氏模量的比,τC是粘弹性测量参照层3的粘度系数与杨氏模量的比。当粘弹性测量参照层3的粘度系数和杨氏模量已知时,从(式5)获得组织的粘度系数和杨氏模量的比τ,并且,通过使用(式2),从在以上的描述中获得的应力振幅σ0和测量的被检体4的应变振幅,计算被检体4的杨氏模量E。此外,从被检体4的杨氏模量E和粘度系数与杨氏模量的比τ,计算被检体4的粘度系数。
此外,将参照图3A、图3B和图3C描述获得被检体4中的组织的粘度系数与杨氏模量的比τ的另一方法。
当在施加周期性压力的同时测量应变和应力并且在曲线图上绘制应变和应力的值时,示出图3A所示的滞后回线。滞后回线的面积S0是粘度系数与杨氏模量的比τ的函数。因此,可从滞后回线的面积获得粘度。此外,可以使用通过粘弹性测量参照层3和被检体4的应变的周期性变化形成的滞后回线,以避免繁琐的应力的随时间的变化的计算。图3B示出在施加周期性压力时绘制粘弹性测量参照层3和被检体4的应变的曲线图。在这种情况下也形成滞后回线,并且,例如,当如上面描述的那样施加关于时间具有三角波形的周期性压力时,滞后回线的面积S1由下式表示。
[数学式6]
(式6)
当将(式2)和(式3)应用于上式时,建立下式。
[数学式7]
(式7)
这里,T是周期性压力的周期,ε和εC分别是通过测量获得的粘弹性测量参照层3的应变振幅和被检体4的应变振幅,并且,τC是粘弹性测量参照层3的已知的粘度系数与杨氏模量的比。因此,可通过除τ以外测量应变来测量左手侧的滞后回线的面积S1和右手侧的各变量的值。因此,可通过逆求解该式以获得τ的值。
特别地,如果周期T被设为约几秒的短时间并且施加周期性压力,那么可使用以下的近似式。
[数学式8]
(式8)
并且,可以使用在(参考文献1)Nitta,Shiina,Ueno,"HysteresisParameter Imaging of Soft Tissue under Quasi-Static Deformation",2003IEEE International Ultrasonics Symposium Proceedings,pp.1606-1609,2003中描述的滞后参数。
图3C是示出如何获得用于计算滞后参数的基准面积S2的示图。
图3C是以与图3B相同的方式在施加周期性压力的同时绘制粘弹性测量参照层3和被检体4的应变的曲线图。计算由被检体4的应变的最大值和最小值包围的面积S2。可通过如下式所示的那样计算滞后回线的面积与基准面积S2的比来获得滞后参数HP。
[数学式9]
可通过以与上述的方式相同的方式使用滞后参数HP获得τ的值。当使用该滞后参数时,通过基准面积将计算标准化,使得存在测量中的误差影响不大的优点。并且,当通过逆求解该式获得τ的值时,存在一般难以产生多值函数并且过程简单的优点。
在这里的描述中,为了简化描述,示出与施加关于时间具有三角波形的周期性压力的情况有关的等式。但是,按压方法可以是周期性重复按压和释放的方法,并且,对于基于图2B所示的模型的解析式和通过模拟获得的拟合表,可使用根据随时间的变化模式在以上的计算中使用的等式。此外,可以使用通过实际测量粘弹性测量参照层获得的拟合表。
作为本发明的另一实施例的测量方法,能够作为周期性地施加压力的替代以阶段状方式施加压力,并且在时间序列中测量被检体4和粘弹性测量参照层3中的应变分布。瞬时施加压力,然后在时间序列中测量被检体4和粘弹性测量参照层3中的应变分布。即使当压力被瞬时施加时,粘弹性测量参照层3的粘度的影响仍在经过足够的时间之后被去除。但是,由于测量时间太长,因此这是不现实的。因此,以与周期性地施加压力的情况相同的方式在短的时间间隔内测量应变值的随时间的变化。
此时,在图2B所示的模型中,粘弹性测量参照层3的应变的随时间的变化ε′(t)和被检体4中的应变的随时间的变化εC′(t)由下式表示。
[数学式10]
[数学式11]
通过使用(式10),从粘弹性测量参照层3的测量应变以及已知的杨氏模量和粘度系数,获得应力振幅σ0。此外,从测量的被检体4的应变变化获得粘度系数与杨氏模量的比τ。具体地,通过下式计算被检体4的应变的时间微分的对数。
[数学式12]
因此,从时间微分的对数,求得作为在多个不同时间点的被检体4的应变的时间微分的对数的斜率的粘度系数与杨氏模量的比τ的倒数。当比τ的倒数被应用于(式11)时,获得被检体4的杨氏模量,并且,从杨氏模量和比τ计算被检体4的粘度系数。
将参照图4描述用于将杨氏模量和粘度系数的测量(计算)应用于被检体4的粘弹性分布测量(计算)的方法。在图4中,附图标记3指示粘弹性测量参照层,附图标记4指示被检体,附图标记13指示对于测量而施加的压力。压力13被基本上垂直地施加到粘弹性测量参照层3与被检体4之间的边界,并且,在粘弹性测量参照层3与被检体4中的组织中在图4中的垂直方向(压力13的方向)上出现应变。在粘弹性测量参照层3中,沿边界出现εC1、εC2和εC3等的应变分布,并且,在被检体4中,出现ε11、ε12和ε13等的应变分布。通过作为弹性波的超声波,测量这些应变分布的随时间的变化。在图4中的垂直方向上施加压力13,使得从同一应力产生在垂直方向上对准的粘弹性测量参照层3的应变εC1以及被检体4的应变ε11、ε12和ε13等。图2B的模型被应用于粘弹性测量参照层3的应变εC1和被检体4的应变ε11的随时间的变化,使得如上面描述的那样计算应变ε11的位置处的杨氏模量和粘度系数。此外,从粘弹性测量参照层3的应变εC1和应变ε12的随时间的变化,计算应变ε12的位置处的杨氏模量和粘度系数。以相同的方式,计算在垂直方向上与粘弹性测量参照层3的应变εC1对准的位置处的杨氏模量和粘度系数。通过使用在垂直方向上与粘弹性测量参照层3的应变εC1、εC2和εC3等的应变分布的各位置对准的被检体4的应变的值执行相同的方法,使得可计算断面的粘度系数和杨氏模量。
基于上述的原理,其粘度系数和弹性模量已知的粘弹性测量参照层3被设置在转换元件阵列2与被检体4之间,该转换元件阵列2为接收由被检体4反射的弹性波并将弹性波转换成电信号的转换设备,并且,基于由被检体4反射的弹性波,测量当向被检体4和参照层3施加压力时产生的被检体4的应变和参照层3的应变。能够通过使用参照层3的测量应变、被检体4的测量应变以及参照层3的已知的粘度系数和弹性模量计算被检体的粘度系数。
如上所述,向被检体4和参照层3施加的压力是包含压力的释放的周期性压力,使得作为计算单元的粘弹性分布计算单元8可通过使用通过该周期性压力产生的被检体4的应变值的随时间的变化和参照层3的应变值的随时间的变化计算被检体4的粘度系数(η)。
具体地,例如,如上式6~8所示,作为计算单元的粘弹性分布计算单元8可通过使用用通过周期性压力产生的被检体4的应变值和参照层3的应变值绘制的滞后回线的面积或者如上式5描述的那样通过使用被检体4的应变值的随时间的变化与参照层3的应变值的随时间的变化之间的延迟时间,计算被检体的粘度系数(η)。
如上所述,向被检体4和参照层3施加的压力瞬时增加到一定压力,然后,该一定压力被保持一定时间段,使得作为计算单元的粘弹性分布计算单元8可通过使用在向被检体4和参照层3施加的压力瞬时增加到一定压力之后的被检体4的应变值和参照层3的应变值,计算被检体4的粘度系数(η)。
虽然在以上的实施例中使用杨氏模量作为弹性特性并且使用粘度系数作为粘度特性,但是,可以使用诸如刚度系数和压力弹性系数的弹性系数和诸如粘度和粘度系数的粘性系数。此外,这里,将描述使用作为弹性特性的存储弹性模量以及作为粘度特性的损失弹性模量的粘度特性的计算方法。当向具有粘弹性特性的材料施加应力σ(t)时,应变ε(t)的相位相对于应变的相位延迟。
σ(t)=σ0sinωt
ε(t)=ε0sin(ωt-δ)
δ表示相位的延迟。此时,存储弹性模量G′、损失弹性模量G′′和绝对动态弹性模量|G|可由下式表示。
G′=|G|cosδ
G"=|G|sinδ
这里,当粘弹性测量参照层的应变振幅为εC0、被检体的应变振幅为εt0且粘弹性测量参照层的应变与被检体的应变之间的相位差为δt-C时,可通过超声波测量这三个值。粘弹性测量参照层的存储弹性模量GC′、粘弹性测量参照层的损失弹性模量GC′′、绝对动态弹性模量|Gc|和粘弹性测量参照层的应变的相位δC被分别测量,使得这些值变为已知的参数。当使用可通过超声波测量的值和已知的参数时,可通过下式计算被检体的应变的相位δt、被检体的存储弹性模量Gt′、被检体的损失弹性模量Gt′′。
δt=δt-c+δc
换句话说,可通过使用粘弹性测量参照层的弹性特性(存储弹性模量)、粘弹性测量参照层的粘度特性(损失弹性模量)和粘弹性测量参照层和被检体的应变,计算被检体的粘度特性(损失弹性模量)。如上面描述的那样计算的存储弹性模量和损失弹性模量是可通过通常的粘弹性特性测量设备测量的值。例如,可直接将这些值与通过使用提取的组织测量的粘弹性的机械测量结果相比较。
下面,将基于例子1连同以上的实施例的配置要素描述本发明的例子。
例子1
以下,将参照附图描述本发明的例子1。图1是示出根据本发明的超声诊断装置的例子的框图。配置要素将连同它们的操作一起被依次描述。
[粘弹性测量参照层的特定配置例]
粘弹性测量参照层3由如下这样的材料形成,该材料传播但是不吸收或者散射作为弹性波的超声波并且具有保形性、柔性和适当的弹性。由于通过粘弹性测量参照层3发射和接收超声波形成图像,因此超声波的传播性能重要。如果超声波传播由于吸收或散射而衰减,那么发射和接收的效率劣化。因此,粘弹性测量参照层3由对于超声波透明并且具有小的衰减常数以及被检体4中的平均声速附近的传播速度(声速)的均匀材料形成。如果粘弹性测量参照层3的声学阻抗与被检体4明显不同,那么在边界表面上出现反射,并且,用于产生图像的发射和接收的效率劣化,使得粘弹性测量参照层3的声学阻抗与被检体4的声学阻抗的差异不那么大。但是,如果粘弹性测量参照层3的声学阻抗与被检体4的声学阻抗相同,那么如后面描述那样难以检测边界表面,使得粘弹性测量参照层3与被检体4的声学阻抗相差适量的声学阻抗。粘弹性测量参照层3具有保形性、柔性和适当的弹性,以通过利用粘弹性测量参照层3的变形量测量应力(压力)的分布。作为具有以上特性的材料,可以使用诸如聚乙烯醇、聚亚安酯的含水凝胶和橡胶材料。但是,可以使用具有以上特性的任何材料。但是,使用诸如杨氏模量、刚度系数和压力弹性系数的弹性系数和诸如粘度和粘度系数的粘性系数已知的材料。从传播超声波的观点看,粘弹性测量参照层3的厚度应薄。但是,粘弹性测量参照层3应适当地变形以测量应力(压力)的分布。此外,粘弹性测量参照层3必须被置于应变分布被检测到的区域中。因此,粘弹性测量参照层3的厚度为0.1mm~50mm、优选为1mm~10mm。
[产生接收束信号并显示B模式图像]
探测器1以机械或电子的方式执行射束扫描,并且向被检体4发射以及从被检体4接收超声波。转换元件阵列2是在其中布置由来自发射单元5的驱动波形驱动并产生超声波的转换元件的转换设备。发射单元5生成用于转换元件阵列2的各转换元件的驱动波形,并基于来自控制单元11的控制调整各转换元件的驱动定时。通过合成通过调整转换元件的驱动定时从各转换元件产生的超声波而获得的合成超声波形成会聚于预定点的超声波发射束。换句话说,控制单元11控制发射单元5,使得从探测器1发射具有在希望的深度的会聚点以及在希望的方向上的方向性的超声波发射束。控制单元11执行在其中摆动超声波发射束的方向并执行射束扫描的扇形扫描。或者,在转换元件阵列2上被驱动的转换元件被限制以形成发射开口并且发射开口移动,使得执行其中超声波发射射束基本上平行地扫描的线性扫描。
接收单元6包含接收电路和整相相加电路(phasing additioncircuit)。接收电路放大通过转换元件阵列2的各转换元件接收的接收信号,并且通过执行AD转换将接收信号转换成多个时间序列数字信号。整相相加电路是要形成超声波接收束。整相相加电路对于数字化的接收信号提供由控制单元11控制的延迟时间并进行相加(整相相加)以形成接收束信号。接收束信号是如下这样的信号,在该信号中通过转换元件接收的信号的接收定时被调整并且所接收的信号被求和,并且,能够形成关于接收器灵敏度具有方向性并具有一个或更多个收敛点的超声波接收束。控制单元11与超声波发射束同步地执行超声波接收束的射束扫描。作为结果,通过从探测器1将超声波发射束发射到被检体4中并接收从被检体4的内部组织反射的超声反射回声信号在超声波接收束中获得接收束信号。通过控制单元11的射束扫描,获得与对于被检体4中的特定点或特定方向的超声反射回声信号对应的多个接收束信号。接收束信号与射束扫描对应地布置,使得获得与被检体4中的超声反射回声信号的强度对应的断面图像。控制单元11重复射束扫描,使得获得一天中的不同时间的断面图像。
通过接收单元6产生的多个接收束信号被输入到成像单元9中。成像单元9从多个接收束信号形成反映被检体4中的断面结构的灰度级B模式断面图像。成像单元9根据射束扫描布置多个接收束信号并对于接收束信号执行诸如增益校正、滤波处理、包络解调和对数压缩的信号处理。此外,成像单元9可包含用于通过执行诸如轮廓强调(upsharp mask)和图像滤波的图像处理将信号转换成供给到显示单元10的显示信号以用于信号处理之后的图像显示的数字扫描转换电路和用于模拟视频信号的DA转换电路。成像单元9导致显示单元10根据来自控制单元11的显示控制显示B模式断面图像。
显示单元10是诸如CRT和LCD的显示装置。显示单元10显示基于来自成像单元9的显示信号的图像。
[作为用于测量应变的测量单元的应变检测单元的特定配置例子]
通过接收单元6产生的多个接收束信号还被输入到应变分布检测单元7中。作为用于测量应变的测量单元的应变分布检测单元7基于多个接收束信号检测被检体4和粘弹性测量参照层3中的应变分布。
可通过以与已知的颜色流动多普勒方法和组织跟踪方法相同的方式执行在不同的射束扫描时段获得的在相同位置处的接收束信号之间的相关性计算以检测应变分布。
应变分布检测单元7具有用于存储多个接收束信号的存储器电路和相关性计算电路。当控制单元11重复射束扫描时获得的多个接收束信号存储于存储器电路中。通过一个射束扫描获得的一组接收束信号形成与一定时间点的断面图像对应的帧接收束信号数据。不同的帧接收束信号数据与不同时间点的射束扫描对应(各帧接收信号数据包含来自多个不同的位置的接收束信号)。不同的帧接收束信号数据之间的相应位置处的接收束信号是在不同的时间点的射束扫描之间来自同一位置的超声反射回声信号。因此,可通过执行这些接收束信号的相关性计算测量在不同时间点的射束扫描之间产生的相应位置处的位移。通过获取沿射束深度方向的差值从位移分布计算应变分布。
相关性计算电路在其输入部分具有Hilbert变换滤波器。相关性计算电路将与相同位置对应的不同时间点的射束扫描的两个接收束信号转换成解析信号,并执行这些信号的复相关性计算。或者,通过使用通过接收束信号的正交检测获得的I-Q信号,计算与相同位置对应的不同时间点的射束扫描的两个接收束信号之间的时间序列中的瞬时相位延迟,并且,可通过将瞬时相位延迟转换成声速来获得位移。此外,可通过使用在相关技术领域中已知的组合自相关方法等提高应变测量的精度。
作为存储器电路和相关性计算电路的具体配置,可以使用以下描述的配置。
第一例子包括存储多个帧接收束信号数据的存储器单元和与在一次射束扫描中使用的射束位置对应的多个相关性计算电路。从存储于存储器电路中的多个帧接收束信号数据选择与不同的时间点的射束扫描对应的两个帧接收束信号数据,并且,与该两个帧接收束信号数据对应的接收束信号被输入到与各射束位置对应的相关性计算电路中。可通过在相关性计算电路中执行相关性计算,获得对应束位置处的在射束深度方向上的位移量分布。在该配置中,重复执行射束扫描,并且,通过从多个获得的帧接收束信号数据选择帧接收束信号数据,执行相关性计算,使得能够测量任何两个时间点的射束扫描之间的位移。可通过计算深度方向上的应变之间的差异来测量应变量。
这里,使用对于各射束位置执行一维相关性计算的多个相关性电路。但是,能够使用对于射束扫描的不同时间点的帧接收束信号数据执行二维相关性计算的相关电路。
另一例子包括与在射束扫描中使用的各射束位置对应的多个存储器电路和多个相关性计算电路。多个存储器电路中的每一个包含两个FIFO型存储器和用于初次存储相关性计算结果的TMP存储器。当重复执行射束扫描时,首先对于各射束位置在第一FIFO型存储器中存储接收束信号。第一FIFO型存储器的输出被输入到第二FIFO型存储器和相关性计算电路中。此外,第二FIFO型存储器的输出被输入到相关性计算电路中。两个FIFO型存储器中的每一个具有用于存储在与各射束扫描对应的射束位置处获得的接收束信号的存储容量。因此,在重复的射束扫描的某一时间点,存储与两个时间连续的射束扫描对应的射束位置处的两个接收束信号。通过相关性计算电路计算两个接收束信号,并且,获得相应的射束位置处的在射束深度方向上的位移量分布。在重复的射束扫描中的两个时间连续的射束扫描之间产生这里获得的位移量分布。为了获得在多个射束扫描之间产生的位移,累积相加相关性计算电路的计算结果。因此,相关性计算的计算结果与TMP存储器中的数据进行相加,并且,相加结果被存储于TMP存储器中。此外,通过计算深度方向上的位移分布的差异,计算应变分布。该配置适于通过适合于实时处理的流水线型处理来实现应变分布检测单元7。特别地,在该配置中,存储器电路的容量可比上述的例子小。
在上述的两个例子中,从与各射束位置对应的多个相关性计算电路获得在不同的射束扫描之间产生的在射束深度方向上的应变分布。应变分布被布置于各射束位置中,使得获得在不同的射束扫描之间产生的二维断面图像的应变分布。除了上述的两个例子以外,还能够通过使用在相关技术中已知的方法、使用在不同时间点的射束扫描获得的帧接收束信号数据,检测二维断面的应变分布。
在本例子中,特别地,使用不仅包含被检体4而且包含粘弹性测量参照层3的区域以测量应变分布。更具体地,该区域被设定为使得测量被检体4与粘弹性测量参照层3之间的边界部分中的位移量。由此,计算粘弹性测量参照层3的应变分布。为了检测粘弹性测量参照层3中的位移量,切除接收束信号,使得接收束信号包含来自边界部分的超声波的反射回声,并且可执行上述的一系列的相关性计算。但是,可单独地设置以下描述的计算电路。
一般地,在粘弹性测量参照层3内部基本上对于超声波透明,使得超声波的反射回声信号非常小,并且,接收束信号也非常小。因此,该区域中的接收束信号的相关性计算结果包含误差,并且结果可能不稳定。另一方面,由于被检体4与粘弹性测量参照层3之间的声学阻抗的差异,被检体4与粘弹性测量参照层3之间的边界部分在本地产生大的反射回声信号,使得可准确地检测接收束信号及其相关性计算结果。粘弹性测量参照层3由厚度方向上均匀的材料形成,并具有在厚度方向上恒定的应变分布,使得通过使用从边界部分中的位移获得的应变值获得层中的应变分布。
因此,阈值电路被设置在相关性计算电路的输入部分,并且,其强度比一定值小的接收束信号被设为0,使得通过掩蔽各接收束信号的前部部分和与粘弹性测量参照层3的厚度对应的部分来执行相关性计算。因此,可以消除上述的不稳定。
此外,接收束信号的振幅强度比一定值小的部分被确定为在粘弹性测量参照层3的内部,使得直接检测到接收束信号的边界部分,并且,可从射束扫描的帧接收束信号数据之间的差获得边界部分位置不同的射束扫描之间的边界部分中的位移。可通过沿深度方向取得差值从被检体4与粘弹性测量参照层3之间的边界部分的位移分布计算粘弹性测量参照层3的应变分布。
如上所述,从各接收束信号的前部部分计算粘弹性测量参照层3的应变量,并且,从在该前部部分之后的部分计算沿深度方向的被检体4的应变分布。从接收束信号计算的沿深度方向的应变分布被布置于相应的射束位置处,使得测量被检体4的断面应变分布和粘弹性测量参照层3的沿粘弹性测量参照层3与被检体4之间的边界的应变分布。这些应变分布从应变分布检测单元7被输入到作为用于计算粘度系数的单元的粘弹性分布计算单元8。
[作为用于计算粘度系数的单元的粘弹性分布计算单元]
作为用于计算粘度系数的单元的粘弹性分布计算单元8根据上述的原理从粘弹性测量参照层3和被检体4的应变分布以及粘弹性测量参照层3的弹性系数和粘性系数计算被检体4中的弹性模量分布和粘度系数分布。
粘弹性分布计算单元8可通过使用微处理器、存储器、控制微处理器和存储器的控制电路以及传输数据的总线电路形成,或者可实现为诸如通用PC和FPGA的可编程处理电路上的处理软件。
因此,以下将描述粘弹性分布计算单元8中的处理。当执行超声束的线性扫描以测量应变时,与各接收束信号对应的超声束的各位置基本上与施加压力的方向对应。因此,对于各接收束信号获得的粘弹性测量参照层3的应变量和被检体4的沿深度方向的应变分布与在图4中描述的垂直方向上对准的同一列中的应变对应。因此,可对于各接收束信号执行包含应变分布检测单元7的处理的处理。当执行超声束的扇形扫描时或者当执行另一扫描时,通过利用执行扫描时的超声束的位置与施加压力的方向之间的几何关系,对于在图4中描述的垂直方向上对准的各列,产生粘弹性测量参照层3的应变量和被检体4的沿深度方向的应变分布的布置。此时,通过内插计算从超声束的位置与在垂直方向上对准的列之间的位置关系产生被检体4的应变分布的布置。为了简化描述,以下将描述执行线性扫描的例子。但是,可通过由从以上的位置关系用在垂直方向上对准的数据串替代对于各接收束信号产生的粘弹性测量参照层3的应变量和被检体4的沿深度方向的应变分布的布置,对于其它的射束扫描方法执行相同的处理。
对于通过一次射束扫描获得的各帧接收束信号数据,获得在射束扫描的帧时间的断面应变分布数据组。断面应变分布数据组包含从在对应于一次射束扫描的时间的各接收束信号获得的粘弹性测量参照层3的应变量和被检体4的沿深度方向的应变值的串。对于各接收束信号的粘弹性测量参照层3的应变量和被检体4的沿深度方向的应变值的数据组被称为束应变分布数据组。断面应变分布数据组包含与在各时间的各接收束信号对应的多个束应变分布数据组。在断面应变分布数据组中,各束应变分布数据组与射束扫描中的超声束位置对应。可通过重复射束扫描来获得与不同时间的扫描对应的多个断面应变分布数据组。
可通过在周期性地施加压力的同时重复射束扫描来获得在各时间的断面应变分布数据组。
从在完成压力施加的一个循环或多个循环时获得的多个断面应变分布数据组,计算粘弹性测量参照层应变-被检体应变滞后回线的面积S1。对于各超声束位置从相应的束应变分布数据组计算粘弹性测量参照层3的最大应变εcmax和最小应变εcmin。(从在压力施加的一个循环期间获得的各时刻的断面应变分布数据组提取与希望的超声束位置对应的束应变分布数据组,并且,获得数据组中的粘弹性测量参照层3的应变的最大值和最小值。对于各超声束位置执行以上的操作。可通过使用在压力施加的多个循环期间获得的断面应变分布数据组计算压力施加的一个循环中的粘弹性测量参照层3的应变值的平均值,并且可获得一个循环中的最大值和最小值。作为最大值和最小值,可以使用束应变分布数据组中的值,即,在各扫描获得数据的时刻所获得的值。或者,值可关于时间在这些值之间被内插,并且,在实际获得数据的时刻之间的时刻的内插值可被估计为最大值或最小值。)对于各超声束位置,从最大应变εcmax和最小应变εcmin之间的差值计算粘弹性测量参照层应变振幅εc。然后,对于各超声束位置的在深度方向上的各点计算粘弹性测量参照层应变-被检体应变滞后回线的面积S1。首先,从沿提供最小应变εcmin时到提供最大应变εcmax时的时刻的被检体的应变值和粘弹性测量参照层的应变值获得外出面积(outwardarea)。(对于束应变分布数据中的各深度的点提供被检体的应变值,并且,对于点所属于的各束应变分布数据提供粘弹性测量参照层的应变值。因此,获得束位置的各深度处的点的各时刻的一对被检体的应变值和粘弹性测量参照层的应变值。例如,从各时刻的这些值的对获得梯形积分近似或Simpson积分近似和面积S1。)然后,类似地,沿从提供最大应变εcmax时到提供最小应变εcmin时获得返回面积,并计算外出面积与返回面积之间的差,从而获得粘弹性测量参照层应变-被检体应变滞后回线的面积S1。同时,在各点计算被检体的最大应变εmax和最小应变εmin,并且,从它们之间的差计算粘弹性测量参照层应变振幅ε。然后,通过逆求解(式8)计算各点的弹性模量和粘度系数。
可通过重复压力施加的多个循环提高用于计算滞后回线面积的数据的精度。具体而言,通过使用在各点和各时间获得的所有的粘弹性测量参照层应变和被检体应变的对来计算面积。在这种情况下,使用在粘弹性测量参照层应变随时间增加的过程中获得的对来计算外出面积,并且,使用在粘弹性测量参照层应变随时间减小的过程中获得的对来计算返回面积。关于最大应变和最小应变,类似地,可以使用跨着整个测量时间被测量的应变的最大值和最小值。可对于各压力循环计算粘弹性测量参照层应变-被检体应变滞后回线的面积S1,并且,可使用面积S1的平均值。
此外,计算粘弹性测量参照层3的应变量的一部分或束应变分布数据组中的平均应变值,并且,可通过测量平均应变值的随时间的变化来计算压力循环时间T。换句话说,可通过计算在施加周期性压力的过程中随时间改变的平均应变值的最大值和最小值的时间点之间的时间差,获得压力循环时间T。当施加多个压力循环时,可通过计算压力循环时间的平均值,获得压力循环时间T。虽然粘弹性测量参照层3的应变量的一部分或平均被用作平均应变值,但是可以使用被检体4中的特定位置的应变量或被检体4中的应变量的平均。
通过对于超声束位置的深度方向上的所有点以及对于射束扫描中的所有超声束位置执行以上的计算,计算被检体的断面中的弹性模量和粘度系数的分布。
所计算的弹性模量和粘度系数的分布数据被输入到成像单元9中。成像单元9从多个的弹性模量和粘度系数的分布数据形成反映被检体4中的粘弹性的粘弹性断面图像。以与形成B模式断面图像相同的方式,成像单元9执行诸如增益校正和滤波处理的信号处理和诸如轮廓强调和图像滤波的图像处理以形成粘弹性断面图像。此外,成像单元9根据来自控制单元11的显示控制使得显示单元10通过诸如切换、同时并列和重叠合成的样式将粘弹性断面图像连同B模式断面图像一起显示。此外,成像单元9可基于由应变分布检测单元7产生的被检体4中的应变分布数据产生应变断面图像,并且根据来自控制单元11的显示控制使得显示单元10通过诸如切换、同时并列和重叠合成的样式将粘弹性断面图像连同B模式断面图像一起显示。弹性模量和粘度系数可分别被分配给色调、色度和亮度中的任一种并被显示。通过这种显示,能够直观地观察诸如弹性模量和粘度系数的这两种物理性能。并且,当在B模式断面图像上重叠弹性模量和粘度系数时,可设定和重叠反映弹性模量和粘度系数中的一个或两者的透明度。通过这种重叠,能够有效地观察具有明显的弹性模量或粘度系数的区域。
[粘弹性分布测量操作和粘弹性断面图像的显示]
下面,将描述当执行粘弹性分布测量时由控制单元11执行的处理。如上所述,控制单元11执行超声波发射/接收扫描的控制,并导致显示单元10显示B模式断面图像。此时,根据输入,控制单元11可导致应变分布检测单元7操作,并且连同B模式断面图像一起或者作为其替代在显示单元10上显示应变断面图像。
当用户在检查显示单元10上的B模式断面图像的同时在被检体4的预定部分上设置探测器1并且接通控制开关12时,控制单元11开始粘弹性分布计算单元8的处理,并且,粘弹性分布测量操作开始。当用户通过探测器1向被检体4施加周期性压力时,粘弹性分布计算单元8操作并且产生各时间点的断面应变分布数据组。当施加多个压力循环时,用户在循环的开始和循环的结束从控制开关12输入循环周期信号。控制单元11从循环周期信号的输入计算压力周期T,并且将压力周期T输入到粘弹性分布计算单元8中。当完成周期性压力施加时,通过关断控制开关12,粘弹性分布计算单元8中的断面应变分布数据组的产生结束,计算粘弹性测量参照层应变-被检体应变滞后回线的面积S1,并且,计算弹性模量和粘度系数的分布。当完成计算处理时,控制单元导致成像单元9操作并且生成粘弹性断面图像,并且导致显示单元10根据事先输入的显示样式显示粘弹性断面图像。
虽然这里从控制开关12输入压力循环的开始和结束,但是,如果从上述的粘弹性分布计算单元8中的获得的断面应变分布数据组中的应变中的一个的随时间的变化获得压力周期T,那么不必须输入循环周期信号。
此外,虽然在获得所有的断面应变分布数据组之后计算粘弹性测量参照层应变-被检体应变滞后回线的面积S1,但是,能够规定数据获取压力循环的数量并在完成规定数量的压力循环之后开始计算粘弹性测量参照层应变-被检体应变滞后回线的面积S1、弹性模量和粘度系数。在更新要使用的断面应变分布数据组的同时计算滞后回线的面积、弹性模量和粘度系数,使得可随时间计算弹性模量和粘度系数的分布,同时可对于各压力循环更新粘弹性断面图像。
虽然在以上的描述中通过使用粘弹性测量参照层应变-被检体应变滞后回线的面积S1计算弹性模量和粘度系数,但是,可通过使用滞后参数计算弹性模量和粘度系数。
能够计算粘弹性测量参照层应变振幅εC并通过使用(式3)计算与各束应变分布数据对应的束位置处的应力振幅σ0。然后,可从束位置处的粘弹性测量参照层应变和束位置中的深度方向上布置的被检体中的应变的随时间的变化计算各点处的粘弹性测量参照层应变的变化和被检体中的应变变化之间的延迟时间,并且可通过使用(式5)计算弹性模量和粘度系数。(可从被检体的应变值和粘弹性测量参照层的应变值分别达到最大值和/或最小值的时间差获得延迟时间。这里,对于被检体中的各点获得被检体的应变值,并且,使用关于点所属于的束位置获得的粘弹性测量参照层的应变值。)或者,可通过对关于各点改变的粘弹性测量参照层的应变和被检体的应变的时间序列数据执行互相关性计算,获得延迟时间。
虽然在本例子中,用户通过使用手持型的探测器施加周期性压力,但能够通过将诸如微电动机的激励单元附接于手持型探测器或者设置用于按压探测器的单独的按压单元并从控制单元11驱动激励单元或单独的按压单元,自动施加周期性压力。在该配置中,能够精确地控制压力施加的随时间的变化,诸如周期性施加的压力周期T,使得可以提高弹性模量和粘度系数的计算精度。
此外,虽然在说明书中在探测器1中包含粘弹性测量参照层3,但是,粘弹性测量参照层3可形成为与探测器1分开的可附接和可拆卸的耦合器。可根据要被诊断的区域选择耦合器的形状和粘弹性。在这种情况下,在控制单元11中,可根据耦合的类型输入或选择弹性模量和粘度系数,并且,粘弹性分布计算单元8可通过使用选择的值计算被检体的弹性模量和粘度系数的分布。
例子2
在本例子中,不是周期性地施加压力,而是阶梯状地施加压力。与例子1的不同主要是粘弹性分布计算单元8的处理,从而将主要描述粘弹性分布计算单元8的处理,并且仅描述不同之处。
当探测器1的控制开关12被接通时,控制单元11控制粘弹性分布计算单元8,使得与射束扫描的各时间点对应地存储断面应变分布数据组。探测器1通过手动操作或者按压单元以阶梯状方式被瞬时按压到一定压力水平,并且,保持压力,直到完成数据获取。粘弹性分布计算单元8使用(式10)以及此后在不同的时间点获得的各束应变分布数据组的粘弹性测量参照层3的应变的随时间的变化,并计算相应的束位置的应力振幅σ0(当使用式10时,可通过将根据阶梯状压力施加的在特定的时间间隔之后获得的粘弹性测量参照层3的应变值、粘弹性测量参照层3的已知的杨氏模量以及从杨氏模量和已知的粘度系数获得的τC的值应用于式(10),获得应力振幅σ0)。但是,希望通过使用多个点的时间间隔并拟合到(式10)来提高测量精度。
通过使用在不同的时间点获得的多个断面应变分布数据组中的各束应变分布数据组中的被检体应变值,提取相应的束位置处的深度方向上的各位置处的被检体应变的随时间的变化。通过使用(式12)从各点的随时间的变化计算被检体中的各点处的杨氏模量和粘度系数的比τ。首先,为了计算被检体应变的微分,通过将与各断面应变分布数据组对应的根据阶梯状压力施加的各特定的时间间隔的被检体应变的增量除以时间间隔,计算应变的近似时间微分值,并且,根据(式12)从应变的近似时间微分值的随时间的变化的斜率获得1/τ。此外,此时,希望通过拟合从多个时间点的应变的近似微分值获得1/τ。同时,从应变的近似时间微分值的随时间的变化的截距(intercept)获得下式的值。
[数学式13]
然后,从在以上的计算中获得的τ和应力振幅σ0获得杨氏模量E。
通过对于超声束位置的深度方向上的所有点以及对于射束扫描中的所有超声束位置执行以上的计算,计算被检体的断面中的弹性模量和粘度系数的分布。
所计算的弹性模量和粘度系数的分布以与例子1相同的方式被输入到成像单元9中,并且,在显示单元10上显示粘弹性断面图像。
虽然这里通过将被检体应变的增量除以时间间隔获得近似时间微分值,但可以使用高次中心差分、前向差分或后向差分。
在本例子中,针对阶梯状压力施加从被检体和粘弹性测量参照层中的应变的随时间的变化获得被检体的杨氏模量和粘度系数。在这种情况下,关于必要的时间序列时间间隔,可通过使用几个到几十个时间点的数据(包含上述的对于拟合使用该数据的情况)执行足够精确的计算。由此,对于各时间点的断面应变分布数据组的必要数量是几个到几十个,并且,用于产生断面应变分布数据组的帧接收束信号数据包含几个到几十个帧。在一到几秒内完成用于产生这些数据的射束扫描,使得存在在短的时间段内获取用于产生粘弹性分布数据的数据的优点。并且,存在可在短的时间段内获取测量数据而无需执行在例子1中执行的压力施加和压力释放的优点。
虽然已参照示例性实施例说明了本发明,但应理解,本发明不限于公开的示例性实施例。以下的权利要求的范围应被赋予最宽泛的解释以包含所有这样的变更方式和等同的结构和功能。
本申请要求在2011年4月8日提交的日本专利申请No.2011-086512的权益,通过引用将其全部内容并入此。
Claims (5)
1.一种被检体信息获取装置,该被检体信息获取装置用于向被检体发射弹性波、接收从被检体内部反射的弹性波并且获取被检体内部的信息,该被检体信息获取装置包括:
转换设备,被配置为接收弹性波并将弹性波转换成电信号;
参照层,设置在转换设备与被检体之间,并且该参照层的粘度系数和弹性模量是已知的;
测量单元,被配置为通过使用所述电信号测量当向被检体和参照层施加压力时产生的被检体的应变和参照层的应变;和
计算单元,被配置为通过使用参照层的粘度系数和弹性模量、被检体的应变以及参照层的应变,计算被检体的粘度系数。
2.根据权利要求1的被检体信息获取装置,其中,向被检体和参照层施加的压力是包含释放压力的周期性压力施加,并且,
计算单元通过使用通过所述周期性压力施加而产生的被检体的应变值的随时间的变化和参照层的应变值的随时间的变化,计算被检体的粘度系数。
3.根据权利要求2的被检体信息获取装置,其中,计算单元通过使用由通过周期性压力产生的被检体的应变值和参照层的应变值绘制的滞后回线的面积,计算被检体的粘度系数。
4.根据权利要求2的被检体信息获取装置,其中,计算单元通过使用通过所述周期性压力施加产生的被检体的应变值的随时间的变化和参照层的应变值的随时间的变化之间的延迟时间,计算被检体的粘度系数。
5.根据权利要求1的被检体信息获取装置,其中,
向被检体和参照层施加的压力瞬间增加到一定压力,然后该一定压力被保持一定时间段,并且,
计算单元通过使用在施加到被检体和参照层的压力瞬间增加到一定压力之后的被检体的应变值和参照层的应变值,计算被检体的粘度系数。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011086512 | 2011-04-08 | ||
JP2011-086512 | 2011-04-08 | ||
PCT/JP2012/059952 WO2012137983A1 (en) | 2011-04-08 | 2012-04-05 | Subject information acquisition apparatus |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103476344A CN103476344A (zh) | 2013-12-25 |
CN103476344B true CN103476344B (zh) | 2015-09-23 |
Family
ID=46178741
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201280017094.2A Active CN103476344B (zh) | 2011-04-08 | 2012-04-05 | 被检体信息获取装置 |
Country Status (5)
Country | Link |
---|---|
US (2) | US9655589B2 (zh) |
EP (1) | EP2693952B1 (zh) |
JP (1) | JP6039220B2 (zh) |
CN (1) | CN103476344B (zh) |
WO (1) | WO2012137983A1 (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2693952B1 (en) * | 2011-04-08 | 2019-06-12 | Canon Kabushiki Kaisha | Subject information acquisition apparatus |
JP5930611B2 (ja) * | 2011-05-26 | 2016-06-08 | キヤノン株式会社 | 被検体情報取得装置 |
JP5753633B2 (ja) * | 2013-03-05 | 2015-07-22 | 日立アロカメディカル株式会社 | 超音波診断装置、及び送受信方法 |
WO2015059878A1 (ja) * | 2013-10-22 | 2015-04-30 | 国立大学法人東京農工大学 | 粘性係数算出装置、押込試験装置、引張試験装置、粘性係数算出方法およびプログラム |
JP5862914B1 (ja) * | 2014-07-10 | 2016-02-16 | 高周波粘弾性株式会社 | 粘弾性特性測定装置及び粘弾性特性測定方法 |
CN104297104A (zh) * | 2014-09-29 | 2015-01-21 | 安徽农业大学 | 一种高精度变温式液体黏滞系数测量仪 |
WO2018031725A1 (en) * | 2016-08-11 | 2018-02-15 | Mayo Foundation For Medical Education And Research | Loss-angle-based determination of a medium viscoelastic parameter in sub-hertz frequency range with the use of local creep response |
KR102627723B1 (ko) * | 2017-09-08 | 2024-01-23 | 삼성메디슨 주식회사 | 초음파 영상장치 및 그 제어방법 |
US11154277B2 (en) * | 2017-10-31 | 2021-10-26 | Siemens Medical Solutions Usa, Inc. | Tissue viscoelastic estimation from shear velocity in ultrasound medical imaging |
CN109470565B (zh) * | 2018-11-29 | 2021-04-23 | 交通运输部公路科学研究所 | 一种沥青混合料动态泊松比的测试方法 |
WO2023186550A1 (en) * | 2022-03-31 | 2023-10-05 | Koninklijke Philips N.V. | Generating an ultrasound device interaction indicator |
EP4252667A1 (en) * | 2022-03-31 | 2023-10-04 | Koninklijke Philips N.V. | System and method for generating an ultrasound device interaction indicator |
CN116026682B (zh) * | 2023-03-30 | 2023-07-04 | 浙江大学 | 基于qme的快速弹性成像计算方法 |
CN116539732B (zh) * | 2023-06-27 | 2023-10-03 | 北京航空航天大学 | 一种检测体外冲击波加载过程及衰减规律的装置与方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1762180A1 (en) * | 2004-06-09 | 2007-03-14 | Hitachi Medical Corporation | Elastic image display method and ultrasonographic device |
CN101431943A (zh) * | 2006-03-29 | 2009-05-13 | 超声成像公司 | 对粘弹性介质进行成像的方法与装置 |
EP2272435A1 (en) * | 2008-04-25 | 2011-01-12 | Hitachi Medical Corporation | Reference deformable body, ultrasonic diagnostic apparatus, and ultrasonic diagnostic method |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU660034B2 (en) | 1991-05-10 | 1995-06-08 | Board Of Regents, The University Of Texas System | Method for elastographic measurement and imaging |
JP2005066041A (ja) * | 2003-08-25 | 2005-03-17 | Hitachi Medical Corp | 超音波探触子及び超音波診断装置 |
JP4189840B2 (ja) | 2003-10-20 | 2008-12-03 | 独立行政法人産業技術総合研究所 | 超音波を利用した軟組織の粘弾性推定装置およびプログラム |
JP4854212B2 (ja) | 2005-03-31 | 2012-01-18 | 日立アロカメディカル株式会社 | 超音波診断装置 |
JP4793726B2 (ja) | 2006-01-24 | 2011-10-12 | 独立行政法人産業技術総合研究所 | 超音波診断装置 |
US20080267499A1 (en) * | 2007-04-30 | 2008-10-30 | General Electric Company | Method and system for automatic detection of objects in an image |
JP5252985B2 (ja) | 2008-05-02 | 2013-07-31 | 花王株式会社 | 皮膚内部の弾性計測方法 |
US8394026B2 (en) | 2008-11-03 | 2013-03-12 | University Of British Columbia | Method and apparatus for determining viscoelastic parameters in tissue |
JP5422335B2 (ja) | 2009-10-15 | 2014-02-19 | 矢崎総業株式会社 | ヒューズ |
EP2693952B1 (en) * | 2011-04-08 | 2019-06-12 | Canon Kabushiki Kaisha | Subject information acquisition apparatus |
JP5996908B2 (ja) * | 2012-04-02 | 2016-09-21 | 富士フイルム株式会社 | 超音波診断装置及び超音波診断装置の表示方法 |
-
2012
- 2012-04-05 EP EP12724414.3A patent/EP2693952B1/en active Active
- 2012-04-05 WO PCT/JP2012/059952 patent/WO2012137983A1/en active Application Filing
- 2012-04-05 CN CN201280017094.2A patent/CN103476344B/zh active Active
- 2012-04-05 US US14/110,125 patent/US9655589B2/en active Active
- 2012-04-06 JP JP2012087762A patent/JP6039220B2/ja active Active
-
2017
- 2017-04-18 US US15/490,776 patent/US10639011B2/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1762180A1 (en) * | 2004-06-09 | 2007-03-14 | Hitachi Medical Corporation | Elastic image display method and ultrasonographic device |
CN101431943A (zh) * | 2006-03-29 | 2009-05-13 | 超声成像公司 | 对粘弹性介质进行成像的方法与装置 |
EP2272435A1 (en) * | 2008-04-25 | 2011-01-12 | Hitachi Medical Corporation | Reference deformable body, ultrasonic diagnostic apparatus, and ultrasonic diagnostic method |
Also Published As
Publication number | Publication date |
---|---|
JP6039220B2 (ja) | 2016-12-07 |
US20140221833A1 (en) | 2014-08-07 |
EP2693952A1 (en) | 2014-02-12 |
WO2012137983A1 (en) | 2012-10-11 |
US9655589B2 (en) | 2017-05-23 |
EP2693952B1 (en) | 2019-06-12 |
CN103476344A (zh) | 2013-12-25 |
US20170215847A1 (en) | 2017-08-03 |
US10639011B2 (en) | 2020-05-05 |
JP2012223565A (ja) | 2012-11-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103476344B (zh) | 被检体信息获取装置 | |
US20230263504A1 (en) | Shear wave elastrography method and apparatus for imaging an anisotropic medium | |
KR101398948B1 (ko) | 진폭―위상 변조된 초음파를 이용한 점탄성 측정 | |
US9986973B2 (en) | Method for shear wave ultrasound vibrometry with interleaved push and detection pulses | |
US5474070A (en) | Method and apparatus for elastographic measurement and imaging | |
US7753847B2 (en) | Ultrasound vibrometry | |
US7736315B2 (en) | Method and apparatus providing improved ultrasonic strain measurements of soft tissue | |
US7601122B2 (en) | Ultrasonic elastography with angular compounding | |
US7785259B2 (en) | Detection of motion in vibro-acoustography | |
CN109077754B (zh) | 一种测量组织力学特性参数的方法及设备 | |
EP2881041A1 (en) | Apparatus and method for ultrasonic diagnosis | |
US10136876B2 (en) | System and method for local estimation of nonlinear tissue elasticity with acoustic radiation force | |
JP5411699B2 (ja) | 超音波撮像装置 | |
KR20140129246A (ko) | 다중의 어퍼처 초음파를 사용한 물질 강성의 결정 | |
JP5075830B2 (ja) | 超音波診断装置 | |
JP5930611B2 (ja) | 被検体情報取得装置 | |
US11540809B2 (en) | Ultrasonic diagnostic apparatus and method for propagation speed analysis of shear wave and elastic modulus measurement of a tissue | |
JP6492230B2 (ja) | スペクトル解析装置、スペクトル解析方法及び超音波撮像装置 | |
WO2021212494A1 (zh) | 粘弹性测量方法和超声测量*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |