Ultrasonic diagnostic apparatus

Abstract

An ultrasonic diagnostic apparatus, including a transceiver for sequentially generating multiple received signals, a color flow mapping signal processing unit that generates blood flow velocity data based on the multiple received signals, and a residual image processing unit that includes an aliasing determination unit and a persistence calculation unit and performs residual image processing on blood flow velocity data, an aliasing determination unit that adds the blood flow velocity data V.sub.current for the most recent frame and the blood flow velocity data V.sub.out-1 for the frame preceding the most recent frame, performs aliasing determination on the frame that was determined by a frame before the most recent frame based on the aliasing determination results of the blood flow velocity data, and a persistence calculation unit that performs a persistence calculation after performing a correction to the blood flow velocity data in which aliasing has occurred.

Claims

1. An ultrasonic diagnostic apparatus comprising: a transceiver that transmits ultrasonic waves by driving a probe repeatedly, receives echoes reflected off a subject, and successively generates a plurality of received signals; a color flow mapping signal processing unit that successively generates blood flow velocity data of blood flow in the subject for each frame of a plurality of frames based on the plurality of received signals; a residual image processing unit that successively receives the each frame of the plurality of frames and performs residual image processing on blood flow velocity of the each frame; the residual image processing unit further comprising; an aliasing determination unit that performs a determination of whether or not aliasing occurs in a first blood flow velocity data for a most recent frame of the plurality of frames as received from the color flow mapping signal processing unit and in a second blood flow velocity data for a frame of the plurality of frames immediately preceding the most recent frame that has previously been subjected to the residual image processing; and a persistence memory unit that stores the second blood flow velocity data for the frame immediately preceding the most recent frame that has been subjected to the residual image processing of the residual processing unit and a first aliasing determination result that has been determined for the second blood flow velocity data for the frame immediately preceding the most recent frame, wherein the determination by the aliasing determination unit is performed based on the first blood flow velocity data, the second blood flow velocity data, and whether the first aliasing determination result that has been previously determined for the frame immediately preceding the most recent frame indicates presence of aliasing in the blood flow velocity data for the frame immediately preceding the most recent frame that has been subjected to the residual image processing, the aliasing determination unit outputting a second aliasing determination result that is an aliasing determination result of the first blood flow velocity data, and a third aliasing determination result that is an aliasing determination result of the second blood flow velocity data, and a persistence calculation unit that, based on the second and third aliasing determination results and a predetermined persistence coefficient for adjusting the residual image effect, performs a persistence calculation using the first blood flow velocity data and the second blood flow velocity data, outputs a persistence calculation result as a third blood flow velocity data for the most recent frame that is the blood flow velocity data that is subjected to residual image processing, determines whether or not aliasing occurs in the third blood flow velocity data, outputs data which is a fourth aliasing determination result, and outputs the third blood flow velocity data, which is a result of the persistence calculation, and the fourth aliasing determination result to the persistence memory unit.

2. The ultrasonic diagnostic apparatus according to claim 1, wherein the persistence calculation unit performs the persistence calculation after correcting the blood flow velocity data of the most recent frame and a residual image processed blood flow velocity data of the frame immediately preceding the most recent frame based on the second and third aliasing determination results.

3. The ultrasonic diagnostic apparatus according to claim 1, wherein the residual image processing unit includes a persistence coefficient determination unit that sets the persistence coefficients based on the second and third aliasing determination results and the first blood flow velocity data.

4. The ultrasonic diagnostic apparatus according to claim 1, the residual image processing unit further comprising: a frame memory unit that stores the first blood flow velocity data, and wherein an aliasing determination unit performs aliasing determination on the first blood flow velocity data and the second blood flow velocity data using the first aliasing determination result, which are obtained from the frame memory unit and the persistence memory unit respectively.

5. The ultrasonic diagnostic apparatus according to claim 1, wherein the aliasing determination unit determines whether or not aliasing occurs in the first blood flow velocity data and the second blood flow velocity data based on a result of comparing magnitudes of the first blood flow velocity data and second blood flow velocity data with a predetermined threshold value, and the first aliasing determination result.

6. The ultrasonic diagnostic apparatus according to claim 1, the residual image processing unit further comprising: a first persistence coefficient determination unit that sets a first persistence coefficient based on the second and third aliasing determination results and the first blood flow velocity data; a second persistence coefficient determination unit that sets a second persistence coefficient based on the second and third aliasing determination results, and the second blood flow velocity data; a first persistence calculation unit that performs a persistence calculation using the first blood flow velocity data and second blood flow velocity data based on the second and third aliasing determination results and the first persistence coefficient; a second persistence calculation unit that performs a persistence calculation using the first blood flow velocity data and second blood flow velocity data based on the second and third aliasing determination results and the second persistence coefficient; a maximum value selection unit that compares absolute values of calculation results outputted from the first and second persistence calculation units, whereby a larger calculation result is outputted as the third blood flow velocity data, determines whether or not aliasing occurs in the third blood flow velocity data, and outputs a fourth aliasing determination result that is an aliasing determination result thereof.

7. The ultrasonic diagnostic apparatus according to claim 1, the residual image processing unit further comprising: a first persistence coefficient determination unit that sets a first persistence coefficient based on the second and third aliasing determination results and the first blood flow velocity data; a second persistence coefficient determination unit that sets a second persistence coefficient based on the second and third aliasing determination results, and the second blood flow velocity data, and a persistence coefficient selection unit that compares an absolute value of the first blood flow velocity data with an absolute value of the second blood flow velocity data, and selects the first persistence coefficient when the absolute value of the first blood flow velocity data is larger and otherwise selects the second persistence coefficient, wherein the persistence calculation unit performs a persistence calculation using the first blood flow velocity data and the second blood flow velocity data based on the second and third aliasing determination results and the selected persistence coefficient.

Description

BRIEF DESCRIPTION OF THE DRAWING

(1) The advantages and features provided by one or more embodiments of the invention will become more fully understood from the detailed description given hereinbelow and the appended drawings which are given by way of illustration only, and thus are not intended as a definition of the limits of the present invention.

(2) FIG. 1 is a block diagram showing an example configuration of an ultrasonic diagnostic apparatus according to the first embodiment of the present invention;

(3) FIG. 2 is a schematic diagram illustrating the blood flow velocity and blood flow direction obtained based on the Doppler Effect;

(4) FIG. 3A shows a complex plane V.sub.current and V.sub.out-1 of condition (10);

(5) FIG. 3B shows a complex plane V.sub.current and V.sub.out-1 of condition (11);

(6) FIG. 3C shows a complex plane V.sub.current and V.sub.out-1 of condition (12);

(7) FIG. 3D shows a complex plane V.sub.current and V.sub.out-1 of condition (13);

(8) FIG. 3E shows a complex plane V.sub.current and V.sub.out-1 of condition (14);

(9) FIG. 3F shows a complex plane V.sub.current and V.sub.out-1 of condition (15);

(10) FIG. 4 is a graph showing an example of the relationship between the reference index and the persistence coefficient;

(11) FIG. 5 is a block diagram showing an example configuration of an ultrasonic diagnostic apparatus according to the second embodiment of the present invention;

(12) FIG. 6A is a graph showing an example of the relationship between the first reference index and the first persistence coefficient;

(13) FIG. 6B is a graph showing an example of the correspondence relationship between the second reference index and the second persistence coefficient;

(14) FIG. 7 is a block diagram showing an example configuration of an ultrasonic diagnostic apparatus according to the third embodiment of the present invention;

(15) FIG. 8 is a block diagram showing an example configuration of an existing ultrasonic diagnostic apparatus.

DETAILED DESCRIPTION OF EMBODIMENTS

(16) Hereinafter, one or more embodiments of the present invention will be described with reference to the drawings. However, the scope of the invention is not limited to the disclosed embodiments.

First Embodiment

Example Configuration of the Ultrasound Diagnostic Apparatus

(17) Hereinafter, an ultrasonic diagnostic apparatus according to the first embodiment of the present invention will be described with reference to the drawings. FIG. 1 is a block diagram showing an example configuration of an ultrasonic diagnostic apparatus according to the first embodiment of the present invention.

(18) The ultrasonic diagnostic apparatus 11 shown in FIG. 1 includes a probe 101, an ultrasonic transceiver unit 102, a CFM signal processing unit 103, a residual image processing unit 115, a tomographic image signal processing unit 111, a CFM DSC unit 110, a tomographic image DSC unit 112, an image combining unit 113, and a display 114. Among these configurations, it is possible to use a general-purpose probe and a display device for the probe 101 and display 114, which is not necessary for the ultrasonic diagnostic apparatus 11 to include the probe 101 and display 114.

(19) The ultrasonic transceiver unit 102 generates a drive signal for driving the probe 101 and outputs it to the probe 101. In response to this, the probe 101 transmits ultrasonic waves toward the subject. The probe 101 receives the reflected echo of the transmitted ultrasonic waves reflected off of the subject and generates reception signals. The ultrasonic transceiver unit 102 corresponds to the ultrasonic transceiver unit of the present invention.

(20) More specifically, the probe 101 includes multiple piezoelectric elements, and ultrasonic waves are transmitted from each piezoelectric element. An ultrasonic beam is formed by ultrasonic waves transmitted from multiple piezoelectric elements. The ultrasonic transceiver unit 102 drives the probe 101 while performing delay control of each piezoelectric element so as to scan the subject with multiple ultrasonic beams.

(21) Reflected echos from the subject are received by each piezoelectric element. The ultrasonic transceiver unit 102 has the probe 101 generate a reception signal that corresponds to the transmitted ultrasonic beam by controlling the delay of each piezoelectric element. When the ultrasonic beam scans the subject once, data for one frame is obtained. The probe 101 can sequentially generate reception signals of several frames to several tens of frames per second by scanning the subject several times to several tens of times per second.

(22) The ultrasonic diagnostic apparatus 11 according to the first embodiment generates a B-mode tomographic image and a color flow mapping image, combines them, and displays them on the display 114. Therefore, the transmission and reception of ultrasonic waves with the ultrasonic transceiver unit 102 is performed for generating both the B-mode tomographic image and color flow mapping image. The number of frames per second of the B-mode tomographic image and the number of frames per second of the color flow mapping image may be the same or different. When the number of frames is the same, transmission and reception of ultrasonic waves for B-mode tomographic image generation and transmission and reception of ultrasonic waves for color flow mapping image generation may be alternately repeated.

(23) When generating a B-mode tomographic image, the ultrasonic transceiver unit 102 performs transmission and reception suitable for generating a B-mode tomographic image and outputs the obtained reception signal to the tomographic image signal processing unit 111. On the other hand, when generating a color flow mapping tomographic image, the ultrasonic transceiver unit 102 performs transmission and reception suitable for generating a color flow mapping tomographic image and outputs the obtained reception signal to the CFM signal processing unit 103. Generally, when generating a color flow mapping tomographic image, in order to obtain a stable color flow mapping tomographic image, the ultrasonic transceiver unit 102 performs transmission and reception of ultrasonic waves on the same acoustic line multiple times.

(24) The CFM signal processing unit 103, after performing orthogonal detection processing, MTI filtering, and auto-correlation processing on the received signal, calculates blood flow velocity and blood flow power, and then performs noise filtering to eliminate system or acoustic noise, and generates CFM frame data. CFM frame data includes at least blood flow velocity data and may also include blood flow power data and dispersion data of the blood flow velocity. The CFM signal processing unit 103 repeats the process for each received signal forming each frame. CFM frame data generated by the CFM signal processing unit 103 is output to the residual image processing unit 115 for each frame. The CFM signal processing unit 103 corresponds to the color flow mapping signal processing unit of the present invention.

(25) The residual image processing unit 115, using the persistence coefficient, performs residual image processing on the CFM frame data for each frame. In the ultrasonic diagnostic apparatus 11 according to the first embodiment, the residual image processing unit 115 sets the persistence coefficient according to the blood flow velocity of the most recent frame. In other words, the persistence coefficient is not constant and is a dynamic value based on the blood flow velocity of the most recent frame. As a result, the persistence coefficient can be changed according to the blood flow velocity, and the residual image effect can be adjusted. However, blood flow velocity that can be measured by Doppler shift is restricted by the pulse repetition frequency (PRF). As a result, aliasing, which is a phenomenon in which the blood flow velocity corresponding to the frequency change exceeding ±PRF/2 is observed as blood flow velocity in the opposite direction, occurs and it becomes difficult to accurately evaluate the blood flow velocity. Therefore, the residual image processing unit 115 determines whether or not aliasing is occurring, and if it is determined that aliasing is occurring, corrects the blood flow velocity data, thereby, canceling the effect of the aliasing. Thus, in the ultrasonic diagnostic apparatus 11 according to the first embodiment, accurate blood flow velocity can be assessed regardless of the occurrence of aliasing.

(26) In the ultrasonic diagnostic apparatus 11 according to the first embodiment, the residual image processing unit 115 uses the blood flow velocity data of the most recent frame and the blood flow velocity data of the immediately preceding frame to determine whether or not aliasing has occurred. As a configuration for performing such processing, the residual image processing unit 115 includes a frame memory unit (first memory unit) 104, an aliasing determination unit 105, a persistence coefficient determination unit 106, a persistence coefficient reference memory unit (third memory unit) 107, a persistence calculation unit 108, and a persistence memory unit (second memory unit) 109.

(27) The frame memory unit 104 stores the CFM frame data of the most recent (currently being scanned) frame. The persistence memory unit 109 stores the CFM frame data and aliasing determination result outputted by the persistence calculation unit 108 for the frame immediately preceding the most recent frame. The CFM frame data stored in the persistence memory unit 109 is subjected to the persistence calculation, that is, residual image processing. In the following, blood flow velocity data of the CFM frame data stored in the frame memory unit 104 will be referred to as V.sub.current, blood flow velocity data of CFM frame data stored in the persistence memory unit 109 will be referred to as V.sub.out-1, and the aliasing determination result stored in the persistence memory unit 109 for the frame before the most recent frame will be referred to as R.sub.aliasing-1.

(28) The aliasing determination unit 105 performs aliasing determination on both the blood flow velocity data of CFM frame data from the frame memory unit 104, V.sub.current, and the blood flow velocity data of CFM frame data from the persistence memory unit 109, V.sub.out-1. More specifically, the aliasing determination unit 105 compares the blood flow velocity data V.sub.current and blood flow velocity data V.sub.out-1 with a plurality of threshold values respectively, so as to determine whether or not aliasing has occurred in the blood flow velocity data for the most recent frame V.sub.current and the blood flow velocity data for the frame immediately preceding the most recent frame V.sub.out-1, and outputs the determination results to the persistence coefficient determination unit 106 and the persistence calculation unit 108.

(29) The persistence coefficient determination unit 106, based on the two determination results from the aliasing determination unit 105 (the determination results of V.sub.current and V.sub.out-1) and the blood flow velocity data V.sub.current read from the frame memory unit 104, creates a reference index to the persistence coefficient reference memory unit 107. The persistence coefficient determination unit 106 accesses the persistence coefficient reference memory unit 107, reads the persistence coefficients mapped to the reference index, and sets them in the persistence calculation unit 108. In the persistence coefficient reference memory unit 107, a reference table of persistence coefficients corresponding to the values of blood flow velocity is stored in advance. The reference table includes two or more different persistence coefficients mapped according to the values of blood flow velocity.

(30) The persistence calculation unit 108 performs the persistence calculation on the most recent blood flow velocity data V.sub.current based on the persistence coefficient set by the persistence coefficient determination unit 106 and two determination results from the aliasing determination unit 105. The persistence calculation is an operation that yields a residual image by performing a weighted addition on images of the current and of a plurality of consecutive frames and displaying the results.

(31) Details regarding the processing performed by the aliasing determination unit 105, the persistence coefficient determination unit 106, and the persistence calculation unit 108 in the residual image processing unit 115 will be described later.

(32) The CFM DSC unit 110 converts the coordinates of the blood flow velocity data output from the persistence calculation unit 108 and outputs it to the image combining unit 113.

(33) The tomographic image signal processing unit 111 performs dynamic filtering on the received signal to remove unnecessary noise and then performs envelope detection processing and dynamic range compression processing, and outputs tomographic image frame data to the tomographic image DSC unit 112. The tomographic image DSC unit 112 converts the coordinates of the tomographic image frame data from the tomographic image signal processing unit 111, and outputs it to the image combining unit 113.

(34) The image combining unit 113 combines frame data outputted from the CFM DSC unit 110 and the tomographic image DSC unit 112 in a way such that the frame data are combined for each pixel or for each corresponding measurement point to generate combined image frame data. Specifically, when blood flow velocity is zero, tomographic image frame data is displayed. Otherwise, two pieces of data for each pixel or data for each corresponding measurement point are combined such that the CFM frame data is displayed. It converts the data into color information according to the velocity and direction of blood flow and outputs it to the display 114. The display 114 displays data received from the image combining unit 113.

Processing of the Aliasing Determination Unit

(35) Next, the aliasing determination method in the aliasing determination unit 105 and the determination method of the persistence coefficient in the persistence coefficient determination unit 106 will be described in more detail.

(36) The aliasing determination unit 105 performs aliasing determination on the blood flow velocity data of the most recent CFM frame data V.sub.current, which is read from the frame memory unit 104, and the residual image processed blood flow velocity data of CFM frame data for the frame immediately preceding the most recent frame V.sub.out-1, which is read from the persistence memory unit 109. The aliasing determination unit 105 uses the read values of V.sub.current and V.sub.out-1, and determines whether aliasing has occurred in V.sub.current or V.sub.out-1.

(37) FIG. 2 is a schematic diagram showing the velocity and direction of blood flow obtained based on the Doppler effect. FIG. 2 is a diagram showing a complex plane. Blood flow velocity can be expressed as a vector having a phase angle on the complex plane. That is, blood flow velocity increasing or decreasing corresponds to a vector on the complex plane, shown in FIG. 2, rotating around the origin.

(38) Vectors V.sub.current and V.sub.out-1 shown in FIG. 2 are the blood flow velocity data V.sub.current and V.sub.out-1, of a certain frame, being generated by the CFM signal processing unit 103 and included in the CFM frame data, which are represented on a complex plane. As described above, in the complex plane shown in FIG. 2, the magnitude of the blood flow velocity data V.sub.current and V.sub.out-1 corresponds to the rotation angle from the +X-axis (i.e., phase angle).

(39) In the complex plane shown in FIG. 2, when the blood flow velocity is a positive value, the vector corresponding to the blood flow velocity is located in the first or second quadrant, and when the blood flow velocity is a negative value, the vector corresponding to the blood flow velocity is located in the third or fourth quadrant. Blood flow velocity with a positive value means that the blood flow is moving in a direction approaching the probe 101, whereas blood flow velocity with a negative value means that the blood flow is moving in a direction away from the probe 101.

(40) In the complex plane, when the magnitude of the velocity vector corresponding to a certain blood flow velocity gradually increases from 0, the vector corresponding to the blood flow velocity data moves in the direction in which the rotation angle from the +X axis increases. Specifically, it moves in order from the first quadrant to the second quadrant, in the direction from the +X axis to the −X axis. When the blood flow velocity exceeds a predetermined value +V.sub.max or −V.sub.max (velocity corresponding to the frequency variation corresponding to ±PRF/2), crosses the −X axis and moves from the second quadrant to the third quadrant, aliasing occurs. In this case, the blood flow velocity, which should be a positive value, is measured as a negative value. When the magnitude of the velocity vector gradually decreases from 0, the vector corresponding to the blood flow velocity data moves in the direction opposite to the above.

(41) Based on the above, the aliasing determination unit 105 determines whether or not aliasing has occurred in V.sub.current and V.sub.out-1 by comparing the magnitude of the blood flow velocity data V.sub.current of the most recent frame as described below and the magnitude of the blood flow velocity data V.sub.out-1 of the frame immediately preceding the most recent frame to a predetermined threshold Vth.

(42) Here, as a method of aliasing determination, the method disclosed in Japanese Patent No. 5652395, the aliasing determination method adopted by the present invention will be described after the items that need to be improved (items that need improvement) in the method have been cited.

Existing Aliasing Determination Method (Method Disclosed in Japanese Patent No. 5652395)

(43) Table 1 below shows the aliasing determination results obtained by using the existing aliasing method. In Table 1, the threshold value V.sub.th is the maximum value of the change in blood flow velocity assumed in the time interval between adjacent frames.

(44) TABLE-US-00001 TABLE 1 Conditions Aliasing (0) V.sub.out−1 > 0 and V.sub.current < −V.sub.th Confirmed (1) V.sub.out−1 < 0 and V.sub.current > V.sub.th Confirmed (2) V.sub.out−1 < −V.sub.th and V.sub.current > 0 Confirmed (3) V.sub.out−1 > V.sub.th and V.sub.current < 0 Confirmed (4) Other than (0) to (3) Not confirmed

(45) As shown in condition (0) of Table 1, when the blood flow velocity data of the frame immediately preceding the most recent frame V.sub.out-1 is positive and the blood flow velocity data of the most recent frame V.sub.current is smaller than −V.sub.th, assuming that aliasing has not occurred, it is unnatural if the blood flow that was in the positive direction (direction approaching the probe 101) in the frame immediately preceding the most recent frame suddenly changes to high-velocity blood flow in the reversed negative direction (direction away from the probe 101). That is, it is presumed that the velocity of blood flow that was positive in the frame immediately preceding the most recent frame increases in the most recent frame, and the blood flow velocity corresponds to +PRF/2, that is, a value greater than the blood flow velocity in which aliasing has occurred. Therefore, when condition (0) is satisfied, aliasing is determined to have occurred in V.sub.current. Condition (1) is a case where the sign of condition (0) is reversed. In the following, the blood flow velocity at which aliasing occurs is specified as V.sub.max or −V.sub.max.

(46) As shown in condition (2), when V.sub.out-1 is smaller than −V.sub.th, and V.sub.current becomes a positive value, assuming that aliasing has not occurred, it is unnatural if the high-velocity blood flow that was in the negative direction in the frame immediately preceding the most recent frame suddenly changes to blood flow in the reversed positive direction in the most recent frame. That is, it is presumed that with aliasing having occurred in V.sub.out-1 and having a value larger than the +PRF/2, the blood flow velocity has slowed down in the most recent frame in the same positive direction. Accordingly, it is determined that aliasing has occurred. Condition (3) is a case where the sign of condition (2) is reversed.

(47) If none of conditions (0) to (3) are satisfied, it is determined that aliasing has not occurred.

Items That Need Improvement in Existing Aliasing Determination Method

(48) However, in the existing aliasing determination method described above, there are improvements that need to be made as described below. The first item that needs improvement is that, even if aliasing is determined to have occurred, it is not specified which of V.sub.current or V.sub.out-1 aliasing has occurred in.

(49) Specifically, for example, when V.sub.out-1>V.sub.th and V.sub.current<−V.sub.th, thereby fulfilling condition (0) or condition (3) in Table 1 and aliasing is determined to have occurred, it may not be possible to specify which of V.sub.current or V.sub.out-1 aliasing occurred in. That is, if condition (0) is true, aliasing has occurred in V current, and if condition (3) is true, aliasing has occurred in V.sub.out-1. In such a case, it is not possible to specify which one aliasing occurred in. Also, even when V.sub.out-1<−V.sub.th and V.sub.current<V.sub.th, corresponding to condition (1) or condition (2), there may be instances where it is not possible to specify which one is aliasing in a similar manner.

(50) The second item that needs improvement is that, in some cases, aliasing is determined not to have occurred even though aliasing has actually occurred. Specifically, for example, when V.sub.out-1<−V.sub.th and V.sub.current<−V.sub.th, thereby fulfilling condition (4) in Table 1, aliasing is determined not to have occurred, and in fact, the, even when aliasing is determined to have occurred in both V.sub.current and V.sub.out-1, V.sub.out-1<−V.sub.th and V.sub.current<−V.sub.th can be satisfied. It is the same when V.sub.out-1>V.sub.th and V.sub.current>V.sub.th.

(51) The aliasing determination method of the present invention to be described below is to improve on such existing aliasing determination methods.

Aliasing Determination Method of the Aliasing Determination Unit of the Present Invention

(52) In the present invention, the aliasing determination unit 105 performs aliasing determination based on the magnitude of blood flow velocity data of the most recent frame V.sub.current, the magnitude of blood flow velocity data of the frame immediately preceding the most recent frame V.sub.out-1, and the aliasing determination result R.sub.aliasing-1 of the frame immediately preceding the most recent frame. Table 2 below shows aliasing determination results from the aliasing determination unit 105. In Table 2, “aliasing present” means the aliasing determination result of the frame immediately preceding the most recent frame R.sub.aliasing-1 determined that aliasing has occurred, and “aliasing absent” means the aliasing determination result of the frame immediately preceding the most recent frame R.sub.aliasing-1 determined that aliasing has not occurred.

(53) TABLE-US-00002 TABLE 2 Aliasing in Aliasing in Conditions V.sub.current V.sub.out−1 (10) V.sub.out−1 < −V.sub.th and V.sub.current < −V.sub.th Confirmed Confirmed and Aliasing confirmed (11) V.sub.out−1 > V.sub.th and V.sub.current > V.sub.th and Confirmed Confirmed Aliasing confirmed (12) V.sub.out−1 > 0 and V.sub.current < −V.sub.th and Confirmed Not Aliasing not confirmed confirmed (13) V.sub.out−1 < 0 and V.sub.current > Vth and Confirmed Not Aliasing not confirmed confirmed (14) V.sub.out−1 < −V.sub.th and V.sub.current > 0 and Not Confirmed Aliasingconfirmed confirmed (15) V.sub.out−1 > V.sub.th and V.sub.current < 0 and Not Confirmed Aliasing confirmed confirmed (16) Other than (10) to (15) Not Not confirmed confirmed

(54) As shown in condition (10) of Table 2, in the case where V.sub.out-1<−V.sub.th and V.sub.current<−V.sub.th, as described above, in practice, either aliasing has occurred in both V.sub.current and V.sub.out-1, or aliasing has occurred in neither. Therefore, the aliasing determination unit 105 refers to R.sub.aliasing-1 and when the aliasing determination result for the frame immediately preceding the most recent frame R.sub.aliasing-1 is “aliasing present” as shown in condition (10), it determines that aliasing has also occurred in V.sub.current and V.sub.out-1. The reason for this is as follows. Since the aliasing determination result R.sub.aliasing-1 for the frame immediately preceding the most recent frame being “aliasing present” means that aliasing has occurred in V.sub.out-1, V.sub.out-1 is larger than +PRF/2. In this case, supposing that aliasing has not occurred in V.sub.current, it is unnatural if the extremely high-velocity blood flow with a velocity that exceeds +PRF/2 in the positive direction in the frame immediately preceding the most recent frame suddenly changes to a high-velocity blood flow in the reversed negative direction in the most recent frame. For this reason, even in the most recent frame, it is estimated to be an extremely high-velocity flow that exceeds a value corresponding to that of the positive directional +PRF/2. Shown in condition (11), it is the same when V.sub.out-1<−V.sub.th and V.sub.current<−V.sub.th.

(55) FIG. 3A is a diagram showing V.sub.current and V.sub.out-1 of condition (10) on the same complex plane shown in FIG. 2. FIG. 3B is a diagram showing V.sub.current and V.sub.out-1 of condition (11) on the same complex plane shown in FIG. 2.

(56) As shown in condition (12), in the case where V.sub.out-1>0 and V.sub.current<−V.sub.th, and the aliasing determination result R.sub.aliasing-1 is “aliasing absent” because aliasing has not occurred in V.sub.out-1, for the same reason as in condition (0) of Table 1, the aliasing determination unit 105 determines that aliasing has not occurred in V.sub.out-1 and has occurred in V.sub.current. Condition (13) is a case where the sign of condition (12) is reversed.

(57) FIG. 3C is a diagram showing V.sub.current and V.sub.out-1 of condition (12) on the same complex plane shown in FIG. 2. FIG. 3D is a diagram showing V.sub.current and V.sub.out-1 of condition (13) on the same complex plane shown in FIG. 2.

(58) As shown in condition (14), in the case where V.sub.out-1 is smaller than −V.sub.th, V.sub.current becomes positive, and in the case where the aliasing determination result R.sub.aliasing-1 is “aliasing present”, aliasing has occurred in V.sub.out-1. In this case, for the same reason as in condition (2) of Table 1, the aliasing determination unit 105 determines that aliasing has occurred in V.sub.out-1 and has not occurred in V.sub.current. Condition (15) is a case where the sign of condition (14) is reversed.

(59) FIG. 3E is a diagram showing V.sub.current and V.sub.out-1 of condition (14) on the same complex plane shown in FIG. 2. FIG. 3F is a diagram showing V.sub.current and V.sub.out-1 of condition (15) on the same complex plane shown in FIG. 2.

(60) Here, since, in the case where V.sub.out-1>V.sub.th and V.sub.current<−V.sub.th, if the aliasing determination result R.sub.aliasing-1 is “aliasing absent”, condition (12) is fulfilled, and if the aliasing determination result R.sub.aliasing-1 is “aliasing present”, condition (15) is fulfilled, the presence or absence of aliasing can be determined for each value of V.sub.current and V.sub.out-1. Also, even when V.sub.out-1<−V.sub.th and V.sub.current>V.sub.th, similarly fulfilling condition (13) or condition (14), the presence or absence of aliasing in V.sub.current and V.sub.out-1 can be specified.

(61) If none of conditions (10) to (15) are satisfied, the aliasing determination unit 105 determines that aliasing has not occurred (condition (16)).

(62) Thus, even if whether or not aliasing has occurred can not be accurately determined with the existing aliasing determination method, it is possible to accurately determine whether or not aliasing has occurred with the aliasing determination method of the aliasing determination unit 105 of the present invention.

Processing by the Persistence Coefficient Determination Unit

(63) Next, the method by which the persistence coefficient determination unit 106 sets persistence numbers to the persistence calculation unit 108 will be described.

(64) The persistence coefficient determination unit 106 creates a reference index for the persistence coefficient reference memory unit 107 based on the aliasing determination results, outputted from the aliasing determination unit 105, of the blood flow velocity data V.sub.current for the most recent frame and blood flow velocity data V.sub.out-1 for the frame immediately preceding the most recent frame, and the absolute value of V.sub.current read from the frame memory unit 104. Table 3 below shows the reference index created by the persistence coefficient determination unit 106.

(65) TABLE-US-00003 TABLE 3 Aliasing in Aliasing in Reference Condition V.sub.current V.sub.out−1 Indexes (Idx) (10) Confirmed Confirmed V.sub.max (11) Confirmed Confirmed V.sub.max (12) Confirmed Not confirmed V.sub.max (13) Confirmed Not confirmed V.sub.max (14) Not confirmed Confirmed Abs (V.sub.current) (15) Not confirmed Confirmed Abs (V.sub.current) (16) Not confirmed Not confirmed Abs (V.sub.current)

(66) When aliasing has occurred in V.sub.current, it is possible that the blood flow velocity V.sub.current has a magnitude such that it exceeds the absolute value of V.sub.max or −V.sub.max. For this reason, the reference index becomes V.sub.max. Otherwise, it becomes the absolute value of V.sub.current, Abs(V.sub.current).

(67) In the persistence coefficient reference memory unit 107, a reference table composed of persistence coefficients mapped to the reference index is stored. The persistence coefficient determination unit 106 accesses the persistence coefficient reference memory unit 107, reads the persistence coefficients mapped to the created reference index, and outputs them to the persistence calculation unit 108.

(68) FIG. 4 is a graph showing an example of the relationship between the reference index and the persistence coefficient. In FIG. 4, the horizontal axis represents the reference index, and the vertical axis shows the persistence coefficient. As shown in Table 3, the reference index is V.sub.max or the absolute value of V.sub.current, Abs(V.sub.current). When the absolute value of V.sub.current is lower than the threshold V.sub.th, the persistence coefficient C.sub.persistence, which decreases steadily with increasing V.sub.current, is mapped to V.sub.current. That is, when the absolute value of V.sub.current is lower than the threshold V.sub.th, a persistence coefficient C.sub.persistence, which is proportional and different from the blood flow velocity of the most recent frame V.sub.current, is mapped to V.sub.current.

(69) As described above, the persistence calculation is a calculation that yields a residual image effect by performing weighted addition on images of the current and of a plurality of consecutive frames and utilizing the display of the results. Details regarding the persistence calculation will be described later.

Processing by the Persistence Calculation Unit 108

(70) The persistence calculation unit 108 performs the persistence calculation on the most recent blood flow velocity data V.sub.current based on the persistence coefficient set by the persistence coefficient determination unit 106 and two determination results obtained from the aliasing determination unit 105.

(71) When the determination result from the aliasing determination unit 105 for V.sub.current is true, i.e., when aliasing is determined to have occurred in V.sub.current, the persistence calculation unit 108, before performing the persistence calculation, performs the following correction to V.sub.current. Specifically, when V.sub.current<0, the persistence calculation unit 108 corrects V.sub.current by adding the repetition frequency PRF to V.sub.current (V.sub.current+PRF), and when V.sub.current>0, V.sub.current is corrected by subtracting PRF from V.sub.current (V.sub.current−PRF).

(72) Similarly, when the determination result from the aliasing determination unit 105 for V.sub.out-1 is true, i.e., when aliasing is determined to have occurred in V.sub.out-1, the persistence calculation unit 108, before performing the persistence calculation, performs a correction to V.sub.out-1 similar to that described above for V.sub.current. Specifically, when V.sub.out-1<0, the persistence calculation unit 108 corrects V.sub.out-1 to V.sub.out-1+PRF. When V.sub.out-1>0, the persistence calculation unit 108 corrects V.sub.out-1 to V.sub.out-1+PRF.

(73) The persistence calculation unit 108, using the blood flow velocity data V.sub.current and V.sub.out-1, performs the persistence calculation with the following equation (1).
V.sub.out=(1−C.sub.persistence)×V.sub.current+C.sub.persistence×V.sub.out-1.  (1)

(74) In equation (1), V.sub.out is the blood flow velocity data that was subjected to the persistence calculation (residual image processing) in the most recent frame, and C.sub.persistence is the persistence coefficient (0≤C persistence≤1) set by the persistence coefficient determination unit 106.

(75) When upon performing the persistence calculation using equation (1) as described above, V.sub.out>+PRF/2, the persistence calculation unit 108 corrects V.sub.out to V.sub.out−PRF, and when V.sub.out is <−PRF/2, V.sub.out is corrected to V.sub.out+PRF. The persistence calculation unit 108 outputs the corrected V.sub.out to the CFM DSC unit 110 and the persistence memory unit 109.

(76) By such processing, even if aliasing has occurred in V.sub.current, V.sub.out-1 or V.sub.out, V.sub.out that has the effect of aliasing canceled can be calculated.

(77) The persistence calculation by the persistence calculation unit 108 described above is performed for each pixel or each measurement point of one frame of blood flow velocity data (CFM frame data).

(78) Thus, by using a persistence coefficient set by the persistence coefficient determination unit 106, the persistence calculation unit 108 can obtain the following results by performing the persistence calculation.

(79) That is, in the case where the magnitude of blood flow velocity V.sub.current of the most recent frame is small, because the persistence coefficient C.sub.persistence set by the persistence coefficient determination unit 106 increases, the weight of the blood flow velocity V.sub.out-1 of the immediately preceding frame increases when the persistence calculation unit 108 performs the persistence calculation such as in equation (1). As a result, when V.sub.current when is small, V.sub.out, that strongly reflects V.sub.out-1 and is the persistence calculation result of the persistence calculation, is calculated. Because of this, using V.sub.out, the change in the color flow mapping image generated by the CFM DSC unit 110 and the image combining unit 113 becomes smooth and black spots are less likely to occur.

(80) Also, in the case where V.sub.current is large, because the persistence coefficient C.sub.persistence set by the persistence coefficient determination unit 106 decreases, the weight of V.sub.out-1 decreases when the persistence calculation unit 108 performs the persistence calculation such as in equation (1). As a result, in the case where V.sub.current is large, with the influence of V.sub.out-1 being small, V.sub.out is calculated, and it is possible to realize a color flow mapping image reflecting the sharp increase in blood flow velocity in real time.

(81) Moreover, the persistence coefficient determination unit 106 monotonically reduces the persistence coefficient C.sub.persistence with increases in persistence. Therefore, in the case where the blood flow velocity increases with time, the persistence coefficient determination unit 106 decreases the persistence coefficient C.sub.persistence with time, so the residual image effect from the persistence calculation of the persistence calculation unit 108 decreases and the color flow mapping image changes abruptly. Conversely, in the case where the blood flow velocity decreases with time, the persistence coefficient determination unit 106 increases the persistence coefficient C.sub.persistence with time, so the residual image effect from the persistence calculation of the persistence calculation unit 108 increases and the color flow mapping image changes slowly.

(82) Back to the description of the processing of persistence calculation unit 108. The persistence calculation unit 108, using the persistence calculation result V.sub.out, determines whether or not aliasing has occurred in the last frame in blood flow velocity data V.sub.out that had been subjected to residual image processing after performing the persistence calculation using equation (1). Here, the persistence calculation unit 108 determines that aliasing has occurred in V.sub.out when V.sub.out>+PRF/2 or V.sub.out<−PRF/2; otherwise (−PRF/2≤V.sub.out≤PRF/2), it determines that aliasing has not occurred.

(83) The persistence calculation unit 108 then outputs the aliasing determination result R.sub.aliasing of the blood flow velocity data V.sub.out of the most recent frame to the persistence memory unit 109. The aliasing determination result R.sub.aliasing that is outputted from the persistence calculation unit 108 is stored in the persistence memory unit 109, and at the time of processing in the next frame, the aliasing determination result R.sub.aliasing-1 of the frame immediately preceding the most recent frame is outputted to the aliasing determination unit 105.

Operation and Results

(84) As described above, in the ultrasonic diagnostic apparatus 11 according to the first embodiment of the present invention, the aliasing determination unit 105 adds the blood flow velocity V.sub.current of the most recent frame and the blood flow velocity V.sub.out-1 of the frame preceding the most recent frame, and based on the aliasing determination result R.sub.aliasing-1 of the blood flow velocity data that was determined for the frame preceding the most recent frame, with the frame having been subjected to residual image processing, performs the aliasing determination. Because of this, whether or not aliasing has occurred in the blood flow velocity data V.sub.current of the most recent frame and the blood flow velocity data V.sub.out-1 of the frame immediately preceding the most recent frame can be accurately determined.

(85) Further, in the ultrasonic diagnostic apparatus 11 according to the first embodiment of the present invention, based on the result of the aliasing determination of such an aliasing determination method, the persistence calculation unit 108 is performs the persistence calculation using V.sub.current and V.sub.out-1 and outputs the result. Therefore, it is possible to generate and display a tomographic image that accurately reflects blood flow velocity.

(86) Furthermore, in the ultrasonic diagnostic apparatus 11 according to the first embodiment of the present invention, based on the blood flow velocity data V.sub.out of the most recent frame and has been subjected to residual image processing, which is the result from the persistence calculation, the persistence calculation unit 108 determines whether or not aliasing has occurred in V.sub.out and outputs that result R.sub.aliasing to the persistence memory unit 109. Then, in the processing of the next frame, the persistence memory unit 109 outputs V.sub.out inputted in the processing of the previous frame as V.sub.out-1 of the next frame, and outputs R.sub.aliasing inputted in the processing of the previous frame as R.sub.aliasing-1.

(87) Then, in the ultrasonic diagnostic apparatus 11 according to a first embodiment of the present invention, when aliasing is determined to have occurred in the blood flow velocity data V.sub.current for the most recent frame or the blood flow velocity data V.sub.out-1 for the frame immediately preceding the most recent frame, the persistence calculation unit 108 performs corrections on the blood flow velocity data in which aliasing occurred. Furthermore, the persistence calculation unit 108 determines whether or not aliasing has occurred for the most recent frame in the blood flow velocity data V.sub.out that has been subjected to residual image processing, and performs a correction to V.sub.out when aliasing is determined to have occurred. Therefore, even if aliasing has occurred in V.sub.current, V.sub.out-1 or V.sub.out, V.sub.out that has the effect of aliasing canceled can be calculated.

Second Embodiment

(88) The following describes the ultrasonic diagnostic apparatus 12 according to the second embodiment of the present invention. FIG. 5 is a block diagram showing an example configuration of an ultrasonic diagnostic apparatus according to the second embodiment of the present invention.

(89) The ultrasonic diagnostic apparatus 12 shown in FIG. 5 is different from the first embodiment of the ultrasonic diagnostic apparatus 11 described above in that the residual image processing unit has two persistence coefficient determination units, two persistence coefficient reference memory units and two persistence calculation units. In the following, a description will be given only for the differences from the first embodiment.

(90) As shown in FIG. 5, the residual image processing unit 115′ includes a frame memory unit 104, an aliasing determination unit 105, a first persistence coefficient determination unit 106A, a second persistence coefficient determination unit 106B, a first persistence coefficient reference memory unit 107A, a second persistence coefficient reference memory unit 107B, a first persistence calculation unit 108A, a second persistence calculation unit 108B, and a persistence memory unit 109.

(91) The first persistence coefficient determination unit 106A creates the first reference index for the first persistence coefficient reference memory unit 107A based on the two aliasing determination results outputted from the aliasing determination unit 105 and the absolute value of V.sub.current read from the frame memory unit 104. Information regarding the first reference index will be described later.

(92) The first persistence coefficient determination unit 106A accesses the first persistence coefficient reference memory unit 107A, reads out the first persistence coefficient mapped to the first reference index, and sets them to the first persistence calculation unit 108A. In the first persistence coefficient reference memory unit 107A, a first reference table containing first persistence coefficients mapped to the magnitude of blood flow velocity data is stored in advance. The first reference table includes two or more different persistence coefficients that are mapped to the magnitude of blood flow velocity data.

(93) In contrast, the second persistence coefficient determination unit 106B creates a second reference index for the second persistence coefficient reference memory unit 107B based on the two aliasing determination results outputted from the aliasing determination unit 105 and the absolute value of V.sub.out-1 read from the persistence memory unit 109. Creation of the second reference index in the second persistence coefficient determination unit 106B will be described later.

(94) The second persistence coefficient determination unit 106B accesses the second persistence coefficient reference memory unit 107B, reads the second persistence coefficients mapped to the second reference index, and sets them in the second persistence calculation unit 108B. The second persistence coefficient reference memory unit 107B, a second reference table containing second persistence coefficients mapped to the magnitude of blood flow velocity data is stored in advance. Although the second reference table also includes persistence coefficients of two or more different values that are mapped according to the values of blood flow velocity, as explained in detail below, the values are different from the first persistence coefficient and second persistence coefficient that are mapped to the same blood flow velocity value.

(95) Similar to the persistence calculation unit 108 of the first embodiment, after correcting for aliasing, the first persistence calculation unit 108A performs the persistence calculation on V.sub.current and V.sub.out-1 using equation (1) based on the persistence coefficient set by the first persistence coefficient determination unit 106A and the two determination results from the aliasing determination unit 105.

(96) Similar to the persistence calculation unit 108 of the first embodiment, after correcting for aliasing, the second persistence calculation unit 108B performs the persistence calculation on V.sub.current and V.sub.out-1 using equation (1) based on the persistence coefficient set by the second persistence coefficient determination unit 106B and the two determination results from the aliasing determination unit 105.

(97) Calculations of the first persistence calculation unit 108A and second persistence calculation unit 108B are the same except for their differing set persistence coefficients.

(98) The first persistence calculation unit 108A and the second persistence calculation unit 108B output their respective results for the persistence calculation to the maximum value selection unit 116.

(99) Based on the persistence calculation results inputted from the first persistence calculation unit 108A and the second persistence calculation unit 108B, for each pixel or for each value of the corresponding measurement point, the maximum value selection unit 116 compares the absolute value of the blood flow velocity data which is the result of the persistence calculation, the larger value is selected, and blood flow velocity data V.sub.out subjected to residual image processing is generated for the most recent frame.

(100) The maximum value selection unit 116 corrects V.sub.out for aliasing, as performed by the persistence calculation unit 108 of the first embodiment.

(101) Then, the maximum value selection unit 116 outputs the corrected V.sub.out to the CFM DSC unit 110 and the persistence memory unit 109.

(102) Furthermore, the maximum value selection unit 116, as performed by the persistence calculation unit 108 of the first embodiment, determines whether or not aliasing has occurred in the most recent frame in the blood flow velocity data V.sub.out that are subjected to residual image processing. The maximum value selection unit 116 then outputs the aliasing determination result R.sub.aliasing of the blood flow velocity data V.sub.out in the most recent frame to the persistence memory unit 109.

(103) Next, a method of determining the first and second persistence coefficients will be described in detail. As described above, the first persistence coefficient determination unit 106A creates the first reference index for persistence coefficient reference memory unit 107A based on the two folded determination results outputted from the aliasing determination unit 105 and the absolute value of V.sub.current read from the frame memory unit 104.

(104) The first reference index created in Table 4 is shown.

(105) TABLE-US-00004 TABLE 4 Reference Reference Indexes Indexes in First in Second Persistence Persistence Coefficient Coefficient Determininion Determininion Aliasing in Aliasing in Unit Unit Condition V.sub.current V.sub.out−1 (Idx1) (Idx2) (10) Confirmed Confirmed V.sub.max V.sub.max (11) Confirmed Confirmed V.sub.max V.sub.max (12) Confirmed Not V.sub.max Abs (V.sub.out−1) confirmed (13) Confirmed Not V.sub.max Abs (V.sub.out−1) confirmed (14) Not Confirmed Abs (V.sub.current) V.sub.max confirmed (15) Not Confirmed Abs (V.sub.current) V.sub.max confirmed (16) Not Not Abs (V.sub.current) Abs (V.sub.out−1) confirmed confirmed

(106) When aliasing has occurred in V.sub.current, it is possible that V.sub.current exceeds the absolute value of V.sub.max or −V.sub.max. Therefore, when aliasing has occurred in V.sub.current, the first reference index becomes V.sub.max. Otherwise, it becomes the absolute value of V.sub.current, Abs(V.sub.current).

(107) FIG. 6A is a graph showing an example of the relationship between the first reference index and the first persistence coefficient. In FIG. 6A, the horizontal axis represents the first reference index and the vertical axis shows the first persistence coefficient. As shown in Table 4, the first reference index is V.sub.max or the absolute value of V.sub.current, Abs(V.sub.current). When the absolute value of V.sub.current is lower than the threshold V.sub.th, the first persistence coefficient C.sub.persistence, which decreases steadily with increasing V.sub.current, is mapped to V.sub.current. That is, when the absolute value of V.sub.current is lower than the threshold V.sub.th, a persistence coefficient C.sub.persistence, which is proportional and different from the blood flow velocity of the most recent frame V.sub.current, is mapped to V.sub.current.

(108) On the other hand, the second persistence coefficient determination unit 106B creates a second reference index for the persistence coefficient reference memory unit 107B based on the two aliasing determination results outputted by the aliasing determination unit 105 and the absolute value of V.sub.out-1 outputted by the persistence memory unit 109. Refer to Table 4 for the second reference index.

(109) When aliasing has occurred in V.sub.out-1, it is possible that V.sub.out-1 exceeds the absolute value of V.sub.max or −V.sub.max. Therefore, when aliasing has occurred in V.sub.out-1, the second reference index becomes V.sub.max. Otherwise, it becomes the absolute value of V.sub.out-1, Abs(V.sub.out-1).

(110) FIG. 6B is a graph showing an example of the correspondence relationship between the second reference index and the second persistence coefficient. In FIG. 6B, the horizontal axis represents the second reference index and the vertical axis shows the second persistence coefficient. As shown in Table 4, the second reference index is V.sub.max or the absolute value of V.sub.out-1, Abs(V.sub.out-1). When the absolute value of V.sub.out-1 is lower than the threshold V.sub.th, the second persistence coefficient C.sub.persistence, which decreases steadily with increasing V.sub.current, is mapped to V.sub.current. That is, when the absolute value of V.sub.out-1 is lower than the threshold V.sub.th, a persistence coefficient C.sub.persistence, which is proportional and different from the blood flow velocity of the frame before the lastest frame V.sub.out-1, is mapped to V.sub.out-1.

(111) As shown in FIGS. 6A and 6B, the first persistence coefficient is mapped to the blood flow velocity of the most recent frame V.sub.current, and the value of the first persistence coefficient is as small as the magnitude of the absolute value of V.sub.current. The greater the first persistence factor, since the operation more considering the blood flow velocity of the previous frame, the first persistence calculation unit 108A reduces residual image effect, quickly changing the blood flow velocity performing an operation to be. In contrast, the second persistence coefficient, blood flow velocity V of the previous frame.sub.out-1 mapped to, and, V.sub.out-1 for a larger value as the absolute value is large, the second persistence calculation unit 108B increases the residual image effect, reduce operation changes in blood flow velocity.

(112) The first persistence calculation unit 108A as described above, to reduce the residual image effect, quickly for performing the operation of changing the blood flow velocity, monotonously the first persistence coefficient as the first reference index is increased reduced, it can realize a color flow mapping image reflecting the sharp increase in blood flow velocity in real time as the blood flow velocity becomes higher.

(113) When the second persistence calculation unit 108B adds color to the color tone or gradation in the blood flow image as the blood flow velocity increases to display an image with an enhanced residual image effect, a dark display will be shown as a residual image for an unnecessarily long time when the blood flow velocity is low. For example, when the probe is moved, the blood flow looks as if it had a tail. Because of this, the second persistence coefficient is increased monotonically as the second reference index is increased, and it becomes possible to suppress the residual image effect as the blood flow velocity decreases.

Third Embodiment

(114) The following describes the ultrasonic diagnostic apparatus 13 according to the third embodiment of the present invention. FIG. 7 is a block diagram showing an example configuration of an ultrasonic diagnostic apparatus according to the third embodiment of the present invention.

(115) The ultrasonic diagnostic apparatus 13 shown in FIG. 7 is different from the first embodiment of the ultrasonic diagnostic apparatus 11 described above in that the residual image processing unit has two persistence coefficient determination units and two persistence coefficient reference memory units, and the two persistence coefficient determination units have a persistence coefficient selection unit that selects from the set persistence coefficients. In the following, a description will be given only of the differences from the first and second embodiments.

(116) As shown in FIG. 7, the residual image processing unit 115″ includes a frame memory unit 104, an aliasing determination unit 105, a first persistence coefficient determination unit 106A, a second persistence coefficient determination unit 106B, a persistence coefficient selection unit 117, a first persistence coefficient reference memory unit 107A, a second persistence coefficient reference memory unit 107B, a persistence calculation unit 108, and a persistence memory unit 109.

(117) The persistence coefficient selection unit 117 reads the blood flow velocity data V.sub.current for the most recent frame from the frame memory unit 104 and the blood flow velocity data V.sub.out-1 for the frame preceding the most recent frame from the persistence memory unit 109, compares the absolute values of V.sub.current and V.sub.out-1, and based on the comparison results, selects one of either the first persistence coefficient set by the first persistence coefficient determination unit 106A or the second persistence coefficient set by the second persistence coefficient determination unit 106B to set in the persistence calculation unit 108.

(118) Specifically, the persistence coefficient selection unit 117, in the case where the absolute value of V.sub.current is larger than the absolute value of V.sub.out-1, selects the first persistence coefficient; otherwise, it selects the second persistence.

(119) Therefore, as in the second embodiment, the first persistence coefficient is mapped to the blood flow velocity V.sub.current of the most recent frame and is as small as the magnitude of the absolute value of V.sub.current. The greater the first persistence factor, since the operation more considering the blood flow velocity of the previous frame, persistence calculation unit 108, in the case of using the first persistence coefficient, reduces residual image effect, it is possible to perform an operation to quickly change the blood flow velocity. In contrast, the second persistence coefficient, blood flow velocity V of the previous frame.sub.out-1 mapped to, and, V.sub.out-1 for a larger value as the absolute value is large, the second persistence calculation unit 108 in the case of using the stance factor enhances the residual image effect, it is possible to reduce operation changes in blood flow velocity.

(120) Having described various embodiments with reference to the accompanying drawings, the present invention is not limited to these examples. Within the scope of the claims, various modifications or amendments conceivable by those skilled in the art are also included in the technical scope of the present invention. Further, without departing from the scope of the disclosure, it may be any combination of the components in the above embodiment.

(121) In the above embodiment, although processing was performed using the blood flow velocity data included in CFM frame data, in the case where CFM frame data includes data not other than blood flow velocity data V.sub.current or V.sub.out-1 (such as blood flow power data or dispersion data of the blood flow velocity), the persistence calculation unit 108, even with respect to data other than blood flow velocity data, uses the blood flow velocity data for the most recent frame and blood flow velocity data for the frame preceding the most recent frame and the persistence coefficient C.sub.persistence and performs a persistence calculation in a similar way as in the above equation (1), and may obtain residual image processed data.

(122) Further, in the above embodiment, it has been performed persistence processing using the blood flow velocity data of the most recent frame and a previous frame, two front, or three or more previous frames of the blood flow velocity data it may be subjected to a persistence processing also be used. Further, not being limited to equation (1), persistence processing may be performed using other calculation equations.

(123) Each functional block employed in the description of the above embodiments are typically implemented as an LSI configured by an integrated circuit having an input terminal and an output terminal. These may be implemented individually as single chips, or may be integrated into one chip including part or all. Although it is referred to here as the LSI, depending on the degree of integration, it may also be called an IC (Integrated Circuit), system LSI, or super LSI.

(124) Further, the method of circuit integration is not limited to the LSI, and may be realized using a dedicated circuit or a general-purpose processor. After manufacturing of the LSI, a programmable FPGA (Field Programmable Gate Array), Reconfigurable Processor (Reconfigurable Processor) connections or settings of circuit cells in the LSI may be utilized.

(125) Furthermore, if, by way of advancements in semiconductor technology or a derivative, an integrated circuit technology comes out to replace LSI's, as a matter of course, it may be used to integrate functional blocks.

(126) The present invention is suitable for a ultrasonic diagnostic apparatus capable of displaying the blood flow velocity.

(127) According to an embodiment, there is provided an ultrasonic diagnostic apparatus which is capable of accurately determining the occurrence of aliasing, and when the aliasing occurs, is capable of adequately correcting for the effect.

(128) Although embodiments of the present invention have been described and illustrated in detail, it is clearly understood that the same is by way of illustration and example only and not limitation, the scope of the present invention should be interpreted by terms of the appended claims.