JP4726407B2

JP4726407B2 - Ultrasonic diagnostic equipment

Info

Publication number: JP4726407B2
Application number: JP2003298673A
Authority: JP
Inventors: 俊一郎谷川; 宏章片岡
Original assignee: GE Medical Systems Global Technology Co LLC
Current assignee: GE Medical Systems Global Technology Co LLC
Priority date: 2003-08-22
Filing date: 2003-08-22
Publication date: 2011-07-20
Anticipated expiration: 2023-08-22
Also published as: JP2005065936A

Description

本発明は、空間フィルタリング（ｆｉｌｔｅｒｉｎｇ）方法および空間フィルタ（ｆｉｌｔｅｒ）並びに超音波診断装置に関し、特に、フロー（ｆｌｏｗ）画像について空間フィルタリングを行う方法および空間フィルタ、並びに、そのような空間フィルタを有する超音波診断装置に関する。 The present invention relates to a spatial filtering method, a spatial filter, and an ultrasonic diagnostic apparatus, and more particularly, to a method and a spatial filter for performing spatial filtering on a flow image, and an ultrasonic having such a spatial filter. The present invention relates to an ultrasonic diagnostic apparatus.

超音波診断装置では、超音波エコー（ｅｃｈｏ）のドップラシフト（Ｄｏｐｐｌｅｒｓｈｉｆｔ）を利用して、血流等に関するフロー画像を撮影することが行われる。フロー画像を得るために、エコー受信信号にはＭＴＩ（ｍｏｖｉｎｇｔａｒｇｅｔｉｎｄｉｃａｔｉｏｎ）処理が行われ、さらに自己相関演算が行われる（例えば、特許文献１参照）。
特開２００１−１２８９７６号公報（第５−６頁、図１，７） In the ultrasonic diagnostic apparatus, a flow image relating to blood flow or the like is taken using a Doppler shift of ultrasonic echo (echo). In order to obtain a flow image, the echo reception signal is subjected to MTI (moving target indication) processing and further subjected to autocorrelation calculation (see, for example, Patent Document 1).
JP 2001-128976 A (page 5-6, FIGS. 1 and 7)

フロー画像には、本来フローがあるべきところがあたかもフローが無いように表示される現象すなわちいわゆる黒抜けが生じることがある。これはデータ処理の過程で何らかの異常によりデータが欠落することに由来する。 In a flow image, there may be a phenomenon that a flow should be originally displayed as if there is no flow, that is, a so-called black spot. This is because data is lost due to some abnormality in the process of data processing.

そこで、本発明の課題は、フロー画像の黒抜けを補正する空間フィルタリング方法および空間フィルタ、並びに、そのような空間フィルタを有する超音波診断装置を実現することである。 Accordingly, an object of the present invention is to realize a spatial filtering method and a spatial filter that correct black spots in a flow image, and an ultrasonic diagnostic apparatus having such a spatial filter.

（１）上記の課題を解決するためのひとつの観点での発明は、フロー画像について、予め定められた条件に当てはまる画素値を持つ画素の位置および連続画素数を求め、前記連続画素数が予め定められた限度以内のときは画素が存在する領域を補正領域と判定し、前記補正領域の画素値を周囲の画素値から補間した画素値で置き換える、ことを特徴とする空間フィルタリング方法である。 (1) In one aspect of the invention for solving the above-described problem, a flow image is obtained by obtaining a position of a pixel having a pixel value that satisfies a predetermined condition and the number of continuous pixels, and the number of continuous pixels is determined in advance. A spatial filtering method is characterized in that when it is within a predetermined limit, an area in which a pixel exists is determined as a correction area, and a pixel value in the correction area is replaced with a pixel value interpolated from surrounding pixel values.

（２）上記の課題を解決するための他の観点での発明は、フロー画像について、予め定められた条件に当てはまる画素値を持つ画素の位置および連続画素数を求める探索手段と、前記連続画素数が予め定められた限度以内のときは画素が存在する領域を補正領域と判定する判定手段と、前記補正領域の画素値を周囲の画素値から補間した画素値で置き換える補正手段と、を具備することを特徴とする空間フィルタである。 (2) According to another aspect of the invention for solving the above-described problem, a search unit for obtaining a position of a pixel having a pixel value satisfying a predetermined condition and the number of continuous pixels in the flow image, and the continuous pixel A determination unit that determines an area in which a pixel is present as a correction area when the number is within a predetermined limit; and a correction unit that replaces a pixel value of the correction area with a pixel value interpolated from surrounding pixel values. It is a spatial filter characterized by doing.

（３）上記の課題を解決するための他の観点での発明は、対象に照射した超音波のエコーに基づいてフロー画像を撮影する撮影手段と、前記フロー画像を処理する空間フィルタと、を有する超音波診断装置であって、前記空間フィルタは、フロー画像について、予め定められた条件に当てはまる画素値を持つ画素の位置および連続画素数を求める探索手段と、前記連続画素数が予め定められた限度以内のときは画素が存在する領域を補正領域と判定する判定手段と、前記補正領域の画素値を周囲の画素値から補間した画素値で置き換える補正手段と、を具備することを特徴とする超音波診断装置である。 (3) In another aspect of the invention for solving the above-described problem, an imaging unit that captures a flow image based on an echo of an ultrasonic wave applied to a target, and a spatial filter that processes the flow image, An ultrasonic diagnostic apparatus, wherein the spatial filter includes a search unit that obtains the position of a pixel having a pixel value that satisfies a predetermined condition and the number of continuous pixels in the flow image, and the number of continuous pixels is predetermined. A determination unit that determines a region in which a pixel is present as a correction region, and a correction unit that replaces a pixel value of the correction region with a pixel value interpolated from surrounding pixel values. This is an ultrasonic diagnostic apparatus.

前記条件は閾値によって規定されることが、条件設定が容易な点で好ましい。前記限度は可変であることが、黒抜け補正を適切に行う点で好ましい。 The condition is preferably defined by a threshold value in terms of easy condition setting. It is preferable that the limit is variable in terms of appropriately performing blackout correction.

本発明では、フロー画像について、予め定められた条件に当てはまる画素値を持つ画素の位置および連続画素数を求め、連続画素数が予め定められた限度以内のときは画素が存在する領域を補正領域と判定し、補正領域の画素値を周囲の画素値から補間した画素値で置き換えるので、フロー画像の黒抜けを補正する空間フィルタリング方法および空間フィルタ、並びに、そのような空間フィルタを有する超音波診断装置を実現することができる。本発明による黒抜け補正は、黒抜け領域の判定を行わない通常の空間フィルタリングのように、空間分解能やコントラスト（ｃｏｎｔｒａｓｔ）を低下させることがなく、また、画像輪郭の拡大を生じるということもない。 In the present invention, for a flow image, the position of a pixel having a pixel value that satisfies a predetermined condition and the number of continuous pixels are obtained, and when the number of continuous pixels is within a predetermined limit, the region where the pixel exists is a correction region. And the pixel value in the correction area is replaced with the pixel value interpolated from the surrounding pixel values, so that the spatial filtering method and the spatial filter for correcting the black omission in the flow image, and the ultrasonic diagnosis having such a spatial filter An apparatus can be realized. The blackout correction according to the present invention does not reduce the spatial resolution and contrast and does not cause the enlargement of the image contour unlike the normal spatial filtering that does not determine the blackout area. .

以下、図面を参照して発明を実施するための最良の形態を詳細に説明する。なお、本発明は、発明を実施するための最良の形態に限定されるものではない。図１に超音波診断装置のブロック（ｂｌｏｃｋ）図を示す。本装置は発明を実施するための最良の形態の一例である。本装置の構成によって、超音波診断装置に関する本発明を実施するための最良の形態の一例が示される。 The best mode for carrying out the invention will be described below in detail with reference to the drawings. Note that the present invention is not limited to the best mode for carrying out the invention. FIG. 1 shows a block diagram of the ultrasonic diagnostic apparatus. This apparatus is an example of the best mode for carrying out the invention. An example of the best mode for carrying out the present invention relating to an ultrasonic diagnostic apparatus is shown by the configuration of the present apparatus.

図１に示すように、本装置は、超音波プローブ２（ｐｒｏｂｅ）を有する。超音波プローブ２は、図示しない複数の超音波トランスデューサ（ｔｒａｎｓｄｕｃｅｒ）のアレイ（ａｒｒａｙ）を有する。個々の超音波トランスデューサは例えばＰＺＴ（チタン（Ｔｉ）酸ジルコン（Ｚｒ）酸鉛）セラミックス（ｃｅｒａｍｉｃｓ）等の圧電材料によって構成される。超音波プローブ２は、操作者により撮影の対象４に当接して使用される。 As shown in FIG. 1, this apparatus has an ultrasonic probe 2 (probe). The ultrasonic probe 2 has an array of a plurality of ultrasonic transducers (not shown). Each ultrasonic transducer is made of a piezoelectric material such as PZT (titanium (Ti) zirconate (Zr) acid) ceramics. The ultrasonic probe 2 is used in contact with an object 4 to be photographed by an operator.

超音波プローブ２は送受信部６に接続されている。送受信部６は、超音波プローブ２に駆動信号を与えて超音波を送波させる。送受信部６は、また、超音波プローブ２が受波したエコー信号を受信する。 The ultrasonic probe 2 is connected to the transmission / reception unit 6. The transmission / reception unit 6 sends a drive signal to the ultrasonic probe 2 to transmit ultrasonic waves. The transmission / reception unit 6 also receives an echo signal received by the ultrasonic probe 2.

送受信部６のブロック図を図２に示す。同図に示すように、送受信部６は送波タイミング（ｔｉｍｉｎｇ）発生ユニット（ｕｎｉｔ）６０２を有する。送波タイミング発生ユニット６０２は、送波タイミング信号を周期的に発生して送波ビームフォーマ（ｂｅａｍｆｏｒｍｅｒ）６０４に入力する。送波タイミング信号の周期は後述の制御部１８により制御される。 A block diagram of the transceiver 6 is shown in FIG. As shown in the figure, the transmission / reception unit 6 includes a transmission timing generation unit (unit) 602. The transmission timing generation unit 602 periodically generates a transmission timing signal and inputs the transmission timing signal to a transmission beamformer 604. The period of the transmission timing signal is controlled by the control unit 18 described later.

送波ビームフォーマ６０４は、送波のビームフォーミング（ｂｅａｍｆｏｒｍｉｎｇ）を行うもので、送波タイミング信号に基づき、所定の方位の超音波ビームを形成するためのビームフォーミング信号を生じる。ビームフォーミング信号は、方位に対応した時間差が付与された複数の駆動信号からなる。ビームフォーミングは後述の制御部１８によって制御される。送波ビームフォーマ６０４は、送波ビームフォーミング信号を送受切換ユニット６０６に入力する。 The transmission beamformer 604 performs beamforming of the transmission, and generates a beamforming signal for forming an ultrasonic beam having a predetermined direction based on the transmission timing signal. The beam forming signal is composed of a plurality of drive signals to which time differences corresponding to directions are given. Beam forming is controlled by the control unit 18 described later. The transmission beam former 604 inputs the transmission beam forming signal to the transmission / reception switching unit 606.

送受切換ユニット６０６は、ビームフォーミング信号を超音波トランスデューサアレイに入力する。超音波トランスデューサアレイにおいて、送波アパーチャ（ａｐｅｒｔｕｒｅ）を構成する複数の超音波トランスデューサは、駆動信号の時間差に対応した位相差を持つ超音波をそれぞれ発生する。それら超音波の波面合成により、所定方位の音線に沿った超音波ビームが形成される。 The transmission / reception switching unit 606 inputs a beamforming signal to the ultrasonic transducer array. In the ultrasonic transducer array, a plurality of ultrasonic transducers constituting a transmission aperture each generate an ultrasonic wave having a phase difference corresponding to a time difference of drive signals. An ultrasonic beam along a sound ray in a predetermined direction is formed by the wavefront synthesis of these ultrasonic waves.

送受切換ユニット６０６には受波ビームフォーマ６１０が接続されている。送受切換ユニット６０６は、超音波トランスデューサアレイ中の受波アパーチャが受波した複数のエコー信号を受波ビームフォーマ６１０に入力する。受波ビームフォーマ６１０は、送波の音線に対応した受波のビームフォーミングを行うもので、複数の受波エコーに時間差を付与して位相を調整し、次いでそれら加算して所定方位の音線に沿ったエコー受信信号を形成する。受波のビームフォーミングは後述の制御部１８により制御される。 A reception beam former 610 is connected to the transmission / reception switching unit 606. The transmission / reception switching unit 606 inputs a plurality of echo signals received by the reception aperture in the ultrasonic transducer array to the reception beam former 610. The receiving beam former 610 performs received beam forming corresponding to the sound ray of the transmitted wave, adjusts the phase by giving a time difference to a plurality of received echoes, and then adds them to generate sound in a predetermined direction. An echo reception signal along the line is formed. The beam forming of the received wave is controlled by the control unit 18 described later.

超音波ビームの送波は、送波タイミング発生ユニット６０２が発生する送波タイミング信号により、所定の時間間隔で繰り返し行われる。それに合わせて、送波ビームフォーマ６０４および受波ビームフォーマ６１０により、音線の方位が所定量ずつ変更される。それによって、対象４の内部が、音線によって順次に走査される。このような構成の送受信部６は、例えば図３に示すような走査を行う。すなわち、放射点２００からｚ方向に延びる音線２０２で扇状の２次元領域２０６をθ方向に走査し、いわゆるセクタスキャン（ｓｅｃｔｏｒｓｃａｎ）を行う。 Transmission of the ultrasonic beam is repeated at predetermined time intervals by a transmission timing signal generated by the transmission timing generation unit 602. In accordance with this, the direction of the sound ray is changed by a predetermined amount by the transmission beam former 604 and the reception beam former 610. Thereby, the inside of the object 4 is sequentially scanned by sound rays. The transceiver unit 6 having such a configuration performs scanning as shown in FIG. 3, for example. That is, the fan-shaped two-dimensional region 206 is scanned in the θ direction by the sound ray 202 extending in the z direction from the radiation point 200, and so-called sector scan is performed.

送波および受波のアパーチャを超音波トランスデューサアレイの一部を用いて形成するときは、このアパーチャをアレイに沿って順次移動させることにより、例えば図４に示すような走査を行うことができる。すなわち、放射点２００からｚ方向に発する音線２０２を直線状の軌跡２０４に沿って平行移動させることにより、矩形状の２次元領域２０６をｘ方向に走査し、いわゆるリニアスキャン（ｌｉｎｅａｒｓｃａｎ）を行う。 When the transmission and reception apertures are formed by using a part of the ultrasonic transducer array, the apertures are sequentially moved along the array to perform scanning as shown in FIG. 4, for example. That is, by moving a sound ray 202 emitted from the radiation point 200 in the z direction along a linear locus 204, a rectangular two-dimensional region 206 is scanned in the x direction, and a so-called linear scan is performed. Do.

なお、超音波トランスデューサアレイが、超音波送波方向に張り出した円弧に沿って形成されたいわゆるコンベックスアレイ（ｃｏｎｖｅｘａｒｒａｙ）である場合は、リニアスキャンと同様な音線走査により、例えば図５に示すように、音線２０２の放射点２００を円弧状の軌跡２０４に沿って移動させ、扇面状の２次元領域２０６をθ方向に走査して、いわゆるコンベックススキャンが行える。 When the ultrasonic transducer array is a so-called convex array formed along an arc extending in the ultrasonic wave transmission direction, for example, as shown in FIG. As described above, the so-called convex scan can be performed by moving the radiation point 200 of the sound ray 202 along the arc-shaped locus 204 and scanning the fan-shaped two-dimensional region 206 in the θ direction.

送受信部６はＢモード（ｍｏｄｅ）処理部１０およびドップラ処理部１２に接続されている。送受信部６から出力される音線ごとのエコー受信信号は、Ｂモード処理部１０およびドップラ処理部１２に入力される。 The transmission / reception unit 6 is connected to a B-mode processing unit 10 and a Doppler processing unit 12. The echo reception signal for each sound ray output from the transmission / reception unit 6 is input to the B-mode processing unit 10 and the Doppler processing unit 12.

Ｂモード処理部１０はＢモード画像を形成するものである。Ｂモード処理部１０は、図６に示すように、対数増幅ユニット１０２と包絡線検波ユニット１０４を備えている。Ｂモード処理部１０は、対数増幅ユニット１０２でエコー受信信号を対数増幅し、包絡線検波ユニット１０４で包絡線検波して音線上の個々の反射点でのエコーの強度を表す信号、すなわちＡスコープ（ｓｃｏｐｅ）信号を得て、このＡスコープ信号の各瞬時の振幅をそれぞれ輝度値として、Ｂモード画像を形成する。 The B-mode processing unit 10 forms a B-mode image. As shown in FIG. 6, the B mode processing unit 10 includes a logarithmic amplification unit 102 and an envelope detection unit 104. The B mode processing unit 10 logarithmically amplifies the echo reception signal by the logarithmic amplification unit 102, envelope detection by the envelope detection unit 104, and a signal representing the echo intensity at each reflection point on the sound ray, that is, an A scope A (scope) signal is obtained, and a B-mode image is formed using the instantaneous amplitude of the A scope signal as a luminance value.

ドップラ処理部１２はドップラ計測データを形成するものである。ドップラ計測データには、後述する速度データ、分散データおよびパワーデータが含まれる。ドップラ処理部１２は、図７に示すように直交検波ユニット１２０、ＭＴＩフィルタ（ｍｏｖｉｎｇｔａｒｇｅｔｉｎｄｉｃａｔｉｏｎｆｉｌｔｅｒ）１２２、自己相関演算ユニット１２４、平均流速演算ユニット１２６、分散演算ユニット１２８およびパワー（ｐｏｗｅｒ）演算ユニット１３０を備えている。 The Doppler processing unit 12 forms Doppler measurement data. The Doppler measurement data includes speed data, dispersion data, and power data, which will be described later. As shown in FIG. 7, the Doppler processing unit 12 includes a quadrature detection unit 120, an MTI filter (moving target indication filter) 122, an autocorrelation calculation unit 124, an average flow velocity calculation unit 126, a dispersion calculation unit 128, and a power calculation unit. 130 is provided.

ドップラ処理部１２は、直交検波ユニット１２０でエコー受信信号を直交検波し、ＭＴＩフィルタ１２２でＭＴＩ処理してエコー信号のドップラシフトを求める。また、自己相関演算ユニット１２４でＭＴＩフィルタ１２２の出力信号について自己相関演算を行い、平均流速演算ユニット１２６で自己相関演算結果から平均流速Ｖを求め、分散演算ユニット１２８で自己相関演算結果から流速の分散Ｔを求め、パワー演算ユニット１３０で自己相関演算結果からドップラ信号のパワーＰＷを求める。 The Doppler processing unit 12 performs quadrature detection of the echo reception signal by the quadrature detection unit 120 and performs MTI processing by the MTI filter 122 to obtain a Doppler shift of the echo signal. Further, the autocorrelation calculation unit 124 performs autocorrelation calculation on the output signal of the MTI filter 122, the average flow velocity calculation unit 126 obtains the average flow velocity V from the autocorrelation calculation result, and the dispersion calculation unit 128 calculates the flow velocity from the autocorrelation calculation result. The variance T is obtained, and the power calculation unit 130 obtains the power PW of the Doppler signal from the autocorrelation calculation result.

これによって、対象４内で移動するエコー源、例えば血液等の平均流速Ｖとその分散Ｔおよびドップラ信号のパワーＰＷを表すそれぞれのデータが音線ごとに得られる。これらデータは、音線上の各点の平均流速、分散およびパワーを示す。以下、平均流速を単に速度という。なお、速度は音線方向の成分として得られる。また、超音波プローブ２に近づく方向と遠ざかる方向とが区別される。なお、エコー源は血液に限るものではなく、例えば血管等に注入されたマイクロバルーン（ｍｉｃｒｏｂａｌｌｏｏｎ）造影剤等であって良い。以下、血液の例で説明するがマイクロバルーン造影剤の場合も同様である。 As a result, an echo source moving within the object 4, for example, an average flow velocity V of blood or the like, its dispersion T, and data representing the power PW of the Doppler signal are obtained for each sound ray. These data show the average flow velocity, dispersion and power at each point on the ray. Hereinafter, the average flow velocity is simply referred to as speed. The speed is obtained as a component in the sound ray direction. Further, a direction approaching the ultrasonic probe 2 is distinguished from a direction moving away. The echo source is not limited to blood, and may be, for example, a micro balloon contrast agent injected into a blood vessel or the like. Hereinafter, an example of blood will be described, but the same applies to the case of a microballoon contrast agent.

Ｂモード処理部１０およびドップラ処理部１２は画像処理部１４に接続されている。画像処理部１４は、Ｂモード処理部１０およびドップラ処理部１２からそれぞれ入力されるデータに基づいて、それぞれＢモード画像およびドップラ画像を生成する。以下、ドップラ画像をフロー画像ともいう。 The B mode processing unit 10 and the Doppler processing unit 12 are connected to the image processing unit 14. The image processing unit 14 generates a B-mode image and a Doppler image based on data input from the B-mode processing unit 10 and the Doppler processing unit 12, respectively. Hereinafter, the Doppler image is also referred to as a flow image.

画像処理部１４は、図８に示すように、バス（ｂｕｓ）１４０によって接続された入力データメモリ（ｄａｔａｍｅｍｏｒｙ）１４２、ディジタル・スキャンコンバータ（ｄｉｇｉｔａｌｓｃａｎｃｏｎｖｅｒｔｅｒ）１４４、画像メモリ１４６およびプロセッサ（ｐｒｏｃｅｓｓｏｒ）１４８を備えている。 As shown in FIG. 8, the image processing unit 14 includes an input data memory 142, a digital scan converter 144, an image memory 146, and a processor connected by a bus 140. 148.

Ｂモード処理部１０およびドップラ処理部１２から音線ごとに入力されたＢモード画像およびドップラ計測データは、入力データメモリ１４２にそれぞれ記憶される。入力データメモリ１４２のデータは、ディジタル・スキャンコンバータ１４４で走査変換されて画像メモリ１４６に記憶される。プロセッサ１４８は、入力データメモリ１４２および画像メモリ１４６のデータについてそれぞれ所定のデータ処理を施す。データ処理には空間フィルタリングが含まれる。空間フィルタリングについては後にあらためて説明する。 The B-mode image and Doppler measurement data input for each sound ray from the B-mode processing unit 10 and the Doppler processing unit 12 are stored in the input data memory 142, respectively. Data in the input data memory 142 is scan-converted by the digital scan converter 144 and stored in the image memory 146. The processor 148 performs predetermined data processing on the data in the input data memory 142 and the image memory 146, respectively. Data processing includes spatial filtering. Spatial filtering will be described later.

画像処理部１４には表示部１６が接続されている。表示部１６は、画像処理部１４から画像信号が与えられ、それに基づいて画像を表示するようになっている。なお、表示部１６は、カラー（ｃｏｌｏｒ）画像が表示可能なグラフィックディスプレー（ｇｒａｐｈｉｃｄｉｓｐｌａｙ）等で構成される。 A display unit 16 is connected to the image processing unit 14. The display unit 16 receives an image signal from the image processing unit 14 and displays an image based on the image signal. Note that the display unit 16 is configured by a graphic display or the like that can display a color image.

以上の送受信部６、Ｂモード処理部１０、ドップラ処理部１２、画像処理部１４および表示部１６には制御部１８が接続されている。制御部１８は、それら各部に制御信号を与えてその動作を制御する。また、被制御の各部から各種の報知信号が入力される。超音波プローブ２、送受信部６、Ｂモード処理部１０、ドップラ処理部１２、画像処理部１４および制御部１８からなる部分は、本発明における撮影手段の一例である。 A control unit 18 is connected to the transmission / reception unit 6, B-mode processing unit 10, Doppler processing unit 12, image processing unit 14, and display unit 16. The control unit 18 gives control signals to these units to control their operation. Also, various notification signals are input from each part to be controlled. A portion including the ultrasonic probe 2, the transmission / reception unit 6, the B-mode processing unit 10, the Doppler processing unit 12, the image processing unit 14, and the control unit 18 is an example of an imaging unit in the present invention.

制御部１８の制御の下で、Ｂモード動作およびドップラモード動作が実行される。制御部１８には操作部２０が接続されている。操作部２０は操作者によって操作され、制御部１８に適宜の指令や情報を入力するようになっている。操作部２０は、例えばキーボード（ｋｅｙｂｏａｒｄ）やポインティングデバイス（ｐｏｉｎｔｉｎｇｄｅｖｉｃｅ）およびその他の操作具を備えた操作パネル（ｐａｎｅｌ）で構成される。 Under the control of the control unit 18, the B mode operation and the Doppler mode operation are executed. An operation unit 20 is connected to the control unit 18. The operation unit 20 is operated by an operator and inputs appropriate commands and information to the control unit 18. The operation unit 20 includes, for example, an operation panel including a keyboard, a pointing device, and other operation tools.

本装置の動作を説明する。操作者は超音波プローブ２を対象４の所望の個所に当接し、操作部２０を操作して、例えばＢモードとドップラモードを併用した撮影動作を行う。これによって、制御部１８による制御の下で、Ｂモード撮影とドップラモード撮影が時分割で行われる。すなわち、例えばドップラモードのスキャンを所定回数行う度にＢモードのスキャンを１回行う割合で、Ｂモードとドップラモードの混合スキャンが行われる。 The operation of this apparatus will be described. The operator abuts the ultrasonic probe 2 on a desired portion of the object 4 and operates the operation unit 20 to perform an imaging operation using both the B mode and the Doppler mode, for example. As a result, under the control of the control unit 18, B-mode shooting and Doppler mode shooting are performed in a time-sharing manner. That is, for example, every time the Doppler mode scan is performed a predetermined number of times, the mixed scan of the B mode and the Doppler mode is performed at a rate of performing the B mode scan once.

Ｂモードにおいては、送受信部６は、超音波プローブ２を通じて音線順次で対象４の内部を走査して逐一そのエコーを受信する。Ｂモード処理部１０は、送受信部６から入力されるエコー受信信号を対数増幅ユニット１０２で対数増幅し包絡線検波ユニット１０４で包絡線検波してＡスコープ信号を求め、それに基づいて音線ごとのＢモード画像を形成する。画像処理部１４は、Ｂモード処理部１０から入力される音線ごとのＢモード画像を入力データメモリ１４２に記憶する。これによって、入力データメモリ１４２内に、Ｂモード画像についての音線データ空間が形成される。 In the B mode, the transmission / reception unit 6 scans the inside of the object 4 in the order of sound rays through the ultrasonic probe 2 and receives the echoes one by one. The B-mode processing unit 10 logarithmically amplifies the echo reception signal input from the transmission / reception unit 6 by the logarithmic amplification unit 102, detects the envelope by the envelope detection unit 104, obtains an A scope signal, and based on the A scope signal, A B-mode image is formed. The image processing unit 14 stores the B-mode image for each sound ray input from the B-mode processing unit 10 in the input data memory 142. As a result, a sound ray data space for the B-mode image is formed in the input data memory 142.

ドップラモードにおいては、送受信部６は超音波プローブ２を通じて音線順次で対象４の内部を走査して逐一そのエコーを受信する。その際、１音線当たり複数回の超音波の送波とエコーの受信が行われる。 In the Doppler mode, the transmission / reception unit 6 scans the inside of the object 4 in the order of sound rays through the ultrasonic probe 2 and receives the echoes one by one. At that time, ultrasonic waves are transmitted and echoes are received a plurality of times per sound ray.

ドップラ処理部１２は、エコー受信信号を直交検波ユニット１２０で直交検波し、ＭＴＩフィルタ１２２でＭＴＩ処理し、自己相関演算ユニット１２４で自己相関を求め、自己相関結果から、平均流速演算ユニット１２６で平均流速を求め、分散演算ユニット１２８で分散を求め、パワー演算ユニット１３０でパワーを求める。 The Doppler processing unit 12 performs quadrature detection on the echo reception signal by the quadrature detection unit 120, performs MTI processing by the MTI filter 122, obtains autocorrelation by the autocorrelation calculation unit 124, and calculates the average by the average flow velocity calculation unit 126 from the autocorrelation result. The flow rate is obtained, the dispersion is obtained by the dispersion operation unit 128, and the power is obtained by the power operation unit 130.

画像処理部１４は、ドップラ処理部１２から入力される音線ごとかつピクセルごとの各ドップラ計測データを入力データメモリ１４２に記憶する。これによって、入力データメモリ１４２内に、各ドップラ計測データについての音線データ空間が成される。 The image processing unit 14 stores the Doppler measurement data for each sound ray and each pixel input from the Doppler processing unit 12 in the input data memory 142. As a result, a sound ray data space for each Doppler measurement data is formed in the input data memory 142.

プロセッサ１４８は、入力データメモリ１４２のＢモード画像およびドップラ計測データをディジタル・スキャンコンバータ１４４でそれぞれ走査変換して画像メモリ１４６に書き込む。その際、ドップラ計測データは、速度に分散を加味したカラーフローマッピング（ＣＦＭ：ｃｏｌｏｒｆｌｏｗｍａｐｐｉｎｇ）画像およびパワードップラ（ＰＤＩ：ｐｏｗｅｒＤｏｐｐｌｅｒｉｍａｇｉｎｇ）画像として書き込まれる。 The processor 148 scan-converts the B-mode image and Doppler measurement data in the input data memory 142 by the digital scan converter 144 and writes them in the image memory 146. At this time, the Doppler measurement data is written as a color flow mapping (CFM) image and a power Doppler imaging (PDI) image in which dispersion is added to the speed.

プロセッサ１４８は、Ｂモード画像、ＣＦＭ画像およびＰＤＩ画像を別々な領域に書き込む。Ｂモード画像は、音線走査面の体内組織の断層像を示すものとなる。ＣＦＭ画像は、音線走査面の血流速度等の２次元分布を示す画像となる。この画像では血流の方向に応じて表示色を異ならせる。また、速度に応じて表示色の輝度を異ならせる。また、分散に応じて所定の色の混色率を高め表示色の純度を変える。これによって、速度の大きさ、方向、分散が視認性良く表示される。ＰＤＩ画像は、音線走査面の血流等の存在を示す画像となる。表示色はＣＦＭ画像に用いる色とは異ならせる。信号強度に応じて表示色の輝度を変える。ＣＦＭ画像およびＰＤＩ画像はいずれもフロー画像となる。 The processor 148 writes the B-mode image, CFM image, and PDI image in separate areas. The B-mode image shows a tomographic image of the body tissue on the sound ray scanning plane. The CFM image is an image showing a two-dimensional distribution such as a blood flow velocity on the sound ray scanning plane. In this image, the display color is varied depending on the direction of blood flow. Further, the brightness of the display color is varied according to the speed. In addition, the purity of the display color is changed by increasing the color mixture ratio of a predetermined color according to the dispersion. Thereby, the magnitude, direction, and dispersion of the speed are displayed with good visibility. The PDI image is an image indicating the presence of blood flow or the like on the sound ray scanning plane. The display color is different from the color used for the CFM image. The brightness of the display color is changed according to the signal intensity. Both the CFM image and the PDI image are flow images.

これらの画像を表示部１６に表示させる場合には、例えばＢモード画像とＣＦＭ画像を重ね合わせて表示する。これにより、体内組織との位置関係が明確な血流速度分布像を観察することができる。また、Ｂモード画像とＰＤＩ画像を重ね合わせて表示する。この場合は、体内組織との位置関係が明確な血管走行状態を観察することができる。 When these images are displayed on the display unit 16, for example, a B-mode image and a CFM image are superimposed and displayed. Thereby, a blood flow velocity distribution image with a clear positional relationship with the body tissue can be observed. In addition, the B-mode image and the PDI image are displayed in a superimposed manner. In this case, it is possible to observe a blood vessel running state in which the positional relationship with the body tissue is clear.

フロー画像について、プロセッサ１４８により、黒抜けを補正するための空間フィルタリングが行われる。空間フィルタリングに着目したプロセッサ１４８の機能ブロック図を図９に示す。同図は空間フィルタの構成を示す。本フィルタの構成によって、空間フィルタに関する本発明を実施するための最良の形態の一例が示される。本フィルタの動作によって、空間フィルタリング方法に関する本発明を実施するための最良の形態の一例が示される。 For the flow image, the processor 148 performs spatial filtering for correcting blackout. A functional block diagram of the processor 148 focusing on spatial filtering is shown in FIG. This figure shows the configuration of the spatial filter. An example of the best mode for carrying out the present invention relating to the spatial filter is shown by the configuration of the present filter. An example of the best mode for carrying out the present invention relating to a spatial filtering method is shown by the operation of the present filter.

同図に示すように、本フィルタは、黒抜け探索ユニット７０２、領域判定ユニット７０４および黒抜け補正ユニット７０６を有する。各ユニットはいずれもプロセッサ１４８の機能の一部である。黒抜け探索ユニット７０２は、入力画像の黒抜けを探索する。領域判定ユニット７０４は、黒抜け探索ユニット７０２の探索結果に基づいて、黒抜け部分が補正領域に該当するか否かを判定する。補正ユニット７０６は、領域判定ユニット７０４の判定結果に基づいて入力画像の黒抜けを補正し、補正後の画像を出力する。出力画像には必要に応じて適宜のスムージング（ｓｍｏｏｔｈｉｎｇ）を施すようにしてもよい。 As shown in the figure, this filter has a black spot search unit 702, an area determination unit 704, and a black spot correction unit 706. Each unit is part of the function of the processor 148. A black spot search unit 702 searches for black spots in the input image. The area determination unit 704 determines whether or not the black hole portion corresponds to the correction area based on the search result of the black hole search unit 702. The correction unit 706 corrects the blackout of the input image based on the determination result of the area determination unit 704, and outputs the corrected image. The output image may be appropriately smoothed as necessary.

黒抜け探索ユニット７０２は本発明における探索手段の一例である。領域判定ユニット７０４は本発明における判定手段の一例である。黒抜け補正ユニット７０６は本発明における補正手段の一例である。 The blackout search unit 702 is an example of search means in the present invention. The area determination unit 704 is an example of determination means in the present invention. The blackout correction unit 706 is an example of correction means in the present invention.

本フィルタの動作のフローチャート（ｆｌｏｗｃｈａｒｔ）を図１０に示す。同図に示すように、本フィルタは、ステージ（ｓｔａｇｅ）７１１で黒抜け探索を行う。黒抜け探索は黒抜け探索ユニット７０２によって行われる。黒抜け探索は、入力画像について、所定の条件に当てはまる画素値を持つ画素の位置（座標）を検出するとともに、そのような画素が連続する数（連続画素数）を求めることにより行われる。 A flow chart of the operation of this filter is shown in FIG. As shown in the figure, this filter performs a black spot search at a stage 711. The black spot search is performed by the black spot search unit 702. The blackout search is performed by detecting the position (coordinates) of a pixel having a pixel value that satisfies a predetermined condition in the input image and obtaining the number of such pixels that are continuous (number of continuous pixels).

所定の条件は閾値によって規定される。閾値を用いることにより条件設定を容易に行うことができる。閾値は単一の値あるいは所定の範囲の両端を規定する２つの閾値である。また、閾値は、画素値そのものについての閾値あるいは画素値の差についての閾値である。 The predetermined condition is defined by a threshold value. Condition setting can be easily performed by using the threshold value. The threshold value is a single value or two threshold values that define both ends of a predetermined range. The threshold value is a threshold value for the pixel value itself or a threshold value for a difference between pixel values.

画素値そのものについての閾値は、例えば、画素値のレンジ（ｒａｎｇｅ）の数％程度と設定される。画素値の差についての閾値は、例えば、減算前の値から見た差の値の減少量が画素値のレンジ（ｒａｎｇｅ）の数％程度と設定される。それらの閾値は可変とするのが、フィルタリングの効果を調節可能とする点で好ましい。 The threshold for the pixel value itself is set to, for example, about several percent of the range of the pixel value. For example, the threshold value for the difference between pixel values is set such that the amount of decrease in the difference value viewed from the value before subtraction is about several percent of the range of the pixel value. It is preferable that these threshold values be variable in that the effect of filtering can be adjusted.

画素値の差について規定された閾値を用いるときは、例えば図１１に示すように、注目画素の座標を（ｘ，ｙ）としたとき、座標（ｘ，ｙ）と座標（ｘ−１，ｙ）の画素値の差および座標（ｘ，ｙ）と座標（ｘ，ｙ−１）の画素値の差がいずれも所定の条件を満足するとき、座標（ｘ，ｙ）を黒抜け領域の始点とする。そして、ここを起点として、同じ条件の画素をｘ，ｙ両方向に探索し、そのような画素が連続して存在する位置および画素数を求める。このような黒抜け探索がフロー画像全体にわたって行われる。その結果、フロー画像における黒抜けの位置および大きさが個々に判明する。 When using a threshold defined for the difference in pixel values, for example, as shown in FIG. 11, when the coordinates of the pixel of interest are (x, y), the coordinates (x, y) and the coordinates (x-1, y ) And the difference between the pixel values of the coordinates (x, y) and the coordinates (x, y-1) satisfy the predetermined condition, the coordinates (x, y) are set as the start point of the blackout region. And Then, starting from here, a pixel having the same condition is searched in both the x and y directions, and the position and the number of pixels where such a pixel continuously exists are obtained. Such blackout search is performed over the entire flow image. As a result, the positions and sizes of black spots in the flow image are individually determined.

次に、ステージ７１３で、領域判定を行う。領域判定は領域判定ユニット７０４によって行われる。領域判定は、黒抜け探索によって求められた連続画素数に基づいて行われる。すなわち、連続画素数が所定限度以内であるときは、黒抜け補正を行うべき領域であると判定し、連続画素数が所定限度を超えるときは、黒抜け補正を行うべき領域でないと判定する。所定限度は、例えば、数十画素程度と設定される。所定限度は可変とするのが、フィルタリングを適切に行う点で好ましい。 Next, in stage 713, area determination is performed. The area determination is performed by the area determination unit 704. The area determination is performed based on the number of continuous pixels obtained by the blackout search. That is, when the number of continuous pixels is within a predetermined limit, it is determined that the region is to perform blackout correction, and when the number of continuous pixels exceeds the predetermined limit, it is determined that the region is not to be subjected to blackout correction. The predetermined limit is set to about several tens of pixels, for example. It is preferable that the predetermined limit is variable in terms of appropriate filtering.

なお、連続画素数が所定限度以内であっても、探索個所がフロー画像が存在する範囲を逸脱したときは、黒抜け補正を行うべき領域でないと判定する。入力画像におけるフロー画像の存在範囲は予め判明している。 Even if the number of continuous pixels is within a predetermined limit, if the search location deviates from the range where the flow image exists, it is determined that the region is not a region where blackout correction should be performed. The existence range of the flow image in the input image is known in advance.

連続画素数が所定限度を超えるということは、その領域が黒抜け部分（すなわち、本来流れが有るはずのところがたまたま黒く抜けたところ）ではなく、実際に流れが無い部分である可能性が高い。そこで、そのような領域は領域判定によって黒抜け補正の対象から除外される。これによって、黒抜け補正のための空間フィルタリングの妥当性を高めることができる。連続画素数が所定限度以内でかつ探索範囲がフロー画像の存在範囲を逸脱したときも同様である。 If the number of continuous pixels exceeds a predetermined limit, it is highly likely that the area is not a black-out portion (that is, a portion where the flow should have flowed by chance happens to be black) but a portion where there is actually no flow. Therefore, such a region is excluded from the target of blackout correction by region determination. As a result, the validity of spatial filtering for blackout correction can be increased. The same applies when the number of continuous pixels is within a predetermined limit and the search range deviates from the existence range of the flow image.

次に、ステージ７１５で、黒抜け補正を行う。黒抜け補正は、黒抜け補正ユニット７０６によって行われる。黒抜け補正ユニット７０６は、領域判定により補正対象と判定された領域について黒抜け補正を行う。黒抜け補正は、黒抜け領域内の画素値を、黒抜け領域の周囲の画素値から補間した画素値で置き換えることによって行う。補間は、単純補間、１次補間、２次補間等、適宜の補間であってよい。 Next, at stage 715, blackout correction is performed. Blackout correction is performed by a blackout correction unit 706. The blackout correction unit 706 performs blackout correction on the area determined as the correction target by the area determination. Blackout correction is performed by replacing pixel values in the blackout area with pixel values interpolated from the pixel values around the blackout area. The interpolation may be an appropriate interpolation such as simple interpolation, primary interpolation, or quadratic interpolation.

図１２に、空間フィルタリングに着目したプロセッサ１４８の他の機能ブロック図を示す。同図は空間フィルタの構成を示す。本フィルタの構成によって、空間フィルタに関する本発明を実施するための最良の形態の一例が示される。本フィルタの動作によって、空間フィルタリング方法に関する本発明を実施するための最良の形態の一例が示される。 FIG. 12 shows another functional block diagram of the processor 148 focusing on spatial filtering. This figure shows the configuration of the spatial filter. An example of the best mode for carrying out the present invention relating to the spatial filter is shown by the configuration of the present filter. An example of the best mode for carrying out the present invention relating to a spatial filtering method is shown by the operation of the present filter.

同図に示すように、本フィルタは、黒抜け探索ユニット８０２、領域判定ユニット８０４、黒抜け補正ユニット８０６および補正確定ユニット８０８を有する。各ユニットはいずれもプロセッサ１４８の機能の一部である。黒抜け探索ユニット８０２は、入力画像の黒抜けを探索する。領域判定ユニット８０４は、黒抜け探索ユニット８０２の探索結果に基づいて、黒抜け部分が補正領域に該当するか否かを判定する。補正ユニット８０６は、黒抜け探索ユニット８０２の探索結果に基づいて入力画像の黒抜けを補正する。補正確定ユニット８０８は、領域判定ユニット８０４の判定結果に基づいて、黒抜け補正ユニット８０６の補正を確定し、補正確定後の画像を出力する。出力画像には必要に応じて適宜のスムージングを施すようにしてもよい。 As shown in the figure, this filter has a black spot search unit 802, an area determination unit 804, a black spot correction unit 806, and a correction confirmation unit 808. Each unit is part of the function of the processor 148. A black spot search unit 802 searches for black spots in the input image. The area determination unit 804 determines whether or not the black missing portion corresponds to the correction area based on the search result of the black missing search unit 802. The correction unit 806 corrects blackout in the input image based on the search result of the blackout search unit 802. The correction confirmation unit 808 confirms the correction of the blackout correction unit 806 based on the determination result of the area determination unit 804, and outputs the image after the correction is determined. The output image may be appropriately smoothed as necessary.

黒抜け探索ユニット８０２は本発明における探索手段の一例である。領域判定ユニット８０４は本発明における判定手段の一例である。黒抜け補正ユニット８０６および補正確定ユニット８０８からなる部分は本発明における補正手段の一例である。 The blackout search unit 802 is an example of search means in the present invention. The area determination unit 804 is an example of determination means in the present invention. A portion including the blackout correction unit 806 and the correction confirmation unit 808 is an example of a correction unit in the present invention.

本フィルタの動作のフローチャートを図１３に示す。同図に示すように、本フィルタは、ステージ９０１で黒抜け探索を行う。黒抜け探索は黒抜け探索ユニット８０２によって行われる。黒抜け探索は、入力画像について、所定の条件に当てはまる画素値を持つ画素の位置（座標）を検出するとともに、そのような画素が連続する数（連続画素数）を求めることにより行われる。所定の条件は前述と同様である。探索の要領も前述と同等である。 FIG. 13 shows a flowchart of the operation of this filter. As shown in the figure, the filter performs a black spot search at stage 901. The blackout search is performed by a blackout search unit 802. The blackout search is performed by detecting the position (coordinates) of a pixel having a pixel value that satisfies a predetermined condition in the input image and obtaining the number of such pixels that are continuous (number of continuous pixels). The predetermined condition is the same as described above. The search procedure is the same as described above.

次に、ステージ９０３で、黒抜け補正を行う。黒抜け補正は、黒抜け補正ユニット８０６によって行われる。黒抜け補正ユニット８０６は、探索ユニット８０２によって探索された黒抜け部分について、探索され次第に黒抜け補正を行う。黒抜け補正の要領は前述と同様である。 Next, in the stage 903, blackout correction is performed. Blackout correction is performed by a blackout correction unit 806. The black missing correction unit 806 performs black missing correction as soon as the black missing portion searched by the search unit 802 is searched. The procedure for blackout correction is the same as described above.

次に、ステージ９０５で、領域判定を行う。領域判定は領域判定ユニット８０４によって行われる。領域判定は、黒抜け探索によって求められた連続画素数に基づいて行われる。すなわち、連続画素数が所定限度以内であるときは、黒抜け補正を行うべき領域であると判定し、連続画素数が所定限度を超えるときは、黒抜け補正を行うべき領域でないと判定する。所定限度は可変とするのが、フィルタリングの効果を調節可能とする点で好ましい。 Next, at stage 905, area determination is performed. The area determination is performed by the area determination unit 804. The area determination is performed based on the number of continuous pixels obtained by the blackout search. That is, when the number of continuous pixels is within a predetermined limit, it is determined that the region is to perform blackout correction, and when the number of continuous pixels exceeds the predetermined limit, it is determined that the region is not to be subjected to blackout correction. It is preferable that the predetermined limit is variable in that the effect of filtering can be adjusted.

連続画素数が所定限度を超えるということは、その領域が黒抜け部分ではなく、実際に流れが無い部分である可能性が高い。そこで、そのような領域は領域判定によって黒抜け補正の対象から除外される。これによって、黒抜け補正のための空間フィルタリングの妥当性を高めることができる。また、連続画素数が所定限度以内でかつ探索範囲がフロー画像の存在範囲を逸脱したときも同様である。 If the number of continuous pixels exceeds a predetermined limit, there is a high possibility that the area is not a black-out portion but a portion that does not actually flow. Therefore, such a region is excluded from the target of blackout correction by region determination. As a result, the validity of spatial filtering for blackout correction can be increased. The same applies when the number of continuous pixels is within a predetermined limit and the search range deviates from the existence range of the flow image.

次に、ステージ９０７で、補正確定を行う。補正確定は補正確定ユニット８０８によって行われる。補正確定は、領域判定により黒抜け補正を行うべきでないと判定された領域について、黒抜け補正を取り消すことによって行われる。補正の取消は、黒抜け探索時に記憶した画素位置および画素値を用いて、黒抜け補正を行うべきでないと判定された領域について元の画素値を復元することにより行われる。 Next, at stage 907, correction is confirmed. The correction confirmation is performed by a correction confirmation unit 808. The correction is determined by canceling the blackout correction for the area determined not to be subjected to the blackout correction by the area determination. The cancellation of the correction is performed by restoring the original pixel value for the area where it is determined that the blackout correction should not be performed, using the pixel position and the pixel value stored during the blackout search.

このように、黒抜けを、探索され次第に全てをいったん補正し、その後に領域判定に基づいて不要な補正を取り消す方法は、前述のように黒抜け補正を領域判定後に行うものよりもフィルタリング速度を高速化できる点で好ましい。 As described above, the method of correcting all black spots as soon as they are searched and then canceling unnecessary corrections based on the area determination has a higher filtering speed than that in which black area correction is performed after area determination as described above. It is preferable in that the speed can be increased.

以上のように、本発明ではフロー画像における黒抜け領域を適切に判定し、その領域の画素値を周辺の画素値から補間した画素値で置き換えるので、複数画素にわたる黒抜けを効果的に補正することができる。このような補正は、黒抜け領域の判定を行わない通常の空間フィルタリングのように、空間分解能やコントラスト（ｃｏｎｔｒａｓｔ）を低下させることがなく、また、画像輪郭の拡大を生じるということもない。 As described above, according to the present invention, the blackout region in the flow image is appropriately determined, and the pixel value in that region is replaced with the pixel value interpolated from the surrounding pixel values, so that the blackout over a plurality of pixels is effectively corrected. be able to. Such correction does not reduce the spatial resolution and contrast, and does not cause enlargement of the image contour, unlike normal spatial filtering that does not determine the blackout region.

超音波診断装置のブロック図である。It is a block diagram of an ultrasonic diagnostic apparatus. 送受信部のブロック図である。It is a block diagram of a transmission / reception part. 音線走査の概念を示す図である。It is a figure which shows the concept of sound ray scanning. 音線走査の概念を示す図である。It is a figure which shows the concept of sound ray scanning. 音線走査の概念を示す図である。It is a figure which shows the concept of sound ray scanning. Ｂモード処理部のブロック図である。It is a block diagram of a B mode processing unit. ドップラ処理部のブロック図である。It is a block diagram of a Doppler processing unit. 画像処理部のブロック図である。It is a block diagram of an image processing part. 空間フィルタのブロック図である。It is a block diagram of a spatial filter. 空間フィルタの動作のフローチャートである。It is a flowchart of operation | movement of a spatial filter. 画素の配列を示す図である。It is a figure which shows the arrangement | sequence of a pixel. 空間フィルタのブロック図である。It is a block diagram of a spatial filter. 空間フィルタの動作のフローチャートである。It is a flowchart of operation | movement of a spatial filter.

Explanation of symbols

２超音波プローブ
４撮影対象
６送受信部
１０Ｂモード処理部
１２ドップラ処理部
１４画像処理部
１６表示部
１８制御部
２０操作部
７０２，８０２黒抜け探索ユニット
７０４，８０４領域判定ユニット
７０６，８０６黒抜け補正ユニット
８０８補正確定ユニット

2 Ultrasonic probe 4 Imaging target 6 Transmission / reception unit 10 B-mode processing unit 12 Doppler processing unit 14 Image processing unit 16 Display unit 18 Control unit 20 Operation units 702 and 802 Blackout search units 704 and 804 Area determination units 706 and 806 Blackout Correction unit 808 Correction confirmation unit

Claims

Ultrasound diagnosis having imaging means for capturing a flow image in which pixels are displayed in black when a value relating to blood flow is small based on an echo of an ultrasound irradiated to a target, and a spatial filter for processing the flow image A device,
The spatial filter is
Search means for obtaining a position of a pixel having a pixel value that meets a condition of a blackout region where a plurality of black pixels are gathered in the flow image and a number of continuous pixels, wherein the coordinates of the pixel of interest in the flow image are (x, y), the difference between the pixel values of the coordinates (x, y) and the coordinates (x-1, y) and the difference of the pixel values of the coordinates (x, y) and the coordinates (x, y-1) are both A search means for searching in both the x and y directions using the pixel value difference condition, with the coordinates (x, y) as a starting point of the black area when the black area condition is satisfied;
In the flow image, when the number of continuous pixels is within a predetermined limit as a condition of the blackout region, a determination unit that determines a region where a pixel exists as a correction region;
Correction means for replacing the pixel value of the correction region with a pixel value interpolated from surrounding pixel values;
The ultrasonic diagnostic apparatus, wherein the spatial filter corrects the blackout region.

Ultrasound diagnosis having imaging means for capturing a flow image in which pixels are displayed in black when a value relating to blood flow is small based on an echo of an ultrasound irradiated to a target, and a spatial filter for processing the flow image A device,
The spatial filter is
Search means for obtaining a position of a pixel having a pixel value that meets a condition of a blackout region where a plurality of black pixels are gathered in the flow image and a number of continuous pixels, wherein the coordinates of the pixel of interest in the flow image are (x, y), the difference between the pixel values of the coordinates (x, y) and the coordinates (x-1, y) and the difference of the pixel values of the coordinates (x, y) and the coordinates (x, y-1) are both A search means for searching in both the x and y directions using the pixel value difference condition, with the coordinates (x, y) as a starting point of the black area when the black area condition is satisfied;
Correction means for replacing the pixel value of the pixel searched by the search means with a pixel value interpolated from surrounding pixel values;
When the number of continuous pixels in the flow image exceeds a predetermined limit as a condition of the blackout region, a determination unit that determines a region where the continuous pixel corrected by the correction unit exists as a restoration region;
Restoring means for restoring the pixel value of the restoration area to the pixel value before correction,
The ultrasonic diagnostic apparatus, wherein the spatial filter corrects the blackout region.

The ultrasonic diagnostic apparatus according to claim 1 or 2,
2. The ultrasonic diagnostic apparatus according to claim 1, wherein the condition of the blackout region is defined by a threshold value.

The ultrasonic diagnostic apparatus according to any one of claims 1 to 3,
The ultrasonic diagnostic apparatus according to claim 1, wherein the predetermined limit is variable.

In the ultrasonic diagnostic apparatus according to any one of claims 1 to 4,
The ultrasonic diagnostic apparatus, wherein the flow image is a color flow mapping image or a power Doppler image.

The ultrasonic diagnostic apparatus according to claim 5,
The ultrasonic diagnostic apparatus, wherein the flow image is displayed so as to overlap with a B-mode image.

The ultrasonic diagnostic apparatus according to any one of claims 1 to 6,
The ultrasonic diagnostic apparatus characterized in that the interpolation of the correction means is simple interpolation, primary interpolation or secondary interpolation.