"造影剤が投入された被検体内の三次元領域を所定期間に亘って超音波で走査し、前記三次元領域に関する超音波データを前記所定期間に亘って取得するデータ取得ユニットと、前記所定期間内の解析期間に亘る前記三次元領域に関する超音波データを用いて、前記解析期間内の各時相における第1のボリュームデータを生成すると共に、前記解析期間についての造影剤時間情報を示す第2のボリュームデータと、前記解析期間の前記三次元領域の各位置における造影剤特徴量を示す第3のボリュームデータと、を生成するボリュームデータ生成ユニットと、前記第2のボリュームデータと前記第3のボリュームデータとを用いて、投影画像を生成する画像生成ユニットと、を具備する超音波診断装置。前記画像生成ユニットは、前記三次元領域の各位置について、前記造影剤特徴量に応じて異なる輝度値を割り当てると共に前記造影剤時間情報に応じて異なる色調を割り当てることで、前記投影画像を生成する請求項1記載の超音波診断装置。前記画像生成ユニットは、前記第2のボリュームデータを用いて、各投影方向につき、前記造影剤時間情報の分散値、前記造影剤時間情報の最大値、前記造影剤時間情報の最小値、前記造影剤時間情報の平均値のうちのいずれかを計算し、前記計算された値に応じて、異なる色調を割り当てる請求項2記載の超音波診断装置。前記画像生成ユニットは、前記第3のボリュームデータを用いて、各投影方向につき、最大或いは最小の輝度値を有するボクセルを選択し、前記選択されたボクセルに対応する前記第2のボリュームデータのボクセルを特定し、前記特定されたボクセルの前記造影剤時間情報に従って、異なる色調を割り当てる請求項2記載の超音波診断装置。前記ボリュームデータ生成ユニットは、前記解析期間内の各時相における第1のボリュームデータを用いて、前記三次元領域内の各位置の輝度値の時間的変化を計算し、前記三次元領域内の各位置の信号値の時間的変化を用いて、前記三次元領域内の各位置における造影剤の到達時刻、輝度値が最大値になるピーク時刻、WashIn時間、WashOut時間のいずれかを前記造影剤時間情報とする前記第2のボリュームデータを生成する請求項1乃至4のうちいずれか一項記載の超音波診断装置。前記投影画像を生成する際に用いる少なくとも一つの投影範囲を設定する設定ユニットをさらに具備し、前記画像生成ユニットは、前記設定された少なくとも一つの投影範囲を用いて、少なくとも一つの前記投影画像を生成する請求項1乃至5のうちいずれか一項記載の超音波診断装置。前記投影画像が複数枚生成された場合には、当該複数の投影画像をマルチビュー表示する表示ユニットをさらに具備する請求項6記載の超音波診断装置。造影剤が投入された被検体内の三次元領域を所定期間に亘って超音波で走査することで取得された前記三次元領域に関する超音波データを記憶する記憶ユニットと、前記所定期間内の解析期間に亘る前記三次元領域に関する超音波データを用いて、前記解析期間内の各時相における第1のボリュームデータを生成すると共に、前記解析期間についての造影剤時間情報を示す第2のボリュームデータと、前記解析期間の前記三次元領域の各位置における造影剤特徴量を示す第3のボリュームデータと、を生成するボリュームデータ生成ユニットと、前記第2のボリュームデータと前記第3のボリュームデータとを用いて、投影画像を生成する画像生成ユニットと、を具備する超音波画像処理装置。前記画像生成ユニットは、前記三次元領域の各位置について、前記造影剤特徴量に応じて異なる輝度値を割り当てると共に前記造影剤時間情報に応じて異なる色調を割り当てることで、前記投影画像を生成する請求項8記載の超音波画像処理装置。前記画像生成ユニットは、前記第2のボリュームデータを用いて、各投影方向につき、前記造影剤時間情報の分散値、前記造影剤時間情報の最大値、前記造影剤時間情報の最小値、前記造影剤時間情報の平均値のうちのいずれかを計算し、前記計算された値に応じて、異なる色調を割り当てる請求項9記載の超音波画像処理装置。前記画像生成ユニットは、前記第3のボリュームデータを用いて、各投影方向につき、、最大或いは最小の輝度値を有するボクセルを選択し、前記選択されたボクセルに対応する前記第2のボリュームデータのボクセルを特定し、前記特定されたボクセルの前記造影剤時間情報に従って、異なる色調を割り当てる請求項9記載の超音波画像処理装置。前記ボリュームデータ生成ユニットは、前記解析期間内の各時相における第1のボリュームデータを用いて、前記三次元領域内の各位置の輝度値の時間的変化を計算し、前記三次元領域内の各位置の信号値の時間的変化を用いて、前記三次元領域内の各位置における造影剤の到達時刻、輝度値が最大値になるピーク時刻、WashIn時間、WashOut時間のいずれかを前記造影剤時間情報とする前記第2のボリュームデータを生成する請求項8乃至11のうちいずれか一項記載の超音波画像処理装置。前記投影画像を生成する際に用いる少なくとも一つの投影範囲を設定する設定ユニットをさらに具備し、前記画像生成ユニットは、前記設定された少なくとも一つの投影範囲を用いて、少なくとも一つの前記投影画像を生成する請求項8乃至12のうちいずれか一項記載の超音波画像処理装置。前記投影画像が複数枚生成された場合には、当該複数の投影画像をマルチビュー表示する表示ユニットをさらに具備する請求項13記載の超音波画像処理装置。造影剤が投入された被検体内の三次元領域を所定期間に亘って撮像し、前記三次元領域に関する画像データを前記所定期間に亘って取得するデータ取得ユニットと、前記所定期間内の解析期間に亘る前記三次元領域に関する画像データを用いて、前記解析期間内の各時相における第1のボリュームデータを生成すると共に、前記解析期間についての造影剤時間情報を示す第2のボリュームデータと、前記解析期間の前記三次元領域の各位置における造影剤特徴量を示す第3のボリュームデータと、を生成するボリュームデータ生成ユニットと、前記第2のボリュームデータと前記第3のボリュームデータとを用いて、投影画像を生成する画像生成ユニットと、を具備する医用画像診断装置。前記画像生成ユニットは、前記三次元領域の各位置について、前記造影剤特徴量に応じて異なる輝度値を割り当てると共に前記造影剤時間情報に応じて異なる色調を割り当てることで、前記投影画像を生成する請求項15記載の医用画像診断装置。前記画像生成ユニットは、前記第2のボリュームデータを用いて、各投影方向につき、前記造影剤時間情報の分散値、前記造影剤時間情報の最大値、前記造影剤時間情報の最小値、前記造影剤時間情報の平均値のうちのいずれかを計算し、前記計算された値に応じて、異なる色調を割り当てる請求項16記載の医用画像診断装置。前記画像生成ユニットは、前記第3のボリュームデータを用いて、各投影方向につき、、最大或いは最小の輝度値を有するボクセルを選択し、前記選択されたボクセルに対応する前記第2のボリュームデータのボクセルを特定し、前記特定されたボクセルの前記造影剤時間情報に従って、異なる色調を割り当てる請求項16記載の医用画像診断装置。前記ボリュームデータ生成ユニットは、前記解析期間内の各時相における第1のボリュームデータを用いて、前記三次元領域内の各位置の輝度値の時間的変化を計算し、前記三次元領域内の各位置の信号値の時間的変化を用いて、前記三次元領域内の各位置における造影剤の到達時刻、輝度値が最大値になるピーク時刻、WashIn時間、WashOut時間のいずれかを前記造影剤時間情報とする前記第2のボリュームデータを生成する請求項15乃至18のうちいずれか一項記載の医用画像診断装置。前記投影画像を生成する際に用いる少なくとも一つの投影範囲を設定する設定ユニットをさらに具備し、前記画像生成ユニットは、前記設定された少なくとも一つの投影範囲を用いて、少なくとも一つの前記投影画像を生成する請求項15乃至19のうちいずれか一項記載の医用画像診断装置。前記投影画像が複数枚生成された場合には、当該複数の投影画像をマルチビュー表示する表示ユニットをさらに具備する請求項20記載の医用画像診断装置。しかしながら、二次元パラメトリックイメージングでは、二次元画像データ断層データを用いているため、血流の空間的連続性の表現に限界がある。例えば、断層像に直交する血管は、実質的に点として表示されてしまう等、奥行き方向の構造を持つ血管等を認識することは困難である。一方、血流動態血流の時間変化を観察する場合、時相の異なる複数のボリュームデータを観察するとすれば、従来の二次元画像断層像の診断時に比べ、観察範囲が1次元増えるため、観察時間が増加し、さらに頭で時間変化情報の構築記憶も曖昧になるという問題がある。従って、時系列のボリュームデータを用いたパラメトリックイメージング四次元パラメトリックイメージングを行うことが効果的と思われる。しかしながら、到達時間等の時間情報の三次元的映像化法については、これまでに臨床の点おいて効果的な方法が提案されていない。これは、二次元パラメトリックイメージングの手法を単純に三次元的に行うだけでは、奥行き方向に関する時間情報を適切に扱うことができないからである。特に、近年、四次元走査が可能な超音波診断装置においても、造影イメージング機能が導入されてきている。従って、四次元パラメトリックイメージングの手法を確立し、ダイナミックな血流動態の三次元的な観察を可能にすることは、臨床的価値も高く、強く望まれるところである。本発明は、上記事情を鑑みてなされたもので、四次元パラメトリックイメージングにより、ダイナミックな血流変化、微細な血管構築等の情報を高い視認性を持って簡単に観察可能な超音波診断装置、超音波画像処理装置及び医用画像診断装置を提供することを目的としている。本発明は、上記目的を達成するため、次のような手段を講じている。一実施形態に係る超音波診断装置は、造影剤が投入された被検体内の三次元領域を所定期間に亘って超音波で走査し、前記三次元領域に関する超音波データを前記所定期間に亘って取得するデータ取得ユニットと、前記所定期間内の解析期間に亘る前記三次元領域に関する超音波データを用いて、前記解析期間内の各時相における第1のボリュームデータを生成すると共に、前記解析期間についての造影剤時間情報を示す第2のボリュームデータと、前記解析期間の前記三次元領域の各位置における造影剤特徴量を示す第3のボリュームデータと、を生成するボリュームデータ生成ユニットと、前記第2のボリュームデータと前記第3のボリュームデータとを用いて、投影画像を生成する画像生成ユニットと、を具備するものである。以下、本発明の実施形態を図面に従って説明する。なお、以下の説明において、略同一の機能及び構成を有する構成要素については、同一符号を付し、重複説明は必要な場合にのみ行う。図1は、本実施形態に係る超音波診断装置1のブロック構成図を示している。同図に示すように、本超音波診断装置10は、超音波プローブ12、入力装置13、モニター14、超音波送信ユニット21、超音波受信ユニット22、Bモード処理ユニット23、ドプラ処理ユニット24、画像生成ユニット25、画像メモリ26、画像合成ユニット27、制御プロセッサCPU28、記憶ユニット29、インターフェースユニット30、ソフトウェア格納ユニット31を具備している。装置本体11に内蔵される超音波送信ユニット21、受信ユニット22等は、集積回路などのハードウェアで構成されることもあるが、ソフトウェア的にモジュール化されたソフトウェアプログラムである場合もある。以下、個々の構成要素の機能について説明する。超音波プローブ12は、超音波送信ユニット21からの駆動信号に基づき超音波を発生し、被検体からの反射波を電気信号に変換する複数の圧電振動子、当該圧電振動子に設けられる整合層、当該圧電振動子から後方への超音波の伝播を防止するバッキング材等を有している。当該超音波プローブ12から被検体Pに超音波が送信されると、当該送信超音波は、体内組織の音響インピーダンスの不連続面で次々と反射され、エコー信号として超音波プローブ12に受信される。このエコー信号の振幅は、反射することになった反射することになった不連続面における音響インピーダンスの差に依存する。また、送信された超音波パルスが、移動している血流や心臓壁等の表面で反射された場合のエコーは、ドプラ効果により移動体の超音波送信方向の速度成分を依存して、周波数偏移を受ける。なお、本超音波装置が具備する超音波プローブ12は、被検体の三次元領域を超音波走査可能なものである。そのため、超音波プローブ12は、振動子をその配列方向の直交方向に沿って機械的に揺動させ、三次元領域を超音波走査する構成、又は二次元的に配列された二次元振動素子を用いて電気的制御により三次元領域を超音波走査する構成等を有する。前者の構成を採用する場合、被検体の三次元的走査は揺動回路揺動機構によって行われるため、検査者はプローブ本体を被検体に接触させるだけで、自動的に複数の二次元断層像を取得することができる。制御された揺動速度から断面間の正確な距離も検知できる。また、後者の構成を採用する場合には、原理的には、従来の二次元断層像を取得するのと同じ時間で、三次元領域を超音波走査することができる。入力装置13は、装置本体11に接続され、オペレータからの各種指示、条件、関心領域ROIの設定指示、種々の画質条件設定指示等を装置本体11にとりこむためのトラックボール、各種スイッチ、ボタン、マウス、キーボード等を有している。モニター14は、画像合成ユニット27らのビデオ信号に基づいて、生体内の形態学的情報や、血流情報を画像として表示する。超音波送信ユニット21は、図示しないトリガ発生回路、遅延回路およびパルサ回路等を有している。パルサ回路では、所定のレート周波数frHz周期;1/fr秒で、送信超音波を形成するためのレートパルスが繰り返し発生される。また、遅延回路では、チャンネル毎に超音波をビーム状に集束し且つ送信指向性を決定するのに必要な遅延時間が、各レートパルスに与えられる。この遅延情報を変化させることで、プローブ振動子面からの送信方向を任意に調整することが可能となる。トリガ発生回路は、このレートパルスに基づくタイミングで、プローブ12に駆動パルスを印加する。超音波受信ユニット22は、図示していないアンプ回路、A/D変換器、加算器等を有している。アンプ回路では、プローブ12を介して取り込まれたエコー信号をチャンネル毎に増幅する。A/D変換器では、増幅されたエコー信号に対し受信指向性を決定するのに必要な遅延時間を与え、その後加算器において加算処理を行う。この加算により、エコー信号の受信指向性に応じた方向からの反射成分が強調され、受信指向性と送信指向性とにより超音波送受信の総合的なビームが形成される。Bモード処理ユニット23は、超音波受信ユニット22からエコー信号を受け取り、対数増幅、包絡線検波処理などを施し、信号強度が輝度の明るさで表現されるデータを生成する。このとき検波周波数を変化させることで、映像化する周波数帯域を変えることができる。また一つの受信データに対して、2つの検波周波数による検波処理を並列に行うことも可能な構成である。それを利用すると一つの受信信号からバブル像と組織像を生成することも可能となる。なお、バブルを映像化する場合は、セカンドハーモニック2次高調波を主に映像化することが多い。本Bモード処理ユニット23において処理を受けたデータは、画像生成ユニット25に出力され、反射波の強度を輝度にて表したBモード画像として再構成される。ドプラ処理ユニット24は、送受信ユニット22から受け取ったエコー信号から速度情報を周波数解析し、ドプラ効果による血流や組織、造影剤エコー成分を抽出し、平均速度、分散、パワー等の血流情報を多点について求める。得られた血流情報は画像生成ユニット24に送られ、平均速度画像、分散画像、パワー画像、これらの組み合わせ画像としてモニター14にカラー表示することができる。画像生成ユニット25は、超音波スキャンの走査線信号列を、テレビなどに代表される一般的なビデオフォーマットの走査線信号列に変換し、表示画像としての超音波診断画像を生成する。画像生成ユニット24は、画像データを格納する記憶メモリを搭載しており、例えば診断の後に操作者が検査中に記録された画像を呼び出すことが可能となっている。また、画像生成ユニット25は、画像処理装置としての機能を持ち、例えばボリュームデータを構築する場合、三次元領域或いは連続する二次元領域を超音波走査することで得られる走査線信号列を空間的に配置し、座標変換、補間処理等を必要に応じて実行することで、ボリュームデータを構築する。画像生成ユニット25は、得られたボリュームデータを用いてボリュームレンダリングや、空間的MIP等の処理、ボリュームデータ内の任意の断層像を切り出しによるMPR処理、マイクロフローイメージング処理、動き補正処理等を実行し、所定の三次元画像を生成する。なお、本画像生成ユニット25における各種画像処理法等は、ソフトウェア的手法或いはハードウェア的手法のいずれであっても構わない。さらに、画像生成ユニット25は、制御プロセッサ28の制御のもと、後述する四次元パラメトリックイメージング機能に従う所定の処理を実行する。画像メモリ26は、複数フレーム或いは複数ボリューム分の超音波データを一時的に記憶する。画像合成ユニット27は、画像生成ユニット25から受け取った画像を種々のパラメータの文字情報や目盛等と共に合成し、ビデオ信号としてモニター14に出力する。制御プロセッサCPU28は、情報処理装置計算機としての機能を持ち、本超音波診断装置本体の動作を制御する。制御プロセッサ28は、記憶ユニット29から各種画像処理法を実現するためのプログラム、後述する四次元パラメトリックイメージング機能を実現するためのプログラムを読み出してソフトウェア格納部31上に展開し、各種処理に関する演算制御等を実行する。記憶ユニット29は、各種スキャンシーケンスを実行するためのプログラム、後述する四次元パラメトリックイメージング機能を実現するための専用プログラム、画像生成、表示処理を実行するための制御プログラム、診断情報患者ID、医師の所見等、診断プロトコル、送受信条件、ボディマーク生成プログラムその他のデータ群が保管されている。また、必要に応じて、画像メモリ26中の画像の保管などにも使用される。記憶ユニット29のデータは、インターフェースユニット30を経由して外部周辺装置へ転送することも可能となっている。インターフェースユニット30は、入力装置13、ネットワーク、新たな外部記憶装置図示せずに関するインターフェースである。当該装置によって得られた超音波画像等のデータや解析結果等は、インターフェースユニット30よって、ネットワークを介して他の装置に転送可能である。四次元パラメトリックイメージング機能次に、本超音波診断装置1が有する四次元パラメトリックイメージング機能について説明する。この機能は、造影イメージングにおいて、所定期間に亘る前記三次元領域に関する超音波データを用いて、所定期間内の各時相における造影剤バブルの空間密度を信号値輝度値として反映する第1のボリュームデータを生成する。また、当該所定期間内の各時相毎の第1のボリュームデータを用いて、三次元領域の各位置における造影剤時間情報を示す第2のボリュームデータと、三次元領域の各位置における造影剤特徴量を示す第3のボリュームデータと、を生成する。そして、生成された第2のボリュームデータ及び第3のボリュームデータを用いて、三次元領域の各位置における造影剤特徴量が輝度値により現され、且つ造影剤時間情報が色調で現される三次元画像を生成し、所定の形態で表示するものである。なお、以下の説明においては、説明を具体的にするために、超音波診断装置を用いて四次元パラメトリックイメージング機能を実現する場合を例とする。しかしながら、当該例に拘泥されず、例えば、例えば医用ワークステーションによる超音波画像処理装置或いは、ビューワを用いて、四次元走査によって事前に取得された超音波データを用いても、行うことができる。係る場合には、例えば、図1の点線内の構成を採用することができる。さらに、造影剤を用いる他の医用画像診断装置X線コンピュータ断層撮像装置、磁気共鳴イメージング装置、X線診断装置等においても、例えば図1の点線内の構成と実質的に同じ構成を実行することで、四次元パラメトリックイメージング機能を実現することが可能である。図2は、本四次元パラメトリックイメージング機能に従う処理四次元パラメトリックイメージング処理の流れを示したフローチャートである。以下、当該フローチャートに示す各ステップにおいて実行される処理の内容について説明する。[患者情報の入力、送受信条件、スキャンシーケンス等の選択:ステップS1]操作ユニット33を介して患者情報の入力、送受信条件被走査領域の大きさを決めるための画角、焦点位置、送信電圧等、被検体の三次元領域を所定期間に亘って超音波走査するためのスキャンシーケンス等の選択が実行されるステップS1。入力、選択された各種情報条件等は、自動的に記憶装置29に記憶される。[所定期間に亘るボリュームデータの収集:ステップS2]次に、送受信制御ユニット31は、被検者の所定部位診断対象とする血管等を含む三次元領域を被走査領域として、リアルタイム三次元超音波走査四次元走査を所定期間に亘って実行し、当該所定期間内の各時相に対応する三次元超音波データを取得するステップS2。[造影剤の空間密度を表す第1のボリュームデータの生成:ステップS3]次に、画像生成ユニット25は、取得された三次元超音波データに対して、造影剤からの信号成分に対応するハーモニック成分の抽出処理、座標変換、補間処理等を施し、例えば図3に示すような各時相kk=1,2,、nに対応する第1のボリュームデータを生成するステップS3。この第1のボリュームデータは、各時相における造影剤の空間密度濃度を信号値輝度値で表すものとなっている。従って、時系列の当該第1のボリュームデータを解析することで、被走査領域である三次元領域内の各位置における造影剤密度の時間的変化を把握することが可能である。[造影剤に関する時間情報を示す第2のボリュームデータの生成:ステップS4]次に、画像生成ユニット25は、所定期間に亘る第1のボリュームデータのうち、解析期間所定期間のうち、四次元パラメトリックイメージングによる解析対象とする期間に対応する時相毎の第1のボリュームデータを用いて、造影剤に関する時間情報を示す第2のボリュームデータを生成するステップS4。ここで、造影剤に関する時間情報とは、各位置における造影剤の到達時刻、信号値又は輝度値のピーク時刻、WashIn時間、WashOut時間等である。なお、解析期間の開始時刻及び終了時刻は、入力装置13からの操作により任意に設定することが可能である。また、例えば、造影剤タイマーのON時刻、バブルを破壊するバーストの終了時刻、ある対象にバブルが入ってきた到達した時刻等を検出し、装置によって自動的に解析期間を開始等するようにしてもよい。また、第2のボリュームデータは、解析期間に対応する複数の第1のボリュームデータと第1の各ボリュームの時相と一致して同じ数だけ作成することを前提とする。しかしながら、これに限定されることなく、例えば第2のボリュームデータを単一のボリュームデータとすることも可能である。以下、本ステップにおいて生成される第2のボリュームを造影剤時間情報の種別毎に説明する。造影剤時間情報を造影剤の到達時刻とする場合まず、画像生成ユニット25は、解析期間の各時相に対応する複数の第1のボリュームデータを用いて、各ボクセルについての解析期間におけるTICTimeIntensityCurveを生成する。図4は、あるボクセルについての解析期間に亘るTICの一例を示した図である。画像生成ユニット25は、生成された各ボクセルについての解析期間におけるTICを用いて各ボクセルの輝度の時間変化を観察し、図4に示すような予め設定された輝度の閾値を始めて越えた時刻Tarrivedを、ボリュームレートとボリューム数とを用いてボクセル毎に計算する。この計算された時刻をボクセルの値とすることで、第2のボリュームデータを生成することができる。また、解析期間の最初の時相に対応する第1のボリュームデータにおいて、輝度値が最初から閾値を越えているボクセルについては、ボクセル値をゼロまたは値なしとしてもよい。なぜなら、この様なボクセルの輝度値は、バブルからの信号に起因するものではなく、組織からの信号に起因するものが多いからである。この様に、造影剤時間情報が造影剤の到達時刻である場合には、解析期間に対応する複数の第1のボリュームデータと同じ数だけ作成してもよい。しかしながら、例えば解析期間内の最後の時相に対応する第2のボリュームデータが一つあれば、各ボクセルについての解析期間における造影剤到達時刻を把握することができる。造影剤時間情報を輝度値がピークになるピーク時刻とする場合まず、画像生成ユニット25は、解析期間の各時相に対応する複数の第1のボリュームデータを用いて、各ボクセルについての解析期間におけるTICTimeIntensityCurveを生成する。次に、画像生成ユニット25は、生成された各ボクセルについての解析期間におけるTICを用いて各ボクセルの輝度の時間変化を観察し、図4に示すように輝度値が最大ピークとなる時刻Tpeakをボクセル毎に計算する。この計算された時刻をボクセルの値とすることで、第2のボリュームデータを生成することができる。この様に、造影剤に関する時間情報が、輝度値がピークになるピーク時刻である場合においても、解析期間に対応する複数の第1のボリュームデータと同じ数だけ作成してもよい。しかしながら、到達時刻は一度求まると更新されることはないのに対し、ピーク時間の場合は、例えば図5に示すように時系列で隣り合うボリューム間で輝度値の大小比較演算fk-1,kただし、k=1,2,、nを実行し、輝度が過去の最大値を超えるたびに、時間情報も更新されることになる。厳密に言えば、輝度がピークを迎え、WashOutする時間まで観察するならば、各ボクセルでピーク時間は一つしかないことになる。従って、例えば解析期間内の最後の時相に対応する第2のボリュームデータが一つあれば、各ボクセルについての解析期間における造影剤の到達時刻を把握することができる。造影剤時間情報を造影剤のWashIn時間とする場合WashInの時間を、図4に示すように、例えばピーク時刻-到達時刻として定義する。係る場合には、解析期間についてのTICを用いて各ボクセルの輝度の時間変化を観察し、輝度値がある閾値を超えた時刻到達時刻と最大ピークとなる時刻とをボクセル毎に特定し、各ボクセルについてのWashIn時間を計算することができる。この計算された時間期間をボクセルの値とすることで、造影剤に関する時間情報を造影剤のWashIn時間とする場合の第2のボリュームデータを、例えば単一のボリュームとして生成することができる。造影剤に関する時間情報を造影剤のWashOut時間とする場合WashOutの時間を、図4に示すように、例えばピークを一度迎えた時刻ピーク時刻からTICが所定の閾値以下になる時刻までの期間として定義する。なお、TICは、フィッティングFitting又は平滑化Smoothingされていることが好ましい。係る場合には、解析期間におけるTICを用いて各ボクセルの輝度の時間変化を観察し、輝度値が最大ピークとなる時刻から所定の閾値以下になる時刻までの時間をボクセル毎に計算する。この計算された時間期間をボクセルの値とすることで、造影剤に関する時間情報を造影剤のWashOut時間とする場合の第2のボリュームデータを、例えば単一のボリュームとして生成することができる。なお、以上述べた各種造影剤に関する時間情報に関する第2のボリュームデータの生成においては、同じボクセルが同位置であるという前提の処理である。しかしながら、実際にはいくらプローブをしっかり保持しようとしても、手が動いてしまったり、患者の呼吸で動いてしまったり、完全に固定しておくことは困難である。そこで、例えば特開2007-330764に記載の動き補正技術を用いて装置が連続するボリューム間の位置ずれを補正することは、実際の検査でこの提案手法を使う場合に、非常に有用である。[各位置の造影剤特徴量を示す第3のボリュームデータの生成:ステップS5]次に、画像生成ユニット25は、解析期間に対応する時相毎の第1のボリュームデータを用いて、解析期間における造影剤特徴量を被走査領域の各位置毎に示す第3のボリュームデータを生成するステップS5。ここで、解析期間における造影剤特徴量とは、解析期間での各位置についての信号値輝度値の最大値、最小値、平均値等である。例えば、血管形態観察のためにバブルをトレースするために、解析期間における造影剤特徴量として、信号値輝度値の最大値を用いる場合を想定する。係る場合には、画像生成ユニット25は、解析期間に対応する時相毎の第1のボリュームデータを用いて、空間的位置が対応する各ボクセルについて、時間方向に最大値を保持temporalmaxholdを実行することで例えば図6に示すように時系列で隣り合うボリューム間で輝度値の最大値大小比較演算Maxk-1,kただし、k=1,2,、nを逐次実行し、輝度が過去の最大値を超えるたびに、時間情報を更新することで、第3のボリュームデータを生成する。なお、この様に信号値輝度値の最大値を用いる場合、第3のボリュームデータは、解析期間の最終時相に対応する単一のデータのみであっても良い。しなしながら、最終的な単一のデータのみを採用した場合、最も遅く染まるパフュージョンに最も早く染まる腫瘍血管/動脈血管が埋もれてしまい、血管は良く見えづらくなる可能性がある。従って、第3のボリュームデータは、解析期間の各時相に対応する複数の第1のボリュームデータと同じ数だけ作成することが好ましい。[第2、第3のボリュームデータを用いた投影画像の生成:ステップS6]次に、画像生成ユニット25は、第2、第3のボリュームデータを用いて投影画像の生成するステップS6。図7は、第2、第3のボリュームデータを用いた投影処理を説明するための図である。また、図8は、当該投影処理において用いられるカラーマップの一例を示した図である。図8のカラーマップは、各画素につき、時間と輝度値明るさとに応じてことなる色調を割り当てるものである。画像生成ユニット25は、第3のボリュームデータに対して、図7に示すように奥行き方向を視線方向とする投影処理を各レイについて実行し、且つ第2のボリュームデータと図8に示すカラーマップとを用いて造影剤時間情報に応じた色調を位置毎に選択することで、三次元領域の各位置における造影剤特徴量が輝度値により現され、造影剤時間情報が色調で現される投影画像四次元パラメトリックイメージを生成する。ここで、投影する範囲において、各視線方向毎にすなわちレイ毎にに時間情報は複数存在することになる言い換えれば、奥行き方向に関する情報が加わることになる。この様な造影剤に関する時間情報に応じた色相の割り当ては、例えば以下の何れかの実施例に係る手法によって表現することができる。分散値表示本実施例に係る手法は、視線方向に沿って複数存在する時間情報を分散値によって表現するものである。例えば、腫瘍などの血管が入り組んだ部分では、奥行き方向に時間情報は分散が大きい。一方、大きな血管が無く、毛細血管パフュージョンがほとんどである部分の場合、時間情報は分散が小さい。すなわち、分散に応じて色を変えることで、染影が均一か不均一か、血管等を多く含むか含まないかなどを表現することが可能である。図9は、ある一つのレイ上に存在する第2のボリュームデータのボクセル値を、視点手前側から奥側にかけてプロットすることで得られるグラフである。本実施例の場合、画像生成ユニット25は、図9に示すように、第2のボリュームデータのボクセルのボクセル値の分散を各レイ毎に計算し、得られた分散値に応じた色調を割り当てる。なお、奥行きの範囲内に例えば他の臓器が入ってしまう場合は、不必要に分散値が大きくなると考えられる。従って、投影範囲の設定は重要であると言える。投影範囲の最適化は、後述する手法によって実現することが可能である。平均値表示本実施例に係る手法は、視線方向に沿って複数存在する時間情報を平均値によって表現するものである。すなわち、画像生成ユニット25は、第2のボリュームデータのボクセルのボクセル値の平均値を各レイ毎に計算し、得られた平均値に応じた輝度値を割り当てる。到達時刻表示本実施例に係る手法は、視線方向に沿って複数存在する時間情報を到達時刻を基準として表現するものである。例えば、画像生成ユニット25は、図10に示すように、到達時刻が最小値ただし、時間情報=0であるものは除くとなるボクセルを見つけ、その到達時刻で色調を決め、そのボクセルに対応する位置における輝度値を、第1のボリュームデータによって決定することで、投影処理を実行する。ピーク時刻表示本実施例に係る手法は、視線方向に沿って複数存在する時間情報をピーク時刻最大値を基準として表現するものである。例えば、画像生成ユニット25は、図11に示すように、ピーク時刻となるボクセルを見つけ、その到達時刻で色調を決め、そのボクセルに対応する位置における輝度値を、第1のボリュームデータによって決定することで、投影処理を実行する。なお、ピーク時刻が複数ある場合には、手前側に位置するボクセルの情報を優先するものとする。その他また、上記各実施例の他、例えば、第3のボリュームデータにおいて奥行き方向に輝度値が最大となるボクセル或いは最小となるボクセルを特定し、当該ボクセルの位置に対応する第2のボリュームでのボクセルの時間情報を用いて、色調を選択することも可能である。[投影画像の表示:ステップS7]生成された投影画像は、画像合成ユニット27において所定の文字情報等と合成され、モニター14に所定の形態で表示されるステップS7。次に、本四次元パラメトリックイメージングの変形例について説明する。変形例1本四次元パラメトリックイメージングにおいては、ステップS6の投影処理において、投影処理の対象とする奥行き方向のデータの厚み投影範囲を調整することが可能である。例えば、ユーザの指示によって、所望する断層像が表示される。画像生成ユニット25は、奥行き方向に表示された断層像を中心とする所定範囲、表示された断層像から画面奥側に所定範囲、表示された断層像から画面手前側に所定範囲、といった具合に、表示された断層像を基準として奥行き方向に関する厚み投影範囲を指定する。なお、この投影範囲は、操作者からのマニュアル操作により、任意に調整することが可能である。この様な構成により、投影範囲を設定可能とすれば、ユーザは非常に簡便に観察範囲を設定でき、描出血管のつながりも、厚みを増やすことで改善される。特に、投影処理において分散表示を採用する場合には、奥行きの範囲内に例えば他の臓器が入ってしまうと、不必要に分散値が大きくなってしまう。従って、本手法を用いて、投影範囲を正確に設定することは効果的である。変形例2本四次元パラメトリックイメージングにおいては、ステップS6の投影処理において、投影処理の対象とするボクセルを造影剤情報に応じて選択する変更/調整することが可能である。例えば、第2のボリュームデータに対して閾値処理を施すことにより、造影剤時間情報のうち時間的に早いもの、遅いボクセルを選択的に表示することも可能である。例えば、ユーザからの指示により造影剤時間情報に所定の閾値を設け、これを用いた閾値処理によって必要のないボクセルを除外し、第2のボリュームデータを投影表示することができる。変形例3本四次元パラメトリックイメージングにおいては、ステップS6の投影処理において、投影処理の対象とする奥行き方向のデータの厚み投影範囲を、第1、第2、第3のボリュームデータの全範囲を網羅するために、複数設定することも可能である。例えば、第1、第2、第3のボリュームデータを、図12に示すように奥行き方向に対して等幅によってサブボリュームデータVs1~Vsnに分割し、各サブボリュームデータを投影範囲としてステップS6の投影処理を実行する。サブボリュームの分割幅は、操作者からのマニュアル操作により、任意に設定することが可能である。また、第1、第2、第3のボリュームデータを、図13に示すように奥行き方向に対してそれぞれが任意の幅を持つサブボリュームデータVs1~Vsnに分割し、各サブボリュームデータを投影範囲としてステップS6の投影処理を実行する。各サブボリュームの幅は、操作者からのマニュアル操作により、任意に設定することが可能である。さらに、第1、第2、第3のボリュームデータを、図14に示すように奥行き方向に対して一部が重複するサブボリュームデータVs1~Vsnに分割し、各サブボリュームデータを投影範囲としてステップS6の投影処理を実行する。なお、同図の斜線部は、隣り合うサブボリュームの重複領域を示している。重複幅は、操作者からのマニュアル操作により、任意の設定することが可能である。これらの場合、投影画像は各サブボリュームに対応した数だけ生成されることになる。従って、ステップS7においては、各サブボリュームに対応する各投影画像を個別に表示、或いはマルチビュー表示される。この様な構成により、第1、第2、第3のボリュームデータの全体について映像化することができる。このため、画像診断に有益な領域を取りこぼすことなく、四次元パラメトリックイメージングによって映像化することができる。効果以上述べた構成によれば、造影イメージングにおいて、所定期間に亘る前記三次元領域に関する超音波データを用いて、所定期間内の各時相における造影剤バブルの空間密度を反映する第1のボリュームデータを生成する。また、当該所定期間内の各時相毎の第1のボリュームデータを用いて、三次元領域の各位置における造影剤時間情報を示す第2のボリュームデータと、三次元領域の各位置における造影剤特徴量を示す第3のボリュームデータと、を生成する。そして、生成された第2のボリュームデータ及び第3のボリュームデータを用いて、三次元領域の各位置における造影剤特徴量が輝度値により現され、且つ造影剤時間情報が色調で現される三次元画像を生成し、所定の形態で表示する。従って、奥行き方向に関する時間情報を適切に扱うことができ、ダイナミックな血流変化、微細な血管構築等の情報を高い視認性を持って簡単に観察することができる。また、時間情報については、到達時刻、ピーク時刻、WashIn時間、WashOut時間等を必要に応じて選択することができる。さらに、また、投影手法については、分散値、到達時刻を基準等を必要に応じて選択することができる。従って、ユーザは状況に応じて所望のパラメータを設定することができ、自由度の高い四次元パラメトリックイメージングを実現することができる。さらに、例えば他の臓器が入らないように投影範囲を調整することで、好適な四次元パラメトリックイメージングを実現することができる。なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。具体的な変形例としては、例えば次のようなものがある。1本実施形態に係る各機能は、当該処理を実行するプログラムをワークステーション等のコンピュータにインストールし、これらをメモリ上で展開することによっても実現することができる。このとき、コンピュータに当該手法を実行させることのできるプログラムは、磁気ディスクフロッピー登録商標ディスク、ハードディスクなど、光ディスクCD-ROM、DVDなど、半導体メモリなどの記録媒体に格納して頒布することも可能である。2上記実施形態においては、ボクセルデータによって構成されるボリュームデータを用いて四次元パラメトリックイメージング処理を実行する例を示した。しかしながら、当該例に拘泥されず、生データによって構成されるボリュームデータを用いて四次元パラメトリックイメージング処理を実行するようにしてもよい。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。"