(57)【要約】      （修正有）【課題】プローブが移動する場合にも正確に変位検出を行うことができる超音波診断装置を提供する。【解決手段】Y方向の動きを検出するためには、X方向の８素子分における時刻の異なる２枚のY方向超音波画像２２が用いられる。Y方向に動くことでラインがずれているので、検出された動きベクトルを用いて、動き補償を伴う変位検出が行われる。これらの処理は、Y方向の３素子分行われる。したがって、この３素子分の結果を集約することで、X方向超音波画像に対して、１スキャンライン分の歪量データが作成される。例えば、超音波画像を撮影するプローブからの信号から、超音波画像を生成し、表示する超音波画像処理装置に適用することができる。【選択図】図３
2015100613.tif 000002
【特許請求の範囲】
【請求項１】
プローブが移動する場合に、第１のフレームと、前記第１のフレームとは異なる時刻のフレームである第２のフレームとの間のずれ量を補正するずれ量補正部と、
前記第１のフレームと前記ずれ量補正部によりずれ量が補正された第２フレームとから検出された位相差に基づいて歪量を算出する歪量算出部と
を備える画像処理装置。
【請求項２】
前記歪量算出部により算出された歪量を、Bモード画像に重畳する重畳部を
さらに備える請求項１に記載の画像処理装置。
【請求項３】
前記第１のフレームおよび前記第２のフレームは、前記Bモード画像が生成される方向に平行な方向以外の方向の画像である
請求項２に記載の画像処理装置。
【請求項４】
前記第１のフレームおよび前記第２のフレームは、前記Bモード画像が生成される方向に平行な方向の画像である
請求項２に記載の画像処理装置。
【請求項５】
前記第１のフレームおよび前記第２のフレームは、前記プローブが移動する方向に平行な方向の画像である
請求項２に記載の画像処理装置。
【請求項６】
前記重畳部は、前記歪量算出部により算出された歪量を、前記歪量が算出された方向に直交する方向のBモード画像に重畳する
請求項２の記載の画像処理装置。
【請求項７】
前記重畳部は、前記歪量算出部により算出された歪量を、前記歪量が算出された方向に平行する方向のBモード画像に重畳する
請求項２の記載の画像処理装置。
【請求項８】
前記歪量算出部により算出された歪量を、被写体に関する内視鏡画像に重畳する重畳部を
さらに備える請求項１に記載の画像処理装置。
【請求項９】
画像処理装置が、
プローブが移動する場合に、第１のフレームと、前記第１のフレームとは異なる時刻のフレームである第２のフレームとの間のずれ量を補正し、
前記第１のフレームと前記ずれ量が補正された第２フレームとから検出された位相差に基づいて歪量を算出する
画像処理方法。
【請求項１０】
プローブと、
画像処理部と
を備え、
前記画像処理部は、プローブが移動する場合に、第１のフレームと、前記第１のフレームとは異なる時刻のフレームである第２のフレームとの間のずれ量を補正するずれ量補正部と、
前記第１のフレームと前記ずれ量補正部によりずれ量が補正された第２フレームとから検出された位相差に基づいて歪量を算出する歪量算出部と
を備える超音波診断装置。
【請求項１１】
前記画像処理部は、前記歪量算出部により算出された歪量を、Bモード画像に重畳する重畳部を
さらに備える請求項１０に記載の超音波診断装置。
【請求項１２】
前記プローブは、第１の方向および第１の方向に直交する第２の方向に、それぞれ複数の圧電素子を有し、
前記第１のフレームおよび前記第２のフレームは、前記Bモード画像が生成される方向に平行な方向以外の方向の画像である
請求項１１に記載の超音波診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、画像処理装置および方法、並びに超音波診断装置に関し、特に、プローブが移動する場合に正確に変位検出を行うことができるようにした画像処理装置および方法、並びに超音波診断装置に関する。
【背景技術】
【０００２】
超音波を用いて組織の硬さを評価するエラストグラフィにおいては、プローブを静止した状態で、深さ方向に振動させることで組織の変位を検出している。
【０００３】
しかしながら、プローブを動かした状態では、組織の変位を正しく検出できない。これは、変位検出で用いる対応点が動きにより移動していまうが、その補正をしていないためである。
【０００４】
なお、音響放射圧インパルス(ARFI)においては、所望の照射位置と実際の照射位置のずれ量の検出を行うものがある（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１０−１７２６９９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
すなわち、プローブを動かした状態では、変位検出で用いる対応点が動きにより移動していまい、フレーム間で対応するラインが揃っていない。このため、プローブを動かした状態では、組織の変位を正しく検出できない。
【０００７】
本開示は、このような状況に鑑みてなされたものであり、プローブが移動する場合に正確に変位検出を行うことができるものである。
【課題を解決するための手段】
【０００８】
本開示の一側面の画像処理装置は、プローブが移動する場合に、第１のフレームと、前記第１のフレームとは異なる時刻のフレームである第２のフレームとの間のずれ量を補正するずれ量補正部と、前記第１のフレームと前記ずれ量補正部によりずれ量が補正された第２フレームとから検出された位相差に基づいて歪量を算出する歪量算出部とを備える。
【０００９】
前記歪量算出部により算出された歪量を、Bモード画像に重畳する重畳部をさらに備えることができる。
【００１０】
前記第１のフレームおよび前記第２のフレームは、前記Bモード画像が生成される方向に平行な方向以外の方向の画像である。
【００１１】
前記第１のフレームおよび前記第２のフレームは、前記Bモード画像が生成される方向に平行な方向の画像である。
【００１２】
前記第１のフレームおよび前記第２のフレームは、前記プローブが移動する方向に平行な方向の画像である。
【００１３】
前記重畳部は、前記歪量算出部により算出された歪量を、前記歪量が算出された方向に直交する方向のBモード画像に重畳することができる。
【００１４】
前記重畳部は、前記歪量算出部により算出された歪量を、前記歪量が算出された方向に平行する方向のBモード画像に重畳することができる。
【００１５】
前記歪量算出部により算出された歪量を、被写体に関する内視鏡画像に重畳する重畳部をさらに備えることができる。
【００１６】
本開示の一側面の画像処理方法は、画像処理装置が、プローブが移動する場合に、第１のフレームと、前記第１のフレームとは異なる時刻のフレームである第２のフレームとの間のずれ量を補正し、前記第１のフレームと前記ずれ量が補正された第２フレームとから検出された位相差に基づいて歪量を算出する。
【００１７】
本開示の一側面の超音波診断装置は、プローブと、画像処理部とを備え、前記画像処理部は、プローブが移動する場合に、第１のフレームと、前記第１のフレームとは異なる時刻のフレームである第２のフレームとの間のずれ量を補正するずれ量補正部と、前記第１のフレームと前記ずれ量補正部によりずれ量が補正された第２フレームとから検出された位相差に基づいて歪量を算出する歪量算出部とを備える。
【００１８】
前記画像処理部は、前記歪量算出部により算出された歪量を、Bモード画像に重畳する重畳部をさらに備えることができる。
【００１９】
前記プローブは、第１の方向および第１の方向に直交する第２の方向に、それぞれ複数の圧電素子を有し、前記第１のフレームおよび前記第２のフレームは、前記Bモード画像が生成される方向に平行な方向以外の方向の画像である。
【００２０】
本開示の一側面においては、プローブが移動する場合に、第１のフレームと、前記第１のフレームとは異なる時刻のフレームである第２のフレームとの間のずれ量が補正される。そして、前記第１のフレームと前記ずれ量が補正された第２フレームとから検出された位相差に基づいて歪量が算出される。
【発明の効果】
【００２１】
本開示によれば、変位検出を行うことができる。特に、プローブが移動する場合に正確に変位検出を行うことができる。
【図面の簡単な説明】
【００２２】
【図１】プローブの軸と動きについて説明する図である。
【図２】本技術において用いられるプローブの構成例を示す図である。
【図３】本技術の概要を説明する図である。
【図４】プローブの偏心について説明する図である。
【図５】本技術を適用した超音波画像処理装置の構成例を示すブロック図である。
【図６】Bモード画像生成部の構成例を示すブロック図である。
【図７】歪画像生成部の構成例を示すブロック図である。
【図８】フレーム間変位算出部の構成例を示すブロック図である。
【図９】IQデータの補間処理の例を示す図である。
【図１０】超音波画像処理の例について説明するフローチャートである。
【図１１】コンピュータの構成例を示すブロック図である。
【発明を実施するための形態】
【００２３】
以下、本開示を実施するための形態（以下実施の形態とする）について説明する。なお、説明は以下の順序で行う。
１．第１の実施の形態（超音波画像処理装置）
２．第２の実施の形態（コンピュータ）
【００２４】
＜第１の実施の形態＞
＜プローブの軸と動きの例＞
まず、図１を参照して、プローブの軸と動きについて説明する。なお、説明を簡単にするため、1Dのプローブが直方体の形状をしていることとする。
【００２５】
プローブ１１の軸として、長軸方向、短軸方向、深さ方向へ、X軸、Y軸、Z軸をそれぞれ割り当て、以下、XYZを用いて説明する。プローブ１１は、複数の圧電素子が配列されて構成される。長軸方向とは、配列する圧電素子の数が短軸方向より多く、圧電素子の全長が長い方向である。短軸方向とは、配列する圧電素子の数が長軸方向より少なく、圧電素子の全長が長軸方向よりも短い方向である。深さ方向とは、超音波出力方向である深さ方向である。
【００２６】
例えば、1Dプローブである場合に、X方向の配列されている圧電素子で受信したRF信号を並べることで超音波画像１２が生成される。
【００２７】
なお、プローブ１１の動きとして、Z方向に対しては、プローブ１１からのエラストグラフィを取得するための振動がある。また、X方向とY方向に対しては、使用者の操作などによる規則的または不規則的な運動が想定される。一例として、Y方向へはスイープ動作による等速移動、X方向へは手振れなどによる不規則な移動などが考えられる。
【００２８】
＜プローブの構成例＞
図２は、本技術において用いられるプローブの構成例を示す図である。なお、本技術においては、図１で示したY方向の動きを検知するために、Y方向に２つ以上の圧電素子（例えば、振動子）３１が配列されている必要がある。すなわち、プローブ２１は、1.5Dまたは2Dアレイのプローブとして構成される。なお、後述するが、使用する画像と重畳する画像の方向によっては、1Dのプローブも用いることができる。
【００２９】
図２の例においては、圧電素子３１が、X方向へ8素子、Y方向へ3素子配列された1.5Dのプローブ２１が示されている。
【００３０】
図１において上述したように、各圧電素子３１で受信したRF信号を並べることで超音波画像が生成される。図２の例においては、X方向に着目すると、８つの圧電素子３１が３列に配列されているので、３枚の超音波画像(以下、X方向超音波画像と称する)２３が生成される。同様に、Y方向に着目すると、３つの圧電素子３１が８列に配列されているので、８枚の超音波画像(以下、Y方向超音波画像と称する)２２が生成される。X方向超音波画像２３は、Bモード画像生成に用いられる。
【００３１】
このように、本技術のプローブにおいては、Y方向にも超音波画像が生成されるので、それを用いてブロックマッチングができて、変位検出で用いる対応点を追跡することができる。
【００３２】
なお、図２の例においては図示されないが、2Dアレイのプローブの場合、圧電素子は、X方向とY方向に同数配列される。
【００３３】
＜本技術の概要＞
図３は、本技術を具体的に説明する図である。図３Ａに示されるように、Y方向の動きがある場合が示されている。Y方向の動きを検出するためには、図３Ｂに示されるように、X方向の８素子分における時刻の異なる２枚（時刻tと時刻t-1の２枚）のY方向超音波画像２２が用いられる。この動き検出には、例えば、ブロックマッチングを使った手法などを適用することができる。これにより、動きベクトルを求めることができる。
【００３４】
図３の例においては、Y方向に３素子配列されているので、素子間ピッチ×２[um]の動きまで検出することが可能である。Y方向により多くの圧電素子３１が並ぶほど、大きな動きの検出が可能となる。
【００３５】
次に、Y方向に動くことで、フレーム間で対応するラインがずれているので、図３Ｂに示されるように、検出された動きベクトルを用いて、動き補償を伴う変位検出が行われる。すなわち、検出された動きベクトルを用いて、t-1の画像に対して動き補償が行われ、動き補償画像が生成される。なお、この際、動きベクトルが素子間ピッチよりも短い場合、補間処理が行われて、小数位置の画素が作られる。よって、素子間ピッチよりも小さな動きも補償可能である。
【００３６】
そして、時刻tの画像と時刻t-1の動き補償画像の間で、位相差検出処理を行う。例えば、位相差検出処理として、自己相関法を用いることが可能である。そして、位相差検出結果より、変位および歪量が検出される。
【００３７】
なお、上述したこれらの処理は、Y方向の３素子分行われる。したがって、この３素子分の結果を集約することで、X方向超音波画像２３に対して、１スキャンライン分の歪量データが作成される。この集約の方法は、平均やメディアンなどがあるが、集約する方法であればどんな方法でもかまわない。また、集約でなくても、３素子分のうちのどれか１素子分（例えば、真ん中の１素子）を選択して用いる方法でもよい。
【００３８】
さらに、動き検出補償においては、多くの素子分を用い、位相差検出処理においてそのうちの数素子分や１素子分を用いることも可能である。
【００３９】
以上の処理を、X方向超音波画像２３の各スキャンラインに対応する３素子を使って繰り返すことで、すべてのスキャンラインの歪量データを生成することが可能である。この歪量データは、図３Ｃに示されるように、X方向超音波画像２３に重畳して表示される。
【００４０】
以上のように、動きの方向に平行な超音波画像を用いて動き検出補償をして、その結果が、動きの方向と直交する超音波画像(Bモード画像）に重畳して表示される。換言するに、Bモード画像が生成される方向に直交する超音波画像を用いて動き検出補償をして、その結果が、Bモード画像に重畳して表示される。
【００４１】
なお、2Dプローブの場合、例えば、X方向に8素子並び、Y方向に8素子並ぶので、Y方向に動きがあるとき、Y方向（平行方向）の画像で動きを検出補償し、Y方向（平行方向）の画像に重畳することもできる。この場合は、8素子のデータを集約することなしに、Y方向の8面分の画像に対してそれぞれ行うことができる。
【００４２】
また、上記説明においてはY方向の動きを検出する例を説明したが、X方向の動きを検出することも可能である。例えば、1Dプローブの場合、X方向に動きがあるとき、X方向（平行方向）の動き検出補償し、X方向（平行方向）の画像に重畳することができる。
【００４３】
さらに、2Dプローブの場合、X方向に動いていて、X方向（平行方向）に動き検出補償して、Y方向（直交方向）に重畳することもできる。
【００４４】
すなわち、本技術においては、プローブが移動する場合に、異なる時刻の２フレームの間のずれ量を補正し、第１のフレームとずれ量が補正された第２フレームとの位相差を検出して、検出された位相差に基づいて、歪量を算出するようにした。
【００４５】
これにより、プローブが動いた状態で、正確な変位検出を行うことができる。したがって、組織の硬さを正しく求めることができる。
【００４６】
なお、本技術は、X軸、Y軸、Z軸に沿ってだけ動くものではなく、斜めに動く場合にも適用することができる。斜めに動く場合には、動きの成分をX方向とY方向に分解して、各方向で動き検出と変位検出を行い、統計処理（例えば、平均化）をした結果から歪量を求めることができる。
【００４７】
さらに、本技術は、プローブ２１が偏心して動く（回転角度を有する動きを行う）場合にも適用可能である。
【００４８】
＜プローブの偏心＞
図４は、プローブの偏心について説明する図である。
【００４９】
プローブ２１の偏心には、プローブ２１から超音波が出力される方向（超音波出力方向）の、ヨー(yaw)方向の回転角度、ピッチ(pitch)方向の回転角度、及び、ロール(roll)方向の回転角度の３種類がある。
【００５０】
プローブ２１が、例えば、直前に、検査対象者に当てられたときの、そのプローブ２１の位置を原点とし、かつ、z軸を重力の方向とする３次元座標系を、プローブ座標系ということとする。このとき、ピッチ(pitch)方向の回転角度は、プローブ座標系のx軸回りの回転角度を表し、ロール(roll)方向の回転角度は、プローブ座標系のy軸回りの回転角度を表し、ヨー(yaw)方向の回転角度は、プローブ座標系のz軸回りの回転角度を表す。
【００５１】
ここで、説明を簡単にするため、プローブ２１が直方体の形状をしていることとする。超音波出力方向がプローブ座標系のz軸と一致している状態のプローブ２１としての直方体の面のうちの、超音波出力方向と平行な１つの面としての正面２１Ａが、例えば、図４に示すように、手前側を向いている場合と、z軸回りに180度回転して向こう側を向いている場合とでは、プローブ２１の見た目の状態は同一であるが、ヨー(yaw)方向の回転角度が異なる。
【００５２】
また、図４の例においては、超音波画像の座標系のx軸及びy軸を図示してあるが、超音波画像の座標系のx軸及びy軸は、それぞれ、超音波画像の横方向及び縦方向を表し、プローブ座標系のx軸及びy軸とは異なる。
【００５３】
なお、ピッチ(pitch)方向の回転は、x方向を２枚使って動き検出可能であり、ロール(roll)方向の回転は、y方向を２枚使って動き検出可能である。ヨー(yaw)方向の回転は、斜め方向の動きと同様に、x方向とy方向に分解して動き検出することで対応可能である。
【００５４】
このようにプローブ２１が偏心して動く（回転角度を有する動きを行う）場合にも、本技術を適用することができる。例えば、アフィンパラメータなどをフレーム間で求めることで、斜めに動く場合などと同様に対応することができる。
【００５５】
＜超音波画像処理装置の構成例＞
図５は、本技術の画像処理装置を適用した超音波画像処理装置の構成例を示す図である。
【００５６】
図５に示される超音波画像処理装置１００は、超音波を用いて検査対象者（対象物）の内部の画像（すなわち、超音波画像）を撮影し、表示する装置である。超音波画像処理装置１００は、例えば、医療用として患者の身体の内部や胎児などの撮影に用いられたり、手術中に臓器や組織に直接当てて用いられたり、工業用として製品の内部の断面などの撮影に用いられたりする。
【００５７】
特に、この超音波画像処理装置１００は、被写体の超音波画像を撮影し、超音波を用いて組織の硬さを評価（検査）するエラストグラフィ画像（歪画像とも称する）を表示する。
【００５８】
超音波画像処理装置１００は、図２のプローブ２１および画像処理部１１１を含むように構成されている。
【００５９】
プローブ２１は、被写体（生体；例えば、皮膚）などに押しつけて、被写体の画像を撮影する部位である。なお、プローブ２１は、1.5Dアレイであってもよいし、2Dアレイであってもよい。
【００６０】
プローブ２１は、超音波を出力するとともに、検査対象者（の体内）で超音波が反射されることにより戻ってくる反射波を受信する超音波プローブであり、超音波送受信部１２１を有する。
【００６１】
超音波送受信部１２１は、例えば、プローブ２１の先端に設けられている。超音波送受信部１２１は、超音波を発生して出力するとともに、検査対象者で反射された超音波の反射波を受信し、その反射波の強度を表す反射波データ(RFデータ)を、画像処理部１１１の複素化部１３１に出力する。
【００６２】
画像処理部１１１は、複素化部１３１、X/Y方向成分分離部１３２、Bモード画像生成部１３３、歪画像生成部１３４、および画像重畳部１３５を含むように構成されている。
【００６３】
複素化部１３１には、プローブ２１において受信されたN×M本のRFデータが入力される。ここで、Nは、X方向の素子数、MはY方向の素子数である。複素化部１３１は、実数のRFデータを複素数のIQデータに変換する。例えば、直交変換やヒルベルト変換が用いられる。複素化部１３１は、IQデータをX/Y方向成分分離部１３２に出力する。
【００６４】
X/Y方向成分分離部１３２は、複素化部１３１からのIQデータのペア（X方向IQデータとY方向IQデータ）を生成する。X方向IQデータは、X方向に並んだIQデータであり、例えば、M/2の位置のN本のIQデータ群である。Y方向IQデータは、Y方向に並んだIQデータであり、例えば、0乃至N-1の各位置におけるM本のIQデータ群である。
【００６５】
X方向IQデータは、Bモード画像生成に用いられる。Y方向IQデータは、動きを考慮した歪量を算出するために用いられる。X/Y方向成分分離部１３２は、X方向IQデータをBモード画像生成部１３３に出力し、Y方向IQデータを歪画像生成部１３４に出力する。
【００６６】
Bモード画像生成部１３３は、X/Y方向成分分離部１３２からのX方向IQデータから、X方向B(Brightness)モード画像（輝度画像）を生成し、生成したX方向Bモード画像を、画像重畳部１３５に出力する。
【００６７】
歪画像生成部１３４は、X/Y方向成分分離部１３２からのY方向IQデータを用いて、Z方向振動に伴う変位を検出し、変位量から歪量を算出し、その歪量を色情報に変換したY方向歪画像（エラストグラフィ画像）を、画像重畳部１３５に出力する。
【００６８】
画像重畳部１３５は、Bモード画像生成部１３３からのX方向のBモード画像と、歪画像生成部１３４からのY方向のエラストグラフィ画像を重畳する。これにより、画像重畳部１３５からは、Bモード画像に重畳されたエラストグラフィ画像（歪画像）が出力され、図示せぬ後段の表示部などに表示される。
【００６９】
＜Bモード画像生成部の構成＞
図６は、Bモード画像生成部の構成例を示すブロック図である。
【００７０】
図６の例において、Bモード画像生成部１３３は、包絡線算出部１４１および対数変換部１４２を含むように構成されている。
【００７１】
包絡線算出部１４１は、X方向IQデータを用いて、包絡線（絶対値）を算出する。対数変換部１４２は、包絡線算出部１４１により算出された包絡線に基づいて、対数変換を行い、包絡線の振幅を輝度に変換することで、X方向Bモード画像（輝度画像）を生成する。対数変換部１４２は、生成したX方向Bモード画像を、画像重畳部１３５に供給する。
【００７２】
＜歪画像生成部の構成＞
図７は、歪画像生成部の構成例を示すブロック図である。なお、歪画像生成部１３４における処理は、0乃至N-1の各位置におけるM本のIQデータ群のそれぞれに対して実行される。
【００７３】
図７の例において、歪画像生成部１３４は、フレームメモリ１５１、フレーム間変位算出部１５２、歪量算出部１５３、および可視化処理部１５４を含むように構成されている。
【００７４】
フレームメモリ１５１およびフレーム間変位算出部１５２には、図５のX/Y方向成分分離部１３２からのY方向IQデータが入力される。
【００７５】
フレームメモリ１５１は、X/Y方向成分分離部１３２からのY方向IQデータ（異なる時刻（例えば、直前）のデータ）を蓄積する。
【００７６】
フレーム間変位算出部１５２は、X/Y方向成分分離部１３２からの時刻tのY方向IQデータ（以下、IQデータ(t)とも称する）と、フレームメモリ１５１からの時刻t-1のY方向IQデータ（以下、IQデータ(t-1)とも称する）とを用いて、Z方向振動に伴う変位を検出する。その際、Y方向に動くことにより発生するラインのずれが検出され、補正されて、Z方向振動に伴う変位が検出される。なお、時刻t-1のデータは、直前のデータを意味するが、これに限らず、異なる時刻のデータであればよい。フレーム間変位算出部１５２は、検出した変位の量を、歪量算出部１５３に供給する。
【００７７】
歪量算出部１５３は、フレーム間変位算出部１５２からの変位量に基づいて歪量を算出する。例えば、歪量は、変位量を微分することで求められる。例えば、微分値の傾きが大きいと大きな歪みがあることを示し、傾きが小さいと歪みが小さいことを示す。歪量算出部１５３は、算出した歪量を可視化処理部１５４に供給する。
【００７８】
可視化処理部１５４は、歪量をカラー情報に変換する。なお、歪みが大きいと、柔らかいということであり、歪みが小さいと硬いということである。例えば、歪量の画面内平均値をGreen、歪量の画面内平均よりも小さな値をBlue、歪量の画面内平均よりも大きな値をRedとして可視化される。可視化処理部１５４は、変換したカラー情報を、エラストグラフィ画像（Y方向歪画像）として、図５の画像重畳部１３５に出力する。
【００７９】
＜フレーム間変位算出部の構成＞
図８は、フレーム間変位算出部の構成例を示すブロック図である。
【００８０】
図８の例において、フレーム間変位算出部１５２は、絶対値／位相分離部１６１−１および１６１−２、動き検出部１６２、動き補償部１６３、および位相差検出部１６４を含むように構成される。
【００８１】
図６のX/Y方向成分分離部１３２からのY方向IQデータ(t)は、絶対値／位相分離部１６１−１および位相差検出部１６４に入力される。図７のフレームメモリ１５１からのY方向IQデータ(t-1)は、絶対値／位相分離部１６１−２および動き補償部１６３に入力される。
【００８２】
絶対値／位相分離部１６１−１は、Y方向IQデータ(t)から絶対値と位相とを分離し、IQデータ(t)の絶対値を取得し、取得した絶対値を、動き検出部１６２に出力する。絶対値／位相分離部１６１−２は、Y方向IQデータ(t-1)から絶対値と位相とを分離し、IQデータ(t-1)の絶対値を取得し、取得した絶対値を、動き検出部１６２に出力する。
【００８３】
動き検出部１６２は、IQデータの絶対値成分を用いて、時刻tと時刻t-1の間の動き量を検出する。例えば、動き量は、ブロックマッチングを用いて、動きベクトルとして検出される。なお、動き検出により位置合わせが行われる。すなわち、位置合わせができれば、動き検出に限定されず、また、動き検出方法もブロックマッチングに限定されない。動き検出部１６２は、検出された動きベクトルを、動き補償部１６３に出力する。
【００８４】
動き補償部１６３は、動き検出部１６２により検出された動きベクトルを用いて、時刻t-1のIQデータに対して補間処理を行い、動き補償後の複素データを生成する。この補間処理について図９を参照して説明する。
【００８５】
図９は、補間処理の例を示す図である。図９の例において、mとm+1は、隣接するIQデータである。移動量wが素子間ピッチ未満のとき、mとm+1の間のIQデータを補間処理によって作ることになる。例えば、移動量wを重み量とした線形補間を、複素数データに適用することで、補間データ（動き補償後の複素データ）を作ることができる。動き補償部１６３は、動き補償後の複素データを、位相差検出部１６４に出力する。
【００８６】
位相差検出部１６４は、時刻tのIQデータと、動き補償後の時刻t-1のIQデータの間の位相差を検出し、変位量を算出する。例えば、位相差検出の方法としては、自己相関法などが用いられる。位相差検出部１６４は、算出した変位量を、歪量算出部１５３に供給する。
【００８７】
＜超音波画像生成処理＞
次に、図１０のフローチャートを参照して、超音波画像処理装置１００の超音波画像処理について説明する。
【００８８】
プローブ２１の超音波送受信部１２１は、超音波を発生して出力するとともに、検査対象者で反射された超音波の反射波を受信し、その反射波の強度を表す反射波データ(RFデータ)を、画像処理部１１１の複素化部１３１に供給する。
【００８９】
ステップＳ１１１において、複素化部１３１は、２次元RFデータを取得する。
【００９０】
ステップＳ１１２において、複素化部１３１は、取得したRFデータの複素化を行う。すなわち、複素化部１３１は、実数のRFデータを複素数のIQデータに変換し、変換されたIQデータを、X/Y方向成分分離部１３２に出力する。
【００９１】
ステップＳ１１３において、X/Y方向成分分離部１３２は、X方向複素数データとY方向複素数データとに分離を行う。X/Y方向成分分離部１３２は、X方向IQデータをBモード画像生成部１３３に出力し、Y方向IQデータを歪画像生成部１３４に出力する。
【００９２】
ステップＳ１１４において、Bモード画像生成部１３３は、X/Y方向成分分離部１３２からのX方向IQデータから、X方向Bモード画像（輝度画像）を生成する。すなわち、包絡線算出部１４１は、X方向IQデータを用いて、包絡線（絶対値）を算出する。対数変換部１４２は、包絡線算出部１４１により算出された包絡線に基づいて、対数変換を行い、包絡線の振幅を輝度に変換することで、X方向Bモード画像（輝度画像）を生成する。
Bモード画像生成部１３３（対数変換部１４２）は、生成したX方向Bモード画像を、画像重畳部１３５に出力する。
【００９３】
歪画像生成部１３４において、図６のX/Y方向成分分離部１３２からのY方向IQデータ(t)は、絶対値／位相分離部１６１−１および位相差検出部１６４に入力される。図７のフレームメモリ１５１からのY方向IQデータ(t-1)は、絶対値／位相分離部１６１−２および動き補償部１６３に入力される。
【００９４】
ステップＳ１１５において、絶対値／位相分離部１６１−１および１６１−２は、Y方向IQデータの絶対値を取得する。すなわち、絶対値／位相分離部１６１−１は、Y方向IQデータ(t)から絶対値と位相とを分離し、IQデータ(t)の絶対値を取得し、取得した絶対値を、動き検出部１６２に出力する。絶対値／位相分離部１６１−２は、Y方向IQデータ(t-1)から絶対値と位相とを分離し、IQデータ(t-1)の絶対値を取得し、取得した絶対値を、動き検出部１６２に出力する。
【００９５】
ステップＳ１１６において、動き検出部１６２は、時刻tと時刻t-nのフレーム間で動き検出を行い、動きベクトルを求める。動き検出部１６２は、検出された動きベクトルを、動き補償部１６３に出力する。
【００９６】
ステップＳ１１７において、動き補償部１６３は、動き検出部１６２により検出された動きベクトルを用いて、時刻t-nのフレームの動き補償を行う。すなわち、動き補償部１６３は、IQデータ(t-n)に対して補間処理を行い、動き補償後の複素データを生成する。動き補償部１６３は、動き補償後の複素データを、位相差検出部１６４に出力する。
【００９７】
ステップＳ１１８において、位相差検出部１６４は、動き補償フレーム（動き補償後の時刻t-1のIQデータ）と、時刻tのIQデータとの間の位相差を検出し、変位量を算出する。位相差検出部１６４は、算出した変位量を、歪量算出部１５３に出力する。
【００９８】
ステップＳ１１９において、歪量算出部１５３は、フレーム間変位算出部１５２からの変位量を歪量に変換する。例えば、歪量は、変位量を微分することで求められる。歪量算出部１５３は、変換により算出した歪量を可視化処理部１５４に供給する。
【００９９】
次に、可視化処理部１５４は、歪量をカラー情報に変換する。可視化処理部１５４は、変換したカラー情報を、エラストグラフィ画像（Y方向歪画像）として、図５の画像重畳部１３５に出力する。
【０１００】
ステップＳ１２０において、画像重畳部１３５は、Bモード画像生成部１３３からのX方向のBモード画像と、歪画像生成部１３４からのY方向のエラストグラフィ画像を重畳する。これにより、画像重畳部１３５からは、Bモード画像に重畳されたエラストグラフィ画像（歪画像）が出力される。
【０１０１】
本技術においては、プローブが移動する場合に、プローブが移動する方向のフレーム間のずれ量を補正した後で、位相差を検出して歪み量を算出するようにした。
【０１０２】
これにより、プローブが動いた状態で、正確な変位検出を行うことができる。したがって、組織の硬さを正しく求めることができる。
【０１０３】
また、プローブを動かしながら、組織の硬さを評価できるため、検査時間の短縮になり、術者（検査対象者）の負担軽減にもつながる。
【０１０４】
なお、プローブの動きとしては、図１を参照して上述したように、使用者の故意の操作による動きだけでなく、使用者の意図せぬ手振れなどの動きも含まれる。すなわち、手振れをしたとしても、正確な変位検出を行うことができる。
【０１０５】
さらに、上記説明においては、エラストグラフィ画像（歪画像）に、Bモード画像を合成して表示する例を説明したが、Bモード画像を合成することなしに表示させてもよい。また、合成する画像としては、Bモード画像だけに限らず、どんな画像であってもよい。具体的には、エラストグラフィ画像に、被写体に関する画像として、被写体のCT画像、被写体のMRI画像、被写体の内視鏡画像（内視鏡により撮影されたビデオ画像または静止画像）などを合成することも可能である。
【０１０６】
なお、本技術は、プローブを故意に動かす場合に限らず、手振れなどで動く場合などにも適用することができる。
【０１０７】
また、本技術は、医療用途および非医療用途のいずれにも用いることが可能である。また、本技術は、人間だけでなく、例えば、動物や植物、人工物など、超音波により被写体の断面の撮影を行う様々な場面に用いることができる。
【０１０８】
上述した一連の処理は、ハードウエアにより実行することもできるし、ソフトウエアにより実行することもできる。一連の処理をソフトウエアにより実行する場合には、そのソフトウエアを構成するプログラムが、コンピュータにインストールされる。ここで、コンピュータには、専用のハードウエアに組み込まれているコンピュータや、各種のプログラムをインストールすることで、各種の機能を実行することが可能な汎用のパーソナルコンピュータなどが含まれる。
【０１０９】
＜第２の実施の形態＞
[コンピュータの構成例]
図１１は、上述した一連の処理をプログラムにより実行するコンピュータのハードウエアの構成例を示すブロック図である。
【０１１０】
コンピュータにおいて、CPU(Central Processing Unit)４０１、ROM(Read Only Memory)４０２、RAM(Random Access Memory)４０３は、バス４０４により相互に接続されている。
【０１１１】
バス４０４には、さらに、入出力インタフェース４０５が接続されている。入出力インタフェース４０５には、入力部４０６、出力部４０７、記憶部４０８、通信部４０９、およびドライブ４１０が接続されている。
【０１１２】
入力部４０６は、キーボード、マウス、マイクロホンなどよりなる。出力部４０７は、ディスプレイ、スピーカなどよりなる。記憶部４０８は、ハードディスクや不揮発性のメモリなどよりなる。通信部４０９は、ネットワークインタフェースなどよりなる。ドライブ４１０は、磁気ディスク、光ディスク、光磁気ディスク、又は半導体メモリなどのリムーバブル記録媒体４１１を駆動する。
【０１１３】
以上のように構成されるコンピュータでは、CPU４０１が、例えば、記憶部４０８に記憶されているプログラムを入出力インタフェース４０５及びバス４０４を介してRAM４０３にロードして実行することにより、上述した一連の処理が行われる。
【０１１４】
コンピュータ（CPU４０１）が実行するプログラムは、例えば、パッケージメディア等としてのリムーバブル記録媒体４１１に記録して提供することができる。また、プログラムは、ローカルエリアネットワーク、インターネット、デジタル放送といった、有線または無線の伝送媒体を介して提供することができる。
【０１１５】
コンピュータでは、プログラムは、リムーバブル記録媒体４１１をドライブ４１０に装着することにより、入出力インタフェース４０５を介して、記憶部４０８にインストールすることができる。また、プログラムは、有線または無線の伝送媒体を介して、通信部４０９で受信し、記憶部４０８にインストールすることができる。その他、プログラムは、ROM４０２や記憶部４０８に、あらかじめインストールしておくことができる。
【０１１６】
なお、コンピュータが実行するプログラムは、本明細書で説明する順序に沿って時系列に処理が行われるプログラムであっても良いし、並列に、あるいは呼び出しが行われたとき等の必要なタイミングで処理が行われるプログラムであっても良い。
【０１１７】
また、本明細書において、システムの用語は、複数の装置、ブロック、手段などにより構成される全体的な装置を意味するものである。
【０１１８】
なお、本開示における実施の形態は、上述した実施の形態に限定されるものではなく、本開示の要旨を逸脱しない範囲において種々の変更が可能である。
【０１１９】
以上、添付図面を参照しながら本開示の好適な実施形態について詳細に説明したが、開示はかかる例に限定されない。本開示の属する技術の分野における通常の知識を有するであれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例また修正例に想到し得ることは明らかであり、これらについても、当然に本開示の技術的範囲に属するものと了解される。
【０１２０】
なお、本技術は以下のような構成も取ることができる。
（１）  プローブが移動する場合に、第１のフレームと、前記第１のフレームとは異なる時刻のフレームである第２のフレームとの間のずれ量を補正するずれ量補正部と、
前記第１のフレームと前記ずれ量補正部によりずれ量が補正された第２フレームとから検出された位相差に基づいて歪量を算出する歪量算出部と
を備える画像処理装置。
（２）  前記歪量算出部により算出された歪量を、Bモード画像に重畳する重畳部を
さらに備える前記（１）に記載の画像処理装置。
（３）  前記第１のフレームおよび前記第２のフレームは、前記Bモード画像が生成される方向に平行な方向以外の方向の画像である
前記（２）に記載の画像処理装置。
（４）  前記第１のフレームおよび前記第２のフレームは、前記Bモード画像が生成される方向に平行な方向の画像である
前記（２）に記載の画像処理装置。
（５）  前記第１のフレームおよび前記第２のフレームは、前記プローブが移動する方向に平行な方向の画像である
前記（２）に記載の画像処理装置。
（６）  前記重畳部は、前記歪量算出部により算出された歪量を、前記歪量が算出された方向に直交する方向のBモード画像に重畳する
前記（２）に記載の画像処理装置。
（７）  前記重畳部は、前記歪量算出部により算出された歪量を、前記歪量が算出された方向に平行する方向のBモード画像に重畳する
前記（２）に記載の画像処理装置。
（８）  前記歪量算出部により算出された歪量を、被写体に関する内視鏡画像に重畳する重畳部を
さらに備える前記（２）に記載の画像処理装置。
（９）  画像処理装置が、
プローブが移動する場合に、第１のフレームと、前記第１のフレームとは異なる時刻のフレームである第２のフレームとの間のずれ量を補正し、
前記第１のフレームと前記ずれ量が補正された第２フレームとから検出された位相差に基づいて歪量を算出する
画像処理方法。
（１０）  プローブと、
画像処理部と
を備え、
前記画像処理部は、プローブが移動する場合に、第１のフレームと、前記第１のフレームとは異なる時刻のフレームである第２のフレームとの間のずれ量を補正するずれ量補正部と、
前記第１のフレームと前記ずれ量補正部によりずれ量が補正された第２フレームとから検出された位相差に基づいて歪量を算出する歪量算出部と
を備える超音波診断装置。
（１１）  前記画像処理部は、前記歪量算出部により算出された歪量を、Bモード画像に重畳する重畳部を
さらに備える前記（１０）に記載の超音波診断装置。
（１２）  前記プローブは、第１の方向および第１の方向に直交する第２の方向に、それぞれ複数の圧電素子を有し、
前記第１のフレームおよび前記第２のフレームは、前記Bモード画像が生成される方向に平行な方向以外の方向の画像である
前記（１１）に記載の超音波診断装置。
【符号の説明】
【０１２１】
２１  プローブ，  ３１  圧電素子，  ２２  X方向超音波画像，  ２３  Y方向超音波画像，  １００  超音波画像処理装置，  １１１  画像処理部，  １２１  超音波送受信部，  １３１  複素化部，  １３２  X/Y方向成分分離部，  １３３  Bモード画像生成部，  １３４  歪画像生成部，  １３５  画像重畳部，  １４１  包絡線算出部，  １４２  対数変換部，  １５１  フレームメモリ，  １５２  フレーム間変位算出部，  １５３  歪量算出部，  １５４  可視化処理部，  １６１−１，１６１−２  絶対値／位相分離部，  １６２  動き検出部，  １６３  動き補償部，  １６４  位相差検出部
【図１】
2015100613.tif 000003
【図２】
2015100613.tif 000004
【図３】
2015100613.tif 000005
【図４】
2015100613.tif 000006
【図５】
2015100613.tif 000007
【図６】
2015100613.tif 000008
【図７】
2015100613.tif 000009
【図８】
2015100613.tif 000010
【図９】
2015100613.tif 000011
【図１０】
2015100613.tif 000012
【図１１】
2015100613.tif 000013