Pulse Doppler filling data determination method and apparatus, and electronic device

Abstract

The embodiment of the invention provides a pulse Doppler filling data determination method and device and electronic equipment. The method comprises the following steps: acquiring reference Doppler data before and / or after a to-be-filled time period, wherein the reference Doppler data is obtained through scanning of Doppler ultrasonic equipment; performing spectrum analysis on the reference Doppler data to obtain corresponding reference spectrum information; and predicting Doppler data to be filled in the time period to be filled based on the reference frequency spectrum information to obtain target Doppler data, target frequency spectrum information corresponding to the target Doppler data being determined based on the reference frequency spectrum information. The target Doppler data obtained through prediction by the scheme is relatively accurate, and the scheme is relatively small in calculation amount and can be suitable for an application scene of real-time scanning.

Description

Technical Field

[0001] This invention relates to the field of ultrasound imaging technology, and more specifically, to a method and apparatus for determining pulsed Doppler filling data, an electronic device, and a storage medium. Background Technology

[0002] Medical imaging technology is a widely used non-invasive or minimally invasive diagnostic technique that can image human tissues and organs through various imaging technologies to obtain ultrasound images of human tissues, assisting users in measuring and diagnosing lesions in human tissues and organs.

[0003] In current clinical ultrasound applications, ultrasound imaging modes can include various modes such as two-dimensional ultrasound imaging (B-mode), color Doppler ultrasound imaging (C-mode), and pulsed Doppler imaging (D-mode). D-mode can be combined with either B-mode or C-mode to form a multiplexed ultrasound imaging system. For example, by combining D-mode with B-mode, during ultrasound scanning, the ultrasound device's display can simultaneously show anatomical images from B-mode and pulsed Doppler images from D-mode.

[0004] Because this technology requires the simultaneous display of images from multiple modes, compared to the simplex mode in D imaging, pulsed Doppler multiplexing requires periodically selecting a period during scanning to perform B-mode or C-mode scans in order to provide the signal data for the corresponding display mode. However, since D imaging requires a continuous one-dimensional signal as input, the Doppler data for the time period using B-mode or C-mode scans needs to be filled in by the ultrasound system to ensure the temporal continuity of the Doppler data.

[0005] In related technologies, data reuse is typically used to fill in the Doppler data for this period. This involves directly reusing Doppler data from a period preceding the time to be filled, followed by signal processing methods to ensure smooth edges. However, this method is demanding, requiring users to continuously debug to find the optimal filling method, and it is also prone to causing a certain degree of spectral leakage. Summary of the Invention

[0006] The present invention was proposed in view of the above-mentioned problems. The present invention provides a method for determining pulse Doppler filling data, a device for determining pulse Doppler filling data, an electronic device, and a storage medium.

[0007] According to one aspect of the present invention, a method for determining pulsed Doppler filling data is provided. The method includes: acquiring reference Doppler data located before and / or after the filling period, wherein the reference Doppler data is obtained by scanning with a Doppler ultrasound device; performing spectral analysis on the reference Doppler data to obtain corresponding reference spectral information; and predicting the Doppler data to be filled within the filling period based on the reference spectral information to obtain target Doppler data, wherein the target spectral information corresponding to the target Doppler data is determined based on the reference spectral information.

[0008] For example, the reference Doppler data includes a sequence of reference data points before the period to be filled and / or a sequence of reference data points after the period to be filled, each sequence of reference data points including one or more Doppler data points; the reference spectral information corresponding to each reference data point sequence includes the spectral components corresponding to the reference data point sequence and the intensity coefficient of each spectral component.

[0009] For example, performing spectral analysis on reference Doppler data to obtain corresponding reference spectral information includes: for each reference data point sequence, performing the following specific analysis operation: dividing the frequency range corresponding to each Doppler data point in the reference data point sequence into a first number of segments to obtain a first number of frequency components; determining a first number of intensity coefficients corresponding one-to-one with the first number of frequency components by minimizing the difference between the signal intensity of each first Doppler data point in the fitted data sequence and the signal intensity of the corresponding second Doppler data point in the reference data point sequence; wherein, the fitted data sequence includes first Doppler data points corresponding one-to-one with the second Doppler data points in the reference data point sequence, and the signal intensity of the first Doppler data points is determined based on the first number of frequency components, the first number of intensity coefficients, and the sampling time of the corresponding second Doppler data points; the first number is an integer greater than or equal to 1.

[0010] For example, the signal strength of each Doppler data point in the reference data point sequence includes a real part and an imaginary part, and a specific analysis operation is performed on the real part and the imaginary part respectively; when performing a specific analysis operation on the real part, the signal strength of the first Doppler data point is the real part of the first Doppler data point, the signal strength of the second Doppler data point is the real part of the second Doppler data point, a first number of frequency components are the first number of frequency components corresponding to the real part, and a first number of intensity coefficients are the first number of intensity coefficients corresponding to the real part; when performing a specific analysis operation on the imaginary part, the signal strength of the first Doppler data point is the imaginary part of the first Doppler data point, the signal strength of the second Doppler data point is the imaginary part of the second Doppler data point, a first number of frequency components are the first number of frequency components corresponding to the imaginary part, and a first number of intensity coefficients are the first number of intensity coefficients corresponding to the imaginary part.

[0011] For example, predicting the Doppler data to be filled within the filling period based on reference spectrum information to obtain target Doppler data includes: predicting the Doppler data to be filled within the filling period based on reference spectrum information and a preset time weighting coefficient to obtain target Doppler data; wherein, the signal strength of the target Doppler data point in the target Doppler data is determined based on the target spectrum information and the sampling time corresponding to the target Doppler data point, and the target spectrum information is obtained by adjusting the reference spectrum information with a preset time weighting coefficient, the time weighting coefficient being related to the sampling time corresponding to the target Doppler data point.

[0012] For example, the reference spectrum information includes spectral components and intensity coefficients for each spectral component; the target Doppler data includes a second number of target Doppler data points, where the second number is an integer greater than or equal to 1; the preset time weighting coefficients include time weighting coefficients corresponding one-to-one with the second number of target Doppler data points; the target spectrum information includes intensity coefficients corresponding one-to-one with the second number of target Doppler data points; and the spectral components in the target spectrum information are consistent with the spectral components in the reference spectrum information. The intensity coefficients corresponding to each target Doppler data point in the target spectrum information are determined based on the product of the corresponding time weighting coefficient and the intensity coefficients in the reference spectrum information. For the reference spectrum information corresponding to a sequence of reference data points before the time period to be filled, the second number of target Doppler data points are arranged in chronological order of sampling time, and the corresponding time weighting coefficients gradually decrease. For the reference spectrum information corresponding to a sequence of reference data points after the time period to be filled, the second number of target Doppler data points are arranged in chronological order of sampling time, and the corresponding time weighting coefficients gradually increase.

[0013] For example, for the reference spectral information corresponding to the reference data point sequence located before the time period to be filled, the time weighting coefficients corresponding to the second number of target Doppler data points are determined based on the following formula:

[0014]

[0015] Where W(s) represents the time weighting coefficient, s represents the s-th target Doppler data point, s∈[1,L], and L represents the second quantity; for the reference spectrum information corresponding to the reference data point sequence located after the time period to be filled, the time weighting coefficient corresponding to the s-th target Doppler data point is equal to 1-W(s).

[0016] For example, the reference Doppler data includes a sequence of reference data points before the period to be filled and a sequence of reference data points after the period to be filled, each sequence of reference data points including one or more Doppler data points; the preset time weighting coefficient includes the time weighting coefficient corresponding to each reference data point sequence; predicting the Doppler data to be filled within the period to be filled based on the reference spectrum information and the preset time weighting coefficient to obtain the target Doppler data includes: for each target Doppler data point in the target Doppler data, for each reference data point sequence, determining the predicted signal intensity corresponding to the reference data point sequence based on the time weighting coefficient corresponding to the reference data point sequence and the reference spectrum information corresponding to the reference data point sequence, and summing the predicted signal intensities corresponding to all reference data point sequences to obtain the signal intensity of the target Doppler data point.

[0017] For example, based on reference spectrum information and preset time weighting coefficients, the Doppler data to be filled within the time period is predicted to obtain the target Doppler data, including: determining the signal strength of each target Doppler data point in the target Doppler data using the following formula:

[0018]

[0019] Among them, g1, g2, g3, ..., g L Let I represent the signal intensity corresponding to the 1st, 2nd, 3rd, ..., Lth target Doppler data points, respectively; let Q represent the real part of each target Doppler data point; let Q represent the imaginary part of each target Doppler data point; and let K represent the intensity coefficient corresponding to the real part of the reference Doppler data point in the reference data point sequence preceding the time period to be filled. K represents the intensity coefficient corresponding to the real part of the reference Doppler data point in the reference data point sequence following the time period to be filled. Q This represents the intensity coefficient corresponding to the imaginary part of the reference Doppler data points in the reference data point sequence preceding the time period to be filled. W1 represents the intensity coefficient corresponding to the imaginary part of the reference Doppler data points in the reference data point sequence after the time period to be filled, and W2 represents the time weight coefficients corresponding to the reference data point sequences before and after the time period to be filled, respectively.

[0020] According to another aspect of the present invention, a pulsed Doppler filling data determination device is also provided. The device includes: an acquisition module for acquiring reference Doppler data located before and / or after the filling time period, wherein the reference Doppler data is obtained by scanning with a Doppler ultrasound device; an analysis module for performing spectral analysis on the reference Doppler data to obtain corresponding reference spectral information; and a prediction module for predicting the Doppler data to be filled within the filling time period based on the reference spectral information to obtain target Doppler data, wherein the target spectral information corresponding to the target Doppler data is determined based on the reference spectral information.

[0021] According to another aspect of the present invention, an electronic device is also provided, including a processor and a memory, wherein a computer program is stored in the memory, and the computer program instructions are executed by the processor to perform the above-described pulse Doppler filling data determination method.

[0022] According to another aspect of the present invention, a storage medium is also provided, which stores a computer program / instructions, which, when executed, are used to perform the above-described pulse Doppler filling data determination method.

[0023] According to embodiments of the present invention, the pulse Doppler filling data determination method, device, electronic device, and storage medium can obtain corresponding reference spectrum information by performing spectral analysis on reference Doppler data located before and / or after the filling time period. Based on the reference spectrum information, the target Doppler data can be obtained by predicting the Doppler data to be filled within the filling time period. During the prediction process, the target spectrum information corresponding to the target Doppler data is determined based on the reference spectrum information. That is, the spectrum situation within the filling time period can be predicted using reference spectrum information before and / or after the filling time period to determine the Doppler data within that period. This solution has low implementation requirements and does not require users to determine the target Doppler data through repeated debugging. Furthermore, this solution obtains the target Doppler data based on reference spectrum information, making the target Doppler data more consistent with the continuity of spectral changes, thus avoiding spectral leakage to a certain extent. Therefore, this solution can more naturally fill in the missing Doppler data during scanning in multiplex ultrasound imaging mode, ensuring the accuracy of the final displayed pulse Doppler image. In addition, this method has low computational complexity and is well-suited for real-time scanning applications.

[0024] The above description is merely an overview of the technical solution of the present invention. In order to better understand the technical means of the present invention and to implement it in accordance with the contents of the specification, and in order to make the above and other objects, features and advantages of the present invention more apparent and understandable, specific embodiments of the present invention are described below. Attached Figure Description

[0025] The above and other objects, features, and advantages of the present invention will become more apparent from the more detailed description of the embodiments of the invention in conjunction with the accompanying drawings. The drawings are provided to further illustrate the embodiments of the invention and form part of the specification. They are used together with the embodiments of the invention to explain the invention and do not constitute a limitation thereof. In the drawings, the same reference numerals generally represent the same parts or steps.

[0026] Figure 1 A schematic flowchart of a pulse Doppler filling data determination method according to an embodiment of the present invention is shown;

[0027] Figure 2 A schematic diagram of Doppler data according to an embodiment of the present invention is shown;

[0028] Figure 3 A schematic block diagram of a pulse Doppler filling data determination device according to an embodiment of the present invention is shown; and

[0029] Figure 4 A schematic block diagram of an electronic device according to an embodiment of the present invention is shown. Detailed Implementation

[0030] To make the objectives, technical solutions, and advantages of the present invention more apparent, exemplary embodiments according to the present invention will be described in detail below with reference to the accompanying drawings. Obviously, the described embodiments are merely a part of the embodiments of the present invention, and not all of the embodiments of the present invention. It should be understood that the present invention is not limited to the exemplary embodiments described herein. Based on the embodiments described in the present invention, all other embodiments obtained by those skilled in the art without inventive effort should fall within the protection scope of the present invention.

[0031] As mentioned above, in related technologies, data reuse is typically used to fill in the Doppler data within the time period to be filled. This involves directly reusing Doppler data from a period preceding the time period to be filled, supplemented by signal processing methods to ensure smooth edges. This method only involves calculating the signal strength of the Doppler data; the signal strength of the Doppler data from a period preceding the time period to be filled is determined directly or after certain adjustments as the signal strength of the Doppler data within the time period to be filled. Since directly setting the signal strength is not robust, users need to continuously debug to find the optimal filling method. Furthermore, this method of directly setting the signal strength is also prone to causing a certain degree of spectral leakage.

[0032] To address the aforementioned technical problems, this invention provides a method for determining pulse Doppler filling data. Figure 1 A schematic flowchart of a pulse Doppler filling data determination method 100 according to an embodiment of the present invention is shown, as follows: Figure 1 As shown, the method 100 may include the following steps S110, S120 and S130.

[0033] Step S110: Obtain reference Doppler data before and / or after the time period to be filled. The reference Doppler data is obtained by scanning with a Doppler ultrasound device.

[0034] For example, a Doppler ultrasound device can be used to scan a target object to obtain Doppler data (i.e., pulsed Doppler data) corresponding to the target object. Those skilled in the art will understand that the Doppler data can be pulsed Doppler blood flow signals. The target object can be any organism or any part of an organism to be tested. The part to be tested can be a biological tissue, such as the coronary arteries or veins of the heart. D-mode imaging can be primarily used to examine moving organs and flowing fluids, such as the heart, blood vessels and the blood flowing within them, or fetal heartbeats, to understand their motion state or measure blood flow velocity and direction. When using D-mode imaging to scan the part to be tested, it can also be combined with B-mode imaging and / or C-mode imaging to form a multiplexed imaging mode. Using multiplexed imaging mode, both the morphology of the part to be tested and the direction and velocity of blood flow can be observed. For example, the T1 duration of the part to be tested can be scanned using D-mode imaging, and then the T2 duration can be scanned using B-mode imaging. After scanning the area under test using the B imaging mode for duration T2, the scan is then performed again using the D imaging mode for duration T3. During the B imaging mode scan, the Doppler data corresponding to the D imaging mode is missing. To ensure the continuity of the Doppler data, the Doppler data within duration T2 needs to be filled. The time period corresponding to duration T2 can represent the time period to be filled. It should be noted that although the above description uses the example of Doppler data loss due to imaging mode switching to illustrate the technical problems in the related art, it should be understood that the pulsed Doppler filling data determination method according to the embodiments of the present invention can be applied to any scenario where Doppler data is missing. For example, if pulsed Doppler imaging is paused during the time period to be filled due to Doppler ultrasound equipment malfunction or other reasons, as long as Doppler data was acquired before and / or after the time period to be filled, the pulsed Doppler filling data determination method according to the embodiments of the present invention can be used to predict the Doppler data within the time period to be filled, in order to fill the missing portion.

[0035] Figure 2 A schematic diagram of Doppler data according to an embodiment of the present invention is shown. Figure 2As shown, time period 210 is the time period to be filled, the period before 210 is the preceding period 220, and the period after 210 is the following period 230. During the time period to be filled 210, the Doppler ultrasound equipment is not operating in D-imaging mode; for example, it may have switched to an imaging mode other than D-imaging mode, or the Doppler ultrasound equipment may have paused pulse Doppler imaging due to malfunction or other reasons. Therefore, the Doppler data within the time period to be filled 210 needs to be filled.

[0036] The reference Doppler data may include first reference Doppler data preceding the period to be filled and / or second reference Doppler data following the period to be filled. The first reference Doppler data may be all or part of the Doppler data collected before the period to be filled. For example, regardless of how much Doppler data was collected before the period to be filled, all Doppler data collected before the period to be filled can be used as the first reference Doppler data. Alternatively, Doppler data within a first preset period can be selected from the Doppler data collected before the period to be filled as the first reference Doppler data. The duration of the first preset period can be set to any suitable duration as needed, and this invention does not limit this. For example, the duration of the first preset period may be in the range of [60, 80] milliseconds (ms), such as 65ms, 70ms, 75ms, etc. Preferably, the first preset period is adjacent to the period to be filled. The second reference Doppler data may be all or part of the Doppler data collected after the period to be filled. For example, regardless of how much Doppler data is collected after the period to be filled, all Doppler data collected after the period to be filled can be used as the second reference Doppler data. Alternatively, Doppler data within a second preset period can be selected from the Doppler data collected after the period to be filled as the second reference Doppler data. The duration of the second preset period can be set to any suitable duration as needed, and this invention does not limit this. For example, the duration of the second preset period can be in the range of [60, 80] ms, such as 65 ms, 70 ms, 75 ms, etc. Preferably, the second preset period is adjacent to the period to be filled.

[0037] Step S120: Perform spectral analysis on the reference Doppler data to obtain the corresponding reference spectral information.

[0038] For example, the reference Doppler data can be processed by performing Fourier transforms or other methods to conduct spectral analysis. Fourier transforms include, but are not limited to, Discrete-time Fourier Transform (DTFT) and Fast Fourier Transform (FFT). By performing spectral analysis on the reference Doppler data, corresponding reference spectral information can be obtained. For example, the reference spectral information may include spectral components (also called spectral elements) and the intensity coefficient of each spectral component. The spectral components can represent the amplitude of the corresponding frequency in the frequency domain for each Doppler data point in the reference Doppler data after performing a Fourier transform.

[0039] Step S130: Based on the reference spectrum information, predict the Doppler data to be filled in the time period to obtain the target Doppler data, wherein the target spectrum information corresponding to the target Doppler data is determined based on the reference spectrum information.

[0040] For example, based on the obtained reference spectral information, the Doppler data to be filled within the time period can be predicted. The prediction result can be used as the target Doppler data. It is understood that the Doppler data to be filled within the time period is not actual data; the Doppler data within the time period is actually blank before filling. After predicting the Doppler data to be filled within the time period, the predicted target Doppler data can be used as the Doppler data within the time period. For example, the target spectral information may include the spectral components and intensity coefficients corresponding to each target Doppler data point in the target Doppler data.

[0041] According to the above technical solution, by performing spectral analysis on reference Doppler data located before and / or after the time period to be filled, corresponding reference spectral information can be obtained. Based on the reference spectral information, the Doppler data to be filled within the time period can be predicted to obtain the target Doppler data. During the prediction process, the target spectral information corresponding to the target Doppler data is determined based on the reference spectral information. In other words, the spectral situation within the time period to be filled can be predicted using reference spectral information before and / or after the time period to be filled, thereby determining the Doppler data within that time period. This solution has low implementation requirements, eliminating the need for users to repeatedly debug to determine the target Doppler data. Furthermore, this solution predicts the target Doppler data based on reference spectral information, ensuring that the target Doppler data better conforms to the continuity of spectral changes, thus avoiding spectral leakage to some extent. Therefore, this solution can more naturally fill in the missing Doppler data during scanning in multiplex ultrasound imaging mode, ensuring the accuracy of the final displayed pulsed Doppler image. In addition, this method has low computational complexity and is well-suited for real-time scanning applications.

[0042] For example, any type of prediction algorithm, such as averaging, exponential smoothing, or linear regression, can be used to predict the target Doppler data. Furthermore, any type of adaptive filter, such as Wiener or Kalman filters, can be used to analyze the reference Doppler data to obtain the target Doppler data.

[0043] For example, the reference Doppler data may include a sequence of reference data points before the period to be filled and / or a sequence of reference data points after the period to be filled, each sequence of reference data points including one or more Doppler data points; the reference spectral information corresponding to each reference data point sequence includes the spectral components corresponding to the reference data point sequence and the intensity coefficient of each spectral component.

[0044] In one embodiment, the reference Doppler data may include a sequence of reference data points preceding the period to be filled (which may be referred to as the first reference data point sequence) and / or a sequence of reference data points following the period to be filled (which may be referred to as the second reference data point sequence). For ease of description and understanding, the first reference data point sequence preceding the period to be filled will be described below as an example. The first reference data point sequence includes one or more Doppler data points. (Referring again...) Figure 2The first reference data point sequence, located before the time period to be filled, may include eight Doppler data points I1, I2, ..., I8. The first reference data point sequence can be represented as [I1, I2, ..., I8]. The reference spectrum information corresponding to the first reference data point sequence may include the spectral components corresponding to that reference data point sequence and the intensity coefficient of each spectral component. The number of Doppler data points included in the first reference data point sequence can be set as needed, and can be any integer greater than or equal to 1. The second reference data point sequence is similar to the first reference data point sequence and will not be described further. The number of Doppler data points included in the first and second reference data point sequences may be the same or different.

[0045] According to the above technical solution, the reference Doppler data may include a sequence of reference data points preceding the period to be filled and / or a sequence of reference data points following the period to be filled. The reference spectral information corresponding to each reference data point sequence includes the spectral components corresponding to that reference data point sequence and the intensity coefficient of each spectral component. This method can determine the corresponding spectral components and intensity coefficients for reference data points preceding and following the period to be filled, thus improving the accuracy of the determined reference spectral information.

[0046] For example, performing spectral analysis on reference Doppler data to obtain corresponding reference spectral information may include: for each reference data point sequence, performing the following specific analysis operation: dividing the frequency range corresponding to each Doppler data point in the reference data point sequence into a first number of segments to obtain a first number of frequency components; determining a first number of intensity coefficients corresponding one-to-one with the first number of frequency components by minimizing the difference between the signal intensity of each first Doppler data point in the fitted data sequence and the signal intensity of the corresponding second Doppler data point in the reference data point sequence; wherein the fitted data sequence includes first Doppler data points corresponding one-to-one with the second Doppler data points in the reference data point sequence, and the signal intensity of the first Doppler data points is determined based on the first number of frequency components, the first number of intensity coefficients, and the sampling time of the corresponding second Doppler data points; the first number is an integer greater than or equal to 1.

[0047] In one embodiment, specific analysis operations can be performed on each reference data point sequence to obtain the reference spectrum information corresponding to that reference data point sequence. The following explanation still uses the first reference data point sequence as an example. The sampling time sequence T of the Doppler data points in the first reference data point sequence can be expressed as (without loss of generality, assuming that data sampling starts from time 0): T = [t1, t2, ..., t...]. N1 ] T =[1 / f s ,2 / fs , ..., N / f s ] T N1 represents the total number of Doppler data points in the first reference data point sequence, where N1 is any integer greater than or equal to 1. According to the Nyquist sampling theorem, the frequency range corresponding to each Doppler data point in the first reference data point sequence can be considered to be [-f...]. s / 2,f s / 2]. This frequency range can be uniformly divided into a first number (M) segments to obtain M frequency components. M is any integer greater than or equal to 1, and this invention does not limit this. Therefore, the spectral component sequence of each Doppler data point in the first reference data point sequence can be represented as [cos(2πf c1 ), cos(2πf c2 ), ..., cos(2πf cM )]. Wherein, cos(2πf c1 ), cos(2πf c2 Each of these can represent a spectral component. The frequency range corresponding to each spectral component is f. cm =-f s / 2+mf s / M, where m∈[1,M]. The spectral components corresponding to the second reference data point sequence may be consistent with those corresponding to the first reference data point sequence, which will not be elaborated further. The intensity coefficients corresponding to the second reference data point sequence may or may not be consistent with those corresponding to the first reference data point sequence, depending on the specific spectral analysis results. The spectral analysis process can be transformed into a data fitting process. The implementation method of data fitting is described below.

[0048] In one embodiment, a fitted data sequence can be constructed, including Doppler data points (referred to as first Doppler data points) that correspond one-to-one with Doppler data points (referred to as second Doppler data points) in the current reference data point sequence. The signal intensity of the first Doppler data points can be determined based on M frequency components, M intensity coefficients, and the sampling time of the corresponding second Doppler data points. The values ​​of the M intensity coefficients can be determined by minimizing the difference between the signal intensity of each first Doppler data point in the fitted data sequence and the signal intensity of the corresponding second Doppler data point in the reference data point sequence. The difference between the signal intensity of the nth first Doppler data point in the fitted data sequence and the signal intensity of the corresponding second Doppler data point can be represented by ω. n Let n ∈ [1, N], where N represents the number of Doppler data points included in the current reference data point sequence. The difference ω can be expressed as... n The strength coefficient k corresponding to the minimum value mThese M strength coefficients can be considered as a sequence K. The sequence K can be expressed as: K = [k1, k2, ... k M ] T For example, ω n It can be calculated using the following formula:

[0049]

[0050] Where n represents the nth first Doppler data point, n∈[1,N]. For the first reference data point sequence, N=N1; for the second reference data point sequence, N=N2. N1 is the number of Doppler data points in the first reference data point sequence, and N2 is the number of Doppler data points in the second reference data point sequence. For ease of description, a spectral analysis matrix C can be defined:

[0051]

[0052] For example, the intensity coefficient sequence K can be expressed as: K = (C T C) -1 C T I.

[0053] For both the first and second reference data point sequences, the corresponding intensity coefficient sequences can be determined using the methods described above.

[0054] According to the above technical solution, for each reference data point sequence, the frequency range corresponding to each Doppler data point in the reference data point sequence can be divided into a first number of segments to obtain a first number of frequency components. Then, by minimizing the difference between the signal intensity of each first Doppler data point in the fitted data sequence and the signal intensity of the corresponding second Doppler data point in the reference data point sequence, a first number of intensity coefficients corresponding one-to-one with the first number of frequency components are determined. The intensity coefficients obtained by this method can minimize the difference between the signal intensity of each first Doppler data point in the fitted data sequence and the signal intensity of the corresponding second Doppler data point in the reference data point sequence, so as to ensure that the signal intensity of the fitted first Doppler data point is closer to the signal intensity of the corresponding second Doppler data point, thereby ensuring the accuracy of the target Doppler data determined based on such intensity coefficients.

[0055] For example, the signal strength of each Doppler data point in the reference data point sequence may include a real part and an imaginary part, and a specific analysis operation is performed on the real part and the imaginary part respectively; when performing a specific analysis operation on the real part, the signal strength of the first Doppler data point is the real part of the first Doppler data point, the signal strength of the second Doppler data point is the real part of the second Doppler data point, a first number of frequency components are the first number of frequency components corresponding to the real part, and a first number of intensity coefficients are the first number of intensity coefficients corresponding to the real part; when performing a specific analysis operation on the imaginary part, the signal strength of the first Doppler data point is the imaginary part of the first Doppler data point, the signal strength of the second Doppler data point is the imaginary part of the second Doppler data point, a first number of frequency components are the first number of frequency components corresponding to the imaginary part, and a first number of intensity coefficients are the first number of intensity coefficients corresponding to the imaginary part.

[0056] In one embodiment, the Doppler data points can be complex points. The sequence of first reference data points preceding the time period to be filled can be represented as Y = [I1 + jQ1, I2 + jQ2, ..., I...]. N +jQ N1 ] T =I+jQ. The sequence of second reference data points following the period to be filled can be represented as: The following explanation uses the first reference data point sequence as an example. The first reference data point sequence I = [I1, I2, ..., I8] can represent the sequence corresponding to the real part of each Doppler data point in the first reference data point sequence. The sequence Q corresponding to the imaginary part of each Doppler data point in the first reference data point sequence can be represented as Q = [Q1, Q2, ..., Q8]. In this embodiment, N1 equals 8. Specific analysis operations are performed on the real and imaginary parts of the Doppler data points respectively. The implementation method of performing specific analysis operations on the real part of the first Doppler data points has been described in detail in the previous embodiment, and will not be repeated here for brevity. The intensity coefficient sequence corresponding to the real part of the Doppler data points in the first reference data point sequence is K = (C T C) -1 C T I. In a similar manner, specific analysis operations can be performed on the imaginary part of the Doppler data points in the first reference data point sequence to obtain the intensity coefficient sequence K corresponding to the imaginary part of the Doppler data points in the first reference data point sequence. Q =(C Q T C Q ) -1 C Q T Q. Wherein, the spectral analysis matrix C corresponds to the imaginary part of the Doppler data points in the first reference data point sequence.Q This can be expressed by the following formula:

[0057]

[0058] The above embodiments mainly focus on performing spectral analysis on the real and imaginary parts of each Doppler data point in the first reference data point sequence before the filling period. Similarly, spectral analysis can be performed on the real and imaginary parts of each Doppler data point in the second reference data point sequence after the filling period to obtain the intensity coefficients corresponding to the real and imaginary parts of each Doppler data point in the second reference data point sequence. and In the above embodiments, the real part of the Doppler data points is represented by a cosine function, and the imaginary part is represented by a sine function. This representation is merely exemplary and does not limit the invention. For example, the real part of the Doppler data points can also be represented by a sine function, and the imaginary part by a cosine function.

[0059] According to the above technical solution, specific analysis operations can be performed on the real and imaginary parts of the Doppler data points respectively, so as to obtain reference spectrum information more comprehensively.

[0060] For example, predicting the Doppler data to be filled within the filling period based on reference spectrum information to obtain target Doppler data may include: predicting the Doppler data to be filled within the filling period based on reference spectrum information and a preset time weighting coefficient to obtain target Doppler data; wherein, the signal strength of the target Doppler data point in the target Doppler data is determined based on the target spectrum information and the sampling time corresponding to the target Doppler data point, and the target spectrum information is determined by adjusting the reference spectrum information with a preset time weighting coefficient, the time weighting coefficient being related to the sampling time corresponding to the target Doppler data point.

[0061] In one embodiment, the target Doppler data can be obtained by predicting the Doppler data to be filled within the filling period based on reference spectrum information and a preset time weighting coefficient. The preset time weighting coefficient can be implemented using any weighting function, such as a simple moving average, an exponential moving average, or a weighted moving average, etc., and this invention is not limited thereto. Based on the reference spectrum information corresponding to the first reference data point sequence and / or the second reference data point sequence and the preset time weighting coefficient, the spectrum information corresponding to the Doppler data to be filled within the filling period (referred to herein as target spectrum information) can be determined. Based on the target spectrum information and the sampling time corresponding to the target Doppler data points, the signal strength of the target Doppler data points in the target Doppler data can be determined. The Doppler data to be filled within the filling period may include one or more Doppler data points, which, after determining their corresponding signal strengths, form the target Doppler data points described herein. The number of target Doppler data points can be set as needed, and can be any integer greater than or equal to 1. The sampling time series T corresponding to the target Doppler data points g It can be represented as: T g = [t1, t2, ... t L ] T =[1 / f s ,2 / f s , ..., L / f s ] T L represents the number of target Doppler data points, which can be any integer greater than or equal to 1. The target spectral information can be determined based on reference spectral information and a preset time weighting coefficient. For example, the target spectral information can be obtained by adjusting the size of the reference spectral information using the preset time weighting coefficient. For instance, the target spectral information can be determined by multiplying the reference spectral information and the preset time weighting coefficient. The time weighting coefficient is related to the sampling time corresponding to the target Doppler data points, thus associating the target spectral information corresponding to the target Doppler data points with the sampling time corresponding to the target Doppler data points. This allows target Doppler data points with different sampling times to have different target spectral information, resulting in target Doppler data points with more temporally smooth and continuous target spectral information.

[0062] According to the above technical solution, the Doppler data to be filled within the time period is predicted based on reference spectral information and a preset time weighting coefficient to obtain the target Doppler data. This method can incorporate a preset time weighting coefficient when obtaining the target spectral information, thus adjusting the reference spectral information using the time weighting coefficient to obtain the target spectral information. This helps to make the obtained target spectral information within the time period to be filled more smooth and continuous in time, resulting in higher accuracy.

[0063] For example, the target Doppler data may include a second number of target Doppler data points, where the second number is an integer greater than or equal to 1. Preset time weighting coefficients include time weighting coefficients corresponding one-to-one with the second number of target Doppler data points. The target spectral information may include intensity coefficients corresponding one-to-one with the second number of target Doppler data points. The spectral components in the target spectral information are consistent with the spectral components in the reference spectral information. The intensity coefficient corresponding to each target Doppler data point in the target spectral information is determined based on the product of the corresponding time weighting coefficient and the intensity coefficient in the reference spectral information. For the reference spectral information corresponding to the reference data point sequence before the time period to be filled, the second number of target Doppler data points are arranged in chronological order of sampling time, and the corresponding time weighting coefficients gradually decrease. For the reference spectral information corresponding to the reference data point sequence after the time period to be filled, the second number of target Doppler data points are arranged in chronological order of sampling time, and the corresponding time weighting coefficients gradually increase.

[0064] In one embodiment, the target Doppler data may include a second number (L) of target Doppler data points. The size of the second number can be arbitrary and can be set as needed; this invention does not limit this. The preset time weighting coefficient may include a time weighting coefficient W(s) corresponding one-to-one with each of the L target Doppler data points. Wherein, s∈[1,L]. The target spectral information may include intensity coefficients corresponding one-to-one with each of the L target Doppler data points. The spectral components in the reference spectral information can be used as the spectral components in the target spectral information. Multiplying any time weighting coefficient by the corresponding intensity coefficient in the reference spectral information determines the intensity coefficient corresponding to the target Doppler data point in the target spectral information. Multiple target Doppler data points can be sequentially represented as d1, d2, d3, ..., d... LFor the reference spectrum information corresponding to the first reference data point sequence before the filling period, the L target Doppler data points are arranged in order of sampling time from front to back, and the corresponding time weight coefficients gradually decrease. That is, for the reference spectrum information corresponding to the first reference data point sequence, the time weight coefficient corresponding to d1 is greater than the time weight coefficient corresponding to d2, the time weight coefficient corresponding to d2 is greater than the time weight coefficient corresponding to d3, ..., the time weight coefficient corresponding to d7 is greater than the time weight coefficient corresponding to d8. For the reference spectrum information corresponding to the second reference data point sequence after the filling period, the L target Doppler data points are arranged in order of sampling time from front to back, and the corresponding time weight coefficients gradually increase. Referring to the above description of "time weight coefficients gradually decreasing," those skilled in the art can understand the meaning of "time weight coefficients gradually increasing," and for the sake of brevity, it will not be elaborated further here.

[0065] According to the above technical solution, for the reference spectrum information corresponding to the first reference data point sequence, the time weighting coefficient gradually decreases with the sampling time of the target Doppler data points, thus gradually reducing the proportion of the reference spectrum information corresponding to the first reference data point sequence in the target spectrum information. Conversely, this gradually increases the proportion of the reference spectrum information corresponding to the second reference data point sequence in the target spectrum information. This solution makes the obtained target Doppler data points more accurate.

[0066] For example, for the reference spectral information corresponding to the reference data point sequence located before the time period to be filled, the time weighting coefficients corresponding to the second number of target Doppler data points are determined based on the following formula:

[0067]

[0068] Where W(s) represents the time weighting coefficient; s represents the s-th target Doppler data point, s∈[1,L]; L represents the second quantity; for the reference spectrum information corresponding to the reference data point sequence located after the time period to be filled, the time weighting coefficient corresponding to the s-th target Doppler data point is equal to 1-W(s).

[0069] It should be noted that the calculation method of formula (4) is only an example and not a limitation of the present invention. For example, W(s) can also be calculated using other functions such as quadratic functions, where s is the independent variable of the function.

[0070] In one embodiment, according to formula (4) It can be seen that if L equals 8, then when s equals 1, W(s) equals 7 / 8; when s equals 4, W(s) equals 1 / 2; ..., when s equals 8, W(s) equals 0. The second number (L) target Doppler data points are arranged in chronological order of sampling time, with the corresponding time weighting coefficients gradually decreasing. Since the preset time weighting coefficients need to be applied to the pulse Doppler data filling process, they can be transformed to obtain the corresponding time weighting coefficient matrix for easier calculation. The time weighting coefficients corresponding to the L target Doppler data points can be represented by the following time weighting coefficient matrix:

[0071]

[0072] Where W1 can represent the time weight coefficient matrix corresponding to the first reference data point sequence, and W2 can represent the time weight coefficient matrix corresponding to the second reference data point sequence.

[0073] According to the above technical solution, the time weighting coefficients corresponding to the second number of target Doppler data points are determined by the above formula, which eliminates the need for complex calculations and has high computational efficiency.

[0074] For example, the reference Doppler data may include a sequence of reference data points before the period to be filled and a sequence of reference data points after the period to be filled, each sequence of reference data points may include one or more Doppler data points; the preset time weighting coefficient may include the time weighting coefficient corresponding to each reference data point sequence; predicting the Doppler data to be filled within the period to be filled based on the reference spectrum information and the preset time weighting coefficient to obtain the target Doppler data includes: for each target Doppler data point in the target Doppler data, for each reference data point sequence, determining the predicted signal intensity corresponding to the reference data point sequence based on the time weighting coefficient corresponding to the reference data point sequence and the reference spectrum information corresponding to the reference data point sequence, and summing the predicted signal intensities corresponding to all reference data point sequences to obtain the signal intensity of the target Doppler data point.

[0075] In one embodiment, the first and second reference data point sequences have been described in detail in the preceding embodiments, and for the sake of brevity, they will not be repeated here. For each target Doppler data point in the target Doppler data, based on the time weighting coefficient corresponding to each reference data point sequence and the reference spectrum information corresponding to that reference data point sequence, the predicted signal intensity corresponding to that reference data point sequence can be determined. By adding the predicted signal intensities corresponding to all reference data point sequences, the signal intensity of the current target Doppler data point can be determined.

[0076] According to the above technical solution, for each target Doppler data point in the target Doppler data, the signal intensity of that target Doppler data point can be obtained by summing the predicted signal intensities corresponding to each reference data point sequence. This solution can combine reference data point sequences before and after the time period to be filled to obtain the signal intensities of multiple target Doppler data points, providing more comprehensive information and helping to further improve the accuracy of the determined target Doppler data points.

[0077] For example, predicting the Doppler data to be filled within the time period based on reference spectrum information and a preset time weighting coefficient to obtain the target Doppler data may include: determining the signal strength of each target Doppler data point in the target Doppler data using the following formula:

[0078]

[0079]

[0080] Among them, g1, g2, g3, ..., g L Let I represent the signal intensity corresponding to the 1st, 2nd, 3rd, ..., Lth target Doppler data points, respectively; let Q represent the real part of each target Doppler data point; let Q represent the imaginary part of each target Doppler data point; and let K represent the intensity coefficient corresponding to the real part of the reference Doppler data point in the reference data point sequence preceding the time period to be filled. K represents the intensity coefficient corresponding to the real part of the reference Doppler data point in the reference data point sequence following the time period to be filled. Q This represents the intensity coefficient corresponding to the imaginary part of the reference Doppler data points in the reference data point sequence preceding the time period to be filled. Let W1 represent the intensity coefficient corresponding to the imaginary part of the reference Doppler data point in the reference data point sequence after the time period to be filled, and W2 represent the time weight coefficients corresponding to the reference data point sequences before and after the time period to be filled, respectively.

[0081] According to the above technical solution, the signal strength of each target Doppler data point in the target Doppler data is determined by the above formula (7). This solution has a small computational load and can ensure the efficiency and accuracy of acquiring target Doppler data.

[0082] The following describes an illustrative process for determining pulsed Doppler filling data according to an embodiment of the present invention. In the following process, it is assumed that the number of Doppler data points before the filling period is the same as the number of Doppler data points after the filling period, both being N.

[0083] Assume the sampling rate of the Doppler data is fs The N complex points located before the time period to be filled are Y = [I1 + jQ1, I2 + jQ2, ..., I...]. N +jQ N ] T =I+jQ, where N complex points are located after the time period to be filled. Suppose the time series corresponding to N complex points is T = [t1, t2, ..., t] N ] T =[1 / f s ,2 / f s , ..., N / f s ] T .

[0084] Step S1: Calculate the spectrum analysis matrix:

[0085]

[0086]

[0087] Step S2: Calculate the spectral analysis results (i.e., the intensity coefficient matrix): K = (C T C) -1 C T I, K Q =(C Q T C Q ) -1 C Q T Q,

[0088] Step S3: Calculate the time weighting function:

[0089] Step S4: Calculate the time weighting coefficient matrix: Calculate the time series within the period to be filled: T g = [t1, t2, ... t L ] T =[1 / f s ,2 / f s , ..., L / f s ] T .

[0090] Step S5: Calculate the target Doppler data point sequence:

[0091] According to another aspect of the present invention, a pulse Doppler filling data determination apparatus is also provided. Figure 3A schematic diagram of a pulsed Doppler fill data determination apparatus 300 according to an embodiment of the present invention is shown, as follows: Figure 3 As shown, the device 300 may include an acquisition module 310, an analysis module 320, and a prediction module 330.

[0092] The acquisition module 310 is used to acquire reference Doppler data located before and / or after the time period to be filled, the reference Doppler data being obtained by scanning with a Doppler ultrasound device.

[0093] Analysis module 320 is used to perform spectral analysis on reference Doppler data to obtain corresponding reference spectral information.

[0094] The prediction module 330 is used to predict the Doppler data to be filled in the time period based on the reference spectrum information to obtain the target Doppler data, wherein the target spectrum information corresponding to the target Doppler data is determined based on the reference spectrum information.

[0095] Those skilled in the art can understand the specific implementation scheme and beneficial effects of the pulse Doppler filling data determination device by reading the above description of the pulse Doppler filling data determination method 100. For the sake of brevity, it will not be described in detail here.

[0096] For example, the reference Doppler data includes a sequence of reference data points before the period to be filled and / or a sequence of reference data points after the period to be filled, each sequence of reference data points including one or more Doppler data points; the reference spectral information corresponding to each reference data point sequence includes the spectral components corresponding to the reference data point sequence and the intensity coefficient of each spectral component.

[0097] For example, the analysis module 320 may include an analysis submodule. The analysis submodule is configured to perform the following specific analysis operations for each reference data point sequence: divide the frequency range corresponding to each Doppler data point in the reference data point sequence into a first number of segments to obtain a first number of frequency components; determine a first number of intensity coefficients corresponding one-to-one with the first number of frequency components by minimizing the difference between the signal intensity of each first Doppler data point in the fitted data sequence and the signal intensity of the corresponding second Doppler data point in the reference data point sequence; wherein the fitted data sequence includes first Doppler data points corresponding one-to-one with the second Doppler data points in the reference data point sequence, and the signal intensity of the first Doppler data points is determined based on the first number of frequency components, the first number of intensity coefficients, and the sampling time of the corresponding second Doppler data points; the first number is an integer greater than or equal to 1.

[0098] For example, the signal strength of each Doppler data point in the reference data point sequence includes a real part and an imaginary part, and a specific analysis operation is performed on the real part and the imaginary part respectively; when the analysis submodule performs a specific analysis operation on the real part, the signal strength of the first Doppler data point is the real part of the first Doppler data point, the signal strength of the second Doppler data point is the real part of the second Doppler data point, the first number of frequency components are the first number of frequency components corresponding to the real part, and the first number of intensity coefficients are the first number of intensity coefficients corresponding to the real part; when the analysis submodule performs a specific analysis operation on the imaginary part, the signal strength of the first Doppler data point is the imaginary part of the first Doppler data point, the signal strength of the second Doppler data point is the imaginary part of the second Doppler data point, the first number of frequency components are the first number of frequency components corresponding to the imaginary part, and the first number of intensity coefficients are the first number of intensity coefficients corresponding to the imaginary part.

[0099] For example, the prediction module 330 may include a prediction submodule. The prediction submodule is used to predict the Doppler data to be filled within the time period based on reference spectrum information and a preset time weighting coefficient, to obtain target Doppler data; wherein the signal strength of the target Doppler data point in the target Doppler data is determined based on the target spectrum information and the sampling time corresponding to the target Doppler data point, and the target spectrum information is determined by adjusting the reference spectrum information with a preset time weighting coefficient, the time weighting coefficient being related to the sampling time corresponding to the target Doppler data point.

[0100] For example, the reference spectrum information includes spectral components and intensity coefficients for each spectral component; the target Doppler data includes a second number of target Doppler data points, where the second number is an integer greater than or equal to 1; the preset time weighting coefficients include time weighting coefficients corresponding one-to-one with the second number of target Doppler data points; the target spectrum information includes intensity coefficients corresponding one-to-one with the second number of target Doppler data points; and the spectral components in the target spectrum information are consistent with the spectral components in the reference spectrum information. The intensity coefficients corresponding to each target Doppler data point in the target spectrum information are determined based on the product of the corresponding time weighting coefficient and the intensity coefficients in the reference spectrum information. For the reference spectrum information corresponding to a sequence of reference data points before the time period to be filled, the second number of target Doppler data points are arranged in chronological order of sampling time, and the corresponding time weighting coefficients gradually decrease. For the reference spectrum information corresponding to a sequence of reference data points after the time period to be filled, the second number of target Doppler data points are arranged in chronological order of sampling time, and the corresponding time weighting coefficients gradually increase.

[0101] For example, for the reference spectral information corresponding to the reference data point sequence located before the time period to be filled, the time weighting coefficients corresponding to the second number of target Doppler data points are determined based on the following formula:

[0102]

[0103] Where W(s) represents the time weighting coefficient, s represents the s-th target Doppler data point, s∈[1,L], and L represents the second quantity; for the reference spectrum information corresponding to the reference data point sequence located after the time period to be filled, the time weighting coefficient corresponding to the s-th target Doppler data point is equal to 1-W(s).

[0104] For example, the reference Doppler data includes a sequence of reference data points before the time period to be filled and a sequence of reference data points after the time period to be filled, each reference data point sequence including one or more Doppler data points. The preset time weighting coefficients include time weighting coefficients corresponding to each reference data point sequence. The prediction module 330 may include a first determining submodule. The first determining submodule is used to, for each target Doppler data point in the target Doppler data, for each reference data point sequence, determine the predicted signal intensity corresponding to the reference data point sequence based on the time weighting coefficients corresponding to the reference data point sequence and the reference spectrum information corresponding to the reference data point sequence, and sum the predicted signal intensities corresponding to all reference data point sequences to obtain the signal intensity of the target Doppler data point.

[0105] For example, the prediction module 330 may include a second determining submodule. The second determining submodule is configured to determine the signal strength of each target Doppler data point in the target Doppler data using the following formula:

[0106]

[0107] Among them, g1, g2, g3, ..., g L Let I represent the signal intensity corresponding to the 1st, 2nd, 3rd, ..., Lth target Doppler data points, respectively; let Q represent the real part of each target Doppler data point; let Q represent the imaginary part of each target Doppler data point; and let K represent the intensity coefficient corresponding to the real part of the reference Doppler data point in the reference data point sequence preceding the time period to be filled. K represents the intensity coefficient corresponding to the real part of the reference Doppler data point in the reference data point sequence following the time period to be filled. Q This represents the intensity coefficient corresponding to the imaginary part of the reference Doppler data points in the reference data point sequence preceding the time period to be filled. W1 represents the intensity coefficient corresponding to the imaginary part of the reference Doppler data points in the reference data point sequence after the time period to be filled, and W2 represents the time weight coefficients corresponding to the reference data point sequences before and after the time period to be filled, respectively.

[0108] According to another aspect of the present invention, an electronic device is also provided. Figure 4 A schematic block diagram of an electronic device according to an embodiment of the present invention is shown. Figure 4 As shown, the electronic device 400 includes a processor 410 and a memory 420. The memory 420 stores a computer program. The computer program instructions are executed by the processor 410 to perform the pulse Doppler filling data determination method described above.

[0109] The aforementioned electronic devices can be ultrasound devices, such as desktop ultrasound diagnostic devices, portable ultrasound diagnostic devices (e.g., laptop ultrasound diagnostic devices or tablet ultrasound diagnostic devices), and also ultrasound imaging workstations.

[0110] According to another aspect of the present invention, a storage medium storing a computer program / instructions is also provided. The storage medium may include, for example, a storage component of a tablet computer, a hard disk of a personal computer, an erasable programmable read-only memory (EPROM), a portable read-only memory (CD-ROM), a USB memory, or any combination of the above storage media. The storage medium may be any combination of one or more computer-readable storage media. The computer program / instructions are used by a processor during runtime to execute the above-described pulse Doppler filling data determination method.

[0111] Those skilled in the art can understand the specific implementation scheme of the above-mentioned electronic device and storage medium by reading the relevant description of the pulse Doppler filling data determination method, and will not be repeated here for the sake of brevity.

[0112] Although exemplary embodiments have been described herein with reference to the accompanying drawings, it should be understood that the above exemplary embodiments are merely illustrative and are not intended to limit the scope of the invention. Various changes and modifications can be made therein by those skilled in the art without departing from the scope and spirit of the invention. All such changes and modifications are intended to be included within the scope of the invention as claimed in the appended claims.

[0113] Those skilled in the art will recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementations should not be considered beyond the scope of this invention.

[0114] In the several embodiments provided by this invention, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For instance, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another device, or some features may be ignored or not executed.

[0115] Numerous specific details are set forth in the specification provided herein. However, it will be understood that embodiments of the invention may be practiced without these specific details. In some instances, well-known methods, structures, and techniques have not been shown in detail so as not to obscure the understanding of this specification.

[0116] Similarly, it should be understood that, in order to streamline the invention and aid in understanding one or more of the various aspects of the invention, features of the invention are sometimes grouped together in a single embodiment, figure, or description thereof in the description of exemplary embodiments of the invention. However, this approach should not be construed as reflecting an intention that the claimed invention requires more features than are expressly recited in each claim. Rather, as reflected in the corresponding claims, its inventive point lies in solving the corresponding technical problem with fewer features than all of those in a single disclosed embodiment. Therefore, the claims following the detailed description are hereby expressly incorporated into that detailed description, wherein each claim itself is a separate embodiment of the invention.

[0117] Those skilled in the art will understand that, apart from the mutual exclusion of features, all features disclosed in this specification (including the accompanying claims, abstract, and drawings) and all processes or units of any method or apparatus so disclosed can be combined in any combination. Unless otherwise expressly stated, each feature disclosed in this specification (including the accompanying claims, abstract, and drawings) may be replaced by an alternative feature that serves the same, equivalent, or similar purpose.

[0118] Furthermore, those skilled in the art will understand that although some embodiments described herein include certain features but not others included in other embodiments, combinations of features from different embodiments are intended to be within the scope of the invention and form different embodiments. For example, in the claims, any of the claimed embodiments can be used in any combination.

[0119] The various component embodiments of the present invention can be implemented in hardware, or as software modules running on one or more processors, or a combination thereof. Those skilled in the art will understand that microprocessors or digital signal processors (DSPs) can be used in practice to implement some or all of the functions of some modules in the pulse Doppler filling data determination device according to embodiments of the present invention. The present invention can also be implemented as a device program (e.g., a computer program and computer program product) for performing part or all of the methods described herein. Such programs implementing the present invention can be stored on a computer-readable medium or can be in the form of one or more signals. Such signals can be downloaded from an Internet website, provided on a carrier signal, or provided in any other form.

[0120] It should be noted that the above embodiments are illustrative of the invention and not restrictive, and that those skilled in the art can devise alternative embodiments without departing from the scope of the appended claims. In the claims, any reference signs placed between parentheses should not be construed as limiting the claims. The word "comprising" does not exclude the presence of elements or steps not listed in the claims. The word "a" or "an" preceding an element does not exclude the presence of a plurality of such elements. The invention can be implemented by means of hardware comprising several different elements and by means of a suitably programmed computer. In the unit claims enumerating several means, several of these means may be embodied by the same item of hardware. The use of the words first, second, and third, etc., does not indicate any order. These words can be interpreted as names.

[0121] The above description is merely a specific embodiment of the present invention or an explanation of that embodiment. The scope of protection of the present invention is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. The scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A method for determining pulse Doppler filling data, characterized in that, The method includes: Acquire reference Doppler data located before and / or after the time period to be filled, the reference Doppler data being obtained by scanning with a Doppler ultrasound device; Spectral analysis is performed on the reference Doppler data to obtain the corresponding reference spectral information; Based on the reference spectrum information, the Doppler data to be filled in the time period is predicted to obtain the target Doppler data, wherein the target spectrum information corresponding to the target Doppler data is determined based on the reference spectrum information.

2. The method as described in claim 1, characterized in that, The reference Doppler data includes a sequence of reference data points preceding the period to be filled and / or a sequence of reference data points following the period to be filled, each sequence of reference data points including one or more Doppler data points; The reference spectrum information corresponding to each reference point data sequence includes the spectral components corresponding to that reference point data sequence and the intensity coefficient of each spectral component.

3. The method as described in claim 2, characterized in that, The step of performing spectral analysis on the reference Doppler data to obtain corresponding reference spectral information includes: For each reference data point sequence, perform the following specific analysis operations: The frequency range corresponding to each Doppler data point in the reference data point sequence is divided into a first number of segments to obtain a first number of frequency components; The first number of intensity coefficients corresponding one-to-one with the first number of frequency components are determined by minimizing the difference between the signal intensity of each first Doppler data point in the fitted data sequence and the signal intensity of the corresponding second Doppler data point in the reference data sequence. The fitted data sequence includes first Doppler data points that correspond one-to-one with the second Doppler data points in the reference data point sequence. The signal intensity of the first Doppler data points is determined based on the first number of frequency components, the first number of intensity coefficients, and the sampling time of the corresponding second Doppler data points. The first number is an integer greater than or equal to 1.

4. The method as described in claim 3, characterized in that, The signal strength of each Doppler data point in the reference data point sequence includes a real part and an imaginary part, and the specific analysis operation is performed on the real part and the imaginary part respectively; When performing the specific analysis operation on the real part, the signal intensity of the first Doppler data point is the real part of the first Doppler data point, the signal intensity of the second Doppler data point is the real part of the second Doppler data point, the first number of frequency components are the first number of frequency components corresponding to the real part, and the first number of intensity coefficients are the first number of intensity coefficients corresponding to the real part. When performing the specific analysis operation on the imaginary part, the signal strength of the first Doppler data point is the imaginary part of the first Doppler data point, the signal strength of the second Doppler data point is the imaginary part of the second Doppler data point, the first number of frequency components are the first number of frequency components corresponding to the imaginary part, and the first number of intensity coefficients are the first number of intensity coefficients corresponding to the imaginary part.

5. The method according to any one of claims 1-4, characterized in that, The step of predicting the Doppler data to be filled within the time period based on the reference spectrum information to obtain the target Doppler data includes: Based on the reference spectrum information and the preset time weighting coefficient, the Doppler data to be filled in the time period is predicted to obtain the target Doppler data. The signal strength of the target Doppler data point in the target Doppler data is determined based on the target spectrum information and the sampling time corresponding to the target Doppler data point. The target spectrum information is determined by adjusting the reference spectrum information with the preset time weighting coefficient, and the time weighting coefficient is related to the sampling time corresponding to the target Doppler data point.

6. The method as described in claim 5, characterized in that, The reference spectrum information includes spectral components and intensity coefficients for each spectral component; the target Doppler data includes a second number of target Doppler data points, where the second number is an integer greater than or equal to 1; the preset time weighting coefficients include time weighting coefficients that correspond one-to-one with the second number of target Doppler data points; the target spectrum information includes intensity coefficients that correspond one-to-one with the second number of target Doppler data points; and the spectral components in the target spectrum information are consistent with the spectral components in the reference spectrum information. The intensity coefficient corresponding to each target Doppler data point in the target spectral information is determined based on the product of the corresponding time weighting coefficient and the intensity coefficient in the reference spectral information. Specifically, for the reference spectrum information corresponding to the reference data point sequence located before the time period to be filled, the second number of target Doppler data points are arranged in order from front to back according to the sampling time, and the corresponding time weight coefficient gradually decreases; for the reference spectrum information corresponding to the reference data point sequence located after the time period to be filled, the second number of target Doppler data points are arranged in order from front to back according to the sampling time, and the corresponding time weight coefficient gradually increases.

7. The method as described in claim 6, characterized in that, For the reference spectral information corresponding to the reference data point sequence located before the time period to be filled, the time weighting coefficients corresponding to the second number of target Doppler data points are determined based on the following formula: Where W(s) represents the time weighting coefficient, s represents the s-th target Doppler data point, s∈[1,L], and L represents the second quantity; For the reference spectrum information corresponding to the reference data point sequence located after the time period to be filled, the time weight coefficient corresponding to the s-th target Doppler data point is equal to 1-W(s).

8. The method as described in claim 5, characterized in that, The reference Doppler data includes a sequence of reference data points before the time period to be filled and a sequence of reference data points after the time period to be filled, each reference data point sequence including one or more Doppler data points; the preset time weighting coefficient includes a time weighting coefficient corresponding to each reference data point sequence. The step of predicting the Doppler data to be filled within the time period based on the reference spectrum information and a preset time weighting coefficient to obtain the target Doppler data includes: For each target Doppler data point in the target Doppler data For each reference data point sequence, based on the time weighting coefficient corresponding to the reference data point sequence and the reference spectrum information corresponding to the reference data point sequence, the predicted signal intensity corresponding to the reference data point sequence is determined, and the predicted signal intensities corresponding to all reference data point sequences are summed to obtain the signal intensity of the target Doppler data point.

9. The method as described in claim 8, characterized in that, The step of predicting the Doppler data to be filled within the time period based on the reference spectrum information and a preset time weighting coefficient to obtain the target Doppler data includes: The signal strength of each target Doppler data point in the target Doppler data is determined by the following formula: Among them, g1, g2, g3, ..., g L Let I represent the signal intensity corresponding to the 1st, 2nd, 3rd, ..., Lth target Doppler data points in the target Doppler data, Q represent the imaginary part of each target Doppler data point, and K represent the intensity coefficient corresponding to the real part of the reference Doppler data point in the reference data point sequence preceding the time period to be filled. K represents the intensity coefficient corresponding to the real part of the reference Doppler data point in the reference data point sequence following the time period to be filled. Q This represents the intensity coefficient corresponding to the imaginary part of the reference Doppler data points in the reference data point sequence preceding the time period to be filled. W1 represents the intensity coefficient corresponding to the imaginary part of the reference Doppler data points in the reference data point sequence located after the time period to be filled, and W2 represents the time weight coefficients corresponding to the reference data point sequences located before and after the time period to be filled, respectively.

10. A pulse Doppler filling data determination device, characterized in that, The device includes: The acquisition module is used to acquire reference Doppler data located before and / or after the time period to be filled, the reference Doppler data being obtained by scanning with a Doppler ultrasound device; The analysis module is used to perform spectral analysis on the reference Doppler data to obtain the corresponding reference spectral information; The prediction module is used to predict the Doppler data to be filled in the time period based on the reference spectrum information to obtain the target Doppler data, wherein the target spectrum information corresponding to the target Doppler data is determined based on the reference spectrum information.

11. An electronic device comprising a processor and a memory, characterized in that, The memory stores a computer program, the computer program instructions of which, when executed by the processor, are used to perform the pulse Doppler filling data determination method as described in any one of claims 1-9.

12. A storage medium storing a computer program / instructions, characterized in that, The computer program / instructions are used to perform the pulse Doppler filling data determination method as described in any one of claims 1-9 when the program is running.

Images