JP2014195512A - 超音波診断装置 - Google Patents

超音波診断装置 Download PDF

Info

Publication number
JP2014195512A
JP2014195512A JP2013071825A JP2013071825A JP2014195512A JP 2014195512 A JP2014195512 A JP 2014195512A JP 2013071825 A JP2013071825 A JP 2013071825A JP 2013071825 A JP2013071825 A JP 2013071825A JP 2014195512 A JP2014195512 A JP 2014195512A
Authority
JP
Japan
Prior art keywords
image frame
frame
luminance value
ultrasonic
image
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2013071825A
Other languages
English (en)
Other versions
JP5956375B2 (ja
Inventor
雅史 野口
Masafumi Noguchi
雅史 野口
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujifilm Corp
Original Assignee
Fujifilm Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fujifilm Corp filed Critical Fujifilm Corp
Priority to JP2013071825A priority Critical patent/JP5956375B2/ja
Publication of JP2014195512A publication Critical patent/JP2014195512A/ja
Application granted granted Critical
Publication of JP5956375B2 publication Critical patent/JP5956375B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

【課題】超音波画像のブレを抑制しつつ同一組織について輝度値をほぼ均一に保つことができる超音波診断装置を提供する。
【解決手段】画像生成部10で生成された所定の画像フレームの輝度値から過去の画像フレームの輝度値を差し引いて画素毎に輝度値の差分ΔIを算出したフレーム相関処理部7は、輝度値の差分ΔIが予め定められた負の閾値T2より小さくなるほど、または、輝度値の差分ΔIが予め定められた正の閾値T1より大きくなるほど所定の画像フレームの重みαを大きく変化させ、且つ、所定の画像フレームの重みαの原点からの変化量が負の閾値T2より小さな領域L3と比べて正の閾値T1より大きな領域L1において小さくなるように所定の画像フレームの重みαを設定し、所定の画像フレームと過去の画像フレームを重み付け加算する。
【選択図】図1

Description

この発明は、超音波診断装置に係り、特に、超音波画像にフレーム相関処理を施す超音波診断装置に関する。
従来から、医療分野において、超音波画像を利用した超音波診断装置が実用化されている。通常、この種の超音波診断装置は、超音波探触子から被検体内に向けて超音波ビームを送信し、被検体からの超音波エコーを超音波探触子で受信して、その受信信号を電気的に処理することにより超音波画像が生成される。
一般に、超音波画像を生成する際には、ノイズを低減すると共に画像を滑らかにするために、フレーム相関処理が行われている。例えば、IIRフィルタ(Infinite Impulse Response Filter)やFIRフィルタ(Finite Impulse Response Filter)などを用いて、現在の画像フレームに対して過去の画像フレームを重み付け加算することによりフレーム相関処理が行われる。
ここで、超音波画像の画像フレーム間で対象が変位している場合、例えば、操作者が超音波探触子を移動させている場合や、心臓などの動きを有する対象を観察する場合には、現在の画像フレームと過去の画像フレームとの間で対象が異なる位置に存在するため、フレーム相関処理を施した超音波画像にブレが生じるといった問題があった。図7に、画像フレーム間で対象が変位している超音波画像にフレーム相関処理を施した一例を示す。このように、超音波画像の全体にブレが生じ、特に、血管内腔Sなどの低輝度な組織では、過去の画像フレームにおいて血管壁Wなどの高輝度な組織が存在していた部分Pにブレが顕著に表れてしまう。
そこで、例えば、特許文献1には、画像フレーム間における対象の変位に起因したブレを抑制しつつフレーム相関処理を行う超音波診断装置が開示されている。この超音波診断装置では、図8に示すように、現在の画像フレームと過去の画像フレームとの輝度値の差ΔIを画素毎に算出し、その輝度値の差ΔIの絶対値が予め定められた閾値Tよりも大きい部分については、画像フレーム間における対象の変位に起因した輝度の変化であると判断して、現在の画像フレームの重みαを急激に正側に大きく変化させる。一方、輝度値の差ΔIの絶対値が閾値Tよりも小さい部分については、ノイズに起因した輝度値の変化であると判断して、過去の画像フレームに基づいた通常の平滑化処理が施される。これにより、過去の画像フレームとの間で対象の変位が生じている部分には現在の画像フレームの輝度を大きく反映させると共に、ノイズが生じている部分には過去の画像フレームとの間で平滑化処理が施されるため、図9に示すように、超音波画像のブレを抑制しつつノイズが除去された超音波画像を生成することができる。
特開2002−34984号公報
しかしながら、画像フレーム間の輝度値の差ΔIにおいて正負の符号が考慮されていないため、図9のように血管壁Wの周囲に血管壁Wと同等の輝度値を有する甲状腺Gなどの組織が存在する場合に、輝度値の差ΔIが負になる部分については明瞭な超音波画像を生成することができるものの、輝度値の差ΔIが正になる部分において輝度が不均一になる部分Rが生じてしまうといった問題があった。
すなわち、血管壁Wや甲状腺Gが存在していた位置に血管内腔Sが移動した場合(輝度値の差ΔIが負になる場合)には、輝度値の差ΔIの絶対値はいずれも閾値Tより大きくなるため、フレーム相関処理では現在の画像フレームにおける血管内腔Sの輝度が一様に大きく反映される。これにより、血管内腔Sにおける画像のブレを抑制することができる。
一方、図10に示すように、血管内腔Sと甲状腺Gの両者が存在していた部分にわたって血管壁Wが移動した場合(輝度値の差ΔIが正になる場合)には、血管壁Wと血管内腔Sとの輝度値の差は大きいのに対し、血管壁Wと甲状腺Gとの輝度値の差は小さいため、血管壁Wの輝度値との差分を算出すると、輝度値の差ΔIの絶対値が閾値Tより大きい部分Waと閾値Tより小さい部分Wbがそれぞれ表れることになる。このため、血管内腔Sが位置していた部分Waについては現在の画像フレームにおける血管壁Wの輝度が大きく反映されるのに対し、甲状腺Gが位置していた部分Wbについては過去の画像フレームとの間で平滑化処理が施されるため、フレーム相関処理後の超音波画像には、血管壁Wの輝度値が不均一な部分R、いわゆるぎらつきが発生してしまう。
この発明は、このような従来の問題点を解消するためになされたもので、超音波画像のブレを抑制しつつ同一組織の輝度値をほぼ均一に保つことができる超音波診断装置を提供することを目的とする。
この発明に係る超音波診断装置は、超音波探触子から被検体内の対象に向けて超音波ビームの送受信を行い、超音波探触子から出力される受信信号に基づいて超音波画像を順次生成すると共に生成された所定の画像フレームに対して過去の画像フレームに基づく画像処理を施す超音波診断装置であって、受信信号に基づいて超音波画像をフレーム毎に順次生成する画像生成部と、画像生成部で生成された所定の画像フレームの輝度値から過去の画像フレームの輝度値を差し引いて画素毎に輝度値の差分を算出し、輝度値の差分に応じて過去の画像フレームに対する所定の画像フレームの重みをそれぞれ設定して、所定の画像フレームと過去の画像フレームを重み付け加算するフレーム相関処理部とを備え、フレーム相関処理部は、輝度値の差分が予め定められた負の閾値より小さくなるほど、または、輝度値の差分が予め定められた正の閾値より大きくなるほど所定の画像フレームの重みを大きく変化させ、且つ、所定の画像フレームの重みの原点からの変化量を負の閾値より小さな領域と比べて正の閾値より大きな領域において小さくするものである。
ここで、正の閾値と負の閾値は、ノイズに起因した輝度値の変化を示す領域と、過去の画像フレームからの対象の変位に起因した輝度値の変化を示す領域との境界部分に定めることができる。
また、フレーム相関処理部は、正の閾値より大きな領域において、過去の画像フレームからの対象の変位に起因して生じる超音波画像のブレを抑制しつつ同一対象の輝度値を均一に保つように、所定の画像フレームの重みを変化させるのが好ましい。また、フレーム相関処理部は、負の閾値より小さな領域において、過去の画像フレームからの対象の変位に起因して生じる超音波画像のブレを抑制するように、所定の画像フレームの重みを変化させるのが好ましい。
さらに、フレーム相関処理部は、正の閾値と負の閾値の間の領域において、所定の画像フレームの重みを、ノイズを抑制するような値で一定に保つのが好ましい。
また、正の閾値と負の閾値は、絶対値が同じ値となるように定められ、フレーム相関処理部は、正の閾値より大きな領域と負の閾値より小さな領域において所定の画像フレームの重みを直線的に変化させ、且つ、その変化の傾きを負の閾値より小さな領域と比べて正の閾値より大きな領域において緩やかにすることができる。
また、過去の画像フレームは、所定の画像フレームとフレーム相関処理する前段階において、過去の画像フレームより以前の画像フレームとの間でフレーム相関処理を施された画像フレームを用いることができる。また、過去の画像フレームは、所定の画像フレームとフレーム相関処理する前段階において、過去の画像フレームより以前の画像フレームとの間でフレーム相関処理が施されていない画像フレームを用いることもできる。
また、過去の画像フレームは、所定の画像フレームとフレーム相関処理する前段階において、過去の画像フレームより以前の画像フレームとの間でフレーム相関処理が施された画像フレームと、過去の画像フレームより以前の画像フレームとの間でフレーム相関処理が施されていない画像フレームとを含むこともできる。
この発明に係る超音波診断装置は、超音波探触子から被検体内の対象に向けて超音波ビームの送受信を行い、超音波探触子から出力される受信信号に基づいて超音波画像を順次生成すると共に生成された所定の画像フレームに対して過去の画像フレームに基づく画像処理を施す超音波診断装置であって、受信信号に基づいて超音波画像をフレーム毎に順次生成する画像生成部と、画像生成部で生成された所定の画像フレームの輝度値から過去の画像フレームの輝度値を差し引いて画素毎に輝度値の差分を算出し、輝度値の差分に応じて過去の画像フレームに対する所定の画像フレームの重みをそれぞれ設定して、所定の画像フレームと過去の画像フレームを重み付け加算するフレーム相関処理部とを備え、フレーム相関処理部は、輝度値の差分が予め定められた負の閾値より小さくなるほど、または、輝度値の差分が予め定められた正の閾値より大きくなるほど所定の画像フレームの重みを大きく変化させ、且つ、所定の画像フレームの重みの原点からの変化量を負の閾値より小さな領域と比べて正の閾値より大きな領域において大きくするものである。
この発明によれば、所定の画像フレームと過去の画像フレームの輝度値の差分に応じて所定の画像フレームの重みを設定する際に所定の画像フレームの重みの原点からの変化量を負の閾値より小さな領域と比べて正の閾値より大きな領域において小さくするので、超音波画像のブレを抑制しつつ同一組織の輝度値をほぼ均一に保つことが可能となる。
この発明の実施の形態1に係る超音波診断装置の構成を示すブロック図である。 輝度値の差分に対する現在の画像フレームの重みと過去の画像フレームの重みの分布を示すグラフである。 画像フレーム間で血管壁が変位する様子を示す図である。 実施の形態1に係る超音波診断装置を用いて超音波画像にフレーム相関処理を施した実施例を示す図である。 実施の形態1の変形例に係る超音波診断装置において用いられる、輝度値の差分に対する現在の画像フレームの重みと過去の画像フレームの重みの分布を示すグラフである。 実施の形態2に係る超音波診断装置において用いられる、輝度値の差分に対する現在の画像フレームの重みと過去の画像フレームの重みの分布を示すグラフである。 従来例において、フレーム相関処理を施した超音波画像を示す図である。 従来例において、輝度値の差分に対する現在の画像フレームの重みと過去の画像フレームの重みの分布を示すグラフである。 図8に示す重みの分布に従ってフレーム相関処理を施した超音波画像を示す図である。 画像フレーム間で血管壁が変位する様子を示す図である。
以下、この発明の実施の形態を添付図面に基づいて説明する。
実施の形態1
図1に、実施の形態1に係る超音波診断装置の構成を示す。超音波診断装置は、超音波探触子1を備え、この超音波探触子1に送信回路2と受信回路3がそれぞれ接続されている。受信回路3には信号処理部4が接続され、この信号処理部4にDSC(Digital Scan Converter)5、画像処理部6、フレーム相関処理部7、表示制御部8および表示部9が順次接続されている。そして、信号処理部4、DSC5および画像処理部6により、画像生成部10が形成されている。
また、受信回路3と信号処理部4には受信データメモリ11が接続されると共に、画像処理部6とフレーム相関処理部7には画像メモリ12が接続されている。さらに、送信回路2、受信回路3、信号処理部4、DSC5、画像処理部6、フレーム相関処理部7および表示制御部8には制御部13が接続され、この制御部13に操作部14と格納部15がそれぞれ接続されている。
超音波探触子1は、複数のトランスデューサ素子が1次元又は2次元に配列されたアレイトランスデューサを有している。これらのトランスデューサ素子は、それぞれ送信回路2から供給される駆動信号に従って超音波を送信すると共に被検体からの超音波エコーを受信して受信信号を出力する。各トランスデューサ素子は、例えば、PZT(チタン酸ジルコン酸鉛)に代表される圧電セラミックや、PVDF(ポリフッ化ビニリデン)に代表される高分子圧電素子、PMN−PT(マグネシウムニオブ酸・チタン酸鉛固溶体)に代表される圧電単結晶等からなる圧電体の両端に電極を形成した振動子によって構成される。
そのような振動子の電極に、パルス状又は連続波の電圧を印加すると、圧電体が伸縮し、それぞれの振動子からパルス状又は連続波の超音波が発生して、それらの超音波の合成により超音波ビームが形成される。また、それぞれの振動子は、伝搬する超音波を受信することにより伸縮して電気信号を発生し、それらの電気信号は、超音波の受信信号として出力される。
送信回路2は、例えば、複数のパルサを含んでおり、制御部13からの制御信号に応じて選択された送信遅延パターンに基づいて、超音波探触子1のアレイトランスデューサから送信される超音波が超音波ビームを形成するように、それぞれの駆動信号の遅延量を調節してアレイトランスデューサに供給する。
受信回路3は、超音波探触子1のアレイトランスデューサから得られた受信信号を増幅してA/D変換した後、制御部13からの制御信号に応じて選択された受信遅延パターンに基づいて設定される音速または音速の分布に従い、各受信信号にそれぞれの遅延を与えて加算することにより、受信フォーカス処理を行う。この受信フォーカス処理により、超音波探触子1で受信した超音波エコーの焦点が絞り込まれたサンプルデータ(音線信号)が生成される。
信号処理部4は、受信回路3で生成されたサンプルデータに対し、超音波の反射位置の深度に応じて距離による減衰の補正を施した後、包絡線検波処理を施すことにより、被検体内の組織に関する断層画像情報であるBモード画像データを生成する。
DSC5は、信号処理部4で生成されたBモード画像データを通常のテレビジョン信号の走査方式に従う画像データに変換(ラスター変換)する。
画像処理部6は、DSC5から入力されるBモード画像データに階調処理等の各種の必要な画像処理を施した後、Bモード画像データをフレーム相関処理部7に出力すると共にBモード画像データをフレーム毎に画像メモリ12に格納する。
フレーム相関処理部7は、画像メモリ12に格納された過去の画像フレームのBモード画像データを読み出して、現在の画像フレームの輝度値から過去の画像フレームの輝度値を差し引いて画素毎に輝度値の差分を算出する。そして、算出された輝度値の差分に応じて過去の画像フレームに対する現在の画像フレームの重みを画素毎に設定し、現在の画像フレームと過去の画像フレームを重み付け加算することにより、フレーム相関処理を施したBモード画像データを生成する。
表示制御部8は、フレーム相関処理部7から入力されたBモード画像データに基づいて、表示部9に超音波診断画像を表示させる。
表示部9は、例えば、LCD等のディスプレイ装置を含んでおり、表示制御部8の制御の下で、超音波診断画像を表示する。
操作部14は、操作者が入力操作を行うためのもので、キーボード、マウス、トラックボール、タッチパネル等から形成することができる。
格納部15は、動作プログラム等を格納するもので、ハードディスク、フレキシブルディスク、MO、MT、RAM、CD−ROM、DVD−ROM、SDカード、CFカード、USBメモリ等の記録メディア、またはサーバ等を用いることができる。
制御部13は、操作者により操作部14から入力された各種の指令信号等に基づいて、超音波診断装置内の各部の制御を行うものである。
なお、信号処理部4、DSC5、画像処理部6、フレーム相関処理部7および表示制御部8は、CPUと、CPUに各種の処理を行わせるための動作プログラムから構成されるが、それらをデジタル回路で構成してもよい。
次に、フレーム相関処理部7において超音波画像に施されるフレーム相関処理について詳細に説明する。
フレーム相関処理は、現在の画像フレームを過去の画像フレームとの相関に基づいて平滑化処理することにより、フレーム間で不均一に発生する熱雑音などのノイズを低減するために行われ、例えば、現在の画像フレームと1つ前の過去の画像フレームとの間でフレーム相関処理を施して超音波画像(相関処理画像)Pを生成する場合には、P=α×F+β×Fn−1で表すことができる。なお、Fは現在の画像フレーム、Fn−1は1つ前の画像フレームをそれぞれ示している。また、αおよびβはそれぞれ相関係数(重み)を示し、α+β=1で表すことができる。
ここで、相関係数αおよびβの値は、図2に示すように、現在の画像フレームの輝度値から過去の画像フレームの輝度値を差し引いて画素毎に算出された輝度値の差分ΔIに応じて変化されるものである。具体的には、輝度値の差分ΔIに対して、正の側と負の側にそれぞれ正の閾値T1と負の閾値T2が予め定められている。この正の閾値T1と負の閾値T2は、フレーム相関処理の目的に応じて定められるもので、例えば、熱雑音などのノイズに起因した輝度値の変化を示す領域と、過去の画像フレームからの対象の変位に起因した輝度値の変化を示す領域との境界部分に定めることができる。すなわち、ノイズに起因して生じた輝度値の差は、正の閾値T1と負の閾値T2で挟まれた領域L1に含まれ、対象の変位に起因して生じた輝度値の差は、正の閾値T1より大きな領域L2と負の閾値T2より小さな領域L3にそれぞれ含まれることになる。ここでは、正の閾値T1と負の閾値T2は、絶対値が同じ値となるように定められているものとする。
フレーム相関処理部7は、正の閾値T1と負の閾値T2の間の領域L1においては、現在の画像フレームの重みαと過去の画像フレームの重みβを、ノイズを抑制するような値で一定に保っている。これにより、輝度値の差分ΔIが領域L1に位置する画素については、過去の画像フレームの輝度値に基づく平滑化処理が施され、ノイズを除去することができる。
一方、フレーム相関処理部7は、輝度値の差分ΔIが正の閾値T1より大きな領域L2と、負の閾値T2より小さな領域L3において、輝度値の差分ΔIが負の閾値T2より小さくなるほど現在の画像フレームの重みαが急激に大きくなるように直線的に変化させると共に、輝度の差分ΔIが正の閾値T1より大きくなるほど現在の画像フレームの重みαが大きくなるように直線的に変化させる。この時、現在の画像フレームの重みαを変化させる傾きは、負の閾値T2より小さな領域L3と比べて、正の閾値T1より大きな領域L2において穏やかとなるように設定されている。すなわち、現在の画像フレームの重みαの原点Oからの変化量は、負の閾値T2より小さな領域L3と比べて、正の閾値T1より大きな領域L2において小さくなるように設定されている。ここで、原点Oは、輝度値の差ΔIの座標軸と重みαおよびβの座標軸が互いにゼロを示す点である。
ここで、フレーム相関処理部7は、負の閾値T2より小さな領域L3において、過去の画像フレームからの対象の変位に起因して生じる超音波画像のブレを抑制するように、現在の画像フレームの重みαを急激に大きくなるように変化させている。なお、領域L3における重みαの変化量(傾き)と重みβの変化量は、画像フレーム間における対象の移動速度、すなわちフレームレートや被検体内における対象の移動速度などに基づいて設定することができる。
また、フレーム相関処理部7は、正の閾値T1より大きな領域L2において、過去の画像フレームからの対象の変位に起因して生じる超音波画像のブレを抑制しつつ同一対象の輝度を均一に保つように、現在の画像フレームの重みαを領域L3に対して緩やかに変化させている。なお、領域L2における重みαの変化量と重みβの変化量は、画像フレーム間における対象の移動速度および超音波画像における輝度の分布などに基づいて設定することができる。
なお、正の閾値T1と負の閾値T2は、ノイズに起因した輝度値の変化を示す領域と、過去の画像フレームからの対象の変位に起因した輝度値の変化を示す領域とを明確に分ける境界部分に設定する必要はなく、ノイズに起因した輝度値の変化を示す領域内、あるいは、過去の画像フレームからの対象の変位に起因した輝度値の変化を示す領域内にある程度侵すように定めてもよい。例えば、甲状腺Gが存在していた位置に血管壁Wが移動した際に生じる輝度値の差ΔIは領域L1内に存在し、且つ、血管内腔Sが存在していた位置に血管壁Wが移動した際に生じる輝度値の差ΔIは領域L2に存在するように、正の閾値T1を定めることができる。
次に、この実施の形態1の動作について説明する。
まず、送信回路2からの駆動信号に従って、超音波探触子1から被検体内の対象に向けて超音波ビームが発せられる。超音波ビームは、被検体内を伝搬して被検体内の対象により超音波エコーとして反射され、この超音波エコーが超音波探触子1で受信される。超音波エコーが超音波探触子1で受信されると、超音波探触子1から受信信号が受信回路3に出力され、受信回路3は入力された受信信号に基づいてサンプルデータを生成する。続いて、サンプルデータが画像生成部10に出力されて、画像生成部10がサンプルデータに基づいて画像データを順次生成する。そして、生成された画像データは、画像処理部10により画像メモリ12にフレーム毎に格納されると共にフレーム相関処理部7に出力される。
フレーム相関処理部7は、画像メモリ12に格納された画像データを読み出し、現在の画像フレームの輝度値から過去の画像フレームの輝度値を差し引いた輝度値の差分ΔIを画素毎に算出する。そして、算出された輝度値の差分ΔIに応じて現在の画像フレームの重みαと過去の画像フレームの重みβをそれぞれ設定して、現在の画像フレームと過去の画像フレームを重み付け加算することにより、超音波画像にフレーム相関処理を施す。
例えば、過去の画像フレームと現在の画像フレームとの間で対象がほとんど変位していない部分、あるいは、図3に示すように甲状腺Gが存在していた位置に同様の輝度値を有する血管壁Wが移動した部分R1など、現在の画像フレームと過去の画像フレームとの間で輝度値の差ΔIが小さい部分(例えば、輝度値の差ΔIが−20dBから20dBの範囲内となる部分)については、正の閾値T1と負の閾値T2の間の領域L1で設定された重みαおよびβによりフレーム相関処理が施される。これにより、画像フレーム間で発生したノイズを除くことができる。
また、現在の画像フレームと過去の画像フレームとの間で輝度値の差ΔIが負の閾値T2より小さな領域L3に含まれる部分、例えば、高輝度な甲状腺Gが存在していた位置に低輝度な血管内腔Sが移動した部分R2、あるいは、高輝度な血管壁Wが存在していた位置に低輝度な血管内腔Sが移動した部分R3など(例えば、輝度値の差ΔIが−20dBより小さい部分)は、領域L3において設定された重みαおよびβによりフレーム相関処理が施される。このとき、領域L3では、重みαが急激に大きく変化されると共に重みβが急激に小さく変化されているため、現在の画像フレームの輝度が大きく反映される。すなわち、血管内腔Sなどの低輝度な部分が移動した画素については、その血管内腔Sの輝度がそのまま反映される。このため、フレーム相関処理後の血管内腔Sの輝度が、過去の画像フレームの甲状腺Gや血管壁Wの輝度で損なわれることはなく、過去の画像フレームからの対象の変位に起因して生じる画像のブレを抑制することができる。
さらに、現在の画像フレームと過去の画像フレームとの間で輝度値の差ΔIが正の閾値T1より大きな領域L2に含まれる部分、例えば、血管内腔Sが存在していた位置に血管壁Wが移動した部分R4、あるいは、血管内腔Sが存在していた位置に甲状腺Gが移動した部分R5など(例えば、輝度値の差ΔIが20dBより大きい部分)は、領域L2において設定された重みαおよびβによりフレーム相関処理が施される。このとき、領域L2では、領域L3と比べて、重みαが正側に緩やかに変化されると共に重みβが負側に緩やかに変化されているため、過去の画像フレームにおいて互いに異なる輝度を有する部分、例えば低輝度な血管内腔Sと高輝度な甲状腺Gをある程度等しい重みβで反映することができる。このため、血管内腔Sと甲状腺Gが存在していた部位にわたって高輝度な血管壁Wが移動した場合でも、血管壁Wを全体的にある程度等しい重みでフレーム相関処理することができ、高輝度な同一対象の輝度を均一に保つことができる。また、血管壁Wなどの高輝度な部分は画像のブレが発生しにくいため、領域L2のように重みαおよびβの変化量を緩やかに設定しても充分にブレを除くことができる。
このようにして、フレーム相関処理が施された超音波画像は、フレーム相関処理部7から表示制御部8に出力され、表示部9に表示される。
本実施の形態によれば、現在の画像フレームと過去の画像フレームの輝度値の差分に応じた重みで画素毎にフレーム相関処理を行うため、ノイズを除去すると共に、超音波画像のブレを抑制しつつ同一対象の輝度を均一に保つことができる。このため、例えば、操作者が超音波探触子を移動させている場合でも、超音波画像全体のブレを抑制しつつ、血管壁Wや甲状腺Gにおける輝度を一定に保つことができ、超音波画像の観察が輝度のぎらつきによって阻害されるのを抑制することができる。
実際に、図2に示したグラフに従って、フレーム相関処理を施した実施例を図4に示す。図8に示したグラフに従って、フレーム相関処理を施した図9の比較例では、血管壁Wにおいて輝度が不均一な領域Rが生じているのに対し、図4の実施例では、その領域Rの輝度が均一化されており、血管壁Wの輝度が全体的にほぼ一定に保たれていることがわかる。また、図4では、画像フレーム間における対象の変位に起因して生じる画像のブレが抑制されており、特に、血管内腔Sではブレがほとんど観察されない。さらに、熱雑音などに起因するノイズも除去されていることがわかる。このことから、本発明により、ノイズを除去すると共に、画像のブレを抑制しつつ同一対象の輝度を均一に保つことができた。
なお、上記の実施の形態では、正の閾値T1と負の閾値T2を絶対値が同じ値となるように設定すると共に現在の画像フレームの重みαを変化させる傾きを領域L3に比べて領域L2において緩やかになるように設定したが、現在の画像フレームの重みαの原点Oからの変化量が領域L3に比べて領域L2において小さければよく、これに限られるものではない。
例えば、図5に示すように、正の閾値T1に対して負の閾値T2を原点Oに近づけると共に、領域L2と領域L3における重みαおよびβを同じ変化量(傾き)で変化させることができる。このように、負の閾値T2を原点Oに近づけることにより、重みαおよびβを変化させる開始位置を領域L2と比べて領域L3において原点Oに近い位置に設定することができる。このため、領域L2と領域L3において重みαおよびβを同じ変化量で変化させても、重みαおよびβの原点Oからの変化量は、領域L3に比べて領域L2において小さくすることができ、上記の実施の形態と同様の効果を得ることができる。
また、上記の実施の形態では、現在の画像フレームと1つ前の過去の画像フレームとの間でフレーム相関処理を施したが、現在の画像フレームとの間でフレーム相関処理を施すことができればよく、これに限られるものではない。例えば、現在の画像フレームと複数の過去の画像フレームとの間でフレーム相関処理を施すことができる。また、現在の画像フレームと、数フレーム離れた過去の画像フレーム(例えば、3つ前の過去の画像フレーム)との間でフレーム相関処理を施すこともできる。
また、過去の画像フレームとしては、現在の画像フレームとフレーム相関処理する前段階において、過去の画像フレームより以前の画像フレームとの間で既にフレーム相関処理が施されているものを用いることができ、例えばIIRフィルタを用いて現在の画像フレームとの間でフレーム相関処理を行うことができる。また、過去の画像フレームは、現在の画像フレームとフレーム相関処理する前段階において、過去の画像フレームより以前の画像フレームとの間でフレーム相関処理が施されていないものを用いることができ、例えばFIRフィルタを用いて現在の画像フレームとの間でフレーム相関処理を行うことができる。さらに、過去の画像フレームは、現在の画像フレームとフレーム相関処理する前段階において、過去の画像フレームより以前の画像フレームとの間でフレーム相関処理が施された画像フレームと、過去の画像フレームより以前の画像フレームとの間でフレーム相関処理が施されていない画像フレームの両者を含むものを用いることもできる。
また、現在の画像フレームとしては、画像メモリ12に格納された所定の画像フレームを用いることができ、この所定の画像フレームより以前の画像フレームが過去の画像フレームとして用いられる。
また、上記の実施の形態では、フレーム相関処理部7は、領域L2およびL3において、現在の画像フレームの重みαを直線的に変化させていたが、輝度値の差分ΔIが負の閾値T2より小さくなるほど重みαを大きく変化させると共に輝度値の差分ΔIが正の閾値T1より大きくなるほど重みαを大きく変化させ、且つ、重みαの原点Oからの変化量を領域L3と比べて領域L2において小さくできればよく、これに限られるものではない。例えば、重みαを段階的に大きく変化させてもよく、重みαを曲線的に大きく変化させてもよい。
実施の形態2
実施の形態1では、同一対象の輝度を均一に保つために、現在の画像フレームの重みαの原点Oからの変化量が領域L3と比べて領域L2において小さくなるように重みαを変化させていたが、重みαの原点Oからの変化量が領域L3と比べて領域L2において大きくなるように重みαを変化させることにより、同一対象の輝度を高めて強調することもできる。
これにより、心臓などの動きの速い組織を観察する超音波画像において、ノイズを除去すると共に、弁や心壁などの高輝度な部分のブレを抑制しつつその位置を明瞭に観察することができる。
例えば、図6に示すように、正の閾値T1と負の閾値T2が絶対値が同じ値となるようにそれぞれ定められ、現在の画像フレームの重みαと過去の画像フレームの重みβが、正の閾値T1と負の閾値T2の間の領域L1において、ノイズを抑制するような値で一定に保たれている。これにより、輝度値の差分ΔIが領域L1に位置する画素については、過去の画像フレームの輝度値に基づいて平滑化処理を施して、ノイズを除くことができる。
また、輝度値の差分ΔIが正の閾値T1より大きな領域L2と、負の閾値T2より小さな領域L3において、輝度値の差分ΔIが負の閾値T2より小さくなるほど現在の画像フレームの重みαが大きくなるように直線的に変化させると共に、輝度の差分ΔIが正の閾値T1より大きくなるほど現在の画像フレームの重みαが急激に大きくなるように直線的に変化させる。この時、現在の画像フレームの重みαを変化させる傾きは、領域L3と比べて領域L2において急になるように設定されている。すなわち、現在の画像フレームの重みαの原点Oからの変化量は、領域L3と比べて領域L2において大きくなるように設定されている。
ここで、負の閾値T2より小さな領域L3では、過去の画像フレームからの対象の変位に起因して生じる超音波画像のブレを抑制するように、現在の画像フレームの重みαを大きく変化させている。これにより、例えば、心壁などの高輝度な部分が存在していた位置に、心臓の内腔などの低輝度な部分が移動した画素については、現在の画像フレームの輝度を大きく反映させることができ、フレーム相関処理後の低輝度な部分が、過去の画像フレームの高輝度な部分に損なわれず、過去の画像フレームからの対象の変位に起因して生じる画像のブレを抑制することができる。
また、正の閾値T1より大きな領域L2では、過去の画像フレームからの対象の変位に起因して生じる超音波画像のブレを抑制しつつ同一対象の輝度を強調するように、現在の画像フレームの重みαを領域L3に対して急激に大きく変化させている。これにより、例えば、心臓の内腔などの低輝度な部分が存在していた位置に、弁や心壁などの高輝度な部分が移動した場合でも、その画素について現在の画像フレームの輝度(弁や心壁などの輝度)を大きく反映させることができ、フレーム相関処理後の高輝度な部分は、過去の画像フレームの低輝度な部分との間でほとんど平滑化されず、その輝度を保つことができる。これにより、過去の画像フレームからの対象の変位に起因して生じる画像のブレを抑制しつつその対象を強調することができる。
本実施の形態によれば、現在の画像フレームと過去の画像フレームの輝度値の差分に応じた重みで画素毎にフレーム相関処理を行うことにより、ノイズを除去すると共に、超音波画像のブレを抑制しつつ同一対象を強調することができる。このため、例えば、心臓などの動きが速い組織を観察する場合でも、心臓の内腔に発生するブレを抑制しつつ、弁や心壁の輝度を強調することができ、動きの速い組織を詳細に観察することができる。
1 超音波探触子、2 送信回路、3 受信回路、4 信号処理部、5 DSC、6 画像処理部、7 フレーム相関処理部、8 表示制御部、9 表示部、10 画像生成部、11 受信データメモリ、12 画像メモリ、13 制御部、14 操作部、15 格納部、α 現在の画像フレームの重み、β 過去の画像フレームの重み、ΔI 輝度値の差分、T1 正の閾値、T2 負の閾値、L1〜L3 領域、W 血管壁、S 血管内腔、G 甲状腺。

Claims (10)

  1. 超音波探触子から被検体内の対象に向けて超音波ビームの送受信を行い、前記超音波探触子から出力される受信信号に基づいて超音波画像を順次生成すると共に生成された所定の画像フレームに対して過去の画像フレームに基づく画像処理を施す超音波診断装置であって、
    前記受信信号に基づいて超音波画像をフレーム毎に順次生成する画像生成部と、
    前記画像生成部で生成された所定の画像フレームの輝度値から過去の画像フレームの輝度値を差し引いて画素毎に輝度値の差分を算出し、前記輝度値の差分に応じて前記過去の画像フレームに対する前記所定の画像フレームの重みをそれぞれ設定して、前記所定の画像フレームと前記過去の画像フレームを重み付け加算するフレーム相関処理部と
    を備え、
    前記フレーム相関処理部は、前記輝度値の差分が予め定められた負の閾値より小さくなるほど、または、前記輝度値の差分が予め定められた正の閾値より大きくなるほど前記所定の画像フレームの重みを大きく変化させ、且つ、前記所定の画像フレームの重みの原点からの変化量を前記負の閾値より小さな領域と比べて前記正の閾値より大きな領域において小さくすることを特徴とする超音波診断装置。
  2. 前記正の閾値と前記負の閾値は、ノイズに起因した輝度値の変化を示す領域と、前記過去の画像フレームからの対象の変位に起因した輝度値の変化を示す領域との境界部分に定められる請求項1に記載の超音波診断装置。
  3. 前記フレーム相関処理部は、前記正の閾値より大きな領域において、前記過去の画像フレームからの対象の変位に起因して生じる超音波画像のブレを抑制しつつ同一対象の輝度値を均一に保つように、前記所定の画像フレームの重みを変化させる請求項1または2に記載の超音波診断装置。
  4. 前記フレーム相関処理部は、前記負の閾値より小さな領域において、前記過去の画像フレームからの対象の変位に起因して生じる超音波画像のブレを抑制するように、前記所定の画像フレームの重みを変化させる請求項1〜3のいずれか一項に記載の超音波診断装置。
  5. 前記フレーム相関処理部は、前記正の閾値と前記負の閾値の間の領域において、前記所定の画像フレームの重みを、ノイズを抑制するような値で一定に保つ請求項1〜4のいずれか一項に記載の超音波診断装置。
  6. 前記正の閾値と前記負の閾値は、絶対値が同じ値となるように定められ、
    前記フレーム相関処理部は、前記正の閾値より大きな領域と前記負の閾値より小さな領域において前記所定の画像フレームの重みを直線的に変化させ、且つ、その変化の傾きを前記負の閾値より小さな領域と比べて前記正の閾値より大きな領域において緩やかにする請求項5に記載の超音波診断装置。
  7. 前記過去の画像フレームは、前記所定の画像フレームとフレーム相関処理する前段階において、前記過去の画像フレームより以前の画像フレームとの間でフレーム相関処理が施された画像フレームである請求項1〜6のいずれか一項に記載の超音波診断装置。
  8. 前記過去の画像フレームは、前記所定の画像フレームとフレーム相関処理する前段階において、前記過去の画像フレームより以前の画像フレームとの間でフレーム相関処理が施されていない画像フレームである請求項1〜6のいずれか一項に記載の超音波診断装置。
  9. 前記過去の画像フレームは、前記所定の画像フレームとフレーム相関処理する前段階において、前記過去の画像フレームより以前の画像フレームとの間でフレーム相関処理が施された画像フレームと、前記過去の画像フレームより以前の画像フレームとの間でフレーム相関処理が施されていない画像フレームとを含む請求項1〜6のいずれか一項に記載の超音波診断装置。
  10. 超音波探触子から被検体内の対象に向けて超音波ビームの送受信を行い、前記超音波探触子から出力される受信信号に基づいて超音波画像を順次生成すると共に生成された所定の画像フレームに対して過去の画像フレームに基づく画像処理を施す超音波診断装置であって、
    前記受信信号に基づいて超音波画像をフレーム毎に順次生成する画像生成部と、
    前記画像生成部で生成された所定の画像フレームの輝度値から過去の画像フレームの輝度値を差し引いて画素毎に輝度値の差分を算出し、前記輝度値の差分に応じて前記過去の画像フレームに対する前記所定の画像フレームの重みをそれぞれ設定して、前記所定の画像フレームと前記過去の画像フレームを重み付け加算するフレーム相関処理部と
    を備え、
    前記フレーム相関処理部は、前記輝度値の差分が予め定められた負の閾値より小さくなるほど、または、前記輝度値の差分が予め定められた正の閾値より大きくなるほど前記所定の画像フレームの重みを大きく変化させ、且つ、前記所定の画像フレームの重みの原点からの変化量を前記負の閾値より小さな領域と比べて前記正の閾値より大きな領域において大きくすることを特徴とする超音波診断装置。
JP2013071825A 2013-03-29 2013-03-29 超音波診断装置 Active JP5956375B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2013071825A JP5956375B2 (ja) 2013-03-29 2013-03-29 超音波診断装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2013071825A JP5956375B2 (ja) 2013-03-29 2013-03-29 超音波診断装置

Publications (2)

Publication Number Publication Date
JP2014195512A true JP2014195512A (ja) 2014-10-16
JP5956375B2 JP5956375B2 (ja) 2016-07-27

Family

ID=52356834

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013071825A Active JP5956375B2 (ja) 2013-03-29 2013-03-29 超音波診断装置

Country Status (1)

Country Link
JP (1) JP5956375B2 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6171246B1 (ja) * 2016-02-09 2017-08-02 本多電子株式会社 超音波画像表示装置及びそのためのプログラムを格納した記録媒体
WO2021020038A1 (ja) 2019-07-26 2021-02-04 富士フイルム株式会社 超音波診断装置および超音波診断装置の制御方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10286256A (ja) * 1997-04-15 1998-10-27 Fujitsu Ltd 超音波診断装置およびフィルタ
JP2002034984A (ja) * 2000-07-21 2002-02-05 Ge Medical Systems Global Technology Co Llc 流れ情報データ平均処理方法および超音波診断装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10286256A (ja) * 1997-04-15 1998-10-27 Fujitsu Ltd 超音波診断装置およびフィルタ
JP2002034984A (ja) * 2000-07-21 2002-02-05 Ge Medical Systems Global Technology Co Llc 流れ情報データ平均処理方法および超音波診断装置

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6171246B1 (ja) * 2016-02-09 2017-08-02 本多電子株式会社 超音波画像表示装置及びそのためのプログラムを格納した記録媒体
WO2017138086A1 (ja) * 2016-02-09 2017-08-17 本多電子株式会社 超音波画像表示装置及び方法、並びにプログラムを格納した記録媒体
WO2021020038A1 (ja) 2019-07-26 2021-02-04 富士フイルム株式会社 超音波診断装置および超音波診断装置の制御方法
US20220125414A1 (en) * 2019-07-26 2022-04-28 Fujifilm Corporation Ultrasound diagnostic apparatus and control method of ultrasound diagnostic apparatus
US11992370B2 (en) 2019-07-26 2024-05-28 Fujifilm Corporation Ultrasound diagnostic apparatus and control method of ultrasound diagnostic apparatus in which an ultrasound probe and a diagnostic apparatus main body are wirelessly connected

Also Published As

Publication number Publication date
JP5956375B2 (ja) 2016-07-27

Similar Documents

Publication Publication Date Title
US20150196278A1 (en) Ultrasound diagnostic apparatus and ultrasound image producing method
US10980515B2 (en) Acoustic wave processing apparatus, signal processing method, and program for acoustic wave processing apparatus
JP6997508B2 (ja) 超音波診断装置及びその制御プログラム
JP2012105950A (ja) 超音波診断装置及びプログラム
JP5663552B2 (ja) 超音波検査装置、超音波検査装置の信号処理方法およびプログラム
JP5956375B2 (ja) 超音波診断装置
JP5325502B2 (ja) 超音波画像形成装置および超音波画像形成方法
JP5777650B2 (ja) 超音波診断装置および超音波画像生成方法
JP5358636B2 (ja) 超音波診断装置および超音波画像生成方法
JP5836197B2 (ja) 超音波診断装置およびデータ処理方法
EP3542346B1 (en) Adaptive ringdown subtraction for coronary and peripheral intravascular ultrasound (ivus)
JP5836241B2 (ja) 超音波検査装置、超音波検査装置の信号処理方法およびプログラム
JP4530834B2 (ja) 超音波画像処理方法及び超音波画像処理装置、並びに、超音波画像処理プログラム
US20160139252A1 (en) Ultrasound diagnostic device, method for generating acoustic ray signal of ultrasound diagnostic device, and program for generating acoustic ray signal of ultrasound diagnostic device
JP2022164443A (ja) 超音波診断装置及び医用画像処理装置
JP5851345B2 (ja) 超音波診断装置およびデータ処理方法
JP5837455B2 (ja) 超音波画像診断装置
JP2020036774A (ja) 超音波診断装置、医用画像処理装置及び医用画像処理プログラム
JP5450488B2 (ja) 超音波診断装置および超音波画像生成方法
JP4664209B2 (ja) 超音波診断装置およびその撮像を実行する超音波イメージングプログラム
JP5917388B2 (ja) 超音波検査装置、超音波検査装置の信号処理方法およびプログラム
US10932759B2 (en) Ultrasound diagnostic apparatus and data processing method
JP2022161534A (ja) 医用画像生成装置、医用画像生成方法、およびプログラム
JP2023104734A (ja) 超音波診断装置、及び画像処理装置
JP4291833B2 (ja) 超音波診断装置および超音波診断画像生成方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20150402

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20160224

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20160329

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160517

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20160607

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160616

R150 Certificate of patent or registration of utility model

Ref document number: 5956375

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250