(57)【特許請求の範囲】
【請求項１】
検査対象に対して超音波を送受信するための超音波探触子を含むユニットと、（n-１）番目のフレームの二次元超音波画像であるｆn-1を保持するための画像保持部と、前記検査対象の動きベクトルを計測するための少なくとも１つの領域を設定するための計測領域設定部と、前記ｆn-1とｎ番目のフレームの二次元超音波画像であるｆnとから前記検査対象の動きベクトルを計測するための動きベクトル計測部と、１番目のフレームの二次元超音波画像であるf1から前記ｆn-1までの複数の二次元超音波画像を加算して累積加算画像を作成する加算画像作成部と、前記加算画像作成部で作成された前記累積加算画像を保持するための加算画像保持部と、前記動きベクトルに基づいて前記累積加算画像を補正するための画像補正部と、前記画像補正部で補正された前記累積加算画像と前記ｆnとの少なくとも一方に重み係数を乗算して処理する重み付け部と、前記加算画像保持部に保持された前記累積加算画像を表示するための表示部とを有し、前記加算画像作成部は、前記重み付け部で処理された前記ｆnおよび前記累積加算画像を加算して新たな累積加算画像を作成して前記加算画像保持部に保持させ、
前記動きベクトル計測部は、最新に取得したフレームの二次元超音波画像ｆnとその１フレーム前の二次元超音波画像ｆn-1との間で動きベクトルを計測し、
前記画像補正部は、計測した動きベクトルに基づいて１フレーム前までの二次元超音波画像である前記ｆ1から前記ｆn-1までを用いて構成した前記累積加算画像を補正し、前記加算画像作成部は補正された前記累積加算画像と前記最新に取得したフレームの二次元超音波画像を加算処理して新たな前記累積加算画像を構成することを特徴とする超音波診断装置。
【請求項２】
前記計測領域設定部は前記ｆn-1上に前記領域を設定することを特徴とする請求項１に記載の超音波診断装置。
【請求項３】
前記加算画像作成部は、前記最新に取得したフレームの二次元超音波画像と補正された前記累積加算画像との双方に一定の重み付けを乗算して加算処理を行なうことを特徴とする請求項１に記載の超音波診断装置。
【請求項４】
前記重み付け部は最新に取得したフレームの二次元超音波画像と前記累積加算画像との双方の画像に一定の重み付け処理をし、前記加算画像作成部は重み付けされた前記双方の画像を加算処理して新たな前記累積加算画像を構成することを特徴とする請求項１に記載の超音波診断装置。
【請求項５】
前記動きベクトルの計測精度を示す相関値に応じて、前記重み係数が自動的に調節されることを特徴とする請求項１に記載の超音波診断装置。
【請求項６】
前記重み係数を変更するためのダイヤルをさらに有することを特徴とする請求項１に記載の超音波診断装置。
【請求項７】
前記表示部は前記ダイアルによって変更された前記重み係数の値を表示することを特徴とする請求項６に記載の超音波診断装置。
【請求項８】
前記加算画像への重み係数を０にするための消去手段を有することを特徴とする請求項１に記載の超音波診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、スペックルノイズを低減し、高コントラスト超音波画像を表示する装置に関する。
【背景技術】
【０００２】
超音波診断装置は、患部を非侵襲に実時間で撮像できるため、外来における初期診断や胎児診断、治療のモニタリングツールとして広く利用されている。近年ではポータブルタイプの超音波診断装置の開発も進み、病院内外問わず、患者の元に装置を持っていく救急処置も可能になっている。また、装置の小型に加えて、超音波診断装置は操作が容易であるため、患者宅と病院とのネットワーク整備を進めることにより、患者本人が装置を操作しながら画像を病院に転送し、遠隔的に診断するという在宅医療の形態も将来的に期待される。【０００３】
超音波画像は撮像時間が短く実時間画像表示が可能な反面、ノイズレベルに対する組織構造を示す信号強度（S/N比）が低く、画像の判別が困難な場合がある。したがって診断の精度向上や医者と患者の情報共有を進めるためには、より客観的な情報を示す判り易い画像表示技術が期待されている。【０００４】
超音波画像には、スペックルと呼ばれる斑点の模様があり、画像の判読を難しくする一つの要因となっている。このスペックルを軽減し、高画質の超音波画像を表示する技術として、空間コンパウンド法がある。超音波二次元画像（B-mode画像）上のスペックルは、検査対象の内部に散在する微小散乱体からの反射信号が互いに干渉して生じる、空間的に定在した信号であり、その空間周波数は検査対象の構造を示す信号に比べて低い。そのため、撮像面をスペックルの大きさと同程度にスライス方向に移動させると、画像上においてスペックルがランダムに変化する一方で、検査対象の構造には大きな変化が生じない。このようにスライス方向に僅かに異なる位置で取得した画像同士を加算処理すると、検査対象の構造からの信号は加算毎に強調されるのに対し、ランダムに変化するスペックルの信号は平滑化されていく。その結果、スペックルと構造とのS/N比が向上し、高いコントラストを呈した画像が構成される。【０００５】
複数の画像を加算処理し、特定の情報を抽出する技術は複数報告されており、その中のいくつかを以下に説明する。
特許文献１の技術では、血流画像（CFM:Color Flow Mapping）を複数枚メモリに保持し、表示される全画像から加算処理に用いるフレームを術者が選択して除き、残った画像を加算処理して画像を構成している。【０００６】
特許文献２に記載の技術は、少量の造影剤投与で安定に効率よく造影剤からの信号を受信し、血管の分布を明瞭に画像化する技術に関するものである。血管情報を画像化するには血流内を流れる造影剤を超音波により圧壊させ、その際に生じる強い非線形信号を受信して画像を構成する。造影剤を圧壊させるには強い超音波強度が必要なため、送波信号の集束領域とそれ以外の領域とでS/Nが大きく異なり、造影剤画像の深さ方向にむらができる課題があった。この課題に対して、特許文献２に記載の技術は、送波信号の集束領域を深さ方向に多段階に変えてマルチショット画像群と呼ばれる複数の画像を構成し、それら画像を加算処理する。また、マルチショット画像群の各フレームの特定領域に重み付けをしてから加算処理する。【０００７】
【特許文献１】特許第3447150号公報
【０００８】
【特許文献２】特開第2002-209898号公報
【０００９】
複数枚の画像を加算もしくは減算処理して特定の情報を強調した画像を構成し、リアルタイムで表示するには、被検体の位置ずれ、処理時間の増大、メモリの増加といった問題が生じる。【００１０】
特許文献１に記載の技術は、検査対象の位置ずれを考慮していないため動きを伴う検査対象への適用が困難であり、また加算処理に用いる画像を術者が選定するため、リアルタイム表示も困難である。特許文献２に記載の技術は、加算するフレーム全てを一旦メモリに保持しなければならず、表示する一枚の加算処理画像を構成するためには、加算処理に用いるフレーム数以上のメモリが必要である問題がある。
【発明の開示】
【００１１】
本発明の目的は、スペックルを低減した高コントラスト超音波画像を表示する超音波診断装置を提供することにある。
【００１２】
本発明は時間変化する情報を抽出することを目的とする。本発明は、重み付けを累積加算画像に対して行い、それにより時間変化する情報の抽出量を術者が制御などして利用する。
【００１３】
上記目的を達成するために、本発明の超音波診断装置では、取り込んだB-mode画像とその直前に取得したB-mode画像との間に生じる検査対象の動きベクトル（体動の情報）を計測し、計測した動きベクトルに基づく画像変形処理を行ないながら、重み係数を乗した画像の加算処理を行ない、検査対象の輪郭構造を強調した画像を表示する。
【００１４】
以下、本発明の超音波診断装置の代表的な構成例について列挙する。
（１）検査対象に対して照射する超音波を送受信するための超音波探触子と、受波信号から検査対象の断層像を示すB-mode画像を構成するための超音波画像構成ユニットと、１フレーム前のB-mode画像fn-1を保持するための画像保持部と、fn-1画像上に検査対象の動きベクトルを計測するための少なくとも１つの領域を設定するための計測領域設定部と、fn-1とfnの間に生じる検査対象の動きベクトルを計測するための動きベクトル計測部と、f1からfn-1を用いて構成される累積加算画像を作成する加算画像作成部と、前記動きベクトルの計測結果に基づいて、f1からfn-1を用いて構成される累積加算画像を変形(補正)させる画像変形(補正)部と、取得画像および累積加算画像との少なくともいずれかに重み係数を乗するための重み付け部と、累積画像と取得画像とを加算処理をするための画像加算部と、次に取り込まれるfn+1への加算処理に用いるために前記画像加算部で構成した累積加算画像を保持するための加算画像保持部と、前記累積加算画像を表示するための表示部とを有する。以下、加算画像を変形させるとは、加算画像を動きベクトルに基づいて補正することを含むものとする。
（２）前記（１）の超音波診断装置において、動きベクトルの計測は取得したB-mode画像fnとその１フレーム前のB-mode画像fn-1との間で行ない、加算処理は取得したB-mode画像fnを基準画像とし、
0005171610.tif 000002
（３）前記（１）の超音波診断装置において、取得画像及び累積加算画像に重み係数を乗した加算処理を行ない、新たな加算画像を構成することを特徴とする。
（４）前記（１）の超音波診断装置において、取得画像及び累積加算画像に乗される重み係数は術者が自由に変更できることを特徴とする。
（５）前記（１）の超音波診断装置において、取得画像及び累積加算画像に乗される重み係数は、動きベクトルの計測の際に算出される相関値により自動的に調整され、例えば検査対象の動きが大きい場合には、取得画像に対する重み係数が０となることを特徴とする。
【００１５】
本発明によれば、スペックルノイズを低減し、組織の輪郭構造を強調した高コントラスト画像の構成が可能であり、視認性の高い超音波画像を表示できる。
【図面の簡単な説明】
【００１６】
【図１】本発明の実施例１の超音波診断装置の構成を示すブロック図。
【図２】実施例１の超音波診断装置において、B-mode画像の取り込みから累積加算画像の表示に至る処理を説明するフローチャート。
【図３】実施例１の超音波診断装置において、動きベクトルの計測手法を説明する図。
【図４】実施例１の超音波診断装置において、取得画像と前フレームまでの累積加算画像とを用いて、新たな累積加算画像を構成するための画像処理を説明する図。
【図５】実施例１の超音波診断装置において、動きベクトルの計測を高速処理するための処理を説明するフローチャート。
【図６】実施例１の超音波診断装置において、動きベクトルの計測手法を説明する図。
【図７】実施例１の超音波診断装置において、累積加算画像を構成する各フレームの重み係数を示す図。
【図８】実施例１の超音波診断装置において、重み係数（α，β）を調節するためのダイヤルを備えた探触子または超音波診断機の例を示す図。
【図９】実施例１の超音波診断装置において、重み係数（α，β）の調節機能とリフレッシュ機能とを備えたダイヤル兼ボタンを備えた探触子または超音波診断機の例を示す図。
【図１０】実施例１の超音波診断装置において、B-mode画像と加算画像とを並べて表示するバイプレーン表示の形態の一例を示す図。
【図１１】実施例２の超音波診断装置の構成を示すブロック図。
【図１２】実施例３の超音波診断装置の構成を示すブロック図。
【図１３】実施例３の超音波診断装置において、５枚の画像で構成されるユニット画像を説明する図。
【図１４】実施例３の超音波診断装置において、５枚のユニット画像および１枚の画像に対する重み係数を示す図。
【発明を実施するための最良の形態】
【００１７】
以下、本発明の実施例について、図面を参照して詳述する。
【実施例１】
【００１８】
図１は、本発明の一実施例になる超音波診断装置の構成を示すブロック図である。【００１９】
本実施例の超音波診断装置では、検査対象への超音波送受信により超音波二次元画像（以下、B-mode画像）を構成し、前記B-mode画像上に検査対象の動きベクトルを計測するための領域を設定し、取り込んだB-mode画像とその１フレーム前の画像間に生じる動きベクトルを前記計測領域毎に計測し、前記動きベクトルに基づいて累積加算画像を変形(補正)し、前記累積加算画像に重み係数を乗して取り込んだB-mode画像に加算し、新たな累積加算画像を構成して実時間で表示部に表示する。【００２０】
まず、超音波画像構成ユニットにおけるB-mode画像の構成から、位置補正を加えた画像加算処理を行なった加算画像の表示に至るまでの装置構成について、図１のブロック図を用いて説明する。【００２１】
超音波探触子（以下、探触子）２は複数の圧電素子が平行に並べられた構造をもつ。送信ビームフォーマ３からＤ/Ａ変換器４を経由して、各圧電素子にアナログ送波信号が送られ、検査対象１に向けて超音波を照射する。各圧電素子から送信される超音波は、送信ビームフォーマ３によって電子的に遅延が掛けられており、所定の深度で集束する。送波信号は、検査対象１内で反射され、再び探触子の各圧電素子で受信される。各圧電素子で受信した反射エコーは、ＴＧＣ（Time Gain Control）部５で、送波の到達深度によって変わる減衰分を補正された後、Ａ/Ｄ変換器６でデジタル信号に変換され、受信ビームフォーマ７に送られる。【００２２】
受信ビームフォーマ７では、焦点位置から各圧電素子までの距離に応じた遅延時間を掛けて加算結果が出力される。この集束超音波を二次元走査させることで、検査対象１の二次元的な反射エコー分布が得られる。受信ビームフォーマ７からは実部と虚部に分けられたＲＦ信号が出力され、包絡線検波部８と計測領域設定部１１に送られる。包絡線検波部８に送られた信号は、ビデオ信号に変換された後、スキャンコンバータ９で走査線間の補間が加えられ、二次元画像データであるB-mode画像が構成される。【００２３】
構成されたB-mode画像は動きベクトル計測部１０に送られる。ここで、画像保持部１１に動きベクトル計測部１０に取り込まれたB-mode画像よりも１フレーム前のB-mode画像が保持されている。構成されたB-mode画像が１フレーム目である場合は、動きベクトル計測部１０では処理されずに通過して画像保持部１１に入力される。画像保持部１１に保持されているB-mode画像は、計測領域設定部１２において動きベクトルを計測するための計測領域が、構造に応じた最適な大きさに設定される。計測領域が設定されたB-mode画像は、動きベクトル計測部１０に送られる。動きベクトル計測部１０では、計測領域設定部１２からのB-mode画像とスキャンコンバータ９から送られるB-mode画像とを用いて、前記計測領域内の動きベクトルが計測される。動きベクトルの計測手法は、相互相関演算または最小二乗法である。画像変形部１３では、前記動きベクトル計測部１０で計測された動きベクトルに基づいて、加算画像保持部１４から読み出された累積加算画像が変形される。重み付け部１６では取得画像および累積加算画像に重み係数が乗され、画像加算部１３にて取得画像と累積加算画像とが加算される。画像加算部１３で構成された累積加算画像は加算画像保持部１４に一旦格納された後、表示部１５にて表示される。加算画像保持部１４に格納された累積加算画像は、次の加算処理の際に重み付け部１６を経て画像加算部１３に送られる。【００２４】
次に、図２のフローチャートに従って、B-mode画像の取得から累積加算画像の表示に至るまでの工程の詳細を説明する。また、図４は画像処理の流れを示すブロック図である。【００２５】
まず、工程１でB-mode画像fnを構成する。fnが最初の画像（n=1）の時は、fnは画像保持部１１に格納され（工程２）、次のB-mode画像が取り込まれる（図４-４１）。fn(n>1)が動きベクトル計測部に送られるのと同時に、画像保存部１１に格納されている１フレーム前の画像fn-1が読み込まれ（工程３）、動きベクトルを計測するための領域が設定される（工程４）。続いて、fn 及びfn-1を用いて、両画像間に生じた動きベクトルが、設定した計測領域毎に計測される（工程５、図４-４２）。
0005171610.tif 000003
工程５で計測された動きベクトルに基づいて変形処理が為される（工程９、図４-４３）。変形処理が為された累積加算画像は一定の重み係数を乗する重み付け処理をする（工程１０）。
0005171610.tif 000004
画像加算部１３にてfnを基準にした加算処理が為される（工程６、図４-４４）。
0005171610.tif 000005
加算画像保持部１４に格納された後（工程７）、表示部１５に表示される（工程８）。【００２６】
工程４の計測領域の設定手法について図３を用いて説明する。
本発明では、B-mode画像（fn-1）２５に計測領域２４を複数個設定し、各計測領域において、最も整合が取れる領域を、取得したB-mode画像（fn）２１から相互相関演算または最小二乗法により探索する。一つ一つの計測領域内での動きは変形を伴わない剛体運動と見なし、各計測領域で求めた個々の動きベクトルを組み合わせることで、変形を含む検査対象全体の動きを計測する。計測領域をB-mode画像（fn-1）２５に設定するのは、後に説明する重み係数を用いた加算処理に用いる領域を、一定にするためである。重み係数は加算画像に対して大きく設定される。一方、計測領域を取得画像（fn）２１に設定すると、加算画像から抽出される領域が、加算領域とスペックルを含む非加算領域を含む。そのため、この領域を加算処理すると、重み係数によってスペックルが強調され、アーチファクトとなってしまう。【００２７】
計測に用いる信号成分は、大別して検査対象の輪郭構造や組織間の境界といった低周波成分と、検査対象に散在する微小散乱体により散乱した超音波が互いに干渉してなる高周波成分のスペックルがある。本発明では、両者を区別せず、画像全体に計測領域を設定して移動量を計算する。B-mode画像を用いて画像間に生じる動きを計測する上で、高周波成分であるスペックルは、動きベクトルの計測精度を向上させるだけでなく、特徴的な構造を持たない組織の実質部位における動きの計測を可能にさせる。計測領域の大きさは、B-mode画像の最小構成要素であるスペックルよりも大きい必要があり、スペックルサイズの約２倍から３倍の大きさに設定する。肝臓や腎臓の腹部組織では、約５×５mm2の計測領域を設定する。【００２８】
計測領域の設定に関しては、組織の部分毎に異なる大きさの計測領域を設定する手法も考えられる。動きベクトルの計測分解能を高くするには、計測領域を小さくする必要があるが、その場合は必然的に領域数が多くなり、計測に要する処理時間が長くなる。したがって、組織の輪郭といった特徴的な構造を含む計測領域は、他の領域と区別して当該領域内の空間周波数に応じて大きく設定する。例えば、図３の中で、肝臓組織２２、肝臓組織内の血管構造２３を含む画像２５の中に計測領域を設定する場合、血管構造を含む計測領域２６の大きさは、特徴的な構造物を含まない周囲の計測領域２４より大きく設定する。また、組織の境界を含む計測領域や変形が大きい計測領域では、周辺の領域よりも相関値が低く、組織の動きベクトルの計測が困難になる。したがって、相関値が低い計測領域では、計測領域をスペックルサイズの２倍を越えない程度に小さく再設定する手法も考えられる。【００２９】
次に、工程５及び工程６の動きベクトルの計測から画像加算処理に関する具体的な手法について、図５及び図６を用いて説明する。動きベクトルの計測は、取り込んだ画像fnと１フレーム前の画像fn-1とを用いて行なう。動きベクトルの計測に関する処理の流れを、図５に示すフローチャートを用いて説明する。まず、画像fn-1に動きベクトルを計測するための計測領域を、前述した手法を用いて設定する（工程１）。実際には画面全体を覆う形で複数の計測領域が設定されるが、ここではその中の一つの計測領域５１に着目し、その大きさを(i,j)とする。次に、計測領域fn-1'５１と最も整合する領域を探索するための探索領域fn'５２をfn上に設定する（工程２）。探索領域fn'５２の大きさは検査対象の動く速さとフレームレートを考慮して可能な限り小さく設定する。フレームレート２０-３０(frames/sec)で、呼吸の影響を受ける肝臓の場合、計測領域fn-1'５１を中心に周囲10ピクセル広い領域(i+20,j+20)の設定で動きベクトルの計測が可能である。次に、計測領域fn-1' ５１及び探索領域fn'５２に対してLow-Pass Filterを適用し（工程３）、１ピクセル間隔で画素を間引くデシメーション処理を行なう（工程４）。Low-Pass Filter処理は、その後に続くデシメーション処理による画素間のエイリアシングを抑制するための処理であり、デシメーション処理により領域探索に要する処理時間を約４倍の短縮が可能である。デシメーション処理の具体的な方法について図６を用いて説明する。デシメーション処理後の探索領域fn'５５内で、同じくデシメーション処理後の計測領域fn-1'５４を１ピクセルずつ走査させ、以下の（式１）または（式２）で定義される相関値ｃが最小値となる位置を探索し、動きベクトルV５７を計測する（工程５）。
0005171610.tif 000006
【００３０】
なお、（式２）の｜｜は絶対値を表す。デシメーション処理後の画像を用いた動きベクトルの計測には、±１ピクセル分の誤差が伴う。この誤差を解消するため、計測領域５１を最初の位置から動きベクトルV５７だけ移動させた計測領域５６と、その位置の周囲１−２ピクセル広い領域をもつ探索領域５５を、画像fn及びfn-1から各々再設定する（工程６）。再設定した計測領域５６及び探索領域５５を用いて、工程５と同様の手法により動きベクトルV2を再計測する（工程７）。以上の処理により、結果的に加算処理の際に補正すべき動きベクトルは、動きベクトルV５７を用いて（（2×V）＋V2）となる。【００３１】
加算処理は、取得画像に対して、その１フレーム前までの画像を用いて構成された累積加算画像を変形させて加算する処理である。加算処理の際に基準となるのは取得画像であるため、
0005171610.tif 000007
動きベクトルの計測が終わると、（n-1）フレームまでの
0005171610.tif 000008
重み係数が乗された後、画像加算部１３に入力され、動きベクトル計測部で計測された動きベクトル(（2×V）＋V2)を用いた変形処理がなされる。
0005171610.tif 000009
（α，β）の重み係数を用いた（式３）に示す処理形式でなされ、新たな累積加算画像が構成される。
0005171610.tif 000010
【００３２】
（式３）を展開し、累積加算画像を構成する各フレームの係数を示すと（式４）になる。
0005171610.tif 000011
【００３３】
（式４）に従って求まる各フレームの重み係数は加算効果を判断する上で重要なパラメータとなる。（α，β）＝（0.95，0.05）,（0.90，0.1）,（0.85，0.15）,（0.8，0.2）、n＝80の時の各フレームの重み係数を図７に示す。α値が小さい場合には（例えばα=0.80）、重み係数は70−80フレームにかけて急速に大きくなるため、輪郭強調やスペックル除去といった加算効果は低い。しかし、70フレーム以前のフレームに関しては重み係数は非常に低いため、過去に生じた計測誤差の残存は少なくなる。一方、α値が大きい場合は（例えばα=0.95）、α値が小さい場合に比べてグラフは水平に近づき、各フレームの重み係数は同等になる。多くのフレームが累積加算画像に影響を与えるため、加算効果は大きいが、過去の１フレームで生じた計測誤差が後の累積加算画像に残存することになる。最適な輪郭強調やスペックル除去の程度は術者又は検査対象によって異なる。したがって、（α，β）は初期値を（α，β）＝（0.85，0.15）として、術者が自由に変更できるようになっている。この初期値は術者が自由に変更できる。操作中のα値の変更の仕方は、図８に示すように、探触子８１に付属するダイヤル８２または診断機８３に設置されたダイヤル８４を回転させる方法が考えられる。その際、診断機の画面上には、その時のα値が表示される。なお、上記重み係数は、累積加算画像及び取得画像の双方に重み係数を乗して加算効果を調節することを目的としており、（式３）に示す形に限定されるものではなく、重み係数の数や乗数を係数として用いるといった様々な係数の形を含む。
【００３４】
累積加算画像に対する重みを制御し、誤差の残存を制御する方法は、いくつか考えられる。第一の方法は、術者が手動でダイヤルを調整する方法である。関心領域を探索するために探触子を大きく動かす際には、α値を低く設定して通常のB-mode画像に近い画像を表示し、関心領域を精査する時にα値を高くする。また、誤差の残存を消去したい時も、一時的にα値を低くすることで計測誤差によるアーチファクトを低減させることができる。第二の方法は、計測する動きベクトルの大きさに応じて自動的にα値を制御する方法である。（式１）、（式２）のｃが予め設定する閾値よりも大きい場合には、自動的にαが低くなり誤差の残存が低減される。つまり第一の方法を自動化した方法である。第三の方法は診断機または探触子に、加算画像をメモリから消去させ、再度１フレーム目から加算処理を行なうためのリフレッシュボタン(消去ボタン)を設ける方法である。これにより術者は任意に誤差の残存を０にすることができる。リフレッシュボタンが押されると、加算処理の際にα=0が入力される。リフレッシュボタンは図８に示すα値を変更するためのダイヤルと別に設けても良いが、装置の複雑化や利便性を考えると、図９に示すようにα値を変更するために探触子８１または診断機８３に備えられたダイヤル９１またはダイヤル９２を押し込むことで画像がリフレッシュされる形態が最も操作しやすい。リフレッシュボタンにより１フレーム目から加算処理をすると、スペックルが急に顕著になり、術者によっては観察しにくいことも考えられる。したがって、リフレッシュボタンによるα値は、術者が自由に設定できるものとする。【００３５】
画像の表示形態は、加算画像を画面全体に表示する形態の他に、図１０に示すようにB-mode画像１０１と加算画像１０２を並べて表示するバイプレーン表示も可能である。また、上記実施例の説明では超音波画像としてB-mode画像を用いたが、実施例１に記す超音波診断機は用いる画像を特定せず、高調波成分を選択的に画像化したTHI（Tissue Harmonic Image）画像や血流を画像化したCFM（Color Flow Mapping）画像を用いてもよい。
【実施例２】
【００３６】
本発明の第二の実施例について説明する。第二の実施例は、第一の実施例に示す動き計測及び加算処理をしている際に、別のメモリで一定枚数で構成される加算画像を構成し、術者がリフレッシュ機能を使用した時に、一定枚数の加算状態から累積加算画像の表示が可能な特徴を持つ。【００３７】
装置構成を図１１に示すブロック図を用いて説明する。B-mode画像の構成、動きベクトルの計測、画像の加算処理といった加算画像の構成に至る処理は実施例１と同じである。第二の実施例では、図１に示した第一の実施例の装置構成に加えて、リフレッシュ用画像保持部１７、動きベクトル保持部１８、画像減算部１９、リフレッシュ用加算画像保持部２０を有している。【００３８】
リフレッシュ用加算画像を構成するフレーム数は搭載可能なメモリ数によるが、ここでは５枚を例として説明する。【００３９】
実施例１と同様に画像保持部１１には、取り込まれた最新の画像よりも１フレーム前の画像が動きベクトル計測用として保持されているが、実施例２では、さらに２フレームから５フレーム前の計４枚の画像がリフレッシュ用画像保持部１７に記憶されている。【００４０】
画像の取り込みが開始され、動きベクトルの計測及び加算処理が進み、
0005171610.tif 000012
画像減算部１９を素通りしてリフレッシュ用加算画像保持部２０に記憶される。尚、動きベクトル計測部１０で計測される各画像の動きベクトル計測結果は、動きベクトル保持部１８に記憶されている。また、この段階でB-mode画像メモリに格納されている画像データは、画像保持部１１にはf5、リフレッシュ用画像保持部１７にはf4、f3、f2、f1が格納されている。【００４１】
次に、新たに画像f6が動きベクトル計測部１０を経て画像加算部１３に取り込まれ、
0005171610.tif 000013
加算画像保持部１４に記憶されると同時に、画像減算部１９に入力される。この時、リフレッシュ用画像保持部１７では画像保持部１１から画像f5が入力されると同時に、画像f1を画像減算部１９に出力する。画像減算部１９では、
0005171610.tif 000014
減算するのは画像f1そのものではなく、画像f1が加算処理されてから画像f6が加算されるまでに、画像f1に対してなされた変形処理を含めた画像f1の情報である。動きベクトル保持部１８に記憶されているのは、画像間に生じた位置ずれの情報であるため、これらを加算すると特定の画像に対してなされた変形の履歴となる。したがって、画像f2から画像f6に対して動きベクトル計測部１０で計測された動きベクトルを全て加算し、その加算結果に基づいて画像f1を変形させると、
0005171610.tif 000015
２フレームから６フレームの計５枚で構成されるリフレッシュ用の
0005171610.tif 000016
に記憶される。次の過程で減算処理されるのは画像f2であるため、画像f1と画像f2の間で求めた動きベクトルは動きベクトル保持部１８から消去される。【００４２】
以上の工程により、リフレッシュ用加算画像保持部２０には、常に取得画像から５フレームの加算画像が保持されている。そしてリフレッシュ機能が作動すると、加算画像保持部１４に保持されている累積加算画像の代わりにリフレッシュ用加算画像が表示部１５に表示され、新たな累積加算画像として加算画像保持部１４に保持される。
【実施例３】
【００４３】
本発明の第三の実施例について説明する。第三の実施例は、取り込んだ画像を累積加算してなる加算画像を複数のユニットに分割してユニット画像を構成し、各ユニット画像に対して重み係数を付けることで、加算効果が高くかつ誤差の残存が低い加算画像を表示することを特徴とする。【００４４】
装置構成を図１２に示す。術者が予め設定する加算枚数を３０枚、ユニットを構成する画像枚数を５枚として、以下に処理工程について説明する。【００４５】
超音波画像構成ユニットから入力された画像fn、及びその１フレーム前の画像fn-1を用いて動きベクトルを計測する工程は実施例１及び実施例２に記載の工程と同じである。但し、実施例３においては、画像加算部１３で５枚分の加算画像が構成されると、ユニット画像としてユニット画像保持部２０３に記憶される。つまり、例えば２６フレーム目の画像f26の加算処理を終えた時点では、ユニット画像保持部２０３には、
0005171610.tif 000017
保持されている（図１３）。【００４６】
画像ｆ26が画像加算部１３に入力されると、ユニット画像保持部２０３に記憶されている５枚のユニット画像は画像変形部２０４に取り込まれ、動きベクトル計測部１０で計測された動きベクトルに基づく変形処理がなされる。変形処理がなされた画像は重み付け部１６に入力され、取得画像ｆ26と共に重み係数が乗され、画像加算部１３にて（式５）に示す加算処理が行われる。加算処理は５枚のユニット画像及び取り込まれたｆ26の各々に所定の重み係数を乗して行なわれる。例えば、加算画像に対する各ユニット画像の重み係数の分布は図１４のグラフになる。
0005171610.tif 000018
【００４７】
３１フレーム目の画像f31が取り込まれると、
0005171610.tif 000019
加算処理に用いられた各ユニット画像は、画像加算部１３にて変形処理がなされた状態で再びユニット画像保持部２０３に記憶される。【００４８】
以上に説明したように、各ユニットに重み係数を乗することにより、加算効果を維持したまま、誤差の残存を自動的に低減することができる。また、動きベクトルの計測結果に応じて、例えば計測誤差が大きなユニットに対する重み係数を低くするなど、特定のユニットに対する重み係数を自動的に制御することも可能である。【００４９】
実施例の説明では加算枚数を３０枚、ユニットを構成する画像枚数を５枚として説明したが、実施例１に記したように検査対象の動きの大きさに応じて重み係数を調整することで、加算効果や誤差の残存を自由に制御することができる。また、ユニットを構成する画像枚数を少なくすることで、重み係数をフレーム単位で制御することが可能になり、加算効果をより自由に制御できるだけでなく、誤差が大きい特定のフレームを加算処理から外すことも可能になる。
【図１】
0005171610.tif 000020
【図２】
0005171610.tif 000021
【図３】
0005171610.tif 000022
【図４】
0005171610.tif 000023
【図５】
0005171610.tif 000024
【図６】
0005171610.tif 000025
【図７】
0005171610.tif 000026
【図８】
0005171610.tif 000027
【図９】
0005171610.tif 000028
【図１０】
0005171610.tif 000029
【図１１】
0005171610.tif 000030
【図１２】
0005171610.tif 000031
【図１３】
0005171610.tif 000032
【図１４】
0005171610.tif 000033