(57)【要約】【課題】  四次元パラメトリックイメージングにより、ダイナミックな血流変化、微細な血管構築等の情報を高い視認性を持って簡単に観察可能な超音波診断装置等を提供する。【解決手段】  造影剤が投入された被検体内の三次元領域を所定期間に亘って超音波で走査し、前記三次元領域に関する超音波データを前記所定期間に亘って取得するデータ取得手段と、前記所定期間内の解析期間に亘る前記三次元領域に関する超音波データを用いて、前記解析期間内の各時相における第１のボリュームデータを生成すると共に、前記解析期間についての造影剤時間情報を示す第２のボリュームデータと、前記解析期間の前記三次元領域の各位置における造影剤特徴量を示す第３のボリュームデータと、を生成するボリュームデータ生成手段と、前記第２のボリュームデータと前記第３のボリュームデータとを用いて、投影画像を生成する画像生成手段と、を具備する。【選択図】図２
2011224354.tif 000002
【特許請求の範囲】
【請求項１】
造影剤が投入された被検体内の三次元領域を所定期間に亘って超音波で走査し、前記三次元領域に関する超音波データを前記所定期間に亘って取得するデータ取得ユニットと、
前記所定期間内の解析期間に亘る前記三次元領域に関する超音波データを用いて、前記解析期間内の各時相における第１のボリュームデータを生成すると共に、前記解析期間についての造影剤時間情報を示す第２のボリュームデータと、前記解析期間の前記三次元領域の各位置における造影剤特徴量を示す第３のボリュームデータと、を生成するボリュームデータ生成ユニットと、
前記第２のボリュームデータと前記第３のボリュームデータとを用いて、投影画像を生成する画像生成ユニットと、
を具備する超音波診断装置。
【請求項２】
前記画像生成ユニットは、前記三次元領域の各位置について、前記造影剤特徴量に応じて異なる輝度値を割り当てると共に前記造影剤時間情報に応じて異なる色調を割り当てることで、前記投影画像を生成する請求項１記載の超音波診断装置。
【請求項３】
前記画像生成ユニットは、
前記第２のボリュームデータを用いて、各投影方向につき、前記造影剤時間情報の分散値、前記造影剤時間情報の最大値、前記造影剤時間情報の最小値、前記造影剤時間情報の平均値のうちのいずれかを計算し、
前記計算された値に応じて、異なる色調を割り当てる請求項２記載の超音波診断装置。
【請求項４】
前記画像生成ユニットは、前記第３のボリュームデータを用いて、各投影方向につき、最大或いは最小の輝度値を有するボクセルを選択し、
前記選択されたボクセルに対応する前記第２のボリュームデータのボクセルを特定し、
前記特定されたボクセルの前記造影剤時間情報に従って、異なる色調を割り当てる請求項２記載の超音波診断装置。
【請求項５】
前記ボリュームデータ生成ユニットは、
前記解析期間内の各時相における第１のボリュームデータを用いて、前記三次元領域内の各位置の輝度値の時間的変化を計算し、
前記三次元領域内の各位置の信号値の時間的変化を用いて、前記三次元領域内の各位置における造影剤の到達時刻、輝度値が最大値になるピーク時刻、Wash-In時間、Wash-Out時間のいずれかを前記造影剤時間情報とする前記第２のボリュームデータを生成する請求項１乃至４のうちいずれか一項記載の超音波診断装置。
【請求項６】
前記投影画像を生成する際に用いる少なくとも一つの投影範囲を設定する設定ユニットをさらに具備し、
前記画像生成ユニットは、前記設定された少なくとも一つの投影範囲を用いて、少なくとも一つの前記投影画像を生成する請求項１乃至５のうちいずれか一項記載の超音波診断装置。
【請求項７】
前記投影画像が複数枚生成された場合には、当該複数の投影画像をマルチビュー表示する表示ユニットをさらに具備する請求項６記載の超音波診断装置。
【請求項８】
造影剤が投入された被検体内の三次元領域を所定期間に亘って超音波で走査することで取得された前記三次元領域に関する超音波データを記憶する記憶ユニットと、
前記所定期間内の解析期間に亘る前記三次元領域に関する超音波データを用いて、前記解析期間内の各時相における第１のボリュームデータを生成すると共に、前記解析期間についての造影剤時間情報を示す第２のボリュームデータと、前記解析期間の前記三次元領域の各位置における造影剤特徴量を示す第３のボリュームデータと、を生成するボリュームデータ生成ユニットと、
前記第２のボリュームデータと前記第３のボリュームデータとを用いて、投影画像を生成する画像生成ユニットと、
を具備する超音波画像処理装置。
【請求項９】
前記画像生成ユニットは、前記三次元領域の各位置について、前記造影剤特徴量に応じて異なる輝度値を割り当てると共に前記造影剤時間情報に応じて異なる色調を割り当てることで、前記投影画像を生成する請求項８記載の超音波画像処理装置。
【請求項１０】
前記画像生成ユニットは、
前記第２のボリュームデータを用いて、各投影方向につき、前記造影剤時間情報の分散値、前記造影剤時間情報の最大値、前記造影剤時間情報の最小値、前記造影剤時間情報の平均値のうちのいずれかを計算し、
前記計算された値に応じて、異なる色調を割り当てる請求項９記載の超音波画像処理装置。
【請求項１１】
前記画像生成ユニットは、前記第３のボリュームデータを用いて、各投影方向につき、、最大或いは最小の輝度値を有するボクセルを選択し、
前記選択されたボクセルに対応する前記第２のボリュームデータのボクセルを特定し、
前記特定されたボクセルの前記造影剤時間情報に従って、異なる色調を割り当てる請求項９記載の超音波画像処理装置。
【請求項１２】
前記ボリュームデータ生成ユニットは、
前記解析期間内の各時相における第１のボリュームデータを用いて、前記三次元領域内の各位置の輝度値の時間的変化を計算し、
前記三次元領域内の各位置の信号値の時間的変化を用いて、前記三次元領域内の各位置における造影剤の到達時刻、輝度値が最大値になるピーク時刻、Wash-In時間、Wash-Out時間のいずれかを前記造影剤時間情報とする前記第２のボリュームデータを生成する請求項８乃至１１のうちいずれか一項記載の超音波画像処理装置。
【請求項１３】
前記投影画像を生成する際に用いる少なくとも一つの投影範囲を設定する設定ユニットをさらに具備し、
前記画像生成ユニットは、前記設定された少なくとも一つの投影範囲を用いて、少なくとも一つの前記投影画像を生成する請求項８乃至１２のうちいずれか一項記載の超音波画像処理装置。
【請求項１４】
前記投影画像が複数枚生成された場合には、当該複数の投影画像をマルチビュー表示する表示ユニットをさらに具備する請求項１３記載の超音波画像処理装置。
【請求項１５】
造影剤が投入された被検体内の三次元領域を所定期間に亘って撮像し、前記三次元領域に関する画像データを前記所定期間に亘って取得するデータ取得ユニットと、
前記所定期間内の解析期間に亘る前記三次元領域に関する画像データを用いて、前記解析期間内の各時相における第１のボリュームデータを生成すると共に、前記解析期間についての造影剤時間情報を示す第２のボリュームデータと、前記解析期間の前記三次元領域の各位置における造影剤特徴量を示す第３のボリュームデータと、を生成するボリュームデータ生成ユニットと、
前記第２のボリュームデータと前記第３のボリュームデータとを用いて、投影画像を生成する画像生成ユニットと、
を具備する医用画像診断装置。
【請求項１６】
前記画像生成ユニットは、前記三次元領域の各位置について、前記造影剤特徴量に応じて異なる輝度値を割り当てると共に前記造影剤時間情報に応じて異なる色調を割り当てることで、前記投影画像を生成する請求項１５記載の医用画像診断装置。
【請求項１７】
前記画像生成ユニットは、
前記第２のボリュームデータを用いて、各投影方向につき、前記造影剤時間情報の分散値、前記造影剤時間情報の最大値、前記造影剤時間情報の最小値、前記造影剤時間情報の平均値のうちのいずれかを計算し、
前記計算された値に応じて、異なる色調を割り当てる請求項１６記載の医用画像診断装置。
【請求項１８】
前記画像生成ユニットは、前記第３のボリュームデータを用いて、各投影方向につき、、最大或いは最小の輝度値を有するボクセルを選択し、
前記選択されたボクセルに対応する前記第２のボリュームデータのボクセルを特定し、
前記特定されたボクセルの前記造影剤時間情報に従って、異なる色調を割り当てる請求項１６記載の医用画像診断装置。
【請求項１９】
前記ボリュームデータ生成ユニットは、
前記解析期間内の各時相における第１のボリュームデータを用いて、前記三次元領域内の各位置の輝度値の時間的変化を計算し、
前記三次元領域内の各位置の信号値の時間的変化を用いて、前記三次元領域内の各位置における造影剤の到達時刻、輝度値が最大値になるピーク時刻、Wash-In時間、Wash-Out時間のいずれかを前記造影剤時間情報とする前記第２のボリュームデータを生成する請求項１５乃至１８のうちいずれか一項記載の医用画像診断装置。
【請求項２０】
前記投影画像を生成する際に用いる少なくとも一つの投影範囲を設定する設定ユニットをさらに具備し、
前記画像生成ユニットは、前記設定された少なくとも一つの投影範囲を用いて、少なくとも一つの前記投影画像を生成する請求項１５乃至１９のうちいずれか一項記載の医用画像診断装置。
【請求項２１】
前記投影画像が複数枚生成された場合には、当該複数の投影画像をマルチビュー表示する表示ユニットをさらに具備する請求項２０記載の医用画像診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波造影剤であるマイクロバブルを用いることで、画像診断の対象とされる臓器等の血流動態の時間的な変化を三次元的に映像化することが可能な超音波診断装置等であって、ダイナミックな血流変化、微細な血管構築等の情報を高い視認性を持って簡単に観察可能な超音波診断装置、超音波画像処理装置及び医用画像診断装置に関する。
【背景技術】
【０００２】
超音波画像診断は、超音波プローブを体表から当てるだけの簡単な操作で心臓の拍動や胎児の動きの様子がリアルタイム表示で得られ、かつ安全性が高いため繰り返して検査が行えるほか、システムの規模がX線、CT、MRIなど他の診断機器に比べて小さく、ベッドサイドへ移動していっての検査も容易に行えるなど簡便である。また、超音波診断はＸ線などのように被曝の影響がなく、産科や在宅医療等においても使用することができる。特に、近年では、揺動プローブ（一次元に配列された超音波振動子列を機械的に揺動させることで、三次元領域を超音波走査可能なプローブ）或いは二次元アレイプローブ（マトリックス状に配列された超音波振動子列により、三次元領域を超音波走査可能なプローブ）を用いて、リアルタイム三次元走査（「四次元走査」ともいう）も可能になっている。
【０００３】
また、近年、静脈投与型の超音波造影剤が製品化されている。超音波画像診断において、この様な造影剤を用いた造影イメージングを行うことで、種々の情報を取得することが可能となっている。例えば、原発性肝癌の場合、動脈相の血流情報や門脈相の血流情報を相補的に用いることで良性悪性の鑑別診断などを行うことができる。また、RFA治療後に造影イメージングを行った場合、焼妁された領域は染影しないことから、効果判定にも応用されている。さらに、心臓および肝臓などの検査で静脈から超音波造影剤を注入して血流信号を増強し、血流動態の評価を行うための撮像方法も提案されている。再還流（replenishment）法、マイクロフローイメージング（Micro Flow Imaging）は、上記血流動態の評価を行うための撮像方法の代表例である。
【０００４】
再還流法は、造影剤気泡が崩壊するという特徴を生かし、（ａ）低音圧照射下でスキャン断面に充満していく気泡の動態を観察し、（ｂ）照射音圧を高音圧に切り替えて、断面内（厳密には照射体積内）の気泡を崩壊させ、（ｃ）再び断面内に流入していく気泡の様子を観察する、という手法である（例えば、特許文献１参照）。
【０００５】
一般に、造影イメージングの特徴は、カラードプラ等では映像化できない微小血流を映像化できる点にある。一方、微小血流内には、少数のバブルしか存在しない。このため、非定常的に染影する。マイクロフローイメージングは、この様な造影イメージング及び造影剤の特徴を利用して、非定常的に映像化されるバブルを時間方向に重畳表示することで、微小血流構造をより明瞭に映像化する手法である（例えば、特許文献２参照）。マイクロフローイメージングを実行する場合、一定時間、患者は息止め、技師等の検者はプローブを固定しなければならないという制約がある。この制約を緩和する技術として、対象領域の動きを画像処理にて補正する技術も提案されている。（例えば、特許文献２参照）。
【０００６】
これらの造影イメージングにおいては、例えば、染影の開始時相（すなわちバブルの到達時間に相当）、染影がピークを迎える時相（ピーク時刻）、バブルの到達時刻−ピーク時刻で定義されるWash-In時間、バブルが静脈に流出する時刻−ピーク時刻で定義されるWash-Out時間等のダイナミック情報を、空間的に（すなわち、各位置において）取得する。これらのダイナミック情報は、動画として観察することが好ましい場合がある。例えば肝臓の場合、腫瘍は動脈血で栄養されており、正常実質は門脈血で栄養されている。このため、動画での観察を行えば、腫瘍部への造影剤到達時刻は周囲に比べて早くなり、腫瘍のほうが早く染影し始める様子を視認することができる。
【０００７】
その一方で、ダイナミック情報の動画観察は、各位置で時間変化を頭の中で覚える必要があるため、必ずしも客観的であるとは言えない。また、観察に或る程度の時間が必要となってしまうという問題もある。これらの問題の解決するものとして、一枚の二次元画像によってダイナミック情報を観察可能とする二次元パラメトリックイメージング技術が提案されている。この技術では、例えば、各ピクセルにおける輝度の時間変化の最大値をピーク時刻とする、或いはノイズ等の影響を受けないような閾値を設け、その閾値を最初に越えた時刻を到達時刻とする等の計算を行い、時間情報を求める。そして、求めた時間情報の値に応じて、位置毎に色調を変えた二次元画像を表示するというものである。上述のMicro Flow Imagingや動き補正の技術とも組み合わせて使うと効果的であり、今後臨床の現場での利用が期待されている。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開平１１−１５５８５８号公報
【特許文献２】特開２００４−３２１６８８号公報
【特許文献３】特願２００６−２６６１
【特許文献４】特許第２７１４３２９号公報
【特許文献５】特許第３４９５７１０号公報
【特許文献６】特開２００９−２３３４０８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００９】
しかしながら、二次元パラメトリックイメージングでは、二次元画像データ（断層データ）を用いているため、血流の空間的連続性の表現に限界がある。例えば、断層像に直交する血管は、実質的に点として表示されてしまう等、奥行き方向の構造を持つ血管等を認識することは困難である。
【００１０】
一方、血流動態（血流の時間変化）を観察する場合、時相の異なる複数のボリュームデータを観察するとすれば、従来の二次元画像（断層像）の診断時に比べ、観察範囲が1次元増えるため、観察時間が増加し、さらに頭で時間変化情報の構築（記憶）も曖昧になるという問題がある。
【００１１】
従って、時系列のボリュームデータを用いたパラメトリックイメージング（四次元パラメトリックイメージング）を行うことが効果的と思われる。しかしながら、到達時間等の時間情報の三次元的映像化法については、これまでに臨床の点おいて効果的な方法が提案されていない。これは、二次元パラメトリックイメージングの手法を単純に三次元的に行うだけでは、奥行き方向に関する時間情報を適切に扱うことができないからである。特に、近年、四次元走査が可能な超音波診断装置においても、造影イメージング機能が導入されてきている。従って、四次元パラメトリックイメージングの手法を確立し、ダイナミックな血流動態の三次元的な観察を可能にすることは、臨床的価値も高く、強く望まれるところである。
【００１２】
本発明は、上記事情を鑑みてなされたもので、四次元パラメトリックイメージングにより、ダイナミックな血流変化、微細な血管構築等の情報を高い視認性を持って簡単に観察可能な超音波診断装置、超音波画像処理装置及び医用画像診断装置を提供することを目的としている。
【課題を解決するための手段】
【００１３】
本発明は、上記目的を達成するため、次のような手段を講じている。
【００１４】
一実施形態に係る超音波診断装置は、造影剤が投入された被検体内の三次元領域を所定期間に亘って超音波で走査し、前記三次元領域に関する超音波データを前記所定期間に亘って取得するデータ取得ユニットと、前記所定期間内の解析期間に亘る前記三次元領域に関する超音波データを用いて、前記解析期間内の各時相における第１のボリュームデータを生成すると共に、前記解析期間についての造影剤時間情報を示す第２のボリュームデータと、前記解析期間の前記三次元領域の各位置における造影剤特徴量を示す第３のボリュームデータと、を生成するボリュームデータ生成ユニットと、前記第２のボリュームデータと前記第３のボリュームデータとを用いて、投影画像を生成する画像生成ユニットと、を具備するものである。
【図面の簡単な説明】
【００１５】
【図１】図１は、本実施形態に係る超音波診断装置１のブロック構成図を示している。
【図２】図２は、本四次元パラメトリックイメージング機能に従う処理（四次元パラメトリックイメージング処理）の流れを示したフローチャートである。
【図３】図３は、所定期間の各時相ｋ（ｋ＝１，２，・・・、ｎ）に対応する第１のボリュームデータを概念的に示した図である。
【図４】図４は、あるボクセルについての解析期間に亘るTICの一例を示した図である。
【図５】図５は、解析期間の各時相に対応する複数の第1のボリュームデータを用いた第２のボリュームデータ生成処理を説明するための図である。
【図６】図６は、解析期間の各時相に対応する複数の第1のボリュームデータを用いた第３のボリュームデータ生成処理を説明するための図である。
【図７】図７は、第２、第３のボリュームデータを用いた投影処理を説明するための図である。
【図８】図８は、第２、第３のボリュームデータを用いた投影画像の生成に用いるカラーマップの一例を示した図である。
【図９】図９は、造影剤時間情報の奥行き方向に関する割り当て処理を説明するための図である。
【図１０】図１０は、造影剤時間情報の奥行き方向に関する割り当て処理を説明するための図である。
【図１１】図１１は、造影剤時間情報の奥行き方向に関する割り当て処理を説明するための図である。
【図１２】図１２は、第１、第２又は第３のボリュームデータを分割して複数の投影領域を生成する場合の一例を説明するための図である。
【図１３】図１３は、第１、第２又は第３のボリュームデータを分割して複数の投影領域を生成する場合の他の例を説明するための図である。
【図１４】図１４は、第１、第２又は第３のボリュームデータを分割して複数の投影領域を生成する場合の他の例を説明するための図である。
【発明を実施するための形態】
【００１６】
以下、本発明の実施形態を図面に従って説明する。なお、以下の説明において、略同一の機能及び構成を有する構成要素については、同一符号を付し、重複説明は必要な場合にのみ行う。
【００１７】
図１は、本実施形態に係る超音波診断装置１のブロック構成図を示している。同図に示すように、本超音波診断装置１０は、超音波プローブ１２、入力装置１３、モニター１４、超音波送信ユニット２１、超音波受信ユニット２２、Ｂモード処理ユニット２３、ドプラ処理ユニット２４、画像生成ユニット２５、画像メモリ２６、画像合成ユニット２７、制御プロセッサ（ＣＰＵ）２８、記憶ユニット２９、インターフェースユニット３０、ソフトウェア格納ユニット３１を具備している。装置本体１１に内蔵される超音波送信ユニット２１、受信ユニット２２等は、集積回路などのハードウェアで構成されることもあるが、ソフトウェア的にモジュール化されたソフトウェアプログラムである場合もある。以下、個々の構成要素の機能について説明する。
【００１８】
超音波プローブ１２は、超音波送信ユニット２１からの駆動信号に基づき超音波を発生し、被検体からの反射波を電気信号に変換する複数の圧電振動子、当該圧電振動子に設けられる整合層、当該圧電振動子から後方への超音波の伝播を防止するバッキング材等を有している。当該超音波プローブ１２から被検体Ｐに超音波が送信されると、当該送信超音波は、体内組織の音響インピーダンスの不連続面で次々と反射され、エコー信号として超音波プローブ１２に受信される。このエコー信号の振幅は、反射することになった反射することになった不連続面における音響インピーダンスの差に依存する。また、送信された超音波パルスが、移動している血流や心臓壁等の表面で反射された場合のエコーは、ドプラ効果により移動体の超音波送信方向の速度成分を依存して、周波数偏移を受ける。
【００１９】
なお、本超音波装置が具備する超音波プローブ１２は、被検体の三次元領域を超音波走査可能なものである。そのため、超音波プローブ１２は、振動子をその配列方向の直交方向に沿って機械的に揺動させ、三次元領域を超音波走査する構成、又は二次元的に配列された二次元振動素子を用いて電気的制御により三次元領域を超音波走査する構成等を有する。前者の構成を採用する場合、被検体の三次元的走査は揺動回路（揺動機構）によって行われるため、検査者はプローブ本体を被検体に接触させるだけで、自動的に複数の二次元断層像を取得することができる。制御された揺動速度から断面間の正確な距離も検知できる。また、後者の構成を採用する場合には、原理的には、従来の二次元断層像を取得するのと同じ時間で、三次元領域を超音波走査することができる。
【００２０】
入力装置１３は、装置本体１１に接続され、オペレータからの各種指示、条件、関心領域（ＲＯＩ）の設定指示、種々の画質条件設定指示等を装置本体１１にとりこむためのトラックボール、各種スイッチ、ボタン、マウス、キーボード等を有している。
【００２１】
モニター１４は、画像合成ユニット２７らのビデオ信号に基づいて、生体内の形態学的情報や、血流情報を画像として表示する。
【００２２】
超音波送信ユニット２１は、図示しないトリガ発生回路、遅延回路およびパルサ回路等を有している。パルサ回路では、所定のレート周波数ｆｒ  Ｈｚ（周期；１／ｆｒ秒）で、送信超音波を形成するためのレートパルスが繰り返し発生される。また、遅延回路では、チャンネル毎に超音波をビーム状に集束し且つ送信指向性を決定するのに必要な遅延時間が、各レートパルスに与えられる。この遅延情報を変化させることで、プローブ振動子面からの送信方向を任意に調整することが可能となる。トリガ発生回路は、このレートパルスに基づくタイミングで、プローブ１２に駆動パルスを印加する。
【００２３】
超音波受信ユニット２２は、図示していないアンプ回路、Ａ／Ｄ変換器、加算器等を有している。アンプ回路では、プローブ１２を介して取り込まれたエコー信号をチャンネル毎に増幅する。Ａ／Ｄ変換器では、増幅されたエコー信号に対し受信指向性を決定するのに必要な遅延時間を与え、その後加算器において加算処理を行う。この加算により、エコー信号の受信指向性に応じた方向からの反射成分が強調され、受信指向性と送信指向性とにより超音波送受信の総合的なビームが形成される。
【００２４】
Ｂモード処理ユニット２３は、超音波受信ユニット２２からエコー信号を受け取り、対数増幅、包絡線検波処理などを施し、信号強度が輝度の明るさで表現されるデータを生成する。このとき検波周波数を変化させることで、映像化する周波数帯域を変えることができる。また一つの受信データに対して、２つの検波周波数による検波処理を並列に行うことも可能な構成である。それを利用すると一つの受信信号からバブル像と組織像を生成することも可能となる。なお、バブルを映像化する場合は、セカンドハーモニック（２次高調波）を主に映像化することが多い。本Ｂモード処理ユニット２３において処理を受けたデータは、画像生成ユニット２５に出力され、反射波の強度を輝度にて表したＢモード画像として再構成される。
【００２５】
ドプラ処理ユニット２４は、送受信ユニット２２から受け取ったエコー信号から速度情報を周波数解析し、ドプラ効果による血流や組織、造影剤エコー成分を抽出し、平均速度、分散、パワー等の血流情報を多点について求める。得られた血流情報は画像生成ユニット２４に送られ、平均速度画像、分散画像、パワー画像、これらの組み合わせ画像としてモニター１４にカラー表示することができる。
【００２６】
画像生成ユニット２５は、超音波スキャンの走査線信号列を、テレビなどに代表される一般的なビデオフォーマットの走査線信号列に変換し、表示画像としての超音波診断画像を生成する。画像生成ユニット２４は、画像データを格納する記憶メモリを搭載しており、例えば診断の後に操作者が検査中に記録された画像を呼び出すことが可能となっている。また、画像生成ユニット２５は、画像処理装置としての機能を持ち、例えばボリュームデータを構築する場合、三次元領域或いは連続する二次元領域を超音波走査することで得られる走査線信号列を空間的に配置し、座標変換、補間処理等を必要に応じて実行することで、ボリュームデータを構築する。画像生成ユニット２５は、得られたボリュームデータを用いてボリュームレンダリングや、空間的ＭＩＰ等の処理、ボリュームデータ内の任意の断層像を切り出しによるＭＰＲ処理、マイクロフローイメージング処理、動き補正処理等を実行し、所定の三次元画像を生成する。なお、本画像生成ユニット２５における各種画像処理法等は、ソフトウェア的手法或いはハードウェア的手法のいずれであっても構わない。
【００２７】
さらに、画像生成ユニット２５は、制御プロセッサ２８の制御のもと、後述する四次元パラメトリックイメージング機能に従う所定の処理を実行する。
【００２８】
画像メモリ２６は、複数フレーム或いは複数ボリューム分の超音波データを一時的に記憶する。
【００２９】
画像合成ユニット２７は、画像生成ユニット２５から受け取った画像を種々のパラメータの文字情報や目盛等と共に合成し、ビデオ信号としてモニター１４に出力する。
【００３０】
制御プロセッサ（ＣＰＵ）２８は、情報処理装置（計算機）としての機能を持ち、本超音波診断装置本体の動作を制御する。制御プロセッサ２８は、記憶ユニット２９から各種画像処理法を実現するためのプログラム、後述する四次元パラメトリックイメージング機能を実現するためのプログラムを読み出してソフトウェア格納部３１上に展開し、各種処理に関する演算・制御等を実行する。
【００３１】
記憶ユニット２９は、各種スキャンシーケンスを実行するためのプログラム、後述する四次元パラメトリックイメージング機能を実現するための専用プログラム、画像生成、表示処理を実行するための制御プログラム、診断情報（患者ＩＤ、医師の所見等）、診断プロトコル、送受信条件、ボディマーク生成プログラムその他のデータ群が保管されている。また、必要に応じて、画像メモリ２６中の画像の保管などにも使用される。記憶ユニット２９のデータは、インターフェースユニット３０を経由して外部周辺装置へ転送することも可能となっている。
【００３２】
インターフェースユニット３０は、入力装置１３、ネットワーク、新たな外部記憶装置（図示せず）に関するインターフェースである。当該装置によって得られた超音波画像等のデータや解析結果等は、インターフェースユニット３０よって、ネットワークを介して他の装置に転送可能である。
【００３３】
（四次元パラメトリックイメージング機能）
次に、本超音波診断装置１が有する四次元パラメトリックイメージング機能について説明する。この機能は、造影イメージングにおいて、所定期間に亘る前記三次元領域に関する超音波データを用いて、所定期間内の各時相における造影剤（バブル）の空間密度を信号値（輝度値）として反映する第１のボリュームデータを生成する。また、当該所定期間内の各時相毎の第１のボリュームデータを用いて、三次元領域の各位置における造影剤時間情報を示す第２のボリュームデータと、三次元領域の各位置における造影剤特徴量を示す第３のボリュームデータと、を生成する。そして、生成された第２のボリュームデータ及び第３のボリュームデータを用いて、三次元領域の各位置における造影剤特徴量が輝度値により現され、且つ造影剤時間情報が色調で現される三次元画像を生成し、所定の形態で表示するものである。
【００３４】
なお、以下の説明においては、説明を具体的にするために、超音波診断装置を用いて四次元パラメトリックイメージング機能を実現する場合を例とする。しかしながら、当該例に拘泥されず、例えば、例えば医用ワークステーションによる超音波画像処理装置（或いは、ビューワ）を用いて、四次元走査によって事前に取得された超音波データを用いても、行うことができる。係る場合には、例えば、図１の点線内の構成を採用することができる。さらに、造影剤を用いる他の医用画像診断装置（Ｘ線コンピュータ断層撮像装置、磁気共鳴イメージング装置、Ｘ線診断装置等）においても、例えば図１の点線内の構成と実質的に同じ構成を実行することで、四次元パラメトリックイメージング機能を実現することが可能である。
【００３５】
図２は、本四次元パラメトリックイメージング機能に従う処理（四次元パラメトリックイメージング処理）の流れを示したフローチャートである。以下、当該フローチャートに示す各ステップにおいて実行される処理の内容について説明する。
【００３６】
［患者情報の入力、送受信条件、スキャンシーケンス等の選択：ステップＳ１］
操作ユニット３３を介して患者情報の入力、送受信条件（被走査領域の大きさを決めるための画角、焦点位置、送信電圧等）、被検体の三次元領域を所定期間に亘って超音波走査するためのスキャンシーケンス等の選択が実行される（ステップＳ１）。入力、選択された各種情報・条件等は、自動的に記憶装置２９に記憶される。
【００３７】
［所定期間に亘るボリュームデータの収集：ステップＳ２］
次に、送受信制御ユニット３１は、被検者の所定部位（診断対象とする血管等）を含む三次元領域を被走査領域として、リアルタイム三次元超音波走査（四次元走査）を所定期間に亘って実行し、当該所定期間内の各時相に対応する三次元超音波データを取得する（ステップＳ２）。
【００３８】
［造影剤の空間密度を表す第１のボリュームデータの生成：ステップＳ３］
次に、画像生成ユニット２５は、取得された三次元超音波データに対して、造影剤からの信号成分に対応するハーモニック成分の抽出処理、座標変換、補間処理等を施し、例えば図３に示すような各時相ｋ（ｋ＝１，２，・・・、ｎ）に対応する第１のボリュームデータを生成する（ステップＳ３）。この第１のボリュームデータは、各時相における造影剤の空間密度（濃度）を信号値（輝度値）で表すものとなっている。従って、時系列の当該第１のボリュームデータを解析することで、被走査領域である三次元領域内の各位置における造影剤密度の時間的変化を把握することが可能である。
【００３９】
［造影剤に関する時間情報を示す第２のボリュームデータの生成：ステップＳ４］
次に、画像生成ユニット２５は、所定期間に亘る第１のボリュームデータのうち、解析期間（所定期間のうち、四次元パラメトリックイメージングによる解析対象とする期間）に対応する時相毎の第１のボリュームデータを用いて、造影剤に関する時間情報を示す第２のボリュームデータを生成する（ステップＳ４）。ここで、造影剤に関する時間情報とは、各位置における造影剤の到達時刻、信号値（又は輝度値）のピーク時刻、Wash-In時間、Wash-Out時間等である。
【００４０】
なお、解析期間の開始時刻及び終了時刻は、入力装置１３からの操作により任意に設定することが可能である。また、例えば、造影剤タイマーのON時刻、バブルを破壊するバーストの終了時刻、ある対象にバブルが入ってきた（到達した）時刻等を検出し、装置によって自動的に解析期間を開始等するようにしてもよい。
【００４１】
また、第２のボリュームデータは、解析期間に対応する複数の第１のボリュームデータと（第1の各ボリュームの時相と一致して）同じ数だけ作成することを前提とする。しかしながら、これに限定されることなく、例えば第２のボリュームデータを単一のボリュームデータとすることも可能である。以下、本ステップにおいて生成される第２のボリュームを造影剤時間情報の種別毎に説明する。
【００４２】
（造影剤時間情報を造影剤の到達時刻とする場合）
まず、画像生成ユニット２５は、解析期間の各時相に対応する複数の第1のボリュームデータを用いて、各ボクセルについての解析期間におけるTIC（Time Intensity Curve）を生成する。
【００４３】
図４は、あるボクセルについての解析期間に亘るTICの一例を示した図である。画像生成ユニット２５は、生成された各ボクセルについての解析期間におけるTICを用いて各ボクセルの輝度の時間変化を観察し、図４に示すような予め設定された輝度の閾値を始めて越えた時刻Ｔａｒｒｉｖｅｄを、ボリュームレートとボリューム数とを用いてボクセル毎に計算する。この計算された時刻をボクセルの値とすることで、第２のボリュームデータを生成することができる。
【００４４】
また、解析期間の最初の時相に対応する第１のボリュームデータにおいて、輝度値が最初から閾値を越えているボクセルについては、ボクセル値をゼロまたは値なしとしてもよい。なぜなら、この様なボクセルの輝度値は、バブルからの信号に起因するものではなく、組織からの信号に起因するものが多いからである。
【００４５】
この様に、造影剤時間情報が造影剤の到達時刻である場合には、解析期間に対応する複数の第１のボリュームデータと同じ数だけ作成してもよい。しかしながら、例えば解析期間内の最後の時相に対応する第２のボリュームデータが一つあれば、各ボクセルについての解析期間における造影剤到達時刻を把握することができる。
【００４６】
（造影剤時間情報を輝度値がピークになるピーク時刻とする場合）
まず、画像生成ユニット２５は、解析期間の各時相に対応する複数の第1のボリュームデータを用いて、各ボクセルについての解析期間におけるTIC（Time Intensity Curve）を生成する。次に、画像生成ユニット２５は、生成された各ボクセルについての解析期間におけるTICを用いて各ボクセルの輝度の時間変化を観察し、図４に示すように輝度値が最大（ピーク）となる時刻Ｔｐｅａｋをボクセル毎に計算する。この計算された時刻をボクセルの値とすることで、第２のボリュームデータを生成することができる。
【００４７】
この様に、造影剤に関する時間情報が、輝度値がピークになるピーク時刻である場合においても、解析期間に対応する複数の第１のボリュームデータと同じ数だけ作成してもよい。しかしながら、到達時刻は一度求まると更新されることはないのに対し、ピーク時間の場合は、例えば図５に示すように時系列で隣り合うボリューム間で輝度値の大小比較演算ｆ（ｋ−１，ｋ）（ただし、ｋ＝１，２，・・・、ｎ）を実行し、輝度が過去の最大値を超えるたびに、時間情報も更新されることになる。厳密に言えば、輝度がピークを迎え、Wash-Outする時間まで観察するならば、各ボクセルでピーク時間は一つしかないことになる。従って、例えば解析期間内の最後の時相に対応する第２のボリュームデータが一つあれば、各ボクセルについての解析期間における造影剤の到達時刻を把握することができる。
【００４８】
（造影剤時間情報を造影剤のWash-In時間とする場合）
Wash-Inの時間を、図４に示すように、例えばピーク時刻−到達時刻として定義する。係る場合には、解析期間についてのTICを用いて各ボクセルの輝度の時間変化を観察し、輝度値がある閾値を超えた時刻（到達時刻）と最大（ピーク）となる時刻とをボクセル毎に特定し、各ボクセルについてのWash-In時間を計算することができる。この計算された時間（期間）をボクセルの値とすることで、造影剤に関する時間情報を造影剤のWash-In時間とする場合の第２のボリュームデータを、例えば単一のボリュームとして生成することができる。
【００４９】
（造影剤に関する時間情報を造影剤のWash-Out時間とする場合）
Wash-Outの時間を、図４に示すように、例えばピークを一度迎えた時刻（ピーク時刻）からTICが所定の閾値以下になる時刻までの期間として定義する。なお、TICは、フィッティング（Fitting）又は平滑化（Smoothing）されていることが好ましい。係る場合には、解析期間におけるTICを用いて各ボクセルの輝度の時間変化を観察し、輝度値が最大（ピーク）となる時刻から所定の閾値以下になる時刻までの時間をボクセル毎に計算する。この計算された時間（期間）をボクセルの値とすることで、造影剤に関する時間情報を造影剤のWash-Out時間とする場合の第２のボリュームデータを、例えば単一のボリュームとして生成することができる。
【００５０】
なお、以上述べた各種造影剤に関する時間情報に関する第２のボリュームデータの生成においては、同じボクセルが同位置であるという前提の処理である。しかしながら、実際にはいくらプローブをしっかり保持しようとしても、手が動いてしまったり、患者の呼吸で動いてしまったり、完全に固定しておくことは困難である。そこで、例えば特開２００７−３３０７６４に記載の動き補正技術を用いて装置が連続するボリューム間の位置ずれを補正することは、実際の検査でこの提案手法を使う場合に、非常に有用である。
【００５１】
［各位置の造影剤特徴量を示す第３のボリュームデータの生成：ステップＳ５］
次に、画像生成ユニット２５は、解析期間に対応する時相毎の第１のボリュームデータを用いて、解析期間における造影剤特徴量を被走査領域の各位置毎に示す第３のボリュームデータを生成する（ステップＳ５）。ここで、解析期間における造影剤特徴量とは、解析期間での各位置についての信号値（輝度値）の最大値、最小値、平均値等である。
【００５２】
例えば、血管形態観察のためにバブルをトレースするために、解析期間における造影剤特徴量として、信号値（輝度値）の最大値を用いる場合を想定する。係る場合には、画像生成ユニット２５は、解析期間に対応する時相毎の第１のボリュームデータを用いて、空間的位置が対応する各ボクセルについて、時間方向に最大値を保持（temporal maxhold）を実行することで（例えば図６に示すように時系列で隣り合うボリューム間で輝度値の最大値（大小）比較演算Ｍａｘ（ｋ−１，ｋ）（ただし、ｋ＝１，２，・・・、ｎ）を逐次実行し、輝度が過去の最大値を超えるたびに、時間情報を更新することで）、第３のボリュームデータを生成する。
【００５３】
なお、この様に信号値（輝度値）の最大値を用いる場合、第３のボリュームデータは、解析期間の最終時相に対応する単一のデータのみであっても良い。しなしながら、最終的な単一のデータのみを採用した場合、最も遅く染まるパフュージョンに最も早く染まる腫瘍血管／動脈血管が埋もれてしまい、血管は良く見えづらくなる可能性がある。従って、第３のボリュームデータは、解析期間の各時相に対応する複数の第１のボリュームデータと同じ数だけ作成することが好ましい。
【００５４】
［第２、第３のボリュームデータを用いた投影画像の生成：ステップＳ６］
次に、画像生成ユニット２５は、第２、第３のボリュームデータを用いて投影画像の生成する（ステップＳ６）。
【００５５】
図７は、第２、第３のボリュームデータを用いた投影処理を説明するための図である。また、図８は、当該投影処理において用いられるカラーマップの一例を示した図である。図８のカラーマップは、各画素につき、時間と輝度値（明るさ）とに応じてことなる色調を割り当てるものである。
【００５６】
画像生成ユニット２５は、第３のボリュームデータに対して、図７に示すように奥行き方向を視線方向とする投影処理を各レイについて実行し、且つ第２のボリュームデータと図８に示すカラーマップとを用いて造影剤時間情報に応じた色調を位置毎に選択することで、三次元領域の各位置における造影剤特徴量が輝度値により現され、造影剤時間情報が色調で現される投影画像（四次元パラメトリックイメージ）を生成する。
【００５７】
ここで、投影する範囲において、各視線方向毎に（すなわちレイ毎に）に時間情報は複数存在することになる（言い換えれば、奥行き方向に関する情報が加わることになる）。この様な造影剤に関する時間情報に応じた色相の割り当ては、例えば以下の何れかの実施例に係る手法によって表現することができる。
【００５８】
（分散値表示）
本実施例に係る手法は、視線方向に沿って複数存在する時間情報を分散値によって表現するものである。例えば、腫瘍などの血管が入り組んだ部分では、奥行き方向に時間情報は分散が大きい。一方、大きな血管が無く、毛細血管（パフュージョン）がほとんどである部分の場合、時間情報は分散が小さい。すなわち、分散に応じて色を変えることで、染影が均一か不均一か、血管等を多く含むか含まないかなどを表現することが可能である。
【００５９】
図９は、ある一つのレイ上に存在する第２のボリュームデータのボクセル値を、視点手前側から奥側にかけてプロットすることで得られるグラフである。本実施例の場合、画像生成ユニット２５は、図９に示すように、第２のボリュームデータのボクセルのボクセル値の分散を各レイ毎に計算し、得られた分散値に応じた色調を割り当てる。
【００６０】
なお、奥行きの範囲内に例えば他の臓器が入ってしまう場合は、不必要に分散値が大きくなると考えられる。従って、投影範囲の設定は重要であると言える。投影範囲の最適化は、後述する手法によって実現することが可能である。
【００６１】
（平均値表示）
本実施例に係る手法は、視線方向に沿って複数存在する時間情報を平均値によって表現するものである。すなわち、画像生成ユニット２５は、第２のボリュームデータのボクセルのボクセル値の平均値を各レイ毎に計算し、得られた平均値に応じた輝度値を割り当てる。
【００６２】
（到達時刻表示）
本実施例に係る手法は、視線方向に沿って複数存在する時間情報を到達時刻を基準として表現するものである。例えば、画像生成ユニット２５は、図１０に示すように、到達時刻が最小値（ただし、時間情報＝０であるものは除く）となるボクセルを見つけ、その到達時刻で色調を決め、そのボクセルに対応する位置における輝度値を、第1のボリュームデータによって決定することで、投影処理を実行する。
【００６３】
（ピーク時刻表示）
本実施例に係る手法は、視線方向に沿って複数存在する時間情報をピーク時刻（最大値）を基準として表現するものである。例えば、画像生成ユニット２５は、図１１に示すように、ピーク時刻となるボクセルを見つけ、その到達時刻で色調を決め、そのボクセルに対応する位置における輝度値を、第1のボリュームデータによって決定することで、投影処理を実行する。なお、ピーク時刻が複数ある場合には、手前側に位置するボクセルの情報を優先するものとする。
【００６４】
（その他）
また、上記各実施例の他、例えば、第３のボリュームデータにおいて奥行き方向に輝度値が最大となるボクセル（或いは最小となるボクセル）を特定し、当該ボクセルの位置に対応する第２のボリュームでのボクセルの時間情報を用いて、色調を選択することも可能である。
【００６５】
［投影画像の表示：ステップＳ７］
生成された投影画像は、画像合成ユニット２７において所定の文字情報等と合成され、モニター１４に所定の形態で表示される（ステップＳ７）。
【００６６】
次に、本四次元パラメトリックイメージングの変形例について説明する。
【００６７】
（変形例１）
本四次元パラメトリックイメージングにおいては、ステップＳ６の投影処理において、投影処理の対象とする奥行き方向のデータの厚み（投影範囲）を調整することが可能である。例えば、ユーザの指示によって、所望する断層像が表示される。画像生成ユニット２５は、奥行き方向に表示された断層像を中心とする所定範囲、表示された断層像から画面奥側に所定範囲、表示された断層像から画面手前側に所定範囲、といった具合に、表示された断層像を基準として奥行き方向に関する厚み（投影範囲）を指定する。なお、この投影範囲は、操作者からのマニュアル操作により、任意に調整することが可能である。
【００６８】
この様な構成により、投影範囲を設定可能とすれば、ユーザは非常に簡便に観察範囲を設定でき、描出血管のつながりも、厚みを増やすことで改善される。特に、投影処理において分散表示を採用する場合には、奥行きの範囲内に例えば他の臓器が入ってしまうと、不必要に分散値が大きくなってしまう。従って、本手法を用いて、投影範囲を正確に設定することは効果的である。
【００６９】
（変形例２）
本四次元パラメトリックイメージングにおいては、ステップＳ６の投影処理において、投影処理の対象とするボクセルを造影剤情報に応じて選択する（変更／調整）することが可能である。例えば、第２のボリュームデータに対して閾値処理を施すことにより、造影剤時間情報のうち時間的に早いもの、遅いボクセルを選択的に表示することも可能である。例えば、ユーザからの指示により造影剤時間情報に所定の閾値を設け、これを用いた閾値処理によって必要のないボクセルを除外し、第２のボリュームデータを投影表示することができる。
【００７０】
（変形例３）
本四次元パラメトリックイメージングにおいては、ステップＳ６の投影処理において、投影処理の対象とする奥行き方向のデータの厚み（投影範囲）を、第１、第２、第３のボリュームデータの全範囲を網羅するために、複数設定することも可能である。
【００７１】
例えば、第１、第２、第３のボリュームデータを、図１２に示すように奥行き方向に対して等幅によってサブボリュームデータＶｓ１〜Ｖｓｎに分割し、各サブボリュームデータを投影範囲としてステップＳ６の投影処理を実行する。サブボリュームの分割幅は、操作者からのマニュアル操作により、任意に設定することが可能である。
【００７２】
また、第１、第２、第３のボリュームデータを、図１３に示すように奥行き方向に対してそれぞれが任意の幅を持つサブボリュームデータＶｓ１〜Ｖｓｎに分割し、各サブボリュームデータを投影範囲としてステップＳ６の投影処理を実行する。各サブボリュームの幅は、操作者からのマニュアル操作により、任意に設定することが可能である。
【００７３】
さらに、第１、第２、第３のボリュームデータを、図１４に示すように奥行き方向に対して一部が重複するサブボリュームデータＶｓ１〜Ｖｓｎに分割し、各サブボリュームデータを投影範囲としてステップＳ６の投影処理を実行する。なお、同図の斜線部は、隣り合うサブボリュームの重複領域を示している。重複幅は、操作者からのマニュアル操作により、任意の設定することが可能である。
【００７４】
これらの場合、投影画像は各サブボリュームに対応した数だけ生成されることになる。従って、ステップＳ７においては、各サブボリュームに対応する各投影画像を個別に表示、或いはマルチビュー表示される。
【００７５】
この様な構成により、第１、第２、第３のボリュームデータの全体について映像化することができる。このため、画像診断に有益な領域を取りこぼすことなく、四次元パラメトリックイメージングによって映像化することができる。
【００７６】
（効果）
以上述べた構成によれば、造影イメージングにおいて、所定期間に亘る前記三次元領域に関する超音波データを用いて、所定期間内の各時相における造影剤（バブル）の空間密度を反映する第１のボリュームデータを生成する。また、当該所定期間内の各時相毎の第１のボリュームデータを用いて、三次元領域の各位置における造影剤時間情報を示す第２のボリュームデータと、三次元領域の各位置における造影剤特徴量を示す第３のボリュームデータと、を生成する。そして、生成された第２のボリュームデータ及び第３のボリュームデータを用いて、三次元領域の各位置における造影剤特徴量が輝度値により現され、且つ造影剤時間情報が色調で現される三次元画像を生成し、所定の形態で表示する。従って、奥行き方向に関する時間情報を適切に扱うことができ、ダイナミックな血流変化、微細な血管構築等の情報を高い視認性を持って簡単に観察することができる。
【００７７】
また、時間情報については、到達時刻、ピーク時刻、Wash-In時間、Wash-Out時間等を必要に応じて選択することができる。さらに、また、投影手法については、分散値、到達時刻を基準等を必要に応じて選択することができる。従って、ユーザは状況に応じて所望のパラメータを設定することができ、自由度の高い四次元パラメトリックイメージングを実現することができる。さらに、例えば他の臓器が入らないように投影範囲を調整することで、好適な四次元パラメトリックイメージングを実現することができる。
【００７８】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。具体的な変形例としては、例えば次のようなものがある。
【００７９】
（１）本実施形態に係る各機能は、当該処理を実行するプログラムをワークステーション等のコンピュータにインストールし、これらをメモリ上で展開することによっても実現することができる。このとき、コンピュータに当該手法を実行させることのできるプログラムは、磁気ディスク（フロッピー（登録商標）ディスク、ハードディスクなど）、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤなど）、半導体メモリなどの記録媒体に格納して頒布することも可能である。
【００８０】
（２）上記実施形態においては、ボクセルデータによって構成されるボリュームデータを用いて四次元パラメトリックイメージング処理を実行する例を示した。しかしながら、当該例に拘泥されず、生データによって構成されるボリュームデータを用いて四次元パラメトリックイメージング処理を実行するようにしてもよい。
【００８１】
また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【００８２】
また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【符号の説明】
【００８３】
１０…超音波診断装置、１１…装置本体、１２…超音波プローブ、１３…入力装置、１４…モニター、２１…超音波送信ユニット、２２…超音波受信ユニット、２３…Ｂモード処理ユニット、２４…ドプラ処理ユニット、２５…画像生成ユニット、２６…画像メモリ、２７…画像合成ユニット、２８…制御プロセッサ（ＣＰＵ）、２９…記憶ユニット、３０…インターフェースユニット、３１…ソフトウェア格納ユニット
【図１】
2011224354.tif 000003
【図２】
2011224354.tif 000004
【図３】
2011224354.tif 000005
【図４】
2011224354.tif 000006
【図５】
2011224354.tif 000007
【図６】
2011224354.tif 000008
【図７】
2011224354.tif 000009
【図８】
2011224354.tif 000010
【図９】
2011224354.tif 000011
【図１０】
2011224354.tif 000012
【図１１】
2011224354.tif 000013
【図１２】
2011224354.tif 000014
【図１３】
2011224354.tif 000015
【図１４】
2011224354.tif 000016