IP Force 特許公報掲載プロジェクト 2022.1.31 β版

知財求人 - 知財ポータルサイト「IP Force」

▶ 一般社団法人メディカル・イメージング・コンソーシアムの特許一覧 ▶ 株式会社MU研究所の特許一覧 ▶ 国立大学法人富山大学の特許一覧

<>
  • 特許-超音波医用システム 図1
  • 特許-超音波医用システム 図2
  • 特許-超音波医用システム 図3
  • 特許-超音波医用システム 図4
  • 特許-超音波医用システム 図5
  • 特許-超音波医用システム 図6
  • 特許-超音波医用システム 図7
  • 特許-超音波医用システム 図8
  • 特許-超音波医用システム 図9
  • 特許-超音波医用システム 図10
  • 特許-超音波医用システム 図11
  • 特許-超音波医用システム 図12
  • 特許-超音波医用システム 図13
  • 特許-超音波医用システム 図14
< >
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2023-08-24
(45)【発行日】2023-09-01
(54)【発明の名称】超音波医用システム
(51)【国際特許分類】
   A61B 8/14 20060101AFI20230825BHJP
【FI】
A61B8/14
【請求項の数】 3
(21)【出願番号】P 2019192437
(22)【出願日】2019-10-23
(65)【公開番号】P2021065395
(43)【公開日】2021-04-30
【審査請求日】2022-10-19
【新規性喪失の例外の表示】特許法第30条第2項適用 平成30年10月29日~31日に、第39回超音波エレクトロニクスの基礎と応用に関するシンポジウムで、本発明の一部である超音波断層像生成方法、超音波断層像生成装置について部分的に発表した。 令和1年6月11日に、Japanese Journal of Applied Physics 58,SGGE06(2019)で、本発明の一部である超音波断層像生成方法、超音波断層像生成装置について部分的に発表した。
(73)【特許権者】
【識別番号】514189468
【氏名又は名称】一般社団法人メディカル・イノベーション・コンソーシアム
(73)【特許権者】
【識別番号】513004353
【氏名又は名称】株式会社MU研究所
(73)【特許権者】
【識別番号】305060567
【氏名又は名称】国立大学法人富山大学
(74)【代理人】
【識別番号】100126000
【弁理士】
【氏名又は名称】岩池 満
(74)【代理人】
【識別番号】100154748
【弁理士】
【氏名又は名称】菅沼 和弘
(74)【代理人】
【識別番号】110002055
【氏名又は名称】弁理士法人iRify国際特許事務所
(72)【発明者】
【氏名】千葉 敏雄
(72)【発明者】
【氏名】山下 紘正
(72)【発明者】
【氏名】望月 剛
(72)【発明者】
【氏名】長谷川 英之
【審査官】冨永 昌彦
(56)【参考文献】
【文献】特開2006-204923(JP,A)
【文献】国際公開第2018/230066(WO,A1)
【文献】特開2015-164516(JP,A)
【文献】特開2015-071028(JP,A)
【文献】米国特許出願公開第2015/0063058(US,A1)
【文献】国際公開第2018/206736(WO,A1)
【文献】国際公開第2007/097108(WO,A1)
(58)【調査した分野】(Int.Cl.,DB名)
A61B 8/00 - 8/15
(57)【特許請求の範囲】
【請求項1】
表示装置としての8Kモニタと、
超音波断層像生成装置と、を備え、
前記超音波断層像生成装置は、
M(2以上の自然数)個の超音波振動子から発せられた超音波のエコーを受信してエコー信号を出力する超音波プローブから出力されるMチャネルのエコー信号における雑音を推定し、受信焦点からのエコーを強調する重み係数を前記Mチャネルのエコー信号における信号対雑音比に応じて算出する推定手段と、
超音波断層像を表すビームフォーマを、前記推定手段にて算出された重み係数を用いて前記Mチャネルのエコー信号から生成する生成手段と、
を含むことを特徴とする超音波医用システム。
【請求項2】
前記推定手段は、
前記Mチャネルのエコー信号に遅延和ビームフォーミングにおける遅延を付与して得られるm番目のエコー信号smまでを累算して得られる累積要素信号umと、エコー信号smに含まれる直流成分yと付加的ノイズに起因するバイアスnとを用いてモデル化したモデル化要素信号Um=m×y+nと、の平均二乗差αが最小になるようにyおよびnの値を設定し、設定したyおよびnを用いて前記平均二乗差αの最小値を算出し、当該最小値と前記設定したyとから前記重み係数を算出し、
前記生成手段は、
前記設定したyに前記重み係数を乗算して、超音波断層像を表すビームフォーマを生成する
ことを特徴とする請求項1に記載の超音波医用システム。
【請求項3】
前記推定手段は、
前記Mチャネルのエコー信号に遅延和ビームフォーミングにおける遅延を付与して得られるm番目のエコー信号smまでに含まれる雑音の積分値nmの二乗平均と、遅延補償後のエコー信号smの平均Y DAS の二乗平均とから前記重み係数を算出し、
前記生成手段は、
前記平均Y DAS に前記重み係数を乗算して、超音波断層像を表すビームフォーマを生成する
ことを特徴とする請求項に記載の超音波医用システム
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、超音波断層像を生成し表示する技術に関する。
【背景技術】
【0002】
超音波診断は、被検者の体内の断層像を非侵襲的に測定する手法であり,広く医療現場に普及している。超音波診断により得られる超音波断層像の空間分解能およびコントラストは,診断精度に直結する重要な要素である。このため、超音波断層像の空間分解能およびコントラストを向上させる技術が種々提案されており、その一例としては非特許文献1に開示の技術が挙げられる。非特許文献1に開示の技術では、複数の超音波振動子からなる配列型超音波振動子により受信された各エコー信号間の相関性(coherence)に基づき,超音波断層像の方位分解能およびコントラストを向上させている。
【先行技術文献】
【非特許文献】
【0003】
【文献】超音波受信信号間の相関性に基づく超音波イメージング手法 P.-C. Li and M.-L. Li, “Adaptive imaging using the generalized coherence factor,” IEEE Trans. Ultrason. Ferroelectr. Freq. Control, vol. 50, no. 2, pp. 128-141, 2003.
【文献】H. Hasegawa and H. Kanai, “Effect of element directivity on adaptive beamforming applied to high-frame-rate ultrasound,” IEEE Trans. Ultrason. Ferroelectr. Freq. Control, vol. 62, no. 3, pp. 511-523, 2015.
【発明の概要】
【発明が解決しようとする課題】
【0004】
前述したように、超音波断層像の空間分解能およびコントラストは,診断精度に直結するため、高ければ高い程好ましい。さらに、空間分解能およびコントラストの向上をより顕在化させることが嘱望されている。
【0005】
本発明は、上述した事情に鑑みてなされたものであり、超音波断層像の空間分解能およびコントラストを向上させ、8Kモニタの採用により、その向上による効果をより顕在化させることを目的とする。
【課題を解決するための手段】
【0006】
上記課題を解決するために、本発明の第1の態様に係る超音波医用システムでは、表示装置としての8Kモニタと、超音波断層像生成装置と、を備え、前記超音波断層像生成装置は、送信方式として、送信時に各振動子への励振パルスの位相を制御し、集束型の送信ビームを形成する集束型送信、又は送信波の位相を制御せずに同時に又は送信ビームが拡散するように位相制御する非集束型送信のいずれかを採用し、受信方式として、M(2以上の自然数)個の超音波振動子からの超音波のエコーを受信してエコー信号を出力する超音波プローブから出力されるMチャネルのエコー信号の位相を制御する適合型ビームフォーミングを採用した。
【0007】
本発明の第2の態様に係る超音波医用システムは、表示装置としての8Kモニタと、超音波断層像生成装置と、を備え、前記超音波断層像生成装置は、M(2以上の自然数)個の超音波振動子から発せられた超音波のエコーを受信してエコー信号を出力する超音波プローブから出力されるMチャネルのエコー信号における雑音を推定し、受信焦点からのエコーを強調する重み係数を前記Mチャネルのエコー信号における信号対雑音比に応じて算出する推定手段と、超音波断層像を表すビームフォーマを、前記推定手段にて算出された重み係数を用いて前記Mチャネルのエコー信号から生成する生成手段と、を含む。
【0008】
ここで、前記推定手段は、前記Mチャネルのエコー信号に遅延和ビームフォーミングにおける遅延を付与して得られるm番目のエコー信号sまでを累算して得られる累積要素信号uと、エコー信号sに含まれる直流成分yと付加的ノイズに起因するバイアスnとを用いてモデル化したモデル化要素信号U=m×y+nと、の平均二乗差αが最小になるようにyおよびnの値を設定し、設定したyおよびnを用いて前記平均二乗差αの最小値を算出し、当該最小値と前記設定したyとから前記重み係数を算出し、前記生成手段は、前記設定したyに前記重み係数を乗算して、超音波断層像を表すビームフォーマを生成するようにしてよい。
【0009】
さらに、前記推定手段は、前記Mチャネルのエコー信号に遅延和ビームフォーミングにおける遅延を付与して得られるm番目のエコー信号sまでに含まれる雑音の積分値nの二乗平均と、遅延補償後のエコー信号sの平均YDASの二乗平均とから前記重み係数を算出し、前記生成手段は、前記平均YDASに前記重み係数を乗算して、超音波断層像を表すビームフォーマを生成するようにしてよい。
【発明の効果】
【0010】
本発明によれば、超音波断層像の空間分解能およびコントラストを向上させ、8Kモニタの採用により、その向上による効果をより顕在化させる技術を提供することができる。
【図面の簡単な説明】
【0011】
図1】本発明の一実施形態による超音波断層像生成装置を含む超音波医用システムの構成例を示すブロック図である。
図2】各方式に対する3つの分解能に係る当該8Kモニタの有用性を示す図である。
図3図3(a)には集束型ビームフォーミング、図3(b)には非集束型ビームフォーミングについて、超音波ビームの様子を示す図である。
図4】受信信号から超音波ビームラインを形成し、その情報を従来モニタ(1K相当)と8Kモニタに表示したときの様子を示す図である。
図5】従来モニタでの表示能力について説明する図である。
図6】8Kモニタでの利点を説明する図である。
図7】理想的な画像(画素が無限小で表現される画像)の枠を示す図である。
図8】フレーム補間の説明図である。
図9】各時の画像の移動の見え方について、従来の1Kモニタと8Kモニタの比較を示す図である。
図10】フレーム内のある画素の輝度が時間とともに変化する様子を示す図である。
図11】画素の大きさによる動きの滑らかさに加えて輝度の量子化の細かさを含めて滑らかな変化となる様子を示す図である。
図12】超音波断層像生成装置の信号処理部が実行する信号処理の流れを示すフローチャートである。
図13】超音波断層像の空間分解能を評価するための点ターゲットをイメージングした図である。
図14】超音波断層像のコントラストを評価するためのファントムを画像化した図である。
【発明を実施するための形態】
【0012】
以下、図面を参照しつつ本発明の実施形態を説明する。
(A.実施形態)
図1は、本発明の一実施形態によれる超音波断層像生成装置20を含む超音波医用システム1の構成例を示す図である。超音波医用システム1は、医療現場において被検者の体内の超音波断層像を非侵襲的に撮像するためのシステムである。図1に示すように、超音波医用システム1は、超音波断層像生成装置20の他に、信号線を介して各々超音波断層像生成装置20に接続された超音波プローブ10、操作装置30および表示装置40、を有する。
【0013】
超音波プローブ10は、複数の超音波振動子よりなる配列型超音波振動子を有する。本実施形態の超音波医用システム1では、超音波プローブ10として、M(2以上の自然数)個の超音波振動子を0.1mm間隔で配列したリニアアレイプローブ(PU-0558:上田日本無線株式会社)が用いられている。複数の超音波素子の各々は、超音波断層像生成装置20による制御の下、被検者の検査部位に向けて超音波を放射するとともに、当該超音波のエコーを受信してエコー信号を出力する。
【0014】
超音波断層像生成装置20は、超音波プローブ10に超音波を送信させるとともに、超音波プローブ10からの出力信号に信号処理を施して画像データを生成する。操作装置30は、例えばマウスなどのポインティングデバイスやキーボードを含む。操作装置30は、超音波医用システム1の利用者(例えば、超音波診断のための各種操作を行う検査技師)に超音波断層像生成装置20に対する各種入力操作を行わせるための装置である。表示装置40は、例えば液晶ディスプレイである。表示装置40は、超音波断層像生成装置20の出力する画像データに応じて画像を表示する。
【0015】
超音波断層像生成装置20は、図1に示すように、制御部200と、送信部210と、受信部220と、信号処理部230と、を有する。図1では詳細な図示を省略したが、超音波断層像生成装置20は、OS(Operating System)などの各種ソフトウェアを記憶した記憶部(例えば、ハードディスク)も有する。
【0016】
制御部200は例えばCPU(Central Processing Unit)である。制御部200は、上記記憶部に記憶されているソフトウェアを実行することにより、超音波断層像生成装置20の制御中枢として機能し、各部の作動制御を行う。より詳細に説明すると、制御部200は、従来と同様のライン毎の取得シーケンスによって超音波断層像が生成されるように、各部の作動制御を行う。
【0017】
送信部210および受信部220には、信号線を介して超音波プローブ10が接続されている。送信部210は、制御部200から与えられる送信データにD/A変換を施して送信信号を生成し、超音波プローブ10が備えるM個の超音波振動子の各々に与える。これにより、超音波プローブ10が備えるM個の超音波振動子の各々は超音波を放射する。受信部220は、超音波プローブ10の複数の超音波振動子の各々から出力されるエコー信号にA/D変換を施し、さらに遅延を付与して遅延補償し、信号処理部230に与える。なお、本実施形態において、受信部220がエコー信号に付与する遅延は、従来の超音波断層像生成手法である遅延和ビームフォーミング(以下、DASビームフォーミング)に準拠した遅延である。
【0018】
超音波プローブ10のm(m=0、1、2・・・M-1)番目の超音波振動子から出力され、受信部220により遅延を付与されたエコー信号をsとすると、超音波プローブ10の受信開口に含まれるM個の超音波振動子により得られるエコー信号は以下の数1に示すベクトルSで表される。遅延補償の後、ベクトルSに含まれる受信焦点からのエコーは受信開口を横切る直流(DC)成分となる。したがって、従来のDASビームフォーミングでは、受信焦点yからのエコーに相当するビームフォーマ(すなわち、超音波断層像を表すビームフォーマ)YDASは、遅延補償後のエコー信号sの平均として以下の数2のように求められていた。
【数1】
【数2】
【0019】
これに対して、信号処理部230は、遅延補償後のエコー信号sに対して、本実施形態の特徴を顕著に示す信号処理(信号対雑音比に基づくビームフォーミング処理)を施して超音波断層像を表すビームフォーマを生成して表示装置40に与える。信号処理部230は、例えばDSP(Digital Signal Processor)であり、図1では詳細な図示を省略したが、信号処理部230には、信号対雑音比に基づくビームフォーミング処理を当該信号処理部230に実行させる信号処理プログラムが予めインストールされている。信号処理部230は、受信部220により遅延を付与された信号に対して、上記信号処理プログラムにしたがって、信号対雑音比ビームフォーミングまたは線形回帰ビームフォーミングを実行する。信号処理部230により実行されるビームフォーミング処理としては、従来技術に係る整相加算受信ビームフォーミングに加えて、本実施形態特有の適応型受信ビームフォーミング、即ち、信号対雑音比ビームフォーミングと線形回帰ビームフォーミングの2種類がある。これら信号対雑音比ビームフォーミングおよび線形回帰ビームフォーミングは、広義にはどちらも信号対雑音比に基づくものであるが、2つの手法を区別するために異なる名称とした。
【0020】
また、本実施形態に係る超音波断層像生成装置20では、送信方式については、例えば、集束型送信と非集束型送信のいずれにも対応可能である。集束型送信とは、送信時に各振動子への励振パルスの位相を制御し、集束型の送信ビームを形成する方法をいう。また、整相加算受信ビームフォーミングとは、各振動子から受信したエコー信号の位相を制御して受信ビームを形成する方法をいう。一方、非集束型送信とは、送信波の位相を制御せずに同時にまたは、送信ビームが拡散するように位相制御する送信方式をいう。適合型ビームファーミングとは、受信信号が各時刻で代々エネルギーとなるように各振動子からのエコー信号の位相を制御する詳細は後述する方法である。
【0021】
本実施形態では、表示装置40として、8Kモニタを採用するが、各方式に対する3つの分解能に係る当該8Kモニタの有用性は、図2にまとめられている。
【0022】
図2の中の数値は、装置の設定条件に応じて変化する範囲を示す。また分解能の仮定として、集束型送信時の分解能は非集束型のそれと比較して2倍高いとする。また、適応型受信ビームフォーミングでは、整相加算受信ビームフォーミングのそれと比較して2倍良くなると仮定する。一般に、分解能の具体的な数値は、超音波周波数、超音波ビームの走査方法などのパラメータにより変化するが、1)集束型+整相加算受信ビームフォーミングの方式では、空間分解が0.1~3mmと仮定し、コントラスト分解能は128送受信チャンネルを仮定して48dBとする。一方、2)非集束型送信+整相加算受信ビームフォーミングの方式では、受信時のみでのビームフォーミングであるので、空間分解能は1)と比較して2倍悪くなるので0.2~6mmと仮定した。さらに、1)2)の整相加算受信ビームフォーマと3)4)の適応型受信ビームフォーマとの比較では、空間分解能が適応型処理をすることで2倍改善されると仮定し、即ち、1)の0.1~3mmが3)では0.05~1.5mmとなる。これらのエコーの分解能を現行の表示系と8Kモニタとを比較し、エコー分解能を十分に表示できているかそれとも不十分であるかを図2は示している。
【0023】
図3(a)には集束型ビームフォーミング、図3(b)には非集束型ビームフォーミングについて、超音波ビームの様子が示されている。
【0024】
前者では、音波の波面が焦点で集束するように超音波パルスの位相が制御されて送信されている。後者では、超音波パルスが同時刻で送信されているために、超音波の波面は平面波となっている。もしくは、送信パルスの位相を制御して、集束とは反対に超音波の波面が円弧状に拡散する場合も含む。受信時には両者ともプローブ内の各アレイ振動子で同時に受信し、受信信号をメモリに記録する。それら記録された受信信号から受信ビームを形成することは共通である。
【0025】
図4には受信信号から超音波ビームラインを形成し、その情報を従来モニタ(1K相当)と8Kモニタに表示したときの様子を示す。
【0026】
従来モニタでのピクセルサイズは図2に示したように0.1~0.6mmと仮定する。この範囲は超音波画像をモニタに表示する際の診断距離(診断装置ではMAG機能に当たる)により変化する。即ち、モニタの200ピクセルに超音波診断距離20mmで表示すると1ピクセル当たり0.1mmとなるが、診断距離12cmを表示する場合には1ピクセル当たり0.6mmとなる。この場合に、8Kモニタではピクセルのサイズが1Kの8分の1となるので、0.0125~0.075mmとなる。
【0027】
図4(a)では従来モニタでの表示時の状態を、図4(b)では8Kモニタでの表示時の状態を示している。従来モニタの1)、4)では超音波ビームの実寸と表示系のピクセルサイズが同じであり、また2)ではビーム幅よりピクセルサイズが小さいためにそれらのビームで検出された反射体からのエコー信号がモニタに表示された画像(ピクセルにより表示された黒点)はビーム幅と等しい。それに対して、3)の場合は超音波ビーム幅がピクセルサイズよりも小さいので、画像は実際のビーム幅よりも大きく表示される。一方、8Kモニタでの表示の場合には表示系のピクセルサイズすべての場合において実際の超音波ビーム幅よりも十分に小さいので、それらのビームから得られた反射体の画像は実際のビームと等しくなる。
【0028】
次に、図5を参照して、従来モニタでの表示能力について説明を行う。
【0029】
生体から得られるエコー強度のダイナミックレンジは70dB近くあることが知られている。つまり、一番小さなエコーから最大のエコーの比が70dB近くなる。一方、従来モニタ明るさの表示能力は暗い輝度から最も明るい輝度までのダイナミックレンジが48dBである。このため、同時に70dBのエコー画像を表示することができない。仮に最大のエコー画像をモニタの最大輝度で表示したとすると、図5(a)の上段に示すように強い輝度にエコーは表示できるが、小さなエコーは画面には表示されない。
【0030】
逆に、図5(a)の下段に示すように最小のエコーをモニタの最小輝度で表示すると、48dBを超える強度のエコー画像はモニタの最大輝度で飽和した画像となってしまう。このような表示方法は、かなり以前にはリニア表示として実用化されていた時代もあったものの、診断としは同時に最大エコーと最小エコーを観察する必要があり、この表示では不十分であった。
【0031】
そこで、採用された方法がモニタの表示ダイナミックレンジに合うように、エコー信号を一度対数圧縮してモニタのダイナミックレンジ内に収める方法であり、この方法は対数圧縮表示、またはログ圧縮表示と呼ばれている(図5(b)参照)。
【0032】
この方法では、同時にすべてのエコー輝度情報を表示できる利点があるが、高輝度のエコー同士のわずかな違いを表現できない問題点があった。
【0033】
以下、図6を参照して、その問題点と8Kモニタでの利点を説明する。
【0034】
図6(a)はリニア表示の場合であり、図5(a)、図5(B)に対応する。この場合ではエコー輝度の低いときのモニタ表示の輝度差ΔAと高輝度での輝度差ΔBとは等しくなる。ただし、欠点として上記のようにダイナミックレンジが狭くエコー情報をカバーできない、それに対して、図6(b)のログ表示の場合には、ダイナミックレンジはカバーできるものの、それらは等しくならずΔA>ΔBとなる。ところが、8Kモニタの場合には、ダイナミックレンジが広くリニア表示でカバーできるので、エコー情報を漏れなく表示でき、かつ、エコー輝度の量子化は均等となる。つまり、ΔA=ΔBとなる。
【0035】
なお、以上の説明では、生体のエコー情報のダイナミックレンジが70dBであるとして説明したが、実際の装置で得られるエコー情報のダイナミックレンジは図2で分類したように、送信の方法や受信ビームフォーミングの方式により異なる。このために、それらの方式で生体情報70dBを十分に検出できない送受信のビームフォーミング、例えば図2の方式1)ではエコー情報は48dBであるので、これを従来モニタの48dBで表現でき、モニタ表現能力としては満足している。ただし、繰り返しとなるが、本来生体情報としては70dBありものが送受信ビームフォーミングの性能として48dBしか得ていないので、このような場合には従来モニタから8Kモニタに替えても効果がないことになる。方式2)の非集束型送信では、さらに検出能力が下がり、42dBとなるので、従来モニタで十分となる。
【0036】
しかし、方式3),4)では、適応型ビームファーミングを用いることにより、検出できるエコー情報のダイナミックレンジが格段と増し、この方式にも様々な方式が検討されているが、ここでは、64~70dBの情報が得られると仮定すると、従来のモニタでは表現できず不十分であり、8Kモニタではそれが可能となる。
【0037】
また、超音波診断装置の大きな特徴の一つに時間分解能がある。
【0038】
CTやMRIなどの医用画像診断装置でも近年、完像時間が格段に短くなってきているが、超音波診断装置はそれらと比較して飛躍的に速くなっている。より具体的には、図2の中の方式1)、3)の集束型送信方式では、画像のフレームレートが30-80fps程度であるが、非集束型送信方式の方式2)、4)ではそれが100~500fpsと圧倒的な速さである。このように、高速に得られる画像を従来のモニタでは60fpsの速度でしか表示できなかった。
【0039】
それに対して、8Kモニタでは120fpsで表示できることになる。これでも先ほどの100~500fpsの全てをカバーしていないと思われるが、人の視覚限界が120fpsと言われる現状では、8Kモニタで十分である。それ以上に早く得られた画像情報は実時間表示ではなく、画像分析等に役立てることが可能である。
【0040】
さらに、複雑な点は単に空間分解と時間分解能が独自にその特徴を発揮するのではなく、人の目の特性では、その両方を同時に感じ取る性質があるために、両者が同時に増すことによる相乗効果がある。
【0041】
図7には理想的な画像(画素が無限小で表現される画像)の枠が示されている。枠の中は黒い画像として表示されるとする。この理想の画像が表示系のフレーム間の時間(F0からF1、またはF1からF2までの時間)に従来モニタのピクセルの半分程度右に平行移動した状況を想定している。また、フレーム間の時間は8Kの120fpsを想定している。この場合従来モニタでは60fpsであるので、F1から次に表示されるフレームはF2となり、F1は表示されないが、ここではもし、単純に従来モニタの空間分解能は変わらず、時間分解能のみ2倍にできたとした場合を例として挙げている。
【0042】
従って、この例では、従来モニタでもF1時相で表示ができているわけであるが、理想の画像の位置が従来モニタのピクセル内での移動であるので結果的には、F0の画像と同じ表示画面となる。F2の時相になって初めて表示画像に変化が起こる。それに対して8K画像ではピクセルが小さいために、F1の時相でもF0のときとは異なる画像が表示されることになり、画像の変化を観測することが可能である。これは静止画でつまりF0の時相の画像では、空間分解能が8Kモニタで上がったとしても、従来モニタと違いが感じ取れないが、画像が時間変化する場合には空間分解能が関係してくることになり、即ち時空間で分解能を考える必要があり、その点でも8Kニターに有用性があると言える。
【0043】
ここで、超音波取得フレームレートとモニタによる画像表示フレームレートの関係をもう少し詳しく説明する。
【0044】
先に図2で説明したように、超音波の送信+受信の組み合わせで、集束型送信を行う方式(図2の1と3)では、超音波取得フレームレート(時間分解能)が30-80fpsである。それを表示するモニタの表示フレームレートが従来のモニタでは、60fpsであり、超音波フレームレートが30-60fpsまでは表示が追いついているが、それ以上の超音波画像を実時間で表示することができない状態であった。そこで、そのような場合には、画像情報を超音波診断装置内のメモリに一旦記録しておき、あとで記録画像を再生する、いわゆるスロモーション画像表示を行っていた。しかし、8Kモニタでは画像表示フレームレートが120fpsであることから、集束型送信方式での画像情報をすべて実時間表示することが可能となった。ところが、非集束型送信方式(図2の2、4)では超音波取得フレームレートが100-500fpsにもなり、もはや8Kモニタといえども実時間表示は不可能である。仮に、それを実現できたとしても、人間の目の能力外であり、人間の目が追い付いていかない。そこで、取得した画像情報はスロモーションにより我々の感知できる8Kモニタの120fpsまで落として表示することとなるが、この場合にも繰り返しになるが、従来モニタの60fpsよりも時間分解能の良い120fpsのスピードで再生が可能であり、よりスムーズなスロモーション画像再生となる。
【0045】
次に、画像補間機能について、従来モニタと8Kモニタとの性能の違いについて説明する。一般に、集束型送信方式では、先に図2に示したように超音波ビームが非集束送信と比較してビーム幅が狭くかつ送信強度が高いなどの利点がある(空間分解能の例で、図2の1と2、または3と4との比較)ものの、前述したように、超音波画像フレームレートがモニタ表示速度よりも遅い問題がある。この場合に、表示系では同じ超音波取得画像を何フレームも表示することとなり、効率が悪い。そこで、画像補間機能を使用することがある。
【0046】
図8には、その補間機能の説明図が示されている。超音波取得フレームF1とF2が取得できたと仮定する。そこで、この2枚の画像からその間の画像IF1,IF2,IF3を生成する。この生成方法には、色々な方法が考えられるが、ここでは一例として線形補間方法を示す。この様に補間機能により生成した画像を8Kモニタの高フレームレートで表示することにより、超音波取得フレームレートと同じ時間で、表示系では補間画像も含めて実時間で表示が可能となり、時間方向でも滑らかな画像表示を実現できる。
【0047】
以上は時間分解に関する効果を説明したが、8Kではそれに空間分解の効果がさらに影響することを次に説明する。
【0048】
図9において、ある点座標(V1,H1)から(V2,H2)まで移動したと仮定し、超音波取得フレームがF1とF2であったと仮定する。さらに、その間の画像IF1、IF2、IF3を、補間生成したと仮定する。これらの画像を、仮に従来の1Kモニタが高速に実時間で表示できたとしても、図9中段の1Kで示すように、補間された点の位置が1Kモニタでは同一の画素内に入るために、F1とIF1では同じ画像となり、かつ同様に、IF2,IF3とF2は同じ画像となるために、せっかく補間機能を使用し高速表示しても、結果として人間の目には、この機能を使用しない場合と同じフレームレートとしか感知されないこととなる。それに対して、同図に8Kで示された画像表示では、補間された画像は8Kモニタの画素が細かいために、F1からF2まで補間画像を含めてすべて異なった画像として表示されるために、時間分解能と空間分解能の両方を同時に感知することが可能となる。
【0049】
さらに、この補間機能の表示では、上述に加えてコントラスト分解能も関係してくることを次に説明する。
【0050】
図10に示すように、ある画素の輝度が、理想的に連続で増加している場合でも、表示モニタでは、決められた階調のどれかに対応して表示することになる。この階調数が多い程、理想の変化に近づく。その点で、従来の1Kモニタでは輝度値を256階調で表現するのに対して、8Kモニタでは16倍の4096階調である。この階調差は、次の結果をもたらす。すなわち、図8と同様に超音波取得フレームF1のある画素値が単調に増加してF2の画素値に変化したと仮定し、その間の画像IF1,IF2,IF3での値を補間機能により生成する場合を考える。この場合に、従来のモニタでは、せっかく補間機能により、異なる輝度値が生成されたとしても、その変化が1階調内である場合には、補間されなかった時と同じ表示になってしまう。しかし、16倍の輝度分解能をもつ8Kモニタでは、画素の輝度変化が補間した画像を含めてすべての画像で異なり、120fpsの速さで我々には感知できるため、あたかも連続に画像が変化しているかのように見える。
【0051】
図11では、補間機能により生成された画像の表示に関して、図9の空間分解能の違いによる表示の違いに加えて、図10の輝度情報の表示の違いを合わせた状態を図示している。このように、図10のところで8Kモニタでは、時間分解能と空間分解能の両方が関係することを述べたが、さらに輝度の変化、すなわちコントラスト分解能も同時に影響するため、1Kと8Kモニタでの補間機能の違いは、人間の目には、時間分解能と空間分解能とコントラスト分解能の3つの要素が個々に独立して感知されるのではなく、同時に感知されることとなり、大きな違いとして認識されることになる。
【0052】
以下では、本実施形態の特徴を顕著に示す信号対雑音比ビームフォーミングおよび線形回帰ビームフォーミングについて詳細に説明する。
【0053】
図12は、信号対雑音比ビームフォーミングおよび線形回帰ビームフォーミングの流れを示すフローチャートである。図2に示すように、両手法には、推定ステップSA100と、推定ステップSA100に後続する生成ステップSA110の2つのステップが含まれる。つまり、信号処理プログラムにしたがって作動している信号処理部230は、図1に示すように、推定ステップSA100を実行する推定手段230a、および生成ステップSA110を実行する生成手段230bとして機能する。
【0054】
信号対雑音比ビームフォーミングにおける推定ステップSA100では、信号処理部230は、受信部220から出力されるMチャネルのエコー信号における信号対雑音比を推定し、受信焦点からのエコーを強調する重み係数(信号対雑音比に応じた重み係数)を算出する。上述したように、受信焦点からのエコーyは、遅延補償後のエコー信号smの直流成分となる。信号対雑音比ビームフォーミングにおける推定ステップSA100では、信号処理部230は、信号成分と雑音成分を遅延補償後のエコー信号smの平均値と分散で推定し、受信焦点からのエコーを強調する重み係数WSNRを以下の数3にしたがって算出する。そして、信号対雑音比ビームフォーミングにおける生成ステップSA110では、信号処理部230は、信号対雑音比ビームフォーミングの出力(すなわち、超音波断層像を表すビームフォーマ)YSNRを以下の数4にしたがって算出して表示装置40に与える。
【数3】
【数4】
【0055】
遅延補償後のエコー信号smの信号対雑音比が非常に高い場合は、数3の分母が非常に小さくなり、その結果WSNRが極端に大きくなってビームフォーマ出力が不安定となる。それを避けるため、安定化パラメータβ(実数)を数5のように導入してもよい。
【数5】
βが0に近いほど数5の分母が小さくなることが回避され、ビームフォーマ出力は安定するが、空間分解能などの改善効果は低下する。安定化パラメータβの値については、ビームフォーマ出力の安定度と空間分解能等の改善効果との兼ね合いで適宜好適な値に定めるようにすればよい。以上が信号対雑音比ビームフォーミングの内容である。
【0056】
次いで、線形回帰ビームフォーミングについて説明する。
線形回帰ビームフォーミングにおける推定ステップSA100では、信号処理部230は、信号対雑音比ビームフォーミングにおける処理とは異なる処理で、受信部220から出力されるMチャネルのエコー信号における雑音を推定し、受信焦点からのエコーを強調する重み係数を算出する。より詳細に説明すると、信号処理部230は、まず、累積要素信号uを以下の数6にしたがって算出する(ただし、u=0)。なお、数6の右辺におけるsはi番目の超音波振動子からの遅延補償後のエコー信号である。
【数6】
【0057】
上述したように、受信焦点からのエコーyは、遅延補償後のエコー信号smの直流成分となる。したがって、累積要素信号uは、以下の数7のような線形関数モデル化される。なお、数7におけるnは付加的ノイズに起因するバイアスである。以下では、数7にしたがってモデル化した信号をモデル化要素信号と呼ぶ。測定された累積要素信号uとモデル化要素信号Uとの平均二乗差αは以下の数8のように定義され、信号処理部230は、数8で定義される平均二乗差αが最小になるように、yおよびnの値(以下、最小二乗推定値)を設定(すなわち、信号対雑音比を推定)する。なお、yおよびnの最小二乗推定値は、数9に示すように、yとnに対するαの偏微分をゼロに設定することによって得られる。
【数7】
【数8】
【数9】
【0058】
次いで、信号処理部230は、まず、数9にしたがって算出した最小二乗推定値YおよびNを、数8におけるyおよびnに代入して平均二乗差αの最小値αminを算出する。そして、信号処理部230は、受信焦点からのエコーを強調する重み係数WLRを以下の数10にしたがって算出し、線形回帰ビームフォーミングにおける推定ステップSA100を終了する。線形回帰ビームフォーミングにおける生成ステップSA110では、線形回帰ビームフォーマの出力(すなわち、超音波断層像を表すビームフォーマ出力)YLRを以下の数11にしたがって算出して表示装置40に与える。
【数10】
【数11】
【0059】
信号対雑音比ビームフォーミングと同様に、遅延補償後のエコー信号smの信号対雑音比が非常に高い場合は、数10の分母が非常に小さくなり、その結果WLRが極端に大きくなってビームフォーマ出力が不安定となる。それを避けるため、安定化パラメータγ(実数)を数12のように導入してもよい。
【数12】
γの値を大きくするほど、ビームフォーマ出力は安定するが、空間分解能などの改善効果は低減する。安定化パラメータγの値についても、前述のβと同様に、ビームフォーマ出力の安定度と空間分解能等の改善効果との兼ね合いで適宜好適な値に定めるようにすればよい。以上が線形回帰ビームフォーミングの内容である。
【0060】
上述した線形回帰ビームフォーミングの推定ステップSA100では、信号対雑音比の推定に最小二乗法が用いられているため、信号対雑音比ビームフォーミングに比べ計算負荷が高い。そこで、線形回帰ビームフォーミングの計算効率を向上させる(すなわち、計算負荷を低減させる)ため、以下のように変形してもよい。
計算効率を向上させた線形回帰ビームフォーミングの推定ステップSA100では、信号処理部230は、m番目の素子による受信信号smに含まれる雑音の積分値nmを数13により算出する。
【数13】
計算効率を向上させた線形回帰ビームフォーミングにおける重み係数WLReは、数13により得られる雑音成分の積分値nmを用いて、以下の数14のように定義される。計算効率を向上させた線形回帰ビームフォーミングの推定ステップSA100では、信号処理部230は、数14にしたがって重み係数WLReを算出する。なお、数14におけるγは数12におけるものと同様に安定化パラメータである。
【数14】
計算効率を向上させた線形回帰ビームフォーミングの生成ステップSA110では、信号処理部230は、超音波断層像を表すビームフォーマ出力YLReを以下の数15にしたがって算出し、表示装置40に与える。
【数15】
【0061】
また、信号対雑音比ビームフォーミング、線形回帰ビームフォーミング、計算効率を向上させた線形回帰ビームフォーミングの何れについても,非特許文献2に示された開口分割処理を組み合わせることで、計算量を低減させてもよい。
【0062】
図13は、超音波断層像の空間分解能を評価するための点ターゲットをイメージングした結果である。より詳細には、図13(a)はDASビームフォーミングにより得られた画像,図13(b)は受信信号間の相関性に基づく手法により得られた画像、図13(c)と(d)はそれぞれ、本実施形態の信号対雑音比ビームフォーミングおよび線形回帰ビームフォーミングにより得られた画像である。図13(a)~(d)の各図において、画像の輝度(白色の強さ)は超音波散乱波の強度を示す。図13(a)~(d)に示す各画像を比較すれば明らかなように、本実施形態の信号対雑音比ビームフォーミングおよび線形回帰ビームフォーミングにより得られた画像(図13(c)、(d))では、白の輝点のサイズが図13(a)および図13(b)に比較して小さくなっている。このことから,本実施形態の信号対雑音比ビームフォーミングおよび線形回帰ビームフォーミングによれば、DASビームフォーミングおよび受信信号管の相関性に基づく手法に比較して、高い空間分解能が得られることが判る。
【0063】
図14は、超音波断層像のコントラストを評価するためのファントム(虚像)を画像化したものである。より詳細に説明すると、図14(a)は、従来のDASビームフォーミングにより得られた画像、図14(b)は、受信信号間の相関性により得られた画像、図14(c)と図14(d)はそれぞれ、本実施形態の信号対雑音比ビームフォーミングおよび線形回帰ビームフォーミングにより得られた画像である。なお、図14(a)~図14(d)の各図においては、中央部の暗い部分は超音波散乱波が発生しない媒質(具体的には、嚢胞模擬部)であり,黒一色で描出されることが望ましい。図14(a)および図14(b)の各画像では,嚢胞模擬部においても白い輝点が発生しており,これらの白い輝点は虚像である。図14(c)では虚像が低減されていることがわかる。さらに、図14(d)に示す画像では嚢胞模擬部において虚像が描出されていない。つまり、本実施形態の信号対雑音比ビームフォーミングおよび線形回帰ビームフォーミングによれば、従来のDASビームフォーミングおよび受信信号管の相関性に基づく手法に比較して、嚢胞模擬部における虚像の描出を抑制することができ、コントラストが向上していることが判る。
【0064】
図14において、虚像の抑圧効果が信号対雑音比ビームフォーミングに比べ線形回帰ビームフォーミングの方が高かった理由は、数6の処理による効果、すなわち、遅延補償後のエコー信号sの積分効果である。積分は低域通過フィルタに対応する。積分操作以外のフィルタを適用することで、線形回帰ビームフォーマの出力をさらに向上させることも可能である。なお、数13の積分処理についても同様に他のフィルタ処理に適宜変更してもよい。
【0065】
以上説明したように、本発明によれば、超音波断層像の空間分解能およびコントラストを、従来のDASビームフォーミングは勿論、超音波受信信号間の相関性に基づく手法と比較しても、さらに向上させることが可能になる。
【0066】
(B.変形)
以上、本発明の一実施形態について説明したが、この実施形態に以下の変形を加えても勿論よい。
(1)上記実施形態では、超音波医用システムへの本発明の適用例を説明したが、物体の非破壊検査など、医療用以外の超音波断層像の生成への適用も可能である。医療用以外の技術分野においても、超音波断層像の空間分解能およびコントラストは,高ければ高い程好ましいからである。
【0067】
(2)線形回帰ビームフォーミングの推定ステップSA100では、Mチャネルのエコー信号の信号対雑音比を最小二乗法により推定したが、尤度を用いる手法など他の手法により信号対雑音比を推定してもよい。要は、M(2以上の自然数)個の超音波振動子から発せられた超音波のエコーを受信してエコー信号を出力する超音波プローブから出力されるMチャネルのエコー信号における雑音を推定し、受信焦点からのエコーを強調する重み係数を前記Mチャネルのエコー信号における信号対雑音比に応じて算出する推定ステップと、超音波断層像を表すビームフォーマを、前記推定ステップにて算出された重み係数を用いて前記Mチャネルのエコー信号から生成する生成ステップと、を含むことを特徴とする超音波断層像生成方法であればよい。
【0068】
(3)上記実施形態では、超音波断層像生成装置20の信号処理部230が推定手段230aおよび生成手段230bとして機能したが、制御部200を推定手段230aおよび生成手段230bとして機能させてもよい。具体的には、受信部220の出力信号を制御部200に与え、制御部200に上記実施形態の信号処理プログラムを実行させるようにすればよい。
【0069】
(4)上記実施形態では、本実施形態の特徴を顕著に示す超音波断層像生成方法を実現する信号処理プログラムが超音波断層像生成装置20に予めインストールされていた。しかし、CPUなどのコンピュータを、M(2以上の自然数)個の超音波振動子を有し、各超音波振動子から発せられた超音波のエコーを受信してエコー信号を出力する超音波プローブから出力されるMチャネルのエコー信号における雑音を推定し、受信焦点からのエコーを強調する重み係数を前記Mチャネルのエコー信号における信号対雑音比に応じて算出する推定手段と、超音波断層像を表すビームフォーマを、前記推定手段により算出された重み係数を用いて前記Mチャネルのエコー信号から生成する生成手段と、して機能させるプログラムを単体で製造し、販売等の配布をしてもよい。上記プログラムの配布態様の具体例としては、インターネットなどの電気通信回線経由のダウンロードにより配布する態様や、CD-ROM(Compact Disk-Read Only Memory)やフラッシュROM(Read Only Memory)などのコンピュータ読み取り可能な記録媒体に書き込んで配布する態様が挙げられる。このようにして配布される上記プログラムにしたがって、コンピュータを作動させることで、当該コンピュータに、本発明の超音波断層像生成方法を実行させることが可能になるからである。
【0070】
(5)上記実施形態では、本実施形態の特徴を顕著に示す超音波断層像生成方法の各ステップを実行する推定手段230aおよび生成手段230bがソフトウェアモジュールとして実現されていた。しかし、M(2以上の自然数)個の超音波振動子を有し、各超音波振動子から発せられた超音波のエコーを受信してエコー信号を出力する超音波プローブから出力されるMチャネルのエコー信号における雑音を推定し、受信焦点からのエコーを強調する重み係数を前記Mチャネルのエコー信号における信号対雑音比に応じて算出する推定手段と、超音波断層像を表すビームフォーマを、前記推定手段により算出された重み係数を用いて前記Mチャネルのエコー信号から生成する生成手段の各々をASICなどの電子回路で構成し、これら電子回路を組み合わせて本発明の超音波断層像生成を構成してもよい。
【符号の説明】
【0071】
1…超音波医用システム、10…超音波プローブ、20…超音波断層像生成装置、30…操作装置、40…表示装置、200…制御部、210…送信部、220…受信部、230…信号処理部、230a…推定手段、230b…生成手段。
図1
図2
図3
図4
図5
図6
図7
図8
図9
図10
図11
図12
図13
図14