(57)【特許請求の範囲】
【請求項１】
生体内の三次元空間から取得された超音波ボリュームデータに対して、レンダリング処理が適用される三次元関心領域を設定する三次元関心領域設定手段と、
前記三次元関心領域内のデータを用いてレンダリング処理を実行し、これにより三次元超音波画像を形成する三次元超音波画像形成手段と、を含み、
前記三次元関心領域は、画像化対象組織と画像化非対象組織とを空間的に分離するためのクリッピング面を有し、
前記三次元関心領域設定手段は、
前記クリッピング面を変形させる変形手段と、
前記クリッピング面それ全体を傾斜させる傾斜手段と、
前記クリッピング面の傾斜角度に応じて前記クリッピング面のサイズを調整するサイズ調整手段と、
を有することを特徴とする超音波ボリュームデータ処理装置。
【請求項２】
請求項１記載の装置において、
前記変形手段は、前記クリッピング面を凸形状又は凹形状に変形させる、ことを特徴とする超音波ボリュームデータ処理装置。
【請求項３】
生体内の三次元空間から取得された超音波ボリュームデータに対して、レンダリング処理が適用される三次元関心領域を設定する三次元関心領域設定手段と、
前記三次元関心領域内のデータを用いてレンダリング処理を実行し、これにより三次元超音波画像を形成する三次元超音波画像形成手段と、を含み、
前記三次元関心領域は、画像化対象組織と画像化非対象組織とを空間的に分離するためのクリッピング面を有し、
前記三次元関心領域設定手段は、
前記クリッピング面を変形させる変形手段と、
前記クリッピング面それ全体を傾斜させる傾斜手段と、
を有し、
前記クリッピング面は、傾斜運動する３つの方向としての第１水平方向、第２水平方向及び垂直方向を有し、
前記変形手段は、ユーザー指定されたパラメータｈに従って前記クリッピング面の代表点についての前記垂直方向の高さを定める、
ことを特徴とする超音波ボリュームデータ処理装置。
【請求項４】
請求項３記載の装置において、
前記第１水平方向において前記代表点の両側の形状が線対称の関係にあり、
前記第２水平方向において前記代表点の両側の形状が線対称の関係にある、
ことを特徴とする超音波ボリュームデータ処理装置。
【請求項５】
請求項１記載の装置において、
前記傾斜手段は、第１傾斜角度θ１及び第２傾斜角度θ２に従って前記クリッピング面の二次元傾斜姿勢を定める、
ことを特徴とする超音波ボリュームデータ処理装置。
【請求項６】
請求項１記載の装置において、
前記クリッピング面はレンダリング開始面である、ことを特徴とする超音波ボリュームデータ処理装置。
【請求項７】
請求項１記載の装置において、
前記サイズ調整手段は、前記クリッピング面の傾斜角度が大きくなればなるほど、前記クリッピング面のサイズを大きくする、ことを特徴とする超音波ボリュームデータ処理装置。
【請求項８】
請求項７記載の装置において、
第１過程において前記変形手段によって変形後のクリッピング面が生成され、次の第２過程において前記傾斜手段によって前記変形後のクリッピング面を傾斜させることにより変形後かつ傾斜後のクリッピング面が生成され、更に第３過程において前記サイズ調整手段によって前記変形後かつ傾斜後のクリッピング面のサイズを調整することにより変形後かつ傾斜後かつサイズ調整後のクリッピング面が生成される、ことを特徴とする超音波ボリュームデータ処理装置。
【請求項９】
請求項８記載の装置において、
前記変形手段は、前記変形後のクリッピング面として第１の三次元形状データを生成し、
前記傾斜手段は、前記第１の三次元形状データの回転変換により前記変形後かつ傾斜後のクリッピング面として第２の三次元形状データを生成し、
前記サイズ調整手段は、前記第２の三次元形状データの拡大変換により前記変形後かつ傾斜後かつサイズ調整後のクリッピング面として第３の三次元形状データを生成し、
前記第３の三次元形状データに基づいて前記レンダリングを実行する際の各レイ上のボクセル演算開始点が規定される、ことを特徴とする超音波ボリュームデータ処理装置。
【請求項１０】
請求項１記載の装置において、
前記超音波ボリュームデータに基づいて互いに直交する第１及び第２の断層画像を形成する断層画像形成手段と、
前記三次元関心領域についての互いに直交する２つの断面を表す第１及び第２のグラフィック画像を形成するグラフィック画像形成手段と、
前記第１の断層画像上に前記第１のグラフィック画像が合成された第１の表示画像を表示し、前記第２の断層画像上に前記第２のグラフィック画像が合成された第２の表示画像を表示する表示手段と、
を含み、
前記三次元関心領域が変更されると、それに連動して、前記第１及び第２のグラフィック画像の内容が変更される、ことを特徴とする超音波ボリュームデータ処理装置。
【請求項１１】
請求項１記載の装置において、
複数の初期パラメータセットを格納した記憶部と、
前記複数の初期パラメータセットの中から特定のパラメータセットを選択する選択手段と、
を含み、
前記三次元関心領域設定手段は、前記特定のパラメータセットに従って初期三次元関心領域を設定する、ことを特徴とする超音波ボリュームデータ処理装置。
【請求項１２】
生体内の三次元空間から取得された超音波ボリュームデータに対してレンダリングが適用される三次元関心領域を設定する機能を備え、超音波ボリュームデータ処理装置において実行されるプログラムであって、
当該プログラムは、前記三次元関心領域に含まれる面であって画像化対象組織と画像化非対象組織とを空間的に分離するためのクリッピング面を生成するクリッピング面生成機能を有し、
前記クリッピング面生成機能は、ユーザー入力に基づいて前記クリッピング面を変形させる変形機能と、ユーザー入力に基づいて前記クリッピング面それ全体を傾斜させる傾斜機能と、
前記クリッピング面の傾斜角度に応じて前記クリッピング面のサイズを調整するサイズ調整機能と、
を有することを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は超音波ボリュームデータ処理装置に関し、特に三次元画像処理が適用される三次元関心領域を設定する技術に関する。
【背景技術】
【０００２】
近時、医療分野において、三次元超音波診断が普及しつつある。例えば、産科においては、母体内の胎児を包含する三次元空間（送受波空間）に対し、超音波が送受波され、これにより、ボリュームデータ（超音波ボリュームデータ）が取得される。そのボリュームデータに対してレンダリング処理を適用することにより、胎児の三次元画像が形成される。三次元画像形成法としては、ボリュームレンダリング法等が周知である。そのようなレンダリング法では、画像化の対象となる空間を通過する複数のレイ（視線）が設定され、各レイごとに開始点から終了点まで逐次的にボクセル演算が実行される。そして、各レイについての最終的なボクセル演算結果が画素値に対応付けられ、複数の画素値の集合として三次元画像が形成される。他の画像処理法として、サーフェイスレンダリング法、積算投影法、等が知られている。
【０００３】
対象組織についての三次元画像を形成する場合には、対象組織に隣接している他の組織（非対象組織）の画像化をできるだけ防止することが望まれる。特にボリュームレンダリング法において各レイ上のボクセル演算が対象組織の手前側に存在する他の組織に対しても適用されてしまうと、最終的に形成される三次元画像上において、対象組織が他の組織の裏側に隠れてしまい、対象組織を観察できなくなるという問題が生じる。例えば、子宮内の胎児を表した三次元画像において、最も観察したい胎児の顔が、その手前にある胎盤に隠れてしまう、という問題がある。
【０００４】
そこで、三次元関心領域（３Ｄ−ＲＯＩ）が利用されている。三次元関心領域は、三次元のデータ処理空間内に存在する部分的な空間であって、レンダリング処理が適用される範囲を限定するための空間である。三次元関心領域内に胎児（正確には胎児データ）が位置し、且つ、胎盤（正確には胎盤データ）が三次元関心領域外に位置するように、三次元データ処理空間内に三次元関心領域を設定すれば、上記問題を解消することができる。三次元データ処理空間において、三次元関心領域は立体的な図形として観念される。通常、三次元関心領域はレンダリング開始面を有する。例えば、胎児の三次元画像を形成する場合、レンダリング開始面ができる限り胎児の顔と胎盤との間に位置するように、当該レンダリング開始面がマニュアル又は自動的に設定される。レンダリング開始面は、組織分離機能を有するので、それをクリッピング面と称することができる。
【０００５】
特許文献１には、三次元関心領域の形状を変形させる技術が開示されている。この技術によれば、三次元関心領域として、立方体形状の特定の角部分を斜めに切り落とした形状を生成できる。特許文献２には、複数の点に基づいて定義された自由曲面を有する三次元関心領域が開示されている。同じく特許文献３には、自由曲面を有する三次元関心領域が開示されている。特許文献３に開示された技術では、自由曲面の形状が、直交する２つの断面上におけるユーザーの座標指定によって定義されている。このため、一方の断面に対する座標指定を他方の断面に反映させるために仮想点という特別な点が用いられている。特許文献４には三次元関心領域を自動的に設定する技術が開示されている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００４−１１３６０３号公報（図４）
【特許文献２】特開２００４−３３６５８号公報（図２）
【特許文献３】特開２００６−６１６９８号公報（図６Ｃ）
【特許文献４】特開２００１−１４５６３１号公報（第００４４段落）
【発明の概要】
【発明が解決しようとする課題】
【０００７】
生体内において、画像化したい組織の形状や画像化したくない組織の形状は色々であるため、三次元関心領域の形状として単なる立方体形状を採用するのは困難である。特に、レンダリング開始面であるクリッピング面を単なる平面とすると、画像化したい組織と画像化したくない組織とを分離することが困難となる。よって、特許文献１に記載された技術は採用し難い。
【０００８】
その一方、ユーザー負担の軽減及び検査時間の短縮のためには、三次元関心領域を簡便かつ迅速に設定できることが望まれる。特に、レンダリング開始面あるいはクリッピング面の位置、形状及び姿勢を簡便かつ迅速に設定できることが望まれる。その観点から見て、特許文献２に記載された技術は、断層画像上において、複数の点を組織表面形状に沿って個別的に位置決める必要があるので、操作上の負担が生じやすい。
【０００９】
本発明者らの研究によれば、画像化したい組織の表面形状及び画像化したくない組織の表面形状は、多くの場合に、およそ凸面状及び凹面状のいずれかである、ということが判明している。また、それらの間の隙間は、多くの場合に、三次元関心領域の中心軸（レンダリング中心線）に対して傾斜している、ということが判明している。このような生体組織の傾向に鑑みると、特許文献３の技術は採用しがたい。同技術では、三次元関心領域の上面（クリッピング面）の形状を自在に変化させることができるものの、それ全体の傾きを変化させることはできないからである。当該上面の四隅の高さは常に互いに同一であり、傾斜した隙間に対して上面を挿入することが非常に困難であると言い得る。もちろん、三次元関心領域それ自体をボリュームデータに対して相対的に傾斜させることもできる。しかし、その場合、クリッピング面をそれ全体として傾斜させることはできるものの、レンダリング方向が変化してしまう、又は、三次元関心領域内に画像化不要な組織が入り込んでしまう、といった他の問題が生じてしまう。
【００１０】
本発明の目的は、レンダリングが適用される三次元関心領域におけるクリッピング面の姿勢を適切に設定できるようにすることにある。あるいは、レンダリングが適用される三次元関心領域におけるクリッピング面の姿勢及び形状をそれぞれ独立して速やかに設定できるようにし、これによりユーザー操作時の負担を軽減し、また検査時間を短縮できるようにすることにある。
【課題を解決するための手段】
【００１１】
後述の実施形態に係る超音波ボリュームデータ処理は、生体内の三次元空間から取得された超音波ボリュームデータに対してレンダリング処理が適用される三次元関心領域を設定する三次元関心領域設定手段と、前記三次元関心領域内のデータを用いてレンダリング処理を実行し、これにより三次元超音波画像を形成する三次元超音波画像形成手段と、を含み、前記三次元関心領域は、画像化対象組織と画像化非対象組織とを空間的に分離するためのクリッピング面を有し、前記三次元関心領域設定手段は、前記クリッピング面を変形させる変形手段と、前記クリッピング面それ全体を傾斜させる傾斜手段と、を有することを特徴とするものである。
【００１２】
上記構成において、超音波ボリュームデータは生体内の三次元空間に対する超音波の送受波によって取得される。三次元関心領域設定手段は超音波ボリュームデータに対して仮想的な三次元関心領域を設定する。当該三次元関心領域はレンダリング処理が適用される領域である。レンダリング処理は、ボリュームレンダリング法に基づく処理であるのが望ましいが、他の処理であってもよい。三次元関心領域は分離面あるいは境界面として機能するクリッピング面を有する。クリッピング面は、特に望ましくは、レンダリング開始面であるが、他の面であってもよい。いずれにしても、クリッピング面は、画像化対象組織と画像化非対象組織とを空間的に分離することを目的とした面である。三次元関心領域設定手段は変形手段と傾斜手段とを備える。それらは望ましくはソフトウエアの機能として実現される。変形手段はクリッピング面の形状を変化させるものである。その場合にクリッピング面が通過する１又は複数の代表点を動かすことにより、クリッピング面を変形させるようにしてもよい。望ましくは、クリッピング面の形状は１つの代表点の高さ位置によって定められる。傾斜手段は、クリッピング面の変形後において、あるいは、変形前において、クリッピング面それ全体を傾斜させる。第１方向にだけ傾斜可能な構成を採用するようにしてもよいが、望ましくは、第１方向及び第２方向の両方向に傾斜可能な構成が採用される。クリッピング面それ全体を傾斜させることができるので、２つの組織表面が対向している状況において、それらの隙間の傾斜方向に合わせてクリッピング面を傾けることが容易となる。クリッピング面が例えば４つの角を有する場合、クリッピング面の二次元傾斜によれば、４つの角の高さを任意に設定できる。予めクリッピング面を変形させておいてからそれを傾斜させるようにしてもよいし、クリッピング面を傾斜させた上でそれを変形させるようにしてもよいし、変形と傾斜とを同時又は交互に行えるようにしてもよい。クリッピング面の傾斜により、そのサイズの補正が必要となる場合には、クリッピング面に対して拡大処理（又は縮小処理）が適用される。
【００１３】
上記の構成によれば、生体内における画像化対象組織と画像化非対象組織との間にクリッピング面を容易かつ迅速に設定できるという利点を得られる。特に、クリッピング面の形状を適切なものにでき、同時に、その姿勢を適切なものにできるから、レンダリング処理結果の品質を向上できる。例えば、胎児の三次元画像を形成する場合に、胎児の顔が胎盤に隠れてしまうといった問題を容易に回避することが可能である。超音波ボリュームデータ処理装置は、超音波診断装置によって構成されてもよいし、超音波診断装置で取得されたデータを処理するコンピュータによって構成されてもよいし、他の装置によって構成されてもよい。
【００１４】
望ましくは、前記変形手段は、前記クリッピング面を凸形状又は凹形状に変形させる。組織間の隙間の形状は、経験則上、概ね凸面又は凹面であることから、クリッピング面の形状として凸形状及び凹形状の両方を選択できるように構成するのが望ましい。クリッピング面の形状として凹凸面のような複雑な形状を採用することも勿論可能である。但し、ユーザーの操作性を考慮するならば、クリッピング面の中心にある代表点を上下に移動させて、それを通過する単純な凸面又は凹面が設定されるように構成するのが望ましい。
【００１５】
望ましくは、前記クリッピング面は、傾斜運動する３つの方向としての第１水平方向、第２水平方向及び垂直方向を有し、前記変形手段は、ユーザー指定されたパラメータｈに従って前記クリッピング面の代表点についての前記垂直方向の高さを定める。３つの方向は直交関係を維持しながら運動するものである。垂直方向は法線方向に相当し、当該方向における代表点の高さの設定により、クリッピング面の変形方向及び曲率を自在に設定できる。
【００１６】
望ましくは、前記第１水平方向において前記代表点の両側の形状が線対称の関係にあり、前記第２水平方向において前記代表点の両側の形状が線対称の関係にある。このような対称性によればユーザーをして三次元関心領域を直感的に容易に認識させることができる。
【００１７】
望ましくは、前記傾斜手段は、第１傾斜角度θ１及び第２傾斜角度θ２に従って前記クリッピング面の二次元傾斜姿勢を定める。２つの方向について独立して傾斜角度を設定できるようにすれば、組織間の隙間が存在する方向にクリッピング面を合わせることが容易となる。
【００１８】
望ましくは、前記クリッピング面の傾斜角度に応じて前記クリッピング面のサイズを調整するサイズ調整手段を有する。傾斜前の状態においてクリッピング面の形状が定義された場合、そこから傾斜を行うと、クリッピング面が三次元関心領域の本体図形（クリッピング面を除いた部分）に合致せずに、両者間に空間的なずれが生じてしまう。そこで、サイズ調整（特に拡大処理）を行えば、両者をサイズ面において整合させることができる。これにより正規の立体形状を維持できる。なお、サイズ調整は、クリッピング面又は本体図形のいずれにも適用可能である。拡大補正の他、縮小補正を採用することもできる。
【００１９】
望ましくは、前記サイズ調整手段は、前記クリッピング面の傾斜角度が大きくなればなるほど、前記クリッピング面のサイズを大きくする。この場合、相似形を維持してクリッピング面のサイズが変化するように構成するのが望ましい。
【００２０】
望ましくは、第１過程において前記変形手段によって変形後のクリッピング面が生成され、次の第２過程において前記傾斜手段によって前記変形後のクリッピング面を傾斜させることにより変形後かつ傾斜後のクリッピング面が生成され、更に第３過程において前記サイズ調整手段によって前記変形後かつ傾斜後のクリッピング面のサイズを調整することにより変形後かつ傾斜後かつサイズ調整後のクリッピング面が生成される。この構成によれば、演算を簡略化できるから、ソフトウエア処理を迅速に行える。望ましくは、前記変形手段は、前記変形後のクリッピング面として第１の三次元形状データを生成し、前記傾斜手段は、前記第１の三次元形状データの回転変換により前記変形後かつ傾斜後のクリッピング面として第２の三次元形状データを生成し、前記サイズ調整手段は、前記第２の三次元形状データの拡大変換により前記変形後かつ傾斜後かつサイズ調整後のクリッピング面として第３の三次元形状データを生成し、前記第３の三次元形状データに基づいて前記レンダリングを実行する際の各レイ上のボクセル演算開始点が規定される。クリッピング面を構成する点群の位置形状データを用いずにクリッピング面を関数で表現することもできる。
【００２１】
望ましくは、前記超音波ボリュームデータに基づいて互いに直交する第１及び第２の断層画像を形成する断層画像形成手段と、前記三次元関心領域についての互いに直交する２つの断面を表す第１及び第２のグラフィック画像を形成するグラフィック画像形成手段と、前記第１の断層画像上に前記第１のグラフィック画像が合成された第１の表示画像を表示し、前記第２の断層画像上に前記第２のグラフィック画像が合成された第２の表示画像を表示する表示手段と、を含み、前記三次元関心領域が変更されると、それに連動して、前記第１及び第２のグラフィック画像の内容が変更される。ユーザーはそのような第１及び第２の表示画像を通じて三次元関心領域の位置、形状等を認識できる。
【００２２】
望ましくは、複数の初期パラメータセットを格納した記憶部と、前記複数の初期パラメータセットの中から特定のパラメータセットを選択する選択手段と、を含み、前記三次元関心領域設定手段は、前記特定のパラメータセットに従って初期三次元関心領域を設定する。複数の初期形状に対応する複数の初期パラメータセットを予め用意しておけば、所望の三次元関心領域を迅速に設定でき、その場合におけるユーザーの負担を少なくできる。過去に設定された三次元関心領域を規定するパラメータセットを追加登録して後日利用できるように構成するのが望ましい。
【００２３】
後述する実施形態に係るプログラムは、生体内の三次元空間から取得された超音波ボリュームデータに対してレンダリングが適用される三次元関心領域を設定する機能を備え、超音波ボリュームデータ処理装置において実行されるプログラムであって、当該プログラムは、前記三次元関心領域に含まれる面であって画像化対象組織と画像化非対象組織とを空間的に分離するためのクリッピング面を生成するクリッピング面生成機能を有し、前記クリッピング面生成機能は、ユーザー入力に基づいて前記クリッピング面を変形させる変形機能と、ユーザー入力に基づいて前記クリッピング面それ全体を傾斜させる傾斜機能と、を有するものである。
【発明の効果】
【００２４】
本発明によれば、レンダリングが適用される三次元関心領域におけるクリッピング面の姿勢を適切に設定できる。あるいは、レンダリングが適用される三次元関心領域におけるクリッピング面の姿勢及び形状をそれぞれ独立して速やかに設定できるようにして、ユーザー操作時の負担を軽減し、また検査時間を短縮できる。
【図面の簡単な説明】
【００２５】
【図１】本発明に係る超音波ボリュームデータ処理装置としての超音波診断装置の全体構成を示すブロック図である。
【図２】表示器に表示される画像の一例を示す図である。
【図３】三次元関心領域を利用したレンダリング処理を説明するための図である。
【図４】記憶部に格納される複数のパラメータセットを示す図である。
【図５】制御部の機能として三次元関心領域の生成処理を説明するための図である。
【図６】傾斜していない変形後のクリップ面を表す概念図である。
【図７】傾斜後のクリップ面を示す概念図である。
【図８】拡大処理後のクリップ面を示す概念図である。
【図９】三次元関心領域の設定方法を説明するためのフローチャートである。
【図１０】三次元関心領域のＹＸ断面形状およびＹＺ断面形状を示す図である。
【図１１】三次元関心領域のＹＸ断面形状の生成処理を示すフローチャートである。
【図１２】三次元関心領域のＹＺ断面形状の生成処理を示すフローチャートである。
【発明を実施するための形態】
【００２６】
図１には、本発明に係る超音波ボリュームデータ処理装置としての超音波診断装置がブロック図として示されている。この超音波診断装置は医療の分野において用いられ、超音波の送受波により生体内の組織の三次元画像を形成する機能を有している。本実施形態において、画像化の対象となる組織は胎児である。もちろん他の組織を画像化するようにしてもよい。
【００２７】
プローブ１０は三次元エコーデータを取り込むための送受波手段であり、本実施形態において、プローブ１０は１Ｄアレイ振動子とそれを機械的に走査する走査機構とを内蔵している。１Ｄアレイ振動子による超音波ビームの電子走査により走査面が形成され、その走査面を機械的に走査することにより、三次元エコーデータ取り込み空間としての三次元空間１２が形成される。プローブ１０に２Ｄアレイ振動子を設け、超音波ビームを二次元的に電子走査するようにしてもよい。このような方式によっても三次元空間１２を形成することが可能である。電子走査方式としては、電子セクタ走査、電子リニア走査等が知られている。プローブ１０は体表面上に当接されるものであるが、もちろんプローブ１０として体腔内に挿入されるものが用いられてもよい。胎児の超音波診断を行う場合には、プローブ１０が母体の腹部表面上に当接され、その状態において超音波の送受波が実行される。
【００２８】
送受信部１６は送信ビームフォーマおよび受信ビームフォーマとして機能するものである。すなわち、送信時において、送受信部１６からプローブ１０へ複数の送信信号が供給され、これにより送信ビームが形成される。受信時において、生体内からの反射波はプローブ１０において受波され、これによりプローブ１０から送受信部１６へ複数の受信信号が出力される。送受信部１６において、複数の受信信号に対する整相加算処理が実行され、これにより整相加算後の受信信号としてビームデータが出力される。ビームデータに対しては図示されていない信号処理モジュールによる各種の信号処理が適用され、そのような信号処理後のビームデータが３Ｄメモリ１８へ格納される。
【００２９】
３Ｄメモリ１８は、三次元空間１２に対応するデータ空間を有し、３Ｄメモリ１８内には三次元空間１２から取得されたエコーデータ集合体としてのボリュームデータが格納される。ボリュームデータは実際には複数本のビームデータの座標変換及び補間処理により構成されるものである。ドプラ情報から成るボリュームデータが構成されるようにしてもよい。各エコーデータに対する座標変換がデータ読み出し時に実行されてもよい。
【００３０】
ボリュームレンダリング部２０は、制御部３６から与えられたレンダリング条件に従って、三次元関心領域内のデータを利用してレンダリング処理を実行し、これにより三次元画像を形成する。その画像データは表示処理部２２へ出力される。ボリュームレンダリング法としては各種の手法が知られており、様々な方法を採用することができる。なお、ボリュームレンダリング法に代えてサーフェイスレンダリング法等の他の画像処理法を適用することも可能である。
【００３１】
断層画像形成部２４，２６，２８は、それぞれＢモード白黒断層画像を形成するモジュールである。本実施形態においては、互いに直交関係にある切断面に対応した３つの断層画像が形成されている。それらの断層画像は一般にトリプレーンと称されている。本実施形態においては、三次元関心領域における中心点（原点）を通過する３つの切断面に対応する３つの断層画像が形成されている。勿論、各切断面を後述するＸ方向、Ｙ方向あるいはＺ方向に変化させるようにしてもよい。断層画像形成部２４，２６，２８により生成された３つの断層画像のデータは表示処理部２２へ出力される。なお、断層画像形成部２４，２６，２８の画像処理にあたって必要な条件は制御部３６から渡されている。
【００３２】
図形イメージ生成部３０，３２，３４は、３つの断層画像に対してオーバーレイ表示されるグラフィックイメージを生成するモジュールである。本実施形態においては、図形イメージ生成部３０が関心領域のＹＸ断面を表すグラフィックイメージを生成しており、図形イメージ生成部３２が三次元関心領域のＹＺ断面を表すグラフィックイメージを生成している。図形イメージ生成部３４は、三次元関心領域のもう一つの断面形状を表すグラフィックイメージを生成している。このように生成された各イメージは表示処理部２２に出力されている。図形イメージ生成部３０，３２，３４におけるイメージ生成にあたって必要な条件およびデータは制御部３６から渡されている。
【００３３】
表示処理部２２は画像合成機能を有し、入力される複数の画像データに基づいて表示画像を構成し、それを表すデータを表示器４２へ出力する。表示器４２に表示される画像例については後に図２を用いて説明する。
【００３４】
制御部３６は、本実施形態において、ＣＰＵおよび動作プログラムによって構成されている。ちなみにボリュームレンダリング部２０、断層画像形成部２４，２６，２８、図形イメージ生成部３０，３２，３４もソフトエアの機能として実現することが可能である。制御部３６には記憶部３８が接続されており、また制御部３６には入力部４０が接続されている。入力部４０は本実施形態において操作パネルにより構成され、その操作パネルはキーボードやトラックボール等を有するものである。ユーザーは入力部４０を用いて三次元関心領域の設定にあたって必要な数値を入力することが可能である。記憶部３８には、後に図４を用いて説明する複数のパラメータセットが予め格納されている。また記憶部３８にはワークエリアが確保されており、そのワークエリア上にはクリップ面を表す形状データが必要に応じて格納される。制御部３６は、本実施形態において三次元関心領域の設定機能を有している。その処理内容については後に図５乃至図１２を用いて詳述する。
【００３５】
図２には、表示器に表示される画像の一例が示されている。符号４４は表示画像を表している。その表示画像４４には、３つの断層画像４６，４８，５０が含まれ、更に三次元画像５２が含まれている。３つの断層画像４６，４８，５０は互いに直交関係を有し、それら全体としてトリプレーンを構成している。具体的に説明すると、断層画像４６はＹＸ断面を表す断層画像であり、その断層画像４６内にはグラフィックイメージとしての断面形状ＳYXが含まれている。この断面形状ＳYXは、三次元関心領域のＹＸ断面を表すものであり、それにはクリッピング面の断面を表すカーブ１１８が含まれている。符号５４は画像化対象となる組織を表しており、符号５６は画像化をしたくない組織を表している。カーブ１１８は組織５４と組織５６の隙間に沿って湾曲かつ傾斜している。
【００３６】
断層画像４８には、三次元関心領域のＹＺ断面を表す断面形状ＳYZが含まれている。それはクリッピング面のＹＺ断面を表すカーブ１２８を有している。カーブ１２８は画像化対象となる組織５４とそうではない組織５６との間の隙間に沿って傾斜および湾曲しており、２つの組織の分離が行われている。図２においては、凸形状を有するカーブ１１８，１２８が表示されていたが、凹形状を有するカーブ１１８，１２８が表示されることもある。
【００３７】
断層画像５０には三次元関心領域のＸＺ断面を表す四角形の断面形状５８が表れている。もちろん、三次元関心領域の位置やサイズはユーザにより任意に可変することが可能である。また後に詳述するように、クリッピング面の形状および姿勢を任意に設定することが可能である。
【００３８】
三次元画像５２は、以上のように設定される三次元関心領域内のデータを用いてレンダリング処理を実行することにより生成された画像である。本実施形態において、クリッピング面がレンダリング開始面に相当しており、図示の例においてはＹ方向に沿って各レイが設定されており、投影視点はＹ方向上方にある。なお、レンダリング処理にあたっては、各レイ上においてボクセルデータがサンプリングされるが、その場合においては各ボクセルデータは周辺に存在する複数のエコーデータを参照することにより補間処理によって生成される。補間処理を行う場合において三次元関心領域外のデータが参照されることもある。
【００３９】
図３には、三次元関心領域を用いたレンダリング処理が概念図として示されている。三次元関心領域Ｖに対してはＹ方向に平行に複数のレイ６４が設定されている。もちろん、平行投影法ではなく他の投影法を適用してもよい。三次元関心領域Ｖにおけるクリッピング面６０は本実施形態において変形可能且つ傾斜可能なものであり、クリッピング面６０は上述したようにレンダリング開始面に相当する。クリッピング面６０と反対側の面が終了面６２であり、それは図３において底面のように表されている。各レイ上においては、所定のサンプリング間隔をもってボクセル演算が逐次的に実行される。各ボクセルごとに入力光量に基づいて出力光量が演算され、その場合においてはオパシティ（不透明度）がパラメータとして利用される。ボクセル演算が各レイに沿って順次進行し、その処理が終了面６２に到達した段階で、当該レイ上のボクセル演算は終了する。あるいは各ボクセルごとに演算される出力光量が最大値（例えば１．０）に到達した場合にもボクセル演算は終了する。その終了時点の出力光量が画素値として用いられ、すなわち仮想的なスクリーン６８上において、レイ６４に対応する画素７０の画素値が当該レイ６４において求められた最終の出力光量により決定される。複数のレイについて同様の処理を行うならば、スクリーン６８上に三次元画像が構築されることになる。
【００４０】
図４には、図１に示した記憶部に格納される複数のパラメータセット７４の一例が示されている。ここで１つのパラメータセット７４に着目すると、当該パラメータセット７４は、三次元関心領域の原点Ｃの座標（Ｘ0，Ｙ0，Ｚ0）と、三次元関心領域のサイズ（ＸW，ＹW，ＺW）と、Ｙ方向におけるオフセット量（Ｙ+）と、クリッピング面の高さｈと、クリッピング面の傾斜角度θX，θZと、で構成されている。もちろん、この構成は一例であり、他の構成を採用することもできる。本実施形態においては、三次元関心領域として、例えば８つの初期形状が定められており、その８つの初期形状に対応した８つのパラメータセット７４が記憶部上に登録されている。よって、ユーザーはパラメータセットの選択により、所望の形状を有する三次元関心領域を瞬時に設定することが可能であり、その上で必要な修正を施すことが可能である。今まで生成した三次元関心領域を記憶部上に登録するようにしてもよい。
【００４１】
図５には、クリッピング面の生成に関する処理の流れが概念図として示されている。符号７６は制御部３６の機能を表しており、符号３８Ａは記憶部３８上のワークエリアを表している。また符号７４Ａはユーザにより選択されたパラメータセットを表している。そのパラメータセット７４Ａを構成する各パラメータを、符号９０で示すように、ユーザーによって任意に変更することが可能である。まず、第１段階においては、符号７８で示すように、スプライン処理により、傾斜していない状態で変形されたクリッピング面が生成される。その具体的な内容については後に図６等を用いて詳述する。その処理にあたってはパラメータセット７４Ａが参照される。スプライン処理７８により生成された傾斜前且つ変形後のクリッピング面は実際には三次元の形状データとして生成されており、それがワークエリア３８Ａ上に一時的に格納される。形状データは符号８０で表されている。
【００４２】
第２段階においては、回転変換処理８２が実行される。具体的にはパラメータセット７４Ａに基づいて、形状データ８０に対して回転変換処理８２を適用することにより、傾斜後のクリッピング面を表す形状データ８４が生成され、それがワークエリア３８Ａ上に格納される。この形状データ８４は過渡的なデータである。
【００４３】
第３段階においては、形状データ８４に対して拡大変換処理８６が適用される。この拡大変換処理８６により、拡大後のクリッピング面を表す形状データ８８が生成され、それがワークエリア３８Ａ上に格納される。この拡大変換処理８６により実際に機能させるクリッピング面が確定することになり、すなわち三次元関心領域の生成が完了する。そして、符号９０で示されるように、いずれかのパラメータがユーザーにより変更された場合、それをトリガとして、再び第１段階から第３段階までの処理が実行され、三次元関心領域が瞬時に更新されることになる。
【００４４】
次に、図６乃至図９を用いて、三次元関心領域の生成処理について、より具体的に説明する。
【００４５】
図６には、初期状態にある三次元関心領域Ｖ１が示されている。三次元関心領域Ｖ１はデータ処理空間上において仮想的に存在するものである。すなわち、以下に説明するクリッピング面１１２以外の本体については、レンダリング処理範囲の条件として数値上存在しているだけであり、実際にそのような形状が生成されているわけではない。ただし、本実施形態の説明に当たっては、その理解を助けるために、三次元関心領域を視覚的に認識できる図形であるものとみなす。
【００４６】
三次元関心領域Ｖ１は原点Ｃを有し、その原点Ｃの空間的位置は上述したパラメータセットに含まれる原点座標（Ｘ0，Ｙ0，Ｚ0）により定められる。三次元関心領域Ｖ１のＸ方向、Ｙ方向およびＺ方向のそれそれのサイズは上述したパラメータセット７４に含まれるサイズ情報（ＸW，ＹW，ＺW）により定められる。それらの情報は各方向における全体幅を定義するものであってもよいし、その半分を定義するものであってもよい。図６において、Ｏは基準点を表しており、それは中心線１００上の点である。基準点Ｏは後に説明するようにクリッピング面１１２を傾斜させる場合における回転中心（不動点）である。Ｙ方向のサイズとは別に、基準点のＹ方向の高さを調整したい場合にはオフセット値Ｙ+の値を可変すればよい。その値の可変により原点Ｃを維持しつつも、その上側のサイズを任意に定めることができる。原点Ｏに対して中心線１００に沿った距離として高さｈが定義されている。基準点Ｏから距離ｈだけ離れた点が代表点Ｐである。ｈが正の値であれば代表点Ｐは基準点Ｏよりも上方に位置し、それとは逆にｈが負の値をとる場合には代表点Ｐは基準点Ｏよりも下方に位置することになる。
【００４７】
ちなみに、図６において、三次元関心領域Ｖ１が有する８個の角がａ１，ａ２，ａ３，ａ４，ａ５，ａ６，ａ７，ａ８で表されている。そして、各辺の中点がａ１２，ａ２３，ａ３４，ａ１４，ａ５６，ａ６７，ａ７８，ａ５８で表されている。
【００４８】
クリッピング面１１２の生成について詳述すると、上述したように、パラメータｈにより代表点Ｐが定められると、代表点Ｐと２つの端点Ｐ１，Ｐ２とに基づくスプライン補間演算により、背骨のような基本線として曲線１０４が生成される。端点Ｐ１は上述した点ａ１２であり、端点Ｐ２は上述した点ａ３４である。曲線１０４は代表点Ｐを挟んでＺ方向両側において線対称の形状を有している。但し、図６においてはクリッピング面１１２の形状をより分かり易く表現方法が採用されている。
【００４９】
以上のように曲線１０４が定められると、辺Ｌ１と辺Ｌ３との間に複数のスプライン曲線１０６，１０８が順次生成される。具体的には、辺Ｌ２，辺Ｌ４に対して等間隔で複数の端点Ｐ３，Ｐ４，Ｐ５，Ｐ６が定められ、同様に曲線１０４上においてもＺ方向に等間隔で複数の通過点が定められる。そして、Ｚ方向の各位置において、２つの端点と１つの通過点とを用いてスプライン補間演算を行うことにより、スプライン曲線を生成することができる。図６においては、端点Ｐ３，Ｐ４を結ぶスプライン補間曲線が符号１０６で表されており、同じく端点Ｐ５，Ｐ６を結ぶスプライン補間曲線が符号１０８で示されている。そのような曲線をＺ方向の各位置において生成することにより、結果として曲線アレイ１１０を構成することができ、それによって傾斜前の状態におけるクリッピング面１１２が構成される。それは実際には三次元形状データからなるものである。ちなみに、図６に示す例においてはＺ方向に複数のスプライン保管曲線が並べられたが、それがＸ方向に生成されるようにしてもよい。
【００５０】
以上のような第１段階の処理により、傾斜していないクリッピング面１１２が生成されると、次に、図７に示すように、傾斜したクリッピング面１１４が生成される。具体的には、パラメータセットを構成する傾斜角度θX，θZが参照され、Ｘ方向およびＺ方向の両方向にクリッピング面が傾けられる。これを、法線１１５の傾斜として理解することもできるが、本実施形態においては法線は演算されていない。三次元の形状データの座標変換により、傾斜後のクリッピング面１１４が構成されている。
【００５１】
法線１１５は基準点Ｏと傾斜後の代表点Ｐ’とを結ぶ線であり、法線１１５のＸ方向の傾斜角度が上述したθXであり、法線１１５のＺ方向の傾斜角度が上述したθZである。クリッピング面１１４は、４つの角ｂ１，ｂ２，ｂ３，ｂ４を有している。図６に示したクリッピング面１１２を単純に傾斜させたため、図７に示す状態では、角ｂ１，ｂ２，ｂ３，ｂ４は、三次元形状本体よりも内側に位置している。このままでは三次元関心領域を構成できない。そこで、本実施形態においては、クリッピング面１１４の相似形を生成する拡大処理が実行されている。相似形が維持されるように拡大処理が適用される。具体的には、２つの水平方向ｕ，ｖ方向にクリッピング面１１４が拡大され、同時に法線１１５の方向（傾斜後の垂直方向）にも、クリッピング面１１４が拡大される。ただし、基準点の位置は不変である。拡大後のクリッピング面も三次元形状データとして構成されている。
【００５２】
図８には、拡大後のクリッピング面１１６が示されている。拡大処理の結果、４つの角ｃ１，ｃ２，ｃ３，ｃ４は４つの垂直辺上に位置している。これにより傾斜後のクリッピング面１１６を含む三次元関心領域Ｖ２が構成される。その三次元関心領域Ｖ２は８つの点ｃ１，ｃ２，ｃ３，ｃ４，ａ５，ａ６，ａ７，ａ８によって囲まれる形状を有している。一方、このような三次元関心領域Ｖ２を２つの断層画像上において確認するために断面形状ＳYX、ＳYZが生成される。その場合においては、２つの端点ｃ１２，ｃ３４を結ぶ断面上において上述同様のスプライン補間演算によりカーブが生成され、それと直交関係にある断面上においても２つの端点ｃ１４，ｃ２３の間におけるスプライン補間演算により上述したカーブが生成される。それらの断面形状の生成に関しては後に図１０乃至図１２を用いて詳述する。
【００５３】
図９には、以上説明した三次元関心領域の生成処理がフローチャートとして示されている。まず、Ｓ１０１では、複数の初期形状の中からユーザーにより実際に使用する初期形状が選択される。具体的には、ユーザーが選択された初期形状に対応するパラメータセットがプログラム実行上において認識されることになる。Ｓ１０２おいては、三次元関心領域の原点座標Ｃ（Ｘ0，Ｙ0，Ｚ0）、サイズ（ＸW，ＹW，ＺW）、オフセットＹ+、クリッピング面の高さｈ、に基づいて、傾斜前のクリッピング面が生成される。そのような生成にあたっては上述したスプライン補間演算が実行される。これによってクリッピング面を表す三次元形状データが生成される。Ｓ１０３においては、傾斜角度をθX，θZに基づいて、傾斜前のクリッピング面を空間的に変換することにより、傾斜後のクリップ面が生成される。そのクリッピング面も三次元形状データにより構成される。Ｓ１０４では、そのような傾斜後のクリッピング面に対して拡大処理が適用される。その場合においては、上述したように相似形を維持しながら形状を拡大する処理が実行される。拡大率については傾斜角度θX、θZから一義的に定めることができる。このような処理により実際に使用するクリッピング面が三次元形状データとして構成されることになる。Ｓ１０５においては以上のように構成されたクリッピング面を有する三次元関心領域が決定される。それは実際には三次元形状データとそこからのレンダリング範囲を規定する数値情報とにより構成されるものである。Ｓ１０６においては以上のように生成されたレンダリング条件に従ってボリュームレンダリング処理が実際に実行される。これにより三次元画像が構成され、その画像データが表示画面上に表示される。Ｓ１０７では、以上のような処理を続行するか否かが判断され、続行する場合にはＳ１０８においてユーザーによりいずれかのパラメータに対する変更があったか否かが判断され、変更有りと判断された場合には、Ｓ１０２以降の各工程が繰り返し実行される。
【００５４】
以上のように、本実施形態によれば、クリッピング面をシンプルな凸面あるいは凹面として構成することができ、しかも、そのようなクリッピング面を二次元傾斜させることができるので、図２に示したように、画像化対象となる組織とそうでない組織との間の隙間が二次元的に傾斜しているような場合においても、それに合わせてクリッピング面の姿勢を適切に設定できるという利点がある。本実施形態においては、クリッピング面の中心点を成す代表点の高さを定めるだけで、その形状を定義することができるから、ユーザーにおいて分かり易く混乱が生じることもない。そのようなシンプルな形状であっても、実際の生体組織の形状への適合性は良好であることが実験により確認されている。本実施形態においては、Ｘ方向およびＹ方向の両方向において、代表点の両側の形状が線対称の関係にあったが、代表点の両側において非対称の形状を採用することも可能である。ちなみに、立方体が有する６面のうちでいずれの面をクリッピング面にしてもよい。また、三次元関心領域をボリュームデータに対して相対的に任意に回転できるように構成するのが望ましい。
【００５５】
次に、図１０乃至図１２を用いて、図２に示した三次元関心領域の断面形状ＳYX，ＳYZの生成方法について説明する。
【００５６】
本実施形態においては、三次元関心領域の生成とその断面形状の生成は、互いに独立した演算処理により実現されている。（Ａ）には断面形状ＳYXの生成過程が示されており、（Ｂ）には断面形状ＳYZの生成過程が示されている。
【００５７】
（Ａ）において、みかけ上の代表点Ｐ１０は、中心線１１０上の高さｈに基づいて、まず奥行方向（Ｚ方向）に、傾斜していないカーブをスプライン補間により生成し、それを傾斜角度θｚ に基づいて回転及び拡大し、それにより生成されるカーブ上において中心線１１０に交わる点を特定することにより、当該交点として決定される。但し、これは一例に過ぎない。符号２１２が初期カーブを表している。この初期カーブ２１２はみかけ上の代表点Ｐ１０と２つの端点ａ１４,ａ２３とを結ぶスプライン補間演算により生成されるものである。補間方法としては各種のものがあるので、必要に応じて補間方法を切り替えてもよい。クリッピング面の生成過程と同様に、Ｘ方向の傾斜角度θXに基づいてカーブ２１２が回転され、その回転後のカーブが符号２１６で表されている。そこにおける２つの端点がａ１４’，ａ２３’であり、移動後の代表点がＰ１０Xで表されている。基準点Ｏと移動後の代表点Ｐ１０Xとを結ぶ直線が、みかけ上の法線２１４である。但し、代表点Ｐ１０Xや法線２１４は実際には演算されない。カーブ２１６に対して拡大処理を適用した後の形状がカーブ２１８である。これにより、２つの端点ｃ１４,ｃ２３が左右の垂直辺上に位置している。このような処理の結果、２つの点ｃ１４,ｃ２３,ａ６７,ａ５８に囲まれるボックス形状として三次元関心領域のＹＸ断面ＳYXが生成される。それは実際にはグラフィックイメージとして断層画面上にオーバーレイ表示される。
【００５８】
ＹＺ断面ＳYZについても上記同様に生成される。具体的には、まず３点Ｐ１２,ａ１２,ａ３４のスプライン補間演算によりカーブ２２２が生成され、Ｚ方向の傾斜角度θZに基づいてそのカーブ２１２が回転変換される。それが符号２２６により表されている。みかけ上の法線２２４上に代表点Ｐ１２Yが表されているが、それらは実際には計算されない。２つの端点ａ１２’,ａ３４’は左右の垂直ライン上から内側に外れている。そこで、それに対して拡大処理を適用した結果がカーブ２２８である。２つの端点ｃ１４,ｃ２３はいずれも左右の垂直ライン上に位置している。これによってＹＺ断面形状ＳYZが生成される。それは断層画像上にオーバーレイ表示される。
【００５９】
図１１には、図１０の（Ａ）に示した断面形状ＳYXの生成過程がフローチャートとして示されている。Ｓ２０１では、ユーザーの選択による初期データとしてパラメータセットが参照される。Ｓ２０２では、そのパラメータセットに基づいて、傾斜していないカーブが生成される。Ｓ２０３では、Ｘ方向の傾斜角度θXに基づいて、傾斜したカーブが生成される。Ｓ２０４では、傾斜角度θXに基づいて、傾斜後のカーブに対する拡大処理が適用される。Ｓ２０５では、拡大処理後のカーブを含むボックス形状として断面形状ＳYXが定められ、それが断層画像上にグラフィックとして表示される。Ｓ２０６では以上の処理を継続するか否かが判断され、Ｓ２０７においていずれかのパラメータについての変更があったと判断された場合、具体的には断面形状ＳYXに関係するパラメータについて変更があったと判断された場合に、Ｓ２０２以降の各工程が繰り返し実行される。
【００６０】
図１２には、図１０（Ｂ）に示した断面形状ＳYZの生成過程がフローチャートとして示されている。この内容は図１１に示した内容と基本的に同一であるが、念のため説明する。
【００６１】
Ｓ３０１では、初期データが参照され、Ｓ３０２ではパラメータセットが参照され、それに基づいて、傾斜していないカーブが生成される。そして、Ｓ３０３では、Ｚ方向の傾斜角度θZに基づいて、傾斜処理が実行されて、傾斜したカーブが生成される。そして、Ｓ３０４では、同じく傾斜角度θZに基づいて、拡大率が定められ、傾斜したカーブに対して拡大処理を適用した結果として拡大後のカーブが生成される。Ｓ３０５では、そのような拡大されたカーブを含む断面形状ＳYZが定められ、それが断層画像上に表示される。Ｓ３０６では、以上の処理を続行するか否かが判断され、続行する場合には、Ｓ３０７において、断面形状ＳYZに影響を与えるパラメータの値が変更されたか否かが判断され、変更があったと判断された場合には、Ｓ３０２以降の各工程が実行されることになる。
【００６２】
上記の実施形態においては、断面形状の生成にあたって２次元空間上での処理が中心となるので、その処理を迅速に行えるという利点が得られる。
【符号の説明】
【００６３】
４６，４８，５０  断層画像、５２  ３次元画像、１１２  傾斜前のクリッピング面、１１４  傾斜後のクリッピング面、１１６  拡大後のクリッピング面。
【図１】
0005437767.tif 000002
【図２】
0005437767.tif 000003
【図３】
0005437767.tif 000004
【図４】
0005437767.tif 000005
【図５】
0005437767.tif 000006
【図６】
0005437767.tif 000007
【図７】
0005437767.tif 000008
【図８】
0005437767.tif 000009
【図９】
0005437767.tif 000010
【図１０】
0005437767.tif 000011
【図１１】
0005437767.tif 000012
【図１２】
0005437767.tif 000013