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; generating two sub-sequences based on the echo data, and calculating third-order cumulants of the two sub-sequences respectively; calculating cross-correlation functions of the two sub-sequences by using the third-order cumulants of the two sub-sequences; and calculating blood flow parameters based on the echo data and the cross-correlation functions. In the application, the echo data of multiple pulses of wall-filtered pixel points in the imaging region is divided into two sub-sequences, the third-order cumulants of the sub-sequences are calculated, and the characteristics that the third-order cumulants are not sensitive to Gaussian color noise are used to eliminate the influence of Gaussian color noise on the accuracy of the blood flow parameter estimated by the autocorrelation method, thereby improving the accuracy of blood flow parameter estimation.

Description

Technical Field

[0001] This application relates to the field of ultrasound technology, and more specifically, to a method and apparatus for determining and displaying blood flow parameters, an electronic device, and a computer-readable storage medium. Background Technology

[0002] Currently, color Doppler imaging systems use autocorrelation methods to estimate blood flow parameters. However, the acquired images contain Gaussian white noise, which becomes Gaussian colored noise after wall filtering. The autocorrelation function of the Gaussian colored noise is no longer the impulse function, which seriously affects the accuracy of blood flow parameter estimation.

[0003] Therefore, improving the accuracy of blood flow parameter estimation is a technical problem that needs to be solved by those skilled in the art. Summary of the Invention

[0004] The purpose of this application is to provide a method, apparatus, electronic device, and computer-readable storage medium for determining and displaying blood flow parameters, thereby improving the accuracy of blood flow parameter estimation.

[0005] To achieve the above objectives, this application provides a method for determining blood flow parameters, comprising:

[0006] Acquire echo data of multiple pulses from pixels after wall filtering in the imaging region;

[0007] Two subsequences are generated based on the echo data, and the third-order cumulants of the two subsequences are calculated respectively;

[0008] The cross-correlation function of the two subsequences is calculated using the third-order cumulants of the two subsequences;

[0009] Blood flow parameters are calculated based on the echo data and the cross-correlation function.

[0010] The generation of two sub-sequences based on the echo data includes:

[0011] A first subsequence is constructed based on the echo data; wherein the first element of the first subsequence is the first element of the echo data, and the length of the first subsequence is the length of the echo data minus one.

[0012] The echo data is delayed by one step to generate a second subsequence; wherein the first element of the second subsequence is the second element of the echo data, and the length of the second subsequence is equal to that of the first subsequence.

[0013] The calculation of the third-order cumulants of the two subsequences includes:

[0014] Determine the first delay parameter, the second delay parameter, and the third delay parameter;

[0015] Based on the first delay parameter, the second delay parameter, and the third delay parameter, the two sub-sequences are delayed respectively to obtain the first delayed sequence, the second delayed sequence, and the third delayed sequence corresponding to the two sub-sequences respectively;

[0016] Perform conjugate operation on at least one delayed sequence corresponding to each of the two subsequences to obtain at least one conjugate operation result corresponding to each of the two subsequences. Use the dot product between the conjugate operation result corresponding to each of the two subsequences and other delayed sequences that have not undergone conjugate operation as the third-order cumulants of the two subsequences.

[0017] Wherein, the third delay parameter is zero, and correspondingly, the step of performing conjugation operations on at least one delay sequence corresponding to each of the two subsequences to obtain at least one conjugation operation result corresponding to each of the two subsequences, and using the dot product between the conjugation operation result corresponding to each of the two subsequences and other delay sequences that have not undergone conjugation operations as the third-order cumulants of the two subsequences, includes:

[0018] Perform conjugate operation on the third delay sequence corresponding to each of the two subsequences to obtain the conjugate operation results corresponding to each of the two subsequences. Use the conjugate operation results corresponding to each of the two subsequences, the dot product between the first delay sequence and the second delay sequence as the third-order cumulants of the two subsequences.

[0019] Accordingly, the step of calculating the cross-correlation function of the two subsequences using the third-order cumulants of the two subsequences includes:

[0020] The dot product between the conjugate of the third-order cumulants of the first subsequence and the third-order cumulants of the second subsequence is used as the cross-correlation function.

[0021] Wherein, the second delay parameter is equal to the third delay parameter, and correspondingly, the step of performing conjugation operation on at least one delay sequence corresponding to each of the two subsequences to obtain at least one conjugation operation result corresponding to each of the two subsequences, and using the dot product result between the conjugation operation result corresponding to each of the two subsequences and other delay sequences that have not undergone conjugation operation as the third-order cumulants of the two subsequences, includes:

[0022] Perform conjugate operations on the first and third delayed sequences corresponding to the two subsequences to obtain the first and second conjugate operation results corresponding to the two subsequences respectively. Use the first conjugate operation result, the second conjugate operation result, and the dot product result between the second delayed sequences corresponding to the two subsequences as the third-order cumulants of the two subsequences respectively.

[0023] Accordingly, the step of calculating the cross-correlation function of the two subsequences using the third-order cumulants of the two subsequences includes:

[0024] The dot product between the conjugate of the third-order cumulant of the second subsequence and the third-order cumulant of the first subsequence is used as the cross-correlation function.

[0025] Wherein, the two subsequences are the first subsequence and the second subsequence, and the calculation of the third-order cumulant of the subsequences includes:

[0026] Determine the first delay parameter, the second delay parameter, and the third delay parameter;

[0027] Based on the first delay parameter, the second delay parameter, and the third delay parameter, the first subsequence is delayed to obtain a first delayed sequence, a second delayed sequence, and a third delayed sequence.

[0028] The second subsequence is delayed based on the first delay parameter to obtain a fourth delayed sequence;

[0029] Perform a conjugate operation on the third delayed sequence, and use the dot product of the conjugate of the third delayed sequence, the first delayed sequence, and the second delayed sequence as the third-order cumulant of the first subsequence;

[0030] The dot product of the conjugate of the third delay sequence, the fourth delay sequence, and the second delay sequence is used as the third-order cumulant of the second subsequence;

[0031] Accordingly, the step of calculating the cross-correlation function of the two subsequences using the third-order cumulants of the two subsequences includes:

[0032] The dot product between the conjugate of the third-order cumulants of the first subsequence and the cross-third-order cumulants of the second subsequence is used as the cross-correlation function.

[0033] The blood flow parameters include any one or a combination of blood flow energy, blood flow velocity, and blood flow variance.

[0034] The calculation of blood flow parameters based on the echo data and the cross-correlation function includes:

[0035] 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;

[0036] And / or, each element in the cross-correlation function is accumulated to obtain a first accumulated value, the first accumulated value is used to perform angle calculation to obtain the pulse frequency, and the blood flow velocity is calculated according to the correspondence between the pulse frequency and the blood flow velocity;

[0037] And / or, summing up each element in the cross-correlation function to obtain a first accumulated value, summing up the absolute values ​​of each element in the cross-correlation function to obtain a second accumulated value, calculating the ratio of the absolute value of the first accumulated value to the second accumulated value, and taking the difference between 1 and the ratio as the blood flow variance.

[0038] To achieve the above objectives, this application provides a method for displaying blood flow parameters, comprising:

[0039] Acquire echo data of multiple pulses from pixels after wall filtering in the imaging region;

[0040] Two subsequences are generated based on the echo data, and the third-order cumulants of the two subsequences are calculated respectively;

[0041] The cross-correlation function of the two subsequences is calculated using the third-order cumulants of the two subsequences;

[0042] Blood flow parameters are calculated based on the echo data and the cross-correlation function;

[0043] The blood flow parameters are displayed.

[0044] To achieve the above objectives, this application provides a device for determining blood flow parameters, comprising:

[0045] The acquisition module is used to acquire echo data of multiple pulses from pixels after wall filtering in the imaging region;

[0046] The first calculation module is used to generate two sub-sequences based on the echo data, and to calculate the third-order cumulants of the two sub-sequences respectively.

[0047] The second calculation module is used to calculate the cross-correlation function of the two subsequences using the third-order cumulants of the two subsequences;

[0048] The third calculation module is used to calculate blood flow parameters based on the echo data and the cross-correlation function.

[0049] To achieve the above objectives, this application provides a blood flow parameter display device, comprising:

[0050] The acquisition module is used to acquire echo data of multiple pulses from pixels after wall filtering in the imaging region;

[0051] The first calculation module is used to generate two sub-sequences based on the echo data, and to calculate the third-order cumulants of the two sub-sequences respectively.

[0052] The second calculation module is used to calculate the cross-correlation function of the two subsequences using the third-order cumulants of the two subsequences;

[0053] The third calculation module is used to calculate blood flow parameters based on the echo data and the cross-correlation function;

[0054] The display module is used to display the blood flow parameters.

[0055] To achieve the above objectives, this application provides an electronic device, comprising:

[0056] Memory, used to store computer programs;

[0057] A processor, used to execute the computer program to implement the steps of the method for determining blood flow parameters as described above.

[0058] To achieve the above objectives, this application provides an electronic device, comprising:

[0059] Memory, used to store computer programs;

[0060] A processor is used to execute the computer program to implement the steps of the blood flow parameter display method described above.

[0061] To achieve the above objectives, this application provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of the method for determining blood flow parameters as described above.

[0062] To achieve the above objectives, this application provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of the blood flow parameter display method described above.

[0063] As can be seen from the above scheme, the method for determining blood flow parameters provided in this application includes: acquiring echo data of multiple pulses of pixels after wall filtering in the imaging region; generating two sub-sequences based on the echo data, and calculating the third-order cumulants of the two sub-sequences respectively; calculating the cross-correlation function of the two sub-sequences using the third-order cumulants of the two sub-sequences; and calculating blood flow parameters based on the echo data and the cross-correlation function.

[0064] The method for determining blood flow parameters provided in this application divides the echo data of the imaging region after wall filtering and multiple pulses into two subsequences, solves the third-order cumulant of each subsequence, and utilizes the characteristic that the third-order cumulant is insensitive to Gaussian colored noise to eliminate the problem of Gaussian colored noise affecting the accuracy of blood flow parameter estimation by autocorrelation method, thereby improving the accuracy of blood flow parameter estimation.

[0065] This application also discloses a method for displaying blood flow parameters, a device for determining blood flow parameters, a device for displaying blood flow parameters, an electronic device, and a computer-readable storage medium, which can achieve the same technical effects as described above.

[0066] It should be understood that the above general description and the following detailed description are merely exemplary and do not limit this application. Attached Figure Description

[0067] 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:

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

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

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

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

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

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

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

[0075] 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.

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

[0077] This application scenario includes an ultrasound device and an electronic device, both connected via a network. The ultrasound device can be of various types, such as B-mode ultrasound, color Doppler ultrasound, endoscopic Doppler ultrasound, or intravascular ultrasound. The electronic device is capable of processing ultrasound data; it can be a terminal device or a server. When it's a terminal device, it can be a personal computer, tablet, or smartphone. In the specific application scenario, the ultrasound device is equipped with an ultrasound probe. The ultrasound device emits ultrasound pulses at a specific frequency towards the object being examined (e.g., a blood vessel). These ultrasound pulses are reflected by the object, generating echo signals, which are received by the ultrasound probe. The ultrasound device processes the echo signals received by the ultrasound probe to obtain echo data of multiple pulses from pixels in the imaging area after wall filtering, and sends this echo data to the electronic device. The electronic device processes the acquired echo data to calculate blood flow parameters.

[0078] In some applications, electronic devices can return the calculated blood flow parameters to ultrasound equipment for display.

[0079] This application discloses a method for determining blood flow parameters, which improves the accuracy of blood flow parameter estimation.

[0080] See Figure 1 The flowchart illustrates a method for determining blood flow parameters according to an exemplary embodiment, wherein the subject executing the method can be an electronic device. Figure 1 As shown, the method includes:

[0081] S101: Acquire echo data of multiple pulses from pixels after wall filtering in the imaging region;

[0082] The CFM (Color Doppler) mode of ultrasound equipment utilizes the Doppler frequency shift generated by multiple ultrasound echoes due to blood flow motion to indirectly describe blood flow motion by estimating the Doppler frequency. In practice, the complex data x acquired by the ultrasound equipment in CFM mode after wall filtering of the imaging area is used. L×P×En Among them, such as Figure 2 As shown, L stands for Line, P for Point, and En for Ensemble. (L, P) represents a pixel in the complex data, and En is the number of pulse repetitions. This method acquires the echo data of multiple pulses for each pixel (L, P). Complex data x after wall filtering of imaging region L×P×En It is a three-dimensional complex matrix, that is, data composed of three dimensions: L, P, and En. The data in the En dimension can be called En-directed data, which is used to estimate the blood flow Doppler frequency and further obtain the blood flow velocity.

[0083] In this embodiment, we assume an additive noise echo model Y(t) = x(t) + n(t), where Y(t) represents the acquired echo data, x(t) represents the actual echo data, and n(t) represents additive noise. That is, additive noise may be present in the acquired echo data.

[0084] Optionally, the ultrasound device acquires real data of the imaging area (which can be an area containing blood vessels) through the ultrasound probe, demodulates the real data to obtain complex data, and then performs wall filtering to filter out data with small motion trends such as blood vessel walls. At this point, the wall-filtered complex data of the imaging area is obtained.

[0085] Optionally, the imaging area can correspond to many pixels. In the specific implementation, one or more pixels can be polled in turn, and the corresponding echo data can be obtained for each polled pixel, and blood flow parameters can be calculated until the blood flow parameters of all pixels in the imaging area have been calculated. Alternatively, echo data can be obtained on a pixel-by-pixel basis, and then the blood flow parameters can be calculated for the entire imaging area.

[0086] It should be noted that since the echo data is obtained based on multiple pulses, each pulse's echo signal corresponds to one data element. Arranging the data elements corresponding to different pulses yields a data sequence. Therefore, the echo data corresponding to all pulses can be a sequence composed of multiple data elements.

[0087] S102: Generate two sub-sequences based on the echo data, and calculate the third-order cumulants of the two sub-sequences respectively;

[0088] In this step, the echo data of multiple pulses for each pixel (L, P) are first processed. Divided into two subsequences The length of a subsequence can be equal to or shorter than the length of the echo data. If shorter, it can be one element shorter, two elements shorter, etc. Furthermore, a subsequence can be defined as a continuous segment of elements in the echo data, or it can be defined as a subsequence by selecting elements at intervals. Optionally, the echo data can be split into two subsequences of the same length but with different elements. For example, for echo data containing N elements, the sequence of N-1 elements (1 to N-1) can be defined as one subsequence, and the sequence of N-1 elements (2 to N) can be defined as the other subsequence.

[0089] As a feasible implementation, generating two sub-sequences based on the echo data includes: 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 length of the echo data minus one; and generating a second sub-sequence by performing a one-step delay on the echo data; wherein the first element of the second sub-sequence is the second element of the echo data, and the length of the second sub-sequence is equal to that of the first sub-sequence. The specific implementation is as follows:

[0090]

[0091] Wherein, ":" indicates that the sequence index increments sequentially, starting with 1 (it should be noted that in some embodiments, it can also increment by other numbers, such as 2, 3, etc.). For the first subsequence, This is the second subsequence.

[0092] Secondly, estimate the third-order cumulants of the two subsequences (the first subsequence and the second subsequence) separately. and The third-order cumulant is the result of the dot product of three sequences.

[0093] Optionally, calculating the third-order cumulants of the two subsequences includes: determining a first delay parameter, a second delay parameter, and a third delay parameter; performing delay processing on the two subsequences based on the first delay parameter, the second delay parameter, and the third delay parameter to obtain a first delay sequence, a second delay sequence, and a third delay sequence corresponding to the two subsequences respectively; performing a conjugation operation on at least one delay sequence corresponding to each of the two subsequences to obtain at least one conjugation operation result corresponding to each of the two subsequences; and using the dot product between the conjugation operation result corresponding to each of the two subsequences and other delay sequences that have not undergone conjugation operation as the third-order cumulants of the two subsequences respectively. In a specific implementation, when calculating the third-order cumulant of a subsequence, at least two delay parameters can be determined first, and the subsequence can be delayed based on the delay parameters to obtain the corresponding delay sequence; the conjugate of at least one delay sequence can be determined, and the conjugate of the determined delay sequence can be multiplied by other delay sequences, and the operation result can be determined as the third-order cumulant of the subsequence. In the embodiments of this application, the method of performing a conjugation operation on one of the delay sequences can be called single conjugation.

[0094] As a feasible implementation, the step of performing conjugation operations on at least one delayed sequence corresponding to each of the two subsequences to obtain at least one conjugation operation result corresponding to each of the two subsequences, and using the dot product between the conjugation operation result corresponding to each of the two subsequences and other delayed sequences that have not undergone conjugation operations as the third-order cumulants of the two subsequences, includes: performing conjugation operations on the third delayed sequence corresponding to each of the two subsequences to obtain the conjugation operation result corresponding to each of the two subsequences, and using the dot product between the conjugation operation result corresponding to each of the two subsequences, the first delayed sequence, and the second delayed sequence as the third-order cumulants of the two subsequences.

[0095] Optionally, when performing the dot product operation, for the first subsequence, the conjugate of the third delay sequence of the first subsequence is determined, and the conjugate of the third delay sequence is multiplied by the first delay sequence and the second delay sequence of the first subsequence to obtain the third-order cumulant of the first subsequence; similarly, for the second subsequence, the conjugate of the third delay sequence of the second subsequence is determined, and the conjugate of the third delay sequence is multiplied by the first delay sequence and the second delay sequence of the second subsequence to obtain the third-order cumulant of the second subsequence.

[0096] The above implementation method is a single conjugate form, and the specific implementation method is as follows:

[0097]

[0098] Among them, τ, γ represents the first delay parameter, the second delay parameter, and the third delay parameter, respectively. These are the first delayed sequence, the second delayed sequence, and the third delayed sequence corresponding to the first subsequence. The conjugate of the third delayed sequence corresponding to the first subsequence. These are the first delayed sequence, the second delayed sequence, and the third delayed sequence corresponding to the second subsequence, respectively. is the conjugate of the third delayed sequence corresponding to the second subsequence, and ⊙ represents the dot product, that is, the multiplication of corresponding elements of the two sequences.

[0099] Preferably, the third delay parameter is zero, i.e., γ = 0. This results in better performance and higher accuracy of the calculated blood flow parameters. The specific implementation is as follows:

[0100]

[0101] Among them, the sequence of dot product results and two sequences The lengths are equal. Because τ, The presence of this reduces the length of the sequence involved in the calculation, effectively improving the efficiency of blood flow parameter computation. To optimize performance, it is possible to make... Furthermore, one of them is zero. To further optimize performance by ensuring long sequences, we can choose τ = 1. Or τ = 0,

[0102] In the above implementation, there are three cases. The first case is that the first delay parameter, the second delay parameter, and the third delay parameter are all equal, i.e. when In this case, the specific implementation method can be represented as:

[0103]

[0104]

[0105]

[0106] Among them, it is necessary to ensure and Consistency in arrangement, which can occur without delay.

[0107] Alternatively, there can be a delay, that is, Compared to no delay, having a delay can reduce the amount of data to be calculated.

[0108] The second case is when the first, second, and third delay parameters are partially equal, i.e. Placing conjugate sequences on sequences with the same delay yields better performance, as shown in the following implementation:

[0109]

[0110] in, The arrangement of the front and back is arbitrary, and Just keep it consistent.

[0111] The third case is when the first delay parameter, the second delay parameter, and the third delay parameter are all unequal, i.e. The performance at this point is not as good as the two situations mentioned above.

[0112] As another feasible implementation, the step of performing conjugation operations on at least one delayed sequence corresponding to each of the two sub-sequences to obtain at least one conjugation operation result corresponding to each of the two sub-sequences, and using the dot product between the conjugation operation result corresponding to each of the two sub-sequences and other delayed sequences that have not undergone conjugation operations as the third-order cumulants of the two sub-sequences, includes: performing conjugation operations on the first delayed sequence and the third delayed sequence corresponding to each of the two sub-sequences to obtain the first conjugation operation result and the second conjugation operation result corresponding to each of the two sub-sequences, and using the dot product between the first conjugation operation result, the second conjugation operation result, and the second delayed sequence corresponding to each of the two sub-sequences as the third-order cumulants of the two sub-sequences. In this embodiment, the method of performing conjugation operations on two delayed sequences can be called double conjugation.

[0113] In this implementation, performance is poor when the first delay parameter, the second delay parameter, and the third delay parameter are all equal or all unequal. Performance is better when the first delay parameter, the second delay parameter, and the third delay parameter are partially equal. The performance is better when the two conjugates are placed on two sequences with different delays. If the second delay parameter is equal to the third delay parameter, i.e. The implementation method is as follows:

[0114]

[0115] in, The arrangement of the front and back is arbitrary, and Just keep it consistent.

[0116] As another feasible implementation, the step of calculating the third-order cumulants of the two subsequences respectively includes: determining a first delay parameter, a second delay parameter, and a third delay parameter; performing delay processing on the first subsequence based on the first delay parameter, the second delay parameter, and the third delay parameter respectively to obtain a first delayed sequence, a second delayed sequence, and a third delayed sequence; performing delay processing on the second subsequence based on the first delay parameter to obtain a fourth delayed sequence; performing a conjugate operation on the third delayed sequence, and using the dot product result between the conjugate of the third delayed sequence, the first delayed sequence, and the second delayed sequence as the third-order cumulant of the first subsequence; and using the dot product result between the conjugate of the third delayed sequence, the fourth delayed sequence, and the second delayed sequence as the cross-third-order cumulant of the second subsequence.

[0117] The above implementation method calculates the mutual tertiary cumulants of two subsequences, and the specific implementation method is as follows:

[0118]

[0119] Among them, the first conjugate and not placed in The performance is better on sequences, τ, The value of γ can be arbitrary, and the sequence can be arranged arbitrarily, ensuring that and As long as the format is consistent, it is acceptable.

[0120] S103: Calculate the cross-correlation function of the two subsequences using the third-order cumulants of the two subsequences;

[0121] In this step, the third-order cumulants of the two subsequences are used. Calculate their cross-correlation to obtain the third-order cumulant cross-correlation function. Optionally, a dot product operation can be performed on the third-order cumulants, and the correlation between the two subsequences can be determined based on the result of the dot product operation.

[0122] For the first conjugate form, the dot product between the conjugate of the third-order cumulants of the first subsequence and the third-order cumulants of the second subsequence is used as the cross-correlation function. The specific implementation is as follows: For the double conjugate form, the dot product between the conjugate of the third-order cumulant of the second subsequence and the third-order cumulant of the first subsequence is used as the cross-correlation function. The specific implementation is as follows: For the cross-third cumulants, the dot product between the conjugate of the third cumulants of the first subsequence and the cross-third cumulants of the second subsequence is used as the cross-correlation function. The specific implementation is as follows:

[0123] S104: Calculate blood flow parameters based on the echo data and the cross-correlation function;

[0124] In this step, echo data from multiple pulses at each pixel (L, P) are used. Cross-correlation function with third-order cumulants Estimate blood flow parameters. Blood flow parameters may include blood flow energy, blood flow velocity, and blood flow variance.

[0125] Optionally, echo data or cross-correlation functions can be calculated separately to obtain the corresponding blood flow parameters. For example, blood flow energy can be calculated based on echo data, blood flow velocity can be calculated based on cross-correlation functions, and blood flow variance can be calculated based on cross-correlation functions. Alternatively, echo data and cross-correlation functions can be calculated together, such as by performing addition, multiplication, or division operations, and the corresponding blood flow parameters can be obtained based on the calculation results.

[0126] The blood flow energy (Power) is calculated as follows: the product of the echo data and the conjugate of the echo data is used as an intermediate sequence, and each element in the intermediate sequence is summed to obtain the blood flow energy. The specific implementation is as follows:

[0127]

[0128] The blood flow velocity (v) is calculated as follows: Each element in the cross-correlation function is summed to obtain a first accumulated value; an angle operation is performed on the first accumulated value to obtain the pulse frequency; and the blood flow velocity is calculated based on the correspondence between the pulse frequency and the blood flow velocity. The specific implementation of the angle calculation is as follows:

[0129] in, This is the first accumulated value;

[0130] It is understandable that, such as Figure 3 As shown, the Doppler frequency of the ultrasonic device is also the echo frequency f. d The relationship with blood flow velocity v is: f d =2vcosθ / λ, where θ represents the angle between the direction of ultrasound and the direction of blood flow velocity, and λ represents the ultrasound carrier wavelength. Therefore, v = f d / (2cosθ / λ). In the above formula, f is obtained after the angle operation angle(). d It can be based on v=f d The blood flow velocity is calculated using / (2cosθ / λ). It should be noted that in actual calculations, since the blood flow direction θ is unknown, the value f is generally used. dThis represents blood flow velocity, which is the blood flow velocity value without blood flow direction estimation. It can be understood as qualitative blood flow.

[0131] The blood flow variance is calculated as follows: Each element in the cross-correlation function is summed to obtain a first accumulated value; the absolute values ​​of each element in the cross-correlation function are summed to obtain a second accumulated value; the ratio of the absolute value of the first accumulated value to the second accumulated value is calculated; and the difference between 1 and the ratio is taken as the blood flow variance. The specific implementation is as follows: in, This is the second accumulated value.

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

[0133] Suppose that the distribution of the random variable x is x ~ N(0, σ). 2 Therefore, the probability density function of x is:

[0134]

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

[0136]

[0137] In the integral formula,

[0138] Combining the two methods, The moment generating function can be obtained as follows:

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

[0140]

[0141]

[0142]

[0143]

[0144] From the relationship between higher-order moments and moment generating functions, we can obtain

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

[0146] By extension, for any integer k, the sum of moments of a Gaussian random variable can be uniformly written as:

[0147]

[0148] The cumulant generating function of a Gaussian random variable x can be directly obtained from the moment generating function:

[0149]

[0150] Its derivatives are:

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

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

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

[0154] From the relationship between higher-order cumulants and cumulant generating functions, we can obtain:

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

[0156] The above results regarding higher-order moments and higher-order cumulants can be generalized to: the second-order moments and second-order cumulants of any zero-mean Gaussian random process are identical, both equal to its variance σ. 2 Its odd-order moments are always zero, but even-order moments are not equal to zero; while higher-order cumulants (third order and above) are always equal to zero. Therefore, in this sense, higher-order cumulants are "blind" to Gaussian random processes, that is, higher-order cumulants are insensitive to Gaussian colored noise.

[0157] The method for determining blood flow parameters provided in this application divides the echo data after multiple pulses of the imaging region wall filtering into two subsequences. The third-order cumulant of each subsequence is calculated. Utilizing the insensitivity of the third-order cumulant to Gaussian noise, the problem of Gaussian noise affecting the accuracy of autocorrelation methods in estimating blood flow parameters is eliminated. Simultaneously, the influence of additive noise in the acquired echo data is effectively removed, and blood flow parameters are calculated based on real echo data, improving the accuracy of blood flow parameter estimation. Furthermore, since there is more Gaussian noise at the vessel edge and in low-velocity blood flow regions, eliminating the influence of Gaussian noise can effectively ensure the stability of the estimation at these regions, reducing the possibility of abrupt changes and improving the coherence and hierarchy of blood flow imaging.

[0158] This application discloses a method for displaying blood flow parameters; for details, please refer to... Figure 4 A flowchart illustrating a method for displaying blood flow parameters according to an exemplary embodiment, as shown below. Figure 4 As shown, it includes:

[0159] S201: Acquire echo data of multiple pulses from pixels after wall filtering in the imaging region;

[0160] S202: Generate two sub-sequences based on the echo data, and calculate the third-order cumulants of the two sub-sequences respectively;

[0161] S203: Calculate the cross-correlation function of the two subsequences using the third-order cumulants of the two subsequences;

[0162] S204: Calculate blood flow parameters based on the echo data and the cross-correlation function;

[0163] S205: Display the blood flow parameters.

[0164] In this step, the blood flow parameters of each pixel are estimated and then displayed. These parameters can be directly displayed on the screen as text, or the blood flow image can be processed based on these parameters, and the processed image can be displayed indirectly.

[0165] It should be noted that the execution entity for the blood flow parameter display method can be an ultrasound device. The ultrasound device acquires echo data, calculates blood flow parameters based on the echo data, and displays the calculated blood flow parameters. Alternatively, the execution entity for the blood flow parameter display method can also be an electronic device. The electronic device acquires echo data from the ultrasound device, calculates blood flow parameters based on the echo data, and displays the calculated blood flow parameters. Of course, the electronic device can also send the calculated blood flow parameters to the ultrasound device for display on the ultrasound device's screen.

[0166] Optionally, based on blood flow parameters, the color and contrast of blood vessels in the ultrasound image can be adjusted, or the ultrasound image can be rendered and then displayed on the screen. For example, for a certain pixel in the ultrasound image, blood flow velocities between (v0, v1) are displayed as orange-red, blood flow velocities between (v1, v2) are displayed as magenta, and blood flow velocities between (v2, v3) are displayed as dark red. According to the blood flow velocity of each pixel, the pixel is rendered with the corresponding color, thus obtaining the entire ultrasound image after color rendering, and the ultrasound image is displayed on the screen of the ultrasound device.

[0167] The blood flow parameter display method provided in this application divides the echo data of the imaging region after wall filtering and multiple pulses into two subsequences, solves the third-order cumulant of each subsequence, and utilizes the characteristic that the third-order cumulant is insensitive to Gaussian noise to eliminate the problem that Gaussian noise affects the accuracy of blood flow parameter estimation by autocorrelation method, thereby improving the accuracy of blood flow parameter estimation and improving the display effect of blood flow.

[0168] It should be noted that the method for calculating blood flow parameters in the blood flow parameter display method provided in this application embodiment can be referred to in conjunction with the method for determining blood flow parameters provided in this application embodiment. That is, the blood flow parameters can be determined by referring to the various embodiments of the blood flow parameter determination method provided in this application embodiment, and then the determined blood flow parameters can be displayed. For the sake of brevity, it will not be described in detail here.

[0169] The following describes a specific implementation method: In CFM mode, the ultrasound device emits ultrasound pulses through the kidney at a specific echo frequency. These ultrasound pulses are reflected by the object being examined, generating echo signals, which are received by the ultrasound probe of the ultrasound device. The ultrasound device processes the echo signals received by the ultrasound probe to obtain echo data of multiple pulses from pixels in the imaging area after wall filtering, and sends this echo data to an electronic device. The electronic device calculates the blood flow parameters for each pixel, including blood flow energy, blood flow velocity, and blood flow variance. The specific process is as follows:

[0170] Echo data based on each pixel Generate two subsequences and

[0171]

[0172] Select delay parameter τ = 1, Or τ = 0, Estimate separately and Third-order cumulative and

[0173]

[0174] Using the third-order cumulants of two subsequences Calculate their cross-correlation function to obtain the third-order cumulant cross-correlation function.

[0175]

[0176] Using echo data of each pixel Cross-correlation function with third-order cumulants Estimate blood flow energy (Power), blood flow velocity (v), and blood flow variance (Variance):

[0177]

[0178]

[0179]

[0180] After the electronic device calculates the blood flow parameters of each pixel, it adjusts the blood vessel color or contrast of each pixel in the ultrasound image based on the blood flow parameters of each pixel, and then returns the adjusted ultrasound image to the ultrasound device for display.

[0181] The following describes a blood flow parameter determination device provided in an embodiment of this application. The blood flow parameter determination device described below and the blood flow parameter determination method described above can be referred to each other.

[0182] See Figure 5 A structural diagram of a blood flow parameter determination device according to an exemplary embodiment is shown, as follows: Figure 5 As shown, it includes:

[0183] The acquisition module 501 is used to acquire echo data of multiple pulses of pixels after wall filtering in the imaging region;

[0184] The first calculation module 502 is used to generate two sub-sequences based on the echo data, and to calculate the third-order cumulants of the two sub-sequences respectively.

[0185] The second calculation module 503 is used to calculate the cross-correlation function of the two subsequences using the third-order cumulants of the two subsequences;

[0186] The third calculation module 504 is used to calculate blood flow parameters based on the echo data and the cross-correlation function.

[0187] The blood flow parameter determination device provided in this application divides the echo data of the imaging region after wall filtering and multiple pulses into two subsequences, solves the third-order cumulant of each subsequence, and utilizes the characteristic that the third-order cumulant is insensitive to Gaussian noise to eliminate the problem that Gaussian noise affects the accuracy of blood flow parameter estimation by autocorrelation method, thereby improving the accuracy of blood flow parameter estimation. It can effectively ensure the stability of estimation in the blood vessel edge and low-velocity blood flow region, reduce the possibility of abrupt changes, and improve the coherence and hierarchy of blood flow imaging.

[0188] Based on the above embodiments, as a preferred embodiment, the first computing module 502 includes:

[0189] A construction unit is used to construct a first subsequence based on the echo data; wherein the first element of the first subsequence is the first element of the echo data, and the length of the first subsequence is the length of the echo data minus one.

[0190] A generation unit is used to perform a one-step delay on the echo data to generate a second subsequence; wherein the first element of the second subsequence is the second element of the echo data, and the length of the second subsequence is equal to that of the first subsequence.

[0191] Based on the above embodiments, as a preferred embodiment, the first computing module 502 includes:

[0192] A determining unit is used to determine the first delay parameter, the second delay parameter, and the third delay parameter;

[0193] The first delay unit is used to perform delay processing on the two sub-sequences based on the first delay parameter, the second delay parameter and the third delay parameter respectively, to obtain the first delay sequence, the second delay sequence and the third delay sequence corresponding to the two sub-sequences respectively;

[0194] The first calculation unit is used to perform conjugate operation on at least one delayed sequence corresponding to each of the two subsequences to obtain at least one conjugate operation result corresponding to each of the two subsequences, and to take the dot product result between the conjugate operation result corresponding to each of the two subsequences and other delayed sequences that have not undergone conjugate operation as the third-order cumulants of the two subsequences.

[0195] Based on the above embodiments, as a preferred implementation, the third delay parameter is zero. Accordingly, the first calculation unit is specifically used to: perform conjugate operation on the third delay sequences corresponding to the two subsequences respectively to obtain the conjugate operation results corresponding to the two subsequences respectively, and use the dot product result between the conjugate operation results corresponding to the two subsequences, the first delay sequence, and the second delay sequence as the third-order cumulants of the two subsequences respectively.

[0196] Accordingly, the second calculation module 503 is specifically used to: use the dot product between the conjugate of the third-order cumulant of the first subsequence and the third-order cumulant of the second subsequence as a cross-correlation function.

[0197] Based on the above embodiments, as a preferred implementation, the second delay parameter is equal to the third delay parameter. Accordingly, the first calculation unit is specifically used to: perform conjugate operations on the first delay sequence and the third delay sequence corresponding to the two sub-sequences respectively to obtain the first conjugate operation result and the second conjugate operation result corresponding to the two sub-sequences respectively, and use the dot product result between the first conjugate operation result, the second conjugate operation result, and the second delay sequence corresponding to the two sub-sequences respectively as the third-order cumulants of the two sub-sequences;

[0198] Accordingly, the second calculation module 503 is specifically used to: use the dot product between the conjugate of the third-order cumulant of the second subsequence and the third-order cumulant of the first subsequence as a cross-correlation function.

[0199] Based on the above embodiments, as a preferred embodiment, the first computing module 502 includes:

[0200] A determining unit is used to determine the first delay parameter, the second delay parameter, and the third delay parameter;

[0201] The second delay unit is used to perform delay processing on the first sub-sequence based on the first delay parameter, the second delay parameter and the third delay parameter respectively, to obtain the first delayed sequence, the second delayed sequence and the third delayed sequence;

[0202] The third delay unit is used to delay the second sub-sequence based on the first delay parameter to obtain a fourth delayed sequence;

[0203] The second calculation unit is used to perform a conjugate operation on the third delayed sequence, and use the dot product result between the conjugate of the third delayed sequence, the first delayed sequence, and the second delayed sequence as the third-order cumulant of the first subsequence.

[0204] The third calculation unit is used to take the dot product result between the conjugate of the third delay sequence, the fourth delay sequence, and the second delay sequence as the mutual third-order cumulant of the second subsequence;

[0205] Accordingly, the second calculation module 503 is specifically used to: take the dot product between the conjugate of the third-order cumulants of the first subsequence and the cross-third-order cumulants of the second subsequence as the cross-correlation function.

[0206] Based on the above embodiments, as a preferred embodiment, the blood flow parameters include any one or a combination of blood flow energy, blood flow velocity, and blood flow variance;

[0207] The third calculation module 504 is specifically used for: taking the echo data and the product of the conjugates of the echo data as an intermediate sequence, accumulating each element in the intermediate sequence to obtain the blood flow energy; and / or, accumulating each element in the cross-correlation function to obtain a first accumulated value, performing angle calculation on the first accumulated value to obtain a pulse frequency, and calculating the blood flow velocity according to the correspondence between the pulse frequency and the blood flow velocity; and / or, accumulating each element in the cross-correlation function to obtain a first accumulated value, accumulating the absolute value of each element in the cross-correlation function to obtain a second accumulated value, calculating the ratio of the absolute value of the first accumulated value to the second accumulated value, and using the difference between 1 and the ratio as the blood flow variance.

[0208] Regarding the apparatus in the above embodiments, the specific manner in which each module performs its operation has been described in detail in the embodiments related to the method, and will not be elaborated upon here.

[0209] The following describes a blood flow parameter display device provided in an embodiment of this application. This device can be referred to in conjunction with the blood flow parameter display method described above.

[0210] See Figure 6 A structural diagram of a blood flow parameter display device according to an exemplary embodiment is shown, as follows: Figure 6 As shown, it includes:

[0211] The acquisition module 601 is used to acquire echo data of multiple pulses of pixels after wall filtering in the imaging region;

[0212] The first calculation module 602 is used to generate two sub-sequences based on the echo data, and to calculate the third-order cumulants of the two sub-sequences respectively.

[0213] The second calculation module 603 is used to calculate the cross-correlation function of the two subsequences using the third-order cumulants of the two subsequences;

[0214] The third calculation module 604 is used to calculate blood flow parameters based on the echo data and the cross-correlation function;

[0215] Display module 605 is used to display the blood flow parameters.

[0216] The blood flow parameter display device provided in this application divides the echo data of the imaging region after wall filtering and multiple pulses into two subsequences, solves the third-order cumulant of each subsequence, and utilizes the characteristic that the third-order cumulant is insensitive to Gaussian colored noise to eliminate the problem that Gaussian colored noise affects the accuracy of blood flow parameter estimation by autocorrelation method, thereby improving the accuracy of blood flow parameter estimation and improving the display effect of blood flow.

[0217] Based on the hardware implementation of the above program modules, and in order to implement the method of the embodiments of this application, the embodiments of this application also provide an electronic device. Figure 7 This is a structural diagram of an electronic device according to an exemplary embodiment, such as... Figure 7 As shown, the electronic device includes:

[0218] Communication interface 1 enables information exchange with other devices, such as network devices;

[0219] Processor 2, connected to communication interface 1, enables information exchange with other devices and, when running a computer program, executes the blood flow parameter determination method provided by one or more of the aforementioned technical solutions. The computer program is stored in memory 3.

[0220] Of course, in practical applications, the various components in an electronic device are coupled together through bus system 4. It can be understood that bus system 4 is used to achieve communication and connection between these components. In addition to the data bus, bus system 4 also includes a power bus, a control bus, and a status signal bus. However, for clarity, in... Figure 7 The general will label all buses as Bus System 4.

[0221] The memory 3 in this embodiment is used to store various types of data to support the operation of the electronic device. Examples of such data include any computer program used to operate on the electronic device.

[0222] It is understood that memory 3 can be volatile memory or non-volatile memory, or both. Non-volatile memory can be read-only memory (ROM), programmable read-only memory (PROM), erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), ferromagnetic random access memory (FRAM), flash memory, magnetic surface memory, optical disc, or compact disc read-only memory (CD-ROM); magnetic surface memory can be disk storage or magnetic tape storage. Volatile memory can be random access memory (RAM), which is used as an external cache. By way of example, but not limitation, many forms of RAM are available, 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 (DDRSDRAM), Enhanced Synchronous Dynamic Random Access Memory (ESDRAM), SyncLink Dynamic Random Access Memory (SLDRAM), and Direct Rambus Random Access Memory (DRRAM).The memory 3 described in the embodiments of this application is intended to include, but is not limited to, these and any other suitable types of memory.

[0223] The methods disclosed in the embodiments of this application can be applied to processor 2, or implemented by processor 2. Processor 2 may be an integrated circuit chip with signal processing capabilities. In the implementation process, each step of the above method can be completed by the integrated logic circuit of the hardware in processor 2 or by instructions in the form of software. The processor 2 may be a general-purpose processor, DSP, or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. Processor 2 can implement or execute the methods, steps and logic block diagrams disclosed in the embodiments of this application. The general-purpose processor may be a microprocessor or any conventional processor, etc. The steps of the methods disclosed in the embodiments of this application can be directly manifested as being executed by a hardware decoding processor, or being executed by a combination of hardware and software modules in the decoding processor. The software modules may be located in a storage medium, which is located in memory 3. Processor 2 reads the program in memory 3 and completes the steps of the aforementioned method in combination with its hardware.

[0224] When processor 2 executes the program, it implements the corresponding processes in the various methods of the embodiments of this application. For the sake of brevity, these will not be described in detail here.

[0225] In an exemplary embodiment, this application also provides a storage medium, namely a computer storage medium, specifically a computer-readable storage medium, such as a memory 3 that stores a computer program, which can be executed by a processor 2 to complete the steps described in the aforementioned method. The computer-readable storage medium may be a memory such as FRAM, ROM, PROM, EPROM, EEPROM, Flash Memory, magnetic surface memory, optical disc, or CD-ROM.

[0226] Those skilled in the art will understand that all or part of the steps of the above method embodiments can be implemented by hardware related to program instructions. The aforementioned program can be stored in a computer-readable storage medium. When the program is executed, it performs the steps of the above method embodiments. The aforementioned storage medium includes various media that can store program code, such as mobile storage devices, ROM, RAM, magnetic disks, or optical disks.

[0227] Alternatively, if the integrated units described above are implemented as software functional modules and sold or used as independent products, they can also be stored in a computer-readable storage medium. Based on this understanding, the technical solutions of the embodiments of this application, or the parts that contribute to the prior art, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause an electronic device (which may be a personal computer, server, or network device, etc.) to execute all or part of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as mobile storage devices, ROM, RAM, magnetic disks, or optical disks.

[0228] The above description is merely a specific embodiment of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the scope of the technology disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A method for determining blood flow parameters, characterized in that, include: Acquire echo data of multiple pulses from pixels after wall filtering in the imaging region; Two subsequences are generated based on the echo data, and the third-order cumulants of the two subsequences are calculated respectively; The cross-correlation function of the two subsequences is calculated using the third-order cumulants of the two subsequences; Blood flow parameters are calculated based on the echo data and the cross-correlation function; The generation of two sub-sequences based on the echo data includes: The echo data is split into two subsequences of the same length but with different elements; Wherein, the blood flow parameters include blood flow velocity and / or blood flow variance; the calculation of blood flow parameters based on the echo data and the cross-correlation function includes: Each element in the cross-correlation function is accumulated to obtain a first accumulated value. An angle operation is performed on the first accumulated value to obtain the pulse frequency. The blood flow velocity is then calculated based on the correspondence between the pulse frequency and the blood flow velocity. And / or, summing up each element in the cross-correlation function to obtain a first accumulated value, summing up the absolute values ​​of each element in the cross-correlation function to obtain a second accumulated value, calculating the ratio of the absolute value of the first accumulated value to the second accumulated value, and taking the difference between 1 and the ratio as the blood flow variance.

2. The method for determining blood flow parameters according to claim 1, characterized in that, The generation of two sub-sequences based on the echo data includes: A first subsequence is constructed based on the echo data; wherein the first element of the first subsequence is the first element of the echo data, and the length of the first subsequence is the length of the echo data minus one. The echo data is delayed by one step to generate a second subsequence; wherein the first element of the second subsequence is the second element of the echo data, and the length of the second subsequence is equal to that of the first subsequence.

3. The method for determining blood flow parameters according to claim 1, characterized in that, The two subsequences are a first subsequence and a second subsequence, respectively. The calculation of the third-order cumulants for the two subsequences includes: Determine the first delay parameter, the second delay parameter, and the third delay parameter; Based on the first delay parameter, the second delay parameter, and the third delay parameter, the two sub-sequences are delayed respectively to obtain the first delayed sequence, the second delayed sequence, and the third delayed sequence corresponding to the two sub-sequences respectively; Perform conjugate operation on at least one delayed sequence corresponding to each of the two subsequences to obtain at least one conjugate operation result corresponding to each of the two subsequences. Use the dot product between the conjugate operation result corresponding to each of the two subsequences and other delayed sequences that have not undergone conjugate operation as the third-order cumulants of the two subsequences.

4. The method for determining blood flow parameters according to claim 3, characterized in that, The third delay parameter is zero. Correspondingly, the step of performing a conjugate operation on at least one delay sequence corresponding to each of the two sub-sequences to obtain at least one conjugate operation result corresponding to each of the two sub-sequences, and using the dot product between the conjugate operation result corresponding to each of the two sub-sequences and other delay sequences that have not undergone conjugate operation as the third-order cumulants of the two sub-sequences, includes: Perform conjugate operation on the third delay sequence corresponding to each of the two subsequences to obtain the conjugate operation results corresponding to each of the two subsequences. Use the conjugate operation results corresponding to each of the two subsequences, the dot product between the first delay sequence and the second delay sequence as the third-order cumulants of the two subsequences. Accordingly, the step of calculating the cross-correlation function of the two subsequences using the third-order cumulants of the two subsequences includes: The dot product between the conjugate of the third-order cumulants of the first subsequence and the third-order cumulants of the second subsequence is used as the cross-correlation function.

5. The method for determining blood flow parameters according to claim 3, characterized in that, The second delay parameter is equal to the third delay parameter. Correspondingly, the step of performing a conjugate operation on at least one delay sequence corresponding to each of the two subsequences to obtain at least one conjugate operation result corresponding to each of the two subsequences, and using the dot product between the conjugate operation result corresponding to each of the two subsequences and other delay sequences that have not undergone conjugate operation as the third-order cumulants of the two subsequences, includes: Perform conjugate operations on the first and third delayed sequences corresponding to the two subsequences to obtain the first and second conjugate operation results corresponding to the two subsequences respectively. Use the first conjugate operation result, the second conjugate operation result, and the dot product result between the second delayed sequences corresponding to the two subsequences as the third-order cumulants of the two subsequences respectively. Accordingly, the step of calculating the cross-correlation function of the two subsequences using the third-order cumulants of the two subsequences includes: The dot product between the conjugate of the third-order cumulant of the second subsequence and the third-order cumulant of the first subsequence is used as the cross-correlation function.

6. The method for determining blood flow parameters according to claim 1, characterized in that, The two subsequences are the first subsequence and the second subsequence, respectively. The calculation of the third-order cumulant of the subsequences includes: Determine the first delay parameter, the second delay parameter, and the third delay parameter; Based on the first delay parameter, the second delay parameter, and the third delay parameter, the first subsequence is delayed to obtain a first delayed sequence, a second delayed sequence, and a third delayed sequence. The second subsequence is delayed based on the first delay parameter to obtain a fourth delayed sequence; Perform a conjugate operation on the third delayed sequence, and use the dot product of the conjugate of the third delayed sequence, the first delayed sequence, and the second delayed sequence as the third-order cumulant of the first subsequence; The dot product of the conjugate of the third delay sequence, the fourth delay sequence, and the second delay sequence is used as the third-order cumulant of the second subsequence; Accordingly, the step of calculating the cross-correlation function of the two subsequences using the third-order cumulants of the two subsequences includes: The dot product between the conjugate of the third-order cumulants of the first subsequence and the cross-third-order cumulants of the second subsequence is used as the cross-correlation function.

7. The method for determining blood flow parameters according to claim 1, characterized in that, The blood flow parameters also include blood flow energy; the calculation of blood flow parameters based on the echo data and the cross-correlation function further 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.

8. A method for displaying blood flow parameters, characterized in that, include: Acquire echo data of multiple pulses from pixels after wall filtering in the imaging region; Two subsequences are generated based on the echo data, and the third-order cumulants of the two subsequences are calculated respectively; The cross-correlation function of the two subsequences is calculated using the third-order cumulants of the two subsequences; Blood flow parameters are calculated based on the echo data and the cross-correlation function; Display the blood flow parameters; The generation of two sub-sequences based on the echo data includes: The echo data is split into two subsequences of the same length but with different elements; Wherein, the blood flow parameters include blood flow velocity and / or blood flow variance; the calculation of blood flow parameters based on the echo data and the cross-correlation function includes: Each element in the cross-correlation function is accumulated to obtain a first accumulated value. An angle operation is performed on the first accumulated value to obtain the pulse frequency. The blood flow velocity is then calculated based on the correspondence between the pulse frequency and the blood flow velocity. And / or, summing up each element in the cross-correlation function to obtain a first accumulated value, summing up the absolute values ​​of each element in the cross-correlation function to obtain a second accumulated value, calculating the ratio of the absolute value of the first accumulated value to the second accumulated value, and taking the difference between 1 and the ratio as the blood flow variance.

9. A device for determining blood flow parameters, characterized in that, include: The acquisition module is used to acquire echo data of multiple pulses from pixels after wall filtering in the imaging region; The first calculation module is used to generate two sub-sequences based on the echo data, and to calculate the third-order cumulants of the two sub-sequences respectively. The second calculation module is used to calculate the cross-correlation function of the two subsequences using the third-order cumulants of the two subsequences; The third calculation module is used to calculate blood flow parameters based on the echo data and the cross-correlation function; Specifically, the first calculation module is used to: split the echo data into two subsequences of the same length but with different elements; The blood flow parameters include blood flow velocity and / or blood flow variance. The third calculation module is specifically used for: accumulating each element in the cross-correlation function to obtain a first accumulated value; performing angle calculations on the first accumulated value to obtain a pulse frequency; calculating the blood flow velocity based on the correspondence between the pulse frequency and the blood flow velocity; and / or, accumulating each element in the cross-correlation function to obtain a first accumulated value; accumulating the absolute values ​​of each element in the cross-correlation function to obtain a second accumulated value; calculating the ratio of the absolute value of the first accumulated value to the second accumulated value; and using the difference between 1 and the ratio as the blood flow variance.

10. A blood flow parameter display device, characterized in that, include: The acquisition module is used to acquire echo data of multiple pulses from pixels after wall filtering in the imaging region; The first calculation module is used to generate two sub-sequences based on the echo data, and to calculate the third-order cumulants of the two sub-sequences respectively. The second calculation module is used to calculate the cross-correlation function of the two subsequences using the third-order cumulants of the two subsequences; The third calculation module is used to calculate blood flow parameters based on the echo data and the cross-correlation function; The display module is used to display the blood flow parameters; Specifically, the first calculation module is used to: split the echo data into two subsequences of the same length but with different elements; The blood flow parameters include blood flow velocity and / or blood flow variance. The third calculation module is specifically used for: accumulating each element in the cross-correlation function to obtain a first accumulated value; performing angle calculations on the first accumulated value to obtain a pulse frequency; calculating the blood flow velocity based on the correspondence between the pulse frequency and the blood flow velocity; and / or, accumulating each element in the cross-correlation function to obtain a first accumulated value; accumulating the absolute values ​​of each element in the cross-correlation function to obtain a second accumulated value; calculating the ratio of the absolute value of the first accumulated value to the second accumulated value; and using the difference between 1 and the ratio as the blood flow variance.

11. An electronic device, characterized in that, include: Memory, used to store computer programs; A processor for executing the computer program to implement the steps of the method as claimed in any one of claims 1 to 8.

12. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, implements the steps of the method as described in any one of claims 1 to 8.

Images