Ultrasonic contrast imaging method, ultrasonic imaging device and storage medium

Abstract

The present application provides an ultrasound contrast imaging method, an ultrasound imaging device, and a storage medium. The method comprises: stimulating an ultrasound probe based on multiple ultrasound pulses to transmit ultrasound waves toward a target medium containing a contrast agent, receiving echoes of the ultrasound waves, and acquiring multiple sets of ultrasound echo data based on the echoes of the ultrasound waves; wherein at least one of the amplitude, phase, and frequency of the multiple ultrasound pulses is different; extracting contrast microbubble signals from the multiple sets of ultrasound echo data to obtain data to be processed; performing singular value decomposition filtering on the processed data to obtain processed data; and obtaining a contrast image based on the processed data. The solution of the present application performs singular value decomposition filtering on the ultrasound echo data containing microbubble signals, and processes the data used for contrast imaging in the non-frequency domain. This can avoid the problem of frequency band overlap that occurs during frequency domain processing, thereby further improving the signal-to-noise ratio and contrast-tissue residual ratio of the contrast image.

Description

[0001] manual Technical Field

[0002] The present application relates to the field of ultrasonic imaging technology, and more specifically to an ultrasonic contrast imaging method, an ultrasonic imaging device, and a storage medium. Background Art

[0003] Penetration (sensitivity) and background carryover are two key metrics for measuring contrast-enhanced ultrasound imaging performance. Further improvements in these metrics will enhance the clinical value of ultrasound imaging. Technically, penetration (sensitivity) can be abstracted as the signal-to-noise ratio (SNR). A higher SNR indicates better penetration and greater sensitivity for contrast agent detection. Background carryover is measured by the contrast-to-tissue carryover ratio (CTR). A higher CTR indicates greater suppression of background carryover. In theory, an uninjected contrast agent image contains only background noise, but this is not always the case. First, tissue itself generates nonlinear components similar to those of the contrast agent, which can be detected alongside the contrast agent's nonlinear echoes. Second, the front-end circuitry inevitably experiences amplitude and phase inconsistencies when transmitting pulse sequences, preventing complete cancellation of tissue components. Furthermore, strongly reflective surfaces such as blood vessel walls, organ capsules, and bones can easily lead to signal saturation, which also manifests as tissue carryover in the image.

[0004] In existing technologies, improving the SNR and CTR of angiographic images primarily relies on pulse sequence methods, including: forward and reverse harmonics (emitting two pulses with a 180° phase difference, then combining the echoes to extract the second harmonic component), amplitude modulation (emitting multiple pulses with different amplitudes, then combining the echoes accordingly to extract the nonlinear fundamental component), and amplitude modulation-phase reversal (emitting multiple pulses with different amplitudes and a 180° phase difference, then combining the echoes accordingly to extract both the nonlinear fundamental and second harmonic components). These methods have two drawbacks: 1) Both tissue and microbubbles generate second harmonic components, and the system noise is mostly white noise, inevitably leading to frequency band overlap between tissue residues, microbubbles, and noise. Subsequent frequency-domain filtering has very limited impact on SNR and CTR. 2) While the nonlinear fundamental can mitigate the frequency band overlap issue to some extent, improving performance compared to the second harmonic, it still cannot completely eliminate broadband residual signals caused by front-end circuit asymmetry and signal saturation. Summary of the Invention

[0005] This application provides an ultrasound contrast-enhanced imaging solution that performs singular value decomposition filtering on ultrasound echo data containing microbubble signals, processing the data used for contrast-enhanced imaging in the non-frequency domain. This improves the signal-to-noise ratio and contrast-to-tissue residual ratio of the contrast-enhanced image. The following briefly describes the ultrasound contrast-enhanced imaging solution proposed in this application, and further details will be described in the detailed description in conjunction with the accompanying drawings.

[0006] On the one hand, the present application provides an ultrasound contrast imaging method, which includes: based on multiple ultrasound pulses, exciting an ultrasound probe to transmit ultrasound to a target medium containing a contrast agent, receiving echoes of the ultrasound, and acquiring multiple groups of ultrasound echo data based on the echoes of the ultrasound; wherein, at least one of the amplitude, phase and frequency of the multiple ultrasound pulses is different; performing contrast microbubble signal extraction on the multiple groups of ultrasound echo data to obtain data to be processed; performing a first matrix processing on the data to be processed to obtain a space-time domain signal matrix; performing singular value decomposition on the space-time domain signal matrix to obtain a singular value decomposition result; performing thresholding processing on the singular value decomposition result to filter out singular values ​​corresponding to tissue residues and / or noise to obtain a processed space-time domain signal matrix; performing a second matrix processing on the processed space-time domain signal matrix to obtain processed data, and the second matrix processing is an inverse process of the first matrix processing; and obtaining a contrast image based on the processed data.

[0007] On the other hand, the present application provides an ultrasound imaging method, which includes: stimulating an ultrasound probe based on multiple ultrasound pulses to transmit ultrasound waves to a target medium containing a contrast agent, receiving echoes of the ultrasound waves, and acquiring multiple sets of ultrasound echo data based on the echoes of the ultrasound waves; wherein at least one of the amplitudes, phases, and frequencies of the multiple ultrasound pulses is different; performing contrast microbubble signal extraction on the multiple sets of ultrasound echo data to obtain data to be processed; performing singular value decomposition on the data to be processed to obtain a singular value decomposition result of a spatial domain signal matrix; performing thresholding on the singular value decomposition result to filter out singular values ​​corresponding to tissue residues and / or noise to obtain processed data; and obtaining a contrast image based on the processed data.

[0008] On the other hand, the present application provides an ultrasound imaging method, which includes: stimulating an ultrasound probe based on multiple ultrasound pulses to transmit ultrasound waves to a target medium containing a contrast agent, receiving echoes of the ultrasound waves, and acquiring multiple sets of ultrasound echo data based on the echoes of the ultrasound waves; wherein at least one of the amplitudes, phases, and frequencies of the multiple ultrasound pulses is different; performing contrast microbubble signal extraction on the multiple sets of ultrasound echo data to obtain data to be processed; performing singular value decomposition filtering on the data to be processed to obtain processed data; and obtaining a contrast image based on the processed data.

[0009] On the other hand, the present application provides an ultrasonic imaging device, which includes an ultrasonic probe, a transmit / receive sequence controller, a processor and a display, wherein: the transmit / receive sequence controller is used to: excite the ultrasonic probe based on multiple ultrasonic pulses to transmit ultrasonic waves to a target medium containing a contrast agent, receive echoes of the ultrasonic waves, and obtain multiple groups of ultrasonic echo data based on the echoes of the ultrasonic waves; wherein at least one of the amplitude, phase and frequency of the multiple ultrasonic pulses is different; the processor is used to: extract contrast microbubble signals from the multiple groups of ultrasonic echo data to obtain data to be processed; perform a first matrix processing on the processed data to obtain a space-time domain signal matrix; perform singular value decomposition on the space-time domain signal matrix to obtain a singular value decomposition result; perform thresholding on the singular value decomposition result to filter out singular values ​​corresponding to tissue residues and / or noise to obtain a processed space-time domain signal matrix; perform a second matrix processing on the processed space-time domain signal matrix to obtain processed data, and the second matrix processing is an inverse process of the first matrix processing; obtain a contrast image based on the processed data; and the display is used to: display the contrast image.

[0010] On the other hand, the present application provides an ultrasonic imaging device, which includes an ultrasonic probe, a transmit / receive sequence controller, a processor and a display, wherein: the transmit / receive sequence controller is used to: excite the ultrasonic probe based on multiple ultrasonic pulses to transmit ultrasonic waves to a target medium containing a contrast agent, receive echoes of the ultrasonic waves, and obtain multiple groups of ultrasonic echo data based on the ultrasonic echoes; wherein at least one of the amplitudes, phases and frequencies of the multiple ultrasonic pulses is different; the processor is used to: extract contrast microbubble signals from the multiple groups of ultrasonic echo data to obtain data to be processed; perform singular value decomposition on the data to be processed to obtain a singular value decomposition result of a spatial domain signal matrix; perform thresholding on the singular value decomposition result to filter out singular values ​​corresponding to tissue residues and / or noise to obtain processed data; obtain a contrast image based on the processed data; and the display is used to: display the contrast image.

[0011] On the other hand, the present application provides an ultrasonic imaging device, which includes an ultrasonic probe, a transmit / receive sequence controller, a processor and a display, wherein: the transmit / receive sequence controller is used to: excite the ultrasonic probe based on multiple ultrasonic pulses to transmit ultrasonic waves to a target medium containing a contrast agent, receive echoes of the ultrasonic waves, and obtain multiple groups of ultrasonic echo data based on the ultrasonic echoes; wherein at least one of the amplitudes, phases and frequencies of the multiple ultrasonic pulses is different; the processor is used to: extract contrast microbubble signals from the multiple groups of ultrasonic echo data to obtain data to be processed; perform singular value decomposition filtering on the data to be processed to obtain processed data; and obtain a contrast image based on the processed data; and the display is used to: display the contrast image.

[0012] In another aspect of the present application, a storage medium is provided, on which a computer program is stored. When the computer program is run, the above-mentioned ultrasound contrast imaging method is executed.

[0013] According to the ultrasound contrast imaging method, ultrasound imaging device and storage medium of the embodiments of the present application, singular value decomposition filtering is performed on the ultrasound echo data containing microbubble signals, and the data used for contrast imaging is processed in the non-frequency domain. This can avoid the frequency band overlap problem that occurs during frequency domain processing, and further improve the signal-to-noise ratio and contrast-to-tissue residual ratio of the contrast image. BRIEF DESCRIPTION OF THE DRAWINGS

[0014] Figure 1 A schematic block diagram of an exemplary ultrasound imaging device for implementing the ultrasound contrast imaging method according to an embodiment of the present application is shown.

[0015] Figure 2 A schematic flowchart of an ultrasound contrast imaging method according to an embodiment of the present application is shown.

[0016] Figure 3 An example diagram showing a contrast-enhanced image obtained without using the ultrasound contrast-enhanced imaging method according to an embodiment of the present application.

[0017] Figure 4 An example diagram of a contrast-enhanced image obtained using the ultrasound contrast-enhanced imaging method according to an embodiment of the present application is shown.

[0018] Figure 5 A schematic flowchart of an ultrasound contrast imaging method according to another embodiment of the present application is shown.

[0019] Figure 6 A schematic flowchart of an ultrasound contrast imaging method according to another embodiment of the present application is shown.

[0020] Figure 7 A schematic block diagram of an ultrasonic imaging device according to an embodiment of the present application is shown.

[0021] Figure 8 A schematic block diagram of an ultrasonic imaging device according to another embodiment of the present application is shown. DETAILED DESCRIPTION

[0022] In order to make the purpose, technical solutions and advantages of the present application more apparent, example embodiments according to the present application will be described in detail below with reference to the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments of the present application, and it should be understood that the present application is not limited to the example embodiments described herein. Based on the embodiments of the present application described in this application, all other embodiments obtained by those skilled in the art without creative work should fall within the scope of protection of this application.

[0023] In the following description, a large number of specific details are provided to provide a more thorough understanding of the present application. However, it will be apparent to those skilled in the art that the present application can be implemented without one or more of these details. In other examples, some technical features well known in the art are not described in order to avoid confusion with the present application.

[0024] It should be understood that the present application can be implemented in different forms and should not be interpreted as being limited to the embodiments set forth herein. On the contrary, providing these embodiments will make the disclosure thorough and complete and will fully convey the scope of the present application to those skilled in the art.

[0025] The purpose of the terms used herein is only to describe specific embodiments and is not intended to limit the present application. When used herein, the singular forms "a", "an", and "the" are also intended to include the plural forms, unless the context clearly indicates otherwise. It should also be understood that the terms "comprising" and / or "comprising", when used in this specification, determine the presence of the features, integers, steps, operations, elements and / or parts, but do not exclude the presence or addition of one or more other features, integers, steps, operations, elements, parts and / or groups. When used herein, the term "and / or" includes any and all combinations of the relevant listed items.

[0026] In order to thoroughly understand the present application, detailed steps and detailed structures will be provided in the following description to illustrate the technical solution proposed by the present application. The preferred embodiments of the present application are described in detail below. However, in addition to these detailed descriptions, the present application may also have other implementation methods.

[0027] First, refer to Figure 1 An exemplary ultrasound imaging device for implementing the ultrasound contrast imaging method according to an embodiment of the present application is described below.

[0028] Figure 1 FIG. 1 is a block diagram of an exemplary ultrasonic imaging device 10 for implementing the ultrasonic contrast imaging method according to an embodiment of the present application. Figure 1As shown, the ultrasonic imaging device 10 may include an ultrasonic probe 100, a transmit / receive selection switch 101, a transmit / receive sequence controller 102, a processor 103, a display 104, and a memory 105. The transmit / receive sequence controller 102 can activate the ultrasonic probe 100 to transmit ultrasonic waves toward a target object (measured object) and can also control the ultrasonic probe 100 to receive ultrasonic echoes returned from the target object, thereby obtaining ultrasonic echo signals / data. The processor 103 processes the ultrasonic echo signals / data to obtain tissue-related parameters and an ultrasonic image of the target object. The ultrasonic images obtained by the processor 103 can be stored in the memory 105 and displayed on the display 104.

[0029] In an embodiment of the present application, the display 104 of the aforementioned ultrasonic imaging device 10 may be a touch screen, a liquid crystal display, etc., or it may be an independent display device such as a liquid crystal display, a television, etc. that is independent of the ultrasonic imaging device 10, or it may be a display screen on an electronic device such as a mobile phone or a tablet computer.

[0030] In the embodiment of the present application, the memory 105 of the aforementioned ultrasonic imaging device 10 may be a flash memory card, a solid-state memory, a hard disk, etc.

[0031] An embodiment of the present application also provides a computer-readable storage medium, which stores multiple program instructions. After the multiple program instructions are called and executed by the processor 103, some or all of the steps or any combination of the steps in the ultrasound contrast imaging method in each embodiment of the present application can be executed.

[0032] In one embodiment, the computer-readable storage medium may be the memory 105 , which may be a non-volatile storage medium such as a flash memory card, a solid-state memory, or a hard disk.

[0033] In the embodiment of the present application, the processor 103 of the aforementioned ultrasound imaging device 10 can be implemented by software, hardware, firmware or a combination thereof, and can use circuits, single or multiple application-specific integrated circuits (ASICs), single or multiple general-purpose integrated circuits, single or multiple microprocessors, single or multiple programmable logic devices, or a combination of the aforementioned circuits or devices, or other suitable circuits or devices, so that the processor 103 can execute the corresponding steps of the ultrasound contrast imaging method in each embodiment.

[0034] The following combination Figures 2 to 6 The ultrasound contrast imaging method of the present application is described in detail. The method can be performed by the aforementioned ultrasound imaging device 10.

[0035] Figure 2FIG. 2 shows a schematic flow chart of an ultrasound contrast imaging method 200 according to an embodiment of the present application. Figure 2 As shown, the ultrasound contrast imaging method 200 includes the following steps:

[0036] In step S210, an ultrasonic probe is excited based on multiple ultrasonic pulses to transmit ultrasonic waves to a target medium containing a contrast agent, echoes of the ultrasonic waves are received, and multiple sets of ultrasonic echo data are acquired based on the echoes of the ultrasonic waves; wherein at least one of the amplitude, phase and frequency of the multiple ultrasonic pulses is different.

[0037] In an embodiment of the present application, an ultrasonic probe is excited to transmit ultrasonic waves to a target medium containing a contrast agent based on multiple ultrasonic pulses with different parameters, where the parameters may include at least one of amplitude, phase, and frequency. Based on this, multiple groups of ultrasonic echo data can be acquired accordingly, and different groups of ultrasonic echo data can correspond to ultrasonic pulses with different parameters. In an embodiment of the present application, the acquired ultrasonic echo data can refer to signal data of all forms and links collected by the ultrasonic imaging platform, including but not limited to: analog signals, digital signals, channel data before beamforming, data after beamforming, data before demodulation, data after demodulation, etc.

[0038] In step S220 , contrast microbubble signals are extracted from multiple sets of ultrasound echo data to obtain data to be processed.

[0039] In the embodiments of the present application, step S210 may correspond to the processing at the transmitter end in the pulse sequence processing method; accordingly, step S220 may perform the processing corresponding to the receiver end in the pulse sequence processing method. This processing is intended to extract nonlinear echo data of the contrast-enhancing microbubble signal, and is therefore referred to as contrast-enhancing microbubble signal extraction. In the embodiments of the present application, contrast-enhancing microbubble signal extraction may include, but is not limited to, processing corresponding to the receiver end, such as positive and negative harmonics, amplitude modulation, and amplitude modulation-phase inversion. Based on the above processing, the data to be processed can be obtained.

[0040] In step S230, a first matrixing process is performed on the data to be processed to obtain a space-time domain signal matrix.

[0041] In the embodiment of the present application, the data to be processed can be processed frame by frame or in batches. Figure 2 In the embodiment described, the method of processing the data to be processed in batches is mainly described, and later in combination with Figure 5In the embodiment described, the method for processing the data to be processed frame by frame is described. In an embodiment of the present application, when processing the data to be processed in batches, the data to be processed of the same number of frames can be obtained each time, or the data to be processed of different numbers of frames can be obtained. In addition, when processing the data to be processed in batches, the data to be processed obtained when two adjacent processes are processed may have data overlap (i.e., the data to be processed obtained when two adjacent processes include partially identical data), or there may be no data overlap (i.e., the data to be processed obtained when two adjacent processes do not include the same data). In other words, when processing the data to be processed in batches, the data to be processed obtained when two adjacent processes are processed may be shared or non-shared. It should be understood that better processing effects can be obtained when data is shared, because the shared data is processed more than once, and more accurate calculation results can be obtained; when data is not shared, the processing speed can be increased because the amount of data for the overall calculation is relatively small. Data sharing or non-sharing can be selected according to specific needs.

[0042] In an embodiment of the present application, when batch processing is performed on the data to be processed, multiple frames of data to be processed can be obtained each time. Since each frame of the data to be processed contains spatial domain information, and the multiple frames of data to be processed contain time domain information, the obtained multiple frames of data to be processed can be matrixed to obtain a spatial-time domain signal matrix for subsequent processing. In addition, since the subsequent processing (step S260) includes a matrixing process that is inverse to the matrixing process, in order to distinguish them from each other, the matrixing process here is referred to as the first matrixing process, and the matrixing process in the subsequent processing is referred to as the second matrixing process.

[0043] As mentioned above, the ultrasonic echo data obtained in step S210 can be various forms of signal data. The following description takes the beamformed baseband data as an example. In the embodiment of the present application, it is assumed that M frames of beamformed baseband data to be processed are obtained in step S220. The M frames of data can be recorded as: k,l (t1),s k,l (t2),…,s k,l (t M ); wherein, the index k∈[1,K] can represent the kth sampling point in the longitudinal direction, and the index l∈[1,L] can represent the lth receiving line in the transverse direction. Then, the M frames of data to be processed are subjected to the first matrix processing to obtain the space-time domain signal matrix S, which can be represented as:

[0044]

[0045] As can be seen from the above formula (1), this step transforms the K×L×M three-dimensional array (M frames of K×L-dimensional data) into a large KL×M-dimensional matrix S. The column dimension of S corresponds to the time domain information of the signal, and the row dimension corresponds to the spatial domain information of the signal, so it is defined as the space-time signal matrix.

[0046] Subsequent processing can be performed based on the obtained space-time domain signal matrix.

[0047] In step S240, singular value decomposition (SVD) is performed on the space-time domain signal matrix to obtain a singular value decomposition result.

[0048] In an embodiment of the present application, singular value decomposition is performed on the space-time domain signal matrix obtained in step S230 to obtain spatial singular vectors and temporal singular vectors containing multiple singular values. Some singular values ​​can be removed together with the corresponding singular vectors through the thresholding processing (setting thresholds based on tissue residues and / or noise) described in subsequent steps to achieve suppression of tissue residues and / or noise in the processed data. The processed data is used for contrast imaging, and the obtained contrast image can achieve improvements in SNR and CTR.

[0049] Now let's continue with the example of the previous step. For the space-time domain signal matrix S obtained in the previous step, the singular value decomposition result obtained by performing singular value decomposition on it can be expressed as:

[0050]

[0051] In the above formula (2), the operator “(·) H " represents the conjugate transpose operation of the vector / matrix; λ1>λ2>…>λ P are P singular values ​​arranged in descending order, where P is the rank of the matrix S; {u1,u2,...,u P} and {v1,v2,...,v P} are the spatial singular vectors and temporal singular vectors corresponding to the above singular values, with dimensions of KL×1 and M×1 respectively, and are orthogonal to each other.

[0052] Subsequent processing can be performed based on the obtained singular value decomposition results.

[0053] In step S250 , thresholding is performed on the singular value decomposition result to filter out singular values ​​corresponding to tissue residues and / or noise, thereby obtaining a processed space-time domain signal matrix.

[0054] As mentioned above, the singular values ​​corresponding to tissue residues and / or noise can be filtered out through thresholding processing to achieve the suppression of tissue residues and noise in the processed data. The processed data is used for contrast imaging, and the obtained contrast image can achieve the improvement of SNR and CTR.

[0055] In one embodiment of the present application, thresholding is performed on the singular value decomposition results, which may include: determining whether all singular values ​​in the singular value decomposition results are within a preset threshold range; filtering out singular values ​​and corresponding singular vectors that are not within the preset threshold range, and reconstructing all remaining singular values ​​and corresponding singular vectors to obtain a processed space-time domain signal matrix; wherein the boundary values ​​of the preset threshold range are respectively a first threshold and a second threshold, the first threshold is less than the second threshold, the first threshold is set based on the singular value corresponding to the noise, and the second threshold is set based on the singular value corresponding to the tissue residue.

[0056] In this embodiment, the singular value decomposition results are subjected to a three-stage thresholding process, i.e., the noise threshold is set as the first threshold and the tissue residue threshold is set as the second threshold. When a singular value in the singular value decomposition result is less than the first threshold, it is considered to be a singular value corresponding to noise and is discarded; when a singular value in the singular value decomposition result is greater than the second threshold, it is considered to be a singular value corresponding to tissue residue and is also discarded; when a singular value in the singular value decomposition result is within the threshold range formed by the two boundary values ​​of the first threshold and the second threshold, that is, the singular value in the singular value decomposition result is greater than or equal to the first threshold and less than or equal to the second threshold, it is considered to be a singular value corresponding to the microbubble signal and is retained, thereby achieving the suppression of tissue residue and noise in the processed data. The retained singular values ​​and corresponding singular vectors are reconstructed to obtain the processed space-time domain signal matrix. Exemplarily, the reconstruction can be linear weighted reconstruction or other reconstruction methods. It should be understood that this thresholding process is only exemplary, and other thresholding processes can also be used, such as setting only a tissue residue threshold or only a noise threshold to achieve the purpose of suppressing only tissue residues or only suppressing noise, or setting the threshold interval more finely to obtain a more precise suppression effect, etc.

[0057] Now, let's continue with the example of the previous step. The three-stage thresholding process is performed on the singular value decomposition result obtained in the previous step, which may include: 1) setting the tissue residue threshold α, when λ p When α is greater than α, it is considered as a singular value corresponding to tissue residue and is discarded; 2) Set the noise threshold β (β < α), when λ p <β is considered as a singular value corresponding to the noise and is also discarded; 3) When α≥λ p When ≥β, it is considered as the singular value corresponding to the microbubble signal and is retained. At this time, the retained singular value and its singular vector are recorded as and Perform linear weighted reconstruction on it, and the processed space-time domain signal matrix can be expressed as:

[0058]

[0059] Subsequent processing can be performed based on the obtained processed space-time domain signal matrix.

[0060] In step S260, a second matrixing process is performed on the processed space-time domain signal matrix to obtain processed data. The second matrixing process is an inverse process of the first matrixing process.

[0061] In the embodiment of the present application, the processed spatial-temporal signal matrix has been filtered to remove tissue residue and / or noise. It can now be restored to a multi-frame data format to prepare for subsequent imaging processing. As previously mentioned, to distinguish it from the first matrixing process in step S230, the matrixing process in step S260 is referred to as the second matrixing process, as the two are inverse processes.

[0062] Now let's continue with the example of the previous step to describe the processed space-time domain signal matrix obtained in the previous step. By performing matrix processing corresponding to step S230 on its column vector, the processed M frames of K×L dimensional data can be obtained, which can be expressed as:

[0063]

[0064] Subsequent imaging processing can be performed based on the obtained multi-frame processed data.

[0065] In step S270 , a contrast image is obtained based on the processed data.

[0066] In one embodiment of the present application, obtaining an angiographic image based on processed data may include: performing signal processing on the processed data to obtain grayscale image data; and generating an angiographic image based on the grayscale image data. Exemplarily, the signal processing may include envelope detection, dynamic range control, brightness compensation, and digital scan conversion (DSC). In another embodiment of the present application, after obtaining the grayscale image data, grayscale atlas mapping and / or pseudo-color atlas mapping may be performed on the grayscale image data to obtain the angiographic image.

[0067] Figure 3 and Figure 4 An example of an ultrasound image obtained without using the ultrasound contrast imaging method according to an embodiment of the present application and an example of an ultrasound image obtained by using the ultrasound contrast imaging method according to an embodiment of the present application are shown respectively. Figure 3 and Figure 4From the comparison, it can be seen that compared with the dog liver contrast imaging image obtained without using the ultrasound contrast imaging method according to the embodiment of the present application, the near-field tissue residue of the dog liver contrast imaging image obtained by using the ultrasound contrast imaging method according to the embodiment of the present application is significantly suppressed, and the effect of enhancing the microbubble signal is achieved under the premise that the background noise remains unchanged, which shows that the ultrasound contrast imaging method according to the embodiment of the present application can simultaneously improve the SNR and CTR of the contrast image.

[0068] Based on the above description, the ultrasound contrast imaging method 200 according to an embodiment of the present application performs singular value decomposition filtering on the ultrasound echo data containing microbubble signals, and processes the data used for contrast imaging in the non-frequency domain, which can avoid the frequency band overlap problem that occurs during frequency domain processing, and further improve the signal-to-noise ratio and contrast-tissue residual ratio of the contrast image.

[0069] The above is an exemplary illustration of an ultrasound contrast imaging method according to an embodiment of the present application. Figure 5 and Figure 6 The ultrasound contrast imaging method according to other embodiments of the present application is described.

[0070] Figure 5 FIG. 5 shows a schematic flow chart of an ultrasound contrast imaging method 500 according to another embodiment of the present application. Figure 5 As shown, the ultrasound contrast imaging method 500 includes the following steps:

[0071] In step S510, an ultrasonic probe is excited based on multiple ultrasonic pulses to transmit ultrasonic waves to a target medium containing a contrast agent, echoes of the ultrasonic waves are received, and multiple sets of ultrasonic echo data are acquired based on the echoes of the ultrasonic waves; wherein at least one of the amplitude, phase and frequency of the multiple ultrasonic pulses is different.

[0072] In step S520 , contrast microbubble signals are extracted from multiple sets of ultrasound echo data to obtain data to be processed.

[0073] In step S530, singular value decomposition is performed on the data to be processed to obtain a singular value decomposition result of the spatial domain signal matrix.

[0074] In step S540 , thresholding is performed on the singular value decomposition result to filter out singular values ​​corresponding to tissue residues and / or noise, thereby obtaining processed data.

[0075] In step S550 , a contrast image is obtained based on the processed data.

[0076] The ultrasound contrast imaging method 500 according to an embodiment of the present application is generally similar to the ultrasound contrast imaging method 200 according to an embodiment of the present application described above, except that the ultrasound contrast imaging method 200 according to an embodiment of the present application described above processes the data to be processed in batches, while the ultrasound contrast imaging method 500 according to an embodiment of the present application processes the data to be processed frame by frame. Due to the frame-by-frame processing, the data to be processed each time only contains spatial domain information and does not contain temporal information. Therefore, it is not necessary to perform the first and second matrixing processes described above. Instead, the processed data is directly subjected to singular value decomposition and thresholding to obtain the processed data for ultrasound imaging. In other words, the ultrasound contrast imaging method 200 according to an embodiment of the present application improves the SNR and CTR of the contrast image through spatial-temporal singular value decomposition filtering, while the ultrasound contrast imaging method 500 according to an embodiment of the present application improves the SNR and CTR of the contrast image through spatial-domain singular value decomposition filtering. In other embodiments, other singular value decomposition filters may also be used to improve the SNR and CTR of the contrast image.

[0077] Figure 6 FIG. 6 is a schematic flow chart of an ultrasound contrast imaging method 600 according to another embodiment of the present application. Figure 6 As shown, the ultrasound contrast imaging method 600 includes the following steps:

[0078] In step S610, an ultrasonic probe is excited based on multiple ultrasonic pulses to transmit ultrasonic waves to a target medium containing a contrast agent, echoes of the ultrasonic waves are received, and multiple sets of ultrasonic echo data are acquired based on the echoes of the ultrasonic waves; wherein at least one of the amplitude, phase and frequency of the multiple ultrasonic pulses is different.

[0079] In step S620 , contrast microbubble signals are extracted from multiple sets of ultrasound echo data to obtain data to be processed.

[0080] In step S630, singular value decomposition filtering is performed on the data to be processed to obtain processed data.

[0081] In step S640 , a contrast image is obtained based on the processed data.

[0082] The ultrasound contrast imaging method 600 according to an embodiment of the present application is substantially similar to the ultrasound contrast imaging methods 200 and 500 according to the embodiments of the present application described above, except that the singular value decomposition filtering of the processed data in the ultrasound contrast imaging method 600 according to the embodiment of the present application is not limited to space-time singular value decomposition filtering and space-domain singular value decomposition filtering, and can also be other singular value decomposition filters. Accordingly, the ultrasound contrast imaging methods 200 and 500 according to the embodiments of the present application described above can serve as two implementations of the ultrasound contrast imaging method 600 according to the embodiment of the present application.

[0083] In one embodiment of the present application, performing singular value decomposition filtering on the processed data in step S630 to obtain processed data may include: performing a first matrix processing on the processed data to obtain a space-time domain signal matrix; performing singular value decomposition on the space-time domain signal matrix to obtain a singular value decomposition result; performing threshold processing on the singular value decomposition result to filter out singular values ​​corresponding to tissue residues and / or noise to obtain a processed space-time domain signal matrix; performing a second matrix processing on the processed space-time domain signal matrix to obtain processed data, where the second matrix processing is an inverse process of the first matrix processing.

[0084] In another embodiment of the present application, performing singular value decomposition and filtering on the processed data in step S630 to obtain processed data may include: performing singular value decomposition on the processed data to obtain a singular value decomposition result of the spatial domain signal matrix; performing thresholding processing on the singular value decomposition result to filter out singular values ​​corresponding to tissue residues and / or noise to obtain processed data.

[0085] Based on the above description, according to the ultrasound contrast imaging methods 500 and 600 of the embodiments of the present application, singular value decomposition filtering is performed on the ultrasound echo data containing microbubble signals, and the data used for contrast imaging is processed in the non-frequency domain, which can avoid the frequency band overlap problem that occurs during frequency domain processing, and further improve the signal-to-noise ratio and contrast-tissue residual ratio of the contrast image.

[0086] The following combination Figure 7 and Figure 8 An ultrasound imaging device according to another aspect of the application is described. Figure 7 FIG. 8 shows a schematic block diagram of an ultrasonic imaging device 700 according to an embodiment of the present application. Figure 7As shown, the ultrasonic imaging device 700 may include a transmit / receive sequence controller 710, an ultrasonic probe 720, a processor 730, and a display 740. The ultrasonic imaging device 700 may be used to implement the ultrasonic contrast imaging methods 200, 500, and 600 described above according to the embodiments of the present application. Those skilled in the art may understand the specific operations of the various components of the ultrasonic imaging device 700 in conjunction with the above description. For the sake of brevity, only the main operations are briefly described here.

[0087] In an embodiment of the present application, when the ultrasonic imaging device 700 is used to implement the ultrasonic contrast imaging method 200 according to the embodiment of the present application, the transmit / receive sequence controller 710 is used to: excite the ultrasonic probe 720 based on multiple ultrasonic pulses to transmit ultrasonic waves to a target medium containing a contrast agent, receive echoes of the ultrasonic waves, and acquire multiple sets of ultrasonic echo data based on the echoes of the ultrasonic waves; wherein at least one of the amplitude, phase, and frequency of the multiple ultrasonic pulses is different; the processor 730 is used to: extract contrast microbubble signals from the multiple sets of ultrasonic echo data to obtain data to be processed; perform a first matrixing process on the data to be processed to obtain a space-time domain signal matrix; perform a singular value decomposition on the space-time domain signal matrix to obtain a singular value decomposition result; perform a thresholding process on the singular value decomposition result to filter out singular values ​​corresponding to tissue residues and / or noise to obtain a processed space-time domain signal matrix; perform a second matrixing process on the processed space-time domain signal matrix to obtain processed data, the second matrixing process being an inverse process of the first matrixing process; and obtain a contrast image based on the processed data; and the display 740 is used to: display the contrast image.

[0088] In an embodiment of the present application, when the ultrasonic imaging device 700 is used to implement the ultrasonic contrast imaging method 500 according to the embodiment of the present application, the transmit / receive sequence controller 710 is used to: stimulate the ultrasonic probe 720 based on multiple ultrasonic pulses to transmit ultrasonic waves to the target medium containing the contrast agent, receive the echo of the ultrasonic waves, and obtain multiple groups of ultrasonic echo data based on the echo of the ultrasonic waves; wherein, at least one of the amplitude, phase and frequency of the multiple ultrasonic pulses is different; the processor 730 is used to: extract contrast microbubble signals from the multiple groups of ultrasonic echo data to obtain data to be processed; perform singular value decomposition on the data to be processed to obtain the singular value decomposition results of the spatial domain signal matrix; perform thresholding on the singular value decomposition results to filter out singular values ​​corresponding to tissue residues and / or noise to obtain processed data; obtain a contrast image based on the processed data; the display 740 is used to: display the contrast image.

[0089] In an embodiment of the present application, when the ultrasonic imaging device 700 is used to implement the ultrasonic contrast imaging method 600 according to the embodiment of the present application, the transmit / receive sequence controller 710 is used to: stimulate the ultrasonic probe 720 based on multiple ultrasonic pulses to transmit ultrasonic waves to a target medium containing a contrast agent, receive the echo of the ultrasonic waves, and obtain multiple groups of ultrasonic echo data based on the echo of the ultrasonic waves; wherein at least one of the amplitude, phase and frequency of the multiple ultrasonic pulses is different; the processor 730 is used to: extract contrast microbubble signals from the multiple groups of ultrasonic echo data to obtain data to be processed; perform singular value decomposition filtering on the data to be processed to obtain processed data; obtain a contrast image based on the processed data; and the display 740 is used to: display the contrast image.

[0090] In an embodiment of the present application, the processor 730 performs thresholding processing on the singular value decomposition result, which may include: determining whether all singular values ​​in the singular value decomposition result are within a preset threshold range; filtering out singular values ​​and corresponding singular vectors that are not within the preset threshold range, and reconstructing all remaining singular values ​​and corresponding singular vectors to obtain a processed space-time domain signal matrix; wherein the boundary values ​​of the preset threshold range are a first threshold and a second threshold, respectively, the first threshold is less than the second threshold, the first threshold is set based on the singular value corresponding to the noise, and the second threshold is set based on the singular value corresponding to the tissue residue.

[0091] In an embodiment of the present application, when performing the first matrix processing on the data to be processed, the processor 730 may obtain the same number of frames of data to be processed each time to perform the first matrix processing.

[0092] In an embodiment of the present application, when performing the first matrix processing on the data to be processed, the data to be processed acquired by the processor 730 twice adjacently may include partially identical data.

[0093] In the embodiment of the present application, when performing the first matrix processing on the data to be processed, the data to be processed acquired by the processor 730 twice adjacently do not include the same data.

[0094] In the embodiment of the present application, the ultrasound echo data may include any one of the following: channel data before beamforming, data after beamforming, data before demodulation, and data after demodulation.

[0095] In an embodiment of the present application, the processor 730 obtains an angiography image based on the processed data, which may include: performing signal processing on the processed data to obtain grayscale image data; and generating an angiography image based on the grayscale image data.

[0096] In an embodiment of the present application, the processor 730 obtains an angiography image based on the processed data, and may further include: after obtaining the grayscale image data, performing grayscale atlas mapping and / or pseudo-color atlas mapping on the grayscale image data to obtain an angiography image.

[0097] In the embodiments of the present application, signal processing may include but is not limited to envelope detection, dynamic range control, brightness compensation, and digital scan conversion.

[0098] Based on the above description, the ultrasonic imaging device according to the embodiment of the present application performs singular value decomposition filtering on the ultrasonic echo data containing microbubble signals, and processes the data used for contrast imaging in the non-frequency domain, which can avoid the frequency band overlap problem that occurs during frequency domain processing, and further improve the signal-to-noise ratio and contrast-tissue residual ratio of the contrast image.

[0099] Figure 8 FIG2 is a schematic block diagram of an ultrasonic imaging device 800 according to another embodiment of the present application. The ultrasonic imaging device 800 includes a memory 810 and a processor 820 .

[0100] The memory 810 stores a program for implementing the corresponding steps of the ultrasound contrast imaging methods 200, 500, and 600 according to the embodiments of the present application. The processor 820 is configured to execute the program stored in the memory 810 to perform the corresponding steps of the ultrasound contrast imaging methods 200, 500, and 600 according to the embodiments of the present application.

[0101] In addition, according to an embodiment of the present application, a storage medium is further provided, on which program instructions are stored, and when the program instructions are executed by a computer or a processor, the corresponding steps of the ultrasound contrast imaging method of the embodiment of the present application are executed. The storage medium may include, for example, a memory card of a smart phone, a storage component of a tablet computer, a hard disk of a personal computer, a read-only memory (ROM), an erasable programmable read-only memory (EPROM), a portable compact disk read-only memory (CD-ROM), a USB memory, or any combination of the above storage media. The computer-readable storage medium may be any combination of one or more computer-readable storage media.

[0102] In addition, according to an embodiment of the present application, a computer program is also provided, which can be stored in a cloud or local storage medium. When the computer program is executed by a computer or processor, it is used to perform the corresponding steps of the ultrasound contrast imaging method of the embodiment of the present application.

[0103] Based on the above description, according to the ultrasound contrast imaging method, ultrasound imaging device and storage medium of the embodiments of the present application, singular value decomposition filtering is performed on the ultrasound echo data containing microbubble signals, and the data used for contrast imaging is processed in the non-frequency domain, which can avoid the frequency band overlap problem that occurs during frequency domain processing, and further improve the signal-to-noise ratio and contrast-tissue residual ratio of the contrast image.

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

[0105] Those skilled in the art will appreciate that the units and algorithm steps of each example 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 performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0106] In the several embodiments provided in this application, 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 example, the division of units is merely a logical function division. In actual implementation, other division methods may be used, such as combining or integrating multiple units or components into another device, or ignoring or not performing some features.

[0107] In the description provided herein, a large number of specific details are described. However, it is understood that the embodiments of the present application can be practiced without these specific details. In some instances, well-known methods, structures, and techniques are not shown in detail so as not to obscure the understanding of this description.

[0108] Similarly, it should be understood that in order to streamline the present application and aid in understanding one or more of the various inventive aspects, in the description of the exemplary embodiments of the present application, the various features of the present application are sometimes grouped together into a single embodiment, figure, or description thereof. However, this approach of the present application should not be interpreted as reflecting the intention that the application claimed for protection requires more features than those explicitly recited in each claim. More precisely, as reflected in the corresponding claims, the inventive point is that the corresponding technical problem can be solved with fewer features than all the features of a single disclosed embodiment. Therefore, the claims following the detailed description are hereby expressly incorporated into the detailed description, with each claim itself serving as a separate embodiment of the present application.

[0109] Those skilled in the art will understand that, except where mutually exclusive, all features disclosed in this specification (including the accompanying claims, abstract, and drawings) and all processes or units of any method or apparatus disclosed herein may be combined in any combination. Unless expressly stated otherwise, each feature disclosed in this specification (including the accompanying claims, abstract, and drawings) may be replaced by an alternative feature providing the same, equivalent, or similar purpose.

[0110] Furthermore, those skilled in the art will appreciate that although some embodiments herein include certain features included in other embodiments but not other features, combinations of features from different embodiments are intended to be within the scope of this application and to form different embodiments. For example, in the claims, any of the claimed embodiments may be used in any combination.

[0111] The various component embodiments of the present application can be implemented in hardware, or in a software module running on one or more processors, or in a combination thereof. Those skilled in the art will appreciate that a microprocessor or digital signal processor (DSP) can be used in practice to implement some or all of the functions of some modules according to the embodiments of the present application. The application can also be implemented as a part or all of a device program (e.g., a computer program and a computer program product) for performing the method described herein. Such a program implementing the present application can be stored on a computer-readable medium, or can have the form of one or more signals. Such a signal can be downloaded from an Internet website, or provided on a carrier signal, or provided in any other form.

[0112] It should be noted that the above embodiments illustrate rather than limit the present application, and that those skilled in the art may devise alternative embodiments without departing from the scope of the appended claims. In the claims, any reference symbols placed between brackets should not be construed as limiting the claims. The present application may be implemented by means of hardware comprising several different elements and by means of a suitably programmed computer. In a unit claim 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 may be interpreted as names.

[0113] The above is merely a description of specific embodiments of the present application, and the scope of protection of the present application is not limited thereto. Any changes or substitutions that can be easily conceived by a person skilled in the art within the technical scope disclosed in this application should be included in the scope of protection of the present application. The scope of protection of the present application should be based on the scope of protection of the claims.

Claims

1. A method for ultrasound contrast imaging, characterized in that: The method comprises: Exciting an ultrasonic probe based on a plurality of ultrasonic pulses to transmit ultrasonic waves toward a target medium containing a contrast agent, receiving echoes of the ultrasonic waves, and acquiring a plurality of sets of ultrasonic echo data based on the echoes of the ultrasonic waves; wherein at least one of amplitude, phase, and frequency of the plurality of ultrasonic pulses is different; Performing contrast microbubble signal extraction on the multiple sets of ultrasound echo data to obtain data to be processed, wherein the contrast microbubble signal extraction includes: extracting nonlinear echo data of the contrast microbubble signal; Performing a first matrix processing on the data to be processed to obtain a space-time domain signal matrix; Performing singular value decomposition on the space-time domain signal matrix to obtain a singular value decomposition result; performing thresholding processing on the singular value decomposition result to filter out singular values ​​corresponding to tissue residues and / or noise, thereby obtaining a processed space-time domain signal matrix; performing a second matrix processing on the processed space-time domain signal matrix to obtain processed data, wherein the second matrix processing and the first matrix processing are inverse processes of each other; A contrast image is obtained based on the processed data.

2. The method according to claim 1, characterized in that The thresholding process of the singular value decomposition result includes: Determine whether all singular values ​​in the singular value decomposition result are within a preset threshold range; Filtering out singular values ​​and corresponding singular vectors that are not within the preset threshold range, and reconstructing all remaining singular values ​​and corresponding singular vectors to obtain the processed space-time domain signal matrix; The boundary values ​​of the preset threshold range are respectively a first threshold and a second threshold, the first threshold is smaller than the second threshold, the first threshold is set based on the singular value corresponding to the noise, and the second threshold is set based on the singular value corresponding to the tissue residue.

3. The method according to claim 1 or 2, characterized in that When performing the first matrix processing on the data to be processed, the first matrix processing is performed on the data to be processed of the same number of frames each time.

4. The method according to claim 3, characterized in that When the first matrix processing is performed on the data to be processed, the data to be processed acquired twice adjacently include some identical data.

5. The method according to claim 3, characterized in that When the first matrix processing is performed on the data to be processed, the data to be processed acquired twice adjacently do not include the same data.

6. The method according to any one of claims 1 to 5, characterized in that The ultrasonic echo data includes any one of the following: channel data before beamforming, data after beamforming, data before demodulation, and data after demodulation.

7. The method according to any one of claims 1 to 6, characterized in that The performing ultrasound contrast imaging based on the processed data includes: performing signal processing on the processed data to obtain grayscale image data; A contrast image is generated based on the grayscale image data.

8. The method according to claim 7, characterized in that The method further comprises: After the grayscale image data is obtained, grayscale atlas mapping and / or pseudo-color atlas mapping is performed on the grayscale image data to obtain an angiographic image.

9. The method according to claim 7 or 8, characterized in that The signal processing includes envelope detection, dynamic range control, brightness compensation and digital scan conversion.

10. A method for ultrasound contrast imaging, characterized in that: The method comprises: Exciting an ultrasonic probe based on a plurality of ultrasonic pulses to transmit ultrasonic waves toward a target medium containing a contrast agent, receiving echoes of the ultrasonic waves, and acquiring a plurality of sets of ultrasonic echo data based on the echoes of the ultrasonic waves; wherein at least one of amplitude, phase, and frequency of the plurality of ultrasonic pulses is different; Performing contrast microbubble signal extraction on the multiple sets of ultrasound echo data to obtain data to be processed, wherein the contrast microbubble signal extraction includes: extracting nonlinear echo data of the contrast microbubble signal; Performing singular value decomposition on the data to be processed to obtain a singular value decomposition result of a spatial domain signal matrix; performing thresholding processing on the singular value decomposition result to filter out singular values ​​corresponding to tissue residues and / or noise to obtain processed data; A contrast image is obtained based on the processed data.

11. The method according to claim 10, characterized in that The thresholding process of the singular value decomposition result includes: Determine whether all singular values ​​in the singular value decomposition result are within a preset threshold range; filtering out singular values ​​and corresponding singular vectors that are not within the preset threshold range, and reconstructing all remaining singular values ​​and corresponding singular vectors to obtain the processed data; The boundary values ​​of the preset threshold range are respectively a first threshold and a second threshold, the first threshold is smaller than the second threshold, the first threshold is set based on the singular value corresponding to the noise, and the second threshold is set based on the singular value corresponding to the tissue residue.

12. An ultrasonic imaging device, characterized in that: The device includes an ultrasound probe, a transmit / receive sequence controller, a processor, and a display, wherein: The transmit / receive sequence controller is used for: The ultrasonic probe is stimulated to transmit ultrasonic waves toward a target medium containing a contrast agent based on a plurality of ultrasonic pulses, echoes of the ultrasonic waves are received, and a plurality of sets of ultrasonic echo data are acquired based on the echoes of the ultrasonic waves; wherein at least one of amplitude, phase, and frequency of the plurality of ultrasonic pulses is different; The processor is configured to: Performing contrast microbubble signal extraction on the multiple sets of ultrasound echo data to obtain data to be processed, wherein the contrast microbubble signal extraction includes: extracting nonlinear echo data of the contrast microbubble signal; Performing a first matrix processing on the data to be processed to obtain a space-time domain signal matrix; Performing singular value decomposition on the space-time domain signal matrix to obtain a singular value decomposition result; performing thresholding processing on the singular value decomposition result to filter out singular values ​​corresponding to tissue residues and / or noise, thereby obtaining a processed space-time domain signal matrix; performing a second matrix processing on the processed space-time domain signal matrix to obtain processed data, wherein the second matrix processing and the first matrix processing are inverse processes of each other; obtaining an angiographic image based on the processed data; The display is used to display the angiography image.

13. The device according to claim 12, characterized in that The processor performs thresholding processing on the singular value decomposition result, comprising: Determine whether all singular values ​​in the singular value decomposition result are within a preset threshold range; Filtering out singular values ​​and corresponding singular vectors that are not within the preset threshold range, and reconstructing all remaining singular values ​​and corresponding singular vectors to obtain the processed space-time domain signal matrix; The boundary values ​​of the preset threshold range are respectively a first threshold and a second threshold, the first threshold is smaller than the second threshold, the first threshold is set based on the singular value corresponding to the noise, and the second threshold is set based on the singular value corresponding to the tissue residue.

14. The device according to claim 12 or 13, characterized in that When performing the first matrix processing on the data to be processed, the processor obtains the same number of frames of the data to be processed each time to perform the first matrix processing.

15. The device according to claim 14, characterized in that When performing the first matrix processing on the data to be processed, the data to be processed acquired by the processor twice consecutively include some identical data.

16. The device according to claim 14, characterized in that When performing the first matrix processing on the data to be processed, the data to be processed acquired by the processor twice consecutively do not include the same data.

17. The device according to any one of claims 12 to 16, characterized in that The ultrasonic echo data includes any one of the following: channel data before beamforming, data after beamforming, data before demodulation, and data after demodulation.

18. The device according to any one of claims 13 to 17, characterized in that The processor performs ultrasound contrast imaging based on the processed data, including: performing signal processing on the processed data to obtain grayscale image data; A contrast image is generated based on the grayscale image data.

19. The device according to claim 18, characterized in that The processor performs ultrasound contrast imaging based on the processed data, further comprising: After the grayscale image data is obtained, grayscale atlas mapping and / or pseudo-color atlas mapping is performed on the grayscale image data to obtain an angiographic image.

20. The device according to claim 18 or 19, characterized in that The signal processing includes envelope detection, dynamic range control, brightness compensation and digital scan conversion.

21. An ultrasonic imaging device, characterized in that: The device includes an ultrasound probe, a transmit / receive sequence controller, a processor, and a display, wherein: The transmit / receive sequence controller is used for: The ultrasonic probe is stimulated to transmit ultrasonic waves toward a target medium containing a contrast agent based on a plurality of ultrasonic pulses, echoes of the ultrasonic waves are received, and a plurality of sets of ultrasonic echo data are acquired based on the echoes of the ultrasonic waves; wherein at least one of amplitude, phase, and frequency of the plurality of ultrasonic pulses is different; The processor is configured to: Performing contrast microbubble signal extraction on the multiple sets of ultrasound echo data to obtain data to be processed, wherein the contrast microbubble signal extraction includes: extracting nonlinear echo data of the contrast microbubble signal; Performing singular value decomposition on the data to be processed to obtain a singular value decomposition result of a spatial domain signal matrix; performing thresholding processing on the singular value decomposition result to filter out singular values ​​corresponding to tissue residues and / or noise to obtain processed data; obtaining an angiographic image based on the processed data; The display is used to display the angiography image.

22. The device according to claim 21, characterized in that The processor performs thresholding processing on the singular value decomposition result, comprising: Determine whether all singular values ​​in the singular value decomposition result are within a preset threshold range; filtering out singular values ​​and corresponding singular vectors that are not within the preset threshold range, and reconstructing all remaining singular values ​​and corresponding singular vectors to obtain the processed data; The boundary values ​​of the preset threshold range are respectively a first threshold and a second threshold, the first threshold is smaller than the second threshold, the first threshold is set based on the singular value corresponding to the noise, and the second threshold is set based on the singular value corresponding to the tissue residue.

23. A storage medium, characterized in that The storage medium stores a computer program, which, when run, executes the ultrasound contrast imaging method according to any one of claims 1 to 11.

Images