JP7469877B2 - Ultrasound diagnostic device, medical image processing device, and medical image processing program - Google Patents

Ultrasound diagnostic device, medical image processing device, and medical image processing program Download PDF

Info

Publication number
JP7469877B2
JP7469877B2 JP2019238036A JP2019238036A JP7469877B2 JP 7469877 B2 JP7469877 B2 JP 7469877B2 JP 2019238036 A JP2019238036 A JP 2019238036A JP 2019238036 A JP2019238036 A JP 2019238036A JP 7469877 B2 JP7469877 B2 JP 7469877B2
Authority
JP
Japan
Prior art keywords
blood flow
region
threshold
threshold value
power
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
JP2019238036A
Other languages
Japanese (ja)
Other versions
JP2021104301A (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 JP2019238036A priority Critical patent/JP7469877B2/en
Publication of JP2021104301A publication Critical patent/JP2021104301A/en
Application granted granted Critical
Publication of JP7469877B2 publication Critical patent/JP7469877B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本明細書および図面に開示の実施形態は、超音波診断装置、医用画像処理装置、および医用画像処理プログラムに関する。 The embodiments disclosed in this specification and the drawings relate to an ultrasound diagnostic device, a medical image processing device, and a medical image processing program.

超音波診装置には、超音波ドプラ画像をカラーで表示させるカラードプラモード(CFM(Color Flow Mapping)モード)での撮影が可能なものがある。CFMモードでは、被検体の形態が描出されたBモード像上に、血流信号(血液からの反射信号)にもとづいて求められた超音波ドプラ画像がカラーで重畳表示される。カラー表示される超音波ドプラ画像は、被検体内の動体において反射することによって周波数がドプラシフトしたドプラ信号にもとづいて求められた血流動態情報を示す画像である。ドプラ信号は、同一場所を所定時間間隔で走査して得られる血流信号の単位時間内の位相シフト量(ドプラシフト量)から求められる。超音波ドプラ画像によれば、血流を可視化することができる。 Some ultrasound diagnostic devices are capable of capturing images in color Doppler mode (CFM (Color Flow Mapping) mode), which displays ultrasonic Doppler images in color. In CFM mode, an ultrasonic Doppler image calculated based on blood flow signals (reflected signals from blood) is superimposed in color on a B-mode image depicting the morphology of the subject. The color ultrasonic Doppler image shows blood flow dynamics information calculated based on Doppler signals whose frequency has been Doppler shifted by reflection from moving objects within the subject. The Doppler signal is calculated from the phase shift amount (Doppler shift amount) per unit time of the blood flow signal obtained by scanning the same location at a specified time interval. The ultrasonic Doppler image makes it possible to visualize blood flow.

しかし、超音波ドプラ画像をBモード画像に重畳する際に、超音波ドプラ画像に描出された血流が、Bモード画像に描出された血管径を大幅にはみ出して表示されるブルーミングあるいはオーバーペイントと呼ばれる現象が生じることがある。ブルーミングは、観察部位によっては心臓の拍動が原因となって生じる。また、下肢では、手で脚をもむように圧迫するミルキング操作が原因となってブルーミングが生じることもある。ブルーミングが生じてしまうと、本来血流が存在しない領域に血流が描出されてしまうため、超音波ドプラ画像における血管の位置の把握および血流動態の把握が著しく困難となり、診断の妨げとなってしまう。 However, when an ultrasound Doppler image is superimposed on a B-mode image, a phenomenon called blooming or overpainting may occur, in which the blood flow depicted in the ultrasound Doppler image extends beyond the diameter of the blood vessel depicted in the B-mode image. Blooming can occur due to the beating of the heart depending on the area being observed. In addition, in the lower limbs, blooming can also occur due to the milking operation, in which the legs are massaged with the hands. When blooming occurs, blood flow is depicted in areas where there is no blood flow, making it extremely difficult to grasp the position of blood vessels and blood flow dynamics in the ultrasound Doppler image, which impedes diagnosis.

特開2014-161554号公報JP 2014-161554 A

本明細書および図面に開示の実施形態が解決しようとする課題の一つは、カラードプラモードのブルーミングを低減することである。ただし、本明細書および図面に開示の実施形態により解決しようとする課題は上記課題に限られない。後述する実施形態に示す各構成による各効果に対応する課題を他の課題として位置づけることもできる。 One of the problems that the embodiments disclosed in this specification and the drawings attempt to solve is to reduce blooming in color Doppler mode. However, the problems that the embodiments disclosed in this specification and the drawings attempt to solve are not limited to the above problem. Problems that correspond 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 determination unit and an image generation unit. The determination unit determines at least one of a power threshold and a variance threshold for displaying ultrasound Doppler image data relating to the subject's blood flow velocity on the display unit based on at least one of the power and variance of the subject's blood flow signal. The image generation unit generates ultrasound Doppler image data relating to the subject's blood flow velocity based on the threshold and displays the data on the display unit.

一実施形態に係る超音波診断装置の一構成例を示すブロック図。FIG. 1 is a block diagram showing an example of the configuration of an ultrasound diagnostic apparatus according to an embodiment. 本実施形態に係るドプラ処理回路の一構成例を示す説明図。FIG. 2 is an explanatory diagram showing an example of the configuration of a Doppler processing circuit according to the present embodiment. 従来のカラードプラ画像の一例を示す説明図。FIG. 1 is an explanatory diagram showing an example of a conventional color Doppler image. (a)はパワー値画像の一例を示す説明図、(b)はパワー値のヒストグラムの一例を示す説明図。FIG. 4A is an explanatory diagram showing an example of a power value image, and FIG. 4B is an explanatory diagram showing an example of a histogram of power values. 本実施形態に係る画像生成回路により生成される超音波ドプラ画像の一例を示す説明図。FIG. 2 is an explanatory diagram showing an example of an ultrasonic Doppler image generated by the image generating circuit according to the embodiment. 閾値に関するユーザの指示を受け付ける様子の一例を示す説明。13 is an explanation showing an example of how a user's instruction regarding a threshold is accepted. ヒストグラムと閾値を示す情報とを同時に表示する場合の例を示す説明図。FIG. 13 is an explanatory diagram showing an example in which a histogram and information indicating a threshold value are displayed simultaneously. (a)は関心領域を深さ方向に3つに分割する場合の分割領域を説明するための図、(b)はスキャン方向に3つに分割する場合の分割領域を説明するための図。FIG. 2A is a diagram for explaining divided regions when a region of interest is divided into three in the depth direction, and FIG. 2B is a diagram for explaining divided regions when a region of interest is divided into three in the scanning direction. (a)は深さ方向に並ぶ3つの分割領域ごとにヒストグラムにもとづいて決定された閾値の一例を示す説明図、(b)は深さと閾値との関係を事前に規定した複数のカーブの一例を示す説明図。FIG. 1A is an explanatory diagram showing an example of thresholds determined based on a histogram for each of three divided regions aligned in the depth direction, and FIG. 1B is an explanatory diagram showing an example of multiple curves in which the relationship between depth and threshold is predefined. 減衰値の測定値に応じて閾値の補正値を決定する場合の例を説明するための図。11A and 11B are diagrams for explaining an example of a case in which a correction value for a threshold is determined in accordance with a measured value of an attenuation value. 関心領域の位置に応じた仮閾値の平均化処理を説明するための図。11A and 11B are diagrams for explaining averaging processing of provisional thresholds according to the position of a region of interest. パワーと分散の両者にもとづいて閾値を決定するためのテーブルの一例を示す説明図。FIG. 13 is an explanatory diagram showing an example of a table for determining a threshold value based on both power and variance.

以下、図面を参照しながら、超音波診断装置、医用画像処理装置、および医用画像処理プログラムの実施形態について詳細に説明する。 Below, embodiments of an ultrasound diagnostic device, a medical image processing device, and a medical image processing program will be described in detail with reference to the drawings.

図1は、一実施形態に係る超音波診断装置10の一構成例を示すブロック図である。超音波診断装置10は、超音波プローブ20、入力インターフェース21、およびディスプレイ22と接続されて用いることができる。なお、超音波診断装置10は、超音波プローブ20、入力インターフェース21、およびディスプレイ22の少なくとも1つを備えてもよい。 FIG. 1 is a block diagram showing an example of the configuration of an ultrasound diagnostic device 10 according to an embodiment. The ultrasound diagnostic device 10 can be used by being connected to an ultrasound probe 20, an input interface 21, and a display 22. Note that the ultrasound diagnostic device 10 may include at least one of the ultrasound probe 20, the input interface 21, and the display 22.

超音波診断装置10は、図1に示すように、送受信回路11、Bモード処理回路12、ドプラ処理回路13、画像生成回路14、画像メモリ15、記憶回路16、ネットワーク接続回路17、および制御回路18を有する。 As shown in FIG. 1, the ultrasound diagnostic device 10 has a transmission/reception circuit 11, a B-mode processing circuit 12, a Doppler processing circuit 13, an image generation circuit 14, an image memory 15, a storage circuit 16, a network connection circuit 17, and a control circuit 18.

送受信回路11は、送信回路および受信回路を有する。送受信回路11は、制御回路18に制御されて、超音波の送受信における送信指向性と受信指向性とを制御する。なお、図1には送受信回路11が超音波診断装置10に設けられる場合の例について示したが、送受信回路11は超音波プローブ20に設けられてもよいし、超音波診断装置10と超音波プローブ20の両方に設けられてもよい。 The transmission/reception circuit 11 has a transmission circuit and a reception circuit. The transmission/reception circuit 11 is controlled by the control circuit 18 to control the transmission directivity and reception directivity in transmitting and receiving ultrasound. Note that while FIG. 1 shows an example in which the transmission/reception circuit 11 is provided in the ultrasound diagnostic device 10, the transmission/reception circuit 11 may be provided in the ultrasound probe 20, or may be provided in both the ultrasound diagnostic device 10 and the ultrasound probe 20.

送信回路は、パルス発生器、送信遅延回路およびパルサ回路などを有し、超音波振動子に駆動信号を供給する。パルス発生器は、所定のレート周波数で、送信超音波を形成するためのレートパルスを繰り返し発生する。送信遅延回路は、超音波振動子から発生される超音波をビーム状に集束して送信指向性を決定するために必要な圧電振動子ごとの遅延時間を、パルス発生器が発生する各レートパルスに対し与える。また、パルサ回路は、レートパルスにもとづくタイミングで、超音波振動子に駆動パルスを印加する。送信遅延回路は、各レートパルスに対し与える遅延時間を変化させることで、圧電振動子面から送信される超音波ビームの送信方向を任意に調整する。 The transmission circuit has a pulse generator, a transmission delay circuit, a pulsar circuit, etc., and supplies a drive signal to the ultrasonic transducer. The pulse generator repeatedly generates rate pulses to form a transmission ultrasonic wave at a predetermined rate frequency. The transmission delay circuit provides each rate pulse generated by the pulse generator with a delay time for each piezoelectric transducer required to focus the ultrasonic waves generated from the ultrasonic transducer into a beam and determine the transmission directivity. In addition, the pulsar circuit applies a drive pulse to the ultrasonic transducer at a timing based on the rate pulse. The transmission delay circuit changes the delay time provided to each rate pulse to arbitrarily adjust the transmission direction of the ultrasonic beam transmitted from the piezoelectric transducer surface.

受信回路は、アンプ回路、A/D変換器、加算器などを有し、超音波振動子が受信したエコー信号を受け、このエコー信号に対して各種処理を行なってエコーデータを生成する。アンプ回路は、エコー信号をチャンネルごとに増幅してゲイン補正処理を行なう。A/D変換器は、ゲイン補正されたエコー信号をA/D変換し、デジタルデータに受信指向性を決定するのに必要な遅延時間を与える。加算器は、A/D変換器によって処理されたエコー信号の加算処理を行なってエコーデータを生成する。加算器の加算処理により、エコー信号の受信指向性に応じた方向からの反射成分が強調される。 The receiving circuit has an amplifier circuit, an A/D converter, an adder, etc., and receives the echo signal received by the ultrasonic transducer and performs various processes on this echo signal to generate echo data. The amplifier circuit amplifies the echo signal for each channel and performs gain correction processing. The A/D converter A/D converts the gain-corrected echo signal and gives the digital data the delay time required to determine the reception directivity. The adder performs addition processing of the echo signal processed by the A/D converter to generate echo data. The addition processing of the adder emphasizes the reflected components from the direction corresponding to the reception directivity of the echo signal.

Bモード処理回路12は、受信回路からエコーデータを受信し、対数増幅、包絡線検波処理などを行なって、信号強度が輝度の明るさで表現されるデータ(Bモードデータ)を生成する。 The B-mode processing circuit 12 receives the echo data from the receiving circuit and performs logarithmic amplification, envelope detection processing, etc. to generate data (B-mode data) in which the signal strength is expressed as luminance brightness.

ドプラ処理回路13は、受信回路から受信したエコーデータの血流信号から速度情報を周波数解析し、ドプラ効果による血流や組織、造影剤エコー成分を抽出し、平均速度、分散、パワーなどの血流動態情報を多点について抽出したデータ(ドプラデータ)を生成する。また、ドプラ処理回路13は、血流信号のパワーおよび速度成分の分散の少なくとも一方にもとづいて、被検体の血流の速度に関する超音波ドプラ画像のデータをディスプレイ22に表示させる際のパワーの閾値および分散の閾値の少なくとも一方を決定する。ドプラ処理回路13の構成の詳細については図2を用いて後述する。 The Doppler processing circuit 13 performs frequency analysis of velocity information from the blood flow signal of the echo data received from the receiving circuit, extracts blood flow, tissue, and contrast agent echo components due to the Doppler effect, and generates data (Doppler data) that extracts blood flow dynamics information such as average velocity, variance, and power for multiple points. In addition, the Doppler processing circuit 13 determines at least one of the power threshold and variance threshold when displaying ultrasonic Doppler image data related to the subject's blood flow velocity on the display 22 based on at least one of the power of the blood flow signal and the variance of the velocity component. The configuration of the Doppler processing circuit 13 will be described in detail later with reference to FIG. 2.

画像生成回路14は、超音波プローブ20が受信したエコー信号にもとづいて超音波画像データを生成し、ディスプレイ22に表示させる。たとえば、画像生成回路14は、Bモード処理回路12から受けたBモードデータから反射波の強度を輝度にて表した2次元Bモード画像データを生成する。また、画像生成回路14は、ドプラ処理回路13から受けたドプラデータおよび閾値にもとづいて、血流動態情報を示す画像としての被検体の血流の速度に関する超音波ドプラ画像のデータを生成する。画像生成回路14は、画像生成部の一例である。 The image generation circuit 14 generates ultrasound image data based on the echo signals received by the ultrasound probe 20 and displays the data on the display 22. For example, the image generation circuit 14 generates two-dimensional B-mode image data that represents the intensity of the reflected waves as brightness from the B-mode data received from the B-mode processing circuit 12. The image generation circuit 14 also generates ultrasound Doppler image data related to the blood flow velocity of the subject as an image showing hemodynamic information based on the Doppler data and threshold value received from the Doppler processing circuit 13. The image generation circuit 14 is an example of an image generation unit.

画像メモリ15は、画像生成回路14が生成したBモード画像データや超音波ドプラ画像データなどの画像データを記憶するメモリである。なお、画像メモリ15は、Bモード処理回路12やドプラ処理回路13が生成したデータを記憶してもよい。 The image memory 15 is a memory that stores image data such as B-mode image data and ultrasonic Doppler image data generated by the image generation circuit 14. The image memory 15 may also store data generated by the B-mode processing circuit 12 and the Doppler processing circuit 13.

記憶回路16は、磁気的もしくは光学的記録媒体または半導体メモリなどの、プロセッサにより読み取り可能な記録媒体を含んだ構成を有する。記憶回路16の記憶媒体内のプログラムおよびデータの一部または全部は電子ネットワークを介した通信によりダウンロードされてもよいし、光ディスクなどの可搬型記憶媒体を介して記憶回路16に与えられてもよい。なお、記憶回路16に記憶される情報の一部または全部は、外部の記憶回路や超音波プローブ20の図示しない記憶回路などの記憶媒体の少なくとも1つに分散されて記憶され、あるいは複製されて記憶されてもよい。 The memory circuit 16 has a configuration including a processor-readable recording medium, such as a magnetic or optical recording medium or a semiconductor memory. Some or all of the programs and data in the storage medium of the memory circuit 16 may be downloaded by communication via an electronic network, or may be provided to the memory circuit 16 via a portable storage medium such as an optical disk. Some or all of the information stored in the memory circuit 16 may be distributed and stored in at least one storage medium, such as an external storage circuit or a storage circuit (not shown) of the ultrasound probe 20, or may be copied and stored.

ネットワーク接続回路17は、ネットワーク100の形態に応じた種々の情報通信用プロトコルを実装する。ネットワーク接続回路17は、この各種プロトコルに従って超音波診断装置10と他の電気機器とを接続する。この接続には、電子ネットワークを介した電気的な接続などを適用することができる。ここで電子ネットワークとは、電気通信技術を利用した情報通信網全般を意味し、病院基幹LAN(Local Area Network)などの無線/有線LANやインターネット網のほか、電話通信回線網、光ファイバ通信ネットワーク、ケーブル通信ネットワークおよび衛星通信ネットワークなどを含む。 The network connection circuit 17 implements various information communication protocols according to the form of the network 100. The network connection circuit 17 connects the ultrasound diagnostic device 10 to other electrical devices according to these various protocols. This connection can be an electrical connection via an electronic network. Here, the electronic network refers to any information communication network that uses electrical communication technology, and includes wireless/wired LANs such as hospital backbone LANs (Local Area Networks) and Internet networks, as well as telephone communication line networks, optical fiber communication networks, cable communication networks, and satellite communication networks.

制御回路18は、超音波診断装置10を統括制御する機能を実現する。 The control circuit 18 provides the function of overall control of the ultrasound diagnostic device 10.

超音波プローブ20は、超音波診断装置10とケーブルを介して着脱自在に接続される。なお、超音波プローブ20は超音波診断装置10と無線接続されてもよい。 The ultrasound probe 20 is detachably connected to the ultrasound diagnostic device 10 via a cable. The ultrasound probe 20 may also be wirelessly connected to the ultrasound diagnostic device 10.

超音波プローブ20としては、スキャン方向(アジマス方向)に複数の超音波振動子が配列されるとともにレンズ方向(エレベーション方向)にも複数の素子が配列された2次元アレイプローブを用いることができる。この種の2次元アレイプローブとしては、たとえば1.5Dアレイプローブ、1.75Dアレイプローブや、2Dアレイプローブなどを用いることができる。 The ultrasonic probe 20 may be a two-dimensional array probe in which multiple ultrasonic transducers are arranged in the scan direction (azimuth direction) and multiple elements are also arranged in the lens direction (elevation direction). Examples of this type of two-dimensional array probe that may be used include a 1.5D array probe, a 1.75D array probe, and a 2D array probe.

なお、超音波プローブ20は、ボリュームデータを取得可能に構成されてもよい。この場合、2次元アレイプローブである超音波プローブ20により被検体を3次元でスキャンしてもよいし、複数の圧電振動子が一列で配置された1次元超音波プローブである超音波プローブ20により被検体を2次元でスキャンするまたはこれら複数の超音波振動子を回転させることで被検体を3次元でスキャンしてもよいし、1次元超音波プローブの複数の圧電振動子を機械的に揺動してもよい。 The ultrasound probe 20 may be configured to acquire volume data. In this case, the subject may be scanned three-dimensionally using the ultrasound probe 20, which is a two-dimensional array probe, or the subject may be scanned two-dimensionally using the ultrasound probe 20, which is a one-dimensional ultrasound probe in which multiple piezoelectric transducers are arranged in a row, or the subject may be scanned three-dimensionally by rotating these multiple ultrasound transducers, or the multiple piezoelectric transducers of the one-dimensional ultrasound probe may be mechanically oscillated.

超音波プローブ20がボリュームデータを取得可能な場合、ユーザは、複数の2次元超音波画像のいずれかをリアルタイムな動画としてまたは静止画として表示させる2次元表示モード(2Dモード)と、リアルタイムに取得されている3次元超音波画像を動画として表示させる4次元表示モード(4Dモード)とのいずれかのモードを選択可能である。 When the ultrasound probe 20 is capable of acquiring volume data, the user can select either a two-dimensional display mode (2D mode) in which any of a number of two-dimensional ultrasound images is displayed as a real-time video or as a still image, or a four-dimensional display mode (4D mode) in which a three-dimensional ultrasound image acquired in real time is displayed as a video.

入力インターフェース21は、たとえばトラックボール、スイッチ、ボタン、マウス、キーボード、操作面へ触れることで入力操作を行なうタッチパッド、光学センサを用いた非接触入力回路、および音声入力回路等などの一般的な入力装置により実現され、ユーザの操作に対応した操作入力信号を制御回路18に出力する。入力インターフェース21は、操作パネルとして構成されてもよい。この場合、操作パネルは、タッチコマンドスクリーンとして機能し、たとえばディスプレイと、このディスプレイの近傍に設けられたタッチ入力回路と、ハードキーとを有する。 The input interface 21 is realized by general input devices such as a trackball, a switch, a button, a mouse, a keyboard, a touchpad that performs input operations by touching the operation surface, a non-contact input circuit using an optical sensor, and a voice input circuit, and outputs an operation input signal corresponding to a user's operation to the control circuit 18. The input interface 21 may be configured as an operation panel. In this case, the operation panel functions as a touch command screen and has, for example, a display, a touch input circuit provided near the display, and hard keys.

ディスプレイ22は、たとえば液晶ディスプレイやOLED(Organic Light Emitting Diode)ディスプレイなどの一般的な表示出力装置により構成され、制御回路18の制御に従って各種情報を表示する。 The display 22 is composed of a general display output device such as a liquid crystal display or an OLED (Organic Light Emitting Diode) display, and displays various information according to the control of the control circuit 18.

また、超音波診断装置10がタブレット型やスマートフォン型の超音波診断装置10である場合は、入力インターフェース21とディスプレイ22とが一体としてタッチパネルを構成してもよい。 In addition, if the ultrasound diagnostic device 10 is a tablet or smartphone type ultrasound diagnostic device 10, the input interface 21 and the display 22 may be integrated to form a touch panel.

また、超音波診断装置10は、医用画像処理装置30および画像サーバ40とネットワーク100を介して互いにデータ送受信可能に接続されてもよい。 The ultrasound diagnostic device 10 may also be connected to the medical image processing device 30 and the image server 40 via the network 100 so that data can be transmitted and received between them.

図2は、本実施形態に係るドプラ処理回路13の一構成例を示す説明図である。ドプラ処理回路13は、記憶回路16に記憶された医用画像処理プログラムを読み出して実行することにより、画像生成回路14と協働してカラードプラモードのブルーミングを低減するための処理を実行するプロセッサである。 Figure 2 is an explanatory diagram showing an example of the configuration of the Doppler processing circuitry 13 according to this embodiment. The Doppler processing circuitry 13 is a processor that reads out and executes a medical image processing program stored in the memory circuitry 16, thereby working with the image generation circuitry 14 to execute processing for reducing blooming in color Doppler mode.

図2に示すように、ドプラ処理回路13のプロセッサは、ウォールフィルタ機能31、自己相関機能32、ヒストグラム生成機能33、閾値決定機能34、補正機能35、およびブランク機能36を実現する。これらの各機能は、それぞれプログラムの形態で記憶回路16に記憶されている。なお、これらの各機能の一部または全部は、プロセッサを用いずハードウエアロジックによって実現されてもよい。 As shown in FIG. 2, the processor of the Doppler processing circuitry 13 realizes a wall filter function 31, an autocorrelation function 32, a histogram generation function 33, a threshold determination function 34, a correction function 35, and a blank function 36. Each of these functions is stored in the memory circuitry 16 in the form of a program. Note that some or all of these functions may be realized by hardware logic without using a processor.

ウォールフィルタ機能31は、血流信号にもとづいて生成されたドプラ信号から、血流以外の臓器等の動きに由来するクラッタ成分を除去する。 The wall filter function 31 removes clutter components that originate from the movement of organs, etc. other than blood flow, from the Doppler signal generated based on the blood flow signal.

自己相関機能32は、ドプラ信号にもとづいてドプラデータ(血流の速度、パワー、分散などの血流動態情報を多点について抽出したデータ)を求める。 The autocorrelation function 32 obtains Doppler data (data extracted from multiple points on hemodynamic information such as blood flow velocity, power, and dispersion) based on the Doppler signal.

ヒストグラム生成機能33は、関心領域における血流信号のパワーのヒストグラムを生成する。本実施形態において、「関心領域」はヒストグラム生成に係る処理対象領域を意味するものとする。関心領域は、血流信号を得るための設定領域内に定められ、当該設定領域と同一の領域としてもよい。関心領域は、血流が存在する領域と血流が存在しない領域の両者を含むように設定される。 The histogram generation function 33 generates a histogram of the power of the blood flow signal in the region of interest. In this embodiment, the "region of interest" refers to the region to be processed for histogram generation. The region of interest is defined within a set region for obtaining the blood flow signal, and may be the same region as the set region. The region of interest is set to include both regions where blood flow exists and regions where blood flow does not exist.

また、ヒストグラム生成機能33は、パワーのヒストグラムに代替的にまたは追加的に、血流の分散成分のヒストグラムを生成してもよい。 The histogram generation function 33 may also generate a histogram of the variance component of blood flow instead of or in addition to the power histogram.

閾値決定機能34は、ヒストグラム生成機能33が生成したヒストグラムにもとづいて、被検体の血流の速度に関する超音波ドプラ画像のデータをディスプレイ22に表示させる際のパワーの閾値および分散の閾値の少なくとも一方を決定する。また、閾値決定機能34は、パワーが低いデータほど、かつ、分散が高いデータほど、ブランク処理において削除対象となるように、パワーと分散の両者に関連付けられた値として閾値を決定してもよい。 The threshold determination function 34 determines at least one of the power threshold and the variance threshold when displaying ultrasonic Doppler image data related to the subject's blood flow velocity on the display 22 based on the histogram generated by the histogram generation function 33. The threshold determination function 34 may also determine the threshold as a value associated with both the power and the variance so that data with lower power and data with higher variance are subject to deletion in the blanking process.

この閾値は、ブランク機能36によるブランク処理で用いられる。閾値決定機能34は、たとえばパワーのヒストグラムの最頻値、中央値、平均値、確率密度関数などに応じて、パワーの閾値を決定することができる。 This threshold is used in blanking processing by the blanking function 36. The threshold determination function 34 can determine the power threshold according to, for example, the mode, median, mean, probability density function, etc. of the power histogram.

補正機能35は、ユーザの指示や、閾値のフレーム間の変化などに応じて、閾値を補正する。閾値決定機能34および補正機能35は、決定部の一例である。 The correction function 35 corrects the threshold value according to user instructions, changes in the threshold value between frames, etc. The threshold value determination function 34 and the correction function 35 are examples of a determination unit.

ブランク機能36は、関心領域内の観測点のそれぞれについて、決定部により決定された閾値にもとづいてデータを削除するブランク処理を行う。たとえば、閾値決定機能34からパワーの閾値を受けると、ブランク機能36は、パワー値が当該閾値以下のデータを削除する。また、閾値決定機能34から分散の閾値を受けると、ブランク機能36は、分散値が当該閾値以上のデータを削除する。 The blanking function 36 performs blanking processing to delete data for each observation point in the region of interest based on the threshold determined by the determination unit. For example, when a power threshold is received from the threshold determination function 34, the blanking function 36 deletes data whose power value is equal to or less than the threshold. Also, when a variance threshold is received from the threshold determination function 34, the blanking function 36 deletes data whose variance value is equal to or greater than the threshold.

また、閾値決定機能34からパワーと分散の両者に関連付けられた値としての閾値を受けると、図12を用いて後述するように、ブランク機能36は、パワーと分散の両方の値にもとづいてブランク処理を行う。 In addition, when the threshold value determination function 34 receives a threshold value associated with both power and variance, the blanking function 36 performs blanking processing based on both the power and variance values, as described below with reference to FIG. 12.

画像生成回路14は、ドプラ処理回路13による閾値にもとづくブランク処理後のドプラデータにもとづいて、被検体の血流の速度に関する超音波ドプラ画像のデータを生成し、ディスプレイ22に表示させる。また、画像生成回路14は、生成した画像データを画像サーバ40に保管させてもよい。 The image generation circuit 14 generates ultrasound Doppler image data relating to the subject's blood flow velocity based on the Doppler data after blank processing based on a threshold value by the Doppler processing circuit 13, and displays the data on the display 22. The image generation circuit 14 may also store the generated image data in the image server 40.

なお、ドプラ処理回路13と画像生成回路14とは、超音波診断装置10と接続された医用画像処理装置30に備えられてもよい。この場合、超音波診断装置10は、血流信号を医用画像処理装置30に与え、医用画像処理装置30のドプラ処理回路13と画像生成回路14は、リアルタイムにまたはポストプロセスで、血流信号のパワーおよび分散の少なくとも一方にもとづいてパワーの閾値および分散の閾値の少なくとも一方を決定し、この閾値にもとづいて被検体の血流の速度に関する超音波ドプラ画像のデータを生成してディスプレイ22または医用画像処理装置30のディスプレイに表示させる。この医用画像処理装置30は、画像処理装置の一例である。 The Doppler processing circuit 13 and the image generation circuit 14 may be provided in a medical image processing device 30 connected to the ultrasound diagnostic device 10. In this case, the ultrasound diagnostic device 10 provides a blood flow signal to the medical image processing device 30, and the Doppler processing circuit 13 and the image generation circuit 14 of the medical image processing device 30 determine at least one of a power threshold and a variance threshold based on at least one of the power and variance of the blood flow signal in real time or in post-processing, and generate ultrasound Doppler image data related to the blood flow velocity of the subject based on this threshold and display it on the display 22 or the display of the medical image processing device 30. This medical image processing device 30 is an example of an image processing device.

以下、図3-図11を参照してパワーの閾値を用いる場合の例を説明する。 Below, we will explain an example of using a power threshold with reference to Figures 3 to 11.

図3は、従来のカラードプラ画像の一例を示す説明図である。図3には、頸動脈のカラードプラ像の一例を示した。従来のブランク処理は、事前に設定した固定値をパワーの閾値に用いる。このため、心臓の拍動やミルキングなどの影響によって図3に示すようなブルーミングが発生してしまう。 Figure 3 is an explanatory diagram showing an example of a conventional color Doppler image. Figure 3 shows an example of a color Doppler image of the carotid artery. Conventional blanking processing uses a preset fixed value as the power threshold. For this reason, blooming as shown in Figure 3 occurs due to the effects of cardiac pulsation, milking, etc.

そこで、本実施形態に係るドプラ処理回路13および画像生成回路14は、ブランク処理の閾値を、表示フレームごとに血流信号にもとづいて適応的に決定する。 Therefore, in this embodiment, the Doppler processing circuit 13 and image generation circuit 14 adaptively determine the blanking threshold for each display frame based on the blood flow signal.

図4(a)はパワー値画像の一例を示す説明図であり、(b)はパワー値のヒストグラムの一例を示す説明図である。また、図5は本実施形態に係る画像生成回路14により生成される超音波ドプラ画像の一例を示す説明図である。なお、図3、図4(a)および図5には、関心領域が図示した画像の全体に一致する場合の例を示した。 Figure 4(a) is an explanatory diagram showing an example of a power value image, and (b) is an explanatory diagram showing an example of a histogram of power values. Also, Figure 5 is an explanatory diagram showing an example of an ultrasonic Doppler image generated by the image generation circuit 14 according to this embodiment. Note that Figures 3, 4(a) and 5 show examples in which the region of interest coincides with the entire image shown.

パワー値のヒストグラムは、拍動などの影響によってフレームごとに変化する。このため、最頻値もまた、フレームごとに変化する。図4(b)に示すように、ヒストグラムのピークはパワー値が小さいところにあり、ピークは血流のない観測点のデータにより構成されると考えられる。このため、ブランク処理のパワー値の閾値を固定値とする場合、固定された閾値をヒストグラムのピークのパワー値が超えてしまうと、ブルーミングが広範囲に発生してしまうことがある。 The histogram of power values changes from frame to frame due to the effects of pulsation, etc. For this reason, the mode also changes from frame to frame. As shown in Figure 4(b), the peak of the histogram is at a point where the power value is small, and the peak is thought to be composed of data from observation points where there is no blood flow. For this reason, if the threshold value of the power value for blanking processing is set to a fixed value, blooming may occur over a wide area if the power value of the histogram peak exceeds the fixed threshold.

一方、閾値決定機能34は、血流信号が強く血流が存在する画素のデータを残しつつ血流信号が弱く血流が無い画素のデータを削除することができるよう、関心領域におけるパワーのヒストグラムの最頻値、中央値、平均値、確率密度関数などに応じて、パワーの閾値を決定する。図3と図5に示す頸動脈のカラードプラ画像を比較して明らかなように、このようにパワーのヒストグラムにもとづいて適応的にパワーの閾値を決定することにより、閾値とヒストグラムのピークとの関係が意図せずに逆転することを未然に防ぐことができ、カラードプラモードのブルーミングを大幅に低減することができる。 On the other hand, the threshold determination function 34 determines the power threshold according to the mode, median, mean, probability density function, etc. of the power histogram in the region of interest so that data of pixels with a weak blood flow signal and no blood flow can be deleted while retaining data of pixels with a strong blood flow signal and where blood flow is present. As is clear from a comparison of the color Doppler images of the carotid artery shown in Figures 3 and 5, by adaptively determining the power threshold in this way based on the power histogram, it is possible to prevent the relationship between the threshold and the histogram peak from being unintentionally reversed, and blooming in color Doppler mode can be significantly reduced.

図6は、閾値に関するユーザの指示を受け付ける様子の一例を示す説明図である。図6に示すように、補正機能35は、ブランク調整受付画像51をディスプレイ22またはタッチコマンドスクリーンとして機能する入力インターフェース21のディスプレイに表示し、ユーザの指示に応じて閾値とヒストグラムのピークとの関係を変更してもよい。図6に示す例では、ユーザは、ブランク調整受付画像51を介してブランク処理の強度を選択する。補正機能35は、選択されたブランク処理の強度が強いほど、多くのデータを消すよう閾値を高く補正する。この場合、ユーザはブランク処理の強度を選択するだけで所望のブルーミング緩和処理結果を得ることができる。 Figure 6 is an explanatory diagram showing an example of how a user's instruction regarding the threshold is accepted. As shown in Figure 6, the correction function 35 may display a blank adjustment acceptance image 51 on the display 22 or the display of the input interface 21 functioning as a touch command screen, and change the relationship between the threshold and the histogram peak in response to a user's instruction. In the example shown in Figure 6, the user selects the strength of the blanking process via the blank adjustment acceptance image 51. The correction function 35 corrects the threshold to be higher so that the more data is erased, the stronger the strength of the selected blanking process. In this case, the user can obtain the desired blooming mitigation processing result simply by selecting the strength of the blanking process.

なお、図6に示すようなヒストグラムのビンの数は、固定ではなくユーザの指示にもとづいて変更可能としてもよい。たとえば、ビンの数を増やす場合、1つのビンに割り当てられるパワー幅が狭くなり、より細かい閾値の設定が可能となる。また、ビンの数を減らし、1つのビンに割り当てられるパワー幅を広くしてもよい。 The number of bins in a histogram such as that shown in FIG. 6 may not be fixed, but may be changeable based on user instructions. For example, when the number of bins is increased, the power range assigned to each bin becomes narrower, allowing for more precise threshold setting. Also, the number of bins may be reduced, and the power range assigned to each bin may be wider.

また、ビンの数はユーザの指示にもとづいて変更される場合に限らず、血流信号を得るための設定領域またはヒストグラム生成に係る処理対象領域である関心領域のサイズ、検査領域(検査部位)、血流信号(速度、パワーまたは分散)の時間的もしくは空間的変動の大きさにもとづいて設定されてもよい。たとえば、関心領域のサイズが大きいほどビンの数を増やすこととしてもよい。また、血流信号の変動(たとえば標準偏差)が大きいほどビンの数を増やすこととしてもよい。 The number of bins may not only be changed based on user instructions, but may also be set based on the size of the region of interest, which is the set region for obtaining a blood flow signal or the region to be processed for generating a histogram, the examination region (examination site), and the magnitude of temporal or spatial fluctuation of the blood flow signal (velocity, power, or variance). For example, the number of bins may be increased as the size of the region of interest increases. Also, the number of bins may be increased as the fluctuation of the blood flow signal (e.g., standard deviation) increases.

図7は、ヒストグラムと閾値を示す情報とを同時に表示する場合の例を示す説明図である。ヒストグラムと閾値を示す情報の両方を含む画像52を、たとえば1フレームごとに、あるいは所定フレームごとに離散的に、表示することにより、ユーザは、ヒストグラムと閾値との関係を直感的に容易に把握することができる。また、補正機能35は、この画像52を介してユーザによる閾値の指定指示やヒストグラムのビン数または各ビン幅の指定指示を、ヒストグラム上で直接受け付けてもよい。 Figure 7 is an explanatory diagram showing an example of simultaneously displaying a histogram and information indicating a threshold value. By displaying an image 52 including both a histogram and information indicating a threshold value, for example, for each frame or discretely for each predetermined number of frames, the user can intuitively and easily grasp the relationship between the histogram and the threshold value. Furthermore, the correction function 35 may directly accept, via this image 52, an instruction from the user to specify a threshold value or an instruction to specify the number of bins or the width of each bin in the histogram on the histogram.

図8(a)は関心領域を深さ方向に3つに分割する場合の分割領域を説明するための図であり、(b)はスキャン方向に3つに分割する場合の分割領域を説明するための図である。 Figure 8 (a) is a diagram for explaining the divided regions when the region of interest is divided into three in the depth direction, and (b) is a diagram for explaining the divided regions when the region of interest is divided into three in the scanning direction.

たとえば、位置が深くなると、減衰し血流信号は減衰して弱まる。また、血管の深さや本数が位置によって異なる場合がある。このため、閾値決定機能34は、関心領域を複数に分割し、分割領域ごとにヒストグラムを作成して閾値を決定してもよい。たとえば、図8(a)および(b)に示す例では、閾値決定機能34は、分割領域61、62、63のそれぞれでヒストグラムを作成して閾値を決定するとよい。この場合、分割領域の位置に適した閾値を決定することができるため、より適切にブルーミングを低減することができる。 For example, as the position becomes deeper, the blood flow signal attenuates and weakens. Also, the depth and number of blood vessels may vary depending on the position. For this reason, the threshold determination function 34 may divide the region of interest into multiple regions and create a histogram for each divided region to determine the threshold. For example, in the example shown in Figures 8 (a) and (b), the threshold determination function 34 may create a histogram for each of the divided regions 61, 62, and 63 to determine the threshold. In this case, a threshold appropriate for the position of the divided region can be determined, so blooming can be more appropriately reduced.

また、関心領域の分割の要否と分割のしかた(分割数および分割方向)は、エコー信号の強度、血流信号の速度、パワー、および分散の少なくとも1つに応じて決定してもよい。たとえば、関心領域や分割領域のパワーの分散が所定の上限値より大きいときは、有用なヒストグラムが得られないと判断して分割し、一つの領域内の分散を小さくする。一方、分散が所定の下限値より小さいときは、血流が存在する領域のみまたは血流が存在しない領域のみとなってしまっている可能性があるため、隣接する領域と結合して一つの領域にするなどして領域の範囲を拡大し、分割数を減らす。 In addition, the necessity for dividing the region of interest and the method of division (number of divisions and division direction) may be determined according to at least one of the strength of the echo signal, the speed, power, and variance of the blood flow signal. For example, when the variance of the power of the region of interest or divided region is greater than a predetermined upper limit, it is determined that a useful histogram cannot be obtained and the region is divided to reduce the variance within each region. On the other hand, when the variance is less than a predetermined lower limit, there is a possibility that only regions with blood flow or only regions without blood flow remain, so the range of the region is expanded by combining adjacent regions into one region, for example, and the number of divisions is reduced.

また、関心領域の分割の要否と分割のしかたは、観察部位に応じて決定してもよい。観察部位の情報は、検査オーダ情報から取得してもよいし、ユーザにより入力インターフェース21を介して与えられてもよい。たとえば、深さ方向に広がりがある観察部位は、関心領域を深さ方向に分割するとよい。 In addition, whether or not to divide the region of interest and how to divide it may be determined depending on the observation site. Information on the observation site may be obtained from the examination order information, or may be provided by the user via the input interface 21. For example, for an observation site that extends in the depth direction, it is advisable to divide the region of interest in the depth direction.

また、関心領域内の血管の走行方向の情報が取得可能な場合は、血管の走行方向に関心領域を分割するとよい。血管の走行方向に関心領域を分割することで、容易に、血流が存在する領域と血流が存在しない領域の両者を各分割領域に含ませることができる。 In addition, if information on the direction in which blood vessels run within the region of interest can be obtained, it is advisable to divide the region of interest in the direction in which the blood vessels run. By dividing the region of interest in the direction in which the blood vessels run, it is easy to include both regions where blood flow exists and regions where blood flow does not exist in each divided region.

図9(a)は、深さ方向に並ぶ3つの分割領域ごとにヒストグラムにもとづいて決定された閾値の一例を示す説明図であり、(b)は深さと閾値との関係を事前に規定した複数のカーブの一例を示す説明図である。 Figure 9 (a) is an explanatory diagram showing an example of thresholds determined based on a histogram for each of three divided regions aligned in the depth direction, and (b) is an explanatory diagram showing an example of multiple curves that predefine the relationship between depth and threshold.

一般に、減衰の影響により、深い位置ほどヒストグラムのピークは小さくなっていく。このため、深い位置ほど閾値は小さくなるものと考えられる。そこで、補正機能35は、分割領域ごとに決定された閾値の分布にもとづいて、分割領域ごとに決定された閾値を補正するとよい。具体的には、あらかじめ深さと閾値との関係を規定した補正テーブル(図9(b)参照)と、分割領域ごとに決定した閾値(図9(a)参照)とを比較し、もっとも近いカーブを選択して閾値を補正するとよい。この場合、補正機能35は、関心領域内の3つの分割領域ごとの閾値のそれぞれを補正して補正後の3つの離散的な閾値を利用してもよいし、分割領域の数とは無関係に、選択したカーブにもとづいてデータの位置ごとに連続的に閾値を決定してもよい。 In general, due to the effect of attenuation, the deeper the position, the smaller the histogram peak. For this reason, it is considered that the deeper the position, the smaller the threshold value. Therefore, the correction function 35 may correct the threshold value determined for each divided region based on the distribution of the threshold value determined for each divided region. Specifically, it is preferable to compare a correction table (see FIG. 9(b)) that predefines the relationship between depth and threshold value with the threshold value determined for each divided region (see FIG. 9(a)), select the closest curve, and correct the threshold value. In this case, the correction function 35 may correct each of the threshold values for the three divided regions in the region of interest and use the three discrete threshold values after correction, or may continuously determine the threshold value for each data position based on the selected curve, regardless of the number of divided regions.

図10は、減衰値の測定値に応じて閾値の補正値を決定する場合の例を説明するための図である。超音波診断装置10が超音波の減衰量を測定する機能を有する場合、補正機能35は、エコー信号または血流信号の減衰測定値やその深さ方向積分値(図10左参照)にもとづいて閾値を補正してもよい。この場合、補正機能35は、減衰値が大きいほど閾値が小さくなるよう補正する。また、超音波診断装置10がSN比を測定する機能を有する場合は、補正機能35はSNのカーブにもとづいて閾値を補正してもよい。この場合、補正機能35は、たとえばSNが低く信号の比率が小さいほど閾値が小さくなるよう補正する。 Figure 10 is a diagram for explaining an example of determining a threshold correction value according to a measured value of attenuation. If the ultrasound diagnostic device 10 has a function for measuring the amount of ultrasonic attenuation, the correction function 35 may correct the threshold based on the measured attenuation value of the echo signal or blood flow signal or its depth-direction integral value (see FIG. 10, left). In this case, the correction function 35 corrects the threshold so that the larger the attenuation value, the smaller the threshold value. Also, if the ultrasound diagnostic device 10 has a function for measuring the S/N ratio, the correction function 35 may correct the threshold based on the S/N curve. In this case, the correction function 35 corrects so that the lower the S/N ratio and the smaller the signal ratio, for example, the smaller the threshold value.

また、閾値決定機能34は、ヒストグラムや閾値を複数フレームの情報にもとづいて作成、決定してもよい。たとえば、ヒストグラムは、複数フレームのデータにもとづいて1つ作成し、この1つのヒストグラムにもとづいて複数フレームに1つの閾値を決定してもよい。また、1フレームごとに1つのヒストグラムを作成して仮に閾値を求め、隣接する複数フレームの仮閾値を平均することにより、複数フレームに1つの閾値を決定してもよい。複数フレームの閾値を平均することによって、閾値の変化が緩やかになり、超音波ドプラ画像をフレームごとに連続表示した場合の滑らかさが向上する。 The threshold determination function 34 may also create and determine a histogram or threshold based on information from multiple frames. For example, a histogram may be created based on data from multiple frames, and a threshold for multiple frames may be determined based on this one histogram. A histogram may also be created for each frame to obtain a provisional threshold, and a threshold for multiple frames may be determined by averaging the provisional thresholds for multiple adjacent frames. Averaging the thresholds for multiple frames makes the change in threshold more gradual, improving smoothness when ultrasound Doppler images are displayed continuously for each frame.

また、閾値決定機能34は、複数フレームの仮閾値間の変化に応じて仮閾値を平均するか仮閾値をそのまま用いるかを決定してもよい。また、平均する場合、閾値決定機能34は平均する隣接フレームの数を当該変化に応じてさらに決定してもよい。 The threshold determination function 34 may also determine whether to average the provisional thresholds or use the provisional thresholds as is, depending on the change between the provisional thresholds of multiple frames. In the case of averaging, the threshold determination function 34 may further determine the number of adjacent frames to average, depending on the change.

たとえば、隣接2フレームのそれぞれでヒストグラムにもとづいて仮閾値を求めた場合を考える。この場合、2つの仮閾値があまり変わらないのであれば、平均してもよいと決定する。 For example, consider the case where a provisional threshold is calculated based on the histogram for each of two adjacent frames. In this case, if the two provisional thresholds are not very different, it is decided that they can be averaged.

一方、2つの仮閾値の差が所定値よりも大きい場合は、フレーム間で血流の状況に変化があったと考えられる。このため、特にブランク処理の強度を高めてブルーミングを極力低減させたい場合は、仮閾値をそのまま用いると決定するとよい。仮閾値の差が大きい場合、平均してしまうと、低い仮閾値の影響で高い仮閾値に対応するフレームの閾値が下がってしまうため、高い仮閾値に対応するフレームにおいてブルーミングを除去することができなくなる場合がある。 On the other hand, if the difference between the two provisional thresholds is greater than a specified value, it is believed that there has been a change in the blood flow conditions between frames. For this reason, it is advisable to decide to use the provisional threshold as is, especially when you want to increase the strength of the blanking process to reduce blooming as much as possible. If the difference between the provisional thresholds is large, averaging will lower the threshold of the frame corresponding to the high provisional threshold due to the influence of the low provisional threshold, and it may not be possible to remove blooming in the frame corresponding to the high provisional threshold.

また、閾値決定機能34は、閾値の変化が正(すなわち増加したか)、負(すなわち減少したか)、に応じて仮閾値を平均するか仮閾値をそのまま用いるかを決定してもよい。たとえば、ブランク処理の強度を高めてブルーミングを極力低減させたい場合であって、変化が負の場合には、平均化することで、過去の高い仮閾値の影響で現在の低い仮閾値を高めに補正することができる。一方、ブランク処理の強度を高めてブルーミングを極力低減させたい場合であって、変化が正の場合には、逆に平均しないことによって過去の低い仮閾値の影響を排除することができる。 The threshold determination function 34 may also determine whether to average the provisional thresholds or use them as is depending on whether the change in the threshold is positive (i.e., an increase) or negative (i.e., a decrease). For example, if it is desired to increase the strength of the blanking process to minimize blooming and the change is negative, averaging can be performed to correct the current low provisional threshold to a higher value due to the influence of the past high provisional threshold. On the other hand, if it is desired to increase the strength of the blanking process to minimize blooming and the change is positive, conversely, not averaging can be performed to eliminate the influence of the past low provisional threshold.

図11は、関心領域の位置に応じた仮閾値の平均化処理を説明するための図である。 Figure 11 is a diagram to explain the averaging process of the provisional threshold according to the position of the region of interest.

閾値決定機能34は、観察部位または関心領域の位置に応じて、仮閾値を平均するか仮閾値をそのまま用いるかを決定するとともに、平均する場合は平均する隣接フレームの数を決定してもよい。 The threshold determination function 34 determines whether to average the provisional thresholds or use the provisional thresholds as is depending on the position of the observed area or region of interest, and may also determine the number of adjacent frames to average if averaging is performed.

たとえば、浅い位置に設定された関心領域71では仮閾値を平均せず、深い位置に設定された関心領域72では、血流信号が小さくSNが小さいことから、浅い位置よりも多いフレーム(たとえば隣接3フレーム)の仮閾値を平均するなど、関心領域の位置に応じて平均化処理を異ならせてもよい。 For example, the averaging process may be different depending on the position of the region of interest. For example, the provisional threshold value is not averaged for a region of interest 71 set at a shallow position, whereas for a region of interest 72 set at a deep position, the blood flow signal is small and the S/N ratio is small, so the provisional threshold values of more frames (e.g., three adjacent frames) are averaged compared to a shallow position.

また、たとえば、観察部位が腹部の場合は隣接2フレームの平均、頸動脈の場合は隣接3フレームの平均、大人の心臓の場合は平均しない、など、観察部位に応じて平均化処理を異ならせてもよい。 Also, the averaging process may be different depending on the observation area, for example, averaging two adjacent frames when the observation area is the abdomen, averaging three adjacent frames when the observation area is the carotid artery, and no averaging when the observation area is the adult heart.

以上、図3-11を参照してパワーの閾値を用いる場合の例について説明したが、上述の通り、本実施形態に係るドプラ処理回路13および画像生成回路14、ならびにこれらを備えた医用画像処理装置30は、分散の閾値を用いることもできる。 The above describes an example of using a power threshold with reference to Figures 3-11, but as mentioned above, the Doppler processing circuit 13 and image generation circuit 14 according to this embodiment, and the medical image processing device 30 equipped with these, can also use a variance threshold.

図12は、パワーと分散の両者にもとづいて閾値を決定するためのテーブルの一例を示す説明図である。 Figure 12 is an explanatory diagram showing an example of a table for determining a threshold value based on both power and variance.

閾値決定機能34は、パワーと分散の両者にもとづいて閾値を決定してもよい。たとえば、パワーおよび分散の両者と、閾値と、をあらかじめ関連付けたテーブルを記憶回路16に記憶させておく場合(図12参照)、閾値決定機能34は、このテーブルを参照してパワーと分散の両者にもとづいて閾値を決定することができる。このテーブルは、図12に示すように、パワーが低いデータほど、かつ、分散が高いデータほど、ブランク処理において削除対象となり非表示となるように設定されるとよい。また、閾値決定機能34は、パワーと分散の少なくとも一方と、血流速度とにもとづいて閾値を決定してもよい。 The threshold determination function 34 may determine the threshold based on both power and variance. For example, if a table that associates both power and variance with a threshold is stored in advance in the memory circuitry 16 (see FIG. 12), the threshold determination function 34 can refer to this table to determine the threshold based on both power and variance. As shown in FIG. 12, this table is preferably set so that data with lower power and higher variance is subject to deletion and hidden in the blanking process. The threshold determination function 34 may also determine the threshold based on at least one of power and variance, and blood flow velocity.

以上説明した少なくとも1つの実施形態によれば、カラードプラモードのブルーミングを低減することができる。 At least one of the embodiments described above can reduce blooming in color Doppler mode.

なお、上記実施形態において、「プロセッサ」という文言は、たとえば、専用または汎用の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))等の回路を意味する。プロセッサがたとえばCPUである場合、プロセッサは記憶回路に保存されたプログラムを読み出して実行することにより、各種機能を実現する。また、プロセッサがたとえばASICである場合、記憶回路にプログラムを保存するかわりに、当該プログラムに相当する機能がプロセッサの回路内に論理回路として直接組み込まれる。この場合、プロセッサは回路内に組み込まれたプログラムを読み出し実行するハードウェア処理により各種機能を実現する。あるいはまた、プロセッサは、ソフトウェア処理とハードウェア処理とを組み合わせて各種機能を実現することもできる。 In the above embodiment, the term "processor" refers to a circuit such as a dedicated or general-purpose CPU (Central Processing Unit), GPU (Graphics Processing Unit), or 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)). When the processor is a CPU, for example, the processor realizes various functions by reading and executing a program stored in a memory circuit. When the processor is an ASIC, for example, instead of storing a program in a memory circuit, a function corresponding to the program is directly built into the processor circuit as a logic circuit. In this case, the processor realizes various functions by hardware processing that reads and executes the program built into the circuit. Alternatively, the processor can realize various functions by combining software processing and hardware processing.

また、上記実施形態では処理回路の単一のプロセッサが各機能を実現する場合の例について示したが、複数の独立したプロセッサを組み合わせて処理回路を構成し、各プロセッサが各機能を実現してもよい。また、プロセッサが複数設けられる場合、プログラムを記憶する記憶回路は、プロセッサごとに個別に設けられてもよいし、1つの記憶回路が全てのプロセッサの機能に対応するプログラムを一括して記憶してもよい。 In addition, in the above embodiment, an example was shown in which a single processor of the processing circuit realizes each function, but a processing circuit may be configured by combining multiple independent processors, and each processor may realize each function. In addition, when multiple processors are provided, a memory circuit for storing programs may be provided separately for each processor, or a single memory circuit may collectively store programs corresponding to the functions of all the processors.

いくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更、実施形態同士の組み合わせを行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。 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.

10 超音波診断装置
13 ドプラ処理回路
14 画像生成回路
16 記憶回路
21 入力インターフェース
22 ディスプレイ
30 医用画像処理装置
34 閾値決定機能
35 補正機能
36 ブランク機能
10 Ultrasound diagnostic device 13 Doppler processing circuit 14 Image generation circuit 16 Memory circuit 21 Input interface 22 Display 30 Medical image processing device 34 Threshold determination function 35 Correction function 36 Blank function

Claims (20)

被検体の血流信号のパワーおよび分散の少なくとも一方にもとづいて、前記被検体の血流の速度に関する超音波ドプラ画像のデータを表示部に表示させる際の前記パワーの閾値および前記分散の閾値の少なくとも一方を決定する決定部と、
前記閾値に基づいて前記超音波ドプラ画像のデータの少なくとも一部を補正するブランク処理を行うブランク処理部と、
前記ブランク処理後のデータに基づいて前記被検体の血流の速度に関する超音波ドプラ画像のデータを生成して表示部に表示させる画像生成部と、
を備え、
前記決定部は、
血流信号を得るための設定領域内に定めた関心領域における前記パワーおよび前記分散の少なくとも一方のヒストグラムにもとづいて前記閾値を決定するとともに、前記ヒストグラムのビンの数を、前記設定領域の大きさ、前記関心領域の大きさ、エコー信号の強度、血流信号の速度、血流信号のパワー、および血流信号の分散の少なくとも1つに応じて変更する、
音波診断装置。
a determination unit that determines at least one of a power threshold value and a variance threshold value when displaying ultrasonic Doppler image data related to a blood flow velocity of the subject on a display unit, based on at least one of a power and a variance of a blood flow signal of the subject;
a blanking processing unit that performs blanking processing to correct at least a part of the data of the ultrasonic Doppler image based on the threshold value;
an image generating unit that generates ultrasonic Doppler image data relating to the blood flow velocity of the subject based on the blanked data and displays the image data on a display unit;
Equipped with
The determination unit is
determining the threshold value based on a histogram of at least one of the power and the variance in a region of interest defined within a set region for obtaining a blood flow signal, and changing the number of bins of the histogram in response to at least one of a size of the set region, a size of the region of interest, an intensity of an echo signal, a velocity of the blood flow signal, a power of the blood flow signal, and a variance of the blood flow signal;
Ultrasound diagnostic equipment.
前記決定部は、
表示フレームごとに、前記ヒストグラムを作成して前記閾値を決定することを繰り返す、
請求項に記載の超音波診断装置。
The determination unit is
Repeating the process of creating the histogram and determining the threshold for each display frame.
The ultrasonic diagnostic apparatus according to claim 1 .
前記決定部は、
前記ヒストグラムの最頻値、中央値、平均値、確率密度関数の少なくとも1つに応じて前記閾値を決定する、
請求項またはに記載の超音波診断装置。
The determination unit is
determining the threshold value according to at least one of a mode, a median, a mean, and a probability density function of the histogram;
3. The ultrasonic diagnostic apparatus according to claim 1 .
前記決定部は、
前記閾値を、ユーザの指示に応じて補正して前記ブランク処理部に与える、
請求項ないしのいずれか1項に記載の超音波診断装置。
The determination unit is
the threshold value is corrected in response to a user instruction and provided to the blanking processing unit;
4. The ultrasonic diagnostic apparatus according to claim 1 .
前記決定部は、
前記ヒストグラムのビンの数を、ユーザの指示に応じて変更する、
請求項ないしのいずれか1項に記載の超音波診断装置。
The determination unit is
Varying the number of bins of the histogram in response to a user instruction.
5. The ultrasonic diagnostic apparatus according to claim 1 .
前記画像生成部は、
前記ヒストグラムと前記閾値を示す情報とを前記表示部に表示させる、
請求項ないしのいずれか1項に記載の超音波診断装置。
The image generating unit includes:
displaying the histogram and information indicating the threshold on the display unit;
6. The ultrasonic diagnostic apparatus according to claim 1 .
前記決定部は、
前記関心領域を複数に分割し、分割領域ごとに、前記ヒストグラムを作成して前記閾値を決定する、
請求項ないしのいずれか1項に記載の超音波診断装置。
The determination unit is
Dividing the region of interest into a plurality of regions, creating the histogram for each divided region, and determining the threshold value;
7. The ultrasonic diagnostic apparatus according to claim 1 .
前記決定部は、
エコー信号の強度、血流信号の速度、パワー、および分散の少なくとも1つに応じて、前記関心領域の分割の要否と、分割数および分割方向とを決定する、
請求項に記載の超音波診断装置。
The determination unit is
determining whether or not the region of interest needs to be divided, and the number and direction of division of the region of interest, depending on at least one of the intensity of the echo signal and the velocity, power, and variance of the blood flow signal;
The ultrasonic diagnostic apparatus according to claim 7 .
前記決定部は、
観察部位に応じて、前記関心領域の分割の要否と、分割数および分割方向とを決定する、
請求項に記載の超音波診断装置。
The determination unit is
determining whether or not the region of interest is to be divided, and the number and direction of division of the region of interest according to the observation site;
The ultrasonic diagnostic apparatus according to claim 7 .
前記決定部は、
前記関心領域内の血管の走行方向に前記関心領域を分割し、分割領域ごとに、前記ヒストグラムを作成して前記閾値を決定する、
請求項ないしのいずれか1項に記載の超音波診断装置。
The determination unit is
Dividing the region of interest in a running direction of a blood vessel within the region of interest, creating the histogram for each divided region, and determining the threshold value;
10. The ultrasonic diagnostic apparatus according to claim 1 .
前記決定部は、
前記分割領域ごとに決定された閾値の分布にもとづいて前記分割領域ごとに決定された閾値を補正する、
請求項ないし10のいずれか1項に記載の超音波診断装置。
The determination unit is
correcting the threshold determined for each of the divided regions based on a distribution of the threshold determined for each of the divided regions;
11. The ultrasonic diagnostic apparatus according to claim 7 .
前記決定部は、
エコー信号または血流信号のSN測定値、またはエコー信号または血流信号の減衰測定値、にもとづいて前記閾値を補正する、
請求項ないし11のいずれか1項に記載の超音波診断装置。
The determination unit is
correcting the threshold value based on an SNR measurement value of the echo signal or blood flow signal, or an attenuation measurement value of the echo signal or blood flow signal;
The ultrasonic diagnostic apparatus according to any one of claims 1 to 11 .
前記決定部は、
前記関心領域のうち深い位置ほどパワーの前記閾値が小さくなるように補正する、
請求項1ないし12のいずれか1項に記載の超音波診断装置。
The determination unit is
correcting the power threshold so that the power threshold becomes smaller as the position of the region of interest becomes deeper;
13. The ultrasonic diagnostic apparatus according to claim 1.
前記決定部は、
複数フレームのデータにもとづいて1つの前記ヒストグラムを作成して前記閾値を決定する、
請求項ないし13のいずれか1項に記載の超音波診断装置。
The determination unit is
creating one histogram based on data of a plurality of frames to determine the threshold value;
14. The ultrasonic diagnostic apparatus according to claim 1 .
前記決定部は、
1フレームごとに前記ヒストグラムを作成して仮に閾値を求め、隣接する複数フレームの仮閾値を平均することにより前記閾値を決定する、
請求項ないし13のいずれか1項に記載の超音波診断装置。
The determination unit is
A histogram is created for each frame to determine a provisional threshold value, and the provisional threshold value is determined by averaging the provisional threshold values of a plurality of adjacent frames.
14. The ultrasonic diagnostic apparatus according to claim 1 .
前記決定部は、
前記複数フレームの仮閾値間の変化に応じて、仮閾値を平均するか仮閾値をそのまま用いるかを決定するとともに、平均する場合は平均する隣接フレームの数を決定する、
請求項15に記載の超音波診断装置。
The determination unit is
determining whether to average the provisional thresholds or to use the provisional thresholds as they are according to the change between the provisional thresholds of the plurality of frames, and, if averaging, determining the number of adjacent frames to be averaged;
The ultrasonic diagnostic apparatus according to claim 15 .
前記決定部は、
観察部位または前記関心領域の位置に応じて、仮閾値を平均するか仮閾値をそのまま用いるかを決定するとともに、平均する場合は平均する隣接フレームの数を決定する、
請求項14または15に記載の超音波診断装置。
The determination unit is
determining whether to average the provisional thresholds or use the provisional thresholds as they are according to the position of the observation site or the region of interest, and, if averaging, determining the number of adjacent frames to average;
16. The ultrasonic diagnostic apparatus according to claim 14 or 15 .
前記決定部は、
血流の速度と、前記パワーおよび前記分散の少なくとも一方と、にもとづいて前記閾値を決定する、
請求項1ないし17のいずれか1項に記載の超音波診断装置。
The determination unit is
determining the threshold value based on a blood flow velocity and at least one of the power and the variance;
18. The ultrasonic diagnostic apparatus according to claim 1.
被検体の血流信号のパワーおよび分散の少なくとも一方にもとづいて、前記被検体の血流の速度に関する超音波ドプラ画像のデータを表示部に表示させる際の前記パワーの閾値および前記分散の閾値の少なくとも一方を決定する決定部と、
前記閾値に基づいて前記超音波ドプラ画像のデータの少なくとも一部を補正するブランク処理を行うブランク処理部と、
前記ブランク処理後のデータに基づいて前記被検体の血流の速度に関する超音波ドプラ画像のデータを表示部に表示させる画像生成部と、
を備え、
前記決定部は、
血流信号を得るための設定領域内に定めた関心領域における前記パワーおよび前記分散の少なくとも一方のヒストグラムにもとづいて前記閾値を決定するとともに、前記ヒストグラムのビンの数を、前記設定領域の大きさ、前記関心領域の大きさ、エコー信号の強度、血流信号の速度、血流信号のパワー、および血流信号の分散の少なくとも1つに応じて変更する、
医用画像処理装置
a determination unit that determines at least one of a power threshold value and a variance threshold value when displaying ultrasonic Doppler image data related to a blood flow velocity of the subject on a display unit, based on at least one of a power and a variance of a blood flow signal of the subject;
a blanking processing unit that performs blanking processing to correct at least a part of the data of the ultrasonic Doppler image based on the threshold value;
an image generating unit that displays ultrasonic Doppler image data relating to the blood flow velocity of the subject on a display unit based on the data after the blanking process;
Equipped with
The determination unit is
determining the threshold value based on a histogram of at least one of the power and the variance in a region of interest defined within a set region for obtaining a blood flow signal, and changing the number of bins of the histogram in response to at least one of a size of the set region, a size of the region of interest, an intensity of an echo signal, a velocity of the blood flow signal, a power of the blood flow signal, and a variance of the blood flow signal;
Medical imaging equipment .
コンピュータに、
被検体の血流信号のパワーおよび分散の少なくとも一方にもとづいて、前記被検体の血流の速度に関する超音波ドプラ画像のデータを表示部に表示させる際の前記パワーの閾値および前記分散の閾値の少なくとも一方を決定するステップと、
前記閾値に基づいて前記超音波ドプラ画像のデータの少なくとも一部を補正するブランク処理を行うステップと、
前記ブランク処理後のデータに基づいて前記被検体の血流の速度に関する超音波ドプラ画像のデータを表示部に表示させるステップと、
を実行させるための医用画像処理プログラムであって、
前記パワーの閾値および前記分散の閾値の少なくとも一方を決定するステップは、
血流信号を得るための設定領域内に定めた関心領域における前記パワーおよび前記分散の少なくとも一方のヒストグラムにもとづいて前記閾値を決定するステップと、
前記ヒストグラムのビンの数を、前記設定領域の大きさ、前記関心領域の大きさ、エコー信号の強度、血流信号の速度、血流信号のパワー、および血流信号の分散の少なくとも1つに応じて変更するステップと、
を含む、
医用画像処理プログラム
On the computer,
determining at least one of a power threshold and a variance threshold when displaying ultrasonic Doppler image data related to the velocity of blood flow in the subject on a display unit, based on at least one of the power and variance of a blood flow signal in the subject;
a step of performing blanking processing to correct at least a part of the data of the ultrasonic Doppler image based on the threshold value;
displaying ultrasonic Doppler image data relating to the blood flow velocity of the subject on a display unit based on the blanked data;
A medical image processing program for executing
The step of determining at least one of the power threshold and the variance threshold comprises:
determining the threshold value based on a histogram of at least one of the power and the variance in a region of interest defined within a set region for obtaining a blood flow signal ;
changing the number of bins of the histogram in response to at least one of the size of the set region, the size of the region of interest, the intensity of the echo signal, the velocity of the blood flow signal, the power of the blood flow signal, and the variance of the blood flow signal ;
including,
A medical image processing program .
JP2019238036A 2019-12-27 2019-12-27 Ultrasound diagnostic device, medical image processing device, and medical image processing program Active JP7469877B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2019238036A JP7469877B2 (en) 2019-12-27 2019-12-27 Ultrasound diagnostic device, medical image processing device, and medical image processing program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2019238036A JP7469877B2 (en) 2019-12-27 2019-12-27 Ultrasound diagnostic device, medical image processing device, and medical image processing program

Publications (2)

Publication Number Publication Date
JP2021104301A JP2021104301A (en) 2021-07-26
JP7469877B2 true JP7469877B2 (en) 2024-04-17

Family

ID=76919527

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019238036A Active JP7469877B2 (en) 2019-12-27 2019-12-27 Ultrasound diagnostic device, medical image processing device, and medical image processing program

Country Status (1)

Country Link
JP (1) JP7469877B2 (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000237192A (en) 1999-02-22 2000-09-05 Ge Yokogawa Medical Systems Ltd B-cfm image display method and ultrasonic diagnostic device
US20020186868A1 (en) 2001-06-12 2002-12-12 Steinar Bjaerum Ultrasound color characteristic mapping
JP2008154891A (en) 2006-12-26 2008-07-10 Ge Medical Systems Global Technology Co Llc Color doppler apparatus and ultrasonic imaging apparatus
JP2009261436A (en) 2008-04-22 2009-11-12 Aloka Co Ltd Ultrasonic diagnostic apparatus
JP2014198240A (en) 2013-03-12 2014-10-23 株式会社東芝 Ultrasonic diagnostic apparatus and ultrasonic image processing program
US20160199022A1 (en) 2015-01-13 2016-07-14 Samsung Medison Co., Ltd. Ultrasound diagnosis apparatus and method of operating the same
JP2019092672A (en) 2017-11-20 2019-06-20 コニカミノルタ株式会社 Ultrasonic signal processing device, ultrasonic diagnostic device, and ultrasonic signal processing method
JP2019181183A (en) 2018-03-30 2019-10-24 キヤノンメディカルシステムズ株式会社 Medical diagnostic apparatus, medical image processing apparatus, and image processing program

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000237192A (en) 1999-02-22 2000-09-05 Ge Yokogawa Medical Systems Ltd B-cfm image display method and ultrasonic diagnostic device
US20020186868A1 (en) 2001-06-12 2002-12-12 Steinar Bjaerum Ultrasound color characteristic mapping
JP2008154891A (en) 2006-12-26 2008-07-10 Ge Medical Systems Global Technology Co Llc Color doppler apparatus and ultrasonic imaging apparatus
JP2009261436A (en) 2008-04-22 2009-11-12 Aloka Co Ltd Ultrasonic diagnostic apparatus
JP2014198240A (en) 2013-03-12 2014-10-23 株式会社東芝 Ultrasonic diagnostic apparatus and ultrasonic image processing program
US20160199022A1 (en) 2015-01-13 2016-07-14 Samsung Medison Co., Ltd. Ultrasound diagnosis apparatus and method of operating the same
JP2019092672A (en) 2017-11-20 2019-06-20 コニカミノルタ株式会社 Ultrasonic signal processing device, ultrasonic diagnostic device, and ultrasonic signal processing method
JP2019181183A (en) 2018-03-30 2019-10-24 キヤノンメディカルシステムズ株式会社 Medical diagnostic apparatus, medical image processing apparatus, and image processing program

Also Published As

Publication number Publication date
JP2021104301A (en) 2021-07-26

Similar Documents

Publication Publication Date Title
US10278670B2 (en) Ultrasound diagnostic apparatus and method of controlling ultrasound diagnostic apparatus
JP6000569B2 (en) Ultrasonic diagnostic apparatus and control program
US10335118B2 (en) Ultrasonic diagnostic apparatus, medical image processing apparatus, and medical image parallel display method
JP6188594B2 (en) Ultrasonic diagnostic apparatus, image processing apparatus, and image processing method
KR102011545B1 (en) Ultrasonic diagnosis apparatus, image processing apparatus, and image processing method
US10835200B2 (en) Ultrasound diagnosis apparatus
JP2017093913A (en) Ultrasonic diagnostic device, signal processing device and analysis program
US8724880B2 (en) Ultrasonic diagnostic apparatus and medical image processing apparatus
CN110575198B (en) Analysis device and analysis method
US20220330920A1 (en) Ultrasonic diagnostic apparatus and medical image processing apparatus
JP7469877B2 (en) Ultrasound diagnostic device, medical image processing device, and medical image processing program
JP5823184B2 (en) Ultrasonic diagnostic apparatus, medical image processing apparatus, and medical image processing program
JP6651405B2 (en) Ultrasound diagnostic apparatus and program
US20210085281A1 (en) Ultrasonic diagnostic apparatus and image processing apparatus
US20200337671A1 (en) Ultrasound diagnosis apparatus and analyzing apparatus
JP7345374B2 (en) Ultrasound diagnostic equipment
JP2012143358A (en) Ultrasonic diagnostic equipment and program
US11690597B2 (en) Ultrasonic diagnostic apparatus
US20230225710A1 (en) Ultrasonic diagnostic apparatus and image processing apparatus
US10709421B2 (en) Ultrasound diagnostic apparatus
JP6793502B2 (en) Ultrasonic diagnostic equipment
US20190298314A1 (en) Ultrasonic diagnostic device, medical image processing device, and medical image processing method
JP2020185122A (en) Analyzing apparatus and ultrasonic diagnostic apparatus

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20220830

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20230329

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230425

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230626

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230905

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20231026

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20231222

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: 20240312

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240405