JPH08154939A - Ultrasonic imaging display method and ultrasonic imaging apparatus - Google Patents
Ultrasonic imaging display method and ultrasonic imaging apparatusInfo
- Publication number
- JPH08154939A JPH08154939A JP6299766A JP29976694A JPH08154939A JP H08154939 A JPH08154939 A JP H08154939A JP 6299766 A JP6299766 A JP 6299766A JP 29976694 A JP29976694 A JP 29976694A JP H08154939 A JPH08154939 A JP H08154939A
- Authority
- JP
- Japan
- Prior art keywords
- data
- signal
- moving object
- ultrasonic imaging
- delay
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 title claims description 11
- 238000013507 mapping Methods 0.000 claims abstract description 8
- 238000005070 sampling Methods 0.000 claims abstract description 6
- 239000006185 dispersion Substances 0.000 claims description 35
- 230000003111 delayed effect Effects 0.000 abstract description 13
- 230000001186 cumulative effect Effects 0.000 abstract description 10
- 210000004369 blood Anatomy 0.000 abstract description 2
- 239000008280 blood Substances 0.000 abstract description 2
- 230000017531 blood circulation Effects 0.000 description 29
- 230000000541 pulsatile effect Effects 0.000 description 16
- 238000010586 diagram Methods 0.000 description 11
- 238000012935 Averaging Methods 0.000 description 8
- 230000005540 biological transmission Effects 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 5
- 239000003086 colorant Substances 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 239000000523 sample Substances 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 210000000601 blood cell Anatomy 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 210000003743 erythrocyte Anatomy 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0609—Display arrangements, e.g. colour displays
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Description
【0001】[0001]
【産業上の利用分野】本発明は超音波イメージング表示
方法及び超音波イメージング装置に関し、特に、血流に
関して拍動流か定常流かを識別可能な状態で表示するこ
とが可能な超音波イメージング表示方法及び超音波イメ
ージング装置に関する。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an ultrasonic imaging display method and an ultrasonic imaging apparatus, and more particularly to an ultrasonic imaging display capable of displaying a pulsatile flow or a steady flow with respect to a blood flow in a distinguishable state. A method and an ultrasonic imaging apparatus.
【0002】[0002]
【従来の技術】超音波を被検体に照射すると、生体組織
を媒体として超音波が伝達されるが、臓器等の組織や病
変部のような周囲の組織との音響インピーダンスの差の
ある所から反射されて、また、血液における赤血球など
からは散乱によって、送波した超音波の一部が戻ってく
る。この反射体や散乱体が視線方向に運動または移動す
る物体であった場合、その反射波の周波数はドプラ効果
によって送信周波数から偏移する。超音波診断装置はこ
の周波数偏移量を測定して運動物体の速度及び移動方向
を表示観察して診断の用に供する装置である。2. Description of the Related Art When an ultrasonic wave is applied to a subject, the ultrasonic wave is transmitted through a living tissue as a medium. However, there is a difference in acoustic impedance from tissues such as organs and surrounding tissues such as lesions. A part of the transmitted ultrasonic wave returns by being reflected and scattered by red blood cells in blood. When the reflector or scatterer is an object that moves or moves in the direction of the line of sight, the frequency of the reflected wave deviates from the transmission frequency due to the Doppler effect. The ultrasonic diagnostic apparatus is an apparatus for measuring the frequency shift amount, displaying and observing the velocity and moving direction of the moving object, and providing the diagnostic object.
【0003】反射体が探触子の方に向かって動いている
時は受信周波数は送信周波数より高く、逆に遠ざかる時
は低くなる。そしてその偏移周波数は反射体の運動速度
に比例する。このドプラ効果を利用して、例えば、心臓
や血管内を流れる血流の方向と速度を知ることができ
る。When the reflector is moving toward the probe, the reception frequency is higher than the transmission frequency, and conversely, the reception frequency is low. The shift frequency is proportional to the moving speed of the reflector. Using this Doppler effect, for example, the direction and velocity of the blood flow in the heart or blood vessel can be known.
【0004】この超音波診断装置において、運動物体の
移動方向を例えば近付く場合に赤,遠ざかる場合に青の
ように色で表示し、その移動速度を輝度で表示するCF
M(Color Flow Mapping)モードを有する装置がある。In this ultrasonic diagnostic apparatus, the moving direction of a moving object is displayed in colors such as red when approaching and blue when moving away, and the moving speed is displayed in brightness.
There is a device having an M (Color Flow Mapping) mode.
【0005】ところで、超音波パルスの反射には固定物
体によるものと運動物体によるものとが混在しており、
特に運動物体に注目する時は固定物体による反射信号を
消去して運動物体の反射信号のみを表示したほうが分り
易いので、その目的のために、運動物体の反射信号のみ
を取り出すMTI(Moving Target Indicator) が用いら
れている。By the way, the reflection of the ultrasonic pulse includes both a fixed object and a moving object.
Especially when paying attention to a moving object, it is easier to understand by removing the reflection signal from the fixed object and displaying only the reflection signal of the moving object. For that purpose, MTI (Moving Target Indicator) that extracts only the reflection signal of the moving object. ) Is used.
【0006】しかし、固定した反射物体に関する情報も
医用超音波診断装置では重要な意味も持っている。それ
は、上記の情報が無いと表示画面中のどこが何を表わし
ているのか特定できないし、また、単に運動物体を表示
しただけではそれが何であるかを判別することが困難で
ある。However, the information about the fixed reflecting object also has important meaning in the medical ultrasonic diagnostic apparatus. It is not possible to specify what is in the display screen and what it does without the above information, and it is difficult to determine what it is by simply displaying a moving object.
【0007】このため、固定物体の表示と運動物体の表
示とを同一画面上に共通に表示し、相互に識別するため
に、血球からの散乱エコーから血流の方向,平均速度と
速度分散値及びパワー値などを演算し、色でそれらの諸
量をリアルタイムでBモード画面上に2次元表示してい
る。このモードをCFMモードといい、これらの機能を
有する装置がCFM超音波診断装置である。Therefore, in order to display the display of a fixed object and the display of a moving object in common on the same screen and distinguish them from each other, the direction of the blood flow, the average velocity and the velocity dispersion value from the scattered echoes from the blood cells are displayed. Also, the power value and the like are calculated, and various quantities thereof are two-dimensionally displayed on the B mode screen in real time by color. This mode is called a CFM mode, and a device having these functions is a CFM ultrasonic diagnostic device.
【0008】図6に従来のCFMモードを有する超音波
診断装置の構成図を示す。この図6の構成において、超
音波探触子1は送波電気信号を超音波に変換して被検体
内に送波し、被検体内から反射されて戻って来た超音波
信号を電気信号に変換する。送受信部2は送信信号を増
幅して超音波探触子1に送り、受波された受信信号を直
交検波等により復調した後に各処理部に送る。すなわ
ち、送受信部2の出力(復調出力(i,q信号等))は
Bモード処理部3及びCFM演算部4にそれぞれ供給さ
れる。FIG. 6 shows a block diagram of a conventional ultrasonic diagnostic apparatus having a CFM mode. In the configuration of FIG. 6, the ultrasonic probe 1 converts the transmitted electric signal into an ultrasonic wave and transmits the ultrasonic wave into the subject, and the ultrasonic signal reflected and returned from the inside of the subject is an electric signal. Convert to. The transmission / reception unit 2 amplifies the transmission signal and sends it to the ultrasonic probe 1, demodulates the received signal received by quadrature detection or the like, and then sends it to each processing unit. That is, the output (demodulation output (i, q signals, etc.)) of the transmission / reception unit 2 is supplied to the B-mode processing unit 3 and the CFM calculation unit 4, respectively.
【0009】Bモード処理部3では増幅,対数圧縮,検
波等を行ってBモード表示のための信号を出力して、画
像処理部5に供給する。また、送受信部2の復調出力は
CFM演算部4に入力される。このCFM演算部4で
は、MTIフィルタ等により移動目標のみの反射波を抽
出し、例えば血流からの反射波を抽出して演算を行い血
流速度を検出する。また、血流速度のばらつきにより、
分散(速度分散)を検出する。The B-mode processing unit 3 performs amplification, logarithmic compression, detection, etc., outputs a signal for B-mode display, and supplies it to the image processing unit 5. The demodulated output of the transmitter / receiver 2 is input to the CFM calculator 4. The CFM calculation unit 4 extracts the reflected wave of only the moving target by using the MTI filter or the like, and extracts the reflected wave from the blood flow, for example, to calculate the blood flow velocity. Also, due to variations in blood flow rate,
Detect dispersion (velocity dispersion).
【0010】これらBモード処理部3及びCFM演算部
4のそれぞれの出力は画像処理部5に供給され、Bモー
ド表示上に運動物体(血流速度,速度分散)に応じて着
色されたイメージを生成し、このイメージの表示をCR
T表示部7に行う。The respective outputs of the B mode processing unit 3 and the CFM calculation unit 4 are supplied to the image processing unit 5, and an image colored according to a moving object (blood flow velocity, velocity dispersion) is displayed on the B mode display. Generate and CR display this image
Perform on the T display unit 7.
【0011】例えば、移動方向で色(色相)を変え(赤
/青)、移動速度に応じて赤/青の輝度を変えるように
して表示を行う。また、速度分散に応じて緑を混ぜるよ
うにした表示を行う。通常、赤は近づく流れを、青は遠
ざかる流れを表し、分散の大きさは緑の輝度の違いよっ
て表している。従って、部位によって、青と緑との混合
が行われ、また、赤と緑との混合が行われる。この場
合、赤,青,緑は光の三原色であるので、混合された色
の変化により、方向,速度,分散の大きさを読みとれ
る。For example, the display is performed by changing the color (hue) in the moving direction (red / blue) and changing the brightness of red / blue according to the moving speed. In addition, the display is made so that green is mixed according to the velocity dispersion. Usually, red represents a flowing flow, blue represents a flowing flow, and the magnitude of dispersion is represented by the difference in the brightness of green. Therefore, depending on the site, blue and green are mixed, and red and green are mixed. In this case, since red, blue, and green are the three primary colors of light, the direction, velocity, and magnitude of dispersion can be read by changing the mixed colors.
【0012】図7にCFM演算部4の内部の回路構成の
一例を示す。この図7において、前述の図6の送受信部
2の復調出力が位相が直交する信号(i信号,q信号)
であった場合、i,q信号それぞれがMTIフィルタ4
00,401によりフィルタリングされて、運動物体か
らの信号成分のみの信号が抽出される。FIG. 7 shows an example of the internal circuit configuration of the CFM calculator 4. In FIG. 7, the demodulated outputs of the transmitter / receiver 2 of FIG. 6 described above are signals whose phases are orthogonal (i signal, q signal)
, The i and q signals are respectively MTI filter 4
The signal of only the signal component from the moving object is extracted after being filtered by 00 and 401.
【0013】MTIフィルタ400,401を通過後の
運動物体からの信号をそれぞれI,Q信号とする。ここ
で、I,Q信号それぞれは遅延部402,403に供給
される。この遅延部402,403の遅延時間はΔtで
あるとする。また、この遅延時間Δtは、送受信部2で
同一音線の信号を受信する間隔(若しくはその整数倍)
の間隔とする。そして、遅延部402,403により遅
延された信号と遅延されていない信号とがそれぞれ掛算
器405〜408で掛算され、更に加算器409と減算
器410で加減算される。The signals from the moving object after passing through the MTI filters 400 and 401 are referred to as I and Q signals, respectively. Here, the I and Q signals are supplied to the delay units 402 and 403, respectively. The delay time of the delay units 402 and 403 is Δt. Further, this delay time Δt is an interval (or an integral multiple thereof) at which the signal of the same sound ray is received by the transmitter / receiver 2.
Interval. Then, the signals delayed by the delay units 402 and 403 and the signals not delayed are respectively multiplied by the multipliers 405 to 408, and further added and subtracted by the adder 409 and the subtractor 410.
【0014】ここで、I信号として、I1 ,I2 ,I
3 ,I4 とΔt毎の4信号を考える。同様にして、Q信
号として、Q1 ,Q2 ,Q3 ,Q4 のΔt毎の4信号を
考える。このような信号が与えられた場合、加算器40
9の出力にはΔt毎に、I1I 2 +Q1 Q2 ,I2 I3
+Q2 Q3 ,I3 I4 +Q3 Q4 が得られる。この信
号をDENと呼ぶことにする。Here, as the I signal, I1, I2, I
Consider 3 signals, I 4 and 4 signals for each Δt. Similarly, four Q signals Q1, Q2, Q3, and Q4 for each Δt are considered. When such a signal is given, the adder 40
The output of 9 is I1I2 + Q1Q2, I2I3 for each Δt.
+ Q2 Q3, I3 I4 + Q3 Q4 are obtained. This signal will be called DEN.
【0015】また、同様に減算器410の出力にはΔt
毎に、I1 Q2 −Q1 I2 ,I2 Q3 −Q2 I3 ,I3
Q4 −Q3 I4 が得られる。この信号をNUMと呼ぶこ
とにする。Similarly, the output of the subtractor 410 has Δt
I1 Q2 -Q1 I2, I2 Q3 -Q2 I3, I3
Q4 -Q3 I4 is obtained. This signal will be called NUM.
【0016】また、パワー算出部404では、2乗加算
処理が行なわれて、Δt毎に、I12 +Q1 2 ,I2 2 +
Q2 2 ,I3 2 +Q3 2 ,I4 2 +Q4 2 が得られる。
この信号をパワー値と呼び、Pと略することにする。Further, the power calculation unit 404 adds squares.
Processing is performed and I12 + Q12, I22+
Q22, I32+ Q32, I42+ Q42Is obtained.
This signal is called a power value and is abbreviated as P.
【0017】そして、遅延加算器411にて、Δt毎の
遅延加算(累加平均)が行なわれて、((I1 I2 +Q
1 Q2 )+(I2 I3 +Q2 Q3 )+(I3 I4 +Q3
Q4))/3が得られる。これをここでは、*DEN*
として略する。Then, the delay adder 411 performs delay addition (cumulative averaging) for each Δt to obtain ((I1 I2 + Q
1 Q2) + (I2 I3 + Q2 Q3) + (I3 I4 + Q3
Q4)) / 3 is obtained. Here it is * DEN *
Abbreviated as.
【0018】同様にして、遅延加算器412にて、Δt
毎の遅延加算(累加平均)が行なわれて、((I1 Q2
−Q1 I2 )+(I2 Q3 −Q2 I3 )+(I3 Q4 −
Q3I4 ))/3が得られる。同様にこれを*NUM*
と略する。Similarly, in the delay adder 412, Δt
Delay addition (cumulative averaging) is performed for each ((I1 Q2
-Q1 I2) + (I2 Q3 -Q2 I3) + (I3 Q4-
Q3I4)) / 3 is obtained. Similarly, set this to * NUM *
Abbreviated.
【0019】また、遅延加算器413にて、Δt毎の遅
延加算(累加平均)が行なわれて、((I1 2 +Q1
2 )+(I2 2 +Q2 2 )+(I3 2 +Q3 2 )+
(I4 2 +Q4 2 ))/4が得られる。これを*P*と
略する。Further, the delay adder 413 performs delay addition (cumulative averaging) for each Δt to obtain ((I1 2 + Q1
2 ) + (I2 2 + Q2 2 ) + (I3 2 + Q3 2 ) +
(I4 2 + Q4 2 )) / 4 is obtained. This is abbreviated as * P *.
【0020】このような*DEN*,*NUM*,*P
*が得られると、平均速度算出部414において、ta
n-1((*NUM*)/(*DEN*))という処理が
実行されて、平均速度*v*が求められる。Such * DEN *, * NUM *, * P
When * is obtained, the average speed calculation unit 414 indicates ta
The process of n -1 ((* NUM *) / (* DEN *)) is executed to obtain the average speed * v *.
【0021】また、分散算出部415において、(1−
((*DEN*2 +*NUM*2 ) 1/2)/*P*)と
いう処理が実行されて、分散*var*が求められる。
このようにして得られた平均速度と分散とを受けた画像
処理部5が、それぞれの信号値に応じて赤/青,緑など
の輝度を変えてカラー表示を行う。Further, in the variance calculation unit 415, (1-
((* DEN *2+ * NUM *2) 1/2) / * P *)
Such processing is executed to obtain the variance * var *.
An image that has been subjected to the average velocity and variance thus obtained
The processing unit 5 uses red / blue, green, etc. according to each signal value.
The color is displayed by changing the brightness of.
【0022】以上のような平均速度,分散と血流との関
係を図8に示す。図8(イ)は血流速度を示しており、
横軸は時間tを示し、縦軸は血流速度vを示している。
尚、ハッチングで示した部分はバラツキを意味してい
る。ここでは、1心拍に相当する時間を示しており、こ
こでは時相I,II,III,IV,Vのそれぞれのタ
イミングで上述の平均速度と分散とを求める処理を実行
している。FIG. 8 shows the relationship between the average velocity and dispersion and the blood flow as described above. FIG. 8A shows the blood flow velocity,
The horizontal axis represents time t, and the vertical axis represents blood flow velocity v.
Incidentally, the hatched portions mean variations. Here, the time corresponding to one heartbeat is shown, and here, the processing for obtaining the above-mentioned average velocity and variance is executed at each timing of the time phases I, II, III, IV, and V.
【0023】従って、各時相I〜Vに応じて平均速度と
分散とが求められ、それぞれ図8(ロ)〜(ヘ)のよう
な速度スペクトラムになる。この各図では、山のピーク
の位置が平均速度、山の広がりが分散を示し、縦軸は強
度を示すスペクトラムを示している。Therefore, the average velocity and the dispersion are obtained according to the respective time phases I to V, and the velocity spectra as shown in FIGS. 8B to 8F are obtained. In each of these figures, the position of the peak of the mountain shows the average velocity, the spread of the mountain shows the dispersion, and the vertical axis shows the spectrum showing the intensity.
【0024】尚、ここでは、各時相の平均速度と分散と
をグラフで示したが、実際には画像処理部5の処理によ
りカラー表示が行なわれる。従って、実際の装置では、
図8(ロ)〜(ヘ)に応じて、平均速度と分散とのカラ
ー表示が刻々と変換するように示される。Here, although the average speed and the dispersion of each time phase are shown in the graph, the color display is actually performed by the processing of the image processing unit 5. So in a real device,
According to FIGS. 8B to 8F, the color display of the average velocity and the variance is shown to be changed every moment.
【0025】[0025]
【発明が解決しようとする課題】以上のようにして平均
速度と分散とを表示する超音波診断装置は、異なった時
相毎に血流の平均速度と分散とを表示するものである。The ultrasonic diagnostic apparatus which displays the average velocity and the dispersion as described above displays the average velocity and the dispersion of the blood flow for each different time phase.
【0026】従って、時間的な拍動性を表示することは
できず、血流が拍動流(血流速度や分散が周期的に変化
する流れ)なのか定常流(血流速度や分散がほぼ一定な
ような流れ)なのかを瞬時に識別することは困難であっ
た。Therefore, it is not possible to display temporal pulsatility, and whether the blood flow is a pulsatile flow (a flow in which the blood flow velocity or dispersion changes periodically) or a steady flow (the blood flow velocity or dispersion is It was difficult to instantly identify whether the flow was almost constant).
【0027】このような場合に、拍動性を識別可能なよ
うに表示するためには、パルスドプラなどによりFFT
で周波数分析を行なって検証する必要が有り、装置が複
雑化する欠点がある。In such a case, in order to display the pulsatility so that the pulsatility can be identified, it is possible to use FFT by pulse Doppler or the like.
It is necessary to perform a frequency analysis in order to verify, and there is a drawback that the device becomes complicated.
【0028】また、CFM演算部4の出力を長時間トレ
ースするような回路を設けることによっても識別が可能
であるが、大規模な回路付加が必要になる問題を有して
いる。Further, the discrimination can be made by providing a circuit for tracing the output of the CFM calculating section 4 for a long time, but there is a problem that a large scale circuit needs to be added.
【0029】以上のような問題のために、簡単な構成で
拍動性を識別可能なように表示することは極めて困難で
あった。本発明は上記の点に鑑みてなされたもので、第
1の目的は、血流に関して拍動流か定常流かを識別可能
な状態で表示することが可能な超音波イメージング表示
方法を実現することである。Due to the above problems, it was extremely difficult to display the pulsatility so as to be distinguishable with a simple structure. The present invention has been made in view of the above points, and a first object thereof is to realize an ultrasonic imaging display method capable of displaying a pulsatile flow or a steady flow with respect to a blood flow in a distinguishable state. That is.
【0030】また、第2の目的は、血流に関して拍動流
か定常流かを識別可能な状態で表示することが可能な超
音波イメージング装置を実現することである。A second object is to realize an ultrasonic imaging apparatus capable of displaying a pulsatile flow or a steady flow with respect to the blood flow in a distinguishable state.
【0031】[0031]
【課題を解決するための手段】前記の課題を解決する第
1の手段は、受波信号より得られた運動物体反射信号中
の同一音線データのサンプリング間隔Δtより大きい間
隔ΔT毎に所定個数の各運動物体反射信号の遅延加算を
行い、遅延加算された運動物体反射信号から流速データ
を求めてカラーフローマッピングによる超音波イメージ
ング表示を行うことを特徴とする超音波イメージング表
示方法である。A first means for solving the above-mentioned problems is a predetermined number for every interval ΔT which is larger than the sampling interval Δt of the same sound ray data in the moving object reflection signal obtained from the received wave signal. The ultrasonic imaging display method is characterized by performing delayed addition of each moving object reflection signal, and calculating flow velocity data from the delayed addition moving object reflection signal to perform ultrasonic imaging display by color flow mapping.
【0032】ここで、運動物体反射信号とは受波信号か
ら検出された運動物体による反射波成分の信号であり、
同一音線のデータはΔt間隔で得られる。また、Δtよ
り大きい間隔ΔTとは異なる時相となるような時間間隔
を意味する。また、流速データとは、血流等の速度等を
示すデータである。Here, the moving object reflection signal is a signal of a reflected wave component from the moving object detected from the received signal,
Data for the same sound ray is obtained at Δt intervals. Further, it means a time interval having a time phase different from the interval ΔT larger than Δt. The flow velocity data is data indicating the velocity of blood flow and the like.
【0033】前記の課題を解決する第2の手段は、受波
信号を分析してBモードデータを得るBモード処理部
と、受波信号を分析して運動物体反射信号を求め、同一
音線データのサンプリング間隔Δtより大きい間隔ΔT
毎に所定個数の運動物体反射信号の遅延加算を行って流
速データを求めるCFM演算部と、前記Bモード処理部
で得られたBモードデータと前記CFM演算部で遅延加
算されて得られた流速データとを用いてカラーフローマ
ッピングによる超音波イメージング表示用のイメージデ
ータを生成する画像処理部とを備えたことを特徴とする
超音波イメージング装置である。The second means for solving the above-mentioned problems is to analyze the received signal to obtain B-mode data, and a B-mode processing unit to analyze the received signal to obtain a moving object reflection signal, and to obtain the same sound line. An interval ΔT larger than the data sampling interval Δt
A CFM calculation unit for delay-adding a predetermined number of moving object reflection signals for each time to obtain flow velocity data, a B mode data obtained by the B mode processing unit, and a flow velocity obtained by delay addition by the CFM calculation unit. And an image processing unit for generating image data for ultrasonic imaging display by color flow mapping using the data and the image processing unit.
【0034】ここで、流速データとは受波信号から得ら
れる血流等の速度や分散等を示すデータのことである。
前記の課題を解決する第3の手段は、受波信号を分析し
てBモードデータを得るBモード処理部と、受波信号を
分析して運動物体反射信号を求め、同一音線データのサ
ンプリング間隔Δtより大きい間隔ΔT毎に所定個数の
運動物体反射信号の遅延加算を行って速度データ及び分
散データを求めるCFM演算部と、前記Bモード処理部
で得られたBモードデータ,前記CFM演算部で遅延加
算されて得られた速度データ及び分散データを用いてカ
ラーフローマッピングによる超音波イメージング表示用
のイメージデータを生成する画像処理部とを備えたこと
を特徴とする超音波イメージング装置である。Here, the flow velocity data is data indicating the velocity and dispersion of the blood flow obtained from the received wave signal.
A third means for solving the above problems is a B-mode processing unit that analyzes a received signal to obtain B-mode data, and a received object signal is analyzed to obtain a moving object reflection signal, and the same sound ray data is sampled. A CFM calculating unit for delay-adding a predetermined number of moving object reflection signals at intervals ΔT larger than the interval Δt to obtain velocity data and dispersion data, B-mode data obtained by the B-mode processing unit, and the CFM calculating unit And an image processing unit for generating image data for ultrasonic imaging display by color flow mapping using the velocity data and the dispersion data obtained by delaying and adding in the ultrasonic imaging apparatus.
【0035】ここで、速度データは血流等の速度を示
し、分散データは血流等の速度のばらつきを示してい
る。Here, the velocity data indicates the velocity of blood flow and the like, and the dispersion data indicates the variation of velocity such as blood flow.
【0036】[0036]
【作用】課題を解決する第1の手段である超音波イメー
ジング表示方法において、Δt間隔で得られた運動物体
反射信号から流速データが求められるが、この際に、Δ
tより大きい間隔であって異なる時相のΔT毎の運動物
体反射信号が遅延加算されたものにより流速データが求
められる。そして、このようにして求められた流速デー
タによりCFMの超音波イメージング表示がなされる。
このようにすることで、拍動流か定常流かを識別可能な
状態で表示することが可能になる。In the ultrasonic imaging display method which is the first means for solving the problem, the flow velocity data is obtained from the moving object reflection signal obtained at the Δt intervals.
The flow velocity data is obtained by delaying and adding the moving object reflection signals for each ΔT of different time phases at intervals larger than t. Then, the CFM ultrasonic imaging display is performed based on the flow velocity data thus obtained.
By doing so, it becomes possible to display the pulsating flow and the steady flow in a distinguishable state.
【0037】課題を解決する第2の手段である超音波イ
メージング装置において、Δt間隔で得られた運動物体
反射信号から流速データが求められるが、この際に、Δ
tより大きい間隔であって異なる時相のΔT毎の運動物
体反射信号が遅延加算されたものにより流速データが求
められる。そして、このようにして求められた流速デー
タによりCFMの超音波イメージング表示がなされる。
このようにすることで、拍動流か定常流かを識別可能な
状態で表示することが可能になる。In the ultrasonic imaging apparatus which is the second means for solving the problem, the flow velocity data is obtained from the moving object reflection signal obtained at Δt intervals.
The flow velocity data is obtained by delaying and adding the moving object reflection signals for each ΔT of different time phases at intervals larger than t. Then, the CFM ultrasonic imaging display is performed based on the flow velocity data thus obtained.
By doing so, it becomes possible to display the pulsating flow and the steady flow in a distinguishable state.
【0038】課題を解決する第3の手段である超音波イ
メージング装置において、Δt間隔で得られた運動物体
反射信号から速度データ及び分散データが求められる
が、この際に、Δtより大きい間隔であって異なる時相
のΔT毎の運動物体反射信号が遅延加算されたものによ
り速度データ及び分散データが求められる。そして、こ
のようにして求められた速度データ及び分散データによ
りCFMの超音波イメージング表示がなされる。このよ
うにすることで、拍動流か定常流かを識別可能な状態で
表示することが可能になる。In the ultrasonic imaging apparatus which is the third means for solving the problem, the velocity data and the dispersion data are obtained from the moving object reflection signal obtained at the Δt intervals. At this time, the intervals are larger than Δt. Then, the velocity data and the dispersion data are obtained by the delayed addition of the moving object reflection signals for each ΔT of different time phases. Then, the CFM ultrasonic imaging display is performed based on the velocity data and the dispersion data thus obtained. By doing so, it becomes possible to display the pulsating flow and the steady flow in a distinguishable state.
【0039】[0039]
【実施例】以下、図面を参照して本発明の実施例を詳細
に説明する。図1は本発明の超音波イメージング表示方
法の一実施例を実現するための装置であり、また、本発
明の超音波イメージング装置の一実施例の概略構成を示
す構成図である。まず、この図1により超音波イメージ
ング装置の概要を説明する。Embodiments of the present invention will be described below in detail with reference to the drawings. 1 is an apparatus for realizing an embodiment of the ultrasonic imaging display method of the present invention, and is a configuration diagram showing a schematic configuration of an embodiment of the ultrasonic imaging apparatus of the present invention. First, the outline of the ultrasonic imaging apparatus will be described with reference to FIG.
【0040】この図1は前述のCFM演算部4の内部の
回路構成の一例を示す構成図である。この図1におい
て、前述の図6の送受信部2の復調出力が位相が直交す
る信号(I信号,Q信号)であった場合、I,Q信号そ
れぞれがMTIフィルタ400,401によりフィルタ
リングされて、運動物体からの信号成分のみの信号(運
動物体反射信号)が抽出される。FIG. 1 is a block diagram showing an example of the internal circuit configuration of the CFM calculator 4 described above. In FIG. 1, when the demodulated output of the transmission / reception unit 2 of FIG. 6 described above is a signal (I signal, Q signal) whose phases are orthogonal to each other, the I and Q signals are filtered by the MTI filters 400 and 401, respectively. A signal (moving object reflection signal) containing only the signal component from the moving object is extracted.
【0041】MTIフィルタ400,401を通過後の
運動物体反射信号をそれぞれI,Q信号とする。ここ
で、I,Q信号それぞれは遅延部402,403に供給
される。この遅延部402,403の遅延時間はΔtで
あるとする。また、この遅延時間Δtは、送受信部2で
同一音線の信号を受信する間隔(若しくはその整数倍)
の間隔とする。そして、遅延部402,403により遅
延された信号と遅延されていない信号とがそれぞれ掛算
器405〜408で掛算され、更に加算器409と減算
器410で加減算される。The moving object reflection signals after passing through the MTI filters 400 and 401 are respectively I and Q signals. Here, the I and Q signals are supplied to the delay units 402 and 403, respectively. The delay time of the delay units 402 and 403 is Δt. Further, this delay time Δt is an interval (or an integral multiple thereof) at which the signal of the same sound ray is received by the transmitter / receiver 2.
Interval. Then, the signals delayed by the delay units 402 and 403 and the signals not delayed are respectively multiplied by the multipliers 405 to 408, and further added and subtracted by the adder 409 and the subtractor 410.
【0042】ここで、I信号として、I1 ,I2 ,I
3 ,I4 とΔt毎の4信号を考える。同様にして、Q信
号として、Q1 ,Q2 ,Q3 ,Q4 のΔt毎の4信号を
考える。このような信号が与えられた場合、加算器40
9の出力にはΔt毎に、I1I 2 +Q1 Q2 ,I2 I3
+Q2 Q3 ,I3 I4 +Q3 Q4 が得られる。この信
号をDENと呼ぶことにする。Here, as the I signal, I1, I2, I
Consider 3 signals, I 4 and 4 signals for each Δt. Similarly, four Q signals Q1, Q2, Q3, and Q4 for each Δt are considered. When such a signal is given, the adder 40
The output of 9 is I1I2 + Q1Q2, I2I3 for each Δt.
+ Q2 Q3, I3 I4 + Q3 Q4 are obtained. This signal will be called DEN.
【0043】また、同様に減算器410の出力にはΔt
毎に、I1 Q2 −Q1 I2 ,I2 Q3 −Q2 I3 ,I3
Q4 −Q3 I4 が得られる。この信号をNUMと呼ぶこ
とにする。Similarly, the output of the subtractor 410 has Δt
I1 Q2 -Q1 I2, I2 Q3 -Q2 I3, I3
Q4 -Q3 I4 is obtained. This signal will be called NUM.
【0044】また、パワー算出部404では、2乗加算
処理が行なわれて、Δt毎に、I1 2 +Q1 2 ,I2 2
+Q2 2 ,I3 2 +Q3 2 ,I4 2 +Q4 2 が得られ
る。この信号をパワー値と呼び、Pと略することにす
る。Further, the power calculation unit 404 adds squares.
Processing is performed and I1 2+ Q12, I22
+ Q22, I32+ Q32, I42+ Q42Is obtained
It This signal is called a power value and is abbreviated as P.
It
【0045】そして、遅延加算器411にて、Δt毎の
遅延加算(累加平均)が行なわれて、((I1 I2 +Q
1 Q2 )+(I2 I3 +Q2 Q3 )+(I3 I4 +Q3
Q4))/3が得られる。同様にこれを*DEN*と略
する。そして、図8において説明した時相I〜Vに対応
して、各時相の*DEN*を、*DEN*I,*DEN
*II,*DEN*III,*DEN*IV,*DEN
*Vのように示す。Then, the delay adder 411 performs delay addition (cumulative averaging) for each Δt to obtain ((I1 I2 + Q
1 Q2) + (I2 I3 + Q2 Q3) + (I3 I4 + Q3
Q4)) / 3 is obtained. Similarly, this is abbreviated as * DEN *. Then, corresponding to the time phases I to V described in FIG. 8, the * DEN * of each time phase is changed to * DEN * I, * DEN.
* II, * DEN * III, * DEN * IV, * DEN
* Indicated as V.
【0046】同様にして、遅延加算器412にて、Δt
毎の遅延加算(累加平均)が行なわれて、((I1 Q2
−Q1 I2 )+(I2 Q3 −Q2 I3 )+(I3 Q4 −
Q3 I4 ))/3が得られる。これをここでは、*NU
M*として略する。Similarly, in the delay adder 412, Δt
Delay addition (cumulative averaging) is performed for each ((I1 Q2
-Q1 I2) + (I2 Q3 -Q2 I3) + (I3 Q4-
Q3 I4)) / 3 is obtained. This is here, * NU
Abbreviated as M *.
【0047】そして、図8において説明した時相I〜V
に対応して、各時相の*NUM*を、*NUM*I,*
NUM*II,*NUM*III,*NUM*IV,*
NUM*Vのように示す。Then, the time phases I to V explained in FIG.
Corresponding to, * NUM * of each time phase, * NUM * I, *
NUM * II, * NUM * III, * NUM * IV, *
It is shown as NUM * V.
【0048】また、遅延加算器413にて、Δt毎の遅
延加算(累加平均)が行なわれて、((I1 2 +Q1
2 )+(I2 2 +Q2 2 )+(I3 2 +Q3 2 )+
(I4 2 +Q4 2 ))/4が得られる。これを*P*と
略する。そして、図8において説明した時相I〜Vに対
応して、各時相の*P*を、*P*I,*P*II,*
P*III,*P*IV,*P*Vのように示す。The delay adder 413 performs delay addition (cumulative averaging) for each Δt to obtain ((I1 2 + Q1
2 ) + (I2 2 + Q2 2 ) + (I3 2 + Q3 2 ) +
(I4 2 + Q4 2 )) / 4 is obtained. This is abbreviated as * P *. Then, corresponding to the time phases I to V described in FIG. 8, the * P * of each time phase is changed to * P * I, * P * II, *.
It is shown as P * III, * P * IV, * P * V.
【0049】ここで、遅延加算器416にて、Δtより
大きい間隔ΔT(異なる時相となるような間隔:ここで
は、ΔTを時相I,II,III,IV,Vの間隔に合
わせた場合を示す)毎の遅延加算(累加平均)が行われ
る。この結果、((*DEN*I)+(*DEN*I
I)+(*DEN*III)+(*DEN*IV)+
(*DEN*V)/5)のような演算が行われる。この
演算結果を、*DEN*Tと呼ぶ。この*DEN*T
は、異なる時相I〜VのデータDENの平均値である。Here, in the delay adder 416, an interval ΔT larger than Δt (intervals having different time phases: here, when ΔT is adjusted to the intervals of the time phases I, II, III, IV, V). (Indicating) is performed for each delay addition (cumulative average). As a result, ((* DEN * I) + (* DEN * I
I) + (* DEN * III) + (* DEN * IV) +
An operation such as (* DEN * V) / 5) is performed. The result of this calculation is called * DEN * T. This * DEN * T
Is an average value of data DEN of different time phases I to V.
【0050】また同様に、遅延加算器417にて、Δt
より大きい間隔ΔT(異なる時相となるような間隔)毎
の遅延加算(累加平均)が行われる。この結果、((*
NUM*I)+(*NUM*II)+(*NUM*II
I)+(*NUM*IV)+(*NUM*V)/5)の
ような演算が行われる。この演算結果を、*NUM*T
と呼ぶ。この*NUM*Tは、異なる時相I〜Vのデー
タNUMの平均値である。Similarly, in the delay adder 417, Δt
Delay addition (cumulative averaging) is performed for each larger interval ΔT (intervals having different time phases). As a result,((*
NUM * I) + (* NUM * II) + (* NUM * II
I) + (* NUM * IV) + (* NUM * V) / 5) is performed. This calculation result is * NUM * T
Call. This * NUM * T is an average value of the data NUM of different time phases I to V.
【0051】更に同様に、遅延加算器418にて、Δt
より大きい間隔ΔT(異なる時相となるような間隔)毎
の遅延加算(累加平均)が行われる。この結果、((*
P*I)+(*P*II)+(*P*III)+(*P
*IV)+(*P*V)/5)のような演算が行われ
る。この演算結果を、*P*Tと呼ぶ。この*P*T
は、異なる時相I〜VのデータPの平均値である。Similarly, in the delay adder 418, Δt
Delay addition (cumulative averaging) is performed for each larger interval ΔT (intervals having different time phases). As a result,((*
P * I) + (* P * II) + (* P * III) + (* P
Calculation such as * IV) + (* P * V) / 5) is performed. The result of this calculation is called * P * T. This * P * T
Is an average value of data P of different time phases I to V.
【0052】このような*DEN*T,*NUM*T,
*P*Tが得られると、平均速度算出部414におい
て、tan-1((*NUM*T)/(*DEN*T))
という処理が実行されて、平均速度*v*Tが求められ
る。Such * DEN * T, * NUM * T,
When * P * T is obtained, the average speed calculation unit 414 calculates tan −1 ((* NUM * T) / (* DEN * T)).
Is executed to obtain the average speed * v * T.
【0053】また、分散算出部415において、(1−
((*DEN*T2 +*NUM*T 2 )1/2 )/*P*
T)という処理が実行されて、分散*var*Tが求め
られる。Further, in the variance calculation unit 415, (1-
((* DEN * T2+ * NUM * T 2)1/2) / * P *
T) is executed to obtain the variance * var * T
Can be
【0054】このようにして得られた平均速度と分散と
を受けた画像処理部5が、それぞれの信号値に応じて赤
/青の色相に輝度を変え、また、この赤/青に緑などの
輝度を変えてた表示を重畳してカラー表示を行う。The image processing unit 5, which receives the average speed and the variance thus obtained, changes the brightness to the hue of red / blue according to the respective signal values, and changes the brightness of the red / blue to green or the like. The color display is performed by superimposing the display in which the luminance of is changed.
【0055】従って、本実施例では、異なった時相(I
〜V)にわたる血流の平均速度と分散との平均値を表示
するものである。このため、時間的な拍動性を表示する
ことができるようになり、血流が拍動流(血流速度や分
散が周期的に変化する流れ)なのか定常流(血流速度や
分散がほぼ一定なような流れ)なのかを瞬時に識別する
ことが可能になる。Therefore, in this embodiment, different time phases (I
.About.V), the average value of the average velocity and variance of the blood flow is displayed. For this reason, it becomes possible to display temporal pulsatility, and whether the blood flow is a pulsatile flow (a blood flow velocity or dispersion changes periodically) or a steady flow (blood flow velocity or dispersion is It is possible to instantly identify whether or not the flow is almost constant).
【0056】尚、上述の実施例では図8に示した時相I
〜VのΔT間隔の5時相のデータを遅延加算したが、こ
れ以外の個数のデータを遅延加算するようにしても構わ
ない。In the above embodiment, the time phase I shown in FIG.
Although the data of 5 time phases at ΔT intervals of ˜V are delayed and added, other numbers of data may be delayed and added.
【0057】図2は拍動性の血流速度の時間変化を示し
ており、横軸は時間tを示し、縦軸は血流速度vを示し
ている。尚、ハッチングで示した部分はバラツキを意味
している。ここでは、1心拍に相当する時間を示してお
り、例えば上述の時相I,II,III,IV,Vに相
当するタイミングで上述の平均速度と分散とを求める処
理を実行している。FIG. 2 shows changes in pulsatile blood flow velocity with time. The horizontal axis represents time t and the vertical axis represents blood flow velocity v. Incidentally, the hatched portions mean variations. Here, the time corresponding to one heartbeat is shown, and for example, the above-described processing for obtaining the average velocity and the variance is executed at the timings corresponding to the above-mentioned time phases I, II, III, IV, and V.
【0058】図3は定常性の血流速度の時間変化を示し
ており、横軸は時間tを示し、縦軸は血流速度vを示し
ている。尚、ハッチングで示した部分はバラツキを意味
している。ここでは、図2と同じ時間を示しており、例
えば上述の時相I,II,III,IV,Vに相当する
タイミングで上述の平均速度と分散とを求める処理を実
行している。FIG. 3 shows the change over time of the steady blood flow velocity. The horizontal axis shows the time t and the vertical axis shows the blood flow velocity v. Incidentally, the hatched portions mean variations. Here, the same time as in FIG. 2 is shown, and for example, the processing for obtaining the above-mentioned average velocity and dispersion is executed at the timings corresponding to the above-mentioned time phases I, II, III, IV, V.
【0059】このような図2,図3のような血流速度に
ついて、上述の演算を実行した場合の平均速度及び分散
はそれぞれ図4,図5に示すようになる。すなわち、拍
動流では、平均速度が広い範囲に分布している。一方、
定常流では、平均速度がほぼ一定の値に集中している。With respect to the blood flow velocity as shown in FIGS. 2 and 3, the average velocity and the variance when the above calculation is executed are as shown in FIGS. 4 and 5, respectively. That is, in the pulsatile flow, the average velocity is distributed in a wide range. on the other hand,
In steady flow, the average velocity is concentrated at a nearly constant value.
【0060】以上の実施例によれば、図4,図5のよう
な特性に応じて着色された表示が行われるので、拍動流
であるか定常流であるかを瞬時に識別することができ
る。尚、図1の構成では、平均速度データと分散データ
との双方でΔT毎の遅延加算を行うようにしたが、いず
れか一方のみで遅延加算を行うような構成も可能であ
る。この場合にも、遅延加算を行ったいずれか一方の流
速データによりCFM表示を行うことで定常流か拍動流
かを識別することが可能であり、回路構成が簡略化でき
る利点も有している。According to the above embodiment, since the display is colored according to the characteristics shown in FIGS. 4 and 5, it is possible to instantly discriminate between the pulsatile flow and the steady flow. it can. In addition, in the configuration of FIG. 1, the delay addition for each ΔT is performed for both the average speed data and the distributed data, but a configuration for performing the delay addition for only one of them is also possible. Also in this case, it is possible to distinguish between the steady flow and the pulsating flow by displaying the CFM based on one of the flow velocity data subjected to the delay addition, and there is also an advantage that the circuit configuration can be simplified. There is.
【0061】以上詳細に説明したように、超音波イメー
ジング表示方法において、Δt間隔で得られた運動物体
反射信号から流速データを求める際に、Δtより大きい
間隔であって異なる時相のΔT毎の運動物体反射信号が
遅延加算されたものにより流速データを求め、このよう
にして求められた流速データによりCFMの超音波イメ
ージング表示を行うようにすることで、拍動流か定常流
かを識別可能な状態で表示することが可能になる。As described in detail above, in the ultrasonic imaging display method, when the flow velocity data is obtained from the moving object reflection signal obtained at the Δt interval, the ΔT is different for each ΔT of the different time phase with the interval larger than Δt. Velocity data can be obtained by delaying and adding moving object reflection signals, and pulsatile flow or steady flow can be identified by performing CFM ultrasonic imaging display based on the thus obtained flow velocity data. It is possible to display in a simple state.
【0062】また、超音波イメージング装置において、
Δt間隔で得られた運動物体反射信号から流速データを
求める際に、Δtより大きい間隔であって異なる時相の
ΔT毎の運動物体反射信号が遅延加算されたものにより
流速データを求め、このようにして求められた流速デー
タによりCFMの超音波イメージング表示を行うように
することで、拍動流か定常流かを識別可能な状態で表示
することが可能になる。In the ultrasonic imaging apparatus,
When the flow velocity data is obtained from the moving object reflection signal obtained at the Δt interval, the flow velocity data is obtained by the delay addition of the moving object reflection signal for each ΔT of an interval larger than Δt and different time phases. By performing the ultrasonic imaging display of the CFM based on the flow velocity data obtained as described above, it becomes possible to display the pulsatile flow and the steady flow in a distinguishable state.
【0063】そして、超音波イメージング装置におい
て、Δt間隔で得られた運動物体反射信号から速度デー
タ及び分散データを求める際に、Δtより大きい間隔で
あって異なる時相のΔT毎の運動物体反射信号が遅延加
算されたものにより速度データ及び分散データを求め、
このようにして求められた速度データ及び分散データに
よりCFMの超音波イメージング表示をすることで、拍
動流か定常流かを識別可能な状態で表示することが可能
になる。Then, in the ultrasonic imaging apparatus, when the velocity data and the dispersion data are obtained from the moving object reflection signal obtained at the Δt intervals, the moving object reflection signal for each ΔT having a different time phase with an interval larger than Δt. The speed data and the distributed data are obtained by the delayed addition of
By performing ultrasonic imaging display of the CFM based on the velocity data and the dispersion data obtained in this way, it becomes possible to display the pulsatile flow and the steady flow in a distinguishable state.
【0064】[0064]
【発明の効果】以上詳細に説明したように、超音波イメ
ージング表示方法において、Δt間隔で得られた運動物
体反射信号から流速データを求める際に、Δtより大き
い間隔であって異なる時相のΔT毎の運動物体反射信号
が遅延加算されたものにより流速データを求め、このよ
うにして求められた流速データによりCFMの超音波イ
メージング表示を行うようにすることで、拍動流か定常
流かを識別可能な状態で表示することが可能になる。As described in detail above, in the ultrasonic imaging display method, when the flow velocity data is obtained from the moving object reflection signal obtained at the Δt interval, the ΔT of the different time phase is obtained at the interval larger than Δt. Velocity data is obtained by delay-adding the moving object reflection signals for each of the signals, and CFM ultrasonic imaging display is performed based on the thus obtained flow velocity data. It becomes possible to display in an identifiable state.
【0065】また、超音波イメージング装置において、
Δt間隔で得られた運動物体反射信号から流速データを
求める際に、Δtより大きい間隔であって異なる時相の
ΔT毎の運動物体反射信号が遅延加算されたものにより
流速データを求め、このようにして求められた流速デー
タによりCFMの超音波イメージング表示を行うように
することで、拍動流か定常流かを識別可能な状態で表示
することが可能になる。In the ultrasonic imaging apparatus,
When the flow velocity data is obtained from the moving object reflection signal obtained at the Δt interval, the flow velocity data is obtained by the delay addition of the moving object reflection signal for each ΔT of an interval larger than Δt and different time phases. By performing the ultrasonic imaging display of the CFM based on the flow velocity data obtained as described above, it becomes possible to display the pulsatile flow and the steady flow in a distinguishable state.
【0066】そして、超音波イメージング装置におい
て、Δt間隔で得られた運動物体反射信号から速度デー
タ及び分散データを求める際に、Δtより大きい間隔で
あって異なる時相のΔT毎の運動物体反射信号が遅延加
算されたものにより速度データ及び分散データを求め、
このようにして求められた速度データ及び分散データに
よりCFMの超音波イメージング表示をすることで、拍
動流か定常流かを識別可能な状態で表示することが可能
になる。Then, in the ultrasonic imaging apparatus, when the velocity data and the dispersion data are obtained from the moving object reflection signals obtained at the Δt intervals, the moving object reflection signals for each ΔT of different time phases with an interval larger than Δt. The speed data and the distributed data are obtained by the delayed addition of
By performing ultrasonic imaging display of the CFM based on the velocity data and the dispersion data obtained in this way, it becomes possible to display the pulsatile flow and the steady flow in a distinguishable state.
【図1】本発明の一実施例の超音波イメージング装置の
ブロック図である。FIG. 1 is a block diagram of an ultrasonic imaging apparatus according to an embodiment of the present invention.
【図2】拍動流の速度スペクトラムを示す特性図であ
る。FIG. 2 is a characteristic diagram showing a velocity spectrum of pulsatile flow.
【図3】定常流の速度スペクトラムを示す特性図であ
る。FIG. 3 is a characteristic diagram showing a velocity spectrum of a steady flow.
【図4】拍動流の平均速度と分散を示す特性図である。FIG. 4 is a characteristic diagram showing average velocity and dispersion of pulsatile flow.
【図5】定常流の平均速度と分散を示す特性図である。FIG. 5 is a characteristic diagram showing the average velocity and dispersion of a steady flow.
【図6】超音波イメージング装置の全体構成を示す構成
図である。FIG. 6 is a configuration diagram showing an overall configuration of an ultrasonic imaging apparatus.
【図7】従来のCFM演算部の回路構成を示す構成図で
ある。FIG. 7 is a configuration diagram showing a circuit configuration of a conventional CFM calculation unit.
【図8】拍動流の速度スペクトラムと平均速度,分散の
様子を示す特性図である。FIG. 8 is a characteristic diagram showing a velocity spectrum of a pulsatile flow, an average velocity, and a state of dispersion.
400,401 MTIフィルタ 402〜404 遅延部 405〜408 乗算器 409 加算器 410 減算器 411〜413 遅延加算器 414 平均速度算出部 415 分散算出部 416〜418 遅延加算部 400, 401 MTI filter 402 to 404 delay unit 405 to 408 multiplier 409 adder 410 subtractor 411 to 413 delay adder 414 average speed calculation unit 415 dispersion calculation unit 416 to 418 delay addition unit
Claims (3)
中の同一音線データのサンプリング間隔Δtより大きい
間隔ΔT毎に所定個数の各運動物体反射信号の遅延加算
を行い、 遅延加算された運動物体反射信号から流速データを求め
てカラーフローマッピングによる超音波イメージング表
示を行うことを特徴とする超音波イメージング表示方
法。1. A delay addition of a predetermined number of moving object reflection signals is performed at every interval ΔT that is larger than the sampling interval Δt of the same sound ray data in the moving object reflection signal obtained from the received signal, and the delay addition is performed. An ultrasonic imaging display method, characterized in that flow velocity data is obtained from a moving object reflection signal and ultrasonic imaging display is performed by color flow mapping.
るBモード処理部と、 受波信号を分析して運動物体反射信号を求め、同一音線
データのサンプリング間隔Δtより大きい間隔ΔT毎に
所定個数の運動物体反射信号の遅延加算を行って流速デ
ータを求めるCFM演算部と、 前記Bモード処理部で得られたBモードデータと前記C
FM演算部で遅延加算されて得られた流速データとを用
いてカラーフローマッピングによる超音波イメージング
表示用のイメージデータを生成する画像処理部とを備え
たことを特徴とする超音波イメージング装置。2. A B-mode processing unit for analyzing received signals to obtain B-mode data, and a received object signal for analyzing moving object reflection signals for each interval ΔT larger than a sampling interval Δt of the same sound ray data. To the C-mode calculation unit for delay-adding a predetermined number of moving object reflection signals to obtain flow velocity data, the B-mode data obtained by the B-mode processing unit, and the C-mode.
An ultrasonic imaging apparatus comprising: an image processing unit that generates image data for ultrasonic imaging display by color flow mapping using flow velocity data obtained by delay addition in the FM calculation unit.
るBモード処理部と、 受波信号を分析して運動物体反射信号を求め、同一音線
データのサンプリング間隔Δtより大きい間隔ΔT毎に
所定個数の運動物体反射信号の遅延加算を行って速度デ
ータ及び分散データを求めるCFM演算部と、 前記Bモード処理部で得られたBモードデータ,前記C
FM演算部で遅延加算されて得られた速度データ及び分
散データを用いてカラーフローマッピングによる超音波
イメージング表示用のイメージデータを生成する画像処
理部とを備えたことを特徴とする超音波イメージング装
置。3. A B-mode processing unit for analyzing a received signal to obtain B-mode data, and a received object signal for analyzing a moving object reflection signal for each interval ΔT larger than a sampling interval Δt of the same sound ray data. A delay time addition of a predetermined number of moving object reflection signals to obtain velocity data and dispersion data; B mode data obtained by the B mode processing unit;
An ultrasonic imaging apparatus, comprising: an image processing unit for generating image data for ultrasonic imaging display by color flow mapping using velocity data and dispersion data obtained by delay addition in the FM calculation unit. .
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP29976694A JP3408647B2 (en) | 1994-12-02 | 1994-12-02 | Ultrasound imaging device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP29976694A JP3408647B2 (en) | 1994-12-02 | 1994-12-02 | Ultrasound imaging device |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH08154939A true JPH08154939A (en) | 1996-06-18 |
JP3408647B2 JP3408647B2 (en) | 2003-05-19 |
Family
ID=17876715
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP29976694A Expired - Fee Related JP3408647B2 (en) | 1994-12-02 | 1994-12-02 | Ultrasound imaging device |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3408647B2 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2020130842A (en) * | 2019-02-25 | 2020-08-31 | キヤノンメディカルシステムズ株式会社 | Ultrasonic image diagnostic apparatus, medical image processing apparatus, and medical image processing program |
-
1994
- 1994-12-02 JP JP29976694A patent/JP3408647B2/en not_active Expired - Fee Related
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2020130842A (en) * | 2019-02-25 | 2020-08-31 | キヤノンメディカルシステムズ株式会社 | Ultrasonic image diagnostic apparatus, medical image processing apparatus, and medical image processing program |
Also Published As
Publication number | Publication date |
---|---|
JP3408647B2 (en) | 2003-05-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5443071A (en) | Quantitative color flow | |
US5445156A (en) | Method for adaptively filtering doppler signals using a complex time domain filter | |
JP5256210B2 (en) | Ultrasonic image processing method and ultrasonic image processing apparatus | |
JP2001178720A (en) | Method and device for visualizing movement in ultrasonic flow imaging using sequential data acquisition | |
US20050054931A1 (en) | Tracking clutter filter for spectral & audio doppler | |
JPH10165400A (en) | Ultrasonic diagnostic device | |
JP3408647B2 (en) | Ultrasound imaging device | |
JP3281435B2 (en) | Ultrasound Doppler diagnostic equipment | |
JP3388043B2 (en) | Ultrasonic imaging display method and ultrasonic imaging apparatus | |
JP2807131B2 (en) | Ultrasound diagnostic equipment | |
JP3847907B2 (en) | Ultrasonic Doppler diagnostic device | |
JP3329589B2 (en) | Ultrasonic imaging display method and ultrasonic diagnostic apparatus | |
JPH02215449A (en) | Ultrasonic diagnostic device | |
JP2703943B2 (en) | Ultrasound diagnostic equipment | |
JP3465197B2 (en) | Ultrasound diagnostic equipment | |
JP3391578B2 (en) | Correlation device and flow information display device | |
JPH0239254B2 (en) | ||
JPH04250148A (en) | Ultrasonic diagnostic device | |
JPH01244738A (en) | Ultrasonic diagnostic device | |
JPH069565B2 (en) | Color flow mapping type ultrasonic diagnostic equipment | |
JP3187169B2 (en) | Ultrasound diagnostic equipment | |
JPH08206113A (en) | Method for processing ultrasonic imaging, and ultrasonic imaging apparatus | |
JPS6096233A (en) | Ultrasonic blood flow measuring apparatus | |
JPS61257631A (en) | Ultrasonic diagnostic apparatus | |
JPS6382634A (en) | Ultrasonic blood flow imaging apparatus |
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 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080314 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090314 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100314 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100314 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100314 Year of fee payment: 7 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100314 Year of fee payment: 7 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100314 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110314 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110314 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120314 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120314 Year of fee payment: 9 |
|
LAPS | Cancellation because of no payment of annual fees |