Ultrasonic Doppler Envelope Determination Method, Device, Computer Equipment, and Storage Medium
By obtaining the power spectral density signal in the ultrasonic spectrum and performing integration processing, determining the maximum frequency point in the forward and reverse directions, the problems of complex calculations and insufficient adaptability of traditional methods are solved, and efficient envelope extraction is achieved.
Patent Information
- Application Number
- CN202210257856.1
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2022-03-16
- Publication Date
- 2025-08-01
- Estimated Expiration
- 2042-03-16
AI Technical Summary
The traditional ultrasonic Doppler spectrum envelope extraction technology has a large calculation volume and high complexity, and cannot accurately handle the positive and reverse blood flow, especially in the presence of reverse blood flow velocity.
By obtaining the power spectral density signals corresponding to each column of signals in the ultrasound spectrum, integrating the forward and reverse maximum frequency points are determined, and connecting these points to form an envelope, the degree of clinging of the waveform is adjusted using the integral curve and the preset threshold.
It realizes the determination of forward and reverse envelopes with simple calculation process and high efficiency, and can accurately obtain the waveform envelope of the ultrasound spectrum to adapt to different blood flow conditions.
Smart Images
Figure CN114642446B_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to the technical field of medical signal processing, and particularly to a method, apparatus, computer device, storage medium, and computer program product for determining an ultrasound Doppler envelope. Background Art
[0002] In recent years, with the rapid development of ultrasound Doppler technology, ultrasound Doppler can non-invasively detect the blood flow condition in human blood vessels, and thus is widely used to diagnose vascular diseases. A large number of studies have shown that vascular diseases can be reflected from the changes in the maximum frequency curve of Doppler signals. Therefore, whether an accurate positive and negative maximum frequency curve (i.e., spectral envelope) of the ultrasound Doppler blood flow signal can be obtained plays an important role in the diagnosis of vascular diseases. The traditional spectral envelope is drawn by experienced ultrasound physicians, which is not only inefficient but also requires a large amount of manual intervention.
[0003] In related technologies, the envelope extraction technology of Doppler spectra mainly includes two types: one is to use image processing technology, such as edge detection of images, etc. to obtain the envelope of the Doppler spectrum. However, this image processing-based method has problems such as large computational amount, high complexity, and low real-time performance. The other is to use the frequency characteristics of the power spectral density integral curve of the ultrasound Doppler blood flow signal to estimate the maximum frequency. Although this digital frequency estimation-based method has a small computational amount and high computational efficiency, it also has certain limitations. If there is a reverse blood flow velocity, it cannot work properly. Summary of the Invention
[0004] Based on this, it is necessary to provide a method, apparatus, computer device, storage medium, and computer program product for determining an ultrasound Doppler envelope that can simultaneously determine positive and negative waveforms in view of the above technical problems.
[0005] The present application provides a method for determining an ultrasound Doppler envelope. The method includes:
[0006] Obtain the power spectral density signal corresponding to each column of signals in the ultrasound map;
[0007] Perform integral processing on the power spectral density signal to obtain the positive maximum frequency point and the negative maximum frequency point in the column signal corresponding to the power spectral density signal;
[0008] Connect the positive maximum frequency point and the negative maximum frequency point corresponding to each column of signals in the ultrasound map to obtain the waveform envelope line corresponding to the ultrasound map.
[0009] Further, the integrating the power spectral density signal to obtain the forward maximum frequency point and the reverse maximum frequency point in the column signal corresponding to the power spectral density signal includes: performing a first integration process on the power spectral density signal based on a first frequency interval corresponding to the power spectral density signal to obtain a corresponding first integration curve; determining the maximum energy point of the power spectral density signal according to the first integration curve; determining two minimum energy points of the power spectral density signal according to the maximum energy point; performing a second integration process on the power spectral density signal in the corresponding interval based on a second frequency interval corresponding to the two minimum energy points to obtain a corresponding second integration curve; and determining the forward maximum frequency point and the reverse maximum frequency point in the column signal corresponding to the power spectral density signal according to the second integration curve.
[0010] Further, the determining the maximum energy point of the power spectral density signal according to the first integration curve includes: obtaining the first and last endpoints of the first integration curve, connecting the two endpoints by a first reference line, and the first reference line is a straight line; obtaining a first intersection point between the first reference line and the first integration curve, and taking the first intersection point as the maximum energy point of the power spectral density signal.
[0011] Further, the determining two minimum energy points of the power spectral density signal according to the maximum energy point includes: dividing the power spectral density signal into two signal intervals according to the maximum energy point; respectively obtaining the point with the minimum ordinate in each signal interval as the minimum energy point of the power spectral density signal.
[0012] Further, the determining the forward maximum frequency point and the reverse maximum frequency point in the column signal corresponding to the power spectral density signal according to the second integration curve includes: obtaining a first endpoint at the head and a second endpoint at the tail of the second integration curve, connecting the first endpoint and the second endpoint by a second reference line, and the second reference line is a straight line; determining a second intersection point, a forward frequency point and a reverse frequency point between the second reference line and the second integration curve, where the forward frequency point is the frequency point corresponding to the maximum forward distance from the second integration curve to the second reference line, and the reverse frequency point is the frequency point corresponding to the maximum reverse distance from the second integration curve to the second reference line; obtaining a first average value difference between the two-side signals divided by the forward frequency point between the first endpoint and the second intersection point, and when the first average value difference is greater than or equal to a preset threshold, determining the forward frequency point as the forward maximum frequency point; obtaining a second average value difference between the two-side signals divided by the reverse frequency point between the second intersection point and the second endpoint, and when the second average value difference is greater than or equal to a preset threshold, determining the reverse frequency point as the reverse maximum frequency point.
[0013] Further, the method further includes: when the first average value difference is less than a preset threshold, moving the forward frequency point in the direction of the second intersection point, and returning to execute the step of obtaining the first average value difference between the two side signals divided by the forward frequency point between the first end point and the second intersection point; when the second average value difference is less than a preset threshold, moving the reverse frequency point in the direction of the second intersection point, and returning to execute the step of obtaining the second average value difference between the two side signals divided by the reverse frequency point between the second intersection point and the second end point.
[0014] The present application further provides an ultrasonic Doppler envelope determination device, and the device includes:
[0015] A power spectral density signal acquisition module, configured to acquire a power spectral density signal corresponding to each column signal in the ultrasonic spectrogram;
[0016] A frequency point determination module, configured to perform an integration process on the power spectral density signal to obtain a forward maximum frequency point and a reverse maximum frequency point in the column signal corresponding to the power spectral density signal;
[0017] An envelope line acquisition module, configured to connect the forward maximum frequency point and the reverse maximum frequency point corresponding to each column signal in the ultrasonic spectrogram to obtain an envelope line corresponding to the ultrasonic spectrogram.
[0018] The present application further provides a computer device. The computer device includes a memory and a processor, the memory stores a computer program, and when the processor executes the computer program, the steps in any one of the above embodiments are implemented.
[0019] The present application further provides a computer-readable storage medium. The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps in any one of the above embodiments are implemented.
[0020] The present application further provides a computer program product. The computer program product includes a computer program, and when the computer program is executed by a processor, the steps in any one of the above embodiments are implemented.
[0021] The above ultrasonic Doppler envelope determination method, device, computer device, storage medium, and computer program product obtain the power spectral density signal corresponding to each column of signals in the ultrasonic spectrogram, perform integral processing on the power spectral density signal to obtain the forward maximum frequency point and the reverse maximum frequency point in the column signal corresponding to the power spectral density signal, and then connect the forward maximum frequency point and the reverse maximum frequency point corresponding to each column of signals in the ultrasonic spectrogram, thereby obtaining the waveform envelope line corresponding to the ultrasonic spectrogram. In this embodiment, by converting the ultrasonic signal into a power spectral density signal and simultaneously determining the forward envelope line and the reverse envelope line corresponding to the ultrasonic spectrogram based on the power spectral density signal, it is possible to obtain a waveform envelope line composed of the forward envelope line and the reverse envelope line, and it has the advantages of simple calculation process and high calculation efficiency. Description of the Drawings
[0022] Figure 1 It is a schematic flowchart of the ultrasonic Doppler envelope determination method in this application.
[0023] Figure 2 It is a schematic diagram of signal conversion in this application.
[0024] Figure 3 It is a schematic flowchart of the step of performing integral processing on the power spectral density signal in this application.
[0025] Figure 4 It is a schematic diagram of the first integral curve in this application.
[0026] Figure 5 It is a schematic diagram of determining the maximum energy point and the minimum energy point in the power spectral density signal in this application.
[0027] Figure 6 It is a schematic diagram of the second integral curve in this application.
[0028] Figure 7 It is a schematic diagram of the waveform envelope line in this application.
[0029] Figure 8 It is a structural block diagram of the ultrasonic Doppler envelope determination device in this application.
[0030] Figure 9 It is an internal structure diagram of the computer device in this application. Detailed Embodiments
[0031] To make the technical solutions and beneficial effects of the present invention more obvious and understandable, the following will be described in detail by listing specific embodiments. Among them, the drawings are not necessarily drawn to scale, and local features can be enlarged or reduced to more clearly show the details of the local features; unless otherwise defined, the technical and scientific terms used herein have the same meanings as those in the technical and scientific fields to which this application belongs.
[0032] This application provides an ultrasonic Doppler envelope determination method. In this embodiment, taking the application of this method to a terminal as an example, it can be understood that this method can also be applied to a server, and can also be applied to a system including a terminal and a server, and is implemented through the interaction between the terminal and the server. Specifically, as Figure 1 shown, the above method includes the following steps:
[0033] Step 102, obtain the power spectral density signal corresponding to each column of signals in the ultrasonic spectrogram.
[0034] Among them, the ultrasonic spectrogram is a spectrogram formed by ultrasonic signals collected by an ultrasonic Doppler probe as a sensor. Each column of signals in the ultrasonic spectrogram includes a number of signal points. For each signal point, it can be uniquely represented in the power spectrum based on the corresponding frequency and energy.
[0035] In this embodiment, as Figure 2 shown, for a number of signal points corresponding to any column of signals n in the ultrasonic spectrogram, according to the frequency and energy respectively corresponding to each signal point in this column of signals, it can be converted into a power spectral density signal S(n) represented by the power spectrum, where the abscissa of the power spectrum is the frequency and the ordinate is the energy.
[0036] Step 104, perform integral processing on the power spectral density signal to obtain the forward maximum frequency point and the reverse maximum frequency point in the column of signals corresponding to the power spectral density signal.
[0037] Among them, the forward maximum frequency point represents the maximum velocity point of the corresponding column of signals, and the reverse maximum frequency point represents the minimum velocity point of the corresponding column of signals. In this embodiment, by performing integral processing on the power spectral density signal S(n) corresponding to each column of signals n in the ultrasonic spectrogram, the forward maximum frequency point and the reverse maximum frequency point in each column of signals are obtained.
[0038] Step 106, connect the forward maximum frequency point and the reverse maximum frequency point corresponding to each column of signals in the ultrasonic spectrogram to obtain the waveform envelope line corresponding to the ultrasonic spectrogram.
[0039] Specifically, by sequentially connecting the forward maximum frequency points corresponding to each column of signals in the ultrasonic spectrogram, the corresponding forward envelope line is obtained. By sequentially connecting the reverse maximum frequency points corresponding to each column of signals in the ultrasonic spectrogram, the corresponding reverse envelope line is obtained. The waveform composed of the forward envelope line and the reverse envelope line is the waveform envelope line corresponding to the ultrasonic spectrogram.
[0040] In the above embodiments, by obtaining the power spectral density signals corresponding to each column of signals in the ultrasonic spectrogram and performing integral processing on the power spectral density signals, the forward maximum frequency point and the reverse maximum frequency point in the column signals corresponding to the power spectral density signals are obtained. Then, the forward maximum frequency points and the reverse maximum frequency points respectively corresponding to each column of signals in the ultrasonic spectrogram are connected, so as to obtain the waveform envelope corresponding to the ultrasonic spectrogram. In this embodiment, by converting the ultrasonic signal into a power spectral density signal and simultaneously determining the forward envelope and the reverse envelope corresponding to the ultrasonic spectrogram based on the power spectral density signal, the waveform envelope composed of the forward envelope and the reverse envelope is obtained, and it has the advantages of simple calculation process and high calculation efficiency.
[0041] In one embodiment, as Figure 3 shown, the steps of performing integral processing on the power spectral density signal to obtain the forward maximum frequency point and the reverse maximum frequency point in the column signal corresponding to the power spectral density signal may specifically include:
[0042] Step 302, performing a first integral processing on the power spectral density signal based on the first frequency interval corresponding to the power spectral density signal to obtain a corresponding first integral curve.
[0043] Wherein, the first frequency interval refers to the frequency interval corresponding to the power spectral density signal. The first integral processing is essentially a process of accumulating the gray levels of the column signals in the ultrasonic spectrogram from low frequency to high frequency. The first integral curve is the integral result obtained by performing the first integral processing on the power spectral density signal.
[0044] In this embodiment, as Figure 4 shown, by performing a first integral processing on the Figure 2 power spectral density signal S(n) therein, the corresponding first integral curve P(n) is obtained. Specifically, the first integral curve P(n) is a discrete data point curve obtained by integrating the power spectral density signal S(n) corresponding to the nth column signal in the ultrasonic spectrogram, that is, a discrete data point curve obtained by integrating the power spectral density signal S(n) as the frequency increases. It can be understood that for different column signals n in the ultrasonic spectrogram, the corresponding power spectral density signals S(n) may be different, and the first integral curves P(n) obtained after integral processing are also different.
[0045] Step 304, determining the maximum energy point of the power spectral density signal according to the first integral curve.
[0046] Since the ordinate of the power spectrum is energy, the maximum energy point is the Figure 2 point with the largest ordinate in the power spectrum as shown. Also, since the first integral curve P(n) is obtained by performing integral processing on the Figure 2 power spectral density signal S(n) as shown, therefore, in this embodiment, byFigure 4 The first integral curve P(n) shown can conveniently determine the maximum energy point of the power spectral density signal S(n).
[0047] Specifically, as Figure 4 shown, first determine the two end points of the first integral curve P(n), and then connect the two end points through the first reference line K1, where the first reference line K1 is a straight line. Further, obtain the first intersection point m1 between the first reference line K1 and the first integral curve P(n), and use the first intersection point m1 as the maximum energy point of the power spectral density signal.
[0048] Step 306: Determine two minimum energy points of the power spectral density signal according to the maximum energy point.
[0049] Since the maximum energy point is the point with the largest ordinate in the power spectrum, the minimum energy point is the point with the smallest ordinate in the power spectrum.
[0050] In this embodiment, the maximum energy point determined based on the above steps, that is, the first intersection point m1, divides the power spectral density signal S(n) into two signal intervals, and respectively obtains the point with the smallest ordinate in each signal interval as the minimum energy point of the power spectral density signal.
[0051] Specifically, as Figure 5 shown, the power spectral density signal S(n) is divided into regions through the position of the first intersection point m1 in the horizontal coordinate direction, obtaining the signal interval from S(1) to S(m1) and the signal interval from S(m1) to S(n) (where S(1) is the signal point corresponding to the coordinate origin, S(m1) is the abscissa corresponding to the m1 signal point, and S(n) is the abscissa corresponding to the signal point n far from the origin). And search for the point with the smallest ordinate in the signal interval from S(1) to S(m1) to obtain lowest_p, and search for the point with the smallest ordinate in the signal interval from S(m1) to S(n) to obtain lowest_r, thereby obtaining two minimum energy points lowest_p and lowest_r.
[0052] Step 308: Perform a second integration process on the power spectral density signal in the corresponding interval based on the second frequency interval corresponding to the two minimum energy points to obtain the corresponding second integral curve.
[0053] In this embodiment, after determining the two minimum energy points through the above steps, the interval composed of the abscissas corresponding to the two minimum energy points respectively is used as the second frequency interval, and the power spectral density signal corresponding to this interval is subjected to a second integration process, thereby obtaining the corresponding second integral curve.
[0054] Specifically, as Figure 5As shown, if the power spectral density signal corresponding to the interval from the lowest energy point lowest_p to lowest_r is S(p_r) (S(p_r) is the signal corresponding to the signal interval intercepted from S(n) from lowest_p to lowest_r), then the second integration process is performed on S(p_r) with the second frequency interval composed of lowest_p to lowest_r as the integration interval, so as to obtain Figure 6 the second integration curve P(m) as shown.
[0055] Step 310, determine the forward maximum frequency point and the reverse maximum frequency point in the column signal corresponding to the power spectral density signal according to the second integration curve.
[0056] Specifically, as Figure 6 shown, first obtain the first endpoint at the head and the second endpoint at the tail of the second integration curve P(m), and then connect the first endpoint and the second endpoint through the second reference line K2, where the second reference line K2 is a straight line. The first endpoint is the endpoint closest to the origin in the second integration curve P(m), and the second endpoint is the endpoint farthest from the origin in the second integration curve P(m).
[0057] Furthermore, determine the second intersection point m2, the forward frequency point locate_p, and the reverse frequency point locate_r between the second reference line K2 and the second integration curve P(m). Among them, the forward frequency point locate_p is the frequency point corresponding to the maximum forward distance from the second integration curve P(m) to the second reference line K2, and the reverse frequency point locate_r is the frequency point corresponding to the maximum reverse distance from the second integration curve P(m) to the second reference line K2.
[0058] Obtain the first average value difference between the two side signals divided by the forward frequency point locate_p between the first endpoint and the second intersection point m2, that is, obtain the first average value of the frequencies corresponding to all signal points between the first endpoint and the forward frequency point locate_p, and obtain the second average value of the frequencies corresponding to all signal points between the forward frequency point locate_p and the second intersection point m2, and calculate the difference between the first average value and the second average value, and this difference is the first average value difference. When the first average value difference is greater than or equal to the preset threshold, determine this forward frequency point locate_p as the forward maximum frequency point. Among them, the range of the preset threshold can be any value between 0 and 255.
[0059] Similarly, obtain the second average value difference between the two sides divided by the reverse frequency point locate_r between the second intersection point m2 and the second endpoint. When the second average value difference is greater than or equal to the preset threshold, determine this reverse frequency point locate_r as the reverse maximum frequency point.
[0060] When the first average value difference is less than the preset threshold, the positive frequency point is moved in the direction of the second intersection point, that is, locate_(p + 1) is determined as the positive frequency point, and the step of obtaining the first average value difference of the two-side signals divided by the positive frequency point between the first end point and the second intersection point is returned for execution. Until the first average value difference is greater than or equal to the preset threshold, the corresponding positive frequency point can be determined as the positive maximum frequency point.
[0061] Similarly, when the second average value difference is less than the preset threshold, the reverse frequency point is moved in the direction of the second intersection point, that is, locate_(r - 1) is determined as the reverse frequency point, and the step of obtaining the second average value difference of the two-side signals divided by the reverse frequency point between the second intersection point and the second end point is returned for execution. Until the second average value difference is greater than or equal to the preset threshold, the corresponding reverse frequency point can be determined as the reverse maximum frequency point.
[0062] By performing the processing as shown in Figure 3 on the power spectral density signals corresponding to each column of signals in the ultrasonic spectrogram respectively, the positive maximum frequency point and the reverse maximum frequency point corresponding to each column of signals can be obtained. By connecting the positive maximum frequency points corresponding to each column of signals in the ultrasonic spectrogram in sequence, the corresponding positive envelope is obtained. By connecting the reverse maximum frequency points corresponding to each column of signals in the ultrasonic spectrogram in sequence, the corresponding reverse envelope is obtained. The waveform composed of the positive envelope and the reverse envelope is the waveform envelope corresponding to the ultrasonic spectrogram (as shown in Figure 7 ).
[0063] In the above embodiment, by performing the first integration processing on the power spectral density signal, the maximum energy point and the minimum energy point are determined to achieve the purpose of filtering and denoising the signal. Furthermore, by performing the second integration processing on the filtered and denoised signal, the positive maximum frequency point and the reverse maximum frequency point of the signal corresponding to the corresponding column in the ultrasonic spectrogram can be determined more accurately. The tightness of the waveform can be flexibly adjusted by means of the preset threshold.
[0064] It should be understood that although the steps in the flowcharts involved in the above-described embodiments are shown in sequence according to the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless otherwise clearly stated in this article, there is no strict order restriction for the execution of these steps, and these steps can be executed in other orders. Moreover, at least a part of the steps in the flowcharts involved in the above-described embodiments may include multiple steps or multiple stages. These steps or stages are not necessarily executed at the same time, but can be executed at different times. The execution order of these steps or stages is not necessarily sequential, but can be executed alternately or in turn with at least a part of other steps or steps or stages in other steps.
[0065] Based on the same inventive concept, an embodiment of the present application further provides an ultrasonic Doppler envelope determination device for implementing the ultrasonic Doppler envelope determination method described above. The solution provided by this device for solving problems is similar to the solution described in the above method. Therefore, the specific limitations in one or more embodiments of the ultrasonic Doppler envelope determination device provided below can refer to the limitations on the ultrasonic Doppler envelope determination method in the above text, and will not be repeated here.
[0066] In one embodiment, as Figure 8 shown, an ultrasonic Doppler envelope determination device is provided, including: a power spectral density signal acquisition module 802, a frequency point determination module 804, and an envelope line acquisition module 806, where:
[0067] The power spectral density signal acquisition module 802 is configured to acquire the power spectral density signal corresponding to each column of signals in the ultrasonic spectrogram;
[0068] The frequency point determination module 804 is configured to perform an integration process on the power spectral density signal to obtain the forward maximum frequency point and the reverse maximum frequency point in the column signal corresponding to the power spectral density signal;
[0069] The envelope line acquisition module 806 is configured to connect the forward maximum frequency point and the reverse maximum frequency point corresponding to each column of signals in the ultrasonic spectrogram to obtain the envelope line corresponding to the ultrasonic spectrogram.
[0070] In one embodiment, the frequency point determination module may further include: a first integration processing unit configured to perform a first integration process on the power spectral density signal based on a first frequency interval corresponding to the power spectral density signal to obtain a corresponding first integration curve; a maximum energy point determination unit configured to determine a maximum energy point of the power spectral density signal according to the first integration curve; a minimum energy point determination unit configured to determine two minimum energy points of the power spectral density signal according to the maximum energy point; a second integration processing unit configured to perform a second integration process on the power spectral density signal in a corresponding interval based on a second frequency interval corresponding to the two minimum energy points to obtain a corresponding second integration curve; and a frequency point determination unit configured to determine a forward maximum frequency point and a reverse maximum frequency point in a column signal corresponding to the power spectral density signal according to the second integration curve.
[0071] In one embodiment, the maximum energy point determination unit is specifically further configured to: obtain two end points at the head and tail of the first integration curve, and connect the two end points by a first reference line, where the first reference line is a straight line; obtain a first intersection point between the first reference line and the first integration curve, and use the first intersection point as the maximum energy point of the power spectral density signal.
[0072] In one embodiment, the minimum energy point determination unit is specifically further configured to: divide the power spectral density signal into two signal intervals according to the maximum energy point; respectively obtain a point with the minimum ordinate in each signal interval as the minimum energy point of the power spectral density signal.
[0073] In one embodiment, the frequency point determination unit is specifically further configured to: obtain a first end point at the head and a second end point at the tail of the second integration curve, and connect the first end point and the second end point by a second reference line, where the second reference line is a straight line; determine a second intersection point, a forward frequency point, and a reverse frequency point between the second reference line and the second integration curve, where the forward frequency point is a frequency point corresponding to the maximum forward distance from the second integration curve to the second reference line, and the reverse frequency point is a frequency point corresponding to the maximum reverse distance from the second integration curve to the second reference line; obtain a first average value difference between two sides of the first end point to the second intersection point divided by the forward frequency point, and when the first average value difference is greater than or equal to a preset threshold, determine the forward frequency point as the forward maximum frequency point; obtain a second average value difference between two sides of the second intersection point to the second end point divided by the reverse frequency point, and when the second average value difference is greater than or equal to a preset threshold, determine the reverse frequency point as the reverse maximum frequency point.
[0074] In one embodiment, the frequency point determination unit is further specifically configured to: when the first average value difference is less than a preset threshold, move the forward frequency point in the direction of the second intersection point, and return to execute the step of obtaining the first average value difference between the two sides divided by the forward frequency point from the first endpoint to the second intersection point; when the second average value difference is less than a preset threshold, move the reverse frequency point in the direction of the second intersection point, and return to execute the step of obtaining the second average value difference between the two sides divided by the reverse frequency point from the second intersection point to the second endpoint.
[0075] Each module in the above ultrasonic Doppler envelope determination device can be implemented in whole or in part by software, hardware, and their combination. Each of the above modules can be embedded in the processor in the computer device in the form of hardware or be independent of it, or can be stored in the memory of the computer device in the form of software, so that the processor can call and execute the operations corresponding to each of the above modules.
[0076] In one embodiment, a computer device is provided. The computer device can be a terminal, and its internal structure diagram can be as Figure 9 shown. The computer device includes a processor, a memory, a communication interface, a display screen, and an input device connected through a system bus. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The communication interface of the computer device is used to communicate with an external terminal in a wired or wireless manner. The wireless manner can be achieved through WIFI, a mobile cellular network, NFC (Near Field Communication), or other technologies. When the computer program is executed by the processor, it implements an ultrasonic Doppler envelope determination method. The display screen of the computer device can be a liquid crystal display screen or an electronic ink display screen. The input device of the computer device can be a touch layer covering the display screen, or a button, a trackball, or a touchpad provided on the housing of the computer device, or an external keyboard, a touchpad, or a mouse, etc.
[0077] Those skilled in the art can understand that Figure 9 the structure shown in
[0078] is only a block diagram of some structures related to the solution of this application, and does not constitute a limitation on the computer device to which the solution of this application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine some components, or have different component arrangements. In one embodiment, a computer device is provided, including a memory and a processor. A computer program is stored in the memory. When the processor executes the computer program, the following steps are implemented:
[0079] Obtain the power spectral density signals corresponding to each column of signals in the ultrasonic spectrogram;
[0080] Perform integral processing on the power spectral density signals to obtain the forward maximum frequency point and the reverse maximum frequency point in the column signals corresponding to the power spectral density signals;
[0081] Connect the forward maximum frequency point and the reverse maximum frequency point corresponding to each column of signals in the ultrasonic spectrogram to obtain the waveform envelope corresponding to the ultrasonic spectrogram.
[0082] In one embodiment, when the processor executes the computer program, the following steps are further implemented: perform first integral processing on the power spectral density signals based on the first frequency interval corresponding to the power spectral density signals to obtain the corresponding first integral curve; determine the maximum energy point of the power spectral density signals according to the first integral curve; determine the two minimum energy points of the power spectral density signals according to the maximum energy point; perform second integral processing on the power spectral density signals in the corresponding interval based on the second frequency interval corresponding to the two minimum energy points to obtain the corresponding second integral curve; determine the forward maximum frequency point and the reverse maximum frequency point in the column signals corresponding to the power spectral density signals according to the second integral curve.
[0083] In one embodiment, when the processor executes the computer program, the following steps are further implemented: obtain the first and last endpoints of the first integral curve, and connect the two endpoints through a first reference line, where the first reference line is a straight line; obtain the first intersection point between the first reference line and the first integral curve, and use the first intersection point as the maximum energy point of the power spectral density signals.
[0084] In one embodiment, when the processor executes the computer program, the following steps are further implemented: divide the power spectral density signals into two signal intervals according to the maximum energy point; respectively obtain the points with the minimum ordinate in each signal interval as the minimum energy points of the power spectral density signals.
[0085] In one embodiment, when the processor executes the computer program, the following steps are further implemented: obtaining a first endpoint at the head and a second endpoint at the tail of the second integral curve, connecting the first endpoint and the second endpoint through a second reference line, where the second reference line is a straight line; determining a second intersection point, a forward frequency point, and a reverse frequency point between the second reference line and the second integral curve, where the forward frequency point is the frequency point corresponding to the maximum forward distance from the second integral curve to the second reference line, and the reverse frequency point is the frequency point corresponding to the maximum reverse distance from the second integral curve to the second reference line; obtaining a first average value difference between two side signals divided by the forward frequency point between the first endpoint and the second intersection point, and when the first average value difference is greater than or equal to a preset threshold, determining the forward frequency point as the maximum forward frequency point; obtaining a second average value difference between two side signals divided by the reverse frequency point between the second intersection point and the second endpoint, and when the second average value difference is greater than or equal to a preset threshold, determining the reverse frequency point as the maximum reverse frequency point.
[0086] In one embodiment, when the processor executes the computer program, the following steps are further implemented: when the first average value difference is less than the preset threshold, moving the forward frequency point in the direction of the second intersection point, and returning to execute the step of obtaining the first average value difference between two side signals divided by the forward frequency point between the first endpoint and the second intersection point; when the second average value difference is less than the preset threshold, moving the reverse frequency point in the direction of the second intersection point, and returning to execute the step of obtaining the second average value difference between two side signals divided by the reverse frequency point between the second intersection point and the second endpoint.
[0087] In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the following steps are implemented:
[0088] Obtaining a power spectral density signal corresponding to each column signal in the ultrasonic spectrogram;
[0089] Performing an integration process on the power spectral density signal to obtain a maximum forward frequency point and a maximum reverse frequency point in the column signal corresponding to the power spectral density signal;
[0090] Connecting the maximum forward frequency point and the maximum reverse frequency point corresponding to each column signal in the ultrasonic spectrogram to obtain a waveform envelope corresponding to the ultrasonic spectrogram.
[0091] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented: performing a first integration process on the power spectral density signal based on a first frequency interval corresponding to the power spectral density signal to obtain a corresponding first integration curve; determining a maximum energy point of the power spectral density signal according to the first integration curve; determining two minimum energy points of the power spectral density signal according to the maximum energy point; performing a second integration process on the power spectral density signal in a corresponding interval based on a second frequency interval corresponding to the two minimum energy points to obtain a corresponding second integration curve; determining a forward maximum frequency point and a reverse maximum frequency point in the column signal corresponding to the power spectral density signal according to the second integration curve.
[0092] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented: obtaining two end points at the head and tail of the first integration curve, and connecting the two end points by a first reference line, where the first reference line is a straight line; obtaining a first intersection point between the first reference line and the first integration curve, and taking the first intersection point as the maximum energy point of the power spectral density signal.
[0093] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented: dividing the power spectral density signal into two signal intervals according to the maximum energy point; respectively obtaining a point with the minimum ordinate in each signal interval as the minimum energy point of the power spectral density signal.
[0094] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented: obtaining a first end point at the head and a second end point at the tail of the second integration curve, and connecting the first end point and the second end point by a second reference line, where the second reference line is a straight line; determining a second intersection point, a forward frequency point and a reverse frequency point between the second reference line and the second integration curve, where the forward frequency point is the frequency point corresponding to the maximum forward distance from the second integration curve to the second reference line, and the reverse frequency point is the frequency point corresponding to the maximum reverse distance from the second integration curve to the second reference line; obtaining a first average value difference between two side signals divided by the forward frequency point between the first end point and the second intersection point, and when the first average value difference is greater than or equal to a preset threshold, determining the forward frequency point as the forward maximum frequency point; obtaining a second average value difference between two side signals divided by the reverse frequency point between the second intersection point and the second end point, and when the second average value difference is greater than or equal to a preset threshold, determining the reverse frequency point as the reverse maximum frequency point.
[0095] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented: when the first average value difference is less than a preset threshold, move the forward frequency point in the direction of the second intersection point, and return to execute the step of obtaining the first average value difference of the two-side signals divided by the forward frequency point between the first endpoint and the second intersection point; when the second average value difference is less than a preset threshold, move the reverse frequency point in the direction of the second intersection point, and return to execute the step of obtaining the second average value difference of the two-side signals divided by the reverse frequency point between the second intersection point and the second endpoint.
[0096] In one embodiment, a computer program product is provided, including a computer program, which when executed by a processor, implements the following steps:
[0097] Obtain the power spectral density signal corresponding to each column signal in the ultrasonic spectrogram;
[0098] Perform integral processing on the power spectral density signal to obtain the forward maximum frequency point and the reverse maximum frequency point in the column signal corresponding to the power spectral density signal;
[0099] Connect the forward maximum frequency point and the reverse maximum frequency point corresponding to each column signal in the ultrasonic spectrogram to obtain the waveform envelope corresponding to the ultrasonic spectrogram.
[0100] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented: perform first integral processing on the power spectral density signal based on the first frequency interval corresponding to the power spectral density signal to obtain a corresponding first integral curve; determine the maximum energy point of the power spectral density signal according to the first integral curve; determine two minimum energy points of the power spectral density signal according to the maximum energy point; perform second integral processing on the power spectral density signal in the corresponding interval based on the second frequency interval corresponding to the two minimum energy points to obtain a corresponding second integral curve; determine the forward maximum frequency point and the reverse maximum frequency point in the column signal corresponding to the power spectral density signal according to the second integral curve.
[0101] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented: obtain the two end points of the first integral curve, connect the two end points through a first reference line, and the first reference line is a straight line; obtain the first intersection point between the first reference line and the first integral curve, and use the first intersection point as the maximum energy point of the power spectral density signal.
[0102] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented: dividing the power spectral density signal into two signal intervals according to the maximum energy point; respectively obtaining the points with the minimum ordinate in each signal interval as the minimum energy points of the power spectral density signal.
[0103] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented: obtaining a first endpoint at the head and a second endpoint at the tail of the second integral curve, connecting the first endpoint and the second endpoint through a second reference line, where the second reference line is a straight line; determining a second intersection point, a forward frequency point, and a reverse frequency point between the second reference line and the second integral curve, where the forward frequency point is the frequency point corresponding to the maximum forward distance from the second integral curve to the second reference line, and the reverse frequency point is the frequency point corresponding to the maximum reverse distance from the second integral curve to the second reference line; obtaining a first average value difference between the two-side signals divided by the forward frequency point between the first endpoint and the second intersection point, and when the first average value difference is greater than or equal to a preset threshold, determining the forward frequency point as the maximum forward frequency point; obtaining a second average value difference between the two-side signals divided by the reverse frequency point between the second intersection point and the second endpoint, and when the second average value difference is greater than or equal to a preset threshold, determining the reverse frequency point as the maximum reverse frequency point.
[0104] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented: when the first average value difference is less than the preset threshold, moving the forward frequency point in the direction of the second intersection point, and returning to execute the step of obtaining the first average value difference between the two-side signals divided by the forward frequency point between the first endpoint and the second intersection point; when the second average value difference is less than the preset threshold, moving the reverse frequency point in the direction of the second intersection point, and returning to execute the step of obtaining the second average value difference between the two-side signals divided by the reverse frequency point between the second intersection point and the second endpoint.
[0105] Those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be completed by instructing relevant hardware through a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above methods. Among them, any reference to a memory, database, or other medium used in the embodiments provided in the present application can include at least one of non-volatile and volatile memories. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetoresistive random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM), etc. The databases involved in the embodiments provided in the present application can include at least one of relational databases and non-relational databases. Non-relational databases can include distributed databases based on blockchain, etc., without limitation. The processors involved in the embodiments provided in the present application can be general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic devices, data processing logics based on quantum computing, etc., without limitation.
[0106] It should be understood that the above embodiments are all exemplary and do not cover all possible implementation manners included in the claims. Without departing from the scope of the present disclosure, various deformations and changes can also be made on the basis of the above embodiments. Similarly, any combination of the technical features of the above embodiments can also be made to form additional embodiments of the present invention that may not be explicitly described. Therefore, the above embodiments only represent several implementation manners of the present invention and do not limit the protection scope of the present invention patent.
Claims
1. An ultrasonic Doppler envelope determination method, characterized in that, The method includes: Obtaining the power spectral density signal corresponding to each column signal in the ultrasonic spectrogram; Performing a first integration process on the power spectral density signal based on a first frequency interval corresponding to the power spectral density signal to obtain a corresponding first integration curve; Determining the maximum energy point of the power spectral density signal according to the first integration curve; Determining two minimum energy points of the power spectral density signal according to the maximum energy point; Performing a second integration process on the power spectral density signal in the corresponding interval based on a second frequency interval corresponding to the two minimum energy points to obtain a corresponding second integration curve; Determining a forward maximum frequency point and a reverse maximum frequency point in the column signal corresponding to the power spectral density signal according to the second integration curve; Connecting the forward maximum frequency point and the reverse maximum frequency point corresponding to each column signal in the ultrasonic spectrogram to obtain a waveform envelope corresponding to the ultrasonic spectrogram.
2. The method according to claim 1, wherein The determining the maximum energy point of the power spectral density signal according to the first integration curve includes: Obtaining the first and last endpoints of the first integration curve, and connecting the two endpoints through a first reference line, where the first reference line is a straight line; Obtaining a first intersection point between the first reference line and the first integration curve, and taking the first intersection point as the maximum energy point of the power spectral density signal.
3. The method according to claim 1, characterized in that, The determining the two minimum energy points of the power spectral density signal according to the maximum energy point includes: Dividing the power spectral density signal into two signal intervals according to the maximum energy point; Respectively obtaining the point with the minimum ordinate in each signal interval as the minimum energy point of the power spectral density signal.
4. The method according to claim 1, wherein The determining the forward maximum frequency point and the reverse maximum frequency point in the column signal corresponding to the power spectral density signal according to the second integration curve includes: Obtaining a first endpoint at the head and a second endpoint at the tail of the second integration curve, and connecting the first endpoint and the second endpoint through a second reference line, where the second reference line is a straight line; Determining a second intersection point, a forward frequency point and a reverse frequency point between the second reference line and the second integration curve, where the forward frequency point is the frequency point corresponding to the maximum positive distance from the second integration curve to the second reference line, and the reverse frequency point is the frequency point corresponding to the maximum negative distance from the second integration curve to the second reference line; Obtaining a first average value difference between the two-side signals divided by the forward frequency point between the first endpoint and the second intersection point, and when the first average value difference is greater than or equal to a preset threshold, determining the forward frequency point as the forward maximum frequency point; Obtaining a second average value difference between the two-side signals divided by the reverse frequency point between the second intersection point and the second endpoint, and when the second average value difference is greater than or equal to a preset threshold, determining the reverse frequency point as the reverse maximum frequency point.
5. The method according to claim 4, wherein The method further includes: When the first average value difference is less than a preset threshold, move the forward frequency point in the direction of the second intersection point, and return to execute the step of obtaining the first average value difference of the two-sided signals corresponding to the two sides divided by the forward frequency point between the first end point and the second intersection point; When the second average value difference is less than a preset threshold, move the reverse frequency point in the direction of the second intersection point, and return to execute the step of obtaining the second average value difference of the two-sided signals corresponding to the two sides divided by the reverse frequency point between the second intersection point and the second end point.
6. An ultrasonic Doppler envelope determination device, characterized in that, The device includes: A power spectral density signal acquisition module, configured to acquire a power spectral density signal corresponding to each column of signals in the ultrasonic spectrogram; A frequency point determination module, configured to perform a first integration process on the power spectral density signal based on a first frequency interval corresponding to the power spectral density signal to obtain a corresponding first integration curve; Determine the maximum energy point of the power spectral density signal according to the first integration curve; Determine two minimum energy points of the power spectral density signal according to the maximum energy point; Perform a second integration process on the power spectral density signal in the corresponding interval based on a second frequency interval corresponding to the two minimum energy points to obtain a corresponding second integration curve; Determine a forward maximum frequency point and a reverse maximum frequency point in the column signal corresponding to the power spectral density signal according to the second integration curve; An envelope acquisition module, configured to connect the forward maximum frequency point and the reverse maximum frequency point corresponding to each column of signals in the ultrasonic spectrogram respectively to obtain an envelope corresponding to the ultrasonic spectrogram.
7. A computer device, comprising a memory and a processor, the memory storing a computer program, characterized in that, When the processor executes the computer program, the steps of the method according to any one of claims 1 to 5 are implemented.
8. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 5 are implemented.
9. A computer program product, comprising a computer program, characterized in that, When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 5 are implemented.