JP4312494B2 - Ultrasonic bone density measuring device, ultrasonic measuring device, and ultrasonic measuring method - Google Patents

Ultrasonic bone density measuring device, ultrasonic measuring device, and ultrasonic measuring method Download PDF

Info

Publication number
JP4312494B2
JP4312494B2 JP2003111647A JP2003111647A JP4312494B2 JP 4312494 B2 JP4312494 B2 JP 4312494B2 JP 2003111647 A JP2003111647 A JP 2003111647A JP 2003111647 A JP2003111647 A JP 2003111647A JP 4312494 B2 JP4312494 B2 JP 4312494B2
Authority
JP
Japan
Prior art keywords
ultrasonic
transfer function
waveform
signal
deriving
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
JP2003111647A
Other languages
Japanese (ja)
Other versions
JP2004344181A (en
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
Priority to JP2003111647A priority Critical patent/JP4312494B2/en
Priority to PCT/JP2004/005455 priority patent/WO2004091409A1/en
Publication of JP2004344181A publication Critical patent/JP2004344181A/en
Application granted granted Critical
Publication of JP4312494B2 publication Critical patent/JP4312494B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

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)

Description

【０００１】
【発明の属する技術分野】
本発明は、骨粗鬆症の診断等に利用される超音波骨密度測定装置、超音波測定装置、超音波骨密度測定方法及び超音波測定方法に関する。
【０００２】
【従来の技術】
骨組織の密度が低下する症状である骨粗鬆症の診断のために、超音波骨密度測定装置が用いられている。超音波は、被検体に含有される骨塩量に応じた速度で被検体内を伝播する。したがって、被検体内における超音波の伝播速度（ＳＯＳ：Speed Of Sound）を測定すれば、被検体である踵の骨塩量に相当する量を測定することができる。超音波骨密度測定装置は、この原理を利用したものであって、被検者の骨中における超音波の伝播速度を測定し、測定された超音波の伝播速度から骨塩密度などの骨特性を定量測定することができる装置である。超音波伝播速度の測定精度を向上させるために、超音波骨密度測定装置による超音波検査は、海綿骨が多い踵骨を被検体として行われることが多い。
【０００３】
特許文献１に記載された超音波骨密度測定装置は、被検体を載せるための基台と、両側から被検体を挟み込んで被検体に密着する一対の測定体とを有している。一方の測定体内には、超音波振動子を有する超音波送波器が配置され、他方の測定体内には、超音波振動子を有する超音波受波器（実質的に超音波送波器と同じもの）が配置されている。超音波送波器から送波された超音波は被検体を通過して超音波受波器によって受波される。また、測定体と被検体との境界面で反射する超音波を測定し、通過した超音波との時間差から被検体の超音波伝播速度を求める。
【０００４】
【特許文献１】
特開平１０−４３１８０号公報
【０００５】
特許文献１に記載されたような超音波骨密度測定装置においては、超音波伝播速度の測定時に、超音波送波器に単パルス信号を与えている。このとき、超音波振動子の周波数帯域が狭いために、超音波送波器から出力される超音波は時間軸に沿って広がった（時間幅の広い）波形を有するものとなる。したがって、超音波骨密度測定装置の時間分解能が悪化してしまう。超音波振動子としてコンポジット振動子を用いれば超音波送波器から出力される超音波を時間幅の短いものとすることが可能であるが、その時間幅短縮効果は十分ではなく、しかもコンポジット振動子自体が高価であって商用の装置に用いることは実用的ではない。また、超音波送波器から出力される超音波は、一部が被検体を挟み込む測定体と被検体との境界面以外の表面で反射し、時間的に遅れて受波器に到達することがあり、その場合には測定精度が悪化してしまう。この反射波（以降反射波と記載）は測定体の形状に依存し、被検体を通過した超音波を受波する場合も、測定体と被検体との境界面で反射した超音波を受波する場合のどちらの場合でも発生し得る。
【０００６】
また、超音波骨密度測定装置を含む超音波測定装置において、ノイズの影響を可能な限り排除して高精度の測定結果を得られるようにすることが望まれている。
【０００７】
そこで、本発明の一つの目的は、送波される超音波の時間幅延長を抑制し且つ反射波の影響を少なくすることで高精度の測定結果を得ることが可能な超音波骨密度測定装置及び超音波骨密度測定方法を提供することである。
【０００８】
また、本発明の別の目的は、ノイズの影響を受けにくい高精度の測定結果を得ることが可能な超音波測定装置及び超音波測定方法を提供することである。
【０００９】
【課題を解決するための手段】
本発明の超音波骨密度測定装置は、被検体を載せるための基台と、少なくともいずれか一方が他方の方向に移動可能に対向した、超音波透過性固体材料からなる一対の測定体と、前記一対の測定体のいずれか一方に設けられた超音波送波器と、前記一対の測定体のいずれか一方に設けられた超音波受波器と、前記超音波送波器から前記超音波受波器に至る測定系の周波数伝達関数を導出するための伝達関数導出手段と、出力信号の所望波形と、前記伝達関数導出手段によって導出された周波数伝達関数とに基づいて、前記超音波受波器が前記所望波形を有する出力信号を出力するときに前記超音波送波器に供給される入力信号の波形を導出するための入力波形導出手段と、前記入力波形導出手段が導出した波形を有する入力信号を前記超音波送波器に供給するための入力波形供給手段とを備えている。
【００１０】
本発明の超音波測定方法は、少なくともいずれか一方が他方の方向に移動可能に対向した、超音波透過性固体材料からなる一対の測定体のいずれか一方に設けられた超音波送波器から、前記一対の測定体のいずれか一方に設けられた超音波受波器に至る測定系の周波数伝達関数を導出する伝達関数導出ステップと、出力信号の所望波形と、前記伝達関数導出ステップで導出された周波数伝達関数とに基づいて、前記超音波受波器が前記所望波形を有する出力信号を出力するときに前記超音波送波器に供給される入力信号の波形を導出する入力波形導出ステップとを備えている。
【００１１】
本発明の超音波測定方法は、少なくともいずれか一方が他方の方向に移動可能に対向しており且つ先端にスタンドオフが設けられた超音波透過性固体材料からなる一対の測定体のいずれか一方に設けられた超音波送波器から、前記一対の測定体のいずれか一方に設けられた超音波受波器に至る測定系の周波数伝達関数を、２つの前記スタンドオフを接触させた状態において導出する伝達関数導出ステップと、出力信号の所望波形と、前記伝達関数導出ステップで導出された周波数伝達関数とに基づいて、前記超音波受波器が前記所望波形を有する出力信号を出力するときに前記超音波送波器に供給される入力信号の波形を導出する入力波形導出ステップと、２つの前記スタンドオフを離隔させるステップと、離隔した２つの前記スタンドオフ間に配置された踵を２つの前記スタンドオフで挟み込むステップと、２つの前記スタンドオフによって踵が挟み込まれた状態において、前記入力波形導出ステップで求められた波形を有する入力信号を前記超音波送波器に供給し、前記超音波受波器からの出力信号を検出するステップとを備えている。本明細書において、「スタンドオフ」とは患者に接触する部分を言う。
【００１２】
この構成によると、周波数伝達関数を用いて超音波受波器から出力される出力信号が所望波形を有するような入力信号の波形を導出しているので、超音波受波器から出力される出力信号の波形を時間幅の狭いものとすることが可能である。また、導出された伝達関数が反射波の影響を組み込んだものとなるので、導出された波形を有する入力信号を超音波送波器に供給したときに超音波受波器から出力される出力信号が反射波成分をほとんど含まなくなる。したがって、高精度の測定結果を得ることが可能となる。
【００１３】
本発明において、各測定体の形状が、前記超音波送波器から送波された超音波の超音波送波器から超音波受波器までの反射波の影響を前記伝達関数に含めて所望波形の出力信号を得るように、被検体の測定に適した形状に決定されていることが好ましい。これによると、測定体がいかなる形状であっても測定体内部での反射波の影響を受けないため、自由に被検体の測定に最も適した形状に決定することができる。
【００１４】
本発明の超音波測定装置は、超音波送波器と、超音波受波器と、前記超音波送波器から前記超音波受波器に至る測定系の周波数伝達関数を導出するための伝達関数導出手段と、出力信号の所望波形と、前記伝達関数導出手段によって導出された周波数伝達関数とに基づいて、前記超音波受波器が前記所望波形を有する出力信号を出力するときに前記超音波送波器に供給される入力信号の波形を導出するための入力波形導出手段と、前記入力波形導出手段が導出した波形を有する入力信号を前記超音波送波器に供給するための入力波形供給手段とを備えている。そして、前記伝達関数導出手段が、超音波を受波した前記超音波受波器からの出力信号にインパルス信号を加えた信号に基づいて、測定系の周波数伝達関数を導出する。
【００１５】
本発明の超音波骨密度測定装置は、別の観点において、超音波送波器と、超音波受波器と、前記超音波送波器から被検査物を経て前記超音波受波器に至る測定系の周波数伝達関数を導出するための伝達関数導出手段と、出力信号の所望波形と、前記伝達関数導出手段によって導出された周波数伝達関数とに基づいて、前記超音波受波器が前記所望波形を有する出力信号を出力するときに前記超音波送波器に供給される入力信号の波形を導出するための入力波形導出手段と、前記入力波形導出手段が導出した波形を有する入力信号を前記超音波送波器に供給するための入力波形供給手段とを備えている。そして、前記伝達関数導出手段が、超音波を受波した前記超音波受波器からの出力信号にインパルス信号を加えた信号に基づいて、測定系の周波数伝達関数を導出する。
【００１６】
本発明の超音波測定方法は、超音波送波器から前記超音波受波器に至る測定系の周波数伝達関数を導出する伝達関数導出ステップと、超音波を受波した前記超音波受波器からの出力信号の所望波形と、前記伝達関数導出ステップで導出された周波数伝達関数とに基づいて、前記超音波受波器が前記所望波形を有する出力信号を出力するときに前記超音波送波器に供給される入力信号の波形を導出する入力波形導出ステップとを備えている。そして、前記伝達関数導出ステップにおいて、出力信号にインパルス信号を加えた信号に基づいて、測定系の周波数伝達関数を導出する。
【００１７】
この構成によると、周波数伝達関数を用いて超音波受波器から出力される出力信号が所望波形を有するような入力信号の波形を導出しているので、超音波受波器から出力される出力信号の波形を時間幅の狭いものとすることが可能である。また、導出された伝達関数が反射波の影響を組み込んだものとなるので、導出された波形を有する入力信号を超音波送波器に供給したときに超音波受波器から出力される出力信号が反射波成分をほとんど含まなくなる。したがって、高精度の測定結果を得ることが可能となる。
【００１８】
また、出力信号とインパルス信号とから成る信号に基づいて測定系の周波数伝達関数を導出するため、ノイズの影響をできる限り排除することができる。これにより、さらに高精度の測定結果を得ることができる。
【００１９】
本発明の超音波測定装置において、前記伝達関数導出手段が、前記超音波送波器への入力信号をフーリエ変換するための第１の変換手段と、前記超音波受波器が超音波を受波したときの出力信号にインパルス信号を付加するためのインパルス信号付加手段と、前記インパルス信号付加手段によってインパルス信号が付加された出力信号をフーリエ変換するための第２の変換手段と、前記第１の変換手段によって得られた変換入力信号と、前記第２の変換手段によって得られた変換出力信号との比を求めるための求比手段とを備えていてよい。これにより、伝達関数を迅速に求めることが可能である。
【００２０】
【発明の実施の形態】
以下、本発明の好適な実施の形態について、図面を参照しつつ説明する。
【００２１】
図１は、本発明の第１の実施の形態による超音波骨密度測定装置を示す一部断面正面図である。図２は、図１に示す超音波骨密度測定装置のII-II線に沿った断面図である。図３は、図１に示す超音波骨密度測定装置のIII-III線に沿った断面図である。
【００２２】
図１〜図３に示す超音波骨密度測定装置１において、測定体２は、これと対向配置された測定体３との間隔を可変するために移動自在にされた円筒状の可動タンク７を有している。測定体２、３内には、超音波送受波器４、５がそれぞれ配置されている。超音波送受波器４、５は、超音波の送波及び受波を共に行うことができる機器である。
【００２３】
可動タンク７は、測定体３側に被検体（超音波骨密度測定装置１は被検体として踵を用いるので、以下において、被検体のことを単に踵ということがある）を押し付けるスタンドオフ８が先端に圧入されることで、内部が密閉状態にされている。密閉された測定体３の内部空間には、超音波送受波器４、５から発生される超音波信号を透過し易い超音波透過物質となる超音波透過性流体（例えば、水）が注入されている。
【００２４】
スタンドオフ８は、測定体３に向かって段々に縮径突出する断面台形を呈しており、可動タンク７の軸方向に直交する２つの反射平面８ａ、８ｂを有している。また、スタンドオフ８は、超音波信号を透過させ易く且つ超音波透過性流体と異なる音響インピーダンスを有するアクリル、エポキシ、ウレタンおよびシリコン等の各種樹脂を成形してなるものである。特に、被検体への押し付けに際して反射平面８ａ、８ｂの変形が少ない硬質材料であるアクリル樹脂は、スタンドオフ８の材料として好ましい。なお、後述するスタンドオフ２５もスタンドオフ８と同じ材料からなる。
【００２５】
測定体２の可動タンク７は、対となる２つのタンク支え１１ａ、１１ｂの間に形成されたＶ溝に支持されていると共に、踵が載せられる基台である診断台９に対して駆動機構１０によって移動させることが可能である。駆動機構１０は、可動タンク７に固設されてその軸方向に延びたラック１０ａに噛み合うピニオン１０ｂを有している。ピニオン１０ｂは、診断台９に立設された支持部材１３のオイルレスブッシュによって回転自在に支持されたダイヤルシャフト１２ａに挿入されることによって固定されている。ダイヤルシャフト１２ａのピニオン１０ｂよりも先端側には、取っ手１２ａが取り付けられたダイアル１２が取り付けられている。そのため、ダイアル１２を正逆回転させると、ピニオン１０ｂの回転がラック１０ａに伝達されて、可動タンク７がスタンドオフ８と共に測定体３に近接又は離隔する方向に直線移動する。このようにして、測定体２と測定体３との間隔を変えることができる。また、駆動機構１０には、ピニオン１０ｂの回転を規制する固定ピン（図示せず）が含まれている。この固定ピンによるピニオン１０ｂの回転規制で、２つの測定体２、３の間隔を一定範囲内に保持できるようになっている。
【００２６】
可動タンク７は、図２及び図３に示すように、その軸方向に並行して延びるタンクガイド１５を有している。タンクガイド１５は、診察台９に沿って配置されたガイドシャフト１６と、ガイドシャフト１６に沿ってＵ字型に曲げられたガイド蓋１８とを主要部品としている。ガイドシャフト１６は、その両端側において診断台９に立設された一対の固定台２２ａ、２２ｂに軸支されている。また、ガイド蓋１８は、可動タンク７の外周に固設されたガイド板１７と一体化されることによって、ガイド板１７と共にガイド筒体を構成している。このガイド筒体は、ガイド板１７とガイド蓋１８とで区画されたガイド孔２３を有している。ガイド孔２３には、ガイドシャフト１６が貫通している。そして、ガイド筒体は、ガイド蓋１８の両端の軸受部でガイドシャフト１６に対して軸方向摺動自在に支持されている。
【００２７】
ガイド孔２３は、ガイドシャフト１６の軸方向に延びると共に、部材１８、１７の側面に開口を有している。ガイド孔２３の軸方向両端部は、シール部材１９ａ、１９ｂ及び密閉板２０ａ、２０ｂがそれぞれこの順にガイドシャフト１６に摺動自在に外嵌されることによって密閉されている。可動タンク７が駆動機構１０によって移動させられる際、シール部材１９ａ、１９ｂと密閉板２０ａ、２０ｂとがガイドシャフト１６に対して摺動するので、可動タンク７は安定して案内される。
【００２８】
ガイド孔２３は、ガイドシャフト１６の軸方向に延びるようにガイド板１７に形成された連絡長孔２１を通して可動タンク７内に連通している。したがって、可動タンク７内の超音波透過性液体は、ガイド孔２３内にも充満する。シール部材１９ａ、１９ｂと密閉板２０ａ、２０ｂとでガイド孔２３が軸方向両端部において密閉されているために、ガイド孔２３内に充満した超音波透過性液体が外部に漏れることがない。
【００２９】
測定体３は、可動タンク７のスタンドオフ８に対向して配されたスタンドオフ２５を有している。スタンドオフ２５は、超音波送受波器４、５で発生した超音波信号を透過し易い超音波透過性固体からなる。スタンドオフ２５は、可動タンク７のスタンドオフ８に向かって徐々に縮径した断面台形の突出形状を呈しており、スタンドオフ８の反射面８ａに平行な反射面２５ａを有している。
【００３０】
測定体３は、スタンドオフ２５の軸心が可動タンク７のスタンドオフ８の軸心と一致し且つ反射面２５ａが反射面８ａと平行となるように、診断台９に立設された固定台２８に固設されている。したがって、駆動機構１０によりスタンドオフ８をスタンドオフ２５に近づく方向に直線移動させることにより、２つのスタンドオフ８、２５で踵を挟み込むことができるので、踵を診断台９上に固定することができる。踵を固定するために必要なスタンドオフ８の移動量は、踵の幅寸法に依存する。
【００３１】
上記のように２つのスタンドオフ８、２５が先端側に向かって縮径した断面台形形状となっているので、反射面８ａ、２５ａの面積は比較的小さい。そのため、表面に凹凸がありその凹凸形状が一定でない踵等が被検体であっても、スタンドオフ８、２５の反射面８ａ、２５ａの全面を被検体に密着させ易くなっている。
【００３２】
超音波送受波器４、５としては、通常、超音波信号の発生および検出を１台で行うことができる超音波トランスデューサと呼ばれる装置が用いられる。超音波送受波器４、５は、相互間の距離Ｌ０が固定された状態で測定体２、３内に配置されている。
【００３３】
超音波送受波器４は、スタンドオフ８の反射平面８ａに向かって超音波信号を送信及び受信可能に透過流体が充満した可動タンク７中に配置されると共に、連結部材２７の先端に取り付けられている。連結部材２７は、その基端がガイドシャフト１６に固定されており、連絡長孔２１から可動タンク７内に突出するように延びている。連絡長孔２１がガイドシャフト１６の軸方向に延びるようにガイド板１７に形成されているために、超音波送受波器４をガイドシャフト１６に固定した状態でスタンドオフ８及び可動タンク７を直線移動させることができるようになっている。
【００３４】
超音波送受波器５は、スタンドオフ２５の反射平面２５ａに向かって超音波信号を送信及び受信可能に、気密状態のスタンドオフ２５内に固定されている。上記のように超音波送受波器４、５を配置することにより、一方の超音波送受波器４、５から発せられ、超音波透過性流体およびスタンドオフ８、２５を透過した超音波信号を他方の超音波送受波器４、５で検出することができる。また、一方の超音波送受波器４、５から発せられ、いずれかの反射面８ａ、８ｂ、２５ａで反射した超音波信号を当該一方の超音波送受波器４、５で受信することも可能となる。
【００３５】
次に、超音波骨密度測定装置１のシステム構成について説明する。超音波骨密度測定装置１のシステム構成を示すブロック図である図４に示すように、超音波骨密度測定装置１において、２つの超音波送受波器４、５には、増幅器３１、Ａ／Ｄコンバータ３２、演算器３３、入力波形発生器３４、及び、増幅器３５が接続されている。
【００３６】
なお、以下の説明では、超音波送受波器５が超音波を送波し、超音波送受波器４が超音波を受波するとして説明するが、これとは逆に、超音波送受波器４が超音波を送波し、超音波送受波器５が超音波を受波するようにしてもよい。また、２つの超音波送受波器４、５のいずれか一方だけが超音波の送波及び受波を両方行うようにしてもよい。
【００３７】
増幅器３１は、超音波送受波器４に接続されている。増幅器３１は、超音波を受波した超音波送受波器４の出力信号の電圧を増幅する。増幅器３１には、Ａ／Ｄコンバータ３２が接続されている。Ａ／Ｄコンバータ３２は、増幅器３１の出力信号を、演算器３３で演算処理を実行できるようにデジタル化する。
【００３８】
演算器３３は、以下に説明するような各種の演算処理を実行する。演算器３３は、伝達関数導出部３３１と、入力波形導出部３３６とを含んでいる。伝達関数導出部３３１は、所望出力波形記憶部３３２と、インパルス信号記憶部３３３と、変換部３３４と、求比部３３５とを含んでいる。入力波形導出部３３６は、変換部３３７と、求比部３３８とを含んでいる。
【００３９】
所望出力波形記憶部３３２は、超音波送受波器４で受信されるのに望ましい波形を有する出力信号ｙ（ｔ）（ｔ：時間）に関する波形情報を記憶している。本実施の形態において、所望出力波形記憶部３３２で記憶されている所望出力信号ｙ（ｔ）の波形情報は、周波数が後述する伝達関数の中心周波数の２倍内である１周期分の正弦波（以下、「ｓｉｎ１波」と称する）である。
【００４０】
インパルス信号記憶部３３３は、超音波送受波器５に供給されるインパルス信号ｉ（ｔ）に関する波形情報を記憶している。
【００４１】
変換部３３４は、インパルス信号記憶部３３３に記憶されているインパルス信号ｉ（ｔ）、及び、超音波送受波器４から出力されるインパルス応答信号ｏ（ｔ）をそれぞれフーリエ変換する。以下の説明において、インパルス信号ｉ（ｔ）及びインパルス応答信号ｏ（ｔ）を変換部３３４によってフーリエ変換した信号を、それぞれ、Ｉ（ω）及びＯ（ω）（ω：角周波数）と記載する。
【００４２】
求比部３３５は、変換部３３４で得られた信号Ｏ（ω）とＩ（ω）との比を求める。この比が、超音波骨密度測定装置１の測定系における伝達関数Ｈ（ω）である。
【００４３】
一方、入力波形導出部３３６の変換部３３７は、所望出力波形記憶部３３２に記憶されているｓｉｎ１波の波形を有する出力信号ｙ（ｔ）をフーリエ変換する。以下の説明において、出力信号ｙ（ｔ）を変換部３３７によってフーリエ変換した信号を、Ｙ（ω）と記載する。
【００４４】
求比部３３８は、変換部３３７で得られた信号Ｙ（ω）と、求比部３３５で得られた伝達関数Ｈ（ω）との比を求める。以下の説明において、この比をＸ（ω）と記載する。求比部３３８で求められたＸ（ω）は変換部３３７に与えられて逆フーリエ変換される。これによって、超音波送受波器５に入力すべき入力信号ｘ（ｔ）に関する波形情報が得られる。得られた波形情報は、演算部３３内に記憶される。
【００４５】
入力波形発生器３４は、入力波形導出部３３６で導出された波形を有する入力信号ｘ（ｔ）を発生させる。発生した入力信号ｘ（ｔ）は、増幅器３５により電圧が増幅されて、超音波送受波器５に供給される。
【００４６】
次に、超音波骨密度測定装置１の全体の処理手順について図５及び図６を参照しつつ説明する。図５は、超音波骨密度測定装置１の全体の処理手順を表したフローチャートである。図６は、入力すべき入力信号ｘ（ｔ）を導出する一連の処理を図示したものである。なお、図６内に描かれた信号の形状はいずれも一例として示すものに過ぎない。
【００４７】
まず、ステップＳ２０１において、スタンドオフ８、２５の反射面８ａ、２５ａを接触させる。スタンドオフ８及び２５の接触時には超音波の透過のため、超音波ゼリーを塗布する。その後、ステップＳ２０３に移行する。ステップＳ２０３では、インパルス信号記憶部３３３に記憶されているインパルス信号ｉ（ｔ）（図６のブロックＢ１１内の（ａ）参照）を超音波送受波器５に与え、超音波送受波器５を励振させる。次にステップＳ２０５に移行して、超音波送受波器５が発生した超音波を超音波送受波器４が受信する。すると、超音波送受波器４は、インパルス応答信号ｏ（ｔ）（図６のブロックＢ１２内の（ａ）参照）を出力する。インパルス応答信号ｏ（ｔ）は、増幅器３１、Ａ／Ｄコンバータ３２を経て演算器３３に与えられる。その後、ステップＳ２０７に移行する。
【００４８】
ステップＳ２０７では、変換部３３４により、インパルス信号ｉ（ｔ）をフーリエ変換する。これにより、信号Ｉ（ω）（図６のブロックＢ１１内の（ｂ）参照）が求められる。その後、ステップＳ２０９に移行する。ステップＳ２０９では、変換部３３４により、ステップＳ２０５で超音波送受波器４が発生したインパルス応答信号ｏ（ｔ）をフーリエ変換する。これにより、信号Ｏ（ω）（図６のブロックＢ１２内の（ｂ）参照）が求められる。その後、ステップＳ２１１に移行する。
【００４９】
ステップＳ２１１では、求比部３３５により、ステップＳ２０７及びステップＳ２０９において求められた信号Ｉ（ω）及び信号Ｏ（ω）から、超音波骨密度測定装置１の測定系の周波数伝達関数Ｈ（ω）が計算される（図６のブロックＢ１３参照）。ここで伝達関数Ｈ（ω）は次式で計算できる。計算後、ステップＳ２１３に移行する。
【００５０】
【数１】

Figure 0004312494
【００５１】
ステップＳ２１３では、入力波形導出部３３６により、超音波送受波器５に励振すべき入力信号ｘ（ｔ）が計算される。具体的には、まず、変換部３３７により、所望出力波形記憶部３３２に記憶されているｓｉｎ１波の波形を有する所望出力信号ｙ（ｔ）（図６のブロックＢ１４内の（ａ）参照）をフーリエ変換し、信号Ｙ（ω）（図６のブロックＢ１４内の（ｂ）参照）を得る。さらに、求比部３３８により、信号Ｙ（ω）とステップＳ２１１において導出された伝達関数Ｈ（ω）との比である信号Ｘ（ω）（図６のＢ１５参照）が求められる。信号Ｘ（ω）は次式のように表される。
【００５２】
【数２】
Figure 0004312494
【００５３】
そして、変換部３３７が信号Ｘ（ω）を逆フーリエ変換する。これによって、超音波送受波器５に入力すべき入力信号ｘ（ｔ）に関する波形情報（図６のブロックＢ１６参照）が得られ、ステップＳ２１５に移行する。
【００５４】
ステップＳ２１５では、ステップＳ２０１において接触させた２つのスタンドオフ８、２５の反射面８ａ、２５ａを再び分離させる。分離させたスタンドオフ８、２５の間に被検体である踵を挿入して、スタンドオフ８、２５の反射面８ａ、２５ａを踵に密着させる。スタンドオフ８及び２５と踵との密着時には超音波の透過のため、超音波ゼリーを塗布する。その後、ステップＳ２１７において、ステップＳ２１３において求められた入力信号ｘ（ｔ）を入力波形発生器３４により発生させる。さらに、増幅器３５により入力信号ｘ（ｔ）の電圧を増幅させて、超音波送受波器５に供給する。
【００５５】
その後、ステップＳ２１９に移行して、超音波送受波器４により、踵を透過した超音波を受信する。その後、ステップＳ２２１で、超音波送受波器４から出力された出力信号に基づいて踵の骨密度を算出する。算出された踵の骨密度は、図示しないディスプレイに表示される。
【００５６】
図７は、超音波骨密度測定装置１での各処理段階で表れる各信号の一例を描いた図である。図７に示すように、インパルス信号（図７のブロックＢ２１内の（ａ）参照）及びインパルス応答信号（図７のブロックＢ２１内の（ｂ）参照）を用いて導出された伝達関数は、測定体２、３及びスタンドオフ８、２５で発生する反射波の影響を組み込んだものとなっているので、この伝達関数を利用して求められた入力すべき入力信号（図７のブロックＢ２２参照）を超音波送受波器５に供給したときに超音波送受波器４から出力された出力信号（図７のブロックＢ２３参照）は、反射波の影響をほとんど受けず、ｓｉｎ１波に近いものとなる。また、超音波送受波器４からの出力信号の時間幅は、高価なコンポジット振動子を用いなくても十分に狭いものとなる。したがって、安価な装置構成によって時間分解能に優れた高精度な測定結果を得ることが可能となる。
【００５７】
また、本実施の形態において、測定体２、３の形状は、前記方法により測定体がいかなる形状であっても測定体内部での反射波の影響を受けないため、自由に被検体である踵の測定に最も適した形状に決定することができる。
【００５８】
また、上述したように伝達関数導出部３３１がフーリエ変換及び求比を行うことで伝達関数を導出しているので、伝達関数を迅速に求めることが可能となっている。
【００５９】
次に、本発明に係る第２の実施の形態による超音波骨密度測定装置について、図面を参照しつつ説明する。本実施の形態の超音波骨密度測定装置は、伝達関数の導出手順において第１の実施の形態の超音波骨密度測定装置と相違している。以下、その相違点を中心に説明する。なお、以下の説明において、第１の実施の形態と同様の部材には同じ符号を付してその説明を省略する。
【００６０】
図８は、本実施の形態による超音波骨密度測定装置のシステム構成を示すブロック図である。演算器３６は、伝達関数導出部３６１、及び、入力波形導出部３３６を含んでいる。伝達関数導出部３６１は、所望出力波形記憶部３３２、インパルス信号記憶部３６３、インパルス信号付加部３６４、変換部３６５、及び、求比部３６６を含んでいる。
【００６１】
インパルス信号記憶部３６３は、超音波送受波器５に供給されるインパルス信号ｉ（ｔ）に関する波形情報、及び、インパルス信号ｉ（ｔ）を超音波送受波器５に供給したときの超音波送受波器４からの出力信号ｏ（ｔ）に付加するためのインパルス信号δ（ｔ）に関する波形情報を記憶している。
【００６２】
インパルス信号付加部３６４は、超音波送受波器４から出力されるインパルス応答信号ｏ（ｔ）にインパルス信号記憶部３６３に記憶されているインパルス信号δ（ｔ）を付加する。インパルス信号付加部３６４によってインパルス信号δ（ｔ）が付加されたインパルス応答信号ｏ（ｔ）を、以下の説明において信号ｋ（ｔ）と記載する。
【００６３】
変換部３６５は、インパルス信号記憶部３３３に記憶されているインパルス信号ｉ（ｔ）、及び、インパルス信号付加部３６４でのインパルス信号δ（ｔ）の付加によって得られた信号ｋ（ｔ）をそれぞれフーリエ変換する。以下の説明において、インパルス信号ｉ（ｔ）及び信号ｋ（ｔ）を変換部３６５によってフーリエ変換した信号を、それぞれ、Ｉ（ω）及びＫ（ω）と記載する。
【００６４】
求比部３６６は、変換部３６５で得られた信号Ｋ（ω）とＩ（ω）との比を求める。この比が、本実施の形態による超音波骨密度測定装置の測定系における伝達関数Ｇ（ω）である。
【００６５】
入力波形導出部３３６での処理は、伝達関数Ｈ（ω）がＧ（ω）と表記される以外は同じであるので、説明を省略する。
【００６６】
次に、本実施の形態に係る超音波骨密度測定装置の全体の処理手順について図９のフローチャートを参照しつつ説明する。図９において、ステップＳ３０７までの各処理手順は、第１の実施の形態のステップＳ２０７までの各処理手順と同じであるため、その説明を省略する。
【００６７】
ステップＳ３０９では、インパルス信号付加部３６４により、ステップＳ３０５で超音波送受波器４が発生したインパルス応答信号ｏ（ｔ）に、インパルス信号記憶部３６３に記憶されているインパルス信号δ（ｔ）を付加する。このとき、インパルス信号δ（ｔ）は、インパルス応答信号ｏ（ｔ）の平坦部分に重ねられる。これにより、信号ｋ（ｔ）が得られる。さらにステップＳ３０９では、変換部３６５により、信号ｋ（ｔ）をフーリエ変換する。これにより、信号Ｋ（ω）が求められる。その後、ステップＳ３１１に移行する。
【００６８】
ステップＳ３１１では、求比部３６５により、ステップＳ３０７及びステップＳ３０９において求められた信号Ｉ（ω）及び信号Ｋ（ω）から、超音波骨密度測定装置の測定系の周波数伝達関数Ｇ（ω）が計算される。ここで伝達関数Ｇ（ω）は次式で計算できる。計算後、ステップＳ３１３に移行する。
【００６９】
【数３】
Figure 0004312494
【００７０】
ステップＳ３１３以降の各処理手順は、第１の実施の形態のステップＳ２１３〜Ｓ２２１の各処理手順と同様である。
【００７１】
ここで、本実施の形態による超音波骨密度測定装置で導出される伝達関数Ｇ（ω）について、第１の実施の形態における伝達関数Ｈ（ω）と比較しつつ説明する。図１０（ａ）、（ｂ）は、それぞれ、第１の実施の形態におけるインパルス応答信号ｏ（ｔ）及び伝達関数Ｈ（ω）の信号波形図である。図１１（ａ）、（ｂ）は、それぞれ、本実施の形態における信号ｋ（ｔ）（インパルス応答信号ｏ（ｔ）にインパルス信号δ（ｔ）を付加した信号）及び伝達関数Ｇ（ω）の信号波形図である。
【００７２】
図１０（ｂ）に示されているように、伝達関数Ｈ（ω）の高周波域及び低周波域には、不要成分であるノイズが重畳している。このようにノイズが重畳した伝達関数Ｈ（ω）に基づいて入力すべき入力信号ｘ（ｔ）を計算すると、求められた入力信号ｘ（ｔ）の全体に、伝達関数Ｈ（ω）の高周波域成分及び低周波域成分が重畳されてしまう。その結果、この入力信号ｘ（ｔ）を超音波送受波器５に与えたときに得られる出力信号がｓｉｎ１波のような波形から大きくずれて、高精度な測定結果が得られなくなる。
【００７３】
一方、図１１（ｂ）に示されている伝達関数Ｇ（ω）は、インパルス信号δ（ｔ）を付加することで、伝達関数Ｇ（ω）の最大値Ｇ1（周波数ｆcに対応する）よりも小さい値Ｇ2以下の部分が埋め込まれて（カットされて）一定値となった、高周波域及び低周波域にノイズを有さない波形となっている。したがって、伝達関数Ｇ（ω）に基づいて算出された入力信号ｚ（ｔ）には、伝達関数Ｇ（ω）の高周波域及び低周波域のノイズ成分が重畳されることがない。そのため、ｚ（ｔ）の波形形状が安定し、実回路による送信が簡易になる。
【００７４】
なお、図１１（ｂ）に示すような一定値以下の領域が埋め込まれた伝達関数Ｇ（ω）は、フーリエ変換後の複素演算によっても得ることができる。しかしながら、この複素演算は注意深く行う必要がある上に、複素演算をする必要上計算量も多く長時間を要する。本実施の形態では、インパルス信号を付加するという簡単な処理によって、複雑なアルゴリズムを必要とせず、上記のような伝達関数Ｇ（ω）を簡易に得ることができる点で利益が大きい。
【００７５】
また、本実施の形態によると、上述した第１の実施の形態と同様の利益、すなわち、反射波の影響をほとんど受けず且つ超音波送受波器４からの出力信号の時間幅が十分に狭いものとなるために安価な装置構成によって時間分解能に優れた高精度な測定結果を得ることが可能であること、及び、伝達関数を迅速に求めることが可能であることが得られる。
【００７６】
以上、本発明の好適な実施の形態について説明したが、本発明は上述の実施の形態に限られるものではなく、特許請求の範囲に記載した限りにおいて様々な設計変更が可能なものである。例えば、上述の実施の形態では、伝達関数を導出する際に、スタンドオフ８、２５の反射面８ａ、２５ａを接触させて導出しているが、スタンドオフ８、２５の反射面８ａ、２５ａを接触させずに超音波を被検体に透過させた状態で伝達関数を導出してもよい。
【００７７】
また、上述の実施の形態では、２つの超音波送受波器を用いているが、どちらか一方を超音波送波器とし、他方を超音波受波器としてもよい。
【００７８】
上述の実施の形態では、所望出力信号ｙ（ｔ）の波形をｓｉｎ１波としているが、それ以外の波形としてもよい。例えば、矩形波などでもよい。上述の実施の形態では、演算を実行する演算器を備えているが、演算器を備えず、演算をパーソナルコンピュータ上で実行してもよい。
【００７９】
上述の実施の形態では、インパルス信号のフーリエ変換Ｉ（ω）で除しているが、十分細いインパルスであればＩ（ω）＝１と見なすこと可能である。この場合、Ｈ（ω）＝Ｏ（ω）となる。
【００８０】
上述の実施の形態では、２つの超音波送受波器を用いて超音波を送受しているが、１つの超音波送受波器のみで超音波を送受してもよい。具体的には、超音波送受波器が送信した超音波の反射波を、同じ超音波送受波器で受信して、伝達関数を導出してもよい。
【００８１】
また、第２の実施の形態による超音波骨密度測定装置は、被検体が骨に限定されるものではなく、超音波測定装置として用いることも可能である。その場合、超音波測定装置は基台や測定体といった部材を備えていなくてもよい。
【００８２】
【発明の効果】
以上説明したように、本発明の超音波骨密度測定装置によると、送波される超音波の時間幅延長を抑制し且つ反射波の影響を少なくすることで高精度の測定結果を得ることが可能である。また、本発明の超音波測定装置によると、ノイズの影響を受けにくい高精度の測定結果を得ることが可能である。
【図面の簡単な説明】
【図１】本発明の第１の実施の形態による超音波骨密度測定装置の全体の構成を示す側面図である。
【図２】図１に示す超音波骨密度測定装置のII-II線に沿った断面図である。
【図３】図１に示す超音波骨密度測定装置のIII-III線に沿った断面図である。
【図４】図１に示す超音波骨密度測定装置のシステム構成を示すブロック図である。
【図５】図１に示す超音波骨密度測定装置による骨密度測定時の処理手順を表したフローチャートである。
【図６】図１に示す超音波骨密度測定装置に入力すべき波形ｘ（ｔ）を導出する一連の処理図である。
【図７】図１に示す超音波骨密度測定装置の出力波形を得るまでの概略図である。
【図８】本発明の第２の実施の形態による超音波骨密度測定装置のシステム構成を示すブロック図である。
【図９】図８に示す超音波骨密度測定装置による骨密度測定時の処理手順を表したフローチャートである。
【図１０】図８に示す超音波骨密度測定装置の出力信号にインパルス信号を付加せず伝達関数を導出する場合の信号波形図である
【図１１】図８に示す超音波骨密度測定装置の出力信号にインパルス信号を付加して伝達関数を導出する場合の信号波形図である。
【符号の説明】
１ 超音波骨密度測定装置
２、３ 測定体
４、５ 超音波送受波器
８、２５ スタンドオフ
９ 診断台（基台）
３３、３６ 演算器
３３１、３６１ 伝達関数導出部（伝達関数導出手段）
３３２ 所望出力波形記憶部
３３３、３６３ インパルス信号記憶部
３３４、３６５ 変換部（第１の変換手段、第２の変換手段）
３３５、３６６ 求比部（求比手段）
３３６ 入力波形導出部（入力波形導出手段）
３６４ インパルス信号付加部（インパルス信号付加手段）[0001]
BACKGROUND OF THE INVENTION
The present invention relates to an ultrasonic bone density measuring device, an ultrasonic measuring device, an ultrasonic bone density measuring method, and an ultrasonic measuring method used for diagnosis of osteoporosis and the like.
[0002]
[Prior art]
An ultrasonic bone density measuring device is used for diagnosis of osteoporosis, which is a symptom in which the density of bone tissue decreases. Ultrasound propagates through the subject at a speed corresponding to the amount of bone mineral contained in the subject. Therefore, if an ultrasonic wave propagation speed (SOS: Speed Of Sound) in the subject is measured, an amount corresponding to the bone mineral content of the subject's eyelid can be measured. The ultrasonic bone density measuring device uses this principle, and measures the propagation speed of ultrasonic waves in the bone of the subject, and the bone characteristics such as bone mineral density from the measured ultrasonic propagation speed. Is a device capable of quantitative measurement. In order to improve the measurement accuracy of the ultrasonic propagation velocity, the ultrasonic inspection by the ultrasonic bone density measuring apparatus is often performed using a rib with many cancellous bones as a subject.
[0003]
The ultrasonic bone density measuring apparatus described in Patent Document 1 includes a base for placing a subject and a pair of measuring bodies that sandwich the subject from both sides and are in close contact with the subject. An ultrasonic transmitter having an ultrasonic transducer is arranged in one measurement body, and an ultrasonic receiver having an ultrasonic transducer (substantially an ultrasonic transmitter and Are the same). The ultrasonic wave transmitted from the ultrasonic wave transmitter passes through the subject and is received by the ultrasonic wave receiver. Further, the ultrasonic wave reflected at the interface between the measurement object and the subject is measured, and the ultrasonic propagation velocity of the subject is obtained from the time difference from the passed ultrasonic wave.
[0004]
[Patent Document 1]
JP-A-10-43180
[0005]
In the ultrasonic bone density measuring apparatus as described in Patent Document 1, a single pulse signal is given to the ultrasonic transmitter when measuring the ultrasonic propagation velocity. At this time, since the frequency band of the ultrasonic transducer is narrow, the ultrasonic wave output from the ultrasonic transmitter has a waveform that spreads along the time axis (wide time width). Therefore, the time resolution of the ultrasonic bone density measuring device is deteriorated. If a composite vibrator is used as the ultrasonic vibrator, it is possible to shorten the time width of the ultrasonic wave output from the ultrasonic transmitter, but the time width shortening effect is not sufficient, and the composite vibration The child itself is expensive and not practical for use in a commercial device. In addition, the ultrasonic wave output from the ultrasonic transmitter partly reflects on the surface other than the boundary surface between the measurement object and the object sandwiching the object, and arrives at the wave receiver with a time delay. In such a case, the measurement accuracy deteriorates. This reflected wave (hereinafter referred to as “reflected wave”) depends on the shape of the measurement object, and even when receiving an ultrasonic wave that has passed through the subject, the ultrasonic wave reflected at the interface between the measurement object and the subject is received. Can occur in either case.
[0006]
Moreover, in an ultrasonic measurement apparatus including an ultrasonic bone density measurement apparatus, it is desired to obtain a highly accurate measurement result by eliminating the influence of noise as much as possible.
[0007]
Accordingly, an object of the present invention is to provide an ultrasonic bone density measuring apparatus capable of obtaining a highly accurate measurement result by suppressing the extension of the time width of transmitted ultrasonic waves and reducing the influence of reflected waves. And an ultrasonic bone density measuring method.
[0008]
Another object of the present invention is to provide an ultrasonic measurement apparatus and an ultrasonic measurement method capable of obtaining a highly accurate measurement result that is not easily affected by noise.
[0009]
[Means for Solving the Problems]
  The ultrasonic bone density measuring apparatus of the present invention includes a base for placing a subject, a pair of measuring bodies made of an ultrasonically transparent solid material, at least one of which is movably opposed to the other, An ultrasonic transmitter provided in one of the pair of measuring bodies, an ultrasonic receiver provided in either of the pair of measuring bodies, and the ultrasonic wave from the ultrasonic transmitter Based on the transfer function deriving means for deriving the frequency transfer function of the measurement system leading to the receiver, the desired waveform of the output signal, and the frequency transfer function derived by the transfer function deriving means, the ultrasonic reception An input waveform deriving unit for deriving a waveform of an input signal supplied to the ultrasonic wave transmitter when the wave generator outputs an output signal having the desired waveform, and a waveform derived by the input waveform deriving unit The input signal having the ultrasonic wave Have an input waveform supply means for supplying to the filterThe
[0010]
  The ultrasonic measurement method of the present invention includes an ultrasonic transmitter provided on one of a pair of measurement bodies made of an ultrasonically transparent solid material, at least one of which is movably opposed to the other. Derived in a transfer function deriving step for deriving a frequency transfer function of a measurement system reaching an ultrasonic wave receiver provided in one of the pair of measuring bodies, a desired waveform of an output signal, and the transfer function deriving step An input waveform deriving step for deriving a waveform of an input signal supplied to the ultrasonic wave transmitter when the ultrasonic wave receiver outputs an output signal having the desired waveform based on the frequency transfer function obtained And withThe
[0011]
  The ultrasonic measurement method of the present invention includes at least one of a pair of measurement bodies made of an ultrasonically transparent solid material, at least one of which is movably opposed to the other direction and provided with a standoff at the tip. The frequency transfer function of the measurement system from the ultrasonic transmitter provided to the ultrasonic receiver provided to one of the pair of measurement bodies in a state where the two standoffs are in contact with each other When the ultrasonic receiver outputs an output signal having the desired waveform based on the transfer function deriving step to be derived, the desired waveform of the output signal, and the frequency transfer function derived in the transfer function deriving step. An input waveform deriving step for deriving a waveform of an input signal supplied to the ultrasonic wave transmitter, a step of separating the two standoffs, and two spaced apart standoffs In the state in which the wrinkles arranged in the two standoffs are sandwiched between the two standoffs and in the state in which the wrinkles are sandwiched between the two standoffs, the ultrasonic wave transmission of the input signal having the waveform obtained in the input waveform deriving step And a step of detecting an output signal from the ultrasonic receiver.TheIn this specification, “stand-off” refers to a portion that comes into contact with a patient.
[0012]
According to this configuration, the output signal output from the ultrasonic receiver is derived using the frequency transfer function so that the output signal output from the ultrasonic receiver has a desired waveform. It is possible to make the waveform of the signal narrow. In addition, since the derived transfer function incorporates the influence of the reflected wave, the output signal output from the ultrasonic receiver when an input signal having the derived waveform is supplied to the ultrasonic transmitter However, almost no reflected wave component is included. Therefore, a highly accurate measurement result can be obtained.
[0013]
  In the present invention, the shape of each measurement object is desired to include the influence of the reflected wave from the ultrasonic transmitter to the ultrasonic receiver of the ultrasonic wave transmitted from the ultrasonic transmitter in the transfer function. It is preferable to determine the shape suitable for the measurement of the subject so that the waveform output signal can be obtained.Yes.According to this, since any shape of the measurement body is not affected by the reflected wave inside the measurement body, it is possible to freely determine the shape most suitable for the measurement of the subject.
[0014]
  The ultrasonic measurement apparatus of the present invention includes an ultrasonic transmitter, an ultrasonic receiver, and a transmission for deriving a frequency transfer function of a measurement system from the ultrasonic transmitter to the ultrasonic receiver. Based on the function derivation means, the desired waveform of the output signal, and the frequency transfer function derived by the transfer function derivation means, the ultrasonic receiver outputs the output signal having the desired waveform. Input waveform deriving means for deriving a waveform of an input signal supplied to the acoustic wave transmitter, and an input waveform for supplying an input signal having a waveform derived by the input waveform deriving means to the ultrasonic wave transmitter Supply means. The transfer function deriving means derives a frequency transfer function of the measurement system based on a signal obtained by adding an impulse signal to the output signal from the ultrasonic wave receiver that has received the ultrasonic wave.The
[0015]
  In another aspect, the ultrasonic bone density measuring apparatus according to the present invention reaches an ultrasonic wave transmitter, an ultrasonic wave receiver, and the ultrasonic wave transmitter from the ultrasonic wave transmitter through the object to be inspected. Based on the transfer function deriving means for deriving the frequency transfer function of the measurement system, the desired waveform of the output signal, and the frequency transfer function derived by the transfer function deriving means, the ultrasonic receiver receives the desired An input waveform deriving means for deriving a waveform of an input signal supplied to the ultrasonic transmitter when outputting an output signal having a waveform; and an input signal having a waveform derived by the input waveform deriving means Input waveform supply means for supplying to the ultrasonic transmitter. The transfer function deriving means derives a frequency transfer function of the measurement system based on a signal obtained by adding an impulse signal to the output signal from the ultrasonic wave receiver that has received the ultrasonic wave.The
[0016]
  The ultrasonic measurement method of the present invention includes a transfer function deriving step for deriving a frequency transfer function of a measurement system from an ultrasonic transmitter to the ultrasonic receiver, and the ultrasonic receiver that receives the ultrasonic wave. When the ultrasonic receiver outputs the output signal having the desired waveform based on the desired waveform of the output signal from the frequency transfer function derived in the transfer function deriving step, the ultrasonic wave transmission An input waveform deriving step for deriving a waveform of an input signal supplied to the device. In the transfer function deriving step, the frequency transfer function of the measurement system is derived based on the signal obtained by adding the impulse signal to the output signal.The
[0017]
According to this configuration, the output signal output from the ultrasonic receiver is derived using the frequency transfer function so that the output signal output from the ultrasonic receiver has a desired waveform. It is possible to make the waveform of the signal narrow. In addition, since the derived transfer function incorporates the influence of the reflected wave, the output signal output from the ultrasonic receiver when an input signal having the derived waveform is supplied to the ultrasonic transmitter However, almost no reflected wave component is included. Therefore, a highly accurate measurement result can be obtained.
[0018]
Further, since the frequency transfer function of the measurement system is derived based on the signal composed of the output signal and the impulse signal, the influence of noise can be eliminated as much as possible. Thereby, it is possible to obtain a more accurate measurement result.
[0019]
  In the ultrasonic measurement apparatus of the present invention, the transfer function deriving means includes first conversion means for Fourier transforming an input signal to the ultrasonic transmitter, and the ultrasonic receiver receives ultrasonic waves. Impulse signal adding means for adding an impulse signal to the output signal when waved, second conversion means for Fourier transforming the output signal to which the impulse signal is added by the impulse signal adding means, and the first And a ratio determining means for determining a ratio between the converted input signal obtained by the converting means and the converted output signal obtained by the second converting means.Yes.Thereby, the transfer function can be obtained quickly.
[0020]
DETAILED DESCRIPTION OF THE INVENTION
Hereinafter, preferred embodiments of the present invention will be described with reference to the drawings.
[0021]
FIG. 1 is a partial sectional front view showing an ultrasonic bone density measuring apparatus according to a first embodiment of the present invention. FIG. 2 is a cross-sectional view taken along line II-II of the ultrasonic bone density measuring apparatus shown in FIG. FIG. 3 is a cross-sectional view taken along line III-III of the ultrasonic bone density measuring apparatus shown in FIG.
[0022]
In the ultrasonic bone density measuring apparatus 1 shown in FIGS. 1 to 3, the measuring body 2 includes a cylindrical movable tank 7 that is movable in order to vary the distance between the measuring body 2 and the measuring body 3 that is disposed to face the measuring body 2. Have. Ultrasonic transducers 4 and 5 are arranged in the measuring bodies 2 and 3, respectively. The ultrasonic transducers 4 and 5 are devices that can both transmit and receive ultrasonic waves.
[0023]
The movable tank 7 has a standoff 8 that presses a subject (the ultrasonic bone density measuring apparatus 1 uses a scissors as a subject to the subject). The inside is sealed by being press-fitted into the tip. An ultrasonic transmission fluid (for example, water) that is an ultrasonic transmission material that easily transmits ultrasonic signals generated from the ultrasonic transducers 4 and 5 is injected into the sealed internal space of the measurement body 3. ing.
[0024]
The standoff 8 has a trapezoidal cross section that gradually projects toward the measuring body 3 and has two reflection planes 8 a and 8 b that are orthogonal to the axial direction of the movable tank 7. Further, the standoff 8 is formed by molding various resins such as acrylic, epoxy, urethane, and silicon that easily transmit ultrasonic signals and have an acoustic impedance different from that of the ultrasonically permeable fluid. In particular, an acrylic resin, which is a hard material with little deformation of the reflection planes 8 a and 8 b when pressed against the subject, is preferable as a material for the standoff 8. Note that a standoff 25 described later is also made of the same material as the standoff 8.
[0025]
The movable tank 7 of the measuring body 2 is supported by a V-groove formed between the two tank supports 11a and 11b as a pair, and has a driving mechanism with respect to the diagnostic table 9 which is a base on which a bag is placed. 10 can be moved. The drive mechanism 10 has a pinion 10b fixed to the movable tank 7 and meshing with a rack 10a extending in the axial direction. The pinion 10 b is fixed by being inserted into a dial shaft 12 a that is rotatably supported by an oilless bush of a support member 13 that is erected on the diagnostic table 9. A dial 12 to which a handle 12a is attached is attached to the distal end side of the dial shaft 12a from the pinion 10b. Therefore, when the dial 12 is rotated forward and backward, the rotation of the pinion 10b is transmitted to the rack 10a, and the movable tank 7 moves linearly in the direction of approaching or separating from the measuring body 3 together with the standoff 8. In this way, the distance between the measurement body 2 and the measurement body 3 can be changed. The drive mechanism 10 includes a fixing pin (not shown) that restricts the rotation of the pinion 10b. By restricting the rotation of the pinion 10b by this fixing pin, the interval between the two measuring bodies 2 and 3 can be kept within a certain range.
[0026]
As shown in FIGS. 2 and 3, the movable tank 7 has a tank guide 15 extending in parallel with the axial direction thereof. The tank guide 15 is mainly composed of a guide shaft 16 disposed along the examination table 9 and a guide lid 18 bent into a U shape along the guide shaft 16. The guide shaft 16 is pivotally supported by a pair of fixed bases 22a and 22b provided upright on the diagnostic base 9 at both ends thereof. The guide lid 18 is integrated with a guide plate 17 fixed on the outer periphery of the movable tank 7, thereby forming a guide cylinder together with the guide plate 17. The guide cylinder has a guide hole 23 defined by a guide plate 17 and a guide lid 18. The guide shaft 16 passes through the guide hole 23. The guide cylinder is supported by the bearing portions at both ends of the guide lid 18 so as to be slidable in the axial direction with respect to the guide shaft 16.
[0027]
The guide hole 23 extends in the axial direction of the guide shaft 16 and has openings on the side surfaces of the members 18 and 17. Both end portions in the axial direction of the guide hole 23 are sealed by externally fitting seal members 19a and 19b and sealing plates 20a and 20b to the guide shaft 16 in this order. When the movable tank 7 is moved by the drive mechanism 10, the seal members 19a and 19b and the sealing plates 20a and 20b slide with respect to the guide shaft 16, so that the movable tank 7 is stably guided.
[0028]
The guide hole 23 communicates with the movable tank 7 through a communication long hole 21 formed in the guide plate 17 so as to extend in the axial direction of the guide shaft 16. Therefore, the ultrasonically transmissive liquid in the movable tank 7 is also filled in the guide hole 23. Since the guide holes 23 are sealed at both ends in the axial direction by the seal members 19a and 19b and the sealing plates 20a and 20b, the ultrasonically permeable liquid filled in the guide holes 23 does not leak to the outside.
[0029]
The measuring body 3 has a standoff 25 arranged to face the standoff 8 of the movable tank 7. The standoff 25 is made of an ultrasonically transparent solid that easily transmits ultrasonic signals generated by the ultrasonic transducers 4 and 5. The standoff 25 has a trapezoidal cross-sectional projection shape with a diameter gradually reduced toward the standoff 8 of the movable tank 7, and has a reflection surface 25 a parallel to the reflection surface 8 a of the standoff 8.
[0030]
The measuring body 3 is a fixed base erected on the diagnostic table 9 such that the axis of the standoff 25 coincides with the axis of the standoff 8 of the movable tank 7 and the reflecting surface 25a is parallel to the reflecting surface 8a. 28 is fixed. Therefore, by moving the standoff 8 linearly in the direction approaching the standoff 25 by the drive mechanism 10, the heel can be sandwiched between the two standoffs 8, 25, so that the heel can be fixed on the diagnostic table 9. it can. The amount of movement of the standoff 8 necessary for fixing the heel depends on the width dimension of the heel.
[0031]
As described above, since the two standoffs 8 and 25 have a trapezoidal shape with a diameter reduced toward the tip, the areas of the reflecting surfaces 8a and 25a are relatively small. For this reason, even if the subject has wrinkles or the like having irregularities on the surface and uneven shapes, the reflective surfaces 8a and 25a of the standoffs 8 and 25 are easily brought into close contact with the subject.
[0032]
As the ultrasonic transducers 4 and 5, an apparatus called an ultrasonic transducer capable of generating and detecting an ultrasonic signal with a single unit is usually used. The ultrasonic transducers 4 and 5 are arranged in the measuring bodies 2 and 3 with the distance L0 between them being fixed.
[0033]
The ultrasonic transducer 4 is disposed in the movable tank 7 filled with a permeable fluid so as to transmit and receive an ultrasonic signal toward the reflection plane 8a of the standoff 8, and is attached to the tip of the connecting member 27. ing. The connecting member 27 has a base end fixed to the guide shaft 16 and extends from the communication slot 21 so as to protrude into the movable tank 7. Since the communication long hole 21 is formed in the guide plate 17 so as to extend in the axial direction of the guide shaft 16, the standoff 8 and the movable tank 7 are straightened with the ultrasonic transducer 4 fixed to the guide shaft 16. It can be moved.
[0034]
The ultrasonic transducer 5 is fixed in the airtight standoff 25 so that an ultrasonic signal can be transmitted and received toward the reflection plane 25a of the standoff 25. By arranging the ultrasonic transducers 4 and 5 as described above, an ultrasonic signal emitted from one ultrasonic transducer 4 and 5 and transmitted through the ultrasonically permeable fluid and the standoffs 8 and 25 is obtained. It can be detected by the other ultrasonic transducers 4 and 5. It is also possible to receive an ultrasonic signal emitted from one of the ultrasonic transducers 4 and 5 and reflected by any one of the reflecting surfaces 8a, 8b and 25a by the one ultrasonic transducer 4 and 5. It becomes.
[0035]
Next, the system configuration of the ultrasonic bone density measuring apparatus 1 will be described. As shown in FIG. 4, which is a block diagram showing the system configuration of the ultrasonic bone density measuring apparatus 1, in the ultrasonic bone density measuring apparatus 1, the two ultrasonic transducers 4 and 5 include an amplifier 31, A / A D converter 32, a calculator 33, an input waveform generator 34, and an amplifier 35 are connected.
[0036]
In the following description, it is assumed that the ultrasonic transducer 5 transmits an ultrasonic wave and the ultrasonic transducer 4 receives an ultrasonic wave. On the contrary, the ultrasonic transducer is used. 4 may transmit ultrasonic waves, and the ultrasonic transducer 5 may receive ultrasonic waves. Further, only one of the two ultrasonic transducers 4 and 5 may perform both ultrasonic transmission and reception.
[0037]
The amplifier 31 is connected to the ultrasonic transducer 4. The amplifier 31 amplifies the voltage of the output signal of the ultrasonic transducer 4 that has received the ultrasonic wave. An A / D converter 32 is connected to the amplifier 31. The A / D converter 32 digitizes the output signal of the amplifier 31 so that the arithmetic unit 33 can execute arithmetic processing.
[0038]
The computing unit 33 executes various types of computing processes as described below. The computing unit 33 includes a transfer function deriving unit 331 and an input waveform deriving unit 336. The transfer function deriving unit 331 includes a desired output waveform storage unit 332, an impulse signal storage unit 333, a conversion unit 334, and a ratio determining unit 335. The input waveform deriving unit 336 includes a converting unit 337 and a ratio determining unit 338.
[0039]
The desired output waveform storage unit 332 stores waveform information related to the output signal y (t) (t: time) having a desired waveform to be received by the ultrasonic transducer 4. In the present embodiment, the waveform information of the desired output signal y (t) stored in the desired output waveform storage unit 332 is a sine wave for one cycle whose frequency is within twice the center frequency of a transfer function described later. (Hereinafter referred to as “sin 1 wave”).
[0040]
The impulse signal storage unit 333 stores waveform information related to the impulse signal i (t) supplied to the ultrasonic transducer 5.
[0041]
The conversion unit 334 performs Fourier transform on the impulse signal i (t) stored in the impulse signal storage unit 333 and the impulse response signal o (t) output from the ultrasonic transducer 4. In the following description, the signals obtained by Fourier transforming the impulse signal i (t) and the impulse response signal o (t) by the conversion unit 334 are described as I (ω) and O (ω) (ω: angular frequency), respectively. .
[0042]
The ratio determining unit 335 calculates a ratio between the signals O (ω) and I (ω) obtained by the converting unit 334. This ratio is the transfer function H (ω) in the measurement system of the ultrasonic bone density measuring apparatus 1.
[0043]
On the other hand, the conversion unit 337 of the input waveform deriving unit 336 performs Fourier transform on the output signal y (t) having a sin 1 wave waveform stored in the desired output waveform storage unit 332. In the following description, a signal obtained by Fourier transforming the output signal y (t) by the conversion unit 337 is referred to as Y (ω).
[0044]
The ratio determining unit 338 calculates a ratio between the signal Y (ω) obtained by the converting unit 337 and the transfer function H (ω) obtained by the ratio determining unit 335. In the following description, this ratio is described as X (ω). X (ω) obtained by the ratio determining unit 338 is given to the converting unit 337 and subjected to inverse Fourier transform. Thereby, waveform information relating to the input signal x (t) to be input to the ultrasonic transducer 5 is obtained. The obtained waveform information is stored in the calculation unit 33.
[0045]
The input waveform generator 34 generates an input signal x (t) having the waveform derived by the input waveform deriving unit 336. The generated input signal x (t) is amplified in voltage by the amplifier 35 and supplied to the ultrasonic transducer 5.
[0046]
Next, an overall processing procedure of the ultrasonic bone density measuring apparatus 1 will be described with reference to FIGS. 5 and 6. FIG. 5 is a flowchart showing the entire processing procedure of the ultrasonic bone density measuring apparatus 1. FIG. 6 illustrates a series of processes for deriving an input signal x (t) to be input. Note that the shape of the signal depicted in FIG. 6 is merely shown as an example.
[0047]
First, in step S201, the reflecting surfaces 8a and 25a of the standoffs 8 and 25 are brought into contact with each other. When the standoffs 8 and 25 are in contact, ultrasonic jelly is applied to transmit ultrasonic waves. Thereafter, the process proceeds to step S203. In step S203, the impulse signal i (t) (see (a) in block B11 of FIG. 6) stored in the impulse signal storage unit 333 is applied to the ultrasonic transducer 5 and the ultrasonic transducer 5 is Excited. Next, the process proceeds to step S <b> 205, and the ultrasonic transducer 4 receives the ultrasonic wave generated by the ultrasonic transducer 5. Then, the ultrasonic transducer 4 outputs an impulse response signal o (t) (see (a) in the block B12 in FIG. 6). The impulse response signal o (t) is given to the calculator 33 through the amplifier 31 and the A / D converter 32. Thereafter, the process proceeds to step S207.
[0048]
In step S207, the conversion unit 334 performs Fourier transform on the impulse signal i (t). Thereby, the signal I (ω) (see (b) in the block B11 in FIG. 6) is obtained. Thereafter, the process proceeds to step S209. In step S209, the conversion unit 334 performs Fourier transform on the impulse response signal o (t) generated by the ultrasonic transducer 4 in step S205. Thereby, the signal O (ω) (see (b) in the block B12 in FIG. 6) is obtained. Thereafter, the process proceeds to step S211.
[0049]
In step S211, the frequency transfer function H (ω) of the measurement system of the ultrasonic bone density measuring apparatus 1 is calculated from the signal I (ω) and the signal O (ω) obtained in step S207 and step S209 by the ratio determining unit 335. Is calculated (see block B13 in FIG. 6). Here, the transfer function H (ω) can be calculated by the following equation. After the calculation, the process proceeds to step S213.
[0050]
[Expression 1]
Figure 0004312494
[0051]
In step S <b> 213, the input waveform deriving unit 336 calculates the input signal x (t) to be excited in the ultrasonic transducer 5. Specifically, first, the desired output signal y (t) (see (a) in the block B14 in FIG. 6) having a sin1 wave waveform stored in the desired output waveform storage unit 332 is converted by the conversion unit 337. Fourier transform is performed to obtain a signal Y (ω) (see (b) in block B14 in FIG. 6). Further, the ratio determining unit 338 determines a signal X (ω) (see B15 in FIG. 6) that is a ratio between the signal Y (ω) and the transfer function H (ω) derived in step S211. The signal X (ω) is expressed as the following equation.
[0052]
[Expression 2]
Figure 0004312494
[0053]
Then, the conversion unit 337 performs an inverse Fourier transform on the signal X (ω). Thereby, waveform information (see block B16 in FIG. 6) regarding the input signal x (t) to be input to the ultrasonic transducer 5 is obtained, and the process proceeds to step S215.
[0054]
In step S215, the reflecting surfaces 8a and 25a of the two standoffs 8 and 25 brought into contact in step S201 are separated again. A scissors as a subject is inserted between the separated standoffs 8 and 25, and the reflecting surfaces 8a and 25a of the standoffs 8 and 25 are brought into close contact with the scissors. When the standoffs 8 and 25 are in close contact with the heel, ultrasonic jelly is applied to transmit ultrasonic waves. Thereafter, in step S217, the input waveform generator 34 generates the input signal x (t) obtained in step S213. Further, the voltage of the input signal x (t) is amplified by the amplifier 35 and supplied to the ultrasonic transducer 5.
[0055]
Thereafter, the process proceeds to step S219, and the ultrasonic transducer 4 receives the ultrasonic wave that has passed through the eyelid. Thereafter, in step S221, the bone density of the heel is calculated based on the output signal output from the ultrasonic transducer 4. The calculated bone density of the heel is displayed on a display (not shown).
[0056]
FIG. 7 is a diagram depicting an example of each signal appearing at each processing stage in the ultrasonic bone density measuring apparatus 1. As shown in FIG. 7, the transfer function derived using the impulse signal (see (a) in block B21 in FIG. 7) and the impulse response signal (see (b) in block B21 in FIG. 7) is measured. Since the influence of the reflected wave generated in the bodies 2 and 3 and the standoffs 8 and 25 is incorporated, the input signal to be input obtained using this transfer function (see block B22 in FIG. 7) The output signal (see block B23 in FIG. 7) output from the ultrasonic transmitter / receiver 4 when the signal is supplied to the ultrasonic transmitter / receiver 5 is hardly affected by the reflected wave and is close to the sin1 wave. . In addition, the time width of the output signal from the ultrasonic transducer 4 is sufficiently narrow without using an expensive composite vibrator. Therefore, it is possible to obtain a highly accurate measurement result with excellent time resolution by an inexpensive apparatus configuration.
[0057]
In the present embodiment, the shape of the measurement bodies 2 and 3 is not affected by the reflected wave inside the measurement body regardless of the shape of the measurement body by the above method. It is possible to determine the most suitable shape for the measurement.
[0058]
In addition, as described above, the transfer function deriving unit 331 derives the transfer function by performing the Fourier transform and the ratio finding, so that the transfer function can be quickly obtained.
[0059]
Next, an ultrasonic bone density measuring apparatus according to a second embodiment of the present invention will be described with reference to the drawings. The ultrasonic bone density measuring apparatus of the present embodiment is different from the ultrasonic bone density measuring apparatus of the first embodiment in the transfer function derivation procedure. Hereinafter, the difference will be mainly described. In the following description, the same members as those in the first embodiment are denoted by the same reference numerals, and the description thereof is omitted.
[0060]
FIG. 8 is a block diagram showing a system configuration of the ultrasonic bone density measuring apparatus according to the present embodiment. The computing unit 36 includes a transfer function deriving unit 361 and an input waveform deriving unit 336. The transfer function deriving unit 361 includes a desired output waveform storage unit 332, an impulse signal storage unit 363, an impulse signal addition unit 364, a conversion unit 365, and a ratio determining unit 366.
[0061]
The impulse signal storage unit 363 transmits waveform information regarding the impulse signal i (t) supplied to the ultrasonic transducer 5 and ultrasonic transmission / reception when the impulse signal i (t) is supplied to the ultrasonic transducer 5. Waveform information relating to the impulse signal δ (t) to be added to the output signal o (t) from the waver 4 is stored.
[0062]
The impulse signal adding unit 364 adds the impulse signal δ (t) stored in the impulse signal storage unit 363 to the impulse response signal o (t) output from the ultrasonic transducer 4. The impulse response signal o (t) to which the impulse signal δ (t) is added by the impulse signal adding unit 364 will be referred to as a signal k (t) in the following description.
[0063]
The conversion unit 365 receives the impulse signal i (t) stored in the impulse signal storage unit 333 and the signal k (t) obtained by adding the impulse signal δ (t) in the impulse signal addition unit 364, respectively. Fourier transform. In the following description, signals obtained by Fourier transforming the impulse signal i (t) and the signal k (t) by the conversion unit 365 are described as I (ω) and K (ω), respectively.
[0064]
The ratio determining unit 366 calculates the ratio between the signal K (ω) and I (ω) obtained by the converting unit 365. This ratio is the transfer function G (ω) in the measurement system of the ultrasonic bone density measuring apparatus according to the present embodiment.
[0065]
The processing in the input waveform deriving unit 336 is the same except that the transfer function H (ω) is expressed as G (ω), and thus the description thereof is omitted.
[0066]
Next, the entire processing procedure of the ultrasonic bone density measuring apparatus according to the present embodiment will be described with reference to the flowchart of FIG. In FIG. 9, each processing procedure up to step S307 is the same as each processing procedure up to step S207 in the first embodiment, and thus description thereof is omitted.
[0067]
In step S309, the impulse signal adding unit 364 adds the impulse signal δ (t) stored in the impulse signal storage unit 363 to the impulse response signal o (t) generated by the ultrasonic transducer 4 in step S305. To do. At this time, the impulse signal δ (t) is superimposed on a flat portion of the impulse response signal o (t). Thereby, the signal k (t) is obtained. In step S309, the conversion unit 365 performs Fourier transform on the signal k (t). Thereby, the signal K (ω) is obtained. Thereafter, the process proceeds to step S311.
[0068]
In step S311, the frequency transfer function G (ω) of the measurement system of the ultrasonic bone density measuring apparatus is obtained from the signal I (ω) and the signal K (ω) obtained in steps S307 and S309 by the ratio determining unit 365. Calculated. Here, the transfer function G (ω) can be calculated by the following equation. After the calculation, the process proceeds to step S313.
[0069]
[Equation 3]
Figure 0004312494
[0070]
Each processing procedure after step S313 is the same as each processing procedure of steps S213 to S221 of the first embodiment.
[0071]
Here, the transfer function G (ω) derived by the ultrasonic bone density measuring apparatus according to the present embodiment will be described in comparison with the transfer function H (ω) in the first embodiment. FIGS. 10A and 10B are signal waveform diagrams of the impulse response signal o (t) and the transfer function H (ω), respectively, in the first embodiment. FIGS. 11A and 11B show the signal k (t) (the signal obtained by adding the impulse signal δ (t) to the impulse response signal o (t)) and the transfer function G (ω) in this embodiment, respectively. FIG.
[0072]
As shown in FIG. 10B, noise that is an unnecessary component is superimposed on the high frequency region and the low frequency region of the transfer function H (ω). When the input signal x (t) to be input is calculated based on the transfer function H (ω) on which noise is superimposed in this way, the high frequency of the transfer function H (ω) is added to the entire input signal x (t) obtained. The band component and the low frequency band component are superimposed. As a result, the output signal obtained when this input signal x (t) is applied to the ultrasonic transducer 5 is greatly deviated from a waveform such as a sin1 wave, and a highly accurate measurement result cannot be obtained.
[0073]
On the other hand, the transfer function G (ω) shown in FIG. 11B is obtained from the maximum value G1 (corresponding to the frequency fc) of the transfer function G (ω) by adding the impulse signal δ (t). Further, a waveform having no noise in the high frequency region and the low frequency region, in which the portion below the small value G2 is embedded (cut) and becomes a constant value, is obtained. Therefore, the high frequency region and low frequency region noise components of the transfer function G (ω) are not superimposed on the input signal z (t) calculated based on the transfer function G (ω). Therefore, the waveform shape of z (t) is stable, and transmission by an actual circuit is simplified.
[0074]
Note that the transfer function G (ω) in which a region of a certain value or less as shown in FIG. 11B is embedded can also be obtained by a complex operation after Fourier transform. However, this complex operation needs to be performed carefully, and requires a large amount of calculation due to the necessity for complex operation. In the present embodiment, a simple process of adding an impulse signal does not require a complicated algorithm, and the transfer function G (ω) as described above can be easily obtained.
[0075]
Further, according to the present embodiment, the same benefits as those of the first embodiment described above, that is, the time width of the output signal from the ultrasonic transducer 4 is sufficiently narrow with little influence of the reflected wave. Therefore, it is possible to obtain a highly accurate measurement result excellent in time resolution by an inexpensive apparatus configuration and to obtain a transfer function quickly.
[0076]
The preferred embodiments of the present invention have been described above. However, the present invention is not limited to the above-described embodiments, and various design changes can be made as long as they are described in the claims. For example, in the above-described embodiment, when the transfer function is derived, the reflection surfaces 8a and 25a of the standoffs 8 and 25 are brought into contact with each other, but the reflection surfaces 8a and 25a of the standoffs 8 and 25 are derived. The transfer function may be derived in a state where ultrasonic waves are transmitted through the subject without contact.
[0077]
In the above-described embodiment, two ultrasonic transducers are used, but one of them may be an ultrasonic transducer and the other may be an ultrasonic receiver.
[0078]
In the above-described embodiment, the waveform of the desired output signal y (t) is the sin 1 wave, but other waveforms may be used. For example, a rectangular wave may be used. In the above-described embodiment, an arithmetic unit that performs an operation is provided, but an arithmetic unit may not be provided, and the operation may be executed on a personal computer.
[0079]
In the above-described embodiment, the impulse signal is divided by the Fourier transform I (ω). However, if the impulse is sufficiently thin, it can be regarded as I (ω) = 1. In this case, H (ω) = O (ω).
[0080]
In the above-described embodiment, ultrasonic waves are transmitted and received using two ultrasonic transducers. However, ultrasonic waves may be transmitted and received using only one ultrasonic transducer. Specifically, the reflected wave of the ultrasonic wave transmitted by the ultrasonic transducer may be received by the same ultrasonic transducer to derive the transfer function.
[0081]
Further, the ultrasonic bone density measuring apparatus according to the second embodiment is not limited to a bone, and can be used as an ultrasonic measuring apparatus. In that case, the ultrasonic measurement device may not include a member such as a base or a measurement body.
[0082]
【The invention's effect】
As described above, according to the ultrasonic bone density measuring device of the present invention, it is possible to obtain a highly accurate measurement result by suppressing the extension of the time width of the transmitted ultrasonic wave and reducing the influence of the reflected wave. Is possible. In addition, according to the ultrasonic measurement apparatus of the present invention, it is possible to obtain a highly accurate measurement result that is not easily affected by noise.
[Brief description of the drawings]
FIG. 1 is a side view showing an overall configuration of an ultrasonic bone density measuring apparatus according to a first embodiment of the present invention.
FIG. 2 is a cross-sectional view taken along line II-II of the ultrasonic bone density measuring device shown in FIG.
3 is a cross-sectional view taken along line III-III of the ultrasonic bone density measuring apparatus shown in FIG.
4 is a block diagram showing a system configuration of the ultrasonic bone density measuring apparatus shown in FIG. 1. FIG.
FIG. 5 is a flowchart showing a processing procedure when bone density is measured by the ultrasonic bone density measuring apparatus shown in FIG. 1;
6 is a series of processing diagrams for deriving a waveform x (t) to be input to the ultrasonic bone density measuring apparatus shown in FIG. 1. FIG.
7 is a schematic view until an output waveform of the ultrasonic bone density measuring device shown in FIG. 1 is obtained. FIG.
FIG. 8 is a block diagram showing a system configuration of an ultrasonic bone density measuring apparatus according to a second embodiment of the present invention.
9 is a flowchart showing a processing procedure when bone density is measured by the ultrasonic bone density measuring apparatus shown in FIG.
10 is a signal waveform diagram in the case of deriving a transfer function without adding an impulse signal to the output signal of the ultrasonic bone density measuring device shown in FIG.
11 is a signal waveform diagram when a transfer function is derived by adding an impulse signal to the output signal of the ultrasonic bone density measuring apparatus shown in FIG. 8. FIG.
[Explanation of symbols]
1 Ultrasonic bone densitometer
2, 3 Measuring object
4, 5 Ultrasonic transducer
8, 25 Standoff
9 Diagnosis table (base)
33, 36 computing unit
331, 361 Transfer function deriving unit (transfer function deriving means)
332 desired output waveform storage unit
333, 363 Impulse signal storage unit
334, 365 conversion unit (first conversion means, second conversion means)
335, 366 Ratio finding unit (ratio finding means)
336 Input waveform deriving unit (input waveform deriving means)
364 Impulse signal adding unit (impulse signal adding means)

Claims (6)

被検体を載せるための基台と、
少なくともいずれか一方が他方の方向に移動可能に対向した、超音波透過性固体材料からなる一対の測定体と、
前記一対の測定体のいずれか一方に設けられた超音波送波器と、
前記一対の測定体のいずれか一方に設けられた超音波受波器と、
前記超音波送波器から前記超音波受波器に至る測定系の周波数伝達関数を導出するための伝達関数導出手段と、
出力信号の所望波形と、前記伝達関数導出手段によって導出された周波数伝達関数とに基づいて、前記超音波受波器が前記所望波形を有する出力信号を出力するときに前記超音波送波器に供給される入力信号の波形を導出するための入力波形導出手段と、
前記入力波形導出手段が導出した波形を有する入力信号を前記超音波送波器に供給するための入力波形供給手段とを備えており、
前記伝達関数導出手段が、
超音波を受波した前記超音波受波器からの出力信号にインパルス信号を加えた信号に基づいて、測定系の周波数伝達関数を導出することを特徴とする超音波骨密度測定装置。A base for placing the subject;
A pair of measuring bodies made of an ultrasonically transparent solid material, at least one of which is movably opposed to the other direction;
An ultrasonic transmitter provided in any one of the pair of measuring bodies;
An ultrasonic receiver provided on one of the pair of measuring bodies;
Transfer function deriving means for deriving a frequency transfer function of a measurement system from the ultrasonic transmitter to the ultrasonic receiver;
Based on the desired waveform of the output signal and the frequency transfer function derived by the transfer function deriving means, when the ultrasonic receiver outputs the output signal having the desired waveform, the ultrasonic transmitter Input waveform deriving means for deriving the waveform of the supplied input signal;
Input waveform supply means for supplying an input signal having a waveform derived by the input waveform deriving means to the ultrasonic wave transmitter ,
The transfer function deriving means,
An ultrasonic bone density measuring apparatus, wherein a frequency transfer function of a measurement system is derived based on a signal obtained by adding an impulse signal to an output signal from the ultrasonic receiver that has received an ultrasonic wave . 各測定体の形状が、前記超音波送波器から送波された超音波の超音波送波器から超音波受波器までの反射波の影響を前記伝達関数に含めて所望波形の出力信号を得るように、被検体の測定に適した形状に決定されていることを特徴とする請求項１に記載の超音波骨密度測定装置。  The shape of each measurement object is an output signal having a desired waveform including the influence of the reflected wave from the ultrasonic transmitter to the ultrasonic receiver of the ultrasonic wave transmitted from the ultrasonic transmitter in the transfer function. The ultrasonic bone density measuring device according to claim 1, wherein the ultrasonic bone density measuring device is determined to have a shape suitable for measurement of a subject. 超音波送波器と、
超音波受波器と、
前記超音波送波器から前記超音波受波器に至る測定系の周波数伝達関数を導出するための伝達関数導出手段と、
出力信号の所望波形と、前記伝達関数導出手段によって導出された周波数伝達関数とに基づいて、前記超音波受波器が前記所望波形を有する出力信号を出力するときに前記超音波送波器に供給される入力信号の波形を導出するための入力波形導出手段と、
前記入力波形導出手段が導出した波形を有する入力信号を前記超音波送波器に供給するための入力波形供給手段とを備えており、
前記伝達関数導出手段が、
超音波を受波した前記超音波受波器からの出力信号にインパルス信号を加えた信号に基づいて、測定系の周波数伝達関数を導出することを特徴とする超音波測定装置。An ultrasonic transmitter,
An ultrasonic receiver,
Transfer function deriving means for deriving a frequency transfer function of a measurement system from the ultrasonic transmitter to the ultrasonic receiver;
Based on the desired waveform of the output signal and the frequency transfer function derived by the transfer function deriving means, when the ultrasonic receiver outputs the output signal having the desired waveform, the ultrasonic transmitter Input waveform deriving means for deriving the waveform of the supplied input signal;
Input waveform supply means for supplying an input signal having a waveform derived by the input waveform deriving means to the ultrasonic wave transmitter,
The transfer function deriving means,
An ultrasonic measurement apparatus for deriving a frequency transfer function of a measurement system based on a signal obtained by adding an impulse signal to an output signal from the ultrasonic receiver that has received an ultrasonic wave. 前記伝達関数導出手段が、
前記超音波送波器への入力信号をフーリエ変換するための第１の変換手段と、
前記超音波受波器が超音波を受波したときの出力信号にインパルス信号を付加するためのインパルス信号付加手段と、
前記インパルス信号付加手段によってインパルス信号が付加された出力信号をフーリエ変換するための第２の変換手段と、
前記第１の変換手段によって得られた変換入力信号と、前記第２の変換手段によって得られた変換出力信号との比を求めるための求比手段とを備えていることを特徴とする請求項３に記載の超音波測定装置。The transfer function deriving means,
First conversion means for Fourier transforming an input signal to the ultrasonic transmitter;
An impulse signal adding means for adding an impulse signal to an output signal when the ultrasonic receiver receives an ultrasonic wave;
Second conversion means for Fourier transforming the output signal to which the impulse signal is added by the impulse signal addition means;
2. A ratio determining means for obtaining a ratio between a converted input signal obtained by the first converting means and a converted output signal obtained by the second converting means. 3. The ultrasonic measurement apparatus according to 3. 超音波送波器と、
超音波受波器と、
前記超音波送波器から被検査物を経て前記超音波受波器に至る測定系の周波数伝達関数を導出するための伝達関数導出手段と、
出力信号の所望波形と、前記伝達関数導出手段によって導出された周波数伝達関数とに基づいて、前記超音波受波器が前記所望波形を有する出力信号を出力するときに前記超音波送波器に供給される入力信号の波形を導出するための入力波形導出手段と、
前記入力波形導出手段が導出した波形を有する入力信号を前記超音波送波器に供給するための入力波形供給手段とを備えており、
前記伝達関数導出手段が、
超音波を受波した前記超音波受波器からの出力信号にインパルス信号を加えた信号に基づいて、測定系の周波数伝達関数を導出することを特徴とする超音波測定装置。An ultrasonic transmitter,
An ultrasonic receiver,
A transfer function deriving means for deriving a frequency transfer function of the measurement system from the ultrasonic wave transmitter to the ultrasonic wave receiver through the inspection object;
Based on the desired waveform of the output signal and the frequency transfer function derived by the transfer function deriving means, when the ultrasonic receiver outputs the output signal having the desired waveform, the ultrasonic transmitter Input waveform deriving means for deriving the waveform of the supplied input signal;
Input waveform supply means for supplying an input signal having a waveform derived by the input waveform deriving means to the ultrasonic wave transmitter,
The transfer function deriving means,
An ultrasonic measurement apparatus for deriving a frequency transfer function of a measurement system based on a signal obtained by adding an impulse signal to an output signal from the ultrasonic receiver that has received an ultrasonic wave. 超音波送波器から前記超音波受波器に至る測定系の周波数伝達関数を導出する伝達関数導出ステップと、A transfer function deriving step for deriving a frequency transfer function of the measurement system from the ultrasonic transmitter to the ultrasonic receiver;
超音波を受波した前記超音波受波器からの出力信号の所望波形と、前記伝達関数導出ステップで導出された周波数伝達関数とに基づいて、前記超音波受波器が前記所望波形を有する出力信号を出力するときに前記超音波送波器に供給される入力信号の波形を導出する入力波形導出ステップとを備えており、Based on the desired waveform of the output signal from the ultrasonic wave receiver that has received the ultrasonic wave and the frequency transfer function derived in the transfer function deriving step, the ultrasonic wave receiver has the desired waveform. An input waveform deriving step for deriving a waveform of an input signal supplied to the ultrasonic wave transmitter when outputting an output signal,
前記伝達関数導出ステップにおいて、In the transfer function derivation step,
出力信号にインパルス信号を加えた信号に基づいて、測定系の周波数伝達関数を導出することを特徴とする超音波測定方法。An ultrasonic measurement method, wherein a frequency transfer function of a measurement system is derived based on a signal obtained by adding an impulse signal to an output signal.
JP2003111647A 2003-04-16 2003-04-16 Ultrasonic bone density measuring device, ultrasonic measuring device, and ultrasonic measuring method Expired - Fee Related JP4312494B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2003111647A JP4312494B2 (en) 2003-04-16 2003-04-16 Ultrasonic bone density measuring device, ultrasonic measuring device, and ultrasonic measuring method
PCT/JP2004/005455 WO2004091409A1 (en) 2003-04-16 2004-04-16 Ultrasonic bone density measurement device, ultrasonic measurement device, ultrasonic bone density measurement method, and ultrasonic measurement method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2003111647A JP4312494B2 (en) 2003-04-16 2003-04-16 Ultrasonic bone density measuring device, ultrasonic measuring device, and ultrasonic measuring method

Publications (2)

Publication Number Publication Date
JP2004344181A JP2004344181A (en) 2004-12-09
JP4312494B2 true JP4312494B2 (en) 2009-08-12

Family

ID=33296004

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2003111647A Expired - Fee Related JP4312494B2 (en) 2003-04-16 2003-04-16 Ultrasonic bone density measuring device, ultrasonic measuring device, and ultrasonic measuring method

Country Status (2)

Country Link
JP (1) JP4312494B2 (en)
WO (1) WO2004091409A1 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102670257A (en) * 2012-05-16 2012-09-19 北京航空航天大学 Heel width measuring device and measuring method suitable for ultrasonic osteoporosis tester

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2575583B2 (en) * 1993-03-16 1997-01-29 アロカ株式会社 Estimation method of coefficient value related to ultrasonic wave propagation characteristics in living tissue
JPH08140971A (en) * 1994-11-17 1996-06-04 Hitachi Medical Corp Ultrasonic diagnostic apparatus
JPH0984788A (en) * 1995-09-26 1997-03-31 Aloka Co Ltd Tissue structure analyzing device
US5947902A (en) * 1996-09-27 1999-09-07 Metra Biosystems, Inc. Apparatus and method for calibration of an ultrasound transmission probe

Also Published As

Publication number Publication date
WO2004091409A1 (en) 2004-10-28
JP2004344181A (en) 2004-12-09

Similar Documents

Publication Publication Date Title
Zhang et al. Estimation of tissue’s elasticity with surface wave speed
US9345448B2 (en) System and method for non-invasive determination of tissue wall viscoelasticity using ultrasound vibrometry
Bernal et al. Material property estimation for tubes and arteries using ultrasound radiation force and analysis of propagating modes
EP3395256B1 (en) Method and apparatus for measuring viscoelastic parameter of viscoelastic medium
Lubinski et al. Lateral displacement estimation using tissue incompressibility
Zhang et al. Noninvasive method for estimation of complex elastic modulus of arterial vessels
Zheng et al. Detection of tissue harmonic motion induced by ultrasonic radiation force using pulse-echo ultrasound and Kalman filter
CA3017627C (en) Ultrasonic probe and ultrasonic detecting device provided with same
AU2016347143B2 (en) Elasticity detection method and device
EP0128635A2 (en) Measuring non-linear parameter of acoustic medium
WO1992007513A1 (en) Method and apparatus for breast imaging
Urban et al. Generalized response of a sphere embedded in a viscoelastic medium excited by an ultrasonic radiation force
Barannik et al. Doppler ultrasound detection of shear waves remotely induced in tissue phantoms and tissue in vitro
Pitre et al. Super-shear evanescent waves for non-contact elastography of soft tissues
Zhang et al. Noninvasive generation and measurement of propagating waves in arterial walls
Tabaru et al. Measurement of elastic properties of tissue by shear wave propagation generated by acoustic radiation force
US20190320997A1 (en) Methods for the use of inherent frequency shifting mechanisms for sensors response reading with continuous wave excitation
Wallace et al. Finite amplitude measurements of the nonlinear parameter B/A for liquid mixtures spanning a range relevant to tissue harmonic mode
Kujawska et al. Determination of nonlinear medium parameter B/A using model assisted variable-length measurement approach
JP4312494B2 (en) Ultrasonic bone density measuring device, ultrasonic measuring device, and ultrasonic measuring method
Baranger et al. The fundamental mechanisms of the Korotkoff sounds generation
US20070016038A1 (en) Ultrasonic method to determine bone parameters
Jang et al. Comparison of techniques for estimating shear-wave velocity in arterial wall using shear-wave elastography-FEM and phantom study
Guzina et al. Viscoelastic characterization of thin tissues using acoustic radiation force and model-based inversion
US20240090875A1 (en) System and Method for Non-Invasive Determination of Bladder Overactivity Using Ultrasound Vibrometry

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20060329

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090120

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090319

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20090512

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20090513

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120522

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

Ref document number: 4312494

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130522

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140522

Year of fee payment: 5

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees