JPH01212543A - 超音波診断装置の信号処理方法 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 〔産業上の利用分野〕 本発明は検体の断層像を得るための超音波診断装置に係
り、詳しくは超音波を検体内に入射し、検体内部の組織
の不連続面からの反射波を受信し、受信した反射波であ
る受信波に基づいて検体の組織の断層像を得るための超
音波診断装置の信号処理方法に関するものである。

〔技術の背景〕

超音波診断装置によって生体等の検体の内部状態を正確
に把握するには、観測によって得られる雑音を含むデー
タから真の情報をいかに取り出すかが重要な問題である
。

音源より放出された超音波が検体中を伝播してから受信
されるまでに種々の変形を受けた波形である入射波をデ
ジタル化して、その時系列なｘ　（ｔ）（ｔ＝０，１．
・・・、ｎ）とし、同様に検査対象である検体内部の反
射係数時系列なＣ（ｔ）　、観測された受信波の時系列
なｙ　（ｔ）とすると ｙ（ｔ）＝Ｃ（ｔ）＊ｘ（ｔ）＝、１９（ｓ）ｘ（ｔ−
ｓ）　・・・・・（１）と表わされる。但し、記号＊は
コンボ８リユージヨン（たたみこみ積分）を表わす。一
般に入射波ｘ（ｔ）は、音源に関する時系列Ｓ　（ｔ）
　、検体内部のフィルタ効果に関する時系列Ｅ　（ｔ）
、受信装置の特性に関する時系列Ｒ（ｔ）によりｘ（ｔ
）　＝　５（ｔ）　＊　Ｅ（ｔ）　＊　Ｒ（ｔ）　　・
・・・・・・・・・・・（２）と表わされる。超音波診
断装置における信号処理とは、検体内部の反射係数時系
列Ｃ（ｔ）を得ることである。そのために、理想的には
入射波ｘ　（ｔ）から波形変形を起す要因をとってイン
パルス（すなわち、瞬間的な単波）に変換する逆フィル
タのオペレータＦ　（ｔ）を求め、オペレータＦ　（ｔ
）と受信波ｙ　（ｔ）とのコンボリューションな行って
Ｃ（ｔ）を得ることになる。すなわち、 Ｃ（ｔ）　＝ｙ（ｔ）　＊　Ｆ（ｔ）′・・・・・・・
・・・・・・・・・・・（３）である。しかし、Ｆ　（
ｔ）を厳密に求めることは不可能であり、近似的にＦ　
（ｔ）を求めることしかできない。その近似方法として
ウィナ−フィルタ演算処理を用いたデコンボリューショ
ン法がある。

ここで、ウィナ−フィルタ演算処理とは、出力波形とイ
ンパルス等の希望波形との誤差の二乗平均が最小になる
ようなオペレータ信号を取り出す演算処理である。

詳しくは、ウィナ−フィルタ演算処理とは、第６図に示
すように、入力波形の自己相関関数ｆ　（ｔ）を求める
演算５０と、入力波形と希望波形との相互相関関数φ（
ｔ”）を求める演算５１、および、ウィナ−の最適フィ
ルタの理論に基づいて、正規方程式 ％式％（４） を解く演算５２によってオペレータ信号Ｆ１（１）を求
める処理からなっている。

また、デコンボリューション法とは、検体内部の反射係
数時系列Ｃ（ｔ）〜がランダム、すなわちＣ（ｔ）の自
己相関関数がインパルスであって、入射波ｘ　（ｔ）が
処理結果と希望波形であるインパルスとの誤差の二乗平
均が最小になる最小位相波形の場合を想定してウィナ−
フィルタ演算処理によってＣ（ｔ）・を得る方法である
。　　。

この場合には、入射波ｘ　（ｔ）の自己相関関数ｆ　（
ｔ）が同時に受信波ｙ　（ｔ）の自己相関関数になり、
更に入射波ｘ　（ｔ）と希望波形であるインパルスとの
□相互相関関数φ（１）が入射波ｘ　（ｔ）と無関係な
定数になる。それ故、入射波ｘ　（ｔ）が未知であって
も、それが最小位相波形であると仮定することによりウ
ィナ−フィルタ演算処理によって（４）式の正規方程式
から変換オペレータＦｌ　（ｔ）を算出することができ
る。

そして、該変換オペレータＦｌ（ｔ）とｙ（１）とのコ
ンボリューションを行うことによって反射係数時系列Ｃ
（ｔ）を求めることができる。

（従来の技術） 従来、超音波診断装置の探触子によって超音波を人体内
に入射し、人体内部の組織の不連続面からの反射波を受
信波として探触子で検出し、この受信波をデコンボリュ
ーション法によって信号処理する方法があった。この方
法は第７図のブロック図に示すように受信波１からその
一部抽出１０１を行ってモデル受信波２とし、このモデ
ル受信波２について、インパルス３を希望波形としてウ
ィナ−フィルタ演算処理１０２を行って、モデル受信波
２を希望波形に変換する変換オペレータ信号４を取り出
し、この変換オペレータ信号４と受信波１とのコンボリ
ューション１０３を行い、反射係数時系列である変換波
５を求めるものである。

（発明が解決しようとする課題） ところが、デコンボリューション法は、技術の背景で述
べたように、入射波が最小位相波形である場合を想定し
た信号処理方法である。入射波が最小位相波形でない場
合には、ウィナ−フィルタ演算処理１０２によって得た
変換オペレータ信号４を使っても受信波１を希望波形に
十分近くする　　−ことができない。したがって、人体
内部の反射係数時系列を十分に再現することができない
ので、超音波診断装置による画像は明瞭なものにならず
画像処理の距離分解能が悪いことになる。

そこで、本発明はデコンボリューション法を適用した場
合に超音波診断装置による画像処理の距離分解能をより
高めることを目的としてなされたものである。

（課題を解決するための手段） 本発明は、超音波診断装置の探触子が検出した受信波を
、インパルスを希望波形とするデコンボリューション法
処理によって変換した変換波に基づいて検体内の検査を
行う超音波診断装置の信号処理方法において、前記受信
波を一旦最小位相波形または遅延インパルスによる変換
波形に変換した後に前記デコンボリューション法処理に
よって変換波を得ることを特徴とする超音波診断装置の
信号処理方法である。

（実施例） 次に本発明の詳細な説明する。

〈第１実施例〉 第１図において、まず、超音波診断装置の発振器によっ
て駆動した探触子によって生体等の検体に超音波を入射
し、検体内部で反射した反射波を探触子が検出して前段
階処理１０４を行って受信波Ａ１を得る。ここで前段階
処理１０４として、探触子が検出した時点での受信波Ａ
１をフィルタにより、ナイキスト周波数以上の成分を遮
断し、アナログ信号である受信波Ａ１をデジタル信号に
変換し、そのデータを記憶する。

また検体の別の部位に超音波を入射した場合の反射受信
波を受信波Ｂ１とし、受信波Ａ１と同様の前段階処理１
０４を行う。

次に該受信波Ａ１に関するデータから一部抽出１０５を
行って受信波Ａ１を代表させたモデル受信波２に相当す
るデータを読み出す。

続いてモデル受信波２に最小位相モデル受信波演算１０
６を行って最小位相波形である最小位相モデル受信波９
を得る。最小位相モデル受信波演算１０６として、この
例では二重Ｗｉｅｎｅｒ−Ｌｅｖｉｎｓｏｎ法による演
算を行っている。すなわち、該演算１０６はモデル受信
波２にインパルス６を希望波形としてウィナ−フィルタ
演算処理１０７を行ってウィナ−フィルタ信号７を取り
出す過程と、該ウィナ−フィルタ信号７にインパルス８
を希望波形としてウィナ−フィルタ演算処理１０８を行
って最小位相モデル受信波９を得る過程とからなってい
る。

更に、該最小位相モデル受信波９、前記モデル受信波２
及び受信波Ｂ１に対して最小位相受信波演算１０９を行
い、最小位相受信波１２を得る。

この最小位相受信波演算１０９はモデル受信波２に対し
て最小位相モデル受信波９を希望波形としてウィナ−フ
ィルタ演算処理１１０を行って最小位相化オペレータ信
号１０を取り出すと′ともに、該受信波Ｂ１と最小位相
化オペレータ信号１０とのコンボリューション１１１を
行い最小位相受信波１２を得る演算からなっている。

その後、最小位相受信波１２に対してデコンボリューシ
ョン法処理１１２を行って変換波１３を得る。ここでデ
コンボリューション法処理１１２とは前記最小位相受信
波１２に対してインパルス３を希望波形としたウィナ−
フィルタ演算処理１０２を行い変換オペレータ信号１１
を取り出すとともに、該変換オペレータ信号１１と最小
位相受信波１２とのコンボリューション１１３を行い変
換波１３を得る処理である。

最後に、該変換波１３の後段階処理１１４を行う。後段
階処理１１４として、デジタル化されている変換波１３
を再度アナログ化し、対数増幅器によりて小振幅信号は
ど大きく増幅し、検波器によって両波整流し、フィルタ
により平滑化し、画像表示回路を経てＣＲＴで画像とす
る。以上の処理で受信波Ｂ１は受信波Ａ１であっても差
支ない。こうして、本実施例では受信波Ｂ１およびＡ１
の各々に対して最小位相化オペレータ信号１０を取り出
すのではなく、モデル受信波２から最小位相化オペレー
タ信号１０を取り出し、該最小位相化オペレータ信号１
０と受信波Ｂ１またはＡ１とのコンボリューション１１
１を行って最小位相受信波１２を得ているので、受信波
Ｂ１およびＡ１の各々から最小位相化オペレータ信号１
０を取り出す場合に比べて処理手数が少なくてすむ。な
お、モデル受信波２としては必ずしも受信波ＡＩの一部
を抽出して代表させたものではなく、受信波Ａ１を理想
化した波形で代表させたものでもよい。

〈第２実施例） 第２図は、第１図が受信波Ｂ１を一旦最小位相波形に変
換した後デコンボリューション法処理を行って変換波を
得るに対し、受信波Ｂ１を−Ｈ遅延インパルスによる変
換波形に変換した後デコンボリューション法処理を行う
ものである。

すなわち、本実施例では、検体内部で反射した反射波を
探触子が検出して前段階−処理１０４を行って受信波Ａ
１を得、次に該受信波Ａ１から一部抽出１０５を行って
受信波Ａ１を代表するモデル受信波２を得ることは第１
図と同様である。

本実施例においては、このモデル受信波２に対して、遅
延インパルス変換演算１３６を行い遅延インパルス変換
受信波２３を得る。遅延インパルス変換演算１３６゛に
おいては、モデル受信波２にウィナ−フィルタ演算処理
１２７を施すのであるが、この場合の希望波形を遅延イ
ンパルス２１とする。モデル受信波２が最小位相波形で
ない場合、希望波形であるインパルスをモデル受信波２
の始動時刻に定めても、インパルスに十分近い出力は得
られない。そこで、希望波形であるインパルスにある遅
延時間を与えてウィナ−フィルタ演算１２７を行うと、
インパルスに近い出力が得られ、遅延時間を変えて演算
をくり返せば最適遅延時間が得られる、この場合、その
モデル受信波２について最適の遅延インパルス変換オペ
レータ信号２２が得られる。

このオペレータ信号２２と受信波Ｂ１とのコンボリュー
ション１３１に・よって遅延インパルス変換受信波２３
が得られる。こうして得られた波形はインパルスに近い
波形であるが、これに更にデコンボリューション法処理
１１２を施すことにより所望の変換波１３を得る。

く第３実施例） 次に第３の実施例を第３図及び第４図に基づいて説明す
る。

第３図は第１図に対してウィナ−フィルタ演算処理の入
力が異っている。即ち第１図は最小位相受信波１２が入
力であるのに対し、第３図では最小位相モデル受信波９
が入力となる。その結果第３図の変換オペレータ信号１
１は第１図の変換オペレータ１１を近似したものとなる
。

第４図は第３図を簡略化したものである。第３の実施例
においても、第１の実施例と同様に、まず検体内部で反
射した反射波を探触子が検出して前段階処理１０４を行
って受信波Ａ１を得る。次に、該受信波Ａ１から一部抽
１１０５を行って受信波Ａ１を代表するモデル受信波２
を得る。、続いて、第１の実施例と同様にモデル受信波
２に最小位相モデル受信波演算１０６を行って最小位相
波形である最小位相モデル受信波９を得る。さらに、前
記モデル受信波２に最小位相モデル受信波９を希望波形
としてウィナ−フィルタ演算処理１１０を行って、最小
位相化オペレータ信号１０を取り出すための最小位相化
オペレータ演算１１５を行う。一方、前記最小位相モデ
ル受信波９、最小位相化オペレータ信号１０及び受信波
１に対してデコンボリューション法処理１１６を行い変
換波１３を得る。ここで、デコンボリューション法処理
１１６とは最小位相モデル受信波９に対してインパルス
３を希望波形とするウィナ−フィルタ演算処理１０２を
行い変換オペレータ信号１１を取り出し、受信波１、最
小位相化オペレータ信号１０及び変換オペレータ信号１
１のコンボリューション１１７を行って変換波１３を得
ることである。最後に該変換波１３に後段階処理１１４
を行う。

第４図の本実施例では、以上述べたように、第１の実施
例と異なり最小位相受信波をあられに得ていない。これ
は、第１の実施例とはコンボリューションの順序と結合
とが異なっているからである。しかし、コンボリューシ
ョンでは、交換剤と結合剤とが成立するので、本実施例
でも表面には表われないが、第１の実施例と同様に、や
はり最小位相受信波を得ていることになる。

〈第４実施例〉 次に第４の実施例を第５図に基づいて説明する。

第４の実施例においても、第１の実施例と同様に、まず
検体内部で反射した反射波を探触子が検出して前段階処
理１０４を行って受信波Ａ１を得る。次に、該受信波Ａ
１から一部抽出１０５を行って受信波Ａ１を代表させた
モデル受信波２を得る。続いて、モデル受信波２に最小
位相モデル受信波演算１１８を行って、最小位相モデル
受信波９及びウィナ−フィルタ信号７を取り出す。この
最小位相モデル受信波演算１１８は、第１及び第２の実
施例の最小位相モデル受信波演算１０６と同じ演算を行
うが、その演算中に生じるウィナ−フィルタ信号７を取
り出す点が異なる。次に、該最小位相モデル受信波９、
モデル受信波２及び受信波Ｂ１より第１の実施例と同様
に最小位相受信波演算１０９を行い最小位相受信波１２
を得る。最後に、前記ウィナ−フィルタ信号７と該最小
位相受信波１２とのコンボリューション１１３を行って
変換波１３を得るためのデコンボリューション法処理１
１９を行う。このデコンボリューション法処理１１９は
変換波１３を得る点では第１及び第２の実施例のデコン
ボリューション処理１１２，１１６と同じであるが、変
換波１３を取り出すためのウィナ−フィルタ演算処理を
行っていない点で異なる。したがって、その分、処理の
手数が減ることになる。これはウィナ−フィルタ信号７
が第１及び第２の実施例における変換オペレータ信号１
１と極めて類似していることを利用したものである。

く第５実施例〉 第５の実施例は、前段階処理１０４または後段階処理１
１４等において、受信波１を該装置固有の特性、検体固
着の特性及び外部環境のうちの少なくとも１つに応じて
補正する操作を施したものである。受信波１を補正する
操作としては次のものがある。

その第１番目は、最良のモデル受信波２を選択する操作
に関するものである。この操作を施す理由は検体内の吸
収によって入射波９反射波のスペクトルは距離に依存す
るからである。そこで、最良の結果を得るためには、検
体の深さや特性に応じた複数個のモデル受信波２または
単一のモデル受信波から複数個のモデルを与える演算方
法を用′意しておき、前段階処理１０４において、検体
内の深さや検体の特性に応じてモデル受信波２を切り換
えるようにする。

こうして、受信波１に対して複数個のモデル受信波２を
用いてデコンボリューション法処理等のウェブレット処
理を行うことができる。

例えば、第８図はヒマシ油中の種々の距離（１５ｍｍ、
２５．５ｍｍ、３６．５ｍｍ、４５．５ｍｍ、５５．０
ｍｍ、６５．０ｍｍ。

７５．０ｍｍ）に反射物としてアクリル板を設置し、超
音波診断装置の探触子から超音波を発射してその反射受
信波を記録したものである。図中りは探触子アクリル板
との距離を表わし、ａ　−ｇはそれぞれの距離に応じた
反射受信波を表わす。各反射受信波の振巾は、最大値が
等しくなるように正規化しである。実際はヒマシ油の距
離依存減衰があるため、ａからｇへ移るに従って振巾が
小さくなつている。ｇのノイズが大きく見えるのはその
ためである。

減衰のみならず、ｇの波形はａに比べて主たる周波数成
分が変化している率がわかる。ゼロクロス周期がａの約
０．３ｐｓであるのに対し、ｇでは約０．４ｐ、ｓであ
る。これをｂをモデル受信波として第１図の処理を葬っ
た結果が第９図である。ｂに近い受信波ではスパイク波
が得られているが、ｇではほとんどスパイク波とは言い
難い。これを以下に述べる方法で処理した結果が第１０
図で、第９図に比べて著しい改善が見られる。

第１１図は第８図のＱ　％　ｇの振巾−周波数特性であ
る。レベルはａの最大振巾成分をＯｄＢにとっである。

これからａの各周波数成分をそれぞれＯｄＢとして５〜
ｇの各周波数成分のａに対する減衰を求めると第１２図
となる。８〜ｇの主たる周　−波数成分的１．１〜４　
ＭＨｚで考えると減衰が距離及び周波数に依存している
ことが判る。従９てこれから各周波数毎に距離Ｉｃｍ当
りの減衰量を求めると第１３図が得られ、ばらつきはあ
るものの１本の推定曲線が求められる。これからａは距
離ゼロにおける仮想波形に対して各周波数成分がどれ程
減衰しているかが求まり、それを第１２図に適用して仮
想波形の各周波数成分なＯｄＢにとったのが、第１４図
である。これから各曲線の傾斜を求め、これらを１ｃｍ
当りの減衰量を求めると、第１５図の吸収係数が得られ
る、あるいは第１４図の各周波数成分の減衰量を１ｃｍ
当りの値に直せば第１６図が得られ、これらから１本の
推定曲線が吸収係数として求められる。吸収係数の計算
に当り、距離は複数距離が用いである。

第１７図は第８図Ｊｉｌ　％　ｇの位相周波数特性であ
る。１．５〜４　ＭＨｚの範囲では何れも略々平行なの
で、この場合は距離毎に位相特性を変える必要は。ない
が、場合によっては１ｃｍ当りの位相周波数特性を求め
る必要があろう。以上の減衰１位相の周波数特性から距
離ゼロの合成仮想波形及び距離毎の合成波形を求めると
第１８図のようになる。第１０図はこうして得られた第
１８図ｆ１ｇをモデル受信波として第１図に示した第１
の実施例の処理な行つた結果である。　　　　　　　　
　　　　　。

その第２番目は、受信波ｌの減衰を補正するゲインコン
トロールを施したものである。ゲインコントロールを施
す理由は、超音波の反射波が検体内で減衰を受けるため
受信波１の振幅が生体内通過距離によって小さくなるか
らである。すなわち、画像の明るさは振幅の大きさによ
って決まるため、この減衰を補正しないと、探触子に近
い検体部分に相当する画像は明るくなり、探触子から遠
い検体部分に相当する画像は暗くなるために本実施例が
必要となる。なお、該ゲインコントロールは、精度を下
げないようにアナログ信号をデジタル信号に変換する前
に行うのが適当である。

その第３番目は、複数個の探触子から得た受信波ｌにつ
いて、それぞれ本発明に係る信号処理方法によって信号
処理を行って得た変換波１３を探触子の位置に関する時
間の遅れを遅延回路によって補正し、同一時刻に揃えた
後にＣＤ　Ｐ　（ＣｏｍｍｏｎＤｅｐｔｈ　Ｐｏ１ｎｔ
）重合演算器で重合することによって変換波１３の振幅
は大となる。一方、変換波１３の前後のノイズは検体内
の異なる経路によって生じているため、位相がまちまち
であり、重合後のノイズの振幅はそれ程大きくならない
。

こうしてＳ−Ｎ比が改善されることになる。

また補正遅延時間によって重合信号反射点までの平均速
度が算出できる。

〈第６実施例） さらに第６の実施例を説明する。本発明に係る信号処理
方法で得られた変換波１３を正負検出器に入力して、正
負の判別を行い、カラー表示回路によって正スパイク画
像と負スパイク画像とを異なる色て表示する。これは、
検体肉組織で超音波が反射する一場合、低密度組織から
高密度組織に進む場合に反射すると、反射波形は入射波
形と同位相であり、逆の場合は、逆位相となることから
である。

こうして本実施例の表示によって組織の性質が表現され
ることになる。

（発明の効果） 本発明は、超音波診断装置が検出した受信波を、−旦最
小位相受信波形に変換した後に、デコンボリューション
法処理を適用して変換波を求めているため、従来に比べ
てより理想的な希望波形に変換することができる。こう
して検体内の組織についての情報をより正確に把握する
ことができるので、その情報によって得られる画像の分
解能が高められることになる。

【図面の簡単な説明】

第１図は本発明の第１の実施例を示すブロック図、第２
図は本発明の第２の実施例を示すブロック図、第３図及
び第４図は本発明の第３の実施例を示すブロック図、第
５図は本発明の第４の実施例を示すブロック図、第６図
はウィナ−フィルタ演算処理を示すブロック図、第７図
は従来の信号処理方法を示すブロック図、第８図はヒマ
シ油中の種々の距離に設けたアクリル板からの反射受信
　−波の正規化波形図、第９図は第８図の波形６をモデ
ル受信波として第１の実施例に係る処理によって変換し
た変換波形図、第１０図は第８図の波形に距離依存の補
正を施した変換波形図、第１１図は第８図の波形の振幅
・周波数特性を示す図、第１２図は第８図のａに対する
ｂ−ｇの減衰の周波数依存性を示す図、第１３図は第１
２図の単位距離あたりの減衰量、第１４図は距離０にお
ける合成仮想波形に対するｆｉ　％　ｇの波形の減衰の
周波数依存性を示す図、第１５図は吸収係数を示す図、
第１６図は第１４図の各周波数成分の単位距離あたりの
減衰量を示す図、第１７図は第８図の８〜ｇの波形の位
相の周波数特性を示す図、第１８図は以上の減衰、位相
の周波数特性から求めた距離の合成仮想波形及び各距離
の合成波形を示す波形図である。 １・・・受信波 ２・・・モデル受信波 ３．６．８・・・インパルス ９・・・最小位相モデル受信波 １０・・・最小位相化オペレータ信号 １１・・・変換オペレータ信号 １２・・・最小位相受信波 １３・・・変換波 ｊＩ４図 ｊｌａ！Ｉ！！！ 第８図 ｑ　　　　　　Ｏ 第１１図 ＦＲＥＱＵＥＮＣＹ（ＭＨ２） 第１２図 ＲＥＦＥＲＦ？ＥＤ　Ｔｏ　　Ｄ−１５，０ＭＭ、　　
ＲＥＦＬＥＣＴＩＯＮＦＲＥＱＵＥＮＣＹ（ＭＨｚ） 第１３図 ＦＲＥＱＵＥＮＣＹ（ＭＨ２） ３１１７図　　　　　　　　８ Ｕ」 の く 工 〇− ＦＲＥＱＵＥＮＣＹ（ＭＨ２１

Claims (4)

【特許請求の範囲】
1. （１）超音波診断装置の探触子が検出した受信波を、イ
   ンパルスを希望波形とするデコンボリューション法処理
   によって変換した変換波に基づいて検体内の検査を行う
   超音波診断装置の信号処理方法において、前記受信波を
   一旦最小位相波形または遅延インパルスによる変換波形
   に変換した後に前記デコンボリューション法処理によっ
   て変換波を得ることを特徴とする超音波診断装置の信号
   処理方法。
2. （２）超音波診断装置の探触子が検出した受信波を代表
   させたモデル受信波を最小位相波形である最小位相モデ
   ル受信波に変換するとともに、ウィナーフィルタ演処理
   によって前記モデル受信波を最小位相モデル受信波に変
   換する最小位相化オペレータ信号を取り出し、該最小位
   相化オペレータ信号と受信波とのコンボリューションを
   行って受信波を最小位相波形である最小位相受信波に一
   旦変換すると共に、前記最小位相モデル受信波に対して
   インパルスを希望波形とするデコンボリューション法処
   理を行って取り出した変換オペレータ信号と前記最小位
   相受信波とのコンボリューションを行って変換波を得る
   請求項１記載の超音波診断装置の信号処理方法。
3. （３）超音波診断装置の探触子が検出した受信波を代表
   させたモデル受信波に、遅延インパルスを希望波形とす
   るウィナーフィルタ演算処理によって遅延インパルス変
   換オペレータ信号を取り出し、これと受信波とのコンボ
   リューションから遅延インパルスに対する最小位相の波
   形である変換受信波に一旦変換した後、この遅延インパ
   ルスによる変換受信波に対してインパルスを希望波形と
   するデコンボリューション法処理によって得た変換オペ
   レータ信号を前記遅延インパルスによる変換受信波にコ
   ンボリューションして変換波を得る請求項１記載の超音
   波診断装置の信号処理方法。
4. （４）超音波診断装置の探触子が検出した受信波を、該
   装置固有の特性、検体固有の特性及び外部環境のうちの
   少なくとも１つに応じて受信波を変換する補正処理を行
   うことを特徴とする請求項１記載の超音波診断装置の信
   号処理方法。