JP2015013099A - 超音波測定装置および超音波測定方法 - Google Patents
超音波測定装置および超音波測定方法 Download PDFInfo
- Publication number
- JP2015013099A JP2015013099A JP2014038976A JP2014038976A JP2015013099A JP 2015013099 A JP2015013099 A JP 2015013099A JP 2014038976 A JP2014038976 A JP 2014038976A JP 2014038976 A JP2014038976 A JP 2014038976A JP 2015013099 A JP2015013099 A JP 2015013099A
- Authority
- JP
- Japan
- Prior art keywords
- blood vessel
- ultrasonic
- determination unit
- signal
- temporal change
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
【課題】より少ない処理負荷でプローブ走査範囲内から血管位置を判定する。
【解決手段】超音波測定では、血管は拍動に伴って常に動くので血管の前壁と後壁を含む走査線上の反射信号の信号強度変化は、血管以外の生体組織のそれにくらべて大きくなる。よって、超音波振動子別に、隣接する2フレーム分の反射波信号が得られる都度、深さ別に信号強度差を算出し、全深さの信号強度差の合計を超音波振動子別に所定時間分積算する。信号強度差の積算値を現わすヒストグラムのピーク値を検索し、ピーク値に対応する超音波振動子を「血管直上にある超音波振動子」として判定する。
【選択図】図4
【解決手段】超音波測定では、血管は拍動に伴って常に動くので血管の前壁と後壁を含む走査線上の反射信号の信号強度変化は、血管以外の生体組織のそれにくらべて大きくなる。よって、超音波振動子別に、隣接する2フレーム分の反射波信号が得られる都度、深さ別に信号強度差を算出し、全深さの信号強度差の合計を超音波振動子別に所定時間分積算する。信号強度差の積算値を現わすヒストグラムのピーク値を検索し、ピーク値に対応する超音波振動子を「血管直上にある超音波振動子」として判定する。
【選択図】図4
Description
本発明は、超音波測定により生体内組織を測定する超音波測定装置等に関する。
超音波測定装置により体内の生体情報を測定することで、血管機能の評価や血管疾患の判断が行われている。
例えば、動脈硬化の指標となる頸動脈のIMT(Intima Media Thickness:内膜中膜複合体厚)を計測することもその1つである。頸動脈の測定では頸動脈を見つけ、測定ポイントを適切に決定しなければならない。通常は、オペレーターが医学的知識に基づいて目的とする頸動脈のおおよその位置に超音波プローブを当て、モニターに表示されるBモード画像を見ながら測定対象とする頸動脈を詳細に探し出し、探し出した頸動脈を計測ポイントとして手動で設定する。
例えば、動脈硬化の指標となる頸動脈のIMT(Intima Media Thickness:内膜中膜複合体厚)を計測することもその1つである。頸動脈の測定では頸動脈を見つけ、測定ポイントを適切に決定しなければならない。通常は、オペレーターが医学的知識に基づいて目的とする頸動脈のおおよその位置に超音波プローブを当て、モニターに表示されるBモード画像を見ながら測定対象とする頸動脈を詳細に探し出し、探し出した頸動脈を計測ポイントとして手動で設定する。
こうした一連の測定操作を速やかに実行し且つ適切に頸動脈を見つけるには熟練を要するとされてきたが、近年、測定操作を補助する機能が考案されるようになった。例えば、特許文献1には、超音波ビームの反射波の受信信号強度を利用して自動で血管を検出する方法が開示されている。
特許文献1に開示されている検出方法では、血管壁からの反射波が最も強いことを前提として血管位置を特定している。しかし、実際には周辺組織など血管以外の組織からも強い反射波を受信する場合がある。また、フォーカスの設定値などによって各深さでの受信強度が変化する場合もある。そのため、一概に反射波の強度のみで血管を特定することは難しい。
また、一般的な画像認識の原理を用いて血管を検出する方法として、テンプレートマッチングやオプティカルフローなどなどが挙げられるが、何れもアルゴリズムが複雑な上、処理負荷も大きいため、小型な装置でしかも短時間に血管を検出しようとする場合には向かない。
具体的には、テンプレートマッチングを用いる検出方法では、予め血管断面のBモード画像などのテンプレート画像を用意しておいて、プローブ走査範囲の画像において最もマッチングする位置を探索する。しかし、フォーカスやプローブの角度などの測定条件によって超音波画像の見え方が異なる場合が多く、血管断面がきれいな円形で映し出されるとは限らない。さらには、それらの見え方や血管の大きさなどの個人差を考慮すると、テンプレートマッチングを用いて血管を自動検出する処理は、えてして複雑で演算負荷が大きくなってしまう傾向にある。
オプティカルフローは、時間的に連続した複数の画像間の相対的な動きによって生じる物体の動きベクトルを求めて、画像の中から動きのある物体を検出する方法である。拍動による血管壁の移動量が、プローブ走査範囲中で最も大きいと考えれば、オプティカルフローを用いて血管の位置を検出することができる。しかし、動きベクトルを求めるにはブロックマッチングや輝度勾配の拘束方程式を用いた勾配法などを用いる必要があり、やはり複雑で演算負荷が大きくなってしまう傾向にある。
複雑で演算負荷が大きいと、必然的に高速なCPUや専用の信号処理回路を必要とするためどうしても装置が大型化しやすく、高価な装置になる。
本発明は、こうした事情を鑑みて、プローブ走査範囲内から血管位置を判定することをより簡単に実現する手法の実現を目的として考案されたものである。
以上の課題を解決するための第1の発明は、生体内組織に対して超音波を照射する制御を行う照射制御部と、前記超音波の反射波信号を受信する受信部と、前記反射波信号に含まれる信号の時間的変化に基づいて、前記生体内組織中の血管の位置を判定する血管位置判定部と、を備えた超音波測定装置である。
別形態として、生体内組織に対して超音波を照射することと、前記超音波の反射波信号を受信することと、前記反射波信号に含まれる信号の時間的変化に基づいて、前記生体内組織中の血管の位置を判定することと、を含む超音波測定方法を構成することができる。
第1の発明及びその別形態によれば、照射範囲の反射波信号の時間的変化を求めることができる。血管とそれ以外の生体組織との動きに着目すれば、血管は拍動に伴って大きく動くが、生体組織の動きは血管に比べて小さい。
よって、照射範囲において、大きな時間的変化を検出している部分に血管が存在すると推測することできる。逆に、小さい時間的変化しか検出できない部分は、血管以外であると推測できる。つまり、血管位置を判定することができる。しかも、判定に係る主たる処理は反射波信号の時間的変化の算出であり、従来の方法に比べて簡単に実現することができる。
よって、照射範囲において、大きな時間的変化を検出している部分に血管が存在すると推測することできる。逆に、小さい時間的変化しか検出できない部分は、血管以外であると推測できる。つまり、血管位置を判定することができる。しかも、判定に係る主たる処理は反射波信号の時間的変化の算出であり、従来の方法に比べて簡単に実現することができる。
第2の発明は、前記血管位置判定部が、前記時間的変化の積算結果が所定の高変化条件を満たす位置に基づいて、前記血管の位置を判定する、第1の発明の超音波測定装置である。
血管とそれ以外の生体組織の動きに着目すれば、血管は拍動に伴い必ず動くので時間的変化の積算結果は血管以外の生体組織のそれに比べてはるかに多くなる。よって、第2の発明によれば、血管位置の判定確度を高めることができる。
また、第3の発明は、更に前記血管位置判定部が、前記時間的変化の発生が所定の周期条件を満たす位置に基づいて、前記血管の位置を判定する、第1又は第2の発明の超音波測定装置である。
血管の動きに伴う時間的変化は拍動と高い相関のある周期性を有する。よって、第3の発明によれば、血管位置を更に的確に識別することが可能であり判定確度をより一層高めることができる。
第4の発明は、前記血管位置判定部が、生体表面上の同一位置における深さ方向の信号の時間的変化を当該同一位置における時間的変化として統計処理することで、前記血管上の生体表面の位置を判定する第1位置判定部を有する、第1〜第3の何れかの発明の超音波測定装置である。
更に第5の発明は、前記第1位置判定部が、前記血管上の生体表面の位置において、前記血管の前壁部の位置の時間的変化と後壁部の位置の時間的変化とを当該表面位置における時間的変化として合算して統計処理する、第4の発明の超音波測定装置である。
第4及び第5の発明によれば、血管上の生体表面の位置を判定することができる。特に、生体表面上の同一位置における深さ方向の信号の時間的変化を当該同一位置における時間的変化として統計処理するため、血管上の生体表面の位置をより適切かつ高速に求めることができる。
第6の発明は、前記第1位置判定部は、生体表面上の同一位置における深さ方向の信号を、前記受信部での受信タイミングが異なる当該信号の組み合わせで自己相関演算した相関値を、当該同一位置における時間的変化を示す値とする、第4の発明の超音波測定装置である。
この第6の発明によれば、生体表面上のある位置における深さ方向の信号を受信タイミングが異なる信号の組み合わせで自己相関演算することで、当該位置における時間的変化をします値を得ることができる。
第7の発明は、前記血管位置判定部が、前記第1位置判定部で判定された判定位置における深さ方向の各位置での信号の時間的変化に基づいて、前記血管の深さ位置を判定する第2の位置判定部を有する、第4〜第6の何れかの発明の超音波測定装置である。
第7の発明によれば、血管直上にあると判定された生体表面位置における血管の深さ位置を更に判定することができる。
なお、適切な判定精度を得るために、第8の発明として、前記血管位置判定部が、20回/秒以上のサンプリング速度で前記時間的変化を判定する、第1〜第7の何れかの発明の超音波測定装置を構成することができる。
図1は、本実施形態における超音波測定装置10のシステム構成例を示す図である。超音波測定装置10は、超音波の反射波を測定することにより被検体2の生体情報を測定する装置である。本実施形態では、血管4の位置を自動的に判定し、生体情報の1つとしてIMT(Intima Media Thickness:血管の内膜中膜複合体厚)を算出する。
超音波測定装置10は、測定結果や操作情報を画像表示するための手段及び操作入力のための手段を兼ねるタッチパネル12と、操作入力をするためのキーボード14と、超音波プローブ16(深触子)と、処理装置30とを備える。処理装置30には、制御基板31が搭載されており、タッチパネル12,キーボード14,超音波プローブ16などの装置各部と信号送受可能に接続されている。
制御基板31には、CPU32や、ASIC(application specific integrated circuit)、各種LSI(Large Scale Integration)の他、ICメモリーやハードディスク等による記憶媒体33と、外部装置とのデータ通信を実現する通信IC34とが搭載されている。処理装置30は、記憶媒体33に記憶されている測定プログラムをCPU32等が実行することにより本実施形態に係る各種機能を実現する。
具体的には、処理装置30の制御により、超音波測定装置10は超音波プローブ16から被検体2へ超音波パルスを発信・照射し、その反射波を受信する。そして受信した反射波を増幅・信号処理することにより、被検体2の生体内構造の位置情報や経時変化などを判定することができる。反射波信号には、いわゆるAモード、Bモード、Mモード、カラードップラーの各モードの画像が含まれる。勿論、これら以外の形式のデータでもよい。超音波を用いた測定(サンプリング)は所定周期で繰り返し実行される。測定単位を「フレーム」と呼称する。本実施形態のサンプリングは20fps(Frames Per Second)以上で行われる。
また、超音波測定装置10は、基準となる反射波信号に関心領域(トラッキングポイント)を設定することで、各関心領域を異なるフレーム間で追跡し変位を算出する、いわゆる「トラッキング処理」を行うことができる。
そして、超音波測定装置10は、超音波測定により得られる反射波信号に基づいてプローブ走査範囲における血管位置の判定と、判定された位置の血管について超音波測定を実行してIMTを算出することができる。
[原理の説明]
さて、生体組織の動きに着目すると、血管は心臓の拍動に伴い周期的に大きく動くが、血管の周りにある他の生体組織の動きは血管に比べてずっと小さい。本実施形態では、この知見に基づいて血管位置の判定を行う。
さて、生体組織の動きに着目すると、血管は心臓の拍動に伴い周期的に大きく動くが、血管の周りにある他の生体組織の動きは血管に比べてずっと小さい。本実施形態では、この知見に基づいて血管位置の判定を行う。
先ず、第1の位置判定として血管直上となる生体表面の位置を判定する。
図2は、超音波プローブ16を被検体2の体表面に当てて超音波測定している状態を簡略的に示す図であり、血管4の短軸方向の断面で示した図である。超音波プローブ16には、複数の超音波振動子18が内蔵されている。図2の例では、一つ一つの超音波振動子18からは1本の超音波ビームが図の上から下に向けて照射される。超音波振動子18が覆う範囲がプローブ走査範囲Asとなる。なお、超音波振動子18は、図2に向かって奥行方向にも複数列配設された、すなわち、平面状に配設された構成であってもよいし、図2に向かって奥行方向には1列のみで左右方向にのみ配列された列状の構成であってもよい。
図2は、超音波プローブ16を被検体2の体表面に当てて超音波測定している状態を簡略的に示す図であり、血管4の短軸方向の断面で示した図である。超音波プローブ16には、複数の超音波振動子18が内蔵されている。図2の例では、一つ一つの超音波振動子18からは1本の超音波ビームが図の上から下に向けて照射される。超音波振動子18が覆う範囲がプローブ走査範囲Asとなる。なお、超音波振動子18は、図2に向かって奥行方向にも複数列配設された、すなわち、平面状に配設された構成であってもよいし、図2に向かって奥行方向には1列のみで左右方向にのみ配列された列状の構成であってもよい。
図2(2)に示すように、血管は心臓の拍動により等方的に収縮拡張を繰り返す。しかし、超音波ビームのビーム方向に直交する面ほど強い反射波を受信できるが、ビーム方向と平行に近いほど反射波を受信し難い。そのため、超音波測定では、図2(1)に示すように、血管4の前壁4f及び後壁4rからの反射波は強く検出されるが、横壁4sからの反射波は弱くなる。言い換えれば、プローブ走査範囲Asに血管4が有れば、その直上にある超音波振動子18の位置における反射波信号には、前壁と後壁に係る強い反射波が現れることになる。
なお、説明の簡明化および理解の容易化のため、血管4の直上の位置を判定することを、血管4の直上の位置の超音波振動子18を判定することとして説明する。
また、反射波信号の受信は、それぞれの超音波振動子18が個別で行うが、超音波振動子18それぞれの受信信号から信号位相差の算出や受信ビームの形成処理等を行うため、複数の超音波振動子18それぞれが反射波信号を受信する必要がある。しかし、説明の簡明化および理解の容易化のため、強い反射波を含む反射波信号の位置を判定することを、強い反射波を含む反射波信号を受信した超音波振動子18を判定することとして説明する。すなわち、「超音波振動子の位置」の判定を、適宜「超音波振動子」の判定として説明する。
また、反射波信号の受信は、それぞれの超音波振動子18が個別で行うが、超音波振動子18それぞれの受信信号から信号位相差の算出や受信ビームの形成処理等を行うため、複数の超音波振動子18それぞれが反射波信号を受信する必要がある。しかし、説明の簡明化および理解の容易化のため、強い反射波を含む反射波信号の位置を判定することを、強い反射波を含む反射波信号を受信した超音波振動子18を判定することとして説明する。すなわち、「超音波振動子の位置」の判定を、適宜「超音波振動子」の判定として説明する。
図3は、血管直上にある超音波振動子18の位置における反射波信号の例を示す図である。図3(1)は第1フレームにおける測定結果を示す「深さ−信号強度グラフ」であり、図3(2)は第2フレームにおける測定結果を示す「深さ−信号強度グラフ」である。図3(3)は第1フレームと第2フレームとの間の「深さ−信号強度グラフ」の差分を示す「フレーム間信号強度差グラフ」である。
前述のように、血管4があればその前壁と後壁に係る強い反射波が検出される。図3(1)及び図3(2)においても、体表面近くの反射波群よりも深い位置に明確に識別できる2つの強い反射波のピークが現れている。そして、第1フレームと第2フレームとの間で深さ毎の信号強度差を求めると、図3(3)のグラフが得られ、フレーム間で血管の前壁及び後壁の動きが明確となる。
図3(3)のグラフから明らかなように、血管以外の生体組織もわずかながらも拍動等の影響で動いているので多少の信号強度差が生じるが、血管(より具体的には血管の前壁や後壁)ほど大きな値が検出されることはない。ましてや、血管直上にない超音波振動子18における反射波信号の信号強度差グラフではこうしたピークは見られない。つまり、拍動に伴う血管の動きは2フレーム間における信号強度の変化に現れるといえる。
本実施形態では、血管の動きに相応しい信号強度の変化が測定されたからといって直ちにその超音波振動子18が血管直上にあるとは見なさず、信号強度の変化を統計処理するとともに、信号強度変化の周期性に基づいて最終的な判定を行う。
図4は、連続する2フレーム間における信号強度の変化の統計処理について説明するための概念図である。図4(1)は、各超音波振動子18における反射波の信号強度を輝度に変換した画像すなわちBモード画像である。図4(2)は、連続する2フレーム間での各超音波振動子における信号強度変化を積算したヒストグラムである。ここで注意する点は、図3(3)の横軸が深さ方向であるのに対して、図4(2)の横軸は超音波振動子の配列順であることである。すなわち、図3(3)は、1つの超音波振動子が受信した結果に基づくグラフである。例えば、この超音波振動子が「Tr1」であるとする。図3(3)の各深さでの値(信号強度差)を積算した値が、図4(2)の「Tr1」における値(積算信号強度変化)となる。
具体的に説明する。ヒストグラムは、連続する2フレーム分の超音波測定が行われる都度、超音波振動子毎に、全ての深さにおける信号強度差の合計を算出することを繰り返し、当該信号強度差の合計を所定時間分(例えば、少なくとも心周期1〜数拍分。数秒程度。)積算することで得られる。
換言すれば、生体表面上の同一位置における深さ方向での信号の時間的変化を当該同一位置における時間的変化として積算する統計処理をした結果である。より具体的には、血管上の生体表面の位置において、血管の前壁部の位置の時間的変化と後壁部の位置の時間的変化とを当該表面位置における時間的変化として合算して統計処理する。
換言すれば、生体表面上の同一位置における深さ方向での信号の時間的変化を当該同一位置における時間的変化として積算する統計処理をした結果である。より具体的には、血管上の生体表面の位置において、血管の前壁部の位置の時間的変化と後壁部の位置の時間的変化とを当該表面位置における時間的変化として合算して統計処理する。
連続する2フレーム分の超音波測定から得られる信号強度差の合計は、血管上にない超音波振動子の合計よりも血管上にある超音波振動子の合計の方が、大きな値を示す。しかも、血管4の中心の直上にある超音波振動子18ほど大きな値を示す。当然、信号強度差の合計を積算したヒストグラムの縦軸の値、すなわち積算信号強度差にもそれは現れる。
よって、ヒストグラムの縦軸の値が所定の高変化条件を満たした超音波振動子18を「血管直上にある超音波振動子」と判定することができる。より具体的には、ヒストグラムの縦軸の値のピークに対応する超音波振動子18を「血管直上にある超音波振動子」と判定する。図4の例では、超音波振動子「Tr1」が該当する。
よって、ヒストグラムの縦軸の値が所定の高変化条件を満たした超音波振動子18を「血管直上にある超音波振動子」と判定することができる。より具体的には、ヒストグラムの縦軸の値のピークに対応する超音波振動子18を「血管直上にある超音波振動子」と判定する。図4の例では、超音波振動子「Tr1」が該当する。
本実施形態では、こうした統計処理を経て判定された超音波振動子について更に信号強度変化の周期性の観点から「血管直上にある超音波振動子」としての妥当性を判断する。
図5は、信号強度の変化の周期性に基づく血管位置の判定の原理について説明するための概念図である。図5(1)は心周期2拍分の血管径の時間変化を示すグラフである。図5(2)は、当該グラフにおける各点間(t1〜t2,t2〜t3,…)におけるフレーム間信号強度差グラフの概略図であって、前壁及び後壁に該当する部分の信号強度差ピークを分かり易く示した図である。なお、図5(2)の破線で描画された2重円は、内円が心拡張期の血管、外円が心収縮期の血管を示している。
統計処理により判定された超音波振動子(図4の例ではTr1)について、判定の過程で既に得られている複数のフレーム間信号強度差グラフ(図3(3))から、フレーム間の信号強度変化の周期を推定する。
先に述べたように、統計処理のために少なくとも心周期1〜数拍分、数秒程度の超音波測定が連続的に行われる。そして、超音波測定のフレームレートは20フレーム/秒であり、既に得られている複数のフレーム間信号強度差グラフに現れている前壁及び後壁部分の信号強度差は、血管が拍動に伴い拡張/縮小する血管壁部の移動速度に相当する。
そこで、既に得られている複数のフレーム間信号強度差グラフ毎に、各深さの信号強度差の合計を算出し、合計のピークが得られるフレーム間信号強度差グラフを検索する。図5の例では、「t2〜t3」と「t6〜t7」におけるフレーム間信号強度差グラフで、ピークが得られる。そして、「t2」と「t6」間の時間差から、ピーク間の周期(Cpeak)を求める。このピーク間周期(Cpeak)が、心周期と近似していれば拍動に伴う血管の動きを捉えたものであると確証を持つことができる。
よって、本実施形態では、図4で示した統計処理により選択された超音波振動子のうち、ピーク間周期(Cpeak)が心周期に近似しているものを最終的に「血管直上にある超音波振動子」と判定する。これで、血管上の生体表面の位置を判定する第1の位置判定が完了したことになる。
血管上の生体表面の位置を判定したならば、血管の深さを判定する第2の位置判定を行う。
図6は、血管の深さ位置を判定する方法の概念図である。図4で示した統計処理により判定された超音波振動子(Tr1)についての超音波測定結果(図3(1)(2)に相当)から、血管の前壁及び後壁に相当する信号強度のピークを検索することで血管軸の深さ(Hc)を特定する。
図6は、血管の深さ位置を判定する方法の概念図である。図4で示した統計処理により判定された超音波振動子(Tr1)についての超音波測定結果(図3(1)(2)に相当)から、血管の前壁及び後壁に相当する信号強度のピークを検索することで血管軸の深さ(Hc)を特定する。
前述の第1の位置判定と、第2の位置判定とを行うことで、プローブ走査範囲Asにおける血管の位置(2次元位置)が特定されたことになる。
[機能構成の説明]
次に、本実施形態を実現するための機能構成について説明する。
図7は、本実施形態の超音波測定装置10の機能構成例を示すブロック図である。超音波測定装置10は、操作入力部100と、送信部102と、受信部104と、処理部200と、画像表示部360と、記憶部500とを備える。
次に、本実施形態を実現するための機能構成について説明する。
図7は、本実施形態の超音波測定装置10の機能構成例を示すブロック図である。超音波測定装置10は、操作入力部100と、送信部102と、受信部104と、処理部200と、画像表示部360と、記憶部500とを備える。
操作入力部100は、オペレーターによる各種操作入力を受け付け、操作入力に応じた操作入力信号を処理部200へ出力する。ボタンスイッチやレバースイッチ、ダイヤルスイッチ、トラックパッド、マウス、などにより実現できる。図1のタッチパネル12やキーボード14がこれに該当する。
超音波プローブ16は、公知の超音波プローブを用いて構成され、送信部102は、パルス電圧に基づいて超音波を照射する。また、受信部104は、送信部102によって照射された超音波が被検体2の生体内で反射した反射波信号を受信し、電気信号に変換して出力する。
処理部200は、例えば、CPUやGPU等のマイクロプロセッサーや、ASIC、ICメモリーなどの電子部品によって実現される。そして、処理部200は、各機能部との間でデータの入出力制御を行い、所定のプログラムや各種データに基づいて各種の演算処理を実行して、血管位置を判定し、被検体2の生体情報を算出する。図1の処理装置30及び制御基板31がこれに該当する。
本実施形態では、処理部200は、超音波測定制御部210と、心拍測定部220と、血管位置判定部230と、生体情報算出部240と、画像生成部260とを有する。
超音波測定制御部210は、照射制御部212と、送受信制御部214と、受信合成部216とを有し、超音波測定を統合的に制御する。本実施形態でのサンプリングレートは20回/秒以上のサンプリング速度とし、その一例として20fpsで測定するものとする。
照射制御部212は、超音波プローブ16から発信する超音波パルスのタイミングを制御し、送信制御信号を送受信制御部214へ出力する。
送受信制御部214は、照射制御部212からの送信制御信号に従ってパルス電圧を発生させて送信部102へ出力する。その際、送信遅延処理を行って各超音波振動子へのパルス電圧の出力タイミングの調整を行うことができる。また、送受信制御部214は、受信部104から出力された反射波信号の増幅やフィルター処理を行って、その結果を受信合成部216へ出力することができる。
受信合成部216は、必要に応じて遅延処理等を行っていわゆる受信信号のフォーカスに係る処理等を実行して反射波信号を生成する。
心拍測定部220は、心周期拍動を測定する。本実施形態では、反射波信号のピークを検索することで心周期拍動を測定するものとするが、外部の血圧測定装置に接続して外部の血圧計測装置から取得した血圧の変動から心周期拍動を測定するとしても良い。
血管位置判定部230は、反射波信号に含まれる照射範囲の各部の信号の時間的変化に基づいて生体内組織中の血管の位置を判定する。具体的には、反射波信号の信号強度の時間的変化の積算結果に基づいて血管の位置を判定して、判定結果を報知する。より具体的には、周期基準判定部232と、第1位置判定部234と、第2位置判定部236と、報知制御部238とを有する。
周期基準判定部232は、反射波信号の信号強度の時間的変化の発生周期を推定し、当該周期が所定の周期条件(例えば、心拍測定部220で測定される心周期に近似することを示す条件)を満たす位置に基づいて血管の位置を判定する(図5参照)。
第1位置判定部234は、血管の直上の生体表面の位置を判定する。本実施形態では、血管上の生体表面の位置において、血管の前壁部の位置の時間的変化と後壁部の位置の時間的変化とを当該表面位置における時間的変化として合算して統計処理する(図3〜図4参照)。
より具体的には、第1位置判定部234は、(1)超音波振動子のフレーム間信号強度差を全深さについて求める処理、(2)全深さのフレーム間信号強度差の合計を超音波振動子別に所定時間積算したヒストグラムを作成する処理、(3)ヒストグラムのピーク値を検索し、検索されたピーク値に対応する超音波振動子を選出する処理、を実行する。
第2位置判定部236は、第1位置判定部234で判定された判定位置における深さ方向の各位置での信号の時間的変化に基づいて血管の深さ位置を判定する(図6参照)。
報知制御部238は、第1位置判定部234及び第2位置判定部236による判定結果を報知する制御を行う。具体的には、第1位置判定部234により選出された超音波振動子に対応する位置を、画像表示部360に表示させる制御を行う。第1位置判定部234により選出された超音波振動子に対応する位置は生体表面(皮膚面)の位置であり、血管の直上位置である。従って、当該位置を表示することで、血管直上の位置を報知することができる。また、第2位置判定部236により判定された深さ位置を画像表示部360に表示させる制御を行う。これにより、血管の深さ位置を報知することができる。なお、「報知」は、表示に限らず、音声出力や、所定の通知信号の出力等でもよい。
生体情報算出部240は、血管位置判定部230にて血管の直上にあると判定された超音波振動子(受信部104)の反射波信号に基づいて目的とする生体情報を算出する。本実施形態ではIMTを算出する。生体情報の算出方法は、公知技術を適宜利用可能である。例えば、内腔内膜境界と中外膜境界を判定して、これらの境界間の距離からIMTを算出する。なお、算出する生体情報の種類は、適宜追加・変更可能である。例えば、血管径を算出し、スティフネスパラメーターβを用いた血圧算出式から血圧を算出するとしてもよい。その場合、図1で示した構成に適宜他の測定装置(例えば、血圧計など)を追加することができる。
画像生成部260は、超音波測定や、生体情報測定に必要な各種操作画面や測定結果を表示するための画像を生成し画像表示部360へ出力する。
画像表示部360は、画像生成部260から入力される画像データを表示する。図1のタッチパネル12がこれに該当する。
記憶部500は、ICメモリーやハードディスク、光学ディスクなどの記憶媒体により実現され、各種プログラムや、処理部200の演算過程のデータなどの各種データを記憶する。図1では、処理装置30の制御基板31に搭載されている記憶媒体33がこれに該当する。なお、処理部200と記憶部500の接続は、装置内の内部バス回路による接続に限らず、LAN(Local Area Network)やインターネットなどの通信回線で実現しても良い。その場合、記憶部500は超音波測定装置10とは別の外部記憶装置により実現されるとしてもよい。
記憶部500は、測定プログラム501と、反射波信号510と、フレーム間信号強度差データ512と、統計処理データ514と、ピーク周期データ518と、血管位置データ520と、生体情報測定結果530とを記憶する。
処理部200は、測定プログラム501を読み出して実行することにより、超音波測定制御部210や、心拍測定部220、血管位置判定部230、生体情報算出部240、画像生成部260等の機能を実現する。なお、これらの機能部を電子回路等のハードウェアで実現する場合には、当該機能を実現させるためのプログラムの一部を省略することができる。例えば、血管位置判定部230をLSI等で実現するならば、血管位置判定部230の機能を実現させるための血管位置判定プログラム502を省略できる。
反射波信号510は、超音波測定により得られた反射波信号のデータであって、超音波測定制御部210によってフレーム毎に生成される。例えば、1つの反射波信号510には、超音波振動子の識別情報(Tr)と測定されたフレーム識別情報(fr)とが対応づけて格納される。図3(1)、図3(2)のグラフのデータに相当する。
フレーム間信号強度差データ512は、同じ超音波振動子の異なるフレーム間における深さ別の信号強度差のデータであって、例えば、連続する2フレーム分の超音波測定が完了する毎に各超音波振動子について生成される。図3(3)のグラフのデータに相当する。1つのフレーム間信号強度差データ512は、対応する超音波振動子の識別情報(Tr)と、比較に使用されたフレームを識別するための2つのフレーム識別情報(fr1,fr2)とが対応づけて格納されている。
統計処理データ514は、所定期間分の反射波信号の信号強度変化の統計処理データであって、所定期間分の超音波測定が完了すると作成される。図4(2)のヒストグラムのデータに相当する。
ピーク周期データ518は、超音波振動子の識別情報(Tr)と対応づけて、フレーム間信号強度の推定周期(Cpeak)を格納する。
血管位置データ520は、血管上の生体表面の位置情報に相当する超音波振動子の識別情報(Tr)と対応づけて、血管の深さ位置(Hc)を格納する。
生体情報測定結果530は、算出された生体情報を時系列に格納する。
生体情報測定結果530は、算出された生体情報を時系列に格納する。
なお、記憶部500は、これら以外にも、各種フラグ、計時用のカウンター値など、血管位置の判定並びに生体情報の算出に必要なデータを適宜記憶することができる。
[処理の流れの説明]
次に、超音波測定装置10の動作について説明する。
図8は、超音波測定装置10の血管位置判定並びに生体情報算出に係る処理の流れを説明するためのフローチャートである。
処理部200は、先ず所定のサンプリング周期での超音波測定を開始する(ステップS2)。サンプリング周期、すなわちフレームレートは20fps以上とすると好適である。以降、記憶部500に逐次、反射波信号510が記憶される。
次に、超音波測定装置10の動作について説明する。
図8は、超音波測定装置10の血管位置判定並びに生体情報算出に係る処理の流れを説明するためのフローチャートである。
処理部200は、先ず所定のサンプリング周期での超音波測定を開始する(ステップS2)。サンプリング周期、すなわちフレームレートは20fps以上とすると好適である。以降、記憶部500に逐次、反射波信号510が記憶される。
また、処理部200は超音波測定の開始に伴って心周期拍動の測定とそれに基づく心周期の判定を開始する(ステップS4)。
次に、連続する複数フレームの超音波測定の測定ができたならば(ステップS6のYES)、処理部200は超音波振動子毎にフレーム間での信号強度差を算出し、記憶部500にフレーム間信号強度差データ512を新たに記憶する(ステップS8)。
次に、新たに得たフレーム間信号強度差に基づいて、超音波振動子別にループAを実行する(ステップS10〜S20)。
ループAでは、処理対象の超音波振動子の全深さにおける新たに算出された信号強度差の合計を算出し(ステップS12)、算出された信号強度差の合計を加算するように統計処理データ514すなわちヒストグラム(図4(2)参照)を更新する(ステップS18)。
ループAでは、処理対象の超音波振動子の全深さにおける新たに算出された信号強度差の合計を算出し(ステップS12)、算出された信号強度差の合計を加算するように統計処理データ514すなわちヒストグラム(図4(2)参照)を更新する(ステップS18)。
フレーム間信号強度差の算出とループAとは、少なくとも心周期1拍〜数拍分に相当する測定が完了するまで繰り返される(ステップS30のNO)。そして所定時間分の測定が完了すると(ステップS30のYES)、処理部200は次に統計処理データ514から、血管の直上にある超音波振動子を1次選出する(ステップS32)。具体的には、統計処理データ514のヒストグラムの縦軸の値のピーク値を検索し、検索されたピーク値に対応する超音波振動子を、高変化条件を満たす位置に相当するとして選出する。
そして、1次選出された超音波振動子毎に信号強度変化の周期を推定し、ピーク周期データ518として記憶する(ステップS34)。次いで、推定された周期に基づいて、ステップS4で測定開始された心周期に近似している超音波振動子を2次選出する(ステップS36)。これにより、「血管直上にある超音波振動子」が最終判定され、血管上の生体表面の位置を判定する第1の位置判定がなされたことになる。処理部200は、当該超音波振動子の識別情報を含む血管位置データ520を記憶部500に記憶させる。また、報知制御部238が、ステップS36で選出された超音波振動子に対応する位置、すなわち第1の位置判定により判定された位置を報知出力する制御を行う(ステップS37)。
次に、処理部200は2次選出された超音波振動子の反射波信号510から当該素子の直下にある血管の深さ位置(Hc)を算出し、対応する血管位置データ520に格納する(ステップS38)。また、報知制御部238が、算出された深さ位置(Hc)を報知出力する制御を行う(ステップS39)。なお、ステップS37の報知を、ステップS38で纏めて行うこととしてもよい。
これで生体情報算出に必要な血管の位置の判定が完了したので、処理部200は2次選出された超音波振動子を用いて生体情報の算出を実行する(ステップS40)。本実施形態ではIMTを算出するので、例えば、2次選出された超音波振動子に対応する反射波信号510から血管の前壁及び後壁の境界を検出してIMTを算出する。勿論、算出する生体情報によっては、ここで改めて超音波測定を行うとしてもよい。例えば、血管壁のトラッキング処理を行うなどである。
以上、本実施形態によれば、超音波振動子毎に異なるフレーム間での信号強度の時間的変化のヒストグラムを作成するといった極めて演算負荷の低い簡単な処理で、プローブ走査範囲内から血管の位置を特定することができる。従って、超音波測定装置10の小型化および低廉化を図ることができる。
なお、本発明を適用可能な形態は上記実施形態に限るものではなく、適宜構成要素の追加・省略・変更を施すことができる。
例えば、第1位置判定部234が血管の直上の生体表面の位置を判定する方式として次のような方式を採用することとしてもよい。すなわち、生体表面上の同一位置における深さ方向の信号を、受信タイミングが異なる当該信号の組み合わせで自己相関演算した相関値を各位置において求める(信号強度の変化の統計処理)。そして、相関値が所定の低レベル条件を満たす位置を血管直上位置として判定する。
より具体的には、ある超音波振動子において、タイミングが異なる2つの反射波信号を得る。例えば、図3(1)及び図3(2)に示した2つのフレームの反射波信号を得る。この2つの反射波信号を自己相関演算することで相関値を算出する。これを、各超音波振動子18について行うことで、信号強度の変化の統計処理を実現する。図10にその結果の一例を示す。
図10(1)は、図4(1)に示したBモード画像と同一の画像であり、図10(2)は、連続する2フレーム間での各超音波振動子における受信信号(反射波信号)の相関値を示すグラフである。2フレーム間の信号に変化(違い)が無いほど、相関値は高い値となり、変化(違い)があるほど、相関値は低い値となる。血管、特に動脈は心臓の拍動により大きく収縮拡張しているため、相関値は低い値となることが分かる。これを、ある閾値と比較することで、所定の低レベル条件を満たす相関値を判定し、対応する超音波振動子の位置を血管の直上位置として判定することができる。図10(2)では、ある閾値以下であり、且つ、最小の相関値を取った振動子Tr1が血管の直上に位置する超音波振動子とされている。また、2フレーム間だけでなく、ある一定の期間に亘って相関値を連続算出(さらには積算)して、所定の低レベル条件を満たす位置を判定することとしてもよい。
例えば、ステップS4や、ステップS34〜S36を省略した構成とすることができる。ステップS34〜S36を省略した場合では、ステップS38ではステップS32で1次選出された超音波振動子について血管の深さ位置を算出すれば良い。また、ステップS38〜S39を省略した構成としてもよい。この場合、ステップS40では1次選出された超音波振動子を用いて生体情報を算出すれば良い。
また、上記実施形態では、血管の深さ位置を2次選出された超音波振動子の反射波信号510から前壁及び後壁に該当するピークを検索することで決定したが、これに統計的な処理を加えることもできる。
例えば、図9に示すように、ステップS38にて、2次選出された超音波振動子に関する信号強度差(図3(3)に相当)を、深さ位置毎に積算した第2のヒストグラムを作成する。そして、当該第2のヒストグラムのピークを検索する。当然、このピークは前壁と後壁それぞれについて求められるので、これら2つのピークの中間に当たる深さを血管の深さ位置(Hc)とすればよい。
例えば、図9に示すように、ステップS38にて、2次選出された超音波振動子に関する信号強度差(図3(3)に相当)を、深さ位置毎に積算した第2のヒストグラムを作成する。そして、当該第2のヒストグラムのピークを検索する。当然、このピークは前壁と後壁それぞれについて求められるので、これら2つのピークの中間に当たる深さを血管の深さ位置(Hc)とすればよい。
2…被検体、4…血管、4f…前壁、4r…後壁、4s…横壁、10…超音波測定装置、12…タッチパネル、14…キーボード、16…超音波プローブ、18…超音波振動子、30…処理装置、31…制御基板、32…CPU、33…記憶媒体、34…通信IC、100…操作入力部、102…送信部、104…受信部、200…処理部、210…超音波測定制御部、212…照射制御部、214…送受信制御部、216…受信合成部、220…心拍測定部、230…血管位置判定部、232…周期基準判定部、234…第1位置判定部、236…第2位置判定部、238…報知制御部、240…生体情報算出部、260…画像生成部、360…画像表示部、500…記憶部、501…測定プログラム、502…血管位置判定プログラム、510…反射波信号、512…フレーム間信号強度差データ、514…統計処理データ、518…ピーク周期データ、520…血管位置データ、530…生体情報測定結果、As…プローブ走査範囲
Claims (9)
- 生体内組織に対して超音波を照射する制御を行う照射制御部と、
前記超音波の反射波信号を受信する受信部と、
前記反射波信号に含まれる信号の時間的変化に基づいて、前記生体内組織中の血管の位置を判定する血管位置判定部と、
を備えた超音波測定装置。 - 前記血管位置判定部は、前記時間的変化の積算結果が所定の高変化条件を満たす位置に基づいて、前記血管の位置を判定する、
請求項1に記載の超音波測定装置。 - 前記血管位置判定部は、前記時間的変化の発生が所定の周期条件を満たす位置に基づいて、前記血管の位置を判定する、
請求項1又は2に記載の超音波測定装置。 - 前記血管位置判定部は、生体表面上の同一位置における深さ方向の信号の時間的変化を当該同一位置における時間的変化として統計処理することで、前記血管上の生体表面の位置を判定する第1位置判定部を有する、
請求項1〜3の何れか一項に記載の超音波測定装置。 - 前記第1位置判定部は、前記血管上の生体表面の位置において、前記血管の前壁部の位置の時間的変化と後壁部の位置の時間的変化とを当該表面位置における時間的変化として合算して統計処理する、
請求項4に記載の超音波測定装置。 - 前記第1位置判定部は、生体表面上の同一位置における深さ方向の信号を、前記受信部での受信タイミングが異なる当該信号の組み合わせで自己相関演算した相関値を、当該同一位置における時間的変化を示す値とする、
請求項4に記載の超音波測定装置。 - 前記血管位置判定部は、前記第1位置判定部で判定された判定位置における深さ方向の各位置での信号の時間的変化に基づいて、前記血管の深さ位置を判定する第2の位置判定部を有する、
請求項4〜6の何れか一項に記載の超音波測定装置。 - 前記血管位置判定部は、20回/秒以上のサンプリング速度で前記時間的変化を判定する、
請求項1〜7の何れか一項に記載の超音波測定装置。 - 生体内組織に対して超音波を照射することと、
前記超音波の反射波信号を受信することと、
前記反射波信号に含まれる信号の時間的変化に基づいて、前記生体内組織中の血管の位置を判定することと、
を含む超音波測定方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014038976A JP2015013099A (ja) | 2013-06-04 | 2014-02-28 | 超音波測定装置および超音波測定方法 |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013117619 | 2013-06-04 | ||
JP2013117619 | 2013-06-04 | ||
JP2014038976A JP2015013099A (ja) | 2013-06-04 | 2014-02-28 | 超音波測定装置および超音波測定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2015013099A true JP2015013099A (ja) | 2015-01-22 |
Family
ID=52435377
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2014038976A Pending JP2015013099A (ja) | 2013-06-04 | 2014-02-28 | 超音波測定装置および超音波測定方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2015013099A (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016193020A (ja) * | 2015-03-31 | 2016-11-17 | セコム株式会社 | 超音波センサ |
CN109475343A (zh) * | 2016-08-01 | 2019-03-15 | 深圳迈瑞生物医疗电子股份有限公司 | 超声弹性测量显示方法及系统 |
-
2014
- 2014-02-28 JP JP2014038976A patent/JP2015013099A/ja active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016193020A (ja) * | 2015-03-31 | 2016-11-17 | セコム株式会社 | 超音波センサ |
CN109475343A (zh) * | 2016-08-01 | 2019-03-15 | 深圳迈瑞生物医疗电子股份有限公司 | 超声弹性测量显示方法及系统 |
US11331076B2 (en) | 2016-08-01 | 2022-05-17 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Method and system for displaying ultrasonic elastic measurement |
CN109475343B (zh) * | 2016-08-01 | 2024-04-02 | 深圳迈瑞生物医疗电子股份有限公司 | 剪切波弹性成像测量显示方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9456804B2 (en) | Ultrasound measurement apparatus and ultrasound measurement method | |
US9797997B2 (en) | Ultrasonic diagnostic system and system and method for ultrasonic imaging | |
Beulen et al. | Toward noninvasive blood pressure assessment in arteries by using ultrasound | |
JP6173886B2 (ja) | 超音波診断装置、画像処理装置及び画像処理方法 | |
JP4189405B2 (ja) | 超音波診断装置 | |
US20150250448A1 (en) | Ultrasonic measurement apparatus and ultrasonic measurement method | |
US20150245820A1 (en) | Ultrasonic measurement apparatus and ultrasonic measurement method | |
Cinthio et al. | Intra-observer variability of longitudinal displacement and intramural shear strain measurements of the arterial wall using ultrasound noninvasively in vivo | |
US20150238169A1 (en) | Ultrasonic measurement apparatus and ultrasonic measurement method | |
EP2910192A1 (en) | Ultrasonic measurement apparatus and ultrasonic measurement method | |
WO2014091999A1 (ja) | 超音波撮像装置及び方法 | |
JP5771306B1 (ja) | 超音波診断装置 | |
JP2007222533A (ja) | 超音波診断装置及び超音波画像処理方法 | |
US20150209004A1 (en) | Ultrasonic measurement apparatus and ultrasonic measurement method | |
JP2018102589A (ja) | 脈波伝播速度測定装置、血圧測定装置、および脈波伝播速度測定方法 | |
JP2021533913A (ja) | 脈波伝播速度測定を実行するためのシステム及び方法 | |
JP2009039277A (ja) | 超音波診断装置 | |
JP2015013099A (ja) | 超音波測定装置および超音波測定方法 | |
JP4745455B2 (ja) | 超音波診断装置、超音波画像処理装置、及び超音波信号処理プログラム | |
JP5851549B2 (ja) | 超音波診断装置 | |
JP2019000151A (ja) | 超音波診断装置及び関心領域設定方法 | |
US11998385B2 (en) | Methods and systems for investigating blood vessel characteristics | |
WO2014129203A1 (ja) | 超音波診断装置および超音波診断装置の制御方法 | |
JP2008183118A (ja) | 超音波診断装置 | |
JP2019000152A (ja) | 超音波診断装置及びモード遷移方法 |