被検体に超音波を送受信する探触子と、前記探触子が受信した超音波信号を整相加算する受信ビームフォーマ部と、前記被検体に造影剤が投与された開始時間を含む情報の指定を受け付ける処理画像選択部と、前記開始時間の情報に基づいて、整相加算された信号に基づく画像データから、前記造影剤で造影された部分を抽出する第1画素選択部と、前記第1画素選択部で抽出された画像データの画素情報に基づいて、更に画像データを抽出する第2画素選択部と、第2画素選択部で抽出された画像データに基づいて血流動態画像を生成する画像構成部と、前記血流動態画像を表示する表示部と、を有する超音波診断装置。前記第2画素選択部は、更に、前記開始時間に基づいて算出された情報に基づいて、画像を抽出することを特徴とする超音波診断装置。前記画像構成部にて構成される血流動態画像に対して、第1のROIで指定する領域のマッピングパラメータの数値を前記表示部に表示し、さらに第2のROI指定する領域のマッピングパラメータの数値又は第1のROIの数値に対する相対的な値を前記表示部に表示する機能を有することを特徴とする請求項1に記載の超音波診断装置。前記処理画像選択部では、指定を受け付けた開始時間に対して、予め設定された時間tstarttpreだけさかのぼった画像データから抽出が開始されることを特徴とする請求項1に記載の超音波診断装置。前記血流動態画像を構成する上で必要となる数値を画素毎に算出するTIC評価部を有し、前記TIC評価部は、閾値輝度Iαから造影前の範囲における最大値Ipremaxを差し引いた輝度変化量ΔIを造影剤による輝度変化量として算出する評価部を有し、前記第1画素選択部は、算出された輝度変化量ΔIに基づいての画素を選択し、前記第2画素選択部は、前記輝度変化量ΔIから算出された時間情報に基づいて着目する画素の抽出をすることを特徴とする請求項2に記載の超音波診断装置。前記画像構成部は、前記時間情報の時間の変化に伴い色が変化する時間カラーバーと、前記カラーバーで特定された色で、画素毎に抽出された前記画像データを色づけした血流動態画像を含むことを特徴とする請求項5に記載の超音波診断装置。前記第2画素選択部は、前記時間カラーバーの色の変化の時間幅を変更可能なことを特徴とする請求項6に記載の超音波診断装置。前記第2画素選択部は、輝度変化量ΔIと画素の関係を示すヒストグラムに基づいて、表示する画素を選択し、ヒストグラム上に表示される描画範囲を受け付け、受け付けた範囲より外の画素のパラメータが0値化され血流動態画像に反映させることを特徴とする請求項5に記載の超音波診断装置。前記第2画素選択部で利用する時間カラーバーは表示画面に表示され、始点と終点を術者が指定することにより、その範囲外の画素のパラメータが0値化され、さらに始点と終点を時間カラーバーの上限と下限とする色の配置が新たに設定され、血流動態画像に即時反映されることを特徴とする請求項6に記載の超音波診断装置。前記第2画素選択部は、輝度変化量ΔIと画素の関係を示すヒストグラムに対応した、輝度の大小に伴い色が変化する輝度のカラーバー基づいて、表示する画素を選択し、前記輝度カラーバーは、ヒストグラムを等分化する数と等分化した各範囲での画素の総数に応じて、輝度カラーバーを区分けする数と各区間の幅を調整したものであり、更に輝度カラーバーにはΔIの大小の他、血管の太さや種類の情報を表示することを特徴とする請求項5に記載の超音波診断装置。表示画面上に設ける設定変更の機能により、前記処理画像選択部で用いる時間幅tstarttpre、前記TIC計測部で算出するパラメータ、前記第1画素選択部で設定する条件式を術者が自由に変更できることを特徴とする請求項1に記載の超音波診断装置。前記表示部では構成した血流動態画像と過去に構成した血流動態画像を比較する場合、比較する2つの血流動態画像上に基準とするROIを設けて各ROIの色と時間の関係を定量的に比較し、基準とするROIのマッピングパラメータと配色を一致させる補間機能を有することを特徴とする請求項1に記載の超音波診断装置。前記画像データは、前記受信ビームフォーマから出力されるRAWデータ含んでいることを特徴とする請求項1に記載の超音波診断装置。超音波探触子から、被検体に超音波を送受信し、前記超音波探触子が受信した超音波信号を、受信ビームフォーマ部で整相加算し、処理画像選択部で、前記被検体に造影剤が投与された開始時間を含む情報の指定を受け付け、第1画素選択部で、前記開始時間の情報に基づいて、整相加算された信号に基づく画像データから、前記造影剤で造影された部分を抽出し、第2画素選択部で、前記第1画素選択部で抽出された画像データの画素情報に基づいて、更に画像データを抽出し、画像構成部は、第2画素選択部で抽出された画像データに基づいて血流動態画像を生成し、表示部で前記血流動態画像を表示する超音波画像表示方法。造影剤の流入時間など血流動態の相対的な違いを色の違いで表現する血流動態画像を表示する機能に関して、着目する血管を選択的に抽出しノイズが少ない血流動態画像を構成する機能と、治療等による血流動態の変化を調べるために、構成した血流画像と過去の血流動態画像との色の違いを比較する際に、基準となる領域の色を統一化する補間処理機能を持つ超音波診断装置及び超音波画像表示方法を提供する。本発明は、処理画像選択部11で範囲指定した画像データを用いて、TIC評価部13にて造影剤の流入開始時間を示すパラメータ及び造影剤による輝度変化量ΔIを算出する。第1画素選択部15で一定の条件を満たす画素を選択的に抽出すると同時に非造影部を除去する。画像構成部16にてtαの値に応じて配色した血流動態画像を構成し、第2画素選択部17にて術者が着目する血管を輝度変化量に基づくヒストグラムを利用して選択的に抽出する。本発明により、着目血管を選択的に抽出した視認性の高い血流動態画像を提示できる。以下、本発明の実施例について図を用いて説明する。実施例になる超音波診断装置のブロック図を図1に示す。実施例になる超音波診断装置は、撮像対象1に対して超音波信号を送受信する探触子2と、探触子2を構成する圧電素子に所望の送受信ビームを形成するための時間遅延を与える送信ビームフォーマ部4及び受信ビームフォーマ部7と、送受信信号をアナログデジタル変換するAD変換器6及びDA変換器3と、超音波信号が生体内部を伝播する過程で生じる振幅減衰を補正するTGCTimeGainController5と、受信したRF信号を検波し画像信号に変換する包落線検波部8と、画像信号から二次元の画像データを構成するスキャンコンバータ9とを持ち、ここまでが一般的に知られる超音波画像を構成する画像データ構成部である。本実施例では更に、スキャンコンバータ9からの画像データを保存する画像メモリ12と、被検体に造影剤が投与された開始時間を含む情報の指定を受け付けるなど、保存した画像データから血流動態画像の構成に利用するデータ範囲を選択する処理画像選択部11と、選択した画像データから血流動態の評価に必要なパラメータを求めるTIC評価部13と、TIC評価部13で算出した数値を保存するパラメータメモリ14と、予め装置内に設定される条件式に基づき画素の抽出及び除去を行う第1画素選択部15と、パラメータメモリ14に保存された数値を用いて血流動態画像を構成する画像構成部16と、画像構成部16で構成した画像に対して、輝度情報と時間情報とを利用して再度画素の抽出を行い、術者が着目する情報を抽出する第2画素選択部17とを有する。前記画像データ構成部で構成される画像データは、臨床現場で汎用的に使用されている白黒画像Bモードや造影剤画像送受信のシーケンスやフィルタ処理等により造影剤からの信号を強調した画像を指し、その構成手法は一般的に既知であるため、ここでは簡単に説明する。探触子2の超音波照射面は、複数の圧電素子が一列に配列した1次元アレイの構成になっており、各素子が超音波の送受信を担う。送信ビームフォーマ部4からの電圧パルスがDA変換器3を介して各圧電素子に入力され、素子の圧電振動によって撮像対象1に向けて超音波が照射される。この時、各圧電素子には所定の時間遅延が電子的に与えられており、各圧電素子から送信された超音波は撮像対象1の内部の所定の位置で焦点を結ぶ。撮像対象1からの反射エコーは各圧電素子で受信され、TGC5で伝播距離に応じた振幅補正がなされる。続いて、受信信号はAD変換器6を介して受信ビームフォーマ部7に送られ、焦点位置から各圧電素子までの距離に応じた遅延時間を掛けて加算結果が出力される整相加算。造影剤からの信号を強調して画像化する手法としては、例えば互いの位相が反転した2つの信号を送信し、その受信信号を加算する方法が良く知られている。受波信号の加算により組織成分を主とする基本周波数成分が抑圧される一方で、造影剤からの信号を主とする高調波成分は増強される。その結果、造影剤を強調した画像が得られる。超音波送受信を、圧電素子列に沿った全ての走査線で行なうことで、撮像対象1の2次元的な反射エコー分布が得られる。受信ビームフォーマ部7からは実部と虚部に分けられたRF信号が出力され、包絡線検波部8に送られる。包絡線検波部8に送られた信号は、ビデオ信号に変換された後、スキャンコンバータ9で走査線間の画素補間が加えられ、二次元の画像データに再構成された後、表示部10に表示される。次に、スキャンコンバータ9からの画像データを用いて血流動態画像を構成するまでの超音波画像表示方法を図2に示すフローチャートに従って説明する。まず術者は、撮像対象1を撮像面内に捉えた状態で探触子2を固定し、造影剤を投与する。造影剤の投与と同時又は投与前後に術者は装置上の画像データの取り込み開始ボタンを押し、検査対象に造影剤が流入する段階から充満する段階までの時系列画像データを画像メモリ12に保存する工程0。次に、取り込んだ画像データから血流動態の評価に用いる画像データを選択する工程1。利用する画像データのセットを画像メモリ12から選択すると図3に示すように選択した画像データ32、血流動態画像の処理の開始を示す開始ボタン31、画像データの時間範囲を示す時間レンジバー33が表示される。時間レンジバー33の単位は初期値を0秒とする秒、又は初期値を1とするフレーム番号が表示される。時間レンジバー33には現画像の時間位置を示す矢印が設けられ、操作者が指定する時間位置の画像データが逐次表示される。時間レンジバー33には処理する画像データの範囲を指定するための開始矢印と終了矢印も設けられており、術者は画像データを見ながら血流動態画像の処理に利用する画像データの開始及び終了の範囲を指定する工程1。この時、血流動態の観察を主たる目的とする上では、一般的には開始矢印は撮像対象1に造影剤が流入する直前又は直後の時間に設定される。終了矢印は、着目する血管への造影剤の流入が完了した時点、又は血管造影を経て組織の造影が確認できた時点が適当である。前者の場合、例えば病変部の血管が動脈支配なのか門脈支配なのかの判断に有効な情報が得られる。後者の場合は、血管だけでなく血管から組織領域への造影経過も色の違いで観察できるため、病変の鑑別や病変の主たる栄養血管の同定に役立つ情報が得られる。特に血管を標的とする治療法を実施する場面では、標的の同定や治療効果の判定に非常に有効であるが、その分処理時間が増大する。超音波診断装置によっては検査対象に充満した撮像面内の造影剤を高音圧超音波によって一時的に消失させる機能を有し、この機能を利用して高音圧超音波の照射後に撮像面内に再度流入する造影剤を観察する場合がある。この場合には、例えば各画像データの輝度の積分値や平均値から高音圧照射が実施された時間に対応する画像データが容易に特定でき、この画像データの位置に開始矢印が自動的に設置される。また画像データ上のROIを利用して空間的な範囲指定が実施し、処理コストの低減を図ることもできる。画像データの範囲を指定した後、開始ボタン31を選択すると工程2で指定された開始時間から更に1秒前の時間が新たな開始時間として自動的に再設定され、当該範囲の画像データがTIC評価部13に送られる。術者が指定する開始時間の直前の1秒間の画像データは、造影剤以外の要因による輝度変動量を算出するために利用される。1秒という数値は1心拍に対応する時間として予め設定され、この間の輝度変動は主に心拍に起因して生じるスペックルパタンの変動と見なせる。算出された輝度変動量はTICに常に内在し、TICから造影剤による輝度変動量を正確に算出するために必要な値となる。なお1秒という設定は術者が自由に変更することが可能であり、例えば患者の心拍等に合わせて適宜調整される。新たに設定される開始時間が取得した画像データの範囲を超える場合には、0秒が開始時間として設定される。TIC評価部13では、造影前の時間範囲tpretstartにおける最大値Ipremax、処理範囲全体tstarttendにおける最大値Imaxと最小値Imin、術者が予め設定する閾値輝度Iα、及び閾値輝度に達した時点の時間tαが画素毎にTICから算出される工程2。TIC計測部13で算出される各パラメータとTICとの関係を図4Aに示す。tpreは前述したように、工程1にて術者が設定した開始時間の1秒前の値をとる。ここでは、血流動態画像を構成する評価値マッピングパラメータを造影剤の流入開始時間を示すtαとして以下説明を続けるが、実際にはTICから計測される最大値や平均値といったTICから求まる様々な統計値がマッピングパラメータの候補となり術者が自由に選択できる。また、Iαは最大値Imaxに術者が予め設定する係数αを掛けた値とする。係数αの値は術者の任意で決められるが、造影剤の流入開始時間をマッピングパラメータとする場合には、経験的にαは0508程度が適当である。適切なtαを算出するには、図4Bに示すように閾値輝度がTICの斜面、つまり造影剤流入の最盛期を捉える必要がある。図中には不適切なαの設定の一例としてα09、α02の場合を示す。双方共に造影剤の流入とは関係しない位置でTICが閾値輝度に達し、不適切なtαを算出する結果となることが判る。ここではα08として説明を続ける。TIC評価部13においてパラメータを算出する際に、TICを35フレーム程度で平均化処理することが望ましい。この処理により画像データに含まれる電気ノイズの影響をTICから排除し、より正確なパラメータの算出が実現する。このとき、画像データのサンプリング数単位時間当たりの画像データ数は維持するのがよい。次に、TIC計測部13で算出したパラメータに基づき、造影剤による輝度変化量ΔIが画素毎に算出される工程3。TICの変動要因は造影剤の流入と組織の移動に分けられる。探触子を固定し、かつ息止め等で呼吸を制御した場合、組織の移動に伴う輝度変動は主に心拍に伴うスペックルの変動と考えられる。したがってIαから造影前段階での輝度変動の最大値Ipremaxを差し引いた値ΔIを造影剤による輝度変化量として算出する。更には、ΔIの値の大きさは血流量を反映しており、動脈、門脈又は腫瘍血管等の血管の種類を判別する上で有効な指標となる。TIC評価部13で算出した各数値とΔIはパラメータメモリ14に保存される。次に、第1画素選択部15で造影部の画素が選択的に抽出される工程4。選択される画素は以下の2つの条件を満たす。条件1:Iα>Ipremax条件2:ΔI>IpremaxImin第1画素選択部15では、造影部の疑いがある画素は残し、確実に非造影部と判断できる画素、又は確実に血流動態画像のノイズとなる画素を除去する。条件1は、閾値輝度が造影前段階の最大輝度を超えることを意味する。この条件を満たさない画素は、造影段階に入っても造影前段階と同じ輝度変動しかないため、造影剤の流入が認められない非造影部と判断する。逆にこの画素を造影部と判断した場合には、図4Bに示すように誤ったマッピングパラメータが算出される可能性が極めて高い。条件2は造影剤による輝度変化量に比べて、心拍や体動による輝度変化量が小さいことを意味する。条件1を満たさない場合、ΔIは負になり条件2を満たすことはない。条件1を満たし、かつ条件2を満たさない画素について説明する。本実施例で用いる画像データは造影剤の信号を強調し、かつ組織からの信号を抑圧した画像を基本とする。その上で条件2を満たさない画素は、造影効果がない、又は極めて低い画素である。条件1を満たしていることから本画素のTICは造影段階にて上昇するものの、その要因は主に造影剤以外の体動等による組織成分の入り込みと考えられることから、非造影部と判断する。上記の条件1、条件2を満たさない画素は非造影部としてマッピングパラメータtα、及びΔIが0値化される。第1画素選択部15における非造影部の除去効果は比較的小さい。第1画素選択部15では確実に非造影部と考えられる画素を除去し、造影部と非造影部の判断が曖昧になる画素は可能な限り抽出して必要な情報を失わないようにする。着目血管の抽出は主に後段の第2画素選択部17で実施され、第1画素選択部15はその補助的な役割を担う。この仕組みにより、血流動態画像の初期画像の視認性が向上し、第2の画素の選択の工程である工程6で実施される術者による画像調整がし易くなる。逆に第1画素選択部15の機能がない場合には、非造影部の画素が持つtαの乱数により前記初期画像の色分布が雑然となる。特に微細な血管の視認性が低下し、第2画素抽出部17における画像調整が困難になるだけでなく、着目すべき血管を見落とす危険性も出てくる。続いて工程5にてそれぞれの画素で求められたマッピングパラメータtαの値に応じて色分けされた血流動態画像51が図5に示すように画像データ32と併せて表示画面34に表示される。ここのtαは、前述したように最大値Imaxに予め設定する係数αをかけた値Iαに達する時間である。tαを用いる理由は、それぞれの画素に対応する着目部位の造影剤の開始時刻の相対的な値を比べるためである。本来は、造影剤開始時刻はtstartであるが、tstartは、図9に示すように、同じ値の輝度値を持つ時間が複数あり、いつが開始時刻なのか決定するのが難しい。その点、計算で求めたtαは、いずれの画素においても、一点であると考えられるので、比較ができる。図5の説明に戻るが、tα0の画素は色づけされない。表示される血流動態画像の色の分布から、着目する組織の造影剤の流入開始時間を相対的に比較できる。その結果、動脈性又は門脈性といった着目病変の支配血流の性質が判り、鑑別診断に有用な情報が得られる。工程6では血流動態画像の描画範囲、色分け、血管抽出といった描画情報の調整が行われる。血流動態画像51にはマッピングパラメータの時間の数値と色との対応関係を示す時間カラーバー52が付属して表示され、時間カラーバーは処理の開始時間tpreを0とする相対的な時間で表示される。時間カラーバー52に設けられる始点と終点の位置を変更すると、図6に示すように時間カラーバー52の最小値と最大値が始点と終点に変更される。同時に、時間と配色の対応も変更される。変更された結果が、対応するそれぞれの画素の配色に反映され、血流動態画像51が更新される。設定範囲から外れる画素のtαは0値化され、血流動態画像から色が除かれる。さらに図6に示すように分割時間を変更することで色分けする時間の区切りを変更することができ、この変更も血流動態画像に随時反映される。始点と終点の設定は表示画面34に表示される数値を直接変更することでも実施できる。表示画面34に設けられた計測ボタン54を選択すると血流動態画像51上に第1のROI55が表示され、始点に対する絶対的な流入時間が数字で表示される。再度計測ボタンを選択すると第2のROI56が表示される。第2のROI56については、第1のROI55に対する相対的な流入時間差が数字で表示される。相対的であることを示すために、第1のROI55から第2のROI56に対して矢印が引かれる。画面中の矢印を選択すると、第2のROI56も始点に対して絶対的な流入時間差が表示され、矢印が消失する。この計測機能により、二つの血管の血流動態の相対的な違いを定量的に評価できる。以上のように、血流動態画像の時間範囲と配色を自由に調整することで、術者が着目する時間範囲の血管を正確に抽出し、さらに観察しやすい形に画像を調整し表示することができる。時間レンジに関する変更は、術者の着目範囲の選別と配色を調整するものである。したがって、この段階では血流動態画像には工程4で除去しきれなかった非造影部の画素を主とするノイズ成分が残留しており、特に腫瘍血管のような比較的細い血管はノイズに埋もれる場合がある。造影部と非造影部の区別は時間情報よりも輝度情報を利用した方が効果的であるため、ΔIに基づく画像調整機能により実施する。図5に示すように表示画像上に設けられた画素選択ボタン53を術者が押すと、第2画素選択部及び画像構成部が機能をし、図7に示すように横軸ΔI、縦軸が画素数のヒストグラムが表示される。ヒストグラムにはΔI=0は含めない。第1画素選択部で条件を満たさない画素は0値化しており、この画素数は他の値の画素数に比べて非常に大きい値となる可能性が高い。更にこの値は血流動態画像を調整する上で不用であるため、ヒストグラムの視認性を高めるためにΔI>0の範囲で表示する。ヒストグラムには更に描画範囲を示す矢印又は色分けされた領域、又はその両方が表示される。術者がヒストグラム上でΔIの描画範囲を指定すると、第2画素選択部が、その範囲外のΔIを持つ画素のマッピングパラメータが即時0値化し、血流動態画像から色が除かれる。非造影部に相当する画素のΔIは造影部に極めて小さいため、この範囲を描画範囲から外すことで効果的にノイズを除去することができる。また、太い血管ほどΔIが大きい傾向もあるため、血管径に応じた選択的な抽出も可能である。描画範囲を指定する際に、図8に示すように、ヒストグラムをΔIの値に応じて例えば図のように3分割し、レンジバーに変換して表示することも可能である。レンジバーの各区間の長さは、分割された各範囲に含まれる画素の総数の比率で決まる。この時にもΔI0の画素数は含めない。レンジバーの軸にはΔIの大小をそのまま表示してもよく、またΔIが小さい方から大きい方にかけて微細血管、大血管と表示することでより判りやすくできる。更に図8に示すように、ΔIの値に応じて濃淡又は色づけしたΔI画像をΔIのレンジバーやヒストグラムと併せて表示し、描画範囲の指定に伴い除去される画素を画像上で確認することで、より正確な調整が可能である。しかしながら、残したい画素が除去されるケース、又はその逆のケースは発生しうる。この問題を解決する機能として、図8に示すようにΔI画像には除去ボタン、抽出ボタン、LPFボタンが設けられている。例えば術者が抽出ボタンを選択した後、常に抽出したい画素又は領域を選択すると、その領域は描画範囲に関わらず常に抽出される。逆に除去ボタンにて選択された画素は常に除去される。LPFボタンを選択すると、ΔI画像に対してLPFが掛けられる。ΔIの分布が空間的に平滑化されるため、描画範囲の変更に伴う画素の抽出及び除去がある程度まとまった範囲で実行される。常に除去する画素を解除する場合、再度除去ボタンを選択して当該画素を選択すれば解除が実行される。抽出画素に関する除去も同様である。LPFに関しては、再度LPFボタンを選択すれば解除される。以上に説明したように、本実施例の超音波診断装置は、被検体1に超音波を送受信する探触子2と、探触子2が受信した超音波信号を整相加算する受信ビームフォーマ部7と、被検体1に造影剤が投与された開始時間を含む情報の指定を受け付ける処理画像選択部11と、開始時間の情報に基づいて、整相加算された信号に基づく画像データから、前記造影剤で造影された部分を抽出する第1画素選択部15と、第1画素選択部15で抽出された画像データの画素情報に基づいて、更に画像データを抽出する第2画素選択部17と、第2画素選択部17で抽出された画像データに基づいて血流動態画像を生成する画像構成部16と、血流動態画像を表示する表示部10とからなる。2段階の画素抽出機能を持ち、第1段階で非造影部を緩い条件下で自動的に除去し、第2段階で術者による時間及び輝度に基づく画像調整が為される。なお、輝度に基づく画像調整の前段階として、時間に基づく調整を行うと、画像の視認性を向上することができ、時間に基づく画像調整が容易になる。また、第1の画像抽出の段階は、緩い条件下にすることにより、ノイズを除去する。ノイズ除去の条件を強くしすぎると、例えば、細い血管の信号までをも除去することになってしまうので、本発明では、例えば、工程4で定めたような、緩い条件とする。輝度情報ΔIに基づく調整では、着目する組織以外の成分を除去することを主目的とする。第1段階での非造影部の除去処理との違いは、着目組織と除去効果を視覚的に観察しながら調整できる点である。ΔIに基づく処理は、画像化しない情報に基づく画像調整であることが重要な点である。例えば、太い血管からの信号を除去し、細い血管を抽出したい場合には、ΔIを図8に示すように、設定することができる。時間tαに基づく調整では、着目する時間範囲を絞り込むことで画像の視認性を向上し、更に単位時間当たりの色数を増やすことで血流動態の細かい違いを大きな色の違いで識別できる。この処理は画像化している情報に基づく画像調整であり、本質的には着目する血流動態の着目範囲及び配色に対する処理である。工程6の説明では時間情報による調整をしてから輝度情報による調整を実施しているが、実際には順番が逆になることもあり、又は交互に実施されることもある。両者の調整は互いに密接に関連するため、着目血管の情報を正確に抽出するには、両者のバランスが重要となる。以上の工程により血流動態画像が構成され、表示部10に表示される工程7。表示画面34に表示される画像は、造影画像と血流動態画像の他、血流動態画像の色を透明化し、造影画像上に重畳表示してもよい。表示画面34には設定変更のボタンが設けられており、これを選択すると図9に示すように血流動態画像の処理工程の順に造影前時間の初期設定と各種パラメータが表示される。また、図4Aに示すTICと各パラメータの関係図も併せて表示される。造影時間や閾値輝度の係数αを変更する場合、表示画面の該当する数字を直接変更する。変更した結果は即時TICの関係図に反映される。ΔIや画素選択の条件式を変更する場合、数式を選択して超音波診断機に付属するキーボードから直接入力してもよく、又は数式上の変更したいパラメータを選択した後、変更後のパラメータ上にドラッグすることでパラメータの変更が実行される。マッピングパラメータはTIC上に丸印で示される。マッピングパラメータを変更する場合、キーボードで直接入力してもよく、またTIC上の丸印を新たにマッピングパラメータとして設定するパラメータ上にドラッグすることで変更が実行される。パラメータのドラッグは一般的な超音波診断装置に付属するリターンボタン及びトラックボール又はそれに変わる機能をもつ装置で実施される。タッチパネルを供えた装置であれば、画面上を指でなぞることで実施される。更に設定変更の画面には初期値減算のONOFFも選択できる。これはTIC評価部13における各種パラメータを算出する際に、各画素のTICの全値から初期値を減算し、TICの初期値を0値化する処理である。この処理の長所について説明する。造影前の段階で輝度値が高い画素は、Imaxが高い値となるためこれに係数を乗じたIαはImaxから値を大きく下げる。一方、TICから初期値を減算するとImaxの値が下がり、同じ係数を乗じた場合でもImaxに対するIαの下げ幅は小さくなる。そのため、初期値減算をONにすることで、輝度第1画素抽出部15の条件1をより満たしやすくなる。この効果は、主に組織成分が強く残留する造影部の画素を積極的に抽出することにある。本実施例では、血管部を選択的に抽出することに重点を置き、更に血管内部は超音波画像上では基本的に低輝度に描出されることから、初期設定はOFFとしている。各種パラメータはパラメータメモリ14に保存してあるため、設定変更と同時に画像の更新が実施される。但し、係数αの変更に関してはマッピングパラメータtαをTICから再度算出する必要があるため、TIC評価部13にて再度パラメータの算出を実行する。但し、予め特定の係数例えばα05060708のマッピングパラメータ値を算出してパラメータメモリ14に保持することで、TIC評価部での再計算を避けることが可能である。処理範囲を変更した場合も同様にTIC評価部13での再計算が必要となる。次に、構成した血流動態画像を過去の血流動態画像と比較する場合の機能について説明する。この機能は血管を標的とした治療の効果判定に有効性を発揮するため、治療前後の血流動態画像を比較することを例に説明する。表示画面34に設けられている画像読み出しボタン100を選択すると、超音波診断装置内に設けられているHDDに保存されている過去の血流動態画像101が表示画面に表示される図10。ここで示す血流動態画像は、造影剤の流入開始時間を例として説明を続ける。治療前後の血流動態画像では時間レンジやΔIに基づく輝度レンジの調整が異なるため、時間カラーバー105とΔIに基づく輝度レンジバー104も併せて表示される。ここでは、過去の血流動態画像を治療前の血流動態画像101とし、既に表示されている血流動態画像51を治療後のものとする。表示画面34に表示される基準設定ボタン107を選択すると、治療前の血流動態画像101及び治療後の血流動態画像51に基準となる領域を設定するための基準ROIが治療前後の血流動態画像上に表示され造影前の基準ROI102、造影後の基準ROI103、同時に時間カラーバー105、時間カラーバー52に基準ROIで指定する領域基準領域の平均時間の位置が矢印で示される。次に、ROIを治療の影響を受けない例えば正常血管の位置に設定する。治療前後の血流動態画像から同じ血管を選択できればよいが、それが困難な場合には、着目する領域に位置関係が近い血管を選択する。血流動態画像の色の違いから各血管の造影開始時間を容易に識別できるため、血管の空間的な位置関係だけでなく、血流動態の類似性を利用して、基準ROI102及び基準ROI103の設定が可能である。次に、基準設定ボタン107を再度選択すると、基準領域の色、つまり流入開始時間の違いから、治療前後の血流動態画像の時間レンジの相対的な時間差が算出され、治療前の血流動態画像の流入開始時間と色の対応関係が、治療後の血流動態画像に一致するように調整される。例えば図11に示すように治療前後で基準領域の時間にtoffsettposttpreの時間差があった場合、治療前の血流動態画像を構成する全画素のマッピングパラメータにtoffsetが加算され、新たなマッピングパラメータTstarttstarttoffsetTpretpretoffsetが再設定される。この補正処理により基準領域の配色、つまりtpostとTpreの配色が治療前後で一致し、病変部の血流動態を色の違いで容易に比較することができる。治療前の血流動態画像に対する色の調整処理の結果は即時に血流動態画像に反映される。そのため術者は画面上で調整結果を確認しながら、基準領域を最適な位置に再設定できる。なお、補正によりカラーバーから外れる時間に関しては、血流動態画像上において無色になる。輝度レンジに関しても、時間レンジと同様の手法で補正される。また、表示される治療前後の画像は、血流動態画像を半透明にした状態で造影画像を重畳した画像でもよい。以上に説明の実施例においては、スキャンコンバータ9からの画像データを利用して解析を行なっているが、受信ビームフォーマ部7から出力されるRAWデータを利用することも可能である。本実施例の画像データは、スキャンコンバータ9からの画像データだけでなく、RAWデータを含んで解釈することができる。図12に示すように、スキャンコンバータ9からの画像データに加え、受信ビームフォーマ部7からのRAWデータを画像メモリ12に蓄積する。図2に示すフローチャートにおいて、画像データの選択工程1は画像データを利用して実施し、その後の数値解析工程2から工程6をRAWデータで実施する。RAWデータは生体内の情報を直接的に反映しており、更に画像データに比べて一般的に深いbit数の信号である。そのため画像データを利用する場合に比べて高感度でかつ高精度な数値解析が可能となり、例えば部位毎のマッピングパラメータ値の違いをより顕著にしたり、ΔIが僅かに異なる血管を区別し易くしたりできる。