Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP7613989B2 - Ultrasound diagnostic equipment - Google Patents
[go: Go Back, main page]

JP7613989B2 - Ultrasound diagnostic equipment - Google Patents

Ultrasound diagnostic equipment Download PDF

Info

Publication number
JP7613989B2
JP7613989B2 JP2021063689A JP2021063689A JP7613989B2 JP 7613989 B2 JP7613989 B2 JP 7613989B2 JP 2021063689 A JP2021063689 A JP 2021063689A JP 2021063689 A JP2021063689 A JP 2021063689A JP 7613989 B2 JP7613989 B2 JP 7613989B2
Authority
JP
Japan
Prior art keywords
signal
blood flow
doppler information
information
filter
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.)
Active
Application number
JP2021063689A
Other languages
Japanese (ja)
Other versions
JP2022158648A (en
Inventor
千尋 柴田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Canon Medical Systems Corp
Original Assignee
Canon Medical Systems Corp
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 Canon Medical Systems Corp filed Critical Canon Medical Systems Corp
Priority to JP2021063689A priority Critical patent/JP7613989B2/en
Publication of JP2022158648A publication Critical patent/JP2022158648A/en
Application granted granted Critical
Publication of JP7613989B2 publication Critical patent/JP7613989B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Description

本明細書及び図面に開示の実施形態は、超音波診断装置に関する。 The embodiments disclosed in this specification and the drawings relate to an ultrasound diagnostic device.

従来、超音波診断装置は、ドプラ(Doppler)効果を利用して、超音波の反射波から血流情報の生成及び表示を行う機能を備える。反射波から得られるドプラ情報には、血流に起因する血流信号成分の他に、被検体の臓器等の固定反射や臓器の呼吸性移動、拍動性移動等に起因するクラッタ信号成分が含まれる。このため、従来の超音波診断装置では、被検体の正確な血流情報を生成するために、ドプラ情報からクラッタ信号成分を除去する技術が利用されている。クラッタ信号成分を除去する技術としては、例えば、ウォールフィルタ(Wall Filter)を用いる手法や、ブランク処理(「ブランキング処理」とも呼ばれる)を用いる手法が知られている。 Conventionally, ultrasound diagnostic devices have the function of generating and displaying blood flow information from reflected ultrasound waves using the Doppler effect. The Doppler information obtained from the reflected waves contains not only blood flow signal components caused by blood flow, but also clutter signal components caused by fixed reflections from the subject's organs and respiratory and pulsatile movements of the organs. For this reason, conventional ultrasound diagnostic devices use a technique for removing clutter signal components from the Doppler information in order to generate accurate blood flow information for the subject. Known techniques for removing clutter signal components include, for example, a technique using a wall filter and a technique using blanking processing (also called "blanking processing").

国際公開第2013/153921号International Publication No. 2013/153921

本明細書及び図面に開示の実施形態が解決しようとする課題の一つは、クラッタ信号成分を適切に除去することである。ただし、本明細書及び図面に開示の実施形態により解決しようとする課題は上記課題に限られない。後述する実施形態に示す各構成による各効果に対応する課題を他の課題として位置づけることもできる。 One of the problems that the embodiments disclosed in this specification and the drawings attempt to solve is to appropriately remove clutter signal components. However, the problems that the embodiments disclosed in this specification and the drawings attempt to solve are not limited to the above problem. Problems corresponding to the effects of each configuration shown in the embodiments described below can also be positioned as other problems.

実施形態に係る超音波診断装置は、信号処理部を備える。信号処理部は、互いに異なるカットオフ値を有する複数のフィルタ処理を超音波の受信信号に対してそれぞれ並列に適用することによって抽出された複数のドプラ情報に基づいて、サンプル点ごとにクラッタ信号成分が除去された血流情報を生成する。 The ultrasound diagnostic device according to the embodiment includes a signal processor. The signal processor generates blood flow information from which clutter signal components have been removed for each sample point based on multiple pieces of Doppler information extracted by applying multiple filter processes having different cutoff values in parallel to the received ultrasound signal.

図1は、第1の実施形態に係る超音波診断装置の構成例を示すブロック図である。FIG. 1 is a block diagram showing an example of the configuration of an ultrasonic diagnostic apparatus according to the first embodiment. 図2は、第1の実施形態に係る信号処理回路の構成例を示すブロック図である。FIG. 2 is a block diagram showing an example of the configuration of a signal processing circuit according to the first embodiment. 図3は、第1の実施形態に係る超音波診断装置の処理手順の一例を示すフローチャートである。FIG. 3 is a flowchart showing an example of a processing procedure of the ultrasonic diagnostic apparatus according to the first embodiment. 図4Aは、第1の実施形態に係るブランク処理機能の処理について説明するための図である。FIG. 4A is a diagram for explaining processing of the blanking function according to the first embodiment. 図4Bは、第1の実施形態に係るブランク処理機能の処理について説明するための図である。FIG. 4B is a diagram for explaining the process of the blanking function according to the first embodiment. 図5Aは、第1の実施形態に係るブランク処理機能の処理について説明するための図である。FIG. 5A is a diagram for explaining processing of the blanking function according to the first embodiment. 図5Bは、第1の実施形態に係るブランク処理機能の処理について説明するための図である。FIG. 5B is a diagram for explaining the process of the blanking function according to the first embodiment. 図6Aは、第1の実施形態に係る超音波診断装置による効果を説明するための図である。FIG. 6A is a diagram for explaining the effects of the ultrasound diagnostic apparatus according to the first embodiment. 図6Bは、第1の実施形態に係る超音波診断装置による効果を説明するための図である。FIG. 6B is a diagram for explaining the effect of the ultrasound diagnostic apparatus according to the first embodiment. 図7Aは、比較例に係る第2の手法における課題を説明するための図である。FIG. 7A is a diagram for explaining a problem with the second method according to the comparative example. 図7Bは、比較例に係る第2の手法における課題を説明するための図である。FIG. 7B is a diagram for explaining a problem with the second method according to the comparative example. 図8は、比較例に係る第2の手法における課題を説明するための図である。FIG. 8 is a diagram for explaining a problem with the second method according to the comparative example. 図9は、第1の実施形態の変形例に係る超音波診断装置の処理手順の一例を示すフローチャートである。FIG. 9 is a flowchart showing an example of a processing procedure of an ultrasound diagnostic apparatus according to a modified example of the first embodiment. 図10は、第2の実施形態に係る信号処理回路の処理について説明するための図である。FIG. 10 is a diagram for explaining the processing of the signal processing circuit according to the second embodiment. 図11は、第2の実施形態に係る信号処理回路の構成例を示すブロック図である。FIG. 11 is a block diagram showing an example of the configuration of a signal processing circuit according to the second embodiment. 図12は、第2の実施形態に係る超音波診断装置の処理手順の一例を示すフローチャートである。FIG. 12 is a flowchart showing an example of a processing procedure of the ultrasound diagnostic apparatus according to the second embodiment.

以下、図面を参照して、実施形態に係る超音波診断装置を説明する。なお、実施形態は、以下の実施形態に限られるものではない。また、一つの実施形態に記載した内容は、原則として他の実施形態にも同様に適用可能である。 The ultrasound diagnostic device according to the embodiment will be described below with reference to the drawings. Note that the embodiment is not limited to the following embodiment. Furthermore, the contents described in one embodiment can, in principle, be similarly applied to other embodiments.

(第1の実施形態)
図1は、第1の実施形態に係る超音波診断装置1の構成例を示すブロック図である。図1に示すように、第1の実施形態に係る超音波診断装置1は、装置本体100と、超音波プローブ101と、入力装置102と、ディスプレイ103とを有する。超音波プローブ101、入力装置102、及びディスプレイ103は、装置本体100と通信可能に接続される。なお、被検体Pは、超音波診断装置1の構成に含まれない。
(First embodiment)
Fig. 1 is a block diagram showing an example of the configuration of an ultrasound diagnostic apparatus 1 according to the first embodiment. As shown in Fig. 1, the ultrasound diagnostic apparatus 1 according to the first embodiment includes an apparatus main body 100, an ultrasound probe 101, an input device 102, and a display 103. The ultrasound probe 101, the input device 102, and the display 103 are communicably connected to the apparatus main body 100. Note that a subject P is not included in the configuration of the ultrasound diagnostic apparatus 1.

超音波プローブ101は、超音波の送受信を行う。例えば、超音波プローブ101は、複数の圧電振動子を有する。これら複数の圧電振動子は、後述する装置本体100が有する送受信回路110から供給される駆動信号に基づき超音波を発生する。また、超音波プローブ101が有する複数の圧電振動子は、被検体Pからの反射波を受信して電気信号に変換する。また、超音波プローブ101は、圧電振動子に設けられる整合層と、圧電振動子から後方への超音波の伝播を防止するバッキング材等を有する。なお、超音波プローブ101は、装置本体100と着脱自在に接続される。 The ultrasonic probe 101 transmits and receives ultrasonic waves. For example, the ultrasonic probe 101 has a plurality of piezoelectric transducers. These piezoelectric transducers generate ultrasonic waves based on a drive signal supplied from a transmission/reception circuit 110 of the device main body 100 described later. The plurality of piezoelectric transducers of the ultrasonic probe 101 receive reflected waves from the subject P and convert them into electrical signals. The ultrasonic probe 101 also has a matching layer provided on the piezoelectric transducers, and a backing material that prevents ultrasonic waves from propagating backward from the piezoelectric transducers. The ultrasonic probe 101 is detachably connected to the device main body 100.

超音波プローブ101から被検体Pに超音波が送信されると、送信された超音波は、被検体Pの体内組織における音響インピーダンスの不連続面で次々と反射され、反射波信号として超音波プローブ101が有する複数の圧電振動子にて受信される。受信される反射波信号の振幅は、超音波が反射される不連続面における音響インピーダンスの差に依存する。なお、送信された超音波パルスが、移動している血流や心臓壁等の表面で反射された場合の反射波信号は、ドプラ効果により、移動体の超音波送信方向に対する速度成分に依存して、周波数偏移を受ける。 When ultrasound waves are transmitted from the ultrasound probe 101 to the subject P, the transmitted ultrasound waves are reflected successively by discontinuous surfaces of acoustic impedance in the tissues of the subject P, and are received as reflected wave signals by the multiple piezoelectric transducers of the ultrasound probe 101. The amplitude of the received reflected wave signal depends on the difference in acoustic impedance at the discontinuous surfaces where the ultrasound waves are reflected. When the transmitted ultrasound pulse is reflected by the surface of a moving blood flow, heart wall, etc., the reflected wave signal undergoes a frequency shift due to the Doppler effect depending on the velocity component of the moving body in the direction of ultrasound transmission.

なお、第1の実施形態に係る超音波プローブ101は、被検体Pを2次元で走査する1Dアレイプローブであっても、被検体Pを3次元で走査するメカニカル4Dプローブや2Dアレイプローブであっても適用可能である。 The ultrasound probe 101 according to the first embodiment can be a 1D array probe that scans the subject P in two dimensions, or a mechanical 4D probe or 2D array probe that scans the subject P in three dimensions.

入力装置102は、マウス、キーボード、ボタン、パネルスイッチ、タッチコマンドスクリーン、フットスイッチ、トラックボール、ジョイスティック等の装置に対応する。入力装置102は、超音波診断装置1の操作者からの各種設定要求を受け付け、装置本体100に対して受け付けた各種設定要求を転送する。 The input device 102 corresponds to devices such as a mouse, a keyboard, a button, a panel switch, a touch command screen, a foot switch, a trackball, a joystick, etc. The input device 102 accepts various setting requests from the operator of the ultrasound diagnostic device 1 and transfers the accepted various setting requests to the device main body 100.

ディスプレイ103は、超音波診断装置1の操作者が入力装置102を用いて各種設定要求を入力するためのGUI(Graphical User Interface)を表示したり、装置本体100において生成された超音波画像データ等を表示したりする。 The display 103 displays a GUI (Graphical User Interface) that enables the operator of the ultrasound diagnostic device 1 to input various setting requests using the input device 102, and displays ultrasound image data generated in the device body 100, etc.

装置本体100は、超音波プローブ101が受信した反射波信号に基づいて超音波画像データを生成する装置である。図1に示す装置本体100により生成される超音波画像データは、2次元の反射波信号に基づいて生成される2次元の超音波画像データであっても、3次元の反射波信号に基づいて生成される3次元の超音波画像データであってもよい。 The device body 100 is a device that generates ultrasound image data based on a reflected wave signal received by the ultrasound probe 101. The ultrasound image data generated by the device body 100 shown in FIG. 1 may be two-dimensional ultrasound image data generated based on a two-dimensional reflected wave signal, or three-dimensional ultrasound image data generated based on a three-dimensional reflected wave signal.

装置本体100は、図1に示すように、送受信回路110と、信号処理回路120と、画像生成回路130と、画像メモリ140と、記憶回路150と、処理回路160とを備える。送受信回路110、信号処理回路120、画像生成回路130、画像メモリ140、記憶回路150、及び処理回路160は、互いに通信可能に接続される。 As shown in FIG. 1, the device main body 100 includes a transmission/reception circuit 110, a signal processing circuit 120, an image generation circuit 130, an image memory 140, a storage circuit 150, and a processing circuit 160. The transmission/reception circuit 110, the signal processing circuit 120, the image generation circuit 130, the image memory 140, the storage circuit 150, and the processing circuit 160 are connected to each other so as to be able to communicate with each other.

送受信回路110は、後述する処理回路160の指示に基づいて、超音波プローブ101が行う超音波送受信(超音波スキャン)を制御する。送受信回路110は、パルス発生器、送信遅延回路、パルサ等を有し、超音波プローブ101に駆動信号を供給する。パルス発生器は、所定の繰り返し周波数で送信超音波を形成するためのレートパルスを繰り返し発生する。また、送信遅延回路は、超音波プローブ101から発生される超音波をビーム状に集束し、かつ送信指向性を決定するために必要な圧電振動子ごとの遅延時間を、パルス発生器が発生する各レートパルスに対し与える。また、パルサは、レートパルスに基づくタイミングで、超音波プローブ101に駆動信号(駆動パルス)を印加する。すなわち、送信遅延回路は、各レートパルスに対し与える遅延時間を変化させることで、圧電振動子面から送信される超音波の送信方向を任意に調整する。 The transmission/reception circuit 110 controls the ultrasonic transmission/reception (ultrasonic scan) performed by the ultrasonic probe 101 based on instructions from the processing circuit 160 described later. The transmission/reception circuit 110 has a pulse generator, a transmission delay circuit, a pulser, etc., and supplies a drive signal to the ultrasonic probe 101. The pulse generator repeatedly generates a rate pulse for forming a transmission ultrasonic wave at a predetermined repetition frequency. The transmission delay circuit also focuses the ultrasonic waves generated from the ultrasonic probe 101 into a beam shape and provides each rate pulse generated by the pulse generator with a delay time for each piezoelectric transducer required to determine the transmission directivity. The pulser also applies a drive signal (drive pulse) to the ultrasonic probe 101 at a timing based on the rate pulse. In other words, the transmission delay circuit arbitrarily adjusts the transmission direction of the ultrasonic waves transmitted from the piezoelectric transducer surface by changing the delay time provided to each rate pulse.

なお、送受信回路110は、後述する処理回路160の指示に基づいて、所定のスキャンシーケンスを実行するために、送信周波数、送信駆動電圧等を瞬時に変更可能な機能を有している。特に、送信駆動電圧の変更は、瞬間にその値を切り替え可能なリニアアンプ型の発信回路、又は、複数の電源ユニットを電気的に切り替える機構によって実現される。 The transmission/reception circuit 110 has the ability to instantly change the transmission frequency, transmission drive voltage, etc., in order to execute a specified scan sequence based on instructions from the processing circuit 160, which will be described later. In particular, the change in transmission drive voltage is achieved by a linear amplifier type transmission circuit that can instantly switch its value, or a mechanism that electrically switches between multiple power supply units.

また、送受信回路110は、アンプ回路、A/D(Analog/Digital)変換器、受信遅延回路、加算器、直交検波回路等を有し、超音波プローブ101が受信した反射波信号に対して各種処理を行って反射波データを生成する。アンプ回路は、反射波信号をチャンネル毎に増幅してゲイン補正処理を行う。A/D変換器は、ゲイン補正された反射波信号をA/D変換する。受信遅延回路は、デジタルデータに受信指向性を決定するのに必要な受信遅延時間を与える。加算器は、受信遅延回路により受信遅延時間が与えられた反射波信号の加算処理を行う。加算器の加算処理により、反射波信号の受信指向性に応じた方向からの反射成分が強調される。 The transmission/reception circuit 110 also has an amplifier circuit, an A/D (Analog/Digital) converter, a reception delay circuit, an adder, a quadrature detection circuit, etc., and performs various processes on the reflected wave signal received by the ultrasound probe 101 to generate reflected wave data. The amplifier circuit amplifies the reflected wave signal for each channel and performs gain correction processing. The A/D converter A/D converts the gain-corrected reflected wave signal. The reception delay circuit gives the digital data a reception delay time required to determine the reception directivity. The adder performs addition processing of the reflected wave signal to which the reception delay time has been given by the reception delay circuit. The addition processing by the adder emphasizes the reflected component from a direction corresponding to the reception directivity of the reflected wave signal.

そして、直交検波回路は、加算器の出力信号をベースバンド帯域の同相信号(I信号、I:In-phase)と直交信号(Q信号、Q:Quadrature-phase)とに変換する。そして、直交検波回路は、I信号及びQ信号(以下、IQ信号と記載する)を反射波データとして、信号処理回路120へ送信する。なお、直交検波回路は、加算器の出力信号を、RF(Radio Frequency)信号に変換した上で、信号処理回路120へ送信してもよい。IQ信号や、RF信号は、位相情報が含まれる信号(受信信号)となる。以下では、送受信回路110が出力する反射波データを、受信信号と記載する場合がある。 The quadrature detection circuit then converts the output signal of the adder into an in-phase signal (I signal, I: In-phase) and a quadrature signal (Q signal, Q: Quadrature-phase) in the baseband. The quadrature detection circuit then transmits the I signal and Q signal (hereinafter referred to as IQ signal) as reflected wave data to the signal processing circuit 120. Note that the quadrature detection circuit may convert the output signal of the adder into an RF (Radio Frequency) signal before transmitting it to the signal processing circuit 120. The IQ signal and RF signal are signals (received signals) that contain phase information. Hereinafter, the reflected wave data output by the transmission/reception circuit 110 may be referred to as a received signal.

なお、送受信回路110は、1回の超音波ビームの送信により得られる各圧電振動子の反射波信号から複数の受信フォーカスの反射波データを生成することができる。すなわち、送受信回路110は、並列同時受信処理を行うことが可能な回路である。なお、第1の実施形態は、送受信回路110が並列同時受信処理を実行できない場合であっても適用可能である。 The transmission/reception circuit 110 can generate reflected wave data for multiple reception foci from the reflected wave signals of each piezoelectric transducer obtained by transmitting an ultrasonic beam once. In other words, the transmission/reception circuit 110 is a circuit capable of performing parallel simultaneous reception processing. The first embodiment can be applied even when the transmission/reception circuit 110 cannot perform parallel simultaneous reception processing.

信号処理回路120は、送受信回路110が反射波信号から生成した反射波データに対して、各種の信号処理を行う。例えば、信号処理回路120は、送受信回路110が生成した反射波データ(IQ信号)を一時的に記憶するバッファを備える。バッファは、数フレーム分のIQ信号、又は、数ボリューム分のIQ信号を記憶する。例えば、バッファは、FIFO(First-In/First-Out)メモリであり、所定フレーム分のIQ信号を記憶する。そして、例えば、バッファは、新たに1フレーム分のIQ信号が送受信回路110にて生成された場合、生成時間が最も古い1フレーム分のIQ信号を破棄して、新たに生成された1フレーム分のI/Q信号を記憶する。そして、信号処理回路120は、バッファから読み出した反射波データ(IQ信号)に対して、対数増幅、包絡線検波処理、対数圧縮などを行って、多点の信号強度が輝度の明るさで表現されるデータ(Bモード情報)を生成する。 The signal processing circuit 120 performs various signal processing on the reflected wave data generated by the transmission/reception circuit 110 from the reflected wave signal. For example, the signal processing circuit 120 includes a buffer that temporarily stores the reflected wave data (IQ signal) generated by the transmission/reception circuit 110. The buffer stores several frames of IQ signals or several volumes of IQ signals. For example, the buffer is a FIFO (First-In/First-Out) memory that stores a predetermined number of frames of IQ signals. Then, for example, when a new frame of IQ signals is generated by the transmission/reception circuit 110, the buffer discards the oldest frame of IQ signals and stores the newly generated frame of I/Q signals. Then, the signal processing circuit 120 performs logarithmic amplification, envelope detection processing, logarithmic compression, etc. on the reflected wave data (IQ signal) read from the buffer to generate data (B-mode information) in which the signal intensity at multiple points is expressed as luminance brightness.

なお、信号処理回路120は、フィルタ処理により、検波周波数を変化させることで、映像化する周波数帯域を変えることができる。この信号処理回路120のフィルタ処理機能を用いることにより、コントラストハーモニックイメージング(CHI:Contrast Harmonic Imaging)や、ティッシュハーモニックイメージング(THI:Tissue Harmonic Imaging)等のハーモニックイメージングを実行可能である。 The signal processing circuit 120 can change the detection frequency through filter processing, thereby changing the frequency band to be visualized. By using the filter processing function of the signal processing circuit 120, it is possible to perform harmonic imaging such as contrast harmonic imaging (CHI) and tissue harmonic imaging (THI).

また、この信号処理回路120のフィルタ処理機能を用いることにより、第1の実施形態に係る超音波診断装置1は、ティッシュハーモニックイメージング(THI:Tissue Harmonic Imaging)を実行可能である。 In addition, by using the filter processing function of this signal processing circuit 120, the ultrasound diagnostic device 1 according to the first embodiment can perform tissue harmonic imaging (THI).

また、CHIやTHIのハーモニックイメージングを行う際、信号処理回路120は、上述したフィルタ処理を用いた方法とは異なる方法により、ハーモニック成分を抽出することができる。ハーモニックイメージングでは、振幅変調(AM:Amplitude Modulation)法や位相変調(PM:Phase Modulation)法、AM法及びPM法を組み合わせたAMPM法と呼ばれる映像法が行われる。AM法、PM法及びAMPM法では、同一の走査線に対して振幅や位相が異なる超音波送信を複数回行う。これにより、送受信回路110は、各走査線で複数の反射波データ(受信信号)を生成し出力する。そして、信号処理回路120は、各走査線の複数の反射波データ(受信信号)を、変調法に応じた加減算処理することで、ハーモニック成分を抽出する。そして、信号処理回路120は、ハーモニック成分の反射波データ(受信信号)に対して包絡線検波処理等を行って、Bモード情報を生成する。 When performing harmonic imaging of CHI or THI, the signal processing circuit 120 can extract harmonic components using a method other than the above-mentioned filter processing method. In harmonic imaging, an imaging method called the AMPM method, which is a combination of the amplitude modulation (AM) method, the phase modulation (PM) method, and the AM method and PM method, is performed. In the AM method, the PM method, and the AMPM method, ultrasonic waves with different amplitudes and phases are transmitted multiple times for the same scanning line. As a result, the transmission/reception circuit 110 generates and outputs multiple reflected wave data (received signals) for each scanning line. Then, the signal processing circuit 120 extracts harmonic components by performing addition and subtraction processing on the multiple reflected wave data (received signals) for each scanning line according to the modulation method. Then, the signal processing circuit 120 performs envelope detection processing or the like on the reflected wave data (received signals) of the harmonic components to generate B-mode information.

信号処理回路120は、バッファから読み出した反射波データを周波数解析することで、走査範囲内にある移動体のドプラ効果に基づく運動情報を抽出したデータ(ドプラ情報)を生成する。具体的には、信号処理回路120は、移動体の運動情報として、平均速度、平均分散値、平均パワー値等を、複数のサンプル点それぞれで推定したドプラ情報を生成する。ここで、移動体とは、例えば、血流や、心壁等の組織、造影剤である。血流の運動情報(血流情報)を生成(抽出)する場合、信号処理回路120は、被検体の臓器等の固定反射や臓器の呼吸性移動、拍動性移動等に起因するクラッタ信号成分をドプラ情報から除去する処理を行う。血流情報は、血流の平均速度、血流の平均分散値、血流の平均パワー値等で表される。 The signal processing circuit 120 performs frequency analysis on the reflected wave data read from the buffer to generate data (Doppler information) that extracts motion information based on the Doppler effect of a moving object within the scanning range. Specifically, the signal processing circuit 120 generates Doppler information that estimates the average velocity, average variance, average power value, etc., as the motion information of the moving object at each of multiple sample points. Here, the moving object is, for example, blood flow, tissue such as the heart wall, or contrast agent. When generating (extracting) motion information of blood flow (blood flow information), the signal processing circuit 120 performs processing to remove clutter signal components caused by fixed reflections of the subject's organs, respiratory movement of the organs, pulsatile movement, etc., from the Doppler information. The blood flow information is expressed as the average velocity, average variance, and average power value of the blood flow.

上記の信号処理回路120の機能を用いて、本実施形態に係る超音波診断装置1は、カラーフローマッピング法(CFM:Color Flow Mapping)とも呼ばれるカラードプラ法を実行可能である。CFM法では、超音波の送受信が複数の走査線上で複数回行われる。そして、CFM法では、同一位置のデータ列に対してMTI(Moving Target Indicator)フィルタを掛けることで、静止している組織、或いは、動きの遅い組織に由来する信号(クラッタ信号)を抑制して、血流に由来する信号を抽出する。そして、CFM法では、この血流信号から血流の速度、血流の分散、血流のパワー等の血流情報を推定する。後述する画像生成回路130は、推定結果の分布を、例えば、2次元でカラー表示した超音波画像データ(ドプラ画像データ)を生成する。そして、ディスプレイ103は、ドプラ画像データを表示する。 Using the above-mentioned function of the signal processing circuit 120, the ultrasound diagnostic device 1 according to this embodiment can perform a color Doppler method, also called a color flow mapping method (CFM). In the CFM method, ultrasonic waves are transmitted and received multiple times on multiple scanning lines. In the CFM method, a moving target indicator (MTI) filter is applied to a data sequence at the same position to suppress signals (clutter signals) originating from stationary tissue or slow-moving tissue and extract signals originating from blood flow. In the CFM method, blood flow information such as the blood flow speed, blood flow dispersion, and blood flow power is estimated from this blood flow signal. The image generation circuit 130, which will be described later, generates ultrasound image data (Doppler image data) that displays the distribution of the estimation results in color, for example, in two dimensions. The display 103 displays the Doppler image data.

MTIフィルタとしては、バタワース型のIIR(Infinite Impulse Response)フィルタや、多項式回帰フィルタ(Polynomial Regression Filter)等、係数が固定されたフィルタや、入力信号に応じて係数を変化させる適応型のMTIフィルタが適用可能である。適応型のMTIフィルタとしては、例えば「Eigenvector Regression Filter」と呼ばれているフィルタを用いる。以下、固有ベクトルを用いた適応型MTIフィルタである「Eigenvector Regression Filter」を、「固有ベクトル型MTIフィルタ」と記載する。 As the MTI filter, a Butterworth-type IIR (Infinite Impulse Response) filter, a polynomial regression filter, or other filter with fixed coefficients, or an adaptive MTI filter that changes coefficients according to the input signal can be used. As an adaptive MTI filter, for example, a filter called an "Eigenvector Regression Filter" is used. Hereinafter, the "Eigenvector Regression Filter", which is an adaptive MTI filter that uses eigenvectors, will be referred to as an "Eigenvector-type MTI filter".

固有ベクトル型MTIフィルタは、相関行列から固有ベクトルを計算し、計算した固有ベクトルから、クラッタ成分抑制処理に用いる係数を計算する。この方法は、主成分分析や、カルーネン・レーベル変換(Karhunen-Loeve transform)、固有空間法で使われている手法を応用したものである。 The eigenvector-type MTI filter calculates eigenvectors from the correlation matrix, and then calculates the coefficients used in the clutter component suppression process from the calculated eigenvectors. This method is an application of the techniques used in principal component analysis, the Karhunen-Loeve transform, and the eigenspace method.

固有ベクトル型MTIフィルタを用いる場合、信号処理回路120は、同一位置(同一サンプル点)の連続した反射波データのデータ列から、走査範囲の相関行列を計算する。例えば、信号処理回路120は、相関行列の固有値及び当該固有値に対応する固有ベクトルを計算する。そして、信号処理回路120は、例えば、各固有値の大きさに基づいて各固有ベクトルを並べた行列のランクを低減した行列を、クラッタ成分を抑制するフィルタ行列として計算する。ここで、信号処理回路120は、例えば、予め設定された値、或いは、操作者が指定した値により、低減される主成分の数、すなわち、ランクカット数の値を決定する。しかし、心臓や血管等、拍動により移動速度が時間により変化する組織が走査範囲内に含まれる場合、ランクカット数の値は、固有値の大きさから適応的に決定されることが好適である。すなわち、信号処理回路120は、相関行列の固有値の大きさに応じて、低減する主成分の数を変更する。本実施形態では、信号処理回路120は、固有値の大きさに応じて、低減するランク数を変更する。 When using an eigenvector type MTI filter, the signal processing circuit 120 calculates the correlation matrix of the scanning range from a data sequence of consecutive reflected wave data at the same position (same sample point). For example, the signal processing circuit 120 calculates the eigenvalues of the correlation matrix and the eigenvectors corresponding to the eigenvalues. Then, the signal processing circuit 120 calculates a matrix in which the rank of the matrix in which the eigenvectors are arranged based on the magnitude of each eigenvalue is reduced as a filter matrix that suppresses clutter components. Here, the signal processing circuit 120 determines the number of principal components to be reduced, that is, the value of the rank cut number, for example, by a preset value or a value specified by the operator. However, when tissues whose moving speed changes over time due to pulsation, such as the heart and blood vessels, are included within the scanning range, it is preferable that the value of the rank cut number is adaptively determined from the magnitude of the eigenvalue. That is, the signal processing circuit 120 changes the number of principal components to be reduced depending on the magnitude of the eigenvalue of the correlation matrix. In this embodiment, the signal processing circuit 120 changes the number of ranks to be reduced depending on the magnitude of the eigenvalue.

信号処理回路120は、フィルタ行列を用いて、同一位置(同一サンプル点)の連続した反射波データのデータ列から、クラッタ成分が抑制され、血流に由来する血流信号が抽出されたデータ列を出力する。信号処理回路120は、出力したデータを用いた自己相関演算等の演算を行って、血流情報を推定し、推定した血流情報を出力する。 The signal processing circuit 120 uses a filter matrix to suppress clutter components from a data string of consecutive reflected wave data from the same position (same sample point) and outputs a data string in which blood flow signals derived from blood flow are extracted. The signal processing circuit 120 performs calculations such as autocorrelation calculations using the output data to estimate blood flow information and output the estimated blood flow information.

画像生成回路130は、信号処理回路120が生成したデータから超音波画像データを生成する。例えば、画像生成回路130は、2次元のBモード情報から反射波の強度を輝度で表した2次元Bモード画像データを生成する。また、画像生成回路130は、2次元のドプラ情報から血流情報が映像化された2次元ドプラ画像データを生成する。2次元ドプラ画像データは、速度画像データ、分散画像データ、パワー画像データ、又は、これらを組み合わせた画像データである。画像生成回路130は、ドプラ画像データとして、血流情報がカラーで表示されるカラードプラ画像データを生成したり、1つの血流情報がグレースケールで表示されるドプラ画像データを生成したりする。 The image generation circuit 130 generates ultrasound image data from the data generated by the signal processing circuit 120. For example, the image generation circuit 130 generates two-dimensional B-mode image data in which the intensity of the reflected wave is expressed as brightness from the two-dimensional B-mode information. The image generation circuit 130 also generates two-dimensional Doppler image data in which blood flow information is visualized from the two-dimensional Doppler information. The two-dimensional Doppler image data is velocity image data, variance image data, power image data, or image data that combines these. The image generation circuit 130 generates color Doppler image data in which blood flow information is displayed in color, or Doppler image data in which one piece of blood flow information is displayed in grayscale, as the Doppler image data.

ここで、画像生成回路130は、一般的には、超音波走査の走査線信号列を、テレビ等に代表されるビデオフォーマットの走査線信号列に変換(スキャンコンバート)し、表示用の超音波画像データを生成する。具体的には、画像生成回路130は、超音波プローブ101による超音波の走査形態に応じて座標変換を行うことで、表示用の超音波画像データを生成する。また、画像生成回路130は、スキャンコンバート以外に、種々の画像処理として、例えば、スキャンコンバート後の複数の画像フレームを用いて、輝度の平均値画像を再生成する画像処理(平滑化処理)や、画像内で微分フィルタを用いる画像処理(エッジ強調処理)等を行う。また、画像生成回路130は、超音波画像データに、種々のパラメータの文字情報、目盛り、ボディーマーク等を合成する。 Here, the image generating circuit 130 generally converts (scan converts) the scan line signal sequence of the ultrasound scan into a scan line signal sequence of a video format, such as that of a television, to generate ultrasound image data for display. Specifically, the image generating circuit 130 generates ultrasound image data for display by performing coordinate conversion according to the ultrasound scanning form of the ultrasound probe 101. In addition to scan conversion, the image generating circuit 130 also performs various image processing, such as image processing (smoothing processing) that regenerates an average brightness image using multiple image frames after scan conversion, and image processing (edge enhancement processing) that uses a differential filter within the image. The image generating circuit 130 also combines text information of various parameters, scales, body marks, etc., with the ultrasound image data.

すなわち、Bモード情報及びドプラ情報は、スキャンコンバート処理前の超音波画像データであり、画像生成回路130が生成するデータは、スキャンコンバート処理後の表示用の超音波画像データである。なお、Bモード情報及びドプラ情報は、生データ(Raw Data)とも呼ばれる。画像生成回路130は、スキャンコンバート処理前の2次元超音波画像データから、表示用の2次元超音波画像データを生成する。 That is, the B-mode information and Doppler information are ultrasound image data before the scan conversion process, and the data generated by the image generation circuit 130 is ultrasound image data for display after the scan conversion process. Note that the B-mode information and Doppler information are also called raw data. The image generation circuit 130 generates two-dimensional ultrasound image data for display from the two-dimensional ultrasound image data before the scan conversion process.

更に、画像生成回路130は、信号処理回路120が生成した3次元のBモード情報に対して座標変換を行うことで、3次元Bモード画像データを生成する。また、画像生成回路130は、信号処理回路120が生成した3次元のドプラ情報に対して座標変換を行うことで、3次元ドプラ画像データを生成する。 Furthermore, the image generating circuit 130 generates three-dimensional B-mode image data by performing coordinate transformation on the three-dimensional B-mode information generated by the signal processing circuit 120. Furthermore, the image generating circuit 130 generates three-dimensional Doppler image data by performing coordinate transformation on the three-dimensional Doppler information generated by the signal processing circuit 120.

更に、画像生成回路130は、ボリュームデータをディスプレイ103にて表示するための各種の2次元画像データを生成するために、ボリュームデータに対してレンダリング処理を行う。画像生成回路130が行うレンダリング処理としては、例えば、断面再構成法(MPR:Multi Planer Reconstruction)を行ってボリュームデータからMPR画像データを生成する処理がある。また、画像生成回路130が行うレンダリング処理としては、例えば、3次元の情報を反映した2次元画像データを生成するボリュームレンダリング(VR:Volume Rendering)処理がある。 Furthermore, the image generation circuit 130 performs rendering processing on the volume data to generate various types of two-dimensional image data for displaying the volume data on the display 103. The rendering processing performed by the image generation circuit 130 includes, for example, processing for performing multi-planar reconstruction (MPR) to generate MPR image data from the volume data. The rendering processing performed by the image generation circuit 130 also includes, for example, volume rendering (VR) processing for generating two-dimensional image data reflecting three-dimensional information.

画像メモリ140は、画像生成回路130が生成した表示用の画像データを記憶するメモリである。また、画像メモリ140は、信号処理回路120や信号処理回路120が生成したデータを記憶することも可能である。画像メモリ140が記憶するBモード情報やドプラ情報は、例えば、診断の後に操作者が呼び出すことが可能となっており、画像生成回路130を経由して表示用の超音波画像データとなる。また、画像メモリ140は、送受信回路110が出力した反射波データを記憶することも可能である。 The image memory 140 is a memory that stores image data for display generated by the image generation circuit 130. The image memory 140 can also store the signal processing circuit 120 and data generated by the signal processing circuit 120. The B-mode information and Doppler information stored in the image memory 140 can be called up by the operator after diagnosis, for example, and becomes ultrasound image data for display via the image generation circuit 130. The image memory 140 can also store reflected wave data output by the transmission/reception circuit 110.

記憶回路150は、超音波送受信、画像処理及び表示処理を行うための制御プログラムや、診断情報(例えば、患者ID、医師の所見等)や、診断プロトコルや各種ボディーマーク等の各種データを記憶する。また、記憶回路150は、必要に応じて、画像メモリ140が記憶する画像データの保管等にも使用される。また、記憶回路150が記憶するデータは、図示しないインターフェースを経由して、外部装置へ転送することができる。また、記憶回路150は、外部装置から図示しないインターフェースを経由して転送されたデータを記憶することも可能である。 The memory circuitry 150 stores various data such as control programs for transmitting and receiving ultrasound, image processing, and display processing, diagnostic information (e.g., patient ID, doctor's findings, etc.), diagnostic protocols, and various body marks. The memory circuitry 150 is also used to store image data stored in the image memory 140 as necessary. The data stored in the memory circuitry 150 can be transferred to an external device via an interface not shown. The memory circuitry 150 can also store data transferred from an external device via an interface not shown.

処理回路160は、超音波診断装置1の処理全体を制御する。具体的には、処理回路160は、入力装置102を介して操作者から入力された各種設定要求や、記憶回路150から読込んだ各種制御プログラム及び各種データに基づき、送受信回路110、信号処理回路120、及び画像生成回路130の処理を制御する。また、処理回路160は、画像メモリ140や記憶回路150が記憶する表示用の超音波画像データをディスプレイ103にて表示するように制御する。 The processing circuitry 160 controls the overall processing of the ultrasound diagnostic device 1. Specifically, the processing circuitry 160 controls the processing of the transmission/reception circuitry 110, the signal processing circuitry 120, and the image generation circuitry 130 based on various setting requests input by the operator via the input device 102 and various control programs and various data read from the memory circuitry 150. The processing circuitry 160 also controls the display 103 to display ultrasound image data for display stored in the image memory 140 and the memory circuitry 150.

ここで、処理回路160が実行する各処理機能は、例えば、コンピュータによって実行可能なプログラムの形態で記憶回路150に記録されている。処理回路160は、各プログラムを記憶回路150から読み出し、実行することで各プログラムに対応する機能を実現するプロセッサである。 Here, each processing function executed by the processing circuitry 160 is recorded in the memory circuitry 150, for example, in the form of a program executable by a computer. The processing circuitry 160 is a processor that realizes the function corresponding to each program by reading each program from the memory circuitry 150 and executing it.

本実施形態にて用いる「プロセッサ」という文言は、例えば、CPU(Central Processing Unit)、GPU(Graphics Processing Unit)、或いは、特定用途向け集積回路(Application Specific Integrated Circuit:ASIC)、プログラマブル論理デバイス(例えば、単純プログラマブル論理デバイス(Simple Programmable Logic Device:SPLD)、複合プログラマブル論理デバイス(Complex Programmable Logic Device:CPLD)、及びフィールドプログラマブルゲートアレイ(Field Programmable Gate Array:FPGA))等の回路を意味する。プロセッサは記憶回路に保存されたプログラムを読み出し実行することで機能を実現する。なお、記憶回路150にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むよう構成しても構わない。この場合、プロセッサは回路内に組み込まれたプログラムを読み出し実行することで機能を実現する。なお、本実施形態の各プロセッサは、プロセッサごとに単一の回路として構成される場合に限らず、複数の独立した回路を組み合わせて1つのプロセッサとして構成し、その機能を実現するようにしてもよい。更に、各図における複数の構成要素を1つのプロセッサへ統合してその機能を実現するようにしてもよい。 The term "processor" used in this embodiment means, for example, a circuit such as a CPU (Central Processing Unit), a GPU (Graphics Processing Unit), an Application Specific Integrated Circuit (ASIC), a programmable logic device (for example, a Simple Programmable Logic Device (SPLD), a Complex Programmable Logic Device (CPLD), and a Field Programmable Gate Array (FPGA)). The processor realizes its function by reading and executing a program stored in a memory circuit. Note that instead of storing a program in the memory circuit 150, the program may be directly built into the circuit of the processor. In this case, the processor realizes its function by reading and executing a program built into the circuit. Note that each processor in this embodiment is not limited to being configured as a single circuit for each processor, but may be configured as a single processor by combining multiple independent circuits to realize its function. Furthermore, multiple components in each figure may be integrated into a single processor to realize its function.

ここで、第1の実施形態に係る超音波診断装置1は、クラッタ信号成分を適切に除去するために、以下に説明する処理機能を備える。すなわち、信号処理回路120は、互いに異なるカットオフ値を有する複数のフィルタ処理を超音波の受信信号に対してそれぞれ並列に適用することによって抽出された複数のドプラ情報に基づいて、サンプル点ごとにクラッタ信号成分が除去された血流情報を生成する。なお、信号処理回路120は、信号処理部の一例である。 The ultrasound diagnostic device 1 according to the first embodiment has the processing functions described below in order to appropriately remove clutter signal components. That is, the signal processing circuit 120 generates blood flow information from which clutter signal components have been removed for each sample point based on multiple pieces of Doppler information extracted by applying multiple filter processes having different cutoff values in parallel to the received ultrasound signal. The signal processing circuit 120 is an example of a signal processing unit.

図2は、第1の実施形態に係る信号処理回路120の構成例を示すブロック図である。図2に示すように、信号処理回路120は、ウォールフィルタ処理機能121A,121B、自己相関演算機能122A,122B、ドプラ情報演算機能123A,123B、ブランク処理機能124、及びパーシスタンス機能125を実行する。 Figure 2 is a block diagram showing an example of the configuration of the signal processing circuit 120 according to the first embodiment. As shown in Figure 2, the signal processing circuit 120 executes wall filter processing functions 121A, 121B, autocorrelation calculation functions 122A, 122B, Doppler information calculation functions 123A, 123B, a blank processing function 124, and a persistence function 125.

なお、信号処理回路120の構成要素であるウォールフィルタ処理機能121A,121B、自己相関演算機能122A,122B、ドプラ情報演算機能123A,123B、ブランク処理機能124、及びパーシスタンス機能125が実行する各処理機能は、例えば、コンピュータによって実行可能なプログラムの形態で記憶回路150に記録されている。処理回路160は、各プログラムを記憶回路150から読み出し、実行することで各プログラムに対応する機能を実現するプロセッサである。なお、信号処理回路120が実行する各処理機能については、図3の説明にて後述する。 The processing functions executed by the wall filter processing functions 121A, 121B, autocorrelation calculation functions 122A, 122B, Doppler information calculation functions 123A, 123B, blank processing function 124, and persistence function 125, which are components of the signal processing circuit 120, are recorded in the memory circuit 150, for example, in the form of a program executable by a computer. The processing circuit 160 is a processor that reads out each program from the memory circuit 150 and executes it to realize the function corresponding to each program. The processing functions executed by the signal processing circuit 120 will be described later in the explanation of FIG. 3.

図3は、第1の実施形態に係る超音波診断装置1の処理手順の一例を示すフローチャートである。図3に示すステップS101~ステップS109の各処理は、例えば、ドプラ画像データの表示を開始する旨の開始指示を操作者から受け付けた場合に開始される。 Figure 3 is a flowchart showing an example of a processing procedure of the ultrasound diagnostic device 1 according to the first embodiment. Each process of steps S101 to S109 shown in Figure 3 is started, for example, when a start instruction to start displaying Doppler image data is received from the operator.

図3に示すように、超音波診断装置1は、超音波スキャンを実行する(ステップS101)。例えば、送受信回路110は、処理回路160からの制御に基づいて、操作者によって設定された関心領域に対して超音波スキャンを実行する。そして、送受信回路110は、フレームごとに、関心領域内の各サンプル点からの反射波データを生成する。 As shown in FIG. 3, the ultrasound diagnostic device 1 executes an ultrasound scan (step S101). For example, the transmission/reception circuitry 110 executes an ultrasound scan on a region of interest set by an operator based on control from the processing circuitry 160. Then, the transmission/reception circuitry 110 generates reflected wave data from each sample point in the region of interest for each frame.

続いて、信号処理回路120は、高カットオフ値のフィルタ処理Aと、低カットオフ値のフィルタ処理Bとをそれぞれ適用する(ステップS102)。ステップS102の処理は、ウォールフィルタ処理機能121A,121Bによって実行される。 Then, the signal processing circuit 120 applies a high cutoff value filter process A and a low cutoff value filter process B (step S102). The process of step S102 is executed by the wall filter processing functions 121A and 121B.

ウォールフィルタ処理機能121A,121Bは、サンプル点ごとに時系列的に受信された超音波の受信信号(反射波データのデータ列)に対してフィルタ処理を実行する。ここで、ウォールフィルタ処理機能121Aは、反射波信号のデータ列に対して、カットオフ値Aを有するフィルタ処理Aを実行する。また、ウォールフィルタ処理機能121Bは、反射波信号のデータ列に対して、カットオフ値Aより小さいカットオフ値Bを有するフィルタ処理Bを実行する。なお、ウォールフィルタ処理機能121Aは、第1フィルタ処理部の一例である。また、ウォールフィルタ処理機能121Bは、第2フィルタ処理部の一例である。また、カットオフ値Aは、第1カットオフ値の一例である。また、カットオフ値Bは、第2カットオフ値の一例である。また、フィルタ処理Aは、第1フィルタ処理の一例である。また、フィルタ処理Bは、第2フィルタ処理の一例である。また、ウォールフィルタ処理機能121A,121Bを区別無く総称する場合には、「ウォールフィルタ処理機能121」と記載する。 The wall filter processing functions 121A and 121B perform filter processing on the ultrasonic reception signal (data sequence of reflected wave data) received in time series for each sample point. Here, the wall filter processing function 121A performs filter processing A having a cutoff value A on the data sequence of the reflected wave signal. Also, the wall filter processing function 121B performs filter processing B having a cutoff value B smaller than the cutoff value A on the data sequence of the reflected wave signal. Note that the wall filter processing function 121A is an example of a first filter processing unit. Also, the wall filter processing function 121B is an example of a second filter processing unit. Also, the cutoff value A is an example of a first cutoff value. Also, the cutoff value B is an example of a second cutoff value. Also, the filter processing A is an example of a first filter processing. Also, the filter processing B is an example of a second filter processing. Also, when the wall filter processing functions 121A and 121B are collectively referred to without distinction, they are described as "wall filter processing function 121".

例えば、ウォールフィルタ処理機能121は、フィルタ処理として、ハイパスフィルタ(例えばIIRフィルタ)を適用する。ウォールフィルタ処理機能121は、第nフレームにおける所定位置(サンプル点)のIIRフィルタ出力データ(ドプラ信号)を得るために、同一位置における、第nフレームの反射波データと、過去の4フレーム分(第(n-4)フレーム~第(n-1)フレーム)の反射波データと、過去の4フレーム分のIIRフィルタ出力データ(ドプラ信号)とを用いる。これらの反射波データは、上述したように、1フレームの走査範囲(第1走査範囲)を形成する複数の走査線それぞれで、1回ずつ超音波送受信を行うことで生成された反射波データである。なお、ウォールフィルタ処理機能121は、ハイパスフィルタに限定されるものではない。例えば、ウォールフィルタ処理機能121Bは、バンドパスフィルタを適用しても良い。 For example, the wall filter processing function 121 applies a high-pass filter (e.g., an IIR filter) as the filter processing. To obtain IIR filter output data (Doppler signal) at a predetermined position (sample point) in the nth frame, the wall filter processing function 121 uses the reflected wave data of the nth frame at the same position, the reflected wave data of the past four frames (frame (n-4) to frame (n-1)), and the IIR filter output data (Doppler signal) of the past four frames. As described above, these reflected wave data are reflected wave data generated by transmitting and receiving ultrasound once for each of the multiple scanning lines that form the scanning range of one frame (first scanning range). Note that the wall filter processing function 121 is not limited to a high-pass filter. For example, the wall filter processing function 121B may apply a band-pass filter.

そして、信号処理回路120は、フィルタ処理Aに基づくドプラ情報Aと、フィルタ処理Bに基づくドプラ情報Bとを算出する(ステップS103)。ステップS103の処理は、自己相関演算機能122A,122B、及びドプラ情報演算機能123A,123Bによって実行される。 Then, the signal processing circuit 120 calculates Doppler information A based on the filter processing A and Doppler information B based on the filter processing B (step S103). The process of step S103 is executed by the autocorrelation calculation functions 122A and 122B and the Doppler information calculation functions 123A and 123B.

自己相関演算機能122A,122Bは、最新フレームのドプラ信号のIQ信号と、1フレーム前のドプラ信号のIQ信号との複素共役をとることで自己相関値を算出する。具体的には、自己相関演算機能122Aは、フィルタ処理Aによって出力されたドプラ信号を用いて自己相関値を算出する。また、自己相関演算機能122Bは、フィルタ処理Bによって出力されたドプラ信号を用いて自己相関値を算出する。なお、自己相関演算機能122A,122Bを区別無く総称する場合には、「自己相関演算機能122」と記載する。 The autocorrelation calculation functions 122A and 122B calculate the autocorrelation value by taking the complex conjugate of the IQ signal of the Doppler signal of the latest frame and the IQ signal of the Doppler signal of the previous frame. Specifically, the autocorrelation calculation function 122A calculates the autocorrelation value using the Doppler signal output by filter processing A. The autocorrelation calculation function 122B calculates the autocorrelation value using the Doppler signal output by filter processing B. Note that when the autocorrelation calculation functions 122A and 122B are referred to collectively without distinction, they are referred to as "autocorrelation calculation function 122."

ドプラ情報演算機能123A,123Bは、自己相関演算機能122が算出した自己相関値から、平均速度、平均分散値、平均パワー値等のドプラ情報を算出する。具体的には、ドプラ情報演算機能123Aは、自己相関演算機能122Aが算出した自己相関値からドプラ情報Aを算出する。また、ドプラ情報演算機能123Bは、自己相関演算機能122Bが算出した自己相関値からドプラ情報Bを算出する。なお、ドプラ情報Aは、第1ドプラ情報の一例である。また、ドプラ情報Bは、第2ドプラ情報の一例である。 The Doppler information calculation functions 123A and 123B calculate Doppler information such as average velocity, average variance, and average power from the autocorrelation values calculated by the autocorrelation calculation function 122. Specifically, the Doppler information calculation function 123A calculates Doppler information A from the autocorrelation values calculated by the autocorrelation calculation function 122A. Furthermore, the Doppler information calculation function 123B calculates Doppler information B from the autocorrelation values calculated by the autocorrelation calculation function 122B. Note that the Doppler information A is an example of first Doppler information. Furthermore, the Doppler information B is an example of second Doppler information.

そして、信号処理回路120は、ドプラ情報Bに基づいて、閾値を設定する(ステップS104)。そして、信号処理回路120は、ドプラ情報Aが閾値未満か否かを判定する(ステップS105)。そして、信号処理回路120は、ドプラ情報Aが閾値未満である場合には(ステップS105,Yes)、血流情報として略ゼロを採用する(ステップS106)。一方、信号処理回路120は、ドプラ情報Aが閾値以上である場合には(ステップS105,No)、血流情報としてドプラ情報Aを採用する(ステップS107)。ステップS104~ステップS107の各処理は、ブランク処理機能124によって実行される。 Then, the signal processing circuit 120 sets a threshold value based on the Doppler information B (step S104). The signal processing circuit 120 then determines whether or not the Doppler information A is less than the threshold value (step S105). If the Doppler information A is less than the threshold value (step S105, Yes), the signal processing circuit 120 adopts approximately zero as the blood flow information (step S106). On the other hand, if the Doppler information A is equal to or greater than the threshold value (step S105, No), the signal processing circuit 120 adopts the Doppler information A as the blood flow information (step S107). Each process from step S104 to step S107 is executed by the blank processing function 124.

ブランク処理機能124は、ドプラ情報Aからドプラ情報Bを部分的に除くブランク処理を実行することで、血流情報を生成する。例えば、ブランク処理機能124は、フィルタ処理Aによって抽出された各サンプル点のドプラ情報Aと、フィルタ処理Bによって抽出された各サンプル点のドプラ情報Bとに基づいて、クラッタ信号成分が除去された血流情報を生成する。なお、ブランク処理機能124は、血流情報生成部の一例である。 The blanking function 124 generates blood flow information by performing blanking processing to partially remove Doppler information B from Doppler information A. For example, the blanking function 124 generates blood flow information from which clutter signal components have been removed, based on the Doppler information A of each sample point extracted by filter processing A and the Doppler information B of each sample point extracted by filter processing B. The blanking function 124 is an example of a blood flow information generating unit.

ここで、図4A、図4B、図5A、及び図5Bを用いて、ブランク処理機能124の処理について説明する。図4A~図5Bは、第1の実施形態に係るブランク処理機能124の処理について説明するための図である。図4A~図5Bにおいて、縦軸はドプラ情報に含まれる各信号成分の強度に対応し、横軸は周波数に対応する。 The processing of the blanking function 124 will now be described with reference to Figures 4A, 4B, 5A, and 5B. Figures 4A to 5B are diagrams for explaining the processing of the blanking function 124 according to the first embodiment. In Figures 4A to 5B, the vertical axis corresponds to the intensity of each signal component included in the Doppler information, and the horizontal axis corresponds to the frequency.

例えば、図4Aに示すように、心臓の動きの大きさに対して十分に血流移動が大きい場合には、受信信号に対して高カットオフ値のハイパスフィルタA(フィルタ処理A)を適用すると、平均速度「Vel-B1」のドプラ情報が抽出される。また、同じ受信信号に対して低カットオフ値のハイパスフィルタB(フィルタ処理B)を適用すると、平均速度「Vel-C1」のドプラ情報が抽出される。ここで、「Vel-B1」は、血流信号が支配的であるため高い値となる。一方、「Vel-C1」は、ハイパスフィルタAの場合と比較して多くのクラッタ信号が残存しているため、クラッタ信号の残存量に依存して低い値となる。この結果、「Vel-B1」と「Vel-C1」の速度値の差は大きくなる。 For example, as shown in FIG. 4A, when the blood flow movement is sufficiently large relative to the magnitude of the heart's movement, applying high-pass filter A (filter processing A) with a high cutoff value to the received signal extracts Doppler information of the average velocity "Vel-B1". Furthermore, applying high-pass filter B (filter processing B) with a low cutoff value to the same received signal extracts Doppler information of the average velocity "Vel-C1". Here, "Vel-B1" is a high value because the blood flow signal is dominant. On the other hand, "Vel-C1" is a low value depending on the amount of remaining clutter signal because more clutter signal remains compared to the case of high-pass filter A. As a result, the difference in the velocity values between "Vel-B1" and "Vel-C1" becomes large.

一方、図4Bに示すように、心臓の動きの大きさに対して血流移動が小さい場合には、受信信号に対して高カットオフ値のハイパスフィルタA(フィルタ処理A)を適用すると、平均速度「Vel-B2」のドプラ情報が抽出される。また、同じ受信信号に対して低カットオフ値のハイパスフィルタB(フィルタ処理B)を適用すると、平均速度「Vel-C2」のドプラ情報が抽出される。ここで、「Vel-B2」及び「Vel-C2」は、いずれもクラッタ信号が支配的である。このため、「Vel-B2」と「Vel-C2」の速度値の差は、図4Aの場合と比較して小さい。 On the other hand, as shown in Figure 4B, when blood flow movement is small relative to the magnitude of cardiac movement, applying high-pass filter A (filter processing A) with a high cutoff value to the received signal extracts Doppler information of average velocity "Vel-B2". Furthermore, applying high-pass filter B (filter processing B) with a low cutoff value to the same received signal extracts Doppler information of average velocity "Vel-C2". Here, clutter signals dominate both "Vel-B2" and "Vel-C2". For this reason, the difference in velocity values between "Vel-B2" and "Vel-C2" is smaller than in the case of Figure 4A.

このように、クラッタ信号に対する血流信号の大きさに依存して、カットオフ値が異なる2つのウォールフィルタから得られる2つのドプラ情報(速度値)の差には違いがある。そこで、ブランク処理機能124は、この2つのドプラ情報(速度値)の差の違いを利用することで、クラッタ信号成分を適切に除去した血流情報を生成する。 As such, depending on the magnitude of the blood flow signal relative to the clutter signal, there is a difference in the difference between the two pieces of Doppler information (velocity values) obtained from two wall filters with different cutoff values. Therefore, the blanking processing function 124 uses the difference between these two pieces of Doppler information (velocity values) to generate blood flow information from which the clutter signal components have been appropriately removed.

具体例を挙げると、ブランク処理機能124は、ドプラ情報Bに基づいて、ブランク処理で用いる閾値を設定する。例えば、ブランク処理機能124は、ドプラ情報演算機能123Bによって算出された平均速度に対して、予め設定された固定値「α」を加算した値を、閾値として設定する。そして、ブランク処理機能124は、設定した閾値を用いた閾値処理により、血流情報を生成する。 As a specific example, the blanking function 124 sets a threshold value to be used in blanking based on the Doppler information B. For example, the blanking function 124 sets a threshold value obtained by adding a preset fixed value "α" to the average velocity calculated by the Doppler information calculation function 123B. Then, the blanking function 124 generates blood flow information by threshold processing using the set threshold value.

ここで、固定値「α」は、操作者によって任意に設定可能な値である。例えば、固定値「α」は、スキャン部位ごとに予め設定されるのが好適である。また、固定値「α」は、「Vel-B1」と「Vel-C1」の速度値の差(図4Aの速度値の差)より小さく、「Vel-B2」と「Vel-C2」の速度値の差(図4Bの速度値の差)より大きい値であるのが好適である。 The fixed value "α" here is a value that can be arbitrarily set by the operator. For example, it is preferable that the fixed value "α" be set in advance for each scan area. It is also preferable that the fixed value "α" is a value that is smaller than the difference in the velocity values between "Vel-B1" and "Vel-C1" (the difference in the velocity values in FIG. 4A) and larger than the difference in the velocity values between "Vel-B2" and "Vel-C2" (the difference in the velocity values in FIG. 4B).

例えば、クラッタ信号に対して血流信号が十分に大きい場合には、ブランク処理機能124は、図4Aのドプラ情報Bに対応する平均速度「Vel-C1」に対して固定値「α」を加算することで、閾値「TH-1」を設定する(図5A)。ここで、固定値「α」は、「Vel-B1」と「Vel-C1」の速度値の差より小さい値であるので、閾値「TH-1」は、平均速度「Vel-B1」より小さくなる。 For example, when the blood flow signal is sufficiently large relative to the clutter signal, the blanking function 124 sets the threshold value "TH-1" by adding a fixed value "α" to the average velocity "Vel-C1" corresponding to the Doppler information B in FIG. 4A (FIG. 5A). Here, the fixed value "α" is a value smaller than the difference in the velocity values between "Vel-B1" and "Vel-C1", so the threshold value "TH-1" is smaller than the average velocity "Vel-B1".

一方、クラッタ信号に対して血流信号が小さい場合には、ブランク処理機能124は、図4Bのドプラ情報Bに対応する平均速度「Vel-C2」に対して固定値「α」を加算することで、閾値「TH-2」を設定する(図5B)。ここで、固定値「α」は、「Vel-B2」と「Vel-C2」の速度値の差より大きい値であるので、閾値「TH-2」は、平均速度「Vel-B2」より大きくなる。 On the other hand, when the blood flow signal is small relative to the clutter signal, the blanking function 124 sets the threshold value "TH-2" by adding a fixed value "α" to the average velocity "Vel-C2" corresponding to the Doppler information B in FIG. 4B (FIG. 5B). Here, the fixed value "α" is a value larger than the difference in the velocity values between "Vel-B2" and "Vel-C2", so the threshold value "TH-2" is larger than the average velocity "Vel-B2".

換言すると、受信信号に含まれるクラッタ信号や血流信号の程度が不明であったとしても、ブランク処理機能124は、2つのウォールフィルタから得られる2つのドプラ情報(速度値)の差の大小を、閾値を用いて判定する。この結果、ブランク処理機能124は、受信信号のうち、クラッタ信号に対する血流信号の大きさを推定することができる。 In other words, even if the level of clutter signals or blood flow signals contained in the received signal is unknown, the blanking function 124 uses a threshold value to determine the magnitude of the difference between the two pieces of Doppler information (velocity values) obtained from the two wall filters. As a result, the blanking function 124 can estimate the magnitude of the blood flow signal relative to the clutter signal in the received signal.

つまり、ブランク処理機能124は、フレームごとに、各サンプル点のドプラ情報Aと閾値とを比較する。そして、ブランク処理機能124は、ドプラ情報Aが閾値未満である場合には、クラッタ信号に対して血流信号が小さく、ドプラ情報に含まれる成分の多くがクラッタ信号であると推定できる。このため、ブランク処理機能124は、ドプラ情報Aが閾値未満であるサンプル点については、血流情報として略ゼロを採用する。一方、ブランク処理機能124は、ドプラ情報Aが閾値以上である場合には、クラッタ信号に対して血流信号が大きく、ドプラ情報に含まれる成分の多くが血流信号であると推定できる。このため、ブランク処理機能124は、ドプラ情報Aが閾値以上であるサンプル点については、血流情報としてドプラ情報Aを採用する。 That is, the blanking function 124 compares the Doppler information A of each sample point with the threshold value for each frame. Then, when the Doppler information A is less than the threshold value, the blanking function 124 can estimate that the blood flow signal is small relative to the clutter signal, and that most of the components included in the Doppler information are clutter signals. For this reason, the blanking function 124 adopts approximately zero as the blood flow information for sample points where the Doppler information A is less than the threshold value. On the other hand, when the Doppler information A is equal to or greater than the threshold value, the blanking function 124 can estimate that the blood flow signal is large relative to the clutter signal, and that most of the components included in the Doppler information are blood flow signals. For this reason, the blanking function 124 adopts the Doppler information A as the blood flow information for sample points where the Doppler information A is equal to or greater than the threshold value.

このように、ブランク処理機能124は、ドプラ情報Aと、ドプラ情報Bに基づく閾値との比較に基づいて、血流情報を生成する。この結果、ブランク処理機能124は、フレームごとに、各サンプル点について適切に閾値が設定されたブランク処理を実行できるので、個人差が大きい生体においてもクラッタ信号成分を適切に除去した血流情報を生成することができる。 In this way, the blanking function 124 generates blood flow information based on a comparison between Doppler information A and a threshold based on Doppler information B. As a result, the blanking function 124 can perform blanking processing with an appropriate threshold set for each sample point for each frame, so that blood flow information from which clutter signal components have been appropriately removed can be generated even for living organisms with large individual differences.

なお、図4A~図5Bの説明では、ブランク処理機能124が、ドプラ情報として平均速度を用いてブランク処理を実行する場合を説明したが、実施形態はこれに限定されるものではない。例えば、ブランク処理機能124が、ドプラ情報として、平均分散値や平均パワー値等を用いてブランク処理を実行することも可能である。また、ブランク処理機能124は、ドプラ情報のうち1種類のパラメータを用いるのではなく、複数種類のパラメータを組み合わせてブランク処理を実行することも可能である。組み合わせの例としては、平均速度と平均分散値、平均速度と平均パワー値、平均分散値と平均パワー値、及び平均速度と平均分散値と平均パワー値の4通りが適用可能である。 Note that in the explanation of Figures 4A to 5B, the blanking function 124 performs blanking using average velocity as Doppler information, but the embodiment is not limited to this. For example, the blanking function 124 can also perform blanking using an average variance value, average power value, or the like as Doppler information. Furthermore, the blanking function 124 can also perform blanking by combining multiple types of parameters rather than using one type of parameter from the Doppler information. As examples of combinations, four combinations are applicable: average velocity and average variance value, average velocity and average power value, average variance value and average power value, and average velocity, average variance value and average power value.

また、ブランク処理は、上述した閾値処理に限定されるものではない。例えば、ブランク処理機能124は、ブランク処理として、差分処理を実行することも可能である。差分処理については、変形例にて後述する。 The blanking process is not limited to the threshold process described above. For example, the blanking process function 124 can also perform differential processing as blanking process. The differential processing will be described later in the modified example.

また、「略ゼロ」という記載は、必ずしも「ゼロ」に限定されるものではなく、僅かな値を有していても良いことを意図したものである。例えば、ドプラ画像データとして画像化される場合、閲覧者が画像上でゼロと区別できない程度の僅かな値であれば、「ゼロ」に代えて採用可能である。 The phrase "almost zero" does not necessarily mean "zero" and is intended to mean that it may have a small value. For example, when visualized as Doppler image data, any small value that is indistinguishable from zero by the viewer on the image can be used in place of "zero."

図3の説明に戻る。信号処理回路120は、平滑化処理を実行する(ステップS108)。例えば、パーシスタンス機能125は、ブランク処理機能124から出力された血流情報に対して、時間方向に平滑化するフィルタ処理を適用する。パーシスタンス機能125は、平滑化処理後の血流情報を、例えば画像生成回路130に出力する。 Returning to the explanation of FIG. 3, the signal processing circuit 120 executes smoothing processing (step S108). For example, the persistence function 125 applies a filter process to the blood flow information output from the blanking processing function 124 to smooth the blood flow information in the time direction. The persistence function 125 outputs the blood flow information after the smoothing process to, for example, the image generation circuit 130.

そして、画像生成回路130は、ドプラ画像データを生成する(ステップS109)。例えば、画像生成回路130は、パーシスタンス機能125から出力された平滑化処理後の血流情報に対してスキャンコンバート等の処理を行って、血流情報に応じた画素値が割り当てられたドプラ画像データを生成する。生成されたドプラ画像データは、画像メモリ140に格納され、必要に応じてディスプレイ103に表示されたり、超音波診断装置1に接続される外部装置に転送されたりする。 Then, the image generation circuit 130 generates Doppler image data (step S109). For example, the image generation circuit 130 performs processing such as scan conversion on the smoothed blood flow information output from the persistence function 125 to generate Doppler image data to which pixel values according to the blood flow information are assigned. The generated Doppler image data is stored in the image memory 140 and, as necessary, is displayed on the display 103 or transferred to an external device connected to the ultrasound diagnostic device 1.

このように、超音波診断装置1は、ステップS101~ステップS109の各処理を実行することで、各フレームのドプラ画像データを動画像として表示する。そして、超音波診断装置1は、ドプラ画像データの表示を終了する旨の終了指示を操作者から受け付けると、ステップS101~ステップS109の各処理の実行を終了する。 In this way, the ultrasound diagnostic device 1 executes each process from step S101 to step S109 to display the Doppler image data of each frame as a moving image. Then, when the ultrasound diagnostic device 1 receives an end instruction from the operator to end the display of the Doppler image data, it ends the execution of each process from step S101 to step S109.

なお、図3にて説明した処理手順の順序はあくまで一例であり、必ずしも図示の順序で実行されるものではない。例えば、図3の処理手順の順序は、処理内容に矛盾が生じない範囲内で任意に変更可能である。 Note that the order of the processing steps described in FIG. 3 is merely an example, and the processing steps are not necessarily executed in the order shown. For example, the order of the processing steps in FIG. 3 can be changed as desired as long as no inconsistencies occur in the processing content.

上述してきたように、第1の実施形態に係る超音波診断装置1は、互いに異なるカットオフ値を有する複数のフィルタ処理を超音波の受信信号に対してそれぞれ並列に適用することによって抽出された複数のドプラ情報に基づいて、サンプル点ごとにクラッタ信号成分が除去された血流情報を生成する。これによれば、超音波診断装置1は、クラッタ信号成分を適切に除去することが可能となる。 As described above, the ultrasound diagnostic device 1 according to the first embodiment generates blood flow information from which clutter signal components have been removed for each sample point, based on multiple pieces of Doppler information extracted by applying multiple filter processes having different cutoff values in parallel to the received ultrasound signal. This enables the ultrasound diagnostic device 1 to appropriately remove the clutter signal components.

図6A及び図6Bを用いて、第1の実施形態に係る超音波診断装置1による効果を説明する。図6A及び図6Bは、第1の実施形態に係る超音波診断装置1による効果を説明するための図である。図6Aには、第1の実施形態に係る超音波診断装置1により得られるドプラ画像データの一例を示す。また、図6Bには、第1の実施形態が適用されない比較例に係る超音波診断装置により得られるドプラ画像データの一例を示す。 The effects of the ultrasound diagnostic device 1 according to the first embodiment will be described using Figures 6A and 6B. Figures 6A and 6B are diagrams for explaining the effects of the ultrasound diagnostic device 1 according to the first embodiment. Figure 6A shows an example of Doppler image data obtained by the ultrasound diagnostic device 1 according to the first embodiment. Also, Figure 6B shows an example of Doppler image data obtained by an ultrasound diagnostic device according to a comparative example to which the first embodiment is not applied.

図6A及び図6Bに示すように、第1の実施形態に係る超音波診断装置1により得られるドプラ画像データでは、関心領域の中央付近に冠動脈血流が明瞭に描出される(図6A)。これに対して、比較例に係る超音波診断装置により得られるドプラ画像データでは、冠動脈血流は描出されず、クラッタ信号のみが描出される。 As shown in Figures 6A and 6B, in the Doppler image data obtained by the ultrasound diagnostic device 1 according to the first embodiment, the coronary artery blood flow is clearly depicted near the center of the region of interest (Figure 6A). In contrast, in the Doppler image data obtained by the ultrasound diagnostic device according to the comparative example, the coronary artery blood flow is not depicted, and only clutter signals are depicted.

ここで、比較例における課題を説明する。第1の実施形態が適用されない比較例としては、大きく下記の2つの手法が考えられるが、いずれの手法においてもクラッタ信号成分を適切に除去することが難しい。 Here, we will explain the problems with the comparative example. The following two methods are broadly considered as comparative examples to which the first embodiment is not applied, but in either method, it is difficult to appropriately remove clutter signal components.

第1の手法は、2つのフィルタ処理を順番に適用する手法である。具体的には、第1の手法では、まず、ドプラ情報に対してローパスフィルタを適用してクラッタ信号成分を積極的に抽出する。そして、抽出したクラッタ信号成分に基づいて必要に応じて元のドプラ情報に対してフィードバック処理をかける。そして、フィードバック処理後のドプラ情報に対してハイパスフィルタを適用することで、血流信号成分とクラッタ信号成分とを分離する。 The first method is a method of applying two filter processes in sequence. Specifically, in the first method, a low-pass filter is first applied to the Doppler information to actively extract clutter signal components. Then, feedback processing is applied to the original Doppler information as necessary based on the extracted clutter signal components. Then, a high-pass filter is applied to the Doppler information after feedback processing to separate the blood flow signal components and the clutter signal components.

しかしながら、第1の手法では、元のドプラ情報に対してローパスフィルタとハイパスフィルタを順に適用するため、フィルタ処理に要する処理時間が2回分かかる上、フィードバック処理にも一定の処理時間がかかっていた。また、フィードバック処理で用いられる係数は経験的に設定せざるを得ず、必ずしも適切な係数を適用できるとは限らない。このように、第1の手法では、クラッタ信号成分を適切に除去することが難しい。 However, in the first method, a low-pass filter and a high-pass filter are applied to the original Doppler information in sequence, so the processing time required for the filter processing is required twice, and the feedback processing also requires a certain amount of processing time. In addition, the coefficients used in the feedback processing must be set empirically, and it is not always possible to apply appropriate coefficients. As such, it is difficult to appropriately remove clutter signal components with the first method.

第2の手法は、ブランク処理を行う手法である。具体的には、第2の手法では、ドプラ情報の速度と閾値とを比較し、閾値未満の速度を略ゼロとすることで、クラッタ信号成分が除去される。 The second method is a method of performing blanking processing. Specifically, in the second method, the velocity of the Doppler information is compared with a threshold, and any velocity below the threshold is set to approximately zero, thereby removing the clutter signal components.

しかしながら、第2の手法では、全フレームの全サンプル点に対して単一の閾値を適用するため、必ずしもクラッタ信号成分を適切に除去することができなかった。例えば、冠動脈血流のように、クラッタ信号成分が大きい撮像環境下では、必ずしもクラッタ信号成分を適切に除去することができない。 However, in the second method, a single threshold is applied to all sample points in all frames, so clutter signal components cannot always be removed appropriately. For example, in an imaging environment with large clutter signal components, such as coronary artery blood flow, the clutter signal components cannot always be removed appropriately.

図7A、図7B、及び図8を用いて、比較例に係る第2の手法における課題を説明する。図7A、図7B、及び図8は、比較例に係る第2の手法における課題を説明するための図である。図7A、図7B、及び図8において、縦軸は信号強度に対応し、横軸は周波数に対応する。 The problems with the second method according to the comparative example will be explained using Figures 7A, 7B, and 8. Figures 7A, 7B, and 8 are diagrams for explaining the problems with the second method according to the comparative example. In Figures 7A, 7B, and 8, the vertical axis corresponds to signal strength, and the horizontal axis corresponds to frequency.

例えば、図7Aに示すように、心臓の動きの大きさに対して血流移動が小さい場合には、クラッタ信号が支配的となる。この場合、ウォールフィルタ(ハイパスフィルタ)を適用すると、フィルタ処理後の信号成分もクラッタ信号が支配的となり、平均速度は「Vel-C0」となる。一方、図7Bに示すように、心臓の動きの大きさに対して十分に血流移動が大きい場合にウォールフィルタを適用すると、フィルタ処理後の信号成分は血流信号が支配的となり、平均速度は「Vel-B0」となる。 For example, as shown in Figure 7A, when blood flow movement is small relative to the magnitude of cardiac movement, clutter signals become dominant. In this case, when a wall filter (high-pass filter) is applied, the signal components after filtering are also dominated by clutter signals, and the average velocity becomes "Vel-C0". On the other hand, as shown in Figure 7B, when a wall filter is applied when blood flow movement is sufficiently large relative to the magnitude of cardiac movement, the signal components after filtering are dominated by blood flow signals, and the average velocity becomes "Vel-B0".

ここで、例えば、平均速度「Vel-C0」が平均速度「Vel-B0」より大きい値となる場合は、平均速度「Vel-C0」と平均速度「Vel-B0」とが近い値となる場合には、単一の閾値を用いるブランク処理では、クラッタ信号成分を適切に除去することができない。 Here, for example, if the average speed "Vel-C0" is greater than the average speed "Vel-B0", and the average speeds "Vel-C0" and "Vel-B0" are close to each other, the clutter signal components cannot be adequately removed by blanking processing using a single threshold.

例えば、図8に示すように、閾値が「Th-A」に設定される場合には、ブランク処理を行っても血流信号もクラッタ信号も残存することとなる。また、閾値が「Th-B」に設定される場合には、ブランク処理によって血流信号もクラッタ信号も除去されることとなる。このように、単一の閾値を用いるブランク処理では、クラッタ信号成分を適切に除去することができない。 For example, as shown in FIG. 8, when the threshold is set to "Th-A," both the blood flow signal and the clutter signal remain even after blanking. When the threshold is set to "Th-B," both the blood flow signal and the clutter signal are removed by blanking. In this way, blanking using a single threshold is unable to adequately remove clutter signal components.

これに対して、第1の実施形態に係る超音波診断装置1は、互いに異なるカットオフ値を有する2つのウォールフィルタを同一の受信信号に対して並列に適用する。このため、超音波診断装置1は、実質1回分の処理時間でウォールフィルタを適用可能である。 In contrast, the ultrasound diagnostic device 1 according to the first embodiment applies two wall filters having different cutoff values in parallel to the same received signal. Therefore, the ultrasound diagnostic device 1 can apply the wall filters in substantially one processing time.

また、超音波診断装置1は、フレームごとに、各サンプル点の受信信号から得られるドプラ情報Bに基づいて、ブランク処理に用いる閾値を設定する。このため、超音波診断装置1は、フレームごとに、各サンプル点の受信信号に含まれるクラッタ信号及び血流信号の割合に応じて適切にブランク処理を実行することが可能である。 In addition, the ultrasound diagnostic device 1 sets a threshold value to be used for blanking processing for each frame based on the Doppler information B obtained from the received signal at each sample point. Therefore, the ultrasound diagnostic device 1 can appropriately perform blanking processing for each frame according to the proportion of clutter signals and blood flow signals contained in the received signal at each sample point.

(第1の実施形態の変形例)
上記の実施形態では、ブランク処理として閾値処理が適用される場合を説明したが、実施形態はこれに限定されるものではない。例えば、ブランク処理機能124は、ブランク処理として、差分処理を実行することも可能である。すなわち、ブランク処理機能124は、ドプラ情報Aから、ドプラ情報Bに基づく減数値を減算する減算処理に基づいて、血流情報を生成する。
(Modification of the first embodiment)
In the above embodiment, a case where threshold processing is applied as blanking processing has been described, but the embodiment is not limited to this. For example, the blanking function 124 can also execute subtraction processing as blanking processing. That is, the blanking function 124 generates blood flow information based on a subtraction process that subtracts a subtraction value based on Doppler information B from Doppler information A.

図9は、第1の実施形態の変形例に係る超音波診断装置1の処理手順の一例を示すフローチャートである。図9に示すステップS201~ステップS210の各処理は、例えば、ドプラ画像データの表示を開始する旨の開始指示を操作者から受け付けた場合に開始される。 Figure 9 is a flowchart showing an example of a processing procedure of the ultrasound diagnostic device 1 according to a modified example of the first embodiment. Each process of steps S201 to S210 shown in Figure 9 is started, for example, when a start instruction to start displaying Doppler image data is received from the operator.

なお、図9の処理手順のうち、ステップS201~ステップS203の各処理は、図3に示したステップS101~ステップS103の各処理と同様であるので説明を省略する。また、ステップS209~ステップS210の各処理は、図3に示したステップS108~ステップS109の各処理と同様であるので説明を省略する。 Note that among the processing steps in FIG. 9, steps S201 to S203 are similar to steps S101 to S103 shown in FIG. 3, and therefore their explanations are omitted. Also, steps S209 to S210 are similar to steps S108 to S109 shown in FIG. 3, and therefore their explanations are omitted.

図9に示すように、信号処理回路120は、ドプラ情報Bに基づいて、減数値を設定する(ステップS204)。そして、信号処理回路120は、ドプラ情報Aから減数値を減算する(ステップS205)。そして、信号処理回路120は、減算後の値がゼロ未満か否かを判定する(ステップS206)。そして、信号処理回路120は、減算後の値がゼロ未満である場合には(ステップS206,Yes)、血流情報として略ゼロを採用する(ステップS207)。一方、信号処理回路120は、減算後の値がゼロ以上である場合には(ステップS206,Yes)、血流情報としてドプラ情報Aを採用する(ステップS208)。 As shown in FIG. 9, the signal processing circuit 120 sets a subtraction value based on Doppler information B (step S204). Then, the signal processing circuit 120 subtracts the subtraction value from Doppler information A (step S205). Then, the signal processing circuit 120 determines whether the value after subtraction is less than zero (step S206). Then, if the value after subtraction is less than zero (step S206, Yes), the signal processing circuit 120 adopts approximately zero as the blood flow information (step S207). On the other hand, if the value after subtraction is equal to or greater than zero (step S206, Yes), the signal processing circuit 120 adopts Doppler information A as the blood flow information (step S208).

例えば、ブランク処理機能124は、ドプラ情報Bに対して予め設定された固定値「α」を加算した値を、減数値として設定する。なお、固定値「α」は、上記の実施形態にて説明した固定値「α」と同様であるので、説明を省略する。 For example, the blanking function 124 sets the subtraction value to a value obtained by adding a preset fixed value "α" to the Doppler information B. Note that the fixed value "α" is the same as the fixed value "α" described in the above embodiment, and therefore will not be described here.

そして、ブランク処理機能124は、ドプラ情報Aから減数値を減算し、減算処理の差分値がゼロ未満であるか否かを判定する。ブランク処理機能124は、減算処理の差分値がゼロ未満であるサンプル点については、血流情報として略ゼロを採用する。一方、ブランク処理機能124は、減算処理の差分値がゼロ以上であるサンプル点については、血流情報としてドプラ情報Aを採用する。 Then, the blanking function 124 subtracts the subtraction value from the Doppler information A and determines whether the difference value of the subtraction process is less than zero. For sample points where the difference value of the subtraction process is less than zero, the blanking function 124 adopts approximately zero as the blood flow information. On the other hand, for sample points where the difference value of the subtraction process is zero or greater, the blanking function 124 adopts the Doppler information A as the blood flow information.

このように、ブランク処理機能124は、ドプラ情報Aから、ドプラ情報Bに基づく減数値を減算する減算処理に基づいて、血流情報を生成する。これにより、ブランク処理機能124は、上記の実施形態にて説明した閾値処理によるブランク処理と同様の処理を行うことができる。 In this way, the blanking function 124 generates blood flow information based on a subtraction process that subtracts a subtraction value based on Doppler information B from Doppler information A. This allows the blanking function 124 to perform processing similar to the blanking process using threshold processing described in the above embodiment.

(第2の実施形態)
上記の実施形態では、2つのウォールフィルタを用いる処理を説明したが、実施形態はこれに限定されるものではない。例えば、本実施形態は、3つ以上のウォールフィルタを用いることも可能である。
Second Embodiment
In the above embodiment, the process using two wall filters has been described, but the embodiment is not limited to this. For example, the present embodiment can also use three or more wall filters.

図10を用いて、第2の実施形態に係る信号処理回路120の処理について説明する。図10は、第2の実施形態に係る信号処理回路120の処理について説明するための図である。図10において、縦軸はドプラ情報に含まれる各信号成分の強度に対応し、横軸は周波数に対応する。 The processing of the signal processing circuit 120 according to the second embodiment will be described with reference to FIG. 10. FIG. 10 is a diagram for explaining the processing of the signal processing circuit 120 according to the second embodiment. In FIG. 10, the vertical axis corresponds to the intensity of each signal component included in the Doppler information, and the horizontal axis corresponds to the frequency.

図10に示すように、血流信号に対してクラッタ信号が大きく無い場合がある。このような場合には、高カットオフ値のハイパスフィルタAにより得られるドプラ情報Aと、低カットオフ値(図10の「低カットオフ値1」)のハイパスフィルタBにより得られるドプラ情報Bは、いずれも血流信号が支配的となる。このため、ドプラ情報Aの平均速度「Vel-B3」と、ドプラ情報Bの平均速度「Vel-C3」との速度値の差は、小さくなる。このため、第1の実施形態に係るブランク処理では、血流信号が支配的であるにもかかわらずドプラ情報Aを除去する(略ゼロにする)ように機能してしまう場合がある。 As shown in FIG. 10, there are cases where the clutter signal is not large relative to the blood flow signal. In such cases, the blood flow signal dominates in both Doppler information A obtained by high-pass filter A with a high cutoff value and Doppler information B obtained by high-pass filter B with a low cutoff value ("low cutoff value 1" in FIG. 10). As a result, the difference in velocity value between the average velocity "Vel-B3" of Doppler information A and the average velocity "Vel-C3" of Doppler information B becomes small. For this reason, the blanking process according to the first embodiment may function to remove Doppler information A (to reduce it to approximately zero) even though the blood flow signal is dominant.

そこで、血流信号に対してクラッタ信号が大きく無いスキャン部位では、3つ目のウォールフィルタを用いるのが好適である。例えば、信号処理回路120は、3つ目のウォールフィルタとして、低カットオフ値1より小さい低カットオフ値2のハイパスフィルタCを適用する。ハイパスフィルタCにより得られるドプラ情報Cの平均速度「Vel-C4」は、クラッタ信号が大きければクラッタ信号に依存してより低い値となるが、クラッタ信号が小さければ平均速度「Vel-B3」や平均速度「Vel-C3」と同程度の値となる。そこで、信号処理回路120は、平均速度「Vel-C4」と、平均速度「Vel-B3」及び平均速度「Vel-C3」のうち少なくとも一方とを比較することで、受信信号に含まれるクラッタ信号及び血流信号の割合を推定する。 Therefore, it is preferable to use a third wall filter in a scan area where the clutter signal is not large relative to the blood flow signal. For example, the signal processing circuit 120 applies a high-pass filter C with a low cutoff value 2, which is smaller than the low cutoff value 1, as the third wall filter. The average velocity "Vel-C4" of the Doppler information C obtained by the high-pass filter C becomes a lower value depending on the clutter signal if the clutter signal is large, but becomes a value similar to the average velocity "Vel-B3" or the average velocity "Vel-C3" if the clutter signal is small. Therefore, the signal processing circuit 120 estimates the ratio of the clutter signal and the blood flow signal contained in the received signal by comparing the average velocity "Vel-C4" with at least one of the average velocity "Vel-B3" and the average velocity "Vel-C3".

図11は、第2の実施形態に係る信号処理回路120の構成例を示すブロック図である。図11に示すように、信号処理回路120は、図2に示した構成に加えて、ウォールフィルタ処理機能121C、自己相関演算機能122C、及びドプラ情報演算機能123Cを備える。なお、第2の実施形態に係る信号処理回路120が実行する各処理機能については、図12の説明にて後述する。 Fig. 11 is a block diagram showing an example of the configuration of the signal processing circuit 120 according to the second embodiment. As shown in Fig. 11, the signal processing circuit 120 includes a wall filter processing function 121C, an autocorrelation calculation function 122C, and a Doppler information calculation function 123C in addition to the configuration shown in Fig. 2. Each processing function executed by the signal processing circuit 120 according to the second embodiment will be described later in the explanation of Fig. 12.

図12は、第2の実施形態に係る超音波診断装置1の処理手順の一例を示すフローチャートである。図12に示すステップS301~ステップS310の各処理は、例えば、ドプラ画像データの表示を開始する旨の開始指示を操作者から受け付けた場合に開始される。 Figure 12 is a flowchart showing an example of a processing procedure of the ultrasound diagnostic device 1 according to the second embodiment. Each process of steps S301 to S310 shown in Figure 12 is started, for example, when a start instruction to start displaying Doppler image data is received from the operator.

図12に示すように、超音波診断装置1は、超音波スキャンを実行する(ステップS101)。なお、ステップS301の処理は、図3のステップS101の処理と同様であるので説明を省略する。 As shown in FIG. 12, the ultrasound diagnostic device 1 executes an ultrasound scan (step S101). Note that the process of step S301 is similar to the process of step S101 in FIG. 3, and therefore a description thereof will be omitted.

続いて、信号処理回路120は、高カットオフ値のフィルタ処理Aと、低カットオフ値1のフィルタ処理Bと、低カットオフ値2のフィルタ処理Cとをそれぞれ適用する(ステップS302)。ステップS302の処理は、ウォールフィルタ処理機能121A,121B,121Cによって実行される。なお、図11に示すウォールフィルタ処理機能121A,121Bの説明は、図2に示したウォールフィルタ処理機能121A,121Bの説明と同様であるので説明を省略する。 Then, the signal processing circuit 120 applies filter processing A with a high cutoff value, filter processing B with a low cutoff value 1, and filter processing C with a low cutoff value 2 (step S302). The processing of step S302 is executed by the wall filter processing functions 121A, 121B, and 121C. Note that the explanation of the wall filter processing functions 121A and 121B shown in FIG. 11 is the same as the explanation of the wall filter processing functions 121A and 121B shown in FIG. 2, so the explanation will be omitted.

ウォールフィルタ処理機能121Cは、カットオフ値B(ウォールフィルタ処理機能121のカットオフ値)より小さいカットオフ値Cを有するフィルタ処理Cを実行する。ウォールフィルタ処理機能121Cは、ウォールフィルタ処理機能121A,121Bと並列でフィルタ処理を実行する。なお、ウォールフィルタ処理機能121Cの構成は、カットオフ値が異なる点を除き、ウォールフィルタ処理機能121A,121Bと同様であるので説明を省略する。また、ウォールフィルタ処理機能121Cは、第3フィルタ処理部の一例である。 The wall filter processing function 121C executes filter processing C having a cutoff value C smaller than the cutoff value B (the cutoff value of the wall filter processing function 121). The wall filter processing function 121C executes filter processing in parallel with the wall filter processing functions 121A and 121B. Note that the configuration of the wall filter processing function 121C is similar to that of the wall filter processing functions 121A and 121B except for the difference in the cutoff value, and therefore a description thereof will be omitted. Also, the wall filter processing function 121C is an example of a third filter processing unit.

そして、信号処理回路120は、フィルタ処理Aに基づくドプラ情報Aと、フィルタ処理Bに基づくドプラ情報B、フィルタ処理Cに基づくドプラ情報Cとを算出する(ステップS303)。ステップS303の処理は、自己相関演算機能122A,122B,122C、及びドプラ情報演算機能123A,123B,123Cによって実行される。
なお、図11に示す自己相関演算機能122A,122B、及びドプラ情報演算機能123A,123Bの説明は、図2に示した自己相関演算機能122A,122B、及びドプラ情報演算機能123A,123Bの説明と同様であるので説明を省略する。
Then, the signal processing circuit 120 calculates Doppler information A based on the filter processing A, Doppler information B based on the filter processing B, and Doppler information C based on the filter processing C (step S303). The process of step S303 is executed by the autocorrelation calculation functions 122A, 122B, and 122C and the Doppler information calculation functions 123A, 123B, and 123C.
The autocorrelation calculation functions 122A, 122B and the Doppler information calculation functions 123A, 123B shown in FIG. 11 are similar to the autocorrelation calculation functions 122A, 122B and the Doppler information calculation functions 123A, 123B shown in FIG. 2, and therefore will not be described again.

自己相関演算機能122Cは、フィルタ処理Cによって出力されたドプラ信号を用いて自己相関値を算出する。なお、自己相関演算機能122Cの構成は、フィルタ処理Cによって出力されたドプラ信号を用いる点を除き、自己相関演算機能122A,122Bと同様であるので説明を省略する。 The autocorrelation calculation function 122C calculates an autocorrelation value using the Doppler signal output by the filter process C. Note that the configuration of the autocorrelation calculation function 122C is similar to that of the autocorrelation calculation functions 122A and 122B, except for the use of the Doppler signal output by the filter process C, so a description thereof will be omitted.

ドプラ情報演算機能123Cは、自己相関演算機能122Cが算出した自己相関値からドプラ情報Cを算出する。なお、ドプラ情報演算機能123Cの構成は、自己相関演算機能122Cが算出した自己相関値を用いる点を除き、ドプラ情報演算機能123A,123Bと同様であるので説明を省略する。また、ドプラ情報Cは、第3ドプラ情報の一例である。 The Doppler information calculation function 123C calculates Doppler information C from the autocorrelation value calculated by the autocorrelation calculation function 122C. Note that the configuration of the Doppler information calculation function 123C is similar to that of the Doppler information calculation functions 123A and 123B, except for the use of the autocorrelation value calculated by the autocorrelation calculation function 122C, and therefore a description thereof will be omitted. Also, the Doppler information C is an example of third Doppler information.

そして、信号処理回路120は、ドプラ情報Bに基づいて、閾値を設定する(ステップS304)。そして、信号処理回路120は、ドプラ情報Aが閾値未満であるか否かを判定する(ステップS305)。そして、信号処理回路120は、ドプラ情報Aが閾値未満である場合には(ステップS305,Yes)、更に、ドプラ情報Bとドプラ情報Cとの差が所定値未満であるか否かを判定する(ステップS306)。そして、信号処理回路120は、ドプラ情報Bとドプラ情報Cとの差が所定値以上である場合には(ステップS305,Yes)、血流情報として略ゼロを採用する(ステップS306)。一方、信号処理回路120は、ドプラ情報Aが閾値以上である場合(ステップS305,No)、及び、ドプラ情報Bとドプラ情報Cとの差が所定値未満である場合には(ステップS306,Yes)、血流情報としてドプラ情報Aを採用する(ステップS308)。ステップS304~ステップS308の各処理は、ブランク処理機能124によって実行される。 Then, the signal processing circuit 120 sets a threshold value based on the Doppler information B (step S304). Then, the signal processing circuit 120 judges whether or not the Doppler information A is less than the threshold value (step S305). Then, when the Doppler information A is less than the threshold value (step S305, Yes), the signal processing circuit 120 further judges whether or not the difference between the Doppler information B and the Doppler information C is less than a predetermined value (step S306). Then, when the difference between the Doppler information B and the Doppler information C is equal to or greater than a predetermined value (step S305, Yes), the signal processing circuit 120 adopts approximately zero as the blood flow information (step S306). On the other hand, when the Doppler information A is equal to or greater than the threshold value (step S305, No) and when the difference between the Doppler information B and the Doppler information C is less than a predetermined value (step S306, Yes), the signal processing circuit 120 adopts the Doppler information A as the blood flow information (step S308). Each process from step S304 to step S308 is performed by the blank processing function 124.

例えば、ブランク処理機能124は、図10における平均速度「Vel-C3」と平均速度「Vel-C4」との差の大小を判定するために、「所定値」を利用する。この所定値は、受信信号に含まれるクラッタ信号の大きさを検出可能に設定される。 For example, the blanking function 124 uses a "predetermined value" to determine the magnitude of the difference between the average speed "Vel-C3" and the average speed "Vel-C4" in FIG. 10. This predetermined value is set so that the magnitude of the clutter signal contained in the received signal can be detected.

そして、ブランク処理機能124は、フレームごとに、各サンプル点のドプラ情報Aと閾値とを比較する。そして、ブランク処理機能124は、ドプラ情報Aが閾値以上である場合には、クラッタ信号に対して血流信号が大きく、ドプラ情報に含まれる成分の多くが血流信号であると推定できる。このため、ブランク処理機能124は、ドプラ情報Aが閾値以上であるサンプル点については、血流情報としてドプラ情報Aを採用する。 The blanking function 124 then compares the Doppler information A of each sample point with a threshold value for each frame. If the Doppler information A is equal to or greater than the threshold value, the blanking function 124 can estimate that the blood flow signal is large relative to the clutter signal, and that most of the components contained in the Doppler information are blood flow signals. For this reason, the blanking function 124 employs the Doppler information A as the blood flow information for sample points where the Doppler information A is equal to or greater than the threshold value.

一方、ブランク処理機能124は、ドプラ情報Aが閾値未満である場合には、更に、ドプラ情報Bとドプラ情報Cとの差と所定値とを比較する。そして、ブランク処理機能124は、ドプラ情報Bとドプラ情報Cとの差が所定値未満である場合には、図10に示したように、血流信号に対してクラッタ信号が小さく、血流信号が支配的であると推定できる。このため、ブランク処理機能124は、ドプラ情報Aが閾値未満であり、かつ、ドプラ情報Bとドプラ情報Cとの差が所定値未満であるサンプル点については、血流情報としてドプラ情報Aを採用する。 On the other hand, if Doppler information A is less than the threshold value, blanking function 124 further compares the difference between Doppler information B and Doppler information C with a predetermined value. Then, if the difference between Doppler information B and Doppler information C is less than the predetermined value, blanking function 124 can estimate that the clutter signal is small relative to the blood flow signal and the blood flow signal is dominant, as shown in FIG. 10. Therefore, for sample points where Doppler information A is less than the threshold value and the difference between Doppler information B and Doppler information C is less than the predetermined value, blanking function 124 adopts Doppler information A as the blood flow information.

また、ブランク処理機能124は、ドプラ情報Bとドプラ情報Cとの差が所定値以上である場合には、血流信号に対してクラッタ信号が大きく、ドプラ情報に含まれる成分の多くがクラッタ信号であると推定できる。このため、ブランク処理機能124は、ドプラ情報Aが閾値未満であり、かつ、ドプラ情報Bとドプラ情報Cとの差が所定値以上であるサンプル点については、血流情報として略ゼロを採用する。 In addition, when the difference between Doppler information B and Doppler information C is equal to or greater than a predetermined value, the blanking function 124 can estimate that the clutter signal is large relative to the blood flow signal, and that most of the components contained in the Doppler information are clutter signals. For this reason, the blanking function 124 adopts approximately zero as the blood flow information for sample points where Doppler information A is less than the threshold value and the difference between Doppler information B and Doppler information C is equal to or greater than a predetermined value.

このように、ブランク処理機能124は、ドプラ情報Aと、ドプラ情報Bと、ドプラ情報Cとに基づいて、血流信号を生成する。 In this way, the blanking processing function 124 generates a blood flow signal based on Doppler information A, Doppler information B, and Doppler information C.

そして、信号処理回路120は、平滑化処理を実行する(ステップS309)。そして、画像生成回路130は、ドプラ画像データを生成する(ステップS310)。なお、ステップS309及びステップS310の各処理は、図3に示したステップS108及びステップS109の各処理と同様であるので説明を省略する。 Then, the signal processing circuit 120 executes a smoothing process (step S309). Then, the image generating circuit 130 generates Doppler image data (step S310). Note that the processes in steps S309 and S310 are similar to those in steps S108 and S109 shown in FIG. 3, and therefore will not be described.

このように、第2の実施形態に係る超音波診断装置1は、3つ以上のウォールフィルタを用いることで、クラッタ信号成分を適切に除去した血流情報を生成することができる。これにより、第2の実施形態に係る超音波診断装置1は、血流信号に対してクラッタ信号が小さいスキャン部位に対しても適用することが可能となる。 In this way, the ultrasound diagnostic device 1 according to the second embodiment can generate blood flow information from which clutter signal components have been appropriately removed by using three or more wall filters. This makes it possible for the ultrasound diagnostic device 1 according to the second embodiment to be applied to scan sites where the clutter signal is small relative to the blood flow signal.

なお、第2の実施形態にて説明した内容はあくまで一例であり、実施形態はこれに限定されるものではない。例えば、第1の実施形態にて説明した変形例や他のバリエーションは、第2の実施形態にも同様に適用可能である。 Note that the content described in the second embodiment is merely an example, and the embodiment is not limited to this. For example, the modified examples and other variations described in the first embodiment can be similarly applied to the second embodiment.

(その他の実施形態)
上述した実施形態以外にも、種々の異なる形態にて実施されてもよい。
Other Embodiments
In addition to the above-described embodiment, the present invention may be embodied in various different forms.

図示した各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部又は一部を、各種の負荷や使用状況等に応じて、任意の単位で機能的又は物理的に分散・統合して構成することができる。更に、各装置にて行なわれる各処理機能は、その全部又は任意の一部が、CPU及び当該CPUにて解析実行されるプログラムにて実現され、或いは、ワイヤードロジックによるハードウェアとして実現され得る。 The components of each device shown in the figure are conceptual and functional, and do not necessarily have to be physically configured as shown. In other words, the specific form of distribution and integration of each device is not limited to that shown, and all or part of it can be functionally or physically distributed and integrated in any unit depending on various loads, usage conditions, etc. Furthermore, all or any part of each processing function performed by each device can be realized by a CPU and a program analyzed and executed by the CPU, or can be realized as hardware using wired logic.

また、上記の実施形態において説明した各処理のうち、自動的に行なわれるものとして説明した処理の全部又は一部を手動的に行なうこともでき、或いは、手動的に行なわれるものとして説明した処理の全部又は一部を公知の方法で自動的に行なうこともできる。この他、上記文書中や図面中で示した処理手順、制御手順、具体的名称、各種のデータやパラメータを含む情報については、特記する場合を除いて任意に変更することができる。 Furthermore, among the processes described in the above embodiments, all or part of the processes described as being performed automatically can be performed manually, or all or part of the processes described as being performed manually can be performed automatically using known methods. In addition, the information including the processing procedures, control procedures, specific names, various data, and parameters shown in the above documents and drawings can be changed as desired unless otherwise specified.

また、上記の実施形態において説明した撮像方法は、予め用意された撮像プログラムをパーソナルコンピュータやワークステーション等のコンピュータで実行することによって実現することができる。この撮像プログラムは、インターネット等のネットワークを介して配布することができる。また、この超音波イメージング方法は、ハードディスク、フレキシブルディスク(FD)、CD-ROM、MO、DVD等のコンピュータで読み取り可能な非一過性の記録媒体に記録され、コンピュータによって記録媒体から読み出されることによって実行することもできる。 The imaging method described in the above embodiment can be realized by executing a prepared imaging program on a computer such as a personal computer or a workstation. This imaging program can be distributed via a network such as the Internet. This ultrasound imaging method can also be recorded on a non-transitory recording medium that can be read by a computer, such as a hard disk, a flexible disk (FD), a CD-ROM, an MO, or a DVD, and executed by being read from the recording medium by a computer.

以上説明した少なくとも1つの実施形態によれば、クラッタ信号成分を適切に除去することができる。 According to at least one of the embodiments described above, clutter signal components can be appropriately removed.

いくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更、実施形態同士の組み合わせを行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。 Although several embodiments have been described, these embodiments are presented as examples and are not intended to limit the scope of the invention. These embodiments can be implemented in various other forms, and various omissions, substitutions, modifications, and combinations of embodiments can be made without departing from the spirit of the invention. These embodiments and their modifications are within the scope of the invention and its equivalents as set forth in the claims, as well as the scope and spirit of the invention.

1 超音波診断装置
120 信号処理回路
121 ウォールフィルタ処理機能
122 自己相関演算機能
123 ドプラ情報演算機能
124 ブランク処理機能
125 パーシスタンス機能
REFERENCE SIGNS LIST 1 Ultrasound diagnostic device 120 Signal processing circuit 121 Wall filter processing function 122 Autocorrelation calculation function 123 Doppler information calculation function 124 Blank processing function 125 Persistence function

Claims (6)

互いに異なるカットオフ値を有する複数のフィルタ処理を超音波の受信信号に対してそれぞれ並列に適用することによって抽出された複数のドプラ情報に基づいて、サンプル点ごとにクラッタ信号成分が除去された血流情報を生成する信号処理部
を備え
前記信号処理部は、
サンプル点ごとに時系列的に受信された前記受信信号に対して、第1カットオフ値を有する第1フィルタ処理を実行することで、各サンプル点の第1ドプラ情報を出力する第1フィルタ処理部と、
前記受信信号に対して、前記第1カットオフ値より小さい第2カットオフ値を有する第2フィルタ処理を実行することで、各サンプル点の第2ドプラ情報を出力する第2フィルタ処理部と、
各サンプル点の第2ドプラ情報に基づいてブランク処理に用いる閾値を前記サンプル点ごとに設定し、設定した前記閾値と前記第1ドプラ情報とを前記サンプル点ごとに比較することで、前記血流情報を生成する血流情報生成部と、
を有する、超音波診断装置。
a signal processing unit that generates blood flow information from which clutter signal components have been removed for each sample point based on a plurality of pieces of Doppler information extracted by applying a plurality of filter processes having mutually different cutoff values in parallel to the received ultrasonic signal ,
The signal processing unit includes:
a first filter processing unit that performs a first filter process having a first cutoff value on the received signal received in time series for each sample point, and outputs first Doppler information for each sample point;
a second filter processing unit that performs a second filter process having a second cutoff value smaller than the first cutoff value on the received signal, and outputs second Doppler information of each sample point;
a blood flow information generating unit that sets a threshold value used for blanking for each sample point based on second Doppler information of the sample point, and generates the blood flow information by comparing the set threshold value with the first Doppler information for each sample point;
An ultrasound diagnostic device comprising :
前記血流情報生成部は、
前記第1ドプラ情報から前記第2ドプラ情報を部分的に除く前記ブランク処理を実行することで、前記血流情報を生成する、
請求項に記載の超音波診断装置。
The blood flow information generating unit
generating the blood flow information by performing the blanking process to partially remove the second Doppler information from the first Doppler information;
The ultrasonic diagnostic apparatus according to claim 1 .
前記血流情報生成部は、
前記第1ドプラ情報が前記閾値以上であるサンプル点については、前記血流情報として前記第1ドプラ情報を採用し、
前記第1ドプラ情報が前記閾値未満であるサンプル点については、前記血流情報として略ゼロを採用する、
請求項に記載の超音波診断装置。
The blood flow information generating unit
For sample points where the first Doppler information is equal to or greater than the threshold, the first Doppler information is adopted as the blood flow information;
For a sample point where the first Doppler information is less than the threshold, approximately zero is adopted as the blood flow information.
The ultrasonic diagnostic apparatus according to claim 1 .
前記血流情報生成部は、前記第2ドプラ情報に対して予め設定された固定値を加算した値を、前記閾値として設定する、
請求項1~3のいずれか一つに記載の超音波診断装置。
the blood flow information generating unit sets a value obtained by adding a preset fixed value to the second Doppler information as the threshold value;
4. The ultrasonic diagnostic apparatus according to claim 1 .
前記第1フィルタ処理部は、前記第1フィルタ処理として、ハイパスフィルタを適用し、
前記第2フィルタ処理部は、前記第2フィルタ処理として、ハイパスフィルタ又はバンドパスフィルタを適用する、
請求項のいずれか一つに記載の超音波診断装置。
The first filtering unit applies a high-pass filter as the first filtering process,
The second filtering unit applies a high-pass filter or a band-pass filter as the second filtering process.
The ultrasonic diagnostic apparatus according to any one of claims 1 to 4 .
前記信号処理部は、
サンプル点ごとに時系列的に受信された超音波の受信信号に対して、第1カットオフ値を有する第1フィルタ処理を実行する第1フィルタ処理部と、
前記受信信号に対して、前記第1カットオフ値より小さい第2カットオフ値を有する第2フィルタ処理を実行する第2フィルタ処理部と、
前記受信信号に対して、前記第2カットオフ値より小さい第3カットオフ値を有する第3フィルタ処理を実行する第3フィルタ処理部と、
前記第1フィルタ処理によって抽出された各サンプル点の第1ドプラ情報と、前記第2フィルタ処理によって抽出された各サンプル点の第2ドプラ情報と、前記第3フィルタ処理によって抽出された各サンプル点の第3ドプラ情報とに基づいて、前記クラッタ信号成分が除去された血流情報を生成する血流情報生成部と
を有する、
請求項1に記載の超音波診断装置。
The signal processing unit includes:
a first filter processing unit that executes a first filter process having a first cutoff value on a reception signal of an ultrasonic wave received in time series for each sample point;
a second filter processing unit that performs a second filter process on the received signal, the second filter process having a second cutoff value smaller than the first cutoff value;
a third filter processing unit that performs a third filter process on the received signal, the third filter process having a third cutoff value smaller than the second cutoff value;
a blood flow information generating unit that generates blood flow information from which the clutter signal components have been removed, based on first Doppler information of each sample point extracted by the first filtering process, second Doppler information of each sample point extracted by the second filtering process, and third Doppler information of each sample point extracted by the third filtering process.
The ultrasonic diagnostic apparatus according to claim 1 .
JP2021063689A 2021-04-02 2021-04-02 Ultrasound diagnostic equipment Active JP7613989B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2021063689A JP7613989B2 (en) 2021-04-02 2021-04-02 Ultrasound diagnostic equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2021063689A JP7613989B2 (en) 2021-04-02 2021-04-02 Ultrasound diagnostic equipment

Publications (2)

Publication Number Publication Date
JP2022158648A JP2022158648A (en) 2022-10-17
JP7613989B2 true JP7613989B2 (en) 2025-01-15

Family

ID=83638528

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021063689A Active JP7613989B2 (en) 2021-04-02 2021-04-02 Ultrasound diagnostic equipment

Country Status (1)

Country Link
JP (1) JP7613989B2 (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005177338A (en) 2003-12-22 2005-07-07 Aloka Co Ltd Ultrasonic diagnostic instrument
JP2008178466A (en) 2007-01-23 2008-08-07 Toshiba Corp Angiography equipment
WO2014142174A1 (en) 2013-03-12 2014-09-18 株式会社 東芝 Ultrasonic diagnostic device and ultrasonic image processing method

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3029315B2 (en) * 1990-03-30 2000-04-04 株式会社東芝 Ultrasound blood flow imaging device
US5111695A (en) * 1990-07-11 1992-05-12 General Electric Company Dynamic phase focus for coherent imaging beam formation
JPH06327673A (en) * 1993-05-27 1994-11-29 Toshiba Corp Ultrasonic diagnostic equipment

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005177338A (en) 2003-12-22 2005-07-07 Aloka Co Ltd Ultrasonic diagnostic instrument
JP2008178466A (en) 2007-01-23 2008-08-07 Toshiba Corp Angiography equipment
WO2014142174A1 (en) 2013-03-12 2014-09-18 株式会社 東芝 Ultrasonic diagnostic device and ultrasonic image processing method

Also Published As

Publication number Publication date
JP2022158648A (en) 2022-10-17

Similar Documents

Publication Publication Date Title
JP6253999B2 (en) Ultrasonic diagnostic apparatus, image processing apparatus, and image processing method
US10603014B2 (en) Ultrasonic diagnostic apparatus, image processing apparatus, and image processing method
CN104883981B (en) Diagnostic ultrasound equipment, image processing apparatus and image processing method
US20140039317A1 (en) Ultrasound diagnosis apparatus and controlling method
JP6580925B2 (en) Ultrasonic diagnostic apparatus and image processing apparatus
US10835200B2 (en) Ultrasound diagnosis apparatus
JP6697609B2 (en) Ultrasonic diagnostic device, image processing device, and image processing method
JP6460707B2 (en) Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, and ultrasonic image processing program
JP7282492B2 (en) Ultrasound diagnostic device, medical image processing device and medical image processing program
US12551192B2 (en) Ultrasonic diagnostic apparatus and medical image processing apparatus
JP7034686B2 (en) Ultrasound diagnostic equipment, medical image processing equipment and their programs
US20170224309A1 (en) Ultrasonic diagnostic device
JP6651405B2 (en) Ultrasound diagnostic apparatus and program
JP7305438B2 (en) Analysis device and program
JP7613989B2 (en) Ultrasound diagnostic equipment
US10709421B2 (en) Ultrasound diagnostic apparatus
JP7693408B2 (en) Ultrasound diagnostic equipment
JP7560240B2 (en) Ultrasound diagnostic device and ultrasound diagnostic program
JP6793502B2 (en) Ultrasonic diagnostic equipment
JP6843591B2 (en) Ultrasonic diagnostic equipment
JP2024121628A (en) Ultrasound diagnostic device and ultrasound diagnostic method
JP2024034087A (en) Ultrasound diagnostic device and blood flow image data generation method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20240301

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20240830

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20240904

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20241031

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20241127

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20241226

R150 Certificate of patent or registration of utility model

Ref document number: 7613989

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150