(57)【要約】【課題】本発明は、様々な血流条件下において血管壁面近傍の血流速度勾配を高精度に推定する技術を提供する。【解決手段】検査対象が反射するエコー信号から、検査対象の血管内の壁面における血流速度勾配の推定値を算出する超音波撮像装置およびその演算方法において、エコー信号から算出される、血管壁面から血管中心へ径方向にとった複数の測定点における、血管壁面に沿った方向の血流速度から血流速度勾配分布を算出し、算出した血流速度勾配分布において所定の範囲を算出し、算出した所定の範囲内の測定点における血流速度および血流速度勾配の値から血管壁面における血流速度勾配の推定値を算出する。【選択図】図１
2018011634.tif 000002
【特許請求の範囲】
【請求項１】
検査対象が反射するエコー信号を受信する受信部と、
前記受信部により受信された前記エコー信号を処理する信号処理部と、を備え、
前記信号処理部は、
前記検査対象の血管壁面から血管中心へ径方向に並んだ複数の測定点における、前記エコー信号から算出される前記血管壁面に沿った方向の血流速度の値から、血流速度勾配分布を算出する速度勾配分布演算部と、
前記血流速度勾配分布において所定の範囲を算出する範囲特定演算部と、
前記所定の範囲内の測定点における血流速度および血流速度勾配の値から、血管壁面における血流速度勾配の推定値を算出する壁面速度勾配演算部と、を備えることを特徴とする超音波撮像装置。
【請求項２】
請求項１に記載の超音波撮像装置において、
前記範囲特定演算部は、前記血流速度勾配分布の形状に基づいて棄却する領域を決定し、前記血流速度勾配分布に対して前記棄却する領域を除いた範囲として、前記所定の範囲を算出することを特徴とする超音波撮像装置。
【請求項３】
請求項２に記載の超音波撮像装置において、
前記範囲特定演算部は、前記血流速度勾配分布において、前記血管壁面から前記血管中心に向かって単調増加から単調減少に転じる点および上に凸から下に凸に転じる変曲点の少なくともいずれかを基準として、前記所定の範囲を算出することを特徴とする超音波撮像装置。
【請求項４】
請求項３に記載の超音波撮像装置において、
前記範囲特定演算部は、前記血管壁面における血流速度勾配の推定に用いる測定点の所定の範囲の境界を、前記血管壁面から前記血管中心に向かって単調増加から単調減少に転じる点とし、
前記壁面速度勾配演算部は、前記境界を含む所定の範囲内の少なくとも1点以上の測定点における血流速度および血流速度勾配の値から、前記血管壁面における血流速度勾配の推定値を算出することを特徴とする超音波撮像装置。
【請求項５】
請求項４に記載の超音波撮像装置において、
前記壁面速度勾配演算部は、前記血管壁面から前記血管中心に向かって単調増加から単調減少に転じる点および前記上に凸から下に凸に転じる変曲点の少なくともいずれかにおける血流速度および血流速度勾配の値から、血管壁面における血流速度勾配の推定値を算出することを特徴とする超音波撮像装置。
【請求項６】
請求項１に記載の超音波撮像装置において、
前記信号処理部は、前記壁面速度勾配演算部が算出した前記血管壁面における血流速度勾配の推定値を用いて、前記検査対象の診断指標となる情報を算出する診断指標演算部を備えることを特徴とする超音波撮像装置。
【請求項７】
請求項６に記載の超音波撮像装置において、
前記診断指標演算部は、前記血管壁面上の血流速度勾配の推定値から壁面せん断応力を算出する壁面せん断応力演算部であることを特徴とする超音波撮像装置。
【請求項８】
請求項６に記載の超音波撮像装置において、
前記診断指標演算部は、前記血管壁面上の血流速度勾配の推定値から血管弾性を算出する血管弾性演算部であることを特徴とする超音波撮像装置。
【請求項９】
請求項６に記載の超音波撮像装置において、
前記診断指標演算部は、前記血管壁面上の血流速度勾配の推定値から血液流量を算出する血液流量演算部であることを特徴とする超音波撮像装置。
【請求項１０】
請求項１に記載の超音波撮像装置において、
前記信号処理部で得られた情報を表示する表示部を備え、
前記表示部は、前記速度勾配分布演算部が算出した前記血流速度勾配分布を表示し、
前記範囲特定演算部が算出した所定の範囲を、前記血流速度勾配分布に重ねて表示することを特徴とする超音波撮像装置。
【請求項１１】
請求項１に記載の超音波撮像装置において、
前記信号処理部は、前記エコー信号から前記検査対象の断層画像を形成する断層画像形成部を備え、
前記壁面速度勾配演算部は、前記断層画像から特定される血管壁面の位置情報と、前記範囲特定演算部が算出した所定の範囲内の測定点における血流速度および血流速度勾配の値から、血管壁面における血流速度勾配の推定値を算出することを特徴とする超音波撮像装置。
【請求項１２】
請求項６に記載の超音波撮像装置において、
前記信号処理部で得られた情報を表示する表示部を備え、
前記信号処理部は、前記エコー信号から前記検査対象の断層画像を形成する断層画像形成部を備え、
前記表示部は、前記検査対象の診断指標となる情報を前記断層画像に重ねて、空間分布図として表示することを特徴とする超音波撮像装置。
【請求項１３】
請求項６に記載の超音波撮像装置において、
前記検査対象の心拍信号情報を入力する入力部と、
前記入力部および前記信号処理部で得られた情報を表示する表示部と、を備え、
前記表示部は、前記検査対象の診断指標となる情報を前記心拍信号情報とともに時系列変化で表示することを特徴とする超音波撮像装置。
【請求項１４】
超音波撮像装置における演算方法であって、
検査対象が反射するエコー信号から前記検査対象の血管内における超音波照射方向の血流速度成分を抽出する工程と、前記血流速度成分に基づいて、前記血管壁面から前記血管中心へ径方向に並んだ複数の測定点における、血管壁面に沿った方向の血流速度分布を算出する工程と、
微分演算により前記血流速度分布から血流速度勾配分布を算出する工程と、
前記血流速度勾配分布において所定の範囲を算出する工程と、
前記所定の範囲内の測定点から1点以上の測定点を選択する工程と、
選択された測定点における血流速度および血流速度勾配の値から血管壁面における血流速度勾配の推定値を算出する工程と、
を備えることを特徴とする演算方法。
【請求項１５】
請求項１４に記載の演算方法において、
前記選択された測定点における血流速度勾配の値を用いて、前記血管壁面近傍における前記血流速度勾配分布に近似モデルを適用し、前記近似モデルにおける血流速度勾配分布の定積分が前記選択された測定点における血流速度の値と等しくなる関係に基づいて、前記選択された測定点における血流速度および血流速度勾配の値から前記血管壁面における血流速度勾配の推定値を算出することを特徴とする演算方法。
【請求項１６】
請求項１５に記載の演算方法において、
前記血管壁面近傍における前記血流速度勾配分布の近似モデルを1次直線とすることを特徴とする演算方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、医療用の超音波撮像装置に関し、血管壁面近傍における血流速度勾配を求める技術に関する。
【背景技術】
【０００２】
先進国における主要な死因の一つは、心不全などの心血管疾患であり、その多くに動脈硬化が関わる。動脈硬化は、血管壁面に対する血流からの刺激の大小により進行速度が変わることが指摘されている。そのため、血管壁面に作用する壁面せん断応力は、動脈硬化早期の診断指標として注目されている。
【０００３】
壁面せん断応力を算出する方法の一つに、超音波撮像装置を用いて、超音波ドプラ法により測定した血流速度分布に基づいて算出する方法がある。この方法では、血流速度分布の空間微分である血流速度勾配を、血管壁面近傍で高精度に求めることが要求される。
【０００４】
超音波ドプラ法では、ドプラ効果を用いてエコー信号から反射体の移動速度情報が得られるが、反射体には血流中の赤血球だけでなく、血管膜や筋肉などの周辺組織も含まれる。血管壁面近傍における血流は、血管壁面との摩擦により周辺組織の動きと同程度に低速となるため、血流のエコー信号が周辺組織等の検出対象外の信号に埋もれてしまい、血管壁面近傍における血流速度の測定精度は低下する。このことから、超音波撮像装置で測定した血流速度分布を単純に微分することでは、血管壁面近傍で血流速度勾配が高精度に得られず、臨床的に動脈硬化の早期診断が可能な精度は実現されていない。
【０００５】
血流速度勾配の測定精度向上を目的とした超音波撮像方法には、非特許文献１に記載の技術がある。この文献には、「II. MATERIALS AND METHODS」の「B. Velocity reconstruction and WSR measurement」中に、「In the first step the two regions that extend for about 5-10% of the diameter from the wall positions towards the vessel lumen are located. The profile measured in these regions is substituted by a line that starts at the wall with velocity 0 and ramps up to join the remaining -measured- profile.」という記載がある。この方法では、血管壁面から血管径の5-10%離れた領域までの血流速度分布を棄却し、棄却した領域に対して補間を行った後、血管壁面近傍の血流速度勾配を推定している。
【先行技術文献】
【非特許文献】
【０００６】
【非特許文献１】Improved Wall Shear Rate Method for Robust Measurements（Stefano Ricci et al., 2014 IEEE International Ultrasonics Symposium Proceedings, 432-435, 2014.）
【発明の概要】
【発明が解決しようとする課題】
【０００７】
非特許文献１では、血管径を基準として血流の測定値を棄却する領域を決定する方法を適用しており、頸動脈をモデルとしたシミュレーションで検証している。しかし、実測した血流速度分布は、検査対象の血管の形状や柔軟性、血管内皮の性状、拍動などにより様々に異なるため、血管径のみを基準として棄却領域を決定する方法では、血流条件次第で、検出対象外の信号に埋もれた血流速度の測定値が十分に棄却されない場合が生じる。その結果、血管壁面近傍の血流速度勾配の推定に信頼性の低い血流速度の測定値を用いる可能性がある。
【０００８】
本発明は、血管の形状、柔軟性、血管内皮の性状、拍動などが異なる様々な血流条件下においても、血管壁面近傍の血流速度勾配を高精度に推定する超音波撮像装置およびその演算方法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
上記課題を解決するため、本発明においては、検査対象が反射するエコー信号を受信する受信部と、前記受信部により受信された前記エコー信号を処理する信号処理部と、を備え、前記信号処理部は、前記検査対象の血管壁面から血管中心へ径方向に並んだ複数の測定点における、前記エコー信号から算出される前記血管壁面に沿った方向の血流速度の値から、血流速度勾配分布を算出する速度勾配分布演算部と、前記血流速度勾配分布において所定の範囲を算出する範囲特定演算部と、前記所定の範囲内の測定点における血流速度および血流速度勾配の値から、血管壁面における血流速度勾配の推定値を算出する壁面速度勾配演算部と、を備えることを特徴とする超音波撮像装置を提供する。
【００１０】
また、上記課題を解決するため、本発明においては、超音波撮像装置における演算方法であって、検査対象が反射するエコー信号から前記検査対象の血管内における超音波照射方向の血流速度成分を抽出する工程と、前記血流速度成分に基づいて、前記血管壁面から前記血管中心へ径方向に並んだ複数の測定点における、血管壁面に沿った方向の血流速度分布を算出する工程と、微分演算により前記血流速度分布から血流速度勾配分布を算出する工程と、前記血流速度勾配分布において所定の範囲を算出する工程と、前記所定の範囲内の測定点から1点以上の測定点を選択する工程と、選択された測定点における血流速度および血流速度勾配の値から前記血管壁面における血流速度勾配の推定値を算出する工程と、を備えることを特徴とする演算方法を提供する。
【発明の効果】
【００１１】
本発明によれば、様々な血流条件下において、血管壁面近傍の血流速度勾配を高精度に推定することができる。
【図面の簡単な説明】
【００１２】
【図１】実施例１における超音波撮像装置の一構成例を示すブロック図。
【図２】実施例１における信号処理部の動作の実施形態を示す計算処理フローを示す図。
【図３】血流速度分布の例および血流速度勾配分布の例を示す図。
【図４】壁面速度勾配演算に用いる血流速度勾配分布の近似モデルと血管壁面の位置を示す図。
【図５】表示画像形成部において形成する所定の範囲の表示画像の一例を示す図。
【図６】表示画像形成部において形成する壁面せん断応力の空間分布の表示画像の一例を示す図。
【図７】表示画像形成部において形成する壁面せん断応力の時系列変化の表示画像の一例を示す図。
【図８】実施例２における信号処理部の動作の実施形態を示す計算処理フローを示す図。
【図９】実施例３における信号処理部の動作の実施形態を示す計算処理フローを示す図。
【図１０】実施例３における血液流量演算の積分範囲を示す図。
【発明を実施するための形態】
【００１３】
以下、本発明の実施例を図面に従い説明する。
【実施例１】
【００１４】
図１は、実施例１の壁面せん断応力測定機能を有する超音波撮像装置の構成例を示すブロック図である。壁面せん断応力は、血管内皮細胞に変化をもたらす血流刺激で、動脈硬化の早期診断における診断指標として臨床研究で注目されているパラメータである。また、壁面せん断応力の大きさは、進行した動脈硬化で発生するプラークの破綻に影響しており、プラーク破綻のリスク予測とも関係が深い。
【００１５】
装置本体１は超音波探触子２を制御しながら超音波画像を生成するものであり、入力部１０、制御部１１、超音波信号を発信する発信部１２、エコー信号を受信する受信部１３、表示部１４、信号処理部１５、およびメモリ１６を備えている。
【００１６】
超音波探触子２は、被検者の生体３に接し、発信部１２で生成された信号に従い、生体３内の血管３０に対し超音波を照射し、受信部１３は血管３０のエコー信号を受信する。超音波探触子２は、スキャン方式に応じて連続波あるいはパルス波を発生する。また超音波探触子２の走査方法により、２次元的な断面を撮像する平面的撮像法、或いは３次元的な領域を撮像する立体的撮像法を適宜選択してよい。
【００１７】
装置本体１の各構成要素の機能を説明する。入力部１０は、超音波撮像装置を操作する医師や技師（以下、まとめて検者という）が制御部１１に対し超音波撮像装置の動作条件を設定するキーボードやポインティングデバイスを備える。また検査に心電図等の外部機器からの情報を利用する場合、外部機器からの情報を取り込む機能も備える。
【００１８】
制御部１１は、入力部１０によって設定された超音波撮像装置の動作条件に基づき発信部１２、受信部１３、表示部１４および信号処理部１５を制御するもので、例えばコンピュータシステムのＣＰＵ（Central Processing Unit）に構築することができる。
【００１９】
発信部１２は、所定の周波数の信号を発生する発振器を備え、超音波探触子２に駆動信号を送る。受信部１３は、図示していないが、受信回路やサンプリング周波数が通常１０ＭＨｚから５０ＭＨｚのＡ/Ｄ（Analog-to-Digital ）コンバーターを含み、そのほかに、超音波探触子２によって受信されたエコー信号に対し整相加算、検波、増幅などの信号処理を行う。この処理には、一般的な超音波撮像装置が備えている低速の速度成分を排除するフィルタ（以下、ウォールフィルタという）を含む。ただし、Ａ／Ｄコンバーターは受信部１３の代わりに信号処理部１５の前段に備えてもよく、その場合は整相加算、検波、増幅、ウォールフィルタなどの信号処理を信号処理部１５が行う。また、受信部１３は、図示していないが、超音波探触子２の受信素子毎、あるいは素子を束ねた開口部ごとのエコー信号を一時的に保存する受信データメモリを有してもよい。
【００２０】
次に、信号処理部１５の詳細な構成要素を説明する。ここでは、ＣＰＵが実行するソフトウェアにより実現する場合を例に説明する。信号処理部１５は、主要な要素として、断層画像形成部１５１、ドプラ速度抽出部１５２、速度分布演算部１５８、速度勾配分布演算部１５３、範囲特定演算部１５４、壁面速度勾配演算部１５５、壁面せん断応力演算部１５６、表示画像形成部１５７をプログラムとして有し、ＣＰＵがプログラムを読み込んで実行することにより、後述する機能を実現する。
【００２１】
断層画像形成部１５１は、受信部１３から出力されるエコー信号から、検査対象の照射領域の２次元的な組織の断層画像、あるいは２次元アレイ探触子やメカニカル探触子を用いることにより、３次元的な組織の断層画像を形成する。ドプラ速度抽出部１５２は、受信部１３から出力されるエコー信号から超音波照射方向の血流速度成分（以下、ドプラ速度という）を抽出する。その際、ドプラ速度は、平面的撮像法を用いることで２次元的な空間分布として、あるいは立体的撮像法を用いることで３次元的な空間分布として得る。
【００２２】
速度分布演算部１５８は、ドプラ速度抽出部１５２の出力したドプラ速度に対し、血管３０の壁面に沿った方向の血流速度（以下、特に方向を示さない場合、血流速度は血管３０の壁面に沿った方向の血流速度を指す）を算出し、血管３０の壁面から血管中心へ径方向に並んだ複数の測定点における血流速度分布を算出する。
【００２３】
速度勾配分布演算部１５３は、速度分布演算部１５８で算出された血流速度分布から微分演算により血流速度勾配分布を算出する。範囲特定演算部１５４は、速度勾配分布演算部１５３の出力した血流速度勾配分布に対し、所定の範囲を算出する。
【００２４】
壁面速度勾配演算部１５５は、範囲特定演算部１５４が算出した所定の範囲内の測定点における血流速度分布および血流速度勾配分布の値を用いて、血管壁面における血流速度勾配（以下、壁面速度勾配という）を算出する。
【００２５】
壁面せん断応力演算部１５６は、検査対象の診断指標となる情報を算出する診断指標演算部を構成する。壁面せん断応力演算部１５６は、壁面速度勾配演算部１５５が算出した壁面速度勾配の値に対し、例えば予め決められた値や、入力部１０から入力される値、もしくはエコー信号から算出される値として与えられる血液の粘性係数を乗じることで血管壁面に作用する壁面せん断応力を算出する。表示画像形成部１５７は、スキャンコンバーターを備え、表示部１４に表示される表示画像を形成する。形成される表示画像としては、断層画像形成部１５１で形成した断層画像や、ドプラ速度抽出部１５２で抽出したドプラ速度、速度分布演算部１５８で算出した血流速度分布、速度勾配分布演算部１５３で算出した血流速度勾配分布、範囲特定演算部１５４で算出した所定の範囲、壁面速度勾配演算部１５５で算出した壁面速度勾配、壁面せん断応力演算部１５６で算出した壁面せん断応力などの情報である。
【００２６】
なお、信号処理部１５の構成要素の一部又は全部の機能は、制御部１１を構成するものと同一のＣＰＵ、あるいは異なるＣＰＵで実行するソフトウェアで実現する他、ＡＳＩＣ （Application Specific Integrated Circuit ）やＦＰＧＡ（Field-Programmable Gate Array）、ＧＰＵ（Graphics Processing Unit）等のハードウェアで実現してもよい。
【００２７】
メモリ１６は、エコー信号、信号処理部１５での演算に必要な情報（検者が入力部１０により指示した情報等）や信号処理部１５の処理結果（断層画像、ドプラ速度、血流速度分布、血流速度勾配分布、血流速度勾配分布における所定の範囲、壁面速度勾配、壁面せん断応力、表示画像等）を記憶する。
【００２８】
以上説明した装置の構成を踏まえ、信号処理部１５の動作の実施形態の一例を、図２に示す計算処理フローを参照して説明する。
〈ステップＳ１〉
ドプラ速度抽出部１５２は、受信部１３が出力したエコー信号を受け取った後、血流情報を抽出するためウォールフィルタを適用し、カラードプラ法を用いてエコー信号からドプラ速度を抽出する。その際、パルス波ドプラ法などの血流速度測定方式を用いてもよい。また、抽出したドプラ速度の空間分布や時間変化に対して、例えばローパスフィルタや、バンドパスフィルタ、平滑化フィルタなどの様々なフィルタを適用し所望の分布や時間変化を得てもよい。
〈ステップＳ２〉
速度分布演算部１５８は、ドプラ速度成分抽出部により抽出されたドプラ速度を受け取った後、血管３０の壁面に沿った方向の血流速度を算出する。その際、まず、二方向以上の角度で求めたドプラ速度から速度ベクトルを算出するベクトルドプラ法を用いて血流の速度ベクトルを算出し、続いて、速度ベクトルの血管３０の壁面に沿った方向の速度成分を求める。なお、血管壁面に沿った方向の血流速度を求める方法は、ベクトルドプラ法に限らず、例えば、抽出されたドプラ速度に、超音波照射角度または血流の流れ方向を考慮することで推定してもよい。その後、算出された血流速度から、血管３０の壁面から血管中心へ径方向に並んだ複数の測定点における血流速度分布を算出する。その際、血管３０の径方向は、断層画像形成部１５１により形成されメモリ１６に保存されたＢモード画像などの断層画像より自動的に計算されてもよいし、検者が入力部１０により指示してもよい。
〈ステップＳ３〉
速度勾配分布演算部１５３は、速度分布演算部１５８により算出された血流速度分布を受け取った後、微分演算により血流速度勾配分布を算出する。その際、算出した血流速度勾配分布に対して、例えばローパスフィルタや、バンドパスフィルタ、平滑化フィルタなどの様々なフィルタを適用し所望の分布を得てもよい。
〈ステップＳ４〉
範囲特定演算部１５４は、速度勾配分布演算部１５３により算出された血流速度勾配分布を受け取った後、まず血流速度勾配分布の極値および変曲点の少なくともいずれかを算出する。この血流速度勾配分布の極値および変曲点の少なくともいずれかの座標を用いて血流速度勾配分布の棄却する領域を決定する根拠について、図３を用いて説明する。
【００２９】
図３（ａ）は血管壁面から血管中心へ径方向（ｙ方向）に並んだ複数の測定点における血流速度の測定値を分布図として表したものであり、図３（ｂ）は図３（ａ）の血流速度分布を微分演算することにより算出される血流速度勾配分布を表したものである。ただし、図中の点線は血管壁面の位置、一点鎖線は血管中心をそれぞれ表す。ここで、図３（ａ）、（ｂ）の分布図は、ウォールフィルタの影響を受けている。
【００３０】
一方、図３（ｃ）は血流速度分布の近似モデルの例、図３（ｄ）は血流速度勾配分布の近似モデルの例である。まず、これらの近似モデルについて説明する。血管を最も簡易的に表す３次元形状は円断面の直管（以下、単に円管という）である。仮に、血流が時間依存性の全く無い定常流れであるとすると、血流速度分布は、数式（１）に示すHagen-Poiseuille流れの式より、２次曲線で表される。
【００３１】
【数１】
2018011634.tif 000003
【００３２】
ただし、uは血流速度、Rは血管半径、ｘは流れ方向の距離、ｒは血管中心からの半径方向の距離を表す。また、ｄｐは流れ方向の２箇所の微小区間ｄｘにおける圧力ｐの変化量であり、μは血液の粘性係数である。
【００３３】
血流速度勾配分布は血流速度分布の空間微分として算出できることから、円管の血流速度勾配分布は、数式（２）に示す式より、１次直線で表される。
【００３４】
【数２】
2018011634.tif 000004
【００３５】
実際の血流では、湾曲などの複雑な血管形状や、拍動などの血流の時間変動が存在するため、血流速度分布はHagen-Poiseuille流れと同じにはならないが、血管壁面近傍の局所では、壁面の滑りなし条件の影響が強く、Hagen-Poiseuille流れと共通する。このため、図３（ｃ）、（ｄ）では、血管壁面近傍の血流速度分布を２次曲線、血管壁面近傍の血流速度勾配分布を１次直線で近似し、血管中心付近の血流速度勾配分布は下に凸で、緩やかに0に収束するものとしてモデル化している。ただし、血流速度分布の近似モデルは２次曲線に限られるものではなく、多項式や指数関数などその他の関数で近似することもできる。
【００３６】
続いて、図３（ａ）および（ｂ）に示すウォールフィルタの影響を受けた血流速度分布および血流速度勾配分布について説明する。ウォールフィルタは、血流と比較して低速の、血管膜や筋肉など周囲組織の動きを除外するためのフィルタである。そのため、特性上、血管壁面から受ける摩擦力により低速となった血管壁面近傍の血流速度の信号を弱め、結果として血流速度が過小に測定される。したがって、ウォールフィルタ適用後の血流速度分布を微分演算することで得られた血流速度勾配分布は、図３（ｂ）に示すように、血管中心から血管壁面へ向かうに従って単調増加から単調減少に転じ、血管壁面における血流速度勾配が０に収束する傾向を示す。
【００３７】
以上に述べた理由から、図３（ｂ）に示す血流速度勾配分布において、少なくとも血管壁面から血管中心に向かって単調増加から単調減少に転じる点（以下、極値aという）から血管壁面までの範囲は、血流速度および血流速度勾配が過小評価されている。また、血流速度勾配分布の変曲点のうち、上に凸から下に凸に転じる点（以下、変曲点ｂという）は、血流速度勾配分布が図３（ｄ）の近似モデルに合う形状となったことを示す指標となる。したがって、血流速度勾配分布の極値aや変曲点ｂは、様々な血流条件下において、ウォールフィルタの影響で血流速度が過小評価された領域、すなわち、周辺組織等の検出対象外の信号に血流速度の測定値が埋もれた領域の指標となる。
〈ステップＳ５〉
範囲特定演算部１５４は、続いて、極値aおよび変曲点ｂの少なくともいずれかを基準として血流速度勾配分布の棄却する領域を決定し、前記血流速度勾配分布から棄却領域を除いた範囲、すなわち血管３０の壁面における血流速度勾配の推定に分布上の値を用いる所定の範囲を算出する。極値aを所定の範囲の境界とした場合、壁面に近い測定点の分布上の値を壁面速度勾配の推定に用いることができる。また、変曲点ｂから血管壁面までをウォールフィルタの影響で過小評価された領域とみなし、変曲点ｂを所定の範囲の境界としてもよい。また、これらの座標を基準として検者により設定される、これらの座標と異なる点を所定の範囲の境界としてもよい。算出される範囲は、着目している血管壁面から血管中心までにおいて定義されてもよいし、血管壁面から血管中心を通って対向する血管壁面までにおいて定義されてもよい。
〈ステップＳ６〉
壁面速度勾配演算部１５５は、速度分布演算部１５８から算出された血流速度分布、速度勾配分布演算部１５３から算出された血流速度勾配分布、範囲特定演算部１５４から算出された所定の範囲を受け取った後、まず所定の範囲内から、壁面速度勾配の演算に用いる測定点を選ぶ。
【００３８】
その際、後述する壁面速度勾配演算の方式によって、選択する測定点は一点でも複数点でもよい。好ましくは、下記の理由より、極値aおよび変曲点ｂの二点の測定点を選択するとよい。極値aを選択することで、所定の範囲内で最も壁面に近い測定点における血流速度および血流速度勾配の値を用いることができ、さらに、変曲点ｂを選択することで、血流速度勾配分布が図３（ｄ）の近似モデルに合う形状となった点における血流速度および血流速度勾配の値を用いることができる。
〈ステップＳ７〉
壁面速度勾配演算部１５５は、続いて、選択された測定点における血流速度および血流速度勾配の値を用いて壁面速度勾配の推定値を算出する。
【００３９】
ここでは、代表的な例として、選択された測定点が一点の場合（ａ）と二点の場合（ｂ）について図４を用いて説明する。図４では、極値aを境界の一方とし、もう一方の境界を血管中心として、極値aから血管中心までの範囲を所定の範囲とする。ただし、所定の範囲の内側は境界を含むものとする。ここでは、図３（ｄ）と同様に、選択された測定点から血管壁面までの血流速度勾配分布を一次直線で近似した場合について説明する。
【００４０】
測定点が一点の場合が図４（ａ）であり、選択された測定点から血管壁面までの血流速度勾配分布を１次直線４１で近似する場合、図１および図２には図示していないが、Ｂモード画像など断層画像形成部１５１で形成した断層画像から血管壁面の位置ｙｗを特定する必要がある。その際、血管壁面の位置ｙｗの特定方法としては、断層画像形成部１５１により形成されメモリ１６に保存されたＢモード画像などの断層画像より自動的に計算されてもよいし、検者が入力部１０により指示してもよい。
血流速度勾配の定積分が血流速度となることから、測定点ｙ０における血流速度は台形ＡＢＣＤの面積となり、数式（３）が成立する。
【００４１】
【数３】
2018011634.tif 000005
【００４２】
ただし、ｕ０は測定点ｙ０における血流速度、ｄｕ／ｄｙは血流速度勾配を表す。
【００４３】
したがって、壁面速度勾配は数式（４）で算出される。
【００４４】
【数４】
2018011634.tif 000006
【００４５】
すなわち、血管壁面近傍の血流速度勾配分布に近似モデル（１次直線）を適用し、血管壁面から選択された測定点までの血流速度勾配分布の定積分が選択された測定点における血流速度と等しくなる関係に基づき、選択された測定点における血流速度および血流速度勾配の値から壁面速度勾配の推定値を算出することができる。
【００４６】
また、測定点が二点の場合が図４（ｂ）であり、各々の測定点から血管壁面までの血流速度勾配分布を１次直線４２で近似する。その際、血管壁面の位置は必ずしも特定しておく必要はない。血管壁面の座標ｙｗを未知数として、測定点ｙ1およびｙ2における血流速度と血流速度勾配の関係は、台形ＡＢＣＤおよび台形ＡＢＥＦより数式（５）で記述される。
【００４７】
【数５】
2018011634.tif 000007
【００４８】
数式（５）を連立方程式として解くと、壁面速度勾配(ｄｕ／ｄｙ)ｙ=ｙｗは算出される。例えば、測定点ｙ１を極値aのｙ座標、測定点ｙ２を変曲点ｂのｙ座標として、数式（５）に代入することで、血管壁面の位置ｙｗが不明であっても壁面速度勾配を求めることが可能である。
〈ステップＳ８〉
壁面せん断応力演算部１５６は、壁面速度勾配演算部１５５から壁面速度勾配の値を受け取った後、壁面せん断応力を算出する。壁面せん断応力τは、数式（６）によって与えられる。
【００４９】
【数６】
2018011634.tif 000008
【００５０】
ここで、μは血液の粘性係数である。
〈ステップＳ９〉
表示画像形成部１５７は、断層画像形成部から血管３０の断層画像、ドプラ速度抽出部１５２からドプラ速度、速度分布演算部１５８から血管３０の壁面に沿った方向の血流速度分布、速度勾配分布演算部１５３から血流速度勾配分布、範囲特定演算部１５４から所定の範囲、壁面速度勾配演算部１５５から壁面速度勾配の値、壁面せん断応力演算部１５６から壁面せん断応力の情報を受け取った後、それらの情報の全部または一部を、予め決められた形式や入力部１０から入力される指示に従い、表示画像として形成する。
【００５１】
本実施例の装置構成においては、表示部１４は速度分布演算部１５８が算出した血流速度分布や速度勾配分布演算部１５３が算出した血流速度勾配分布を表示し、範囲特定演算部１５４が算出した所定の範囲を、血流速度分布や血流速度勾配分布に重ねて表示してもよい。例えば、図５（ａ）または（ｂ）に示すように、血管壁面に沿った方向の血流速度分布や血流速度勾配分布に対し、所定の範囲５１、５２の範囲内または範囲外を色付けする、あるいは単に所定の範囲５１、５２の境界の座標を数値として表示する表示画像を形成する。
【００５２】
また、図６（ａ）または（ｂ）に示すように、複数点について算出した壁面速度勾配または壁面せん断応力の空間分布情報を知覚情報として、血管の断層画像６１、６２と共に検者に提示する。すなわち、断層画像形成部１５１が形成した検査対象の血管の断層画像と、壁面速度勾配または壁面せん断応力の空間分布情報などの検査対象の診断指標となる情報を重ねて、空間分布図として表示する。
【００５３】
その際、知覚情報を、図６（ａ）に示すように壁面せん断応力の大きさを幅で表示してもよいし、図６（ｂ）に示すようにカラーバー６３を利用して壁面せん断応力の大きさを色分けして表示してもよい。このように、Ｂモード画像などの血管の断層画像に、血流速度分布および血流速度勾配分布の少なくとも一つに重ねて表示することは、これらの情報の空間分布に対する検者の理解を助ける。
【００５４】
また、図７に示すように、壁面速度勾配または壁面せん断応力の時系列変化情報７１を、検者に提示してもよい。同図の横軸は時間（ｓ）、縦軸は壁面せん断応力（Ｐａ）である。特に、時系列変化の情報を心拍信号とともに表示することは、壁面せん断応力の時系列変化と拍動の相関に対する検者の理解を助ける。すなわち、入力部１０から検査対象の心拍信号情報を入力するよう指示し、表示部１４は、検査対象の診断指標となる壁面せん断応力などの情報を入力された心拍信号情報とともに時系列変化で表示する。さらに、壁面せん断応力の空間分布情報や時系列変化情報の統計値、例えば最大値、最小値、平均値、中央値などの少なくとも一つを、検者に提示してもよい。
【実施例２】
【００５５】
実施例１においては、検査対象の診断指標となる情報を算出する診断指標演算部として、壁面速度勾配の推定値を用いて壁面せん断応力を算出する壁面せん断応力演算部を備える構成を説明したが、壁面速度勾配は、圧較差による血管弾性測定法の計算にも用いられる。実施例２として、検査対象の診断指標となる情報を算出する診断指標演算部として血管弾性演算部を有する超音波撮像装置の実施例を、実施例１の超音波撮像装置と異なる点を中心に説明する。計算処理フローではステップＳ８以降が異なるため、図２のステップＳ８、Ｓ９の代わりにステップＳ１０、Ｓ１１を挿入したフローチャートを図８に示す。
〈ステップＳ１０〉
実施例１では、壁面速度勾配演算部１５５の出力を壁面せん断応力演算部１５６が受け取り、壁面せん断応力を算出した。実施例２では、図示を省略するが、壁面せん断応力演算部１５６に代え血管弾性演算部を備え、血管弾性演算部は、図１の壁面速度勾配演算部１５５の出力を受け取り、血管弾性を算出する。また、算出された血管弾性はメモリ１６に保存される。血管弾性Ｅは数式（７）によって与えられる。
【００５６】
【数７】
2018011634.tif 000009
【００５７】
ただし、Ｒは血管の内径、ｈは血管の厚さ、Ｌは圧較差算出に用いる２点間の距離、Δｄは前記２点間の内径差、μは血液の粘性係数を示す。これらのパラメータは、例えば予め決められた値や、入力部１０から入力される値、もしくはエコー信号から算出される値のいずれを採用してもよい。
〈ステップＳ１１〉
実施例２では、表示画像形成部１５７は、断層画像形成部１５１から血管３０の断層画像、ドプラ速度抽出部１５２からドプラ速度、速度分布演算部１５８から血管３０の壁面に沿った方向の血流速度分布、速度勾配分布演算部１５３から血流速度勾配分布、範囲特定演算部１５４から所定の範囲、壁面速度勾配演算部１５５から壁面速度勾配の値、血管弾性演算部から血管弾性の情報を受け取った後、それらの情報の全部または一部を、予め決められた形式や入力部１０から入力される指示に従い、表示画像として形成する。表示形式の具体例は、図２のステップＳ９に準じる。
【実施例３】
【００５８】
実施例３は、検査対象の診断指標となる情報を算出する診断指標演算部として血液流量演算部を有する超音波撮像装置の実施例である。実施例３の構成を、実施例１の超音波撮像装置と異なる点を中心に説明する。本実施例では、図１に示した装置構成の壁面せん断応力演算部１５６に代え流量演算部を備える。流量演算部は壁面速度勾配演算部１５５の出力を受け取り、血液流量を算出する。具体的には、壁面速度勾配演算部１５５が算出した壁面速度勾配の値を用いて、流量演算部が血管壁面から測定点ｙ１までの範囲の血流速度分布を前記の近似モデル（２次曲線）で置換し、置換後の流速分布を積分することにより血液流量を算出する。
【００５９】
本実施例の構成によれば、より高精度に血液流量を求めることができる。計算処理フローでは実施例１と内容が異なるステップＳ８以降に関して、図２のステップＳ８、Ｓ９の代わりにステップＳ１２、Ｓ１３、Ｓ１４、Ｓ１５を挿入したフローチャートを図９に示す。
〈ステップＳ１２〉
図示を省略した流量演算部は、壁面速度勾配の数値を受け取った後、血液流量を算出する。まず、壁面速度勾配演算部１５５が算出した壁面速度勾配の値を用いると、血流速度勾配分布の近似モデルは数式（８）によって与えられる。
【００６０】
【数８】
2018011634.tif 000010
【００６１】
〈ステップＳ１３〉
数式（８）を積分することにより、血流速度分布の近似モデルは数式（９）によって与えられる。
【００６２】
【数９】
2018011634.tif 000011
【００６３】
〈ステップＳ１４〉
数式（９）を血管断面全体で積分することにより、血液流量が求められる。図１０は実施例３の構成で血液流量演算の積分範囲を説明するための図であり、１０１、１０３はモデル化した血流速度分布を積分する領域を、１０２は実測した血流速度分布を積分する領域を示す。同図に示すように、モデル化した血流速度分布を積分する領域１０１における血管壁面Ａ近傍の血流速度分布の近似モデルをｕＡ、モデル化した血流速度分布を積分する領域１０３における血管壁面Ｂ近傍の血流速度分布の近似モデルをｕＢ、実測した血流速度分布を積分する領域１０２における血管中央の実測した血流速度分布をｕＣとすると、血液流量は次式で求められる。
【００６４】
【数１０】
2018011634.tif 000012
【００６５】
ただし、ｙｗＡは血管壁面Ａの座標、ｙｗＢは血管壁面Ｂの座標、ｙＡは血管壁面Ａ近傍の測定点、ｙＢは血管壁面Ｂ近傍の測定点をそれぞれ示す。
〈ステップＳ１５〉
実施例３では、表示画像形成部１５７は、断層画像形成部１５１から血管３０の断層画像、ドプラ速度抽出部１５２からドプラ速度、速度分布演算部１５８から血管３０の壁面に沿った方向の血流速度分布、速度勾配分布演算部１５３から血流速度勾配分布、範囲特定演算部１５４から所定の範囲、壁面速度勾配演算部１５５から壁面速度勾配の値、血液流量演算部から血液流量の情報を受け取った後、それらの情報の全部または一部を、予め決められた形式や入力部１０から入力される指示に従い、表示画像として形成する。表示形式の具体例は、ステップＳ９に準じる。
【００６６】
本発明の超音波撮像装置は、上記の実施例に限定されず、適宜要素の追加、削除などが可能である。例えば、各実施例は検査対象の診断指標となる情報を算出する診断指標演算部として、それぞれ壁面せん断応力演算部、血管弾性演算部、あるいは血液流量演算部を備えていたが、これらの二つ、あるいは三つを組み合わせて備える構成とすることも可能である。また、実施例１は動脈硬化の早期診断を例に説明したが、これは本発明の適用対象を動脈に限定するものではなく、下肢静脈などの静脈における血流速度勾配測定にも適用することが可能である。下肢静脈は、血栓や静脈瘤の生じやすい部位であり、これらの診断にも本発明を適用してもよい。
【符号の説明】
【００６７】
１  装置本体
２  超音波探触子
３  生体
１０  入力部
１１  制御部
１２  発信部
１３  受信部
１４  表示部
１５  信号処理部
１６  メモリ
３０  血管
４１、４２  一次直線
５１、５２  所定の範囲
６１、６２  血管の断層画像
６３  カラーバー
７１  時系列変化情報
１０１、１０３  モデル化した血流速度分布を積分する領域
１０２  実測した血流速度分布を積分する領域
１５１  断層画像形成部
１５２  ドプラ速度抽出部
１５３  速度勾配分布演算部
１５４  範囲特定演算部
１５５  壁面速度勾配演算部
１５６  壁面せん断応力演算部
１５７  表示画像形成部
１５８  速度分布演算部
【図１】
2018011634.tif 000013
【図２】
2018011634.tif 000014
【図３】
2018011634.tif 000015
【図４】
2018011634.tif 000016
【図５】
2018011634.tif 000017
【図６】
2018011634.tif 000018
【図７】
2018011634.tif 000019
【図８】
2018011634.tif 000020
【図９】
2018011634.tif 000021
【図１０】
2018011634.tif 000022