JP2014217605A - Video system and image forming method - Google Patents
Video system and image forming method Download PDFInfo
- Publication number
- JP2014217605A JP2014217605A JP2013099248A JP2013099248A JP2014217605A JP 2014217605 A JP2014217605 A JP 2014217605A JP 2013099248 A JP2013099248 A JP 2013099248A JP 2013099248 A JP2013099248 A JP 2013099248A JP 2014217605 A JP2014217605 A JP 2014217605A
- Authority
- JP
- Japan
- Prior art keywords
- frequency
- wave
- component data
- data
- frequency component
- 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.)
- Granted
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
Abstract
Description
本発明は、映像装置及び映像形成方法に関する。 The present invention relates to a video apparatus and a video forming method.
超音波診断装置は、超音波プローブに内蔵された超音波トランスデューサから発生した超音波を被検体内に放射し、被検体組織の音響インピーダンスの差異によって生ずる反射信号を超音波トランスデューサによって受信する。そして、この信号をフーリエ変換等の演算処理を施して映像を形成し、表示部に表示する。
この超音波診断装置は、超音波プローブを被検体の体表に接触させるだけの、非侵襲且つ簡単な操作で、リアルタイムにて二次元映像が容易に観察できる。このため、心臓等の機能検査や各種臓器の形態診断に広く用いられている。特に、超音波診断は、被検体に与える照射障害が殆ど無いことから、胎児の発育診断には不可欠な映像診断法である。
The ultrasonic diagnostic apparatus radiates an ultrasonic wave generated from an ultrasonic transducer built in an ultrasonic probe into a subject, and receives a reflected signal generated by a difference in acoustic impedance of the subject tissue by the ultrasonic transducer. The signal is subjected to arithmetic processing such as Fourier transform to form an image and displayed on the display unit.
This ultrasonic diagnostic apparatus can easily observe a two-dimensional image in real time by a non-invasive and simple operation in which an ultrasonic probe is brought into contact with the body surface of a subject. For this reason, it is widely used for functional tests such as the heart and morphological diagnosis of various organs. In particular, the ultrasonic diagnosis is an image diagnostic method indispensable for the fetal growth diagnosis because there is almost no irradiation obstacle to the subject.
なお、本発明の技術分野に近い先行技術文献を特許文献１及び非特許文献１に示す。
特許文献１には、被検体ボリューム中の超音波散乱媒体の三次元イメージング方法およびシステムの技術内容が開示されている。
非特許文献１には、外部より加えられた低周波振動の位相変化から伝搬速度を推定する、加振映像法の技術内容が開示されている。
Prior art documents close to the technical field of the present invention are shown in Patent Document 1 and Non-Patent Document 1.
Patent Document 1 discloses the technical content of a three-dimensional imaging method and system for an ultrasonic scattering medium in a subject volume.
Non-Patent Document 1 discloses the technical content of an excitation video method that estimates the propagation velocity from the phase change of low-frequency vibration applied from the outside.
ある対象物を映像化するイメージング技術には、走査型電子顕微鏡等の電子線等による近接的な手法と、超音波、電波、光等の波動を用いる遠隔的な手法に大別できる。超音波診断装置は後者の内、コヒーレントな波動を用いたイメージング手法に属する。
この、コヒーレントな波動を用いる遠隔的なイメージング技術では通常コヒーレント波を送受信し、送受信機をスキャンすることによって、ビーム方向（深さ方向）とそれに直交する方向（横方向）の反射係数の映像化が可能である。この映像化において、深さ方向の空間分解能は送受信機の帯域幅、横方向の空間分解能は送受信機の点拡がり関数のサイズによって決定される。つまり、超音波診断装置が出力する映像の横方向の分解能は、超音波プローブ内の超音波トランスデューサの開口幅と出力される超音波の波長に依存する。このため、どんなに装置を高性能化しても、センサの点広がり関数の空間分解能以上に詳細な映像を得ることはできない。そして、超音波トランスデューサの設計上、超音波の波長を短くする技術にも限界がある。
Imaging techniques for imaging a certain object can be broadly divided into a proximity technique using an electron beam such as a scanning electron microscope and a remote technique using waves such as ultrasonic waves, radio waves, and light. Among the latter, the ultrasonic diagnostic apparatus belongs to an imaging technique using coherent waves.
In this remote imaging technology using coherent waves, coherent waves are usually transmitted and received, and the reflection coefficient in the beam direction (depth direction) and the direction orthogonal to it (lateral direction) are visualized by scanning the transceiver. Is possible. In this visualization, the spatial resolution in the depth direction is determined by the bandwidth of the transceiver, and the spatial resolution in the lateral direction is determined by the size of the point spread function of the transceiver. That is, the resolution in the horizontal direction of the image output by the ultrasonic diagnostic apparatus depends on the aperture width of the ultrasonic transducer in the ultrasonic probe and the wavelength of the output ultrasonic wave. For this reason, no matter how high the performance of the device, no more detailed images can be obtained than the spatial resolution of the sensor point spread function. Further, there is a limit to the technology for shortening the wavelength of the ultrasonic wave due to the design of the ultrasonic transducer.
本発明は係る状況に鑑みてなされたものであり、探査用キャリア信号の周波数を操作せずに、極めて簡便な技術要素を追加するだけで、容易に横方向（あるいはセンサのスキャン方向）の空間分解能を向上できる、映像装置及び映像形成方法を提供することを目的とする。 The present invention has been made in view of such a situation, and it is easy to add a very simple technical element without manipulating the frequency of the carrier signal for exploration. An object of the present invention is to provide a video apparatus and a video forming method capable of improving the resolution.
上記課題を解決するために、本発明の映像装置は、第一の周波数の波動を測定対象に送信し、測定対象から反射した波動を受信するトランスデューサと、測定対象に第二の周波数の振動を加える加振器と、トランスデューサが出力する信号を検波する直交検波器とを具備する。更に、直交検波器が出力する信号に基づくデータにフーリエ変換を施して、第一の周波数の成分である基本波成分データと、加振器によって第一の周波数から周波数が増加した増加周波数成分データとを具備する。更に、加振器によって第一の周波数から周波数が減少した減少周波数成分データとを出力する離散フーリエ変換器と、増加周波数成分データ及び減少周波数成分データに補正演算を行う補償演算部とを具備する。更に、補償演算部によって補正された増加周波数成分データ及び減少周波数成分データを基本波成分データに加算する加算器とを具備する。 In order to solve the above problems, an image device of the present invention transmits a wave having a first frequency to a measurement object, receives a wave reflected from the measurement object, and vibrates a second frequency in the measurement object. It includes a vibrator to be added and a quadrature detector for detecting a signal output from the transducer. Furthermore, the fundamental wave component data which is the component of the first frequency by performing Fourier transform on the data based on the signal output from the quadrature detector, and the increased frequency component data whose frequency is increased from the first frequency by the vibrator It comprises. Furthermore, a discrete Fourier transformer that outputs reduced frequency component data whose frequency has been reduced from the first frequency by the vibrator, and a compensation calculation unit that performs a correction operation on the increased frequency component data and the reduced frequency component data are provided. . And an adder for adding the increased frequency component data and the decreased frequency component data corrected by the compensation calculation unit to the fundamental wave component data.
本発明によれば、探査用キャリア信号の周波数を操作せずに、極めて簡便な技術要素を追加するだけで、容易に空間分解能を向上できる、映像装置及び映像形成方法を提供できる。
上記した以外の課題、構成及び効果は、以下の実施形態の説明により明らかにされる。
According to the present invention, it is possible to provide a video apparatus and a video forming method that can easily improve the spatial resolution by adding a very simple technical element without manipulating the frequency of the carrier signal for search.
Problems, configurations, and effects other than those described above will be clarified by the following description of embodiments.
［超音波映像装置１０１の全体構成］
図１Ａ及びＢは、本発明の実施形態である超音波映像装置１０１の全体構成を示すブロック図である。
図１Ａは、超音波映像装置１０１の概略を示すブロック図である。
超音波映像装置１０１は、超音波プローブ１０２と、信号処理装置１０３と、表示装置１０４で構成される。
本実施形態の超音波プローブ１０２は、筐体の横側に加振器１０５が取り付けられている。
信号処理装置１０３は、測定対象１０６に接触する超音波プローブ１０２を通じて超音波の送受信を行うと共に、加振器１０５を駆動する。そして、受信した超音波の信号に基づいて、映像情報を作成し、表示装置１０４に送出する。
表示装置１０４は一般的なディスプレイの他、パソコンも利用できる。
[Overall configuration of ultrasonic imaging apparatus 101]
1A and 1B are block diagrams showing the overall configuration of an ultrasound imaging apparatus 101 that is an embodiment of the present invention.
FIG. 1A is a block diagram illustrating an outline of the ultrasound imaging apparatus 101.
The ultrasound imaging apparatus 101 includes an ultrasound probe 102, a signal processing device 103, and a display device 104.
In the ultrasonic probe 102 of this embodiment, a vibrator 105 is attached to the side of the casing.
The signal processing device 103 transmits and receives ultrasonic waves through the ultrasonic probe 102 that is in contact with the measurement target 106 and drives the vibrator 105. Then, image information is created based on the received ultrasonic signal and sent to the display device 104.
The display device 104 can use a personal computer in addition to a general display.
図１Ｂは、超音波映像装置１０１の機能をより詳細に示すブロック図である。
信号処理装置１０３は、信号処理部１０７と映像作成部１０８を有する。
信号処理部１０７は、超音波プローブ１０２を通じて超音波の送受信を行うと共に、加振器１０５を駆動する。そして、複素成分データ１０９を生成する。
映像作成部１０８は、複素成分データ１０９に基づいて、映像情報を作成し、表示装置１０４に送出する。
FIG. 1B is a block diagram showing the functions of the ultrasound imaging apparatus 101 in more detail.
The signal processing device 103 includes a signal processing unit 107 and a video creation unit 108.
The signal processing unit 107 transmits and receives ultrasonic waves through the ultrasonic probe 102 and drives the vibrator 105. Then, complex component data 109 is generated.
The video creation unit 108 creates video information based on the complex component data 109 and sends it to the display device 104.
［超音波プローブ１０２の構成］
図２は、超音波プローブ１０２の一部断面図である。超音波プローブ１０２には、、超音波トランスデューサ２０１のスキャン方向の延長線上に加振器１０５が取り付けられている。加振器１０５は超音波トランスデューサ２０１が超音波を送受信する方向と平行な方向に振動する。この振動は、映像を捉える対象である人体２０２の皮膚２０３に密着させた際、人体２０２を振動させ、ずり弾性波を励起する。励起したずり弾性波は横波として人体２０２内を伝搬し、人体２０２の体内の超音波散乱体を振動させ、反射超音波にドップラー効果を与える。このドップラー効果によって、超音波トランスデューサ２０１が受信する超音波の周波数は、周波数変調或は位相変調される。
[Configuration of Ultrasonic Probe 102]
FIG. 2 is a partial cross-sectional view of the ultrasonic probe 102. A vibrator 105 is attached to the ultrasonic probe 102 on an extension line in the scanning direction of the ultrasonic transducer 201. The vibrator 105 vibrates in a direction parallel to the direction in which the ultrasonic transducer 201 transmits and receives ultrasonic waves. This vibration vibrates the human body 202 when it is brought into close contact with the skin 203 of the human body 202 that is a target for capturing an image, and excites shear elastic waves. The excited shear elastic wave propagates in the human body 202 as a transverse wave, vibrates an ultrasonic scatterer in the human body 202, and gives a Doppler effect to the reflected ultrasonic wave. Due to the Doppler effect, the frequency of the ultrasonic wave received by the ultrasonic transducer 201 is frequency-modulated or phase-modulated.
［信号処理部１０７の構成］
図３は、信号処理部１０７の機能ブロック図である。
発振器３０１は、超音波プローブ１０２の内部に設けられる複数の超音波トランスデューサ２０１を駆動する、一定周波数の駆動信号を出力する。例えば５ＭＨｚの駆動信号を出力する。
この駆動信号は、超音波トランスデューサ２０１の他にも利用される。
先ず、分周器３０２を通じて加振器１０５の駆動を制御する信号として使われる。分周器３０２によって５ＭＨｚの駆動信号は５００Ｈｚまでダウンコンバートされる。これ以降、この加振器１０５の駆動を制御する信号を加振信号と呼ぶ。また、加振信号の波数を加振波数と呼ぶ。
次に、直交検波に用いられる二つの乗算器のうち、第一乗算器３０３には発振器３０１からの駆動信号が直接入力され、第二乗算器３０４には、駆動信号が９０°移相器３０５を通して入力される。
また、発振器３０１からの駆動信号は、逓倍器３０６を通じて第一Ａ／Ｄ変換器３０７及び第二Ａ／Ｄ変換器３０８のサンプリングクロックとして用いられる。そして、駆動信号はタイミングクロックとしてシーケンサ３０９にも入力される。
このように、同一の駆動信号に基づいて加振信号とサンプリングクロックを作成することは、図４にて後述するＤＦＴ４０１において、サンプリングデータの位相を揃えるために必要な構成である。
[Configuration of Signal Processing Unit 107]
FIG. 3 is a functional block diagram of the signal processing unit 107.
The oscillator 301 outputs a drive signal having a constant frequency that drives a plurality of ultrasonic transducers 201 provided inside the ultrasonic probe 102. For example, a 5 MHz drive signal is output.
This drive signal is used in addition to the ultrasonic transducer 201.
First, it is used as a signal for controlling the drive of the vibrator 105 through the frequency divider 302. The 5 MHz drive signal is down-converted to 500 Hz by the frequency divider 302. Hereinafter, a signal for controlling the driving of the vibrator 105 is referred to as a vibration signal. The wave number of the excitation signal is referred to as the excitation wave number.
Next, of the two multipliers used for quadrature detection, the driving signal from the oscillator 301 is directly input to the first multiplier 303, and the driving signal is input to the second multiplier 304 by the 90 ° phase shifter 305. Is entered through.
The drive signal from the oscillator 301 is used as a sampling clock for the first A / D converter 307 and the second A / D converter 308 through the multiplier 306. The drive signal is also input to the sequencer 309 as a timing clock.
Thus, the creation of the excitation signal and the sampling clock based on the same drive signal is a configuration necessary for aligning the phases of the sampling data in the DFT 401 described later with reference to FIG.
超音波プローブ１０２には、加振器１０５と複数の超音波トランスデューサ２０１が収納されている。複数の超音波トランスデューサ２０１はセレクタ３１０によって選択される。セレクタ３１０はシーケンサ３０９によって制御される。
セレクタ３１０は更に、シーケンサ３０９によって制御される切替スイッチ３１１に接続されている。切替スイッチ３１１は三点の端子を有し、一点は発振器３０１に、もう一点は第一乗算器３０３及び第二乗算器３０４に接続され、もう一点は何処にも接続されない。切替スイッチ３１１は超音波トランスデューサ２０１に対し、超音波の送信と受信を切り替える。
The ultrasonic probe 102 houses a vibrator 105 and a plurality of ultrasonic transducers 201. The plurality of ultrasonic transducers 201 are selected by the selector 310. The selector 310 is controlled by the sequencer 309.
The selector 310 is further connected to a changeover switch 311 controlled by the sequencer 309. The changeover switch 311 has three terminals, one point is connected to the oscillator 301, the other point is connected to the first multiplier 303 and the second multiplier 304, and the other point is not connected anywhere. The changeover switch 311 switches the transmission and reception of ultrasonic waves to the ultrasonic transducer 201.
切替スイッチ３１１を通じて超音波トランスデューサ２０１から出力される受信信号は、直交検波器を構成する第一乗算器３０３及び第二乗算器３０４に入力される。
第一乗算器３０３は、超音波トランスデューサ２０１からの受信信号と発振器３０１からの駆動信号を乗算する。
第二乗算器３０４は超音波トランスデューサ２０１からの受信信号と、９０°移相器３０５によって位相が９０°遅れた駆動信号を乗算する。
第一乗算器３０３の出力信号は第一ＬＰＦ３１２によって不必要な高周波成分を除去した後、第一Ａ／Ｄ変換器３０７によってデジタルデータに変換される。この、第一Ａ／Ｄ変換器３０７の出力データは、実数成分に相当する「In-phaseデータ」である。
第二乗算器３０４の出力信号は第二ＬＰＦ３１３によって不必要な高周波成分を除去した後、第二Ａ／Ｄ変換器３０８によってデジタルデータに変換される。この、第二Ａ／Ｄ変換器３０８の出力データは、虚数成分に相当する「Quadrature phaseデータ」である。
In-phaseデータとQuadrature phaseデータを合わせたデータを、複素成分データ１０９として、図示しないＲＡＭに一旦保持する。
A reception signal output from the ultrasonic transducer 201 through the changeover switch 311 is input to the first multiplier 303 and the second multiplier 304 that constitute the quadrature detector.
The first multiplier 303 multiplies the reception signal from the ultrasonic transducer 201 and the drive signal from the oscillator 301.
The second multiplier 304 multiplies the received signal from the ultrasonic transducer 201 by the drive signal whose phase is delayed by 90 ° by the 90 ° phase shifter 305.
The output signal of the first multiplier 303 is converted into digital data by the first A / D converter 307 after unnecessary high frequency components are removed by the first LPF 312. The output data of the first A / D converter 307 is “In-phase data” corresponding to the real number component.
The output signal of the second multiplier 304 is converted into digital data by the second A / D converter 308 after unnecessary high frequency components are removed by the second LPF 313. The output data of the second A / D converter 308 is “Quadrature phase data” corresponding to the imaginary component.
Data obtained by combining In-phase data and Quadrature phase data is temporarily held in a RAM (not shown) as complex component data 109.
［映像作成部１０８の構成］
図４は、映像作成部１０８の機能ブロック図である。
複素成分データ１０９は、ＤＦＴ（Discrete Fourier Transform: 離散フーリエ変換器）４０１を用いて、三つのデータに展開される。
第一のデータは、基本周波数成分のデータ（基本波成分データ）である、０周波数成分複素ビットマップデータ４０２である。
第二のデータは、基本周波数成分に対して、加振波数だけプラス方向にシフトしたデータ（正の加振周波数成分データ）である、＋ｆｖ周波数成分複素ビットマップデータ４０３である。
第三のデータは、基本周波数成分に対して、加振波数だけマイナス方向にシフトしたデータ（負の加振周波数成分データ）である、−ｆｖ周波数成分複素ビットマップデータ４０４である。
[Configuration of Video Creation Unit 108]
FIG. 4 is a functional block diagram of the video creation unit 108.
The complex component data 109 is developed into three data using a DFT (Discrete Fourier Transform) 401.
The first data is 0-frequency component complex bitmap data 402, which is fundamental frequency component data (fundamental wave component data).
The second data is + fv frequency component complex bitmap data 403, which is data (positive excitation frequency component data) shifted in the positive direction by the number of excitation waves with respect to the fundamental frequency component.
The third data is −fv frequency component complex bitmap data 404, which is data (negative excitation frequency component data) shifted in the minus direction by the number of excitation waves with respect to the fundamental frequency component.
加振波数推定部４０５は、第二のデータ（＋ｆｖ周波数成分複素ビットマップデータ４０３）、及び第三のデータ（−ｆｖ周波数成分複素ビットマップデータ４０４）に含まれている空間変調によりダウンコンバートされた高い波数スペクトルの情報を、復調操作により元の波数帯域へ戻すために、以下の式で求められる配列データ（空間復調項）を作成する。 The excitation wave number estimation unit 405 is down-converted by spatial modulation included in the second data (+ fv frequency component complex bitmap data 403) and the third data (-fv frequency component complex bitmap data 404). In order to return the information of the high wave number spectrum to the original wave number band by the demodulation operation, array data (spatial demodulation term) obtained by the following equation is created.
振動振幅推定部４０６は、＋ｆｖ周波数成分複素ビットマップデータ４０３及び−ｆｖ周波数成分複素ビットマップデータ４０４の振幅を、０周波数成分複素ビットマップデータ４０２に合わせるための係数ｋδ（スカラ値）を作成する。
除算器４０７は、加振波数推定部４０５が出力する配列データである空間復調項に、振動振幅推定部４０６が出力する係数ｋδを除算して、以下の式に示す、配列データである＋ｆｖでの空間復調振幅補償データ及び−ｆｖでの空間復調振幅補償データを作成する。
The vibration amplitude estimation unit 406 creates a coefficient kδ (scalar value) for matching the amplitudes of the + fv frequency component complex bitmap data 403 and the −fv frequency component complex bitmap data 404 with the 0 frequency component complex bitmap data 402. .
The divider 407 divides the coefficient kδ output from the vibration amplitude estimation unit 406 by the spatial demodulation term, which is the array data output from the excitation wave number estimation unit 405, by + fv, which is the array data shown in the following equation. Spatial demodulation amplitude compensation data and spatial demodulation amplitude compensation data at -fv are created.
第三乗算器４０８は、＋ｆｖ周波数成分複素ビットマップデータ４０３と−ｆｖでの空間復調振幅補償データを乗算する。
第四乗算器４０９は、−ｆｖ周波数成分複素ビットマップデータ４０４と＋ｆｖでの空間復調振幅補償データを乗算する。
加算器４１０は、第三乗算器４０８の出力である、−ｆｖでの空間復調振幅補償データで補正された＋ｆｖ周波数成分複素ビットマップデータ４０３と、第四乗算器４０９の出力である、＋ｆｖでの空間復調振幅補償データで補正された−ｆｖ周波数成分複素ビットマップデータ４０４と、０周波数成分複素ビットマップデータ４０２を加算して、ビデオメモリ（以下「ＶＲＡＭ」と略す）４１１に展開する。そして、このＶＲＡＭ４１１に展開されたビットマップデータは、表示装置１０４によって表示される。
なお、加振波数推定部４０５、振動振幅推定部４０６、除算器４０７、第三乗算器４０８及び第四乗算器４０９は、０周波数成分複素ビットマップデータ４０２に＋ｆｖ周波数成分複素ビットマップデータ４０３及び−ｆｖ周波数成分複素ビットマップデータ４０４を加算器４１０で加算するために振幅補償を行う、補償演算部と総称できる。
The third multiplier 408 multiplies the + fv frequency component complex bitmap data 403 and the spatial demodulation amplitude compensation data at −fv.
The fourth multiplier 409 multiplies the -fv frequency component complex bitmap data 404 by the spatial demodulation amplitude compensation data at + fv.
The adder 410 outputs + fv frequency component complex bitmap data 403 corrected by the spatial demodulation amplitude compensation data at −fv, which is the output of the third multiplier 408, and + fv, which is the output of the fourth multiplier 409. The -fv frequency component complex bitmap data 404 corrected with the spatial demodulation amplitude compensation data of the above and the 0 frequency component complex bitmap data 402 are added and developed in a video memory (hereinafter abbreviated as “VRAM”) 411. The bitmap data developed in the VRAM 411 is displayed by the display device 104.
The excitation wave number estimation unit 405, the vibration amplitude estimation unit 406, the divider 407, the third multiplier 408, and the fourth multiplier 409 include the + fv frequency component complex bitmap data 403 and the 0 frequency component complex bitmap data 402. -Fv Frequency component complex bitmap data 404 can be collectively referred to as a compensation operation unit that performs amplitude compensation in order to add the adder 410.
これより、空間復調振幅補償データについて説明する。
加振により横波が平面波としてｘ方向に伝搬するとすれば、伝搬による空間の変位は以下のように表される。
The spatial demodulation amplitude compensation data will now be described.
If the transverse wave propagates in the x direction as a plane wave by excitation, the displacement of the space due to propagation is expressed as follows.
ここで、δは横波の振動振幅、ｆｖは加振周波数、ｋｖは加振波の波数である。
この波により計測対象が振動し、その位置が変化する。ある波長ｋのセンシング波（たとえば超音波）によって対象を計測するとき、その波長での反射係数は位相変調を受けて、以下に表す式に変化する。
Here, δ is the vibration amplitude of the transverse wave, f v is the excitation frequency, and k v is the wave number of the excitation wave.
This wave vibrates the measurement object and changes its position. When an object is measured by a sensing wave (for example, an ultrasonic wave) having a certain wavelength k, the reflection coefficient at that wavelength undergoes phase modulation and changes to the following expression.
その波長帯域での送受信での点広がり関数をｗ（ｘ）とすると、センシング波の送受信により得られる信号は以下のように表される。 If the point spread function for transmission / reception in the wavelength band is w (x), a signal obtained by transmission / reception of the sensing wave is expressed as follows.
ここで、ｋδ≪１という近似を置けば、センシング波の送受信により得られる信号は以下のように表される。 Here, if an approximation of kδ << 1 is placed, a signal obtained by transmitting and receiving a sensing wave is expressed as follows.
この信号は時間的に変動しない項と正負の周波数±ｆｖで表される項に分けられる。周波数０の項は通常の計測により得られる画像であり、以下のように表される。 This signal is divided into a term that does not vary with time and a term that is represented by positive and negative frequencies ± f v . The term of frequency 0 is an image obtained by normal measurement and is expressed as follows.
この信号は点広がり関数により、低い波数の情報しか得られない。
一方、時間方向にｆｖの逆数の整数倍の時間幅でフーリエ変換すれば、センシング波の送受信により得られる信号は以下の式に変化する。
This signal can obtain only low wave number information by a point spread function.
On the other hand, if the Fourier transform in the time width of an integral multiple of the reciprocal of a time direction to f v, signals obtained by transmitting and receiving the sensing wave is changed to the following equation.
±ｆｖ成分の画像は、やはり点広がり関数により低い波数の情報が抽出されるが、反射係数が波数−（±ｋｖ）の正弦的な変動で空間変調（波数領域ではスペクトルのシフト）されており、等価的に高い波数の情報を含んだ、独立な画像となっている。 In the ± fv component image, low wavenumber information is extracted by the point spread function, but the reflection coefficient is spatially modulated (spectrum shift in the wavenumber domain) with a sinusoidal fluctuation of wavenumber− (± kv ). It is an independent image that contains information of high wave number equivalently.
これらを合成した以下の式が、本実施形態における超解像アルゴリズムの画像である。 The following formula combining these is an image of the super-resolution algorithm in the present embodiment.
加振波数ｋｖをうまく選んで、点広がり関数のスペクトル幅と同程度にすれば、波数領域での帯域幅を以下のように三倍程度に広げることができる。 If the excitation wave number kv is selected appropriately and is made to be about the same as the spectrum width of the point spread function, the bandwidth in the wave number region can be expanded to about three times as follows.
［アナログ信号の態様］
図５は、信号処理部１０７における、アナログ信号の態様を示す波形図である。
発振器３０１は、一定周波数（５ＭＨｚ）及び一定振幅の駆動信号５０１を出力し続ける。これに対し、シーケンサ３０９は切替スイッチ３１１を第一ゲート信号５０２及び第二ゲート信号５０３で制御する。第一ゲート信号５０２は、発振器３０１を超音波トランスデューサ２０１へ接続する信号である。第二ゲート信号５０３は、超音波トランスデューサ２０１を第一乗算器３０３及び第二乗算器３０４へ接続する信号である。つまり、第一ゲート信号５０２は超音波を送信するタイミングのゲート信号であり、第二ゲート信号５０３は超音波を受信するタイミングのゲート信号である。
[Analog signal mode]
FIG. 5 is a waveform diagram showing an aspect of an analog signal in the signal processing unit 107.
The oscillator 301 continues to output a drive signal 501 having a constant frequency (5 MHz) and a constant amplitude. On the other hand, the sequencer 309 controls the changeover switch 311 with the first gate signal 502 and the second gate signal 503. The first gate signal 502 is a signal for connecting the oscillator 301 to the ultrasonic transducer 201. The second gate signal 503 is a signal for connecting the ultrasonic transducer 201 to the first multiplier 303 and the second multiplier 304. That is, the first gate signal 502 is a gate signal at a timing for transmitting an ultrasonic wave, and the second gate signal 503 is a gate signal at a timing for receiving an ultrasonic wave.
第一ゲート信号５０２によって、駆動信号５０１に基づく送信バースト信号５０４が、超音波トランスデューサ２０１に印加される。そして、第一ゲート信号５０２がオフ状態になってから遅れて、第二ゲート信号５０３によって、受信信号５０５が超音波トランスデューサ２０１から第一乗算器３０３及び第二乗算器３０４へ出力される。そして、第一乗算器３０３及び第二乗算器３０４から直交検波信号５０６が得られる。
厳密には、第一乗算器３０３と第二乗算器３０４から出力される信号の波形はそれぞれ異なるが、図５ではその詳細を省略している。
A transmission burst signal 504 based on the drive signal 501 is applied to the ultrasonic transducer 201 by the first gate signal 502. The received signal 505 is output from the ultrasonic transducer 201 to the first multiplier 303 and the second multiplier 304 by the second gate signal 503 with a delay after the first gate signal 502 is turned off. Then, the quadrature detection signal 506 is obtained from the first multiplier 303 and the second multiplier 304.
Strictly speaking, the waveforms of the signals output from the first multiplier 303 and the second multiplier 304 are different from each other, but details thereof are omitted in FIG.
改めて、実際の数値例を伴って説明する。
先ず、第一ゲート信号５０２は駆動信号５０１のゼロクロス点から始まり、４波長分の期間だけオン状態になる。したがって、超音波トランスデューサ２０１は４波長分の送信バースト信号５０４を出力する。
次に、第二ゲート信号５０３は第一ゲート信号５０２がオフ状態になってから、映像化したい対象の、最も浅い位置の深度に合わせた時間だけ遅れて、駆動信号５０１のゼロクロス点から所定の期間だけオン状態になる。
人体２０２における超音波の伝搬速度は凡そ１５００ｍ／ｓｅｃである。例えば、２μｓｅｃだけ第二ゲート信号５０３を遅延させるとすると、２μｓｅｃ×１５００ｍ／ｓｅｃ＝３ｍｍとなる。この３ｍｍという値は、超音波が往復３ｍｍの時間差を形成することを意味し、それは片道１．５ｍｍに相当する。つまり、体表から１．５ｍｍの深さから測定を行うこととなる。
第二ゲート信号５０３は、映像化したい対象の、最も深い位置の深度に合わせた時間だけ、オン状態にする。例えば、第二ゲート信号の遅延時間を０．１ｍｓｅｃに設定すると、往復で１５ｃｍ、片道７．５ｃｍの深さになる。この程度の深度にすれば、胎児や臓器の映像化に適したものとなる。
第二ゲート信号５０３はそのまま第一Ａ／Ｄ変換器３０７及び第二Ａ／Ｄ変換器３０８のオン・オフ制御信号としても使用される。
A description will be given again with actual numerical examples.
First, the first gate signal 502 starts from the zero cross point of the drive signal 501 and is turned on for a period of four wavelengths. Therefore, the ultrasonic transducer 201 outputs a transmission burst signal 504 for four wavelengths.
Next, after the first gate signal 502 is turned off, the second gate signal 503 is delayed from the zero cross point of the drive signal 501 by a time corresponding to the depth of the shallowest position of the object to be imaged. Turns on only for a period.
The propagation speed of ultrasonic waves in the human body 202 is approximately 1500 m / sec. For example, if the second gate signal 503 is delayed by 2 μsec, then 2 μsec × 1500 m / sec = 3 mm. This value of 3 mm means that the ultrasonic waves form a time difference of 3 mm for reciprocation, which corresponds to 1.5 mm for one way. That is, measurement is performed from a depth of 1.5 mm from the body surface.
The second gate signal 503 is turned on only for a time corresponding to the depth of the deepest position of the object to be imaged. For example, when the delay time of the second gate signal is set to 0.1 msec, the depth is 15 cm in a round trip and 7.5 cm in one way. This depth is suitable for imaging fetuses and organs.
The second gate signal 503 is directly used as an on / off control signal for the first A / D converter 307 and the second A / D converter 308.
［複素成分データ１０９と複素ビットマップデータの関係］
図６は、映像作成部１０８における、複素成分データ１０９と、ＤＦＴ４０１によって変換された複素ビットマップデータとの関係を示す図である。
先の図５で説明した、第二ゲート信号５０３に基づく複素成分データ１０９は、図６中、Ｎ個のサンプル数の配列データとなる。ここでＮは、第二ゲート信号５０３がオン状態の期間に含まれる、サンプリングパルスの数である。
この配列データを、一つの超音波トランスデューサ２０１について、Ｍ回繰り返して取得する。したがって、ある超音波トランスデューサ２０１について、Ｎ個の配列データがＭ個得られる。図６中、「パルス繰り返し方向」と記されている矢印が、この繰り返しを示す。
[Relationship between Complex Component Data 109 and Complex Bitmap Data]
FIG. 6 is a diagram illustrating a relationship between the complex component data 109 and the complex bitmap data converted by the DFT 401 in the video creation unit 108.
The complex component data 109 based on the second gate signal 503 described with reference to FIG. 5 is array data of N samples in FIG. Here, N is the number of sampling pulses included in the period in which the second gate signal 503 is on.
This array data is repeatedly acquired M times for one ultrasonic transducer 201. Therefore, M array data are obtained for an ultrasonic transducer 201. In FIG. 6, an arrow marked “pulse repetition direction” indicates this repetition.
サンプリングパルスは、第一ゲート信号５０２が形成する送信バースト信号５０４が体内で反射して戻ってきた際の、深度に相当する。
そこで、同じサンプリングパルスのデータを対象に、ＤＦＴ４０１にて周波数領域のデータに変換する。すると、幾つかの周波数成分のデータが得られる。
従来技術の超音波映像装置では、基本周波数である駆動信号５０１の周波数成分（例えば５ＭＨｚ）のデータだけを取り出して、映像化していた。これが０周波数成分複素ビットマップデータ４０２である。
本実施形態の超音波映像装置は、この駆動信号５０１の周波数成分に加え、加振器１０５の振動周波数だけプラス方向とマイナス方向にシフトした周波数成分のデータを取り出す。これが、＋ｆｖ周波数成分複素ビットマップデータ４０３と、−ｆｖ周波数成分複素ビットマップデータ４０４である。図６中、「ドップラー周波数方向」と記されている矢印が、−ｆｖ周波数成分複素ビットマップデータ４０４と、０周波数成分複素ビットマップデータ４０２と、＋ｆｖ周波数成分複素ビットマップデータ４０３の並び、すなわち周波数の相違を示す。
The sampling pulse corresponds to the depth when the transmission burst signal 504 formed by the first gate signal 502 is reflected back inside the body.
Therefore, DFT 401 converts the data of the same sampling pulse into data in the frequency domain. Then, data of several frequency components are obtained.
In the conventional ultrasound imaging apparatus, only the data of the frequency component (for example, 5 MHz) of the drive signal 501 that is the fundamental frequency is extracted and visualized. This is the zero frequency component complex bitmap data 402.
In addition to the frequency component of the drive signal 501, the ultrasonic imaging apparatus according to the present embodiment extracts data of a frequency component shifted in the positive direction and the negative direction by the vibration frequency of the vibrator 105. This is + fv frequency component complex bitmap data 403 and -fv frequency component complex bitmap data 404. In FIG. 6, an arrow marked “Doppler frequency direction” indicates an arrangement of −fv frequency component complex bitmap data 404, 0 frequency component complex bitmap data 402, and + fv frequency component complex bitmap data 403, that is, Indicates the difference in frequency.
こうして、一つの超音波トランスデューサ２０１について、深さ方向に要素Ｎ個の一次元の複素成分データ１０９がＭ個得られ、更にＤＦＴ４０１によって要素Ｎ個の一次元の複素ビットマップデータが三個得られる。この作業を、超音波プローブ１０２の内部にＫ個設けられている超音波トランスデューサ２０１に対して、セレクタ３１０を切り替えて実行する。図６中、「スキャン方向」と記されている矢印が、Ｋ個の超音波トランスデューサ２０１から、セレクタ３１０を選択することによって形成されるデータを示す。
ＤＦＴ４０１によって、時間関数である複素成分データ１０９が複素ビットマップデータに変換される際、測定の繰り返し数であるＭという要素は除去される。
In this way, M one-dimensional complex component data 109 of N elements in the depth direction is obtained for one ultrasonic transducer 201, and three pieces of one-dimensional complex bitmap data of N elements are obtained by DFT 401. . This operation is executed by switching the selector 310 for the K ultrasonic transducers 201 provided inside the ultrasonic probe 102. In FIG. 6, an arrow marked “scan direction” indicates data formed by selecting the selector 310 from the K ultrasonic transducers 201.
When the complex component data 109, which is a time function, is converted into complex bitmap data by the DFT 401, the element M, which is the number of measurement repetitions, is removed.
［複素ビットマップデータの補正演算］
本実施形態の超音波映像装置は、空間分解能を向上させるために、計測対象を超音波の進行方向と平行な方向へ加振し、超音波トランスデューサのスキャン方向に横波を伝搬させる。このように加振しながらそのドップラー周波数成分を計測することで、空間周波数のカバー範囲を拡張することができる。
図７は、映像作成部１０８における、複素ビットマップデータの補正演算を説明するための模式図である。図７は、図４のブロック図とほぼ同じ態様ではあるが、データの態様をより明確に表している。
図８Ａ、Ｂ及びＣは、超音波トランスデューサ２０１の点拡がり関数のスペクトルと、復調時のドップラー成分と、振幅の差を示す図である。
[Complex bitmap data correction]
In order to improve the spatial resolution, the ultrasonic imaging apparatus of the present embodiment vibrates the measurement target in a direction parallel to the traveling direction of the ultrasonic wave, and propagates the transverse wave in the scanning direction of the ultrasonic transducer. By measuring the Doppler frequency component while vibrating in this way, the spatial frequency coverage can be expanded.
FIG. 7 is a schematic diagram for explaining the complex bitmap data correction calculation in the video creation unit 108. 7 is almost the same mode as the block diagram of FIG. 4, but more clearly represents the mode of data.
8A, 8B, and 8C are diagrams showing the difference between the spectrum of the point spread function of the ultrasonic transducer 201, the Doppler component at the time of demodulation, and the amplitude.
図７中、０周波数成分複素ビットマップデータ４０２、＋ｆｖ周波数成分複素ビットマップデータ４０３及び−ｆｖ周波数成分複素ビットマップデータ４０４は、それぞれＮ行Ｋ列の配列データである。Ｎは深度（サンプリングパルスの数）であり、Ｋは横方向座標（超音波トランスデューサ２０１の数）である。
加振波数推定部４０５は、加振信号の波数によって駆動信号５０１の周波数がプラス方向とマイナス方向にシフトするシフト量を推定演算する。何故周波数ではなく「波数」であるのかの理由は、変調の形態の違いに起因する。
In FIG. 7, 0 frequency component complex bitmap data 402, + fv frequency component complex bitmap data 403, and -fv frequency component complex bitmap data 404 are array data of N rows and K columns, respectively. N is the depth (number of sampling pulses), and K is the horizontal coordinate (number of ultrasonic transducers 201).
The excitation wave number estimation unit 405 estimates and calculates the shift amount by which the frequency of the drive signal 501 shifts in the plus direction and the minus direction depending on the wave number of the excitation signal. The reason why it is “wave number” instead of frequency is due to the difference in the form of modulation.
通常の信号処理装置１０３では、時間領域における電気信号の正弦波変調をフーリエ変換して、周波数領域における周波数シフトの情報を得る。
しかし、本実施形態の超音波映像装置が扱う信号は、空間領域における超音波に対して、正弦波で駆動される媒質を用いた空間変調を検出した信号である。この空間領域の信号をフーリエ変換すると、波数領域における波数シフトの情報に変わる。つまり、通常の信号処理装置１０３は電気信号が光の速さに等しいので、周期や周波数を容易に特定できるが、本実施形態の超音波映像装置は、波が伝搬する媒質によって波の伝搬速度がまちまちである。このため、媒質における波の伝搬速度を考慮する必要が生じる。
In the normal signal processing apparatus 103, information on frequency shift in the frequency domain is obtained by performing Fourier transform on the sinusoidal modulation of the electric signal in the time domain.
However, the signal handled by the ultrasonic imaging apparatus of the present embodiment is a signal obtained by detecting spatial modulation using a medium driven by a sine wave with respect to ultrasonic waves in the spatial domain. When the signal in this spatial domain is Fourier-transformed, it changes to information on the wave number shift in the wave number domain. That is, since the normal signal processing apparatus 103 can easily specify the period and frequency because the electric signal is equal to the speed of light, the ultrasonic imaging apparatus of the present embodiment has a wave propagation speed depending on the medium through which the wave propagates. There are various towns. For this reason, it is necessary to consider the wave propagation speed in the medium.
図８Ａは、空間分解能Ｄの超音波トランスデューサ２０１における、点拡がり関数のスペクトルを示す図である。図８Ａ中、横軸は波数ｋｘである。横軸方向の２π／Ｄが、実質的に有効な波数範囲である。
図８Ｂは、復調時のドップラー成分を含めた、点拡がり関数のスペクトルを示す図である。図８Ｂ中、λｖは加振波の波長である。点拡がり関数のスペクトルが超音波トランスデューサ２０１の点拡がり関数と同じ大きさになるように、加振波の波長（すなわち加振波の周波数）を定めることで、三つの点拡がり関数のスペクトルが波数領域上でうまく合わさり、有効な波数範囲を拡大できる。このためには、加振波の波長λｖと、センサの点拡がり関数サイズＤが、ほぼ同じ値であることが求められる。
もし、λｖがＤより極端に大きい（λｖ≫Ｄ）場合、点拡がり関数の帯域同士が余計に重なり、合成した帯域幅が狭くなり、実用性が薄れる。
もし、λｖがＤより極端に小さい（λｖ≪Ｄ）場合、スペクトルが平坦にならない。
FIG. 8A is a diagram showing a spectrum of a point spread function in the ultrasonic transducer 201 with a spatial resolution D. FIG. In FIG. 8A, the horizontal axis represents the wave number kx. 2π / D in the horizontal axis direction is a substantially effective wave number range.
FIG. 8B is a diagram illustrating a spectrum of a point spread function including a Doppler component at the time of demodulation. In FIG. 8B, λv is the wavelength of the excitation wave. By determining the wavelength of the excitation wave (that is, the frequency of the excitation wave) so that the spectrum of the point spread function has the same size as the point spread function of the ultrasonic transducer 201, the spectrum of the three point spread functions becomes the wave number. The effective wave number range can be expanded by combining well on the area. For this purpose, it is required that the wavelength λv of the excitation wave and the point spread function size D of the sensor have substantially the same value.
If λv is extremely larger than D (λv >> D), the bands of the point spread function overlap each other, the combined bandwidth becomes narrow, and the practicality is reduced.
If λv is extremely smaller than D (λv << D), the spectrum does not become flat.
図８Ｃは、０周波数成分複素ビットマップデータ４０２と、＋ｆｖ周波数成分複素ビットマップデータ４０３及び−ｆｖ周波数成分複素ビットマップデータ４０４の点拡がり関数のスペクトルを示す図である。この点広がり関数は、超音波トランスデューサのスキャン方向における指向性の特性をフーリエ変換して得られるものである。そして、この点広がり関数の外形は、ガウス関数で近似モデル化できる。
加振波数推定部４０５が出力する配列データは、＋ｆｖ周波数成分複素ビットマップデータ４０３及び−ｆｖ周波数成分複素ビットマップデータ４０４に対し、周波数シフトを戻して、超音波の基本周波数成分に一致させる。
但し、＋ｆｖ周波数成分複素ビットマップデータ４０３及び−ｆｖ周波数成分複素ビットマップデータ４０４の振幅は、図８Ｃに示すように、０周波数成分複素ビットマップデータ４０２の振幅と比べると、ｋδだけ低下している。ｋは超音波の波数であり、δは加振波の振動振幅である。ｋδが小さければ小さいほど、＋ｆｖ周波数成分複素ビットマップデータ４０３及び−ｆｖ周波数成分複素ビットマップデータ４０４の振幅は低下してしまう。したがって、加振波の振動振幅が小さすぎるとノイズに埋もれてしまう。
帯域合成の際、０周波数成分複素ビットマップデータとこの振幅の差を揃えるのが振動振幅推定部４０６である。
FIG. 8C is a diagram illustrating the spectrum of the point spread function of the 0 frequency component complex bitmap data 402, the + fv frequency component complex bitmap data 403, and the −fv frequency component complex bitmap data 404. This point spread function is obtained by Fourier transforming the directivity characteristic in the scanning direction of the ultrasonic transducer. The outline of this point spread function can be approximated by a Gaussian function.
The array data output by the excitation wave number estimation unit 405 returns the frequency shift to the + fv frequency component complex bitmap data 403 and the −fv frequency component complex bitmap data 404 to match the fundamental frequency component of the ultrasonic wave.
However, the amplitudes of the + fv frequency component complex bitmap data 403 and the −fv frequency component complex bitmap data 404 are lower by kδ than the amplitude of the 0 frequency component complex bitmap data 402, as shown in FIG. 8C. Yes. k is the wave number of the ultrasonic wave, and δ is the vibration amplitude of the excitation wave. The smaller kδ is, the lower the amplitude of the + fv frequency component complex bitmap data 403 and the −fv frequency component complex bitmap data 404 is. Therefore, if the vibration amplitude of the excitation wave is too small, it will be buried in noise.
At the time of band synthesis, the vibration amplitude estimation unit 406 aligns the difference between the amplitude and the zero frequency component complex bitmap data.
図９Ａ、Ｂ及びＣは、シミュレーションを用いた、従来技術と本実施形態の比較結果を示す図である。
図９Ａは、シミュレーションの説明図である。ある空間において、超音波プローブ装置から深さ35mmにターゲットとなる点散乱体を２個、１．２５ｍｍの間隔を空けて並べる。この点散乱体を、超音波プローブ装置を用いて検出する。計測用超音波の中心周波数を５ＭＨｚ、４バーストのパルス波とし、送信パルス波形を１ｎｓの離散時間間隔で作成した。超音波プローブ装置の空間分解能は２ｍｍとなるよう超音波トランスデューサの開口面の音圧分布を正弦波の半波長で模擬した。超音波の減衰は１ｄＢ／ＭＨｚ／ｃｍ、伝搬速度を１５００ｍ／ｓｅｃとした。受信波形は直交検波後１．５ＭＨｚの帯域幅のローパスフィルタを通して複素成分データを得る。超音波パルスの繰り返し周波数を１０ｋＨｚとし、１０ｍｓの間１００個の超音波パルスからの散乱波を用いて、ドップラー信号取得のためのフーリエ変換を行う。さらに、超音波トランスデューサの移動間隔を０．５ｍｍとして１２ｍｍ移動させて映像化するものとした。
図９Ｂは、従来技術によるシミュレーション結果の映像である。
従来技術、すなわち０周波数成分複素ビットマップデータ４０２は、図９Ｂに示すように、超音波プローブ装置の空間分解能より狭い間隔である二つの点が結合して一つの点として映像が形成されてしまう。
図９Ｃは、加振を行った場合のシミュレーション結果の映像である。加振周波数は５００Ｈｚ、波長は２ｍｍ、振動振幅は５μｍとした。
本実施形態の場合は、図９Ｃに示すように、同じ超音波プローブ装置を用いていながら、空間分解能が三倍になったことで、二つの点が明確に分離した映像が形成される。
9A, 9B, and 9C are diagrams showing a comparison result between the related art and the present embodiment using a simulation.
FIG. 9A is an explanatory diagram of the simulation. In a certain space, two point scatterers to be targeted are arranged at a depth of 35 mm from the ultrasonic probe device, with an interval of 1.25 mm. This point scatterer is detected using an ultrasonic probe apparatus. The center frequency of the ultrasonic waves for measurement was 5 MHz and 4 burst pulse waves, and the transmission pulse waveform was created at discrete time intervals of 1 ns. The sound pressure distribution on the opening surface of the ultrasonic transducer was simulated with a half wavelength of a sine wave so that the spatial resolution of the ultrasonic probe device was 2 mm. The attenuation of the ultrasonic wave was 1 dB / MHz / cm, and the propagation speed was 1500 m / sec. The received waveform is subjected to quadrature detection, and complex component data is obtained through a low-pass filter having a bandwidth of 1.5 MHz. The repetition frequency of the ultrasonic pulse is set to 10 kHz, and the Fourier transform for acquiring the Doppler signal is performed using the scattered waves from 100 ultrasonic pulses for 10 ms. Furthermore, the moving distance of the ultrasonic transducer was set to 0.5 mm, and the image was moved by 12 mm for imaging.
FIG. 9B is an image of a simulation result according to the prior art.
As shown in FIG. 9B, in the prior art, that is, the zero frequency component complex bitmap data 402, two points having an interval smaller than the spatial resolution of the ultrasonic probe device are combined to form an image as one point. .
FIG. 9C is an image of a simulation result when excitation is performed. The excitation frequency was 500 Hz, the wavelength was 2 mm, and the vibration amplitude was 5 μm.
In the case of the present embodiment, as shown in FIG. 9C, an image in which the two points are clearly separated is formed by using the same ultrasonic probe device and the spatial resolution being tripled.
以上説明した実施形態には、以下に記す応用例が可能である。
（１）上記の実施形態において、空間分解能を更に向上させるには、二種類以上の周波数で測定対象を加振する。そして、ＤＦＴ４０１で加振周波数に応じた成分を取り出し、補正演算を施す。例えば、上記の実施形態では５００Ｈｚの単一の正弦波で測定対象を加振していたが、１０００Ｈｚも同時に与えて加振させることで、理論的には５倍の解像度を得られる。
In the embodiment described above, the following application examples are possible.
(1) In the above embodiment, in order to further improve the spatial resolution, the measurement object is vibrated at two or more frequencies. Then, a component corresponding to the excitation frequency is extracted by the DFT 401 and a correction operation is performed. For example, in the above-described embodiment, the measurement object is vibrated with a single sine wave of 500 Hz. However, theoretically, a resolution of 5 times can be obtained by simultaneously applying 1000 Hz.
（２）本実施形態では、測定対象を映像化するためのコヒーレントな波動に超音波を用いた、超音波映像装置を開示した。このコヒーレントな波動は超音波に限られない。反射波を受信するセンサの空間分解能を向上させるために、測定対象を適切な周波数と振幅で加振できるのであれば、光や電磁波も利用可能である。 (2) In the present embodiment, an ultrasonic imaging apparatus is disclosed that uses ultrasonic waves for coherent waves for imaging a measurement target. This coherent wave is not limited to ultrasound. In order to improve the spatial resolution of the sensor that receives the reflected wave, light and electromagnetic waves can be used as long as the measurement target can be vibrated with an appropriate frequency and amplitude.
（３）加振波数推定部４０５が出力する空間復調項は、予め算出した値をＲＯＭやハードディスク装置等の不揮発性ストレージに記憶しておいてもよい。振動振幅推定部４０６が与える補正定数ｋδは、人体模型等の模擬被験体を撮影して、表示装置１０４を見ながら、カメラのピントを合わせるように調整してもよい。 (3) The spatial demodulation term output by the excitation wave number estimation unit 405 may store a value calculated in advance in a nonvolatile storage such as a ROM or a hard disk device. The correction constant kδ given by the vibration amplitude estimation unit 406 may be adjusted so that the camera is focused while photographing a simulated subject such as a human body model and viewing the display device 104.
（４）本実施形態では、超音波プローブ１０２の筐体の横に加振器１０５を取り付けていたが、加振器１０５を超音波プローブ１０２と別体で設けてもよい。測定中は、操作者が超音波プローブ１０２と加振器１０５の両方を測定対象である人体にあてがう。加振器１０５から発生し、測定対象である人体に伝搬する横波は、数ｍ／ｓｅｃの速度であるので、超音波プローブ１０２の測定範囲である数ｃｍの幅を０．０１秒程度で伝搬する。リアルタイムでこの横波を追いかけてスキャンすることは、技術的には市販の装置で可能な次元である。また、プローブのスキャンを１００回繰り返しながら計測すると仮定しても、１秒程度で数ｃｍの計測が終わるので、この間できるだけ動かないように手で保持するとしても十分実用面で可能である。すなわち、超音波プローブ１０２と加振器１０５との相対距離を維持する手段として、人間の手を用いて固定しても、実質的に問題となる測定誤差は生じない。 (4) In this embodiment, the vibrator 105 is attached to the side of the casing of the ultrasonic probe 102. However, the vibrator 105 may be provided separately from the ultrasonic probe 102. During the measurement, the operator places both the ultrasonic probe 102 and the vibrator 105 on the human body to be measured. The transverse wave generated from the vibrator 105 and propagating to the human body to be measured has a velocity of several m / sec. Therefore, the width of several cm, which is the measurement range of the ultrasonic probe 102, propagates in about 0.01 seconds. To do. It is technically possible to follow the shear wave and scan in real time with a commercially available device. Further, even if it is assumed that measurement is performed while repeating the scan of the probe 100 times, the measurement of several centimeters is completed in about one second, and even if it is held by hand so that it does not move as much as possible during this time, it is sufficiently practical. That is, even if it fixes using a human hand as a means to maintain the relative distance of the ultrasonic probe 102 and the vibrator 105, the measurement error which becomes a problem substantially does not arise.
本実施形態では、超音波映像装置を開示した。
空間分解能を上げたい方向に、計測対象を加振させる。そして、直交検波で得た複素成分データ１０９にＤＦＴ４０１を行う際、加振周波数成分を取り出し、補正演算を施して基本波成分に加算する。このように比較的簡便なプロセスの追加で、容易に空間分解能を向上できる。
In the present embodiment, an ultrasonic imaging apparatus has been disclosed.
Vibrate the measurement object in the direction you want to increase the spatial resolution. When performing DFT 401 on the complex component data 109 obtained by quadrature detection, an excitation frequency component is extracted, subjected to a correction operation, and added to the fundamental component. Thus, the spatial resolution can be easily improved by adding a relatively simple process.
以上、本発明の実施形態例について説明したが、本発明は上記実施形態例に限定されるものではなく、特許請求の範囲に記載した本発明の要旨を逸脱しない限りにおいて、他の変形例、応用例を含む。
例えば、上記した実施形態例は本発明をわかりやすく説明するために装置及びシステムの構成を詳細且つ具体的に説明したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。また、ある実施例の構成の一部を他の実施例の構成に置き換えることは可能であり、更にはある実施例の構成に他の実施例の構成を加えることも可能である。また、各実施例の構成の一部について、他の構成の追加・削除・置換をすることも可能である。
また、上記の各構成、機能、処理部等は、それらの一部又は全部を、例えば集積回路で設計するなどによりハードウェアで実現してもよい。また、上記の各構成、機能等は、プロセッサがそれぞれの機能を実現するプログラムを解釈し、実行するためのソフトウェアで実現してもよい。各機能を実現するプログラム、テーブル、ファイル等の情報は、メモリや、ハードディスク、ＳＳＤ（Solid State Drive）等の揮発性或は不揮発性のストレージ、または、ＩＣカード、光ディスク等の記録媒体に保持することができる。
また、制御線や情報線は説明上必要と考えられるものを示しており、製品上必ずしもすべての制御線や情報線を示しているとは限らない。実際には殆ど全ての構成が相互に接続されていると考えてもよい。
The embodiment of the present invention has been described above. However, the present invention is not limited to the above-described embodiment, and other modifications may be made without departing from the gist of the present invention described in the claims. Includes application examples.
For example, the above-described exemplary embodiments are detailed and specific descriptions of the configuration of the apparatus and system in order to explain the present invention in an easy-to-understand manner, and are not necessarily limited to those having all the configurations described. . Further, a part of the configuration of one embodiment can be replaced with the configuration of another embodiment, and the configuration of another embodiment can be added to the configuration of one embodiment. Moreover, it is also possible to add, delete, and replace other configurations for a part of the configuration of each embodiment.
Each of the above-described configurations, functions, processing units, and the like may be realized by hardware by designing a part or all of them with, for example, an integrated circuit. Further, each of the above-described configurations, functions, and the like may be realized by software for interpreting and executing a program that realizes each function by the processor. Information such as programs, tables, and files for realizing each function is stored in a memory, a hard disk, a volatile or non-volatile storage such as an SSD (Solid State Drive), or a recording medium such as an IC card or an optical disk. be able to.
In addition, the control lines and information lines are those that are considered necessary for the explanation, and not all the control lines and information lines on the product are necessarily shown. Actually, it may be considered that almost all the components are connected to each other.
１０１…超音波映像装置、１０２…超音波プローブ、１０３…信号処理装置、１０４…表示装置、１０５…加振器、１０６…測定対象、１０７…信号処理部、１０８…映像作成部、１０９…複素成分データ、２０１…超音波トランスデューサ、２０２…人体、２０３…皮膚、３０１…発振器、３０２…分周器、３０３…第一乗算器、３０４…第二乗算器、３０５…移相器、３０６…逓倍器、３０７…第一Ａ／Ｄ変換器、３０８…第二Ａ／Ｄ変換器、３０９…シーケンサ、３１０…セレクタ、３１１…切替スイッチ、３１２…第一ＬＰＦ、３１３…第二ＬＰＦ、４０１…ＤＦＴ、４０２…０周波数成分複素ビットマップデータ、４０３…＋ｆｖ周波数成分複素ビットマップデータ、４０４…−ｆｖ周波数成分複素ビットマップデータ、４０５…加振波数推定部、４０６…振動振幅推定部、４０７…除算器、４０８…第三乗算器、４０９…第四乗算器、４１０…加算器、４１１…ＶＲＡＭ、５０１…駆動信号、５０２…第一ゲート信号、５０３…第二ゲート信号、５０４…送信バースト信号、５０５…受信信号、５０６…直交検波信号 DESCRIPTION OF SYMBOLS 101 ... Ultrasonic imaging device, 102 ... Ultrasonic probe, 103 ... Signal processing apparatus, 104 ... Display apparatus, 105 ... Exciter, 106 ... Measurement object, 107 ... Signal processing part, 108 ... Video production part, 109 ... Complex Component data 201 ... Ultrasonic transducer 202 ... Human body 203 ... Skin 301 ... Oscillator 302 ... Divisor 303 ... First multiplier 304 ... Second multiplier 305 ... Phase shifter 306 ... Multiplication 307 ... first A / D converter, 308 ... second A / D converter, 309 ... sequencer, 310 ... selector, 311 ... switch, 312 ... first LPF, 313 ... second LPF, 401 ... DFT 402 ... 0 frequency component complex bitmap data, 403 ... + fv frequency component complex bitmap data, 404 ...- fv frequency component complex bitmap data, 405 ... Wave number estimation unit, 406 ... vibration amplitude estimation unit, 407 ... divider, 408 ... third multiplier, 409 ... fourth multiplier, 410 ... adder, 411 ... VRAM, 501 ... drive signal, 502 ... first gate Signal 503 ... Second gate signal 504 ... Transmission burst signal 505 ... Reception signal 506 ... Quadrature detection signal
Claims (5)
前記測定対象に第二の周波数の振動を加える加振器と、
前記トランスデューサが出力する信号を検波する直交検波器と、
前記直交検波器が出力する信号に基づくデータにフーリエ変換を施して、前記第一の周波数の成分である基本波成分データと、前記加振器によって前記第一の周波数から周波数が増加した増加周波数成分データと、前記加振器によって前記第一の周波数から周波数が減少した減少周波数成分データとを出力する離散フーリエ変換器と、
前記増加周波数成分データ及び前記減少周波数成分データに補正演算を行う補償演算部と、
前記補償演算部によって補正された前記増加周波数成分データ及び前記減少周波数成分データを前記基本波成分データに加算する加算器と
を具備する映像装置。 A transducer for transmitting a wave of a first frequency to a measurement object and receiving the wave reflected from the measurement object;
An exciter that applies vibrations of a second frequency to the measurement object;
A quadrature detector for detecting a signal output from the transducer;
The fundamental frequency component data that is the component of the first frequency by performing Fourier transform on the data based on the signal output by the quadrature detector, and the increased frequency that is increased from the first frequency by the vibrator A discrete Fourier transformer that outputs component data and reduced frequency component data having a frequency reduced from the first frequency by the vibrator;
A compensation calculation unit for performing a correction calculation on the increased frequency component data and the decreased frequency component data;
An image device comprising: an adder that adds the increased frequency component data and the decreased frequency component data corrected by the compensation calculation unit to the fundamental wave component data.
前記第一の周波数の信号を発生する発振器と、
前記発振器を分周して前記第二の周波数の信号を発生する分周器と
を具備する、請求項１記載の映像装置。 Furthermore,
An oscillator for generating a signal of the first frequency;
The video apparatus according to claim 1, further comprising a frequency divider that divides the oscillator to generate a signal of the second frequency.
前記加振器が発生する振動の波数に相当するデータを推定する加振波数推定部と、
前記加振器が発生する振動の振幅に相当するデータを推定する振動振幅推定部と
を具備する、請求項２記載の映像装置。 The compensation calculation unit
An excitation wave number estimating unit that estimates data corresponding to the wave number of vibration generated by the vibrator;
The video apparatus according to claim 2, further comprising: a vibration amplitude estimation unit that estimates data corresponding to an amplitude of vibration generated by the vibrator.
前記受信信号に直交検波を行う直交検波ステップと、
前記直交検波ステップによって得られる信号に基づくデータにフーリエ変換を施して、前記第一の周波数の成分である基本波成分データと、前記波動受信ステップによって前記第一の周波数から周波数が増加した増加周波数成分データと、前記波動受信ステップによって前記第一の周波数から周波数が減少した減少周波数成分データとを出力する離散フーリエ変換ステップと、
前記増加周波数成分データ及び前記減少周波数成分データに補正演算を行う補正演算ステップと、
前記補正演算ステップによって補正された前記増加周波数成分データ及び前記減少周波数成分データを前記基本波成分データに加算する加算ステップと
を有する映像形成方法。 A wave receiving step of transmitting a wave of a first frequency to a measurement object and applying a vibration of a second frequency to the measurement object, and receiving the wave reflected from the measurement object to obtain a reception signal;
A quadrature detection step for performing quadrature detection on the received signal;
The fundamental frequency component data which is a component of the first frequency by performing Fourier transform on the data based on the signal obtained by the quadrature detection step, and an increased frequency whose frequency is increased from the first frequency by the wave receiving step Discrete Fourier transform step of outputting component data and reduced frequency component data in which the frequency is reduced from the first frequency by the wave receiving step;
A correction calculation step for performing a correction calculation on the increased frequency component data and the decreased frequency component data;
An image forming method comprising: an addition step of adding the increased frequency component data and the decreased frequency component data corrected by the correction calculation step to the fundamental wave component data.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2013099248A JP6161028B2 (en) | 2013-05-09 | 2013-05-09 | Image apparatus and image forming method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2013099248A JP6161028B2 (en) | 2013-05-09 | 2013-05-09 | Image apparatus and image forming method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2014217605A true JP2014217605A (en) | 2014-11-20 |
| JP6161028B2 JP6161028B2 (en) | 2017-07-12 |
Family
ID=51936654
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2013099248A Active JP6161028B2 (en) | 2013-05-09 | 2013-05-09 | Image apparatus and image forming method |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6161028B2 (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH02119851A (en) * | 1988-10-29 | 1990-05-07 | Shimadzu Corp | Ultrasonic inspection device |
| JPH0367165A (en) * | 1989-08-04 | 1991-03-22 | Mitsubishi Electric Corp | Ultrasonic flaw detector |
| JP2002502030A (en) * | 1998-02-03 | 2002-01-22 | ザ トラスティーズ オブ ザ スティーブンス インスティテュート オブ テクノロジー | Method and apparatus for acoustic detection and localization of structural defects or ice on structures |
| JP2013512026A (en) * | 2009-11-25 | 2013-04-11 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Ultrasonic shear wave imaging with focused scanline beamforming |
-
2013
- 2013-05-09 JP JP2013099248A patent/JP6161028B2/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH02119851A (en) * | 1988-10-29 | 1990-05-07 | Shimadzu Corp | Ultrasonic inspection device |
| JPH0367165A (en) * | 1989-08-04 | 1991-03-22 | Mitsubishi Electric Corp | Ultrasonic flaw detector |
| JP2002502030A (en) * | 1998-02-03 | 2002-01-22 | ザ トラスティーズ オブ ザ スティーブンス インスティテュート オブ テクノロジー | Method and apparatus for acoustic detection and localization of structural defects or ice on structures |
| JP2013512026A (en) * | 2009-11-25 | 2013-04-11 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Ultrasonic shear wave imaging with focused scanline beamforming |
Also Published As
| Publication number | Publication date |
|---|---|
| JP6161028B2 (en) | 2017-07-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4931910B2 (en) | Ultrasonic imaging device | |
| CN107106120B (en) | Method for ultrasound elastography by continuous vibration of an ultrasound transducer | |
| JP5373308B2 (en) | Ultrasonic imaging apparatus and ultrasonic imaging method | |
| CN110063749B (en) | Ultrasonic measurement device, ultrasonic imaging device, and ultrasonic measurement method | |
| WO2013051943A1 (en) | Nonlinear imaging with dual band pulse complexes | |
| US20110319756A1 (en) | Method for Ultrasound Vibrometry Using Orthogonal Basis Functions | |
| JP3510025B2 (en) | Ultrasound imaging device | |
| JP6733478B2 (en) | Ultrasonic signal processing device, ultrasonic signal processing method, and ultrasonic diagnostic device | |
| Urban et al. | Modulation of ultrasound to produce multifrequency radiation force | |
| KR102257942B1 (en) | Frequency sweep for acoustic radiation force impulse | |
| Nomura et al. | Feasibility of low-frequency ultrasound imaging using pulse compressed parametric ultrasound | |
| CN107049361A (en) | It is imaged using the speed of sound of shearing wave | |
| Noumura et al. | Feasibility of low-frequency directive sound source with high range resolution using pulse compression technique | |
| JP6498183B2 (en) | Ultrasound imaging system | |
| CN107569254B (en) | Ultrasonic signal processing device, ultrasonic signal processing method, and ultrasonic diagnostic device | |
| KR20160090059A (en) | Ultrasound probe, ultrasound diagnostic apparatus having the same and ultrasound signal generating method | |
| JP5513976B2 (en) | Ultrasonic diagnostic equipment | |
| JP6161028B2 (en) | Image apparatus and image forming method | |
| Archer et al. | Longitudinal motion of focused shear wave beams in soft elastic media | |
| JP5260897B2 (en) | Ultrasonic diagnostic equipment | |
| JP2010110503A (en) | Ultrasonic diagnosing device | |
| JP2004041324A (en) | Ultrasonograph | |
| JP2011240006A (en) | Ultrasonograph | |
| JP2011229557A (en) | Ultrasonograph | |
| JP2011234846A (en) | Ultrasonograph |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20160502 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20170313 |
|
| 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: 20170509 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20170606 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6161028 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |