JP5497821B2

JP5497821B2 - 流体流速検出装置及びプログラム

Info

Publication number: JP5497821B2
Application number: JP2012060973A
Authority: JP
Inventors: 稔小野; 慶一板谷; 孝岡田
Original assignee: University of Tokyo NUC; Hitachi Aloka Medical Ltd
Current assignee: Hitachi Ltd; University of Tokyo NUC
Priority date: 2012-03-16
Filing date: 2012-03-16
Publication date: 2014-05-21
Anticipated expiration: 2032-03-16
Also published as: CN104168835A; EP2826425B1; US9267828B2; US20150013471A1; WO2013136573A1; EP2826425A1; JP2013192643A; EP2826425A4; CN104168835B

Description

本発明は、血管や心臓等、生体内部の管腔状組織内を流れる流体の流速を検出する流体流速検出装置及びプログラムに関する。

従来、超音波ビームを生体に向けて送受して得た超音波エコー信号から血流のドプラ信号を検出し，その信号から血管や心臓等の内部を流れる血液の血流速度を検出する方法として、例えば特許文献１に開示するものが知られている。この特許文献１に開示された方法では、血流を基本流成分と渦流成分とにわけて演算している。

また一方で、非特許文献１に開示された方法は、心臓の左心室（ＬＶ）内の血流速度を検出する例が開示されている。この方法は具体的には次のようなものである。すなわち図５に示すように、公知の超音波診断法により２次元カラードプラ画像１００と、連続的なＢモードの画像２００とを得ておき、２次元カラードプラの画像１００からは、超音波ビーム方向の速度（Radial Velocities）Ｖrの情報１０１を得る。一方、Ｂモードの画像２００からはスペックル・トラッキング等の広く知られた方法により左心室の壁の運動速度を検出する。

次に連続の方程式

を用いて、超音波ビーム方向に直交する角度方向の速度Ｖθの情報１０２を得る。
ここで

は、

を簡略に示したものである（以下同じ）。

具体的には、図６に示すように、超音波ビーム源Ｔを中心とした２次元極座標系（ｒ，θ）において、ある動径位置ｒにおける左心室の一方側の壁（例えば後壁：posterior wall）の運動速度をＶ_θ ^-（ｒ）、他方側の壁（例えば前壁：anterior wall）の運動速度をＶ_θ ⁺（ｒ）とする。

このとき、ある動径位置ｒ，角度位置θの点（以下、注目位置と呼ぶ）における、超音波ビーム方向に直交する角度方向の速度Ｖ_θは、（１）式を積分することによって得られる。この積分には、後壁側から積分していく第１の積分結果と、前壁側から積分していく第２の積分結果とがある。つまり、後壁の位置θ-から動径ｒが一定の線に沿って（角度方向に）注目位置θまで積分し、積分定数として後壁の運動速度Ｖ_θ ^-（ｒ）を用いる第１の積分は、

となる。また、前壁の位置θ⁺から動径ｒが一定の線に沿って（角度方向に）注目位置θまで積分し、積分定数として前壁の運動速度Ｖ_θ ⁺（ｒ）を用いる第２の積分は、

となる。

これは結局、一階の偏微分方程式に対して、後壁側と前壁側とで２つの境界条件が設定できる結果、余分な解が得られていることに相当する。そこでこれらを重みづけ平均することにより、解の誤差を低減する。すなわち、注目位置（ｒ，θ）における、超音波ビーム方向に直交する角度方向の速度Ｖ_θを、上記第１の積分結果と第２の積分結果とを重みｗで、重みづけ平均して、

とする。

ここにｗは重みであり、技術常識として、（２），（３）式のそれぞれの解V_θが一様に等しくなるように選択する。具体的には、積分経路に沿った位置を表すｒ及びθの関数として、

で定められる。また（５）式で分母が「０」に近くなり、発散する場合は

と、近似的な式で置換える。ここで∂_θV_θは（１）で得られるものである。

特許第４６９８６８６号公報

DamienGarcia, et. al., "Two-Dimensional Intraventricular Flow Mapping by DigitalProcessing Conventional Color-Doppler Echocardiography Images", IEEETransactions on Medical Imaging, Vol.29, No.10, October 2010, 1701-1712

しかしながら非特許文献１に開示の技術では、（５）または（６）式でV_θ（ｒ，θ+）を用いる場合にはθ^-側の壁側における重みが過大ないし過小となって壁の運動方向に直交する向きに流体の速度成分が大きくなってしまったり、壁の運動速度以上に、流体が壁から遠ざかる速度が大きくなったりするなど、不自然な結果が得られることになっており、医学的にみて演算結果の妥当性が低いという問題があった。

本発明は上記実情に鑑みて為されたもので、医学的見地から妥当な評価を得ることのできる流体流速検出装置及びプログラムを提供することを、その目的の一つとする。

上記従来例の問題点を解決するための本発明は、流体流速検出装置であって、生体内部の管腔状組織内を流れる流体に対し、複数の角度方向位置で超音波ビームを送受波し、超音波ビームごとにエコーデータ列を出力する動作を、繰返し行う超音波探触子と、前記エコーデータ列から、前記流体の前記超音波ビーム方向成分に係る流速を演算する第１演算手段と、前記繰返し得られたエコーデータ列から、前記管腔状組織の複数の時点でのＢモード画像を得て、当該複数の時点でのＢモード画像に基づいて、前記管腔状組織の前記複数の角度方向での一方側の壁と、他方側の壁との運動速度を検出する検出手段と、前記流体の前記超音波ビーム方向に直交する成分に係る流速と、前記超音波ビーム方向成分に係る流速の情報とに関連する連続の方程式を、管腔状組織の一方側から、流速を求めたい注目位置（ｒ，θ）まで積分し、当該積分における境界条件として、前記検出手段により検出した当該一方側の壁の運動速度を用いて、第１の積分結果として、前記流体の前記超音波ビーム方向に直交する成分に係る流速の計算値Ｖ_θ ^-を得るとともに、前記流体の前記超音波ビーム方向に直交する成分に係る流速と、前記超音波ビーム方向成分に係る流速の情報とに関連する連続の方程式を、管腔状組織の他方側から、流速を求めたい注目位置まで積分し、当該積分における境界条件として、前記検出手段により検出した当該他方側の壁の運動速度を用いて、第２の積分結果として、前記流体の前記超音波ビーム方向に直交する成分に係る流速の計算値Ｖ_θ ⁺を得て、これら第１の積分結果Ｖ^- _θと、第２の積分結果Ｖ⁺ _θと、重みｗとを用いて、Ｖ_θ（ｒ，θ）=ｗ・Ｖ_θ ^-＋（１−ｗ）・Ｖ_θ ⁺を演算し、前記流体の前記超音波ビーム方向に直交する成分に係る流速の計算値Ｖ_θ（ｒ，θ）を演算する第２の演算手段と、を含み、前記重みｗは、少なくとも前記一方側の壁から予め定めた距離までは、当該一方側の壁からの距離ｄに比例し、前記他方側の壁から予め定めた距離までは、当該他方側の壁からの距離ｄ′に比例する値であることとしたものである。また前記重みｗは、一方側の壁からの距離ｄに比例する値としてもよい。

また本発明の一態様に係るプログラムは、生体内部の管腔状組織内を流れる流体に対し、複数の角度方向位置で超音波ビームを送受波し、超音波ビームごとにエコーデータ列を出力する動作を、繰返し行う超音波探触子に接続されたコンピュータを、前記エコーデータ列から、前記流体の前記超音波ビーム方向成分に係る流速を演算する第１演算手段と、前記繰返し得られたエコーデータ列から、前記管腔状組織の複数の時点でのＢモード画像を得て、当該複数の時点でのＢモード画像に基づいて、前記管腔状組織の前記複数の角度方向での一方側の壁と、他方側の壁との運動速度を検出する検出手段と、前記流体の前記超音波ビーム方向に直交する成分に係る流速と、前記超音波ビーム方向成分に係る流速の情報とに関連する連続の方程式を、管腔状組織の一方側から、流速を求めたい注目位置（ｒ，θ）まで積分し、当該積分における境界条件として、前記検出手段により検出した当該一方側の壁の運動速度を用いて、第１の積分結果として、前記流体の前記超音波ビーム方向に直交する成分に係る流速の計算値Ｖ_θ ^-を得るとともに、前記流体の前記超音波ビーム方向に直交する成分に係る流速と、前記超音波ビーム方向成分に係る流速の情報とに関連する連続の方程式を、管腔状組織の他方側から、流速を求めたい注目位置まで積分し、当該積分における境界条件として、前記検出手段により検出した当該他方側の壁の運動速度を用いて、第２の積分結果として、前記流体の前記超音波ビーム方向に直交する成分に係る流速の計算値Ｖ⁺ _θを得て、これら第１の積分結果Ｖ^- _θと、第２の積分結果Ｖ_θ ⁺と、重みｗとを用いて、Ｖ_θ（ｒ，θ）=ｗ・Ｖ_θ ^-＋（１−ｗ）・Ｖ_θ ⁺を演算し、前記流体の前記超音波ビーム方向に直交する成分に係る流速の計算値Ｖθ（ｒ，θ）を演算する第２の演算手段であって、前記重みｗは、少なくとも前記一方側の壁から予め定めた距離までは、当該一方側の壁からの距離ｄに比例し、前記他方側の壁から予め定めた距離までは、当該他方側の壁からの距離ｄ′に比例する値である第２の演算手段と、として機能させることとしたものである。

本発明によると、医学的見地から妥当な評価を得ることができる。

本発明の実施の形態に係る流体流速検出装置の構成例を表すブロック図である。本発明の実施の形態に係る流体流速検出装置による超音波ビームの送受波の方向の例を表す説明図である。本発明の実施の形態に係る流体流速検出装置の動作例を表すフローチャート図である。本発明の実施の形態に係る流体流速検出装置が用いる重みの例を表す説明図である。従来の血流速度を得る方法の例を表す概要図である。対象となる組織の壁面の運動速度の例を表す説明図である。

本発明の実施の形態について図面を参照しながら説明する。本実施の形態の流体流速検出装置は、図１に例示するように、超音波探触子１と、送受信部２と、送受信制御部３と、演算部４と、表示部５とを含んで構成される。ここで演算部４は、断層画像形成部４１と、カラードプラ演算部４２と、壁運動速度演算部４３と、速度ベクトル演算部４４と、重み設定部３４と、表示画像形成部４５とを含んで構成される。これら各部は、ＣＰＵ等のプログラム制御デバイスによりソフトウエア的に実現されてもよい。この場合、プログラム制御デバイスは、記憶部等コンピュータ可読な記録媒体に格納されたプログラムに従って動作し、上記各部の構成をソフトウエア的に実現する。

超音波探触子１は、例えばセクタ走査型またはコンベックス走査型のプローブである。送受信部２は、送受信制御部３の制御を受けて、この超音波探触子１を介して体表側から人体等の生体内部にある血管や心臓など、測定対象となる流体を内部に含んだ管腔状組織（以下、対象管腔状組織と呼ぶ）内を流れる流体に対し、複数の角度方向位置で超音波ビームＥを送受波する。
送受信制御部３の動作により、図２に例示するように、体表より超音波探触子１側にある仮想的な中心点を中心として、仮想的な面内（観測面内）で、所定角度範囲（θa＜θ＜θb）内で予め定められた複数の角度方向θ1，θ2，…，θn（θa≦θ1＜θ2＜…＜θn≦θb）に対して、それぞれ超音波ビームＥ（θj）（ｊ＝１，２，…，ｎ）が送信される。送信された超音波ビームは、体内の組織や流体の組成物（血液ならば例えば赤血球）で反射し、超音波探触子１へ戻る。送受信部２は、超音波探触子１を介して、一定の時間ごとに、複数回、到来する超音波受信信号データ（エコーデータ）ｅi（θj）＝｛（θj，ｔ1），（θj，ｔ2），…，（θj，ｔm）｝（ｊ＝１，２，…，ｎ、また｛＊｝は、＊の配列（データアレイ）を表す）を取得して、超音波ビームごとにエコーデータ列を得る。そしてこの送受信部２は、超音波ビームごとのエコーデータ列を、演算部４の断層画像形成部４１と、カラードプラ演算部４２とに出力する。送受信制御部３は、以上の動作を、予め定めたタイミングごとに繰返し行うよう送受信部２を制御する。そして送受信部２はｉ回目の繰返しのときに当該タイミングで取得したエコーデータ列ｅi（θ1）…ｅi（θn）（ｉ＝１，２，…）を出力する。

本実施の形態では、断層画像形成部４１が、送受信部が出力するエコーデータ列から、Ｂモード画像を形成する。またカラードプラ演算部４２は、送受信部が出力するエコーデータ列から、流体の超音波ビーム方向成分に係る流速を演算して、２次元カラードプラの画像を形成する。また壁運動速度演算部４３は、断層画像形成部から対象管腔状組織の複数の時点でのＢモード画像を受け入れ、当該複数の時点でのＢモード画像に基づいて、対象管腔状組織の一方側の壁と、他方側の壁との運動速度を検出する。

そして速度ベクトル演算部４４は、流体の超音波ビーム方向に直交する成分に係る流速と、超音波ビーム方向成分に係る流速の情報とに関連する連続の方程式を、対象管腔状組織の一方側から、流速を求めたい注目位置（ｒ，θ）まで積分し、当該積分における境界条件として、先に検出した当該一方側の壁の運動速度を用いて、第１の積分結果として、流体の超音波ビーム方向に直交する成分に係る流速の計算値Ｖ_θ ^-を得るとともに、流体の超音波ビーム方向に直交する成分に係る流速と、超音波ビーム方向成分に係る流速の情報とに関連する連続の方程式を、対象管腔状組織の他方側から、流速を求めたい注目位置まで積分し、当該積分における境界条件として、検出手段により検出した当該他方側の壁の運動速度を用いて、第２の積分結果として、流体の超音波ビーム方向（ｒ軸方向）に直交する成分に係る流速の計算値Ｖ⁺ _θを得る。

速度ベクトル演算部４４は、これら第１の積分結果Ｖ^- _θと、第２の積分結果Ｖ⁺ _θと、重みｗとを用いて、
Ｖ_θ（ｒ，θ）=ｗ・Ｖ_θ ^-＋（１−ｗ）・Ｖ_θ ⁺
を演算し、これより流体の超音波ビーム方向に直交する成分に係る流速の計算値Ｖ_θ（ｒ，θ）を得て出力する。

以下、これらの演算部４の各部の動作をより具体的に説明する。
断層画像形成部４１は、送受信部が出力するエコーデータ列ｅi（θ1）…ｅi（θn）（ｉ＝１，２，…）を用いて、Ｂモード画像を形成する。
カラードプラ演算部４２は、送受信部が出力するエコーデータ列ｅi（θ1）…ｅi（θn）（ｉ＝１，２，…）のうち、予め定めた条件で選択されるｐ回目のエコーデータ列ｅp（θ1）…ｅp（θn）を用いて、２次元カラードプラの画像を形成する。そしてカラードプラ演算部４２は、この２次元カラードプラの画像から対象管腔状組織内の座標（ｒ，θ）における流体の超音波ビーム方向成分（ｒ成分）に係る流速Ｖｒ（ｒ，θ）（Radial
Velocities）を演算する。この演算方法は広く知られているので、ここでの説明は省略する。なお、条件は例えばｐ＝１としてもよいし、測定開始から予め定めた時間Ｔが経過した時点でのエコーデータ列を用いる条件であってもよい。

壁運動速度演算部４３は形成されたＢモード画像から、スペックル・トラッキング等の方法によって、超音波探触子１（超音波の放射源）からの距離ｒにある、対象管腔状組織の一方側で流体に接する壁面（（ｒ，θ^-）の位置にあるものとする）の運動速度のうち、超音波ビーム方向に直交した（θ軸）方向の成分Ｖθ_-（ｒ）を検出する。

また、この壁運動速度演算部４３は、形成されたＢモード画像から、同様の方法で、超音波探触子１（超音波の放射源）からの距離ｒにある、対象管腔状組織の他方側で流体に接する壁面（（ｒ，θ⁺）の位置にあるものとする）の運動速度のうち、超音波ビーム方向に直交した（θ軸）方向の成分Ｖ_θ ⁺（ｒ）を検出する。

この壁運動速度演算部４３が処理するＢモード画像は、ＸＹ直交座標系上に描画されている。従ってスペックル・トラッキング等によって得られる対象管腔状組織の壁の運動速度は、ＸＹ直交座標系における値となる。壁運動速度演算部４３は、このＸＹ直交座標系における値を、超音波探触子１（超音波の放射源）の位置を原点とした、極座標系（ｒ，θ）に変換して出力する。

速度ベクトル演算部４４は、カラードプラ演算部４２が出力する２次元カラードプラ画像（超音波ビーム方向の血液の流速情報Ｖr（ｒ，θ）と、壁運動速度演算部４３が出力するＶ_θ ^-（ｒ）及びＶ_θ ⁺（ｒ）とを用いた連続の方程式

により、超音波ビーム方向に直交する角度方向の流体の速度Ｖ_θを得る。ここで

は、

を簡略に示したものである（以下同じ）。

このとき、ある動径位置ｒ，角度位置θの点（以下、注目位置と呼ぶ）における、超音波ビーム方向に直交する角度方向の速度Ｖ_θは、（１）式を積分することによって得られる。この積分には、対象管腔状組織の一方側の壁の位置θ^-（ｒ）から積分していく第１の積分結果と、対象管腔状組織の他方側の壁の位置θ⁺（ｒ）から積分していく第２の積分結果とがある。つまり、一方側の壁の位置θ^-（ｒ）から動径ｒが一定の線に沿って（角度方向に）注目位置θ（ｒ）まで積分し、積分定数として一方側の壁の運動速度Ｖ_θ ^-（ｒ）を用いる第１の積分は、

となる。また、他方側の壁の位置の位置θ⁺（ｒ）から動径ｒが一定の線に沿って（角度方向に）注目位置θ（ｒ）まで積分し、積分定数として前壁の運動速度Ｖ_θ ⁺（ｒ）を用いる第２の積分は、

となる。

速度ベクトル演算部４４は、これら（２），（３）式の積分を、予め定めた複数の動径方向位置ｒ1，ｒ2…ごとに演算して、各動径方向位置における第１，第２の積分結果Ｖ_θ ^-（ｒ1），Ｖ_θ ⁺（ｒ1），Ｖ_θ ^-（ｒ2），Ｖ_θ ⁺（ｒ2），…を得る。

これは結局、一階の偏微分方程式に対して、一方側及び他方側の壁の２箇所で合計２つの境界条件が設定できる結果、余分な解が得られていることに相当する。そこでこれらを重みづけ平均することにより、解の誤差を低減する。すなわち速度ベクトル演算部４４は、上記第１の積分結果Ｖ_θ ^-（ｒ）と第２の積分結果Ｖ_θ ⁺（ｒ）と（ｒ＝ｒ1，ｒ2…）を、重み設定部３４が出力する重みｗで重みづけ平均して、

を演算し、注目位置（ｒ，θ）における、超音波ビーム方向に直交する角度方向の速度Ｖ_θ（ｒ，θ）を得る。

これにより速度ベクトル演算部４４は、対象管腔状組織内（ｒ，θ）における流体の流速を、超音波ビーム方向の成分Ｖｒ（ｒ，θ）と、超音波ビーム方向に直交する成分Ｖ_θ（ｒ，θ）とに分けて取得することとなる。

重み設定部３４は、少なくとも一方側の壁から予め定めた距離までは、当該一方側の壁からの距離ｄに比例し、他方側の壁から予め定めた距離までは、当該他方側の壁からの距離ｄ′に比例する値であるように重みを設定する。一例として、本実施の形態のこの重み設定部３４が設定する重みｗは、次のようになる。まず、動径方向座標ｒi（ｒiはｒ1，ｒ2…のいずれか）における、対象管腔状組織の一方側の壁から他方側の壁までの積分経路（動径方向座標ｒ＝一定の経路）上の距離Ｄ（ｒi）は、Ｄ（ｒi）＝ｒi・｜θ⁺（ｒi）−θ^-（ｒi）｜で得られる（｜＊｜は、＊の絶対値を演算することを意味する）。また一方側の壁から注目位置（ｒi，θ）までの距離ｄは、ｄ＝ｒi・｜θ−θ^-（ｒi）｜で得られる。そこで重み設定部３４は、重みｗを、ｗ＝ｄ／Ｄ（ｒi）＝｜θ−θ^-（ｒi）｜／｜θ⁺（ｒi）−θ^-（ｒi）｜とする。これにより設定される重みは、注目位置の角度方向座標θに関して一次（距離ｄに比例）する値となる。この比例関係は動径方向位置ｒによらずに成立する。なお、このように、重みｗを、一方側からの壁からの距離ｄに比例するものとすると、当然にｗ＝ｄ／Ｄ（ｒi）＝（Ｄ（ｒi）−ｒ′）／（ｒi）である（ここでｒ′は、他方側から注目位置までの積分経路上の距離）ので、他方側の壁からの距離ｄ′にも比例することとなる。

この重みは、（５）または（６）式により設定される重みとは異なり、積分経路上の各点で第１，第２の積分結果を一致させるものではない。しかしながら一見技術常識に反するこの重みの設定により、両側の壁近傍における流体の流速は、壁の運動速度を自然に反映したものとなることがわかった。
表示画像形成部４５は、断層画像形成部４１が出力するＢモード画像に、速度ベクトル演算部４４にて演算した、（ｒ，θ）座標の各点での速度ベクトルを表す画像を合成して、表示用画像を形成して出力する。

本実施の形態の流体流速検出装置は、以上の構成を有してなり、次のように動作する。具体的に以下の例では、心臓の左心室（ＬＶ）内の血流速度を検出する場合を例とする。この例では、超音波探触子１は、体表側から人体等の生体内部にある心臓の左心室内を流れる血液に対し、複数の角度方向位置で超音波ビームＥを送波する。そして超音波探触子１は、送波した超音波ビームごとに得られるエコーデータ列ｅ（θ1）…ｅ（θn）を送受信部２を介して演算部４に対して出力する。超音波探触子１は、この動作を、予め定めたタイミングごと（例えば間隔Δｔのタイミングごと）に繰返し行う。ｋ番目の繰返しで得られるエコーデータ列を、ｅk（θ1）…ｅk（θn）（ｋ＝１，２，…）とする。

速度ベクトル演算部４４は、図３に例示するように、超音波探触子１が出力するエコーデータ列ｅk（θ1）…ｅk（θn）（ｋ＝１，２，…）のうち、予め定めた条件で選択されるｐ回目のエコーデータ列ｅp（θ1）…ｅp（θn）から演算される２次元カラードプラの画像を得る（Ｓ１）。また壁運動速度演算部４３は、超音波探触子１が出力するエコーデータ列から、複数のＢモードの画像を得る（Ｓ２）。そして２次元カラードプラの画像から、複数の予め定められた点（規定点と呼ぶ）（ｒi，θj）（ここでｉ＝１，２，…、ｊ＝１，２…）での血液の超音波ビーム方向の速度（Radial Velocities）Ｖｒ（ｒi，θj）を得る（Ｓ３）。

一方、壁運動速度演算部４３は、処理Ｓ２で得た複数のＢモードの画像の各々からスペックル・トラッキング等の広く知られた方法により対象管腔状組織である左心室のｒ＝ｒi（ｉ＝１，２，…）における後壁（posterior wall）の位置θ^-（ｒi）と、前壁（anterior wall）の位置θ⁺（ｒi）とを求める。この際、スペックル・トラッキング等の方法では画像のＸ，Ｙ座標で得られる位置の情報を、超音波源である超音波探触子１の位置を原点とした極座標に変換して、これらの位置を求める。

また壁運動速度演算部４３は、この各壁の位置の情報と、複数のＢモードの画像のそれぞれが得られたタイミングの情報（例えばΔｔ間隔）を用いて、動径位置ｒi（ｉ＝１，２，…）のそれぞれにおける左心室の一方側の壁（例えば後壁：posterior wall）の運動速度Ｖ_θ ^-（ｒi）と、他方側の壁（例えば前壁：anterior wall）の運動速度Ｖ_θ ⁺（ｒi）とを求める（Ｓ４）。

次に重み設定部３４は、規定点のある動径方向位置ｒ＝ｒi（ｉ＝１，２，…）のそれぞれにおいて、（１）式の積分経路の一方端（後壁側）の位置θ^-（ｒi）と他方端（前壁側）の位置θ⁺（ｒi）とを得る。そして積分経路上の規定点（ｒi，θj）（ｊ＝１，２，…）において、一方端からの距離ｄ＝ｒi・｜θ（ｒi）−θ^-（ｒi）｜と、積分経路全体の長さＤ＝ｒi・｜θ⁺（ｒi）−θ^-（ｒi）｜とを用いて、重みｗ＝ｄ／Ｄ（ｒi）＝｜θ（ｒi）−θ^-（ｒi）｜／｜θ⁺（ｒi）−θ^-（ｒi）｜を演算する（Ｓ５）。

また速度ベクトル演算部４４は、この重みｗを用いて（４）式を演算し、各規定点（ｒi，θj）（ｉ＝１，２，…、ｊ＝１，２…）での、超音波ビーム方向に直交する角度方向の血液の速度Ｖθ（ｒi，θj）を得る（Ｓ６）。

これにより、各規定点（ｒi，θj）（ｉ＝１，２，…、ｊ＝１，２…）での、血液の速度Ｖが、超音波ビーム方向の速度成分Ｖｒ（ｒi，θj）と、超音波ビーム方向に直交する角度方向の速度成分Ｖθ（ｒi，θj）とにわけて得られたこととなる。速度ベクトル演算部４４は、この得られた各規定点での速度Ｖの値を出力する（Ｓ７）。

速度ベクトル演算部４４は、例えば以上の処理を、処理Ｓ１においてｐ＝１，２，…と予め定めた値ずつインクリメントしながら繰返し実行する。そうすると、各繰返しで得られるエコーデータ列ｅk（θ1）…ｅk（θn）（ｋ＝１，２，…）に基づいて、各エコーデータ列が得られた時点でのＢモード画像と、血液の速度Ｖとが得られることとなる。速度ベクトル演算部４４は、ｋ番目のエコーデータ列に基づいて各規定点で得られた速度Ｖk（ｒi，θj）について、ｋ番目のエコーデータ列に基づいて得られたＢモード画像上で各規定点（ｒi，θj）の位置を根元とし、この速度ベクトルＶk（ｒi，θj）の長さ及び方向を表す矢印図形を描画してもよい。このためには、演算部４の表示画像形成部４５は、断層画像形成部４１が出力するＢモード画像に、速度ベクトル演算部４４にて演算した、（ｒ，θ）座標の各点での速度ベクトルを表す画像を合成して、表示用画像を形成して出力する。つまり、断層画像形成部４１が出力するＢモード画像上で各規定点に対応する位置（Ｘ，Ｙ座標系上の位置）Ｐ（ｘ，ｙ）を演算するとともに、速度ベクトルＶk（ｒi，θj）もまた、Ｘ，Ｙ座標系上での各成分（Ｘ軸方向成分とＹ軸方向成分）の値Ｖk_x（ｘ，ｙ）（Ｘ軸方向成分）及びＶk_y（ｘ，ｙ）（Ｙ軸方向成分）に変換する。そして表示画像形成部４５は、予め定めた係数α（拡大縮小率）を乗じて、位置（α・Ｖk_x（ｘ，ｙ），α・Ｖk_y（ｘ，ｙ））を得ておき、規定点（ｘ，ｙ）からここで得た位置（α・Ｖk_x（ｘ，ｙ），α・Ｖk_y（ｘ，ｙ））までの矢印図形を描画する。

これを各繰返しで得られるエコーデータ列ｅk（θ1）…ｅk（θn）（ｋ＝１，２，…）に基づいて得られた画像や値ごとに繰返すと、心臓の動きとともに心臓内における血液の流れがアニメーションのコマとして描画されることとなる。表示画像形成部４５は、各コマの画像を連続的に置換えながら表示することで、アニメーションとして表示してもよいし、各コマの画像をマトリクス状に配列して表示してもよい。

図４に、（４）式を演算する際の重みについて、従来の重みの例と、本実施の形態の重みの例とを比較して示す。この図４に示されるように従来例の重み（図中破線にて示す）は、動径方向位置（図中のｒ1，ｒ2，ｒ3）により異なり、また、角度方向に沿って壁近傍でも必ずしも積分経路全体の長さのうち、壁からの距離の割合に対して比例的に変化する重みとなっていなかった。これに対して、本実施の形態の重み（図中実線で示す）は、少なくとも各壁から予め定めた範囲（積分範囲の両端から予め定めた範囲）では動径方向の位置によらず、少なくとも壁近傍において、積分経路全体の長さのうち、壁からの距離の割合に対して比例的に変化する重みとした。

この本実施の形態の重みは、技術常識的には、一方側の壁からの積分結果と、他方側の壁からの積分結果とを積分経路上の各点で一致させないため採用しがたい重みであったが、本実施の形態では、これを敢えて採用したものである。そしてこの本実施の形態におけるような重みによれば、医学的見地から妥当な評価を得ることができることが見出された。

なお、本実施の形態において、重みは、各壁から予め定めた距離ｌまでｗ＝ｄ／Ｄ（ｒi）＝｜θ（ｒi）−θ^-（ｒi）｜／｜θ⁺（ｒi）−θ^-（ｒi）｜にて演算された値であれば、各壁から予め定めた距離ｌより離れた位置では、この重みに従わなくてもよい。例えば、ｒ｜θ−θ⁺｜とｒ｜θ−θ^-｜との双方が、予め定めた距離ｌより大きい位置では（５）式（またはそれを近似した（６）式）とし、ｒ｜θ−θ⁺｜とｒ｜θ−θ^-｜とのいずれか一方が予め定めた距離ｌより小さい位置では重みｗ＝ｄ／Ｄ（ｒi）＝｜θ（ｒi）−θ^-（ｒi）｜／｜θ⁺（ｒi）−θ^-（ｒi）｜としてもよい。すなわち

としてもよい。

なお、ｒ｜θ−θ⁺｜とｒ｜θ−θ^-｜との双方が、予め定めた距離ｌより大きい範囲における重みは、（５）ないし（６）式に限られず、｜θ−θ^-｜に関して実験的に定められる単調減少関数とすることとしてもよい。

１超音波探触子、２送受信部、３送受信制御部、４演算部、５表示部、３４重み設定部、４１断層画像形成部、４２カラードプラ演算部、４３壁運動速度演算部、４４速度ベクトル演算部、４５表示画像形成部。

Claims

生体内部の管腔状組織内を流れる流体に対し、複数の角度方向位置で超音波ビームを送受波し、超音波ビームごとにエコーデータ列を出力する動作を、繰返し行う超音波探触子と、
前記エコーデータ列から、前記流体の前記超音波ビーム方向成分に係る流速を演算する第１演算手段と、
前記繰返し得られたエコーデータ列から、前記管腔状組織の複数の時点でのＢモード画像を得て、当該複数の時点でのＢモード画像に基づいて、前記管腔状組織の前記複数の角度方向での一方側の壁と、他方側の壁との運動速度を２次元座標系内で検出する検出手段と、
前記流体の前記超音波ビーム方向に直交する成分に係る流速と、前記超音波ビーム方向成分に係る流速の情報とに関連する連続の方程式を、管腔状組織の一方側から、流速を求めたい注目位置（ｒ，θ）まで積分し、当該積分における境界条件として、前記検出手段により検出した当該一方側の壁の運動速度を用いて、第１の積分結果として、前記流体の前記超音波ビーム方向に直交する成分に係る流速の計算値Ｖ^- _θを得るとともに、
前記流体の前記超音波ビーム方向に直交する成分に係る流速と、前記超音波ビーム方向成分に係る流速の情報とに関連する連続の方程式を、管腔状組織の他方側から、流速を求めたい注目位置まで積分し、当該積分における境界条件として、前記検出手段により検出した当該他方側の壁の運動速度を用いて、第２の積分結果として、前記流体の前記超音波ビーム方向に直交する成分に係る流速の計算値Ｖ⁺ _θを得て、
これら第１の積分結果Ｖ^- _θと、第２の積分結果Ｖ⁺ _θと、重みｗとを用いて、
Ｖ_θ（ｒ，θ）=ｗ・Ｖ^- _θ＋（１−ｗ）・Ｖ⁺ _θ
を演算し、前記流体の前記超音波ビーム方向に直交する成分に係る流速の計算値Ｖθ（ｒ，θ）を演算する第２の演算手段と、
を含み、
前記重みｗは、少なくとも前記一方側の壁から予め定めた距離までは、当該一方側の壁からの距離ｄに比例し、前記他方側の壁から予め定めた距離までは、当該他方側の壁からの距離ｄ′に比例する値である流体流速検出装置。
請求項１記載の流体流速検出装置であって、
前記重みｗは、一方側の壁からの距離ｄに比例する値である流体流速検出装置。
生体内部の管腔状組織内を流れる流体に対し、複数の角度方向位置で超音波ビームを送受波し、超音波ビームごとにエコーデータ列を出力する動作を、繰返し行う超音波探触子に接続されたコンピュータを、
前記エコーデータ列から、前記流体の前記超音波ビーム方向成分に係る流速を演算する第１演算手段と、
前記繰返し得られたエコーデータ列から、前記管腔状組織の複数の時点でのＢモード画像を得て、当該複数の時点でのＢモード画像に基づいて、前記管腔状組織の前記複数の角度方向での一方側の壁と、他方側の壁との運動速度を２次元座標系内で検出する検出手段と、
前記流体の前記超音波ビーム方向に直交する成分に係る流速と、前記超音波ビーム方向成分に係る流速の情報とに関連する連続の方程式を、管腔状組織の一方側から、流速を求めたい注目位置（ｒ，θ）まで積分し、当該積分における境界条件として、前記検出手段により検出した当該一方側の壁の運動速度を用いて、第１の積分結果として、前記流体の前記超音波ビーム方向に直交する成分に係る流速の計算値Ｖ^- _θを得るとともに、
前記流体の前記超音波ビーム方向に直交する成分に係る流速と、前記超音波ビーム方向成分に係る流速の情報とに関連する連続の方程式を、管腔状組織の他方側から、流速を求めたい注目位置まで積分し、当該積分における境界条件として、前記検出手段により検出した当該他方側の壁の運動速度を用いて、第２の積分結果として、前記流体の前記超音波ビーム方向に直交する成分に係る流速の計算値Ｖ⁺ _θを得て、
これら第１の積分結果Ｖ_θ ^-と、第２の積分結果Ｖ_θ ⁺と、重みｗとを用いて、
Ｖ_θ（ｒ，θ）=ｗ・Ｖ_θ ^-＋（１−ｗ）・Ｖ_θ ⁺
を演算し、前記流体の前記超音波ビーム方向に直交する成分に係る流速の計算値Ｖ_θ（ｒ，θ）を演算する第２の演算手段であって、
前記重みｗは、少なくとも前記一方側の壁から予め定めた距離までは、当該一方側の壁からの距離ｄに比例し、前記他方側の壁から予め定めた距離までは、当該他方側の壁からの距離ｄ′に比例する値である第２の演算手段と、
として機能させるプログラム。