Visual tracking system and method for dynamic distribution of ultrasonic contrast agent

Abstract

The invention relates to the technical field of ultrasonic contrast, in particular to a visual tracking system and method for dynamic distribution of an ultrasonic contrast agent, and the method comprises the steps: obtaining sound pressure and microbubble breakage data of ultrasonic contrast, analyzing the influence of the sound pressure on microbubble breakage, and obtaining a breakage influence coefficient; sound pressure and microbubble movement data of ultrasonic contrast are obtained, the influence of the sound pressure on microbubble movement is analyzed, and a movement influence coefficient is obtained; according to the breaking influence coefficient and the movement influence coefficient, the sound pressure is dynamically regulated and controlled, and regulated and controlled sound pressure is obtained; analyzing the vulnerability coefficient of the current microbubble according to the damage influence of sound waves on the microbubble within the scanning time; the regulation and control sound pressure is corrected through the vulnerability coefficient of the current microbubble, and optimized regulation and control sound pressure is obtained; adjusting the tracking area according to the current microbubble fracture rate and the vulnerability coefficient of the current microbubble; according to the method, the service life of the microbubbles is prolonged, stable development is maintained, non-specific rupture is reduced, biological risks are reduced, and the stability of a tracking area is effectively guaranteed.

Description

Technical Field

[0001] The present invention relates to the technical field of ultrasound contrast imaging, and in particular to a visualization tracking system and method for the dynamic distribution of an ultrasound contrast agent. Background Art

[0002] In the field of ultrasound contrast imaging, microbubbles have been widely used as ultrasound contrast agents in clinical scenarios such as cardiovascular imaging, tumor diagnosis, and treatment evaluation. The dynamic behavior of microbubbles in the acoustic field (such as rupture, migration, and distribution) directly affects imaging quality and treatment efficacy.

[0003] In the paper "Dynamic Table Movement Contrast Agent Bolus Tracking in 3D-MRA of Lower Extremity Arteries," the application of dynamic table movement contrast agent tracking in lower extremity arterial MRA and its application in the diagnosis of vascular lesions were explored. Methods: 17 patients with suspected pelvic and lower extremity vascular lesions underwent 3D-CEMRA, and the scans were compared with surgical findings. Care-bolus scanning was combined with three-segment dynamic table movement scanning. Results: Satisfactory 3D-CEMRA imaging of the abdominal, pelvic, and lower extremity arteries was achieved in all 17 patients, with clear vessels and distinct lesions.

[0004] Current ultrasound contrast imaging technology suffers from rough sound pressure control, inaccurate prediction of microbubble behavior, and lack of damage accumulation management, resulting in unstable imaging quality, low treatment efficiency, and potential biological risks. Summary of the Invention

[0005] The purpose of the present invention is to provide an integrated and efficient method for oxidative removal of hypophosphorous acid in chemical plating wastewater. The technical problem solved by the present invention is that the current ultrasound contrast imaging technology has unstable imaging quality, low treatment efficiency and potential biological risks due to the extensive sound pressure control, inaccurate prediction of microbubble behavior and lack of damage accumulation management.

[0006] The purpose of the present invention can be achieved through the following technical solutions:

[0007] A method for visually tracking the dynamic distribution of an ultrasound contrast agent comprises the following steps:

[0008] During ultrasound scanning, the acoustic pressure and microbubble fragmentation data of ultrasound contrast imaging are obtained, the effect of acoustic pressure on microbubble fragmentation is analyzed, and the fragmentation influence coefficient is obtained;

[0009] During ultrasound scanning, the acoustic pressure and microbubble movement data of ultrasound contrast imaging are obtained, the influence of acoustic pressure on microbubble movement is analyzed, and the movement influence coefficient is obtained;

[0010] According to the crushing influence coefficient and the movement influence coefficient, the sound pressure is dynamically regulated to obtain the regulated sound pressure;

[0011] Then, according to the damage effect of the acoustic wave on the microbubbles during the scanning time, the vulnerability coefficient of the current microbubbles is analyzed;

[0012] The regulated sound pressure is corrected by the current microbubble vulnerability coefficient to obtain the optimized regulated sound pressure;

[0013] Then, the tracking area is adjusted according to the current microbubble rupture rate and the current microbubble vulnerability coefficient.

[0014] As a further solution of the present invention: the process of obtaining the crushing influence coefficient is:

[0015] Obtain the total number of initially injected microbubbles and the number of ruptured microbubbles, and calculate the percentage of the number of ruptured microbubbles to the total number of initially injected microbubbles to obtain the microbubble rupture rate;

[0016] The microbubble rupture rate and sound pressure present an S-shaped curve. The influence of sound pressure on microbubble rupture is analyzed and the fragmentation influence coefficient is obtained.

[0017] As a further solution of the present invention: the process of obtaining the mobile influence coefficient is:

[0018] Acquire microbubble movement data, analyze the effect of sound pressure on microbubble movement, and obtain the movement influence coefficient.

[0019] As a further solution of the present invention, the specific process of dynamically regulating the sound pressure according to the fragmentation influence coefficient and the movement influence coefficient is as follows:

[0020] Based on the fragmentation influence coefficient and the movement influence coefficient, the microbubble concentration and velocity in the future time window are predicted, and the regulated sound pressure is calculated through the PID control model; wherein, the PID control model is:

[0021]

[0022] Among them, K p is the proportional gain, K i is the integral gain, K d is the differential gain, e(t) is the error signal, and P(t) is the regulated sound pressure.

[0023] As a further solution of the present invention: the expression of e(t) is: e(t)=|C(t)-C target |+|V(t)-V target |.

[0024] As a further solution of the present invention, according to the damage effect of the acoustic wave on the microbubbles during the scanning time, the process of analyzing the vulnerability coefficient of the current microbubbles is as follows:

[0025] The sound pressure is normalized, and the historical data is segmented by time window. The damage caused by a single pulse to the microbubble is calculated to obtain the single pulse damage. The decay rate of the microbubble damage over time is also obtained. All single pulse damage and decay rates within the scanning period are comprehensively analyzed to obtain the cumulative damage. The cumulative damage is mapped to obtain the vulnerability coefficient of the current microbubble.

[0026] As a further solution of the present invention, a comprehensive analysis of all single pulse damage and attenuation rates during the scanning period is performed to obtain the cumulative damage process:

[0027] The damage Dj of a single pulse to the microbubble is multiplied by the decay rate of the microbubble damage over time to obtain the single damage. The cumulative damage is obtained by summing up all the single damages within the scanning time.

[0028] As a further solution of the present invention: the calculation process of optimizing and controlling the sound pressure is:

[0029] The current microbubble fragility coefficient is multiplied by the control sound pressure to obtain the optimized control sound pressure.

[0030] As a further solution of the present invention, the process of adjusting the tracking area according to the current microbubble rupture rate and the current microbubble vulnerability coefficient is as follows:

[0031] Obtain the inter-frame change rate of the second harmonic signal intensity within the ROI; calculate the adjustment value of the tracking area using the ROI width adjustment formula;

[0032]

[0033] Wherein, DR is the adjustment value of the tracking region, ΔI is the inter-frame change rate, and W0 is the current ROI width.

[0034] A visual tracking system for dynamic distribution of ultrasound contrast agents, the system comprising:

[0035] Acoustic pressure first impact module: During ultrasound scanning, the acoustic pressure and microbubble fragmentation data of ultrasound contrast imaging are obtained, the influence of acoustic pressure on microbubble fragmentation is analyzed, and the fragmentation influence coefficient is obtained;

[0036] Acoustic pressure second impact module: During ultrasound scanning, it obtains the acoustic pressure and microbubble movement data of ultrasound contrast imaging, analyzes the impact of acoustic pressure on microbubble movement, and obtains the movement influence coefficient;

[0037] Optimization control module: dynamically adjusts the sound pressure according to the crushing influence coefficient and the movement influence coefficient to obtain the regulated sound pressure;

[0038] Then, according to the damage effect of the acoustic wave on the microbubbles during the scanning time, the vulnerability coefficient of the current microbubbles is analyzed;

[0039] The regulated sound pressure is corrected by the current microbubble vulnerability coefficient to obtain the optimized regulated sound pressure;

[0040] Then, the tracking area is adjusted according to the current microbubble rupture rate and the current microbubble vulnerability coefficient.

[0041] Beneficial effects of the present invention:

[0042] The present invention obtains the acoustic pressure and microbubble fragmentation data of ultrasound contrast imaging during the ultrasonic scanning process, analyzes the influence of the acoustic pressure on the microbubble rupture, and obtains the fragmentation influence coefficient; obtains the acoustic pressure and microbubble movement data of ultrasound contrast imaging during the ultrasonic scanning process, analyzes the influence of the acoustic pressure on the microbubble movement, and obtains the movement influence coefficient; dynamically regulates the acoustic pressure according to the fragmentation influence coefficient and the movement influence coefficient to obtain the regulated acoustic pressure; then analyzes the current microbubble vulnerability coefficient according to the damage effect of the acoustic wave on the microbubble during the scanning time; corrects the regulated acoustic pressure according to the current microbubble vulnerability coefficient to obtain the optimized regulated acoustic pressure; and then adjusts the regulated acoustic pressure according to the regulated acoustic pressure according to the regulated acoustic pressure. The current microbubble rupture rate and the current microbubble vulnerability coefficient are used to adjust the tracking area. During the visual tracking of the dynamic distribution of ultrasound contrast agents, the present invention analyzes the influence of the sound pressure distribution on microbubble fragmentation and the influence on microbubble movement, comprehensively analyzes the influence of microbubble fragmentation and microbubble movement, and combines the damage of sound pressure to microbubbles in the historical process, so as to accurately control the sound pressure during the scanning process. This extends the life of microbubbles, maintains stable imaging, reduces nonspecific rupture, reduces biological risks, and effectively ensures the stability of the tracking area, which will greatly improve the imaging quality in ultrasound contrast imaging. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0044] Figure 1 It is a flow chart of a method for visually tracking the dynamic distribution of an ultrasound contrast agent according to the present invention;

[0045] Figure 2 It is a structural diagram of a visualization tracking system for dynamic distribution of ultrasound contrast agents according to the present invention. DETAILED DESCRIPTION

[0046] In order to enable those skilled in the art to better understand the solutions of the present invention, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the embodiments described are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts should fall within the scope of protection of the present invention.

[0047] Example 1

[0048] like Figure 1 As shown, an embodiment of the present invention provides a method for visually tracking the dynamic distribution of an ultrasound contrast agent, which specifically includes the following steps:

[0049] Step 1: During ultrasound scanning, obtain the acoustic pressure and microbubble breakup data of ultrasound contrast imaging, analyze the effect of acoustic pressure on microbubble breakage, and obtain the breakup influence coefficient;

[0050] In some embodiments, during ultrasound scanning, the acoustic pressure of ultrasound contrast imaging is obtained by a hydrophone, and the total number of initially injected microbubbles (calibrated in vitro) and the number of ruptured microbubbles (by counting signal dips or noise events) are obtained. The percentage of the number of ruptured microbubbles to the total number of initially injected microbubbles is calculated to obtain the microbubble rupture rate.

[0051] The microbubble rupture rate and sound pressure usually present an S-shaped curve. The influence of sound pressure on microbubble rupture is analyzed by formula (1) to obtain the fragmentation influence coefficient;

[0052]

[0053] Where R max is the maximum rupture rate; P O is the sound pressure corresponding to the half-maximum rupture rate (critical threshold), P is the sound pressure of contrast-enhanced ultrasound; K1 is the fragmentation influence coefficient, and R is the goodness of fit;

[0054] For example, the physical meanings of the parameters in the calculation formula of the crushing influence coefficient are:

[0055]

[0056]

[0057] Step 2: During ultrasound scanning, obtain the acoustic pressure and microbubble movement data of ultrasound contrast imaging, analyze the effect of acoustic pressure on microbubble movement, and obtain the movement influence coefficient;

[0058] In some embodiments, during ultrasound scanning, the acoustic pressure of ultrasound contrast imaging is obtained by a hydrophone, and microbubble movement data (optical flow method or particle image velocimetry) is obtained. The influence of the acoustic pressure on the movement of microbubbles is analyzed by formula (2) to obtain the movement influence coefficient;

[0059] ΔV=K2*P norm +A formula (2)

[0060] Where, ΔV is the change in microbubble migration velocity, K2 is the migration influence coefficient, Pnorm is the normalized sound pressure of contrast-enhanced ultrasound, and A is the intercept;

[0061] Step 3: Dynamically control the sound pressure according to the crushing influence coefficient and the movement influence coefficient to obtain the controlled sound pressure;

[0062] Then, according to the damage effect of the acoustic wave on the microbubbles during the scanning time, the vulnerability coefficient of the current microbubbles is analyzed;

[0063] The regulated sound pressure is corrected by the current microbubble vulnerability coefficient to obtain the optimized regulated sound pressure;

[0064] Then, the tracking area is adjusted according to the current microbubble rupture rate and the current microbubble vulnerability coefficient;

[0065] In step 3, first, specifically, the sound pressure is dynamically regulated according to the crushing influence coefficient and the movement influence coefficient. The specific process of obtaining the regulated sound pressure is:

[0066] Obtain the crushing influence coefficient and the movement influence coefficient, and construct the objective function J = a1*|C(t)-C target |+a2*|V(t)-V target |, where J is the objective function, which aims to accurately control the concentration and velocity of microbubbles while ensuring the survival rate of microbubbles (avoiding excessive rupture) to maximize the imaging quality or treatment effect. C(t) is the microbubble concentration, V(t) is the microbubble velocity, and C target is the preset target value of microbubble concentration, V target is the preset target value of microbubble velocity, a1 and a2 are weight coefficients, balancing the concentration and velocity priorities;

[0067] And set the constraints, the sound pressure range P min ≤P≤P max ; Rupture rate safety threshold R≤R threshold ;

[0068] Based on the fragmentation influence coefficient and the movement influence coefficient, the microbubble concentration and velocity in the future time window are predicted, and the regulated sound pressure is calculated through the PID control model; wherein, the PID control model is:

[0069]

[0070] Where Kp is the proportional gain, Ki is the integral gain, Kd is the differential gain, and e(t) = |C(t)-C target |+|V(t)-V target |, is the error signal, P(t) is the regulated sound pressure;

[0071] Secondly, specifically, according to the damage effect of the acoustic wave on the microbubbles during the scanning time, the process of analyzing the current microbubble vulnerability coefficient is as follows:

[0072] The acoustic pressure is normalized, and the historical data is segmented by time window. The damage caused by a single pulse to the microbubble is calculated to obtain the single pulse damage. The decay rate of the microbubble damage over time is also obtained. All single pulse damage and decay rates within the scanning period are comprehensively analyzed to obtain the cumulative damage. The cumulative damage is mapped to obtain the vulnerability coefficient of the current microbubble.

[0073] For example, the calculation process of the damage of a single pulse to microbubbles to obtain the single pulse damage is:

[0074] The damage Dj of a single pulse to microbubbles is related to its sound pressure Pj and duration tj. j =k*P j α , where K and α are both influence coefficients;

[0075] Comprehensively analyzing all single pulse damage and attenuation rates during the scanning period, the cumulative damage process is obtained as follows:

[0076] The damage Dj of a single pulse to the microbubble is calculated by multiplying the decay rate of the microbubble damage over time to obtain the single damage. The cumulative damage is obtained by summing up all the single damages within the scanning time.

[0077] The decay rate of microbubble damage over time is obtained by:

[0078] Obtain the same pulse and the rupture rate corresponding to the pulse, record the first rupture rate of the first pulse as the first rupture rate, and the second rupture rate of the second pulse as the second rupture rate. Calculate the difference between the first rupture rate and the natural rupture rate when there is no pulse to obtain the first rupture rate difference. Calculate the difference between the first rupture rate and the natural rupture rate when there is no pulse to obtain the second rupture rate difference. Calculate the ratio of the second rupture rate difference to the first rupture rate difference to obtain the residual damage effect after damage attenuation.

[0079] According to the damage model D = 1 + e -λΔt , the decay rate λ is calculated by fitting the experimental data through nonlinear regression, where D is the residual damage effect after damage decay, and Δt is the interval between the first pulse and the second pulse;

[0080] The cumulative damage mapping is used to obtain the vulnerability coefficient of the current microbubble. The mapping formula is:

[0081]

[0082] Among them, X(t) is the vulnerability coefficient of the current microbubble, D threshold is the damage threshold of microbubble rupture, D(t) is the cumulative damage;

[0083] Thirdly, specifically, the calculation process for optimizing and controlling the sound pressure is:

[0084] The current microbubble vulnerability coefficient is multiplied by the control sound pressure to obtain the optimized control sound pressure;

[0085] Fourthly, specifically, the process of adjusting the tracking area according to the current microbubble rupture rate and the current microbubble vulnerability coefficient is:

[0086] Obtain the inter-frame rate of change of the second harmonic (2fo) signal intensity within the ROI; calculate the adjustment value of the tracking area using the ROI width adjustment formula;

[0087]

[0088] Where DR is the adjustment value of the tracking area, ΔI is the inter-frame change rate, and W0 is the current ROI width.

[0089] The technical solution of the embodiment of the present invention is as follows: during the ultrasonic scanning process, the acoustic pressure and microbubble fragmentation data of the ultrasonic contrast imaging are obtained, the influence of the acoustic pressure on the microbubble burst is analyzed, and the fragmentation influence coefficient is obtained; during the ultrasonic scanning process, the acoustic pressure and microbubble movement data of the ultrasonic contrast imaging are obtained, the influence of the acoustic pressure on the microbubble movement is analyzed, and the movement influence coefficient is obtained; according to the fragmentation influence coefficient and the movement influence coefficient, the acoustic pressure is dynamically regulated to obtain the regulated acoustic pressure; then, according to the damage effect of the acoustic wave on the microbubble during the scanning time, the current microbubble vulnerability coefficient is analyzed; the regulated acoustic pressure is corrected by the current microbubble vulnerability coefficient to obtain the optimized regulated acoustic pressure. pressure; and then adjust the tracking area according to the current microbubble rupture rate and the current microbubble vulnerability coefficient; in the process of visual tracking of the dynamic distribution of ultrasound contrast agent, the present invention analyzes the influence of acoustic pressure distribution on microbubble fragmentation and the influence on microbubble movement, comprehensively analyzes the influence of microbubble fragmentation and microbubble movement, and combines the damage of acoustic pressure to microbubbles in the historical process, so as to accurately control the acoustic pressure during the scanning process; thereby extending the life of microbubbles, maintaining stable imaging, reducing nonspecific rupture, reducing biological risks, and effectively ensuring the stability of the tracking area, which will greatly improve the imaging quality in ultrasound contrast imaging.

[0090] Example 2

[0091] like Figure 2 As shown, an embodiment of the present invention provides a visualization tracking system for the dynamic distribution of ultrasound contrast agents, which specifically includes the following modules:

[0092] Acoustic pressure first impact module: During ultrasound scanning, the acoustic pressure and microbubble fragmentation data of ultrasound contrast imaging are obtained, the influence of acoustic pressure on microbubble fragmentation is analyzed, and the fragmentation influence coefficient is obtained;

[0093] More specifically, during the ultrasound scanning process, the acoustic pressure of ultrasound contrast imaging is obtained by a hydrophone, as well as the total number of initially injected microbubbles and the number of ruptured microbubbles. The percentage of the number of ruptured microbubbles to the total number of initially injected microbubbles is calculated to obtain the microbubble rupture rate.

[0094] The microbubble rupture rate and sound pressure usually present an S-shaped curve. The influence of sound pressure on microbubble rupture is analyzed by formula (1) to obtain the fragmentation influence coefficient;

[0095]

[0096] Where R max is the maximum rupture rate; P O is the sound pressure corresponding to the half-maximum rupture rate (critical threshold), P is the sound pressure of contrast-enhanced ultrasound; K1 is the fragmentation influence coefficient, and R is the goodness of fit;

[0097] Acoustic pressure second impact module: During ultrasound scanning, it obtains the acoustic pressure and microbubble movement data of ultrasound contrast imaging, analyzes the impact of acoustic pressure on microbubble movement, and obtains the movement influence coefficient;

[0098] In more detail, during the ultrasound scanning process, the acoustic pressure of ultrasound contrast imaging is obtained by a hydrophone, and the microbubble movement data (optical flow method or particle image velocimetry) is obtained. The influence of the acoustic pressure on the microbubble movement is analyzed by formula (2) to obtain the movement influence coefficient;

[0099] ΔV=K2*P norm +A formula (2)

[0100] Where, ΔV is the change in microbubble migration velocity, K2 is the migration influence coefficient, Pnorm is the normalized sound pressure of contrast-enhanced ultrasound, and A is the intercept;

[0101] Optimization control module: dynamically adjusts the sound pressure according to the crushing influence coefficient and the movement influence coefficient to obtain the regulated sound pressure;

[0102] Then, according to the damage effect of the acoustic wave on the microbubbles during the scanning time, the vulnerability coefficient of the current microbubbles is analyzed;

[0103] The regulated sound pressure is corrected by the current microbubble vulnerability coefficient to obtain the optimized regulated sound pressure;

[0104] Then, the tracking area is adjusted according to the current microbubble rupture rate and the current microbubble vulnerability coefficient;

[0105] In more detail, first, specifically, the specific process of dynamically regulating the sound pressure according to the fragmentation influence coefficient and the movement influence coefficient is:

[0106] Obtain the crushing influence coefficient and the movement influence coefficient, and construct the objective function J = a1*|C(t)-C target |+a2*|V(t)-V target |, where J is the objective function, which aims to accurately control the microbubble concentration and velocity under the premise of satisfying the microbubble survival rate (avoiding excessive rupture) to maximize the imaging quality or treatment effect, C(t) is the microbubble concentration, V(t) is the microbubble velocity, and C target is the preset target value of microbubble concentration, V target is the preset target value of microbubble velocity, a1 and a2 are weight coefficients, balancing the concentration and velocity priorities;

[0107] And set the constraints, the sound pressure range P min ≤P≤P max ; Rupture rate safety threshold R≤R threshold ;

[0108] Based on the fragmentation influence coefficient and the movement influence coefficient, the microbubble concentration and velocity in the future time window are predicted, and the optimized sound pressure is calculated through the PID control model; wherein, the PID control model is:

[0109]

[0110] Where Kp is the proportional gain, Ki is the integral gain, Kd is the differential gain, and e(t) = |C(t)-C target |+|V(t)-V target |, is the error signal;

[0111] Secondly, specifically, according to the damage effect of the acoustic wave on the microbubbles during the scanning time, the process of analyzing the current microbubble vulnerability coefficient is as follows:

[0112] The acoustic pressure is normalized, and the historical data is segmented by time window. The damage caused by a single pulse to the microbubble is calculated to obtain the single pulse damage. The decay rate of the microbubble damage over time is also obtained. All single pulse damage and decay rates within the scanning period are comprehensively analyzed to obtain the cumulative damage. The cumulative damage is mapped to obtain the vulnerability coefficient of the current microbubble.

[0113] Thirdly, specifically, the calculation process for optimizing and controlling the sound pressure is:

[0114] The current microbubble vulnerability coefficient is multiplied by the control sound pressure to obtain the optimized control sound pressure;

[0115] Fourthly, specifically, the process of adjusting the tracking area according to the current microbubble rupture rate and the current microbubble vulnerability coefficient is:

[0116] Obtain the inter-frame change rate of the second harmonic (2fo) signal intensity within the ROI; calculate the adjustment value of the tracking area using the ROI width adjustment formula;

[0117]

[0118] Wherein, DR is the adjustment value of the tracking region, ΔI is the inter-frame change rate, and W0 is the current ROI width.

[0119] The above is a detailed description of an embodiment of the present invention. However, the content described is only a preferred embodiment of the present invention and should not be considered to limit the scope of the present invention. All equivalent changes and improvements made within the scope of the present invention should still fall within the scope of the patent coverage of the present invention.

Claims

1. A method for visually tracking the dynamic distribution of an ultrasound contrast agent, characterized in that: The following steps are involved: During ultrasound scanning, the acoustic pressure and microbubble fragmentation data of ultrasound contrast imaging are obtained, the effect of acoustic pressure on microbubble fragmentation is analyzed, and the fragmentation influence coefficient is obtained; During ultrasound scanning, the acoustic pressure and microbubble movement data of ultrasound contrast imaging are obtained, the influence of acoustic pressure on microbubble movement is analyzed, and the movement influence coefficient is obtained; According to the crushing influence coefficient and the movement influence coefficient, the sound pressure is dynamically regulated to obtain the regulated sound pressure; Then, according to the damage effect of the acoustic wave on the microbubbles during the scanning time, the vulnerability coefficient of the current microbubbles is analyzed; The regulated sound pressure is corrected by the current microbubble vulnerability coefficient to obtain the optimized regulated sound pressure; Then, the tracking area is adjusted according to the current microbubble rupture rate and the current microbubble vulnerability coefficient.

2. The method for visually tracking the dynamic distribution of an ultrasound contrast agent according to claim 1, characterized in that: The process of obtaining the crushing influence coefficient is as follows: Obtain the total number of initially injected microbubbles and the number of ruptured microbubbles, and calculate the percentage of the number of ruptured microbubbles to the total number of initially injected microbubbles to obtain the microbubble rupture rate; The microbubble rupture rate and sound pressure present an S-shaped curve. The influence of sound pressure on microbubble rupture is analyzed and the fragmentation influence coefficient is obtained.

3. The method for visually tracking the dynamic distribution of an ultrasound contrast agent according to claim 1, characterized in that: The process of obtaining the mobile influence coefficient is as follows: Acquire microbubble movement data, analyze the effect of sound pressure on microbubble movement, and obtain the movement influence coefficient.

4. The method for visually tracking the dynamic distribution of an ultrasound contrast agent according to claim 1, characterized in that: According to the fragmentation influence coefficient and the movement influence coefficient, the specific process of dynamically controlling the sound pressure is as follows: According to the breakup influence coefficient and the movement influence coefficient, the microbubble concentration and velocity in the future time window are predicted, and the regulated sound pressure is calculated through the PID control model; wherein, the PID control model is: Among them, K p is the proportional gain, K i is the integral gain, K d is the differential gain, e(t) is the error signal, and P(t) is the regulated sound pressure.

5. The method for visually tracking the dynamic distribution of an ultrasound contrast agent according to claim 4, characterized in that: The expression of e(t) is: e(t)=|C(t)-C target |+|V(t)-V target |.

6. The method for visually tracking the dynamic distribution of an ultrasound contrast agent according to claim 1, characterized in that: According to the damage effect of acoustic waves on microbubbles during the scanning time, the process of analyzing the current microbubble vulnerability coefficient is as follows: The sound pressure is normalized, and the historical data is segmented by time window. The damage caused by a single pulse to the microbubble is calculated to obtain the single pulse damage. The decay rate of the microbubble damage over time is also obtained. All single pulse damage and decay rates within the scanning period are comprehensively analyzed to obtain the cumulative damage. The cumulative damage is mapped to obtain the vulnerability coefficient of the current microbubble.

7. The method for visually tracking the dynamic distribution of an ultrasound contrast agent according to claim 6, characterized in that: Comprehensively analyzing all single pulse damage and attenuation rates during the scanning period, the cumulative damage process is obtained as follows: The damage Dj of a single pulse to the microbubble is multiplied by the decay rate of the microbubble damage over time to obtain the single damage. The cumulative damage is obtained by summing up all the single damages within the scanning time.

8. The method for visually tracking the dynamic distribution of an ultrasound contrast agent according to claim 1, characterized in that: The calculation process for optimizing and controlling the sound pressure is: The current microbubble fragility coefficient is multiplied by the control sound pressure to obtain the optimized control sound pressure.

9. The method for visually tracking the dynamic distribution of an ultrasound contrast agent according to claim 8, characterized in that: According to the current microbubble rupture rate and the current microbubble vulnerability coefficient, the process of adjusting the tracking area is as follows: Obtain the inter-frame change rate of the second harmonic signal intensity within the ROI; calculate the adjustment value of the tracking area using the ROI width adjustment formula; Wherein, DR is the adjustment value of the tracking region, ΔI is the inter-frame change rate, and W0 is the current ROI width.

10. A visualization tracking system for the dynamic distribution of ultrasound contrast agents, characterized in that: The system is used to execute the method according to any one of claims 1 to 9, and the system comprises: Acoustic pressure first impact module: During ultrasound scanning, the acoustic pressure and microbubble fragmentation data of ultrasound contrast imaging are obtained, the influence of acoustic pressure on microbubble fragmentation is analyzed, and the fragmentation influence coefficient is obtained; Acoustic pressure second impact module: During ultrasound scanning, it obtains the acoustic pressure and microbubble movement data of ultrasound contrast imaging, analyzes the impact of acoustic pressure on microbubble movement, and obtains the movement influence coefficient; Optimization control module: dynamically adjusts the sound pressure according to the crushing influence coefficient and the movement influence coefficient to obtain the regulated sound pressure; Then, according to the damage effect of the acoustic wave on the microbubbles during the scanning time, the vulnerability coefficient of the current microbubbles is analyzed; The regulated sound pressure is corrected by the current microbubble vulnerability coefficient to obtain the optimized regulated sound pressure; Then, the tracking area is adjusted according to the current microbubble rupture rate and the current microbubble vulnerability coefficient.

Images