CN111265246B - Ultrasonic color imaging processing method and device - Google Patents
Ultrasonic color imaging processing method and device Download PDFInfo
- Publication number
- CN111265246B CN111265246B CN202010097056.9A CN202010097056A CN111265246B CN 111265246 B CN111265246 B CN 111265246B CN 202010097056 A CN202010097056 A CN 202010097056A CN 111265246 B CN111265246 B CN 111265246B
- Authority
- CN
- China
- Prior art keywords
- color
- data
- blood flow
- image data
- black
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/488—Diagnostic techniques involving Doppler signals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
Abstract
The disclosure relates to the technical field of ultrasonic imaging, and particularly provides an ultrasonic color imaging processing method and device. The method comprises the following steps: acquiring color echo signals, and acquiring color characteristic data of each sampling position according to the color echo signals; performing DSC digital scanning transformation according to the color characteristic data to obtain color image data of each pixel point of the image display area; determining data of pixel points meeting a preset threshold condition as blood flow data, wherein the preset threshold condition is determined according to the color image data; and obtaining blood flow velocity data of corresponding pixel points according to the blood flow data, and carrying out color coding imaging on the image according to the blood flow velocity data. According to the method, DSC conversion is carried out on all color characteristic data parameters of the sampling position, threshold judgment is carried out according to the parameters after space conversion, the speed estimation matrix is refined from the sampling position to the pixel point, threshold judgment is more comprehensive and accurate, and imaging accuracy is higher.
Description
Technical Field
The disclosure relates to the technical field of ultrasonic imaging, in particular to an ultrasonic color imaging processing method and device.
Background
Ultrasonic imaging has become one of the most widely used medical devices in clinical practice because of its advantages such as non-invasiveness, real-time performance, and convenient operation, and color doppler imaging is more significant for ultrasonic imaging. The color Doppler imaging superimposes color blood flow on a black-and-white image, has a two-dimensional ultrasonic structural image, provides dynamic information of the blood flow, and has a great effect on diagnosis of vascular diseases.
In ultrasound imaging, the moving medium on which the imaging effective information depends is blood flow, but other movements than blood flow may also produce doppler phenomena, such as tissue movement and the like. Due to the influence of factors such as respiration, scattering, electronic components and the like, echo signals received by the actual transducer contain various data, such as common noise, blood flow movement, non-blood flow movement and the like. And common noise and non-blood flow motion are noise which is not needed for imaging, so that the filtering of the common noise and the non-blood flow motion, and the maximum representation of the blood flow motion are the key points of color signal processing.
Disclosure of Invention
In order to improve the accuracy of ultrasonic color imaging, the disclosure provides an ultrasonic color imaging processing method and device.
In a first aspect, the present disclosure provides an ultrasound color imaging processing method, including:
acquiring color echo signals, and acquiring color characteristic data of each sampling position according to the color echo signals;
performing DSC digital scanning conversion according to the color characteristic data to obtain color image data of each pixel point of the image display area;
determining data of pixel points meeting a preset threshold condition as blood flow data, wherein the preset threshold condition is determined according to the color image data;
and obtaining blood flow velocity data of corresponding pixel points according to the blood flow data, and carrying out color coding and display on the image according to the blood flow velocity data.
In some embodiments, the determining that the data of the pixel points satisfying the preset threshold condition is blood flow data includes:
dividing an image display area into a plurality of sub-areas, wherein each sub-area comprises a plurality of pixel points;
according to the color image data of the pixel points in the sub-area, calculating to obtain the color image data of the sub-area;
determining the sub-area meeting the preset threshold condition as a color area;
and determining the data of the pixel points meeting the preset threshold condition as blood flow data according to the color image data of each pixel point in the color area.
In some embodiments, the method further comprises:
acquiring black and white echo signals, and processing the black and white echo signals to obtain black and white characteristic data of each sampling position;
performing DSC digital scanning conversion according to the black-and-white characteristic data to obtain black-and-white image data of each pixel point of an image display area, wherein the black-and-white image data comprises gray values;
dividing an image display area into a plurality of sub-areas, wherein each sub-area comprises a plurality of pixel points;
calculating black-and-white image data of the sub-region according to the black-and-white image data of the pixel points in the sub-region;
determining the sub-area meeting the preset threshold condition as a color area; the preset threshold condition is determined according to the color image data and the black and white image data;
and determining the data meeting the preset threshold condition as blood flow data according to the color image data and the black and white image data of each pixel point in the color area.
In some embodiments, before the obtaining the color feature data of each sampling position according to the color echo signals, the method further includes:
carrying out wall filtering processing on the color echo signals;
the color image data comprises a first energy value before wall filtering of each pixel point, a second energy value after wall filtering of each pixel point and a variance value.
In some embodiments, the preset threshold condition comprises:
the first energy value is less than or equal to a first preset energy threshold value;
the ratio of the first energy value to the second energy value is less than or equal to a preset energy ratio threshold;
the second energy value is greater than or equal to a second preset energy threshold;
the variance value is less than or equal to a preset variance threshold;
the gray value is less than or equal to a preset gray value;
the determining that the data of the pixel points meeting the preset threshold condition is blood flow data includes:
and when the data of the pixel point meets all conditions, determining that the data of the pixel point meets the preset threshold condition.
In some embodiments, the obtaining color feature data for each sampling location from the color echo signals includes:
carrying out autocorrelation processing on the color echo signals subjected to wall filtering processing to obtain the characteristic intermediate parameters of each sampling position;
obtaining the first energy value for each sample position before wall filtering; the color feature data comprises the feature intermediate parameters and the first energy value for each sampling position before the wall filtering.
In some embodiments, the performing DSC digital scan conversion according to the color feature data to obtain color image data of each pixel point of an image display area includes:
and calculating to obtain a second energy value and a variance value of each pixel point according to the characteristic intermediate parameters after DSC digital scanning conversion.
In some embodiments, after said color-coded imaging of an image from said blood flow velocity data, further comprising:
obtaining blood flow speed data of a previous frame image of a current frame image;
judging whether the difference value of the blood flow speed data of the current frame image and the blood flow speed data of the previous frame image is larger than a first preset threshold value or not;
if not, time correlation processing is carried out on the current frame image.
In a second aspect, the present disclosure provides an ultrasound color imaging processing apparatus comprising:
the acquisition module is used for acquiring color echo signals and acquiring color characteristic data of each sampling position according to the color echo signals;
the space transformation module is used for carrying out DSC digital scanning transformation according to the color characteristic data to obtain color image data of each pixel point of the image display area;
the determining module is used for determining the data of the pixel points meeting a preset threshold condition as blood flow data, wherein the preset threshold condition is determined according to the color image data;
and the coding module is used for obtaining blood flow velocity data of corresponding pixel points according to the blood flow data and carrying out color coding imaging on the image according to the blood flow velocity data.
In some embodiments, the determining module is specifically configured to:
dividing an image display area into a plurality of sub-areas, wherein each sub-area comprises a plurality of pixel points;
according to the color image data of the pixel points in the sub-area, the color image data of the sub-area is obtained through calculation;
determining the sub-area meeting the preset threshold condition as a color area;
and determining the data of the pixel points meeting the preset threshold condition as blood flow data according to the color image data of each pixel point in the color area.
In some embodiments, the obtaining module is further configured to obtain a black-and-white echo signal, and process the black-and-white echo signal to obtain black-and-white feature data of each sampling location;
the space transformation module is further used for carrying out DSC digital scanning transformation according to the black-and-white characteristic data to obtain black-and-white image data of each pixel point of an image display area, and the black-and-white image data comprises a gray value;
the determining module is further configured to divide the image display area into a plurality of sub-areas, where each sub-area includes a plurality of pixel points;
calculating black-and-white image data of the sub-region according to the black-and-white image data of the pixel points in the sub-region;
determining the sub-area meeting the preset threshold condition as a color area; the preset threshold condition is determined according to the color image data and the black and white image data;
and determining the data meeting the preset threshold condition as blood flow data according to the color image data and the black-and-white image data of each pixel point in the color area.
In some embodiments, the apparatus further comprises:
the wall filtering module is used for carrying out wall filtering processing on the color echo signals; the color image data comprises a first energy value before wall filtering of each pixel point, a second energy value after wall filtering of each pixel point and a variance value;
the autocorrelation module is used for performing autocorrelation processing on the color echo signals after wall filtering processing to obtain the characteristic intermediate parameters of each sampling position; the color feature data includes the feature intermediate parameters and the first energy value for each sampling position before the wall filtering.
In a third aspect, the present disclosure provides an ultrasonic color imaging apparatus comprising:
a processor; and
a memory communicatively coupled to the processor and storing computer readable instructions executable by the processor, the processor performing the method according to any of the embodiments of the first aspect when the computer readable instructions are read.
In a fourth aspect, the present disclosure provides a storage medium storing computer-readable instructions for causing a computer to perform the method according to any one of the embodiments of the first aspect.
The ultrasonic color imaging processing method in the embodiment of the disclosure comprises the steps of acquiring color echo signals, obtaining color feature data of each sampling position according to the color echo signals, inputting the color feature data into a DSC (differential scanning calorimetry) for conversion to obtain color image data of each pixel point, then performing threshold judgment on the color image data of each pixel point, and determining data meeting a preset threshold condition as blood flow data, thereby removing data interference representing tissue motion or noise and improving imaging accuracy. And obtaining the blood flow velocity estimation value of the corresponding pixel point according to the screened blood flow data, and further carrying out color coding and display to obtain the ultrasonic Doppler image. According to the method, DSC conversion is performed on all color characteristic data parameters of the sampling position, threshold judgment is performed according to the parameters after space conversion, the speed estimation matrix is refined from the sampling position to pixel points, the threshold judgment is more comprehensive and accurate, and the noise removal effect is better.
In the ultrasonic color imaging processing method in the embodiment of the disclosure, when the threshold is determined, the image display area is firstly divided into a plurality of sub-areas, each sub-area includes a plurality of pixel points, then data representing the sub-area is obtained by calculation according to the color image data of all the pixel points in the sub-area, then the threshold is determined according to the data of the sub-area, and the sub-area meeting the preset threshold condition is determined to be a color area, so that the noise of non-blood flow can be roughly screened out. And then, for the color area, the threshold judgment of pixel points one by one is carried out on the area to obtain more accurate blood flow data, and meanwhile, the calculated amount is greatly reduced by dividing the area and combining the judgment of pixel points one by one.
Drawings
In order to more clearly illustrate the embodiments of the present disclosure or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present disclosure, and other drawings can be obtained by those skilled in the art without creative efforts.
Fig. 1 is a schematic diagram of color signal processing in the prior art.
Fig. 2 is a flow chart of an ultrasound color imaging processing method in some embodiments according to the present disclosure.
FIG. 3 is a flow chart of an ultrasound color imaging processing method according to further embodiments of the present disclosure.
FIG. 4 is a flow chart of a method of ultrasound color imaging processing in accordance with one embodiment of the present disclosure.
Fig. 5 is a block diagram of an ultrasound color imaging processing apparatus according to some embodiments of the present disclosure.
FIG. 6 is a schematic diagram of a computer system suitable for use in practicing the disclosed methods.
Detailed Description
The technical solutions of the present disclosure will be described clearly and completely with reference to the accompanying drawings, and it is to be understood that the described embodiments are only some embodiments of the present disclosure, but not all embodiments. All other embodiments, which can be derived by one of ordinary skill in the art from the embodiments disclosed herein without making any creative effort, shall fall within the protection scope of the present disclosure. In addition, technical features involved in different embodiments of the present disclosure described below may be combined with each other as long as they do not conflict with each other.
The ultrasonic color Doppler imaging is to superpose color blood flow on a black-and-white image, wherein black-and-white signals can be imaged according to the intensity of reflection signals at different depths of tissues, the higher the intensity of the reflection signals is, the larger the gray value of the image is, and otherwise, the smaller the gray value of the image is. The color image signal calculates the frequency shift of the ultrasonic signal caused by the Doppler effect, if the medium moves towards the transducer, the echo frequency will be increased, if the medium moves away from the transducer, the echo frequency will be decreased, and the blood flow velocity is obtained by the frequency shift and is represented by different colors (such as red and blue) and shades, so as to obtain a real-time color blood flow image. And fusing the color blood flow image and the black and white image together for display to form a color Doppler image.
In the prior art, when the ultrasound imaging system performs color signal processing, the processing flow may be as shown in fig. 1, for example. The color echo signals enter an autocorrelation module after being processed by demodulation, wall filtering and the like, the autocorrelation module obtains three intermediate parameters A, B and C of color characteristics of the sampling position through autocorrelation operation, and parameters such as speed, energy or variance of the sampling position can be obtained through calculation according to the intermediate parameters A, B and C. The threshold judgment module is used for carrying out threshold judgment on the data of the sampling position, the data of the sampling point meeting the threshold condition can be identified as blood flow data, the data of the sampling point not meeting the threshold condition represents noise, and the noise such as tissue motion and the like can be removed through the threshold judgment. Calculating according to the screened blood flow data to obtain blood flow velocity data of the sampling position, performing two-dimensional smoothing processing, entering a DSC (Digital scanning Converter) for spatial transformation, converting the sampling point data into pixel point data of an image display area, and finally performing frame averaging processing, color coding and the like to obtain a color image for display.
In the scheme in the prior art, the threshold judgment module performs screening, denoising and speed estimation according to the color characteristic value of the sampling point, and because the sampling point data amount is small and the fineness is low, the accuracy of screening and denoising the noise is poor, and after the obtained blood flow velocity data of the sampling point enters DSC transformation, the data accuracy of each pixel point in the display area is further reduced, so that the imaging accuracy is low.
Based on the above-mentioned defects in the prior art, the present disclosure provides an ultrasound color imaging processing method, so as to improve the accuracy of ultrasound color imaging. A method in some embodiments according to the present disclosure is illustrated in fig. 2.
As shown in fig. 2, in some embodiments, the present disclosure provides an ultrasound color imaging processing method comprising:
s201, acquiring color echo signals, and acquiring color characteristic data of each sampling position according to the color echo signals.
Specifically, in COLOR ultrasound imaging, the received digital signals of each channel may be subjected to a beamformer, filtering, quadrature demodulation, and the like, and finally a complex signal including a black-and-white echo signal (BW) and a COLOR echo signal (COLOR) is obtained. After the color echo signals are acquired, color feature data of each sampling position can be obtained through calculation according to the color echo signals.
For example, in one embodiment, three characteristic intermediate parameters a, B, C for each sampling position may be obtained by an autocorrelation algorithm, and the characteristic intermediate parameters are used to calculate speed values, energy values, and variance values of the sampling points, which are described in detail below. The characteristic intermediate parameter is the color characteristic data of the sampling position.
In some embodiments, before performing the autocorrelation operation on the color echo signals, the signals may be first subjected to a wall filtering process, so that most of the low-frequency tissue signals can be filtered out, and the interference can be reduced.
S202, DSC digital scanning conversion is carried out according to the color characteristic data to obtain color image data of each pixel point of the image display area.
Specifically, the DSC (Digital Scan Converter) may convert the data of the sampling position into data values of the entire display area through coordinate conversion and interpolation according to the size and shape of the display area. The data processing process of the DSC refers to the DSC operating principle in the related art, and the DSC is not improved in the present disclosure, and is not described in detail herein.
And carrying out DSC digital scanning conversion on the color characteristic data of each sampling position to obtain color image data of each pixel point of the whole image display area. The color feature data may include one or more of feature intermediate parameters, black and white feature data, wall filtering initial energy, and other parameters, and after DSC conversion, the color image data corresponding to each pixel point of the image is obtained through calculation, for example, the energy value, variance value, black and white gray value, initial energy value, and the like of each pixel point.
In an exemplary embodiment, the color feature data includes feature intermediate parameters a, B, and C of each sampling location, and the DSC transform is performed, and an energy value and a variance value of each pixel point are obtained through calculation according to the feature intermediate parameters a, B, and C, and the energy value and the variance value of each pixel point are the color image data of the point.
S203, determining the data of the pixel points meeting the preset threshold condition as blood flow data, wherein the preset threshold condition is determined according to the color image data.
Specifically, in step S202, color image data for each pixel in the display region is calculated, and it is necessary to screen out data representing blood flow motion by removing noise data such as tissue motion from these data.
First, a preset threshold condition is determined according to color image data, and the condition is used as a judgment condition for judging the threshold of the data, for example: confirming that the energy value after wall filtering is smaller than a preset threshold value, namely confirming that the noise data does not belong to blood flow movement, and then confirming that the preset threshold value condition is that the energy value is larger than or equal to the preset energy threshold value; for another example: the larger the variance value of the pixel point is, the stronger the randomness of the representative signal is, and the more likely the representative signal is noise, the condition that the variance value is smaller than or equal to the preset variance threshold value can be determined as the preset threshold value. It is to be understood that the determination of the preset threshold condition can be based on various parameters and combinations thereof, and is not limited to the above examples, and will be described in detail hereinafter, which is not detailed herein.
According to the determined preset threshold condition, threshold judgment can be carried out on data of each pixel point in the display area, the preset threshold condition is met, the data are represented as blood flow data, the corresponding points are color areas, and on the contrary, the corresponding points are noise data, and removal is carried out.
And S204, obtaining blood flow velocity data of the corresponding pixel points according to the blood flow data, and carrying out color coding imaging on the image according to the blood flow velocity data.
Specifically, in step S203, the noise data is subjected to threshold value determination and filtered, and for the pixel points determined to be blood flow motion, blood flow velocity data indicating an estimated value of the velocity of the blood flow motion at the position is calculated from the blood flow data at the pixel points. And carrying out color coding imaging on the image according to the blood flow velocity data, and further superposing the image on a black-and-white image for displaying.
In view of the above, in the ultrasonic color imaging processing method in some embodiments of the present disclosure, the color feature data of the sampling point is input to the DSC to perform spatial transformation, so as to obtain the color image data of each pixel point. And then carry out threshold value judgement according to the color image data of pixel, rather than carrying out threshold value judgement according to the color characteristic data of sampling point, the blood flow velocity estimation matrix is the pixel by sampling position is meticulous, and threshold value judgement is more comprehensive meticulous, and the noise removal effect is better.
In some embodiments, as can be seen from the foregoing, the present disclosure refines the threshold into the pixel point data of the display area, greatly improves the denoising effect, and correspondingly increases the computation amount for threshold analysis and calculation. In order to reduce the amount of computation of data analysis processing while improving the imaging accuracy, as shown in fig. 3, in step S203, the disclosed method includes:
s301, dividing the image display area into a plurality of sub-areas, wherein each sub-area comprises a plurality of pixel points.
Specifically, the display area is firstly divided into a plurality of sub-areas, and the division of the sub-areas may be determined according to the size and shape of the display area, for example, in one embodiment, each sub-area is a pixel area of 6 × 6, that is, each sub-area contains 36 pixel points in total. It is understood that, for the division of the sub-area, one skilled in the art can select the sub-area according to the size, shape, operation precision, etc. of the specific display area, and the above-mentioned exemplary embodiments do not limit the present disclosure.
S302, according to the color image data of the pixel points in the sub-area, the color image data of the sub-area is obtained through calculation.
And calculating color image data representing the sub-region data according to the color image data of the pixel points in each sub-region.
In an exemplary embodiment, the average value of the color image data of all the pixels in the sub-region may be used as the color image data of the sub-region. For example, the energy average and variance average of all pixel points in the sub-region are calculated, and the average represents the energy value and variance value of the sub-region.
And S303, determining the sub-area meeting the preset threshold condition as a color area.
Specifically, the method in step S203 may be referred to determine the preset threshold condition, which is not described herein again. And judging the threshold of each sub-area according to a preset threshold condition, and when the preset threshold condition is met, indicating that the sub-area contains color, namely determining that the sub-area is a color area. On the contrary, if the condition is not satisfied with the preset threshold, it indicates that no color is contained in the sub-region, and the color value is 0.
By carrying out threshold judgment on the divided sub-regions, the noise of the non-blood flow region can be quickly removed, and the sub-region containing blood flow, namely the color region, is obtained.
S304, according to the color image data of each pixel point in the color area, determining the data of the pixel points meeting the preset threshold condition as blood flow data.
For the color region obtained after the screening, the pixel points in the region are subjected to threshold judgment one by one, and the judgment process is only referred to the step S203, which is not described herein again. And removing the noise data which do not meet the preset threshold condition in the color area, thereby obtaining accurate blood flow data.
In the embodiment, the threshold value judgment is carried out by dividing the region and combining the pixel points one by one, so that the imaging precision is improved, and the calculation amount of data analysis and calculation is greatly reduced.
Fig. 4 shows an ultrasound color imaging method in an embodiment of the present disclosure, and the ultrasound color imaging method of the present disclosure is specifically described in conjunction with the embodiment of fig. 4. As shown in fig. 4, in the present embodiment, the method of the present disclosure includes:
and S01, beam synthesis line data. For example, the ultrasonic probe receives digital echo signals of respective channels, and synthesizes line data by using a beamformer.
S02, the data combined by the beamformer is subjected to, for example, dynamic filtering and quadrature demodulation to obtain a complex signal including a black-and-white echo signal (BW) and a COLOR echo signal (COLOR).
And S03, processing the black-and-white echo signal (BW) to obtain black-and-white characteristic data of each sampling position.
In this embodiment, the black-and-white echo signal enters the black-and-white signal processing module to be processed, and then undergoes envelope detection, compression conversion and spatial preprocessing in sequence to obtain black-and-white feature data at each sampling position, where the black-and-white feature data includes, for example, a black-and-white gray scale value. The black and white signal processing can be realized by the method in the related art, and the details of the disclosure are not repeated.
And S04, performing wall filtering processing on the color echo signals, filtering most low-frequency tissue signals, and reducing signal interference.
And S05, acquiring a first energy value of each sampling position before wall filtering. The first energy value represents the pre-filter energy D for each sample position in the color echo signal.
And S06, performing autocorrelation processing on the color echo signals after wall filtering processing to obtain characteristic intermediate parameters of each sampling position.
Specifically, let the number of times of repeating transmission for each scanning process be N, if I i,j,k And Q i,j,k Respectively representing the spatial sampling sequence values of the orthogonal signal after the K-time emission of the jth sampling point on the ith scanning line in the scanning plane after the orthogonal demodulation, and then the initial energy value D of the sampling position i,j The expression of (a) is as follows:
if with X i,j,k And Y i,j,k Respectively representing the values of the sampling sequence after wall filtering, and let A i,j And B i,j Representing the real and imaginary parts, respectively, of their estimated autocorrelation function by C i,j Representing an estimate of power, the expression is:
through the above autocorrelation operation, three characteristic intermediate parameters a, B, and C of each sampling position are obtained, the characteristic intermediate parameters are used to calculate the velocity value, the energy value, and the variance value of the sampling point, and the calculation manner is described in detail below.
And S07, performing two-dimensional smoothing processing on the characteristic intermediate parameter and the first energy value. The two-dimensional smoothing carries out some related smoothing processing on the data in the scanning line direction, and the data precision is improved.
And S08, carrying out DSC digital scanning conversion according to the characteristic intermediate parameter, the first energy value and the black-and-white image data to obtain the full characteristic parameter of each pixel point of the image display area.
Specifically, the DSC can convert data of the sampling positions into data values of the entire display area through coordinate conversion and interpolation according to the size and shape of the display area. In this embodiment, spatial transformation is performed jointly with the characteristic intermediate parameters (a, B, and C), the first energy value D, and the black-and-white image data to obtain a full characteristic parameter of each pixel point in the display area, including the first energy value of each pixel point, the black-and-white image data, and the characteristic intermediate parameters.
And S09, determining the data of the pixel points meeting the preset threshold condition as blood flow data.
Specifically, in the present embodiment, the data threshold determination may be divided into three parts:
1. and calculating the pixel point parameters one by one.
In this embodiment, a velocity value V, a variance value σ, and a power value P of each pixel are calculated according to the characteristic intermediate parameters a, B, and C, where the power value P represents a second energy value after wall filtering of the pixel. The specific calculation formula is as follows:
V i,j ＝Farctan(-B i,j /A i,j )
P i,j ＝C i,j
wherein T represents the time interval between two scanning lines, theta is the included angle between the scanning line and blood flow, and omega 0 Is a constant.
And calculating to obtain a variance value sigma and a second energy value P of each pixel point, and combining the first energy value D and the black-and-white image data to jointly serve as parameters for subsequent threshold judgment, wherein each parameter contributes to noise suppression differently.
2. A preset threshold condition is determined.
In this embodiment, the preset threshold condition includes:
1) The first energy value is less than or equal to a first preset energy threshold.
The first energy value D represents the initial energy value before wall filtering, and pixels with too high initial energy are generally considered to be strong reflection noise.
2) The ratio of the first energy value to the second energy value is less than or equal to a preset energy ratio threshold.
Pixels with too high a first energy value D/second energy value P represent noise below the cutoff frequency.
3) The second energy value is greater than or equal to a second preset energy threshold.
The second energy value P represents the power of the pixel point after wall filtering, and according to the characteristic of the high-pass filter, the power of the pass band is higher at this time, and random noise with a certain amplitude exists in the whole frequency band, so the point where the second energy value P is smaller than a certain threshold value should be noise.
4) The variance value is less than or equal to a preset variance threshold.
The variance value σ represents the variance of the signal, and the larger the variance value is, the stronger the randomness of the representative signal is, and the easier it is to be noise.
5) The gray scale value of the black and white image data is less than or equal to a preset gray scale value.
Black and white grey values are an indication of the reflected energy from tissue, and for strong reflection points with higher grey values it is unlikely to be liquid, i.e. not likely to be blood flow, and should therefore be filtered out as noise.
In this embodiment, the five determination conditions are used as the preset threshold condition, and when the data of a pixel meets all the conditions, it is indicated that the pixel is blood flow motion, and the data of the pixel is blood flow data. For the value of each preset value, those skilled in the art can perform corresponding preset according to different implementation scenarios and use environments, which is not limited by the present disclosure.
3. And (6) judging a threshold value.
In the present embodiment, in view of improving the imaging accuracy and reducing the computation amount of data analysis processing, the method of area division and pixel-by-pixel combination determination described in the embodiment of fig. 3 is specifically:
firstly, dividing a display area into a plurality of sub-areas of 6-by-6 pixels;
and calculating a first energy value D, a second energy value P, a variance value sigma and an average value of gray values of the pixel points in the sub-region, and taking the average value of the parameters as the image data of the sub-region.
And judging whether the data of each subarea simultaneously meets the preset threshold condition, if so, indicating that the subarea contains the color, and determining that the subarea is a color area. If not, the sub-area does not contain color, and the color value is 0.
And performing threshold judgment on each pixel point data in the color area one by one, judging whether the pixel point simultaneously meets the preset threshold condition, and if so, indicating that the pixel point data is blood flow data. If not, the pixel point is not blood flow motion.
And S10, estimating the speed.
The velocity estimation in the ultrasound imaging is based on different parameter displays provided by different display modes, such as a color display mode, an energy display mode, a variance velocity display mode, and the like, and any display mode can be calculated in step S09.
And S11, carrying out color coding imaging on the image according to the data of the pixel points representing the blood flow motion and displaying the image. For example, the direction of blood flow moving towards the transducer is indicated in red, the direction of blood flow moving away from the transducer is indicated in blue, while the speed of blood flow is indicated by the brightness of the color.
In some embodiments, to improve the continuity of the ultrasound color image, after the velocity estimation, the ultrasound image is subjected to a time-dependent process, i.e., a frame-dependent process, based on the blood flow velocity. The time-dependent processing includes:
acquiring blood flow velocity data of a previous frame image of a current frame image;
judging whether the difference value of the blood flow speed data of the current frame image and the blood flow speed data of the previous frame image is larger than a first preset threshold value or not;
and if not, performing time correlation processing on the current frame image and the previous frame image.
Specifically, the blood flow velocity data of the image of the previous frame of the current frame image may be obtained from the stored previous frame data, and the difference between the blood flow velocity data of the image of the previous frame and the blood flow velocity data of the image of the current frame may be compared with a first preset threshold value, where the first preset threshold value represents the maximum variation value of the two adjacent frames of images.
If the difference value of the two adjacent frames of image data is greater than the first preset threshold value, it indicates that the blood flow has shifted or the scanning position has been switched, and the two frames of images are unrelated without frame correlation processing. And if the difference value of the two adjacent frames of image data is less than or equal to a first preset threshold value, the two frames of images before and after are related, and the images are continuous images, and the current image is subjected to frame related processing.
When the current image is subjected to frame correlation processing, firstly, color data of N frames of images in front of the current image are obtained, and according to the blood flow speed of the N frames of images and the weight of each frame of image, the blood flow speed of the current image is obtained by processing pixel points one by one, and the time correlation processing is completed. And carrying out color coding and displaying according to the processed blood flow velocity data.
In the embodiment, in the time correlation processing, the change of the blood flow is fully considered, and whether the blood flow is shifted or not is judged according to the difference value of two adjacent frames of images, so that the time correlation processing is more accurate.
In a second aspect, as shown in fig. 5, the present disclosure provides an ultrasound color imaging processing apparatus, comprising:
the acquisition module 100 is configured to acquire color echo signals and obtain color feature data of each sampling location according to the color echo signals;
the space transformation module 200 is configured to perform DSC digital scanning transformation according to the color feature data to obtain color image data of each pixel point in the image display area;
a determining module 300, configured to determine that data of a pixel point meeting a preset threshold condition is blood flow data, where the preset threshold condition is determined according to color image data;
and the encoding module 400 is configured to obtain blood flow velocity data of corresponding pixel points according to the blood flow data, and perform color encoding imaging on the image according to the blood flow velocity data.
In some embodiments, the determining module 300 is specifically configured to:
dividing an image display area into a plurality of sub-areas, wherein each sub-area comprises a plurality of pixel points;
according to the color image data of the pixel points in the sub-area, calculating to obtain the color image data of the sub-area;
determining a sub-area meeting a preset threshold condition as a color area;
and determining the data of the pixel points meeting the preset threshold condition as blood flow data according to the color image data of each pixel point in the color area.
In some embodiments, the obtaining module 100 is further configured to obtain a black-and-white echo signal, and process the black-and-white echo signal to obtain black-and-white feature data of each sampling position;
the spatial transformation module 200 is further configured to perform DSC digital scanning transformation according to the black-and-white feature data to obtain black-and-white image data of each pixel point in the image display area, where the black-and-white image data includes a gray value;
the determining module 300 is further configured to divide the image display area into a plurality of sub-areas, where each sub-area includes a plurality of pixel points;
calculating black-and-white image data of the sub-region according to the black-and-white image data of the pixel points in the sub-region;
determining a sub-area meeting a preset threshold condition as a color area; the preset threshold condition is determined according to the color image data and the black and white image data;
and determining the data meeting the preset threshold condition as blood flow data according to the color image data and the black-and-white image data of each pixel point in the color area.
In some embodiments, the apparatus further comprises:
the wall filtering module is used for carrying out wall filtering processing on the color echo signals; the color image data comprises a first energy value before wall filtering of each pixel point, a second energy value after wall filtering of each pixel point and a variance value;
the autocorrelation module is used for carrying out autocorrelation processing on the color echo signals after wall filtering processing to obtain characteristic intermediate parameters of each sampling position; the color feature data includes feature intermediate parameters and a first energy value for each sampling location prior to wall filtering.
In a third aspect, the present disclosure provides an ultrasonic color imaging apparatus comprising:
a processor; and
a memory communicatively coupled to the processor and storing computer readable instructions executable by the processor, the processor performing a method according to any of the embodiments described above when the computer readable instructions are read.
In a fourth aspect, the present disclosure provides a storage medium storing computer-readable instructions for causing a computer to perform the method according to any of the above embodiments.
Specifically, fig. 6 shows a schematic structural diagram of a computer system 600 suitable for implementing the method or processor of the present disclosure, and the electronic device and the storage medium provided in the third and fourth aspects are implemented by the system shown in fig. 6.
As shown in fig. 6, the computer system 600 includes a Central Processing Unit (CPU) 601 that can perform various appropriate actions and processes according to a program stored in a Read Only Memory (ROM) 602 or a program loaded from a storage section 608 into a Random Access Memory (RAM) 603. In the RAM 603, various programs and data necessary for the operation of the system 600 are also stored. The CPU 601, ROM 602, and RAM 603 are connected to each other via a bus 604. An input/output (I/O) interface 605 is also connected to bus 604.
The following components are connected to the I/O interface 605: an input portion 606 including a keyboard, a mouse, and the like; an output portion 607 including a display such as a Cathode Ray Tube (CRT), a Liquid Crystal Display (LCD), and the like, and a speaker; a storage section 608 including a hard disk and the like; and a communication section 609 including a network interface card such as a LAN card, a modem, or the like. The communication section 609 performs communication processing via a network such as the internet. The driver 610 is also connected to the I/O interface 605 as needed. A removable medium 611 such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, or the like is mounted on the drive 610 as necessary, so that a computer program read out therefrom is mounted in the storage section 608 as necessary.
In particular, the above method processes may be implemented as a computer software program according to embodiments of the present disclosure. For example, embodiments of the present disclosure include a computer program product comprising a computer program tangibly embodied on a machine-readable medium, the computer program comprising program code for performing the above-described method. In such an embodiment, the computer program may be downloaded and installed from a network through the communication section 609 and/or installed from the removable medium 611.
The flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present disclosure. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems that perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
It should be understood that the above embodiments are only examples for clarity of description, and are not limiting. Other variations and modifications will be apparent to persons skilled in the art in light of the above description. And are neither required nor exhaustive of all embodiments. And obvious variations or modifications of the present disclosure may be made without departing from the scope of the present disclosure.
Claims (12)
1. An ultrasonic color imaging processing method, comprising:
acquiring color echo signals, and acquiring color characteristic data of each sampling position according to the color echo signals;
performing DSC digital scanning conversion according to the color characteristic data to obtain color image data of each pixel point of the image display area;
determining data of pixel points meeting a preset threshold condition as blood flow data, wherein the preset threshold condition is determined according to the color image data;
obtaining blood flow velocity data of corresponding pixel points according to the blood flow data, and carrying out color coding imaging on an image according to the blood flow velocity data; the determining that the data of the pixel points meeting the preset threshold condition is blood flow data includes:
dividing an image display area into a plurality of sub-areas, wherein each sub-area comprises a plurality of pixel points;
according to the color image data of the pixel points in the sub-area, calculating to obtain the color image data of the sub-area;
determining the sub-area meeting the preset threshold condition as a color area;
and determining the data of the pixel points meeting the preset threshold condition as blood flow data according to the color image data of each pixel point in the color area.
2. The method of claim 1, further comprising:
acquiring black and white echo signals, and processing the black and white echo signals to obtain black and white characteristic data of each sampling position;
performing DSC digital scanning conversion according to the black-and-white characteristic data to obtain black-and-white image data of each pixel point of an image display area, wherein the black-and-white image data comprises gray values;
dividing an image display area into a plurality of sub-areas, wherein each sub-area comprises a plurality of pixel points;
calculating black-and-white image data of the subarea according to the black-and-white image data of the pixel points in the subarea;
determining the sub-area meeting the preset threshold condition as a color area; the preset threshold condition is determined according to the color image data and the black and white image data;
and determining the data meeting the preset threshold condition as blood flow data according to the color image data and the black and white image data of each pixel point in the color area.
3. The method of claim 2, further comprising, prior to said deriving color characterization data for each sampling location from said color echo signals:
performing wall filtering processing on the color echo signals;
the color image data comprises a first energy value before wall filtering of each pixel point, a second energy value after wall filtering of each pixel point and a variance value.
4. The method of claim 3, wherein the preset threshold condition comprises:
the first energy value is less than or equal to a first preset energy threshold value;
the ratio of the first energy value to the second energy value is less than or equal to a preset energy ratio threshold;
the second energy value is greater than or equal to a second preset energy threshold;
the variance value is less than or equal to a preset variance threshold;
the gray value is less than or equal to a preset gray value;
the determining that the data of the pixel points meeting the preset threshold condition is blood flow data includes:
and when the data of the pixel point meets all conditions, determining that the data of the pixel point meets the preset threshold condition.
5. The method of claim 3, wherein said deriving color characterization data for each sampling location from the color echo signals comprises:
carrying out autocorrelation processing on the color echo signals subjected to wall filtering processing to obtain characteristic intermediate parameters of each sampling position;
obtaining the first energy value for each sample position before wall filtering; the color feature data includes the feature intermediate parameters and the first energy value for each sampling position before the wall filtering.
6. The method according to claim 5, wherein the performing the DSC digital scan conversion according to the color feature data to obtain color image data of each pixel point of an image display area comprises:
and calculating to obtain a second energy value and a variance value of each pixel point according to the characteristic intermediate parameters after DSC digital scanning conversion.
7. The method of claim 1, further comprising, after said color-coded imaging an image from said blood flow velocity data:
acquiring blood flow velocity data of a previous frame image of a current frame image;
judging whether the difference value of the blood flow speed data of the current frame image and the blood flow speed data of the previous frame image is larger than a first preset threshold value or not;
if not, time correlation processing is carried out on the current frame image.
8. An ultrasonic color imaging processing apparatus, comprising:
the acquisition module is used for acquiring color echo signals and acquiring color characteristic data of each sampling position according to the color echo signals;
the space transformation module is used for carrying out DSC digital scanning transformation according to the color characteristic data to obtain color image data of each pixel point of the image display area;
the determining module is used for determining that the data of the pixel points meeting a preset threshold condition is blood flow data, and the preset threshold condition is determined according to the color image data;
the encoding module is used for obtaining blood flow velocity data of corresponding pixel points according to the blood flow data and carrying out color encoding imaging on the image according to the blood flow velocity data; the determining module is specifically configured to:
dividing an image display area into a plurality of sub-areas, wherein each sub-area comprises a plurality of pixel points;
according to the color image data of the pixel points in the sub-area, calculating to obtain the color image data of the sub-area;
determining the sub-area meeting the preset threshold condition as a color area;
and determining the data of the pixel points meeting the preset threshold condition as blood flow data according to the color image data of each pixel point in the color area.
9. The apparatus of claim 8,
the acquisition module is further used for acquiring black and white echo signals and processing the black and white echo signals to obtain black and white characteristic data of each sampling position;
the space transformation module is further used for carrying out DSC digital scanning transformation according to the black-and-white characteristic data to obtain black-and-white image data of each pixel point of an image display area, and the black-and-white image data comprises a gray value;
the determining module is further configured to divide the image display area into a plurality of sub-areas, where each sub-area includes a plurality of pixel points;
calculating black-and-white image data of the sub-region according to the black-and-white image data of the pixel points in the sub-region;
determining the sub-area meeting the preset threshold condition as a color area; the preset threshold condition is determined according to the color image data and the black and white image data;
and determining the data meeting the preset threshold condition as blood flow data according to the color image data and the black-and-white image data of each pixel point in the color area.
10. The apparatus of claim 8, further comprising:
the wall filtering module is used for carrying out wall filtering processing on the color echo signals; the color image data comprises a first energy value before wall filtering of each pixel point, a second energy value after wall filtering of each pixel point and a variance value;
the autocorrelation module is used for carrying out autocorrelation processing on the color echo signals after wall filtering processing to obtain characteristic intermediate parameters of each sampling position; the color feature data includes the feature intermediate parameters and the first energy value for each sampling position before the wall filtering.
11. An ultrasonic color imaging apparatus, comprising:
a processor; and
a memory, communicatively coupled to the processor, storing computer readable instructions executable by the processor, the processor performing the method of any of claims 1 to 7 when the computer readable instructions are read.
12. A storage medium having stored thereon computer-readable instructions for causing a computer to perform the method of any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010097056.9A CN111265246B (en) | 2020-02-17 | 2020-02-17 | Ultrasonic color imaging processing method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010097056.9A CN111265246B (en) | 2020-02-17 | 2020-02-17 | Ultrasonic color imaging processing method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111265246A CN111265246A (en) | 2020-06-12 |
| CN111265246B true CN111265246B (en) | 2023-02-24 |
Family
ID=70991142
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010097056.9A Active CN111265246B (en) | 2020-02-17 | 2020-02-17 | Ultrasonic color imaging processing method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111265246B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2021253230A1 (en) * | 2020-06-16 | 2021-12-23 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic imaging method, ultrasonic imaging system and computer storage medium |
| CN112826535B (en) * | 2020-12-31 | 2022-09-09 | 青岛海信医疗设备股份有限公司 | Method, device and equipment for automatically positioning blood vessel in ultrasonic imaging |
| CN113096056B (en) * | 2021-04-06 | 2022-04-12 | 全景恒升（北京）科学技术有限公司 | Intravascular image fusion method based on region complementation |
| CN114170239A (en) * | 2021-11-11 | 2022-03-11 | 苏州国科昂卓医疗科技有限公司 | Large-view-field self-adaptive wall filtering method and system based on hierarchical search |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0690952A (en) * | 1992-06-16 | 1994-04-05 | Toshiba Corp | Ultrasonic diagnostic device |
| CN101416885A (en) * | 2007-10-22 | 2009-04-29 | 深圳迈瑞生物医疗电子股份有限公司 | Frame equilibration circuit for ultrasonic imaging system |
| CN102805649A (en) * | 2011-06-03 | 2012-12-05 | 深圳迈瑞生物医疗电子股份有限公司 | Color ultrasonic imaging method and device |
| CN105590315A (en) * | 2015-12-04 | 2016-05-18 | 南京云石医疗科技有限公司 | Ultrasonic color Doppler image post-processing method |
-
2020
- 2020-02-17 CN CN202010097056.9A patent/CN111265246B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0690952A (en) * | 1992-06-16 | 1994-04-05 | Toshiba Corp | Ultrasonic diagnostic device |
| CN101416885A (en) * | 2007-10-22 | 2009-04-29 | 深圳迈瑞生物医疗电子股份有限公司 | Frame equilibration circuit for ultrasonic imaging system |
| CN102805649A (en) * | 2011-06-03 | 2012-12-05 | 深圳迈瑞生物医疗电子股份有限公司 | Color ultrasonic imaging method and device |
| CN105590315A (en) * | 2015-12-04 | 2016-05-18 | 南京云石医疗科技有限公司 | Ultrasonic color Doppler image post-processing method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111265246A (en) | 2020-06-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111265246B (en) | Ultrasonic color imaging processing method and device | |
| JP4810048B2 (en) | Method and apparatus for flash suppression in two-dimensional ultrasound imaging | |
| US6277075B1 (en) | Method and apparatus for visualization of motion in ultrasound flow imaging using continuous data acquisition | |
| JP4757307B2 (en) | Ultrasonic image processing device | |
| US20030174890A1 (en) | Image processing device and ultrasonic diagnostic device | |
| US8202220B2 (en) | Ultrasound diagnostic apparatus | |
| US20060184021A1 (en) | Method of improving the quality of a three-dimensional ultrasound doppler image | |
| JP4627366B2 (en) | Method and apparatus for motion visualization in ultrasonic flow imaging using packet data acquisition | |
| US20080077011A1 (en) | Ultrasonic apparatus | |
| US9039620B2 (en) | Ultrasound diagnostic apparatus | |
| EP1050761B1 (en) | Color imaging display method and apparatus for color doppler imaging system | |
| CN105073020A (en) | Ultrasonic diagnostic device and ultrasonic image processing method | |
| US8500646B2 (en) | Color Doppler mode image processing in an ultrasound system | |
| EP1796037B1 (en) | Image processing system and method for controlling gains for color flow images | |
| JP2020524552A (en) | Method and system for compound ultrasound imaging | |
| JP4347860B2 (en) | Ultrasonic diagnostic equipment | |
| JP3802462B2 (en) | Ultrasonic diagnostic equipment | |
| US6500125B1 (en) | Ultrasound b/color priority threshold calculation | |
| JP3662835B2 (en) | Ultrasonic diagnostic equipment | |
| JP2006051202A (en) | Image processing device, image processing method, and program | |
| JP2723453B2 (en) | Ultrasound blood flow observation device | |
| JP4616614B2 (en) | Ultrasonic diagnostic equipment | |
| JP5134757B2 (en) | Image processing apparatus, image processing method, and ultrasonic diagnostic apparatus | |
| CN117522722A (en) | Ultrasonic Doppler blood flow image processing method and related equipment | |
| JP2005318921A (en) | Ultrasonic diagnostic equipment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |