(57)【要約】      （修正有）【課題】ノイズの影響を低減し、正確な性状測定を行うことのできる超音波診断装置を提供する。【解決手段】生体組織へ超音波を送信するため超音波プローブ１３を駆動する送信部１４と、生体組織からの超音波反射波を、前記超音波プローブ１３を用いて受信する受信部１５と、受信信号を位相検波し、前記生体組織に設定された複数の測定対象位置から位相検波信号を生成する位相検波部１７と、前記複数の測定対象位置の複数領域を前記生体組織に設定し、各領域内の位相検波信号の空間平均処理を行って平均位相検波信号を求め、前記各領域内の測定対象位置における位相検波信号の大きさに基づく前記平均位相検波信号の重心位置を求め、前記平均位相検波信号および前記重心位置に基づいて、前記複数の測定対象位置の位置変位量、前記複数の測定対象位置の２つの間の厚さ変化量、歪み量、弾性特性のうち少なくとも一つを求める演算部１９とを備える超音波診断装置。【選択図】図２
2007020724.tif 000002
【特許請求の範囲】
【請求項１】
生体組織へ超音波送信波を送信するための超音波プローブを駆動する送信部と、
前記超音波送信波が前記生体組織において反射することにより得られる超音波反射波を、前記超音波プローブを用いて受信し、受信信号を生成する受信部と、
前記受信信号を位相検波し、前記生体組織に設定された複数の測定対象位置にそれぞれ対応した複数の位相検波信号を生成する位相検波部と、
前記複数の測定対象位置の少なくとも２つをそれぞれ含む複数の領域を前記生体組織に設定し、各領域内の位相検波信号の空間平均処理を行うことにより平均位相検波信号を求め、前記各領域内の測定対象位置における位相検波信号の大きさに基づいて前記平均位相検波信号の重心位置を求め、前記平均位相検波信号および前記重心位置に基づいて、前記複数の測定対象位置の位置変位量、前記複数の測定対象位置から任意に選ばれる２つの間の厚さ変化量、歪み量、弾性特性のうち少なくとも一つを求める演算部と、
を備える超音波診断装置。
【請求項２】
生体組織へ超音波送信波を送信するための超音波プローブを駆動する送信部と、
前記超音波送信波が前記生体組織において反射することにより得られる超音波反射波を、前記超音波プローブを用いて受信し、受信信号を生成する受信部と、
前記受信信号を位相検波し、前記生体組織に設定された複数の測定対象位置にそれぞれ対応した複数の位相検波信号を生成する位相検波部と、
前記複数の位相検波信号の時間相関演算を行うことにより、複数の相関信号を求め、前記複数の測定対象位置の少なくとも２つをそれぞれ含む複数の領域を前記生体組織に設定し、各領域内の相関信号の空間平均処理を行うことにより平均相関信号を求め、前記各領域内の測定対象位置における相関信号の大きさに基づいて前記平均相関信号の重心位置を求め、前記平均相関信号および前記重心位置に基づいて、前記複数の測定対象位置の位置変位量、前記複数の測定対象位置から任意に選ばれる２つの間の厚さ変化量、歪み量、弾性特性のうち少なくとも一つを求める演算部と、
を備える超音波診断装置。
【請求項３】
前記演算部は、所定の重み付け係数を用い、前記各領域内の平均位相検波信号および重心位置の少なくとも一方を求める請求項１または２に記載の超音波診断装置。
【請求項４】
前記演算部は、前記複数の領域から選ばれる任意の２つ領域の重心位置の差分から求められる厚さを用い、前記２つの領域間の歪み量および弾性特性の少なくとも１つを求める請求項１から３のいずれかに記載の超音波診断装置。
【請求項５】
前記各領域内の複数の測定対象位置は同一音響線上に位置している請求項１から４のいずれかに記載の超音波診断装置。
【請求項６】
超音波診断装置の制御部による超音波診断装置の制御方法であって、
超音波プローブを駆動し、超音波を送信するステップと、
前記超音波が生体組織において反射することにより得られる超音波反射波を前記超音波プローブにより受信し、受信信号を生成するステップと、
前記受信信号を位相検波し、前記生体組織の複数の測定対象位置にそれぞれ対応した複数の位相検波信号を生成するステップと、
前記複数の測定対象位置の少なくとも２つをそれぞれ含む複数の領域を前記生体組織に設定し、各領域内の位相検波信号の空間平均処理を行うことにより平均位相検波信号を求めるステップと、
前記領域内の測定対象位置における位相検波信号の大きさに基づいて前記平均位相検波信号の重心位置を求めるステップと、
前記平均位相検波信号および前記重心位置に基づいて、前記複数の測定対象位置の位置変位量、前記複数の測定対象位置から任意に選ばれる２つの間の厚さ変化量、歪み量、弾性特性のうち少なくとも一つを求めるステップと、
を包含する超音波診断装置の制御方法。
【請求項７】
超音波診断装置の制御部による超音波診断装置の制御方法であって、
超音波プローブを駆動し、超音波を送信するステップと、
前記超音波が生体組織において反射することにより得られる超音波反射波を前記超音波プローブにより受信し、受信信号を生成するステップと、
前記受信信号を位相検波し、前記生体組織の複数の測定対象位置にそれぞれ対応した複数の位相検波信号を生成するステップと、
前記複数の位相検波信号の時間相関演算を行うことにより、複数の相関信号を求めるステップと、
前記複数の測定対象位置の少なくとも２つをそれぞれ含む複数の領域を前記生体組織に設定し、各領域内の相関信号の空間平均処理を行うことにより平均相関信号を求めるステップと、
前記領域内の測定対象位置における相関信号の大きさに基づいて前記平均相関信号の重心位置を求めるステップと、
前記平均相関信号および前記重心位置に基づいて、前記複数の測定対象位置の位置変位量、前記複数の測定対象位置から任意に選ばれる２つの間の厚さ変化量、歪み量、弾性特性のうち少なくとも一つを求めるステップと、
を包含する超音波診断装置の制御方法。
【請求項８】
前記重心位置を求めるステップは、所定の重み付け係数を用い、前記各領域内の平均位相検波信号または平均相関信号と重心位置との少なくとも一方を求める請求項６または７に記載の超音波診断装置の制御方法。
【請求項９】
前記位置変位量、厚さ変化量、歪み量、弾性特性のうち少なくとも一つを求めるステップは、前記複数の領域から選ばれる任意の２つ領域の重心位置の差分から求められる厚さを用い、前記２つの領域間の歪み量および弾性特性の少なくとも１つを求める請求項６から８のいずれかに記載の超音波診断装置の制御方法。
【請求項１０】
前記各領域内の複数の測定対象位置は同一音響線上に位置している請求項６から９のいずれかに記載の超音波診断装置の制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波診断装置に関し、特に、生体内の組織の性状診断を行う超音波診断装置および超音波診断装置の制御方法に関する。
【背景技術】
【０００２】
近年、心筋梗塞や脳梗塞などの循環器系疾病を患う人々が増加してきており、このような疾病の予防および治療を行うことが大きな課題となっている。
【０００３】
心筋梗塞や脳梗塞の発病には、動脈硬化が深く関係している。具体的には、動脈壁に粥腫が形成されたり、高血圧等の種々の要因によって動脈の新しい細胞が作られなくなったりすると、動脈は弾力性を失い、硬く、脆くなる。そして、粥腫が形成された部分において血管が閉塞したり、粥腫を覆う血管組織が破裂することにより粥腫が血管内へ流出し、別の部分において動脈を閉塞させたり、動脈が硬化した部分が破裂したりすることによって、これらの疾病が引き起こされる。このため、動脈硬化を早期に診断することがこれらの疾病予防や治療には重要となる。
【０００４】
従来、動脈硬化病変の診断は、血管カテーテルを用いて血管内部の様子を直接観察することによって行われていた。しかし、この診断には、血管カテーテルを血管に挿入する必要があるため、被験者への負荷が大きいという問題があった。このため、血管カテーテルによる観察は、動脈硬化病変が存在していることが確かである被験者に対して、その場所を特定するために用いられ、例えば、健康管理のための検査として、この方法が用いられることはなかった。
【０００５】
動脈硬化の一因であるコレステロール値を測定したり、血圧値を測定したりすることは、被験者への負担が少なく、容易に行うことのできる検査である。しかし、これらの値は、動脈硬化の度合いを直接示すものではない。
【０００６】
また、動脈硬化を早期に診断して、動脈硬化の治療薬を被験者に対して投与することができれば、動脈硬化の治療に効果を発揮する。しかし、動脈硬化が進行してしまうと、治療薬によって動脈硬化の進展を抑制することはできても、硬化した動脈を完全に回復させることは難しいと言われている。
【０００７】
こうした理由から、被験者への負担が少なく、動脈硬化が進行する前に早期段階で診断する診断方法あるいは診断装置が求められている。
【０００８】
一方、被験者への負担が少ない非侵襲の医療診断装置として、超音波診断装置やＸ線診断装置が従来用いられている。超音波やＸ線を体外から照射することによって、被験者に苦痛を与えることなく、体内の形状情報、あるいは形状の時間変化情報を得ることができる。体内の測定対象物の形状の時間変化情報（運動情報）が得られると、測定対象物の性状情報を求めることができる。つまり、生体内の血管の弾性特性を求めることができ、動脈硬化の度合いを直接知ることが可能となる。
【０００９】
特に超音波診断は、Ｘ線診断と比較した場合、被験者に超音波プローブをあてるだけで測定できるので、被験者への造影剤投与が不要である点やＸ線被爆のおそれがない点で優れている。
【００１０】
また、近年のエレクトロニクス技術の進歩によって、超音波診断装置の測定精度を飛躍的に向上させることも可能になってきた。これに伴って、生体組織の微小運動を計測する超音波診断装置の開発が進んでいる。例えば、特許文献１は、制約付き最小二乗法を用いて超音波エコー信号の振幅と位相を解析することにより、測定対象を高精度でトラッキングする技術を開示している。この技術を位相差トラッキング法と呼ぶ。この技術によれば、血管運動の振幅数ミクロンで数百Ｈｚまでの速い振動成分を高精度に計測できるため、血管壁の厚さ変化や歪みを数ミクロンのオーダーで高精度な計測をすることが可能になると報告されている。
【００１１】
このような高精度な計測手法を用いることにより、動脈壁の弾性特性の二次元分布を詳細に測定することが可能となる。例えば非特許文献１では、頸動脈血管壁の弾性率の二次元分布の様子をＢモード断層像に重ねて表示した一例を示している。
【００１２】
動脈壁の硬さ度合いは一様ではなく、ある分布を持って存在していることが知られおり、動脈硬化症の診断においては、動脈壁の硬さの分布を正確に把握することが重要である。非特許文献１に開示された方法によれば、動脈の硬化度合いを示す特徴量である弾性率が二次元で表示されるため、動脈壁の硬化した部分をより正確に特定することが可能となる。
【特許文献１】特開平１０−５２２６号公報
【非特許文献１】Hiroshi Kanai et al, “Elasticity Imaging of Atheroma With Transcutaneous Ultrasound Preliminary Study”, Circulation, Vol.107, p.3018-3021, 2003.
【発明の開示】
【発明が解決しようとする課題】
【００１３】
しかしながら、特許文献１および非特許文献１などに開示された技術によって、数ミクロンオーダーの計測が可能となるにつれてノイズの影響が問題となる。具体的には、超音波エコー信号には測定環境に応じて様々なノイズを含んでいるため、得られる位置変位量が不安定となる。その結果、性状測定値を正確に求めることが困難となる。
【００１４】
本発明はこのような課題を解決し、ノイズの影響を低減し、正確な性状測定を行うことのできる超音波診断装置を提供することを目的とする。
【課題を解決するための手段】
【００１５】
本発明の超音波診断装置は、生体組織へ超音波送信波を送信するための超音波プローブを駆動する送信部と、前記超音波送信波が前記生体組織において反射することにより得られる超音波反射波を、前記超音波プローブを用いて受信し、受信信号を生成する受信部と、前記受信信号を位相検波し、前記生体組織に設定された複数の測定対象位置にそれぞれ対応した複数の位相検波信号を生成する位相検波部と、前記複数の測定対象位置の少なくとも２つをそれぞれ含む複数の領域を前記生体組織に設定し、各領域内の位相検波信号の空間平均処理を行うことにより平均位相検波信号を求め、前記各領域内の測定対象位置における位相検波信号の大きさに基づいて前記平均位相検波信号の重心位置を求め、前記平均位相検波信号および前記重心位置に基づいて、前記複数の測定対象位置の位置変位量、前記複数の測定対象位置から任意に選ばれる２つの間の厚さ変化量、歪み量、弾性特性のうち少なくとも一つを求める演算部とを備える。
【００１６】
また、本発明の他の超音波診断装置は、生体組織へ超音波送信波を送信するための超音波プローブを駆動する送信部と、前記超音波送信波が前記生体組織において反射することにより得られる超音波反射波を、前記超音波プローブを用いて受信し、受信信号を生成する受信部と、前記受信信号を位相検波し、前記生体組織に設定された複数の測定対象位置にそれぞれ対応した複数の位相検波信号を生成する位相検波部と、前記複数の位相検波信号の時間相関演算を行うことにより、複数の相関信号を求め、前記複数の測定対象位置の少なくとも２つをそれぞれ含む複数の領域を前記生体組織に設定し、各領域内の相関信号の空間平均処理を行うことにより平均相関信号を求め、前記各領域内の測定対象位置における相関信号の大きさに基づいて前記平均相関信号の重心位置を求め、前記平均相関信号および前記重心位置に基づいて、前記複数の測定対象位置の位置変位量、前記複数の測定対象位置から任意に選ばれる２つの間の厚さ変化量、歪み量、弾性特性のうち少なくとも一つを求める演算部とを備える。
【００１７】
ある好ましい実施形態において、前記演算部は、所定の重み付け係数を用い、前記各領域内の平均位相検波信号および重心位置の少なくとも一方を求める。
【００１８】
ある好ましい実施形態において、前記演算部は、前記複数の領域から選ばれる任意の２つ領域の重心位置の差分から求められる厚さを用い、前記２つの領域間の歪み量および弾性特性の少なくとも１つを求める。
【００１９】
ある好ましい実施形態において、前記各領域内の複数の測定対象位置は同一音響線上に位置している。
【００２０】
本発明の超音波診断装置の制御方法は、超音波診断装置の制御部による超音波診断装置の制御方法であって、超音波プローブを駆動し、超音波を送信するステップと、前記超音波が生体組織において反射することにより得られる超音波反射波を前記超音波プローブにより受信し、受信信号を生成するステップと、前記受信信号を位相検波し、前記生体組織の複数の測定対象位置にそれぞれ対応した複数の位相検波信号を生成するステップと、前記複数の測定対象位置の少なくとも２つをそれぞれ含む複数の領域を前記生体組織に設定し、各領域内の位相検波信号の空間平均処理を行うことにより平均位相検波信号を求めるステップと、前記領域内の測定対象位置における位相検波信号の大きさに基づいて前記平均位相検波信号の重心位置を求めるステップと、前記平均位相検波信号および前記重心位置に基づいて、前記複数の測定対象位置の位置変位量、前記複数の測定対象位置から任意に選ばれる２つの間の厚さ変化量、歪み量、弾性特性のうち少なくとも一つを求めるステップとを包含する。
【００２１】
また、本発明の他の超音波診断装置の制御方法は、超音波診断装置の制御部による超音波診断装置の制御方法であって、超音波診断装置の制御部による超音波診断装置の制御方法であって、超音波プローブを駆動し、超音波を送信するステップと、前記超音波が生体組織において反射することにより得られる超音波反射波を前記超音波プローブにより受信し、受信信号を生成するステップと、前記受信信号を位相検波し、前記生体組織の複数の測定対象位置にそれぞれ対応した複数の位相検波信号を生成するステップと、前記複数の位相検波信号の時間相関演算を行うことにより、複数の相関信号を求めるステップと、前記複数の測定対象位置の少なくとも２つをそれぞれ含む複数の領域を前記生体組織に設定し、各領域内の相関信号の空間平均処理を行うことにより平均相関信号を求めるステップと、前記領域内の測定対象位置における相関信号の大きさに基づいて前記平均相関信号の重心位置を求めるステップと、前記平均相関信号および前記重心位置に基づいて、前記複数の測定対象位置の位置変位量、前記複数の測定対象位置から任意に選ばれる２つの間の厚さ変化量、歪み量、弾性特性のうち少なくとも一つを求めるステップとを包含する。
【００２２】
ある好ましい実施形態において、前記重心位置を求めるステップは、所定の重み付け係数を用い、前記各領域内の平均位相検波信号または平均相関信号と重心位置との少なくとも一方を求める。
【００２３】
ある好ましい実施形態において、前記位置変位量、厚さ変化量、歪み量、弾性特性のうち少なくとも一つを求めるステップは、前記複数の領域から選ばれる任意の２つ領域の重心位置の差分から求められる厚さを用い、前記２つの領域間の歪み量および弾性特性の少なくとも１つを求める。
【００２４】
ある好ましい実施形態において、前記各領域内の複数の測定対象位置は同一音響線上に位置している。
【発明の効果】
【００２５】
本発明の超音波診断装置によれば、各測定対象位置から得られる位相検波信号または相関信号を空間平均処理することによってノイズの影響を低減するとともに、各測定対象位置から得られる位相検波信号または相関信号の重心位置を用いて、測定対象とする組織の位置を正確に追跡し、また、弾性率を求めるための任意の２つの測定位置間の厚さを正確に求めることが可能となる。したがって、被検体の性状特性、特に組織の歪み量や弾性率を正確に求めることができる。
【発明を実施するための最良の形態】
【００２６】
本発明の超音波診断装置は、測定時におけるノイズの影響を低減するために、被検体の各測定対象位置から得られる情報を空間平均処理する。具体的には、特許文献１および非特許文献１などに開示された位相差トラッキング法により、各測定対象位置から得られる位相検波信号、または、位相検波信号を時間軸上で相関演算処理することにより得られる相関信号の空間平均処理を行う。これにより、位相検波信号あるいは相関信号に重畳しているノイズの影響が低減され、各測定対象位置の位置変位量が安定する。
【００２７】
このとき、空間平均処理を行う範囲内には、大きさの異なる複数の位相検波信号（または相関信号、以下、位相検波信号に言及するが相関信号も同様である）がある。これら位相検波信号の加算平均処理を行い平均位相検波信号を求める場合、平均位相検波信号の角度におよぼす影響は、各位相検波信号の大きさに依存する。すなわち、平均処理範囲内において、相対的に小さい位相検波信号は、位相検波信号の角度に与える影響は小さい。
【００２８】
このため、空間平均処理の対象となった複数の位相検波信号からなる位相検波信号群の位置、あるいは、複数の位相検波信号から求められた平均位相検波信号の位置を各位相検波信号の大きさに応じて決定する。つまり、複数の位相検波信号の重心位置を求め、求めた位置において平均位相検波信号が得られたものとする。例えば、等間隔に並んだ連続した５点（Ｐ１、Ｐ２、Ｐ３、Ｐ４、Ｐ５）からなる群において、この群内の位相検波信号データを空間平均処理する場合、群の位置あるいは得られた位相検波信号の位置は５点の中央にある点Ｐ３とするのが一般的である。しかし、本発明では、Ｐ１のデータとＰ２のデータが他の３点のデータに比べて極端に小さい場合、データ平均値に及ぶＰ１のデータとＰ２のデータの影響は極めて小さいため、群あるいは平均位相検波信号の位置は中央のＰ３ではなく、重心のＰ４とする。これによって、測定対象とする組織の位置を正確に追跡することが可能となる。
【００２９】
このように各測定対象位置から得られる位相検波信号または相関信号を空間平均処理することによってノイズの影響を低減するとともに、各測定対象位置から得られる位相検波信号または相関信号の重心位置を用いて、測定対象とする組織の位置を正確に追跡し、また、弾性率を求めるための任意の２つの測定位置間の距離（厚さ）を正確に求めることが可能となる。したがって、被検体の性状特性、特に組織の歪み量や弾性率を正確に求めることができる。
【００３０】
以下、本発明による超音波診断装置の実施形態を図面を用いて詳細に説明する。
【００３１】
図１は、本実施形態の超音波診断装置１１を用いて血管壁組織の性状診断を行うための構成を示すブロック図である。超音波診断装置１１に接続された超音波プローブ１３は被験者の体表２に密着するよう支持され、血管外組織１の内部へ超音波を送信する。送信された超音波は血管３や血液５にて反射、散乱し、その一部が超音波プローブ１３へ戻り、超音波反射波（エコー）として受信される。超音波診断装置１１は、受信信号の解析、演算を行い、血管前壁４の形状情報や運動情報を求める。また、超音波診断装置１１には血圧計１２が接続されており、血圧計１２が測定した被験者の血圧データは超音波診断装置１１へと入力される。超音波診断装置１１は、例えば特許文献１に開示されている方法にしたがって、検波信号の振幅および位相の両方を用い、制約付き最小二乗法によって対象の瞬時の位置を決定する。高精度な（位置変化量の測定精度は±０．２ミクロン）位相トラッキングを行うことにより、血管前壁４における微小部位の厚さや厚さ変化量の時間変化の様子を十分な精度で測定することができる。さらに、血圧計１２から得た血圧データを用いることで、血管前壁４における微小部位の弾性特性を求めることができる。また、超音波診断装置１１には心電計２２が接続されており、心電計２２が測定した心電波形は超音波診断装置１１へと入力され、データ取得やデータリセットのタイミングを決定するトリガ信号として使用される。心電計２２は他の生体信号検出手段である心音計や脈波計と置き換えることも可能であり、心電波形の替わりに心音波形や脈波波形をトリガ信号として用いることも可能である。
【００３２】
以下、超音波診断装置１１の構成および動作を詳細に説明する。図２は、超音波診断装置１１の構成を示すブロック図である。超音波診断装置１１は、送信部１４、受信部１５、遅延時間制御部１６、位相検波部１７、フィルタ部１８、演算部１９、演算データ記憶部２０および表示部２１を備えている。また、コンピュータなどからなる制御部３０がこれら各部を制御する。
【００３３】
送信部１４は、所定の駆動パルス信号を生成し、超音波プローブ１３に出力する。駆動パルス信号により超音波プローブ１３から送信される超音波送信波は、血管３等の生体組織において反射、散乱し、生じた超音波反射波が超音波プローブ１３で受信される。超音波を発生させる駆動パルスの周波数は、時間軸上で隣接している前後の超音波パルスが重ならないように、測定対象の深さと超音波の音速とを考慮して決定される。
【００３４】
受信部１５は超音波プローブ１３を用いて超音波反射波を検出し、検出によって得られた信号を増幅することにより、受信信号を生成する。受信部１５はＡ／Ｄ変換部を含み、受信信号をさらにデジタル信号に変換する。送信部１４および受信部１５は電子部品などを用いて構成される。
【００３５】
遅延時間制御部１６は送信部１４および受信部１５に接続されており、送信部１４から超音波プローブ１３の超音波振動子群に与える駆動パルス信号の遅延時間を制御する。これにより、超音波プローブ１３から送信される超音波送信波の超音波ビームの音響線の方向や焦点深度を変化させる。また、超音波プローブ１３によって受信され、受信部１５によって増幅された受信信号の遅延時間を制御することにより、開口径を変化させたり、焦点位置を変化させたりすることができる。遅延時間制御部１６の出力は位相検波部１７に入力される。
【００３６】
位相検波部１７は、遅延時間制御部１６で遅延制御された受信信号を位相検波し、実部信号と虚部信号とからなる位相検波信号を作成する。位相検波信号はフィルタ部１８に入力される。フィルタ部１８は、高周波成分、測定対象以外からの反射成分およびノイズ成分等を除去する。位相検波部１７およびフィルタ部１８はソフトウエアによってもハードウエアによっても構成することができる。
【００３７】
位相検波信号の実部信号および虚部信号は、演算部１９へ入力される。図３は演算部１９の構成を詳細に示すブロック図である。演算部１９は、形状測定値演算部３１および性状測定値演算部３２を含む。演算部１９は、ソフトウエアによってもハードウエアによっても構成することができる。
【００３８】
形状測定値演算部３１は、受信信号の実部信号および虚部信号に基づいて、複数の測定対象位置における生体組織の運動速度をそれぞれ求め、運動速度を積分することによって、位置変位量（位置の時間変位量）を求める。求めた位置変位量から各測定対象位置間における生体組織の厚さ変化量（伸縮量）を求める。また、心電計２２から一心周期に関する情報を受け取り、一心周期における厚さ変化量の最大値と最小値との差分である最大厚さ変化量と厚さの最大値を求める。
【００３９】
性状測定値演算部３２は、最大厚さ変化量と厚さの最大値を受け取り、生体組織の歪みを求める。さらに血圧計１２から得られる血圧データを用いて、各測定対象位置間の組織の弾性特性を求める。
【００４０】
このようにして求められた生体組織の最大厚さ変化量、歪み量、あるいは弾性特性は、計測領域に対応してマッピングされ、形状測定値または性状測定値の空間分布を示す心周期毎の空間分布フレームとして表示部２１へ出力される。
【００４１】
図４および図５を参照して、演算部１９におけるこれらの演算をさらに詳しく説明する。図４は、生体６０を伝播する超音波ビーム６７を模式的に示しており、図では血管以外の生体組織６２、血管壁６３および血液６４が示されている。生体６０の表面に配置された超音波プローブ１３から送信した超音波送信波は、生体６０中を進行する。超音波送信波は、ある有限の幅を持つ超音波ビーム６７として生体６０中を伝播し、その過程において生体組織６２、血管壁６３そして血液６４等によって反射または散乱し、反射または散乱された超音波の一部が超音波プローブ１３へ戻り、超音波反射波として受信される。超音波反射波は離散的な時系列信号ｒｋ（ｔ）として検出され、超音波プローブ１３に近い組織から得られる反射の時系列信号ほど、時間軸上で原点近くに位置する。超音波ビーム６７の幅（ビーム径）は、遅延時間を変化させることにより制御することができる。
【００４２】
超音波ビームの中心軸である音響線６６上に位置する複数の測定対象位置Ｐｎ（Ｐ１、Ｐ２、Ｐ３、Ｐｋ・・・Ｐｎ、ｎは３以上の自然数）は、測定開始時には、ある一定間隔Ｌで超音波プローブ１３に近い順にＰ１、Ｐ２、Ｐ３、Ｐｋ・・・Ｐｎと配列している。反射波信号ｒｋ（ｔ）を位相検波部１７において予め定めた検波周波数を用いて直交検波し、実部信号および虚部信号からなる位相検波信号Ｒｋ（ｔ）を求め、フィルタ部１８を通過させる。演算部１９の形状測定値演算部３１では、位相検波信号Ｒｋ（ｔ）と微小時間Δｔ後の位相検波信号Ｒｋ（ｔ＋Δｔ）において振幅は変化せず、位相および反射位置のみが変化するという制約のもとで、位相検波信号Ｒｋ（ｔ）とＲｋ（ｔ＋Δｔ）との相関演算を行い、位置ｋの相関信号Ｑｋ（ｔ）を以下の式から求める。
【００４３】
Ｑｋ（ｔ）＝Ｒｋ（ｔ＋Δｔ）×Ｒｋ*（ｔ）    （＊は複素共役を示す）
【００４４】
以下において詳細に説明するように本発明では、ここで求めた位相検波信号Ｒｋ（ｔ）または相関信号Ｑｋ（ｔ）に対して空間平均処理を行う。その後、相関信号Ｑｋ（ｔ）の実部と虚部とが成す角度θｋ（ｔ）を求め、角度θｋ（ｔ）から測定対象位置Ｐｋの時間変化の様子、すなわち位置Ｐｋの位置変位量ｈｋ（ｔ）を求めることができる。
【００４５】
図５は、測定対象位置Ｐｎと弾性特性演算の対象組織Ｔｎとの関係を示している。対象組織Ｔｋは、隣接する測定対象位置ＰｋとＰｋ＋１とに挟まれた範囲に厚さＬを有して位置している。ｎ個の測定対象位置Ｐ１・・・・Ｐｎからは（ｎ−１）個の対象組織Ｔ１・・・・Ｔｎ−１を設けることができる。
【００４６】
対象組織Ｔｋの伸縮量である厚さ変化量Ｈｋ（ｔ）は、測定対象位置ＰｋとＰｋ＋１の位置変位量ｈｋ（ｔ）とｈｋ＋１（ｔ）とから、Ｈｋ（ｔ）＝ｈｋ＋１（ｔ）−ｈｋ（ｔ）として求められる。
【００４７】
血管壁６４の組織Ｔｋの厚さの変化は、心拍による血圧の変化に応じて生じ、心周期におよそ同期して繰り返される。したがって、弾性特性も心周期に同期して一心周期毎の数値を求めることが好適である。一心周期内の厚さ変化量Ｈｋ（ｔ）から最大値と最小値とを抽出し、最大値と最小値との差分を最大厚さ変化量ΔＨｋとする。また、血圧の最大値と最小値との差分を脈圧Δｐとする。対象組織の厚さの最大値（あるいは初期値）をＨｍとしたとき、ＨｍはＬ×｛（ｋ＋１）−ｋ｝＝Ｌなので、歪み量Ｓｋおよび弾性特性χｋはそれぞれ次式で求めることができる。
【００４８】
Ｓｋ＝ΔＨｋ／Ｈｍ＝ΔＨｋ／Ｌ
χｋ＝Δｐ／Ｓｋ＝Δｐ・Ｈｍ／ΔＨｋ
【００４９】
測定対象位置Ｐｎの数やその間隔は、測定の目的や測定対象物である生体組織の特性に応じて任意に設定できる。
【００５０】
次に空間平均処理を詳細に説明する。上述したように、測定された反射波信号ｒｋ（ｔ）は測定環境に応じて様々なノイズを含んでおり、ここから求めた位置変位量ｈｋ（ｔ）や厚さ変化量Ｈｋ（ｔ）もノイズを含んでしまうことがある。厚さの最大値Ｈｍは一心周期内の厚さ変化量Ｈｋ（ｔ）の最大値と最小値の差分であるため、ノイズの影響を大きく受けてしまい、結果として得られる歪み量Ｓｋや弾性特性χｋは誤った値となってしまう。
【００５１】
こうしたノイズによる影響を低減するために、超音波診断装置１１は、信号の空間平均処理を行う。反射波信号ｒｋ（ｔ）が含むノイズを除去するためには、位相検波信号Ｒｋ（ｔ）や相関信号Ｑｋ（ｔ）の空間平均処理が有効である。以下の説明では、図６および図７を参照して、相関信号Ｑｋ（ｔ）の空間平均処理を行う例を説明する。図６は演算部１９における信号処理を示すフローチャートである。演算部１９は、フィルタ部１８から出力された位相検波信号Ｒｋ（ｔ）を相関演算処理し、相関信号Ｑｋ（ｔ）を生成する（ステップＳ１）。さらに、相関信号Ｑｋ（ｔ）を空間平均処理し、平均相関信号Ｑ’ｋ（ｔ）を生成する（ステップＳ２）。
【００５２】
この空間平均処理を図７を用いて詳細に説明する。ある時刻ｔにおいて、少なくとも２つの測定対象位置を含む領域を設定する。好ましくは、２以上１０以下の測定対象位置を含む。この領域は同一音響線上の複数の測定対象位置を含むことが好ましい。例えば、位置Ｐｋの相関信号Ｑｋ（ｔ）を含む同一音響線上の複数の位置を含む領域を設定する。この領域に含まれる複数の位置における複数この相関信号を相関信号群とする。例えば、図７に示すように、位置Ｐｋを中心とした連続する５つの位置Ｐｋ−２、Ｐｋ−１、Ｐｋ、Ｐｋ＋１、Ｐｋ＋２の相関信号Ｑｋ−２（ｔ）、Ｑｋ−１（ｔ）、Ｑｋ（ｔ）、Ｑｋ＋１（ｔ）、Ｑｋ＋２（ｔ）を相関信号群とする。次にこれら信号群の空間平均処理、もしくは加算処理を行い、これを位置Ｐｋの平均相関信号Ｑ’ｋ（ｔ）とする。相関信号Ｑｋ（ｔ）は複素数であり、平均処理や加算処理はベクトル演算となる。測定対象位置の組み合わせが異なるように領域をシフトさせて設定し、この空間平均処理をすべての位置Ｐｋについて求める。例えば平均相関信号Ｑ’ｋ＋１（ｔ）は相関信号群Ｑｋ−１（ｔ）、Ｑｋ（ｔ）、Ｑｋ＋１（ｔ）、Ｑｋ＋２（ｔ）、Ｑｋ＋３（ｔ）から求められる。
【００５３】
以上のように求めた平均相関信号Ｑ’ｋ（ｔ）の位相角を用いて位相トラッキングを行い、位置変位量ｈ’ｋ（ｔ）を演算する（ステップＳ３）。そして、位置Ｐｋの位置変位量ｈ’ｋ（ｔ）と位置Ｐｋ＋ｄの位置変位量ｈ’ｋ＋ｄ（ｔ）との差分から厚さ変化量Ｈｋ（ｔ）と最大厚さ変化量ΔＨｋとを求め（ステップＳ４）、続いて歪み量Ｓｋを求める（ステップＳ５）。
【００５４】
ここで歪み量Ｓｋは、一心周期内の最大厚さ変化量ΔＨｋと対象組織の厚さの初期値Ｈｍから求められる。一般的にはＨｍは位置がｄ離れている二点間の距離初期値なので、Ｌ×｛（ｋ＋ｄ）−ｋ｝＝Ｌ×ｄというある一定の数となる。しかし、本発明では、相関信号群を構成する各相関信号の大きさに応じて平均相関信号の位置座標初期値を別途求め、求めた位置座標初期値を用いて二点間の距離初期値を算出する。具体的には、平均相関信号の位置座標初期値は、相関信号群を構成する各相関信号の位置の重心とする。なぜならば、空間平均処理を求めるベクトル演算において、各相関信号が結果に与える影響は、その大きさ（絶対値）に依存するからである。より具体的には、座標初期値Ｚｋ−２、Ｚｋ−１、Ｚｋ、Ｚｋ＋１、Ｚｋ＋２を有する位置Ｐｋ−２、Ｐｋ−１、Ｐｋ、Ｐｋ＋１、Ｐｋ＋２の相関信号Ｑｋ−２（ｔ）、Ｑｋ−１（ｔ）、Ｑｋ（ｔ）、Ｑｋ＋１（ｔ）、Ｑｋ＋２（ｔ）からなる相関信号群から求めた平均相関信号Ｑ’ｋ（ｔ）の位置座標初期値Ｚ’ｋは、次式で求める。
【００５５】
【数１】
2007020724.tif 000003
【００５６】
上式を用い、すべての位置Ｐｋにおける平均相関信号Ｑ’ｋ（ｔ）の位置座標初期値Ｚ’ｋを求める（ステップＳ７）。すなわち、例えば平均相関信号Ｑ’ｋ＋ｄ（ｔ）の位置座標初期値Ｚ’ｋ＋ｄは相関信号群Ｑｋ＋ｄ−２（ｔ）、Ｑｋ＋ｄ−１（ｔ）、Ｑｋ＋ｄ（ｔ）、Ｑｋ＋ｄ＋１（ｔ）、Ｑｋ＋ｄ＋２（ｔ）、から求められる。求めた位置座標初期値Ｚ’ｋは、小数を含む形式のままでも良いし、四捨五入した整数値としてもかまわないが、四捨五入すると精度は低くなる。
【００５７】
続いて、ステップＳ７にて求めた平均相関信号の位置座標初期値Ｚ’ｋを用いて、歪み量や弾性特性演算の対象組織の厚さ初期値Ｈｍ’ｋを求める（ステップＳ８）。ステップＳ４においては、位置がｄ離れた二点間の厚さ変化量を求めており、ステップＳ８においても位置がｄ離れた二点間の厚さの初期値を求める。位置Ｐｋの位置座標初期値Ｚ’ｋと位置Ｐｋ＋ｄの位置座標初期値Ｚ’ｋ＋ｄとの差分から、厚さ初期値Ｈｍ’ｋ＝Ｚ’ｋ＋ｄ−Ｚ’ｋとして求められる。
【００５８】
そして、最大厚さ変化量ΔＨｋと厚さ初期値Ｈｍ’ｋとから歪み量Ｓｋを求め（ステップＳ５）、続いて血圧計１２で測定された脈圧Δｐと歪み量Ｓｋとから弾性特性χｋが演算される（ステップＳ６）。
【００５９】
以上、詳述したように、相関信号群の空間平均処理により平均相関信号を求める場合には、平均相関信号を求めた複数の相関信号の重心位置を求め、平均相関信号の位置座標初期値とすることにより、正確な歪み量や弾性特性を求めることが可能となる。本発明は、相関信号群が組織の境界部（例えば血管中膜と血管外膜、血液と血管内膜、など）を跨いでいるときに特に効果的である。それは、組織が異なると音響インピーダンスが異なり、相関信号群内の各相関信号の大きさにばらつきが顕著になるためである。
【００６０】
なお、相関信号群内の相関信号を単純に平均するだけでなく重み付け平均を行えば、より正確に歪み量や弾性特性を求めることが可能となる。重み付けに用いる係数は、相関信号群の中心にピークが位置するように設けられたガウス関数やｓｙｎｃ関数を用いることが好適であるが、例えば連続した５点に対して１：２：３：２：１となる係数を積算するような重み付けでも効果的である。またこの場合、位置の重心を求める過程にも上記重み付けを導入することが好適であるが、平均相関信号の演算と、位置の重心演算のどちらか一方のみに重み付けを行ってもよい。
【００６１】
さらに、ステップＳ７で求める位置座標初期値は、ステップＳ１で求めた相関信号Ｑｋ（ｔ）をもとに演算されるが、相関信号Ｑｋ（ｔ）の大きさのばらつき分布と位相検波信号Ｒｋ（ｔ）の大きさのばらつき分布とがほぼ等しいことから、図６において点線で示すように位置座標初期値を位相検波信号Ｒｋ（ｔ）をもとに演算してもよい。
【００６２】
また、上述の例では相関信号に空間平均処理を施しているが、位相検波信号に空間平均処理を施しても同様な効果が得られる。図８は、位相検波信号を空間平均処理するフローチャートを示している。フィルタ部１８から出力された位相検波信号Ｒｋ（ｔ）は、まず空間平均処理され、平均位相検波信号Ｒ’ｋ（ｔ）となり（ステップＳ１１）、平均位相検波信号Ｒ’ｋ（ｔ）が相関演算処理され相関信号Ｑｋ（ｔ）となる（ステップＳ１２）。この相関信号Ｑｋ（ｔ）の位相角を用いて位相トラッキングを行い、位置変位量ｈｋ（ｔ）を演算する（ステップＳ１３）。そして、位置Ｐｋの位置変位量ｈｋ（ｔ）と位置Ｐｋ＋ｄの位置変位量ｈｋ＋ｄ（ｔ）との差分から厚さ変化量Ｈｋ（ｔ）と最大厚さ変化量ΔＨｋとを求める（ステップＳ１４）。一方、位相検波信号Ｒｋ（ｔ）をもとに平均位相検波信号Ｒ’ｋ（ｔ）の位置座標初期値Ｚ’ｋを求め（ステップＳ１７）、歪み量や弾性特性演算の対象組織の厚さ初期値Ｈｍ’ｋを求める（ステップＳ１８）。そして、ステップＳ１４で求めた最大厚さ変化量ΔＨｋとステップＳ１８で求めた厚さ初期値Ｈｍ’ｋとを用いて歪み量Ｓｋを演算し（ステップＳ１５）、血圧計１２で測定された脈圧Δｐと歪み量Ｓｋとから弾性特性χｋが演算される（ステップＳ１６）。
【００６３】
また、本実施形態では血管壁組織の性状診断を行う場合を例として説明を行ったが、本発明の特徴は、音響インピーダンスが変化する境界部での組織性状診断をより正確に行える点にある。このため、本発明の超音波診断装置は、血管壁組織以外の種々の生体組織の性状測定を好適に行うことができる。さらに、生体以外のゲルやゴムの性状を調べる場合にも好適である。
【産業上の利用可能性】
【００６４】
本発明の超音波診断装置は、生体組織、特に血管壁の形状および性状特性を好適に測定できる。また、医師が生体組織の形状および性状特性を正確に診断することが可能な超音波診断装置として有用である。
【図面の簡単な説明】
【００６５】
【図１】本発明による超音波診断装置を用いて血管壁組織性状の診断を行うための構成を示すブロック図である。
【図２】本発明による超音波診断装置の実施形態の構成を示すブロック図である。
【図３】図２に示す超音波診断装置の演算部の構成を詳細に示すブロック図である。る。
【図４】血管壁を伝播する超音波ビームと測定対象位置とを模式的に示す図である。
【図５】測定対象位置と弾性特性を求める対象組織との関係を示す図である。
【図６】本発明による超音波診断装置の実施形態における演算部の信号処理を説明するフローチャートである。
【図７】測定対象位置と平均相関信号および位置座標初期値との関係を示す図である。
【図８】本発明による超音波診断装置の実施形態における演算部の信号処理を説明するフローチャートである。
【符号の説明】
【００６６】
１  血管外組織
２  体表
３  血管
４  血管前壁
５  血液
１１  超音波診断装置
１２  血圧計
１３  超音波プローブ
１４  送信部
１５  受信部
１６  遅延時間制御部
１７  位相検波部
１８  フィルタ部
１９  演算部
２０  演算データ記憶部
２１  表示部
２２  心電計
３１  形状測定値演算部
３２  性状特性値演算部
６０  生体
６２  生体組織
６３  血管壁
６４  血液
６６  音響線
６７  超音波ビーム
【図１】
2007020724.tif 000004
【図２】
2007020724.tif 000005
【図３】
2007020724.tif 000006
【図４】
2007020724.tif 000007
【図５】
2007020724.tif 000008
【図６】
2007020724.tif 000009
【図７】
2007020724.tif 000010
【図８】
2007020724.tif 000011