Blood flow parameter determination and display method, apparatus, and electronic device and storage medium

Abstract

The application discloses a blood flow parameter determination and display method and device, an electronic device and a computer readable storage medium. The blood flow parameter determination method comprises the following steps: acquiring echo data of multiple pulses of wall-filtered pixel points in an imaging region; constructing four sub-sequences based on the echo data, and calculating fourth-order cumulants of the echo data based on the four sub-sequences; and calculating the blood flow parameter based on the echo data and the fourth-order cumulants. In the application, the four sub-sequences are constructed based on the echo data of multiple pulses of wall-filtered pixel points in the imaging region, the fourth-order cumulants of the echo data are calculated by using the constructed sub-sequences, and the problem that the autocorrelation method is affected by the accuracy of the blood flow parameter estimation in the additive Gaussian color noise and multiplicative noise background is eliminated by using the characteristics that the fourth-order cumulants are not sensitive to the additive Gaussian color noise and multiplicative noise, so that the accuracy of the blood flow parameter estimation is improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of ultrasound, in particular to a blood flow parameter determination and display method and device, an electronic device and a computer readable storage medium. BACKGROUND

[0002] Currently, a self-correlation method is used in a color blood flow imaging system to estimate blood flow parameters. The method assumes that only additive noise exists in echo signals without multiplicative noise. However, multiplicative noise exists in actual complex signals, which affects the accuracy of blood flow parameter estimation. Meanwhile, additive Gaussian white noise becomes additive Gaussian colored noise after wall filtering. The autocorrelation function of Gaussian colored noise is no longer an impulse function, which further affects the accuracy of the estimation result.

[0003] Therefore, how to improve the accuracy of blood flow parameter estimation is a technical problem to be solved by those skilled in the art. SUMMARY

[0004] The present application aims to provide a blood flow parameter determination and display method and device, an electronic device and a computer readable storage medium, which improves the accuracy of blood flow parameter estimation.

[0005] To achieve the above-mentioned purpose, the present application provides a blood flow parameter determination method, comprising:

[0006] Obtaining echo data of multiple pulses of wall-filtered pixel points in an imaging region;

[0007] Constructing four sub-sequences based on the echo data, and calculating fourth-order cumulants of the echo data based on the four sub-sequences;

[0008] Calculating blood flow parameters based on the echo data and the fourth-order cumulants.

[0009] Among them, the four sub-sequences are constructed based on the echo data, comprising:

[0010] Delay processing the echo data according to a preset delay rule to construct four sub-sequences.

[0011] Among them, the delay processing the echo data according to a preset delay rule to construct four sub-sequences, comprising:

[0012] Determining a first delay parameter, a second delay parameter and a third delay parameter; wherein the difference between the sum of the second delay parameter and the third delay parameter and the first delay parameter is 1.

[0013] constructing a first sub-sequence based on the echo data; wherein a first element of the first sub-sequence is a first element of the echo data, and a length of the first sub-sequence is a difference between a length of the echo data and a target maximum value, the target maximum value being a maximum value of the first delay parameter and the second delay parameter;

[0014] delay processing the echo data based on the first delay parameter to construct a second sub-sequence; wherein a length of the second sub-sequence is a difference between the length of the echo data and a first target difference value, the first target difference value being a difference between the target maximum value and the first delay parameter;

[0015] delay processing the echo data based on the second delay parameter to construct a third sub-sequence; wherein a length of the third sub-sequence is a difference between the length of the echo data and a second target difference value, the second target difference value being a difference between the target maximum value and the second delay parameter;

[0016] delay processing the echo data based on the third delay parameter to construct a fourth sub-sequence; wherein a length of the fourth sub-sequence is a difference between the length of the echo data and a third target difference value, the third target difference value being a difference between the target maximum value and the third delay parameter.

[0017] The four sub-sequences include the first sub-sequence, the second sub-sequence, the third sub-sequence, and the fourth sub-sequence.

[0018] The calculating the fourth-order cumulant of the echo data based on the four sub-sequences includes:

[0019] calculating a first point multiplication result between a conjugate of the first sub-sequence, a conjugate of the second sub-sequence, the third sub-sequence, and the fourth sub-sequence, and accumulating each element in the first point multiplication result to obtain a first accumulated value;

[0020] calculating a second point multiplication result between the conjugate of the first sub-sequence and the third sub-sequence, and accumulating each element in the second point multiplication result to obtain a second accumulated value;

[0021] calculating a third point multiplication result between the conjugate of the second sub-sequence and the fourth sub-sequence, and accumulating each element in the third point multiplication result to obtain a third accumulated value;

[0022] calculating a fourth point multiplication result between the conjugate of the first sub-sequence and the fourth sub-sequence, and accumulating each element in the fourth point multiplication result to obtain a fourth accumulated value;

[0023] conjugate of the second sub-sequence, a fifth point multiplication result between the third sub-sequence, and accumulating an absolute value of each element in the fifth point multiplication result to obtain a fifth accumulated value;

[0024] calculating a first product between the second accumulated value and the third accumulated value, and a second product between the fourth accumulated value and the fifth accumulated value;

[0025] taking a difference between the first accumulated value and the first product, the second product as the fourth-order cumulant of the echo data.

[0026] wherein the blood flow parameter comprises any one or a combination of blood flow energy, blood flow velocity and blood flow variance;

[0027] calculating the blood flow parameter based on the echo data and the fourth-order cumulant comprises:

[0028] taking a product between the echo data and a conjugate of the echo data as an intermediate sequence, and accumulating each element in the intermediate sequence to obtain the blood flow energy;

[0029] and / or, performing an angle operation on a negative of the fourth-order cumulant to obtain a pulse frequency, and calculating the blood flow velocity according to a corresponding relationship between the pulse frequency and the blood flow velocity;

[0030] and / or, calculating a fourth-order absolute value cumulant based on the four sub-sequences, calculating a ratio between an absolute value of the fourth-order cumulant and the fourth-order absolute value cumulant, and taking a difference between 1 and the ratio as the blood flow variance.

[0031] wherein the calculating the fourth-order absolute value cumulant based on the four sub-sequences comprises:

[0032] calculating a first point multiplication result between a conjugate of the first sub-sequence, a conjugate of the second sub-sequence, the third sub-sequence and the fourth sub-sequence, and accumulating an absolute value of each element in the first point multiplication result to obtain a sixth accumulated value;

[0033] calculating a second point multiplication result between a conjugate of the first sub-sequence and the third sub-sequence, and accumulating an absolute value of each element in the second point multiplication result to obtain a seventh accumulated value;

[0034] calculating a third point multiplication result between a conjugate of the second sub-sequence and the fourth sub-sequence, and accumulating an absolute value of each element in the third point multiplication result to obtain an eighth accumulated value;

[0035] calculating a fourth point multiplication result between a conjugate of the first sub-sequence and the fourth sub-sequence, and accumulating an absolute value of each element in the fourth point multiplication result to obtain a ninth accumulated value;

[0036] calculating a fifth point multiplication result between the second subsequence and the third subsequence, accumulating the absolute value of each element in the fifth point multiplication result to obtain a tenth accumulated value;

[0037] calculating a third product between the seventh accumulated value and the eighth accumulated value, and a fourth product between the ninth accumulated value and the tenth accumulated value;

[0038] accumulating the sum of the sixth accumulated value, the third product and the fourth product as a fourth-order absolute value accumulation.

[0039] To achieve the above object, the present application provides a blood flow parameter display method, comprising:

[0040] acquiring echo data of multiple pulses of a pixel point after wall filtering in an imaging region;

[0041] constructing four sub-sequences based on the echo data, and calculating a fourth-order accumulation of the echo data based on the four sub-sequences;

[0042] calculating a blood flow parameter based on the echo data and the fourth-order accumulation;

[0043] displaying the blood flow parameter.

[0044] To achieve the above object, the present application provides a blood flow parameter determination device, comprising:

[0045] an acquisition module, configured to acquire echo data of multiple pulses of a pixel point after wall filtering in an imaging region;

[0046] a first calculation module, configured to construct four sub-sequences based on the echo data, and calculate a fourth-order accumulation of the echo data based on the four sub-sequences;

[0047] a second calculation module, configured to calculate a blood flow parameter based on the echo data and the fourth-order accumulation.

[0048] To achieve the above object, the present application provides a blood flow parameter display device, comprising:

[0049] an acquisition module, configured to acquire echo data of multiple pulses of a pixel point after wall filtering in an imaging region;

[0050] a first calculation module, configured to construct four sub-sequences based on the echo data, and calculate a fourth-order accumulation of the echo data based on the four sub-sequences;

[0051] a second calculation module, configured to calculate a blood flow parameter based on the echo data and the fourth-order accumulation.

[0052] a display module configured to display the blood flow parameter.

[0053] To achieve the above object, the present application provides an electronic device, comprising:

[0054] a memory configured to store a computer program;

[0055] a processor configured to execute the computer program to implement the steps of the blood flow parameter determination method.

[0056] To achieve the above object, the present application provides an electronic device, comprising:

[0057] a memory configured to store a computer program;

[0058] a processor configured to execute the computer program to implement the steps of the blood flow parameter display method.

[0059] To achieve the above object, the present application provides a computer readable storage medium, wherein the computer readable storage medium stores a computer program, and the computer program is executed by a processor to implement the steps of the blood flow parameter determination method.

[0060] To achieve the above object, the present application provides a computer readable storage medium, wherein the computer readable storage medium stores a computer program, and the computer program is executed by a processor to implement the steps of the blood flow parameter display method.

[0061] According to the above scheme, the blood flow parameter determination method provided by the present application comprises: acquiring echo data of multiple pulses of a pixel point after wall filtering in an imaging region; constructing four sub-sequences based on the echo data, and calculating fourth-order cumulants of the echo data based on the four sub-sequences; and calculating a blood flow parameter based on the echo data and the fourth-order cumulants.

[0062] The blood flow parameter determination method provided by the present application constructs four sub-sequences based on echo data of multiple pulses after wall filtering in an imaging region, calculates fourth-order cumulants of the echo data by using the constructed sub-sequences, and eliminates the problem of affecting the accuracy of blood flow parameter estimation by using the characteristics of the fourth-order cumulants that are not sensitive to additive Gaussian color noise and multiplicative noise, thereby improving the accuracy of blood flow parameter estimation.

[0063] The present application also discloses a blood flow parameter display method, a blood flow parameter determination device, a blood flow parameter display device, an electronic device and a computer readable storage medium, which can also achieve the above technical effects.

[0064] It should be understood that the above general description and the following detailed description are only exemplary and cannot limit the present application. Attached Figure Description

[0065] To more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort. The drawings are used to provide a further understanding of this disclosure and constitute a part of the specification. They are used together with the following detailed description to explain this disclosure, but do not constitute a limitation of this disclosure. In the drawings:

[0066] Figure 1 This is a flowchart illustrating a method for determining blood flow parameters according to an exemplary embodiment;

[0067] Figure 2 This is a schematic diagram illustrating a representation format of echo data according to an exemplary embodiment;

[0068] Figure 3 This is a schematic diagram illustrating the relationship between Doppler frequency and blood flow according to an exemplary embodiment;

[0069] Figure 4 This is a flowchart illustrating a method for displaying blood flow parameters according to an exemplary embodiment;

[0070] Figure 5 This is a structural diagram of a blood flow parameter determination device according to an exemplary embodiment;

[0071] Figure 6 This is a structural diagram of a blood flow parameter display device according to an exemplary embodiment;

[0072] Figure 7 This is a structural diagram of an electronic device according to an exemplary embodiment. Detailed Implementation

[0073] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of this application, and not all of them. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application. Furthermore, in the embodiments of this application, "first," "second," etc., are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence.

[0074] The method for determining blood flow parameters provided in this application embodiment can be applied to the following application scenarios:

[0075] The application scenario includes an ultrasound device and an electronic device, both of which can be connected to a network. The ultrasound device can be various types of devices, such as a B-ultrasound, a color Doppler ultrasound, an endoscopic Doppler ultrasound, an intravascular ultrasound, etc. The electronic device is a device capable of processing ultrasound data, which can be a terminal device or a server. When it is a terminal device, it can be a personal computer, a tablet computer, a smartphone, etc. In a specific application scenario, the ultrasound device is configured with an ultrasound probe. The ultrasound device transmits an ultrasound wave pulse to a detected object (such as a blood vessel region) at a certain frequency through the ultrasound probe. The ultrasound wave pulse generates a return signal through reflection of the detected object and is received by the ultrasound probe. The ultrasound device processes the return signal received by the ultrasound probe to obtain return data of multiple pulses of the imaged region wall-filtered pixel points, and sends the return data to the electronic device. The electronic device processes the obtained return data to calculate a blood flow parameter.

[0076] In some application scenarios, the electronic device can return the calculated blood flow parameter to the ultrasound device for display.

[0077] Embodiments of the present application disclose a blood flow parameter determination method, which improves the accuracy of blood flow parameter estimation. The execution subject of the method can be an electronic device.

[0078] Referring to Figure 1 , a flowchart of a blood flow parameter determination method according to an exemplary embodiment is shown, as Figure 1 shown, comprising:

[0079] S101: Obtain return data of multiple pulses of imaged region wall-filtered pixel points;

[0080] The ultrasound device CFM (color Doppler) mode utilizes the characteristic that blood flow motion causes multiple ultrasound returns to produce Doppler frequency shift, and indirectly describes blood flow motion by estimating Doppler frequency. In a specific implementation, the imaged region wall-filtered complex data x L×P×En is acquired in the CFM mode of the ultrasound device, as Figure 2 shown, L is Line (line), P is Point (point), En is Ensemble (snap shot), (L, P) represents a pixel point in the complex data, En is the number of pulse repetitions, and return data of multiple pulses of each pixel point (L, P) is obtained The imaged region wall-filtered complex data x L×P×En is a three-dimensional complex matrix, that is, it is data composed of three dimensions of L, P and En. The data in the En dimension can be referred to as En-direction data, which is used to estimate the blood flow Doppler frequency and further obtain the blood flow velocity.

[0081] In the embodiment, an additive and multiplicative noise echo model Y(t) = T(t)x(t) + n(t) is assumed, Y(t) represents the acquired echo data, x(t) represents the real echo data, n(t) represents the additive noise, and T(t) represents the multiplicative noise. That is, additive noise and multiplicative noise can exist in the acquired echo data.

[0082] Optionally, the ultrasound device acquires real data of the imaging region through the ultrasound probe, obtains complex data after demodulating the real data, and then performs wall filtering processing to filter out data with small motion trend such as blood vessel walls, to obtain wall-filtered complex data of the imaging region.

[0083] Optionally, the imaging region can correspond to a plurality of pixel points. In the specific implementation process, a single or multiple pixel points can be polled, and the corresponding echo data is acquired and the blood flow parameter is calculated until all the pixel points in the imaging region complete the blood flow parameter calculation. In addition, the echo data can also be acquired in units of pixels, and then the blood flow parameter is calculated as a whole object of the imaging region.

[0084] It should be noted that the echo data is data obtained based on multiple pulses, and the echo signal of each pulse can correspond to a data element. Arranging the data elements corresponding to different pulses can obtain a data sequence. Therefore, the echo data corresponding to all pulses can be a sequence composed of multiple data elements.

[0085] S102: four sub-sequences are constructed based on the echo data, and fourth-order cumulants of the echo data are calculated based on the four sub-sequences;

[0086] In this step, first, four sub-sequences are constructed based on the echo data The four sub-sequences are constructed according to a certain delay rule, that is, the echo data is delayed according to a predetermined delay rule to construct four sub-sequences (first sub-sequence, second sub-sequence, third sub-sequence, and fourth sub-sequence). The length of the sub-sequence can be equal to or shorter than the length of the echo data. For the case of being shorter, it can be shorter than the length of the echo data by 1 element, 2 elements, etc. In addition, the sub-sequence can be a continuous segment of elements in the echo data, or the elements can be selected at intervals and the elements at intervals can be determined as the sub-sequence.

[0087] Optionally, a sequence of a certain length is selected from the echo data to form the first sub-sequence, the second sub-sequence, the third sub-sequence, and the fourth sub-sequence. The lengths of the four sequences are the same and less than the length of the echo data. For example, for echo data containing N elements, a sequence formed by 1 to N-3 elements is determined as the first sub-sequence, a sequence formed by 2 to N-2 elements is determined as the second sub-sequence, a sequence formed by 3 to N-1 elements is determined as the third sub-sequence, and a sequence formed by 4 to N elements is determined as the fourth sub-sequence.

[0088] As a feasible implementation, the echo data is processed according to the preset delay rule to construct four sub-sequences, including: determining a first delay parameter, a second delay parameter, and a third delay parameter; wherein the difference between the sum of the second delay parameter and the third delay parameter and the first delay parameter is 1; constructing a first sub-sequence based on the echo data; wherein the first element of the first sub-sequence is the first element of the echo data, and the length of the first sub-sequence is the difference between the length of the echo data and a target maximum value, the target maximum value being the maximum value of the second delay parameter and the third delay parameter; processing the echo data according to the first delay parameter to construct a second sub-sequence; wherein the length of the second sub-sequence is the difference between the length of the echo data and a first target difference value, the first target difference value being the difference between the target maximum value and the first delay parameter; processing the echo data according to the second delay parameter to construct a third sub-sequence; wherein the length of the third sub-sequence is the difference between the length of the echo data and a second target difference value, the second target difference value being the difference between the target maximum value and the second delay parameter; processing the echo data according to the third delay parameter to construct a fourth sub-sequence; wherein the length of the fourth sub-sequence is the difference between the length of the echo data and a third target difference value, the third target difference value being the difference between the target maximum value and the third delay parameter.

[0089] The specific implementation is as follows:

[0090]

[0091] wherein “:” represents that the sequence index is incremented by 1 (it should be noted that in some embodiments, the sequence index can also be incremented by other numbers, for example, incremented by 2, incremented by 3, etc.), τ, γ are the first delay parameter, the second delay parameter, and the third delay parameter, respectively, is the target maximum value, is the first target difference value, is a second target difference value, is a third target difference value, are respectively a first sub-sequence, a second sub-sequence, a third sub-sequence, and a fourth sub-sequence.

[0092] Embodiments of the present application need to determine a fourth-order cumulant, i.e. a point multiplication result of four sequences, to determine blood flow velocity and other blood flow parameters. Assuming that an estimated value of the fourth-order cumulant obtained is then a coefficient of a Doppler frequency corresponding to the blood flow velocity is Let It can be ensured that the estimated value of the blood flow velocity will not be amplified or reduced. Therefore, preferably, a sum of the second delay parameter and the third delay parameter is different from the first delay parameter, i.e. At this time, the performance is better, i.e. the accuracy of the calculated blood flow parameter is higher. Specifically, the blood flow velocity of the kidney is slow, and the velocity value range is wide; the blood flow velocity of the carotid artery is large, and the velocity value range is relatively concentrated; considering the characteristics of the kidney blood flow and the carotid artery blood flow, they can be respectively set as or This setting can make the sequence the longest and retain as much data as possible, so that the estimated value obtained is optimal.

[0093] Further, according to the four sub-sequences constructed the fourth-order cumulant of the echo data is calculated. First, the conjugate of two delay sub-sequences is determined, the determined conjugate of the delay sub-sequences is multiplied with other delay sub-sequences, and the fourth-order cumulant is obtained based on the multiplication result. The process of obtaining the fourth-order cumulant based on the multiplication result can be: a first accumulated value is obtained by accumulating the multiplication result, then the conjugate of each delay sub-sequence is multiplied with each delay sub-sequence which has not been conjugated, four accumulated values, i.e. a second accumulated value, a third accumulated value, a fourth accumulated value, and a fifth accumulated value, are obtained by accumulating the multiplication result respectively, then the accumulated values corresponding to the conjugate of different delay sub-sequences are multiplied, two multiplication results are obtained, and finally the difference between the first accumulated value and the two multiplication results is determined as the fourth-order cumulant.

[0094] As a feasible implementation, the calculating the fourth-order cumulant of the echo data based on the four sub-sequences comprises: calculating a first point multiplication result between a conjugate of the first sub-sequence, a conjugate of the second sub-sequence, the third sub-sequence and the fourth sub-sequence, accumulating each element in the first point multiplication result to obtain a first accumulated value; calculating a second point multiplication result between the conjugate of the first sub-sequence and the third sub-sequence, accumulating each element in the second point multiplication result to obtain a second accumulated value; calculating a third point multiplication result between the conjugate of the second sub-sequence and the fourth sub-sequence, accumulating each element in the third point multiplication result to obtain a third accumulated value; calculating a fourth point multiplication result between the conjugate of the first sub-sequence and the fourth sub-sequence, accumulating each element in the fourth point multiplication result to obtain a fourth accumulated value; calculating a fifth point multiplication result between the conjugate of the second sub-sequence and the third sub-sequence, accumulating each element in the fifth point multiplication result to obtain a fifth accumulated value; calculating a first product between the second accumulated value and the third accumulated value, and a second product between the fourth accumulated value and the fifth accumulated value; taking a difference between the first accumulated value and the first product and the second product as the fourth-order cumulant of the echo data. The specific implementation is as follows:

[0095]

[0096] wherein, ⊙ represents point multiplication, that is, multiplication of elements corresponding to positions of two sequences, is a conjugate of the first sub-sequence, is a conjugate of the second sub-sequence, is a first point multiplication result, is a first accumulated value, is a second point multiplication result, is a second accumulated value, is a third point multiplication result, is a third accumulated value, is a fourth point multiplication result, is a fourth accumulated value, is a fifth point multiplication result, is a fifth accumulated value, is a first product, is a second product.

[0097] S103: calculating a blood flow parameter based on the echo data and the fourth-order cumulant;

[0098] In this step, the blood flow parameter is estimated by using the echo data of multiple pulses of each pixel point (L, P) and the fourth-order cumulant The blood flow parameter can include blood flow energy, blood flow velocity and blood flow variance, etc.

[0099] Optionally, the echo data or the cross-correlation function can be calculated alone to obtain the corresponding blood flow parameters, for example, the blood flow energy can be calculated based on the echo data, the blood flow velocity can be calculated based on the cross-correlation function, and the blood flow variance can be calculated based on the cross-correlation function. The echo data and the cross-correlation function can also be operated together, for example: addition, multiplication, division, etc. The corresponding blood flow parameters are obtained based on the operation results.

[0100] The calculation method of the blood flow energy Power is: taking the product between the echo data and the conjugate of the echo data as an intermediate sequence, and accumulating each element in the intermediate sequence to obtain the blood flow energy. The specific implementation method is as follows:

[0101]

[0102] The calculation method of the blood flow velocity Velocity (i.e. v) is: performing angle operation on the negative number of the fourth-order cumulant to obtain the pulse frequency, and calculating the blood flow velocity according to the corresponding relationship between the pulse frequency and the blood flow velocity. The specific implementation method is It can be understood that, as Figure 3 shown, the Doppler frequency of the ultrasonic equipment, that is, the echo frequency f d is related to the blood flow velocity v as: f d = 2vcosθ / λ, where θ represents the angle between the ultrasonic wave and the blood flow velocity direction, and λ represents the ultrasonic carrier wavelength, and v = f d / (2cosθ / λ) can be obtained. In the above formula, f d is obtained after angle operation, and since the blood flow direction θ is unknown, the value f d is generally used to represent the blood flow velocity Velocity, that is, the blood flow velocity value without blood flow direction estimation, which can be understood as qualitative blood flow.

[0103] The calculation method of the blood flow variance Variance is: calculating the fourth-order absolute value cumulant based on the four sub-sequences, calculating the ratio of the absolute value of the fourth-order cumulant to the fourth-order absolute value cumulant, and taking the difference between 1 and the ratio as the blood flow variance. The specific implementation method is as follows:

[0104]

[0105] wherein, is the fourth-order absolute value cumulant.

[0106] In the calculation of the fourth-order absolute value cumulant, the conjugate of two of the delay sub-sequences is first determined, the determined conjugate of the delay sub-sequences is point-multiplied with the other delay sub-sequences, the absolute values of the point multiplication results are accumulated to obtain a sixth accumulated value, the conjugate of each delay sub-sequence is respectively point-multiplied with each delay sub-sequence which is not conjugated, the absolute values of the point multiplication results are respectively accumulated to obtain four accumulated values, i.e., a seventh accumulated value, an eighth accumulated value, a ninth accumulated value and a tenth accumulated value, the accumulated values corresponding to the conjugate of different delay sub-sequences are point-multiplied to obtain two point multiplication results, and finally the sum of the sixth accumulated value and the two point multiplication results is determined as the fourth-order absolute value cumulant.

[0107] As a feasible implementation, the fourth-order absolute value cumulant is calculated based on the four sub-sequences, which includes: calculating the conjugate of the first sub-sequence, the conjugate of the second sub-sequence, the third sub-sequence and the fourth sub-sequence, accumulating the absolute values of each element in the first point multiplication result to obtain a sixth accumulated value; calculating the conjugate of the first sub-sequence and the third sub-sequence, accumulating the absolute values of each element in the second point multiplication result to obtain a seventh accumulated value; calculating the conjugate of the second sub-sequence and the fourth sub-sequence, accumulating the absolute values of each element in the third point multiplication result to obtain an eighth accumulated value; calculating the conjugate of the first sub-sequence and the fourth sub-sequence, accumulating the absolute values of each element in the fourth point multiplication result to obtain a ninth accumulated value; calculating the conjugate of the second sub-sequence and the third sub-sequence, accumulating the absolute values of each element in the fifth point multiplication result to obtain a tenth accumulated value; calculating the third product between the seventh accumulated value and the eighth accumulated value, and the fourth product between the ninth accumulated value and the tenth accumulated value; and taking the sum of the sixth accumulated value and the third product and the fourth product as the fourth-order absolute value cumulant. The specific implementation is as follows:

[0108]

[0109] wherein, is the sixth accumulated value,

[0110] is the seventh accumulated value,

[0111] is the eighth accumulated value,

[0112] is the ninth accumulated value,

[0113] is the tenth accumulated value,

[0114] The third product is

[0115] The fourth product is

[0116] The following explains why high-order cumulants are insensitive to Gaussian colored noise:

[0117] Suppose the distribution of random variable x is represented as x ~ N(0, σ 2 )S, so the probability density function of x is:

[0118] Therefore, the moment generating function of the Gaussian function is:

[0119]

[0120] In the integral formula Joint two formulas, let The moment generating function is:

[0121]

[0122] The derivatives of Φ(ω) are:

[0123]

[0124]

[0125]

[0126]

[0127] From the relationship between high-order moments and moment generating functions, we can get

[0128] m1 = 0, m2 = σ 2 , m3 = 0, m4 = 3σ 4

[0129] By extension, for any integer k, the sum of moments of a Gaussian random variable is written as:

[0130]

[0131] From the moment generating function, we can directly get the cumulant generating function of the Gaussian random variable x:

[0132] The derivatives are:

[0133] ψ 1 (ω) = -σ 2 ω

[0134] ψ 2(ω) = -σ 2

[0135] ψ k (ω)≡0,k=3,4,…

[0136] From the relationship between high-order cumulants and cumulant generating function, we have:

[0137] c1=0,c2=σ 2 ,m3≡0,k=3,4,…

[0138] The above results of high-order moments and high-order cumulants can be extended to: the second-order moments and the second-order cumulants of any zero-mean Gaussian random process are the same, both equal to its variance σ 2 ; the odd-order moments are always zero, but the even-order moments are not equal to zero; and the high-order cumulants (third-order and above) are always equal to zero. Therefore, in this sense, the high-order cumulants are "blind" to the Gaussian random process, that is, the high-order cumulants are not sensitive to the Gaussian color noise.

[0139] The method for determining blood flow parameters provided by the embodiments of the present application constructs four sub-sequences based on the echo data of multiple pulses of the wall-filtered pixel points in the imaging region, calculates the fourth-order cumulants of the echo data by using the constructed sub-sequences, and eliminates the problem of affecting the accuracy of blood flow parameter estimation by using the characteristics that the fourth-order cumulants are not sensitive to additive Gaussian color noise and multiplicative noise, thereby improving the accuracy of blood flow parameter estimation. Further, since there is more Gaussian color noise at the edges of blood vessels and in low-speed blood flow regions, eliminating the influence of Gaussian color noise can effectively ensure the stability of estimation at the edges of blood vessels and in low-speed blood flow regions, reduce the possibility of sudden state occurrence, and improve the continuity and hierarchy of blood flow imaging. In addition, since the embodiments of the present application eliminate the influence of multiplicative noise, the robustness of estimating blood flow parameters in the case of small samples is higher.

[0140] The embodiments of the present application disclose a blood flow parameter display method, specifically, referring to Figure 4 , a flow chart of a blood flow parameter display method according to an example embodiment is shown, as Figure 4 shown, comprising:

[0141] S201: acquiring echo data of multiple pulses of wall-filtered pixel points in an imaging region;

[0142] S202: constructing four sub-sequences based on the echo data, and calculating fourth-order cumulants of the echo data based on the four sub-sequences;

[0143] S203: calculating blood flow parameters based on the echo data and the fourth-order cumulants;

[0144] S204: displaying the blood flow parameters.

[0145] In this step, the blood flow parameter of each pixel point is estimated and displayed. The blood flow parameter can be directly displayed in the display screen in the form of text, or the blood flow image can be processed based on the blood flow parameter and the processed blood flow image can be displayed to indirectly display the blood flow parameter.

[0146] It should be noted that the execution subject of the blood flow parameter display method can be an ultrasound device. The ultrasound device acquires echo data, calculates the blood flow parameter based on the echo data, and displays the blood flow parameter based on the calculated blood flow parameter. In addition, the execution subject of the blood flow parameter display method can also be an electronic device. The electronic device acquires echo data from the ultrasound device, calculates the blood flow parameter based on the echo data, and displays the blood flow parameter based on the calculated blood flow parameter. Of course, the electronic device can also send the calculated blood flow parameter to the ultrasound device for display on the display screen of the ultrasound device.

[0147] Optionally, based on the blood flow parameter, the color, contrast, etc. of the blood vessel in the ultrasound image are adjusted, or the ultrasound image is rendered, and then displayed in the display screen. For example, for a certain pixel point in the ultrasound image, the blood flow velocity between (v0, v1) is displayed in orange, the blood flow velocity between (v1, v2) is displayed in magenta, and the blood flow velocity between (v2, v3) is displayed in deep red. According to the blood flow velocity of each pixel point, the pixel point is rendered into the corresponding color, and then the entire ultrasound image after color rendering is obtained, and the ultrasound image is displayed on the display screen of the ultrasound device.

[0148] The blood flow parameter display method provided by the embodiments of the present application constructs four sub-sequences based on the echo data of multiple pulses after wall filtering of the imaging region, calculates the fourth-order cumulant of the echo data using the constructed sub-sequences, and uses the characteristic that the fourth-order cumulant is not sensitive to additive Gaussian color noise and multiplicative noise to eliminate the problem that the autocorrelation method is affected by the accuracy of the blood flow parameter in the additive Gaussian color noise and multiplicative noise background, thereby improving the accuracy of the blood flow parameter estimation and improving the display effect of the blood flow.

[0149] It should be noted that the method for calculating the blood flow parameter in the blood flow parameter display method provided by the embodiments of the present application can be mutually referred to with the method for determining the blood flow parameter provided by the embodiments of the present application. That is, the blood flow parameter can be determined by referring to each embodiment of the method for determining the blood flow parameter provided by the embodiments of the present application, and then the determined blood flow parameter is displayed. For the sake of brevity, it will not be repeated here.

[0150] A specific embodiment is introduced below. In the CFM mode, the ultrasound device transmits ultrasound pulses at a certain echo frequency through the kidney. The ultrasound pulses generate echo signals through reflection of the detected object, and are received by the ultrasound probe of the ultrasound device. After processing the echo signals received by the ultrasound probe, the ultrasound device obtains echo data of multiple pulses of the wall-filtered pixel points in the imaging region, and sends the echo data to the electronic device. The blood flow parameters of each pixel point, including blood flow energy, blood flow velocity and blood flow variance, are calculated in the electronic device, and the specific process is as follows:

[0151] Selecting a delay parameter Or Based on echo data Constructing four sub-sequences

[0152]

[0153] According to the four constructed sub-sequences Calculate the fourth-order cumulant of the echo data

[0154]

[0155] Using the echo data of each pixel point And the fourth-order cumulant Estimate the blood flow energy Power, blood flow velocity v and blood flow variance Variance:

[0156]

[0157]

[0158]

[0159]

[0160] After the electronic device calculates the blood flow parameters of each pixel point, the blood vessel color or contrast of each pixel point in the ultrasound image is adjusted based on the blood flow parameters of each pixel point, and the adjusted ultrasound image is returned to the ultrasound device for display.

[0161] A blood flow parameter determination device provided by an embodiment of the present application is introduced below. The blood flow parameter determination device described below can be referred to in conjunction with the blood flow parameter determination method described above.

[0162] Referring to Figure 5 , a structural diagram of a blood flow parameter determination device according to an exemplary embodiment is shown, as shown in Figure 5 , comprising:​

[0163] The acquisition module 501 is configured to acquire echo data of multiple pulses of a pixel point after wall filtering in an imaging region;

[0164] The first calculation module 502 is configured to construct four sub-sequences based on the echo data, and calculate fourth-order cumulants of the echo data based on the four sub-sequences.

[0165] The second calculation module 503 is configured to calculate a blood flow parameter based on the echo data and the fourth-order cumulants.

[0166] The blood flow parameter determination device provided by the embodiments of the present application constructs four sub-sequences based on echo data of multiple pulses after wall filtering in an imaging region, calculates fourth-order cumulants of the echo data by using the constructed sub-sequences, and eliminates the problem of affecting the accuracy of blood flow parameter estimation by using the characteristics that the fourth-order cumulants are not sensitive to additive Gaussian color noise and multiplicative noise, thereby improving the accuracy of blood flow parameter estimation, effectively ensuring the stability of estimation in a blood vessel edge and a low-speed blood flow region, reducing the possibility of sudden state, and improving the continuity and hierarchy of blood flow imaging. In addition, since the embodiments of the present application eliminate the influence of multiplicative noise, the robustness of estimating the blood flow parameter in a small sample is higher.

[0167] On the basis of the above-mentioned embodiments, as a preferred implementation manner, the first calculation module 502 comprises:

[0168] The construction unit is configured to perform delay processing on the echo data according to a preset delay rule, so as to construct four sub-sequences.

[0169] On the basis of the above-mentioned embodiments, as a preferred implementation form, the constructing unit is specifically configured to: determine a first delay parameter, a second delay parameter, and a third delay parameter; wherein a difference between the first delay parameter and a sum of the second delay parameter and the third delay parameter is 1; construct a first subsequence based on the echo data; wherein a first element of the first subsequence is a first element of the echo data, and a length of the first subsequence is a difference between a length of the echo data and a target maximum value, the target maximum value being a maximum value of the second delay parameter and the third delay parameter; perform delay processing on the echo data based on the first delay parameter to construct a second subsequence; wherein a length of the second subsequence is a difference between the length of the echo data and a first target difference value, the first target difference value being a difference between the target maximum value and the first delay parameter; perform delay processing on the echo data based on the second delay parameter to construct a third subsequence; wherein a length of the third subsequence is a difference between the length of the echo data and a second target difference value, the second target difference value being a difference between the target maximum value and the second delay parameter; and perform delay processing on the echo data based on the third delay parameter to construct a fourth subsequence; wherein a length of the fourth subsequence is a difference between the length of the echo data and a third target difference value, the third target difference value being a difference between the target maximum value and the third delay parameter.

[0170] On the basis of the above-mentioned embodiments, as a preferred implementation form, the first calculating module 502 comprises:

[0171] The calculating unit is configured to: calculate a first point multiplication result between a conjugate of the first subsequence, a conjugate of the second subsequence, the third subsequence, and the fourth subsequence, accumulate each element in the first point multiplication result to obtain a first accumulated value; calculate a second point multiplication result between the conjugate of the first subsequence and the third subsequence, accumulate each element in the second point multiplication result to obtain a second accumulated value; calculate a third point multiplication result between the conjugate of the second subsequence and the fourth subsequence, accumulate each element in the third point multiplication result to obtain a third accumulated value; calculate a fourth point multiplication result between the conjugate of the first subsequence and the fourth subsequence, accumulate each element in the fourth point multiplication result to obtain a fourth accumulated value; calculate a fifth point multiplication result between the conjugate of the second subsequence and the third subsequence, accumulate each element in the fifth point multiplication result to obtain a fifth accumulated value; calculate a first product between the first accumulated value and a difference between the second accumulated value and the third accumulated value, and a second product between the fourth accumulated value and the fifth accumulated value; and take a difference between the first accumulated value and the first product and the second product as the fourth-order cumulant of the echo data.

[0172] In the above embodiment, as a preferred implementation, the blood flow parameter comprises any one or a combination of blood flow energy, blood flow velocity and blood flow variance.

[0173] The second calculation module 503 is specifically configured to: take a product between the echo data and a conjugate of the echo data as an intermediate sequence, accumulate each element in the intermediate sequence to obtain the blood flow energy; and / or, perform angle operation on a negative number of the fourth-order cumulant to obtain a pulse frequency, and calculate the blood flow velocity according to a corresponding relationship between the pulse frequency and the blood flow velocity; and / or, calculate a fourth-order absolute value cumulant based on the four sub-sequences, calculate a ratio between an absolute value of the fourth-order cumulant and the fourth-order absolute value cumulant, and take a difference between 1 and the ratio as the blood flow variance.

[0174] In the above embodiment, as a preferred implementation, the second calculation module 403 is specifically configured to: calculate a first point multiplication result between a conjugate of the first sub-sequence, a conjugate of the second sub-sequence, the third sub-sequence and the fourth sub-sequence, accumulate an absolute value of each element in the first point multiplication result to obtain a sixth accumulated value; calculate a second point multiplication result between the conjugate of the first sub-sequence and the third sub-sequence, accumulate an absolute value of each element in the second point multiplication result to obtain a seventh accumulated value; calculate a third point multiplication result between the conjugate of the second sub-sequence and the fourth sub-sequence, accumulate an absolute value of each element in the third point multiplication result to obtain an eighth accumulated value; calculate a fourth point multiplication result between the conjugate of the first sub-sequence and the fourth sub-sequence, accumulate an absolute value of each element in the fourth point multiplication result to obtain a ninth accumulated value; calculate a fifth point multiplication result between the conjugate of the second sub-sequence and the third sub-sequence, accumulate an absolute value of each element in the fifth point multiplication result to obtain a tenth accumulated value; calculate a third product between the seventh accumulated value and the eighth accumulated value, and a fourth product between the ninth accumulated value and the tenth accumulated value; and take a sum of the sixth accumulated value, the third product and the fourth product as a fourth-order absolute value cumulant.

[0175] As to the apparatus in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiments of the method, and will not be described in detail here.

[0176] The following describes a blood flow parameter display apparatus provided by the embodiments of the present application. The blood flow parameter display apparatus described below can be referred to the blood flow parameter display method described above.

[0177] Reference is made to Figure 6, a structural diagram of a blood flow parameter display device according to an example embodiment is shown, as shown in Figure 6 , including:

[0178] The acquisition module 601 is configured to acquire echo data of multiple pulses of a wall-filtered pixel point in an imaging region;

[0179] The first calculation module 602 is configured to construct four sub-sequences based on the echo data, and calculate a fourth-order cumulant of the echo data based on the four sub-sequences;

[0180] The second calculation module 603 is configured to calculate a blood flow parameter based on the echo data and the fourth-order cumulant;

[0181] The display module 604 is configured to display the blood flow parameter.

[0182] The blood flow parameter display device provided by the example embodiment is based on the echo data of multiple pulses of a wall-filtered pixel point in an imaging region to construct four sub-sequences, and uses the constructed sub-sequences to calculate a fourth-order cumulant of the echo data. The fourth-order cumulant is not sensitive to additive Gaussian color noise and multiplicative noise, which eliminates the problem of affecting the accuracy of blood flow parameter estimation by using autocorrelation methods in the background of additive Gaussian color noise and multiplicative noise, improves the accuracy of blood flow parameter estimation, effectively ensures the stability of estimation in the vascular edge and low-speed blood flow area, reduces the possibility of sudden state, improves the continuity and hierarchy of blood flow imaging, and improves the display effect of blood flow.

[0183] Based on the hardware implementation of the above program modules, and in order to implement the method of the example embodiment, the example embodiment further provides an electronic device, Figure 7 , a structural diagram of an electronic device according to an example embodiment is shown, as shown in Figure 7 , the electronic device includes:

[0184] The communication interface 1 can interact with other devices such as network devices and the like.

[0185] The processor 2 is connected with the communication interface 1 to realize information interaction with other devices, and is used to run a computer program to execute the blood flow parameter determination method provided by one or more technical solutions. The computer program is stored on the memory 3.

[0186] Of course, in actual application, each component in the electronic device is coupled together through a bus system 4. It can be understood that the bus system 4 is used to realize the connection and communication between the components. The bus system 4 includes a data bus, a power bus, a control bus and a state signal bus. However, in order to clearly illustrate, all kinds of buses are marked as the bus system 4 in Figure 7 .

[0187] The memory 3 in the embodiments of the present application is used to store various types of data to support the operation of the electronic device. Examples of these data include: any computer programs used to operate on the electronic device.

[0188] It can be understood that the memory 3 can be a volatile memory or a non-volatile memory, and can also include both volatile and non-volatile memories. Among them, the non-volatile memory can be a read-only memory (ROM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), an electrically erasable programmable read-only memory (EEPROM), a ferromagnetic random access memory (FRAM), a flash memory, a magnetic surface memory, an optical disc, or a compact disc read-only memory (CD-ROM). The magnetic surface memory can be a disk memory or a tape memory. The volatile memory can be a random access memory (RAM) used as an external cache. By way of example but not limitation, many forms of RAM can be used, such as static random access memory (SRAM), synchronous static random access memory (SSRAM), dynamic random access memory (DRAM), synchronous dynamic random access memory (SDRAM), double data rate synchronous dynamic random access memory (DDR SDRAM), enhanced synchronous dynamic random access memory (ESDRAM), sync link dynamic random access memory (SLDRAM), and direct rambus random access memory (DRRAM).The memory 3 described in the embodiments of the present application is intended to include, but not limited to, these and any other suitable types of memory.

[0189] The method disclosed in the embodiments of the present application can be applied in the processor 2 or implemented by the processor 2. The processor 2 can be an integrated circuit chip with processing capability. In the implementation process, each step of the above method can be completed by the integrated logic circuit of hardware in the processor 2 or the instruction in the form of software. The processor 2 described above can be a general processor, a DSP, or other programmable logic device, discrete gate or transistor logic device, discrete hardware component, etc. The processor 2 can implement or execute the disclosed methods, steps and logic block diagrams in the embodiments of the present application. The general processor can be a microprocessor or any conventional processor, etc. In combination with the steps of the method disclosed in the embodiments of the present application, the execution can be directly embodied as hardware decoding processor or executed by the combination of hardware and software modules in the decoding processor. The software module can be located in the storage medium, which is located in the memory 3. The processor 2 reads the program in the memory 3 and combines the hardware to complete the steps of the above method.

[0190] The processor 2 implements the corresponding flow in each method of the embodiments of the present application when executing the program. For brevity, it will not be repeated here.

[0191] In the exemplary embodiments, the embodiments of the present application also provide a storage medium, i.e. a computer storage medium, specifically a computer readable storage medium, for example, including the memory 3 storing the computer program, the above computer program can be executed by the processor 2 to complete the steps of the above method. The computer readable storage medium can be FRAM, ROM, PROM, EPROM, EEPROM, Flash Memory, magnetic surface memory, optical disc, or CD-ROM, etc.

[0192] Those skilled in the art can understand that all or part of the steps of the above method embodiments can be completed by program instruction related hardware, and the above program can be stored in a computer readable storage medium. The program is executed to perform the steps of the above method embodiments; and the above storage medium includes mobile storage device, ROM, RAM, magnetic disc or optical disc, etc. various program code storage media.

[0193] Alternatively, the above-mentioned integrated units of the present application, if realized in the form of software function modules and sold or used as independent products, can also be stored in a computer-readable storage medium. Based on such understanding, the technical solutions of the embodiments of the present application can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a plurality of instructions for causing an electronic device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the methods described in the embodiments of the present application. The aforementioned storage medium includes: mobile storage devices, ROM, RAM, magnetic disks or optical disks, and various media that can store program codes.

[0194] The above is only a specific implementation of the present application, but the protection scope of the present application is not limited thereto, any person skilled in the art can easily think of changes or replacements within the technical scope disclosed in the present application, which should be covered within the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.

Claims

1. A method of determining a blood flow parameter, characterized by, The method comprises: acquiring echo data of multiple pulses of a wall-filtered pixel point in an imaging region; constructing four sub-sequences based on the echo data, and calculating fourth-order cumulants of the echo data based on the four sub-sequences; calculating a blood flow parameter based on the echo data and the fourth-order cumulants; wherein the blood flow parameter comprises a blood flow velocity and / or a blood flow variance; the calculation of the blood flow parameter based on the echo data and the fourth-order cumulants comprises: performing an angle operation on a negative number of the fourth-order cumulants to obtain a pulse frequency, and calculating the blood flow velocity according to a corresponding relationship between the pulse frequency and the blood flow velocity; and / or, calculating fourth-order absolute value cumulants based on the four sub-sequences, calculating a ratio of an absolute value of the fourth-order cumulants to the fourth-order absolute value cumulants, and taking a difference between 1 and the ratio as the blood flow variance.

2. The method of claim 1, wherein the blood flow parameter is determined by: the construction of the four sub-sequences based on the echo data comprises: delaying the echo data according to a preset delay rule to construct the four sub-sequences.

3. The method of claim 2, wherein the blood flow parameter is determined by: the delaying of the echo data according to the preset delay rule to construct the four sub-sequences comprises: determining a first delay parameter, a second delay parameter, and a third delay parameter; wherein a difference between the sum of the second delay parameter and the third delay parameter and the first delay parameter is 1; constructing a first sub-sequence based on the echo data; wherein a first element of the first sub-sequence is a first element of the echo data, a length of the first sub-sequence is a difference between a length of the echo data and a target maximum value, and the target maximum value is a maximum value of the second delay parameter and the third delay parameter; delaying the echo data according to the first delay parameter to construct a second sub-sequence; wherein a length of the second sub-sequence is a difference between the length of the echo data and a first target difference value, and the first target difference value is a difference between the target maximum value and the first delay parameter; delaying the echo data according to the second delay parameter to construct a third sub-sequence; wherein a length of the third sub-sequence is a difference between the length of the echo data and a second target difference value, and the second target difference value is a difference between the target maximum value and the second delay parameter; delaying the echo data according to the third delay parameter to construct a fourth sub-sequence; wherein a length of the fourth sub-sequence is a difference between the length of the echo data and a third target difference value, and the third target difference value is a difference between the target maximum value and the third delay parameter.

4. The method of claim 1, wherein the blood flow parameter is determined by: ###00003### The four sub-sequences comprise a first sub-sequence, a second sub-sequence, a third sub-sequence, and a fourth sub-sequence. ​ the calculation of the fourth-order cumulants of the echo data based on the four sub-sequences comprises: calculating a conjugate of the first sub-sequence, a conjugate of the second sub-sequence, a first point multiplication result among the third sub-sequence, the fourth sub-sequence, and accumulating each element in the first point multiplication result to obtain a first accumulated value; Calculate the second dot product between the conjugate of the first subsequence and the third subsequence, and sum each element in the second dot product to obtain a second accumulated value; Calculate the third dot product between the conjugate of the second subsequence and the fourth subsequence, and sum each element in the third dot product to obtain the third accumulated value; Calculate the conjugate of the first subsequence and the fourth dot product between the fourth subsequences, and sum each element in the fourth dot product to obtain the fourth accumulated value; Calculate the fifth dot product between the conjugate of the second subsequence and the third subsequence, and sum each element in the fifth dot product to obtain the fifth accumulated value; Calculate the first product between the second accumulated value and the third accumulated value, and the second product between the fourth accumulated value and the fifth accumulated value; The difference between the first accumulated value and the first product and the second product is taken as the fourth-order cumulant of the echo data.

5. The method of claim 1, wherein the blood flow parameter is determined by: ###00005### The blood flow parameters also include blood flow energy; ​ The calculation of blood flow parameters based on the echo data and the fourth-order cumulative quantity also includes: The blood flow energy is obtained by accumulating the product of the echo data and the conjugate of the echo data as an intermediate sequence.

6. The method of determining a blood flow parameter according to claim 5, wherein, The four subsequences include a first subsequence, a second subsequence, a third subsequence, and a fourth subsequence; The calculation of the fourth-order absolute value cumulant based on the four subsequences includes: Calculate the first dot product result between the conjugate of the first subsequence, the conjugate of the second subsequence, the third subsequence, and the fourth subsequence, and sum the absolute values ​​of each element in the first dot product result to obtain the sixth accumulated value; Calculate the second dot product between the conjugate of the first subsequence and the third subsequence, and sum the absolute values ​​of each element in the second dot product to obtain a seventh accumulated value; Calculate the conjugate of the second subsequence and the third dot product between the fourth subsequence, and sum the absolute values ​​of each element in the third dot product to obtain the eighth accumulated value; Calculate the conjugate of the first subsequence and the fourth dot product between the fourth subsequences, and sum the absolute values ​​of each element in the fourth dot product to obtain the ninth accumulated value; Calculate the fifth dot product between the conjugate of the second subsequence and the third subsequence, and sum the absolute values ​​of each element in the fifth dot product to obtain the tenth accumulated value; Calculate the third product between the seventh accumulated value and the eighth accumulated value, and the fourth product between the ninth accumulated value and the tenth accumulated value; The sum of the sixth accumulated value, the third product, and the fourth product is taken as the fourth-order absolute value accumulator.

7. A blood flow parameter display method characterized by comprising: include: Acquire echo data of multiple pulses from pixels after wall filtering in the imaging region; Four subsequences are constructed based on the echo data, and the fourth-order cumulant of the echo data is calculated based on the four subsequences. Blood flow parameters are calculated based on the echo data and the fourth-order cumulative volume. Display the blood flow parameters; The blood flow parameter comprises a blood flow velocity and / or a blood flow variance. The blood flow parameter is calculated based on the echo data and the fourth-order cumulant, comprising: The negative of the fourth-order cumulant is subjected to an angle operation to obtain a pulse frequency, and the blood flow velocity is calculated according to a corresponding relationship between the pulse frequency and the blood flow velocity. And / or, a fourth-order absolute value cumulant is calculated based on the four sub-sequences, a ratio of an absolute value of the fourth-order cumulant to the fourth-order absolute value cumulant is calculated, and a difference between 1 and the ratio is taken as the blood flow variance.

8. An apparatus for determining a blood flow parameter, characterized in that Comprising: An acquisition module is configured to acquire echo data of multiple pulses of a pixel point after wall filtering in an imaging region; A first calculation module is configured to construct four sub-sequences based on the echo data, and calculate a fourth-order cumulant of the echo data based on the four sub-sequences; A second calculation module is configured to calculate a blood flow parameter based on the echo data and the fourth-order cumulant; The blood flow parameter comprises a blood flow velocity and / or a blood flow variance. The second calculation module is specifically configured to: subject the negative of the fourth-order cumulant to an angle operation to obtain a pulse frequency, and calculate the blood flow velocity according to a corresponding relationship between the pulse frequency and the blood flow velocity; and / or, calculate a fourth-order absolute value cumulant based on the four sub-sequences, calculate a ratio of an absolute value of the fourth-order cumulant to the fourth-order absolute value cumulant, and take a difference between 1 and the ratio as the blood flow variance.

9. A blood flow parameter display device, characterized by comprising: Comprising: An acquisition module is configured to acquire echo data of multiple pulses of a pixel point after wall filtering in an imaging region; A first calculation module is configured to construct four sub-sequences based on the echo data, and calculate a fourth-order cumulant of the echo data based on the four sub-sequences; A second calculation module is configured to calculate a blood flow parameter based on the echo data and the fourth-order cumulant; A display module is configured to display the blood flow parameter; The blood flow parameter comprises a blood flow velocity and / or a blood flow variance. The second calculation module is specifically configured to: subject the negative of the fourth-order cumulant to an angle operation to obtain a pulse frequency, and calculate the blood flow velocity according to a corresponding relationship between the pulse frequency and the blood flow velocity; and / or, calculate a fourth-order absolute value cumulant based on the four sub-sequences, calculate a ratio of an absolute value of the fourth-order cumulant to the fourth-order absolute value cumulant, and take a difference between 1 and the ratio as the blood flow variance.

10. An electronic device, comprising: Comprising: A memory is configured to store a computer program; A processor is configured to implement steps of the method in any one of claims 1 to 7 when the computer program is executed.

11. A computer readable storage medium, characterized in that, The computer program is stored on the computer readable storage medium, and the computer program is executed by the processor to implement steps of the method in any one of claims 1 to 7.

Images