Abstract:
Disclosed is an ultrasound color Doppler image system. The ultrasound color Doppler image system includes a calculation unit that calculates a mean value associated with an I/Q signal corresponding to a pixel of a color image in an ultrasound image and generates a multiplication value using the calculated mean value; a comparison unit that compares the generated multiplication value with the mean value; and a masking unit that performs masking of the pixel based on a comparison result, and, the calculation unit calculates the mean value of the I/Q signal for each frame, selects a reference mean value based on scales of the calculated mean values, and generates the multiplication value by multiplying the selected reference mean value and a scale factor.

Description:
CROSS-REFERENCE TO RELATED APPLICATION 
       [0001]    This application claims the benefit of Korean Patent Application No. 10-2009-0118439, filed on Dec. 2, 2009, in the Korean Intellectual Property Office, the disclosure of which is incorporated herein by reference. 
       BACKGROUND 
       [0002]    1. Field 
         [0003]    The present invention relates to an ultrasound color Doppler image system, and more specifically, to an ultrasound color Doppler image system that filters a clutter signal from a color Doppler image. 
         [0004]    2. Description of the Related Art 
         [0005]    An ultrasound color Doppler image system may emit an ultrasound to a human body and may measure a Doppler shift frequency of the ultrasound reflected from a blood flow, thereby detecting a distribution of the blood flow in a real time. 
         [0006]    Although the emitted ultrasound is well focused on the blood flow, a part of the ultrasound may be emitted in an undesired direction, and thus, a reflection signal reflected from the blood flow may be returned together with an undesired signal reflected from other tissue besides the blood flow. In this instance, the reflection signal from the blood flow is referred to as a Doppler signal, and the undesired signal reflected from other tissue is referred to as a clutter signal. 
         [0007]    Usually, when the ultrasound is well focused, most of ultrasound energy may be transmitted to a focal point, and a relatively small portion of the ultrasound energy may leak to the outside, and thus, a desired signal from the focal point may be mostly returned. However a reflectivity of blood flow is relatively small compared with a reflectivity of peripheral tissue, such as a vessel wall, a muscle, and the like, and thus, in general, the clutter signal is greater than the Doppler signal even though a relatively small amount of ultrasound energy generally leaks. 
         [0008]    When there is little movement of the peripheral tissue, a Doppler shift frequency is relatively close to zero, and thus, it is generally well-known that the clutter signal is easily eliminated by a high-pass clutter filter. Also, when movement of the peripheral tissue becomes faster, the Doppler shift frequency of the clutter signal may increase. Accordingly, the clutter signal may not be completely eliminated with a general clutter filter, and thus, an aberrant color may appear on a screen. The fast movement of the tissue may be displayed as though the fast movement is the blood flow. The aberrant color is an artifact displayed on the screen, and the artifact may appear as an instantaneous flash on the screen, thereby being referred to as a flash artifact. 
       SUMMARY 
       [0009]    An aspect of the present invention provides an apparatus and a method that filters a clutter signal to eliminate a flash artifact component from a color Doppler image. 
         [0010]    Another aspect of the present invention provides an apparatus and a method that eliminates a clutter signal from both a wideband frequency, and a narrowband frequency. 
         [0011]    According to exemplary embodiments of the present invention, there may be provided an ultrasound color Doppler image system, the system including a calculation unit to calculate a mean value associated with an I/Q signal corresponding to a pixel of a color image in an ultrasound image, and to generate a multiplication value using the calculated mean value, a comparison unit to compare the generated multiplication value with the mean value, and a masking unit to perform masking of the pixel based on a comparison result. 
         [0012]    The calculation unit may calculate the mean value of the I/Q signal for each frame, may select a reference mean value based on scales of the calculated mean values, and may generate the multiplication value by multiplying the selected reference mean value and a scale factor. 
         [0013]    According to exemplary embodiments of the present invention, there may be provided a method of filtering a clutter signal in an ultrasound color Doppler image system, the method including calculating a mean value associated with an I/Q signal corresponding to a pixel of a color image in an ultrasound image, generating a multiplication value using the calculated mean value, comparing the generated multiplication value with the mean value, and performing a masking of the pixel based on a comparison result. 
         [0014]    The calculating may calculate the mean value of the I/Q signal for each frame. The generating of the multiplication value may select a reference mean value based on scales of calculated mean values; and may generate the multiplication value by multiplying the selected reference mean value and a scale factor. 
         [0015]    Additional aspects and/or advantages will be set forth in part in the description which follows and, in part, will be apparent from the description, or may be learned by practice of the embodiments. 
       EFFECT 
       [0016]    According to exemplary embodiments of the present invention, a clutter signal is filtered, and thus, a flash artifact component is eliminated from a color Doppler image. 
         [0017]    According to exemplary embodiments of the present invention, a clutter signal is eliminated from both a wideband frequency, and a narrowband frequency. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0018]    These and/or other aspects and advantages will become apparent and more readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings of which: 
           [0019]      FIG. 1  is a block diagram illustrating an ultrasound color Doppler image system according to an exemplary embodiment of the present invention; 
           [0020]      FIGS. 2 through 4  are diagrams illustrating an example of an phase and a magnitude of two successive signals in an ensemble I/Q signal according to an exemplary embodiment of the present invention; 
           [0021]      FIG. 5  is a diagram illustrating an example of masking a pixel in a color image according to an exemplary embodiment of the present invention; 
           [0022]      FIG. 6  is a diagram illustrating an image before eliminating a flash artifact from a color Doppler image and an image after eliminating the flash artifact from the color Doppler image according to an exemplary embodiment of the present invention; and 
           [0023]      FIG. 7  is a flowchart illustrating a method of filtering a clutter signal in an ultrasound color Doppler image system according to an exemplary embodiment of the present invention. 
       
    
    
     DETAILED DESCRIPTION 
       [0024]    Reference will now be made in detail to example embodiments, examples of which are illustrated in the accompanying drawings, wherein like reference numerals refer to the like elements throughout. An ultrasound diagnosis apparatus utilizing a touch interaction is described below to explain the present disclosure by referring to the figures. 
         [0025]      FIG. 1  illustrates an ultrasound color Doppler image system according to an exemplary embodiment of the present invention. 
         [0026]    Referring to  FIG. 1 , an ultrasound color Doppler image system  100  may include a calculation unit  110 , a comparison unit  120 , a masking unit  130 , and a controller  140 . 
         [0027]    The calculation unit  110  may calculate a mean phase and a mean magnitude of a first signal corresponding to a pixel of a color image in an ultrasound image, and may multiply, by a predetermined scale factor, the calculated mean phase and the calculated mean magnitude. Here, the color image may be an image before a clutter signal is eliminated, namely, an image including a blood flow and a peripheral tissue, such as a vessel wall, muscle, and the like. Also, the first signal may be an In-phase/Quadrature-phase signal (I/Q signal) per pixel or an I/Q signal per sample among an I/Q signal based on an ensemble unit or I/Q signal based on a frame unit. 
         [0028]    In this instance, the calculation unit  110  may calculate the mean phase and the mean magnitude of the first signal (an I/Q signal) using a second signal (another I/Q signal) that has a temporal continuity with the first signal (the I/Q signal). As an example, when the calculation unit  110  calculates the mean phase and the mean magnitude of the first signal based on Equation 1 may be as given below: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       MeanPhase 
                        
                       
                         ( 
                         
                           u 
                           , 
                           v 
                         
                         ) 
                       
                     
                     = 
                     
                       
                         ( 
                         
                           
                             ∑ 
                             
                               k 
                               = 
                               0 
                             
                             
                               EL 
                               - 
                               1 
                             
                           
                            
                           
                              
                             
                               
                                 tan 
                                 
                                   - 
                                   1 
                                 
                               
                                
                               
                                 ( 
                                 
                                   
                                     E 
                                      
                                     
                                       ( 
                                       
                                         k 
                                         , 
                                         u 
                                         , 
                                         v 
                                       
                                       ) 
                                     
                                   
                                   × 
                                   
                                     
                                       E 
                                        
                                       
                                         ( 
                                         
                                           
                                             k 
                                             + 
                                             1 
                                           
                                           , 
                                           u 
                                           , 
                                           v 
                                         
                                         ) 
                                       
                                     
                                     * 
                                   
                                 
                                 ) 
                               
                             
                              
                           
                         
                         ) 
                       
                       n 
                     
                   
                    
                   
                     
 
                   
                    
                   
                     
                       MeanMagn 
                        
                       
                         ( 
                         
                           u 
                           , 
                           v 
                         
                         ) 
                       
                     
                     = 
                     
                       
                         ( 
                         
                           
                             ∑ 
                             
                               k 
                               = 
                               0 
                             
                             
                               EL 
                               - 
                               1 
                             
                           
                            
                           
                             ( 
                             
                               
                                  
                                 
                                   E 
                                    
                                   
                                     ( 
                                     
                                       k 
                                       , 
                                       u 
                                       , 
                                       v 
                                     
                                     ) 
                                   
                                 
                                  
                               
                               × 
                               
                                  
                                 
                                   
                                     E 
                                      
                                     
                                       ( 
                                       
                                         
                                           k 
                                           + 
                                           1 
                                         
                                         , 
                                         u 
                                         , 
                                         v 
                                       
                                       ) 
                                     
                                   
                                   * 
                                 
                                  
                               
                             
                             ) 
                           
                         
                         ) 
                       
                       n 
                     
                   
                 
               
               
                 
                   [ 
                   
                     Equation 
                      
                     
                         
                     
                      
                     1 
                   
                   ] 
                 
               
             
           
         
       
     
         [0029]    Here, MeanPhase denotes the mean phase, MeanMagn denotes the mean magnitude, u denotes a scanline number, v denotes a sample number or a pixel number, k denotes an ensemble number, EL denotes a length of ensemble, * denotes a conjugated complex number, and E(k,u,v) denotes a modulated I/Q signal that is expressed in a form of a complex number (I+jQ). 
         [0030]    The calculation unit  110  may generate a first multiplication value and a second multiplication value by multiplying a predetermined scale factor and each of a maximum value for the calculated mean phase and a maximum value for the calculated mean magnitude. Here, the scale factor may be set in advance depending on a user or a system environment. 
         [0031]    Specifically, the calculation unit  110  may generate the first multiplication value by multiplying the maximum value for the mean phase, namely, a maximum mean phase, by the scale factor. Also, the calculation unit  110  may generate the second multiplication value by multiplying the maximum value for the mean magnitude, namely a maximum mean magnitude, by the scale factor. 
         [0032]    The calculation unit  110  may calculate the maximum value for the calculated mean phase and the maximum value for the calculated mean magnitude based on an example, such as Equation 2 as given below. The calculation unit  110  may compare mean phases with each other to calculate the maximum value for the mean phase, and may compare mean magnitudes with each other to calculate the maximum value for the mean magnitude. 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       MaxPhase 
                       
                         init 
                         = 
                         0 
                       
                     
                     = 
                     
                       
                         
                           Max 
                           V 
                         
                         
                           v 
                           = 
                           0 
                         
                       
                        
                       
                         ( 
                         
                           
                             
                               Max 
                               U 
                             
                             
                               u 
                               = 
                               0 
                             
                           
                            
                           
                             ( 
                             
                               
                                 MeanPhase 
                                  
                                 
                                     
                                 
                                  
                                 
                                   ( 
                                   
                                     u 
                                     , 
                                     v 
                                   
                                   ) 
                                 
                               
                               , 
                               MaxPhase 
                             
                             ) 
                           
                         
                         ) 
                       
                     
                   
                    
                   
                     
 
                   
                    
                   
                     
                       MaxMagn 
                       
                         init 
                         = 
                         0 
                       
                     
                     = 
                     
                       
                         
                           Max 
                           V 
                         
                         
                           v 
                           = 
                           0 
                         
                       
                        
                       
                         ( 
                         
                           
                             
                               Max 
                               U 
                             
                             
                               u 
                               = 
                               0 
                             
                           
                            
                           
                             ( 
                             
                               
                                 MeanMagn 
                                  
                                 
                                     
                                 
                                  
                                 
                                   ( 
                                   
                                     u 
                                     , 
                                     v 
                                   
                                   ) 
                                 
                               
                               , 
                               MaxMagn 
                             
                             ) 
                           
                         
                         ) 
                       
                     
                   
                 
               
               
                 
                   [ 
                   
                     Equation 
                      
                     
                         
                     
                      
                     2 
                   
                   ] 
                 
               
             
           
         
       
     
         [0033]    Here, MaxPhase denotes the maximum value for the mean phase, MaxMagn denotes the maximum value for the mean magnitude, u denotes a scanline number, v denotes a sample number or a pixel number, U denotes a maximum value for the scaneline, V denotes a maximum value for the sample number or the pixel number, MeanPhase denotes a mean phase, and MeanMagn denotes a mean magnitude. 
         [0034]    For reference, the ultrasound color Doppler image system  100  may perform modulation based on a pixel unit with respect to all pixels in an I/Q signal based on an ensemble unit, before performing calculation by the calculation unit  110 . Accordingly, the first signal may be modulated signal. 
         [0035]    The comparison unit  120  may compare the generated first multiplication value, the generated second multiplication value with the calculated mean magnitude and the calculated mean phase, respectively. Specifically, the comparison unit  120  may compare the generated first multiplication value with the calculated mean phase, and may compare the generated second multiplication value with the calculated mean magnitude 
         [0036]    The masking unit  130  may perform masking of the corresponding pixel of the color image based on a comparison result. 
         [0037]    As an example, when the first multiplication value is less than the mean phase and the second multiplication value is greater than the mean magnitude according to the comparison result (first multiplication value&lt;mean phase &amp; second multiplication value&gt;mean magnitude), the masking unit  130  may perform masking of the corresponding pixel of the color image. 
         [0038]    Also, when the first multiplication value is greater than the mean phase, and the second multiplication value is greater than the mean magnitude according to the comparison result (first multiplication value&gt;mean phase &amp; second multiplication value&gt;mean magnitude), the masking unit  130  may perform masking of the corresponding pixel of the color image. 
         [0039]    Also, when the first multiplication value is less than the mean phase, and the second multiplication value is less than the mean magnitude according to the comparison result (first multiplication value&lt;mean phase &amp; second multiplication value&lt;mean magnitude), the masking unit  130  may perform masking of the corresponding pixel of the color image. 
         [0040]    Conversely, when the first multiplication value is greater than the mean phase and the second multiplication value is less than the mean magnitude according to the comparison result (first multiplication value&gt;mean phase &amp; second multiplication value&lt;mean magnitude), the masking unit  130  may not perform masking of the corresponding pixel of the color image. 
         [0041]    This may be simply shown by Table 1 as given below: 
         [0000]    
       
         
               
               
               
               
               
             
               
               
               
               
               
             
           
               
                   
                 TABLE 1 
               
               
                   
                   
               
               
                   
                 Comparison between first 
                   
                 Comparison between second 
                   
               
               
                   
                 multiplication value with 
                   
                 multiplication value with 
               
               
                   
                 mean phase 
                   
                 mean magnitude 
                 Making value 
               
               
                   
                   
               
             
             
               
                   
               
             
          
           
               
                 1 
                 mean phase &lt; scale 
                 &amp; 
                 mean magnitude &gt; scale 
                 0 (Masking) 
               
               
                   
                 factor × max mean phase 
                   
                 factor × max mean magnitude 
               
               
                 2 
                 mean phase &gt; scale 
                 &amp; 
                 mean magnitude &lt; scale 
                 1 (No Masking) 
               
               
                   
                 factor × max mean phase 
                   
                 factor × max mean magnitude 
               
               
                 3 
                 mean phase &gt; scale 
                 &amp; 
                 mean magnitude &gt; scale 
                 0 (Masking) 
               
               
                   
                 factor × max mean phase 
                   
                 factor × max mean magnitude 
               
               
                 4 
                 mean phase &lt; scale 
                 &amp; 
                 mean magnitude &lt; scale 
                 0 (Masking) 
               
               
                   
                 factor × max mean phase 
                   
                 factor × max mean magnitude 
               
               
                   
               
             
          
         
       
     
         [0042]    As shown in Table 1, the ultrasound color Doppler image system  100  may set a masking value as “0” or “1” depending on the comparison result between the first multiplication value and the mean phase and the comparison result between the second multiplication value and the mean magnitude. Accordingly, the masking unit  130  may perform masking of the corresponding pixel of the color image, when the masking value is “0”. Conversely, the masking unit  130  may not perform masking of the corresponding pixel of the color image, when the masking value is “1”. 
         [0043]    The controller  140  may generally control an operation of the ultrasound color Doppler image system  100 , such as operations of the calculation unit  110 , the comparison unit  120 , the masking unit  130 , and the like. 
         [0044]      FIGS. 2 through 4  are diagrams illustrating an example of a phase and a magnitude of two successive signals in an ensemble I/Q signal according to an exemplary embodiment of the present invention. 
         [0045]    A graph of modulated I/Q signals, namely, a graph of E( 0 ), E( 1 ), E( 2 ), E( 3 ), . . . , E(EL−2), and E(EL−1), is illustrated in  FIGS. 2 through 4 , a number of the modulated I/Q signals being same as an ensemble length that is a number of ultrasound beams to be used for calculating a color image of a single pixel. A phase and a magnitude of two temporally successive I/Q signals may be calculated using the two signals. 
         [0046]    First, as illustrated in  FIG. 2 , the ultrasound Doppler image system  100  according to an exemplary embodiment of the present invention may calculate a phase and a magnitude of E( 0 ) and E( 1 ) using the two temporally successive signals, namely, E( 0 ) and E( 1 ). 
         [0047]    Specifically, the ultrasound color Doppler image system  100  may multiply the two temporally successive signals, namely, E( 0 ) and E( 1 ), (E( 0 )×E( 1 )*), and may calculate a phase 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       ( 
                       
                         
                           C 
                            
                           
                             ( 
                             0 
                             ) 
                           
                         
                         = 
                         
                            
                           
                             
                               
                                 ∠tan 
                                 
                                   - 
                                   1 
                                 
                               
                                
                               
                                 ( 
                                 
                                   b 
                                   
                                     ? 
                                   
                                 
                                 ) 
                               
                             
                             - 
                             
                               
                                 ∠tan 
                                 
                                   - 
                                   1 
                                 
                               
                                
                               
                                 ( 
                                 
                                   d 
                                   c 
                                 
                                 ) 
                               
                             
                           
                            
                         
                       
                       ) 
                     
                      
                     
                       
 
                     
                      
                     
                       
                         ? 
                       
                        
                       
                         indicates text missing or illegible when filed 
                       
                     
                   
                    
                   
                       
                   
                 
               
               
                 
                     
                 
               
             
           
         
       
     
         [0000]    and a magnitude (M( 0 )=√{square root over (a 2 +b 2 )}×√{square root over (c 2 +d 2 )}) of the two signals (E( 0 ) and E( 1 )) from a multiplication result 
         [0000]    
       
         
           
             
               ( 
               
                 
                   
                     
                       a 
                       2 
                     
                     + 
                     
                       b 
                       2 
                     
                   
                 
                 × 
                 
                   
                     
                       c 
                       2 
                     
                     + 
                     
                       d 
                       2 
                     
                   
                 
                  
                 
                    
                   
                     ∠ 
                      
                     
                        
                       
                         
                           
                             tan 
                             
                               - 
                               1 
                             
                           
                            
                           
                             ( 
                             
                               b 
                               g 
                             
                             ) 
                           
                         
                         - 
                         
                           
                             tan 
                             
                               - 
                               1 
                             
                           
                            
                           
                             ( 
                             
                               d 
                               c 
                             
                             ) 
                           
                         
                       
                        
                     
                   
                 
               
               ) 
             
             . 
           
         
       
     
         [0048]    As illustrated in  FIG. 3 , the ultrasound color Doppler image system  100  may calculate a phase and a magnitude of temporally successive two signals, namely, E( 1 ) and E( 2 ) using E( 1 ) and E( 2 ). 
         [0049]    Specifically, the ultrasound color Doppler image system  100  may multiply the two temporally successive signals, namely, E( 1 ) and E( 2 ), (E( 1 )×E( 2 *), and may calculate a phase 
         [0000]    
       
         
           
             ( 
             
               
                 C 
                  
                 
                   ( 
                   1 
                   ) 
                 
               
               = 
               
                  
                 
                   
                     
                       ∠tan 
                       
                         - 
                         1 
                       
                     
                      
                     
                       ( 
                       
                         
                           ? 
                         
                          
                         
                           ? 
                         
                       
                       ) 
                     
                   
                   - 
                   
                     
                       ∠tan 
                       
                         - 
                         1 
                       
                     
                      
                     
                       ( 
                       
                         d 
                         c 
                       
                       ) 
                     
                   
                 
                  
               
             
             ) 
           
         
       
       
         
           
             
               
                 ? 
               
                
               
                 indicates text missing or illegible when filed 
               
             
              
             
                 
             
           
         
       
     
         [0000]    and a magnitude (M( 1 )=√{square root over (e 2 +f 2 )}×√{square root over (c 2 +d 2 )}) of the two signals (E( 1 ) and E( 2 )) from a multiplication result 
         [0000]    
       
         
           
             
               
                 
                   ( 
                   
                     
                       
                         
                            
                           2 
                         
                         + 
                         
                           f 
                           2 
                         
                       
                     
                     × 
                     
                       
                         
                           c 
                           2 
                         
                         + 
                         
                           d 
                           2 
                         
                       
                     
                      
                     
                       ? 
                     
                   
                   ) 
                 
                 . 
                 
                   
 
                 
                  
                 
                   ? 
                 
               
                
               
                 indicates text missing or illegible when filed 
               
             
              
             
                 
             
           
         
       
     
         [0050]    The ultrasound color Doppler image system  100  may continuously calculate a phase and a magnitude of each of two other signals, such as E( 2 ) and E( 3 ), E( 3 ) and E( 4 ), . . . , and E(EL−3) and E(EL−2) in the same manner as described with reference to  FIGS. 2 and 3 , and when a phase and a magnitude of E(EL−2) and E(EL−1) is calculated based on a method of  FIG. 4 , phases and magnitudes of each pixel or each sample in an ensemble I/Q signal, namely the color image, are completely calculated. 
         [0051]    As illustrated in  FIG. 4 , the ultrasound color Doppler image system  100  may multiply two temporally successive signals, namely, E(EL−2) and E(EL−1), ((E(EL−2)×E(EL−1)*), and may calculate a phase 
         [0000]    
       
         
           
             ( 
             
               
                 C 
                  
                 
                   ( 
                   
                     EL 
                     - 
                     2 
                   
                   ) 
                 
               
               = 
               
                  
                 
                   
                     
                       ∠tan 
                       
                         - 
                         1 
                       
                     
                      
                     
                       ( 
                       
                         n 
                         m 
                       
                       ) 
                     
                   
                   - 
                   
                     
                       ∠tan 
                       
                         - 
                         1 
                       
                     
                      
                     
                       ( 
                       
                         ? 
                       
                       ) 
                     
                   
                 
                  
               
             
             ) 
           
         
       
       
         
           
             
               
                 ? 
               
                
               
                 indicates text missing or illegible when filed 
               
             
              
             
                 
             
           
         
       
     
         [0000]    and a magnitude (M(EL−2)=√{square root over (m 2 +n 2 )}×√{square root over (k 2 +h 2 )}) of the two signals (E(EL−2), E(EL−1)) from a multiplication result 
         [0000]    
       
         
           
             
               
                 
                   ( 
                   
                     
                       
                         
                           m 
                           2 
                         
                         + 
                         
                           n 
                           2 
                         
                       
                     
                     × 
                     
                       
                         
                           k 
                           2 
                         
                         + 
                         
                           
                             h 
                             2 
                           
                            
                           
                             ? 
                           
                         
                       
                     
                      
                     
                       ? 
                     
                   
                   ) 
                 
                 . 
                 
                   
 
                 
                  
                 
                   ? 
                 
               
                
               
                 indicates text missing or illegible when filed 
               
             
              
             
                 
             
           
         
       
     
         [0052]    As described above, the ultrasound color Doppler image system  100  may calculate phases and a magnitudes of each pixel or each sample in the ensemble I/Q signal, namely, the color image, using two temporally successive signals. Accordingly, the ultrasound color Doppler image system  100  may calculate the mean phase and the mean magnitude as illustrated in  FIG. 3  based on the phases and magnitudes calculated in  FIGS. 2 through 4 . 
         [0000]    
       
         
           
             
               
                 
                   
                     MeanPhase 
                     = 
                     
                       
                         ( 
                         
                           
                             ∑ 
                             
                               k 
                               = 
                               0 
                             
                             
                               EL 
                               - 
                               1 
                             
                           
                            
                           
                             C 
                              
                             
                               ( 
                               k 
                               ) 
                             
                           
                         
                         ) 
                       
                       n 
                     
                   
                    
                   
                     
 
                   
                    
                   
                     MeanMagnitude 
                     = 
                     
                       
                         
                           ( 
                           
                             
                               ∑ 
                               
                                 k 
                                 = 
                                 0 
                               
                               
                                 EL 
                                 - 
                                 1 
                               
                             
                              
                             
                               C 
                                
                               
                                 ( 
                                 k 
                                 ) 
                               
                             
                           
                           ) 
                         
                         n 
                       
                       . 
                     
                   
                 
               
               
                 
                   [ 
                   
                     Equation 
                      
                     
                         
                     
                      
                     3 
                   
                   ] 
                 
               
             
           
         
       
     
         [0053]    Here, Meanphase denotes a mean phase, MeanMagnitude denotes a mean magnitude, n denotes a total number of phases and a total number of magnitudes, 
         [0000]    
       
         
           
             
               ∑ 
               
                 k 
                 = 
                 0 
               
               
                 EL 
                 - 
                 1 
               
             
              
             
               C 
                
               
                 ( 
                 k 
                 ) 
               
             
           
         
       
     
         [0000]    denotes a total sum of phases, and 
         [0000]    
       
         
           
             
               ∑ 
               
                 k 
                 = 
                 0 
               
               
                 EL 
                 - 
                 1 
               
             
              
             
               M 
                
               
                 ( 
                 k 
                 ) 
               
             
           
         
       
     
         [0000]    denotes a total sum of magnitudes. 
         [0054]      FIG. 5  illustrates an example of masking a pixel in a color image according to an exemplary embodiment of the present invention. For reference, in  FIG. 5 , each square block denotes a pixel, and a number in the square block denotes a pixel number. Also, numbers “0” and “1” in a right side of a black block and a white block denote a masking value. 
         [0055]    As illustrated in  FIG. 5 , a masking value corresponding to pixel numbers  1 ,  2 ,  4 ,  5 ,  9 ,  11 ,  12 ,  13 , and  14  is “0”, and a masking value corresponding to pixel numbers  3 ,  6 ,  7 ,  8 ,  10 ,  15 , and  16  is “1”. 
         [0056]    Accordingly, the ultrasound color Doppler image system  100  may perform masking of pixels masking values of which are “0”, namely, the pixels corresponding to the pixel numbers  1 ,  2 ,  4 ,  5 ,  9 ,  11 ,  12 ,  13 , and  14 . The corresponding pixels are marked in oblique lines and the corresponding pixels may be marked in black when actually embodied. 
         [0057]    Conversely, the ultrasound color Doppler image system may not perform masking of pixels masking values of which are “1”, namely, the pixels corresponding to the pixel numbers  3 ,  6 ,  7 ,  8 ,  10 ,  15 , and  16 . 
         [0058]    As described above, the ultrasound color Doppler image system  100  may set the masking value through a process, such as calculation, comparison, and the like, with respect to each pixel of the color image, and may perform masking of a corresponding pixel according to a predetermined masking value. Accordingly, a flash artifact component may be effectively eliminated from a color Doppler image according to an exemplary embodiment of the present invention. 
         [0059]    As described above, the ultrasound color Doppler image system may effectively eliminate the flash artifact component from the color Doppler image, thereby providing a more accurate Doppler image as illustrated in  FIG. 6 . For reference,  FIG. 6  illustrates an image (left image) before eliminating a flash artifact from a color Doppler image and an image (right image) after eliminating the flash artifact from the color Doppler image according to an exemplary embodiment of the present invention 
         [0060]      FIG. 7  illustrates a method of filtering a clutter signal in an ultrasound color Doppler image system according to an exemplary embodiment of the present invention. 
         [0061]    Referring to  FIG. 7 , the ultrasound color Doppler image system  100  may calculate a mean phase and a mean magnitude of a first signal corresponding to a pixel of a color image in the ultrasound image in operation  5710 . Here, the color image may be an image before a clutter signal is eliminated, namely, an image including a blood flow and a peripheral tissue, such as a vessel wall, a muscle, and the like. Also, the first signal may be an I/Q signal per pixel or an I/Q signal per sample among an I/Q signal based on an ensemble unit. 
         [0062]    In this instance, the ultrasound color Doppler image system may calculate a mean phase and a mean magnitude of the first signal using the first signal and a second signal that has a temporal continuity with the first signal. 
         [0063]    Subsequently, the ultrasound color Doppler image system  100  may generate a first multiplication value and a second multiplication value by multiplying each of a maximum value for the calculated mean phase and a maximum value for the calculated mean magnitude by a predetermined scale factor. Here, the scale factor may be set by a user in advance, and may include positive numbers, such as 1, 1.5, 2, 2.5, and the like, as coefficients for scaling. 
         [0064]    Specifically, the ultrasound color Doppler image system  100  may generate the first multiplication value by multiplying the scale factor and the maximum value for the calculated mean phase, namely a maximum mean phase (scale factor×maximum mean phase). Also, the ultrasound color Doppler image system may generate the second multiplication value by multiplying the scale factor and the maximum value for the calculated mean magnitude, namely a maximum mean magnitude (scale factor×maximum mean magnitude). 
         [0065]    Subsequently, the ultrasound color Doppler image system  100  may compare the generated first multiplication value and the generated second multiplication value with each of the mean phase and the mean magnitude. 
         [0066]    Specifically, the ultrasound color Doppler image system  100  may compare the first multiplication value with the mean phase, and may compare the second multiplication value with the mean magnitude. 
         [0067]    Subsequently, the ultrasound color Doppler image system may perform masking of the corresponding pixel of the color image based on the comparison result in operation S 740 . 
         [0068]    As an example, when the first multiplication value is less than the mean phase and the second multiplication value is greater than the mean magnitude according to the comparison result (first multiplication value&lt;mean phase &amp; second multiplication value&gt;mean magnitude), the ultrasound color Doppler image system  100  may perform masking of the corresponding pixel of the color image. 
         [0069]    When the first multiplication value is greater than the mean phase, and the second multiplication value is greater than the mean magnitude according to the comparison result (first multiplication value&gt;mean phase &amp; second multiplication value&gt;mean magnitude), the ultrasound color Doppler image system  100  may perform masking of the corresponding pixel of the color image. 
         [0070]    When the first multiplication value is less than the mean phase, and the second multiplication value is less than the mean magnitude according to the comparison result (first multiplication value&lt;mean phase &amp; second multiplication value&lt;mean magnitude), the ultrasound color Doppler image system may perform masking of the corresponding pixel of the color image. 
         [0071]    Conversely, when the first multiplication value is greater than the mean phase, and the second multiplication value is less than the mean magnitude according to the comparison result (first multiplication value&gt;mean phase &amp; second multiplication&gt;mean magnitude), the ultrasound color Doppler image system  100  may not perform masking of the corresponding pixel of the color image. 
         [0072]    As described above, the ultrasound color Doppler image system  100  may determine whether to perform the masking of each pixel in the color image through a calculation and a comparison, and may perform masking of a corresponding pixel. Accordingly, example embodiments may effectively eliminate a flash artifact component from the color Doppler image by filtering a clutter signal from a Doppler signal. 
         [0073]    The exemplary embodiments of the present invention include computer-readable media including program instructions to implement various operations embodied by a computer. The media may also include, alone or in combination with the program instructions, data files, data structures, tables, and the like. The media and program instructions may be those specially designed and constructed for the purposes of the present invention, or they may be of the kind well known and available to those having skill in the computer software arts. Examples of computer-readable media include magnetic media such as hard disks, floppy disks, and magnetic tape; optical media such as CD ROM disks; magneto-optical media such as optical disks; and hardware devices that are specially configured to store and perform program instructions, such as read-only memory devices (ROM) and random access memory (RAM). Examples of program instructions include both machine code, such as produced by a compiler, and files containing higher level code that may be executed by the computer using an interpreter. The described hardware devices may be configured to act as one or more software modules in order to perform the operations of the above-described embodiments of the present invention, or vice versa. 
         [0074]    Although a few example embodiments have been shown and described, it would be appreciated by those skilled in the art that changes may be made in these example embodiments without departing from the principles and spirit of the invention, the scope of which is defined in the claims and their equivalents.