(57)【要約】【課題】超音波の音速制約により生じるドプラスペクトル画像の画質の劣化を抑える。【解決手段】実施系形態に係る超音波診断装置において、設定部は、複数の観測部位を設定する。距離判定部は、前記複数の観測部位のうち少なくとも１つの観測部位の走査線上における深さと所定の閾値とを比較する。スキャン切替部は、前記少なくとも１つの観測部位の走査線上における深さが前記閾値を下回っていた場合に、前記複数の観測部位それぞれに対して１回ずつ交互に超音波を送受信する第１のスキャンを行い、前記少なくとも１つの観測部位の走査線上における深さが前記閾値を上回っていた場合に、前記複数の観測部位のうち少なくとも１つの観測部位については複数回超音波を送受信して、前記複数の観測部位それぞれに対して交互に超音波を送受信する第２のスキャンを行うようにスキャン方式を切り替える。【選択図】図１
2012139489.tif 000002
【特許請求の範囲】
【請求項１】
複数の観測部位を設定する設定部と、
前記複数の観測部位のうち少なくとも１つの観測部位の走査線上における深さと所定の閾値とを比較する距離判定部と、
前記少なくとも１つの観測部位の走査線上における深さが前記閾値を下回っていた場合に、前記複数の観測部位それぞれに対して１回ずつ交互に超音波を送受信する第１のスキャンを行い、前記少なくとも１つの観測部位の走査線上における深さが前記閾値を上回っていた場合に、前記複数の観測部位のうち少なくとも１つの観測部位については複数回超音波を送受信して、前記複数の観測部位それぞれに対して交互に超音波を送受信する第２のスキャンを行うようにスキャン方式を切り替えるスキャン切替部と、
前記第１のスキャン又は前記第２のスキャンにより受信された反射波データに基づいて、前記複数の観測部位それぞれにおける移動速度の経時的な変化を示すドプラスペクトラム画像をそれぞれ生成する画像生成部と、
前記ドプラスペクトル画像を表示する表示部と、
を備える、超音波診断装置。
【請求項２】
前記距離判定部は、前記複数の観測部位のうち少なくとも２つの観測部位の走査線上における深さの合計と所定の閾値とを比較し、
前記スキャン切替部は、前記少なくとも２つの観測部位の走査線上における深さの合計が前記閾値を下回っていた場合に、前記第１のスキャンを行い、前記少なくとも２つの観測部位の走査線上における深さの合計が前記閾値を上回っていた場合に、前記第２のスキャンを行うようにスキャン方式を切り替える、
請求項１に記載の超音波診断装置。
【請求項３】
前記ドプラスペクトラム画像により示される移動速度とから得られる計測値を算出する計測値算出部と、
前記計測値を表示部に表示させる計測値表示部と、
をさらに備える、請求項１又は２に記載の超音波診断装置。
【請求項４】
前記距離判定部は、診断部位又は患者情報に基づいて前記閾値を設定して前記深さの判定を行う、請求項１、２又は３に記載の超音波診断装置。
【請求項５】
前記スキャン切替部は、前記ドプラスペクトル画像の速度レンジが所定の速度閾値を下回った場合に、前記第２のスキャンを行うようにスキャン方式を切り替える、請求項１〜４のいずれか一つに記載の超音波診断装置。
【請求項６】
前記スキャン切替部は、前記ドプラスペクトル画像に折り返しが生じているか否かを検出し、折り返しが生じていることを検出した場合に、前記第２のスキャンを行うようにスキャン方式を切り替える、請求項１〜４のいずれか一つに記載の超音波診断装置。
【請求項７】
超音波診断装置の制御方法であって、
前記超音波診断装置の制御部が、
複数の観測部位を設定し、
前記複数の観測部位のうち少なくとも１つの観測部位の走査線上における深さと所定の閾値とを比較し、
前記少なくとも１つの観測部位の走査線上における深さが前記閾値を下回っていた場合に、前記複数の観測部位それぞれに対して１回ずつ交互に超音波を送受信する第１のスキャンを行い、前記少なくとも１つの観測部位の走査線上における深さが前記閾値を上回っていた場合に、前記複数の観測部位のうち少なくとも１つの観測部位については複数回超音波を送受信して、前記複数の観測部位それぞれに対して交互に超音波を送受信する第２のスキャンを行うようにスキャン方式を切り替え、
前記第１のスキャン又は前記第２のスキャンにより受信された反射波データに基づいて、前記複数の観測部位それぞれにおける移動速度の経時的な変化を示すドプラスペクトラム画像をそれぞれ生成し、
前記ドプラスペクトル画像を表示部に表示させる、
ことを含む、制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、超音波診断装置及びその制御方法に関する。
【背景技術】
【０００２】
従来、Ｂモード画像などの血管像上に血流情報の観察部位としてレンジゲートを設定し、そのレンジゲートにおける血流速度の経時的な変化を示すドプラスペクトラム画像を表示する超音波診断装置が知られている。また、かかる超音波診断装置によって、複数箇所に設定されたレンジゲートそれぞれにおけるドプラスペクトル画像を表示するデュアルドプラ技術も知られている。
【０００３】
ここで、デュアルドプラ技術で用いられるスキャン方式として、インターリーブスキャン及びセグメントスキャンがある。インターリーブスキャンは、複数箇所に設定されたレンジゲートそれぞれに対して１回ずつ交互に超音波を送受信することで、各レンジゲートにおける血流情報を取得する方式である。また、セグメントスキャンは、複数箇所に設定されたレンジゲートそれぞれに対して複数回ずつ交互に超音波を送受信することで、各レンジゲートにおける血流情報を取得する方式である。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開平９−２０６３０３号公報
【特許文献２】特開平６−７３５２号公報
【特許文献３】特開２００８−９２９８１号公報
【特許文献４】特開平１１−９４９３２号公報
【特許文献５】特開平６−７３４８号公報
【特許文献６】特開２００９−１３６４４６号公報
【特許文献７】特開２００７−２０２６１７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、上述した従来技術では、超音波の音速制約によって、良好なドプラスペクトル画像が得られない場合があった。例えば、インターリーブスキャンでは、速度レンジが制約されて低くなるため、被検体の深部に流れる速い血流に関するドプラスペクトラム画像に折り返し現象が発生しやすい。また、セグメントスキャンでは、１つのレンジゲートに連続して超音波を送受信している間は、その他のレンジゲートには超音波が送受信されない。そのため、各レンジゲートにおけるドプラスペクトラム画像に周期的なデータの欠落が発生し、この欠落によって画像の劣化が生じる場合があった。
【課題を解決するための手段】
【０００６】
実施形態に係る超音波診断装置は、設定部と、距離判定部と、スキャン切替部と、画像生成部と、表示部とを備える。設定部は、複数の観測部位を設定する。距離判定部は、前記複数の観測部位のうち少なくとも１つの観測部位の走査線上における深さと所定の閾値とを比較する。スキャン切替部は、前記少なくとも１つの観測部位の走査線上における深さが前記閾値を下回っていた場合に、前記複数の観測部位それぞれに対して１回ずつ交互に超音波を送受信する第１のスキャンを行い、前記少なくとも１つの観測部位の走査線上における深さが前記閾値を上回っていた場合に、前記複数の観測部位のうち少なくとも１つの観測部位については複数回超音波を送受信して、前記複数の観測部位それぞれに対して交互に超音波を送受信する第２のスキャンを行うようにスキャン方式を切り替える。画像生成部は、前記第１のスキャン又は前記第２のスキャンにより受信された反射波データに基づいて、前記複数の観測部位それぞれにおける移動速度の経時的な変化を示すドプラスペクトラム画像をそれぞれ生成する。表示部は、前記ドプラスペクトル画像を表示する。
【図面の簡単な説明】
【０００７】
【図１】図１は、本実施形態に係る超音波診断装置の構成を示すブロック図である。
【図２】図２は、本実施形態に係る制御部の機能構成を示すブロック図である。
【図３Ａ】図３Ａは、本実施形態に係る超音波診断装置におけるシングルドプラモードを説明するための図（１）である。
【図３Ｂ】図３Ｂは、本実施形態に係る超音波診断装置におけるシングルドプラモードを説明するための図（２）である。
【図４Ａ】図４Ａは、本実施形態に係る超音波診断装置におけるシングルドプラモードを説明するための図（１）である。
【図４Ｂ】図４Ｂは、本実施形態に係る超音波診断装置におけるシングルドプラモードを説明するための図（２）である。
【図５Ａ】図５Ａは、本実施形態に係る超音波診断装置におけるデュアルドプラモードを説明するための図（１）である。
【図５Ｂ】図５Ｂは、本実施形態に係る超音波診断装置におけるデュアルドプラモードを説明するための図（２）である。
【図６Ａ】図６Ａは、本実施形態に係る超音波診断装置におけるデュアルドプラモードを説明するための図（１）である。
【図６Ｂ】図６Ｂは、本実施形態に係る超音波診断装置におけるデュアルドプラモードを説明するための図（２）である。
【図７Ａ】図７Ａは、本実施形態に係る超音波診断装置におけるデュアルドプラモードを説明するための図（１）である。
【図７Ｂ】図７Ｂは、本実施形態に係る超音波診断装置におけるデュアルドプラモードを説明するための図（２）である。
【図８】図８は、本実施形態に係る距離判定部による距離の判定を説明するための図である。
【図９】図９は、本実施形態に係るインターリーブスキャンのシーケンスを示す図である。
【図１０】図１０は、本実施形態に係るインターリーブスキャンにおける処理の流れを示す図である。
【図１１】図１１は、左室流入血流ピーク速度及び僧帽弁輪移動速度が選択された場合のインターリーブスキャンのシーケンスを示す図である。
【図１２】図１２は、左室流入血流ピーク速度及び僧帽弁輪移動速度が選択された場合のインターリーブスキャンにおける処理の流れを示す図である。
【図１３】図１３は、本実施形態に係るセグメントスキャンのシーケンスを示す図である。
【図１４】図１４は、本実施形態に係るセグメントスキャンにおける処理の流れを示す図である。
【図１５】図１５は、超音波の音速制約を説明するための図である。
【図１６】図１６は、本実施形態に係る計測値表示部による計測値の表示の一例を示す図である。
【図１７】図１７は、本実施形態に係る計測値算出部による計測値の算出の一例を示す図である。
【図１８】図１８は、本実施形態に係る計測値表示部による計測値の表示の一例を示す図である。
【図１９Ａ】図１９Ａは、本実施形態に係る計測値算出部による計測値の算出の一例を示す図（１）である。
【図１９Ｂ】図１９Ｂは、本実施形態に係る計測値算出部による計測値の算出の一例を示す図（２）である。
【図１９Ｃ】図１９Ｃは、本実施形態に係る計測値算出部による計測値の算出の一例を示す図（３）である。
【図２０】図２０は、本実施形態に係る超音波診断装置によるＢ／Ｄ同時スキャンの処理手順を示すフローチャートである。
【図２１】図２１は、本実施形態に係る超音波診断装置による自動計測処理の処理手順を示すフローチャートである。
【発明を実施するための形態】
【０００８】
まず、本実施形態に係る超音波診断装置の構成について説明する。図１は、本実施形態に係る超音波診断装置１００の構成を示すブロック図である。図１に示すように、本実施形態に係る超音波診断装置１００は、超音波プローブ１と、表示部２と、入力部３と、装置本体１０とを有する。
【０００９】
超音波プローブ１は、複数の圧電振動子を有し、これら複数の圧電振動子は、後述する装置本体１０が有する送信部１１から供給される駆動信号に基づき超音波を発生する。また、超音波プローブ１は、被検体Ｐからの反射波を受信して電気信号に変換する。また、超音波プローブ１は、圧電振動子に設けられる整合層及び音響レンズ、圧電振動子から後方への超音波の伝播を防止するバッキング材などを有する。かかる超音波プローブ１は、装置本体１０と着脱自在に接続される。
【００１０】
超音波プローブ１から被検体Ｐに超音波が送信されると、送信された超音波は、被検体Ｐの体内組織における音響インピーダンスの不連続面で次々と反射され、反射波信号として超音波プローブ１が有する複数の圧電振動子にて受信される。受信される反射波信号の振幅は、超音波が反射される不連続面における音響インピーダンスの差に依存する。なお、送信された超音波パルスが血流や心臓壁などの移動体の表面で反射された場合の反射波信号は、ドプラ効果により、移動体の超音波送信方向に対する速度成分に依存して周波数偏移（ドプラ偏移）を受ける。
【００１１】
なお、本実施形態は、複数の圧電振動子が一列で配置された１次元超音波プローブである超音波プローブ１により、被検体Ｐを２次元でスキャンする場合であっても、１次元超音波プローブの複数の圧電振動子を機械的に揺動する超音波プローブ１や複数の圧電振動子が格子状に２次元で配置された２次元超音波プローブである超音波プローブ１により、被検体Ｐを３次元でスキャンする場合であっても、適用可能である。
【００１２】
入力部３は、マウス、キーボード、ボタン、パネルスイッチ、タッチコマンドスクリーン、フットスイッチ、トラックボールなどを有し、超音波診断装置１００の操作者から各種要求を受け付け、受け付けた各種要求を装置本体１０に転送する。
【００１３】
例えば、操作者は、入力部３が有するトラックボールを用いて、Ｂモード画像などの血管像上に血流情報の観察部位を示すレンジゲートの設定を行なう。また、例えば、操作者は、入力部３が有するパネルスイッチなどを用いて、Ｂモード画像とドプラスペクトル画像とを表示するＢ／Ｄ同時スキャンの開始要求及び終了要求を行う。
【００１４】
表示部２は、超音波診断装置１００の操作者が入力部３を用いて各種要求を入力するためのＧＵＩ（Graphical  User  Interface）を表示したり、装置本体１０において生成された超音波画像などを表示したりする。
【００１５】
装置本体１０は、超音波プローブ１が受信した反射波に基づいて超音波画像を生成する。具体的には、装置本体１０は、送信部１１と、受信部１２と、Ｂモード処理部１３と、ドプラ処理部１４と、画像生成部１５と、画像メモリ１６と、制御部１７と、内部記憶部１８とを有する。
【００１６】
送信部１１は、トリガ発生回路、送信遅延回路及びパルサ回路などを有し、超音波プローブ１に駆動信号を供給する。パルサ回路は、所定の繰り返し周波数（ＰＲＦ：Pulse  Repetition  Frequency）の送信超音波を形成するためのレートパルスを繰り返し発生する。なお、ＰＲＦは、レート周波数とも呼ばれる。送信遅延回路は、パルサ回路が発生する各レートパルスに対して、超音波プローブ１から発生される超音波をビーム状に集束して送信指向性を決定するために必要な圧電振動子ごとの送信遅延時間を与える。また、トリガ発生回路は、レートパルスに基づくタイミングで、超音波プローブ１に駆動信号（駆動パルス）を印加する。すなわち、送信遅延回路は、各レートパルスに対し与える送信遅延時間を変化させることで、圧電振動子面からの送信方向を任意に調整する。
【００１７】
なお、送信部１１は、後述する制御部１７の指示に基づいて、所定のスキャンシーケンスを実行するために、送信周波数、送信駆動電圧などを瞬時に変更可能な機能を有している。特に、送信駆動電圧の変更は、瞬間にその値を切り替え可能なリニアアンプ型の発信回路、又は、複数の電源ユニットを電気的に切り替える機構によって実現される。
【００１８】
ここで、送信遅延時間は、超音波ビームの送信フォーカスの音響レンズからの位置（深さ）によって決定される。そして、送信部１１は、送信遅延時間を用いることで、超音波の送信における送信指向性を制御する。
【００１９】
受信部１２は、アンプ回路、Ａ／Ｄ変換器、受信遅延回路、加算器などを有し、超音波プローブ１が受信した反射波信号に対して各種処理を行なって反射波データを生成する。アンプ回路は、反射波信号をチャンネルごとに増幅してゲイン補正処理を行なう。Ａ／Ｄ変換器は、ゲイン補正された反射波信号をＡ／Ｄ変換する。受信遅延回路は、デジタルデータに受信指向性を決定するのに必要な受信遅延時間を与える。加算器は、受信遅延回路により受信遅延時間が与えられた反射波信号の加算処理を行なって反射波データを生成する。加算器の加算処理により、反射波信号の受信指向性に応じた方向からの反射成分が強調される。
【００２０】
ここで、受信遅延時間は、超音波ビームの受信フォーカスの音響レンズからの位置（深さ）によって決定される。そして、受信部１２は、受信遅延時間を用いることで、超音波の受信における受信指向性を制御する。
【００２１】
また、本実施形態にかかる超音波プローブ１は、送信フォーカス及び受信フォーカスの位置に応じて、送受信に用いる圧電振動子（送信用口径及び受信用口径）を変更することが可能である。例えば、近い位置からの反射波信号を受信する際には、強い受信フォーカスをかけるために、受信する振動子の数を少なくしておき、中央部分の圧電振動子で受信した反射波信号のみが超音波画像の生成に用いられるように、小さな受信用口径が受信条件として決定される。また、遠い位置からの反射波信号を受信する際には、圧電振動子の口径が大きいほど受信フォーカスを強くできるので、距離に応じて受信用口径を大きくするように受信条件が決定される。
【００２２】
Ｂモード処理部１３は、受信部１２によって生成された反射波データに対して対数増幅、包絡線検波処理などを行うことで、信号強度が輝度の明るさで表現されるデータ（Ｂモードデータ）を生成する。
【００２３】
ドプラ処理部１４は、受信部１２によって生成された反射波データから速度情報を周波数解析することでドプラ偏移を抽出し、ドプラ偏移を用いることで、ドプラ効果による血流や組織、造影剤エコー成分を抽出し、平均速度、分散、パワーなどの移動体情報を多点について抽出したデータ（ドプラデータ）を生成する。
【００２４】
なお、本実施形態に係るＢモード処理部１３及びドプラ処理部１４は、２次元の反射波データ及び３次元の反射波データの両方について処理可能である。
【００２５】
画像生成部１５は、Ｂモード処理部１３及びドプラ処理部１４によって生成されたデータから超音波画像を生成する。すなわち、画像生成部１５は、Ｂモード処理部１３によって生成されたＢモードデータから、反射波の強度を輝度で表したＢモード画像を生成する。または、画像生成部１５は、Ｂモード処理部１３によって生成された所定のスキャンラインにおけるＢモードデータから、所定のスキャンラインにおける反射波強度の時系列に沿った変化を輝度にて表したＭモード画像を生成する。
【００２６】
また、画像生成部１５は、ドプラ処理部１４によって生成されたドプラデータから、移動体情報（血流情報や組織の移動情報）を表す平均速度画像、分散画像、パワー画像、又は、これらの組み合わせ画像としてのカラードプラ画像を生成する。さらに、画像生成部１５は、ドプラ処理部１４によって生成されたドプラデータから、移動体の速度情報（血流の速度情報や組織の速度情報）を時系列に沿ってプロットしたドプラスペクトラム画像を生成する。
【００２７】
画像メモリ１６は、画像生成部１５が生成した超音波画像を記憶するメモリである。また、画像メモリ１６は、Ｂモード処理部１３やドプラ処理部１４によって生成されたデータを記憶することも可能である。
【００２８】
内部記憶部１８は、超音波送受信、画像処理及び表示処理を行なうための制御プログラムや、診断情報（例えば、患者ＩＤ、医師の所見など）、診断プロトコル、各種ボディーマークなどの各種データを記憶する。また、内部記憶部１８は、必要に応じて、画像メモリ１６によって記憶される画像の保管などにも使用される。また、内部記憶部１８によって記憶されるデータは、図示しないインターフェースを経由して、外部の周辺装置へ転送することができる。
【００２９】
制御部１７は、超音波診断装置１００の処理全体を制御する。具体的には、制御部１７は、入力部３を介して操作者から入力された各種要求や、内部記憶部１８から読込んだ各種制御プログラム及び各種データに基づいて、送信部１１、受信部１２、Ｂモード処理部１３、ドプラ処理部１４、画像生成部１５の処理を制御する。また、制御部１７は、画像メモリ１６によって記憶される超音波画像や、画像生成部１５によって行われる各種処理を指定するためのＧＵＩなどを表示部２に表示するよう制御する。
【００３０】
以上、本実施形態に係る超音波診断装置１００の構成について説明した。このような構成のもと、本実施形態に係る超音波診断装置１００では、制御部１７が、血流情報の観測部位として設定された少なくとも２つのレンジゲートについて、第１のレンジゲートから超音波プローブまでの距離と第２のレンジゲートから超音波プローブまでの距離との合計長が閾値未満であるか否かを判定する。また、制御部１７は、距離の合計長が閾値未満であると判定された場合にはインターリーブスキャンを行い、距離の合計長が閾値以上であると判定された場合にはセグメントスキャンを行うようにスキャン方式を切り替える。また、画像生成部１５が、セグメントスキャン又はインターリーブスキャンにより受信された反射波データに基づいて、第１のレンジゲートにおける血流速度の経時的な変化を示す第１のドプラスペクトラム画像と第２のレンジゲートにおける血流速度の経時的な変化を示す第２のドプラスペクトル画像とをそれぞれ生成する。そして、表示部２が、画像生成部１５により生成された第１のドプラスペクトラム画像及び第２のドプラスペクトル画像を表示する。
【００３１】
すなわち、本実施形態に係る超音波診断装置１００は、複数箇所に設定されたレンジゲートそれぞれにおけるドプラスペクトル画像を表示する場合に、各レンジゲートの深さの合計に応じて、インターリーブスキャンとセグメントスキャンとを自動的に切り替える。ここで、インターリーブスキャンは、第１のレンジゲート及び第２のレンジゲートそれぞれに対して１回ずつ交互に超音波を送受信する方式である。また、セグメントスキャンは、第１のレンジゲート及び第２のレンジゲートそれぞれに対して複数回ずつ交互に超音波を送受信する方式である。
【００３２】
以下では、かかる超音波診断装置１００について詳細に説明する。なお、本実施形態では、超音波診断装置１００は、Ｂモード画像の血管像上に血流情報の観察部位として２つのレンジゲートを設定し、各レンジゲートにおけるドプラスペクトル画像をそれぞれ表示する。このように、２つのレンジゲートそれぞれにおけるドプラスペクトル画像を表示する表示モードを以下ではデュアルドプラモードと呼ぶ。なお、超音波診断装置１００は、２つのレンジゲートにおけるドプラスペクトル画像を１つずつ表示することも可能である。このように、２つのレンジゲートにおけるドプラスペクトル画像を１つずつ表示する表示モードを以下ではシングルドプラモードと呼ぶ。
【００３３】
また、超音波診断装置１００は、診断対象の臓器や診断の種類に応じて各種のアプリケーションを実行することができる。本実施形態では、超音波診断装置１００が、心臓の診断用のアプリケーションと頚動脈の診断用のアプリケーションとを実行する場合について説明する。さらに、超音波診断装置１００は、診断部位に応じて、ドプラスペクトル画像の表示態様を切り替えることができる。本実施形態では、診断部位が心臓の左室流入血流（Left  Ventricular  Inflow：ＬＶＩ）及び左室流出血流（Left  Ventricular  Outflow：ＬＶＯ）である場合と、心臓の左室流入血流ピーク速度（Ｅ）及び僧帽弁輪移動速度（ｅ’）である場合と、頚動脈の総頚動脈（Common  Carotid  Artery：ＣＣＡ）及び内頚動脈（Internal  Carotid  Artery：ＩＣＡ）である場合とについて説明する。
【００３４】
次に、本実施形態に係る制御部１７について詳細に説明する。図２は、本実施形態に係る制御部１７の機能構成を示すブロック図である。図２に示すように、制御部１７は、表示制御部１７ａと、設定部１７ｆと、距離判定部１７ｂと、スキャン切替部１７ｃと、計測値算出部１７ｄと、計測値表示部１７ｅとを有する。
【００３５】
表示制御部１７ａは、入力部３を介して操作者から各種要求を受け付け、受け付けた各種要求に応じて、画像メモリ１６によって記憶される超音波画像や、画像生成部１５によって行われる各種処理を指定するためのＧＵＩなどを表示部２に表示させる。また、超音波診断装置１００は、入力部３が有するタッチコマンドスクリーンを介して、上述した表示モード、アプリケーション及び診断部位を選択する操作を操作者から受け付ける。
【００３６】
例えば、表示制御部１７ａは、タッチスクリーン上に「Ｄｕａｌ  Ｄｏｐｐｌｅｒ」ボタン、「ＰＷＤ１」ボタン、「ＰＷＤ２」ボタンを表示させる。「Ｄｕａｌ  Ｄｏｐｐｌｅｒ」ボタンは、シングルモード又はデュアルモードの選択と、診断部位の選択とを操作者から受け付けるためのボタンである。この「Ｄｕａｌ  Ｄｏｐｐｌｅｒ」は、操作者によって押下されるたびに、「Ｄｕａｌ  Ｄｏｐｐｌｅｒ（ｏｆｆ）」、「Ｄｕａｌ  Ｄｏｐｐｌｅｒ（ＬＶＩ／ＬＶＯ）」、「Ｄｕａｌ  Ｄｏｐｐｌｅｒ（Ｅ／ｅ’）」の順で表示を切り替える。
【００３７】
また、「ＰＷＤ１」ボタン及び「ＰＷＤ２」ボタンは、２つのレンジゲートのうちいずれかを選択する操作を操作者から受け付けるためのボタンである。これら「ＰＷＤ１」ボタン及び「ＰＷＤ２」ボタンは、「Ｄｕａｌ  Ｄｏｐｐｌｅｒ」ボタンが「Ｄｕａｌ  Ｄｏｐｐｌｅｒ（ｏｆｆ）」の状態では「ＰＷＤ１」及び「ＰＷＤ２」と表示され、「Ｄｕａｌ  Ｄｏｐｐｌｅｒ」ボタンが「Ｄｕａｌ  Ｄｏｐｐｌｅｒ（ＬＶＩ／ＬＶＯ）」の状態では「ＰＷＤ１（ＬＶＩ）」及び「ＰＷＤ２（ＬＶＯ）」と表示され、「Ｄｕａｌ  Ｄｏｐｐｌｅｒ」ボタンが「Ｄｕａｌ  Ｄｏｐｐｌｅｒ（Ｅ／ｅ’）」の状態では「ＰＷＤ１（Ｅ）」及び「ＰＷＤ２（ｅ’）」と表示される。
【００３８】
そして、例えば、表示制御部１７ａは、操作者からＢ／Ｄ同時スキャンの開始要求を受け付けた場合には、画像生成部１５によって生成されたＢモード画像及びドプラスペクトラム画像を表示部２に表示させる。また、表示制御部１７ａは、表示部２に表示されたＢモード画像上に、超音波の送受信方向を示す２本のスキャンラインを表示する。また、表示制御部１７ａは、各スキャンライン上にレンジゲートを表示する。なお、表示制御部１７ａは、入力部３が有するトラックボールを介して操作者から受け付けた操作に応じて、各スキャンラインを走査方向に移動させたり、各レンジゲートの位置をスキャンラインに沿って移動させたりする。
【００３９】
ここで、表示制御部１７ａは、操作者によって選択された表示モードやアプリケーション、診断部位に応じて、Ｂモード画像上に表示するスキャンライン及びレンジゲートの位置や、ドプラスペクトラム画像の種類を変化させる。例えば、スキャンライン及びレンジゲートの位置は、あらかじめアプリケーション及び診断部位ごとに定義されたプリセット情報に基づいて決められる。
【００４０】
図３Ａ、３Ｂ、４Ａ及び４Ｂは、本実施形態に係る超音波診断装置１００におけるシングルドプラモードを説明するための図である。図３Ａ、３Ｂ、４Ａ及び４Ｂは、心臓の診断用のアプリケーションが選択され、診断部位として左室流入血流及び左室流出血流が選択された場合を示している。なお、図３Ａ及び４Ａは、表示部２が有する表示領域を示しており、図３Ｂ及び４Ｂは、タッチコマンドスクリーンを示している。
【００４１】
図３Ａ、３Ｂ、４Ａ及び４Ｂに示すように、心臓の診断用のアプリケーションが選択され、診断部位として左室流入血流及び左室流出血流が選択された場合には、表示制御部１７ａは、表示部２にＢモード画像３１を表示し、そのＢモード画像３１上に２本のスキャンラインＰＷＤ１及びＰＷＤ２を表示する。また、表示制御部１７ａは、スキャンラインＰＷＤ１上にレンジゲートＲＧ１を表示し、スキャンラインＰＷＤ２上にレンジゲートＲＧ２を表示する。このとき、表示制御部１７ａは、左室流入血流の位置にレンジゲートＲＧ１が配置され、左室流出血流の位置にレンジゲートＲＧ２が配置されるように、スキャンラインＰＷＤ１及びＰＷＤ２並びにレンジゲートＲＧ１及びＲＧ２を表示する。
【００４２】
そして、例えば、図３Ａ及び３Ｂに示すように、タッチコマンドスクリーンに表示された「Ｄｕａｌ  Ｄｏｐｐｌｅｒ」ボタンが「Ｄｕａｌ  Ｄｏｐｐｌｅｒ（ｏｆｆ）」の状態で、「ＰＷＤ１」ボタンが押下された場合には、表示制御部１７ａは、表示部２が有する表示領域上に、スキャンラインＰＷＤ１上に設定されたレンジゲートＲＧ１におけるドプラスペクトル画像３２を表示させる。なお、この状態では、表示制御部１７ａは、スキャンラインＰＷＤ１及びレンジゲートＲＧ１に対する操作を受け付け可能な状態にする。
【００４３】
また、例えば、図４Ａ及び４Ｂに示すように、タッチコマンドスクリーンに表示された「Ｄｕａｌ  Ｄｏｐｐｌｅｒ」ボタンが「Ｄｕａｌ  Ｄｏｐｐｌｅｒ（ｏｆｆ）」の状態で、「ＰＷＤ２」ボタンが選択された場合には、表示制御部１７ａは、表示部２が有する表示領域上に、スキャンラインＰＷＤ２上に設定されたレンジゲートＲＧ２におけるドプラスペクトル画像４２を表示させる。なお、この状態では、表示制御部１７ａは、スキャンラインＰＷＤ２及びレンジゲートＲＧ２に対する操作を受け付け可能な状態にする。
【００４４】
図５Ａ、５Ｂ、６Ａ、６Ｂ、７Ａ、７Ｂ及び７Ｃは、本実施形態に係る超音波診断装置１００におけるデュアルドプラモードを説明するための図である。図５Ａ及び５Ｂは、心臓の診断用のアプリケーションが選択され、診断部位として左室流入血流及び左室流出血流が選択された場合を示している。また、図６Ａ及び６Ｂは、心臓の診断用のアプリケーションが選択され、診断部位として左室流入血流ピーク速度及び僧帽弁輪移動速度が選択された場合を示している。また、図７Ａ、７Ｂ及び７Ｃは、頚動脈の診断用のアプリケーションが選択され、診断部位として総頚動脈及び内頚動脈が選択された場合を示している。
【００４５】
図５Ａ及び５Ｂに示すように、心臓の診断用のアプリケーションが選択され、診断部位として左室流入血流及び左室流出血流が選択された場合には、表示制御部１７ａは、表示部２に心臓のＢモード画像５１を表示し、そのＢモード画像５１上に２本のスキャンラインＰＷＤ１及びＰＷＤ２を表示する。また、表示制御部１７ａは、スキャンラインＰＷＤ１上にレンジゲートＲＧ１を表示し、スキャンラインＰＷＤ２上にレンジゲートＲＧ２を表示する。このとき、表示制御部１７ａは、左室流入血流の位置にレンジゲートＲＧ１が配置され、左室流出血流の位置にレンジゲートＲＧ２が配置されるように、スキャンラインＰＷＤ１及びＰＷＤ２並びにレンジゲートＲＧ１及びＲＧ２を表示する。
【００４６】
そして、例えば、図５Ａ及び５Ｂに示すように、タッチコマンドスクリーンに表示された「Ｄｕａｌ  Ｄｏｐｐｌｅｒ」ボタンが「Ｄｕａｌ  Ｄｏｐｐｌｅｒ（ＬＶＩ／ＬＶＯ）」の状態である場合には、表示制御部１７ａは、表示部２が有する表示領域上に、レンジゲートＲＧ１における正側の速度成分を示すドプラスペクトル画像５２と、レンジゲートＲＧ２における負側の速度成分を示すドプラスペクトル画像５３とを上下に配置して表示させる。なお、この状態で、「ＰＷＤ１」ボタンが押下された場合には、表示制御部１７ａは、スキャンラインＰＷＤ１及びレンジゲートＲＧ１に対する操作を受け付け可能な状態にする。一方、「ＰＷＤ２」ボタンが押下された場合には、表示制御部１７ａは、スキャンラインＰＷＤ２及びレンジゲートＲＧ２に対する操作を受け付け可能な状態にする。
【００４７】
また、図６Ａ及び６Ｂに示すように、心臓の診断用のアプリケーションが選択され、診断部位として左室流入血流ピーク速度及び僧帽弁輪移動速度が選択された場合には、表示制御部１７ａは、表示部２に心臓のＢモード画像６１を表示し、そのＢモード画像６１上に２本のスキャンラインＰＷＤ１及びＰＷＤ２を表示する。また、表示制御部１７ａは、スキャンラインＰＷＤ１上にレンジゲートＲＧ１を表示し、スキャンラインＰＷＤ２上にレンジゲートＲＧ２を表示する。このとき、表示制御部１７ａは、左室流入血流の位置にレンジゲートＲＧ１が配置され、僧帽弁輪の位置にレンジゲートＲＧ２が配置されるように、スキャンラインＰＷＤ１及びＰＷＤ２並びにレンジゲートＲＧ１及びＲＧ２を表示する。
【００４８】
そして、例えば、図６Ａ及び６Ｂに示すように、タッチコマンドスクリーンに表示された「Ｄｕａｌ  Ｄｏｐｐｌｅｒ」ボタンが「Ｄｕａｌ  Ｄｏｐｐｌｅｒ（Ｅ／ｅ’）」の状態である場合には、表示制御部１７ａは、表示部２が有する表示領域上に、レンジゲートＲＧ１における左室流入血流ピーク速度のドプラスペクトル画像６２と、僧帽弁輪移動速度のドプラスペクトル画像６３とを上下に配置して表示させる。なお、この状態で、「ＰＷＤ１」ボタンが押下された場合には、表示制御部１７ａは、スキャンラインＰＷＤ１及びレンジゲートＲＧ１に対する操作を受け付け可能な状態にする。一方、「ＰＷＤ２」ボタンが押下された場合には、表示制御部１７ａは、スキャンラインＰＷＤ２及びレンジゲートＲＧ２に対する操作を受け付け可能な状態にする。
【００４９】
また、例えば、図７Ａ及び７Ｂに示すように、頚動脈の診断用のアプリケーションが選択され、診断部位として総頚動脈及び内頚動脈が選択された場合には、表示制御部１７ａは、表示部２に頚動脈のＢモード画像７１を表示し、そのＢモード画像７１上に２本のスキャンラインＰＷＤ１及びＰＷＤ２を表示する。また、表示制御部１７ａは、スキャンラインＰＷＤ１上にレンジゲートＲＧ１を表示し、スキャンラインＰＷＤ２上にレンジゲートＲＧ２を表示する。このとき、表示制御部１７ａは、総頚動脈の位置にレンジゲートＲＧ１が配置され、内頚動脈の位置にレンジゲートＲＧ２が配置されるように、スキャンラインＰＷＤ１及びＰＷＤ２並びにレンジゲートＲＧ１及びＲＧ２を表示する。
【００５０】
そして、例えば、図７Ａ及び７Ｂに示すように、タッチコマンドスクリーンに表示された「Ｄｕａｌ  Ｄｏｐｐｌｅｒ」ボタンが「Ｄｕａｌ  Ｄｏｐｐｌｅｒ（ＣＣＡ／ＩＣＡ）」の状態である場合には、表示制御部１７ａは、表示部２が有する表示領域上に、レンジゲートＲＧ１における総頚動脈のドプラスペクトル画像７２と、内頚動脈のドプラスペクトル画像７３とを上下に配置して表示させる。なお、この状態で、「ＰＷＤ１」ボタンが押下された場合には、表示制御部１７ａは、スキャンラインＰＷＤ１及びレンジゲートＲＧ１に対する操作を受け付け可能な状態にする。一方、「ＰＷＤ２」ボタンが押下された場合には、表示制御部１７ａは、スキャンラインＰＷＤ２及びレンジゲートＲＧ２に対する操作を受け付け可能な状態にする。
【００５１】
図２の説明にもどって、設定部１７ｆは、複数の観測部位を設定する。本実施形態では、設定部１７ｆは、表示制御部１７ａによって表示部２に表示されたレンジゲートの位置に基づいて、観測部位を設定する。具体的には、設定部１７ｆは、表示部２に表示されたＢモード画像上でレンジゲートが位置付けられた箇所を、観測部位として設定する。
【００５２】
距離判定部１７ｂは、複数の観測部位のうち少なくとも１つの観測部位の走査線上における深さと所定の閾値とを比較する。本実施形態では、距離判定部１７ｂは、複数の観測部位のうち少なくとも２つの観測部位の走査線上における深さの合計と所定の閾値とを比較する。
【００５３】
具体的には、距離判定部１７ｂは、血流情報の観測部位として設定された少なくとも２つのレンジゲートについて、第１のレンジゲートから超音波プローブまでの距離と第２のレンジゲートから超音波プローブまでの距離との合計長が閾値未満であるか否かを判定する。
【００５４】
図８は、本実施形態に係る距離判定部１７ｂによる距離の判定を説明するための図である。図８に示すように、例えば、Ｂモード画像８１上に設定された２本のスキャンラインＰＷＤ１及びＰＷＤ２が設定されており、スキャンラインＰＷＤ１上にレンジゲートＲＧ１が設定され、スキャンラインＰＷＤ２上にレンジゲートＲＧ２が設定されていたとする。この場合には、距離判定部１７ｂは、超音波プローブ１のプローブ原点８０からレンジゲートＲＧ１までの距離Ｒ１と、超音波プローブ１の原点８０からレンジゲートＲＧ２までの距離Ｒ２とをそれぞれ算出する。そして、距離判定部１７ｂは、算出した距離Ｒ１と距離Ｒ２との合計長を算出し、その合計長が所定の閾値未満であるか否かを判定する。
【００５５】
ここで、本実施形態では、距離判定部１７ｂは、診断部位に基づいて閾値を設定して距離の合計長の判定を行う。例えば、距離判定部１７ｂは、診断部位が心臓の左室流入血流及び左室流出血流である場合には、インターリーブスキャンを行った場合でも折り返しを生じないレンジゲートの深さの２倍の値を閾値として設定する。ここで、折り返しを生じないレンジゲートの深さは、例えば、あらかじめ実験的にレンジゲートの深さを少しずつ深くしながらインターリーブスキャンを行い、ドプラスペクトル画像に折り返しが発生した時点でのレンジゲートの深さより浅くすればよい。この深さから求められる閾値は、例えば、あらかじめ操作者によって所定の記憶部に格納される。そして、距離判定部１７ｂは、記憶部に格納されている閾値を取得して、距離の合計長の判定を行う。なお、診断部位が心臓の左室流入血流ピーク速度及び僧帽弁輪移動速度である場合についても、同様に閾値を設定することができる。
【００５６】
また、例えば、距離判定部１７ｂは、診断部位が頚動脈の総頚動脈及び内頚動脈である場合には、設定し得るレンジゲートの深さの最大値の２倍より大きな値を閾値として設定する。これにより、頚動脈の総頚動脈及び内頚動脈の診断が行われる場合には、各レンジゲートから超音波プローブまでの距離の合計長が閾値以上となることがなくなるので、常にインターリーブスキャンでデータが収集されることになる。一般的に、頚動脈は体表から浅い位置にあるので、インターリーブスキャンでデータを収集しても、折り返し現象が発生する可能性が低く、十分な画質のドプラ画像が得られる。
【００５７】
なお、ここでは、診断部位に基づいて閾値を設定する場合について説明したが、例えば、距離判定部１７ｂは、患者情報に基づいて閾値を設定してもよい。例えば、距離判定部１７ｂは、診断が行われる際に操作者によって超音波診断装置１００に入力された患者の性別や年齢に基づいて、閾値を設定する。例えば、ドプラの速度レンジは、加齢とともに低下することが知られている。そこで、例えば、距離判定部１７ｂは、患者の年齢が上がるにつれて値が低くなるように、閾値を設定する。
【００５８】
図２の説明にもどって、スキャン切替部１７ｃは、少なくとも１つの観測部位の走査線上における深さが閾値を下回っていた場合に、複数の観測部位それぞれに対して１回ずつ交互に超音波を送受信する第１のスキャンを行い、少なくとも１つの観測部位の走査線上における深さが閾値を上回っていた場合に、複数の観測部位のうち少なくとも１つの観測部位については複数回超音波を送受信して、複数の観測部位それぞれに対して交互に超音波を送受信する第２のスキャンを行うようにスキャン方式を切り替える。
【００５９】
本実施形態では、スキャン切替部１７ｃは、少なくとも２つの観測部位の走査線上における深さの合計が閾値を下回っていた場合に、第１のスキャンを行い、少なくとも２つの観測部位の走査線上における深さの合計が閾値を上回っていた場合に、第２のスキャンを行うようにスキャン方式を切り替える。
【００６０】
なお、第２のスキャンは、例えば、複数の観測部位それぞれに対して同じ数だけ超音波を送受信するものでもよいし、複数の観測部位それぞれに対して異なる数だけ超音波を送受信するものでもよい。例えば、第２のスキャンは、複数の観測部位のうち１つ又は複数の観測部位については１回ずつ超音波を送受信し、他の観測部位については複数回超音波を送受信してもよい。
【００６１】
ここで、各観測部位に何回ずつ送受信を行うかは、例えば、観測部位の深さによって決められる。一般的に、一方の観測部位の深さが深いほど、その観測部位に対して複数回超音波を送受信するのにかかる時間が長くなるため、他方の観測部位に関するドプラ波形に生じる空隙が大きくなってしまう。そこで、例えば、深い位置にある観測部位については、浅い位置にある観測部位よりも送受信の回数を減らす。
【００６２】
また、送受信の回数は、例えば、必要とされる測定精度によって決められてもよい。例えば、高精度の測定が必要な観測部位や、ＳＮ比の悪い観測部位については、送受信の回数を増やすようにする。また、送受信の回数は、例えば、流速によって決められてもよい。例えば、ある観測部位の流速が低い場合には、その観測部位については１回ずつ超音波を送受信し、他の観測部位については複数回超音波を送受信するようにする。これにより、流速が低い観測部位については、時間間隔を空けて超音波を送受信することで、低い流速の検出を行うことができる。
【００６３】
具体的には、スキャン切替部１７ｃは、距離判定部１７ｂにより距離の合計長が閾値未満であると判定された場合にインターリーブスキャンを行い、距離判定部１７ｂにより距離の合計長が閾値以上であると判定された場合にセグメントスキャンを行うようにスキャン方式を切り替える。ここで、インターリーブスキャン及びセグメントスキャンについて具体的に説明する。なお、ここでは、図８に示したレンジゲートＲＧ１及びＲＧ２からデータを収集する場合について説明する。
【００６４】
まず、インターリーブスキャンについて説明する。図９は、本実施形態に係るインターリーブスキャンのシーケンスを示す図である。図９において、横軸は時間を示している。また、Ｔｘは、超音波プローブ１から送信される超音波のＰＲＦと送信タイミングとを示している。また、Ｒｘは、超音波プローブ１により反射波が受信されるタイミングを示している。また、Ｄ１は、ＲＧ１におけるドプラスペクトル画像用のデータがサンプリングされるタイミングを示している。また、Ｄ２は、ＲＧ２におけるドプラスペクトル画像用のデータがサンプリングされるタイミングを示している。
【００６５】
インターリーブスキャンでは、レンジゲートＲＧ１及びレンジゲートＲＧ２それぞれに対して１回ずつ交互に超音波が送受信される。例えば、図９に示すように、インターリーブスキャンでは、スキャンラインＰＷＤ１に沿ってＰＲＦが８ｋＨｚの超音波が送信され、スキャンラインＰＷＤ２に沿ってＰＲＦが４ｋＨｚの超音波が送信される。ここで、スキャンラインＰＷＤ１に対する送信とスキャンラインＰＷＤ２に対する送信とは、１回ずつ交互に行われる。
【００６６】
また、インターリーブスキャンでは、例えば、レンジゲートＲＧ１の反射波とレンジゲートＲＧ２の反射波とが交互に受信される。そして、例えば、ＲＧ１におけるドプラスペクトル画像用のデータとＲＧ２におけるドプラスペクトル画像用のデータとが、それぞれ２．７ｋＨｚの周期でサンプリングされる。なお、インターリーブスキャンでは、レンジゲートＲＧ１及びＲＧ２の位置に応じて、各レンジゲートからのデータ収集が最短で可能となるＰＲＦが設定される。
【００６７】
図１０は、本実施形態に係るインターリーブスキャンにおける処理の流れを示す図である。図１０に示すように、インターリーブスキャンでは、ドプラ処理部１４が、レンジゲートＲＧ１からの反射波データに対して、ウォールフィルター（Wall  Filter）、高速フーリエ変換（Fast  Fourier  Transformation：ＦＦＴ）、後処理（Ｐｏｓｔ処理）を順に施すことで、レンジゲートＲＧ１における血流の流速を示すドプラデータを生成する。
【００６８】
一方、ドプラ処理部１４は、レンジゲートＲＧ２からの反射波データに対して、ウォールフィルター（Wall  Filter）、ＦＦＴ（Fast  Fourier  Transformation）、後処理（Post処理）を順に施すことで、レンジゲートＲＧ２における血流の流速を示すドプラデータを生成する。そして、画像生成部１５が、ドプラ処理部１４によって生成された各ドプラデータから、レンジゲートＲＧ１におけるドプラスペクトラム画像とレンジゲートＲＧ２におけるドプラスペクトル画像とをそれぞれ生成して表示部２に表示させる（Ｄｕａｌ−Ｄ表示）。
【００６９】
なお、診断部位として心臓の左室流入血流ピーク速度及び僧帽弁輪移動速度が選択された場合のインターリーブスキャンでは、僧帽弁輪移動速度のドプラスペクトル画像は組織ドプラになるので、シーケンスが図９に示したものと少し異なる。図１１は、左室流入血流ピーク速度及び僧帽弁輪移動速度が選択された場合のインターリーブスキャンのシーケンスを示す図である。なお、ここでは、左室流入血流の位置にレンジゲートＲＧ１が配置され、僧帽弁輪の位置にレンジゲートＲＧ２が配置されているとする。図１１において、横軸は時間を示している。また、Ｔｘ、Ｒｘ、Ｄ１及びＤ２の意味は図９と同様である。また、Ｄ３は、ＲＧ２における僧帽弁輪移動速度のドプラスペクトル画像用のデータがサンプリングされるタイミングを示している。
【００７０】
左室流入血流ピーク速度及び僧帽弁輪移動速度が選択された場合のインターリーブスキャンでは、レンジゲートＲＧ１及びレンジゲートＲＧ２それぞれに対して複数回ずつ交互に超音波が送受信される。例えば、図１１に示すように、インターリーブスキャンでは、スキャンラインＰＷＤ１に沿ってＰＲＦが５ｋＨｚの超音波が送信され、スキャンラインＰＷＤ２に沿ってＰＲＦが４ｋＨｚの超音波が送信される。ここで、スキャンラインＰＷＤ１に対する送信とスキャンラインＰＷＤ２に対する送信とは、１回ずつ交互に行われる。
【００７１】
また、インターリーブスキャンでは、例えば、レンジゲートＲＧ１の反射波とレンジゲートＲＧ２の反射波とが交互に受信される。そして、例えば、ＲＧ１におけるドプラスペクトル画像用のデータとＲＧ２におけるドプラスペクトル画像用のデータとが、それぞれ２．２ｋＨｚの周期でサンプリングされる。また、ＲＧ２における僧帽弁輪移動速度のドプラスペクトル画像用のデータは、例えば、１．１ｋＨｚの周期に間引いて収集が行われる。これは、組織の移動速度は血流の流速に比べて遅いためである。なお、インターリーブスキャンでは、レンジゲートＲＧ１及びＲＧ２の位置に応じて、各レンジゲートからのデータ収集が最短で可能となるＰＲＦが設定される。
【００７２】
図１２は、左室流入血流ピーク速度及び僧帽弁輪移動速度が選択された場合のインターリーブスキャンにおける処理の流れを示す図である。図１２に示すように、左室流入血流ピーク速度及び僧帽弁輪移動速度が選択された場合のインターリーブスキャンでは、ドプラ処理部１４が、レンジゲートＲＧ２からの反射波データに対して、ウォールフィルター（Wall  Filter）を施す前に、ローパスフィルタ（Low  Pass  Filter：ＬＰＦ）とスケーリング（Ｓｃａｌｉｎｇ）とを施す。これにより、僧帽弁輪移動速度のドプラスペクトル画像用のデータが間引かれる。
【００７３】
次に、セグメントスキャンについて説明する。図１３は、本実施形態に係るセグメントスキャンのシーケンスを示す図である。図１３において、横軸は時間を示している。また、Ｔｘ及びＲｘの意味は図９と同じである。また、Ｄ１は、ＲＧ１におけるドプラスペクトル画像用のデータがサンプリングされるタイミングを示している。また、Ｄ２は、ＲＧ２におけるドプラスペクトル画像用のデータがサンプリングされるタイミングを示している。また、Ｄ３は、ＲＧ１におけるドプラスペクトル画像用のデータに関する信号処理を示している。また、Ｄ４は、ＲＧ２におけるドプラスペクトル画像用のデータに関する信号処理を示している。
【００７４】
セグメントスキャンでは、レンジゲートＲＧ１及びレンジゲートＲＧ２それぞれに対して複数回ずつ交互に超音波が送受信される。例えば、図１３に示すように、セグメントスキャンでは、スキャンラインＰＷＤ１に沿ってＰＲＦが５ｋＨｚの超音波が複数回連続して送信され、スキャンラインＰＷＤ２に沿ってＰＲＦが４ｋＨｚの超音波が複数回連続して送信される。ここで、スキャンラインＰＷＤ１に対する送信とスキャンラインＰＷＤ２に対する送信とは、複数回ずつ交互に行われる。
【００７５】
また、セグメントスキャンでは、例えば、レンジゲートＲＧ１の反射波が複数回連続して受信され、レンジゲートＲＧ２の反射波が複数回連続して受信される。ここで、レンジゲートＲＧ１からの反射波の受信とレンジゲートＲＧ２からの反射波の受信とは、複数回ずつ交互に行われる。そして、例えば、ＲＧ１におけるドプラスペクトル画像用のデータとＲＧ２におけるドプラスペクトル画像用のデータとが、複数回の反射波データからなるセグメントの単位で交互にサンプリングされる。なお、セグメントスキャンでは、レンジゲートＲＧ１及びＲＧ２の位置に応じて、各レンジゲートからのデータ収集が最短で可能となるＰＲＦが設定される。
【００７６】
ここで、セグメントスキャンでは、レンジゲートＲＧ１に連続して超音波を送受信している間は、レンジゲートＲＧ２には超音波が送受信されないため、各レンジゲートにおけるドプラスペクトラム画像に周期的なデータの欠落が発生する。そこで、本実施形態では、各レンジゲートにおけるドプラスペクトラム画像に周期的なデータの欠落区間に補間データが挿入される。これにより、データの欠落によって生じる画像の劣化を抑えることができる。
【００７７】
図１４は、本実施形態に係るセグメントスキャンにおける処理の流れを示す図である。図１４に示すように、セグメントスキャンでは、ドプラ処理部１４が、レンジゲートＲＧ１及びＲＧ２それぞれからの反射データに対してウォールフィルター（Wall  Filter）、高速フーリエ変換（Fast  Fourier  Transformation：ＦＦＴ）を順に施す。このとき、ウォールフィルター及び高速フーリエ変換は時分割処理となる。なお、ドプラ処理部１４は、時分割処理でウォールフィルター及び高速フーリエ変換を行うのではなく、図１０に示した流れと同様に、レンジゲートＲＧ１及びＲＧ２それぞれからの反射波データに対して、別々にウォールフィルター及び高速フーリエ変換を施してもよい。
【００７８】
そして、ドプラ処理部１４は、高速フーリエ変換が施されたレンジゲートＲＧ１からのデータに対して、パラメータ同定処理、補間データ生成処理、後処理（Ｐｏｓｔ処理）を施すことで、レンジゲートＲＧ１におけるドプラスペクトラム画像に生じているデータの欠落区間に補間データを補填する。同様に、ドプラ処理部１４は、高速フーリエ変換が施されたレンジゲートＲＧ２からのデータに対しても、パラメータ同定処理、補間データ生成処理、後処理（Ｐｏｓｔ処理）を施すことで、レンジゲートＲＧ２におけるドプラスペクトラム画像に生じているデータの欠落区間に補間データを補填する。そして、画像生成部１５が、ドプラ処理部１４によって生成された各ドプラデータから、レンジゲートＲＧ１におけるドプラスペクトラム画像とレンジゲートＲＧ２におけるドプラスペクトル画像とをそれぞれ生成して表示部２に表示させる（Ｄｕａｌ−Ｄ表示）。
【００７９】
ここで、上述したインターリーブスキャン又はセグメントスキャンのいずれかを単独で行った場合には、従来のように、超音波の音速制約によって、良好なドプラスペクトル画像が得られない可能性がある。図１５は、超音波の音速制約を説明するための図である。図１５に示すように、超音波診断装置で用いられる超音波には、視野深度と、ＰＲＦと、ドプラの速度レンジとの間にトレードオフが生じる。
【００８０】
図１５に示す関係からも分かるように、ＰＲＦが小さくなると、視野深度は深くなるのに対して、ドプラの速度レンジは低くなる。インターリーブスキャンでは、レンジゲートの数に応じてＰＲＦが小さくなってしまうので、このような音速制約からドプラスペクトラム画像の速度レンジが小さくなってしまい、折り返し現象が発生しやすい。このため、インターリーブスキャンでは、例えば、深い位置に設定されたレンジゲートにおける速い血流を診断することが困難である。
【００８１】
これに対し、本実施形態では、スキャン切替部１７ｃが、距離判定部１７ｂによって距離の合計長が閾値未満であると判定された場合にはインターリーブスキャンを行い、距離判定部１７ｂにより距離の合計長が閾値以上であると判定された場合にはセグメントスキャンを行うようにスキャン方式を切り替える。すなわち、本実施形態では、深い位置に設定されたレンジゲートにおける血流を診断する場合には、自動的にスキャン方式がインターリーブスキャンからセグメントスキャンに切り替わる。したがって、本実施形態によれば、深い位置に設定されたレンジゲートにおける速い血流についても画質のよいドプラスペクトル画像を得ることができる。
【００８２】
図２の説明にもどって、計測値算出部１７ｄは、画像生成部１５により生成された第１のドプラスペクトラム画像により示される移動速度と第２のドプラスペクトラム画像により示される移動速度とから得られる計測値を算出する。
【００８３】
例えば、計測値算出部１７ｄは、診断部位が心臓の左室流入血流及び左室流出血流である場合には、左室流入血流の位置に設定されたレンジゲートＲＧ１におけるドプラスペクトル画像により示される血流速度と、左室流出血流の位置に設定されたレンジゲートＲＧ２におけるドプラスペクトル画像により示される血流速度とから、各種計測値を算出する。例えば、計測値算出部１７ｄは、Ｍｉｔｒａｌ系の計測値として、Ｅｖｅｌ、Ａｖｅｌ、Ｅ／Ａ（Ｅｖｅｌ／Ａｖｅｌ）、ＤｃＴなどの計測値を算出する。また、計測値算出部１７ｄは、Ａｏｒｔｉｃ系の計測値として、ＶＴＩ、ＶＰ、ＰＰＧ、ＭＰＧなどの計測値を算出する。
【００８４】
さらに、計測値算出部１７ｄは、左室流入血流及び左室流出血流に関する計測値として、ＩＲＴ（Isovolumetric  Relaxation  Time）、ＩＣＴ（Isovolumetric  Contraction  Time）、Ｔ．Ｉｎｄｅｘなどの計測値を算出する。図１７は、本実施形態に係る計測値算出部１７ｄによる計測値の算出の一例を示す図である。例えば、図１７に示すように、拡張期心室流入血流速度波形の終了から開始までの時間をａ、駆出時間（ＥＴ：Ejection  Time）をｂ、心電図のＲ波から心室流入血流速度波形の開始までの時間をｃ、心電図のＲ波から左室駆出血流速度波形の終了までの時間をｄとすると、ＩＲＴ、ＩＣＴは、それぞれ以下に示す式により算出される。
【００８５】
ＩＲＴ＝ｃ−ｄ
ＩＣＴ＝ａ−ｂ−ＩＲＴ
Ｔ．Ｉｎｄｅｘ＝（ａ−ｂ）／ｂ
【００８６】
また、計測値算出部１７ｄは、診断部位が心臓の左室流入血流ピーク速度及び僧帽弁輪移動速度である場合には、左室流入血流の位置に設定されたレンジゲートＲＧ１におけるドプラスペクトル画像により示される血流速度と、僧帽弁輪の位置に設定されたレンジゲートＲＧ２におけるドプラスペクトル画像により示される僧帽弁輪移動速度とから、各種計測値を算出する。例えば、計測値算出部１７ｄは、ＥＰＶ、ｅ’、ｅ’／Ｅなどの計測値を算出する。ここで、ＥＰＶは、左室流入血流速度波形におけるＥ波のピーク速度である。また、ｅ’は、僧帽弁輪移動速度のピーク値である。
【００８７】
また、計測値算出部１７ｄは、診断部位が頚動脈の総頚動脈及び内頚動脈である場合には、Ｂモード画像に基づいて各種計測値を算出する。図１９Ａ、１９Ｂ及び１９Ｃは、本実施形態に係る計測値算出部１７ｄによる計測値の算出の一例を示す図である。例えば、図１９Ａに示すように、計測値算出部１７ｄは、レンジゲートＲＧ１とレンジゲートＲＧ２との間の距離Ｌを算出する。また、例えば、図１９Ｂに示すように、計測値算出部１７ｄは、頚動脈の上下壁厚ｈ１及びｈ２、頚動脈の内径Ｄを算出する。
【００８８】
さらに、計測値算出部１７ｄは、診断部位が頚動脈の総頚動脈及び内頚動脈である場合には、総頚動脈の位置に設定されたレンジゲートＲＧ１におけるドプラスペクトル画像により示される血流速度と、内頚動脈の位置に設置されたレンジゲートＲＧ２おけるドプラスペクトル画像により示される血流速度とから、各種計測値を算出する。例えば、図１９Ｃに示すように、計測値算出部１７ｄは、ＣＣＡＶｅｌ、ＩＣＡＶｅｌ、Ｔ１などの計測値を算出する。ここで、Ｃｃａｖｅｌは、ＣＣＡの最高速度であり、Ｉｃａｖｅｌは、ＩＣＡの最高速度である。また、Ｔ１は、ＣＣＡピークとＩＣＡピークとの時間差である。また、計測値算出部１７ｄは、脈波速度Ｃから動脈硬化度Ｅを算出してもよい。例えば、動脈硬化度Ｅは、部位ごとにあらかじめ決められたプリセット値をρとすると、以下に示す式（１）により算出される。
【００８９】
【数１】
2012139489.tif 000003
【００９０】
図２の説明にもどって、計測値表示部１７ｅは、計測値算出部１７ｄにより算出された計測値を表示部２に表示させる。
【００９１】
図１６は、本実施形態に係る計測値表示部１７ｅによる計測値の表示の一例を示す図である。例えば、図１６に示すように、計測値表示部１７ｅは、診断部位が心臓の左室流入血流及び左室流出血流である場合には、Ｍｉｔｒａｌ系の計測値を表示するための表示領域１６１に、計測値算出部１７ｄにより算出されたＥｖｅｌ、Ａｖｅｌ、Ｅ／Ａ（Ｅｖｅｌ／Ａｖｅｌ）、ＤｃＴなどの計測値を表示する。また、計測値算出部１７ｄは、Ａｏｒｔｉｃ系の計測値を表示するための表示領域１６２に、計測値算出部１７ｄによって算出されたＶＴＩ、ＶＰ、ＰＰＧ、ＭＰＧなどの計測値を表示する。さらに、計測値表示部１７ｅは、左室流入血流及び左室流出血流に関する計測値を表示するための表示領域１６３に、計測値算出部１７ｄによって算出されたＩＲＴ、ＩＣＴ、Ｔ．Ｉｎｄｅｘなどの計測値を表示する。
【００９２】
図１８は、本実施形態に係る計測値表示部１７ｅによる計測値の表示の一例を示す図である。例えば、図１８に示すように、計測値表示部１７ｅは、診断部位が心臓の左室流入血流ピーク速度及び僧帽弁輪移動速度である場合には、左室流入血流ピーク速度及び僧帽弁輪移動速度に関する計測値を表示するための表示領域１８１に、計測値算出部１７ｄによって算出されたＥＰＶ、ｅ’、ｅ’／Ｅなどの計測値を出力する。
【００９３】
また、例えば、計測値表示部１７ｅは、診断部位が頚動脈の総頚動脈及び内頚動脈である場合には、図１９Ｃに示したように、ＣＣＡＶｅｌ、ＩＣＡＶｅｌなどの計測値を表示部２に表示する。
【００９４】
次に、本実施形態に係る超音波診断装置１００によるＢ／Ｄ同時スキャンの処理手順について説明する。図２０は、本実施形態に係る超音波診断装置１００によるＢ／Ｄ同時スキャンの処理手順を示すフローチャートである。
【００９５】
図２０に示すように、本実施形態に係る超音波診断装置１００では、制御部１７が、操作者からＢ／Ｄ同時スキャンの開始要求を受け付けたか否かを判定する（ステップＳ１０１）。そして、Ｂ／Ｄ同時スキャンの開始要求が受け付けられた場合には（ステップＳ１０１，Ｙｅｓ）、表示制御部１７ａが、画像生成部１５によって生成されたＢモード画像を表示部２に表示する（ステップＳ１０２）。
【００９６】
その後、表示制御部１７ａは、操作者によって診断用のアプリケーションが選択されるまで待機する（ステップＳ１０３，Ｎｏ）。そして、アプリケーションが選択された場合には（ステップＳ１０３，Ｙｅｓ）、表示制御部１７ａは、操作者によって診断部位が選択されるまで待機する（ステップＳ１０４，Ｎｏ）。
【００９７】
そして、診断部位が選択された場合には（ステップＳ１０４，Ｙｅｓ）、距離判定部１７ｂが、スキャン方式の切り換え判定に用いられる閾値を設定する（ステップＳ１０５）。その後、距離判定部１７ｂは、レンジゲートＲＧ１及びレンジゲートＲＧ２が設定されるまで待機する（ステップＳ１０６，Ｎｏ）。
【００９８】
そして、レンジゲートＲＧ１及びレンジゲートＲＧ２が設定された場合には（ステップＳ１０６，Ｙｅｓ）、距離判定部１７ｂは、超音波プローブ１からレンジゲートＲＧ１までの距離Ｒ１と超音波プローブ１からレンジゲートＲＧ２までの距離Ｒ２との合計長を算出する（ステップＳ１０７）。その後、距離判定部１７ｂは、算出した距離の合計長が閾値未満であるか否かを判定する（ステップＳ１０８）。
【００９９】
ここで、距離の合計長が閾値未満であった場合には（ステップＳ１０８，Ｙｅｓ）、スキャン切替部１７ｃが、スキャン方式をインターリーブスキャンに切り替える（ステップＳ１０９）。一方、距離の合計長が閾値以上であった場合には（ステップＳ１０８，Ｎｏ）、スキャン切替部１７ｃは、スキャン方式をセグメントスキャンに切り替える（ステップＳ１１０）。
【０１００】
続いて、表示制御部１７ａが、操作者によってデュアルドプラモードが選択されている場合には（ステップＳ１１１，Ｙｅｓ）、レンジゲートＲＧ１及びレンジゲートＲＧ２それぞれにおけるドプラスペクトル画像を表示部２に表示させる（ステップＳ１１２）。一方、表示制御部１７ａが、操作者によってデュアルドプラモードが選択されていない場合には（ステップＳ１１１，Ｎｏ）、レンジゲートＲＧ１又はレンジゲートＲＧ２おけるドプラスペクトル画像を表示部２に表示させる（ステップＳ１１３）。
【０１０１】
その後、操作者によってレンジゲートが変更された場合には（ステップＳ１１４，Ｙｅｓ）、制御部１７が、ステップＳ１０７に制御を戻す。こうして、制御部１７は、レンジゲートが変更されている間は、上述したスキャンの切り替えに関する処理を繰り返す。
【０１０２】
また、制御部１７は、レンジゲートが変更されずに（ステップＳ１１４，Ｎｏ）、操作者からＢ／Ｄ同時スキャンの終了要求も受け付けなかった場合には（ステップＳ１１５，Ｎｏ）、ステップＳ１０３に制御を戻す。こうして、操作者からＢ／Ｄ同時スキャンの終了要求を受け付けるまでの間は、制御部１７は、ステップＳ１０３〜Ｓ１１４の処理を繰り返す。そして、制御部１７は、操作者からＢ／Ｄ同時スキャンの終了要求を受け付けた場合には（ステップＳ１１５，Ｙｅｓ）、Ｂ／Ｄ同時スキャンに関する処理を終了する。
【０１０３】
次に、本実施形態に係る超音波診断装置１００による自動計測処理の処理手順について説明する。図２１は、本実施形態に係る超音波診断装置１００による自動計測処理の処理手順を示すフローチャートである。
【０１０４】
図２１に示すように、本実施形態に係る超音波診断装置１００では、制御部１７が、操作者からフリーズ要求を受け付けたか否かを判定する（ステップＳ２０１）。そして、フリーズ要求が受け付けられた場合には（ステップＳ２０１，Ｙｅｓ）、表示制御部１７ａが、Ｂモード画像及びドプラスペクトル画像をフリーズ（停止）する（ステップＳ２０２）。
【０１０５】
続いて、計測値算出部１７ｄが、画像生成部１５により生成された各ドプラスペクトラム画像により示される移動速度から得られる計測値を算出する（ステップＳ２０３）。そして、計測値表示部１７ｅが、計測値算出部１７ｄにより算出された計測値を表示部２に表示させる（ステップＳ２０４）。
【０１０６】
上述したように、本実施形態に係る超音波診断装置１００は、距離判定部１７ｂと、スキャン切替部１７ｃと、画像生成部１５と、表示部２とを有する。距離判定部１７ｂは、血流情報の観測部位として設定された少なくとも２つのレンジゲートについて、第１のレンジゲートから超音波プローブまでの距離と第２のレンジゲートから超音波プローブまでの距離との合計長が閾値未満であるか否かを判定する。スキャン切替部１７ｃは、距離の合計長が閾値未満であると判定された場合にはインターリーブスキャンを行い、距離の合計長が閾値以上であると判定された場合にはセグメントスキャンを行うようにスキャン方式を切り替える。画像生成部１５は、セグメントスキャン又はインターリーブスキャンにより受信された反射波データに基づいて、第１のレンジゲートにおける血流速度の経時的な変化を示す第１のドプラスペクトラム画像と第２のレンジゲートにおける血流速度の経時的な変化を示す第２のドプラスペクトル画像とをそれぞれ生成する。表示部２は、画像生成部１５により生成された第１のドプラスペクトラム画像及び第２のドプラスペクトル画像を表示する。
【０１０７】
このように、本実施形態に係る超音波診断装置１００は、複数箇所に設定されたレンジゲートそれぞれにおけるドプラスペクトル画像を表示する場合に、各レンジゲートの深さの合計に応じて、インターリーブスキャンとセグメントスキャンとを自動的に切り替える。これにより、深い位置に設定されたレンジゲートにおける血流を診断する場合には、自動的にスキャン方式がインターリーブスキャンからセグメントスキャンに切り替わる。したがって、本実施形態によれば、深い位置に設定されたレンジゲートにおける速い血流についても画質のよいドプラスペクトル画像を得ることができる。すなわち、本実施形態によれば、超音波の音速制約により生じるドプラスペクトル画像の画質の劣化を抑えることができる。
【０１０８】
なお、上記実施形態では、スキャン切替部１７ｃは、各レンジゲートから超音波プローブ１までの距離の合計長が閾値以上であると判定された場合に、スキャン方式をセグメントスキャンに切り替えることとした。例えば、これに加えて、スキャン切替部１７ｃは、画像生成部１５により生成された第１のドプラスペクトル画像又は第２のドプラスペクトル画像の速度レンジが所定の速度閾値を下回った場合にも、セグメントスキャンを行うようにスキャン方式を切り替えるようにしてもよい。これにより、超音波の音速制約により生じるドプラスペクトル画像の画質の劣化をより確実に抑えることができる。
【０１０９】
また、上記実施形態では、距離判定部１７ｂが、複数の観測部位のうち少なくとも２つの観測部位の走査線上における深さの合計と所定の閾値とを比較することとしたが、実施形態はこれに限られない。
【０１１０】
例えば、距離判定部１７ｂは、３つ以上の観測部位の走査線上における深さの合計と所定の閾値とを比較してもよい。この場合には、スキャン切替部１７ｃは、３つ以上の観測部位の走査線上における深さの合計が前記閾値を下回っていた場合に、インターリーブスキャンを行い、前記少なくとも２つの観測部位の走査線上における深さの合計が前記閾値を上回っていた場合に、セグメントスキャンを行うようにスキャン方式を切り替える。
【０１１１】
また、例えば、距離判定部１７ｂは、複数の観測部位のうち、いずれか１つの観測部の走査線上における深さを閾値と比較してもよい。例えば、距離判定部１７ｂは、複数の観測部位のうち、基準とすべき観測部位を指定する操作を操作者から受け付け、操作者によって指定された観測部位の走査線上における深さを閾値と比較する。この場合には、スキャン切替部１７ｃは、操作者によって指定された観測部位の深さが閾値を下回っていた場合に、インターリーブスキャンを行い、操作者によって指定された観測部位の深さが閾値を上回っていた場合に、セグメントスキャンを行うようにスキャン方式を切り替える。
【０１１２】
さらに、例えば、距離判定部１７ｂは、１つの観測部位を基準とするのではなく、複数の観測部位それぞれを閾値と比較してもよい。この場合には、スキャン切替部１７ｃは、複数の観測部位のうち、少なくとも１つの観測部位の深さが閾値を下回っていた場合に、インターリーブスキャンを行う。また、スキャン切替部１７ｃは、複数の観測部位のうち、少なくとも１つの観測部位の深さが閾値を上回っていた場合に、セグメントスキャンを行うようにスキャン方式を切り替える。
【０１１３】
また、上記実施形態では、スキャン切替部１７ｃが、閾値に基づいてスキャン方式を切り替えることとしたが、実施形態はこれに限られない。
【０１１４】
例えば、スキャン切替部１７ｃは、閾値に基づくスキャン方式の切り替えを行うとともに、さらに、ドプラスペクトル画像に折り返しが生じているか否かを検出し、折り返しが生じていることを検出した場合に、セグメントスキャンを行うようにスキャン方式を切り替えてもよい。
【０１１５】
この場合には、例えば、スキャン切替部１７ｃは、スキャン中に、所定の時間間隔で、ドプラスペクトル画像に折り返しが生じているか否かを検出する。ここで、折り返しを検出する方法としては、各種の方法を用いることができる。
【０１１６】
例えば、スキャン切替部１７ｃは、ドプラ処理部１４によって生成されたドプラデータに基づいて、血流速度の最大値の経時的な変化をトレースすることで、血流速度の最大値のトレース波形を検出する。このトレース波形は、ドプラスペクトラム画像の辺縁部をトレースした波形となる。さらに、スキャン切替部１７ｃは、検出したトレース波形に基づいて各速度の頻度を求め、速度の頻度分布を表すヒストグラムを作成する。そして、スキャン切替部１７ｃは、ヒストグラムから上限値ＵＬと下限値ＬＬを求め、絶対値｜ＵＬ−ＬＬ｜が第１の閾値より大きい値であり、かつ、｜ＵＬ｜又は｜ＬＬ｜のいずれかが第２の閾値より大きい場合に、ドプラスペクトル画像に折り返しが生じていると判定する。ここで、第１の閾値は、ノイズなどを判定するための値である。また、第２の閾値は、第１の閾値より大きい値であり、例えば、ナイキスト周波数（ＰＲＦの１／２）の値である。
【０１１７】
そして、スキャン切替部１７ｃは、折り返しが生じていることを検出した場合に、検出した時点でスキャン方式をセグメントスキャンに切り替える。なお、スキャン切替部１７ｃは、ドプラスペクトル画像に折り返しが生じていることを検出した場合に、ただちにスキャン方式を切り替えるのではなく、距離判定部１７ｂによって用いられる閾値を、折り返しを検出した時点での観測部位の深さ（少なくとも１つの観測部位の深さ又は複数の観測部位の深さの合計）より小さい値に変更してもよい。このように閾値を変更した場合には、距離判定部１７ｂが観測部位の深さと閾値とを比較した際に、観測部位の深さが閾値を上回ることになり、その結果、スキャン切替部１７ｃによってスキャン方式がセグメントスキャンに切り替えられる。
【０１１８】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【０１１９】
１００  超音波診断装置
２  表示部
１５  画像生成部
１７  制御部
１７ｂ  距離判定部
１７ｃ  スキャン切替部
１７ｆ  設定部
【図１】
2012139489.tif 000004
【図２】
2012139489.tif 000005
【図８】
2012139489.tif 000006
【図９】
2012139489.tif 000007
【図１０】
2012139489.tif 000008
【図１１】
2012139489.tif 000009
【図１２】
2012139489.tif 000010
【図１３】
2012139489.tif 000011
【図１４】
2012139489.tif 000012
【図１５】
2012139489.tif 000013
【図２０】
2012139489.tif 000014
【図２１】
2012139489.tif 000015
【図３Ａ】
2012139489.tif 000016
【図３Ｂ】
2012139489.tif 000017
【図４Ａ】
2012139489.tif 000018
【図４Ｂ】
2012139489.tif 000019
【図５Ａ】
2012139489.tif 000020
【図５Ｂ】
2012139489.tif 000021
【図６Ａ】
2012139489.tif 000022
【図６Ｂ】
2012139489.tif 000023
【図７Ａ】
2012139489.tif 000024
【図７Ｂ】
2012139489.tif 000025
【図１６】
2012139489.tif 000026
【図１７】
2012139489.tif 000027
【図１８】
2012139489.tif 000028
【図１９Ａ】
2012139489.tif 000029
【図１９Ｂ】
2012139489.tif 000030
【図１９Ｃ】
2012139489.tif 000031