JP4113377B2 - 超音波診断装置 - Google Patents
超音波診断装置 Download PDFInfo
- Publication number
- JP4113377B2 JP4113377B2 JP2002134888A JP2002134888A JP4113377B2 JP 4113377 B2 JP4113377 B2 JP 4113377B2 JP 2002134888 A JP2002134888 A JP 2002134888A JP 2002134888 A JP2002134888 A JP 2002134888A JP 4113377 B2 JP4113377 B2 JP 4113377B2
- Authority
- JP
- Japan
- Prior art keywords
- displacement
- phase difference
- displacement amount
- amount
- diagnostic apparatus
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、超音波診断装置に係り、特に対象組織の変位量を測定する装置に関する。
【０００２】
【従来の技術】
対象組織の変位量を測定するには、対象組織に超音波を繰り返し送受信して、受信信号の位相を比較し、その位相差から変位量を求める方法が用いられる。図６は、上記方法を説明するため、時間とともに変位する対象組織１１に対して、探触子２２から超音波１３を繰り返し送信し、対象組織１１から超音波１５が反射され探触子２２に戻る様子を示したものである。（ａ），（ｂ）は、それぞれ時刻Ｔ＝ｔ1，ｔ2における対象組織１１の変化状態を示す。対象組織１１は時間経過とともに探触子２２側に変位し、それに応じて対象組織１１の表面と探触子２２の間隔が時間経過とともに小さくなる。
【０００３】
そこで、時刻Ｔ＝ｔ1において超音波１３を送信し、対象組織１１の表面から反射された超音波１５を探触子２２が受信した信号波形を＃１とする。同様に時刻Ｔ＝ｔ2に送信した超音波に対応する受信信号波形を＃２とする。
【０００４】
図７は、横軸に時間、縦軸に振幅をとり、対象組織１１の表面に相当する受信信号波形＃１，＃２を示したものである。それぞれの受信信号波形についての時間軸の原点は、それぞれの超音波の送信時刻を基準に取った。図７に示すように、対象組織１１の変位の大きさに応じ、これらの受信信号＃１，＃２の相互間に位相差が生ずる。例えば受信信号＃１，＃２のそれぞれにつき、その波形の立上り側のゼロクロス点の相互間の差で位相差Δθ21を求めることができる。位相差Δθ21の大きさは、図６から理解できるように、対象組織１１と探触子２２との距離、つまり対象組織１１の変位量に基づくものである。ここで、２πの位相差が超音波の１波長の長さに相当するので、求められた位相差から変位量を算出できる。したがって、対象組織に超音波を繰り返し送受信して、受信信号の位相を比較し、その位相差から刻々と変化する変位量の時間変化を求めることができる。
【０００５】
【発明が解決しようとする課題】
しかし、この方法では、変位量が大きくなり、受信信号相互間の位相差が±π以上となると、以下に述べる位相差に関する公知の「折り返し」の問題が起こり、測定可能な最大変位量に制限が出る。
【０００６】
図８は、位相差に関する折り返しの問題を説明するため、受信信号相互間の位相差が±π以上となる２個の受信信号＃１と＃３を、図７と同じ表し方で示したものである。図８において、受信信号＃１の立上り側ゼロクロス点を点Ｒとすると、信号波形の１波長は２πに相当するので、受信信号＃１を基準とした位相差の識別は、点Ｒを中心とした±πの位相差の範囲に対応する点Ｓと点Ｔの範囲内でのみ有効に行うことができる。いま、対象組織の変位量が大きく、その変位量に対応した受信信号＃３の変化が、受信信号＃１との相互間の位相差Δθ31がπを超えているときは、受信信号＃１に関する位相差の有効識別範囲の点Ｓと点Ｔの範囲外に、受信信号＃３の真のゼロクロス点Ｕがくる。このとき、受信信号＃１に関する位相差の有効識別範囲の点Ｓと点Ｔの範囲内でも、受信信号＃３は他のゼロクロス点Ｖをもつ。したがって、いまの場合、受信信号＃１のゼロクロス点Ｒから位相差Δθ31だけ離れた真のゼロクロス点Ｕでなく、位相差［２π−Δθ31］だけ離れた他のゼロクロス点Ｖを位相差の算出に用いてしまう。
【０００７】
このように、位相差が±πを超えると、観測する位相差の折り返しが起こるために、位相差から変位量を求める方法は、位相差に関し±πの範囲の測定に制限される。したがって、測定可能な最大変位量は、超音波の波長をλとして、±λ／２の範囲に制限され、これを超える大きな変位量を測定することができない。
【０００８】
超音波の周波数を低くして、波長λを大きくすれば大きな変位量の測定が可能になるが、変位量の測定精度が低くなる。
【０００９】
一方、超音波の送信する間隔を短くし、変位量の測定時刻ｔ1，ｔ2の間隔を短くすれば、送信間隔内の対象組織の変位量が小さくなるので、位相差を±πの範囲に納めるようにすることができる。しかし、超音波のパルス繰り返し時間（Ｐｕｌｓｅ Ｒｅｐｅｔｉｔｉｏｎ Ｔｉｍｅ：ＰＲＴ）を短くすることになり、周知のように、対象組織の診断深さが浅くなる制限が出る。
【００１０】
また、対象組織の測定点、例えば血管の前壁にトラッキングゲートを設定し、血管壁の変位を逐次検出するトラッキング機能において、位相検出を用いて血管前壁の変位の逐次検出を行うときは同様の問題が起こることがある。すなわち、この場合のトラッキングゲートの変位計測可能な窓幅は、超音波の波長をλとして、±λ／２であるので、血管前壁が±λ／２以上変位するときは、トラッキングゲートはその変位に追随できないか、あるいは誤った動作をする。
【００１１】
本発明の目的は、かかる従来技術の課題を解決し、変位量の大きい対象組織の変位を精度よく測定する超音波診断装置を提供することである。他の目的は、対象組織の診断深さを変えることなく、変位量の大きい対象組織の変位を精度よく測定する超音波診断装置を提供することである。
【００１２】
【課題を解決するための手段】
上記目的を達成するため、本発明に係る超音波診断装置は、対象組織に対し、送受信シーケンスに従って、第１周波数の超音波を繰り返し送受信して複数の第１受信信号を出力し、前記第１周波数より低い第２周波数の超音波を繰り返し送受信して複数の第２受信信号を出力する送受信手段と、前記複数の第１受信信号についての相互間の第１位相差と、前記複数の第２受信信号についての相互間の第２位相差とに基づいて、対象組織の変位量を演算する変位量演算手段と、を備える。
【００１３】
この構成により、短い波長の第１周波数の超音波による計測と、長い波長の第２周波数の超音波による計測の両者の利点を生かして、変位量の大きい対象組織の変位を精度よく測定することができる。変位量は、位置、速度などの測定量を含む。
【００１４】
ここで、前記変位量演算手段は、前記複数の第１受信信号の相互間の波形比較により前記第１位相差を演算する第１位相差演算部と、前記複数の第２受信信号の相互間の波形比較により前記第２位相差を演算する第２位相差演算部と、前記第２位相差に基づいて、前記第１位相差に関する折り返しの有無および回数を判定する判定部と、前記判定部の判定結果に従い、前記第１位相差に基づいて前記対象組織の変位量を演算する変位量演算部と、を備える。
【００１５】
この構成により、第１周波数の超音波についてその位相差の折り返しにより識別できない部分があるか否か、あるときは何回折り返しがあるかを判別する。その結果を用いることで、第１周波数の超音波のＰＲＴを変えることなく、第１周波数の超音波について位相差の折り返しで制限される測定可能な最大変位量を超える分を、第２周波数の超音波を用いた位相差の測定で補うことができる。したがって対象組織の診断深さを変えることなく、変位量の大きい対象組織の変位を精度よく測定することができる。
【００１６】
望ましくは、前記送受信シーケンスは、前記第１周波数の超音波と前記第２周波数の超音波とを交互に送信する送受信シーケンスである。
【００１７】
この構成により、第１位相差と第２位相差とを相互に近接した時刻で測定できるので、それらの測定値の同時性がよく、対象組織の変位量を正確に演算できる。また、対象組織の変位量を短時間で演算できる。
【００１８】
望ましくは、前記折り返し回数は、前記第１位相差の１回の折り返し当たりの単位変位量で前記第２位相差に相当する参照変位量を除することにより算出される整数である。
【００１９】
望ましくは、前記変位量演算部は、前記第１位相差から演算される見かけの変位量と、前記折り返し回数に前記単位変位量を乗じて得られる折り返し変位量とに基づいて前記対象組織の変位量を演算する。
【００２０】
ここで見かけ変位量は、位相差の折り返しが生じている状態での第１位相差に基づいて演算される見かけの変位量である。したがって、見かけ変位量と、折り返しの生じた大きな変位量すなわち折り返し変位量とから、対象組織の変位量を得ることができる。
【００２１】
望ましくは、前記見かけ変位量は、前記位相差の折り返しが生じない変位量観測範囲である変位計測窓を、変位方向に前記折り返し変位量分だけシフトさせて演算される。
【００２２】
ここで変位計測窓は、長さの次元では、第１超音波の波長をλとして、±λ／２の幅を有するものとできる。また、±λ／２の幅に対応する時間幅または位相幅であってもよい。
【００２３】
望ましくは、前記見かけ変位量は、前記変位計測窓のシフトにあわせ、上記変位計測窓内における変位前の第１受信信号の変位量測定のための特徴点の位置を、変位方向に前記折り返し変位量分だけシフトさせ、上記シフトさせた変位計測窓内における変位前後の前記特徴点の位置の差から演算される。
【００２４】
ここで特徴点の位置は、第１受信信号のゼロクロス点の位置とすることができる。また、第１受信信号の他の点、例えばピーク点の位置等を特徴点の位置とすることもできる。位置は、距離についての位置の他、時間、位相についての位置でもよい。
【００２５】
【発明の実施の形態】
以下に図面を用い、本発明に係る実施の形態につき詳細に説明する。図１は、本発明の原理説明図である。（ａ）は、縦軸に振幅、横軸に時間を取り、変位前の対象組織に対し第１周波数の超音波を送信して得られた第１受信信号１０２ａを実線で、変位後の対象組織に対し第１周波数の超音波を送信して得られた第１受信信号１０２ｂを破線で示してある。同様に、（ｂ）は、第１周波数より低い第２周波数の超音波を変位前の対象組織に対し送信して得られた第２受信信号１０４ａを実線で、変位後の対象組織に対し第２周波数の超音波を送信して得られた第２受信信号１０４ｂを破線で示してある。（ａ）と（ｂ）の時間軸の原点は、同じ位置として配列した。
【００２６】
すでに図７で説明したように、対象組織の変位量は、変位前の受信信号のゼロクロス点と変位後の受信信号のゼロクロス点との相互間の変化から求められる。図１（ａ）の第１受信信号１０２ａ，１０２ｂの相互間においては、変位前後でゼロクロス点が点Ｒから点Ｕへ変化し、図１（ｂ）の第２受信信号１０４ａ，１０４ｂの相互間においては、変位前後でゼロクロス点が点Ｐから点Ｑへ変化している。ここで、第１受信信号１０２ａ，１０２ｂ相互間におけるゼロクロス点の変化は±πを超えているので、図８で説明した位相差の折り返し問題が起こっている。したがって、真のゼロクロス点はＵであるが、従来技術では、位相差の折り返しによって現れる他のゼロクロス点Ｖを観測してしまい、ゼロクロス点Ｖに基づいて位相差の算出を行ってしまう。
【００２７】
第１周波数の超音波の変位量計測においては、位相差の折り返しが起こらない範囲、つまり±λ／２の範囲で高精度の計測ができる。したがって、±λ／２の範囲が、変位量を有効に観測できる変位計測窓である。
【００２８】
これに対し、第２周波数は第１周波数より周波数が低く、（ｂ）に示すように第２受信信号におけるゼロクロス点の変化について位相差の折り返しが起こっていない。そこで、変位前後の２つの第２受信信号のゼロクロス点の変化、すなわち点Ｐと点Ｑとの相互間の位相差から変位量を演算することができる。この変位量を参照変位量Ｄと呼ぶことにする。参照変位量は、全体変位量を表しているが、周波数の低い第２周波数の超音波を用いて算出するので、第１周波数の超音波を用いて算出する値に比べて、精度が低下する。
【００２９】
図１（ａ）において、上記のように点Ｒと点Ｕとの間の位相差はπを超え、位相差の折り返しが起こっているので、位相差の折り返しの回数を求める。位相差の折り返しは、±πを単位として起こる。変位量の次元で換算すると±λ／２を単位として起こる。したがって、参照変位量Ｄをλ／２で除した整数をｍとすると、ｍが位相差の折り返し回数になる。図示した例では、λは第１周波数の超音波の波長であるので、（ａ）と（ｂ）を比較して、位相差の折り返し回数ｍは１である。
【００３０】
いま、位相差の折り返し回数をｍとしてｍ＊（±λ／２）に相当する変位量を、折り返し回数に相当する変位量という意味で折り返しで生じた変位量あるいは単に折り返し変位量ということにする。
【００３１】
全変位量から折り返しで生じた変位量を差し引いた変位量の大きさをΔｄと表すとすると、Δｄは±λ／２の範囲の大きさ、すなわち第１受信信号の変位計測窓内の大きさであり、この部分は高精度の変位量計測ができる。Δｄは高精度ではあるが、上記のように全変位量を表さず、±λ／２の範囲の成分のみを表す、いわば見かけの変位量である。
【００３２】
見かけ変位量Δｄは、変位計測窓をｍ＊（λ／２）だけ変位方向にシフトすることで計測できる。つまり、図１（ａ）において、変位前の変位計測窓１０６は、点Ｔと点Ｓの範囲である。この点Ｔと点Ｓの範囲の変位計測窓１０６内では、変位後の第１受信信号のゼロクロス点は点Ｖとして観測される。すでに述べたように点Ｖは、変位量を示す真のゼロクロス点ではない。そこで、点Ｔと点Ｓの範囲の変位計測窓１０６をｍ＊（λ／２）だけ変位方向にシフトする。
【００３３】
図１において、変位方向は時間軸で進み方向（右方向）であり、ｍ＝１であるから、変位計測窓１０６をλ／２だけ時間軸に沿い進み方向にシフトする。シフト後の変位計測窓１０８内で、変位後の第１受信信号のゼロクロス点は点Ｕとして観測される。すなわち、変位計測窓をｍ＊（λ／２）だけ変位方向にシフトすることで、シフト後の変位計測窓１０８内で、真のゼロクロス点Ｕを検出できる。
【００３４】
検出された真のゼロクロス点Ｕを用いて、見かけ変位量Δｄを算出するには、変位前の第１受信信号のゼロクロス点Ｒの位置も変位計測窓のシフトとともにシフトすればよい。図１（ａ）の場合、ゼロクロス点Ｒの位置を、変位方向にλ／２だけシフトすると、点Ｓの位置になる。そこで、図７の説明に従い、点Ｓと点Ｕの間の位相差から見かけ変位量Δｄを求めることができる。
【００３５】
したがって、対象組織の変位量は、折り返しで生じた変位量ｍ＊（λ／２）に、それを超えた±λ／２の範囲の見かけ変位量Δｄを加算して求めることができる。図１の場合ではｍ＝１であるから、対象組織の変位量は、（λ／２）＋Δｄとして求められる。このように、変位量の大きい場合、周波数の高い第１受信信号で位相差の折り返しが起こっても、周波数の低い第２受信信号の参照変位量Ｄを用い、その折り返し回数ｍを算出して補うことで、生体組織の変位量を精度よく測定できる。
【００３６】
なお、点Ｖは変位量を示す真のゼロクロス点ではないが、点Ｖと点Ｒの位相差は、（λ／２）−Δｄの値に相当するので、点Ｔと点Ｓの範囲の変位計測窓をシフトせずに、その範囲で検出されるゼロクロス点Ｖを用いて（λ／２）−Δｄを求め、この値から−Δｄを演算することもできる。この場合、Δｄに付される符号が正負逆になる。
【００３７】
また、第１受信信号の変位前のゼロクロス点Ｒを変位計測窓のシフトとともにシフトする他に、変位計測窓の基準位置、例えば時間軸で最も遅れ側の点（図１では変位計測窓１０６，１０８のそれぞれ左端）を基準位置と定め、その基準位置からのゼロクロス点の位置の変化からΔｄを求めることもできる。すなわち、変位前の変位計測窓１０６の左端から変位前のゼロクロス点Ｒの位置までの量と、変位後の変位計測窓１０８の左端と変位後の真のゼロクロス点Ｕ点の位置までの量を比較して、Δｄを求めることもできる。
【００３８】
図２は、超音波診断装置２０のブロック図である。図２において、探触子２２は、超音波の送波およびエコーの受波を行う超音波探触子で、例えば、被検者の体表面上に当接して用いられる。探触子２２内には図示されていないアレイ振動子が設けられる。アレイ振動子は複数の振動素子からなり、アレイ振動子により超音波ビームが形成される。この超音波ビームは電子走査され、その電子走査方式としては、例えば電子リニア走査や電子セクタ走査を用いることができる。
【００３９】
送受信部２４は、送信回路２６と受信回路２８とを備え、探触子２２と接続される。送受信部は、後述する送受信制御回路３２の制御の下で、超音波の送受信により、断層画像用の超音波ビームと、これとは別に変位量計測用の２つの異なる周波数の超音波ビームを形成し、受信信号を出力する回路である。より詳しくは、送信回路２６は、アレイ振動子の各チャネルごとに遅延された送信信号を供給する回路である。受信回路２８は、探触子２２からのエコー信号を増幅し、各チャネル間の受信信号の位相差を調整する整相加算等の処理を行い、受信信号としてＢモード信号処理部３６と変位量演算部３８に出力する回路である。
【００４０】
送受信制御回路３２の制御の下で、探触子２２から断層画像用の超音波ビームが生体内に走査される走査面１０の様子を図２中に示す。後述するように、断層画像用の超音波のエコー信号は信号処理および表示処理がなされ、走査面１０に断層画像として表示される。ユーザは、この断層画像を見ながら、対象組織１１の変位量計測個所、例えば血管の前壁点Ａにつき、点Ａを通過する方位を設定できる。その設定に基づき送受信制御回路３２の制御の下で、探触子２２からその方位に第１周波数の超音波１２の送受信と、第１周波数より低い第２周波数の超音波１４の送受信が交互に行われる。例えば第１周波数を８ＭＨｚ、第２周波数を４ＭＨｚとすることができる。
【００４１】
図３は、送受信制御回路３２により制御される第１周波数の超音波と第２周波数の超音波の送信タイミングチャートである。図３は横軸に時間をとり、上段に同期クロックを、中段に第１周波数の超音波の送信タイミングを、下段に第２周波数の超音波の送信タイミングを示した。このように、第１周波数の超音波の送信と第２周波数の超音波の送信は、一定の送信間隔で交互に行われる。
【００４２】
再び図２に戻り、Ｂモード信号処理部３６は、断層画像用の超音波について受信回路２８から出力される整相加算後の受信信号に基づき、エコー信号の包絡振幅を抽出する検波、包絡振幅信号の対数圧縮等の処理を行う回路である。信号処理の結果は、表示処理部４２に出力される。
【００４３】
変位量演算部３８は、変位量計測用の超音波について受信回路２８から受取った２つの異なる周波数の超音波の受信信号から、対象組織１１の変位量を演算する機能を有する。具体的には、後述するユーザにより設定されたトラッキングゲートによって、対象組織１１である血管前壁の点Ａの位置をトラッキングし、点Ａの変位量を演算する機能を有する。トラッキング機能は、後述する制御部３０の制御の下において、設定されたトラッキングゲート内で、位相を検出し、血管前壁の点Ａの変位を逐次検出する機能である。変位量を演算する機能のさらに詳細については後に詳述する。変位量演算の結果は表示処理部４２に出力される。
【００４４】
トラッキングゲートの位相検出に基づく変位計測窓の範囲は、超音波の波長をλとして±λ／２なので、この範囲を超えるときは、位相差の折り返しの問題が起こる。したがって、制御部３０の制御の下で、図１で説明したと同様に、生じた位相差の折り返し回数をｍとして、ｍ＊（λ／２）だけ変位方向に変位計測窓をシフトさせることで、±λ／２の範囲を超える変位量の場合でも誤動作することなくトラッキング機能を行うことができる。
【００４５】
表示処理部４２は、Ｂモード信号処理部３６の出力からＢモード断層画像を形成し、変位量演算部３８の演算結果を示す数値やグラフの表示に必要な処理を行う回路である。変位量に関するグラフとしては、例えば変位量の時間変化を示すグラフ等である。また、これらの画像を合成する機能も有している回路である。表示処理部４２は、いわゆるディジタルスキャンコンバータ（ＤＳＣ）や各種の画像処理回路によって構成することができる。形成されたＢモード断層画像と変位量に関する画像は、表示器４４に出力される。
【００４６】
表示器４４は、表示処理部４２により形成されたＢモード断層画像や変位量の演算結果を示す数値やグラフを表示する。また、例えば血管前壁の変位量の他に、血管後壁の変位量を計測し、その結果をあわせて表示することもできる。さらに、血管前壁の変位量と血管後壁の変位量から血管の直径を演算し、その結果を表示することもできる。
【００４７】
入力部４６は、後述の制御部３０と接続され、ユーザにより、超音波診断装置２０の動作モードの選択や必要なデータの入力が行われるスイッチ、キーボード、トラックボール等である。例えば、ユーザは、入力部４６を用いて、動作モードとして変位計測を選択できる。また、ユーザは、表示器４４の画面を見ながらトラックボール等を操作して断層画像上で変位量計測を行う点Ａを通過する所定の方位を設定でき、点Ａの位置に変位量計測のためのトラッキングゲートを設定できる。
【００４８】
制御部３０は、超音波診断装置２０の各構成と接続され、それぞれの制御を行うコントローラである。制御部３０は、上記の送受信部２４を制御する送受信制御回路３２を備える。
【００４９】
制御部３０は、変位量計測のための超音波ビームの方位を設定する制御を行う機能を有する。すなわち、入力部４６から入力された所定の方位に、第１周波数の超音波と第１周波数より低い第２周波数の超音波の方位を設定する制御を行う機能を有する。また、制御部３０は、変位量計測のためのトラッキングを制御する機能を有する。すなわち、入力部４６からなされたトラッキングゲートの設定に基づき、血管前壁の点Ａの変位を逐次検出するトラッキングを制御する機能を有する。
【００５０】
図４は、変位量演算部３８の詳細ブロック図である。変位量演算部３８は、受信回路２８が繰り返し出力する複数の第１受信信号５０を記憶する第１受信信号メモリ５４、複数の第２受信信号５２を記憶する第２受信信号メモリ５６、および変位量演算回路５８を備え、変位前後の２つの第１受信信号と変位前後の２つの第２受信信号の計４個の受信信号に基づき、対象組織の変位量６０を演算し出力する機能を有する。
【００５１】
第１受信信号メモリ５４は、受信回路２８から受取った複数の第１受信信号５０を記憶する半導体メモリ等の記憶装置である。同様に第２受信信号メモリ５６は、受信回路２８から受取った複数の第２受信信号５２を記憶する記憶装置である。各受信信号の受信順序の対応付けは、例えばメモリのアドレス等に対応付けて管理される。
【００５２】
変位量演算回路５８の中の参照変位量演算器６２は、第２受信信号メモリ５６から変位前後の２つの第２受信信号を読み出し、図７で説明した方法に従い、その信号相互間の位相差である第２位相差を求め、その第２位相差に基づき参照変位量Ｄを演算する機能を有する。
【００５３】
変位量判定器６４は、演算された参照変位量Ｄに基づき、図１で説明した方法に従い、位相差の折り返しの有無とその回数を判定する機能を有する。すなわち、参照変位量Ｄを第１周波数の超音波の波長をλとしてλ／２で除したときに得られる整数部分ｍが、位相差の折り返しの回数である。判定式の形で示すと、次式を満たす整数ｍが位相差の折り返し回数である。ｍ＝０のときは位相差の折り返しが生じていない。
【数１】
ｍ＊（λ／２）≦Ｄ≦（ｍ＋１）＊（λ／２） ・・・・（１）
【００５４】
見かけ変位量演算器６６は、図１で説明した方法に従い、判定された位相差の折り返し回数ｍに従って、第１受信信号の変位計測窓をシフトし、見かけ変位量Δｄを演算する機能を有する。図１で説明した方法を再言すれば、まず変位前後の２つの第１受信信号を第１受信信号メモリ５４から読出した後、第１受信信号の変位計測窓を、折り返し変位量ｍ＊（λ／２）だけ変位方向にシフトし、その状態で変位後の第１受信信号のゼロクロス点Ｕの位置を求める。そして、変位前の第１受信信号のゼロクロス点Ｒの位置を、折り返し変位量ｍ＊（λ／２）だけ変位方向にシフトした点Ｓと点Ｕとの相互間の位相差から見かけ変位量Δｄが算出される。
【００５５】
対象組織変位量演算器６８は、位相差の折り返し回数ｍに従って、対象組織の変位量を演算する機能を有する。すなわち、ｍが０以外の整数のときは、見かけ変位量Δｄと、折り返し変位量ｍ＊（λ／２）とを加算し、ｍ＊（λ／２）＋Δｄを演算し、その値を対象組織の変位量６０として表示処理部４２に出力する機能を有する。また、ｍが０のときは、変位前後の２つの第１受信信号を第１受信信号メモリ５４から読出し、図７で説明した方法に従い、変位前後の第１受信信号の位相差から求めた変位量をそのまま対象組織の変位量として表示処理部４２に出力する機能を有する。
【００５６】
図５は、対象組織の変位量演算の手順をフローチャートで示したものである。
【００５７】
まず受信した第１受信信号を第１受信信号メモリに順次記憶する（Ｓ１０）。同様に受信した第２受信信号を第２受信信号メモリに順次記憶する（Ｓ１２）。
【００５８】
次に第２受信信号メモリから変位前後の２つの第２受信信号を読出し、参照変位量演算器において、その相互間の位相差から参照変位量Ｄを演算する（Ｓ１４）。
【００５９】
そして、変位量判定器で、参照変位量Ｄを用いて、第１周波数の位相差の折り返しの有無と、その回数ｍを判定する（Ｓ１６）。折り返しの回数ｍは、式（１）が成立するｍを算出することで求めることができる。
【００６０】
次に、見かけ変位量演算器は、第１受信信号メモリから、変位前後の２つの第１受信信号を読出す（Ｓ１８）。そして、ｍが０以外の整数のときは、図１で説明した方法に従い、変位前の第１受信信号の変位計測窓をｍ＊（λ／２）だけ変位方向にシフトさせた上で、変位前後の第１受信信号の相互間の位相差を求め、その位相差における見かけ変位量Δｄを演算する（Ｓ２０）。
【００６１】
そして、ｍ＊（λ／２）＋Δｄを演算する（Ｓ２２）。この値を、対象組織の変位量とすることができる。
【００６２】
上記の加算、乗算等の四則演算処理は、演算処理装置を用いて行うことができる。また、演算処理装置に代わり、数式化したメモリを用い、必要なパラメータ、例えば折り返し回数ｍ、λ、参照変位量Ｄ等をメモリに入力して数式を処理し、同じ結果を得ることもできる。
【００６３】
【発明の効果】
本発明に係る超音波診断装置によれば、対象組織の診断深さを変えることなく、変位量の大きい対象組織の変位を精度よく測定することができる。
【図面の簡単な説明】
【図１】 本発明の原理説明図である。
【図２】 本発明に係る実施の形態の超音波診断装置のブロック図である。
【図３】 第１周波数の超音波と第２周波数の超音波の送信タイミングチャートである。
【図４】 変位量演算部の詳細ブロック図である。
【図５】 対象組織の変位量演算のフローチャートである。
【図６】 位相差から変位量を求める原理の説明図である。
【図７】 ２つの受信信号相互間の位相差の関係を示す図である。
【図８】 位相差に関する折り返しの問題を説明する図である。
【符号の説明】
１０ 走査面、１１ 対象組織、１２ 第１周波数の超音波、１４ 第２周波数の超音波、２０ 超音波診断装置、２２ 探触子、２４ 送受信部、２６ 送信回路、２８ 受信回路、３０ 制御部、３２ 送受信制御回路、３６ Ｂモード信号処理部、３８ 変位量演算部、５４ 第１受信信号メモリ、５６ 第２受信信号メモリ、５８ 変位量演算回路、６２ 参照変位量演算器、６４ 変位量判定器、６６ 見かけ変位量演算器、６８ 対象組織変位量演算器、１０６，１０８ 変位計測窓。
【発明の属する技術分野】
本発明は、超音波診断装置に係り、特に対象組織の変位量を測定する装置に関する。
【０００２】
【従来の技術】
対象組織の変位量を測定するには、対象組織に超音波を繰り返し送受信して、受信信号の位相を比較し、その位相差から変位量を求める方法が用いられる。図６は、上記方法を説明するため、時間とともに変位する対象組織１１に対して、探触子２２から超音波１３を繰り返し送信し、対象組織１１から超音波１５が反射され探触子２２に戻る様子を示したものである。（ａ），（ｂ）は、それぞれ時刻Ｔ＝ｔ1，ｔ2における対象組織１１の変化状態を示す。対象組織１１は時間経過とともに探触子２２側に変位し、それに応じて対象組織１１の表面と探触子２２の間隔が時間経過とともに小さくなる。
【０００３】
そこで、時刻Ｔ＝ｔ1において超音波１３を送信し、対象組織１１の表面から反射された超音波１５を探触子２２が受信した信号波形を＃１とする。同様に時刻Ｔ＝ｔ2に送信した超音波に対応する受信信号波形を＃２とする。
【０００４】
図７は、横軸に時間、縦軸に振幅をとり、対象組織１１の表面に相当する受信信号波形＃１，＃２を示したものである。それぞれの受信信号波形についての時間軸の原点は、それぞれの超音波の送信時刻を基準に取った。図７に示すように、対象組織１１の変位の大きさに応じ、これらの受信信号＃１，＃２の相互間に位相差が生ずる。例えば受信信号＃１，＃２のそれぞれにつき、その波形の立上り側のゼロクロス点の相互間の差で位相差Δθ21を求めることができる。位相差Δθ21の大きさは、図６から理解できるように、対象組織１１と探触子２２との距離、つまり対象組織１１の変位量に基づくものである。ここで、２πの位相差が超音波の１波長の長さに相当するので、求められた位相差から変位量を算出できる。したがって、対象組織に超音波を繰り返し送受信して、受信信号の位相を比較し、その位相差から刻々と変化する変位量の時間変化を求めることができる。
【０００５】
【発明が解決しようとする課題】
しかし、この方法では、変位量が大きくなり、受信信号相互間の位相差が±π以上となると、以下に述べる位相差に関する公知の「折り返し」の問題が起こり、測定可能な最大変位量に制限が出る。
【０００６】
図８は、位相差に関する折り返しの問題を説明するため、受信信号相互間の位相差が±π以上となる２個の受信信号＃１と＃３を、図７と同じ表し方で示したものである。図８において、受信信号＃１の立上り側ゼロクロス点を点Ｒとすると、信号波形の１波長は２πに相当するので、受信信号＃１を基準とした位相差の識別は、点Ｒを中心とした±πの位相差の範囲に対応する点Ｓと点Ｔの範囲内でのみ有効に行うことができる。いま、対象組織の変位量が大きく、その変位量に対応した受信信号＃３の変化が、受信信号＃１との相互間の位相差Δθ31がπを超えているときは、受信信号＃１に関する位相差の有効識別範囲の点Ｓと点Ｔの範囲外に、受信信号＃３の真のゼロクロス点Ｕがくる。このとき、受信信号＃１に関する位相差の有効識別範囲の点Ｓと点Ｔの範囲内でも、受信信号＃３は他のゼロクロス点Ｖをもつ。したがって、いまの場合、受信信号＃１のゼロクロス点Ｒから位相差Δθ31だけ離れた真のゼロクロス点Ｕでなく、位相差［２π−Δθ31］だけ離れた他のゼロクロス点Ｖを位相差の算出に用いてしまう。
【０００７】
このように、位相差が±πを超えると、観測する位相差の折り返しが起こるために、位相差から変位量を求める方法は、位相差に関し±πの範囲の測定に制限される。したがって、測定可能な最大変位量は、超音波の波長をλとして、±λ／２の範囲に制限され、これを超える大きな変位量を測定することができない。
【０００８】
超音波の周波数を低くして、波長λを大きくすれば大きな変位量の測定が可能になるが、変位量の測定精度が低くなる。
【０００９】
一方、超音波の送信する間隔を短くし、変位量の測定時刻ｔ1，ｔ2の間隔を短くすれば、送信間隔内の対象組織の変位量が小さくなるので、位相差を±πの範囲に納めるようにすることができる。しかし、超音波のパルス繰り返し時間（Ｐｕｌｓｅ Ｒｅｐｅｔｉｔｉｏｎ Ｔｉｍｅ：ＰＲＴ）を短くすることになり、周知のように、対象組織の診断深さが浅くなる制限が出る。
【００１０】
また、対象組織の測定点、例えば血管の前壁にトラッキングゲートを設定し、血管壁の変位を逐次検出するトラッキング機能において、位相検出を用いて血管前壁の変位の逐次検出を行うときは同様の問題が起こることがある。すなわち、この場合のトラッキングゲートの変位計測可能な窓幅は、超音波の波長をλとして、±λ／２であるので、血管前壁が±λ／２以上変位するときは、トラッキングゲートはその変位に追随できないか、あるいは誤った動作をする。
【００１１】
本発明の目的は、かかる従来技術の課題を解決し、変位量の大きい対象組織の変位を精度よく測定する超音波診断装置を提供することである。他の目的は、対象組織の診断深さを変えることなく、変位量の大きい対象組織の変位を精度よく測定する超音波診断装置を提供することである。
【００１２】
【課題を解決するための手段】
上記目的を達成するため、本発明に係る超音波診断装置は、対象組織に対し、送受信シーケンスに従って、第１周波数の超音波を繰り返し送受信して複数の第１受信信号を出力し、前記第１周波数より低い第２周波数の超音波を繰り返し送受信して複数の第２受信信号を出力する送受信手段と、前記複数の第１受信信号についての相互間の第１位相差と、前記複数の第２受信信号についての相互間の第２位相差とに基づいて、対象組織の変位量を演算する変位量演算手段と、を備える。
【００１３】
この構成により、短い波長の第１周波数の超音波による計測と、長い波長の第２周波数の超音波による計測の両者の利点を生かして、変位量の大きい対象組織の変位を精度よく測定することができる。変位量は、位置、速度などの測定量を含む。
【００１４】
ここで、前記変位量演算手段は、前記複数の第１受信信号の相互間の波形比較により前記第１位相差を演算する第１位相差演算部と、前記複数の第２受信信号の相互間の波形比較により前記第２位相差を演算する第２位相差演算部と、前記第２位相差に基づいて、前記第１位相差に関する折り返しの有無および回数を判定する判定部と、前記判定部の判定結果に従い、前記第１位相差に基づいて前記対象組織の変位量を演算する変位量演算部と、を備える。
【００１５】
この構成により、第１周波数の超音波についてその位相差の折り返しにより識別できない部分があるか否か、あるときは何回折り返しがあるかを判別する。その結果を用いることで、第１周波数の超音波のＰＲＴを変えることなく、第１周波数の超音波について位相差の折り返しで制限される測定可能な最大変位量を超える分を、第２周波数の超音波を用いた位相差の測定で補うことができる。したがって対象組織の診断深さを変えることなく、変位量の大きい対象組織の変位を精度よく測定することができる。
【００１６】
望ましくは、前記送受信シーケンスは、前記第１周波数の超音波と前記第２周波数の超音波とを交互に送信する送受信シーケンスである。
【００１７】
この構成により、第１位相差と第２位相差とを相互に近接した時刻で測定できるので、それらの測定値の同時性がよく、対象組織の変位量を正確に演算できる。また、対象組織の変位量を短時間で演算できる。
【００１８】
望ましくは、前記折り返し回数は、前記第１位相差の１回の折り返し当たりの単位変位量で前記第２位相差に相当する参照変位量を除することにより算出される整数である。
【００１９】
望ましくは、前記変位量演算部は、前記第１位相差から演算される見かけの変位量と、前記折り返し回数に前記単位変位量を乗じて得られる折り返し変位量とに基づいて前記対象組織の変位量を演算する。
【００２０】
ここで見かけ変位量は、位相差の折り返しが生じている状態での第１位相差に基づいて演算される見かけの変位量である。したがって、見かけ変位量と、折り返しの生じた大きな変位量すなわち折り返し変位量とから、対象組織の変位量を得ることができる。
【００２１】
望ましくは、前記見かけ変位量は、前記位相差の折り返しが生じない変位量観測範囲である変位計測窓を、変位方向に前記折り返し変位量分だけシフトさせて演算される。
【００２２】
ここで変位計測窓は、長さの次元では、第１超音波の波長をλとして、±λ／２の幅を有するものとできる。また、±λ／２の幅に対応する時間幅または位相幅であってもよい。
【００２３】
望ましくは、前記見かけ変位量は、前記変位計測窓のシフトにあわせ、上記変位計測窓内における変位前の第１受信信号の変位量測定のための特徴点の位置を、変位方向に前記折り返し変位量分だけシフトさせ、上記シフトさせた変位計測窓内における変位前後の前記特徴点の位置の差から演算される。
【００２４】
ここで特徴点の位置は、第１受信信号のゼロクロス点の位置とすることができる。また、第１受信信号の他の点、例えばピーク点の位置等を特徴点の位置とすることもできる。位置は、距離についての位置の他、時間、位相についての位置でもよい。
【００２５】
【発明の実施の形態】
以下に図面を用い、本発明に係る実施の形態につき詳細に説明する。図１は、本発明の原理説明図である。（ａ）は、縦軸に振幅、横軸に時間を取り、変位前の対象組織に対し第１周波数の超音波を送信して得られた第１受信信号１０２ａを実線で、変位後の対象組織に対し第１周波数の超音波を送信して得られた第１受信信号１０２ｂを破線で示してある。同様に、（ｂ）は、第１周波数より低い第２周波数の超音波を変位前の対象組織に対し送信して得られた第２受信信号１０４ａを実線で、変位後の対象組織に対し第２周波数の超音波を送信して得られた第２受信信号１０４ｂを破線で示してある。（ａ）と（ｂ）の時間軸の原点は、同じ位置として配列した。
【００２６】
すでに図７で説明したように、対象組織の変位量は、変位前の受信信号のゼロクロス点と変位後の受信信号のゼロクロス点との相互間の変化から求められる。図１（ａ）の第１受信信号１０２ａ，１０２ｂの相互間においては、変位前後でゼロクロス点が点Ｒから点Ｕへ変化し、図１（ｂ）の第２受信信号１０４ａ，１０４ｂの相互間においては、変位前後でゼロクロス点が点Ｐから点Ｑへ変化している。ここで、第１受信信号１０２ａ，１０２ｂ相互間におけるゼロクロス点の変化は±πを超えているので、図８で説明した位相差の折り返し問題が起こっている。したがって、真のゼロクロス点はＵであるが、従来技術では、位相差の折り返しによって現れる他のゼロクロス点Ｖを観測してしまい、ゼロクロス点Ｖに基づいて位相差の算出を行ってしまう。
【００２７】
第１周波数の超音波の変位量計測においては、位相差の折り返しが起こらない範囲、つまり±λ／２の範囲で高精度の計測ができる。したがって、±λ／２の範囲が、変位量を有効に観測できる変位計測窓である。
【００２８】
これに対し、第２周波数は第１周波数より周波数が低く、（ｂ）に示すように第２受信信号におけるゼロクロス点の変化について位相差の折り返しが起こっていない。そこで、変位前後の２つの第２受信信号のゼロクロス点の変化、すなわち点Ｐと点Ｑとの相互間の位相差から変位量を演算することができる。この変位量を参照変位量Ｄと呼ぶことにする。参照変位量は、全体変位量を表しているが、周波数の低い第２周波数の超音波を用いて算出するので、第１周波数の超音波を用いて算出する値に比べて、精度が低下する。
【００２９】
図１（ａ）において、上記のように点Ｒと点Ｕとの間の位相差はπを超え、位相差の折り返しが起こっているので、位相差の折り返しの回数を求める。位相差の折り返しは、±πを単位として起こる。変位量の次元で換算すると±λ／２を単位として起こる。したがって、参照変位量Ｄをλ／２で除した整数をｍとすると、ｍが位相差の折り返し回数になる。図示した例では、λは第１周波数の超音波の波長であるので、（ａ）と（ｂ）を比較して、位相差の折り返し回数ｍは１である。
【００３０】
いま、位相差の折り返し回数をｍとしてｍ＊（±λ／２）に相当する変位量を、折り返し回数に相当する変位量という意味で折り返しで生じた変位量あるいは単に折り返し変位量ということにする。
【００３１】
全変位量から折り返しで生じた変位量を差し引いた変位量の大きさをΔｄと表すとすると、Δｄは±λ／２の範囲の大きさ、すなわち第１受信信号の変位計測窓内の大きさであり、この部分は高精度の変位量計測ができる。Δｄは高精度ではあるが、上記のように全変位量を表さず、±λ／２の範囲の成分のみを表す、いわば見かけの変位量である。
【００３２】
見かけ変位量Δｄは、変位計測窓をｍ＊（λ／２）だけ変位方向にシフトすることで計測できる。つまり、図１（ａ）において、変位前の変位計測窓１０６は、点Ｔと点Ｓの範囲である。この点Ｔと点Ｓの範囲の変位計測窓１０６内では、変位後の第１受信信号のゼロクロス点は点Ｖとして観測される。すでに述べたように点Ｖは、変位量を示す真のゼロクロス点ではない。そこで、点Ｔと点Ｓの範囲の変位計測窓１０６をｍ＊（λ／２）だけ変位方向にシフトする。
【００３３】
図１において、変位方向は時間軸で進み方向（右方向）であり、ｍ＝１であるから、変位計測窓１０６をλ／２だけ時間軸に沿い進み方向にシフトする。シフト後の変位計測窓１０８内で、変位後の第１受信信号のゼロクロス点は点Ｕとして観測される。すなわち、変位計測窓をｍ＊（λ／２）だけ変位方向にシフトすることで、シフト後の変位計測窓１０８内で、真のゼロクロス点Ｕを検出できる。
【００３４】
検出された真のゼロクロス点Ｕを用いて、見かけ変位量Δｄを算出するには、変位前の第１受信信号のゼロクロス点Ｒの位置も変位計測窓のシフトとともにシフトすればよい。図１（ａ）の場合、ゼロクロス点Ｒの位置を、変位方向にλ／２だけシフトすると、点Ｓの位置になる。そこで、図７の説明に従い、点Ｓと点Ｕの間の位相差から見かけ変位量Δｄを求めることができる。
【００３５】
したがって、対象組織の変位量は、折り返しで生じた変位量ｍ＊（λ／２）に、それを超えた±λ／２の範囲の見かけ変位量Δｄを加算して求めることができる。図１の場合ではｍ＝１であるから、対象組織の変位量は、（λ／２）＋Δｄとして求められる。このように、変位量の大きい場合、周波数の高い第１受信信号で位相差の折り返しが起こっても、周波数の低い第２受信信号の参照変位量Ｄを用い、その折り返し回数ｍを算出して補うことで、生体組織の変位量を精度よく測定できる。
【００３６】
なお、点Ｖは変位量を示す真のゼロクロス点ではないが、点Ｖと点Ｒの位相差は、（λ／２）−Δｄの値に相当するので、点Ｔと点Ｓの範囲の変位計測窓をシフトせずに、その範囲で検出されるゼロクロス点Ｖを用いて（λ／２）−Δｄを求め、この値から−Δｄを演算することもできる。この場合、Δｄに付される符号が正負逆になる。
【００３７】
また、第１受信信号の変位前のゼロクロス点Ｒを変位計測窓のシフトとともにシフトする他に、変位計測窓の基準位置、例えば時間軸で最も遅れ側の点（図１では変位計測窓１０６，１０８のそれぞれ左端）を基準位置と定め、その基準位置からのゼロクロス点の位置の変化からΔｄを求めることもできる。すなわち、変位前の変位計測窓１０６の左端から変位前のゼロクロス点Ｒの位置までの量と、変位後の変位計測窓１０８の左端と変位後の真のゼロクロス点Ｕ点の位置までの量を比較して、Δｄを求めることもできる。
【００３８】
図２は、超音波診断装置２０のブロック図である。図２において、探触子２２は、超音波の送波およびエコーの受波を行う超音波探触子で、例えば、被検者の体表面上に当接して用いられる。探触子２２内には図示されていないアレイ振動子が設けられる。アレイ振動子は複数の振動素子からなり、アレイ振動子により超音波ビームが形成される。この超音波ビームは電子走査され、その電子走査方式としては、例えば電子リニア走査や電子セクタ走査を用いることができる。
【００３９】
送受信部２４は、送信回路２６と受信回路２８とを備え、探触子２２と接続される。送受信部は、後述する送受信制御回路３２の制御の下で、超音波の送受信により、断層画像用の超音波ビームと、これとは別に変位量計測用の２つの異なる周波数の超音波ビームを形成し、受信信号を出力する回路である。より詳しくは、送信回路２６は、アレイ振動子の各チャネルごとに遅延された送信信号を供給する回路である。受信回路２８は、探触子２２からのエコー信号を増幅し、各チャネル間の受信信号の位相差を調整する整相加算等の処理を行い、受信信号としてＢモード信号処理部３６と変位量演算部３８に出力する回路である。
【００４０】
送受信制御回路３２の制御の下で、探触子２２から断層画像用の超音波ビームが生体内に走査される走査面１０の様子を図２中に示す。後述するように、断層画像用の超音波のエコー信号は信号処理および表示処理がなされ、走査面１０に断層画像として表示される。ユーザは、この断層画像を見ながら、対象組織１１の変位量計測個所、例えば血管の前壁点Ａにつき、点Ａを通過する方位を設定できる。その設定に基づき送受信制御回路３２の制御の下で、探触子２２からその方位に第１周波数の超音波１２の送受信と、第１周波数より低い第２周波数の超音波１４の送受信が交互に行われる。例えば第１周波数を８ＭＨｚ、第２周波数を４ＭＨｚとすることができる。
【００４１】
図３は、送受信制御回路３２により制御される第１周波数の超音波と第２周波数の超音波の送信タイミングチャートである。図３は横軸に時間をとり、上段に同期クロックを、中段に第１周波数の超音波の送信タイミングを、下段に第２周波数の超音波の送信タイミングを示した。このように、第１周波数の超音波の送信と第２周波数の超音波の送信は、一定の送信間隔で交互に行われる。
【００４２】
再び図２に戻り、Ｂモード信号処理部３６は、断層画像用の超音波について受信回路２８から出力される整相加算後の受信信号に基づき、エコー信号の包絡振幅を抽出する検波、包絡振幅信号の対数圧縮等の処理を行う回路である。信号処理の結果は、表示処理部４２に出力される。
【００４３】
変位量演算部３８は、変位量計測用の超音波について受信回路２８から受取った２つの異なる周波数の超音波の受信信号から、対象組織１１の変位量を演算する機能を有する。具体的には、後述するユーザにより設定されたトラッキングゲートによって、対象組織１１である血管前壁の点Ａの位置をトラッキングし、点Ａの変位量を演算する機能を有する。トラッキング機能は、後述する制御部３０の制御の下において、設定されたトラッキングゲート内で、位相を検出し、血管前壁の点Ａの変位を逐次検出する機能である。変位量を演算する機能のさらに詳細については後に詳述する。変位量演算の結果は表示処理部４２に出力される。
【００４４】
トラッキングゲートの位相検出に基づく変位計測窓の範囲は、超音波の波長をλとして±λ／２なので、この範囲を超えるときは、位相差の折り返しの問題が起こる。したがって、制御部３０の制御の下で、図１で説明したと同様に、生じた位相差の折り返し回数をｍとして、ｍ＊（λ／２）だけ変位方向に変位計測窓をシフトさせることで、±λ／２の範囲を超える変位量の場合でも誤動作することなくトラッキング機能を行うことができる。
【００４５】
表示処理部４２は、Ｂモード信号処理部３６の出力からＢモード断層画像を形成し、変位量演算部３８の演算結果を示す数値やグラフの表示に必要な処理を行う回路である。変位量に関するグラフとしては、例えば変位量の時間変化を示すグラフ等である。また、これらの画像を合成する機能も有している回路である。表示処理部４２は、いわゆるディジタルスキャンコンバータ（ＤＳＣ）や各種の画像処理回路によって構成することができる。形成されたＢモード断層画像と変位量に関する画像は、表示器４４に出力される。
【００４６】
表示器４４は、表示処理部４２により形成されたＢモード断層画像や変位量の演算結果を示す数値やグラフを表示する。また、例えば血管前壁の変位量の他に、血管後壁の変位量を計測し、その結果をあわせて表示することもできる。さらに、血管前壁の変位量と血管後壁の変位量から血管の直径を演算し、その結果を表示することもできる。
【００４７】
入力部４６は、後述の制御部３０と接続され、ユーザにより、超音波診断装置２０の動作モードの選択や必要なデータの入力が行われるスイッチ、キーボード、トラックボール等である。例えば、ユーザは、入力部４６を用いて、動作モードとして変位計測を選択できる。また、ユーザは、表示器４４の画面を見ながらトラックボール等を操作して断層画像上で変位量計測を行う点Ａを通過する所定の方位を設定でき、点Ａの位置に変位量計測のためのトラッキングゲートを設定できる。
【００４８】
制御部３０は、超音波診断装置２０の各構成と接続され、それぞれの制御を行うコントローラである。制御部３０は、上記の送受信部２４を制御する送受信制御回路３２を備える。
【００４９】
制御部３０は、変位量計測のための超音波ビームの方位を設定する制御を行う機能を有する。すなわち、入力部４６から入力された所定の方位に、第１周波数の超音波と第１周波数より低い第２周波数の超音波の方位を設定する制御を行う機能を有する。また、制御部３０は、変位量計測のためのトラッキングを制御する機能を有する。すなわち、入力部４６からなされたトラッキングゲートの設定に基づき、血管前壁の点Ａの変位を逐次検出するトラッキングを制御する機能を有する。
【００５０】
図４は、変位量演算部３８の詳細ブロック図である。変位量演算部３８は、受信回路２８が繰り返し出力する複数の第１受信信号５０を記憶する第１受信信号メモリ５４、複数の第２受信信号５２を記憶する第２受信信号メモリ５６、および変位量演算回路５８を備え、変位前後の２つの第１受信信号と変位前後の２つの第２受信信号の計４個の受信信号に基づき、対象組織の変位量６０を演算し出力する機能を有する。
【００５１】
第１受信信号メモリ５４は、受信回路２８から受取った複数の第１受信信号５０を記憶する半導体メモリ等の記憶装置である。同様に第２受信信号メモリ５６は、受信回路２８から受取った複数の第２受信信号５２を記憶する記憶装置である。各受信信号の受信順序の対応付けは、例えばメモリのアドレス等に対応付けて管理される。
【００５２】
変位量演算回路５８の中の参照変位量演算器６２は、第２受信信号メモリ５６から変位前後の２つの第２受信信号を読み出し、図７で説明した方法に従い、その信号相互間の位相差である第２位相差を求め、その第２位相差に基づき参照変位量Ｄを演算する機能を有する。
【００５３】
変位量判定器６４は、演算された参照変位量Ｄに基づき、図１で説明した方法に従い、位相差の折り返しの有無とその回数を判定する機能を有する。すなわち、参照変位量Ｄを第１周波数の超音波の波長をλとしてλ／２で除したときに得られる整数部分ｍが、位相差の折り返しの回数である。判定式の形で示すと、次式を満たす整数ｍが位相差の折り返し回数である。ｍ＝０のときは位相差の折り返しが生じていない。
【数１】
ｍ＊（λ／２）≦Ｄ≦（ｍ＋１）＊（λ／２） ・・・・（１）
【００５４】
見かけ変位量演算器６６は、図１で説明した方法に従い、判定された位相差の折り返し回数ｍに従って、第１受信信号の変位計測窓をシフトし、見かけ変位量Δｄを演算する機能を有する。図１で説明した方法を再言すれば、まず変位前後の２つの第１受信信号を第１受信信号メモリ５４から読出した後、第１受信信号の変位計測窓を、折り返し変位量ｍ＊（λ／２）だけ変位方向にシフトし、その状態で変位後の第１受信信号のゼロクロス点Ｕの位置を求める。そして、変位前の第１受信信号のゼロクロス点Ｒの位置を、折り返し変位量ｍ＊（λ／２）だけ変位方向にシフトした点Ｓと点Ｕとの相互間の位相差から見かけ変位量Δｄが算出される。
【００５５】
対象組織変位量演算器６８は、位相差の折り返し回数ｍに従って、対象組織の変位量を演算する機能を有する。すなわち、ｍが０以外の整数のときは、見かけ変位量Δｄと、折り返し変位量ｍ＊（λ／２）とを加算し、ｍ＊（λ／２）＋Δｄを演算し、その値を対象組織の変位量６０として表示処理部４２に出力する機能を有する。また、ｍが０のときは、変位前後の２つの第１受信信号を第１受信信号メモリ５４から読出し、図７で説明した方法に従い、変位前後の第１受信信号の位相差から求めた変位量をそのまま対象組織の変位量として表示処理部４２に出力する機能を有する。
【００５６】
図５は、対象組織の変位量演算の手順をフローチャートで示したものである。
【００５７】
まず受信した第１受信信号を第１受信信号メモリに順次記憶する（Ｓ１０）。同様に受信した第２受信信号を第２受信信号メモリに順次記憶する（Ｓ１２）。
【００５８】
次に第２受信信号メモリから変位前後の２つの第２受信信号を読出し、参照変位量演算器において、その相互間の位相差から参照変位量Ｄを演算する（Ｓ１４）。
【００５９】
そして、変位量判定器で、参照変位量Ｄを用いて、第１周波数の位相差の折り返しの有無と、その回数ｍを判定する（Ｓ１６）。折り返しの回数ｍは、式（１）が成立するｍを算出することで求めることができる。
【００６０】
次に、見かけ変位量演算器は、第１受信信号メモリから、変位前後の２つの第１受信信号を読出す（Ｓ１８）。そして、ｍが０以外の整数のときは、図１で説明した方法に従い、変位前の第１受信信号の変位計測窓をｍ＊（λ／２）だけ変位方向にシフトさせた上で、変位前後の第１受信信号の相互間の位相差を求め、その位相差における見かけ変位量Δｄを演算する（Ｓ２０）。
【００６１】
そして、ｍ＊（λ／２）＋Δｄを演算する（Ｓ２２）。この値を、対象組織の変位量とすることができる。
【００６２】
上記の加算、乗算等の四則演算処理は、演算処理装置を用いて行うことができる。また、演算処理装置に代わり、数式化したメモリを用い、必要なパラメータ、例えば折り返し回数ｍ、λ、参照変位量Ｄ等をメモリに入力して数式を処理し、同じ結果を得ることもできる。
【００６３】
【発明の効果】
本発明に係る超音波診断装置によれば、対象組織の診断深さを変えることなく、変位量の大きい対象組織の変位を精度よく測定することができる。
【図面の簡単な説明】
【図１】 本発明の原理説明図である。
【図２】 本発明に係る実施の形態の超音波診断装置のブロック図である。
【図３】 第１周波数の超音波と第２周波数の超音波の送信タイミングチャートである。
【図４】 変位量演算部の詳細ブロック図である。
【図５】 対象組織の変位量演算のフローチャートである。
【図６】 位相差から変位量を求める原理の説明図である。
【図７】 ２つの受信信号相互間の位相差の関係を示す図である。
【図８】 位相差に関する折り返しの問題を説明する図である。
【符号の説明】
１０ 走査面、１１ 対象組織、１２ 第１周波数の超音波、１４ 第２周波数の超音波、２０ 超音波診断装置、２２ 探触子、２４ 送受信部、２６ 送信回路、２８ 受信回路、３０ 制御部、３２ 送受信制御回路、３６ Ｂモード信号処理部、３８ 変位量演算部、５４ 第１受信信号メモリ、５６ 第２受信信号メモリ、５８ 変位量演算回路、６２ 参照変位量演算器、６４ 変位量判定器、６６ 見かけ変位量演算器、６８ 対象組織変位量演算器、１０６，１０８ 変位計測窓。
Claims (6)
- 対象組織に対し、送受信シーケンスに従って、第１周波数の超音波を繰り返し送受信して複数の第１受信信号を出力し、前記第１周波数より低い第２周波数の超音波を繰り返し送受信して複数の第２受信信号を出力する送受信手段と、
前記複数の第１受信信号についての相互間の第１位相差と、前記複数の第２受信信号についての相互間の第２位相差とに基づいて、対象組織の変位量を演算する変位量演算手段と、
を備え、
前記変位量演算手段は、
前記複数の第１受信信号の相互間の波形比較により前記第１位相差を演算する第１位相差演算部と、
前記複数の第２受信信号の相互間の波形比較により前記第２位相差を演算する第２位相差演算部と、
前記第２位相差に基づいて、前記第１位相差に関する折り返しの有無および回数を判定する判定部と、
前記判定部の判定結果に従い、前記第１位相差に基づいて前記対象組織の変位量を演算する変位量演算部と、
を備えることを特徴とする超音波診断装置。 - 請求項１に記載の超音波診断装置において、
前記送受信シーケンスは、
前記第１周波数の超音波と前記第２周波数の超音波とを交互に送信する送受信シーケンスであることを特徴とする超音波診断装置。 - 請求項１に記載の超音波診断装置において、
前記折り返し回数は、
前記第１位相差の１回の折り返し当たりの単位変位量で前記第２位相差に相当する参照変位量を除することにより算出される整数であることを特徴とする超音波診断装置。 - 請求項３に記載の超音波診断装置において、
前記変位量演算部は、
前記第１位相差から演算される見かけ変位量と、前記折り返し回数に前記単位変位量を乗じて得られる折り返し変位量とに基づいて前記対象組織の変位量を演算することを特徴とする超音波診断装置。 - 請求項４に記載の超音波診断装置において、
前記見かけ変位量は、
前記位相差の折り返しが生じない変位量観測範囲である変位計測窓を、変位方向に前記折り返し変位量分だけシフトさせて演算されることを特徴とする超音波診断装置。 - 請求項５に記載の超音波診断装置において、
前記見かけ変位量は、
前記変位計測窓のシフトにあわせ、上記変位計測窓内における変位前の第１受信信号の変位量測定のための特徴点の位置を、変位方向に前記折り返し変位量分だけシフトさせ、上記シフトさせた変位計測窓内における変位前後の前記特徴点の位置の差から演算されることを特徴とする超音波診断装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002134888A JP4113377B2 (ja) | 2002-05-10 | 2002-05-10 | 超音波診断装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002134888A JP4113377B2 (ja) | 2002-05-10 | 2002-05-10 | 超音波診断装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2003325520A JP2003325520A (ja) | 2003-11-18 |
| JP4113377B2 true JP4113377B2 (ja) | 2008-07-09 |
Family
ID=29697356
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002134888A Expired - Fee Related JP4113377B2 (ja) | 2002-05-10 | 2002-05-10 | 超音波診断装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4113377B2 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11103218B2 (en) | 2015-11-02 | 2021-08-31 | Fujifilm Corporation | Ultrasound diagnostic apparatus and control method of ultrasound diagnostic apparatus |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4810876B2 (ja) * | 2005-04-27 | 2011-11-09 | トヨタ自動車株式会社 | ボルト軸力測定装置 |
| CN113124948B (zh) * | 2021-05-20 | 2022-08-30 | 中国计量大学 | 一种基于fpga与互相关法的高精度时差测量方法 |
-
2002
- 2002-05-10 JP JP2002134888A patent/JP4113377B2/ja not_active Expired - Fee Related
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11103218B2 (en) | 2015-11-02 | 2021-08-31 | Fujifilm Corporation | Ultrasound diagnostic apparatus and control method of ultrasound diagnostic apparatus |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2003325520A (ja) | 2003-11-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP1176910B1 (en) | Method and apparatus for automatic vessel tracking in ultrasound imaging | |
| US9895138B2 (en) | Ultrasonic diagnostic apparatus | |
| US8469887B2 (en) | Method and apparatus for flow parameter imaging | |
| US5535747A (en) | Ultrasonic equipment | |
| EP1967866A1 (en) | Ultrasonic apparatus and ultrasonic diagnostic method | |
| JP3381968B2 (ja) | 超音波診断装置 | |
| JP2001128976A (ja) | Ｐｒｆ調節方法および装置並びに超音波撮影装置 | |
| WO2019187649A1 (ja) | 超音波診断装置および超音波診断装置の制御方法 | |
| JP2023155494A (ja) | 超音波診断装置 | |
| US20070073152A1 (en) | Systems and methods for acquiring images simultaneously | |
| US20230346334A1 (en) | Ultrasound diagnostic apparatus and control method of ultrasound diagnostic apparatus | |
| JP4113377B2 (ja) | 超音波診断装置 | |
| US7371219B2 (en) | Ultrasound diagnosis apparatus operable in doppler mode | |
| JP5455567B2 (ja) | 超音波診断装置 | |
| JP3688562B2 (ja) | 超音波診断装置 | |
| JP6494784B2 (ja) | 超音波診断装置および超音波診断装置の制御方法 | |
| JP2002224063A (ja) | 脈波伝播速度計測装置及び超音波診断装置 | |
| JP6778649B2 (ja) | 超音波撮像装置及び血流速度算出方法 | |
| JPH09131346A (ja) | ２次元ドプラ超音波診断装置 | |
| JP7334494B2 (ja) | 超音波診断装置、超音波診断装置の制御方法、及び、超音波診断装置の制御プログラム | |
| US20220022849A1 (en) | Ultrasound diagnostic apparatus, control method of ultrasound diagnostic apparatus, and processor for ultrasound diagnostic apparatus | |
| JP2001212140A (ja) | 超音波診断装置 | |
| JP2801450B2 (ja) | 超音波組織変位計測装置 | |
| JP2803308B2 (ja) | 超音波診断装置 | |
| JPH05317312A (ja) | 超音波診断装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050223 |
|
| RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20050223 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20070710 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070724 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070914 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20080401 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080411 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110418 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130418 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140418 Year of fee payment: 6 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |