(57)【要約】      （修正有）【課題】スペクトル拡散広瞬時帯域幅コヒーレント符号化波形を用いる合成アパーチャ超音波イメージングのための技術、システム及びデバイスを提供する。【解決手段】異なる周波数帯に対応し、対応する位相による固有の周波数を有する互いに直交する複数の個別の直交符号化波形から形成された複合波形を合成することと、１つ以上の送信位置からターゲットに向けて、複合波形に基づく音響波を送信することと、１つ以上の受信位置において、送信された音響波形に対応する、ターゲットの少なくとも一部から戻る音響エネルギを受信することとを含み、送受信位置は、それぞれ、ターゲットに対するトランスデューサのアレイの空間位置と、アレイのビーム位相中心位置との一方又は両方を含み、送信された音響波形及び応答音響波形は、拡大された実効アパーチャを生成する。【選択図】図１Ｃ
2018183651.tif 000002
【特許請求の範囲】
【請求項１】
音響イメージングデバイスにおいて音響波形を生成する方法であって、
１つ以上の波形シンセサイザにおいて、ターゲットに向けて送信される１つ以上の複合波形を合成することであって、前記複合波形は、それぞれが互いに直交し及び異なる周波数帯に対応する個別の直交符号化波形から形成され、前記各個別の直交符号化波形は、対応する位相による固有の周波数を有することと、
前記ターゲットに対する１つ以上の送信位置から、複数の音響波形からなる１つ以上の複合音響波形を送信することであって、前記送信することは、アレイの１つ以上のトランスデューサ素子を選択し、前記１つ以上の複合音響波形のそれぞれの前記複数の個別の直交符号化波形を、前記１つ以上の複合音響波形のそれぞれの複数の対応する音響波形にトランスデュースすることを含むことと、
前記ターゲットに対する１つ以上の受信位置において、前記送信される音響波形に対応するターゲットの少なくとも一部から戻る応答音響波形を受信することであって、前記受信することは、前記アレイのトランスデューサ素子の少なくとも幾つかを選択して、前記応答音響波形を受信することを含むことと、を含み、
前記送信位置及び前記受信位置は、それぞれ、前記ターゲットに対するトランスデューサ素子のアレイの空間位置と、前記アレイのビーム位相中心位置との一方又は両方を含み、
前記送信される音響波形及び応答音響波形は、前記音響イメージングデバイスの拡大された実効アパーチャを生成する、方法。
【請求項２】
前記音響波形をターゲットに送信する際、前記アレイのトランスデューサ素子を制御して、前記ターゲットに対する前記複合波形の向きを変更し、前記ターゲットが、イメージング期間に亘って、異なる波形の向きで前記音響波形を受信するようにする請求項１記載の方法。
【請求項３】
前記トランスデューサ素子の制御は、前記ターゲットに対して、前記複数の空間位置に沿って前記アレイを平行移動することによって、前記ターゲットに対する前記複合波形の向きを変更することを含む請求項２記載の方法。
【請求項４】
前記トランスデューサ素子の制御は、前記アレイの１つ以上のトランスデューサ素子における前記送信される音響波形のビーム位相中心位置を変更することによって、前記ターゲットに対する前記複合波形の向きを変更することを含む請求項２記載の方法。
【請求項５】
前記複数の個別の直交符号化波形の各波形は、それぞれが個別に振幅重み付け及び位相重み付けされた複数の振幅及び複数の位相を含む請求項１記載の方法。
【請求項６】
前記複合波形の前記個別の直交符号化波形を合成することは、周波数帯を選択することと、各直交符号化波形の１つ以上の振幅、時間−帯域幅積パラメータ及び位相パラメータを判定することとを含む請求項１記載の方法。
【請求項７】
前記位相パラメータは、一組の擬似乱数又は一組の決定論的な数から決定される請求項６記載の方法。
【請求項８】
前記複数の音響波形の送信は、前記複数の位置の少なくとも１つの位置から前記音響波形を順次的に又は同時に送信することを含む請求項１記載の方法。
【請求項９】
前記個別の直交符号化波形は、コヒーレント波形を含む請求項１記載の方法。
【請求項１０】
前記複合波形に基づいて無線周波数（ＲＦ）波形を形成することを更に含み、
前記送信される音響波形は、前記アレイの１つ以上のトランスデューサ素子においてＲＦベースの複合波形をトランスデュースすることによって生成される請求項１記載の方法。
【請求項１１】
ＲＦベースの複合波形を増幅することを更に含む請求項１０記載の方法。
【請求項１２】
前記受信した応答音響波形を増幅することを更に含む請求項１記載の方法。
【請求項１３】
前記受信した応答音響波形を、前記１つ以上の複合波形に対応する１つ以上の受信複合波形として、アナログフォーマットからデジタルフォーマットに変換することを更に含み、各受信複合波形は、前記ターゲットの情報を含み、
前記情報は、前記受信複合波形の対応する周波数帯に関連する振幅及び位相を含む請求項１記載の方法。
【請求項１４】
前記受信複合波形の対応する周波数帯の少なくとも１つの周波数帯について、前記振幅及び位相の少なくとも１つは、それぞれ個別に、振幅重み付け又は位相重み付けされる請求項１３記載の方法。
【請求項１５】
前記受信した応答音響波形を処理して、前記ターゲットの少なくとも一部の画像を生成することを更に含む請求項１記載の方法。
【請求項１６】
前記受信した応答音響波形を、前記１つ以上の複合波形に対応する１つ以上の受信複合波形として、アナログフォーマットからデジタルフォーマットに変換することを更に含み、各受信複合波形は、前記ターゲットの情報を含み、
前記１つ以上の複合波形及び前記対応する１つ以上の受信複合波形をデータブロックのメモリマップ内に保存することを更に含み、各データブロックは、各サンプル点についての前記複合波形の受信した応答音響波形、前記対応する個別の直交符号化波形、及び前記サンプル点についての前記１つ以上のトランスデューサ素子の対応する位置データを保存する請求項１５記載の方法。
【請求項１７】
前記処理は、
前記保存された受信複合波形の軸範囲圧縮を実行することと、
周波数領域処理技術又は時間領域処理技術の１つ以上を用いて、前記実効アパーチャを含む軸範囲圧縮データの各保存されたブロックを処理することによって合成アパーチャ画像を形成することと、を含む請求項１５記載の方法。
【請求項１８】
波形発生器に接続された１つ以上の波形シンセサイザを有する波形生成ユニットであって、前記波形発生器が提供する波形情報に応じて前記１つ以上の波形シンセサイザによって生成される異なる周波数帯に対応する複数の個別の直交符号化波形からなる複合波形を合成し、前記個別の直交符号化波形は、それぞれ互いに直交し及び異なる周波数帯に対応し、前記各個別の直交符号化波形は、対応する位相による固有の周波数を有する波形生成ユニットと、
送信モードと受信モードの間で切り換えられる送信／受信スイッチングユニットと、
前記送信／受信スイッチングユニットと通信して、ターゲットに対する１つ以上の送信位置から複数の音響波形からなる複合音響波形を送信し、前記複合音響波形の送信される音響波形は、前記複合波形の合成された個別の直交符号化波形に基づいており、及び前記ターゲットに対する１つ以上の受信位置において、前記ターゲットの少なくとも一部から戻る前記複数の送信された音響波形に対応する応答音響波形を受信し、前記送信される音響波形及び応答音響波形は、合成アパーチャ音響波形イメージングシステムの拡大された実効アパーチャを生成し、前記送信位置及び前記受信位置は、それぞれ、空間位置と、ビーム位相中心位置との一方又は両方を含むトランスデューサ素子のアレイと、
前記トランスデューサ素子のアレイと通信して、前記複数の個別の直交符号化波形を前記複数の対応する音響波形にトランスデュースする前記アレイの１つ以上のトランスデューサ素子を選択し、及び前記応答音響波形を受信する前記アレイの１つ以上のトランスデューサ素子を選択する多重化ユニットと、
前記ターゲットの情報を提供する、前記トランスデューサ素子のアレイによって受信された前記応答音響波形を、アナログフォーマットからデジタルフォーマットに変換するアナログ／デジタル（Ａ／Ｄ）変換器のアレイと、
前記波形生成ユニット及びＡ／Ｄ変換器のアレイと通信し、データを保存するメモリユニットと、前記メモリユニットに接続されて、前記情報をデータとして処理する処理ユニットとを備えるコントローラユニットと、
前記コントローラユニットと通信するユーザインタフェースユニットとを備える合成アパーチャ音響波形イメージングシステム。
【請求項１９】
前記保存されるデータは、前記受信した応答音響波形の前記デジタルフォーマット、前記対応する合成された複合波形、及び前記１つ以上の送信位置及び受信位置における前記トランスデューサ素子のアレイの対応する位置データを含む請求項１８記載のシステム。
【請求項２０】
前記波形生成ユニットは、前記送信／受信スイッチングユニットと前記１つ以上の波形シンセサイザとの間に配設され、前記複合波形を変更する１つ以上の増幅器を更に備える請求項１８記載のシステム。
【請求項２１】
前記送信／受信スイッチングユニットと前記Ａ／Ｄコンバータのアレイとの間に配設され、前記受信した応答音響波形を変更する１つ以上の前置増幅器のアレイを更に備える請求項１８記載のシステム。
【請求項２２】
前記処理ユニットは、デジタル信号プロセッサを備える請求項１８記載のシステム。
【請求項２３】
前記コントローラユニットは、前記音響波形画像処理システムの要素の少なくとも１つと時間的に同期するマスタクロックを更に備える請求項１８記載のシステム。
【請求項２４】
前記コントローラユニットは、前記情報を処理して、前記ターゲットの少なくとも一部の画像を生成する請求項１８記載のシステム。
【請求項２５】
前記ユーザインタフェースユニットは、前記画像を表示するディスプレイと、前記システムの動作のための動作モードを含むユーザ入力データを受信するユーザ入力端子とを備える請求項１８記載のシステム。
【請求項２６】
前記動作モードは、前記応答音響波形から取得される少なくとも１つの測定された特性の少なくとも１つの特徴に基づいて、画像色分けを可能にする生物組織のイメージングのための少なくとも１つの人工組織染色（Artificial Tissue Staining：ＡＴＳ）モードを含む請求項２５記載のシステム。
【請求項２７】
前記動作モードは、前記応答音響波形から取得される１つ以上の測定された特性の少なくとも１つの特徴に基づいて、１つ以上のアルゴリズム分類器を用いて生物組織タイプを分類して生物組織をイメージングする少なくとも１つのコンピュータ支援診断（Computer Aided Diagnostic-Mode：ＣＡＤ）モードを含む請求項２５記載のシステム。
【請求項２８】
前記ディスプレイは、前記分類された生物組織タイプに基づいて、前記生物組織の色分けされた画像を表示する請求項２５記載のシステム。
【請求項２９】
前記音響波形をターゲットに送信する際、前記コントローラユニットは、前記トランスデューサ素子のアレイを制御して、前記ターゲットに対する前記複合波形の向きを変更し、前記ターゲットが、イメージング期間に亘って、異なる波形向きで前記複合音響波形を受信するようにする請求項１８記載のシステム。
【請求項３０】
前記トランスデューサ素子のアレイは、前記送信位置に沿って１次元、２次元又は３次元で移動して、前記複数の音響波形を送信し、及び受信位置に沿って１次元、２次元又は３次元で移動して、前記応答音響波形を受信するように動作する請求項１８記載のシステム。
【請求項３１】
前記アレイのトランスデューサ素子の少なくとも１つは、前記送信アレイの他のトランスデューサ素子から独立して、１次元、２次元又は３次元で移動することができる請求項１８記載のシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
関連出願への相互参照
この特許文献は、２０１３年９月１３日に出願された、米国仮出願番号第６１／８７７，８８４号、発明の名称「COHERENT SPREAD SPECTRUM CODED WAVEFORMS IN SYNTHETIC APERTURE IMAGE FORMATION」の優先権の利益を主張し、この文献は、参照によって本文献の一部として援用される。
技術分野
この特許文献は、音響画像形成（acoustic image formation）に関する。
【背景技術】
【０００２】
音響イメージングは、媒体を伝播する音波の特性を用いて、視覚的な画像を生成するイメージング方式である。高周波音響イメージングは、数十年に亘って、様々な生医学的分野において、動物及び人間の内部構造及び機能を視覚化するイメージング方式として使用されている。生医学的なイメージングにおいて使用される高周波音響波は、例えば、１〜２０ＭＨｚ又はこれより高い異なる周波数で動作することができ、超音波と呼ばれることも多い。従来の超音波イメージングの技術は、不十分な空間的分解能及び組織弁別を含む幾つかの要因のために、望ましい画質が得られず、多くの臨床的適応又は応用の用途が制限されることがある。
【発明の概要】
【課題を解決するための手段】
【０００３】
コヒーレントな、スペクトル拡散、広瞬時帯域幅、周波数符号化及び／又は位相符号化音響波形を用いる合成アパーチャ超音波イメージングのための技術、システム及びデバイスを開示する。ここに開示する技術、システム及び装置を用いて、生物組織の１次元（１Ｄ）、２次元（２Ｄ）及び／又は３次元（３Ｄ）超音波画像を形成することができる。
【０００４】
本開示に記述する主題は、以下の特徴の１つ以上を提供することができ、多くの用途に適用できる。例えば、ここに開示する技術を使用することによって、日常的なプライマリケアのスクリーニングにおいて、悪性腫瘍を含む初期疾患並びに進行期癌を特定及び検出することができ、これにより、潜在的に、診断が困難な無症候性患者の生存率を高めることができる。また、認定放射線科医がここに開示する技術を使用することによって、外科生検又は切除侵襲を行う前に、腫瘍が良性か悪性かを診断でき、不必要な生検を減らし、患者の生存率を向上させることができる。また、ここに開示する技術は、細針生検器具に統合されて、医療的処置において、非侵襲的診断を検証することができ、これによって、生検処置の侵襲度を低めることができる。また、ここに開示する技術は、手術前にコンピュータ断層撮影（Computed Tomography：ＣＴ）Ｘ線撮像に統合され、及び／又は手術中に最低侵襲外科高解像度ビデオ測定器に統合され、ＣＴ画像、光学画像及び超音波画像を融合することができ、これによって、外科医は、健康な組織に損傷を与えることなく疾患のある組織を特定し、診断して、外科的に切除することができる。ここに開示する技術は、専門的な外科器具に統合され、超音波画像を他のデータに融合することができ、外科医は、近傍の構造への不必要なダメージを最小化しながら、解剖学的関心領域を特定し、処置することができる。ここに開示する技術により、例えば、密封放射性線源及びカテーテルの適切な位置への挿入を正確に案内することによって、悪性腫瘍の小線源治療の時間を短縮できる。同様に、ここに開示する技術により、疾患の高用量の局所的な医薬投与治療を補助することができる。
【図面の簡単な説明】
【０００５】
【図１Ａ】超音波イメージングにおいて実際のアパーチャビームを形成するトランスデューサアレイによって生成される超音波ビームを示す図である。
【図１Ｂ】位相中心が連続的な位置にあり、超音波イメージングにおける合成アパーチャビームを形成するトランスデューサアレイによって生成される超音波ビームを示す図である。
【図１Ｃ】複数の送信及び／又は受信位置から合成アパーチャビームを形成するトランスデューサアレイによって生成される例示的な合成超音波送信及び／又は受信ビームを示す図である。
【図１Ｄ】複数の送信及び／又は受信位置から合成アパーチャビームを形成するトランスデューサアレイによって生成される例示的な合成超音波送信及び／又は受信ビームを示す図である。
【図１Ｅ】ここに開示する技術の合成超音波ビームを生成するための例示的なトランスデューサアレイを示す図である。
【図２Ａ】スペクトル拡散符号化波形を用いた例示的な合成アパーチャ超音波システムのブロック図である。
【図２Ｂ】スペクトル拡散符号化波形を用いた他の例示的な合成アパーチャ超音波システムのブロック図である。
【図２Ｃ】ここに開示する技術の例示的な合成アパーチャイメージング技術を示す図である。
【図２Ｄ】ここに開示する技術の例示的な合成アパーチャイメージング技術を示す図である。
【図２Ｅ】スペクトル拡散符号化波形を用いた例示的な合成アパーチャ超音波画像処理システムの動作を示す図である。
【図３】複数の波形から構成される例示的なスペクトル拡散、広瞬時帯域幅、周波数及び／又は位相符号化複合波形のグラフ図である。
【図４】例示的なスペクトル拡散符号化波形の曖昧度関数の特徴を示す図である。
【図５Ａ】ビーム操縦、動的フォーカシング、及び整形の例を示す図である。
【図５Ｂ】ビーム操縦、動的フォーカシング、及び整形の例を示す図である。
【図５Ｃ】ビーム操縦、動的フォーカシング、及び整形の例を示す図である。
【図６】超音波画像を生成するための例示的な合成アパーチャ超音波信号処理技術のブロック図である。
【図７】例示的な合成アパーチャによりサンプリングされた受信エコーデータ及びメモリへの保存を示す図である。
【図８】広帯域、スペクトル拡散、雑音状、コヒーレント合成アパーチャ画像形成のための例示的な波数アルゴリズムのブロック図である。
【図９】ストルト変換処理の例示的な実施例を示す図である。
【図１０】例示的な反復、閉ループ、合成アパーチャ画像自動フォーカスサブルーチンのブロック図である。
【発明を実施するための形態】
【０００６】
超音波イメージングは、２つの媒体（例えば、生物組織構造）の間の境界において一部が反射し、一部が透過する物理的弾性媒体内に音響波形（パルス）を出射することによって実行できる。反射は、２つの媒体の間の（例えば、２つの異なる生物組織タイプ間の界面における）音響インピーダンスの差に依存する。例えば、送信された音響波形の音響エネルギの一部は、界面において散乱してトランスデューサに戻り、受信及び処理されて情報が抽出され、残りのエネルギは、次の媒体に伝播する。幾つかの例では、散乱中心として機能する反射性媒体に含まれる２つ以上のインピーダンスの結果として反射の散乱が生じることがある。更に、例えば、音響エネルギは、媒体の特性及び／又は音響波の性質に基づいて屈折し、回折し、遅延し及び／又は減衰することがある。
【０００７】
幾つかの既存の超音波イメージングシステムでは、典型的なトランスデューサは、圧電素子のアレイを使用して、ターゲット領域（例えば、生体組織）に向けて超音波パルスを送信し、内部の散乱構造から戻る反射された超音波信号（エコー）を受信する。このトランスデューサアレイは、イメージングシステムのアパーチャとして機能する。超音波パルスは、シーケンスパルスとして、平面又は立体において電子的に操縦及びフォーカシングでき、ターゲットの画像を形成するために用いられる応答エコーの１Ｄ、２Ｄ及び／又は３Ｄマップを生成するために使用することができる。超音波パルスを操縦及びフォーカシングするプロセスは、ビーム整形とも呼ばれる。幾つかの具体例では、トランスデューサアレイにおいて送信及び受信される超音波パルス及び応答エコーは、フェーズドアレイとして動作するアレイの各トランスデューサにおいて、個別に時間的に遅延させることができる。
【０００８】
従来の実際のアパーチャ超音波イメージングシステムでは、画像の品質は、超音波システムのトランスデューサで発生する音場に直接的に依存する。図１Ａは、例えば、単一の電気音響トランスデューサ素子又は複数の電気音響トランスデューサ素子のアレイとして構成することができる幅ｗＲのトランスデューサアレイ１１０によって生成され、関心ボリューム（volume of interest：ＶＯＩ）１１５に向けられる超音波ビーム１１１を形成するための実際のアパーチャを形成する実際のアパーチャ超音波ビーム１１１を表す図である。例えば、画像は、通常、１軸の画像線毎に順次的に取得され（すなわち、スライス毎のターゲット領域範囲の走査）、イメージングの間、フレームレートに制約が生じ、これは、例えば、動くターゲットのイメージングを含む様々なリアルタイム超音波画像化アプリケーションにおいて、問題となることがある。
【０００９】
従来の実際のアパーチャ超音波イメージングにおける制約の問題を解決するために、合成アパーチャ超音波イメージング（synthetic aperture ultrasound imaging）を用いて、超音波画像の品質を改善することができる。「合成アパーチャ」とは、１つ以上のより小さい実際のアパーチャ（サブアパーチャ）を連続的に使用してＶＯＩを調べ、その位相中心が、特定及び／又は任意の形状の既知の１次元（１Ｄ）、２次元（２Ｄ）及び／又は、３次元（３Ｄ）経路に沿って動かされ、より大きな実効（仮想）アパーチャを実現して画像を取得する概念である。合成アパーチャは、電気音響トランスデューサ（例えば、トランスデューサアレイ）の空間位置を、連続的なビーム送信及び／又は受信位置に機械的に変更することによって、又は電気音響トランスデューサアレイ上の連続するビーム送信及び／又は受信位置の位相中心を電子的に変更することによって、或いはこれらの組合せによって形成することができる。合成アパーチャベースのイメージングは、当初、航空機によって上空から地上の関心領域を走査して、地上の広い領域をイメージングするレーダシステムで使用されたものである。超音波イメージングにおける合成アパーチャのフォーカシングは、超音波送信要素からＶＯＩの位置への及びこの位置から超音波受信素子に戻るまでの幾何学的な距離に基づく。超音波イメージングでは、合成アパーチャを使用することによって、全方向から、複数の送受信位置のそれぞれで記録された応答エコー（例えば、モノスタティック（mono-static）及びバイスタティック（bi-static）エコー）の受信された振幅及び位相データを解析することによって、ターゲット領域内の点へのフォーカシングが可能になり、これによって、全領域に関する情報を提供することができる。応答エコーの方向は、１つの受信チャンネルのみからでは判定することができないので、多くの受信チャンネルを用いて、応答エコーに含まれる情報を判定し、この情報をチャンネルの一部又は全部に亘って処理し、ターゲット領域の画像を生成するために使用される情報を最終的にレンダリングする。
【００１０】
合成アパーチャアレイは、１つ以上の実際のビームアパーチャサブアレイを含んでいてもよく、その位相中心は、図１Ｂに示すように、サンプリング位置間を移動する。例えば、トランスデューサアレイは、複数の実際のアパーチャサブアレイを含んでいてもよく、これらは、組み合わされて全体のアレイを構成し、サブアレイの１つ以上の位相中心は、（例えば、電子的及び／又は機械的に）サンプリング位置間を移動する。
【００１１】
合成アパーチャ超音波技術の１つの例では、単一又は複数のトランスデューサ素子を用いて、ＶＯＩを含む領域に亘る複数の位置において発散波面を送信し、画像領域全体をカバーする実効アパーチャを形成することができる。図１Ｂは、位相中心が連続的な位置を移動し、超音波イメージングにおける合成アパーチャビームを形成するトランスデューサアレイによって生成される超音波ビームを表している。この図に示すように、任意であるが既知の、開いた又は閉じた１Ｄ、２Ｄ又は３Ｄ経路Ｗに沿って、Ｎ個の連続的な位置（位置ａ、ｂ、ｃ…Ｎ）に位相中心が機械的及び／又は電子的に位置決めされる１つ以上のトランスデューサアレイ１１０によって複数の超音波１２１ａ、１２１ｂ…１２１Ｎが生成され、この経路は、例えば、直線、弧、円、螺旋又は定義された如何なる曲線経路等であってもよく、これにより、ＶＯＩ１２５の超音波イメージングのための合成アパーチャが形成される。実効アパーチャ内のトランスデューサ１１０のトランスデューサ素子の全部又は一部について受信されたモノスタティック及びバイスタティックエコー信号は、送信毎にサンプリングされる。例えば、個別の受信超音波信号データを用いて、個別のフォーカシングされていない送信に基づいて、低解像度画像を生成することができ、続いてこれを処理して、個別の受信超音波モノスタティック及びバイスタティックエコー信号から、２次元（２Ｄ）及び／又は３次元（３Ｄ）のフォーカシングされた画像を生成することができる。例えば、合成アパーチャ超音波イメージングを用いることによって、例えば、同等な性能の実際のビームシステムに比べて、超音波イメージングにおけるシステムの複雑性とコストを低減することができる。
【００１２】
音響パルスを生成するために用いられる波形のタイプも、超音波イメージングにおいて生成される画像の品質に影響を及ぼすことがある。従来の幾つかの超音波映像技術は、反射された信号から振幅情報のみを使用する。例えば、１つのパルスを出射し、反射信号を連続的にサンプリングしてもよい。生物組織では、音響速度は、略一定（例えば、骨を除いて、１０％未満の範囲内）であると考えられ、ここで、波形の出射から反射信号の受信までの時間は、波形が組織構造を伝播する距離（例えば、反射構造の深さ）に依存する。したがって、反射信号を複数の時間間隔でサンプリングし、複数の深さから反射する反射信号を受信することができる。また、異なる深さの異なる組織は、異なる量のエネルギで入射波形を部分的に反射することがあり、したがって、異なる媒体からの反射信号は、異なる振幅を有することがある。対応する超音波画像は、深さに基づいて生成することができる。したがって、新しい波形を出射するまでの時間は、画像に望まれる最大の深さに依存することができる。パルス式モノクロマチック及び／又は狭瞬時帯域幅波形を使用する超音波映像技術では、画像処理及び生成の解像度が低いという問題が生じることがある。一方、（例えば、周波数及び／又は位相によって）符号化されたスペクトル拡散、広瞬時帯域幅特徴を有する波形は、超音波イメージングのリアルタイム制御を可能にし、より高い品質の画像を生成することができる。
【００１３】
ここでは、合成アパーチャ超音波（synthetic aperture ultrasound：ＳＡＵ）イメージングにおいて使用されるコヒーレントなスペクトル拡散広瞬時帯域幅符号化波形（spread-spectrum, instantaneous-wideband, coded waveforms）を生成し、送信し、受信し及び処理する技術、システム及びデバイスを開示する。
【００１４】
ここに開示するＳＡＵイメージング技術は、既存の超音波映像技術に対して改善された画質、コントラスト及び解像度を提供することができ、及び組織の弁別及び分類を可能にすることができる。更に、ここに開示する技術において採用される例示的なスペクトル拡散広瞬時帯域幅符号化波形は、既存の医療超音波デバイスに存在するようなハードウェア設計上の制約を受けない。
【００１５】
ここに開示するＳＡＵ技術の機能のコヒーレント波形を使用することによって、エコー応答の一部又は全部を、選択された基準信号、例えば、送信波形に複素相関させることができる。このようなコヒーレントな複素相関によって、画像及び信号のアーチファクトを低減することができ、及び低い信号対雑音比及び干渉の存在下でデータを抽出することができる。
【００１６】
ここに開示するＳＡＵ技術の実現例におけるスペクトル拡散信号の使用によって、意図的で明確な振幅及び位相周波数コンテンツを有する音響波形の決定的な設計が可能になる。例えば、スペクトル拡散複合音響波形の各周波数成分の振幅及び／又は位相を明確に定義することによって、信号及び情報処理技術を用いて、例えば、数学的な限界に近い最大限の量の情報をエコー応答から抽出することができる。
【００１７】
ここに開示するＳＡＵ技術の実現例における瞬間的な、コヒーレントな、広帯域のスペクトル拡散符号化波形の使用によって、各送信−受信間隔の間に、利用可能な全ての情報を捕捉でき、これによって、例えば、生体の生物学的試料の不均質で動的な性質による、及び収集プロセスの動作によって誘導されるアーチファクトによる反射信号の劣化を最小化することができる。更に、ここに開示する技術に基づく信号及び情報処理方法を用いて、例えば、基本的な物理的パラメータ（例えば、バルク率、密度、減衰、音響インピーダンス、振幅反射、グループ遅延等）を抽出し、ＶＯＩの組織を弁別及び分化することができる。例えば、ここに開示するＳＡＵ技術の幾つかの信号及び情報処理方法は、ＶＯＩの組織を弁別及び／又は分類するための、受信した周波数及び角度依存広帯域、スペクトル拡散、合成アパーチャ受信信号エコーに対する数学的逆算法、並びに専門的システム技術、例えば、決定論的な、サポートベクトルネットワーク及びニューラルネットワーク技術を含む。
【００１８】
ここに開示するＳＡＵ技術の実現例における波形の各周波数成分の明確な振幅及び／又は位相符号化は、複数の利益を提供することができる。例えば、振幅符号化により、トランスデューサアレイ及び音響伝播チャンネルの周波数分散特性の明確な補償が可能になる。各周波数成分の振幅及び／又は位相符号化によって、広瞬時帯域幅波形の決定的なビーム整形及び操縦が可能になる。例示的な送信信号の各周波数成分の明確な振幅及び位相符号化によって、ピーク−平均パワー比（peak-to-average power ratio：ＰＡＰＲ）を最小化でき及び音響パワーを広帯域に広げることができ、例えば、有害な生物学的影響を最小にすることができる。例えば、スペクトル拡散信号の各周波数成分の振幅及び／又は位相を明確に定義することによって、同時に送信でき、互いの干渉が最小化された複数の波形を構築することができ、これによって、信号及び情報処理技術を用いて、個別の送信波形に関連する受信信号を復元することができる。更に、ここに開示するＳＡＵ技術の符号化されたスペクトル拡散音響波形によって、これらの波形の特定の曖昧性（ambiguity properties）のために、動き補償を行うことができる。
【００１９】
一側面として、音響イメージングデバイスにおいて音響波形を生成する方法を開示する。この方法は、１つ以上の波形シンセサイザにおいて、音響イメージングデバイスのトランスデューサアレイの１つ以上の空間位置及び／又はトランスデューサアレイの１つ以上のビーム位相中心位置からターゲットに向けて送信される複数の個別の符号化波形から形成される１つ以上の複合波形を合成することを含む。複合波形の個別の符号化波形は、互い直交し、異なる周波数帯にあり、個別の相互に直交する符号化波形のそれぞれは、対応する位相による固有の周波数を有する。この方法は、ターゲットに対する１つ以上の送信位置から、複数の音響波形からなる１つ以上の複合音響波形を送信することを含み、送信は、アレイの１つ以上のトランスデューサ素子を選択し、１つ以上の複合音響波形のそれぞれの複数の個別の直交符号化波形を、１つ以上の複合音響波形のそれぞれの複数の対応する音響波形にトランスデュースすることを含む。更に、この方法は、ターゲットに対する１つ以上の受信位置において、送信される音響波形に対応するターゲットの少なくとも一部から戻る応答音響波形を受信することを含み、受信は、アレイのトランスデューサ素子の少なくとも幾つかを選択して、応答音響波形を受信することを含む。受信位置は、ターゲットに対するトランスデューサ素子のアレイの空間位置と、応答音響波形を受信するアレイのビーム位相中心位置との一方又は両方を含む。送信される音響波形及び応答音響波形は、音響イメージングデバイスの拡大された実効アパーチャを生成する。
【００２０】
幾つかの具体例では、方法は、例えば、音響波形をターゲットに送信する際、アレイのトランスデューサ素子を制御して、ターゲットに対する複合波形の向きを変更し、ターゲットが、イメージング期間に亘って、異なる波形の向きで音響波形を受信するようにする。例えば、ターゲットに対する複合波形の向きの変更は、同じ又は異なる空間位置から異なる複合波形を送信すること、異なる空間位置から同じ又は異なる複合波形を送信すること、及びトランスデューサ素子のアレイ上の異なるビーム位相中心位置から同じ又は異なる複合波形を送信することを含む。
【００２１】
幾つかの具体例においては、方法は、例えば、受信した応答音響波形を、１つ以上の複合波形に対応する１つ以上の受信複合波形として、アナログフォーマットからデジタルフォーマットに変換することを更に含み、各受信複合波形は、ターゲットの情報を含み、情報は、受信複合波形の対応する周波数帯に関連する振幅及び位相を含む。また、幾つかの具体例においては、方法は、例えば、受信した（１つ以上の受信複合波形の）応答音響波形を処理して、ターゲットの少なくとも一部の画像（例えば、２Ｄ及び／又は３Ｄ画像）を生成することを更に含む。
【００２２】
一側面として、合成アパーチャ超音波画像処理システムを開示する。システムは、波形発生器に接続された１つ以上の波形シンセサイザを有する波形生成ユニットを含む。波形生成ユニットは、波形発生器が提供する波形情報に応じて１つ以上の波形シンセサイザによって生成される異なる周波数帯に対応する複数の個別の直交符号化波形からなる複合波形を合成し、個別の直交符号化波形は、それぞれ互いに直交し及び異なる周波数帯に対応し、各個別の直交符号化波形は、対応する位相による固有の周波数を有する。システムは、送信モードと受信モードの間で切り換えられる送信／受信スイッチングユニットを含む。また、システムは、送信／受信スイッチングユニットと通信するトランスデューサ素子のアレイを含む。トランスデューサ素子のアレイは、ターゲットに対する１つ以上の送信位置から複数の音響波形からなる複合音響波形を送信するように構成されており、複合音響波形の送信される音響波形は、複合波形の合成された個別の直交符号化波形に基づいている。更に、トランスデューサ素子のアレイは、例えば、ターゲットに対する１つ以上の受信位置において、ターゲットの少なくとも一部から戻る複数の送信された音響波形に対応する応答音響波形を受信するように構成されている。送信される音響波形及び応答音響波形は、合成アパーチャ音響波形イメージングシステムの拡大された実効アパーチャを生成する。それぞれの波形を送信及び受信するための送信位置及び受信位置は、それぞれ、ターゲットに対するトランスデューサ素子のアレイの空間位置及び音響波形を送信及び／又は受信するためのアレイの上のビーム位相中心位置の一方又は両方を含む。システムは、更に、トランスデューサ素子のアレイと通信して、複数の個別の直交符号化波形を複数の対応する音響波形にトランスデュースするアレイの１つ以上のトランスデューサ素子を選択し、及び応答音響波形を受信するアレイの１つ以上のトランスデューサ素子を選択する多重化ユニットを含む。幾つかの実現例では、システムは、例えば、ターゲットの情報を提供する、トランスデューサ素子のアレイによって受信された応答音響波形を、アナログフォーマットからデジタルフォーマットに変換するアナログ／デジタル（Ａ／Ｄ）変換器のアレイを含む。また、システムは、波形生成ユニット及びトランスデューサ素子のアレイ（例えば、Ａ／Ｄコンバータのアレイを介してもよい。）と通信するコントローラユニットを含み、コントローラユニットは、データを保存するメモリユニットと、メモリユニットに接続されて、情報をデータとして処理する処理ユニットとを含む。システムは、コントローラユニットと通信するユーザインタフェースユニットを含むことができる。システムの幾つかの具体例では、例えば、コントローラユニットは、処理されたデータから、ターゲットの少なくとも一部の画像を生成するように構成される。
【００２３】
図１Ｃ及び図１Ｄは、複数の送信位置から合成アパーチャビームを形成するトランスデューサアレイによって生成される例示的な複合超音波ビームを示している。図に示すように、トランスデューサアレイ１３０は、複数の実際のアパーチャサブアレイＳｕｂ１、Ｓｕｂ２…ＳｕｂＮを含み、各サブアレイは、個別のトランスデューサ素子（例えば、１〜１６個、３２個、６４個等の素子）を含む。トランスデューサ素子の一部又は全部は、超音波イメージングのための合成アパーチャを形成する１以上の個別の相互に直交する符号化された音響波形１３１ａ、１３１ｂ…１３１ｎの複合音響波形を（例えば、順次的に、同時に又はランダムに）送信でき、この複合音響波形は、複数のサブアレイ位相中心位置からターゲット関心ボリューム（volume of interest：ＶＯＩ）１３５に送信され、アレイ１３０を形成する。また、アレイ１３０を形成するトランスデューサ素子の一部又は全部は、（個別の相互に直交する符号化された音響波形１３１ａ、１３１ｂ…１３１ｎに基づいて形成される）送信される音響波形に対応する応答音響波形を受信し、ここで、受信される音響波形は、ＶＯＩ１３５の少なくとも一部から（例えば、反射、屈折、回折、遅延及び／又は減衰されて）後方散乱して戻る。これによって、受信される個別の音響波形は、送信された複合音響波形に対応する１つ以上の受信複合波形を形成する。複合音響波形は、個別の音響波形を生成するために用いられる複数のスペクトル拡散広瞬時帯域幅符号化波形から形成される複合的な合成波形に基づいて生成される。個別の複合音響波形１３１ａ、１３１ｂ…１３１ｎは、トランスデューサアレイ１３０のサブアレイの１つ以上によってトランスデュースすることができる。トランスデューサアレイ１３０は、既知の経路に沿った複数の物理的位置（例えば、図１Ｃに示す位置１から位置２）及び／又は複数のビーム操縦位置（例えば、図１Ｄに示す位置１から位置２）に配置することができ、位相中心は、例えば、合成アパーチャを形成する連続的な位置に機械的に、電子的に、又は機械的且つ電子的に位置決めされる。
【００２４】
幾つかの実現例では、例えば、１つ以上のサブアレイに亘るトランスデューサ素子の様々な組合せから構成される様々なアレイセグメントを用いて、直交符号化音響波形１３１ａ、１３１ｂ…１３１ｎを生成することができる。図１Ｅは、ここに開示する技術の複合超音波ビームを生成するための例示的なトランスデューサアレイを示している。図１Ｅの例では、トランスデューサアレイ１３０は、４つのトランスデューサセグメント（例えば、トランスデューサセグメント１、２、３、４）に構成された６４個の個別のトランスデューサ素子を含む。この例では、（例えば、４つのトランスデューサセグメントの１つ以上に亘るトランスデューサ素子を含む）６４個の個別のトランスデューサ素子の何れかを含む１つ以上のサブアレイは、個別の直交する符号化された音響波形１３１ａ、１３１ｂ  …１３１ｎを（例えば、順次的に、同時に又はランダムに）送信することができる。サブアレイは、１つのトランスデューサセグメント内の又は複数のトランスデューサセグメントに亘る個別のトランスデューサ素子の組合せを含むことができる。例えば、サブアレイ１は、トランスデューサセグメント１のトランスデューサ素子２、３、６、７、１０、１１、１４、１５及びトランスデューサセグメント３のトランスデューサ素子２、３を含み、サブアレイ２は、トランスデューサセグメント２のトランスデューサ素子１、２、３、５、６、７、９、１０、１１を含み、サブアレイ３は、トランスデューサセグメント３のトランスデューサ素子９、１０、１１、１２、１３、１４、１５、１６及びトランスデューサセグメント４のトランスデューサ素子９、１３を含み、サブアレイ４は、トランスデューサセグメント４のトランスデューサ素子５、６、７、９、１０、１１、１３、１４、１５を含む。図２Ａ及び図２Ｂを用いて後述するように、サブアレイの構成は、波形発生器とトランスデューサアレイの間に接続されたスイッチング素子（例えば、マルチプレクサユニット）を用いて生成することができる。
【００２５】
図２Ａは、拡張された実効（合成）アパーチャに亘って強化された波形特性を有する音響波形を生成できる例示的な合成アパーチャ超音波（Synthetic Aperture Ultrasound：ＳＡＵ）システム２００のブロック図である。ＳＡＵシステム２００によって生成される複合音響波形の強化された波形特性は、スペクトル拡散、広瞬時帯域幅、コヒーレント、擬似ランダム雑音特性及び符号化を含む。ＳＡＵシステム２００は、多くのシステム設計の１つによって構成できる。ある例では、ＳＡＵシステム２００は、時間的同期のためのマスタクロック２０１を含むことができる。マスタクロック２０１は、互いに同期されて動作するシステムコントローラ２０２並びにＳＡＵシステム２００の他のモジュールにインタフェースすることができる。システムコントローラ２０２は、処理ユニット、例えば、ＲＩＳＣベース又は他のタイプのＣＰＵアーキテクチャの中央演算処理装置（central processing unit：ＣＰＵ）を含むことができる。システムコントローラ２０２は、少なくとも１つの入出力（Ｉ／Ｏ）ユニット及び／又はメモリユニットを含んでいてもよく、これらは、処理ユニットと通信して、システムコントローラ２０２の様々な機能をサポートする。例えば、処理ユニットは、システム制御バス、例えば、データ及び制御バス２０３に接続することができる。システムコントローラ２０２は、例えば、パーソナルコンピュータ（personal computer：ＰＣ）、ラップトップコンピュータ、タブレット及びモバイル通信デバイスアーキテクチャ等の様々なデータ処理アーキテクチャの１つとして実現できる。
【００２６】
メモリユニットは、命令、ソフトウェア、値、画像、及び処理ユニットによって処理又は参照される他のデータ等の他の情報及びデータを保存できる。メモリユニットのストレージ機能を実現するために、ランダムアクセスメモリ（Random Access Memory：ＲＡＭ）デバイス、読取専用メモリ（Read Only Memory：ＲＯＭ）デバイス、フラッシュメモリデバイス及び他の適切な様々なタイプのストレージ媒体を用いることできる。メモリユニットは、事前記憶波形（pre-stored waveforms）、係数データ及び情報を保存でき、これらを用いて、例えば、スペクトル拡散、広瞬時帯域幅、コヒーレント、擬似ランダム雑音、周波数符号化波形及び／又は位相符号化波形等の波形を生成することができる。メモリユニットは、受信及び処理された波形から得られたデータ及び情報を保存でき、これを用いて、新たな波形を生成及び送信することができる。メモリユニットは、システム制御バス、例えば、データ及び制御バス２０３に接続してもよい。
【００２７】
入出力ユニットは、外部インタフェース、データストレージのソース及び／又はディスプレイデバイスに接続してもよい。入出力ユニットは、システム制御バス、例えば、データ及び制御バス２０３に接続してもよい。以下に限定されるわけではないが、汎用シリアルバス（Universal Serial Bus：ＵＳＢ）、ＩＥＥＥ１３９４（ファイアワイヤ：FireWire）、ＩＥＥＥ８０２．１１１（ブルートゥース：Bluetooth(登録商標)）、無線ローカルエリアネットワーク（Wireless Local Area Network：ＷＬＡＮ）、無線パーソナルエリアネットワーク（Wireless Personal Area Network：ＷＰＡＮ）、無線ワイドエリアネットワーク（Wireless Wide Area Network：ＷＷＡＮ）、ＩＥＥＥ８０２．１６（Worldwide Interoperability for Microwave Access：ＷｉＭＡＸ）及びパラレルインタフェース等の一般的なデータ通信規格に互換性がある様々なタイプの有線又は無線インタフェースを用いて入出力ユニットを実現することができる。入出力ユニットは、外部インタフェース、データストレージのソース又はディスプレイデバイスに接続し、データ及び情報を受信又は出力してもよく、これらのデータ及び情報は、プロセッサユニットで処理し、メモリユニットに保存し、又は出力ユニットにおいて表示してもよい。
【００２８】
システムコントローラ２０２は、例えば、データ及び制御バス２０３への接続を介して、ＳＡＵシステム２００の全てのモジュールを制御できる。例えば、データ及び制御バス２０３は、システムコントローラ２０２を１つ以上の付属デジタルシグナルプロセッサ、例えば、デジタルシグナルプロセッサ２０４に接続し、その機能を制御して、波形を処理してもよい。デジタルシグナルプロセッサ２０４は、１つ以上のプロセッサ、例えば、以下に限定されるものではないが、特定用途向けＩＣ（application-specific integrated circuit：ＡＳＩＣ）、フィールドプログラマブルゲートアレイ（field-programmable gate array：ＦＰＧＡ）、デジタルシグナルプロセッサ（digital signal processors：ＤＳＰ）、ＡｓＡＰ（asynchronous array of simple processors）及び他のタイプのデータ処理アーキテクチャを含んでいてもよい。また、データ及び制御バス２０３は、システムコントローラ２０２及びデジタルシグナルプロセッサ２０４をユーザインタフェースのためのモジュールを有する１つ以上のディスプレイユニット、例えば、モジュールユーザインタフェース２０６を有するディスプレイ２０５に接続し、ユーザ又はオペレータに情報を提供し、及びユーザ又はオペレータから入力／命令を受け取ってもよい。ディスプレイ２０５は、視覚的ディスプレイとして、多くの適切なディスプレイユニット、例えば、以下に限定されるものではないが、陰極線管（cathode ray tube：ＣＲＴ）、発光ダイオード（light emitting diode：ＬＥＤ）、液晶ディスプレイ（liquid crystal display：ＬＣＤ）モニタ及び／又はスクリーンを含むことができる。また、ディスプレイ２０５は、様々なタイプのディスプレイ、スピーカ又はプリントインタフェースを含むことができる。他の具体例では、ディスプレイ２０５は、他のアウトプット装置、例えば、トナー、液体インクジェット、固体インク、昇華式又はインクレス（感熱式又は紫外線式）プリンタ、及び様々なタイプのオーディオ信号再生装置を含んでいてもよい。ユーザインタフェース２０６は、適切な多くのインタフェース、例えば、様々なタイプのキーボード、マウス、音声コマンド、タッチパッド、ブレインマシンインタフェース（brain-machine interface）装置を含むことができる。
【００２９】
ＳＡＵシステム２００は、システムコントローラ２０２によって制御されて、１つ以上のデジタル波形を生成する波形発生器２０７を備えていてもよい。波形シンセサイザ及びビームコントローラのアレイの少なくとも１つの素子、例えば、この具体例では、波形シンセサイザ及びビームコントローラ２０８として示されている素子によって、１つ以上のデジタル波形をアナログ電気信号（例えば、アナログ波形）として生成することができる。波形発生器２０７は、関数発生器及び任意波形発生器（arbitrary waveform generator：ＡＷＧ）の少なくとも１つであってもよい。例えば、波形発生器２０７をＡＷＧとして構成し、任意のデジタル波形を生成し、これを波形シンセサイザ及びビームコントローラ２０８によって、個々のアナログ波形及び／又は複合アナログ波形として合成してもよい。また、波形発生器２０７は、デジタル波形の生成に使用される事前記憶波形、係数データ及び情報を保存する少なくとも１つのメモリユニットを含むことができる。
【００３０】
図２Ａに示すＳＡＵシステム２００は、Ｉ個のアレイ素子を含む波形シンセサイザ及びビームコントローラ２０８を備える。一具体例では、波形シンセサイザ及びビームコントローラ２０８は、Ｉ個のアレイ式波形シンセサイザの各ラインに少なくとも１つの波形シンセサイザ要素を含むように構成できる。他の具体例では、波形シンセサイザ及びビームコントローラ２０８は、Ｉ個のアレイ式ビームコントローラの各ラインに少なくとも１つのビームコントローラ要素を含むことができる。他の具体例では、波形シンセサイザ及びビームコントローラ２０８は、Ｉ個のアレイ式波形シンセサイザ及びビームコントローラの各ラインに少なくとも１つの波形シンセサイザ要素及びビームコントローラ要素を含むことができる。波形シンセサイザ及びビームコントローラ２０８は、電気信号、例えば、無線周波数（radio frequency：ＲＦ）波形を生成するための位相ロックループシステムを含んでいてもよい。例示的なＲＦ波形は、波形シンセサイザ及びビームコントローラ２０８のアレイ素子によって生成された個々の波形から、波形シンセサイザ及びビームコントローラ２０８によって合成され、例えば、個々のＲＦ波形は、波形シンセサイザ及びビームコントローラ２０８の他のアレイ素子によって生成される他の全ての個々の波形と実質的に同時に１つのアレイ素子によって生成してもよい。個々のＲＦ波形のそれぞれは、周波数成分又は「チップ」とも呼ばれる特定の周波数帯毎に定義され、個々の波形のそれぞれの波形特性は、波形発生器２０７によって決定され、チップに対応する少なくとも１つの振幅値及び少なくとも１つの位相値を含んでいてもよい。波形発生器２０７は、波形シンセサイザ及びビームコントローラ２０８に命令を発し、個々の波形のそれぞれの特性に関する情報を含む波形データを波形シンセサイザ及びビームコントローラ２０８に供給し、個々の直交するＲＦ波形を生成し、波形シンセサイザ及びビームコントローラ２０８は、これらを合成して、複合ＲＦ波形を生成することができる。
【００３１】
波形シンセサイザ及びビームコントローラ２０８によって生成される個々のＲＦ波形及び／又は複合ＲＦ波形は、Ｉ個の増幅器のアレイを含む出力増幅器２０９によって、例えば、波形を増幅し及び／又は波形の位相をシフトさせることによって変更される。幾つかの具体例では、出力増幅器２０９は、線形アンプとして構成される。出力増幅器２０９は、トランスデューサドライバとして用いることができる。個々のＲＦ波形及び／又は複合ＲＦ波形は、例えば、Ｎ極双投送信／受信スイッチ（N-pole double-throw transmit/receive switch）である送信／受信（Transmit/Receive：Ｔ／Ｒ）スイッチ２１０に供給される。Ｔ／Ｒスイッチ２１０は、ＳＡＵシステム２００のトランスデューサモジュールに接続してもよい。Ｔ／Ｒスイッチ２１０は、例えば、Ｎ極多重化スイッチを含むことによって、多重化ユニットとして動作できる。生成されたＲＦ波形、例えば、ターゲット媒体に送信される複合ＲＦ波形及び／又は少なくとも１つの個々のＲＦ波形は、トランスデューサモジュールによって、例えば、音響波にトランスデュースできる。図２Ａに示す例では、トランスデューサモジュールは、トランスデューサ素子のアレイ、例えば、Ｘ個のトランスデューサ素子からなるトランスデューサアレイ２１１として構成される。例えば、トランスデュースされた音響波は、音響波形バーストの形式で出射することができる。例えば、トランスデューサアレイ２１１の選択されたアレイ素子は、波形発生器２０７によって決定され、空間的に組み合われて複合音響波形を形成する個別の直交波形に対応する２つ以上の個別の直交音響波形を生成（例えば、トランスデュース（変換））してもよい。更なる例として、トランスデューサアレイ２１１の選択されたアレイ素子は、波形発生器２０７によって決定された複合波形に対応する１つ以上の複合音響波形を生成（例えば、トランスデュース（変換））してもよい。
【００３２】
例えば、Ｔ／Ｒスイッチ２１０が送信モードに構成されている場合、例示的な、トランスデュースされ、送信された複合音響波形は、ターゲット、例えば、生物組織に対するトランスデューサアレイ２１１の複数の位置からターゲット領域に向けて送信することができ、このトランスデュースされ、送信された音響波形は、空間的に結合された音響波形を形成する。送信された複合音響波形は、ターゲット媒体内に伝播し、ターゲット媒体は、例えば、送信された音響波形の一部を透過し、一部を反射する１つ以上の不均質な媒体を有することができる。例えば、音響波形が送信された後、Ｔ／Ｒスイッチ２１０は、受信モードに構成することができる。ターゲットによって（少なくとも部分的に）反射された例示的な複合音響波形は、トランスデューサアレイ２１１によって受信することができ、これは、応答音響波形と呼ばれる。幾つかの具体例では、トランスデューサアレイ２１１のＸ個のアレイ素子の選択されたアレイ素子は、個別の直交波形（例えば、周波数チップ）に対応する応答音響波形を受信し、これをアナログＲＦ波形に変換するように構成することができる。他の例では、トランスデューサアレイ２１１の選択されたトランスデューサ素子は、例えば、例示的なトランスデューサアレイ２１１の制御ロジックと通信するシステムコントローラ２０２によって判定される選択制御信号に基づいて、送信された複合波形に対応する応答音響波形を受信し、これを複合アナログＲＦ波形に変換するように構成してもよい。
【００３３】
幾つかの実現例では、例えば、トランスデューサアレイ２１１は、ＳＡＵシステム２００を用いた超音波画像化を行っている間、トランスデューサアレイ２１１を空間的に移動させ、合成アパーチャを生成することによって、データサンプリング／超音波スキャン位置の１次元、２次元及び／又は３次元に沿って、ビーム位相中心を機械的に平行移動させることができる。更なる例では、トランスデューサアレイ２１１は、静止していてもよく、ＳＡＵシステム２００を用いた超音波画像化を行っている間に、システムコントローラ２０２によって、データサンプリング／超音波スキャン位置として、Ｘ個のトランスデューサ素子の部分を順次的に又はランダムにアドレス指定し、合成アパーチャを生成することによって、静止したトランスデューサアレイ２１１に沿って、ビーム位相中心を１次元、２次元及び／又は３次元に電子的に平行移動させてもよい。更なる例として、超音波画像化を行っている間に、ＳＡＵシステム２００は、位相中心を、データサンプリング／超音波スキャン位置の１次元、２次元及び／又は３次元に機械的及び電子的に平行移動させて合成アパーチャを生成してもよい。
【００３４】
幾つかの実現例では、例えば、個別のＲＦ波形及び／又は複合ＲＦ波形を、例えば、時間的に個別に又は同時に、別々の送信トランスデューサアレイに直接的に渡し、別々の受信トランスデューサアレイを用いて、関心ボリュームから散乱した音響波形を受信することができる。図２Ｂは、ＳＡＵシステム２００ａを表すブロック図として、ここに開示する技術の合成アパーチャ超音波システムの他の例示的な実施形態を示している。ＳＡＵシステム２００ａは、ＳＡＵシステム２００と同様に動作して、拡張された実効（合成）アパーチャに亘って、改良された波形特性を有する音響波形を生成するものであり、（１）（例えば、波形発生器２０７によって決定される個別の直交符号化波形に対応する）Ｉ個の個別の波形をトランスデュースし、例えば、順次的に、又は同時に、又は時間的にランダムに１つ以上の音響波形を送信するＸ個のトランスデューサ素子からなる送信トランスデューサアレイ２１１ａを含み、更に、（２）Ｙ個のトランスデューサ素子からなる別の受信トランスデューサアレイ２１１ｂを含み、これは、応答音響波形を受信し、それぞれをアナログＲＦ波形に変換するように構成することができる。幾つかの具体例では、Ｘは、Ｙに等しく、他の具体例では、Ｘは、Ｙに等しくない。幾つかの実現例では、ＳＡＵシステム２００ａは、出力増幅器２０９と送信トランスデューサアレイ２１１ａの間に配設されたＴ／Ｒスイッチ２１０ａを含むことができ、Ｔ／Ｒスイッチ２１０ａは、例えば、Ｎ極多重化スイッチを含むことによって、多重化ユニットとして動作することができる。また、幾つかの実現例では、ＳＡＵシステム２００ａは、受信トランスデューサアレイ２１１ｂと前置増幅器モジュール２１２の間に配置されたＴ／Ｒスイッチ２１０ｂを含むことができ、Ｔ／Ｒスイッチ２１０ｂは、例えば、Ｎ極多重化スイッチを含むことによって、多重化ユニットとして動作することができる。例えば、送信トランスデューサアレイ２１１ａは、Ｉ個の個別の直交符号化波形を、波形発生器２０７によって決定されるＲＦ波形に対応する１つ以上の音響波形にトランスデュースすることができる。トランスデュースされた音響波形は、音響波形バーストの形式で出射してもよい。例えば、トランスデューサアレイ２１１ａの選択されたアレイ素子は、波形発生器２０７によって決定される個別の直交波形に対応する２以上の個別の直交音響波形を生成（例えば、トランスデュース）してもよく、これらを空間的に結合して、例えば、順次的に、又は同時に、又は時間的にランダムに送信される１つ以上の音響波形である、１以上の複合音響波形を形成してもよい。更なる例として、トランスデューサアレイ２１１ａの選択されたアレイ素子は、波形発生器２０７によって決定される複合波形に対応する１つ以上の複合音響波形を生成（例えば、トランスデュース）してもよく、これらは、順次的に、又は同時に、又は時間的にランダムに送信される。
【００３５】
例示的なＳＡＵシステム２００ａでは、送信トランスデューサアレイ２１１ａ及び受信トランスデューサアレイ２１１ｂは、ＳＡＵシステム２００ａを用いた超音波画像化を行っている間に、データサンプリング／超音波スキャン位置の１次元、２次元及び／又は３次元に、位相中心の一部又は全部を機械的、電子的に又は機械的且つ電子的に平行移動させて合成アパーチャを生成することができる。幾つかの実現例では、例えば、送信トランスデューサアレイ２１１ａのＸ個のトランスデューサ素子は、位相中心の一部又は全部を１次元、２次元及び／又は３次元に機械的に、電子的に、又は機械的且つ電子的に一斉に平行移動させ、他の実現例では、例えば、送信トランスデューサアレイ２１１ａのＸ個のトランスデューサ素子の位相中心の１つ以上は、アレイ２１１ａの他の素子から独立して、１次元、２次元及び／又は３次元に機械的に、電子的に、又は機械的且つ電子的に平行移動させてもよい。幾つかの実現例では、例えば、トランスデューサアレイ２１１ａ、２１１ｂのＸ個又はＹ個のトランスデューサ素子は、それぞれ、放射音響ビーム及び受信ビームを、１つ又は２つの角度次元（例えば、方位角及び／又は高度角）で、（例えば、機械的に、電子的に、又は機械的且つ電子的に）走査することができる。同様に、例えば、幾つかの実現例では、受信トランスデューサアレイ２１１ｂのＹ個のトランスデューサ素子の位相中心は、（例えば、機械的に、電子的に、又は機械的且つ電子的に）１次元、２次元及び／又は３次元に一斉に平行移動させてもよく、他の実現例では、受信トランスデューサアレイ２１１ｂのＹ個のトランスデューサ素子の１つ以上は、アレイ２１１ｂの他の素子から独立して、１次元、２次元及び／又は３次元に平行移動させてもよい。例えば、一実施形態においては、トランスデューサアレイ２１１又は２１１ａのＸ個のトランスデューサ素子のそれぞれは、Ｉ個の出力増幅器２０９の１つに対応し、例えば、Ｘ＝Ｉである。これに代えて、例えば、構成要素の数を減らすために、全てのＸ個のトランスデューサ素子から送信トランスデューサ素子の幾つかグループを構成し、これらを互いに多重化し、Ｎ極多重化スイッチを介して、Ｉ個より少ない増幅器２０９と通信させてもよい。
【００３６】
図２Ａ及び図２Ｂに示すように、個別の受信（アナログ）ＲＦ波形は、前置増幅器モジュール２１２によって修正してもよく、前置増幅器モジュール２１２は、Ｊ個以下の数の増幅器を含み、幾つかの例では、例えば、アレイ内の各増幅器は、それぞれ、トランスデューサアレイ２１１又は受信トランスデューサアレイ２１１ｂのＸ又はＹ個のトランスデューサ素子に対応していてもよい。これに代えて、例えば、構成要素の数を減らすために、全ての受信トランスデューサ素子から受信トランスデューサ素子の幾つかグループを構成し、これらを互いに多重化し、Ｎ極多重化スイッチを介して、前置増幅器と通信させてもよい。例えば、前置増幅器モジュール２１２の対応する増幅器のそれぞれ又は少なくとも幾つかの増幅器は、個別の受信（アナログ）ＲＦ波形を増幅し及び／又は位相をシフトさせることができる。幾つかの具体例では、前置増幅器モジュール２１２のＪ個の増幅器のアレイは、線形増幅器として構成される。幾つかの具体例では、前置増幅器モジュール２１２は、受信波形に対して他の信号調整技術を実行するように実現してもよい。増幅及び／又は信号調整の後、例えば、Ｊ個のＡ／Ｄコンバータのアレイを含むデジタル−アナログ（Ａ／Ｄ）変換器モジュール２１３によって、個別の受信波形をアナログフォーマットからデジタルフォーマットに変換できる。Ａ／Ｄ変換器モジュール２１３は、最下位ビット（least significant bit：ＬＳＢ）ジッタ、スプリアスフリーダイナミックレンジ（spurious-free dynamic range：ＳＦＤＲ）及び波形依存性を有し、例示的な波形を適切にデコードすることができるＡ／Ｄコンバータを含むことができる。個別の受信波形の変換されたデジタル表現は、ターゲット媒体の代表的な画像を生成及び形成するように、プロセッサ、例えば、デジタルシグナルプロセッサ２０４によって処理することができる。
【００３７】
ＳＡＵシステム２００は、多くの動作モードの１つで動作することができる。ある例では、マスタクロック２０１は、例えば、波形発生器２０７、波形シンセサイザ２０８及びＤＳＰ２０４を含むＳＡＵシステム２００のモジュールを同期させるためのタイムベースを提供することができる。マスタクロック２０１は、例示的な波形の位相符号化を可能にするように、低位相雑音クロックとして構成することができる。オペレータは、ユーザインタフェース２０６によって、合成アパーチャモードの動作を選択することができる。ＳＡＵシステム２００の幾つかの実現例では、合成アパーチャモードは、以下に限定されるわけではないが、合成アパーチャストリップ走査（synthetic aperture strip scan：ＳＡＳＳ）モード及び例示的な合成アパーチャスポットライト（synthetic aperture spotlight：ＳＡＳｐｌ）モードを含むことができる。
【００３８】
図２Ｃ及び図２Ｄは、それぞれ、ここに開示する技術の例示的なＳＡＳＳモード及び例示的なＳＡＳｐｌモードを示している。図２Ｃに示すように、例えば、ＳＡＳＳモードは、トランスデューサアレイ２１１の位相中心の電子的、機械的、又は電子的且つ機械的な１次元、２次元及び／又は３次元の平行移動によって、ターゲットの周囲の複数の位置から送信される複数の音響波形に基づいて、合成アパーチャを生成するように動作可能である。図２Ｄに示すように、例えば、ＳＡＳｐｌモードは、送信及び／又は受信ビームを操作してＶＯＩに中心合わせするように動作可能であり、したがって、例えば、ＳＡＳＳモードより大きな合成アパーチャを生成し、これによって、横方向の画像解像度を高めることができる。ＳＡＳＳ及びＳＡＳｐｌの例示的な動作モードは、２Ｄ平面スライス及び／又は３Ｄ立体断層撮影レンダリングを含み、オペレータは、例えば、ＳＡＵシステム２００を使用して、これを操作し、矢状面、冠状面、横断面又は任意のあらゆる平面を提供することができる。オペレータは、実際のアパーチャの動作モードを選択することができる。
【００３９】
例えば、ユーザに対して用意され、ユーザインタフェース２０６で選択できる幾つかの例示的な動作モードは、従来のＡ−モード（例えば、１Ｄ深さのみの画像）、従来のＢ−モード（例えば、２Ｄ平面画像−横断方向対深さ方向）、従来のＣ−モード（例えば、選択された深さにおける２Ｄ平面画像）、従来のＤ−モード（例えば、ドップラモード）、及び従来の又は新しい動作モードの何れか１つ以上と組み合わされた統合型外科的治療モードとしての高強度超音波（High Intensity Focused Ultrasound：ＨＩＦＵ）を含む。例示的なドップラモードは、カラードップラ（例えば、色分けされたドップラ及びＢ−モード画像の重ね合せ）、連続ドップラ（例えば、１Ｄドップラプロファイル対深さ）、パルス波ドップラ（例えば、選択されたボリュームのためのドップラ対時間）、及び２重／３重ドップラ（例えば、従来のＢ−モード、従来のＣ−モード又は色ドップラ及びパルス波ドップラの重ね合せ）を含む。例示的なＳＡＵシステム２００は、更に、スペクトル拡散、広瞬時帯域幅、周波数及び／又は位相符号化波形を生成できる新しい動作モードを実行することができる。例えば、ユーザは、従来のＢ−モードに類似するが、画質が大幅に向上した（例えば、解像度、コントラスト比等がより高い）高解像度２Ｄ画像モードを選択することができ、又はユーザは、ＣＴ及び磁気共鳴映像法（Magnetic Resonance Imaging：ＭＲＩ）と同様の手法でユーザが選択可能な２Ｄ画像として表示することができる立体画像を生成する高解像度３Ｄイメージングモードを選択することができる。更に、例えば、ユーザは、微視的組織学的研究のための組織染色に類似する組織弁別を補助するための画像色分けが、Ｂ−モード、Ｃ−モード、Ｄ−モード又は他のモードと組み合わされた例示的なＡＴＳモード（人工組織染色（Artificial Tissue Staining）モード）及び組織タイプを弁別及び特定する例示的なＣＡＤモード（コンピュータ支援診断（Computer Aided Diagnostic）モード）を選択することができる。ＡＴＳモードは、ターゲット領域からの応答エコー波形、例示的な送信されたスペクトル拡散、広瞬時帯域幅、符号化音響波形からの、例えば、応答エコーから取得される複数の測定された特性の１つ以上に基づいて、画像処理において画像色分けのための特徴を採用することができる。ＣＡＤモードは、分類器（アルゴリズム）を用いて、例えば、ターゲット領域からの応答エコー、例えば、例示的なスペクトル拡散、広瞬時帯域幅、符号化、角度ダイバーシチ、モノスタティック及びバイスタティック、合成アパーチャ音響波形からの応答エコーの測定された特性の特徴に基づいて、組織タイプを分類することができる。特徴的特性は、異なるインピーダンス、（波長の関数としての）波形反射、グループ遅延等を含むことができる。ＣＡＤモードで採用できる幾つかの例示的な分類器は、決定論的分類器、確率的分類器（例えば、ベイズ及びサポートベクトルネットワーク分類器）、及びニューラルネットワーク分類器等を含むことができる。
【００４０】
図２Ｅは、合成アパーチャ超音波イメージングのためにＳＡＵシステム２００又は２００ａを動作させる例示的な動作方法２５０を表している。タイムエポック毎に、方法２５０は、プロセス２５１を実行して、ユーザ定義の動作モードをチェックすることによって開始できる。例えば、動作モードは、ユーザがユーザインタフェース２０６を用いて選択してもよく、他のエンティティによって選択してもよく、ＳＡＵシステム２００又は２００ａの内部で選択してもよい。選択された動作モードに基づき、システムコントローラ２０２は、例えば、プロセス２５２において、波形発生器２０７が、波形シンセサイザ及びビームコントローラ２０８の１つ以上の要素に対し、望ましい広帯域複合ＲＦ波形を構成する各周波数チップの周波数、振幅及び位相を定義するデジタルメッセージ（データ）を発行するよう命令する。プロセス２５１は、方法２５０の間のどこで行ってもよく、複数のインスタンスで実現してもよい。方法２５０は、波形シンセサイザ及びビーム整形器、例えば、波形シンセサイザ及びビームコントローラ２０８に波形データ（例えば、例示的なデジタルメッセージ／データ）を発行するプロセス２５２を含む。発行される波形データは、直交周波数及び／又は位相符号化波形として合成される望ましい周波数チップの周波数、振幅及び位相情報を含むことができ、各符号化波形は、異なる周波数帯に対応する。方法２５０は、ＳＡＵシステム２００又は２００ａの送信モードをイネーブルにするプロセス２５３を含む。例えば、ＳＡＵシステム２００をインプリメントする場合、プロセス２５３は、システムコントローラ２０２を用いて、Ｎ極双投式Ｔ／Ｒスイッチ、例えば、Ｔ／Ｒスイッチ２１０を送信位置にすることによって、ＳＡＵシステム２００を送信モードに切り換える。プロセス２５３は、方法２５０の間、複数のインスタンスで実行してもよく、例えば、後述するプロセス２５８において応答エコー信号を受信した後に実行してもよい。方法２５０は、定義済みの周波数チップに対応する個別のアナログＲＦ波形を生成するプロセス２５４を含む。例えば、波形シンセサイザ及びビームコントローラ２０８のアレイの各素子は、波形発生器２０７からのデジタルメッセージ／データを、コヒーレントなアナログ広帯域複合波形を構成することができる個別のアナログ波形に変換することができる。幾つかの実現例では、方法２５０は、コヒーレントなアナログ広帯域複合波形を構成することができる個別のアナログ波形を増幅するプロセス２５５を含むことができる。例えば、各アナログ波形及び／又は広帯域複合波形は、出力増幅器２０９のアレイ素子によって増幅してもよい。増幅されたアナログ広帯域複合波形は、Ｔ／Ｒスイッチ２１０を通過し、例えば、超音波プローブ内のトランスデューサアレイ２１１の各アレイ素子を励起する。方法２５０は、複合アナログ波形を、走査されるボリュームを介して伝播することができる音響波形にトランスデュースするプロセス２５６を含む。例えば、トランスデューサアレイ２１１の各素子は、波形シンセサイザ及びビームコントローラ２０８で生成された、周波数チップに対応する個別のアナログ波形のそれぞれから、広帯域複合音響波形を構成する音響波形を提供することができる。トランスデューサアレイ２１１は、ターゲット媒体、例えば、調査対象の生物組織ボリューム内に伝播する音響ビームを生成することができる。プロセス２５６の終了後、方法２５０は、ＳＡＵシステム２００又は２００ａの受信モードをイネーブルにするプロセス２５７を実行することができる。例えば、ＳＡＵシステム２００をインプリメントする場合、プロセス２５７は、システムコントローラ２０２を用いて、Ｎ極双投式Ｔ／Ｒスイッチ２１０を受信位置にすることによって、ＳＡＵシステム２００を受信モードに切り換える。方法２５０は、応答音響波形を受信するプロセス２５８を含み、応答音響波形は、１つ以上の応答音響波形の形式であってもよい。プロセス２５８は、応答音響波形エコーを、例えば、生成された個別の波形の周波数チップに対応する個別の受信アナログ波形にトランスデュースすることを含む。例えば、応答音響波形は、トランスデューサアレイ２１１に戻り、受信される。トランスデューサアレイ２１１の各素子は、受信した音響波形をアナログ信号（波形）に変換することができる。
【００４１】
方法２５０の幾つかの具体例では、トランスデューサアレイ２１１は、ターゲットに対して他の位置に平行移動させることができる。プロセス２５３〜２５８は、ターゲットに関するトランスデューサアレイ２１１の複数の位置のそれぞれについて繰り返してもよく、これにより、図２Ｃに例示する合成アパーチャを形成することができる。方法２５０の幾つかの具体例では、トランスデューサアレイ２１１を用いて、例えば、トランスデューサアレイ２１１の対応するトランスデューサ素子の１つ以上の位置であることができる複数の位置のうちの１つ以上の位置から送信される音響波形を操縦し、図２Ｄに示すような合成アパーチャを形成することができる。例えば、送信された音響波形は、ターゲットに向かう直接的な経路内の複数の位置のうちの１つ以上の位置において、生成された合成アパーチャ画像に基づいて操縦してもよい。また、例えば、システムコントローラ２０２は、複数の符号化波形の各波形について、複数の振幅及び複数の位相を生成することができ、これらは、それぞれ個別に振幅重み付け及び個別に位相重み付けされ、これによって、ターゲットに向かう位置から音響波形の操縦を行うことができる。
【００４２】
幾つかの実現例では、方法２５０は、個別の受信アナログ波形を増幅するプロセス２５９を含むことができる。例えば、各受信アナログ波形は、前置増幅器モジュール２１２内の低雑音前置増幅器要素によって増幅することができる。方法２５０は、個別のアナログ波形をデジタル波形データに変換するプロセス２６０を含む。例えば、各受信され（及び増幅された）アナログ波形信号は、Ａ／Ｄ変換器モジュール２１３の各Ａ／Ｄ要素によってデジタルワードに変換できる。デジタルフォーマットデータは、信号処理のためにデジタルシグナルプロセッサ２０４に送信してもよい。
【００４３】
方法２５０は、デジタル波形データを、合成アパーチャ画像及びターゲット媒体を表す画像フレームに処理するプロセス２６１を含む。プロセス２６１については、図６を用いて後に更に詳細に説明する。プロセス２６１は、また、デジタル波形データを、個別の複合受信アナログ波形を表す複合デジタル信号に複合することを含むことができる。例えば、デジタルシグナルプロセッサ２０４は、合成アパーチャについて、トランスデューサアレイ素子のそれぞれによって受信された広帯域複合音響波形を構成する周波数チップの各振幅及び位相を検出することができる。デジタルシグナルプロセッサ２０４は、受信ビームを形成し、ビームの各解像度要素の振幅及ドップラ成分を分離し、先にオペレータによって選択されているモードに関連する画像フレームを形成することができる。デジタルシグナルプロセッサ２０４によって形成される画像フレームは、ディスプレイ２０５によってユーザに表示することができる。
【００４４】
ＳＡＵシステム２００は、超音波イメージングのためのスペクトル拡散、広瞬時帯域幅（例えば、比帯域が最大１００％以上）、コヒーレント、擬似ランダム雑音（ＰＲＮ）、周波数及び／又は位相符号化波形を生成するように実現することができる。このような波形の実施形態は、無限にある。図３は、複数の個別の波形（例えば、周波数チップ）からなる生成された複合波形３００の例示的なプロットの一例を示している。幾つかの実現例では、複合波形３００の個別の波形は、周波数チップ毎にパルスのシーケンスを含み、シーケンス又はコード期間（Ｔ）の後に繰り返されるＰＲＮ波形であってもよく、例えば、このシーケンスは、周波数チップのセット内の他のどのシーケンスに対しても、又は有意に異なる時間フレームにおける同じシーケンスに対しても、非常に低い相関性を有し、又は狭帯域干渉又は熱雑音を有する。例えば、ＳＡＵシステム２００は、送信側及び受信側で、例示的なＰＲＮ波形の正確に同じシーケンスを生成することができ、したがって、（送信信号シーケンスに基づく）受信信号シーケンスは、ターゲットの音響画像を生成する信号処理において、高い相関性を示すことができる。
【００４５】
図３に示すように、複合波形３００の例示的な個別の波形又はチップ３０１は、例えば、図２Ｅ内のプロセス２５６について説明したように、時間フレームｔ０から開始される送信期間Ｔの間に送信された周波数チップｆＮ−２に対応する。図３に示すように、送信期間Ｔには、受信時間間隔ＴＲが続き、ここで、図２Ｅ内のプロセス２５８について説明したように、応答音響波形エコーが受信される。送信期間Ｔ及び受信時間間隔ＴＲは、フレーム期間Ｔｆを構成し、これは、後続する時間フレーム（ｔ１、ｔ２、ｔ３…）において繰り返すことができる。
【００４６】
例示的な複合波形３００は、時間領域を複素数として表す以下の式（１）として表される波形Ｗのための式によって表すことができる。
【００４７】
【数１】
2018183651.tif 000003
【００４８】
Ｗは、Ｍ個の個別の直交波形（例えば、直交周波数チップ）から構成され、ｊ＝−√−１である。式（１）において、ｎは、複合波形Ｗ内の周波数チップの数を表し、ｋは、時間（例えば、フレーム）内の変化の数を表し、Ｔは、チップ期間又は符号化されたシーケンスの期間を表し、そしてｆ０は、基本チップ周波数を表し、ｆ０＝１／ＮＴであり、Ｎｆ０は、最大周波数であり、（Ｍ−Ｎ＋１）ｆ０は、最小周波数である。例えば、周波数チップの数ｎは、Ｎ−Ｍ＋１からＮまでの正整数のシーケンスを表す。波形繰り返し周波数は、１／Ｔｆであり、Ｔｆは、フレーム又はエポックの継続時間であり、０≦ｘ≦ｔｆの場合、Ｕ（ｘ）＝１である。Φｎｋは、ｋ番目のタイムエポックにおけるｎ番目のチップの周波数チップ位相項を表し、Ａｎはｎ番目のチップの振幅である。例えば、周波数チップ位相項Φｎｋは、擬似乱数位相項であってもよく、ここで、擬似ランダム的にスクランブルされた開始位相Φｎｋは、集合｛Ｉｎｋ２π／Ｎ｝内の乱数であり、ここで、Ｉｎｋは、Ｎを大きな数として、数列Ｉ＝０、１、２、３…Ｎから置き換えなしで選択されるランダムな正整数のシーケンスである。別の例では、周波数チップ位相項Φｎｋは、所望の程度で統計的に互いに直交する波形Ｗｓ（ｔ）の集合を生成する多くの数学的手法の１つを用いて選択することができる。付加的な位相項であるＣｎは、０から２πの間の数である。例えば、周波数チップ位相擬似乱数Φｎｋを、システムコントローラ２０２及び／又は波形発生器２０７のメモリユニット内の例示的なデータベースに予め保存してもよい。
【００４９】
複合波形Ｗは、個別の、実質的に直交する符号化波形（例えば、周波数チップ）を合成することによって形成することができ、ここで、個別の符号化波形は、異なる周波数帯に対応し、符号化波形は、周波数符号化波形又は位相符号化波形の少なくとも１つ、例えば、波形シンセサイザ２０８において合成される符号化波形を含む。複合符号化波形は、周波数チップの搬送周波数を定義する２つ以上の周波数を選択し（例えば、最小周波数及び最大周波数を選択することを含む。）、及び周波数チップのＡｎ振幅値を決定することによって、周波数符号化波形として合成することができる。周波数符号化波形の合成は、符号化波形の各波形の時間−帯域幅積（Ｍｆ０Ｔ）パラメータを決定することを含んでいてもよい。幾つかの実現例では、特定のタイムエポックの間に、特定の周波数チップのための振幅を、その周波数チップのための単一の値として決定し、その特定の周波数チップのための後続するタイムエポックにおいてその単一の値を繰り返してもよい。他の実現例では、特定のタイムエポックの間に、特定の周波数チップのための振幅を、その周波数チップのための単一の値として決定し、その特定の周波数チップのための後続するタイムエポックにおいては、異なる単一の値を割り当ててもよい。更に他の実現例では、特定のタイムエポックの間に、特定の周波数チップのための振幅が、その周波数チップのための複数の振幅値を含むように決定してもよく、ここで、特定の周波数のための後続するタイムエポックにおいて、Ａｎの複数の値を繰り返してもよく、変化させてもよい。最大周波数（Ｎｆ０）から最小周波数（（Ｍ−Ｎ１）ｆ０）までの周波数範囲並びに個別の波形振幅項の集合の選択では、多くの既知のコードシーケンス（例えば、プッシュシーケンス（pushing sequences）、バーカコード（Barker Codes）等）の１つを利用してもよく、例えば、擬似乱数コード上の数値検索を利用してもよく、又は最小曖昧サイドローブ（minimum ambiguity sidelobes）のための他の如何なるコードを利用してもよい。
【００５０】
これに加えて又はこれに代えて、複合符号化波形は、個別の符号化された直交波形（例えば、周波数チップ）の各波形の個別の波形位相項（Φｎｋ）を決定することによって位相符号化波形として合成してもよい。例えば、複合波形Ｗの変化を提供するために、位相Φｎｋは、送信期間Ｔ内の周波数チップのための１つ以上の位相値を含むことができる。幾つかの実現例では、特定のタイムエポックの間に、特定の周波数チップのための位相Φｎｋを、その周波数チップのための単一の値として決定し、その特定の周波数チップのための後続するタイムエポックにおいてその単一の値を繰り返してもよい。他の実現例では、特定のタイムエポックの間に、特定の周波数チップのための位相Φｎｋを、この周波数チップのための単一の値として決定し、その特定の周波数チップのための後続するタイムエポックにおいては、異なる単一の値を割り当ててもよい。更に他の実現例では、特定のタイムエポックの間に、特定の周波数チップのための位相Φｎｋが、その周波数チップのために複数の値を含むように決定してもよく、ここで、特定の周波数のための後続するタイムエポックにおいて、Φｎｋの複数の値を繰り返してもよく、変化させてもよい。例えば、第１のタイムエポック（ｔ０）の波形３０１は、例えば、送信期間Ｔの前半部分のためのその位相シフトとして、第１の位相ΦＡを含むことができ、例えば、送信期間Ｔの後半部分のためのその位相シフトとして、第２の位相ΦＢを含むことができる。次のタイムエポック（ｔ１）の波形３０１は、その前半及び後半の位相シフトとして、例示的な位相ΦＡ及びΦＢを繰り返してもよく、或いは、他の位相シフトシーケンス（例えば、ΦＡ、ΦＢ、ΦＣ、又は例えば、ΦＢ及びΦＡ、又は他の構成）を含んでもよい。周波数符号化波形の合成は、符号化波形の各波形の時間−帯域幅積（Ｍｆ０Ｔ）パラメータを決定することを含んでいてもよい。
【００５１】
図３に示すように、例示的な送信複合波形Ｗは、周波数範囲ｆＮ−Ｍ＋１〜ｆＮに亘って完全に広がり、互いに直交するＭ個の個別の波形から構成さる。パラメータＮは、Ｗに広瞬時帯域幅を与えるために、大きな数を選択することができる。特定のケースでは、最低周波数がｆＮ−Ｍ＋１＝１／Ｔの場合、Ｗは、この周波数範囲に含むことができる如何なる広帯域の波形を表していてもよい。Ｍ個の個別の波形のうちのどの波形も、間隔Ｔの間に、１つ以上の位相（例えば、Φｎｋ）を１つの波形に符号化してもよい。更に、Ｍ個の個別の波形のうちのどの波形も、１つの波形に符号化された複数の振幅を含むことができる。これは、振幅重み付けと位相重み付けによって実現してもよい。
【００５２】
式（１）によって表される個別の相互に直交する波形チップのファミリは、コヒーレントな、擬似ランダム雑音、周波数符号化及び／又は位相符号化、スペクトル拡散複合波形を形成することができる。個別の波形は、パラメータの選択によって、所望の程度で統計的に互いに直交させることができる。例えば、式（２）に表す所与の波形の曖昧度関数の遅延及び周波数サイドローブレベルは、この波形の直交の程度を表している。波形の特定のパラメータを決定することによって、医療用超音波画像の解像度を大幅に向上させることができる。例えば、医療用超音波画像の解像度に影響するパラメータは、結合された軸の固有の範囲（例えば、ドップラ分解能）及び波形のスペックル低減能力を決定する時間−帯域幅積（Ｍｆ０Ｔ）パラメータ、並びに統計的な直交の程度を決定し、これにより、例えば、波形が生物組織の不均質な媒体内で機能する程度を決定する個別の波形位相項（Φｎｋ）を含む。例えば、サイドローブが小さいほど、直交性が高く、解像度が高く（雑音が少なく）なる。個別の波形位相項（Φｎｋ）の集合の選択では、（例えば、バーカ（Barker）、フランク（Frank）、ゴレイ（Golay）等を含む）多くの既知のコードシーケンスの１つを利用してもよく、例えば、擬似乱数コード上の数値検索を利用してもよく、又は最小曖昧サイドローブ（minimum ambiguity sidelobes）のための他の如何なる符号を利用してもよい。
【００５３】
幾つかの実現例では、式（１）によって表される複合波形３００は、例えば、単一の広帯域、コヒーレント、周波数及び／又は位相符号化波形であってもよい。例えば、パラメータの選択によって、単一の波形が他の全ての信号波形又はターゲット媒体に存在する雑音信号に対して統計的に直交するようにしてもよい。
【００５４】
ｎ番目のチップの振幅であるパラメータＡｎ及び付加的な位相項Ｃｎは、組合せて、トランスデューサアレイ２１１の個別の素子を励起して、Ｗの周波数範囲に亘って望ましい振幅及び位相特徴を有する送信音響波形を生成するアナログ信号のプリエンファシス（pre-emphasis）を提供することができる。送信波形のプリエンファシスは、周波数の関数としてのトランスデューサ素子の一定ではない振幅及び位相応答及び組織層の非一様な伝播特性の両方を補償することができる。例えば、プリエンファシス項は、一定の（例えば、フラットな）振幅を有するチップ及び既知の位相対周波数特性を有する音響波形を提供することができる。このような振幅対周波数が一定の音響波形は、「白色（white）」波形と呼ぶことができる。一方、プリエンファシスが提供されない場合、送信された音響波形は、トランスデューサの周波数応答を複製することができ、このような波形は、「色付き（colored）」波形と呼ばれる。受信波形のデエンファシス（de-emphasis）によって、ターゲット媒体のボリューム、例えば、生物組織ボリュームの反射特性を判定することができる。
【００５５】
式（１）によって表される複合波形Ｗは、２つ以上の個別の相互に直交する符号化波形であり、これは、チップとも呼ばれる。複合波形の個別の相互に直交する符号化波形のそれぞれは、固有の周波数を有し、各固有の周波数には、対応する特定の位相が関連付けられている。幾つかの実現例では、個別の相互に直交する符号化波形は、振幅符号化及び位相符号化でき、固有の周波数波形のそれぞれは、固有の周波数のそれぞれに関連する対応する特定の位相及び振幅を含む。幾つかの実現例では、例えば、複合波形Ｗの個別の相互に直交する符号化波形は、ターゲットに向けて順次的に又は同時に送信することができ、他の幾つかの実現例では、ターゲットに向けてランダムに送信することができる。
【００５６】
１つの実例では、複合波形Ｗ１は、互いに直交する５つの個別の符号化波形、すなわち、対応する特定の位相φ１を有する周波数ｆ１を有する第１の波形、対応する特定の位相φ２を有する周波数ｆ２を有する第２の波形、対応する特定の位相φ３を有する周波数ｆ３を有する第３の波形、対応する特定の位相φ４を有する周波数ｆ４を有する第４の波形、及び対応する特定の位相φ５を有する周波数ｆ５を有する第５の波形を含み、これは、以下のように表すことができる。
Ｗ１＝ｆ１，φ１＋ｆ２，φ２＋ｆ３，φ３＋ｆ４，φ４＋ｆ５，φ５
同様に、例えば、複合波形Ｗ２は、５つの個別の直交符号化波形、すなわち、対応する特定の位相φ６を有する周波数ｆ１を有する第１の波形、対応する特定の位相φ７を有する周波数ｆ２を有する第２の波形、対応する特定の位相φ８を有する周波数ｆ３を有する第３の波形、対応する特定の位相φ９を有する周波数ｆ４を有する第４の波形、及び対応する特定の位相φ１０を有する周波数ｆ５を有する第５の波形を含むことができ、これは、以下のように表すことができる。
Ｗ２＝ｆ１，φ６＋ｆ２，φ７＋ｆ３，φ８＋ｆ４，φ９＋ｆ５，φ１０
同様に、例えば、５つの個別の直交符号化波形を有する複合波形Ｗ３は、対応する特定の位相φ６を有する周波数ｆ６を有する第１の波形、対応する特定の位相φ７を有する周波数ｆ７を有する第２の波形、対応する特定の位相φ８を有する周波数ｆ８を有する第３の波形、対応する特定の位相φ９を有する周波数ｆ９を有する第４の波形、及び対応する特定の位相φ１０を有する周波数ｆ１０を有する第５の波を有する波形を含むことができ、これは、以下のように表すことができる。
Ｗ３＝ｆ６，φ６＋ｆ７，φ７＋ｆ８，φ８＋ｆ９，φ９＋ｆ１０，φ１０
同様に、例えば、５つの個別の直交符号化波形を有する複合波形Ｗ４は、対応する特定の位相φ１を有する周波数ｆ１を有する第１の波形（例えば、Ｗ１の第１の波形と同じ）、対応する特定の位相φ２を有する周波数ｆ２を有する第２の波形（例えば、Ｗ１の第２の波形と同じ）、対応する特定の位相φ１１を有する周波数ｆ８を有する第３の波形、対応する特定の位相φ８を有する周波数ｆ３を有する第４の波形（例えば、Ｗ２の第３の波形と同じ）及び対応する特定の位相φ１０を有する周波数ｆ１０を有する第５の波形（例えば、Ｗ３の第５の波形と同じ）を含むことができ、これは、以下のように表すことができる。
Ｗ４＝ｆ１，φ１＋ｆ２，φ２＋ｆ８，φ１１＋ｆ３，φ８＋ｆ１０，φ１０
同様に、例えば、５つの個別の直交符号化波形を有する複合波形Ｗ５は、対応する特定の位相φ１２を有する周波数ｆ１を有する第１の波形、対応する特定の位相φ１２を有する周波数ｆ２を有する第２の波形、対応する特定の位相φ１２を有する周波数ｆ８を有する第３の波形、対応する特定の位相φ１２を有する周波数ｆ３を有する第４の波形、及び対応する特定の位相φ１０を有する周波数ｆ１０を有する第５の波形（例えば、Ｗ３の第５の波形と同じ）含むことができ、これは、以下のように表すことができる。
Ｗ５＝ｆ１，φ１２＋ｆ２，φ１２＋ｆ８，φ１２＋ｆ３，φ１２＋ｆ１０，φ１０
これらの例示的な複合波形（Ｗ１、Ｗ２、Ｗ３、Ｗ４、Ｗ５）の全ては、互いに直交させることができ、或いは、必要に応じて、低い相互相関性を有するように設計することができる。
【００５７】
他の具体例では、複合波形Ｗ６は、互いに直交する５つの個別の符号化波形、すなわち、対応する特定の位相φ１及び振幅Ａ１を有する周波数ｆ１を有する第１の波形、対応する特定の位相φ２及び振幅Ａ２を有する周波数ｆ２を有する第２の波形、対応する特定の位相φ３及び振幅Ａ３を有する周波数ｆ３を有する第３の波形、対応する特定の位相φ４及び振幅Ａ４を有する周波数ｆ４を有する第４の波形、及び対応する特定の位相φ５及び振幅Ａ５を有する周波数ｆ５を有する第５の波形を含み、これは、以下のように表すことができる。
Ｗ６＝（ｆ１，φ１，Ａ１）＋（ｆ２，φ２，Ａ２）＋（ｆ３，φ３，Ａ３）＋（ｆ４，φ４，Ａ４）＋（ｆ５，φ５，Ａ５）
同様に、例えば、５つの個別の直交符号化波形を有する複合波形Ｗ７は、対応する特定の位相φ６及び振幅Ａ６を有する周波数ｆ１を有する第１の波形、対応する特定の位相φ７及び振幅Ａ７を有する周波数ｆ２を有する第２の波形、対応する特定の位相φ８及び振幅Ａ８を有する周波数ｆ３を有する第３の波形、対応する特定の位相φ９及び振幅Ａ９を有する周波数ｆ４を有する第４の波形、及び対応する特定の位相φ１０及び振幅Ａ１０を有する周波数ｆ５を有する第５の波形、を含むことができ、これは、以下のように表すことができる。
Ｗ７＝（ｆ１，φ６，Ａ６）＋（ｆ２，φ７，Ａ７）＋（ｆ３，φ８，Ａ８）＋（ｆ４，φ９，Ａ９）＋（ｆ５，φ１０，Ａ１０）
同様に、例えば、５つの個別の直交符号化波形を有する複合波形Ｗ８は、対応する特定の位相φ６及び振幅Ａ６を有する周波数ｆ６を有する第１の波形、対応する特定の位相φ７及び振幅Ａ７を有する周波数ｆ７を有する第２の波形、特定の位相φ８及び振幅Ａ８を有する周波数ｆ８を有する第３の波形、対応する特定の位相φ９及び振幅Ａ９を有する周波数ｆ９を有する第４の波形、及び対応する特定の位相φ１０及び振幅Ａ１０を有する周波数ｆ１０を有する第５の波形を含むことができ、これは、以下のように表すことができる。
Ｗ８＝（ｆ６，φ６，Ａ６）＋（ｆ７，φ７，Ａ７）＋（ｆ８，φ８，Ａ８）＋（ｆ９，φ９，Ａ９）＋（ｆ１０，φ１０，Ａ１０）．
同様に、例えば、５つの個別の直交符号化波形を有する複合波形Ｗ９は、対応する特定の位相φ１及び振幅Ａ１を有する周波数ｆ１を有する第１の波形（例えば、Ｗ６の第１の波形と同じ）、対応する特定の位相φ２及び振幅Ａ２を有する周波数ｆ２を有する第２の波形（例えば、Ｗ６の第２の波形と同じ）、対応する特定の位相φ１１及び振幅Ａ１１を有する周波数ｆ８を有する第３の波形、対応する特定の位相φ８及び振幅Ａ８を有する周波数ｆ３を有する第４の波形（例えば、Ｗ７の第３の波形と同じ）、及び対応する特定の位相φ１０及び振幅Ａ１０を有する周波数ｆ１０を有する第５の波形（例えば、Ｗ８の第５の波形と同じ）を含むことができ、これは、以下のように表すことができる。
Ｗ９＝（ｆ１，φ１，Ａ１）＋（ｆ２，φ２，Ａ２）＋（ｆ８，φ１１，Ａ１１）＋（ｆ３，φ８，Ａ８）＋（ｆ１０，φ１０，Ａ１０）
これらの例示的な複合波形（Ｗ６、Ｗ７、Ｗ８、Ｗ９）の全ては、互いに直交させることができ、或いは、必要に応じて、低い相互相関性を有するように設計することができる。
【００５８】
単一周波数モード（例えば、従来のＡ−モード、Ｂ−モード及びＣ−モード）は、これらのモノクロマチックな性質のために、プリエンファシスを必要としないことが検証された。このような単一周波数波形は、例えば、生物学的に安全な音響強度制約を確実にするために、振幅制御を必要とすることがある。
【００５９】
各チップの位相がランダムである場合、送信波形Ｗは、ランダム雑音のような特徴を有することがある。各チップの位相（Φｎｋ＋Ｃｎ）が一意的に決定され、反復可能で、（図２Ａに示すように）マスタクロックに同期される場合、送信波形Ｗは、擬似ランダム雑音として分類することができる。このような擬似ランダム雑音波形は、コヒーレントであり、これにより、コヒーレントな受信機を実現することができる。
【００６０】
広瞬時帯域幅擬似ランダム雑音波形の画像処理の利点は、適切な波形選択による、従来の医療用超音波画像に共通するスペックル、例えば、相互干渉波形によって生成される、ランダムな強度パターンであるスペックル／スペックルパターンの低減及び潜在的除去を含むことができる。この例示的なスペックルの低減は、スペックルが観察されない広帯域ガウス雑音状白色光によって照射されたシーンと、同じシーンで強いスペックルを示す狭帯域レーザ照射との比較に類似する。
【００６１】
コヒーレントな擬似ランダム雑音周波数及び位相符号化波形の信号処理の利点は、タイムサイドローブ及びドップラサイドローブが非常に小さい波形を含む。例えば、曖昧度関数Ａ（τ，ν）は、ドップラシフト（τ）又は伝播遅延（ν）の影響のために受信機内の整合フィルタによって処理される受信波形の歪みを示す２次元表現であってもよい。詳しくは、例示的な曖昧度関数Ａ（τ，ν）は、式（２）によって、波形特性及び受信機の特徴によってのみ定義され、シナリオに依存しない。曖昧度関数Ａ（τ，ν）は、以下のように定義される。
【００６２】
【数２】
2018183651.tif 000004
【００６３】
ここで、０≦ｔ≦Ｔの場合、
【００６４】
【数３】
2018183651.tif 000005
【００６５】
であり、その他の場合、Ｘｋ（ｔ）＝０である。
【００６６】
式（１）によって表されるタイプについては、以下の式が得られる。
【００６７】
【数４】
2018183651.tif 000006
【００６８】
ここで、Δｔ＝τ−ｔ、Δｆ＝ν−（ｆｎ−ｆｍ）及びΔΦ＝Φｎ−Φｍであり、これにより、式（４）に示す完全曖昧式（complete ambiguity equation）が得られる。
【００６９】
【数５】
2018183651.tif 000007
【００７０】
ここで、ｎ及びｍは、Ｎ−Ｍ＋１からＮまでの正整数のシーケンスである。
【００７１】
図４は、波形Ｗについての式（１）によって表される擬似ランダム雑音周波数符号化複合波形４０１の例示的な曖昧度関数の特徴を表している。例示的な周波数符号化複合波形４０１は、１２８のコード長を有する。図４に示すように、この曖昧度関数のサイドローブ４０２は、チップ間の位相相互作用に起因し、遅延及び周波数の両方において、Ｎ２の関数であるピークを下回るプラトーレベル（plateau level）を有する。
【００７２】
選択される特定のコード（Ｉｎｋ）に応じて、多くの複合波形Ｗが可能であることが検証された。しかしながら、定義される全ての波形毎にサイドローブ性能が保証できるわけではなく、したがって、可能なコードの集合の数値検索によって判定される遅延及び周波数の両方においてサイドローブが十分に小さくなるコードのみを使用する必要がある。
【００７３】
例えば、医療用超音波アプリケーションでは、伝播媒体としての生体組織は、不均質である。伝播媒体の不均質性によって、差動時間遅延が生じることがあり、また、生体組織は、非意図的な動きによってドップラ効果を引き起こすことがある。超音波振動子アレイも、（例えば、物理的寸法の制約のために）超音波ビームの軸外部分に好ましくないサイドローブ及びグレーティングローブを有することがあり、これにより、メインローブの応答に望ましくない時間遅延及びドップラ応答が加わることがある。低曖昧度関数サイドローブを示す波形は、差動時間遅延、動きによって引き起こされるドップラ効果及びトランスデューササイドローブ効果からの干渉を低減することによって、ＳＡＵフォーカシング及びターゲットコントラストを著しく向上させることができる。
【００７４】
コヒーレントな擬似ランダム雑音、周波数及び／又は位相符号化波形は、より高次のクロス範囲フォーカシング技術の採用を可能にし、これは、サイズが限られた超音波振動子アレイ、例えば、医学的な超音波振動子アレイの横方向の解像度を向上させる。
【００７５】
例えば、各生物組織タイプ及び各病変組織タイプは、周波数、モノスタティック及びバイスタティック角度、並びに空間的形状の関数として、それら自身の固有の超音波エコー応答信号を示すことがある。従来のエラストグラフモード（Elastograph-Mode：Ｅモード）方式を用いた場合、例えば、測定誤差のために、例えば、重なっている不均質な媒体を介する超音波伝播を正確に特徴付けることができないために、このような特性を利用して組織を分類することが困難であることがある。例示的なＳＡＵシステム２００によって生成される例示的な波形、例えば、広瞬時帯域幅、コヒーレント、擬似ランダム雑音、周波数及び／又は位相符号化波形によれば、検査中のターゲットボリュームのエコー応答から物理的組織の特徴を判定することによって組織の弁別が可能である。受信エコーの測定された特徴から得られる特徴データに分類器、例えば、ベイズ推論分類器等を適用することによって、ターゲットボリューム内に観察される組織タイプを自動的に分類することができ、コンピュータ支援診断モード（Computer Aided Diagnostic-Mode：ＣＡＤモード）が実現される。
【００７６】
生来的に画質を実質的に低下させ、個々のオペレータの技能に依存する従来のＥモードとは異なり、式（１）によって表される例示的な波形は、ＡＴＳ及び／又はＣＡＤモードにおいて、生来的に画質を向上さることができるとともに、結果として生成される画像において、組織毎に色付けを行うことができる。この利点によって、ユーザの技能に関する要求が緩和され、病変の縁を識別でき、診断を向上させることができる。
【００７７】
更に、（図２Ａに示すように）送信側に配置された波形シンセサイザ２０８及び受信側に配置されたデジタルシグナルプロセッサ２０４によって、ビーム制御機能（例えば、ビーム操縦、動的ビームフォーカシング及びビーム整形）を実現することができる。図５Ａ〜図５Ｃは、フェーズドアレイの各素子間の差動時間遅延又はこれと等価な位相シフト、及び振幅重み付けを導入することによる、これらのデジタル電子機能の基本を示している。図５Ａに示すように、差動位相シフトは、各音線（ｒ１、ｒ２…ｒｉ…）がｉ番目の素子から焦点（ｐ）に至る距離（ｄ）における差動変化を補償することができる。角度（φ）は、焦点（ｐ）がトランスデューサアレイ２１１のターゲット媒体内のターゲットに向けての直接的な狙い／整列がｚ軸方向に沿っていないために形成される。更に、各要素に差動振幅重み付け適用することによって、ビーム形状を制御し、サイドローブ及びグレーティングローブを抑制することができる。また、波形発生器２０７は、例示的な波形内の１つ以上のチップについて、位相遅延をプリエンコードして、トランスデューサアレイ２１１の各トランスデューサ素子から送信される１つ以上のチップの位相を遅延させることができる。図５Ｂ及び図５Ｃは、この特徴の例示的な結果を示している。１つ以上のチップのための例示的な位相遅延値は、デジタルシグナルプロセッサ２０４及び／又はシステムコントローラ２０２に通信することができ、これによって、受信複合波形の信号処理に位相遅延値を組み込むことができる。
【００７８】
狭瞬時帯域幅超音波デバイスの場合、この機能は、各素子を駆動する複合アナログ信号に位相シフト及び振幅減衰を導入することによって達成される。なお、ＳＡＵシステム２００によって生成される例示的なスペクトル拡散、広瞬時帯域幅、周波数及び位相符号化波形では、式（５）に示すように、波形の個々のチップ（Ｗｉ）は、Ｘ個の素子の全てについて、アレイ素子（ｉ）毎に、周波数（ｎ）の関数として、個別に振幅重み付け（Ｂｎｉ）及び位相重み付け（Ｄｎｉ）される。
【００７９】
【数６】
2018183651.tif 000008
【００８０】
必要な各チップの振幅及び位相重みは、送信時に、システムコントローラ２０２によって算出し、波形発生器２０７に指示として供給してもよい。そして、波形発生器２０７は、デジタルワード（実数成分及び虚数成分）を波形シンセサイザ及びビームコントローラ２０８に供給し、波形シンセサイザ及びビームコントローラ２０８は、アナログ駆動信号を生成し、このアナログ駆動信号は、増幅器２０９によって増幅され、トランスデューサアレイ２１１のアレイの各要素に供給される。
【００８１】
合成アパーチャ超音波信号処理は、受信した合成アパーチャ超音波信号データを収集し、送信波形の複製を利用する一連のアルゴリズムを用いてこれらのデータを処理し、合成アパーチャ超音波画像を生成することを含む。例えば、ここに開示する技術の合成アパーチャ超音波信号処理は、図４に示すように、例えば、「画鋲（thumb-tack）型」曖昧度関数を利用することによって、生体の不均質な組織の伝播効果を打ち消す。受信時には、Ａ／Ｄ変換器モジュール２１３が受信アナログ信号をデジタルワードに変換することができる。Ａ／Ｄ変換器モジュール２１３は、アレイ素子毎に実際のアレイ素子の各サンプル位置について、各チップの振幅及び位相情報を表すデジタルワードをデジタルシグナルプロセッサ２０４に供給することができ、このデータは、合成アパーチャ全体に亘って収集及び保存される。
【００８２】
合成アパーチャ波形の軸方向分解能（∂ａ）は、表現∂ａ＝ｃ／Ｂによって与えられ、Ｂは、波形の帯域幅であり、ｃは、音速である。実際のアパーチャの横方向分解能（∂ｌ）は、表現∂ｌ＝ｄλｃ／ｗＲによって近似される。例えば、ＳＡＵシステム２００によって形成される合成アパーチャのための横方向分解能は、∂ｌ＝ｄλｃ／２Ｗによって近似され、ここで、λｃは、送信波形の中心周波数の波長である。
【００８３】
図６は、デジタル信号処理を用いて超音波画像を生成する例示的な合成アパーチャ超音波信号処理技術６００のブロック図である。この技術は、例示的なコヒーレント、広帯域、瞬間的、スペクトル拡散、符号化、雑音状超音波波形を処理することによって実現してもよい。例示的な合成アパーチャ超音波信号処理技術は、ＳＡＵシステム２００のデジタルシグナルプロセッサ２０４を用いて、それぞれ図２Ｃ及び図２Ｄに示す例示的なＳＡＳＳ及びＳＡＳｐｌモードで実現してもよい。図６に示すように、ＳＡＵ信号処理技術６００は、各サンプリング位置について、送信超音波パルス及び応答エコー信号の各周波数チップの振幅及び位相情報を表すデジタルワードをＳＡＵシステム２００のメモリユニットに保存するプロセス６１０を含む。幾つかの実現例では、例えば、ＳＡＵ信号処理技術６００は、合成アパーチャの形成の前に軸範囲圧縮（axial range compression）を実行するプロセス６２０を含むことができ、これは、例えば、メモリを節約できる利点を有する。ＳＡＵ信号処理技術６００は、合成アパーチャ全体を含む軸範囲圧縮されたエコーデータの保存されたブロックを処理することによって、合成アパーチャ画像を形成するプロセス６３０を含むことができる。ＳＡＵ信号処理技術６００は、ＳＡ画像データを処理して、表示のための超音波画像を形成するプロセス６４０を含むことができる。
【００８４】
図７は、例示的な合成アパーチャによりサンプリングされた受信エコーデータ及びメモリへの保存を示している。この例示的なＳＡＵデータ収集シナリオでは、単にプローブと呼ぶこともできるトランスデューササブアレイの位相中心は、より大きなアレイの一部として、任意であるが既知の開いた又は閉じた３Ｄ経路に沿って示されるＳＡ座標系において機械的に又は電子的に平行移動され、通常のケースでは、各ｕ，ｖ，ｗデータ点において及び各θ、φビーム角について、生のエコー応答データをサンプリングする。サンプル点毎のデジタル化された受信エコーデータ及びその３Ｄ位置は、完全な走査毎にメモリブロックに保存される。図６のプロセス６１０に示すように、受信され、デジタル化されたエコー応答は、収集され、メモリに保存される。図７は、保存されたデータブロックの例示的な合成アパーチャ受信データメモリマップの概念を示している。収集され、デジタル化された受信エコーデータは、データのブロックとして、ブロック毎に処理することができる。
【００８５】
例えば、ＳＡＵシステム２００のデジタルシグナルプロセッサ２０４を用いて軸範囲圧縮を実行するプロセス６２０は、以下の例のように実現してもよい。図８は、例示的な広帯域、スペクトル拡散、雑音状、コヒーレント合成アパーチャ画像形成のための波数アルゴリズム（wave-number algorithm）のブロック図である。この例では、データ収集サンプル点毎の受信信号Ｙｉ（τ，ν）を、時間ｔ０によってシフトされた送信デジタル波形の複素共役複製Ｗｉ＊（τ）に乗算してもよい。この乗算演算は、例えば、通常、平行してＪ回繰り返すことができ、ここで、図８に示すように、演算毎に波形複製を先の演算から整数増分（integral increment）ｊ・ｔ０だけ時間的にシフトさせ、図７に示すように、例示的なメモリマップ内のサンプリングされたデータ点に基づいて、メモリに保存する。このステップは、サンプリングされたデータ点毎に繰り返され、図２Ｃ及び図２Ｄに示す軸範囲スライスｊについての圧縮された軸範囲データが生成される。
【００８６】
例えば、ＳＡＵシステム２００のデジタルシグナルプロセッサ２０４を用いて合成アパーチャ画像を形成するプロセス６３０は、幾つか異なるアルゴリズムのセットの何れを用いて実現してもよい。例えば、これらの合成アパーチャアルゴリズムのセットは、周波数領域アルゴリズム又は時間領域アルゴリズムとして分類することができる。例えば、周波数領域タイプのアルゴリズムは、波数（wave-number）アルゴリズム、チャープスケーリング（chirp-scaling）アルゴリズム、及びスケーリングされた逆フーリエ変換アルゴリズム等を含むことができる。また、例えば、時間領域タイプのアルゴリズムは、明示的整合フィルタ（explicit matched-filter）アルゴリズム、及び逆投影（back-projection）アルゴリズム等を含むことができる。これらの例示的なアルゴリズムの全ては、合成アパーチャ、コヒーレント、広瞬時帯域幅、スペクトル拡散、医療用超音波画像の形成に適用できるが、例えば、波数アルゴリズムが概念的に最も理解しやすいため、図８の例では、周波数領域波数アルゴリズムを示している。
【００８７】
プロセス６３０は、図８に示すように、生成された圧縮軸範囲データの窓フィルタリングステップを含むことができる。例えば、有限期間信号レコードの端部における不連続性に起因するスペクトル漏洩を最小化するために、窓関数フィルタリングを適用して、画質を改善してもよい。この例では、軸範囲スライスｊ毎に、生成される積に窓関数、例えば、矩形（rectangular）、ハン（Hann）、ハミング（Hamming）、テューキー（Tukey）、カイザーベッセル（Kaiser Bessel）、ドルフ−チェビシェフ（Dolph-Chebyshev）、ブラックマン−ハリス（Blackman-Harris）等の窓関数を乗算する。これらの例示的な窓関数は、端部においてゼロとなるように傾斜されている。窓関数の選択では、スペクトル漏洩によって誘導される画像サイドローブの抑制と、画像解像度との間でトレードオフが生じる。例えば、窓をゼロに向けて滑らかに傾斜させると、望ましくない画像サイドローブが低減されるが、解像度が犠牲となる。逆に、例えば、端部におけるゼロに向けて関数をより「急峻」にすれば、解像度はより高くなるが、画像サイドローブも大きくなる。ある例では、矩形窓関数によって、解像度が最高になり及びサイドローブも最大になる。これに代えて、例えば、ブラックマン−ハリス窓関数では、例えば、矩形窓関数に比べて、サイドローブが小さくなり（例えば、−７０ｄＢ以下）、解像度も低くなる。プロセス６３０では、適用される特定の１又は複数の窓関数の選択は、試行錯誤によって決定されるので、例えば、ＳＡＵシステム２００のオペレータに任されている。
【００８８】
プロセス６３０は、図８に示すように窓関数フィルタ処理されたデータのフーリエ変換ステップを含むことができる。幾つかの具体例では、合成アパーチャ全体について先に収集され、メモリに保存されている、デジタル化され、範囲圧縮された受信エコーデータのブロックに２次元離散フーリエ変換（Discrete Fourier TransformＤＦＴ）を適用することができる。この処理で使用できる多くの代替の離散フーリエ変換アルゴリズムがある。如何なるＤＦＴでも十分であるが、例示的なＤＦＴクラスのうち、クーリ−テューキー（Cooley-Tukey）アルゴリズムを適用してもよい。プロセス６３０において、適用されるフーリエ変換アルゴリズムの選択は、ＳＡＵシステム２００のオペレータに任せてもよい。
【００８９】
プロセス６３０は、図８に示すように、フーリエ変換ステップの後に、ストルト（Stolt）変換ステップを含むことができる。ストルト変換は、フーリエ変換されたデータの合成アパーチャの全体に亘って収集された横方向の圧縮のために使用することができる。例えば、波数アルゴリズムは、空間周波数領域の補間による座標変換（ストルトマッピング）によって、イメージングシステムデータを基本的に逆にする。フーリエ変換ステップにおいて波数領域に変換された圧縮エコーデータは、合成焦点範囲において、ターゲットに関して整合フィルタリングが適用され、これに続いて、１ステップの座標変換が行われる。ストルト変換は、図１Ｂ及び図２Ｄに示すように、合成アパーチャに沿って、各点Ｕにおいて収集された複数の極性データセットサンプルを単一の画像に変換し、図９に示すように、原点においてフォーカシングされた新たなＸ−Ｙ座標系にこれをマッピングする。
【００９０】
図９は、極座標から合成アパーチャに亘って収集された複数のサンプル画像を、所望の焦点を原点とする新しい座標系において、１つの解像度がより高い画像に変換するストルト変換プロセスの例示的な具体例を示す図である。説明を目的として、この例の合成アパーチャは、Ｕ軸に沿って平行移動する１次元トランスデューサプローブを含むと仮定する。通常の変換に続いて、伝播の方向における時間的波数（temporal wavenumber）をＫ＝２πｆ／ｃと定義する。空間的波数（spatial wavenumber）は、ＫＵ＝ＫｓｉｎΦと定義され、ここで角度Φは、以下の式によって与えられる。
【００９１】
【数７】
2018183651.tif 000009
【００９２】
これはＸ−Ｙ座標系の原点に中心合わせ（例えば、フォーカシング）されている。空間的波数ＫＸ及びＫＹは、それぞれＸ軸及びＹ軸の方向におけるものである。図９に示すこの特定の幾何学的構成例では、ストルトマッピングは、以下のように表される。
【００９３】
【数８】
2018183651.tif 000010
【００９４】
図８に示すように、プロセス６３０は、ストルト変換ステップの後に、逆フーリエ変換ステップを含むことができる。先のストルト変換ステップの結果に２次元逆離散フーリエ変換を適用してもよい。例えば、逆ＤＦＴは、１／Ｎ係数及び指数の符号が逆である点を除いて、基本的にＤＦＴと同じアルゴリズムであるため、ＤＦＴと同様に他の多くの逆ＤＦＴをこのステップに適用することができる。例えば、一般的に用いられるＤＦＴは、クーリ−テューキー（Cooley-Tukey）アルゴリズムを含むが、如何なる逆ＤＦＴアルゴリズムでも十分であり、如何なるＤＦＴアルゴリズムを適用してもよい。プロセス６３０では、適用される逆フーリエ変換アルゴリズムの選択は、ＳＡＵシステム２００のオペレータに任せてもよい。
【００９５】
例えば、プロセス６３０の幾つかの具体例では、オートフォーカスステップを実行してもよい。例えば、合成アパーチャ開ループ形成画像（synthetic aperture open-loop formed image）のフォーカシングは、結果の画像の点像分布関数（Point Spread Function：ＰＳＦ）のサイドローブを低減する周知の幾つかの合成アパーチャレーダ（Synthetic Aperture Radar：ＳＡＲ）及び合成アパーチャソナー（Synthetic Aperture Sonar：ＳＡＳ）オートフォーカスアルゴリズムを用いて任意に向上させることができる。受信合成アパーチャデータは、ブロック処理されているため、画像における焦点から遠い部分は、任意の位相誤差のためにフォーカシングが不十分であることがある。１つ以上のオートフォーカス技術を適用して、少なくとも部分的な修正を行うことによって、画質を改善してもよい。
【００９６】
合成アパーチャ画像形成においては、上述したような開ループ合成アパーチャデジタル処理アルゴリズムは、媒体の均質性、定常性及びシーンに対する正確な送受信トランスデューサ−サンプリング位置を仮定することができる。例えば、生きた哺乳類の被験者では、所望のレベルの均質性、定常性及びトランスデューサ位置は、完全には実現できないことがある。したがって、幾つかの例示的な実現例では、一組の開ループ合成アパーチャアルゴリズムを適用し、続いて、包括的又は一連の局所的なオートフォーカスアルゴリズムを適用し、合成アパーチャ画像を形成し、波形のコヒーレントな、スペクトル拡散特性を利用して、アパーチャサンプリング点間でコヒーレント性を有さないＰＳＦのサイドローブを低減することによって、画質を改善することが望ましい場合がある。例えば、開ループ合成アパーチャアルゴリズムは、トランスデューササンプル位置及び組織パラメータについて利用できる最良の推定を使用し、続いて、オートフォーカスアルゴリズムが最適画質を検索して、画像形成パラメータを微調整し、例えば、トランスデューサ位置、不均質組織特性及び非コヒーレントＰＳＦサイドローブにおける推定の不正確性を軽減する。オートフォーカス技術は、生データのコヒーレント性及び冗長性に依存するので、例えば、画質の大幅な改善が実現する場合、受信エコー信号を、ナイキスト制限を超えて時間的及び空間的にオーバーサンプリングすることが望ましいこともある。例えば、時間的サンプリング間隔は、１／（２Ｎｆ０）未満であってもよく空間的サンプリング間隔は、中心周波数又はｃ／［（Ｎ＋Ｍ）ｆ０］において、２分の１波長未満であってもよい。
【００９７】
図１０は、例示的なコヒーレント、反復的、閉ループ、合成アパーチャ画像オートフォーカスサブルーチンのブロック図である。例示的なコヒーレント、反復的、閉ループ、オートフォーカス画像補正プロセスを実行することによって、デジタル信号処理の負荷が大きくなることを代償にして、軸方向−横方向の合成アパーチャ位相履歴における好ましくない変化についての画質補償を向上させることができる。
【００９８】
図１０の例示的な反復的、閉ループ、合成アパーチャ画像オートフォーカスサブルーチンでは、まず、生の受信エコーデータに基づく位相誤差を推定し、続いて、推定された／非コヒーレントな位相誤差を除去する。例えば、再構築された画像における様々な位相誤差の影響を示す表１に示すように、位相誤差の性質及び大きさに応じて、位相誤差は、幾何学的性質、解像度、コントラスト及び低減された信号対雑音比（ＳＮＲ）に関して、画質を大幅に低下させることがある。
【００９９】
【表１】
2018183651.tif 000011
【０１００】
多くのオートフォーカスアルゴリズムが知られているが、これらは、一般的に、非パラメトリック、モデルベース又はこれらの２つのタイプの組合せにグループ化することができる。例えば、一般的に用いられる非パラメトリックアルゴリズムは、様々なシーンに亘って高次の位相誤差を取り除く能力を示す位相傾斜アルゴリズム（phase-gradient algorithm）を含む。更に、位相傾斜アルゴリズムを強化した幾つかのアルゴリズムも開発されている。例えば、多くの例のうちの２つとして、最尤推定器である固有ベクトルアルゴリズム（eigenvector algorithm）及び位相誤差の変動を最小化する重み付け最小二乗アルゴリズムがある。
【０１０１】
モデルベースの、オートフォーカスアルゴリズムは、あらゆる物理的測定システムに存在する、系統的で、位置に依存し、緩やかに変化する位相誤差のモデルを採用する。例えば、トランスデューサ走査アセンブリの機械的構成部品のためのデバイスに依存するサンプリング位置誤差モデルが開発されている。このようなモデルが与えられると、位相誤差が推定され、幾つかの所定の画質測定基準に基づいて、最良の画像が得られるまで、位相誤差補正が繰り返される。
【０１０２】
非パラメトリック及びモデルベースのハイブリッド型の手法の一例では、走査システムの誤差モデルを用いて、各下位画像に現れる位相誤差の位置が不変となるように、画像を下位画像に分割する。次に、非パラメトリックオートフォーカスアルゴリズム、例えば、位相傾斜アルゴリズムを各下位画像に個別に適用することができる。そして、例えば、個別の下位画像を再び組み合わせ、完全なオートフォーカス画像を形成する。
【０１０３】
式５のコヒーレントな、広帯域波形の利点の１つは、受信信号がサブバンドに分割され、位相誤差推定のために、最低周波数（例えば、最長波長）サブバンドを最初に選択できる点である。この波長が最も長いサブバンドは、例えば、周波数がより高い（波長がより短い）サブバンドと比較して、不均質な組織及びサンプリング位置の不確実性に起因する位相誤差の影響が最小であるという効果がある。例えば、選択されたオートフォーカスアルゴリズム、例えば、位相傾斜アルゴリズムによって予め定義された画質のレベルを達成する際には、存在する誤差の推定を更に高精度化するために、必要に応じて、選択された徐々に波長が短くなるサブバンドデータを使用してもよい。最後に、例えば、位相補正されたサブバンドデータを再構築して、完全なオートフォーカス画像を形成する。
【０１０４】
ここに開示する技術の幾つかの適用例及び用途は、上述したシステム、方法及びデバイスの特徴を利用して実現することができる。ここに開示する技術の臨床的用途について、幾つかの具体例を説明する。
【０１０５】
１つの例示的な適用例では、結果として生成される画質、並びに例示的なスペクトル拡散超音波デバイスのＡＴＳ及びＣＡＤモードによって、プライマリケア担当医は、定期検診スクリーニングプロトコルにこの様式を組み込み、初期悪性腫瘍（例えば、ステージ０又は１）及び進行期癌を発見することができる。この適用例の結果として、このデバイスは、潜在的に、例えば、胃癌、膵臓癌、膀胱癌等の悪性腫瘍に罹患する、診断が困難な無症候性患者の生存率を高めることができる。
【０１０６】
他の例示的な適用例では、結果として生成される画質、並びに例示的なスペクトル拡散超音波デバイスのＡＴＳ及びＣＡＤモードによって、認定放射線科医は、外科生検又は切除侵襲を行う前に、腫瘍が良性か悪性かを診断できる。この適用例の結果として、放射線科医が初期悪性腫瘍（例えば、ステージ０又は１）を特定及び診断できるため、患者の生存率を潜在的に向上させることができる。更に、不必要な生検及びこれに伴う治療が困難な又は致命的な合併症、例えば、メチシリン耐性黄色ブドウ球菌（ＭＲＳＡブドウ状球菌）感染症に至る危険性を潜在的に回避できる。
【０１０７】
他の例示的な適用例では、例示的なスペクトル拡散超音波デバイスによって生成される２Ｄ又は３Ｄ画質、及びオプションとして、（例えば、３Ｄ画像の連続的な蓄積から導出することができる）４Ｄのイメージング能力を、細針生検及び他の医療処置に使用することができる。例えば、例示的なスペクトル拡散超音波デバイス（例えば、デバイスのトランスデューサプローブ）を例示的な細針生検器具に統合することによって、非常に小さな初期ステージ（例えば、ステージ０又は１）の腫瘍を細針生検して、非侵襲的診断を検証することができる。この適用例の結果として、外科医は、切開生検法を回避でき、治療が困難で致命的な合併症の危険性を回避でき、これは、患者にとって明らかに有益である。
【０１０８】
他の例示的な適用例では、このデバイスのスペクトル拡散トランスデューサプローブを最小侵襲外科的高解像度ビデオ測定器に統合することによって、光学画像と超音波画像とを融合することができる。このスペクトル拡散超音波デバイスの改善された２Ｄ又は３Ｄ画質及びオプションとしての４Ｄのイメージング能力、並びにＡＴＳ及びＣＡＤモードを用いて、このような融合されたビデオ及び超音波画像により、外科医は、健康な組織を過度に切除することなく、疾患のある組織を特定し、外科的に切除することができる。
【０１０９】
他の例示的な適用例では、この装置のスペクトル拡散、２Ｄ又は３Ｄ高解像度イメージング動作モードと、このデバイスの最小侵襲ＨＩＦＵ動作モードの統合は、正確な最小侵襲外科的治療のオプションを提供する。このスペクトル拡散超音波デバイスの改善された２Ｄ又は３Ｄ画質オプションとしての４Ｄのイメージング能力、並びにＡＴＳ及びＣＡＤモードを用いて、このような超音波画像により、外科医は、疾患のある組織を特定し、健康な組織を過度に破壊することなく、急速な温度上昇によって、疾患のある組織を外科的に破壊することができる。
【０１１０】
他の例示的な適用例では、このスペクトル拡散超音波デバイスの改善された２Ｄ又は３Ｄ画質オプションとしての４Ｄのイメージング能力、並びにＡＴＳを用いて、例示的なスペクトル拡散超音波デバイスは、カテーテル及び密封放射性線源の適切な位置への挿入を正確に案内することによって、悪性腫瘍の小線源治療の時間を短縮できる。このスペクトル拡散超音波デバイスの小線源治療への適用は、小さく、特定困難な腫瘍及びそのマージンの処置に特に有用である。
【０１１１】
他の例示的な適用例では、このスペクトル拡散超音波デバイスの改善された２Ｄ又は３Ｄ画質オプションとしての４Ｄのイメージング能力、並びにＡＴＳを用いて、例示的なスペクトル拡散超音波デバイスは、カテーテル及び医薬の適切な位置への挿入を正確に案内することによって、疾患の高用量の局所的な医薬投与治療を効果的に実現することができる。このスペクトル拡散超音波デバイスの医薬投薬治療への適用は、小さく、特定困難な腫瘍の処置に特に有用である。
【０１１２】
実施例
以下の実施例は、本技術の幾つか実施形態の実例である。本技術の他の例示的な実施形態は、以下に列挙する実施例の前に又は以下に列挙する実施例の後に、説明されたものであってもよい。
【０１１３】
本技術１つの例（実施例１）では、音響イメージングデバイスにおいて音響波形を生成する方法は、１つ以上の波形シンセサイザにおいて、ターゲットに向けて送信される１つ以上の複合波形を合成することであって、複合波形は、それぞれが互いに直交し及び異なる周波数帯に対応する個別の直交符号化波形から形成され、各個別の直交符号化波形は、対応する位相による固有の周波数を有することと、ターゲットに対する１つ以上の送信位置から、複数の音響波形からなる１つ以上の複合音響波形を送信することであって、送信は、アレイの１つ以上のトランスデューサ素子を選択し、１つ以上の複合音響波形のそれぞれの複数の個別の直交符号化波形を、１つ以上の複合音響波形のそれぞれの複数の対応する音響波形にトランスデュースすることを含むことと、ターゲットに対する１つ以上の受信位置において、送信される音響波形に対応するターゲットの少なくとも一部から戻る応答音響波形を受信することであって、受信は、アレイのトランスデューサ素子の少なくとも幾つかを選択して、応答音響波形を受信することを含むこととを含み、送信位置及び受信位置は、それぞれ、ターゲットに対するトランスデューサ素子のアレイの空間位置と、アレイのビーム位相中心位置との一方又は両方を含み、送信される音響波形及び応答音響波形は、音響イメージングデバイスの拡大された実効アパーチャを生成する。
【０１１４】
実施例２は、実施例１の方法を含み、音響波形をターゲットに送信する際、アレイのトランスデューサ素子を制御して、ターゲットに対する複合波形の向きを変更し、ターゲットが、イメージング期間に亘って、異なる波形の向きで音響波形を受信するようにする。
【０１１５】
実施例３は、実施例２の方法を含み、トランスデューサ素子の制御は、ターゲットに対して、複数の空間位置に沿ってアレイを平行移動することによって、ターゲットに対する複合波形の向きを変更することを含む。
【０１１６】
実施例４は、実施例２の方法を含み、トランスデューサ素子の制御は、アレイの１つ以上のトランスデューサ素子における送信される音響波形のビーム位相中心位置を変更することによって、ターゲットに対する複合波形の向きを変更することを含む。
【０１１７】
実施例５は、実施例１の方法を含み、複数の個別の直交符号化波形の各波形は、それぞれが個別に振幅重み付け及び位相重み付けされた複数の振幅及び複数の位相を含む。
【０１１８】
実施例６は、実施例１の方法を含み、複合波形の個別の直交符号化波形を合成することは、周波数帯を選択することと、各直交符号化波形の１つ以上の振幅、時間−帯域幅積パラメータ及び位相パラメータを判定することとを含む。
【０１１９】
実施例７は、実施例６の方法を含み、位相パラメータは、一組の擬似乱数又は一組の決定論的な数から決定される。
【０１２０】
実施例８は、実施例１の方法を含み、複数の音響波形の送信は、複数の位置の少なくとも１つの位置から音響波形を順次的に又は同時に送信することを含む。
【０１２１】
実施例９は、実施例１の方法を含み、個別の直交符号化波形は、コヒーレント波形を含む。
【０１２２】
実施例１０は、実施例１の方法を含み、複合波形に基づいて無線周波数（ＲＦ）波形を形成することを更に含み、送信される音響波形は、アレイの１つ以上のトランスデューサ素子においてＲＦベースの複合波形をトランスデュースすることによって生成される。
【０１２３】
実施例１１は、実施例１０の方法を含み、ＲＦベースの複合波形を増幅することを更に含む。
【０１２４】
実施例１２は、実施例１の方法を含み、受信した応答音響波形を増幅することを更に含む。
【０１２５】
実施例１３は、実施例１の方法を含み、受信した応答音響波形を、１つ以上の複合波形に対応する１つ以上の受信複合波形として、アナログフォーマットからデジタルフォーマットに変換することを更に含み、各受信複合波形は、ターゲットの情報を含み、情報は、受信複合波形の対応する周波数帯に関連する振幅及び位相を含む。
【０１２６】
実施例１４は、実施例１３の方法を含み、受信複合波形の対応する周波数帯の少なくとも１つの周波数帯について、振幅及び位相の少なくとも１つは、それぞれ個別に、振幅重み付け又は位相重み付けされる。
【０１２７】
実施例１５は、実施例１の方法を含み、受信した応答音響波形を処理して、ターゲットの少なくとも一部の画像を生成することを更に含む。
【０１２８】
実施例１６は、実施例１５の方法を含み、受信した応答音響波形を、１つ以上の複合波形に対応する１つ以上の受信複合波形として、アナログフォーマットからデジタルフォーマットに変換することを更に含み、各受信複合波形は、ターゲットの情報を含み、１つ以上の複合波形及び対応する１つ以上の受信複合波形をデータブロックのメモリマップ内に保存することを更に含み、各データブロックは、各サンプル点についての複合波形の受信した応答音響波形、対応する個別の直交符号化波形、及びサンプル点についての１つ以上のトランスデューサ素子の対応する位置データを保存する。
【０１２９】
実施例１７は、実施例１５の方法を含み、処理は、保存された受信複合波形の軸範囲圧縮を実行することと、周波数領域処理技術又は時間領域処理技術の１つ以上を用いて、実効アパーチャを含む軸範囲圧縮データの各保存されたブロックを処理することによって合成アパーチャ画像を形成することとを含む。
【０１３０】
本技術１つの例（実施例１８）では、合成アパーチャ音響波形イメージングシステムは、波形発生器に接続された１つ以上の波形シンセサイザを有する波形生成ユニットであって、波形発生器が提供する波形情報に応じて１つ以上の波形シンセサイザによって生成される異なる周波数帯に対応する複数の個別の直交符号化波形からなる複合波形を合成し、個別の直交符号化波形は、それぞれ互いに直交し及び異なる周波数帯に対応し、各個別の直交符号化波形は、対応する位相による固有の周波数を有する波形生成ユニットと、送信モードと受信モードの間で切り換えられる送信／受信スイッチングユニットと、送信／受信スイッチングユニットと通信して、ターゲットに対する１つ以上の送信位置から複数の音響波形からなる複合音響波形を送信し、複合音響波形の送信される音響波形は、複合波形の合成された個別の直交符号化波形に基づいており、及びターゲットに対する１つ以上の受信位置において、ターゲットの少なくとも一部から戻る複数の送信された音響波形に対応する応答音響波形を受信し、送信される音響波形及び応答音響波形は、合成アパーチャ音響波形イメージングシステムの拡大された実効アパーチャを生成し、送信位置及び受信位置は、それぞれ、空間位置と、ビーム位相中心位置との一方又は両方を含むトランスデューサ素子のアレイと、トランスデューサ素子のアレイと通信して、複数の個別の直交符号化波形を複数の対応する音響波形にトランスデュースするアレイの１つ以上のトランスデューサ素子を選択し、及び応答音響波形を受信するアレイの１つ以上のトランスデューサ素子を選択する多重化ユニットと、ターゲットの情報を提供する、トランスデューサ素子のアレイによって受信された応答音響波形を、アナログフォーマットからデジタルフォーマットに変換するアナログ／デジタル（Ａ／Ｄ）変換器のアレイと、波形生成ユニット及びＡ／Ｄ変換器のアレイと通信し、データを保存するメモリユニットと、メモリユニットに接続されて、情報をデータとして処理する処理ユニットとを含むコントローラユニットと、コントローラユニットと通信するユーザインタフェースユニットとを含む。
【０１３１】
実施例１９は、実施例１８のシステムを含み、保存されるデータは、受信した応答音響波形のデジタルフォーマット、対応する合成された複合波形、及び１つ以上の送信位置及び受信位置におけるトランスデューサ素子のアレイの対応する位置データを含む。
【０１３２】
実施例２０は、実施例１８のシステムを含み、波形生成ユニットは、送信／受信スイッチングユニットと１つ以上の波形シンセサイザとの間に配設され、複合波形を変更する１つ以上の増幅器を更に含む。
【０１３３】
実施例２１は、実施例１８のシステムを含み、送信／受信スイッチングユニットとＡ／Ｄコンバータのアレイとの間に配設され、受信した応答音響波形を変更する１つ以上の前置増幅器のアレイを更に含む。
【０１３４】
実施例２２は、実施例１８のシステムを含み、処理ユニットは、デジタル信号プロセッサを含む。
【０１３５】
実施例２３は、実施例１８のシステムを含み、コントローラユニットは、音響波形画像処理システムの要素の少なくとも１つと時間的に同期するマスタクロックを更に含む。
【０１３６】
実施例２４は、実施例１８のシステムを含み、コントローラユニットは、情報を処理して、ターゲットの少なくとも一部の画像を生成する。
【０１３７】
実施例２５は、実施例１８のシステムを含み、ユーザインタフェースユニットは、画像を表示するディスプレイと、システムの動作のための動作モードを含むユーザ入力データを受信するユーザ入力端子とを含む。
【０１３８】
実施例２６は、実施例２５のシステムを含み、動作モードは、応答音響波形から取得される少なくとも１つの測定された特性の少なくとも１つの特徴に基づいて、画像色分けを可能にする生物組織のイメージングのための少なくとも１つの人工組織染色（Artificial Tissue Staining：ＡＴＳ）モードを含む。
【０１３９】
実施例２７は、実施例２５のシステムを含み、動作モードは、応答音響波形から取得される１つ以上の測定された特性の少なくとも１つの特徴に基づいて、１つ以上のアルゴリズム分類器を用いて生物組織タイプを分類して生物組織をイメージングする少なくとも１つのコンピュータ支援診断（Computer Aided Diagnostic-Mode：ＣＡＤ）モードを含む。
【０１４０】
実施例２８は、実施例２５のシステムを含み、ディスプレイは、分類された生物組織タイプに基づいて、生物組織の色分けされた画像を表示する。
【０１４１】
実施例２９は、実施例１８のシステムを含み、音響波形をターゲットに送信する際、コントローラユニットは、トランスデューサ素子のアレイを制御して、ターゲットに対する複合波形の向きを変更し、ターゲットが、イメージング期間に亘って、異なる波形向きで複合音響波形を受信するようにする。
【０１４２】
実施例３０は、実施例１８のシステムを含み、トランスデューサ素子のアレイは、送信位置に沿って１次元、２次元又は３次元で移動して、複数の音響波形を送信し、及び受信位置に沿って１次元、２次元又は３次元で移動して、応答音響波形を受信するように動作する。
【０１４３】
実施例３１は、実施例１８のシステムを含み、アレイのトランスデューサ素子の少なくとも１つは、送信アレイの他のトランスデューサ素子から独立して、１次元、２次元又は３次元で移動することができる。
【０１４４】
本技術の１つの例（実施例３２）では、音響波形から画像を生成する方法は、異なる周波数帯に対応する複数の符号化波形を結合し、異なる周波数帯における個別の直交波信号を含む複合波形を生成することであって、符号化波形は、対応する位相及び振幅を有する固有の周波数を含むことと、複合波形を用いて、ターゲットに対する第１の空間位置からターゲットに向かって送信される音響波を生成することであって、音響波は、複合波形の符号化波形に対応する個別の音響波信号を含むことと、音響波をターゲットに送信することであって、個別の音響波信号は、互いに向きを変化させ、ターゲットが、イメージング期間に亘って、異なる波形向きで個別の音響波信号を受信するようにすることと、音響波がターゲットに送信された後に、ターゲットの少なくとも一部から応答音響信号を受信することと、ターゲットに対する少なくとも第２の位置から、結合するステップ、生成するステップ及び送信するステップを繰り返し、結合、生成及び送信のステップは、複数の位置について繰り返されて、合成アパーチャを形成することと、複数の位置からの受信した応答音響信号を、それぞれがターゲットの情報を含む対応するデジタル複合波形に変換することと、受信した複合波形を処理して、ターゲットの少なくとも一部の画像を生成することとを含む。
【０１４５】
実施例３３は、実施例３２の方法を含み、受信複合波形をリアルタイムで処理して、合成アパーチャ画像を生成することを更に含む。
【０１４６】
実施例３４は、実施例３３の方法を含み、生成された合成アパーチャ画像に基づいて、ターゲットに向かう直接的な経路内の複数の位置の１つ以上の位置において、送信された音響波の方向を操縦することを更に含む。
【０１４７】
本技術の１つの例（実施例３５）では、音響波形から画像を生成する方法において、異なる周波数帯に対応する複数の符号化波形を結合し、異なる周波数帯における個別の直交波信号を含む複合波形を生成することであって、符号化波形は、対応する位相及び振幅を有する固有の周波数を含むことと、ターゲットに向けて送信される、複合波形の符号化波形に対応する個別の音響波信号から構成される音響波を生成することと、音響波をターゲットに送信することであって、個別の音響波信号は、ターゲットに対する１つ以上の空間位置におけるビーム位相中心位置の第１の組から送信され、及び個別の音響波信号は、互いに向きを変化させ、ターゲットがイメージング期間に亘って、異なる波形向きで個別の音響波信号を受信するようにすることと、音響波がターゲットに送信された後に、ターゲットの少なくとも一部から応答音響信号を受信することと、ターゲットに対する１つ以上の空間位置におけるターゲットに対するビーム位相中心位置の少なくとも第２の組から、結合するステップ、生成するステップ及び送信するステップを繰り返し、合成アパーチャを形成することと、受信した応答音響信号を、それぞれがターゲットの情報を含む対応するデジタル複合波形に変換することと、受信した複合波形を処理して、ターゲットの少なくとも一部の画像を生成することとを含む。
【０１４８】
実施例３６は、実施例３５の方法を含み、受信複合波形をリアルタイムで処理して、合成アパーチャ画像を生成することを更に含む。
【０１４９】
実施例３７は、実施例３６の方法を含み、生成された合成アパーチャ画像に基づいて、送信された音響波の方向を操縦することを更に含む。
【０１５０】
本技術の１つの例（実施例３８）では、合成アパーチャ音響波形イメージングシステムは、波形発生器に接続された１つ以上の波形シンセサイザを有する波形生成ユニットであって、波形発生器が提供する波形情報に応じて１つ以上の波形シンセサイザによって生成される異なる周波数帯に対応する複数の個別の直交符号化波形からなる複合波形を合成し、個別の直交符号化波形は、それぞれ互いに直交し及び異なる周波数帯に対応し、各個別の直交符号化波形は、対応する位相による固有の周波数を有する波形生成ユニットと、波形生成ユニットと通信して、ターゲットに対する１つ以上の送信位置から複数の音響波形からなる複合音響波形を送信するトランスデューサ素子の送信アレイであって、複合音響波形の送信される音響波形は、複合波形の合成された個別の直交符号化波形に基づいており、送信位置は、送信アレイの空間位置及び送信アレイのトランスデューサ素子のビーム位相中心位置の一方又は両方を含む送信アレイと、波形生成ユニットと通信して、ターゲットに対する１つ以上の受信位置において、ターゲットの少なくとも一部から戻る送信された音響波形に対応する応答音響波形を受信するトランスデューサ素子の受信アレイであって、送信される音響波形及び受信される音響波形は、合成アパーチャ音響波形画像処理システムの拡大された実効アパーチャを生成し、送信位置及び受信位置は、それぞれ、空間位置と、ビーム位相中心位置との一方又は両方を含む受信アレイと、それぞれ、送信アレイ及び受信アレイと通信して、複数の個別の直交符号化波形を複数の対応する音響波形にトランスデュースする送信アレイのトランスデューサ素子の１つを選択し、及び応答音響波形を受信する受信アレイの１つ以上のトランスデューサ素子を選択する第１の多重化ユニット及び第２の多重化ユニットと、ターゲットの情報を提供する、トランスデューサ素子の受信アレイによって受信された応答音響波形を、アナログフォーマットからデジタルフォーマットに変換するアナログ／デジタル（Ａ／Ｄ）変換器のアレイと、波形生成ユニット及びＡ／Ｄ変換器のアレイと通信し、データを保存するメモリユニットと、メモリユニットに接続されて、情報をデータとして処理する処理ユニットとを含むコントローラユニットと、コントローラユニットと通信するユーザインタフェースユニットとを含む。
【０１５１】
実施例３９は、実施例３８のシステムを含み、保存されるデータは、受信した応答音響波形のデジタルフォーマット、対応する合成された複合波形、及び１つ以上の送信位置及び受信位置のそれぞれにおける送信アレイ及び受信アレイの対応する位置データを含む。
【０１５２】
実施例４０は、実施例３８のシステムを含み、波形生成ユニットは、送信／受信スイッチングユニットと１つ以上の波形シンセサイザとの間に配設され、複合波形を変更する１つ以上の増幅器を更に含む。
【０１５３】
実施例４１は、実施例３８のシステムを含み、受信アレイと、Ａ／Ｄコンバータのアレイとの間に配設され、受信した応答音響波形を変更する１つ以上の前置増幅器のアレイを更に含む。
【０１５４】
実施例４２は、実施例３８のシステムを含み、処理ユニットは、デジタル信号プロセッサを含む。
【０１５５】
実施例４３は、実施例３８のシステムを含み、コントローラユニットは、音響波形画像処理システムの要素の少なくとも１つと時間的に同期するマスタクロックを更に含む。
【０１５６】
実施例４４は、実施例３８のシステムを含み、コントローラユニットは、情報を処理して、ターゲットの少なくとも一部の画像を生成する。
【０１５７】
実施例４５は、実施例３８のシステムを含み、ユーザインタフェースユニットは、画像を表示するディスプレイと、システムの動作のための動作モードを含むユーザ入力データを受信するユーザ入力端子とを含む。
【０１５８】
実施例４６は、実施例４５のシステムを含み、動作モードは、応答音響波形から取得される少なくとも１つの測定された特性の少なくとも１つの特徴に基づいて、画像色分けを可能にする生物組織のイメージングのための少なくとも１つの人工組織染色（Artificial Tissue Staining：ＡＴＳ）モードを含む。
【０１５９】
実施例４７は、実施例４５のシステムを含み、動作モードは、応答音響波形から取得される１つ以上の測定された特性の少なくとも１つの特徴に基づいて、１つ以上の分類器を用いて生物組織タイプを分類して生物組織をイメージングする少なくとも１つのコンピュータ支援診断（Computer Aided Diagnostic-Mode：ＣＡＤ）モードを含む。
【０１６０】
実施例４８は、実施例４５のシステムを含み、ディスプレイは、分類された生物組織タイプに基づいて、生物組織の色分けされた画像を表示する。
【０１６１】
実施例４９は、実施例３８のシステムを含み、音響波形をターゲットに送信する際、コントローラユニットは、送信アレイを制御して、ターゲットに対する複合波形の向きを変更し、ターゲットが、イメージング期間に亘って、異なる波形向きで複合音響波形を受信するようにする。
【０１６２】
実施例５０は、実施例３８のシステムを含み、トランスデューサ素子の送信アレイは、複数の位置に沿って１次元、２次元又は３次元で移動して、複数の音響波形を送信するように動作する。
【０１６３】
実施例５１は、実施例３８のシステムを含み、送信アレイの１つ以上のトランスデューサ素子は、送信アレイの他のトランスデューサ素子から独立して、１次元、２次元又は３次元で移動することができる。
【０１６４】
実施例５２は、実施例３８のシステムを含み、トランスデューサ素子の受信アレイは、複数の位置に沿って１次元、２次元又は３次元で移動して、応答音響波形を受信するように動作する。
【０１６５】
実施例５３は、実施例３８のシステムを含み、受信アレイの１つ以上のトランスデューサ素子は、送信アレイの他のトランスデューサ素子から独立して、１次元、２次元又は３次元で移動することができる。
【０１６６】
実施例５４は、実施例３８のシステムを含み、送信アレイのトランスデューサ素子の数は、受信アレイのトランスデューサ素子の数より大きい。
【０１６７】
本明細書に開示した主題の具体例及び機能的動作、例えば様々なモジュールは、デジタル電子回路で実現してもよく、本明細書に開示した構造及びこれらの構造的な均等物を含むコンピュータソフトウェア、ファームウェア又はハードウェアで実現してもよく、これらの１つ以上の組合せで実現してもよい。ここに開示した主題の具体例は、１つ以上のコンピュータプログラム製品、例えば、実体がある不揮発性の媒体内に符号化され、データ処理装置によって実行され、又はデータ処理装置の動作を制御するコンピュータプログラム命令の１つ以上のモジュールとして実現することもできる。コンピュータが読取可能な媒体は、機械可読のストレージデバイス、機械可読のストレージ基板、メモリデバイス、機械可読の伝播信号に作用する組成物又はこれらの１つ以上の組合せであってもよい。用語「データ処理装置」は、データを処理するための全ての装置、デバイス及び機械を包含し、一例としてプログラミング可能なプロセッサ、コンピュータ、複数のプロセッサ又はコンピュータがこれに含まれる。装置は、ハードウェアに加えて、当該コンピュータプログラムの実行環境を作成するコード、例えば、プロセッサファームウェアを構成するコード、プロトコルスタック、データベース管理システム、オペレーティングシステム又はこれらの１つ以上の組合せを含むことができる。
【０１６８】
コンピュータプログラム（プログラム、ソフトウェア、ソフトウェアアプリケーション、スクリプト又はコードとも呼ばれる。）は、コンパイラ言語又はインタープリタ言語を含む如何なる形式のプログラミング言語で書いてもよく、例えば、スタンドアロンプログラムとして、若しくはモジュール、コンポーネント、サブルーチン又は演算環境での使用に適する他のユニットとして、如何なる形式で展開してもよい。コンピュータプログラムは、必ずしもファイルシステム内のファイルに対応していなくてもよい。プログラムは、他のプログラム又はデータを含むファイル（例えば、マークアップ言語文書内に保存された１つ以上のスクリプト）の一部に保存してもよく、当該プログラムに専用の単一のファイルに保存してもよく、連携する複数のファイル（例えば、モジュール、サブプログラム又はコードの一部を保存する１つ以上のファイル）に保存してもよい。コンピュータプログラムは、１つのコンピュータ上で実行されるように展開してもよく、１つの場所に設けられた又は複数の場所に亘って分散され、通信ネットワークによって相互接続された複数のコンピュータ上で実行されるように展開してもよい。
【０１６９】
本明細書に開示したプロセス及びロジックフローは、入力データを処理し、出力を生成することによって機能を実現する１つ以上のコンピュータプログラムを実行する１つ以上のプログラミング可能なプロセッサによって実現してもよい。プロセス及びロジックフローは、例えば、フィールドプログラマブルゲートアレイ（field programmable gate array：ＦＰＧＡ）又は特定用途向け集積回路（application specific integrated circuit：ＡＳＩＣ）等の専用論理回路等として実現できる装置によって実行してもよい。
【０１７０】
コンピュータプログラムの実行に適するプロセッサには、一例として、汎用マイクロプロセッサ及び専用マイクロプロセッサの両方、例えば、デジタルシグナルプロセッサ（digital signal processor：ＤＳＰ）並びにあらゆる種類のデジタルコンピュータの１つ以上のプロセッサの何れかを含ませてもよい。プロセッサは、通常、読出専用メモリ若しくはランダムアクセスメモリ、又はこれらの両方から命令及びデータを受け取る。コンピュータの基本的な要素は、命令を実行するプロセッサと、命令及びデータを保存する１つ以上のメモリデバイスである。また、コンピュータは、通常、データを保存するための１つ以上の大容量記憶装置、例えば、磁気ディスク、光磁気ディスク又は光ディスクを含み、若しくは、大容量記憶装置からデータを受信し、大容量記憶装置にデータを送信し、又はこの両方の動作を行うように大容量記憶装置に動作的に接続されている。但し、コンピュータは、必ずしもこのような装置を有する必要はない。コンピュータプログラム命令及びデータの格納に適するコンピュータ可読媒体には、一例として挙げれば、半導体記憶デバイス、例えば、ＥＰＲＯＭ、ＥＥＰＲＯＭ及びフラッシュメモリデバイスを含む全ての形式の不揮発性メモリが含まれる。プロセッサ及びメモリは、専用論理回路によって補ってもよく、専用論理回路に組み込んでもよい。
【０１７１】
本明細書は、多くの詳細事項を含んでいるが、これらの詳細事項は、特許請求している又は特許請求することができる本発明の範囲を限定するものとは解釈されず、本発明の特定の実施形態の特定の特徴の記述として解釈される。本明細書において、別個の実施形態の文脈で開示した幾つかの特徴を組み合わせて、単一の実施形態として実現してもよい。逆に、単一の実施形態の文脈で開示した様々な特徴は、複数の実施形態に別個に具現化してもよく、適切な如何なる部分的組合せとして具現化してもよい。更に、以上では、幾つかの特徴を、ある組合せで機能するものと説明しているが、初期的には、そのように特許請求している場合であっても、特許請求された組合せからの１つ以上の特徴は、幾つかの場合、組合せから除外でき、特許請求された組合せは、部分的組合せ又は部分的な組合せの変形に変更してもよい。
【０１７２】
同様に、図面では、動作を特定の順序で示しているが、このような動作は、所望の結果を達成するために、図示した特定の順序又は順次的な順序で行う必要はなく、また、図示した全ての動作を行う必要もない。ある特定の状況では、多重タスキング及び並列処理が有利であることもある。更に、上述した実施形態における様々なシステムの構成要素の分離は、全ての実施形態においてこのような分離が必要であることを意図してはいない。
【０１７３】
幾つかの具体例及び実施例のみを開示したが、この特許文献に記述し例示した内容に基づいて、他の具体例、拡張例及び変形例を想到することができる
【図１Ａ】
2018183651.tif 000012
【図１Ｂ】
2018183651.tif 000013
【図１Ｃ】
2018183651.tif 000014
【図１Ｄ】
2018183651.tif 000015
【図１Ｅ】
2018183651.tif 000016
【図２Ａ】
2018183651.tif 000017
【図２Ｂ】
2018183651.tif 000018
【図２Ｃ】
2018183651.tif 000019
【図２Ｄ】
2018183651.tif 000020
【図２Ｅ】
2018183651.tif 000021
【図３】
2018183651.tif 000022
【図４】
2018183651.tif 000023
【図５Ａ】
2018183651.tif 000024
【図５Ｂ】
2018183651.tif 000025
【図５Ｃ】
2018183651.tif 000026
【図６】
2018183651.tif 000027
【図７】
2018183651.tif 000028
【図８】
2018183651.tif 000029
【図９】
2018183651.tif 000030
【図１０】
2018183651.tif 000031