JP4616614B2 - Ultrasonic diagnostic equipment - Google Patents

Ultrasonic diagnostic equipment Download PDF

Info

Publication number
JP4616614B2
JP4616614B2 JP2004306926A JP2004306926A JP4616614B2 JP 4616614 B2 JP4616614 B2 JP 4616614B2 JP 2004306926 A JP2004306926 A JP 2004306926A JP 2004306926 A JP2004306926 A JP 2004306926A JP 4616614 B2 JP4616614 B2 JP 4616614B2
Authority
JP
Japan
Prior art keywords
sample point
value
component
evaluation
blood flow
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP2004306926A
Other languages
Japanese (ja)
Other versions
JP2006116070A (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.)
Hitachi Ltd
Original Assignee
Aloka Co Ltd
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 Aloka Co Ltd filed Critical Aloka Co Ltd
Priority to JP2004306926A priority Critical patent/JP4616614B2/en
Publication of JP2006116070A publication Critical patent/JP2006116070A/en
Application granted granted Critical
Publication of JP4616614B2 publication Critical patent/JP4616614B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本発明は超音波診断装置に関し、特に速度信号に含まれる不要成分の除去に関する。   The present invention relates to an ultrasonic diagnostic apparatus, and more particularly to removal of unnecessary components contained in a velocity signal.

血流画像を形成する超音波診断装置においては、一般に、超音波の送受波により得られた受信信号が、直交検波処理回路を介して、MTIフィルタ(Moving Target Indicator Filter)などによって構成されるウォールモーションフィルタに入力される。そのフィルタにおいては、受信信号に含まれるクラッタ成分(心臓壁などの低域運動体成分)が除去される。フィルタを通過した信号は自己相関回路及び逆正接演算回路などを含む速度演算部に入力され、受信信号が速度信号に変換される。ここで、速度信号は各サンプル点ごとの速度値によって構成されるものである。速度信号はノイズ除去回路などを介して画像形成部へ出力される。ノイズ除去回路については後述する。画像形成部はデジタルスキャンコンバータなどを含み、速度信号に基づいてカラーの二次元血流画像を形成する。上記の自己相関回路は、一般に、同じビームアドレス上において複数の超音波ビームを形成して得られる複数のビームデータ間で自己相関演算を行う回路であり、具体的には、ビーム上の各サンプル点ごとに、時間的に隣接するビーム間で自己相関演算を行って速度値を求めるものである。   In an ultrasonic diagnostic apparatus that forms a blood flow image, generally, a received signal obtained by transmitting and receiving an ultrasonic wave is configured by an MTI filter (Moving Target Indicator Filter) or the like via an orthogonal detection processing circuit. Input to the motion filter. In the filter, clutter components (low-frequency moving body components such as a heart wall) included in the received signal are removed. The signal that has passed through the filter is input to a speed calculation unit including an autocorrelation circuit and an arctangent calculation circuit, and the received signal is converted into a speed signal. Here, the speed signal is constituted by a speed value for each sample point. The speed signal is output to the image forming unit via a noise removal circuit or the like. The noise removal circuit will be described later. The image forming unit includes a digital scan converter and the like, and forms a color two-dimensional blood flow image based on the speed signal. The above autocorrelation circuit is generally a circuit that performs autocorrelation calculation between a plurality of beam data obtained by forming a plurality of ultrasonic beams on the same beam address, specifically, each sample on the beam. For each point, a speed value is obtained by performing autocorrelation calculation between temporally adjacent beams.

上記の各サンプル点の速度値からなる速度信号には、ウォールモーションフィルタで除去できなかった残留クラッタ成分が含まれ、また、ノイズ除去回路で除去できなかった残留ノイズ成分も含まれる。残留クラッタ成分及び残留ノイズ成分は残留不要成分(広義の残留ノイズ成分)として定義される。そのような残留不要成分は、血流画像の画質を低下させる要因となるため、その効果的な除去が望まれる。   The velocity signal composed of the velocity values of the respective sample points includes a residual clutter component that cannot be removed by the wall motion filter, and also includes a residual noise component that cannot be removed by the noise removal circuit. The residual clutter component and the residual noise component are defined as a residual unnecessary component (residual noise component in a broad sense). Such residual unnecessary components cause a decrease in the image quality of the blood flow image, and therefore effective removal thereof is desired.

下記の特許文献1には、MTIフィルタ、相関演算器、セレクタなどを有する超音波診断装置が開示されている。当該特許文献1には、MTIフィルタの出力信号のレベルに基づいて、あるいは、MTIフィルタの入力信号と出力信号の差に基づいて、自己相関演算により得られた速度信号を選択的に通過させるセレクタの動作を制御することが記載されている。下記の特許文献2には、自己相関演算器の出力における実数部成分及び虚数部成分の組み合わせに基づいて、不要信号(ノイズ成分及びクラッタ成分)を除去する超音波診断装置が記載されている。下記の特許文献3には、パワー、流速値、及び、分散の組み合わせに基づいて、ノイズ成分とクラッタ成分とを弁別、除去する超音波診断装置が記載されている。しかし、特許文献1乃至3には、空間的な画素値のランダム性(バラツキ)を考慮すること、あるいは、空間的な相関演算については記載されていない。   Patent Document 1 below discloses an ultrasonic diagnostic apparatus having an MTI filter, a correlation calculator, a selector, and the like. The Patent Document 1 discloses a selector that selectively passes a speed signal obtained by autocorrelation based on the level of the output signal of the MTI filter or based on the difference between the input signal and the output signal of the MTI filter. It is described that the operation is controlled. Patent Document 2 below describes an ultrasonic diagnostic apparatus that removes unnecessary signals (noise component and clutter component) based on a combination of a real part component and an imaginary part component in the output of an autocorrelation calculator. Patent Document 3 below describes an ultrasonic diagnostic apparatus that discriminates and removes noise components and clutter components based on a combination of power, flow velocity value, and dispersion. However, Patent Documents 1 to 3 do not describe consideration of spatial pixel value randomness (variation) or spatial correlation calculation.

下記の特許文献4には、二値化処理後の窓内の画素値の合計及び画素パターンに基づいて、窓内の注目画素についての処理内容を判定する超音波診断装置が記載されている。下記の特許文献5には、孤立画素あるいは少数個連続した画素をノイズとみなして除去する超音波診断装置が開示されている。特許文献4及び5には、上記特許文献1乃至3と同様に、空間的な画素値のランダム性を考慮すること、あるいは、空間的な相関演算については記載されていない。   Patent Document 4 below describes an ultrasonic diagnostic apparatus that determines the processing content of a pixel of interest in a window based on the sum of pixel values in the window after binarization processing and the pixel pattern. Patent Document 5 below discloses an ultrasonic diagnostic apparatus that removes isolated pixels or a small number of consecutive pixels as noise. Patent Documents 4 and 5 do not describe consideration of spatial pixel value randomness or spatial correlation calculation, similar to Patent Documents 1 to 3.

特開2003−250802号公報JP 2003-250802 A 特開平7−178092号公報Japanese Patent Laid-Open No. 7-178092 特開平4−218143号公報JP-A-4-218143 特開平7−108043号公報JP 7-108043 A 特開平8−280683号公報JP-A-8-280683

血流画像の画質を高めるためには、クラッタ成分やノイズ成分、つまり不要成分を効果的に除去することが求められ、そのために、クラッタ成分除去手段及びノイズ成分除去手段が設けられるが、不要成分は血流成分と同じような速度値、パワーなどを有する場合もあり、従来方法では必ずしも効果的に不要信号を除去できない場合がある。特に、クラッタ成分及びノイズ成分の除去処理後に残留する残留不要成分を更に除去することが望まれるが、上記各文献にはそのような後処理を行うための構成は開示されていない。   In order to improve the image quality of a blood flow image, it is required to effectively remove clutter components and noise components, that is, unnecessary components. For this purpose, clutter component removing means and noise component removing means are provided. May have the same velocity value and power as the blood flow component, and the conventional method may not always effectively remove unnecessary signals. In particular, it is desirable to further remove residual unnecessary components that remain after the removal processing of the clutter component and the noise component, but the above-mentioned documents do not disclose a configuration for performing such post-processing.

本発明の目的は、速度信号中に含まれる不要成分を効果的に除去することにある。   An object of the present invention is to effectively remove unnecessary components contained in a speed signal.

本発明の他の目的は、クラッタ成分又はノイズ成分の除去処理後においてなお残留する残留不要成分を除去することにある。   Another object of the present invention is to remove residual unnecessary components that remain after the removal processing of clutter components or noise components.

本発明の他の目的は、不要成分又は残留不要成分をその空間的な性質を利用して効果的に除去することにある。   Another object of the present invention is to effectively remove unnecessary components or residual unnecessary components by utilizing their spatial properties.

(1)本発明は、超音波を送受波して得られた受信信号に基づいて、各サンプル点の速度値を演算する速度演算手段と、前記各サンプル点の速度値を三値化処理し、前記各サンプル点ごとに属性値を求める三値化手段と、前記三値化処理後の空間的な属性値配列に基づいて、前記各サンプル点を基準として設定される注目領域について空間的なランダム性を評価するためのランダム性評価演算を実行し、前記各サンプル点ごとに評価値を求める評価手段と、前記各サンプル点ごとの評価値に基づいて、前記各サンプル点の速度値が血流成分であるか不要成分であるかを判定し、不要成分を除去する処理を実行する除去手段と、前記不要信号を除去する処理を実行した後に、血流成分に基づいて血流画像を形成する血流画像形成手段と、を含むことを特徴とする。 (1) In the present invention, based on a received signal obtained by transmitting and receiving an ultrasonic wave, speed calculation means for calculating a speed value at each sample point, and a ternary process for the speed value at each sample point. , A ternary means for obtaining an attribute value for each sample point, and a spatial region for a region of interest set on the basis of each sample point based on the spatial attribute value array after the ternarization process Randomness evaluation calculation for evaluating randomness is performed, evaluation means for obtaining an evaluation value for each sample point, and the velocity value of each sample point is calculated based on the evaluation value for each sample point. After determining whether it is a flow component or an unnecessary component and executing a process for removing the unnecessary component and executing the process for removing the unnecessary signal, a blood flow image is formed based on the blood flow component Blood flow image forming means And wherein the door.

上記構成によれば、超音波の送受波により得られた受信信号に基づいて各サンプル点の速度値が演算され、その速度値が三値化処理によりいずれかの属性値に規格化される。そして、空間的な属性値配列(属性値集合)における属性値のランダム性が評価値として評価され、その評価値に基づいて不要成分(又は残留不要成分)が判別、除去される。つまり、比較的大きなランダム性をもった信号成分を不要成分として認識し、それが血流画像形成に先立って取り除かれる。   According to the above configuration, the velocity value of each sample point is calculated based on the received signal obtained by the transmission / reception of the ultrasonic wave, and the velocity value is normalized to any attribute value by the ternary processing. Then, the randomness of the attribute values in the spatial attribute value array (attribute value set) is evaluated as an evaluation value, and unnecessary components (or residual unnecessary components) are discriminated and removed based on the evaluation values. That is, a signal component having a relatively large randomness is recognized as an unnecessary component and is removed prior to blood flow image formation.

一般に、血流部については、それを空間的に見ると、いわゆる折り返し部を除いて、個々の速度値の大きさや符号がある程度連続する傾向にあり、その一方、不要成分は空間的に散在し、その大きさや符号はランダム性を呈する。特に、クラッタ成分除去後且つノイズ成分除去後の速度信号については、そのことが顕著に認められる。その場合、単に注目するサンプル点の速度値を注目するだけでは、そのような空間的な傾向を認識できないが、本発明によれば、注目するサンプル点を基準として注目領域について空間的なランダム性を評価し、その評価結果を不要成分(又は残留不要成分)の判別に利用することができる。   Generally, regarding the blood flow part, when viewed spatially, the magnitude and sign of individual velocity values tend to be somewhat continuous except for the so-called folded part, while unnecessary components are spatially scattered. The size and code exhibit randomness. This is particularly noticeable for the speed signal after the clutter component is removed and after the noise component is removed. In such a case, it is not possible to recognize such a spatial tendency simply by focusing on the velocity value of the sample point of interest. However, according to the present invention, spatial randomness of the region of interest based on the sample point of interest is used. The evaluation result can be used for discrimination of unnecessary components (or residual unnecessary components).

上記の三値化処理から除去処理までの一連の処理は、クラッタ成分除去処理及びノイズ成分除去処理がなされた後に行われるのが特に望ましい。つまり、残留不要成分に対する処理であるのが望ましい。ただし、クラッタ成分除去処理及びノイズ成分除去処理の一方又は両方がなされない段階において適用しても、あるいは、それらの処理と並行して適用しても、一定の効果が得られる。クラッタ成分除去処理は公知のMTIフィルタを用いて低速運動体成分を除去する処理であるのが望ましい。ノイズ成分除去処理としては各種の公知方法をあげることができ、例えば、サンプル点のパワー値に対する閾値処理を含む処理などがある。なお、ノイズ成分除去処理の中に上記の一連の処理を組み込んでもよい。   It is particularly desirable that the series of processes from the ternarization process to the removal process be performed after the clutter component removal process and the noise component removal process are performed. That is, it is desirable that the process be performed on residual unnecessary components. However, even if it is applied at a stage where one or both of the clutter component removal processing and the noise component removal processing are not performed, or is applied in parallel with these processing, a certain effect can be obtained. The clutter component removal processing is desirably processing for removing low-speed moving body components using a known MTI filter. Various known methods can be used as the noise component removal process, and examples include a process including a threshold process for the power value of the sample point. Note that the series of processes described above may be incorporated in the noise component removal process.

上記構成においては、空間的なランダム性の評価に先立って三値化処理が施されており、すなわち速度値が規格化されているので、不要信号の空間的性質に着目して、空間的なランダム性の評価演算を簡略化でき、また評価精度を高められるという利点がある。変形例としては、三値化処理に代替する二値化処理や四値以上の規格化処理をあげることができるが、不要成分を空間的に判別する際には三値化処理が望ましい。上記の注目領域は、一次元領域あるいは二次元領域であるのが望ましいが、三次元領域であってもよい。   In the above configuration, the ternarization processing is performed prior to the evaluation of the spatial randomness, that is, the speed value is standardized. There are advantages that the randomness evaluation calculation can be simplified and the evaluation accuracy can be improved. As a modification, a binarization process that substitutes for the ternarization process and a standardization process of four or more values can be given, but the ternarization process is desirable when spatially discriminating unnecessary components. The region of interest is preferably a one-dimensional region or a two-dimensional region, but may be a three-dimensional region.

望ましくは、前記三値化処理では、前記各サンプル点の速度値が第1条件を満たす場合に正属性値が付与され、前記各サンプル点の速度値が第2条件を満たす場合にゼロ属性値が付与され、前記各サンプル点の速度値が第3条件を満たす場合に負属性値が付与される。すなわち、速度軸が3つの区間に区分され、サンプル点の速度値が属する区分が判定される。ここで、第2条件は、速度値がゼロであるという条件であってもよいし、速度値ゼロを挟んだ正負極性にまたがる一定区間に速度値が入るという条件であってもよい。上記の一連の処理の前段でクラッタ成分除去処理及びノイズ成分除去処理を行う場合には不要信号がかなり除去されている筈であり、このため第2条件として、速度値がゼロであるという条件を採用することができる。この三値化処理は、空間的なランダム性を評価するための前処理に相当するものである。   Preferably, in the ternarization process, a positive attribute value is given when the velocity value of each sample point satisfies the first condition, and a zero attribute value when the velocity value of each sample point satisfies the second condition. And a negative attribute value is assigned when the velocity value of each sample point satisfies the third condition. That is, the velocity axis is divided into three sections, and the division to which the velocity value of the sample point belongs is determined. Here, the second condition may be a condition that the speed value is zero, or may be a condition that the speed value enters a certain section that spans the positive and negative polarity across the speed value zero. When the clutter component removal process and the noise component removal process are performed in the preceding stage of the above-described series of processes, the unnecessary signal should be considerably removed. Therefore, the second condition is that the speed value is zero. Can be adopted. This ternary processing corresponds to preprocessing for evaluating spatial randomness.

望ましくは、前記評価手段は、前記注目領域とそれを空間的にシフトさせた参照領域との間で、前記ランダム性評価演算を実行する。注目領域と参照領域の形態(大きさ及び形状)は一致させておくのが望ましく、両者間のシフト方向及びシフト量は適宜設定することができる。例えば、シフト量はサンプル点1個であってもよい。注目領域と参照領域とが互いに異なるビームあるいはフレームに存在してもよい。望ましくは、前記ランダム性評価演算は相関演算に相当し、前記評価値は相関値に相当する。   Desirably, the evaluation means performs the randomness evaluation calculation between the region of interest and a reference region obtained by spatially shifting the region of interest. It is desirable that the forms (size and shape) of the attention area and the reference area coincide with each other, and the shift direction and the shift amount between them can be appropriately set. For example, the shift amount may be one sample point. The attention area and the reference area may exist in different beams or frames. Preferably, the randomness evaluation calculation corresponds to a correlation calculation, and the evaluation value corresponds to a correlation value.

望ましくは、前記注目領域はビーム方向に沿った一次元領域であり、前記参照領域は前記注目領域をビーム方向にシフトさせた領域である。一次元領域の場合には、1つのビームに対応する1つの受信信号(時系列信号)を用いてビーム上の各サンプル点についてランダム性の評価演算を行えるので簡便である。望ましくは、前記注目領域はフレーム上の二次元領域であり、前記参照領域は前記注目領域をフレーム上でシフトさせた領域である。シフト方向は水平方向及び垂直方向の一方又は両方である。   Preferably, the region of interest is a one-dimensional region along the beam direction, and the reference region is a region obtained by shifting the region of interest in the beam direction. In the case of a one-dimensional region, it is convenient because a randomness evaluation calculation can be performed for each sample point on the beam using one received signal (time-series signal) corresponding to one beam. Preferably, the attention area is a two-dimensional area on the frame, and the reference area is an area obtained by shifting the attention area on the frame. The shift direction is one or both of the horizontal direction and the vertical direction.

望ましくは、前記除去手段は、前記各サンプル点ごとに、それが有する速度値とそれについて演算された評価値との組み合わせに基づいて、前記各サンプル点の速度値が血流成分であるか不要成分であるかを判定する。この構成によれば、血流成分と不要成分との間における速度値の差及び評価値の差を利用して、両者を高精度に弁別できる。更に、第3の参照情報を参酌してもよい。望ましくは、前記除去手段は、前記速度値が第1の判定条件を満たし、且つ、前記評価値が第2の判定条件を満たす場合に不要成分であると判定する。   Preferably, for each of the sample points, the removing means is based on a combination of the velocity value of the sample point and the evaluation value calculated for the sample point, or whether the velocity value of each sample point is a blood flow component or not. Determine if it is a component. According to this configuration, the difference between the velocity value and the evaluation value between the blood flow component and the unnecessary component can be used to distinguish the two with high accuracy. Further, the third reference information may be taken into consideration. Preferably, the removing unit determines that the velocity value is an unnecessary component when the velocity value satisfies a first determination condition and the evaluation value satisfies a second determination condition.

(2)本発明は、超音波を送受波して得られた受信信号に含まれるクラッタ成分を除去するクラッタ成分除去手段と、前記クラッタ成分の除去後における受信信号に基づいて、各サンプル点の速度値を演算する速度演算手段と、前記各サンプル点ごとの速度値によって構成される速度信号に含まれるノイズ成分を除去するためのノイズ成分除去手段と、前記各サンプル点の速度値に基づいて、前記各サンプル点ごとに属性値を求める属性値演算手段と、前記属性値演算後の空間的な属性値配列に基づいて、前記各サンプル点を基準として設定される注目領域について空間的なランダム性を評価するためのランダム性評価演算を実行し、前記各サンプル点ごとに評価値を求める評価手段と、前記各サンプル点ごとの速度値及び評価値に基づいて、前記各サンプル点の速度値が血流成分であるか残留不要成分であるかを判定し、残留不要成分を除去する処理を実行する残留不要成分除去手段と、前記残留不要成分を除去する処理を実行した後に、血流成分に基づいて血流画像を形成する血流画像形成手段と、を含むことを特徴とする。 (2) The present invention provides a clutter component removing means for removing a clutter component included in a received signal obtained by transmitting / receiving ultrasonic waves, and the received signal after the removal of the clutter component based on the received signal of each sample point. Based on speed calculation means for calculating a speed value, noise component removal means for removing a noise component included in a speed signal constituted by the speed value for each sample point, and the speed value of each sample point , An attribute value calculation means for obtaining an attribute value for each sample point, and a spatial random for a region of interest set based on each sample point based on a spatial attribute value array after the attribute value calculation Based on the evaluation means for calculating the evaluation value for each sample point, the speed value and the evaluation value for each sample point Determining whether the velocity value of each sample point is a blood flow component or a residual unnecessary component, and executing a process of removing the residual unnecessary component; and a process of removing the residual unnecessary component Blood flow image forming means for forming a blood flow image based on a blood flow component after execution.

上記構成によれば、まず、クラッタ成分除去手段によって受信信号中のクラッタ成分に対する除去処理がなされ、次に、ノイズ成分除去手段によってホワイトノイズなどのノイズ成分に対する除去処理がなされる(その場合、残留クラッタも除去処理の対象となる)。加えて、それらの処理によっても残存してしまう残留不要成分(残留ノイズと残留クラッタの混在体とみなせる)を効果的に除去できる。   According to the above configuration, the clutter component removing unit firstly removes the clutter component in the received signal, and then the noise component removing unit removes the noise component such as white noise (in this case, the residual component). Clutter is also subject to removal). In addition, it is possible to effectively remove residual unnecessary components (which can be regarded as a mixture of residual noise and residual clutter) that remain even after such treatment.

個々のサンプル点に着目すると、残留不要成分は血流成分と同じような速度値及び符号を有する場合も多く、速度値のみを参酌するだけでは、血流成分と残留不要成分とを弁別するのが困難となるが、空間的にある程度の広がりをもって見れば、血流成分と残留不要成分との間にランダム性の評価値において違いが生じる。そこで、それを血流成分と残留不要成分との弁別に利用できる。ここで、評価値のみを基準として弁別を行うことも場合によっては可能であるが、弁別精度を確保するために、評価値に加えて速度値自体を参酌するのが望ましい。もちろん、それ以外の第3の情報を参酌することもできる。   Focusing on the individual sample points, the residual unnecessary component often has the same velocity value and sign as the blood flow component, and the blood flow component and the residual unnecessary component are discriminated only by considering only the velocity value. However, if it is viewed with a certain spatial extent, there is a difference in the evaluation value of randomness between the blood flow component and the residual unnecessary component. Therefore, it can be used for discrimination between blood flow components and residual unnecessary components. Here, although it is possible in some cases to discriminate on the basis of only the evaluation value, it is desirable to consider the speed value itself in addition to the evaluation value in order to ensure discrimination accuracy. Of course, other third information can be taken into consideration.

望ましくは、前記属性値演算手段は、前記各サンプル点の速度値が第1条件を満たす場合に正属性値を付与し、前記各サンプル点の速度値が第2条件を満たす場合にゼロ属性値を付与し、前記各サンプル点の速度値が第3条件を満たす場合に負属性値を付与する三値化処理手段である。望ましくは、前記評価手段は、前記注目領域とそれを空間的にシフトさせた参照領域との間で、相関演算に相当する前記ランダム性評価演算を実行する。   Preferably, the attribute value calculation means gives a positive attribute value when the velocity value of each sample point satisfies a first condition, and a zero attribute value when the velocity value of each sample point satisfies a second condition. And ternarization processing means for assigning a negative attribute value when the velocity value of each sample point satisfies the third condition. Preferably, the evaluation means executes the randomness evaluation calculation corresponding to a correlation calculation between the attention area and a reference area obtained by spatially shifting the attention area.

以上説明したように、本発明によれば、血流画像の画質を向上できる。特に、速度信号中に含まれる不要成分あるいは残留不要信号を効果的に除去できる。   As described above, according to the present invention, the image quality of a blood flow image can be improved. In particular, unnecessary components or residual unnecessary signals included in the speed signal can be effectively removed.

以下、本発明の好適な実施形態を図面に基づいて説明する。   DESCRIPTION OF EXEMPLARY EMBODIMENTS Hereinafter, preferred embodiments of the invention will be described with reference to the drawings.

図1には、本発明に係る超音波診断装置の好適な実施形態が示されており、図1はその全体構成を示すブロック図である。この超音波診断装置は超音波の送受波により得られた受信信号に基づいてドプラ情報を抽出し、それに基づいて二次元(あるいは三次元)の血流画像を形成する機能を備えた装置である。   FIG. 1 shows a preferred embodiment of an ultrasonic diagnostic apparatus according to the present invention, and FIG. 1 is a block diagram showing the overall configuration thereof. This ultrasonic diagnostic apparatus is an apparatus having a function of extracting Doppler information based on a received signal obtained by transmission / reception of an ultrasonic wave and forming a two-dimensional (or three-dimensional) blood flow image based on the extracted signal. .

プローブ10は超音波の送受波を行う送受波器であり、本実施形態においてはプローブ10内に複数の振動素子からなるアレイ振動子が設けられている。アレイ振動子によって超音波ビームが形成され、その超音波ビームが電子走査されて走査面が形成される。電子走査方式としては電子リニア走査、電子セクタ走査などをあげることができる。血流画像の形成のために、一般に各ビームアドレスごとに複数回の超音波パルスの送信が行われる。プローブ10上にいわゆる2Dアレイ振動子を設け、これによって三次元データ取込領域を形成するようにしてもよい。プローブ10は一般に体表面上に当接して用いられるものであるが、それが体腔内に挿入されるものであってもよい。   The probe 10 is a transducer that transmits and receives ultrasonic waves. In this embodiment, an array transducer including a plurality of transducer elements is provided in the probe 10. An ultrasonic beam is formed by the array transducer, and the ultrasonic beam is electronically scanned to form a scanning surface. Examples of the electronic scanning method include electronic linear scanning and electronic sector scanning. In order to form a blood flow image, transmission of ultrasonic pulses is generally performed a plurality of times for each beam address. A so-called 2D array transducer may be provided on the probe 10 to form a three-dimensional data capturing area. The probe 10 is generally used in contact with the body surface, but it may be inserted into a body cavity.

送受信部12は送信ビームフォーマー及び受信ビームフォーマーとして機能する。すなわち送受信部12は、送信時において、複数の送信信号を所定の遅延関係をもって複数の振動素子に供給する。これによって送信ビームが形成される。一方、体内からの反射波は複数の振動素子にて受波され、それによって得られた複数の受信信号が送受信部12に入力される。送受信部12はそれらの受信信号に対して整相加算処理を適用し、これによって整相加算後の受信信号を得る。すなわちこのような整相加算処理により電子的に受信ビームが形成される。送受信部12は受信信号をサンプリングするA/D変換部を有する。   The transmission / reception unit 12 functions as a transmission beam former and a reception beam former. That is, the transmission / reception unit 12 supplies a plurality of transmission signals to the plurality of vibration elements with a predetermined delay relationship at the time of transmission. As a result, a transmission beam is formed. On the other hand, reflected waves from the body are received by a plurality of vibration elements, and a plurality of reception signals obtained thereby are input to the transmission / reception unit 12. The transmission / reception unit 12 applies a phasing addition process to these reception signals, thereby obtaining a reception signal after phasing addition. That is, a reception beam is electronically formed by such a phasing addition process. The transmission / reception unit 12 includes an A / D conversion unit that samples a received signal.

直交検波部14は、入力される受信信号を複素信号に変換する機能を有する。本実施形態において、直交検波部14は受信信号に対して参照信号を混合する一対のミキサ回路を備えている。もちろん、受信信号を複素信号に変換できる限りにおいて各種の回路構成を採用することができる。直交検波後の受信信号(複素信号)はMTIフィルタ部16に入力される。このMTIフィルタ部16は実数部用のMTIフィルタ及び虚数部用のMTIフィルタによって構成され、受信信号に含まれる低速運動体成分(クラッタ)を除去する。すなわち、このMTIフィルタ部16は複素ウォールモーションフィルタとして機能するものである。ちなみに、血流ではなく心臓壁などの運動組織を画像化する場合にはこのMTIフィルタ16を経由しない受信信号が後段の速度演算部18へ入力される。   The quadrature detection unit 14 has a function of converting an input received signal into a complex signal. In the present embodiment, the quadrature detection unit 14 includes a pair of mixer circuits that mix a reference signal with a received signal. Of course, various circuit configurations can be employed as long as the received signal can be converted into a complex signal. The received signal (complex signal) after quadrature detection is input to the MTI filter unit 16. The MTI filter unit 16 includes an MTI filter for a real part and an MTI filter for an imaginary part, and removes a low-speed moving body component (clutter) included in the received signal. That is, the MTI filter unit 16 functions as a complex wall motion filter. Incidentally, when imaging not a blood flow but a moving tissue such as a heart wall, a received signal that does not pass through the MTI filter 16 is input to the speed calculation unit 18 at the subsequent stage.

速度演算部18は自己相関演算器及び速度演算器を有する。自己相関演算器は公知のように同じビームアドレス上において得られた時間的に隣接する複数の受信信号間において自己相関演算を実行する。それにより得られた自己相関結果は速度演算器に入力され、その演算結果から血流についての速度値(平均流速)が得られる。具体的には各サンプル点ごとに速度値が得られ、それらは速度信号を構成するものである。速度演算部18が速度値の演算機能の他、パワーを演算する機能、分散値(各サンプル点ごとの時間軸上の速度分散値)を演算する機能などを具備していてもよい。   The speed calculation unit 18 includes an autocorrelation calculator and a speed calculator. As is well known, the autocorrelation calculator performs an autocorrelation calculation between a plurality of temporally adjacent received signals obtained on the same beam address. The autocorrelation result obtained thereby is input to a velocity calculator, and a velocity value (average flow velocity) for the blood flow is obtained from the calculation result. Specifically, velocity values are obtained for each sample point, and they constitute a velocity signal. In addition to the speed value calculation function, the speed calculation unit 18 may include a function for calculating power, a function for calculating a variance value (speed variance value on the time axis for each sample point), and the like.

ノイズ処理部20は、速度信号中に含まれるホワイトノイズ上のノイズを除去するための回路である。MTIフィルタ部16によって除去仕切れなかった残留クラッタ成分についてもノイズ処理部20によってその除去が図られている。ノイズ処理部20は例えばパワー値と一定の閾値とを比較してノイズレベルの信号を除去する回路として構成され、あるいは速度値、パワー値、分散値などの複数の情報の組合せにしたがって不要信号を除去する回路として構成される。例えば、上記の特許文献1〜特許文献5などに掲げた回路を利用してもよい。   The noise processing unit 20 is a circuit for removing noise on the white noise included in the speed signal. The residual clutter component that has not been removed by the MTI filter unit 16 is also removed by the noise processing unit 20. The noise processing unit 20 is configured, for example, as a circuit that removes a noise level signal by comparing a power value with a certain threshold value, or an unnecessary signal according to a combination of a plurality of information such as a speed value, a power value, and a variance value. It is configured as a circuit to be removed. For example, you may utilize the circuit hung up in said patent document 1-patent document 5.

本来的には、上記のMTIフィルタ部16及びノイズ処理部20によって受信信号中に含まれるクラッタ成分やノイズ成分はすべて除去されるはずであるが、実際にはノイズ処理部20の後段においても若干ながら不要成分が存在する。その不要成分(残留不要成分)は、MTIフィルタ部16及びノイズ処理部20によって除去仕切れなかった残留クラッタ成分、及び、ノイズ処理部20によって除去仕切れなかったホワイトノイズ状の残留ノイズ成分、などを含むものである。そのような残留不要成分が速度信号中に含まれると血流画像を形成した場合においてその画質を低下させることになるため、その効果的な除去が望まれる。そのために本実施形態においては残留不要成分処理部22が設けられている。   Originally, all of the clutter components and noise components included in the received signal should be removed by the MTI filter unit 16 and the noise processing unit 20 described above. However, there are unnecessary components. The unnecessary components (residual unnecessary components) include residual clutter components that are not removed by the MTI filter unit 16 and the noise processing unit 20, white noise-like residual noise components that are not removed by the noise processing unit 20, and the like. It is a waste. If such a residual unnecessary component is included in the velocity signal, the image quality is deteriorated when a blood flow image is formed. Therefore, effective removal thereof is desired. Therefore, the residual unnecessary component processing unit 22 is provided in the present embodiment.

それについて具体的に説明すると、ノイズ処理部20から出力される速度値velは三値化部28及び除去部32に出力されている。三値化部28は、本実施形態において、vel>0の場合にはvel’=+1とし、vel<0の場合にはvel’=−1とし、vel=0の場合には属性値vel’=0とする。すなわち、後の相関値演算の前処理として速度値velを規格化処理し、各サンプル点ごとに属性値vel’を求めるものである。このような前処理により後の相関演算を簡略化でき、しかも、血流成分と不要信号成分とをそれぞれの空間的な性質に基づいて効果的に弁別することが可能となる。   More specifically, the velocity value vel output from the noise processing unit 20 is output to the ternarization unit 28 and the removal unit 32. In the present embodiment, the ternary unit 28 sets vel ′ = + 1 when vel> 0, sets vel ′ = − 1 when vel <0, and sets the attribute value vel ′ when vel = 0. = 0. That is, the velocity value vel is standardized as a pre-processing for the subsequent correlation value calculation, and the attribute value vel 'is obtained for each sample point. Such pre-processing can simplify subsequent correlation calculations, and can effectively discriminate blood flow components and unwanted signal components based on their spatial properties.

相関値演算部30は、走査面に対応するフレーム上において注目するサンプル点を中心として注目領域を設定し、また、その注目領域を一定方向にシフトした同形状の領域として参照領域を設定する。そして、注目領域と参照領域との間において後に詳述する相関値演算を実行し、その演算結果である相関値Rを評価値として除去部32へ出力する。   The correlation value calculation unit 30 sets a region of interest around a sample point of interest on a frame corresponding to the scanning plane, and sets a reference region as a region of the same shape that is shifted in a certain direction. Then, correlation value calculation, which will be described in detail later, is performed between the attention area and the reference area, and the correlation value R that is the calculation result is output to the removal unit 32 as an evaluation value.

フレーム上における各サンプル点ごとに相関値Rが求められ、それが除去部32へ出力されることになる。このような相関値の演算を行うために三値化部28と相関値演算部30との間に必要に応じてラインメモリあるいはフレームメモリが設けられる。相関値演算部30の具体的な作用については後に詳述する。   A correlation value R is obtained for each sample point on the frame and is output to the removal unit 32. In order to perform such a correlation value calculation, a line memory or a frame memory is provided between the ternarization unit 28 and the correlation value calculation unit 30 as necessary. The specific operation of the correlation value calculation unit 30 will be described in detail later.

除去部32は、各サンプル点の速度値ごとに、それが血流成分であるか、残留不要成分であるかの判定を行って、残留不要成分を除去する処理を実行する。具体的には、後に詳述するように、速度値velに対して設定されている閾値Vthと相関値Rに対して設定されている閾値Rthとを用いて、二次元の閾値判定を行うことにより血流成分と残留不要成分を弁別し、その弁別結果に基づいて残留不要成分の除去を行っている。残留不要成分の除去を行った場合、必要に応じてダミーデータが挿入される。そのようなダミーデータの挿入は以下に説明する表示処理部24の作用として行わせてもよい。 The removal unit 32 determines whether the velocity value of each sample point is a blood flow component or a residual unnecessary component, and executes a process of removing the residual unnecessary component. Specifically, as will be described in detail later, the two-dimensional threshold determination is performed using the threshold V th set for the velocity value vel and the threshold R th set for the correlation value R. By doing so, the blood flow component and the residual unnecessary component are discriminated, and the residual unnecessary component is removed based on the discrimination result. When residual unnecessary components are removed, dummy data is inserted as necessary. Such insertion of dummy data may be performed as an operation of the display processing unit 24 described below.

いずれにしても、MTIフィルタ部16及びノイズ処理部20の作用によってある程度のクラッタ成分除去及びノイズ成分除去が達成されているが、上記のような残留不要成分処理部22を更に設けることにより、残留不要成分を効果的に除去して高品質の二次元血流画像を形成することが可能となる。本発明者による実験によれば、二次元血流画像の画質をかなり改善できることが確認されている。すなわち、本実施形態においては血流成分と残留不要成分の空間的な性質の違いを利用したため、従来のクラッタ処理及びノイズ処理とは異なる観点から不要成分の除去を行ったため、単に二段階の不要成分除去を行う場合よりも、より効果的な不要成分除去を達成できている。   In any case, a certain amount of clutter component removal and noise component removal has been achieved by the action of the MTI filter unit 16 and the noise processing unit 20, but by providing the residual unnecessary component processing unit 22 as described above, It is possible to effectively remove unnecessary components and form a high-quality two-dimensional blood flow image. According to experiments by the present inventors, it has been confirmed that the image quality of a two-dimensional blood flow image can be significantly improved. That is, in this embodiment, since the difference in the spatial characteristics of the blood flow component and the residual unnecessary component is used, unnecessary components are removed from a viewpoint different from the conventional clutter processing and noise processing, so that only two steps are unnecessary. More effective removal of unnecessary components can be achieved than when components are removed.

表示処理部24は本実施形態においてデジタルスキャンコンバータ(DSC)を含むものである。そのDSCは座標変換、補間処理機能などを有している。表示処理部24によって残留不要成分処理部22から出力される速度信号に基づいて二次元血流画像が形成される。この場合においては必要に応じて残留不要信号の除去によって消失したデータが補間値によって補われる。表示部26には血流画像が表示される。   The display processing unit 24 includes a digital scan converter (DSC) in this embodiment. The DSC has coordinate conversion and interpolation processing functions. A two-dimensional blood flow image is formed based on the velocity signal output from the residual unnecessary component processing unit 22 by the display processing unit 24. In this case, the data lost by the removal of the residual unnecessary signal is supplemented by the interpolation value as necessary. A blood flow image is displayed on the display unit 26.

図1には図示されていないが、二次元血流画像と一緒に二次元断層画像を表示するため、二次元断層画像形成部が設けられる。上記実施形態においては、MTIフィルタ部16及びノイズ処理部20の後段に残留不要成分処理部22が設けられていたが、MTIフィルタ部16及びノイズ処理部20の一方あるいは両方を設けない場合においても、当該構成を採用することができ、そのような場合によっても一定の効果を得ることができる。   Although not shown in FIG. 1, a two-dimensional tomographic image forming unit is provided to display a two-dimensional tomographic image together with the two-dimensional blood flow image. In the above embodiment, the residual unnecessary component processing unit 22 is provided after the MTI filter unit 16 and the noise processing unit 20, but even when one or both of the MTI filter unit 16 and the noise processing unit 20 are not provided. The configuration can be adopted, and a certain effect can be obtained even in such a case.

以下に相関値演算部30の作用について詳述する。上記のように相関値演算部30は、三値化処理後の属性値vel’の集合(本実施形態では一次元配列)に基づいてランダム性の評価値として相関値(あるいはそれに相当する値)を演算する。具体的には、図2におけるW1は、注目するサンプル点jを中心としてビーム上で上下均等に広がった一次元の注目領域を示しており、その垂直方向(ビーム方向)のサイズはNである。図2におけるW1’は、注目領域W1を垂直方向へnだけシフトさせることによって定義される参照領域を示している。ここでnは1以上の整数(ピクセル数)である。nは例えば1であるが、その値を状況に応じて可変設定するようにしてもよい。注目領域W1と参照領域W1’との間で以下に説明する相関値(評価値)の演算が実行される。ここで、注目領域W1と参照領域W1’との間における相関係数(自己相関係数)をRj(n)とすると、以下の(1)式によってそれが定義される。

Figure 0004616614
The operation of the correlation value calculation unit 30 will be described in detail below. As described above, the correlation value calculation unit 30 uses the correlation value (or a value corresponding to it) as an evaluation value of randomness based on the set of attribute values vel ′ after ternarization processing (one-dimensional array in this embodiment). Is calculated. Specifically, W1 in FIG. 2 indicates a one-dimensional region of interest that spreads evenly up and down on the beam around the sample point j of interest, and the size in the vertical direction (beam direction) is N. . W1 ′ in FIG. 2 indicates a reference area defined by shifting the attention area W1 by n in the vertical direction. Here, n is an integer (number of pixels) of 1 or more. For example, n is 1, but the value may be variably set according to the situation. A correlation value (evaluation value) described below is executed between the attention area W1 and the reference area W1 ′. Here, if the correlation coefficient (autocorrelation coefficient) between the attention area W1 and the reference area W1 ′ is R j (n), it is defined by the following equation (1).
Figure 0004616614

残留不要成分についての速度値の空間分布は上述したようにランダム性が強いので、注目領域W1と参照領域W1’との間における相関係数Rj(n)は、血流成分のそれに比べて小さくなる。その性質を利用して、個々のサンプル点の相関係数Rj(n)が閾値Rth以下となる場合には、残留不要成分である可能性が大きいことを認識できる。但し、いわゆる信号処理上の折り返し現象によりモザイク状となった血流部が不要成分であるとして不必要に除外されてしまうことを防止するため、注目するサンプル点の速度値velの絶対値がある閾値Vthよりも高い場合には除去対象から除外するようにする。上記の2つの条件の組み合わせが図3に示されており、ここで横軸は速度値を示しており、縦軸は相関値を示している。ハッチングが付された領域100が残留不要成分として除去される範囲を示している。 Since the spatial distribution of the velocity values for the residual unnecessary component is highly random as described above, the correlation coefficient R j (n) between the attention area W1 and the reference area W1 ′ is larger than that of the blood flow component. Get smaller. Using this property, it can be recognized that there is a high possibility of being a residual unnecessary component when the correlation coefficient R j (n) of each sample point is equal to or less than the threshold value R th . However, there is an absolute value of the velocity value vel of the sample point of interest in order to prevent the mosaic blood flow part from being unnecessarily excluded as an unnecessary component due to the so-called folding phenomenon in signal processing. If it is higher than the threshold value V th, it is excluded from the removal target. A combination of the above two conditions is shown in FIG. 3, where the horizontal axis indicates the velocity value and the vertical axis indicates the correlation value. A hatched region 100 is shown as a range in which residual unnecessary components are removed.

以上の処理がフレームを構成する各サンプル点(あるいはピクセル点)ごとに実行され、フレーム内に含まれる残留不要成分が効果的に取り除かれた情報を利用して血流画像を形成することが可能となる。ちなみに、残留不要成分を取り除いた後に別のダミーデータを補間処理によって補填するようにしてもよい。そのような補間処理は、デジタルスキャンコンバータにおける変換処理に委ねることもできる。   The above processing is executed for each sample point (or pixel point) constituting the frame, and a blood flow image can be formed by using information from which residual unnecessary components contained in the frame are effectively removed. It becomes. Incidentally, after removing unnecessary residual components, another dummy data may be supplemented by interpolation processing. Such an interpolation process can be left to the conversion process in the digital scan converter.

上記の(1)式に代えて、以下の(2)式又は(3)式を用いることもできる。それらの式は厳密には相関値を求める演算式ではないが、それと等価な値を求めることができる式である。

Figure 0004616614
Instead of the above formula (1), the following formula (2) or formula (3) can also be used. Strictly speaking, these expressions are not arithmetic expressions for obtaining correlation values, but are expressions that can obtain equivalent values.
Figure 0004616614

注目領域W1と参照領域W1’との間で相関関係を求めることが可能であれば上記以外の計算式を利用することができる。更に、上記の(2)式及び(3)式におけるRj(0)以下の(4)式に示すものとしてもよい。

Figure 0004616614
If it is possible to obtain a correlation between the attention area W1 and the reference area W1 ′, a calculation formula other than the above can be used. Furthermore, it is good also as what is shown to (4) Formula below Rj (0) in said (2) Formula and (3) Formula.
Figure 0004616614

上記の(2)式を用いた場合には、図3に示した判定条件と同じものを利用することができ、上記の(3)式を用いた場合には、ランダム性が大きい信号であるほどRj(n)が大きくなるため、図4に示すような判定条件が利用される。図4において、横軸は速度値を示しており、縦軸は係数を示している。係数Rj(n)が閾値Rth以上で、且つ、速度値velの絶対値が閾値Vth以下の場合には不要成分であると判定される。それがハッチングが付された領域102として示されている。その領域102に該当するサンプル点の速度値は除去され(ゼロに修正され)、それ以外の領域に該当するサンプル点の速度値はそのまま保存される。 When the above equation (2) is used, the same judgment conditions as those shown in FIG. 3 can be used. When the above equation (3) is used, the signal is highly random. Since R j (n) increases, determination conditions as shown in FIG. 4 are used. In FIG. 4, the horizontal axis indicates the speed value, and the vertical axis indicates the coefficient. When the coefficient R j (n) is equal to or greater than the threshold value R th and the absolute value of the velocity value vel is equal to or less than the threshold value V th, it is determined as an unnecessary component. This is shown as hatched area 102. The velocity value of the sample point corresponding to the region 102 is removed (corrected to zero), and the velocity value of the sample point corresponding to the other region is stored as it is.

上記においては、注目領域及び参照領域がビーム方向に沿った一次元の領域であったが、ビーム走査方向に一次元の注目領域及び参照領域を設定することもできるし、以下に説明するように二次元の注目領域及び参照領域を設定することもできる。図5においては、注目するサンプル点Pを中心として注目領域W2が定義され、それを水平方向へm且つ垂直方向へnだけシフトした領域として同一形態の参照領域W2’が定義されている。ここで、m,nはいずれも1以上の整数である、なお、Bはビームライン(音線)を示している。図6においては、注目領域W3とそれに対して水平方向にmだけシフトした参照領域W3’とが定義され、図7においては、注目領域W4とそれに対して垂直方向にnだけシフトした参照領域W4’とが定義されている。ここで、mは1以上の整数であり、nも1以上の整数である。このように二次元の領域のシフトを利用して相関値を演算することもできる。   In the above, the attention area and the reference area are one-dimensional areas along the beam direction. However, a one-dimensional attention area and a reference area can be set in the beam scanning direction, and will be described below. A two-dimensional attention area and reference area can also be set. In FIG. 5, a region of interest W2 is defined around the sample point P of interest, and a reference region W2 'having the same form is defined as a region shifted by m in the horizontal direction and n in the vertical direction. Here, m and n are both integers of 1 or more, and B indicates a beam line (sound ray). In FIG. 6, the attention area W3 and a reference area W3 ′ shifted by m in the horizontal direction are defined, and in FIG. 7, the attention area W4 and the reference area W4 shifted by n in the vertical direction are defined. 'And are defined. Here, m is an integer of 1 or more, and n is an integer of 1 or more. In this way, the correlation value can be calculated using the shift of the two-dimensional region.

ここで、図5の場合における相関係数Ri,j(m,n)は、領域の大きさをM×Nとして、以下の(5)式のように求められる。これは上記(1)式に対応したものである。

Figure 0004616614
Here, the correlation coefficient R i, j (m, n) in the case of FIG. 5 is obtained as in the following equation (5), where the size of the region is M × N. This corresponds to the above equation (1).
Figure 0004616614

また、同様に上記(2)式及び(3)式に対応する式は以下の通りである。

Figure 0004616614
Similarly, the expressions corresponding to the expressions (2) and (3) are as follows.
Figure 0004616614

この場合においても、(7)式及び(8)式に含まれるRi,j(0,0)を以下のように置き換えることができる。

Figure 0004616614
Even in this case, R i, j (0,0) included in the equations (7) and (8) can be replaced as follows.
Figure 0004616614

以上説明したように、本実施形態によれば、MTIフィルタ部16及びノイズ処理部20によって除去仕切れなかった残留不要成分を効果的に除去して二次元血流画像の画質を高めることが可能である。すなわち、残留不要成分についての空間的な性質に着目してそれを血流成分と精度良く弁別して残留不要成分を速度信号中から精度良く除外することができる。ちなみに、三値化部28における規格化処理においては、三値化処理に代えて二値化処理や四値化処理などを採用することもできるが、残留不要成分の弁別という観点及び演算の簡略化の観点から言えば、三値化処理を行うのが望ましい。また相関値演算部30において採用する演算式としては空間的なランダム性を評価できる限りにおいて上記で掲げた計算式以外の計算式を採用することもできる。   As described above, according to the present embodiment, it is possible to effectively remove residual unnecessary components that have not been removed by the MTI filter unit 16 and the noise processing unit 20 and to improve the image quality of the two-dimensional blood flow image. is there. That is, focusing on the spatial properties of the residual unnecessary component, it can be accurately distinguished from the blood flow component, and the residual unnecessary component can be accurately excluded from the speed signal. Incidentally, in the standardization processing in the ternarization unit 28, binarization processing, quaternization processing, etc. can be adopted instead of the ternarization processing. However, from the viewpoint of discriminating residual unnecessary components, the calculation is simplified. From the viewpoint of conversion, it is desirable to perform ternary processing. In addition, as a calculation formula employed in the correlation value calculation unit 30, calculation formulas other than the calculation formulas described above may be employed as long as spatial randomness can be evaluated.

本発明に係る超音波診断装置の全体構成を示すブロック図である。1 is a block diagram showing an overall configuration of an ultrasonic diagnostic apparatus according to the present invention. 一次元の注目領域及び参照領域を示す図である。It is a figure which shows a one-dimensional attention area and a reference area. 不要信号判定条件の一例を示す図である。It is a figure which shows an example of an unnecessary signal determination condition. 不要信号判定条件の他の例を示す図である。It is a figure which shows the other example of unnecessary signal determination conditions. 水平方向及び垂直方向にシフトした関係にある二次元の注目領域及び参照領域を示す図である。It is a figure which shows the two-dimensional attention area | region and reference area | region which have the relationship shifted to the horizontal direction and the vertical direction. 水平方向にシフトした関係にある二次元の注目領域及び参照領域を示す図である。It is a figure which shows the two-dimensional attention area | region and reference area | region which have the relationship shifted to the horizontal direction. 垂直方向にシフトした関係にある二次元の注目領域及び参照領域を示す図である。It is a figure which shows the two-dimensional attention area | region and reference area | region which have the relationship shifted to the perpendicular direction.

符号の説明Explanation of symbols

16 MTIフィルタ部、20 ノイズ処理部、22 残留不要成分処理部、28
三値化部、30 相関値演算部、32 除去部。
16 MTI filter section, 20 noise processing section, 22 residual unnecessary component processing section, 28
Trinization unit, 30 correlation value calculation unit, 32 removal unit.

Claims (10)

超音波を送受波して得られた受信信号に基づいて、各サンプル点の速度値を演算する速度演算手段と、
前記各サンプル点の速度値を三値化処理し、前記各サンプル点ごとに属性値を求める三値化手段と、
前記三値化処理後の空間的な属性値配列に基づいて、前記各サンプル点を基準として設定される注目領域について空間的なランダム性を評価するためのランダム性評価演算を実行し、前記各サンプル点ごとに評価値を求める評価手段と、
前記各サンプル点ごとの評価値に基づいて、前記各サンプル点の速度値が血流成分であるか不要成分であるかを判定し、不要成分を除去する処理を実行する除去手段と、
前記不要信号を除去する処理を実行した後に、血流成分に基づいて血流画像を形成する血流画像形成手段と、
を含み、
前記評価手段は、前記注目領域とそれを空間的にシフトさせた参照領域との間で、前記ランダム性評価演算を実行することを特徴とする超音波診断装置。
Based on a received signal obtained by transmitting and receiving ultrasonic waves, speed calculating means for calculating a speed value of each sample point;
A ternarization means for ternarizing the velocity value of each sample point and obtaining an attribute value for each sample point;
Based on the spatial attribute value array after the ternarization processing, a randomness evaluation calculation is performed for evaluating spatial randomness for a region of interest set with reference to each sample point, and An evaluation means for obtaining an evaluation value for each sample point;
Based on the evaluation value for each sample point, a removal unit that determines whether the velocity value of each sample point is a blood flow component or an unnecessary component, and executes processing for removing the unnecessary component;
Blood flow image forming means for forming a blood flow image based on a blood flow component after executing the process of removing the unnecessary signal;
Only including,
The ultrasonic diagnostic apparatus , wherein the evaluation unit executes the randomness evaluation calculation between the region of interest and a reference region obtained by spatially shifting the region of interest .
請求項1記載の装置において、
前記三値化処理では、前記各サンプル点の速度値が第1条件を満たす場合に正属性値が付与され、前記各サンプル点の速度値が第2条件を満たす場合にゼロ属性値が付与され、前記各サンプル点の速度値が第3条件を満たす場合に負属性値が付与されることを特徴とする超音波診断装置。
The apparatus of claim 1.
In the ternary processing, a positive attribute value is assigned when the velocity value of each sample point satisfies the first condition, and a zero attribute value is assigned when the velocity value of each sample point satisfies the second condition. The ultrasonic diagnostic apparatus is characterized in that a negative attribute value is given when the velocity value of each sample point satisfies the third condition.
請求項1記載の装置において、
前記注目領域からの前記参照領域のシフト量を可変設定する手段を含むことを特徴とする超音波診断装置。
The apparatus of claim 1.
An ultrasonic diagnostic apparatus comprising: means for variably setting a shift amount of the reference area from the attention area .
超音波を送受波して得られた受信信号に基づいて、各サンプル点の速度値を演算する速度演算手段と、
前記各サンプル点の速度値を三値化処理し、前記各サンプル点ごとに属性値を求める三値化手段と、
前記三値化処理後の空間的な属性値配列に基づいて、前記各サンプル点を基準として設定される注目領域について空間的なランダム性を評価するためのランダム性評価演算を実行し、前記各サンプル点ごとに評価値を求める評価手段と、
前記各サンプル点ごとの評価値に基づいて、前記各サンプル点の速度値が血流成分であるか不要成分であるかを判定し、不要成分を除去する処理を実行する除去手段と、
前記不要信号を除去する処理を実行した後に、血流成分に基づいて血流画像を形成する血流画像形成手段と、
を含み、
前記ランダム性評価演算は相関演算に相当し、
前記評価値は相関値に相当することを特徴とする超音波診断装置。
Based on the received signal obtained by transmitting and receiving ultrasonic waves, speed calculation means for calculating the speed value of each sample point,
A ternarization means for ternarizing the velocity value of each sample point and obtaining an attribute value for each sample point;
Based on the spatial attribute value array after the ternarization processing, a randomness evaluation calculation is performed for evaluating spatial randomness for a region of interest set with each sample point as a reference, An evaluation means for obtaining an evaluation value for each sample point;
Based on the evaluation value for each sample point, a removal unit that determines whether the velocity value of each sample point is a blood flow component or an unnecessary component, and executes processing for removing the unnecessary component;
Blood flow image forming means for forming a blood flow image based on a blood flow component after executing the process of removing the unnecessary signal;
Including
The randomness evaluation calculation corresponds to a correlation calculation,
The ultrasonic diagnostic apparatus, wherein the evaluation value corresponds to a correlation value.
請求項3記載の装置において、
前記注目領域はビーム方向に沿った一次元領域であり、
前記参照領域は前記注目領域をビーム方向にシフトさせた領域であることを特徴とする超音波診断装置。
The apparatus of claim 3.
The region of interest is a one-dimensional region along the beam direction;
The ultrasonic diagnostic apparatus, wherein the reference region is a region obtained by shifting the region of interest in a beam direction.
請求項3記載の装置において、
前記注目領域はフレーム上の二次元領域であり、
前記参照領域は前記注目領域をフレーム上でシフトさせた領域であることを特徴とする超音波診断装置。
The apparatus of claim 3.
The region of interest is a two-dimensional region on a frame;
The ultrasonic diagnostic apparatus, wherein the reference region is a region obtained by shifting the region of interest on a frame.
超音波を送受波して得られた受信信号に基づいて、各サンプル点の速度値を演算する速度演算手段と、
前記各サンプル点の速度値を三値化処理し、前記各サンプル点ごとに属性値を求める三値化手段と、
前記三値化処理後の空間的な属性値配列に基づいて、前記各サンプル点を基準として設定される注目領域について空間的なランダム性を評価するためのランダム性評価演算を実行し、前記各サンプル点ごとに評価値を求める評価手段と、
前記各サンプル点ごとの評価値に基づいて、前記各サンプル点の速度値が血流成分であるか不要成分であるかを判定し、不要成分を除去する処理を実行する除去手段と、
前記不要信号を除去する処理を実行した後に、血流成分に基づいて血流画像を形成する血流画像形成手段と、
を含み、
前記除去手段は、前記各サンプル点ごとに、それが有する速度値とそれについて演算された評価値との組み合わせに基づいて、前記各サンプル点の速度値が血流成分であるか不要成分であるかを判定することを特徴とする超音波診断装置。
Based on a received signal obtained by transmitting and receiving ultrasonic waves, speed calculating means for calculating a speed value of each sample point;
A ternarization means for ternarizing the velocity value of each sample point and obtaining an attribute value for each sample point;
Based on the spatial attribute value array after the ternarization processing, a randomness evaluation calculation is performed for evaluating spatial randomness for a region of interest set with reference to each sample point, and An evaluation means for obtaining an evaluation value for each sample point;
Based on the evaluation value for each sample point, a removal unit that determines whether the velocity value of each sample point is a blood flow component or an unnecessary component, and executes processing for removing the unnecessary component;
Blood flow image forming means for forming a blood flow image based on a blood flow component after executing the process of removing the unnecessary signal;
Including
The removing means is a blood flow component or an unnecessary component for each sample point based on a combination of a velocity value of the sample point and an evaluation value calculated for the sample point. An ultrasonic diagnostic apparatus characterized by determining whether or not.
請求項7記載の装置において、
前記除去手段は、前記速度値が第1の判定条件を満たし、且つ、前記評価値が第2の判定条件を満たす場合に不要成分であると判定することを特徴とする超音波診断装置。
The apparatus of claim 7.
The ultrasonic diagnostic apparatus according to claim 1, wherein the removing unit determines that the velocity value satisfies a first determination condition and the evaluation value is an unnecessary component when the evaluation value satisfies a second determination condition.
超音波を送受波して得られた受信信号に含まれるクラッタ成分を除去するクラッタ成分除去手段と、
前記クラッタ成分の除去後における受信信号に基づいて、各サンプル点の速度値を演算する速度演算手段と、
前記各サンプル点ごとの速度値によって構成される速度信号に含まれるノイズ成分を除去するためのノイズ成分除去手段と、
前記各サンプル点の速度値に基づいて、前記各サンプル点ごとに属性値を求める属性値演算手段と、
前記属性値演算後の空間的な属性値配列に基づいて、前記各サンプル点を基準として設定される注目領域について空間的なランダム性を評価するためのランダム性評価演算を実行し、前記各サンプル点ごとに評価値を求める評価手段と、
前記各サンプル点ごとの速度値及び評価値に基づいて、前記各サンプル点の速度値が血流成分であるか残留不要成分であるかを判定し、残留不要成分を除去する処理を実行する残留不要成分除去手段と、
前記残留不要成分を除去する処理を実行した後に、血流成分に基づいて血流画像を形成する血流画像形成手段と、
を含み、
前記評価手段は、前記注目領域とそれを空間的にシフトさせた参照領域との間で、相関演算に相当する前記ランダム性評価演算を実行することを特徴とする超音波診断装置。
A clutter component removing means for removing clutter components contained in a received signal obtained by transmitting and receiving ultrasonic waves;
Based on the received signal after removal of the clutter component, speed calculation means for calculating the speed value of each sample point;
Noise component removing means for removing a noise component contained in a velocity signal constituted by velocity values for the respective sample points;
Attribute value calculation means for obtaining an attribute value for each sample point based on the velocity value of each sample point;
Based on the spatial attribute value array after the attribute value calculation, a randomness evaluation calculation is performed to evaluate spatial randomness for a region of interest set based on each sample point, and each sample An evaluation means for obtaining an evaluation value for each point;
Based on the velocity value and the evaluation value for each sample point, it is determined whether the velocity value of each sample point is a blood flow component or a residual unnecessary component, and a residual is executed to remove the residual unnecessary component Unnecessary component removing means;
Blood flow image forming means for forming a blood flow image based on a blood flow component after executing the process of removing the residual unnecessary component;
Only including,
The ultrasonic diagnostic apparatus , wherein the evaluation unit performs the randomness evaluation calculation corresponding to a correlation calculation between the attention area and a reference area obtained by spatially shifting the attention area .
請求項9記載の装置において、
前記属性値演算手段は、前記各サンプル点の速度値が第1条件を満たす場合に正属性値を付与し、前記各サンプル点の速度値が第2条件を満たす場合にゼロ属性値を付与し、前記各サンプル点の速度値が第3条件を満たす場合に負属性値を付与する三値化処理手段であることを特徴とする超音波診断装置。
The apparatus of claim 9.
The attribute value calculation means gives a positive attribute value when the velocity value of each sample point satisfies the first condition, and gives a zero attribute value when the velocity value of each sample point satisfies the second condition. The ultrasonic diagnostic apparatus is a ternary processing means for giving a negative attribute value when the velocity value of each sample point satisfies the third condition.
JP2004306926A 2004-10-21 2004-10-21 Ultrasonic diagnostic equipment Expired - Fee Related JP4616614B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2004306926A JP4616614B2 (en) 2004-10-21 2004-10-21 Ultrasonic diagnostic equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2004306926A JP4616614B2 (en) 2004-10-21 2004-10-21 Ultrasonic diagnostic equipment

Publications (2)

Publication Number Publication Date
JP2006116070A JP2006116070A (en) 2006-05-11
JP4616614B2 true JP4616614B2 (en) 2011-01-19

Family

ID=36534504

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2004306926A Expired - Fee Related JP4616614B2 (en) 2004-10-21 2004-10-21 Ultrasonic diagnostic equipment

Country Status (1)

Country Link
JP (1) JP4616614B2 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4291833B2 (en) 2006-06-23 2009-07-08 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ultrasonic diagnostic apparatus and ultrasonic diagnostic image generation method
JP5366678B2 (en) * 2009-06-25 2013-12-11 株式会社東芝 Three-dimensional ultrasonic diagnostic apparatus and program
JP6580915B2 (en) * 2015-09-14 2019-09-25 キヤノンメディカルシステムズ株式会社 Ultrasonic diagnostic apparatus and signal processing apparatus

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07108043A (en) * 1993-10-13 1995-04-25 Aloka Co Ltd Ultrasonic blood flow observing device

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07108043A (en) * 1993-10-13 1995-04-25 Aloka Co Ltd Ultrasonic blood flow observing device

Also Published As

Publication number Publication date
JP2006116070A (en) 2006-05-11

Similar Documents

Publication Publication Date Title
US8864671B2 (en) Methods and systems for color flow imaging
KR100742467B1 (en) Method and apparatus for visualization of motion in ultrasound flow imaging using packet data acquisition
JP5160825B2 (en) Ultrasonic diagnostic apparatus and image processing program
CN111265246B (en) Ultrasonic color imaging processing method and device
CN102105107A (en) Signal processing apparatus, ultrasonic apparatus and method for detecting a unique region such as a reflector of high reflectivity
JP2006116149A (en) Ultrasonic diagnostic apparatus and control program thereof
KR101051555B1 (en) Ultrasonic Imaging Apparatus and Method for Forming an Improved 3D Ultrasound Image
EP1050761B1 (en) Color imaging display method and apparatus for color doppler imaging system
EP1796037B1 (en) Image processing system and method for controlling gains for color flow images
US20150173717A1 (en) Ultrasonic diagnostic device
EP2103954A1 (en) Adaptive clutter signal filtering in an ultrasound system
Lovstakken et al. Blood flow imaging-a new real-time, flow imaging technique
JP4616614B2 (en) Ultrasonic diagnostic equipment
JP4778325B2 (en) Ultrasonic diagnostic equipment
KR20080060625A (en) Ultrasound diagnostic system and method for acquiring ultrasound images based on motion of a target object
EP1279048A1 (en) Ultrasound system with b/color priority threshold calculation
JP3685737B2 (en) Ultrasonic diagnostic equipment
KR100998443B1 (en) Ultrasound system and method for processing ultrasound data considering scan conversion
JP3662835B2 (en) Ultrasonic diagnostic equipment
JP2723453B2 (en) Ultrasound blood flow observation device
JP2013174473A (en) Underwater detection device and program
JP5999997B2 (en) Ultrasonic diagnostic equipment
JP5936857B2 (en) Ultrasonic diagnostic apparatus and image processing program thereof
JP3802450B2 (en) Ultrasonic diagnostic equipment
JP7471257B2 (en) ULTRASONIC IMAGING APPARATUS AND METHOD FOR GENERATING COLOR DOPPLER IMAGE

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20070913

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20100628

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20100727

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20100909

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

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20101022

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

Free format text: PAYMENT UNTIL: 20131029

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees