超音波診断装置
【課題】血管と器具を両方とも同時に明瞭な三次元画像として可視化することのできる超音波診断装置を提供する。
【解決手段】器具が挿入された血管に超音波を送受波して得られたボクセルデータ群に対して設定された各レイについて、深さ方向にボクセル演算を逐次実行する。ボクセル演算を行うにあたって第１レンダリング手段によって、第１開始点から前記ボクセル演算の逐次実行結果により、第１三次元画像を構成する画素値を求める。前記各レイについて判定手段を用いて第２開始点を判定する。第２開始点が判定されたら第２レンダリング手段によって、前記器具が反映された第２三次元画像を構成する画素値を求める。画像を合成することによって器具が挿入された血管を三次元投影画像として認識することができる。
【解決手段】器具が挿入された血管に超音波を送受波して得られたボクセルデータ群に対して設定された各レイについて、深さ方向にボクセル演算を逐次実行する。ボクセル演算を行うにあたって第１レンダリング手段によって、第１開始点から前記ボクセル演算の逐次実行結果により、第１三次元画像を構成する画素値を求める。前記各レイについて判定手段を用いて第２開始点を判定する。第２開始点が判定されたら第２レンダリング手段によって、前記器具が反映された第２三次元画像を構成する画素値を求める。画像を合成することによって器具が挿入された血管を三次元投影画像として認識することができる。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は超音波診断装置に関し、特に、三次元超音波画像を形成する超音波診断装置に関する。
【背景技術】
【０００２】
血液透析を行う場合に、透析患者の体から血液を取り出して体内に戻すためには、前腕部の血管が使用されることが多い。血液透析を長く続けていると、前腕部の血管が痛んで狭窄が発生することがある。
【０００３】
狭窄が発生して内径が細くなった血管を治療するために、血管内にバルーンやステントを挿入し血管の内径を拡げる治療が行われている。この治療を行う際には、カテーテルの先端部を体外（例えば肘や足の付け根）から体内の血管に挿入し、狭窄した血管の位置までカテーテルの先端部を移動させ、位置決めを行う必要がある。カテーテルの先端部の位置及び経路の確認は、従来、Ｘ線撮影装置を用いた診断方法が用いられてきた。しかし、Ｘ線を用いると被曝の問題があるので、生体への影響の少ない超音波診断装置を適用することが研究されている。
【０００４】
特許文献１には、例えば前腕部を診断対象として前腕部の血管を診断するために、超音波を利用したボリュームレンダリング法に基づいた画像形成技術が記されている。特許文献２には、血管内に挿入されたカテーテルの位置を観察することのできる医用装置が示されている。この医用装置においては、三次元画像を得るために２台のＸ線撮影装置等を用いて、被検者の２方向からの投影画像を形成している。また、特許文献３には、血管に挿入するカテーテルの三次元画像を形成する手術支援方法が記されている。この方法においては、ユーザが手動で手術器具の形状及び初期位置を設定した上で血管内部の三次元画像が形成される。
【０００５】
【特許文献１】特開２００２−３３６２５３号公報
【特許文献２】特開２００２−１１９５０７号公報
【特許文献３】特開平１０−３３４２２０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
超音波診断装置を用いて三次元画像を取得するためには、特許文献１記載の装置を用いることにより血管の立体投影画像を得ることができる。血管の画像を明瞭に表示するためには、血管に対して適切なレンダリング条件が設定される。血管を表示対象としてレンダリング条件を最適化すると、血管の内部に何らかの器具が存在していたとしても、その器具の画像化は困難である。
【０００７】
一方において、血管内の器具を画像化することだけを目的にするのであれば、レンダリングパラメータである不透明度（オパシティ）の設定を一律で変更することによって画像化は可能である。カテーテルのような器具は一般的には高反射材であるため高レベルのエコー信号が得られる。しかし、不透明度のパラメータを一律に変化させることにより器具を明瞭に画像化しようすると、今度は血管を画像化することが困難になってしまう。
【０００８】
血管内の器具の状態を画像化することができれば、治療対象である血管と器具との位置を把握することができるのであるが、従来の超音波診断装置を用いて不透明度のパラメータを調整しても、血管と血管内に挿入された器具とを同時に明瞭に可視化することはできなかった。
【０００９】
本発明の目的は、血管と器具を両方とも同時に明瞭な三次元画像として可視化することのできる超音波診断装置を提供することにある。
【課題を解決するための手段】
【００１０】
本発明は、器具が挿入された血管に超音波を送受波して得られたボクセルデータ群に対して設定された各レイについて、深さ方向にボクセル演算を逐次実行する超音波診断装置において、前記各レイについて、第１開始点から前記ボクセル演算の逐次実行結果により、第１三次元画像を構成する画素値を求める第１レンダリング手段と、前記各レイについて第２開始点を判定する判定手段と、複数のレイの中で前記第２開始点が判定されたレイについて、前記第２開始点からのボクセル演算の逐次実行結果により、前記器具が反映された第２三次元画像を構成する画素値を求める第２レンダリング手段と、前記第１三次元画像と前記第２三次元画像とを合成し合成三次元画像を形成する合成処理手段と、前記合成三次元画像を表示する表示手段と、を有することを特徴とする。
【００１１】
上記構成によれば、ボクセルデータ群に対して設定された各レイ上において、第１開始点からのボクセル演算の逐次実行によって各々の画素値が求められる。すなわち、第１レンダリング処理によって各々の画素値が求められ第１三次元画像が形成される。また、ボクセルデータ群に対して設定された各レイ上において、第２開始点が判定された場合には、第２開始点より奥側のレイ上のデータに対して第２三次元画像の形成のための第２レンダリング処理が実行される。
【００１２】
第１レンダリング手段によっては第１三次元画像を形成することが可能となり、第２レンダリング手段によっては第２開始点を判定した後の複数のボクセルデータを用いて、血管内に挿入された器具を反映させた第２三次元画像を形成することが可能となる。つまり、第２開始点が存在する場合には、第１三次元画像と第２三次元画像との２つの画像を形成して、合成三次元画像を形成することにより、生体組織である血管と血管内に挿入された器具とのそれぞれを明瞭に三次元画像化することができる。
【００１３】
第２開始点の判定方法については、組織と器具のそれぞれの画像を明瞭に画像表示できる限りにおいて、様々な手法を適用することができる。すなわち本発明はレンダリングの段階的な切換えを特徴事項の一つとしており、その切換えの条件は各種の用途に応じて適宜定めればよい。
【００１４】
なお、三次元空間内において任意の視点位置からレイを設定してもよいし、レイと超音波ビームとが同一方向でもよい。ボクセル演算としては、不透明度（オパシティ）を利用したボリュームレンダリング法による透過光量演算を行ってもよいし、その他の三次元画像構成のための演算を用いてもよい。
【００１５】
望ましくは、前記判定手段は、前記血管内かつ前記器具の手前側の深さ位置を前記第２開始点として判定することを特徴とする。上記構成によれば、各レイの方向に対して、器具の手前側に第２開始点を設定することにより、血管内の器具に対する第２レンダリング処理を行うことができる。
【００１６】
望ましくは、前記判定手段は、前記各レイについての前記第１開始点からのボクセル演算の逐次実行結果に基づいて、前記第２開始点の判定を行うことを特徴とする。上記構成によれば、ボクセルデータ群に対して設定された各レイ上において、第１開始点からのボクセル演算の逐次実行結果から所定の条件を満たすことにより第２開始点が設定される。ボクセル演算を行う上でのボリュームレンダリングの演算条件は変更することができるので、例えば可視化対象である血管の位置でレンダリング処理を終了させて、その終了点をもって第２開始点とすることもできる。
【００１７】
望ましくは、前記判定手段は、前記各レイについての前記第１開始点からのボクセル演算の逐次実行結果に基づいて、前記第２開始点の判定を行うことを特徴とする。上記構成によれば、ボクセル演算によって逐次算出される値に基づいて、第２開始点の判定が行われるので、超音波のエコー強度、すなわち生体の組織内部の状態を反映した第２開始点の判定を行うことができる。
【００１８】
望ましくは、前記判定手段は、前記各レイについて逐次実行される前記第１開始点からのボクセル演算で用いられるパラメータの積算値に基づいて、前記第２開始点の判定を行うことを特徴とする。上記構成によれば、ボクセル演算で用いられるパラメータの積算値に基づいて、第２開始点の判定が行われるので、パラメータに設定される値に応じて、生体組織の状態に適合した柔軟な第２開始点の設定を行うことができる。ここでパラメータとしては、例えば、不透明度、透明度を示す値などを利用できる。
【００１９】
望ましくは、前記各レイ上のボクセルデータ列に対して輝度反転処理を施す前処理手段を含み、前記判定手段、前記第１レンダリング手段及び前記第２レンダリング手段の内の少なくとも１つは前記輝度反転処理が施されたボクセルデータ列に基づいて処理を実行することを特徴とする。
【００２０】
上記構成によれば、輝度反転処理が施されたボクセルデータ列が処理対象となる。超音波の受信信号強度を反転して画像形成処理を行い、表示手段で合成三次元画像を表示することによって、本来であれば低レベル信号が高レベル信号として置換処理され、逆に本来の高レベル信号が低レベル信号として置換処理される。従って、血液から得られる低レベルのエコー信号を高レベルのエコー信号として扱うことにより、血液が流れる部分についての信号処理することができる。
【００２１】
更に、器具からは一般的には高レベル信号が得られる場合が多いが、血管内に挿入された器具についても同様に高レベル信号が得られる。よって、通常の信号処置であれば器具が強調されるような信号処理になる。しかし、器具を主体としたボクセルデータ列についても輝度反転処理を行うことによって、器具を強調しないような処理を行うことも可能となる。このように、反転処理をする信号と反転処理をしない信号とを任意に選択することによって、三次元画像の表示の際に強調する部分や強調しない部分を選択して処理することが可能となる。
【００２２】
望ましくは、超音波を送受波して得られた受信信号を直交検波し、血液の移動に伴って発生するドプラ偏移周波数成分に応じて血流の存在位置を検出するカラードプラ処理手段を更に含み、前記判定手段は、前記カラードプラ処理手段が検出する血流位置情報に基づいて、前記第２開始点の判定を行うことを特徴とする。上記構成によれば、血液が生体内を移動することに伴って発生するドプラ偏移周波数に応じて、血液の存在位置を検出することができるので、確実に血流が存在する部分に第２レンダリングの第２開始点を定めることができる。
【００２３】
望ましくは、前記合成処理手段は、前記血管を主体的に表した第１三次元画像上に、前記器具を主体的に表した第２三次元画像を重ね合わせて前記合成三次元画像を形成することを特徴とする。上記構成によれば、第１三次元画像と第２三次元画像とを合成して三次元画像を形成しても、元々の第１三次元画像と第２三次元画像の画像情報を視覚的に識別することができる。第１三次元画像は生体組織を主体とした画像であり、第２三次元画像は器具を主体とした画像であるので、合成三次元画像においても器具の位置を明瞭に把握することができる。
【００２４】
望ましくは、前記合成処理手段は、前記第１三次元画像と前記第２三次元画像とを視覚的に異ならせつつ合成処理を行うことを特徴とする。上記構成によれば、第１三次元画像と第２三次元画像との識別が容易になる。視覚的に異ならせる合成処理としては、例えば着色表示、輝度を変動させる点滅表示あるいは輪郭強調表示などが挙げられる。
【発明の効果】
【００２５】
以上説明したように、本発明によれば、血管と器具を両方とも同時に明瞭な三次元画像として可視化することができる。
【発明を実施するための最良の形態】
【００２６】
以下、本発明の好適な実施形態を図面に基づいて説明する。図１は、人体の前腕部１０に対して超音波の送受信器を当接している状態を模式的に示した図である。本実施形態では、血液透析を受けている透析患者を超音波診断の対象としており、前腕部１０において狭窄が発生している血管１４に対して器具１６を挿入し、治療を行う場面を例として説明する。
【００２７】
この器具１６は、血管１４に挿入される手術用の器具であり、具体的にはカテーテルを示している。探触子本体１２はアレイ振動子（図示せず）を備える電子走査方式の超音波プローブであり、探触子本体１２の送受波面側には水袋１８が装着されている。探触子本体１２は、後述する機械走査機構２２（図示せず）によって保持されており、機械走査機構２２は、探触子本体１２と水袋１８とを合体させることによりプローブユニット３０として構成される。
【００２８】
図２は、本発明に係る超音波診断装置の概略的な全体構成を示すブロック図である。本実施形態における超音波診断装置は、大別してプローブユニット３０と装置本体４０とから構成される。まず、プローブユニット３０について更に詳述する。プローブユニット３０は、可動部２０と機械走査機構２２とから構成されている。可動部２０は、前述したように探触子本体１２と水袋１８とから構成されている。探触子本体１２の電子走査方式は、本実施形態においては電子リニア走査方式が用いられている。なお、電子セクタ走査方式やその他の電子走査方式を用いてもよい。水袋１８の外表面は、柔軟な膜状部材で形成されており、その内部には例えば蒸留水などが充填されている。
【００２９】
探触子本体１２を前腕部１０に当接させると、水袋１８が前腕部１０の曲面形状に適合して柔軟に変形し、アレイ振動子から送信される超音波を前腕部１０に送受信する超音波伝達経路が形成される。本実施形態においては、１本の超音波ビームの形成に応じて１つのエコーデータ列が取得される。エコーデータ列は、プローブユニット３０から装置本体４０側に出力される。超音波ビームを電子走査することによって、前腕の長手方向に対して垂直な方向に走査面１９（図１参照）が形成される。走査面上のエコーデータアレイを装置本体４０において取り込むことにより、１枚の走査面上での２次元断層情報を得ることができる。
【００３０】
機械走査機構２２は、機械走査制御部２４からの信号に基づいて位置決め走査される機構である。探触子本体１２を電子走査することによって走査面が形成されるが、機械走査機構２２は、その走査面と直交する方向に可動部２０を機械走査する。本実施形態においては、水袋１８を前腕部１０に当接した状態を保ちながら、機械走査機構２２によって可動部２０を腕の長手方向に沿って移動させる。そうすると、前腕の長手方向に渡って互いに平行な複数の走査面が順番に形成される。可動部２０を前腕の長手方向に移動させることによって、複数のエコーデータアレイからなるエコーデータ群を取得することが可能となり、前腕部１０の内部の３次元空間情報を得ることができる。
【００３１】
次に、装置本体４０の構成について説明する。装置本体４０は、キーボードやトラックボールなどからなる入力部４２を備えている。この入力部４２は、三次元超音波画像を形成する上で必要な条件設定やパラメータ設定、あるいは三次元超音波画像の色彩調整を行うための数値入力機能を有する。入力部４２での入力内容は主制御部４４にて読み取られ、主制御部４４は、本装置全体の動作制御を行う。主制御部４４によって動作制御される機械走査制御部２４は、可動部２０を制御するための信号を機械走査機構２２に出力する。
【００３２】
送受信部４６は、送信ビームフォーミング機能及び受信ビームフォーミング機能を備えている。つまり、送受信部４６は、アレイ振動子を構成する各振動素子に対して、それぞれの遅延時間を設定して超音波ビームを形成する送信ビームフォーミング機能を有し、受信信号の信号受信タイミングを調整し、複数の受信信号の位相を整合させる整相加算機能を含む受信ビームフォーミング機能を有している。なお、送受信部４６は、受信信号の強度を補正するダイナミックレンジ変換回路などの信号処理回路を備えている。送受信部４６から出力されるエコーデータ列は、三次元画像処理部４８に伝送される。ちなみに、送受信部４６から出力されるエコーデータ列は、三次元画像処理部４８に出力するのと同時に、図示されていないカラードプラ処理部に対して出力してもよい。
【００３３】
三次元画像処理部４８は、複数の三次元画像を形成する機能を有し、本実施形態においては、反転処理部５０、第１レンダリング開始点判定部５２、第１レンダリング処理部５４、第２レンダリング開始点判定部５６、第２レンダリング処理部５８を備えている。ちなみに、本実施形態においては、ボリュームレンダリングの視線であるレイは、超音波ビームの進行方向と同一の方向に設定される。従って、三次元画像処理部４８に入力されたエコーデータ列は、そのままボクセル演算の対象であるボクセルデータ列とみなして処理することができる。
【００３４】
反転処理部５０は、入力されるデータの輝度を反転して出力する機能を有する。反転処理部５０においては、エコーデータ列を構成する各エコーデータの信号レベルを反転変換する。この反転変換を行うと低レベルのエコーデータは高レベルのエコーデータに置換処理される。故に、反転処理部５０は、超音波に対する反射強度の弱い生体組織、つまり前腕部においては血液の流れる部分からのエコー信号の検出を容易にする手段となる。
【００３５】
第１レンダリング開始点判定部５２は、第１レンダリング処理部５４で行われるボクセル演算の第１開始点を判定する機能を有する。その第１開始点の判定には、反転処理部５０から出力される反転エコーデータ列が用いられる。第１開始点の判定条件は、例えば、当該超音波診断装置のオペレータが入力部４２を操作して行うＲＯＩの領域設定を利用することにより定められる。ちなみに、第１開始点の判定条件としては、特開２００２−３３６２５３号公報に記されている被検体表面とカップリング体との境界面を検出する方法を用いることもできる。
【００３６】
第１レンダリング処理部５４は、前述の第１レンダリング開始点判定部５２で第１開始点が判定されたことを信号処理の起点としてボクセル演算の逐次演算を行い、あるレイ上の三次元画像の画素値を求める機能を有する。ボクセル演算については、例えば特開平１０−３３５３８号公報に記されている公知のボクセル演算法を用いることができる。
【００３７】
ボクセル演算法の一例としては、ボクセルの不透明度とエコー値との積によって表される発光光量と、ボクセルの透明度とそのボクセルの入力光との積によって表される減衰光量とを加算することにより、当該ボクセルからの出力光量を求める演算が挙げられる。ボクセルへの入力光量は、そのボクセルに近接する一つ前のボクセルの出力光量と同値であるため、逐次演算を行うことによってレイ方向の画素の輝度値を求めることができる。但し、演算対象のボクセルがレイ上の最終ボクセルに到達したり、レイ上の複数のボクセルのオパシティ積算値が１になったりした場合には、逐次演算を終了するという終了条件が付加される。
【００３８】
第２レンダリング開始点判定部５６は、第２レンダリング処理部５８で行われるボクセル演算の第２開始点を判定する機能を有する。その第２開始点の判定には、本実施形態においては、第１レンダリング処理部５４で行われるボクセル演算の逐次実行結果が利用される。すなわち、第１レンダリング処理部５４で行われるボクセル演算による透過光量の値が、所定の閾値を越える場合に、そのボクセル演算の演算対象であるエコーデータを第２開始点とする。
【００３９】
本実施形態にかかる超音波診断装置は、第１レンダリング処理の結果に基づいて、第２レンダリング処理を実行することを特徴事項の一つとしており、その第２レンダリング処理は第２開始点から開始される。第２開始点は、器具を主体とした第２三次元画像を取得するためのレンダリング処理の再スタートポイントあるいはリセットポイントになる。血管内に挿入された器具を描出する第２三次元画像を形成するためには、血管壁の部分、更に望ましくは血管内を流れる血流部分に入り込んだ部分に第２開始点を設定する必要がある。
【００４０】
血管壁と血液のそれぞれのエコー強度が大きく変化するという物理現象を利用して、血管の内部に入り込む位置を第２開始点として検出することができる。血管の内部に入り込んだ位置である第２開始点から、超音波画像を形成するレンダリング処理を変えることは、血管内に存在する血液以外の異物を視覚化する上で非常に有効な方法であり、特に、三次元画像化の場合には、三次元空間内での器具の位置を認識するために有効である。
【００４１】
第２開始点の判定方法には、前述したような透過光量の値が所定の閾値を越えることを利用する方法以外にも、いくつかの他の判定方法を用いることもできる。
【００４２】
第２開始点の他の判定方法としては、例えば、ボクセル演算の逐次演算と同時に並列実行されるオパシティの積算演算を用いることができる。すなわち、ボクセル演算によるオパシティの積算値が所定の閾値を超える場合には、そのボクセル演算の演算対象であるエコーデータを第２開始点とすることができる。
【００４３】
また、第２開始点の他の判定方法としては、前述の特開２００２−３３６２５３号公報に記されている、被検体表面とカップリング体との境界面を検出する方法と同様の手法を用いて、第２開始点を判定することができる。その同様の手法とは、すなわち関心領域を被検体表面を包含しない位置に予め設定しておき、その被検体表面を包含しない関心領域内におけるデータを輝度反転させた後に第１レンダリング処理部５４にて処理する手法である。高強度の反転エコーデータが連続してＮ個（一例として５個）続いた場合に、そのＮ個目のエコーデータを第２開始点とすることができる。
【００４４】
更に、第２開始点の他の判定方法としては、連続するデータ列の微分値を用いることもできる。すなわち、ボクセル演算による画像再構築を行う場合にレイの方向に沿って、連続するデータ列の構成要素である複数のデータを対象に微分を行い、微分値が急激に変化するデータに相当する位置を第２開始点とすることができる。
【００４５】
以上に記したいくつかの検出方法のいずれかの方法を用いて、第２開始点を検出することができる。
【００４６】
第２レンダリング処理部５８は、前述の第２レンダリング開始点判定部５６で判定された第２開始点からボクセル演算の逐次演算を行って、あるレイ上の三次元画像の画素値を求める機能を有する。第２レンダリング処理部５８には、反転処理される前のそのままのエコーデータが入力される。第２レンダリング処理部５８で行われるボクセル演算は、第１レンダリング処理部５４で行われるボクセル演算と同じ方法であってよいし、そのボクセル演算とは異なる三次元画像の構築演算であってもよい。
【００４７】
次に、三次元画像合成部６０について説明する。三次元画像合成部６０は、２種類の三次元画像を重ね合わせて、１つの合成三次元画像を形成する機能を有する。三次元画像合成部６０は、第１三次元画像メモリ６２、第２三次元画像メモリ６４などの構成を備えている。第１三次元画像メモリ６２は、第１レンダリング処理部５４で求められる複数の画素値をマッピングして構成される第１三次元画像を記憶するメモリである。第２三次元画像メモリ６４は、第２レンダリング処理部５８で求められる複数の画素値をマッピングして構成される第２三次元画像を記憶するメモリである。これらのメモリの中の第１三次元画像と第２三次元画像とを合成することにより、１つの合成三次元画像が形成される。
【００４８】
表示部６８は、ＣＲＴやＬＣＤで構成されており、合成三次元画像をオペレータが観察できるように表示する。
【００４９】
図３は、三次元画像処理部４８で行われる第１及び第２レンダリング処理の各演算処理に関する説明図である。図３には、可動部２０を前腕部１０に当接させて得られる三次元画像取り込み領域が模式的に示してある。
【００５０】
可動部２０を移動させることで互いに平行な複数の走査面が形成される。最前面の走査面７０上には、血管１４の断面と血管に挿入された器具１６が含まれる。また、走査面７０上には血管１４を中心にした三次元画像を形成するための関心領域（ＲＯＩ）８０の設定範囲が示されている。また、走査面７０上で可動部２０に近接した部分には、前腕部１０の体表面と水袋１８とが当接する境界部分８２が示されている。ちなみに、ＲＯＩ８０の設定形状は、Ｙ軸方向に伸長した三次元の直方体形状である。
【００５１】
図３においては、１本の超音波ビームに応じて取得されるボクセルデータ列を、走査面上のＺ軸に平行な１本の直線として表している。例えば、ボクセルデータ列Ａ７２は、血管１４が存在しない位置で取得されるボクセルデータ列である。
【００５２】
ボクセルデータ列Ａ７２は、ＲＯＩ８０の設定領域との関係において第１開始点が判定される。その第１開始点の位置は、図３に示す走査面７０においてボクセルデータ列Ａ７２を示す直線とＲＯＩ８０の上端面とが交差する点Ｐ１と定めることができる。
【００５３】
ボクセルデータ列Ａ７２に対するボクセル演算は、第１開始点である点Ｐ１の位置から、ＲＯＩ８０の下端面と交差する点Ｐ２までの範囲のデータ列に関して演算される。累積的にボクセル演算を実行している途中において、演算の終了条件を満たした場合には演算が中止され、そのときの演算結果がボクセルデータ列Ａ７２に関する画素値となる。または、ＲＯＩ８０の下端面上の点Ｐ２に至っても演算の終了条件が最後まで満たさない場合には、最終ボクセルでの演算結果が画素値となる。つまり、第１のレンダリング処理を行うにあたって、演算の終了条件を満足する／しないに関わらず、１つのボクセルデータ列に対して１つの画素値が求められる。その画素値は、図３に示すビーム対応画素値Ａ（８６）となる。第１のレンダリング処理を実行している最中に、第２開始点の判定条件を満足する点が検出されれば、第２のレンダリング処理に切り替わるが、ボクセルデータ列Ａ７２は、血管が存在しない部分で取得されるボクセルデータ列であるので、第２開始点は検出されない。
【００５４】
次に、血管１４と器具１６とが存在する位置で取得されるボクセルデータ列Ｂ８８に関する演算処理について説明する。ボクセルデータ列Ｂ８８は、前述のボクセルデータ列Ａ７２と同様に、ＲＯＩ８０の設定領域との関係において第１開始点が判定される。その第１開始点の位置は、図３に示す走査面７０においてボクセルデータ列Ｂ８８を示す直線とＲＯＩ８０の上端面とが交差する点Ｐ３と定めることができる。ボクセルデータ列Ｂ８８に対しては、第１開始点であるＰ３が判定された後から、第１レンダリング処理が行われる。
【００５５】
第１レンダリング処理を順に実行していくと、血管１４の内側に入った位置におけるボクセルデータの演算が行われる。血液からは弱いエコーデータが得られるが、反転処理によって高いエコーデータに変換されている。従って、第１レンダリング処理の出力光量値がある一定の閾値を越えるという判定条件を用いて、この条件を満たす点が判定される。なお、第１レンダリング処理の出力光量値がある一定の閾値を越えたという判定条件を満たした時点をもって、第１レンダリング処理を中止することも可能である。その場合には、第１レンダリング処理の終了点である第１終了点Ｐ４は、前述の第２開始点Ｐ５と同一の点となる。第１開始点Ｐ３から第１終了点Ｐ４までの間のボクセルデータに対して、第１レンダリング処理を行うことによって、図３に示すビーム対応画素値Ｂ１（９０）が求められる。
【００５６】
次に、第１レンダリング処理の第１終了点Ｐ４に基づいて、第２レンダリング処理の開始点である第２開始点Ｐ５が設定される。第２開始点Ｐ５は、血管１４の内側に入った位置に設定される。第２開始点Ｐ５の設定に伴って第２レンダリング処理が実行される。第２レンダリング処理では、反転されていないボクセルデータ列Ｂ８８が演算処理の対象となる。第２レンダリング処理は、血管内の器具を画像化するための処理である。器具の画像化のためには、器具を反映した高レベルのエコーデータにレンダリングパラメータが設定される。
【００５７】
第２レンダリング処理を順に実行していくと、器具１６の表面から強いレベルのエコーデータが検出される。ボクセルデータ列Ｂ８８の上であって、器具１６の表面に交差する位置を第２レンダリング処理の終了点Ｐ６として設定することができる。第２開始点Ｐ５から第２終了点Ｐ６までの間のボクセルデータに対して、第２レンダリング処理を行うことによって、図３に示すビーム対応画素値Ｂ２（９２）が求められる。
【００５８】
なお、一般的には、器具の表面からは高レベルのエコー信号が得られるので、器具１６の表面上において第２レンダリング処理が終了するように設定することができる。但し、ボリュームレンダリングのパラメータ設定の如何によっては、更に深さ方向に存在するボクセルデータの演算を累積することもできる。
【００５９】
以上説明したように、２段階のレンダリング処理を行うことにより、第１レンダリング処理でビーム対応画素値Ａ（８６）とビーム対応画素値Ｂ１（９０）が求められる。そして、第２レンダリング処理によりビーム対応画素値Ｂ２（９２）が求められる。例示した２つのボクセルデータ列Ａ７２とボクセルデータ列Ｂ８８のみならず、ＲＯＩ８０内に存在する全てのボクセルデータ列に対して第１レンダリング処理を行うことにより、第１三次元画像９４が形成される。また、第２開始点が判定されたボクセルデータ列に対して第２レンダリング処理を行うことにより、第２三次元画像９６が形成される。
【００６０】
第１三次元画像９４は、血管が存在しない生体組織に相当するボクセルデータ列と、血管が存在する部分については血管よりも手前側の生体組織に相当するボクセルデータ列に基づいて形成される。つまり、第１三次元画像９４は血管の内部の情報を含まず、生体組織を主体的に表した三次元投影画像となる。また、第２三次元画像９６は、血管あるいはその血管の裏側に位置する生体組織に相当するボクセルデータ列に基づいて形成される。よって、第２三次元画像９６では、器具投影部分９８を認識することができる。
【００６１】
図４は、図３に示す走査面７０内において、第１レンダリング処理と第２レンダリング処理が適応される領域を区別して示した図である。図４には、図３に示す走査面７０内のＲＯＩ８０の領域だけが示してある。
【００６２】
ＲＯＩ８０の上端には、第１レンダリング処理が開始される第１開始点群１６２が直線として示してある。前述の点Ｐ３も、第１開始点群１６２の上に存在する。また、ＲＯＩ８０の下端には、第１終了点群（１６４Ａ、１６４Ｂ）が２つの線分として示してある。血管１４が存在する部分は、第１レンダリング処理が血管内で終了するようにレンダリング条件が設定されているので、血管に対する第１終了点群１６６は血管壁の形状に対応した円弧形状となる。第１のレンダリング処理が施される領域は、第１開始点群１６２と第１終了点群（１６４Ａ，１６４Ｂ，１６６）とに挟まれる領域であるので、図５において図示する逆凹型の斜線領域１６０となる。
【００６３】
本実施形態においては、第２レンダリング処理は、血管が存在する位置から得られるボクセルデータ列に対して実行される。よって、血管内の第１終了点群である１６６が第２開始点群１６８となる。第２レンダリング処理は、レイ１７０の深さ方向に沿って、第２開始点群１６８から開始されるので、第２レンダリング処理が施される領域は、図５において図示する斜線領域１７２となる。ちなみに、器具１６の表面の位置で第２レンダリング処理を中止することもできる。
【００６４】
図５は、ある一つのエコーデータ列を信号処理することで生成する各種信号を並べて示した図である。図５の最上段に示す図５（ａ）は、エコーデータ列の強度レベルを示している。図５（ａ）の横軸は時間軸であり、縦軸はエコーデータの強度レベルを示す。このエコーデータ列は、血管内に存在する器具に向かって１本の超音波ビームを送信することで得られたエコーデータ列を表している。図３に示したボクセルデータ列Ｂ８８がこのエコーデータ列に相当する。
【００６５】
エコーデータ列１００は、そのデータ列の中間部分１０２において、レベルが急激に低下したり（１０４ａ，１０４ｂ）、一時的にピークを示したり（１０６）している。レベルが急激に低下しているのは、この部分のエコーデータが血管内の低いエコー強度を反映しているためであり、レベルが一時的にピークを示すのは、このピーク部分のエコーデータが器具によるエコー強度を反映しているためである。なお、中間部分１０２に含まれず時間軸上で前後に広がるエコーデータは、生体組織によって得られるエコーデータである。
【００６６】
図５（ｂ）は、エコーデータ列１００を反転処理して得られた反転エコーデータ列１０８を示すグラフである。図５（ｂ）の横軸と縦軸の定義は図５（ａ）と同一である。反転処理を行うことにより、エコーデータ列１００において低レベルである部分（１０４ａ，１０４ｂ）が、高レベルのデータに変換される。
【００６７】
ここで、図５の説明を一時的に離れて、図６を用いてエコーデータの反転処理特性の一例を示す。図６は、エコーデータの反転処理特性の一例を示す図である。横軸はエコー信号入力レベルを示し、縦軸はエコー信号出力レベルを示している。この曲線１２４に従った変換処理を行うと、低レベルの入力信号は高レベルの出力信号として変換され、逆に高レベルの入力信号は低レベルの出力信号として変換される。なお、グラフは非線形特性を示しているので、入力信号のレベルが低い程、高いレベルに強調された反転信号が得られるようになっている。また、ある閾値を上回る入力信号については、出力信号として変換されないカットオフ特性も有している。
【００６８】
図５に戻って、反転エコーデータ列１０８の信号処理について説明する。本実施形態に係る超音波診断装置においては、この反転エコーデータ列１０８について、第１レンダリング処理を施している。第１レンダリング処理部５４でボクセル演算を実行するのに伴って、各ボクセルデータからの出力光量値が求められ、それと平行して同時に不透明度（オパシティ）の積算値も求められる。
【００６９】
図５（ｃ）は第１レンダリング処理によるオパシティの積算値（Σαi）の累積的な変
化を示し、図５（ｄ）は第１レンダリング処理による出力光量値（Ｃ１out）の累積的な変化を示している。本実施形態においては、出力光量値（Ｃ１out）の値が一定の閾値を超えた場合に、第２レンダリング処理が開始される。よって、図５（ｄ）において出力光量値を表す曲線１１２と、閾値Ｃthを表す水平線との交差点Ｔが存在することが、第２レンダリング処理の開始を意味する。本実施形態においては、第１レンダリング処理は、第２レンダリング処理が開始されたことをもって、その演算処理を中断する。但し、第２レンダリング処理の開始に関わらず、第１レンダリング処理の演算を終了条件を満たすまで継続してもよい。第１レンダリング処理は第２レンダリング処理とは別個独立に実行されるからである。図５（ｄ）においては、第１レンダリング処理の中断又は継続が任意であることを破線によって示している。ボクセル演算はオパシティの積算演算と表裏一体に実行されるので、図５（ｃ）においてもその積算値を表す曲線の一部を破線で示している。
【００７０】
図５（ｅ）は、第２レンダリング処理の処理対象となるエコーデータ列１１８を示している。このエコーデータ列１１８は、図５（ａ）に示すエコーデータ列１００の後半部、つまり時間軸上にて一点鎖線１１６で示す時点以降で得られたエコーデータ列と同値である。第２レンダリング処理に関して注目すべきことは、処理対象のデータ列は反転エコーデータ列ではないことと、第２レンダリング処理のための第２のオパシティが新たにゼロから積算されることである。第２レンダリング処理を行う場合に、反転エコーデータ列を用いると、最初から高レベルのデータを処理することになり、器具の画像化に適さない。第２レンダリング処理は、ある判定条件を満たした場合に開始されるので、判定条件を満足した時点からオパシティの積算演算が新たに開始される。第２レンダリング処理を行うことによって、エコーデータ列１１８の中でピークを示している器具のエコーデータ１２０に基づいて器具に対応する画素値が求められる。
【００７１】
図５（ｆ）には、第２レンダリング処理による出力光量値Ｃ２out（１２２）の変化が示されている。第２レンダリング処理では、器具を明瞭に画像化することを優先してレンダリングパラメータを設定することが好ましい。従って、この出力光量値Ｃ２outについても、器具または血管内に相当する時間的な範囲での出力光量値の変化は大きいが、それ以降は飽和傾向になるため時間的な変化は少なくなる。
【００７２】
図７は、他の実施形態に係る超音波診断装置の構成ブロック図である。図２の構成ブロック図と比較して、三次元画像処理部１３０の構成が異なっている。よって、図７において、図２に示した構成と同様の構成には同一符号を付しその説明を省略する。
【００７３】
三次元画像処理部１３０においては、第２レンダリング開始点判定部５６に対して、反転されていないエコーデータ列が入力される。第２レンダリング開始点判定部５６では、低レベルのエコーデータが連続してＭ個（一例として５個）続いた場合に、そのＭ個目のエコーデータを第２開始点とする判定処理が行われる。低レベルのエコーデータが検出された後に、引き続いてＭ個のデータが連続しなければ判定条件が成立しないので、血管内に確実に入り込んだ位置に第２開始点を設定することができる。判定条件を満たす第２開始点が検出されると、その情報は、第１レンダリング処理部５４に伝達される。この情報伝達の方向は、図２に示した機能ブロック図で示した情報伝達の方向と逆である。図７に示す第１レンダリング処理部５４においては、第２開始点の検出をもって、その検出の直前まで実行していたボクセル演算を中止する。一方、本実施形態によると、第２開始点が検出されることをもって、生体組織ではなく血管内でのレンダリング処理に移行していることを確実に判断できる。よって、三次元画像処理部１３０によれば、第２開始点よりも更に深い位置から得られるエコーデータの情報に基づいて、第１レンダリング処理が実質的に処理不要になる位置を検出することにより、第１レンダリング処理を強制的に中止することができる。
【００７４】
また、他の実施形態に係る超音波診断装置においては、図７に示す超音波診断装置の装置本体４０が備えているカラードプラ処理部（図示せず）を用いて、以下に示す手法により第２開始点を検出することもできる。
【００７５】
カラードプラ処理部は、送受信部４６が出力する受信信号を入力して、その分析結果を第２レンダリング開始点判定部５６に対して出力する。カラードプラ処理部においては、受信信号を直交検波した後にＭＴＩフィルタにより低いドプラ偏移周波数を除去し、自己相関演算することでドプラ偏移周波数の平均周波数とパワーを検出することができる。このようにドプラ偏移周波数の平均周波数とパワーに基づいて血流のイメージング表示を行う技術は、いわゆるカラーフロー又はパワーフローと呼ばれる公知の技術である。この公知技術を用いれば、超音波断層画像上で血流の存在する部分と存在しない部分の境界を判定することができる。その判定結果、つまりカラードプラ処理部で求めた各レイ上での血流の存在位置を示す境界情報を第２レンダリング開始点判定部５６に伝送することによって、第２レンダリング開始点判定部５６において第２開始点を定めることができる。この判定方法によれば、血流の存在する血管内部の位置に確実に第２開始点を設定することができる利点がある。
【００７６】
カラードプラ処理部にて行われる処理は、カラーフロー処理よりもパワーフロー処理の方が好ましい。パワーフロー処理は、カラーフロー処理に比べてドプラアングルの依存性が少ないため、血管の走行を連続的に捉えられる。なお、一般的に、器具が血管内を移動する速度は血液が血管内を流れる速度に比べて無視できる程度に小さいので、器具の運動によってドプラ偏移周波数が影響を受けることはない。
【００７７】
図８は、他の実施形態に係る三次元画像合成部の構成ブロック図である。図８に示す三次元画像合成部１４０は、図２に示す三次元画像合成部６０の変形例である。三次元画像合成部１４０は、三次元画像合成部６０と置換して用いることが可能であり、２枚の三次元画像を合成する際に、各画像に輝度の重み付けを付加する機能を有している。三次元画像合成部１４０は、第１表示処理部１４２と第２表示処理部１４４を備える点において、前述の三次元画像合成部６０と相違する。第１表示処理部１４２は、第１三次元画像の輝度と色彩を変換あるいは調整する機能を有しており、その機能を使用するために、例えば輝度変換係数Ｊ（設定可能範囲：−１≦Ｊ≦１）のパラメータが活用される。また、第２表示処理部１４４は、第２三次元画像の輝度と色彩を変換あるいは調整する機能を有しており、その機能を使用するために、例えば輝度変換係数Ｋ（設定可能範囲：−１≦Ｋ≦１）のパラメータが活用される。ちなみに、輝度変換係数Ｊ及びＫは、入力部４２からのオペレータの操作入力信号に基づいて、主制御部４４から出力され、三次元画像合成部１４０に対して入力される係数である。
【００７８】
輝度変換係数Ｊ及びＫのそれぞれの機能は同一である。すなわち、輝度変換係数が１（Ｊ＝１，Ｋ＝１）である場合には三次元画像の輝度調整を一切行わず、輝度変換係数が−１（Ｊ＝−１，Ｋ＝−１）の場合には輝度反転処理を行う。そして、輝度変換係数が１から−１の間の小数値が指令された場合には輝度のグラデーション調整が行われる。
【００７９】
このような輝度反転及び調整機能を有することにより、第１三次元画像と第２三次元画像とは、それぞれが独立に輝度調整を行うことができる。よって、血流を明瞭に画像化するためにレンダリング処理前にエコーデータの反転処理が行われる場合であっても、合成三次元画像においては生体組織と器具とを明暗により、あるいは色彩の違いによって識別することが可能となり、双方とも明瞭に三次元画像として形成することができる。
【００８０】
図９は、三次元画像合成部１４０を用いて実行される画像合成の説明図である。図２にて説明した三次元画像合成部６０、あるいは図８を用いて説明した三次元画像合成部１４０のどちらにおいても、第１三次元画像と第２三次元画像との画像合成が行われる。画像合成は、第１三次元画像の上に第２三次元画像を重ね合わせる、いわゆるオーバーレイ処理によって実行される。図９に示す第１三次元画像１４６には、生体組織内で屈曲した状態の血管１５２が示されている。この第１三次元画像１４６は、三次元のレンダリング処理を施して形成された三次元投影画像である。血管１５２の部分が白く表示されているのは、輝度反転したエコーデータ群によって第１三次元画像１４６が形成されるためである。第２三次元画像１４８には、血管に挿入された器具１５４の形状が示されている。この第２三次元画像１４８も三次元投影画像である。器具１５４の部分が白く表示されているのは、輝度反転していないエコーデータ群によって第２三次元画像１４８が形成されるためである。図９に示した状態のまま、血管１５２と器具１５４の双方が白く表示された状態で２枚の三次元画像を合成処理すると、血管内の器具１５４の位置を識別することが難しくなる。識別を容易にするために、図８を用いて説明した輝度変換係数ＪとＫをＪ＝１，Ｋ＝−１と設定することにより、第２三次元画像１４８だけについて反転処理が実行される。そうすると、器具１５４の部分が黒く変換されるため、合成三次元画像１５０においては、血管１５２の立体形状が明確に表示されると同時に、血管内の器具１５６の位置も明確に把握することができる。なお、第２三次元画像に対して着色表示を施せば、器具１５４の部分を生体組織と異なる色に表示することが可能となり、より明瞭に両者を識別することができる。
【図面の簡単な説明】
【００８１】
【図１】前腕部に超音波の送受信器を当接している状態を模式的に示した図である。
【図２】本発明に係る超音波診断装置の概略的な全体構成を示すブロック図である。
【図３】三次元画像処理部で行われる第１及び第２レンダリング処理の各演算処理に関する説明図である。
【図４】第１及び第２レンダリング処理の領域の区別を視覚的に示した図である。
【図５】ある一つのエコーデータ列の２信号処理についての説明図である。
【図６】エコーデータの反転処理特性の一例を示す図である。
【図７】他の実施形態に係る超音波診断装置の構成ブロック図である。
【図８】他の実施形態に係る三次元画像合成部の構成ブロック図である。
【図９】三次元画像合成部において実行される画像合成の説明図である。
【符号の説明】
【００８２】
１４ 血管、１６ 器具、２０ 可動部、７０ 走査面、７２ ボクセルデータＡ、８０ 関心領域（ＲＯＩ）、８２ 境界部分、８６ ビーム対応画素値Ａ、９０ ビーム対応画素値Ｂ１、９２ ビーム対応画素値Ｂ２、９４ 第１三次元画像、 ９６ 第２三次元画像、９８ 器具投影部分、Ｐ１，Ｐ３ 第１開始点、Ｐ２，Ｐ４ 第１終了点、Ｐ５ 第２開始点、Ｐ６ 第２終了点。
【技術分野】
【０００１】
本発明は超音波診断装置に関し、特に、三次元超音波画像を形成する超音波診断装置に関する。
【背景技術】
【０００２】
血液透析を行う場合に、透析患者の体から血液を取り出して体内に戻すためには、前腕部の血管が使用されることが多い。血液透析を長く続けていると、前腕部の血管が痛んで狭窄が発生することがある。
【０００３】
狭窄が発生して内径が細くなった血管を治療するために、血管内にバルーンやステントを挿入し血管の内径を拡げる治療が行われている。この治療を行う際には、カテーテルの先端部を体外（例えば肘や足の付け根）から体内の血管に挿入し、狭窄した血管の位置までカテーテルの先端部を移動させ、位置決めを行う必要がある。カテーテルの先端部の位置及び経路の確認は、従来、Ｘ線撮影装置を用いた診断方法が用いられてきた。しかし、Ｘ線を用いると被曝の問題があるので、生体への影響の少ない超音波診断装置を適用することが研究されている。
【０００４】
特許文献１には、例えば前腕部を診断対象として前腕部の血管を診断するために、超音波を利用したボリュームレンダリング法に基づいた画像形成技術が記されている。特許文献２には、血管内に挿入されたカテーテルの位置を観察することのできる医用装置が示されている。この医用装置においては、三次元画像を得るために２台のＸ線撮影装置等を用いて、被検者の２方向からの投影画像を形成している。また、特許文献３には、血管に挿入するカテーテルの三次元画像を形成する手術支援方法が記されている。この方法においては、ユーザが手動で手術器具の形状及び初期位置を設定した上で血管内部の三次元画像が形成される。
【０００５】
【特許文献１】特開２００２−３３６２５３号公報
【特許文献２】特開２００２−１１９５０７号公報
【特許文献３】特開平１０−３３４２２０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
超音波診断装置を用いて三次元画像を取得するためには、特許文献１記載の装置を用いることにより血管の立体投影画像を得ることができる。血管の画像を明瞭に表示するためには、血管に対して適切なレンダリング条件が設定される。血管を表示対象としてレンダリング条件を最適化すると、血管の内部に何らかの器具が存在していたとしても、その器具の画像化は困難である。
【０００７】
一方において、血管内の器具を画像化することだけを目的にするのであれば、レンダリングパラメータである不透明度（オパシティ）の設定を一律で変更することによって画像化は可能である。カテーテルのような器具は一般的には高反射材であるため高レベルのエコー信号が得られる。しかし、不透明度のパラメータを一律に変化させることにより器具を明瞭に画像化しようすると、今度は血管を画像化することが困難になってしまう。
【０００８】
血管内の器具の状態を画像化することができれば、治療対象である血管と器具との位置を把握することができるのであるが、従来の超音波診断装置を用いて不透明度のパラメータを調整しても、血管と血管内に挿入された器具とを同時に明瞭に可視化することはできなかった。
【０００９】
本発明の目的は、血管と器具を両方とも同時に明瞭な三次元画像として可視化することのできる超音波診断装置を提供することにある。
【課題を解決するための手段】
【００１０】
本発明は、器具が挿入された血管に超音波を送受波して得られたボクセルデータ群に対して設定された各レイについて、深さ方向にボクセル演算を逐次実行する超音波診断装置において、前記各レイについて、第１開始点から前記ボクセル演算の逐次実行結果により、第１三次元画像を構成する画素値を求める第１レンダリング手段と、前記各レイについて第２開始点を判定する判定手段と、複数のレイの中で前記第２開始点が判定されたレイについて、前記第２開始点からのボクセル演算の逐次実行結果により、前記器具が反映された第２三次元画像を構成する画素値を求める第２レンダリング手段と、前記第１三次元画像と前記第２三次元画像とを合成し合成三次元画像を形成する合成処理手段と、前記合成三次元画像を表示する表示手段と、を有することを特徴とする。
【００１１】
上記構成によれば、ボクセルデータ群に対して設定された各レイ上において、第１開始点からのボクセル演算の逐次実行によって各々の画素値が求められる。すなわち、第１レンダリング処理によって各々の画素値が求められ第１三次元画像が形成される。また、ボクセルデータ群に対して設定された各レイ上において、第２開始点が判定された場合には、第２開始点より奥側のレイ上のデータに対して第２三次元画像の形成のための第２レンダリング処理が実行される。
【００１２】
第１レンダリング手段によっては第１三次元画像を形成することが可能となり、第２レンダリング手段によっては第２開始点を判定した後の複数のボクセルデータを用いて、血管内に挿入された器具を反映させた第２三次元画像を形成することが可能となる。つまり、第２開始点が存在する場合には、第１三次元画像と第２三次元画像との２つの画像を形成して、合成三次元画像を形成することにより、生体組織である血管と血管内に挿入された器具とのそれぞれを明瞭に三次元画像化することができる。
【００１３】
第２開始点の判定方法については、組織と器具のそれぞれの画像を明瞭に画像表示できる限りにおいて、様々な手法を適用することができる。すなわち本発明はレンダリングの段階的な切換えを特徴事項の一つとしており、その切換えの条件は各種の用途に応じて適宜定めればよい。
【００１４】
なお、三次元空間内において任意の視点位置からレイを設定してもよいし、レイと超音波ビームとが同一方向でもよい。ボクセル演算としては、不透明度（オパシティ）を利用したボリュームレンダリング法による透過光量演算を行ってもよいし、その他の三次元画像構成のための演算を用いてもよい。
【００１５】
望ましくは、前記判定手段は、前記血管内かつ前記器具の手前側の深さ位置を前記第２開始点として判定することを特徴とする。上記構成によれば、各レイの方向に対して、器具の手前側に第２開始点を設定することにより、血管内の器具に対する第２レンダリング処理を行うことができる。
【００１６】
望ましくは、前記判定手段は、前記各レイについての前記第１開始点からのボクセル演算の逐次実行結果に基づいて、前記第２開始点の判定を行うことを特徴とする。上記構成によれば、ボクセルデータ群に対して設定された各レイ上において、第１開始点からのボクセル演算の逐次実行結果から所定の条件を満たすことにより第２開始点が設定される。ボクセル演算を行う上でのボリュームレンダリングの演算条件は変更することができるので、例えば可視化対象である血管の位置でレンダリング処理を終了させて、その終了点をもって第２開始点とすることもできる。
【００１７】
望ましくは、前記判定手段は、前記各レイについての前記第１開始点からのボクセル演算の逐次実行結果に基づいて、前記第２開始点の判定を行うことを特徴とする。上記構成によれば、ボクセル演算によって逐次算出される値に基づいて、第２開始点の判定が行われるので、超音波のエコー強度、すなわち生体の組織内部の状態を反映した第２開始点の判定を行うことができる。
【００１８】
望ましくは、前記判定手段は、前記各レイについて逐次実行される前記第１開始点からのボクセル演算で用いられるパラメータの積算値に基づいて、前記第２開始点の判定を行うことを特徴とする。上記構成によれば、ボクセル演算で用いられるパラメータの積算値に基づいて、第２開始点の判定が行われるので、パラメータに設定される値に応じて、生体組織の状態に適合した柔軟な第２開始点の設定を行うことができる。ここでパラメータとしては、例えば、不透明度、透明度を示す値などを利用できる。
【００１９】
望ましくは、前記各レイ上のボクセルデータ列に対して輝度反転処理を施す前処理手段を含み、前記判定手段、前記第１レンダリング手段及び前記第２レンダリング手段の内の少なくとも１つは前記輝度反転処理が施されたボクセルデータ列に基づいて処理を実行することを特徴とする。
【００２０】
上記構成によれば、輝度反転処理が施されたボクセルデータ列が処理対象となる。超音波の受信信号強度を反転して画像形成処理を行い、表示手段で合成三次元画像を表示することによって、本来であれば低レベル信号が高レベル信号として置換処理され、逆に本来の高レベル信号が低レベル信号として置換処理される。従って、血液から得られる低レベルのエコー信号を高レベルのエコー信号として扱うことにより、血液が流れる部分についての信号処理することができる。
【００２１】
更に、器具からは一般的には高レベル信号が得られる場合が多いが、血管内に挿入された器具についても同様に高レベル信号が得られる。よって、通常の信号処置であれば器具が強調されるような信号処理になる。しかし、器具を主体としたボクセルデータ列についても輝度反転処理を行うことによって、器具を強調しないような処理を行うことも可能となる。このように、反転処理をする信号と反転処理をしない信号とを任意に選択することによって、三次元画像の表示の際に強調する部分や強調しない部分を選択して処理することが可能となる。
【００２２】
望ましくは、超音波を送受波して得られた受信信号を直交検波し、血液の移動に伴って発生するドプラ偏移周波数成分に応じて血流の存在位置を検出するカラードプラ処理手段を更に含み、前記判定手段は、前記カラードプラ処理手段が検出する血流位置情報に基づいて、前記第２開始点の判定を行うことを特徴とする。上記構成によれば、血液が生体内を移動することに伴って発生するドプラ偏移周波数に応じて、血液の存在位置を検出することができるので、確実に血流が存在する部分に第２レンダリングの第２開始点を定めることができる。
【００２３】
望ましくは、前記合成処理手段は、前記血管を主体的に表した第１三次元画像上に、前記器具を主体的に表した第２三次元画像を重ね合わせて前記合成三次元画像を形成することを特徴とする。上記構成によれば、第１三次元画像と第２三次元画像とを合成して三次元画像を形成しても、元々の第１三次元画像と第２三次元画像の画像情報を視覚的に識別することができる。第１三次元画像は生体組織を主体とした画像であり、第２三次元画像は器具を主体とした画像であるので、合成三次元画像においても器具の位置を明瞭に把握することができる。
【００２４】
望ましくは、前記合成処理手段は、前記第１三次元画像と前記第２三次元画像とを視覚的に異ならせつつ合成処理を行うことを特徴とする。上記構成によれば、第１三次元画像と第２三次元画像との識別が容易になる。視覚的に異ならせる合成処理としては、例えば着色表示、輝度を変動させる点滅表示あるいは輪郭強調表示などが挙げられる。
【発明の効果】
【００２５】
以上説明したように、本発明によれば、血管と器具を両方とも同時に明瞭な三次元画像として可視化することができる。
【発明を実施するための最良の形態】
【００２６】
以下、本発明の好適な実施形態を図面に基づいて説明する。図１は、人体の前腕部１０に対して超音波の送受信器を当接している状態を模式的に示した図である。本実施形態では、血液透析を受けている透析患者を超音波診断の対象としており、前腕部１０において狭窄が発生している血管１４に対して器具１６を挿入し、治療を行う場面を例として説明する。
【００２７】
この器具１６は、血管１４に挿入される手術用の器具であり、具体的にはカテーテルを示している。探触子本体１２はアレイ振動子（図示せず）を備える電子走査方式の超音波プローブであり、探触子本体１２の送受波面側には水袋１８が装着されている。探触子本体１２は、後述する機械走査機構２２（図示せず）によって保持されており、機械走査機構２２は、探触子本体１２と水袋１８とを合体させることによりプローブユニット３０として構成される。
【００２８】
図２は、本発明に係る超音波診断装置の概略的な全体構成を示すブロック図である。本実施形態における超音波診断装置は、大別してプローブユニット３０と装置本体４０とから構成される。まず、プローブユニット３０について更に詳述する。プローブユニット３０は、可動部２０と機械走査機構２２とから構成されている。可動部２０は、前述したように探触子本体１２と水袋１８とから構成されている。探触子本体１２の電子走査方式は、本実施形態においては電子リニア走査方式が用いられている。なお、電子セクタ走査方式やその他の電子走査方式を用いてもよい。水袋１８の外表面は、柔軟な膜状部材で形成されており、その内部には例えば蒸留水などが充填されている。
【００２９】
探触子本体１２を前腕部１０に当接させると、水袋１８が前腕部１０の曲面形状に適合して柔軟に変形し、アレイ振動子から送信される超音波を前腕部１０に送受信する超音波伝達経路が形成される。本実施形態においては、１本の超音波ビームの形成に応じて１つのエコーデータ列が取得される。エコーデータ列は、プローブユニット３０から装置本体４０側に出力される。超音波ビームを電子走査することによって、前腕の長手方向に対して垂直な方向に走査面１９（図１参照）が形成される。走査面上のエコーデータアレイを装置本体４０において取り込むことにより、１枚の走査面上での２次元断層情報を得ることができる。
【００３０】
機械走査機構２２は、機械走査制御部２４からの信号に基づいて位置決め走査される機構である。探触子本体１２を電子走査することによって走査面が形成されるが、機械走査機構２２は、その走査面と直交する方向に可動部２０を機械走査する。本実施形態においては、水袋１８を前腕部１０に当接した状態を保ちながら、機械走査機構２２によって可動部２０を腕の長手方向に沿って移動させる。そうすると、前腕の長手方向に渡って互いに平行な複数の走査面が順番に形成される。可動部２０を前腕の長手方向に移動させることによって、複数のエコーデータアレイからなるエコーデータ群を取得することが可能となり、前腕部１０の内部の３次元空間情報を得ることができる。
【００３１】
次に、装置本体４０の構成について説明する。装置本体４０は、キーボードやトラックボールなどからなる入力部４２を備えている。この入力部４２は、三次元超音波画像を形成する上で必要な条件設定やパラメータ設定、あるいは三次元超音波画像の色彩調整を行うための数値入力機能を有する。入力部４２での入力内容は主制御部４４にて読み取られ、主制御部４４は、本装置全体の動作制御を行う。主制御部４４によって動作制御される機械走査制御部２４は、可動部２０を制御するための信号を機械走査機構２２に出力する。
【００３２】
送受信部４６は、送信ビームフォーミング機能及び受信ビームフォーミング機能を備えている。つまり、送受信部４６は、アレイ振動子を構成する各振動素子に対して、それぞれの遅延時間を設定して超音波ビームを形成する送信ビームフォーミング機能を有し、受信信号の信号受信タイミングを調整し、複数の受信信号の位相を整合させる整相加算機能を含む受信ビームフォーミング機能を有している。なお、送受信部４６は、受信信号の強度を補正するダイナミックレンジ変換回路などの信号処理回路を備えている。送受信部４６から出力されるエコーデータ列は、三次元画像処理部４８に伝送される。ちなみに、送受信部４６から出力されるエコーデータ列は、三次元画像処理部４８に出力するのと同時に、図示されていないカラードプラ処理部に対して出力してもよい。
【００３３】
三次元画像処理部４８は、複数の三次元画像を形成する機能を有し、本実施形態においては、反転処理部５０、第１レンダリング開始点判定部５２、第１レンダリング処理部５４、第２レンダリング開始点判定部５６、第２レンダリング処理部５８を備えている。ちなみに、本実施形態においては、ボリュームレンダリングの視線であるレイは、超音波ビームの進行方向と同一の方向に設定される。従って、三次元画像処理部４８に入力されたエコーデータ列は、そのままボクセル演算の対象であるボクセルデータ列とみなして処理することができる。
【００３４】
反転処理部５０は、入力されるデータの輝度を反転して出力する機能を有する。反転処理部５０においては、エコーデータ列を構成する各エコーデータの信号レベルを反転変換する。この反転変換を行うと低レベルのエコーデータは高レベルのエコーデータに置換処理される。故に、反転処理部５０は、超音波に対する反射強度の弱い生体組織、つまり前腕部においては血液の流れる部分からのエコー信号の検出を容易にする手段となる。
【００３５】
第１レンダリング開始点判定部５２は、第１レンダリング処理部５４で行われるボクセル演算の第１開始点を判定する機能を有する。その第１開始点の判定には、反転処理部５０から出力される反転エコーデータ列が用いられる。第１開始点の判定条件は、例えば、当該超音波診断装置のオペレータが入力部４２を操作して行うＲＯＩの領域設定を利用することにより定められる。ちなみに、第１開始点の判定条件としては、特開２００２−３３６２５３号公報に記されている被検体表面とカップリング体との境界面を検出する方法を用いることもできる。
【００３６】
第１レンダリング処理部５４は、前述の第１レンダリング開始点判定部５２で第１開始点が判定されたことを信号処理の起点としてボクセル演算の逐次演算を行い、あるレイ上の三次元画像の画素値を求める機能を有する。ボクセル演算については、例えば特開平１０−３３５３８号公報に記されている公知のボクセル演算法を用いることができる。
【００３７】
ボクセル演算法の一例としては、ボクセルの不透明度とエコー値との積によって表される発光光量と、ボクセルの透明度とそのボクセルの入力光との積によって表される減衰光量とを加算することにより、当該ボクセルからの出力光量を求める演算が挙げられる。ボクセルへの入力光量は、そのボクセルに近接する一つ前のボクセルの出力光量と同値であるため、逐次演算を行うことによってレイ方向の画素の輝度値を求めることができる。但し、演算対象のボクセルがレイ上の最終ボクセルに到達したり、レイ上の複数のボクセルのオパシティ積算値が１になったりした場合には、逐次演算を終了するという終了条件が付加される。
【００３８】
第２レンダリング開始点判定部５６は、第２レンダリング処理部５８で行われるボクセル演算の第２開始点を判定する機能を有する。その第２開始点の判定には、本実施形態においては、第１レンダリング処理部５４で行われるボクセル演算の逐次実行結果が利用される。すなわち、第１レンダリング処理部５４で行われるボクセル演算による透過光量の値が、所定の閾値を越える場合に、そのボクセル演算の演算対象であるエコーデータを第２開始点とする。
【００３９】
本実施形態にかかる超音波診断装置は、第１レンダリング処理の結果に基づいて、第２レンダリング処理を実行することを特徴事項の一つとしており、その第２レンダリング処理は第２開始点から開始される。第２開始点は、器具を主体とした第２三次元画像を取得するためのレンダリング処理の再スタートポイントあるいはリセットポイントになる。血管内に挿入された器具を描出する第２三次元画像を形成するためには、血管壁の部分、更に望ましくは血管内を流れる血流部分に入り込んだ部分に第２開始点を設定する必要がある。
【００４０】
血管壁と血液のそれぞれのエコー強度が大きく変化するという物理現象を利用して、血管の内部に入り込む位置を第２開始点として検出することができる。血管の内部に入り込んだ位置である第２開始点から、超音波画像を形成するレンダリング処理を変えることは、血管内に存在する血液以外の異物を視覚化する上で非常に有効な方法であり、特に、三次元画像化の場合には、三次元空間内での器具の位置を認識するために有効である。
【００４１】
第２開始点の判定方法には、前述したような透過光量の値が所定の閾値を越えることを利用する方法以外にも、いくつかの他の判定方法を用いることもできる。
【００４２】
第２開始点の他の判定方法としては、例えば、ボクセル演算の逐次演算と同時に並列実行されるオパシティの積算演算を用いることができる。すなわち、ボクセル演算によるオパシティの積算値が所定の閾値を超える場合には、そのボクセル演算の演算対象であるエコーデータを第２開始点とすることができる。
【００４３】
また、第２開始点の他の判定方法としては、前述の特開２００２−３３６２５３号公報に記されている、被検体表面とカップリング体との境界面を検出する方法と同様の手法を用いて、第２開始点を判定することができる。その同様の手法とは、すなわち関心領域を被検体表面を包含しない位置に予め設定しておき、その被検体表面を包含しない関心領域内におけるデータを輝度反転させた後に第１レンダリング処理部５４にて処理する手法である。高強度の反転エコーデータが連続してＮ個（一例として５個）続いた場合に、そのＮ個目のエコーデータを第２開始点とすることができる。
【００４４】
更に、第２開始点の他の判定方法としては、連続するデータ列の微分値を用いることもできる。すなわち、ボクセル演算による画像再構築を行う場合にレイの方向に沿って、連続するデータ列の構成要素である複数のデータを対象に微分を行い、微分値が急激に変化するデータに相当する位置を第２開始点とすることができる。
【００４５】
以上に記したいくつかの検出方法のいずれかの方法を用いて、第２開始点を検出することができる。
【００４６】
第２レンダリング処理部５８は、前述の第２レンダリング開始点判定部５６で判定された第２開始点からボクセル演算の逐次演算を行って、あるレイ上の三次元画像の画素値を求める機能を有する。第２レンダリング処理部５８には、反転処理される前のそのままのエコーデータが入力される。第２レンダリング処理部５８で行われるボクセル演算は、第１レンダリング処理部５４で行われるボクセル演算と同じ方法であってよいし、そのボクセル演算とは異なる三次元画像の構築演算であってもよい。
【００４７】
次に、三次元画像合成部６０について説明する。三次元画像合成部６０は、２種類の三次元画像を重ね合わせて、１つの合成三次元画像を形成する機能を有する。三次元画像合成部６０は、第１三次元画像メモリ６２、第２三次元画像メモリ６４などの構成を備えている。第１三次元画像メモリ６２は、第１レンダリング処理部５４で求められる複数の画素値をマッピングして構成される第１三次元画像を記憶するメモリである。第２三次元画像メモリ６４は、第２レンダリング処理部５８で求められる複数の画素値をマッピングして構成される第２三次元画像を記憶するメモリである。これらのメモリの中の第１三次元画像と第２三次元画像とを合成することにより、１つの合成三次元画像が形成される。
【００４８】
表示部６８は、ＣＲＴやＬＣＤで構成されており、合成三次元画像をオペレータが観察できるように表示する。
【００４９】
図３は、三次元画像処理部４８で行われる第１及び第２レンダリング処理の各演算処理に関する説明図である。図３には、可動部２０を前腕部１０に当接させて得られる三次元画像取り込み領域が模式的に示してある。
【００５０】
可動部２０を移動させることで互いに平行な複数の走査面が形成される。最前面の走査面７０上には、血管１４の断面と血管に挿入された器具１６が含まれる。また、走査面７０上には血管１４を中心にした三次元画像を形成するための関心領域（ＲＯＩ）８０の設定範囲が示されている。また、走査面７０上で可動部２０に近接した部分には、前腕部１０の体表面と水袋１８とが当接する境界部分８２が示されている。ちなみに、ＲＯＩ８０の設定形状は、Ｙ軸方向に伸長した三次元の直方体形状である。
【００５１】
図３においては、１本の超音波ビームに応じて取得されるボクセルデータ列を、走査面上のＺ軸に平行な１本の直線として表している。例えば、ボクセルデータ列Ａ７２は、血管１４が存在しない位置で取得されるボクセルデータ列である。
【００５２】
ボクセルデータ列Ａ７２は、ＲＯＩ８０の設定領域との関係において第１開始点が判定される。その第１開始点の位置は、図３に示す走査面７０においてボクセルデータ列Ａ７２を示す直線とＲＯＩ８０の上端面とが交差する点Ｐ１と定めることができる。
【００５３】
ボクセルデータ列Ａ７２に対するボクセル演算は、第１開始点である点Ｐ１の位置から、ＲＯＩ８０の下端面と交差する点Ｐ２までの範囲のデータ列に関して演算される。累積的にボクセル演算を実行している途中において、演算の終了条件を満たした場合には演算が中止され、そのときの演算結果がボクセルデータ列Ａ７２に関する画素値となる。または、ＲＯＩ８０の下端面上の点Ｐ２に至っても演算の終了条件が最後まで満たさない場合には、最終ボクセルでの演算結果が画素値となる。つまり、第１のレンダリング処理を行うにあたって、演算の終了条件を満足する／しないに関わらず、１つのボクセルデータ列に対して１つの画素値が求められる。その画素値は、図３に示すビーム対応画素値Ａ（８６）となる。第１のレンダリング処理を実行している最中に、第２開始点の判定条件を満足する点が検出されれば、第２のレンダリング処理に切り替わるが、ボクセルデータ列Ａ７２は、血管が存在しない部分で取得されるボクセルデータ列であるので、第２開始点は検出されない。
【００５４】
次に、血管１４と器具１６とが存在する位置で取得されるボクセルデータ列Ｂ８８に関する演算処理について説明する。ボクセルデータ列Ｂ８８は、前述のボクセルデータ列Ａ７２と同様に、ＲＯＩ８０の設定領域との関係において第１開始点が判定される。その第１開始点の位置は、図３に示す走査面７０においてボクセルデータ列Ｂ８８を示す直線とＲＯＩ８０の上端面とが交差する点Ｐ３と定めることができる。ボクセルデータ列Ｂ８８に対しては、第１開始点であるＰ３が判定された後から、第１レンダリング処理が行われる。
【００５５】
第１レンダリング処理を順に実行していくと、血管１４の内側に入った位置におけるボクセルデータの演算が行われる。血液からは弱いエコーデータが得られるが、反転処理によって高いエコーデータに変換されている。従って、第１レンダリング処理の出力光量値がある一定の閾値を越えるという判定条件を用いて、この条件を満たす点が判定される。なお、第１レンダリング処理の出力光量値がある一定の閾値を越えたという判定条件を満たした時点をもって、第１レンダリング処理を中止することも可能である。その場合には、第１レンダリング処理の終了点である第１終了点Ｐ４は、前述の第２開始点Ｐ５と同一の点となる。第１開始点Ｐ３から第１終了点Ｐ４までの間のボクセルデータに対して、第１レンダリング処理を行うことによって、図３に示すビーム対応画素値Ｂ１（９０）が求められる。
【００５６】
次に、第１レンダリング処理の第１終了点Ｐ４に基づいて、第２レンダリング処理の開始点である第２開始点Ｐ５が設定される。第２開始点Ｐ５は、血管１４の内側に入った位置に設定される。第２開始点Ｐ５の設定に伴って第２レンダリング処理が実行される。第２レンダリング処理では、反転されていないボクセルデータ列Ｂ８８が演算処理の対象となる。第２レンダリング処理は、血管内の器具を画像化するための処理である。器具の画像化のためには、器具を反映した高レベルのエコーデータにレンダリングパラメータが設定される。
【００５７】
第２レンダリング処理を順に実行していくと、器具１６の表面から強いレベルのエコーデータが検出される。ボクセルデータ列Ｂ８８の上であって、器具１６の表面に交差する位置を第２レンダリング処理の終了点Ｐ６として設定することができる。第２開始点Ｐ５から第２終了点Ｐ６までの間のボクセルデータに対して、第２レンダリング処理を行うことによって、図３に示すビーム対応画素値Ｂ２（９２）が求められる。
【００５８】
なお、一般的には、器具の表面からは高レベルのエコー信号が得られるので、器具１６の表面上において第２レンダリング処理が終了するように設定することができる。但し、ボリュームレンダリングのパラメータ設定の如何によっては、更に深さ方向に存在するボクセルデータの演算を累積することもできる。
【００５９】
以上説明したように、２段階のレンダリング処理を行うことにより、第１レンダリング処理でビーム対応画素値Ａ（８６）とビーム対応画素値Ｂ１（９０）が求められる。そして、第２レンダリング処理によりビーム対応画素値Ｂ２（９２）が求められる。例示した２つのボクセルデータ列Ａ７２とボクセルデータ列Ｂ８８のみならず、ＲＯＩ８０内に存在する全てのボクセルデータ列に対して第１レンダリング処理を行うことにより、第１三次元画像９４が形成される。また、第２開始点が判定されたボクセルデータ列に対して第２レンダリング処理を行うことにより、第２三次元画像９６が形成される。
【００６０】
第１三次元画像９４は、血管が存在しない生体組織に相当するボクセルデータ列と、血管が存在する部分については血管よりも手前側の生体組織に相当するボクセルデータ列に基づいて形成される。つまり、第１三次元画像９４は血管の内部の情報を含まず、生体組織を主体的に表した三次元投影画像となる。また、第２三次元画像９６は、血管あるいはその血管の裏側に位置する生体組織に相当するボクセルデータ列に基づいて形成される。よって、第２三次元画像９６では、器具投影部分９８を認識することができる。
【００６１】
図４は、図３に示す走査面７０内において、第１レンダリング処理と第２レンダリング処理が適応される領域を区別して示した図である。図４には、図３に示す走査面７０内のＲＯＩ８０の領域だけが示してある。
【００６２】
ＲＯＩ８０の上端には、第１レンダリング処理が開始される第１開始点群１６２が直線として示してある。前述の点Ｐ３も、第１開始点群１６２の上に存在する。また、ＲＯＩ８０の下端には、第１終了点群（１６４Ａ、１６４Ｂ）が２つの線分として示してある。血管１４が存在する部分は、第１レンダリング処理が血管内で終了するようにレンダリング条件が設定されているので、血管に対する第１終了点群１６６は血管壁の形状に対応した円弧形状となる。第１のレンダリング処理が施される領域は、第１開始点群１６２と第１終了点群（１６４Ａ，１６４Ｂ，１６６）とに挟まれる領域であるので、図５において図示する逆凹型の斜線領域１６０となる。
【００６３】
本実施形態においては、第２レンダリング処理は、血管が存在する位置から得られるボクセルデータ列に対して実行される。よって、血管内の第１終了点群である１６６が第２開始点群１６８となる。第２レンダリング処理は、レイ１７０の深さ方向に沿って、第２開始点群１６８から開始されるので、第２レンダリング処理が施される領域は、図５において図示する斜線領域１７２となる。ちなみに、器具１６の表面の位置で第２レンダリング処理を中止することもできる。
【００６４】
図５は、ある一つのエコーデータ列を信号処理することで生成する各種信号を並べて示した図である。図５の最上段に示す図５（ａ）は、エコーデータ列の強度レベルを示している。図５（ａ）の横軸は時間軸であり、縦軸はエコーデータの強度レベルを示す。このエコーデータ列は、血管内に存在する器具に向かって１本の超音波ビームを送信することで得られたエコーデータ列を表している。図３に示したボクセルデータ列Ｂ８８がこのエコーデータ列に相当する。
【００６５】
エコーデータ列１００は、そのデータ列の中間部分１０２において、レベルが急激に低下したり（１０４ａ，１０４ｂ）、一時的にピークを示したり（１０６）している。レベルが急激に低下しているのは、この部分のエコーデータが血管内の低いエコー強度を反映しているためであり、レベルが一時的にピークを示すのは、このピーク部分のエコーデータが器具によるエコー強度を反映しているためである。なお、中間部分１０２に含まれず時間軸上で前後に広がるエコーデータは、生体組織によって得られるエコーデータである。
【００６６】
図５（ｂ）は、エコーデータ列１００を反転処理して得られた反転エコーデータ列１０８を示すグラフである。図５（ｂ）の横軸と縦軸の定義は図５（ａ）と同一である。反転処理を行うことにより、エコーデータ列１００において低レベルである部分（１０４ａ，１０４ｂ）が、高レベルのデータに変換される。
【００６７】
ここで、図５の説明を一時的に離れて、図６を用いてエコーデータの反転処理特性の一例を示す。図６は、エコーデータの反転処理特性の一例を示す図である。横軸はエコー信号入力レベルを示し、縦軸はエコー信号出力レベルを示している。この曲線１２４に従った変換処理を行うと、低レベルの入力信号は高レベルの出力信号として変換され、逆に高レベルの入力信号は低レベルの出力信号として変換される。なお、グラフは非線形特性を示しているので、入力信号のレベルが低い程、高いレベルに強調された反転信号が得られるようになっている。また、ある閾値を上回る入力信号については、出力信号として変換されないカットオフ特性も有している。
【００６８】
図５に戻って、反転エコーデータ列１０８の信号処理について説明する。本実施形態に係る超音波診断装置においては、この反転エコーデータ列１０８について、第１レンダリング処理を施している。第１レンダリング処理部５４でボクセル演算を実行するのに伴って、各ボクセルデータからの出力光量値が求められ、それと平行して同時に不透明度（オパシティ）の積算値も求められる。
【００６９】
図５（ｃ）は第１レンダリング処理によるオパシティの積算値（Σαi）の累積的な変
化を示し、図５（ｄ）は第１レンダリング処理による出力光量値（Ｃ１out）の累積的な変化を示している。本実施形態においては、出力光量値（Ｃ１out）の値が一定の閾値を超えた場合に、第２レンダリング処理が開始される。よって、図５（ｄ）において出力光量値を表す曲線１１２と、閾値Ｃthを表す水平線との交差点Ｔが存在することが、第２レンダリング処理の開始を意味する。本実施形態においては、第１レンダリング処理は、第２レンダリング処理が開始されたことをもって、その演算処理を中断する。但し、第２レンダリング処理の開始に関わらず、第１レンダリング処理の演算を終了条件を満たすまで継続してもよい。第１レンダリング処理は第２レンダリング処理とは別個独立に実行されるからである。図５（ｄ）においては、第１レンダリング処理の中断又は継続が任意であることを破線によって示している。ボクセル演算はオパシティの積算演算と表裏一体に実行されるので、図５（ｃ）においてもその積算値を表す曲線の一部を破線で示している。
【００７０】
図５（ｅ）は、第２レンダリング処理の処理対象となるエコーデータ列１１８を示している。このエコーデータ列１１８は、図５（ａ）に示すエコーデータ列１００の後半部、つまり時間軸上にて一点鎖線１１６で示す時点以降で得られたエコーデータ列と同値である。第２レンダリング処理に関して注目すべきことは、処理対象のデータ列は反転エコーデータ列ではないことと、第２レンダリング処理のための第２のオパシティが新たにゼロから積算されることである。第２レンダリング処理を行う場合に、反転エコーデータ列を用いると、最初から高レベルのデータを処理することになり、器具の画像化に適さない。第２レンダリング処理は、ある判定条件を満たした場合に開始されるので、判定条件を満足した時点からオパシティの積算演算が新たに開始される。第２レンダリング処理を行うことによって、エコーデータ列１１８の中でピークを示している器具のエコーデータ１２０に基づいて器具に対応する画素値が求められる。
【００７１】
図５（ｆ）には、第２レンダリング処理による出力光量値Ｃ２out（１２２）の変化が示されている。第２レンダリング処理では、器具を明瞭に画像化することを優先してレンダリングパラメータを設定することが好ましい。従って、この出力光量値Ｃ２outについても、器具または血管内に相当する時間的な範囲での出力光量値の変化は大きいが、それ以降は飽和傾向になるため時間的な変化は少なくなる。
【００７２】
図７は、他の実施形態に係る超音波診断装置の構成ブロック図である。図２の構成ブロック図と比較して、三次元画像処理部１３０の構成が異なっている。よって、図７において、図２に示した構成と同様の構成には同一符号を付しその説明を省略する。
【００７３】
三次元画像処理部１３０においては、第２レンダリング開始点判定部５６に対して、反転されていないエコーデータ列が入力される。第２レンダリング開始点判定部５６では、低レベルのエコーデータが連続してＭ個（一例として５個）続いた場合に、そのＭ個目のエコーデータを第２開始点とする判定処理が行われる。低レベルのエコーデータが検出された後に、引き続いてＭ個のデータが連続しなければ判定条件が成立しないので、血管内に確実に入り込んだ位置に第２開始点を設定することができる。判定条件を満たす第２開始点が検出されると、その情報は、第１レンダリング処理部５４に伝達される。この情報伝達の方向は、図２に示した機能ブロック図で示した情報伝達の方向と逆である。図７に示す第１レンダリング処理部５４においては、第２開始点の検出をもって、その検出の直前まで実行していたボクセル演算を中止する。一方、本実施形態によると、第２開始点が検出されることをもって、生体組織ではなく血管内でのレンダリング処理に移行していることを確実に判断できる。よって、三次元画像処理部１３０によれば、第２開始点よりも更に深い位置から得られるエコーデータの情報に基づいて、第１レンダリング処理が実質的に処理不要になる位置を検出することにより、第１レンダリング処理を強制的に中止することができる。
【００７４】
また、他の実施形態に係る超音波診断装置においては、図７に示す超音波診断装置の装置本体４０が備えているカラードプラ処理部（図示せず）を用いて、以下に示す手法により第２開始点を検出することもできる。
【００７５】
カラードプラ処理部は、送受信部４６が出力する受信信号を入力して、その分析結果を第２レンダリング開始点判定部５６に対して出力する。カラードプラ処理部においては、受信信号を直交検波した後にＭＴＩフィルタにより低いドプラ偏移周波数を除去し、自己相関演算することでドプラ偏移周波数の平均周波数とパワーを検出することができる。このようにドプラ偏移周波数の平均周波数とパワーに基づいて血流のイメージング表示を行う技術は、いわゆるカラーフロー又はパワーフローと呼ばれる公知の技術である。この公知技術を用いれば、超音波断層画像上で血流の存在する部分と存在しない部分の境界を判定することができる。その判定結果、つまりカラードプラ処理部で求めた各レイ上での血流の存在位置を示す境界情報を第２レンダリング開始点判定部５６に伝送することによって、第２レンダリング開始点判定部５６において第２開始点を定めることができる。この判定方法によれば、血流の存在する血管内部の位置に確実に第２開始点を設定することができる利点がある。
【００７６】
カラードプラ処理部にて行われる処理は、カラーフロー処理よりもパワーフロー処理の方が好ましい。パワーフロー処理は、カラーフロー処理に比べてドプラアングルの依存性が少ないため、血管の走行を連続的に捉えられる。なお、一般的に、器具が血管内を移動する速度は血液が血管内を流れる速度に比べて無視できる程度に小さいので、器具の運動によってドプラ偏移周波数が影響を受けることはない。
【００７７】
図８は、他の実施形態に係る三次元画像合成部の構成ブロック図である。図８に示す三次元画像合成部１４０は、図２に示す三次元画像合成部６０の変形例である。三次元画像合成部１４０は、三次元画像合成部６０と置換して用いることが可能であり、２枚の三次元画像を合成する際に、各画像に輝度の重み付けを付加する機能を有している。三次元画像合成部１４０は、第１表示処理部１４２と第２表示処理部１４４を備える点において、前述の三次元画像合成部６０と相違する。第１表示処理部１４２は、第１三次元画像の輝度と色彩を変換あるいは調整する機能を有しており、その機能を使用するために、例えば輝度変換係数Ｊ（設定可能範囲：−１≦Ｊ≦１）のパラメータが活用される。また、第２表示処理部１４４は、第２三次元画像の輝度と色彩を変換あるいは調整する機能を有しており、その機能を使用するために、例えば輝度変換係数Ｋ（設定可能範囲：−１≦Ｋ≦１）のパラメータが活用される。ちなみに、輝度変換係数Ｊ及びＫは、入力部４２からのオペレータの操作入力信号に基づいて、主制御部４４から出力され、三次元画像合成部１４０に対して入力される係数である。
【００７８】
輝度変換係数Ｊ及びＫのそれぞれの機能は同一である。すなわち、輝度変換係数が１（Ｊ＝１，Ｋ＝１）である場合には三次元画像の輝度調整を一切行わず、輝度変換係数が−１（Ｊ＝−１，Ｋ＝−１）の場合には輝度反転処理を行う。そして、輝度変換係数が１から−１の間の小数値が指令された場合には輝度のグラデーション調整が行われる。
【００７９】
このような輝度反転及び調整機能を有することにより、第１三次元画像と第２三次元画像とは、それぞれが独立に輝度調整を行うことができる。よって、血流を明瞭に画像化するためにレンダリング処理前にエコーデータの反転処理が行われる場合であっても、合成三次元画像においては生体組織と器具とを明暗により、あるいは色彩の違いによって識別することが可能となり、双方とも明瞭に三次元画像として形成することができる。
【００８０】
図９は、三次元画像合成部１４０を用いて実行される画像合成の説明図である。図２にて説明した三次元画像合成部６０、あるいは図８を用いて説明した三次元画像合成部１４０のどちらにおいても、第１三次元画像と第２三次元画像との画像合成が行われる。画像合成は、第１三次元画像の上に第２三次元画像を重ね合わせる、いわゆるオーバーレイ処理によって実行される。図９に示す第１三次元画像１４６には、生体組織内で屈曲した状態の血管１５２が示されている。この第１三次元画像１４６は、三次元のレンダリング処理を施して形成された三次元投影画像である。血管１５２の部分が白く表示されているのは、輝度反転したエコーデータ群によって第１三次元画像１４６が形成されるためである。第２三次元画像１４８には、血管に挿入された器具１５４の形状が示されている。この第２三次元画像１４８も三次元投影画像である。器具１５４の部分が白く表示されているのは、輝度反転していないエコーデータ群によって第２三次元画像１４８が形成されるためである。図９に示した状態のまま、血管１５２と器具１５４の双方が白く表示された状態で２枚の三次元画像を合成処理すると、血管内の器具１５４の位置を識別することが難しくなる。識別を容易にするために、図８を用いて説明した輝度変換係数ＪとＫをＪ＝１，Ｋ＝−１と設定することにより、第２三次元画像１４８だけについて反転処理が実行される。そうすると、器具１５４の部分が黒く変換されるため、合成三次元画像１５０においては、血管１５２の立体形状が明確に表示されると同時に、血管内の器具１５６の位置も明確に把握することができる。なお、第２三次元画像に対して着色表示を施せば、器具１５４の部分を生体組織と異なる色に表示することが可能となり、より明瞭に両者を識別することができる。
【図面の簡単な説明】
【００８１】
【図１】前腕部に超音波の送受信器を当接している状態を模式的に示した図である。
【図２】本発明に係る超音波診断装置の概略的な全体構成を示すブロック図である。
【図３】三次元画像処理部で行われる第１及び第２レンダリング処理の各演算処理に関する説明図である。
【図４】第１及び第２レンダリング処理の領域の区別を視覚的に示した図である。
【図５】ある一つのエコーデータ列の２信号処理についての説明図である。
【図６】エコーデータの反転処理特性の一例を示す図である。
【図７】他の実施形態に係る超音波診断装置の構成ブロック図である。
【図８】他の実施形態に係る三次元画像合成部の構成ブロック図である。
【図９】三次元画像合成部において実行される画像合成の説明図である。
【符号の説明】
【００８２】
１４ 血管、１６ 器具、２０ 可動部、７０ 走査面、７２ ボクセルデータＡ、８０ 関心領域（ＲＯＩ）、８２ 境界部分、８６ ビーム対応画素値Ａ、９０ ビーム対応画素値Ｂ１、９２ ビーム対応画素値Ｂ２、９４ 第１三次元画像、 ９６ 第２三次元画像、９８ 器具投影部分、Ｐ１，Ｐ３ 第１開始点、Ｐ２，Ｐ４ 第１終了点、Ｐ５ 第２開始点、Ｐ６ 第２終了点。
【特許請求の範囲】
【請求項１】
器具が挿入された血管に超音波を送受波して得られたボクセルデータ群に対して設定された各レイについて、深さ方向にボクセル演算を逐次実行する超音波診断装置において、
前記各レイについて、第１開始点から前記ボクセル演算の逐次実行結果により、第１三次元画像を構成する画素値を求める第１レンダリング手段と、
前記各レイについて第２開始点を判定する判定手段と、
複数のレイの中で前記第２開始点が判定されたレイについて、前記第２開始点からのボクセル演算の逐次実行結果により、前記器具が反映された第２三次元画像を構成する画素値を求める第２レンダリング手段と、
前記第１三次元画像と前記第２三次元画像とを合成し合成三次元画像を形成する合成処理手段と、
前記合成三次元画像を表示する表示手段と、
を有することを特徴とする超音波診断装置。
【請求項２】
請求項１記載の超音波診断装置において、
前記判定手段は、前記血管内かつ前記器具の手前側の深さ位置を前記第２開始点として判定することを特徴とする超音波診断装置。
【請求項３】
請求項１記載の超音波診断装置において、
前記第２開始点は、前記第１開始点からのボクセル演算の終了条件を満たす第１終了点に基づいて、前記血管内に設定されることを特徴とする超音波診断装置。
【請求項４】
請求項１記載の超音波診断装置において、
前記判定手段は、前記各レイについての前記第１開始点からのボクセル演算の逐次実行結果に基づいて、前記第２開始点の判定を行うことを特徴とする超音波診断装置。
【請求項５】
請求項１記載の超音波診断装置において、
前記判定手段は、前記各レイについて逐次実行される前記第１開始点からのボクセル演算で用いられるパラメータの積算値に基づいて、前記第２開始点の判定を行うことを特徴とする超音波診断装置。
【請求項６】
請求項１記載の超音波診断装置において、
前記各レイ上のボクセルデータ列に対して輝度反転処理を施す前処理手段を含み、
前記判定手段、前記第１レンダリング手段及び前記第２レンダリング手段の内の少なくとも１つは前記輝度反転処理が施されたボクセルデータ列に基づいて処理を実行することを特徴とする超音波診断装置。
【請求項７】
請求項１記載の超音波診断装置において、
超音波を送受波して得られた受信信号を直交検波し、血液の移動に伴って発生するドプラ偏移周波数成分に応じて血流の存在位置を検出するカラードプラ処理手段を更に含み、
前記判定手段は、前記カラードプラ処理手段が検出する血流位置情報に基づいて、前記第２開始点の判定を行うことを特徴とする超音波診断装置。
【請求項８】
請求項１記載の超音波診断装置において、
前記合成処理手段は、前記血管を主体的に表した第１三次元画像上に、前記器具を主体的に表した第２三次元画像を重ね合わせて前記合成三次元画像を形成することを特徴とする超音波診断装置。
【請求項９】
請求項１記載の超音波診断装置において、
前記合成処理手段は、前記第１三次元画像と前記第２三次元画像とを視覚的に異ならせつつ合成処理を行うことを特徴とする超音波診断装置。
【請求項１】
器具が挿入された血管に超音波を送受波して得られたボクセルデータ群に対して設定された各レイについて、深さ方向にボクセル演算を逐次実行する超音波診断装置において、
前記各レイについて、第１開始点から前記ボクセル演算の逐次実行結果により、第１三次元画像を構成する画素値を求める第１レンダリング手段と、
前記各レイについて第２開始点を判定する判定手段と、
複数のレイの中で前記第２開始点が判定されたレイについて、前記第２開始点からのボクセル演算の逐次実行結果により、前記器具が反映された第２三次元画像を構成する画素値を求める第２レンダリング手段と、
前記第１三次元画像と前記第２三次元画像とを合成し合成三次元画像を形成する合成処理手段と、
前記合成三次元画像を表示する表示手段と、
を有することを特徴とする超音波診断装置。
【請求項２】
請求項１記載の超音波診断装置において、
前記判定手段は、前記血管内かつ前記器具の手前側の深さ位置を前記第２開始点として判定することを特徴とする超音波診断装置。
【請求項３】
請求項１記載の超音波診断装置において、
前記第２開始点は、前記第１開始点からのボクセル演算の終了条件を満たす第１終了点に基づいて、前記血管内に設定されることを特徴とする超音波診断装置。
【請求項４】
請求項１記載の超音波診断装置において、
前記判定手段は、前記各レイについての前記第１開始点からのボクセル演算の逐次実行結果に基づいて、前記第２開始点の判定を行うことを特徴とする超音波診断装置。
【請求項５】
請求項１記載の超音波診断装置において、
前記判定手段は、前記各レイについて逐次実行される前記第１開始点からのボクセル演算で用いられるパラメータの積算値に基づいて、前記第２開始点の判定を行うことを特徴とする超音波診断装置。
【請求項６】
請求項１記載の超音波診断装置において、
前記各レイ上のボクセルデータ列に対して輝度反転処理を施す前処理手段を含み、
前記判定手段、前記第１レンダリング手段及び前記第２レンダリング手段の内の少なくとも１つは前記輝度反転処理が施されたボクセルデータ列に基づいて処理を実行することを特徴とする超音波診断装置。
【請求項７】
請求項１記載の超音波診断装置において、
超音波を送受波して得られた受信信号を直交検波し、血液の移動に伴って発生するドプラ偏移周波数成分に応じて血流の存在位置を検出するカラードプラ処理手段を更に含み、
前記判定手段は、前記カラードプラ処理手段が検出する血流位置情報に基づいて、前記第２開始点の判定を行うことを特徴とする超音波診断装置。
【請求項８】
請求項１記載の超音波診断装置において、
前記合成処理手段は、前記血管を主体的に表した第１三次元画像上に、前記器具を主体的に表した第２三次元画像を重ね合わせて前記合成三次元画像を形成することを特徴とする超音波診断装置。
【請求項９】
請求項１記載の超音波診断装置において、
前記合成処理手段は、前記第１三次元画像と前記第２三次元画像とを視覚的に異ならせつつ合成処理を行うことを特徴とする超音波診断装置。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【公開番号】特開２００８−７９８２１（Ｐ２００８−７９８２１Ａ）
【公開日】平成２０年４月１０日（２００８．４．１０）
【国際特許分類】
【出願番号】特願２００６−２６３１６６（Ｐ２００６−２６３１６６）
【出願日】平成１８年９月２７日（２００６．９．２７）
【出願人】（３９００２９７９１）アロカ株式会社 (899)
【Ｆターム（参考）】
【公開日】平成２０年４月１０日（２００８．４．１０）
【国際特許分類】
【出願日】平成１８年９月２７日（２００６．９．２７）
【出願人】（３９００２９７９１）アロカ株式会社 (899)
【Ｆターム（参考）】
[ Back to top ]