(57)【要約】【課題】三次元画像に関する視点の位置を基準として速度情報を表示する。【解決手段】３Ｄデータ生成部４０は、データ空間内の各座標ごとに、３つの３Ｄプローブ１１，１２，１３から得られるエコーデータを合成処理して三次元の速度ベクトルを得る。レイ設定部５２は、視点からデータ空間内を通る複数のレイを設定する。ボクセル値算出部５０は、設定された各レイに沿って並ぶ複数のボクセルに関して、各ボクセルごとに速度値と輝度値を算出する。各ボクセルの速度値として、そのボクセルに対応した速度ベクトルのレイ方向成分が算出される。各レイに沿って並ぶ複数のボクセルの各々についての速度値と輝度値が算出されると、レンダリング処理部６０において、速度値と輝度値を利用して、各レイについてのレンダリング演算が実行される。【選択図】図１
2012045253.tif 000002
【特許請求の範囲】
【請求項１】
三次元空間に対して互いに異なる複数の原点方向から超音波を送受する送受波部と、
前記送受波部を制御することにより、前記各原点方向ごとに超音波ビームを立体的に走査して各原点方向ごとに三次元的にエコーデータを収集する送受信処理部と、
前記三次元空間に対応したデータ空間内の各座標ごとに、前記複数の原点方向から得られるその座標に対応した複数のエコーデータに基づいて、その座標における三次元の速度ベクトル情報を得る三次元情報生成部と、
前記データ空間に対して設定された視点側からデータ空間内を通る複数のレイを設定するレイ設定部と、
前記各レイごとに、そのレイに沿って並ぶ複数のボクセルに関して、各ボクセルに対応した座標の前記速度ベクトル情報に基づいて、各ボクセルごとにそのレイの方向の速度値を算出するボクセル値算出部と、
前記各レイごとに、そのレイに沿って並ぶ複数のボクセルに関する複数の速度値に基づいてレンダリング演算を実行するレンダリング処理部と、
前記レンダリング演算により前記複数のレイから得られる複数の画素値に基づいて前記三次元空間に関する表示画像を形成する画像形成部と、
を有する、
ことを特徴とする超音波診断装置。
【請求項２】
請求項１に記載の超音波診断装置において、
前記レンダリング処理部は、前記各レイごとに、前記速度値が正の値となる複数のボクセルを対象とした正のレンダリング演算と、前記速度値が負の値となる複数のボクセルを対象とした負のレンダリング演算と、を実行し、
前記画像形成部は、前記各レイごとに、前記正のレンダリング演算により得られる結果値と前記負のレンダリング演算により得られる結果値に基づいて画素値を決定する、
ことを特徴とする超音波診断装置。
【請求項３】
請求項２に記載の超音波診断装置において、
前記各ボクセルに対応した座標のエコーデータに基づいて各ボクセルごとに輝度値が決定され、
前記レンダリング処理部は、前記各レイごとに、前記正のレンダリング演算と前記負のレンダリング演算に加えて、そのレイに沿って並ぶ複数のボクセルに関する複数の輝度値に基づいた輝度のレンダリング演算を実行し、
前記画像形成部は、前記各レイごとに、前記正のレンダリング演算と前記負のレンダリング演算により得られる色相と、前記輝度のレンダリング演算により得られる輝度と、に基づいて画素値を決定する、
ことを特徴とする超音波診断装置。
【請求項４】
請求項３に記載の超音波診断装置において、
前記レンダリング処理部は、前記各レイごとに、レンダリング演算の進行と共に、そのレイに沿って並ぶ複数のボクセルに関する複数の速度値と複数の輝度値を複合的に累積処理し、その累積処理の結果が予め設定された閾値に達した時点をレンダリング演算の終了時点とする、
ことを特徴とする超音波診断装置。
【請求項５】
請求項１から４のいずれか１項に記載の超音波診断装置において、
前記送受波部は、三次元空間に対して互いに異なる複数の原点方向から超音波を送受する複数の三次元プローブを備える、
ことを特徴とする超音波診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、三次元空間内から得られる速度情報に基づいて表示画像を形成する超音波診断装置に関する。
【背景技術】
【０００２】
三次元空間内からエコー情報を収集して三次元の超音波画像を形成する装置が知られている。例えば、特許文献１には、三次元の生体組織画像と血流画像とを同時に表示する超音波診断装置が示されている。この装置により、例えば、生体組織部分を白で表示させて血流部分を赤く表示させ、さらに、生体組織部分と血流部分の前後関係をより明確に表現した画像を形成することができる。
【０００３】
また、例えば、正方向の血流に対して暖色系の色相を割り当て、負方向の血流に対して寒色系の色相を割り当てることにより、血流の方向を表現したカラードプラ画像が知られており、特許文献２には、カラードプラの情報を三次元画像に反映させる技術が示されている。
【０００４】
上述した基礎的な技術に支えられ、近年においては、三次元画像を形成する際の視点を超音波プローブに固定せず、視点の位置を任意に変更できる装置も登場している。この装置によれば、例えば医師等のユーザが診断状況等に応じた所望の位置に視点を設定し、診断対象の三次元画像を観察することが可能になる。
【０００５】
ところが、診断対象が血流でありその方向を色で表現する場合において、例えば、超音波ビームの方向を基準として、固定的に、正方向の血流を暖色系の色相として負方向の血流を寒色系の色相とすると、視点の位置が超音波プローブ側にある場合には、視点に向かってくる血流が暖色系となり視点から遠ざかる血流が寒色系となるため、カラードプラと同様な色の表現が維持されるものの、視点の位置が血流を挟んで超音波プローブの反対側に設定されてしまうと、視点に向かってくる血流が寒色系となり視点から遠ざかる血流が暖色系となるため、カラードプラにおける色の表現と逆になってしまう。
【０００６】
つまり、超音波ビームの方向を基準として速度情報を表示してしまうと、視点の位置の変更に対して十分に対応しきれない。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特許第３９７７７７９号公報
【特許文献２】特開２００５−２７８９８８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
上述した背景技術に鑑み、本願の発明者は、視点の位置を任意に変更して三次元画像を形成する技術について研究開発を重ねてきた。
【０００９】
本発明は、その研究開発の過程において成されたものであり、その目的は、三次元画像に関する視点の位置を基準として速度情報を表示する技術を提供することにある。
【課題を解決するための手段】
【００１０】
上記目的にかなう好適な超音波診断装置は、三次元空間に対して互いに異なる複数の原点方向から超音波を送受する送受波部と、前記送受波部を制御することにより、前記各原点方向ごとに超音波ビームを立体的に走査して各原点方向ごとに三次元的にエコーデータを収集する送受信処理部と、前記三次元空間に対応したデータ空間内の各座標ごとに、前記複数の原点方向から得られるその座標に対応した複数のエコーデータに基づいて、その座標における三次元の速度ベクトル情報を得る三次元情報生成部と、前記データ空間に対して設定された視点側からデータ空間内を通る複数のレイを設定するレイ設定部と、前記各レイごとに、そのレイに沿って並ぶ複数のボクセルに関して、各ボクセルに対応した座標の前記速度ベクトル情報に基づいて、各ボクセルごとにそのレイの方向の速度値を算出するボクセル値算出部と、前記各レイごとに、そのレイに沿って並ぶ複数のボクセルに関する複数の速度値に基づいてレンダリング演算を実行するレンダリング処理部と、前記レンダリング演算により前記複数のレイから得られる複数の画素値に基づいて前記三次元空間に関する表示画像を形成する画像形成部と、を有することを特徴とする。
【００１１】
上記構成によれば、三次元の速度ベクトル情報に基づいて各ボクセルごとにレイの方向の速度値が算出される。そのため、例えば、視点の位置の変更に伴ってレイの方向が変化した場合においても、変化後のレイの方向に応じた速度値が算出され、その速度値に基づいて表示画像が形成される。このように、例えば任意の位置に設定される視点を基準として速度値を得て表示画像を形成することが可能になる。
【００１２】
望ましい具体例において、前記レンダリング処理部は、前記各レイごとに、前記速度値が正の値となる複数のボクセルを対象とした正のレンダリング演算と、前記速度値が負の値となる複数のボクセルを対象とした負のレンダリング演算と、を実行し、前記画像形成部は、前記各レイごとに、前記正のレンダリング演算により得られる結果値と前記負のレンダリング演算により得られる結果値に基づいて画素値を決定する、ことを特徴とする。
【００１３】
望ましい具体例において、前記各ボクセルに対応した座標のエコーデータに基づいて各ボクセルごとに輝度値が決定され、前記レンダリング処理部は、前記各レイごとに、前記正のレンダリング演算と前記負のレンダリング演算に加えて、そのレイに沿って並ぶ複数のボクセルに関する複数の輝度値に基づいた輝度のレンダリング演算を実行し、前記画像形成部は、前記各レイごとに、前記正のレンダリング演算と前記負のレンダリング演算により得られる色相と、前記輝度のレンダリング演算により得られる輝度と、に基づいて画素値を決定する、ことを特徴とする。
【００１４】
望ましい具体例において、前記レンダリング処理部は、前記各レイごとに、レンダリング演算の進行と共に、そのレイに沿って並ぶ複数のボクセルに関する複数の速度値と複数の輝度値を複合的に累積処理し、その累積処理の結果が予め設定された閾値に達した時点をレンダリング演算の終了時点とする、ことを特徴とする。
【００１５】
望ましい具体例において、前記送受波部は、三次元空間に対して互いに異なる複数の原点方向から超音波を送受する複数の三次元プローブを備える、ことを特徴とする。
【発明の効果】
【００１６】
本発明により、三次元画像に関する視点の位置を基準として速度情報を表示することが可能になる。
【図面の簡単な説明】
【００１７】
【図１】本発明の実施において好適な超音波診断装置の全体構成を示す図である。
【図２】エコーデータの合成処理を説明するための図である。
【図３】ボリュームレンダリングの原理を説明するための図である。
【図４】レンダリング処理部における演算を説明するための図である。
【図５】レンダリング演算の結果に基づいた３原色の信号値を示す図である。
【発明を実施するための形態】
【００１８】
以下に本発明の好適な実施形態を説明する。
【００１９】
図１は、本発明の実施において好適な超音波診断装置の全体構成を示す図である。３Ｄプローブ１１は、血流などの対象物を含む三次元空間に対して超音波を送受波する超音波プローブである。３Ｄプローブ１１は、超音波を送受する複数の振動素子を備えており、複数の振動素子がビームフォーマ２１によって送信制御されて送信ビームが形成される。また、複数の振動素子が三次元空間内から得られる超音波を受波し、これにより得られた信号がビームフォーマ２１へ出力され、ビームフォーマ２１が受信ビームを形成する。
【００２０】
３Ｄプローブ１１は、超音波ビーム（送信ビームと受信ビーム）を三次元空間内において走査して立体的にエコーデータを収集する三次元プローブである。例えば、一次元的に配列された複数の振動素子（１Ｄアレイ振動子）によって電子的に形成される走査面を機械的に動かすことにより超音波ビームが立体的に走査される。また、二次元的に配列された複数の振動素子（２Ｄアレイ振動子）を電子的に制御して超音波ビームを立体的に走査してもよい。３Ｄプローブ１１を利用して三次元空間内から収集されたエコーデータは３Ｄメモリ３１に記憶される。
【００２１】
図１の超音波診断装置は、３Ｄプローブ１１に加えて、３Ｄプローブ１２，１３を備えている。３Ｄプローブ１２，１３も三次元プローブであり、例えば３Ｄプローブ１１と同じ構成のプローブが利用される。但し、３つの３Ｄプローブ１１，１２，１３は、互いに異なる方向から三次元空間に対して超音波を送受するように配置される。
【００２２】
３Ｄプローブ１２は、ビームフォーマ２２により制御される。これにより超音波ビームが三次元空間内において立体的に走査されてエコーデータが収集される。３Ｄプローブ１２を利用して三次元空間内から収集されたエコーデータは３Ｄメモリ３２に記憶される。
【００２３】
同様に、３Ｄプローブ１３は、ビームフォーマ２３により制御される。これにより、超音波ビームが三次元空間内において立体的に走査されてエコーデータが収集される。３Ｄプローブ１３を利用して三次元空間内から収集されたエコーデータは３Ｄメモリ３３に記憶される。
【００２４】
３つの３Ｄプローブ１１，１２，１３により三次元空間内からエコーデータが収集されると、３Ｄデータ生成部４０は、３Ｄメモリ３１，３２，３３に記憶されたそれらのエコーデータを合成処理することにより、三次元空間に対応したデータ空間内において、各座標ごとに輝度や速度ベクトルを算出する。
【００２５】
図２は、エコーデータの合成処理を説明するための図であり、図２には、図１を利用して説明した３Ｄプローブ１１，１２，１３が図示されている。図２に示すように、３つの３Ｄプローブ１１，１２，１３は、三次元空間１００に対して互いに異なる方向から超音波を送受するように配置されている。例えば、３つの３Ｄプローブ１１，１２，１３の互いの相対的な配置関係が既知とされ、３Ｄプローブ１２を基準として三次元空間１００の座標系が設定される。
【００２６】
３Ｄデータ生成部４０（図１）は、三次元空間１００に対応したデータ空間内の各座標ごとに、３つの３Ｄプローブ１１，１２，１３から得られるエコーデータを合成処理して三次元の速度ベクトルを得る。そこで、図２に示す三次元空間１００内の座標Ｃにおける合成処理について説明する。
【００２７】
まず、３Ｄプローブ１１から座標Ｃに向けられた超音波ビームＢ１に沿って得られるエコーデータに含まれるドプラ情報に基づいて、超音波ビームＢ１と同じ方向の速度ベクトル成分Ｖ１が決定される。３Ｄプローブ１１についてのエコーデータは、３Ｄメモリ３１（図１）から読み出される。また、３Ｄプローブ１２から座標Ｃに向けられた超音波ビームＢ２に沿って得られるエコーデータに含まれるドプラ情報に基づいて、超音波ビームＢ２と同じ方向の速度ベクトル成分Ｖ２が決定される。３Ｄプローブ１２についてのエコーデータは、３Ｄメモリ３２（図１）から読み出される。さらに、３Ｄプローブ１３から座標Ｃに向けられた超音波ビームＢ３に沿って得られるエコーデータに含まれるドプラ情報に基づいて、超音波ビームＢ３と同じ方向の速度ベクトル成分Ｖ３が決定される。３Ｄプローブ１３についてのエコーデータは、３Ｄメモリ３３（図１）から読み出される。そして、速度ベクトル成分Ｖ１と速度ベクトル成分Ｖ２と速度ベクトル成分Ｖ３をベクトル的に加算することにより、三次元空間１００内の座標Ｃにおける三次元の速度ベクトルＶが得られる。
【００２８】
また、３Ｄデータ生成部４０は、三次元空間１００に対応したデータ空間内の各座標ごとに、３つの３Ｄプローブ１１，１２，１３から得られるエコーデータに基づいて輝度（エコーデータの大きさ）を得る。なお、輝度はスカラー値であるため、例えば３つの３Ｄプローブ１１，１２，１３のうちのいずれか１つのプローブから得られるエコーデータのみにより座標Ｃにおける輝度が決定されてもよいし、複数のプローブから得られるエコーデータの平均値などから座標Ｃにおける輝度が決定されてもよい。
【００２９】
こうして、３Ｄデータ生成部４０は、三次元空間１００内（データ空間内）の全座標について、各座標ごとに三次元の速度ベクトルと輝度を決定する。全座標についての速度ベクトルと輝度は、例えば座標値をアドレスとして３Ｄデータ生成部４０が備えるメモリに記憶される。
【００３０】
図１の超音波診断装置は、３Ｄデータ生成部４０において得られた三次元空間１００内（データ空間内）の各座標の速度ベクトルと輝度に基づいて、ボリュームレンダリングの手法を利用して表示画像を形成する。そこで、そのボリュームレンダリングの原理について説明する。
【００３１】
図３は、ボリュームレンダリングの原理を説明するための図である。図３には、三次元空間１００（データ空間）がＸＹＺの直交座標系で示されている。三次元空間１００内の全座標について、各座標ごとにデータ（三次元の速度ベクトルと輝度）が決定されていることは先に説明したとおりである。
【００３２】
ボリュームレンダリングにおいては、通常、三次元空間１００（データ空間）の外側に仮想的に視点ＶＰが設定され、三次元空間１００を間に介して、視点ＶＰと反対側に二次元平面としてのスクリーン１１０が仮想的に設定される。その視点ＶＰを基準として複数のレイ（透視線）１１２が定義される。レイ１１２は、三次元空間１００を貫通するように設定される。そして、三次元空間１００内において、レイ１１２上には複数のボクセル（黒丸）からなるボクセル列が対応付けられる。
【００３３】
なお、各ボクセルのデータ（ボクセル値）は、そのボクセルに対応した座標の速度ベクトルや輝度から算出される。各ボクセルの位置に合致する座標が無ければ、そのボクセルの近傍にある座標の速度ベクトルや輝度を用いて、補間処理などにより、そのボクセルのボクセル値が算出される。本実施形態におけるボクセル値については後に説明する。
【００３４】
ボリュームレンダリングにおいては、視点ＶＰ側からレイ１１２に沿って並ぶ複数のボクセル値を対象として、レンダリング演算（ボクセル演算）が逐次的に実行される。そして、最終のレンダリング演算の結果として画素値が決定される。その画素値がスクリーン上における当該レイ１１２に対応する画素Ｐにマッピングされる。
【００３５】
図３において、スクリーン１１０は、Ｘｓ，Ｙｓの座標系で示されている。つまり、スクリーン１１０内の各座標はＸｓ，Ｙｓの座標値で定義される。このスクリーン１１０内の各座標ごとに１本のレイ１１２が設定されており、上述にように各レイ１１２ごとに得られた画素値をスクリーン１１０上にマッピングすることにより、スクリーン１１０上に三次元画像が形成される。なお、複数のレイ１１２が互いに平行の場合もあるが、複数のレイ１１２が互いに非平行となる場合もある。
【００３６】
レンダリング演算の好適な例では、ｉ番目のボクセルのデータ（ボクセル値）をｅｉとし、その際のオパシティ（不透明度）をαｉ（但し０≦αｉ≦１）とし、ｉ番目の演算結果をＣＯＵＴｉとし、ｉ番目の演算における入力値（（ｉ−１）番目の演算結果）をＣＩＮｉとすると、次式で示すことができる。もちろん、次式以外の公知の演算が用いられてもよい。
【００３７】
【数１】
2012045253.tif 000003
【００３８】
例えば数１式を利用してレイ１１２に沿って逐次的にレンダリング演算を行う場合においては、それと並行して各オパシティαｉが積算され、その積算値が所定値（例えば１）以上となった時点において、そのレイ１１２についての演算が終了する。また、最終のボクセル値についての演算が終了した場合にもそのレイ１１２についての演算は終了する。さらに、本実施形態においては、後に詳述する終了条件に基づいてレイ１１２についての演算を終了させてもよい。
【００３９】
なお、本実施形態に係る超音波診断装置は、上記のようなボリュームレンダリングを利用することが特に好適であるものの、各レイ１１２に沿って逐次演算を行う他のレンダリング処理が適用されてもよい。例えば、積算法、最大値や最小値の検出法などが適用されてもよい。
【００４０】
図１に戻り、レイ設定部５２は、視点ＶＰ（図２）からデータ空間内を通る複数のレイ（図２の符号１１２）を設定する。視点ＶＰは、任意の位置に設定される。例えば、図示しない操作デバイスを利用してユーザが所望の位置に視点ＶＰを設定する。また、スクリーン（図２の符号１１０）は、データ空間を間に挟むように、視点ＶＰの反対側に設定される。レイ設定部５２は、設定された視点ＶＰからデータ空間内を通りスクリーンに達する複数のレイを設定する。
【００４１】
複数のレイが設定されると、ボクセル値算出部５０は、設定された各レイに沿って並ぶ複数のボクセルに関して、各ボクセルごとにボクセル値を算出する。ボクセル値算出部５０は、ボクセル値として各ボクセルごとに速度値と輝度値を算出する。
【００４２】
まず、各ボクセルの速度値は、そのボクセルに対応した座標の速度ベクトルに基づいて算出される。複数のボクセルは三次元空間内（データ空間内）にあるため（図３参照）、三次元空間内における各ボクセルの位置に対応した座標の速度ベクトル（図２）が参照される。そして、各ボクセルの速度値として、そのボクセルに対応した速度ベクトルのレイ方向成分が算出される。つまり、対象となるボクセルが載るレイとそのボクセルに対応した速度ベクトルとの間の角度をθとし、その速度ベクトルの大きさをＶとすると、そのボクセルの速度値は、速度値＝Ｖ×ｃｏｓθと算出される。なお、視点に向かう方向を速度の正方向とし、視点から遠ざかる方向を速度の負方向とする。
【００４３】
また、各ボクセルの輝度値は、そのボクセルに対応した座標の輝度（エコーデータの大きさ）に基づいて算出される。つまり、三次元空間内における各ボクセルの位置に対応した座標の輝度（エコーデータの大きさ）が参照され、そのエコーデータの大きさがそのボクセルの輝度値とされる。
【００４４】
こうして、ボクセル値算出部５０により、各レイに沿って並ぶ複数のボクセルの各々についての速度値と輝度値が算出されると、レンダリング処理部６０において、ボクセル値として算出された速度値と輝度値を利用して、各レイについてのレンダリング演算が実行される。
【００４５】
図４は、レンダリング処理部における演算を説明するための図である。レンダリング処理部（図１の符号６０）は、ボクセル値として算出された速度値と輝度値を利用して、各レイごとにレンダリング演算を実行する。
【００４６】
図４（Ｉ）は、１本のレイについてのレンダリング演算を示しており、１本のレイについて、輝度値を対象としたレンダリング演算（輝度レンダリング）と、正の速度値を対象としたレンダリング演算（正レンダリング）と、負の速度値を対象としたレンダリング演算（負レンダリング）が実行される。
【００４７】
輝度レンダリングにおいては、レイに沿って並ぶ複数のボクセルについて、視点に近い側のボクセル１から順に、ボクセル値として輝度値を利用して、レンダリング演算（ボクセル演算）が実行される。
【００４８】
正レンダリングにおいては、レイに沿って並ぶ複数のボクセルについて、視点に近い側のボクセル１から順に、ボクセル値として速度値を利用して、レンダリング演算（ボクセル演算）が実行される。但し、正レンダリングにおいては、正の速度値つまり視点に向かう方向の速度値のみが演算の対象とされ、速度値がゼロまたは負となるボクセルについては演算の対象から除かれる。なお、速度値がゼロまたは負となるボクセルのボクセル値を０として演算を進めることにより、結果的に演算結果に反映させないようにしてもよい。
【００４９】
一方、負レンダリングにおいては、負の速度値つまり視点から遠ざかる方向の速度値のみが演算の対象とされ、速度値がゼロまたは正となるボクセルについては演算の対象から除かれる。なお、速度値がゼロまたは正となるボクセルのボクセル値を０として演算を進めることにより、結果的に演算結果に反映させないようにしてもよい。
【００５０】
図４（ＩＩ）は、（Ｉ）に示したレンダリング演算の終了判定を示しており、（Ｉ）のレンダリング演算と並行して実行される処理を示している。この終了判定においては、レイに沿って並ぶ複数のボクセルに関する速度値と輝度値が複合的に累積処理され、累積処理の結果が予め設定された閾値に達した時点がレンダリング演算の終了時点とされる。
【００５１】
具体的には、例えば、各ボクセルに関する速度値と輝度値の絶対値が利用される。そして、各ボクセルに関する輝度値の絶対値と速度値の絶対値が累積処理の対象とされ、（Ｉ）に示すレンダリング演算に並行するように、視点に近い側のボクセル１から順に絶対値の加算（累積処理）が進行される。図４（ＩＩ）に示す例においては、ボクセル１の輝度値の絶対値と正の速度値の絶対値が加算され、さらにその加算結果に対して、ボクセル２の輝度値の絶対値と正の速度値の絶対値が加算され、さらにその加算結果に対して、ボクセル３の輝度値の絶対値と負の速度値の絶対値が加算される。こうして、累積処理が進められ、ボクセルｉにおいて累積処理の結果が予め設定された閾値に達すると、（Ｉ）に示す各レンダリング演算がボクセルｉの段階で終了され、その段階における各レンダリング演算の値が演算結果とされる。
【００５２】
図１に戻り、レンダリング処理部６０は、複数のレイの各々について３つのレンダリング演算を実行し、各レイごとに演算結果（Ｅ）と演算結果（＋Ｖ）と演算結果（−Ｖ）を得る（図４参照）。そして、画像形成部７０は、各レイごとに得られる３つの演算結果に基づいて、スクリーン上におけるそのレイに対応した画素Ｐ（図３参照）の画素値を決定する。各画素Ｐの画素値は、ＲＧＢ（赤緑青）の３原色に対応したＲ信号値とＧ信号値とＢ信号値により決定される。これらＲ信号値とＧ信号値とＢ信号値は、各レイごとに３つのレンダリング演算の結果として得られる演算結果（Ｅ）と演算結果（＋Ｖ）と演算結果（−Ｖ）に基づいて決定される。
【００５３】
図５は、レンダリング演算の結果に基づいた３原色の信号値を示す図である。演算結果（＋Ｖ）は、正の速度値を対象としたレンダリング演算の演算結果であり、演算結果（−Ｖ）は、負の速度値を対象としたレンダリング演算の演算結果であり、演算結果（Ｅ）は輝度値を対象としたレンダリング演算の演算結果である（図４参照）。
【００５４】
演算結果（＋Ｖ）と演算結果（−Ｖ）が共にゼロの場合には、Ｒ信号値とＧ信号値とＢ信号値の全ての信号値を演算結果（Ｅ）とする。一方、演算結果（＋Ｖ）と演算結果（−Ｖ）の少なくとも一方がゼロでない場合には、Ｒ信号値を演算結果（＋Ｖ）としてＢ信号値を演算結果（−Ｖ）とし、Ｇ信号値は０または０以外の固定値とする。
【００５５】
これにより、演算結果（＋Ｖ）つまり視点に向かってくる正の速度値に関する演算結果が赤（Ｒ）を基調とした暖色系で画素に反映され、演算結果（−Ｖ）つまり視点から遠ざかる負の速度値に関する演算結果が青（Ｂ）を基調とした寒色系で画素に反映される。
【００５６】
図１に戻り、画像形成部７０は、各レイごとに得られる３つの演算結果に基づいて、スクリーン上におけるそのレイに対応した画素の画素値（Ｒ信号値，Ｇ信号値，Ｂ信号値）を決定する（図５参照）。こうして、仮想的なスクリーン上における全ての画素の画素値が決定され、このスクリーン上に形成される画像が三次元空間に関する表示画像としてモニタなどに表示される。なお、図１に示した超音波診断装置内の各部は、制御部８０により制御される。
【００５７】
以上、本発明の好適な実施形態を説明したが、上述した実施形態によれば、三次元の速度ベクトルに基づいて各ボクセルごとにレイの方向の速度値が算出され、その算出において、視点に向かう方向が正の速度値とされ視点から遠ざかる方向が負の速度値される。つまり、レイの起点となる視点ＶＰが任意の位置に変更されたとしても、常に視点を基準として速度の正負が決定されている。これにより、例えば、視点に向かってくる血流等が赤を基調とした暖色系で表現され、視点から遠ざかる血流等が青を基調とした寒色系で表現され、視点の位置に関わらずカラードプラと同様な色の表現を維持することができる。
【００５８】
なお、上述した実施形態は、あらゆる点で単なる例示にすぎず、本発明の範囲を限定するものではない。本発明は、その本質を逸脱しない範囲で各種の変形形態を包含する。例えば、図１では、３つの３Ｄプローブ１１，１２，１３を利用した実施形態を示したが、４つ以上の複数の３Ｄプローブが利用されてもよい。その場合には、例えばそれら複数の３Ｄプローブの中から３つの３Ｄプローブが選択的に利用されて、図２以降に説明した処理が実現される。
【符号の説明】
【００５９】
１１，１２，１３  ３Ｄプローブ、２１，２２，２３  ビームフォーマ、４０  ３Ｄデータ生成部、５０  ボクセル値算出部、５２  レイ設定部、６０  レンダリング処理部。
【図１】
2012045253.tif 000004
【図２】
2012045253.tif 000005
【図３】
2012045253.tif 000006
【図４】
2012045253.tif 000007
【図５】
2012045253.tif 000008