(57)【要約】  本発明による超音波診断システムは、被検体の周囲に配置され、超音波の送受信を行う複数の素子と、少なくとも１つの素子が超音波を送信し、少なくとも一部の素子が散乱波を受信するように制御する制御部と、前記散乱波を受信した素子から得たデータである測定データを収集する収集部と、前記被検体の少なくとも一部が含まれる撮像領域を分割した複数の分割領域のそれぞれの、所定の素子から送信した超音波が分轄領域における被検体で散乱して複数の素子の各々で受信されるまでの時間である到達時間を成分とした第１要素と、前記測定データを成分とした第２要素とに基づき、それぞれの分割領域における散乱音圧強度を算出する計算部と、各分割領域の散乱音圧強度を画素値に変換し、散乱画像を作成する画像作成部と、を備える。
2019082892.tif 000002
【特許請求の範囲】
【請求項１】
被検体の周囲に配置され、超音波の送信及び受信の少なくともいずれか一方を行う複数の素子と、
前記複数の素子の少なくとも１つが超音波を送信し、前記複数の素子の全部または一部が、前記超音波が前記被検体で散乱した散乱波を受信するように、前記複数の素子を制御する制御部と、
前記散乱波を受信した素子から得たデータである測定データを収集するデータ収集部と、
前記被検体の全部または一部が含まれる撮像領域を分割した分割領域のそれぞれの、前記複数の素子のうち所定の素子から送信した超音波が前記分割領域における前記被検体で散乱して前記複数の素子の全部または一部の各々で受信されるまでの時間である到達時間を成分とした第１の要素と、前記測定データを成分とした第２の要素に基づき、前記分割領域のそれぞれにおける前記散乱波の音圧の強度である散乱音圧強度を算出する計算部と、
前記分割領域のそれぞれにおける前記散乱音圧強度を画素値に変換した画像である散乱画像を作成する画像作成部と、
を備える超音波診断システム。
【請求項２】
前記分割領域は、前記撮像領域を格子状に分割した領域であり、
前記第１の要素は前記到達時間を成分とした行列の逆行列であり、
前記第２の要素は前記測定データを成分としたベクトルであり、
前記計算部は、前記第１の要素と、前記第２の要素の積から、前記散乱音圧強度を算出することを特徴とする請求項１に記載の超音波診断システム。
【請求項３】
前記分割領域を構成する格子状に分割した領域の縦の個数と横の個数の積と、受信素子数と前記測定データの収集における時間軸方向のデータサンプリング点数の積が、それぞれ、前記行列の列数と行数と一致することを特徴とする請求項２に記載の超音波診断システム。
【請求項４】
前記制御部は、第１素子が超音波を送信した後、第２素子が超音波を送信するように制御し、
前記データ収集部は、前記第１素子が送信した超音波に対応する散乱波を受信した素子から第１測定データを収集し、前記第２素子が送信した超音波に対応する散乱波を受信した素子から第２測定データを収集し、
前記計算部は、前記第１素子を超音波の送信素子とした場合の第１逆行列と、前記第１測定データを並べたベクトルとの積から第１散乱音圧強度を算出し、前記第２素子を超音波の送信素子とした場合の第２逆行列と、前記第２測定データを並べたベクトルとの積から第２散乱音圧強度を算出し、前記第１散乱音圧強度と前記第２散乱音圧強度とを合成することを特徴とする請求項２又は３に記載の超音波診断システム。
【請求項５】
前記行列は、ランクが、前記格子状に分割した領域の縦の個数と横の個数の積に等しいことを特徴とする請求項３に記載の超音波診断システム。
【請求項６】
前記到達時間は、乳房内部の前記超音波の音速と乳房外部の前記超音波の音速との差分が発生することを踏まえて計算されることを特徴とする請求項５に記載の超音波診断システム。
【請求項７】
前記被検体からみて、受信素子は送信素子側に配置されていることを特徴とする請求項５に記載の超音波診断システム。
【請求項８】
前記散乱画像を所定時間毎に作成し、作成した散乱画像における画素値の変化が所定値以上となる部分を抽出することを特徴とする請求項１に記載の超音波診断システム。
【請求項９】
被検体の周囲に配置された複数の素子のいずれか１つから超音波を送信し、前記複数の素子の全部または一部で、前記超音波が前記被検体で散乱した散乱波を受信する工程と、
前記散乱波を受信した素子から得たデータである測定データを収集する工程と、
前記被検体の全部または一部が含まれる撮像領域を分割した分割領域のそれぞれの、前記複数の素子のうち所定の素子から送信した超音波が前記分割領域における前記被検体で散乱して前記複数の素子の全部または一部の各々で受信されるまでの時間である到達時間を成分とした第１の要素と、前記測定データを成分とした第２の要素に基づき、前記分割領域のそれぞれにおける前記散乱波の音圧の強度である散乱音圧強度を算出する工程と、
前記分割領域のそれぞれにおける前記散乱音圧強度を画素値に変換した画像である散乱画像を作成する工程と、
を備える超音波診断方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波を照射して被検体の断層画像を作成する超音波診断システム及び超音波診断方法に関する。
【背景技術】
【０００２】
非侵襲性である超音波による診断システムは、生体を直接切開して観察する外科手術の必要がないため、被検体内部の情報を診断する技術として医療分野で広く用いられている。【０００３】
超音波診断の一手法である超音波ＣＴ（Computed  Tomography）は、超音波を被検体に照射し、反射超音波や透過超音波を用いて被検体の断層画像を作成するものであり、近年の研究により、乳がんの検出に有用性があることが示されている。超音波ＣＴは、例えば、超音波の送受信を行う多数の素子をリング状に配置したリング型アレイトランスデューサを使用し、断層像を作成する。【０００４】
従来の断層像作成方法の１種である開口合成法では、まず、１つの素子から超音波を送信し、エコー信号を全素子で受信して、第１軸が受信素子番号、第２軸がエコー信号到達時間を示す２次元データ（フレームデータ）を生成する。超音波を送信する素子を順に変えていくことで、リング型アレイトランスデューサの素子数分のフレームデータが生成される。【０００５】
図１２ａに示すように、送信素子Ｅｍから断層像の１画素に対応する着目箇所ＰＩまでの距離ＬＴＸ、この着目箇所ＰＩから受信素子Ｅｎまでの距離ＬＲＸ、及び音速ｃから、エコー信号到達時間ｔ＝（ＬＴＸ＋ＬＲＸ）／ｃが求まる。図１２ｂに示すように、送信素子Ｅｍのフレームデータのうち、受信素子Ｅｎ、時間ｔのエコーデータが、着目箇所ＰＩに対応する。【０００６】
１つのフレームデータには、受信素子毎に、着目箇所ＰＩに対応するエコーデータが含まれる。リング型アレイトランスデューサがＮ個の素子からなる場合、受信素子はＮ個となるため、１つのフレームデータには着目箇所ＰＩに対応するエコーデータがＮ個含まれる。フレームデータはＮ個あるため、着目箇所ＰＩに対応する１画素の輝度は、Ｎ×Ｎ個のエコーデータの合成となる。このようにして各画素の輝度を算出し、画像を作成していた。【０００７】
上述したように、従来は、１素子から送信した超音波のエコー信号を全素子で受信することを、素子数分繰り返し行うため、超音波の送信回数が多く、計測に時間がかかっていた。また、大量のデータを取得するため、計算機へのデータ転送に時間がかかっていた。【０００８】
【特許文献１】国際公開第２０１７／０５１９０３号
【発明の概要】
【０００９】
本発明は、上記従来の実状に鑑みてなされたものであり、被検体の計測及びデータ転送に要する時間を短縮できる超音波診断システム及び超音波診断方法を提供することを目的とする。
【００１０】
本発明による超音波診断システムは、被検体の周囲に配置され、超音波の送信及び受信の少なくともいずれか一方を行う複数の素子と、前記複数の素子のいずれか１つが超音波を送信し、前記複数の素子の全部または一部が、前記超音波が前記被検体で散乱した散乱波を受信するように、前記複数の素子を制御する制御部と、前記散乱波を受信した素子から得たデータである測定データを収集するデータ収集部と、前記被検体の全部または一部が含まれる撮像領域を分割した分割領域のそれぞれの、前記複数の素子のうち所定の素子から送信した超音波が前記分割領域における前記被検体で散乱して前記複数の素子の全部または一部の各々で受信されるまでの時間である到達時間を成分とした第１の要素と、前記測定データを成分とした第２の要素に基づき、前記分割領域のそれぞれにおける前記散乱波の音圧の強度である散乱音圧強度を算出する計算部と、前記分割領域のそれぞれにおける前記散乱音圧強度を画素値に変換した画像である散乱画像を作成する画像作成部と、を備えるものである。
【００１１】
本発明の一態様によれば、前記分割領域は、前記撮像領域を格子状に分割した領域であり、前記第１の要素は前記到達時間を成分とした行列の逆行列であり、前記第２の要素は前記測定データを成分としたベクトルであり、前記計算部は、前記第１の要素と、前記第２の要素の積から、前記散乱音圧強度を算出する。
【００１２】
本発明の一態様によれば、前記分割領域を構成する格子状に分割した領域の縦の個数と横の個数の積と、受信素子数と前記測定データの収集における時間軸方向のデータサンプリング点数の積が、それぞれ、前記行列の列数と行数と一致する。
【００１３】
本発明の一態様によれば、前記制御部は、第１素子が超音波を送信した後、第２素子が超音波を送信するように制御し、前記データ収集部は、前記第１素子が送信した超音波に対応する散乱波を受信した素子から第１測定データを収集し、前記第２素子が送信した超音波に対応する散乱波を受信した素子から第２測定データを収集し、前記計算部は、前記第１素子を超音波の送信素子とした場合の第１逆行列と、前記第１測定データを並べたベクトルとの積から第１散乱音圧強度を算出し、前記第２素子を超音波の送信素子とした場合の第２逆行列と、前記第２測定データを並べたベクトルとの積から第２散乱音圧強度を算出し、前記第１散乱音圧強度と前記第２散乱音圧強度とを合成する。
【００１４】
本発明の一態様によれば、前記行列は、ランクが、前記格子状に分割した領域の縦の個数と横の個数の積に等しい。
【００１５】
本発明の一態様によれば、前記到達時間は、乳房内部の前記超音波の音速と乳房外部の前記超音波の音速との差分が発生することを踏まえて計算される。
【００１６】
本発明の一態様によれば、前記被検体からみて、受信素子は送信素子側に配置されている。
【００１７】
本発明の一態様によれば、前記散乱画像を所定時間毎に作成し、作成した散乱画像における画素値の変化が所定値以上となる部分を抽出する。
【００１８】
本発明による超音波診断方法は、被検体の周囲に配置された複数の素子のいずれか１つから超音波を送信し、前記複数の素子の全部または一部で、前記超音波が前記被検体で散乱した散乱波を受信する工程と、前記散乱波を受信した素子から得たデータである測定データを収集する工程と、前記被検体の全部または一部が含まれる撮像領域を分割した分割領域のそれぞれの、前記複数の素子のうち所定の素子から送信した超音波が前記分割領域における前記被検体で散乱して前記複数の素子の全部または一部の各々で受信されるまでの時間である到達時間を成分とした第１の要素と、前記測定データを成分とした第２の要素に基づき、前記分割領域のそれぞれにおける前記散乱波の音圧の強度である散乱音圧強度を算出する工程と、前記分割領域のそれぞれにおける前記散乱音圧強度を画素値に変換した画像である散乱画像を作成する工程と、を備えるものである。
【発明の効果】
【００１９】
本発明によれば、被検体の計測及びデータ転送に要する時間を短縮できる。
【図面の簡単な説明】
【００２０】
【図１】本発明の実施形態における超音波診断システムの概略構成図である。
【図２】図１のＡ−Ａ線断面図である。
【図３】演算装置の機能ブロック図である。
【図４】散乱波の受信例を示す図である。
【図５】散乱波の受信例を示す図である。
【図６】関心領域のピクセル分割例を示す図である。
【図７】図７ａ〜図７ｄは測定行列の作成方法を説明する図である。
【図８】図８ａはファントムを示す図であり、図８ｂ〜図８ｆは比較例による画像作成方法を適用したシミュレーション結果を示す図である。
【図９】図９ａはファントムを示す図であり、図９ｂ〜図９ｆは比較例による画像作成方法を適用したシミュレーション結果を示す図である。
【図１０】図１０ａ、図１０ｄはファントムを示す図であり、図１０ｂ、図１０ｅは実施形態による画像作成方法を適用したシミュレーション結果を示す図である。
【図１１】図１１ａは受信開口制限の例を示す図であり、図１１ｂは透過波と散乱波の伝搬時間の例を示す図である。
【図１２】図１２ａ、図１２ｂは従来の開口合成法でのデータ取得方法を説明する図である。
【図１３】図１３ａ、図１３ｂは評価用モデルを示す図である。
【図１４】図１４ａ、図１４ｂは実施形態による再構成画像を示す図である。
【図１５】図１５ａ、図１５ｂは開口合成法による再構成画像を示す図である。
【図１６】図１６ａは信号対雑音比の解析結果を示すグラフであり、図１６ｂは解像度の解析結果を示すグラフである。
【発明を実施するための形態】
【００２１】
以下、図面を参照して本発明についてさらに詳細に説明する。本発明の実施形態に係る超音波診断システムは、人体等の被検体に超音波を照射し、受信したエコー信号を用いて散乱画像（散乱音圧強度のマップ）を作成する。医師は、作成された散乱画像を確認することで、悪性腫瘍等の病変を診断することができる。
【００２２】
図１に示すように、本実施形態に係る超音波診断システム１０は、リングアレイＲと、スイッチ回路１１０と、送受信回路１２０と、演算装置１３０と、画像表示装置１４０とを備えている。
【００２３】
リングアレイＲは、複数の振動子が組み合わさって構成される、好ましくは直径８０〜５００ｍｍ、より好ましくは直径１００〜３００ｍｍのリング型形状の振動子である。また、リングアレイＲは、直径を可変とする構成をとることもできる。本実施形態では一例として、４つの凹面型振動子Ｐ０１〜Ｐ０４を組み合わせたリング形状の振動子を用いる。
【００２４】
例えば、凹面型振動子Ｐ０１〜Ｐ０４が、それぞれ２５６個の短冊形圧電素子Ｅ（以下、単に「素子Ｅ」とも呼ぶ。）を有する場合、リングアレイＲは１０２４個の素子Ｅから構成されることになる。凹面型振動子Ｐ０１〜Ｐ０４に設けられる素子Ｅの数は限定されず、好ましくは１〜１０００個、より好ましくは１００〜５００個である。
【００２５】
各素子Ｅは、電気的信号と超音波信号とを相互変換する機能を有する。素子Ｅは被検体Ｔに超音波を送信し、被検体Ｔで散乱（反射）される散乱波（前方散乱波、側方散乱波、後方散乱波）を受信し、電気的信号を測定データとして形成する。
【００２６】
本実施形態では、各素子Ｅが、超音波の送信及び受信の両方の機能を備えるものとして説明するが、これに限定されない。例えば、超音波の送信機能及び受信機能のうちいずれか一方のみを有する送信素子又は受信素子を使用し、複数の送信素子及び複数の受信素子をリング状に配置してもよい。また、送信及び受信の両方の機能を備える素子と、送信素子と、受信素子とが混在する構成であってもよい。
【００２７】
図２は、図１のＡ−Ａ線断面図である。例えば、リングアレイＲは、穴の開いたベッドの下に、ベッドの穴と挿入部ＳＰとが重畳するように設置される。被験者はベッドの穴から、撮像対象となる身体の部位（被検体Ｔ）を挿入部ＳＰに挿入する。
【００２８】
被検体Ｔを挿入するための挿入部ＳＰは、リングアレイＲの中央に設けられている。リングアレイＲの複数の素子Ｅは、リングに沿って挿入部ＳＰの周囲に等間隔で設けられている。リングアレイＲの内周側には、音響レンズと呼ばれる凸面レンズが表面に取り付けられている。このような表面加工をリングアレイＲの内周側に施すことで、各素子Ｅが送信する超音波を、リングアレイＲを含む平面内に収束させることができる。
【００２９】
本実施形態では、各素子Ｅを等間隔にリング状に配置しているが、リングアレイＲの形状は円形に限定されず、例えば、六角形、正方形、三角形など任意の多角形、少なくとも一部に曲線や円弧を含む形状、その他任意の形状、または、これらの形状の一部（例えば、半円や円弧）であってもよい。すなわち、リングアレイＲは、アレイＲと一般化することができる。また、アレイＲを構成する各素子Ｅの配置は、被検体Ｔの周囲を断続的に少なくとも９０度またはそれ以上囲むような配置であれば好ましいものの、これらに限定されるものではない。
【００３０】
リングアレイＲはスイッチ回路１１０を介して送受信回路１２０に接続されている。送受信回路１２０（制御部）は、リングアレイＲの素子Ｅに制御信号（電気的信号）を送信し、超音波の送受信を制御する。例えば、送受信回路１２０は、素子Ｅに対して、送信する超音波の周波数や大きさ、波の種類（連続波やパルス波等）等を指示する。
【００３１】
スイッチ回路１１０は、リングアレイＲの複数の素子Ｅの各々に接続されており、送受信回路１２０からの信号を任意の素子Ｅに伝達し、素子Ｅを駆動させ、信号の送受信を行わせる。例えば、スイッチ回路１１０が、送受信回路１２０からの制御信号を供給する素子Ｅを切り替えることで、複数の素子Ｅのいずれか１つを、超音波を送信する送信素子として機能させ、複数（例えば全て）の素子Ｅで散乱波を受信させる。
【００３２】
全ての素子Ｅを同時に駆動して測定データを収集してもよいし、リングアレイＲの複数の素子Ｅをいくつかのグループに分けて、グループ単位で順に測定データを収集してもよい。グループの切り替えを数μｓ〜ｍｓオーダ以下で行うことで、ほぼリアルタイムに測定データの収集を行うことができる。
【００３３】
リングアレイＲは、ステッピングモータ等により上下動可能に設置されている。リングアレイＲを上下動させて、被検体Ｔの全体のデータ収集が行われる。
【００３４】
演算装置１３０は、例えばＣＰＵ、記憶部（ＲＡＭ、ＲＯＭ、ハードディスク等）、通信部等を備えたコンピュータにより構成されている。記憶部に格納されたプログラムが実行されることで、図３に示すような、送信素子決定部１３１、データ収集部１３２、計算部１３３、画像作成部１３４等の機能が実現され、行列データ格納領域１３５及び測定データ格納領域１３６が記憶部に確保される。各部による処理については後述する。
【００３５】
次に、本実施形態による散乱画像の作成方法について説明する。図４に示すように、１つの点散乱体ＰＳ（被検体Ｔの１点）に着目し、１つの送信素子ＥＴから送信された超音波が、この点散乱体ＰＳで散乱し、１つの受信素子ＥＲで受信される場合を考える。このとき、受信素子ＥＲでの測定データには、この点散乱体ＰＳの影響による固有パターンが含まれる。
【００３６】
本実施形態では、被検体Ｔが、多数の点散乱体で構成されていることを前提に、被検体Ｔを構成する特定の部分が超音波を反射し散乱させるものとする。図５に示すように、多数の点散乱体ＰＳが存在する場合、受信素子ＥＲでの測定データには、各点散乱体ＰＳの影響による固有パターンの線形結合が含まれる。本実施形態では、測定データから、各固有パターンを切り分けて判別し、各点散乱体ＰＳの散乱波の音圧の強度である散乱音圧強度を算出する。
【００３７】
本実施形態では、被検体の前部又は一部が含まれる撮像領域を分割した分割領域を設定する。具体的には、図６に示すように、撮像領域（関心領域）Ｒを格子状に分割し、複数のピクセル領域を設定する。例えば、撮像領域ＲをＭ２（＝Ｍ×Ｍ）個のピクセル領域に分割する。図６は、Ｍ＝１１とし、撮像領域Ｒを１２１個のピクセル領域Ｐ１〜Ｐ１２１に分割した例を示す。
【００３８】
本実施形態では、分割領域のそれぞれの散乱波の散乱音圧強度を、所定の素子から送信した超音波が分割領域における被検体で散乱して複数の素子の各々で受信されるまでの信号到達時間と、測定データに基づき、算出する。
【００３９】
各ピクセル領域の散乱音圧強度は１次元ベクトルｘに展開することができる。ピクセル領域がＭ２個ある場合、Ｍ２×１のベクトルｘとなる。ある特定の１つの送信素子ＥＴから超音波を送信し、Ｎ個の受信素子ＥＲを使用し、一素子あたりのサンプリング数をＮｔとした場合の測定データは、Ｎｔ×Ｎ×１のベクトルｙで表現できる。
【００４０】
ここで、適切な測定行列Ｇを使用すると、Ｇｘ＝ｙとして数式化することができる。Ｇの逆行列Ｇ−１を用いて、ｘ＝Ｇ−１ｙからｘが求まる。すなわち、逆行列Ｇ−１と、Ｎ個の受信素子ＥＲでの測定データとから、各ピクセル領域の散乱音圧強度を算出することができる。
【００４１】
なお、逆行列以外にも、いわゆる疑似逆行列を用いる方法も有効である。これはＧ−１Ｇ＝Ｅ（Ｅは単位行列）となるＧ−１以外にも、Ｇ’Ｇ＝Ｅ’、Ｅ’は定められた条件において、トレースの絶対値和が最小となるようなＧ’を用いる方法である。Ｇ’は、最小二乗法や拘束条件付きの変分法などにより求めることが出来る。この手法は、逆行列の発散に起因して処理後のノイズが大きくなることを防ぐ観点から優れている。同じように、Ｈ＝Ｇ＋λＥに対する逆行列Ｈ’を計算する手法なども有効である。
【００４２】
次に、測定行列Ｇの構造について説明する。送信素子ＥＴを固定し、各ピクセル領域を１つの点散乱体とみなし、ｉ番目の点散乱体で散乱された超音波を複数の受信素子で受信する場合の信号到達時間をｉ番目の列ベクトルとすることで、測定行列Ｇを構築する。
【００４３】
例えば、図７ａに示すように、送信素子ＥＴから送信した超音波が、１番目のピクセル領域Ｐ１の位置にある点散乱体で散乱し、複数（例えば全て）の受信素子で受信されるまでの信号到達時間を、受信素子の配置順に並べた列ベクトルｃ１が、測定行列Ｇの１番目の列ベクトルとなる。
【００４４】
図７ｂに示すように、送信素子ＥＴから送信した超音波が、２番目のピクセル領域Ｐ２の位置にある点散乱体で散乱し、複数の受信素子で受信されるまでの信号到達時間を、受信素子の配置順に並べた列ベクトルｃ２が、測定行列Ｇの２番目の列ベクトルとなる。
【００４５】
図７ｃに示すように、送信素子ＥＴから送信した超音波が、６１番目のピクセル領域Ｐ５６の位置にある点散乱体で散乱し、複数の受信素子で受信されるまでの信号到達時間を、受信素子の配置順に並べた列ベクトルｃ５６が、測定行列Ｇの６１番目の列ベクトルとなる。
【００４６】
図７ｄに示すように、送信素子ＥＴから送信した超音波が、１２１番目のピクセル領域Ｐ１２１の位置にある点散乱体で散乱し、複数の受信素子で受信されるまでの信号到達時間を、受信素子の配置順に並べた列ベクトルｃ１２１が、測定行列Ｇの１２１番目の列ベクトルとなる。
【００４７】
このようにして測定行列Ｇが構築される。図示しない計算機により測定行列Ｇの逆行列Ｇ−１が計算され、行列データ格納領域１３５に格納される。
【００４８】
送信素子決定部１３１は、この測定行列Ｇで送信素子として仮定した素子Ｅから超音波が送信されるように、送受信回路１２０に指示する。
【００４９】
データ収集部１３２は、スイッチ回路１１０及び送受信回路１２０を介して、複数の素子により得られたデータである測定データ（受信データ）を収集（受信又は取得することを含む）する。測定データは、測定データ格納領域１３６に格納される。
【００５０】
計算部１３３は、行列データ格納領域１３５に格納されている逆行列Ｇ−１と、測定データ格納領域１３６に格納されている測定データを成分とするベクトルｙとの積から、各ピクセル領域の散乱音圧強度ｘを算出する。
【００５１】
画像作成部１３４は、各ピクセル領域の散乱音圧強度を画素値に変換し、２次元アレイ状に配置することで、撮像領域の散乱画像（散乱音圧強度のマップ）を作成する。作成された散乱画像は、画像表示装置１４０に表示される。
【００５２】
逆行列Ｇ−１を求めることが可能となるのは、行列Ｇのランク（階数）がピクセル数Ｍ２と等しい必要がある。ランクとは、この行列がもつ異なる固有値の数と言い換えることもできる。逆にいうと、行列Ｇがランク＝Ｍ２を満たす条件で、求めるべきピクセル数と受信素子数の対応関係を定めて、その条件にてデータ取得を行うように構成する。予め計算したＧ-1に対して、取得データｙを乗算した時に、解が望ましくない場合（著しいアーチファクトが存在する場合など）には、データ取得後により小さなピクセル数（Ｍ２）２に対して構成したＧ２に対して計算されたＧ２-1にデータを乗算して画像再構成を行うように構成することも可能である。（Ｍ２＜Ｍ）
【００５３】
このように、本実施形態によれば、単一の送信素子から超音波を送信し、複数の受信素子で受信したエコーデータを並べたベクトルと、予め準備しておいた測定行列Ｇの逆行列Ｇ−１との積から、撮像領域（関心領域）の各ピクセル領域の散乱音圧強度を求め、散乱画像を作成する。
【００５４】
そのため、１素子から送信した超音波のエコー信号を全素子で受信することを、送信素子を切り替えながら繰り返し行う従来の開口合成法と比較して、被検体の計測に要する時間を短縮できる。また、取得するデータ量を削減し、データ転送に要する時間を短縮できる。
【００５５】
上記実施形態では、送信素子を１つだけとする例について説明したが、これは、ピクセル領域の数が少ない場合や、ノイズが極めて少ない場合に好適である。ピクセル領域の数が多い場合や、ノイズが大きい場合は、複数の送信素子を切り替えて順に超音波を送信し、測定データを収集することが好ましい。この場合、送信素子毎に逆行列Ｇ−１を予め準備しておく。逆行列Ｇ−１と測定データのベクトルとの積を計算し、送信素子数分の散乱音圧強度ｘを求める。複数の散乱音圧強度ｘを合成することで、信号対雑音比を改善することができる。
【００５６】
上記実施形態において逆行列Ｇ−１が解けない場合は、最小二乗法や、ペナルティタームを正則化法により解くことで、散乱音圧強度ｘを求めることができる。特にノイズｎの影響が無視できない場合には、Ｇｘ＋ｎ＝ｙとなり、ｘ＝Ｇ-1（ｙ-ｎ）となるが、一般的にはｎを特定できないためである。
【００５７】
シミュレーション
２種類のファントムに対し、上記実施形態による散乱画像作成方法を適用したシミュレーションを行った。また、比較例として、開口合成法を適用したシミュレーションを行った。ファントムは、図８ａ、図１０ａに示す格子状に配置された４９点（７×７）の点散乱体と、図９ａ、図１０ｄに示す６４×６４ピクセルのＳｈｅｐｐ−Ｌｏｇａｎファントムとした。シミュレーション条件を以下の表１に示す。
【００５８】
【表１】
2019082892.tif 000003
【００５９】
比較例による開口合成法を適用した結果を図８、図９に示す。図８は４９点の離散的な散乱体をモデル化した場合を示し、図９は構造をモデル化した場合を示す。図８ａ、図９ａが正解（元のモデル）を示す。
【００６０】
図８ｂ〜図８ｆ、図９ｂ〜図９ｆは、それぞれ、送受信素子数を１２８個、６４個、３２個、１６個、８個とした場合の開口合成撮像の結果を示す。図８、図９に示すように、送受信素子数の減少に伴い、本来輝度が存在しない画素においてノイズが生じていることが確認できた。特に図８では、送受信素子数が１２８個、６４個、３２個の場合はノイズがランダムに分布しているのみであるが、送受信素子数が１６個、８個の場合、特定のパターンとしてアーチファクトを形成しており、誤診等の原因となりうる。
【００６１】
図９では、送受信素子数１２８個のときから、腫瘍と乳腺の分離が困難となるようなコントラストの低下が生じており、送受信素子数が３２個、１６個、８個の画像においては、様々なアーチファクトの発生により画像の視認性が著しく低下した。
【００６２】
一方、図１０は本実施形態による方法を適用した結果であり、図１０ａ、図１０ｄが正解（元のモデル）を示し、図１０ｂ、図１０ｅが送受信素子数を８個とした場合の復元画像を示す。図１０から、元のモデルを完全に復元できていることが確認された。
【００６３】
送信素子数が増えた場合は、送信条件数の数だけ部分行列を縦にならべて行列Ｇを形成した。他にも、送信素子数を１個、受信素子数を１６個とした場合においても、同様の結果が得られることが確認された。更に、送信素子数を１個とし、受信素子数を３２個、６４個、１２８個とした場合も、受信素子数が１６個の場合と同様な再構成結果が得られることを確認した。
【００６４】
これらの結果から、
［１］リングアレイを構成する素子数を維持したまま、送信条件を１に削減可能
［２］リングアレイを構成する素子数を１２８から８に削減し、送信条件は素子数の分だけ実施する
という２つの使途があることが確認できる。
【００６５】
［１］の場合は、開口合成法と比較して、撮像速度が１２８倍、取得データ量は１／１２８になる。上記の例では素子数を１２８個としているが、例えば２０４８素子、２５６送信条件のシーケンスでは、撮像速度が２５６倍、取得データ量は１／２５６になる。
【００６６】
［２］の場合は、リングアレイを構成する素子数の削減や、マルチプレクサの回路規模低減、素子とマルチプレクサ間のケーブル本数の低減（いずれも１／１６）というコスト削減や装置サイズ縮小の効果がある。
【００６７】
一断面の撮像に要する時間は、直径×２／音速×撮像条件数となるので、直径２００ｍｍ、音速１５００ｍ／ｓ、代表的な開口合成法の撮像条件数２５６回の場合、一断面の撮像に約７０ｍｓかかる。１０００断面撮像する場合、モータでの移動速度が十分に小さいと、７０秒程度要することになる。本発明の方法を適用すると、１０００断面撮像しても、撮像時間は０．２３秒となる。このように、撮像が高速化され、データ量も削減された分、撮像枚数を増やすことに自由度を使うことも可能となる。
【００６８】
データ量に関しては、素子数２０４８、送信条件２５６、直径２０ｃｍのリングでサンプリング周波数が４０ＭＨｚ、ＡＤ変換器の出力が２Ｂｙｔｅの開口合成法の場合、深さ方向（超音波伝搬方法）のサンプル点数が２００ｅ−３×２／１５００×４０ｅ６＝１ｅ４であり、１フレーム当たり、２×１ｅ４×２０４８×２５６＝１０ＧＢとなる。断面数が１０００の場合は、１つのボリュームデータあたり１０ＴＢとなる。このため、断面数を削減するか、エコーデータではなく、画像データ（約２ＧＢ／体積）で保存することを余儀なくされる。１つのエコーデータに多様なアプリケーション処理を加えることで診断能力の向上が期待されるが、毎回１０ＴＢのデータを保存することは現実的ではなくなってしまっている。一方、本発明により、大よそ２桁から３桁でデータ削減が可能となることにより、エコーデータ自体の保存も可能となり、画像化する前のエコーデータの状態で過去データとの比較が可能となるなど、新しい使い道の発展にもつながる。
【００６９】
ボリュームでの撮像が高速化されると、乳腺撮像の主要なアプリケーションの一つであるエラストグラフィにおいても、大きな改善が可能となる。エラストグラフィでは一断面の１ライン上での加圧前後での相互相関によって歪を抽出、その分布を可視化することで病変を検出する技術である。加圧により、ラインやスライスがずれると相互相関精度が低下する。通常は撮像スライス位置を固定して、ラインのずれのみ相互相関対象に含めることで、相関エラーを提言しているが、ボリュームでの撮像が高速化されれば、ボリューム内でのエコーライン間の相関が可能とあり、相関エラーの低減が可能となる。
【００７０】
ここまでの説明では、Ｇ行列を生成する際にインパルス応答を前提とした説明を行った。実際にはトランスデューサは共振周波数を中心としたバンドパスフィルタであるので、帯域幅は有限である。この場合においても、Ｇ行列の生成時にインパルス応答に対応した波形を繰り込むことで演算が可能となる。
【００７１】
また、ここまでの説明では、音速不均質が与える影響に関して議論を行っていない。すなわち、音波の音速は、乳房の内部を透過する場合と、そうでない場合、また１つの経路中で乳房内部を透過する部分と外部を透過する部分の比の変化とで不均一となることから、両音波の音速の違いを踏まえることが好ましい。例えば、散乱画像を作成する前に乳房の内外を検出する二値化処理を行い、積分伝搬時間から乳房内の平均音速を算出することが可能である。この音速分布を活用することで、音速不均質に対応したＧ行列の修正をすることは、再構成アルゴリズムのロバスト性向上に有効である。
【００７２】
Ｇ行列の生成における留意点としては、透過波の分離の観点もある。送信に用いる周波数によるが、一般的には散乱波より透過波の方が強いことが多い。図１１ｂに示すように、透過波が散乱エコー信号の上に重なっていると、本発明の方法により再構成を行う際に影響を受け得る。そこで、図１１ａに示すように、被検体（撮像領域）からみて送信素子側に配置された一部の素子のみを用いて受信を行うように受信開口制限を行うことで、透過波の影響を限定的にすることが可能となる。図示した例では、受信開口を全素子の３／８に制限することで、本発明の手法が実現可能となることが確認されている。
【００７３】
ここまでの実施例ではリング内の広い領域で等間隔に画素を設定する方法に関して説明を行った。リング状のアレイを用いた乳腺組織の撮像においては、乳房内の領域と、それを取り囲む水の領域に大別される。当然、計測においては水の領域に関しては撮像を行う必要がない。そこで撮像領域を乳房が存在する領域のみに設定することで、画素数が削減され、それに応じて、必要な送信条件数、サンプリング周波数を低下させることが、本発明によって可能となる。リングアレイが体幹部に近い（ベッドの天板に近い）場合では、水の領域でデータを取らないことによるデータ削減効果は少ないが、体幹部から離れ、乳房の先端に近づいていくと撮像面内での乳房断面積は低下し、本発明で最適化された撮像条件における取得データ量の削減効果は大きくなる。この時に、水の中に気泡などが存在することによって、ノイズが存在すると、エラーとなる可能性があるので、一条件に関して、時間を変えて撮像を行い、経時変化しない成分のみを抽出する前処理を適用することで、水中に浮遊する散乱体に起因するエコー信号の除去を行うことがノイズ低減のために有効となる。
【００７４】
次に本発明が成立する条件、主要なパラメータの設定方法に関して、補足的な説明を行う。本発明のポイントは、計測エリアをグリッド化することで、離散化して、離散的なモデルとして取り扱うために、行列での表現が可能とすることである。この時、グリッドサイズが重要である。散乱波を取得して、画像化する従来法の代表例である開口合成法と比較して、グリッドサイズに対して画質がどのように変化するかを説明する。
【００７５】
図１３に２つの評価用モデルを示す。図１３ａは点の配列からなる構造を示し、図１３ｂは乳房断層像を模擬した構造を示す。中心周波数２ＭＨｚ、空間を画素サイズ０．２ｍｍ（波長の約８分の１）でグリッドした時に、グリッド内に存在する散乱体の散乱量の総和をグリッド中心に存在するものとして取り扱う。この時、グリッド内の散乱体の空間分布を、グリッド中央に存在するδ関数に置き換える近似操作が、画像再構成プロセスにおけるアーチファクトの生成や、解像度の低下の原因となる。典型的な結果として、図１４ａ、図１４ｂに本発明による結果、図１５ａ、図１５ｂに比較例として公知の開口合成法を用いた結果を示す。図１４ａ、図１４ｂでは、解像度は維持されているが、ノイズの増加が顕著である。一方、図１５ａ、図１５ｂでは、ノイズの増加は顕著ではないが、解像度の低下が顕著である。
【００７６】
そこで、グリッド内の散乱体位置誤差（設定位置と、グリッド中央の距離）を横軸にして、２つの画質評価因子に関する解析を行った。すなわち、信号対雑音比（ＳＮＲ）と、解像度（半値幅）を縦軸にとり、開口合成を用いた場合と、本発明の場合に関して、図１６ａ、図１６ｂにプロットした。結果として、開口合成の場合は、誤差が増えると解像度は低下する（分解能＝分解可能な下限サイズ［ｍ］が増大する）が、ＳＮＲには大きな変化がない。一方、本発明では解像度の変化は小さいがＳＮＲの低下が大きい。
【００７７】
解像度の低下は、画像がぼやけるが、観察対象が消失するわけではない。（もちろん、程度による。）一方、ＳＮＲがある値を下回ると、そもそも観察対象を視認することが出来なくなる。従来、逆行列の演算が実用的に計算可能な範囲で考えると、グリッドサイズが大きくなり、位置誤差が増え、ＳＮＲが急激に劣化するため、本発明のような手法が検討されることは無かった。本発明は、十分なグリッド数があれば対象を離散的に取り扱えるという着想から、発明に至ったものである。この結果で示すように、本発明が効果を発揮するのは、開口合成法に比べると、より厳しい条件となっているが、効果を発揮する条件においては、より高い性能を実現することが可能となる。
【００７８】
また、別の実施例として、乳房内の温度変化計測に本発明を用いる例について説明を行う。がん細胞は他の正常組織を構成する細胞に比べ、増殖速度が速く、早い増殖を実現するために代謝が活発であることが知られている。このため、既存のＰＥＴ（ポジトロン断層法）では、放射性崩壊する同位体でラベリングした糖を投与して、代謝が活発な部位から放射されるγ線源の空間分布を可視化することによって、代謝が活発な部位＝腫瘤を検出する診断法が臨床現場で広く用いられている。高いコントラストで病変検出が可能な技術であるが、放射性同位体薬剤を生成するための加速器などの付帯設備が高額、大型であることや、体内被曝を伴うことが検診での利用を妨げている。
【００７９】
代謝計測をより簡易に行う方法として、赤外線カメラを用いて、代謝に伴う温度上昇部位を検出する手法に関しても、検討の歴史は長く、かつては臨床装置も販売されていた。（近年、ＡＩ技術の発展に伴い、再度開発を検討するグループもある。）赤外線で計測する場合には、生体中の水分が熱輻射に伴う赤外線を吸収してしまうため、計測可能な領域は体表に近い浅部に限られることが課題であった。
【００８０】
リングアレイを用いた計測では、温度変化に伴う音速分布の変化を計測することが可能である。実際に水分や脂肪の音速の温度依存性は、Ｘｍ／ｓ／ｋ程度であり、リングアレイでも検出可能である。体内深部も原理的には計測可能である点が赤外線カメラを用いた場合に比べた長所である。ただし、音速の温度依存性から計測可能なのは、絶対温度ではなく、温度変化であるので、乳房に温度変化を与え、代謝量が多い部位と、他の部位で、温度変化速度の違いからコントラストを得るのが有効である。本発明によらない、通常のリングエコー撮像のシーケンスでは、１ボリュームのデータを取得するのに、５分から１０分の時間を要する。腫瘤で発生した熱は、熱拡散方程式に従い拡散し、平衡状態においては、熱源以外の周囲の組織との温度差は小さくなってしまう。この平衡状態に達するまでの温度変化を効率的に取得し、かつ付加的な設備を設けずに温度変化を与えるには、リングアレイが格納された水槽中に、被検者が乳房を入れた直後の温度の非平衡過程を観察することが望ましい。
【００８１】
この目的おいて、本発明による、送信条件数の削減は効果的である。上記の１ボリュームのデータ取得時間５分に対して、過渡的な温度変化時間が５〜１０分程度であり、この間に過渡的な温度変化を観察するには、撮像速度を５〜１０倍高速化する必要がある。本発明による高速化手法を用い、散乱画像を所定時間毎に連続して作成し、経時的な画像変化から、変化の大きい領域（画素値や対応するＲＦデータが所定値以上変化する部分）を抽出し、それを温度変化の大小に換算することで、代謝による熱抵抗の相違を可視化することが可能となる。
【００８２】
本発明を特定の態様を用いて詳細に説明したが、本発明の意図と範囲を離れることなく様々な変更が可能であることは当業者に明らかである。
本出願は、２０１７年１０月２４日付で出願された日本特許出願２０１７−２０５３４３に基づいており、その全体が引用により援用される。
【符号の説明】
【００８３】
１０  超音波診断システム
１１０  スイッチ回路
１２０  送受信回路
１３０  演算装置
１４０  画像表示装置
【図１】
2019082892.tif 000004
【図２】
2019082892.tif 000005
【図３】
2019082892.tif 000006
【図４】
2019082892.tif 000007
【図５】
2019082892.tif 000008
【図６】
2019082892.tif 000009
【図７】
2019082892.tif 000010
【図８】
2019082892.tif 000011
【図９】
2019082892.tif 000012
【図１０】
2019082892.tif 000013
【図１１】
2019082892.tif 000014
【図１２】
2019082892.tif 000015
【図１３】
2019082892.tif 000016
【図１４】
2019082892.tif 000017
【図１５】
2019082892.tif 000018
【図１６】
2019082892.tif 000019
【国際調査報告】
2019082892.tif 000020
2019082892.tif 000021