JP6689568B2 - 超音波測定装置及び超音波測定方法 - Google Patents

Description

本発明は、超音波測定装置及び超音波測定方法に関し、特に、サンプル内の反射体に対する画像を測定する超音波測定装置及び超音波測定方法に関する。

超音波測定装置は、サンプル内のフォーカシング地点に超音波信号を送信し、サンプル内の反射体によって生成されたアコースティックフィールドを受信し、受信されたアコースティックフィールドを用いてサンプル内の反射体画像を測定（生成）する装置である。
さらに、超音波測定装置は、リアルタイムでサンプルに対する画像を処理することができ、人体内部の組織検査のような医療用又は非破壊検査用などの様々な産業分野において利用されている。

超音波測定装置は、複数のトランスデューサー素子を含み得、複数のトランスデューサー素子の送信時期を調整してサンプル内の特定の地点にビームフォーカシングすることができる。
超音波測定装置は、ビームフォーカシングを行って超音波信号を送信し、反射体によって反射又は散乱されたアコースティックフィールドを用いて反射体に対する画像を測定（生成）する。

しかし、サンプル内の測定しようとする反射体以外の他の反射体がさらに含まれる場合には、他の反射体からのアコースティックフィールドが測定しようとする反射体からのアコースティックフィールドと干渉することがある。
そのため、測定しようとする反射体を正確に画像測定（生成）することが困難であるという問題があった。

本発明は上記従来の超音波測定装置における問題点に鑑みてなされたものであって、本発明の目的は、反射体の正確な画像生成ができる超音波測定方法及び超音波測定装置を提供することにある。

上記目的を達成するためになされた本発明による超音波測定方法は、実際フォーカシング地点に超音波信号をビームフォーカシングするステップと、ビームフォーカシングに対応して、反射体から反射されるアコースティックフィールドを受信してトランスデューサー素子のチャネルごとのデータを含む受信データを取得するステップと、前記受信データに基づいて複数の画像フレームを生成し、前記複数の画像フレームに基づいて前記トランスデューサー素子のチャネルごとのデータに動き補正を行うステップと、前記実際フォーカシング地点と前記反射体の周辺領域の合成フォーカシング地点との間の第１ディレイに基づいて前記合成フォーカシング地点に合成フォーカシング方法を適用し、前記合成フォーカシング地点に対応する少なくとも１つの合成データを取得するステップと、前記受信データ及び前記少なくとも１つの合成データを用いて前記反射体の画像を生成するステップと、を有し、前記動き補正を行うステップは、前記受信データに基づいて前記反射体に対する複数の画像フレームを取得するステップと、前記複数の画像フレームを互いに比較し、前記反射体の動きを検出するステップと、前記反射体の動きに基づいて前記トランスデューサー素子のチャネルごとのデータのディレイを調整することにより、前記動き補正を行うステップと、を含み、前記合成データを取得するステップは、前記動き補正が行われたトランスデューサー素子のチャネルごとのデータに基づいて前記合成データを取得することを特徴とする。

前記合成データを取得するステップは、前記受信データに前記合成フォーカシング方法を用いて合成データセットを生成するステップを含み、前記反射体の画像を生成するステップは、前記合成データセットを用いて前記反射体の画像を生成するステップを含むことが好ましい。
前記合成データを取得するステップは、前記反射体の周辺領域に前記合成フォーカシング地点を決定するステップと、前記合成フォーカシング地点と前記実際フォーカシング地点との間の第１ディレイを予測するステップと、を含むことが好ましい。
前記合成データを取得するステップは、前記第１ディレイに基づいて少なくとも１つのトランスデューサー素子それぞれにおけるディレイを予測するステップと、前記予測された前記少なくとも１つのトランスデューサー素子それぞれにおけるディレイに基づいて、前記受信データをトランスデューサー素子チャネルごとにディレイ補正するステップと、をさらに含むことが好ましい。
前記合成データを取得するステップは、前記反射体の周辺領域内の複数の合成フォーカシング地点を決定するステップと、前記複数の合成フォーカシング地点それぞれに対してディレイ補正を行って、前記複数のフォーカシング地点に対応する合成データセットを取得するステップをさらに含むことが好ましい。
前記反射体の画像を生成するステップは、前記単一の受信データに対して収差補正するステップをさらに含むことが好ましい。

前記反射体の画像を生成するステップは、前記少なくとも１つの合成データ及び前記受信データを用いて前記反射体の画像を生成することが好ましい。
前記受信データに、フィルタリング又は復調の少なくとも１つを行うステップをさらに有することが好ましい。
前記合成フォーカシング方法又は前記動き補正の少なくとも１つを用いてディレイエラーを再算出するステップと、前記合成フォーカシングを繰り返し処理するステップと、をさらに有することが好ましい。

上記目的を達成するためになされた本発明による超音波測定装置は、実際フォーカシング地点に超音波信号をビームフォーカシングし、ビームフォーカシングに対応して、反射体から反射されるアコースティックフィールドを受信してトランスデューサー素子のチャネルごとのデータを含む受信データを取得するトランスデューサーと、前記受信データに基づいて複数の画像フレームを生成し、前記複数の画像フレームに基づいて前記トランスデューサー素子のチャネルごとのデータに動き補正を行う動き補正部と、前記実際フォーカシング地点と前記反射体の周辺領域の合成フォーカシング地点との間のディレイに基づいて前記合成フォーカシング地点に合成フォーカシング方法を適用し、前記合成フォーカシング地点に対応する少なくとも１つの合成データを取得する合成データ取得部と、前記受信データ及び前記少なくとも１つの合成データを用いて、前記反射体の画像を生成する画像生成部と、を備え、前記動き補正部は、前記受信データに基づいて前記反射体に対する複数の画像フレームを取得し、前記複数の画像フレームを互いに比較し、前記反射体の動きを検出し、前記反射体の動きに基づいて前記トランスデューサー素子のチャネルごとのデータのディレイを調整することにより、前記動き補正を行い、前記合成データ取得部は、前記動き補正が行われたトランスデューサー素子のチャネルごとのデータに基づいて前記合成データを取得することを特徴とする。

前記合成フォーカシング地点それぞれに対する合成データをチャネルごとに合算して単一の受信データを生成し、前記単一の受信データに対して収差補正するデータ合算及び収差補正部をさらに備えることが好ましい。

本発明に係る超音波測定装置及び超音波測定方法によれば、相対的に少ない実際のフォーカシング地点についてのみビームフォーカシングを行っても、合成フォーカシング方法を用いて測定地点周辺領域の相対的に多くのフォーカシング地点に対する受信データを取得することができる。そのため、演算時間及び測定時間を著しく縮小することができる、という効果がある。

本発明の実施形態に係る超音波測定装置の概略を示すブロック図である。 本発明の実施形態に係る超音波測定装置の概略を示すブロック図である。 サンプル内に単一反射体が存在する場合を説明するための概念図である。 サンプル内の複数反射体が含まれた場合のアコースティックフィールドの概念図である。 本発明の実施形態に係る超音波測定装置で、比較例による観測対象である第１反射体以外の他の反射体が含まれたサンプルに対する超音波測定方法を説明するためのフローチャートである。 図３に示す超音波測定方法を説明するための概念図である。 図３に示す超音波測定方法を説明するための概念図である。 受信データ及び受信データのチャネルごとの合算を説明するための図である。 受信データ及び受信データのチャネルごとの合算を説明するための図である。 受信データのチャネルごとの合算による他反射体の影響除去を説明するための概念図である。 受信データのチャネルごとの合算による他反射体の影響除去を説明するための概念図である。 本発明による超音波測定装置の他の実施形態に係る超音波測定方法を説明するためのフローチャートである。 様々な実施形態に係るビームフォーカシング方法を説明するための概念図である。 様々な実施形態に係るビームフォーカシング方法を説明するための概念図である。 フォーカシング地点にビームフォーカシングを行った構成を説明するための概念図である。 フォーカシング地点にビームフォーカシングを行った構成を説明するための概念図である。 フォーカシング地点にビームフォーカシングを行った構成を説明するための概念図である。 本発明の一実施形態に係るビームフォーカシング方法を説明するためのフローチャートである。 本発明の複数フォーカシング地点に対するビームフォーカシング方法を説明するためのフローチャートである。 本発明の一実施形態に係る実際フォーカシング受信データを説明するための図である。 サンプル内のフォーカシング地点を説明するための概念図である。 サンプル内のフォーカシング地点を説明するための概念図である。 本発明の一実施形態に係る合成データ取得方法を説明するためのフローチャートである。 本発明の他の実施形態に係る複数の合成フォーカシング地点に対する合成データ取得方法を説明するためのフローチャートである。 本発明の一実施形態に係る合成データの例示である。 本発明の他の実施形態に係る受信データ測定時のサンプル内反射体が移動する場合、反射体動き補正を行う超音波測定方法を説明するためのフローチャートである。 本発明の他の実施形態に係る動き補正を行う超音波測定方法を説明するためのフローチャートである。 サンプル内の反射体の動きを説明するための概念図である。 本発明の実施形態に係る複数の画像フレームの概念図である。 動きベクトル検出及びディレイ補正工程を説明するための概念図である。 チャネルごとにディレイを補正する工程を説明するための図である。 本発明の一実施形態に係る動き補正を行う超音波測定装置の概略構成を示すブロック図である。 収差補正工程を説明するための概念図である。 収差補正工程を説明するための概念図である。 収差補正部の相互相関に基づく収差補正を説明するための概念図である。 相互相関に基づく収差補正以前の画像写真である。 相互相関に基づく収差補正以後の画像写真である。 本発明のさらに他の実施形態に係る超音波測定方法を説明するためのフローチャートである。

次に、本発明に係る超音波測定装置及び超音波測定方法を実施するための形態の具体例を図面を参照しながら説明する。

図１Ａ及び１Ｂは、本発明の実施形態に係る超音波測定装置の概略を示すブロック図である。
図１Ａに示すように、超音波測定装置１８は、プロセッサ３０、表示部２０、媒質コンテナ４０、送受信部５０、増幅器及びフィルタ７０、及びアナログ−デジタル変換部８０を備える。
送受信部５０は、少なくとも１つのトランスデューサー素子（５１〜５４）を含む。

プロセッサ３０は、超音波測定装置全般の動作を制御する。
プロセッサ３０は、格納部（図示せず）に格納されたプログラム又はアルゴリズムに基づいて超音波測定装置全般の動作を制御する。
プロセッサ３０は、送受信部５０の超音波信号の送信時期及び送信地点を決定する。
詳細には、プロセッサ３０は、超音波信号のビームフォーカシング地点を決定する。

プロセッサ３０は、決定されたビームフォーカシング地点に対応して少なくとも１つのトランスデューサー素子（５１〜５４）の各駆動時期を決定する。
プロセッサ３０は、少なくとも１つのトランスデューサー素子（５１〜５４）の駆動時期に対する電気的な信号を生成し、生成された駆動時期に対する電気的な信号を少なくとも１つのトランスデューサー素子（５１〜５４）のそれぞれに出力する。
例えば、プロセッサ３０は、超音波信号がサンプル６０上の少なくとも１つの地点にビームフォーカスされるよう、少なくとも１つのトランスデューサー素子（５１〜５４）のそれぞれの送信時期を決定する。

少なくとも１つのトランスデューサー素子（５１〜５４）のそれぞれは、電気的な信号を超音波信号に変換する。又は、少なくとも１つのトランスデューサー素子（５１〜５４）それぞれは、受信されたアコースティックフィールドを電気的な信号に逆変換し得る。
少なくとも１つのトランスデューサー素子（５１〜５４）のそれぞれは、プロセッサ３０から受信された電気的な信号に基づいて各駆動時期に超音波信号を送信する。少なくとも１つのトランスデューサー素子（５１〜５４）からの超音波信号は互いに干渉してもよく、特定の地点にビームフォーカスされてもよい。

少なくとも１つのトランスデューサー素子（５１〜５４）は、受信されたアコースティックフィールドを電気的な信号に変換する。一方、ここで、アコースティックフィールドは、サンプル６０内の反射体によって反射又は散乱された媒質上のフィールドである。
変換された電気的な信号は、増幅器及びフィルタ７０によって予め設定された利得に増幅されてフィルタリングされる。
さらに、フィルタリングされた電気的な信号は、アナログ−デジタル変換器８０によってデジタル信号に変換される。本実施形態では、受信されたアコースティックフィールドが変換された電気的な信号は受信データと命名するか、又は、変換されたデジタル信号を受信データとして命名してもよい。
受信データは、少なくとも１つのトランスデューサー素子（５１〜５４）のチャネルごとのデータを含み、これに関連しては後で、表１を参照して詳細に説明する。

少なくとも１つのトランスデューサー素子（５１〜５４）は、一定期間の間にアコースティックフィールドを受信する。
時間の流れに伴って、少なくとも１つのトランスデューサー素子（５１〜５４）のそれぞれに受信されるアコースティックフィールドも変更される。
受信データも一定期間の間の時間によるデータであってもよい。例えば、受信データは、アナログ−デジタル変換器８０のサンプリング期間による時系列的なデータであってもよく、これに関連しては後で、表１を参照して詳細に説明する。

プロセッサ３０は、受信データが入力され、入力された受信データに基づいてサンプル６０に対する画像を生成する。
プロセッサ３０は、生成された画像を表示部２０が表示するように制御してもよい。
一方、プロセッサ３０は、サンプル６０の少なくとも１つの実際フォーカシング地点に超音波信号をビームフォーカシングするように少なくとも１つのトランスデューサー素子（５１〜５４）の送信時期を決定する。
さらに、プロセッサ３０は、少なくとも１つのトランスデューサー素子（５１〜５４）のチャネルごとに受信データを取得する。

プロセッサ３０は、合成フォーカシング方法に基づいて、サンプル６０内の反射体周辺領域の複数の合成フォーカシング地点それぞれに対する合成データを取得する。
ここで、複数の合成データを合成データセットと命名する。
プロセッサ３０は、合成データに基づいて反射体の画像を取得する。上述した工程については詳細に後述する。
媒質コンテナ４０は、超音波信号が移動するための媒質を入れるコンテナであってもよい。
一方、媒質が準備された測定環境として、媒質コンテナ４０は超音波測定装置に含まれないこともある。

図１Ｂは、本発明の他の実施形態に係る超音波測定装置の概略を示すブロック図である。
図１Ｂの実施形態に係る超音波測定装置は、図１Ａの実施形態に係る超音波測定装置とは対照的に受信部９０をさらに備える。
受信部９０は、例えば、少なくとも１つのチャネルを有するハイドロフォン（ｈｙｄｒｏｐｈｏｎｅ）又は追加トランスデューサー素子などで実現され得る。

図２Ａは、サンプル内に単一反射体が存在する場合を説明するための概念図である。
図２Ａに示していないが、少なくとも１つのトランスデューサー素子（ｃｈ１〜ｃｈ５）は、サンプル１０に超音波信号を送信する。
サンプル１０は超音波信号を受信し、これに対応してアコースティックフィールドを生成する。
一方、図２Ａに示すように、サンプル１０内には単一の反射体２０１が含まれてもよい。
反射体２０１は、超音波信号に対応したり、又はビームフォーカシング地点からの超音波信号に対応してアコースティックフィールド２１０を放出する。

少なくとも１つのトランスデューサー素子（ｃｈ１〜ｃｈ５）は、アコースティックフィールド２１０を受信し、これを電気的な信号に変換する。
超音波測定装置は、変換された信号に基づいて反射体２０１の位置を測定し得る。
図２Ａに示すように、単一の反射体２０１によって生成されたアコースティックフィールド２１０は、実質的に球面波の形態に形成される。これはホイヘンスの原理によって反射体２０１を中心に同心円形態のアコースティックフィールドが形成されることから起因する。
単一の反射体２０１が実質的な球面波の形態であるアコースティックフィールド２１０を放出することによって、超音波測定装置は球面波の収差（ａｂｅｒｒａｔｉｏｎ）補正を行って単一の反射体２０１に対する画面を生成する。超音波測定装置が収差補正する構成については後述する。

一方、図２Ｂは、サンプル内の複数反射体が含まれた場合のアコースティックフィールドの概念図である。
図２Ｂに示すように反射体（２０１、２８０、２９０）がサンプル１０内に存在する。
少なくとも１つのトランスデューサー素子（ｃｈ１〜ｃｈ５）は、サンプル１０に超音波信号を送信する。
反射体２０１は、超音波信号を受信してアコースティックフィールド２１０を生成して放出する。反射体２８０は、超音波信号を受信してアコースティックフィールド２２０を生成して放出する。反射体２９０は、超音波信号を受信してアコースティックフィールド２３０を生成して放出する。
複数のアコースティックフィールド（２１０、２２０、２３０）は媒質内で互いに干渉する。

図２Ｂに示すように、複数の反射体（２０１、２８０、２９０）それぞれによって生成されたアコースティックフィールド（２１０、２２０、２３０）は、実質的に球面波の形態に形成される。
ただし、球面波形態のアコースティックフィールド（２１０、２２０、２３０）は互いに干渉され得る。そのため、複数の球面波が干渉された形態のアコースティックフィールドが生成され、少なくとも１つのトランスデューサー素子（ｃｈ１〜ｃｈ５）は、干渉された形態のアコースティックフィールドを受信して電気的な信号に変換する。
トランスデューサー素子（ｃｈ１〜ｃｈ５）に球面波の形態と相違の大きいアコースティックフィールドが入力されることによって、超音波測定装置は、特定の反射体２０１に対する正確な画像を生成することができない。

図３は、本発明の実施形態に係る超音波測定装置で、比較例による観測対象である第１反射体以外の他の反射体が含まれたサンプルに対する超音波測定方法を説明するためのフローチャートである。
図３に示す比較例は、図４Ａ〜４Ｃを参照して詳細に説明することにする。

まず、図３を参照すると、ステップＳ３１０において、超音波測定装置は測定地点を決定する。
ここで、測定地点は、画像生成を行うための地点である。例えば、図４Ａに示すように、サンプル１０内には複数の反射体（４０１、４８０、４９０）が含まれる。
超音波測定方法は、複数の反射体（４０１、４８０、４９０）のうちから第１反射体４０１を測定地点に決定するとする。超音波測定方法は、後述する方法に基づいて第１反射体４０１以外の残りの反射体（４８０、４９０）の測定時に及ぼす影響を抑制したり最小化することができる。

再び図３を参照すると、ステップＳ３２０において、超音波測定方法は、測定地点周辺領域のｎ地点にビームフォーカシングを行う。
超音波測定装置は、複数のトランスデューサー素子を含み得る。超音波測定装置は、複数のトランスデューサー素子それぞれの超音波信号の送信時期を調整する。
超音波測定装置は、複数のトランスデューサー素子それぞれの送信時期を含むトランスデューサー制御信号を出力する。
そのため、超音波測定装置は、サンプル１０内の特定地点に超音波信号を集中するようにするビームフォーカシングを行う。

例えば、図４Ｂに示す第１地点４１１にビームフォーカシングを行う場合に、超音波測定装置は第２トランスデューサー素子ｃｈ２の超音波信号の送信時期を相対的に遅くに調整して、第５トランスデューサー素子ｃｈ５の超音波信号の送信時期を相対的に早く調整する。
例えば、第１地点４１１から第２トランスデューサー素子ｃｈ２までの距離が最短であり、第１地点４１１から第５トランスデューサー素子ｃｈ５までの距離が最長であると仮定する。
第２トランスデューサー素子ｃｈ２からの超音波信号が相対的に遅れた時期に送信され、第５トランスデューサー素子ｃｈ５からの超音波信号が相対的に早い時期に送信されることで、複数の超音波信号が第１地点４１１に同時に集中する。超音波測定装置は、残りのトランスデューサー素子の超音波信号の送信時期を決定して全てのトランスデューサー素子からの超音波信号が第１地点４１１に同時に集中するように制御する。

他の実施形態に係る超音波測定方法として、超音波信号が第２地点（図示せず）から発生するようにビームフォーカシングを行ってもよい。
例えば、超音波測定装置は、第２地点（図示せず）をトランスデューサーの後方の一地点として決定してもよく、第２地点（図示せず）から超音波が拡散されるようにビームフォーカシングを行ってもよい。これについて図７Ｂを参照して詳細に説明する。
それと共に、スキャットラインごとに深度方向にフォーカスされる位置が互いに異なってもよい。

第１地点４１１にビームフォーカシングを行った後、第１反射体４０１によってアコースティックフィールドが生成される。
ステップＳ３３０において、少なくとも１つのトランスデューサー素子（ｃｈ１〜ｃｈ５）はアコースティックフィールドを受信し、受信されたアコースティックフィールドを電気的な信号に変換する。
超音波測定装置は、受信されたアコースティックフィールドに対する受信データを格納する。超音波測定装置は、受信データをチャネルごとにｃｈ１〜ｃｈ５に格納してもよい。
超音波測定装置は、受信データを一定期間の間に格納してもよい。例えば、表１は、第１地点４１１に対応した格納された受信データの例示である。

表１の受信データは、予め設定された期間でサンプリングされて格納されたものとして、「ｔ」が指示する１〜６はサンプリングナンバーである。
格納された受信データは、チャネルごと（ｃｈ１〜ｃｈ５）に区分されて格納される。

ステップＳ３４０において、超音波測定方法は「ｎ」が「Ｎ」を超過するか否かを判断する。
例えば、「Ｎ」が“２５”と設定されてもよい。「ｎ」が「Ｎ」を超過しない場合、超音波測定方法はステップＳ３５０において「ｎ」の値を“１”増加させ、これに対応してステップＳ３６０においてフォーカシング地点を移動させる。
超音波測定装置は、第２地点４１２に対してビームフォーカシングを行う。超音波測定装置は、第２地点４１２に対応して少なくとも１つのトランスデューサー素子（ｃｈ１〜ｃｈ５）の送信時期を再調整し、第２地点４１２に対しビームフォーカシングを行う。

超音波測定装置は、測定地点４０１の周辺地域にフォーカシング地点（４１１〜４３５）を決定してもよく、例えば、図４Ｂに示すように、フォーカシング地点（４１１〜４３５）がそれぞれ予め設定された間隔に離隔して全体的に長方形の格子形態を有するようビームフォーカシング地点を決定してもよい。
すなわち、一部のビームフォーカシング地点は、第１反射体４０１よりトランスデューサーに近く配置し、一部のビームフォーカシング地点は第１反射体４０１よりトランスデューサーから遠く配置し、又は、一部のビームフォーカシング地点は。第１反射体４０１とトランスデューサーから同じ距離に配置する。
一方、図４Ｂでは、２次元の地点に２５個のフォーカシング地点に対してビームフォーカシングするように図に示しているが、これは一例に過ぎず、他の実施形態に係る超音波測定装置は、３次元の例を挙げて１２５個フォーカシング地点に対してビームフォーカシングを行ってもよい。

超音波測定装置は、決定された全てのフォーカシング地点に対してビームフォーカシングを行い、各フォーカシング地点に対してチャネルごとのデータを格納する。
例えば、超音波測定装置は、図４Ｃに示すような受信データを格納してもよい。
受信データは、フォーカシング地点ごと（ｎ＝１、２、．．．、２５）に格納される。
各フォーカシング地点ごとの受信データは、チャネルごとに格納される。
第１フォーカシング地点４１１に対する受信データは、第１チャネル〜第５チャネル（４４１〜４４５）ごとに格納される。
さらに、第１チャネル〜第５チャネル（４４１〜４４５）ごとのデータそれぞれは、表１のように予め設定された時間の間に格納される。第２フォーカシング地点４１１に対する受信データも第１チャネル〜第５チャネル（４５１〜４５５）ごとに格納され、第２５フォーカシング地点４３５に対する受信データも第１チャネル〜第５チャネル（４６１〜４６５）ごとに格納される。

全てのフォーカシング地点（４１１〜４３５）に対して測定が完了すれば（ｎ＞Ｎ）、ステップＳ３７０において超音波測定装置はチャネルごとに受信データを合算する。
例えば、図４Ｄに示すように、超音波測定装置は、チャネルごとに受信データを合算する。
図４Ｄに示すように、超音波測定装置は、複数のフォーカシング地点（４１１〜４３５）に対する受信データの第１チャネルデータ（４４１、４５１、．．．、４６１）を合算し、合算された第１チャネルデータ４７１を算出する。

さらに、超音波測定装置は、複数のフォーカシング地点に対する受信データをチャネルごとに合算してチャネルごとの合算データ（４７１〜４７５）を生成する。
そのため、超音波測定装置は、単一の受信データを生成する。
ここで、単一の受信データとは、複数フォーカシング地点に対する受信データがチャネルごとに合算されて１つの単一の受信データが生成されることを示す。
チャネルごとのデータ合算を行うために、観測地点（第１反射体４０１）以外の他の反射体（４８０、４９０）によるアコースティックフィールドの影響が相殺され得る。
これについて、図４Ｅ及び４Ｆを参照して説明することにする。

ステップＳ３８０において、超音波測定装置は収差補正をさらに行う。
ステップＳ３９０において、超音波測定装置は、チャネルごとの合算データに基づいて画像を生成する。

図４Ｅ及び４Ｆは、チャネルごとのデータ合算によって他の反射体による影響が相殺される構成を説明するための概念図である。
図４Ｅは、フォーカシング地点が第１地点４１１である場合に対する概念図である。
第１反射体４０１及び第２反射体４８０のそれぞれからトランスデューサー素子に受信される時間はそれぞれ「ｔ１」及び「ｔ２」である。
さらに、図４Ｆは、フォーカシング地点が第２地点４１２である場合に対する概念図である。
第１反射体４０１及び第２反射体４８０のそれぞれからトランスデューサー素子に受信される時間はそれぞれ「ｔ３」及び「ｔ４」である。ここで、“ｔ１−ｔ２”と“ｔ３−ｔ４”は互いに異なるが、チャネルごとのデータ合算時の観測地点ではない他の反射体からの信号はディレイ差で相殺される。
一方、図３に示す比較例による超音波測定方法では、複数のフォーカシング地点に対してビームフォーカシングするための測定時間が増加する。

図５は、本発明による超音波測定装置の他の実施形態に係る超音波測定方法を説明するためのフローチャートである。
ステップＳ５１０において、超音波測定装置は、少なくとも１つの実際フォーカシング地点に超音波信号をビームフォーカシングする。
特に、本実施形態に係る超音波測定方法では、トランスデューサー素子の個数以下の個数の実際フォーカシング地点にビームフォーカシングを行う。
図６Ａ及び６Ｂは、様々な実施形態に係るビームフォーカシング方法を説明するための概念図である。
図６Ａ及び６Ｂの実施形態で、サンプル１０内の第１地点６０１が測定地点に決定されたと仮定する。

本実施形態に係る超音波測定方法は、第１地点６０１の周辺領域に複数のフォーカシング地点（６１１〜６１５）を決定する。
超音波測定装置は、フロントフォーカシング方法に基づいて、トランスデューサー素子（ｃｈ１〜ｃｈ５）の前面に複数のフォーカシング地点（６１１〜６１５）を決定する。
超音波測定装置は、トランスデューサー素子の個数である“５つ”より以下の個数でフォーカシング地点（６１１〜６１５）の個数を決定する。
他の実施形態に係る超音波測定方法では、リアフォーカシング方法に基づいてトランスデューサー素子（ｃｈ１〜ｃｈ５）の後面に複数のフォーカシング地点（６２１〜６２５）を決定する。

図７Ａ〜７Ｃは、フォーカシング地点にビームフォーカシングを行った構成を説明するための概念図である。
図７Ａは、スキャンライン（ｓｃａｎ ｌｉｎｅ）ごとのフォーカシング方法に基づいたビームフォーカシング実行結果であり、図７Ｂはリアフォーカシング方法に基づいたビームフォーカシング実行結果であり、図７Ｃはフロントフォーカシング方法に基づいたビームフォーカシング実行結果である。

図８は、本発明の実施形態に係るビームフォーカシング方法を説明するためのフローチャートである。
ステップＳ８０５において、超音波測定装置は、例えば、図６Ａ及び図６Ｂに示す測定地点６０１のように測定地点を決定する。
ステップＳ８１０において、超音波測定装置は、少なくとも１つの実際フォーカシング地点を決定する。
ステップＳ８２０において、超音波測定装置は、少なくとも１つのトランスデューサーの送信時期を決定、すなわち、トランスデューサー素子間の駆動ディレイを決定する。
ステップＳ８３０において、超音波測定装置は、少なくとも１つのトランスデューサー素子を決定された送信時期に超音波信号を送信するように制御し、そのため、特定実際フォーカシング地点にビームフォーカシングを行う。

図９は、本発明の実施形態による複数フォーカシング地点に対するビームフォーカシング方法を説明するためのフローチャートである。
ステップＳ９１０において、超音波測定装置は測定地点を決定する。
ステップＳ９２０において、超音波測定装置は、少なくとも１つの実際フォーカシング地点を決定する。例えば、超音波測定装置は、図６Ａ又は６Ｂに示すように、フロントフォーカシング方式又はリアフォーカシング方式に基づいて少なくとも１つの実際フォーカシング地点を決定する。
ステップＳ９３０において、超音波測定方法は、「ｋ」地点にビームフォーカシングする。
これに対応して、ステップＳ９４０において、アコースティックフィールドを受信して格納する。
「ｋ」は「Ｋ」値になるまでフォーカシング地点を移動させながらビームフォーカシングを行い、「ｋ」が「Ｋ」値以下であれば（ステップＳ９５０）、超音波測定装置はステップＳ９６０で「ｋ」値を１ずつ増加させ、ステップＳ９７０でフォーカシング地点を移動させる。

再び図５を参照すると、ステップＳ５２０において、超音波測定装置は、フォーカシング地点それぞれに対応するアコースティックフィールドを受信して格納する。
ステップＳ５３０において、超音波測定装置は、受信データそれぞれのディレイを調整して合成フォーカシングを行う。
合成フォーカシングは、サンプル上の第１実際フォーカシング地点に対応した受信データを用いて、第２合成フォーカシング地点に対応する合成データを生成する方法である。
詳細には、合成フォーカシングは、第１実際フォーカシング地点にビームフォーカシングを行った後、受信されたチャネルごとのデータそれぞれのディレイを調整し、第２合成フォーカシング地点にビームフォーカシングを行い、これに対応して受信されたデータを生成する方法である。

例えば、超音波測定装置は、第１フォーカシング地点に実際ビームフォーカシングを実行して、これに対応する第１受信データを取得する。ここで、第１実際フォーカシング地点は複数であってもよく、第１受信データも複数であってもよい。
超音波測定装置は、実際ビームフォーカシングを行った第１実際フォーカシング地点と合成フォーカシングを行う第２合成フォーカシング地点との間の距離、及びディレイのうち少なくとも１つを判断する。
さらに、超音波測定装置は、第１実際フォーカシング地点と第２合成フォーカシング地点との間のディレイに基づいて、トランスデューサー素子チャネルごとのディレイを予測する。

超音波測定装置は、予測されたトランスデューサー素子チャネルごとのディレイに基づいて、第１受信データのチャネルごとのデータをディレイ調整する。
そのため、超音波測定装置は、第２合成フォーカシング地点にビームフォーカシングを行った場合に受信される受信データの第２受信データを取得する。
すなわち、超音波測定装置は、第２合成フォーカシング地点に実際にビームフォーカシングすることなく、第２合成フォーカシング地点に対応した受信データを取得する。

ステップＳ５４０において、超音波測定装置は、測定地点周辺領域の複数の合成フォーカシング地点に対する合成データを取得する。
合成データを取得する方法の詳細は後述する。
ステップＳ５５０において、超音波測定装置は、複数の合成データのチャネルごとのデータを合算する。
チャネルごとのデータ合算にしたがって、サンプル内に複数の反射体がある場合にも、測定地点以外の他の反射体の影響が抑制される。
他の実施形態で、超音波測定装置は、合成データ及び受信データに対してチャネルごとのデータを合算してもよい。

ステップＳ５６０において、超音波測定装置は、チャネルごとの合算データに基づいて測定地点に対する画像を生成する。
図５に示していないが、超音波測定装置は、収差補正を行ってチャネルごとの合算データに基づいて測定地点に対する画像を生成してもよい。
又、超音波測定装置は、受信データからハーモニック（ｈａｒｍｏｎｉｃ）成分を抽出して用いてもよい。
さらに、超音波測定装置は、受信データに対してフィルタリング又は復調を行って、受信端ディレイ補正又は合成フォーカシング方法を適用してもよい。
詳細には、受信データは、フィルタリングを用いてバンドパスフィルタリングされたデータに加工されてハーモニック成分を抽出してもよく、パルスインバージョン（ｐｕｌｓｅ ｉｎｖｅｒｓｉｏｎ）のような方式を用いてハーモニック成分が抽出されてもよく、復調によりベースバンド信号に変換されてもよい。

図１０は、本発明の一実施形態に係る実際フォーカシング受信データを説明するための図である。
例えば、超音波測定装置は、図１０に示すような実際フォーカシング受信データを取得する。
ここで、実際フォーカシング受信データは、超音波測定装置が実際に超音波信号をサンプルに送信し、これに対応してアコースティックフィールドを受信した場合のデータを示す。
実際フォーカシング受信データは、図１０に示すように、チャネルごとに受信データを含む。

図１０に示す実際フォーカシング受信データは、例えば、図６Ａに示すように、５個のフォーカシング地点（６１１〜６１５）にビームフォーカシングを行って取得された受信データであり得る。
第１地点６１１に対するビームフォーカシングの結果、受信された第１データは、第１チャネル受信データ乃至第５チャネル受信データ（１００１〜１００５）を含む。
第２地点６１２に対するビームフォーカシングの結果、受信された第２データは、第１チャネル受信データ乃至第５チャネル受信データ（１０１１〜１０１５）を含む。
フォーカシングの結果、受信された第５データは、第１チャネル受信データ乃至第５チャネル受信データ（１０２１〜１０２５）を含む。

図１１Ａ及び１１Ｂは、サンプル内の合成フォーカシング地点を説明するための概念図であり、図１２は、本発明の一実施形態に係る合成データ取得方法を説明するためのフローチャートである。
図１１Ａ及び１１Ｂは、図１２を参照してより詳しく説明する。
まず、図１２を参照すると、ステップＳ１２１０において、超音波測定装置は合成フォーカシング地点を決定する。
例えば、超音波測定装置は、図１１Ａに示すように、合成フォーカシング地点（１１０１〜１１２０）を決定する。超音波測定装置は、測定地点６０１の周辺領域に合成フォーカシング地点を決定する。
例えば、超音波測定装置は、フォーカシング地点（１１０１〜１１２０）がそれぞれ予め設定された間隔で離隔して全体的に長方形の格子形態を有するように合成フォーカシング地点を決定してもよい。

ステップＳ１２２０において、超音波測定装置は、合成フォーカシング地点（１１０１〜１１２０）それぞれに対応する受信データ取得のために格納された受信データのディレイを調整し、ステップＳ１２３０において、測定地点周辺領域の合成データを取得する。
例えば、図１１Ｂに示すように、超音波測定装置は、実際フォーカシング地点（６１１〜６１５）と合成フォーカシング地点１１０１との間のディレイがそれぞれＡ１、Ａ２、Ａ３、Ａ４、Ａ５であることを取得する。
超音波測定装置は、合成フォーカシング地点１１０１と実際フォーカシング地点（６１１〜６１５）との間で取得されたディレイ、及び少なくとも１つのトランスデューサー素子それぞれにおけるディレイに基づいて合成フォーカシング方法を行う。

詳細には、合成フォーカシング地点１１０１に対応する合成データは、図１０に示すような実際フォーカシング受信データに基づいて取得される。
例えば、トランスデューサー素子から実際フォーカシング地点（６１１〜６１５）までのフォーカシング時間を「ｔ１」と仮定する。さらに、上述したように、実際フォーカシング地点（６１１〜６１５）と合成フォーカシング地点１１０１との間のディレイはそれぞれＡ１、Ａ２、Ａ３、Ａ４、Ａ５である。

超音波測定装置は、第１実際フォーカシング地点に対応する受信データに“ｔ１−Ａ１”のディレイを適用させる。さらに、超音波測定装置は、第２実際フォーカシング地点から第５実際フォーカシング地点に対応する受信データそれぞれに“ｔ１−Ａ２”〜“ｔ１−Ａ５”のディレイを適用させる。
超音波測定装置は、ディレイが適用された第１実際フォーカシング地点に対応する受信データ乃至ディレイが適用された第５実際フォーカシング地点に対応する受信データを合算し、合成フォーカシング地点１１０１に対応する合成データを取得する。

一方、上述した“ｔ１−Ａ１”のディレイは、合成フォーカシング地点１１０１が実際フォーカシング地点（６１１〜６１５）より相対的にトランスデューサー素子に近い場合の例示である。
実際フォーカシング地点（６１１〜６１５）が合成フォーカシング地点１１０１より相対的にトランスデューサー素子に近い場合、ディレイは“ｔ１＋Ａ１”に決定される。
超音波測定装置は、上述したように合成データを取得し、そのため、図４Ｃに示すような、実際フォーカシング受信データと同一のデータを取得することができる。

すなわち、超音波測定装置は、合成フォーカシング方法を受信データに適用し、図４Ｂに示すように実際フォーカシングを行うことなく、図４Ｃに示すような受信データを取得することができる。
上述したように、本発明の一実施形態に係る超音波測定方法は、相対的に少ない実際のフォーカシング地点についてのみビームフォーカシングを行っても、合成フォーカシング方法を用いて測定地点周辺領域の相対的に多くのフォーカシング地点に対する受信データを取得することができる。
そのため、演算時間及び測定時間を著しく縮小することができる。

図１３は、本発明の他の実施形態に係る複数の合成フォーカシング地点に対する合成データを取得する方法を説明するためのフローチャートである。
図１３に示す実施形態は、超音波測定装置が実際フォーカシング地点からの受信データを格納したと仮定する。
ステップＳ１３１０において、超音波測定装置は複数の合成フォーカシング地点を決定する。
例えば、超音波測定装置は、測定地点近所領域で複数のフォーカシング地点を決定する。

ステップＳ１３２０において、超音波測定装置は、合成フォーカシング地点と実際フォーカシング地点との間のディレイを予測する。
超音波測定装置は、「ｍ」合成フォーカシング地点と実際フォーカシング地点との間の距離を予測する。
超音波測定装置は、「ｍ」合成フォーカシング地点と少なくとも１つのトランスデューサー素子までの測定時間と、実際フォーカシング地点と少なくとも１つのトランスデューサー素子までの測定時間との間の差に基づいて、トランスデューサー素子チャネルごとにディレイを予測する。そのため、トランスデューサー素子チャネルごとのディレイを予測する。

ステップＳ１３３０において、超音波測定装置は、予測されたディレイに基づいて受信データのディレイを補正する。
合わせて、ステップＳ１３４０において、超音波測定装置は、「ｍ」合成フォーカス地点に対応する合成データを取得して格納する。
超音波測定装置は、「ｍ」が「Ｍ」を超過するまで上述した工程を繰り返す（ステップＳ１３５０、Ｓ１３６０）。ここで、「Ｍ」は合成フォーカシング地点の個数である。
超音波測定装置は、ステップＳ１３７０において、測定地点周辺領域の複数フォーカシング地点に対応する合成データを取得する。

図１４は、本発明の一実施形態に係る合成データの例示である。
図１４に示すように、合成データは、例えば、２５個であってもよく、それぞれの合成データはチャネルごとのデータを含んでもよい。
例えば、第１合成データは、第１チャネルデータ乃至第５チャネルデータ（１４０１〜１４０５）を含み、第２合成データは第１チャネルデータ乃至第５チャネルデータ（１４１１〜１４１５）を含み、第３合成データは第１チャネルデータ乃至第５チャネルデータを含む。
合成データは、実際受信データを含んでもよい。例えば、超音波測定装置は、実際受信データ５個及び合成データ２０個を生成して２５個の受信データを生成してもよい。
又は、超音波測定装置は、実際受信データ５個に基づいて合成データ２５個を生成して２５個の受信データを生成してもよい。

図１５は、本発明の他の実施形態に係る超音波測定方法を説明するためのフローチャートである。
図１５は、特に受信データ測定時のサンプル内反射体が移動する場合、反射体動き補正を行う超音波測定方法を説明するためのフローチャートである。
ステップＳ１５１０において、超音波測定装置はサンプル上の複数地点に超音波信号ビームフォーカシングを行う。

例えば、超音波測定装置は、図６Ａ又は６Ｂに示すような、フロントフォーカシング方法又はリアフォーカシング方法に基づいてビームフォーカシングを行う。
超音波測定装置は、サンプル上の複数地点に送信側ビームフォーカシングを行う。
超音波測定装置は、トランスデューサー素子の個数よりも小さい個数のフォーカシング地点にビームフォーカシングを行う。

ステップＳ１５２０において、超音波測定装置は、サンプル上の複数地点に対するビームフォーカシングのそれぞれに対応して受信されるアコースティックフィールドを受信する。
超音波測定装置は、アコースティックフィールドを受信して受信データを取得及び格納する。
受信データは、トランスデューサー素子のチャネルごとのデータを含んでもよく、予め設定された時間の間測定されてもよい。

ステップＳ１５３０において、超音波測定装置は、チャネルごとのデータに対して動き補正を行う。
詳細には、超音波測定装置は、受信データに基づいて画像を生成する。超音波測定装置は、時間の流れに応じて複数の画像フレームを生成する。超音波測定装置は、複数の画像フレームのうち互いに隣接する画像フレームを比較して反射体の動きを検出する。
超音波測定装置は、反射体動きに対応してチャネルごとのデータを補正する。

例えば、超音波測定装置は、反射体が動く以前の位置と反射体が動いた以後の位置との間のディレイを予測する。
超音波測定装置は、予測されたディレイに基づいてチャネルごとのデータを補正して動き補正を行う。
動きベクトル検出は、画面次元で行われてもよく、動き補正はチャネルごとのデータのディレイを調整することによって行われてもよい。

ステップＳ１５４０において、超音波測定装置は合成データを取得する。
例えば、超音波測定装置はフォーカシング地点を決定し、合成フォーカシング地点と実際フォーカシング地点との間のディレイを予測する。
超音波測定装置は、合成フォーカシング地点と実際フォーカシング地点との間のディレイに基づいてチャネルごとのデータを補正し、そのため、合成フォーカシング地点に対応する合成データを取得する。

ステップＳ１５５０において、超音波測定装置は、合成データ又は実際フォーカシング受信データ及び合成データ両者に基づいて画像を生成する。
超音波測定装置は、チャネルごとにチャネルごとのデータを合算する。
上述した構成によって測定地点以外の他の反射体の影響を予防することができる。

また、超音波測定装置は、収差補正（ａｂｅｒｒａｔｉｏｎ ｃｏｒｒｅｃｔｉｏｎ）を行ってもよい。
超音波測定装置は、収差補正された合算データに基づいて画像を生成する。
また、超音波測定装置は、合成フォーカシング又は動き補正によってディレイエラーを再算出してもよく、合成フォーカシングを再び行うことによって上記のステップが繰り返すように制御する。

図１６は、本発明の他の実施形態に係る動き補正を行う超音波測定方法を説明するためのフローチャートである。
図１６に示す動き補正の方法は、図１７Ａ〜１７Ｃと図１８を参照してより詳しく説明する。
図１６に示す動き補正の方法の実施形態で、サンプル１０内の反射体１７００は、図１７Ａに示すように移動するものと仮定する。

ステップＳ１６１０において、超音波測定装置は、サンプル上の複数地点に超音波信号ビームフォーカシングを行う。
例えば、超音波測定装置は、図６Ａ又は図６Ｂに示すような、フロントフォーカシング方法又はリアフォーカシング方法に基づいてビームフォーカシングを行う。
ステップＳ１６２０において、超音波測定装置は、サンプル上の複数地点に対するビームフォーカシングそれぞれに対応して受信されるアコースティックフィールドを受信する。
超音波測定装置は、アコースティックフィールドを受信して受信データを取得及び格納を行う。受信データはトランスデューサー素子のチャネルごとのデータを含んでもよく、予め設定された時間の間測定されてもよい。

ステップＳ１６３０において、超音波測定装置は、受信データに基づいて画像フレームを生成する。
上述したように、受信データは予め設定された時間の間に測定され、そのために時間の流れに応じて格納される。
超音波測定装置は、時間の流れにより格納された受信データに基づいて予め設定された時間の間の複数の画像フレームを生成する。
例えば、超音波測定装置は、図１７Ｂに示すように、第１、２、３画像フレーム（１７１０、１７２０、１７３０）を生成する。

超音波測定装置は、第１、２、３画像フレーム（１７１０、１７２０、１７３０）での反射体の動きを検出する。
例えば、超音波測定装置は、第１画像フレーム１７１０での反射体１７１１と第２画像フレーム１７２０での反射体１７２１との間の距離が閾値以上であることを確認する。
ステップＳ１６４０において、超音波測定装置は、隣接する画像フレームの反射体位置の間の距離が閾値以上であるかを判断して動きベクトルを抽出する。
さらに、超音波測定装置は、第３画像フレーム１７３０での反射体１７３１も第１画像フレーム１７１０と比較して動きを確認する。

例えば、超音波測定装置は、第３画像フレーム１７３０と第２画像フレーム１７２０とを比較し、動きベクトルが「０」であることを確認する。
これによって、超音波測定装置は、第３画像フレーム１７３０の反射体１７３１と第１画像フレーム１７１０の反射体１７１１との位置間の動きベクトルと、第２画像フレーム１７２０の反射体１７２１と第１画像フレーム１７１０の反射体１７１１との位置間の動きベクトルとが同一であることを確認できる。

一方、図に示していないが、第４画像フレームで反射体が再び動いてもよい。
このような場合、超音波測定装置は、第４画像フレームと第３画像フレームとを比較して動きベクトルを検出する。
超音波測定装置は、第４画像フレームでの反射体と第１画像フレームでの反射体との間の動きベクトルを取得する。
詳細には、超音波測定装置は、第３フレームと第１フレームとの間の動きベクトルに、第４フレームと第３フレームとの間の動きベクトルを合算し、第４フレームと第１フレームとの間の動きベクトルを取得する。

ステップＳ１６５０において、超音波測定装置は、動きベクトルに基づいて受信データの動きを補正する。
例えば、超音波測定装置は、図１７Ｃに示すように、移動以前地点１７４１と移動以後地点１７４２との間のディレイ「Δｔ」を予測する。
超音波測定装置は、移動以前地点１７４１と移動以後地点１７４２との間のディレイ「Δｔ」に基づいて、トランスデューサー素子チャネルごとのディレイを予測することができる。超音波測定装置は、予測されたトランスデューサー素子チャネルごとのディレイに基づいて、受信データをチャネルごとにディレイ補正を行う。それにより、受信データは動き補正される。

図１８は、チャネルごとにディレイを補正する工程を説明するための図である。
図１８は、例えば、第１地点にビームフォーカシングを行った結果に応じて取得された受信データである。
受信データはｔ＝１、２、．．．、Ｔである場合に対して格納される。ここで、「ｔ」はサンプリング区間に対する時間インデックスである。さらに、各時間ごとの受信データは、チャネルごとのデータを含む。
例えば、ｔ＝１である場合の受信データは、第１チャネル受信データ乃至第５チャネル受信データ（１８０１〜１８０５）を含む。ｔ＝２である場合の受信データは、第１チャネル受信データ乃至第５チャネル受信データ（１８１１〜１８１５）を含んでもよい。ｔ＝Ｔである場合の受信データは、第１チャネル受信データ乃至第５チャネル受信データ（１８２１〜１８２５）を含む。

一方、例えば、図１７Ａ及び１７Ｂに示す実施形態では、第１サンプリング区間と第２サンプリング区間との間でのみ反射体が移動した。さらに、第２サンプリング区間以後の第Ｔサンプリング区間まで反射体が停止した状況を想定する。
この場合、超音波測定装置は、第２サンプリング区間に対応する受信データから第Ｔサンプリング区間に対応する受信データのチャネルごとのディレイ補正を行う。
詳細には、超音波測定装置は、各サンプリング区間でチャネルごとのディレイを予測し、予測されたチャネルごとのディレイに基づいてチャネルごとのディレイ補正を行う。
そのため、超音波測定装置は、第２サンプリング区間に対応するデータ乃至第Ｔサンプリング区間に対応する受信データが、反射体が初期位置に継続して位置したような受信データになるよう補正する。

例えば、超音波測定装置は、図１７Ａに示すように、反射体１７００が初期位置に続けて位置するように受信データを補正する。
それにより、超音波測定装置は、図１８に示すように、第２サンプリング区間のチャネルごとのデータを動き補正した、動き補正された受信データ（１８３１〜１８３５）と、第Ｔサンプリング区間のチャネルごとのデータを動き補正した、動き補正された受信データ（１８４１〜１８４５）を生成する。
従って、第２サンプリング区間の受信データと第Ｔサンプリング区間の受信データとの間の受信データも動き補正されたことを、当業者は容易に理解するのであろう。

再び図１６を参照すると、ステップＳ１６６０において、超音波測定装置は、動き補正された受信データに基づいて合成データを取得する。
また、ステップＳ１６７０において、超音波測定装置は、合成データに基づいて画像を生成する。
超音波測定装置は、チャネルごとにチャネルごとのデータを合算する。
上述した構成によって測定地点以外の他の反射体の影響を予防することができる。
超音波測定装置は、収差補正を行ってもよい。超音波測定装置は、収差補正された合算データに基づいて画像を生成してもよい。
上述した工程によって、サンプル内に動いている反射体が存在する場合にも動きを補正することによって明確な画像生成を行うことができる。

図１９は、本発明の一実施形態に係る動き補正を行う超音波測定装置の概略構成を示すブロック図である。
図１９に示すように、超音波測定装置は、超音波信号送受信部１９１０、動き補正部１９２０、合成データ取得部１９３０、データ合算及び収差補正部１９４０、画像生成部１９５０、動き推定部１９６０、及び表示部１９７０を備える。

超音波信号送受信部１９１０は、サンプルの少なくとも１つの実際フォーカシング地点に超音波信号を送信し、ビームフォーカシングに対応して反射体が生成するアコースティックフィールドを受信する。

画像生成部１９５０は、受信したアコースティックフィールドに対応する受信データに基づいて複数の画像フレームを生成する。
受信データは、上述したように時系列的なデータであってもよい。例えば、受信データは表１のように、サンプリングナンバー及びトランスデューサー素子チャネルごとのデータを含み得る。そのため、画像生成部１９５０は、時間の流れに応じた複数の画像フレームを生成することができる。

動き補正部１９２０は、生成された複数の画像フレーム間の比較によって動く反射体を補正する。
詳細には、動き推定部１９６０で推定された動きベクトルに基づいて動き補正部１９２０は受信データをチャネルごとに動き補正を行う。
合成データ取得部１９３０は、動き補正された受信データに合成フォーカシング方法を適用する。そのため、合成データ取得部１９３０は、少なくとも１つの合成データを取得する。

データ合算及び収差補正部１９４０は、少なくとも１つの合成データをチャネルごとに合算して単一の受信データを取得する。
さらに、データ合算及び収差補正部１９４０は、単一の受信データに収差補正を行う。
画像生成部１９５０は、収差補正された単一の受信データに基づいて画像を生成する。
表示部１９７０は、生成された画像を表示する。

図２０Ａ及び２０Ｂは、収差補正工程を説明するための概念図である。
図２０Ａに示すように、（ｘ、ｚ）座標を有する反射体２００１からアコースティックフィールド２０４１が放出される。
反射体２００１は、トランスデューサー素子アレイ２００２から“Ｒ”だけ離隔される。ここで、アコースティックフィールド２０４１は球面波の形態であってもよく、これによってトランスデューサー素子それぞれに到達する時間が異なってもよい。

トランスデューサー素子は、受信したアコースティックフィールド２０４１を電気的な信号２００３に変換して出力する。
収差補正部２００４は、受信した電気的な信号２００３のディレイを調整して収差補正された電気的な信号２０１０を出力する。
信号合算部２００５は、収差補正された電気的な信号２０１０を合算して画像信号２００６を生成する。

一方、図２０Ｂに示すように、音速度が一定ではない媒質２０３０が配置されている場合、アコースティックフィールド２０４２がトランスデューサー素子のそれぞれに球面波の形態ではない不規則な時間に到達する。
そのため、変換された電気的な信号２０１３に対して収差補正部２００４が収差補正しても収差補正された信号２０４０の間にディレイが発生することがある。
そのため、画像信号２０６０は、反射体２００１に対する歪曲された画像信号となる。

本発明の一実施形態に係る収差補正部２００４は、中央チャネル信号を他のチャネル信号と相互相関（ｃｒｏｓｓ ｃｏｒｒｅｌａｔｉｏｎ）して他のチャネル信号が中央チャネル信号と整列するように収差補正してもよい。
例えば、収差補正部２００４は、相互相関して図２０Ｃに示すように、他のチャネル信号のディレイ「Δｔ」を調整して中央チャネル信号と整列させる。
図２０Ｄ及び２０Ｅは、本発明の実施形態に係る相互相関に基づいた収差補正以前と以後の画像写真である。

図２１は、本発明のさらに他の実施形態に係る超音波測定方法を説明するためのフローチャートである。
超音波測定装置は、ステップＳ２１１０において、少なくとも１つの実際フォーカシング地点に対応する受信データに合成フォーカシング方法を適用し、反射体周辺領域の合成フォーカシング地点に対応する少なくとも１つの合成データを取得する。
超音波測定装置は、ステップＳ２１２０において、少なくとも１つの合成データを用いて反射体の画像を生成する。

以上で説明した装置は、ハードウェア構成要素、ソフトウェア構成要素、及び／又はハードウェア構成要素及びソフトウェア構成要素の組み合わせで実現され得る。
例えば、プロセッサ、コントローラ、ＡＬＵ（ａｒｉｔｈｍｅｔｉｃ ｌｏｇｉｃ ｕｎｉｔ）、デジタル信号プロセッサ（ｄｉｇｉｔａｌ ｓｉｇｎａｌ ｐｒｏｃｅｓｓｏｒ）、マイクロコンピュータ、ＦＰＡ（ｆｉｅｌｄ ｐｒｏｇｒａｍｍａｂｌｅ ａｒｒａｙ）、ＰＬＵ（ｐｒｏｇｒａｍｍａｂｌｅ ｌｏｇｉｃ ｕｎｉｔ）、マイクロプロセッサー、又は、命令（ｉｎｓｔｒｕｃｔｉｏｎ）を実行して応答できる異なる装置のように、１つ以上の汎用コンピュータ又は特殊目的のコンピュータを用いて実現してもよい。

処理装置は、オペレーションシステム（ＯＳ）及びオペレーションシステム上で行われる１つ以上のソフトウェアアプリケーションを実行し得る。また、処理装置は、ソフトウェアの実行に応答してデータに対しアクセス、格納、操作、処理、及び生成を行い得る。
理解の便宜のために、処理装置は１つの使用されるものとして説明した場合もあるが、当該の技術分野で通常の知識を有する者であるならば、処理装置が複数の処理要素（ｐｒｏｃｅｓｓｉｎｇ ｅｌｅｍｅｎｔ）及び／又は複数の類型の処理要素を含んでいることが理解されよう。例えば、処理装置は、複数のプロセッサ又は１つのプロセッサ及び１つのコントローラを含んでもよい。また、並列プロセッサ（ｐａｒａｌｌｅｌ ｐｒｏｃｅｓｓｏｒ）のような、他の処理構成も可能である。

ソフトウェアは、コンピュータプログラム、コード、命令、又はこのうちの１つ以上の組み合わせを含んでもよく、所望の通りに動作するよう処理装置を構成したり、独立的又は結合的に処理装置を命令してもよい。
ソフトウェア及び／又はデータは、処理装置によって解釈されたり処理装置に命令又はデータを提供するためどのような類型の機械、構成要素、物理的装置、仮想装置、コンピュータ格納媒体又は装置、送信される信号波に永久的又は一時的に具体化できる。
ソフトウェアは、ネットワークに接続されたコンピュータシステム上に分散し、分散された方法で格納されたり実行されてもよい。ソフトウェア及びデータは１つ以上のコンピュータで読取可能な記録媒体に格納されてもよい。

上述の本発明の実施形態に係る方法は、多様なコンピュータ手段を介して様々な処理を実行することができるプログラム命令の形態で実現され、コンピュータ読取可能な記録媒体に記録されてもよい。
コンピュータ読取可能な記録媒体は、プログラム命令、データファイル、データ構造などのうち１つ又はその組み合わせを含んでもよい。記録媒体に記録されるプログラム命令は、本発明の目的のために特別に設計されて構成されたものでもよく、コンピュータソフトウェア分野の技術を有する当業者にとって公知のものであり、使用可能なものであってもよい。

コンピュータ読取可能な記録媒体の例としては、ハードディスク、フロッピー（登録商標）ディスク及び磁気テープのような磁気媒体、ＣＤ−ＲＯＭ、ＤＶＤのような光記録媒体、光ディスクのような光磁気媒体、及びＲＯＭ、ＲＡＭ、フラッシュメモリなどのようなプログラム命令を保存して実行するように特別に構成されたハードウェア装置が含まれてもよい。
プログラム命令の例には、コンパイラによって作られるような機械語コードだけでなく、インタープリタなどを用いてコンピュータによって実行できる高級言語コードが含まれる。前記したハードウェア装置は、本発明の動作を行うために１つ以上のソフトウェアモジュールとして作動するように構成されてもよく、その逆も同様である。

尚、本発明は、上述の実施形態に限られるものではない。本発明の技術的範囲から逸脱しない範囲内で多様に変更実施することが可能である。

１０、６０ サンプル
１８ 超音波測定装置
２０ 表示部
３０ プロセッサ
４０ 媒質コンテナ
５０ 送受信部
５１〜５４ トランスデューサー素子
７０ 増幅器及びフィルタ
８０ アナログ−デジタル変換部
９０ 受信部
２０１、２８０、２９０ 反射体
２１０、２２０、２３０ アコースティックフィールド

Claims (11)

1. 実際フォーカシング地点に超音波信号をビームフォーカシングするステップと、
   ビームフォーカシングに対応して、反射体から反射されるアコースティックフィールドを受信してトランスデューサー素子のチャネルごとのデータを含む受信データを取得するステップと、
   前記受信データに基づいて複数の画像フレームを生成し、前記複数の画像フレームに基づいて前記トランスデューサー素子のチャネルごとのデータに動き補正を行うステップと、
   前記実際フォーカシング地点と前記反射体の周辺領域の合成フォーカシング地点との間の第１ディレイに基づいて前記合成フォーカシング地点に合成フォーカシング方法を適用し、前記合成フォーカシング地点に対応する少なくとも１つの合成データを取得するステップと、
   前記受信データ及び前記少なくとも１つの合成データを用いて前記反射体の画像を生成するステップと、を有し、
   前記動き補正を行うステップは、前記受信データに基づいて前記反射体に対する複数の画像フレームを取得するステップと、
   前記複数の画像フレームを互いに比較し、前記反射体の動きを検出するステップと、
   前記反射体の動きに基づいて前記トランスデューサー素子のチャネルごとのデータのディレイを調整することにより、前記動き補正を行うステップと、を含み、
   前記合成データを取得するステップは、前記動き補正が行われたトランスデューサー素子のチャネルごとのデータに基づいて前記合成データを取得することを特徴とする超音波測定方法。
2. 前記合成データを取得するステップは、前記受信データに前記合成フォーカシング方法を用いて合成データセットを生成するステップを含み、
   前記反射体の画像を生成するステップは、前記合成データセットを用いて前記反射体の画像を生成するステップを含むことを特徴とする請求項１に記載の超音波測定方法。
3. 前記合成データを取得するステップは、前記反射体の周辺領域に前記合成フォーカシング地点を決定するステップと、
   前記合成フォーカシング地点と前記実際フォーカシング地点との間の第１ディレイを予測するステップと、を含むことを特徴とする請求項１に記載の超音波測定方法。
4. 前記合成データを取得するステップは、前記第１ディレイに基づいて少なくとも１つのトランスデューサー素子それぞれにおけるディレイを予測するステップと、
   前記予測された前記少なくとも１つのトランスデューサー素子それぞれにおけるディレイに基づいて、前記受信データをトランスデューサー素子チャネルごとにディレイ補正するステップと、をさらに含むことを特徴とする請求項３に記載の超音波測定方法。
5. 前記合成データを取得するステップは、
   前記反射体の周辺領域内の複数の合成フォーカシング地点を決定するステップと、
   前記複数の合成フォーカシング地点それぞれに対してディレイ補正を行って、前記複数のフォーカシング地点に対応する合成データセットを取得するステップをさらに含むことを特徴とする請求項４に記載の超音波測定方法。
6. 前記反射体の画像を生成するステップは、単一の前記受信データに対して収差補正するステップをさらに含むことを特徴とする請求項１に記載の超音波測定方法。
7. 前記反射体の画像を生成するステップは、前記少なくとも１つの合成データ及び前記受信データを用いて前記反射体の画像を生成することを特徴とする請求項１に記載の超音波測定方法。
8. 前記受信データに、フィルタリング又は復調の少なくとも１つを行うステップをさらに有することを特徴とする請求項１に記載の超音波測定方法。
9. 前記合成フォーカシング方法又は前記動き補正の少なくとも１つを用いてディレイエラーを再算出するステップと、
   前記合成フォーカシングを繰り返し処理するステップと、をさらに有することを特徴とする請求項１に記載の超音波測定方法。
10. 実際フォーカシング地点に超音波信号をビームフォーカシングし、ビームフォーカシングに対応して、反射体から反射されるアコースティックフィールドを受信してトランスデューサー素子のチャネルごとのデータを含む受信データを取得するトランスデューサーと、
    前記受信データに基づいて複数の画像フレームを生成し、前記複数の画像フレームに基づいて前記トランスデューサー素子のチャネルごとのデータに動き補正を行う動き補正部と、
    前記実際フォーカシング地点と前記反射体の周辺領域の合成フォーカシング地点との間のディレイに基づいて前記合成フォーカシング地点に合成フォーカシング方法を適用し、前記合成フォーカシング地点に対応する少なくとも１つの合成データを取得する合成データ取得部と、
    前記受信データ及び前記少なくとも１つの合成データを用いて、前記反射体の画像を生成する画像生成部と、を備え、
    前記動き補正部は、前記受信データに基づいて前記反射体に対する複数の画像フレームを取得し、
    前記複数の画像フレームを互いに比較し、前記反射体の動きを検出し、
    前記反射体の動きに基づいて前記トランスデューサー素子のチャネルごとのデータのディレイを調整することにより、前記動き補正を行い、
    前記合成データ取得部は、前記動き補正が行われたトランスデューサー素子のチャネルごとのデータに基づいて前記合成データを取得することを特徴とする超音波測定装置。
11. 前記合成フォーカシング地点それぞれに対する合成データをチャネルごとに合算して単一の受信データを生成し、前記単一の受信データに対して収差補正するデータ合算及び収差補正部をさらに備えることを特徴とする請求項１０に記載の超音波測定装置。

Images