JP2824673B2 - Ultrasound diagnostic equipment - Google Patents

Ultrasound diagnostic equipment

Info

Publication number
JP2824673B2
JP2824673B2 JP1241327A JP24132789A JP2824673B2 JP 2824673 B2 JP2824673 B2 JP 2824673B2 JP 1241327 A JP1241327 A JP 1241327A JP 24132789 A JP24132789 A JP 24132789A JP 2824673 B2 JP2824673 B2 JP 2824673B2
Authority
JP
Japan
Prior art keywords
sound ray
signal
element group
image
received
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 - Lifetime
Application number
JP1241327A
Other languages
Japanese (ja)
Other versions
JPH03103246A (en
Inventor
一博 河野
芳克 野田
Original Assignee
ジーイー横河メディカルシステム株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by ジーイー横河メディカルシステム株式会社 filed Critical ジーイー横河メディカルシステム株式会社
Priority to JP1241327A priority Critical patent/JP2824673B2/en
Publication of JPH03103246A publication Critical patent/JPH03103246A/en
Application granted granted Critical
Publication of JP2824673B2 publication Critical patent/JP2824673B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Description

【発明の詳細な説明】 (産業上の利用分野) 本発明は超音波診断装置のBモード画像のスペックル
ノイズの実時間低減法に関する。
Description: TECHNICAL FIELD The present invention relates to a method for real-time reduction of speckle noise in a B-mode image of an ultrasonic diagnostic apparatus.

(従来の技術) パルス式超音波診断装置は超音波パルスを被検体内に
放射し、超音波の減衰や反射の度合が組織やその病変部
によって異なることを利用して、被検体からの反射波を
分析して診断する装置である。しかしこの超音波診断装
置で得られたBモード画像にはスペックルと呼ばれる斑
紋状の模様が現れ、画質劣化の主原因となっている。ス
ペックルの発生の原因は生体内に入射される超音波パル
スが位相の揃ったものであるために、被検体内部に存在
する多数の散乱体からの反射波はトランスデューサの受
波面上で干渉を起こしてしまう。その結果、検出された
エコー信号の振幅にゆらぎを生じ、これがBモード画像
上の斑紋状の模様、即ち、スペックルとなって現れる。
画像に現れたスペックルパタンは多数の散乱体からの反
射波の干渉の結果生じたものであり、組織の微細な構造
を表すものではない。従って、受信信号をBモード画像
として表示するためには、スペックルを減少させる必要
がある。
(Prior art) A pulse-type ultrasonic diagnostic apparatus radiates an ultrasonic pulse into a subject, and utilizes the fact that the degree of attenuation and reflection of the ultrasonic wave differs depending on a tissue or a lesion thereof, and is reflected from the subject. It is a device that analyzes and diagnoses waves. However, a mottled pattern called speckle appears in the B-mode image obtained by the ultrasonic diagnostic apparatus, which is a main cause of image quality deterioration. The speckles are generated because the ultrasonic pulses incident on the living body have the same phase, and the reflected waves from a large number of scatterers inside the subject cause interference on the receiving surface of the transducer. Will wake up. As a result, the amplitude of the detected echo signal fluctuates, which appears as a mottled pattern on the B-mode image, that is, speckle.
The speckle pattern that appears in the image is a result of interference of reflected waves from a large number of scatterers and does not represent the fine structure of the tissue. Therefore, in order to display a received signal as a B-mode image, it is necessary to reduce speckle.

(発明が解決しようとする課題) 上記のスペックル減少のために、スペックルに関して
相関の小さい複数個のBモード画像を重ね合わせる方法
を取っており、この方法としてトランスデューサの位置
変化(空間コンパウンド法)、或るいは入射する超音波
パルスの中心周波数変化(周波数コンパウンド法)によ
って得ている。
(Problem to be Solved by the Invention) In order to reduce the above-mentioned speckle, a method of superimposing a plurality of B-mode images having a small correlation with respect to the speckle is adopted. ) Or by changing the center frequency of the incident ultrasonic pulse (frequency compound method).

空間コンパウンド法では、受信アレイトランスデュー
サの受信エレメントの位置を変えることによりスペック
ルに関して相関の小さい複数の信号を得、これをBモー
ド画像化した後、加算して重ね合わせて表示することに
よりスペックルノズルを低減している。この方法では、
アレイトランスデューサによる送受信を異なるエレメン
トで行い、これに増幅,検波、整相加算等の処理を施こ
してBモード画像にし、この行程を必要な回数(最小限
2回)繰り返して加算して画像表示するため、時間が掛
かり、実時間処理が困難であるという問題があった。
又、相関の小さい複数個のBモード画像を得るために
は、複数個の画像を同一領域から得ることができず、ス
ペックルノイズの低減される領域が第5図の斜線で示さ
れる領域に狹まってしまうという問題もあった。図にお
いて、1はアレイトランスデューサで、複数のエレメン
トによって等価的に図のA,B,C3点から照射されている場
合のビームの照射範囲を示している。この場合にスペッ
クルノイズの低減される領域は斜線の範囲のみである。
In the spatial compound method, a plurality of signals having a small correlation with respect to speckles are obtained by changing the position of a receiving element of a receiving array transducer, and these signals are formed into a B-mode image, and then added and superimposed and displayed. Nozzles are reduced. in this way,
Transmission and reception by the array transducer are performed by different elements, and processing such as amplification, detection, and phasing addition is performed to produce a B-mode image. This process is repeated a required number of times (minimum two times) and added to display an image. Therefore, there is a problem that it takes time and real-time processing is difficult.
Further, in order to obtain a plurality of B-mode images having a small correlation, a plurality of images cannot be obtained from the same region, and the region in which the speckle noise is reduced corresponds to the region shown by oblique lines in FIG. There was also a problem that it narrowed. In the drawing, reference numeral 1 denotes an array transducer, which shows a beam irradiation range in a case where irradiation is performed from a plurality of elements equivalently from points A, B, and C3 in the drawing. In this case, the region in which the speckle noise is reduced is only the range of the oblique line.

本発明は上記の点に鑑みてなされたもので、その目的
は、スペックルの低減を高速に実行する超音波診断装置
を提供することにある。
The present invention has been made in view of the above points, and an object of the present invention is to provide an ultrasonic diagnostic apparatus that executes speckle reduction at high speed.

(課題を解決するための手段) 前記の課題を解決する本発明の超音波診断装置は、被
検体における一領域からの超音波におるエコー信号を異
なる複数の超音波ビームとして受信する受信手段と、前
記受信された複数の信号を対数増幅する増幅手段と、前
記対数増幅された複数の信号についての絶対値を演算す
る演算手段と、前記演算により得られた複数の絶対値を
加算する加算手段と、前記加算手段の出力信号を包絡線
検波する検波手段と、前記検波された信号に基づいて前
記被検体の画像を生成する画像生成手段とを備えたこと
を特徴とする。
(Means for Solving the Problems) An ultrasonic diagnostic apparatus of the present invention for solving the above-mentioned problems includes a receiving unit for receiving echo signals in an ultrasonic wave from one region in a subject as a plurality of different ultrasonic beams. Amplifying means for logarithmically amplifying the plurality of received signals, calculating means for calculating absolute values of the plurality of logarithmically amplified signals, and adding means for adding a plurality of absolute values obtained by the calculation A detection unit that performs envelope detection on an output signal of the addition unit; and an image generation unit that generates an image of the subject based on the detected signal.

(作用) アレイトランスデューサの複数個のエレメント群を主
音線用エレメント群として送信する。主音線用エレメン
ト群でエコー信号を受信すると共に、主音線用エレメン
ト群とは位置又は開口の異なる副音線用エレメント群で
も受信し、別個に増幅した後、それぞれの絶対値を求め
て加算し、スペックルノイズを低減する。
(Operation) A plurality of element groups of the array transducer are transmitted as main sound ray element groups. The echo signal is received by the main sound ray element group, and also the sub sound ray element group having a different position or opening from the main sound ray element group is received, amplified separately, and the absolute value of each is calculated and added. Reduce speckle noise.

(実施例) 以下、図面を参照して本発明の実施例を詳細に説明す
る。
(Example) Hereinafter, an example of the present invention will be described in detail with reference to the drawings.

第1図は本発明の方法の一実施例のアレイトランスデ
ューサの送受エレメントの配分を示す図である。図にお
いて、1は複数のエレメントが一列に配置されたN個の
エレメントを有するアレイトランスデューサ、2はi番
目のエレメントからi+n−1番目のエレメントまでn
個のエレメントにより送信し受信する主音線用エレメン
ト群、3はi+a番目からi+a+n−1番目までn個
のエレメントから成る受信専用の副音線用エレメント群
である。主音線用エレメント群2から送信し、被検体の
目標からのエコー信号を受信する。同時に副音線用エレ
メント群3により同じエコー信号を受信する。これが終
ればそれぞれ隣のエレメントに移って送受信を行う。
FIG. 1 is a diagram showing the distribution of transmitting and receiving elements of an array transducer according to one embodiment of the method of the present invention. In the figure, 1 is an array transducer having N elements in which a plurality of elements are arranged in a line, and 2 is n from the i-th element to the (i + n-1) -th element.
The main sound ray element group 3 transmitted and received by the number of elements is a reception-only sub sound ray element group composed of n elements from the (i + a) th element to the (i + a + n-1) th element. The echo signal is transmitted from the main sound ray element group 2 and received from the target of the subject. At the same time, the same echo signal is received by the sub sound ray element group 3. When this is completed, the process moves to the next element and performs transmission and reception.

第2図は第1図に示した、アレイトランスデューサ1
で受信した2信号を処理する回路のブロック図である。
図において、第1図と同等の部分には同一の符号を付し
てある。この図では主音線用エレメント群2は1番目か
らn−1番目のエレメント、副音線用エレメント群3は
2番目からn番目までエレメントで構成される例を示し
ている。図中、11は主音線用エレメント群2からの信号
(以下主音線信号という)を対数圧縮増幅する対数増幅
器、12は副音線用エレメント群3からの信号(以下副音
線信号という)を対数圧縮増幅する対数増幅器である。
13は対数増幅器11の出力の主音線信号の位相情報を除去
するため、その絶対値である振幅値を求めるための絶対
値演算回路、14は同じく副音線信号の振幅値を得るため
の絶対値演算回路である。15は絶対値演算回路13,14の
出力を加算する加算器、16は加算器15の出力信号を検波
して、その包絡線信号を出力する包絡線検波器、17は包
絡線検波器16で検波された信号を所要の信号処理をする
信号処理回路である。全エレメントのデータが集まった
段階で画像再構成回路18で画像再構成される。
FIG. 2 shows the array transducer 1 shown in FIG.
FIG. 3 is a block diagram of a circuit that processes two signals received by the control unit.
In the figure, the same parts as those in FIG. 1 are denoted by the same reference numerals. This figure shows an example in which the main sound ray element group 2 is composed of the first to n-1st elements, and the sub sound ray element group 3 is composed of the second to n-th elements. In the figure, reference numeral 11 denotes a logarithmic amplifier for logarithmically compressing and amplifying a signal from the main sound ray element group 2 (hereinafter referred to as main sound ray signal), and reference numeral 12 denotes a signal from the sub sound ray element group 3 (hereinafter referred to as sub sound ray signal). It is a logarithmic amplifier that performs logarithmic compression amplification.
13 is an absolute value calculation circuit for obtaining an amplitude value which is an absolute value of the phase information of the main sound ray signal output from the logarithmic amplifier 11, and 14 is also an absolute value calculation circuit for obtaining an amplitude value of the sub sound ray signal. It is a value calculation circuit. 15 is an adder that adds the outputs of the absolute value calculation circuits 13 and 14, 16 is an envelope detector that detects the output signal of the adder 15 and outputs its envelope signal, and 17 is an envelope detector 16. This is a signal processing circuit that performs required signal processing on the detected signal. When the data of all the elements is collected, the image is reconstructed by the image reconstruction circuit 18.

次に上記のように構成された装置の動作を第3図のフ
ローチャートを参照して説明することにより、本実施例
の方法を説明する。
Next, the operation of the apparatus configured as described above will be described with reference to the flowchart of FIG. 3 to explain the method of the present embodiment.

ステップ1 アレイトランスデューサ1からi番目〜(i+n−
1)番目のn個のエレメント群を送信用に選択する。
Step 1 i-th to (i + n-)
1) The n-th element group is selected for transmission.

ステップ2 i=1として送信用エレメント群を1番目からn番目
のn個のエレメントとする。
Step 2 Assuming that i = 1, the transmission element group is set to the first to n-th n elements.

ステップ3 n個の送信用エレメント群から送信する。Step 3 Transmission is performed from n transmission element groups.

(i〜i+n−1) ステップ4 送信用エレメント群と同じ主音線用エレメント群2
(i〜i+n−1)で受信する。(i+1〜i+n)番
目の副音線用エレメント群3でも受信する。
(I to i + n-1) Step 4 Same main sound ray element group 2 as the transmission element group
(I to i + n-1). The (i + 1-i + n) -th sub sound ray element group 3 is also received.

ステップ5 主音線信号を対数増幅器11で対数圧縮増幅して受信信
号の広いダイナミックレンジを画像表示に適する範囲の
レベルに圧縮する。副音線信号を対数増幅器12で対数圧
縮増幅する。
Step 5 The main sound ray signal is logarithmically compressed and amplified by the logarithmic amplifier 11 to compress a wide dynamic range of the received signal to a level suitable for image display. The auxiliary sound ray signal is logarithmically compressed and amplified by the logarithmic amplifier 12.

ステップ6 対数圧縮された主音線信号及び副音線信号から位相情
報を消去して位相差による信号間の相違をなくし、振幅
の最大値のみを得るために、主音線信号を絶対値演算回
路13において絶対値演算を行って絶対値である振幅値を
算出する。副音線信号を絶対値演算回路14において絶対
値演算を行って振幅値を算出する。
Step 6 In order to eliminate the phase information from the logarithmically compressed main sound ray signal and sub sound ray signal to eliminate the difference between the signals due to the phase difference and obtain only the maximum value of the amplitude, the main sound ray signal is subjected to an absolute value calculation circuit 13. , An absolute value calculation is performed to calculate an amplitude value which is an absolute value. The absolute value calculation circuit 14 performs an absolute value calculation on the auxiliary sound ray signal to calculate an amplitude value.

ステップ7 絶対値演算回路13と絶対値演算回路14との出力を加算
器15により加算し、主音線信号と副音線信号の絶対値の
和を求める。
Step 7 The outputs of the absolute value calculation circuit 13 and the absolute value calculation circuit 14 are added by the adder 15 to obtain the sum of the absolute values of the main sound ray signal and the sub sound ray signal.

ステップ8 加算器15の出力を包絡線検波器16により検波する。Step 8 The output of the adder 15 is detected by the envelope detector 16.

ステップ9 検波された信号を通常の超音波診断装置の信号処理回
路17において、AD変換,記憶装置に格納等の信号処理を
行う。
Step 9 The detected signal is subjected to signal processing such as AD conversion and storage in a storage device in a signal processing circuit 17 of a normal ultrasonic diagnostic apparatus.

ステップ10 アレイトランスデューサ1の全エレメントについて送
受信が終ったかをチェックする。終っていなければ、ス
テップ11に進む。終っていれば、ステップ12に進む。
Step 10 It is checked whether transmission / reception has been completed for all elements of the array transducer 1. If not, go to step 11. If so, go to step 12.

ステップ11 エレメント番号iに1を加えて、次のエレメント群を
選択し、ステップ3に戻る。
Step 11 One is added to the element number i, the next element group is selected, and the process returns to step 3.

ステップ12 全エレメントの送受信が終了した段階で補間を行いB
モード画像の再構成を行う。
Step 12 Interpolation is performed when all elements have been transmitted and received.
Reconstruct a mode image.

上記に説明したように、同一の送信信号に対して、2
組の伝播経路の異なる受信信号を同時に受信することに
より、スペックルについて相関の小さい信号を実時間
で、同一領域で得ることができ、相関の小さい信号を重
ね合わせることで、実時間で、イメージ領域を減少させ
ることなく、スペックルの低減したBモード画像が得ら
れるようになる。
As described above, for the same transmission signal, 2
By simultaneously receiving received signals having different sets of propagation paths, a signal having a small correlation with respect to speckles can be obtained in the same region in real time. A B-mode image with reduced speckles can be obtained without reducing the area.

尚、本発明は上記実施例に限定されるものではない。 The present invention is not limited to the above embodiment.

実施例では副音線用エレメント群3を主音線用エレメ
ント群2の隣りのエレメントから選んだが、適当数離隔
するようにしてもよい。但し、それがためにスペックル
低減領域を狹めないように留意する必要がある。
In the embodiment, the sub sound ray element group 3 is selected from the elements adjacent to the main sound ray element group 2, but may be separated by an appropriate number. However, care must be taken not to narrow the speckle reduction region because of that.

又、実施例では音響経路を変化させた2音線によりス
ペックル低減をしていたが、開口の大きさを変えた2音
線によるようにしてもよい。その場合のエレメント群の
選択は例えば次のようにして選ぶ。
Further, in the embodiment, speckle reduction is performed by using two sound rays whose acoustic paths are changed. However, two sound rays whose opening size is changed may be used. In this case, the element group is selected as follows, for example.

主音線用エレメント群…… i+a番目〜i+a+n−1番目のn個 副音線用エレメント群…… i番目〜i+2a+n−1番目のn+2a個 ここで、i=1〜N−(2a+n−1) この実施例では副音線のエレメントの数を主音線の両
側に等しい数だけ増加させた場合を示したが、左右の増
加数を等しくする必要はなく、異なる数を選んでも差し
支えない。
Primary sound ray element group ... i + a-th to i + a + n-1 n-th auxiliary sound ray element group ... i-th to i + 2a + n-1-th n + 2a where i = 1 to N- (2a + n-1) In the embodiment, the case where the number of elements of the auxiliary sound ray is increased by the same number on both sides of the main sound ray has been described. However, it is not necessary to make the left and right increase numbers equal, and a different number may be selected.

更に、主音線1個に対し副音線を複数個選んで行うこ
ともできる。
Further, it is also possible to select a plurality of sub sound rays for one main sound ray.

(発明の効果) 以上詳細に説明したように本発明によればスペックル
の低減を高速に実行することができる。
(Effects of the Invention) As described in detail above, according to the present invention, speckle reduction can be performed at high speed.

【図面の簡単な説明】[Brief description of the drawings]

第1図は本発明の方法の一実施例の送受信方法を示す
図、 第2図は第1図の方法で得られた信号の処理回路を示す
図、 第3図は本発明の方法のフローチャート、 第4図は従来の空間コンパウンド法によるスペックル低
減領域を示す図である。 1……アレイトランスデューサ 2……主音線用エレメント群 3……副音線用エレメント群 11,12……対数増幅器 13,14……絶対値演算回路 15……加算器 16,21,22……包絡線検波器 17……信号処理回路、18……画像再構成回路
FIG. 1 is a diagram showing a transmission / reception method according to an embodiment of the method of the present invention, FIG. 2 is a diagram showing a circuit for processing a signal obtained by the method of FIG. 1, and FIG. 3 is a flowchart of the method of the present invention. FIG. 4 is a diagram showing a speckle reduction region by the conventional spatial compound method. 1 Array transducer 2 Primary sound ray element group 3 Secondary sound ray element group 11,12 Logarithmic amplifier 13,14 Absolute value calculation circuit 15 Adder 16,21,22 Envelope detector 17: Signal processing circuit, 18: Image reconstruction circuit

───────────────────────────────────────────────────── フロントページの続き (58)調査した分野(Int.Cl.6,DB名) A61B 8/14 G01N 29/06,29/22──────────────────────────────────────────────────続 き Continued on the front page (58) Field surveyed (Int.Cl. 6 , DB name) A61B 8/14 G01N 29 / 06,29 / 22

Claims (1)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】被検体における一領域からの超音波による
エコー信号を異なる複数の超音波ビームとして受信する
受信手段と、 前記受信された複数の信号を対数増幅する増幅手段と、 前記対数増幅された複数の信号についての絶対値を演算
する演算手段と、 前記演算により得られた複数の絶対値を加算する加算手
段と、 前記加算手段の出力信号を包絡線検波する検波手段と、 前記検波された信号に基づいて前記被検体の画像を生成
する画像生成手段とを備えたことを特徴とする超音波診
断装置。
A receiving means for receiving echo signals from an area of the subject as a plurality of ultrasonic beams; an amplifying means for logarithmically amplifying the received signals; and a logarithmically amplified signal. Calculating means for calculating absolute values of the plurality of signals, adding means for adding the plurality of absolute values obtained by the calculation, detecting means for performing envelope detection on an output signal of the adding means, And an image generating means for generating an image of the subject based on the received signal.
JP1241327A 1989-09-18 1989-09-18 Ultrasound diagnostic equipment Expired - Lifetime JP2824673B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP1241327A JP2824673B2 (en) 1989-09-18 1989-09-18 Ultrasound diagnostic equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP1241327A JP2824673B2 (en) 1989-09-18 1989-09-18 Ultrasound diagnostic equipment

Publications (2)

Publication Number Publication Date
JPH03103246A JPH03103246A (en) 1991-04-30
JP2824673B2 true JP2824673B2 (en) 1998-11-11

Family

ID=17072649

Family Applications (1)

Application Number Title Priority Date Filing Date
JP1241327A Expired - Lifetime JP2824673B2 (en) 1989-09-18 1989-09-18 Ultrasound diagnostic equipment

Country Status (1)

Country Link
JP (1) JP2824673B2 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4590256B2 (en) * 2004-05-20 2010-12-01 富士フイルム株式会社 Ultrasonic imaging apparatus, ultrasonic image processing method, and ultrasonic image processing program
JP4644621B2 (en) * 2006-03-29 2011-03-02 株式会社日立エンジニアリング・アンド・サービス Ultrasonic inspection method and ultrasonic inspection apparatus
JP4862471B2 (en) * 2006-04-24 2012-01-25 Jfeスチール株式会社 Bath roll speed detecting device, bath roll speed control method and apparatus

Also Published As

Publication number Publication date
JPH03103246A (en) 1991-04-30

Similar Documents

Publication Publication Date Title
JP2777197B2 (en) Ultrasound diagnostic equipment
US5172343A (en) Aberration correction using beam data from a phased array ultrasonic scanner
Trahey et al. Experimental results with a real-time adaptive ultrasonic imaging system for viewing through distorting media
JP4620261B2 (en) Ultrasound imaging with higher order nonlinear components
US11372094B2 (en) Reverberation artifact cancellation in ultrasonic diagnostic images
KR20080039446A (en) Ultrasound imaging system and method for flow imaging using real-time spatial compounding
US11650300B2 (en) Ultrasound system and method for suppressing noise using per-channel weighting
JPS6134333B2 (en)
US7291108B2 (en) Ultrasonic transmission/reception apparatus for generating an image based on ultrasonic echoes and vibro-acoustic sounds
US4679565A (en) Ultrasonic diagnostic apparatus using non-linear parameters of an organ
US5187687A (en) Production of images
JPH0467856A (en) Ultrasonic imager by bistatic transmission/reception
US6432054B1 (en) Medical ultrasonic imaging with adaptive synthesis and compounding
US5476098A (en) Partially coherent imaging for large-aperture phased arrays
JP2824673B2 (en) Ultrasound diagnostic equipment
Hazard et al. Effects of motion on a synthetic aperture beamformer for real-time 3D ultrasound
JPH0614934B2 (en) Ultrasonic diagnostic equipment
JP4445255B2 (en) Method and apparatus for ultrasonic speckle reduction using broadband frequency synthesis with harmonics generated in tissue
JP2007021192A (en) Ultrasonic diagnostic apparatus and ultrasonic transmission/reception condition optimizing program
JP3078569B2 (en) Ultrasound diagnostic equipment
JPS5869538A (en) Ultrasonic diagnostic apparatus
JP2831719B2 (en) Ultrasound diagnostic equipment
JPS6214294B2 (en)
JP2002301071A (en) Ultrasonic imaging method and apparatus
JP4154043B2 (en) Ultrasonic imaging device

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080911

Year of fee payment: 10

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090911

Year of fee payment: 11

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090911

Year of fee payment: 11

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090911

Year of fee payment: 11

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100911

Year of fee payment: 12

EXPY Cancellation because of completion of term
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100911

Year of fee payment: 12