(57)【要約】【課題】超音波画像処理装置において、深さに依存した色付け条件を定める関数に従ってカラー三次元画像を生成する場合に、関数に含まれるパラメータを設定する際のユーザーの負担を軽減する。【解決手段】ボリュームレンダリングの結果として生じる深さマップ２００に基づいて、関数に含まれる深さパラメータが自動的に決定される。具体的には、深さマップ２００における、平均値、最小値（最浅点の深さ）又は最大値（最深点の深さ）に基づいて、深さパラメータが自動的に決定される。ユーザー指定点の深さに基づいてパラメータが自動的に決定されてもよい。関数に含まれる範囲パラメータ（混合率を規定するパラメータ）が自動的に決定されてもよい。【選択図】図２０
2018000780.tif 000002
【特許請求の範囲】
【請求項１】
生体内の三次元空間から取得されたボリュームデータに対してボリュームレンダリングを適用する手段であって、前記ボリュームレンダリングの結果として輝度マップ及び深さマップを生成するレンダリング手段と、
深さに応じた色付け条件を定める関数に従って、前記輝度マップ及び前記深さマップからカラー三次元画像を生成するカラー処理手段と、
前記深さマップに基づいて、前記関数に含まれる少なくとも１つのパラメータを決定するパラメータ決定手段と、
を含むことを特徴とする超音波画像処理装置。
【請求項２】
請求項１記載の装置において、
前記カラー処理手段は、
前記カラー三次元画像における画素ごとに、当該画素の輝度に基づいてメインカラー及びサブカラーを決定する手段と、
前記関数に対して当該画素に対応する深さを与えることにより決定される混合割合に基づいて、前記メインカラー及び前記サブカラーを混合することにより、前記画素のカラーを決定する手段と、
を含み、
前記パラメータ決定手段は、
前記深さマップに基づいて代表深さを特定する手段と、
前記代表深さに基づいて、前記メインカラーと前記サブカラーとの混合に関わる深さパラメータを決定する手段と、
を含むことを特徴とする超音波画像処理装置。
【請求項３】
請求項２記載の装置において、
前記代表深さには、前記深さマップの全部又は一部についての平均値、最小値及び最大値の内の少なくとも１つが含まれる、
ことを特徴とする超音波画像処理装置。
【請求項４】
請求項２記載の装置において、
前記カラー三次元画像上において又は前記ボリュームデータから生成される断層画像上においてユーザーにより指定点を指定するための手段を含み、
前記代表深さには、前記深さマップにおける前記指定点に対応する深さが含まれる、
ことを特徴とする超音波画像処理装置。
【請求項５】
請求項２記載の装置において、
前記パラメータ決定手段は、前記メインカラーと前記サブカラーとの混合が生じる深さ範囲に関わる範囲パラメータを決定する手段を含む、
ことを特徴とする超音波画像処理装置。
【請求項６】
請求項１記載の装置において、
前記ボリュームレンダリングでは、前記ボリュームデータに対して複数のレイが設定され、レイごとに終了条件が満たされるまで当該レイに沿って出力光量演算が繰り返し実行され、
前記深さマップを構成する各深さは前記各レイにおいて終了条件が満たされた時点での深さである、
ことを特徴とする超音波画像処理装置。
【請求項７】
超音波画像処理装置に超音波画像処理方法を実行させるためのプログラムであって、
前記超音波画像処理方法は、
生体内の三次元空間から取得されたボリュームデータに対してボリュームレンダリングを適用する手段であって、前記ボリュームレンダリングの結果として輝度マップ及び深さマップを生成する工程と、
深さに依存した色付け条件を定める関数に従って、前記輝度マップ及び前記深さマップからカラー三次元画像を生成する工程と、
前記深さマップに基づいて前記関数に含まれる少なくとも１つのパラメータを決定する工程と、
を含むことを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は超音波画像処理装置に関し、特に、超音波ボリュームデータに基づいてカラー三次元画像を形成する装置に関する。
【背景技術】
【０００２】
超音波画像処理装置として、超音波診断装置や、超音波診断装置で取得された超音波ボリュームデータ（以下、単に「ボリュームデータ」という。）を処理する情報処理装置があげられる。超音波画像処理装置においては、ボリュームデータに対してボリュームレンダリングを適用することにより、三次元超音波画像（対象組織を三次元表現した画像）が形成される。
【０００３】
ボリュームレンダリングの実行に際しては、ボリュームデータに対して、視点から出る複数のレイ（視線）が設定される。個々のレイごとに、つまりカラー三次元画像を構成する個々の画素ごとに、レンダリング演算が実行されて画素値が求められる。より詳しくは、開始点から奥行き方向に沿って、ボクセル単位で、オパシティ（不透明度）を利用した出力光量演算（ボクセル演算）が繰り返し実行される。所定の終了条件が満たされた時点での出力光量つまり演算終了点での出力光量が画素値（輝度）を構成する。所定の終了条件としては、オパシティ累積値が１又はそれに近い所定値に到達したことを内容とする第１の終了条件、レイ上の最後のサンプル点まで出力光量演算が到達したことを内容とする第２の終了条件、等が知られている。オパシティ累積値は、一般に、出力光量演算ごとにそこで利用されたオパシティを順次加算することによって求められるものである。出力光量演算式としては各種の式が知られている。複数のレイについて求められた複数の輝度は輝度マップとして観念され、それが三次元超音波画像を構成する。
【０００４】
上記のような一般的なボリュームレンダリング法とは異なる特別なボリュームレンダリング法として、輝度マップと共に深さマップ（距離マップ）を生成し、それらのマップに基づいてカラー三次元画像を形成するボリュームレンダリング法が知られている。その方法では、レイ単位でのレンダリング演算に際し、終了点での出力光量（輝度）と、終了点の深さと、が特定される。これによって、複数のレイに対応する複数の輝度からなる輝度マップと、複数のレイに対応する複数の深さからなる深さマップ（距離マップ）と、が生成される。深さマップは、典型的には、対象組織の三次元形態を表す。例えば、終了点の深さに応じてカラー（例えばRＧＢ値の組み合わせ）を変化させることにより、カラー三次元画像上において奥行き感を増大させることが可能である。
【０００５】
特許文献１、２は、三次元超音波画像の形成に際して、終了点の深さに応じた色付けを提案するものである。特許文献３には、ボリュームデータに対して設定される三次元関心領域が示されている。その関心領域はクリッピング面としてのレンダリング開始面（以下「開始面」と称する。）を有する。その開始面は変形可能な面である（更に高さ方向への並行移動や傾斜運動も可能な面である）。例えば、胎児と子宮壁（胎盤を含む）との間の隙間の形状に合わせて開始面を変形させれば、子宮壁を画像化対象から除外しつつ、胎児だけをレンダリング対象とすることが可能である。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特許第５１１７４９０号公報
【特許文献２】特許第５５２５９３０号公報
【特許文献３】特開２０１１−８３４３９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
上記のように、輝度マップ及び深さマップに基づくカラー三次元画像の生成に際しては、一般に、輝度に基づいて決定されたメインカラー及びサブカラーの混合割合を終了点の深さに応じて変化させる「関数」が利用される。その関数は、混色関数、重み付け合成関数又はブレンド関数として機能するものである。その他の関数が利用されることもある。
【０００８】
しかしながら、上記関数に含まれるパラメータをユーザーが設定、調整して所望のカラー三次元画像を得るのには一般に熟練を要する。結果として所望のカラー三次元画像が得られたとしても、それまでにユーザーに大きな負担が生じてしまう。
【０００９】
本発明の目的は、深さに応じて変化する色付け条件を定める関数に従ってカラー三次元画像を生成する場合に、その関数に含まれるパラメータをユーザーに格別な負担をかけることなく決定できるようにすることにある。あるいは、状況に応じてパラメータが最適化されるようにすることにある。
【課題を解決するための手段】
【００１０】
本発明に係る超音波画像処理装置は、生体内の三次元空間から取得されたボリュームデータに対してボリュームレンダリングを適用する手段であって、前記ボリュームレンダリングの結果として輝度マップ及び深さマップを生成するレンダリング手段と、深さに応じた色付け条件を定める関数に従って、前記輝度マップ及び前記深さマップからカラー三次元画像を生成するカラー処理手段と、前記深さマップに基づいて、前記関数に含まれる少なくとも１つのパラメータを決定するパラメータ決定手段と、を含むことを特徴とする。
【００１１】
上記構成によれば、深さマップに基づいて、関数に含まれる１又は複数のパラメータが決定される。よって、パラメータをユーザーが直接的に指定しなければならない場合に比べて、ユーザーの負担を軽減できる。あるいは、対象組織の形態その他の状況に応じてパラメータを定められるという利点が得られる。複数のパラメータの全部が自動的に決定されてもよいが、それらの一部が自動的（又は半自動的）に決定されてもよい。望ましくは、関数の内容又は作用を示すガイド像が表示される。ガイド像の観察を通じて、自動的に決定されたパラメータを確認できるように構成するのが望ましい。
【００１２】
ボリュームレンダリング法としては各種の方法を利用することが可能である。但し、通常のボリュームレンダリング法では輝度マップしか生成されないため、輝度マップと共に深さマップ（距離マップ）が生成されるように、公知のボリュームレンダリング法を修正又は拡張して利用するのが望ましい。望ましくは、対象組織は子宮内の胎児であり、特に望ましくは、胎児の顔面である。
【００１３】
望ましくは、前記カラー処理手段は、前記カラー三次元画像における画素ごとに、当該画素の輝度に基づいてメインカラー及びサブカラーを決定する手段と、前記関数に対して当該画素に対応する深さを与えることにより決定される混合割合に基づいて、前記メインカラー及び前記サブカラーを混合することにより、前記画素のカラーを決定する手段と、を含み、前記パラメータ決定手段は、前記深さマップに基づいて代表深さを特定する手段と、前記代表深さに基づいて、前記メインカラーと前記サブカラーとの混合に関わる深さパラメータを決定する手段と、を含む。
【００１４】
メインカラーは、望ましくは、視点に近い組織を表現するためのカラーであり、一方、サブカラーは、メインカラーに混合されるカラーであって、望ましくは、視点から遠い組織を表現する（あるいは混色により奥行き感を出す）ためのカラーである。深さに応じて、メインカラーに対するサブカラーの混合割合を変化させることにより、立体感や奥行き感を増大できる。
【００１５】
望ましくは、前記代表深さには、前記深さマップの全部又は一部についての平均値、最小値及び最大値の内の少なくとも１つが含まれる。望ましくは、前記カラー三次元画像又は前記ボリュームデータから生成される断層画像上においてユーザーにより指定点を指定するための手段を含み、前記代表深さには、前記深さマップにおける前記指定点に対応する深さが含まれる。望ましくは、前記パラメータ決定手段は、前記メインカラーと前記サブカラーとの混合に関わる範囲パラメータを決定する手段を含む。例えば、深さのヒストグラムに基づいて範囲パラメータが決定されてもよいし、深さパラメータに基づいて範囲パラメータが決定されてもよい。
【００１６】
望ましくは、前記ボリュームレンダリングでは、前記ボリュームデータに対して複数のレイが設定され、レイごとに終了条件が満たされるまで当該レイに沿って出力光量演算が繰り返し実行され、前記深さマップを構成する各深さは前記各レイにおいて終了条件が満たされた時点での深さである。深さは、開始点から終了点までの距離、基準面（基準点）から終了点までの距離、等である。
【００１７】
本発明に係る超音波画像処理方法は、生体内の三次元空間から取得されたボリュームデータに対してボリュームレンダリングを適用する工程であって、前記ボリュームレンダリングの結果として輝度マップ及び深さマップを生成する工程と、深さに依存した色付け条件を定める関数に従って、前記輝度マップ及び前記深さマップからカラー三次元画像を生成する工程と、前記深さマップに基づいて前記関数に含まれる少なくとも１つのパラメータを決定する工程と、を含むことを特徴とする。超音波画像処理方法は超音波画像処理装置において実行されるプログラムによって構成されてもよい。そのプログラムはネットワークを介して又は記憶媒体を介して超音波画像処理装置へインストールされる。
【発明の効果】
【００１８】
本発明によれば、深さに応じて変化する色付け条件を定める関数に従ってカラー三次元画像を生成する場合に、その関数に含まれるパラメータをユーザーの格別な負担なく決定できる。あるいは、状況に応じて最適なパラメータを決定できる。
【図面の簡単な説明】
【００１９】
【図１】本発明に係る超音波画像処理装置の実施形態を示すブロック図である。
【図２】輝度マップ及び距離マップの生成を示す図である。
【図３】基準面を有しない従来例を示す図である。
【図４】基準面の第１例を示す図である。
【図５】基準面の第２例を示す図である。
【図６】基準面の第３例を示す図である。
【図７】輝度マップ及び深さマップの生成を説明するためのフローチャートである。
【図８】カラー処理部の構成例を示す図である。
【図９】混色関数（重み付け合成関数）を示す図である。
【図１０】表示処理部の動作例を示すフローチャートである。
【図１１】ガイド像の第１例を示す図である。
【図１２】ガイド像の第２例を示す図である。
【図１３】ガイド像の第３例を示す図である。
【図１４】ガイド像の第４例を示す図である。
【図１５】ガイド像の第５例を示す図である。
【図１６】ガイド像の第６例を示す図である。
【図１７】ガイド像の第７例を示す図である。
【図１８】ガイド像の第８例を示す図である。
【図１９】パラメータ決定方法を説明するための図である。
【図２０】深さパラメータ決定方法の第１例を示す図である。
【図２１】深さパラメータ決定方法の第２例を示す図である。
【図２２】最浅点の深さと最深点の深さとに基づく深さパラメータの決定を示す図である。
【図２３】深さパラメータ決定方法の第３例を示す図である。
【図２４】深さパラメータ決定方法の第４例を示す図である。
【図２５】混色関数の変化を示す図である。
【図２６】範囲パラメータ決定方法の第１例を示す図である。
【図２７】範囲パラメータ決定方法の第１例における混色関数の変化を示す図である。
【図２８】範囲パラメータ決定方法の第２例を示す図である。
【図２９】範囲パラメータ決定方法の第２例における混色関数の変化を示す図である。
【図３０】カラー演算部の変形例を示す図である。
【発明を実施するための形態】
【００２０】
以下、本発明の好適な実施形態を図面に基づいて説明する。
【００２１】
（１）装置構成
図１には、本発明に係る超音波画像処理装置の好適な実施形態が示されている。この超音波画像処理装置は、本実施形態において超音波診断装置である。超音波画像処理装置が、画像処理を実行する情報処理装置（コンピュータ）によって構成されてもよい。その場合、超音波診断装置で取得されたボリュームデータが情報処理装置へ渡される。
【００２２】
図１に示す超音波診断装置は医療機関に設置され、生体に対する超音波の送受波により超音波画像を形成する装置である。本実施形態においては、以下に詳述するように、超音波画像として、カラー三次元画像が生成されている。それは、例えば、子宮内の胎児の三次元形態を輝度及び色相の変化によって立体的に表現した画像である。
【００２３】
図１において、３Ｄプローブ１０は、生体内の三次元空間に対して超音波の送受波を行うことによりボリュームデータを取得する送受波器である。３Ｄプローブ１０は、例えば、２Ｄアレイ振動子を有する。２Ｄアレイ振動子は二次元配列された複数の振動素子からなるものである。２Ｄアレイ振動子によれば、超音波ビームを電子的に二次元走査することが可能である。電子走査方式としては、電子セクタ走査方式等が知られている。１Ｄアレイ振動子を機械的に走査することによってボリュームデータを取得してもよい。一般に、１つのボリュームデータは複数の走査面データ（フレームデータ）により構成され、各走査面データは複数のビームデータにより構成される。各ビームデータは深さ方向に並ぶ複数のエコーデータにより構成される。
【００２４】
送受信部１２は電子回路であり、それは送信ビームフォーマー及び受信ビームフォーマーとして機能する。すなわち、送信時において、送受信部１２から３Ｄプローブ１０へ複数の送信信号が供給され、これにより送信ビームが形成される。受信時において、生体内からの反射波が３Ｄプローブ１０によって受波されると、３Ｄプローブ１０から複数の受信信号が送受信部１２へ出力される。送受信部１２では複数の受信信号を整相加算し、これにより受信ビームに相当するビームデータを生成する。
【００２５】
送受信部１２からビームデータ処理部１４へ複数のビームデータが順次出力される。ビームデータ処理部１４は、検波器、対数変換器等を備える電子回路である。そこで処理された各ビームデータが座標変換部１６へ送られる。
【００２６】
座標変換部１６は、座標変換機能、補間機能、等を有する電子回路である。座標変換部１６は、各ビームデータを構成する個々のエコーデータを三次元記憶空間にマッピングする処理を実行する。実際には、三次元座標ごとに補間データを生成しそれをマッピングしている。走査面に対応するフレームデータ単位で座標変換を実行し、それを繰り返すことによりボリュームデータを生成するようにしてもよい。その場合には二次元スキャンコンバータが利用される。いずれにしても、本実施形態においては、座標変換部１６が有するメモリ１８内にボリュームデータが格納される。なお、そのようなボリュームデータから切り出される面データに基づいて二次元断層画像（Ｂモード画像）が形成される。この場合、面データは、後述するレイ群中の特定のレイ列に相当する。この場合、面データは深さ軸を有する。
【００２７】
ボリュームレンダリング部２０は、三次元関心領域設定手段、レンダリング手段、等として機能する。それは電子回路又はソフトウエア機能として実現される。三次元関心領域（3D-ROI）は、ボリュームデータの中で画像化したい部分を切り出すための領域であり、逆に言えば、画像化したくない部分を除外した領域である。三次元関心領域の形態や位置が自動的に又はユーザー入力により定められる。例えば、子宮内の胎児の三次元画像を形成したい場合、胎児と子宮壁（胎盤）との間の隙間（羊水中）に、クリッピング面としての開始面が設定されるように、三次元関心領域が設定される。その場合において、視点に対してボリュームデータを相対的に回転させれば、それに伴って三次元関心領域に対してボリュームデータが相対的に回転する。クリッピング面つまり開始面は、デフォルト状態では平面であるが、上記特許文献３に記載されているように、それは例えば凸面又は凹面に変形可能な面である。更に、その高さ、傾きを変更することも可能である。
【００２８】
本実施形態においては、ボリュームレンダリング部２０が、後述する基準面の設定手段として機能する。基準面は、典型的には、開始面よりも視点側に設けられる平面であり、後述するレイ群に対して直交する平面である。基準面については後に図３乃至図６を用いて詳述する。
【００２９】
ボリュームレンダリング部２０は、ボリュームレンダリング法に基づいて、三次元関心領域つまり対象データに対してレイ群を設定し、個々のレイごとにレンダリング演算を実行する。複数のレイは本実施形態において互いに平行であり、それらは視点（又はそれに相当する面）から奥行き方向に延びる。複数のレイは三次元画像を構成する複数の画素に対応する。ボリュームレンダリング部２０は、レイ単位で、開始点（レンダリング開始点）から終了点（レンダリング終了条件が満たされた点）までにわたって、一定のピッチで、出力光量演算（ボクセル演算）を繰り返し実行する。より正確に言えば、所定の終了条件が満たされるまで、演算点を順次シフトさせながら、オパシティ（不透明度）を利用した出力光量演算が繰り返し実行される。
【００３０】
所定の終了条件には、第１の終了条件と第２の終了条件とがある。出力光量演算で参照するオパシティの累積値が所定値（例えば０．９８）を超える場合に、第１の終了条件が満たされたと判断される。一方、出力光量演算が終了面まで到達した場合に第２の終了条件が満たされたと判断される。オパシティ関数を適宜設定しておけば、例えば、胎児表面又はその近傍で、つまり胎児の表層で、第１の終了条件が満たされる確率が高まる。その結果、専ら胎児表面が表された三次元画像を形成することが可能である。所定の終了条件が満たされた時点での出力光量、つまり終了点での出力光量、がそのレイに対応する輝度とされる。その輝度がメモリ２２上にマッピングされる。複数のレイに対応する複数の輝度のマッピング結果として、メモリ２２上に輝度マップが構成される。
【００３１】
一方、出力光量演算に際しては、上記のように、オパシティが順次加算されており、その累積値が所定値に到達して第１の終了条件が満たされた時点で、あるいは、レイ上の演算点が三次元関心領域の終了面まで到達して第２の終了条件が満たされた時点で、終了点の深さが特定される。本実施形態では、その深さとして、基準点（基準面）から終了点までの長さが演算される。具体的には、終了条件が満たされるまでの出力光量演算回数に対し、レイ上でのサンプル点ピッチを乗算することにより、開始面から終了点までの第１の距離が求められる。一方、当該レイ上において又はその延長線上において、基準面から開始面までの長さが第２の距離として求められる。第１の距離と第２の距離とを加算することにより、カラー処理で参照される深さ（基準面を基準として規格化された深さ）が求められる。その深さがメモリ２４上にマッピングされる。複数のレイに対応する複数の深さにより、メモリ２４上に深さマップが構成される。それは、輝度マップと同様、三次元画像を構成する画素配列と同じ配列を有するものである。本実施形態では、上記のように、深さマップを構成する個々の深さが基準面からの深さつまり規格化された距離である。但し、開始面から終了点までの距離あるいは他の距離を上記深さとしてメモリ２４上にマッピングするようにしてもよい。ボリュームレンダリングについては後に図２及び図７を用いて詳述する。
【００３２】
表示処理部３０は、色付け機能、画像合成機能、等の各種の表示処理機能を備えている。その内で、カラー三次元画像を生成する部分がカラー処理部３２として示されている。また、それに関連する機能としてガイド像生成部２６及びパラメータ決定部２８が示されている。表示処理部３０はプロセッサ又は電子回路により構成される。制御部３６を構成するＣＰＵ及びプログラムによって表示処理部３０が構成されてもよい。
【００３３】
カラー処理部３２は、カラー三次元画像を構成する画素ごとに、輝度Ｉと深さZとから、画素値としてのカラーを決定するものである。例えば、輝度Ｉに応じて決定される暖色系のカラーと、同じく輝度Ｉに応じて決定される寒色系のカラーと、が深さＺに応じて決定される色付け条件（混合率）に従って混合され、その混合結果として当該画素のカラーが決定される。その例において、前者の暖色系のカラーはメインカラーであり、つまり近距離表現用カラーであり、それは例えばオレンジ系のカラーである。後者の寒色系のカラーはサブカラーであり、つまり遠距離表現用カラーであり、例えばブルー系のカラーである。各画素に対応する深さは、典型的には、対象組織の表層の深さを表しているので、浅い位置に存在する表層については暖色系カラーで表現され、深い位置に存在する表層については寒色系カラーの混合割合が高められたカラーで表現される。もちろん、上記であげた具体的なカラーは例示である。
【００３４】
本実施形態では、カラー処理部３２が、輝度に応じてメインカラー及びサブカラーを決定する２つの変換テーブル（変換関数）を有しており、また、メインカラー及びサブカラーを深さに応じて重み付け合成する合成器を有している。その重み付け合成条件は混色関数によって規定される。混色関数には、深さ方向における混色範囲の中間深さを示す深さパラメータと、深さ方向における混色範囲の大きさを示す範囲パラメータと、が含まれる。但し、それらは例示である。それらに代わる複数のパラメータが混色関数に含まれていてもよい。本実施形態において、混色関数に含まれる個々のパラメータは、カラー三次元画像及び断層画像を観察しながらユーザーによって指定され、あるいは、深さマップ等に基づいて自動的に決定される。カラー処理部３２の構成及び作用については後に図８及び図９を用いて詳述する。
【００３５】
ガイド像生成部２６は、断層画像と共に表示されるガイド像を生成するモジュールである。ガイド像は、後述するように、グラフィック要素として構成され、つまり図形、色見本等として構成され、混色関数の内容又は作用の一部を視覚的に表現した像である。ガイド像も深さ軸（第１の深さ軸）を有する。断層画像は深さ軸（第２の深さ軸）を有する。第２の深さ軸はレイ群に対して平行であり、それは通常、断層画像の左側辺又は右側辺に一致する。第１の深さ軸が第２の深さ軸に平行になるように、ガイド像が断層画像の近傍に又は断層画像の上に表示される。第１の深さ軸が後述するマーカーの運動軸であってもよい。混色関数に含まれる深さパラメータ及び範囲パラメータをユーザーが設定する場合にガイド像が参照される。これによりパラメータ設定を支援できる。カラー三次元画像を観察する場合にもガイド像が参照される。これによりカラー三次元画像として表示された胎児の形態の把握を促進できる。ガイド像については後に図１１乃至図１８を用いて詳述する。
【００３６】
パラメータ決定部２８は、混色関数に含まれる深さパラメータ（具体的には後に示すZs）及び範囲パラメータ（具体的には後に示すΔZ）の内の前者又は両者を自動的に（又は半自動的に）設定するモジュールである。この構成により、ユーザーの負担を軽減でき、また、カラー三次元画像の色合いを最適化することが容易となる。パラメータ自動設定については後に図１９乃至図２９を用いて詳述する。
【００３７】
制御部３６はＣＰＵ及びプログラムによって構成され、制御部３６により、図１に示されている個々の構成の動作が制御される。操作パネル３８は入力装置を構成し、それを利用して三次元関心領域をマニュアル設定したり混色関数に含まれるパラメータを設定したりすることが可能である。クリッピング面としての開始面をマニュアルで設定する場合にも操作パネル３８が利用される。本実施形態では、三次元関心領域よりも視点側に基準面が自動的に設定されているが、その基準面をマニュアルで設定できるように構成してもよい。表示器３４には、カラー三次元画像、断層画像、ガイド像の他、必要に応じて、任意断層画像等が表示される。表示器３４にトリプレーン画像が表示されてもよい。
【００３８】
（２）輝度マップ及び深さマップの生成
図２には、輝度マップ及び深さマップの生成方法が模式的に示されている。但し、図２には一般的な生成方法が示されており、そこには基準面を利用した深さマップの生成方法は示されていない。
【００３９】
三次元関心領域４０に対して互いに並行に複数のレイ４２が設定される。Ｘ−Ｙ−Ｚは直交座標系を表している。Ｚ方向が深さ方向である。符号４１で示す矢印は視点から出る投影方向を示している。その方向に対して複数のレイ４２が平行に設定される。なお、ボリュームデータに対して任意の位置に視点を設定することが可能である。
【００４０】
図示の例では、符号４０Ａがクリッピング面としての開始面を示している。それは平面として描かれている。実際には開始面４０Ａの形状等を自在に変更することが可能である。個々のレイ４２において、開始面４０Ａ上の開始点から奥行き方向（図示の例ではＺ方向）へ出力光量演算が繰り返し実行される。その演算の都度、オパシティが累積加算される。上述した終了条件が満たされた時点で、出力光量演算が止まり、つまり、レンダリング演算が終了する。これによって、レイ４２ごとに終了点４４が特定される。その終了点４４は例えば胎児の表層に相当する。レンダリング演算終了時点での出力光量が輝度としてメモリ上にマッピングされる。複数のレイ４２に対応する複数の輝度のマッピング結果として輝度マップ４６が構成される。それは白黒のボリュームレンダリング画像に相当する。
【００４１】
一方、レイ単位で、終了点４４の深さがメモリ上にマッピングされる。複数のレイ４２に対応する複数の深さのマッピング結果として深さマップ４８が構成される。グレースケール５０に示されているように、図２においては、最小深さから最大深さまでが輝度の変化として表現されている。
【００４２】
輝度マップ４６及び深さマップ４８は、三次元画像を構成する二次元画素配列と同じ二次元画素配列を有する。輝度マップ４６及び深さマップ４８に基づいて、カラー三次元画像が形成される。例えば、上記のように、個々の画素ごとに、輝度に対応する２つのカラーが、深さに応じて重み付け合成され、これにより当該画素のカラーが決定される。これが全画素にわたって繰り返される。深さマップ４８を構成する各深さは、以下に説明するように、基準点から終了点までの距離であるのが望ましい。但し、開始点から終了点までの距離であってもよい。それらの距離は、いずれも視点側から見た対象組織の深さを示すものである。
【００４３】
（３）基準面に基づく深さの演算
図３には基準面を有しない従来例が比較例として示されている。符号５４は胎児を示しており、符号５６は胎盤を示している。符号５７は羊水を示している。符号５２はボリュームデータを示し、符号５８は三次元関心領域を示している（図３には三次元関心領域の断面が現れている）。三次元関心領域５８は、クリッピング面として機能する開始面６０、底面（下面）としての終了面６５、及び、４つの側面（図３には２つの側面６２，６４が現れている）を有するボックス状の形態を有する。開始面６０の凹形変形により、開始面６０が胎児５４と胎盤５６との間の隙間に差し込まれている。このように三次元関心領域５８の形態が定められた上で、それに対して設定される個々のレイごとにレンダリング演算が実行される。
【００４４】
図３において、個々の下向き矢印６３はレンダリング演算経路を示しており、矢印６３の先端が終了点を示している。なお、図３には、複数の下向き矢印６３の中で胎児表層まで到達していないものが含まれているが、それは作図上の都合によるものである。図３に示す例では、基準面が採用されていない。つまり、各レイにおいて、開始面から終了点までの長さが深さとして演算され、その深さがそのまま色付け処理で参照される。それ故、開始面６０が湾曲していることに起因して、深さマップにおいて胎児表面の形態とは関係のない深さ変化が生じてしまう。それは無用な色相変化を引き起こし得るものである。
【００４５】
図４には基準面の第１例が示されている。なお、図３に示した要素と同様の要素には同一符号を付し、その説明を省略する。このことは、後に説明する図５及び図６についても同様である。
【００４６】
図４においては、開始面６０よりも視点に近い側に基準面６６が設定されている。その基準面６６は、複数のレイに直交する平面である。デフォルト状態の開始面６０を基準面と定めてもよい。基準面６６は三次元関心領域５８の全体を覆っている。具体的には、三次元関心領域５８が有する４つの側面の上辺に連接する平面として基準面が構成されている。但し、視点により近い位置に基準面が設定されてもよい。基準面６６が対象組織を横切らないように当該基準面６６の位置を定めるのが望ましい。
【００４７】
個々のレイにおいては、従来同様にレンダリング演算が実行されるが、深さの演算に限っては、基準面６６が基準として利用される。すなわち、各レイにおいて、基準面６６から終了点までの規格化された長さが深さとして演算される。実際には、開始面６０上の開始点から終了点までの第１の距離が演算され、一方、基準面６６から開始面６０までの第２の距離が演算される。そして、第１の距離と第２の距離とを加算することにより、その加算結果として、マッピングされる深さが演算される。もっとも、基準面から終了点までの距離が直接的に求められてもよい。本実施形態では、レンダリング演算範囲とは別に、深さが定義されるようにしたので、開始面の形態に依存しない、規格化された深さをカラー処理で利用できる。
【００４８】
図５には基準面の第２例が示されている。この第２例では、胎盤５６Ａの形態を考慮して、三次元関心領域５８Ａの開始面６０Ａが凸面を構成しており、これにより、開始面６０Ａから見て視点とは反対側に基準面７０が設定されている。但し、基準面７０は、胎児５４にかかっておらず、各レイにおいて演算される深さは、いずれもプラスの深さである（符号７２を参照）。もっとも、図５に示す第２例では、基準面７０が対象組織にかかり易いので、基準面７０がより上方に位置するように、例えば、基準面７０が開始面６０Ａの頂点を通過するように、基準面設定条件を変更するのが望ましい。
【００４９】
図６には基準面の第３例が示されている。この例では、胎盤５６Ｂがかなり複雑な形態を有しており、それに応じて、三次元関心領域５８Ｂの開始面６０Ｂも複雑な形態を有している。開始面６０Ｂの視点側に基準面７４が設定されている。それは平面であるが、各レイに対しては傾斜している。その傾斜角度は小さく、傾斜を原因とする深さ変化量はかなり小さい。しかも一様に傾斜しているので不自然な深さ変化や局所的な深さ変化は生じていない。このように、場合によっては、傾斜した基準面７４を利用することも可能である。但し、傾斜に起因して見かけ上の深さ増減が生じるので、できる限り、各レイに直交する基準面（上記第１例として示した基準面）を採用するのが望ましい。
【００５０】
（４）ボリュームレンダリング
図７に輝度マップ及び深さマップの生成に関する動作の一例が示されている。Ｓ１０では、レイ番号を示すｋが初期化され、Ｓ１２では、距離を示すｄが初期化される。Ｓ１４では、ｋ番目のレイ上において距離ｄにある注目点について、エコー値（ボクセル値）が演算される。例えば、注目点の近傍に存在する８個のエコー値に基づいて三次元線形補間演算により注目点のエコー値が算出される。Ｓ１６では、注目点でのエコー値、それに対応したオパシティ、等に基づく出力光量演算（ボクセル演算）が実行される。これにより注目点での出力光量が求まり、その出力光量が次の注目点での入力光量とされる。このようなアルゴリズムとして各種のものが知られている。Ｓ１６では、オパシティの累積加算も実行される。
【００５１】
Ｓ１８では、所定の終了条件が満たされたか否かが判断される。満たされていない場合、Ｓ２８においてｄがインクリメントされた上で、Ｓ１４からの各工程が繰り返し実行される。一方、Ｓ１８において、終了条件が満たされたと判断された場合、Ｓ１９において、最終的な出力光量が輝度としてメモリ上にマッピングされる。
【００５２】
Ｓ２０では、基準面から開始面までの距離Δｄが計算される。距離Δｄは開始面の位置及び形態に依存する。Ｓ２２では、終了点の距離であるｄと、基準面と開始面との間の距離Δｄと、が加算され、その結果として、規格化された深さＺが求まる。Ｓ２４では深さＺがメモリ上にマッピングされる。Ｓ２６では、ｋが最大値に到達したか否かが判断され、ｎｏであればＳ３０においてｋがインクリメントされた上で、Ｓ１２以降の各工程が繰り返し実行される。
【００５３】
この後、輝度マップ及び距離マップに基づいてカラーの三次元画像が形成される。但し、三次元画像形成に先立って、少なくとも深さマップに対してはＳ２７で平滑化処理が施される。
【００５４】
（５）カラー処理
図８には、図１に示したカラー処理部３２の構成例が示されている。カラー処理部３２は、２つのテーブル８４，８６と、重み付け合成器８８と、を有している。テーブル８４は、フロント側変換テーブルであり、輝度Ｉに応じてメインカラー（ＲＧＢ値）Ｃｆを決定するルックアップテーブルである。一般的には、メインカラーとして、温かみを感じられるオレンジ系のカラーが採用される。つまり、近距離に存在する組織がそのようなカラーによって表現される。テーブル８６は、バック側変換テーブルであり、輝度Ｉに応じてサブカラー（ＲＧＢ値）Ｃｂを決定するルックアップテーブルである。一般的には、サブカラーとして、メインカラーへの混色によって遠近感を増大できるブルー系のカラーが採用される。もっとも、上記のいずれのカラーも例示に過ぎない。
【００５５】
重み付け合成器８８は、輝度Ｉに応じて決定されたカラーCｆ及びカラーＣｂを、深さＺに応じて重み付け合成するモジュールである。深さＺに応じた混色割合は、以下に説明する混色関数（符号８８Ａ参照）によって決定される。重み付け合成器８８から、重み付け合成後のカラー（ＲＧＢ値）Ｃが出力される。もっとも、カラーＣｆに対して重み１．０が与えられる場合にはカラーＣｆがそのままカラーＣとして出力され、同様に、カラーＣｂに対して重み１．０が与えられる場合にはカラーＣｂがそのままカラーＣとして出力される。本実施形態においては、混色関数に含まれる２つのパラメータZs、ΔZを操作することにより、色付け条件が変更される。
【００５６】
図９には混色関数の一例が示されている。横軸は深さＺを示している。縦軸は重みを示しており、具体的には、左側の縦軸はメインカラー（フロントカラー）Cｆに与えられるフロント重みを示しており、右側の縦軸はサブカラー（バックカラー）Cｂに与えられるバック重みを示している。第１の混色関数１００によってメインカラーＣｆに与えられる重みが決定される。同様に、第２の混色関数１０２によってサブカラーＣｂに与えられる重みが決定される。２つの重みの総和は１．０である。よって、一方の混色関数から他方の混色関数が自動的に定まるので、一方の混色関数だけを利用するようにしてもよい。以下においては、必要に応じて２つの混色関数１００，１０２の内で、メインカラー用の混色関数１００だけを明示する。
【００５７】
図９において、横軸方向に３つの区間が設定されている。左側の区間（Ｚ：０−Ｚ１）においてはメインカラーが１００％採用される。中央の区間（Ｚ：Ｚ１―Ｚ２）においてはメインカラーとサブカラーとが深さ応じた混合割合で混合される。当該区間は深さ方向における混色区間あるいは混合範囲と言い得る。右側の区間（Ｚ：Ｚ２−Ｚmax）においてはサブカラーが１００％採用される。混合区間における中間深さが深さパラメータZsであり、混合区間の深さ方向の範囲が範囲パラメータΔZである。本実施形態では、それらのパラメータによって、混色関数１００，１０２の形態つまり混合作用が規定される。例えば、対象組織の表層における平均的な深さ付近に、あるいは、特に形態観察を行いたい部位の深さ付近に、深さパラメータZsが設定される。また、例えば、対象組織の形態の変化の度合いに応じて、範囲パラメータΔZが設定される。
【００５８】
深さパラメータZsを増減すると、混色関数１００，１０２が横軸方向にシフトすることから、深さパラメータZsをシフトパラメータと称することもできる。範囲パラメータΔZを増減すると、カラー三次元画像におけるメインカラー及びサブカラーの割合と、混合カラーの割合とが変化することから、範囲パラメータΔZをコントラストパラメータと称することもできる。なお、他のパラメータによって混色関数が定義されてもよい。また、曲線によって混色関数が定義されてもよい。
【００５９】
カラー三次元画像を観察するだけで、上記パラメータをマニュアルで設定することは難しく、その設定を支援することが望まれる。そこで、本実施形態ではガイド像が生成されており、それが三次元画像及び断層画像と共に表示される。また、パラメータ設定の際のユーザー負担を軽減するためパラメータを自動的に決定する仕組みが設けられている。
【００６０】
それらの概要を図１０に基づいて説明する。図１０には図１に示した表示処理部の動作がフローチャートとして示されている。Ｓ４０ではモードが選択される。パラメータをマニュアルで設定をする場合にはＳ４２が実行され、パラメータを自動的に（又は半自動的に）設定する場合にはＳ５０が実行される。Ｓ４２では、画素ごとに、輝度に応じてフロントカラー及びバックカラーが決定され、それらが深さに応じて重み付け合成される。これにより、当該画素のカラーが決定される。重み付け合成の条件（つまり色付け条件）は混色関数により定められる。重み付け合成に先立って、Ｓ４４において、ユーザーにより、深さパラメータZs及び範囲パラメータΔZが指定され、それらが受け付けられる。それらのパラメータZs,ΔZによって混色関数が定義されることになる。Ｓ４６において、混色関数を象徴する像として、混色関数の内容又は作用の少なくとも一部又は一面を表すガイド像が表示される。それを参照しながらパラメータ設定を行える。ガイド像に影響を与えるパラメータが変更された場合、その変更がガイド像にリアルタイムで反映される。Ｓ４８では、本処理を継続するか否かが判定され、本処理を継続する場合にはＳ４０以降の各工程が繰り返し実行される。
【００６１】
Ｓ４０において、パラメータの自動設定が選択された場合、Ｓ５０において、深さマップに基づいて深さパラメータZsが自動的に決定される。その上で、Ｓ５２において範囲パラメータΔZが例えば深さパラメータZsに基づいて自動的に決定される。但し、範囲パラメータΔZがマニュアルで決定されてもよい。Ｓ５６では、上記Ｓ４６と同様に、決定された深さパラメータZsと範囲パラメータΔZの内の両方又は一方がガイド像の内容に反映される。Ｓ５４では、画素ごとに、輝度に応じてフロントカラー及びバックカラーが決定され、それらが深さに応じて重み付け合成される。これにより当該画素のカラーが決定される。この場合も、重み付け合成の条件が混色関数により定められる。
【００６２】
（６）ガイド像
図１１にはガイド像の第１例が示されている。表示画面１０４上には白黒の断層画像１０６及びカラー三次元画像１０８が表示されている。カラー三次元画像１０８は、三次元関心領域内のボリュームデータに対してボリュームレンダリングを適用し、更にカラー処理を適用することにより生成された画像である。ここでは胎児の顔が立体的に表現されている。断層画像１０６は、三次元空間（ボリュームデータ）の垂直断面を示し、それはレイ列に相当する。この例では、断層画像１０６において、横軸がＸ軸であり、縦軸がＺ軸である。断層画像１０６における左側の垂直辺及び右側の垂直辺はいずれもＺ軸に相当する。断層画像１０６内にはボックスとして三次元関心領域１１０が表示されている。
【００６３】
断層画像１０６の左側端部には、具体的には、三次元関心領域１１０の左側垂直辺の近傍には、ガイド像１１２が表示されている。ガイド像１１２は、混色関数の内容又は作用を示す像である。この第１例では、ガイド像が、深さパラメータZs（つまり混色範囲の中間レベル１２０）を示すマーカー１１４を有している。符号１１６は開始面レベルを表し、符号１１８は終了面レベルを示している。マーカー１１４はそれらの間を垂直方向につまり深さ軸方向に運動するグラフィック要素である。図示の例では、マーカー１１４が三角形の形態を有しているが、それは例示に過ぎない。
【００６４】
ガイド像１１２は運動軸としての深さ軸を有し、その深さ軸は断層画像が有する深さ軸と平行である。操作パネル３８における摘み１２２を操作すると、符号１２６で示すように、深さパラメータZsが変化し、それに伴ってマーカー１１４が上下方向に運動する。同時にカラー三次元画像１０８のカラー態様も変化する。深さ軸上のマーカー１１４によって深さパラメータZsを認識できる。本実施形態によれば、断層画像１０６上に表示されている組織構造を参照しつつ、深さパラメータZsを最適化することが可能である。なお、操作パネル３８には範囲パラメータΔZを調整するための摘み１２４も含まれる。図１１に示す第１例では、範囲パラメータΔZを変更しても、ガイド像１１２の表示態様は変化しないが、カラー三次元画像１０８のカラー表現態様は変化する。この第１例によると、混色関数の内容（特にそこに含まれる主要なパラメータ）及び混色関数が作用する深さ（混色範囲の中間レベル）を断層画像との対比において認識することが可能である。
【００６５】
図１２には、ガイド像の第２例が示されている。なお、図１１に示した要素と同様の要素には同一符号を付しその説明を省略する。このことは図１３以後の各図においても同様である。
【００６６】
図１２において、ガイド像１２８は、深さ方向に伸長した帯状の形態を有する。この例では、ガイド像１２８の上端は開始面レベル１１６に一致しており、ガイド像１２８の下端は終了面レベル１１８に一致している。ガイド像１２８は、大別して、倒立した台形状の上部１２６ａと、台形状の下部１２６ｂと、からなる。それらの間が境界線１２６ｃであり、それは共通斜辺をなしている。境界線１２６ｃは混色関数（フロント重みを規定する混色関数）における混色範囲を表している。境界線１２６ｃの深さ方向の幅が混色範囲の深さ方向のサイズに対応しており、境界線１２６ｃの傾きが混合率つまり勾配に対応している。つまり、ガイド像１２８は混色関数の形態それ自体を表現している。
【００６７】
上部１２６ａはフロント側変換テーブルにより生成されるメインカラー群の内の代表的なメインカラーによって着色されている。下部１２６ｂはバック側変換テーブルにより生成されるサブカラー群の内の代表的なサブカラーによって着色されている。よって、上部１２６ａ及び下部１２６ｂはそれぞれ色見本として機能する。
【００６８】
符号１２６で示されているように、摘み１２２を操作して深さパラメータZｓを変化させると、境界線１２６ｃがその傾斜角度を維持したまま深さ方向へ運動する。一方、符号１３４で示されているように、摘み１２４を操作して範囲パラメータΔZを変化させると、境界線１２６ｃの中間深さが維持されたまま、境界線１２６ｃの両端間の深さ方向のサイズが変化し、同時に傾斜角度が変化する。この第２例によれば、深さ軸上において、２つのパラメータを具体的に認識することが可能である。換言すれば、この第２例によれば、混合関数それ自体の具体的な内容及びそれによる作用の全体を直感的に認識し易い。
【００６９】
図１３には、ガイド像の第３例が示されている。ガイド像１３６はグラフ状の形態を有し、具体的には、横軸１３６ａ，縦軸１３６ｂ及び折れ線１３６ｃを有する。横軸１３６
ａはフロント重み（メインカラーに与える重み）の大きさを示している。縦軸１３６ｂは深さＺを示している。折れ線１３６ｃはメインカラー用の混色関数を示している。サブカラー用の混色関数は、メインカラー用の混色関数から容易に導ける。第４例によれば、混色関数の内容を具体的に認識できる。
【００７０】
図１４には、ガイド像の第４例が示されている。ガイド像１３８は深さ方向に伸長した帯状の形態を有する。ガイド像の上端１３８ａから下端１３８ｂにかけて深さ方向に色が変化する色変化見本が構成されている。実際には、混色範囲（符号１３０で示される深さZ1からが符号１３２で示される深さZ2の間）において、色変化（グラデーション）が生じている。深さZ1よりも上側ではメインカラーが表示され、深さZ2よりも下側ではサブカラーが表示される。このような深さ方向に変化する色変化見本により、混色関数の作用を具体的に認識することができる。
【００７１】
図１５には、ガイド像の第５例が示されている。ガイド像１４０は、色変化見本１４２と、２つのマーカー１４４，１４６と、により構成されている。色変化見本１４２は、深さ方向に伸長した帯状の形態を有し、その上辺１４２ａが開始面深さを表しており、その下辺１４２ｂが終了面深さを表している。色変化見本１４２は、上記のガイド像１３８と同様に、深さ方向に変化するカラーによって構成されている。符号１３０及び符号１３２は、混色範囲の両端の深さZ1，Z2を示している。それらの深さZ1，Z2が２つのマーカー１４４，１４６によって特定されている。各マーカー１４４，１４６は深さZ1，Z2の変化に伴って深さ方向に運動するグラフィック要素である。その形態として三角形以外の形態（例えば単なる水平ライン）を採用してもよい。各マーカー１４４，１４６の位置を調整することにより混色関数を設定することが可能である。その場合、深さパラメータZs及び範囲パラメータΔZが深さZ1，Z2によって定められてもよいし、深さZ1，Z2がパラメータとして混色関数に直接的に与えられてもよい。各マーカー１４４，１４６をポインティングデバイスによって直接的に動かすようにしてもよいし、摘み操作によって動かすようにしてもよい。
【００７２】
図１６には、ガイド像の第６例が示されている。表示画面１４８上にはいわゆるトリプレーンを構成する直交関係にある３つの白黒断層画像１５０，１５２，１５４が表示されている。更に、カラー三次元画像１５６が表示されている。断層画像１５０はＸ軸及びＺ軸を有し、第１垂直断面としてのＡ面を表すものである。断層画像１５２はＹ軸及びＺ軸を有し、第２垂直断面としてのＢ面を表すものである。断層画像１５４はＸ軸及びＹ軸を有し、水平断面としてのＣ面を表すものである。断層画像１５０内には三次元関心領域の近傍にガイド像１５８が表示されている。断層画像１５２内には三次元関心領域の近傍にガイド像１６０が表示されている。このように複数の断層画像と共に複数のガイド像を表示するようにしてもよい。
【００７３】
なお、ガイド像１５８，１６０は、それぞれ、図１２に示したガイド像と同じ態様を有している。もちろん、ガイド像１５８，１６０として、他の態様を採用してもよいし、それらの態様を互いに異ならせてもよい。
【００７４】
図１７には、ガイド像の第７例が示されている。断層画像上には三次元関心領域１６２が示されており、それはクリッピング面として機能する上面１６２Ａを有する。それとは別に、規格化された深さを演算するための基準面１６４が設定されている。ガイド像１６６は左側部分１６６Ｌ及び右側部分１６６Ｒを含む。左側部分１６６Ｌは深さパラメータZsを示すマーカー１６８を有している。深さパラメータZsを示す水平ライン１７４も表示されている。右側部分１６６Ｒは混色範囲の両端の深さZ1,Z2を示す２つのマーカー１７０，１７２を有している。
【００７５】
このようにガイド像を構成する複数の要素を分散的に表示することにより、視認性や見栄えを良好にできる。更に、水平ライン１７４を表示すれば、組織構造と深さパラメータZsとの関係を具体的に特定できるので、深さパラメータZsの調整がし易くなる。
【００７６】
図１８には、ガイド像の第８例が示されている。三次元関心領域１７６の上面１７６Ａが凸面を構成している。この例では基準面が採用されていない。つまり、上面１７６Ａから終了点までの距離が深さとして演算されている。表示エリア内における個々の開始点の位置が揃わなくなることから（湾曲することから）、水平方向の各位置での深さパラメータZsの位置も画面上で湾曲することになる。これを踏まえ、ガイド像１７８には、マーカー１８０を含む左側部分１７８Ｌ、及び、マーカー１８２，１８４を含む右側部分１７８Ｒに加えて、深さパラメータZsを表す曲線１８６、深さZ1を表す曲線１８８及び深さZ2を表す曲線１９０が含まれている。それらの曲線１８６，１８８，１９０は、断層画像を横断している。この第８例によれば、各水平位置において、深さパラメータZs、深さZ1及び深さZ2を正しく認識することが可能である。
【００７７】
以上のように、断層画像と共に（それに重畳して又はその近傍に）ガイド像を表示すれば、断層画像上に表示された組織との関係において、関数の内容又は作用を直感的に認識することが可能である。よって、そのようなガイド像を見ながらパラメータの設定を行えば、ユーザーの負担や混乱を軽減できる。ガイド像はカラー三次元画像を観察する場合においても参考になるものである。
【００７８】
（７）パラメータ自動決定
既に説明したように、本実施形態においては、必要に応じて、混色関数に含まれるパラメータが自動的に決定される。図１９において、符号１９４で示されている深さパラメータZs（以下、単に「深さZs」ともいう。）は、符号１９２で示されている深さマップから決定される。その場合、以下に説明するように、深さマップの全部又は一部における、平均値、最浅点の深さ、最深点の深さ、ユーザー指定点の深さ等として、深さパラメータZsが決定される。符号１９６で示されている範囲パラメータΔZ（以下、単に「範囲ΔZ」ともいう。）がマニュアルで設定されてもよいが、望ましくは、それは、深さパラメータZsから自動的に決定され、あるいは、深さマップから自動的に決定される。三次元画像上での指定点に対応する深さに基づいて範囲パラメータΔZが決定されてもよい。
【００７９】
深さパラメータZs及び範囲パラメータΔZが決定されると、それらによって混色関数１９８が定義される。混色関数１９８は、既に詳述したように、画素の輝度から並列的に求められるメインカラー及びサブカラーを、画素に対応する深さに応じて混合する関数である。以下にパラメータ決定方法の具体例を説明する。
【００８０】
図２０には、深さパラメータ決定方法の第１例が示されている。カラー三次元画像２００内にはユーザー設定により又は自動的に参照エリア２０６が設定されている。図示の例では、垂直中心線２０４及び水平中心線２０２の交点（中心点）を中心として広がる矩形のエリアとして参照エリア２０６が設定されている。参照エリア２０６は、エコーが不安定となる、あるいは、エコーが消失し易い周囲部分を除外したものである。カラー三次元画像２００の全体を参照エリア２０６として設定してもよいが、観察したい部位が支配的に含まれ、観察したくない部位が除外されるように、参照エリアのサイズや位置を調整するのが望ましい。
【００８１】
この第１例では、深さマップ上に参照エリア２０６と同じ参照エリアが設定され、その参照エリア内において深さの平均値が演算される。参照エリア２０６はそのために設定されるエリアであるとも言える。深さマップにおける参照エリア内の平均値がそのまま深さパラメータZsとして設定される。なお、範囲パラメータについては、マニュアルで設定されてもよいし、後述するように、深さパラメータZsに連動する値としてもよい。
【００８２】
参照エリア内の平均値は、画像化したい組織（表層）の平均的な深さと言えるから、そこを深さパラメータZsとすれば、混色範囲を最適化できる可能性を高められる自動設定誤、必要に応じて、マニュアルで微調整してもよい。
【００８３】
図２１には、深さパラメータ決定方法の第２例が示されている。第１例と同様に、カラー三次元画像２００に対して参照エリア２０６が設定されている。この参照エリア２０６と同じ参照エリアが深さマップ上に設定される。そして、その参照エリア内における最も浅い点（最浅点）Ａの深さZaと、その参照エリア内における最も深い点（最深点）Ｂの深さZbと、が自動的に特定される。その上で、例えば次の(1)式により、深さパラメータZsが決定される（符号２１０参照）。但しαは係数である。
【００８４】
Zs=(Zb-Za)×α+Za  ・・・(1)
【００８５】
この構成によれば、参照エリア内における深さの広がりを考慮して、深さパラメータZsを自動的に決定できる。決定された深さパラメータZsはガイド像を通じて確認し得る。
【００８６】
図２２には上記(1)式を構成する要素が示されている。符号２１２は断層画像上における三次元関心領域を示している。それは胎児の頭部２１４を包含している。最浅点Ａの深さZaと、最深点Ｂの深さZbとの間を一定の割合αで内分する深さとして深さZsが決定される。αを固定値とせずにパラメータとし、ユーザーによって又は対象組織によってそれを可変できるようにしてもよい。
【００８７】
図２３には、深さパラメータ決定方法の第3例が示されている。カラー三次元画像２００上には参照エリア２０６が設定されており、また、カラー三次元画像２００上にはユーザー指定された指定点Ａが設定されている。指定点Ａが断層画像上において指定されてもよい。指定点Ａは、例えば、胎児において最も浅いと思われる部位上に指定される。深さマップから指定点に対応する深さZaが特定される。また、深さマップにおける参照エリア（上記参照エリア２０６と同じ参照エリア）内における最深点Ｂが特定され、その深さZbが特定される。そして、上記(1)式に従って、深さパラメータZsが決定される（符号２１６参照）。最深点をユーザーが指定し、最浅点を自動的に決定するようにしてもよい。
【００８８】
図２４には、深さパラメータ決定方法の第４例が示されている。カラー三次元画像２００において最浅点Ａ及び最深点Ｂを自動的に特定することにより、深さマップから最浅点Ａの深さZa及び最深点Ｂの深さZbが決定される。その上で、上記(1)式に従って、深さパラメータZsが決定される（符号２１８参照）。
【００８９】
図２５には、深さパラメータZsの変化に伴う混色関数の変化が示されている。範囲パラメータΔZを維持したまま、深さパラメータZsを変化させると、符号２２０ａ、符号２２０ｂ、符号２２０ｃで示すように、混色関数が平行移動（シフト）する。その場合、混色範囲における重み変化率つまり勾配は不変である。もっとも、以下に説明するように、範囲パラメータΔZを深さパラメータZs等に応じて変化させるようにしてもよい。
【００９０】
図２６には、範囲パラメータ決定方法の第１例が示されている。横軸は深さＺを示しており、縦軸は範囲ΔZを示している。その最大値はΔZmaxである。Zs0は深さパラメータZsの初期値又は標準値を示しており、Zs1は現在の深さパラメータZsを示している。このような前提の下、以下の(2)式に従って、現在の範囲パラメータΔZ１が決定される（符号２２４参照）。以下のβ１は係数である。
【００９１】
ΔZ1=ΔZmax−β1×｜Zs1−Zs0｜  ・・・(2)
【００９２】
上記(2)式において、β1は係数である。但し、混色範囲の下限Z1及び上限Z2が０≦Z1，Z2≦MaxZの条件を満たすことが上記(2)式の前提となる。上記第１例において、範囲パラメータΔZmaxは任意の値を取り得るようにし、これに応じて係数β１を調整した上で(2)式を用いてΔZ1を決定するようにしてもよい。
【００９３】
上記第1例によると、深さパラメータZsがその初期値又は標準値からずれればずれるほど、範囲パラメータΔZは小さくなる。つまり、混色変化が生じる深さ範囲が小さくなり、同時に、その範囲内での変化率が大きくなる。その変化の様子が図２７に示されている。深さパラメータZsが深さゼロ又は深さMaxZに近付くに従って、符号２２６Ａ，２２６Ｂ，２２６Ｃで示すように、混色範囲での変化率が急峻となっている。
【００９４】
図２８には、範囲パラメータ決定方法の第２例が示されている。横軸は深さＺを示しており、縦軸は範囲ΔZを示している。関数２２８に示すように、深さパラメータZsが小さくなるに従って、範囲パラメータΔZが増大されている。その変化率を定数β２とすると、以下の(３)式によって、範囲パラメータΔZが決定される（符号２３０参照）。
【００９５】
ΔZ=β2（MaxZ−Zs）
【００９６】
上記の範囲パラメータΔZの変化に伴う混色関数の変化が図２９に示されている。横軸が深さＺを示しており、縦軸がフロント重みを示している。符号２３０Ａ、２３０Ｂ、２３０Ｃで示すように、深さパラメータZsの増大に伴い、混色関数が左側へシフトしつつ、混色範囲内での変化率が緩和されている。
【００９７】
以上示したパラメータ決定方法は例示に過ぎず、上記以外の手法を利用して、深さパラメータZs及び範囲パラメータΔZが決定されてもよい、また、それらのパラメータに代えて他のパラメータ（例えば深さZ1,Z2）が自動的に決定されてもよい。混色関数として曲線からなる関数（例えばロジステック曲線）を利用してもよい。いずれにしても、混色関数を規定する複数のパラメータの内で少なくとも１つのパラメータを自動的に決定すればユーザーの負担を軽減できる。
【００９８】
パラメータ自動調整を行うタイミングとして、幾つかのタイミングをあげることができる。例えば、混色関数の調整を行うスイッチがオンされた時点の直後の画像更新タイミングでパラメータ自動調整が実行されてもよい。その後において、例えば、混色関数に関わる何等かの操作が行われたタイミング、画像が更新されたタイミング、一定の時間間隔をおいたタイミング、で自動的にパラメータ自動調整が実行されてもよい。
【００９９】
図３０には、カラー処理部３２の変形例が示されている。カラー処理部３２は、フロント側変換テーブル２３２、バック側変換テーブル２３４及び重み付け合成器２３６を有している。重み付け合成器２３６には、上記同様に、深さパラメータZs及び範囲パラメータΔZ等のパラメータ群が与えられる。フロント側変換テーブル２３２及びバック側変換テーブル２３４はそれぞれルックアップテーブルとして構成される。フロント側変換テーブル２３２は、入力される輝度Ｉ及び深さＺに基づいてメインカラーを決定するテーブルである。バック側変換テーブル２３４は、入力される輝度Ｉ及び深さＺに基づいてサブカラーを決定するテーブルである。それらのカラーが重み付け合成器２３６において合成されて画素のカラーが決定される。カラー処理部３２において、３つ以上のテーブルによって決定される３つのカラーを合成するようにしてもよい。
【符号の説明】
【０１００】
２０  ボリュームレンダリング部、２６  ガイド像生成部、２８  パラメータ決定部、３０  表示処理部、３２  カラー処理部、４６  輝度マップ、４８  深さマップ（距離マップ）。
【図１】
2018000780.tif 000003
【図２】
2018000780.tif 000004
【図３】
2018000780.tif 000005
【図４】
2018000780.tif 000006
【図５】
2018000780.tif 000007
【図６】
2018000780.tif 000008
【図７】
2018000780.tif 000009
【図８】
2018000780.tif 000010
【図９】
2018000780.tif 000011
【図１０】
2018000780.tif 000012
【図１１】
2018000780.tif 000013
【図１２】
2018000780.tif 000014
【図１３】
2018000780.tif 000015
【図１４】
2018000780.tif 000016
【図１５】
2018000780.tif 000017
【図１６】
2018000780.tif 000018
【図１７】
2018000780.tif 000019
【図１８】
2018000780.tif 000020
【図１９】
2018000780.tif 000021
【図２０】
2018000780.tif 000022
【図２１】
2018000780.tif 000023
【図２２】
2018000780.tif 000024
【図２３】
2018000780.tif 000025
【図２４】
2018000780.tif 000026
【図２５】
2018000780.tif 000027
【図２６】
2018000780.tif 000028
【図２７】
2018000780.tif 000029
【図２８】
2018000780.tif 000030
【図２９】
2018000780.tif 000031
【図３０】
2018000780.tif 000032