JP3400247B2 - Ultrasound diagnostic equipment - Google Patents
Ultrasound diagnostic equipmentInfo
- Publication number
- JP3400247B2 JP3400247B2 JP16276896A JP16276896A JP3400247B2 JP 3400247 B2 JP3400247 B2 JP 3400247B2 JP 16276896 A JP16276896 A JP 16276896A JP 16276896 A JP16276896 A JP 16276896A JP 3400247 B2 JP3400247 B2 JP 3400247B2
- Authority
- JP
- Japan
- Prior art keywords
- unit
- storage unit
- ultrasonic
- scanning method
- section
- 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
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Description
【0001】[0001]
【発明の属する技術分野】本発明は、医療診断に用いる
超音波診断装置に関し、特に超音波断層像を高フレーム
レートで書き替えることができる超音波診断装置に関す
るものである。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an ultrasonic diagnostic apparatus used for medical diagnosis, and more particularly to an ultrasonic diagnostic apparatus capable of rewriting an ultrasonic tomographic image at a high frame rate.
【0002】[0002]
【従来の技術】従来文献である『医用超音波機器ハンド
ブック』pp.129−161(コロナ社,1985年
発行)には、リアルタイムの超音波ビーム走査方法とし
てBモード走査とMモード走査が紹介されている。図2
は、Bモード走査での超音波ビーム形成とMモード走査
での超音波ビーム形成を示す説明図である。図2(a)
はBモード走査、図2(b)はMモード走査である。B
モード走査では、超音波ビームを移動させながら撮像対
象の2次元断層像を得る。簡単のために、図2(a)で
は4本の超音波ビームで1画像を構成している。探触子
はA,B,C,Dの位置に超音波ビームを順次形成す
る。ビームが形成される順にビーム1,ビーム2・・・
と表す。超音波送受の間隔をTとおくと、例えばAの位
置では4T間隔で新しいビームが構成される。撮像対象
の2次元断層像も4T間隔で更新される。一方、Mモー
ド走査は超音波ビームを1つの位置、例えばAに固定す
る。探触子は、Aの位置でビーム1,ビーム2・・・と
T間隔で超音波ビームを形成する。なお、図2(a)の
超音波ビーム移動方法はリニア走査であるが、この他に
も口径上の1つの固定位置から扇方向にビームを送出す
るセクタ走査、口径から離れた1つの固定位置から扇方
向にビームを送出するコンベックス走査等がある。2. Description of the Related Art "Medical Ultrasonic Equipment Handbook" pp. 129-161 (Corona Publishing, 1985) introduces B-mode scanning and M-mode scanning as real-time ultrasonic beam scanning methods. Figure 2
FIG. 3 is an explanatory diagram showing ultrasonic beam formation in B mode scanning and ultrasonic beam formation in M mode scanning. Figure 2 (a)
Shows B mode scanning, and FIG. 2B shows M mode scanning. B
In mode scanning, a two-dimensional tomographic image of an imaging target is obtained while moving the ultrasonic beam. For the sake of simplicity, in FIG. 2A, four ultrasonic beams form one image. The probe sequentially forms ultrasonic beams at positions A, B, C and D. Beam 1, beam 2 ...
Express. If the ultrasonic transmission / reception interval is set to T, new beams are formed at 4T intervals at the position A, for example. The two-dimensional tomographic image of the imaging target is also updated at 4T intervals. On the other hand, the M mode scanning fixes the ultrasonic beam at one position, for example, A. The probe forms an ultrasonic beam with a beam 1, a beam 2, ... Although the ultrasonic beam moving method in FIG. 2A is linear scanning, in addition to this, sector scanning for sending a beam in a fan direction from one fixed position on the aperture and one fixed position away from the aperture are also used. There is convex scanning and the like in which a beam is sent in a fan direction from.
【0003】[0003]
【発明が解決しようとする課題】前述のBモード走査と
Mモード走査は、それぞれ異なる利点を有している。す
なわち、Bモード走査には撮像対象の2次元断層像を構
成できる利点がある。しかし、1つの超音波ビーム位置
に着目すると、超音波ビームの更新間隔が4Tであるた
め、その間の撮像対象の変化がわからないという問題が
ある。一方、Mモード走査には1つの超音波ビーム位置
での撮像対象の変化を、超音波の送受信間隔Tで観察で
きる利点がある。しかしながら、撮像対象の2次元断層
像を構成できないという問題がある。これら2つの利点
を併せ持つ走査方法が得られるならば、最も望ましい。
そこで本発明の目的は、これら従来の課題を解決し、撮
像対象の2次元断層像を構成することができるととも
に、1つの超音波ビーム位置における超音波ビームの更
新間隔を超音波送受信間隔に等しくして、高フレームレ
ートで断層像を書き替えることができる超音波診断装置
を提供することにある。The B-mode scanning and the M-mode scanning described above have different advantages. That is, the B-mode scanning has an advantage that a two-dimensional tomographic image of an imaging target can be formed. However, when focusing on one ultrasonic beam position, there is a problem that the change of the imaging target during that period cannot be recognized because the ultrasonic beam update interval is 4T. On the other hand, the M-mode scanning has an advantage that the change of the imaging target at one ultrasonic beam position can be observed at the ultrasonic transmission / reception interval T. However, there is a problem that a two-dimensional tomographic image of the imaging target cannot be constructed. It would be most desirable if a scanning method could be obtained that combines these two advantages.
Therefore, an object of the present invention is to solve these conventional problems, to construct a two-dimensional tomographic image of an imaging target, and to make the ultrasonic beam update interval at one ultrasonic beam position equal to the ultrasonic wave transmission / reception interval. Then, an ultrasonic diagnostic apparatus capable of rewriting a tomographic image at a high frame rate is provided.
【0004】[0004]
【課題を解決するための手段】上記目的を達成するた
め、本発明の超音波診断装置は、被検体内に超音波パル
スを送受信する配列素子からなる探触子(図1の1)
と、該探触子(1)中の超音波を送受信する素子を選択
する口径選択部(同2)と、該口径選択部(2)によっ
て選択された探触子送信素子(1)に送波パルスを与え
る送波部(同3)と、該口径選択部(2)によって選択
された探触子受信素子(1)からの信号を並列にアナロ
グ・ディジタル変換する並列A/D変換部(同4)と、
該並列A/D変換部(4)の出力に遅延を与えるディジ
タル遅延部(同5)と、該ディジタル遅延部(5)の出
力信号を加算し、超音波ビームを形成するディジタル加
算部(同6)とから構成される超音波診断装置におい
て、該口径選択部(2)に対し2通りの超音波ビーム走
査方法を指示する走査方法選択部(同7)を有し,該走
査方法選択部(7)が第1の超音波ビーム走査方法(例
えば、Bモード)を選択した場合には、該ディジタル加
算部(6)の出力が切換部(同8)により第1の記憶部
(同9)に格納され、該走査方法選択部(7)が第2の
超音波ビーム走査方法(例えば、Mモード)を選択した
場合には、該ディジタル加算部(6)の出力が切換部
(8)により第2の記憶部(同10)に格納され、該第
1の記憶部(9)のデータと該第2の記憶部(10)の
データとの相関演算を実施すると共に、該第2の記憶部
(10)のデータを用いて該第1の記憶部(9)のデー
タを補間する相関補間部(同11)を有し、該相関補間
部(11)の出力を第3の記憶部(同12)に格納する
ことを特徴としている。これにより、撮像対象の2次元
断層像を構成することができるとともに、1つの超音波
ビーム位置における超音波ビームの更新間隔を超音波送
受信間隔に等しくすることができる。すなわち、2次元
断層像を高フレームレートで書き替え可能となる。In order to achieve the above object, the ultrasonic diagnostic apparatus of the present invention is a probe (1 in FIG. 1) including an array element for transmitting and receiving ultrasonic pulses in a subject.
And an aperture selection section (2) for selecting an element for transmitting and receiving ultrasonic waves in the probe (1), and a probe transmission element (1) selected by the aperture selection section (2). And a parallel A / D converter for parallel-analog-to-digital converting signals from the probe receiving element (1) selected by the aperture selection section (2), Same as 4),
A digital delay unit (5) that gives a delay to the output of the parallel A / D conversion unit (4) and a digital addition unit (the same as the digital addition unit that forms an ultrasonic beam by adding the output signals of the digital delay unit (5)). And a scanning method selection unit (7) for instructing the aperture selection unit (2) of two ultrasonic beam scanning methods. When (7) selects the first ultrasonic beam scanning method (for example, B mode), the output of the digital addition section (6) is changed by the switching section (8) to the first storage section (9). ) And the scanning method selection unit (7) selects the second ultrasonic beam scanning method (for example, M mode), the output of the digital addition unit (6) is the switching unit (8). Stored in the second storage unit (same as 10) by the Data and the data of the second storage unit (10) are correlated, and the data of the first storage unit (9) is interpolated using the data of the second storage unit (10). It is characterized in that it has a correlation interpolator (11) and stores the output of the correlation interpolator (11) in the third storage (12). Thereby, a two-dimensional tomographic image of the imaging target can be formed, and the update interval of the ultrasonic beam at one ultrasonic beam position can be made equal to the ultrasonic transmission / reception interval. That is, the two-dimensional tomographic image can be rewritten at a high frame rate.
【0005】[0005]
【発明の実施の形態】以下、本発明の実施例を、図面に
より説明する。図1は、本発明の第1の実施例における
超音波診断装置の構成図である。図1において、1は被
検体内に超音波パルスを送信し、反射波を受信する探触
子、2は探触子中の超音波を送受信する素子を選択する
口径選択部、3は探触子の送信素子に送波パルスを与え
る送波部、4は探触子の受信素子からの信号を並列にア
ナログ・ディジタル変換する並列A/D変換部、5は並
列A/D変換部4の出力に遅延を与えるディジタル遅延
部、6はディジタル遅延部5の出力信号を加算して超音
波ビームを形成するディジタル加算部、7は口径選択部
2に対して2通りの超音波ビーム走査方法を指示する走
査方法選択部、8はディジタル加算部6の出力を第1ま
たは第2の記憶部に切替え接続する切換部、9は第1の
記憶部、10は第2の記憶部、11は第1の記憶部9の
データと第2の記憶部10のデータとの相関演算を行う
とともに、第2の記憶部10のデータを用いて第1の記
憶部9のデータを補間する相関補間部、12は相関補間
部11の出力を格納する第3の記憶部である。BEST MODE FOR CARRYING OUT THE INVENTION Embodiments of the present invention will be described below with reference to the drawings. FIG. 1 is a configuration diagram of an ultrasonic diagnostic apparatus according to the first embodiment of the present invention. In FIG. 1, reference numeral 1 is a probe that transmits an ultrasonic pulse into the subject and receives a reflected wave, 2 is a diameter selection unit that selects an element for transmitting and receiving ultrasonic waves in the probe, and 3 is a probe. A transmitting unit for giving a transmission pulse to the transmitting element of the child, a parallel A / D converting section for analog-digital converting signals from the receiving element of the probe in parallel, and 5 for the parallel A / D converting section 4. A digital delay unit for delaying the output, 6 is a digital adding unit for adding the output signals of the digital delay unit 5 to form an ultrasonic beam, and 7 is a ultrasonic selecting unit for the aperture selecting unit 2. A scanning method selecting section for instructing, 8 a switching section for switching and connecting the output of the digital adding section 6 to the first or second storage section, 9 a first storage section, 10 a second storage section, and 11 a second storage section. The correlation calculation between the data in the first storage unit 9 and the data in the second storage unit 10 is performed. Utotomoni, first correlation interpolation unit for interpolating the data storage unit 9 with the data of the second storage unit 10, 12 denotes a third storage unit for storing the output of the correlation interpolation section 11.
【0006】走査方法選択部7は、最初に口径選択部2
に対しBモード走査を指定する。ここでは、図2(a)
を用いて動作を説明しているので、断層像1枚を構成す
るために必要な超音波送受信回数を4とする。図2
(a)の位置Aに送波ビームを形成するように、口径選
択部2で選択された送信素子群に対して送波部3から打
ち出しパルス(パルス電圧)を与える。撮像対象からの
反射信号を口径選択部2で選択した受信素子群が受ける
と、並列A/D変換部4により素子ごとにA/D変換さ
れる。なお、口径選択部2が選択する送信素子群と受信
素子群が一致する必要はない。例えば、小さい口径の送
信素子で送出して、大きい口径の受信素子で受信して
も、また大きい口径の送信素子で送出して、小さい口径
の受信素子で受信しても、差し支えない。A/D変換さ
れた各信号は、ディジタル遅延部5でフォーカス用遅延
が与えられ、その後にディジタル加算部6で加算されて
位置Aに受波ビームが形成される。ディジタル加算部6
の出力は、切換部8により第1の記憶部9に接続され
る。以上の操作により、第1の記憶部9に超音波ビーム
位置Aでの反射信号が格納される。口径選択部2が送受
信素子群を順次切り換えて、上述と同様の動作を繰り返
すことにより、位置B,C,Dに超音波ビームを形成
し、その反射信号を第1の記憶部9に格納する。これに
より、断層像1枚を構成するために必要なデータが第1
の記憶部9に蓄えられる。The scanning method selection unit 7 firstly includes the aperture selection unit 2.
For B mode scanning. Here, FIG. 2 (a)
Since the operation is described by using, the number of ultrasonic wave transmission / reception required to form one tomographic image is set to 4. Figure 2
In order to form the transmitted beam at the position A in (a), the wave-transmitting unit 3 applies a pulse (pulse voltage) to the transmitting element group selected by the aperture selecting unit 2. When the reception element group selected by the aperture selection unit 2 receives the reflection signal from the imaging target, the parallel A / D conversion unit 4 performs A / D conversion for each element. The transmitting element group and the receiving element group selected by the aperture selection unit 2 do not have to match. For example, there is no problem even if a small-diameter transmitting element transmits and a large-aperture receiving element receives, or a large-aperture transmitting element transmits and a small-aperture receiving element receives. The A / D-converted signals are given a focusing delay by the digital delay unit 5, and then added by the digital addition unit 6 to form a received beam at the position A. Digital adder 6
The output of is connected to the first storage unit 9 by the switching unit 8. By the above operation, the reflection signal at the ultrasonic beam position A is stored in the first storage unit 9. The aperture selection unit 2 sequentially switches the transmission / reception element group and repeats the same operation as described above to form ultrasonic beams at the positions B, C, and D, and store the reflected signal in the first storage unit 9. . As a result, the data necessary to compose one tomographic image is the first
It is stored in the storage unit 9 of.
【0007】図3は、図1において、Bモード走査時に
第1の記憶部に格納されたデータの図である。装置の使
用者は、超音波ビーム位置A,B,C,Dでの反射信号
を繰り返し第1の記憶部9に格納して、自分が撮像対象
を観察したい時間分の断層像構成に必要な全データを第
1の記憶部9に蓄える。ここでは、A→B→C→D→A
→B→C→D→A→B→C→Dと3枚の断層像構成に必
要なデータを格納したとする。このときの第1の記憶部
9でのデータ配置は、図3に示すように、A〜Dの反射
データが3回ずつ配列される。ここでは、格納された反
射信号をs1〜s12とする。1つの反射信号のデータ
数をNとおけば、第1の記憶部9の総データ数は12N
である。例えば、位置Aでの反射信号はs1,s5,s
9,位置Bでの反射信号はs2,s6,s10である。
またs1,s2,s3,s4を輝度変調して2次元表示
すると1枚目の断層像,s5,s6,s7,s8を輝度
変調して2次元表示すると2枚目の断層像が得られる。
この断層像は、一般にBモード画像と呼ばれる。超音波
送受の間隔をTとおけば、断層像間の同一位置における
反射信号(例えば位置Aでは、s1,s5,s9)の時
間間隔は4Tである。時間間隔が4Tでは、その間の変
化がわからないので、本発明では、相関を用いたデータ
補間によりこの断層像間の同一位置における反射信号の
時間間隔をTとして、断層像の書き替えを高速に行い、
高フレームレートにするものである。FIG. 3 is a diagram of data stored in the first storage section during B mode scanning in FIG. The user of the apparatus repeatedly stores the reflection signals at the ultrasonic beam positions A, B, C, and D in the first storage unit 9 so as to configure a tomographic image corresponding to the time he / she wants to observe the imaging target. All data are stored in the first storage unit 9. Here, A → B → C → D → A
It is assumed that the data necessary for constructing three tomographic images in the order of B, C, D, A, B, C, and D are stored. In the data arrangement in the first storage unit 9 at this time, as shown in FIG. 3, reflection data A to D are arranged three times each. Here, the stored reflection signals are s1 to s12. If the number of data of one reflected signal is N, the total number of data of the first storage unit 9 is 12N.
Is. For example, the reflection signals at position A are s1, s5, s
9, the reflection signals at the position B are s2, s6, and s10.
When s1, s2, s3, and s4 are brightness-modulated and two-dimensionally displayed, a first tomographic image is obtained, and when s5, s6, s7, and s8 are brightness-modulated and two-dimensionally displayed, a second tomographic image is obtained.
This tomographic image is generally called a B-mode image. Assuming that the ultrasonic transmission / reception interval is T, the time interval between the reflection signals (for example, s1, s5, and s9 at the position A) between the tomographic images is 4T. When the time interval is 4T, the change between them is unknown. Therefore, in the present invention, the time interval of the reflection signal at the same position between the tomographic images is set to T by the data interpolation using the correlation, and the tomographic image is rewritten at high speed. ,
High frame rate.
【0008】図4は、図3に示す第1の記憶部のデータ
から作成されたBモード行列の図である。第1の記憶部
9に12Nのデータが格納されると、次に相関補間部1
1が第1の記憶部9に格納された反射信号データを、超
音波ビームの位置ごとに行列に変換する。位置Aでの変
換方法を、図4に示す。s1,s5,s9をゼロづめに
より9行N列の行列に変換する。変換後の行列は、再び
第1の記憶部9に保存される。この行列を、以降の説明
ではBモード行列と呼ぶことにする。Bモード行列は、
位置Aでの反射信号を行方向に時間間隔Tで配列したも
のである。ゼロづめされた行は、データの存在しない時
刻の反射信号である。同様に位置B,C,Dにおいて
も、Bモード行列を作成する。例えば、位置BのBモー
ド行列では、図4においてs1,s5,s9がs2,s
6,s10となる。一般に、1つの反射信号のデータ数
がN個、格納した断層像がL枚、断層像間の同一位置に
おける反射信号の時間間隔がKTのとき、Bモード行列
の大きさはK・(L−1)+1行N列となる。次に、走
査方法選択部7は、口径選択部2に対しMモード走査を
指定する。この動作は、上記第1の記憶部9での行列変
換と並列に実行できる。すなわち、Bモード走査を行っ
た後、第1の記憶部9の行列変換を行うのと並行してM
モード走査を行う。あるいは、上記第1の記憶部9での
行列変換より先に実行することもできる。すなわち、B
モード走査を行った後、Mモード走査を行い、次にBモ
ード行列を作成してもよい。FIG. 4 is a diagram of a B-mode matrix created from the data in the first storage section shown in FIG. When 12N data is stored in the first storage unit 9, the correlation interpolation unit 1
1 converts the reflection signal data stored in the first storage unit 9 into a matrix for each position of the ultrasonic beam. The conversion method at the position A is shown in FIG. s1, s5, s9 are converted into a matrix of 9 rows and N columns by zero-padding. The converted matrix is stored again in the first storage unit 9. This matrix will be referred to as a B-mode matrix in the following description. The B-mode matrix is
The reflection signals at the position A are arranged in the row direction at time intervals T. The zero-padded row is the reflected signal at the time when there is no data. Similarly, at the positions B, C and D, a B mode matrix is created. For example, in the B-mode matrix at position B, s1, s5, s9 are s2, s in FIG.
6 and s10. Generally, when the number of data of one reflection signal is N, the number of stored tomographic images is L, and the time interval of the reflection signals at the same position between the tomographic images is KT, the size of the B-mode matrix is K · (L− 1) +1 row and N column. Next, the scanning method selection unit 7 specifies the M mode scanning for the aperture selection unit 2. This operation can be executed in parallel with the matrix conversion in the first storage unit 9. That is, after the B mode scanning is performed, the matrix conversion of the first storage unit 9 is performed in parallel with the M conversion.
Perform mode scan. Alternatively, it may be executed before the matrix conversion in the first storage unit 9. That is, B
After performing the mode scan, the M mode scan may be performed, and then the B mode matrix may be created.
【0009】図5は、Mモード走査時に第2の記憶部に
格納されたデータの説明図である。図2(b)の位置A
に送波ビームを形成するように、口径選択部2で選択し
た送信素子群に送波部3が打ち出しパルスを与える。こ
れ以降、並列A/D変換部4,ディジタル遅延部5,お
よびディジタル加算部6の動作は、Bモード走査のとき
と全く同様である。ディジタル加算部6の出力は、切換
部8により第2の記憶部10に接続され、第2の記憶部
10に超音波ビーム位置Aでの反射信号が格納される。
なお、反射信号のデータ数は、Bモード走査のときと同
じNとする。Mモード走査のため、口径選択部2は送受
信素子群を固定して同様の動作を繰り返し、図4に示す
位置AにおけるBモード行列の第1行と最終行の時間間
隔(ここでは8T)よりも十分に長い時間にわたり、位
置Aでの反射信号を第2の記憶部10に格納する。ここ
では、160本の反射信号を格納したとする。このとき
の第2の記憶部10でのデータ配置は、図5に示すよう
になる。ここでは、格納された反射信号をr1〜r16
0とする。従って、第2の記憶部10の総データ数は1
60Nである。r1〜r160を輝度変調して2次元表
示すると、位置Aにおいて、一般にMモード画像と呼ば
れる断層像が構成できる。図5における隣接反射信号の
時間間隔はTである。FIG. 5 is an explanatory diagram of data stored in the second storage section during M mode scanning. Position A in FIG. 2B
The transmitting unit 3 gives a launch pulse to the transmitting element group selected by the aperture selecting unit 2 so as to form a transmitting beam. After that, the operations of the parallel A / D converter 4, the digital delay unit 5, and the digital adder 6 are exactly the same as those in the B mode scanning. The output of the digital addition unit 6 is connected to the second storage unit 10 by the switching unit 8, and the reflection signal at the ultrasonic beam position A is stored in the second storage unit 10.
The number of data of the reflection signal is N, which is the same as in the B mode scanning. Because of the M mode scanning, the aperture selecting unit 2 fixes the transmitting and receiving element group and repeats the same operation, and from the time interval (here, 8T) between the first row and the last row of the B mode matrix at the position A shown in FIG. The reflected signal at the position A is stored in the second storage unit 10 for a sufficiently long time. Here, it is assumed that 160 reflected signals are stored. The data arrangement in the second storage unit 10 at this time is as shown in FIG. Here, the stored reflection signals are r1 to r16.
Set to 0. Therefore, the total number of data in the second storage unit 10 is 1
It is 60N. When r1 to r160 are brightness-modulated and two-dimensionally displayed, a tomographic image generally called an M-mode image can be formed at the position A. The time interval between adjacent reflection signals in FIG. 5 is T.
【0010】図6は、図5における第2の記憶部のデー
タから作成したMモード行列の図である。第2の記憶部
10に160Nのデータが格納されると、次に相関補間
部11が第2の記憶部10に格納された反射信号データ
を160行N列の行列に変換する。変換後の行列は、再
び第2の記憶部10に保存される。この行列を、これ以
降の説明ではMモード行列と呼ぶ。Mモード行列は、B
モード行列と同様に位置Aでの反射信号を行方向に時間
間隔Tで並べたものである。Bモード行列との相違は、
Tごとの全時刻でデータが存在するために、ゼロづめの
行が存在しない点である。一般に、1つの反射信号のデ
ータ数がN個、1つの超音波ビーム位置において格納し
たMモード走査の反射信号数がK個のとき、Mモード行
列の大きさはK行N列である。位置AにおけるMモード
行列を作成した後、位置B,C,Dでの反射信号160
本を同様に第2の記憶部10に格納し、それぞれの位置
でのMモード行列を作成する。これにより、断層像の高
フレームレート化に必要な2つの行列が第1の記憶部9
と第2の記憶部10に用意される。ここで、相関補間部
11が2つの行列の相関をとる。位置Aを例にとり、B
モード行列の要素をBij、Mモード行列の要素をMi
jとおく。相関補間部11は、次式(1)により行列間
の相関係数を計算する。FIG. 6 is a diagram of an M-mode matrix created from the data in the second storage section in FIG. When 160N data is stored in the second storage unit 10, the correlation interpolation unit 11 next converts the reflected signal data stored in the second storage unit 10 into a matrix of 160 rows and N columns. The converted matrix is again stored in the second storage unit 10. This matrix will be referred to as an M-mode matrix in the following description. The M-mode matrix is B
Similar to the mode matrix, the reflection signals at position A are arranged in the row direction at time intervals T. The difference from the B-mode matrix is
The point is that there is no zero-padded row because data exists at all times for each T. Generally, when the number of data of one reflection signal is N and the number of reflection signals of M mode scanning stored at one ultrasonic beam position is K, the size of the M mode matrix is K rows and N columns. After creating the M-mode matrix at position A, the reflected signals 160 at positions B, C, D
The book is similarly stored in the second storage unit 10 and the M-mode matrix at each position is created. As a result, two matrices necessary for increasing the frame rate of the tomographic image are stored in the first storage unit 9.
And prepared in the second storage unit 10. Here, the correlation interpolation unit 11 calculates the correlation between the two matrices. Taking position A as an example, B
The elements of the mode matrix are Bij, and the elements of the M mode matrix are Mi.
Let j. The correlation interpolator 11 calculates the correlation coefficient between matrices by the following equation (1).
【数1】
ここで、aは0≦a≦151の整数であり、各aにつき
相関係数を計算する。なお、式(1)における前のΣは
図4の1行〜9行の値の積分、後のΣは図4の1列〜N
列の値の積分、Bijは図4のBモード行列の要素であ
り、M(i+a)jは図6のMモード行列の要素であって、M
モードが160行存在するため、a=0〜151にする
ことにより、約18倍の大きさを有するMモード行列の
上をBモード行列が1行ずつ移動しながら、両要素が掛
算され、移動毎の要素積の総和が求められる。[Equation 1] Here, a is an integer of 0 ≦ a ≦ 151, and the correlation coefficient is calculated for each a. In equation (1), Σ before is the integration of the values in the first row to ninth row in FIG. 4, and Σ is the first column to N in FIG.
The integration of the column values, B ij, is the element of the B-mode matrix of FIG. 4, M (i + a) j is the element of the M-mode matrix of FIG.
Since there are 160 rows of modes, by setting a = 0 to 151, both elements are multiplied and moved while the B-mode matrix moves row by row on the M-mode matrix having a size of about 18 times. The sum of element products for each is calculated.
【0011】図7は、Bモード行列とMモード行列の相
関計算を示す図であり、図8は、Bモード行列をMモー
ド行列で補間する方法を示す図である。前式(1)の相
関演算は、図7に模式的に示されている。先ず、図7に
示すように、Mモード行列とBモード行列を列番号が一
致するように重ね合わせ、次にBモード行列を行方向に
1つずつ移動しながら、移動毎に重なった部分の要素積
の総和を求めている。例えばa=100のときに、前式
(1)の相関演算結果が最大になったとする。このと
き、図8に示す方法で、Mモード行列を用いてBモード
行列を補間するのである。図8では、Bモード行列ゼロ
づめ行の各要素を、相関係数が最大となった個所におけ
る乗算相手のMモード行列の各要素で入れ換えている。
すなわち、s1の次のゼロづめの行をr102に、s1
の2行下のゼロづめの行をr103に、s5の前のゼロ
づめの行をr104に、それぞれ置き換える。なお、図
8では、補間後のBモード行列の第1,5,9行をs
1,s5,s9としたが、(s1+r101)/2,
(s5+r105)/2,(s9+r109)/2とす
ることも可能である。補間後のBモード行列は、第3の
記憶部12に保存される。また、相関係数が最大となっ
た個所における乗算相手のMモード行列各要素を用いる
のではなく、相関係数最大の個所と、2番目に大きい個
所の2つのMモード行列の各要素を用いることも可能で
ある。これは、例えば相関係数最大個所がa=100、
2番目に大きい個所がa=50のときには、補間後のB
モード行列第2行を(r102+r52)/2,第3行
を(r103+r53)/2とすることである。FIG. 7 is a diagram showing the correlation calculation of the B-mode matrix and the M-mode matrix, and FIG. 8 is a diagram showing the method of interpolating the B-mode matrix with the M-mode matrix. The correlation calculation of the above equation (1) is schematically shown in FIG. First, as shown in FIG. 7, the M-mode matrix and the B-mode matrix are superposed so that the column numbers are the same, and then the B-mode matrix is moved one by one in the row direction. Seeking the sum of element products. For example, when a = 100, it is assumed that the correlation calculation result of the above equation (1) is maximized. At this time, the B-mode matrix is interpolated using the M-mode matrix by the method shown in FIG. In FIG. 8, each element of the zero-setting row of the B-mode matrix is replaced with each element of the M-mode matrix of the multiplication partner at the position where the correlation coefficient is maximum.
That is, the next zero-padded row of s1 is set to r102 and s1
The zero padding row two rows below is replaced with r103, and the zero padding row before s5 is replaced with r104. In FIG. 8, the first, fifth, and ninth rows of the B-mode matrix after interpolation are s
1, s5, s9, but (s1 + r101) / 2
It is also possible to set (s5 + r105) / 2, (s9 + r109) / 2. The interpolated B-mode matrix is stored in the third storage unit 12. Also, instead of using each element of the M-mode matrix of the multiplication partner at the location where the correlation coefficient is maximum, each element of the two M-mode matrices at the location having the maximum correlation coefficient and the second largest location is used. It is also possible. This means that, for example, the maximum correlation coefficient is a = 100,
When the second largest part is a = 50, B after interpolation
The second row of the mode matrix is (r102 + r52) / 2, and the third row is (r103 + r53) / 2.
【0012】以上により、位置AでのBモード行列の補
間が完了した。同様の方法で位置B,C,DにおけるB
モード行列も補間する。その結果、図3における1枚目
の断層像と2枚目の断層像の間、2枚目の断層像と3枚
目の断層像の間にそれぞれ新たに3枚の断層像が挿入さ
れて、合計9枚の断層像を構成するのに必要なデータが
第3の記憶部12に格納された。これらの9枚の断層像
においては、断層像間の同一位置における反射信号の時
間間隔がTとなっているので、輝度変調して表示するこ
とにより高フレームレートの断層像を得ることができ
る。また、前式(1)では、列要素を一致させBモード
行列とMモード行列の相関を計算したが、他の方法とし
て列方向にもBモード行列を移動して相関を計算するこ
とも可能である。この場合の計算式を、次式(2)に示
す。By the above, the interpolation of the B-mode matrix at the position A is completed. B at positions B, C, D in a similar manner
The mode matrix is also interpolated. As a result, three new tomographic images are newly inserted between the first tomographic image and the second tomographic image in FIG. 3 between the second tomographic image and the third tomographic image, respectively. The data necessary for forming a total of 9 tomographic images are stored in the third storage unit 12. In these nine tomographic images, since the time interval of the reflection signal at the same position between the tomographic images is T, it is possible to obtain a tomographic image with a high frame rate by displaying with luminance modulation. Further, in the above formula (1), the column elements are matched and the correlation between the B-mode matrix and the M-mode matrix is calculated, but as another method, the B-mode matrix can be moved in the column direction to calculate the correlation. Is. The calculation formula in this case is shown in the following formula (2).
【数2】
ここで、aの範囲は前式(1)と同様であり、bは−N
+1≦b≦N−1の整数である。すなわち、先ず同じ行
番号のMモード行列の左側にBモード行列を1列だけ重
ね合った位置に配列し、順次1列ずつ右側に移動させ
る。そして、Mモード行列の右側にBモード行列を1列
だけ重ね合っただけの位置まで移動させ、それぞれの移
動毎に要素積の総和を求める。Bモード行列で乗算の対
象となる最初の列番号がmax(1,1−b)であり、
最後の列番号がmin(N,N−b)である。次に、前
式(1)と同じように、Bモード行列を行方向に1行ず
つ移動しながら、上述の列方向の移動を行って移動毎の
要素積の総和を求めていく。この計算方式は、撮像対象
の深さ方向の平行移動を考えない場合に用いる。なお前
式(2)を用いる場合も、Bモード行列ゼロづめ行の各
要素に新たに代入する値は、あくまでも相関係数が最大
となった個所における乗算相手のMモード行列各要素で
ある。[Equation 2] Here, the range of a is the same as that of the previous expression (1), and b is -N.
It is an integer of + 1 ≦ b ≦ N−1. That is, first, the B-mode matrix is arranged on the left side of the M-mode matrix of the same row number at a position where the B-mode matrix is overlapped by one column, and is sequentially moved to the right by one column. Then, the B-mode matrix is moved to the right side of the M-mode matrix to a position where only one column is overlapped, and the sum of element products is calculated for each movement. The first column number to be multiplied in the B-mode matrix is max (1,1-b),
The last column number is min (N, N-b). Next, as in the previous expression (1), while moving the B-mode matrix row by row in the row direction, the above-mentioned column direction movement is performed to obtain the sum of element products for each movement. This calculation method is used when the parallel movement of the imaging target in the depth direction is not considered. Also in the case of using the above equation (2), the value newly assigned to each element of the zero-setting row of the B-mode matrix is each element of the M-mode matrix of the multiplication partner at the position where the correlation coefficient becomes the maximum.
【0013】以上述べた補間方法は、撮像物体の各超音
波ビームの位置における反射信号が、周期的に同じ変化
を繰り返すことを仮定している。つまり、最初に作成し
たBモード行列と同じ反射信号パターンが、十分に長い
時間Mモード走査を行い、収集した反射信号から作成し
たMモード行列の中に現れることを前提にしている。生
体での代表的な周期としては、心拍,呼吸がある。The above-described interpolation method assumes that the reflection signals at the positions of the ultrasonic beams of the imaged object repeat the same change periodically. In other words, it is premised that the same reflection signal pattern as the B-mode matrix created first appears in the M-mode matrix created from the reflection signals acquired by performing M-mode scanning for a sufficiently long time. Heartbeat and respiration are typical cycles in the living body.
【0014】図9は、本発明の第2の実施例における超
音波診断装置の構成図である。図9において、1は探触
子、2は口径選択部、3は送波部、7は走査方法選択
部、8は切換部、9は第1の記憶部、10は第2の記憶
部、11は相関補間部、12は第3の記憶部、13はア
ナログ遅延部、14はアナログ加算部、15はA/D変
換部15である。図1の診断装置では、探触子1の全素
子信号がディジタル化され、以降の信号処理がディジタ
ルで行われるのに対して、図9の診断装置では、加算に
よる超音波ビームの形成まではアナログ信号処理で行
い、その先をディジタル信号処理とするものである。こ
の方式では、補間処理に移る前まではアナログ信号のま
ま処理を行うが、補間処理はディジタル処理を行うの
で、精度的には前述の方式と何等変わらない。従って、
高フレームレートの断層像を求めるアルゴリズムは、図
1の診断装置と全く同様である。以上の実施例では、ビ
ームの移動方法としてリニア走査を例にとって説明した
が、その他にもセクタ走査,コンベックス走査があり、
これらにも本発明が適用できることは明らかである。な
お、セクタ走査とは、同位置の送受信素子から扇方向に
超音波ビームを送出する方式であり、コンベックス走査
とは、円弧状に並んだ送受信素子を移動させながら扇方
向にビームを送出する方式である。また、第1の記憶部
9、第2の記憶部10、第3の記憶部12を共通のディ
ジタルメモリとして、アドレスを分割して使うことも可
能である。なお、高フレームレート化が実現される本発
明の効果としては、断層画像を形成するラスタ数で決ま
り、例えば1枚の画像が128本のラスタで形成される
場合には、128倍の書き替え速度となり、またラスタ
数が64本、300本の場合には、それぞれ64倍、3
00倍の効果を得ることができる。FIG. 9 is a block diagram of an ultrasonic diagnostic apparatus according to the second embodiment of the present invention. In FIG. 9, 1 is a probe, 2 is a diameter selection unit, 3 is a wave transmission unit, 7 is a scanning method selection unit, 8 is a switching unit, 9 is a first storage unit, 10 is a second storage unit, Reference numeral 11 is a correlation interpolation unit, 12 is a third storage unit, 13 is an analog delay unit, 14 is an analog addition unit, and 15 is an A / D conversion unit 15. In the diagnostic device of FIG. 1, all the element signals of the probe 1 are digitized and the subsequent signal processing is digitally performed, whereas in the diagnostic device of FIG. 9, the ultrasonic beam is not formed by addition. This is performed by analog signal processing, and then the digital signal processing. In this method, the analog signal is processed as it is before the interpolation processing is started. However, since the interpolation processing is digital processing, there is no difference in accuracy from the above method. Therefore,
The algorithm for obtaining a high frame rate tomographic image is exactly the same as that of the diagnostic apparatus in FIG. In the above embodiments, linear scanning is taken as an example of the beam moving method, but other than that, there are sector scanning and convex scanning.
It is obvious that the present invention can be applied to these as well. Note that sector scanning is a method of transmitting ultrasonic beams in the fan direction from transmitting / receiving elements at the same position, and convex scanning is a method of transmitting beams in the fan direction while moving the transmitting / receiving elements arranged in an arc shape. Is. It is also possible to divide the address and use the first storage unit 9, the second storage unit 10, and the third storage unit 12 as a common digital memory. The effect of the present invention for realizing a high frame rate is determined by the number of rasters forming a tomographic image. For example, when one image is formed by 128 rasters, rewriting by 128 times is performed. If the number of rasters is 64 and 300, 64 times and 3 respectively.
The effect of 00 times can be obtained.
【0015】[0015]
【発明の効果】以上説明したように、本発明によれば、
撮像対象の断層像を構成する超音波診断装置において、
Mモード走査で得られた信号とBモード走査で得られた
信号の相関演算から、Bモード走査で得られた信号を補
間することにより、超音波送受信間隔で画面の書きかわ
る高フレームレートの超音波診断装置を実現することが
できる。As described above, according to the present invention,
In an ultrasonic diagnostic apparatus that forms a tomographic image of an imaging target,
By interpolating the signal obtained in the B mode scanning from the correlation calculation of the signal obtained in the M mode scanning and the signal obtained in the B mode scanning, the screen is rewritten at the ultrasonic transmission / reception intervals to achieve a high frame rate. A sound wave diagnostic apparatus can be realized.
【図1】本発明の第1の実施例を示す超音波診断装置の
構成図である。FIG. 1 is a configuration diagram of an ultrasonic diagnostic apparatus showing a first embodiment of the present invention.
【図2】Bモード走査での超音波ビーム形成とMモード
走査での超音波ビーム形成を表す図である。FIG. 2 is a diagram showing ultrasonic beam formation in B mode scanning and ultrasonic beam formation in M mode scanning.
【図3】Bモード走査時に第1の記憶部に格納されたデ
ータを示す図である。FIG. 3 is a diagram showing data stored in a first storage unit during B mode scanning.
【図4】図3における第1の記憶部のデータから作成さ
れたBモード行列の図である。FIG. 4 is a diagram of a B-mode matrix created from the data in the first storage unit in FIG.
【図5】Mモード走査時に第2の記憶部に格納されたデ
ータを示す図である。FIG. 5 is a diagram showing data stored in a second storage unit during M mode scanning.
【図6】図5における第2の記憶部のデータから作成さ
れたMモード行列の図である。FIG. 6 is a diagram of an M-mode matrix created from the data in the second storage unit in FIG.
【図7】Bモード行列とMモード行列の相関計算を表す
図である。FIG. 7 is a diagram showing a correlation calculation of a B-mode matrix and an M-mode matrix.
【図8】Bモード行列をMモード行列で補間する方法を
表す図である。FIG. 8 is a diagram illustrating a method of interpolating a B-mode matrix with an M-mode matrix.
【図9】本発明の第2の実施例を示す超音波診断装置の
構成図である。FIG. 9 is a configuration diagram of an ultrasonic diagnostic apparatus showing a second embodiment of the present invention.
1…探触子、2…口径選択部、3…送波部、4…並列A
/D変換部、5…ディジタル遅延部、6…ディジタル加
算部、7…走査方法選択部、8…切換部、9…第1の記
憶部、10…第2の記憶部、11…相関補間部、12…
第3の記憶部、13…アナログ遅延部、14…アナログ
加算部、15…A/D変換部、s1〜s12…Bモード
走査時の反射信号、r1〜r160…Mモード走査時の
反射信号。1 ... Probe, 2 ... Caliber selection part, 3 ... Wave transmission part, 4 ... Parallel A
/ D conversion unit, 5 ... Digital delay unit, 6 ... Digital addition unit, 7 ... Scan method selection unit, 8 ... Switching unit, 9 ... First storage unit, 10 ... Second storage unit, 11 ... Correlation interpolation unit , 12 ...
Third storage unit, 13 ... Analog delay unit, 14 ... Analog addition unit, 15 ... A / D conversion unit, s1 to s12 ... Reflection signal in B mode scanning, r1 to r160 ... Reflection signal in M mode scanning.
フロントページの続き (56)参考文献 特開 昭56−8041(JP,A) 特開 平2−98344(JP,A) 特開 平2−257943(JP,A) 特開 平5−118821(JP,A) (58)調査した分野(Int.Cl.7,DB名) A61B 8/00 Continuation of front page (56) Reference JP-A-56-8041 (JP, A) JP-A-2-98344 (JP, A) JP-A-2-257943 (JP, A) JP-A-5-118821 (JP , A) (58) Fields investigated (Int.Cl. 7 , DB name) A61B 8/00
Claims (2)
素子からなる探触子と、該探触子中の超音波を送受信す
る素子を選択する口径選択部と、該口径選択部によって
選択された探触子送信素子に送波パルスを与える送波部
と、該口径選択部によって選択された探触子受信素子か
らの信号を並列にアナログ・ディジタル変換する並列A
/D変換部と、該並列A/D変換部の出力に遅延を与え
るディジタル遅延部と、該ディジタル遅延部の出力信号
を加算して超音波ビームを形成するディジタル加算部と
を備えた超音波診断装置において、 該口径選択部に対し、超音波ビームを移動させながら撮
像対象の2次元断層像を得るBモード走査方法、あるい
は該超音波ビームを1つの位置に固定して走査するMモ
ード走査方法を指示する走査方法選択部と、 該走査方法選択部が該Bモード走査方法を選択した場合
には、該ディジタル加算部の出力を格納する第1の記憶
部と、 該走査方法選択部が該Mモード走査方法を選択した場合
には、該ディジタル加算部の出力を格納する第2の記憶
部と、 該第1の記憶部のデータと該第2の記憶部のデータとの
相関演算を実施すると共に、該第2の記憶部のデータを
用いて該第1の記憶部のデータを補間する相関補間部
と、 該相関補間部の出力を格納する第3の記憶部とを有する
ことを特徴とする超音波診断装置。1. A probe comprising an array element for transmitting and receiving ultrasonic pulses in a subject, a caliber selecting section for selecting an element for transmitting and receiving ultrasonic waves in the probe, and selection by the caliber selecting section. And a parallel transmission circuit for analog-digital converting signals from the probe receiving element selected by the aperture selecting section in parallel to the transmitting section for applying a transmission pulse to the selected probe transmitting element.
An ultrasonic wave including an A / D conversion unit, a digital delay unit that delays an output of the parallel A / D conversion unit, and a digital addition unit that adds an output signal of the digital delay unit to form an ultrasonic beam. In the diagnostic device, a B-mode scanning method for obtaining a two-dimensional tomographic image of an imaging object while moving the ultrasonic beam with respect to the aperture selection unit, or an M-mode scanning for scanning the ultrasonic beam while fixing it at one position A scanning method selecting section for instructing a method, a first storage section for storing an output of the digital adding section when the scanning method selecting section selects the B-mode scanning method, and the scanning method selecting section. When the M-mode scanning method is selected, the correlation calculation between the second storage unit that stores the output of the digital addition unit and the data of the first storage unit and the data of the second storage unit is performed. With the implementation An ultrasonic diagnostic device comprising: a correlation interpolating unit that interpolates data of the first storage unit using data of the second storage unit; and a third storage unit that stores an output of the correlation interpolating unit. apparatus.
素子からなる探触子と、該探触子中の超音波を送受信す
る素子を選択する口径選択部と、該口径選択部によって
選択された探触子送信素子に送波パルスを与える送波部
と、該口径選択部によって選択された探触子受信素子か
らの信号に遅延を与えるアナログ遅延部と、該アナログ
遅延部の出力信号を加算して超音波ビームを形成するア
ナログ加算部と、該アナログ加算部の出力をアナログ・
ディジタル変換するA/D変換部とを備えた超音波診断
装置において、 該口径選択部に対し、超音波ビームを移動させながら撮
像対象の2次元断層像を得るBモード走査方法、あるい
は該超音波ビームを1つの位置に固定して走査するMモ
ード走査方法を指示する走査方法選択部と、 該走査方法選択部が該Bモード走査方法を選択した場合
には、該A/D変換部の出力を格納する第1の記憶部
と、 該走査方法選択部が該Mモード走査方法を選択した場合
には、該A/D変換部の出力を格納する第2の記憶部
と、 該第1の記憶部のデータと該第2の記憶部のデータとの
相関演算を実施すると共に、該第2の記憶部のデータを
用いて該第1の記憶部のデータを補間する相関補間部
と、 該相関補間部の出力を格納する第3の記憶部とを有する
ことを特徴とする超音波診断装置。2. A probe comprising an array element for transmitting and receiving ultrasonic pulses within a subject, a caliber selecting section for selecting an element for transmitting and receiving ultrasonic waves in the probe, and selection by the caliber selecting section. Transmitting section for giving a transmission pulse to the probe transmitting element thus selected, an analog delay section for delaying a signal from the probe receiving element selected by the aperture selecting section, and an output signal of the analog delay section Is added to form an ultrasonic beam, and the output of the analog addition unit
In an ultrasonic diagnostic apparatus including an A / D conversion unit for digital conversion, a B-mode scanning method for obtaining a two-dimensional tomographic image of an imaging target while moving an ultrasonic beam to the aperture selection unit, or the ultrasonic wave A scanning method selecting section for instructing an M mode scanning method for fixing the beam at one position and scanning, and an output of the A / D converting section when the scanning method selecting section selects the B mode scanning method. A first storage unit that stores the output, a second storage unit that stores the output of the A / D conversion unit when the scanning method selection unit selects the M-mode scanning method, and A correlation interpolator that performs a correlation operation between the data in the storage unit and the data in the second storage unit, and interpolates the data in the first storage unit using the data in the second storage unit; And a third storage unit for storing the output of the correlation interpolator. Ultrasonic diagnostic apparatus according to claim.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP16276896A JP3400247B2 (en) | 1996-06-24 | 1996-06-24 | Ultrasound diagnostic equipment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP16276896A JP3400247B2 (en) | 1996-06-24 | 1996-06-24 | Ultrasound diagnostic equipment |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH105214A JPH105214A (en) | 1998-01-13 |
JP3400247B2 true JP3400247B2 (en) | 2003-04-28 |
Family
ID=15760861
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP16276896A Expired - Fee Related JP3400247B2 (en) | 1996-06-24 | 1996-06-24 | Ultrasound diagnostic equipment |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3400247B2 (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6553140B2 (en) * | 2013-01-22 | 2019-07-31 | キヤノンメディカルシステムズ株式会社 | Ultrasonic diagnostic apparatus, image processing apparatus, and image processing method |
CN104936531B (en) | 2013-01-22 | 2017-06-09 | 东芝医疗系统株式会社 | Diagnostic ultrasound equipment, image processing apparatus and image processing method |
-
1996
- 1996-06-24 JP JP16276896A patent/JP3400247B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JPH105214A (en) | 1998-01-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1620015B1 (en) | Efficient ultrasound system for two-dimentional c-scan imaging and related method thereof | |
JP2777197B2 (en) | Ultrasound diagnostic equipment | |
KR101303626B1 (en) | Diagnosis system for diagnosing subject, medical image system for providing diagnosis image of subject and method for displaying diagnosis image of subject | |
EP1972280A1 (en) | Ultrasound system and method for forming ultrasound images | |
JP3144819B2 (en) | Ultrasound diagnostic equipment | |
JP2007222390A (en) | Ultrasonograph | |
JP5346245B2 (en) | Ultrasonic diagnostic apparatus, ultrasonic imaging method and program | |
JP2010511420A (en) | Method and apparatus for multi-line color flow and vascular ultrasound imaging | |
JP3400247B2 (en) | Ultrasound diagnostic equipment | |
JP3884370B2 (en) | Ultrasonic diagnostic equipment | |
JP3544722B2 (en) | Ultrasound diagnostic equipment | |
JP3806229B2 (en) | Ultrasonic diagnostic equipment | |
JP4418052B2 (en) | Ultrasonic beam scanning method and apparatus and ultrasonic imaging apparatus | |
JPH0324865B2 (en) | ||
JP2639391B2 (en) | Fluid map image forming device | |
JP3046424B2 (en) | Ultrasound diagnostic equipment | |
JP2007215630A (en) | Ultrasonograph | |
JPH0678922A (en) | Ultrasonic diagnostic device | |
JP4565822B2 (en) | Ultrasonic diagnostic equipment | |
JPH08173431A (en) | Ultrasonic diagnostic device | |
JPH0399651A (en) | Ultrasonic diagnosis apparatus | |
JPH06205773A (en) | Ultrasonic diagnostic system | |
JP3620953B2 (en) | Ultrasonic diagnostic equipment | |
JP2009112491A (en) | Ultrasonic diagnostic system | |
JPH03173548A (en) | Ultrasonic diagnostic device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080221 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090221 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090221 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100221 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100221 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110221 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110221 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120221 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130221 Year of fee payment: 10 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140221 Year of fee payment: 11 |
|
LAPS | Cancellation because of no payment of annual fees |