WO2021232192A1 - Procédé et appareil d'imagerie de contraste ultrasonore, et support de stockage - Google Patents

Procédé et appareil d'imagerie de contraste ultrasonore, et support de stockage Download PDF

Info

Publication number
WO2021232192A1
WO2021232192A1 PCT/CN2020/090855 CN2020090855W WO2021232192A1 WO 2021232192 A1 WO2021232192 A1 WO 2021232192A1 CN 2020090855 W CN2020090855 W CN 2020090855W WO 2021232192 A1 WO2021232192 A1 WO 2021232192A1
Authority
WO
WIPO (PCT)
Prior art keywords
data
contrast
frame
difference
adjacent frames
Prior art date
Application number
PCT/CN2020/090855
Other languages
English (en)
Chinese (zh)
Inventor
章希睿
桑茂栋
朱磊
何绪金
Original Assignee
深圳迈瑞生物医疗电子股份有限公司
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 深圳迈瑞生物医疗电子股份有限公司 filed Critical 深圳迈瑞生物医疗电子股份有限公司
Priority to CN202080006417.2A priority Critical patent/CN113164160B/zh
Priority to PCT/CN2020/090855 priority patent/WO2021232192A1/fr
Publication of WO2021232192A1 publication Critical patent/WO2021232192A1/fr

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques
    • A61B8/481—Diagnostic techniques involving the use of contrast agent, e.g. microbubbles introduced into the bloodstream
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts

Definitions

  • This application relates to the field of ultrasound contrast imaging technology, and more specifically to an ultrasound contrast imaging method, device and storage medium.
  • tissue residue there are three factors that can cause tissue residue.
  • One is that the tissue itself will also produce a nonlinear component, which is detected together with the nonlinear echo of the contrast agent; the other is that the linear component of the tissue cannot be cancelled after the amplitude of each pulse in the nonlinear imaging method is modulated due to the signal saturation. , Which mostly appear on strong reflective surfaces such as blood vessel walls, organ capsules, and bones; third, movement during scanning (probe displacement, patient respiration, etc.), which will cause the loss of correlation between pulse echoes in nonlinear imaging methods. Furthermore, the condition of tissue signal cancellation is not satisfied.
  • Tissue residues can interfere with the user's observation of the contrast agent signal, thereby affecting the accuracy of the diagnosis result.
  • the existing solutions start from the front-end hardware on the one hand, and study how to further improve the symmetry of the transmitted waveform, but the improvement is extremely difficult; on the other hand, it uses new imaging methods (sub-harmonics, etc.) or signal processing methods (filtering). Etc.)
  • the tissue residue is suppressed as much as possible, but in addition to the non-linear components generated by the tissue itself, the tissue residue and the contrast agent signal caused by other reasons cannot be distinguished in the frequency domain, and the suppression ability is limited.
  • the present application provides a contrast-enhanced ultrasound imaging solution, which can effectively suppress tissue residue in the contrast-enhanced data.
  • a contrast-enhanced ultrasound imaging method includes: controlling an ultrasound probe to transmit ultrasound to a target object, receiving an echo of the ultrasound, and obtaining an ultrasound echo signal based on the echo of the ultrasound; Acquire at least two frames of contrast signal data based on the ultrasound echo signal, wherein each frame of contrast signal data is used to generate a frame of contrast image; calculate the information difference between the contrast signal data of two adjacent frames, and calculate the information difference
  • the relationship with one frame of contrast signal data in the two adjacent frames of contrast signal data based on the relationship, the one frame of contrast signal data is divided into regions, and the contrast signal data in different regions is differently suppressed Factor processing; and generating a contrast image based on the processed one frame of contrast signal data.
  • the method includes: controlling an ultrasound probe to transmit ultrasound to a target object, receiving echoes of the ultrasound, and obtaining ultrasound echo signals based on the echo of the ultrasound.
  • Acquire contrast signal data based on the ultrasound echo signal and generate contrast images based on the contrast signal data; acquire at least two frames of contrast images, calculate the information difference between two adjacent frames of contrast images, and calculate the information difference And the relationship with one of the two adjacent frames of contrast images; and based on the relationship, the one frame of contrast images is divided into regions, and different regions of the one frame of contrast images are used.
  • Inhibitors are processed.
  • a contrast-enhanced ultrasound imaging device comprising: an ultrasound probe, a transmission/reception sequence controller, and a processor, wherein the transmission/reception sequence controller is used to: excite the ultrasound probe Transmit ultrasonic waves to a target object, receive echoes of the ultrasonic waves, and obtain ultrasonic echo signals based on the echoes of the ultrasonic waves; the processor is used to execute the above-mentioned contrast-enhanced ultrasound imaging method.
  • a contrast-enhanced ultrasound imaging method includes: acquiring at least two frames of contrast data; calculating the information difference between two adjacent frames of contrast data, and calculating the difference between the information and the phase.
  • the relationship between one frame of contrast data in two adjacent frames of contrast data classifying the one frame of contrast data based on the relationship, and suppressing at least part of the one frame of contrast data based on the result of the classification deal with.
  • a contrast-enhanced ultrasound imaging method includes: controlling an ultrasound probe to transmit ultrasound to a target object, receiving echoes of the ultrasound, and obtaining ultrasound echo signals based on the echo of the ultrasound.
  • At least part of the contrast signal data is subjected to suppression processing; and a contrast image is generated based on the processed one frame of contrast signal data.
  • a contrast-enhanced ultrasound imaging method includes: controlling an ultrasound probe to transmit ultrasound to a target object, receiving echo of the ultrasound, and obtaining ultrasound echo signals based on the echo of the ultrasound.
  • a contrast-enhanced ultrasound imaging device which includes an ultrasound probe, a transmission/reception sequence controller, and a processor, wherein: the transmission/reception sequence controller is used to: The target object emits ultrasonic waves, receives the echoes of the ultrasonic waves, and obtains ultrasonic echo signals based on the echoes of the ultrasonic waves; the processor is configured to: obtain at least two frames of contrast data based on the ultrasonic echo signals; The information difference between the contrast data of the frame is calculated, and the relationship between the information difference and one frame of the contrast data of the two adjacent frames of contrast data is calculated; the one frame of contrast data is classified based on the relationship, and Based on the result of the classification, at least part of the data of the one frame of contrast data is suppressed.
  • a contrast-enhanced ultrasound imaging device which includes an ultrasound probe, a transmission/reception sequence controller, and a processor, wherein: the transmission/reception sequence controller is used to: The target object emits ultrasonic waves, receives the echoes of the ultrasonic waves, and obtains ultrasonic echo signals based on the echoes of the ultrasonic waves; the processor is used to execute the above-mentioned contrast-enhanced ultrasound imaging method.
  • a storage medium is provided, and a computer program is stored on the storage medium, and the computer program executes the above-mentioned ultrasound contrast imaging method during operation.
  • the contrast-enhanced ultrasound imaging method, device, and storage medium excavate the difference between the residual tissue and the contrast agent signal from the perspective of time domain difference (inter-frame difference), and can compare the residual part of the tissue and the contrast agent signal in the contrast data.
  • the parts are effectively distinguished, so that the residual part of the tissue can be effectively suppressed and the effect of the contrast agent strength is not affected as much as possible.
  • Fig. 1 shows a schematic diagram of residual tissue in a contrast image.
  • Fig. 2 shows a schematic block diagram of an ultrasound contrast imaging apparatus according to an embodiment of the present application.
  • Fig. 3 shows a schematic flowchart of a contrast-enhanced ultrasound imaging method according to an embodiment of the present application.
  • Fig. 4 shows double real-time images of human thyroid contrast without injection of contrast agent.
  • FIG. 5A shows the comparison result of the signal amplitude of the first category data points in one frame of contrast data before and after the inter-frame difference.
  • FIG. 5B shows the comparison result of the signal amplitude of the second type of data points in one frame of contrast data before and after the inter-frame difference.
  • FIG. 5C shows the comparison result of the signal amplitude of the third category data points in one frame of contrast data before and after the inter-frame difference.
  • FIG. 6A shows an example of a double real-time image of human thyroid contrast when no contrast agent is injected and is not processed by the ultrasound contrast imaging method according to an embodiment of the present application.
  • FIG. 6B shows an example of a double real-time image of human thyroid contrast after being processed by the ultrasound contrast imaging method according to the embodiment of the present application when no contrast agent is injected.
  • FIG. 7A shows an example of a canine liver contrast double real-time image after injection of a contrast agent and not processed by the ultrasound contrast imaging method according to an embodiment of the present application.
  • FIG. 7B shows an example of a double real-time image of canine liver contrast after injection of a contrast agent and processed by the ultrasound contrast imaging method according to an embodiment of the present application.
  • Fig. 8 shows a schematic flowchart of a contrast-enhanced ultrasound imaging method according to another embodiment of the present application.
  • Fig. 9 shows a schematic flowchart of a contrast-enhanced ultrasound imaging method according to still another embodiment of the present application.
  • Fig. 10 shows a schematic flowchart of a contrast-enhanced ultrasound imaging method according to another embodiment of the present application.
  • Fig. 11 shows a schematic flow chart of a contrast-enhanced ultrasound imaging method according to still another embodiment of the present application.
  • Fig. 12 shows a schematic block diagram of an ultrasound contrast imaging apparatus according to another embodiment of the present application.
  • the present application provides a contrast-enhanced ultrasound imaging solution, which excavates the difference between the tissue residue and the contrast agent signal from the perspective of time domain difference (inter-frame difference), which can effectively suppress the tissue residue in the contrast data.
  • time domain difference inter-frame difference
  • Fig. 2 is a schematic diagram of an ultrasound contrast imaging device in an embodiment of the application.
  • the ultrasound contrast imaging apparatus 100 may include a probe 110, a transmitting circuit 112, a receiving circuit 114, a beam combining circuit 116, a processor 118, and a display 120.
  • the transmitting circuit 112 and the receiving circuit 114 can be connected to the ultrasonic probe 110 through the transmitting/receiving selection switch 122.
  • the transmitting circuit 112 can excite the probe 110 to transmit ultrasonic waves to the target tissue containing the contrast agent; the receiving circuit 114 can receive the ultrasonic echo returned from the target tissue through the probe 110 to obtain the ultrasonic echo signal/data; the ultrasonic echo signal/data After the beam combining circuit 116 performs beam combining processing, it is sent to the processor 118.
  • the processor 118 processes the ultrasound echo signal/data to obtain an ultrasound image and a contrast image of the target tissue.
  • the ultrasound image and the contrast image obtained by the processor 118 may be stored in the memory 124. These ultrasound images and contrast images can be displayed on the display 120.
  • contrast-enhanced ultrasound imaging apparatus 100 shown in FIG. 2 are only schematic, and may include more or fewer components. The present invention is not limited to this.
  • the contrast-enhanced ultrasound imaging apparatus 100 shown in FIG. 2 may be used to perform the following contrast-enhanced ultrasound imaging method according to the embodiment of the present invention.
  • FIG. 3 shows a contrast-enhanced ultrasound imaging method 300 according to an embodiment of the present application.
  • the contrast-enhanced ultrasound imaging method 300 may include the following steps:
  • step S310 at least two frames of contrast data are acquired.
  • step S320 the information difference between two adjacent frames of contrast data is calculated, and the relationship between the information difference and one frame of contrast data in the two adjacent frames of contrast data is calculated.
  • step S330 one frame of contrast data is classified based on the relationship, and at least part of the one frame of contrast data is suppressed based on the classification result.
  • the at least two frames of contrast data acquired in step S310 may be at least two frames of contrast signal data, wherein each frame of contrast signal data is used to generate one frame of contrast image; in addition, the contrast image acquired in step S310
  • the at least two frames of contrast data may also be at least two frames of contrast image data.
  • the at least two frames of contrast data obtained in step S310 may be generated in real time, or may be from any source.
  • the data with strong tissue residue in one frame of contrast data can be Points and data points with weak tissue residue (or even no tissue residue) are clearly distinguished (described in further detail by examples below), so that a frame of angiographic data can be classified based on this relationship (for example, classified as data with strong tissue residue) Points, data points with weak tissue residues, data points without tissue residues, etc.), and based on the classification results, at least part of the data of a frame of contrast data (such as data points with strong tissue residue and data points with weak tissue residue) are suppressed Processing, so as to achieve the effect of effectively suppressing the residual part of the tissue without affecting the intensity of the contrast agent as much as possible.
  • the contrast data of two adjacent frames includes the contrast data of the previous frame (for example, called the nth frame of contrast data D n , where n is a natural number) and the contrast data of the next frame (for example, called the nth + One frame of contrast data D n+1 .
  • calculating the information difference between two adjacent frames of contrast data in step S320 may include: calculating the contrast data of the next frame and the contrast data of the previous frame The result of the forward difference is the result of D n+1- D n .
  • the relationship between the calculated information difference in step S320 and one frame of contrast data in two adjacent frames of contrast data may include: before calculation The data relationship between the difference result (ie D n+1- D n ) and the contrast data of the previous frame (ie D n ), that is, the data relationship between (D n+1- D n ) and D n is calculated.
  • calculating the information difference between two adjacent frames of contrast data in step S320 may include: calculating the backward difference result between the contrast data of the previous frame and the contrast data of the next frame, that is, D n- D n+ results 1.
  • the information calculated in step S320 and the relationship between the contrast difference data of the contrast data of two adjacent one may include: calculating the result of the difference (i.e., D n -D n + 1)
  • the data relationship with the next frame of contrast data ie D n+1
  • D n+1 the data relationship between (D n- D n+1 ) and D n+1 is calculated. Therefore, in step S320, the difference between the calculated information and the adjacent
  • the relationship between one frame of contrast data in the two frames of contrast data can also be described as: the relationship between the difference in calculation information and the difference of one frame of contrast data in two adjacent frames of contrast data.
  • one frame of contrast data when it is the aforementioned D n or D n+1 , can be a complex analysis containing information such as amplitude and phase.
  • the signal is, for example, an I ⁇ J-dimensional complex matrix, where I and J can represent the number of vertical dots and the number of horizontal lines, respectively.
  • the aforementioned information difference between two adjacent frames of contrast data may include amplitude difference, phase difference or frequency difference.
  • one frame of contrast data (whether it is the aforementioned D n or D n+1 ) can be an analytical signal containing at least amplitude information, for example, an I ⁇ J dimensional matrix, Among them, I and J can respectively represent the number of vertical dots and the number of horizontal lines. Based on this, the aforementioned difference in information between two adjacent frames of contrast data may include a difference in amplitude.
  • step S320 when the information difference between two adjacent frames of contrast data is a magnitude difference, in one embodiment, in step S320, the information difference between two adjacent frames of contrast data is calculated, and the information difference is calculated as compared with the adjacent one.
  • the relationship between one frame of contrast data in two frames of contrast data may include: calculating the difference data between two adjacent frames of contrast data, and calculating the respective amplitude data of the difference data and one frame of contrast data; the amplitude based on the difference data
  • the data and the amplitude data of one frame of contrast data are calculated to calculate the amplitude ratio of the difference data to one frame of contrast data, and the result of the relationship between the difference data and one frame of contrast data is obtained.
  • the ratio may be calculated after logarithmic compression is performed on the amplitude data to reduce the amount of calculation. For example, if the ratio (result of the relationship) is expressed as R, then R can be expressed as the following formula:
  • abs represents the amplitude (absolute value)
  • log10 represents the logarithmic operation with 10 as the base
  • symbol “./” represents the point-to-point division operation of the matrix
  • R is the I ⁇ J-dimensional real matrix
  • I and J are such as The foregoing. It should be understood that this formula is only exemplary.
  • the information difference between the contrast data of two adjacent frames is the amplitude difference
  • the information difference between the contrast data of the two adjacent frames is calculated in step S320, and the information difference is calculated as compared with that of the adjacent two frames.
  • the relationship between one frame of contrast data in the contrast data may include: calculating the difference data between two adjacent frames of contrast data, calculating the ratio between the difference data and one frame of contrast signal data; solving the amplitude of the contrast value to obtain the difference data and The amplitude ratio of one frame of contrast data is used as the result of the relationship between the difference data and one frame of contrast data.
  • the ratio between the difference data and the differenced frame of contrast data is calculated first, and then the amplitude is solved for the ratio, and finally the amplitude ratio between the difference data and the differenced frame of contrast data can be obtained.
  • the difference between the difference data and the differenced frame of contrast signal data can also be calculated, and then based on the difference and The relationship of the preset threshold value determines the result of the relationship between the difference data and the differenced one frame of contrast data.
  • step S320 when the information difference between two adjacent frames of contrast data is a phase difference, in one embodiment, in step S320, the information difference between two adjacent frames of contrast data is calculated and calculated
  • the relationship between the information difference and one frame of contrast data in two adjacent frames of contrast data may include: calculating the difference data between two adjacent frames of contrast data, and calculating the respective phase data of the difference data and one frame of contrast data; Based on the phase data of the differential data and the phase data of one frame of contrast data, the phase difference between the differential data and one frame of contrast data is calculated, and the result of the relationship between the differential data and one frame of contrast data is obtained.
  • first calculate the phase data of the difference data and the differenced frame of contrast data then calculate the difference between the two, and use the calculated phase difference value as the difference data and the differenced frame of contrast data. The result of the relationship between the data.
  • the information difference between two adjacent frames of contrast data is a phase difference
  • the ratio between the difference data and the differenced frame of contrast data is calculated first, and then the phase is solved for the ratio, and finally the phase difference between the difference data and the differenced frame of contrast data can be obtained.
  • Value as the result of the relationship between the difference data and the differenced one frame of contrast data.
  • one frame of contrast data when the contrast data is contrast image data, one frame of contrast data (whether it is the aforementioned D n or D n+1 ) may be a data set containing gray information, for example, V ⁇ H A dimensional matrix, where V and H can represent the number of vertical pixels and the number of horizontal pixels, respectively.
  • calculating the information difference between two adjacent frames of contrast data in step S320, and calculating the relationship between the information difference and one frame of contrast data in the adjacent two frames of contrast data may include : Calculate the forward difference result between the next frame of contrast image and the previous frame; calculate the gray ratio between the forward difference result and the previous frame of contrast image.
  • calculating the information difference between two adjacent frames of contrast data and calculating the relationship between the information difference and one frame of contrast data in the two adjacent frames of contrast data in step S320 may include: calculating The backward difference result of the previous frame of contrast image and the next frame of contrast image; calculate the gray ratio of the backward difference result and the next frame of contrast image.
  • the above exemplarily describes the calculation method of the relationship between the information difference between two adjacent frames of contrast data and one frame of contrast data in the adjacent two frames of contrast data.
  • the relationship between the information difference between two frames of angiographic data and the differenced one frame of angiographic data in two adjacent frames of angiographic data can reduce the data points with strong tissue residue and weak (or even no tissue residue) in one frame of angiographic data.
  • the data points of tissue residues are clearly distinguished. This will be described below in conjunction with FIG. 4 and FIGS. 5A to 5C.
  • Fig. 4 shows a double real-time image of human thyroid angiography when no contrast agent is injected. As shown in Figure 4, it shows that a total of 100 frames of contrast beam synthesis of the human thyroid gland are collected under the condition of no contrast agent.
  • One of the corresponding contrast double real-time images of the I/Q data is shown in Figure 4 The left half of is the tissue image, and the right half is the contrast image. As shown in the right half of FIG. 4, the tissue residues are mainly concentrated at the data points included in the boxes 410 and 420, and there is only noise at the remaining data points.
  • a frame of contrast data (contrast signal data or contrast image data) can be classified.
  • the data points included in the box 410 are mainly represented as cortex that is likely to cause signal saturation.
  • the layer which is the most obvious data point location area of tissue residue, is generally the near field area, and the relative displacement with the probe is small, and it is far away from the organ to be scanned, and is not easily affected by the patient’s breathing; after the frame difference, the data point The signal amplitude in the location area will be significantly reduced, as shown in Figure 5A.
  • the data points included in the frame 420 (for example, called the second category data points) mainly include tissue residues caused by strong reflective surfaces such as blood vessel walls and organ envelopes.
  • the data points included in the box 430 (for example, referred to as the third category data points) mainly include inter-frame noise, which is generally Gaussian-white noise, and generally the far-field area; after the inter-frame difference, the overall variance will become larger. , As shown in Figure 5C.
  • the embodiments of the present application are based on this, and the strong residual tissue in a frame of angiography data can be Data points and data points with weak tissue residue (or even no tissue residue) are clearly distinguished, so that a frame of contrast data can be classified based on the relationship, and at least part of the frame of contrast data can be suppressed based on the classification result Therefore, the residual part of the tissue can be effectively suppressed and the effect of the intensity of the contrast agent is not affected as much as possible.
  • one frame of contrast data can be expressed as a matrix, so the difference data of two adjacent frames of contrast data can also be expressed as a matrix.
  • the difference data of two adjacent frames of contrast data can be expressed as the first A matrix
  • a frame of contrast data to be compared with the first matrix is expressed as a second matrix
  • the step S330 to classify a frame of contrast data based on the relationship may include: determining the data points in the first matrix and The numerical interval to which the relationship result of the data point at the corresponding position in the second matrix belongs; and the data point at the corresponding position in the second matrix is classified as belonging to the corresponding category based on the numerical interval to which the relationship result belongs.
  • data points of different categories belonging to different numerical intervals in the second matrix can be divided by presetting a number of numerical intervals.
  • the numerical interval may be determined based on one or more preset thresholds. Obviously, the greater the number of preset thresholds, the more numerical intervals are divided, and the finer the classification of data points in a frame of angiographic data.
  • a numerical interval based on a preset threshold T (where the value range of T is, for example, 0 to 1) is used as an example for description.
  • the preset numerical interval may include a first numerical interval and a second numerical interval, the first numerical interval is, for example, (0, T], and the second numerical interval is, for example, (T, 1). (0,T) corresponds to the first type of data point, and the second numerical interval (T,1) corresponds to the second type of data point.
  • the data point in the first matrix (the difference data of two adjacent frames of angiographic data) and the first The relationship result (for example, the aforementioned amplitude ratio R) of the data points at the corresponding positions in the two matrices (one frame of contrast data that is differenced in two adjacent frames of contrast data) is within the first numerical interval, that is, R ⁇ (0, T], the data points at the corresponding positions in the second matrix are divided into data points of the first category; when determining the relationship between the data points in the first matrix and the data points at the corresponding positions in the second matrix (for example, When the aforementioned amplitude ratio R) is within the second numerical interval, that is, R ⁇ (T,1], the data point at the corresponding position in the second matrix is divided into the second type of data point.
  • the signal amplitude will be significantly reduced compared to the signal amplitude of the differenced frame of contrast data.
  • the ratio between the amplitudes is relatively small, that is, the R value is relatively small; as described above in conjunction with Figure 4 and Figure 5B, after the inter-frame difference, the signal amplitude is slightly weaker or close to the signal amplitude of the differenced frame of contrast data , That is, the ratio between the signal amplitude of the differential data and the signal amplitude of one frame of contrast data to be differentiated is relatively large, that is, the R value is relatively large.
  • the first category data points corresponding to R ⁇ (0, T) are compared with the second category data points corresponding to R ⁇ (T, 1], and the tissue residue at the first category data points is stronger than (more than ) Tissue residues at data points of the second category.
  • the numerical interval based on a preset threshold T may also include a third numerical interval, that is, a numerical interval (1,+ ⁇ ) ⁇ (- ⁇ ,0).
  • a third numerical interval that is, a numerical interval (1,+ ⁇ ) ⁇ (- ⁇ ,0).
  • the signal amplitude is stronger than the signal amplitude of the differenced one frame of contrast data, that is, R is greater than 1, which means the third
  • the category data points are noise data points (R less than 0 is considered to be a manifestation of the randomness of noise).
  • performing suppression processing on at least part of the data in one frame of contrast data based on the classification result in step S330 may include: determining, based on the classification result, data that needs to be suppressed in one frame of contrast data; Different types of data in the data that need to be suppressed are processed with different suppression factors.
  • different types of data can be processed differently, for example, the data points with strong residual tissues described in the foregoing (such as the data points in the box 410 in FIG. 4 or the foregoing are classified as the first
  • a weaker suppression factor can be used for processing (moderate suppression ), in order to appropriately suppress the weak tissue residue in the contrast data without affecting the intensity of the contrast agent, where the suppression factor can be expressed as F, exemplarily, For positive decimal R, it tends to 1 after opening the root sign, in order to achieve a more conservative inhibitory effect.
  • the suppression process may not be required, or the suppression factor F It is equal to 1, that is, no processing is performed. Based on this, the effect of effectively suppressing tissue residue without affecting the intensity of the contrast agent can be achieved. It should be understood that the different types of data points represented by the above boxes 410-430 are only an example, and are not used to limit the present invention.
  • F can also be expressed as a matrix.
  • the point-to-point multiplication of the two matrices can realize the tissue residue printing processing of the angiographic data. As mentioned earlier, F tends to 0, which means that the inhibitory effect is gradually strengthened; F tends to 1, which means that the inhibitory effect is weakened, or even no treatment.
  • the foregoing numerical interval divided by a single threshold, the calculation method of the suppression factor, and the action mode of the suppression factor are only exemplary, and multiple thresholds may be used to divide the numerical interval according to the teachings of the present application, and the suppression factor may also be other.
  • the calculation method (such as other calculation methods other than square and root), and the action method of the suppression factor can also be other methods (such as not point-to-point multiplication, but addition, division, etc.).
  • the contrast-enhanced ultrasound imaging method 300 excavates the difference between the residual tissue and the contrast agent signal from the perspective of time domain difference (inter-frame difference), and can compare the residual part of the tissue and the contrast agent in the contrast data.
  • the signal part is effectively distinguished, so that the residual part of the tissue can be effectively suppressed and the effect of the intensity of the contrast agent is not affected as much as possible.
  • 6A to 7B can be used to check the tissue residue suppression effect obtained by the ultrasound contrast imaging method according to the embodiment of the present application.
  • two scenarios are considered: the residual suppression effect when the contrast agent is not injected; the residual suppression effect after the injection of the contrast agent and the degree of influence on the intensity of the contrast agent.
  • 6A shows an example of a human thyroid contrast dual real-time image processed by the ultrasound contrast imaging method according to the embodiment of the present application when no contrast agent is injected
  • FIG. An example of double real-time images of human thyroid angiography processed by the contrast-enhanced ultrasound imaging method. As shown in FIGS.
  • FIG. 7A shows an example of a canine liver contrast double real-time image after injection of a contrast agent and not processed by the ultrasound contrast imaging method according to an embodiment of the present application.
  • FIG. 7B shows an example of a canine liver contrast double real-time image after injection of a contrast agent and processed by the ultrasound contrast imaging method according to an embodiment of the present application.
  • the near-field residual is significantly suppressed, the contrast agent signal is slightly lost, and the background noise is unchanged, and the desired effect is also achieved.
  • the above exemplarily shows the ultrasound contrast imaging method 300 according to an embodiment of the present invention.
  • a contrast-enhanced ultrasound imaging method according to other embodiments of the present invention will be described with reference to FIGS. 8 to 12.
  • FIG. 8 shows a schematic flowchart of a contrast-enhanced ultrasound imaging method 800 according to another embodiment of the present application. As shown in FIG. 8, the contrast-enhanced ultrasound imaging method 800 may include the following steps:
  • step S810 the ultrasonic probe is controlled to transmit ultrasonic waves to the target object, receive ultrasonic echoes, and obtain ultrasonic echo signals based on the ultrasonic echoes.
  • step S820 at least two frames of contrast signal data are acquired based on the ultrasound echo signal, wherein each frame of contrast signal data is used to generate one frame of contrast image.
  • step S830 the information difference between two adjacent frames of contrast signal data is calculated, and the relationship between the information difference and one frame of contrast signal data in the two adjacent frames of contrast signal data is calculated.
  • step S840 a frame of contrast signal data is divided into regions based on the relationship, and the contrast signal data of different regions is processed with different suppression factors.
  • step S850 a contrast image is generated based on the processed one frame of contrast signal data.
  • the contrast-enhanced ultrasound imaging method 800 at least two frames of contrast signal data are acquired in real-time, and based on the difference in information between two adjacent frames of contrast signal data and the difference in the contrast signal data of two adjacent frames
  • one frame of contrast signal data is divided into regions, the contrast signal data of different regions are processed with different suppression factors, and finally the contrast image is generated based on the processed contrast signal data.
  • the contrast-enhanced ultrasound imaging method 800 in the contrast-enhanced ultrasound imaging method 800 according to the embodiment of the present application, it is also based on the difference in information between two adjacent frames of contrast data.
  • the relationship between the differenced one-frame contrast data in the two-frame contrast data performs tissue residue suppression processing on the differenced one-frame contrast data.
  • the difference is that in the contrast-enhanced ultrasound imaging method 800 according to the embodiment of the present application, it is based on the difference between the information difference between two adjacent frames of contrast signal data and the differenced one frame of contrast signal data.
  • the frame contrast signal data is divided into regions, and contrast images are generated based on the contrast signal data obtained after different regions are processed by different suppression factors.
  • contrast-enhanced ultrasound imaging method 800 of the embodiment of the present application differs from the contrast-enhanced ultrasound imaging method 300 of the embodiment of the present application, and the remaining similarities are only briefly summarized.
  • one frame of contrast signal data can be organized
  • the area with strong residual such as the contrast signal data corresponding to the area in the box 410 shown in Figure 4 above
  • the area with weak tissue residual such as the contrast signal data corresponding to the area in the box 420 shown in Figure 4 above
  • tissue residue such as the contrast signal data corresponding to the area in the box 430 shown in Figure 4 above
  • a frame of contrast signal data can be partitioned based on the relationship (for example, divided into tissue residue Strong areas, areas with weak tissue residues, areas without tissue residues, etc.), and based on the results of region division, different areas of a frame of contrast signal data are suppressed by different suppression factors, so that the generated contrast image can be organized
  • the residual part is effectively suppressed and does not affect the effect of the contrast agent intensity.
  • the contrast signal data of two adjacent frames includes the contrast signal data of the previous frame and the contrast signal data of the next frame.
  • the calculation of the information difference between the contrast signal data of the two adjacent frames in step S830 may include : Calculate the forward difference result between the contrast signal data of the next frame and the contrast signal data of the previous frame; the relationship between the calculated information difference in step S830 and one frame of contrast signal data in two adjacent frames of contrast signal data includes: Calculate the relationship between the forward difference result and the contrast signal data of the previous frame.
  • the contrast signal data of two adjacent frames includes the contrast signal data of the previous frame and the contrast signal data of the next frame.
  • the calculation of the information difference between the contrast signal data of the two adjacent frames in step S830 includes: Calculate the backward difference result between the contrast signal data of the previous frame and the contrast signal data of the next frame; the relationship between the calculated information difference in step S830 and one frame of contrast signal data in two adjacent frames of contrast signal data includes: calculation The backward difference result is related to the size of the contrast signal data of the next frame.
  • the contrast signal data acquired based on the ultrasound echo signal in step S820 may be the contrast signal data obtained after at least beam synthesis processing is performed on the ultrasound echo signal.
  • the information difference includes amplitude difference, phase difference or frequency difference.
  • the contrast signal data acquired based on the ultrasound echo signal in step S820 is the contrast signal data obtained after at least performing analog-to-digital conversion processing on the ultrasound echo signal.
  • the difference in information includes the difference in amplitude.
  • step S830 the information difference between two adjacent frames of contrast signal data is calculated, and the relationship between the information difference and one frame of contrast signal data in the two adjacent frames of contrast signal data is calculated, It may include: calculating the difference data between two adjacent frames of contrast signal data, performing envelope detection on the difference data and the differenced frame of contrast signal data, respectively, to obtain the respective amplitudes of the difference data and the differenced frame of contrast signal data Data; based on the amplitude data of the differential data and the amplitude data of the differenced one frame of contrast signal data, calculate the amplitude ratio of the difference data and the differenced one frame of contrast signal data to obtain the difference between the difference data and the differenced one frame of contrast signal data The result of the relationship.
  • logarithmic compression may also be performed on the result obtained after envelope detection, and the area division of a frame of contrast signal data to be differentiated based on the amplitude ratio may be based on the logarithmically compressed amplitude ratio.
  • step S830 the information difference between two adjacent frames of contrast signal data is calculated, and the information difference is calculated between the difference of one frame of contrast signal data in the adjacent two frames of contrast signal data
  • the relationship may include: calculating the difference data between two adjacent frames of contrast signal data, calculating the ratio of the difference data and the differenced one frame of contrast signal data; performing envelope detection on the comparison value, and obtaining the difference data and the differenced one
  • the amplitude ratio of the frame contrast signal data is taken as the result of the relationship between the difference data and the differenced one frame of contrast signal data.
  • step S830 the information difference between the contrast signal data of two adjacent frames is calculated, and the difference between the information difference and the contrast signal data of one frame of the contrast signal data of the adjacent two frames is calculated.
  • the relationship may include: calculating the difference data between two adjacent frames of contrast signal data, and calculating the respective phase data of the difference data and the differenced one frame of contrast signal data; the phase data based on the difference data and the differenced one frame of contrast data.
  • the phase data of the signal data calculates the phase difference value between the difference data and the differenced frame of contrast signal data, and obtains the relationship result between the difference data and the differenced frame of contrast signal data.
  • step S830 the information difference between two adjacent frames of contrast signal data is calculated, and the information difference is calculated between the difference of one frame of contrast signal data in the adjacent two frames of contrast signal data
  • the relationship may include: calculating the difference data between two adjacent frames of contrast signal data, calculating the ratio between the difference data and the differenced one frame of contrast signal data; performing the argument detection on the comparison value, and obtaining the difference data and the differenced one
  • the phase difference value of the frame contrast signal data is the result of the relationship between the difference data and the differenced one frame of contrast signal data.
  • the difference data of two adjacent frames of contrast signal data may include a first matrix
  • the differenced frame of contrast signal data may include a second matrix.
  • one frame of contrast signal data is determined based on the relationship.
  • the area division may include: determining the numerical interval to which the relationship result between the data point in the first matrix and the data point at the corresponding position in the second matrix belongs; and dividing the corresponding position in the second matrix based on the numerical interval to which the relationship result belongs.
  • the data points at are divided into corresponding areas.
  • data points in different regions belonging to different numerical intervals in the second matrix can be divided by pre-setting several numerical intervals.
  • the numerical interval is determined based on one or more preset thresholds. Obviously, the greater the number of preset thresholds, the more numerical intervals are divided, and the finer the area division of one frame of contrast signal data.
  • a numerical interval based on a preset threshold T (where the value range of T is, for example, 0 to 1) is used as an example for description.
  • the preset numerical interval may include a first numerical interval and a second numerical interval. The first numerical interval is, for example, (0, T], and the second numerical interval is, for example, (T, 1).
  • the relationship result of the data points belongs to the first numerical interval, that is, when R ⁇ (0,T], the data point at the corresponding position in the second matrix is divided into the data point of the first region, For example, the data points in the near-field residual area shown in the figure.
  • the relationship result between any data point in the first matrix and the data point at the corresponding position in the second matrix belongs to the second value
  • the interval that is, when R ⁇ (T,1]
  • the data points at the corresponding positions in the second matrix are divided into data points in the second area, such as the data points in the midfield residual area shown in the figure.
  • the first area There are more tissue residues in the second region than in the second region.
  • the numerical interval based on a preset threshold T may also include a third numerical interval, that is, a numerical interval (1,+ ⁇ ) ⁇ (- ⁇ ,0).
  • a third numerical interval that is, a numerical interval (1,+ ⁇ ) ⁇ (- ⁇ ,0).
  • the suppression processing can be performed in a targeted manner.
  • the residual tissue in the near-field residual area is the most obvious, the residual tissue in the mid-field residual area is weak, and the noise area is basically free of tissue residual. Therefore, in the embodiment of the present application, the data points in the near-field residual area are processed by the first suppression factor, and the data points in the midfield residual area are processed by the second suppression factor.
  • the suppression effect of the first suppression factor is greater than Inhibition of the second inhibitor.
  • the second inhibitory factor can be expressed as F2, exemplarily, For positive decimal R, it tends to 1 after opening the root sign, in order to achieve a more conservative inhibitory effect.
  • the data points in the noise area are processed by a third suppression factor, and the suppression effect of the third suppression factor is smaller than that of the second suppression factor.
  • the differenced one frame of contrast signal data can be expressed as a second matrix
  • the data points in different regions in the second matrix can be processed by different suppression factors. Therefore, the suppression factor can also be expressed as a matrix.
  • the second matrix performs point-to-point multiplication to realize the tissue residue printing processing of the angiographic data.
  • the inhibitory factor tends to 0, which means that the inhibitory effect is gradually strengthened; the inhibitory factor tends to 1, it means that the inhibitory effect is weakened, or even no treatment.
  • the foregoing numerical interval divided by a single threshold, the calculation method of the suppression factor, and the action mode of the suppression factor are only exemplary, and multiple thresholds may be used to divide the numerical interval according to the teachings of the present application, and the suppression factor may also be other.
  • the calculation method (such as other calculation methods other than square and root), and the action method of the suppression factor can also be other methods (such as not point-to-point multiplication, but addition, division, etc.).
  • the contrast-enhanced ultrasound imaging method 700 excavates the difference between the residual tissue and the contrast agent signal from the perspective of time domain difference (inter-frame difference), and can compare the residual tissue area in the contrast signal data with the contrast agent signal.
  • the agent signal area is effectively distinguished, so that the residual area of the tissue can be effectively suppressed and the effect of the contrast agent intensity is not affected as much as possible.
  • FIG. 9 shows a schematic flowchart of a contrast-enhanced ultrasound imaging method 900 according to another embodiment of the present application.
  • the contrast-enhanced ultrasound imaging method 900 may include the following steps:
  • step S910 the ultrasonic probe is controlled to transmit ultrasonic waves to the target object, receive ultrasonic echoes, and obtain ultrasonic echo signals based on the ultrasonic echoes.
  • step S920 the contrast signal data is acquired based on the ultrasound echo signal, and a contrast image is generated based on the contrast signal data.
  • step S930 at least two frames of contrast images are acquired, the information difference between two adjacent frames of contrast images is calculated, and the relationship between the information difference and one of the two adjacent frames of contrast images is calculated.
  • step S940 a frame of contrast image is divided into regions based on the relationship, and different regions of a frame of contrast image are processed with different suppression factors.
  • the ultrasound contrast imaging method 900 at least two frames of contrast image data are acquired in real-time, and based on the difference in information between two adjacent frames of contrast image data and the difference between the two adjacent frames of contrast image data For the relationship between one frame of contrast image data, one frame of contrast image data is divided into regions, and the contrast image data of different regions is processed with different suppression factors.
  • the contrast-enhanced ultrasound imaging method 900 is also based on the difference in information between two adjacent frames of contrast data.
  • the relationship between one frame of contrast data in the two frames of contrast data performs tissue residue suppression processing on one frame of contrast data.
  • the difference is that in the ultrasound contrast imaging method 800 according to the embodiment of the present application, it is based on the difference between the information difference between two adjacent frames of contrast image data and the difference between one frame of contrast image data in two adjacent frames of contrast image data.
  • the relationship divides a frame of contrast image data into regions, instead of dividing a frame of contrast signal based on the information difference between two adjacent frames of contrast signal data and the relationship between one frame of contrast signal data in two adjacent frames of contrast signal data
  • the data is divided into regions.
  • one frame of contrast image data can be organized Areas with strong residues (such as the area within the box 410 shown in Figure 4 above) and areas with weak tissue residues (such as the area within the box 420 shown in Figure 4 above) and regions with basically no tissue residues (such as the above
  • the area within the frame 430 shown in FIG. 4 is clearly distinguished, so that a frame of contrast image data can be partitioned based on the relationship (for example, divided into areas with strong tissue residue, areas with weak tissue residue, and no tissue residue).
  • different areas of a frame of contrast image data are suppressed by different suppression factors, so that the residual part of the processed contrast image can be effectively suppressed without affecting the intensity of the contrast agent. Effect.
  • two adjacent frames of contrast images include the previous frame of contrast images and the next frame of contrast images
  • the calculation of the information difference between the two adjacent frames of contrast images in step S930 may include: calculating the latter The result of the forward difference between a frame of contrast image and the previous frame of contrast image
  • the relationship between the calculated information difference in step S930 and one frame of contrast image in two adjacent frames of contrast image may include: the calculation of the forward difference result and the previous frame The gray scale ratio of a contrast image.
  • two adjacent frames of contrast image include the previous frame of contrast image and the next frame of contrast image
  • the calculation of the information difference between the two adjacent frames of contrast image in step S930 may include: calculating the previous frame The result of the backward difference between one frame of contrast image and the next frame of contrast image
  • the relationship between the calculated information difference in step S930 and one of the two adjacent frames of contrast image may include: the result of calculating the backward difference and the latter one The gray ratio of the contrast image of the frame.
  • the difference in information includes a first matrix
  • a frame of contrasted image to be differentiated includes a second matrix.
  • dividing the frame of contrasted image into regions may include: determining the difference in the first matrix The numerical interval to which the gray scale ratio of the data point and the data point at the corresponding position in the second matrix belongs; and the data point at the corresponding position in the second matrix is divided into the corresponding area based on the numerical interval to which the gray scale ratio belongs.
  • the numerical interval is determined based on one or more preset thresholds.
  • the numerical interval includes a first numerical interval and a second numerical interval determined based on a preset threshold; the first numerical interval corresponds to the first region, the second numerical interval corresponds to the second region, and More tissue residues in the second area than in the second area; when the gray scale ratio of any data point in the first matrix to the data point at the corresponding position in the second matrix belongs to the first numerical interval, the first area The data points in the second matrix are processed by the first suppression factor; when the grayscale ratio of any data point in the first matrix to the data point at the corresponding position in the second matrix belongs to the second numerical interval, the second area is The data points of using the second inhibitory factor for processing, the inhibitory effect of the first inhibitory factor is greater than the inhibitory effect of the second inhibitory factor.
  • the first area is the near field residual area
  • the second area is the midfield residual area
  • the numerical interval further includes a third numerical interval corresponding to the noise area; when the grayscale ratio of any data point in the first matrix to the data point at the corresponding position in the second matrix belongs to the third In the numerical interval, the data points in the noise area are processed by the third suppression factor, and the suppression effect of the third suppression factor is less than the suppression effect of the second suppression factor.
  • the third suppression factor is equal to 1.
  • the contrast-enhanced ultrasound imaging method 900 excavates the difference between the residual tissue and the contrast agent signal from the perspective of time domain difference (inter-frame difference), and can compare the residual tissue area in the contrast image data with the contrast agent signal.
  • the agent image area is effectively distinguished, so that the residual area of the tissue can be effectively suppressed and the effect of the contrast agent strength is not affected as much as possible.
  • FIG. 10 is a schematic flowchart of a contrast-enhanced ultrasound imaging method 1000 according to another embodiment of the present application. As shown in FIG. 10, the contrast-enhanced ultrasound imaging method 100 may include the following steps:
  • step S1010 the ultrasonic probe is controlled to transmit ultrasonic waves to the target object, receive echoes of the ultrasonic waves, and obtain ultrasonic echo signals based on the echoes of the ultrasonic waves.
  • step S1020 at least two frames of contrast signal data are acquired based on the ultrasound echo signal, wherein each frame of contrast signal data is used to generate one frame of contrast image.
  • step S1030 the information difference between two adjacent frames of contrast signal data is calculated, and the relationship between the information difference and one frame of contrast signal data in the two adjacent frames of contrast signal data is calculated.
  • step S1040 the one frame of contrast signal data is classified based on the relationship, and at least part of the data of the one frame of contrast signal data is suppressed based on the classification result.
  • step S1050 a contrast image is generated based on the processed one frame of contrast signal data.
  • the contrast-enhanced ultrasound imaging method 1000 according to the embodiment of the present application is similar to the foregoing contrast-enhanced ultrasound imaging method 300 according to the embodiment of the present application as a whole, except that the contrast-enhanced ultrasound imaging method 1000 according to the embodiment of the present application is real-time in steps S1010-S1020.
  • Two frames of contrast data are acquired, and the acquired contrast data is contrast signal data used to generate a contrast image.
  • the contrast signal data acquired in steps S1010-S1020 may be baseband data, and the information difference between the contrast signal data of two adjacent frames includes amplitude difference, phase difference or frequency difference.
  • the contrast-enhanced ultrasound imaging method 1000 excavates the difference between the residual tissue and the contrast agent signal from the perspective of time domain difference (inter-frame difference), and can compare the residual tissue area in the contrast signal data with the contrast agent signal.
  • the agent signal area is effectively distinguished, so that the residual area of the tissue can be effectively suppressed and the effect of the contrast agent intensity is not affected as much as possible.
  • FIG. 11 shows a schematic flowchart of a contrast-enhanced ultrasound imaging method 1100 according to another embodiment of the present application. As shown in FIG. 10, the contrast-enhanced ultrasound imaging method 1100 may include the following steps:
  • step S1110 the ultrasonic probe is controlled to transmit ultrasonic waves to the target object, receive ultrasonic echoes, and obtain ultrasonic echo signals based on the ultrasonic echoes.
  • step S1120 the contrast signal data is acquired based on the ultrasound echo signal, and a contrast image is generated based on the contrast signal data.
  • step S1130 the information difference between two adjacent contrast images is calculated, and the relationship between the information difference and one of the two adjacent contrast images is calculated.
  • step S1140 the one-frame contrast image is classified based on the relationship, and at least part of the one-frame contrast image is suppressed based on the classification result.
  • the contrast-enhanced ultrasound imaging method 1100 according to the embodiment of the present application is similar to the foregoing contrast-enhanced ultrasound imaging method 300 according to the embodiment of the present application as a whole, except that the contrast-enhanced ultrasound imaging method 1100 according to the embodiment of the present application acquires in real time in steps S1110-S1120 At least two contrast images.
  • calculating the information difference between two adjacent frames of contrast images in step S1130 may include: calculating the gray difference result between two adjacent frames of contrast images; calculating the information difference in step S1130
  • the relationship with one frame of contrast images in two adjacent frames of contrast images may include: calculating the gray level ratio between the gray difference result and one frame of contrast images in two adjacent frames of contrast images.
  • the contrast-enhanced ultrasound imaging method 1100 excavates the difference between the residual tissue and the contrast agent signal from the perspective of time domain difference (inter-frame difference), and can compare the residual tissue area in the contrast image data with the contrast agent signal.
  • the agent image area is effectively distinguished, so that the residual area of the tissue can be effectively suppressed and the effect of the contrast agent strength is not affected as much as possible.
  • an embodiment of the present invention also provides a contrast-enhanced ultrasound imaging apparatus 100, which can be used to implement the above-mentioned contrast-enhanced imaging method.
  • the ultrasound contrast imaging apparatus 100 may include a probe 110, a transmitting circuit 112, a receiving circuit 114, a beam combining circuit 116, a processor 118 and a display 120, a transmitting/receiving selection switch 122, and some or all of the components in the memory 124, and the correlation of each component
  • the description can refer to the above.
  • the transmitting circuit 112 is used to excite the ultrasonic probe 110 to transmit ultrasonic waves to the target object.
  • the receiving circuit 114 is used to control the ultrasonic probe 110 to receive the ultrasonic echo returned from the target object to obtain the ultrasonic echo signal.
  • the processor 118 may be used to execute the above-mentioned contrast-enhanced ultrasound imaging methods 300, 800, 900, 1000, and 1100. Hereinafter, only the main functions of the ultrasound contrast imaging apparatus 100 are described, and the details that have been described above are omitted.
  • the processor 118 calculates the information difference between two adjacent frames of contrast data, and calculates the relationship between the information difference and one frame of contrast data in the two adjacent frames of contrast data; The relationship classifies the one frame of contrast data, and performs suppression processing on at least part of the one frame of contrast data based on the classification result.
  • the difference in information between two adjacent frames of contrast data includes an amplitude difference, a phase difference, or a frequency difference.
  • the processor 118 calculates the difference data between the two adjacent frames of contrast data, and calculates the difference data and the one.
  • the amplitude data of each frame of contrast data based on the amplitude data of the differential data and the amplitude data of the one frame of contrast data, the amplitude ratio of the differential data and the one frame of contrast data is calculated to obtain the difference data and the The result of the relationship between one frame of contrast data.
  • the processor 118 calculates the difference data between the two adjacent frames of contrast data, and compares the difference data with the one frame. Calculating the ratio of the contrast data; and solving the amplitude of the ratio to obtain the ratio of the amplitude of the difference data to the one frame of contrast data, which is used as the result of the relationship between the difference data and the one frame of contrast data.
  • the processor 118 calculates the difference data between the two adjacent frames of contrast data, and calculates the difference data and the one The phase data of each frame of contrast data; the phase difference between the difference data and the one frame of contrast data is calculated based on the phase data of the differential data and the phase data of the one frame of contrast data, to obtain the difference data and the phase data Describe the results of the relationship between one frame of contrast data.
  • the processor 118 calculates the difference data between the two adjacent frames of contrast data, and compares the difference data with the one frame. Calculating the ratio of the contrast data; and solving the phase of the ratio to obtain the phase difference value between the difference data and the one frame of contrast data as a result of the relationship between the difference data and the one frame of contrast data.
  • the difference data includes a first matrix
  • the one frame of contrast data includes a second matrix
  • the processor 118 is configured to determine the data points in the first matrix and the corresponding data points in the second matrix.
  • the numerical interval to which the relationship result of the data point at the position belongs, and the data point at the corresponding position in the second matrix is classified as belonging to the corresponding category based on the numerical interval to which the relationship result belongs.
  • the processor 118 is configured to determine, based on the result of the classification, data that needs to be suppressed in the one frame of angiographic data; and use different suppression for different types of data in the data that needs to be suppressed. Factor processing.
  • the numerical interval is determined based on one or more preset thresholds.
  • the numerical interval includes a first numerical interval and a second numerical interval determined based on the preset threshold; the first numerical interval corresponds to a first type of data point, and the second numerical interval corresponds to a first numerical interval.
  • Two types of data points, and the tissue residues at the first type data points are stronger than the tissue residues at the second type data points; when any data point in the first matrix and the second matrix When the result of the relationship between the data points at the corresponding positions belongs to the first numerical interval, the data points at the corresponding positions in the second matrix are divided into data points of the first category; when any of the data points in the first matrix When the result of the relationship between a data point and the data point at the corresponding position in the second matrix belongs to the second numerical interval, the data point at the corresponding position in the second matrix is divided into a second type of data point ; Using a first inhibitor to process the first category of data points, using a second inhibitor to process the second category of data points, the inhibitory effect of the first inhibitory factor is greater than that of the second
  • the processor 118 calculates the information difference between two adjacent frames of contrast data, and calculates the relationship between the information difference and one frame of contrast data in the two adjacent frames of contrast data; The relationship divides the one frame of contrast data into regions, and uses different suppression factors to process the contrast data in different regions.
  • the contrast data may be contrast signal data used to generate one frame of contrast image, or one frame of contrast image data.
  • FIG. 12 shows a schematic block diagram of an ultrasound contrast imaging apparatus 1200 according to an embodiment of the present application.
  • the ultrasound contrast imaging apparatus 1200 includes a memory 1210 and a processor 1220.
  • the memory 1210 stores programs for implementing corresponding steps in the ultrasound contrast imaging methods 300, 800, 900, 1000, and 1100 according to the embodiments of the present application.
  • the processor 1220 is configured to run a program stored in the memory 1210 to execute corresponding steps of the ultrasound contrast imaging methods 300, 800, 900, 1000, and 1100 according to the embodiments of the present application.
  • a storage medium on which program instructions are stored, and when the program instructions are run by a computer or a processor, they are used to execute the ultrasound contrast imaging The corresponding steps of the method.
  • 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), 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.
  • a computer program is also provided, and the computer program can be stored in a cloud or a local storage medium.
  • the computer program is run by a computer or a processor, it is used to execute the corresponding steps of the ultrasound contrast imaging method of the embodiment of the present application.
  • the contrast-enhanced ultrasound imaging method, device, and storage medium excavate the difference between the tissue residue and the contrast agent signal from the perspective of time domain difference (inter-frame difference), and can detect the tissue residue in the contrast data.
  • the part is effectively distinguished from the signal part of the contrast agent, so that the residual part of the tissue can be effectively suppressed and the effect of the intensity of the contrast agent is not affected as much as possible.
  • the disclosed device and method may be implemented in other ways.
  • the device embodiments described above are merely illustrative, for example, the division of the units is only a logical function division, and there may be other divisions in actual implementation, for example, multiple units or components may be combined or It can be integrated into another device, or some features can be ignored or not implemented.
  • the various component embodiments of the present application may be implemented by hardware, or by software modules running on one or more processors, or by a combination of them.
  • a microprocessor or a digital signal processor (DSP) may be used in practice to implement some or all of the functions of some modules according to the embodiments of the present application.
  • This application can also be implemented as a device program (for example, a computer program and a computer program product) for executing part or all of the methods described herein.
  • Such a program for implementing the present application may be stored on a computer readable medium, or may 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.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Veterinary Medicine (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Hematology (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

Procédé et appareil d'imagerie de contraste ultrasonore, et support de stockage. Le procédé consiste : à acquérir au moins deux trames de données de contraste (S310); à calculer une différence d'informations entre deux trames adjacentes de données de contraste, et à calculer une relation entre la différence d'informations et une trame de données de contraste parmi les deux trames adjacentes de données de contraste (S320); et à classifier la trame de données de contraste en fonction de la relation, et à effectuer un traitement de suppression sur au moins certaines données de la trame de données de contraste en fonction d'un résultat de classification (S330). Le procédé permet d'explorer la différence entre des résidus de tissu et un signal d'agent de contraste du point de vue d'une différence de domaine temporel (une différence entre trames), de sorte qu'une partie résidu de tissu des données de contraste et une partie de son signal d'agent de contraste puissent être distinguées efficacement l'une de l'autre, la partie résidu de tissu pouvant être ainsi supprimée efficacement, et son influence sur l'effet de l'intensité d'un agent de contraste étant réduite dans la plus grande mesure possible.
PCT/CN2020/090855 2020-05-18 2020-05-18 Procédé et appareil d'imagerie de contraste ultrasonore, et support de stockage WO2021232192A1 (fr)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202080006417.2A CN113164160B (zh) 2020-05-18 2020-05-18 超声造影成像方法、装置和存储介质
PCT/CN2020/090855 WO2021232192A1 (fr) 2020-05-18 2020-05-18 Procédé et appareil d'imagerie de contraste ultrasonore, et support de stockage

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2020/090855 WO2021232192A1 (fr) 2020-05-18 2020-05-18 Procédé et appareil d'imagerie de contraste ultrasonore, et support de stockage

Publications (1)

Publication Number Publication Date
WO2021232192A1 true WO2021232192A1 (fr) 2021-11-25

Family

ID=76879265

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2020/090855 WO2021232192A1 (fr) 2020-05-18 2020-05-18 Procédé et appareil d'imagerie de contraste ultrasonore, et support de stockage

Country Status (2)

Country Link
CN (1) CN113164160B (fr)
WO (1) WO2021232192A1 (fr)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106028947A (zh) * 2015-05-07 2016-10-12 深圳迈瑞生物医疗电子股份有限公司 血流成像方法及系统
CN107361793A (zh) * 2017-07-18 2017-11-21 深圳开立生物医疗科技股份有限公司 超声波成像方法、系统及超声成像设备
CN108135566A (zh) * 2016-06-30 2018-06-08 北京深迈瑞医疗电子技术研究院有限公司 一种增强超声造影图像的方法、系统及超声造影成像设备
CN111012381A (zh) * 2019-12-30 2020-04-17 无锡祥生医疗科技股份有限公司 超声造影方法、装置及存储介质
US10646203B2 (en) * 2015-10-08 2020-05-12 Samsung Medison Co., Ltd. Ultrasound diagnosis method and apparatus for analyzing contrast enhanced ultrasound image

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106061396B (zh) * 2015-04-08 2019-05-21 深圳迈瑞生物医疗电子股份有限公司 一种超声造影成像方法及装置
US20190154821A1 (en) * 2017-11-21 2019-05-23 The Charles Stark Draper Laboratory, Inc. System and method for imaging and localization of contrast-enhanced features in the presence of accumulating contrast agent in a body
CN110458836A (zh) * 2019-08-16 2019-11-15 深圳开立生物医疗科技股份有限公司 一种超声造影成像方法、装置和设备及可读存储介质

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106028947A (zh) * 2015-05-07 2016-10-12 深圳迈瑞生物医疗电子股份有限公司 血流成像方法及系统
US10646203B2 (en) * 2015-10-08 2020-05-12 Samsung Medison Co., Ltd. Ultrasound diagnosis method and apparatus for analyzing contrast enhanced ultrasound image
CN108135566A (zh) * 2016-06-30 2018-06-08 北京深迈瑞医疗电子技术研究院有限公司 一种增强超声造影图像的方法、系统及超声造影成像设备
CN107361793A (zh) * 2017-07-18 2017-11-21 深圳开立生物医疗科技股份有限公司 超声波成像方法、系统及超声成像设备
CN111012381A (zh) * 2019-12-30 2020-04-17 无锡祥生医疗科技股份有限公司 超声造影方法、装置及存储介质

Also Published As

Publication number Publication date
CN113164160A (zh) 2021-07-23
CN113164160B (zh) 2022-09-16

Similar Documents

Publication Publication Date Title
US11635514B2 (en) Imaging methods and apparatuses for performing shear wave elastography imaging
US6368277B1 (en) Dynamic measurement of parameters within a sequence of images
US9002080B2 (en) Singular value filter for imaging or detection
US11000263B2 (en) Ultrasound diagnostic apparatus, image processing device, and image processing method
KR102286299B1 (ko) 초음파 이미지 디스플레이 장치 및 초음파 이미지를 디스플레이하기 위한 방법
CN103202714B (zh) 超声波诊断装置、医用图像处理装置以及医用图像处理方法
US20120154400A1 (en) Method of reducing noise in a volume-rendered image
US20150359507A1 (en) Ultrasound diagnosis apparatus and ultrasound image processing method
JP2011125690A (ja) 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム
US20240050062A1 (en) Analyzing apparatus and analyzing method
Hu et al. Acoustic shadow detection: study and statistics of B-mode and radiofrequency data
US20060030777A1 (en) T-statistic method for suppressing artifacts in blood vessel ultrasonic imaging
US10537305B2 (en) Detecting amniotic fluid position based on shear wave propagation
Magud et al. Medical ultrasound image speckle noise reduction by adaptive median filter
US8801616B2 (en) Ultrasonic diagnostic apparatus and ultrasonic image processing method
WO2021232192A1 (fr) Procédé et appareil d'imagerie de contraste ultrasonore, et support de stockage
CN112890866A (zh) 超声成像方法、系统及计算机可读存储介质
JP2019181183A (ja) 医用診断装置、医用画像処理装置、及び画像処理プログラム
JP2020110296A (ja) 超音波撮像装置及び画像処理方法
US20220087650A1 (en) Ultrasound imaging system and method for generating an enhanced image
WO2020047805A1 (fr) Procédé et dispositif d'imagerie élastique d'onde de cisaillement et support de stockage informatique
JP7109986B2 (ja) 画像処理装置及び超音波診断装置
JP2020054815A (ja) 解析装置及び解析プログラム
KR100656861B1 (ko) 볼륨 데이터를 렌더링하는 영상 처리 시스템 및 방법
US20230037641A1 (en) Elastography method, system and storage medium

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 20936826

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 20936826

Country of ref document: EP

Kind code of ref document: A1