(57)【要約】【課題】反射体の正確な画像生成ができる超音波測定方法及び超音波測定装置を提供する。【解決手段】本発明による超音波測定方法は、実際フォーカシング地点に対応する受信データに合成フォーカシング方法を適用し、反射体周辺領域の合成フォーカシング地点に対応する少なくとも１つの合成データを取得するステップと、前記少なくとも１つの合成データを用いて前記反射体の画像を生成するステップと、を有する。【選択図】      図３
2015150423.tif 000002
【特許請求の範囲】
【請求項１】
実際フォーカシング地点に対応する受信データに合成フォーカシング方法を適用し、反射体周辺領域の合成フォーカシング地点に対応する少なくとも１つの合成データを取得するステップと、
前記少なくとも１つの合成データを用いて前記反射体の画像を生成するステップと、を有することを特徴とする超音波測定方法。
【請求項２】
前記合成データを取得するステップ以前に、前記実際フォーカシング地点に超音波信号をビームフォーカシングするステップと、
前記超音波信号に対応して反射体が生成するアコースティックフィールドを受信して前記受信データを取得するステップと、をさらに有することを特徴とする請求項１に記載の超音波測定方法。
【請求項３】
前記合成データを取得するステップは、前記受信データに前記合成フォーカシング方法を用いて合成データセットを生成するステップを含み、
前記反射体の画像を生成するステップは、前記合成データセットを用いて前記反射体の画像を生成するステップを含むことを特徴とする請求項１に記載の超音波測定方法。
【請求項４】
前記合成データを取得するステップは、前記反射体の周辺領域に前記合成フォーカシング地点を決定するステップと、
前記合成フォーカシング地点と前記実際フォーカシング地点との間の第１ディレイを予測するステップと、を含むことを特徴とする請求項１に記載の超音波測定方法。
【請求項５】
前記合成データを取得するステップは、前記第１ディレイに基づいて少なくとも１つのトランスデューサー素子それぞれにおけるディレイを予測するステップと、
前記予測された前記少なくとも１つのトランスデューサー素子それぞれにおけるディレイに基づいて、前記受信データをトランスデューサー素子チャネルごとにディレイ補正するステップと、をさらに含むことを特徴とする請求項４に記載の超音波測定方法。
【請求項６】
前記合成データを取得するステップは、
前記反射体の周辺領域内の複数の合成フォーカシング地点を決定するステップと、
前記複数の合成フォーカシング地点それぞれに対してディレイ補正を行って、前記複数のフォーカシング地点に対応する合成データセットを取得するステップをさらに含むことを特徴とする請求項５に記載の超音波測定方法。
【請求項７】
前記反射体の画像を生成するステップは、前記合成フォーカシング地点それぞれに対する合成データをチャネルごとに合算して単一の受信データを生成し、単一の受信データを用いて前記反射体の画像を生成するステップを含むことを特徴とする請求項１に記載の超音波測定方法。
【請求項８】
前記反射体の画像を生成するステップは、前記単一の受信データに対して収差補正するステップをさらに含むことを特徴とする請求項７に記載の超音波測定方法。
【請求項９】
前記受信データに対して動き補正を行うステップをさらに有することを特徴とする請求項１に記載の超音波測定方法。
【請求項１０】
前記受信データに対して動き補正を行うステップは、前記受信データに基づいて前記反射体の複数の画像フレームを取得するステップと、
前記複数の画像フレームを相互に比較し、前記反射体の動きベクトルを検出するステップとを含むことを特徴とする請求項９に記載の超音波測定方法。
【請求項１１】
前記受信データに対して動き補正を行うステップは、前記検出された動きベクトルに基づいて、動き以前の反射体の位置と動き以後の反射体の位置との間の第２ディレイを予測するステップと、
前記予測された第２ディレイに基づいて、少なくとも１つのトランスデューサー素子それぞれにおけるディレイを予測するステップと、
前記予測された前記少なくとも１つのトランスデューサー素子それぞれにおけるディレイに基づいて、前記受信データをトランスデューサー素子チャネルごとにディレイ補正するステップと、をさらに含むことを特徴とする請求項１０に記載の超音波測定方法。
【請求項１２】
前記受信データをトランスデューサー素子チャネルごとにディレイ補正するステップは、最初の画像フレームでの反射体の位置から前記動きベクトルが検出された画像フレームに対応するサンプリング区間の受信データをトランスデューサー素子チャネルごとにディレイ補正することを特徴とする請求項１１に記載の超音波測定方法。
【請求項１３】
前記反射体の画像を生成するステップは、前記少なくとも１つの合成データ及び前記受信データを用いて前記反射体の画像を生成することを特徴とする請求項１に記載の超音波測定方法。
【請求項１４】
前記受信データに、フィルタリング又は復調の少なくとも１つを行うステップをさらに有することを特徴とする請求項１に記載の超音波測定方法。
【請求項１５】
前記合成フォーカシング方法又は前記動き補正の少なくとも１つを用いてディレイエラーを再算出するステップと、
前記合成フォーカシングを繰り返し処理するステップと、をさらに有することを特徴とする請求項９に記載の超音波測定方法。
【請求項１６】
実際フォーカシング地点に対応する受信データに合成フォーカシング方法を適用し、反射体周辺領域の合成フォーカシング地点に対応する少なくとも１つの合成データを取得する合成データ取得部と、
前記少なくとも１つの合成データを用いて、前記反射体の画像を生成する画像生成部と、を備えることを特徴とする超音波測定装置。
【請求項１７】
前記合成フォーカシング地点それぞれに対する合成データをチャネルごとに合算して単一の受信データを生成し、前記単一の受信データに対して収差補正するデータ合算及び収差補正部をさらに備えることを特徴とする請求項１６に記載の超音波測定装置。
【請求項１８】
前記画像生成部は、前記合成データを生成する以前に前記受信データに基づいて複数の画像フレームを生成し、
前記複数の画像フレームに基づいて前記反射体の動きを補正する動き補正部をさらに備えることを特徴とする請求項１６に記載の超音波測定装置。
【請求項１９】
前記複数の画像フレームを相互に比較し、前記反射体の動きベクトルを検出する動き推定部をさらに備えることを特徴とする請求項１８に記載の超音波測定装置。
【請求項２０】
実際フォーカシング地点に超音波信号をビームフォーカシングし、ビームフォーカシングに対応して反射体が生成するアコースティックフィールドを受信する少なくとも１つのトランスデューサー素子と、
前記少なくとも１つのトランスデューサー素子チャネルごとに受信データを取得し、前記受信データに合成フォーカシング方法を適用して前記反射体周辺領域の複数の合成フォーカシング地点それぞれに対する合成データを取得し、前記合成データに基づいて前記反射体の画像を取得するプロセッサと、を備えることを特徴とする超音波測定装置。
【請求項２１】
実際フォーカシング地点に超音波信号を送信し、ビームフォーカシングに対応して反射体が生成するアコースティックフィールドを受信する超音波信号送受信部と、
前記超音波信号送受信部で受信した受信データに基づいて複数の画像フレームを生成する画像生成部と、
前記生成された複数の画像フレームに基づいて前記反射体の動きを推定する動き推定部と、
前記推定された動きに基づいて前記受信データに対して動き補正を行う動き補正部と、を備えることを特徴とする超音波測定装置。
【請求項２２】
前記動き補正された受信データに対して合成フォーカシング方法を適用し、前記反射体周辺領域の複数の地点に対する合成データを生成する複数の合成フォーカシング生成部と、
前記合成データをチャネルごとに合算して単一の受信データを生成し、収差補正を行うデータ合算及び収差補正部と、をさらに備え、
前記画像生成部は、収差補正された単一の受信データを用いて前記反射体の画像を生成することを特徴とする請求項２１に記載の超音波測定装置。
【請求項２３】
実際フォーカシング地点に超音波信号をビームフォーカシングするようにトランスデューサー制御信号を出力するステップと、
前記ビームフォーカシングに対応して反射体が生成するアコースティックフィールドに対する受信データが入力されるステップと、
合成フォーカシング方法に基づいて、前記反射体周辺領域の複数の合成フォーカシング地点それぞれに対する合成データを取得するステップと、
前記合成データに基づいて前記反射体の画像を取得するステップと、を有することを特徴とする超音波測定方法。
【請求項２４】
超音波信号を実際ビームフォーカシングするようにトランスデューサー制御信号を出力するステップと、
ビームフォーカシングに対応して反射体が生成するアコースティックフィールドに対する受信データが入力されるステップと、
受信データからハーモニック成分を抽出するステップと、
前記ハーモニック成分を抽出した受信データに合成フォーカシング方法を適用して、前記反射体周辺領域の複数の合成フォーカシング地点それぞれに対する合成データを取得するステップと、
前記合成データに基づいて前記反射体の画像を取得するステップと、を有することを特徴とする超音波測定方法。
【請求項２５】
測定地点の周辺の実際フォーカシング地点を決定するステップと、
超音波トランスデューサーを用いて前記実際フォーカシング地点に実際ビームフォーカシングを行うステップと、
前記実際ビームフォーカシングに応答して超音波データを取得するステップと、
前記測定地点周辺の仮想フォーカシング地点を決定するステップと、
前記仮想フォーカシング地点に対応して、それぞれの超音波トランスデューサーに対して取得された超音波データのディレイを調整することによって合成データセットを算出するステップと、
それぞれの超音波トランスデューサーに対する合成データセットを合わせることによって、単一トランスデューサーデータセットを取得するステップと、
前記単一トランスデューサーデータセットに基づいて画像を生成するステップと、を有し、
前記測定地点は、前記仮想フォーカシング地点によって仮想的に取り囲まれることを特徴とする超音波測定方法。
【請求項２６】
前記実際フォーカシング地点と対応する仮想フォーカシング地点との間の第１ディレイを予測するステップをさらに有することを特徴とする請求項２５に記載の超音波測定方法。
【請求項２７】
前記第１ディレイに基づいて前記超音波トランスデューサーそれぞれに関する第２ディレイを予測するステップと、
前記超音波トランスデューサーそれぞれに対して、前記第２ディレイに基づいて取得された超音波データに対するディレイ補正を行うステップと、をさらに有することを特徴とする請求項２６に記載の超音波測定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波測定装置及び超音波測定方法に関し、特に、サンプル内の反射体に対する画像を測定する超音波測定装置及び超音波測定方法に関する。
【背景技術】
【０００２】
超音波測定装置は、サンプル内のフォーカシング地点に超音波信号を送信し、サンプル内の反射体によって生成されたアコースティックフィールドを受信し、受信されたアコースティックフィールドを用いてサンプル内の反射体画像を測定（生成）する装置である。
さらに、超音波測定装置は、リアルタイムでサンプルに対する画像を処理することができ、人体内部の組織検査のような医療用又は非破壊検査用などの様々な産業分野において利用されている。
【０００３】
超音波測定装置は、複数のトランスデューサー素子を含み得、複数のトランスデューサー素子の送信時期を調整してサンプル内の特定の地点にビームフォーカシングすることができる。
超音波測定装置は、ビームフォーカシングを行って超音波信号を送信し、反射体によって反射又は散乱されたアコースティックフィールドを用いて反射体に対する画像を測定（生成）する。
【０００４】
しかし、サンプル内の測定しようとする反射体以外の他の反射体がさらに含まれる場合には、他の反射体からのアコースティックフィールドが測定しようとする反射体からのアコースティックフィールドと干渉することがある。
そのため、測定しようとする反射体を正確に画像測定（生成）することが困難であるという問題があった。
【発明の概要】
【発明が解決しようとする課題】
【０００５】
本発明は上記従来の超音波測定装置における問題点に鑑みてなされたものであって、本発明の目的は、反射体の正確な画像生成ができる超音波測定方法及び超音波測定装置を提供することにある。
【課題を解決するための手段】
【０００６】
上記目的を達成するためになされた本発明による超音波測定方法は、実際フォーカシング地点に対応する受信データに合成フォーカシング方法を適用し、反射体周辺領域の合成フォーカシング地点に対応する少なくとも１つの合成データを取得するステップと、前記少なくとも１つの合成データを用いて前記反射体の画像を生成するステップと、を有することを特徴とする。
【０００７】
前記合成データを取得するステップ以前に、前記実際フォーカシング地点に超音波信号をビームフォーカシングするステップと、前記超音波信号に対応して反射体が生成するアコースティックフィールドを受信して前記受信データを取得するステップと、をさらに有することが好ましい。
前記合成データを取得するステップは、前記受信データに前記合成フォーカシング方法を用いて合成データセットを生成するステップを含み、前記反射体の画像を生成するステップは、前記合成データセットを用いて前記反射体の画像を生成するステップを含むことが好ましい。
前記合成データを取得するステップは、前記反射体の周辺領域に前記合成フォーカシング地点を決定するステップと、前記合成フォーカシング地点と前記実際フォーカシング地点との間の第１ディレイを予測するステップと、を含むことが好ましい。
前記合成データを取得するステップは、前記第１ディレイに基づいて少なくとも１つのトランスデューサー素子それぞれにおけるディレイを予測するステップと、前記予測された前記少なくとも１つのトランスデューサー素子それぞれにおけるディレイに基づいて、前記受信データをトランスデューサー素子チャネルごとにディレイ補正するステップと、をさらに含むことが好ましい。
前記合成データを取得するステップは、前記反射体の周辺領域内の複数の合成フォーカシング地点を決定するステップと、前記複数の合成フォーカシング地点それぞれに対してディレイ補正を行って、前記複数のフォーカシング地点に対応する合成データセットを取得するステップをさらに含むことが好ましい。
前記反射体の画像を生成するステップは、前記合成フォーカシング地点それぞれに対する合成データをチャネルごとに合算して単一の受信データを生成し、単一の受信データを用いて前記反射体の画像を生成するステップを含むことが好ましい。
前記反射体の画像を生成するステップは、前記単一の受信データに対して収差補正するステップをさらに含むことが好ましい。
【０００８】
前記受信データに対して動き補正を行うステップをさらに有することが好ましい。
前記受信データに対して動き補正を行うステップは、前記受信データに基づいて前記反射体の複数の画像フレームを取得するステップと、前記複数の画像フレームを相互に比較し、前記反射体の動きベクトルを検出するステップとを含むことが好ましい。
前記受信データに対して動き補正を行うステップは、前記検出された動きベクトルに基づいて、動き以前の反射体の位置と動き以後の反射体の位置との間の第２ディレイを予測するステップと、前記予測された第２ディレイに基づいて、少なくとも１つのトランスデューサー素子それぞれにおけるディレイを予測するステップと、前記予測された前記少なくとも１つのトランスデューサー素子それぞれにおけるディレイに基づいて、前記受信データをトランスデューサー素子チャネルごとにディレイ補正するステップと、をさらに含むことが好ましい。
前記受信データをトランスデューサー素子チャネルごとにディレイ補正するステップは、最初の画像フレームでの反射体の位置から前記動きベクトルが検出された画像フレームに対応するサンプリング区間の受信データをトランスデューサー素子チャネルごとにディレイ補正することが好ましい。
前記反射体の画像を生成するステップは、前記少なくとも１つの合成データ及び前記受信データを用いて前記反射体の画像を生成することが好ましい。
前記受信データに、フィルタリング又は復調の少なくとも１つを行うステップをさらに有することが好ましい。
前記合成フォーカシング方法又は前記動き補正の少なくとも１つを用いてディレイエラーを再算出するステップと、前記合成フォーカシングを繰り返し処理するステップと、をさらに有することが好ましい。
【０００９】
上記目的を達成するためになされた本発明による超音波測定装置は、実際フォーカシング地点に対応する受信データに合成フォーカシング方法を適用し、反射体周辺領域の合成フォーカシング地点に対応する少なくとも１つの合成データを取得する合成データ取得部と、前記少なくとも１つの合成データを用いて、前記反射体の画像を生成する画像生成部と、を備えることを特徴とする。
【００１０】
前記合成フォーカシング地点それぞれに対する合成データをチャネルごとに合算して単一の受信データを生成し、前記単一の受信データに対して収差補正するデータ合算及び収差補正部をさらに備えることが好ましい。
前記画像生成部は、前記合成データを生成する以前に前記受信データに基づいて複数の画像フレームを生成し、前記複数の画像フレームに基づいて前記反射体の動きを補正する動き補正部をさらに備えることが好ましい。
前記複数の画像フレームを相互に比較し、前記反射体の動きベクトルを検出する動き推定部をさらに備えることが好ましい。
【００１１】
また、上記目的を達成するためになされた本発明による超音波測定装置は、実際フォーカシング地点に超音波信号をビームフォーカシングし、ビームフォーカシングに対応して反射体が生成するアコースティックフィールドを受信する少なくとも１つのトランスデューサー素子と、前記少なくとも１つのトランスデューサー素子チャネルごとに受信データを取得し、前記受信データに合成フォーカシング方法を適用して前記反射体周辺領域の複数の合成フォーカシング地点それぞれに対する合成データを取得し、前記合成データに基づいて前記反射体の画像を取得するプロセッサと、を備えることを特徴とする。
【００１２】
また、上記目的を達成するためになされた本発明による超音波測定装置は、実際フォーカシング地点に超音波信号を送信し、ビームフォーカシングに対応して反射体が生成するアコースティックフィールドを受信する超音波信号送受信部と、前記超音波信号送受信部で受信した受信データに基づいて複数の画像フレームを生成する画像生成部と、前記生成された複数の画像フレームに基づいて前記反射体の動きを推定する動き推定部と、前記推定された動きに基づいて前記受信データに対して動き補正を行う動き補正部と、を備えることを特徴とする。
【００１３】
前記動き補正された受信データに対して合成フォーカシング方法を適用し、前記反射体周辺領域の複数の地点に対する合成データを生成する複数の合成フォーカシング生成部と、前記合成データをチャネルごとに合算して単一の受信データを生成し、収差補正を行うデータ合算及び収差補正部と、をさらに備え、前記画像生成部は、収差補正された単一の受信データを用いて前記反射体の画像を生成することが好ましい。
【００１４】
また、上記目的を達成するためになされた本発明による超音波測定方法は、実際フォーカシング地点に超音波信号をビームフォーカシングするようにトランスデューサー制御信号を出力するステップと、前記ビームフォーカシングに対応して反射体が生成するアコースティックフィールドに対する受信データが入力されるステップと、合成フォーカシング方法に基づいて、前記反射体周辺領域の複数の合成フォーカシング地点それぞれに対する合成データを取得するステップと、前記合成データに基づいて前記反射体の画像を取得するステップと、を有することを特徴とする。
【００１５】
また、上記目的を達成するためになされた本発明による超音波測定方法は、超音波信号を実際ビームフォーカシングするようにトランスデューサー制御信号を出力するステップと、ビームフォーカシングに対応して反射体が生成するアコースティックフィールドに対する受信データが入力されるステップと、受信データからハーモニック成分を抽出するステップと、前記ハーモニック成分を抽出した受信データに合成フォーカシング方法を適用して、前記反射体周辺領域の複数の合成フォーカシング地点それぞれに対する合成データを取得するステップと、前記合成データに基づいて前記反射体の画像を取得するステップと、を有することを特徴とする。
【００１６】
また、上記目的を達成するためになされた本発明による超音波測定方法は、測定地点の周辺の実際フォーカシング地点を決定するステップと、超音波トランスデューサーを用いて前記実際フォーカシング地点に実際ビームフォーカシングを行うステップと、前記実際ビームフォーカシングに応答して超音波データを取得するステップと、前記測定地点周辺の仮想フォーカシング地点を決定するステップと、前記仮想フォーカシング地点に対応して、それぞれの超音波トランスデューサーに対して取得された超音波データのディレイを調整することによって合成データセットを算出するステップと、それぞれの超音波トランスデューサーに対する合成データセットを合わせることによって、単一トランスデューサーデータセットを取得するステップと、前記単一トランスデューサーデータセットに基づいて画像を生成するステップと、を有し、前記測定地点は、前記仮想フォーカシング地点によって仮想的に取り囲まれることを特徴とする。
【００１７】
前記実際フォーカシング地点と対応する仮想フォーカシング地点との間の第１ディレイを予測するステップをさらに有することが好ましい。
前記第１ディレイに基づいて前記超音波トランスデューサーそれぞれに関する第２ディレイを予測するステップと、前記超音波トランスデューサーそれぞれに対して、前記第２ディレイに基づいて取得された超音波データに対するディレイ補正を行うステップと、をさらに有することが好ましい。
【発明の効果】
【００１８】
本発明に係る超音波測定装置及び超音波測定方法によれば、相対的に少ない実際のフォーカシング地点についてのみビームフォーカシングを行っても、合成フォーカシング方法を用いて測定地点周辺領域の相対的に多くのフォーカシング地点に対する受信データを取得することができる。そのため、演算時間及び測定時間を著しく縮小することができる、という効果がある。
【図面の簡単な説明】
【００１９】
【図１Ａ】本発明の実施形態に係る超音波測定装置の概略を示すブロック図である。
【図１Ｂ】本発明の実施形態に係る超音波測定装置の概略を示すブロック図である。
【図２Ａ】サンプル内に単一反射体が存在する場合を説明するための概念図である。
【図２Ｂ】サンプル内の複数反射体が含まれた場合のアコースティックフィールドの概念図である。
【図３】本発明の実施形態に係る超音波測定装置で、比較例による観測対象である第１反射体以外の他の反射体が含まれたサンプルに対する超音波測定方法を説明するためのフローチャートである。
【図４Ａ】図３に示す超音波測定方法を説明するための概念図である。
【図４Ｂ】図３に示す超音波測定方法を説明するための概念図である。
【図４Ｃ】受信データ及び受信データのチャネルごとの合算を説明するための図である。
【図４Ｄ】受信データ及び受信データのチャネルごとの合算を説明するための図である。
【図４Ｅ】受信データのチャネルごとの合算による他反射体の影響除去を説明するための概念図である。
【図４Ｆ】受信データのチャネルごとの合算による他反射体の影響除去を説明するための概念図である。
【図５】本発明による超音波測定装置の他の実施形態に係る超音波測定方法を説明するためのフローチャートである。
【図６Ａ】様々な実施形態に係るビームフォーカシング方法を説明するための概念図である。
【図６Ｂ】様々な実施形態に係るビームフォーカシング方法を説明するための概念図である。
【図７Ａ】フォーカシング地点にビームフォーカシングを行った構成を説明するための概念図である。
【図７Ｂ】フォーカシング地点にビームフォーカシングを行った構成を説明するための概念図である。
【図７Ｃ】フォーカシング地点にビームフォーカシングを行った構成を説明するための概念図である。
【図８】本発明の一実施形態に係るビームフォーカシング方法を説明するためのフローチャートである。
【図９】本発明の複数フォーカシング地点に対するビームフォーカシング方法を説明するためのフローチャートである。
【図１０】本発明の一実施形態に係る実際フォーカシング受信データを説明するための図である。
【図１１Ａ】サンプル内のフォーカシング地点を説明するための概念図である。
【図１１Ｂ】サンプル内のフォーカシング地点を説明するための概念図である。
【図１２】本発明の一実施形態に係る合成データ取得方法を説明するためのフローチャートである。
【図１３】本発明の他の実施形態に係る複数の合成フォーカシング地点に対する合成データ取得方法を説明するためのフローチャートである。
【図１４】本発明の一実施形態に係る合成データの例示である。
【図１５】本発明の他の実施形態に係る受信データ測定時のサンプル内反射体が移動する場合、反射体動き補正を行う超音波測定方法を説明するためのフローチャートである。
【図１６】本発明の他の実施形態に係る動き補正を行う超音波測定方法を説明するためのフローチャートである。
【図１７Ａ】サンプル内の反射体の動きを説明するための概念図である。
【図１７Ｂ】本発明の実施形態に係る複数の画像フレームの概念図である。
【図１７Ｃ】動きベクトル検出及びディレイ補正工程を説明するための概念図である。
【図１８】チャネルごとにディレイを補正する工程を説明するための図である。
【図１９】本発明の一実施形態に係る動き補正を行う超音波測定装置の概略構成を示すブロック図である。
【図２０Ａ】収差補正工程を説明するための概念図である。
【図２０Ｂ】収差補正工程を説明するための概念図である。
【図２０Ｃ】収差補正部の相互相関に基づく収差補正を説明するための概念図である。
【図２０Ｄ】相互相関に基づく収差補正以前の画像写真である。
【図２０Ｅ】相互相関に基づく収差補正以後の画像写真である。
【図２１】本発明のさらに他の実施形態に係る超音波測定方法を説明するためのフローチャートである。
【発明を実施するための形態】
【００２０】
次に、本発明に係る超音波測定装置及び超音波測定方法を実施するための形態の具体例を図面を参照しながら説明する。
【００２１】
図１Ａ及び１Ｂは、本発明の実施形態に係る超音波測定装置の概略を示すブロック図である。
図１Ａに示すように、超音波測定装置１８は、プロセッサ３０、表示部２０、媒質コンテナ４０、送受信部５０、増幅器及びフィルタ７０、及びアナログ−デジタル変換部８０を備える。
送受信部５０は、少なくとも１つのトランスデューサー素子（５１〜５４）を含む。
【００２２】
プロセッサ３０は、超音波測定装置全般の動作を制御する。
プロセッサ３０は、格納部（図示せず）に格納されたプログラム又はアルゴリズムに基づいて超音波測定装置全般の動作を制御する。
プロセッサ３０は、送受信部５０の超音波信号の送信時期及び送信地点を決定する。
詳細には、プロセッサ３０は、超音波信号のビームフォーカシング地点を決定する。
【００２３】
プロセッサ３０は、決定されたビームフォーカシング地点に対応して少なくとも１つのトランスデューサー素子（５１〜５４）の各駆動時期を決定する。
プロセッサ３０は、少なくとも１つのトランスデューサー素子（５１〜５４）の駆動時期に対する電気的な信号を生成し、生成された駆動時期に対する電気的な信号を少なくとも１つのトランスデューサー素子（５１〜５４）のそれぞれに出力する。
例えば、プロセッサ３０は、超音波信号がサンプル６０上の少なくとも１つの地点にビームフォーカスされるよう、少なくとも１つのトランスデューサー素子（５１〜５４）のそれぞれの送信時期を決定する。
【００２４】
少なくとも１つのトランスデューサー素子（５１〜５４）のそれぞれは、電気的な信号を超音波信号に変換する。又は、少なくとも１つのトランスデューサー素子（５１〜５４）それぞれは、受信されたアコースティックフィールドを電気的な信号に逆変換し得る。
少なくとも１つのトランスデューサー素子（５１〜５４）のそれぞれは、プロセッサ３０から受信された電気的な信号に基づいて各駆動時期に超音波信号を送信する。少なくとも１つのトランスデューサー素子（５１〜５４）からの超音波信号は互いに干渉してもよく、特定の地点にビームフォーカスされてもよい。
【００２５】
少なくとも１つのトランスデューサー素子（５１〜５４）は、受信されたアコースティックフィールドを電気的な信号に変換する。一方、ここで、アコースティックフィールドは、サンプル６０内の反射体によって反射又は散乱された媒質上のフィールドである。
変換された電気的な信号は、増幅器及びフィルタ７０によって予め設定された利得に増幅されてフィルタリングされる。
さらに、フィルタリングされた電気的な信号は、アナログ−デジタル変換器８０によってデジタル信号に変換される。本実施形態では、受信されたアコースティックフィールドが変換された電気的な信号は受信データと命名するか、又は、変換されたデジタル信号を受信データとして命名してもよい。
受信データは、少なくとも１つのトランスデューサー素子（５１〜５４）のチャネルごとのデータを含み、これに関連しては後で、表１を参照して詳細に説明する。
【００２６】
少なくとも１つのトランスデューサー素子（５１〜５４）は、一定期間の間にアコースティックフィールドを受信する。
時間の流れに伴って、少なくとも１つのトランスデューサー素子（５１〜５４）のそれぞれに受信されるアコースティックフィールドも変更される。
受信データも一定期間の間の時間によるデータであってもよい。例えば、受信データは、アナログ−デジタル変換器８０のサンプリング期間による時系列的なデータであってもよく、これに関連しては後で、表１を参照して詳細に説明する。
【００２７】
プロセッサ３０は、受信データが入力され、入力された受信データに基づいてサンプル６０に対する画像を生成する。
プロセッサ３０は、生成された画像を表示部２０が表示するように制御してもよい。
一方、プロセッサ３０は、サンプル６０の少なくとも１つの実際フォーカシング地点に超音波信号をビームフォーカシングするように少なくとも１つのトランスデューサー素子（５１〜５４）の送信時期を決定する。
さらに、プロセッサ３０は、少なくとも１つのトランスデューサー素子（５１〜５４）のチャネルごとに受信データを取得する。
【００２８】
プロセッサ３０は、合成フォーカシング方法に基づいて、サンプル６０内の反射体周辺領域の複数の合成フォーカシング地点それぞれに対する合成データを取得する。
ここで、複数の合成データを合成データセットと命名する。
プロセッサ３０は、合成データに基づいて反射体の画像を取得する。上述した工程については詳細に後述する。
媒質コンテナ４０は、超音波信号が移動するための媒質を入れるコンテナであってもよい。
一方、媒質が準備された測定環境として、媒質コンテナ４０は超音波測定装置に含まれないこともある。
【００２９】
図１Ｂは、本発明の他の実施形態に係る超音波測定装置の概略を示すブロック図である。
図１Ｂの実施形態に係る超音波測定装置は、図１Ａの実施形態に係る超音波測定装置とは対照的に受信部９０をさらに備える。
受信部９０は、例えば、少なくとも１つのチャネルを有するハイドロフォン（ｈｙｄｒｏｐｈｏｎｅ）又は追加トランスデューサー素子などで実現され得る。
【００３０】
図２Ａは、サンプル内に単一反射体が存在する場合を説明するための概念図である。
図２Ａに示していないが、少なくとも１つのトランスデューサー素子（ｃｈ１〜ｃｈ５）は、サンプル１０に超音波信号を送信する。
サンプル１０は超音波信号を受信し、これに対応してアコースティックフィールドを生成する。
一方、図２Ａに示すように、サンプル１０内には単一の反射体２０１が含まれてもよい。
反射体２０１は、超音波信号に対応したり、又はビームフォーカシング地点からの超音波信号に対応してアコースティックフィールド２１０を放出する。
【００３１】
少なくとも１つのトランスデューサー素子（ｃｈ１〜ｃｈ５）は、アコースティックフィールド２１０を受信し、これを電気的な信号に変換する。
超音波測定装置は、変換された信号に基づいて反射体２０１の位置を測定し得る。
図２Ａに示すように、単一の反射体２０１によって生成されたアコースティックフィールド２１０は、実質的に球面波の形態に形成される。これはホイヘンスの原理によって反射体２０１を中心に同心円形態のアコースティックフィールドが形成されることから起因する。
単一の反射体２０１が実質的な球面波の形態であるアコースティックフィールド２１０を放出することによって、超音波測定装置は球面波の収差（ａｂｅｒｒａｔｉｏｎ）補正を行って単一の反射体２０１に対する画面を生成する。超音波測定装置が収差補正する構成については後述する。
【００３２】
一方、図２Ｂは、サンプル内の複数反射体が含まれた場合のアコースティックフィールドの概念図である。
図２Ｂに示すように反射体（２０１、２８０、２９０）がサンプル１０内に存在する。
少なくとも１つのトランスデューサー素子（ｃｈ１〜ｃｈ５）は、サンプル１０に超音波信号を送信する。
反射体２０１は、超音波信号を受信してアコースティックフィールド２１０を生成して放出する。反射体２８０は、超音波信号を受信してアコースティックフィールド２２０を生成して放出する。反射体２９０は、超音波信号を受信してアコースティックフィールド２３０を生成して放出する。
複数のアコースティックフィールド（２１０、２２０、２３０）は媒質内で互いに干渉する。
【００３３】
図２Ｂに示すように、複数の反射体（２０１、２８０、２９０）それぞれによって生成されたアコースティックフィールド（２１０、２２０、２３０）は、実質的に球面波の形態に形成される。
ただし、球面波形態のアコースティックフィールド（２１０、２２０、２３０）は互いに干渉され得る。そのため、複数の球面波が干渉された形態のアコースティックフィールドが生成され、少なくとも１つのトランスデューサー素子（ｃｈ１〜ｃｈ５）は、干渉された形態のアコースティックフィールドを受信して電気的な信号に変換する。
トランスデューサー素子（ｃｈ１〜ｃｈ５）に球面波の形態と相違の大きいアコースティックフィールドが入力されることによって、超音波測定装置は、特定の反射体２０１に対する正確な画像を生成することができない。
【００３４】
図３は、本発明の実施形態に係る超音波測定装置で、比較例による観測対象である第１反射体以外の他の反射体が含まれたサンプルに対する超音波測定方法を説明するためのフローチャートである。
図３に示す比較例は、図４Ａ〜４Ｃを参照して詳細に説明することにする。
【００３５】
まず、図３を参照すると、ステップＳ３１０において、超音波測定装置は測定地点を決定する。
ここで、測定地点は、画像生成を行うための地点である。例えば、図４Ａに示すように、サンプル１０内には複数の反射体（４０１、４８０、４９０）が含まれる。
超音波測定方法は、複数の反射体（４０１、４８０、４９０）のうちから第１反射体４０１を測定地点に決定するとする。超音波測定方法は、後述する方法に基づいて第１反射体４０１以外の残りの反射体（４８０、４９０）の測定時に及ぼす影響を抑制したり最小化することができる。
【００３６】
再び図３を参照すると、ステップＳ３２０において、超音波測定方法は、測定地点周辺領域のｎ地点にビームフォーカシングを行う。
超音波測定装置は、複数のトランスデューサー素子を含み得る。超音波測定装置は、複数のトランスデューサー素子それぞれの超音波信号の送信時期を調整する。
超音波測定装置は、複数のトランスデューサー素子それぞれの送信時期を含むトランスデューサー制御信号を出力する。
そのため、超音波測定装置は、サンプル１０内の特定地点に超音波信号を集中するようにするビームフォーカシングを行う。
【００３７】
例えば、図４Ｂに示す第１地点４１１にビームフォーカシングを行う場合に、超音波測定装置は第２トランスデューサー素子ｃｈ２の超音波信号の送信時期を相対的に遅くに調整して、第５トランスデューサー素子ｃｈ５の超音波信号の送信時期を相対的に早く調整する。
例えば、第１地点４１１から第２トランスデューサー素子ｃｈ２までの距離が最短であり、第１地点４１１から第５トランスデューサー素子ｃｈ５までの距離が最長であると仮定する。
第２トランスデューサー素子ｃｈ２からの超音波信号が相対的に遅れた時期に送信され、第５トランスデューサー素子ｃｈ５からの超音波信号が相対的に早い時期に送信されることで、複数の超音波信号が第１地点４１１に同時に集中する。超音波測定装置は、残りのトランスデューサー素子の超音波信号の送信時期を決定して全てのトランスデューサー素子からの超音波信号が第１地点４１１に同時に集中するように制御する。
【００３８】
他の実施形態に係る超音波測定方法として、超音波信号が第２地点（図示せず）から発生するようにビームフォーカシングを行ってもよい。
例えば、超音波測定装置は、第２地点（図示せず）をトランスデューサーの後方の一地点として決定してもよく、第２地点（図示せず）から超音波が拡散されるようにビームフォーカシングを行ってもよい。これについて図７Ｂを参照して詳細に説明する。
それと共に、スキャットラインごとに深度方向にフォーカスされる位置が互いに異なってもよい。
【００３９】
第１地点４１１にビームフォーカシングを行った後、第１反射体４０１によってアコースティックフィールドが生成される。
ステップＳ３３０において、少なくとも１つのトランスデューサー素子（ｃｈ１〜ｃｈ５）はアコースティックフィールドを受信し、受信されたアコースティックフィールドを電気的な信号に変換する。
超音波測定装置は、受信されたアコースティックフィールドに対する受信データを格納する。超音波測定装置は、受信データをチャネルごとにｃｈ１〜ｃｈ５に格納してもよい。
超音波測定装置は、受信データを一定期間の間に格納してもよい。例えば、表１は、第１地点４１１に対応した格納された受信データの例示である。
【００４０】
【表１】
2015150423.tif 000003
【００４１】
表１の受信データは、予め設定された期間でサンプリングされて格納されたものとして、「ｔ」が指示する１〜６はサンプリングナンバーである。
格納された受信データは、チャネルごと（ｃｈ１〜ｃｈ５）に区分されて格納される。
【００４２】
ステップＳ３４０において、超音波測定方法は「ｎ」が「Ｎ」を超過するか否かを判断する。
例えば、「Ｎ」が“２５”と設定されてもよい。「ｎ」が「Ｎ」を超過しない場合、超音波測定方法はステップＳ３５０において「ｎ」の値を“１”増加させ、これに対応してステップＳ３６０においてフォーカシング地点を移動させる。
超音波測定装置は、第２地点４１２に対してビームフォーカシングを行う。超音波測定装置は、第２地点４１２に対応して少なくとも１つのトランスデューサー素子（ｃｈ１〜ｃｈ５）の送信時期を再調整し、第２地点４１２に対しビームフォーカシングを行う。
【００４３】
超音波測定装置は、測定地点４０１の周辺地域にフォーカシング地点（４１１〜４３５）を決定してもよく、例えば、図４Ｂに示すように、フォーカシング地点（４１１〜４３５）がそれぞれ予め設定された間隔に離隔して全体的に長方形の格子形態を有するようビームフォーカシング地点を決定してもよい。
すなわち、一部のビームフォーカシング地点は、第１反射体４０１よりトランスデューサーに近く配置し、一部のビームフォーカシング地点は第１反射体４０１よりトランスデューサーから遠く配置し、又は、一部のビームフォーカシング地点は。第１反射体４０１とトランスデューサーから同じ距離に配置する。
一方、図４Ｂでは、２次元の地点に２５個のフォーカシング地点に対してビームフォーカシングするように図に示しているが、これは一例に過ぎず、他の実施形態に係る超音波測定装置は、３次元の例を挙げて１２５個フォーカシング地点に対してビームフォーカシングを行ってもよい。
【００４４】
超音波測定装置は、決定された全てのフォーカシング地点に対してビームフォーカシングを行い、各フォーカシング地点に対してチャネルごとのデータを格納する。
例えば、超音波測定装置は、図４Ｃに示すような受信データを格納してもよい。
受信データは、フォーカシング地点ごと（ｎ＝１、２、．．．、２５）に格納される。
各フォーカシング地点ごとの受信データは、チャネルごとに格納される。
第１フォーカシング地点４１１に対する受信データは、第１チャネル〜第５チャネル（４４１〜４４５）ごとに格納される。
さらに、第１チャネル〜第５チャネル（４４１〜４４５）ごとのデータそれぞれは、表１のように予め設定された時間の間に格納される。第２フォーカシング地点４１１に対する受信データも第１チャネル〜第５チャネル（４５１〜４５５）ごとに格納され、第２５フォーカシング地点４３５に対する受信データも第１チャネル〜第５チャネル（４６１〜４６５）ごとに格納される。
【００４５】
全てのフォーカシング地点（４１１〜４３５）に対して測定が完了すれば（ｎ＞Ｎ）、ステップＳ３７０において超音波測定装置はチャネルごとに受信データを合算する。
例えば、図４Ｄに示すように、超音波測定装置は、チャネルごとに受信データを合算する。
図４Ｄに示すように、超音波測定装置は、複数のフォーカシング地点（４１１〜４３５）に対する受信データの第１チャネルデータ（４４１、４５１、．．．、４６１）を合算し、合算された第１チャネルデータ４７１を算出する。
【００４６】
さらに、超音波測定装置は、複数のフォーカシング地点に対する受信データをチャネルごとに合算してチャネルごとの合算データ（４７１〜４７５）を生成する。
そのため、超音波測定装置は、単一の受信データを生成する。
ここで、単一の受信データとは、複数フォーカシング地点に対する受信データがチャネルごとに合算されて１つの単一の受信データが生成されることを示す。
チャネルごとのデータ合算を行うために、観測地点（第１反射体４０１）以外の他の反射体（４８０、４９０）によるアコースティックフィールドの影響が相殺され得る。
これについて、図４Ｅ及び４Ｆを参照して説明することにする。
【００４７】
ステップＳ３８０において、超音波測定装置は収差補正をさらに行う。
ステップＳ３９０において、超音波測定装置は、チャネルごとの合算データに基づいて画像を生成する。
【００４８】
図４Ｅ及び４Ｆは、チャネルごとのデータ合算によって他の反射体による影響が相殺される構成を説明するための概念図である。
図４Ｅは、フォーカシング地点が第１地点４１１である場合に対する概念図である。
第１反射体４０１及び第２反射体４８０のそれぞれからトランスデューサー素子に受信される時間はそれぞれ「ｔ１」及び「ｔ２」である。
さらに、図４Ｆは、フォーカシング地点が第２地点４１２である場合に対する概念図である。
第１反射体４０１及び第２反射体４８０のそれぞれからトランスデューサー素子に受信される時間はそれぞれ「ｔ３」及び「ｔ４」である。ここで、“ｔ１−ｔ２”と“ｔ３−ｔ４”は互いに異なるが、チャネルごとのデータ合算時の観測地点ではない他の反射体からの信号はディレイ差で相殺される。
一方、図３に示す比較例による超音波測定方法では、複数のフォーカシング地点に対してビームフォーカシングするための測定時間が増加する。
【００４９】
図５は、本発明による超音波測定装置の他の実施形態に係る超音波測定方法を説明するためのフローチャートである。
ステップＳ５１０において、超音波測定装置は、少なくとも１つの実際フォーカシング地点に超音波信号をビームフォーカシングする。
特に、本実施形態に係る超音波測定方法では、トランスデューサー素子の個数以下の個数の実際フォーカシング地点にビームフォーカシングを行う。
図６Ａ及び６Ｂは、様々な実施形態に係るビームフォーカシング方法を説明するための概念図である。
図６Ａ及び６Ｂの実施形態で、サンプル１０内の第１地点６０１が測定地点に決定されたと仮定する。
【００５０】
本実施形態に係る超音波測定方法は、第１地点６０１の周辺領域に複数のフォーカシング地点（６１１〜６１５）を決定する。
超音波測定装置は、フロントフォーカシング方法に基づいて、トランスデューサー素子（ｃｈ１〜ｃｈ５）の前面に複数のフォーカシング地点（６１１〜６１５）を決定する。
超音波測定装置は、トランスデューサー素子の個数である“５つ”より以下の個数でフォーカシング地点（６１１〜６１５）の個数を決定する。
他の実施形態に係る超音波測定方法では、リアフォーカシング方法に基づいてトランスデューサー素子（ｃｈ１〜ｃｈ５）の後面に複数のフォーカシング地点（６２１〜６２５）を決定する。
【００５１】
図７Ａ〜７Ｃは、フォーカシング地点にビームフォーカシングを行った構成を説明するための概念図である。
図７Ａは、スキャンライン（ｓｃａｎ  ｌｉｎｅ）ごとのフォーカシング方法に基づいたビームフォーカシング実行結果であり、図７Ｂはリアフォーカシング方法に基づいたビームフォーカシング実行結果であり、図７Ｃはフロントフォーカシング方法に基づいたビームフォーカシング実行結果である。
【００５２】
図８は、本発明の実施形態に係るビームフォーカシング方法を説明するためのフローチャートである。
ステップＳ８０５において、超音波測定装置は、例えば、図６Ａ及び図６Ｂに示す測定地点６０１のように測定地点を決定する。
ステップＳ８１０において、超音波測定装置は、少なくとも１つの実際フォーカシング地点を決定する。
ステップＳ８２０において、超音波測定装置は、少なくとも１つのトランスデューサーの送信時期を決定、すなわち、トランスデューサー素子間の駆動ディレイを決定する。
ステップＳ８３０において、超音波測定装置は、少なくとも１つのトランスデューサー素子を決定された送信時期に超音波信号を送信するように制御し、そのため、特定実際フォーカシング地点にビームフォーカシングを行う。
【００５３】
図９は、本発明の実施形態による複数フォーカシング地点に対するビームフォーカシング方法を説明するためのフローチャートである。
ステップＳ９１０において、超音波測定装置は測定地点を決定する。
ステップＳ９２０において、超音波測定装置は、少なくとも１つの実際フォーカシング地点を決定する。例えば、超音波測定装置は、図６Ａ又は６Ｂに示すように、フロントフォーカシング方式又はリアフォーカシング方式に基づいて少なくとも１つの実際フォーカシング地点を決定する。
ステップＳ９３０において、超音波測定方法は、「ｋ」地点にビームフォーカシングする。
これに対応して、ステップＳ９４０において、アコースティックフィールドを受信して格納する。
「ｋ」は「Ｋ」値になるまでフォーカシング地点を移動させながらビームフォーカシングを行い、「ｋ」が「Ｋ」値以下であれば（ステップＳ９５０）、超音波測定装置はステップＳ９６０で「ｋ」値を１ずつ増加させ、ステップＳ９７０でフォーカシング地点を移動させる。
【００５４】
再び図５を参照すると、ステップＳ５２０において、超音波測定装置は、フォーカシング地点それぞれに対応するアコースティックフィールドを受信して格納する。
ステップＳ５３０において、超音波測定装置は、受信データそれぞれのディレイを調整して合成フォーカシングを行う。
合成フォーカシングは、サンプル上の第１実際フォーカシング地点に対応した受信データを用いて、第２合成フォーカシング地点に対応する合成データを生成する方法である。
詳細には、合成フォーカシングは、第１実際フォーカシング地点にビームフォーカシングを行った後、受信されたチャネルごとのデータそれぞれのディレイを調整し、第２合成フォーカシング地点にビームフォーカシングを行い、これに対応して受信されたデータを生成する方法である。
【００５５】
例えば、超音波測定装置は、第１フォーカシング地点に実際ビームフォーカシングを実行して、これに対応する第１受信データを取得する。ここで、第１実際フォーカシング地点は複数であってもよく、第１受信データも複数であってもよい。
超音波測定装置は、実際ビームフォーカシングを行った第１実際フォーカシング地点と合成フォーカシングを行う第２合成フォーカシング地点との間の距離、及びディレイのうち少なくとも１つを判断する。
さらに、超音波測定装置は、第１実際フォーカシング地点と第２合成フォーカシング地点との間のディレイに基づいて、トランスデューサー素子チャネルごとのディレイを予測する。
【００５６】
超音波測定装置は、予測されたトランスデューサー素子チャネルごとのディレイに基づいて、第１受信データのチャネルごとのデータをディレイ調整する。
そのため、超音波測定装置は、第２合成フォーカシング地点にビームフォーカシングを行った場合に受信される受信データの第２受信データを取得する。
すなわち、超音波測定装置は、第２合成フォーカシング地点に実際にビームフォーカシングすることなく、第２合成フォーカシング地点に対応した受信データを取得する。
【００５７】
ステップＳ５４０において、超音波測定装置は、測定地点周辺領域の複数の合成フォーカシング地点に対する合成データを取得する。
合成データを取得する方法の詳細は後述する。
ステップＳ５５０において、超音波測定装置は、複数の合成データのチャネルごとのデータを合算する。
チャネルごとのデータ合算にしたがって、サンプル内に複数の反射体がある場合にも、測定地点以外の他の反射体の影響が抑制される。
他の実施形態で、超音波測定装置は、合成データ及び受信データに対してチャネルごとのデータを合算してもよい。
【００５８】
ステップＳ５６０において、超音波測定装置は、チャネルごとの合算データに基づいて測定地点に対する画像を生成する。
図５に示していないが、超音波測定装置は、収差補正を行ってチャネルごとの合算データに基づいて測定地点に対する画像を生成してもよい。
又、超音波測定装置は、受信データからハーモニック（ｈａｒｍｏｎｉｃ）成分を抽出して用いてもよい。
さらに、超音波測定装置は、受信データに対してフィルタリング又は復調を行って、受信端ディレイ補正又は合成フォーカシング方法を適用してもよい。
詳細には、受信データは、フィルタリングを用いてバンドパスフィルタリングされたデータに加工されてハーモニック成分を抽出してもよく、パルスインバージョン（ｐｕｌｓｅ  ｉｎｖｅｒｓｉｏｎ）のような方式を用いてハーモニック成分が抽出されてもよく、復調によりベースバンド信号に変換されてもよい。
【００５９】
図１０は、本発明の一実施形態に係る実際フォーカシング受信データを説明するための図である。
例えば、超音波測定装置は、図１０に示すような実際フォーカシング受信データを取得する。
ここで、実際フォーカシング受信データは、超音波測定装置が実際に超音波信号をサンプルに送信し、これに対応してアコースティックフィールドを受信した場合のデータを示す。
実際フォーカシング受信データは、図１０に示すように、チャネルごとに受信データを含む。
【００６０】
図１０に示す実際フォーカシング受信データは、例えば、図６Ａに示すように、５個のフォーカシング地点（６１１〜６１５）にビームフォーカシングを行って取得された受信データであり得る。
第１地点６１１に対するビームフォーカシングの結果、受信された第１データは、第１チャネル受信データ乃至第５チャネル受信データ（１００１〜１００５）を含む。
第２地点６１２に対するビームフォーカシングの結果、受信された第２データは、第１チャネル受信データ乃至第５チャネル受信データ（１０１１〜１０１５）を含む。
フォーカシングの結果、受信された第５データは、第１チャネル受信データ乃至第５チャネル受信データ（１０２１〜１０２５）を含む。
【００６１】
図１１Ａ及び１１Ｂは、サンプル内の合成フォーカシング地点を説明するための概念図であり、図１２は、本発明の一実施形態に係る合成データ取得方法を説明するためのフローチャートである。
図１１Ａ及び１１Ｂは、図１２を参照してより詳しく説明する。
まず、図１２を参照すると、ステップＳ１２１０において、超音波測定装置は合成フォーカシング地点を決定する。
例えば、超音波測定装置は、図１１Ａに示すように、合成フォーカシング地点（１１０１〜１１２０）を決定する。超音波測定装置は、測定地点６０１の周辺領域に合成フォーカシング地点を決定する。
例えば、超音波測定装置は、フォーカシング地点（１１０１〜１１２０）がそれぞれ予め設定された間隔で離隔して全体的に長方形の格子形態を有するように合成フォーカシング地点を決定してもよい。
【００６２】
ステップＳ１２２０において、超音波測定装置は、合成フォーカシング地点（１１０１〜１１２０）それぞれに対応する受信データ取得のために格納された受信データのディレイを調整し、ステップＳ１２３０において、測定地点周辺領域の合成データを取得する。
例えば、図１１Ｂに示すように、超音波測定装置は、実際フォーカシング地点（６１１〜６１５）と合成フォーカシング地点１１０１との間のディレイがそれぞれＡ１、Ａ２、Ａ３、Ａ４、Ａ５であることを取得する。
超音波測定装置は、合成フォーカシング地点１１０１と実際フォーカシング地点（６１１〜６１５）との間で取得されたディレイ、及び少なくとも１つのトランスデューサー素子それぞれにおけるディレイに基づいて合成フォーカシング方法を行う。
【００６３】
詳細には、合成フォーカシング地点１１０１に対応する合成データは、図１０に示すような実際フォーカシング受信データに基づいて取得される。
例えば、トランスデューサー素子から実際フォーカシング地点（６１１〜６１５）までのフォーカシング時間を「ｔ１」と仮定する。さらに、上述したように、実際フォーカシング地点（６１１〜６１５）と合成フォーカシング地点１１０１との間のディレイはそれぞれＡ１、Ａ２、Ａ３、Ａ４、Ａ５である。
【００６４】
超音波測定装置は、第１実際フォーカシング地点に対応する受信データに“ｔ１−Ａ１”のディレイを適用させる。さらに、超音波測定装置は、第２実際フォーカシング地点から第５実際フォーカシング地点に対応する受信データそれぞれに“ｔ１−Ａ２”〜“ｔ１−Ａ５”のディレイを適用させる。
超音波測定装置は、ディレイが適用された第１実際フォーカシング地点に対応する受信データ乃至ディレイが適用された第５実際フォーカシング地点に対応する受信データを合算し、合成フォーカシング地点１１０１に対応する合成データを取得する。
【００６５】
一方、上述した“ｔ１−Ａ１”のディレイは、合成フォーカシング地点１１０１が実際フォーカシング地点（６１１〜６１５）より相対的にトランスデューサー素子に近い場合の例示である。
実際フォーカシング地点（６１１〜６１５）が合成フォーカシング地点１１０１より相対的にトランスデューサー素子に近い場合、ディレイは“ｔ１＋Ａ１”に決定される。
超音波測定装置は、上述したように合成データを取得し、そのため、図４Ｃに示すような、実際フォーカシング受信データと同一のデータを取得することができる。
【００６６】
すなわち、超音波測定装置は、合成フォーカシング方法を受信データに適用し、図４Ｂに示すように実際フォーカシングを行うことなく、図４Ｃに示すような受信データを取得することができる。
上述したように、本発明の一実施形態に係る超音波測定方法は、相対的に少ない実際のフォーカシング地点についてのみビームフォーカシングを行っても、合成フォーカシング方法を用いて測定地点周辺領域の相対的に多くのフォーカシング地点に対する受信データを取得することができる。
そのため、演算時間及び測定時間を著しく縮小することができる。
【００６７】
図１３は、本発明の他の実施形態に係る複数の合成フォーカシング地点に対する合成データを取得する方法を説明するためのフローチャートである。
図１３に示す実施形態は、超音波測定装置が実際フォーカシング地点からの受信データを格納したと仮定する。
ステップＳ１３１０において、超音波測定装置は複数の合成フォーカシング地点を決定する。
例えば、超音波測定装置は、測定地点近所領域で複数のフォーカシング地点を決定する。
【００６８】
ステップＳ１３２０において、超音波測定装置は、合成フォーカシング地点と実際フォーカシング地点との間のディレイを予測する。
超音波測定装置は、「ｍ」合成フォーカシング地点と実際フォーカシング地点との間の距離を予測する。
超音波測定装置は、「ｍ」合成フォーカシング地点と少なくとも１つのトランスデューサー素子までの測定時間と、実際フォーカシング地点と少なくとも１つのトランスデューサー素子までの測定時間との間の差に基づいて、トランスデューサー素子チャネルごとにディレイを予測する。そのため、トランスデューサー素子チャネルごとのディレイを予測する。
【００６９】
ステップＳ１３３０において、超音波測定装置は、予測されたディレイに基づいて受信データのディレイを補正する。
合わせて、ステップＳ１３４０において、超音波測定装置は、「ｍ」合成フォーカス地点に対応する合成データを取得して格納する。
超音波測定装置は、「ｍ」が「Ｍ」を超過するまで上述した工程を繰り返す（ステップＳ１３５０、Ｓ１３６０）。ここで、「Ｍ」は合成フォーカシング地点の個数である。
超音波測定装置は、ステップＳ１３７０において、測定地点周辺領域の複数フォーカシング地点に対応する合成データを取得する。
【００７０】
図１４は、本発明の一実施形態に係る合成データの例示である。
図１４に示すように、合成データは、例えば、２５個であってもよく、それぞれの合成データはチャネルごとのデータを含んでもよい。
例えば、第１合成データは、第１チャネルデータ乃至第５チャネルデータ（１４０１〜１４０５）を含み、第２合成データは第１チャネルデータ乃至第５チャネルデータ（１４１１〜１４１５）を含み、第３合成データは第１チャネルデータ乃至第５チャネルデータを含む。
合成データは、実際受信データを含んでもよい。例えば、超音波測定装置は、実際受信データ５個及び合成データ２０個を生成して２５個の受信データを生成してもよい。
又は、超音波測定装置は、実際受信データ５個に基づいて合成データ２５個を生成して２５個の受信データを生成してもよい。
【００７１】
図１５は、本発明の他の実施形態に係る超音波測定方法を説明するためのフローチャートである。
図１５は、特に受信データ測定時のサンプル内反射体が移動する場合、反射体動き補正を行う超音波測定方法を説明するためのフローチャートである。
ステップＳ１５１０において、超音波測定装置はサンプル上の複数地点に超音波信号ビームフォーカシングを行う。
【００７２】
例えば、超音波測定装置は、図６Ａ又は６Ｂに示すような、フロントフォーカシング方法又はリアフォーカシング方法に基づいてビームフォーカシングを行う。
超音波測定装置は、サンプル上の複数地点に送信側ビームフォーカシングを行う。
超音波測定装置は、トランスデューサー素子の個数よりも小さい個数のフォーカシング地点にビームフォーカシングを行う。
【００７３】
ステップＳ１５２０において、超音波測定装置は、サンプル上の複数地点に対するビームフォーカシングのそれぞれに対応して受信されるアコースティックフィールドを受信する。
超音波測定装置は、アコースティックフィールドを受信して受信データを取得及び格納する。
受信データは、トランスデューサー素子のチャネルごとのデータを含んでもよく、予め設定された時間の間測定されてもよい。
【００７４】
ステップＳ１５３０において、超音波測定装置は、チャネルごとのデータに対して動き補正を行う。
詳細には、超音波測定装置は、受信データに基づいて画像を生成する。超音波測定装置は、時間の流れに応じて複数の画像フレームを生成する。超音波測定装置は、複数の画像フレームのうち互いに隣接する画像フレームを比較して反射体の動きを検出する。
超音波測定装置は、反射体動きに対応してチャネルごとのデータを補正する。
【００７５】
例えば、超音波測定装置は、反射体が動く以前の位置と反射体が動いた以後の位置との間のディレイを予測する。
超音波測定装置は、予測されたディレイに基づいてチャネルごとのデータを補正して動き補正を行う。
動きベクトル検出は、画面次元で行われてもよく、動き補正はチャネルごとのデータのディレイを調整することによって行われてもよい。
【００７６】
ステップＳ１５４０において、超音波測定装置は合成データを取得する。
例えば、超音波測定装置はフォーカシング地点を決定し、合成フォーカシング地点と実際フォーカシング地点との間のディレイを予測する。
超音波測定装置は、合成フォーカシング地点と実際フォーカシング地点との間のディレイに基づいてチャネルごとのデータを補正し、そのため、合成フォーカシング地点に対応する合成データを取得する。
【００７７】
ステップＳ１５５０において、超音波測定装置は、合成データ又は実際フォーカシング受信データ及び合成データ両者に基づいて画像を生成する。
超音波測定装置は、チャネルごとにチャネルごとのデータを合算する。
上述した構成によって測定地点以外の他の反射体の影響を予防することができる。
【００７８】
また、超音波測定装置は、収差補正（ａｂｅｒｒａｔｉｏｎ  ｃｏｒｒｅｃｔｉｏｎ）を行ってもよい。
超音波測定装置は、収差補正された合算データに基づいて画像を生成する。
また、超音波測定装置は、合成フォーカシング又は動き補正によってディレイエラーを再算出してもよく、合成フォーカシングを再び行うことによって上記のステップが繰り返すように制御する。
【００７９】
図１６は、本発明の他の実施形態に係る動き補正を行う超音波測定方法を説明するためのフローチャートである。
図１６に示す動き補正の方法は、図１７Ａ〜１７Ｃと図１８を参照してより詳しく説明する。
図１６に示す動き補正の方法の実施形態で、サンプル１０内の反射体１７００は、図１７Ａに示すように移動するものと仮定する。
【００８０】
ステップＳ１６１０において、超音波測定装置は、サンプル上の複数地点に超音波信号ビームフォーカシングを行う。
例えば、超音波測定装置は、図６Ａ又は図６Ｂに示すような、フロントフォーカシング方法又はリアフォーカシング方法に基づいてビームフォーカシングを行う。
ステップＳ１６２０において、超音波測定装置は、サンプル上の複数地点に対するビームフォーカシングそれぞれに対応して受信されるアコースティックフィールドを受信する。
超音波測定装置は、アコースティックフィールドを受信して受信データを取得及び格納を行う。受信データはトランスデューサー素子のチャネルごとのデータを含んでもよく、予め設定された時間の間測定されてもよい。
【００８１】
ステップＳ１６３０において、超音波測定装置は、受信データに基づいて画像フレームを生成する。
上述したように、受信データは予め設定された時間の間に測定され、そのために時間の流れに応じて格納される。
超音波測定装置は、時間の流れにより格納された受信データに基づいて予め設定された時間の間の複数の画像フレームを生成する。
例えば、超音波測定装置は、図１７Ｂに示すように、第１、２、３画像フレーム（１７１０、１７２０、１７３０）を生成する。
【００８２】
超音波測定装置は、第１、２、３画像フレーム（１７１０、１７２０、１７３０）での反射体の動きを検出する。
例えば、超音波測定装置は、第１画像フレーム１７１０での反射体１７１１と第２画像フレーム１７２０での反射体１７２１との間の距離が閾値以上であることを確認する。
ステップＳ１６４０において、超音波測定装置は、隣接する画像フレームの反射体位置の間の距離が閾値以上であるかを判断して動きベクトルを抽出する。
さらに、超音波測定装置は、第３画像フレーム１７３０での反射体１７３１も第１画像フレーム１７１０と比較して動きを確認する。
【００８３】
例えば、超音波測定装置は、第３画像フレーム１７３０と第２画像フレーム１７２０とを比較し、動きベクトルが「０」であることを確認する。
これによって、超音波測定装置は、第３画像フレーム１７３０の反射体１７３１と第１画像フレーム１７１０の反射体１７１１との位置間の動きベクトルと、第２画像フレーム１７２０の反射体１７２１と第１画像フレーム１７１０の反射体１７１１との位置間の動きベクトルとが同一であることを確認できる。
【００８４】
一方、図に示していないが、第４画像フレームで反射体が再び動いてもよい。
このような場合、超音波測定装置は、第４画像フレームと第３画像フレームとを比較して動きベクトルを検出する。
超音波測定装置は、第４画像フレームでの反射体と第１画像フレームでの反射体との間の動きベクトルを取得する。
詳細には、超音波測定装置は、第３フレームと第１フレームとの間の動きベクトルに、第４フレームと第３フレームとの間の動きベクトルを合算し、第４フレームと第１フレームとの間の動きベクトルを取得する。
【００８５】
ステップＳ１６５０において、超音波測定装置は、動きベクトルに基づいて受信データの動きを補正する。
例えば、超音波測定装置は、図１７Ｃに示すように、移動以前地点１７４１と移動以後地点１７４２との間のディレイ「Δｔ」を予測する。
超音波測定装置は、移動以前地点１７４１と移動以後地点１７４２との間のディレイ「Δｔ」に基づいて、トランスデューサー素子チャネルごとのディレイを予測することができる。超音波測定装置は、予測されたトランスデューサー素子チャネルごとのディレイに基づいて、受信データをチャネルごとにディレイ補正を行う。それにより、受信データは動き補正される。
【００８６】
図１８は、チャネルごとにディレイを補正する工程を説明するための図である。
図１８は、例えば、第１地点にビームフォーカシングを行った結果に応じて取得された受信データである。
受信データはｔ＝１、２、．．．、Ｔである場合に対して格納される。ここで、「ｔ」はサンプリング区間に対する時間インデックスである。さらに、各時間ごとの受信データは、チャネルごとのデータを含む。
例えば、ｔ＝１である場合の受信データは、第１チャネル受信データ乃至第５チャネル受信データ（１８０１〜１８０５）を含む。ｔ＝２である場合の受信データは、第１チャネル受信データ乃至第５チャネル受信データ（１８１１〜１８１５）を含んでもよい。ｔ＝Ｔである場合の受信データは、第１チャネル受信データ乃至第５チャネル受信データ（１８２１〜１８２５）を含む。
【００８７】
一方、例えば、図１７Ａ及び１７Ｂに示す実施形態では、第１サンプリング区間と第２サンプリング区間との間でのみ反射体が移動した。さらに、第２サンプリング区間以後の第Ｔサンプリング区間まで反射体が停止した状況を想定する。
この場合、超音波測定装置は、第２サンプリング区間に対応する受信データから第Ｔサンプリング区間に対応する受信データのチャネルごとのディレイ補正を行う。
詳細には、超音波測定装置は、各サンプリング区間でチャネルごとのディレイを予測し、予測されたチャネルごとのディレイに基づいてチャネルごとのディレイ補正を行う。
そのため、超音波測定装置は、第２サンプリング区間に対応するデータ乃至第Ｔサンプリング区間に対応する受信データが、反射体が初期位置に継続して位置したような受信データになるよう補正する。
【００８８】
例えば、超音波測定装置は、図１７Ａに示すように、反射体１７００が初期位置に続けて位置するように受信データを補正する。
それにより、超音波測定装置は、図１８に示すように、第２サンプリング区間のチャネルごとのデータを動き補正した、動き補正された受信データ（１８３１〜１８３５）と、第Ｔサンプリング区間のチャネルごとのデータを動き補正した、動き補正された受信データ（１８４１〜１８４５）を生成する。
従って、第２サンプリング区間の受信データと第Ｔサンプリング区間の受信データとの間の受信データも動き補正されたことを、当業者は容易に理解するのであろう。
【００８９】
再び図１６を参照すると、ステップＳ１６６０において、超音波測定装置は、動き補正された受信データに基づいて合成データを取得する。
また、ステップＳ１６７０において、超音波測定装置は、合成データに基づいて画像を生成する。
超音波測定装置は、チャネルごとにチャネルごとのデータを合算する。
上述した構成によって測定地点以外の他の反射体の影響を予防することができる。
超音波測定装置は、収差補正を行ってもよい。超音波測定装置は、収差補正された合算データに基づいて画像を生成してもよい。
上述した工程によって、サンプル内に動いている反射体が存在する場合にも動きを補正することによって明確な画像生成を行うことができる。
【００９０】
図１９は、本発明の一実施形態に係る動き補正を行う超音波測定装置の概略構成を示すブロック図である。
図１９に示すように、超音波測定装置は、超音波信号送受信部１９１０、動き補正部１９２０、合成データ取得部１９３０、データ合算及び収差補正部１９４０、画像生成部１９５０、動き推定部１９６０、及び表示部１９７０を備える。
【００９１】
超音波信号送受信部１９１０は、サンプルの少なくとも１つの実際フォーカシング地点に超音波信号を送信し、ビームフォーカシングに対応して反射体が生成するアコースティックフィールドを受信する。
【００９２】
画像生成部１９５０は、受信したアコースティックフィールドに対応する受信データに基づいて複数の画像フレームを生成する。
受信データは、上述したように時系列的なデータであってもよい。例えば、受信データは表１のように、サンプリングナンバー及びトランスデューサー素子チャネルごとのデータを含み得る。そのため、画像生成部１９５０は、時間の流れに応じた複数の画像フレームを生成することができる。
【００９３】
動き補正部１９２０は、生成された複数の画像フレーム間の比較によって動く反射体を補正する。
詳細には、動き推定部１９６０で推定された動きベクトルに基づいて動き補正部１９２０は受信データをチャネルごとに動き補正を行う。
合成データ取得部１９３０は、動き補正された受信データに合成フォーカシング方法を適用する。そのため、合成データ取得部１９３０は、少なくとも１つの合成データを取得する。
【００９４】
データ合算及び収差補正部１９４０は、少なくとも１つの合成データをチャネルごとに合算して単一の受信データを取得する。
さらに、データ合算及び収差補正部１９４０は、単一の受信データに収差補正を行う。
画像生成部１９５０は、収差補正された単一の受信データに基づいて画像を生成する。
表示部１９７０は、生成された画像を表示する。
【００９５】
図２０Ａ及び２０Ｂは、収差補正工程を説明するための概念図である。
図２０Ａに示すように、（ｘ、ｚ）座標を有する反射体２００１からアコースティックフィールド２０４１が放出される。
反射体２００１は、トランスデューサー素子アレイ２００２から“Ｒ”だけ離隔される。ここで、アコースティックフィールド２０４１は球面波の形態であってもよく、これによってトランスデューサー素子それぞれに到達する時間が異なってもよい。
【００９６】
トランスデューサー素子は、受信したアコースティックフィールド２０４１を電気的な信号２００３に変換して出力する。
収差補正部２００４は、受信した電気的な信号２００３のディレイを調整して収差補正された電気的な信号２０１０を出力する。
信号合算部２００５は、収差補正された電気的な信号２０１０を合算して画像信号２００６を生成する。
【００９７】
一方、図２０Ｂに示すように、音速度が一定ではない媒質２０３０が配置されている場合、アコースティックフィールド２０４２がトランスデューサー素子のそれぞれに球面波の形態ではない不規則な時間に到達する。
そのため、変換された電気的な信号２０１３に対して収差補正部２００４が収差補正しても収差補正された信号２０４０の間にディレイが発生することがある。
そのため、画像信号２０６０は、反射体２００１に対する歪曲された画像信号となる。
【００９８】
本発明の一実施形態に係る収差補正部２００４は、中央チャネル信号を他のチャネル信号と相互相関（ｃｒｏｓｓ  ｃｏｒｒｅｌａｔｉｏｎ）して他のチャネル信号が中央チャネル信号と整列するように収差補正してもよい。
例えば、収差補正部２００４は、相互相関して図２０Ｃに示すように、他のチャネル信号のディレイ「Δｔ」を調整して中央チャネル信号と整列させる。
図２０Ｄ及び２０Ｅは、本発明の実施形態に係る相互相関に基づいた収差補正以前と以後の画像写真である。
【００９９】
図２１は、本発明のさらに他の実施形態に係る超音波測定方法を説明するためのフローチャートである。
超音波測定装置は、ステップＳ２１１０において、少なくとも１つの実際フォーカシング地点に対応する受信データに合成フォーカシング方法を適用し、反射体周辺領域の合成フォーカシング地点に対応する少なくとも１つの合成データを取得する。
超音波測定装置は、ステップＳ２１２０において、少なくとも１つの合成データを用いて反射体の画像を生成する。
【０１００】
以上で説明した装置は、ハードウェア構成要素、ソフトウェア構成要素、及び／又はハードウェア構成要素及びソフトウェア構成要素の組み合わせで実現され得る。
例えば、プロセッサ、コントローラ、ＡＬＵ（ａｒｉｔｈｍｅｔｉｃ  ｌｏｇｉｃ  ｕｎｉｔ）、デジタル信号プロセッサ（ｄｉｇｉｔａｌ  ｓｉｇｎａｌ  ｐｒｏｃｅｓｓｏｒ）、マイクロコンピュータ、ＦＰＡ（ｆｉｅｌｄ  ｐｒｏｇｒａｍｍａｂｌｅ  ａｒｒａｙ）、ＰＬＵ（ｐｒｏｇｒａｍｍａｂｌｅ  ｌｏｇｉｃ  ｕｎｉｔ）、マイクロプロセッサー、又は、命令（ｉｎｓｔｒｕｃｔｉｏｎ）を実行して応答できる異なる装置のように、１つ以上の汎用コンピュータ又は特殊目的のコンピュータを用いて実現してもよい。
【０１０１】
処理装置は、オペレーションシステム（ＯＳ）及びオペレーションシステム上で行われる１つ以上のソフトウェアアプリケーションを実行し得る。また、処理装置は、ソフトウェアの実行に応答してデータに対しアクセス、格納、操作、処理、及び生成を行い得る。
理解の便宜のために、処理装置は１つの使用されるものとして説明した場合もあるが、当該の技術分野で通常の知識を有する者であるならば、処理装置が複数の処理要素（ｐｒｏｃｅｓｓｉｎｇ  ｅｌｅｍｅｎｔ）及び／又は複数の類型の処理要素を含んでいることが理解されよう。例えば、処理装置は、複数のプロセッサ又は１つのプロセッサ及び１つのコントローラを含んでもよい。また、並列プロセッサ（ｐａｒａｌｌｅｌ  ｐｒｏｃｅｓｓｏｒ）のような、他の処理構成も可能である。
【０１０２】
ソフトウェアは、コンピュータプログラム、コード、命令、又はこのうちの１つ以上の組み合わせを含んでもよく、所望の通りに動作するよう処理装置を構成したり、独立的又は結合的に処理装置を命令してもよい。
ソフトウェア及び／又はデータは、処理装置によって解釈されたり処理装置に命令又はデータを提供するためどのような類型の機械、構成要素、物理的装置、仮想装置、コンピュータ格納媒体又は装置、送信される信号波に永久的又は一時的に具体化できる。
ソフトウェアは、ネットワークに接続されたコンピュータシステム上に分散し、分散された方法で格納されたり実行されてもよい。ソフトウェア及びデータは１つ以上のコンピュータで読取可能な記録媒体に格納されてもよい。
【０１０３】
上述の本発明の実施形態に係る方法は、多様なコンピュータ手段を介して様々な処理を実行することができるプログラム命令の形態で実現され、コンピュータ読取可能な記録媒体に記録されてもよい。
コンピュータ読取可能な記録媒体は、プログラム命令、データファイル、データ構造などのうち１つ又はその組み合わせを含んでもよい。記録媒体に記録されるプログラム命令は、本発明の目的のために特別に設計されて構成されたものでもよく、コンピュータソフトウェア分野の技術を有する当業者にとって公知のものであり、使用可能なものであってもよい。
【０１０４】
コンピュータ読取可能な記録媒体の例としては、ハードディスク、フロッピー（登録商標）ディスク及び磁気テープのような磁気媒体、ＣＤ−ＲＯＭ、ＤＶＤのような光記録媒体、光ディスクのような光磁気媒体、及びＲＯＭ、ＲＡＭ、フラッシュメモリなどのようなプログラム命令を保存して実行するように特別に構成されたハードウェア装置が含まれてもよい。
プログラム命令の例には、コンパイラによって作られるような機械語コードだけでなく、インタープリタなどを用いてコンピュータによって実行できる高級言語コードが含まれる。前記したハードウェア装置は、本発明の動作を行うために１つ以上のソフトウェアモジュールとして作動するように構成されてもよく、その逆も同様である。
【０１０５】
尚、本発明は、上述の実施形態に限られるものではない。本発明の技術的範囲から逸脱しない範囲内で多様に変更実施することが可能である。
【符号の説明】
【０１０６】
１０、６０    サンプル
１８    超音波測定装置
２０    表示部
３０    プロセッサ
４０    媒質コンテナ
５０    送受信部
５１〜５４    トランスデューサー素子
７０    増幅器及びフィルタ
８０    アナログ−デジタル変換部
９０    受信部
２０１、２８０、２９０    反射体
２１０、２２０、２３０    アコースティックフィールド
【図１Ａ】
2015150423.tif 000004
【図１Ｂ】
2015150423.tif 000005
【図２Ａ】
2015150423.tif 000006
【図２Ｂ】
2015150423.tif 000007
【図３】
2015150423.tif 000008
【図４Ａ】
2015150423.tif 000009
【図４Ｂ】
2015150423.tif 000010
【図４Ｃ】
2015150423.tif 000011
【図４Ｄ】
2015150423.tif 000012
【図４Ｅ】
2015150423.tif 000013
【図４Ｆ】
2015150423.tif 000014
【図５】
2015150423.tif 000015
【図６Ａ】
2015150423.tif 000016
【図６Ｂ】
2015150423.tif 000017
【図７Ａ】
2015150423.tif 000018
【図７Ｂ】
2015150423.tif 000019
【図７Ｃ】
2015150423.tif 000020
【図８】
2015150423.tif 000021
【図９】
2015150423.tif 000022
【図１０】
2015150423.tif 000023
【図１１Ａ】
2015150423.tif 000024
【図１１Ｂ】
2015150423.tif 000025
【図１２】
2015150423.tif 000026
【図１３】
2015150423.tif 000027
【図１４】
2015150423.tif 000028
【図１５】
2015150423.tif 000029
【図１６】
2015150423.tif 000030
【図１７Ａ】
2015150423.tif 000031
【図１７Ｂ】
2015150423.tif 000032
【図１７Ｃ】
2015150423.tif 000033
【図１８】
2015150423.tif 000034
【図１９】
2015150423.tif 000035
【図２０Ａ】
2015150423.tif 000036
【図２０Ｂ】
2015150423.tif 000037
【図２０Ｃ】
2015150423.tif 000038
【図２０Ｄ】
2015150423.tif 000039
【図２０Ｅ】
2015150423.tif 000040
【図２１】
2015150423.tif 000041