JP6222829B2 - 超音波診断装置、画像処理装置及び画像処理方法 - Google Patents

Description

本発明の実施形態は、超音波診断装置、画像処理装置及び画像処理方法に関する。

近年、静脈投与型の超音波造影剤が製品化され、「造影エコー法」が行なわれている。以下、超音波造影剤を省略して造影剤と記載する場合がある。造影エコー法は、例えば、心臓や肝臓等の検査において、静脈から造影剤を注入して血流信号を増強し、血流動態の評価を行うことを目的としている。造影剤の多くは、微小気泡（マイクロバブル）が反射源として機能する。例えば、近年、日本で発売されたソナゾイド（登録商標）と呼ばれる第二世代の超音波造影剤は、リン脂質によりフッ化炭素（perfluorobutane）ガスを内包した微小気泡である。造影エコー法では、微小気泡を破壊させない程度の中低音圧の送信超音波を用いることで、造影剤の還流の様子を安定して観察することができる。

上記造影剤の投与後の診断部位（例えば、肝臓がん）を超音波走査すると、医師等の操作者は、血流により還流する造影剤の流入から流出までの信号強度の上昇と減少とを観察することができる。また、かかる信号強度の経時変化の違いから、腫瘤性病変の良性／悪性の鑑別診断、或いは、「びまん性」の疾患等の診断を行なう研究も行なわれている。

造影剤の還流動態を示す信号強度の経時変化は、単なる形態情報とは異なり、通常、動画像をリアルタイムで、或いは、記録後に、読影される必要がある。従って、造影剤の還流動態の読影のための所要時間は、一般的に長くなる。そこで、通常は動画で観察される造影剤の流入時刻情報を、一枚の静止画像上にマッピングする手法が提案されている。かかる手法は、造影剤の信号のピーク時間の違いを、異なる色相で表現した静止画像を生成表示するものである。この静止画像を参照することで、読影者は、診断部位の断層面内の各所での流入時刻を容易に把握することができる。また、特定領域で造影剤が停滞する時間（流入開始から流出終了までの時間）の違いを、異なる色相で表現した静止画像を生成表示する手法も提案されている。

ところで、腫瘍血管の走行は、正常な血管に比べて複雑であるため、行き場を失ったマイクロバブルが腫瘍に停滞したり、更に、停滞したマイクロバブルが逆流したりする現象が観察されている。このような腫瘍血管内でのマイクロバブルの挙動は、造影超音波撮影を行った腫瘍マウスで実際に観察されている。すなわち、生体イメージングが可能な造影超音波撮影でマイクロバブルの挙動を評価できれば、造影エコー法は、腫瘍血管の異常の評価にも応用できる可能性がある。

また、近年、治験が行なわれている抗がん剤である血管新生阻害剤は、腫瘍を栄養する血管を破壊して、腫瘍血管の断片化や狭小化を起こすことが、病理組織学的な観察で確認されている。造影超音波撮影で血管新生阻害剤により断片化された血管内で造影超音波観察が停滞する様子を映像化したり定量化したりできれば、造影エコー法は、治療効果判定に応用することも期待できる。

しかし、信号強度の変化、すなわち、超音波画像の輝度の変化は、撮影条件や、観測領域よって変化する。例えば、輝度の変化は、造影剤の種類や、観察領域内の血管の性質及び血管周囲の組織の性質により変化する。一方、上記の静止画像は、撮影条件や観測領域に関わらず、観察される絶対的な特徴量（例えば、絶対時間や絶対輝度）に基づいて、造影剤流入時刻が決定され、造影剤流入時刻に基づいて信号強度の経時変化の解析が行なわれることで、生成表示される。

特開２００１−２６９３４１号公報 特開２００２−２３８９０１号公報 特開２０１１−２５４９６３号公報

本発明が解決しようとする課題は、造影剤の還流動態を客観的な基準により解析することができる超音波診断装置、画像処理装置及び画像処理方法を提供することである。

実施形態の超音波診断装置は、輝度変化情報生成部と、解析部と、制御部とを備える。輝度変化情報生成部は、造影剤が投与された被検体を超音波走査して収集された時系列データから、超音波走査領域内に設定された解析領域における輝度の時間変化を示す輝度変化情報を生成する。解析部は、前記輝度変化情報に基づいて、前記解析領域における造影剤の還流動態を時間について正規化したパラメータを取得する。制御部は、前記パラメータを画像又は文字のうち少なくともいずれか一方の形式で表示部に表示させる。

図１は、本実施形態に係る超音波診断装置の構成例を示すブロック図である。 図２は、解析領域の一例を示す図（１）である。 図３は、解析領域の一例を示す図（２）である。 図４は、解析領域の一例を示す図（３）である。 図５は、解析部を説明するための図（１）である。 図６は、解析部を説明するための図（２）である。 図７は、解析部を説明するための図（３）である。 図８は、解析部を説明するための図（４）である。 図９は、変化画像生成部を説明するための図（１）である。 図１０は、変化画像生成部を説明するための図（２）である。 図１１は、変化画像生成部を説明するための図（３）である。 図１２は、本実施形態に係る超音波診断装置が処理の一例を示すフローチャートである。 図１３は、本実施形態に係る変形例を説明するための図（１）である。 図１４は、本実施形態に係る変形例を説明するための図（２）である。

以下、添付図面を参照して、超音波診断装置の実施形態を詳細に説明する。

（実施形態）
まず、本実施形態に係る超音波診断装置の構成について説明する。図１は、本実施形態に係る超音波診断装置の構成例を示すブロック図である。図１に例示するように、第１の実施形態に係る超音波診断装置は、超音波プローブ１と、モニタ２と、入力装置３と、装置本体１０とを有する。

超音波プローブ１は、複数の圧電振動子を有し、これら複数の圧電振動子は、後述する装置本体１０が有する送受信部１１から供給される駆動信号に基づき超音波を発生する。また、超音波プローブ１は、被検体Ｐからの反射波を受信して電気信号に変換する。また、超音波プローブ１は、圧電振動子に設けられる整合層と、圧電振動子から後方への超音波の伝播を防止するバッキング材等を有する。なお、超音波プローブ１は、装置本体１０と着脱自在に接続される。

超音波プローブ１から被検体Ｐに超音波が送信されると、送信された超音波は、被検体Ｐの体内組織における音響インピーダンスの不連続面で次々と反射され、反射波信号として超音波プローブ１が有する複数の圧電振動子にて受信される。受信される反射波信号の振幅は、超音波が反射される不連続面における音響インピーダンスの差に依存する。なお、送信された超音波パルスが、移動している血流や心臓壁等の表面で反射された場合の反射波信号は、ドプラ効果により、移動体の超音波送信方向に対する速度成分に依存して、周波数偏移を受ける。

例えば、装置本体１０は、２次元走査用の超音波プローブ１として、複数の圧電振動子が一列で配置された１Ｄアレイプローブと接続される。或いは、例えば、装置本体１０は、３次元走査用の超音波プローブ１として、メカニカル４Ｄプローブや２Ｄアレイプローブと接続される。メカニカル４Ｄプローブは、１Ｄアレイプローブのように一列で配列された複数の圧電振動子を用いて２次元走査が可能であるとともに、複数の圧電振動子を所定の角度（揺動角度）で揺動させることで３次元走査が可能である。また、２Ｄアレイプローブは、マトリックス状に配置された複数の圧電振動子により３次元走査が可能であるとともに、超音波を集束して送信することで２次元走査が可能である。

本実施形態は、超音波プローブ１により、被検体Ｐが２次元走査される場合であっても、３次元走査される場合であっても適用可能である。

入力装置３は、マウス、キーボード、ボタン、パネルスイッチ、タッチコマンドスクリーン、フットスイッチ、トラックボール、ジョイスティック等を有し、超音波診断装置の操作者からの各種設定要求を受け付け、装置本体１０に対して受け付けた各種設定要求を転送する。例えば、入力装置３は、操作者から超音波造影剤の還流動態を解析するための解析領域の設定を受け付ける。なお、本実施形態で設定される解析領域については、後に詳述する。

モニタ２は、超音波診断装置の操作者が入力装置３を用いて各種設定要求を入力するためのＧＵＩ（Graphical User Interface）を表示したり、装置本体１０において生成された超音波画像等を表示したりする。

装置本体１０は、超音波プローブ１が受信した反射波信号に基づいて超音波画像データを生成する装置である。図１に示す装置本体１０は、超音波プローブ１が受信した２次元の反射波データに基づいて２次元の超音波画像データを生成可能な装置である。また、図１に示す装置本体１０は、超音波プローブ１が受信した３次元の反射波データに基づいて３次元の超音波画像データを生成可能な装置である。以下、３次元の超音波画像データを「ボリュームデータ」と記載する場合がある。

装置本体１０は、図１に示すように、送受信部１１と、Ｂモード処理部１２と、ドプラ処理部１３と、画像生成部１４と、画像処理部１５と、画像メモリ１６と、内部記憶部１７と、制御部１８とを有する。

送受信部１１は、パルス発生器、送信遅延部、パルサ等を有し、超音波プローブ１に駆動信号を供給する。パルス発生器は、所定のレート周波数で、送信超音波を形成するためのレートパルスを繰り返し発生する。また、送信遅延部は、超音波プローブ１から発生される超音波をビーム状に集束し、かつ送信指向性を決定するために必要な圧電振動子ごとの遅延時間を、パルス発生器が発生する各レートパルスに対し与える。また、パルサは、レートパルスに基づくタイミングで、超音波プローブ１に駆動信号（駆動パルス）を印加する。すなわち、送信遅延部は、各レートパルスに対し与える遅延時間を変化させることで、圧電振動子面から送信される超音波の送信方向を任意に調整する。

なお、送受信部１１は、後述する制御部１８の指示に基づいて、所定のスキャンシーケンスを実行するために、送信周波数、送信駆動電圧等を瞬時に変更可能な機能を有している。特に、送信駆動電圧の変更は、瞬間にその値を切り替え可能なリニアアンプ型の発信回路、又は、複数の電源ユニットを電気的に切り替える機構によって実現される。

また、送受信部１１は、プリアンプ、Ａ／Ｄ（Analog／Digital）変換器、受信遅延部、加算器等を有し、超音波プローブ１が受信した反射波信号に対して各種処理を行って反射波データを生成する。プリアンプは、反射波信号をチャネルごとに増幅する。Ａ／Ｄ変換器は、増幅された反射波信号をＡ／Ｄ変換する。受信遅延部は、受信指向性を決定するために必要な遅延時間を与える。加算器は、受信遅延部によって処理された反射波信号の加算処理を行なって反射波データを生成する。加算器の加算処理により、反射波信号の受信指向性に応じた方向からの反射成分が強調され、受信指向性と送信指向性とにより超音波送受信の総合的なビームが形成される。

送受信部１１は、被検体Ｐを２次元走査する場合、超音波プローブ１から２次元の超音波ビームを送信させる。そして、送受信部１１は、超音波プローブ１が受信した２次元の反射波信号から２次元の反射波データを生成する。また、送受信部１１は、被検体Ｐを３次元走査する場合、超音波プローブ１から３次元の超音波ビームを送信させる。そして、送受信部１１は、超音波プローブ１が受信した３次元の反射波信号から３次元の反射波データを生成する。

なお、送受信部１１からの出力信号の形態は、ＲＦ（Radio Frequency）信号と呼ばれる位相情報が含まれる信号である場合や、包絡線検波処理後の振幅情報である場合等、種々の形態が選択可能である。

Ｂモード処理部１２は、送受信部１１から反射波データを受信し、対数増幅、包絡線検波処理等を行なって、信号強度が輝度の明るさで表現されるデータ（Ｂモードデータ）を生成する。

また、Ｂモード処理部１２は、フィルタ処理により、検波周波数を変化させることで、映像化する周波数帯域を変えることができる。このＢモード処理部１２の機能を用いることにより、造影エコー法、例えば、コントラストハーモニックイメージング（ＣＨＩ：Contrast Harmonic Imaging）を実行可能である。すなわち、Ｂモード処理部１２は、超音波造影剤が注入された被検体Ｐの反射波データから、微小気泡（マイクロバブル）を反射源とする反射波データ（高調波データ又は分周波データ）と、被検体Ｐ内の組織を反射源とする反射波データ（基本波データ）とを分離することができる。これにより、Ｂモード処理部１２は、被検体Ｐの反射波データから高調波データ又は分周波データを抽出して、造影画像データを生成するためのＢモードデータを生成することができる。造影画像データを生成するためのＢモードデータは、造影剤を反射源とする反射波の信号強度を輝度で表わしたデータとなる。また、Ｂモード処理部１２は、被検体Ｐの反射波データから基本波データを抽出して、組織画像データを生成するためのＢモードデータを生成することができる。

なお、ＣＨＩを行なう際、Ｂモード処理部１２は、上述したフィルタ処理を用いた方法とは異なる方法により、ハーモニック成分（高調波成分）を抽出することができる。ハーモニックイメージングでは、振幅変調（ＡＭ：Amplitude Modulation）法や位相変調（ＰＭ：Phase Modulation）法、ＡＭ法及びＰＭ法を組み合わせたＡＭＰＭ法と呼ばれる映像法が行なわれる。ＡＭ法、ＰＭ法及びＡＭＰＭ法では、同一の走査線に対して振幅や位相が異なる超音波送信を複数回（複数レート）行なう。これにより、送受信部１１は、各走査線で複数の反射波データを生成し出力する。そして、Ｂモード処理部１２は、各走査線の複数の反射波データを、変調法に応じた加減算処理することで、高調波成分を抽出する。そして、Ｂモード処理部１２は、高調波成分の反射波データに対して包絡線検波処理等を行なって、Ｂモードデータを生成する。

例えば、ＰＭ法が行なわれる場合、送受信部１１は、制御部１８が設定したスキャンシーケンスにより、例えば（−１，１）のように、位相極性を反転させた同一振幅の超音波を、各走査線で２回送信させる。そして、送受信部１１は、「−１」の送信による反射波データと、「１」の送信による反射波データとを生成し、Ｂモード処理部１２は、これら２つの反射波データを加算する。これにより、基本波成分が除去され、２次高調波成分が主に残存した信号が生成される。そして、Ｂモード処理部１２は、この信号に対して包絡線検波処理等を行なって、ＣＨＩのＢモードデータ（造影画像データを生成するためのＢモードデータ）を生成する。ＣＨＩのＢモードデータは、造影剤を反射源とする反射波の信号強度を輝度で表わしたデータとなる。また、ＣＨＩでＰＭ法が行なわれる場合、Ｂモード処理部１２は、例えば、「１」の送信による反射波データをフィルタ処理することで、組織画像データを生成するためのＢモードデータを生成することができる。

ドプラ処理部１３は、送受信部１１から受信した反射波データから速度情報を周波数解析し、ドプラ効果による血流や組織、造影剤エコー成分を抽出し、速度、分散、パワー等の移動体情報を多点について抽出したデータ（ドプラデータ）を生成する。

なお、本実施形態に係るＢモード処理部１２及びドプラ処理部１３は、２次元の反射波データ及び３次元の反射波データの両方について処理可能である。すなわち、Ｂモード処理部１２は、２次元の反射波データから２次元のＢモードデータを生成し、３次元の反射波データから３次元のＢモードデータを生成する。また、ドプラ処理部１３は、２次元の反射波データから２次元のドプラデータを生成し、３次元の反射波データから３次元のドプラデータを生成する。

画像生成部１４は、Ｂモード処理部１２及びドプラ処理部１３が生成したデータから超音波画像データを生成する。すなわち、画像生成部１４は、Ｂモード処理部１２が生成した２次元のＢモードデータから反射波の強度を輝度で表した２次元Ｂモード画像データを生成する。また、画像生成部１４は、ドプラ処理部１３が生成した２次元のドプラデータから移動体情報を表す２次元ドプラ画像データを生成する。２次元ドプラ画像データは、速度画像、分散画像、パワー画像、又は、これらを組み合わせた画像である。

ここで、画像生成部１４は、一般的には、超音波走査の走査線信号列を、テレビ等に代表されるビデオフォーマットの走査線信号列に変換（スキャンコンバート）し、表示用の超音波画像データを生成する。具体的には、画像生成部１４は、超音波プローブ１による超音波の走査形態に応じて座標変換を行なうことで、表示用の超音波画像データを生成する。また、画像生成部１４は、スキャンコンバート以外に種々の画像処理として、例えば、スキャンコンバート後の複数の画像フレームを用いて、輝度の平均値画像を再生成する画像処理（平滑化処理）や、画像内で微分フィルタを用いる画像処理（エッジ強調処理）等を行なう。また、画像生成部１４は、超音波画像データに、種々のパラメータの文字情報、目盛り、ボディーマーク等を合成する。

すなわち、Ｂモードデータ及びドプラデータは、スキャンコンバート処理前の超音波画像データであり、画像生成部１４が生成するデータは、スキャンコンバート処理後の表示用の超音波画像データである。なお、Ｂモードデータ及びドプラデータは、生データ（Raw Data）とも呼ばれる。

更に、画像生成部１４は、Ｂモード処理部１２が生成した３次元のＢモードデータに対して座標変換を行なうことで、３次元Ｂモード画像データを生成する。また、画像生成部１４は、ドプラ処理部１３が生成した３次元のドプラデータに対して座標変換を行なうことで、３次元ドプラ画像データを生成する。すなわち、画像生成部１４は、「３次元のＢモード画像データや３次元ドプラ画像データ」を「３次元超音波画像データ（ボリュームデータ）」として生成する。

更に、画像生成部１４は、ボリュームデータをモニタ２にて表示するための各種の２次元画像データを生成するために、ボリュームデータに対してレンダリング処理を行なう。画像生成部１４が行なうレンダリング処理としては、断面再構成法（ＭＰＲ：Multi Planer Reconstruction）を行なってボリュームデータからＭＰＲ画像データを生成する処理がある。また、画像生成部１４が行なうレンダリング処理としては、ボリュームデータに対して「Curved MPR」を行なう処理や、ボリュームデータに対して「Maximum Intensity Projection」を行なう処理がある。また、画像生成部１４が行なうレンダリング処理としては、３次元の情報を反映した２次元画像データを生成するボリュームレンダリング（ＶＲ：Volume Rendering）処理がある。

画像メモリ１６は、画像生成部１４が生成した表示用の画像データを記憶するメモリである。また、画像メモリ１６は、Ｂモード処理部１２やドプラ処理部１３が生成したデータを記憶することも可能である。画像メモリ１６が記憶する表示用の画像データは、例えば、診断の後に操作者が呼び出すことが可能である。また、画像メモリ１６が記憶するＢモードデータやドプラデータも、例えば、診断の後に操作者が呼び出すことが可能となっており、画像生成部１４を経由して表示用の超音波画像データとなる。また、画像メモリ１６は、送受信部１１から出力されたデータを記憶することも可能である。

画像処理部１５は、コンピュータ支援診断（Computer-Aided Diagnosis：ＣＡＤ）を行なうために、装置本体１０に設置される。画像処理部１５は、画像メモリ１６に格納されたデータを取得して、診断支援のための画像処理を行なう。そして、画像処理部１５は、画像処理結果を、画像メモリ１６や後述する内部記憶部１７に格納する。なお、画像処理部１５が行なう処理については、後に詳述する。

内部記憶部１７は、超音波送受信、画像処理及び表示処理を行なうための制御プログラムや、診断情報（例えば、患者ＩＤ、医師の所見等）や、診断プロトコルや各種ボディーマーク等の各種データを記憶する。また、内部記憶部１７は、必要に応じて、画像メモリ１６が記憶する画像データの保管等にも使用される。また、内部記憶部１７が記憶するデータは、図示しないインターフェースを経由して、外部装置へ転送することができる。なお、外部装置は、例えば、各種医用画像診断装置、画像診断を行なう医師が使用するＰＣ（Personal Computer）や、ＣＤやＤＶＤ等の記憶媒体、プリンター等である。

制御部１８は、超音波診断装置の処理全体を制御する。具体的には、制御部１８は、入力装置３を介して操作者から入力された各種設定要求や、内部記憶部１７から読込んだ各種制御プログラム及び各種データに基づき、送受信部１１、Ｂモード処理部１２、ドプラ処理部１３、画像生成部１４及び画像処理部１５の処理を制御する。また、制御部１８は、画像メモリ１６や内部記憶部１７が記憶する画像データをモニタ２にて表示するように制御する。

以上、本実施形態に係る超音波診断装置の全体構成について説明した。かかる構成のもと、本実施形態に係る超音波診断装置は、造影剤の還流動態の解析を目的とする造影エコー法を行なう。そして、本実施形態に係る超音波診断装置は、超音波造影剤が投与された被検体Ｐを超音波走査して収集された時系列データから、超音波走査領域内に設定された解析領域における造影剤の還流動態を客観的な基準により解析することができる画像データを生成し表示する。

かかる画像データを生成するために、本実施形態に係る画像処理部１５は、図１に示すように、輝度変化情報生成部１５１と、解析部１５２と、変化画像生成部１５３とを有する。

図１に示す輝度変化情報生成部１５１は、造影剤が投与された被検体Ｐを超音波走査して収集された時系列データから、超音波走査領域内に設定された解析領域における輝度の時間変化を示す輝度変化情報を生成する。具体的には、輝度変化情報生成部１５１は、輝度変化情報として、解析領域における輝度の時間変化を示す曲線である輝度変化曲線を生成する。なお、輝度変化情報生成部１５１は、輝度変化曲線を再現可能な情報であるならば、任意の形態で、輝度変化情報を生成可能である。上記の時系列データは、造影期間中に画像生成部１４が時系列に沿って生成した複数の２次元又は３次元の造影画像データである。或いは、上記の時系列データは、造影期間中にＢモード処理部１２が時系列に沿って抽出した複数の２次元又は３次元の高周波データ（高調波成分）である。或いは、上記の時系列データは、造影期間中にＢモード処理部１２が時系列に沿って造影画像データ用に生成した複数の２次元又は３次元のＢモードデータである。

すなわち、造影撮影が２次元の超音波走査領域で行なわれる場合、輝度変化情報生成部１５１は、被検体Ｐを２次元走査して収集された時系列データから、２次元走査領域内に設定された２次元の解析領域における輝度変化曲線を生成する。また、造影撮影が３次元の超音波走査領域で行なわれる場合、輝度変化情報生成部１５１は、被検体Ｐを３次元走査して収集された時系列データから、３次元走査領域内に設定された３次元の解析領域、又は、２次元の解析領域における輝度変化曲線を生成する。

以下では、輝度変化情報生成部１５１が、被検体Ｐを２次元走査して時系列に沿って収集された複数の造影画像データから、２次元走査領域内に設定された２次元の解析領域における輝度変化曲線を生成する場合について説明する。

ここで、本実施形態に係る輝度変化情報生成部１５１は、複数の輝度変化曲線を生成する。例えば、輝度変化情報生成部１５１は、超音波走査領域内に設定された複数の解析領域それぞれにおける複数の輝度変化曲線を生成する。或いは、輝度変化情報生成部１５１は、異なる複数の時期に同一の超音波走査領域で行なわれた超音波走査により収集された複数の時系列データそれぞれから該領域内に設定された少なくとも１つの同一の解析領域それぞれにおける複数の輝度変化曲線を生成する。図２、図３及び図４は、解析領域の一例を示す図である。なお、以下の説明において、超音波プローブ１の位置は、解析領域の設定前後で、同一の位置に固定される。

例えば、操作者は、図２に示すように、造影前のＢモード画像データ（組織画像データ）に描出された肝臓の腫瘍部位に解析領域１００を設定し、肝臓の門脈に解析領域１０１を設定し、腎臓内に解析領域１０２を設定する。解析領域１０１は、腫瘍部位を還流する血流の動態を、肝臓全体を還流する血流の動態と比較するために設定される。また、通常、腎臓が造影剤により染影された後に肝臓は、染影される。このため、解析領域１０２は、肝臓全体を還流する血流の動態を、腎臓全体を還流する血流の動態と比較するために設定される。

輝度変化情報生成部１５１は、解析領域１００〜１０２の設定後、時系列に沿って収集された複数の造影画像データそれぞれから、解析領域１００内の平均輝度、解析領域１０１内の平均輝度及び解析領域１０２内の平均輝度を算出する。これにより、輝度変化情報生成部１５１は、３つの輝度変化曲線を生成する。

或いは、例えば、操作者は、図３の左図に示すように、血管新生阻害剤を用いた治療前に、造影前のＢモード画像データに解析領域１００を設定する。輝度変化情報生成部１５１は、解析領域１００の設定後、時系列に沿って収集された複数の造影画像データそれぞれから、解析領域１００内の平均輝度を算出して、解析領域１００の輝度変化曲線を生成する。

更に、例えば、操作者は、図３の右図に示すように、血管新生阻害剤を用いた治療後に、造影前のＢモード画像データに解析領域１００と同一の位置に解析領域１００’を設定する。輝度変化情報生成部１５１は、解析領域１００’の設定後、時系列に沿って収集された複数の造影画像データそれぞれから、解析領域１００’内の平均輝度を算出して、解析領域１００’の輝度変化曲線を生成する。解析領域１００の輝度変化曲線は、治療前の輝度変化曲線となり、解析領域１００’の輝度変化曲線は、治療後の輝度変化曲線となる。これにより、輝度変化情報生成部１５１は、２つの輝度変化曲線を生成する。

或いは、例えば、操作者は、図４の左図に示すように、造影前のＢモード画像データ（組織画像データ）に解析領域１００を設定して、最初に、造影剤Ａを用いた造影撮影を行なう。輝度変化情報生成部１５１は、解析領域１００の設定後、時系列に沿って収集された複数の造影画像データそれぞれから、解析領域１００内の平均輝度を算出して、解析領域１００における造影剤Ａの輝度変化曲線を生成する。

更に、例えば、所定期間（例えば、１０分）が経過した後、操作者は、図４の右図に示すように、造影剤Ａとは異なる種類の造影剤Ｂを用いた造影撮影を行なう。輝度変化情報生成部１５１は、造影剤Ｂの投与後、時系列に沿って収集された複数の造影画像データそれぞれから、解析領域１００内の平均輝度を算出して、解析領域１００における造影剤Ｂの輝度変化曲線を生成する。これにより、輝度変化情報生成部１５１は、２つの輝度変化曲線を生成する。

なお、図３及び図４では、収集時期が異なる２つの時系列データそれぞれから、１つの同一解析領域の輝度変化曲線が生成される場合について説明した。ただし、本実施形態は、収集時期が異なる２つの時系列データそれぞれから、複数の同一解析領域それぞれの輝度変化曲線が生成される場合であっても良い。また、本実施形態は、収集時期が異なる時系列データが３つ以上である場合であっても良い。

図１に示す解析部１５２は、輝度変化情報に基づいて、解析領域における造影剤の還流動態を時間について正規化したパラメータを取得する。ここで、解析部１５２は、解析領域における造影剤の還流動態を、輝度、又は、輝度及び時間について正規化したパラメータを取得可能である。本実施形態では、解析部１５２が、輝度変化情報に基づいて、解析領域における造影剤の還流動態を輝度及び時間について正規化したパラメータを取得する場合について説明する。換言すると、解析部１５２は、輝度変化曲線の形状を解析して、解析領域における造影剤の還流動態が正規化されたパラメータを取得する。具体的には、解析部１５２は、時間軸、又は、輝度軸及び時間軸の正規化を行なって輝度変化曲線から正規化曲線を生成する。本実施形態では、解析部１５２は、輝度軸及び時間軸の正規化を行なって輝度変化曲線から正規化曲線を生成する。例えば、正規化曲線を生成するために、解析部１５２は、輝度変化曲線において、輝度が最大値となる最大点と、最大点となる前に最大値に第１割合を乗算した第１乗算値となる第１点と、最大点となった後に最大値に第２割合を乗算した第２乗算値となる第２点とを取得する。第１割合及び第２割合は、初期的に設定されたり、操作者により事前に予め設定されたりする。また、第１割合及び第２割合は、操作者により任意に変更可能である。

以下、図２に例示した解析領域１００〜１０２の輝度変化曲線を用いて解析部１５２が行なう処理について、図５〜図８を用いて説明する。図５〜図８は、解析部を説明するための図である。

図５では、解析領域１００の輝度変化曲線を曲線Ｃ０（一点鎖線）とし、解析領域１０１の輝度変化曲線を曲線Ｃ１（二点鎖線）とし、解析領域１０２の輝度変化曲線を曲線Ｃ２（実線）として示している。なお、図５に例示する輝度変化曲線は、輝度変化情報生成部１５１が、数理モデルを用いて、解析領域内の平均輝度の時系列データから生成した近似曲線である。ここで、以下では、第１割合及び第２割合が共に「５０％」として設定されている場合について説明する。なお、本実施形態は、第１割合及び第２割合が異なる割合（例えば、２０％及び３０％）として設定される場合であっても良い。

解析部１５２は、図５に示すように、曲線Ｃ０を解析して、最大点「時間：ｔ０ｍａｘ、輝度：Ｉ０ｍａｘ」を取得する。そして、解析部１５２は、最大輝度の半分の値「Ｉ０ｍａｘ／２」を算出する。そして、解析部１５２は、図５に示すように、曲線Ｃ０において、最大時間以前に「Ｉ０ｍａｘ／２」となる第１点「時間：ｔ０ｓ、輝度：Ｉ０ｍａｘ／２」を取得する。そして、解析部１５２は、図５に示すように、曲線Ｃ０において、最大時間以後に「Ｉ０ｍａｘ／２」となる第２点「時間：ｔ０ｅ、輝度：Ｉ０ｍａｘ／２」を取得する。

同様の処理により、解析部１５２は、図５に示すように、曲線Ｃ１を解析して、最大点「時間：ｔ１ｍａｘ、輝度：Ｉ１ｍａｘ」、第１点「時間：ｔ１ｓ、輝度：Ｉ１ｍａｘ／２」、第２点「時間：ｔ１ｅ、輝度：Ｉ１ｍａｘ／２」を取得する。また、同様の処理により、解析部１５２は、図５に示すように、曲線Ｃ２を解析して、最大点「時間：ｔ２ｍａｘ、輝度：Ｉ２ｍａｘ」、第１点「時間：ｔ２ｓ、輝度：Ｉ２ｍａｘ／２」、第２点「時間：ｔ２ｅ、輝度：Ｉ２ｍａｘ／２」を取得する。

ここで、解析部１５２は、「最大点の時間」を、解析領域において造影剤が最大に流入した時間「最大時間」とする。また、解析部１５２は、「第１点の時間」を、解析領域において造影剤流入が開始した時間とみなし、この時間を血流動態の解析を開始する時間「開始時間」とする。すなわち、解析部１５２は、輝度変化曲線の時間軸の逆方向で、輝度の最大値から一定の割合（第１割合）まで減少するのに要する時間から、開始時間を設定する。換言すると、解析部１５２は、客観的な同一の基準（第１割合）を用いることで、解析対象となる輝度変化曲線の形状に応じた閾値（第１乗算値）を算出して、開始時間を設定する。かかる開始時間は、最大時間が決定された後、過去に遡って設定される時間、すなわち、「Retrospective」に設定される時間である。

また、解析部１５２は、「第２点の時間」を、解析領域から造影剤の流出が終了した時間とみなし、この時間を血流動態の解析を終了する時間「終了時間」とする。すなわち、解析部１５２は、輝度変化曲線の時間軸の順方向で、輝度の最大値から一定の割合（第２割合）まで減少するのに要する時間から、終了時間を設定する。換言すると、解析部１５２は、客観的な同一の基準（第２割合）を用いることで、解析対象となる輝度変化曲線の形状に応じた閾値（第２乗算値）を算出して、終了時間を設定する。かかる終了時間は、最大時間が決定された時点で予想される時間、すなわち、「Prospective」に設定される時間である。

そして、解析部１５２は、これら３つの点から選択された少なくとも２点を用いて、輝度変化曲線を正規化した正規化曲線を生成する。そして、本実施形態では、解析部１５２は、生成した正規化曲線から正規化されたパラメータを取得する。ここで、造影剤流入に関するパラメータを取得する場合、解析部１５２は、第１点及び最大点を用いて、正規化曲線を生成する。また、造影剤流出に関するパラメータを取得する場合、解析部１５２は、最大点及び第２点を用いて、正規化曲線を生成する。また、造影剤流入及び造影剤流出に関するパラメータを取得する場合、解析部１５２は、第１点、最大点及び第２点を用いて、正規化曲線を生成する。

本実施形態では、複数の輝度変化曲線が生成されることから、解析部１５２は、複数の輝度変化曲線それぞれから正規化曲線を生成する。そして、本実施形態では、解析部１５２は、生成した複数の正規化曲線それぞれからパラメータを取得する。以下、輝度軸及び時間軸の正規化により複数の輝度変化曲線それぞれから正規化曲線を生成する方法の一例について説明する。

まず、造影剤流入に関するパラメータを取得する場合について説明する。かかる場合、解析部１５２は、第１点が各輝度変化曲線で同一の正規化第１点にプロットされ、最大点が各輝度変化曲線で同一の正規化最大点にプロットされる正規化時間軸及び正規化輝度軸を設定することで、複数の輝度変化曲線それぞれから複数の正規化曲線を生成する。

具体的には、解析部１５２は、各輝度変化曲線の第１点と最大点との輝度幅及び時間軸を取得する。そして、解析部１５２は、取得した各輝度幅が一定値となるように、各輝度変化曲線の輝度軸のスケールを変更する。また、解析部１５２は、取得した各時間幅が一定値となるように、各輝度変化曲線の時間軸のスケールを変更する。そして、解析部１５２は、スケール変更後の輝度軸及び時間軸において、各輝度変化曲線の第１点を同一座標の正規化第１点に設定し、各輝度変化曲線の最大点を同一座標の正規化最大点に設定する。これにより、解析部１５２は、正規化時間軸及び正規化輝度軸を設定する。そして、解析部１５２は、各輝度変化曲線における第１点から最大点までの曲線を構成する各点を正規化時間軸及び正規化輝度軸において再プロットすることで、複数の輝度変化曲線それぞれから複数の正規化曲線を生成する。

例えば、解析部１５２は、図５に示す曲線Ｃ０、Ｃ１及びＣ２それぞれから、「Ｉ０ｍａｘ／２」、「Ｉ１ｍａｘ／２」及び「Ｉ２ｍａｘ／２」を取得する。また、例えば、解析部１５２は、図５に示す曲線Ｃ０、Ｃ１及びＣ２それぞれから、「ｔ０ｍａｘ−ｔ０ｓ＝ｔ０ｒ」、「ｔ１ｍａｘ−ｔ１ｓ＝ｔ１ｒ」及び「ｔ２ｍａｘ−ｔ２ｓ＝ｔ２ｒ」を取得する。そして、例えば、解析部１５２は、図６に示すように、「Ｉ０ｍａｘ／２、Ｉ１ｍａｘ／２、Ｉ２ｍａｘ／２」を「５０」とする。そして、例えば、解析部１５２は、図６に示すように、「ｔ０ｍａｘ−ｔ０ｓ＝ｔ０ｒ、ｔ１ｍａｘ−ｔ１ｓ＝ｔ１ｒ、ｔ２ｍａｘ−ｔ２ｓ＝ｔ２ｒ」を「１００」とする。これにより、解析部１５２は、正規化時間軸及び正規化輝度軸のスケールを決定する。

そして、解析部１５２は、例えば、曲線Ｃ０〜曲線Ｃ２の第１点が、正規化第１点「正規化時間：−１００、正規化輝度：５０」となり、曲線Ｃ０〜曲線Ｃ２の最大点が、正規化最大点「正規化時間：０、正規化輝度：１００」となるように、正規化時間軸及び正規化輝度軸の座標を決定する。これにより、解析部１５２は、正規化時間軸及び正規化輝度軸の設定を完了する。そして、解析部１５２は、曲線Ｃ０の第１点から最大点までの曲線を構成する各点を正規化時間軸及び正規化輝度軸において再プロットすることで、図６に示す正規化曲線ＮＣ０（ｉｎ）を生成する。同様に、解析部１５２は、曲線Ｃ１から、図６に示す正規化曲線ＮＣ１（ｉｎ）を生成する。同様に、解析部１５２は、曲線Ｃ２から、図６に示す正規化曲線ＮＣ２（ｉｎ）を生成する。

次に、造影剤流出に関するパラメータを取得する場合について説明する。かかる場合、解析部１５２は、最大点が各輝度変化曲線で同一の正規化最大点にプロットされ、第２点が各輝度変化曲線で同一の正規化第２点にプロットされる正規化時間軸及び正規化輝度軸を設定することで、複数の輝度変化曲線それぞれから複数の正規化曲線を生成する。

具体的には、解析部１５２は、各輝度変化曲線の最大点と第２点との輝度幅及び時間軸を取得する。そして、解析部１５２は、取得した各輝度幅が一定値となるように、各輝度変化曲線の輝度軸のスケールを変更する。また、解析部１５２は、取得した各時間幅が一定値となるように、各輝度変化曲線の時間軸のスケールを変更する。そして、解析部１５２は、スケール変更後の輝度軸及び時間軸において、各輝度変化曲線の最大点を同一座標の正規化最大点に設定し、各輝度変化曲線の第２点を同一座標の正規化第２点に設定する。これにより、解析部１５２は、正規化時間軸及び正規化輝度軸を設定する。そして、解析部１５２は、各輝度変化曲線における最大点から第２点までの曲線を構成する各点を正規化時間軸及び正規化輝度軸において再プロットすることで、複数の輝度変化曲線それぞれから複数の正規化曲線を生成する。

例えば、解析部１５２は、図５に示す曲線Ｃ０、Ｃ１及びＣ２それぞれから、「Ｉ０ｍａｘ／２」、「Ｉ１ｍａｘ／２」及び「Ｉ２ｍａｘ／２」を取得する。本実施形態では、第１割合及び第２割合が同一の割合であることから、輝度変化曲線ごとに、最大点と第２点との輝度幅は、最大点と第１点との輝度幅と同じ値となる。また、例えば、解析部１５２は、図５に示す曲線Ｃ０、Ｃ１及びＣ２それぞれから、「ｔ０ｅ−ｔ０ｍａｘ＝ｔ０ｐ」、「ｔ１ｅ−ｔ１ｍａｘ＝ｔ１ｐ」及び「ｔ２ｅ−ｔ２ｍａｘ＝ｔ２ｐ」を取得する。そして、例えば、解析部１５２は、図７に示すように、「Ｉ０ｍａｘ／２、Ｉ１ｍａｘ／２、Ｉ２ｍａｘ／２」を「５０」とする。そして、例えば、解析部１５２は、図７に示すように、「ｔ０ｅ−ｔ０ｍａｘ＝ｔ０ｐ、ｔ１ｅ−ｔ１ｍａｘ＝ｔ１ｐ、ｔ２ｅ−ｔ２ｍａｘ＝ｔ２ｐ」を「１００」とする。これにより、解析部１５２は、正規化時間軸及び正規化輝度軸のスケールを決定する。

そして、解析部１５２は、例えば、曲線Ｃ０〜曲線Ｃ２の最大点が、正規化最大点「正規化時間：０、正規化輝度：１００」となり、曲線Ｃ０〜曲線Ｃ２の第２点が、正規化第２点「正規化時間：１００、正規化輝度：５０」となるように、正規化時間軸及び正規化輝度軸の座標を決定する。これにより、解析部１５２は、正規化時間軸及び正規化輝度軸の設定を完了する。そして、解析部１５２は、曲線Ｃ０の最大点から第２点までの曲線を構成する各点を正規化時間軸及び正規化輝度軸において再プロットすることで、図７に示す正規化曲線ＮＣ０（ｏｕｔ）を生成する。同様に、解析部１５２は、曲線Ｃ１から、図７に示す正規化曲線ＮＣ１（ｏｕｔ）を生成する。同様に、解析部１５２は、曲線Ｃ７から、図６に示す正規化曲線ＮＣ２（ｏｕｔ）を生成する。

次に、造影剤流入及び造影剤流出に関するパラメータを取得する場合について説明する。かかる場合、解析部１５２は、第１点、最大点及び第２点それぞれが各輝度変化曲線で正規化第１点、正規化最大点及び正規化第２点それぞれにプロットされる正規化時間軸及び正規化輝度軸を設定することで、複数の輝度変化曲線それぞれから複数の正規化曲線を生成する。

具体的には、解析部１５２は、各輝度変化曲線の第１点と最大点との輝度幅（第１輝度幅）及び時間幅（第１時間軸）を取得する。また、解析部１５２は、各輝度変化曲線の最大点と第２点との輝度幅（第２輝度幅）及び時間幅（第２時間軸）を取得する。そして、解析部１５２は、各輝度変化曲線の第１輝度幅が一定値（ｄＩ１）となり、各輝度変化曲線の第２輝度幅が一定値（ｄＩ２）となるように、各輝度変化曲線の輝度軸のスケールを変更する。ここで、解析部１５２は、「ｄＩ１：ｄＩ２＝第１割合：第２割合」とする。また、解析部１５２は、各輝度変化曲線の第１時間幅が一定値（ｄＴ１）となり、各輝度変化曲線の第２時間幅が一定値（ｄＴ２）となるように、各輝度変化曲線の時間軸のスケールを変更する。ここで、解析部１５２は、「ｄＴ１：ｄＴ２＝第１割合：第２割合」とする。

そして、解析部１５２は、スケール変更後の輝度軸及び時間軸において、各輝度変化曲線の第１点を同一座標の正規化第１点に設定し、各輝度変化曲線の最大点を同一座標の正規化最大点に設定し、各輝度変化曲線の第２点を同一座標の正規化第２点に設定する。例えば、第１割合が「２０％」であり、第２割合が「３０％」である場合、正規化第１点の座標は、「正規化時間：−１００、正規化輝度：２０」に設定され、正規化最大点の座標は、「正規化時間：０、正規化輝度：１００」に設定され、正規化第２点の座標は、「正規化時間：１５０、正規化輝度：３０」に設定される。

これにより、解析部１５２は、正規化時間軸及び正規化輝度軸を設定する。そして、解析部１５２は、各輝度変化曲線における第１点から最大点を経由して第２点までに至る曲線を構成する各点を正規化時間軸及び正規化輝度軸において再プロットすることで、複数の輝度変化曲線それぞれから複数の正規化曲線を生成する。

本実施形態では第１割合及び第２割合が同じ「５０％」であるので、解析部１５２は、曲線Ｃ０から生成した正規化曲線ＮＣ０（ｉｎ）と正規化曲線ＮＣ０（ｏｕｔ）とを合わせることで、図８に示す正規化曲線ＮＣ０を生成する。同様に、解析部１５２は、曲線Ｃ１から生成した正規化曲線ＮＣ１（ｉｎ）と正規化曲線ＮＣ１（ｏｕｔ）とを合わせることで、図８に示す正規化曲線ＮＣ１を生成する。同様に、解析部１５２は、曲線Ｃ２から生成した正規化曲線ＮＣ２（ｉｎ）と正規化曲線ＮＣ２（ｏｕｔ）とを合わせることで、図８に示す正規化曲線ＮＣ２を生成する。

解析部１５２は、上述した正規化曲線から、正規化されたパラメータを取得する。例えば、解析部１５２は、正規化曲線において、正規化輝度が「８０」となる正規化時間や、正規化時間が「５０」となる正規化輝度を、正規化パラメータとして取得する。

そして、制御部１８は、パラメータ（正規化パラメータ）を画像又は文字のうちいずれか一方の形式でモニタ２に表示させる。パラメータの表示形態については、様々な形態を実行可能であるが、本実施形態では、パラメータの表示が画像の形式で行われる場合について説明する。具体的には、以下では、画像（画像形式）による表示形態の１つとして、正規化曲線から得られるパラメータを用いたパラメトリックイメージングが実行される場合について説明する。なお、パラメータの文字の形式での表示形態、並びに、パラメトリックイメージング以外のパラメータの画像の形式での表示形態については、後に詳述する。

画像による表示形態の１つとしてパラメトリックイメージングが設定されている場合、図１に示す変化画像生成部１５３は、制御部１８の指示により、以下の処理を行なう。すなわち、変化画像生成部１５３は、パラメータの値に応じて色調を変化させた変化画像データを生成する。そして、制御部１８は、画像による表示形態の１つとして、変化画像データをモニタ２に表示させる。本実施形態では、画像による表示形態の１つとして変化画像データの生成表示が設定されている。このことから、変化画像生成部１５３は、複数の正規化曲線それぞれで取得されたパラメータを用いて、変化画像データを生成する。以下、図９〜図１１を用いて、変化画像生成部１５３が生成する変化画像データについて説明する。図９〜図１１は、変化画像生成部を説明するための図である。

変化画像生成部１５３は、造影剤流入、又は、造影剤流出に関するパラメータを画像化する場合、正規化時間軸における正規化時間に応じて異なる色調が対応付けられた対応マップ（時間カラーマップ）を用いて、変化画像データを生成する。時間カラーマップは、例えば、内部記憶部１７に予め記憶されている。図９は、図７に示す正規化曲線ＮＣ０（ｏｕｔ）、ＮＣ１（ｏｕｔ）及びＮＣ２（ｏｕｔ）を用いて、造影剤流出に関するパラメータとして、正規化時間を画像化する場合の一例である。

例えば、制御部１８は、図９の上図に示すように、正規化曲線ＮＣ０（ｏｕｔ）、ＮＣ１（ｏｕｔ）及びＮＣ２（ｏｕｔ）をモニタ２に表示させる。また、制御部１８は、更に、任意の正規化輝度を操作者が設定可能なスライドバーＢ１を表示させる。スライドバーＢ１は、図９の上図に示すように、正規化時間軸に平行し、正規化輝度軸に直交した線となる。また、制御部１８は、図９の上図に示すように、正規化時間軸上に、正規化時間軸のスケールと同じスケールで時間カラーマップを表示させる。なお、時間カラーマップが表示される位置及びスケールは、任意に変更可能である。

そして、操作者は、例えば、図９の上図に示すように、正規化輝度「８０」の位置にスライドバーＢ１を移動する。解析部１５２は、正規化輝度「８０」となる正規化時間を、ＮＣ０（ｏｕｔ）、ＮＣ１（ｏｕｔ）及びＮＣ２（ｏｕｔ）それぞれから取得する。そして、解析部１５２は、ＮＣ０（ｏｕｔ）で取得した正規化時間を解析領域１００のパラメータとし、ＮＣ１（ｏｕｔ）で取得した正規化時間を解析領域１０１のパラメータとし、ＮＣ２（ｏｕｔ）で取得した正規化時間を解析領域１０２のパラメータとして変化画像生成部１５３に通知する。

変化画像生成部１５３は、図９の下図に示すように、ＮＣ０（ｏｕｔ）で取得された正規化時間に対応する色調を時間カラーマップから取得し、取得した色調により、超音波画像データの解析領域１００を着色する。また、変化画像生成部１５３は、図９の下図に示すように、ＮＣ１（ｏｕｔ）で取得された正規化時間に対応する色調を時間カラーマップから取得し、取得した色調により、超音波画像データの解析領域１０１を着色する。また、変化画像生成部１５３は、図９の下図に示すように、ＮＣ２（ｏｕｔ）で取得された正規化時間に対応する色調を時間カラーマップから取得し、取得した色調により、超音波画像データの解析領域１０２を着色する。なお、時間カラーマップから取得した色調により着色される超音波画像データは、例えば、解析領域１００〜１０２が設定された超音波画像データである。

そして、制御部１８は、図９の下図に示す変化画像データをモニタ２に表示させる。この変化画像データは、造影剤流出過程で、造影剤の存在量が最大量から「最大量の所定割合まで減少する流出時間を各解析領域で正規化して画像化したデータとなる。

なお、操作者がスライドバーＢ１を移動するに応じて、解析部１５２は、各解析領域のパラメータを更新して取得し、変化画像生成部１５３は、変化画像データを更新して生成する。なお、正規化輝度の設定は、操作者が数値を入力する場合等、任意の方法により行なっても良い。また、本実施形態は、例えば、自動的に、正規化輝度の値が変更されることで、変化画像データが動画として生成表示される場合であっても良い。

なお、造影剤流入に関するパラメータ（正規化時間）を画像化する場合は、例えば、図６に示す正規化曲線ＮＣ０（ｉｎ）、ＮＣ１（ｉｎ）及びＮＣ２（ｉｎ）を用いて、同様に行なわれる。かかる場合に生成表示される変化画像データは、造影剤流出過程で、造影剤の存在量が最大量の所定割合から最大量まで増加する流入時間を各解析領域で正規化して画像化したデータとなる。

また、変化画像生成部１５３は、造影剤流入、又は、造影剤流出に関するパラメータを画像化する場合、正規化輝度軸における正規化輝度に応じて異なる色調が対応付けられた対応マップ（輝度カラーマップ）を用いて、変化画像データを生成する。輝度カラーマップは、例えば、内部記憶部１７に予め記憶されている。図１０は、図７に示す正規化曲線ＮＣ０（ｏｕｔ）、ＮＣ１（ｏｕｔ）及びＮＣ２（ｏｕｔ）を用いて、造影剤流出に関するパラメータとして、正規化輝度を画像化する場合の一例である。

例えば、制御部１８は、図１０の上図に示すように、正規化曲線ＮＣ０（ｏｕｔ）、ＮＣ１（ｏｕｔ）及びＮＣ２（ｏｕｔ）をモニタ２に表示させる。また、制御部１８は、更に、任意の正規化時間を操作者が設定可能なスライドバーＢ２を表示させる。スライドバーＢ２は、図１０の上図に示すように、正規化輝度軸に平行し、正規化時間軸に直交した線となる。また、制御部１８は、図１０の上図に示すように、正規化輝度軸上に、正規化輝度軸のスケールと同じスケールで輝度カラーマップを表示させる。なお、輝度カラーマップが表示される位置及びスケールは、任意に変更可能である。

そして、操作者は、例えば、図１０の上図に示すように、正規化時間「６０」の位置にスライドバーＢ２を移動する。解析部１５２は、正規化時間「６０」となる正規化輝度を、ＮＣ０（ｏｕｔ）、ＮＣ１（ｏｕｔ）及びＮＣ２（ｏｕｔ）それぞれから取得する。そして、解析部１５２は、ＮＣ０（ｏｕｔ）で取得した正規化輝度を解析領域１００のパラメータとし、ＮＣ１（ｏｕｔ）で取得した正規化輝度を解析領域１０１のパラメータとし、ＮＣ２（ｏｕｔ）で取得した正規化輝度を解析領域１０２のパラメータとして変化画像生成部１５３に通知する。

変化画像生成部１５３は、図１０の下図に示すように、ＮＣ０（ｏｕｔ）で取得された正規化輝度に対応する色調を輝度カラーマップから取得し、取得した色調により、超音波画像データの解析領域１００を着色する。また、変化画像生成部１５３は、図１０の下図に示すように、ＮＣ１（ｏｕｔ）で取得された正規化輝度に対応する色調を輝度カラーマップから取得し、取得した色調により、超音波画像データの解析領域１０１を着色する。また、変化画像生成部１５３は、図１０の下図に示すように、ＮＣ２（ｏｕｔ）で取得された正規化輝度に対応する色調を輝度カラーマップから取得し、取得した色調により、超音波画像データの解析領域１０２を着色する。なお、輝度カラーマップから取得した色調により着色される超音波画像データは、例えば、解析領域１００〜１０２が設定された超音波画像データである。

そして、制御部１８は、図１０の下図に示す変化画像データをモニタ２に表示させる。この変化画像データは、造影剤流出過程を正規化した時間軸上の同一時点で、各解析領域から流出した造影剤の流出量を正規化して画像化したデータとなる。

なお、操作者がスライドバーＢ２を移動するに応じて、解析部１５２は、各解析領域のパラメータを更新して取得し、変化画像生成部１５３は、変化画像データを更新して生成する。なお、正規化時間の設定は、操作者が数値を入力する場合等、任意の方法により行なっても良い。また、本実施形態は、例えば、自動的に、正規化時間の値が変更されることで、変化画像データが動画として生成表示される場合であっても良い。

なお、造影剤流入に関するパラメータ（正規化輝度）を画像化する場合は、例えば、図６に示す正規化曲線ＮＣ０（ｉｎ）、ＮＣ１（ｉｎ）及びＮＣ２（ｉｎ）を用いて、同様に行なわれる。かかる場合に生成表示される変化画像データは、造影剤流入過程を正規化した時間軸上の同一時点で、各解析領域へ流入した造影剤の流入量を正規化して画像化したデータとなる。

また、変化画像生成部１５３は、造影剤流入及び造影剤流出に関するパラメータを画像化する場合、第１対応マップ（第１時間カラーマップ）と第２対応マップ（第２時間カラーマップ）とを混同した第３対応マップを用いて、変化画像データを生成する。ここで、第１時間カラーマップは、正規化時間軸において、正規化最大点の正規化最大時間以前の正規化時間に応じて第１色相で異なる色調が対応付けられたマップである。また、第２時間カラーマップは、正規化時間軸において、正規化最大時間以後の正規化時間に応じて第２色相で異なる色調が対応付けられたマップである。例えば、第１時間カラーマップは、青色系のカラーマップであり、第２時間カラーマップは、赤色系のカラーマップである。第１時間カラーマップ及び第２時間カラーマップは、例えば、内部記憶部１７に予め記憶されている。図１１は、図８に示す正規化曲線ＮＣ０、ＮＣ１及びＮＣ２を用いて、造影剤流入及び造影剤流出に関するパラメータとして、正規化時間を画像化する場合の一例である。

例えば、制御部１８は、図１１に示すように、正規化曲線ＮＣ０、ＮＣ１及びＮＣ２をモニタ２に表示させる。また、制御部１８は、更に、任意の正規化輝度を操作者が設定可能なスライドバーＢ３を表示させる。スライドバーＢ３は、図１１に示すように、正規化時間軸に平行し、正規化輝度軸に直交した線となる。また、制御部１８は、図１１に示すように、正規化時間軸上に、正規化時間軸のスケールと同じスケールで第１時間カラーマップ及び第２時間カラーマップを表示させる。図１１では、正規化最大時間は「０」であり、第１時間カラーマップは「−１００〜０」の正規化時間軸上にスケーリングされて表示され、第２時間カラーマップは、「０〜１００」の正規化時間軸上にスケーリングされて表示されている。なお、第１時間カラーマップ及び第２時間カラーマップが表示される位置及びスケールは、任意に変更可能である。

そして、操作者は、例えば、図１１に示すように、正規化輝度「６５」の位置にスライドバーＢ３を移動する。解析部１５２は、正規化輝度「６５」となる２つの正規化時間（マイナスの正規化時間及びプラスの正規化時間）を、ＮＣ０、ＮＣ１及びＮＣ２それぞれから取得する。そして、解析部１５２は、ＮＣ０で取得した２つ正規化時間を解析領域１００のパラメータとし、ＮＣ１で取得した２つ正規化時間を解析領域１０１のパラメータとし、ＮＣ２で取得した２つ正規化時間を解析領域１０２のパラメータとして変化画像生成部１５３に通知する。

変化画像生成部１５３は、図１１に示すように、ＮＣ０で取得されたマイナスの正規化輝度に対応する色調を第１時間カラーマップから取得し、ＮＣ０で取得されたプラスの正規化輝度に対応する色調を第２時間カラーマップから取得する。そして、変化画像生成部１５３は、図１１に示すように、取得した２つの色調を混合した色調により、超音波画像データの解析領域１００を着色する。

変化画像生成部１５３は、ＮＣ１で取得された２つの正規化輝度に対して、同様の色調取得処理を行ない、図１１に示すように、取得した２つの色調を混合した色調により、超音波画像データの解析領域１０１を着色する。また、変化画像生成部１５３は、ＮＣ２で取得された２つの正規化輝度に対して、同様の色調取得処理を行ない、図１１に示すように、取得した２つの色調を混合した色調により、超音波画像データの解析領域１０２を着色する。

そして、制御部１８は、図１１を用いて生成された変化画像データをモニタ２に表示させる。この変化画像データは、「造影剤の存在量が最大量から最大量の所定割合まで減少する流出時間」と「造影剤の存在量が最大量の所定割合から最大量まで増加する流入時間」とを各解析領域で正規化して、これらの正規化時間を同時に画像化したデータとなる。

なお、操作者がスライドバーＢ３を移動するに応じて、解析部１５２は、各解析領域のパラメータを更新して取得し、変化画像生成部１５３は、変化画像データを更新して生成する。なお、正規化時間の設定は、操作者が数値を入力する場合等、任意の方法により行なっても良い。また、本実施形態は、例えば、自動的に、正規化時間の値が変更されることで、変化画像データが動画として生成表示される場合であっても良い。また、本実施形態は、第１時間カラーマップ及び第２時間カラーマップを混合した２次元の時間カラーマップを用いる場合であっても良い。また、本実施形態は、２つの時間カラーマップを混合するのではなく、単に、正規化時間幅の値に応じた時間カラーマップを用いる場合であっても良い。

また、変化画像生成部１５３は、造影剤流入及び造影剤流出に関するパラメータを画像化する場合、以下の処理により変化画像データを生成しても良い。すなわち、変化画像生成部１５３は、第１輝度カラーマップと第２輝度カラーマップとを用いて、変化画像データを生成する。第１輝度カラーマップは、正規化輝度軸において、正規化最大点の正規化最大時間以前の正規化輝度に応じて第１色相で異なる色調が対応付けられた第１対応マップである。また、第２輝度カラーマップは、正規化輝度軸において、正規化最大時間以後の正規化輝度に応じて第２色相で異なる色調が対応付けられた第２対応マップである。

かかる場合、解析部１５２は、指定された２つの正規化時間「−Ｔ、＋Ｔ」に対応する２つの正規化輝度を各正規化曲線から取得する。そして、変化画像生成部１５３は、「−Ｔ」の正規化輝度に対応する色調を第１輝度カラーマップから取得し、「＋Ｔ」の正規化輝度に対応する色調を第２輝度カラーマップから取得し、取得した２つの色調を混合する。これにより、変化画像生成部１５３は、変化画像データを生成する。なお、上記の処理は、第１輝度カラーマップ及び第２輝度カラーマップを混合した２次元の輝度カラーマップを用いる場合であっても良い。また、本実施形態は、２つの輝度カラーマップを混合するのではなく、単に、正規化輝度幅の値に応じた輝度カラーマップを用いる場合であっても良い。

なお、図３や図４の設定が行なわれる場合は、異なる２つの時期の時系列データそれぞれから同一の解析領域における１つの輝度変化曲線が生成され、２つの正規化曲線が生成される。かかる場合は、変化画像生成部１５３は、同一の超音波画像データを２つ並列し、各超音波画像データの解析領域を、各正規化曲線から取得された正規化パラメータに応じた色調により着色する。

また、異なる３つの時期の時系列データそれぞれから同一の解析領域における１つの輝度変化曲線が生成され、３つの正規化曲線が生成される場合がある。かかる場合、変化画像生成部１５３は、同一の超音波画像データを３つ並列し、各超音波画像データの解析領域を、各正規化曲線から取得された正規化パラメータに応じた色調により着色する。

また、異なる２つの時期の時系列データそれぞれから２つの同一の解析領域それぞれで輝度変化曲線が生成され、２つの時期それぞれで２つの正規化曲線が生成される場合がある。かかる場合、変化画像生成部１５３は、同一の超音波画像データを２つ並列し、各超音波画像データの２つの解析領域それぞれを、各正規化曲線から取得された２つの正規化時間それぞれに応じた色調により着色する。

また、異なる２つの時期の時系列データそれぞれから１つ又は複数の同一の解析領域で輝度変化曲線が生成される場合、変化画像生成部１５３は各正規化曲線で得られた正規化パラメータの比に応じて色調を変化させることで、１つの変化画像データを生成しても良い。

なお、本実施形態は、変化画像データを参照した操作者が、正規化パラメータの値に応じて着色された解析領域を指定したことにともない、当該当該解析領域内、又は、当該解析領域近傍に、正規化パラメータの値を表示させても良い。また、本実施形態は、正規化パラメータの値に応じて解析領域を着色するとともに、当該解析領域内、又は、当該解析領域近傍に、正規化パラメータの値が文字で描画された超音波画像データを、変化画像データとして生成表示する場合であっても良い。また、本実施形態は、解析領域の着色を行なわずに、解析領域内、又は、解析領域近傍に、正規化パラメータの値が文字で描画された超音波画像データを、変化画像データとして生成表示する場合であっても良い。

次に、図１２を用いて、本実施形態に係る超音波診断装置が行なう処理の一例について説明する。図１２は、本実施形態に係る超音波診断装置が処理の一例を示すフローチャートである。なお、図１２は、解析領域の設定及び造影画像データ群の収集が完了し、輝度変化曲線の生成が開始された後の処理を示すフローチャートであり、パラメータの表示形態として変化画像データが設定されている場合の処理を示すフローチャートである。

図１２に例示するように、本実施形態に係る超音波診断装置の解析部１５２は、画像メモリ１６に、複数の輝度変化曲線が格納されたか否かを判定する（ステップＳ１０１）。ここで、複数の輝度変化曲線が格納されていない場合（ステップＳ１０１否定）、解析部１５２は、複数の輝度変化曲線が格納されるまで待機する。

一方、複数の輝度変化曲線が格納された場合（ステップＳ１０１肯定）、解析部１５２は、形状特徴の解析を行なって、複数の輝度変化曲線それぞれから正規化曲線を生成する（ステップＳ１０２）。そして、解析部１５２は、複数の正規化曲線それぞれから、正規化パラメータを取得する（ステップＳ１０３）。

そして、変化画像生成部１５３は、取得されたパラメータの値に応じた色調を対応マップから取得して、変化画像データを生成する（ステップＳ１０４）。そして、制御部１８の制御により、モニタ２は、変化画像データを表示し（ステップＳ１０５）、処理を終了する。

上述したように、本実施形態では、解析対象となる輝度変化曲線の形状特徴を解析して、正規化曲線を生成する。すなわち、本実施形態では、解析対象となる輝度変化曲線が如何なる条件（例えば、時系列データの撮影条件や解析領域の位置）で生成される場合であっても、客観的な同一の基準（最大輝度、第１割合及び第２割合）を用いて輝度変化曲線から正規化曲線を生成する。そして、本実施形態では、正規化曲線から、造影剤の流入量や流出量が正規化されたパラメータ、造影剤の流入時間や流出時間が正規化されたパラメータを取得する。

そして、本実施形態では、正規化されたパラメータにより、血流動態に関するパラメトリックイメージングを行なう。すなわち、本実施形態は、輝度変化曲線から得られる絶対値をパラメータとして用いた従来のパラメトリックイメージングではなく、正規化曲線から得られる相対値をパラメータとして用いたパラメトリックイメージングを行なう。従って、本実施形態では、造影剤の還流動態を客観的な基準により解析することができる。また、本実施形態では、造影剤の流入過程だけでなく、造影剤の流出過程も、正規化されたパラメータにより画像化することができる。

また、本実施形態では、正規化曲線を用いたパラメトリックイメージングを行なうことで、異なる解析領域の造影剤の還流動態を相対的に比較できる。例えば、本実施形態では、医師は、図９〜図１０等を用いて説明した変化画像データを参照して、基準となる組織（門脈や腎臓）の造影剤の還流動態と、腫瘍部位の造影剤の還流動態とを比較して、腫瘍部位の鑑別診断や腫瘍血管の異常の程度を判定することができる。

また、本実施形態では、正規化曲線を用いたパラメトリックイメージングを行なうことで、治療前後の同一解析領域における造影剤の還流動態を相対的に比較できる。例えば、本実施形態では、医師は、図３に例示する解析領域を設定することで生成された変化画像データを参照して、血管新生阻害剤による治療効果を判定することができる。

また、本実施形態では、正規化曲線を用いたパラメトリックイメージングを行なうことで、同一解析領域において異なる特性を有する複数の造影剤それぞれの還流動態を相対的に比較できる。例えば、本実施形態では、医師は、図４に例示する解析領域を設定することで生成された変化画像データを参照して、クッパー細胞に取り込まれやすい造影剤Ａの還流動態と、クッパー細胞に取り込まれにくい造影剤Ｂの還流動態とを比較して、腫瘍部位の鑑別診断や腫瘍血管の異常の程度を判定することができる。

（変形例）
上記の実施形態に係る超音波診断は、上述した処理以外にも、様々な変形例により行なうことが可能である。以下では、上記の実施形態の様々な変形例ついて説明する。なお、以下に説明する様々な変形例の処理は、上記の実施形態で説明した処理と、任意の形態で組み合わせることも可能である。

例えば、上記では、解析領域における造影剤の還流動態を輝度及び時間について正規化したパラメータを取得する場合について説明した。すなわち、上記では、輝度変化曲線の正規化を、時間軸及び輝度軸双方に対して行なう場合について説明した。しかし、解析部１５２は、解析領域における造影剤の還流動態を時間について正規化したパラメータを取得する場合であっても良い。すなわち、本実施形態は、輝度軸については正規化処理を行なわず、時間軸に対して正規化処理を行なって正規化時間軸を設定することで、正規化曲線を生成する場合であっても良い。かかる場合、解析部１５２は、輝度については実データのままとし、時間軸を正規化曲線にスケーリングすることで、輝度変化曲線から正規化曲線を生成する。そして、解析部１５２は、指定された正規化時間に対応する輝度（絶対輝度）を取得し、変化画像生成部１５３は、取得された輝度に応じて色調を変化させた変化画像データを生成する。

また、解析部１５２は、操作者が指定した場合、解析領域における造影剤の還流動態を輝度について正規化したパラメータを取得する場合であっても良い。すなわち、本実施形態は、時間軸については正規化処理を行なわず、輝度軸に対して正規化処理を行なって正規化輝度軸を設定することで、正規化曲線を生成する場合であっても良い。かかる場合、解析部１５２は、時間については実データのままとし、輝度軸を正規化曲線にスケーリングすることで、輝度変化曲線から正規化曲線を生成する。そして、解析部１５２は、指定された正規化輝度に対応する時間（絶対時間）を取得し、変化画像生成部１５３は、取得された時間に応じて色調を変化させた変化画像データを生成する。

また、上記実施形態の変形例として、解析部１５２が正規化曲線そのものをパラメータとして取得し、制御部１８が、画像での表示形態の１つとして、正規化曲線をモニタ２に表示させる場合であっても良い。正規化曲線は、造影剤の還流動態が正規化された曲線であることから、操作者は、正規化曲線そのものを参照しても、造影剤の還流動態を客観的な基準により解析することができる。そこで、例えば、解析部１５２は、複数の正規化曲線を生成した場合、これら複数の正規化曲線をパラメータとして制御部１８に出力する。そして、制御部１８は、複数の正規化曲線をモニタ２に表示させる。

この変形例では、図６〜図８それぞれに例示したグラフが、モニタ２に表示される。或いは、パラメータとして表示される正規化曲線は、上記の変形例で説明したように、一方の軸のみが正規化された曲線であっても良い。なお、上記実施形態の変形例として、解析部１５２は、正規化曲線をパラメータとして取得するとともに、正規化曲線からパラメータを取得する場合であっても良い。例えば、正規化曲線と、上記の実施形態で説明した変化画像データとは、同時にパラメータとして表示される場合であっても良い。

また、上記実施形態の変形例として、解析部１５２が、正規化曲線から取得した少なくとも１つの値をパラメータとして制御部１８に出力し、制御部１８が、上記の少なくとも１つの値をテーブルとして、又は、グラフとしてモニタ２に表示させても良い。この変形例では、解析部１５２は、従来、輝度変化曲線（近似曲線）から得ていたパラメータに対応するパラメータを正規化曲線から取得する。なお、この変形例で用いられる正規化曲線は、２軸が正規化された曲線でも、一方の軸のみが正規化された曲線であっても良い。

以下、従来、解析領域の輝度変化曲線から得られる代表的なパラメータについて説明する。従来のパラメータとしては、例えば、輝度の最大値（最大輝度）、輝度が最大値となるまでに要する時間（最大輝度時間）、平均通過時間（ＭＴＴ：Mean Transit Time）がある。ＭＴＴは、造影剤が流入して輝度が「最大輝度の５０％」となった時点から、最大輝度となった後に、造影剤が流出して輝度が「最大輝度の５０％」となる時点までの時間である。

また、従来のパラメータとしては、例えば、造影剤流入中に、輝度が「最大輝度の５０％」となった時点での輝度変化曲線の微分値、すなわち、傾き（Slope）がある。また、従来のパラメータとしては、例えば、輝度変化曲線の輝度を『造影剤流入時点から最大輝度時間までの積分期間で積分した面積値「Area Wash in」』、輝度変化曲線の輝度を『最大輝度時間から造影剤流出時点までの積分期間で積分した面積値「Area Wash out」』、輝度変化曲線の輝度を『造影剤流入時点から造影剤流出時点までの積分期間で積分した面積値「Area Under Curve」』がある。「Area Wash in」は、造影剤の流入期間内に解析領域に存在した造影剤の総量を示す値となる。「Area Wash out」は、造影剤の流出期間内に解析領域に存在した造影剤の総量を示す値となる。「Area Under Curve」は、造影剤の流入時点から流出時点までに解析領域に存在した造影剤の総量を示す値となる。

以下、解析部１５２が、図１１に示す３つの正規化曲線を用いて解析領域１００、２００及び３００ぞれぞれでの「造影剤の還流動態を客観的に評価できる代表的な正規化パラメータ」を取得する場合について説明する。図１３及び図１４は、変形例を説明するための図である。例えば、従来のＭＴＴに対応する正規化パラメータを得るために、解析部１５２は、正規化輝度の値が立ち上がって正規化最大輝度「１００」の６５％となった時点から、正規化輝度の値が立ち下がって正規化最大輝度「１００」の６５％となった時点までの時間（正規化時間）を取得する。この時間を、解析部１５２は、「６５％」での正規化平均通過時間（ｎＭＴＴ＠６５％）として取得する。解析部１５２は、解析領域１００、２００及び３００ぞれぞれの「ｎＭＴＴ＠６５％」を取得する。なお、正規化平均通過時間の算出に用いる割合は、６５％以外、任意の値に変更可能である。

また、例えば、従来の「最大輝度の５０％」となった時点での「Slope」に対応する正規化パラメータを得るために、解析部１５２は、造影剤流入中に正規化最大輝度「１００」の６５％となった時間での正規化曲線の傾きを、「ｎＳｌｏｐｅ＠６５％」として取得する。解析部１５２は、解析領域１００、２００及び３００ぞれぞれの「ｎＳｌｏｐｅ＠６５％」を取得する。

また、例えば、従来の「Area Under Curve」に対応する正規化パラメータを得るために、解析部１５２は、正規化曲線の正規化輝度を正規化時間「−１００〜１００」まで積分した面積値を、「ｎＡｒｅａ」として取得する。解析部１５２は、解析領域１００、２００及び３００ぞれぞれの「ｎＡｒｅａ」を取得する。なお、解析部１５２は、正規化曲線の正規化輝度を正規化時間「−１００〜０」まで積分した面積値した値を、「Area Wash in」に対応する正規化パラメータとして取得しても良い。また、解析部１５２は、正規化曲線の正規化輝度を正規化時間「０〜１００」まで積分した面積値した値を、「Area Wash out」に対応する正規化パラメータとして取得しても良い。

そして、例えば、制御部１８は、図１３に示すように、解析領域１００、２００及び３００ぞれぞれの「ｎＭＴＴ＠６５％」と、解析領域１００、２００及び３００ぞれぞれの「ｎＳｌｏｐｅ＠６５％」と、解析領域１００、２００及び３００ぞれぞれの「ｎＡｒｅａ」とを、テーブルに変換し、モニタ２に表示させる。テーブル形式の表示形態は、文字形式での表示形態の一例である。或いは、例えば、制御部１８は、図１４に示すように、解析領域１００、２００及び３００ぞれぞれの「ｎＭＴＴ＠６５％」を、棒グラフに変換し、モニタ２に表示させる。また、図示していないが、制御部１８は、解析領域１００、２００及び３００ぞれぞれの他の正規化パラメータも、棒グラフに変換し、モニタ２に表示させる。棒グラフ形式の表示形態は、画像形式での表示形態の一例である。上記の変形例を行なうことでも、造影剤の還流動態を客観的な基準により解析することができる。

なお、上記の変形例では、解析部１５２が、正規化曲線の時間軸上で１時点での傾きを正規化パラメータとして取得する場合について説明した。しかし、解析部１５２は、正規化曲線の時間軸上で複数時点での傾きを正規化パラメータとして取得する場合であっても良い。すなわち、上記の変形例は、解析部１５２が、正規化曲線の各正規化時間での微分値を、正規化曲線として算出しても良い。かかる場合、制御部１８は、各正規化時間での微分値をテーブルとして表示させる。或いは、制御部１８は、各正規化時間での微分値をプロットしたグラフを生成し、このグラフを表示させる。

また、本実施形態の変形例として、輝度変化曲線が１つ生成される場合であっても良い。かかる場合、解析部１５２は、１つの輝度変化曲線から上述した正規化曲線を生成する。そして、上記の実施形態や変形例で説明したように、制御部１８は、様々な形式でパラメータを表示させる。例えば、制御部１８は、１つの正規化曲線から生成された変化画像データをモニタ２に表示させる。これによっても、造影剤の還流動態を客観的な基準により解析することができる。また、上記の画像処理方法は、造影剤の還流動態を客観的な基準により解析可能とすることから、異なる被検体それぞれで解析領域が設定される場合でも適用可能である。

例えば、解析部１５２は、患者Ａの肝臓の腫瘍部位に設定された解析領域の輝度変化曲線から正規化曲線Ａを生成する。また、例えば、患者Ｂの肝臓の腫瘍部位に設定された解析領域の輝度変化曲線から正規化曲線Ｂを生成する。なお、両者の腫瘍部位は、解剖学的に略同じ部位であることが好適である。そして、例えば、変化画像生成部１５３は、正規化曲線Ａの変化画像データＡを生成し、正規化曲線Ｂの変化画像データＢを生成する。或いは、例えば、解析部１５２は、正規化曲線ＡのｎＭＴＴ（Ａ）を算出し、正規化曲線ＢのｎＭＴＴ（Ｂ）を算出する。仮に、患者Ａと患者Ｂとで肝臓がんの進行度が異なる場合、正規化パラメータの値が異なる可能性が高い。すなわち、患者Ａと患者Ｂとで、肝臓がんの進行度が異なる場合、変化画像データＡと変化画像データＢとでは、色調のパターンが異なり、ｎＭＴＴ（Ａ）とｎＭＴＴ（Ｂ）とでは、値が異なる。このようなことから、医師は、例えば、変化画像データＡと変化画像データＢとを比較することで、肝臓がんの進行度の違いを判断することができる。

また、上記の方法により、例えば、肝臓がんの進行度の異なる複数の患者それぞれの正規化パラメータを収集してデータベース化することが可能となる。かかる場合、医師は、新たな肝臓がんの患者Ｃの正規化パラメータが得られた際には、データベースを参照して、患者Ｃの進行度を鑑別することができる。

また、上記では、造影期間中の時系列データが収集後に生成された輝度変化曲線を用いる場合について説明した。しかし、上記の実施形態の変形例として、造影期間中の時系列データの収集中に、輝度変化曲線をリアルタイムで生成する場合であっても良い。すなわち、本実施形態は、輝度変化曲線にて最大点が取得された時点から、少なくとも、造影剤流入に関する正規化パラメータの画像化等が、リアルタイムで行なわれる場合であっても良い。

なお、本実施形態及び変形例で説明した画像処理方法は、超音波診断装置とは独立に設置された画像処理装置により行なわれる場合であっても良い。かかる画像処理装置は、造影剤が投与された被検体Ｐを超音波走査して収集された時系列データを取得することで、本実施形態で説明した画像処理方法を行なうことができる。或いは、かかる画像処理装置は、輝度変化曲線を取得することで、本実施形態で説明した画像処理方法を行なう場合であっても良い。

また、図示した各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部又は一部を、各種の負荷や使用状況等に応じて、任意の単位で機能的又は物理的に分散・統合して構成することができる。更に、各装置にて行なわれる各処理機能は、その全部又は任意の一部が、ＣＰＵ及び当該ＣＰＵにて解析実行されるプログラムにて実現され、或いは、ワイヤードロジックによるハードウェアとして実現され得る。

また、本実施形態及び変形例で説明した画像処理方法は、あらかじめ用意された画像処理プログラムをパーソナルコンピュータやワークステーション等のコンピュータで実行することによって実現することができる。この画像処理プログラムは、インターネット等のネットワークを介して配布することができる。また、この制御プログラムは、ハードディスク、フレキシブルディスク（ＦＤ）、ＣＤ−ＲＯＭ、ＭＯ、ＤＶＤ、ＵＳＢメモリ及びＳＤカードメモリ等のＦｌａｓｈメモリ等のコンピュータで読み取り可能な非一時的な記録媒体に記録され、コンピュータによって非一時的な記録媒体から読み出されることによって実行することもできる。

以上、説明したとおり、本実施形態及び変形例によれば、造影剤の還流動態を客観的な基準により解析することができる。

本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。

１５ 画像処理部
１５１ 輝度変化情報生成部
１５２ 解析部
１５３ 変化画像生成部
１８ 制御部

Claims (14)

1. 造影剤が投与された被検体を超音波走査して収集された時系列データから、超音波走査領域内に設定された解析領域における輝度の時間変化を示す輝度変化情報を生成する輝度変化情報生成部と、
   前記輝度変化情報に基づいて、前記解析領域における造影剤の還流動態を時間について正規化したパラメータを取得する解析部と、
   前記パラメータを画像又は文字のうち少なくともいずれか一方の形式で表示部に表示させる制御部と、
   を備える、超音波診断装置。
2. 前記解析部は、前記解析領域における造影剤の還流動態を、輝度、又は、輝度及び時間について正規化したパラメータを取得する、請求項１に記載の超音波診断装置。
3. 前記パラメータの値に応じて色調を変化させた変化画像データを生成する変化画像生成部、
   を更に備え、
   前記制御部は、前記画像による表示形態の１つとして、前記変化画像データを表示部に表示させる、請求項２に記載の超音波診断装置。
4. 前記輝度変化情報生成部は、前記輝度変化情報として、前記解析領域における輝度の時間変化を示す曲線である輝度変化曲線を生成し、
   前記解析部は、時間軸、又は、輝度軸及び時間軸の正規化を行なって前記輝度変化曲線から正規化曲線を生成し、当該正規化曲線を前記パラメータとして取得する処理、及び、当該正規化曲線から前記パラメータを取得する処理の少なくとも１つの処理を行なう、請求項３に記載の超音波診断装置。
5. 前記解析部は、前記輝度変化曲線において、輝度が最大値となる最大点と、前記最大点となる前に前記最大値に第１割合を乗算した第１乗算値となる第１点と、前記最大点となった後に前記最大値に第２割合を乗算した第２乗算値となる第２点とから選択された少なくとも２点を用いて、前記正規化曲線を生成する、請求項４に記載の超音波診断装置。
6. 前記輝度変化情報生成部は、超音波走査領域内に設定された複数の解析領域それぞれにおける複数の輝度変化曲線、又は、異なる複数の時期に同一の超音波走査領域で行なわれた超音波走査により収集された複数の時系列データそれぞれから該領域内に設定された少なくとも１つの同一の解析領域それぞれにおける複数の輝度変化曲線を生成し、
   前記解析部は、前記複数の輝度変化曲線それぞれから前記正規化曲線を生成し、生成した複数の正規化曲線それぞれを前記パラメータとして取得する処理、及び、当該複数の正規化曲線それぞれから前記パラメータを取得する処理の少なくとも１つの処理を行ない、
   前記変化画像生成部は、前記画像による表示形態の１つとして設定されている場合、前記複数の正規化曲線それぞれで取得された前記パラメータを用いて、前記変化画像データを生成する、請求項５に記載の超音波診断装置。
7. 前記解析部は、
   造影剤流入に関するパラメータを取得する場合、前記第１点が各輝度変化曲線で同一の正規化第１点にプロットされ、前記最大点が各輝度変化曲線で同一の正規化最大点にプロットされる正規化時間軸及び正規化輝度軸を設定することで、前記複数の輝度変化曲線それぞれから複数の正規化曲線を生成し、
   造影剤流出に関するパラメータを取得する場合、前記最大点が各輝度変化曲線で同一の正規化最大点にプロットされ、前記第２点が各輝度変化曲線で同一の正規化第２点にプロットされる正規化時間軸及び正規化輝度軸を設定することで、前記複数の輝度変化曲線それぞれから複数の正規化曲線を生成し、
   造影剤流入及び造影剤流出に関するパラメータを取得する場合、前記第１点、前記最大点及び前記第２点それぞれが各輝度変化曲線で前記正規化第１点、前記正規化最大点及び前記正規化第２点それぞれにプロットされる正規化時間軸及び正規化輝度軸を設定することで、前記複数の輝度変化曲線それぞれから複数の正規化曲線を生成する、請求項６に記載の超音波診断装置。
8. 前記変化画像生成部は、造影剤流入、又は、造影剤流出に関するパラメータを画像化する場合、前記正規化時間軸における正規化時間に応じて異なる色調が対応付けられた対応マップを用いて、前記変化画像データを生成する、請求項７に記載の超音波診断装置。
9. 前記変化画像生成部は、造影剤流入、又は、造影剤流出に関するパラメータを画像化する場合、前記正規化輝度軸における正規化輝度に応じて異なる色調が対応付けられた対応マップを用いて、前記変化画像データを生成する、請求項７に記載の超音波診断装置。
10. 前記変化画像生成部は、造影剤流入及び造影剤流出に関するパラメータを画像化する場合、前記正規化時間軸において、前記正規化最大点の正規化最大時間以前の正規化時間に応じて第１色相で異なる色調が対応付けられた第１対応マップと、前記正規化最大時間以後の正規化時間に応じて第２色相で異なる色調が対応付けられた第２対応マップとを用いて、前記変化画像データを生成する、請求項７に記載の超音波診断装置。
11. 前記解析部は、前記正規化曲線から取得した少なくとも１つの値を前記パラメータとして前記制御部に出力し、
    前記制御部は、前記少なくとも１つの値をテーブルとして、又は、グラフとして前記表示部に表示させる、請求項４〜６のいずれか１つに記載の超音波診断装置。
12. 前記解析部は、前記正規化曲線の時間軸上での傾きを前記パラメータとして取得する、請求項１１に記載の超音波診断装置。
13. 造影剤が投与された被検体を超音波走査して収集された時系列データから、超音波走査領域内に設定された解析領域における輝度の時間変化を示す輝度変化情報を生成する輝度変化情報生成部と、
    前記輝度変化情報に基づいて、前記解析領域における造影剤の還流動態を時間について正規化したパラメータを取得する解析部と、
    前記パラメータを画像又は文字のうち少なくともいずれか一方の形式で表示部に表示させる制御部と、
    を備える、画像処理装置。
14. 輝度変化情報生成部が、造影剤が投与された被検体を超音波走査して収集された時系列データから、超音波走査領域内に設定された解析領域における輝度の時間変化を示す輝度変化情報を生成し、
    解析部が、前記輝度変化情報に基づいて、前記解析領域における造影剤の還流動態を時間について正規化したパラメータを取得し、
    制御部が、前記パラメータを画像又は文字のうち少なくともいずれか一方の形式で表示部に表示させる、
    ことを含む、画像処理方法。