JP5654198B2 - 超音波診断装置 - Google Patents

超音波診断装置 Download PDF

Info

Publication number
JP5654198B2
JP5654198B2 JP2008156770A JP2008156770A JP5654198B2 JP 5654198 B2 JP5654198 B2 JP 5654198B2 JP 2008156770 A JP2008156770 A JP 2008156770A JP 2008156770 A JP2008156770 A JP 2008156770A JP 5654198 B2 JP5654198 B2 JP 5654198B2
Authority
JP
Japan
Prior art keywords
signal
ultrasonic diagnostic
diagnostic apparatus
target position
continuous wave
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
Application number
JP2008156770A
Other languages
English (en)
Other versions
JP2009297351A (ja
Inventor
国田 正徳
正徳 国田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Hitachi Aloka Medical Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hitachi Aloka Medical Ltd filed Critical Hitachi Aloka Medical Ltd
Priority to JP2008156770A priority Critical patent/JP5654198B2/ja
Priority to AT09007364T priority patent/ATE536806T1/de
Priority to EP09007364A priority patent/EP2135556B1/en
Priority to US12/480,874 priority patent/US8343055B2/en
Priority to CN200910147985XA priority patent/CN101606851B/zh
Publication of JP2009297351A publication Critical patent/JP2009297351A/ja
Application granted granted Critical
Publication of JP5654198B2 publication Critical patent/JP5654198B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023Details of receivers
    • G01S7/5203Details of receivers for non-pulse systems, e.g. CW systems
    • G01S7/52031Extracting wanted echo signals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8959Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using coded signals for correlation purposes
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06Measuring blood flow
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/488Diagnostic techniques involving Doppler signals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/02Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
    • G01S15/50Systems of measurement, based on relative movement of the target
    • G01S15/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • G01S15/586Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves and based upon the Doppler effect resulting from movement of targets

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Description

本発明は、超音波診断装置に関し、特に、変調された連続波を利用する超音波診断装置に関する。
超音波診断装置の連続波を利用した技術として、連続波ドプラが知られている。連続波ドプラでは、例えば、数MHzの正弦波として構成される送信波が生体内へ連続的に放射され、生体内からの反射波が連続的に受波される。反射波には、生体内における運動体(例えば血流)によるドプラシフト情報が含まれる。そこで、そのドプラシフト情報を抽出して周波数解析することにより、運動体の速度情報を反映させたドプラ波形などを形成することができる。
連続波を利用した連続波ドプラは、パルス波を利用したパルスドプラに比べて一般に高速の速度計測の面で優れている。こうした事情などから、本願発明者は、連続波ドプラに関する研究を重ねてきた。その成果の一つとして、特許文献1において、周波数変調処理を施した連続波ドプラ(FMCWドプラ)に関する技術を提案している。
一方、連続波ドプラでは、連続波を利用していることにより位置計測が困難である。例えば、従来の一般的な連続波ドプラの装置(FMCWドプラを利用しない装置)では、位置計測を行うことができなかった。これに対し、本願発明者は、特許文献2において、FMCWドプラにより生体内組織の速度に加えて生体内組織の位置を計測することができる技術を提案している。
特開2005−253949号公報 特開2006−14916号公報
特許文献1や特許文献2に記載されたFMCWドプラの技術は、それまでにない超音波診断の可能性を秘めた画期的な技術である。本願発明者は、この画期的な技術の改良についてさらに研究を重ねてきた。
本発明は、このような背景において成されたものであり、その目的は、連続波を利用して目標位置からの生体内情報を抽出する改良技術を提供することにある。
上記目的を達成するために、本発明の好適な態様である超音波診断装置は、周期的な信号列に基づいてデジタル変調処理された連続波の送信信号を出力する送信信号処理部と、送信信号に対応した超音波を生体に送波して生体からの反射波を受波することにより受信信号を得る送受波部と、実質的に送信信号に等しい波形の参照信号を用いて受信信号に対して復調処理を施すことにより復調信号を得る受信信号処理部と、復調信号から生体内情報を抽出する生体内情報抽出部と、を有し、生体内の目標位置から得られる受信信号の周期的な信号列と参照信号の周期的な信号列との間の相関関係を調整して復調処理を施すことにより目標位置からの生体内情報を選択的に抽出する、ことを特徴とする。
上記態様では、生体内の目標位置から得られる受信信号の周期的な信号列と参照信号の周期的な信号列との間の相関関係を調整しているため、例えば、参照信号との相関が比較的大きい信号成分として目標位置からの受信信号を抽出することができる。さらに、その受信信号から例えばバンドパスフィルタあるいはローパスフィルタなどを用いてドプラ情報などを抽出することにより、目標位置からの生体内情報を選択的に抽出することが可能になる。なお、上記態様において、デジタル変調処理として、例えば、周波数シフトキーイングや位相シフトキーイングや振幅シフトキーイングなどを利用することができる。また、参照信号と送信信号は、完全に等しい波形であることが望ましい。但し、参照信号と送信信号は、実質的に等しい波形とみなせる程度の対応関係であってもよい。
望ましい態様において、前記超音波診断装置は、生体内の目標位置の深さに応じた遅延処理を施して受信信号と参照信号との間の遅延関係を調整することにより前記相関関係を調整する、ことを特徴とする。
望ましい態様において、前記超音波診断装置は、前記遅延処理により、目標位置から得られる受信信号の周期的な信号列と参照信号の周期的な信号列との間の相関を強める、ことを特徴とする。
望ましい態様において、前記超音波診断装置は、目標位置の深さに応じた遅延量だけ参照信号を遅延処理することにより、目標位置から得られる受信信号の信号列パターンと参照信号の信号列パターンとを互いに一致させる、ことを特徴とする。
望ましい態様において、前記送信信号処理部は、周期的な信号列に基づいた周波数シフトキーイングにより周波数を変化させて形成された連続波の送信信号を出力する、ことを特徴とする。
望ましい態様において、前記送信信号処理部は、位相を連続的に変化させつつ周波数を変化させた連続波の送信信号を出力する、ことを特徴とする。
望ましい態様において、前記送信信号処理部は、前記周期的な信号列の各信号ごとにその信号の期間T内において位相を+π/2または−π/2だけ変化させた連続波の送信信号を出力する、ことを特徴とする。
望ましい態様において、前記生体内情報抽出部は、生体内情報として、復調信号に含まれる直流信号成分に対応したドプラ信号成分を抽出する、ことを特徴とする。
望ましい態様において、前記超音波診断装置は、互いに異なる複数の遅延量に基づいて参照信号を遅延処理することにより、複数の目標位置に対応した複数の遅延参照信号を形成する遅延処理部をさらに有し、前記受信信号処理部は、複数の遅延参照信号を用いて受信信号に対して復調処理を施すことにより複数の目標位置に対応した複数の復調信号を形成し、前記生体内情報抽出部は、複数の復調信号に基づいて、生体内の深さ方向に並んだ複数の目標位置からの生体内情報を抽出する、ことを特徴とする。
本発明により、デジタル変調処理された連続波を利用して目標位置からのドプラ情報を抽出することが可能になる。
以下、本発明の好適な実施形態を説明する。
図1には、本発明に係る超音波診断装置の好適な実施形態が示されており、図1はその全体構成を示す機能ブロック図である。送信用振動子10は生体内へ送信波を連続的に送波し、また、受信用振動子12は生体内からの反射波を連続的に受波する。このように、送信および受信がそれぞれ異なる振動子で行われて、いわゆる連続波ドプラ法による送受信が実行される。本実施形態において利用される連続波は、デジタル変調された連続波であり、FSK変調器20によって形成される。
FSK変調器20は、パターン発生器24から供給される周期的な信号列に基づいた周波数シフトキーイング(FSK)により、RF波発振器22から供給されるRF波に対してデジタル変調処理を施して連続波を発生する。周波数シフトキーイング(FSK)により形成される連続波の波形については後に説明する。FSK変調器20は、デジタル変調された連続波を電力増幅器14に出力する。
電力増幅器14は、デジタル変調された連続波を電力増幅して送信用振動子10に供給する。そして、デジタル変調された連続波に対応する送信波が送信用振動子10から送波され、生体内からの反射波が連続的に受信用振動子12によって受波される。
前置増幅器16は、受信用振動子12から供給される受波信号に対して低雑音増幅等の受信処理を施し、受信RF信号を形成して受信ミキサ30へ出力する。受信ミキサ30は受信RF信号に対して直交検波を施して複素ベースバンド信号を生成する回路であり、2つのミキサ32,34で構成される。各ミキサは受信RF信号を所定の参照信号と混合する回路である。
受信ミキサ30の各ミキサに供給される参照信号は、デジタル変調された連続波(送信信号)に基づいて生成される。つまり、FSK変調器20から出力される連続波が遅延回路26Iと遅延回路26Qにおいて遅延処理され、遅延回路26Iにおいて遅延処理された連続波がミキサ32に供給され、遅延回路26Qにおいて遅延処理された連続波がミキサ34に供給される。
遅延回路26Iと遅延回路26Qは、目標位置の深さに応じた遅延量だけ連続波に遅延処理を施し、遅延された参照信号を出力する。遅延回路26Iと遅延回路26Qは、各々、例えばn段のシフトレジスタによって形成することができる。この場合、シフトレジスタのn段のタップから目標位置の深さに応じた遅延量のタップが選択され、選択されたタップから目標位置の深さに応じた参照信号(遅延処理された連続波)が出力される。
なお、遅延回路26Iと遅延回路26Qは、互いに連続波の位相をπ/2だけずらして遅延処理を行う。その結果、ミキサ32から同相信号成分(I信号成分)が出力され、ミキサ34から直交信号成分(Q信号成分)が出力される。そして、受信ミキサ30の後段に設けられるLPF(ローパスフィルタ)36,38によって、同相信号成分および直交信号成分の各々の高周波数成分がカットされて検波後の必要な帯域のみの復調信号が抽出される。
後の原理説明で詳述するが、各ミキサで実行される受信RF信号と参照信号との混合処理の結果である受信ミキサ出力信号(復調信号)には、目標位置からの受信信号成分が多く含まれている。LPF36,38において、その目標位置からの受信信号成分に含まれている直流信号成分(移動体が存在すればドプラ信号成分が含まれる)が抽出される。
FFT回路(高速フーリエ変換回路)40,42は、復調信号(同相信号成分および直交信号成分)の各々に対してFFT演算を実行する。その結果、FFT回路40,42において復調信号が周波数スペクトラムに変換される。なお、FFT回路40,42から出力される周波数スペクトラムは、回路の設定条件などにより周波数分解能δfの周波数スペクトラムデータとして出力される。
ドプラ情報解析部44は、周波数スペクトラムに変換された復調信号から、その極性を含んだドプラ情報を抽出する。その際、予め遅延回路26I,26Qによって、生体内の目標位置の深さに応じて参照信号と受信信号との間の遅延関係が調整されているため、目標位置からのドプラ情報が選択的に抽出される。遅延関係の調整と目標位置からのドプラ情報の抽出との関連については、後の原理説明において詳述する。ドプラ情報解析部44は、生体内の各深さ(各位置)ごとにドプラ情報を抽出して、例えば、超音波ビーム(音線)上の各深さごとに生体内組織の速度を算出し、リアルタイムで出力する。なお、超音波ビームを走査させて二次元的あるいは三次元的に生体内組織の各位置の速度を算出してもよい。
表示処理部46は、生体組織の深さ(位置)ごとの速度に基づいて、例えばドプラ波形や、深さ速度の情報を含むグラフなどを形成し、形成したドプラ波形やグラフなどを表示部48にリアルタイムで表示させる。なお、図1に示す超音波診断装置内の各部は、システム制御部50によって制御される。つまり、システム制御部50は、送信制御や受信制御や表示制御などを行う。
以上、概説したように、本実施形態では、デジタル変調処理された連続波に対応した超音波を送受波して受信信号を得て、生体内の目標位置の深さに応じて参照信号と受信信号との間の遅延関係を調整し、目標位置からの受信信号と参照信号との間の相関を強めて復調処理を施すことにより、目標位置からの生体内情報としてドプラ情報を選択的に抽出する。そこで、目標位置からのドプラ情報が選択的に抽出される原理について詳述する。
図2は、周波数シフトキーイング(FSK)により形成される連続波の送信信号を説明するための図である。図2(I)には、パターン発生器(図1の符号24)から出力される周期的な信号列の一例が示されている。パターン発生器は、例えば図2(I)に示すような、ランダムに値を変化させた2値符号(擬似ランダム信号)を発生する。周期的な信号列である擬似ランダム信号の符号系列としては、PN(Pseudo Noise)系列、M系列、Gorey系列など、パルス圧縮などで実用化されている符号系列を用いればよい。
図2(II)には、FSK変調器(図1の符号20)において形成される変調された連続波(送信信号)が示されている。FSK変調器は、図2(I)の2値符号に基づいた周波数シフトキーイングにより、RF波(搬送波)に対してデジタル変調処理を施して連続波の送信信号を形成する。FSK変調器は、例えば、2値符号が「1」のビット期間において周波数f1とし、2値符号が「0」のビット期間において周波数f2とすることにより、図2(II)の送信信号を形成する。
こうして、例えば図2(II)の送信信号に対応した連続波の超音波が送信用振動子(図1の符号10)から出力され、受信用振動子(図1の符号12)を介して生体内から受信信号が得られる。
図3は、FSK変調器から出力される連続波の位相の変化を説明するための図である。なお、図3では、搬送波の位相は常に一定として0(ゼロ)で表現し、周波数f1あるいは周波数f2に対応して位相回転した量を表現している。周波数シフトキーイング(FSK)によりデジタル変調処理された信号を復調処理する際には、その信号の周波数を弁別すればよい。例えば、図2(II)の連続波であれば、周波数f1と周波数f2を弁別することにより、図2(II)の信号から図2(I)の2値符号を復調することができる。単に周波数を弁別するだけであれば、符号に対応した各周波数の位相は特に決めなくてもよい。例えば、符号の切り替わりの時点で位相が不連続に変化しても、単に周波数を弁別するだけであれば問題はない。しかし、符号変化時に位相が不連続となると、その分だけ周波数帯域が広がってしまう。そこで、符号変化時に位相が連続的に変化するように各周波数と1ビットの時間長Tを決めておけば、より狭い帯域でFSK変調および復調を行うことが可能となる。
図3は、その条件を満足したときの連続波の位相変化を描いたものである。図3に示す連続波の位相変化は、2値符号の符号変化時に位相が連続的に変化し、2値符号が「1」のビット期間(周波数f1の期間)において連続波の位相が+φだけ変化し(φだけ進み)、2値符号が「0」のビット期間(周波数f2の期間)において連続波の位相が−φだけ変化する(φだけ遅れる)場合の連続波の位相変化である。
図3に示す例の連続波では、符号がどのように変化しようとも、連続波の位相が連続的に変化し、必ず図3のグラフ内の破線の増加直線または破線の減少直線上を移動するという特徴がある。さらに、符号に対応した周波数の間隔をある程度小さく設定しておけば、狭帯域のFSK連続波を得ることができる。周波数間隔が狭くなると1ビットで回転する位相も小さくなる。この位相回転量がπ/2となるように周波数間隔を設定すると、つまり1ビットの期間T内において位相を+π/2または−π/2だけ変化させた連続波とすることにより、狭帯域でありながら、符号誤り率を小さく保った状態でデジタル信号の送受信をすることができる。このときのFSK方式は特別にMSK(Minimum Shift Keying)方式と呼ばれる。
図4は、FSK変調器から出力される連続波の周波数スペクトラムを説明するための図である。図4には、横軸を周波数として縦軸に各周波数成分の電力を示したグラフが示されており、グラフ内には、周波数シフトキーイングにより周波数f0の搬送波(RF波)から周波数をf1,f2だけ変化させたFSK信号(FSKの連続波)の周波数スペクトラムが変調度により変化する様子が示されている。なお、FSKの変調度mは次式のように定義される。
Figure 0005654198
上式で表現される変調度mは、周波数差(f1とf2の差)とビットレートfb(bit/s.)との比であり、アナログFM変調の変調度βと類似の意味を持っている。変調度mを大きくするということは、周波数差を大きく設定することに相当し、その結果、周波数スペクトラムは広がってゆく。したがって、例えば超音波プローブなどの周波数帯域を考慮しながら変調度mを設定することが望ましい。
本実施形態では、FSK変調器(図1の符号20)において形成された送信信号に対して遅延処理を施して参照信号を形成し、受信ミキサ(図1の符号30)においてその参照信号を用いて受信信号に対してミキサ処理(参照信号と受信信号の乗算)が行われる。このミキサ処理において、遅延処理された参照信号の位相に対応する深さ(目標位置の深さ)からの受信信号と参照信号との間の相関が強められて最大となり、その他の深さからの受信信号と参照信号との間の相関が極端に小さくなる。つまり、本実施形態の超音波診断装置は、位置選択性の機能を備えている。
図5から図7は、本実施形態における位置選択性を説明するための図であり、図5から図7の各図には、目標位置からの受信信号(受信波)と参照信号(参照波)の各位相の時間変化の様子と、目標位置からの受信信号と参照信号との間の位相差の時間変化の様子が図示されている。図5から図7の各図において、位相φは、擬似ランダム信号の1ビットの期間Tにおける連続波の位相の変化である。そして、各図において位相差の時間変化については、φ=π/2の場合における変化も図示されている。
乗算器(図1の受信ミキサ30)は、受信信号と参照信号の両者の位相差に応じた結果(例えば電圧)を出力する。例えば、両者の位相差がπ/2,3π/2,−π/2,−3π/2の場合に乗算器の出力値が「0(ゼロ)」であり、両者の位相差が0(ゼロ),2πの場合に乗算器の出力値が「+1」であり、両者の位相差がπ,−πの場合に乗算器の出力値が「−1」となる。
なお、図5から図7の例においては、7ビットの周期で変化する擬似ランダムパターンを用いているが、装置を具現化する際には、例えば100〜200ビット程度の周期で変化する擬似ランダムパターンを用いることが望ましい。
図5は、参照信号の位相が目標位置からの受信信号と一致している場合の様子を示している。参照信号の位相が目標位置からの受信信号の位相と一致している場合には、両者の位相差は常に0(ゼロ)となる。両者の位相差が0の場合の乗算器(図1の受信ミキサ30)の出力を「+1」とすると、両者の位相が一致している場合に7ビットの期間内において得られる乗算器出力の合計値は「+7」となる。
これに対し、図6,図7は、参照信号の位相が目標位置からの受信信号と一致していない場合の様子を示している。
図6は、目標位置からの受信信号(実線)に対して、参照信号(破線)の位相がTだけ遅れている場合の様子を示している。この場合の両者の位相差は、図6の下段に示すとおりであり、φ=π/2の場合において、乗算器(図1の受信ミキサ30)の出力が常に「0(ゼロ)」となり、7ビットの期間内において得られる乗算器出力の合計値も「0(ゼロ)」となる。
また、図7は、目標位置からの受信信号(実線)に対して、参照信号(破線)の位相が2Tだけ遅れている場合の様子を示している。この場合の両者の位相差は、図7の下段に示すとおりであり、φ=π/2の場合において、乗算器(図1の受信ミキサ30)の出力は「+1」と「−1」をランダムに繰り返し、7ビットの期間内において得られる乗算器出力の合計値が「+1」となる。
このように、参照信号の位相が目標位置からの受信信号の位相と一致している場合には乗算器出力の合計値が「+7」と極端に大きくなり、両者の位相が一致していない場合には乗算器出力の合計値が「0」や「+1」と極端に小さくなる。
図8は、参照波の位相と乗算器出力との関係を示す図である。図8に示すグラフは、図5から図7の例に対応しており、図8に示すグラフの縦軸は、7ビットの期間内において得られる乗算器出力の合計値である。また、図8に示すグラフの横軸は、受信信号に対する参照波(参照信号)の位相、つまり、受信信号と参照信号の位相差を示している。なお、図8の横軸は、φ=π/2を基準とした場合の(φ=π/2を「1」とした場合の)位相の相対的な大きさを示している。
参照信号と受信信号の位相差が0(ゼロ)の場合には、図5を利用して説明したように、乗算器出力の合計値は「+7」となる。したがって、図8において、参照波の位相が0の場合における乗算器出力値が「+7」となっている。また、図6,図7を利用して説明したように、位相差がφの場合における乗算器出力値は「0」となり、位相差が2φの場合における乗算器出力値は「+1」となる。したがって、図8において、参照波の位相が1の場合における乗算器出力値は「0」となり、参照波の位相が2の場合における乗算器出力値は「+1」となっている。
図8には、参照波の位相を3,4,5,6,7,・・・とさらに変化させた場合の乗算器出力値も示している。乗算器出力値(7ビットの期間内の合計値)は、位相差が存在する場合に比べて、位相差が0(ゼロ)の場合に極端に大きな値「+7」を示していることがわかる。また、参照波の位相が相対値で7だけずれると、7ビットの周期の擬似ランダム信号が1周期分だけずれるため、乗算器出力値は、位相差が0の場合と同じく極端に大きな値「+7」を示す。
以上のことから、本実施形態においては、目標までの遅延時間に対応するように参照信号の位相を調整することにより、その目標の深さのみの反射波電力とドプラ情報を選択的に検出することができる。また、その選択性は、擬似ランダムパターンの長さが長くなるほどシャープになる。そして、その選択性は、パターン長を大きくするほど、パルス波ドプラ(PWドプラ)のレンジゲートと同様の特性に近づいてくる。つまり、連続波の長所を維持しつつ、パルス波のメリットを得ることも可能になる。
図9は、本実施形態における復調信号を説明するための図であり、図9には、ミキサ処理により得られる復調信号の周波数スペクトラムが示されている。図9の復調信号は、相関が最大の場合における参照信号と受信信号の乗算結果に相当する。つまり、目標位置からの受信信号と、目標位置の深さに位相を合わせた参照信号との間の乗算結果が、図9の復調信号となる。
図9に示す復調信号には、直流信号成分と、RF信号の周波数f0の2倍の高調波成分が含まれている。なお、図2(II)に示したように、FSKによりデジタル変調された連続波の周波数は、例えば周波数f1と周波数f2をランダムに繰り返す。周波数f1と周波数f2の差が小さく、二つの周波数に対応したスペクトラムが周波数f0のスペクトラムとして重なって検出できる程度の場合には、図9のような周波数スペクトラムとなる。
ドプラ信号は、直流信号成分と周波数f0の2倍の高調波成分に付着した形で出現する。なお、LPF(図1の符号36,38)において、高調波成分がカットされて直流信号成分のみが抽出されるため、FFT回路(図1の符号40,42)においては、図9に示す直流信号成分と周波数f0の2倍の高調波成分のうち、直流信号成分の周波数スペクトラムのみが形成される。そして、ドプラ情報解析部(図1の符号44)において、図9に示す直流信号成分の周波数スペクトラムからドプラ信号が抽出され、ドプラシフト量などに基づいて、目標位置に存在する血流の流速などが算出される。受信ミキサ(図1の符号30)において、直交検波を施しているため、流速の極性を判断することもできる。直流信号成分の周波数スペクトラムからクラッタ信号を抽出して、目標位置に存在する血管壁の位置などを算出してもよい。
なお、超音波が生体内を伝搬する際の減衰特性は、周波数依存性減衰(FDA)と呼ばれ、生体内の軟部組織では、ほぼ周波数と伝搬距離に比例することが知られている。FDAは、復調信号の周波数スペクトラムの形に影響を与える。つまり、体表付近から得られ復調信号の周波数スペクトラムに比べて、深い位置から得られる復調信号の周波数スペクトラムほど減衰が大きい。そこで、FFT回路(図1の符号40,42)から出力される周波数スペクトラムに対して、生体内における周波数依存性減衰の補償処理を施すようにしてもよい。例えば、深い位置から得られる信号ほど大きく増幅し、位置(深さ)に応じた減衰量の変化をキャンセルする。受信ミキサ(図1の符号30)の出力に対して周波数依存性減衰の補償処理を施してもよい。
図10は、パターン発生器(図1の符号24)の具体例を示す図である。図10の回路構成例では、水晶発振器(OSC)から出力される搬送波(RF信号)を矩形信号に変換してからn分周し、n分周後の信号(データ)をシフトレジスタに入力している。そして、シフトレジスタの途中から引き出したn1〜n3の出力とシフトレジスタの最終段からの出力とを利用して、排他論理和をとることにより、PNパターンを得ることができる。PNパターンの長さは、シフトレジスタがn段の場合は、2n―1ビットとなる。例えば、n=10であれば1023ビットのPNパターンが容易に構成できる。この場合、受信信号と参照信号の位相がぴったり一致したときの乗算器出力は、1ビット分の1023倍となる。しかし、他の部分の出力はせいぜい数倍なので、選択能力を著しく向上させることができる。図10の回路は、デジタル回路で構成できるので、容易にIC化できるという利点もある。
パターン発生器から出力される周期的な信号列に基づいて、FSK変調器(図1の符号20)において周波数シフトキーイング(FSK)された連続波が形成される。FSK変調器は、例えば、予め用意された周波数f1と周波数f2の搬送波を利用し、パターン発生器から出力される周期的な信号列、例えば図10の回路によって得られるPNパターンに応じて、周波数f1と周波数f2の搬送波を適宜選択し、選択した搬送波を時間軸上において互いに接続するようにして、FSKされた連続波を形成してもよい。
図11は、遅延回路(図1の符号26I,26Q)と受信ミキサ(図1の符号30)の変形例を説明するための図である。図11の回路では、n段のシフトレジスタによって参照信号が遅延処理され、シフトレジスタのn1からnnまでの各タップから、互いに遅延時間の異なる複数の遅延参照信号が同時に出力される。そして、n1からnnまでの各タップに対応したn個のミキサにおいて、受信信号と各遅延参照信号の乗算処理が行われ、複数の遅延参照信号に対応した複数の乗算結果が同時に(並列的に)出力される。生体内の深さ方向に並んだ複数の目標位置と、シフトレジスタのn1からnnまでの各タップとを対応付けておくことにより、複数の目標位置からのドプラ情報を同時に取得することができる。例えば、超音波ビーム上における全ての位置からの組織情報を同時に検出するようにしてもよい。
以上、本発明の好適な実施形態といくつかの変形形態を説明したが、上述した本発明の好適な実施形態等は、あらゆる点で単なる例示にすぎず、本発明の範囲を限定するものではない。本発明は、その本質を逸脱しない範囲で各種の変形形態を包含する。
例えば、上述した実施形態においては、連続波をデジタル変調する際に周波数シフトキーイング(FSK)を利用している。このFSKに換えて、デジタル変調方式として当業者において明らかな位相シフトキーイング(PSK)や振幅シフトキーイング(ASK)などを利用してもよい。なお、デジタル変調された連続波のデータをメモリなどに記憶しておき、このメモリから読み出されるデータに基づいて、当該連続波を生成してもよい。
本発明に係る超音波診断装置の全体構成を示す機能ブロック図である。 FSKにより形成される連続波の送信信号を説明するための図である。 FSK変調器から出力される連続波の位相変化を説明するための図である。 FSK変調器から出力される連続波の周波数スペクトラムを説明するための図である。 参照信号の位相が受信信号と一致している場合の様子を示す図である。 参照信号の位相がTだけ遅れている場合の様子を示す図である。 参照信号の位相が2Tだけ遅れている場合の様子を示す図である。 参照波の位相と乗算器出力との関係を示す図である。 ミキサ処理により得られる復調信号の周波数スペクトラムを示す図である。 パターン発生器の具体例を示す図である。 遅延回路と受信ミキサの変形例を説明するための図である。
符号の説明
20 FSK変調器、22 RF波発振器、24 パターン発生器、26I,26Q 遅延回路、40,42 FFT回路、44 ドプラ情報解析部。

Claims (7)

  1. 周期的な信号列に基づいてデジタル変調処理された連続波の送信信号を出力する送信信号処理部と、
    送信信号に対応した超音波を生体に送波して生体からの反射波を受波することにより受信信号を得る送受波部と、
    実質的に送信信号に等しい波形の参照信号を用いて受信信号に対して復調処理を施すことにより復調信号を得る受信信号処理部と、
    復調信号から生体内情報を抽出する生体内情報抽出部と、
    を有し、
    生体内の目標位置の深さに応じた遅延処理により前記目標位置から得られる受信信号の周期的な信号列と参照信号の周期的な信号列との間の相関関係を調整して復調処理を施すことにより、前記目標位置からの生体内情報を選択的に抽出し、
    前記送信信号処理部は、周期的な信号列に基づいた周波数シフトキーイングにより周波数を変化させて形成された連続波の送信信号を出力するにあたり、位相を連続的に変化させつつ周波数を変化させた連続波の送信信号を出力する、
    ことを特徴とする超音波診断装置。
  2. 請求項1に記載の超音波診断装置において、
    生体内の目標位置の深さに応じた遅延処理を施して受信信号と参照信号との間の遅延関係を調整することにより前記相関関係を調整する、
    ことを特徴とする超音波診断装置。
  3. 請求項2に記載の超音波診断装置において、
    前記遅延処理により、目標位置から得られる受信信号の周期的な信号列と参照信号の周期的な信号列との間の相関を強める、
    ことを特徴とする超音波診断装置。
  4. 請求項3に記載の超音波診断装置において、
    目標位置の深さに応じた遅延量だけ参照信号を遅延処理することにより、目標位置から得られる受信信号の信号列パターンと参照信号の信号列パターンとを互いに一致させる、
    ことを特徴とする超音波診断装置。
  5. 請求項に記載の超音波診断装置において、
    前記送信信号処理部は、前記周期的な信号列の各信号ごとにその信号の期間T内において位相を+π/2または−π/2だけ変化させた連続波の送信信号を出力する、
    ことを特徴とする超音波診断装置。
  6. 請求項1からのいずれか1項に記載の超音波診断装置において、
    前記生体内情報抽出部は、生体内情報として、復調信号に含まれる直流信号成分に対応したドプラ信号成分を抽出する、
    ことを特徴とする超音波診断装置。
  7. 請求項1からのいずれか1項に記載の超音波診断装置において、
    互いに異なる複数の遅延量に基づいて参照信号を遅延処理することにより、複数の目標位置に対応した複数の遅延参照信号を形成する遅延処理部をさらに有し、
    前記受信信号処理部は、複数の遅延参照信号を用いて受信信号に対して復調処理を施すことにより複数の目標位置に対応した複数の復調信号を形成し、
    前記生体内情報抽出部は、複数の復調信号に基づいて、生体内の深さ方向に並んだ複数の目標位置からの生体内情報を抽出する、
    ことを特徴とする超音波診断装置。
JP2008156770A 2008-06-16 2008-06-16 超音波診断装置 Expired - Fee Related JP5654198B2 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
JP2008156770A JP5654198B2 (ja) 2008-06-16 2008-06-16 超音波診断装置
AT09007364T ATE536806T1 (de) 2008-06-16 2009-06-03 Diagnostisches ultraschallgerät
EP09007364A EP2135556B1 (en) 2008-06-16 2009-06-03 Ultrasound diagnostic apparatus
US12/480,874 US8343055B2 (en) 2008-06-16 2009-06-09 Ultrasound diagnostic apparatus
CN200910147985XA CN101606851B (zh) 2008-06-16 2009-06-12 超声波诊断装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008156770A JP5654198B2 (ja) 2008-06-16 2008-06-16 超音波診断装置

Publications (2)

Publication Number Publication Date
JP2009297351A JP2009297351A (ja) 2009-12-24
JP5654198B2 true JP5654198B2 (ja) 2015-01-14

Family

ID=41059705

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008156770A Expired - Fee Related JP5654198B2 (ja) 2008-06-16 2008-06-16 超音波診断装置

Country Status (5)

Country Link
US (1) US8343055B2 (ja)
EP (1) EP2135556B1 (ja)
JP (1) JP5654198B2 (ja)
CN (1) CN101606851B (ja)
AT (1) ATE536806T1 (ja)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8100832B2 (en) * 2007-04-27 2012-01-24 Hitachi Aloka Medical, Ltd. Ultrasound diagnostic apparatus
JP5629052B2 (ja) * 2008-06-03 2014-11-19 日立アロカメディカル株式会社 超音波診断装置
JP5654198B2 (ja) 2008-06-16 2015-01-14 日立アロカメディカル株式会社 超音波診断装置
JP6116853B2 (ja) * 2011-11-30 2017-04-19 東芝メディカルシステムズ株式会社 超音波診断装置及び画像処理方法
US8737555B2 (en) * 2011-12-22 2014-05-27 Landis+Gyr Technologies, Llc Digital signal processing for PLC communications having communication frequencies
CN103549975B (zh) * 2013-10-14 2016-08-10 华南理工大学 基于跳频超声信号分离的结石检测装置及其方法
US10270491B2 (en) 2017-08-31 2019-04-23 Landis+Gyr Technologies, Llc Power-line communication systems AMD methods having location-extendable collector for end-point data
US10340980B1 (en) 2018-05-07 2019-07-02 Landis+Gyr Technologies, Llc Time synchronization apparatuses and methods for power-distribution systems and the like
US11255948B1 (en) * 2019-06-08 2022-02-22 Robert Dickerman Determining frequency spectra and polarity of frequency for imbalanced quadrature signals including for Doppler system target speed and direction
FR3101955B1 (fr) * 2019-10-09 2021-10-29 Bh Tech Dispositif et procede de mesure de distance par ultrasons

Family Cites Families (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4167879A (en) * 1978-06-06 1979-09-18 Panametrics, Inc. Method and apparatus for examining a solid
US4176351A (en) * 1978-08-18 1979-11-27 Raytheon Company Method of operating a continuous wave radar
FR2447041A1 (fr) * 1979-01-19 1980-08-14 Inst Nat Sante Rech Med Perfectionnements aux velocimetres doppler a bruit pseudo-aleatoires
US4578677A (en) * 1983-09-23 1986-03-25 The United States Of America As Represented By The Secretary Of The Navy Range doppler coupling magnifier
CN86203861U (zh) 1986-05-31 1987-11-11 西安交通大学 一种扩展连续波多普勒血流仪功能的装置
US5022400A (en) * 1988-02-08 1991-06-11 Walters Glenn A Large time bandwidth echographic signal processor
US4918706A (en) * 1988-12-28 1990-04-17 Sperry Marine Inc. Spread spectrum long loop receiver
US5224482A (en) * 1991-04-08 1993-07-06 Kabushiki Kaisha Toshiba Ultrasound high velocity flow correlation measurement using coded pulses
JPH04351093A (ja) * 1991-05-28 1992-12-04 Mitsubishi Heavy Ind Ltd 超音波テレメーター
JPH08204766A (ja) * 1995-01-25 1996-08-09 Katsuyoshi Azeyanagi 時間制限波形の周波数分析検出方式
US6179781B1 (en) * 1999-03-31 2001-01-30 Acuson Corporation Medical diagnostic ultrasound method and apparatus for improving doppler processing
US6491631B2 (en) * 2001-01-11 2002-12-10 General Electric Company Harmonic golay-coded excitation with differential pulsing for diagnostic ultrasound imaging
JP2002136522A (ja) * 2000-11-02 2002-05-14 Japan Science & Technology Corp 超音波測定装置
US7094204B2 (en) * 2002-08-23 2006-08-22 Siemens Medical Solutions Usa, Inc. Coded excitation imaging for use with bipolar, unipolar and other waveforms
US6953434B2 (en) * 2002-09-24 2005-10-11 Ge Medical Systems Global Technology Company, Llc Method and apparatus to enhance ultrasound contrast imaging using stepped-chirp waveforms
US6960169B2 (en) * 2003-05-19 2005-11-01 Siemens Medical Solutions Usa, Inc. Spread spectrum coding for ultrasound contrast agent imaging
JP4688466B2 (ja) * 2003-10-10 2011-05-25 アロカ株式会社 超音波診断装置
WO2005058167A2 (ja) * 2003-12-02 2005-06-30 Hitachi Medical Corp 超音波診断装置
JP4490128B2 (ja) * 2004-01-29 2010-06-23 アロカ株式会社 超音波診断装置
JP4652731B2 (ja) 2004-07-01 2011-03-16 アロカ株式会社 超音波診断装置
US7338450B2 (en) * 2004-08-27 2008-03-04 General Electric Company Method and apparatus for performing CW doppler ultrasound utilizing a 2D matrix array
EP1769747A4 (en) * 2004-10-20 2008-04-09 Toshiba Kk ULTRASOUND DOPPLER DIAGNOSTIC DEVICE
JP2006288974A (ja) * 2005-04-14 2006-10-26 Kanazawa Inst Of Technology 超音波診断装置
US7887487B2 (en) * 2005-07-11 2011-02-15 Siemens Medical Solutions Usa, Inc. Ultrasound diagnostic flow imaging with coded excitation
JP2007330541A (ja) 2006-06-15 2007-12-27 Kanazawa Inst Of Technology 超音波診断装置
US8100832B2 (en) * 2007-04-27 2012-01-24 Hitachi Aloka Medical, Ltd. Ultrasound diagnostic apparatus
JP5654198B2 (ja) 2008-06-16 2015-01-14 日立アロカメディカル株式会社 超音波診断装置

Also Published As

Publication number Publication date
ATE536806T1 (de) 2011-12-15
JP2009297351A (ja) 2009-12-24
CN101606851A (zh) 2009-12-23
EP2135556B1 (en) 2011-12-14
US20090312636A1 (en) 2009-12-17
CN101606851B (zh) 2012-10-03
US8343055B2 (en) 2013-01-01
EP2135556A1 (en) 2009-12-23

Similar Documents

Publication Publication Date Title
JP5629052B2 (ja) 超音波診断装置
JP5654198B2 (ja) 超音波診断装置
JP5459963B2 (ja) 超音波診断装置
JP5308748B2 (ja) 超音波診断装置
JP2010279531A (ja) 超音波診断装置
JP5325503B2 (ja) 超音波診断装置
JP2010110503A (ja) 超音波診断装置
JP2009297444A (ja) 超音波診断装置
JP2011217898A (ja) 超音波診断装置
JP2010269004A (ja) 超音波診断装置
JP2009261749A (ja) 超音波診断装置
JP2011036599A (ja) 超音波診断装置
JP2011240006A (ja) 超音波診断装置
JP2011104241A (ja) 超音波診断装置
JP2010194109A (ja) 超音波診断装置
JP2010075289A (ja) 超音波診断装置
JP2011011041A (ja) 超音波診断装置
JP5202138B2 (ja) 超音波診断装置
JP2011217897A (ja) 超音波診断装置
JP2010162082A (ja) 超音波診断装置
JP2011217875A (ja) 超音波診断装置
JP2009254526A (ja) 超音波診断装置
JP2011036506A (ja) 超音波診断装置
JP2011062386A (ja) 超音波診断装置
JP2009247365A (ja) 超音波診断装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20110329

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20121228

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130108

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130301

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20130618

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130917

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20130925

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20141020

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20141120

R150 Certificate of patent or registration of utility model

Ref document number: 5654198

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

LAPS Cancellation because of no payment of annual fees