CN101865931B - 声速测定装置以及声速测定方法 - Google Patents

声速测定装置以及声速测定方法 Download PDF

Info

Publication number
CN101865931B
CN101865931B CN201010163281.4A CN201010163281A CN101865931B CN 101865931 B CN101865931 B CN 101865931B CN 201010163281 A CN201010163281 A CN 201010163281A CN 101865931 B CN101865931 B CN 101865931B
Authority
CN
China
Prior art keywords
sound
ripple
travel
time
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.)
Expired - Fee Related
Application number
CN201010163281.4A
Other languages
English (en)
Other versions
CN101865931A (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.)
Furuno Electric Co Ltd
Original Assignee
Furuno Electric 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 Furuno Electric Co Ltd filed Critical Furuno Electric Co Ltd
Publication of CN101865931A publication Critical patent/CN101865931A/zh
Application granted granted Critical
Publication of CN101865931B publication Critical patent/CN101865931B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0875—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of bone

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Rheumatology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Orthopedic Medicine & Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明提供一种可以正确且稳定地求出声速的声速测定装置以及声速测定方法。骨强度诊断装置具备波发送专用振子、多个振子、假定传播时间计算部、波形累计部、声速导出部。波发送专用振子向皮质骨的表面发送超声波。振子接收来自皮质骨的超声波，并输出与接收的超声波对应的波形信号。假定传播时间计算部根据该皮质骨内的声速的假定值和该皮质骨的表面形状计算波发送专用振子发送的超声波在皮质骨的表面附近传播之后到达各振子的传播时间。波形累计部求出将多个振子输出的各波形信号移动传播时间之后累计所得的波形的包络线的振幅。声速导出部根据包络线的振幅求出皮质骨内的声速。

Description

声速测定装置以及声速测定方法
技术领域
本发明涉及声速测定装置以及声速测定方法。详细而言，涉及在采用超声波从外部测定被测体内部的声速的声速测定装置以及声速测定方法中，用于可以进行正确且稳定的声速测定的技术。
背景技术
专利文献1公开了一种采用超声波来进行骨评价的超声波诊断装置。该超声波诊断装置放射超声波脉冲，并接收在皮质骨内传播之后从该皮质骨放出的超声波。然后，超声波诊断装置测定从波发送到波接收的时间，并根据已知的传播路径求出皮质骨内的声速。
物体中的声速根据该物体的弹性的性质而发生变化，所以通过测定骨中的声速能够作为骨强度(骨的健全性)的指标。专利文献1的超声波诊断装置将如上述那样所求出的声速作为骨评价值。
专利文献1：特表2003-517328号公报
如专利文献1那样，在对超声波脉冲的从波发送到波接收的间隔进行定时器测量的结构的情况下，不得不适当地从振子输出的波形信号中检测出接收脉冲的峰值。这里，若在理想的条件(例如，置于水中的对象物的测定等)下，则检测接收脉冲是比较容易的。但是，在以生物体为对象的实际的测定中，因为在超声波的传播路径中混在有脂肪成分和水分所以产生多余回波、或来自皮质骨背面的反射等不需要的波形重叠等，产生各种噪声。
这点，对于专利文献1的结构，一旦在波形信号中包含噪声，则难以从波形信号中检测必要的波接收脉冲，不能稳定地测量超声波脉冲的从波发送到波接收的时间。因此，即使在实际现场要检查骨强度，也难以获得正确的结果，留下了改良的余地。
发明内容
本发明鉴于以上的情况而做，其主要的目的在于提供一种能够正确且稳定地求出声速的声速测定装置以及声速测定方法。
本发明要解决的课题如以上所述，下面说明用于解决该课题的手段及其效果。
根据本发明的第1观点，提供一种以下结构的声速测定装置。即，该声速测定装置具备波发送部、多个波接收部、假定传播时间计算部、妥当性指标值计算部和声速导出部。所述波发送部向被测体的表面发送超声波。所述波接收部接收来自所述被测体的超声波，并输出与该接收到的超声波相应的波形信号。所述假定传播时间计算部根据该被测体内的声速的假定值和该被测体的表面形状，计算所述波发送部发送的超声波在所述被测体的表面或该被测体的内部传播之后到达各所述波接收部的传播时间。所述妥当性指标值计算部根据所述多个波接收部中的至少两个波接收部输出的波形信号，求出成为所述传播时间的妥当性的指标的妥当性指标值。所述声速导出部根据所述妥当性指标值求出所述被测体内的声速。
即，通过判断传播时间的妥当性，能够判断在计算该传播时间时所假定的声速的值是否正确，所以由此能够求出被测体内的声速。这里，在判断传播时间的妥当性时，不需要严密地检测接收波形的峰值，所以即使在接收波形含有噪声也可以适当地求出声速。此外，根据来自多个波接收部的波形信号来判断传播时间的妥当性，所以可以减弱噪声的影响。
在所述声速测定装置中，优选所述妥当性指标值计算部将所述多个波接收部中的至少两个波接收部输出的波形信号分别移动与所述传播时间对应的时间，求出将该移动后的波形信号彼此累计所得的累计波形，根据该累计波形求出所述妥当性指标值。
即，累计多个波接收部输出的波形信号彼此时，在波形的峰值位置一致的情况下，最加强振幅。通过利用该性质，可以仅强调希望的峰值，并减少其他噪声的影响。在上述结构中，在假定传播时间计算部所假定的声速正确时，累计后的波形信号的峰值位置一致，所以加强振幅。因此，即使在波接收部输出的波形信号中含有噪声的情况下，通过调查累计波形，也能够判断声速的假定值是否正确。
在所述声速测定装置中，优选所述妥当性指标值计算部将所述多个波接收部中的至少两个波接收部输出的波形信号分别移动与所述传播时间对应的时间，求出将该移动后的波形信号彼此相乘所得的乘法波形，根据该乘法波形求出所述妥当性指标值。
即，在将多个波接收部输出的波形信号彼此相乘时，在波形的峰值位置一致时，最加强振幅。通过利用该性质，可以仅强调希望的峰值，并减少其他噪声的影响。在上述结构中，在假定传播时间计算部所假定的声速正确的情况下，相乘所得的波形信号的峰值位置一致，所以加强振幅。因此，即使在波接收部输出的波形信号中含有噪声的情况下，通过调查乘法波形，也能判断声速的假定值是否正确。
所述声速测定装置优选以下那样的结构。即，所述假定传播时间计算部采用多个值计算所述传播时间作为所述声速的假定值。所述妥当性指标值计算部根据所述多个声速的假定值，求出各个情况的妥当性指标值。所述声速导出部根据所述多个声速的假定值各个情况的妥当性指标值，求出所述被测体内的声速。
由此，试行多个声速的假定值，在传播时间最妥当时，可以采用计算该传播时间时声速的假定值作为被测体内的声速的测定值。
在所述声速测定装置中，优选所述假定传播时间计算部根据预先测定的被测体表面的形状，计算所述传播时间。
由此，能够考虑被测体表面的实际形状从而计算假定传播时间，所以即使在例如被测体表面是曲面状的情况下，也能够精度良好的求出声速。
所述声速测定装置优选以下那样的结构。即，该声速测定装置具备形状检测部。所述多个波接收部中至少一部分构成为能够向被测体的表面发送超声波。所述形状检测部根据从超声波被发送开始到在所述被测体的表面反射之后到达各所述波接收部为止的时间，检测所述被测体表面的形状。所述假定传播时间计算部根据所述形状检测部检测的所述被测体表面的形状，计算所述传播时间。
由此，利用用于测定声速的结构的一部分，能够检测被测体表面的形状，所以可以简化装置从而削减费用。
所述声速测定装置可以以下那样地构成。即，所述多个波接收部构成为以大致等间隔排列为1列的线性阵列，并且所述多个波接收部中的至少任意一个波接收部作为所述波发送部发挥功能。
由此，能够省略波发送专用的结构，所以可以简化装置。此外，从接近的多个波接收部错开定时，通过调整使该定时错开的时间，能够以任意的角度形成超声波束。由此，能够对被测体表面以最合适的角度照射超声波。此外，从波接收部向被测体发送超声波，根据反射波到达各波接收部的时间，检测该被测体的表面形状的结构的情况下，能够检测该表面形状的仅是波接收部的正下方的位置。因此，在分别构成波发送部和波接收部的情况下，波发送部正下方的表面形状不能检测。对于这点，通过波接收部兼具波发送部，从而能够检测作为波发送部的波接收部的正下方的表面形状。由此，可以正确地计算传播时间。
根据本发明的第2观点，是一种利用具备了根据所接收到的超声波来输出波形信号的多个波接收部的声速测定装置的声速测定方法，提供一种以下那样的声速测定方法。即，该声速测定方法包括波发送工序、波接收工序、假定传播时间计算工序、妥当性指标值计算工序。在所述波发送工序中，向被测体的表面发送超声波。在所述波接收工序中，由各所述波接收部接收来自该被测体的超声波。在所述假定传播时间计算工序中，根据该被测体内的声速的假定值和该被测体的表面形状，计算在所述波发送工序发送出的超声波在所述被测体的表面或该被测体的内部传播之后到达各所述波接收部为止的传播时间。在所述妥当性指标值计算工序中，根据所述多个波接收部中的至少两个波接收部输出的波形信号，求出成为所述传播时间的妥当性的指标的妥当性指标值。而且，使所述声速的假定值变化，并且反复进行所述假定传播时间计算工序和所述妥当性指标值计算工序，由此求出多个声速的假定值各自情况的妥当性指标值，根据所述各自情况的妥当性指标值，求出所述被测体内的声速。
即，通过判断传播时间的妥当性，能够判断在计算该传播时间时所假定的声速的值是否正确，所以由此能够求出被测体内的声速。这里，在判断传播时间的妥当性时，不需要严密地检测接收波形的峰值，所以即使在接收波形中含有噪声也可以适当地求出声速。此外，因为根据来自多个波接收部的波形信号判断传播时间的妥当性，所以可以减弱噪声的影响。而且，试行多个声速的假定值，在传播时间最妥当时，可以采用计算该传播时间时的声速的假定值作为被测体内的声速的测定值。
在所述声速测定方法中，优选采用以下的方法。即，在所述妥当性指标值计算工序中，将所述多个波接收部中的至少两个波接收部输出的波形信号分别移动与所述传播时间对应的时间，求出将该移动后的波形信号彼此进行累计所得累计波形，根据该累计波形来求出所述妥当性指标值。
即，在累计多个波接收部输出的波形信号彼此时，波形的峰值位置一致时最加强振幅。通过利用该性质，可以仅强调希望的峰值并减少其他噪声的影响。在上述结构中，在假定传播时间计算工序所假定的声速正确时，被累计的波形信号的峰值位置一致，所以加强振幅。因此，即使在波接收部输出的波形信号中含有噪声的情况下，通过观察累计波形的振幅，也能够判断声速的假定值是否正确。
在所述声速测定方法中可以采用以下的方法。即，在所述妥当性指标值计算工序中，将所述多个波接收部中的至少两个波接收部输出的波形信号分别移动与所述传播时间对应的时间，求出将该移动后的波形信号彼此相乘所得的乘法波形，根据该乘法波形来求出所述妥当性指标值。
即，在将多个波接收部输出的波形信号彼此相乘时，在波形的峰值位置一致时，最加强振幅。通过利用该性质，可以仅强调希望的峰值并减少其他噪声的影响。在上述结构中，在假定传播时间计算工序所假定的声速正确的情况下，被相乘的波形信号的峰值位置一致，所以加强振幅。因此，即使在波接收部输出的波形信号中含有噪声的情况下，通过观察乘法波形的振幅，也可以判断声速的假定值是否正确。
在所述声速测定方法中优选采用以下的方法。即，该声速测定方法包括检测所述被测体的表面形状的形状检测工序。而且，在所述假定传播时间计算工序中，根据所述表面形状来计算所述传播时间。
由此，可以考考被测体表面的实际形状来计算假定传播时间，所以即使在例如被测体表面是曲面状的情况下，也可以精度良好地求出声速。
附图说明
图1是本发明的一实施方式涉及的骨强度诊断装置的示意性剖面图和功能模块图。
图2是说明超声波的多个传播路径的概念图。
图3是例示各振子输出的波形信号的曲线图。
图4是本实施方式的声速测定方法的流程图。
图5(a)是表示由阵列振子发送平面波的样子的图。(b)是表示阵列振子发送的平面波在皮质骨的表面或背面反射的样子的图。
图6(a)是放大了接收表面反射波的振子组的附近的示意图。(b)是说明到达构成振子组的两个振子的表面反射波的传播路径的差的示意图。
图7是说明表面折射波的传播路径的计算方法的示意图。
图8是表示假定传播时间计算部求出的t-x曲线的图。
图9是表示对各振子的波形信号实施了窗口函数的样子的图。
图10是表示移动了各振子的波形信号的样子的图。
图11是表示波形累计部求出的累计波形的样子的图。
图12是表示使声速假定值变化从而描画t-x曲线的样子的图。
图13是表示边使声速假定值变化边求出包络线的振幅从而描画的结果的例的图。
图14是表示由本实施方式的骨强度诊断装置测定的牛胫骨皮质骨的声速的结果的图。
图15是变形例涉及的超声波收发器的示意性剖面图。
图16是表示对表面传播波来自表面的反射波和来自背面的反射波重合的样子的图。
图17是在图13的曲线图中追加描画了累计波形成为最大振幅的时间的曲线图。
符号说明
1骨强度诊断装置(声速测定装置)
2超声波收发器
3装置主体
10皮质骨
21波发送专用振子(波发送部)
22阵列振子
24振子(波接收部)
35运算部
40形状检测部
50声速检测部
51假定传播时间计算部
52波形累计部(妥当性指标值计算部)
53声速导出部
具体实施方式
下面，参照附图说明发明的实施方式。图1是作为本发明的一实施方式涉及的声速测定装置的骨强度诊断装置1的示意性剖面图和功能模块图。
作为声速测定装置的骨强度诊断装置1例如是对胫骨等长管状骨的皮质骨诊断骨强度的装置(其中，诊断对象不限于此)。若具体说明，则骨一般由皮质骨10和存在于皮质骨10的内侧的网状松质骨12构成。此外，皮质骨10的周围被肌肉和脂肪等软组织11覆盖。本实施方式的骨强度诊断装置1构成为从软组织11的外侧对皮质骨10放射超声波，测定该皮质骨10中的声速(骨声速)。
在图1的左侧示出用与骨的长边方向正交的平面切断人体的胫部分的剖面。如图1所示，皮质骨10表面的轮廓形状成为向径方向(与骨的长边方向正交的方向)鼓起的平稳的曲线状。因此，在以下的说明中，有将长管状的骨比作圆柱、将在该剖面内超声波沿皮质骨10表面那样传播的方向称为圆周方向的情况。另一方面，虽然图示省略，但是在用与骨的长边方向平行的平面切断的剖面中，皮质骨10的表面轮廓成为大致直线状。
可是，人的骨具有各向异性，所以长轴方向的声速与圆周方向的声速不同。因此，分别测定圆周方向和长轴方向的声速并调查各向异性构造在诊断骨强度上是重要的。因此，本实施方式的骨强度诊断装置1构成为，如测定骨的圆周方向的声速的情况那样，即使在测定对象的剖面轮廓是曲线状的情况下，也能够正确地测定声速。
以下，具体地说明骨强度诊断装置1的结构。如图1所示，骨强度诊断装置1由超声波收发器2和装置主体3构成。
超声波收发器2进行超声波的波发送和波接收。该超声波收发器2具备：与测定部位的软组织11的表面抵接的抵接面2a、波发送专用振子21、阵列振子22、隔音件23。阵列振子22由排成1列而排列的多个振子24构成。此外，波发送专用振子21和隔音件23配置为沿阵列振子22的排列方向而排列。另外，作为在本实施方式使用的振子，采用若给予电信号则其表面振动从而发射超声波并且若在其表面接收超声波则生成并输出电信号的部件。
波发送专用振子(波发送部)21构成为设置为其表面相对于抵接面2a倾斜，能够从抵接面2a沿倾斜方向发送超声波。作为波发送专用振子21，采用放射的超声波的指向性弱(超声波的角度范围广)的部件。
阵列振子22具有多个振子(波接收部)24。另外，在以下的说明中，在需要区别多个振子24的情况下，从靠近波发送专用振子21的一侧顺次在符号的末尾附加小写的字母，表记为振子24a、振子24b、振子24c、……。该振子24配置为按照与抵接面2a平行的方式以等间隔排列为1列。此外，各振子24构成为可以发送和接收超声波。
隔音件23形成为板状，配置在波发送专用振子21和阵列振子22之间。该隔音件23用于防止以下的情况，即：从波发送专用振子21发送的超声波在超声波收发器2内传播之后直接到达阵列振子22。另外，作为隔音件23的材料，例如，可以采用软木、合成橡胶、多孔材料(例如发泡树脂)等具有吸音作用的材料。
在采用该超声波收发器2实际进行超声波的收发时，在测定部位的皮肤表面(即，软组织11的外侧表面)涂敷超声波胶(jelly)，使该皮肤表面抵接抵接面2a。然后，通过波发送专用振子21或阵列振子22，发送超声波。由此，超声波经由软组织11到达测定对象即皮质骨10。然后，从皮质骨10返回来的超声波被阵列振子22接收。另外，所述超声波胶防止在软组织11和抵接面2a之间产生间隙，并且调整抵接面2a和软组织11之间的声阻抗，抑制从波发送专用振子21或阵列振子22发送的超声波在软组织11的表面反射。
下面，说明装置主体3。装置主体3构成为通过电缆与超声波收发器2连接，能够在与该超声波收发器2之间收发信号。具体而言，该装置主体3具备超声波控制部30、发送电路31、发送切换部32、多个接收电路33、收发分离部34和运算部35。
发送电路31构成为生成用于使波发送专用振子21或阵列振子22振动从而产生超声波的电脉冲信号，并且将该电脉冲信号发送到发送切换部32。电脉冲振动的中心频率是例如1～10MHz左右。另外，也可以代替电脉冲信号而采用例如线性调频(chirp)信号。
另外，在由阵列振子22产生超声波时，发送电路31构成为能够分别对多个振子24生成任意定时的电脉冲信号。此外，超声波控制部30构成为连接于发送电路31，并将用于从多个振子24发送超声波的控制信号发送给发送电路31。由此，成为能够按照从多个振子24一齐或在个别的定时发送超声波的方式进行控制。
发送切换部32构成为进行切换将从发送电路31发送的电脉冲信号发送到波发送专用振子21或阵列振子22的任一个。即，通过发送切换部32，选择进行超声波的发送的振子。
多个接收电路33分别与构成阵列振子22的多个振子24对应地连接。各接收电路33构成为接收一个振子24通过接收超声波而输出的电信号，并生成对该电信号实施了放大处理、滤波器处理、数字变换处理等的数字信号，之后发送给运算部35。另外，从阵列振子22直接输出的信号是模拟的波形信号，向运算部35发送的信号是信号处理后的数字的波形信号，但在以下的说明中有不区别这些而仅称为「波形信号」的情况。
收发分离部34连接于阵列振子22和所述发送电路31以及所述接收电路33之间。该收发分离部34用于防止从发送电路31向阵列振子22发送的电信号(电脉冲信号)直接流入接收电路33，并且防止从阵列振子22向接收电路33发送的电信号流入发送电路31侧。
下面，说明由波发送专用振子21发送超声波时的样子。在由波发送专用振子21发送超声波时，通过发送切换部32，决定波发送专用振子21作为进行超声波的波发送的振子。然后，若来自发送电路31的脉冲信号发送到波发送专用振子21，则该波发送专用振子21对皮质骨10沿倾斜方向发送被调制为脉冲状的超声波。
从波发送专用振子21发送的超声波经由多个传播路径被阵列振子22接收。若被阵列振子22接收了超声波，则来自各振子24的波形信号发送到运算部35。
对于传播从波发送专用振子21发送的超声波的多个路径，参照图2来进行说明。图2是示意性表示到达某振子24的超声波的多个传播路径的图。另外，虽然在图2中仅示出到达一个振子24的超声波的例，但是实际上，对全部振子24的每一个，超声波可以由多个路径到达。
如图2所示，直接波是从波发送专用振子21发送的超声波沿软组织11的表面传播，从而直接到达振子24的波。来自表面的反射波是从波发送专用振子21发送的超声波在软组织11内传播并在皮质骨10的表面(软组织11和皮质骨10的边界面)反射从而到达振子24的波。此外，来自背面的反射波是在软组织11中传播的超声波射入到皮质骨10内后在该皮质骨10内传播并在皮质骨10的背面(皮质骨10和松质骨12的边界面)反射之后再次放射到软组织11中从而到达振子24的波。
此外还有如下的情况，即：从波发送专用振子21发送的超声波在软组织11中传播，到达皮质骨10，在该皮质骨10表面附近传播，并再次放射到软组织11中从而到达振子24。在本说明书中，将经过这样的传播路径而接收的超声波称为表面传播波。在表面传播波中有两个传播路径，即泄露表面波和表面折射波。
若超声波以临界角到达皮质骨10的表面，则在皮质骨10表面产生表面波。该表面波向软组织11侧沿规定方向(射出角成为临界角的方向)放射泄露波，并且沿皮质骨10的表面传播。将此时通过振子24接收的泄露波称为泄露表面波。另一方面，若超声波以比临界角小的角度到达皮质骨10的表面，则在皮质骨10表面折射。此时，在在超声波以接近临界角的射入角射入时，在皮质骨10内在该皮质骨10的表面附近传播之后，在皮质骨10表面向阵列振子22侧折射从而放射到软组织11中。在本说明书中将此时接收的超声波称为表面折射波。表面折射波仅在皮质骨10的剖面轮廓形状是曲线状的情况下产生。
另外，根据骨的形状、波接收的振子24的位置以及从波发送专用振子21发送的超声波的角度等条件，会有上述多种超声波中一部分未产生的情况、以及即使产生也未被振子24接收的情况。但是，在本实施方式中，如前所述，因为采用了指向性广的波发送专用振子21，所以能够使超声波以临界角或接近临界角的角度到达皮质骨10。由此，能够确实地产生泄露表面波或表面折射波，并且能够使泄露表面波或表面折射波到达多个振子24中的至少任意一个。
下面，参照图3说明各个振子24生成的信号波形。图3是表示波发送专用振子21发送了超声波之后各振子24接收超声波并输出的波形信号的波形的曲线图。对于各振子24的波形，横轴t表示时间，各波形信号的纵轴表示该波形信号的振幅。此外，曲线图的x轴表示从距离波发送专用振子21最近的振子即振子24a到各振子24b、24c、……的距离。
另外，伴随波接收的振子24远离波发送专用振子21，波形信号的振幅变小。图3所示的波形的曲线是对各波形信号实施适当的增益从而调整为各波形信号的振幅大致相同后的曲线。
如图3所示，在各振子24的波形信号中含有多个峰值。在图3中，用虚线连接各振子24的波形信号中所含的峰值的位置，该峰值表示基于直接波、来自表面的反射波、来自背面的反射波、以及表面传播波的任一个的峰值。另外，在该曲线图中，因为难以区别泄露表面波和表面折射波，所以将二者统一表示为表面传播波。
由图3的曲线图可知，若有直接波比表面传播波先到达振子24的情况，则也有后到达的情况。此外，有来自背面的反射波和来自表面的反射波几乎与表面传播波同时到达振子24的情况，在该情况下，峰值重合，难以区别。这样，即使波形信号中含有峰值，也难以判别其是从哪个路径传播来的超声波的峰值。
此外，因为软组织11是脂肪和水分的混合体，所以在实际的测定中，检测出多余回波。在该情况下，在图3的波形中进一步加入噪声，所以稳定地检测出峰值变得越发困难。
这里，在以往的声速测定装置中，通过测定从发送超声波到接收表面传播波为止的时间间隔，从而求出骨声速。因此，即使在以往的声速测定装置中也需要检测表面传播波的峰值，但是如上所述，稳定地检测出表面传播波的峰值是困难的。因此，在以往的方法中，未能稳定地测定骨声速。
因此，在本实施方式中，通过对多个振子24输出的波形信号的波形进行累计，从而实现了抗噪声强的声速测定方法。
即，在多个波形的峰值位置一致的情况下，累计这些波形后，该峰值位置的振幅变大。通过利用该性质，能够仅强调希望的峰值，使噪声变得不显著。但是，因为波到达各个振子24的定时各异，所以要加强基于特定波的峰值，需要在按照使该峰值一致的方式使各波形的相位移动的基础上累计波形彼此。
以下，对本实施方式的声速测定方法进行具体地说明。图4表示本实施方式涉及的声速测定方法的流程图。本实施方式的声速测定方法含有形状检测工序、波发送工序、波接收工序、假定传播时间计算工序、波形累计工序和声速导出工序。
在形状检测工序中，测定皮质骨10的表面形状(S101)。表面形状也可以用例如X射线等来测定，但在本实施方式中，通过骨强度诊断装置1来测定皮质骨10的表面形状。另外，后面叙述该表面形状测定的详细情况。
在波发送工序中，通过波发送专用振子21对皮质骨10发送超声波(S102)。
在波接收工序中，由各振子24接收从波发送专用振子21发送的超声波(S103)，获得图3所示的波形信号。如图3所示，各振子24接收的表面传播波的峰值伴随该振子24的位置远离波发送专用振子21而较晚地被检测。该延迟量相当于从波发送专用振子21发送超声波开始到各振子24接收表面传播波为止的时间(传播时间)。
因此，通过使各波形信号提前表面传播波的传播时间(使各波形信号向图3的左侧偏移)，能够使表面传播波的峰值的相位在各振子24之间一致。但是，因为骨声速是未知的，所以不能预知传播时间。
因此在本实施方式中，在假定传播时间计算工序中，假定骨声速，并根据该假定的骨声速(声速假定值)、和在形状检测工序检测出的皮质骨10的表面形状，求出对各振子24的传播时间的假定值(假定传播时间)(S104)。
接着，在波形累计工序(妥当性指标值计算工序)中，在使各波形信号偏移与所述假定传播时间对应的时间的基础上，求出对波形信号进行了累计的累计波形(S105)。此时，作为表示假定传播时间的妥当性的妥当性指标值，求出所述累计波形的包络线(envelope)的振幅(详细情况后述)。
然后，在规定的范围内，使声速假定值接连变化，从而求出对各个声速假定值的累计波形(S104～S106的循环)。于是，在声速假定值与实际的骨声速一致时，各波形的表面传播波的峰值的相位一致，可获得峰值彼此最加强的累计波形。
因此，在声速导出工序中，比较针对多个声速假定值而求出的累计波形的包络线的振幅彼此，检测出累计波形的振幅变得最大的时刻，将该时刻的声速假定值作为骨声速的测定值(S107)。
根据以上的方法，峰值的相位不一致的波(即，表面传播波的峰值以外的噪声)不会给予累计波形的振幅大的影响，所以仅强调表面传播波(若是表面传播波则使焦点匹配)，就能够精度良好地求出骨声速。
下面，说明用于通过本实施方式的骨强度诊断装置1来实现上述声速测定方法的结构。
图1所示的骨强度诊断装置1的装置主体3具备的运算部35由CPU、RAM、ROM等硬件、和所述ROM中所存储的程序等软件构成。而且，该运算部35构成为通过所述硬件和软件协作，作为形状检测部40、声速检测部50等发挥功能。
首先，说明形状检测部40。形状检测部40进行声速的导出所需要的骨表面形状的检测。
即，在根据假定的骨声速(声速假定值)求出假定传播时间时，需要与表面传播波的传播路径相关的信息。这里，为了求出传播路径，不得不利用某些方法取得皮质骨10表面的形状(皮质骨10的剖面轮廓形状)。因此在本实施方式中，在由声速检测部50求出骨声速之前，进行形状检测部40进行的皮质骨10的表面形状的检测。
在由形状检测部40进行骨表面形状的检测时，首先，从阵列振子22进行超声波的发送。参照图5说明由阵列振子22发送超声波的样子。图5(a)是表示由阵列振子22发送超声波的样子的图，图5(b)是表示阵列振子22发送的超声波在皮质骨10的表面或背面反射的样子的图。
在由阵列振子22发送超声波时，通过发送切换部32，决定阵列振子22作为进行超声波的波发送的振子。而且，来自发送电路31的脉冲信号送给阵列振子22，构成该阵列振子22的多个振子24同时对骨发送同相位的超声波。
通过多个振子24同时进行波发送，能够产生图5(a)所示的平面波。该平面波是与抵接面2a平行的波，沿与抵接面2a正交的方向在软组织11中前进。平面波在皮质骨10的表面和背面如图5(b)那样反射，被振子24接收。
各振子24接收超声波后，向运算部35发送与各振子24接收的超声波对应的波形信号。运算部35构成为通过作为形状检测部40发挥功能，从而检测振子24接收反射波时的角度和时间，并基于此来求出皮质骨10的表面形状。
具体而言，形状检测部40由到来方向检测部41、到达时间检测部42、表面反射点检测部43、骨表面线检测部44构成。
首先，说明到来方向检测部41。到来方向检测部41将多个振子24中相邻的2个振子作为1组来决定振子组25，并检测到达各振子组25的超声波的到来方向。另外，在关于以下的形状检测部40的说明中，有将从阵列振子22发送的平面波在皮质骨10的表面反射之后被接收的超声波称为表面反射波，相同地将在皮质骨10的背面反射之后被接收的超声波称为背面反射波的情况。此外，在需要区别各振子组25的情况下，从距波发送专用振子21近的一侧开始依次在符号的末尾附加大写字母，表记为振子组25A、振子组25B、……。
以下，参照图6具体地进行说明。图6(a)是扩大了接收表面反射波的振子组25A的附近的示意图，图6(b)是说明到达构成振子组的2个振子24a、24b的表面反射波的传播路径的差的示意图。在某振子组25中，表面反射波对相邻的2个振子24的到来方向近似。例如在图6中，可以看做表面反射波分别以到来角度θa到达构成振子组25A的振子24a和振子24b。这里，为了求出到来角度θa，进行以下那样的运算。
首先，到来方向检测部41测定构成振子组25A的2个振子24a、24b检测表面反射波的峰值的时间差Δt。另外，如前所述，从阵列振子22发送平面波时产生表面反射波和背面反射波，但是因为必定比背面反射波先接收表面反射波，所以能够适当地检测表面反射波的峰值。
接着，基于该时间差Δt，求出表面反射波对该振子组25A的到来角度θa。如图6(b)所示，设振子24a和振子24b间的间隔为W时，表面反射波传播比振子24b长Wsinθa的距离到达振子24a。这里，设软组织中的声速为SOSsoft时，
SOSsoftΔt＝Wsinθa
因此，到来角度θa能够由下式求出。
θa＝arc sin(SOSsoftΔt/W)
到来方向检测部41对其他振子组25也同样地求出到来角度。另外，在本实施方式中作为软组织11中的声速SOSsoft采用了由经验获得的值，但也可以利用实测值。
下面，说明到达时间检测部42。到达时间检测部42求出从由阵列振子22发送了超声波之后开始到表面反射波到达振子组25为止的到达时间Ta。在本实施方式中，将从由阵列振子22发送了超声波之后开始到表面反射波分别到达构成振子组25的2个振子24为止的时间的平均值作为到达时间Ta。另外，不限于平均值，也可以将例如表面反射波到达任一方的振子24为止的时间原样作为到达时间Ta来采用。
下面，说明表面反射点检测部43。表面反射点检测部43根据到来角度θa和到达时间Ta，检测到达各振子组25的表面反射波的反射点Ra。
这里，在图6所示的平面内，将阵列振子22排列的方向设为x轴，将与x轴正交的方向设为y轴。而且，设从振子组25A到反射点Ra的x轴方向的距离为X，设y轴方向的距离为Y。由图6可知，表面反射波的传播距离La成为
La＝Y+Y/cosθa。
另一方面，利用到达时间Ta和软组织11中的声速SOSsoft时，
La＝SOSsoft×Ta
所以表示反射点Ra的位置的距离X，Y能够由下式求出。
Y＝SOSsoft×Ta×cosθ/(1+cosθ)
X＝Y×tanθ＝SOSsoft×Ta×sinθ/(1+cosθ)
这样，根据平面波的到来角度θa和到达时间Ta，能够计算反射点Ra的位置。而且，表面反射点检测部43针对其他振子组25也可以同样地求出反射点。
骨表面线检测部44通过用直线或曲线连接表面反射点检测部43求出的多个反射点，从而检测骨表面线。反射点是皮质骨10表面上的点，所以骨表面线表示皮质骨10的表面形状。
如以上那样，通过形状检测部40能够获得皮质骨10的表面形状(骨表面线)。
下面，说明声速检测部50。声速检测部50构成为检测皮质骨10中的声速(骨声速)。
在由声速检测部50检测骨声速之前，首先，进行基于形状检测部40的骨表面线的检测(形状检测工序)。然后，由波发送专用振子21进行超声波的波发送(波发送工序)，由阵列振子22接收返回的超声波，对运算部35发送波形信号(波接收工序)。而且，根据各振子24的波形信号，进行基于声速检测部50的骨声速的导出。
具体而言，声速检测部50由假定传播时间计算部51、波形累计部52、声速导出部53构成。
假定传播时间计算部51进行与所述假定传播时间计算工序对应的运算处理。即，假定传播时间计算部51假定骨声速，并根据该假定的骨声速来求出假定传播时间。
假定传播时间计算部51为了求出假定传播时间，首先，求出表面传播波的传播路径。即，若已知皮质骨10的表面形状、骨声速、软组织11中的声速，则通过公知的斯涅耳(Snell)法则，能够唯一地求出表面传播波从波发送专用振子21到各振子24传播的路径。
这里，骨形状利用形状检测部40检测的骨表面线。作为骨声速，从基于骨声速的经验值而设定的规定范围内假定适当的值(声速假定值)来用于计算。作为软组织11中的声速，在本实施方式中利用经验值，但也可以利用预先测定的值。此外，软组织11中的声速可以通过波发送专用振子21或阵列振子22发送超声波并由阵列振子22来进行波接收从而进行测定。
另外，如上述说明的那样，在各振子24接收的表面传播波中，由传播路径的不同而存在泄露表面波和表面折射波的2种。在本实施方式中，作为计算表面折射波的传播路径来说明。
以下，参照图7进行说明。图7是说明表面折射波的传播路径的计算方法的示意图。
假定传播时间计算部51为了决定表面折射波的传播路径，首先，决定来自波发送专用振子21的超声波射入皮质骨10的内部的射入点Pin。
另外，射入点Pin处于从阵列振子22的正下方偏离的位置，所以在从阵列振子22发送平面波从而检测表面反射波的方法(形状检测部40执行的皮质骨表面形状的检测)中，不能取得该射入点Pin附近的表面形状。因此在本实施方式中，根据形状检测部40求出的骨表面线来预测射入点Pin附近的形状。在图7中，用虚线示出所预测的骨表面线。
下面，假定传播时间计算部51求出超声波对该预测的线射入的位置(射入点Pin)。超声波可以在比临界角θc小的任意的角度射入皮质骨10的内部。但是，表面折射波被振子24接收，限于超声波以与临界角θc接近的射入角射入到皮质骨10内的情况。因此在本实施方式中，与超声波以临界角θc射入皮质骨10内的情况近似地求出射入位置。临界角θc由骨声速和软组织中的声速和来决定，所以通过假定骨声速，能够唯一地求出射入点Pin。
接着，假定传播时间计算部51在骨表面线上设定任意的点Pout。从射入点Pin到点Pout的直线距离是超声波在皮质骨10内传播的距离。
而且，假定传播时间计算部51求出超声波从Pout向软组织11侧放射时的角度。超声波放射时的角度可以由斯涅耳法则唯一地求出。具体而言，设声速假定值为SOSbone、设软组织中的声速为SOSsoft时，超声波从皮质骨10向软组织11侧射出时的折射角β1和射入角β2的关系是
sinβ1/SOSsoft＝sinβ2/SOSbone。
射入角β2可以由直线PinPout的斜率、和骨表面线的形状求出。因此，通过上述式可以求出折射角β1。
这里，如图7所示，将距波发送专用振子21最近的振子24a的位置作为原点，取阵列振子22排列的方向为x轴。因此已知折射角β1和骨表面线的形状，所以确定从点Pout放射的超声波与x轴相交的点Px。通过以上说明的运算处理，能够求出从波发送专用振子21到点Px的表面折射波的传播路径。
另外，在求出射入点Pin时，上述那样地与超声波以临界角θc射入皮质骨10内近似，是为了容易计算，所计算的传播路径是近似的值。另一方面，也可以不用临界角θc近似向皮质骨10的射入角，而严密地计算传播路径。即，若已知皮质骨10的表面形状，则通过使向皮质骨10的射入角变化并两次采用斯涅耳法则，能够严密地计算从波发送专用振子21到点Px的传播路径。在该情况下，射入点Pin通过点Px变化。但是，这样的传播路径的严密的计算与用临界角θc近似射入角的计算方法相比，花费处理时间。在本实施方式中，构成为可以选择近似的计算方法和严密的计算方法的任一个来计算传播路径。
而且，假定传播时间计算部51根据所求出的传播路径，求出从由波发送专用振子21发送超声波之后开始到表面折射波到达点Px为止的传播时间。即，成为：
(到点Px的传播时间)＝(在软组织中传播的距离)×SOSsoft+(在皮质骨中传播的距离)×SOSbone。
假定传播时间计算部51边改变Pout的位置边反复上述传播时间的计算。由此，成为边改变点Px的位置，边求出到各个点Px的表面折射波的传播时间，所以能够求出图8所示的t-x曲线。另外，图8的纵轴x是距振子24a的距离，横轴t是从波发送专用振子21发送超声波之后开始到表面折射波到达所述x轴为止的传播时间。
最后，假定传播时间计算部51对各振子24求出表面折射波到达为止的时间(假定传播时间)。即，已知各振子24的x坐标(距振子24a的距离)，所以通过参照t-x曲线，能够求出各振子24的假定传播时间。
下面，说明波形累计部52。波形累计部(妥当性指标值计算部)52进行与所述波形累计工序对应的运算处理。即，波形累计部52构成为分别将各振子24输出的波形信号的波形移动(偏移)所述假定传播时间之后进行累计，求出累计波形。
以下，具体地进行说明。首先，波形累计部52在累计波形信号之前，为了消除表面折射波以外的波而对各波形信号实施适当的窗口函数(图9)。另外，为了决定窗口函数的范围，必须某种程度地预先知道骨声速或传播时间。但是，即使通过窗口函数不能全部消除表面折射波以外的波，通过波形信号彼此的累计能够减弱所残留的波的影响，所以不需要按照完全消除表面折射波以外的波的方式严密地决定窗口函数。因此，只要根据骨声速和软组织11中的声速的经验值等并考虑适当的余量来决定窗口函数就足够。
此外，波形信号的振幅伴随振子24远离波发送专用振子21而变弱。因此，波形累计部52对各振子24的波形信号施加适当的增益，按照各波形信号的振幅成为相同程度的方式进行调整。这里，所述增益的施加方法可以根据由窗口函数获得的波形信号的最大振幅来确定，也可以假定适当的指数函数的衰减来决定。
接着，波形累计部52使乘以所述窗口函数和增益而调整后的各波形信号移动，以使时间提前与各振子24的假定传播时间对应的时间(即，向图9的左侧偏移)。使各波形信号移动后，例如成为图10所示那样。另外，在图10中，各波形信号所含的表面折射波的峰值的相位一致。这样移动波形信号的结果峰值一致，是声速假定值与实际的骨声速一致的情况。
而且，波形累计部52对移动后的波形信号彼此进行累计，求出图11那样的累计波形。
最后，波形累计部52生成所述累计波形的包络线，并求出该包络线的振幅(参照图11)。如后述那样，该包络线的振幅成为偏移后的表面传播波的峰值是否相互一致(进一步说，假定传播时间是否与实际的骨声速一致)的指标。即，包络线的振幅是表示假定传播时间的妥当性的妥当性指标值。
下面，说明声速导出部53。声速导出部53构成为调用假定传播时间计算部51和波形累计部52的功能，反复求出所述包络线的振幅的处理(循环处理)。
在声速导出部53中进行的上述循环处理，假定传播时间计算部51每当被调用时从而与上次不同的声速假定值来计算假定传播时间。即，假定传播时间计算部51在每次所述循环中，采用新的声速假定值来求出新的t-x曲线。例如图12中示出使声速假定值变化而求出的多个t-x曲线。如图12的曲线图所示，声速假定值不同时，t-x曲线的倾斜、以及该t-x曲线与t坐标轴的交点等发生变化。因此，若声速假定值不同，则累计波形信号时，使各波形信号移动的量也不同。
使各波形信号移动的量变化时，累计波形也变化。因此，通过使声速假定值变化，所述包络线的振幅变化。声速导出部53通过边在规定范围内使声速假定值变化边反复所述循环，从而针对该规定范围内全部情况求出包络线的振幅值。在图13中示出这样边使声速假定值变化边求出包络线的振幅而描画的结果的例。
在图13的例中，假定声速为SOStrue时包络线的振幅成为最大。即，在该SOStrue时，认为各波形信号所含的表面折射波的峰值最加强(如图10那样，表面折射波的峰值的相位一致)。因此，声速导出部53将该SOStrue采用为声速的测定值。如以上那样，声速导出部53进行与声速导出工序对应的运算处理而导出声速。
下面，说明通过采用了上述骨强度诊断装置1的声速测定方法实际测定骨声速的実施例。本申请发明者为了确认上述声速测定方法的有效性，通过上述骨强度诊断装置1测定了牛胫骨皮质骨的声速。此外，为了评价其结果是否正确，将该牛胫骨切出并使其为块(ブロツク)状，用透过法测定长轴方向和圆周方向的声速，并进行了与骨强度诊断装置1的测定值的比较。图14中示出该结果。
首先，说明通过透过法测定声速的结果。已知牛的骨由骨单位(osteon)构成的Haversian构造、和层状的Plexiform构造构成。另外，由显微镜的观察结果，确认了在皮质骨的外侧多存在Plexiform构造，而在内侧多存在Haversian构造。因此，在透过法中，分别针对Plexiform构造和Haversian构造测定了圆周方向和长轴方向的声速。
另外，长管状的骨受到长边方向的负荷较多，所以一般在长轴方向被强化。因此，一般地，与圆周方向相比，长轴方向的一方声速大。如图14所示，在基于透过法的测定中也获得了长轴方向的声速大的结果。
下面，说明由本实施方式的骨强度诊断装置1测定声速的结果。如图14所示，由本实施方式的骨强度诊断装置1获得的长轴方向的声速与用透过法测定的结果几乎一致。此外，关于圆周方向的声速，基于本实施方式的骨强度诊断装置1的测定结果是Plexiform构造和Haversian构造的平均值。因此，可以确认通过本实施方式的骨强度诊断装置1可以适当地进行骨声速的测定。另外，一般认为人的骨大部分由Haversian构造构成。因此认为，对于人的骨，通过骨强度诊断装置1也可以进行适当的骨声速的测定。
如以上说明的那样，本实施方式的骨强度诊断装置1具备波发送专用振子21、多个振子24、假定传播时间计算部51、波形累计部52、声速导出部53。波发送专用振子21向皮质骨10的表面发送超声波。振子24接收来自皮质骨10的超声波，并输出与该接收的超声波相应的波形信号。假定传播时间计算部51根据该皮质骨10内的声速的假定值、该皮质骨10的表面形状计算波发送专用振子21发送的超声波在皮质骨10的表面附近传播之后到达各振子24为止的传播时间。波形累计部52根据多个振子24中至少2个振子输出的波形信号，求出成为所述传播时间的妥当性的指标的累计波形的包络线的振幅。声速导出部53根据累计波形的包络线的振幅，求出皮质骨10内的声速。
即，通过判断假定传播时间的妥当性，能够判断在计算该假定传播时间时所假定的声速的值是否正确，所以由此能够求出皮质骨10内的声速。这里，在判断假定传播时间的妥当性时，没有必要严密地检测接收波形的峰值的位置，所以即使在接收波形中含有噪声也能够适当地求出声速。此外，因为根据来自多个振子24的波形信号判断假定传播时间的妥当性，所以能够减弱噪声的影响。
此外，在本实施方式的骨强度诊断装置1中，波形累计部52使多个振子24中至少2个振子输出的波形信号分别移动与所述假定传播时间对应的时间，并求出累计该移动后的波形信号彼此的累计波形，求出该累计波形的包络线的振幅。
即，累计多个振子24波接收部输出的波形信号彼此时，波形的峰值位置一致时，最加强振幅。通过利用该性质，可以仅强调希望的峰值，并减少其他噪声的影响。在上述结构中，在假定传播时间计算部51假定的声速正确的情况下，被累计的波形信号的峰值位置一致，所以加强振幅。因此，即使在振子24输出的波形信号中含有噪声的情况下，通过观察累计波形，也能够判断声速的假定值是否正确。
此外，本实施方式的骨强度诊断装置1如下的构成。即，假定传播时间计算部51采用多个声速假定值来计算传播时间。波形累计部52根据多个声速假定值，针对各个情况的累计波形求出包络线的振幅。而且，声速导出部53通过比较多个声速假定值各个情况的包络线的振幅，从而求出皮质骨10内的声速。
由此，试行多个声速假定值，能够采用峰值彼此最加强的(峰值位置最一致)情况的声速假定值作为骨声速的测定值。
此外，在本实施方式的骨强度诊断装置1中，假定传播时间计算部51根据预先测定的皮质骨10表面的形状计算传播时间。
由此，可以考虑皮质骨10表面的实际形状来计算传播时间，所以即使例如皮质骨10表面是曲面状时，也可以精度良好地求出声速。
此外，本实施方式的骨强度诊断装置1具备形状检测部40。此外，多个振子24构成为能够通过向皮质骨10的表面同时发送超声波来产生平面波。形状检测部40根据平面波在皮质骨10的表面反射之后到达各振子24时的到达角度，检测皮质骨10表面的形状。而且，假定传播时间计算部51根据形状检测部40检测的皮质骨10表面的形状，计算所述传播时间。
由此，利用用于测定声速的结构的一部分，能够检测皮质骨10表面的形状，所以能够简化装置从而削减费用。
此外，本实施方式的声速测定方法包括波发送工序(S102)、波接收工序(S103)、假定传播时间计算工序(S104)、波形累计工序(S105)。在波发送工序中，向皮质骨10的表面发送超声波。在波接收工序中，由各振子24接收来自该皮质骨10的超声波。在假定传播时间计算工序中，根据该皮质骨10内的声速的假定值、该皮质骨10的表面形状计算在波发送工序发送的超声波在皮质骨10的表面附近传播之后到达各振子24为止的传播时间。在波形累计工序中，根据多个振子24中至少2个振子输出的波形信号，求出成为所述传播时间的妥当性的指标的累计波形的包络线的振幅。而且，通过使声速假定值变化并反复假定传播时间计算工序和所述波形累计工序(S104～S106的循环)，求出多个声速假定值各自情况的累计波形，通过比较所述各个情况的累计波形的包络线的振幅，求出皮质骨10内的声速(S107)。
即，通过判断传播时间的妥当性，能够判断计算该传播时间时所假定的声速的值是否正确，所以由此能够求出皮质骨10内的声速。这里，在判断传播时间的妥当性时，因为不需要严密地检测接收波形的峰值，所以即使接收波形中含有噪声也能够适当地求出声速。此外，因为根据来自多个振子24的波形信号判断传播时间的妥当性，所以能够减弱噪声的影响。而且，实行多个声速的假定值，在传播时间最妥当时，能够采用计算该传播时间时的声速的假定值作为被测体内的声速的测定值。
此外，本实施方式的声速测定方法采用以下的方法。即，在波形累计工序中，使所述多个波接收部中至少2个波接收部输出的波形信号移动与所述传播时间对应的时间，求出对该移动后的波形信号彼此进行累计后的累计波形，求出该累计波形的包络线的振幅。
即，累计多个振子24输出的波形信号彼此时，在波形的峰值位置一致时最加强振幅。通过利用该性质，可以仅强调希望的峰值，并减少其他噪声的影响。在上述结构中，在假定到达时间计算工序所假定的声速正确时，所累计的波形信号的峰值位置一致，所以加强振幅。因此，即使在振子24输出的波形信号中含有噪声时，通过调查累计波形，能够判断声速的假定值是否正确。
此外，本实施方式的声速测定工序包括检测皮质骨10的表面形状的形状检测工序(S101)。而且，在假定传播时间计算工序中，根据所述表面形状计算传播时间。
由此，可以考虑皮质骨10表面的实际形状来计算假定传播时间，所以即使在例如皮质骨10表面是曲面状的情况下，也可以精度良好地求出声速。
下面，说明上述实施方式的骨强度诊断装置1的变形例。
该变形例涉及的骨强度诊断装置具备波形乘法部作为妥当性指标值计算部。即，在上述实施方式中，作为妥当性指标值计算部的波形累计部52使各振子24的波形信号根据假定传播时间移动，并累计移动后的波形信号彼此从而求出累计波形。对于该点，本变形例的波形乘法部构成为根据假定传播时间使各振子24的波形信号移动，将移动后的波形信号彼此相乘从而求出乘法波形，并根据该乘法波形求出妥当性指标值。
如以上说明的那样，在该变形例涉及的声速测定装置中，波形乘法部使多个振子24中至少2个振子输出的波形信号分别移动与假定传播时间对应的时间，并求出将该移动后的波形信号彼此相乘所得的乘法波形，根据该乘法波形来求出妥当性指标值。
此外，在采用了该变形例涉及的声速测定装置的声速测定方法中，作为妥当性指标值计算工序，代替上述实施方式的波形累计工序而进行波形乘法工序。即，在波形乘法工序中，将多个波接收部中至少2个波接收部输出的波形信号分别移动与所述传播时间对应的时间，并求出将该移动后的波形信号彼此相乘所得的乘法波形，根据该乘法波形来求出妥当性指标值。
即，多个振子24波接收部输出的波形信号彼此相乘时，波形的峰值位置一致时最加强振幅。通过利用该性质，可以仅强调希望的峰值，并减少其他噪声的影响。在上述结构中，假定传播时间计算部51所假定的声速正确时，相乘所得的波形信号的峰值位置一致，所以加强振幅。因此，即使在振子24输出的波形信号中含有噪声时，通过观察乘法波形，也能够判断声速的假定值是否正确。
下面，说明上述实施方式的骨强度诊断装置1的其他变形例。
该变形例涉及的骨强度诊断装置1具备的超声波收发器90的构造如图15所示。该超声波收发器90省略波发送专用振子，仅具备多个振子24以大致等间隔排列为1列的阵列振子22。这里，阵列振子22设置为延长到上述实施方式中配置波发送专用振子21的位置。
这里，在上述实施方式的骨强度诊断装置1中，具有以下这样的课题。即，在上述实施方式中，在检测骨表面线时，超声波对皮质骨10的射入点Pin附近的形状未能检测。这时因为在由形状检测部40检测皮质骨表面的形状时，仅能检测阵列振子22的正下方的形状。因此，在计算表面传播波的传播路径时，根据由骨表面线预测的线求出了射入点Pin的位置，但是在该方法中，不能求出正确的射入点Pin的位置，所以也涉及导出的声速的误差。
因此，本变形例省略波发送专用振子，进行声速测定时的超声波的发送也由阵列振子(线性阵列(linear array))22来进行。例如，采用振子24a和振子24b作为上述实施方式中的波发送专用振子21的代替。
即，通过相邻的振子24错开定时而发送超声波，从而能够向倾斜方向发送超声波束。据此，阵列振子22不仅能兼具波发送专用振子21的功能，还能够以任意的角度发送超声波，所以能够对皮质骨10以适当的角度(临界角或临界角附近的角度)照射超声波。由此，效率良好，还能降低由指向性引起的声速误差。
而且，根据上述那样地构成，对皮质骨10射入超声波的射入点的位置总是阵列振子22的正下方的位置，所以根据基于形状检测部40的骨表面线的检测结果，能够精度良好地计算该射入点的位置。此外，本变形例的超声波收发器90如图15所示跨整个宽度配置阵列振子22，所以测定时能够把握骨形状的整体，测定变得容易进行。
如以上说明的那样，在该变形例涉及的声速测定装置具备的超声波收发器中，多个振子24构成以等间隔排列为1列的阵列振子22，并且多个振子中至少任一个兼具上述实施方式中的波发送专用振子21的功能。
由此，能够省略波发送专用的结构，所以能够简化装置。此外，通过使从相邻振子24发送超声波的定时错开，并调整错开该定时的间隔，从而能够以任意的角度形成超声波束。由此，能够对皮质骨10的表面以最合适的角度照射超声波。此外，在由阵列振子22产生平面波从而检测皮质骨10的表面形状的结构的情况下，能检测该表面形状的位置仅是阵列振子22的正下方的位置。因此，在分别构成波发送专用振子21和阵列振子22的情况下，不能检测向皮质骨表面的射入点附近的形状。对于这点，通过振子24兼具作为波发送专用振子21的功能，从而能够检测向皮质骨表面的射入点附近的形状。由此，能够正确地计算传播时间。
以上，说明了本发明优选的实施方式及其变形例，但是上述结构可以例如以下那样变更。
在上述实施方式中，说明了全部累计多个振子24输出的波形信号，但是不需要一定累计全部波形信号。例如，在预先判断出表面传播波不会到达某振子24的情况下，也可以不累计来自该振子24的波形信号。
此外，例如，直接波、来自表面的反射波、来自背面的反射波和表面传播波的峰值有成为图16所示的关系的情况。如图16所示，表面传播波和来自表面的反射波从某位置起连续地分离，表面传播波先到达。与此相对，来自背面的反射波从某位置起逐渐接近表面传播波从而重叠。因此，在图16的曲线图中在用虚线包围的部分最良好地含有表面传播波的信息(在该部分，表面传播波的峰值与其他波的峰值不重叠)。因此，在累计波形信号时，若仅选择包括用该虚线包围的部分的波形信号来累计，则可以进行精度良好的声速测定。
另外，表面传播波和来自表面的反射波分离的位置能够根据皮质骨的表面形状进行预测，所以可以利用在选择进行累计的波形信号时。另一方面，表面传播波和来自背面的反射波重叠的位置难以预测。但是，对来自表面传播波与来自表面的反射波分离的位置的全部波形信号进行累计(对来自背面的反射波重叠的波形信号进行累计)，根据本发明的声速测定方法，也可以抑制噪声波形的影响，所以能够适当地导出声速。
也可以变更为在超声波收发器侧具备上述运算部。此外，不限于超声波收发器和装置主体分别的结构，也可以使超声波收发器和装置主体一体化。
在上述实施方式中，假定传播时间计算部51是求出折射表面波的传播路径以及假定传播时间的结构，但是也可以构成为求出泄露表面波的传播路径和假定传播时间。即使在该情况下，通过累计波形信号彼此，也能够使焦点与该泄露表面波的峰值一致，所以能够适当地测定声速。
形状检测部40根据表面反射波检测皮质骨10的表面形状时，也可以根据背面反射波检测皮质骨10的背面形状。在该情况下，根据表面形状和背面形状，能够求出皮质骨厚度。所求出的皮质骨厚度可以用作骨的健全性的指标的一种。
此外，基于形状检测部40的被测体的表面形状的检测方法，也可以代替从多个振子24同时进行波发送从而产生平面波的方法，而采用例如以下的方法。即，从阵列振子22的一端开始，由一个或多个振子24依次进行超声波的发送。然后，根据从进行了各个波发送之后开始到接收表面反射波为止的时间，求出皮质骨10的表面形状。
根据上述形状检测方法，因为从多个振子24依次发送超声波，所以在测定中多少花费时间，但是即使在被测体的表面是凹状的情况下，也可以正确地检测表面形状。另一方面，如实施方式中说明的那样，在从多个振子24同时进行波发送从而产生平面波的方法中，能够缩短测定时间，但是在被测体的表面是凹状的情况下，不能测定表面形状。胫骨剖面的形状各种各样，在一部分中包括凹形状的情况也有，所以在上述实施方式的骨强度诊断装置1中，有不能正确地求出皮质骨10的表面形状的情况。
因此，特别优选构成为用产生平面波的上述实施方式的方法测定皮质骨10的表面形状，在产生错误时(在皮质骨10的剖面形状中包括凹部时)，采用从阵列振子22的一端的振子24开始依次发送超声波的上述方法。
在上述实施方式中，波形累计部52构成为求出累计波形的包络线的振幅作为妥当性指标值，但不限于此。例如，即使是将累计波形的振幅本身作为妥当性指标值的结构，也能够检测波形信号的峰值是否加强(即，声速假定值是否正确)。此外例如，也可以构成为将累计波形的累计值(面积)作为妥当性指标值。
此外，在上述实施方式和变形例中，构成为使波形信号移动之后求出累计波形或乘法波形，之后判断假定传播时间的妥当性，但是不限于该结构。即，也可以不使各振子24输出的波形信号移动，而与假定传播时间曲线(t-x曲线)直接比较从而决定声速。在该情况下，可以原样利用波形信号，也可以进行取出振幅的处理。
例如，将图3所示的波形信号的数据看做二维图像，与图8所示的t-x曲线比较，求出表示一致度的指标(妥当性指标值)。作为上述那样将波形信号和t-x曲线直接比较的方法，可以利用以往的数字图像处理中所用的特征提取法(一般化半变换等)、以及图案识别技术(利用了神经网络的技术等)。而且，使声速假定值变化，并且求出对图12所示那样的各种t-x曲线的一致度，采用一致度最高的t-x曲线的声速假定值作为声速的测定值。
即使是以上的方法，也可以与使波形信号移动之后进行累计或乘法的方法同样地抑制噪声的影响从而进行稳定的声速测定。即，即使在各个波形信号中含有噪声，也可以利用来自多个振子24的波形信号，进行上述特征提取或图案识别等，所以可以降低噪声的影响，良好地求出妥当性指标值。此外，因为是观察与t-x曲线的一致度的结构，所以即使不由波形信号严密地求出峰值位置，也可以判断假定传播时间是否妥当。
对于软组织中的声速，也可以构成为根据软组织的厚度而变更。另外，一般而言，厚的软组织含有脂肪多，所以软组织越厚，该软组织中的声速越延迟。
在上述实施方式中，在各振子24输出的波形信号中乘以了窗口函数，但是也可以省略使用窗口函数的运算。
此外，也可以省略所述形状检测工序，例如假定皮质骨10的剖面轮廓形状是直线状，进行假定传播时间的计算。但是，在骨的圆周方向测定声速时，皮质骨10表面的剖面轮廓是曲线状，所以若如上述那样假定传播路径，则误差变大。因此，优选检测皮质骨表面的形状并根据该表面形状来进行假定传播时间的计算。
如图13的曲线图所示，累计波形的包络线的振幅取得极大值的声速假定值有产生多个的情况。在该情况下，作为采用哪个作为测定值的判断材料，可以辅助性地使用累计波形的振幅成为最大的时间。图17是对图13的曲线图追加描画了累计波形的振幅成为最大的时间的曲线图。在有两个累计波形的包络线的振幅取极大值的声速假定值(如图17的曲线图所示)的情况下，可以认为一个是表面传播波的峰值加强而包络线的振幅变大的情况，另一个是例如来自背面的反射波的峰值加强而包络线的振幅变大的情况。如图17那样描画累计波形的振幅成为最大的时间(即，检测出加强的峰值的时间)时，可以容易地确认切换了被增强了的峰值。这也可以作为焦点是否与表面传播波一致的指标。
此外，本发明的声速测定方法和声速测定装置还能够广泛地使用于测定骨以外的声速的用途。尤其，本发明的声速测定方法和声速测定装置在噪声多的环境下且在测定表面形状是非直线状的被测体的声速时极为有效。

Claims (11)

1.一种声速测定装置，其特征在于，具备：
波发送部，其向被测体的表面发送超声波；
多个波接收部，接收来自所述被测体的超声波，并输出与该接收到的超声波对应的波形信号；
假定传播时间计算部，其根据该被测体内的声速的假定值和该被测体的表面形状，计算所述波发送部发送出的超声波在所述被测体的表面或该被测体的内部传播之后到达各所述波接收部为止的传播时间；
妥当性指标值计算部，其根据所述多个波接收部中的至少两个波接收部输出的波形信号，求出成为所述传播时间的妥当性的指标的妥当性指标值；和
声速导出部，其根据所述妥当性指标值，求出所述被测体内的声速。
2.根据权利要求1所述的声速测定装置，其特征在于，
所述妥当性指标值计算部使所述多个波接收部中的至少两个波接收部输出的波形信号分别移动与所述传播时间对应的时间，求出对该移动后的波形信号彼此进行累计所得的累计波形，并根据该累计波形来求出所述妥当性指标值。
3.根据权利要求1所述的声速测定装置，其特征在于，
所述妥当性指标值计算部使所述多个波接收部中的至少两个波接收部输出的波形信号分别移动与所述传播时间对应的时间，求出将该移动后的波形信号彼此相乘所得的乘法波形，并根据该乘法波形来求出所述妥当性指标值。
4.根据权利要求1～3中任意一项所述的声速测定装置，其特征在于，
所述假定传播时间计算部采用多个值来计算所述传播时间作为所述声速的假定值，
所述妥当性指标值计算部根据多个声速的假定值，求出各种情况下的妥当性指标值，
所述声速导出部根据所述多个声速的假定值各自情况的妥当性指标值，求出所述被测体内的声速。
5.根据权利要求1～3中任意一项所述的声速测定装置，其特征在于，
所述假定传播时间计算部根据预先测定的被测体表面的形状，计算所述传播时间。
6.根据权利要求1～3中任意一项所述的声速测定装置，其特征在于，
具备形状检测部，
所述多个波接收部中的至少一部分构成为能够向被测体的表面发送超声波，
所述形状检测部根据从超声波被发送之后开始到在所述被测体的表面反射之后到达各所述波接收部为止的时间，检测所述被测体表面的形状，
所述假定传播时间计算部根据所述形状检测部检测出的所述被测体表面的形状，计算所述传播时间。
7.根据权利要求1～3中任意一项所述的声速测定装置，其特征在于，
所述多个波接收部构成为以大致等间隔排列为一列的线性阵列，并且所述多个波接收部中的至少任意一个波接收部作为所述波发送部而发挥功能。
8.一种声速测定方法，利用具备了根据接收到的超声波来输出波形信号的多个波接收部的声速测定装置，其特征在于，包括：
波发送工序，向被测体的表面发送超声波；
波接收工序，由各所述波接收部接收来自该被测体的超声波；
假定传播时间计算工序，根据该被测体内的声速的假定值和该被测体的表面形状，计算在所述波发送工序发送出的超声波在所述被测体的表面或该被测体的内部传播之后到达各所述波接收部为止的传播时间；和
妥当性指标值计算工序，根据所述多个波接收部中的至少两个波接收部输出的波形信号，求出成为所述传播时间的妥当性的指标的妥当性指标值，
使所述声速的假定值变化，并反复进行所述假定传播时间计算工序和所述妥当性指标值计算工序，由此求出多个声速的假定值各自情况的妥当性指标值，根据所述各自情况的妥当性指标值，求出所述被测体内的声速。
9.根据权利要求8所述声速测定方法，其特征在于，
在所述妥当性指标值计算工序中，将所述多个波接收部中的至少两个波接收部输出的波形信号移动与所述传播时间对应的时间，求出将该移动后的波形信号彼此进行累计所得的累计波形，根据该累计波形求出所述妥当性指标值。
10.根据权利要求8所述声速测定方法，其特征在于，
在所述妥当性指标值计算工序中，将所述多个波接收部中的至少两个波接收部输出的波形信号移动与所述传播时间对应的时间，求出将该移动后的波形信号彼此相乘所得的乘法波形，根据该乘法波形求出所述妥当性指标值。
11.根据权利要求8～10中任意一项所述声速测定方法，其特征在于，
包括形状检测工序，检测所述被测体的表面形状，
在所述假定传播时间计算工序中，根据所述表面形状计算所述传播时间。
CN201010163281.4A 2009-04-14 2010-04-14 声速测定装置以及声速测定方法 Expired - Fee Related CN101865931B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2009098414A JP5280927B2 (ja) 2009-04-14 2009-04-14 音速測定装置及び音速測定方法
JP2009-098414 2009-04-14

Publications (2)

Publication Number Publication Date
CN101865931A CN101865931A (zh) 2010-10-20
CN101865931B true CN101865931B (zh) 2013-03-13

Family

ID=42315830

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201010163281.4A Expired - Fee Related CN101865931B (zh) 2009-04-14 2010-04-14 声速测定装置以及声速测定方法

Country Status (5)

Country Link
US (1) US8256295B2 (zh)
EP (1) EP2241260B1 (zh)
JP (1) JP5280927B2 (zh)
CN (1) CN101865931B (zh)
AT (1) ATE533406T1 (zh)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5584300B2 (ja) * 2010-08-31 2014-09-03 株式会社日立メディコ 画像再構成方法及び装置
JP2012189352A (ja) * 2011-03-09 2012-10-04 Ihi Inspection & Instrumentation Co Ltd 表面を伝播する超音波の音速測定装置と方法
WO2013008746A1 (ja) * 2011-07-14 2013-01-17 株式会社メガチップス 直線検出装置および直線検出方法
JP2013088118A (ja) * 2011-10-13 2013-05-13 Ihi Inspection & Instrumentation Co Ltd ガイド波を用いた検査方法
EP2781192A4 (en) * 2011-11-17 2015-08-12 Hitachi Aloka Medical Ltd AN ECHOGRAPHIC DIAGNOSTIC DEVICE AND AN EPOGRAPHIC IMAGE GENERATION METHOD
WO2013175867A1 (ja) * 2012-05-21 2013-11-28 古野電気株式会社 伝播速度測定装置、伝播速度測定プログラム、及び伝播速度測定方法
WO2014038569A1 (ja) * 2012-09-10 2014-03-13 古野電気株式会社 厚み測定装置及び厚み測定方法
CN104837411B (zh) * 2012-12-06 2017-03-29 古野电气株式会社 测量装置及测量方法
US9912415B2 (en) * 2013-11-12 2018-03-06 Qualcomm Incorporated Fast service discovery and pairing using ultrasonic communication
JP2015173879A (ja) * 2014-03-17 2015-10-05 古野電気株式会社 音速測定装置、音速測定方法、及び音速測定プログラム
JP2015232519A (ja) * 2014-06-10 2015-12-24 アズビル株式会社 クランプオン式超音波流量計及び流量の計測方法
CN104483385B (zh) * 2014-12-05 2017-02-22 中国航空工业集团公司北京航空材料研究院 一种各向异性材料纵波声速的测量方法
CN108348221B (zh) * 2015-09-01 2021-02-19 戴尔菲纳斯医疗科技公司 使用超声波波形断层成像的组织成像和分析
CN108601584B (zh) * 2015-12-04 2021-05-28 古野电气株式会社 超声波解析装置、超声波解析方法及存储介质
CN106525979A (zh) * 2016-11-16 2017-03-22 中国矿业大学(北京) 基于波形包络的各向同性介质超声波速度自动计算方法
WO2018234472A1 (en) * 2017-06-21 2018-12-27 Charité - Universitätsmedizin Berlin SYSTEM, METHOD AND COMPUTER PROGRAM PRODUCT FOR DETERMINING THE CHARACTERISTICS OF A CORTICAL BONE
CN113109446B (zh) * 2021-04-15 2022-11-29 复旦大学 一种超声断层成像方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4279157A (en) * 1978-06-22 1981-07-21 U.S. Philips Corporation Method of and device for determining the internal structure of a body by means of acoustic beams
CN1732853A (zh) * 2004-02-23 2006-02-15 计算机逻辑公司 超声波骨头评估的装置及方法
CN1766530A (zh) * 2005-11-29 2006-05-03 澳诺制药有限公司 用双向逼近法测定骨中声速的方法
CN1846628A (zh) * 2006-04-29 2006-10-18 澳诺制药有限公司 用时间极值法测定超声波在骨骼中传播速度的方法

Family Cites Families (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5426979A (en) * 1990-06-04 1995-06-27 Medicano Systems Ltd. Frequency spectrum apparatus for determining mechanical properties
JPH06109709A (ja) * 1992-08-31 1994-04-22 Medicano Syst Ltd 中実物の力学的性質判定装置とその方法、並びに中実物の走査及び力学的性質判定装置とその方法
JP2575583B2 (ja) * 1993-03-16 1997-01-29 アロカ株式会社 生体組織中の超音波伝搬特性に関する係数値の推定方法
JP2840040B2 (ja) * 1994-12-22 1998-12-24 アロカ株式会社 組織内音速測定方法
CA2233852C (en) * 1995-10-04 2004-03-16 Sunlight Ultrasound Technologies Limited Ultrasonic device for determining bone characteristics
EP0821913A4 (en) * 1996-02-21 1999-04-28 Seikisui Chemical Co Ltd METHOD AND OSTEOPOROSIS DIAGNOSTIC DEVICE
JP2003517328A (ja) * 1998-03-03 2003-05-27 サンライト メディカル リミテッド 骨内の音響速度決定
JP4688262B2 (ja) * 2000-07-27 2011-05-25 アロカ株式会社 超音波診断装置
FR2839877B1 (fr) * 2002-05-27 2005-04-08 Centre Nat Rech Scient Procede, sonde et appareil pour evaluer de maniere non-invasive une duree de parcours ou une vitesse d'ultra-sons le long d'une interface, notamment osseuse
JP2005043164A (ja) * 2003-07-25 2005-02-17 Toshiba Corp 音波伝播時間測定装置
DE10338940B3 (de) * 2003-08-22 2005-02-10 Fresenius Medical Care Deutschland Gmbh Verfahren zum Messen einer Signallaufzeit in Blut und Vorrichtung zur Anwendung des Verfahrens
GB2424276B (en) * 2005-03-17 2008-09-17 Furuno Electric Co Ultrasonic bone evaluation apparatus
US7753848B2 (en) * 2005-11-11 2010-07-13 Brainlab Ag Determining speed-of-sound factors in ultrasound images of a body
JP5280647B2 (ja) * 2007-05-29 2013-09-04 古野電気株式会社 超音波を用いた骨強度診断装置及び超音波を用いた骨強度診断装置が作動する方法
JP4994182B2 (ja) 2007-10-17 2012-08-08 アルパイン株式会社 オーディオ装置
JP2009153945A (ja) * 2007-12-28 2009-07-16 Furuno Electric Co Ltd 骨厚み測定装置
JP5374086B2 (ja) * 2008-07-25 2013-12-25 古野電気株式会社 骨強度診断装置及び骨強度測定方法
JP5235112B2 (ja) * 2008-07-25 2013-07-10 古野電気株式会社 形状検出装置、形状検出方法、及び形状検出装置を用いた骨強度診断装置
JP2010099452A (ja) * 2008-09-25 2010-05-06 Fujifilm Corp 超音波診断装置及び超音波診断方法
JP5174604B2 (ja) * 2008-09-30 2013-04-03 富士フイルム株式会社 超音波信号処理装置及び方法
JP2010234013A (ja) * 2009-03-31 2010-10-21 Fujifilm Corp 超音波診断装置及び超音波診断方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4279157A (en) * 1978-06-22 1981-07-21 U.S. Philips Corporation Method of and device for determining the internal structure of a body by means of acoustic beams
CN1732853A (zh) * 2004-02-23 2006-02-15 计算机逻辑公司 超声波骨头评估的装置及方法
CN1766530A (zh) * 2005-11-29 2006-05-03 澳诺制药有限公司 用双向逼近法测定骨中声速的方法
CN1846628A (zh) * 2006-04-29 2006-10-18 澳诺制药有限公司 用时间极值法测定超声波在骨骼中传播速度的方法

Also Published As

Publication number Publication date
EP2241260B1 (en) 2011-11-16
CN101865931A (zh) 2010-10-20
JP2010246692A (ja) 2010-11-04
US20100257935A1 (en) 2010-10-14
ATE533406T1 (de) 2011-12-15
JP5280927B2 (ja) 2013-09-04
EP2241260A1 (en) 2010-10-20
US8256295B2 (en) 2012-09-04

Similar Documents

Publication Publication Date Title
CN101865931B (zh) 声速测定装置以及声速测定方法
JP5847941B2 (ja) 波形トラッキング装置、超音波診断装置、及び波形トラッキング方法
EP1194920B1 (en) Recursive ultrasound imaging
CN104203112B (zh) 超声波诊断装置以及弹性评价方法
CN103181780B (zh) 信号分析方法、超声波图像分析方法以及超声波成像系统
CN103462643B (zh) 一种剪切波速度测量方法、装置和系统
EP2881041A1 (en) Apparatus and method for ultrasonic diagnosis
US20210356573A1 (en) Two dimension and three dimension imaging based on speed changes of sound/ultrasound
US9465063B2 (en) Method and system for the estimation and cancellation of multipath delay of electromagnetic signals, in particular SSR replies
JP6134323B2 (ja) 厚み測定装置及び厚み測定方法
US6322507B1 (en) Ultrasonic apparatus and method for evaluation of bone tissue
CN105877783B (zh) 二维剪切波弹性成像方法和装置
WO2008146513A1 (ja) 骨強度診断装置及び骨強度診断方法
JP5801956B2 (ja) 伝播速度測定装置、伝播速度測定プログラム、及び伝播速度測定方法
CN105116371B (zh) 一种基于连续发射调频信号的目标定位方法与装置
CN109589138A (zh) 一种剪切波波速计算方法及弹性成像设备
CN104837411B (zh) 测量装置及测量方法
KR101055580B1 (ko) Ｂｃ-모드 영상을 형성하는 초음파 시스템 및 방법
Granja et al. Ultrasound modulation and codification for localization systems
JP2015173879A (ja) 音速測定装置、音速測定方法、及び音速測定プログラム
JP6633854B2 (ja) 超音波測定装置
CN104297756B (zh) 水下旋转目标缩比声纳成像方法及系统
CN104619264B (zh) 厚度测量装置及厚度测量方法
IL142447A (en) Ultrasonic apparatus and method for evaluation of bone tissue

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130313

Termination date: 20180414