Method for removing artifacts from Doppler ultrasound images.

Abstract

A computer-implemented noise reduction method for removing artifacts from Doppler-type ultrasound images includes the steps of: acquiring a plurality (n) of frames comprising ultrasound images; identifying pixels containing Doppler activation; generating a first image data set for each frame of the plurality of frames by classifying the pixels of the corresponding frames into categories and associating at least a first or second category with each pixel position in the corresponding frame to identify, for each frame, pixels containing Doppler activation; generating at least duration sequences associated with frame portions containing identification data ordered according to the frame sequence, such that adjacent identification data in the sequence refer to consecutive frames, wherein each identification data assumes a first value if the corresponding frame portion contains a Doppler activation pixel, or a second value different from the first value if the corresponding frame portion does not contain a Doppler activation pixel; calculating the length of each duration subsequence having consecutive identification data with a first value in the duration sequence; automatically calculating a reference threshold; and rejecting Doppler pixels if the Doppler pixels belong to a duration sequence below the threshold.

Description

[Technical Field]

[0001] The present invention relates to a method for removing artifacts contained in Doppler images from ultrasound examinations. [Background technology]

[0002] Ultrasound imaging technology is commonly used to generate diagnostic images of the internal features of objects, such as the structure of the human body.

[0003] Ultrasound imaging is a technique that reproduces anatomical cross-sections of the human body corresponding to the scanning plane on a two-dimensional image. The mechanism by which information contained in the echo signal is converted into a two-dimensional image is complex and consists of several steps. Some of these depend on the propagation of ultrasound in biological tissue, while others depend on the equipment technology or how the operator handles them. Ultrasound is generated and properly focused by a transducer. Ultrasound propagates within tissue at a generally constant speed, with slight variations depending on density or acoustic impedance. Contact between ultrasound and various anatomical structures initiates various physical phenomena, such as reflection, dispersion, diffraction, and absorption. The first two of these phenomena generate echo signals that return to the transducer, which are then properly processed to create an ultrasound image. Basically, the information contained in the original signal is not sufficient to create an ultrasound image. However, it must be integrated with other essential information: firstly, the depth of the signal source, which is established as a function of the time elapsed between beam emission and the return of the echo signal; and secondly, the signal direction along the lines that constitute the scanning plane. In other words, in order to generate an ultrasound image, these three pieces of information need to be obtained for each point in the scanning area.

[0004] Ultrasound imaging devices are high-performance instruments capable of handling several imaging modes. A simplified block diagram of an ultrasound imaging device is shown below: - A pulse generator, responsible for producing electrical signals; - A piezoelectric transducer driven to perform both the functions of a source and a receiver. As a source, the transducer is responsible for transmitting an ultrasonic acoustic signal generated by converting electrical pulses into waves. Acoustic wave pulses traveling through tissue are partially reflected by materials with different acoustic impedances. Thereafter, as a receiver, the transducer is responsible for detecting the reflected pulses and converting them into radio frequency (RF) electrical signals; - A TGC module that acts on the received signal by amplifying it proportionally to the depth of the echo signal generation, and by amplifying the received signal to a signal return time that compensates for the gain, in order to accurately represent the difference in acoustic impedance of the surfaces constituting the discontinuity. In other words, a signal reaching the probe from a surface in a deep sound field is weaker than an echo signal emitted from a reflector located on a shallow surface with the same characteristics; - A demodulator that converts the amplified RF signal into a single peak representing the signal from each reflector the ultrasound hit. This allows for the identification of all points where reflections occur and the quantification of the echo amplitude; - A scanning converter that converts a time-continuous input signal, samples it at specified intervals, and converts it into a separate numerical format encoded by signal intensity with 16, 64, and 256 gray levels according to a format, which can be displayed on a monitor so that the image is acquired along lines that make up the scanning field, and also displayed along horizontal lines on the monitor. This is realized by two main macrophases: in the first input phase, the scanning converter converts the RF signal from analog to digital and stores the data in binary format in a memory matrix; in the second output phase, it performs a digital-to-analog conversion on the data stored in the memory matrix; - A display / storage system that typically allows tracking to be viewed via a monitor and then stored in a suitable storage medium for later viewing.

[0005] A specific application of diagnostic ultrasound imaging is the detection and display of fluid velocity using Doppler measurements. This is a physical phenomenon in which the frequency of sound waves striking a moving body fluctuates, and it is directly proportional to the speed of the body's movement itself.

[0006] Therefore, the Doppler effect is based on measuring the frequency variation between the input beam and the reflected beam from a moving body (in medical applications, the moving body often represents red blood cells, and thus blood).

[0007] Therefore, the Doppler effect is the principle behind many ultrasound techniques used to investigate the movement of biological tissues. The variation in the reflected echo frequency is related to the velocity and direction of the reflecting body. The received echo signal is compared to a reference value to establish the fluid velocity through the region. The fluid velocity can be established by measuring the shift in the Doppler frequency of the echo signal received from the structure under inspection.

[0008] There are two basic Doppler acquisition systems: - Continuous waves (CW) enable the investigation of fluid patterns and motion; - A pulsed wave (PW) also provides information about the distance between the reflecting surface and the transducer.

[0009] Color Doppler operates on the same principle, integrating flow information with real-time two-dimensional images. Through conversion, approaching flows are displayed in red, and receding flows in blue. As known from prior art, flow velocity can be visualized through coloring, where different shades and color intensities represent flow rate and direction within the grayscale image. For example, if disturbance is present at a vascular bifurcation, an alternating pattern of blue and red sections will be present. Thus, it becomes possible to distinguish the flow direction relative to the probe. These systems further assist in evaluating the volume of flow itself and the duration of laminar flow or disturbance.

[0010] The importance of obtaining the most accurate results possible from ultrasound machines is evident from the fact that they are sometimes used by doctors for accurate diagnoses.

[0011] This highlights two main categories of challenges in these systems: on the one hand, the experience of a physician or operator using an ultrasound machine that must be operated with parameters set to obtain effective results, and on the other hand, the system itself.

[0012] To achieve the most even distribution of image brightness by amplifying signals from deeper layers, maximizing contrast, and avoiding saturation, the parameters of an ultrasound machine that a physician or technician can set include, for example, gain control and ultrasound focus settings. In an ultrasound machine, there are commands that allow the physician or technician to adjust the scale settings used by the TGC according to their individual working habits, thereby improving the depiction of deeper layers.

[0013] This parameter, understood as the signal-to-noise ratio, directly affects the signal background. The signal-to-noise ratio is generated by the electrical circuits in the generation phase, the reception of the ultrasound signal, and processing, and it appears in the image as a false representation of the flow signal in the vascular region by offsetting the characteristics of the echo structure, which increases with excessive signal amplification or gain values. An example of this false representation is the extreme artifact shown in Figure 5e. Focusing settings allow for changes to the number of active focus points, the shape of the ultrasound beam, and the thickness. If not set correctly, this affects the spatial resolution of the ultrasound image. Another setting that can be configured is the imaging frequency. Increasing its value allows for optimization of the resolution, while decreasing its value allows for increased beam penetration. Another parameter that the physician or operator can act on concerns the compression curve, which allows for changes to the correspondence between signal amplitude and the gray level displayed on the monitor. Using a linear intensity / amplitude relationship, the grayscale is directly proportional to the difference in acoustic impedance of the examined tissue. In some cases, this is not the optimal setting for the diagnostic need, and therefore the physician or operator may represent several acoustic impedance ranges, which are larger traces than others of less importance, by, for example, compressing higher intensity echoes and allowing lower intensity echoes to be represented with a higher grayscale number. A sensitive parameter in clinical examination is the "pulse repetition rate" (or "flow sensitivity"), which controls the system's performance in acquiring flow at varying speeds. If the pulse frequency is too high, sensitivity to slow flow decreases, resulting in a corresponding loss of signal.

[0014] Other important categories belonging to this system, or limitations derived from this technique, primarily due to the interaction between ultrasound and biological structures, will be highlighted in the following description of possible artifacts that may appear in the results of ultrasound examinations. Significant challenges at the instrument's origin, and / or challenges relating to the physician's or technician's experience, are reflected in the results produced by the ultrasound machine itself, both in the use phase and the configuration phase, and are usually identified by the name of artifact. An artifact refers to erroneous or distorted information generated by the ultrasound machine, or by the interaction of ultrasound with tissue overlaid with noise in the Doppler signal. In detail, the physician's or technician's experience in using the ultrasound machine directly leads to greater or lesser variations in the amount of artifact, depending, for example, on probe positioning, movement made by the probe itself, and the speed at which such movement is performed.

[0015] Artifacts seen in Doppler ultrasound imaging, for example in color Doppler ultrasound images, can, from the perspective of flow information, cause confusion or misjudgment considering that this artifact can be defined as all colored pixels and it does not accurately represent the vascular distribution. As mentioned above, three main factors contribute to this problem. That is, improper configuration of equipment and inappropriate signal acquisition due to human error, anatomical factors, and technical constraints. For example, incorrect gain settings, wall filter settings, or velocity scale settings can result in losses of clinically relevant information such as tracking of distortions that display situations that can vary greatly from the presence or absence of flow in blood vessels, flow direction and velocity, or the actual physiological situation. When referring to artifacts due to inappropriate acquisition, motion artifacts are particularly relevant. These include, for example, those due to improper irradiation angle setting and those due to overly fast signal acquisition which is the most frequent error causing flash artifacts. Regarding the issue caused by the irradiation angle which is the angle between the operator's hand and the probe, when the angle is greater than 60°, the amplitude of the spectral curve is gradually reduced and the velocity calculation becomes gradually less reliable. On the other hand, when the angle approaches 90°, even if the flow can generate a signal with a low amplitude, the signal is not recorded. The following main artifacts were found in color Doppler ultrasound according to the categories defined above. Artifacts due to poor configuration

[0016] Doppler gain setting error: Appropriate gain setting is important for accurate reproduction of flow characteristics. If the gain setting is too low, some relevant information may be lost, so the gain is frequently adjusted to maximize the visualization of tracking. In contrast, too high a gain degrades the envelope signal, interrupts its reproduction on the screen, and mimics the spectral broadening that can manifest, for example, turbulent flow.

[0017] Improper angle setting: Regarding the issue of the irradiation angle, it is possible to set correction parameters in the ultrasonic machine. However, if the angle is completely incorrect, this correction will also be in vain. The occurrence of this type of artifact may also depend on using a converter with a very high frequency or no gain adjustment.

[0018] Improper filter setting: The filter phase is designed to remove low-frequency Doppler signals coming from echoes of slow soft tissues. The cut-off frequency of this filter can be selected by the operator. If the setting is too high, diagnostically relevant velocity information may be lost.

[0019] Artifact due to spectral dispersion: Spectral dispersion may occur due to excessive system gain or changes in grayscale sensitivity. Artifacts based on biological structures

[0020] Flash artifact: This artifact occurs as a sudden color burst that extends to a somewhat wide scanning field area. The color coding is a completely artificial result and can occur due to too fast movement of the converter, or due to heart movement, or venous pulsations that cause slight movement of the reflecting surface.

[0021] Pseudo-flow artifact: Pseudo-flow represents actual fluid motion that is not blood. These types of artifacts may be observed as a result of excessive amplification of the Doppler signal, overly low color coding (which may result in signals appearing within cystic structures, for example), adjacent moving structures, respiratory dynamics. It can also occur due to the presence of mirror artifacts or fluctuations. Artifacts that rely on technical limitations

[0022] Directional ambiguity: Directional ambiguity can occur when the ultrasound beam obstructs a blood vessel at a 90° angle. The detected Doppler signal appears as an up-and-down tracking at the baseline of the spectrum. Furthermore, at high gain settings, directional ambiguity is more pronounced, and tracking becomes even less precise. When color Doppler images are generated using a sector-type transducer, flow perpendicular to the beam is generally present along small segments of blood vessels parallel to the transducer surface, while the challenge becomes more pronounced when using a linear probe.

[0023] Side lobe artifact: An electrically focused array transducer directs the primary beam to the critical area to be inspected. However, due to the array elements, weak secondary ultrasonic lobes in space can target areas unrelated to the primary beam. The precise location of these lobes depends on the transducer's design. If these secondary lobes hit a highly reflective surface (such as bone), the echoes returning to the transducer may be detected on the screen along with the primary beam echoes.

[0024] Random noise artifacts: In Doppler ultrasound machines, as with all electrical circuits, noise is proportional to the gain. Random noise occurs through the appearance of flash artifacts or the appearance / disappearance of color regions, especially when the gain is set too high, whenever there is interaction between the probe and tissue.

[0025] Sparkling artifacts: This type of artifact is visible in highly reflective structures and manifests itself as a fluctuating color mosaic, related to the signal characteristics in the presence of background noise. Its appearance is strictly dependent on the ultrasonic mechanical configuration and is generated by a narrowband of inherent electronic noise called phase jitter.

[0026] In the prior art, several signal processing techniques are known that are used to filter out undesirable signals within an ultrasound imaging system that cause one or more artifacts as described above. However, these solutions, performed as part of the image display generation process within an ultrasound machine, either as an upstream processing phase and / or filtering phase, exhibit self-destructive effects in more or less related ways, such as loss of instrument sensitivity and / or relevant information for diagnostic purposes.

[0027] Furthermore, the removal of potentially different types of artifacts relies on various combinations of ultrasonic machine parameter configurations, available processing, and enhanced filtering based on the ultrasonic machine model. [Overview of the project] [Problems that the invention aims to solve]

[0028] The object of the present invention is to solve at least some of the shortcomings of the aforementioned instrument or human factors by acting on the ultrasonic image generated by the ultrasonic machine, instead of acting on the signal information and / or image generation by onboard filtering and / or processing blocks. This result can be achieved through a denoising step in the image generated by the ultrasonic machine acting on 2D reproduction. This analysis aims to identify any artifacts present in the image in order to remove artifacts and / or minimize the impact of suboptimal settings and / or operator experience and / or technical limitations, thereby providing an artifact-free image or an image with significantly reduced artifacts. [Means for solving the problem]

[0029] Therefore, noise reduction is the analysis of Doppler ultrasound images, and through the steps of the method described below, the goal is to identify the presence of any possible artifacts in the image that can be superimposed on the actual signal or positioned alongside the actual signal to remove artifacts, thereby obtaining an image free from distortions that could interfere with or mislead the diagnostic phase. Other technical terms included in this invention that are worth defining and that help relate their use within the description are: - Doppler activation, i.e., a set of pixels that retain Doppler artifacts such as Doppler signals and colored pixels; - A set of pixels that hold non-Doppler ultrasonic signals, such as grayscale pixels, belonging to the outer region from which ultrasonic signals, i.e., Doppler signals, are accumulated; - A set of pixels that hold the Doppler signal, i.e., the actual Doppler signal indicating the presence of blood vessel distribution.

[0030] The computer-based noise reduction method of the present invention begins by decomposing the video into individual frames using standard video segment decomposition techniques, and includes the following steps: - Steps to obtain n frames from ultrasound images; - For each of the n frames, the steps are to identify a first pixel category representing a first pixel in the nth frame containing an ultrasonic signal, and a second pixel category representing a second pixel in the nth frame containing Doppler activation; - For each of the multiple frames mentioned above, the steps include: generating a first set of image data to classify the pixels of the corresponding frames and to associate at least a first or second category with the location of each pixel in the corresponding frame, in order to identify pixels containing Doppler activation for each frame; - A step of generating at least a persistent sequence associated with a frame portion containing identification data ordered according to a frame sequence such that adjacent identification data points to adjacent frames, wherein each identification data is assumed to have a first value if the corresponding frame portion contains a Doppler-activated pixel, or a second value different from the first value if the corresponding frame portion does not contain a Doppler-activated pixel; - A step of calculating the length of each persistence subsequence, which in the above persistence sequence has consecutive identification data having a first value; - A step of calculating a reference threshold based on the length distribution of the persistence subsequences described above, which refer to the identification data of all n frames; -For each of the n frames: If frame n corresponds to one of the above persistence subsequences having a length longer than the above threshold so as to represent a pixel with an actual Doppler signal, then the above pixel position contains information belonging to the above second category; or If frame n corresponds to another of the above persistence subsequences having a length shorter than the threshold so as to represent an artifact, then the information that the above pixel position belongs to the above first category is, A step of identifying a second set of image data (SEED) to associate with each pixel position of the corresponding frame.

[0031] This method is based on the persistence and identification data of the Doppler signal across n frames. While this identification data is considered an artifact even if it belongs to the first value indicating the presence of Doppler activation, its persistence is relatively weak across n frames. Because the identification data belongs to the first value indicating the presence of Doppler activation, the data is more persistent across n frames.

[0032] The step of generating a temporally sustained sequence of Doppler activation is applied to a portion of each frame. According to one embodiment, such a region is a single pixel in the frame, and the more pixels present in the frame of the video format, the more sustained sequences are generated. That is, the pixel matrix that makes up the frame is, for example, an 800x566 format. Each sustained sequence contains identification data for each frame of the video. According to another embodiment, the portion of the frame that the sustained sequence points to is a connected component of pixels with Doppler activation.

[0033] In this second event, it is necessary to define the above frame portion by performing the step of running a search algorithm for connecting elements, in particular for pixels of the first category in frame n and frame n+1. The step of generating a persistence sequence here includes the steps of comparing a first parameter of the first connecting element in frame n with a second connecting element in frame n+1, and generating a second persistence sequence for the second connecting element based on the above comparison step in order to associate the second connecting element with a first persistence sequence in the first connecting element, or to realize tracking of the connecting element between frames n and n+1. The above parameters are preferably the centroid and / or a parameter representing the overlap of the first and second connecting elements, and / or a parameter representing the similarity of the shapes of the first and second connecting elements, and / or a parameter representing the size or dimensions of the first and second connecting elements.

[0034] Thus, the connecting components enable tracking or monitoring of the Doppler activation region, and processing accuracy is improved in the presence of Doppler signal regions that change during Doppler data acquisition.

[0035] In a preferred embodiment where the frame portion is pixels, the method includes the step of extending a second pixel dataset with other pixels containing the Doppler signal from which pixels representing artifacts have been pre-removed, based on a second image dataset (SEED) for the nth frame, by performing a segment decomposition algorithm in at least one nth frame. In particular, the contiguous element analysis is also performed in the second dataset for each frame, referencing pixels representing the actual Doppler signal. For each frame, the contiguous elements of the video frame and the contiguous elements of the pixels in the actual Doppler signal pixels of the actual Doppler signal are superimposed to define, for example, a 90% superimposition threshold. Contiguous elements of video frames that pass the superimposition test and do not function are retained and, for example, modified. Pixels in the Doppler-activated contiguous elements of video frames that did not pass the superimposition test are therefore considered artifacts.

[0036] In detail, it has been found that the second set of image data tends to underestimate the range of representative pixels with actual Doppler signals (i.e., Doppler signals with pixels containing artifacts excluded). To improve the accuracy of processing the original frame, a segment decomposition algorithm is applied to extend the second set of image data to merge additional pixels with previously removed Doppler signals, as shown in the previous paragraph. These pixels are considered to represent artifacts.

[0037] According to a preferred embodiment, each acquired frame can be represented by a three-dimensional matrix, where the first dimension represents the number of pixels on the vertical axis, the second dimension represents the number of pixels on the horizontal axis, and finally the third dimension represents the number of channels (R, G, B), thereby indicating each pixel with three values.

[0038] According to embodiments of the present invention, the identification of colored pixels in step 2 of the noise reduction method, which represent pixels containing Doppler signals, can be achieved by dividing the R, G, and B values ​​of the pixels constituting the frame into clusters, and by selecting pixels that are not classified into clusters that identify the light-colored scale of the ultrasonic machine.

[0039] According to another embodiment of the present invention, the identification of colored pixels in step 2 of the noise reduction method can also be achieved by selecting pixels in which the intensity difference between the R and G, G and B, or R and B channels is greater than or equal to a predetermined threshold.

[0040] According to a preferred embodiment, the threshold is calculated as the 90th percentile of the distribution. This is because the objective is to obtain the minimum activation length value of the Doppler signal, thereby assuming that activation sequence lengths longer than this threshold correspond to the actual signal, while activation sequence lengths shorter than the threshold are assumed to be artifacts.

[0041] According to another preferred embodiment, when the minimum activation length of the Doppler signal is derived by tracking the temporal duration of connected components contained within a closed and isolated Doppler region, i.e., a frame, the threshold is calculated as the 98th percentile.

[0042] According to another embodiment, the embodiment threshold in both frame portions is calculated as the sum of the mean values ​​of representative Doppler activation subsequence lengths over a logarithmic scale, with a 2x standard deviation, in order to compensate for the asymmetry of the distribution.

[0043] In the step of determining the actual signal location, according to a preferred embodiment, the obvious result in the denoised frame may include Doppler signal pixels with the color channel [R, G, B] values ​​of the original frame. Meanwhile, other pixels would assume the average values ​​in the three channels, which are classified as grayscale.

[0044] According to a preferred embodiment in the growth step region, a morphological collapse method is used, followed by a morphological expansion method, with a kernel diameter or side of 5 pixels and a threshold of 10%.

[0045] According to a preferred embodiment, there is a step of displaying the nth frame, in which pixels corresponding to the actual Doppler signal are retained, and pixels corresponding to artifacts are modified to present a predetermined grayscale value so that they are not displayed, for example, as Doppler-activated pixels.

[0046] Preferred embodiments of the present invention are described below for illustrative purposes only, with reference to the accompanying drawings. [Brief explanation of the drawing]

[0047] [Figure 1] This figure shows an example of a frame extracted from an ECO-PLA type video. [Figure 2] The image on the left shows an example of identifying colored pixels, and the image on the right shows the corresponding pixel set, such as a binary mask. [Figure 3] The image on the left shows an example of a seed, and the image on the right shows a related frame. [Figure 4] This diagram shows four substeps in step 7 of the method described in the present invention, starting from the top left and proceeding clockwise, allowing us to see the identification of connecting components, the result of the dilation operation, the calculation of colored pixels, and the final result of the frame in the last execution of the steps forming the method described in the present invention. [Figure 5a] On the left is an example of a frame extracted from an ultrasound image, and on the right is a diagram showing the corresponding frame at the end of a noise reduction process performed by the method of the present invention for several types of artifacts detectable in ultrasound images. [Figure 5b]On the left is an example of a frame extracted from an ultrasound image, and on the right is a diagram showing the corresponding frame at the end of a noise reduction process performed by the method of the present invention for several types of artifacts detectable in ultrasound images. [Figure 5c] On the left is an example of a frame extracted from an ultrasound image, and on the right is a diagram showing the corresponding frame at the end of a noise reduction process performed by the method of the present invention for several types of artifacts detectable in ultrasound images. [Figure 5d] On the left is an example of a frame extracted from an ultrasound image, and on the right is a diagram showing the corresponding frame at the end of a noise reduction process performed by the method of the present invention for several types of artifacts detectable in ultrasound images. [Figure 5e] On the left is an example of a frame extracted from an ultrasound image, and on the right is a diagram showing the corresponding frame at the end of a noise reduction process performed by the method of the present invention for several types of artifacts detectable in ultrasound images. [Figure 6] This is a schematic diagram of multiple frames in a temporal sequence, each accompanied by a corresponding image dataset associated with it. [Figure 7] This figure shows an example of temporal fluctuations in signal and noise between two consecutive frames in an ECO-PLA type video. [Figure 8] This is a flowchart for noise reduction processing. [Modes for carrying out the invention]

[0048] Figure 1 shows an example of a single frame belonging to an Echo Doppler image, where it is possible to identify the actual Doppler signal present in the upper left image region, and artifacts present in the dark region of the lower part of the image, in a special case showing a sac.

[0049] In the following description, steps comprising the present invention's method for analyzing ultrasonic images, designed to remove multiple artifacts that may be located on the side of and / or overlap with the actual signal, will be analyzed in more detail, as illustrated by an example of a preferred embodiment. Acquisition of Echo Cupra video input

[0050] In this step, the ultrasound image is acquired in the form of frame packets. Alternatively, the image can be acquired directly, which will be decompressed by the algorithm itself, and the analysis of a single frame performed by the algorithm will result in an evaluation of the temporal duration through the sequence of all frames. Each frame can be represented as a three-dimensional matrix, where the first dimension represents the number of pixels on the vertical axis, the second dimension represents the number of pixels on the horizontal axis, and the third dimension represents the number of color channels (R, G, B). Thus, the value of each pixel can be represented by three values ​​in the range of 0 to 255. Pixel identification including Doppler signals

[0051] To separate the actual Doppler signal from the noise caused by artifacts, it is necessary to separate all Doppler activations (actual + artifact / noise) within each frame of every frame included in the Echo-Doppler image. The result of this operation is two sets of data identification for each frame: the first set contains Doppler activations (colored pixels), and the second set contains the ultrasonic signal (grayscale or monochrome).

[0052] The integration of the two sets identified above for all frames contained in the ultrasound image can be represented through a matrix of dimensions of n frames × n pixels on the vertical axis × n pixels on the horizontal axis.

[0053] In detail, in this step, colored pixels, i.e., pixels containing Doppler activation, are identified for each frame, and then an image dataset is generated in which the pixel value of Doppler activation (colored pixels) is 1, and the pixel value of other pixels is 0 (grayscale or monochrome pixels).

[0054] Colored pixels can be identified by selecting those pixels, for example, by ensuring that the ratio in the difference between the intensities of the R and G, G and B, or R and B channels is greater than or equal to a predetermined threshold. In a preferred embodiment, this threshold is identified as |RG|≧30V|RB|≧30V|GB|≧30. This is due to the demands of computational speed and the fact that it is not necessary to assign known colors to the colored pixels holding the Doppler signal in order to identify which of the two datasets the pixel under inspection belongs to. An alternative method for identification through differences in channel intensity is machine learning.

[0055] According to this method, all pixels in a frame are divided into separate sets by color, and then merged into two subsets: the first containing all the color sets seen in the Doppler spectrum (colored pixels), and the second containing the color sets specific to grayscale ultrasound. The basis of these operations is deductive knowledge of the specific color tones of the Doppler signal representation in a given image. This assignment of color labels for each pixel can be done by uncontrolled clustering operations or with the help of a classifier / neural network trained to extract the dominant colors in the image.

[0056] In one implementation, this method includes the following: (1) Initialize the distribution of pixel values ​​[R, G, B] for uncontrolled clustering methods such as the K-method; (2) Once pixels are grouped into clusters, the centroids are calculated and defined as the three central [R, G, B] of the cluster, i.e., the three most representative of the cluster; (3) Assign a color label to each centroid by extending it to all pixels that make up a given cluster. This can be done by calculating the distance from a known color label to each centroid in the color space (ciede2000) and by selecting a label for one that is at a shorter distance from the centroid. (4) All pixels with color labels other than gray values ​​(including white and black) are set to the same set, and pixels with grayscale values ​​are set to a separate set.

[0057] These methods require computational power, and in the case of uncontrolled methods, they also require deductive knowledge of color to which Doppler signals can be assumed in the image for label generation. As a result, there is a dependency on the color characteristics to which Doppler signals can be assumed in the image. Consequently, there is a dependency on the characteristics and settings of the ultrasonic machine used, as well as a dependency on the type of Doppler, color Doppler, rather than power Doppler. This method is preferable, for example, when color information should be preserved in artifact analysis, for example, when it is possible to observe the filtering of pixels indicating the direction of approach (conventionally red) due to color information.

[0058] An example of the colored pixel identification step using an intensity difference identification method is shown in Figure 2. The left side of the image shows the frame under inspection, while the right side of the image shows the corresponding element belonging to the image dataset, where the bright areas corresponding to the colored pixels identified in the frame and the dark areas of the other pixels can be seen. Identification of activation sequences (second category)

[0059] The basic assumption of the method of the present invention is that Doppler artifacts, or noise, have temporal persistence, i.e., a period shown with respect to a frame, which is shorter than the actual signal for a given Doppler image under inspection.

[0060] To quantify the temporal persistence of the Doppler signal, the "Doppler activation" length of each pixel is calculated through two operations.

[0061] According to the first embodiment, the input to this step is a binary image matrix with size n (frames) × n pixels (vertical axis) × n pixels (horizontal axis) created in the previous step. First, for each pixel of the matrix, an activation vector is added and defined as a vector of length n (equal to the number of frames in the image) and radix [0, 1] (a value that can be estimated by each pixel depending on whether it is colored or not), and each vector element is valued for each frame using the corresponding pixel value of the frame. The above values ​​correspond to 0 for the first category of information, i.e., pixels without Doppler activation, and to 1 or "non-zero" for the second category, i.e., pixels with Doppler activation (Figure 6).

[0062] Subsequently, for each vector, the length of the activation sequence is calculated (a sequence of non-zero values ​​- second category), the sequence (a sequence of non-zero values ​​- second category), and the number of frames for which each sequence sustains (or disappears) the Doppler without interruption.

[0063] Following the application of these operations, the results obtained can be described as a matrix of size n pixels on the vertical axis × n pixels on the horizontal axis, where each matrix element is a persistence vector of different lengths of 1. The reason for the different lengths of the persistence vectors is the number of activations / deactivations detected at each pixel. Each vector element is a pair of values ​​(count, value), where the value refers to the sequence value (0, 1), while the count refers to the number of frames in the sequence.

[0064] For example, if we assume that the video is divided into 20 frames and each pixel has coordinates x1 and y1, then its activation vector is defined as follows, as described above: [0,0,0,1,1,0,0,0,0,0,1,1,1,1,1,0,0,1,1,1]

[0065] The corresponding duration vector will have a length of 6 (equal to the number of activation / inactivation Siemens in the activation vector) and will consist of the following pairs (count, value): [(3, 0), (2, 1), (5, 0), (5, 1), (2, 0), (3, 1)]

[0066] According to another preferred embodiment, the temporal persistence of a Doppler signal can be achieved not by using closed, isolated Doppler regions, i.e., single pixels, but by tracking the temporal persistence of connecting components of Doppler regions within a frame.

[0067] Within a frame belonging to a Doppler-type ultrasound image, it is possible to identify areas defined by grayscale (or monochrome) pixels within the important region targeted by the ultrasound scan. These are called non-Doppler pixels and represent areas defined by Doppler-type signals and colored pixels, while these are called Doppler pixels and represent Doppler-type signals. Doppler pixels can be further distinguished according to Doppler pixels that represent actual Doppler activation, and a standard has been established for actual Doppler pixels and Doppler pixels that represent artifacts, which are called artifact Doppler pixels.

[0068] The connection components signify the identification of distinct objects present in an image. Each of these objects possesses characteristics formed by a set of pixels that satisfy the same adjacency relationship, called a connection.

[0069] It is possible to identify connecting elements within a frame by algorithms that allow for the identification of the characteristics of pixels within a frame, labeling them, and determining whether they belong to a single set representing an object with a specific shape. These algorithms can distinguish pixels belonging to the category of actual Doppler pixels from pixels belonging to the category of artificial Doppler pixels based on criteria such as morphology, centroid distance, or overlap. This family of algorithms allows for the tracking of connecting elements in subsequent frames, once they have been identified in a given frame, so that objects belonging to the same connecting elements can be distinguished even if the objects are deformed or decomposed in subsequent frames. Additionally, when connecting elements are no longer visible in subsequent frames, it is possible to identify them by filtering out artifacts.

[0070] The fundamental need arises from the various implementations in Doppler ultrasound imaging. In fact, in the case of still images, the probe does not move, and important regions and blood vessels always cover the same areas along the sequence of frames that make up the image. In this case, temporal persistence can be determined pixel by pixel by considering the same column and row values ​​in the binary image matrix of all frames in the image. Conversely, when creating images in which the probe moves during the acquisition phase, the background is no longer stationary along the sequence of frames that make up the image, and a single pixel does not maintain a fixed correspondence between consecutive frames. In this case, temporal persistence can be implemented by considering the connecting elements that exist along the sequence of frames that make up the image. Once all Doppler activations are separated (considering both the actual signal and artifact / noise signals), within each frame, connecting elements are identified through two main steps, as described in the previous paragraph, “Pixel Identification Including Doppler Signals.”

[0071] In the first initial setup phase, all connection components (targets) located in the first non-empty binary mask associated with the corresponding video frames are identified.

[0072] For this mask, the centroid of the rectangular bounding box, or the centroid of the connected component itself, is calculated for each identified component, and the vector V0 associated with this component is initialized using the following information: - Unique identifier - The centroid of the rectangular bounding box, or the centroid of the connecting element itself; - Coordinates of the rectangle that defines the constituent elements within the current mask; - A vector V1 of equal size to the total length of the frame, which, when present in the associated frame, will contain a representation of the coordinates (columns and rows) occupied by the object, wherein in its initial setup phase, the coordinates (columns and rows) occupied by the object in the current frame are stored as the first element of the vector V1; - Boolean information B to indicate whether the above subject is still being tracked or not. V0 This is initially set to a value indicating that the target has been tracked. Examples of values ​​include T, 1, Y, S, etc. - A progressive I that indicates the number of frames in which the object is no longer visible after its initial appearance. V0 ; - A vector V2 with dimensions equal to the total length of the frame that makes up the image, where each element indicates the presence of an object in the associated frame.

[0073] Referring to Figure 7, assuming the image on the left represents the first non-empty binary mask, the initialization phase will detect four components within the mask: S(Tn), N1(Tn), N2(Tn), and N3(Tn). For each of these components, a vector is initialized, and vector V S(Tn) , V N1(Tn) , V N2(Tn) , V N3(Tn) This includes the information mentioned above.

[0074] After being identified by the initial setting phase, for each binary mask, instead of tracking the objects present in the current mask Tn+1 as one of the moving objects identified in the previous mask Tn, it is determined whether they represent new objects. To achieve this result, the following is required: - Identifying the positions of all objects in the current mask Tn+1; - For each object belonging to mask Tn+1, determining the center of gravity and the coordinates of the rectangle defining the boundaries of the components inside mask Tn+1; - For all possible pairs of objects between Tn and Tn+1, the distance between the center of gravity and the overlapping values between the coordinates is calculated.

[0075] Once the above values are calculated, if at least one of the following conditions is met, the object OTn+1 located in mask Tn+1 shall be identified as an object in the previous mask Tn. - Object OTn exists, and for the pair (OT n , OT n+1 ), its overlapping value is at least equal to a predetermined threshold Ts; - Object OT n exists, and for the pair (OT n , OT n+1 ), the distance between its coordinates is at most equal to a predetermined threshold Td.

[0076] Using additional comparison criteria, it is possible to determine whether an object in a given mask can be identified as the same object as in the previous mask, such as shape, features, and similarity of overlapping object masks. According to an embodiment of the algorithm, the difference in pixel size between object OT n and object OT n+1 is less than a given threshold Tm.

[0077] If object OT n+1 is identified as the object tracked in the previous mask, the associated initial setting vector is updated using the characteristics of object OT n+1 .: - Vector V1 is the object OT in the current framen+1 It is updated by adding a new element that includes the coordinates (columns and rows) occupied by the existing element; - The vector V2 is updated by setting the value of the element corresponding to the frame under inspection.

[0078] If multiple objects in the current mask Tn+1 are identified as tracked objects in the previous mask, the V1 and V2 vectors are updated using their respective properties.

[0079] Target OT n+1 If a new target is identified as not being tracked by the previous mask, it is tracked by defining a new V0 vector in the same way as described during the initial setup phase.

[0080] No longer a binary mask T n+1 Binary mask T that does not have a match in n+1 The objects present in are also updated as follows: - Boolean property B that identifies whether the target was tracked or not. V0 However, it is set to a value that indicates the subject is no longer being tracked. Examples of values ​​include F, 0, and N. - A progressive I representing the number of frames since the subject appeared, when it is no longer being tracked. V0 This increases by 1. The tolerance threshold allows objects that have been lost for a given number of frames to be tracked for a given number of frames.

[0081] According to the embodiment, the tolerance threshold is defined to consider the possibility that an object may disappear for a given number of frames and then reappear in subsequent frames. This threshold allows the event to be handled such that an object that has temporarily disappeared is considered as a single tracked object rather than two separate objects. According to the embodiment, the threshold is 0, and therefore throughout the frame sequence, the events of object disappearance and reappearance are treated as two separate objects.

[0082] Figure 7 shows an example of the temporal variation of signal and noise between two consecutive frames.

[0083] Previous explanation regarding the identification of connection components, and binary mask T n Assuming that is the first non-empty mask, then: -During the initial setup phase, the four vectors V01, V02, V03, and V04 are defined in association with the target S(Tn), N1(Tn), N2(Tn), and N3(Tn), respectively, and are set to values ​​with the respective characteristics described above; -In the subsequent mask Tn+1, the signal degrades, resulting in the appearance of new targets [S1(Tn+1), S2(Tn+1), S3(Tn+1)] and the disappearance of others [N2(Tn), N3(Tn)] compared to the situation present in Tn; -At the end of tracking the subsequent mask Tn+1: • Objects S1(Tn+1), S2(Tn+1), and S3(Tn+1) will be tracked as connected components of object S(Tn), and their properties will be represented by a vector V initialized at time T. 01 It will be updated; • Object N1(Tn+1) is tracked as a connected component of object N1(Tn), and its properties are determined by a vector V initialized at time T. 02 It will be updated; • Subjects N2(Tn) and N3(Tn) in mask Tn that no longer have a fit in mask Tn+1 will be updated to no longer be tracked. Calculation of reference thresholds in activation sequence length allocation

[0084] The purpose of this step is to obtain the minimum Doppler activation length value, so that activations longer than the threshold are considered to belong to the actual Doppler signal, while activations shorter than the threshold are associated with artifacts. This step begins with the n-pixel dimensional set of the duration (or duration) vector calculated in the previous step. However, in this step, only the activation length value is considered, and the deactivation length is ignored.

[0085] We will proceed using this example, which was defined in the previous step, where we found the persistence vector for pixels x1 and y1 along the 20 frames that make up the ultrasound image: [(3, 0), (2, 1), (5, 0), (5, 1), (2, 0), (3, 1)]

[0086] Only the activation length is considered, i.e. [(2, 1), (5, 1), (3, 1)] And from these, it is possible to define a set of periods as follows: [2, 5, 3].

[0087] This procedure is performed for all pixels that make up the frame. The reference threshold is calculated for all pixels in all activation continuation distributions as (mean + 1.282 × standard deviation) for the 90th percentile value or normal distribution.

[0088] The threshold is calculated as the 98th percentile if the activation length of the smallest Doppler signal is obtained by tracking the temporal duration of closed, isolated Doppler regions within the frame. Identifying the actual location of the signal seed

[0089] Once a threshold is defined, the actual signal coordinate position needs to be determined for each frame, and the Doppler signal is seeded.

[0090] For each pixel, a pre-calculated activation sequence vector is taken into consideration and compared to a threshold as follows: - Sequences with a length exceeding a threshold are kept "on" (in activation vector sequences, the pixel value remains 1); - Shorter activation sequences are "switched off" (by setting the pixel value to 0 for those activation sequences, even though the threshold for "switching on" is not reached and they are considered to belong to the actual Doppler signal).

[0091] Proceeding with the example defined in the previous step, and assuming that the previous step occurred at a threshold of 4, the activation length vector for pixels x1 and y1 is: We can confirm [(3,0), (2,1), (5,0), (5,1), (2,0), (3,1)]. Only the following sequences have a continuation that exceeds the threshold: (5, 1).

[0092] Therefore, the activation vector is initially calculated as follows: [0,0,0,1,1,0,0,0,0,0,1,1,1,1,1,0,0,1,1,1] At the end of this step, the result will be as defined above. [0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,0,0,0,0,0].

[0093] At the end of this step, a matrix will be generated with dimensions of n frames × n pixels on the vertical axis and n pixels on the horizontal axis. Here, for each frame, each pixel will have a value of 1 if the associated activation sequence is at least equal to a threshold, and 0 otherwise. Identifying the actual signal location

[0094] This step takes into account a matrix of dimensions of n frames × n pixels on the vertical axis and n pixels on the horizontal axis, as defined in the pixel identification step that includes the Doppler signal. Here, for each frame of the matrix, the actual Doppler signal in the corresponding original frame of the video is separated (a denoised frame is obtained) using the image dataset defined in the previous step, such as a binary mask. For each frame, the [R, G, B] values ​​of the pixels belonging to the original frame are then preserved as the actual signal by taking the three relevant [R, G, B] channel values ​​corresponding to each coordinate in the original frame, considering the pixel coordinates with a value of 1 in the image dataset obtained in the previous step, or by multiplying by the binary mask up to the original frame.

[0095] A possible denoising frame view is that the Doppler signal pixels retain the [R, G, B] channel values ​​of the original frame and are therefore colored, while the other pixels assume the average value of the three [R, G, B] channel pixels in the original frame and are therefore grayscale.

[0096] Figure 3 shows an example of a signal seed (second category pixel) for a given frame obtained by saving only the activation sequences of pixels with a length equal to at least the threshold (Figure 3a) on the left. Figure 3b shows a frame obtained by applying the seed to the original frame in Figure 1. It is possible to see the colored pixels corresponding to the actual signal, while artifacts are masked in grayscale. Segment decomposition and expansion

[0097] The previous step of locating the actual signal returns pixel coordinates containing the actual Doppler signal. However, the returned actual signal region may not be accurately assumed and may be underestimated. For example, removing artifacts superimposed on the signal may also result in the removal of pixels belonging to the actual signal. Other possible examples of pixels belonging to the removed signal include pixels around blood vessels being thinned in subsequence frames due to probe movement by the operator, or due to oscillations in the signal itself.

[0098] To accurately identify the edges of the signal region identified in the seed positioning process, starting from each seed (such as the signal region / coordinates recognized at the end of the previous step), connected Doppler activation components from the original video that overlap each seed are attached to the signal region / coordinates. This is performed for each frame according to the following operation of the first segment decomposition algorithm. For example: a. The corresponding activated image dataset is selected, where a value of 1 indicates a colored pixel and a value of 0 indicates something else (results obtained in step 2 of this method). Note that from this image dataset, pixels corresponding to artifacts have not yet been "switched off" (Figure 2b). b. Morphological collapse behavior is applied to the image dataset to separate the main connecting components of the actual signal from artifact-related components as much as possible. In this embodiment, collapse is applied using a circular kernel with a 5-pixel diameter (Figure 4a). c. In the image dataset at point b, connected elements are identified and labeled. That is, each macroscopic element formed by adjacent and connected n pixels > 1 (Figure 4a - connected elements with their colors). Connected elements with associated colors. d. The connecting components in Figure 4a are selected based on an image dataset in which only activation sequences of pixels with a length equal to at least a threshold are stored (second category of pixels - Figure 3a). For example, based on data obtained by overlaying the data from Figure 3a onto Figure 4e, only the brightest regions indicated by arrows are retained, along with the percentage of pixels overlaid with the activation pixels from Figure 3a that are above a predetermined threshold. e. For each connected component selected based on the pixels of the second category (Figure 3a), a morphological expansion operation is applied to the same extent as the previous collapse application (and in reverse, to restore the collapse). In this embodiment (Figure 4b), the expansion is applied using a circular kernel with a 5-pixel diameter.

[0099] When the enlarged image needs to be displayed on the screen, the following steps are performed, starting from Figure 4b: For each extended connection component, count how many pixels remain after filtering (pixels containing the actual signal from the previous step); In the new denoised frame, if the number of colored pixels is greater than or equal to a given threshold, the color is extended to all components, and the pixels are assigned the corresponding [R, G, B] values ​​from the original frame. In this embodiment, the optimal threshold is set to 10%.

[0100] According to an alternative embodiment, the region in Figure 3a is extended as follows (not shown): The image data in Figure 3a (second category) is overlaid on the associated original frame (Figure 1) to define the seed; The segment decomposition algorithm is applied to the original frame based on a seed, such as the region that grows to identify the largest region, which is the seed in Figure 3a.

[0101] This alternative method also allows for obtaining a larger area representing the Doppler signal in the original frame, because artifacts, particularly motion artifacts, have been removed beforehand, and the underestimated Doppler region of the extension has been identified, as shown in Figure 3a.

[0102] Finally, Figure 5 shows several embodiments of the method described in the present invention for removing artifacts from an ultrasound image frame by placing the original frame side by side with its denoised frame.

[0103] The left side shows the frame before denoising was applied, and the right side shows it afterward. Note that colored pixels, considered "artifacts," have been replaced with framescale pixels. These five images were selected because they contain different types of artifacts, they depict different anatomical objects internally, and they have different characteristics. In its original size, each frame is 800 x 566 pixels. Overall, colored pixels have been identified to show intensity differences between R, G, and B channels at least 30 points. The selected activation threshold is 90% of the total activation length for each pixel in the image.

[0104] Figure 5a: Example of flash artifact removal. This video consists of 113 frames and was recorded at 57 fps; the algorithm sets a receptive threshold for 17 consecutive activation frames.

[0105] Figure 5b: Example of flash artifact removal. This video consists of 70 frames and is recorded at 13 fps; the algorithm sets a receptive threshold for 11 consecutive activated frames.

[0106] Figure 5c: Example of flash artifact removal. This video consists of 138 frames and was recorded at 57 fps; the algorithm sets a receptive threshold for 16 consecutive activation frames.

[0107] Figure 5d: Example of pseudo-flow artifact removal. This video consists of 499 frames and was recorded at 57 fps; the algorithm sets a reception threshold for 19 consecutive activation frames.

[0108] Figure 5e: An example of random and extreme artifact removal. This video consists of 229 frames and was recorded at 57 fps; the algorithm sets a receptive threshold for 14 consecutive activation frames.

[0109] According to the present invention, n-frames with corresponding pixels are superimposed, and the j-th pixel at position p×q (height×width) corresponds to the relevant pixel. This is because each frame has the same dimensions, and is based on, for example, a machine that performs ultrasound imaging, and its stability is high.

[0110] Figure 8 shows a flowchart of the artifact removal method for Doppler images, where the aforementioned step sequence can be seen, which makes it possible to obtain images with no artifacts or with significantly reduced artifacts.

Claims

1. A computer-based noise reduction method for removing artifacts from Doppler ultrasound images, The steps include: acquiring n frames that make up the ultrasound image, For each of the n frames, the steps include identifying a first pixel category representing a first pixel in the nth frame containing an ultrasonic signal, and a second pixel category representing a second pixel in the nth frame containing Doppler activation. For each of the n frames, the steps include: generating a first set of image data to associate at least a first or second category with the position of each pixel in the corresponding frame, so as to classify the pixels of the corresponding frame and identify pixels containing Doppler activation for each frame; A step of generating at least a persistent sequence relating to a frame portion containing identification data ordered according to a frame sequence such that adjacent identification data in the sequence point to consecutive frames, wherein each identification data is assumed to have a first value if the corresponding frame portion contains a Doppler-activated pixel, or a second value different from the first value if the corresponding frame portion does not contain a Doppler-activated pixel; The steps include calculating the length of each persistence subsequence, which comprises consecutive identification data having the first value, in the persistence sequence associated with the frame portion, A step of calculating a reference threshold based on the length distribution of the persistence subsequences, which refers to the identification data of all n frames, For each of the n frames, If frame n corresponds to one of the persistence subsequences having a length longer than the reference threshold so as to represent a pixel with an actual Doppler signal, then the information that the position of the pixel belongs to the second category is, or If frame n corresponds to another of the persistence subsequences having a length shorter than the reference threshold to represent an artifact, then the information that the pixel position belongs to the first category is The steps include identifying a second set of image data (SEED) in order to associate it with each pixel position of the corresponding frame, Noise reduction methods, including those mentioned above.

2. The noise reduction method according to claim 1, wherein the frame portion is a pixel of the frame, and the generation step includes the step of generating a plurality of sustained sequences for each pixel of a plurality of frame formats.

3. The denoising method according to claim 1, comprising the step of searching for connected elements having pixels of the first category in frame n and frame n+1, and executing an algorithm that defines the frame portion of the frame, wherein the step of generating a persistence sequence in order to obtain a trace of the connected elements between frame n and frame n+1 includes the step of comparing a first parameter in a first connected element of frame n and a second parameter in a second connected element of frame n+1 to associate the second connected element with a first persistence sequence in the first connected element, or generating a second persistence sequence in the second connected element based on the comparison step, wherein the parameter is a centroid and / or a parameter representing the overlap of the first and second connected elements and / or a parameter representing the similarity of the shapes of the first and second connected elements and / or a parameter representing the size or dimensions of the first and second connected elements.

4. The noise reduction method according to claim 2, wherein the step of identifying the second set (SEED) is performed based on a connection component search algorithm to identify a region representing the actual Doppler signal, and includes the step of performing a segment decomposition algorithm based on the search for connection components in pixels of the second category in at least one nth frame of the ultrasound image, comprising the connection components of the second category and connection components representing the actual Doppler signal based on overlap criteria, and deactivating the connection components of the second category that do not satisfy the overlap criteria in order to obtain a region enhancement effect of the connection components representing the actual Doppler signal.

5. The noise reduction method according to claim 1, wherein the step of identifying the second pixel category is performed for each Doppler active color present in the frame according to a plurality of subcategories.

6. The denoising method according to claim 1, wherein the step of identifying the second category of pixels is carried out by clustering.

7. The noise reduction method according to claim 1, wherein the step of identifying the second category of pixels is performed by selecting via the intensity differences of the R, G, and B channels.

8. The denoising method according to claim 2, wherein the step of calculating a threshold is performed by calculating the 90th percentile in the length distribution of the subsequence, and the continuous identification data of the subsequence is based on the sum of means of subsequences having a first representative Doppler activation value or having continuous identification data of subsequences having a first representative Doppler activation value, with a sum of means of twice the standard deviation in length.

9. The denoising method according to claim 3, wherein the step of calculating the threshold is by calculating the 98th percentile in the activation length distribution, or based on the sum of means of the persistent subsequences having a first representative value of Doppler activation, with a length of twice the standard deviation.

10. The noise reduction method according to claim 1, wherein the calculation of the reference threshold is performed by color labeling of each pixel obtained from the previous step and by using a classifier trained to extract the main colors seen in the image.

11. The noise reduction method according to claim 1, wherein the noise-reduced frame has R, G, and B channel values ​​equal to the values ​​of the original frame when color is represented, and is an average value when grayscale is represented.

12. The denoising method according to claim 4, wherein the region expansion is performed by a morphological collapse method and by expansion with a kernel having a diameter or side length of 5 pixels and a threshold of 10%.

13. The noise reduction method according to claim 1, wherein the frame can be represented as a three-dimensional matrix.

14. The method according to any one of claims 1 to 13, comprising the step of displaying the nth frame, where pixels corresponding to the actual Doppler signal are retained and pixels corresponding to artifacts are modified.

Images