(57)【要約】【課題】より高精度な血管機能評価を可能とする超音波診断装置を提供する。【解決手段】エコートラッキング処理部２２は、受信信号に基づいて血管壁の位置をトラッキングする。血管径変化波形生成部２４は、トラッキング処理部２２により特定される血管１６の前壁１８および後壁１９の位置に基づいて血管径を測定し、時々刻々と変化する血管径の値を表す血管径変化波形を生成する。血管径特性曲線生成部２８は、血管径変化波形に対して心拍変動成分を除去する処理を行うことで血管径変化の特性曲線を生成する。【選択図】        図１
2004290408.tif 000002
【特許請求の範囲】
【請求項１】
血管内の血液の流れを制御することによって生じた時間的変化状態を計測する超音波診断装置であって、
血管に対して超音波の送受波を行うことで受信信号を出力する送受波手段と、
前記時間的変化状態において得られた前記受信信号に基づいて、血管の形態および血液の流れの少なくとも一方についての時間的変化を表す変化波形を測定する変化波形測定手段と、
前記変化波形に対して局所変動成分の除去処理を行うことで、前記変化波形の大局的な時間変化に対応した診断特性を生成する診断特性生成手段と、
を有することを特徴とする超音波診断装置。
【請求項２】
請求項１に記載の超音波診断装置であって、
前記局所変動成分は心拍変動成分であることを特徴とする超音波診断装置。
【請求項３】
請求項２に記載の超音波診断装置であって、
前記局所変動成分の除去処理は、前記変化波形において一心拍期間内ごとの代表値を求め、これらの代表値を利用した前記変化波形の曲線近似処理であることを特徴とする超音波診断装置。
【請求項４】
請求項３に記載の超音波診断装置であって、
前記診断特性生成手段は、心電図波形を利用して前記一心拍期間を特定することを特徴とする超音波診断装置。
【請求項５】
請求項４に記載の超音波診断装置であって、
前記代表値は、前記変化波形の一心拍期間内における最大値、最小値または平均値であることを特徴とする超音波診断装置。
【請求項６】
請求項１から５のいずれか１項に記載の超音波診断装置であって、
前記診断特性を表示する表示手段をさらに有することを特徴とする超音波診断装置。
【請求項７】
血管内の血液の流れを制御することによって生じた血管径の時間的変化状態を計測する超音波診断装置であって、
血管に対して超音波の送受波を行うことで受信信号を出力する送受波手段と、
前記受信信号に基づいて前記血管の壁部の位置をトラッキングするエコートラッキング処理手段と、
前記血管の壁部の位置に基づいて前記血管の血管径を計測し、前記時間的変化状態における前記血管径の時間的変化を表す血管径変化波形を生成する血管径変化波形生成手段と、
前記血管径変化波形に対して心拍変動成分の除去処理を行うことで、血管径変化の特性曲線を生成する血管径特性曲線生成手段と、
を有することを特徴とする超音波診断装置。
【請求項８】
請求項７に記載の超音波診断装置であって、
前記特性曲線に基づいて、血管径の変化率を演算する血管径変化率演算手段をさらに有することを特徴とする超音波診断装置。
【請求項９】
血管内の血液の流れを制御することによって形成された血流速度の時間的変化状態を計測する超音波診断装置であって、
血管に対して超音波の送受波を行うことで受信信号を出力する送受波手段と、
前記受信信号に基づいて前記血管の壁部の位置をトラッキングするエコートラッキング処理手段と、
前記血管の壁部の位置に基づいて血流速度計測用のサンプルゲートを設定し、
前記受信信号に基づいて前記サンプルゲート内の血流速度を計測する血流速度計測手段と、
前記時間的変化状態において得られた前記血流速度に基づいて、血流速度の時間的変化を表す血流速度変化波形を生成する血流速度変化波形生成手段と、
前記血流速度変化波形に対して心拍変動成分の除去処理を行うことで、血流速度変化の特性曲線を生成する血流速度特性曲線生成手段と、
を有することを特徴とする超音波診断装置。
【請求項１０】
請求項９に記載の超音波診断装置であって、
前記特性曲線に基づいて、血流速度の変化率を演算する血流速度変化率演算手段をさらに有することを特徴とする超音波診断装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は超音波診断装置に関し、特に血管機能評価を行う超音波診断装置に関する。
【０００２】
【従来の技術】
血流反応性血管拡張（以下、ＦＭＤ：Ｆｌｏｗ  Ｍｅｄｉａｔｅｄ  Ｄｉｌａｔａｔｉｏｎ）による血管機能計測は、血管の内皮機能の評価や動脈硬化の評価などに有用とされている。その計測手順は次のとおりである。まず被検査者の安静時における上腕動脈の血管径を計測し、次に被検査者の前腕部分をカフなどで５分程度駆血する。一般に駆血状態において上腕動脈は安静時よりもその血管径が小さくなる。その後駆血を解除すると上腕動脈の血管径が拡張し、徐々に安静時の血管径に戻る。そして、駆血解除から安静時状態に戻るまでの血管径の変化の様子から上腕動脈の内皮機能が評価される。
【０００３】
ＦＭＤによる血管計測において、血管径の計測には超音波診断装置が利用されている。つまり、超音波診断装置により安静時、駆血時、駆血解除後の上腕動脈のＢモード断層画像をＶＴＲに録画しておき、録画されたＢモード断層画像から、所定時刻における血管径を手動で計測していた。このため、血管径の変化情報は断片的であり、また、血管径の計測に多大な時間を要するにもかかわらず、血管径の計測精度も低いものであった。
【０００４】
一方、超音波診断装置には、血管壁の位置を追従検出することで、血管径の変化波形を取得するものがある（例えば、特許文献１参照）。この超音波装置を利用することにより、血管径の変化をリアルタイムで、しかも高精度に計測することが可能になる。
【０００５】
【特許文献１】
特開２００１−２９９７５２号公報
【０００６】
【発明が解決しようとする課題】
上述の超音波診断装置を、ＦＭＤによる計測に利用することで、その計測精度を高めることが期待できる。ところが、血管径は脈動の影響で一心拍内においても局所的に変化している。したがって、ＦＭＤによる測定において、一心拍内の局所的変動を考慮しつつ、駆血解除後の血管径の拡張収縮変化を計測できれば、より高精度な測定を実現することができる。
【０００７】
そこで本発明は、より高精度な血管機能評価を可能とする超音波診断装置を提供することを目的とする。
【０００８】
【課題を解決するための手段】
（１）上記目的を達成するために、本発明に係る超音波診断装置は、血管内の血液の流れを制御することによって形成された時間的変化状態を計測する超音波診断装置であって、血管に対して超音波の送受波を行うことで受信信号を出力する送受波手段と、前記時間的変化状態において得られた前記受信信号に基づいて、血管の形態および血液の流れの少なくとも一方についての時間的変化を表す変化波形を測定する変化波形測定手段と、前記変化波形に対して局所変動成分の除去処理を行うことで、前記変化波形の大局的な時間変化に対応した診断特性を生成する診断特性生成手段とを有するものとする。
【０００９】
上記構成において、血管内の血液の流れの制御とは、例えば、カフなどにより駆血を行った後に駆血を解除することであり、あるいは、血管拡張剤を患者に投与することなどである。また診断特性とは、例えば、変化波形に対して局所変動成分の除去処理を施して形成される波形や、形成された波形から求められる各種評価値などである。上記構成によれば、変化波形に対して局所変動成分の除去処理が施され、変化波形における局所的な変動が実質的にキャンセルされて診断特性が生成されるため、局所的な変動に左右されずに高精度な血管機能評価が可能となる。
【００１０】
望ましくは、前記局所変動成分は心拍変動成分とする。ここで心拍変動成分とは、脈動による局所的な変動成分である。上記構成によれば、変化波形における心拍変動成分が実質的にキャンセルされて診断特性が生成されるため、心拍による局所的な変動に左右されない診断が容易に可能になる。
【００１１】
望ましくは、前記局所変動成分の除去処理は、前記変化波形において一心拍期間内ごとの代表値を求め、これら代表値を利用した前記変化波形の曲線近似処理とする。さらに望ましくは、前記診断特性生成手段は、心電図波形を利用して前記一心拍期間を特定するものとする。さらに望ましくは、前記代表値は、前記変化波形の一心拍期間内における最大値、最小値または平均値とする。さらに望ましくは、前記診断特性を表示する表示手段を有するものとする。
【００１２】
（２）また、上記目的を達成するために、本発明に係る超音波診断装置は、血管内の血液の流れを制御することによって形成された血管径の時間的変化状態を計測する超音波診断装置であって、血管に対して超音波の送受波を行うことで受信信号を出力する送受波手段と、前記受信信号に基づいて前記血管の壁部の位置をトラッキングするエコートラッキング処理手段と、前記血管の壁部の位置に基づいて前記血管の血管径を計測し、前記時間的変化状態における前記血管径の時間的変化を表す血管径変化波形を生成する血管径変化波形生成手段と、前記血管径変化波形に対して心拍変動成分の除去処理を行うことで、血管径変化の特性曲線を生成する血管径特性曲線生成手段とを有するものとする。
【００１３】
上記構成によれば、血管壁部の変位に追従して血管径変化波形が計測され、また、血管径変化波形における心拍変動成分が実質的にキャンセルされて特性曲線が生成されるため、リアルタイムでしかも心拍による局所的な変動に左右されずに血管径の広域的な変化を計測することができる。
【００１４】
望ましくは、前記特性曲線に基づいて、血管径の変化率を演算する血管径変化率演算手段をさらに有するものとする。上記構成によれば、リアルタイムで高精度な血管径の変化率の自動評価が可能になる。
【００１５】
（３）また、上記目的を達成するために、本発明に係る超音波診断装置は、血管内の血液の流れを制御することによって形成された血流速度の時間的変化状態を計測する超音波診断装置であって、血管に対して超音波の送受波を行うことで受信信号を出力する送受波手段と、前記受信信号に基づいて前記血管の壁部の位置をトラッキングするエコートラッキング処理手段と、前記血管の壁部の位置に基づいて血流速度計測用のサンプルゲートを設定し、前記受信信号に基づいて前記サンプルゲート内の血流速度を計測する血流速度計測手段と、前記時間的変化状態において得られた前記血流速度に基づいて、血流速度の時間的変化を表す血流速度変化波形を生成する血流速度変化波形生成手段と、前記血流速度変化波形に対して心拍変動成分の除去処理を行うことで、血流速度変化の特性曲線を生成する血流速度特性曲線生成手段とを有するものとする。
【００１６】
上記構成によれば、血流速度計測用のサンプルゲートを、血管径上の所望の位置に血管壁部の位置に追従して設定することができる。また、血流速度変化波形における心拍変動成分が実質的にキャンセルされて特性曲線が生成される。したがって、リアルタイムでしかも心拍による局所的な変動に左右されずに、所望の位置における血流速度の大局的な時間変化を測定することができる。
【００１７】
望ましくは、前記特性曲線に基づいて、血流速度の変化率を演算する血流速度変化率演算手段をさらに有するものとする。上記構成によれば、リアルタイムで高精度な血流速度の変化率の自動評価が可能になる。
【００１８】
【発明の実施の形態】
以下、本発明の好適な実施の形態を図面に基づいて説明する。
【００１９】
図１には、本発明に係る超音波診断装置の好適な実施形態が示されており、図１はその全体構成を示すブロック図である。送受信部１２は、プローブ１０を介して患者１４内の血管１６に対して超音波の送受波を行うことで受信信号を取得し、断層画像形成部２０およびエコートラッキング処理部２２に出力する。断層画像形成部２０は受信信号に基づいて血管１６の断層画像を形成する。形成された断層画像は表示画像形成部４２において表示処理が施されてディスプレイ４４に表示される。
【００２０】
エコートラッキング処理部２２は、受信信号に基づいて血管壁の位置をトラッキングする。血管壁は比較的強いエコーを生じるため、受信信号のエコーデータレベルを所定の閾値と比較することによって血管壁が検出される。エコートラッキング処理部２２は、測定ライン１７上においてプローブ１０側に位置する前壁１８とプローブ１０から遠方側に位置する後壁１９とを、血管１６の変動に伴う位置変化に追従して検出する。
【００２１】
血管径変化波形生成部２４は、エコートラッキング処理部２２により特定される血管１６の前壁１８および後壁１９の位置に基づいて、前壁１８と後壁１９との間の距離である血管径を測定する。血管径変化波形生成部２４は時々刻々と変化する血管径の値を表す血管径変化波形を生成して表示画像形成部４２およびメモリ（１）２６へ出力する。表示画像形成部４２へ出力された血管径変化波形はディスプレイ４４にリアルタイム表示される。メモリ（１）２６に出力された血管径変化波形は血管径特性曲線生成部２８で利用される。
【００２２】
血管径特性曲線生成部２８は、血管径変化波形に対して心拍変動成分を除去する処理を行うことで血管径変化の特性曲線を生成する。血管径は脈動の影響で一心拍内において変化している。血管機能評価、例えばＦＭＤによる測定では駆血解除から安静時状態に戻るまでの血管径の拡張収縮の様子などを調べている。一般的に駆血解除から安静時状態に戻るまでの時間は数分程度である。つまり、血管径は心拍による局所的な拡張収縮変化を伴いつつ、数分程度に亘る大局的な収縮変化を行っている。このため、ＦＭＤによる測定において駆血解除から安静時状態に戻るまでの血管径の収縮の様子を精度よく調べるためには、心拍による局所的変動成分の影響を考慮しつつ大局的な血管径の収縮変化を計測する必要がある。
【００２３】
血管径特性曲線生成部２８は、心拍変動成分を除去するために血管径変化波形における一心拍期間内ごとの代表値を求め、この代表値を利用して血管径変化の特性曲線を生成する。血管径特性曲線生成部２８は心電図波形を利用して一心拍期間を抽出している。つまり、心臓左室の収縮初期時に発せられる心電図波形中のＲ波を利用し、Ｒ波の発生から次のＲ波の発生までの期間における血管径変化波形の部分波形をメモリ（１）２６から取得して、この期間における代表値を演算する。代表値としては期間内、つまり一心拍内の最大値や最小値、あるいは平均値などであり、代表値の演算は各心拍ごとに実施される。血管径特性曲線生成部２８は、次々に演算される代表値を曲線で結んだ特性曲線を生成し、表示画像形成部４２および評価値演算部４０に出力する。表示画像形成部４２へ出力された特性曲線はディスプレイ４４に表示される。
【００２４】
血流速度計測部３２は、血管１６の壁部の位置に基づいて血流速度計測用のサンプルゲート３３を設定し、受信信号に基づいてサンプルゲート３３内の血流速度を計測する。サンプルゲート３３は、エコートラッキング処理部２２により特定される血管１６の前壁１８および後壁１９の位置に基づいて設定され、例えば、前壁１８および後壁１９の中点を基準に設定される。サンプルゲート３３はユーザにより所望の位置に設定されてもよい。血流速度計測部３２は、サンプルゲート３３内の受信信号を参照することで、その受信信号からドプラ情報を抽出して速度情報を取得し、サンプルゲート３３内における血流速度を演算する。血流速度とは、サンプルゲート３３内における血流の平均速度、最大速度、最小速度などである。計測された血流速度は血流速度変化波形生成部３４に出力される。
【００２５】
血流速度変化波形生成部３４は、血流速度計測部３２から出力される時々刻々と変化する血流速度に基づいて血流速度変化波形を生成して表示画像形成部４２およびメモリ（２）３６へ出力する。表示画像形成部４２へ出力された血流速度変化波形はディスプレイ４４にリアルタイム表示される。メモリ（２）３６に出力された血流速度変化波形は血流速度特性曲線生成部３８で利用される。
【００２６】
血流速度特性曲線生成部３８は、血流速度変化波形に対して心拍変動成分を除去する処理を行うことで、血流速度変化の特性曲線を生成する。血管径と同様に血流速度も脈動の影響で一心拍内において変化している。ＦＭＤによる測定では、駆血解除から安静時状態に戻るまでの血流速度の増加減少の様子などを調べている。一般的に駆血解除から安静時状態に戻るまでの時間は数分程度である。つまり、血流速度は心拍による局所的な増加減少を伴いつつ、数分程度に亘る大局的な増加減少傾向を示す。このため、ＦＭＤによる測定において駆血解除から安静時状態に戻るまでの血流速度の増加減少傾向を精度よく調べるためには、心拍による局所的変動成分の影響を考慮しつつ大局的な増加減少傾向を計測する必要がある。
【００２７】
血流速度特性曲線生成部３８は、心拍変動成分を除去するために血流速度変化波形における一心拍期間内ごとの代表値を求め、この代表値を利用して血流速度変化の特性曲線を生成する。血流速度特性曲線生成部３８は心電図波形を利用して一心拍期間を抽出している。つまり、心臓左室の収縮初期時に発せられる心電図波形中のＲ波を利用し、Ｒ波の発生から次のＲ波の発生までの期間における血流速度変化波形の部分波形をメモリ（２）３６から取得して、この期間における代表値を演算する。代表値としては期間内、つまり一心拍内の最大値、最小値あるいは平均値であり、代表値の演算は各心拍ごとに実施される。血流速度特性曲線生成部３８は、次々に演算される代表値を曲線で結んだ特性曲線を生成し、表示画像形成部４２および評価値演算部４０に出力する。表示画像形成部４２へ出力された特性曲線はディスプレイ４４に表示される。
【００２８】
評価値演算部４０は、血管径特性曲線生成部２８から出力される血管径変化の特性曲線や血流速度特性曲線生成部３８から出力される血流速度変化の特性曲線に基づいて、血管機能評価に必要な各種評価値を演算する。評価値は、例えば、ＦＭＤによる測定における駆血解除から安静時状態に戻るまでの時間、血管径の変化率や血流速度の変化率、血管径と血流速度から求まる血流量やその変化などである。評価値演算部４０で演算された各種評価値は表示画像形成部４２に出力されディスプレイ４４に表示される。
【００２９】
図１に示す超音波診断装置内の各部は、制御部５０により動作制御される。例えば、ユーザ操作による各種指示は制御部５０を介して各部に伝えられる。
【００３０】
次に、図１の超音波診断装置を利用した血管機能評価について説明する。なお、以下の説明において図１に記載した部分については図１の符号を利用する。
【００３１】
図２は、図１の超音波診断装置を利用して得られる血管径変化波形および血流速度変化波形を示す図であり、縦軸は血管径または血流速度を、横軸は計測時刻を示している。図２に示す血管径変化波形６０および血流速度変化波形６２は、ＦＭＤの検査手順によって得られる。そこで、ＦＭＤの検査手順を説明する。
【００３２】
まず、患者の安静時状態をベースラインデータとして、上腕動脈の血管径および血流速度を計測する。つまり、プローブ１０を患者１４の上腕部に当接させて超音波診断を行い、血管径変化波形生成部２４により生成される血管径変化波形６０、および、血流速度変化波形生成部３４により生成される血流速度変化波形６２を取得する。安静時状態の波形期間は図２の期間Ｉである。
【００３３】
次に、患者の前腕部をカフなどで駆血して血管径変化波形６０、および、血流速度変化波形６２を取得する。駆血状態の波形期間は図２の期間ＩＩである。期間ＩＩの波形が示すように、血管径は駆血直後から徐々に減少した後、ほぼ一定値を維持する。一方、血流速度は駆血直後から徐々に減少した後、ほぼ一定値を維持する。駆血は５分程度行われる。
【００３４】
次に、カフを開放して血管径変化波形６０および血流速度変化波形６２を取得する。カフ開放後の波形期間は図２の期間ＩＩＩである。期間ＩＩＩの波形が示すように、血管径はカフ開放直後拡大し、徐々に減少した後、安静時状態の血管径に落ち着いていく。一方、血流速度はカフ開放直後から徐々に増加した後、安静時状態の血流速度に戻っていく。
【００３５】
上記手順により、図２に示す血管径変化波形６０および血流速度変化波形６２が取得される。ＦＭＤによる血管機能評価では、駆血解除から安静時状態に戻るまでの血管径の拡大縮小の様子や血流速度の増加減少の様子などを評価するが、血管径変化波形６０や血流速度変化波形６２の心拍による局所的変動成分の影響を考慮しつつ大局的な時間変化を評価する必要がある。
【００３６】
図３は、血管径変化波形６０および血流速度変化波形６２の部分波形を示す図であり、図３には最大径による血管径変化の特性曲線６４、最大血流速度による血流速度変化の特性曲線６６および心電波形７０も示されている。
【００３７】
図３に示す血管径変化波形６０は、図２の血管径変化波形６０の減少傾向部分、例えば、期間ＩＩＩの後半部分に相当する部分波形であり、一心拍期間内での局所的な増減を伴いつつ開放直後増加し、その後減少傾向を示している。したがって、血管径変化波形６０の大局的な時間変化を評価するためには心拍による局所的変動成分の影響を考慮する必要がある。そこで、血管径特性曲線生成部２８は血管径変化波形６０における一心拍期間内ごとの最大値を求めて、最大径による血管径変化の特性曲線６４を生成する。つまり、血管径特性曲線生成部２８は心電波形７０の連続する二つのＲ波の発生期間であるＲ１からＲ２の期間内において、血管径変化波形６０の最大値Ａ１を算出する。同様にして次のＲ波の発生期間内において最大値Ａ２を、さらに次のＲ波の発生期間内において最大値Ａ３を順次算出する。このようにして次々に算出される最大値Ａ１からＡ５に対して、例えば、スプライン補間を利用して、Ａ１からＡ５を結ぶ曲線である、最大径による血管径変化の特性曲線６４を生成する。なお、血管径変化波形６０における一心拍期間内での局所的な増減は図２の期間ＩＩＩのみではなくその他の期間においても同様である。したがって、図２に示す血管径変化波形６０の全域において、最大径による血管径変化の特性曲線６４が生成される。
【００３８】
一方、図３に示す血流速度変化波形６２は、図２の血流速度変化波形６２の減少傾向部分、例えば、期間ＩＩＩの後半部分に相当する部分波形であり、一心拍期間内での局所的な増減を伴いつつ開放直後増加し、その後減少傾向を示している。したがって、血流速度変化波形６２の大局的な時間傾向を評価するためには心拍による局所的変動成分の影響を考慮する必要がある。そこで、血流速度特性曲線生成部３８は血流速度変化波形６２における一心拍期間内ごとの最大値を求めて、最大血流速度による血流速度変化の特性曲線６６を生成する。つまり、血流速度特性曲線生成部３８は心電波形７０の連続する二つのＲ波の発生期間であるＲ１からＲ２の期間内において、血流速度変化波形６２の最大値Ｂ１を算出する。同様にして次のＲ波の発生期間内において最大値Ｂ２を、さらに次のＲ波の発生期間内において最大値Ｂ３を順次算出する。このようにして次々に算出される最大値Ｂ１からＢ５を結ぶ曲線として、最大血流速度による血流速度変化の特性曲線６６を生成する。なお、血流速度変化波形６２における一心拍期間内での局所的な増減は図２の期間ＩＩＩのみではなくその他の期間においても同様である。したがって、図２に示す血流速度変化波形６２の全域において、最大血流速度による血流速度変化の特性曲線６６が生成される。
【００３９】
以上のようにして得られた血管径変化の特性曲線や血流速度変化の特性曲線に基づいて、評価値演算部４０は血管機能評価のための各種評価値を算出する。評価値の一つに、血管径が駆血解除されてから安静時状態に戻るまでに要する時間が挙げられる。つまり、血管径変化の特性曲線に基づいて、ＦＭＤの検査手順における安静時状態の血管径値Ｄ１（図２参照）を計測し、駆血解除後（図２における期間ＩＩＩの開始時点）から血管径が血管径値Ｄ１に戻るまでの時間を測定する。
【００４０】
別の評価値として血管径の変化率が挙げられる。つまり、血管径変化の特性曲線に基づいて、ユーザにより設定された駆血解除後の任意時刻における血管径値Ｄ３（図２参照）を計測して（Ｄ３−Ｄ２）／Ｄ１×１００（パーセント）を算出する。血管径の変化率として特性曲線の傾きなどを算出してもよい。
【００４１】
この他にも、血流速度の変化率、血管径と血流速度から求まる血流量やその変化なども評価値として算出可能である。いずれの評価値も心拍成分が除去された特性曲線に基づいて算出されるため、心拍成分による影響を抑えた血管機能評価が可能になる。
【００４２】
【発明の効果】
以上説明したように、本発明に係る超音波診断装置により、より高精度な血管機能評価が可能となる。
【図面の簡単な説明】
【図１】本発明に係る超音波診断装置の全体構成を示すブロック図である。
【図２】血管径変化波形および血流速度変化波形を示す図である。
【図３】血管径変化波形および血流速度変化波形の部分拡大波形を示す図である。
【符号の説明】
２２  エコートラッキング処理部、２４  血管径変化波形生成部、２８  血管径特性曲線生成部、３２  血流速度計測部、３４  血流速度変化波形生成部、３８  血流速度特性曲線生成部。
【図１】
2004290408.tif 000003
【図２】
2004290408.tif 000004
【図３】
2004290408.tif 000005