JPH08280672A - 超音波骨評価装置 - Google Patents
超音波骨評価装置Info
- Publication number
- JPH08280672A JPH08280672A JP8563395A JP8563395A JPH08280672A JP H08280672 A JPH08280672 A JP H08280672A JP 8563395 A JP8563395 A JP 8563395A JP 8563395 A JP8563395 A JP 8563395A JP H08280672 A JPH08280672 A JP H08280672A
- Authority
- JP
- Japan
- Prior art keywords
- bone
- received signal
- spectrum
- peak
- evaluation value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
ができる超音波骨評価装置を提供する。 【構成】 第1フーリエ変換部20は、骨を透過した超
音波の受波信号をフーリエ変換することにより受波信号
スペクトルを求める。第2フーリエ変換部22は、受波
信号スペクトルに対して更にフーリエ変換を施す。受波
信号スペクトルに周期的なピークが存在する場合、第2
フーリエ変換部22によるフーリエ変換結果である二次
スペクトルに、その周期的なピークの繰り返し周波数に
対応するピークが現れる。そこで、評価値算出部24
は、二次スペクトルに現れるピークの位置に基づき骨の
評価値を算出する。この値が大きいことは、受波信号ス
ペクトルにおけるピーク繰り返しの周波数が大きいこ
と、すなわち受波信号スペクトルのピークの繰り返し周
期が小さいことを示し、これは骨梁間隔が粗であること
を示す。
Description
価・診断を行う超音波骨評価装置に関する。
粗鬆症・骨軟化症などの骨の疾患をもつ人が増加してお
り、その診断や予防の手法が要望されている。
来X線を利用した骨塩量測定装置があったが、近年これ
に加えて、生体に超音波を透過させて骨組織の評価・診
断を行う超音波骨評価装置が提案されている。このよう
な骨評価装置としては、例えば米国特許3,847,1
41号で提案されている装置がある。この装置は、人の
踵骨の評価を行う装置であり、踵に超音波を透過させて
踵骨を透過する超音波の音速を求め、この音速を踵骨の
評価値としている。
って求められる評価値は、骨中の音速や骨塩量など、い
わば巨視的な値であった。
側にある海綿骨とから成っており、このうち海綿骨は、
図3に示すように、比較的弾性波を吸収する硬い骨梁1
10と、骨梁110の間に存在し、軟らかくて弾性波の
伝達性が比較的よい骨髄120とから構成されている。
加齢などにより骨の量が減少してくると、骨梁110が
細くなったり消滅したりしてくるため、骨梁110同士
の間隔が大きくなると考えられるが、このような骨梁間
隔などの骨の微視的な構造は、上記骨中の音速や骨塩量
からは直接的に知ることができない。
な評価を行うために、骨の微視的構造を直接的に評価す
ることができる装置が要望されている。
ものであり、骨の微視的構造に関する評価値を求めるこ
とができる新たなタイプの超音波骨評価装置を提供する
ことを目的とする。
めに、本発明に係る超音波骨評価装置は、被検体に対し
て超音波を送波する送波手段と、前記被検体を透過した
超音波を受波する受波手段と、受波信号を周波数分析し
て受波信号スペクトルを求める周波数分析手段と、求め
られた受波信号スペクトルにおけるピークの周期性に基
づいて被検体内の骨に関する評価値を演算する評価値演
算手段とを有することを特徴とする。
前記周波数分析手段が前記受波信号をフーリエ変換する
ことにより前記受波信号スペクトルを求め、前記評価値
演算手段が前記受波信号スペクトルに更にフーリエ変換
を施し、この結果得られる二次スペクトルに現れるピー
クの位置に基づいて前記評価値を求めることを特徴とす
る。
超音波を、図3に示すような海綿骨骨を含む被検体に照
射し、透過超音波を受波した場合、得られる受波信号の
スペクトルは図4に示すようになる。図4から分かるよ
うに、受波信号のスペクトルには図中Rに示す周期的な
ピークが観察される。
の微視的構造を反映したものと考えられる。すなわち、
海綿骨は、超音波に対して一種のフィルタとして作用
し、このフィルタとしての効果は、海綿骨の微視的な構
造、特に骨梁同士の間隔を反映する。骨梁と骨髄の超音
波伝搬特性からいって、骨梁間隔に合った波長の超音波
ほど透過しやすいので、受波信号のスペクトルにおい
て、海綿骨を透過しやすい波長に対応する周波数成分の
強度が高くなるからである。従って、海綿骨内の骨梁の
間隔が狭い場合は、比較的波長の短い、すなわち比較的
高周波の成分がよく透過されるため、受波信号のスペク
トルには高い周波数を基本周波数とする高調波成分が強
く現れ、これらが周期的なピークを形成する。このた
め、骨梁が密な場合は、受波信号のスペクトルのピーク
の周期が大きくなる。これに対し、海綿骨の骨梁間隔が
大きくなると、よく透過される超音波の波長が大きく
(すなわち、周波数が小さく)なるため、受波信号スペ
クトルに現れるピークの周期が小さくなる。
ルの周期性と骨の微視的構造との関係に注目し、受波信
号スペクトルにおけるピークの繰り返しの周期性を求
め、これに基づいて骨の微視的構造に関する評価値を算
出する。
段から放射された超音波は、骨を透過して受波手段によ
って検出される。その際の受波信号が、周波数分析手段
によって分析され、受波信号のスペクトルが求められ
る。評価値算出手段は、その受波信号スペクトルのピー
クの周期性を検出し、この周期性に基づいて骨の評価値
を算出する。
信号をフーリエ変換することにより得られた受波信号ス
ペクトルに対し更にフーリエ変換を施し、二次スペクト
ルを求める。受波信号スペクトルに周期性が存在する
と、その周期性に対応したピークが二次スペクトルに現
れるので、この二次スペクトルにおいてピークが現れる
位置に基づいて骨の評価値を定める。
施例を図面に基づいて説明する。
概略構成を示すブロック図である。図において、互いに
対向する一対の振動子10a,10bは、図示しない走
査ユニットに固定されている。被検体を検査する際に
は、振動子10a,10bの間に被検体(例えば、踵)
100を配置して、超音波の送受波を行う。このとき、
図示しない走査制御部を通じて走査ユニット全体を駆動
することにより、一対の振動子10a,10bを、被検
体100に対して順次走査することができる。
ユニット12a,12bに接続されている。送受信ユニ
ット12a,12bは、それぞれ振動子10a,10b
に対して送信パルスを供給すると共に、振動子10a,
10bからの受波信号を取り込む。本実施例では、振動
子10a,10bはどちらも送受両用であり、適宜使い
分けることができる。なお、本実施例においては、送波
側の振動子からは、多くの周波数成分を含んだインパル
ス状の超音波パルスが送波される。
御部14からの制御に従い、送信パルスの発生や受波信
号の取り込みを行う。また、送受信ユニット12a,1
2bは、A/D変換器16に接続されており、送受信ユ
ニット12a又は12bによって振動子10a又は10
bから取り込まれた受波信号は、A/D変換器16でデ
ジタル化された後、計測制御部14を介して演算処理部
18に入力される。
であり、図に示すように、演算処理部18は、第1フー
リエ変換部20、第2フーリエ変換部22及び評価値算
出部24を有している。
た受波信号は、まず第1フーリエ変換部20に入力され
る。第1フーリエ変換部20は、受波信号をフーリエ変
換することにより、受波信号のスペクトルを求める。求
められた受波信号スペクトルのデータは、第2フーリエ
変換部22に入力される。第2フーリエ変換部22は、
受波信号スペクトルのデータに対して更にフーリエ変換
を施す。すなわち、第2フーリエ変換部22は、受波信
号スペクトルを周波数の関数と見て、その受波信号のス
ペクトル自体を更に周波数分析する。受波信号スペクト
ルに周期的なピークが存在する場合は、この受波信号ス
ペクトル自体の周波数分析により、その周期性が検出さ
れる。
を有する受波信号スペクトルをフーリエ変換したときに
得られるスペクトルの例である。図5において、縦軸は
強度であり、横軸は受波信号スペクトルに現れる繰り返
し成分の周波数である。以下、この受波信号スペクトル
のフーリエ変換結果のことを二次スペクトルと呼ぶ。
トルには特徴的なピークAが現れる。このピークAは、
受波信号スペクトルにおけるピークの繰り返しによるも
のであると考えられる。従って、二次スペクトルのピー
クから受波信号スペクトルの周期性を検出することがで
きる。
は、第2フーリエ変換部22で得られた二次スペクトル
のデータを受け取り、この二次スペクトルにおけるピー
クの位置Tを検出し、これを骨に関する評価値として出
力する。この評価値が大きいということは、受波信号ス
ペクトルにおけるピーク繰り返しの周波数が大きいこ
と、すなわち受波信号スペクトルのピークの繰り返し周
期が小さいことを示し、これは骨梁間隔が粗であること
を示す。逆に、この評価値が小さいことは、骨梁間隔が
密であることを示す。
スペクトルを更に周波数分析することにより、受波信号
スペクトルに含まれる繰り返し成分の周期性(すなわ
ち、ピークの間隔)に関する評価値を得ることができ
る。受波信号スペクトルのピーク間隔は骨梁間隔を反映
しているので、その評価値は、骨の微視的構造における
一つの特徴量である骨梁間隔を評価するための指標値と
なる。
ピークの繰り返しに注目したものであったが、このピー
クの繰り返しは、1MHzを超えるかなり高い周波数帯
域にならないと現れない。ところが、骨を含む生体自体
は高い周波数成分に対する透過性が小さいので、この繰
り返し部分のレベルは、透過性の大きい比較的周波数の
低い成分の受波信号レベル(図4のスペクトルのメイン
ローブM)に比べてかなり低くなってしまう。このた
め、ピークの繰り返し部分は、回路自体のノイズの影響
を受けやすい。そこで、本実施例では、被検体に対して
何回も超音波パルスを透過させて受波信号データを採取
し、それらを加算処理(あるいは加算平均処理)するこ
とにより、ランダムに発生するノイズの影響を低減す
る。
返しは、メインローブより遥かに高い周波数帯に現れる
ので、ピークの周期性の検出にはメインローブ付近のデ
ータは基本的に必要ない。そこで、ノイズの影響を低減
する方法としては、メインローブ付近の受波信号レベル
が飽和してもよいものとして超音波パルスの送信パワー
を大きくし、ピークの繰り返しが現れる高周波数の領域
の受波レベルを大きくする方法もある。
号スペクトルのピークの周期性を検出するために、受波
信号スペクトルをフーリエ変換して二次スペクトルを求
めたが、これに限らず、受波信号スペクトル自体からピ
ークを検出してピーク同士の間隔を求め、そのピーク間
隔を評価値として採用することもできる。
受波信号スペクトルに現れるピークの繰り返しの周期性
を検出することにより、骨の微視的構造を示す特徴量の
一つである骨梁間隔に関する評価値を得ることができ
る。
示すブロック図である。
の機能ブロック図である。
号のスペクトルを示す図である。
して得られる二次スペクトルを示す図である。
ット、14 計測制御部、16 A/D変換器、18
演算処理部、20 第1フーリエ変換部、22第2フー
リエ変換部、24 評価値算出部。
Claims (2)
- 【請求項1】 被検体に対して超音波を送波する送波手
段と、 前記被検体を透過した超音波を受波する受波手段と、 受波信号を周波数分析して受波信号スペクトルを求める
周波数分析手段と、 求められた受波信号スペクトルにおけるピークの周期性
に基づいて被検体内の骨に関する評価値を演算する評価
値演算手段と、 を有することを特徴とする超音波骨評価装置。 - 【請求項2】 請求項1に記載の超音波骨評価装置にお
いて、 前記周波数分析手段は、前記受波信号をフーリエ変換す
ることにより前記受波信号スペクトルを求め、 前記評価値演算手段は、前記受波信号スペクトルに更に
フーリエ変換を施し、この結果得られる二次スペクトル
に現れるピークの位置に基づいて前記評価値を求めるこ
とを特徴とする超音波骨評価装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP08563395A JP3472376B2 (ja) | 1995-04-11 | 1995-04-11 | 超音波骨評価装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP08563395A JP3472376B2 (ja) | 1995-04-11 | 1995-04-11 | 超音波骨評価装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08280672A true JPH08280672A (ja) | 1996-10-29 |
| JP3472376B2 JP3472376B2 (ja) | 2003-12-02 |
Family
ID=13864245
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP08563395A Expired - Lifetime JP3472376B2 (ja) | 1995-04-11 | 1995-04-11 | 超音波骨評価装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3472376B2 (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2019505285A (ja) * | 2016-01-15 | 2019-02-28 | バタフライ ネットワーク,インコーポレイテッド | 超音波信号処理回路並びに関連する装置及び方法 |
| JP2020508461A (ja) * | 2017-02-23 | 2020-03-19 | サザン リサーチ インスティチュート | スペクトル領域及び時間領域を使用する超音波検査システム |
-
1995
- 1995-04-11 JP JP08563395A patent/JP3472376B2/ja not_active Expired - Lifetime
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2019505285A (ja) * | 2016-01-15 | 2019-02-28 | バタフライ ネットワーク,インコーポレイテッド | 超音波信号処理回路並びに関連する装置及び方法 |
| JP2020508461A (ja) * | 2017-02-23 | 2020-03-19 | サザン リサーチ インスティチュート | スペクトル領域及び時間領域を使用する超音波検査システム |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3472376B2 (ja) | 2003-12-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4338026B2 (ja) | 骨の非侵襲的検査のための方法および装置 | |
| Droin et al. | Velocity dispersion of acoustic waves in cancellous bone | |
| JP3432204B2 (ja) | 超音波診断装置 | |
| US7785259B2 (en) | Detection of motion in vibro-acoustography | |
| AU689570B2 (en) | Method and device for estimating and characterising bone properties | |
| JP2005510283A6 (ja) | 骨の非侵襲的検査のための方法および装置 | |
| JP2883290B2 (ja) | 超音波骨評価装置 | |
| JPH06504453A (ja) | 乳房組織の画像を生成するための方法および装置 | |
| WO2008036513A2 (en) | Ultrasonic method and apparatus for assessment of bone | |
| Gerig et al. | Statistics of ultrasonic scatterer size estimation with a reference phantom | |
| CN108024795A (zh) | 频谱解析装置、频谱解析方法以及超声波摄像装置 | |
| JPH0713631B2 (ja) | 媒体の超音波エコーグラフィック検査方法及び装置 | |
| JP2596701B2 (ja) | 骨評価装置 | |
| US20110046517A1 (en) | Method and system for monitoring skeletal defects | |
| KR20180037351A (ko) | 초음파 비선형 변수를 이용한 골구조 예측 장치 및 방법 | |
| JP3472376B2 (ja) | 超音波骨評価装置 | |
| KR101840349B1 (ko) | 초음파 합주파수 성분을 이용한 골밀도 예측 장치 및 방법 | |
| KR102364309B1 (ko) | 초음파 감쇠계수 및 후방산란계수를 이용한 골밀도 및 골구조 예측방법 | |
| WO2021010836A1 (en) | Method and system for using wave analysis for speed of sound measurement | |
| US20060253025A1 (en) | Ultrasonic Bone Assessment Apparatus and Method | |
| CN113109825A (zh) | 基于Radon变换的长骨相控超声信号定征与骨质评价系统 | |
| WO2001077707A1 (en) | Ultrasound contrast imaging with multiple-pulse excitation waveforms | |
| KR102303922B1 (ko) | 초음파 감쇠계수 및 위상속도를 이용한 골밀도 및 골구조 예측방법 | |
| JP3317988B2 (ja) | 超音波による骨診断装置 | |
| JPH0348789A (ja) | Cwドプラ装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080912 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100912 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100912 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120912 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130912 Year of fee payment: 10 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| EXPY | Cancellation because of completion of term |