JPH09173334A - 超音波診断装置 - Google Patents
超音波診断装置Info
- Publication number
- JPH09173334A JPH09173334A JP7341798A JP34179895A JPH09173334A JP H09173334 A JPH09173334 A JP H09173334A JP 7341798 A JP7341798 A JP 7341798A JP 34179895 A JP34179895 A JP 34179895A JP H09173334 A JPH09173334 A JP H09173334A
- Authority
- JP
- Japan
- Prior art keywords
- spectrum
- ultrasonic
- filter
- diagnostic apparatus
- output
- 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
Classifications
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
(57)【要約】
【課題】本発明の目的は、従来では深層部で充分に確保
することが困難であった信号帯域幅を確保することで距
離分解能を向上させ、画質の良い超音波診断装置を提供
することである。 【解決手段】本発明による超音波診断装置は、超音波プ
ローブ１と、超音波プローブ１を駆動して被検体中に超
音波を送信する送信系２と、超音波の被検体による反射
波を超音波プローブ１で受信して電気信号として取り出
す受信系３と、受信系３から出力される受信信号を深さ
毎に所定のスペクトラム波形に変形するスペクトラム処
理部７と、スペクトラム処理部７の出力について超音波
の画像を生成するための処理を行う画像生成手段４，５
と、画像を表示する表示ユニット６とを具備する。
することが困難であった信号帯域幅を確保することで距
離分解能を向上させ、画質の良い超音波診断装置を提供
することである。 【解決手段】本発明による超音波診断装置は、超音波プ
ローブ１と、超音波プローブ１を駆動して被検体中に超
音波を送信する送信系２と、超音波の被検体による反射
波を超音波プローブ１で受信して電気信号として取り出
す受信系３と、受信系３から出力される受信信号を深さ
毎に所定のスペクトラム波形に変形するスペクトラム処
理部７と、スペクトラム処理部７の出力について超音波
の画像を生成するための処理を行う画像生成手段４，５
と、画像を表示する表示ユニット６とを具備する。
Description
【０００１】
【発明の属する技術分野】本発明は、被検体からの超音
波の受信信号のスペクトラムを変形して画質の向上を図
った超音波診断装置に関する。
波の受信信号のスペクトラムを変形して画質の向上を図
った超音波診断装置に関する。
【０００２】
【従来の技術】近年、超音波の生体内での周波数依存減
衰に伴う受信信号の劣化を補正し、画質の向上を図った
超音波診断装置として、特公昭５９−４１７３１号に記
載されているものがある。
衰に伴う受信信号の劣化を補正し、画質の向上を図った
超音波診断装置として、特公昭５９−４１７３１号に記
載されているものがある。
【０００３】
【発明が解決しようとする課題】一般的に、高空間分解
能であれば高画質が得られるが、空間分解能のうち方位
方向の分解能は中心周波数に比例してよくなり、距離方
向の分解能は帯域幅に比例してよくなることが知られて
いる、ここで、高画質を得るために広帯域駆動された超
音波の受信信号は、生体内での周波数依存減衰により、
被検体の深層部になればなるほどその中心周波数は減少
し、また狭帯域化してしまう（図８参照）。
能であれば高画質が得られるが、空間分解能のうち方位
方向の分解能は中心周波数に比例してよくなり、距離方
向の分解能は帯域幅に比例してよくなることが知られて
いる、ここで、高画質を得るために広帯域駆動された超
音波の受信信号は、生体内での周波数依存減衰により、
被検体の深層部になればなるほどその中心周波数は減少
し、また狭帯域化してしまう（図８参照）。
【０００４】この問題に対して、上記の従来例は、受信
信号の高周波成分を深さに応じて逐次強調するように補
正する手段を有する。ここでこの従来例にも以下の問題
がある。
信号の高周波成分を深さに応じて逐次強調するように補
正する手段を有する。ここでこの従来例にも以下の問題
がある。
【０００５】受信信号の周波数特性は、プローブ単体の
周波数特性と生体の周波数依存減衰特性とによって決ま
る。ここで、生体の周波数依存減衰は深さに従って著し
く変化していること、およびプローブ単体の周波数特性
が必ずしも単純な形状をしていないことにより、超音波
診断装置内で実際に得られる受信信号の周波数特性（以
下、スペクトラムと称する）は、深さやプローブの設定
によって非常に複雑な形状となる。一方、従来例ではこ
のような複雑な変化を受けた受信信号に対して、単純に
高周波成分のみが減衰しているという仮定に基づいて、
その高周波成分を補正しようとしているために、必ずし
も最適な補正の効果が得られるとは限らない。
周波数特性と生体の周波数依存減衰特性とによって決ま
る。ここで、生体の周波数依存減衰は深さに従って著し
く変化していること、およびプローブ単体の周波数特性
が必ずしも単純な形状をしていないことにより、超音波
診断装置内で実際に得られる受信信号の周波数特性（以
下、スペクトラムと称する）は、深さやプローブの設定
によって非常に複雑な形状となる。一方、従来例ではこ
のような複雑な変化を受けた受信信号に対して、単純に
高周波成分のみが減衰しているという仮定に基づいて、
その高周波成分を補正しようとしているために、必ずし
も最適な補正の効果が得られるとは限らない。
【０００６】本発明はこのような従来の超音波生体減衰
補正法の問題点に鑑みて創案されたものである。すなわ
ち、本発明の目的は、生体内で複雑な変化を受けた受信
信号のスペクトラムを、あらかじめ規定した最適なスペ
クトラムになるように変形することで、従来よりも最適
な補正を実施して画質の良い超音波診断装置を提供する
ことにある。本発明の前記ならびにその他の目的と新規
な特徴は、本明細書の記述及び図面によって明らかとな
るであろう。
補正法の問題点に鑑みて創案されたものである。すなわ
ち、本発明の目的は、生体内で複雑な変化を受けた受信
信号のスペクトラムを、あらかじめ規定した最適なスペ
クトラムになるように変形することで、従来よりも最適
な補正を実施して画質の良い超音波診断装置を提供する
ことにある。本発明の前記ならびにその他の目的と新規
な特徴は、本明細書の記述及び図面によって明らかとな
るであろう。
【０００７】
【課題を解決するための手段】本発明による超音波診断
装置は、超音波プローブと、前記超音波プローブを駆動
して被検体中に超音波を送信する送信手段と、前記超音
波の前記被検体による反射波を前記超音波プローブで受
信して電気信号として取り出す受信手段と、前記受信手
段から出力される受信信号を深さ毎に所定のスペクトラ
ム波形に変形するスペクトラム変形手段と、前記スペク
トラム変形手段の出力について超音波の画像を生成する
ための処理を行う画像生成手段と、前記画像を表示する
手段とを具備したことを要部とする。
装置は、超音波プローブと、前記超音波プローブを駆動
して被検体中に超音波を送信する送信手段と、前記超音
波の前記被検体による反射波を前記超音波プローブで受
信して電気信号として取り出す受信手段と、前記受信手
段から出力される受信信号を深さ毎に所定のスペクトラ
ム波形に変形するスペクトラム変形手段と、前記スペク
トラム変形手段の出力について超音波の画像を生成する
ための処理を行う画像生成手段と、前記画像を表示する
手段とを具備したことを要部とする。
【０００８】ここで、スペクトラム変形手段は、深さ毎
にあらかじめ規定された所定のスペクトラム形状になる
ように受信信号スペクトラムを変形させる。従って、単
純に高周波成分を強調する従来に比べて、最適な補正を
実施することが可能となり、深部でも画質が向上され
る。
にあらかじめ規定された所定のスペクトラム形状になる
ように受信信号スペクトラムを変形させる。従って、単
純に高周波成分を強調する従来に比べて、最適な補正を
実施することが可能となり、深部でも画質が向上され
る。
【０００９】
【発明の実施の形態】以下、本発明の実施の形態を図面
を参照して説明する。図１に本実施の形態による超音波
診断装置の構成を示す。システムコントローラー０は本
システム全体の動作に必要とされる全ての制御を司る。
超音波プローブ１は、複数の微小圧電素子が先端に配列
されているアレイ型プローブであり、セクタスキャン
用、リニアスキャン用、他のスキャン用のいずれのタイ
プであってもよい。
を参照して説明する。図１に本実施の形態による超音波
診断装置の構成を示す。システムコントローラー０は本
システム全体の動作に必要とされる全ての制御を司る。
超音波プローブ１は、複数の微小圧電素子が先端に配列
されているアレイ型プローブであり、セクタスキャン
用、リニアスキャン用、他のスキャン用のいずれのタイ
プであってもよい。
【００１０】送信系２は、クロック発生器、レートパル
ス発生器、送信遅延回路、パルサを有する。クロック発
生器はクロックパルスを発生する。レートパルス発生器
は、クロックパルスを分周し、例えば５ＫＨｚのレート
パルスを発生する。このレートパルスはチャンネル数分
に分配され、送信遅延回路に送られる。送信遅延回路
は、チャンネル毎に個々に設けられる。送信遅延回路
は、超音波をビーム状に集束し、且つ超音波ビームを予
定の方向に振るために必要な遅延時間をチャンネル毎に
与える。パルサは、チャンネル毎に個々に設けられる。
パルサは、同じチャンネルに対応する送信遅延回路から
レートパルスをトリガ信号として入力し、高周波の電圧
パルスを対応するチャンネルの圧電素子に出力する。電
圧パルスが供給された圧電素子は、超音波パルスを発生
する。各圧電素子から発生された超音波は、遅延時間の
差異に対応する向きに超音波ビームとして集束する。
ス発生器、送信遅延回路、パルサを有する。クロック発
生器はクロックパルスを発生する。レートパルス発生器
は、クロックパルスを分周し、例えば５ＫＨｚのレート
パルスを発生する。このレートパルスはチャンネル数分
に分配され、送信遅延回路に送られる。送信遅延回路
は、チャンネル毎に個々に設けられる。送信遅延回路
は、超音波をビーム状に集束し、且つ超音波ビームを予
定の方向に振るために必要な遅延時間をチャンネル毎に
与える。パルサは、チャンネル毎に個々に設けられる。
パルサは、同じチャンネルに対応する送信遅延回路から
レートパルスをトリガ信号として入力し、高周波の電圧
パルスを対応するチャンネルの圧電素子に出力する。電
圧パルスが供給された圧電素子は、超音波パルスを発生
する。各圧電素子から発生された超音波は、遅延時間の
差異に対応する向きに超音波ビームとして集束する。
【００１１】超音波は、被検体内の音響インピーダンス
の境界で反射する。反射波は超音波プローブ１で電気信
号に変換される。超音波プローブ１で変換された電気信
号は、受信系３に取り込まれる。受信系３としては、特
公平６−１４９３４号に記載される構成例のようなディ
ジタルシームフォーマが好適であり、この場合、ディジ
タル化（深さ方向にサンプリング）された受信信号を出
力する。受信系３は、プリアンプ、受信遅延回路、加算
器を含む。超音波プローブ１で変換された電気信号は、
プリアンプを介して受信遅延回路に送り込まれる。受信
遅延回路では受信波を集束し、且つ予定の方向からの反
射波成分を強調するために必要な遅延時間がチャンネル
毎に与えられる。遅延時間が与えられた電気信号は加算
器で加算される。これにより受信信号が生成される。
の境界で反射する。反射波は超音波プローブ１で電気信
号に変換される。超音波プローブ１で変換された電気信
号は、受信系３に取り込まれる。受信系３としては、特
公平６−１４９３４号に記載される構成例のようなディ
ジタルシームフォーマが好適であり、この場合、ディジ
タル化（深さ方向にサンプリング）された受信信号を出
力する。受信系３は、プリアンプ、受信遅延回路、加算
器を含む。超音波プローブ１で変換された電気信号は、
プリアンプを介して受信遅延回路に送り込まれる。受信
遅延回路では受信波を集束し、且つ予定の方向からの反
射波成分を強調するために必要な遅延時間がチャンネル
毎に与えられる。遅延時間が与えられた電気信号は加算
器で加算される。これにより受信信号が生成される。
【００１２】受信系３で生成された受信信号は、スペク
トラム処理制御部８による制御を受けるスペクトラム処
理部７を介して、検波器４に送られる。検波器４は、ス
ペクトラム処理部７からの出力信号を検波する。検波器
４から出力される検波信号は、ログ圧縮器５を介して表
示ユニット６に送られ、組織断層像（Ｂモード画像）と
して表示される。
トラム処理制御部８による制御を受けるスペクトラム処
理部７を介して、検波器４に送られる。検波器４は、ス
ペクトラム処理部７からの出力信号を検波する。検波器
４から出力される検波信号は、ログ圧縮器５を介して表
示ユニット６に送られ、組織断層像（Ｂモード画像）と
して表示される。
【００１３】スペクトラム処理部７と、スペクトラム処
理制御部８とは、本発明の特徴的な構成要素であり、以
下に、これらの詳細を説明する。ここで、以下で取り扱
うスペクトラム特性について簡単にふれておく。図７に
示すように、縦軸が振幅を示し、各々の最大値で規格化
した［ｄＢ］表示となっている。横軸は周波数を示す。
また、ｄは深さを示す。図２（ａ）は、スペクトラム処
理部７の１つの構成例を示す。この構成の場合、一般的
なＦＩＲ(finite impulse responcs) 型ディジタルフィ
ルタ７−１によってスペクトラム形状の変形を実現する
のが好適である。以下、このフィルタをエコーフィルタ
と称する。
理制御部８とは、本発明の特徴的な構成要素であり、以
下に、これらの詳細を説明する。ここで、以下で取り扱
うスペクトラム特性について簡単にふれておく。図７に
示すように、縦軸が振幅を示し、各々の最大値で規格化
した［ｄＢ］表示となっている。横軸は周波数を示す。
また、ｄは深さを示す。図２（ａ）は、スペクトラム処
理部７の１つの構成例を示す。この構成の場合、一般的
なＦＩＲ(finite impulse responcs) 型ディジタルフィ
ルタ７−１によってスペクトラム形状の変形を実現する
のが好適である。以下、このフィルタをエコーフィルタ
と称する。
【００１４】エコーフィルタ７−１は、深さ方向にサン
プリングされた受信信号ｆ（ｄ）を順次ラッチするＮ−
１個のデータラッチ７−１−Ｄ1 ，７−１−Ｄ2 ，・・
・７−１−ＤN-1 と、Ｎ個の係数乗算器７−１−Ｍ1 ，
７−１−Ｍ2 ，・・・７−１−ＭN と、加算機７−１−
Ａとから構成される。
プリングされた受信信号ｆ（ｄ）を順次ラッチするＮ−
１個のデータラッチ７−１−Ｄ1 ，７−１−Ｄ2 ，・・
・７−１−ＤN-1 と、Ｎ個の係数乗算器７−１−Ｍ1 ，
７−１−Ｍ2 ，・・・７−１−ＭN と、加算機７−１−
Ａとから構成される。
【００１５】本エコーフィルタの構成の場合、スペクト
ラム処理制御部８は、図２（ｂ）に示すような係数記憶
部８−１として実現するのが好適である。係数記憶部８
−１は、ＲＯＭやＲＡＭといった記憶装置で構成され、
システムコントローラー０の制御に基づき、受信信号の
深さｄに応じて、各深さ毎にあらかじめ記憶されていた
任意のＮ個の係数列ｈ１（ｋ，ｄ）を発生する。ここ
で、システムコントローラー０は、特に図示しないが、
例えば、操作卓上に設置された係数制御スイッチ乃至プ
ローブ情報テーブル（超音波プローブの種類を表す情報
のテーブル）を記憶したＲＯＭといった係数列変更手段
と接続されており、これらの係数列変更手段の設定に対
応して上述の係数記憶部８−１の制御を実施するように
する。すなわち、係数制御スイッチの設定や、動作可能
となっている超音波プローブの種類に対応して、係数記
憶部８−１に記憶されていた任意のＮ個の係数列ｈｊが
ｈ２（ｋ，ｄ），ｈ３（ｋ，ｄ），…ｈ１（ｋ，ｄ）の
中から選択される。これらの係数列の特性は、以下に示
すフィルタ特性算出法によって決定されている。
ラム処理制御部８は、図２（ｂ）に示すような係数記憶
部８−１として実現するのが好適である。係数記憶部８
−１は、ＲＯＭやＲＡＭといった記憶装置で構成され、
システムコントローラー０の制御に基づき、受信信号の
深さｄに応じて、各深さ毎にあらかじめ記憶されていた
任意のＮ個の係数列ｈ１（ｋ，ｄ）を発生する。ここ
で、システムコントローラー０は、特に図示しないが、
例えば、操作卓上に設置された係数制御スイッチ乃至プ
ローブ情報テーブル（超音波プローブの種類を表す情報
のテーブル）を記憶したＲＯＭといった係数列変更手段
と接続されており、これらの係数列変更手段の設定に対
応して上述の係数記憶部８−１の制御を実施するように
する。すなわち、係数制御スイッチの設定や、動作可能
となっている超音波プローブの種類に対応して、係数記
憶部８−１に記憶されていた任意のＮ個の係数列ｈｊが
ｈ２（ｋ，ｄ），ｈ３（ｋ，ｄ），…ｈ１（ｋ，ｄ）の
中から選択される。これらの係数列の特性は、以下に示
すフィルタ特性算出法によって決定されている。
【００１６】以上のような構成で、（１）式に示した係
数の畳み込み演算を実施する（図９参照）。これにより
本構成では、時間軸上の演算によって所望のスペクトラ
ム変形を受けた信号ｇが得られる。
数の畳み込み演算を実施する（図９参照）。これにより
本構成では、時間軸上の演算によって所望のスペクトラ
ム変形を受けた信号ｇが得られる。
【００１７】
【数１】
【００１８】ここで、最適な所定のスペクトラム波形に
変形するためのフィルタ特性算出法の一例について以下
に説明する。基本的な考え方は、従来の「ある決まった
タイプ（ここでいうタイプとは低域阻止域とか帯域通過
域、ないし従来例に記した高域強調型といったものをさ
す）のフィルタ特性を与え、結果として信号のスペクト
ラムが得られる」というものではなく、「フィルタの出
力において得られるスペクトラムを仕様として与え、そ
のスペクトラムを得るためのフィルタ特性を、出力スペ
クトラム期待値と、入力スペクトラムとの両方を用いて
決定する」というものである。具体例としては、式
（２）に示したように出力スペクトラム期待値と入力ス
ペクトラムの比としてフィルタ特性を算出する。
変形するためのフィルタ特性算出法の一例について以下
に説明する。基本的な考え方は、従来の「ある決まった
タイプ（ここでいうタイプとは低域阻止域とか帯域通過
域、ないし従来例に記した高域強調型といったものをさ
す）のフィルタ特性を与え、結果として信号のスペクト
ラムが得られる」というものではなく、「フィルタの出
力において得られるスペクトラムを仕様として与え、そ
のスペクトラムを得るためのフィルタ特性を、出力スペ
クトラム期待値と、入力スペクトラムとの両方を用いて
決定する」というものである。具体例としては、式
（２）に示したように出力スペクトラム期待値と入力ス
ペクトラムの比としてフィルタ特性を算出する。
【００１９】
【数２】
【００２０】入力スペクトラムとして、各々周波数特性
の異なる超音波プローブ毎に、平均的な被検体の減衰特
性を与えた時に各深さ毎に得られるスペクトラムを用い
るのが好適である。このような入力スペクトラムは、例
えば実際に平均的な被検体や、被検体を模擬したファン
トムから受信信号を得て規定する方法が考えられる。ま
た、超音波プローブの周波数特性と生体減衰の深さ応答
を式で与えて、計算機で演算して規定する方法も考えら
れる。
の異なる超音波プローブ毎に、平均的な被検体の減衰特
性を与えた時に各深さ毎に得られるスペクトラムを用い
るのが好適である。このような入力スペクトラムは、例
えば実際に平均的な被検体や、被検体を模擬したファン
トムから受信信号を得て規定する方法が考えられる。ま
た、超音波プローブの周波数特性と生体減衰の深さ応答
を式で与えて、計算機で演算して規定する方法も考えら
れる。
【００２１】本発明では、フィルタの出力において得ら
れるスペクトラムを仕様として与えるので、この出力ス
ペクトラム期待値の仕様が、得られる画質を左右する。
これは期待される出力スペクトラムを上手に与えること
で、最適な画質を得ることが可能であることを意味して
いる。本発明では、最適な所定のスペクトラム期待値の
一例として、（３）式で与えられるガウス関数を用いる
ことを特徴とする。ガウス関数は、（４）式に示すよう
に、周波数軸においても時間軸においてもガウス型の形
状が維持されることが知られている。
れるスペクトラムを仕様として与えるので、この出力ス
ペクトラム期待値の仕様が、得られる画質を左右する。
これは期待される出力スペクトラムを上手に与えること
で、最適な画質を得ることが可能であることを意味して
いる。本発明では、最適な所定のスペクトラム期待値の
一例として、（３）式で与えられるガウス関数を用いる
ことを特徴とする。ガウス関数は、（４）式に示すよう
に、周波数軸においても時間軸においてもガウス型の形
状が維持されることが知られている。
【００２２】
【数３】
【００２３】従って、フィルタの出力において得られる
スペクトラムをガウス型にすれば、フィルタ出力結果で
得られる時間特性もガウス型となり、フィルタのインパ
ルス応答によって時間特性すなわち距離分解能が劣化す
ることが軽減されるので、高画質化の観点から好適な関
数となり得る。ガウス関数を用いた出力スペクトラム期
待値の一例を、図１０に示した。
スペクトラムをガウス型にすれば、フィルタ出力結果で
得られる時間特性もガウス型となり、フィルタのインパ
ルス応答によって時間特性すなわち距離分解能が劣化す
ることが軽減されるので、高画質化の観点から好適な関
数となり得る。ガウス関数を用いた出力スペクトラム期
待値の一例を、図１０に示した。
【００２４】ここで、超音波プローブの周波数特性が必
ずしも単峰性でないことや、生体減衰の分散性のため
に、上述のようにして決定されるフィルタ特性は、一般
的な単峰性を有することはまれで、複雑な多峰性を有す
る特性が深さ毎に刻々と変化することも本発明の特徴で
ある。図８に示す受信信号スペクトラムと、図１０に示
した出力スペクトラム期待値とから得られたフィルタの
スペクトラム特性の一例を図１１に示す。
ずしも単峰性でないことや、生体減衰の分散性のため
に、上述のようにして決定されるフィルタ特性は、一般
的な単峰性を有することはまれで、複雑な多峰性を有す
る特性が深さ毎に刻々と変化することも本発明の特徴で
ある。図８に示す受信信号スペクトラムと、図１０に示
した出力スペクトラム期待値とから得られたフィルタの
スペクトラム特性の一例を図１１に示す。
【００２５】また、このような複雑な周波数特性をフィ
ルタで実現するためには、アナログフィルタで構成する
のは望ましくなく、ディジタルフィルタで実現すること
もまた、本発明の特徴である。
ルタで実現するためには、アナログフィルタで構成する
のは望ましくなく、ディジタルフィルタで実現すること
もまた、本発明の特徴である。
【００２６】フィルタをＦＩＲ型で構成するのが好適で
ある理由は、フィルタ特性を深さ毎に刻々と変化させた
際のフィルタの過渡応答が、ＩＩＲ型よりも小さく抑え
ることが容易なためである。
ある理由は、フィルタ特性を深さ毎に刻々と変化させた
際のフィルタの過渡応答が、ＩＩＲ型よりも小さく抑え
ることが容易なためである。
【００２７】次に、（２）式で算出されたフィルタの周
波数特性Ｆ（ｗ，ｄ）（ｗ＝１，２，…Ｍ）から、Ｎ個
の係数例ｈ（ｋ，ｄ）を得るには、Ｆ（ｗ，ｄ）を逆フ
ーリエ変換した後、必要なＮ個のインパルス系列を係数
列として取り出せば良い。但し、この方法で得られるイ
ンパルス系列は、Ｆ（ｗ，ｄ）が中心対称となるとは限
らないため、直線位相特性が得られにくい。そこで、直
線位相特性を確保して、波形歪みを極力軽減するため
に、Ｆ（ｗ，ｄ）の中心から片側の周波数特性のみが有
効となるように設定し、この特性を中心対称に折り返し
たＦ’（ｗ，ｄ）を逆フーリエ変換した後、必要なＮ個
のインパルス系列を係数列として取り出す方法も考えら
れる。図８の受信信号スペクトラムに対して、このよう
にして得られた出力スペクトラムの一例を、受信信号ス
ペクトラムと共に図１２（ａ），（ｂ）に示した。
波数特性Ｆ（ｗ，ｄ）（ｗ＝１，２，…Ｍ）から、Ｎ個
の係数例ｈ（ｋ，ｄ）を得るには、Ｆ（ｗ，ｄ）を逆フ
ーリエ変換した後、必要なＮ個のインパルス系列を係数
列として取り出せば良い。但し、この方法で得られるイ
ンパルス系列は、Ｆ（ｗ，ｄ）が中心対称となるとは限
らないため、直線位相特性が得られにくい。そこで、直
線位相特性を確保して、波形歪みを極力軽減するため
に、Ｆ（ｗ，ｄ）の中心から片側の周波数特性のみが有
効となるように設定し、この特性を中心対称に折り返し
たＦ’（ｗ，ｄ）を逆フーリエ変換した後、必要なＮ個
のインパルス系列を係数列として取り出す方法も考えら
れる。図８の受信信号スペクトラムに対して、このよう
にして得られた出力スペクトラムの一例を、受信信号ス
ペクトラムと共に図１２（ａ），（ｂ）に示した。
【００２８】もちろん、このような前処理をせずに、Ｆ
（ｗ，ｄ）をそのまま逆フーリエ変換した後に取り出し
たＮ個のインパルス系列を用いてフィルタ処理を実施す
る方法も考えられる。このようにして得られる係数列
は、一般に複素数となるので、図６に示すような複素型
のフィルタ構成にすれば、本方法の実現が可能である。
上記のＦ’（ｗ，ｄ）を用いた係数列では、中心から片
側しか出力スペクトラム期待値と合致しないが、本複素
型のフィルタ構成では両側の出力スペクトラム期待特性
を再現するので、信号帯域外の不要なノイズ成分をより
多く除去することが出来、Ｓ／Ｎの面で有利となる。所
望のスペクトラム変形を受けた信号ｇ（ｄ）を（５）式
に示す。
（ｗ，ｄ）をそのまま逆フーリエ変換した後に取り出し
たＮ個のインパルス系列を用いてフィルタ処理を実施す
る方法も考えられる。このようにして得られる係数列
は、一般に複素数となるので、図６に示すような複素型
のフィルタ構成にすれば、本方法の実現が可能である。
上記のＦ’（ｗ，ｄ）を用いた係数列では、中心から片
側しか出力スペクトラム期待値と合致しないが、本複素
型のフィルタ構成では両側の出力スペクトラム期待特性
を再現するので、信号帯域外の不要なノイズ成分をより
多く除去することが出来、Ｓ／Ｎの面で有利となる。所
望のスペクトラム変形を受けた信号ｇ（ｄ）を（５）式
に示す。
【００２９】
【数４】
【００３０】以上の構成例では、時間軸の演算によって
スペクトラムの変形を実施しているが、これを周波数軸
上で実施してもよい。以下に、周波数軸上での構成例を
示す。
スペクトラムの変形を実施しているが、これを周波数軸
上で実施してもよい。以下に、周波数軸上での構成例を
示す。
【００３１】図３（ａ）は、スペクトラム処理部７の他
の構成例を示す。この構成の場合、周波数イコライザ７
−２によってスペクトラム形状の変形を実現する。以
下、この周波数イコライザを、エコーイコライザと称す
る。
の構成例を示す。この構成の場合、周波数イコライザ７
−２によってスペクトラム形状の変形を実現する。以
下、この周波数イコライザを、エコーイコライザと称す
る。
【００３２】エコーイコライザ７−２は、深さ方向にサ
ンプリングされた受信信号ｆ（ｄ）を一時記憶する信号
メモリ７−２−Ｍ1 と、この信号メモリ７−２−Ｍ1 に
記憶された受信信号ｆ（ｄ）をＭ個の塊毎に順次読み出
して、フーリエ変換処理を行うＦＦＴ処理部７−２−Ｆ
と、ＦＦＴ処理結果のＦ（ｗ，ｄ）に対してＭ個の各周
波数に応じて必要なゲインを与えるイコライザ７−２−
Ｅと、イコライザ処理結果のＭ個の周波数データＧ
（ｗ，ｄ）を一時記憶する信号メモリ７−２−Ｍ2と、
この信号メモリ７−２−Ｍ2 に記憶された周波数データ
Ｇ（ｗ，ｄ）をＭ個の塊毎に順次読み出して、逆フーリ
エ変換処理を行う逆ＦＦＴ処理部７−２−Ｉとから構成
される。尚、本エコーイコライザ７−２の構成は、演算
の高速化のために、ＡＳＩＣで実現されるのが好適であ
る。
ンプリングされた受信信号ｆ（ｄ）を一時記憶する信号
メモリ７−２−Ｍ1 と、この信号メモリ７−２−Ｍ1 に
記憶された受信信号ｆ（ｄ）をＭ個の塊毎に順次読み出
して、フーリエ変換処理を行うＦＦＴ処理部７−２−Ｆ
と、ＦＦＴ処理結果のＦ（ｗ，ｄ）に対してＭ個の各周
波数に応じて必要なゲインを与えるイコライザ７−２−
Ｅと、イコライザ処理結果のＭ個の周波数データＧ
（ｗ，ｄ）を一時記憶する信号メモリ７−２−Ｍ2と、
この信号メモリ７−２−Ｍ2 に記憶された周波数データ
Ｇ（ｗ，ｄ）をＭ個の塊毎に順次読み出して、逆フーリ
エ変換処理を行う逆ＦＦＴ処理部７−２−Ｉとから構成
される。尚、本エコーイコライザ７−２の構成は、演算
の高速化のために、ＡＳＩＣで実現されるのが好適であ
る。
【００３３】このようなエコーイコライザ構成の場合、
スペクトラム処理制御部８は、図３（ｂ）に示すような
周波数特性記憶部８−２として実現するのが好適であ
る。周波数特性記憶部８−２は、ＲＯＭやＲＡＭといっ
た記憶装置で構成され、システムコントローラー０の制
御に基づき、受信信号の深さｄに応じて、各深さ毎にあ
らかじめ記憶されていた任意のＭ個の周波数特性列Ｈ１
（ｗ，ｄ）を発生する。
スペクトラム処理制御部８は、図３（ｂ）に示すような
周波数特性記憶部８−２として実現するのが好適であ
る。周波数特性記憶部８−２は、ＲＯＭやＲＡＭといっ
た記憶装置で構成され、システムコントローラー０の制
御に基づき、受信信号の深さｄに応じて、各深さ毎にあ
らかじめ記憶されていた任意のＭ個の周波数特性列Ｈ１
（ｗ，ｄ）を発生する。
【００３４】ここで、上述と同様の制御を実施すること
により、係数制御スイッチの設定や、動作可能となって
いる超音波プローブの種別に対応して周波数特性記憶部
８−２に記憶されていた任意のＭ個の周波数特性列Ｈｊ
が、Ｈ２（ｗ，ｄ），Ｈ３（ｗ，ｄ），…Ｈ１（ｗ，
ｄ）の中から選択される。
により、係数制御スイッチの設定や、動作可能となって
いる超音波プローブの種別に対応して周波数特性記憶部
８−２に記憶されていた任意のＭ個の周波数特性列Ｈｊ
が、Ｈ２（ｗ，ｄ），Ｈ３（ｗ，ｄ），…Ｈ１（ｗ，
ｄ）の中から選択される。
【００３５】以上のような構成で、（６）式に示した周
波数特性の乗算を実施する。これにより本構成では、周
波数軸上の演算によって所望のスペクトラム変形を受け
た信号ｇ（ｄ）が得られる。
波数特性の乗算を実施する。これにより本構成では、周
波数軸上の演算によって所望のスペクトラム変形を受け
た信号ｇ（ｄ）が得られる。
【００３６】
【数５】
【００３７】本構成においても、上述の時間軸での構成
で説明したのと基本的に同じ原理と特徴を持って、スペ
クトラムの変形が達成される。上述の時間軸での構成の
説明における「フィルタ」を、「イコライザ」に読み替
えれば理解される。本構成は、上記時間軸の構成よりも
複雑となるが、スペクトラム上で直接的に変形を行うこ
とが出来ることが上記時間軸の構成に対する利点であ
る。
で説明したのと基本的に同じ原理と特徴を持って、スペ
クトラムの変形が達成される。上述の時間軸での構成の
説明における「フィルタ」を、「イコライザ」に読み替
えれば理解される。本構成は、上記時間軸の構成よりも
複雑となるが、スペクトラム上で直接的に変形を行うこ
とが出来ることが上記時間軸の構成に対する利点であ
る。
【００３８】これまでの構成は、基本的にあらかじめ規
定した平均的な受信信号スペクトラムによって得られた
スペクトラム変形特性を用いている。従って、個体差に
よって平均からずれた被検体の場合には、出力スペクト
ラムも最適ではなくなる。この問題の対策の一つとし
て、同一の超音波プローブにおいて複数のスペクトラム
変形特性を有し、上述の実施の形態で説明したような、
操作卓に設置された係数制御スイッチを用いて、検査者
が最適な変形特性を選択する方法が考えられる。
定した平均的な受信信号スペクトラムによって得られた
スペクトラム変形特性を用いている。従って、個体差に
よって平均からずれた被検体の場合には、出力スペクト
ラムも最適ではなくなる。この問題の対策の一つとし
て、同一の超音波プローブにおいて複数のスペクトラム
変形特性を有し、上述の実施の形態で説明したような、
操作卓に設置された係数制御スイッチを用いて、検査者
が最適な変形特性を選択する方法が考えられる。
【００３９】また、別の対策としては、適応型の処理に
よって逐次受信信号のスペクトラムを演算し、各被検体
に応じて最適な変形特性を自動的に与える構成も考えら
れる。本構成は、上述したエコーフィルタ構成や、エコ
ーイコライザ構成のどちらでも実施可能であるが、以下
では、エコーイコライザ構成を例にとって説明する。
よって逐次受信信号のスペクトラムを演算し、各被検体
に応じて最適な変形特性を自動的に与える構成も考えら
れる。本構成は、上述したエコーフィルタ構成や、エコ
ーイコライザ構成のどちらでも実施可能であるが、以下
では、エコーイコライザ構成を例にとって説明する。
【００４０】図４（ａ）は、適応型処理を実施するため
のスペクトラム処理部７の別の構成を示す。基本的構成
と動作は図３で説明したものと重複する部分は省略す
る。本構成で特徴となるのは、周波数特性演算部７−３
と、周波数特性一時記憶部８−３であり、この動作を説
明する。図４（ｂ）は、周波数特性記憶部８−２の構成
を示す。
のスペクトラム処理部７の別の構成を示す。基本的構成
と動作は図３で説明したものと重複する部分は省略す
る。本構成で特徴となるのは、周波数特性演算部７−３
と、周波数特性一時記憶部８−３であり、この動作を説
明する。図４（ｂ）は、周波数特性記憶部８−２の構成
を示す。
【００４１】周波数特性演算部７−３は、Ｆ（ｗ，ｄ）
と、Ｈｊ（ｗ，ｄ）とを入力とし、最適なスペクトラム
変形特性Ｈ’（ｗ，ｄ）を（８）式に従って演算し出力
する。演算されたスペクトラム変形特性Ｈ’（ｗ，ｄ）
を（７）式に示す。
と、Ｈｊ（ｗ，ｄ）とを入力とし、最適なスペクトラム
変形特性Ｈ’（ｗ，ｄ）を（８）式に従って演算し出力
する。演算されたスペクトラム変形特性Ｈ’（ｗ，ｄ）
を（７）式に示す。
【００４２】
【数６】
【００４３】尚、本構成では、周波数特性記憶部８−２
に記憶されているＨｊ（ｗ，ｄ）が、出力スペクトラム
期待値であることが、これまでの説明と異なる。 Ｈ’（ｗ，ｄ）＝Ｈｊ（ｗ，ｄ）／Ｆ（ｗ，ｄ） …（８） 周波数特性演算部７−３で演算されたスペクトラム変形
特性Ｈ’（ｗ，ｄ）は、周波数特性一時記憶部８−３に
一時記憶され、イコライザ７−２−Ｅは周波数特性一時
記憶部８−３から最適なスペクトラム変形特性である
Ｈ’（ｗ，ｄ）を読み出して、スペクトラムの変形を実
施する。ここで、必要に応じて（被検体や対象部位の変
化があった時が好適）、周波数特性演算部７−３での演
算を逐次実施して、周波数特性記憶部８−２に随時記憶
し直せば、適応型の処理が実現される。
に記憶されているＨｊ（ｗ，ｄ）が、出力スペクトラム
期待値であることが、これまでの説明と異なる。 Ｈ’（ｗ，ｄ）＝Ｈｊ（ｗ，ｄ）／Ｆ（ｗ，ｄ） …（８） 周波数特性演算部７−３で演算されたスペクトラム変形
特性Ｈ’（ｗ，ｄ）は、周波数特性一時記憶部８−３に
一時記憶され、イコライザ７−２−Ｅは周波数特性一時
記憶部８−３から最適なスペクトラム変形特性である
Ｈ’（ｗ，ｄ）を読み出して、スペクトラムの変形を実
施する。ここで、必要に応じて（被検体や対象部位の変
化があった時が好適）、周波数特性演算部７−３での演
算を逐次実施して、周波数特性記憶部８−２に随時記憶
し直せば、適応型の処理が実現される。
【００４４】以上説明したように本発明によれば、各深
さにおいて受信信号を最適なスペクトラムに変形して空
間分解能の向上をはかったので、各深さで解像度の優れ
た最適な画像が得られる。従って、従来から画質劣化が
大きかった深部領域はもちろん、浅部領域においても従
来より解像度の優れた画像を得ることが可能となり、有
効性の高い診断結果を得ることができる超音波診断装置
を提供することが出来る。
さにおいて受信信号を最適なスペクトラムに変形して空
間分解能の向上をはかったので、各深さで解像度の優れ
た最適な画像が得られる。従って、従来から画質劣化が
大きかった深部領域はもちろん、浅部領域においても従
来より解像度の優れた画像を得ることが可能となり、有
効性の高い診断結果を得ることができる超音波診断装置
を提供することが出来る。
【００４５】本発明の構成は上述してきたような実施の
形態だけにとどまらず、その要旨を逸脱しない範囲にお
いて様々な変形例が得られることは言うまでもない。例
えば従来から用いられているアナログビームフォーマを
有する超音波診断装置にも、本発明は適応される。この
場合、上述の実施の形態では図１の受信系３に含まれて
いたＡ／Ｄ変換器が、後段のブロックに含まれるため、
スペクトラム処理部７の入力はアナログの信号となる。
従って、スペクトラム処理部７の入力側にＡ／Ｄ変換器
を具備し、スペクトラム処理部７として上記で説明した
種々の構成のいずれかを用いることで本発明が実現され
る。本変形例における超音波診断装置の２種の構成例
を、Ａ／Ｄ変換器の位置に着目して、図５（ａ），
（ｂ）に示す。
形態だけにとどまらず、その要旨を逸脱しない範囲にお
いて様々な変形例が得られることは言うまでもない。例
えば従来から用いられているアナログビームフォーマを
有する超音波診断装置にも、本発明は適応される。この
場合、上述の実施の形態では図１の受信系３に含まれて
いたＡ／Ｄ変換器が、後段のブロックに含まれるため、
スペクトラム処理部７の入力はアナログの信号となる。
従って、スペクトラム処理部７の入力側にＡ／Ｄ変換器
を具備し、スペクトラム処理部７として上記で説明した
種々の構成のいずれかを用いることで本発明が実現され
る。本変形例における超音波診断装置の２種の構成例
を、Ａ／Ｄ変換器の位置に着目して、図５（ａ），
（ｂ）に示す。
【００４６】
【発明の効果】本発明による超音波診断装置は、超音波
プローブと、前記超音波プローブを駆動して被検体中に
超音波を送信する送信手段と、前記超音波の前記被検体
による反射波を前記超音波プローブで受信して電気信号
として取り出す受信手段と、前記受信手段から出力され
る受信信号を深さ毎に所定のスペクトラム波形に変形す
るスペクトラム変形手段と、前記スペクトラム変形手段
の出力について超音波の画像を生成するための処理を行
う画像生成手段と、前記画像を表示する手段とを具備し
たことを要部とする。
プローブと、前記超音波プローブを駆動して被検体中に
超音波を送信する送信手段と、前記超音波の前記被検体
による反射波を前記超音波プローブで受信して電気信号
として取り出す受信手段と、前記受信手段から出力され
る受信信号を深さ毎に所定のスペクトラム波形に変形す
るスペクトラム変形手段と、前記スペクトラム変形手段
の出力について超音波の画像を生成するための処理を行
う画像生成手段と、前記画像を表示する手段とを具備し
たことを要部とする。
【００４７】ここで、スペクトラム変形手段は、深さ毎
にあらかじめ規定された所定のスペクトラム形状になる
ように受信信号スペクトラムを変形させる。従って、単
純に高周波成分を強調する従来に比べて、最適な補正を
実施することが可能となり、深部でも画質が向上され
る。
にあらかじめ規定された所定のスペクトラム形状になる
ように受信信号スペクトラムを変形させる。従って、単
純に高周波成分を強調する従来に比べて、最適な補正を
実施することが可能となり、深部でも画質が向上され
る。
【図１】本発明の一実施の形態に係る超音波診断装置の
構成を示す図。
構成を示す図。
【図２】図１のスぺクトラム処理部の一構成例を示す
図。
図。
【図３】図１のスペクトラム処理部の他の構成を示す
図。
図。
【図４】図１のスペクトラム処理部のさらに他の構成を
示す図。
示す図。
【図５】本発明の他の実施の形態に係る超音波診断装置
の構成を示す図。
の構成を示す図。
【図６】複索型のフィルタの一構成例を示す図。
【図７】スペクトラム特性の説明図。
【図８】代表的な深さでの受信信号スぺクトラム特性例
を示す図。
を示す図。
【図９】本発明におけるスペクトル変形手段の特性算出
法の一例の説明図。
法の一例の説明図。
【図１０】本発明における出力スペクトラムの期待値の
一例を示す図。
一例を示す図。
【図１１】代表的な深さでのスペクトル変形手段の特性
の一例を示す図。
の一例を示す図。
【図１２】出力スペクトラムの一例を受信信号スペクト
ラムと共に示す図。
ラムと共に示す図。
１…プローブ、 ２…送信系、 ３…受信系、 ４…検波器、 ５…ログ圧縮器、 ６…表示ユニット、 ７…スペクトラム処理部、 ８…スペクトラム処理制御部、 ０…システムコントローラー、 ７−１…エコーフィルタ部、 ７−２…エコーイコライザ部、 ７−３…周波数特性演算部、 ８−１…係数記憶部、 ８−２…周波数特性記憶部、 ８−３…周波数特性一時記憶部。
─────────────────────────────────────────────────────
【手続補正書】
【提出日】平成８年８月１２日
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】００１５
【補正方法】変更
【補正内容】
【００１５】本エコーフィルタの構成の場合、スペクト
ラム処理制御部８は、図２（ｂ）に示すような係数記憶
部８−１として実現するのが好適である。係数記憶部８
−１は、ＲＯＭやＲＡＭといった記憶装置で構成され、
システムコントローラー０の制御に基づき、受信信号の
深さｄに応じて、各深さ毎にあらかじめ記憶されていた
任意のＮ個の係数列ｈｊ（ｋ，ｄ）を発生する。ここ
で、システムコントローラー０は、特に図示しないが、
例えば、操作卓上に設置された係数制御スイッチ乃至プ
ローブ情報テーブル（超音波プローブの種類を表す情報
のテーブル）を記憶したＲＯＭといった係数列変更手段
と接続されており、これらの係数列変更手段の設定に対
応して上述の係数記憶部８−１の制御を実施するように
する。すなわち、係数制御スイッチの設定や、動作可能
となっている超音波プローブの種類に対応して、係数記
憶部８−１に記憶されていた任意のＮ個の係数列ｈｊが
ｈ１（ｋ，ｄ），ｈ２（ｋ，ｄ），…ｈｌ（ｋ，ｄ）の
中から選択される。これらの係数列の特性は、以下に示
すフィルタ特性算出法によって決定されている。
ラム処理制御部８は、図２（ｂ）に示すような係数記憶
部８−１として実現するのが好適である。係数記憶部８
−１は、ＲＯＭやＲＡＭといった記憶装置で構成され、
システムコントローラー０の制御に基づき、受信信号の
深さｄに応じて、各深さ毎にあらかじめ記憶されていた
任意のＮ個の係数列ｈｊ（ｋ，ｄ）を発生する。ここ
で、システムコントローラー０は、特に図示しないが、
例えば、操作卓上に設置された係数制御スイッチ乃至プ
ローブ情報テーブル（超音波プローブの種類を表す情報
のテーブル）を記憶したＲＯＭといった係数列変更手段
と接続されており、これらの係数列変更手段の設定に対
応して上述の係数記憶部８−１の制御を実施するように
する。すなわち、係数制御スイッチの設定や、動作可能
となっている超音波プローブの種類に対応して、係数記
憶部８−１に記憶されていた任意のＮ個の係数列ｈｊが
ｈ１（ｋ，ｄ），ｈ２（ｋ，ｄ），…ｈｌ（ｋ，ｄ）の
中から選択される。これらの係数列の特性は、以下に示
すフィルタ特性算出法によって決定されている。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００３３
【補正方法】変更
【補正内容】
【００３３】このようなエコーイコライザ構成の場合、
スペクトラム処理制御部８は、図３（ｂ）に示すような
周波数特性記憶部８−２として実現するのが好適であ
る。周波数特性記憶部８−２は、ＲＯＭやＲＡＭといっ
た記憶装置で構成され、システムコントローラー０の制
御に基づき、受信信号の深さｄに応じて、各深さ毎にあ
らかじめ記憶されていた任意のＭ個の周波数特性列Ｈｊ
（ｗ，ｄ）を発生する。
スペクトラム処理制御部８は、図３（ｂ）に示すような
周波数特性記憶部８−２として実現するのが好適であ
る。周波数特性記憶部８−２は、ＲＯＭやＲＡＭといっ
た記憶装置で構成され、システムコントローラー０の制
御に基づき、受信信号の深さｄに応じて、各深さ毎にあ
らかじめ記憶されていた任意のＭ個の周波数特性列Ｈｊ
（ｗ，ｄ）を発生する。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００３４
【補正方法】変更
【補正内容】
【００３４】ここで、上述と同様の制御を実施すること
により、係数制御スイッチの設定や、動作可能となって
いる超音波プローブの種別に対応して周波数特性記憶部
８−２に記憶されていた任意のＭ個の周波数特性列Ｈｊ
が、Ｈ１（ｗ，ｄ），Ｈ２（ｗ，ｄ），…Ｈｌ（ｗ，
ｄ）の中から選択される。
により、係数制御スイッチの設定や、動作可能となって
いる超音波プローブの種別に対応して周波数特性記憶部
８−２に記憶されていた任意のＭ個の周波数特性列Ｈｊ
が、Ｈ１（ｗ，ｄ），Ｈ２（ｗ，ｄ），…Ｈｌ（ｗ，
ｄ）の中から選択される。
【手続補正４】
【補正対象書類名】図面
【補正対象項目名】図６
【補正方法】変更
【補正内容】
【図６】
【手続補正５】
【補正対象書類名】図面
【補正対象項目名】図１２
【補正方法】変更
【補正内容】
【図１２】
Claims (7)
- 【請求項１】 超音波プローブと、 前記超音波プローブを駆動して被検体中に超音波を送信
する送信手段と、 前記超音波の前記被検体による反射波を前記超音波プロ
ーブで受信して電気信号として取り出す受信手段と、 前記受信手段から出力される受信信号を深さ毎に所定の
スペクトラム波形に変形するスペクトラム変形手段と、 前記スペクトラム変形手段の出力について超音波の画像
を生成するための処理を行う画像生成手段と、 前記画像を表示する手段とを具備したことを特徴とする
超音波診断装置。 - 【請求項２】 前記スペクトラム変形手段は、ディジタ
ルフィルタで構成されることを特徴とする請求項１記載
の超音波診断装置。 - 【請求項３】 前記ディジタルフィルタの周波数特性が
フィルタ入力での平均的な信号スペクトラムとフィルタ
出力での所定の信号スペクトラム期待値とにより決定さ
れることを特徴とする請求項１記載の超音波診断装置。 - 【請求項４】 前記ディジタルフィルタの周波数特性が
多峰性であることを特徴とする請求項２記載の超音波診
断装置。 - 【請求項５】 前記ディジタルフィルタの周波数特性を
前記超音波プローブ毎に切り換える手段をさらに備える
ことを特徴とする請求項２記載の超音波診断装置。 - 【請求項６】 前記前記フィルタ入力での信号スペクト
ラムを逐次演算するスペクトラム演算手段と、前記スペ
クトラム演算手段により得られる信号スペクトラムとフ
ィルタ出力での所定の信号スペクトラム期待値とにより
前記ディジタルフィルタの周波数特性を逐次決定する手
段とをさらに備えることを特徴とする請求項２記載の超
音波診断装置。 - 【請求項７】 前記前記フィルタ出力での所定の信号ス
ペクトラム期待値としてガウス関数を用いることを特徴
とする請求項３又は請求項６記載の超音波診断装置。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP34179895A JP3682104B2 (ja) | 1995-12-27 | 1995-12-27 | 超音波診断装置 |
| US08/777,281 US5911160A (en) | 1995-12-27 | 1996-12-27 | Ultrasound diagnostic apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP34179895A JP3682104B2 (ja) | 1995-12-27 | 1995-12-27 | 超音波診断装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH09173334A true JPH09173334A (ja) | 1997-07-08 |
| JP3682104B2 JP3682104B2 (ja) | 2005-08-10 |
Family
ID=18348846
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP34179895A Expired - Lifetime JP3682104B2 (ja) | 1995-12-27 | 1995-12-27 | 超音波診断装置 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US5911160A (ja) |
| JP (1) | JP3682104B2 (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2013166059A (ja) * | 2010-11-11 | 2013-08-29 | Olympus Medical Systems Corp | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
| WO2014115377A1 (ja) | 2013-01-23 | 2014-07-31 | オリンパスメディカルシステムズ株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6016700A (en) * | 1997-10-29 | 2000-01-25 | Krautkramer Branson, Inc. | Signal processing apparatus for ultrasonic flaw detection system |
| US6142942A (en) * | 1999-03-22 | 2000-11-07 | Agilent Technologies, Inc. | Ultrasound imaging system and method employing an adaptive filter |
| JP2010125118A (ja) * | 2008-11-28 | 2010-06-10 | Toshiba Corp | 超音波診断装置及び超音波診断装置制御方法 |
| GB0910704D0 (en) * | 2009-06-22 | 2009-08-05 | Univ Leeds | A novel tomographic sensing system for high conductivity multiphase flow measurement |
| DE102014201129A1 (de) * | 2014-01-22 | 2015-07-23 | Siemens Aktiengesellschaft | Ultraschallprüfvorrichtung und Verfahren zur Ultraschallprüfung |
Family Cites Families (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4202215A (en) * | 1978-10-26 | 1980-05-13 | Kurt Orban Company, Inc. | Sonic pulse-echo method and apparatus for determining attenuation coefficients |
| JPS5941731A (ja) * | 1982-09-02 | 1984-03-08 | Mitsui Petrochem Ind Ltd | 倉庫内の結露防止方法 |
| JPS59171540A (ja) * | 1983-03-18 | 1984-09-28 | 富士通株式会社 | スペクトル形状利用測定方法 |
| JPH0613027B2 (ja) * | 1985-06-26 | 1994-02-23 | 富士通株式会社 | 超音波媒体特性値測定装置 |
| US5315687A (en) * | 1986-03-07 | 1994-05-24 | Adler Research Associates | Side fed superlattice for the production of linear predictor and filter coefficients |
| US5265217A (en) * | 1987-03-03 | 1993-11-23 | Adler Research Associates | Optimal parametric signal processor for least square finite impulse response filtering |
| US4881549A (en) * | 1987-05-29 | 1989-11-21 | Marquette Electronics | Apparatus and method for obtaining ultrasonic backcatter measurement from tissue |
| US5121364A (en) * | 1991-08-07 | 1992-06-09 | General Electric Company | Time frequency control filter for an ultrasonic imaging system |
| US5349525A (en) * | 1993-01-08 | 1994-09-20 | General Electric Company | Color flow imaging system utilizing a frequency domain wall filter |
| US5349524A (en) * | 1993-01-08 | 1994-09-20 | General Electric Company | Color flow imaging system utilizing a time domain adaptive wall filter |
| US5396285A (en) * | 1993-05-07 | 1995-03-07 | Acuson Corporation | Ultrasound imaging method and apparatus with dynamic non-linear filtering |
| US5609155A (en) * | 1995-04-26 | 1997-03-11 | Acuson Corporation | Energy weighted parameter spatial/temporal filter |
| US5524626A (en) * | 1995-06-30 | 1996-06-11 | Siemens Medical Systems, Inc. | System and method for determining local attenuation for ultrasonic imaging |
-
1995
- 1995-12-27 JP JP34179895A patent/JP3682104B2/ja not_active Expired - Lifetime
-
1996
- 1996-12-27 US US08/777,281 patent/US5911160A/en not_active Expired - Lifetime
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2013166059A (ja) * | 2010-11-11 | 2013-08-29 | Olympus Medical Systems Corp | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
| WO2014115377A1 (ja) | 2013-01-23 | 2014-07-31 | オリンパスメディカルシステムズ株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
| US9360550B2 (en) | 2013-01-23 | 2016-06-07 | Olympus Corporation | Ultrasonic observation apparatus, operation method of the same, and computer readable recording medium |
Also Published As
| Publication number | Publication date |
|---|---|
| US5911160A (en) | 1999-06-08 |
| JP3682104B2 (ja) | 2005-08-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ward et al. | Nonlinear propagation applied to the improvement of resolution in diagnostic medical ultrasound | |
| CN1120695C (zh) | 超声波谱反衬成像 | |
| US7749167B2 (en) | Apparatus and method for ultrasonic diagnostic imaging | |
| US5016641A (en) | Spectral interpolation of ultrasound Doppler signal | |
| JP4796105B2 (ja) | 超音波診断装置 | |
| JP5322522B2 (ja) | 超音波診断装置 | |
| EP0813182B1 (en) | Ultrasonic diagnostic apparatus which performs complex vector phase correction to echo signals obtained in time sequence from each ultrasonic transducer | |
| US20060058677A1 (en) | Ultrasonograph | |
| US7666142B2 (en) | Ultrasound doppler diagnostic apparatus and image data generating method | |
| JP3645347B2 (ja) | 超音波診断装置及び遅延時間最適化方法 | |
| US20030060712A1 (en) | Ultrasonic diagnosis apparatus and ultrasound imaging method | |
| CN1515229B (zh) | 超声波诊断装置 | |
| US20100137715A1 (en) | Ultrasonic imaging apparatus and control method for ultrasonic imaging apparatus | |
| JPH09173334A (ja) | 超音波診断装置 | |
| US8652049B2 (en) | Ultrasonic diagnostic apparatus | |
| JP4728466B2 (ja) | 物質をイメージングする方法及びイメージング・システム | |
| JP2004073620A (ja) | 超音波診断装置 | |
| JP2005245479A (ja) | 超音波診断装置 | |
| JP3887040B2 (ja) | 超音波診断装置 | |
| JP2012005789A (ja) | 超音波診断装置 | |
| JPH0678922A (ja) | 超音波診断装置 | |
| JPH08173431A (ja) | 超音波診断装置 | |
| JP2000139911A (ja) | 音響走査線補間方法および装置と超音波診断装置 | |
| JP4698073B2 (ja) | 超音波診断装置 | |
| JP2002301072A (ja) | 超音波撮像方法及び装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20041207 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050126 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050222 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050420 |
|
| 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: 20050517 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050520 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090527 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090527 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100527 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100527 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110527 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110527 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120527 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120527 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130527 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130527 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140527 Year of fee payment: 9 |
|
| EXPY | Cancellation because of completion of term |