(57)【特許請求の範囲】
【請求項１】
被検体における同一領域を少なくとも一部重複させるように超音波スキャンを複数回行い、受信エコーに基づいて複数フレームの画像を生成する超音波診断装置であって、
前記複数フレームの画像のそれぞれに対して、所定の変換処理により多重解像度解析を行う解析手段と、
前記解析手段による多重解像度解析により取得した各解像度の各係数毎に、前記複数フレームの画像間の対応する係数同士でそれぞれフィルタ演算処理を行うフィルタ演算手段と、
前記フィルタ演算手段によるフィルタ演算処理結果に対して、前記所定の変換処理の逆変換処理を施して１フレームの画像を生成する画像生成手段と、
を具備することを特徴とする超音波診断装置。
【請求項２】
前記解析手段により行われる前記所定の変換処理とは、ウェーブレット変換、ステーショナリウェーブレット変換、及びピラミッド変換のうちの何れかによる変換処理であること特徴とする請求項１に記載の超音波診断装置。
【請求項３】
前記フィルタ演算処理手段による前記フィルタ演算処理とは、前記多重解像度解析により取得した前記係数のうち低周波側の係数については平均値を取る演算であり、高周波側の係数については絶対値の大きい値を出力する演算処理であることを特徴とする請求項１に記載の超音波診断装置。
【請求項４】
前記フィルタ演算処理手段による前記フィルタ演算処理とは、全ての前記係数について、絶対値の大きい値を出力する演算処理であることを特徴とする請求項１に記載の超音波診断装置。
【請求項５】
前記フィルタ演算処理手段は、前記複数フレームの画像における所定の位置の係数の出力値を決定する際に、当該所定の位置の近傍の係数も利用して前記フィルタ演算処理を行うことを特徴とする請求項１に記載の超音波診断装置。
【請求項６】
前記解析手段による前記多重解像度解析を実施した後であって、前記フィルタ演算処理手段による前記フィルタ演算処理を実施する前に、前記複数フレームの画像における互いに重複する領域の辺縁部において、前記解析手段による前記所定の変換処理に起因して発生するエッジ成分を除去するエッジ成分除去手段を含むことを特徴とする請求項１に記載の超音波診断装置。
【請求項７】
前記フィルタ演算手段により前記フィルタ演算処理を行う前及び／又は後に、各々の前記画像毎に非線形拡散フィルタ処理を施す非線形拡散フィルタ手段を含むことを特徴とする請求項１に記載の超音波診断装置。
【請求項８】
前記画像生成手段による前記逆変換処理を行う前に、前記フィルタ演算処理後の各ウェーブレット係数に対して所定の係数を乗算する画像処理手段を含むことを特徴とする請求項２に記載の超音波診断装置。
【請求項９】
前記画像を所定の表示手段に表示させる為の表示座標へ座標変換する座標変換手段を含み、
前記解析手段、前記フィルタ演算処理手段、及び前記画像生成手段による処理は、前記座標変換手段による座標変換処理前又は座標変換処理後に行うことを特徴とする請求項１乃至請求項８のうち何れか一つに記載の超音波診断装置。
【請求項１０】
前記超音波スキャンは、
被検体の同一部位に対して複数の異なる方向から超音波ビームの送受信を行なうスキャン、
被検体上で超音波プローブを移動させながら画像を取得していき、これら取得した画像に基づいて前記超音波プローブの動きを検出して前記画像の位置合わせを行なうスキャン、
被検体における同一の領域に対して複数の異なる位置から３次元スキャンを実施するスキャン、
超音波ビームの送信周波数及び／又は受信周波数を異ならせて被検体に対して同一の方向から複数回の超音波ビームの送受信を行うスキャン、及び
被検体に対して同一の方向から送信フォーカスを変えて複数回の超音波ビームの送受信を行うスキャン、
を含むことを特徴とする請求項９に記載の超音波診断装置。
【請求項１１】
被検体における同一領域を少なくとも一部重複させるように超音波スキャンを複数回行い、受信エコーに基づいて複数フレームの画像を生成する超音波診断装置としてコンピュータを機能させるプログラムであって、
前記複数フレームの画像のそれぞれに対して、所定の変換処理により多重解像度解析を行う解析機能と、
前記解析機能による多重解像度解析により取得した各解像度の各係数毎に、前記複数フレームの画像間の対応する係数同士でそれぞれフィルタ演算処理を行うフィルタ演算機能と、
前記フィルタ演算機能によるフィルタ演算処理結果に対して、前記所定の変換処理の逆変換処理を施して１フレームの画像を生成する画像生成機能と、
をコンピュータに実現させることを特徴とするプログラム。
【請求項１２】
前記解析機能により行われる前記所定の変換処理とは、ウェーブレット変換、ステーショナリウェーブレット変換、及びピラミッド変換のうちの何れかによる変換処理であること特徴とする請求項１１に記載のプログラム。
【請求項１３】
前記フィルタ演算処理機能による前記フィルタ演算処理とは、前記多重解像度解析により取得した前記係数のうち低周波側の係数については平均値を取る演算であり、高周波側の係数については絶対値の大きい値を出力する演算処理であることを特徴とする請求項１１に記載のプログラム。
【請求項１４】
前記フィルタ演算処理機能による前記フィルタ演算処理とは、全ての前記係数について、絶対値の大きい値を出力する演算処理であることを特徴とする請求項１１に記載のプログラム。
【請求項１５】
前記フィルタ演算処理機能は、前記複数フレームの画像における所定の位置の係数の出力値を決定する際に、当該所定の位置の近傍の係数も利用して前記フィルタ演算処理を行う機能であることを特徴とする請求項１１に記載のプログラム。
【請求項１６】
前記解析機能による前記多重解像度解析を実施した後であって、前記フィルタ演算処理機能による前記フィルタ演算処理を実施する前に、前記複数フレームの画像における互いに重複する領域の辺縁部において、前記解析機能による前記所定の変換処理に起因して発生するエッジ成分を除去するエッジ成分除去機能を更にコンピュータに実現させることを特徴とする請求項１１に記載のプログラム。
【請求項１７】
前記フィルタ演算機能により前記フィルタ演算処理を行う前及び／又は後に、各々の前記画像毎に非線形拡散フィルタ処理を施す非線形拡散フィルタ機能を更にコンピュータに実現させることを特徴とする請求項１１に記載のプログラム。
【請求項１８】
前記画像生成機能による前記逆変換処理を行う前に、前記フィルタ演算処理後の各ウェーブレット係数に対して所定の係数を乗算する画像処理機能をコンピュータに更に実現させることを特徴とする請求項１２に記載のプログラム。
【請求項１９】
前記画像を所定の表示手段に表示させる為の表示座標へ座標変換する座標変換機能をコンピュータに実現させ、
前記解析機能、前記フィルタ演算処理機能、及び前記画像生成機能による処理は、前記座標変換機能による座標変換処理前又は座標変換処理後に行うことを特徴とする請求項１１乃至請求項１８のうち何れか一つに記載のプログラム。
【請求項２０】
前記超音波スキャンは、
被検体の同一部位に対して複数の異なる方向から超音波ビームの送受信を行なうスキャン、
被検体上で超音波プローブを移動させながら画像を取得していき、これら取得した画像に基づいて前記超音波プローブの動きを検出して前記画像の位置合わせを行なうスキャン、
被検体における同一の領域に対して複数の異なる位置から３次元スキャンを実施するスキャン、
超音波ビームの送信周波数及び／又は受信周波数を異ならせて被検体に対して同一の方向から複数回の超音波ビームの送受信を行うスキャン、及び
被検体に対して同一の方向から送信フォーカスを変えて複数回の超音波ビームの送受信を行うスキャン、
を含むことを特徴とする請求項１９に記載のプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波診断装置及びプログラムに関し、特に超音波スキャンして取得した複数フレームの画像を合成処理する超音波診断装置及びプログラムに関する。
【背景技術】
【０００２】
超音波診断装置は、超音波プローブに内蔵された圧電振動子から発生する超音波を被検体内に放射し、被検体組織の音響インピーダンスの差異によって生ずる反射信号を前記圧電振動子によって受信してモニタ上に表示するものである。この超音波診断装置によれば、超音波プローブを体表に接触させるだけの簡単な操作でリアルタイムの画像データを容易に得ることができる為、臓器の機能診断や形態診断に広く用いられている。
【０００３】
被検体の組織或いは血球からの反射波により生体情報を得る超音波診断法は、超音波パルス反射法と超音波ドプラ法の２つの大きな技術開発により急速な進歩を遂げ、前記技術を用いて得られるＢモード画像データとカラードプラ画像データは、今日の超音波診断において不可欠なものとなっている。
【０００４】
ところで、超音波診断装置によって得られる超音波画像データにおいては、被検体内の複数の反射体からの反射超音波のランダムな干渉によって生ずる粒状性のノイズ（所謂、スペックルノイズ）が発生し、このスペックルノイズは装置のシステムノイズと共に超音波診断の妨げとなることが知られている。
【０００５】
このようなスペックルノイズを低減する方法としては、例えばコンパウンドスキャン法を挙げることができる。コンパウンドスキャン法は更に空間コンパウンドスキャンと周波数コンパウンドスキャンに分類される。
【０００６】
前記空間コンパウンドスキャンは、被検体の同一部位に対して複数の異なる方向から超音波送受波を行ない、得られた複数の画像データを加算合成して表示用の画像データを生成する技術である。他方、前記周波数コンパウンドスキャンは、被検体の同一部位に対し異なる超音波周波数によって複数枚の画像データを収集し、これらの画像データを加算合成して表示用画像データを生成する技術である。
【０００７】
このように画像データを合成処理する技術として、例えば特許文献１には、次のような技術が開示されている。すなわち、特許文献１には、個々の連続するイメージフレームを複数のサブイメージ領域に分割し、該サブイメージ領域の局部移動ベクトルを推定評価する段階と、推定評価された局部移動ベクトルを基に包括的イメージ移動を推定評価する段階と、推定評価された包括的なイメージ移動を基に複合されたＸＦＯＶイメージを表示する段階とより成る、複合された超音波イメージを発生させる方法が開示されている。
【０００８】
この特許文献１に開示された方法によれば、拡大された視野（ＸＦＯＶ）イメージとして表される、大きな複合超音波イメージを発生させることができる。
【０００９】
また、前記スペックルノイズを低減する方法としては、例えば特許文献２に次のような技術が開示されている。すなわち、特許文献２には、被検体の複数方向に対して超音波送受波を行なって画像データを生成する画像データ生成手段と、前記画像データに対しウェーブレット変換を行なって各解像度レベルにおける低周波ウェーブレット変換係数と複数の高周波ウェーブレット変換係数を算出するウェーブレット変換手段と、前記高周波ウェーブレット変換係数に基づいて前記画像データのエッジ情報を検出するエッジ情報検出手段と、検出された前記エッジ情報に基づいて前記高周波ウェーブレット変換係数を重み付け演算する係数処理手段と、重み付け演算された前記高周波ウェーブレット変換係数を用いたウェーブレット逆変換によってエッジ強調された画像データを生成するウェーブレット逆変換手段と、前記エッジ強調された画像データを表示する表示手段を備えた超音波診断装置が開示されている。
【００１０】
この特許文献２に開示された超音波診断装置によれば、超音波画像データにおいて効果的なスペックルノイズの低減とエッジ強調が可能となる。
【００１１】
ところで、超音波診断装置で組織像を表示する際に、超音波スキャンして得られた画像データを合成して表示する際の合成方法としては、主として次の３種類の合成方法が知られている。
【００１２】
（合成方法１）空間コンパウンドスキャン
空間コンパウンドスキャンでは、図８に示すように、超音波プローブ１０１を固定して保持し、フレーム毎に異なる方向からスキャンを行う。このように異なる方向から超音波スキャンを行うと、同じ場所からのエコーであっても干渉条件が異なる為、スペックルの相関が少なくなる。従って、重複する領域（第１スキャン画像１０３と第２スキャン画像１０５との重複領域１０７）について平均値演算又は最大値演算を行うことにより、スペックルを低減させることができる。
【００１３】
また、超音波ビームに垂直な組織境界は鏡面反射により、より明瞭な境界エコーを発生する為、多数の方向から超音波ビームを照射することで、いろいろな角度の組織境界からのエコーを明瞭に表示させることが可能となる。
【００１４】
（合成方法２）パノラミックビュー
パノラミックビューとは、特許文献１で示されるような映像法であり、図９に示すように超音波プローブ１０１をスキャン断面の方向に動かし、得られた画像データから超音波プローブ１０１の動きを検出して、画像データを並行移動且つ回転移動して位置合わせを行い合成処理することで、より広い視野を表示する技術である。現在は最新の画像が上になるような上書きによる合成処理が一般に行われている。
【００１５】
（合成方法３）コンビネーションフォーカス
コンビネーションフォーカスとは、図１０に示すように、異なる方向からではなく同一の方向から、送信フォーカスを変えて複数回の超音波ビームの送受信を行い、例えば近距離の領域の画像１１１Ａに関しては近距離フォーカスの超音波ビーム１１１により取得し、遠距離の領域の画像１１３Ａに関しては遠距離フォーカスの超音波ビーム１１３により取得し、それらの画像を合成する技術である。
【特許文献１】特開平８−２８０６８８号公報
【特許文献２】特開２００５−２９６３３１号公報
【発明の開示】
【発明が解決しようとする課題】
【００１６】
しかしながら、上述した画像データの合成方法は、それぞれ次のような問題を抱えている。
【００１７】
まず、空間コンパウンドスキャンに関しては、画像データを合成する際に、図８に示すように当該合成に係る画像データ（第１スキャン画像１０３と第２スキャン画像１０５と）の間で、重複する領域１０７と重複しない領域１０９とが生じる。ここで、重複しない領域１０９に関しては、当該合成に係る画像データをそのまま表示することになる。他方、重複する領域１０７における画像合成方法としては、例えば平均値演算や最大値演算等を挙げることができる。
【００１８】
ここで、最大値演算を採る場合には、重複する領域１０７と重複しない領域１０９とのつなぎ目（境界）における表示は自然な表示となるが、一つの画像データに階調の高いアーティファクトがあると、該アーティファクトがそのまま表示されてしまう。一方、平均値演算を採る場合には、重複する領域１０７と重複しない領域１０９とのつなぎ目（境界）における表示が不自然に見えてしまう。更に、最大値演算を採る場合であっても、平均値演算を採る場合であっても、当該合成処理後の画像データは、どこかぼけた表示であるという印象を観察者に与える画像データとなるという問題がある。
【００１９】
また、パノラミックビューによる場合は、重複する領域が生じた場合には、上述したように当該領域の画像データを上書き処理を行う。従って、合成処理に起因する画像のぼけは生じない。しかしながら、重複する領域の画像データが存在するにも関わらず、それを利用しないということは、画像データを合成処理することによるメリットを充分に享受しているとは言い難い。換言すれば、重複する領域の画像データを利用することで、上述した空間コンパウンドスキャンにより得られる効果を得ることができる。しかしながら、空間コンパウンドスキャンの抱える問題も同時に抱えてしまう為、いずれにしろ空間コンパウンドスキャンの抱える上述した問題を解決しなくてはならない。
【００２０】
さらには、例えばパノラミックビューによる場合のように、異なる位置からスキャンして取得した画像データには、それぞれ取得した位置によって、当該画像データ中に鮮明な領域と不鮮明な領域とが存在する。例えば図１１に示すように、関心領域１３１の右側に不鮮明領域が存在する第１画像１３３と、関心領域１３１の左側に不鮮明領域が存在する第２画像１３５とについて、平均値演算や最大値演算を実行して合成処理を行うと、不鮮明な領域が拡大した画像データしか得ることができない。
【００２１】
従って、このような場合には、当該合成処理に係る画像データ中の不鮮明な領域は用いずに、鮮明な領域のみを用いて合成処理を行うことで、鮮明な領域のみを有する画像データ１３７を得ることが好ましい。しかしながら、現在、このような合成処理を自動的に行う技術は開示されていない。
【００２２】
そして、コンビネーションフォーカスによる場合は、画像データを分割して合成する為、つなぎ目（境界）が、当該合成処理後の画像データ上で見えてしまうという問題がある。つなぎ目（境界）について少しずつ重みを変えて加算するという方法もあるが、重みを変えて加算する方法よりも分解能向上に寄与するような合成処理方法がより好ましい。
【００２３】
なお、特許文献１及び特許文献２に開示されている技術は、当然ながら図８乃至図１１を参照して説明した上述の問題を解決するものではない。
【００２４】
本発明は、前記の事情に鑑みてなされたものであり、合成処理により作成した画像データにおいて、当該合成処理に係る画像データ同士の重複する領域と重複しない領域との境界が不連続に見えないこと、高輝度のアーティファクトが強調されてしまわないこと、当該合成処理に係る画像データに不鮮明な領域を有する画像データが含まれていた場合であってもこれに起因する不鮮明な領域が生じないこと、を実現した超音波画像診断装置及びプログラムを提供することを目的とする。
【課題を解決するための手段】
【００２５】
上記の目的を達成するために、本発明の第１の態様による超音波診断装置は、被検体における同一領域を少なくとも一部重複させるように超音波スキャンを複数回行い、受信エコーに基づいて複数フレームの画像を生成する超音波診断装置であって、前記複数フレームの画像のそれぞれに対して、所定の変換処理により多重解像度解析を行う解析手段と、前記解析手段による多重解像度解析により取得した各解像度の各係数毎に、前記複数フレームの画像間の対応する係数同士でそれぞれフィルタ演算処理を行うフィルタ演算手段と、前記フィルタ演算手段によるフィルタ演算処理結果に対して、前記所定の変換処理の逆変換処理を施して１フレームの画像を生成する画像生成手段と、を具備することを特徴とする。
【００２６】
上記の目的を達成するために、本発明の第２の態様によるプログラムは、被検体における同一領域を少なくとも一部重複させるように超音波スキャンを複数回行い、受信エコーに基づいて複数フレームの画像を生成する超音波診断装置としてコンピュータを機能させるプログラムであって、前記複数フレームの画像のそれぞれに対して、所定の変換処理により多重解像度解析を行う解析機能と、前記解析機能による多重解像度解析により取得した各解像度の各係数毎に、前記複数フレームの画像間の対応する係数同士でそれぞれフィルタ演算処理を行うフィルタ演算機能と、前記フィルタ演算機能によるフィルタ演算処理結果に対して、前記所定の変換処理の逆変換処理を施して１フレームの画像を生成する画像生成機能と、をコンピュータに実現させることを特徴とする。
【発明の効果】
【００２７】
本発明によれば、合成処理により作成した画像データにおいて、当該合成処理に係る画像データ同士の重複する領域と重複しない領域との境界が不連続に見えないこと、高輝度のアーティファクトが強調されてしまわないこと、当該合成処理に係る画像データに不鮮明な領域を有する画像データが含まれていた場合であってもこれに起因する不鮮明な領域が生じないこと、を実現した超音波画像診断装置及びプログラムを提供することができる。
【発明を実施するための最良の形態】
【００２８】
以下、図面を参照して本発明の一実施形態を説明する。
【００２９】
図１は、本一実施形態に係る超音波診断装置の構成例を示す図である。同図に示すように本一実施形態に係る超音波診断装置１は、超音波プローブ２と、送信回路４と、受信回路６と、Ｂモード処理手段８と、カラードプラ処理手段１０と、座標変換手段１２と、画像合成手段１４と、モニタ１６と、制御手段１８と、を具備する。
【００３０】
前記超音波プローブ２は、被検体との間で照射・反射される超音波信号の送受波を担うデバイス（探触子）であり、電気／機械可逆的変換素子としての圧電セラミック等の圧電素子で形成されている。この超音波プローブ２は、例えばアレイ状に配列される複数の圧電素子を先端部に装備したフェーズドアレイタイプのもので構成される。これにより、超音波プローブ２は、供給されるパルス駆動電圧を超音波パルス信号に変換して被検体のスキャン領域内の所望方向に送信し、且つ被検体から反射してきた超音波信号をこれに対応する電圧のエコー信号に変換する。
【００３１】
前記送信回路４は、前記制御手段１８による制御信号に基づいて、送信チャンネル毎に所定の送信遅延時間が付与されたタイミングで、前記超音波プローブ２の各圧電素子に駆動信号を送信し、これにより前記超音波プローブ２の各圧電素子から被検体内に向けて超音波信号を送波させる。
【００３２】
前記受信回路６は、超音波信号の送波に応じて被検体内の音響インピーダンスの不整合面で反射され、組織内の散乱体によって散乱された成分等を含む超音波エコー信号を、前記超音波プローブ２の各圧電素子を介してそれに対応する電圧量のエコー信号として受信し、そのエコー信号に受信遅延・加算処理を施し、その処理が施された受信信号を前記Ｂモード処理手段８及びカラードプラ処理手段１０にそれぞれ出力する。
【００３３】
前記Ｂモード処理手段８は、前記受信回路６からの出力信号に包絡線検波を行い、これで検波された信号を被検体内の組織の形態情報を担うものとして、前記座標変換手段１２に出力する。
【００３４】
前記カラードプラ処理手段１０は、前記受信回路６からの出力信号に基づいて速度情報を周波数解析し、その解析結果を被検体内の血流又は組織の移動速度情報を担う信号として、前記座標変換手段１２に出力する。
【００３５】
前記座標変換手段１２は、前記Ｂモード処理手段８及び前記カラードプラ処理手段１０からの出力信号を、表示座標へ座標変換して前記画像合成手段１４へ出力する。
【００３６】
前記画像合成手段１４は、前記座標変換手段１２からの信号を受信し、制御手段１８による制御で、Ｂモード像及びＣＦＭ（カラー・フロー・マッピング）像に関する各種画像処理及び合成処理（詳細は後述）を行い、更にこれらの画像に基づいて各種の定量解析や計測を行いその結果を示す情報を画像上に付加する等の画像処理を行い、その結果得られた画像信号をＴＶ用のスキャン信号に変換してモニタ１６に出力する。
【００３７】
前記制御手段１８は、図示しない内部バスに接続されたＣＰＵ（プロセッサ）、メモリ（ＲＡＭ／ＲＯＭ）、ハードディスク装置、リムーバブルメディア（ＣＤ−ＲＯＭ、フロッピー（登録商標）ディスク、メモリカード等）の駆動装置、及びその他の周辺装置を有するコンピュータとしての機能を有し、検査時に予めプログラムされた手順に従って当該超音波診断装置１全体の動作を制御する。
【００３８】
以上説明した構成により、前記モニタ１６上には、超音波画像（Ｂモード像及びＣＦＭ像に関する各種の画像や、その計測、解析結果に関する情報等を含む）が表示される。
【００３９】
なお、空間コンパウンドスキャンにおいてはフレーム毎に前記超音波プローブ２による斜めスキャンの方向（角度）を変化させるのであるが、フレーム毎にスキャンを変化させるのはＢモード像のみとし、ＣＦＭ像はフレーム毎に同じスキャンとする。つまり、ここでは空間コンパウンドスキャン及び合成処理を行うのはＢモードのみとする。以下、Ｂモード像の処理についてのみ詳細に説明する。
【００４０】
以下、本第１実施形態に係る超音波診断装置１の主な特徴の一つである前記画像合成手段１４による合成処理について詳細に説明する。
【００４１】
まず、前記画像合成手段１４による合成処理についての入出力関係は、例えば図２に示すような入出力関係となる。すなわち、Ｎフレーム（あるいはＮフレーム以上）の演算によって画像を合成し、入力と同じフレーム数の画像データを出力する。より具体的には、図２に示すように画像データＡ１と画像データＢ１とを合成して画像データＦ１を生成し、画像データＢ１と画像データＡ２とを合成して画像データＦ２を生成し、画像データＡ２と画像データＢ２とを合成して画像データＦ３を生成する。
【００４２】
なお、図２においては、説明の便宜上２方向のコンパウンド処理の場合を例に示しているが、本一実施形態に係る超音波診断装置によりＮ方向のコンパウンド処理を行うことが可能であるのは勿論である。
【００４３】
次に、前記画像合成手段１４による合成処理の流れについて、図３に示すフローチャートを参照して説明する。ここでは、（ｋ＋１）フレームの画像データを合成処理することを例に説明する。なお、それぞれのステップにおける処理の詳細は、図４を参照して後に詳述する。
【００４４】
まず、前記画像合成手段１４は、最新の画像データ（ｎフレーム目の画像）が入力される（ステップＳ１）と、Ｌレベルの離散的２次元ウェーブレット変換を行う（ステップＳ２）。
【００４５】
続いて、前記画像合成手段１４は、合成処理に係る画像同士の重複する領域に存在するエッジ成分を含むエッジ領域を除去する（ステップＳ３）。通常、画像の辺縁部においては、ウェーブレット変換を行うと高周波側にエッジ成分が発生する。そして、合成処理後の画像においては、当然ながら前記画像の辺縁部は、もはや辺縁部ではなくなる。従って、前記辺縁部におけるエッジ成分は不要になる。この為、ステップＳ３においてエッジ成分の除去を行う。
【００４６】
その後、前記画像合成手段１４は、後述する所定の演算処理（(k+1)枚のフレームの各ウェーブレット係数毎に演算を行う処理）を行い（ステップＳ４）、該演算処理結果に対して２次元離散逆ウェーブレット変換を行って（ステップＳ５）、合成処理後の画像データであるA(0)を得る。
【００４７】
以下、本一実施形態における画像データの合成処理の原理を、図４を参照して詳細に説明する。ここでは、画像ａ及び画像ｂの２つの画像を合成処理する場合を例に説明する。具体的には、前記画像合成手段１４は、次のようにして画像ａ及び画像ｂの多重解像度解析（ここではレベル２の２次元離散ウェーブレット変換）を行う（ステップＳ２）。
【００４８】
すなわち、前記画像合成手段１４は、まずレベル１の２次元離散ウェーブレット変換によって、A1(Approximation), H1(Horizontal detail), V1(Vertical detail), D1(Diagonal detail)のウェーブレット係数に分解する。続いて、前記A1についてのみ、更に２次元離散ウェーブレット変換を行う。これにより、図４に示すようにA2,H2,V2,D2を得る。これらの処理を、画像ａ及び画像ｂの双方についてそれぞれ行う。
【００４９】
続いて、図４において楕円形状５１で囲って示す領域であるエッジ領域を除去する（当該エッジ領域中の高周波側の係数を０にする）処理を行う（ステップＳ３）。
【００５０】
なお、このエッジ領域は、合成処理後の画像においては辺縁部とならないエッジ領域である。換言すれば、このエッジ領域は、合成処理前の或る１枚の画像にとっては当該画像の辺縁部であるが、合成処理後の画像においては辺縁部とはならない領域であり、この領域は当該合成処理前に既知の領域である。同様に、ウェーブレット変換を行った場合に、その画像の辺縁部が影響を与える範囲の領域も当該合成処理前に既知の領域である。
【００５１】
以上説明したステップＳ１乃至ステップＳ３の処理によれば、当然ながら同じ画像からは同じ結果を得ることができる。従って、現在の結果を順次保存しておくことで、毎回再度演算することなく図３において破線で囲った枠内に示す(n-k)フレームから(n-1)フレームまでの離散的２次元ウェーブレット変換の値を得ることができる。
【００５２】
そして、前記ステップＳ４における演算処理及び前記ステップＳ５における２次元離散逆ウェーブレット変換は、例えば次のような演算を行う処理である。
【００５３】
ここで、演算関数をfunc(a1,a2,・・・,aN, M, Level) で示し、離散的２次元逆ウェーブレット変換をidwt2で示す。但し、a1,a2,・・・,aNは入力データを示し、M=0はApproximationを示し、M=1はHorizontal detailを示し、M=2はVertical detailを示し、M=3はDiagonal detailを示し、Levelはウェーブレット展開数を示す。
【００５４】
A(L)=func(A(L,n),A(L,n-1),・・・,A(L,n-k),0,L)
H(L)=func(H(L,n),H(L,n-1),・・・,H(L,n-k),1,L)
V(L)=func(V(L,n),V(L,n-1),・・・,V(L,n-k),2,L)
D(L)=func(D(L,n),D(L,n-1),・・・,D(L,n-k),3,L)
H(L-1)=func(H(L-1,n),H(L-1,n-1),・・・,H(L-1,n-k),1,L-1)
V(L-1)=func(V(L-1,n),V(L-1,n-1),・・・,V(L-1,n-k),2,L-1)
D(L-1)=func(D(L-1,n),D(L-1,n-1),・・・,D(L-1,n-k),3,L-1)
・・・・・・・・
H(1)=func(H(L,n),H(L,n-1),・・・,H(L,n-k),1,1)
V(1)=func(V(L,n),V(L,n-1),・・・,V(L,n-k),2,1)
D(1)=func(D(L,n),D(L,n-1),・・・,D(L,n-k),3,1)
A(L-1)=idwt2(A(L),H(L),V(L),D(L))
A(L-2)=idwt2(A(L-1),H(L-1),V(L-1),D(L-1))
・・・・・・・
A(0)=idwt2(A(1),H(1),V(1),D(1))
ここで、前記A(0)が合成処理により最終的に作成された画像データである。
【００５５】
なお、上述した演算関数の一例としては、すべてのLevelにおいて、
func(a1,a2,・・・,aN, M, Level) = Mean(a1,a2,・・・,aN)    M=0
func(a1,a2,・・・,aN, M, Level) = AbsMax(a1,a2,・・・,aN)  M=1,2,3
を挙げることができる。
【００５６】
但し、Mean(a,b) は、aとbとの平均値を計算する演算であるが、a=0の場合はbの値を出力し、b=0の場合はaの値を出力する演算である。また、AbsMax(a,b) はa及びbのうち何れか絶対値の大きい方を出力する演算である。
【００５７】
なお、より具体的には、前記ステップＳ４においては、例えば次のような演算を行う。
【００５８】
A2 = Mean(A2a,A2b)
H2 = AbsMax(H2a,H2b)
V2 = AbsMax(V2a,V2b)
D2 = AbsMax(D2a,D2b)
H1 = AbsMax(H1a,H1b)
V1 = AbsMax(V1a,V1b)
D1 = AbsMax(D1a,D1b)
そして、上述した演算を行った後、前記画像合成手段１４は、前記の係数A2,H2,V2,D2から２次元離散逆ウェーブレット変換を行ってA1を算出し、A1,H1,V1,D1から２次元離散逆ウェーブレット変換を行って（ステップＳ５）、合成処理後の画像データであるA(0)を得る。
【００５９】
以下、前記ステップＳ４における演算処理を行うことで得られる効果について説明する。
【００６０】
まず、前記Approximationの係数については平均値演算を行うことで、平均輝度が変化しないという効果、及び画像の大まかな構造は合成処理に係る各々の画像間の平均値で示されることになるという効果を得ることができる。
【００６１】
また、前記H,V,Dのdetail係数については最大値演算を行うことで、それぞれの画像のうち最大の分解能を持つ画像が選択され、各画素毎に最大の分解能を持つ画像が自動的に選択される。従って、全体としては重ね合わせの境界も見えずに自然な重ね合わせとして見え、各画素においては重ね合わせた画像の中の最大の分解能を持つ画像が表示されるので、分解能が向上して見えるという効果を得ることができる。
【００６２】
つまり、画像データ同士の重複する領域と重複しない領域との境界が不連続に見えないこと、高輝度のアーティファクトが強調されてしまわないこと、画像データに不鮮明な領域を有する画像データが含まれていた場合であってもこれに起因する不鮮明な領域が生じないこと、という目的を達成することができる。
【００６３】
このような演算により、コンパウンドスキャンで期待される鏡面反射による画像の鮮明化が達成される。ここで、更にスペックル除去という効果を得る為には、ウェーブレット変換を行うことを利用して、ウェーブレット縮退(Ｗａｖｅｌｅｔ  Ｓｈｒｉｎｋａｇｅ)という画像処理技術を併用すればよい。レベルの低いdetail係数はノイズである場合が多いので、ウェーブレット縮退を行うことで、スペックルを除去又は低減することができる。
【００６４】
なお、この場合の前記func関数としては、例えば次の関数を挙げることができる。
【００６５】
func(a1,a2,・・・,aN, M, Level) = AbsMax(a1,a2,・・・,aN) * c(M, Level)    M=1,2,3
ここで、c(M, Level) はウェーブレット係数の重みを示し、例えばLevelが１のような小さい場合に０に近い値を取り、Levelが大きい場合は１に近い値を取る。また、c(M, Level) の値は、Mによって指定されるHorizontal, Vertical, Diagonalの場合でそれぞれ値を変えても良い。さらには、c(M, Level) の値は、組織の辺縁等のエッジを強調する為に、１より大きい値を取っても良い。
【００６６】
なお、以上説明した前記画像合成手段１４による処理は、例えば図５に示す互いに重複するような複数回の３次元スキャンにより取得した画像データに対しても適用できる。
【００６７】
図５は、３次元スキャンをしながら前記超音波プローブ２を動かしたり、別のウィンドウから同一の組織を観察する様子、つまり一旦前記超音波プローブ２を被検体から離して当該被検体の別の位置へ置きかえて再び画像を取得するというプロセスを繰り返す様子を示す図である。
【００６８】
この方法によれば、被検体における同一の領域に対して異なる位置から複数回３次元スキャンを実施して、取得した画像の位置合わせを３次元的に行うことにより、完全に同一断面又は同一ボクセルのデータを得ることができる。この場合、同一の領域について複数のボクセルデータが存在し、複数のボクセルデータ又は切り出された複数の２次元画像データを合成処理することになる。そして、このような方法を採る場合にも、上述した前記画像合成手段１４による合成処理を適用することができることは勿論である。
【００６９】
ところで、前記画像合成手段１４による上述した一連の合成処理は、プログラム化することで、或いはプログラム化した後当該プログラムを記憶媒体に読み込むことによって、超音波診断装置１とは独立したソフトウェア製品単体としての販売、配布も容易になり、また本一実施形態に係る技術を他のハードウェア上で利用することも可能となる。
【００７０】
以上説明したように、本一実施形態によれば、合成処理により作成した画像データにおいて、当該合成処理に係る画像データ同士の重複する領域と重複しない領域との境界が不連続に見えないこと、高輝度のアーティファクトが強調されてしまわないこと、当該合成処理に係る画像データに不鮮明な領域を有する画像データが含まれていた場合であってもこれに起因する不鮮明な領域が生じないこと、を実現した超音波画像診断装置及びプログラムを提供することができる。
【００７１】
ところで、上述した前記画像合成手段１４による合成処理によって得られた画像データにおいては、当該合成処理前の画像データ同士の重複する領域と重複しない領域との境界が自然であり且つ分解能が向上する。分解能が向上する理由は、当該合成処理後の画像データは、合成処理前の画像データにおける高周波成分のうち最も大きい成分を持つ為、最も解像度が高い状態で合成処理されるからである。
【００７２】
以上、一実施形態に基づいて本発明を説明したが、本発明は上述した実施形態に限定されるものではなく、本発明の要旨の範囲内で種々の変形及び応用が可能なことは勿論である。
【００７３】
（第１変形例）
合成処理における演算方法としては上述した演算方法以外にも、例えば全ての場合においてMean演算を用いる方法、又は全ての場合においてAbsMax演算を用いる方法を採ってもよい。
【００７４】
全ての場合においてMean演算を用いる方法は、ウェーブレット変換を行わずに単にMeanを取る処理と基本的には同等の処理であるが、本一実施形態においては上述したようにウェーブレット係数上で重複する側のエッジ領域除去を行っている為、重複する領域における境界が不連続に見えることがない。
【００７５】
全ての場合においてAbsMax演算を用いる方法は、ウェーブレット変換を行わずに単にMaxを取る処理とは異なり、各画素毎に最大の分解能を持つ画像が出力されることになる。但し、平均輝度が上昇したり、高輝度のアーティファクトが残る場合がある。
【００７６】
（第２変形例）
合成処理における上述した演算方法では、各画素毎の係数のみを用いて合成処理を行ったが、当該画素の近傍の画素の係数を用いて合成処理を行っても勿論良い。
【００７７】
例えば、２枚の画像の位置(x,y)におけるAbsMaxによる計算方法として、１枚目の画像における位置(x,y)を含む周囲５×５点の平均値a1を計算し、且つ２枚目の画像における位置(x,y)を含む周囲５×５点の平均値a2を計算し、前記a1及び前記a2のうち絶対値の大きい方の値を位置(x,y)における出力値とする。この処理により、ノイズが強調されて見えることを低減することができる。
【００７８】
（第３変形例）
合成処理における上述した演算方法では、多重解像度解析に離散的ウェーブレット変換を用いた。一般的な離散的ウェーブレット変換では、当該画像データに対してＬＰＦ(Ｌｏｗ  Ｐａｓｓ  Ｆｉｌｔｅｒ)とＨＰＦ(Ｈｉｇｈ  Ｐａｓｓ  Ｆｉｌｔｅｒ)を掛け、それぞれの出力を半分にダウンサンプリングする。この為、情報量は当該変換前と変換後とでは変わらない。
【００７９】
しかしながら、ウェーブレット縮退で係数を処理した後に、逆ウェーブレット変換を行うと、当該画像上にブロック状のアーティファクトが見える場合がある。これを防ぐ為に、離散的ウェーブレット変換の代わりに、例えばステーショナリ(ｓｔａｔｉｏｎａｒｙ)ウェーブレット変換を行っても良い。
【００８０】
このステーショナリウェーブレット変換では、ダウンサンプリングを行わない。この為、２次元画像を処理する場合には、１回の変換によって情報量が４倍になってしまう。しかしながら、ウェーブレット縮退を行った後に逆変換を行った場合であっても、当該画像上に上述したブロック状のアーティファクトは見られないという効果を得ることができる。
【００８１】
（第４変形例）
合成処理における上述した演算方法では、多重解像度解析に離散的ウェーブレット変換を用いた。しかしながら、多重解像度変換にウェーブレット変換の代わりに、例えばラプラシアン(Ｌａｐｌａｃｉａｎ) ピラミッド、ＲｏＬ(Ｒａｔｉｏ  ｏｆ  Ｌｏｗ  ｐａｓｓ)ピラミッド、又はグラディエント(Ｇｒａｄｉｅｎｔ)ピラミッド、といった各種のピラミッド変換(Ｐｙｒａｍｉｄ  ｔｒａｎｓｆｏｒｍ)を用いても良い。
【００８２】
ウェーブレット変換とピラミッド変換との相違点は、ウェーブレット変換が直交変換であるのに対して、ピラミッド変換は必ずしも直交変換ではない点である。しかしながら、ウェーブレット変換及びピラミッド変換の双方共、それらの変換によって多重解像度解析が可能である点では共通する。
【００８３】
（第５変形例）
図６に示すようにスペックル低減やエッジを強調する為の非線形拡散フィルタ処理（ステップＳ１１）を、前記ステップＳ４における演算処理の前或いは前記ステップＳ４における演算処理の後の何れか一方又は両方に挿入しても良い。なお、図６においては、説明の便宜上、前記ステップＳ４の前後双方に非線形拡散フィルタ処理（ステップＳ１１）を挿入して示している。
【００８４】
すなわち、例えば特開２００６−１１６３０７に開示されているような、多重解像度解析と非線形拡散フィルタとを組み合わせた処理を、合成処理前に挿入することで、スペックルが少なく且つエッジが強調された見やすい合成画像を得ることができる。このような処理を行う場合には、多重解像度解析は既に行われているので、この処理に掛かる時間を節約することができる。
【００８５】
（第６変形例）
前記一実施形態においてはコンパウンドスキャンを用いる場合を例に説明したが、コンパウンドスキャンに限らず、例えばパノラミックビュー又は３次元スキャンを行いながら超音波プローブ２を動かしたり、別のウィンドウから同一の組織を観察する場合であっても同様に適用できる。
【００８６】
なお、パノラミックビューの場合には、画像から動き量を推定して、当該画像を並行移動及び／又は回転移動させる必要がある。詳細は、前記特許文献１に開示されている。別の位置から３次元スキャンを行う場合も、３次元的に並行移動及び／又は回転移動させることで、断面を所望の断面に合わせることができる。このようにして断面を一致させた後は、前記一実施形態において説明した画像の合成処理を行うことができる。
【００８７】
（第７変形例）
コンビネーションフォーカス法による合成処理に、前記一実施形態を適用してもよい。この場合、合成処理に係る画像の境界をオーバーラップさせて、前記一実施形態で説明した画像合成手段１４による合成処理を行う。
【００８８】
図７は、２段のコンビネーションフォーカス法の概念を示す図である。近距離フォーカスの超音波ビーム１１１により取得した画像５１Ａと、遠距離フォーカスの超音波ビーム５３による画像５３Ａとの境界をオーバーラップさせてオーバーラップ領域５５Ａを作り、このオーバーラップ領域５５Ａについて、前記一実施形態において説明した合成処理方法を用いる。
【００８９】
これにより、合成処理後の画像における分解能の向上且つスペックルの低減が期待できる。この場合、前記画像合成手段１４による処理と、前記座標変換手段１２による処理との順序を逆にしてもよい。
【００９０】
（第８変形例）
周波数コンパウンドスキャンによる合成処理に、前記一実施形態を適用してもよい。なお、基本波画像と２次高調波画像との合成も周波数コンパウンドの一種である。周波数コンパウンドスキャンによれば、通常は、各々のスキャン範囲は一致するので、コンビネーションフォーカス法と同様にして前記一実施形態を適用することができる。
【００９１】
さらに、上述した実施形態には種々の段階の発明が含まれており、開示される複数の構成要件の適当な組み合わせにより種々の発明が抽出され得る。例えば、実施形態に示される全構成要件からいくつかの構成要件が削除されても、発明が解決しようとする課題の欄で述べた課題が解決でき、発明の効果の欄で述べられている効果が得られる場合には、この構成要件が削除された構成も発明として抽出され得る。
【図面の簡単な説明】
【００９２】
【図１】本発明の一実施形態に係る超音波診断装置の構成例を示す図。
【図２】画像合成手段の入出力を示す図。
【図３】画像合成手段による画像データの合成処理のフローチャートを示す図。
【図４】本発明の一実施形態における画像データの合成処理の原理を示す図。
【図５】３次元スキャンの概念を示す図。
【図６】スペックル低減やエッジを強調するフィルタ処理を挿入する場合の、画像合成手段による画像データの合成処理のフローチャートを示す図。
【図７】２段のコンビネーションフォーカス法の概念を示す図。
【図８】空間コンパウンドスキャンの概念を示す図。
【図９】パノラミックビューの概念を示す図。
【図１０】コンビネーションフォーカスの概念を示す図。
【図１１】従来技術が抱える問題点を示す図。
【符号の説明】
【００９３】
１…超音波診断装置、  ２…超音波プローブ、  ４…送信回路、  ６…受信回路、  ８…Ｂモード処理手段、  １０…カラードプラ処理手段、  １２…座標変換手段、  １４…画像合成手段、  １６…モニタ、  １８…制御手段。
【図１】
0005002397.tif 000002
【図２】
0005002397.tif 000003
【図３】
0005002397.tif 000004
【図４】
0005002397.tif 000005
【図５】
0005002397.tif 000006
【図６】
0005002397.tif 000007
【図７】
0005002397.tif 000008
【図８】
0005002397.tif 000009
【図９】
0005002397.tif 000010
【図１０】
0005002397.tif 000011
【図１１】
0005002397.tif 000012