JPS5919852A - 超音波生体組織特性測定方式 - Google Patents
超音波生体組織特性測定方式Info
- Publication number
- JPS5919852A JPS5919852A JP57129902A JP12990282A JPS5919852A JP S5919852 A JPS5919852 A JP S5919852A JP 57129902 A JP57129902 A JP 57129902A JP 12990282 A JP12990282 A JP 12990282A JP S5919852 A JPS5919852 A JP S5919852A
- Authority
- JP
- Japan
- Prior art keywords
- frequency
- tissue
- value
- biological tissue
- spectrum
- 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.)
- Pending
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0858—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/11—Analysing solids by measuring attenuation of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/34—Generating the ultrasonic, sonic or infrasonic waves, e.g. electronic circuits specially adapted therefor
- G01N29/341—Generating the ultrasonic, sonic or infrasonic waves, e.g. electronic circuits specially adapted therefor with time characteristics
- G01N29/343—Generating the ultrasonic, sonic or infrasonic waves, e.g. electronic circuits specially adapted therefor with time characteristics pulse waves, e.g. particular sequence of pulses, bursts
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/34—Generating the ultrasonic, sonic or infrasonic waves, e.g. electronic circuits specially adapted therefor
- G01N29/348—Generating the ultrasonic, sonic or infrasonic waves, e.g. electronic circuits specially adapted therefor with frequency characteristics, e.g. single frequency signals, chirp signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/895—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques characterised by the transmitted frequency spectrum
- G01S15/8954—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques characterised by the transmitted frequency spectrum using a broad-band spectrum
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/0289—Internal structure, e.g. defects, grain size, texture
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Radar, Positioning & Navigation (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Immunology (AREA)
- Chemical & Material Sciences (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
（発明の技術分野〕
本発明は超音波の反射波を用いて生体組織の如く減衰定
数が周波数に比例（比例定数β）シ、且つ反射係数が各
種生体組織特性値により定まる周波数の関数である媒体
において、β等の各種特性値を測定する方式にかかわり
、特にパワー・スペクトル測定値から理論式・実験式に
回帰することによって、各種特性値を求める方式に関す
る。
数が周波数に比例（比例定数β）シ、且つ反射係数が各
種生体組織特性値により定まる周波数の関数である媒体
において、β等の各種特性値を測定する方式にかかわり
、特にパワー・スペクトル測定値から理論式・実験式に
回帰することによって、各種特性値を求める方式に関す
る。
従来、ｆｌ音波の透過における減衰定数が周波数ｆに比
例し、その比例定数βは組織特性を示す一つの量である
こと、又反射係数は周波数ｆのｎ乗に比例し、この墓指
数ｎは組織特性を示すもう一つの量であることが、実験
的に知られている。
例し、その比例定数βは組織特性を示す一つの量である
こと、又反射係数は周波数ｆのｎ乗に比例し、この墓指
数ｎは組織特性を示すもう一つの量であることが、実験
的に知られている。
又９本発明者の上田等は反射係数が。
６．己。−ｃｄ１′ノ２
の形であることを理論的に示している。ここに。
ｂ、ｄは組織特性と関連を持つ量である。
このように反射係数がｆの関数である場合について、パ
ワー・スペクトルの形状から生体組織特性値を求める方
式は一般的なものはなく９本発明者の三輪等によるエネ
ルギー・レシオ法により。
ワー・スペクトルの形状から生体組織特性値を求める方
式は一般的なものはなく９本発明者の三輪等によるエネ
ルギー・レシオ法により。
該関数形がｒのｎ乗である場合についてｎ、βを求める
方式が特許出願されている。
方式が特許出願されている。
この方法は有効であるが、少なくとも３つの狭帯域のエ
ネルギー成分に着目するため、所謂スカロソピングとし
て知られているスペクトルの局所的凸凹による誤差を招
き易い。このため、有効帯域内で多数の３周波数の組に
ついて計算し、得られたｎ、βの集合について統計処理
を施す必要があり、計算量が多い欠点があった。
ネルギー成分に着目するため、所謂スカロソピングとし
て知られているスペクトルの局所的凸凹による誤差を招
き易い。このため、有効帯域内で多数の３周波数の組に
ついて計算し、得られたｎ、βの集合について統計処理
を施す必要があり、計算量が多い欠点があった。
本発明の目的は、超音波反射波の受信信号から生体の伝
達関数を求め、生体組織特性値を求めるに当って。
達関数を求め、生体組織特性値を求めるに当って。
■組織領域端界の不連続変化を避けたスペクトルの形状
を求める方法と。
を求める方法と。
■その形状の理論・実験から推定される関数形が指数関
数項と非指数関数項との積で表わされるのに対し、実測
値からの回帰が容易な方式を与え。
数項と非指数関数項との積で表わされるのに対し、実測
値からの回帰が容易な方式を与え。
■その回帰によって該関数形に含まれるパラメータを決
定し、更に生体組織特性値を求める方式を提供すること
にある。
定し、更に生体組織特性値を求める方式を提供すること
にある。
本発明は、生体組織伝達関数の周波数応答スペクトルを
正規化して、その絶対値には関係なく。
正規化して、その絶対値には関係なく。
その形状のみに着目するとともに、理論または実験より
推定される関数形が指数関数項と非指数関数項との積で
与えられる場合、非指数関数項で除して対数をとること
によって得られる測定値から回帰することにより、指数
関数項の指数部に含まれる生体組織特性値に関連するパ
ラメータを求めるようにしたものである。
推定される関数形が指数関数項と非指数関数項との積で
与えられる場合、非指数関数項で除して対数をとること
によって得られる測定値から回帰することにより、指数
関数項の指数部に含まれる生体組織特性値に関連するパ
ラメータを求めるようにしたものである。
生体組織に第１図に示す如く生体表面から生体深部へ（
方向Ｚに）トランスデユーサ−〇によって超音波パルス
（中心周波数ｆｏ、帯域２Ω）を送信し、そのパルス波
形が生体内容部位に音速Ｃの速さで順次到達し、夫々か
ら逐次反射波が音速Ｃで逆方向に伝播し、再びトランス
デユーサ−■によって受（言されるものとする。
方向Ｚに）トランスデユーサ−〇によって超音波パルス
（中心周波数ｆｏ、帯域２Ω）を送信し、そのパルス波
形が生体内容部位に音速Ｃの速さで順次到達し、夫々か
ら逐次反射波が音速Ｃで逆方向に伝播し、再びトランス
デユーサ−■によって受（言されるものとする。
今、深さ２における組織の特性を測定するものとす−る
。生体は２の深さまでにｉ種の組織領域から構成されて
いるとする。又、音速Ｃはほぼ一定と仮定する。
。生体は２の深さまでにｉ種の組織領域から構成されて
いるとする。又、音速Ｃはほぼ一定と仮定する。
体表■、即ち２＝０から送信された音圧パルスは、各領
域ｉの中では周波数ｆに比例した減衰定数αｉの減衰を
受ける。この比例定数をａｉとすると。
域ｉの中では周波数ｆに比例した減衰定数αｉの減衰を
受ける。この比例定数をａｉとすると。
αｉ＝ａ＋＋βｒ＊ｆ （ａｉは定数）となり、ａ
ｉは組織ｉに対応した特性を示すパラメータで減衰傾斜
と略称される。
ｉは組織ｉに対応した特性を示すパラメータで減衰傾斜
と略称される。
又、各領域の中では周波数ｆの関数であるパワー反射係
数ｒ（ｆ）を示す。
数ｒ（ｆ）を示す。
又、各領域相互の境界では大きく音響インピーダンスが
変化し、且つその表面が略滑らかである場合が多く、こ
の境界１例えば領域ｉ−１から領域ｉに通過するときは
１段階的な若干の透過１員失を生じる。この透過率をτ
ｊとする。
変化し、且つその表面が略滑らかである場合が多く、こ
の境界１例えば領域ｉ−１から領域ｉに通過するときは
１段階的な若干の透過１員失を生じる。この透過率をτ
ｊとする。
深さ２からの反射波が逆方向に領域ｉからｉ−１に通過
するときも同様で、この透過率をτ′ｉとする。τｉと
τ′ｉは周波数特性を持たないと考えてよい。
するときも同様で、この透過率をτ′ｉとする。τｉと
τ′ｉは周波数特性を持たないと考えてよい。
さて、トランスデユーサ■は送信後、順次深い所からの
反射波を連続的に受信するが１その中。
反射波を連続的に受信するが１その中。
深さ２に対応する反射波が送信後ｔ ＝ ２ ｚ ／
Ｃの時刻に受信されるので、その時点付近のある微小時
間幅τ内の反射波形を分析することで、深さ２における
組織特性を求めることができる。
Ｃの時刻に受信されるので、その時点付近のある微小時
間幅τ内の反射波形を分析することで、深さ２における
組織特性を求めることができる。
深さ２からの反射波受信信号のパワー・スペクトルＥ
ｒ ｆｆ）を求めることは、よく知られたＦＦＴ等の周
波数分析装置で容易に実現される。
ｒ ｆｆ）を求めることは、よく知られたＦＦＴ等の周
波数分析装置で容易に実現される。
このＥ ｒ （ｆ）はトランスデユーサの周波数特性。
深さ２におりるビームの収束の周波数特性、その他計測
系等に起因する伝達関数と、生体中の透過・反射等に起
因する生体組織の伝達関数との積の自乗で与えられる。
系等に起因する伝達関数と、生体中の透過・反射等に起
因する生体組織の伝達関数との積の自乗で与えられる。
今、水中などで深さ２の位置に標準反射体をおいて、そ
の受信波のパワー・スペクトルＥ ｏ ［ｆｌを求める
と、Ｅｏ（ｆ）は計測系の伝達関数の自乗を示すと考え
られる。ここで、Ｅｒ（ｆ）をＥ ｏ ｆｆｌで除して
正規化すると、生体組織固有の伝達関数の自乗が得られ
る。この実測された生体組織伝達関数のパワー・スペク
トルをＲ（ｆｌで示す。
の受信波のパワー・スペクトルＥ ｏ ［ｆｌを求める
と、Ｅｏ（ｆ）は計測系の伝達関数の自乗を示すと考え
られる。ここで、Ｅｒ（ｆ）をＥ ｏ ｆｆｌで除して
正規化すると、生体組織固有の伝達関数の自乗が得られ
る。この実測された生体組織伝達関数のパワー・スペク
トルをＲ（ｆｌで示す。
一方において、Ｒ（ｆｌは上述の説明から１次式で表わ
されるはずである。
されるはずである。
ｋ：ｆに依存しない定数
１１：領域ｉ内の通過長
以下、Ｒ（ｆｌと（１）式の行とを比較することによっ
て９組織特性を求める方式を詳細に説明する。
て９組織特性を求める方式を詳細に説明する。
一つのトランスデユーサではその有効帯域幅内のｆにつ
いてのみしかＲ（［１を求めることができないので、十
分広いｆの範囲でＲ（ｆ）を求めるためには複数のトラ
ンスデユーサで夫々の帯域幅について夫々のＲ（ｆ）を
求めてつなげればよい。
いてのみしかＲ（［１を求めることができないので、十
分広いｆの範囲でＲ（ｆ）を求めるためには複数のトラ
ンスデユーサで夫々の帯域幅について夫々のＲ（ｆ）を
求めてつなげればよい。
又、実際の測定に当っては、近接する反射体からの反射
波の干渉、その他によってスペクトルの形状は局所的な
凸凹を示し、所謂スペクトル・スカロソピングの現象を
呈し、゛著しくスペクトル形状の測定に誤差をもたらす
。この防止のためにはその測定点の前後・左右・上下環
を測定したり。
波の干渉、その他によってスペクトルの形状は局所的な
凸凹を示し、所謂スペクトル・スカロソピングの現象を
呈し、゛著しくスペクトル形状の測定に誤差をもたらす
。この防止のためにはその測定点の前後・左右・上下環
を測定したり。
複数回測定したりして、それらを時間的・空間的に平均
する等の統計的処理が必要になる。
する等の統計的処理が必要になる。
さて、（１）式においである特定周波数ｆｏ（例えば（
１）式の最大値を与える周波数ｆｍ等）におけるＲ（ｆ
ｏ）で（１）式を規格化すると。
１）式の最大値を与える周波数ｆｍ等）におけるＲ（ｆ
ｏ）で（１）式を規格化すると。
となり、に、τｉ、τ′ｉの因子を除去することができ
る。未知数であるτｉ、τ′ｉが消去されることの効果
は大きい。
る。未知数であるτｉ、τ′ｉが消去されることの効果
は大きい。
実測値でもｆｏ及びｆｏにおけるＲ（ｆｏ）を求めるこ
とができるので、同様にＲ（ｆｏ）で規格化するを求め
ることができる。即ち、τｉ、τ′ｉの影響を受けるこ
となく、Ｐ（ｆ）とＱ Ｉｆ）が等しいとしてＱ （ｆ
）に含まれるパラメータを回帰により決定することがで
きる。
とができるので、同様にＲ（ｆｏ）で規格化するを求め
ることができる。即ち、τｉ、τ′ｉの影響を受けるこ
となく、Ｐ（ｆ）とＱ Ｉｆ）が等しいとしてＱ （ｆ
）に含まれるパラメータを回帰により決定することがで
きる。
次にｒ（［１の関数形としてはその近似度により種々の
ものがあり得るが１次の実験式 ％式％ （４） （ａ：定数、ｎ：組織により定まる定数）及び、上田等
による理論式 ｂニドランスデューサ寸法１曲率に関係するｆの関数（
既知） ｂ＝組織の平均的微細構造に依存する定数σに：組織の
値送受方向における空間的自己相関距離 Ｃ：音速 の場合について説明する。
ものがあり得るが１次の実験式 ％式％ （４） （ａ：定数、ｎ：組織により定まる定数）及び、上田等
による理論式 ｂニドランスデューサ寸法１曲率に関係するｆの関数（
既知） ｂ＝組織の平均的微細構造に依存する定数σに：組織の
値送受方向における空間的自己相関距離 Ｃ：音速 の場合について説明する。
（４）式、（５）式は一見著しく異なった形をしている
が、実用されるある周波数範囲では両者はほぼ一致する
。
が、実用されるある周波数範囲では両者はほぼ一致する
。
（４）式、又は（５）式を（２）式に代入すると、
Ｑ４．Ｑ５が得られる。
Ｑ４．Ｑ５が得られる。
Ｑ ４ ＝ Ａ４・ｆηｅｘｐ［（−４，ｇρｄＪｆ］
（６１Ｑ５−）＼、、ｂ’・ｆ’：ｅｘｐ［ニ
ー４（２ゴｐｄｔ）ｆ−Ｃし乙２）ぜ４シ＝ｆｒ”ノ
（７）（６）式、（７）式から判るように、Ｑは指数関
数の項と非指数関数の項との積で与えられることが特徴
的である。両式を夫々非指数関数項で除し、対数をとる
と。
（６１Ｑ５−）＼、、ｂ’・ｆ’：ｅｘｐ［ニ
ー４（２ゴｐｄｔ）ｆ−Ｃし乙２）ぜ４シ＝ｆｒ”ノ
（７）（６）式、（７）式から判るように、Ｑは指数関
数の項と非指数関数の項との積で与えられることが特徴
的である。両式を夫々非指数関数項で除し、対数をとる
と。
Ａ Ｑ４／ Ａ４− ’７？−Ｊｎｆ ＝ −４（Σβ
ｉ Ｊｉ ） ・ｆ （８１Ｊ！５ＱｒＪ／Ａｒ
、−、ｌｎｌ；ｒ”＝ −４（Ｚβ；ｌｂ） ・ ｆ
−（？）白９）ここでＱ４．Ｑ５の代りに実測値Ｐを
代入した（８）式（９）式の左辺を計算し、（８）式又
は（９）式が最も精度良く成立するように回帰計算する
ことにより、パラメータｎ、’Ｅｏ；、ｉｉ、 、ｙ
、等を求めることができる。
ｉ Ｊｉ ） ・ｆ （８１Ｊ！５ＱｒＪ／Ａｒ
、−、ｌｎｌ；ｒ”＝ −４（Ｚβ；ｌｂ） ・ ｆ
−（？）白９）ここでＱ４．Ｑ５の代りに実測値Ｐを
代入した（８）式（９）式の左辺を計算し、（８）式又
は（９）式が最も精度良く成立するように回帰計算する
ことにより、パラメータｎ、’Ｅｏ；、ｉｉ、 、ｙ
、等を求めることができる。
（８）式は左辺に未知数ｎがあるが、ｎを種々に変えて
右辺が一次式となるｎを決定すればよい。
右辺が一次式となるｎを決定すればよい。
夫々の深さＺにおけるΣβ７７ｔを求めることによって
、差分によって各ｌにおけるβｉを求めることができる
。βが深さ２の連続関数の時は、この線積分は Ｉβｄ
ｚとなるので、微分によりβｆＺｌを求めたり、Ｘ線Ｃ
Ｔのアルゴリスムで各方向での線積分から計算すること
ができる。
、差分によって各ｌにおけるβｉを求めることができる
。βが深さ２の連続関数の時は、この線積分は Ｉβｄ
ｚとなるので、微分によりβｆＺｌを求めたり、Ｘ線Ｃ
Ｔのアルゴリスムで各方向での線積分から計算すること
ができる。
このようにして、ｎ、β、ｌを一つの測定方向走査線上
で求めることができる。この走査線をある面内で走査す
ることにより、その面上のｎ、βｔの分布像をもとめる
こともできる。
で求めることができる。この走査線をある面内で走査す
ることにより、その面上のｎ、βｔの分布像をもとめる
こともできる。
以上原理を述べたが、実際に装置化するに当っては、あ
る深さ２でのサンプル点からのパラメータを抽出する技
術はドツプラー測定等で広く行われており、その抽出波
形をフーリエ分析し、パワー・スペクトルを求める装置
は、ＤＦＴ （デジタル・フーリエ変換装置）等で広く
知られているものが容易に利用できる。その出力は計算
機に容易に入力され、それ以降のすべての計算はプログ
ラムによって遂行できる。以下システムの概略につき、
第２図で説明する。
る深さ２でのサンプル点からのパラメータを抽出する技
術はドツプラー測定等で広く行われており、その抽出波
形をフーリエ分析し、パワー・スペクトルを求める装置
は、ＤＦＴ （デジタル・フーリエ変換装置）等で広く
知られているものが容易に利用できる。その出力は計算
機に容易に入力され、それ以降のすべての計算はプログ
ラムによって遂行できる。以下システムの概略につき、
第２図で説明する。
タイミング制御部１からは（イ）に示す超音波送信同期
信号Ｐｄが送信回路２に送られ、超音波送信に必要なだ
けの電力を持ったパルスで超音波振動子３が駆動され、
超音波が生体組織内（または標準反射体を入れた媒質内
）に送出される。
信号Ｐｄが送信回路２に送られ、超音波送信に必要なだ
けの電力を持ったパルスで超音波振動子３が駆動され、
超音波が生体組織内（または標準反射体を入れた媒質内
）に送出される。
生体組織（または標準反射体）４がらの反射波は再び振
動子３によって受信され、受信回路５によって適当なレ
ベルまで増幅されて（ロ）に示す受信信号Ｖｒとなって
データ収集部６に送られる。
動子３によって受信され、受信回路５によって適当なレ
ベルまで増幅されて（ロ）に示す受信信号Ｖｒとなって
データ収集部６に送られる。
タイミング制御部１がらは、振動子表面から測定したい
反射部位までの距離に対応した時間ＴＩだけＰｄから遅
れたタイミングで、データ収集部６に対して（ハ）に示
すゲートパルスｐｇを送り。
反射部位までの距離に対応した時間ＴＩだけＰｄから遅
れたタイミングで、データ収集部６に対して（ハ）に示
すゲートパルスｐｇを送り。
希望する部位の反射信号をデータとして取り込む。
なお、Ｐｇの幅τは測定したい範囲に対応して定められ
る。収集されたデータは例えば（ニ）に示す如き波形に
対応している。
る。収集されたデータは例えば（ニ）に示す如き波形に
対応している。
収集されたデータは周波数解析部７に送られ。
周波数解析結果がデータ・メモリ８に送られる。
周波数解析結果としては５例えば（Ａ）に示すような生
体組織からの反射波スペクトルや１ （Ｂ）に示すよう
な標準反射体からの反射波スペクトルが得られる。
体組織からの反射波スペクトルや１ （Ｂ）に示すよう
な標準反射体からの反射波スペクトルが得られる。
演算処理部９は、データ・メモリ８に保持されている周
波数解析結果に対して前記のごとき種々の計算を行い、
希望する結果を得る。
波数解析結果に対して前記のごとき種々の計算を行い、
希望する結果を得る。
具体的には、 （Ａ）で示されるスペクトルを（Ｂ）で
示されるスペクトルで正規化し、 （Ｃ）で示される如
きスペクトル（（１）式に対応する）を得て（３）式の
ｒｍを求め、更に（８）式、（９）式の左辺を計算すれ
ばよい。しかる後に回帰計算を行い、ｎβ 、σを決定
すればよい。各深さ２のβ がらβｉも求めることがで
きる。
示されるスペクトルで正規化し、 （Ｃ）で示される如
きスペクトル（（１）式に対応する）を得て（３）式の
ｒｍを求め、更に（８）式、（９）式の左辺を計算すれ
ばよい。しかる後に回帰計算を行い、ｎβ 、σを決定
すればよい。各深さ２のβ がらβｉも求めることがで
きる。
■の構成は上記の演算が実現できる構成であればどの様
なものでもよく２例えばマイクロ・プロセッサ、ＲＡＭ
、ＲＯＭ、Ｉ１０ポート等から成るマイクロコンピュー
タでもよい。
なものでもよく２例えばマイクロ・プロセッサ、ＲＡＭ
、ＲＯＭ、Ｉ１０ポート等から成るマイクロコンピュー
タでもよい。
本発明によれば。
■正規化によって計測系に影響されない生体組織特性関
数を抽出し、さらにある周波数の値で再正規化すること
によって１組織領域境界の不連続透過の影響を除くこと
ができ。
数を抽出し、さらにある周波数の値で再正規化すること
によって１組織領域境界の不連続透過の影響を除くこと
ができ。
■推定関数を非指数関数項で除して対数を採ることによ
り１回帰計算が容易になり。
り１回帰計算が容易になり。
■得られたパラメータから組織特性値を求めることが可
能となる効果を有する。
能となる効果を有する。
第１図は超音波走査方向Ｚ軸におりる生体組織構造の断
面を示す模式図であり、■はトランスデユーサ、■は体
表、ｉは領域ｉ−１と領域ｉとの境界を示し、１は夫々
の領域中の超音波行路長を示す。 第２図は本発明の一実施例装置のブロック図であり、１
はタイミング制御部、２は送信回路、３は超音波振動子
、４は標準反射体、５は受信回路。 ６はデータ収集部２７は周波数解析部、８はデータメモ
リ、９は演算処理部である。 第３図は第２図における各部の波形を示す図。 第４図は周波数特性の関係を示す図である。
面を示す模式図であり、■はトランスデユーサ、■は体
表、ｉは領域ｉ−１と領域ｉとの境界を示し、１は夫々
の領域中の超音波行路長を示す。 第２図は本発明の一実施例装置のブロック図であり、１
はタイミング制御部、２は送信回路、３は超音波振動子
、４は標準反射体、５は受信回路。 ６はデータ収集部２７は周波数解析部、８はデータメモ
リ、９は演算処理部である。 第３図は第２図における各部の波形を示す図。 第４図は周波数特性の関係を示す図である。
Claims (1)
- 【特許請求の範囲】 （１１超音波パルスを生体中に送信し、その反射波を受
信し２その受信信号を分析することにより生体組織特性
を測定する方式において。 少なくとも次のステップを有することを特徴とする超音
波生体組織特性測定方式。 イ）分析すべき部位からの反射波受信信号を周波数分析
し、そのパワー・スペクトルを求めるステップ。 口）該部位に標準反射体を有する場合の送受測定系のパ
ワー・スペクトルで、上記イのステップで求めたパワー
・スペクトルを正規化し、生体の組織伝達関数Ｒｆｆ）
を求めるステップ。 ハ）特定の周波数ｆｏでの値Ｒ（ｆｏ）でＲ（［１を正
規化したＰ ｆｆｌを求めるステップ。 二）一方において、理論及び／又は実験から予測される
生体組織伝達関数を同じ（ｆｏでの値で正規化した関数
Ｑ ｆｆ）が Ｑ（ｆｌ ＝ Ａ （［１＊ ｅβ（チノで表わされ
ることを利用して、上記ハのステップで得られたＰ（「
）をＡ （ｆｌで除し。 Ｐ ｔｆｌ ／ Ａ （ｆｌ を求めるステップ。 ホ）上記二のステップで得た値の対数をとるステップ。 へ）上記ホのステップで得られた周波数の関数としての
値を予測式のＢ ｆｆｌと等値であるとして。 Ａ （ｆｌ 、 Ｂ （ｆｌに含まれる各種のパラメ
ータの最適値を求める回帰ステップ。 （２）上記二のステップにおけるＡ （ｆ）が周波数ｆ
の０乗（ｎは生体組織特性を示す定数）に比例する式で
表わされ、かつＢ （ｆｌがｆの１次式で表わされ、ｆ
の１次の係数が減衰定数の周波数傾斜（減衰傾斜）の線
積分と関連した量であることを特徴とする特許請求の範
囲第１項に記載の生体組織特性測定方式。 （３１１記二のステップにおりるＢ （ｆｌが周波数ｆ
の２次式で表わされ、定数項、１次項の係数、２次項の
係数がそれぞれ、平均的微細構造、減衰（頃斜線積分、
走査方向の自己相関距離に関連した量であることを特徴
とする特許請求の範囲第１項に記載の生体組織特性測定
方式。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP57129902A JPS5919852A (ja) | 1982-07-26 | 1982-07-26 | 超音波生体組織特性測定方式 |
| US06/516,555 US4512195A (en) | 1982-07-26 | 1983-07-25 | Ultrasonic living body tissue characterization method |
| EP83304306A EP0100234B1 (en) | 1982-07-26 | 1983-07-26 | Ultrasonic measurement of characteristic values of a medium |
| DE8383304306T DE3376880D1 (en) | 1982-07-26 | 1983-07-26 | Ultrasonic measurement of characteristic values of a medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP57129902A JPS5919852A (ja) | 1982-07-26 | 1982-07-26 | 超音波生体組織特性測定方式 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPS5919852A true JPS5919852A (ja) | 1984-02-01 |
Family
ID=15021207
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP57129902A Pending JPS5919852A (ja) | 1982-07-26 | 1982-07-26 | 超音波生体組織特性測定方式 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US4512195A (ja) |
| EP (1) | EP0100234B1 (ja) |
| JP (1) | JPS5919852A (ja) |
| DE (1) | DE3376880D1 (ja) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1985002682A1 (en) * | 1983-12-14 | 1985-06-20 | Fujitsu Limited | Apparatus for measuring ultrasonic characteristic values of a medium |
| JPS61226653A (ja) * | 1985-03-29 | 1986-10-08 | エヌ・ベー・フイリツプス・フルーイランペンフアブリケン | 超音波エコーグラフィ検査方法及び装置 |
| JPS62123354A (ja) * | 1985-03-29 | 1987-06-04 | エヌ・ベ−・フイリツプス・フル−イランペンフアブリケン | 超音波エコーグラフィによる物体の検査方法及び装置 |
| JPS62123355A (ja) * | 1985-03-29 | 1987-06-04 | エヌ・ベ−・フイリツプス・フル−イランペンフアブリケン | 媒体の超音波エコーグラフィック検査方法及び装置 |
| WO1993003674A1 (en) * | 1991-08-13 | 1993-03-04 | Daikin Industries, Ltd. | Inspection apparatus |
| JP2011505213A (ja) * | 2007-12-03 | 2011-02-24 | ユッカ・ユルヴェリン | 超音波技法を用いて材料の厚みを測定するための方法 |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS60116345A (ja) * | 1983-11-30 | 1985-06-22 | 富士通株式会社 | 超音波診断装置 |
| FR2556844B1 (fr) * | 1983-12-14 | 1987-11-13 | Labo Electronique Physique | Appareil d'examen de milieux par echographie ultrasonore |
| DE3504210A1 (de) * | 1985-02-07 | 1986-08-07 | Kraftwerk Union AG, 4330 Mülheim | Einrichtung zum bestimmen von oberflaechenrissen |
| US4982339A (en) * | 1985-11-18 | 1991-01-01 | The United States Of America As Represented By Department Of Health And Human Service | High speed texture discriminator for ultrasonic imaging |
| US4817015A (en) * | 1985-11-18 | 1989-03-28 | The United States Government As Represented By The Secretary Of The Health And Human Services | High speed texture discriminator for ultrasonic imaging |
| US4836210A (en) * | 1987-06-15 | 1989-06-06 | Fujitsu Limited | Ultrasonic analyzer |
| FR2619448B1 (fr) * | 1987-08-14 | 1990-01-19 | Edap Int | Procede et dispositif de caracterisation tissulaire par reflexion d'impulsions ultrasonores a large bande de frequences, transposition du spectre de frequence des echos dans une gamme audible et diagnostic par ecoute |
| US4855911A (en) * | 1987-11-16 | 1989-08-08 | Massachusetts Institute Of Technology | Ultrasonic tissue characterization |
| US5417215A (en) * | 1994-02-04 | 1995-05-23 | Long Island Jewish Medical Center | Method of tissue characterization by ultrasound |
| AU1454699A (en) | 1997-11-14 | 1999-06-07 | Colorado Seminary | Ultrasonic system for grading meat |
| US8600497B1 (en) * | 2006-03-31 | 2013-12-03 | Pacesetter, Inc. | Systems and methods to monitor and treat heart failure conditions |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4047172A (en) * | 1972-05-30 | 1977-09-06 | General Electric Company | Signal processor for pulse-echo system providing interference level compensation and feed-forward normalization |
| US3896411A (en) * | 1974-02-19 | 1975-07-22 | Westinghouse Electric Corp | Reverberation condition adaptive sonar receiving system and method |
| JPS57550A (en) * | 1980-06-03 | 1982-01-05 | Fujitsu Ltd | Identification systm for specimen |
| JPS57179745A (en) * | 1981-04-30 | 1982-11-05 | Fujitsu Ltd | Method and device for measuring material property by ultrasonic wave |
| US4470303A (en) * | 1982-09-20 | 1984-09-11 | General Electric Company | Quantitative volume backscatter imaging |
-
1982
- 1982-07-26 JP JP57129902A patent/JPS5919852A/ja active Pending
-
1983
- 1983-07-25 US US06/516,555 patent/US4512195A/en not_active Expired - Fee Related
- 1983-07-26 EP EP83304306A patent/EP0100234B1/en not_active Expired
- 1983-07-26 DE DE8383304306T patent/DE3376880D1/de not_active Expired
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1985002682A1 (en) * | 1983-12-14 | 1985-06-20 | Fujitsu Limited | Apparatus for measuring ultrasonic characteristic values of a medium |
| JPS61226653A (ja) * | 1985-03-29 | 1986-10-08 | エヌ・ベー・フイリツプス・フルーイランペンフアブリケン | 超音波エコーグラフィ検査方法及び装置 |
| JPS62123354A (ja) * | 1985-03-29 | 1987-06-04 | エヌ・ベ−・フイリツプス・フル−イランペンフアブリケン | 超音波エコーグラフィによる物体の検査方法及び装置 |
| JPS62123355A (ja) * | 1985-03-29 | 1987-06-04 | エヌ・ベ−・フイリツプス・フル−イランペンフアブリケン | 媒体の超音波エコーグラフィック検査方法及び装置 |
| WO1993003674A1 (en) * | 1991-08-13 | 1993-03-04 | Daikin Industries, Ltd. | Inspection apparatus |
| US5450849A (en) * | 1991-08-13 | 1995-09-19 | Daikin Industries, Ltd. | Location apparatus |
| JP2011505213A (ja) * | 2007-12-03 | 2011-02-24 | ユッカ・ユルヴェリン | 超音波技法を用いて材料の厚みを測定するための方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| DE3376880D1 (en) | 1988-07-07 |
| EP0100234B1 (en) | 1988-06-01 |
| EP0100234A3 (en) | 1985-05-22 |
| EP0100234A2 (en) | 1984-02-08 |
| US4512195A (en) | 1985-04-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP0123427B1 (en) | Ultrasonic medium characterization | |
| JPS5919852A (ja) | 超音波生体組織特性測定方式 | |
| US6634233B2 (en) | Method for determining the wall thickness and the speed of sound in a tube from reflected and transmitted ultrasound pulses | |
| Saniie et al. | Quantitative grain size evaluation using ultrasonic backscattered echoes | |
| US4754760A (en) | Ultrasonic pulse temperature determination method and apparatus | |
| JP2849827B2 (ja) | パルスドップラ信号のフーリエプロセッシング方法および装置 | |
| US4257278A (en) | Quantitative volume blood flow measurement by an ultrasound imaging system featuring a Doppler modality | |
| EP0008517B1 (en) | Duplex ultrasonic imaging system with repetitive excitation of common transducer in doppler modality | |
| US5010885A (en) | Ultrasonic echograph with controllable phase coherence | |
| CA1130440A (en) | Sonic pulse-echo method and apparatus for determining attenuation coefficients | |
| US5097836A (en) | Untrasound diagnostic equipment for calculating and displaying integrated backscatter or scattering coefficients by using scattering power or scattering power spectrum of blood | |
| Jensen et al. | Ultrasound fields in an attenuating medium | |
| Wilson et al. | Frequency domain processing for ultrasonic attenuation measurement in liver | |
| US4511984A (en) | Ultrasound diagnostic apparatus | |
| JP2006506190A (ja) | 診断信号処理方法及びシステム | |
| JPH0713631B2 (ja) | 媒体の超音波エコーグラフィック検査方法及び装置 | |
| WO2014055410A1 (en) | Shear wave attenuation from k-space analysis system | |
| US6142944A (en) | Doppler motion detection with automatic angle correction | |
| Hoskins et al. | Finite beam-width ray model for geometric spectral broadening | |
| Nowak et al. | A comparison of methods for the determination of sound velocity in biological materials: A case study | |
| US20170156701A1 (en) | System and method for gradient-based k-space search for shear wave velocity dispersion estimation | |
| JP2669204B2 (ja) | 探査装置 | |
| Lacefield et al. | Spatial coherence analysis applied to aberration correction using a two-dimensional array system | |
| US4546772A (en) | Method and means for determining ultrasonic wave attenuation in tissue using phase locked loop | |
| JPS5937766B2 (ja) | インピ−ダンス分布の測定方法 |